常用单元的刚度矩阵

.实用文档.

. rurrur22)(2

由于各点在圆周方向上无位移，因而剪应变rv和rv均为零。将应变写成向量的形式，那么rwzuzwrururzzr

根据上式，可推导出几何方程)(eB

其中几何矩阵ijjikiikjkkjjiikkjkjiijkjjkzrzrzrrrrrzrNrzrNrzrNzzzB0000),(0),(0),(00021

3.弹性方程和弹性矩阵[D]

依照广义虎克定律，同样可以写出在轴对称中应力和应变之间的弹性方程，其形式为

)(1zrruE

)(1zruE

)(1rzzuE

rzrzEr)1(2

所以弹性方程为D

式中应力矩阵Trzzr .实用文档.

. 弹性矩阵221000010101)21)(1(ED

4.单元刚度矩阵)(ek

与平面问题相同，仍用虚功原理来建立单元刚度矩阵，其积分式为

dVBDBkVTe)(

在柱面坐标系中，drdzdV2

将drdzdV2代入dVBDBkVTe)(，那么rdrdzBDBkTe2)(

即为轴对称问题求单元刚度矩阵的积分式。

与弹性力学平面问题的三角形单元不同，在轴对称问题中，几何矩阵[B]内有的元素〔如rzrNi),(等〕是坐标r、z的函数，不是常量。因此，乘积BDBT不能简单地从式rdrdzBDBkTe2)(的积分号中提出。如果对该乘积逐项求积分，将是一个繁重的工作。一般采用近似的方法：用三角形形心的坐标值代替几何矩阵[B]内的r和z的值。用B表示在形心),(zr处计算出的矩阵[B]。其中

3)(,3)(kjikjizzzzrrrr

只要单元尺寸不太大，经过这样处理引起的误差也不大。被积函数又成为常数，可以提出到积分号外面：.实用文档.

. rBDBrdrdzBDBkTTe22)(式中——三角形的面积。

由式rBDBrdrdzBDBkTTe22)(可以看出，两轴对称的三角形单元，当形状、大小及方位完全相同而位置不同时，其刚度矩阵也不相同。距离主轴线越远的单元，其刚度越大。这与平面问题不一样。

二、等参数的刚度矩阵

对一些由曲线轮廓的复杂结构，如果采用直角边单元进行离散，由于用直线代替了曲线，除非网格划分得很细，否那么不能获得较高的精度；对另一些应力随坐标急剧变化的结构，采用简单的常应力单元离散时，也必须划分成大量的微小单元，以保证足够的精度。为此引入一种高精度的单元——等参数单元。它既能简化复杂单元划分的工作，又能在满足同样精度的要求时，大大减少使用的单元数。目前流行的大程序中较常用，它成功地解决了许多二维和三维的弹性力学问题。

为导出等参数单元的刚度矩阵，首先要建立根据每个单元的形状确定的自然坐标系，然后将位移模式和形状函数都写成自然坐标的函数。

一个单元在自然坐标系内的点余元整体坐标系内的点成一一对应的关系。通过映射，可以将整体坐标系中的图形转化为自然坐标系中的相应徒刑。例如可以将整体坐标系中的一个任意四边形〔实际单元〕映射到自然坐标系中成为一.实用文档.

. 个正方形〔根本单元〕。同样也可以将任意四面体、六面体〔包括直边和曲边的〕分别映射成正四面体和正六面体。

这里只介绍较简单的一种平面问题的情况，将整体坐标系中的一个任意四边形映射成自然坐标系中的一个正方体，并导出单元刚度矩阵。其它种单元的映射，可依次原理进行。不再表达。

1. 位移模式和形状函数

图4-2中的任意四边形单元上，作连接对边中点的直线，取其交点为原点，这两条直线分别为和轴，并令四条边上的和值分别为1，建立一个新的坐标系，称之为该单元的自然坐标系。原坐标系XOY称为整体坐标系。在整体坐标系中，自然坐标系非正交，它由任意四边形的形状所确定。

图4-19

如果将自然坐标系改画成直角坐标系，那么图4-19〔a〕中的任意四边形单元就成为图4-19〔b〕所示的正方形。

上述两个四边形的点〔包括顶点〕一一对应，即它们之间相互映射。因此，需要写出整体坐标X、Y和自然坐标、之间的坐标转换式，即

87654321YX *

四边形四个顶点的坐标值在XOY坐标系中分别为44332211,,,,,,,YXYXYXYX：在o坐标系中相应为1,1,1,1,1,1,1,1。将有关数据代入*中的第一式，那么有 .实用文档.

. 43214432134321243211,,XXXX

求解上述方程组得：

4,44,443214432134321243211XXXXXXXXXXXXXXXX

坐标变换方程*成为

4321111141XXXXX同理

4321111141YYYYY当引入函数,iN后，坐标变换方程成为

iiiiiiYNYXNX,,4141

式中iiiN1141,

变量、的正负号由相应节点的坐标值ii、决定。例如当i=4时，1,144，因此，411,4N。

下面再来研究函数,iN的特性。

对节点11,1YX，相应的自然坐标值为〔-1，-1〕。从式iiiN1141,中很容易看出，除N1=1外，N2=N3=N4=0。对其余各节点也一样。总而言之，对节点i(i=1,2,3,4)，除Ni=1外，其余三个N值均为零。同时，不难看出

1,,,,4321NNNN，即四个节点的Ni函数之.实用文档.

. 和等于1。

函数,iN具备上章所介绍的形状函数应满足的条件，可作为本单元的形状函数。

采用,iN做形状函数，其位移模式为

iiiiiivNvuNu,,,4141

比照iiiiiiYNYXNX,,4141和iiiiiivNvuNu,,,4141可以看出：在这种实际单元〔任意四边形〕中，坐标变换式和位移模式不仅采用了相同的形状函数,iN，而且具有相同的数学模型。这种性质的实际单元称为等参数单元。

对用节点位移值ui〔或vi等〕求单元内某一点位移量u(或v等)的插值公式iiiuNu,41，只要将u(或v等)换成X(或Y等)，便成为利用节点值Xi(或Yi等)求相应点坐标X(或Y等)的插值公式。相反也是这样。

2.几何矩阵[B]

由于几何矩阵[B]通过对位移求偏导数而得出，所以首先必须利用复合函数求导的规那么得出下述公式

yuXuJuuYYuXXuYYuXXuuu或写成 .实用文档.

. 式中YXYXJ，此式称为雅可比矩阵。

为了将几何矩阵[B]写成变量、的函数，必须将式yuXuJuu改写成

uuJyuXu1，同理，vvJyvXv1

从表示单元内各点位移与其应变关系的几何方程可知：

TYvXvYuuvuXYYX01101000000100

将式uuJyuXu1和vvJyvXv1合并，那么

TTvvuuJJYvXvYuXu1100

对单元〔e〕，任意一点的位移u，v对自然坐标，的偏导数可利用上式求出，写成矩阵形式为：

epTNvvuu

式中Tevuvuvuvu11111111

pppppppppNNNNNNNNN4321432100000000 .实用文档.

. 对于i=1,2,3,4 TiiipNNN

将TTvvuuJJYvXvYuXu1100和epTNvvuu代入TYvXvYuuvuXYYX01101000000100，那么可得出表示在整体坐标系中位移和应变关系的几何方程：eeeB

式中的几何矩阵[B]是自然坐标，的函数：

pNJJB1100011010000001

也可利用TiiipNNN求得的ipN以及iiiiiiYNYXNX,,4141和YXYXJ求出J，TppppYYYYXXXXNNNNJ432143214321。

3.单元刚度矩阵ek

设单元板厚为t，根据虚功方程有：tdABDBkATe，此式中几何矩阵[B]和弹性矩阵[D]都已求出。因为几何矩阵[B]中的变量是自然坐标，，所以也要用自然坐标表示微分面积dA。