PFC5.0代码：含三种矿物组成的岩石或类岩石材料GBM单轴压缩2d算例代码，仅供学习与提升

PFC5.0代码，含三种矿物组成的岩石或者类岩石材料，GBM，单轴压缩2d，算例代码仅供学习以及提升

打开PFC5.0的建模界面，突然想把花岗岩里的石英、长石、云母做成颗粒组合。先整点暴力的——直接拿球体颗粒拼成矿物晶粒，三种矿物用不同颜色标记，看着就像早餐麦片里的混合谷物。

PFC5.0代码，含三种矿物组成的岩石或者类岩石材料，GBM，单轴压缩2d，算例代码仅供学习以及提升

创建矿物组别时记得用group命令打标签：

group quartz range position-x 0 0.5 ; 左边区域是石英 group feldspar range position-x 0.5 1.0 group mica range position-y 0.8 1.0 ; 顶部云母层

这个范围选择可以玩出花样，比如用圆形区域划分矿物分布更符合真实岩石结构。颗粒替换操作是GBM建模的灵魂：

ball replace 1000 group quartz clump 1 ; 石英颗粒替换为clump1 ball replace 800 group feldspar clump 2 ball replace 200 group mica clump 3

clump的惯性参数要特别注意，石英的转动惯量得设大些才不会乱滚。接触模型选线性接触带摩擦：

contact cmat default model linear ... property kn 1e8 ks 1e8 fric 0.5

不同矿物间的接触参数可以搞差异化，比如云母与长石接触的摩擦系数调低到0.3，模拟层状结构的滑移。墙体加载时用伺服控制保持恒定应变率：

wall servo on velocity 1e-5 ; 上下墙以1e-5 m/s压缩

实时监测应力曲线的代码要嵌入measure对象：

measure create name 'stress' type stress measure create name 'strain' type axial

跑着跑着发现云母颗粒开始成片飞溅——这可不是bug，是真实的片状矿物破坏特征。用高帧率录像功能捕捉裂纹扩展：

image sequence interval 1000 ; 每1000步存一帧 plot bitmap size 1024 768 ; 高清画质记录脆性断裂

最后用Python脚本处理结果时发现个彩蛋：把应力-应变曲线的下降段做成动态条形图，居然和重金属乐队音量波形神似。完整代码包已扔在GitHub，记得把颗粒数量调小到5000以下，除非你想让电脑风扇表演直升机起飞。