离散元软件PFC入门指导讲述.ppt

PFC2D颗粒流软件培训中仿科技有限公司 主要内容 第一部分 PFC2D颗粒流程序简介 第二部分 有限差分法基础介绍 第三部分 离散元法基础介绍 第四部分 PFC2D的使用 第一部分 PFC2D颗粒流程序简介 1 、理论背景 2、颗粒流方法的基本假设 3、颗粒流方法的特点 4、可选特性 5、应用领域 6、求解步骤 作为离散元的一种，二维颗粒流程序(Particle Follow Code PFC2D)数值模拟新技术，其理论基础是Cundall [1979]提出的离散单元法，用于颗粒材料力学性态分析，如颗粒团粒体的稳定、变形及本构关系，专门用于模拟固体力学大变形问题。它通过圆形(或异型)离散单元来模拟颗粒介质的运动及其相互作用。由平面内的平动和转动运动方程来确定每一时刻颗粒的位置和速度。作为研究颗粒介质特性的一种工具，它采用有代表性的数百个至上万个颗粒单元，通过数值模拟实验可以得到颗粒介质本构模型。 PFC2D (Particle Follow Code 2 Dimension)即二维颗粒流程序，是通过离散单元方法来模拟圆形颗粒介质的运动及其相互作用。最初，这种方法是研究颗粒介质特性的一种工具，它采用数值方法将物体分为有代表性的数百个颗粒单元，期望利用这种局部的模拟结果来研究边值间题连续计算的本构模型。以下两种因素促使PFC2D方法产生变革与发展:(1)通过现场实验来得到颗粒介质本构模型相当困难:(2)随着微机功能的逐步增强，用颗粒模型模拟整个问题成为可能，一些本构特性可以在模型中自动形成。因此，PFC2D便成为用来模拟固体力学和颗粒流问题的一种有效手段。 2、颗粒流方法的基本假设 颗粒流方法在模拟过程中作了如下假设: 1)颗粒单元为刚性体; 2)接触发生在很小的范围内，即点接触; 3)接触特性为柔性接触，接触处允许有一定的“重叠”量; 4) “重叠”量的大小与接触力有关，与颗粒大小相比，“重叠”量很小; 5)接触处有特殊的连接强度; 6)颗粒单元为圆盘形(或球形)。 其中，颗粒为刚性体的假设，对于模拟介质运动为只沿相互接触面的表面发生的问题非常重要，比如象砂土或粮食这种颗粒组合体材料，利用这种假设在总体上来讲是比较恰当的，因为这种材料的变形是来自于颗粒刚性体间的滑动和转动以及接触面处的张开和闭锁，而不是来自于每个刚性颗粒本身的变形，对于这种特殊材料，没有必要采用非常精确的数值模型，来得到对材料特性的近似。 3、颗粒流方法的特点 PFC2D可以直接模拟圆形颗粒的运动和相互作用问题。颗料可以代表材料中的个别颗粒，例如砂粒，也可以代表粘结在一起的固体材料，例如混凝土或岩石。当粘结以渐进的方式破坏时，它能够破裂。粘结在一起的集合体可以是各向同性，也可以被分成一些离散的区域或块体。这类物理系统可以用处理角状块体的离散单元程序UDEC和3DEC来模拟。 PFC2D有三个优点： 第一、它有潜在的高效率。因为圆形物体间的接触探测比角状物体间的更简单。 第二、对可以模拟的位移大小实质上没有限制。 第三、由于它们是由粘结的粒子组成，块体可以破裂，不象UDEC和3DEC模拟的块体不能破裂。 用PFC2D模拟块体化系统的缺点是，块体的边界不是平的，用户必须接受不平的边界以换取PFC2D提供的优点。 PFC2D中几何特征、物理特性和解题条件的说明不如FLAC和UDEC程序那样直截了当。 例如用连续介质程序，创建网格、设置初始压力、设置固定或自由边界。在象PFC2D这样的颗粒程序中，由于没有唯一的方法在一个指定的空间内组合大量的粒子，粒子紧密结合的状态一般不能预先指定。必须跟踪类似于物体压实的过程，直到获得要求的孔隙率。 由于颗粒相对位置变化产生接触力，初始应力状态的确定与初始压密有关。由于边界不是由平面组成，边界条件的设定比连续介质程序更复杂。 当要求满足有实验室实际测试的模拟物体的力学特性时，出现了更大的困难。在某种程度上，这是一个反复试验的过程，因为目前还没有完善的理论可以根据微观特性来预见宏观特性。 然而，给出一些准则应该有助于模型与原型的匹配，如哪些因素对力学行为的某些方面产生影响，哪些将不产生影响。应该意识到，由于受现有知识的限制，这样的模拟很难。然而，用PFC2D进行试验，对固体力学，特别是对断裂力学和损伤力学，可以获得一些基本认识。 PFC2D能模拟任意大小圆形粒子集合体的动态力学行为。 粒子生成器根据粒子的指定分布规律自动概率地生成。粒子半径按均匀分布或按高斯分布规律分布。 初始孔隙度一般比较高，但通过控制