第三章 疲劳理论基础
疲劳，这个词在工程圈里听得太多了。说白了，就是材料在反复加载下，慢慢「累坏了」的过程。我刚开始做电机轴系仿真那会儿，总觉得静强度够了就行，结果被一个扭振疲劳案例狠狠教育了一回——轴没断，但裂纹已经悄悄长出来了。从那以后，我对疲劳分析再也不敢马虎。
这一章，咱们把疲劳的底牌翻出来看看。我会结合自己做过的项目，把那些书本上没写的坑，一个一个指给你。
3.1 疲劳破坏机理：裂纹是怎么来的？
疲劳破坏，不是一下子发生的。它分三个阶段：裂纹萌生、裂纹扩展、瞬间断裂。
裂纹萌生，通常发生在应力集中区。比如轴肩、键槽、油孔这些地方。微观上，是材料内部的滑移带反复运动，慢慢挤出了微小的裂纹。我见过一个案例，电机轴在键槽根部出现了肉眼都看不见的微裂纹，但用显微镜一看，滑移带已经很明显了。
裂纹扩展，是裂纹慢慢长大的过程。每加载一次，裂纹尖端就往前推进一点点。这个阶段占了疲劳寿命的大部分时间。你想想看，如果能在扩展阶段早期发现裂纹，很多事故是可以避免的。
瞬间断裂，是最后的结局。当裂纹长到临界尺寸，剩下的截面撑不住了，「啪」一下就断了。断口上通常能看到两个区域：一个是光滑的疲劳区，一个是粗糙的瞬断区。我每次做失效分析，第一件事就是看断口形貌——它能告诉我很多信息。
核心要点：疲劳破坏的三个阶段——萌生、扩展、断裂。其中扩展阶段占寿命的80%以上。
注意：高周疲劳和低周疲劳的破坏机理不同。高周疲劳（N>10⁴）以弹性变形为主，低周疲劳（N<10⁴）伴随明显的塑性变形。电机轴系扭振通常属于高周疲劳范畴。
3.2 S-N曲线与疲劳极限
S-N曲线，是疲劳分析的「圣经」。横轴是循环次数N，纵轴是应力幅S。它描述的是：在某个应力水平下，材料能扛多少次循环。
我习惯把S-N曲线分成三段来看：
- 低周区（N<10⁴）：应力高，寿命短，塑性变形明显
- 高周区（10⁴
- 无限寿命区（N>10⁷）：应力低于疲劳极限，理论上不会破坏
疲劳极限，是S-N曲线上的一个拐点。对于钢材，通常10⁷次循环后曲线就平了，这个水平对应的应力就是疲劳极限。但要注意，铝合金、铸铁这些材料没有明显的疲劳极限，曲线会一直往下走。我做电机轴材料选型时，一般优先选有疲劳极限的钢材，设计起来心里有底。
S-N曲线的数学表达式，最常用的是Basquin公式：
σ_a = σ_f' * (2N_f)^b其中：
- σ_a — 应力幅
- σ_f' — 疲劳强度系数
- N_f — 循环次数
- b — 疲劳强度指数（通常-0.05到-0.12之间）
这个公式在Abaqus里可以直接用。我一般会在材料模块里定义好S-N曲线数据，然后交给软件去算。
个人经验：做电机轴系疲劳分析时，别直接用标准S-N曲线。最好根据实际表面加工状态、尺寸效应、应力比做修正。我曾经因为没考虑表面粗糙度，算出来的寿命比实际高了3倍，差点出大事。
3.3 Miner线性累积损伤理论
实际工况下，电机轴的载荷不是恒幅的。今天满负荷，明天半负荷，后天可能还来个冲击。这时候，Miner线性累积损伤理论就派上用场了。
Miner理论的核心思想很简单：每个应力循环都会消耗一点「寿命」，消耗完了就坏了。数学表达是：
D = Σ (n_i / N_i)其中：
- D — 总损伤（D=1时破坏）
- n_i — 第i级应力下的实际循环次数
- N_i — 第i级应力下的允许循环次数（从S-N曲线查得）
举个例子：假设某轴在100MPa下能扛10⁶次，实际跑了5×10⁵次；在80MPa下能扛10⁷次，实际跑了2×10⁶次。那么总损伤D = 5×10⁵/10⁶ + 2×10⁶/10⁷ = 0.5 + 0.2 = 0.7。还没到1，理论上还能继续用。
但我要提醒你，Miner理论有个大前提——损伤是线性累加的，不考虑加载顺序。实际上，先高应力后低应力，和先低应力后高应力，结果是不一样的。我曾经做过对比试验，发现加载顺序对寿命的影响能有20%~30%。所以，Miner理论只能算是个工程近似。
实用建议：在Abaqus里做多轴疲劳分析时，我一般用Miner理论配合雨流计数法。先对时域载荷做雨流计数，得到各级应力幅的循环次数，再代入Miner公式算损伤。这套流程在Fe-Safe和Abaqus里都有现成工具。
3.4 疲劳寿命评估方法
疲劳寿命评估，方法不少。我按自己的使用频率排个序：
- 名义应力法：最经典的方法。用S-N曲线直接算寿命。适合高周疲劳、应力水平不高的情况。电机轴系扭振分析，我90%的情况用这个方法。
- 局部应力-应变法：用ε-N曲线，考虑局部塑性。适合低周疲劳、应力集中严重的区域。比如轴肩根部，我会用这个方法做校核。
- 断裂力学法：基于Paris公式，计算裂纹扩展寿命。适合有初始缺陷的情况。我一般用在焊缝、铸造缺陷的评估上。
- 频域法：把时域载荷转换成功率谱密度，在频域里算疲劳。适合随机振动工况。电机轴系扭振，如果载荷是宽频随机信号，我会用这个方法。
在Abaqus里做疲劳分析，我通常走这个流程：
1. 建立有限元模型，施加边界条件
2. 做静力分析或瞬态分析，提取应力结果
3. 定义材料S-N曲线（或ε-N曲线）
4. 选择疲劳分析方法（名义应力法/局部应变法）
5. 设置载荷谱（恒幅/变幅/随机）
6. 提交计算，查看寿命云图
7. 对危险区域做局部细化，重新校核这个流程看起来简单，但每一步都有坑。比如第二步，应力提取的位置很关键。我习惯在表面节点提取应力，因为疲劳裂纹通常从表面萌生。再比如第五步，载荷谱的准确性直接影响结果。我曾经因为载荷谱少了一个冲击工况，算出来的寿命比实际高了5倍。
避坑指南：做疲劳分析时，千万别忽略平均应力的影响。高平均应力会显著降低疲劳寿命。在Abaqus里，可以用Goodman、Gerber或Soderberg公式做修正。我个人偏好Goodman修正，简单实用，工程上够用。
知识体系总览
下面这张图，是我自己整理的疲劳理论基础框架。你可以把它当作本章的「地图」：
这张图把四个核心知识点串起来了。从破坏机理到S-N曲线，再到累积损伤理论，最后落到评估方法上。你可以在脑子里过一遍：裂纹怎么来的？S-N曲线怎么用？损伤怎么累加？用什么方法算寿命？这四个问题搞清楚了，疲劳分析的基本功就算打牢了。
好了，这一章的内容就到这里。疲劳理论看起来枯燥，但每一条公式、每一个概念，背后都是无数工程师用失败换来的经验。我写这些，也是希望你能少走一些我当年走过的弯路。
公众号:蓝海资料掘金营，微信deep3321