1. 热结构耦合分析概述

在工程结构中，温度变化几乎从不孤立存在。发动机缸体在冷启动与全负荷之间经历上百摄氏度的温差，电子封装在回流焊冷却后残留翘曲，管道在输送高温介质时因热膨胀受限而产生轴向推力。这些现象的共同点是：温度场先改变结构的应变状态，应变若受到约束又反过来影响应力，而应力与变形在多数情况下并不会显著改变温度场本身。理解这种“单向主导、弱反馈”的物理链条，是进入热结构耦合分析的第一步。

热应力产生的机理

热应力的根源可以归结为一句判断：热变形受到约束，才会产生应力。若一个各向同性线弹性体可以自由膨胀，温度由 T₀ 升至 T，其自由热应变为 ε_th = α(T − T₀)，其中 α 为线膨胀系数。此时物体内部不产生应力，只是整体尺寸发生变化。一旦这种自由变形被外部边界（如两端固定的管道）或内部约束（如温度分布不均匀导致相邻材料互相牵制）所限制，约束应变 ε_mech = ε_total − ε_th 就非零，应力随之出现：σ = D·(ε − ε_th)，D 为弹性矩阵。

由此可以区分两类典型来源。第一类是外部约束型：构件整体温度均匀升高，但支座、焊缝或相邻构件阻止其伸长。一根两端固定的钢杆，α = 1.2×10⁻⁵ /℃，温升 ΔT = 100 ℃，若完全约束，则热应力 σ = −EαΔT ≈ −2.0×10⁵ × 1.2×10⁻⁵ × 100 = −240 MPa，已接近普通碳钢屈服强度，这就是管道热应力必须设置膨胀节的直接原因。第二类是内部约束型：构件温度分布不均，表层与芯部、厚壁与薄壁之间膨胀量不同，即使没有外部固定也会产生自平衡应力。淬火冷却时表面先收缩、芯部后收缩，表面受拉、芯部受压，正是残余应力的来源。

需要特别强调参照约定：本文及后续图示统一采用拉应力为正、压应力为负，温度以无应力参考温度 T_ref 为基准，ΔT = T − T_ref。当 ΔT 为正且变形被完全约束时，应力为压应力（负值）；当 ΔT 为负（冷却）且收缩被约束时，应力为拉应力（正值）。这一符号约定在后续所有算例中保持一致。

耦合场分类与求解策略

热与结构的相互作用按信息传递方向可分为两类。顺序耦合（弱耦合）：先求解温度场，再把温度作为体载荷施加到结构模型上求解应力，温度场不受结构变形影响。绝大多数工程热应力问题属于此类，因为热变形对温度分布的反馈通常很小。直接耦合（强耦合）：温度与位移在同一方程组中同时求解，适用于压电热、热-流-固、摩擦生热等双向影响显著的场景，计算代价高但物理更完整。

在 ANSYS 体系中，顺序耦合又有多种实现路径。经典 ANSYS 可用热单元求解温度后通过 ETCHG 转换为结构单元，或用物理环境法（PHYSICS）在同一模型中切换；Workbench 平台则通过“热分析→结构分析”的数据链接自动传递温度结果，操作直观、适合参数化与多方案对比。选择哪种路径，取决于模型规模、是否需要参数扫描、以及是否涉及非线性材料。

典型工程应用场景

热结构耦合分析的应用面很宽，但可以按“温度来源—约束形式—关注结果”三个维度快速归类。发动机缸体与缸盖：燃气侧高温、冷却水侧对流，温度梯度导致缸盖鼻梁区产生压应力、冷却后转为拉应力，关注疲劳寿命与密封面变形。电子封装：回流焊冷却过程中芯片、基板、焊球热膨胀系数不匹配，关注焊点塑性应变与翘曲量。高温管道与换热器：管壁内外温差与端部约束共同作用，关注轴向应力是否超过许用值、是否需要膨胀节。涡轮盘与叶片：离心载荷叠加高温蠕变，属于热-结构-材料非线性耦合。焊接与增材制造：移动热源引起瞬态温度场，冷却后形成残余应力与变形，通常需要顺序耦合的瞬态分析。

这些场景的共同分析流程可以概括为：建立几何与网格 → 施加温度边界（温度、对流、热流、辐射）求解温度场 → 将温度作为体载荷映射到结构模型 → 施加位移约束与机械载荷 → 求解应力与变形 → 用许用应力、疲劳曲线或变形公差判据评估结果。其中最容易出错的两个环节是参考温度的设定和温度载荷的映射精度：参考温度若与工艺实际不符，热应变整体偏移；热网格与结构网格不匹配时，映射误差会直接污染应力结果。