第5章轴对称问题的有限单元法.ppt

6.1 轴对称问题的有限单元法 4.2 轴对称体的离散化 4.3 三角形环单元的位移函数 4.4 三角形环单元的应变和应力 4.5 三角形环单元的单元刚度矩阵 第4章 轴对称体的有限单元法[专题2] 中南大学 * 在工程问题中经常会遇到一些实际问题，如旋转机械中的盘、轴、承力环及支撑圈等，它们都有一个对称轴，而整个物体是通过轴的一个平面上某个图形绕此轴旋转而成的回转体，称之为轴对称体。 ?如果轴对称体的载荷是轴对称的，其约束也是轴对称的，则在载荷作用下产生的位移、应变和应力必然是轴对称的，这种结构的应力分析问题称轴对称问题。(结构轴对称＋载荷轴对称＋约束轴对称) ?如果轴对称体的载荷是复杂载荷，如弯矩、扭矩等，此时载荷不是对称的，则在载荷作用下产生的位移、应变和应力也不再是轴对称的。这种结构的应力分析问题称轴对称体问题(结构轴对称＋载荷或约束非轴对称) 在轴对称问题中，通常采用圆柱坐标(r,?,Z)，对称轴为Z轴，半 径方向为r轴，其正方向如图所示，以Z轴为正向的右螺旋转动方向表示?的正向。 图4－1 三角形环单元 空间轴对称问题虽属三维问题，但由于几何形状的轴对称性，在轴对称载荷作用下，所产生的位移、应变和应力与 ? 无关，只是 r和Z的函数，任一点的位移只有两个方向的分量，即沿r方向的径向位移和沿Z方向的轴向位移。由于轴对称， ?方向的位移?等于0，因此轴对称问题是准二维问题，它可以按平面问题处理，但与平面问题不尽相同。 从左式可见，周向（ ? 方向）应变? ?是由径向位移u而引起的，这是与平面问题的重要区别 物理方程 对于各向同性材料，当不考虑温度变化时，轴对称问题的应力与应变之间的物理方程为 由于轴对称问题的位移和应力仅与坐标r,z有关，因此结构离散化只在子午面内进行。离散轴对称体时，采用的单元是环单元，这些环单元与rz平面（子午面）正交的截面可以有不同的形状，例如3结点三角形、6结点三角形或其它形状，各单元在rz平面内形成网格，如图6－1所示。当用3结点三角形环单元进行网格离散时，每一个三角形环单元都有3条棱边，这3条棱边是3个圆周，称为结点圆，而它们与rz面的交点是三角形的3个顶点i,j,m，称为结点。单元结点是圆环形的铰链，三角形环单元之间用这些铰链互相连接传力。 轴对称问题进行计算时，只需取出一个截面进行网格划分和分析，但应注意到单元是环状的，所有的结点载荷都应理解为作用在单元结点所在圆周上。 在轴对称问题分析中，由于回转体承受轴对称载荷，其变形是轴对称的，所以周向位移w=0，r,z方向的位移分量u,v只与坐标位置r,z有关，因此可以依照平面问题，取线性位移函数，即其广义坐标形式为： 1.单元的应变 利用式（4－2）及（4－6）有 从式（4－10）、（4－12）可以看出，单元中的应变分量?x、?y、?xy是常量，周向应变??不是常量，不仅与结点坐标有关，而且与单元内各点的位置(r,z)若有关，同时矩阵B中包含了1/r项，它将给计算带来麻烦。 2.初应变 在轴对称问题中，物体的温度T=T(r,z)是坐标r,z的函数，与坐标?无关。单元内的初应变一般说来不是均匀的，但当单元的数目划分得 较多，即单元的尺寸很小时，可以用一个平均温度来表示。若单元的3个结点温度升高值分别为Ti,Tj,Tm，则单元的平均温度升高值为 对于轴对称问题，三角形环单元可以利用第3章导出的单元刚度矩阵的普遍式写出，如下： 第4章 轴对称体的有限单元法[专题2] 中南大学