第六讲三维实体单元和轴对称实体.ppt

6 三维实体单元和轴对称实体单元 一、三维弹性力学问题几何方程和物理方程 * * 1、几何方程 应变 ~ 位移关系： 矩阵形式几何方程： 2、物理方程 线弹性应力~应变关系： 其中应力分量： 各向同性弹性系数矩阵 为： （对称） 二、三维简单四面体单元 1、位移模式和形函数 四面体是最简单的三维单元，简单四面体单元把四个顶点作为节点，节点编号i，j，m，l。 每个节点3个自由度，单元共12个自由度。单元节点位移和单元节点力分量共12个。 简单四面体单元 因此，单元内每个坐标轴方向的位移设为三个坐标的完全一次多项式，以满足完备性条件： 用3节点三角形单元中同样的插值过程，进行广义坐标替换，得到用节点位移作广义坐标的位移模式： 用矩阵表示为： 其中形函数矩阵 表示为： 为三阶单位矩阵。 各节点形函数计算公式如下： 由行列式性质不难证明，该单元的形函数满足形函数性质！ 上式中 V为四面体的体积，由几何知识： （节点坐标行列式） 由于该单元是线性位移模式，单元的任一三角形界面在变形后仍然保持为一平面，界面位移由其上的3个节点位移决定。因此相邻两单元在变形过程中，三角形交界面上位移能保持协调，因此是协调单元，从而单元的收敛性得到保证。 2、应变矩阵和单元刚度矩阵 由几何方程得到单元内应变~节点位移关系： 应变矩阵 的任一个子块为： （k=i, j, m, l) 显然，应变矩阵是常数矩阵，该单元的应变、应力是常数，称为常应变单元，与平面三角形单元相似。 由单元刚度矩阵的通式得到简单四面体单元的刚度矩阵： 刚度矩阵按节点分块后，任一个子块为： 3、单元等效节点力计算 三维四节点四面体单元等效节点力计算公式导出原理和平面单元相同。 1）体力 2）面力 如果体力或面力是均匀力，可以把单元上的体力合力或一个三角形面上面力合力均匀分配到单元或一个面的所有节点上，即体力合力每个节点1/4，面力合力每个节点1/3。 4、单元评价 四节点四边形单元的优点是简单，适应性好。 由于采用线性位移模式，单元上应力、应变是常数，精度低，收敛速度慢，一般的单元网格使结构刚度明显偏大 。要得到一定准确度的结果，需要的单元和节点特别多，效益低。结构分析中应尽量避免使用该单元。 进行三维分析，一般采用高精度实体单元。 三、空间轴对称问题简介 工程中有大量的回转体结构，几何形状对称于中心轴。如果回转体所受到的载荷和约束条件也对称于中心轴，则其变形和应力也对称于此轴，这类问题在力学上称为空间轴对称问题。 典型的轴对称问题的例子：高速旋转的气轮机转子、轮盘，受均匀内压的旋转体高压厚壁容器（包括火炮身管），活塞，汽阀等各种密封件结构。 轴对称问题用圆柱坐标系 描述。Z轴为对称轴（回转中心轴），过中心轴的平截面（子午面）为r-z面。 所有力学量都与 无关，只是r和z的函数。因此轴对称问题只需在（r-z面）内描述，因为回转体的任何一个子午面都是圆周方向的对称面。因此，轴对称问题在数学上是二维问题。 轴对称问题的回转体上一点的集中力实际是沿过该点的圆周上的线分布力，表面（一段边界）上的分布力也是绕对称轴回转面上的轴对称分布力。 在轴对称问题中，弹性体内任意一点上，不存在切向位移，轴对称问题的位移有2个分量：沿径向r的位移u和沿轴向z的位移w。 轴对称问题有4个应变分量：三个子午面内的应变分量： ，一个圆周方向的应变分量 。轴对称问题的应变~位移关系——几何方程如下： 轴对称问题的应力分量有4个： 四、轴对称问题的简单三角形单元 1、轴对称问题离散化的几何意义 结构离散 对于轴对称问题的离散,通常在子午面roz上进行,其形状常为三角形和四边形。实际上，子午面上的每个三角形(或四边形)单元表示的是一个绕z轴一周的三棱(或四棱)环单元。因此,有限元轴对称问题的离散就是将连续体离散成由有限个圆环组成的离散体,单元与单元间通过环线(称为节线)相连接,作用于单元上的载荷,也作用于节线上。实际分析时，考虑到轴对称问题位移与周向无关，故可只需取一个截面,按平面情况进行分析。 2、位移模式和形函数 单元位移多项式与简单平面三角形单元相同，是完全一次多项式。通过插值得到节点位移表示的位移模式： 与平面问题相似，单元有6个自由度，每个节点有r，z方向2个位移。单元的节点位移列阵为： 形函数矩阵 与平面问题的简单三角形单元完全相同。 各形函数也完全相同，只是坐标不同： 轴对称问题三结点三角形单元的位移函数形式与平面问题三结点三角形单元基本相同， 主要区别在于结点位移分量不同。 3、应变与应力 4个应变分量中，面内（子午面） 3个应变分量为常量，环向应变不是常应变，