Cohesive讲解.doc

复合材料模型建模与分析1．Cohesive 单元建模方法1.1 几何模型使用内聚力模型( cohesivezone )模拟裂纹的产生和扩展，需要在预计产生裂纹的区域加入 cohesive 层建立 cohesive 层的方法主要有：方法一、建立完整的结构(如图1(a)所示)，然后在上面切割出一个薄层来模拟 cohesive单元，用这种方法建立的 cohesive 单元与其他单元公用节点，并以此传递力和位移 方法二、分别建立 cohesive 层和其他结构部件的实体模型，通过“ tie ”绑定约束，使得cohesive单元两侧的单元位移和应力协调，如图 1 (b)所示a)cohesive单元与其他单元公用节点(b)独立的网格通过“tie ”绑定图 1. 建模方法上述两种方法都可以用来模拟复合材料的分层失效，第一种方法划分网格比较复杂；第二种方法赋材料属性简单，划分网格也方便，但是装配及“ tie ”很繁琐；因此在实际建模中我们应根据实际结构选取较简单的方法1.2 材料属性应用cohesive单元模拟复合材料失效，包括两种模型：一种是基于 traction-separation描述；另一种是基于连续体描述。

其中基于 tractio n-separatio n 描述的方法应用更加广泛而在基于traction-separation 描述的方法中，最常用的本构模型为图 2所示的双线性本构模型它给出了材料达到强度极限前的线弹性段和材料达到强度极限后的刚度线性 降低软化阶段注意图中纵坐标为应力，而横坐标为位移，因此线弹性段的斜率代表的实 际是cohesive单元的刚度曲线下的面积即为材料断裂时的能量释放率因此在定义 cohesive的力学性能时，实际就是要确定上述本构模型的具体形状：包括刚度、极限强度、 以及临界断裂能量释放率，或者最终失效时单元的位移常用的定义方法是给定上述参数 中的前三项，也就确定了 cohesive的本构模型Cohesive单元可理解为一种准二维单元， 可以将它看作被一个厚度隔开的两个面，这两个面分别和其他实体单元连接 Cohesive单元只考虑面外的力，包括法向的正应力以及 XZ, YZ两个方向的剪应力下文对cohesive单元的参数进行阐述，并介绍参数的选择方法图2.双线性本构模型1.2.1cohesive 单元的刚度基于traction-separation 模型的界面单元的刚度可以通过一个简单杆的变形公式来理解PLAE(1)其中L为杆长，E为弹性刚度，A为初始截面积，P为载荷。

公式（1）又可以写成其中S P A为名义应力，K EL为材料的刚度为了更好的理解K,我们把K E L写成:k E E!空（3）L 1 L这里我们用L来代替1，其中L可以理解为建模厚度，即建模时 cohesive in terface 的几 何厚度；L为实际厚度，即cohesiveinterface 的真实厚度，这个厚度在cohesivesection中定义E L可以理解为几何刚度，即模型中 cohesive in terface 所具有的刚度；乎为cohesive in terface 的真实刚度当L为1时，计算界面刚度就采用几何刚度 E； L，当L为0.001时，计算时界面刚度变为1000E L o举个小例子，如果界面的实际厚度为 0.01，而在建模时就是按照这个厚度建立的，在定义material-section 时又specify这层的厚度为 0.01，实际上就等于把界面刚度提高了 2个数量级，模拟结果当然是不对的，这时定义section时应采用默认厚度1ABAQU在cohesive建模中使用了很“人性化”的设计，实际问题中界面可能很薄，有的只有0.001mm 甚至更小有些问题 cohesive单元的in terface 还可能是0厚度（比如 crack问题），而相对来说整体模型也许很大，如果不引入这两个厚度，我们就要在很大 的模型中去创建这个很小的界面这是一个很麻烦的事情。

引入这两个厚度，在建模时我们 就可以用有限的厚度来代替这个很小的界面厚度，只要在 section中定义这个L就好了注：以上大部分内容来自仿真论坛：再议 cohesive应用中对于一些参数的理解）F面举例来说明cohesive单元刚度的设置过程，以 ABAQUS6.为例:进入property 界面，点击Material — Creat ,在弹出的EditMaterial 对话框中，可以 编辑新创建的cohesive材料的名称，然后点击 Mechanical — Elasticity — Elastic —Traction ，在空格中输入相应的刚度图 3.cohesive 单元刚度的定义1.2.2 损伤准则1.2.2.1 初始损伤准则 初始损伤对应于材料开始退化，当应力或应变满足于定义的初始临界损伤准则，则此时退化开始Abaqus的 Damagefortractionseparationlaws 中包括： QuadeDamageMaxeDamage QuadsDamage MaxsDamage MaxpeDamage MaxpsDamag六种初始损伤准贝U, 其中前四种用于一般复合材料分层模拟，后两种主要是在扩展有限元法模拟不连续体（比 如 crack 问题）问题时使用。

使用图2所示的双线本构模型，其中：t0、t0及tt0分别代表纯I型、纯U型或纯川破坏 的最大名义应力， n0、 s0， t0 代表相应的最大名义应变，当定义界面单元的初始厚度为 1时，则名义应变等于与之相对应的相对位移 n , s及tQuadeDamag为二次名义应变准则：当名义应变比的平方和等于 1时，损伤开始MaxeDamag为最大名义应变准则：当任何一个名义应变的比值达到 1时，损伤开始QuadsDamag为二次名义应力准则：当各个方向的名义应变比的平方和等于 1时，损伤开始MaxsDamag为最大名义应力准则：当任何一个名义应力比值达到 1时，损伤开始图 4. 初始损伤准则定义EditMaterial 对话框中，点击 Mecha ni cal — DamageforTract ion Separati on Laws ，然后根据自己的需要点击相应的损伤准则其中最常用是 QuadsDamage1.2.2.2 损伤演化规律选择了初始损伤准则之后，然后点击Suboptions — DamageEvolution，窗口如图5所示 其中Type包括Displacement和Energy，Displacement为基于位移的损伤演化规律，而 Energy 为基于能量的损伤演化规律。

Softening 中包括 Linear ， Exponential 及 Tabular 三种刚度退化方式 DamageEvolution中的所有的选项都是用来确定单元达到强度极限以后的刚度降阶方式一般常用：以能量来控制单元的退化，即 Type— En ergy ;线性软化模型，即 Softening —Linear ， Degradation —Maximum；Mixedmodebehavior—BK， Modemixratio — Energy, 并选中 Power图 5. 损伤演化规律定义1.3cohesive 单元界面属性还是在Property 界面中，点击 Section — Create，在弹出的EditSection 对话框中， 选择 Other —Cohesive图 6. 定义材料的界面属性在 EditSection 对话框中，在 material 的下拉菜单中选择刚才创建的 cohesive 材料， 也可以点击右侧的 create 创建一组新的材料； Response 选择 tractionseparation Initialthickness 为前文提到的 L ，默认值为 1，也可以在 specify 中指定一个特定的值。

1.4 将所创建的界面属性赋予几何实体点击Assign — Section，然后在视图中选中要赋的几何实体，点击左下角的 Done,则弹出如下窗口，在窗口是Section中下拉选中所创建的Cohesive截面，点击0K操作完成图 7. 给实体赋截面属性1.5cohesive 单元网格划分Cohesive 单元网格的划分与其他单元基本一致，但是以下几点不同与其他单元，划分网格 时应特别注意一、 网格密度， cohesive 单元的网格尺寸不能太大，通常需要比较精细的网格，不然 容易引起收敛性问题，甚至无法继续计算二、 必须使用sweep （扫掠）戈扮网格的方法，并且扫掠的方向垂直于 cohesive面，即沿着 cohesive 单元的厚度方向单元种类的选择图 8.cohesive 单元种类选择在单元库中选择 cohesive ，可以在 Viscosity ，specify 中指定一粘性系数，来改善收敛性，但是粘性系数的设置不能太大，不然会影响计算结果，我们一般设置为 0.001 ；Eleme ntdeletion :用于设置单元的删除情况，一般选 yes，即当单元完全失效时被删除； maxdegradation :一般设置为1，即当SDEG=时，认为单元失效。

2.cohesive 单元在复合材料分层分析中的应用为了验证商用有限元软件 ABAQU中的cohesive单元在复合材料分层计算时的有效性， 我们通过其与一实验值的对比验证了其计算的准确性一 DCB试验件，长150mm宽20mm单臂厚度1.98mm预置55mn长的初始裂纹，如图 9 所示材料属性为 E11=150GPa， E22=E33=11GPa， G12=G13=6.0GPa， G23=3.7GPa， 12=0.25，13 =0.25， 23 =0.45 ; cohesive 单元的材料属性为 K=1X 105MPa/mm 界面强度 T=15MPa临界能量释放率 G=0.268KJ/m2悬臂梁一端固支，一端施加位移载荷a) 侧视图(b) 俯视图图9.DCB几何模型Abaqus和实验［1］得到的力位移曲线如图10所示，从图中可以看出，数值模拟的力位 移曲线与实验得到的力位移曲线吻合的很好，数值模拟得到的最大力为 65.8N，而实验得到的最大力为62.52N，数值模拟结果略高于实验结果由此，我们可以得到有限元软件 ABAQU中的cohesive单元可以有效的模拟复合材料层合板的分层计算得到的变形过程的应力及位移云图如图 11、12所示a) abaqus 计算值(b)实验值图 10. 实验及数值模拟结果( a) t=0.25s ( b) t=0.5s( c ) t=0.75s ( d ) t=1.0s图 11. 变形过程中应力云图( a) t=0.25s ( b) t=0.5s( c ) t=0.75s ( d ) t=1.0s图 12. 变形过程中位移云图[1]Camanho,P.P.,Davila,C.G.,2002.Mixed-modedecohesionfiniteelementsforthesimulationofdelaminationincompositematerials.NASA/TM-2002-211737.。