有限元分析理论基础.docx
- 文档编号:15184687
- 上传时间:2023-07-02
- 格式:DOCX
- 页数:23
- 大小:163.65KB
有限元分析理论基础.docx
《有限元分析理论基础.docx》由会员分享,可在线阅读,更多相关《有限元分析理论基础.docx(23页珍藏版)》请在冰点文库上搜索。
有限元分析理论基础
有限元理论基础
有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。
釆用不同的权函数和插值函数形式,便构成不同的有限元方法。
4.加权余量法:
是指采用使余量的加权函数为零求得微分方程近似解的方法称为加权余量法。
(WeightedresidualmethodWRM)是一种直接从所需求解的微分方程及边界条件出发,寻求边值问题近似解的数学方法。
加权余量法是求解微分方程近似解的一种有效的方法。
设问题的控制微分方程为:
在V域内厶(")-八0(5.1.1)
在S边界上〃(“)-&=0(5.1.2)
式中:
L、B——分别为微分方程和边界条件中的微分算子;
f、g——为与未知函数u无关的己知函数域值;
u——为问题待求的未知函数
当弄!
J用力u权余•肚法求近丁以解首先在求耳军域上理立一个T式閑数H一般兵升如下形式:
仁土CN=NC(5.1.3)
Tm
式中:
c{彳寺定系数.也可称为广义坐标;
N:
取白完备函冬攵*S线.性无关的基函孕攵°
由于〃一般只圮彳守求函缨攵U的近1以耳岂因u匕将式(51.3)代入式(511)牙口式(51.2)后将诃•不誉斯兄,昔迅:
|R]=L(flb—f在V域内
\rb=B(^~g在S边界上("14)
城然&、尽反映了r式函竽攵与实解之问的偏差.它丁门分另U称
做内召卩牙口边界余覺。
若在域\'内引入内部权函数硏,在边界S上引入边界权函数WB则可理立11个消除余甘的条件.一般可农示为:
L兀WB1RBdS=0(/=L2.L,〃)(51-5)
•V•S
不同的权函数幵;和jrR反映了不同的消除余•眩的准则。
从上式可以得到求解待定系数矩阵C的代数方程组」一经解得待定系数.由式(5.1.3)即可得所需求解边值问题的进似解」由于试函数〃的不同■余址A和Rb可冇如下三种T/r况.
依”匕加权余•竝法可分为:
1.内部法
试函数满尺边界条件,也即《=B©-g=°此时消除氽呈的条件成为:
1;%尺"=°072丄E)(5.16)
2.边界法
试函数满足控制r方程.也即R■砂CO此时消除氽址的条件为:
“詁/£二°(心L2Lji)(5.1.7)
•s
3.混合法
试函数不测兄控制方程和边界条件.此时用式(51.5)来消除氽显、
限元分析理论基础
有限元分析概念
有限元法:
把求解区域看作由许多小的在节点处相互连接的单元(子域)所构成,其模型给出基本方程的分片(子域)近似解,由于单元(子域)可以被分割成各种形状和大小不同的尺寸,所以它能很好地适应复杂的几何形状、复杂的材料特性和复杂的边界条件
有限元模型:
它是真实系统理想化的数学抽象。
由一些简单形状的单元组成,单元之间通过节点连接,并承受一定载荷。
有限元分析:
是利用数学近似的方法对真实物理系统(几何和载荷工况)进行模拟。
并利用简单而又相互作用的元素,即单元,就可以用有限数量的未知量去逼近无限未知量的真实系统。
线弹性有限元是以理想弹性体为研究对象的,所考虑的变形建立在小变形假设的基础上。
在这类问题中,材料的应力与应变呈线性关系,满足广义胡克定律;应力与应变也是线性关系,线弹性问题可归结为求解线性方程问题,所以只需要较少的计算时间。
如果采用高效的代数方程组求解方法,也有助于降低有限元分析的时间。
线弹性有限元一般包括线弹性静力学分析与线弹性动力学分析两方面。
非线性问题与线弹性问题的区别:
1)非线性问题的方程是非线性的,一般需要迭代求解;
2)非线性问题不能采用叠加原理;
3)非线性问题不总有一致解,有时甚至没有解。
有限元求解非线性问题可分为以下三类:
1)材料非线性问题
材料的应力和应变是非线性的,但应力与应变却很微小,此时应变与位移呈线性关系,这类问题属于材料的非线性问题。
由于从理论上还不能提供能普遍接受的本构关系,所以,一般材料的应力与应变之间的非线性关系要基于试验数据,有时非线性材料特性可用数学模型进行模拟,尽管这些模型总有他们的局限性。
在工程实际中较为重要的材料非线性问题有:
非线性弹性(包括分段线弹性)、弹塑性、粘塑性及蠕变等。
2)几何非线性问题
几何非线性问题是由于位移之间存在非线性关系引起的。
当物体的位移较大时,应变与位移的关系是非线性关系。
研究这类问题一般都是假定材料的应力和应变呈线性关系。
它包括大位移大应变及大位移小应变问题。
如结构的弹性屈曲问题属于大位移小应变问题,橡胶部件形成过程为大应变问题。
3)非线性边界问题
在加工、密封、撞击等问题中,接触和摩擦的作用不可忽视,接触边界属于高度非线性边界。
平时遇到的一些接触问题,如齿轮传动、冲压成型、轧制成型、橡胶减振器、紧配合装配等,当一个结构与另一个结构或外部边界相接触时通常要考虑非线性边界条件。
实际的非线性可能同时出现上述两种或三种非线性问题。
有限元理论基础
有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。
采用不同的权函数和插值函数形式,便构成不同的有限元方法。
1・加权余量法:
是指采用使余量的加权函数为零求得微分方程近似解的方法称为加权余量法。
(WeightedresidualmethodWRM)是一种直接从所需求解的微分方程及边界条件出发,寻求边值问题近似解的数学方法O加权余量法是求解微分方程近似解的一种有效的方法。
设问题的控制微分方程为:
在V域内L(u)-f=O(5.1.1)
在S边界上Bg-g=0(5.1.2)
式中:
L、B——分别为微分方程和边界条件中的微分算子;
f、9——为与未知函数U无关的己知函数域值;
U——为问题待求的未知函数
当牙!
J用力口权余试法求近T以解日寸,甘先在求解喊上理立一个壬式函数〃一般兵冇如帀形式:
z?
o=VC\Nt=AV(5」3)
式中:
G钳定系数.她可称为广义坐标;
N、J队白完备函数共的线J生无关白勺基函数二
由于R一般只足彳寺求函掬;u的近似解.因D匕将•式(5.1.3)代入式(5.11)芽口式(5.12)后将•得•不至U满足.若2:
、Ri=L(M_.f在V域内
\RB=B(ft^-g在S边界上
显然耳、為反映了试《函数与為卖稱之川1白勺偏差,它T门分另州尔
做内占卩才口边界余Tib,
若在域¥內引入内部权函数为,在边界S上引入边界权函数WB则可建立11个消除氽量的条件,一般可农示为:
|r,^r+|rs)A^zs=o(/=l2,l,〃)er
♦p*s
不同的权函数〃;和%反映了不同的消除余-址的准贝!
J。
从上式可以得到求解诗定系数矩阵C的代数方程组。
一经解得待定系数,由式(5.13)即可得所需求稱边值问題的近似解〉由于试函数"的不同.余壁&和凫可方如下三种‘侨况,
依此加权余崔•法可分为:
1.内部法
试函数满兄边界条件.她即&7衢甘Q此时消除余壁的条件成为:
I:
麻"卩=0Q=1、2.L")(5.1.6)
2.边界法
试函数满足控制方程,电卯只=乂$7-/=0此时消除金量的条件为:
曲阳洽=0(/=1.2XJi)(S-i-7)
3.混合法
试函数不濒叉控制方程和边界条件.此时用式(5.1.5)来消除余蚤°
混合法对于试函数的选取最方便,但在相同精度条件下,工作量最大。
对内部法和边界法必须使基函数事先满足一定条件,这对复杂结构分析往往有一定困难,但试函数一经建立,其工作量较小。
无论采用何种方法,在建立试函数时均应注意以下几点:
(1)试函数应由完备函数集的子集构成。
已被采用过的试函数有幕级数、三角级数、样条函数、贝赛尔函数、切比雪夫和勒让德多项式等等。
(2)试函数应具有直到比消除余量的加权积分表达式中最高阶导数低一阶的导数连续性。
(3)试函数应与问题的解析解或问题的特解相关联。
若计算问题具有对称性,应充分利用它。
显然,任何独立的完全函数集都可以作为权函数。
按照对权函数的不同选择得到不同的加权余量计算方法,主要有:
配点法、子域法、最小二乘法、力矩法和伽辽金法。
其中伽辽金法的精度最高。
下而以内部法为例,介绍按权函数分类时加权余蚤的五种基本方法。
对内部法来说,消除余蚤的统一格式是:
0%亦=0(/=1,2.L屮)
1.子域法(SiibdoinninMethod)
此法首先M各求解域丫划分成n个子域%,衣宓个子域内令权函数等于1.而在子域之外取权函数为零,也艮P:
卩(豐)
”|0代外)
如果在各个子域里分别选取试函数.那么它的求禅在形式上将
2.配点法(CollocationMethod)
子域法是令余呈在一个子域上的总和为冬。
而配点法是使
余•豪在扌旨定的n个点上等于寺,这些点称为配点。
0匕法的杈
函数为:
Wh=&P_P)
ooo
篦羽诚£(犹拉克)函数,它的定义为:
P、Pj—分别代表求稱域内任一点和色己点。
由于此法只在配点上偌证余量为零,因此不需耍作积分计算,
所以是最简单的加权余蚤法
3.最小二乘法(LeastSquareMethod)
本法立旌立无丫更在涎个求解域上余-处的平方和取极力、来建立消除余量的条件。
若记余泌平方和为1(C),即2(C)=f砖如=IR]R,dV
则极值条件为:
2厂(箸)'巳沙=°
由此可见,本法权函数为:
昭=淫0=L21,«)
4.伽辽金法Method)
本法是使余涅'与每一个基函数正交,也即以基函数作为权函数
Wh=N((/=1,2,Lji)
当•£式函数田包含注个完备函数集时,用本法必可衣得耕确解
5.矩法(MethodofMoment)
本法与伽辽金法相似,也是用完备函数集作权函数。
但本法的权函数与伽辽金法又冇区别,它与f式函数夭关。
消除余•昼的条件是从岑开始的各阶夕巨为冬.因此
对一^维问题,W=兀nU=l,2,L.n)
对二维丿诃题叫=“》八o’」=12l卫)
其余类推
这五种基本方,去在待定系数尺够多(牙尔做高阶近似)时,其榜皮彼此相近。
但对低阶近似(口较刀寸忻况下.后三科的精皮要高于前两种。
基本方法芈例
为说明上述基本棋念.以图所示等截面悬皆梁.淡満跨均布荷
载作用,求悬臂端B的竖向位移心为例,说明基本方法的应用。
图示采•的扌空制方程为:
怙=。
其边界条件为:
(x=0)
若取试函数为:
c(x5+lx4-14/2.v3+26l3x2)(*)
(x=0
不难验证其'萌足边界条件,也叩^=0o而控制方程的内部余苛R?
为:
A;=E/c(120.v+24/)-<7
因此本问题屈内郃法。
下而•步别用基本方法进行求解。
子域法解
由于f式函数仅一个待定常'数,因此只需取一个子域(芳于全域)即可.消除余•昼的条件为:
[呼(12O.y+24/)-必二0
由此可解得:
「—q
S4EI1
代回(巧式可得•:
忆点法解
=0
同上所述.只需选一个配点来建立消除余呈的条件」若令:
R】
可得:
G=Q峻用
11457/丄5752
若令:
RiL=°
贝呵
c-q、应)(恂
144512方归助
可见不同的色己点结果是不一样的。
最小二乘法解
此时消除余量的条件为:
r為簧扛二\[EIci,120%+24/I- 伽辽金法解 U匕时,N、=x5+/x4-14/2x3+26ft2 消除余呈的条件为: 仙也=() 矩法解 由于只冇一个待定常数.因此消除余量条件只需零次矩即可.此时显然与子域法完全相同。 本例各方法的精度比较 本问越的粉•确解由三養位移计算可■停为: _ql4_0.1257/4 _SEI~EI 由U匕可彳孚,上述各方法对本T列计算的i吴差依次为: ・33.3%: 1.75%(22.2%);13.9%;0.96%;-33.3%上面22.2%为式(林)结巣。 2、虚功原理 ——平衡方程和几何方程的等效积分“弱”形式 虚功原理包含虚位移原理和虚应力原理,是虚位移原理和虚应力原理的总称。 他们都可以认为是与某些控制方程相等效的积分“弱”形式。 虚功原理: 变形体中任意满足平衡的力系在任意满足协调条件的变形状态上作的虚功等于零,即体系外力的虚功与内力的虚功之和等于零。 虚位移原理是平衡方程和力的边界条件的等效积分的“弱”形式; 虚应力原理是几何方程和位移边界条件的等效积分“弱”形式。 虚位移原理的力学意义: 如果力系是平衡的,则它们在虚位移和虚应变上所作的功的总和为零。 反之,如果力系在虚位移(及虚应变)上所作的功的和等于零,则它们一定满足平衡方程。 所以,虚位移原理表述了力系平衡的必要而充分条件。 一般而言,虚位移原理不仅可以适用于线弹性问题,而且可以用于非线性弹性及弹塑性等非线性问题。 虚应力原理的力学意义: 如果位移是协调的,则虚应力和虚边界约束反力在他们上面所作的功的总和为零。 反之,如果上述虚力系在他们上面所作的功的和为零,则它们一定是满足协调的。 所以,虚应力原理表述了位移协调的必要而充分条件。 虚应力原理可以应用于线弹性以及非线性弹性等不同的力学问题。 但是必须指出,无论是虚位移原理还是虚应力原理,他们所依赖的几何方程和平衡方程都是基于小变形理论的,他们不能直接应用于基于大变形理论的力学问题。 3、最小总势能法 应变能: 作用在物体上的外载荷会引起物体变形,变形期间外力所做的功以弹性能的形式储存在物体中,即为应变能。 由n个单元和m个节点组成的物体的总势能为总应变能和外力所做功的差: nm r-11-1 最小势能原理: 对于一个稳定的系统,相对于平衡位置发生的位移总会使系统的总势能最小,BP: 斜貓*遗汕=。 ,曰,2,3,……,n 有限元法的收敛性 有限元法是一种数值分析方法,因此应考虑收敛性问题。 有限元法的收敛性是指: 当网格逐渐加密时,有限元解答的序列收敛到精确解;或者当单元尺寸固定时,每个单元的自由度数越多,有限元的解答就越趋近于精确解。 有限元的收敛条件包括如下四个方面: 1)单元内,位移函数必须连续。 多项式是单值连续函数,因此选择多项式作为位移函数,在单元内的连续性能够保证。 2)在单元内,位移函数必须包括常应变项。 每个单元的应变状态总可以分解为不依赖于单元内各点位置的常应变和由各点位置决定的变量应变。 当单元的尺寸足够小时,单元中各点的应变趋于相等,单元的变形比较均匀,因而常应变就成为应变的主要部分。 为反映单元的应变状态,单元位移函数必须包括常应变项。 3)在单元内,位移函数必须包括刚体位移项。 一般情况下,单元内任一点的位移包括形变位移和刚体位移两部分。 形变位移与物体形状及体积的改变相联系,因而产生应变;刚体位移只改变物体位置,不改变物体的形状和体积,即刚体位移是不产生变形的位移。 空间一个物体包括三个平动位移和三个转动位移,共有六个刚体位移分量。 由于一个单元牵连在另一些单元上,其他单元发生变形时必将带动单元做刚体位移,由此可见,为模拟一个单元的真实位移,假定的单元位移函数必须包括刚体位移项。 4)位移函数在相邻单元的公共边界上必须协调。 对一般单元而言,协调性是指相邻单元在公共节点处有相同的位移,而且沿单元边界也有相同的位移,也就是说,要保证不发生单元的相互脱离开裂和相互侵入重叠。 要做到这一点,就要求函数在公共边界上能由公共节点的函数值唯一确定。 对一般单元,协调性保证了相邻单元边界位移的连续性。 但是,在板壳的相邻单元之间,还要求位移的一阶导数连续,只有这样,才能保证结构的应变能是有界量。 总的说来,协调性是指在相邻单元的公共边界上满足连续性条件。 前三条又叫完备性条件,满足完备条件的单元叫完备单元;第四条是协调性要求,满足协调性的单元叫协调单元;否则称为非协调单元。 完备性要求是收敛的必要条件,四条全部满足,构成收敛的充分必要条件。 在实际应用中,要使选择的位移函数全部满足完备性和协调性要求是比较困难的,在某些情况下可以放松对协调性的要求。 需要指出的是,有时非协调单元比与它对应的协调单元还要好,其原因在于近似解的性质。 假定位移函数就相当于给单元施加了约束条件,使单元变形服从所加约束,这样的替代结构比真实结构更刚一些。 但是,这种近似结构由于允许单元分离、重叠,使单元的刚度变软了,或者形成了(例如板单元在单元之间的绕度连续,而转角不连续时,刚节点变为较接点)对于非协调单元,上述两种影响有误差相消的可能,因此利用非协调单元有时也会得到很好的结果。 在工程实践中,非协调元必须通过“小片试验后”才能使用。 应力的单元平均或节点平均处理方法 最简单的处理应力结果的方法是取相邻单元或围绕节点各单元应力的平均值。 ・1•取相邻单元应力的平均值 这种方法最常用于3节点三角形单元中。 这种最简单而又相当实用的单元得到的应力解在单元内是常数。 可以将其看作是单元内应力的平均值,或是单元形心处的应力。 由于应力近似解总是在精确解上下振荡,可以取相邻单元应力的平均值作为此两个单元合成的较大四边形单元形心处的应力。 如2单元的情况下,取平均应力可以采用算术平均, 即平均应力二(单元1的应力+单元2的应力)12. 也可以釆用精确一些的面积加权平均, 即平均应力二[单元1应力X单元1的面积+单元2应力X单元2面积]/(单元1面积+单元2面积) 当相邻两单元面积相差不大时,两者的结果基本相同。 在单元划分时应避免相邻两单元的面积相差太多,从而使求解的误差相 近。 一般而言,3节点三角形单元的最佳应力点是单元的中心点,此点的应力具有1阶的精度。 ・2•取围绕节点各单元应力的平均值 首先计算围绕该节点(i)周围的相关单元在该节点出的应力值, 1加 然后以他们的平均值作为该节点的最后应力值即6=丄工bmW 其中,仃m是围绕在i节点周围的全部单元。 取平均值时也可进行面积加权。 有限元法求解问题的基本步骤 1•结构离散化 对整个结构进行离散化,将其分割成若干个单元,单元间彼此通过节点相连; 2•求出各单元的刚度矩阵[KF [K]⑹是由单元节点位移量{①严)求单元节点力向量{F}(e)的转移矩 阵,其关系式为: {F}(4[K](e){O}(e) 3•集成总体刚度矩阵[K]并写出总体平衡方程: 总体刚度矩阵[K]是由整体节点位移向量{①}求整体节点力向量的转移矩阵,其关系式为{F}=[K]{O},此即为总体平衡方程。 4•引入支撑条件,求出各节点的位移 节点的支撑条件有两种: 一种是节点n沿某个方向的位移为零,另一种是节点n沿某个方向的位移为一给定值。 5•求出各单元内的应力和应变。 对于有限元方法,其基本思路和解题步骤可归纳为: (1)建立积分方程,根据变分原理或方程余量与权函数正交化原理,建立与微分方程初边值问题等价的积分表达式,这是有限元法的岀发点。 (2)区域单元剖分,根据求解区域的形状及实际问题的物理特点, 将区域剖分为若干相互连接、不重叠的单元。 区域单元划分是釆用有限 元方法的前期准备工作,这部分工作量比较大,除了给计算单元和节点进行编号和确定相互之间的关系之外,还要表示节点的位置坐标,同时还需要列出自然边界和本质边界的节点序号和相应的边界值。 (3)确定单元基函数,根据单元中节点数目及对近似解精度的要求,选择满足一定插值条件的插值函数作为单元基函数。 有限元方法中的基函数是在单元中选取的,由于各单元具有规则的几何形状,在选取基函数时可遵循一定的法则。 (4)单元分析: 将各个单元中的求解函数用单元基函数的线性组合表达式进行逼近;再将近似函数代入积分方程,并对单元区域进行积分,可获得含有待定系数(即单元中各节点的参数值)的代数方程组,称为单元有限元方程。 (5)总体合成: 在得出单元有限元方程之后,将区域中所有单元有限元方程按一定法则进行累加,形成总体有限元方程。 (6)边界条件的处理: 一般边界条件有三种形式,分为本质边界条件(狄里克雷边界条件)、自然边界条件(黎曼边界条件)、混合边界条件(柯西边界条件)。 对于自然边界条件,一般在积分表达式中可自动得到满足。 对于本质边界条件和混合边界条件,需按一定法则对总体有限元方程进行修正满足。 (7)解有限元方程: 根据边界条件修正的总体有限元方程组,是含所有待定未知量的封闭方程组,采用适当的数值计算方法求解,可求得各节点的函数值。 单元刚度矩阵的特性 单元刚度矩阵无论在局部坐标系中还是在整体坐标系中都具有相同的三个特性: 1)对称性 由材料力学中的位移互等定理可知,对一个构件,作用在点j的力引起点i的绕度等于有同样大小而作用于点i的力引起的点j的绕度,即kij(e)=kji(e),表明单元刚度矩阵是一个对称矩阵。 2)奇异性 无逆阵的矩阵就叫做奇异矩阵,其行列式的值为0,即1以叫=(),这一点可以从例题直接得到验证。 其物理意义是引入支撑条件之前,单元可 平移。 3)分块性 有前面所讲的内容可以看出,矩阵[k®]可以用虚线分成四块,因此可写成如下的分块形式, 式中际泸)—局部坐标系中单元(e)按局部码标记的节点hi、n之间的刚度子矩阵 刚架结构中非节点载荷的处理的方法 在刚架结构以及其他较复杂的结构上,他们所受的载荷可以直接作用在节点上,又可以不直接作用在节点上而作用于单元节点间的其他位置上。 后一种情况下的载荷称为非节点载荷。 有限元分析时,总体刚度方程中所用到的力向量是节点力向量。 因此在进行整体分析前应当进行载荷的移植,将作用于单元上的力移植到节点上。 移植时按静力等效的原则进行。 处理非节点载荷一般可直接在整体坐标系内进行,其过程为: 1)将各杆单元看成一根两端固定的梁,分别求出两个固定端的约束反力。 其结果可直接利用材料力学的公式求得; 2)将各固定端的约束反力变号,按节点进行集成,获得各节点的等效载荷 总体刚度矩阵的集成法 使用刚度矩阵获得的方法获得总体刚度矩阵。 在此将其扩展到由整体坐标系中的单元刚度矩阵的子矩阵集成总体刚度矩阵。 步骤如下: 1)对一个有n个节点的结构,将总体刚度矩阵[K]划分为nXn各子区间,然后按节点总码的顺序进行编号; 2)将整体坐标系中单元刚度矩阵的各子矩阵根据其下标的两个总码对号入座,写在总体刚度矩阵相应的子区间; 3)同一子区间内的子矩阵相加,成为总体刚度矩阵中的相应的子矩阵。 总体刚度矩阵的特性 1)对称性: 因为由此特性,在计算机中只需存储其上三角部
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 有限元分析 理论基础