3、材料疲劳性能曲线：S-N曲线、P-S-N曲线、e-N曲线的绘制方法、数学表达及工程应用

材料疲劳性能曲线是疲劳设计与寿命评估的核心基础，它将外加循环应力（或应变）与材料达到失效时的循环次数（寿命）之间的定量关系以图形和数学形式表达出来。工程中常用的三类曲线——S-N曲线（应力-寿命）、P-S-N曲线（概率-应力-寿命）和e-N曲线（应变-寿命）——分别适用于不同的设计场景：S-N曲线用于高周疲劳（寿命>10⁴~10⁵次）的无限寿命设计；P-S-N曲线用于考虑分散性的可靠度设计；e-N曲线则用于低周疲劳（寿命<10⁴~10⁵次）的有限寿命设计，尤其是存在塑性应变的场合。本章将系统讲解这三类曲线的数据获取、绘制方法、数学拟合公式及其在工程实践中的典型应用。

疲劳性能曲线体系与工程应用映射 S-N曲线 应力-寿命（高周疲劳） N > 10⁴~10⁵ P-S-N曲线 概率-应力-寿命（可靠度） 含存活率p（如p=95%/99%） e-N曲线 应变-寿命（低周疲劳） N < 10⁴~10⁵，含塑性应变 数据获取：成组法（定应力水平，多根试样） + 升降法（测定条件疲劳极限） • 成组法：每个应力水平至少5~8根试样，取对数寿命均值，用于S-N曲线中段 • 升降法：相邻应力级差Δσ≈5%σ₀，配对数据点用于确定疲劳极限（10⁷次） • 应变控制：引伸计标距段应变幅恒定，记录半寿命稳定滞后环，用于e-N曲线 数学表达（常用形式） ① 幂函数式：σᵃ·N = C ② 指数式：σ = σ₀·e^(−kN) ③ 三参数式：σ = σₑ(1 + A/N^α) ④ Coffin-Manson：εₐ = (σ′f/E)(2N)^b + ε′f(2N)^c 工程应用场景 • 无限寿命设计：σ < σₑ（疲劳极限） • 有限寿命设计：Miner线性累积损伤 • 可靠度设计：P-S-N给出给定存活率下的许用应力 • 缺口/尺寸修正：修正系数Kσ、ε、β

上图展示了三类疲劳性能曲线在数据获取、数学表达和工程应用上的整体框架。S-N曲线是基础，P-S-N曲线是在S-N基础上引入概率统计处理，而e-N曲线则从应力控制转向应变控制，适用于存在明显塑性变形的低周疲劳问题。下面分别展开论述。

3.1 S-N曲线：高周疲劳设计的基础

3.1.1 核心原理与绘制方法

S-N曲线（Wöhler曲线）描述的是在给定应力比R（最小应力与最大应力之比）或平均应力σm下，施加的循环应力幅σa（或最大应力σmax）与试样达到失效的循环次数N之间的关系。其典型特征是在双对数坐标（lgσa - lgN）下，中寿命区（10³~10⁶次）呈现近似线性下降关系，而在长寿命区（N>10⁶~10⁷次）出现水平渐近线，对应的应力值即为疲劳极限σₑ（对于钢材等体心立方金属）。

绘制S-N曲线的标准流程如下：

  1. 试样制备与试验方案设计：采用光滑圆棒试样（直径6~10mm），表面磨削至Ra≤0.2μm。至少选取5~7个应力水平，最高应力水平应使寿命在10⁴~10⁵次，最低应力水平应接近预期疲劳极限。每个应力水平采用成组法测试5~8根试样。
  2. 试验执行：在旋转弯曲疲劳试验机或电磁谐振试验机上施加恒定幅值循环载荷，频率通常为20~200Hz，记录每根试样的失效循环次数Nf。若试样在10⁷次未失效，则记为“越出”，作为疲劳极限的初步估计。
  3. 数据整理：对每个应力水平下的寿命取对数，计算对数寿命均值lgN̄和标准差。剔除异常数据（如断口显示有明显夹杂物缺陷的试样）。
  4. 曲线拟合：在双对数坐标下，对中寿命区数据采用最小二乘法线性拟合，得到S-N曲线的斜线段；结合升降法测得的疲劳极限，形成完整的S-N曲线。
经验提醒：S-N曲线通常以应力幅σa或最大应力σmax为纵坐标，但必须标注应力比R或平均应力σm。同一材料在不同R下的S-N曲线差异显著，工程中常采用Goodman或Gerber公式进行平均应力修正，将非对称循环等效为对称循环（R=-1）下的应力幅。

3.1.2 数学表达

工程中最常用的S-N曲线数学表达式为幂函数式：

σᵃ · N = C

取对数后变为线性形式：

lgN = lgC - a·lgσ

其中a为斜率参数（通常为3~10，取决于材料），C为材料常数。对于中长寿命区，也可采用三参数式：

σ = σₑ · (1 + A / N^α)

该式能同时描述斜线段和水平渐近线，便于计算机程序直接调用。典型钢材的S-N曲线参数示例：45钢调质态，σb=750MPa，R=-1时，a≈8.5，C≈10³⁵·⁵（单位：MPa^a·cycle）。

3.1.3 工程应用与案例拆解

案例：汽车传动半轴疲劳寿命校核

某商用车传动半轴材料为40Cr，调质处理，σb=980MPa，σₑ=420MPa（R=-1）。设计寿命要求为30万公里，对应循环次数约2×10⁶次。实测半轴在最大扭矩工况下危险截面应力幅σa=380MPa，平均应力σm=120MPa。

分析过程：

  1. 采用Goodman公式将非对称循环等效为对称循环：σa_eq = σa / (1 - σm/σb) = 380 / (1 - 120/980) ≈ 433MPa。
  2. 查该材料S-N曲线（R=-1），当σa=433MPa时，对应寿命N≈8×10⁵次，小于要求的2×10⁶次。
  3. 结论：该半轴疲劳寿命不足，需改进设计。

改进措施与结果：将半轴花键根部圆角从R1.5mm增大至R3mm，表面滚压强化处理（表面残余压应力-400MPa），使σa降至320MPa。重新校核：σa_eq=320/(1-120/980)≈365MPa，对应寿命N≈5×10⁶次，满足要求。台架试验验证寿命达到2.3×10⁶次，与S-N曲线预测吻合良好。

风险提示：S-N曲线基于光滑试样获得，实际零件存在缺口、尺寸效应和表面状态差异，必须引入疲劳缺口系数Kf、尺寸系数ε和表面系数β进行修正。修正后的许用应力为：σ许用 = σₑ·ε·β / Kf。忽略修正将导致严重过估寿命。

3.2 P-S-N曲线：考虑分散性的可靠度设计

3.2.1 核心原理与绘制方法

疲劳寿命数据具有显著的分散性——即使在相同应力水平下，同一批试样的寿命可相差3~5倍甚至更多。P-S-N曲线以存活率p（可靠度）为参数，给出在给定存活率下应力与寿命的关系。例如p=50%对应中值S-N曲线，p=95%表示有95%的试样寿命大于该曲线给出的N值。

绘制P-S-N曲线的关键步骤：

  1. 多应力水平成组试验：每个应力水平至少需要10~15根试样，以获取可靠的寿命分布统计量。
  2. 寿命分布检验：对每个应力水平下的寿命取对数，检验其是否服从正态分布（通常lgN服从正态分布或Weibull分布）。采用Shapiro-Wilk检验或K-S检验。
  3. 计算各应力水平下的对数寿命均值与标准差：lgN̄p = lgN̄ + up·s，其中up为标准正态偏量（p=50%时up=0；p=95%时up=1.645；p=99%时up=2.326）。
  4. 拟合各存活率下的S-N曲线：将不同应力水平下相同p值的lgNp点连接，用最小二乘法拟合直线，得到一组平行或近似平行的P-S-N曲线族。

3.2.2 数学表达与工程应用

P-S-N曲线的数学表达通常为：

lgNp = ap + bp·lgσ

其中ap和bp为与存活率p相关的材料常数。工程中更常用的是给定存活率下的疲劳极限σₑp，例如σₑ95表示95%存活率下的疲劳极限。

存活率p标准正态偏量up疲劳极限σₑp (MPa)应用场景
50%0420一般机械设计
90%1.282398汽车零部件
95%1.645388航空结构件
99%2.326372关键安全部件

案例：风电叶片螺栓连接疲劳可靠度设计

某2MW风电叶片根部T型螺栓连接，材料为42CrMoA，σb=1080MPa。螺栓承受拉-拉疲劳载荷，应力幅σa=180MPa，平均应力σm=300MPa。设计要求20年寿命（约10⁷次循环），可靠度不低于99%。

分析过程：

  1. 查该材料P-S-N曲线（R=0.1），p=99%时，10⁷次对应的许用应力幅σa99=165MPa。
  2. 考虑螺纹缺口系数Kf=3.2，尺寸系数ε=0.8，表面系数β=0.9，则实际许用应力幅：σ许用 = 165×0.8×0.9/3.2 ≈ 37MPa。
  3. 实际应力幅180MPa远大于37MPa，疲劳强度严重不足。

改进方案：采用预紧力控制使螺栓始终处于压应力状态（σm=-200MPa），此时应力幅降至σa=80MPa，且压应力状态显著提高疲劳强度。重新校核：查P-S-N曲线（R=-2），p=99%时10⁷次许用应力幅为95MPa，考虑修正后σ许用=95×0.8×0.9/3.2≈21.4MPa，仍不足。最终方案改为增加螺栓数量（从24根增至36根），使单根螺栓应力幅降至50MPa，同时采用冷滚压螺纹工艺（表面系数提高至1.2），σ许用=95×0.8×1.2/3.2≈28.5MPa，仍不足。最后通过增大螺栓规格（M36→M42）并优化预紧力，使应力幅降至25MPa，满足要求。

复盘要点：P-S-N曲线揭示了疲劳强度的概率本质。在关键安全部件设计中，必须明确给出目标可靠度，并据此选择相应的P-S-N曲线。盲目采用中值S-N曲线（p=50%）进行设计，实际可靠度仅为50%，这在航空航天、风电、压力容器等领域是不可接受的。

3.3 e-N曲线：低周疲劳与应变寿命法

3.3.1 核心原理与绘制方法

当循环应力超过材料屈服强度时，塑性应变成为主导因素，此时应力-寿命关系（S-N曲线）不再适用，必须采用应变-寿命（e-N）曲线。e-N曲线描述总应变幅εa（弹性应变幅εea与塑性应变幅εpa之和）与失效寿命2Nf（反向次数）之间的关系。

e-N曲线的核心是Coffin-Manson方程：

εa = (σ′f / E) · (2Nf)^b + ε′f · (2Nf)^c

其中：σ′f为疲劳强度系数（约等于真实断裂强度σf），b为疲劳强度指数（-0.05~-0.12），ε′f为疲劳延性系数（约等于真实断裂延性εf），c为疲劳延性指数（-0.5~-0.7），E为弹性模量。第一项为弹性应变幅贡献，第二项为塑性应变幅贡献。

绘制e-N曲线的试验方法：

  1. 应变控制疲劳试验：采用轴向应变控制，应变比Rε=-1（完全反向），应变幅从0.2%到2%选取5~7个水平。
  2. 滞后环记录：在每个应变水平下，记录半寿命时的稳定滞后环，分解出弹性应变幅εea和塑性应变幅εpa。
  3. 双对数线性拟合：分别对εea-2Nf和εpa-2Nf数据进行双对数线性回归，得到Coffin-Manson方程中的四个材料常数。

3.3.2 工程应用与案例拆解

案例：压力容器接管根部热疲劳寿命评估

某化工压力容器接管根部，材料为16MnR（Q345R），E=206GPa，σ′f=930MPa，b=-0.09，ε′f=0.26，c=-0.56。设备启停过程中，接管根部因温差产生应变幅εa=0.35%（其中弹性应变幅0.15%，塑性应变幅0.20%）。每天启停一次，设计寿命15年（约5500次）。

计算过程：

εa = (σ′f/E)(2Nf)^b + ε′f(2Nf)^c
0.0035 = (930/206000)(2Nf)^(-0.09) + 0.26(2Nf)^(-0.56)

采用数值迭代求解：设2Nf=10000，则弹性项=(930/206000)×10000^(-0.09)=0.00451×0.437=0.00197；塑性项=0.26×10000^(-0.56)=0.26×0.00575=0.00150；总应变幅=0.00347≈0.0035。因此2Nf≈10000次，即Nf≈5000次。

结论与措施：计算寿命5000次略低于设计要求的5500次，存在风险。改进措施：在接管根部增加过渡圆角（R10→R25），降低应力集中系数，使应变幅降至0.28%。重新计算：2Nf≈25000次，Nf≈12500次，满足要求。同时建议将启停速率从10℃/h降至5℃/h，进一步降低热应力。

常见误区：e-N曲线中的寿命2Nf是反向次数（一个循环包含两次反向），而工程中通常关心循环次数Nf。两者相差2倍，在计算中务必注意。另外，e-N曲线对平均应变敏感，当存在非零平均应变时，需采用Morrow修正或Smith-Watson-Topper（SWT）参数进行修正。

3.4 三类曲线的对比与工程选型

对比维度S-N曲线P-S-N曲线e-N曲线
控制变量应力幅σa应力幅σa（含概率）应变幅εa
适用寿命范围N>10⁴~10⁵（高周）N>10⁴~10⁵（高周）N<10⁴~10⁵（低周）
是否考虑塑性否（弹性为主）否（弹性为主）是（弹塑性）
分散性处理中值曲线存活率p参数化中值曲线（可扩展为P-e-N）
数据获取成本中（5~8根/级）高（10~15根/级）高（应变控制设备昂贵）
典型应用轴类、齿轮、弹簧航空结构、风电、压力容器热疲劳、焊接接头、缺口件

工程选型决策流程：首先判断设计寿命所处的循环次数范围——若N>10⁵次且应力水平低于屈服强度，选用S-N曲线；若涉及可靠度要求，升级为P-S-N曲线；若N<10⁴次或局部存在塑性应变（如缺口根部），必须采用e-N曲线。实际工程中，许多关键部件（如涡轮叶片、焊接结构）同时存在高周和低周疲劳损伤，需要采用Miner线性累积损伤法则结合两类曲线进行双轨评估。

实战建议：在缺乏试验数据时，可采用通用斜率法（Universal Slopes Method）估算e-N曲线参数：b≈-0.12，c≈-0.6，σ′f≈3.5σb，ε′f≈εf。对于钢材，εf≈ln(1/(1-ZA))，其中ZA为断面收缩率。该方法在初步设计阶段具有足够的工程精度。

综上所述，S-N曲线、P-S-N曲线和e-N曲线构成了材料疲劳性能表征的完整体系。工程技术人员应根据部件的载荷特征、寿命要求和可靠度目标，合理选择曲线类型，并严格遵循标准试验方法获取数据，正确应用数学表达式进行寿命预测。在实际工程中，三类曲线往往需要配合使用——例如，对同一个焊接结构，焊缝处采用e-N曲线评估低周疲劳，母材区采用S-N曲线评估高周疲劳，再通过累积损伤法则综合评估整体寿命。掌握这三类曲线的内在联系和工程适用条件，是进行准确疲劳寿命预测的前提。

← 上一章 📖 返回目录 下一章 →