PFC2D颗粒流软件中文说明Andy墨攻1高级材料主要内容第一部分 PFC2D颗粒流程序简介 第二部分 有限差分法基础介绍第三部分 离散元法基础介绍第四部分 PFC2D的使用2高级材料第一部分 PFC2D颗粒流程序简介1 、理论背景2、颗粒流方法的基本假设3、颗粒流方法的特点4、可选特性5、应用领域6、求解步骤3高级材料 作为离散元的一种，二维颗粒流程序(Particle Follow Code PFC2D)数值模拟新技术，其理论基础是Cundall 1979提出的离散单元法，用于颗粒材料力学性态分析，如颗粒团粒体的稳定、变形及本构关系，专门用于模拟固体力学大变形问题它通过圆形(或异型)离散单元来模拟颗粒介质的运动及其相互作用由平面内的平动和转动运动方程来确定每一时刻颗粒的位置和速度作为研究颗粒介质特性的一种工具，它采用有代表性的数百个至上万个颗粒单元，通过数值模拟实验可以得到颗粒介质本构模型1.理论背景理论背景4高级材料 PFC2D (Particle Follow Code 2 Dimension)即二维颗粒流程序，是通过离散单元方法来模拟圆形颗粒介质的运动及其相互作用最初，这种方法是研究颗粒介质特性的一种工具，它采用数值方法将物体分为有代表性的数百个颗粒单元，期望利用这种局部的模拟结果来研究边值间题连续计算的本构模型。
以下两种因素促使PFC2D方法产生变革与发展:(1)通过现场实验来得到颗粒介质本构模型相当困难:(2)随着微机功能的逐步增强，用颗粒模型模拟整个问题成为可能，一些本构特性可以在模型中自动形成因此，PFC2D便成为用来模拟固体力学和颗粒流问题的一种有效手段 5高级材料2、颗粒流方法的基本假设 颗粒流方法在模拟过程中作了如下假设:1)颗粒单元为刚性体;2)接触发生在很小的范围内，即点接触;3)接触特性为柔性接触，接触处允许有一定的“重叠”量;4) “重叠”量的大小与接触力有关，与颗粒大小相比，“重叠”量很小;5)接触处有特殊的连接强度;6)颗粒单元为圆盘形(或球形)6高级材料 其中，颗粒为刚性体的假设，对于模拟介质运动为只沿相互接触面的表面发生的问题非常重要，比如象砂土或粮食这种颗粒组合体材料，利用这种假设在总体上来讲是比较恰当的，因为这种材料的变形是来自于颗粒刚性体间的滑动和转动以及接触面处的张开和闭锁，而不是来自于每个刚性颗粒本身的变形，对于这种特殊材料，没有必要采用非常精确的数值模型，来得到对材料特性的近似7高级材料3、颗粒流方法的特点 PFC2D可以直接模拟圆形颗粒的运动和相互作用问题。
颗料可以代表材料中的个别颗粒，例如砂粒，也可以代表粘结在一起的固体材料，例如混凝土或岩石当粘结以渐进的方式破坏时，它能够破裂粘结在一起的集合体可以是各向同性，也可以被分成一些离散的区域或块体这类物理系统可以用处理角状块体的离散单元程序UDEC和3DEC来模拟8高级材料 PFC2D有三个优点： 第一、它有潜在的高效率因为圆形物体间的接触探测比角状物体间的更简单 第二、对可以模拟的位移大小实质上没有限制 第三、由于它们是由粘结的粒子组成，块体可以破裂，不象UDEC和3DEC模拟的块体不能破裂 用PFC2D模拟块体化系统的缺点是，块体的边界不是平的，用户必须接受不平的边界以换取PFC2D提供的优点9高级材料 PFC2D中几何特征、物理特性和解题条件的说明不如FLAC和UDEC程序那样直截了当 例如用连续介质程序，创建网格、设置初始压力、设置固定或自由边界在象PFC2D这样的颗粒程序中，由于没有唯一的方法在一个指定的空间内组合大量的粒子，粒子紧密结合的状态一般不能预先指定必须跟踪类似于物体压实的过程，直到获得要求的孔隙率 由于颗粒相对位置变化产生接触力，初始应力状态的确定与初始压密有关由于边界不是由平面组成，边界条件的设定比连续介质程序更复杂。
10高级材料 当要求满足有实验室实际测试的模拟物体的力学特性时，出现了更大的困难在某种程度上，这是一个反复试验的过程，因为目前还没有完善的理论可以根据微观特性来预见宏观特性 然而，给出一些准则应该有助于模型与原型的匹配，如哪些因素对力学行为的某些方面产生影响，哪些将不产生影响应该意识到，由于受现有知识的限制，这样的模拟很难然而，用PFC2D进行试验，对固体力学，特别是对断裂力学和损伤力学，可以获得一些基本认识 11高级材料 PFC2D能模拟任意大小圆形粒子集合体的动态力学行为 粒子生成器根据粒子的指定分布规律自动概率地生成粒子半径按均匀分布或按高斯分布规律分布 初始孔隙度一般比较高，但通过控制粒子半径的扩大可以获得密度压实在任何阶段任何因素都可以改变半径所以不需反复试验就可以获得指定孔隙度的压实状态 12高级材料 属性与各个粒子或接触有关，而不是与“类型号”有关 因此，可以指定属性和半径的连续变化梯度节理生成器”用来修改沿指定轨迹线的接触特性假定这些线叠加在颗粒集合体上用这种方法，模型可以被成组的弱面，如岩石节理切割 粒子颜色也是一种属性，用户可以指定各种标记方案13高级材料 PFC2D模型中为了保证数据长期不漂移，用双精度数据存储坐标和半径。
接触的相对位移直接根据坐标而不是位移增量计算接触性质由下列单元组成： 1）线性弹簧或简化的Hertz-Mindlin准则； 2）库仑滑块； 3）粘结类型：粘结接触可承受拉力，粘结存在有限的抗拉和抗剪强度 可设定两种类型的粘结，接触粘结和平行粘结这两种类型粘结对应两种可能的物理接触：接触粘结再现了作用在接触点一个很小区域上的附着作用；平行粘结再现了粒子接触后浇注其它材料的作用（如水泥灌浆）平行粘结中附加材料的有效刚度具有接触点的刚度 14高级材料 块体逻辑支持附属粒子组或块体的创建，促进了程序的推广普及块体内粒子可以任意程度的重叠，作为刚性体具有可变形边界的每一个块体，可作为一般形状的超级粒子通过指定墙的速度、混合的粒子速度、施加外力和重力来给系统加载扩展的FISH库”提供了在集合体内设置指定应力场或施加应力边界条件的函数时步计算是自动的，包括因为Hertz接触模型刚度变化的影响模拟过程中，根据每个粒子周围接触数目和瞬间刚度值，时步也在变化基于估计的粒子数，单元映射策略采用最佳的单元数目，自动调整单元的外部尺寸来适应粒子缺失和指定的新对象单元映射方案支持接触探测算法以保证求解时间随粒子数目线性增加，而不是二次方增加。
15高级材料 类似于FLAC，PFC提供了局部无粘性阻尼这种阻尼形式有以下优点： 1）对于匀速运动，体力接近于零，只有加速运动时才有阻尼； 2）阻尼系数是无因次的； 3）因阻尼系数不随频率变化，集合体中具有不同自然周期的区域被同等阻尼，采用同样的阻尼系数 PFC2D可以在半静态模式下运行以保证迅速收敛到静态解，或者在完全动态模式下运行 PFC2D包含功能强大的内嵌式程序语言FISH，允许用户定义新的变量和函数使数值模型适合用户的特殊需求例如，用户可以定义特殊材料的模型和性质、加载方式、实验条件的伺服控制、模拟的顺序以及绘图和打印用户定义的变量等 16高级材料4、可选特性 1）热学分析2）并行处理技术3）能写用户定义接触模型4）用户写C+程序的C+编程 17高级材料 热学选项用来模拟材料内热量的瞬间流动和热诱导位移和力的顺序发展热学模型可以独立运行或耦合到力学模型通过修改粒子半径和平行粘结承受的力，产生热应变来解释粒子和粘结材料的受热 用户定义的接触本构模型可以用C+语言来编写，并编译成动态链接库文件，一旦需要就可以加载 用户写的C+程序选项允许用户用C+语言写自己的程序，创建可执行的PFC2D个人版本。
这个选项可以用来代替FISH函数，大大提高运行的速度 并行处理技术允许将一个PFC2D模型分成几个部分，每个部分可以在单独的处理器上平行运行与一个PFC2D模型在一个处理器上运行相比，平行处理在内存容量和计算速度方面得到大大提高18高级材料5、应用领域 PFC2D既可解决静态问题也可解决动态问题，既可用于参数预测，也可用于在原始资料详细情况下的实际模拟PFC2D 模拟试验可以代替室内试验在岩石与土体中开挖问题的研究与设计方面，实测资料相对较少，关于初始应力、不连续性等问题也只能部分了解而在松散介质流动问题中，影响流动介质不规律分布的影响因素很难定量描述因此，应用PFC2D 初步研究影响整个系统的一些参数的特性，对整个系统的特性有所了解后，就可以方便地设计模型模拟整个过程 19高级材料 PFC2D可以模拟颗粒间的相互作用问题、大变形问题、断裂问题等，适用于以下领域：（1）在槽、管、料斗、筒仓中松散物体的流动问题； （2）矿山冒落法开采中的岩体断裂、坍塌、破碎和岩块的流动问题； （3）铸模中粉料的压实问题； （4）由粘结粒子组成物体的碰撞及其动态破坏； （5）梁结构的地震响应及垮塌； （6）颗粒材料的基本特性研究，如屈服、流动、体积变化等； （7）固体的基本特性研究，如累积破坏、断裂。
20高级材料6、求解步骤 1)定义模拟对象 根据模拟意图定义模型的详细程序，假如只对某一力学机制的不同解释作出判断时，可以建立一个比较粗略的模型，只要在模型中能体现要解释的机制即可，对所模拟问题影响不大的特性可以忽略21高级材料2)建立力学模型的基本概念 首先对分析对象在一定初始特性形成初步概念为此，应先提出一些问题，如系统是否将变为不稳定系统、问题变形的大小、主要力学特性是否非线性、是否需要定义介质的不连续性、系统边界是实际边界还是无限边界、系统结构有无对称性等 综合以上内容来描述模型的大致特征，包括颗粒单元的设计、接触类型的选择、边界条件的确定以及初始平衡状态的分析22高级材料 3)构造并运行简化模型 在建立实际工程模型之前，先构造并运行一系列简化的测试模型，可以提高解题效率通过这种前期简化模型的运行，可对力学系统的概念有更深入的了解，有时在分析简化模型的结果后(例如所选的接触类型是否有代表性、边界条件对模型结果的影响程度等)，还需将第二步加以修改23高级材料 4)补充模拟问题的数据资料 模拟实际工程问题需要大量简化模型运行的结果，对于地质力学来说包括: a)几何特性，如地下开挖酮室的形状、地形地貌、坝体形状、岩土结构等； b)地质构造位置，如断层、节理、层面等； c)材料特性，如弹/塑性、后破坏特性等； d)初始条件，如原位应力状态、孔隙压力、饱和度等； e)外荷载，如冲击荷载、开挖应力等。
因为一些实际工程性质的不确定性(特别是应力状态、变形和强度特性)，所以必须选择合理的参数研究范围第三步简化模型的运行有助于这项选择，从而为更进一步的试验提供资料24高级材料 5)模拟运行的进一步准备 a)合理确定每一时步所需时间，若运行时间过长，很难得到有意义的结论，所以应该考虑在多台计算机上同时运行 b)模型的运行状态应及时保存，以便在后续运行中调用其结果例如如果分析中有多次加卸荷过程，要能方便地退回到每一过程，并改变参数后可以继续运行 c)在程序中应设有足够的监控点(如参数变化 处、不平衡等)，对中间模拟结果随时作出比较分析，并分析颗粒流动状态25高级材料 6)运行计算模型 在模型正式运行之前先运行一些检验模型，然后暂停，根据一些特性参数的试验或理论计算结果来检查模拟结果是否合理，当确定模型运行正确无误时，连接所有的数据文件进行计算 7)解释结果 计算结果与实测结果进行分析比较图形应集中反应要分析的区域如应力集中区，各种计算结果应能方便地输出分析26高级材料第二部分 有限差分法基础介绍 连续介质三维快速拉格朗日有限差分计算方法( FLAC3D) 是近20年来逐步成熟完善起来的一种新型数值计算方法,它基于显式差分法来求解运动方程和动力方程,可模拟岩土或其他材料的三维力学行为。
其求解时首先将计算区域离散化,分成若干三维单元,单元之间由节点联结,节点受荷载作用后,其平衡方程(运动方程。