PFC5.0单轴压缩、级配，双轴压缩，巴西劈裂，直剪试验
在颗粒流模拟领域，PFC5.0的几个典型试验就像武林绝学里的基本功。咱们先从单轴压缩这个"扎马步"说起。做这个试验时，颗粒群的生成特别有意思——就像把不同大小的鹅卵石倒进玻璃罐，关键要让它们自然堆积形成初始孔隙结构。下面这段代码展示了如何用Fish语言生成级配颗粒：
;生成粒径分布 def gradation loop local n (1,5) d = 0.5 + 0.2*n ;粒径从0.7到1.5mm command ball distribute porosity 0.3 ... diameter [d] ... box -10 10 -10 10 -10 10 end_command end_loop end gradation这里用循环生成了五组不同粒径的颗粒，0.3的孔隙率参数不是随便定的——这是通过反复试算发现既能保证计算效率又能反映真实岩土特性的平衡值。注意直径增量设为0.2mm，过大的步长会导致级配曲线出现"台阶"，影响试样均匀性。
双轴压缩试验就像给颗粒系统戴上紧箍咒。在代码层面，需要同时控制x和y方向的围压。有经验的工程师都知道，下面这个伺服控制循环是成败关键：
def servoed_walls loop n (1,1000) wall servo on ;激活墙体伺服 solve age 50 ;分步计算 if math.abs(wall.force(1)-target_stress) < tol then break end_loop end servoed_walls这个控制逻辑看似简单，但有个魔鬼细节——循环次数上限设为1000次，其实是个经验值。在实际操作中，遇到高孔隙率试样时可能要调整到2000次以上才能收敛。另外，50个计算步长作为分段求解间隔，可以有效避免数值震荡。
巴西劈裂试验的代码实现堪称暴力美学。核心在于如何用两个刚性板模拟加载过程：
;生成加载板 wall generate id 1 plane (0,5) (0,-5) (0,0) wall generate id 2 plane (0,-5) (0,5) (0,10) ;位移控制加载 wall attr yvelocity 0.1 id 1 wall attr yvelocity -0.1 id 2这里y方向加载速度设置为0.1m/s看似偏大，其实在准静态分析中需要通过试算验证。有个小技巧是把速度设为负值实现反向加载，这比用正方向配合坐标变换更直观。监测裂纹发展时，建议同时跟踪接触力链和配位数变化，这样能更早发现破坏前兆。
直剪试验的难点在于边界条件处理。下面这段接触力监测代码经常被新手忽视：
fish define shear_force local sum_f = 0.0 contact.list c if c.pos(2) > shear_plane_y sum_f = sum_f + c.force.x endif end_loop global shear_force = sum_f end这个自定义函数实时统计剪切面上的水平接触力分量。特别注意用contact.list遍历所有接触比用ball.contact效率高40%左右，在处理十万级接触量时这个优化至关重要。当剪切位移达到试样尺寸的15%时，建议开始降低加载速度以避免数值发散。
这些试验看似基础，但就像做菜的火候掌控，参数设置的微妙差异会导致结果天差地别。有个项目曾因把伺服控制的容差参数tol从1e3调到1e4，结果峰值强度预测值偏差了18%。所以每次看到论文里漂亮的应力-应变曲线，都应该想到背后可能藏着上百次的参数调校。