《常用单元的刚度矩阵》由会员分享,可在线阅读,更多相关《常用单元的刚度矩阵(9页珍藏版)》请在金锄头文库上搜索。
1、.由于各点在圆周方向上无位移 ,因而剪应变vr和vr均为urrur零。将应变写成向量的形式,则wzzrzuwzr根据上式,可推导出几何方程B(e)zjkN (r,z)ir10其中几何矩阵B2 rkj00rkjzjkzkjNj(r,z)r0rik00rikzkizijNk(r,z)r0rji0 0 rjizij3.弹性方程和弹性矩阵D依照广义虎克定律 ,同样可以写出在轴对称中应力和应变之间的弹性方程,其形式为所以弹性方程为D式中应力矩阵rzrzT0110E弹性矩阵D10(1)(12)1200024.单元刚度矩阵k(e)与平面问题相同,仍用虚功原理来建立单元刚度矩阵 ,其积分式为在柱面坐标系中,d
2、V 2drdz1 / 9.将dV 2drdz代入k(e)BTDBdV,则k(e) 2BTDBrdrdzV即为轴对称问题求单元刚度矩阵的积分式。与弹性力学平面问题的三角形单元不同 ,在轴对称问题中,几何矩阵B有的元素如Ni(r,z)等是坐标 r、z 的函数,不r是 常 量 。 因 此 , 乘 积BTDB不 能 简 单 地 从 式k(e) 2BTDBrdrdz的积分号中提出。如果对该乘积逐项求积分,将是一个繁重的工作。一般采用近似的方法:用三角形形心的坐标值代替几何矩阵B的 r 和 z 的值。 用B表示在形心(r,z)处计算出的矩阵B。其中只要单元尺寸不太大,经过这样处理引起的误差也不大。被 积
3、函 数 又 成 为 常 数 , 可 以 提 出 到 积 分 号 外 面 :k(e) 2BDBrdrdz 2BDBr式中三角形的面积。TT由式k(e) 2BDBrdrdz 2BDBr可以看出,两轴对TT称的三角形单元,当形状、 大小及方位完全相同而位置不同时,其刚度矩阵也不相同。距离主轴线越远的单元,其刚度越大。这与平面问题不一样。二、等参数的刚度矩阵对一些由曲线轮廓的复杂结构 ,如果采用直角边单元进行离散,由于用直线代替了曲线,除非网格划分得很细,否则不能获得较高的精度;对另一些应力随坐标急剧变化的结构 ,采用简单的常应力单元离散时 ,也必须划分成大量的微小单2 / 9.元,以保证足够的精度。
4、 为此引入一种高精度的单元等参数单元。 它既能简化复杂单元划分的工作,又能在满足同样精度的要求时,大大减少使用的单元数。 目前流行的大程序中较常用,它成功地解决了许多二维和三维的弹性力学问题。为导出等参数单元的刚度矩阵 ,首先要建立根据每个单元的形状确定的自然坐标系 ,然后将位移模式和形状函数都写成自然坐标的函数。一个单元在自然坐标系的点余元整体坐标系的点成一一对应的关系。 通过映射,可以将整体坐标系中的图形转化为自然坐标系中的相应徒刑。例如可以将整体坐标系中的一个任意四边形实际单元映射到自然坐标系中成为一个正方形基本单元。同样也可以将任意四面体、六面体包括直边和曲边的分别映射成正四面体和正六
5、面体。这里只介绍较简单的一种平面问题的情况 ,将整体坐标系中的一个任意四边形映射成自然坐标系中的一个正方体 ,并导出单元刚度矩阵。其它种单元的映射,可依次原理进行。不再叙述。1. 位移模式和形状函数图 4-2 中的任意四边形单元上,作连接对边中点的直线,取其交点为原点,这两条直线分别为和轴,并令四条边上的和值分别为1,建立一个新的坐标系,称之为该单元的自然坐标系。原坐标系 XOY 称为整体坐标系。在整体坐标系中,3 / 9.自然坐标系非正交,它由任意四边形的形状所确定。图图 4-194-19如果将自然坐标系改画成直角坐标系,那么图 4-19a 中的任意四边形单元就成为图 4-19b 所示的正方
6、形。上述两个四边形的点 包括顶点一一对应,即它们之间相互映射。因此,需要写出整体坐标 X、Y 和自然坐标、之间的坐标转换式,即X 1234Y 5678*四边形四个顶点的坐标值在XOY 坐标系中分别为X1,Y1,X2,Y2,X3,Y3,X4,Y4: 在o坐 标 系 中 相 应 为1,1,1,1,1,1,1,1。将有关数据代入*中的第一式,则有求解上述方程组得:坐标变换方程*成为X 11X11X21X31X44同理Y 11Y11Y21Y31Y44当引入函数Ni,后,坐标变换方程成为式中Ni,11i1i4变量、的正负号由相应节点的坐标值i、i决定。例如当 i=4 时, 1,1,因此,N,11。444
7、4下面再来研究函数Ni,的特性。对节点1X1,Y1,相应的自然坐标值为-1,-1。从式4 / 9.Ni,11i1i中很容易看出,除4N1=1 外,N2=N3=N4=0。对其余各节点也一样。总而言之,对节点 i,除 Ni=1外,其余三个 N 值均为零。同时,不难看出N1, N2, N3, N4,1,即四个节点的 Ni函数之和等于 1。函数Ni,具备上章所介绍的形状函数应满足的条件,可作为本单元的形状函数。采用Ni,做形状函数,其位移模式为对比X Ni,Xii144Y Ni,Yii1和u Ni,ui,v Ni,vi可以看出:i1i144在这种实际单元任意四边形中,坐标变换式和位移模式不仅采用了相同
8、的形状函数Ni,而且具有相同的数学模型。这种性质的实际单元称为等参数单元。对用节点位移值 ui或 vi等求单元某一点位移量 u的插值公式u Ni,ui,只要将 u换成 X,i14便成为利用节点值 Xi求相应点坐标 X的插值公式。相反也是这样。2.几何矩阵B由于几何矩阵B通过对位移求偏导数而得出,所以首先必须利用复合函数求导的规则得出下述公式5 / 9.X式中J XY ,此式称为雅可比矩阵。Y为了将几何矩阵B写成变量、的函数,必须将式u u XuJu改写成y u v u v XX11uJu,同理,vJvy y 从表示单元各点位移与其应变关系的几何方程可知:u v u v XX11将式uJu和vJ
9、v合并,则y y 对单元e,任意一点的位移 u,v 对自然坐标,的偏导数可利用上式求出,写成矩阵形式为:eu1式中v1u1v1u1v1u1v1Tuvv T对于 i=1,2,3,4将uuuXuYTNipNiNiTvXTJ1v Y00uJ1和入vv e Np 代6 / 9. X 0Y01000uu 0001Yv0110XuYvXv YT,则可得出表示在整体坐标系中位移和应变关系的几何方程:eBee式中的几何矩阵B是自然坐标,的函数:Ni也可利用NipNiT求得的Nip以及X Ni,Xii144和Y Ni,Yii1XJ XY X1求出J,JN1pN2pN3pN4pYY1 X2Y2X3Y3X4Y4T。
10、3.单元刚度矩阵ke设单元板厚为 t,根据虚功方程有:keBTDBtdA,此式A中几何矩阵B和弹性矩阵D都已求出。因为几何矩阵B中的变量是自然坐标,,所以也要用自然坐标表示微分面积dA。在实际单元中任取一点 p,其整体坐标位 X、Y,其相应的自然坐标为,。过 p 点做,的等值线,同时做 d, d的等值线,围成一小块微分面积 dA,如图 4-20a 所示。为便于分析,将四边形 pqrs 放大,如图 4-20 b 所示。 实际上,d,d取得很小,因此该四边形可视为平行四边形。 若相邻的两边用, b表示,则两者的乘积等于该平行四边形的面积 dA。向量a图图 4-204-207 / 9.若a axi
11、ayj,b bxi byjbxby,则axdA axi ayjbxi byjay为了求出ax,ay,bx,by的值,要先写出 a 和 b 两端节点 p、q、s的坐标值。点 p:Xp点 q:Xq点 s:Xs X,Yp Y, X d,Yq Y d, X, d,Ys Y, d利用泰勒技术展开并略去高阶项,可得对Y d,Y, d,也可写出相应的展开式。利用式Xd可得:XX, d X,dX d, X,XYd,ay YqYpd,XYbx Xs Xpd,by YsYpdax Xq Xpax将此式代入式dA axi ayjbxi byjayaxayXbxXbybxby得到:dA Ydd,简写为dA J ddY
12、单元刚度矩阵为:ke11T BDBtJdd,这个积分可以采用数值方法,用高斯1 1求积分公式很方便的求出,在此不作介绍。8 / 9.例:求如图所示四边形的雅可比矩阵。解: 求雅可比矩阵可在整体坐标系中进行,也可以在实际单元的局部坐标系中进行。为便于计算 ,本例在局部坐标系中进行。对单元1:将四个节点的自然坐标值-1,-1、 1,-1、 1,1、 -1,1 代入下式:Ni,11i1i,再将所得到的Ni,值及四个节点4实际单元在局部坐标系的坐标值-3,-2、 3,-2、 3,2、 -3,2代入下式计算:X Ni,Xii144,则X 3,Y 2Y Ni,Yii1X雅克比矩阵为J XY 30。Y02对单元2:四个节点 在局 部坐标系 中的 坐标值分 别为 -1,-3/4、1,-3/4、 1,5/4、 -1,1/4。因此,因而雅可比矩阵为也可利用JN1pN2pN3pN4pX1Y1X2Y2X3Y3X4Y4T求雅可比矩阵 ,其结果与上相同,同学们可自行验证。图图9 / 9