基于VCCT的裂纹扩展模拟

基于VCCT的裂纹扩展模拟虚拟裂纹闭合技术(VCCT)最初用于计算裂纹体的能量释放率。因此被广泛用于层合复合材料的界面裂纹扩展模拟，并假定裂纹扩展总是沿着预先定义的路径，特别是在界面处。基于VCCT的裂纹扩展模拟当前可用的线性单元如plane182和solid185。基于VCCT的裂纹扩展模拟包含下述假定：,裂纹扩展沿着预先定义的路径,路径通过界面单元来定义,分析为准静态分析，不考虑瞬态效应,材料为线弹性材料，可以是各向同性，正交各向异性，各向性。裂纹可以位于一种材料或者两种材料的界面。断裂准则基于采用VCCT方法计算的能量释放率。可采用多种断裂准则或自定义的准则。同一分析中可定义多条裂纹。VCCT裂纹扩展模拟使用：.界面单元INTER202(2D)和INTER205(3D)・CINT命令计算能量释放率•CGROW命令定义裂纹扩展集，断裂准则，裂纹扩展路径和求解控制参数。VCCT裂纹扩展模拟过程基于VCCT的裂纹扩展模拟假定为准静态模拟。下面为进行模拟的主要步骤：Stepi:建立预先定义裂纹路径的有限元模型Step2:进行能量释放率的计算

Step3:进行裂纹扩展计算裂纹扩展模拟为非线性结构分析，这里详述了一些特点，特别是裂纹扩展的分析细节。Stepl:建立预先定义裂纹路径的有限元模型标准的非线性求解过程需要建立有限元模型，有正确的求解控制设置，载荷和边界条件。预先定义的裂纹路径离散为界面单元，并分为一个单元组，如下图所示：图12.1采用界面单元离散裂纹路径界面单元可以通过CZMESH命令划分或者能生成界面单元的第三方工具划分。MPC约束单元选项(KEYOPT(2)=1)在裂纹扩展前把潜在的裂纹面绑定在一起。当满足断裂准则时，MPC约束随后释放，从而裂纹扩展。在二维问题中，裂纹尖端后的一个界面单元如果在一个指定的子步满足断裂准则则可能张开。在三维问题中，裂纹前沿后的所有界面单元如果满足断裂准则可能张开。裂纹尖端/前沿周围的单元尺寸影响能量释放率的计算精度。当程序采用修正算法,可能不能产生精确的结果。改为沿着裂纹扩展路径使用相同单元尺寸的网格。Step2：进行能量释放率计算基于VCCT的裂纹扩展模拟，必须首先进行能量释放率的计算。计算能量释放率，采用CINT,TYPE,VCCT命令，随后使用CINT命令指定其它选项比如裂纹尖端节点组和裂纹面/边的法向。VCCT计算采用下述假定：•当裂纹增加一个小量时，释放的应变能等于裂纹闭合相同的小量所需的能量•当裂纹扩展一个小量时，裂纹尖端（前沿）位置的裂尖场（变形）不变。当裂纹扩展接近边界或者两条裂纹彼此接近时假定不在适用。因此，使用VCCT计算要仔细检查分析结果。Step3进行裂纹扩展计算裂纹扩展计算在应力计算之后，solution步进行。为了进行裂纹扩展计算，必须先定义裂纹扩展集合，然后指定裂纹路径，断裂准则，裂纹扩展求解控制。求解命令CGROW定义裂纹扩展计算所有必需的参数。进行裂纹扩展计算步骤如下：Step3a:初始裂纹扩展集Step3b:指定裂纹路径Step3c：指定裂纹计算的ID和断裂准则Step3d：指定裂纹扩展的求解控制Step3a:初始裂纹扩展集定义裂纹扩展集，使用CGROW,NEW,n命令，其中n是裂纹扩展集的编号Step3b:指定裂纹路径定义裂纹路径，采用CGROW,PATH,cmname其中cmname是界面单元组的名称。Step3c：指定裂纹计算的ID和断裂准则指定裂纹计算ID,通过CGROW,CID,n命令，其中n是采用VCCT计算能量释放率的裂纹计算(CINT)的ID。(CINT命令定义的参数和断裂参数计算一致)对于简单的断裂准则，比如临界能量释放率，可以通过命令CGROW,FCOPTION,GC,指定，其中value为临界能量释放率。对于一些更复杂的断裂准则，可以通过材料数据表定义断裂准则。采用CGROW，FCOPTION，MTAB，matid命令，其中matid是材料表的材料ID号。有多种断裂准则可以使用，比如linear,bilinear,B-K,修正B-K，PowerLaw，和用户自定义的准则。更多信息请参阅TB，CGCR命令和Fracturecriterial.对于每个裂纹扩展集，只可以定义一条断裂准则和一个单元组。可以采用不同的断裂准则来定义多个裂纹扩展集。多条裂纹可以同时扩展或者彼此独立。当多条裂纹位于同一界面时，也可以合并为一条裂纹，如下图所示：图12.2裂纹扩展和合并f f+ —►Cracks也可以在各自的断裂扩展集中，对同一裂纹定义不同的断裂准则。裂纹可以基于

不同的准则（根据哪条准则满足）扩展，并且彼此独立。这种方法对于比较断裂机理很有帮助。使用CGROW命令定义求解控制参数如下:指定求解控制参数使用CGROW命令:断裂准则系数（fc）CGROW,FCRAT,value,Value为该系数初始时间步长（裂纹扩展开始时）CGROW,DTIME,value,Value为初始时间步长为避免过度预测承载能匕力，指定一个小的初始时间步长随后裂纹扩展的最小时间步长CGROW,DTMIN,value,value为最小时间步长大小随后裂纹扩展的最大时间步长CGROW,DTMAX,value,value为最大时间步长大小裂纹扩展前沿节点允许的最大裂纹扩展量CGROW,STOP,CEMX,value,value为允许的最大的裂纹扩展量裂纹扩展模拟很耗时，当达到感兴趣的指定裂纹扩展量时，使用该命令来终止分析当裂纹迅速扩展时（比如，裂纹扩展不稳定），使用较小的DTMAX和DTMIN来容许载荷重新平衡。当裂纹不再增长时，指定的时间步长控制被忽略，结果依赖于标准时间步长控制。示例:裂纹扩展集定义下面的输入示例定义一个裂纹扩展集：CGROW,NEW,1CGROW,CPATH,cpath1CGROW,FCOPTION,MTAB,5CGROW,DTIME,1.0e-4CGROW,DTMIN,1.0e-5CGROW,DTMAX,2.0e-3...12.1.2裂纹扩展在裂纹扩展模拟中，一个关心的量是裂纹扩展量。VCCT方法测量裂纹扩展是基于已经张开的界面单元的长度，如下述方程和图片所示：图12.3二维和三维裂纹扩展对于二维问题，裂纹扩展是当前已经张开的界面单元的长度之和（a）o对于三维问题，裂纹扩展在每个裂纹前沿节点进行测量，为沿着裂纹扩展方向的界面单元边长的和（b>裂纹扩展量（CEXT）是裂纹求解结果的一部分，和裂纹计算的ID号一致，可以和能量释放率一样通过P0ST1和POST26后处理命令（比如PRCINT,PLCINT,和CISOL）进行后处理。12.1.2断裂准则为建立裂纹扩展,必须定义裂纹开始和随后裂纹扩展的断裂准则。对于线弹性断裂力学，断裂准则通常假定为三种断裂模式的临界能量释放率的函数。表达为：对于一些模型可能需要其它的参数。当断裂准则满足时，发生断裂，表述为：其中“为断裂准则比率。推荐值为。.95至IJL05,默认为L0。可以使用下列断裂准则：临界能量释放率准则线性断裂准则双线性断裂准则B-K断裂准则修正B-K断裂准则指数断裂准则自定义断裂准则用户自定义选项需要提供子程序来定义你自己的断裂准则。12.13L临界能量释放率准则临界能量释放率准则使用总的能量释放率(GT)作为断裂准则。总的能量释放率是三种模式的能量释放率的和，表述为：其中为临界能量释放率。对于I型断裂模式，断裂准则简化为：示例示例12.2线性准则输入能量释放率准则是最简单的断裂准则，适用于2D和3D的断裂扩展模拟。示例12.1临界能量释放率输入gtcval=10.0CGROW,FCOPTION,GTC,gtcval12.1.3.2.线性断裂准则线性选项假设断裂准则是三种模式的能量释放率的线性函数。表述为：其中可，」，G「：分别为I型，II型，川型断裂模式的临界能量释放率。这三个值通过命令丁8口人丁人输入，如下：ConstantTBDATAInputCommentsGfC1ModeI的临界能量释放率，G【＞0GfiC2rcModeII的临界能量释放率，缶【【＞0GIIIC3ModeIII的临界能量释放率,G【U＞0g1c=10.0g2c=20.0g3c=25.0TB,CGCR,1,,,LINEARTBDATA,1,g1c,g2c,g3c三个参数不能同时为零。如果其中一个设为零，相应的项被忽略。当三个临界能量释放率相等，则线性断裂准则简化为临界能量释放率准则。线性断裂准则适用于当三种断裂模式的临界能量释放率明显存在时的三维混合断裂模式。12.1.3.3.双线性断裂准则双线性断裂选项假设断裂准则是I型和II型断裂模式的能量释放率的线性函数，表达式为：其中可，G］：分别为I型，II型断裂模式的临界能量释放率，£和Z为材料常数。四个值都能通过TBDATA命令定义为温度的函数。如下所示：ConstantTBDATAInputConindentsCjcriticalMode1energy-refease「ate』=口球CriticalMode]]■energ^-reieaserate..匕”>0C33G白口 Example12.3双向性准则输入g1c=10.0g2c=20.0x=2y=2TB,CGCR,1,,,BILINEARTBDATA,1,g1c,g2c,x,y双线性断裂准则适用于二维混合断裂模式的模拟。B-K断裂准则B-K选项表述为：f= 13n冏一时［智J其中G；,Gj；,分别为I型，11型断裂模式的临界能量释放率,n为材料常数。三个值都能够通过TBDATA命令定义为温度的函数。如下所示：ConstantTBDATAinputCommeots嘴ciCribcdMode[erfirgy'-rde5se「畛"'>0时C2CribcslMo*ITerergy-reba^raberl-=，11>0nn>oB-K准则本来用于复合材料界面断裂，适用于三维混合断裂模式的模拟。Example12.4B-K准则输入g1c=10.0g2c=20.0h=2TB,CGCR,1,,,BKTBDATA,1,g1c,g2c,h修正B-K断裂准则修正B-K选项，表述为:其中G；,G*5二匚分别为I型，11型，111型断裂模式的临界能量释放率,n为材料常数。四个值都能够通过TBDATA命令定义为温度的函数。如下所示：constantTBDATAinputcomrnerts-iclCribcdM3de]Erergy-releaserate,⑶>0师C2CriMaiMade]]ere「gy-rel635urate」"">0斯C3oiocaiMode]]]tnergy-rdeastrate,>o-1C4q>0Gi：二G%时，修正B-K准则简化为B-K准则。修正B-K准则用于复合材料界面断裂，考虑明显的II型和III型临界能量释放率，适用于三维混合断裂模式的模拟。Example12.5修正B-K准则输入glc=10.0g2c=20.0g3c=25.0h=2TB,CGCR,1,,,MBKTBDATA,1,g1c,g2c,g3c,h指数断裂准则指数准则选项假设断裂准则是三种模式的能量释放率的指数函数。表述为：r / \n2 >f=,_5l .ku-nC p.0Jk^ll) kulll>其中G；，G］：，Gj；：分别为I型，II型，III型断裂模式的临界能量释放率。口.， , 是指数，为常数。六个量都可以通过TBDATA命令定义为温度的函数，如下所示：ConstantTBDATA[routComments*elGLCriticalMode]energy-releasermte,1>0战C2.qCCriticalMode]Iererqv-releaserater- >0C3CriticalModeinenergy-relaaserater'布>0niSJr.>0C5n2>0飞csn3>0三个临界能量释放率不能同时为零。如果其中一个设为零，相应的项被忽略。当指数吗， 设置为1时，指数准则简化为线性断裂准则。指数准则适用于当三种断裂模式的临界能量释放率明显存在时的三维混合断裂模式。示例12.6指数准则输入glc=10.0g2c二20。g3c=25.0nl=2n2=2n3=3TB,CGCR,L,,POWERLAWTBDATA,1,g1c,g2c,g3c,n1,n2,n3用户自定义断裂准则一个自定义的断裂准则可以表述为：其中断裂准则是三种断裂模式的能量释放率和材料常数的函数。所有的值通过TBDATA命令输入。你必须提供一个子程序。下面为一个子程序定义线性断裂准则的例子：*deck,user_cgfcritoptimizeSUBROUTINEuser_cgfcrit(cgi,cid,kct,& nprop,prop,fcscl,& var1,var2,var3,var4)c*****************************************************************cc***primaryfunction:c computefacturecriterionforcrackgrowthc userfracturecriterionexamplec***notice:thisroutinecontainssasiconfidentialinformation***c#include"impcom.inc"#include"ansysdef.inc"cc inputargumentsc ===============ccgi (int,sc,in)CGROWsetidcid(int,sc,in)CINTIDtobeused

cidckct(int,sc,in)Currentcracktipnodecnprop(int,sc,in)numberofpropertiescprop(dp,ar(*),in)propertyarrayccOutputargumentscfcscl(dp,sc,cfcscl(dp,sc,ou)ccfracturecriterionareturnvalueofoneorbiggerindicatesfractureccMisc.argumentsccvar1(,,)notusedcvar2(,,)notusedcvar3(,,)notusedcvar4(,,)notusedcc*****************************************************************cc***subroutines/functionc***get_cgfpar:APItoaccessfracturedatac***wrinqr:ansysstandardiofunctionc***external get_cgfparexternal wrinqrinteger wrinqrc***argumentcINTEGER cgi,cid,kct,npropdoubleprecisionfcscl,& var1,var2,var3,var4doubleprecisionprop(nprop)cc***localvariablecinteger debugflag,iottinteger nndoubleprecisiong1c,g2c,g3c,g1,g2,g3doubleprecisiongs(4),da(1)cc***localparametersDOUBLEPRECISIONZERO,ONEparameter(ZERO=0.0d0,ONE=1.0d0)cc*****************************************************************c***initializationfcscl=ZEROc***retrieveenergy-releaseratesc***forcrackcidandcracktipnodekctc***gs(1:3)willbereturnedasG1,G2,G3c***getenergy-releaseratesnn=3gs(1:nn)=ZEROcallget_cgfpar('GS',cid,kct,0,nn,gs(1))c***getcrackextensionnn=1da(1) =ZEROcallget_cgfpar('DA',cid,kct,0,nn,da(1))c***energy-releaseratesg1=abs(gs(1))g2=abs(gs(2))g3=abs(gs(3))c***inputpropertyfromTBDATA,1,c1,c2,c3g1c = prop(1)g2c = prop(2)g3c = prop(3)c***linearfracturecriterionfcscl=ZEROif(g1c.gt.TINY)fcscl=fcscl+g1/g1cif(g2c.gt.TINY)fcscl=fcscl+g2/g2cif(g3c.gt.TINY)fcscl=fcscl+g3/g3cc***userdebugoutputdebugflag=1if(debugflag.gt.0)theniott=wrinqr(WR_OUTPUT)write(iott,1000)cgi,cid,kct,da(1),fcscl,gs(1:3)1000format(5x,'userfracturecriterion:'/=',i5/&5x,'crackgrowthsetID=',i5/&5x,'crackID =',i5/&5x,'cracktipnode =',i5/&5x,'crackextension =',g11.5/&5x,'calculatedfractureparameter=',g11.5/&5x,'energy-releaseratesGs(1:3)=',3g12.5)endifreturnend12.1.4裂纹扩展模拟例题该例子采用双悬臂梁,一端有线裂纹。在梁端裂纹上下施加相反方向的相同位移，来张开裂纹，如图所示：图12.4双悬臂梁的裂纹扩展下图显示了有限元网格：

采用PLANE182单元，打开enhancedstrain选项(keyopt(l)=2)来建立头体模型。INTER202单元来建立裂纹路径。假设为平面应变条件。在垂直方向采用6个单元，水平方向200个单元。下图显示了预测的荷载-位移曲线。图12.6双悬臂梁荷载-位移曲线图12.6双悬臂梁荷载-位移曲线DisplacementA(mm)2;U2£在裂纹开始扩展前，反力随着位移增大并迅速达到峰值。然后在裂纹扩展的初始阶段，反力迅速下降，随着裂纹增长下降速度变慢。结果和文献结果符合的很好。下面给出最大主应力的云图。图12.7双悬臂梁的云图FU7TMD,1MCEiflLSQLBTIOrSTTP-LSUB=4TIbE=lSI 【际｝PowerGraphicsEF^CET-lAWRES-WatDtK-.9口OSNSMT-^&331SSMX-262.3C657.5^0957.5^09116.051LJ5.iO6H-1□15^5-5,0943XT-45,9143下面给出该双悬臂梁的裂纹扩展模拟的输入文件：/BATCH/TITLE,,CRACKGROWTHSIMULATIONOFADCBPROBLEM-2DPLANESTRAIN/PREP7DIS1=0.9DIS2=12.0N1=1000N2=1000N3=10DL=100DH=3A0=30NEL=200NEH=6TOLER=0.1E-5ET,1,182!*2D4-NODESTRUCTURALSOLIDELEMENTET,1,182KEYOPT,1,1,2!*ENHANCESTRAINFORMULATIONKEYOPT,1,1,2KEYOPT,1,3,2!*PLANESTRAINET,2,182KEYOPT,2,1,2KEYOPT,1,3,2!*PLANESTRAINET,2,182KEYOPT,2,1,2KEYOPT,2,3,2ET,3,202ET,3,202!*2D4-NODECOHESIVEZONEELEMENT!KEYOPT,3,2,2!*ELEMENTFREEOPTIONKEYOPT,3,3,2!*PLANESTRAINMP,EX,1,1.353E5!*E11=135.3GPAMP,EY,1,9.0E3MP,EY,1,9.0E3!*E22=9.0GPAMP,EZ,1,9.0E3MP,EZ,1,9.0E3!*E33=9.0GPAMP,GXY,1,5.2E3!*G12=5.2GPAMP,PRXY,1,0.24MP,PRXZ,1,0.24MP,PRYZ,1,0.46G1C=0.28G1C=0.28!*CRITICALENERGY-RELEASERATEG2C=0.80G3C=0.80TB,CGCR,1,,3,LINEAR!*LINEARFRACTURECRITERIONTBDATA,1,G1C,G2C,G3C!FEMODELRECTNG,0,DL,DH/2RECTNG,0,DL,0,-DH/2LSEL,S,LINE,,2,8,2LESIZE,ALL,DH/NEHLSEL,INVELESIZE,ALL,,,NELALLSEL,ALLTYPE,1MAT,1LOCAL,11,0,0,0,0ESYS,11AMESH,2CSYS,0TYPE,2ESYS,11AMESH,1!*DEFINEAREAS!*DEFINELINEDIVISION!*MESHAREA2!*MESHAREA1CSYS,0NSEL,S,LOC,X,A0-TOLER,DLNUMMRG,NODESESLNTYPE,3MAT,5CZMESH,,,1,Y,0,ALLSEL,ALLNSEL,S,LOC,X,DLD,ALL,ALLNSEL,ALL!ESEL,S,ENAME,,202CM,CPATH,ELEM!*GENERATEINTERFACEELEMENTS!*APPLYCONSTRAINTS!*SELECTINTERFACEELEMENTTO!*DEFINECRACKGROWTHPATHNSLENLISTNSEL,S,LOC,X,A0NSEL,R,LOC,Y,0NLISTESLNELISTCM,CRACK1,NODECM,CRACK1,NODE!*DEFINECRACKTIPNODECOMPONENTNLISTALLSFINISH/SOLURESC,,NONEESEL,S,TYPE,,2NSLE,SNSEL,R,LOC,XNSEL,R,LOC,Y,DH/2NSEL,R,LOC,Y,DH/2!*APPLYDISPLACEMENTLOADINGONTOPD,ALL,UY,DIS1NSEL,ALLESEL,ALLESEL,S,TYPE,,1NSLE,SNSEL,R,LOC,XNSEL,R,LOC,Y,-DH/2NSEL,R,LOC,Y,-DH/2!*APPLYDISPLACEMENTLOADINGONBOTTOMD,ALL,UY,-DIS1NSEL,ALLESEL,ALLAUTOTS,ONTIME,1CINT,NEW,1!*CRACKIDCIN