PFC GBM岩石矿物多组分模型构建与力学模拟分析

PFC(6.0)模拟:GBM模型（grain- based model ，参考文献：Potyondy. D. 2010)， pb-sj或pb-pb 单轴。 模拟花岗岩等矿物晶体岩石，多种矿物晶体模型，其中矿物种类 数量分布可以自定义。可以监测sj裂纹，和各矿物内裂纹。

颗粒流代码（PFC）是一种基于离散元方法的数值模拟工具，广泛应用于岩石力学和材料科学领域。本文介绍的代码实现了一个基于晶粒基模型（GBM）的岩石矿物多组分模拟，能够模拟由多种矿物组成的岩石材料的力学行为和破坏过程。

模型概述

该GBM模型模拟了包含四种不同矿物的岩石材料：石英（shiying）、长石（jiachangshi）、斜长石（xiechangshi）和云母（yunmu）。每种矿物具有不同的物理特性和力学参数，通过分层建模的方式构建了一个真实的岩石微观结构模型。

模型构建流程

第一阶段：初始颗粒生成与平衡

模型首先在指定区域内生成四种不同粒径分布的矿物颗粒：

ball distribute porosity 0.05 resolution 1.0 number-bins 4 bin 1 radius 3e-3 2.5e-3 volume-fraction 0.15 group 'shiying' bin 2 radius 2.5e-3 2e-3 volume-fraction 0.20 group 'jiachangshi' bin 3 radius 2e-3 1.5e-3 volume-fraction 0.60 group 'xiechangshi' bin 4 radius 1.5e-3 1e-3 volume-fraction 0.05 group 'yunmu'

这一阶段设置了颗粒的基本物理属性（密度2500 kg/m³，阻尼系数0.7），并通过多次循环计算使系统达到力学平衡状态。

第二阶段：Voronoi结构生成

基于初始颗粒分布，代码通过Voronoi算法生成多边形晶粒结构：

geom set 'rock' rblock construct from-geometry 'rock' voronoi

这一步骤将离散的颗粒转化为连续的晶粒结构，为后续的矿物组分分配奠定基础。每个Voronoi单元根据其内部颗粒的矿物类型被标记为相应的矿物组分。

第三阶段：几何结构清理与重建

在完成矿物组分分配后，系统清理临时几何结构并重新生成更精细的颗粒体系。这一阶段包括：

• 删除临时几何集合
• 重新生成边界墙
• 以更低孔隙率（0.02）生成更小粒径的颗粒体系
• 根据几何距离将颗粒重新归类到相应的矿物组

第四阶段：接触模型与力学参数设置

这是模型构建中最关键的部分，为不同矿物和矿物界面设置了复杂的接触力学行为：

矿物内部接触（晶内接触）

对于同种矿物颗粒间的接触，采用Flatjoint模型模拟晶内力学行为：

contact model flatjoint range group 'shiying' match 2 contact method bond gap 0.5e-4 range group 'shiying' match 2 contact method deform emod @emod_a krat @emod_a range group 'shiying' match 2

矿物界面接触（晶界接触）

对于不同矿物颗粒间的接触，采用Smoothjoint模型模拟晶界行为：

contact model smoothjoint range group 'jiemian' contact property sj_kn 1e8 sj_kn 1e8 sj_fa 30 sj_ten 1.0e6 sj_coh 1.0e6 sj_fric 0.1 range group 'jiemian'

矿物力学参数

四种矿物被赋予了不同的力学特性：

• 石英：弹性模量20.0 GPa，抗拉强度10.0 MPa，内聚力10.0 MPa
• 长石：弹性模量15.0 GPa，抗拉强度8.0 MPa，内聚力8.0 MPa
• 斜长石：弹性模量10.0 GPa，抗拉强度6.0 MPa，内聚力6.0 MPa
• 云母：弹性模量5.0 GPa，抗拉强度4.0 MPa，内聚力4.0 MPa

这种参数设置反映了真实岩石中不同矿物的强度特性，石英通常是最坚硬的矿物，而云母则相对较弱。

第五阶段：单轴压缩试验模拟

模型设置了标准的单轴压缩试验环境：

wall attribute velocity-y -5 range id 3 wall attribute velocity-y 5 range id 1

通过上下墙体的相对运动对试样施加轴向荷载，同时监测系统的力学响应：

• 记录轴向应力和应变历史
• 跟踪裂纹的产生和发展
• 设置停止条件（应力峰值下降至30%时停止计算）

裂纹追踪与断裂分析

模型集成了先进的裂纹追踪系统，能够记录和分析材料的破坏过程：

fish define add_crack(entries) local contact = entries(1) local mode = entries(2) // 裂纹创建和分类逻辑 end

该系统能够：

• 区分拉伸破坏和剪切破坏
• 记录裂纹的位置、方向和尺寸
• 统计裂纹数量和发展趋势
• 可视化裂纹分布模式

技术特点与优势

1. 多组分真实性：模拟了四种不同矿物，更接近真实岩石的矿物组成

1. 微观结构准确性：通过Voronoi算法生成合理的晶粒结构

1. 力学行为差异化：不同矿物和界面具有不同的力学参数

1. 破坏过程可视化：完整的裂纹追踪和断裂分析功能

1. 计算稳定性：多阶段平衡过程确保模型在加载前的稳定性

应用价值

该GBM模型为岩石力学研究提供了强大的数值实验平台，可用于：

• 研究多矿物岩石的宏观力学行为与微观结构关系
• 分析裂纹在不同矿物和界面间的扩展规律
• 预测岩石材料的强度和破坏模式
• 优化岩石工程设计和风险评估

这种基于微观结构的建模方法能够更准确地反映真实岩石材料的力学行为，为深部资源开采、地质灾害防治和地下工程建设提供重要的理论支撑和技术手段。

PFC(6.0)模拟:GBM模型（grain- based model ，参考文献：Potyondy. D. 2010)， pb-sj或pb-pb 单轴。 模拟花岗岩等矿物晶体岩石，多种矿物晶体模型，其中矿物种类 数量分布可以自定义。可以监测sj裂纹，和各矿物内裂纹。