. .
???2?(r?u)?2?ru?
2?rr由于各点在圆周方向上无位移，因而剪应变vr?和v?r均为
??u???r???r??u??????????零。将应变写成向量的形式，则?????????rw??
??z????z???rz????u?w???????z?r??根据上式，可推导出几何方程?????B????(e)
?zjk?N(r,z)?i?r1?0其中几何矩阵?B??2??rkj????00rkjzjkzkjNj(r,z)r0rik00rikzkizijNk(r,z)r0rji0??0??rji? zij?????3.弹性方程和弹性矩阵[D]
依照广义虎克定律，同样可以写出在轴对称中应力和应变之间的弹性方程，其形式为
1??r?u(?z???)? E1???????u(?r??z)?
E1?z???z?u(?r???)?
E2(1??)rrz??rz
E?r?所以弹性方程为?????D???? 式中应力矩阵??????r???z?rz?T
. . .
. .
??0??1?????1???0??E弹性矩阵?D?? ????1??0(1??)(1?2?)?1?2??00?0?2??4.单元刚度矩阵?k?(e)
与平面问题相同，仍用虚功原理来建立单元刚度矩阵，其积分式为
?k?(e)???B?T?D??B?dV
V
在柱面坐标系中，dV?2?drdz
将dV?2?drdz代入?k?(e)???B?T?D??B?dV，则?k?(e)?2????B?T?D??B?rdrdz
V即为轴对称问题求单元刚度矩阵的积分式。
与弹性力学平面问题的三角形单元不同，在轴对称问题中，几何矩阵[B]有的元素（如Ni(r,z)等）是坐标r、z的函
r数，不是常量。因此，乘积?B?T?D??B?不能简单地从式
?k?(e)?2????B?T?D??B?rdrdz的积分号中提出。如果对该乘积逐项求
积分，将是一个繁重的工作。一般采用近似的方法：用三角形形心的坐标值代替几何矩阵[B]的r和z的值。用?B?表示在形心(r,z)处计算出的矩阵[B]。其中
r?(ri?rj?rk)3,z?(zi?zj?zk)3
只要单元尺寸不太大，经过这样处理引起的误差也不大。被积函数又成为常数，可以提出到积分号外面：
. . .
. .
?k?(e)?2??B??D??B???rdrdz?2??B??D??B?r?式中?——三角形的面积。
TT由式?k?(e)?2?B???D??B???rdrdz?2??B??D??B?r?可以看出，两轴对
TT称的三角形单元，当形状、大小及方位完全相同而位置不同时，其刚度矩阵也不相同。距离主轴线越远的单元，其刚度越大。这与平面问题不一样。
二、等参数的刚度矩阵
对一些由曲线轮廓的复杂结构，如果采用直角边单元进行离散，由于用直线代替了曲线，除非网格划分得很细，否则不能获得较高的精度；对另一些应力随坐标急剧变化的结构，采用简单的常应力单元离散时，也必须划分成大量的微小单元，以保证足够的精度。为此引入一种高精度的单元——等参数单元。它既能简化复杂单元划分的工作，又能在满足同样精度的要求时，大大减少使用的单元数。目前流行的大程序中较常用，它成功地解决了许多二维和三维的弹性力学问题。
为导出等参数单元的刚度矩阵，首先要建立根据每个单元的形状确定的自然坐标系，然后将位移模式和形状函数都写成自然坐标的函数。
一个单元在自然坐标系的点余元整体坐标系的点成一一对应的关系。通过映射，可以将整体坐标系中的图形转化为自然坐标系中的相应徒刑。例如可以将整体坐标系中的一个任意四边形（实际单元）映射到自然坐标系中成为一个正
. . .
. .
方形（基本单元）。同样也可以将任意四面体、六面体（包括直边和曲边的）分别映射成正四面体和正六面体。
这里只介绍较简单的一种平面问题的情况，将整体坐标系中的一个任意四边形映射成自然坐标系中的一个正方体，并导出单元刚度矩阵。其它种单元的映射，可依次原理进行。不再叙述。
1. 位移模式和形状函数
图4-2中的任意四边形单元上，作连接对边中点的直线，取其交点为原点，这两条直线分别为?和?轴，并令四条边上的?和?值分别为?1，建立一个新的坐标系，称之为该单元的自然坐标系。原坐标系XOY称为整体坐标系。在整体坐标系中，自然坐标系非正交，它由任意四边形的形状所确定。
图4-19
如果将自然坐标系改画成直角坐标系，那么图4-19（a）中的任意四边形单元就成为图4-19（b）所示的正方形。
上述两个四边形的点（包括顶点）一一对应，即它们之间相互映射。因此，需要写出整体坐标X、Y和自然坐标?、?之间的坐标转换式，即
X??1??2???3???4??Y??5??6???7???8?? *
四边形四个顶点的坐标值在XOY坐标系中分别为?X1,Y1?,?X2,Y2?,?X3,Y3?,?X4,Y4?：在?o?坐标系中相应为
??1,?1?,?1,?1?,?1,1?,??1,1?。将有关数据代入*中的第一式，则有
. . .
. .
X1??1??2??3??4,X2??1??2??3??4X3??1??2??3??4,X4??1??2??3??4
求解上述方程组得：
X1?X2?X3?X4?X1?X2?X3?X4,?2?44?X1?X2?X3?X4X?X2?X3?X4?3?,?4?144?1?
坐标变换方程*成为
X?1??1????????X1??1????????X2??1????????X3??1????????X4?4同理
Y?1??1????????Y1??1????????Y2??1????????Y3??1????????Y4?4当引入函数Ni??,??后，坐标变换方程成为
X??Ni??,??Xii?144
Y??Ni??,??Yii?1式中Ni??,???1?1??i???1??i??
4变量?、?的正负号由相应节点的坐标值?i、?i决定。例如当i=4时，???1,??1，因此，N??,????1????1???。
4444下面再来研究函数Ni??,??的特性。
对节点1?X1,Y1?，相应的自然坐标值为（-1，-1）。从式
Ni??,???1?1??i???1??i??中很容易看出，除4N1=1外，N2=N3=N4=0。
对其余各节点也一样。总而言之，对节点i(i=1,2,3,4)，除Ni=1外，其余三个N值均为零。同时，不难看出
N1??,???N2??,???N3??,???N4??,???1，即四个节点的
Ni函数之
和等于1。
. . .
常用单元的刚度矩阵
