2. 热弹性理论基础

热弹性理论要回答的核心问题只有一个：当温度改变时，材料内部为什么会产生应力，以及这种应力如何用可计算的量表达。理解这条主线，后续的热冲击、热疲劳、热屈曲才有落脚点。本节从热应变的定义出发，建立本构方程，再讨论热膨胀系数随温度变化时的工程处理方法。

一、热应变：自由膨胀与约束膨胀

先约定符号：以拉应变为正、压应变为负，温度升高为正。设材料在参考温度 T₀ 时长度为 L₀，当温度均匀升至 T 时，若试件完全不受约束，其长度变为 L，则定义线热应变

ε_th = (L − L₀)/L₀ = α·(T − T₀)

其中 α 为线热膨胀系数，单位 1/K 或 1/℃。这里的关键在于：热应变是一种“自由应变”，它本身不产生应力。一根可以自由伸缩的杆，无论升到多高温度，内部应力都为零。只有当膨胀受到约束、或相邻部分膨胀量不一致时，热应变才被“挡住”，转化为弹性应变，进而产生热应力。

因此热应力的本质是：总应变中扣除热应变后剩下的弹性部分，才对应应力。这一句话是全部热弹性本构关系的出发点。

二、本构方程：热弹性胡克定律

在一维情形下，总应变 ε 由弹性应变与热应变叠加：

ε = σ/E + α·(T − T₀)

反解得到应力表达式：

σ = E·[ε − α·(T − T₀)]

这个式子说明：如果总应变 ε 被约束为零（完全约束杆），则 σ = −E·α·ΔT，负号表示升温时受压。以常见钢材为例，取 E = 210 GPa、α = 12×10⁻⁶ /K，若 ΔT = 100 K 且完全约束，则 σ = −210×10⁹ × 12×10⁻⁶ × 100 ≈ −252 MPa。这个量级已接近普通碳钢屈服强度，说明不大的温升在强约束下就能产生可观应力，这正是工程中必须重视热应力的原因。

推广到三维各向同性热弹性体，本构方程为：

ε_ij = (1+ν)/E · σ_ij − (ν/E)·σ_kk·δ_ij + α·ΔT·δ_ij

其物理含义是：正应变由三向应力引起的弹性变形与各向同性热膨胀共同组成，而剪应变不受温度影响（各向同性材料热膨胀不产生剪应变）。工程上常把上式改写为应力形式，用于有限元求解：

σ_ij = 2G·ε_ij + λ·ε_kk·δ_ij − (3λ + 2G)·α·ΔT·δ_ij

其中 G 为剪切模量，λ 为拉梅常数。可以看到温度项以“等效初应变”的形式进入，这也是有限元热应力分析中把温度载荷转化为等效节点力的理论依据。

三、热膨胀系数与温度相关性

实际材料的 α 并非常数。以奥氏体不锈钢为例，20～100 ℃ 时 α 约为 16.5×10⁻⁶ /K，而 20～600 ℃ 的平均值可升至约 18.5×10⁻⁶ /K；铝合金、钛合金、陶瓷的 α 随温度变化规律各不相同。若仍按常温常数计算，在高温工况下热应变可能被低估 10% 以上。

工程处理有两条路径。第一条是平均热膨胀系数法：直接查取从参考温度 T₀ 到目标温度 T 的平均系数 ᾱ(T)，用 ε_th = ᾱ(T)·(T − T₀) 计算总热应变。这种方法简单，适合手算和初步评估，但要求参考温度与数据表一致。第二条是瞬时系数积分法：若已知瞬时系数 α(T)，则

ε_th = ∫_T₀^T α(T′) dT′

在有限元中通常直接输入随温度变化的 α(T) 表或热应变—温度曲线，由求解器积分。需要注意：平均系数与瞬时系数不可混用，混用会引入系统性误差。当温度跨度大、且材料经历多次升降温时，还应考虑参考温度的选取——通常取无应力状态对应的温度作为 T₀。

四、约束程度与热应力的关系

实际结构很少是完全自由或完全约束，多数处于部分约束状态。设杆件两端由刚度为 k 的弹性约束连接，自由热伸长量为 δ_th = α·ΔT·L，则实际被约束住的变形为 δ_th·k/(k + EA/L)，对应应力为

σ = −E·α·ΔT · [k/(k + EA/L)]

括号内即约束度，取值 0～1。k→0 时应力为零（自由膨胀），k→∞ 时趋于完全约束值 −EαΔT。这一关系说明：降低约束刚度是控制热应力的有效手段，工程中的膨胀节、滑动支座、柔性连接正是基于此原理。

五、工程处理要点与常见误区

把上述理论落到实际计算，需要把握几个要点。第一，明确参考温度。热应变是相对量，参考温度对应无应力状态，装配应力、残余应力都会影响这一基准。第二，区分平均系数与瞬时系数，大温差工况优先使用积分形式或随温度变化的材料数据。第三，约束度决定应力上限，完全约束公式 σ = −EαΔT 是保守估算，实际结构应结合边界刚度折减。第四，多材料连接处要分别计算自由伸长量，两者之差才是被约束的变形量，这是异种材料接头热应力分析的关键。

常见误区包括：把热应变直接乘以弹性模量当作热应力（忽略了约束条件）；在温度跨度大时仍用常温 α；以及把平均系数当作瞬时系数代入积分。避免这些错误，热弹性计算的结果才具备工程可信度。

小结本节逻辑链：温度变化产生自由热应变 → 约束使自由热应变无法完全释放 → 剩余部分成为弹性应变 → 通过胡克定律转化为热应力。热膨胀系数的温度相关性通过平均系数或积分方式进入计算，约束程度则决定热应力的实际大小。掌握这条链条，就掌握了热弹性理论最基本的分析框架。