4. 热循环下的弹性与塑性行为

热循环环境下，材料承受的应力并非来自外加载荷，而是由温度变化与约束条件共同“逼”出来的。理解热应力起源、约束分类、屈服判据、包辛格效应以及循环塑性本构模型的适用边界，是判断热循环退化路径的基础。本节从物理机制出发，逐步过渡到工程分析中可直接使用的判据与模型选择方法。

热应力的起源：自由应变被约束“截留”

温度变化使材料产生自由热应变 ε_th = α·ΔT（α 为线膨胀系数，ΔT 为温差）。若材料完全自由，该应变不产生应力；一旦变形受到约束，自由热应变与实际发生的总应变之间出现差值，弹性部分被“截留”为应力：σ = E·(ε_total − ε_th − ε_p)。这就是热应力的本质——约束把热膨胀/收缩转化为弹性应变。

以某奥氏体不锈钢管件为例（示例参数，非实测）：α = 17×10⁻⁶ /℃，E = 180 GPa，ΔT = 200℃。若两端完全固定，自由热应变约 3.4×10⁻³，弹性模量对应热应力约 612 MPa，远超其室温屈服强度约 250 MPa。这说明在强约束下，热循环几乎必然进入塑性。实际工程中，约束往往介于完全自由与完全固定之间，因此需要先对约束条件分类。

约束条件分类：决定热应力上限

工程中常按约束来源与刚度分为三类：

• 外部机械约束：构件被支座、法兰、焊缝或相邻结构限制，位移或转角被强制约束，热应变几乎全部转化为应力，最容易进入塑性。
• 内部自约束：同一构件内温度梯度导致不同区域膨胀量不同，如厚壁管内外壁、焊接热影响区与母材之间，产生自平衡应力。
• 材料不匹配约束：异种材料连接（如不锈钢—碳钢）因 α 与 E 差异，在界面附近形成局部高应力。

约束刚度越高、ΔT 越大、α 越大，热应力越大。分析时应先判断约束类型，再决定是否需要考虑塑性。若弹性计算应力已超过屈服强度，就必须引入塑性行为，否则会严重高估应力、低估实际变形。

屈服判据：从单轴到多轴

单轴条件下，当 |σ| ≥ σ_y 时材料屈服。热循环构件多处于多轴应力状态，需采用等效应力判据。工程分析中常用 von Mises 等效应力：

σ_eq = sqrt(0.5·[(σ1−σ2)² + (σ2−σ3)² + (σ3−σ1)²])
屈服条件：σ_eq ≥ σ_y(T)

注意 σ_y 随温度升高而下降。例如某低合金钢室温 σ_y = 350 MPa，400℃ 时降至约 220 MPa（示例值）。热循环中温度与应力同步变化，必须使用对应温度下的屈服强度。若采用 Tresca 判据，则以最大剪应力为准，结果偏保守。对于循环加载，还需区分首次屈服与循环稳定后的屈服，后者受包辛格效应影响。

包辛格效应：反向屈服提前

材料在某一方向塑性变形后，反向加载时的屈服强度会降低，这就是包辛格效应。例如拉伸至塑性应变 0.5% 后卸载再反向压缩，压缩屈服可能从 −350 MPa 降至 −280 MPa（示例）。在热循环中，升温压缩、降温拉伸交替出现，包辛格效应使反向屈服更早发生，循环塑性应变幅增大，从而加速疲劳与棘轮效应。

若本构模型采用各向同性硬化，将无法描述这一现象，会低估反向塑性。因此涉及明显反向加载的热循环分析，应优先选择随动硬化或混合硬化模型。

循环塑性本构模型及其适用边界

热循环分析中常用的循环塑性本构模型主要有三类，其适用边界差异明显：

选择模型时，先判断是否存在反向塑性。若热循环中应力幅超过屈服强度且反向加载明显，各向同性硬化会低估塑性应变，导致疲劳寿命高估。随动硬化中的非线性模型（如 Chaboche）可同时描述包辛格效应与循环稳定，适合多数热循环退化分析。混合硬化则用于同时存在循环硬化和棘轮的场景。

在热循环退化分析中的落地要点

实际分析中，温度场与力学场往往解耦或弱耦合。以管道热循环为例，先做瞬态热分析得到温度历史，再映射到结构分析中。关键输出包括：等效塑性应变幅 Δε_p/2、棘轮应变率、以及每个循环的应力范围。若 Δε_p/2 稳定，说明进入弹性安定；若持续累积，则发生棘轮，需调整约束或材料。包辛格效应会放大反向塑性，因此本构参数中的随动硬化模量 C 和背应力饱和值需从循环试验提取，不能直接沿用单调拉伸数据。

最后，热循环退化分析不是一次计算，而是“判据—模型—参数—验证”的闭环。先用屈服判据判断是否进入塑性，再根据反向加载程度选择本构，最后用循环稳定后的应力应变响应验证模型是否合理。只有本构适用边界与实际问题匹配，退化预测才具有工程可信度。