PFC直剪试验：当岩石遇上浆液的暴力美学

原创于 2026-02-25 18:15:00 发布 · 174 阅读

4 ·

本内容遵循CC 4.0 BY-SA版权协议

GEO检测

PFC案例11，浆-岩充填结构面二维直剪，包含代码源文件、代码解释、曲线分析

实验室的直剪仪见过吧？今天咱们玩点更刺激的——用代码在虚拟世界复现岩体结构面的剪切破坏。打开案例包，先瞄一眼核心参数：法向应力2MPa，剪切速率0.1m/s，胶结强度30MPa...嗯，这配置是要搞事情啊。

先看模型搭建部分，这波操作堪称暴力美学：

; 剪切盒搭建
wall generate box ...
wall attribute density 2600 young 1e9 

; 颗粒生成
ball distribute ... radius 0.008 0.012
ball attribute density 2600 young 3e9 

; 胶结设置
contact method bond
contact property ... pb_ten 3e7 pb_coh 3e7

这里有个魔鬼细节——胶结参数设置的单位是帕斯卡不是兆帕！当年在这栽过跟头的我，看到数值3e7时DNA动了。颗粒young模量设3e9相当于花岗岩量级，wall用1e9模拟较软的围岩，这种刚度差异会让接触面更易产生滑动。

伺服控制这段代码最有意思：

[servoFlag = 0]
while math.abs(axz - target_stress)/target_stress > 0.02
    wall.force(..., component.z)
    wall.vel(..., component.z) *= 0.9
end

这段看似简单的循环藏着玄机：当实际应力与目标偏差超过2%时，自动调整墙体速度。就像用脚踩弹簧，踩轻了加点力，踩重了收点劲，直到压力稳定在目标值。调试时记得给速度乘衰减系数0.9，否则系统会像喝了假酒一样震荡。

PFC案例11，浆-岩充填结构面二维直剪，包含代码源文件、代码解释、曲线分析

剪切过程的数据记录堪称灵魂画手：

fish define record_data
    global disp = wall.disp(..., component.x)
    global shear_force = wall.force(..., component.x)
    file.write(disp, shear_force)
end
[extract.step = 1000]

每1000步记录一次位移和剪力，这采样频率就像高速摄像机，既不会漏掉关键帧，也不会被海量数据撑爆内存。偷偷说，把step改到500能看到更丝滑的曲线，当然得做好等咖啡的觉悟。