📚 课程目录第3章材料热物理与热力学基础

3. 材料热物理与热力学基础

热循环环境下，材料性能退化的根源往往不在单一力学载荷，而在于温度场与应力场之间的耦合。要建立可靠的热-力耦合分析，第一步不是写方程，而是把材料的热物理参数体系理清楚：哪些参数控制温度分布，哪些参数控制热应变，哪些参数在温度变化时发生非线性突变。本节围绕五个核心参数——热膨胀系数、热导率、比热容、热扩散率和相变潜热——说明它们的物理含义、温度依赖性以及如何组织成可用的分析参数体系。

一、五个参数的物理分工

可以把热-力耦合问题拆成两条链。第一条是“热链”：热量如何传入、如何在材料内部重新分布。这条链由热导率 λ、比热容 c_p 和密度 ρ 共同决定，三者组合出热扩散率 α = λ/(ρc_p)，它直接决定温度场随时间的推进速度。第二条是“力链”：温度变化如何转化为应变和应力。这条链由热膨胀系数 α_L 控制，在受约束条件下产生热应力 σ = Eα_LΔT/(1−ν)（弹性、单轴约束的简化形式）。相变潜热 L 则是一个“隐藏项”，它不改变稳态温度分布，却会在相变区间吸收或释放大量热量，使温度曲线出现平台或拐点，从而改变瞬态热应力的时间历程。

需要特别强调参照约定：本文热膨胀系数均指线膨胀系数，定义为 α_L = (1/L₀)(dL/dT)，单位 1/K 或 10⁻⁶/K；热应变 ε_th = α_L(T)(T − T_ref)，其中 T_ref 为无应力参考温度，全文与图示保持一致。若改用体膨胀系数，需明确 α_V ≈ 3α_L（各向同性材料）。

二、温度依赖性：为什么不能只取一个常数

工程手册常给出室温单值，但热循环跨度可能从 −60 ℃ 到 1000 ℃ 以上，此时参数随温度的变化不可忽略。典型规律如下：

热膨胀系数 α_L(T)：多数金属随温度升高而增大，且在相变点附近出现阶跃；例如某低合金钢在 20 ℃ 约 11.5×10⁻⁶/K，到 600 ℃ 可升至约 14×10⁻⁶/K。若用室温值外推高温，热应变会被系统性低估。
热导率 λ(T)：金属中电子导热占主导，温度升高使电子散射增强，λ 通常下降；陶瓷与玻璃则可能随温度升高而略升。同一材料在 20 ℃ 与 800 ℃ 的 λ 可相差 30% 以上。
比热容 c_p(T)：低温段随 T³ 规律上升，接近德拜温度后趋于平缓；相变时出现尖峰，因为潜热被折算进等效比热。
热扩散率 α(T)：由 λ、ρ、c_p 组合而成，金属通常随温度升高而下降，是瞬态热分析中决定“热穿透深度”的关键量。
相变潜热 L：在相变区间内，可用等效比热 c_eff = c_p + L/(ΔT_相变) 处理，把潜热“摊平”到相变温区，便于数值求解。

三、参数体系如何组织成分析输入

建立热-力耦合分析时，建议按“温度场求解—热应变计算—应力求解”的顺序组织参数，而不是把所有参数堆在一张表里。温度场求解需要 λ(T)、ρ(T)、c_p(T)（或直接给 α(T)）以及相变潜热；热应变计算需要 α_L(T) 与参考温度 T_ref；应力求解还需要弹性模量 E(T)、泊松比 ν(T) 以及必要的屈服与蠕变参数。下面这张图给出参数在两条链中的位置与流向，可作为建模前的参数清单核对图。

四、教学示例：从参数表到热应力估算

下面用一个简化示例说明参数如何真正参与计算。假设某耐热合金构件在 20 ℃ 装配（T_ref = 20 ℃），服役时局部升至 620 ℃，构件受刚性约束无法自由膨胀。已知该合金在 20～620 ℃ 区间的平均线膨胀系数 ᾱ_L = 13.8×10⁻⁶/K，620 ℃ 时 E = 155 GPa，ν = 0.31。则：

ΔT = 620 − 20 = 600 K
ε_th = ᾱ_L · ΔT = 13.8e-6 × 600 = 8.28e-3
σ = E · ε_th / (1 − ν)
  = 155e9 × 8.28e-3 / (1 − 0.31)
  ≈ 1.86e9 Pa ≈ 1860 MPa

这个量级已远超多数耐热合金 620 ℃ 的屈服强度（示例中假设约 400～500 MPa），说明刚性约束下热应力足以驱动塑性变形甚至开裂。若改用室温 α_L = 11.5×10⁻⁶/K 估算，ε_th 仅 6.9×10⁻³，σ 约 1550 MPa，低估约 17%。这正是必须使用温度相关参数的直接理由。

注意：上述数值为教学示例中的合理假设，并非某具体牌号的实测数据。实际分析中应使用材料手册或实测的 α_L(T)、E(T)、ν(T) 曲线，并明确 T_ref 的选取依据。若构件存在温度梯度，还需先由热链求出 T(x,t)，再逐点计算热应变，不能只用单一 ΔT。

五、相变潜热的处理要点

当材料在热循环中发生固-固相变、熔化或凝固时，潜热会使温度曲线在相变温区出现“平台”。数值处理上有两种常用方式：一是等效比热法，把 L 折算进 c_p，在相变区间设置一个高峰值，实现简单但需控制时间步长以免跳过相变区；二是焓法，直接以焓 H(T) 作为状态变量，H 在相变区连续但斜率突变，数值稳定性更好。无论哪种方式，都必须保证相变温区 ΔT_相变与潜热 L 的取值来自同一数据源，避免混用不同测试条件的数据。

把上述参数整理成分析输入时，建议按温度节点逐点给出，而不是只给一个区间平均值。例如在 20、100、200、400、600、800 ℃ 六个节点上分别列出 λ、c_p、ρ、α_L、E、ν，中间温度由线性或样条插值。这样既保留了温度依赖性，又便于在有限元软件中直接以表格形式输入。参数体系的完整性，最终决定了热-力耦合分析结果是否可信。