数值模拟与仿真基础——4 轴对称问题的有限元法.pptVIP

数值模拟与仿真技术基础 麻向军 华南理工大学 聚合物新型成型装备国家工程研究中心 聚合物成型加工工程教育部重点实验室 轴对称问题的有限元法 5.1 轴对称体的离散化 5.2 位移函数 5.3 单元的应变和应力 5.4 单元刚度矩阵 5.5 等效节点载荷 5.6 结构的总体刚度矩阵 5.7 应力计算 5.1 轴对称体的离散化 轴对称体问题的位移和应力仅与坐标 r、z 有关，因此，结构离散化只在子午面内进行。 在子午面内虽然是三角形单元，实际上代表一个 “三棱环体”单元，称为三角形环单元。每一个三角形环单元都有三条棱边，这三条棱边是 3 个圆周，它们与 rz 面的交点是三角形的 3 个顶点 i,j,m。 3 条棱边称为节点圆，三个顶点 i,j,m 称为节点。单元节点是圆环形的铰链，并相互连接和传力。 5.2 位移函数 取线性位移函数： 写成矩阵形式为： 5.3 单元的应变和应力 式中：S —— 应力矩阵应力子矩阵 5.4 单元刚度矩阵 式中子矩阵 Kts 为2 ? 2阶矩阵。与平面问题不同的是矩阵B内含有元素。 子刚度矩阵近似计算 为了简化计算，引入三角形形心坐标： 子刚度矩阵的精确积分 将矩阵 Bt（Bs）分成与坐标 r、z 无关的常值部分 和与坐标 r 有关的变值部分 Bt’（Bs’），即： 对于给的三角形单元而言， 是个常量矩阵，与坐标 r、 z 无关。 式中： 5.5 等效节点载荷 P——分布在圆周上的轴对称载荷， r——集中力作用点的半径。 故等效节点载荷也是节点圆周上的分布力。 对该式可采用数值积分，即可求得温度变化引起的等效节点载荷。 如果采用近似计算， 即取 则有 谢谢大家！ z r O i i j j m m 可求出式中的待定系数： 式中：A——三角形环单元的截面积 系数 同样可以得到平面问题类同的位移函数表达式： 上式中形函数 1、单元的应变 单元中的应变分量 ?r 、 ?z 和?rz为常量，而周向应变不是常量。不仅与节点坐标有关，而且还与单元内各点的位置有关，同时包含的项给计算带来了麻烦。 2、单元的应力 应用普遍公式： 应用载荷移置的普遍公式： 1、集中力： 2 .体积力 3 .离心力 4 .变温应力 由于温度变化，体内将产生应变。可被作为初应变来处理。假如温度改变为T，产生的初应变为： 将各矩阵代入 5、表面力 α i j q m z r O * *