2. 粘接界面力学基础

在Abaqus中模拟粘接结构，核心难点并不在胶层几何本身，而在于如何用连续介质单元或界面单元去“等效”一层极薄、却承担载荷传递与破坏的胶层。内聚力模型（Cohesive Zone Model, CZM）正是把胶层的损伤与断裂压缩到一条牵引—分离（Traction–Separation）本构关系上的方法。理解这条曲线上的每一个参数，是后续建模能否收敛、结果是否可信的前提。

一、内聚力模型：把断裂过程压缩成一条本构曲线

CZM的基本假设是：界面在完全开裂前仍能传递应力，应力大小取决于上下表面的相对分离量。以最常见的双线性本构为例，横轴为分离量δ，纵轴为牵引力t。曲线分为三段：上升段（δ从0到δ⁰）为线弹性，斜率即界面刚度K；峰值点（δ⁰, t⁰）为损伤起始点，t⁰即界面强度；下降段（δ⁰到δᶠ）为损伤演化段，刚度随损伤变量D从0增至1而退化；δᶠ处牵引力归零，界面完全脱粘。曲线下的面积就是临界能量释放率G_c。

需要强调：CZM不是真实胶层的微观物理，而是一种“等效”本构。它把胶层厚度、塑性、微裂纹萌生与聚合等复杂过程，折算为界面刚度、强度和断裂能三个可标定的宏观参数。因此参数取值必须与胶层实际厚度、加载模式相匹配，不能随意套用文献值。

二、断裂力学参数：G_c、t⁰与K的物理含义

临界能量释放率G_c是CZM中最具物理意义的参数，它表示形成单位面积新裂纹面所需的能量，单位常用N/mm（即mJ/mm²）。G_c可由断裂力学试验获得：Mode I用双悬臂梁（DCB）试验，Mode II用端部缺口弯曲（ENF）试验，混合模式用MMB试验。对于同一种胶，G_IC通常明显小于G_IIC，说明界面抗剥离能力弱于抗剪切能力。

界面强度t⁰（Mode I为t_n⁰，Mode II为t_s⁰）是损伤起始的应力门槛。它并不等于胶的拉伸强度，而是界面在特定模式下的等效强度。界面刚度K则是牵引—分离曲线的初始斜率，单位MPa/mm。K越大，界面越“硬”，初始响应越接近刚性连接；K过小则界面过软，会引入非物理的柔度。

三、损伤起始与演化准则：从何时开始坏，到怎样坏完

损伤起始准则回答“什么时候开始损伤”。Abaqus中常用二次应力准则（QUADS）或二次应变准则（QUADSCRT），其形式为：

(⟨t_n⟩/t_n⁰)² + (t_s/t_s⁰)² + (t_t/t_t⁰)² = 1

其中⟨t_n⟩表示只考虑拉伸正应力，压缩不贡献损伤。当左侧达到1时，D从0开始增长。损伤演化准则回答“损伤如何发展”。最常用的是基于能量的线性或指数演化，即用断裂能G_c控制软化段形状。以线性演化为例，等效分离量δᶠ由G_c和t⁰共同决定：

δᶠ = 2G_c / t⁰   （双线性、纯Mode I）

损伤变量D在演化段按D = (δ - δ⁰)/(δᶠ - δ⁰)线性增长，单元刚度按(1-D)退化。当D=1时，界面完全失效，牵引力为零。混合模式下还需引入Benzeggagh–Kenane（BK）准则描述G_IC与G_IIC的混合关系：

G_c = G_IC + (G_IIC - G_IC) · (G_II/G_T)^η

其中η为BK指数，通常取1.5～2.5，需由混合模式试验标定。

四、参数标定与建模中的实际取值

以一个典型结构胶为例：胶层厚度0.2 mm，弹性模量E=2000 MPa，Mode I断裂能G_IC=0.5 N/mm，界面强度t_n⁰=20 MPa。按双线性模型计算：

K = E/t = 2000/0.2 = 10000 MPa/mm
δ⁰ = t⁰/K = 20/10000 = 0.002 mm
δᶠ = 2G_IC/t⁰ = 2×0.5/20 = 0.05 mm

可见δ⁰远小于δᶠ，软化段很长，这正是胶层“脆而韧”的典型特征。在Abaqus中，这些参数分别填入Cohesive Section的Elastic（K）、Damage Initiation（t⁰）和Damage Evolution（G_c）中。若采用牵引—分离本构，需选择“Traction Separation”而非“Continuum”；若用连续单元模拟胶层，则需在材料中定义Damage for Traction Separation Laws。

五、网格与粘性正则化的影响

CZM对网格尺寸敏感。软化段内单元尺寸过大，会导致损伤局部化无法分辨，载荷—位移曲线出现虚假的“跳变”。经验上，沿预期裂纹扩展路径的单元尺寸应小于δᶠ的1/3～1/5。以上例δᶠ=0.05 mm为例，单元尺寸宜取0.01 mm左右。若模型规模不允许，可引入粘性正则化（Viscosity）系数，典型取1e-5～1e-3，帮助收敛，但过大会人为提高断裂能，需做敏感性检查。

此外，损伤起始准则中的二次应力准则适用于混合模式，若只关心Mode I，可简化为最大正应力准则。演化准则中，线性与指数软化对最终G_c无影响，但会影响收敛性和峰值后响应，指数软化通常更易收敛。

理解上述参数的物理意义后，才能在Abaqus中做出有依据的取舍：K决定初始柔度，t⁰决定何时开始损伤，G_c决定最终承载能力，而δᶠ与网格共同决定计算能否稳定走完软化段。这四者相互制约，任何一项脱离实际胶层行为，都会让脱粘仿真失去预测价值。