3. Abaqus粘接建模技术路线
在Abaqus中模拟粘接界面的脱粘,可选的技术路线不止一条。内聚力单元(Cohesive Element)、内聚力接触(Cohesive Surface/Contact)、虚拟裂纹闭合技术(VCCT)和扩展有限元法(XFEM)是工程与科研中最常用的四种。它们并非互相替代的“四个按钮”,而是分别对应不同的物理假设、损伤描述能力和网格策略。选错路线,往往不是算不出来,而是算出来的结果在物理上站不住脚。
四种方法的核心差异
内聚力单元在粘接界面处插入一层零厚度或极薄厚度的实体单元,通过牵引-分离本构(Traction-Separation Law)描述界面从完好到损伤再到完全失效的过程。它天然支持损伤起始、演化与最终脱粘,能输出界面刚度退化、损伤变量SDEG等状态量,适合需要观察脱粘萌生与扩展全过程的场景。代价是必须在几何上为界面预留单元层,网格划分和界面位置需要提前确定。
内聚力接触把内聚力本构定义在接触对(Surface-to-Surface)上,不需要为界面单独建单元。它保留了内聚力模型的损伤演化能力,同时避免了界面单元带来的网格约束,特别适合界面位置明确但几何复杂、难以插入单元层的模型,例如曲面粘接、多材料装配体。其本质仍是“界面有强度、有断裂能”,只是承载本构的载体从单元变成了接触面。
VCCT基于断裂力学能量释放率准则,假设裂纹始终沿预先定义的界面路径扩展。它不模拟损伤萌生,而是直接判断裂纹尖端能量释放率是否达到临界值。因此VCCT适合界面已经存在初始裂纹、需要计算裂纹扩展驱动力和扩展速率的场景,例如分层、脱粘扩展分析。它计算效率高,但对初始裂纹和扩展路径有强依赖。
XFEM允许裂纹在单元内部独立于网格扩展,不需要预先定义裂纹路径。用于粘接脱粘时,它适合裂纹路径不确定、可能偏离界面或需要在连续体中自然萌生裂纹的场景。但XFEM对界面粘接的“界面属性”描述不如内聚力模型直接,通常需要结合损伤准则或接触属性使用,且计算收敛性和结果解释难度较高。
选择依据:从问题出发,而不是从方法出发
选择技术路线的第一步,是明确你要回答的问题。如果问题是“粘接接头在拉伸载荷下何时开始脱粘、脱粘如何扩展”,那么需要的是能描述损伤萌生与演化的模型,内聚力单元或内聚力接触是首选。如果问题是“已有分层在给定载荷下的能量释放率是多少、裂纹会不会继续扩展”,VCCT更直接。如果问题是“裂纹可能从粘接界面附近任意位置萌生并偏离界面”,XFEM的网格无关性更有优势。
第二步是检查几何与网格条件。界面位置是否明确?能否在几何上插入一层单元?如果界面是复杂曲面或装配体中的接触面,插入单元层代价很高,内聚力接触更现实。如果界面路径完全未知,VCCT的预设路径假设就不成立,XFEM或内聚力接触配合损伤准则更合适。
第三步是评估计算成本与收敛性。内聚力单元和接触都需要定义损伤起始准则(如二次应力准则QUADS)和损伤演化(基于能量或位移),参数包括界面强度、断裂能、损伤演化类型。VCCT不需要损伤演化参数,但需要临界能量释放率GIC、GIIC。XFEM需要定义损伤起始与演化,且富集单元带来的收敛问题更突出。在工程实践中,如果界面参数已知且界面位置明确,内聚力单元往往是最稳妥的起点;如果界面复杂,内聚力接触是更灵活的选择。
一个可参考的选择决策表
| 判断条件 | 内聚力单元 | 内聚力接触 | VCCT | XFEM |
|---|---|---|---|---|
| 界面位置是否明确 | 必须明确 | 必须明确 | 必须明确 | 可不明确 |
| 是否需要损伤萌生 | 支持 | 支持 | 不支持 | 支持 |
| 是否需要预设裂纹路径 | 是(界面层) | 是(接触面) | 是 | 否 |
| 是否适合复杂曲面 | 较难 | 适合 | 较难 | 适合 |
| 关键输入参数 | 界面强度、断裂能 | 界面强度、断裂能 | GIC、GIIC | 损伤准则、断裂能 |
| 典型输出 | SDEG、脱粘长度 | CSDMG、脱粘长度 | G_I、G_II | PHILSM、裂纹路径 |
参数设置中的关键差异
内聚力单元和接触都需要在材料属性中定义牵引-分离本构。以双线性本构为例,需要输入界面刚度K、最大牵引力(界面强度)和断裂能Gc。损伤起始准则常用QUADS(二次应力准则),损伤演化常用基于能量的线性或指数软化。VCCT则需要在接触属性或裂纹面定义中指定Bonded接触和VCCT准则,输入GIC、GIIC以及混合模式准则(如BK准则)。XFEM需要在材料中定义损伤起始(如最大主应力准则)和损伤演化(基于能量),并在相互作用中指定XFEM裂纹域。
一个常见的误区是:把内聚力单元的界面刚度设得过高,导致收敛困难;或者把断裂能设得过低,导致脱粘过早发生。合理的做法是先用文献或实验数据确定界面强度和断裂能,再根据网格尺寸调整界面刚度,使牵引-分离曲线在数值上稳定。对于内聚力接触,还需要控制接触稳定化参数,避免刚接触时出现刚体位移。
最终,技术路线的选择不是非此即彼。在复杂模型中,可以组合使用:例如在主要粘接界面使用内聚力接触,在局部关键区域插入内聚力单元,或者用VCCT判断已有分层的扩展趋势,再用内聚力模型模拟新脱粘的萌生。理解每种方法的物理假设和数值特性,才能让仿真结果真正服务于工程判断。