2. 热传导与热应力基础

热力耦合仿真的第一步，不是急着在 Abaqus 里点“Coupled Temperature-Displacement”，而是先把四个物理量之间的关系理清楚：热流怎么走、材料受热膨胀多少、膨胀被约束后产生多大应力、这些量在 Abaqus 里分别由哪些参数和边界条件控制。本节围绕这条主线展开，并给出可直接对照的示例参数。

一、傅里叶定律：热流的方向与大小

傅里叶定律描述导热的基本规律：热流密度矢量 q 与温度梯度 ∇T 成正比、方向相反，即

q = -k · ∇T        （一维形式：q_x = -k · dT/dx）

其中 k 为导热系数（W/(m·K)），负号表示热量自发从高温流向低温。三维各向同性材料中，Abaqus 只需输入一个 k；各向异性材料则需给出 k11、k22、k33 及剪切项。稳态无内热源时，控制方程退化为 ∇·(k∇T)=0；瞬态问题还要引入密度 ρ 与比热 c，即 ρc·∂T/∂t = ∇·(k∇T) + Q。

举例：一块 100 mm 厚的钢板，两侧温度分别为 200 ℃ 与 20 ℃，k=45 W/(m·K)，一维稳态下温度梯度为 (20−200)/0.1 = −1800 K/m，热流密度 q = −45×(−1800) = 81000 W/m²，方向由高温侧指向低温侧。这个量级说明：温差不大时，厚壁件内部的热流依然可观，不能凭直觉忽略。

二、热膨胀系数与热应变

热应变描述“温度变化引起的自由变形”，其定义为

ε_th = α · (T − T_ref)

α 为热膨胀系数（1/K），T_ref 为应力自由参考温度。Abaqus 中 α 在材料属性 Expansion 下定义，参考温度在 Expansion → Reference temperature 或分析步的初始温度场中给定。注意两点：一是 α 随温度变化时应以表格形式输入，Abaqus 会做插值；二是 T_ref 必须与装配时的“无应力温度”一致，否则整个模型会凭空产生初始热应变。

示例：铝合金 α=23×10⁻⁶ /K，T_ref=20 ℃，升温到 120 ℃，则 ε_th = 23×10⁻⁶×100 = 2.3×10⁻³，即 0.23%。一根 500 mm 长的杆若两端自由，伸长量 ΔL = 0.23%×500 = 1.15 mm。这个数值在精密装配中往往就是间隙与过盈的分界线。

三、热应力是怎么产生的

关键判断只有一句：热应变本身不产生应力，只有当热变形被约束时才产生热应力。约束有三种典型来源：外部约束（如两端固支）、内部约束（材料 α 不一致或温度分布不均）、以及结构超静定带来的自约束。

一维全约束杆的经典公式为 σ = −E·α·ΔT/(1−ν)（平面应变）或 σ = −E·α·ΔT（单轴）。以钢为例，E=210 GPa，α=12×10⁻⁶ /K，ΔT=100 ℃，单轴全约束下 σ = −210×10³×12×10⁻⁶×100 ≈ −252 MPa，已接近普通碳钢屈服强度。这说明：几百摄氏度的温差在强约束下足以让构件屈服，热应力绝不是小量。

四、Abaqus 中的基本设置

在 Abaqus/Standard 中做热力耦合，常用两条路径：一是 Coupled Temperature-Displacement 单元（如 C3D8T），温度与位移同时求解；二是先做 Heat Transfer 分析得到温度场，再以 Predefined Field 导入做 Static, General 热应力分析。前者适合强耦合、接触传热问题，后者计算成本低，适合热应力主导、温度场不受变形影响的场景。

关键设置清单如下：

• 材料：Conductivity（k）、Density（ρ）、Specific Heat（c）、Elastic（E、ν）、Expansion（α 与参考温度）。瞬态分析缺 ρ 或 c 会直接报错。
• 分析步：Coupled Temp-Displacement 需选 Transient 或 Steady-state；瞬态要设初始增量步、最小/最大增量步与总时间。
• 边界条件：温度侧用 Temperature 或 Surface Film Condition（对流换热系数 h 与环境温度）；位移侧按实际约束施加。
• 初始温度场：用 Predefined Field → Temperature 指定，默认 0 ℃ 是常见错误来源，务必改为实际装配温度。
• 网格：热梯度大的区域（如热源附近、薄壁过渡处）需局部加密；C3D8T 单元在弯曲主导区域易剪切自锁，可改用 C3D10T 或开启增强沙漏控制。

下面这张图把“温度场 → 热应变 → 约束 → 热应力”的因果链，以及 Abaqus 中对应的设置入口放在同一张关系图里，便于对照检查。

五、一个可复现的示例与结果判读

取 200 mm×20 mm×20 mm 的 45 钢方杆，一端固定、另一端自由，初始温度 20 ℃，整体升温至 120 ℃。材料参数：E=210 GPa，ν=0.3，α=12×10⁻⁶ /K，k=45 W/(m·K)，ρ=7850 kg/m³，c=480 J/(kg·K)。

若只施加温度场、不施加位移约束，Abaqus 结果中杆会整体伸长 ΔL = 12×10⁻⁶×100×200 = 0.24 mm，Mises 应力接近 0，验证“自由膨胀无热应力”。若把自由端也改为固支，则杆处于全约束状态，理论热应力 σ = −E·α·ΔT = −210×10³×12×10⁻⁶×100 = −252 MPa，Abaqus 中沿轴向的 S11 应接近该值（受端部效应影响，端部附近会有应力集中）。

判读时重点看三处：一是 NT11 温度场是否与预期分布一致；二是 U 位移是否等于自由膨胀量；三是 S 应力是否与 σ=−E·α·ΔT 量级吻合。若应力远小于理论值，多半是约束没加对或参考温度设错；若应力远大于理论值，则要检查是否误加了多余约束或接触定义过刚。

把以上四步串起来：傅里叶定律决定温度怎么分布，α 与 T_ref 决定自由热应变多大，约束决定这些应变有多少被“憋”成应力，Abaqus 中的材料、分析步、边界条件与初始温度场则分别对应这四个环节。任何一步设置错位，结果都会偏离物理真实。