PFC2D(ParticleFollowCode2Dimension)即二维颗粒流程序，是通过离散单元方法来模拟圆形颗粒介质的运动及其相互作用。最初，这种方法是研究颗粒介质特性的一种工具，它采用数值方法将物体分为有代表性的数百个颗粒单元，期望利用这种局部的模拟结果来研究边值间题连续计算的本构模型。以下两种因素促使PFC2D方法产生变革与发展:(1)通过现场实验来得到颗粒介质本构模型相当困难:(2)随着微机功能的逐步增强，用颗粒模型模拟整个问题成为可能，一些本构特性可以在模型中自动形成。因此，PFC2D便成为用来模拟固体力学和颗粒流问题的一种有效手段。颗粒流方法的基本假设：颗粒流方法在模拟过程中作了如下假设:1)颗粒单元为刚性体;2)接触发生在很小的范围内，即点接触;3)接触特性为柔性接触，接触处允许有一定的“重叠”量;4)“重叠”量的大小与接触力有关，与颗粒大小相比，“重叠”量很小;5)接触处有特殊的连接强度;6)颗粒单元为圆盘形(或球形)。3、颗粒流方法的特点：PFC2D可以直接模拟圆形颗粒的运动和相互作用问题。颗料可以代表材料中的个别颗粒，例如砂粒，也可以代表粘结在一起的固体材料，例如混凝土或岩石。当粘结以渐进的方式破坏时，它能够破裂。粘结在一起的集合体可以是各向同性，也可以被分成一些离散的区域或块体。这类物理系统可以用处理角状块体的离散单元程序UDEC和3DEC来模拟。PFC2D有三个优点：第一、它有潜在的高效率。因为圆形物体间的接触探测比角状物体间的更简单。第二、对可以模拟的位移大小实质上没有限制。第三、由于它们是由粘结的粒子组成，块体可以破裂，不象UDEC和3DEC模拟的块体不能破裂。用PFC2D模拟块体化系统的缺点是，块体的边界不是平的，用户必须接受不平的边界以换取PFC2D提供的优点。PFC2D能模拟任意大小圆形粒子集合体的动态力学行为。粒子生成器根据粒子的指定分布规律自动概率地生成。粒子半径按均匀分布或按高斯分布规律分布。初始孔隙度一般比较高，但通过控制粒子半径的扩大可以获得密度压实。在任何阶段任何因素都可以改变半径。所以不需反复试验就可以获得指定孔隙度的压实状态。属性与各个粒子或接触有关，而不是与“类型号”有关。因此，可以指定属性和半径的连续变化梯度。“节理生成器”用来修改沿指定轨迹线的接触特性。假定这些线叠加在颗粒集合体上。用这种方法，模型可以被成组的弱面，如岩石节理切割。粒子颜色也是一种属性，用户可以指定各种标记方案。PFC2D模型中为了保证数据长期不漂移，用双精度数据存储坐标和半径。接触的相对位移直接根据坐标而不是位移增量计算。接触性质由下列单元组成：1）线性弹簧或简化的Hertz-Mindlin准则；2）库仑滑块；3）粘结类型：粘结接触可承受拉力，粘结存在有限的抗拉和抗剪强度。可设定两种类型的粘结，接触粘结和平行粘结。这两种类型粘结对应两种可能的物理接触：①接触粘结再现了作用在接触点一个很小区域上的附着作用；②平行粘结再现了粒子接触后浇注其它材料的作用（如水泥灌浆）。平行粘结中附加材料的有效刚度具有接触点的刚度。块体逻辑支持附属粒子组或块体的创建，促进了程序的推广普及。块体内粒子可以任意程度的重叠，作为刚性体具有可变形边界的每一个块体，可作为一般形状的超级粒子。通过指定墙的速度、混合的粒子速度、施加外力和重力来给系统加载。“扩展的FISH库”提供了在集合体内设置指定应力场或施加应力边界条件的函数。时步计算是自动的，包括因为Hertz接触模型刚度变化的影响。模拟过程中，根据每个粒子周围接触数目和瞬间刚度值，时步也在变化。基于估计的粒子数，单元映射策略采用最佳的单元数目，自动调整单元的外部尺寸来适应粒子缺失和指定的新对象。单元映射方案支持接触探测算法以保证求解时间随粒子数目线性增加，而不是二次方增加。类似于FLAC，PFC提供了局部无粘性阻尼。这种阻尼形式有以下优点：1）对于匀速运动，体力接近于零，只有加速运动时才有阻尼；2）阻尼系数是无因次的；3）因阻尼系数不随频率变化，集合体中具有不同自然周期的区域被同等阻尼，采用同样的阻尼系数。PFC2D可以在半静态模式下运行以保证迅速收敛到静态解，或者在完全动态模式下运行。PFC2D包含功能强大的内嵌式程序语言FISH，允许用户定义新的变量和函数使数值模型适合用户的特殊需求。例如，用户可以定义特殊材料的模型和性质、加载方式、实验条件的伺服控制、模拟的顺序以及绘图和打印用户定义的变量等。FISH语言简介FISH是一种内置于Itasca软件内的编程语言，使用FISH用户可以定义新的变量和函数，从而使得这些函数被用来扩展Itasca软件的用法或者增加用户自定义的特性。例如，用户可以绘制(PLOT