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

式中子矩阵[krs]为3×3的矩阵 : 5) 等效结点荷载 体积力与表面力的计算公式与平面三角形单元公式 相似，可以采用静力等效原则简化计算。 崇尚实践 知行并重 * ? BIPT ? BIPT ? BIPT 第4章 轴对称问题与空间 问题有限元方法 张 洪 伟 轴对称问题的有限元法 空间问题的有限元法 1 2 内容提要 * 一. 轴对称问题的定义 工程中有一类结构，它们的几何形状、约束条件及作用的荷载都对称于某一固定轴(可视为子午面内平面物体绕轴旋转一周的结果)，其力学分析称为轴对称问题。 1.离散化 由于可视为子午面内平面物体绕轴旋转一周的结果，因此轴对称问题分析可在子午面内划分单元，实际是取子午面内图形绕对称轴旋转所得“圆环形单元”对物体进行离散。因此可用的单元与平面问题一样。 2.应力和应变 对轴对称问题进行分析一般取柱坐标系，对称轴为Z轴，径向为r 轴，环向为θ轴。 x y z 对称面上任一点p只在该对称面上发生位移，即所有的应力、应变和位移只是z和r的函数，而与坐标θ无关。那么轴对称问题就可转化为二维平面问题来进行研究。但因与平面问题有区别，常称二维半问题。 轴对称问题的有限元法 　　如图所示的受均布内压作用的长圆筒，通过z轴的一个纵截面就是对称面。由于对称性，轴对问题共有4个应力分量： 其中　 表示沿半径方向的正应力，称为径向应力； 表示沿 方向的正应力，称为环向应力或切向应力； 表示沿ｚ方向的正应力，称为轴向应力；　 表示在圆柱面上沿ｚ方向作用的剪应力。 　　同样，轴对称问题共有4个应变分量： 其中　表示沿半径方向的正应变，称为径向正应变；　表示沿　方向的正应变，称为环向正应变或切向正应变；　表示沿ｚ方向的正应变，称为轴向正应变；　表示沿ｒ和ｚ方向的剪应变。 　　在轴对称问题中，弹性体内任意一点上，不存在切向位移，只存在径向位移 u 和轴向位移 w ，两个位移分量表示为: x y z 基本方程 1.平衡方程 2.几何方程 3.物理方程 轴对称单元的特点(与平面三角形单元的区别) 轴对称单元为圆环体，单元与单元间为节圆相连接； 节点力与节点载荷是施加于节圆上的均布力； 单元边界是一回转面； 应变分量 中出现了 ，即应变不是常量；且应变矩阵在r=0时，存在奇异点，需特殊处理，通常用该单元的形心坐标替代节点坐标。 轴对称结构 * 轴对称问题的有限元法 1.离散化 由于可视为子午面内平面物体绕轴旋转一周的结果，因此轴对称问题分析可在子午面内划分单元，实际是取子午面内图形绕对称轴旋转所得“圆环形单元”对物体进行离散。因此可用的单元与平面问题一样。 2.单元分析 参照平面问题的三角形单元位移函数，轴对称问题的三结点三角形单元位移函数取为， 其中： 形函数： 用矩阵表示的单元位移为： 单元应变： 将单元位移函数带入几何方程得： 其中， 用几何矩阵表示单元的应变： 　　由于　是坐标 r、z 的函数，　分量在单元中不为常量，其它三个应变分量在单元中仍为常量。应变矩阵 不再是常数，轴对称三角形单元内的应变也不全为常量。 单元应力： 由弹性矩阵 D 和应变矩阵 B 可以得到应力矩阵 S，并计算出单元内的应力分量： 其中： 式中： 　　由于应变矩阵中的元素不是常量，单元刚度矩阵需要通过积分得到，为简化计算可以用三角形单元形心位置的坐标 　　 代替 B 矩阵中的变量，将单元中的r和z近似地当作常量，并且分别等于 。 　　经过简化，就可以把各个单元近似地当作常应变单元。 3.单元刚度矩阵 有了单元应力场和应变场，可以利用虚位移原理或最小势能原理建立单元刚度矩阵，注意体积积分是沿单元的整个圆环求体积积分 单元刚度矩阵的分块矩阵为， 　　由于几何矩阵中的元素不是常量，单元刚度矩阵需要通过积分得到，为简化计算可以用三角形单元形心位置的坐标 　　代替 B 矩阵中的变量　　。 实践证明采用近似积分也能达到一定的精度，具体对于三角形环单元用形心处坐标代替应变矩阵中的坐标变量。 应变矩阵变成： 其中： 单元刚度矩阵的近似表达式为： 单元刚度矩阵的分块矩阵近似表达式为： 4.总刚度矩阵集成 　　求出了每一个三角形单元的刚度矩阵后，按照平面问题介绍的总刚矩阵的集成方法，就可以得到结构的总刚矩阵。 5.等效节点载荷计算 　　计算轴对称问题的等效节点载荷与平面问题的有所不同，因为周对称结构的子午面上的一个节点是关于对称轴中心对称的圆环，故当计算集中力、表面力和体积力时，应在整个环上积分。 （1）集中力移置 集中力为 　　--集中力作用点的径向坐标。 （2）体积力移置 若体积力为自重，则单位体积的