DIN 2093 疲劳裂纹扩展：Paris 公式与剩余寿命

1. 疲劳裂纹扩展的基本概念

碟形弹簧在交变载荷下，即使初始应力未超过材料的屈服强度，表面或内部的微小缺陷（如夹杂、腐蚀坑、加工刀痕）也可能在循环应力作用下逐渐扩展，最终达到临界尺寸发生失稳断裂。这一过程称为疲劳裂纹扩展，其速率由 Paris 公式 描述。

疲劳裂纹扩展分析用于： - 损伤容限设计：确定在定期检查间隔内，允许的缺陷大小。 - 剩余寿命预测：已知当前裂纹尺寸，估算还能承受多少循环。 - 材料与工艺选择：比较不同材料、表面处理对裂纹扩展抗力的影响。

2. Paris 公式

在稳定扩展阶段（裂纹扩展速率 $da/dN$ 与应力强度因子幅 $\Delta K$ 在双对数坐标下呈线性关系），Paris 公式为：

• $a$ — 裂纹深度或长度（mm 或 m）；
• $N$ — 应力循环次数；
• $da/dN$ — 每次循环裂纹扩展量（mm/cycle 或 m/cycle）；
• $\Delta K$ — 应力强度因子幅值（MPa·√m），$\Delta K = K_{max} - K_{min}$；
• $C$ — Paris 常数（与材料、应力比、环境有关）；
• $m$ — Paris 指数，对于大多数金属材料，$m \approx 3 \sim 4$。

该公式表明，裂纹扩展速率由 $\Delta K$ 主导，且对 $\Delta K$ 的变化十分敏感（因幂指数 $m$ 较大）。

3. 应力强度因子幅 $\Delta K$ 的计算

对于碟形弹簧的典型裂纹形态（表面半椭圆裂纹、边缘角裂纹），$\Delta K$ 可表达为：

• $\Delta \sigma$ — 危险点的应力幅（MPa），对于疲劳问题，取该点在工作循环中最大应力与最小应力的差值。例如，碟簧 I 点（下表面外缘）承受拉应力，则 $\Delta \sigma = \sigma_{I,max} - \sigma_{I,min}$。
• $Y$ — 几何修正因子（无量纲），取决于裂纹形状、构件尺寸和边界条件。常用表面半椭圆裂纹的 $Y \approx 0.65 \sim 1.2$，具体可查断裂力学手册。
• $a$ — 裂纹深度（mm 或 m）。

注意：若危险点处于压应力循环（如 OM 点始终受压），裂纹面会闭合，裂纹扩展驱动力 $\Delta K$ 极小，裂纹往往不扩展或扩展极慢。因此，疲劳裂纹扩展分析应重点关注拉应力区，如碟簧的 I 点、II 点等可能受拉的位置。

4. 材料常数 $C$ 与 $m$

对于高强度弹簧钢（如 50CrV4、C75S），在室温大气环境下，Paris 常数典型值：

注：$C$ 值强烈依赖于 $\Delta K$ 的单位（MPa·√m 或 N/mm³/²），使用时应确保单位统一。此处 $C$ 按 $\Delta K$ 以 MPa·√m 计，$a$ 以 mm 计。

应力比 $R$ 的影响：Paris 公式未显含平均应力，但 $R = \sigma_{min}/\sigma_{max}$ 会影响裂纹闭合程度，进而影响有效 $\Delta K_{eff}$。设计时可保守使用 $R = 0$ 的数据，或选用 Forman 修正公式。

5. 剩余寿命积分

从初始裂纹尺寸 $a_0$ 扩展到临界裂纹尺寸 $a_c$ 所需的循环数 $N_f$，由 Paris 公式积分得到：

分离变量并积分：

若假设几何因子 $Y$ 在裂纹扩展过程中为常数（或取平均值），则可直接积分。对于 $m \neq 2$，积分结果为：

当 $m = 2$ 时（极少），则有对数形式。多数钢材 $m \approx 3$，则积分式变为：

参数说明： - $a_0$ — 初始裂纹尺寸（mm），取决于无损检测能力或假设的初始缺陷大小（如 0.1 mm）； - $a_c$ — 临界裂纹尺寸（mm），由断裂韧性确定：$a_c = \frac{1}{\pi} \left( \frac{K_{IC}}{Y \cdot \sigma_{max}} \right)^2$； - $\Delta \sigma$ — 应力幅（MPa），取工作循环中危险点的应力变化范围。

6. 在碟形弹簧中的应用步骤

1. 确定危险点：通常为拉应力区，如 I 点（下外缘）或 II 点（上外缘），这些位置在循环中可能承受拉应力或拉‑压缩循环，裂纹易于扩展。
2. 应力分析：通过弹性或弹塑性 FEA，计算对应最小和最大压缩量下的应力，得到 $\sigma_{max}, \sigma_{min}$ 及 $\Delta \sigma$。
3. 确定初始裂纹：根据无损检测灵敏度或保守假定（如 $a_0 = 0.1\ \text{mm}$）。
4. 确定临界裂纹：由材料 $K_{IC}$ 和 $\sigma_{max}$ 计算 $a_c$。
5. 选择 $Y$ 因子：根据裂纹形态和位置查表或计算。例如碟簧内缘表面裂纹，可近似取 $Y \approx 1.12$（表面浅裂纹修正）。
6. 计算剩余寿命 $N_f$：代入上述积分式。
7. 安全系数：疲劳裂纹扩展寿命预测分散性大，建议取安全系数 $S_{crack} \ge 2.0$，即要求 $N_f \ge N_{service} \cdot S_{crack}$。

7. 计算示例

已知： - 碟簧 I 点：$\sigma_{max} = 800\ \text{MPa}$，$\sigma_{min} = 100\ \text{MPa}$ → $\Delta \sigma = 700\ \text{MPa}$，$R \approx 0.125$ - 材料 50CrV4，$K_{IC} = 65\ \text{MPa·√m}$，Paris 常数 $C = 3.0 \times 10^{-12}$（mm/cycle, MPa·√m 单位），$m = 3.2$ - 表面半椭圆裂纹，深度方向为 $a$，长度方向较长，简化取 $Y = 0.72$（经查表） - 初始裂纹 $a_0 = 0.1\ \text{mm}$，临界裂纹 $a_c$ 由 $K_{IC} = Y \cdot \sigma_{max} \cdot \sqrt{\pi a_c}$ 求得：

计算剩余寿命（$m=3.2 \neq 2$）：

分步计算： - $Y \cdot \Delta \sigma = 0.72 \times 700 = 504\ \text{MPa}$ - $(504)^{3.2} = \exp(3.2 \ln 504) = \exp(3.2 \times 6.2226) \approx \exp(19.912) \approx 4.42\times10^8$ - $\pi^{1.6} = \exp(1.6 \ln \pi) = \exp(1.6 \times 1.1447) \approx \exp(1.8315) \approx 6.24$ - 分母部分：$1.2 \times 3.0\times10^{-12} \times 4.42\times10^8 \times 6.24 \approx 1.2 \times 3.0\times10^{-12} \times 2.76\times10^9 \approx 1.2 \times 8.28\times10^{-3} = 9.94\times10^{-3}$ - 分子 $2 / (9.94\times10^{-3}) \approx 201.2$ - 括号项：$1/0.1^{0.6} - 1/4.05^{0.6}$
$0.1^{0.6} = \exp(0.6 \ln 0.1) = \exp(0.6 \times (-2.3026)) = \exp(-1.3816) \approx 0.251$

，倒数为
$4.05^{0.6} = \exp(0.6 \ln 4.05) = \exp(0.6 \times 1.3987) \approx \exp(0.8392) \approx 2.315$

，倒数为
差值 $= 3.98 - 0.432 = 3.548$

次循环。

解读：预期裂纹扩展寿命仅约 714 次，表明该应力水平过高，不适合长期使用。若要求 10^5 次寿命，必须大幅降低应力幅。

8. 设计建议与局限

• 压应力区无需裂纹扩展评估：OM 点压应力循环裂纹扩展极慢，通常不控制寿命；重点评估 I 点等拉应力区。
• Paris 公式适用范围：用于应力强度因子幅 $\Delta K$ 高于门槛值 $\Delta K_{th}$ 且低于快速断裂区的情况。对于弹簧钢，$\Delta K_{th} \approx 3 \sim 6\ \text{MPa·√m}$，低于此值裂纹不扩展。
• 变幅载荷：实际载荷谱若为变幅，需使用等效 $\Delta K$ 或采用条带法，不能用单一 $\Delta \sigma$。
• 安全系数：由于裂纹扩展数据分散性大，设计应采用较高的安全系数（≥ 2），并结合定期无损检测。
• 试验验证：关键件建议进行裂纹扩展速率试验，确定准确的 $C$ 和 $m$。

总结：Paris 公式是碟形弹簧损伤容限设计的核心工具，通过积分可预测从初始裂纹到临界裂纹的剩余寿命。正确选择危险点应力幅、几何因子和材料常数，并运用断裂力学准则，可有效保证碟簧在有限寿命内的安全运行。