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 3clump的惯性参数要特别注意石英的转动惯量得设大些才不会乱滚。接触模型选线性接触带摩擦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以下除非你想让电脑风扇表演直升机起飞。