第五章平面问题有限元分析ppt课件.ppt

《第五章平面问题有限元分析ppt课件.ppt》由会员分享，可在线阅读，更多相关《第五章平面问题有限元分析ppt课件.ppt（128页珍藏版）》请在三一办公上搜索。

1、第5章平面问题的有限元分析,5-1引言,在有限元法中，把单元与单元之间设置的相互连接点，称为结点。一般用号码1、2、进行结点编号。结点可为铰接、固接或其它形式的连接。结点的设置、性质及数目等均视所研究问题的性质、描绘变形状态的需要和计算精度的要求等而定。在有限元法中引进结点概念是至关重要的。有了结点，才可将实际连续体看成是仅在结点处相互连接的单元集合组成的离散型结构，从而可使研究的对象转化成可以使用电子计算机计算的教学模型。由单元、结点、结点连线构成的集合称为有限元模型。它是有限元分析与计算的对象。,5-1-1 结构离散化,1、单元划分类型 单元类型：三角形单元、四边形单元 单元数目：根据计算

2、精度要求来确定结点设置：使单元的的结点编号尽量靠近有限元模型：由单元、结点、结点连线构成的集合,2、 离散化时应注意的问题 相邻单元的尺寸尽可能相近；同一单元的最大和最小尺寸之比尽可能接近1，不宜过大；应使全部单元中相关结点的结点编码差值为最小；单元划分时内部每个结点所连接的单元数应尽量相等,连续介质的离散,对于二维连续介质，以图所示的建筑在岩石基础上的支墩坝为例，用有限单元法进行分析的步骤如下： (1)用虚拟的直线把原介质分割成有限个三角形单元，这些直线是单元的边界，几条直线的交点称为结点。 (2)假定各单元在结点上互相铰接，结点位移是基本的未知量。 (3)选择位移函数。 (4)通过位移函数

3、，用结点位移唯一地表示单元内任一点的应变；再利用广义虎克定律，用结点位移可唯一地表示单元内任一点的应力。,(5)利用能量原理，找到与单元内部应力状态等效的结点力，再利用单元应力与结点位移的关系，建立等效结点力与结点位移的关系。 (6)将每一单元所承受的荷载，按静力等效原则移置到结点上。 (7)在每一结点建立用结点位移表示的静力平衡方程，得到一个线性方程组：解出这个方程组，求出结点位移，然后可求得每个单元的应力。,连续介质的有限单元分析包含三个基本方面：介质的离散化、单元特性计算以及单元组合体的结构分析。,3、位移函数,在选择多项式时，为了使有限单元法的计算精度和收敛性得到保障，还需要满足完备性

4、和连续性的要求。为了使位移模式尽可能地反映物体中的真实位移形态，它应满足下列条件：（1）位移模式必须能反映单元的刚体位移；（2）位移模式必须能反映单元的常量应变；（3）位移模式应尽可能地反映位移的连续性。弹性力学平面问题一般选择多项式函数作为位移函数。,5-1-2 平面问题的总势能表达式,将第二章的原理用于平面问题单元分析时，其总势能表达式为,式中相对第二章所讲的总势能表达式多了一项 ，它是单元结点力的外力势能。从整体分析角度，由于结点平衡单元结点力是不出现的，因此有的数上没有这一项，但从最小势能原理作为单元分析和整体分析的推导来说，应该包括这一项。,(5-1),5-2常应变三角形单元,三角形

5、单元是一种简单方便、对边界适应性强的单元，由于以三角形的三个顶点作为结点，因此又成为三结点三角形单元。这种单元的计算精度较低，使用的时候必须进行精细的网格划分，但他仍然是一种常用的单元,5-2-1 单元结点位移和结点力,图5-2为任一的典型单元。单元局部坐标结点编号记作1，2，3逆时针进行标记，其对应的整体结点编号记作 ，。由于平面问题局部坐标和整体坐标是一致的，因此没有坐标转换问题，故也可只标记整体编号以便形成点位向量，如图所示单元每个结点有两个位移分量，称作结点位移，记作,图5-2 常应变三角形单元,（5-2）,将三个结点位移按结点编号排在一起称为单元结点位移记作：,（5-3）,与结点位移

6、相对应，每个结点上受有2个其他单元对它作用的力，称为结点力，记作,将三个结点力按结点编号排在一起称为单元结点力记作：,（5-4）,（5-5）,单元上作用的体积力记为,（5-6）,若单元的边界是物体边界，并且该边界有表面力的话，该表面力记作：,（5-7）,体积力和表面力表达式找中矩阵元素均是沿坐标方向的分布荷载集度。,5-2-2 用面积坐标建立单元位移场,1、面积坐标的定义,设P为三角形单元中的一点，与三个顶点 ， ， 相连，则可将 分割3小块，如图5-3a所示，分别记这三个小三角形面积为,图5-3a 面积坐标示意,对于三角形单元，完全可以同杆系单元一样，在直角坐标下采用广义坐标法建立形函数及单

7、元位移场。但下面将引入的面积坐标除了也可确定点的位置外，对三角形单元的分析具有许多优越性，因此首先介绍面积坐标。,图5-3b 面积坐标示意,则显然存在如下恒等关系,等坐标轴线,有图5-3b可见，P点位置除了可用直角表示外，也可以用 ， ， 中的任意两个来确定。,等坐标轴线,等坐标轴线,（5-8）,若令,（5-9）,那么P点位置就可用量纲参数 ， ， 中的两个来确定，因此 可以作为确定点位置的一种坐标方法。因此 是根据面积来定义的，故称 ， ， 为面积坐标。,由上述定义可见:面积坐标是一种固定于单元的局部坐标，它具有如下性质： Li+Lj+Lk=1； 当P在结点 时， ；当P在结点 所对的边线上

8、时 ； 当P在与 边平行的直线上时， ， ，因此称这些直线为等坐标线，如图5-3b所示。,若以Li和Lj构成面积坐标系，则单元面积坐标系中的图形是个直角边为1的等腰直角三角形，如图5-4所示。三个结点的面积坐标分别为:,图5-4 面积坐标下的单元,结点 ：结点 ：结点 ：单元的三条边线方程为：单元形心处的坐标为,2、面积坐标与直角坐标之间的关系,单元面积 为,设P点坐标(x,y)，则有数学可知，三块小面积为,因此有,（5-10）,若将式（5-10）写成矩阵关系，则有,（5-11）,由式（5-11）可求得,（5-13）,（5-12）,即,式（5-10）和式（5-13）即为两种坐标系间的变换关系,

9、建立了坐标系间的转换关系，则按求导法则可得,（5-13）,利用数学知识 可以证明,（5-15）,（5-14）,式中 ， ， 为面积坐标的幂， 为 边的长度， 为平面问题板厚。这些公式将在今后研究三角形单元等参单元计算中应用。,3、常应变三角形单元的位移场,由形函数的性质可知，本结点为1，它结点为零，在单元内任一点全部形函数和均为1.而面积坐标的性质正好与形函数的性质相同，所以常应变三角形单元的形函数可取面积坐标，即,由此可得形函数矩阵为,其中,（5-15）,（5-18）,（5-17）,为二阶单位矩阵,有了形函数矩阵，则单元内任意一点的位移可表示为：,附、广义坐标法,位移函数,对三角形单元，假定

10、单元内的位移分量是坐标的线性函数,则 有,利用线性代数中解方程组的克来姆法则, 可解出待定常数,式中行列式:,A为ijm的面积,只要A不为0,则可由上式解出:,式中:,为了书写方便，可将上式记为:,整理上式后可得:,同理:,式中:,将三角形单元的位移函数用矩阵表示：,或,例题5-1 试求图5-5所示等腰直角三角形单元的形函数矩阵N。,图5-5 等腰直角三角形单元,应变,根据单元的位移场函数式（5-20），由几何方程可以得到单元 的应变场表达式：,（5-21）,式中微分算子矩阵A为,（5-22）,5-2-3 基于最小势能原理的单元分析,其中，矩阵称为几何矩阵。矩阵可以表示为分块矩阵的形式这里，,

11、由此位移模式可得单元中任意一点的应变矩阵,（5-25）,将式（5-19）代入式（5-24）可得,（5-24）,（5-23）,可见，由于采用了线性位移函数，应变矩阵是常数矩阵。因而单元中的应力及应变式是常数，这就是把这种单元称为常应变单元的原因。,若单元中存在初应变（由温度改变或是收缩的因素引起）,（5-27）,（5-26）,应力场,由物理方程及式（5-27），可以得到单元的应力场表达式： （5-28）其中为应力矩阵，称为弹性矩阵，对于平面应力问题，,将应力矩阵表示为分块矩阵的形式，有： （5-29）其中：,对于平面应变问题，只需将换为，换为 ，则（5-30）变为：,（5-30）,（5-31）,

12、单元刚度矩阵,（1）单元的应变能 （5-32）（2）单元上外力的势能,（5-34）,（5-33）,由式（5-32）和式（5-34）得到单元的总势能为：,利用最小势能原理，取结点位移的变分，得到：,当为内部单元时,（5-35）,（5-36）,当单元在边界处时,（5-37）,（5-38）,若记：,或,（5-39）,则式（5-36）可改写为,（5-41）,（5-40）,上式即为单元刚度方程式。,单元刚度矩阵的显式表示,对于平面应力问题，其刚度矩阵的显式：,平面应变问题,（5-42）,等效结点荷载,等效结点荷载矩阵可由式（5-40）计算获得，但是对一些简单荷载情况也不按式（5-40）计算，而是由静力等

13、效原则直接取得。因此作用在弹性体上的外力，需要移置到相应的结点上成为结点荷载。荷载移置要满足静力等效原则。,、,静力等效原则：指原荷载与结点荷载在任何虚位移上的虚功都相等。在一定的位移模式下这样的移置结果是唯一的，而且总能符合通常理解的对刚体而言的静力等效原则。,分布边界力的等效结点荷载,ij边上均布力px,ij边上三角形荷载px,分布体积力的等效结点荷载,例题5-2 试求例题5-1等腰直角三角形单元的刚度矩阵；若在单元边界点 处作用有竖向荷载 ，试求等效结点荷载矩阵 ，设 。,解：算例5-1已求出了形函数N1，N2和N3于是由式（5-24）有,当 时，平面应力问题与平面应变问题弹性矩阵D彼此

14、相等，即,由式（5-29）可得,因此可得：,又由（5-40），单元厚度为t=1:,3.2.6整体刚度矩阵,得到了单元刚度矩阵后，需要将一系列的将单元组成一个整体结构，然后根据结点载荷平衡的原则进行分析，得到整体刚度矩阵。整体分析包括以下4个步骤：（1）建立整体刚度矩阵，（2）根据支承条件修改整体刚度矩阵，（3）解方程组，求出结点的位移，（4）根据结点位移，求出单元的应变和应力。由单元刚度矩阵得到整体刚度矩阵的基本方法是刚度集成法，即整体刚度矩阵是单元刚度矩阵的集成。,整体刚度矩阵的集成,以二单元为例 进行集装,整体刚度矩阵,其中 K为二阶矩阵。,由上表可以看出，整体刚度矩阵也具有对称性。同时它

15、是一个稀疏矩阵，即其中有大量的零元素，并且非零元素都集中于主对角线附近呈带状。和单元刚度矩阵一样，由于位移函数中包含刚体位移，所以整体刚度矩阵也是一个奇异矩阵。必须要排除刚体位移后，才能变为正定矩阵。,边界条件的处理,边界的约束情况,(1)基础支承结构,(2)具有对称轴的结构,（3）具有给定位移边界的结构,边界条件的处理方法,(1)直接代入法 按结点位移已知和待定重新组合方程,对角元素改1法,只能用于给定零位移。,对角元素乘大数法,5-2-6有限元的收敛性,解的收敛性也可理解为一个问题的解的精度，较粗的分，影响一个实际问题的解的精度可分成三个方面： 实际物理问题理想化力学模型有限元求解方法（解

16、法）数字截断。 此处仅讨论解法。绪论中提到，有限元作为一种数值方法可以认为是里兹法（弹性力学解）的一种特殊形式，不同之处在于有限元法的形函数（在弹力称试探函数）是定义于单元（子域）而不是全域。里兹法的收敛条件是要求试函数具有完全性和连续性。那么它在有限元法中又是如何具体体现的？ 可从两方面： 严格的数学论证（用变分原理）； 物理方面，也就是单元模型问题（三角形单元的位移摸式）。,、位移模式的收敛性条件,（1）完备性 条件,位移函数中必须包含刚体位移；位移函数必须能反映单元的常应变状态。,（2）协调性条件,位移函数应在单元内连续，在单元之间边界上要协调。,不满协调条件的单元不一定不收敛，满足一定

17、条件的一些非协调单元不仅收敛而且收敛速度比协调单元还快、精度更高。因此，对任何单元完备性条件必须满足，如再满足协调条件的单元称为协调元。,用能量变分法可以证明，有限元位移法得到的位移解，总体上不大于真正解。即解具有下限性质。 位移解下限性质可以这样理解：所划分单元是原来连续体的一部分，具有无限多个自由度，在假定了单元的位移函数后，自由度限制为只有以结点位移表示的有限自由度，即位移函数对单元的变形进行了约束和限制，使单元的刚度较实际加强了，因此连续体刚度随之加强，故求得的近似解总体上（而不是每一点）将小于精度解。 又由此可知，杆件结构原本只在结点约束，且位移模式与假设者相同（精确方法：力法，位移

18、法），故有限元法对于杆件结构所求得的是精确解。,、位移法解的下限性质,位移法解的下限性质证明,5-3 平面矩形单元,矩形单元也是常用的单元之一，由于采用了比常应变三角形单元更高次数的位移模式，故可以更好地反映弹性体的位移状态和应力状态。,(a)直角坐标,(b)自然坐标,图5-9 矩形双线性单元,图5-9a所示为有四个结点的平面矩形单元，共有8个结点位移参数。可以直接在直角坐标下进行分析。为使分析过程简洁明了，这里引入一个量纲一的正则坐标(自然坐标),在正则坐标系下原矩形单元映射为边长为2的正方形单元如图5-9b。,5-43,5-3-1 用正则坐标建立单元位移场,广义坐标法,和杆系单元一样，根据

19、结点位移条件，可设单元的位移函数为,（5-44）,式中 为广义坐标，根据结点位移求出，代会式（5-44）整理后便可获得以形函数表达的位移场。,将四个结点的局部坐标值代入（5-44）式，可以列出八个结点位移分量方程，从而得到两组关于 的四元联立方程，即,分别解除八个位移参数，代回式（5-44），按结点位移分类合并可得：,附 形函数求解过程,其中：,其中：,写成统一表达式：,需要说明的是，位移函数中二次 而不是取平方项 或是 ，基于如下考虑：,（1）位移应该只与物体受力和变形有关，和坐标系的选取无关，即所选位移函数应该具有“坐标不变性”； （2） 如果取平方项，则沿单元边线（另一坐标为常数）位移变

20、化将是坐标的二次曲线 而此边线两相邻的单元只有两个结点具有相同的 值，有它们只能保证直线的位移变化，不可能保证两单元具有相同的二次曲线位移变化，也即不能保证相邻单元位移协调。因此，从收敛性可知，取坐标平方项是不合适的。,试凑法确定形函数,在 坐标系下，单元的四条边界线的方程分别为,根据形函数应具有的性质：本点处形函数为1，它点处形函数为0，例如,（5-45）,（a）,由式（5-45）可知下列函数且自动满足形函数它点为零的性质,（b）,代入本点坐标且令其等于1可求得,（c）,将式(c)代回式(b)且引入如下记号,则形函数可写作,可以验证广义坐标法和试凑法的结果完全相同。形函数的图形如图5-10

21、。,（5-46）,（5-47）,图5-10 双线性单元形函数 示意图,单元位移函数,有了形函数，和三角形单元一样，可由单元结点位移用下式构造单元位移场,（5-48）,式中形函数矩阵 为,（5-49）,单元协调性,从式（5-44）或式（5-47）和式（5-48）可见位移函数中包含常数项和一次项，并且从式（5-47）可验证 ，故完备性条件满足。在单元边界上 ，位移是按线性变化的，在相邻单元边界的结点位移协调的情况下，边界位移也是协调的。因此这种单元式收敛的。,5-3-2 应变和应力矩阵,采用与三角形单元完全相同的分析方法和步骤可得一下结果：,记为：,应变矩阵,这里，B矩阵称为几何矩阵。B矩阵可以表

22、示为分块矩阵的形式,（5-50）,其中,（5-51）,式中 、 为 、 方向单元边长的一半。,应力矩阵,由物理方程及式几何方程，可以得到单元的应力场表达式，,其中为应力矩阵，D称为弹性矩阵，对于平面应力问题，,将应力矩阵表示为分块矩阵的形式,其中：,对于平面应变问题，只需将E换为，换为。,（5-52）,（5-53）,5-3-3单元刚度矩阵,和三角形单元一样，可以根据最小势能原理导出结点位移向量和结点力向量之间关系，即单元的刚度矩阵 ，可以将其写成分块的形式。,其中如果单元厚度t为常数，则得到,（5-55）,（5-54）,对于平面应力问题,（5-56）,对于平面应变问题，只需将E换为，换为。,单

23、元等效结点荷载矩阵,（5-57）,式（5-57）中后一项只有单元处于边界且受表面力时才有。和三角形单元一样，一些简单情况的等效荷载也可由静力等效来得到。,5-5 平面等参单元,对于线性边界问题，矩形和三角形单元不仅精度低，而且必然产生离散误差。如果通过缩小单元尺寸来改善计算结果，又将使未知量数目剧增。为了提高精度减少未知量，本节介绍一种工程有限元分析中广泛应用的“等参数单元（简称等参元）”。,等参元的基本思想：首先建立规整形状单元（称为母单元）的形函数，然后利用它做两件事：（1）根据坐标映射用母单元形函数和实际单元的结点坐标确定所划分单元的几何形状，这个实际划分单元称为“子单元”；（2）利用母

24、单元形函数和单元结点位移建立子单元的位移场。 解决了这两个问题后，利用最小势能原理，进行一定的数学推导，就可以建立等参元的单元刚度方程。,在平面问题的有限单元法中，最简单也是常用的单元是具有三个结点的三角形单元，其次是具有四结点的矩形单元。这两种单元形状简单、规整，单元各基本矩阵如N、B、k、R等的求解比较容易，具有显式表示。并且三角形单元具有适应性强的优点，能应用于曲折的几何边界，分布不均匀的材料类型和梯度不等的应力区域，但它的精度较低。而矩形单元具有精度较高、形状规整、便于实现计算自动化等优点，但它的适应性差，不便用于曲线边界和非正交的直线边界。对于材料不均匀的结构和应力梯度不等的区域也难

25、以布置大小不等的网格。,5-5-1 基本概念,现以任意四边形等参元为例，说明等参元的有关基本概念，为进一步介绍打好基础。,1、实际单元几何形状的描述，图形变换,图5-14b为任一四结点四边形单元（子单元）。在直角坐标下要用一个解析的式子描述单元上的任意一点，也即确定其几何形状不太容易，为此选用自然坐标表示四结点正方形单元做为参照（母单元），见图5-14a。通过映射来实现从母单元到子单元之间的图形变换。其实质是建立两个单元上所有点的一一对应关系。,(a)母单元及试凑直线方程示意,(b)等参元,图5-14 四结点等参元示意,从母单元映射到任意子单元的映射关系可借子单元的结点坐标表达为,（5-58）

26、,若记单元任一点的坐标矩阵（简称坐标矩阵），结点坐标矩阵和单元结点坐标矩阵分别为,则式（5-58）可用矩阵形式表示如下,（5-59）,式中： 为母单元的形函数矩阵。,由式（5-58）或式（5-59）可得如下结论由母单元形函数的性质可知，通过映射可得到子单元结点直角坐标母单元中某一坐标线 代入式（5-58），可得它是以 为参变量的子单元整体xy坐标系中直线的参数方程。也即映射结果是子单元的一条斜直线。母单元的正交坐标轴 映射到子单元上，得到一个斜角坐标（以后的高阶单元为曲线坐标）轴，记作 轴，因此，现在子单元有两种坐标 ， 是整体坐标，而 是固定于单元的局部坐标。当母单元形函数确定后，再由各种具

27、体问题实际单元划分所确定的子单元结点坐标，由式（5-59）关系，可映射得到所有实际单元，因此，关键是建立母单元的形函数。由于母单元直线经式（5-59）映射后仍为直线，而两个相邻子单元的公共边线有相同的结点，故映射后的两子单元相邻边线处处吻合，不会产生分离和重叠。,（5-60）,2、建立单元位移场,由于 坐标既是母元的直角坐标，又是等参元的曲线（现在是斜交直线）坐标。母元的形函数是 的函数，因此，等参元的位移场亦可仿母元的位移场建立如下：,（5-61b）,（5-61a）,或,式中： 为母元的形函数； 为子元中任意一点的位移； 为子元的结点位移。,3、几点说明,等参元的位移分布与母元的位移分布即使

28、在结点位移相同的情况下也可能是不同的。单元位移场与单元几何形状都是用相同的形函数构造的，并用相同个数的结点参数来描述，因此称这类单元为等参数单元。若由结点坐标插值构造单元几何形状所用的形函数比结点位移插值构造单元位移场的形函数阶次低，并且所用的结点参数个数少，则称为亚参元。反之，若阶次高，结点个数多，则称为超参元 。,4、两类坐标系之间导数的变换关系,形函数矩阵N只是局部坐标,的显函数,为求形函数对整体坐标x,y的偏导数，必须用复合函数求导公式:,或写成：,（5-62）,（5-63）,将式（5-58）代入式（5-63）则有,由式（5-64）可得,（5-64）,（5-65）,式中,称为雅可比矩阵

29、。而,（5-66）,称为雅可比行列式。基于式（5-65）和 式（5-66）对于平面问题的微分算子矩阵为,（5-67）,式中,（5-68）,（5-69）,雅可比矩阵计算示例,单元a,单元b,单元c,1、八结点曲边四边形等参元,5-5-2 几种常用单元描述和位移模式,图5-15给出了八结点曲边四边形等参元母元与子元示意图，四边形等参元四结点曾加到八结点，一方面可以进一步提高精度，另一面又可拟合曲边界线，因而此单元用途比较广。,母元的形函数,八结点四边形等参元母元是八结点正方形单元（图5-15a），除四个角顶结点外，还取四边中点作为新结点。根据形函数性质，可直接由试凑法得到八个形函数。下面介绍由低阶

30、母元形函数按统一方法生成高阶母元形函数的步骤。,(a)母单元及试凑法方程示意图,(b)等参元,图5-15 八结点等参元示意,(1)根据形函数性质，用试凑法建立边界上新增结点函数。对八结点等参元为,（2）计算低阶母元函数在新增结点处的值，对八结点等参元分别为：,（3）对低阶母元角顶结点形函数进行“它点为零的修正”。例如八结点等参元的N1为：,不难验证直接试凑和统一方法所得的结果完全相同。八结点等单元的全部形函数为：,（5-70）,或是写成如下形式,（5-71）,四边形十二结点等参元 参考教材119页【例题5-7】,位移场,按照等参元建立位移场思路，八结点四边形等参元的位移模式可由八结点母元形函数

31、构造得到:,或,（5-73a）,（5-73b）,补充：拉氏族单元（矩形单元）,拉格朗日插值：,例：边长为2,其他形函数自行推出。,结点的等效荷载,2、六结点三角形等参元,图5-17给出了用面积坐标表示的六结点三角形母元及六结点三角形等参元示意。,(a)母单元及试凑法方程示意图,(b)等参元,图5-15 六结点三角形等参元示意,母元的形函数,在面积坐标下，由直接试凑或统一方法可得六结点三角形等参元形函数为,（5-74b）,（5-74a）,角结点,边中点,三角形十结点等参元参考120页【例题5-8】。,图形变换,根据等参元定义，子单元形状可由下式确定,或,式中 又子单元的六个结点坐标分量所组成。,

32、（5-75）,位移场,根据等参元定义，六结点三角形等参元的位移模式，可由母元形函数构造而成。即,或,（5-76）,补充：索氏单元（三角形单元）,索氏插值：函数构造就是选择能盖着其他点的直线的乘积。,例：,二次单元,其他形函数自行推出。,三次或三次以上单元,例：,其他形函数自行推出。,5-5-3 等参元单元特性分析,建立了上述三个单元的位移场后，可按统一的单元分析步骤进行推导，得出各自的单元刚度矩阵。,1、分析的基础,单元应变场,式中,（5-77）,（5-78）,（5-79）,其中,（5-80）,单元应力场,单元总势能,将位移、应变和应力代入最小势能原理表达式，可得,（5-81）,（5-82）,

33、2、微面积及微弧长的计算,为由积分得到单元刚度方程，尚需解决规则坐标下dA,dl 映射后的计算问题。,dA的计算,图5-19 微面积示意,（5-83）,dl的计算,如图所示 界面,因此在此界上,（5-84b）,（5-84a）,对一般情况,（5-85）,3、单元刚度特性公式,将式（5-83）和式（5-84）代入式（5-82）经整理可得,（5-86）,式中,（5-87）,对正方形母元,式中,（5-88）,（5-89）,或由结点体积力，表面力的值插值构成。,三角形类母元,形函数为,若令,（5-90）,（5-91）,则可将形函数转换成 的二元函数。因此,（5-92）,除此之外还需修改积分限，从图5-2

34、1可见,图5-21 积分限变化示意,（5-93）,对于边界表力的等效荷载的微分弧长应视不同边界而异。,（5-94a）,（5-94b）,不论对 还是对 积分，其积分限均为0到1。,参考例题教材124页【例题5-9】,5-5-4 数值积分,在求解刚度矩阵和结点荷载时，需计算如 的积分。但一般是很复杂的，通常难以用显式表示其积分，一般都用数值积分方法计算积分值，即在单元内选出某些点，称为积分点，求出被积函数 在这些点的值，然后根据这些数值求出积分值。数值积分有两类方法，一类方法积分点是等间距的，如辛普生方法；另一类方法积分点是不等间距的，如高斯方法。,一维高斯积分公式,和 是根据计算精度最高而选定的，积分点 应是勒让德多项式 的根。加权系数 按下式计算。,方形域,二维及三维高斯积分公式,先令 保持常数，计算沿 方向的积分,再沿 方向积分,对三重积分有,三角形域,采用面积坐标L1，L2，L3，在三角形域上的高斯积分公式可写出,这里n积分点个数，L1i，L2i，L3i是积分点i的面积坐标，Hi是加数系数，可按教材127页表5-2选取。,附 例题,结束,The end,