有限元分析基础理论与方法.ppt

《有限元分析基础理论与方法.ppt》由会员分享，可在线阅读，更多相关《有限元分析基础理论与方法.ppt（47页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、二有限元分析基础理论与方法山东科技大学机电学院 2012.10 李学艺第一页，编辑于星期五：十点 十二分。2弹性力学的基本理论弹性力学的基本理论(补充)1.弹性力学中的基本概念弹性力学中的基本概念 (1)体力体力是分布在物体全部体积内的力,作用在物体的每一个质点上，如重力、运动物体的惯性力和磁力等。为了表明物体内某点P所受体力的大小和方向,在这一点取物体的一小部分，它包含着P点，且体积为V。(2)面力面力是分布在物体表面上的力,例如流体压力和接触力。(3)应力物体受外力作用。或由于温度改变，其内部将发生应力。第二页，编辑于星期五：十点 十二分。3弹性力学的基本理论弹性力学的基本理论(补充)(4

2、)应变物体在受到外力和温度的作用下将发生变形。为研究物体内部一点P的变形情况,从P点处取出一个平行六面微元体开始研究。由于平行六面微元体三个棱的边长为无穷小量,所以在物体变形后,仍然是直边,但是三个边的长度和边与边之间的夹角将发生变化。各边的每单位长度的伸长或缩短量称为线应变,用表示；边与边之间的直角的改变称为切应变切应变,用 表示。第三页，编辑于星期五：十点 十二分。4弹性力学的基本理论弹性力学的基本理论(补充)(5)主应力如果过弹性体内任一点P的某一截面上的切应力等于零,则该截面上的正应力称为该点的主应力,主应力作用的这一截面称为过点P的一个应力主平面,主平面的法线方向(即主应力的方向)称

3、为P点的应力主向。(6)主应变给定应变状态下,弹性体内任意一点也存在着三个相互垂直的应变主轴,三应变主轴之间的三个直角变形后仍为直角(切应变为零),沿三个应变主轴有三个主应变,用1、2和3表示。(7)位移物体受力变形过程中,其内部各点发生的位置变化称为位移。位移为矢量。第四页，编辑于星期五：十点 十二分。5弹性力学的基本理论弹性力学的基本理论(补充)2.弹性力学中的基本方程弹性力学中的基本方程 (1)平衡微方程平衡微方程反映微分体内力与外力的关系。(2)几何方程几何方程反映弹性体位移与应变的关系。*弹性体在变形中应保持连续,即假定假定弹性体由许多微小平行六面分体构成,微小平行六面体间不应有空隙

4、,也不应重叠,其应变满足变形协调方程。(3)物理方程物理方程反映应力与应变的关系。第五页，编辑于星期五：十点 十二分。6弹性力学的基本理论弹性力学的基本理论(补充)(4)弹性力学问题的解法弹性力学问题的解法位移法位移法:取位移分量u,v,w为基本未知量，先利用位移表示的平衡微分方程和边界条件求解位移，再利用几何方程求应变，利用物理方程求应力。应力法应力法:取应力分量为基本未知量。混合法混合法：同时取部分位移分量和应力分量作为基本未知量。求解上述方程可知，弹性力学问题在数学上是由偏微分方程及其边界条件描述的,微分方程的建立和求解比较复杂,只有在物体形状和受力较简单的情况下才能获得精确解。为了避免

5、直接建立和求解这些微分方程的困难,提出了弹性力学微分方程的等价表达形式及其建立方法,即能量原理。第六页，编辑于星期五：十点 十二分。7弹性力学的基本理论弹性力学的基本理论(补充)3.弹性力学中的能量原理弹性力学中的能量原理 (1)应变能应变能弹性体在受到外载荷作用发生变形的过程中,将把克服内力所做的功作为应变能存储在弹性体内部,当外力去除后,应变能做功,使弹性体恢复原状。(2)虚位移原理虚位移原理将理论力学中刚体的虚功原理之推广到弹性体上,由于外力会使弹性体变形,因此,如果假定不存在热能和动能的改变,根据能量守恒定律,当处于平衡状态的弹性体发生体系所允许的任意微小位移时,外力在虚位移上所做的功

6、等于虚位移发生时引起的弹性体的应变能增量,称为虚位移原理。第七页，编辑于星期五：十点 十二分。8弹性力学的基本理论弹性力学的基本理论(补充)(3)最小势能原理最小势能原理最小势能原理可方便地建立弹性体基本未知量位移与外力之间的关系。一般弹性问题的最小势能原理可表述为:在满足位移边界条件的所有可能位移中,其中真实位移使系统的势能取最小值。(4)里兹法里兹法最小势能原理是利用求取积分问题的最小值将弹性力学的偏微分方程的求解化为线性方程组的求解,而积分方程极值问题的求解也不易实现。作为一种变分问题,可以利用里兹法进行求解。在1908年,里兹提出了一种“泛函变分的近似计算法”,它简化了求解过程。其基本

7、思德是假定位移函数为级数形式(富氏级数、幂级数等),且满足有关边界条件,但级数中包含一些待定的参数ai(i=1,2,)。将之代入积分方程,使其成为待定参数ai的函数。这样,将泛函求极值的变分问题转化为函数求极值的问题。第八页，编辑于星期五：十点 十二分。9弹性力学的基本理论弹性力学的基本理论(补充)4.两类平面问题两类平面问题任何一个弹性体都是空间物体,一般的外力都是空间力系,因此,任何一个实际的弹性力学问题都是空间问题，但当所研究的弹性体有特殊的形状并承受特定的载荷时，就可以把空间问题简化为近似的平面问题。这样分析和计算的工作量将大大地减少，而所得到的结果却仍然能满足工程上对精度的要求。根据

8、弹性体内的应力状态不同,平面问题可分为平面应力问题和平面应变问题两种。在有限元方法中,利用平面应力和平面应变单元来分析弹性力学的平面问题,大大简化了有限元建模的复杂性,提高了计算效率。(1)平面应力问题；(2)平面应变问题。第九页，编辑于星期五：十点 十二分。10弹性力学的基本理论弹性力学的基本理论(补充)(1)平面应力问题若物体的某一方向的尺寸较另外两个方向的尺寸小得多,如一很薄的等厚平板,仅受平行于板面的沿厚度方向均匀分布的面力,且体力也平行于板面并且不沿厚度变化,此类问题可按平面应力问题考虑，此时垂直板面方向无任何应力。(2)平面应变问题若物体的长度远大于横向尺寸，且仅受平行于横截面、沿

9、长度不变的外力作用，则按平面应变问题处理。第十页，编辑于星期五：十点 十二分。11有限元分析基础理论与方法有限元分析基础理论与方法（一）有限元分析方法的基本概念（一）有限元分析方法的基本概念有限元法分析计算的思路和步骤可归纳如下。1.分析模型的定义根据工程实际结构的形状、尺寸以及受载特点确定分析问题的力学模型。2.连续体的离散化将连续区域离散为由各种单元组成的计算模型,此过程称为有限元网格划分。离散后单元与单元之间利用节点相互连接起来，有限元法中分析的结构已不是原有的物体或结构物,而是由众多单元以一定方式连接而成的模型。这样，用有限元分析计算所获得的结果是近似的。如果划分单元类型合理且数目非常

10、多,则所获得的结果就与实际情况较接近。第十一页，编辑于星期五：十点 十二分。12有限元分析基础理论与方法有限元分析基础理论与方法3.单元特性分析(1)选择位移模式在有限元法中,一般都是选择节点位移作为未知量的所谓“位移法”进行求解。当采用位移法时,可把单元中的一些物理量如位移、应变和应力等由节点位移来表示。此时对单元中位移的分布采用一些能逼近原函数的近似函数予以描述。通常,有限元法中将位移表示为单元坐标变量的简单函数。这种函数称为位移模式或位移函数。通常假定位移函数为多项式。(2)分析单元的力学特性根据单元的材料性质、形状、尺寸、节点数目、位置及其含义等，找出单元节点力和节点位移的关系式，这是

11、单元分析中的关键一步。此时需要应用弹性力学中的几何第十二页，编辑于星期五：十点 十二分。13有限元分析基础理论与方法有限元分析基础理论与方法方程和物理方程来建立力和位移的方程式，从而导出单元刚度矩阵，这是有限元法的基本步骤之一。设节点载荷列阵用Fe表示，节点位移列阵用qe表示：则单元的力和位移的方程式为 Fe=kqe（1）式中 k单元刚度矩阵。(3)计算等效节点力物体被离散后，假定力是通过节点从一个单元传递到另一个单元。但是，对于实际的连续体，力是从单元的公共边界传递到另一个单元的。因而，需要把作用在单元边界上的表面力、体积力或集中力都等效地移到节点上，即用等效的节点力来代替所有作用在单元上的

12、力。第十三页，编辑于星期五：十点 十二分。14有限元分析基础理论与方法有限元分析基础理论与方法4.单元组装由上述获得的单元刚度矩阵和单元等效节点载荷列阵组装成表示整个结构的结构刚度矩阵和结构载荷列阵，从而建立起整个结构已知量(总节点载荷)和整个物体未知量(总节点位移)的关系式。设总刚度矩阵为、总载荷列阵为为F，整个连续体的节点位移为q，三者之间构成整体平衡方程式，即 F=kq（2）5.求解平衡方程式考虑边界条件和初始条件，求解上述联立方程组，获得节点的位移值。可根据方程的具体特点来选择合适的计算方法。第十四页，编辑于星期五：十点 十二分。15有限元分析基础理论与方法有限元分析基础理论与方法（二

13、）结构静力的有限元分析（二）结构静力的有限元分析1.单元特性推导单元特性推导有限元分析的基本步骤之一就是要导出所划分单元的刚度矩阵,一般来说,建立刚度矩阵的方法有:直接法;虚功原理法;能量变分原理方法。直接方法直接方法是直接应用物理概念来建立单元的有限元方程和分析单元特性的一种方法。该方法仅用于简单形状的单元,如梁单元。但它可以帮助理解有限元法的物理概念。虚功原理虚功原理是理论力学的一个基本原理,既可解决线形力学问题,也可用于一切非线性力学问题。最小势能原理只是虚功原理对弹性体导出的一种表述形式。能量变分原理能量变分原理方法是从势能的泛函表达式出发进行变分求极值的结果。能量变分原理方法的应用范

14、围可以方便第十五页，编辑于星期五：十点 十二分。16有限元分析基础理论与方法有限元分析基础理论与方法地扩大到机械结构位移以外的其他不含非线性的领域。如求解热传导、电磁场、流体力学等连续性问题。本课程简要介绍用虚功原理虚功原理法推导单元刚度矩阵的过程,以平面问题中的三角形单元为例进行。(1)位移函数设定对于弹性力学的平面问题，一个节点具有2个自由度,即节点有沿x轴及y轴的2个位移分量,故图1所示三角形单元有6个自由度，即6个位移分量。设三角形单元内的位移函数为:d=(u(x,y)v(x,y)T，当单元很小时,单元内一点的位移可以通过节点的位移插值来表示。可假设单元内位移为x,y的线性函数,且具有

15、6个待定系数a1a6。即第十六页，编辑于星期五：十点 十二分。17有限元分析基础理论与方法有限元分析基础理论与方法 图1 三角形单元a)三角形单元节点位移 b)三角形单元节点力第十七页，编辑于星期五：十点 十二分。18有限元分析基础理论与方法有限元分析基础理论与方法 （3）或写成矩阵形式 （3a）第十八页，编辑于星期五：十点 十二分。19有限元分析基础理论与方法有限元分析基础理论与方法在i、j、k三点有 （3b）或写成矩阵 （3c）第十九页，编辑于星期五：十点 十二分。20有限元分析基础理论与方法有限元分析基础理论与方法为了能用单元节点位移qe表示单元内某点位移d，即表达成节点位移插值函数的形

16、式，应从上式中解出a=C-1qe。其中C逆阵为 式中A三角形面积第二十页，编辑于星期五：十点 十二分。21有限元分析基础理论与方法有限元分析基础理论与方法 为使面积A不为负值，图1种的i,j,k的顺序按逆时针方向标注。将a=C-1qe=代入到式(3a)中，得第二十一页，编辑于星期五：十点 十二分。22有限元分析基础理论与方法有限元分析基础理论与方法 相乘后得(4)第二十二页，编辑于星期五：十点 十二分。23有限元分析基础理论与方法有限元分析基础理论与方法或写成(4a)可简写为d=Nqe (4b)式中 (4c)式中，Ni、Nj、Nk是坐标的连续函数，反映单元内位移分布状态，称为位移的形态函数。N

17、称为形函数矩阵。Ni、Nj、Nk为第二十三页，编辑于星期五：十点 十二分。24有限元分析基础理论与方法有限元分析基础理论与方法(4d)插值函数的性质如下：（a）插值函数在节点上的值有 (5)上式表明，在节点i上，根据式(4a)，应有u(xi,yi)=ui，其它两个形函数同样如此。（b）在单元的任一点处，三个形函数之和等于1,即第二十四页，编辑于星期五：十点 十二分。25有限元分析基础理论与方法有限元分析基础理论与方法(6)因为若单元发生刚体位移,如x方向有刚体位移u0，则单元内(包括节点上)到处应有位移u0，即ui=uj=uk=u0，由式(4b)有因此，必然要求Ni(x,y)+Nj(x,y)+

18、Nk(x,y)=1。若函数插值不满足此要求，则不能反映单元的刚体位移，用以求解则得不到正确的结果。（c）单元插值函数是线性的，在单元内部及边界上位移也是线性的，可由节点的位移唯一确定。由于公共边界上的位移完全由两个公共节点位移所确定，所以相邻单元的位移是连续的。第二十五页，编辑于星期五：十点 十二分。26有限元分析基础理论与方法有限元分析基础理论与方法为了能保证有限元法得到正确解,构造出的单元位移函数必须满足一定的条件,以使得当单元划分越来越细,网格越来越密时,所得的解答能收敛于实际问题的精确解。条件1构造出的位移函数必须在单元内连续,且在相邻单元的公共边界上协调。后者意味着变形后相邻单元之间

19、既不互相脱离(开裂),也不相互侵入(重叠)。条件2位移函数中必须包含单元的刚体位移的常数项常数项，以反映刚形位移。条件3位移函数中必须含有单元常应变单元常应变的一次项。当单元的尺寸越来越小时，单元内各点应变趋于相等。当单元尺寸无限小时,应变逼近于常量,即单元处于常应变状态。因此位移函数中就应包括一次项,因为一次项求导后为常数。第二十六页，编辑于星期五：十点 十二分。27有限元分析基础理论与方法有限元分析基础理论与方法(2)求应变由弹性力学知，或写成矩阵形式第二十七页，编辑于星期五：十点 十二分。2836(7)有限元分析基础理论与方法有限元分析基础理论与方法第二十八页，编辑于星期五：十点 十二分

20、。29有限元分析基础理论与方法有限元分析基础理论与方法(3)求应力根据虎克定律，平面问题的应力式中，对于平面问题有(8)(9)第二十九页，编辑于星期五：十点 十二分。30有限元分析基础理论与方法有限元分析基础理论与方法(4)求单元的刚度矩阵根据虚功原理，当结构受载荷作用处于平衡状态时，在任意给出的节点虚位移下，外力（节点力）Fe及内力所做的虚功之和应等于零，即(10)现给单元节点任意虚位移qe=(ui vi uj vj uk vk)T,则单元内各点将产生相应的虚位移u，v和虚应变x,y,z，它们都为坐标x,y的函数。可分别按式(4b)和式(7)求得(11)(12)第三十页，编辑于星期五：十点

21、十二分。31有限元分析基础理论与方法有限元分析基础理论与方法求单元节点力的虚功：(13)再求内力虚功：(14)(15)式中，V单元体积。将式(12)和式(8)代入到式(14)中，得到式中,qeT和qe可视为常值,将其移出积分号之外,即第三十一页，编辑于星期五：十点 十二分。32有限元分析基础理论与方法有限元分析基础理论与方法(15a)(16)(17)考虑到虚位移的任意性,可以将式(16)中两边的qeT消除,得或式中,将式(13)和式(15a)代入到虚功方程,得(17a)(17b)将前述的B和D 矩阵代入到式(17b)中，可得到平面应第三十二页，编辑于星期五：十点 十二分。33有限元分析基础理论

22、与方法有限元分析基础理论与方法(18)式中,分块阵的表达式为式中,t单元厚度。力问题三角形单元刚度矩阵的分块表达式为(18a)将式(18)代入到式(17a)中并展开得第三十三页，编辑于星期五：十点 十二分。34有限元分析基础理论与方法有限元分析基础理论与方法(19)2.总刚度矩阵的形成总刚度矩阵的形成通过单元分析得到了单元特性方程Fe=Keqe,根据此式还不能求出qe。因为Fe是单元之间的作用力,属于内力,是未知的，而通常知道的是结构所受的外力。由于内力成对出现,它们大小相等,方向相反,因此,如果将每个单元的特性方程叠加,便能消除这些成对的内力,最终只剩下已知的外力,而单元的刚度则集成为整体刚

23、度。这也就是由单元刚度形成总刚的目的。第三十四页，编辑于星期五：十点 十二分。35有限元分析基础理论与方法有限元分析基础理论与方法单元分析时已对单元的每一个节点建立了平衡方程。如式(19)中的第一式表明单元中节点i上的节点力等于该单元三个节点i、j、k的位移在节点i上所引起的节点力的叠加。这是一个单元对它的一个节点的节点力。而整体结构是由多个单元所组成的,即一个节点往往为几个单元所共有,则这个节点上的节点力就应该是共有这个节点的几个单元的所有节点位移在该节点上引起的节点力叠加总和。结构平衡时,每个节点也是平衡的。设作用在节点i上的载荷为Ri,则节点i处的平衡方程为。(20)将式(19)的第一式

24、代入到式(20)中,得第三十五页，编辑于星期五：十点 十二分。36有限元分析基础理论与方法有限元分析基础理论与方法对于结构中的所有节点，则有(21)(22)式中n节点总数。将式(22)记为式中，q=(q1，q2，qn)T，是由所有节点的位移分量组成的列阵；R=(R1，R2，Rn)T，是由所有作用在节点上的载荷组成的列阵；K是要求的总刚度矩阵，表达式为(23)第三十六页，编辑于星期五：十点 十二分。37有限元分析基础理论与方法有限元分析基础理论与方法式(23a)为整个结构的平衡方程，称为有限元方程或刚度方程。(23a)3.载荷转移载荷转移在式(23)中载荷列阵R的元素为作用在节点上的外载荷，是集

25、中力。但结构上的载荷除了集中力之外，往往还有分布的面力和体力，不可能只作用在节点上，另外，即使是集中力也不一定刚好作用在节点上，因此，需要将这些载荷转换为节点载荷,即将载荷向节点移置。载荷移置遵循能量等效原则，即原载荷与移置载荷产生的节点载荷在虚位移上所做的功相等。对于给定位移的函数，第三十七页，编辑于星期五：十点 十二分。38有限元分析基础理论与方法有限元分析基础理论与方法这种移置的结果是唯一的。在线形位移函数情况下，也可以按静力等效原则进行移置。载荷移置是在结构的局部区域内进行的。根据圣维南原理,这种移置可能在局部产生误差,但不会影响整个结构的力学特性。(1)集中力的移置集中力的移置是面力

26、和体力的基础。如图2所示,设平面单元e中某一点(x,y)作用集中力Pc设Pc移置后产生的等效节点载荷为第三十八页，编辑于星期五：十点 十二分。39有限元分析基础理论与方法有限元分析基础理论与方法如果节点发生虚位移qe,则单元内任一点位移为集中力Pc所做的虚功为第三十九页，编辑于星期五：十点 十二分。40有限元分析基础理论与方法有限元分析基础理论与方法等效节点载荷所做的虚功为根据能量等效原则,有由于虚位移是任意的，可以从上式两边同时约去，故有或由上两式可见，载荷移置的结果仅与单元形函数有关，(24)(24a)第四十页，编辑于星期五：十点 十二分。41有限元分析基础理论与方法有限元分析基础理论与方

27、法当形函数确定后，移置的结果是唯一的。(2)面力的移置设厚度为t的平面单元单位面积上作用的面力为Ps(psx,psy)T,如图3所示。并将微元面积dA=tdl上的面力PsdA视为集中力，利用式(24)并积分，可得与面力等效的移置节点载荷为 或(25)(25a)第四十一页，编辑于星期五：十点 十二分。42有限元分析基础理论与方法有限元分析基础理论与方法根据形函数的特点，在ij边上有Nk=0，所以 。因此在ij边上作用的面力只能移置到该边的两个节点上。(3)体力的移置设单元体积内作用的体力为Pv(pvx,pvy)，若将微元体tdxdy上的体力Pv tdxdy视为集中力，则利用式(24)并积分，可得

28、与体力等效的移置节点载荷为或(26)(26a)第四十二页，编辑于星期五：十点 十二分。43有限元分析基础理论与方法有限元分析基础理论与方法一般情况下,由叠加原理,一个单元上总的节点载荷应为上述三种载荷移置结果之和,即若是一个节点与多个单元相连,则节点载荷应为所有相关单元向该节点移置的载荷的叠加。因此,对于整个结构而言,有式中，ne结构的单元数量。(28)(27)第四十三页，编辑于星期五：十点 十二分。44有限元分析基础理论与方法有限元分析基础理论与方法4.约束处理约束处理在没有计入约束条件之前,结构在外力作用下,会产生刚体运动,总刚度矩阵是奇异的,因而,无法求解平衡方程获得节点位移值。故必须对

29、结构施加足够的几何约束,消除结构刚体运动,使得为非奇异阵。如果约束是位移边界条件为零,则将总刚度矩阵中的零位移分量所对应行和列的主对角元素置1,其他元素置为0。在节点载荷列阵中，将零位移分量所对应的节点载荷也变为0。如果节点位移是大于零的已知值,则将该位移分量所对应的主对角元素置为大数,再将载荷列阵中对应的分量置为大数乘以已知的节点位移,其余各元素保持不变。第四十四页，编辑于星期五：十点 十二分。45有限元分析基础理论与方法有限元分析基础理论与方法5.求解线性方程组求解线性方程组对矩阵K、R进行约束处理之后，式(23)变为式中。经约束处理后的总刚度矩阵和载荷列阵。式(29)是一个关于节点位移分

30、量的线性方程组，利用适当方法，如高斯消元法、高斯赛德尔迭代法或波前法等就可求出所有节点位移。6.计算其他物理量计算其他物理量在解式(29)求出基本未知量节点位移分量后，便可利用式(4)、式(7)式(8)求出单元上任一点的位移及单元应变和应力。(29)第四十五页，编辑于星期五：十点 十二分。46有限元分析基础理论与方法有限元分析基础理论与方法7.计算结果处理计算结果处理由于在相邻单元的公共边界上应力和应变存在不连续而发生突变，为使计算结果接近实际值,通常需要对单元应力、应变进行平均化处理。处理的方法有两种:一是绕节点平均;另一种是两单元之间平均。绕节点平均就是将环绕节点的所有单元的应力加以平均,并用该平均值表示该节点的应力。两单元平均就是将两个相邻单元的常量应力加以平均,并用该平均值表示公共边中点的应力。8.结果显示结果显示将计算得到的各种物理量以一定方式显示出来,研究结果的合理性和可靠性,评估结构性能和设计方案的优劣,做出相应的改进措施。目前数据显示方式有等值线、变形图、矢量图、二维或三维曲线以及动画等。第四十六页，编辑于星期五：十点 十二分。47思考思考如图所示，三角形单元的 边作用有均布的法向载荷q，试求其等效节点载荷 （单元厚度为t）。第四十七页，编辑于星期五：十点 十二分。