这次咱们来拆解PFC二维浆岩直剪案例。这个案例有意思的地方在于它展示了颗粒材料与刚性墙体接触面的剪切行为，咱们边看代码边分析剪切曲线的门道
PFC案例9，浆-岩二维直剪，包含代码源文件、代码解释、曲线分析
先看模型搭建的关键代码段：
;生成浆体颗粒 ball distribute ... ;创建上下剪切盒 wall generate id 1 vertices 0 0 1 0 1 1 0 1 wall generate id 2 vertices 0 0.2 1 0.2 ;设置接触参数 contact cmat ...这里通过wall命令生成两个墙体，上墙体保持固定，下墙体施加水平速度模拟剪切。颗粒生成时要注意粒径分布，太均匀的话剪切曲线会不真实。我一般习惯用高斯分布，参数设置成dmin=0.1e-3, dmax=0.3e-3比较合适。
剪切过程的控制逻辑是重点：
[shear_velocity = 0.1] while mechanical age < 10 wall attr velocity-x [shear_velocity] id 2 step 500 ... fish history shear_force = wall.force.x(id2) fish history displacement = wall.pos.x(id2) end这里采用位移控制法，每次循环移动下墙体并记录剪切力。注意速度参数shear_velocity别设太大，建议分阶段调整——初始阶段0.01，峰值后可以提到0.1，这样能更好捕捉软化段。
跑完模拟后咱们来看特征曲线（配图示意）。典型曲线会有四个阶段：
- 弹性段：力随位移线性增长
- 屈服段：曲线斜率开始减小
- 峰值段：达到最大剪切力
- 残余段：力值震荡趋于稳定
有个坑要注意：当剪切盒间隙过小时会出现颗粒被压碎的情况。这时候在contact cmat里调整kn/ks比例，把颗粒刚度调高些，或者在初始化时增大wall的gap值。
PFC案例9，浆-岩二维直剪，包含代码源文件、代码解释、曲线分析
裂纹发展监测可以用这个技巧：
fish define crack_count local cnt = 0 loop foreach bp ball.list if ball.contact.count(bp) < 2 cnt = cnt +1 endif endloop crack_count = cnt end这个简易计数器通过统计接触数小于2的颗粒数量，间接反映剪切带形成情况。实际跑起来会发现裂纹数量在峰值后突然增加，和应力跌落点基本对应。
最后给个实用建议：想得到漂亮的剪切曲线，别忘了在剪切前做充分的平衡计算。我习惯在施加剪切前先跑2000步静力松弛，这样初始接触力稳定，曲线毛刺会少很多。参数敏感性方面，摩擦系数对残余强度影响最直接，把颗粒摩擦角从30度调到35度，残余段力值能提高约18%。
下次可以试试在剪切面预置弱面，观察剪切路径的变化，那又是另一个有意思的故事了。代码文件已打包，需要调试的可以直接改参数里的标记处。