桥梁软件应用有限元简介

1、第二章第二章 结构分析的有限元法结构分析的有限元法 2.1 2.1 有限元法发展简况有限元法发展简况 利用定义在三角形区域上的分片连续函数和最小位能原理St.Venant扭转问题的近似解有限元法的研究现代有限元法1943Courant应用数学家、物理学家、工程师1960Tumer、Clough第一次用三角形单元平面应力问题解答提出了有限单元法的名称各种非线性问题多物理场耦合问题多尺度问题商品化有限元软件20世纪70年代国外几何非线性：因几何变形引起结构刚度改变材料非线性：弹性（超弹和多线性弹性）、粘弹性、非弹性状态非线性：接触问题2.1 有限元法发展简况 固体力学流体力学传热学电磁学学科应用力

2、学计算结构优化计算功能计算技术纯粹数值技术前、后处理技术的高度智能化和与CAD的集成化2.2 2.2 有限元法的基本思路及其求解步骤有限元法的基本思路及其求解步骤 经典的解析法 从连续体的微分方程入手，寻求满足微分方程和定解条件的适合全域的解析解，一旦得到解析解，就可知道域内任意点的解大多数问题，特别是实际问题 很难甚至无法用解析法得到问题的解析解在整个求解域上满足控制方程在边界上满足边界条件的场函数寻找很困难很困难有限元法有限元法单元节点有限元模型2.2 有限元法的基本思路及其求解步骤 有限元法基本思路抛弃寻找一个满足整个求解域的场函数的思路把求解域划分成有限个四边形单元对每一个单元通过插值

3、的方法，用其节点上的位移建立该单元的位移函数123每个单元都有与其对应的位移函数表达式用全部单元域之和代替整个求解域，用全部单元的位移函数之和代替满足整个求解域的位移函数4对单元进行力学特性分析，建立单元节点力与单元节点位移的关系，并将结构的外载荷等效移植到节点上，再在节点上建立力的平衡方程，求解后得到节点上的位移，继而得到各个单元的应力5以以节点位移节点位移为未知量，通过为未知量，通过求解求解力的平衡方程力的平衡方程获得节点获得节点位移，然后按位移，然后按单元单元计算应力计算应力2.2 有限元法的基本思路及其求解步骤 有限元法求解步骤1 离散化离散化将结构（求解域）划分为有限个单元，让全部单

4、元的集合与原结构近似等价划分单元时，二者在几何形体上越逼近越好，特别是在位移和应力急剧变化的地方2 选择单元位移函数选择单元位移函数在有限元法中，需要用单元节点位移通过插值方法建立单元位移函数（单元位移模式），即用单元节点位移来描述单元位移。单元位移函数的合理与否，直接关系到有限元分析的计算精度、效率和收敛性。通常取为多项式形式2.2 有限元法的基本思路及其求解步骤 3 单元特性分析单元特性分析（1）依照应变与位移之间的几何关系，根据所选择的单元位移函数，建立单元应变与单元节点位移之间的关系式。据此式，在求出节点位移后，可以求得单元应变。（2）依照物理关系（胡克定律），建立单元应力与单元节点位

5、移之间的关系式。据此式，在求出节点位移后，可以求得单元应力。（3）根据虚位移原理或最小势能原理，建立单元刚度方程，即单元节点力与单元节点位移之间的关系式。此步骤核心是计算单元刚度矩阵。4 外载荷处理外载荷处理将外载荷（体力、面力等）等效移植到节点上。2.2 有限元法的基本思路及其求解步骤 5 建立节点上的力平衡方程建立节点上的力平衡方程按照有限元法的统一格式，形成如下形式的以节点位移为未知量的代数方程组 KF K由各个单元的刚度矩阵组装成的总体刚度矩阵 待求的节点位移列阵 F按节点编号顺序形成的节点载荷列阵6 处理边界条件、解算节点位移处理边界条件、解算节点位移(2.1)按照实际位移边界条件，

6、对式(2.1)进行整理，解之，可得单元节点位移。有了节点位移，即可根据单元特性分析中建立的关系式，求应力、应变、内力等。有了节点位移，即可根据单元特性分析中建立的关系式，求应力、应变、内力等。后处理后处理：对所选应力、应变等，以：对所选应力、应变等，以彩色云图彩色云图或或图表图表的形式显示的形式显示计算结果。计算结果。2.3 2.3 有限元程序的结构简介有限元程序的结构简介 对一个题目或一个实际工程问题进行有限元分析，大体上对一个题目或一个实际工程问题进行有限元分析，大体上分分3个个主要步骤主要步骤有限元建模有限元求解计算结果分析与整理前处理求解器后处理程程序序结结构构2.3 有限元程序的结构

7、简介 前处理几何模型的建立定义约束条件网格剖分确定材料参数和载荷有限元建模形成有限元分析所需用的有限元计算数据形成有限元分析所需用的有限元计算数据可视化可视化 有限元模型有限元模型前处理中，可以用图形显示所建立的几何模型、单元网格、约束条件等求解器2.3 有限元程序的结构简介 有限元程序的核心部分主要完成有限元模型的力学计算，即根据前处理形成的有限元计算数据，完成以下工作：计算单元刚度矩阵计算节点载荷组装总体刚度矩阵将载荷等效简化到节点上形成总体有限元平衡方程求解节点位移计算应力、应变、内力等2.3 有限元程序的结构简介 后处理根据计算者的要求对计算结果进行检查、分析、整理、打印输出等进行数据

8、检索响应量合成绘制变形图、应力图、应变图、曲线图等可视化可视化的方式分析、观察计算结果的方式分析、观察计算结果计算者计算者进行有限元分析的工作量主要体进行有限元分析的工作量主要体现在现在前处理前处理和和后处理后处理方面方面2.4 2.4 算例算例 2.4.1 平面三角形常应变单元任意区域三角形单元网 格剖分示意图mvmujvjuiviuijm典型三角形单元yxo单元内任意点(x,y)的位移uv、坐标x和y的函数建立单元位移函数通过插值方法建立，即用单元的节点位移来表示单元内任意点的位移1. 单元位移函数2.4 算例 mvmujvjuiviuijm典型三角形单元yxo单元位移函数选用坐标x和y的

9、一次多项式123456uxyvxy(1)123456、待定系数123123123iiijjjmmmuxyuxyuxy(2)未知量求解(2)(3)得到123456、123456、456456456iiijjjmmmvxyvxyvxy(3)2.4 算例 123123123iiijjjmmmuxyuxyuxy(2)1121iijjmmxyDxyAxy123112111121111121iiijjjiijjmmmmmiijjiijjmmmmiijjiijjmmmmuxyuxyaua ua uDAuxyuyuybub ub uDAuyxuxucuc uc uDAxu456456456iiijjjmmmv

10、xyvxyvxy(3)1121iijjmmxyDxyAxy456112111121111121iiijjjiijjmmmmmiijjiijjmmmmiijjiijjmmmmvxyvxyava va vDAvxyvyvybvb vb vDAvyxvxvcvc vc vDAxv2.4 算例 ijmmjjmiimmijjiax yx yax yx yax yx yijmjmimijbyybyybyyijmjmimijcxxcxxcxx 123456uxyvxy(1)代入123456、将求得的得到用单元的节点位移表示的单元位移函数iijjmmiijjmmuN uN uN uvN vN vN v(4)式

11、中是单元形状函数，简称形函数121212iiiijjjjmmmmNab xc yANab xc yANab xc yA（5）,ijmN NN是常数，取决于单元的三个节点坐标, , ,iiima b cc，返回P242.4 算例 11211111.()( 1).() 1.()()222iijjmmjmmjimmiijjiijmxyDxyAxyADx yx yx yx yx yx yaaa 1.()( 1).()1.()()()()()()()0.jmmjimmiijjijmimijijmjjmmiijjjmjmmmiiimimjjmiimjmimjiijmijjix yx yx yx yx yx

12、 yxyyxyyx yyb xb xb xb xb xb xb xb xb xb xbxxb xxbbb xb cbcxbcb c 1()2ijjiAbcb c三角形单元的面积A单元位移函数表达式iijjmmiijjmmuN uN uN uvN vN vN v(4)2.4 算例 000000iijeijmijmjmmuvuNNNuuNNNvvuv eeuN写成矩阵形式简写为其中， eu表示单元内任意点处位移的单元位移函数列阵 000000ijmijmNNNNNNN为形函数矩阵返回P282.4 算例 形函数的性质在节点上形函数的值是式（6）表示形函数Ni在其自身节点上的值等于1，在其他节点上的值

13、等于0，即1(,)( , ,)0ijjijjiN xyi j mji（6）( ,)1(,)0(,)0iiiijjimmN x yN xyN xy，( ,)0(,)1(,)0jiijjjjmmNx yNxyNxy，( ,)0(,)0(,)1miimjjmmmNx yNxyNxy，1单元中任意一点上的各个形函数之和等于1，即22.4 算例 ( , )( , )( , )1ijmN x yNx yNx y由121212iiiijjjjmmmmNab xc yANab xc yANab xc yA（5） ( , )( , )( , )11122212ijmiiijjjmmmijmijmijmN x y

14、Nx yNx yab xc yab xc yab xc yAAAaaabbbxcccyAijmjmimijbyybyybyyijmjmimijcxxcxxcxx 1()2ijmAaaa120.0.21AxyA2.4 算例 小结小结（1 1）本节的三角形单元，形函数是）本节的三角形单元，形函数是线性线性的，为的，为x x、y y的一的一次函数；次函数；（2 2）在单元内部和各条单元边上，位移也是）在单元内部和各条单元边上，位移也是线性线性的，可的，可由两个节点的位移由两个节点的位移唯一唯一确定；确定；（3 3）相邻单元的公共节点的节点位移是相等的，因此，）相邻单元的公共节点的节点位移是相等的，因

15、此，能保证相邻单元在公共边界上以及单元内部的能保证相邻单元在公共边界上以及单元内部的位移连续位移连续性性。2.4 算例 单元位移函数确定后，根据几何方程xyxyuxvyuvyx求得单元内任意点处的应变，即单元应变 000000ijimixjejimyjxyjjiimmmmuNNNuvxxxxuNNNvvyyyyuvNNNNNNuyxyxyxyxv（7）iijjmmiijjmmuN uN uN uvN vN vN v(4)2. 单元应变和单元应力2.4 算例 由121212iiiijjjjmmmmNab xc yANab xc yANab xc yA（5）形函数对坐标变量求偏导111222111

16、222jimijmjimijmNNNbbbxAxAxANNNcccyAyAyA，（8），式（8）代入式（7）中，得到 0001000(9)2iixijmjeeyijmjiijjmmxymmuvbbbucccBvAcbcbcbuv =返回P28由式（9）得到 0001000(10)2ijmijmijmiijjmmbbbBcccBBBAcbcbcb2.4 算例 单元应变矩阵（几何矩阵）000111000222ijmiijjmmiijjmmbbbBcBcBcAAAcbcbcb，（11）分块矩阵参数,; ,ijmijmb b bc c c由单元的节点坐标确定，因此，它们取决于单元形状，当单元的节点坐标

17、确定后，它们都是常量，所以，3节点三角形单元的应变矩阵B是常数矩阵2.4 算例 根据物理方程 xeeeeyxyDDBS（12）其中 000200101011002ED为平面应力（平面应变）问题的弹性矩阵平面应力问题00,EE平面应变问题002,11EE 0000000200000002(1)111111222222iijjmmiijjmmiijjmmbcbcbcESDBbcbcbcAcbcbcb（13）应力矩阵S也是常数矩阵单元应力矩阵返回P28（1 1）3 3节点平面三角形单元，节点平面三角形单元，应变矩阵应变矩阵和和应力矩阵应力矩阵都为都为常数矩阵常数矩阵；（2 2） 3节点平面三角形单元

18、，各点的应变和应力都是相节点平面三角形单元，各点的应变和应力都是相同的，且是常数，所以同的，且是常数，所以3节点三角形单元是节点三角形单元是常应变单元常应变单元，也是也是常应力单元常应力单元；（3 3）采用）采用3 3节点三角形单元时，在应力变化剧烈或应力节点三角形单元时，在应力变化剧烈或应力梯度较大的部位，单元划分应适当加密。梯度较大的部位，单元划分应适当加密。2.4 算例 小结小结3. 单元刚度矩阵2.4 算例 分析结构在载荷作用下产生变形和应力，于是在各单元之间就产生相互作用。实际上，各单元之间的相互作用是通过相邻边界上（即，单元的边，实际是面）的分布力而产生的。按照有限元方法，结构离散

19、化为一个个单元后，单元之间的相互作用就由单元的节点力来实现，即用单元节点力等效代替相邻边界上的相互作用力，这样，节点力就与单元应力相关，而单元应力与节点位移相关，因此，单元节点力与单元节点位移相关。建立单元节点力与单元节点位移之间的关系 eR表示单元节点力mvmujvjuiviuijm (14)Teiijjmmuvuvuu (15)TeixiyjxjymxmyRRRRRRR单元节点位移由（9）式2.4 算例 (9)eeB= eeDB（12）由（12）式把一个单元作为分析对象时，可以把节点力看作外力。单元节点力和单元节点位移之间的关系可由虚位移原理导出在外力作用下，处于平衡状态的变形体，当发生约

20、束允许的任意微小的虚位移时，外力在虚位移上所做的虚功等于整个体积内的应力在虚应变上所做的虚功。推导1 令单元的节点虚位移为 *eTiijjmmuvuvuv eeuN由P18*(16)eeuN eeB由(9)式 *(17)eeB2.4 算例 (15)TeixiyjxjymxmyRRRRRRR *eTiijjmmuvuvuv2 节点力在虚位移上所做的虚功为 *(18)Teeiixiiyjjxjjymmxmmyu Rv Ru Rv Ru Rv RR3 单元应力在虚应变上所做的虚功为 *(19)eTeeVdV单元体积 *(17)eeB将 eeDB（12）代入(19)式 *eeeTeeVTTeeTeeV

21、VdVBDBdVBDBdV由于节点虚位移是任意的 *eTeTeVBDB dV4 建立单元的虚功方程为 *(20)eTTeeTeeVRBDB dV由单元的虚功方程 *(20)eTTeeTeeVRBDB dV2.4 算例 任意性相等 (21)eTeeVRBD B dV节点力和节点位移之间的关系式就建立起来了 (22)eeTVKBDB dV令 (23)eeeRK则(21)式变为单元刚度方程注意：这里，节点力不是结构上的外载荷，而是按虚位移原理把单元边界上的分布力近似等效到单元节点上的一种节点力。节点力在实际结构中是不存在的总结：式（22）、（23）是由三角形常应变单元推导得到的，但是，这两式及其推导

22、过程所基于的原理和方法具有普遍性。原则上说（22）式是位移有限元分析中普遍适用的单元刚度矩阵表达式，对于不同单元，只是其中的具体计算细节不同。三角形常应变单元的刚度矩阵分析2.4 算例 0001000(10)2ijmijmiijjmmbbbBcccAcbcbcb一般情况，单元应变矩阵B是坐标的函数矩阵。 000200101011002ED这里，三角形常应变单元，B是常数矩阵。如果材料是线性的、匀质的，矩阵D也是常数矩阵。单元厚度t是常量,则dV=tdxdy,因此，三角形常应变单元的刚度矩阵可以写成： (24)eeTTTTVAAKBDB dVBDB tdxdyBDB tdxdyBDB tA=将单

23、元刚度矩阵写成分块形式： 222221111112222221111112222221124(1)eTiiiiiiijijijijimimimimiiiiiiijijijijimimimimjijijiKBDB tAbcbccbbbccbccbbbccbccbcbbccbcbbcccbbcbbcccbbb bc cb cEtA222221111222221111112222221111122222jijjjjjjjmjmjmjmjijijijijjjjjjjmjmjmjmm im im im imjmjmjmjmmc bbcb cc bb bc cb cc bc bb cc cb bb cc

24、bcbc bb cc cb bb bc cb cc bb bc cb cc bbc22212111111222222immmmm im im im imjmjmjmjmmmmmmiiijimjijjjmmimjmmcb cc bc bb cc cb bc bb cc cb bc bb ccbkkkkkkkkk（25）2.4 算例 21122, ,114(1)22rsrsrsrsrsrsrsrsrsb bc cb cc bEtkr si j mAc bb cc cb b2.4 算例 （1 1）对称性对称性单元刚度矩阵是单元刚度矩阵是对称矩阵对称矩阵；（2 2）奇异性奇异性单元刚度矩阵是单元刚度矩

25、阵是奇异矩阵奇异矩阵，它不存在逆矩，它不存在逆矩阵阵；（3 3）主元恒正主元恒正单元刚度矩阵对角元素的数值恒大于单元刚度矩阵对角元素的数值恒大于0 0，可由（可由（2525）式看出；）式看出；（4 4）单元刚度矩阵的元素具有明确的物理意义，如图所示）单元刚度矩阵的元素具有明确的物理意义，如图所示单元，其刚度矩阵第一列的单元，其刚度矩阵第一列的6 6个元素个元素k ki1i1( (i i=1,2,3,4,5,6)=1,2,3,4,5,6)的物理意义是，当单元的的物理意义是，当单元的第第1 1个节点位移个节点位移( (节点节点i i的的u ui i) )为为1 1，而而其他节点位移其他节点位移全为

26、全为0 0时，需要在时，需要在6 6个节点位移方向上施加个节点位移方向上施加的的节点力节点力的的大小大小。单元刚度矩阵的特性单元刚度矩阵的特性mvmujvjuiviuijm2.4 算例 4. 单元等效节点载荷应用虚位移原理进行载荷等效移植：移植后的节点载荷和移植前的载荷在约束允许的任意虚位移上所做的功相等。单元的节点虚位移为 *eTiijjmmuvuvuv*(16)eeuN单元的虚位移为划分单元时，一般都将作用有集中力的地方划分为节点，集中力即可直接施加到节点上。以下说明体积力和表面力向节点移植：F2.4 算例 (1)体积力等效移植体积力等效移植 令单位体积的力为： (26)Txyggg单位体

27、积力在x轴和y轴方向的分量令单元体积力等效的移植到单元节点上的等效节点载荷为： (27)eTgigxigyjgxjgymgxmgyFFFFFFF由虚位移原理， egF g和在虚位移上所做的虚功相等，即 *(28)TTeeegAFug tdxdy单元面积单元厚度(16)式代入(28)式： *TTTTeeeeeTgAAAFug tdxdyNg tdxdyNg tdxdy2.4 算例 000000ixiyTeTjxTijmgxyAAAijmjymxmyN gN gN gNNNFNg tdxdytggdxdytdxdyNNNN gN gN g(29)特殊地，体积力是重力，且重力方向为负y方向，单元的单

28、位体积力是 0TTxyggg= 00000000100113001TeTTijmgxyAAijmixiyijxAAjjymxmmyNNNFNg tdxdytggdxdyNNNN gN gNN gtdxdytdxdytANN gN gNN g (30)2.4 算例 (2)表面力等效移植表面力等效移植 2.4 算例 工程问题中，表面力一般都垂直于其作用面，所以，在有限元法中，要求定义的表面力垂直于其作用面，这样，可以将表面力分解到沿坐标轴方向。令q为表面力矢量，则它可以表示为： (31)Txyqqq表面力在x轴和y轴方向的分量令表面力等效的移植到单元节点上的等效节点载荷为： (32)eTqiqxi

29、qyjqxjqymqxmqyFFFFFFF由虚位移原理， eqF q和在虚位移上所做的虚功相等，即 *(33)eTTeeeqSFuq tds表面力作用的边界*(16)eeuN单元的虚位移为代入(33)式得到： *eeeTTTTeeeeeTqSSSFuq tdsNq tdsNq tds 000000eeeixiyTeTjxTijmqxySSSijmjymxmyN qN qN qNNNFNq tdsqqtdstdsNNNN qN qN q(34)2.4 算例 2.4 算例 xy0ijmq沿单元边界均匀分布的表面力mi边上有垂直于边界的均匀分布的表面力q,将分解为x和y方向的均布力qx和qy，这样，

30、mi边上的表面力可以表示为： 11sincos11mijxyimjyybqqllqqqqqxxcll边的边长为mil 002eTTxqxyxylyqtlFNtdsqqqqq(35)2.4 算例 5. 总体平衡方程的建立厚度为t的正方形板，左边固定，在其右上角分别作用有x、y方向的集中力F1x，F1y，建立其总体有限元平衡方程。(1)划分单元，所有节点总体编号1xF1yF121234(2)对号入座12324123( ) i( ) i( ) j( )m( )m单元号局部编码整体编码i14j32m2312( ) j单元刚度矩阵，其分块形式为2.4 算例 11111111112222222222iii

31、jimjijjjmmimjmmiiijimjijjjmmimjmmkkkkkkkkkkkkkkkkkkkk整体刚度矩阵，其分块形式为 1111112131411212221222324112122313233342224142434400iiimijmimmjjmjjmjijijmmjjjmmmiijimiikkkKKKKkkkkkkKKKKKKKKKkkkkkkKKKKkkk单元号局部编码整体编码i14j32m23122.4 算例 (3)总体平衡方程 1111111212222112122332224400iiimijmimmjjmjjmjijijmmjjjmmmiijimiikkkUFkk

32、kkkkUFUFkkkkkkUFkkk 1xF1yF121234单元号局部编码整体编码i14j32m2312 KF总体平衡方程 K整体刚度矩阵 整体节点位移 F整体节点力2.4 算例 1111111212222112122332224400iiimijmimmjjmjjmjijijmmjjjmmmiijimiikkkUFkkkkkkUFUFkkkkkkUFkkk 11111121121222211212233222344400iiimijmimmjjmjjmjijijmmjjjmmmiijimiiuvkkkUukkkkkkUvKUuUkkkkkkvUkkkuv， 111222333444xyxyxyxyFFFFFFFFFFFFF，1xF1y