ABAQUS中Cohesive单元建模方法

1、复合材料模型建模与分析1.CohesiV元建模方法几何模型使用内聚力模型(cohesivezone)模拟裂纹的产生和扩展，需要在预计产生裂纹的区域加入cohesive层。建立cohesive层的方法主要有：方法一、建立完整的结构(如图1(a)所示)，然后在上面切割出一个薄层来模拟cohesive单元，用这种方法建立的cohesive单元与其他单元公用节点，并以此传递力和位移。方法二、分别建立cohesive层和其他结构部件的实体模型，通过“tie"绑定约束，使得cohesive单元两侧的单元位移和应力协调，如图1(b)所示。(a)cohesive单元与其他单元公用节点(b)独立的网格

2、通过“tie"绑定图1建模方法上述两种方法都可以用来模拟复合材料的分层失效，第一种方法划分网格比较复杂；第二种方法赋材料属性简单，划分网格也方便，但是装配及“tie"很繁琐；因此在实际建模中我们应根据实际结构选取较简单的方法。材料属性应用cohesive单元模拟复合材料失效，包括两种模型：一种是基于traction-separation描述；另一种是基于连续体描述。其中基于traction-separation描述的方法应用更加广泛。而在基于traction-separation描述的方法中，最常用的本构模型为图2所示的双线性本构模型。它给出了材料达到强度极限前的线弹性段和

3、材料达到强度极限后的刚度线性降低软化阶段。注意图中纵坐标为应力，而横坐标为位移，因此线弹性段的斜率代表的实际是cohesive单元的刚度。曲线下的面积即为材料断裂时的能量释放率。因此在定义cohesive的力学性能时，实际就是要确定上述本构模型的具体形状：包括刚度、极限强度、以及临界断裂能量释放率，或者最终失效时单元的位移。常用的定义方法是给定上述参数中的前三项，也就确定了cohesive的本构模型。Cohesive单元可理解为一种准二维单元，可以将它看作被一个厚度隔开的两个面，这两个面分别和其他实体单元连接。Cohesive单元只考虑面外的力，包括法向的正应力以及XZ,YZ两个方向的剪应力。

4、下文对cohesive单元的参数进行阐述，并介绍参数的选择方法。tradiohCohesive单元的刚度基于traction-separation模型的界面单元的刚度可以通过一个简单杆的变形公式来理解(1)PL5=AE其中L为杆长，E为弹性刚度，A为初始截面积，P为载荷。公式（1）又可以写成(2)5=丄K其中S=pA为名义应力，K=EL为材料的刚度。为了更好的理解K,我们把K=EL写成:(3)K=E=EL=ELL1这里我们用L来代替1,其中L可以理解为建模厚度，即建模时cohesiveinterface的几何厚度；L'为实际厚度，即cohesiveinterface的真实厚度，这个厚度

5、在cohesivesection中定义。EL可以理解为几何刚度，即模型中cohesiveinterface所具有的刚度；=厂为cohesiveinterface的真实刚度。当L为1时，计算界面刚度就采用几何刚度EL，当L为时，计算时界面刚度变为1000EL。举个小例子，如果界面的实际厚度为，而在建模时就是按照这个厚度建立的，在定义material-section时又specify这层的厚度为，实际上就等于把界面刚度提高了2个数量级，模拟结果当然是不对的，这时定义section时应采用默认厚度1。ABAQUS在cohesive建模中使用了很“人性化”的设计，实际问题中界面可能很薄，有的只有0.0

6、01mm，甚至更小。有些问题cohesive单元的interface还可能是0厚度（比如crack问题），而相对来说整体模型也许很大，如果不引入这两个厚度，我们就要在很大的模型中去创建这个很小的界面这是一个很麻烦的事情。引入这两个厚度，在建模时我们就可以用有限的厚度来代替这个很小的界面厚度，只要在section中定义这个L就好了。（注：以上大部分内容来自仿真论坛：再议cohesive应用中对于一些参数的理解）一个解释“另外有个我的经验公式：大体上energy>*(damageinitiation)A2/(stiffness)这个公式不难理解，就是锐角三角形的总面积大于一条侧边下的面积，将

7、traction-separationlaw画成图线你就一目了然了。不过根据不同的法则，会稍微有些区别的。”一-以上的话引自dava的个人空间，这里我想解释下这个不等式有些新手可能一下还看不明白。damageinitiation为开始破坏时的应力，即三角形的高；stiffness为刚度，也就是斜率，即tanq；所以侧边三角形的底边为damageinitiation/stiffness,*(damageinitiation)A2/(stiffness)即为侧边下的三角形面积。实际上能量还要大于这个侧边下三角形的面积很多,因为斜率一般都很大。关于材料参数定义cohesive的材料时，要填入材料的参

8、数，这些材料参数是材料固有的特性，与几何没有关系，所以放心大胆的填入吧。材料参数是由试验得到的，如果不能做实验(多数情况如此)，就去查国际上相关的文献吧，数据甚至比你自己做试验都要详细，在填入数据时要注意单位的统一。再说句，断裂能为单位面积上的能量，如你的单位选取N(力的单位)和M(长度单位)，那么能量的单位为N/M。下面举例来说明cohesive单元刚度的设置过程，以为例：进入property界面，点击MaterialTCreat,在弹出的EditMaterial对话框中，可以编辑新创建的cohesive材料的名称，然后点击MechanicalTElasticityTElasticTTrac

9、tion,在空格中输入相应的刚度。P威刮于亡丸山宙笈:|U何-2：.口刚剧：回®O令画匕山汗二二山代1234A*1:%¥百站yE-j.XiaH.d*L-L国oIjtkl半AmaLa.11oiix|Arihlycics£Jois0电gidtf-U-FLiyJri-:£c<-l4Ca<nLL4LT11直Jiold.：Iidil-7JirL：Mui*：:|llflJaL-L.|Jul:pyL-L.|f：TtcprAyJi£a.J.t3*匕口尸匸匸口丄1卫34厶IlTlIE>!Ut-jtlVi»flp4Tt*imI*li1i

10、i-iklgrtLE>*ifiL-ivipsL*fci*Jrcafn*1FAbgh'i-FlajriweHij苗iLissasHr回m已：。心曲.耻es“Eb医皿也爾.barlB*iSH-dQ5丢LDLF4BiiLjnnMjiit&UcasMi-jlyds隔压)AdipiLTi1rJriISC旷HCibllligSrnI.xa宝宅-P門丄血+.木"叵jFi.liUtdtlViwij4rtViml*li*l«ri»1SactLorilrifiL«Cn(iLt»丄EtjpSpiriIa*tw»TaiIcPIuj-iw

11、eHijUfig卿中FQq侖tl|乍连|Hl底口曲临日前心<0匚眉|*'1F'rtpBrliSFdafaulti*i23理HHuoltsPcjult:rr-sj-triro怔L:価血-1Ttfl：|tE-L1吁:lidtla：la飾通-【ffii.npSiEtLlyili晋uatLiicuI-|JAxhlycai9£Mt(S)嗒3旳¥"切计Trit書1Cci-»XKi.Li,sai图3.cohesive单元刚度的定义损伤准则初始损伤准则初始损伤对应于材料开始退化，当应力或应变满足于定义的初始临界损伤准则，则此时退化开始。Abaqus

12、的Damagefortractionseparationlaws中包括：QuadeDamage、MaxeDamage、QuadsDamage、MaxsDamage、MaxpeDamage、MaxpsDamage六种初始损伤准则，其中前四种用于一般复合材料分层模拟，后两种主要是在扩展有限元法模拟不连续体（比如crack问题）问题时使用。使用图2所示的双线本构模型，其中：to、to及to分别代表纯I型、纯II型或纯III破nst坏的最大名义应力，80、80，80代表相应的最大名义应变，当定义界面单元的初始厚度为nst1时，则名义应变等于与之相对应的相对位移§，5及5。nstQuadeDa

13、mage为二次名义应变准则：当名义应变比的平方和等于1时，损伤开始。8nMaxeDamage为最大名义应变准则：当任何一个名义应变的比值达到1时，损伤开始。maxQuadsDamage为二次名义应力准则：当各个方向的名义应变比的平方和等于1时，损伤开始。2MaxsDamage为最大名义应力准则当任何一个名义应力比值达到1时，损伤开始。maxKIfl-i虫3aLTi#tp£ftIS*?孚d<al丄讹U.姐鬥C-：*pi£l.t4&E5I.0L軒4!曲11JftitlF-a£1415TllSTilS；加Ip电iEspirlFskf帕lit占4-*LtJp

14、DLi口丄23410n昌"：如障出二门:fanOtM出二I116!：|loUl-lMFirt：|Firb-lMl|1010staijc-jjT田irepim厶择hn*a.yii;i.Ka£Jotsd)“AlxitLidtyFtx扎Ca-«E；«csil;imi图4.初始损伤准则定义EditMaterial对话框中，点击MechanicalTDamageforTractionSeparationLaws,然后根据自己的需要点击相应的损伤准则。其中最常用是QuadsDamage。损伤演化规律选择了初始损伤准则之后，然后点击SuboptionsTDamageE

15、volution,窗口如图5所示。其中Type包括Displacement和Energy,Displacement为基于位移的损伤演化规律，而Energy为基于能量的损伤演化规律。Softening中包括Linear，Exponential及Tabular三种刚度退化方式DamageEvolution中的所有的选项都是用来确定单元达到强度极限以后的刚度降阶方式。一般常用：以能量来控制单元的退化，即TypeTEnergy；线性软化模型，即SofteningTLinear，DegradationTMaximum；MixedmodebehaviorTBK，ModemixratioTEnergy,并选

16、中Power。flM團JjJLfcTnbmtITairl*lat«r-itlJr-afiL-iCnjrati-XteaQratlEjuLBE片？:uJH中（*，AHQdLtlattrialEh-aT-M-Lfi+tnJhEv|e“.llsiitL”Ex£1i：3：latM-a酋，”b|£iba.LX0g打4i±jg荡sLiyu-iAroMEthHnr二i特ivudris3AJita闻%AdijltvitsrFmrCDIEnVIJ«5CTipU.iWf*X+ThlHl-3k.U.lChlijthrWQaadtDui_iIfi-iIJkliiHdf

17、libtirb£fl+ldVflnlithL-hLUtsfLl4sl<9fFiarsEl*Ct3CIiut4pli-n-rl0*图5损伤演化规律定义1.3Cohesive单元界面属性还是在Property界面中，点击SectionTCreate，在弹出的EditSection对话框中，选择OtherTCohesive。EditSectionLitygoryTyre<)Soli!G-asketOShellCoheEiveBe=dirijVciMisticA.COUEtlCirLfirLitirLtrfaceName：COHESIVE>?iOtherCntinuH.C：

18、=LTLcelMaterial:cohesiveCreate：.Response:TracticmSeparationVTrtit.ioltKizI'h-lc：gg:；Use：匚imde：tiixLtN：=diie：COKESITETjpe:Cokesive；Usenodalcoordinate：sCreateSec±ioiL图6.定义材料的界面属性在EditSection对话框中，在material的下拉菜单中选择刚才创建的cohesive材料,也可以点击右侧的create创建一组新的材料；Response选择tractionseparation。Initialthickn

19、ess为前文提到的L，默认值为1,也可以在specify中指定一个特定的值。将所创建的界面属性赋予几何实体点击AssignTSection,然后在视图中选中要赋的几何实体，点击左下角的Done，则弹出如下窗口，在窗口是Section中下拉选中所创建的Cohesive截面，点击OK,操作完成。图7.给实体赋截面属性Cohesive单元网格划分Cohesive单元网格的划分与其他单元基本一致，但是以下几点不同与其他单元，划分网格时应特别注意。 网格密度，cohesive单元的网格尺寸不能太大，通常需要比较精细的网格，不然容易引起收敛性问题，甚至无法继续计算。 必须使用sweep（扫掠）划分网格的方

20、法，并且扫掠的方向垂直于cohesive面，即沿着cohesive单元的厚度方向。 单元种类的选择Ridd.Eiaxultx“血&Vgiv:I弓掘PhAilx(3)占U詔*IT田lie逐xEdyiri由事Jibsdfi曲AiijltvitpFracCirr-BtutiMie営醍T人ElwwlXIWm-k和isrdtfrd0FMrlLcLL匐laid*TnVimpNark.Ad*ptmtj<J*i1ut-iTi-J-tFI-hj-itt也)p科IE.于(弋住因心看創忒固3施西3,01尬動富EE",UkKAbbnJhrw'Itat/nndA1?3©石Cio

21、ilgd7#ipwr«itars5"Ji血h出i"Ae-zuLbi;ESaeIxmilltofld.:IrmI-Iv|Pljget:Ob«Uy®>»ri：llart-D：那I时【l.Fwl"I匚s：LIOE3R(#4itirt«ftrstar1-13*0-htTeEvl4：£i«h.gl<wt卅闻qfar和53旳esalici:iFUslr：口mtrolsVrrstQarfeij>*4onEr图单元种类选择在单元库中选择cohesive,可以在Viscosity,specify中指定一粘性系数，来改善收敛性,但是粘性系数的设置不能太大，不然会影响计算结果，我们一般设置为；Elementdeletion：用于设置单元的删除情况，一般选yes,即当单元完全失效时被删除；maxdegradation：般设置为1,即当SDEG=1时，认为单