有限元方法 第五章 平面三角形单元.ppt

第五章 平面三角形单元 本章内容结束 例如，图3-11(a)所示受纯弯曲的梁，其结构对于x、y轴都是几何对称的，而所受的载荷则是对于x轴对称，对于x轴反对称。可知，梁的应力和变形也将具有同样的对称特性，所以只需取1/4梁进行计算即可。取分离体如图3-11(b)所示，对于其它部分结构对此分离体的影响，可以作相应的处理，即对处于y轴对称面内各节点的x方向位移都设置为零，而对于在x轴反对称面上的各节点的x方向位移也都设置为零。这些条件就等价于在图3-11(b)中相应节点位置处施加约束，图中o点y方向施加的约束是为了消除刚体位移。 　　节点的布置是与单元的划分互相联系的。通常，集中载荷的作用点、分布载荷强度的突变点，分布载荷与自由边界的分界点、支承点等都应该取为节点。并且，当物体是由不 二. 节点的选择及单元的划分 图3-11 　　节点的多少及其分布的疏密程度（即单元的大小），一般要根据所要求的计算精度等方面来综合考虑。从计算结果的精度上讲，当然是单元越小越好，但计算所需要的时间也要大大增加。另外，在微机上进行有限元分析时，还要考虑计算机的容量。因此，在保证计算精度的前提下，应力求采用较少的单元。为了减少单元 (a) (b) 图3-12 同的材料组成时，厚度不同或材料不同的部分，也应该划分为不同的单元。 　　在进行节点编号时，应该注意要尽量使同一单元的相邻节点的号码差尽可能地小，以便最大限度地缩小刚度矩阵的带宽，节省存储、提高计算效率。如前所述，平面问题的半带宽为 B =2 (d+1) ，在划分单元时，对于应力变化梯度较大的部位单元可小一些，而在应力变化比较平缓的区域可以划分得粗一些。 　　还有一点值得注意的是，单元各边的长度不要相差太大，以免出现过大的计算误差或出现病态矩阵。例如，图3-12所示的(a)、(b)两种单元划分，虽然都是同样的四个节点，但(a)的划分方式显然要比(b)的方式好。 三. 节点的编号 例如，对图5-3所示的单元 jm和ijn ，具有公共边ij。 这样，不论按哪个单元来计算，根据（5-11）式，公共边ij上的位移均由下式表示 图5-3 由(5-23)式可知，在ij边上 式中 Ni , Nj 的表达形式如（5-23）式所示。 (i) 利用形函数的这一性质可以证明，相邻单元的位移分别进行线性插值之后，在其公共边上将是连续的。 由此可见，在公共边上的位移 u、v 将完全由公共边上的两个 节点i、j 的位移所确定，因而 相邻单元的位移是保持连续的。 为了在以后讨论问题中能够比较 方便地确定单元中任意一点处的 形函数数值，这里引入面积坐标的概念。 在图5-4所示的三角形单元ijm中， 任意一点P(x , y)的位置可 以用 以下三个比值来确定 图5-4 式中 ?为三角形单元ijm的面积，?i 、?j 、?m 分别是三角形Pjm、Pmi、Pij的面积。这三个比值就叫做P点的面积坐标。 (5-24) 显然这三个面积坐标并不是完全独立的，由于 所以有： 而三角形pjm的面积为： 故有： 类似地有 (5-25) (5-26) 由此可见，前述的三角形常应变单元中的形函数Ni 、Nj 、Nm 就是面积坐标Li 、Lj 、Lm 。 根据面积坐标的定义，我们不难发现，在平行jm边的直线上的所有各点，都有相同的坐标Li ，并且该坐标就等于“该直线至jm边的距离”与“节点i至jm边的距离”之比，图5-4中给出了Li 的一些等值线。 容易看出，单元三个节点的面积坐标分别为 节点 i： Li =1 Lj =0 Lm =0 节点 j： Li =0 Lj =1 Lm =0 节点m： Li =0 Lj =0 Lm =1 不难验证，面积坐标与直角坐标之间存在以下变换关系： (5-27) 当面积坐标的函数对直角坐标求导时，可利用下列公式： (5-28) 一. 单元刚度矩阵 为了推导单元的节点力和节点位移之间的关系，可应用虚位移原理对图5-2中的单元e进行分析。单元e是在等效节点力的作用下处于平衡的，而这种节点力可采用列阵表示为 (a) 假设在单元e中发生有虚位移，则相应的三个节点i、j、m 的虚位移为 且假设单元内各点的虚位移为{f *}，并具有与真实位移相同的位移模式。 §5-4 刚度矩阵 故有 (c) 参照（5-13）式，单元内的虚应变{? *}为 于是，作用在单元体上的外力在虚位移上所做的功可写为 (d) (f) 而单元内的