📚 课程目录第2章断裂力学基础与Abaqus实现

2. 断裂力学基础与Abaqus实现

线弹性断裂力学（LEFM）的核心任务，是把裂纹尖端附近那个理论上趋于无穷的应力场，用一组有限、可比较、可测量的参数来刻画。Abaqus 并不直接"算应力"，而是围绕这些参数提供了一套完整的求解与后处理工具。理解参数定义与 Abaqus 实现之间的对应关系，是后续所有裂纹扩展仿真的前提。

一、裂纹尖端的应力场与应力强度因子

对 I 型（张开型）裂纹，取裂纹尖端为原点、裂纹面沿 x 轴、裂纹前沿沿 z 轴，极坐标 (r, θ) 下的应力场可写为：

σ_ij(r,θ) = K_I / sqrt(2πr) · f_ij(θ) + 高阶项

其中 K_I = lim(r→0) [ σ_yy · sqrt(2πr) ]   （θ = 0 方向）

K_I 即 I 型应力强度因子，单位 MPa·mm^1/2 或 MPa·m^1/2。它把"载荷大小"和"裂纹几何"合并成一个量：K_I 越大，尖端场越强。当 K_I 达到材料的断裂韧度 K_IC 时，裂纹起裂。II 型（滑开）与 III 型（撕开）分别对应 K_II、K_III，混合模式下三者共同决定起裂。

需要强调：K 的解析解只对少数规则几何存在（如无限大板中心穿透裂纹 K_I = σ√(πa)）。工程结构几乎都依赖数值方法，这正是 Abaqus 围线积分的用武之地。

二、能量释放率 G 与 J 积分

能量释放率 G 定义为裂纹扩展单位面积所释放的弹性能：G = −∂Π/∂A，单位 mJ/mm²（即 N/mm）。在线弹性、无塑性耗散的条件下，G 与 K 存在一一对应关系：

平面应力：G = K² / E
平面应变：G = K² (1 − ν²) / E

J 积分是 Rice 提出的围线积分，沿包围裂纹尖端的任意逆时针路径 Γ 计算：

J = ∫_Γ ( W dy − T_i · ∂u_i/∂x ds )

其中 W 为应变能密度，T_i 为路径上的面力，u_i 为位移。J 的关键性质是路径无关性——只要路径包围尖端且不穿过其他塑性区，J 值相同。在线弹性条件下 J = G；在弹塑性条件下 J 仍可作为扩展驱动力，这是它比 K 更通用的原因。Abaqus 的围线积分正是基于这一路径无关性，用多条同心围线分别计算并相互校核。

三、Abaqus 中的围线积分设置

Abaqus/Standard 通过 *CONTOUR INTEGRAL 关键字或 CAE 的 History Output Request 激活围线积分。关键设置包括：

裂纹前沿定义：需指定裂纹尖端节点与裂纹扩展方向。在 2D 中，裂纹前沿是一个点；在 3D 中是一条线，需沿前沿逐点定义。
围线数量：默认 5 条，建议取 5～10 条。第一条围线应紧贴尖端但避开高度畸变单元；最后一条应远离塑性区。
奇异性单元：裂纹尖端周围通常采用 1/4 节点位移的二次单元（CPS6、CPE8 等），使应变场具有 1/√r 奇异性，从而更准确地捕捉 K 和 J。
对称性：若模型具有对称性，可只建半模型，在对称面上施加法向约束，围线积分仍可正常计算。

一个典型的 2D 平板中心裂纹模型（宽 100 mm、高 200 mm、中心裂纹长 2a = 20 mm、远场拉应力 100 MPa、E = 210 GPa、ν = 0.3）在 Abaqus 中的操作流程如下：

# 1. 创建部件与材料
#    材料：E=210000 MPa, ν=0.3, 弹性
# 2. 装配并划分网格
#    裂纹尖端附近使用 CPS6 三角形单元，1/4 节点
# 3. 定义分析步
#    *STATIC, 施加 100 MPa 拉应力
# 4. 定义裂纹
#    *CONTOUR INTEGRAL, CRACK NAME=CRACK1, CONTOURS=8
#    裂纹尖端节点集：TIP_NODE
#    裂纹扩展方向：X 方向
# 5. 提交计算
# 6. 后处理查看 History Output：
#    J 积分、K_I、K_II、G

计算完成后，在 History Output 中会看到 J 和 K 随分析步的变化。对于线弹性问题，8 条围线的 J 值应基本一致（差异通常小于 2%～5%）。若第一条围线结果明显偏离，说明尖端网格过粗或单元畸变严重；若最后一条围线偏离，说明围线已进入边界或塑性区，需调整围线半径。

四、结果判读与常见问题

以示例参数估算：无限大板中心裂纹 K_I = σ√(πa) = 100 × √(π × 10) ≈ 560 MPa·mm^1/2。有限宽板需乘以修正系数 F(a/W)，当 a/W = 0.1 时 F ≈ 1.01，故 K_I ≈ 566 MPa·mm^1/2。Abaqus 围线积分结果应接近该值，偏差主要来自网格密度和围线位置。

实用检查清单：

围线是否全部位于弹性区内？若材料已屈服，J 仍可用，但 K 不再严格有效。
多条围线的 J 值是否收敛？取第 3～6 条的平均值通常最稳定。
裂纹尖端单元是否为 1/4 节点奇异单元？普通单元会低估 K。
2D 问题注意平面应力与平面应变的区别，E′ 不同，G 与 K 的换算关系不同。
3D 裂纹前沿各点的 K 沿厚度方向可能不同，需逐点提取并绘制分布曲线。

常见误区：把 J 积分当作"只能用于弹性"的量。实际上 J 在弹塑性条件下仍是有效的扩展驱动力参数，只是此时 J ≠ G，且需要保证围线不穿过卸载区。Abaqus 在弹塑性分析中同样输出 J，但应结合塑性区大小判断其有效性。

掌握 K、G、J 三者的定义、换算关系与 Abaqus 围线积分的设置方法后，就可以进一步讨论裂纹起裂判据与扩展方向，这也是下一阶段混合模式断裂分析的基础。