ABAQUS基于扩展有限元与内聚力行为的微细观结构断裂损伤力学分析

树脂砂广泛应用于铸造以及现代 3D 打印行业。原砂粒度及其分布、树脂加入量、树脂砂粘结桥有效横截面积，以及空隙分布等构成其微细观结构特征，这些特征与相关性能决定着树脂砂宏观力学行为。

1、研究目的：研究树脂砂微细观结构对宏观力学行为的作用规律，以及树脂砂粘结桥在载荷作用下的破坏过程与断裂机理，可以较深入的认识树脂砂结构常温断裂机理，并指导铸造工程实践，具有重要的理论和实践价值。

2、研究方法：针对原砂与树脂间弱界面的脱粘与树脂基体的失效这两种失效模式，构建了基于断裂力学三维应力微细观结构单胞模型。基于界面接触的内聚力行为方法构建三维微细观结构单胞模型，刻画了粘结砂粒/树脂界面失效，采用位移附加函数的 XFEM（Extended Finite Element Method，扩展有限元方法）捕捉树脂基体损伤与裂纹扩展。研究了树脂砂在拉伸载荷作用下粘结桥的破坏形式。研究结果见下面分解：

二、三维树脂砂微细观结构断裂有限元模型仿真

接下来，笔者为大家带来三维树脂砂微细观结构断裂有限元模型的构建过程。

第一步：Part模块，有限元模型的创建。

1、通过切分的方式创建1/8模型  

同样的方法创建另外一般大小不一样的1/8模型，如下：

2、在装配模块通过布尔运算创建PART3。PART包含两个PART1和两个PART2

3、在装配模块通过布尔运算创建PART4。PART4包含两个PART3 ,至此完成三维模型的创建。

第二步：Property模块，材料的创建，截面创建与赋予（两种材料，sand 与resin）

第三步：Assembly模块

只需装配PART4,默认装配即可

第四步：Step模块、step设置与变量输出设置

第五步：Interation模块、Interation设置（选用通用接触，也可采用面面接触等），建立参考点，耦合加载面。

1、创建接触属性，采用内聚力行为方法模拟截面失效，XFEM的定义

砂粒与树脂截面创建内聚力行为

2、耦合，创建一次性边界条件

XFEM的定义

选择树脂区域为可能损伤区域，材料已设定了XFEM的损伤起始判据与损伤演化方式，此处选择富集区。

第六步：Load模块、边界条件与载荷的设置，采用smooth step加载

采用对称边界条件，如下图所示三个面，右端面加载位移

第七步：Mesh模块、三个实体的网格划分

通过边布种子与总体布种子方法，网格划分结果如下图：

 第八步：求解模块

采用4个CPU计算，总时长约48h

第九步：后处理结果（载荷过程中裂纹的起始与演化）

分析：图所示为三维应力单胞模型在位移载荷作用下断裂过程的最大主应力云图。由图可知，初始加载时，最大应力位于树脂桥中心区域，如图（a）所示；

随着载荷增加，砂粒/树脂粘结界面应力不断增加，导致损伤出现的驱动力为 I 型的界面拉伸应力；在拉伸载荷的持续作用下，损伤不断积累，当损伤度达到 1 时，界面出现破坏，产生初始裂纹；位移载荷继续增加，界面初始裂纹沿砂粒/树脂粘结界面扩展，裂纹扩展一定程度后，砂粒曲面结构特殊性导致树脂砂结构特征应力由界面拉伸应力转变为界面局部剪切应力，界面裂纹扩展速度下降，如图（b）所示；

继续加载时，当界面裂纹面在粘结球面扩展一定程度后，树脂颈横截面积最小位置外缘应力达到树脂拉伸强度，局部出现树脂基体初始裂纹，继续加载，裂纹沿着树脂颈环向扩展，形成环状裂纹面（基体裂纹），如图（c）所示；

界面裂纹沿着砂粒/树脂粘结界面继续向前扩展，扩展速度较慢，基体裂纹向树脂颈内部扩展，裂纹尖端为一环状裂纹线。基体裂纹开始扩展后，结构的特殊性导致基体特征应力不是由加载方向拉伸应力主导，存在部分剪应力，且基体裂纹始终沿着裂纹富集区最大最应力方向扩展，即：裂纹不是一直沿着与加载垂直方向扩展，而是有一定的偏转，偏转角度与结构位置有关。随着载荷的继续增加，基体裂纹扩展至砂粒/树脂粘结界面时，界面裂纹与基体裂纹相遇，形成完整的树脂颈断裂面，粘结桥破坏，如图（d）所示。

以上就是笔者为大家带来的三维树脂砂微细观结构断裂有限元模型的构建过程。希望对学习型仿真秀工程师有所帮助，此外，本文相关模型文件可以在仿真秀官网或APP的文章附件下载，欢迎大家关注我在仿真秀专栏-Wind12306。

此外，受仿真秀平台的邀请，前不久笔者在仿真秀平台发布了一套原创视频教程，希望对裂纹学习者学习知识和增长技能有所帮助，该视频教程为笔者多年积累的经验和案例。

ABAQUS 裂纹入门提高18讲—案例实操精讲 多种裂纹模拟方法、实现形式

以下是课程安排：