2. 疲劳寿命预测的基本模型
疲劳寿命预测不是单一公式，而是一组按损伤机制和寿命区间分工的模型体系。工程上最常用的三条主线是：S-N 曲线（应力—寿命，适用于高周疲劳）、ε-N 曲线（应变—寿命，适用于低周疲劳与局部塑性显著的情形）、裂纹扩展公式（如 Paris 律，适用于已存在裂纹或缺陷的剩余寿命评估）。理解它们的适用边界，比记住公式本身更重要——因为不确定性评估的第一步，就是判断"该用哪个模型"，以及"这个模型的哪些输入本身就不确定"。
一、S-N 曲线：高周疲劳的基线模型
S-N 曲线描述恒幅应力幅 Sa 与失效循环数 Nf 之间的关系。在双对数坐标下，中高寿命段通常近似为直线：
σ_a = σ_f' · (2N_f)^b （Basquin 形式）
lg N_f = A - m · lg σ_a （工程幂律形式）其中 σf' 为疲劳强度系数，b 为疲劳强度指数（钢材常在 −0.05 ~ −0.12），m 为 S-N 曲线斜率。对钢材，曲线在约 106~107 次处常出现水平段，对应疲劳极限；铝合金等则无明确极限，需指定条件疲劳强度（如 107 次对应的应力）。
实际构件承受的是变幅载荷，需借助Miner 线性累积损伤把各级载荷折算为等效损伤：
D = Σ (n_i / N_i)
D ≥ 1 判为失效（工程上常取 0.3~0.7 作为考虑不确定性的临界值）在不确定性评估中，S-N 曲线的角色是提供"应力→寿命"的映射骨架。它的不确定性来源非常具体：材料批次导致的 σf'、b 分散；表面粗糙度、缺口系数 Kt、尺寸效应带来的修正系数偏差；以及 Miner 准则本身忽略载荷顺序的模型误差。因此工程上常把 S-N 曲线写成 P-S-N 曲线族（存活率 P = 50%、90%、99%），用对数寿命的标准差 σlgN 表征分散性。
二、ε-N 曲线：低周疲劳与局部应变法
当名义应力接近或超过屈服强度、循环寿命低于约 104~105 次时，塑性应变主导损伤，必须用应变寿命模型。最经典的是 Manson-Coffin 与 Basquin 合成的总应变方程：
Δε/2 = Δε_e/2 + Δε_p/2
= (σ_f'/E)(2N_f)^b + ε_f'(2N_f)^c其中 εf' 为疲劳延性系数，c 为疲劳延性指数（常在 −0.5 ~ −0.7），E 为弹性模量。第一项是弹性应变贡献，第二项是塑性应变贡献；两条直线在双对数坐标下的交点对应过渡寿命 Nt，是判断"该用 S-N 还是 ε-N"的量化依据。
ε-N 方法通常与Neuber 缺口修正配合使用，把名义弹性应力换算为缺口根部局部弹塑性应变：
Δσ · Δε = (K_t · ΔS)² / E在不确定性评估中，ε-N 曲线的角色是把不确定性从"载荷"传递到"局部应变"。循环应力—应变曲线（Ramberg-Osgood 参数 K'、n'）的分散、Neuber 修正的近似误差、缺口根部残余应力，都会显著放大寿命预测的方差。低周疲劳的寿命分散通常比高周疲劳更大，对数寿命标准差可达 0.2~0.4。
三、裂纹扩展公式：损伤容限与剩余寿命
当构件已存在可检裂纹或初始缺陷时，寿命由裂纹从初始尺寸 a0 扩展到临界尺寸 ac 的循环数决定。Paris 律给出稳定扩展阶段的关系：
da/dN = C · (ΔK)^m
ΔK = Y · Δσ · √(πa)其中 ΔK 为应力强度因子幅，Y 为几何修正因子，C、m 为材料常数（钢材 m 常在 2.5~4）。积分可得剩余寿命：
N_p = ∫[a_0→a_c] da / [C·(Y·Δσ·√(πa))^m]完整描述还需考虑门槛值 ΔKth（低于则不扩展）和断裂韧性 KIC（决定 ac）。在不确定性评估中，裂纹扩展公式的角色最为特殊：它把不确定性集中在初始裂纹尺寸 a0 和材料常数 C、m 上。由于 a0 通常只能通过无损检测给出概率分布（检出概率 POD 曲线），而 m 的微小偏差会经指数放大，裂纹扩展寿命的对数分散往往比 S-N 方法更大。
四、三类模型在不确定性评估中的角色对比
把三者放在同一框架下看，它们的差别不只是公式，而是不确定性的入口位置和放大机制不同。下表给出工程上常用的对照关系（数值为常见范围，属示例性参考，实际应以材料数据手册或试验为准）：
| 对比项 | S-N 曲线 | ε-N 曲线 | 裂纹扩展公式 |
|---|---|---|---|
| 适用寿命区间 | 高周，N > 10⁴~10⁵ | 低周，N < 10⁴~10⁵ | 含裂纹阶段，任意 N |
| 控制参量 | 应力幅 σ_a | 总应变幅 Δε/2 | 应力强度因子幅 ΔK |
| 主要材料常数 | σ_f'、b、疲劳极限 | σ_f'、b、ε_f'、c、K'、n' | C、m、ΔK_th、K_IC |
| 主导不确定源 | 材料分散、表面/缺口修正 | 循环本构分散、Neuber 近似 | 初始裂纹 a₀、C、m |
| 对数寿命标准差（示例） | 0.1~0.2 | 0.2~0.4 | 0.3~0.6 或更大 |
| 典型评估输出 | P-S-N 曲线、可靠度 R(N) | 局部应变—寿命分布 | 剩余寿命分布、检修间隔 |
五、实际评估中的模型选择与参数处理
以一个示例场景说明：某轴类零件材料为 45 钢，工作应力幅约 180 MPa，循环数约 5×105，表面精车，无宏观裂纹。判断流程如下：
- 先判区间：5×105 次属高周疲劳，优先用 S-N 曲线；若局部缺口处应力集中系数 Kt = 2.5，缺口根部名义应力已达 450 MPa，接近屈服，则应改用 ε-N 局部应变法复核。
- 再判缺陷：若超声检测发现 0.5 mm 当量裂纹，则无论寿命区间如何，都应以裂纹扩展公式计算剩余寿命，S-N 结果仅作参考。
- 参数取值：S-N 曲线取 P=50% 与 P=99% 两条，分别对应中值寿命与保守寿命；σlgN 取 0.15（示例假设），用于构造对数正态分布。
- 传递计算：对 σf'、b、Kt 各设定分布（如正态或对数正态），用蒙特卡洛抽样 104 次，得到 N_f 的分布直方图，进而读出 P_f 与 R(N)。
- 结果判读：若 P=99% 对应的寿命低于设计循环数，则需调整材料、表面处理或降低应力幅；若裂纹扩展寿命的 95% 置信下限低于检修间隔，则必须缩短检查周期。
综上，S-N、ε-N 与裂纹扩展公式分别覆盖高周、低周与含裂纹三类场景，它们在不确定性评估中的共同角色是提供寿命对输入参数的响应关系；差别在于不确定源清单、放大机制与分散程度。实际评估时，先按寿命区间和缺陷状态选模型，再列出该模型的关键随机输入，最后用蒙特卡洛或一次二阶矩方法把输入分布传递为寿命分布，才能得到可用于设计决策的可靠度指标。