常用单元的刚度矩阵

1、 2兀(r+u)一2兀rur9由于各点在圆周方向上无位移，因而剪应变v和v均为0rdu零。将应变写成向量的形式，则=r80=8zYrzdrurdwdzdudw+一、dzdr一根据上式，可推导出几何方程S=|bUe)其中几何矩阵b=丄2AZjkN(r,z)ir0ZkjN(r,z)jZijN(r,z)krkjrkjzjkrikrikZkirjirjizij3.弹性方程和弹性矩阵D依照广义虎克定律，同样可以写出在轴对称中应力和应变之间的弹性方程，其形式为+0)+Q)+0)rrz2(1+卩)Trz所以弹性方程为0暑怡式中应力矩阵QUttT口0zrz弹性矩阵dLe（1+卩）（1一2卩）1|LX01|LX

2、01|LX000012|LX24.单元刚度矩阵kL）与平面问题相同，仍用虚功原理来建立单元刚度矩阵，其积分式为kle）二JbtdBJdVV在柱面坐标系中，dv二2兀drdz将dV二2.drdz代入以）JbdBdV，则订）=2兀JJbDBdrdzV即为轴对称问题求单元刚度矩阵的积分式。与弹性力学平面问题的三角形单元不同，在轴对称问题中，几何矩阵B内有的元素（如N（r，z）等）是坐标r、z的函r数，不是常量。因此，乘积blblBl不能简单地从式二2川B卜dlBldrdz的积分号中提出。如果对该乘积逐项求积分，将是一个繁重的工作。一般采用近似的方法：用三角形形心的坐标值代替几何矩阵B内的r和z的值。

3、用胃1表示在形心（衣）处计算出的矩阵B。其中-（r+r+r）_（z.+z.+z）ijkijk3，z只要单元尺寸不太大，经过这样处理引起的误差也不大。被积函数又成为常数，可以提出到积分号外面：Ik1e）=2兀h】b由式kle）二2兀B1d11rdrdz二2兀B】dJb可以看出，两轴对称的三角形单元，当形状、大小及方位完全相同而位置不同时，其刚度矩阵也不相同。距离主轴线越远的单元，其刚度越大。这与平面问题不一样。二、等参数的刚度矩阵对一些由曲线轮廓的复杂结构，如果采用直角边单元进行离散，由于用直线代替了曲线，除非网格划分得很细，否则不能获得较高的精度；对另一些应力随坐标急剧变化的结构，采用简单的常

4、应力单元离散时，也必须划分成大量的微小单元，以保证足够的精度。为此引入一种高精度的单元等参数单元。它既能简化复杂单元划分的工作，又能在满足同样精度的要求时，大大减少使用的单元数。目前流行的大程序中较常用，它成功地解决了许多二维和三维的弹性力学问题。为导出等参数单元的刚度矩阵，首先要建立根据每个单元的形状确定的自然坐标系，然后将位移模式和形状函数都写成自然坐标的函数。一个单元在自然坐标系内的点余元整体坐标系内的点成一一对应的关系。通过映射，可以将整体坐标系中的图形转化为自然坐标系中的相应徒刑。例如可以将整体坐标系中的一个任意四边形（实际单元）映射到自然坐标系中成为一个正方形（基本单元）。同样也可

5、以将任意四面体、六面体（包括直边和曲边的）分别映射成正四面体和正六面体。这里只介绍较简单的一种平面问题的情况，将整体坐标系中的一个任意四边形映射成自然坐标系中的一个正方体，并导出单元刚度矩阵。其它种单元的映射，可依次原理进行。不再叙述。1.位移模式和形状函数图4-2中的任意四边形单元上，作连接对边中点的直线，取其交点为原点，这两条直线分别为g和耳轴，并令四条边上的g和耳值分别为土1,建立一个新的坐标系，称之为该单元的自然坐标系。原坐标系XOY称为整体坐标系。在整体坐标系中，自然坐标系非正交，它由任意四边形的形状所确定。图4-19如果将自然坐标系改画成直角坐标系，那么图4-19（a）中的任意四边

6、形单元就成为图4-19（b）所示的正方形。上述两个四边形的点（包括顶点）一一对应，即它们之间相互映射。因此，需要写出整体坐标X、Y和自然坐标g、耳之间的坐标转换式，即X二a+ag+aq+agq*1234Y=a+ag+aq+agq5678四边形四个顶点的坐标值在XOY坐标系中分别为（X,Y）,（X,Y）,（X,Y）,（X,Y）：在gon坐标系中相应为11223344（-1,-1）,Cl）6,1）（一1,1）。将有关数据代入*中的第一式，则有 1234 X=aex一ex+a,X=a+aex一exTOC o 1-5 h z1123421234X=x+x+x+x,X=xex+xex3123441234

7、求解上述方程组得：X+X+X+XX+X+XXa二一123,a二12341424XX+X+XXX+XXa二1234,a二一12343444坐标变换方程*成为二-K1-g-q+gq)X+(1+g-q-gq)X+(1+g+q+gq)X+(1-g-q-gq)X41234同理Y二K1gq+gq)Y+(1+gqgq)Y+G+g+q+gq)Y+(1gqgq)Y41234当引入函数n(g,n)后，坐标变换方程成为iX=工N/q)Xiii=1Y=N(g,n)Yiii=1式中N(g,q)=1(1+gg)(i4ii变量g、耳的正负号由相应节点的坐标值g、耳决定。例如ii当i=4时，g=1,n=1，因此，n(g,n)

8、=(1g)(+n)。4444下面再来研究函数N(g,耳)的特性。i对节点1(X,Y)，相应的自然坐标值为(-1，-1)。从式11N(g,n)=1(1+gg)(+nn)中很容易看出，除N1=1外，i4ii1N2=N3=N4=0。对其余各节点也一样。总而言之，对节点i(i=1,2,3,4)，除N.=1夕卜，其余三个N值均为零。同时，不难看出n(g,n)+n(g,n)+n(g,n)+n(g,耳)=1，即四个节点的斗函数之和等于1。函数N，丿具备上章所介绍的形状函数应满足的条件，i可作为本单元的形状函数。采用N匕做形状函数，其位移模式为iu=XN,v）u,v=工N，耳丄iiiiTOC o 1-5 h

9、zi=1i=1x=Xn毛,n）x对比i=1ii和u=XN电,心,v=XN（g“可以看出：4iiiiy=乙n（g,nYi=1i=1iii=1在这种实际单元（任意四边形）中，坐标变换式和位移模式不仅采用了相同的形状函数N匕，而且具有相同的数学i模型。这种性质的实际单元称为等参数单元。对用节点位移值u.（或V.等）求单元内某一点位移量u（或vii等）的插值公式u=XN（g,心，只要将u（或V等）换成X（或Yiii=1等），便成为利用节点值X.（或Y.等）求相应点坐标X（或Y等）的插值公式。相反也是这样。2.几何矩阵B由于几何矩阵B通过对位移求偏导数而得出，所以首先必须利用复合函数求导的规则得出下述公

10、式dududXdudYdu何=SdXdgdYdg或写成鬲=J1dududXdudYldu、耳十dXdndYdvaxQY式中J=正鬲，此式称为雅可比矩阵。axQY为了将几何矩阵B写成变量耳的函数，必须将式QuQuaxQu=JI1=JaxQu改写成lan)一从表示单元内各点位移与其应变关系的几何方程可知：aax爲oaaY0ax001auagauaYavaxQuQuQvQv将式axQu=JI1=JI1鬲Qv合并，则一lan)一lan)QvaxQvTJliaYJ_o0JliQuQuaauaYav对单元（e）,任意一点的位移u，v对自然坐标g,耳的偏导数可利用上式求出，写成矩阵形式为：vu叫讥时vuii

11、2pipvuvTkrlokInloL.10Lt)tv3p4p对于i=l,2,3,4InL竺竺!iP况眄J将dudYdvdXdvdYJ01d0ddX=oddYdYddX0000011J110J11dududvdqdg和dudg代入色空色r，则可得出表示dYdXdYJ在整体坐标系中位移和应变关系的几何方程：bh)=Be)|/h)式中的几何矩阵B是自然坐标g耳的函数:也可利用陰J110竺JT求得的n以及dqJiPX=工N匕q)x匸1ii和Y=N毛,qiii=11P2P3P4PXXXX1234YYYY1234dXdYJ=勇求出/，dXdYdqdq3.单元刚度矩阵kq)设单元板厚为t,根据虚功方程有：肛

12、)“国blBjdA，A此式中几何矩阵B和弹性矩阵D都已求出。因为几何矩阵B中的变量是自然坐标g,q，所以也要用自然坐标表示微分面积dA。在实际单元中任取一点p,其整体坐标位X、Y,其相应x 的自然坐标为g,耳。过p点做g,耳的等值线，同时做g+dg,n+dn的等值线，围成一小块微分面积dA,如图4-20(a)所示。为便于分析，将四边形pqrs放大，如图4-20(b)所示。实际上，dg，dn取得很小，因此该四边形可视为平行四边形。若相邻的两边用向量表示，则两者的乘积等于该平行四边形的面积dA。图4-20dA=|a|bsin0=|axb|若a=ai+aj,b=bi+bjxyxy()()abJi+a

13、bi+bjA=xxxyxyabyydA二为了求出a，a,b，b的值，要先写出a和b两端节点p、q、sxyxy的坐标值。点p：x二x(g,n),y二y(g,n)pp点q：X二X(g+dg,n),Y二Y(g+dg,n)qq点s：x二x(g,n+dn)y二y(g,n+dn)ss利用泰勒技术展开并略去高阶项，可得x(g+dg,n)=x(g,n)+|gdgx(g,n+dn)=x(g,n)+Xdn的对y(g+dg,n)y(g,n+dn)，也可写出相应的展开式。利用式x-dg,n)=X,n)+2dg可得：,n+dn)=x,n)+空dnon TOC o 1-5 h zOXdYa=XX=dg,a=YY=dqxq

14、pOgyqp術OXgOYqb=XX=dg,b=YY=dqab5i+aj，丿bi+bj彳=xxxyxyab得到：将此式代入式dA=yy简写为dA=JdgdqxspOgyspOqadA=xaybxbyOXOYOgOgdgdq，OXOY71OqOq单元刚度矩阵为：订e）=JfIbbBJJLgdq，这个积分可以采用“数值方法”用11高斯求积分公式很方便的求出，在此不作介绍。例：求如图所示四边形的雅可比矩阵。解：求雅可比矩阵可在整体坐标系中进行，也可以在实际单元的局部坐标系中进行。为便于计算，本例在局部坐标系中进行。对单元（1）：将四个节点的自然坐标值（-1，-1）、（1，-1）、（1，1）（-1，1）代入下式：N（g,q）=1G+gg）G+qq），四个节点i4iii实际单元在局部坐标系的坐标值（-3，-2）、（3，-2）、（3，2）（-3，2）代入下式计算：X=工Nm)Xi=i1，贝X二3g,Y二2qY=工N,n)Yii