4. Abaqus断裂仿真建模准备
断裂仿真的建模准备与常规静力分析有本质区别：裂纹尖端附近存在高应变梯度，材料本构需要包含损伤起始与演化，输出请求必须能捕捉损伤变量与能量释放过程。本节以“带预制裂纹的CT试样（紧凑拉伸试样）在拉伸载荷下的裂纹扩展”为教学示例，完整走通从部件到历史输出的建模准备链路。示例中的几何尺寸、材料参数均为合理假设值，用于说明操作逻辑，不代表任何实测结果。
一、部件与裂纹几何的建立
断裂仿真中，裂纹本身不是“画出来”的，而是通过几何分区或网格不连续来定义。对于CT试样，典型做法是：在Part模块中先建立二维平面（或三维）试样轮廓，再在裂纹预期扩展路径上切分（Partition）出一条窄带或一条线，作为后续指派裂纹面或seam的几何基础。示例几何：宽度W=50 mm，厚度B=12.5 mm，初始裂纹长度a₀=20 mm（即a₀/W=0.4），加载孔中心距裂纹面12.5 mm。裂纹尖端位置必须精确落在几何点上，否则后续轮廓积分或XFEM富集区会偏移。
若采用XFEM（扩展有限元），裂纹面通过“Special → Crack → Assign seam”指派，裂纹尖端通过“Crack → Edit → Specify crack front”定义；若采用VCCT或内聚单元，则需在裂纹路径上布置一层零厚度或极薄的内聚单元层。两种路线对几何准备的要求不同：XFEM允许裂纹穿过单元内部，网格无需与裂纹面对齐；内聚单元法则要求网格边界与裂纹路径严格一致。选择哪条路线，应在Part阶段就确定，因为它直接决定后续网格划分策略。
二、材料本构：损伤起始与演化
断裂仿真的材料定义至少包含三部分：弹性、塑性（可选）、损伤。以铝合金为例，弹性模量E=71000 MPa，泊松比ν=0.33。损伤采用Abaqus内置的Ductile Damage（延性损伤）模型：损伤起始由等效塑性应变阈值控制，示例取ε₀=0.05；损伤演化采用基于能量的线性软化，断裂能Gf=10 N/mm。在Material → Damage for Ductile Metals → Ductile Damage中，需同时设置“Damage Initiation”和“Damage Evolution”两个子选项，缺一不可——只设起始不设演化，材料会在达到阈值后瞬间失去承载能力，导致收敛困难且结果失真。
若采用XFEM，还需在材料中激活“Damage Evolution”下的“Maximum Principal Strain”或“Energy”准则，并在Section中勾选“XFEM”选项。注意：XFEM与Ductile Damage可以共存，但损伤演化参数需与XFEM的富集区行为协调，否则会出现裂纹不扩展或扩展路径异常。
三、截面指派与装配
截面类型取决于几何维度：二维平面应力问题选“Solid, Homogeneous”，并指定厚度（示例B=12.5 mm）；三维问题选“Solid, Homogeneous”并指派给实体区域。截面指派后，务必在Assembly模块中检查实例是否独立（Independent）或依赖（Dependent）。对于断裂仿真，若需在装配后对裂纹路径进行局部网格细化，建议使用Independent实例，以便在Mesh模块中对实例直接划分而不影响原始Part。
装配阶段还需定义加载孔与夹具的耦合关系。示例中，上加载孔通过“Coupling”约束耦合到参考点RP-1，下加载孔耦合到RP-2，载荷施加在参考点上。这样做的目的是避免孔边应力集中干扰裂纹尖端场，同时便于施加位移或力控制。
四、分析步与求解控制
断裂仿真通常采用“Static, General”分析步，但必须开启几何非线性（Nlgeom=On），因为裂纹扩展伴随大位移和单元畸变。增量步设置是关键：初始增量建议取0.01，最小增量1E-6，最大增量0.1。若采用损伤演化，需在“Other → General Solution Controls”中适当放宽收敛容差，或改用“Automatic Stabilization”添加少量阻尼，以帮助损伤软化阶段的收敛。
对于裂纹扩展这类强非线性问题，建议分两个分析步：第一步施加预载荷（如0.1 mm位移），让接触和耦合关系稳定；第二步施加主载荷（如1.0 mm位移），触发损伤起始与演化。这样可以在第一步中检查模型是否正常，避免直接进入损伤阶段后因其他问题导致不收敛而难以排查。
五、输出请求与历史输出设置
断裂仿真的输出设置直接决定能否提取裂纹扩展阻力曲线（R曲线）和损伤演化数据。场输出（Field Output）中必须勾选：S（应力）、E（应变）、PE（塑性应变）、DAMAGET（损伤起始准则值）、DAMAGEC（损伤演化变量）、STATUS（单元状态，用于XFEM或内聚单元）、U（位移）。输出频率建议每5个增量步一次，避免文件过大。
历史输出（History Output）是断裂仿真的核心。需在裂纹尖端附近定义历史输出区域（Set），输出：裂纹尖端开口位移CTOD、J积分（通过Contour Integral定义）、裂纹扩展长度（通过XFEM的Crack Length或内聚单元的SDEG）。示例中，在裂纹尖端前方沿路径定义一组单元集，输出其DAMAGEC和STATUS，用于判断裂纹何时到达该位置。此外，参考点RP-1的位移和反力必须输出，用于绘制载荷-位移曲线。
六、建模准备的检查清单
在提交计算前，建议按以下顺序逐项核对，避免因准备疏漏导致计算中断或结果不可用：
- 几何：裂纹尖端位置是否精确落在几何点上；XFEM的seam和crack front是否指派正确；内聚单元层是否与裂纹路径对齐。
- 材料：损伤起始与演化是否同时定义；断裂能单位是否与几何单位制一致（示例为N/mm）；XFEM选项是否在Section中激活。
- 截面：厚度是否与几何维度匹配；Independent/Dependent实例选择是否合理。
- 装配：加载孔耦合参考点是否建立；边界条件是否约束了刚体位移。
- 分析步：Nlgeom是否开启；增量步范围是否合理；是否分步加载。
- 输出：DAMAGET、DAMAGEC、STATUS是否勾选；历史输出区域是否覆盖裂纹路径；参考点反力是否输出。
完成上述准备后，模型即可提交计算。断裂仿真的建模准备不是一次性动作，而是一个需要根据试算结果反复调整的过程——尤其是损伤演化参数和网格密度，往往需要2~3轮试算才能获得稳定的裂纹扩展响应。建议先用较粗网格和较大断裂能试算，确认裂纹能正常扩展后，再逐步细化网格并调整参数至目标精度。