Abaqus非均质煤岩体水力压裂模拟，案例加模型文件分享

大家可以先在知网搜索这篇论文：非均质弹塑性煤体水压致裂裂纹形态研究-中国知网 (cnki.net)。

模拟裂纹扩展的方法有很多，但本人认为FDEM方法是模拟缝网形态的最好方法，上面这篇论文使用Abaqus通过二次开发实现FDEM方法，相关方法星辰北极星大佬早就使用过了，这里我们主要解决两个问题：

1、材料非均质性的实现方法

2、多射孔多注入点的集中流量注入实现

模拟的基本模型很简单，办公室的工作站只够带动的二维的；因为考虑了地层的非均质性，所以计算不易收敛，这个二维模型为100w自由度，40核工作站算一步10分钟左右，单步增量时长1e-5到-7的样子，所以还请大家斟酌后建模。具体模型如下图：15m*15m的地层中有一簇射孔，详细尺寸图中已经给出啦。

这一簇射孔在际水力压裂中是一个分隔段，我们需要将中间四个注液点耦合到一个节点上，给这个节点设置集中注液流量，被耦合的四个点根据孔隙压力自己平衡每个点的注液流量。想实现上述目标需要用到管单元FP3D2，这个单元是不支持CAE操作的，关键字如下：

1）单元类型及参数定义

*Element, type=fp3d2, elset=elem-pipe

**

**单元编号**

**

*fluid pipe section, elset=elem-pipe, material=mat-flpipe

*fluid pipe flow loss, type=blasius, laminar flow transition = 1.0

0.031415, 0.2, 0., 0., 0., 0., 0.

2）材料定义

*material, name=mat-flpipe

*density, porefluid

1000

*viscosity, definition=newtonian

0.001

在CAE操作中先要建立一个“一入四出”的二维线单元，之后将线单元四个出口与地层的四个注入点绑定，之后的单元定义以及材料属性赋值操作就都要在inp文件中进行啦，每个参数啥意思帮助文档里有的哦。之后给我们的管单元的“一入”设置集中注入，他就会把流量分进“四出”后进入地层。

至于如何实现材料非均质性的描述，本文中使用的是Weibull分布，对模型中的Cohesive单元进行随机赋值，如何实现材料的随机赋值参照大佬，大佬的脚本讲解细致，可以自己学会的。

最终的模拟结果如下，效果还是不错的，想要模拟复杂缝网，考虑材料的非均质性必不可少！

付费压缩文件里有CAE文件，inp文件以及随机赋值的脚本。

说明：inp文件才是最终提交计算的文件，CAE里的模型是不能直接提交计算的

Python脚本的运行要在CAE中完成

CAE是2021版本建的模型

▼ 扫码获取案例完整内容 ▼

如果您有类似分析等项目仿真需求，可扫码添加客服回复【企服】咨询相关合作。

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删