第五章.等参数单元

母单元 首先，根据形函数的定义，在局部坐标中，建立起几何形 状简单且规整的单元，我们称之为母单元。
1. 一维母单元 采用局部坐标ξ，单元为直线段，即。具体形式如下： 1） 线性单元（2结点）
1 2 1 2

1 -1 0 (a) 线 性 单 元
2 1
N1
N2

2） 二次单元（3结点）
(8-14)
其中， N是用局部坐标表示的形函数，(x,y)是结点i 的整体坐标，上式即为平面坐标变换公式。
返回
图5-4表示了一维单元的坐标变换。原来的直线状的母单 元分别变换成了直线、二次曲线和三次曲线状的子单元，这是 因为变换式中的形函数Ni分别是ξ的一次、二次和三次函数。
y

3 1 2
1 -1
(8-19)
返回
其中，[J]-1是[J]的逆阵
y 1 J x
3. 三维母单元 三维母单元是坐标系中的2×2×2正六面体
1 1 1 1 1 1
如图5-3所示，坐标原点在单元形心上，单元边界是六个平面。 单元结点在角点及各边的等分点上。 1) 线性单元（8结点） 5 8
13
5
16
15 14
8

6
1
这正方形单元的位移模式是：
而其中形函数为：
由图（b）可知
• 假如图 (a)中的任意四边形单元能用上式的位移 模式及形函数进行计算，则前面所提的位移连续 性条件就可以得到满足，所以问题归结为：如何 将任意四边形单元的整体坐标(x，y)，变换成正 方形单元的局部坐标( , )。
根据形函数的两条性质：
2
图5-2 以上形函数也可以合并表示为 1 1
Ni
0 0
二维母单元
（i=1, 2, 3, 4）
4
其中
0 i
0 i

2) 二次单元（8结点）
角点：
Ni 1 4
4
7
3 6
8
o

1 1
0 0
0
0 1
如果有任意四边形单元，如图(a)所示就可以克服矩形单元 之不足，但是这种单元的位移模式如何能否满足前面所述的 条体则是本节要解决的问题。
在图（a)中的任意四边形单元上，作连接对边中点的直线，称 之为 及 ，取其交点为原点，并令四边上的坐标值分别为 1，就得出一新坐标系，称之为单元的局部坐标系。 将局部坐标系改画成直角坐标系，则图 (a)中的任意四边形单元 就变成图 (b)所示的正方形单元。
16
N4
9 1 1 16
图5-1 一维母单元
2. 二维母单元 二维母单元是平面中的2×2正方形
1 1 1 1
如图5-2所示，坐标原点在单位形心上。单元边界是四
条直线：
1
，
1 。为保证用形函数定义的未知量在相
邻单元之间的连续性，单元结点数目应与形函数阶次相适应 。因此，对于线性、二次和三次形函数，单元每边的结点数
N1

N3 1
2
1
2
N2
1
2
1 -1
2 0
3 1

3） 三次单元（4结点）
N1 N3
(b) 二 次 单 元
1 9
16 9 1
2
2
1
N2
1 9
16
2
1
1 3
2
0
3 1

0
x
(a) 线 性 单 元
(b) 二 次 单 元
图5-4 一维单元的平面坐标变换
返回
图5-5表示了二维单元的平面坐标变换。母单元是正方形， 子单元则分别变换成任意四边形和曲边四边形。而且相邻子单 元在公共边上的整体坐标是连续的。以二次单元为例，两个相 邻单公共边界上都是二次曲线（抛物线），而在三个公共结点 上具有相同的坐标。因此，整个公共边界都有相同的坐标，即 相邻单元是连续的。
返回
而曲线坐标系 , , 则只适用于单个独立的子单元，所以称 为局部坐标。整体坐标在整体分析中采用，局部坐标则在单 元分析中采用。 现在讨论两类坐标系中有关偏导数的关系，以二维坐标 为例：根据复合函数的求导法则，有


x x x x
i 1,
2
i 1;
0 0
1 1 1
(8-12)
3) 三次单元（32结点） 角点： 典型边中点：
Ni
1 1 1 9 64
1
0 0 0
2

2

2
19
(8-13)
i
Ni 9 64
从图形变换的角度看， , , 和 x , y , z 可以分别看成是母单 元和子单元这两个不同单元的坐标系，它们都是直角坐标系。 而从另一角度看， , , 和 x , y , z 又可以看成是同一单元（子单 元）的两种不同的坐标系。 x , y , z 是子单元的直角坐标系，而 , , 可看成是子单元的曲线坐标系。可以看出 x , y , z 始终扮演 同一角色，即子单元的直角坐标；而 , , 则扮演两种角色， 它既是母单元的直角坐标，又是子单元的曲线坐标。 在有限元分析中，两者的作用是不同的。直角坐标系在 x , y , z 整个结构的所有子单元中共同采用，所以称为整体坐标。

y y y y
(8-16)
上式可写成矩阵形式

x J y
(8-17)
返回
其中：[J]称为雅可比（Jacobi）矩阵
第五章 等参数单元
第一节 位移模式和形函数
第二节 等参元的概念 第三节 平面等参元 第三节 空间等参元
返回
第一节
位移模式和形函数
一、位移模式 在前面几章中已阐明位移模式就是：单元内任意 一点的位移，被表述为其坐标的函数。在平面问题 的单元中，任一点的位移分量可用下列多项式表示：
为了使有限元的解能够收敛于精确解，任何单 元的位移模式都必须满足以下三个条件： （1）位移模式中必须包括反应刚体位移的常数 项。 （2）位移模式中必须包括反应常应变的线性位 移项。 （3）位移模式必须能保证单元之间位移的连续 性。
分别为两个、三个和四个。除四个交点外，其他结点位于各
边的二分点或三分点上。
返回
1) 线性单元（4结点）
1 1
4

4
3
N1 N3
N2 N4
1 1
4
o

1 1
4
1 1
4
1 (a) 线 性 单 元
7 4
6
17 18 9
7 20
19
1
2
3
(a) 线 性 单 元பைடு நூலகம்
2 三维母单元

10
12
4 3
11
(b) 二 次 单 元
图5-3
返回
2) 二次单元（20结点） 角点： 典型边中点：
i 0 ,
Ni 1 4
Ni
1 8
1 1 1
0 0 0
0
0 0 2
显然在四个节点处，上式所示的关系无疑是成立的现在 要证明在四条边上，这关系也是正确的。以1—2边为例， 在此边上局部坐标 ＝-1，代入，得
或改写成：
同样可得：
等号左边 (x,y)是整体坐标，等号右边 ( )是局部坐标，因此上式被称为坐标 , 变换式 。
通常称局部坐标的正方形单元为母单元或基本单元，称整体 坐标的任意四边形单元为子单元或实际单元。

y
4
7
3
6
1
3 6 1

8
o

4
1
7 8 1
1
5
2
0
5 1
2
（a）母单元
（b）子单元
x
返回
图5-5 二维单元的平面坐标变换
=1 5
13
16
8
15
y
5 13 16 15 17 14 7 9 1 12 19 10 3 8
6 =-1
18 17
14
=-1 7 20
1
（i=1, 2, 3, 4） 边中点：
Ni Ni 1 2 1 2
2 5 (b) 二 次 单 元
1 1
2 0
（i = 5, 6）
(8-9)
1 1
2 0
（i = 7, 8）
3) 三次单元（12结点） 角点：
Ni 1 32
1 1 9
描述位移和描述坐标都采用相同形函数 ，所以这种单元称为等 参数单元。
1. 平面坐标变换 在整体坐标系中，子单元内任一点的坐标用形函数表示如下
N , x N , y
i i
x y
i i
N 1 , x 1 N 2 , x 2 N 1 , y 1 N 2 , y 2
1 3
,
2
i 1,
i 1;
1 1 9 1 1
0 0 0
返回
第二节
等参数元的概念
在平面问题的有限元中，最简单的单元是三节点的三角形单 元，由于这种单元中的应变及应力是常数，而通常计算对象 的应力场又往往随坐标而急剧变化，所以在应用常应变的三 角形单元时，必须划分大量的微小单元，才能得到较好的计 算精度，因而用三节点三角形单元算题时，往往节点数最多， 原始输入数据庞大。四节点的矩形单元能够比三角形单元更 好的反映实际应力变化，但它不能适府曲线边界和非直角的 直线边界，也不便随意改变大小。所以上述的两种单元都有 其不足之处。