1. 热力耦合仿真工程导论
在工程结构中，温度很少是一个孤立变量。发动机缸体受热后膨胀，却受到螺栓与缸盖的约束；电子封装在回流焊冷却过程中，因芯片、基板与焊球的热膨胀系数不同而产生翘曲与应力；刹车盘在摩擦生热后表面温度骤升，又被内部低温区域约束。这些现象的共同点是：温度场会改变结构的变形与应力，而变形与接触状态又会反过来改变传热路径。当这种相互作用不可忽略时，就必须做热力耦合分析，而不能先算温度、再把温度当载荷单向加到结构上。
一、什么是热力耦合
热力耦合（Thermo-Mechanical Coupling）指温度场与力学场之间存在双向或单向的相互影响。其物理基础可归纳为三条：
- 热应变：温度变化 ΔT 引起热应变 ε_th = α·ΔT，α 为热膨胀系数。若结构受约束，热应变无法自由释放，就转化为应力。
- 材料参数随温度变化：弹性模量 E、屈服强度 σ_y、导热系数 k、比热 c、密度 ρ 通常都是温度的函数，高温下 E 和 σ_y 明显下降。
- 几何与接触状态影响传热：接触压力越大，接触热导越高；间隙、翘曲会改变辐射与对流路径。这就是力学场对温度场的反向作用。
按耦合强度，工程上通常分为两类。若力学变形对温度场影响很小（如多数稳态热应力问题），采用顺序耦合：先做传热分析得到温度场，再把温度作为预定义场导入应力分析。若变形、接触、摩擦生热显著改变热边界（如刹车、大变形接触、金属成形），则需完全耦合，在同一求解步中同时求解温度与位移自由度。
二、工程意义与典型场景
热力耦合分析的价值在于回答三类问题：结构在温度载荷下是否安全（应力、塑性应变、蠕变是否超限）；热变形是否影响功能（间隙、密封、光学对准、装配精度）；热-力循环是否导致寿命问题（热疲劳、焊点失效、热机械疲劳裂纹）。典型应用包括：
| 场景 | 关注量 | 耦合方式 |
|---|---|---|
| 电子封装回流焊/热循环 | 焊球塑性应变、翘曲 | 顺序耦合为主 |
| 刹车盘摩擦制动 | 表面温度、热应力、热裂纹 | 完全耦合（摩擦生热） |
| 发动机缸盖热负荷 | 热应力、低周疲劳 | 顺序耦合 |
| 金属热成形/淬火 | 相变、残余应力、变形 | 完全耦合 |
三、Abaqus 的求解能力
Abaqus 提供完整的温度-位移耦合求解体系。核心单元类型为耦合温度-位移单元（如 C3D8T、C3D10T、CAX4T 等），每个节点同时具有温度自由度（11）和位移自由度（1、2、3）。分析类型主要有：
- Coupled Temperature-Displacement（完全耦合）：适用于摩擦生热、大变形接触、热-力强耦合。可选用隐式（Standard）或显式（Explicit）求解。
- Heat Transfer + Static/Visco：顺序耦合的常用组合，先算温度再导入应力，计算成本低。
- Steady-State / Transient：稳态用于长期热平衡，瞬态用于升温、冷却、热冲击过程。
材料方面，Abaqus 支持随温度变化的弹性、塑性、蠕变、超弹性等本构，并可通过表格输入 E(T)、σ_y(T)、α(T)、k(T)、c(T)。接触方面，可定义接触热导（Gap Conductance）、摩擦生热分配系数、辐射视角因子等，使热边界随力学状态动态变化。
四、课程实战路线图
本课程围绕“能建模、能算准、能判断”三条主线展开，从基础概念逐步推进到复杂工况。整体路线可概括为：热力耦合基础与单元选择 → 材料参数随温度定义 → 稳态/瞬态传热分析 → 顺序耦合与完全耦合实现 → 接触热导与摩擦生热 → 热应力与热疲劳评估 → 复杂工况综合案例。下图给出课程的知识结构与求解路径关系。
学习本课程时，建议始终带着一个具体工程对象来练习，例如一块受热约束的金属板、一个含焊球的封装结构或一个摩擦副。每学完一个模块，就回到该对象上增加一层物理：先只算温度，再加入热膨胀，再引入温度相关材料，最后打开接触与摩擦生热。这样能清楚看到每一层物理对结果的影响量级，也能判断何时必须用完全耦合、何时顺序耦合已经足够。