岩土工程新技术

PFC2D 为二维颗粒流程序，用离散单元法模拟圆形颗粒介质运动与相互作用；将对象划分为数百个颗粒单元，本构性质由模型自动涌现。

与有限差/有限元/边界元等不同，离散元将含弱面的岩体视为块体集合，可发生平动、转动与分离。

颗粒流特点：圆形颗粒可直接模拟；块体可由多颗粒胶结；粒径可均匀或高斯分布；可模拟节理弱面；接触模型含线性弹簧、Hertz-Mindlin 与库仑滑块。

PFC2D 与 UDEC/3DEC 相比，圆形颗粒接触计算更简单，可有效模拟大变形与块体分离；与 FLAC 类似采用显式时步与局部非粘性阻尼。

阻尼优点：均匀运动时体 force 近零、阻尼系数无量纲且不随频率变；缺点：几何与边界条件设定较 FLAC/UDEC 困难，孔隙率与初应力常需试算。

PFC2D 用 ANSI C/C++ 编写，16 MB 内存可模拟约 3 万颗粒；内置 FISH 语言用于边界条件、伺服控制与参数分析。

颗粒流基本假设：颗粒刚体平动与转动；接触力由重叠量与接触模型决定；采用显式中心差分积分。

解题途径：建立几何与颗粒体系→赋予接触模型与微观参数→设置边界与加载→确定时步→运行并监测不平衡力→解释结果。

模拟前须补充几何、地质结构、材料性质、初始应力与外部荷载等简化模型数据，参数范围须考虑工程不确定性。

运行准备：确定合适时步；频繁保存状态便于重启；设监测点跟踪不平衡力与参数变化。