Ansys热变形仿真：热固耦合里的温度路径和约束刚度是一对互相伤害的变量

热变形分析的计算框架很清晰：先算温度场（稳态或瞬态），再把温度场作为热载荷喂给结构分析，算热应变和应力。Ansys Workbench里两条系统拖在一起，一条线连过去，看起来一气呵成。

但”先算温度场”这一步有个隐藏的抉择——你用稳态温度场还是瞬态温度场？稳态认为热源和散热已经达到平衡，整个结构有一个固定的温度分布。瞬态考虑升温过程——冷态启动→温度沿时间和空间演化→达到稳态。

对于排气歧管，稳态和瞬态的应力差别有多大？这个项目里，差别接近50%。

稳态和瞬态的温度分布不是同一张图

排气歧管在发动机冷启动后的前120秒是急剧升温阶段——排气温度从常温跳到约850°C（全负荷），热气流以约150m/s的速度冲刷歧管的内壁面。外壁面通过自然对流和辐射向发动机舱散热。

如果只做稳态热分析，假设歧管壁已经达到热平衡——内外壁温差按照壁面厚度(4mm)+导热系数(40W/m·K)+换热系数(内壁350W/m²·K,外壁15W/m²·K)算出来约85°C。这个温差在稳态是对的。

但在瞬态，前60秒内，内壁面温度已经爬到了620°C，外壁面才到380°C——瞬态内外壁温差约240°C，是稳态的近3倍。这个240°C的瞬态温差造成的热弯曲力矩远大于稳态的85°C。

Ansys的瞬态热分析给出了在40秒时歧管内壁温度650°C、外壁380°C的分布——对应热弯曲应力约280MPa。而用稳态温度场算出来的最大热应力只有195MPa。差了85MPa。对于铸铁歧管（HT250，屈服强度250MPa），195MPa在许用范围、280MPa超了。

冷启动热冲击才是排气歧管热变形分析的最危险工况——不是稳态全负荷。

约束超静定效应

热变形仿真的另一个坑是约束。热膨胀本身不产生应力——只在膨胀被约束的时候才产生热应力。排气歧管通过法兰螺栓固定在缸盖上，螺栓约束了歧管的自由热膨胀。

这个约束是弹性的（螺栓有刚度），不是刚性的。如果歧管法兰和缸盖之间的每个螺栓都设成固定约束（全自由度约束），歧管在热膨胀下无处释放——管段的压缩热应力被系统性地高估。

正确的做法是把缸盖也建入模型中——歧管螺栓连接到缸盖模型上，缸盖本身在受热膨胀（铝合金缸盖的热膨胀系数~21×10⁻⁶/K，高于铸铁歧管的~12×10⁻⁶/K）。缸盖的热膨胀方向与歧管膨胀方向相反，两个膨胀互相部分抵消。

加上缸盖模型后，歧管法兰处的热应力从350MPa降到了265MPa。不加缸盖就把约束从弹性变成了刚性——多算了约85MPa的虚假约束应力。

蠕变松弛在高周次热循环中的累积

排气歧管在整个寿命期内经历几千次冷-热-冷循环。铸铁在500°C以上有明显的蠕变——即恒载下的持续变形。每次热循环产生的塑性应变在下次循环中不再能恢复——这就是热机械疲劳（TMF）的损伤机理。

这个项目在Ansys中做了简单的粘塑性分析——用Norton蠕变模型（ε_dot = A·σ^n·exp(-Q/RT)），铸铁的蠕变参数A=2.3×10⁻⁸ (MPa⁻n·h⁻¹)、n=3.2、Q=285kJ/mol。跑了10个热循环，歧管分支处的累积塑性应变已经达到了0.8%。在几千次循环后，这个量会累积到几个百分点——对应低周疲劳寿命的快速消耗。

Ansys热变形仿真不是”把温度场算出来再贴到结构上”这么简单。温度场是瞬态还是稳态、约束从哪一级开始建、是否考虑材料的热蠕变松弛——三个选择差出上百MPa的应力、改变了安全与否的判定边界。