ABAQUS中的断裂力学及裂纹分析总结Word文档格式.docx

ABAQUS中的断裂力学及裂纹分析总结Word文档格式.docx

《ABAQUS中的断裂力学及裂纹分析总结Word文档格式.docx》由会员分享，可在线阅读，更多相关《ABAQUS中的断裂力学及裂纹分析总结Word文档格式.docx（8页珍藏版）》请在冰点文库上搜索。

随后在工程中发现了越来越多的LEFM无法解释的问题。

cohesivefracturemechnics开始引起更多的关注。

在研究以混凝土为代表的quassi-brittlematerial时，cohesivefracturemechnics提供了非常好的结果，所以在70年代到90年代，cohesivefracturemechnics被大量应用于混凝土研究中。

目前比较常用的方法主要fictitiouscrackapproach和effective-elasticcrackapproach或是称为equivalent-elasticcrackapproach.其中fictitiouscrackapproach只考虑了Dugdale-Barenblattenergymechanism而effective-elasticcrackapproach只考虑了基於LEFM的Griffith-Irwinenergydissipationmechanism，但作了一些修正。

做裂纹ABAQUS有几种常见方法。

最简单的是用debond命令,定义*FRACTURECRITERION,TYPE=XXX,

参数。

。

**

*DEBOND,SLAVE=XXX,MASTER=XXX,timeincrement=XX

0,1,

……

......

time,0

要想看到开裂特别注意需要在指定的开裂路径上定义一个*Nset，然后在*INITIALCONDITIONS,TYPE=CONTACT中定义

master,slave,及指定的Nset

这种方法用途其实较为有限。

例子如图

debondexample

另一种方法，在interaction模块，special,定义crackseam,网格最好细化，用collapseelement模拟singularity.这种方法可以计算J积分，应力强度因子等常用的断裂力学参数。

裂尖及奇异性定义:

在interaction-special,先定义crack,定义好裂尖及方向,然后在singularity选择:

midsidenodeparameter:

输入0.25,然后选Collapsedelementside,duplicate

nodes，8节点单元对应（1/r）+（1/r^1/2）奇异性。

这里midsidenodeparameter选0.25对应裂尖collapse成1/4节点单元。

如果midsidenodes不移动到1/4处,则对应（1/r）奇异性,适合perfectplasticity的情况.网格划分:

裂尖网格划分有一些技巧需要注意，partition后先处理最外面的正方形，先在对角线和边上布点，记住要点constraint,然后选第三个选项donotallowthenumberofelementstochange不准seed变化，密度可以自己调整.最里面靠近圆的正方形可以只在对角线上布点.也可以进一步分割内圆及在圆周上布点.里面裂尖周围的内圆选freemesh,elementtype选cps6或者cpe6，外面四边形选sweepmesh,elementtype选cps8或者cpe8,记住把quad下那个缩减积分的勾去掉。

mesh

vonmises

2Dcrack

vonmisesatcracktip

这种方法的几个值得注意的问题，见不少朋友问过。

主要是对断裂力学的理解问题。

1（为什么设置理想弹塑性（epp）分析的时候得到的xx,yy方向或者最大应力值Sxx,Syy会超过材料的屈服强度Sy呢,这分析结果可能吗,

这是因为在ABAQUS中对应等于材料的屈服强度的是vonMises等效应力Se=Sy，因此在平面应变的条件下，xx方向的应力Sxx=Sy*pi/SRQT（3）>

Sy,而Syy=Sy*（2+pi）/SRQT（3）,大概是3倍的屈服应力。

所以得到大于材料的屈服强度的xx及yy方向应力是正常的。

2（为什么设置collapseelement的时候对弹性分析在中间就一个点而要把单元边上的中点移到1/4处，但弹塑性分析却要在中间设置一圈点并且保持单元边上的中点位置不变呢,这个其实不是随便定的，在有限元中分析裂纹时，对弹性分析需要模拟裂尖1/SQRT（r）的奇异性，这样在把单元边上的中点移到1/4处后计算出来的等参单元拉格郎日型函数对应的ufield正好包含1/SQRT（r）项，事实上这一方法在断裂力学的数值模拟发展史上是很巧妙的一个发现，至今仍然被广泛采用。

至于理想弹塑性分析需要模拟裂尖1/r的奇异性,

这样大家都知道在把单元边上的点放在到1/2处后计算出来的正常的等参单元拉格郎日型函数对应的ufield包含1/r项,可以模拟弹塑性分析需要的裂尖1/r的奇异性。

所以在看似动手点几下就能实现的分析模式后面有很清楚漂亮的理论作支持。

也可以使用python控制seam尺寸,然后移动partition和网格，我也没尝试过。

但有一些学者有类似的结果:

FRANC/FAM-Asoftwaresystemforthepredictionof

crackpropagation.In:

JournalofStructuralEngineering26,No.1,1999,p

p.39-48.

再不就是用一些专业的断裂力学软件如zencrack，感兴趣的可以自己看看:

还有就是比较新的cohesiveelement单元。

我仔细读了ABAQUScohesiveelement的理论帮助，个人意见ABAQUS的cohesiveelement采用的是广泛应用于混凝土的类似fictitiouscrack的方法。

只考虑了Dugdale-Barenblattenergymechanism。

这其中softeninglaw的影响是非常重要的。

但ABAUQS似乎只提供了linear或者exponential的softeninglaw，复杂的本构关系还需要另想办法。

至於基于Griffith-Irwinenergydissipationmechanism的J-integral值可以在LEFM分析中单独算。

（ABAQUS用的是SuoZhigang和Hutchinson在1990一篇论文中提出的方法）目前cohesivefracturemechnics已经被应用于各种材料。

不过在使用到纳米或者更小数量级的研究中碰到了不少问题，可能需要结合位错和分子动力学的一些理论。

现有的cohesiveelement单元需要定义damageinitiation和evolution的准则，softening准则目前好像只有linear和exponential，但对一般材料也够用了。

然后通过设置后处理displaygroup可以看到裂纹扩展情况。

裂纹扩展不是ABAQUS的强项，目前比较方便的只能用cohesiveelement，我做过几个模型效果还可以，但对应的参数需要一定的实验数据支持，否则做出来了也不知道对不对。

要注意geometricthickness和constitutivethickness;

materialstiffness和interfacialstiffness的区别以及厚度与精度的影响。

Cohesive

element的核心主要是TS-Law，无论里面的数据如何选取，厚度如何变化，cohesiveelement的表现取决于TS-Law的定义。

具体dava的popo10已经给过详细的解释的讨论，可以搜索他们的帖子。

如图是我做的3dcohesiveelementinterfacialcrack的例子。

3Dinterfacialcohesiveelmentcrack

cohesiveelementlayerfailure

另外自己可以编写cohesiveUel,以便更加灵活的定义cohesiveelement的T-Sl

aw,也有不少人做过,基本流程如下，源程序抱歉我不方便贴上来。

CSUBROUTINEUEL（RHS,AMATRX,SVARS,ENERGY,NDOFEL,NRHS,NSVARS,

PROPS,NPROPS,COORDS,MCRD,NNODE,U,DU,V,A,JTYPE,TIME,DTIME,

2KSTEP,KINC,JELEM,PARAMS,NDLOAD,JDLTYP,ADLMAG,PREDEF,NPREDF,

3LFLAGS,MLVARX,DDLMAG,MDLOAD,PNEWDT,JPROPS,NJPROP,PERIOD）C

INCLUDE'

ABA_PARAM.INC'

C

DIMENSIONRHS（MLVARX,*）,AMATRX（NDOFEL,NDOFEL）,PROPS（*）,

1SVARS（*）,ENERGY（8）,COORDS（MCRD,NNODE）,U（NDOFEL）,

2DU（MLVARX,*）,V（NDOFEL）,A（NDOFEL）,TIME

（2）,PARAMS（*），

3JDLTYP（MDLOAD,*）,ADLMAG（MDLOAD,*）,DDLMAG（MDLOAD,*）,

4PREDEF（2,NPREDF,NNODE）,LFLAGS（*）,JPROPS（*）

DIMENSIONC_COOR（3,8）,R_COOR（3,8）,R_F（24）

DIMENSIONR_MATRX（24,24）,ANGLE（3,9）,R_DUX（4）,R_DUY（4）,R_DUZ（4）

DIMENSIONTN（4）,TT_2（4）,TT_3（4）,STFN（4）,STFT_2（4）,STFT_3（4）

DIMENSIONAX（3,2）,AY（3,2）,AZ（3,2）,AB（3,2）

DIMENSIONPMID（3,2）,R_U（3,8）,R_DU（3,8）

CINITIALIZELOCALVARIABLES

CINITIALIZETHESTIFFNESSMATRIXANDRESIDUALFORCEARRAYCCALCULATETHECURRENTELEMENTGEOMETRY

CDEFINEROTATIONANGLESBETWEENLOCALANDGLOBALSYSTEMCCOMPUTELOCALELEMENALCOORDINATESYSTEM

CCOMPUTECOORDINATESINLOCALCOORDINATESYSTEMCCOMPUTERELATIVEDISP.INLOCALCOOR.SYSTEMCCOMPUTELOCALSTIFFNESSANDNODALTRACTIONSC---------------------------------------------------------------

CCOMPUTETHELOCALSTIFFNESSANDNODALTRACTIONSFORNODEPAIR1AND5,2AND6,3AND7,4AND8

C---------------------------------------------------------------CCHECKIFCRITICALDISP.ISREACHED

CCHECKENERGYFRACTURECRITERION

C-----------------------------------------------------------------------------

CCALCULATETHERESIDUALFORCEANDTHESTIFFNESSMATRIXCENDOFMAINCODE

CSUBROUTINES

感觉掌握好以上这些方法这些对ABAQUS的断裂分析应该算是比较完善了，个人见解

难免有偏差，不足之处欢迎牛人补充。