第五章轴对称问题有限元法.pptVIP

有限元理论与应用 径向变形将引起周向应变，即 轴对称问题中，每个点具有四个应力分量 二、单元分析 1.位移函数 2.单元应变 3.单元应力 三、单元刚度矩阵 四、总刚集成 求出每个三角形环单元的刚度矩阵后，即可按照第二章介绍的总体刚度矩阵的集成方法，得出结构的总刚矩阵。 * 第五章 轴对称问题有限元法 第一节、轴对称问题的定义和特点 一、轴对称问题的定义 当分析结构同时满足以下三个条件时，可认为是轴对称问题。 （1）几何形状轴对称 要求结构是相对对称轴的旋转体。 （2）边界条件轴对称 载荷和位移约束条件具有轴对称性。 （3）材料轴对称 各向同性；各向异性 正交各向异性：材料主轴应与结构的径向、切向和轴向一致。 二、轴对称问题的应力应变特点 轴对称问题的特点是结构的位移、应变和应力都呈轴对称分布。 通常采用柱坐标系（ ），并以 轴为对称轴。结构 中的位移、应变和应力与角度 无关。 可以取出结构的任一子午面进行分析，从而将三维问题转化为二维问题 和四个应变分量 它们的关系为 为轴对称问题的弹性矩阵。 轴对称结构的几何模型是一个表示子午面形状的平面图形，与平面问题相比，轴对称问题的应力与应变分量各多一个。 一、结构离散 本身是一个三维结构，由于形状和载荷的特殊性，其网格划分仅在任一子午面上进行，为平面网格。 几何模型：一个表示子午面形状的平面图形，用相应的轴对称实体单元划分。 三节点三角形环单元 第二节 轴对称问题有限元法 选择线性位移函数 其中，Ni是形函数，其表达式为 将节点i,j,m的坐标值和位移值带入 选择线性位移函数，将节点i,j,m的坐标值和位移值带入 由此可见，周向应变分量 随 而改变，不是常量。 选择线性位移函数，将节点i,j,m的坐标值和位移值带入 用虚位移原理来推导三角形环单元的单元刚度矩阵。单元等效节点力所作 虚功等于三角环单元中的应力所作的虚功。 单元虚应变为： 考虑到虚位移的任意性，将两边的 同时消去， 此时节点力是指整个三角形环单元上的力 把单元中随位置变化而不断变化的r和z用单元截面的形心坐标来近似， 计算轴对称问题的等效节点载荷与平面问题有所不同，因为轴对称结构 的子午面上的一个节点是一个关于对称轴中心对称的圆环，故当计算集 中力、表面力和体积力时，应在整个环上积分。 1.集中力的移置 设三角形环单元内任意一点（r，z）处作用有集中外载荷Pc，它在r和z 方向上的分量分别为Pcr 和Pcz，用矩阵表示为 五、等效节点载荷的计算 Pc移置后对应的节点载荷列阵为 根据虚位移原理，等效节点载荷与原载荷在虚位移上作的虚功应相等，即 2.表面力的移置 O z r Ps l j i m n(r,z) s ds 单元的jm边作用有均布载荷Ps，其方向以压向单元边界为正。 表面力的矩阵表示为 由等效节点载荷与原载荷在虚位移上作的虚功应相等，同时，虚位移的任意性， 3.体积力的移置 工程中经常遇到的轴对称问题，其体积力一般有两种：一种是重力； 另一种是由于绕z轴旋转产生的惯性离心力。 （1）重力 设单元容重为v，则单元体积力为 单元体积力列阵为