应力-寿命（S-N）方法全解析
从Wöhler曲线到Miner损伤累积
疲劳分析用于预测结构在循环载荷下的使用寿命（载荷循环次数）。应力-寿命（S-N）方法是最经典的单轴疲劳分析方法，由Wöhler首次提出。本文系统梳理S-N方法的核心内容：S-N曲线、循环计数、等效名义应力、平均应力修正与损伤累积模型。
一、S-N曲线（S-N Curve）
S-N曲线由Wöhler首次提出，定义了应力与破坏循环次数之间的关系。材料的S-N曲线（以及其他疲劳性能）通常通过实验获得——即完全反转的旋转弯曲试验（fully reversed rotating bending tests）。
由于试验结果通常存在较大的分散性，还应提供数据的统计特征描述：根据曲线的标准误差，用存活概率（certainty of survival）对S-N曲线进行修正——可靠性要求越高，所需的存活概率越大。
图1：试验得到的S-N数据
当S-N试验数据以双对数坐标绘制——交变名义应力幅值Sa或应力范围SR对破坏循环数N——S与N的关系可以用直线段描述。通常采用一段式或两段式理想化。
图2：双对数坐标下的一段式S-N曲线
第1段的S-N关系为：
式中：S为名义应力范围；Nf为疲劳破坏循环数；b₁为第一疲劳强度指数；S₁为疲劳强度系数。
S-N方法基于弹性循环加载，因此S-N曲线在寿命轴上应限制在大于1000次循环的范围内，以确保不发生明显的塑性变形——这通常被称为高周疲劳（high-cycle fatigue）。
S-N曲线数据可查阅各大材料手册，比如MMPDS等。
二、循环计数（Cycle Counting）
循环计数用于从随机载荷序列中提取离散的"等效"常幅循环。可以将其理解为不断变化的应力-应变随时间信号：循环计数统计应力-应变滞后环的数量，并记录它们的范围/均值或最大/最小值。
雨流计数法（Rainflow）是应用最广泛的循环计数方法。它要求将应力时间历程重新排列，使其只包含峰谷值，并以绝对值最大的峰值或谷值作为起点。然后取三个连续应力点S₁、S₂、S₃，定义两个连续范围：
只有当ΔS₁≤ΔS₂时，才从S₁到S₂提取一个循环。一旦提取了一个循环，构成该循环的两个点即被丢弃，剩余的点重新连接；重复此过程，直到剩余数据点耗尽。
图3：雨流计数法确定循环
三、等效名义应力（Equivalent Nominal Stress）
由于S-N理论处理的是单轴应力，因此需要在每个计算点、每个时间步将应力分量合成为一个综合值，作为施加到S-N曲线上的等效名义应力。
提供了多种应力组合类型，默认类型为"绝对最大主应力"（Absolute maximum principle stress）。其中：脆性材料推荐使用"绝对最大主应力"，延性材料推荐使用"带符号的von Mises应力"（Signed von Mises stress）。带符号参数的符号取自最大绝对主应力值的符号。
四、平均应力修正（Mean Stress Correction）
通常，S-N曲线来自完全反转循环加载的标准试验。然而，实际疲劳载荷往往不是完全反转的，法向平均应力对构件疲劳性能有显著影响：拉伸平均应力是有害的，而压缩平均应力是有利的。平均应力修正用于考虑非零平均应力的影响。
在Haigh坐标系中，考虑平均应力影响时广泛使用Gerber抛物线和Goodman直线，可表示为：
Gerber：
Goodman：
式中：
Sm（平均应力）：
Sr（应力范围）：
Se：平均应力修正后的应力范围（对应应力范围Sr和平均应力Sm）
Su：极限强度
Gerber方法对正、负平均应力修正一视同仁——平均应力总会加速疲劳破坏；而Goodman方法忽略负的平均应力。两种方法对压缩平均应力均给出保守结果。其中：脆性材料推荐Goodman方法，延性材料推荐Gerber方法。对于Goodman方法，若拉伸平均应力大于UTS（极限抗拉强度），损伤将大于1.0；对于Gerber方法，无论拉伸还是压缩，若平均应力大于UTS，损伤也将大于1.0。
Haigh图描述了给定破坏循环数下应力幅值与平均应力的不同组合：
图4：Haigh图与平均应力修正方法
影响平均应力影响的参数可在Fatigue Module对话框和Assign Material对话框中定义。
1.GERBER2
改进了Gerber方法，通过忽略负平均应力的影响来改进修正。
2.SODERBERG
与GOODMAN略有不同：平均应力用屈服应力而非极限抗拉强度进行归一化：
式中：Se为等效应力幅值；Sa为应力幅值；Sm为平均应力；Sy为屈服应力。
3.FKM方法
若平均应力修正只指定一个斜率，则平均应力修正对应的平均应力敏感度值M等于第2区间的斜率。依据FKM指南，Haigh图按应力比的值划分为四个区间。修正后的值随后用于选择损伤与寿命计算阶段的S-N曲线。
图5：FKM应力修正对话框（1斜率）
注：下面的FKM方程说明了修正应力幅值（SeA）的计算方法。损伤计算中实际使用的应力是修正应力范围，即：
这些方程适用于用户输入的S-N曲线（默认情况下，任何用户自定义S-N曲线均按应力比R=1.0输入）。
HyperLife提供两种FKM修正选项，通过在Assign Material对话框中将FKM MSS设置为1slope（单斜率）或4slopes（四斜率）来激活。
（1）单斜率（1slope）情形
仅定义单一斜率，且SN模块的平均应力修正设置为FKM时：
区间1（R>1.0）：
区间2（−∞≤R≤0.0）：
区间3（0.0<R<0.5）：
区间4（R≥0.5）：
（2）四斜率（4slopes）情形
定义全部四个斜率，且平均应力修正设置为FKM时，各区间对应的平均应力敏感度值即为控制四个区间的斜率：
区间1（R>1.0）：
区间2（−∞≤R≤0.0）：
区间3（0.0<R<0.5）：
区间4（R≥0.5）：
式中：SeA为平均应力修正后的应力幅值（持久极限/Endurance stress）；Sm为平均应力；Sa为应力幅值；M为区间2输入的斜率（单斜率情形）；M₁~M₄为各区间的斜率（四斜率情形）；Se为完全反转疲劳强度（持久极限）。
五、损伤累积模型（Damage Accumulation Model）
采用Palmgren-Miner线性损伤累加法则。当满足下式时预测失效：
Nif：材料在应力幅值与平均应力组合i下，由其S-N曲线得到的疲劳寿命（破坏循环数）；
ni：载荷水平i下的应力循环数；
Di：ni次载荷循环下的累积损伤。
线性损伤累加法则没有考虑循环疲劳载荷下载荷顺序对损伤累积的影响；然而，它已被证明在许多应用中表现良好。
好了，这期的分享就到这里啦！如果内容对你有帮助，欢迎点赞、转发、收藏。你们的每一次鼓励，都是我持续输出的动力。感谢大家的支持，我们下期见！
【免责声明】本文内容仅供技术交流与学习参考，部分素材来源于公开文献或网络，版权归原作者所有。本文所述观点系基于笔者对资料的个人理解与归纳，不代表原作立场。若涉及版权问题或认为引用不当，请及时告知，我们将在核实后第一时间处理(删除或替换)。感谢您的理解与支持！
#疲劳分析#S-N曲线#HyperLife#CAE仿真#结构强度
