ANSYS热仿真：耦合分析与热应力计算方法

热仿真与温度仿真的根本区别在于”耦合”——温度场如何影响结构变形，结构变形又如何反作用于传热。ANSYS中的热-结构耦合分析分为直接耦合和顺序耦合两种路径，选择哪种取决于耦合强度和计算效率要求。

两种耦合策略的物理本质

顺序耦合（Sequential / One-Way）

先做热分析得到温度场T(x,y,z)，然后将温度场作为体载荷导入结构分析，计算热应力。单向耦合假设结构变形不影响温度场——这在大多数工程场景中是合理的，因为热应变通常很小（金属的线膨胀系数α≈10-15×10⁻⁶/°C，100°C温升对应0.1-0.15%应变）。

ANSYS Workbench中的操作流程：

1. 拖拽Steady-State Thermal到项目视图
2. 拖拽Static Structural到Thermal的Solution上，形成数据链接
3. 在Static Structural中，温度场自动作为 Imported Body Temperature载入
4. 添加固定约束和外力载荷
5. 求解热应力

这种方法的优势是计算效率高（两次独立求解），劣势是温度场和应力场不同步，不适用于大变形或接触状态随温度显著变化的场景。

直接耦合（Direct / Two-Way）

在同一个求解器中同时求解温度场和位移场，每一步迭代中温度和位移相互更新。ANSYS用耦合场单元（如PLANE223, SOLID226, SOLID227）实现。

直接耦合适用于：

• 热应变导致几何变化影响传热路径（如间隙闭合）
• 接触状态随温度变化（热接触传导）
• 摩擦生热（高速滑动界面）

计算量通常是顺序耦合的3-10倍，应优先尝试顺序耦合，仅在结果不收敛或物理现象明显双向耦合时切换。

热应力的物理来源

热应力的本质是温度变化导致材料自由膨胀或收缩，但几何约束阻止了自由变形，从而产生内应力。对于各向同性材料，自由热应变为：

ε_thermal = α * ΔT

其中α是线膨胀系数（1/°C），ΔT是温差。如果结构完全自由膨胀，热应变为零（不受约束的热膨胀不产生应力）。应力仅来自约束：

σ = E * α * ΔT * (约束因子)

完全约束时约束因子=1，部分约束时0<约束因子<1。

典型材料的热膨胀系数和100°C温升下的完全约束热应力：

不锈钢304的热应力最高——高膨胀系数+高弹性模量的组合使其在温度梯度下极易产生热应力开裂。

材料属性的温度依赖

真实材料的热物性和力学性能都随温度变化，ANSYS支持定义温度相关的材料曲线：

导热系数的温度依赖

多数金属的导热系数随温度升高而降低。例如不锈钢304：

• 20°C: 16.2 W/m·K
• 200°C: 17.3 W/m·K（略微上升）
• 500°C: 21.0 W/m·K（持续上升）

碳钢则相反，导热系数随温度下降明显：

• 20°C: 60.5 W/m·K
• 500°C: 38.5 W/m·K
• 800°C: 25.0 W/m·K

屈服强度的温度退化

高温下金属屈服强度显著下降，这是热应力分析中必须考虑的：

如果不考虑屈服强度退化，常温下安全系数2.0的零件在500°C时实际安全系数可能降到1.0以下。

实际案例：管道热应力分析

以蒸汽管道为例说明热-结构耦合分析流程：

问题描述

• 管道：Φ219×8 mm，材料20#钢
• 介质温度：350°C
• 环境温度：20°C
• 管道长度：12 m，两端固定
• 保温层：50 mm厚岩棉

热分析

边界条件：

• 内壁：恒温350°C
• 保温层外壁：对流h=8 W/m²·K，T_amb=20°C

稳态温度场结果：

• 管壁内壁温度：350°C
• 管壁外壁温度：348°C（8mm壁厚，金属导热好，温降很小）
• 保温层外壁温度：52°C

结构分析

将温度场导入Static Structural：

• 材料属性：E=186 GPa（350°C值），α=12.5×10⁻⁶/°C，ν=0.3
• 约束：两端面全固定
• 体载荷：Imported Body Temperature（来自热分析）
• 附加载荷：内压4.0 MPa

热应力计算：

• 轴向热应力σ_axial = EαΔT = 186e312.5e-6330 = 767 MPa

这个数值远超材料屈服强度（20#钢在350°C下σ_y≈155 MPa），说明两端固定约束下管道必然发生屈服。实际工程中通过设置膨胀节或U型弯来释放热膨胀，消除过高的轴向热应力。

修改设计——将一端改为滑动支座（允许轴向位移）后重新计算：

• 轴向热应力降至12 MPa（仅由摩擦力产生）
• 最大应力出现在弯头处，约85 MPa，安全系数1.82

这个案例说明，热应力的核心不是计算精度，而是结构设计——通过合理的约束布置和补偿装置，可以从根本上消除热应力风险。

接触面热传导

装配体中两个零件的接触面存在接触热阻，ANSYS用Thermal Contact Conductance（TCC）描述：

q = TCC * (T_1 – T_2)

TCC的典型值：

在热-结构耦合分析中，接触压力会改变实际接触面积，从而影响TCC。ANSYS的纯Thermal分析中TCC是常数，但在耦合场分析中可以用TCC命令的PRESSURE选项定义TCC随接触压力变化的关系，这在螺栓连接热分析中很重要——螺栓预紧力不足会导致TCC偏低，界面温降增大。