第6章 轴对称有限元分析.pptVIP

轴对称问题是空间问题的一种特殊情况，在实际工程中存在大量的轴对称问题。 飞轮、回转类的压力容器、发动机汽缸套、烟囱及受内压的球壳等，无限大、半无限大的弹性体受集中载荷作用时也可以处理为轴对称问题。 单元位移函数 参照弹性平面问题有限元法，单元位移函数为 几何方程 尽管点的位移发生在平面内，但是，对于垂直于平面的线元素却存在着伸缩的可能，因此，轴对称问题的环向应变不为零。 几何方程 对于周向应变，尽管不存在周向位移，但由于A点发生径向位移后，它与轴的距离变为，从而导致产生周向的变形，如图所示，则产生周向应变为 几何方程 矩阵表示为 几何方程 式是[B] 是几何矩阵（应变转换矩阵）。 物理方程 由于圆柱坐标也是正交坐标，相应的物理方程为： 物理方程 D为轴对称问题的弹性矩阵： 刚度方程（平衡方程的等价形式） 其中 等效节点载荷 对轴对称问题，按静力等效原则把轴对称外加载荷和体积力转化到单元的节点上，其载荷仍然是轴对称载荷，是分布在圆周节点上的线载荷，把它们集中在一起称为轴对称问题的节点载荷。 体积力 惯性力 面力 线载荷 体积力 惯性力 面力 线载荷 其中 (4-42) 应力公式形式上与前面轴对称三角形单元相应公式相同，即 (4-43) 但是 (4-44) 其中常数 (4-45) 当r=0时，即对称轴上有 ，因此，可以用 来代替 单元刚度矩阵具有以下的形式 其中每一个子矩阵是 (4-47) 以消除(4-44)式中的奇异项。 (4-46) 而 (4-48) 式中常数A1、A2、A3按(4-45)式计算。 等效结点力可以按照以下公式计算 (1) 集中力 设在单元径向坐标为rc的某点作用有集中力 移置到结点上的等效结点力，可按下式计算 式中(Ni)c是形函数Ni在载荷作用点的值。 (4-49) (2) 体积力 设单元体积力是 ，则作用在单元各结点上的等效结点力分别按下式计算 式中r用(4-36)式的第一式代入即可。 (4-50) (3) 表面力 设单元的某边界上作用有表面力 ，其中 和 分别是单元表面力在作用边外法线方向和切线方向的投影。于是，此边上各结点的等效结点力分别按下式计算： 式中r可用(4-36)式的第一式代入。比如， 是对应 的那条边，那么 (4-51) (4) 温度应力 如果考虑温度变化产生的初应变 则任意结点上的等效结点力是 (4-52) 计算应力的公式(5-52)(浙大有限元P65)式要改写为 (4-53) 1. 什么叫轴对称问题？其位移有何特点？ 2. 试定义具有线性位移模式的三角形环单元。 3. 三角形环单元的应力和应变的特点。 作业： 整体刚度矩阵 于是(f)式便可以写成与平面问题相同的标准形式 这就是求解结点位移的平衡方程组。 (6-6) (g) 整体刚度矩阵也可以写成分块形式 其中子矩阵为 和平面问题一样，整体刚度矩阵 是对称的带状稀疏阵，在消除刚体位移后，它是正定的。 式右边的载荷列阵展开的形式为 其中 (b) (a) 与平面问题一样，等效结点力也是由作用在环形单元上的集中力、表面力和体积力分别移置到结点上而得到的。移置的原则也是根据这些力和等效结点力在任意虚位移上所作的虚功相等，即 等效节点载荷 式中 rc 为集中载荷 作用点的径向坐标。将上节(b)式代入上式并考虑到 ，上式可以化成 式中右边第一项是环形单元上的集中力 移置到结点的等效结点力，第二项是环形单元上表面力 的等效结点力，第三项是环形单元体积力 的等效结点力。 (c) 采用平面问题中相同的符号： 集中力的等效结点力 表面力的等效结点力 体积力的等效结点力 (6-9) (6-7) (6-8) 集中力等效节点力 表面力等效节点力 体积力等效节点力 (i) (ii) (iii) 直角坐标 于是(c)式可以改写成 再将上式代入(a)式，等效载荷列阵可写成 将(6-8)、(6-9)式和(ii)、(iii)式比较可见，在轴对称情况下积分号后的被积函数比平面问题的多一个变量r，所以虽然也是采用线性位移模式，但是不能象平面问题那