Abaqus中Cohesive单元建模：从本构到失效的完整指南

1. 为什么我们需要Cohesive单元？从“胶水”失效说起

大家好，我是老张，在CAE仿真这个行当里摸爬滚打了十几年，复合材料、胶接结构的失效分析做了没有一百也有八十个。今天想和大家聊聊Abaqus里一个既让人爱又让人恨的家伙——Cohesive单元。很多朋友一听到“内聚力模型”、“损伤演化”这些词就头大，觉得这是高深的理论，离实际应用很远。其实不然，你可以把它想象成结构里的“胶水”。我们做设计的时候，最怕的就是这个“胶水”在关键时刻掉链子，比如飞机机翼的蒙皮和桁条脱胶了，或者手机屏幕的玻璃和边框开胶了。Cohesive单元，就是我们在电脑里模拟这种“胶水”粘接和失效过程的最有力工具。

在Abaqus里模拟界面行为，主要有两种思路：一种是定义接触属性，另一种就是插入Cohesive单元。接触属性更偏向于描述两个已经分离的表面之间的相互作用，比如摩擦、硬接触。而Cohesive单元则不同，它被设计用来模拟那些原本粘结在一起，但在载荷下会逐渐产生损伤、最终完全分离的界面。它就像是在两个实体之间预先铺好的一层“胶水”单元，这层单元有自己的力学行为规则。为什么我更喜欢用Cohesive单元呢？因为它能更物理地描述从完好到损伤起始，再到损伤扩展直至完全失效的全过程，这对于评估结构的损伤容限和剩余强度至关重要。

原始文章提到了Cohesive单元通常被假设为零厚度，这是一个非常重要的概念，也是新手容易困惑的地方。你可能会问，一个没有厚度的单元，它的应变、位移怎么定义？其实，这是一种数学和力学上的简化。我们把界面的力学响应（应力与相对位移的关系）凝聚在这个“几何厚度为零”的面上。它的方向性也很关键，通常定义一个法向（垂直于界面）和两个切向（在界面平面内）。所有的本构关系都建立在这个局部坐标系下。理解了它是个“零厚度的胶水层”，后面的刚度、强度设置就容易入手了。

2. Cohesive单元的核心：理解它的“脾气”（本构与损伤）

想把Cohesive单元用对、用好，你必须摸清它的“脾气”，也就是它的力学行为规律。这个规律主要由三部分组成：弹性本构、损伤起始准则和损伤演化定律。这三部分环环相扣，共同决定了你这层“胶水”在受力时如何变形、何时开始“开裂”、以及“开裂”后如何彻底失效。

2.1 弹性本构：胶水的“初始刚度”

在“胶水”还没开裂之前，它处于纯弹性阶段。这个阶段，应力和位移（或应变）是成正比的，比例系数就是刚度。这是你需要输入的第一个关键参数。对于Cohesive单元，刚度通常按方向定义：一个法向刚度（Kn）和两个剪切刚度（Ks, Kt）。在Abaqus的材料模块中，你会在“Cohesive”或“Traction”相关的选项中找到它们。

这里有个非常关键的坑，我踩过不止一次：刚度值不能乱设。理论上，这个刚度应该足够大，以确保在界面失效前，Cohesive单元本身的弹性变形远小于相邻实体材料的变形，从而不会不适当地“软化”整体结构。但也不能设得太大，否则会导致严重的数值收敛问题。一个经验法则是，让Cohesive单元的刚度与相邻实体材料的等效模量除以一个