基于某VCCT地裂纹扩展模拟.docx

基于某VCCT地裂纹扩展模拟.docx

《基于某VCCT地裂纹扩展模拟.docx》由会员分享，可在线阅读，更多相关《基于某VCCT地裂纹扩展模拟.docx（24页珍藏版）》请在冰豆网上搜索。

基于某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,1

CGROW,CPATH,cpath1

CGROW,FCOPTION,MTAB,5

CGROW,DTIME,1.0e-4

CGROW,DTMIN,1.0e-5

CGROW,DTMAX,2.0e-3

...

12.1.2裂纹扩展

在裂纹扩展模拟中，一个关心的量是裂纹扩展量。

VCCT方法测量裂纹扩展是基于已经张开的界面单元的长度，如下述方程和图片所示：

图12.3二维和三维裂纹扩展

对于二维问题，裂纹扩展是当前已经张开的界面单元的长度之和（a）。

对于三维问题，裂纹扩展在每个裂纹前沿节点进行测量，为沿着裂纹扩展方向的界面单元边长的和（b）.

裂纹扩展量

（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）作为断裂准则。

总的能量释放率是三种模式的能量释放率的和，表述为：

其中

为临界能量释放率。

对于I型断裂模式，断裂准则简化为：

能量释放率准则是最简单的断裂准则，适用于2D和3D的断裂扩展模拟。

示例12.1临界能量释放率输入

gtcval=10.0

CGROW,FCOPTION,GTC,gtcval

12.1.3.2.线性断裂准则

线性选项假设断裂准则是三种模式的能量释放率的线性函数。

表述为：

其中

，

，

分别为I型，II型，III型断裂模式的临界能量释放率。

这三个值通过命令TBDATA输入，如下：

Constant

TBDATAInput

Comments

C1

ModeI的临界能量释放率, 

>0

C2

ModeII的临界能量释放率, 

>0

C3

ModeIII的临界能量释放率, 

>0

示例12.2线性准则输入

g1c=10.0

g2c=20.0

g3c=25.0

TB,CGCR,1,,,LINEAR

TBDATA,1,g1c,g2c,g3c

三个参数不能同时为零。

如果其中一个设为零，相应的项被忽略。

当三个临界能量释放率相等，则线性断裂准则简化为临界能量释放率准则。

线性断裂准则适用于当三种断裂模式的临界能量释放率明显存在时的三维混合断裂模式。

12.1.3.3.双线性断裂准则

双线性断裂选项假设断裂准则是I型和II型断裂模式的能量释放率的线性函数，表达式为：

其中

，

分别为I型，II型断裂模式的临界能量释放率，ξ和ζ为材料常数。

四个值都能通过TBDATA命令定义为温度的函数。

如下所示：

Example 12.3  双向性准则输入

g1c=10.0

g2c=20.0

x=2

y=2

TB,CGCR,1,,,BILINEAR

TBDATA,1,g1c,g2c,x,y

双线性断裂准则适用于二维混合断裂模式的模拟。

12.1.3.4.B-K断裂准则

B-K选项表述为：

其中

，

，分别为I型，II型断裂模式的临界能量释放率,η为材料常数。

三个值都能够通过TBDATA命令定义为温度的函数。

如下所示：

B-K准则本来用于复合材料界面断裂，适用于三维混合断裂模式的模拟。

Example 12.4  B-K准则输入

g1c=10.0

g2c=20.0

h=2

TB,CGCR,1,,,BK

TBDATA,1,g1c,g2c,h

12.1.3.5.修正B-K断裂准则

修正B-K选项,表述为:

其中

，

，

分别为I型，II型，III型断裂模式的临界能量释放率,η为材料常数。

四个值都能够通过TBDATA命令定义为温度的函数。

如下所示：

=

时，修正B-K准则简化为B-K准则。

修正B-K准则用于复合材料界面断裂，考虑明显的II型和III型临界能量释放率，适用于三维混合断裂模式的模拟。

Example 12.5  修正B-K准则输入

g1c=10.0

g2c=20.0

g3c=25.0

h=2

TB,CGCR,1,,,MBK

TBDATA,1,g1c,g2c,g3c,h

12.1.3.6.指数断裂准则

指数准则选项假设断裂准则是三种模式的能量释放率的指数函数。

表述为：

其中

，

，

分别为I型，II型，III型断裂模式的临界能量释放率。

是指数，为常数。

六个量都可以通过TBDATA命令定义为温度的函数，如下所示:

三个临界能量释放率不能同时为零。

如果其中一个设为零，相应的项被忽略。

当指数

设置为1时，指数准则简化为线性断裂准则。

指数准则适用于当三种断裂模式的临界能量释放率明显存在时的三维混合断裂模式。

示例12.6指数准则输入

g1c=10.0

g2c=20.0

g3c=25.0

n1=2

n2=2

n3=3

TB,CGCR,1,,,POWERLAW

TBDATA,1,g1c,g2c,g3c,n1,n2,n3

12.1.3.7.用户自定义断裂准则

一个自定义的断裂准则可以表述为：

其中断裂准则是三种断裂模式的能量释放率和材料常数的函数。

所有的值通过TBDATA命令输入。

你必须提供一个子程序。

下面为一个子程序定义线性断裂准则的例子：

*deck,user_cgfcritoptimize

SUBROUTINEuser_cgfcrit（cgi,cid,kct,

&nprop,prop,fcscl,

&var1,var2,var3,var4）

c*****************************************************************

c

c***primaryfunction:

ccomputefacturecriterionforcrackgrowth

cuserfracturecriterionexample

c***notice:

thisroutinecontainssasiconfidentialinformation***

c

#include"impcom.inc"

#include"ansysdef.inc"

c

cinputarguments

c===============

ccgi（int,sc,in）CGROWsetid

ccid（int,sc,in）CINTIDtobeused

ckct（int,sc,in）Currentcracktipnode

cnprop（int,sc,in）numberofproperties

cprop（dp,ar（*）,in）propertyarray

c

cOutputarguments

c===============

cfcscl（dp,sc,ou）fracturecriterion

careturnvalueofoneorbigger

cindicatesfracture

c

cMisc.arguments

c===============

cvar1（,,）notused

cvar2（,,）notused

cvar3（,,）notused

cvar4（,,）notused

c

c*****************************************************************

c

c***subroutines/function

c***get_cgfpar:

APItoaccessfracturedata

c***wrinqr:

ansysstandardiofunction

c***

externalget_cgfpar

externalwrinqr

integerwrinqr

c***argument

c

INTEGERcgi,cid,kct,nprop

doubleprecisionfcscl,

&var1,var2,var3,var4

doubleprecisionprop（nprop）

c

c***localvariable

c

integerdebugflag,iott

integernn

doubleprecisiong1c,g2c,g3c,g1,g2,g3

doubleprecisiongs（4）,da

（1）

c

c***localparameters

DOUBLEPRECISIONZERO,ONE

parameter（ZERO=0.0d0,

&ONE=1.0d0）

c

c*****************************************************************

c***initialization

fcscl=ZERO

c***retrieveenergy-releaserates

c***forcrackcidandcracktipnodekct

c***gs（1:

3）willbereturnedasG1,G2,G3

c***getenergy-releaserates

nn=3

gs（1:

nn）=ZERO

callget_cgfpar（'GS',cid,kct,0,nn,gs

（1））

c***getcrackextension

nn=1

da

（1）=ZERO

callget_cgfpar（'DA',cid,kct,0,nn,da

（1））

c***energy-releaserates

g1=abs（gs

（1））

g2=abs（gs

（2））

g3=abs（gs（3））

c***inputpropertyfromTBDATA,1,c1,c2,c3

g1c=prop

（1）

g2c=prop

（2）

g3c=prop（3）

c***linearfracturecriterion

fcscl=ZERO

if（g1c.gt.TINY）fcscl=fcscl+g1/g1c

if（g2c.gt.TINY）fcscl=fcscl+g2/g2c

if（g3c.gt.TINY）fcscl=fcscl+g3/g3c

c***userdebugoutput

debugflag=1

if（debugflag.gt.0）then

iott=wrinqr（WR_OUTPUT）

write（iott,1000）cgi,cid,kct,da

（1）,fcscl,gs（1:

3）

1000format（5x,'userfracturecriterion:

'/

&5x,'crackgrowthsetID=',i5/

&5x,'crackID=',i5/

&5x,'cracktipnode=',i5/

&5x,'crackextension=',g11.5/

&5x,'calculatedfractureparameter=',g11.5/

&5x,'energy-releaseratesGs（1:

3）=',3g12.5）

endif

return

end

12.1.4裂纹扩展模拟例题

该例子采用双悬臂梁,一端有线裂纹。

在梁端裂纹上下施加相反方向的相同位移，来张开裂纹，如图所示：

图12.4双悬臂梁的裂纹扩展

下图显示了有限元网格：

采用PLANE182单元，打开enhancedstrain选项（keyopt

（1）=2）来建立实体模型。

INTER202单元来建立裂纹路径。

假设为平面应变条件。

在垂直方向采用6个单元，水平方向200个单元。

下图显示了预测的荷载-位移曲线。

图12.6双悬臂梁荷载-位移曲线

在裂纹开始扩展前，反力随着位移增大并迅速达到峰值。

然后在裂纹扩展的初始阶段，反力迅速下降，随着裂纹增长下降速度变慢。

结果和文献结果符合的很好。

下面给出最大主应力的云图。

图12.7双悬臂梁的云图

下面给出该双悬臂梁的裂纹扩展模拟的输入文件：

/BATCH

/TITLE,,CRACKGROWTHSIMULATIONOFADCBPROBLEM-2DPLANESTRAIN

/PREP7

DIS1=0.9

DIS2=12.0

N1=1000

N2=1000

N3=10

DL=100

DH=3

A0=30

NEL=200

NEH=6

TOLER=0.1E-5

ET,1,182!

*2D4-NODESTRUCTURALSOLIDELEMENT

KEYOPT,1,1,2!

*ENHANCESTRAINFORMULATION

KEYOPT,1,3,2!

*PLANESTRAIN

ET,2,182

KEYOPT,2,1,2

KEYOPT,2,3,2

ET,3,202!

*2D4-NODECOHESIVEZONEELEMENT

!

KEYOPT,3,2,2!

*ELEMENTFREEOPTION

KEYOPT,3,3,2!

*PLANESTRAIN

MP,EX,1,1.353E5!

*E11=135.3GPA

MP,EY,1,9.0E3!

*E22=9.0GPA

MP,EZ,1,9.0E3!

*E33=9.0GPA

MP,GXY,1,5.2E3!

*G12=5.2GPA

MP,PRXY,1,0.24

MP,PRXZ,1,0.24

MP,PRYZ,1,0.46

G1C=0.28!

*CRITICALENERGY-RELEASERATE

G2C=0.80

G3C=0.80

TB,CGCR,1,,3,LINEAR!

*LINEARFRACTURECRITERION

TBDATA,1,G1C,G2C,G3C

!

FEMODEL

RECTNG,0,DL,DH/2!

*DEFINEAREAS

RECTNG,0,DL,0,-DH/2

LSEL,S,LINE,,2,8,2!

*DEFINELINEDIVISION

LESIZE,ALL,DH/NEH

LSEL,INVE

LESIZE,ALL,,,NEL

ALLSEL,ALL

TYPE,1!

*MESHAREA2

MAT,1

LOCAL,11,0,0,0,0

ESYS,11

AMESH,2

CSYS,0

TYPE,2!

*MESHAREA1

ESYS,11

AMESH,1

CSYS,0

NSEL,S,LOC,X,A0-TOLER,DL

NUMMRG,NODES

ESLN

TYPE,3

MAT,5

CZMESH,,,1,Y,0,!

*GENERATEINTERFACEELEMENTS

ALLSEL,ALL

NSEL,S,LOC,X,DL!

*APPLYCONSTRAINTS

D,ALL,ALL

NSEL,ALL

!

ESEL,S,ENAME,,202!

*SELECTINTERFACEELEMENTTO

CM,CPATH,ELEM!

*DEFINECRACKGROWTHPATH

NSLE

NLIST

NSEL,S,LOC,X,A0

NSEL,R,LOC,Y,0

NLIST

ESLN

ELIST

CM,CRACK1,NODE!

*DEFINECRACKTIPNODECOMPONENT

NLIST

ALLS

FINISH

/SOLU