基于ANSYS的冷罩外壳有限元热应力分析.pdfVIP

第 12期 张仁俭等 ：基于 ANSYS的冷罩外壳有限元热应力分析 太阳辐射 2 热应力耦合分析 2．1 创建几何模型 该几何模型的建立是依据某 型高速飞行器的冷 罩的尺寸及实际结构进行很少的必要的简化而建立 的，造型去掉了连接 的形式。而冷罩的内壳与高速 压力低—～——— = ≥呈 三 制冷剂 飞行器的头部通过绝热层紧贴 ，又 由于高速飞行器 的头部的刚性足够大，因此 ，在温度变化 的情况下 ， 图 1 冷罩工作示意图 冷罩内壳的变形是可以忽略 的。因此 ，笔者可以不 考虑冷罩的内壳变形 ，而 只对冷罩 的外壳进行分 1 热应力方程的建立 析 。 物体由于热膨胀只产生线应变 ，剪切应变为零。 外壳受到隔梁 的支撑作用 ，隔梁与 内外壳相连 这种 由于热变形产生的应变可以看作是物体的初应 接，故在该模型 中包括了隔梁和 内罩。内壳跟外壳 变。计算热应力时只需计算 出热变形引起的初应变 通过隔梁连接在一起。隔梁在 内外壳之间是点状均 。 ，求得相应 的初应变引起 的等效结点载荷 P (简 匀分布的，由于视角关系，这里只能看到其中一部分 称温度载荷)，然后按通常求解应力一样解得 由于 的隔梁。其建模的结果如图2所示。 热变形引起的结点位移口，然后可 以由。求得热应力 。 也可以将热变形引起的等效结点载荷JP 与其 隔梁 它载荷项合在一起，求得包括热应力在 内的综合应 力。计算应力时应包括初应变项 = D(8— 0) (1) 式中：s。是温度变化引起温度应变 ，它现在是作为 初应变出现在应力应变关系式 中，对于三维问题是 占。= ( 一咖。)[1，1，1，0，0，0] (2) 式中： 是材料的热膨胀系数 (1／~C)；咖。是结构 的 初始温度场。 虚位移原理的表达式为 f(占 一8u )d—f8u于ds=0 (3) J JS 将式(1)代人式(3)中，则可得到包含温度应变 图2 三维模型图 在内，用以求解热应力问题的最小位能原理，它的泛 函表达式如下 2．2 材料参数的确定 查阅材料手册 ，铝合金的密度为2700kg／m’，铝 ： fo(Ds一占D占。一ud力一frurdr 合金的导热系数取 204W／(m ·oI=)，铝合金的比热 (4) 容为880J(kg·c=【)，铝的弹性模量为7X10加；泊松 将求解域 进行有限元离散，从 an．=0将可 比为 0．3；导热系数为 205；线膨胀系数为 2．4× 得到有限元求解方程 10～ 。 KA =P (5) 2．3 单元类型和网格的划分 与不包含温度应变的有限元求解方程相区别的