4. 内聚力单元建模实战

内聚力（Cohesive，简称COH）单元是Abaqus中模拟粘接界面脱粘的核心工具。它的本质是一层“零厚度或薄厚度”的界面单元，通过牵引-分离（Traction-Separation）本构描述界面从完好到损伤再到完全脱粘的过程。本节从零开始，完整走一遍COH单元建模的关键操作链：创建单元、定义几何厚度、设置本构、指派材料与截面、划分网格。每一步都直接影响后续能否收敛、能否得到合理的脱粘响应。

一、创建COH单元：两种几何来源

在Abaqus/Standard中，COH单元通常使用COH2D4（二维四节点）或COH3D8（三维八节点）等单元类型。创建方式主要有两种：

• 独立建层法：在Part模块中单独绘制一层几何（二维为一条线，三维为一个面），划分网格后通过Element Type指定为COH单元。这种方式几何清晰，适合界面位置明确、需要单独控制网格的情形。
• 共享节点法：不单独建层，而是让上下两个部件在界面处共享节点，再通过Mesh Edit → Create Cohesive Elements在已有网格基础上插入COH单元。这种方式不需要额外几何，但要求上下部件在界面处网格完全匹配。

实际工程中，若粘接层本身有物理厚度（如结构胶层0.2 mm），推荐独立建层；若模拟的是理想界面（如复合材料层间），共享节点法更简洁。无论哪种方式，COH单元的节点顺序必须正确——二维单元沿厚度方向的两个节点对，三维单元沿法向的两个面，节点编号顺序错误会导致法向反向，牵引力方向随之出错。

二、几何厚度：本构厚度 vs 物理厚度

这是COH单元最容易出错的地方。Abaqus中COH单元有两个“厚度”概念：

• 几何厚度（Geometric thickness）：Part中实际绘制的厚度。若独立建层时画了0.2 mm的胶层，几何厚度就是0.2 mm；若画的是零厚度线/面，几何厚度为0。
• 本构厚度（Constitutive thickness）：在Section定义中通过Initial thickness指定，默认值为1.0。它用于将名义应变转换为分离量：δ = ε × T₀。

关键原则：牵引-分离本构中的分离量δ和刚度K，都是基于本构厚度定义的。若几何厚度为0（零厚度界面），必须将Initial thickness设为1.0，此时δ直接等于位移跳跃；若几何厚度为0.2 mm，则应将Initial thickness设为0.2，使δ = ε × 0.2，与物理胶层一致。若两者不匹配，会导致刚度、损伤起始位移和失效应变全部错位。

三、本构定义：双线性牵引-分离

最常用的本构是双线性（Bilinear）牵引-分离模型，需要输入三个核心参数：

• 弹性刚度K：初始线弹性斜率，单位MPa/mm。对于结构胶，K通常取10³~10⁵ MPa/mm量级；K过大导致收敛困难，过小则界面“过软”。
• 损伤起始强度t₀：法向或切向的峰值牵引力，单位MPa。达到该值即开始损伤。
• 断裂能G_C：单位面积脱粘所需能量，单位N/mm（=mJ/mm²）。由G_C和t₀可反算失效应变：δ_f = 2G_C/t₀。

在Abaqus中，通过Material → Mechanical → Damage for Traction Separation Laws → Quads Damage定义损伤起始，通过Damage Evolution → Energy定义损伤演化。若需要混合模式，还需指定模式混合准则（如BK准则）及法向/切向强度比。

*Material, name=Adhesive
*Elastic, type=TRACTION
 1000., 1000., 1000.
*Damage Initiation, criterion=QUADS
 30., 20., 20.
*Damage Evolution, type=ENERGY, mixed mode behavior=BK, power=2.0
 0.5, 0.8, 0.8

上述示例中，法向强度30 MPa，切向20 MPa，法向断裂能0.5 N/mm，切向0.8 N/mm，BK指数2.0。这些数值属于教学示例，实际应通过试验或文献标定。

四、指派材料与截面

创建Section时选择Cohesive类型，指定材料名称和Initial thickness。然后将Section指派给COH单元区域。注意：COH单元的截面方向（法向）由单元节点顺序决定，而非由Section指定。若发现法向反了，需回到Mesh模块用Element Type或Edit Mesh翻转单元法向。

五、网格划分要点

COH单元的网格质量直接决定收敛性：

• 沿厚度方向只划分一层单元。COH单元本身代表界面，多层划分没有物理意义且会引入额外刚度。
• 面内网格尺寸应足够细以捕捉损伤前沿，通常与相邻实体单元尺寸协调，建议界面附近实体单元尺寸不超过COH单元长度的2~3倍。
• 单元形状尽量规则，避免大长宽比。二维COH2D4建议长宽比不超过5:1。
• 网格匹配：若采用共享节点法，上下部件在界面处的网格必须完全一致，否则无法插入COH单元。

六、建模检查清单

完成上述步骤后，建议按以下清单逐项核对，避免常见错误：