弹性力学第6章--用有限单元法求平面问题课件.ppt

1、第五节第五节 单元的结点力列阵与劲度矩阵单元的结点力列阵与劲度矩阵第四节第四节 单元的应变列阵和应力列阵单元的应变列阵和应力列阵 第三节第三节 单元的位移模式与解答的收敛性单元的位移模式与解答的收敛性第二节第二节 有限单元法的概念有限单元法的概念第一节第一节 基本量及基本方程的矩阵表示基本量及基本方程的矩阵表示概述概述第六节第六节 荷载向结点移置荷载向结点移置 单元的结点荷载列阵单元的结点荷载列阵例题例题第十一节第十一节 应用变分原理导出有限单元法的基本方程应用变分原理导出有限单元法的基本方程第十节第十节 计算实例计算实例第九节第九节 计算成果的整理计算成果的整理第八节第八节 解题的具体步骤解

2、题的具体步骤 单元的划分单元的划分第七节第七节 结构的整体分析结点平衡方程组结构的整体分析结点平衡方程组1、有限元法有限元法（Finite Element Method）FEM2、FEM的特点的特点 概述概述(1)(1)具有具有通用性和灵活性通用性和灵活性。首先将连续体变换为离散化结构，然后再利用分片插值技术与虚功原首先将连续体变换为离散化结构，然后再利用分片插值技术与虚功原理或变分方法进行求解。理或变分方法进行求解。简称简称FEM，是弹性力学的一种是弹性力学的一种近似解法。近似解法。(2)(2)对同一类问题对同一类问题,可以编制出可以编制出通用程序通用程序,应用计算机进行应用计算机进行计算。

3、计算。(3)(3)只要适当加密网格，就可以达到工程要求的精度。只要适当加密网格，就可以达到工程要求的精度。3、FEM简史简史 19431943年柯朗年柯朗(CurrantCurrant)第一次提出了第一次提出了FEMFEM的概念。的概念。FEMFEM是上世纪中期才出现，并得到迅速发展和广泛应用的一是上世纪中期才出现，并得到迅速发展和广泛应用的一种数值解法。种数值解法。简史19561956年年,特纳特纳(Tunner)等人提出了等人提出了FEMFEM。上世纪上世纪5050年代年代,平面问题的平面问题的FEMFEM建立建立,并应用于工程问题。并应用于工程问题。19601960年提出了年提出了FEM

4、FEM的名称。的名称。上世纪上世纪6060年代后年代后,FEM,FEM应用于各种力学问题和非线性问题应用于各种力学问题和非线性问题,并得到迅速发展。并得到迅速发展。3、FEM简史简史 19431943年柯朗年柯朗(Currant)第一次提出了第一次提出了FEMFEM的概念。的概念。FEMFEM是上世纪中期才出现，并得到迅速发展和广泛应用的一是上世纪中期才出现，并得到迅速发展和广泛应用的一种数值解法。种数值解法。上世纪上世纪7070年代后年代后,FEM,FEM被引入我国被引入我国,并很快地得到应用和发展。并很快地得到应用和发展。5 5、本章介绍平面问题的本章介绍平面问题的FEMFEM4、FEMF

5、EM的主要导出方法的主要导出方法 应用静力方法或变分方法导出。应用静力方法或变分方法导出。仅叙述按位移求解的方法。仅叙述按位移求解的方法。且一般都以平面应力问题来表示。且一般都以平面应力问题来表示。6-1 基本量和基本方程的矩阵表示基本量和基本方程的矩阵表示 本章无特别指明本章无特别指明,均表示为均表示为平面应力问题平面应力问题的公式。的公式。采用采用矩阵表示矩阵表示,可使公式统一、简洁可使公式统一、简洁,且便于编制程序。且便于编制程序。()Txyfff(,),(,)Tu x y v x yd()Txyxy。Txyyx)()Tiijju v uv()TixiyjxjyFFFFF1 1、基本物理

6、量的矩阵表示、基本物理量的矩阵表示()Txyfff体力体力:位移函数位移函数:应变应变:应力应力:结点位移列阵结点位移列阵:结点力列阵结点力列阵:面力面力:基本物理量(2)物理方程物理方程:(b)D21010(c)1001/2ED 2、FEM中应用的方程中应用的方程 (a)Tuvuvxyxy(1)几何方程几何方程:应用的方程其中,D为弹性矩阵,对于平面应力问题是:(3)虚功方程虚功方程:()()TTAdxdyt*F为结点虚位移及对应的虚应变。*,*其中，在FEM中用结点的平衡方程代替平衡微分方程,后者不再列出。(3)(3)整体分析整体分析。6 6-2-2 有限单元法的概念有限单元法的概念 FE

7、MFEM的概念的概念，可以简述为：，可以简述为：FEM的概念(1)(1)将连续体变换为离散化结构将连续体变换为离散化结构(结构的离散化结构的离散化)；(2)(2)单元分析；单元分析；FEMFEM的分析过程：的分析过程：该方法的理论基础是该方法的理论基础是分片插值技术分片插值技术与与变分原理变分原理。采用有限自由度的离散单元组合体模型去描述实际具有采用有限自由度的离散单元组合体模型去描述实际具有无限自由度的考察体，是一种在力学模型上进行近似的数值无限自由度的考察体，是一种在力学模型上进行近似的数值计算方法。计算方法。结构力学研究的对象是离散化结构。如桁架，各单元(杆件)之间除结点铰结外，没有其他

8、联系(图(a)。弹力研究的对象，是连续体(图(b)。结构离散化1.结构离散化结构离散化FEMFEM的分析过程的分析过程(1)(1)将连续体变换为离散化结构(图(c)：即将连续体划分为有限多个、有限大小的单元，并使这些单元仅在一些结点处用绞连结起来，构成所谓“离散化结构”。(c)深 梁（离 散 化 结 构）图(c)与图(a)相比,两者都是离散化结构；区别是，桁架的单元是杆件，而图(c)的单元是三角形块体（注意:三角形单元内部仍是连续体）。结构离散化 比如：将深梁划分为许多三角形单元，这些单元仅在角点用铰连接起来。(c)深 梁（离 散 化 结 构）将连续体变换为离散化结构(图(c)：即将连续体划分

9、为有限多个、有限大小的单元，并使这些单元仅在一些结点处用绞连结起来，构成所谓“离散化结构”。1.结构离散化结构离散化FEMFEM的分析过程的分析过程(1)(1)2.2.单元分析单元分析 求解方法 每个三角形单元仍然假定为连续的、均匀的、各向同性的完全弹性体。因单元内部仍是连续体，应按弹性力学方法进行分析。取各结点位移 为基本未知量，然后对每个单元,分别求出各物理量,并均用 来表示。()(1,2,)Tiiiu viFEMFEM的分析过程的分析过程(2)(2)(1,2,)ii(1)应用插值公式,由单元结点位移 ,求单元的位移函数Tmjie)(,),(,)Tu x y v x yd.该插值公式称为单

10、元的位移模式位移模式,记为ed.单元分析的主要内容：单元分析的主要内容：(2)应用几何方程,由单元的位移函数d,求出单元的应变单元的应变eB.求解方法(4)应用虚功方程,由单元的应力 ,求出单元的结点力单元的结点力表示为(3)应用物理方程,由单元的应变 ,求出单元的应力单元的应力.eS(eeijmFF FFk.其中，为结点对单元的作用力，作用于单元，称为结点力，以正标向为正。(TiixiyFFF2.2.单元分析单元分析FEMFEM的分析过程的分析过程(2)(2)(1)应用插值公式,由单元结点位移 ,求单元的位移函数Tmjie)(,),(,)Tux y vx yd.该插值公式称为单元的位移模式位

11、移模式,记为ed.单元分析的主要内容：单元分析的主要内容：(2)应用几何方程,由单元的位移函数d,求出单元的应变单元的应变eB.(5)将每一单元中的各种外荷载,按虚功等效原则移置到结点上,化为结点荷载 (.eeLLiLjLmFFFF 为已知值,是用结点位移表示的值。各单位移置到i 结点上的结点荷载iF求解方法LiF3.3.整体分析整体分析,iF,FLi),2,1(,ieLieiFFe各单元对i 结点的结点力作用于结点i上的力有：FEMFEM的分析过程的分析过程(3)(3)其中,表示对围绕i 结点的单元求和；通过求解联立方程，得出各结点位移值，从而求出各单元的应变和应力。求解方法 3.3.整体分

12、析整体分析 2.2.对单元进行分析对单元进行分析 1.1.将连续体变换为离散化结构将连续体变换为离散化结构归纳起来，归纳起来，FEMFEM分析的主要步骤分析的主要步骤：（1 1）单元的位移模式）单元的位移模式（2 2）单元的应变列阵）单元的应变列阵（4 4）单元的结点力列阵单元的结点力列阵（5 5）单元的等效结点荷载列阵）单元的等效结点荷载列阵建立结点平衡方程组，求解各结点的位移。建立结点平衡方程组，求解各结点的位移。（3 3）单元的应力列阵）单元的应力列阵思考题1.有限单元法求解问题的基本步骤是什么？2.试说明单元分析的主要内容。复复 习习1 1、基本物理量与基本方程的矩阵表示、基本物理量与

13、基本方程的矩阵表示(2)物理方程物理方程:Tuvuvxyxy(1)几何方程几何方程:(3)虚功方程虚功方程:D()()TTAdxdyt*F 3.3.整体分析整体分析 2.2.对单元进行分析对单元进行分析 1.1.将连续体变换为离散化结构将连续体变换为离散化结构2.2.FEM分析的主要步骤：分析的主要步骤：位移模式位移模式 应变列阵应变列阵结点力列阵结点力列阵 等效结点荷载列阵等效结点荷载列阵应力列阵应力列阵 应用插值公式,可由 求出位移 。首先首先,必须解决由单元的结点位移 来求出单元的位移函数 FEM是取结点位移 为基本未知数的。问题是如何求应变、应力。该插值公式表示了单元中位移的分布形式,

14、因此称为位移模式。Tmjie(i(,)(,)Tu x y v x yd。e6-3 单元的位移模式与解答的收敛性单元的位移模式与解答的收敛性 位移模式d插值公式(a)在结点 应等于结点位移值 。由此可求出 在结点三角形单元中,可以假定位移分量只是坐标的线性函数,也就是假定：123456,auxyvxy。),(,mjiyxii,(,)iiu v i j m。61三角形单元 其中 包含。及,iiiivuyx61将式(a)按未知数 归纳为:,iivu,biijjmmi ijjm muNuN uN uvNvN vN v。三角形单元或用矩阵表示为:123456,auxyvxy。000 000iiijmji

15、jmjmmuvNNNuuNNNvvuv edN.c()2,(,)iiiiNa bx cyAi j m11,(,)11jjjjiiimmmmxyyxabci j mxyyx N 称为形函数矩阵形函数矩阵,其非零元素为其中，A为ijm的面积(图示坐标系中,i,j,m按逆时针编号),有：imjxyoi11121iijjmmxyAxyxy。三结点三角形单元的位移模式，略去了2次以上的项，因而其误差量级是 且其中只包含了x,y的1次项，所以在单元中Ni的分布如图(a)所示,u,v的分布如图(b)、(c)所示。);(2xo 三角形单元 所以当单元趋于很小时，即 时，为了使FEM之解逼近于真解。则为了保证F

16、EM收敛性,位移模式应满足下列条件：FEM中以后的一系列工作，都是以位移模式为基础的。0,yx收敛性条件 所以当单元趋于很小时，即 时，为了使FEM之解逼近于真解。则为了保证FEM收敛性,位移模式应满足下列条件：0,yx 因为当单元尺寸趋于0时，单元中的位移和应变都趋近于基本量刚体位移和常量位移。(1)位移模式必须能反映单元的刚体位移。(2)位移模式必须能反映单元的常量应变。xxyvyyxu22,22353564353521,00 xvvyuu可见刚体位移项在式(a)中均已反映。而刚体位移形式(P17(2-9)式)为，将式(a)写成对式(a)求应变,得：2635,.xyxy 可见常量应变也已反

17、映。(3)位移模式应尽可能反映位移的连续性。ij,收敛性条件即应尽可能反映原连续体的位移连续性。在三角形单元内部，位移为连续；在两单元边界ij 上，之间均为线性变化，也为连续。收敛性条件 所以当单元趋于很小时，即 时，为了使FEM之解逼近于真解。则为了保证FEM收敛性,位移模式应满足下列条件：0,yx(1)位移模式必须能反映单元的刚体位移。(2)位移模式必须能反映单元的常量应变。(1)和(2)是必要条件，而加上(3)就为充分条件。为了保证为了保证FEM的收敛性：的收敛性：imjxyoi思考题1.应用泰勒级数公式来选取位移模式，为什么必须从低次项开始选取？2.试考虑：将结构力学解法引入到求解连续

18、体的问题时，位移模式的建立是一个关键性工作，它使得单元(连续体)内部的分析工作都有可能进行了。6-4 6-4 单元的应变列阵和应力列阵单元的应变列阵和应力列阵 应用几何方程，求出单元的应变列阵：。mmjjiimmjjiivNvNvNvuNuNuNu,()/2(,)iiiiNabxc yAi j m,位移函数其中，单元中的位移函数单元中的位移函数用位移模式表示为()0001000 2TiiijmjijmjiijjmmmmuvvuxyxyuvbbbucccvAcbcbcbuveB.a应变,应力应用几何方程，求出单元的应变列阵：()0001000 2Tiiijmjijmjiijjmmmmuvvuxy

19、xyuvbbbucccvAcbcbcbuveB.a)(),(bmjiBBBB)(),(0021cmjibccbAiiii。iB其中,B 称为应变矩阵应变矩阵，用分块矩阵表示，,(d)eeDDBS再应用物理方程，求出单元的应力列阵：应变,应力其中,S称为应力转换矩阵应力转换矩阵，写成分块形式为,(d)eeDDBS 对于线性位移模式，求导后得到的应变和应力，均成为常量，因此，称为常应变（应力）单元。应变和应力的误差量级是 ,其精度比位移低一阶，且相邻单元的应力是跳跃式的。(),(e)ijmSS SS2(,)(f)2(1)1122iiiiiiiibcEbci j mAcbSDB.()ox再应用物理方

20、程，求出单元的应力列阵：思考题 如果在位移模式中取到泰勒级数中的二如果在位移模式中取到泰勒级数中的二次幂项，略去次幂项，略去 高阶小量，试考虑位移、高阶小量，试考虑位移、应变和应力的误差量级。应变和应力的误差量级。3x6-5 6-5 单元的结点力列阵与劲度矩阵单元的结点力列阵与劲度矩阵 现在来考虑其中一个单元：模型,结点力 在FEM中，首先将连续体变换为离散化结构的模型。(2)单元与周围的单元在边界上已没有联系,只在结点i,j,m互相联系。(1)将作用于单元上的各种外荷载,按静力等效原则移置到结点上去,化为等效结点荷载。故单元内已没有外荷载。假想将单元与结点i 切开，则：),(,)(mjiFF

21、TiyixiF),(,)(mjiFFTiyixiF其数值与 相同，而方向相反。iF以沿正坐标向为正。对单元而言,这是作用在其上的“外力外力”。(1)结点作用于单元上的力，称为结点力，(2)单元作用于结点的力，为：;)(TmjieFFFF().Txyxy 而其内部有应力作用，考察已与结点切开后的单元i,j,m，则此单元上作用有外力,即结点力,应用虚功方程，求单元的结点力：假设发生一组结点虚位移 则单元内任一点(x,y)的虚位移为 单元内任一点(x,y)的虚应变为 代入虚功方程：在单元中，外力(结点力 )在虚位移(结点虚位移)上的虚功，等于应力 在虚应变 上的虚功，即：,)(e*,)(e*Nd,)

22、(e*B eF)()(*e)(*()().ae TTAdxdyt*e*F虚功方程,)()()(TTeTeTBB*代入其中 与x,y无关，故式(a)成为e)(*()()e Te TAdxdyt*e*TFB.因为 是独立的任意的虚位移，虚功方程对任意的 均应满足，故e)(*e)(*式式(b)是是由应力求结点力的一般公式由应力求结点力的一般公式。TAdxdyteFB.b式(c)是由结点位移求结点力的一般公式，k k称为单元的劲度劲度矩阵矩阵.其中,再将应力公式代入上式，得单元劲度矩阵 ceTeeAdxdytFB DBk dTAdxdytkB DB式(b)是由应力求结点力的一般公式。TAdxdyteF

23、B.b对于三角形单元，B B矩阵内均为常数，有 etATkB DB代入B B，D D，得出k k如书中(6-37)及(6-38)所示。(6-37)iiijimjijjjmmimjmmkkkkkkkkkk22 21122=,114(1)22 ,;,.(6-38)rsrsrsrslnrsrsrsrsrsrsb bc cb cc bEtkAc bb cc cb br si j m l nx yk对于三角形单元，B B矩阵内均为常数，有 etATkB DB代入B B，D D，即得平面应力问题中三结点三角形单元的刚(劲)度矩阵，可写成如下分块矩阵分块矩阵的形式：其中，(1)k k是66的方阵,k k中元

24、素 表示仅在单元结点s沿n方向产生单位位移时引起结点r沿l方向的结点力。(2)由反力互等定理，所以k k是对称矩阵，以对角线为对称轴。,kkrssr单元劲度矩阵单元劲度矩阵k k的性质的性质：(3)当单元作刚体平移时，如 三角形内不产生应力和应变，结点力也为0。1,ijmuuu(4)由(3)可导出行列式 (即k为为奇异矩阵奇异矩阵)。(5)k k的元素与 单元的形状和方位等有关，但与单元的大小和刚体的平动以及作 度转动无关。即,k k中每一行（或列）元素之和为0（其中第1、3、5元素之和(对应x向)或2、4、6元素之和(对应y向)也为0）。,tE0knlnrsk例题例题 某等腰直角三角形单元i

25、jm如图所示,已知在所选取的坐标系中,单元结点坐标分别为:0000ijmijmxaxxyyay 应用11,(,)11jjiimmyxbci j myx 11121iijjmmxyAxyxy。可得200 .2ijmijmbabbaaccacaA 应用公式可得该单元的应力转换矩阵应力转换矩阵为()ijmSS SS2(,)2(1)1122iiiiiiiibcEbci j mAcbSDB.2(1)1001 0011 (e)1111002222EaS.可得该单元的应力转换矩阵应力转换矩阵为210010011 (e)(1)1111002222EaS.应用教材式(6-37)及式(6-38)可得该单元的单元刚

26、度矩阵单元刚度矩阵为2110211022 (f)0012(1)1131222111312222Etk对称 现考察结点力与单元中的应力之间的关系。为了简单起见，假定只有结点i发生位移ui，如右图(a)所示。由上面的单元刚度单元刚度矩阵矩阵得相应的结点力为：210012(1)1001,TixiyjxjymxmyTiTFFFFFFEtuF其中，。相应的结点位移及结点力如图所示22(1)iEtuF 另一方面，由于发生了位移ui,则根据上面得到的应力转换矩阵，可得应力分量为210(1)210,TxyxyTiTEuaFta如上图(b)中单元的两直角面所示。根据单元的平衡条件,还可得出斜面上的应力,如图所示

27、。若将这三个面上的应力分别按静力等效原则积铢累寸到结点上去，可以得到图(a)中相同的结点力。思考题思考题试求出书中例题(P117)的位移模式。()()eTTLixLiyLjxLjyLmxLmyFFFFFFLLiLjLmFFFF .6 6-6-6 荷载向结点移置荷载向结点移置,单元的结点荷载列阵单元的结点荷载列阵 在FEM中,须将作用于单元上的外荷载向结点移置，化为等效结点荷载等效结点荷载，(2)变形体静力等效原则变形体静力等效原则:即在任意的虚位移上，使原荷载与移置荷载的虚功相等。1、等效原则、等效原则 刚体静力等效原则只从运动效应来考虑,得出移置荷载不是唯一的解;变形体的静力等效原则考虑了变

28、形效应,在一定的位移模式下,其结果是唯一的,且满足了前者条件。所以在FEM中，采用变形体的静力等效原则。(1)刚体静力等效原则刚体静力等效原则:使原荷载与移置荷载的主矢量以及对同一点的主矩也相同。移置原则 假设发生一组结点虚位移 ,则点的虚位移为 使移置荷载的虚功等于原荷载的虚功：原荷载 作用于单元中任一点(x,y),为单位厚度上的作用力;移置荷载 作用于结点ijm。,)(TPyPxffPf(),TeLLiLjLmFFFF().e*dN()()().e TeTe Ttt*TLPPFdfN fe)(*集中力,面力,体力2 2、集中力的移置公式集中力的移置公式 对于任意的虚位移 ,虚功方程都必须满

29、足,得：aetTLPFN f.e)(*3、单元边界单元边界 上面力上面力 的移置公式的移置公式 应用式(a),将 代之为 并在边界 上积分,得:SftPf,dstfs beSdstTLFN f.应用式(a)，将 代之为 并对单元域A 积分,得 ceAdxdytTLFN f,dxdytftPf4 4、单元内体力、单元内体力f 的移置公式的移置公式 说明：说明：当位移模式当位移模式为线性函数时，由虚功为线性函数时，由虚功方程得出的移置荷载，方程得出的移置荷载，与按刚体静力等效原则与按刚体静力等效原则得出的结点荷载相同。得出的结点荷载相同。思考题1.试导出书中例题(P119)的荷载移置公式。在单元分

30、析中，从单元的结点位移求位移分布求应变求应力求结点力，为单元的内力分析；外荷载移置到结点荷载，为单元的外力分析。6 6-7-7 结构的整体分析结构的整体分析 结点平衡方程组结点平衡方程组 iFLiF 假设将结点i与周围的单元切开，则围绕i结点的每个单元对i 结点有结点力()的作用,也有外荷载移置的结点荷载()的作用。下面考虑整体分析整体分析。对某一个单元ijm，其中 是对围绕i 结点的单元求和。结点结点i 的平衡条件的平衡条件为 ,(1,2,)aiLieeinFFe,iinnni j mFk 结点平衡条件其中,ijm是单元结点的局部编号;i=1,2,n是其整体编号。代入式(a)，可表示为 ,(

31、)(1,2,)binnLieni j meinkF。将式(b)按整体结点编号排列,得整个结构的平衡方程组。对某一个单元ijm，其中 是对围绕i 结点的单元求和。结点结点i 的平衡条件的平衡条件为 ,(1,2,)aiLieeinFFe,iinnni j mFk .cLKF 其中，分别为整体结点位移列阵，整体结点荷载列阵和整体劲度矩阵。12(),(),TT12nLLLLnFFFFK 考虑结构的约束条件后,从式(c)求出 ,就可以求出各单元的位移和应力。结点平衡方程组例例2例例1 列出图示结构i 结点的平衡条件。（见书中P.121）,()innLien i j me k F1、有限单元法的具体计算步

32、骤、有限单元法的具体计算步骤6 6-8-8 解题的具体步骤解题的具体步骤 单元的划分单元的划分 1、划分单元网格，对单元和结点编号。2、选定直角坐标系，按程序要求填写和输入有关信息。单元内的ijm的局部编号应按书中规定的右手规则编号。否则会使三角形的面积出现负号等问题。3、使用已编好的程序进行上机计算。事先须将有限单元法的公式，计算方法和步骤都编入程序。4、对成果进行整理、分析整理、分析。对第1步和第4步的工作，也尽可能由计算机完成，以减少人工的工作量。如自动划分网格，整理成果等。2 2、单元划分注意事项、单元划分注意事项（8）结构具有凹槽或孔洞等应力集中处等。（1）单元大小问题；（2）单元在

33、不同部位的合理布置问题；（3）三角形三个内角最好较接近；（4）利用对称性和反对称性；（5）厚度突变之处和材料不同之处；（6）载荷作用（集中力或突变分布载荷）处；（7）水利闸坝工程问题；在FEM中，位移的精度较高，其误差量级是 ，即与单元尺度的二次幂成正比。应力的误差量级是，即与单元的大小成正比。6 6-9-9 计算成果的整理计算成果的整理 )(2xo()ox 三结点三角形单元的应力的成果，不但应力的精度较低，而且还产生了所谓应力的波动性应力的波动性。对于结点位移的成果，不需整理就可以直接采用。应力的波动性在三结点三角形单元中较为显著。这是由于计算出的应力的精度较低。假设单元的应力成果为 ，其中

34、 为真解，为误差。则由于在结点都列出了平衡方程并令其满足，从而使相邻的单元的应力趋近于 。这就产生了应力的波动性。为提高应力精度，解决其波动性问题，可采取以下应力成果整理方法：1)两相邻单元平均法;2)绕结点平均法。一般地讲，两相邻单元平均法的精度较好，因为它涉及的区域范围较小。在面力边界线附近，求得的应力误差较大。可采用向外插值的方法（例抛物线插值）来解决。为了提高应力精度，可采用两种方法：一是加密网格，减少单元的尺寸，以提高应力的精度。二是可以采用较多结点的单元，并使位移模式中包含一些高幂次的项，从而提高位移和应力的精度。应力的波动性在三结点三角形单元中较为显著。这是由于计算出的应力的精度

35、较低。假设单元的应力成果为 ，其中 为真解，为误差。则由于在结点都列出了平衡方程并令其满足，从而使相邻的单元的应力趋近于 。这就产生了应力的波动性。为提高应力精度，解决其波动性问题，可采取以下应力成果整理方法：1)两相邻单元平均法;2)绕结点平均法。一般地讲，两相邻单元平均法的精度较好，因为它涉及的区域范围较小。书中应用三结点三角形单元，计算了下列例题：6 6-10-10 计算实例计算实例 1.楔形体受自重及齐顶水压力。2.简支梁受均布荷载。3.圆孔附近的应力集中。在整理应力成果时,读者应注意,应用三角形单元时，(1)采用两单元平均法和绕结点平均法两单元平均法和绕结点平均法的应力成果比较接近，

36、但前者的精度略好于后者。(2)边界面的应力,宜采用向外插值向外插值的方法求出。在FEMFEM中，将连续体变换为离散化结构之后，有两种导出FEM公式的主要方法，即静力法与变分原理法。6 6-11-11 应用变分原理导出有限单元法基本方程应用变分原理导出有限单元法基本方程 (2)建立单元位移模式,求出单元中的位移分布，beB.ceSi aedN1.按静力方法导出按静力方法导出FEMFEM公式公式(1)取结点位移为基本未知数；(3)由几何方程求出单元的应变，(4)由物理方程求出单元的应力，eesAtdstdxdyteTTLPFN fN fN f fLKF.deeFkeLF(5)由虚功方程求出单元的结

37、点力，(6)由虚功方程求出单元的结点荷载 ，(7)建立结点平衡方程组，按结构力学方法导出FEM公式（1）变分原理中的极小势能原理是。minVUEP）（g)().(,21hdxdytdsttVdxdytUATsTTAfdfdfdPT2.按变分方法导出按变分方法导出FEMFEM公式公式 保留上述（1）-（4）步骤，然后应用极小势能原理导出FEM基本方程。按变分法导出FEM公式对于平面问题，对于连续体,变分的宗量是位移函数 变分方程 可表示为总势能 对 的导数等于0，即PEvu,.,vu）（g,0uEP。0vEP）（i变分宗量由 变换成（2）将经典变分原理应用到离散化结构，则).,2,1(nii,e

38、ePPEE,eeUU。eeVV）（jvu,总势能、形变势能和外力势能，可以用单元的势能之和来表示其中 为三角形单元的面积。应用前面记号,)()(21)(2121 eAeTTeeAeTeeAeeeeedxdytdxdytdxdytUUDBBDB TeA。eeTeUk)(21）（k内力势能为其中 为三角形单元的受面力边界。引用前面记号,)()()(eATSPTPTTeeATSPTPTeeeedxdytdsttdxdytdsttVVfdfdfNfdfdfds。eeLTeVF)(）（leeLTeTePVUE)()(21Fk）（m外力势能为 总势能为故总势能极小值条件 变换为（3）对于离散化结构，泛函数

39、 的宗量变 换为 PE)(i)(),2,1(0nniEiP。.)(,2)(ccaababaaaTT),2,1(nii则式(n)成为引用矩阵运算公式，)(),2,1(0)(oniEieTeP。,001,)(ieeLmLjLimjieeLeePEFFFFFFFk 其中 代入式(o)，得出与结构力学方法导出的相同方程，)(),2,1(pnieLiei。FF 从物理意义上讲，将连续体的经典变分原 理(g)或(i)应用到离散化结构，成为式(p)。比较物理意义：凡是与微分方程对应的变分原理存在的任何问题，均可应用变分法导出FEM。式(p)表示总势能在所有结点处的极值条件。式(g)表示总势能的整体极值条件;作业作业(06)P141习题6-2,6-4,6-6第四章例题本章结束本章结束