打开PFC2D的操作界面，先别急着敲代码。咱们今天要搞的这个二维岩石单轴压缩模型，核心在于怎么让颗粒乖乖排好队再被压碎。直接上硬菜，看看这个模型的骨架结构

PFC2D5.0颗粒流『二维岩石单轴压缩』完整代码（附能量分析） 该代码包括： （1）完整代码及部分注释，可根据理解自行修改参数，点击运行即可得到结果，无需调试，可以直接使用，也可供参考学习； （2）应力应变曲线，完整能量分析，岩石裂隙分析，科研常规分析信息俱全。 代码部分详细信息： （1）岩石尺寸50mm×100mm，可修改； （2）平行粘结本构模型，模拟岩石胶结； （3）生成颗粒，平衡，伺服，卸载，加载，实验步骤完整，逻辑正确。 附图：应力应变曲线图，能量分析图，岩石裂隙统计图。

;=======模型初始化======= model new model large-strain on model domain extent -60 60 -60 60

这个域范围设置得比岩石尺寸大两倍是有讲究的——给颗粒变形留足空间。接下来生成岩石试件的关键操作：

;=======颗粒生成======= ball generate id 1 box -25 25 -50 50 radius 0.8 1.2 ball attribute density 2600 damp 0.7

看到没？0.8到1.2的随机半径让颗粒分布更接近真实岩石结构。密度2600kg/m³对应典型花岗岩参数，阻尼系数0.7能有效控制震荡。这时候运行代码会看到颗粒像撒豆子似的随机分布，别慌，后面有伺服控制来收拾残局。

粘结设置是重头戏：

;=======平行粘结======= contact cmat default model linearpbond property pb_ten 1e8 pb_coh 1e8 pb_fa 30

这里的1e8帕相当于100MPa抗拉强度，30度摩擦角是典型值。有意思的是，如果把这个抗拉强度改成5e7，模拟的岩石会像苏打饼干一样脆，各位可以自己试试。

加载阶段的操作看似简单实则暗藏玄机：

;=======加载控制======= wall id 100 name top wall generate plane 0 50 1 50 wall id 101 name bottom wall generate plane 0 -50 1 -50 wall attr yvelocity -0.1 target y -50 ; 加载速率0.1m/s

这个加载速率要是设得太猛，颗粒直接飞散就像爆米花。建议新手先调成-0.01慢慢观察破坏过程。当看到应力曲线出现断崖式下跌，说明岩石大哥正式躺平了。

能量分析模块才是科研狗的最爱：

;=======能量追踪======= history mech energy history unbalanced

跑完模拟直接导出energy.txt，用Python画个能量演化曲线，三根线分别对应应变能、动能和耗散能。记得看动能那个小突起——那是破裂开始的信号，比肉眼观察裂纹早得多。

PFC2D5.0颗粒流『二维岩石单轴压缩』完整代码（附能量分析） 该代码包括： （1）完整代码及部分注释，可根据理解自行修改参数，点击运行即可得到结果，无需调试，可以直接使用，也可供参考学习； （2）应力应变曲线，完整能量分析，岩石裂隙分析，科研常规分析信息俱全。 代码部分详细信息： （1）岩石尺寸50mm×100mm，可修改； （2）平行粘结本构模型，模拟岩石胶结； （3）生成颗粒，平衡，伺服，卸载，加载，实验步骤完整，逻辑正确。 附图：应力应变曲线图，能量分析图，岩石裂隙统计图。

裂隙统计有骚操作：

;=======裂纹计数器======= fish define crack_count local cnt = 0 loop foreach local c contact.list if contact.prop(c,"pb_ten") == 0.0 cnt = cnt +1 endif endloop crack_count = cnt end

这个fish函数实时统计粘结断裂数量，配合history记录，能精确捕捉每个破裂时刻。把裂纹数和应力曲线叠在一起，马上就能看出哪个应变值对应主破裂发生。

参数调整方面，墙裂建议把加载速率从0.1调到0.01试试，会发现应力-应变曲线的峰后跌落更平缓——这模拟的是不同应变率下的脆-延转变。有次手滑把颗粒数翻倍，结果显卡哼唧了半小时，所以模型尺寸别随便放大。

最后吐槽下PFC的能量平衡校验：当unbalanced值保持在1e-5以下才能算可靠结果。有次忘了设阻尼，看着unbalanced像窜天猴似的飙升，整个模型直接表演粒子大逃杀，场面堪比宇宙大爆炸。

完整代码跑完会自动输出三张关键图：应力应变曲线揭示强度特性，能量图展现破坏机理，裂纹统计图堪比岩石的CT扫描。把这些分析组合起来，妥妥能水一篇二区论文——当然，我是说严谨的科研成果。