有限元-第4讲-轴对称问题与空间问题有限元法

第4章轴对称问题与空间问题有限元方法张洪伟轴对称问题的有限元法空间问题的有限元法12内容提要3一.轴对称问题的定义

工程中有一类结构，它们的几何形状、约束条件及作用的荷载都对称于某一固定轴(可视为子午面内平面物体绕轴旋转一周的结果)，其力学分析称为轴对称问题。1.离散化

由于可视为子午面内平面物体绕轴旋转一周的结果，因此轴对称问题分析可在子午面内划分单元，实际是取子午面内图形绕对称轴旋转所得“圆环形单元”对物体进行离散。因此可用的单元与平面问题一样。2.应力和应变

对轴对称问题进行分析一般取柱坐标系，对称轴为Z轴，径向为r轴，环向为θ轴。xyz

对称面上任一点p只在该对称面上发生位移，即所有的应力、应变和位移只是z和r的函数，而与坐标θ无关。那么轴对称问题就可转化为二维平面问题来进行研究。但因与平面问题有区别，常称二维半问题。轴对称问题的有限元法

如图所示的受均布内压作用的长圆筒，通过z轴的一个纵截面就是对称面。由于对称性，轴对问题共有4个应力分量：其中表示沿半径方向的正应力，称为径向应力；表示沿方向的正应力，称为环向应力或切向应力；表示沿ｚ方向的正应力，称为轴向应力；表示在圆柱面上沿ｚ方向作用的剪应力。

同样，轴对称问题共有4个应变分量：

其中表示沿半径方向的正应变，称为径向正应变；表示沿方向的正应变，称为环向正应变或切向正应变；表示沿ｚ方向的正应变，称为轴向正应变；表示沿ｒ和ｚ方向的剪应变。

在轴对称问题中，弹性体内任意一点上，不存在切向位移，只存在径向位移u

和轴向位移w

，两个位移分量表示为:xyz基本方程1.平衡方程2.几何方程3.物理方程

轴对称单元的特点(与平面三角形单元的区别)轴对称单元为圆环体，单元与单元间为节圆相连接；节点力与节点载荷是施加于节圆上的均布力；单元边界是一回转面；应变分量中出现了，即应变不是常量；且应变矩阵在r=0时，存在奇异点，需特殊处理，通常用该单元的形心坐标替代节点坐标。轴对称结构9轴对称问题的有限元法1.离散化

由于可视为子午面内平面物体绕轴旋转一周的结果，因此轴对称问题分析可在子午面内划分单元，实际是取子午面内图形绕对称轴旋转所得“圆环形单元”对物体进行离散。因此可用的单元与平面问题一样。2.单元分析

参照平面问题的三角形单元位移函数，轴对称问题的三结点三角形单元位移函数取为，其中：形函数：用矩阵表示的单元位移为：单元应变：

将单元位移函数带入几何方程得：其中，用几何矩阵表示单元的应变：

由于是坐标r、z

的函数，分量在单元中不为常量，其它三个应变分量在单元中仍为常量。应变矩阵不再是常数，轴对称三角形单元内的应变也不全为常量。单元应力：由弹性矩阵D

和应变矩阵B

可以得到应力矩阵S，并计算出单元内的应力分量：其中：式中：

由于应变矩阵中的元素不是常量，单元刚度矩阵需要通过积分得到，为简化计算可以用三角形单元形心位置的坐标代替B

矩阵中的变量，将单元中的r和z近似地当作常量，并且分别等于。

经过简化，就可以把各个单元近似地当作常应变单元。3.单元刚度矩阵有了单元应力场和应变场，可以利用虚位移原理或最小势能原理建立单元刚度矩阵，注意体积积分是沿单元的整个圆环求体积积分单元刚度矩阵的分块矩阵为，

由于几何矩阵中的元素不是常量，单元刚度矩阵需要通过积分得到，为简化计算可以用三角形单元形心位置的坐标代替B

矩阵中的变量。

实践证明采用近似积分也能达到一定的精度，具体对于三角形环单元用形心处坐标代替应变矩阵中的坐标变量。应变矩阵变成：其中：单元刚度矩阵的近似表达式为：单元刚度矩阵的分块矩阵近似表达式为：4.总刚度矩阵集成

求出了每一个三角形单元的刚度矩阵后，按照平面问题介绍的总刚矩阵的集成方法，就可以得到结构的总刚矩阵。5.等效节点载荷计算

计算轴对称问题的等效节点载荷与平面问题的有所不同，因为周对称结构的子午面上的一个节点是关于对称轴中心对称的圆环，故当计算集中力、表面力和体积力时，应在整个环上积分。（1）集中力移置集中力为

--集中力作用点的径向坐标。（2）体积力移置若体积力为自重，则单位体积的力为ρ为密度度如果单元离开对称轴较远（）可认为将1/3的自重移置到每个节点上。单元自重移置到节点i，j，m上的等效节点载荷为(i，j，m)若体积力为惯性离心力，则单位体积的力为单元离心力移置到节点i，j，m上的等效节点力(i，j，m)(3)分布面力移置

如果物体的形状、尺寸和边界条件不具备某种特殊性，不能简化为平面问题和轴对称问题处理，则必须作为空间问题来分析。与平面问题和轴对称问题相比，空间问题的计算要复杂的多。

（1）单元为块体形状。常用单元：四面体单元、长方体单元、直边六面体单元、曲边六面体单元、轴对称单元。（2）结点位移3个分量。（3）基本方程比平面问题多。3个平衡方程，6个几何方程，6个物理方程。

空间问题的有限元法4结点四面体单元：是空间问题最简单的单元，也是常应变、常应力单元，可以类似平面问题三结点三角形单元进行分析。8结点长方体单元：可以类似平面四结点矩形单元进行分析。8结点直边六面体单元：可以类似平面四结点任意四边形等参元分析。20结点曲边六面体单元：等参单元,可以类似平面八结点曲边四边形等参元进行分析。轴对称单元：一平面单元绕一对称轴旋转形成的空间问题。只需在rz平面划分网格，就像平面问题xy平面中的网格一样，这样这类空间问题可以得到简化。（环向位移等于零）1.四结点四面体单元分析1）位移函数单元结点位移向量单元位移函数xzyijmn将节点坐标代入位移函数

由此可解出代定常数α1~α4再代回到式单元位移函数的第1式，可得位移函数u:编号约定：当沿i,j,m的方向转动时，n在大拇指所指的方向采用同样的方法，可得单元位移：2）单元应变将单元中位移代入上式3）单元应力弹性矩阵[D]:代入单元应变计算公式，整理后: