📚 课程目录第4章热机械耦合控制方程

4. 热机械耦合控制方程

热应力问题的本质，是温度场与位移场之间通过材料本构和几何约束相互影响。要判断一个工程问题该用“先算温度、再算应力”的弱耦合方式，还是必须把两个场放进同一个方程组联立求解，关键不在于问题听起来复杂不复杂，而在于耦合项是否足以改变温度场或位移场本身。本节从控制方程出发，把力平衡、能量守恒、耦合项写清楚，再给出弱耦合与强耦合的判别依据和选择方法。

一、力平衡方程：位移场的控制方程

在准静态、小变形假设下，忽略惯性项和体积力，力平衡（动量守恒）可写为应力张量的散度为零：

∇·σ = 0

σ = C : (ε - ε_th)
ε = ½(∇u + ∇uᵀ)
ε_th = α(T - T_ref) I

其中 σ 为柯西应力张量，C 为弹性刚度张量，ε 为总应变，ε_th 为热应变，u 为位移矢量，α 为线膨胀系数，T_ref 为无应力参考温度，I 为单位张量。把本构代入平衡方程，得到以位移为未知量的形式：

∇·[C : (½(∇u + ∇uᵀ) - α(T - T_ref)I)] = 0

可以看到，温度 T 以“已知热载荷”的身份出现在位移方程中。若 T 不依赖 u，则位移方程是线性的，温度场可以独立求出后代入。这正是弱耦合的数学根源。

二、能量守恒方程：温度场的控制方程

瞬态热传导的能量守恒方程为：

ρ c_p ∂T/∂t = ∇·(k ∇T) + q_vol + Q_mech

其中 ρ 为密度，c_p 为定压比热容，k 为导热系数，q_vol 为外部体积热源，Q_mech 为力学耗散或变形引起的热源项。在热弹性问题中，Q_mech 主要来自两类效应：一是热弹效应（体积变形引起的可逆温度变化），二是非弹性耗散（塑性、黏弹性、损伤等不可逆过程产热）。

若忽略 Q_mech，温度方程与位移完全解耦，温度场可单独求解。若 Q_mech 不可忽略，温度方程中就出现了位移或应变率，两个场必须联立。这就是强耦合的数学标志。

三、耦合项从哪里来

热机械耦合的耦合项主要有三条通道：

热膨胀通道（T → u）：温度变化产生热应变，直接进入力平衡方程。这是最基本的单向影响，几乎所有热应力问题都存在。
热弹效应通道（u → T）：体积应变率引起温度变化，在快速加载、高约束或大体积构件中才显著。例如快速压缩或冲击下，材料被绝热压缩而升温。
耗散产热通道（非弹性应变率 → T）：塑性功、黏性耗散、摩擦生热等转化为热量。在循环载荷、高速成形、刹车盘摩擦等场景中，这一项往往是强耦合的主导因素。

判断耦合强弱，本质上就是判断后两条通道产生的温度变化，是否足以反过来改变热应变和应力分布。

四、弱耦合与强耦合的判别

工程上常用两个无量纲量来辅助判断。第一个是热弹耦合参数，反映热弹效应强度：

Π_te = E α² T_ref / (ρ c_p)

Π_te 越小，热弹效应越弱。对大多数金属，Π_te 在 10⁻³～10⁻² 量级，常规热应力分析中可忽略。第二个是耗散生热与传导散热之比，可用类似 Peclet 数的形式估计：

Π_diss = (β σ_y ε̇_p L²) / (k ΔT)

其中 β 为 Taylor-Quinney 系数（塑性功转热比例，常取 0.9 左右），σ_y 为屈服应力，ε̇_p 为塑性应变率，L 为特征长度，ΔT 为允许温升。Π_diss 越大，局部温升越显著，越需要强耦合。

判别经验：若温度变化引起的热应力远小于材料屈服强度，且非弹性应变率很低，弱耦合足够；若局部温升超过几十摄氏度、或应变率高于 10⁻¹ s⁻¹、或存在明显塑性耗散，应优先考虑强耦合。

五、弱耦合与强耦合的选择与实现

弱耦合（顺序耦合）的典型流程是：先做瞬态热分析得到 T(x,t)，再把温度场作为热载荷逐时间步映射到力学模型，求解位移和应力。它的优点是计算量小、两个物理场可分别用成熟求解器、便于调试。缺点是忽略了力学对温度的反馈，当耗散生热显著时会低估温度，进而低估热应力和软化。

强耦合（直接耦合）则把位移和温度作为同一方程组的未知量同时求解，典型形式为：

[ K_uu   K_uT ] { Δu }   { F_u }
[ K_Tu   K_TT ] { ΔT } = { F_T }

其中 K_uT 表示温度对力平衡的贡献（热膨胀），K_Tu 表示变形/耗散对能量方程的贡献。强耦合精度高，但矩阵规模大、非线性强、收敛难度高，通常需要 Newton-Raphson 迭代和一致的切线刚度。

选择时可参考以下判断路径：先估计 Π_te 和 Π_diss；若两者都很小，选弱耦合；若 Π_diss 较大或温度与应力强烈相互影响，选强耦合；若介于两者之间，可先用弱耦合试算，检查温升和塑性耗散占比，再决定是否升级为强耦合。

六、一个可操作的判别流程

实际分析时，可以按以下顺序快速判断：

先做纯热分析，得到温度场和最大温升 ΔT_max。
用 ΔT_max 估算热应力 σ_th ≈ EαΔT_max/(1−ν)，与屈服强度比较。
若 σ_th 远低于屈服强度，且无显著非弹性变形，选弱耦合。
若 σ_th 接近或超过屈服强度，检查塑性区应变率和循环次数，估算耗散生热。
若耗散温升与外部温升同量级，或温度变化会显著改变材料参数（如 k、α、σ_y 随 T 变化），选强耦合。
若仍不确定，先用弱耦合试算，提取塑性耗散功率和温升，再决定是否升级。

注意：材料参数随温度变化本身并不构成强耦合。只要温度场不依赖位移，仍可用弱耦合，只需在力学求解时按温度更新材料参数。真正的强耦合来自“位移/应变率反过来改变温度场”的反馈通道。

七、选择建议与常见误区

对于大多数常规热应力问题，如热交换器、电子封装回流、管道热膨胀，弱耦合已经足够，计算效率高且结果可靠。对于高速成形、刹车盘摩擦、激光冲击、焊接熔池附近、循环塑性棘轮等场景，耗散产热和热弹效应不可忽略，应使用强耦合或至少做一次弱耦合-强耦合对比验证。

常见误区有两个：一是把所有热应力问题都做成强耦合，导致计算资源浪费和收敛困难；二是明明存在显著塑性耗散却仍用弱耦合，低估温度后进一步低估热软化和应力松弛，得到偏危险的结论。判别时抓住“温度场是否被力学量显著改变”这一核心，就能在精度和成本之间做出合理选择。