PFC 离散元数值模拟仿真 技术与应用

近几年，随着计算能力的提高和算法的优化，离散元仿真技术得到了快速发展，并在学术界产生了大量研究成果。在 PFC 离散元计算中无需给定材料的宏观本构关系和对应的参数，这些传统的参数和力学特性在程序中可以自动得到。据调查，运用 PFC 离散元仿真技术工具近几年发表的论文主要集中在以下几个方面：

1. 材料特性的模拟：研究不同材料（如岩石、土壤、混凝土等）在各种加载条件下的力学行为，包括弹性、塑性、断裂和破坏。

2. 颗粒形状和级配的影响：探讨颗粒的形状、大小分布（级配）对颗粒系统整体性能的影响，以及如何通过调整这些参数来模拟实际材料。

3. 接触模型和参数标定：开发和优化颗粒间的接触模型，以及如何准确标定模型参数以反映实际材料特性。

4. 边界条件和加载方式：研究不同的边界条件和加载方式对仿真结果的影响，以及如何设置这些条件以模拟实际工程问题。

5. 多尺度和多物理场耦合：探索 PFC 与其他数值方法（如有限元方法）的耦合，以实现多尺度或多物理场（如热-力耦合）的仿真。

6. 工程应用案例：在实际工程问题中的应用，如隧道开挖、边坡稳定性分析、地下结构设计等。

7. 环境因素和可持续性：探讨环境因素（如湿度、温度变化）对材料行为的影响，以及 PFC 在评估材料可持续性方面的应用。

PFC 离散元数值模拟仿真技术与应用

理论基础及 PFC 入门
1 岩土工程数值模拟方法概述

1.1 基于网格的模拟方法：

有限元、有限差分、大变形处理 CEL、ALE、XFEM

1.2 基于点的模拟方法：

离散单元法 DEM、光滑粒子流方法 SPH、物质点法 MPM

1.3 基于块体的模拟方法

2 离散元与 PFC 软件操作

2.1 离散元的基本原理（计算原理、宏观参量与微观参量的关系）

2.2 PFC 软件界面操作

2.3 文件系统

2.4 显示控制

2.5 帮助文档的使用

FISH 、N PYTHON 语言及COMMAND 命令
3 PFC 软件的计算控制方法

3.1 PFC 计算控制的语言逻辑

3.2 FISH 语言（基本语法、函数定义与调用、创建模型、控制模拟过程、

处理模拟结果、FISH Callback 操作等）

3.3 COMMAND 命令（命令结构、创建模型、状态监测与绘图、控制模拟

过程、求解控制、状态查询、与 FISH 语言的混合使用等）

3.4 PYTHON 语言（基本语法、Numpy 库的使用、接口的使用等）

离散元模拟方法
4 离散元模拟方法

4.1 离散元数值试样的生成方法

4.1.1 单元试样模型生成方法

4.1.2 边值问题（场地）模型生成方法

4.1.3 连续—非连续耦合模型生成方法

4.1.4 复杂颗粒形状的模拟方法（Rblock 方法、Clump 方法）

4.2 接触模型选择与参数标定

4.2.1 离散元接触模型的选择原则—12 个内置模型

4.2.2 接触模型参数的标定方法与参数意义—以胶结颗粒材料（岩石、胶结砂土等）为例，讲授参数标定步骤

4.3 其他问题

4.3.1 模型边界条件施加方法（达到初始平衡状态、开挖类模拟、填筑类模拟、加载类模拟、周期性边界、应力伺服）

4.3.2 各种阻尼的选择（粘滞阻尼、局部阻尼、滞回接触模型）

4.3.3 时步与时步缩放（静力、动力问题时步及相关命令）

4.3.4 试样尺寸、颗粒数量、级配选择

4.3.5 并行计算

土体单元试验模拟
5 土体单元试验模拟方法

5.1 常规三轴剪切试验模拟(命令流+FISH)

5.1.1 建模方法与注意事项

5.1.2 模拟结果分析

5.1.3 模拟结果可视化

5.2 真三轴剪切模拟(命令流+FISH)

5.2.1 真三轴加载路径的模拟

5.2.2 真三轴强度准则

5.2.3 微观结构演变过程

5.3 不排水三轴剪切模拟(命令流+FISH)

5.4 循环三轴剪切的模拟(命令流+FISH)

5.5 颗粒破碎过程模拟(命令流+FISH)

5.6 岩石（胶结颗粒）材料的剪切过程模拟

5.7 离散元模拟与弹塑性本构模型

工程实例分析
6 工程实例分析

6.1 活动门试验模拟(命令流+FISH)

6.1.1 试样级配控制

6.1.2 应力状态控制

6.1.3 孔隙比的控制

6.1.4 活动门加载的实现

6.2 盾构隧道掌子面稳定性(命令流+FISH)

6.2.1 主动失稳模式

6.2.2 被动失稳模式

6.3 节理岩体中的硐室开挖稳定性(命令流+FISH)

6.3.1 节理裂隙岩体的生成

6.3.2 初始应力状态控制

6.3.3 开挖模拟

PFC3D 与 FLAC3D 耦合模拟与分析
7 离散 — 连续域耦合模拟

7.1 离散—连续耦合模拟方法

 与 FLAC3D 中一维结构单元耦合

 与 FLAC3D 中二维壳结构单元或三维实体单元的面的耦合

 与 FLAC3D 中三维实体单元的耦合（实例）

7.2 离散—连续域参数匹配

7.3 基于离散—连续域耦合的三轴剪切试验模拟(命令流+FISH)

实例 操作 ：二维壳结构单元耦合 （壳单元模拟橡胶膜-创建耦合墙-施加应力边界等向压缩-剪切模拟）

7.4 基于离散—连续域耦合的地基承载力分析(命令流+FISH)

实例操作 ： 基于 n Punch indentation 案例的修改与实现

PFC- -CFD 耦合模拟与分析
8 流固耦合分析

8.1 颗粒与流体相互作用理论（CFD 模块概况、体积平均粗网格法、颗粒与流体相互作用计算）

8.2 流固耦合框架

 CFD 网格、流体域边界设置、网格导入、网格流体参数设置

 孔隙率计算

 耦合时间间隔、耦合时步、网格与颗粒尺寸

 耦合步骤

8.3 实例操作分析(命令流+FISH)

8.3.1 单向耦合

8.3.2 孔隙介质中 Darcy 流模拟（Fipy 应用）

8.3.3 与 FLAC3D 的渗流耦合模拟