abaqus在岩土工程中的应用_ABAQUS软件中扩展有限元（XFEM）的应用

机械设备的破坏断裂和裂纹扩展是密切相关的，造成裂纹扩展的原因复杂多样，但是主要归结于裂纹尖端在非常小的区域内的应力、应变场的变化以及微观结构的损伤。但由于三维断裂问题是复杂的不连续问题，给断裂的数值模拟造成一定的困难。

---

传统的数值模拟方法采用的是基于连续介质理论的有限单元法进行计算，该方法需要对裂纹尖端网格加密、或引入奇异单元，造成操作复杂等问题。由于上述方法在计算时都有一定的困难。因此在传统有限元法基础之上发展起来了扩展有限元法(XFEM)的裂纹计算方法，该方法是在有限元位移函数中加入了跳跃函数及裂尖增强函数。应用XFEM方法简化模型建模，不需要再对裂纹几何进行网格划分，裂纹可以在单元边的内部通过，且XFEM可不设置裂纹扩展面，模型在加载过程中在最危险点处开始断裂，可以更准确的模拟真实的裂纹扩展情况。

---

本文就基于ABAQUS平台，应用XFEM方法对楔形预裂纹的扩展过程进行分析，希望对大家的断裂数值模拟提供思路。

1. 几何模型的的建立

简单模型使用ABAQUS的CAE模块可以很快捷的进行建模，但对于较为复杂的结构，ABAQUS的CAE 模块就较为繁琐。因此，在建模过程中可采用其它的前处理软件代替ABAQUS的这部分功能，从而减少工作量，本文应用solidworks三维软件建模。

2.材料和截面属性创建

设置楔形预裂纹模型材料，输入弹性模量为及泊松比为。在定义属性时可设定损伤黏度系数为1×10e-5来提高计算的收敛性。接着创建截面属性，先把材料属性赋予截面属性，再把截面属性赋予几何模型。

3.网格划分

可使用ABAQUS自带网格划分模块或使用 HYPERMESH 等网格处理软件，网格单元的质量和数量对求解结果和求解过程影响较大，如果结构单元全部由等边三角形、正方形、正四面体、立方六面体等单元构成，则求解精度可接近实际值，但实际中很难做到，因此我们就要进行模型简化或其他单元格形式。

4.分析步与预裂纹的定义

建立分析步，在研究裂纹扩展速度以及裂纹扩展路径时需打开几何非线性，在提取应力强度因子时需关闭几何非线性。适当增加最大增量步数，减小最小增量步大小，可降低计算的不收敛性。在相互作用模块设置类型为XFEM的裂纹，在研究裂纹扩展速度以及裂纹扩展路径时需允许裂纹生长，在提取应力强度因子时不允许裂纹生长。

5.载荷、边界条件定义及输出设置

在相互作用模块创建耦合约束，将螺栓孔内表面与中心点建立耦合，在耦合点加载集中作用力，将模型下表面固定。在提取应力强度因子时，历程输出中将作用域选定为裂纹，类型为应力强度因子。

6.求解

提交工作，可查看有无错误，并可监控整个计算过程，对于较复杂模型，计算时间会较长。

7.后处理

Abaqus求解完成后可以查看Mises应力云图及裂纹扩展路径等需要的数据或云图。