4. S-N曲线与疲劳寿命预测
各位同学，大家好。今天我们来聊聊复合材料疲劳分析里最基础、也最核心的一块——S-N曲线。说白了，S-N曲线就是材料疲劳性能的“身份证”。你拿到一个复合材料结构，想知道它能在交变载荷下撑多久，就得靠它。
我个人习惯，在讲理论之前，先跟大家说说这东西到底怎么来的。S-N曲线，全称是应力-寿命曲线。横坐标是循环次数N（通常取对数），纵坐标是应力幅值S。一条曲线，就描述了材料在不同应力水平下能活多久。
核心要点：S-N曲线是疲劳分析的基石。没有它，后面所有的寿命预测都是空中楼阁。
4.1 S-N曲线的获取方法
获取S-N曲线，主要有两条路：一条是老老实实做试验，另一条是用经验公式估算。我两种都干过，各有各的坑。
4.1.1 试验法
试验法最直接。拿一批标准试件，在不同的应力水平下做疲劳试验，记录下每个试件断裂时的循环次数。然后把数据点画在图上，拟合出一条曲线。
具体步骤是这样的：
- 确定应力水平：一般取5-8个不同的应力幅值。最高应力水平通常取静强度的70%-80%，最低取30%-40%。
- 进行疲劳试验：每个应力水平下至少做3-5个试件。为什么？因为疲劳数据分散性大，你懂的。
- 记录失效循环数：记录每个试件断裂时的N值。
- 数据拟合：用最小二乘法拟合S-N曲线。常用的拟合公式是幂函数形式：S^m * N = C。
我记得有一次做碳纤维复合材料的疲劳试验，数据点散得一塌糊涂。同一个应力水平下，有的试件撑了10万次，有的才1万次就断了。后来发现是试件加工质量不一致，边缘有毛刺。嗯，这里要注意，试件质量直接影响试验结果。
我的经验：做疲劳试验前，一定要先做静力试验，摸清材料的静强度。这样选应力水平才有依据。另外，试件数量别省，至少每个应力水平5个，否则数据没法看。
4.1.2 经验公式法
试验法虽然准，但费时费力。有时候项目周期紧，没时间做全套试验，就得靠经验公式估算。
常用的经验公式有：
- Basquin公式：S = S_f * (2N)^b。这是最经典的公式，适用于高周疲劳。S_f是疲劳强度系数，b是疲劳强度指数。
- Manson-Coffin公式：这个主要针对低周疲劳，考虑塑性应变。公式是：ε_a = σ_f/E * (2N)^b + ε_f * (2N)^c。
- 复合材料专用公式：比如Hashin-Rotem公式，考虑了纤维和基体的不同失效模式。
你想想看，经验公式的好处是快，但缺点也很明显——精度有限。不同材料体系、不同铺层顺序，公式里的参数都得重新标定。说白了，经验公式只能给你一个大概的参考，不能完全替代试验。
警告：千万别直接拿文献里的参数套用。我曾经吃过这个亏，用碳纤维/环氧树脂的参数去算玻璃纤维/聚酯树脂的疲劳寿命，结果差了整整一个数量级。不同材料体系，参数必须重新标定。
4.2 应力比与平均应力的影响
实际工程中的疲劳载荷，很少是标准的对称循环（R=-1）。更多时候是拉-拉疲劳（R>0）或者拉-压疲劳（R<0）。这时候，应力比R和平均应力σ_m的影响就出来了。
应力比R的定义很简单：R = σ_min / σ_max。R=-1是对称循环，R=0是脉动循环，R=0.1是典型的拉-拉疲劳。
平均应力σ_m = (σ_max + σ_min) / 2。平均应力越大，疲劳寿命越短。为什么？因为平均应力相当于给材料施加了一个恒定的拉伸载荷，加速了损伤累积。
处理平均应力影响，常用的方法有：
- Goodman修正：σ_a = σ_{-1} * (1 - σ_m/σ_b)。σ_{-1}是对称循环下的疲劳极限，σ_b是静强度。
- Gerber修正：σ_a = σ_{-1} * (1 - (σ_m/σ_b)^2)。这个比Goodman更保守一些。
- Soderberg修正：σ_a = σ_{-1} * (1 - σ_m/σ_y)。σ_y是屈服强度，这个最保守。
我个人习惯，在复合材料里常用Goodman修正。但要注意，复合材料的疲劳行为比金属复杂得多，这些修正公式只能作为近似。我在做风电叶片疲劳分析时，就发现Goodman修正对碳纤维复合材料偏保守，对玻璃纤维复合材料又偏危险。所以，最好还是用试验数据来验证。
关键点：应力比和平均应力对疲劳寿命的影响不可忽视。在工程应用中，一定要根据实际载荷谱来选择合适的修正方法。
4.3 疲劳寿命预测模型
有了S-N曲线，知道了载荷谱，接下来就是预测寿命了。这里最经典的就是Miner线性累积损伤理论。当然，也有更复杂的非线性模型。
4.3.1 Miner线性累积损伤
Miner理论的核心思想很简单：每个循环造成的损伤是独立的，可以线性叠加。公式是：D = Σ (n_i / N_i)。其中n_i是第i级应力水平下的实际循环数，N_i是该应力水平下对应的疲劳寿命。
当D=1时，材料就失效了。
Miner理论的优点是简单、好用。但缺点也很明显——它忽略了载荷顺序效应。你想想看，先高应力后低应力，和先低应力后高应力，损伤累积过程能一样吗？实际上，高应力先作用会产生微裂纹，后续低应力反而会加速裂纹扩展。Miner理论完全没考虑这个。
我记得有一次做直升机旋翼桨叶的疲劳分析，用Miner理论算出来的寿命是5000小时，但实际台架试验只撑了3500小时。后来分析发现，是载荷谱中的高应力峰值造成了额外的损伤，Miner理论低估了。
避坑指南：我曾经在风电叶片疲劳分析中，直接用Miner理论算寿命，结果跟实测差了30%。后来改用修正的Miner理论（比如考虑门槛值），精度才上来。所以，Miner理论可以用，但一定要结合工程经验进行修正。
4.3.2 非线性累积损伤模型
为了克服Miner理论的不足，学者们提出了各种非线性累积损伤模型。常用的有：
- Corten-Dolan模型：考虑了损伤的加速效应。公式是：D = Σ (n_i / N_i)^a，其中a是材料常数，通常大于1。
- Manson-Halford模型：基于损伤曲线法，考虑了载荷顺序效应。公式比较复杂，但精度更高。
- 基于剩余强度的模型：认为材料的剩余强度随着循环次数增加而下降，当剩余强度等于最大应力时，材料失效。
非线性模型虽然精度高，但参数多，计算复杂。在工程应用中，我一般这样选：
- 如果载荷谱比较简单（比如等幅疲劳），用Miner理论就够了。
- 如果载荷谱复杂（比如变幅疲劳、随机载荷），建议用Corten-Dolan模型。
- 如果对精度要求极高（比如航空结构），那就得上Manson-Halford模型。
总结：没有完美的模型，只有最适合的模型。工程应用的关键，是在精度和效率之间找到平衡点。
4.4 本章知识体系
为了让大家更直观地理解本章的知识结构，我画了一张图：
这张图把本章的核心内容串起来了。从S-N曲线的获取，到应力比和平均应力的修正，再到寿命预测模型的选择，每一步都环环相扣。你想想看，少了任何一环，疲劳分析都不完整。
最后说一句：疲劳分析不是纯理论，更不是纯试验。它是理论和经验的结合。多积累工程案例，多对比试验数据，你才能慢慢找到感觉。