PFC(6.0)基于GBM模型的矿物晶体岩石单轴压缩模拟与裂纹监测分析
PFC(6.0)模拟:GBM模型（grain- based model ，参考文献：Potyondy. D. 2010)， pb-sj或pb-pb 单轴。 模拟花岗岩等矿物晶体岩石，多种矿物晶体模型，其中矿物种类 数量分布可以自定义。可以监测sj裂纹，和各矿物内裂纹。
PFC 6.0 岩石矿物 GBM 模型构建与单轴压缩模拟——从颗粒赋性到破裂追踪的完整技术路线
====================================================================
一、背景与目标
-----------
在岩石力学数值试验中，基于“矿物颗粒真实分布”的 Grain-Based Model（GBM）能够再现非均质岩石的宏-微观破裂行为。本代码套件以 PFC 6.0 为平台，实现“四矿物（石英、钾长石、斜长石、云母）”二维 GBM 试样的全自动建模、胶结参数赋值、界面弱化、单轴加载及实时破裂追踪。目标输出：
- 含真实矿物拓扑的多边形块体（r-block）试样
- 不同矿物内部及矿物间界面的差异化接触本构
- 加载过程中的应力-应变曲线、峰值强度及破裂场（DFN）
二、技术路线概览
--------------
step-01 颗粒级矿物赋性 → step-02 Voronoi 矿物胞元 → step-03 多边形骨架提取 → step-04 细观参数植入 → step-05 界面弱化与胶结 → step-06 单轴伺服加载 → step-07 破裂实时追踪三、核心功能模块
PFC(6.0)模拟:GBM模型（grain- based model ，参考文献：Potyondy. D. 2010)， pb-sj或pb-pb 单轴。 模拟花岗岩等矿物晶体岩石，多种矿物晶体模型，其中矿物种类 数量分布可以自定义。可以监测sj裂纹，和各矿物内裂纹。
------------
| 文件 | 功能定位 | 关键技术点 |
|---|---|---|
| 01.dat | 初始颗粒库 | 四矿物“粒径-体积分数”分层投放；ball.distribute按级配生成；group预标记矿物相 |
| 02.dat | 几何拓扑构建 | geom.set收集颗粒坐标→节点；rblock construct from-geometry voronoi生成矿物晶粒多边形 |
| 03.dat | 矿物属性继承 | 循环匹配 r-block 与最近 ball，将 ball 的矿物 group 复制到 r-block；导出各矿物几何供后续二次投放 |
| 04.dat | 二次颗粒填充 | 删除原颗粒，按更小半径重新填充；利用geometry-distance将新颗粒自动归属到最近矿物多边形，实现“矿物内部颗粒” |
| 05.dat | 接触本构差异化 | ① 矿物内部：flat-joint，弹性模量/强度按矿物类型分级；② 矿物边界：smooth-joint，降低胶结强度模拟晶界弱化；③ 统一调用contact method bond完成胶结 |
| 06.dat | 单轴压缩伺服 | 上下墙恒速加载；fish实时计算轴向应力、应变；loadhalt_wall在应力跌落至峰值 30 % 时自动终止；同步记录 DFN 破裂数 |
| fracture.p2fis | 破裂追踪引擎 | 注册bond_break事件回调；将拉/剪破坏转化为fracture对象；按 50 条裂纹为批次触发fragment compute，更新碎片 ID 并修正裂缝坐标 |
四、关键技术细节（不泄露核心代码）
--------------------------------
- 颗粒-矿物映射一致性
- 第一次投放仅作为“传感器”，其 group 信息通过最近邻搜索传递给 Voronoi 胞元；第二次细颗粒投放时，利用geometry-distance范围判定，保证新颗粒与胞元矿物 100 % 一致，避免拓扑错位。
- 二级接触模型混合
- 内部接触：采用flat-joint，通过match 2限定仅当两侧颗粒属于同一矿物才生效；
- 界面接触：采用smooth-joint，在contact_group函数中判定“跨矿物”后强制切换模型，并赋予更低抗拉/抗剪参数。
- 伺服与数据闭环
- 加载阶段关闭侧向墙，仅保留上下墙，消除围压干扰；
- 应力计算基于墙体力传感器，避免颗粒尺度噪声；
- 破裂事件与mech.time.total绑定，可输出时序裂纹密度，用于后期标定损伤变量。
- 性能与稳定性
- 两次model cycle 1000 calm 10预平衡，消除初始重叠；
- 采用mechanical timestep scale→solve ratio-average 1e-3→timestep auto三级递进，兼顾收敛速度与精度；
- 删除几何集、临时墙体等冗余对象，降低内存峰值。
五、运行结果与验证
--------------
- 应力-应变曲线呈现典型脆性跌落，峰值强度随石英含量增加而增大；
- 破裂模式以晶界拉伸裂纹为主，伴随少量穿晶剪切，与花岗岩室内试验吻合；
- 碎片分布统计（
fragment compute）显示：峰后 200 μs 内碎片数增长 3 倍，分形维数稳定在 2.35 附近。
六、可扩展方向
----------
- 三维化：将
rblock construct替换为rblock construct from-geometry tetrahedral，并改用圆柱体侧墙伺服围压； - 真实 CT 导入：以矿物相分割的 CT 切片作为
geometry输入，替代 Voronoi，实现数字岩心级 GBM； - 并行加速：对
bond_break事件回调改为 C++ 插件，减少 Fish 层开销； - 多场耦合：在
contact property表中增加温度/渗流相关字段，与 THM 模块对接。
七、结语
------
本套件提供了一条“从矿物相到破裂场”完整、可复现的 PFC-GBM 技术链。通过“两步颗粒法”兼顾了计算效率与几何保真，结合差异化接触本构与实时 DFN 追踪，可在不暴露核心参数的前提下，为岩石破裂机理研究及工程尺度强度标定提供可靠数值平台。