📚 课程目录第1章断裂力学基础与Abaqus仿真概览

1. 断裂力学基础与Abaqus仿真概览

断裂力学要回答的核心问题只有一个：含裂纹的构件在给定载荷下是否安全，裂纹会不会扩展、沿什么方向扩展、扩展多快。围绕这个问题，线弹性断裂力学（LEFM）建立了三个彼此等价的核心参量——应力强度因子 K、能量释放率 G 和 J 积分。理解它们的定义、适用条件和相互换算关系，是后续在 Abaqus 中正确建模、正确提取结果、正确判断裂纹行为的前提。

一、三个核心参量的物理含义

应力强度因子 K 描述裂纹尖端附近应力场的奇异性强度。以 I 型（张开型）裂纹为例，裂纹尖端附近的应力场可写为 σ_ij = K_I / √(2πr) · f_ij(θ)，其中 r 为到裂尖的距离，θ 为相对裂纹面的角度。K_I 的量纲是应力×√长度，常用单位为 MPa·√mm 或 MPa·√m。K 的数值取决于远场应力、裂纹尺寸与几何构型，典型形式为 K_I = Y·σ·√(πa)，Y 为几何修正因子。当 K_I 达到材料的断裂韧度 K_IC 时，裂纹发生失稳扩展。

能量释放率 G 从能量角度描述裂纹扩展单位面积所释放的弹性能，G = −∂Π/∂A，Π 为系统总势能，A 为裂纹面积。对线弹性材料，G 与 K 存在直接换算关系：平面应力状态下 G = K²/E，平面应变状态下 G = K²(1−ν²)/E。这一关系说明 K 与 G 并非独立判据，而是同一物理过程的两种表述。

J 积分 是围绕裂纹尖端的一条与路径无关的围线积分：J = ∫_Γ (W dy − T_i ∂u_i/∂x ds)，其中 W 为应变能密度，T_i 为面力分量，u_i 为位移分量。J 积分的价值在于它同时适用于线弹性和弹塑性材料，在小变形、单调加载条件下，J 积分可作为弹塑性断裂的控制参量，其临界值 J_IC 用于判断延性材料的起裂。

三个参量的适用边界要分清：K 与 G 仅适用于线弹性或小范围屈服；J 积分可扩展到弹塑性，但要求单调加载、无卸载，且裂纹尖端塑性区尺寸远小于构件特征尺寸。

二、Abaqus 断裂仿真模块的整体架构

Abaqus 并不存在一个名为“断裂模块”的独立入口，断裂仿真能力分散在多个功能组件中，按求解方式可分为三大类：

围线积分法（Contour Integral）：在 Interaction 模块中定义裂纹尖端（Crack Front），通过 History Output 请求 K_I、K_II、K_III、J 积分和能量释放率。这是最经典、最直接的方法，适用于二维和三维、线弹性与弹塑性。
VCCT（虚拟裂纹闭合技术）：基于 Irwin 裂纹闭合假设，通过节点力和节点位移差计算能量释放率，适用于界面裂纹和复合材料分层问题，在 Interaction 模块中定义 Debond 接触或直接使用 VCCT 关键字。
XFEM（扩展有限元）：通过富集函数描述裂纹，裂纹可以穿过单元内部，无需网格与裂纹面共形，适合模拟裂纹沿任意路径扩展。在 Interaction 模块中定义 XFEM 裂纹，配合损伤起始与演化准则使用。

此外，Cohesive 单元（内聚力单元）用于模拟界面脱粘和胶层断裂，属于另一种独立的建模路径。选择哪种方法，取决于材料行为（脆性/延性）、裂纹路径是否已知、是否需要自动扩展以及计算成本。

三、从参量到仿真：一个可执行的入门算例

下面用一个带中心穿透裂纹的有限宽板拉伸算例，把上述参量与实际操作串起来。板宽 W = 100 mm，半裂纹长度 a = 10 mm，板厚 t = 5 mm，材料弹性模量 E = 210000 MPa，泊松比 ν = 0.3，远场拉伸应力 σ = 100 MPa。几何修正因子按有限宽板公式取 Y ≈ 1.12（近似值，用于教学演示）。

理论估算：
K_I = Y · σ · √(πa)
    = 1.12 × 100 × √(π × 10)
    ≈ 1.12 × 100 × 5.605
    ≈ 627.8 MPa·√mm

平面应力下：
G = K² / E = 627.8² / 210000 ≈ 1.877 N/mm
J = G ≈ 1.877 N/mm

在 Abaqus 中实现这一算例的关键操作如下：

在 Part 模块建立二维平面应力板，用 Partition 在裂纹位置切分，使裂纹面成为独立边。
在 Property 模块赋予线弹性材料，E = 210000 MPa，ν = 0.3。
在 Assembly 模块创建实例，Step 模块建立 Static, General 分析步，打开几何非线性（NLGEOM = ON）并非必需，但建议开启以观察大变形影响。
在 Interaction 模块选择 Special → Crack → Assign Seam，指定裂纹面；再选择 Special → Crack → Create，定义裂纹尖端（Crack Front），并指定裂纹扩展方向 q 向量。
在 History Output 中请求 Contour Integral 输出，围线数量建议取 5～10 条，观察 J 积分是否路径无关。
在 Load 模块施加 100 MPa 拉应力，边界条件约束板端刚体位移。
提交计算后，在 Visualization 模块查看 History Output，读取第 3～5 条围线的 J 积分值，若数值趋于稳定，说明路径无关性满足。

围线积分结果必须检查路径无关性。如果不同围线的 J 值差异超过 5%，通常说明裂纹尖端网格过粗、围线穿过塑性区或接触面，需要细化裂尖网格或调整围线半径。

四、结果判读与常见误区

计算完成后，把 Abaqus 输出的 K_I 与理论估算值 627.8 MPa·√mm 对比。若仿真值在 600～650 MPa·√mm 范围内，说明建模合理；若偏差超过 10%，应检查裂纹尖端网格质量、围线定义和载荷施加方式。对于弹塑性材料，K 不再适用，应直接读取 J 积分，并与材料的 J_IC 比较。