疲劳有限元分析

机械零件很少是”一下压垮”的，大多数是在循环载荷下慢慢萌生裂纹、扩展、最后疲劳断裂。我这个项目客户的一批带圆孔平板，要在交变载荷下长期服役，他想知道疲劳寿命分布和裂纹最可能从哪开始。我做的 疲劳有限元分析 把 S-N 曲线和寿命云图结合起来给了答案。

疲劳的物理起点是循环塑性。哪怕名义应力低于屈服，局部应力集中处也会反复进入微小塑性，每一循环积累一点损伤，直到裂纹萌生。所以疲劳分析的第一步是输入材料 S-N 曲线：我给了钢的 S-N 数据，高周区应力幅随循环次数 N 增加而下降，最后趋近一个 endurance limit（耐久极限）S_e ≈ 220 MPa——意思是在 10⁶ 次循环以下、应力幅低于这个值，理论上可无限寿命。

第二步是把静力应力的 von Mises 场作为疲劳输入。我在 ANSYS nCode DesignLife 里做，载荷用的是恒幅正弦载荷（简单但能说明问题）。结果用 log₁₀(N_f) 云图表示：圆孔边缘应力集中最明显，疲劳寿命最低（红色区域 N_f 只有 10⁴ 量级），远离孔的区域蓝色（>10⁷）。这个分布和客户后期疲劳试验的裂纹起始位置完全一致，验证了模型。

判读逻辑要讲清：疲劳寿命云图不是”应力图”，而是”到失效还差多少循环”的图。红色=危险，蓝色=安全。客户据此把孔边做了一道加强筋，最低寿命从 10⁴ 拉到 10⁵ 量级，设计就稳了。

工程里载荷往往不是恒幅。真实工况是变幅或随机载荷，这时要先做雨流计数（rainflow counting）把复杂载荷谱拆成一系列恒幅循环，再用 Miner 线性累积损伤准则 Σ(n_i/N_i)=1 判失效。我这个项目作为演示用了恒幅，但报告里给了变幅的处理流程。

平均应力修正也是个关键点。上面 S-N 曲线是基于对称循环（平均应力为 0）的，实际载荷常有平均拉应力，会显著降低疲劳寿命。修正方法有三：Goodman（保守，直线）、Gerber（抛物线，较准）、Soderberg（最保守）。我这个项目材料是钢材，用 Goodman 修正比较稳。这里踩过坑：一开始忘了做平均应力修正，算出来的寿命比试验高了一个数量级，加上修正后才对上。

踩坑记录：第一，应力集中系数。圆孔的理论应力集中因子 K_t 如果不从应力场直接读、而是用经验公式近似，会引入误差，我直接用 FEA 结果提取孔边峰值应力。第二，网格。寿命云图对孔边网格极敏感，孔边必须加密到能分辨应力梯度，我孔边单元尺寸 0.2 mm。第三，表面状态。疲劳对表面粗糙度、残余应力极敏感，S-N 曲线若取自抛光试样，实际机加工件寿命要打折，我加了表面系数。第四，低周 vs 高周。我这个项目是高周疲劳（应力低、循环多），若是低周（塑性应变主导）要用基于应变的 Coffin-Manson 模型，方法完全不同。

从我的工程经验看，疲劳有限元分析 的核心输入是”载荷谱”和”S-N/疲劳参数”，这两个错了后面全错。我给客户的建议是：先拿到真实载荷谱再算，平均应力修正和表面系数一个都不能省。更多疲劳分析，见 [疲劳有限元分析](https://www.keyanxueshu.com/category/fea/)；需要变幅载荷分析，[疲劳有限元分析](https://www.keyanxueshu.com/) 上有入口。