企业官网建设流程全解析

简介：面向ABAQUS复合材料失效分析的Cohesive单元建模方法文档，适合进行裂纹扩展、分层失效与断裂模拟的仿真工程师和科研人员。文档系统梳理了两类几何建模方式：整体切割薄层公用节点、独立网格通过tie绑定，并对比了各自在网格划分与赋材属性方面的适用性。随后重点讲解基于traction-separation的双线性本构模型，覆盖刚度、极限强度、临界断裂能量释放率等核心参数，以及单位统一、section实际厚度与建模厚度换算等易错细节，帮助规避界面刚度过高导致的失真结果。全包共1个doc文件，大小1.6MB，已有2229人学习，可作为ABAQUS中cohesive单元从理论到参数设置的入门参考。 开头就从我自己的经历说起吧。做仿真分析这几年，只要有同事问"两个件之间的胶层怎么模拟开裂""复合材料层合板怎么算分层"，我第一个想到的准是ABAQUS里的Cohesive单元。这东西在断裂力学和界面失效分析里太常用了，做胶接接头、涂层层间、电池极片卷绕、岩石裂纹扩展的工程问题几乎都绕不开它。今天我把这套Cohesive单元的建模思路、参数设置和踩过的坑一起整理出来，给正要上手或者已经在参数里挣扎的朋友做个参考。

这套内容的适用面很广，但主要帮助三类人：一是刚接触ABAQUS、想学界面开裂模拟的研究生；二是做结构强度校核、要评估胶接或层间失效的工程师；三是哪怕已经建过Cohesive模型，但是被不收敛、刚度穿透、单元删除显示异常这些问题折腾过，想看看别人是怎么处理的。不夸张地说，搞懂了Cohesive单元的参数逻辑，你对"仿真里怎么让材料裂开"这件事就算真正入了门。

1. 先搞明白Cohesive单元到底在算什么

1.1 它和普通实体单元的本质区别

普通实体单元是用应力-应变关系描述材料的变形，而我们常说的Cohesive单元本质上是描述"两个面之间的分离行为"——它没有实际的厚度刚度意义，而是用界面力和界面相对位移之间的关系来控制裂纹的开裂过程。

你可以把Cohesive单元想象成一层双面胶：厚度方向被拉的时候，它先给一个阻力（刚度），拉力超过某个极限（损伤起始），双面胶开始被慢慢撕开（损伤演化），直到完全失效脱离（刚度退化到零）。这种思路非常适合模拟粘接界面、复合材料层间、以及任何"裂纹路径已知或基本已知"的断裂问题。

1.2 牵引-分离法则（Traction-Separation Law）背后的直觉

ABAQUS里最常用的Cohesive单元本构是牵引-分离法则，也叫T-S法则。它描述的是界面牵引力（单位面积上的力）与相对位移（法向和切向的分离量）之间的关系。典型曲线分成两个阶段：线弹性上升段和损伤软化段。曲线下的面积就是断裂能，代表使界面完全失效所需要消耗的能量。

这个本构的优点是对参数要求相对精简：你只需要给出线弹性阶段的刚度K、起裂强度、以及断裂能。相比断裂力学里的应力强度因子，它对工程人员更友好，因为你不需要知道初始裂纹的精确位置，只要把Cohesive单元放在潜在开裂路径上，它会自己"感受"应力并决定在哪里起裂。当然这也有个前提：你得知道裂纹大概会在哪个界面上扩展，否则Cohesive单元放错了位置也白搭。

2. 三种建模方式与选型考量

2.1 基于面的Cohesive接触

第一种方式是在两个面之间定义surface-based cohesive behavior，本质上是一种含cohesive本构的接触属性，不需要在网格里插入Cohesive单元。这种方式适合模拟初始无厚度、且两个界面已经存在接触的贴合面，比如热固化成型的复合材料层板、没有胶层厚度的装配接触面。

它的优点在于建模简单，不会额外增加单元数量，也不会带来网格奇异，收敛性相对较好。但要注意，如果界面本身有真实的胶层厚度和材料属性（比如结构胶0.2mm厚），用基于面的cohesive就不够精确，因为它不包含真实的胶层体积效应，只能传递界面力。

2.2 嵌入有限厚度粘聚力层

第二种方式是直接建立有厚度的Cohesive单元层，比如把胶层几何画出来，给一层cohesive网格。这种方式最直观，参数也方便和实验数据对照——胶层的弹性模量、拉伸强度、断裂能直接对应T-S法则里的K、起裂强度和断裂能。

但这有个常见的模型刚度过高问题：如果你把Cohesive单元的实际厚度建得比胶层大很多，比如为了方便网格划分把0.1mm画成1mm，那么刚度K就会相差十倍，计算出来的开裂荷载完全不对。所以很多人在搭模型时要么用真实厚度，要么用截面属性里的constitutive thickness来修正几何厚度和材料厚度的差异。constitutive thickness是ABAQUS计算T-S法则位移时实际用的厚度，务必和真实物理厚度对应。

2.3 零厚度Cohesive层（offset节点法）

第三种是零厚度方法，也就是在相邻实体单元之间插入一层几何厚度为0的Cohesive单元。实现办法是在网格划分阶段把两个实体部分共节点区域整体偏移复制一层节点，用offset方法创建cohesive单元。这样做的好处是：你不需要建立薄层几何体，也不用担心厚胶层影响整体应力分布，适合模拟复合材料层板之间的层间界裂。

麻烦的地方在于网格生成和节点编号处理。手动操作时比较繁琐，建议用ABAQUS内置的Mesh→Edit→Insert Cohesive Cracks或者在Part模块里通过孤立网格修改来完成。另外，零厚度单元的初始穿透问题也得留意：如果两侧实体的节点位置并不完全重合，插入的cohesive单元在第一步计算时就可能出现contact穿透，导致应力异常。

为了让选型更直观，我把三种方式的适用场景做成了一个对比表，这在项目初期定方案时非常有用。

3. 关键参数设置与实操要点

3.1 刚度K的取值逻辑

Cohesive单元的初始刚度不能拍脑袋给。在T-S法则里，K = 弹性模量 / 有效厚度，但这里有个陷阱：K太小会让界面在未损伤时就产生较大相对位移，相当于人为把界面变"软"，整体结构刚度下降；K太大又容易引起数值病态，收敛变得困难。

我的实践经验是，初始刚度应尽量大，但不至于引起数值振荡。对于2D模型，一般取周围实体材料刚度E的10~100倍量级作为初始K值；对于3D问题则可以按K = E / t来估算，其中t取胶层真实厚度。如果发现收敛困难，再把K往下降一个数量级试，同时注意对比拉伸——位移曲线在起裂前的斜率是否和实验吻合。我只见过有人为了收敛把K压到和实体刚度一个量级甚至更小，结果曲线斜率明显不对，这种情况宁可调网格也不能乱牺牲刚度。

3.2 损伤起始准则：从起裂到失效

Cohesive单元的损伤起始判据有Maximum Stress、Maximum Strain、Quadratic Stress和Quadratic Strain四种常用准则。对于工程实践，我默认推荐Quadratic Stress，因为它同时考虑法向和切向应力的耦合影响，和胶接界面的失效机理更接近。Maximum Stress更偏保守，适合快速估算起裂荷载；而Quadratic Strain在实际工程中使用较少，除非你获得的实验数据是用应变形式给出的。

需要注意法向和切向强度取值不能相差太多，否则会掩盖某种失效模式。比如层合板分层问题，I型断裂能（法向拉开）和II型断裂能（切向滑移）往往差距很大，你要在损伤演化参数里单独定义。这里没有捷径，材料和工况决定一切。

3.3 粘性正则化这味"止痛药"怎么用

非线性计算不收敛的时候，很多人第一反应就是把coefficient of viscosity（粘性系数）打开，调大一点，问题确实能收敛，但这是有代价的：粘性正则化会给界面软化段增加一个"粘性阻力"，相当于人为延迟了损伤演化，峰值力可能会变得偏大，开裂位移被拉长。

我的建议是，先不加粘性正则化，优先检查网格质量和载荷步增量。如果实在要加，使用一个很小的值，比如1e-6或1e-5量级，只要能帮助收敛就可以，不要贪图省事直接上0.01。这跟用药一样，剂量要拿捏好，千万别整成常规操作。我之前处理过一个胶层剥离案例，粘性系数从1e-6改到1e-4，峰值力涨了接近10%，这个差异在工程判断里足以影响结论。

4. 网格划分与前处理技巧

4.1 Cohesive层的网格长宽比注意事项

Cohesive单元本身就是"薄层"，所以网格通常沿界面法向只有一个单元厚度，而面内尺寸取决于实体网格的大小。这导致长宽比可能做到几十甚至上百，在普通实体单元里是禁忌，但对Cohesive单元来说问题不大。不过面内尺寸不要取得太悬殊，特别是损伤扩展区域，网格尺寸越均匀，裂纹扩展路径越稳定。

在隐藏层间裂纹扩展的模型中，Cohesive层沿界面方向的网格尺寸最好和上下实体单元在界面处的尺寸对齐，否则会因节点不匹配造成应力传递不连续，最终影响损伤区域判断。遇到稳定裂纹扩展评价，建议沿裂纹扩展方向细化网格并进行网格敏感性验证，否则你无法确定裂纹扩展阻力到底是材料属性还是网格粗带来的假象。

4.2 单元删除、输出变量与后处理显示

计算完成后想看到界面开裂效果，需要设置单元删除。在Damage演化那里勾选at element deletion，允许单元在刚度退化到一定程度后被移除。同时在Field Output里输出SDEG（总刚度损伤标量）、STATUS（单元状态）和COORD（用于判断变形几何），这样后处理才能看到开裂过程。

另外要注意，输出结果默认是积分点值，Cohesive单元在厚度方向通常只用一个积分点，所以你看到的损伤云图可能显得"像一层薄壳"。这时可以开启结果可视化里的color mapping，让damage2025等外部变量参与显示，或者把SDEG切换为节点平均，预设的展示效果会好很多。操作很简单，但很多初学者都卡在"模型裂了，云图里却看不到裂纹"这一步上。

5. 常见问题与排查技巧实录

5.1 常见报错与收敛问题排查

我在实际项目里遇到最多的问题基本上集中在三类：一是单元过度变形导致负特征值；二是接触穿透或初始过盈；三是损伤穿透后单元失效导致载荷无法传递。

• 遇到负特征值，优先检查Cohesive单元是否被异常的压应力挤到畸变，同时看是不是在完全软化后还继续承载。一个常见解法是在Interaction里给失效后的实体面补一个通用接触，避免实体单元之间的穿透。
• 初始过盈，多半出现在零厚度模型里。用Surface Contact检测一下，或者干脆把零厚度单元做成极小厚度（0.001mm），能帮解决大部分初始接触异常。
• 如果加载到一半时，载荷-位移曲线突然掉得厉害，先不要急着怀疑材料参数。建议把Cohesive单元单个拿出来做单胞拉伸验证，单独确认它的起裂和断裂行为。这个"单元级标定"能省你到处试错的时间。

5.2 参数敏感性排查表

为了在排查时头脑清晰，我常把问题现象、可能原因和优先排查顺序列成一个表，这样即使隔了很久再回头看模型也知道从哪下手。

5.3 几个值得记住的实操心得

第一，Cohesive单元不适合放在单元尺寸差异过大的区域，尤其是损伤扩展路径上，网格必须过渡平滑。如果裂纹要走的方向是斜着的，Cohesive单元原则上只能在预先设置的界面方向扩展，它并不能模拟任意方向的裂纹转向，这时候要考虑XFEM或XFEM+Cohesive混合方案。很多人拿着Cohesive当万能裂缝神器用，结果计算完裂纹方向不对，那不是参数问题，是方法选择问题。

第二，如果做显式分析（Explicit），时间步长会被Cohesive单元的高刚度严重拖慢。一种常用手段是适当减小初始刚度K或给材料加小量密度缩放，但要对照动态效应是否被放大。隐式分析（Standard）虽然不怕时间步，但收敛性对参数更敏感。同一套参数在Standard里能算，拿到Explicit就得重新试，这点我希望各位提前有心理准备。

第三，当整个模型中包含多个Cohesive界面时，尽量对每个界面单独设置截面属性和单元集，方便在后续调试中单独关掉某个界面来排查问题。避免把所有cohesive混在一个set里，否则后处理能把你折磨疯，尤其当损伤区域不止一个时，哪个界面先裂甚至都分不出来。

6. 一个完整小案例：胶接接头拉伸

6.1 建模思路与关键步骤

就拿最常见的单搭接胶接接头来说，这个案例最适合用来检验Cohesive单元参数设得是否合理。模型由两块铝板和一个胶层组成，铜板用实体单元C3D8R，胶层用COH3D8（三维cohesive单元，厚度方向一层）。

操作流程是这样的：先在Part里建立两块铝板几何，中间留出0.2mm间隙；在Property模块定义一个cohesive材料，设置弹性属性Knn=Kss=Ktt，损伤起始用Quadratic Stress，损伤演化按能量类型选BK准则并给出GIc、GIIc和指数eta=1.2（BK准则指数）；然后在Assembly中装配，其中胶层用独立Part或者直接在实体网格间通过offset节点生成零厚度cohesive单元。如果需要体现真实胶层厚度，就采用有限厚度法建立0.2mm厚度的cohesive层，然后定义截面时constitutive thickness设为0.2。

边界条件这块没什么特别，一端固定，另一端给一个位移加载，每秒0.5mm，分50步加载到10mm。为了让结果更容易收敛，我把胶层两端边缘做了局部细化，让Cohesive单元长度在1mm左右保持规则。

6.2 参数和结果怎么看

给出一个推荐初始参数集（单位制统一为mm、N、MPa体系）：Knn=Kss=Ktt=1e6 N/mm³，法向强度3MPa，切向强度5MPa，GIc=0.4 N/mm，GIIc=0.8 N/mm，BK指数设为1.2。这一组数偏中性，不同胶种可以在此基础上调整。

计算完以后，先看载荷-位移曲线：起裂前是明显的线性段，曲线斜率对应界面刚度；第一个峰值就是起裂荷载，随后进入软化段，如果曲线出现断崖式下跌说明断裂能给的偏低或网格太粗；如果曲线缓慢下降，说明断裂能相对偏高。我会再用后处理SDEG确认损伤的起始位置：应该从搭接区端部开始，沿胶层逐渐向中心推进，这个趋势和实际胶接接头失效模式一致。如果损伤先出现在中间，那基本就是参数或者边界条件有问题。

我可以很明确地讲，照着上面这套流程做一遍，Cohesive单元的基础应用你基本就能拿下。后面再遇到复合材料层合板低速冲击、岩石裂纹扩展、电池卷绕极组分层这类问题，思考路径都是一样的：找到潜在开裂界面，选对建模方式，标定好那三个核心参数（K、强度、断裂能），剩下的就是耐心调试网格和收敛性。

最后再分享一个我自己的小习惯：每次新建Cohesive模型时，都会先跑一个单胞拉伸的简单模型，验证参数在单个单元上的响应符合预期，再放进完整结构里计算。这看起来多花十分钟，但能帮你省下后面数不清的排错时间。毕竟仿真里"算得出来"和"算得对"是两回事，Cohesive单元这种对参数灵敏得不得了的模型，一定要逼自己养成先单元级标定、再整机计算的好习惯。

本文还有配套的精品资源，点击获取