基于VCCT的裂纹扩展模拟

基于VCCT的裂纹扩展模拟虚拟裂纹闭合技术(VCCT)最初用于计算裂纹体的能量释放率。因此被广泛用于层合复合材料的界面裂纹扩展模拟，并假定裂纹扩展总是沿着预先定义的路径，特别是在界面处。基于VCCT的裂纹扩展模拟,当前可用的线性单元如plane182和solid185。基于VCCT的裂纹扩展模拟包含下述假定：·裂纹扩展沿着预先定义的路径·路径通过界面单元来定义·分析为准静态分析，不考虑瞬态效应·材料为线弹性材料，可以是各向同性，正交各向异性,各向性。裂纹可以位于一种材料或者两种材料的界面。断裂准则基于采用VCCT方法计算的能量释放率。可采用多种断裂准则或自定义的准则。同一分析中可定义多条裂纹。VCCT裂纹扩展模拟使用：·界面单元INTER202（2D）和INTER205（3D）·CINT命令计算能量释放率·CGROW命令定义裂纹扩展集，断裂准则，裂纹扩展路径和求解控制参数。12.1.1VCCT裂纹扩展模拟过程基于VCCT的裂纹扩展模拟假定为准静态模拟。下面为进行模拟的主要步骤：Step1：建立预先定义裂纹路径的有限元模型Step2：进行能量释放率的计算Step3：进行裂纹扩展计算裂纹扩展模拟为非线性结构分析，这里详述了一些特点，特别是裂纹扩展的分析细节。12.1.1.1Step1：建立预先定义裂纹路径的有限元模型标准的非线性求解过程需要建立有限元模型，有正确的求解控制设置，载荷和边界条件。预先定义的裂纹路径离散为界面单元，并分为一个单元组，如下图所示：图12.1采用界面单元离散裂纹路径界面单元可以通过CZMESH命令划分或者能生成界面单元的第三方工具划分。MPC约束单元选项(KEYOPT(2)=1)在裂纹扩展前把潜在的裂纹面绑定在一起。当满足断裂准则时，MPC约束随后释放，从而裂纹扩展。在二维问题中，裂纹尖端后的一个界面单元如果在一个指定的子步满足断裂准则则可能张开。在三维问题中，裂纹前沿后的所有界面单元如果满足断裂准则可能张开。裂纹尖端/前沿周围的单元尺寸影响能量释放率的计算精度。当程序采用修正算法,可能不能产生精确的结果。改为沿着裂纹扩展路径使用相同单元尺寸的网格。12.1.1.2Step2：进行能量释放率计算基于VCCT的裂纹扩展模拟，必须首先进行能量释放率的计算。计算能量释放率，采用CINT,TYPE,VCCT命令，随后使用CINT命令指定其它选项比如裂纹尖端节点组和裂纹面/边的法向。VCCT计算采用下述假定：●当裂纹增加一个小量时，释放的应变能等于裂纹闭合相同的小量所需的能量●当裂纹扩展一个小量时，裂纹尖端（前沿）位置的裂尖场（变形）不变。当裂纹扩展接近边界或者两条裂纹彼此接近时假定不在适用。因此，使用VCCT计算要仔细检查分析结果。12.1.1.3.Step3进行裂纹扩展计算裂纹扩展计算在应力计算之后，solution步进行。为了进行裂纹扩展计算，必须先定义裂纹扩展集合，然后指定裂纹路径，断裂准则，裂纹扩展求解控制。求解命令CGROW定义裂纹扩展计算所有必需的参数。进行裂纹扩展计算步骤如下：Step3a：初始裂纹扩展集Step3b：指定裂纹路径Step3c：指定裂纹计算的ID和断裂准则Step3d：指定裂纹扩展的求解控制12.1.1.3.1Step3a：初始裂纹扩展集定义裂纹扩展集，使用CGROW,NEW,n命令，其中n是裂纹扩展集的编号12.1.1.3.1Step3b：指定裂纹路径定义裂纹路径，采用CGROW,PATH,cmname,其中cmname是界面单元组的名称。12.1.1.3.1Step3c：指定裂纹计算的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裂纹扩展和合并也可以在各自的断裂扩展集中，对同一裂纹定义不同的断裂准则。裂纹可以基于不同的准则（根据哪条准则满足）扩展，并且彼此独立。这种方法对于比较断裂机理很有帮助。使用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来容许载荷重新平衡。当裂纹不再增长时，指定的时间步长控制被忽略，结果依赖于标准时间步长控制。12.1.1.4.示例:裂纹扩展集定义下面的输入示例定义一个裂纹扩展集：CGROW,NEW,1CGROW,CPATH,cpath1CGROW,FCOPTION,MTAB,5CGROW,DTIME,1.0e-4CGROW,DTMIN,1.0e-5CGROW,DTMAX,2.0e-3...12.1.2裂纹扩展在裂纹扩展模拟中，一个关心的量是裂纹扩展量。VCCT方法测量裂纹扩展是基于已经张开的界面单元的长度，如下述方程和图片所示：∆a=∑∆i图12.3二维和三维裂纹扩展对于二维问题，裂纹扩展是当前已经张开的界面单元的长度之和（a）。对于三维问题，裂纹扩展在每个裂纹前沿节点进行测量，为沿着裂纹扩展方向的界面单元边长的和（b).裂纹扩展量∆a（CEXT)是裂纹求解结果的一部分，和裂纹计算的ID号一致，可以和能量释放率一样通过POST1和POST26后处理命令（比如PRCINT,PLCINT,和CISOL)进行后处理。12.1.2断裂准则为建立裂纹扩展,必须定义裂纹开始和随后裂纹扩展的断裂准则。对于线弹性断裂力学，断裂准则通常假定为三种断裂模式的临界能量释放率的函数。表达为：=对于一些模型可能需要其它的参数。当断裂准则满足时，发生断裂，表述为：其中为断裂准则比率。推荐值为0.95到1.05，默认为1.0。可以使用下列断裂准则：·临界能量释放率准则·线性断裂准则·双线性断裂准则·B-K断裂准则·修正B-K断裂准则·指数断裂准则·自定义断裂准则用户自定义选项需要提供子程序来定义你自己的断裂准则。12.1.3.1.临界能量释放率准则临界能量释放率准则使用总的能量释放率（GT）作为断裂准则。总的能量释放率是三种模式的能量释放率的和，表述为：=GTGTCT=++其中T为临界能量释放率。对于I型断裂模式，断裂准则简化为：=GIGIC能量释放率准则是最简单的断裂准则，适用于2D和3D的断裂扩展模拟。示例12.1临界能量释放率输入gtcval=10.0CGROW,FCOPTION,GTC,gtcval12.1.3.2.线性断裂准则线性选项假设断裂准则是三种模式的能量释放率的线性函数。表述为：=++其中，，分别为I型，II型，III型断裂模式的临界能量释放率。这三个值通过命令TBDATA输入，如下：ConstantTBDATAInputCommentsC1ModeI的临界能量释放率,0C2ModeII的临界能量释放率,0C3ModeIII的临界能量释放率,0示例12.2线性准则输入g1c=10.0g2c=20.0g3c=25.0TB,CGCR,1,,,LINEARTBDATA,1,g1c,g2c,g3c三个参数不能同时为零。如果其中一个设为零，相应的项被忽略。当三个临界能量释放率相等，则线性断裂准则简化为临界能量释放率准则。线性断裂准则适用于当三种断裂模式的临界能量释放率明显存在时的三维混合断裂模式。12.1.3.3.双线性断裂准则双线性断裂选项假设断裂准则是I型和II型断裂模式的能量释放率的线性函数，表达式为：其中，分别为I型，II型断裂模式的临界能量释放率，ξ和ζ为材料常数。四个值都能通过TBDATA命令定义为温度的函数。如下所示：Example12.3双向性准则输入g1c=10.0g2c=20.0x=2y=2TB,CGCR,1,,,BILINEARTBDATA,1,g1c,g2c,x,y双线性断裂准则适用于二维混合断裂模式的模拟。12.1.3.4.B-K断裂准则B-K选项表述为：其中，，分别为I型，II型断裂模式的临界能量释放率,η为材料常数。三个值都能够通过TBDATA命令定义为温度的函数。如下所示：B-K准则本来用于复合材料界面断裂，适用于三维混合断裂模式的模拟。Example12.4B-K准则输入g1c=10.0g2c=20.0h=2TB,CGCR,1,,,BKTBDATA,1,g1c,g2c,h12.1.3.5.修正B-K断裂准则修正B-K选项,表述为:其中，，分别为I型，II型，III型断裂模式的临界能量释放率,η为材料常数。四个值都能够通过TBDATA命令定义为温度的函数。如下所示：=时，修正B-K准则简化为B-K准则。修正B-K准则用于复合材料界面断裂，考虑明显的II型和III型临界能量释放率，适用于三维混合断裂模式的模拟。Example12.5修正B-K准则输入g1c=10.0g2c=20.0g3c=25.0h=2TB,CGCR,1,,,MBKTBDATA,1,g1c,g2c,g3c,h12.1.3.6.指数断裂准则指数准则选项假设断裂准则是三种模式的能量释放率的指数函数。表述为：其中，，分别为I型，II型，III型断裂模式的临界能量释放率。n1,n2,n3是指数，为常数。六个量都可以通过TBDATA命令定义为温度的函数，如下所示:三个临界能量释放率不能同时为零。如果其中一个设为零，相应的项被忽略。当指数n1,n2,n3设置为1时，指数准则简化为线性断裂准则。指数准则适用于当三种断裂模式的临界能量释放率明显存在时的三维混合断裂模式。示例12.6指数准则输入g1c=10.0g2c=20.0g3c=25.0n1=2n2=2n3=3TB,CGCR,1,,,POWERLAWTBDATA,1,g1c,g2c,g3c,n1,n2,n312.1.3.7.用户自定义断裂准则一个自定义的断裂准则可以表述为：=T其中断裂准则是三种断裂模式的能量释放率和材料常数的函数。所有的值通过TBDATA命令输入。你必须提供一个子程序。下面为一个子程序定义线性断裂准则的例子：*deck,user_cgfcritoptimizeSUBROUTINEuser_cgfcrit(cgi,cid,kct,&nprop,prop,fcscl,&var1,var2,var3,var4)c*****************************************************************cc***primaryfunction:ccomputefacturecriterionforcrackgrowthcuserfracturecriterionexamplec***notice:thisroutinecon