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 试样的全自动建模、胶结参数赋值、界面弱化、单轴加载及实时破裂追踪。目标输出：

1. 含真实矿物拓扑的多边形块体（r-block）试样
2. 不同矿物内部及矿物间界面的差异化接触本构
3. 加载过程中的应力-应变曲线、峰值强度及破裂场（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 并修正裂缝坐标

四、关键技术细节（不泄露核心代码）

--------------------------------

1. 颗粒-矿物映射一致性
   - 第一次投放仅作为“传感器”，其 group 信息通过最近邻搜索传递给 Voronoi 胞元；第二次细颗粒投放时，利用 geometry-distance 范围判定，保证新颗粒与胞元矿物 100 % 一致，避免拓扑错位。

1. 二级接触模型混合
   - 内部接触：采用 flat-joint，通过 match 2 限定仅当两侧颗粒属于同一矿物才生效；
   - 界面接触：采用 smooth-joint，在 contact_group 函数中判定“跨矿物”后强制切换模型，并赋予更低抗拉/抗剪参数。

1. 伺服与数据闭环
   - 加载阶段关闭侧向墙，仅保留上下墙，消除围压干扰；
   - 应力计算基于墙体力传感器，避免颗粒尺度噪声；
   - 破裂事件与 mech.time.total 绑定，可输出时序裂纹密度，用于后期标定损伤变量。

1. 性能与稳定性
   - 两次 model cycle 1000 calm 10 预平衡，消除初始重叠；
   - 采用 mechanical timestep scale → solve ratio-average 1e-3 → timestep auto 三级递进，兼顾收敛速度与精度；
   - 删除几何集、临时墙体等冗余对象，降低内存峰值。

五、运行结果与验证

--------------

• 应力-应变曲线呈现典型脆性跌落，峰值强度随石英含量增加而增大；
• 破裂模式以晶界拉伸裂纹为主，伴随少量穿晶剪切，与花岗岩室内试验吻合；
• 碎片分布统计（fragment compute）显示：峰后 200 μs 内碎片数增长 3 倍，分形维数稳定在 2.35 附近。

六、可扩展方向

----------

1. 三维化：将 rblock construct 替换为 rblock construct from-geometry tetrahedral，并改用圆柱体侧墙伺服围压；
2. 真实 CT 导入：以矿物相分割的 CT 切片作为 geometry 输入，替代 Voronoi，实现数字岩心级 GBM；
3. 并行加速：对 bond_break 事件回调改为 C++ 插件，减少 Fish 层开销；
4. 多场耦合：在 contact property 表中增加温度/渗流相关字段，与 THM 模块对接。

七、结语

------

本套件提供了一条“从矿物相到破裂场”完整、可复现的 PFC-GBM 技术链。通过“两步颗粒法”兼顾了计算效率与几何保真，结合差异化接触本构与实时 DFN 追踪，可在不暴露核心参数的前提下，为岩石破裂机理研究及工程尺度强度标定提供可靠数值平台。