内容发布更新时间 : 2024/6/7 16:07:50星期一 下面是文章的全部内容请认真阅读。

武汉理工大学有限元分析基础教案

?u?a1?a2x?a3y?a4xy ??v?a5?a6x?a7y?a8xy由于我们选择的是多项式插值。所以对于2元函数多项式项元的选择，可以用帕斯卡三角形确定。如左图。

从该三角形可以看出，2元多项式的二次项有三项，即x2、y2、xy。如果选取前面两项中的任何一项，都会造成位移模式偏惠与你所选择的那个方向。从而使位移出现不对称的情况。只有增加xy这一项，才能避免这种情况出现。也就满足了我们材料各向同性的要求。

选取局部坐标系，将四个节点的坐标代入该公式，并求解待定系数，最后可以得出和前面三角形单元类似的关系：

?Ni??????0其中的Ni?0NiNj00NjNm00NpNm0?u??ui?0??i???????N?????

Np???v??v??p??p?1?x??y?1?x??y?1?1?N?1?1? ???????? j4?a??b?4?a??b?1?x??y?1?x??y??1???1?? Nm??1???1?? 4?a??b?4?a??b? Np?在上面公式中，代表节点坐标的是什么？

位移模式的讨论。 1，应变和刚体位移 由几何方程得

?x??y??xy??u???a1?a2x?a3y?a4xy??a2?a4y ?x?x?v???a5?a6x?a7y?a8xy??a2?a8x?y?y?v?u?????a5?a6x?a7y?a8xy???a1?a2x?a3y?a4xy??a3?a4x?a6?a8y ?x?y?x?y同三角形单元位移模式一样，常应变和刚体位移包括在位移模式之中。同时我们可以看到，应变中还有随x、y的变量应变。所以在矩形单元不再是常应变或常应力单元了。

2，在局部坐标系中，当x=±a（y=±b）时，位移模式是x（或y）的线性函数（将x=±a（y=±b）代入位移表达式中即可）。这表明在边界上位移按线性规律变化，在公共边界上的位移是连续的，满足收敛性的充分条件。 二，矩形单元的刚度矩阵

??x??????y?????xy?e???e?????x??0??????y??0?0b?y0b?y0??b?y?0??ui????b?y?e????u?1???0??a?x?0??a?x?0a?x0a?x??????????v?y??4ab???a?xb?ya?x??b?y???vm????a?x???b?y???a?x?b?y?????x??应力向量

???e??D??B???e?。分别代入平面应力和平面应变的弹性矩阵，就可以得出两类问题的应

武汉理工大学教务处制 46

武汉理工大学有限元分析基础教案

力解。

刚度矩阵?K??eT相乘后逐项积分即得单元刚度??B??A?D??B?tdxdy代入[B]和[D]矩阵，

矩阵，见教材中的具体计算公式。需要说明的是，由于[B]不再是常数矩阵（还有x、y），所

以积分运算较三角形单元复杂些。具体见如下矩阵。

矩形单元的载荷向节点移置的方法同三角形单元总刚矩阵组集的原则和方法、边界条件

武汉理工大学教务处制 47

武汉理工大学有限元分析基础教案

的处理等也同三角形单元。其中不同的位置仅是单元节点位移和节点载荷的数目，不再是6个而是8个。 三，讨论

矩形单元的优缺点：

1）其位移模式推导出的应力、应变不再是常量，分布更接近实际物体中的分布。从而使这种单元的计算精度更高。

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?来确

AjAiALi? Lj? Lm?m

AAAA——三角形的面积。

Ai——是三角形pjm的面积。 Aj——是三角形pim的面积。 Am——是三角形pij的面积。

当P在图中虚线上任一点是，Li是相同的（三角形面积为底×高的一半），在该三角形的三

武汉理工大学教务处制 48

武汉理工大学有限元分析基础教案

个顶点，分别有： i，Li=1、Lm=Lj=0 j，Lj=1、Lm=Li=0 m，Lm=1、Li=Lj=0

而且有Li?Lj?Lm?1，将节点的面积坐标和前面我们推导出的形函数Ni进行比较，可以知道，形函数实际就是面积坐标。 三，六节点三角形单元的位移模式

?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??N?2LL?? ii?i所以

2??同理可得Nj?2Lj?Lj???1?1??? Nm?2Lm?Lm??

2?2??对于N1应满足j、m节点处为零，所以构造函数?LjLm， 令?LjLm=1（在节点1处应为1），可以求得β=4 所以N1?4LjLm，同理得N2?4LiLm

N3?4LiLj

综合以上推导，六节点三角形单元的形函数如下：

?1??N?2LL??ii?i?2??? 1、2、3、i、j、m循环 N1?4LjLm??1?ai?bix?ciy??Li?2A?四，单元刚度矩阵

武汉理工大学教务处制 49

武汉理工大学有限元分析基础教案

推导过程类似于三角形单元的刚度矩阵，不同之处在于在对形函数求偏导时，要运用复合函数求导的方法（因为插值基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???? ii?ii2??x2Ai???????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°，试分析其齿的受力状况。

武汉理工大学教务处制 50