疲劳裂纹扩展分析方法
1. 疲劳失效概述
疲劳是机械设备最常见的失效模式之一，约占失效案例的80%以上。疲劳裂纹扩展分析是损伤容限设计的重要组成部分。
2. 疲劳裂纹扩展的三个阶段
2.1 阶段I：裂纹萌生
- 微观短裂纹阶段
- 长度：几个晶粒尺度
- 扩展速率很低
2.2 阶段II：稳定扩展
- 宏观裂纹稳定扩展
- Paris公式适用
- 占据大部分寿命
2.3 阶段III：失稳扩展
- 快速扩展
- 达到断裂韧性$K_{IC}$
- 最终断裂
3. Paris公式
Paris-Erdogan公式描述了阶段II的裂纹扩展速率：
$$ \frac{da}{dN} = C(\Delta K)^m $$
其中：
- $a$：裂纹长度
- $N$：循环次数
- $\Delta K$：应力强度因子范围
- $C, m$：材料常数
3.1 应力强度因子范围
$$ \Delta K = K_{max} - K_{min} = \Delta \sigma \sqrt{\pi a} \cdot Y $$
其中：
- $\Delta \sigma$：应力范围
- $Y$：几何修正因子
4. Walker公式
考虑应力比$R$的影响：
$$ \frac{da}{dN} = C[\Delta K(1-R)^{1-\gamma}]^m $$
其中$\gamma$为材料常数，对于钢材通常取0.5～0.7。
5. 疲劳寿命预测
5.1 积分法
对Paris公式积分：
$$ N = \int_{a_0}^{a_f} \frac{da}{C(\Delta K)^m} $$
其中：
- $a_0$：初始裂纹长度
- $a_f$：最终裂纹长度
5.2 数值积分
实际工程中，$\Delta K$随$a$变化，需采用数值积分：
$$ N \approx \sum_{i=1}^{n} \frac{\Delta a_i}{(da/dN)_i} $$
6. 影响因素
6.1 应力比效应
- $R = \sigma_{min}/\sigma_{max}$
- $R$越大，扩展速率越快
6.2 载荷谱
- 恒幅载荷
- 变幅载荷：需考虑载荷相互作用效应
6.3 环境因素
- 腐蚀环境
- 高温蠕变
- 低温脆化
7. 工程应用实例
7.1 压力容器焊缝裂纹
已知条件：
- 材料：Q345R
- 初始裂纹：$a_0 = 3$ mm
- 设计压力波动：$\Delta P = 2$ MPa
- Paris常数：$C = 6.9 \times 10^{-12}$, $m = 3.0$（MPa·m^{1/2}单位）
计算步骤：
-
计算应力范围： $$ \Delta \sigma = \frac{PR}{t} $$
-
计算$\Delta K$： $$ \Delta K = \Delta \sigma \sqrt{\pi a} \cdot Y $$
-
代入Paris公式求解$N$
7.2 结果分析
- 预测疲劳寿命
- 确定检测周期
- 评估是否需要修复
8. BS 7910标准
8.1 疲劳评估流程
- 确定初始缺陷尺寸
- 确定应力历程
- 计算应力强度因子
- 选择裂纹扩展模型
- 积分求解寿命
- 失效评定图（FAD）校核
8.2 安全系数
- 裂纹尺寸：安全系数2
- 寿命：安全系数3～10
9. 有限元分析
9.1 FEA在疲劳分析中的作用
- 计算复杂几何的$K$值
- 考虑多轴应力状态
- 模拟裂纹扩展路径
9.2 常用软件
- ANSYS SMART Crack
- ABAQUS XFEM
- FRANC3D
- Zencrack
10. 总结
疲劳裂纹扩展分析是现代损伤容限设计的核心内容，结合断裂力学理论和有限元分析，可以有效预测结构的剩余寿命，为维护决策提供科学依据。
相关阅读：LBB技术