0302空间问题四面体单元解析.docVIP

03-02空间问题的周围体单元解析 03-02空间问题的周围体单元解析 PAGE / NUMPAGES 03-02空间问题的周围体单元解析 第三章 轴对称、三维和高次单元 §3-2 空间问题的周围体单元 空间问题的有限单元法，和平面问题及轴对称问题的有限单元法的原理和解析过程完好同样。由于空间问题应采用三维坐标系，因此单元的自由度、刚度矩阵的元素个数，方程组内方程个数等要较平面问题和轴对称问题多，因此空间问题的规模一般比轴对称问题和平面问题大得多。它要求计算机的内存大，且计算时间长，花销高。这些问题都给三维有限单元法的详尽运用带来好多困难。 和平面问题同样，空间有限单元法采用单元 也是多种多样的，其中最简单的是四节点周围体 单元。采用周围体单元和线性位移模式来办理空 间问题，可以看作平面问题中三角形单元的实行。 在采用周围体单元失散化后的空间结构物 中，一系列不相互重叠的周围体之间仅在节点处 以空间铰相互连接。四节点周围体单元仅在四个 极点处取为节点， 其编号为 i,j,m,p 。每个单元的计算简图如图 3-7 所示。 在位移法中，取节点位移为基本未知量，四 节点周围体单元共有十二个自由度 ( 位移重量 ) ，其节点位移列阵为  图 3-7 空间周围体单元 i e j m p ui vi wi u j v j w j um vm wm u p vp wp T 其子矩阵 相应的节点力列阵为  i ui vi wi (i,j,m) F e Fi F j Fm Fp T 1 其子矩阵 Fi U i Vi Wi T 一、单元法位移函数 结构中各点的位移是坐标 x、 y、 z 的函数。当单元足够小时，单元内各点的位移可用简单的线性多项式来近似描述，即 u v w  1 2 x 3 y 4 z 5 6 x 7 y 8 z (3-49) 0 10 x 11 y 12 z 式中 1， 2，, ， 12 是十二个待定系数，它们可由单元的节点位移和坐标确定。假设节 点 i,j,m,p 的坐标分别为 ( x i y i z ) 、( x j y j z j ) 、( xm y m zm ) 、 ( x p y p z ) ， i p 将它们代入 (3-49) 式的第一式可得各个节点在 x 方向的位移 ui 1 2 xi 3 yi 4 zi u j 1 2 x j 3 y j 4 z j (3-50) um 2 xm 3 ym 4 zm 1 u p 1 2 x p 3 y p 4 z p 解上述线性方程组，可获得 1 ， 2 ， 3 ， 4 ，再代入 (3-50) 式，得 u 1 [( ai bi x ci y d i z)ui (a j b j x c j y d j z)u j (3-51) 6V ( am bm x cm y d m z)um ( ap bp x c p y d p z)u p ] 其中 V 为周围体 ijmp 的体积， a ,b i ，, ， c ,d p 为系数。 i p 1 xi yi zi V 1 x j y j zj (3-52) 1 xm ym zm 1 x p y p zp 2 x j y j zj 1 y j zj ai xm ym zm bi 1 ym zm x p y p z p 1 yp zp x j 1 z j x j y j 1 ci xm 1 zm di xm ym 1 (i,j,m,p) (3-53) x p 1 zp x p y p 1 为了使周围体的体积 V 不致为负值，单元四个节点的标号 i,j,m,p 必定依照必然的次序： 在右手坐标系中，要使得右手螺旋在依照 i → j → m的转向转动时，向 p 的方向前进，象图 3-1 中单元那样。 用同样方法，可以得出其他二个位移重量： v 1 [( ai bi x ci y di z)vi (a j b j x c j y d j z)v j (3-54) 6V (am bm x cm y dm z)vm (a p bp x cp y d p z)v p ] w 1 [( ai bi x ci y di z)wi ( a j b j x c j y d j z) w j 6V (3-55) (am bm x cm y dm z)wm (a p bp x cp y d p z)wp ] 综合表达式 (3-51) 、 (3-54) 及 (3-55) ，可以将位移重量表示成为 f u v w T [ N ] e (3-56) IN i I