有限元第八讲 三维实体单元和轴对称实体.pptVIP

7 三维实体单元和轴对称实体单元 一、三维弹性力学问题几何方程和物理方程 * * 1、几何方程 应变 ~ 位移关系： 矩阵形式几何方程： 2、物理方程 线弹性应力~应变关系： 其中应力分量： 各向同性弹性系数矩阵 为： （对称） 二、三维简单四面体单元 1、位移模式和形函数 四面体是最简单的三维单元，简单四面体单元把四个顶点作为节点，节点编号i，j，m，l。 每个节点3个自由度，单元共12个自由度。单元节点位移和单元节点力分量共12个。 简单四面体单元 因此，单元内每个坐标轴方向的位移设为三个坐标的完全一次多项式，以满足完备性条件： 用3节点三角形单元中同样的插值过程，进行广义坐标替换，得到用节点位移作广义坐标的位移模式： 用矩阵表示为： 其中形函数矩阵 表示为： 为三阶单位矩阵。 各节点形函数计算公式如下： 由行列式性质不难证明，该单元的形函数满足形函数性质！ 上式中 V为四面体的体积，由几何知识： （节点坐标行列式） 由于该单元是线性位移模式，单元的任一三角形界面在变形后仍然保持为一平面，界面位移由其上的3个节点位移决定。因此相邻两单元在变形过程中，三角形交界面上位移能保持协调，因此是协调单元，从而单元的收敛性得到保证。 2、应变矩阵和单元刚度矩阵 由几何方程得到单元内应变~节点位移关系： 应变矩阵 的任一个子块为： （k=i