ABAQUS中的断裂力学及裂纹分析总结

1ABAQUS中的断裂力学及裂纹分析总结（转自simwe）（1）做裂纹ABAQUS有几种常见方法。最简单的是用debond命令,定义*FRACTURECRITERION,TYPE=XXX,参数。。。***DEBOND,SLAVE=XXX,MASTER=XXX,timeincrement=XX0,1,……......time,0要想看到开裂特别注意需要在指定的开裂路径上定义一个*Nset，然后在*INITIALCONDITIONS,TYPE=CONTACT中定义master,slave,及指定的Nset这种方法用途其实较为有限。（2）另一种方法，在interaction模块，special,定义crackseam,网格最好细化，用collapseelement模拟singularity.这种方法可以计算J积分，应力强度因子等常用的断裂力学参数.裂尖及奇异性定义:在interaction-special,先定义crack,定义好裂尖及方向,然后在singularity选择：midsidenodeparameter:输入0.25,然后选Collapsedelementside,duplicatenodes，8节点单元对应(1/r)+(1/r^1/2)奇异性。这里midsidenodeparameter选0.25对应裂尖collapse成1/4节点单元。如果midsidenodes不移动到1/4处,则对应(1/r)奇异性,适合perfectplasticity的情况.网格划分:裂尖网格划分有一些技巧需要注意，partition后先处理最外面的正方形，先在对角线和边上2布点，记住要点constraint,然后选第三个选项donotallowthenumberofelementstochange不准seed变化，密度可以自己调整.最里面靠近圆的正方形可以只在对角线上布点.也可以进一步分割内圆及在圆周上布点.里面裂尖周围的内圆选freemesh,elementtype选cps6或者cpe6，外面四边形选sweepmesh,elementtype选cps8或者cpe8,记住把quad下那个缩减积分的勾去掉。这种方法的几个值得注意的问题，见不少朋友问过。主要是对断裂力学的理解问题。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的奇异性。所以在看似动手点几下就能实现的分析模式后面有很清楚漂亮的理论作支持。还有就是比较新的cohesiveelement单元。我仔细读了ABAQUScohesiveelement的理论帮助，个人意见ABAQUS的cohesiveelement采用的是广泛应用于混凝土的类似fictitiouscrack的方法。只考虑了Dugdale-Barenblattenergymechanism。这其中softeninglaw的影响是非常重要的。但ABAUQS似乎只提供了linear或者exponential的softeninglaw，复杂的本构关系还需要另想办法。至於基于Griffith-Irwinenergydissipationmechanism的J-integral值可以在LEFM分析中单独算。(ABAQUS用的是SuoZhigang和Hutchinson在1990一篇论文中提出的方法)目前cohesivefracturemechnics已经被应用于各种材料。不过在使用到纳米或者更小数量级的研究中碰到了不少问题，可能需要结合位错和分子动力学3的一些理论。现有的cohesiveelement单元需要定义damageinitiation和evolution的准则，softening准则目前好像只有linear和exponential，但对一般材料也够用了。然后通过设置后处理displaygroup可以看到裂纹扩展情况。裂纹扩展不是ABAQUS的强项，目前比较方便的只能用cohesiveelement，我做过几个模型效果还可以，但对应的参数需要一定的实验数据支持，否则做出来了也不知道对不对。要注意geometricthickness和constitutivethickness；materialstiffness和interfacialstiffness的区别以及厚度与精度的影响。Cohesiveelement的核心主要是TS-Law，无论里面的数据如何选取，厚度如何变化，cohesiveelement的表现取决于TS-Law的定义。具体dava的popo10已经给过详细的解释的讨论，可以搜索他们的帖子。如图是我做的3dcohesiveelementinterfacialcrack的例子。