断裂力学 (Fracture Mechanics) — CAE术语解说

很好的问题。普通的强度计算是基于材料完好无缺的前提，检查应力是否在许可值以下。但实际的结构物中总是存在焊接缺陷、疲劳裂纹、腐蚀坑等各种不连续性。

完全正确。断裂力学（Fracture Mechanics）是以"裂纹已经存在"为出发点，定量评估该裂纹是否保持稳定、何时会导致不稳定破裂的学科。例如，当飞机机身发现微小裂纹时，可以判断"在下次检查前，该裂纹是否会扩展到危险尺寸"。这是损伤容限设计的基础。

完全不同。应力集中系数 $K_t$ 是无量纲倍数，表示缺口或孔附近的名义应力增大的倍数。而应力强度因子 $K$ 的单位是 $\mathrm{MPa}\sqrt{\mathrm{m}}$，表示裂纹尖端应力场的"强度"本身。其定义式为：

它代表裂纹的变形模式。正好乘此机会讲讲裂纹模式。

是的，有3种。这对实务工作很重要，需要牢记。

基本上是这样。这是线性弹性断裂力学（LEFM）的破裂判定准则。但这个判定有很重要的前提条件。为了讲清楚，我们来讨论LEFM与EPFM的区别。

LEFM是 Linear Elastic Fracture Mechanics（线性弹性断裂力学），EPFM是 Elastic-Plastic Fracture Mechanics（弹塑性断裂力学）。简单来说，根据裂纹尖端的塑性区大小来判断使用哪一种。

LEFM在"小规模屈服（Small-Scale Yielding: SSY）"条件下有效，即裂纹尖端塑性区尺寸 $r_p$ 相对于裂纹长度 $a$ 或试件尺寸很小。一般用Irwin塑性区尺寸估算公式：

大致如此。但即使是同一材料，薄板会因为平面应力状态导致塑性区增大，低温时可能变成脆性破裂而适用LEFM。压力容器设计规范（ASME BPVC第XI卷或BS 7910）中使用FAD（失效评估图）来同时覆盖LEFM和EPFM两种情况。

直观地说，就是裂纹扩展单位面积时释放的能量。由Rice（1968）定义，是沿环绕裂纹尖端任意路径 $\Gamma$ 的路径积分：

这正是J积分的妙处——它路径独立（path-independent）。无论在裂纹尖端附近还是远处的路径上计算，都得到相同的值。因此在FEM中，即使裂纹尖端特异场的精度有所不足，也能在离裂纹尖端较远的地方精确计算J值。Abaqus、ANSYS等都支持"contour integral"功能，可以自动用多条路径计算J并检验收敛性。

在小规模屈服条件下，两者有这样的关系：

传统FEM处理裂纹时，必须让网格单元边界与裂纹面重合。由于裂纹尖端存在 $1/\sqrt{r}$ 的应力奇异性，需要在尖端部署特异单元（将中间节点移动到1/4点位置）。但每次裂纹扩展都要重新划分网格。在3D问题中这非常繁琐。

完全正确。这是XFEM的最大优势。XFEM在传统FEM形状函数中添加富集函数（enrichment function）：

Abaqus 6.9以后的版本标准支持XFEM，裂纹插入和扩展相对容易设置。ANSYS Mechanical也有Smart Crack Growth功能。开源软件中，Code_Aster支持X-FEM。不过XFEM并非万能，对于3D复杂裂纹分岔或涉及接触的情况仍有困难。这时可以考虑Cohesive Zone Model（CZM）或Phase-field方法。

最重要的是裂纹尖端的网格划分。传统FEM中，需要在裂纹尖端周围建立同心圆放射状的"蜘蛛网网格"，尖端节点合并为一点，采用特异单元。最内层单元尺寸应为裂纹长度的1/20以下。J积分的轮廓线选5～10条，检查外层轮廓的数值是否已收敛。

首先计算每个荷载循环的应力强度因子范围 $\Delta K = K_{\max} - K_{\min}$，然后用Paris法则求裂纹扩展速率：

可以这样总结：