#ABAQUS围线积分+网格重划分--模拟裂纹任意路径扩展

总体分析结果.jpg

前面的帖子已经介绍过在ABAQUS中模拟裂纹都有哪些技术,事实上,目前各大商业软件对于裂纹扩展的模拟还是有很多很多限制的,更不要说模拟裂纹的任意路径扩展了,截止目前,分析裂纹任意路径扩展的方法主要包括以下几种方法:
1 使用ABAQUS(或者ls-dyna等)自带材料损伤本构,达到失效的单元会被删除;
   例如:brittle cracking、ductile damage等
2 使用ABAQUS自带的扩展有限单元法xfem模拟裂纹任意路径扩展,裂纹可以穿过单元;
   例如:基于LEFM或者粘性片段法的xfem
3 在实体单元间批量插入cohesive单元模拟裂纹的任意路径扩展;
   例如:在所有实体单元间批量插入cohesive单元
4 使用围线积分+网格重新划分模拟裂纹的任意路径扩展;
   例如:自己编程实现围线积分+网格重新划分,或者franc2d/3d,zencrack3d,ALOF,adapcrack等
********
前三种其实目前用的比较多,各种例子视频讲解都有,但是对于第四种方法,往往由于使用起来太复杂被大家所抛弃,下面我们就重点讲解一下这种方法,为大家提供一种思路,有兴趣的可以自己使用python二次开发来完成。

首先,第四种模拟裂纹任意路径扩展的方法的思想主要是:
初始裂纹在时间t0达到断裂韧性,裂纹发生扩展,删除原始模型网格,裂纹往前扩展一个增量,重新划分网格,把上一步的应力场使用*map solution 映射到新模型上,载荷施加t0之后的载荷段,再计算应力强度因子,如果在t1时刻裂纹尖端的应力强度因子又达到了断裂韧性,裂纹会进一步扩展......后面不断判断循环进行下去。

1.jpg

我们采用这种思想通过python二次开发实现了裂纹沿着任意路径的自动扩展,而且将程序得到的结果与理论解析解进行了对比,结果完全一致(这里不再给出)。
下面我们将采用该程序模拟一个单裂纹的任意扩展:

2.jpg

当我们只考虑crack1裂纹扩展的时候(这时候我们采用的载荷及约束的位置是不同于上图的,主要是展示一下裂纹扩展而已),结果如下(只象征性地给出裂纹扩展过程中的几个阶段云图):

mesh1.jpg
mesh5.jpg
mesh10.jpg
mesh15.jpg
mesh19.jpg
r1.jpg
r5.jpg
r10.jpg
r15.jpg
r19.jpg

裂纹扩展过程中相关数据的变化:

裂纹扩展过程中相关数据变化1.jpg
裂纹扩展过程中相关数据变化2.jpg
裂纹扩展过程中相关数据变化3.jpg

当我们只考虑crack1和crack2裂纹扩展的时候,结果如下(只象征性地给出裂纹扩展过程中的几个阶段云图):

mm1.jpg
mm3.jpg
mm6.jpg
mm9.jpg
mm13.jpg
rr1.jpg
rr3.jpg
rr6.jpg
rr13.jpg

当我们只考虑crack1、crack2、crack3、crack4裂纹扩展的时候,结果如下(只象征性地给出裂纹扩展过程中的几个阶段云图):

mmmm1.jpg
mmmm3.jpg
mmmm4.jpg
mmmm5.jpg
mmmn6.jpg
rrrr1.jpg
rrrr3.jpg
rrrr4.jpg
rrrr5.jpg
rrrr6.jpg

几种不同情况下的裂纹扩展与预期是一致的,而且此方法可以输出裂纹扩展过程中的裂纹尖端的应力强度因子和能量释放率,对于分析脆性材料中的裂纹任意路径扩展而言是一种非常不错的办法。

裂纹扩展45次.jpg

ABAQUS断裂模拟收徒 ,保证快速学会各种ABAQUS断裂模拟方法    1200/人(将享有各种插件以及程序,价值3000+、专门定制视频、全程亲自教学、各种模型调试及解答问题等等,倾囊相教)

如果有需要或者需要讨论的可以加qq1057593923,希望此帖子可以给大家提供一定的帮助。

#ABAQUS围线积分+网格重划分--模拟裂纹任意路径扩展的评论4条

#ABAQUS围线积分+网格重划分--模拟裂纹任意路径扩展的相关案例教程

1. ABAQUS裂纹体的建立 在ABAQUS中,根据裂纹尖端的形状可以建立两种类型的裂纹体:尖锐(Sharp)型裂纹和钝形(Blunt)裂纹,如图1所示。对于尖锐型裂纹,裂纹尖端存在奇异性;而钝形裂纹尖端可以看做是一个具有给定缺口半径的缺口,因此裂纹尖端不存在奇异性,可以按照常规有限元的建模方式来建立。在ABAQUS中,两种类型的裂纹均可以进行应力强度因子分析。 图1 尖锐型裂纹和钝形裂纹 对于
Embedded crack model 通过强不连续性方法(Strong Discontinuity Approach)将 interface 模型内嵌在实体单元中。在连续体单元中,形函数矩阵发生改变,导致在积分点处出现“断裂”。相比起用 interface单元连接实体单元的方式模拟开裂,Embedded crack model 在建模上更加方便,可以直接采用有限元软件里的网格划分功能,不需要另
工程实际中经常遇见材料开裂问题。对于受到外力作用的结构,材料的承载力随着结构的变形发生改变,比如延性材料。因此可以采用连续体塑性力学框架对结构进行分析模拟。但是对于脆性材料,连续体塑性力学尽管可以模拟出结构宏观的力-位移关系,但是很难模拟出材料的开裂破坏。 为了采用有限元模拟材料的开裂,cohesive单元常被预设在连续体单元之间,通过cohesive塑性或者损伤本构模拟出材料强度下降的过程。但是
In model: Model-1, XFEM crack growth interaction "Int-2"references a stationary XFEM crack "Crack-1". Delete or suppressthe interaction, or change the crack to allow growth.输入文件未生成,作业未提交分析 想请教一下,这个应该怎
The output variables listed below are available in Abaqus/Explicit. Mechanical analysis–nodal quantities CFORCE Field: yes History: no .fil: no Contact normal force (CNORMF) and frictional shear force
博士/博士生
影响力
粉丝
内容
获赞
收藏
    项目客服
    培训客服
    4 12