断裂接触模型fracture contact-model的细微变化

1 引言

昨天回顾过去的笔记【离散断裂网络(DFN)[P3]: fracture contact-model】时运行了一小段程序，发现原来的代码已经运行不下去，原来的这篇笔记写于2020年6月18日，距今差不多已经过了2年，在这两年中，PFC命令和fracture contact-model模型都发生了不少变化。这个笔记仅基于这一小段代码，记录了这些变化之处。

2 命令变化

(1) model large-strain on 

PFC和3DEC默认的计算模式是大应变【小应变和大应变】，因此这个命令可以有也可以没有；此外，旧格式的largestrain仍然可以使用。

(2) contact cmat default model

新的接触模型都以contact cmat开头，发现没有contact，直接以cmat default model的命令仍然有效。

相关文章，在仿真秀官网搜索：

FLAC3D 7.0 新特性简介(P4)---软件集成

"contact cmat default model linear"

离散断裂网络(DFN)[P4]: 创建一个合成岩体SRM

生成复杂颗粒形状的三种方法

一个阶段性的小结: 离散断裂网络DFN

PFC 之 attribute & property

岩石力学---从物理试验到数值试验

(3) model clean

新的格式必须加上model, 只用clean会产生错误信息。

(4) fracture generate fracture-count

离散断裂网络命令已经完全改变，旧格式的dfn generate nfrac不能再继续使用。

相关文章，在仿真秀官网搜索：

离散断裂网络Discrete Fracture Network (DFN）[P1]

离散断裂网络 (DFN) [P2]: fracture generate

岩石破碎(Rock Fragmentation)的数值模拟(block fragment compute)

简化和优化DFN---fracture combine命令

岩体统计模型---对数正态分布(log-normal distribution)

离散断裂网络[DFN.]FISH函数分类与解析(UDEC和3DEC)

离散断裂网络(DFN)[P4]: 创建一个合成岩体SRM

Fracman读取FLAC3D生成的离散断裂网络(DFN)模型

(5) fracture property

fracture property的所有关键字，例如sj_kn, sj_ks等必须加引号，否则会产生错误，这是一个强制性的规则。

fracture property 'sj_kn' 1e5 'sj_ks' 1e-4 'sj_coh' 0

(6) fracture contact-model model

与fracture property的规则类似，安装模型时，模型名称必须加引号，否则会产生错误，这是一个强制性的规则。另外，旧的格式useaperture已经不能再使用，必须使用use-aperture。

fracture contact-model model 'smoothjoint' install use-aperture

3 模型改进

现在，可以使用的fracture contact-model接触模型增加到7个，包括：linear, linear contact bond, linear parallel bond, hertz, hysteretic, smooth joint以及flat joint property tables。此外，除了原来的6个关键字activate, deactivate, distance f, install, model s, use-aperture，又新增了3个关键字：apply-groups，apply-properties和type。模型的工作机理是：(1)使用CMAT分配一个缺省的线性接触模型; (2)产生ball和fracture; (3) fracture contact-model命令在所有与断裂相交的接触中安装指定的模型。受接触的影响，球的尺寸影响着断裂的尺寸和分布。下图所示的分别为直径0.2,0.6和1.0时相应的断裂，可以看出，随着球半径的增加，平均的断裂长度在不断减小。

参考文档:

(1) Slip on a Fault

(2) Creation of a Synthetic Rock Mass (SRM) Specimen