保姆级教程：用PFC模拟岩石巴西劈裂试验，从建模到结果分析全流程

离散元仿真实战：PFC岩石巴西劈裂试验全流程解析

岩土工程领域的数值模拟技术正经历从宏观连续体到微观离散介质的范式转移。作为离散元方法（DEM）的代表性工具，PFC（Particle Flow Code）通过颗粒集合体的力学行为模拟，为研究人员提供了观察岩石破裂过程的"数值显微镜"。本文将聚焦巴西劈裂试验这一经典岩石力学测试方法，通过PFC实现从模型构建、参数设定到结果分析的全流程实战演示。

1. 试验原理与模型构建基础

巴西劈裂试验最初由巴西学者提出，现已成为测定岩石抗拉强度的标准方法。其核心原理是通过径向压缩圆盘试样，在圆心区域产生均匀拉应力场。与传统单轴压缩试验相比，该试验能更准确反映岩石的抗拉特性——这一在岩体破坏中起关键作用的参数。

1.1 离散元模型初始化

在PFC中构建圆形试样时，需特别注意以下参数设置：

new def par sample_radius=0.4 rdmin=0.006 rdmax=0.009 poro=0.12 end @par

表：关键初始化参数说明

颗粒生成阶段采用ball distribute命令时，通过range annulus限定颗粒分布区域，确保形成规则的圆形试样边界：

ball distribute porosity @poro radius [rdmin] [rdmax] range annulus... center 0 0 radius 0 @sample_radius

提示：resolution参数控制边界墙体的离散精度，值越小圆形越光滑，但会相应增加计算量

1.2 材料参数与接触模型

岩石材料的力学特性通过线性接触模型定义：

cmat default model linear method deform ball property emod 100e6 kratio 1.5 fric 0.5

• emod：弹性模量(Pa)
• kratio：法向/切向刚度比
• fric：摩擦系数

初始平衡阶段使用calm命令消除系统初始不平衡力：

cycle 2000 calm 50

2. 伺服控制预压技术详解

预压阶段通过伺服控制机制实现试样均匀压缩，这是确保试验结果准确性的关键步骤。伺服控制的本质是通过反馈调节维持目标应力状态。

2.1 伺服算法实现

核心伺服函数包含三个主要部分：

def calStress sumForce=0 loop foreach ct contact.list("ball-facet") sumForce+=contact.force.normal(ct) endloop wsrr=sumForce/(2*math.pi*wlr) ; 环向应力计算 end

应力计算采用环向应力公式，其中wlr为实时更新的试样半径。速度调节算法通过应力偏差控制墙体运动速度：

def sevro_wall calStress if global.step>time_record then getg time_record=global.step+sevro_freq endif rvel=gr*(wsrr-trr) ; 速度计算 ... ; 速度分配 end

伺服控制关键参数

• servo_factor：伺服稳定性系数(0.1-1.0)
• sevro_freq：伺服调用频率(步数)

2.2 边界条件处理

墙体顶点速度分配采用矢量运算：

loop foreach vt wall.vertexlist(wp) direct=math.unit(wall.vertex.pos(vt)) vel_vector=direct*rvel wall.vertex.vel(vt)=vel_vector endloop

注意：顶点速度方向始终沿径向，确保试样均匀压缩

3. 胶结模型与加载设置

预压完成后需建立颗粒间胶结，模拟岩石内部结构。PFC中可通过平行粘结模型实现：

3.1 胶结参数设定

contact method bond gap 0.001 contact property pb_ten 1e6 pb_coh 1e6 pb_fa 30

胶结参数影响分析

• 抗拉强度(pb_ten)决定裂纹萌生阈值
• 内聚力(pb_coh)影响裂纹扩展路径
• 摩擦角(pb_fa)控制破坏后行为

3.2 加载系统配置

巴西劈裂采用上下加载板方式：

wall create id 1 vertices [-wlr] [wlr] [wlr] [wlr] wall create id 2 vertices [-wlr] [-wlr] [wlr] [-wlr]

位移控制加载通过设定加载板速度实现：

[strainRate=1e-1] wall attribute yvel [strainRate*wlr*2] range id 2 wall attribute yvel [-strainRate*wlr*2] range id 1

4. 结果分析与后处理

试验数据的实时采集通过FISH函数实现，主要监测两个关键指标：

4.1 力学响应计算

def computer_strain disp=wall.disp.y(wpup)-wall.disp.y(wpdown) stress=(wall.force.contact.y(wpup)-wall.force.contact.y(wpdown))/(math.pi*2*wlr) end

巴西劈裂强度公式$$ \sigma_t = \frac{2P_{max}}{\pi Dt} $$ 其中：

• $P_{max}$：峰值载荷
• $D$：试样直径
• $t$：试样厚度

4.2 破裂过程观测

PFC提供多种裂纹可视化方法：

1. 接触力链显示

   plot create view 1 plot item create ball colorby contactforce

2. 粘结断裂监测

   history bond break

3. 能量演化分析

   history energy

典型破坏模式表现为：

• 初始微裂纹在圆心处萌生
• 主裂纹沿加载轴向扩展
• 最终形成贯通破裂面

5. 常见问题排查指南

在实际操作中常遇到以下典型问题：

试样不均匀变形

• 检查伺服控制参数servo_factor
• 验证颗粒尺寸分布范围
• 确保边界条件对称性

过早破裂

• 提高胶结强度参数
• 降低加载速率strainRate
• 检查初始孔隙率设置

数值震荡

• 调整阻尼系数ball attribute damp
• 减小时间步set dt
• 增加伺服频率sevro_freq

6. 单轴压缩与巴西劈裂对比

两种试验在PFC实现中的主要差异：

表：试验配置对比

实现巴西劈裂时需要特别注意：

1. 圆形试样生成精度
2. 径向伺服控制算法
3. 加载板接触处理
4. 裂纹扩展方向判定

实际项目中，我们常通过调整resolution参数来平衡计算效率和边界精度。当研究各向异性岩石时，还需考虑颗粒排列方向性的影响。