2. 疲劳基础理论与关键参数

在Abaqus中做疲劳与耐久性仿真，第一步不是点“提交作业”，而是先把载荷谱翻译成材料能“听懂”的语言。材料对循环载荷的响应，最终都落在几条曲线和几个参数上：S-N曲线、ε-N曲线、应力比R、平均应力σ_m、疲劳极限σ_e，以及把不同幅值循环拼在一起的损伤累积法则。理解这些量之间的换算关系，才能正确设置Abaqus的疲劳分析模块（如fe-safe接口或直接循环加载）并解释寿命结果。

一、S-N曲线：高周疲劳的“寿命地图”

S-N曲线描述应力幅S_a与失效循环数N_f的关系，通常画在双对数坐标上。对于钢等有明确疲劳极限的材料，曲线在约10⁶～10⁷周次后趋于水平，该水平对应的应力幅即疲劳极限σ_e；铝合金等无疲劳极限材料则持续下降，工程上常取10⁷或5×10⁸周次对应的条件疲劳极限。

在Abaqus/fe-safe中录入S-N曲线时，需要给出若干(S_a, N_f)数据点，并指定插值方式（对数-对数线性插值最常用）。例如某低合金钢的示例数据（教学假设值，非实测）：

S_a (MPa)   N_f (cycles)
  400        1.0e4
  300        5.0e4
  220        2.0e5
  180        1.0e6
  160        5.0e6
  150        1.0e7   ← 疲劳极限附近

注意：S-N曲线通常对应对称循环（R=-1）或某一特定应力比。若实际工况的R不同，必须做平均应力修正，否则寿命会被严重高估或低估。

二、ε-N曲线：低周疲劳与局部应变法

当循环应力接近或超过屈服强度，塑性应变不可忽略，此时用应变幅ε_a比应力幅更合适。ε-N曲线由Basquin（弹性项）和Coffin-Manson（塑性项）叠加而成：

ε_a = σ'_f/E · (2N_f)^b + ε'_f · (2N_f)^c

其中σ'_f为疲劳强度系数，b为疲劳强度指数（约-0.05～-0.12），ε'_f为疲劳延性系数，c为疲劳延性指数（约-0.5～-0.7）。在Abaqus中做低周疲劳时，通常先用弹塑性分析得到危险点的应变历程，再用ε-N曲线和Morrow、SWT等平均应力修正模型计算寿命。ε-N曲线适用于总寿命低于约10⁴～10⁵周次的低周疲劳问题。

三、应力比、平均应力与应力幅

一个循环由最大应力σ_max和最小应力σ_min定义。三个基本参数：

• 应力幅 σ_a = (σ_max - σ_min)/2
• 平均应力 σ_m = (σ_max + σ_min)/2
• 应力比 R = σ_min/σ_max

拉伸平均应力会降低疲劳寿命，压缩平均应力则往往有利。因此当R≠-1时，需要把实际循环等效到R=-1的S-N曲线上。常用修正模型包括：

• Goodman：σ_a/σ_-1 + σ_m/σ_b = 1，偏保守，适合脆性材料；
• Gerber：σ_a/σ_-1 + (σ_m/σ_b)² = 1，适合延性材料；
• Soderberg：用屈服强度σ_y代替σ_b，最保守；
• Morrow：用真实断裂应力σ_f修正，常用于应变-寿命法。

在Abaqus/fe-safe中，这些修正通过“Mean Stress Correction”选项选择。若工况中平均应力随载荷块变化，应逐块修正，而不是用全局平均值。

四、疲劳极限与条件疲劳极限

疲劳极限σ_e是材料在无限次循环下不失效的最大应力幅。对钢材，若S-N曲线在10⁷周次后水平，该水平即σ_e；对铝合金、镁合金等，曲线无水平段，只能定义条件疲劳极限，例如10⁷周次对应的应力幅。在Abaqus/fe-safe中，若材料数据只给到10⁷，软件默认将之后视为无限寿命，此时应确认该假设是否与设计规范一致。

疲劳极限还受应力比、表面质量、尺寸效应、加载频率等影响。工程上常用修正系数：σ_e,实际 = σ_e · k_表面 · k_尺寸 · k_载荷。这些系数应在材料数据录入前完成，而不是在Abaqus后处理中临时调整。

五、Miner线性累积损伤理论及其适用条件

实际载荷谱通常包含多个幅值级别。Miner理论假设：每个循环造成的损伤独立、可线性叠加，总损伤D = Σ(n_i/N_fi)。当D≥1时判定失效。例如某零件经历三个载荷块（教学示例）：

Miner理论的适用条件需要特别注意：

• 载荷顺序无关假设：实际中高-低顺序比低-高顺序更危险，Miner无法区分；
• 无过载迟滞：单个大过载可能引入残余压应力，延缓后续裂纹扩展，Miner会低估寿命；
• 损伤线性：小载荷循环在Miner中仍贡献损伤，但实际可能低于疲劳极限而不造成损伤，需设置截止限；
• 无交互作用：不同幅值之间的相互影响被忽略。

六、在Abaqus中落地这些参数

典型流程是：先做静力或弹塑性分析得到危险点应力/应变场，再把结果导入fe-safe或使用Abaqus内置的疲劳模块。材料卡片中录入S-N或ε-N曲线，选择平均应力修正模型，定义载荷谱和应力比，设置Miner累积与截止限，最后提取寿命云图和损伤分布。判断结果时，重点看三个量：危险点寿命是否满足设计目标、损伤是否集中在预期位置、以及D值是否接近1但未超过。若D远小于1而实际工况更严酷，应检查载荷谱是否被过度简化。