ANSYS 裂纹仿真：疲劳裂纹扩展与断裂韧性评估

本项目针对 ASTM E399 标准紧凑拉伸（CT）试样开展断裂力学仿真。材料为某中碳钢，试样厚度 25 mm，裂纹初始长度 15 mm，外载按疲劳循环加载。模型在 ANSYS Mechanical 中建立二维平面应变模型，裂纹尖端采用奇异单元（1/4 节点）并计算应力强度因子 K_I。

案例分析

裂纹尖端应力场显示，在 150 kN 外载下 K_I 达到 28.4 MPa·m⁰·⁵，塑性区尺寸约 2.5 mm。材料的平面应变断裂韧性 K_IC 为 62 MPa·m⁰·⁵，因此当前裂纹长度下的安全系数约为 2.18。疲劳裂纹扩展按 Paris 定律计算，参数 C=2.6×10⁻¹²、指数 m=3.1；从初始 15 mm 扩展到临界裂纹 40 mm 约需 10⁷ 循环，检查间隔建议取 4.0×10⁶ 循环。

模型细节、验证与工程意义

模型在 ANSYS Mechanical 2023 R1 中建立，裂纹尖端采用 1/4 节点奇异单元（PLANE183），其余区域为六面体网格。载荷按疲劳试验机正弦波形施加，R=0.1。验证时将 K_I 与 ASTM E399 标定公式对比，偏差 <3%；Paris 曲线与实测 da/dN 数据对比，寿命预测偏差在 2 倍分散带内。模型可用于剩余寿命评估与检查周期制定。

总结教训

第一个教训是奇异单元方向必须指向裂纹尖端，否则 K_I 计算会偏差 10% 以上。第二个教训是塑性区尺寸需满足小范围屈服条件（<2.5% 韧带），本项目 2.5 mm 塑性区在 25 mm 韧带内满足。第三个教训是 Paris 定律参数对寿命预测极敏感，C 值 2 倍变化可使寿命估算相差 3–5 倍，必须按实测 da/dN 曲线标定。

项目经验

1. ANSYS 裂纹仿真建议用奇异单元处理裂尖 1/√r 应力奇异性。
2. 标准紧凑拉伸试样厚度 25 mm，初始裂纹 15 mm。
3. 150 kN 外载下 K_I = 28.4 MPa·m⁰·⁵。
4. 材料断裂韧性 K_IC = 62 MPa·m⁰·⁵，安全系数约 2.18。
5. 裂纹尖端塑性区尺寸约 2.5 mm。
6. Paris 定律参数：C=2.6×10⁻¹²，m=3.1。
7. 临界裂纹长度 40 mm，对应剩余寿命约 10⁷ 循环。
8. 建议检查间隔：4.0×10⁶ 循环。
9. 奇异单元方向必须指向裂尖，否则 K_I 偏差 >10%。
10. 小范围屈服条件要求塑性区 <2.5% 韧带宽度。
11. da/dN 曲线必须实测标定，不能直接用文献值。
12. 裂纹尖端奇异单元采用 1/4 节点 PLANE183。
13. 载荷比 R=0.1，符合疲劳试验标准。
14. 本项目方法可迁移到压力容器、轨道交通部件的裂纹寿命评估。

补充经验

1. 三维模型可进一步考虑厚度方向的约束效应与表面裂纹形貌。
2. 过载迟滞效应需用 Wheeler/Willenborg 模型修正 Paris 预测。
3. 裂纹尖端钝化会延缓启裂，Paris 区预测需用有效裂纹长度修正。
4. 腐蚀环境会显著降低 K_IC，疲劳寿命需按环境修正。

复现参数速查（案例图硬参数）

更多相关仿真经验，可参考本站的 ANSYS 仿真案例栏目。