1. 热结构耦合仿真概述

热结构耦合仿真研究的是温度场与位移场之间的相互作用：温度变化引起结构热膨胀或收缩，当这种自由变形受到约束时，结构内部便产生热应力；而应力与变形反过来又可能改变接触状态、间隙与传热路径，从而影响温度分布。理解这条"温度→应变→应力→再影响传热"的闭环，是选择耦合方式与求解器的前提。

1.1 热应力的产生机理

热应力的根源不是温度本身，而是热变形受到约束。对无外部约束、可自由膨胀的各向同性线弹性体，温度由 T₀ 均匀升至 T 时，自由热应变分量为

ε_th = α · ΔT        （ΔT = T − T₀，α 为线膨胀系数）
σ = D : (ε − ε_th)   （D 为弹性刚度张量）

可见只有当总应变 ε 与自由热应变 ε_th 不一致时，应力才非零。约束来源主要有三类：外部约束（如两端固支的管道、被螺栓固定的板）、内部约束（构件内部温度梯度导致相邻材料互相牵制，如厚壁圆筒内外壁温差）、材料不匹配（不同 α 的材料连接，如陶瓷涂层与金属基体、电子封装中芯片与基板）。

以一根两端完全固定的钢杆为例：取 α = 1.2×10⁻⁵ /℃，E = 210 GPa，温升 ΔT = 100 ℃。若允许自由伸长，应变为 1.2×10⁻³；因被完全约束，实际应变为零，故弹性压应变为 −1.2×10⁻³，对应压应力 σ = E·α·ΔT ≈ −252 MPa。这一量级已接近常用结构钢的屈服强度，说明热应力在工程中绝不可忽略。若杆一端自由，则应力为零——这正是"约束决定热应力"的直接体现。

还需注意两类非线性来源：一是材料参数随温度变化（E、α、屈服强度、导热系数 k、比热 c 均随 T 改变），二是高温下进入塑性或蠕变后应力松弛。因此热应力分析常需与温度相关的材料数据，并区分弹性、弹塑性与蠕变本构。

1.2 耦合类型：顺序耦合与完全耦合

按温度场与位移场的求解次序，可分为两大类。

顺序耦合（Sequential / Weakly Coupled）：先解温度场，再把温度作为已知载荷施加到结构分析中求解应力。其成立条件是"单向耦合"——温度显著影响应力，而变形对温度分布的影响可忽略。典型场景是热传导以固体导热为主、接触间隙变化不大、无显著流体对流反馈的问题，如热压成型模具、发动机缸盖热负荷、焊接后的残余应力评估。顺序耦合的优点是两个物理场可分别用最合适的网格与求解器，计算成本低、易调试；缺点是忽略了变形对传热的反馈。

完全耦合（Fully Coupled / Strongly Coupled）：在同一增量步内同时求解温度与位移自由度，温度与变形相互影响。适用于接触热传导随间隙变化、大变形改变辐射角系数、摩擦生热、形状记忆合金相变等强耦合问题。代价是每步需求解更大的耦合方程组，收敛更困难，计算量显著增加。

选择判据可归纳为：若变形引起的传热路径变化对温度场的影响小于工程允许误差（例如间隙变化小于 5%、温升预测偏差可接受），优先顺序耦合；若接触开合、大变形或生热与变形强相关，则必须完全耦合。

1.3 Abaqus 求解器与实现路径

Abaqus 提供多条实现路径，需按耦合强度与问题类型选择：

• Abaqus/Standard + 顺序耦合：先用 *HEAT TRANSFER（或 *COUPLED TEMPERATURE-DISPLACEMENT 的稳态/瞬态传热步）算温度，再以 *TEMPERATURE 读入温度场做 *STATIC 应力分析。适合大多数常规热应力问题。
• Abaqus/Standard 完全耦合：使用 *COUPLED TEMPERATURE-DISPLACEMENT 单元（如 C3D8T、C3D10T、CAX4T），在同一分析中同时求解。适合接触传热、摩擦生热等强耦合。
• Abaqus/Explicit 热-力耦合：用 *DYNAMIC TEMPERATURE-DISPLACEMENT, EXPLICIT，适合高速冲击、大变形、短时瞬态生热问题，无需迭代求解耦合方程。
• Abaqus/CFD 与共轭传热：涉及流体对流换热时，可用 CFD 求解流场温度，再与结构做顺序或协同耦合。

关键设置包括：定义随温度变化的 conductivity、specific heat、density、elastic、expansion；设置热边界（*FILM 对流、*RADIATION 辐射、*BOUNDARY 指定温度）；在完全耦合中为接触对指定 *GAP CONDUCTANCE 与 *GAP HEAT GENERATION。初始温度用 *INITIAL CONDITIONS, TYPE=TEMPERATURE 给定，参考温度（零热应变温度）由 *EXPANSION 中的 ZERO 参数或分析步的 temperature 设定，二者必须与建模意图一致，否则热应变基准会出错。

1.4 典型工程应用场景

热结构耦合在以下场景中最为常见：电子封装与功率器件（芯片-基板-焊点因 α 不匹配产生翘曲与焊点疲劳）、航空发动机热端部件（涡轮叶片、燃烧室在温度梯度下的低周疲劳与蠕变）、汽车制动与排气系统（制动盘热裂纹、排气歧管热变形）、核电与压力容器（启停热冲击下的热应力与断裂评定）、增材制造（逐层熔凝产生的温度梯度与残余应力、变形控制）、焊接与热处理（焊接残余应力与变形预测）。这些场景的共同点是：温度场不均匀、存在约束或材料不匹配、且热应力与变形直接影响寿命或装配精度。

1.5 建模前的判断清单

为便于实际选型，可将上述内容整理为一份判断清单，在建模前逐项确认：

需要强调的是，顺序耦合与完全耦合并非互斥，工程中常先用顺序耦合快速评估温度场与热应力量级，再对关键部位或关键工况切换为完全耦合进行精细分析。无论采用哪种方式，都必须保证参考温度、初始温度与热边界条件三者一致，并检查材料参数的温度相关性是否覆盖实际工作温区，否则再精细的网格也无法给出可信的热应力结果。