4. 围线积分与网格收敛性分析

在断裂力学仿真中，应力强度因子 K 与 J 积分是判断裂纹是否起裂、扩展的核心参量。Abaqus 通过围线积分（contour integral）在裂纹尖端周围构造若干条环绕路径，对每条路径分别计算 J 积分，并进一步换算为 K。围线积分结果是否可信，取决于两件事：积分区域与路径设置是否合理，以及网格密度是否已经收敛。本节以一个标准紧凑拉伸（CT）试样为例，完整走一遍从建模、设置围线积分、加密网格到判读收敛性的过程。

一、围线积分的物理含义与 Abaqus 实现方式

J 积分定义为围绕裂纹尖端的一条逆时针闭合路径上的能量线积分，其值在理想弹性条件下与路径无关，这正是它可作为断裂参量的根本原因。Abaqus 并不真的让你手画一条路径，而是采用“域积分”等价形式：以裂纹尖端为圆心，按等比例缩放生成一组同心环状单元域，每个域对应一条等效积分路径，称为一个围线（contour）。默认输出 5 条围线，可通过关键字或界面修改。

要让围线积分正常工作，模型必须满足几个硬性条件：裂纹尖端必须是奇异单元或退化的四边形/六面体单元（collapsed element），裂纹面必须真实分开，尖端附近网格需要围绕尖端呈放射状布置。Abaqus/CAE 中通过 Interaction → Special → Crack → Contour integral 定义裂纹前沿，选择裂纹尖端节点与裂纹面，再指定围线数量。若使用 XFEM，则无需预置裂纹面，但本节采用传统尖锐裂纹建模，以便直接观察网格收敛行为。

二、CT 试样几何与材料设置

采用标准 CT 试样，宽度 W=50 mm，厚度 B=12.5 mm，初始裂纹长度 a=25 mm（a/W=0.5），销孔直径 12.5 mm。材料取线弹性，弹性模量 E=210 GPa，泊松比 ν=0.3，用于先验证 K 的收敛性；随后切换为弹塑性（双线性等向强化，屈服 350 MPa，切线模量 2 GPa），观察 J 积分在大变形下的收敛特征。

加载方式采用位移控制：在上销孔施加竖直向上 0.5 mm 的位移，下销孔完全固定。位移控制比力控制更稳定，便于在接近起裂时仍能获得收敛解。输出请求中必须勾选 Contour integral，并指定输出到 .dat 与 .odb。

三、围线积分区域与路径选择

围线数量并非越多越好。第 1 条围线紧贴裂纹尖端，受奇异单元畸变影响最大，数值往往偏大或偏小；最后几条围线若延伸到远场或跨越销孔、边界，又会因路径无关性被破坏而失真。工程上通常取中间若干条围线的平均值作为报告值，并检查首末围线是否偏离均值超过 5%。

在 Abaqus 中，围线半径由程序根据尖端周围单元尺寸自动按比例缩放，用户可通过 *CONTOUR INTEGRAL 关键字中的 RADIUS 参数控制起始半径。建议起始半径取裂纹尖端最小单元边长的 2～3 倍，围线数量取 5～8 条，使最外圈仍落在 K 主导区内（一般要求 r < a/10）。

四、网格密度对 K 与 J 的影响

围线积分对尖端网格极为敏感。若尖端单元过大，奇异场无法被正确捕捉，K 会明显偏低；若单元过小但过渡区太陡，又会因刚度矩阵病态导致 J 值震荡。合理的做法是围绕尖端做放射状网格，径向尺寸按等比递增，环向单元数保持 16～24 个，径向第一层单元尺寸约为 a/100～a/50。

下面给出一组示例性的收敛研究数据（属于教学假设，非实测结果）。固定 a/W=0.5、位移 0.5 mm，仅改变尖端最小单元尺寸 h，读取第 3 条围线的 K 值：

判据：当相邻两级网格的 K 变化小于 2% 时，可认为已进入收敛区间。上例中 M4 与 M5 之间仅差 0.45%，说明 h≈0.06 mm 已足够；继续加密到 M5 收益很小，计算成本却翻倍。工程上通常取 M4 作为报告网格。

五、J 积分与弹塑性收敛的差异

切换到弹塑性材料后，J 积分不再与路径严格无关，因为塑性区尺寸随载荷增大而扩展。此时应关注：塑性区是否被最外圈围线完全包围。若最外圈仍落在塑性区内，J 会随路径号明显漂移。建议在 .dat 中查看各围线 J 值，若 C1～C5 的离散度超过 8%，说明围线区域不足，需要增大围线半径或增加围线数量。

此外，弹塑性下网格收敛比线弹性更慢。示例中当 h 从 0.125 mm 降到 0.0625 mm 时，J 仍变化约 4%，直到 h=0.03 mm 才降到 1.5% 以内。因此弹塑性断裂分析通常要求尖端单元比线弹性更细，且需要保证裂纹尖端塑性区内有至少 3～4 层单元。

六、收敛性研究的操作闭环

完整流程可归纳为：建立 CT 模型 → 定义裂纹与围线积分（5～8 条）→ 提交线弹性计算 → 提取各围线 K → 按 h 逐级加密并记录 → 判断相邻级变化是否小于 2% → 切换弹塑性 → 检查各围线 J 离散度 → 确认塑性区被包围 → 输出报告网格与报告围线均值。每一步都应在 .dat 或 .odb 中留下可复查的数据，而不是凭感觉判断“差不多了”。

最终交付物建议包含一张收敛曲线（K 或 J 对 h 的双对数图）和一张围线分布表（各路径的 K/J 值及相对偏差）。这两份材料能直接说明你的网格与围线设置是否可信，也是后续裂纹扩展分析的基础。