有限元分析教案
矩形单元的载荷向节点移置的方法同三角形单元总刚矩阵组集的原则和方法、边界条件的处理等也同三角形单元。其中不同的位置仅是单元节点位移和节点载荷的数目，不再是6个而是8个。
三，讨论
矩形单元的优缺点：
1）其位移模式推导出的应力、应变不再是
46
有限元分析教案
常量，分布更接近实际物体中的分布。从而使这种单元的计算精度更高。
2）对于斜边性边界和曲线边界的拟合性差，且不便于在不同部位采用不同大小的单元。实际应用中常需要三角形单元在大小不同的矩形单元之间进行过渡。见下图。
§2.17 六节点的三角形单元
从以上对矩形单元的推导，可以看出，增加单元节点的数目可以提高单元位移模式的插值次数，进而可以提高单元计算精度。沿用这样的思路，我们对三角形单元也可以采用相同的办法。不过增加的节点是在每条边的中点。这样就得到了一个六节点的三角形单元。如图所示。
一，位移模式
由于这样的单元有6个节点，所以节点位移向量的个数就是2×6=12个，那么在位移模式中可以有12个待定系数。参看帕斯卡三角形，我们选用以下的位移模式：
?u?a1?a2x?a3y?a4x2?a5xy?a6y2 ?22?v?a7?a8x?a9y?a10x?a11xy?a12y将每个单元公共边上的方程代入，可以求得一个二次的抛物线方程，而公共边上有三个公共节点，所以可以唯一确定这条抛物线。这表明该插值模式满足位移连续的收敛性充分条件。
仍然可以按照前面三节点三角形单元的方法，分别将6个节点位移代入，然后解联立方程组，求12个待定系数。这样的计算，由于待定系数过多，计算过程也过于繁杂。所以实际中是采用面积坐标的形式进行计算。 二，面积坐标
三角形三个顶点i、j、m，p为三角形中任意一点，
其在三角形中的位置，可以用Li定。其中：
?LjLm?来确
AjAiAL?Li? j Lm?m
AAAA——三角形的面积。
Ai——是三角形pjm的面积。 Aj——是三角形pim的面积。 Am——是三角形pij的面积。
当P在图中虚线上任一点是，Li是相同的（三角形面积为底×高的一半），在该三角形的三个顶点，分别有： i，Li=1、Lm=Lj=0 j，Lj=1、Lm=Li=0 m，Lm=1、Li=Lj=0
而且有Li?Lj?Lm?1，将节点的面积坐标和前面我们推导出的形函数Ni进行比较，可以知道，形函数实际就是面积坐标。 三，六节点三角形单元的位移模式
47
有限元分析教案
?u?Niui?Njuj?Nmum?N1u1?N2u2?N3u3? v?Nv?Nv?Nv?Nv?Nv?Nviijjmm112233?推求如下：
Ni在i节点处应为1，在j、3节点处应为零，面积坐标Li虽在i节点处为1，但在2、3节点处却为1/2，所以还是用Li作为Ni就不行了。为此我们构造一个函数?Li?Li?该函数是否满足形函数的性质。 由于Li在2、3节点处，Li???1??，考察2?1所以很显然上式等于零。 2??1??=1，可以计算出来） 2?在i节点处要为1，那么β=2。（令?Li?Li?1??? 所以Ni?2Li?Li?2??同理可得Nj?2Lj?Lj???1?1??N?2LL? ? ?mm?m2?2??对于N1应满足j、m节点处为零，所以构造函数?LjLm， 令?LjLm=1（在节点1处应为1），可以求得β=4 所以N1?4LjLm，同理得N2?4LiLm
N3?4LiLj
综合以上推导，六节点三角形单元的形函数如下：
?1??N?2LL???iii?2??? 1、2、3、i、j、m循环 N1?4LjLm??1?ai?bix?ciy??Li?2A?四，单元刚度矩阵
推导过程类似于三角形单元的刚度矩阵，不同之处在于在对形函数求偏导时，要运用复合函数求导的方法（因为插值基N是面积坐标的函数，面积坐标才是x、y的函数） 如：?x??u???Niui?Njuj?Nmum?N1u1?N2u2?N3u3? ?x?x?Ni?Ni?Li????x?Li?x?Li??1????1????2LL??a?bx?cy???? iiiii??2???x?2A???48
有限元分析教案
?N1?N1?Lj?N1?Lm ???x?Lj?x?Lm?x在求得[B][D]矩阵后，运用
?K?e???A?B?T?D??B?tdxdy就可以求得单刚矩阵。但是由于
[B]矩阵的复杂性，所以求解[K]e也不再是很容易的一项工作了。必须运用一些相关的积分公式。具体的计算公式可以见教材的相关内容。 五，讨论
在单元数目相同的情况下，六节点的三角形单元计算精度远比三节点三角形单元高，也比矩形单元高。但由于一个节点的相关节点数目，在6节点的三角形单元中大大增加，所以总刚度矩阵的的带宽也较三节点三角形单元宽的多。所需的内存也相应增加。
从理论上说来，运用这种增加节点，改善单元计算精度的方法，可以不断的运用下去。如4节点矩形单元可以变为8节点矩形单元等等，但是实际运用中，我们可以看到，节点数目的增加，导致单元计算复杂性的大大增加，有时可能求解不出单刚计算公式（在现有的计算分析理论上）。所以在有限元程序中，再高次数的单元运用就比较少了。6节点三角形单元为二次单元，而三节点三角形单元为一次单元。
§2.18 平面问题的计算实例
例1 一直齿齿轮，齿数为20，模数是3，齿厚20mm，压力角为20°，试分析其齿的受力状况。
49
有限元分析教案
例2 变厚度圆筒的压力容器，受有内压，试用有限元法计算圆筒内外壁的应力。
50