符号	名称	单位
stress_MPa	应力	MPa
temp_C	温度	°C
time_hours	时间	h

F-SSHT-K107 蠕变应变累积与寿命评估

1. 蠕变三阶段曲线

典型蠕变曲线由三个阶段组成，反映了材料在高温恒载下应变随时间的变化规律：

第一阶段（减速蠕变 / 初始蠕变）：
加载瞬间产生弹性应变，随后蠕变速率随时间逐渐降低。微观上，位错运动受到增殖位错和亚晶界阻碍，产生应变硬化效应，使材料“越来越难”蠕变。

第二阶段（稳态蠕变 / 最小蠕变速率）：
应变硬化与热回复（位错攀移湮灭）达到动态平衡，蠕变速率保持恒定且最小。这是碟簧服役寿命中最长的阶段，也是 Norton-Bailey 公式直接描述的对象。

第三阶段（加速蠕变 / 最终断裂）：
微观孔洞聚合成微裂纹，有效承载面积减小，真实应力急剧升高，蠕变加速直至断裂。一旦进入第三阶段，碟簧寿命已接近终点，必须立即更换。

2. 各阶段本构方程

2.1 第一阶段：时间硬化模型

\boxed{\varepsilon_I(t) = \varepsilon_0 + A_1 \cdot \sigma^n \cdot t^m}

$\varepsilon_0$ ：初始弹性应变 ( $\sigma/E$ )
$$A_1, n, m$$ ：材料常数（ $$0 < m < 1$$ ，例如 $m \approx 0.3 \sim 0.5$ ）
$\sigma$ ：恒定应力 (MPa)
$$t$$ ：时间 (h)

蠕变速率随时间递减： $\dot{\varepsilon}_I = A_1 \cdot \sigma^n \cdot m \cdot t^{m-1}$ ，因 $$(m-1) < 0$$ ，速率下降。

2.2 第二阶段：稳态蠕变（Norton-Bailey）

\boxed{\dot{\varepsilon}_{II} = A_2 \cdot \sigma^n \cdot \exp\left( -\frac{Q}{R \cdot T} \right)}

\boxed{\varepsilon_{II}(t) = \dot{\varepsilon}_{II} \cdot t}

这是蠕变寿命计算的核心。 $\dot{\varepsilon}_{II}$ 越小，碟簧在高温下维持预紧力的能力越强。

2.3 第三阶段：损伤耦合加速蠕变

第三阶段无法用简单幂律描述，通常引入连续损伤力学，将损伤变量 $$D$$ 与蠕变速率耦合：

\dot{\varepsilon}_{III} = \frac{\dot{\varepsilon}_{II}}{(1 - D)^\phi}

\dot{D} = \frac{A_3 \cdot \sigma^\chi}{(1 - D)^\psi}

$$D$$ ：损伤变量， $$D=0$$ 表示完好， $$D=1$$ 表示断裂。
$\phi, \psi, \chi$ ：材料损伤指数，通常通过试验拟合。

当 $D \to 1$ 时，蠕变速率趋于无穷，对应宏观裂纹形成。该阶段在工程中只用于剩余寿命预警，不作为正常设计状态。

3. 蠕变应变阈值与更换准则

对于碟形弹簧，蠕变应变直接影响其几何尺寸和力学性能：

蠕变应变 ≤ 0.5%：自由高度损失极小，刚度变化可忽略，状态安全。
0.5% < 蠕变应变 ≤ 1.0%：需加强监测，评估预紧力是否仍满足功能要求。若预紧力降额已接近下限，应准备备件。
蠕变应变 > 1.0%：自由锥高 $$h_0$$ 显著减小，展平力下降可能超过 10%，且材料已积累大量微观损伤，建议立即更换。

工程阈值： $\varepsilon_c > 1\%$ 是碟簧因蠕变丧失弹性补偿能力的临界点，必须作为维护计划的强制更换指标。

4. 累积蠕变应变计算

当碟簧在变温变载工况下工作时，各阶段的蠕变应变需按时间分段累加，并考虑应力松弛的耦合效应。

分段累加公式（适用于温度和应力分段变化）：

\boxed{\varepsilon_{c,total} = \sum_{i=1}^{N} \dot{\varepsilon}_{II}(T_i, \sigma_i) \cdot \Delta t_i}

$(T_i, \sigma_i)$ ：第 $$i$$ 段时间区间内的温度和 OM 点应力。
$\Delta t_i$ ：第 $$i$$ 段的持续时间 (h)。
$$N$$ ：时间分段总数。

应力松弛耦合修正：
碟形弹簧在恒定位移下工作，蠕变导致应力不断下降（松弛）。此时需将 Norton-Bailey 公式代入松弛方程迭代求解：

\frac{d\sigma}{dt} = -E \cdot \dot{\varepsilon}_c(\sigma, T)

逐步数值积分，同时计算蠕变应变累积和应力衰减，直至达到设计寿命或蠕变应变超限。

5. 计算示例

已知：H13 碟簧，工作温度 500 °C (773 K)，OM 点初始应力 $\sigma = 600\ \text{MPa}$ 。
Norton 参数： $A_2 = 8.0 \times 10^{-20}$ ， $$n = 5.5$$ ， $Q = 320,000\ \text{J/mol}$ 。
设计要求：运行 20,000 小时，蠕变应变 ≤ 1%。

步骤1：计算稳态蠕变速率

\dot{\varepsilon}_{II} = 8.0 \times 10^{-20} \cdot (600)^{5.5} \cdot \exp\left( -\frac{320000}{8.314 \times 773} \right)

$600^{5.5} = \exp(5.5 \cdot \ln 600) \approx \exp(35.2) \approx 1.9 \times 10^{15}$
$\frac{Q}{RT} = \frac{320000}{8.314 \times 773} \approx 49.7$
$\exp(-49.7) \approx 2.5 \times 10^{-22}$ $\dot{\varepsilon}_{II} \approx 8.0 \times 10^{-20} \times 1.9 \times 10^{15} \times 2.5 \times 10^{-22} \approx 3.8 \times 10^{-26}\ \text{s}^{-1}$

换算为每年（ $3.15 \times 10^7$  s）蠕变量： $\approx 1.2 \times 10^{-18}$ ，20,000 小时 ≈ 2.28 年，累积应变 完全可忽略。

步骤2：修正应力效应
若应力升高至 900 MPa（因碟簧过载），则：

900^{5.5} = \exp(5.5 \cdot \ln 900) \approx \exp(37.4) \approx 1.7 \times 10^{16}

速率增大 约 9 倍，但仍极低。
但若同时温度升至 600 °C (873 K)， $\frac{Q}{RT} \approx 44.1$ ， $\exp(-44.1) \approx 7.0 \times 10^{-20}$ ，速率将暴增 6 个数量级以上，20,000 小时内的蠕变应变可能超过 5%，必须更换材料为 Inconel 718。

总结：蠕变应变累积遵循三阶段规律。稳态蠕变阶段决定碟簧长时寿命，Norton-Bailey 公式是量化核心。工程中设定 1% 蠕变应变为强制更换阈值，并结合变温变载分段累积法进行全寿命评估。

蠕变应变累积

公式表达式

参数列表

需要计算该公式？

详细计算指南