3. 线弹性断裂力学参数提取实战
本节以中心裂纹拉伸试样（CCT，Center Cracked Tension）为对象，完整走一遍 Abaqus/Standard 中轮廓积分（Contour Integral）的建模、提交与结果提取流程，并把数值结果与线弹性断裂力学理论解逐项对比。轮廓积分法在 Abaqus 中同时给出应力强度因子 K 与 J 积分，是线弹性断裂力学参数提取最常用的手段。
一、问题定义与理论解
取一块宽 2W、高 2H、厚 t 的薄板，中心贯穿裂纹长度为 2a，两端施加均布拉伸应力 σ。本节示例参数取：半宽 W = 50 mm，半高 H = 100 mm，半裂纹长 a = 10 mm，板厚 t = 5 mm，弹性模量 E = 210000 MPa，泊松比 ν = 0.3，远端均布应力 σ = 100 MPa。平面应力假设。
对无限大板中心裂纹，应力强度因子理论解为 KI = σ√(πa)。有限宽板需引入修正系数 F(α)，其中 α = a/W：
K_I = σ · √(πa) · F(α)
F(α) = √(sec(πα/2)) （Isida 近似，适用于 α 较小）
本例：α = 10/50 = 0.2
F(0.2) = √(sec(π·0.2/2)) = √(sec(0.31416)) ≈ √(1.0515) ≈ 1.0254
K_I = 100 × √(π × 10) × 1.0254 ≈ 100 × 5.60499 × 1.0254 ≈ 574.7 MPa·√mm
平面应力下 J 与 K 的关系：
J = K_I² / E = 574.7² / 210000 ≈ 1.573 N/mm注意单位一致性：Abaqus 中若长度用 mm、力用 N，则应力为 MPa，K 的单位为 MPa·√mm，J 的单位为 N/mm（即 mJ/mm²）。这是后续对比时最容易出错的地方。
二、建模要点与轮廓积分设置
由于对称性，可只建 1/4 模型以节省计算量，但轮廓积分要求裂纹面与积分路径完整，1/4 模型需要在对称面上施加对称边界条件，并保证裂纹尖端位于对称轴上。为便于教学演示，本节采用全模型，二维平面应力单元 CPE4R，裂纹尖端附近使用退化四边形单元（Collapsed Element）或直接使用常规单元配合足够细密的网格。
关键操作顺序如下：
- 在 Part 模块中建立带中心裂纹的平板，裂纹用一条线表示，两侧面为自由面。
- 在 Property 模块定义线弹性材料，E = 210000 MPa，ν = 0.3，并创建截面赋给平板。
- 在 Assembly 模块创建独立实例，进入 Step 模块建立 Static, General 分析步，打开几何非线性开关（NLGEOM = ON）对线弹性问题影响很小，但建议保持默认关闭以匹配理论解。
- 在 Interaction 模块选择 Special → Crack → Assign Seam，将裂纹线定义为 seam，使裂纹面在网格中分离。
- 在 History Output Request 中勾选 J-integral，并指定积分路径数量（Number of contours），一般取 5～10 条。
- 在 Load 模块施加远端均布压力或拉力，本例在上下端面施加 100 MPa 拉应力。
- 在 Mesh 模块对裂纹尖端区域加密，建议尖端附近单元尺寸约为 a/20 或更小。
三、轮廓积分路径与结果提取
Abaqus 的轮廓积分默认以裂纹尖端为圆心，自动生成一系列围绕尖端的环形路径。每条路径对应一个 J 积分值，同时可输出 KI、KII、KIII。对于纯 I 型问题，KII 与 KIII 应接近零，可作为结果合理性的辅助判断。
提交作业后，在 Visualization 模块中通过 Result → History Output 查看 J 积分随轮廓编号的变化。若各条轮廓的 J 值趋于稳定，说明路径无关性得到满足，结果可信。随后用 J 反算 K：平面应力下 K = √(J·E)，平面应变下 K = √(J·E/(1−ν²))。本例为平面应力，故 K = √(1.573 × 210000) ≈ 574.7 MPa·√mm，与理论解一致。
四、结果判读与常见问题
从上面的收敛曲线可以看到，第 1 条轮廓的 J 值略高于理论值，这是因为最靠近裂纹尖端的路径受应力奇异性影响较大，数值误差更明显；第 3 条轮廓之后 J 值基本稳定在 1.573 N/mm 附近，与理论解吻合良好。工程实践中建议取第 3～5 条轮廓的平均值作为报告值。
若结果与理论解偏差超过 5%，可按以下顺序排查：
- 单位制不一致：检查长度、力、应力是否统一为 mm-N-MPa 体系，K 与 J 的单位换算是否用错。
- 平面应力/平面应变设置错误：CPE4R 为平面应力，CPE4 为平面应变，两者 J-K 换算关系不同。
- 裂纹尖端网格过粗：尖端单元尺寸应小于 a/20，且退化单元中间节点应移到 1/4 位置。
- seam 未正确指派：裂纹面未分离会导致 J 积分路径穿过材料，结果严重偏小。
- 载荷或边界条件不对称：纯 I 型问题中 KII 应接近零，若 KII 明显非零，说明模型存在不对称或约束过强。
五、可直接复用的参数与检查清单
下表汇总了本例的关键输入参数与推荐取值，可作为同类 CCT 试样仿真的起点。实际使用时只需替换几何尺寸、材料参数与载荷，其余设置基本可沿用。
| 项目 | 本例取值 | 说明 |
|---|---|---|
| 试样半宽 W | 50 mm | 裂纹长度与宽度比 α = a/W |
| 半裂纹长 a | 10 mm | α = 0.2，属于小范围屈服 |
| 板厚 t | 5 mm | 平面应力假设 |
| 弹性模量 E | 210000 MPa | 钢类材料典型值 |
| 泊松比 ν | 0.3 | — |
| 远端应力 σ | 100 MPa | 远低于屈服强度 |
| 单元类型 | CPE4R | 平面应力四边形 |
| 积分路径数 | 7 | 取第 3～5 条平均 |
| 尖端单元尺寸 | ≈ 0.5 mm | 约为 a/20 |
| KI 理论值 | 574.7 MPa·√mm | 含有限宽修正 |
| J 理论值 | 1.573 N/mm | J = K²/E |
完成本节练习后，读者应能独立完成一个 CCT 试样的轮廓积分建模，正确设置 seam、积分路径与输出请求，并依据路径无关性判断结果是否可信。下一类问题会把同样的流程迁移到含缺口或双边裂纹的试样上，届时只需调整几何与边界条件，轮廓积分的操作逻辑完全一致。