注：今天偶然间见到颗粒流的概念，以前一直不了解，今天查了查，贴在这里，以备以后可以温故知新。本文内容源自浙江大学罗永先生的博士论文，使得吾辈能花较少的时间看到广博的知识，在此特向其表示感谢！岩土工程数值计算总体上可以分为两大类:一类是基于连续介质力学理论的方法，如有限元法(FEM)和快速拉格朗日法 (FLAC(1tasea，2002))等;另一类是不连续介质力学的方法，如离散元法 UDEC(1tasca，2000)、3DEC(Itasea，1998)、PFC(Itasea，2002)和块体理论DDA(石根华，1988)等。离散元方法按其用途又可以分为宏观离散元方法和细观离散元方法，前者主要针对解决规模相对较大的不连续面，如断层节理结构与基础之间的结合面等引起的问题(UDEC，3DEC)，后者则着重于数目众多具有不连续特性的接触面或点，如破碎岩体中的破裂面、砂土中的接触面(点)和材料中颗粒之间的接触面(点)等。PFC(Particle Flow Code)是在著名学者 Peter Cundall主持下采用细观离散元理论(又称为粒子流理论)开发的一种数值计算平台，可以广泛地应用于研究细观结构控制问题。目前，PFC在世界上的应用并不广泛，成果报道也主要集中在PFC国际会议论文集中。
颗粒流PFCZD (Particle Flow Code in 2 Dimensions)平台数值模拟单元有两种:颗粒圆筒和颗粒(disc or particle)，主要用于平面应力和平面应变的特殊情况;颗粒流PFC3D(Particle Flow Code in 3 Dimensions)的数值模拟单元是三维球体颗粒(granular)，主要用于三维受力分析。Cundall(2002)博士认为PFC在描述岩土体介质特殊特性方面有着其他常用数值方法不可比拟的优势，主要表现在如下方面:(l)能自动模拟介质基本特性随应力环境的变化;(2)能实现岩土体对历史应力一应变记忆特性的模拟(屈服面变化Kaiser效等);(3)反映剪胀及其对历史应力等的依赖性;(4)自动反映介质的连续非线行应力一应变关系屈服强度和此后的应变软化或硬化过程;(5)能描述循环加载条件下的滞后效应;(6)描述中间应力增大时介质特性的脆性一塑性转化;(7)能考虑增量刚度对中间应力和应力历史的依赖性;(8)能反映应力一应变路径引起的刚度和强度的各向异性问题;(9)描述了强度包线的非线性特征;(10)介质材料微裂缝的自然产生过程;PFC的基本特点粒子流属于不连续介质力学的一种方法，这里的粒子并不直接与介质中是否存在颗粒状物质有关，只是用来描述介质特性的一种方式。比如PFC既可以用来描述具有颗粒物质的粗粒花岗岩一类的介质，也可以用来研究非晶质材料的特性。粒子流模型主要反映了颗粒集合体的力学行为，在粒子流模型内，离散的粒子被认为是刚性的，粒子之间的接触方式和力学特征可以不同，但符合基本的牛顿运动定律，即当粒子间的静力平衡被破坏时，粒子产生运动。粒子间的接触方式和接触强度是最关键的一个环节，决定了粒子集合体即介质的基本力学特性以及具体的承受张剪压力和保持静力平衡的基本能力。粒子集合体的各种复杂力学特性，比如其非线性特征和破坏特征都是通过粒子间的基本状态体现出来的。粒子间的接触出现破坏标志着粒子集合体由线性到非线性力学特征的开始。因此，在利用粒子流方法求解有关问题时不需要定义介质的本构关系。介质在复杂应力状态下的应力一应变关系将由其内部粒子间接触变化(如裂纹扩展)的情况所决定。粒子间的接触关系可以处理成非豁结和钻结两种方式。当豁结强度在没有达到破坏时，勃结介质对粒子集合体基本特性的影响所起的作用就开始起到重要作用，当薪结强度达到极限破坏强度时，颗粒之间发生相对位移，滑动和转动随着勃结介质作用消失而发生，这也是粒子流所研究的对象不局限于颗粒状介质的一个基本原因。颗粒流最基本的特征有:(1)允许粒子发生有限位移和转动粒子间可以完全脱离;(2)在计算过程中能够自动辩识新的接触。PFC的基本思想是采用介质最基本单元一粒子和最基本的力学关系一牛顿第二定律来描述介质的复杂力学行为，故是一种本质性和根本性的描述。该数值计算理论在应用环节的思路和方法，因为其基本思想的不同，很大程度上不同于其他连续和非连续力学理论方法。这些差别主要体现在如下几个方面:(l)模型介质的宏观基本物理力学特征不可能通过直接赋值的形式实现，只有粒子的几何特性和粒子间接触的细(微)观力学参数可以赋值，粒子的几何参数包括介质颗粒大小和分布(土体的颗粒级配和岩石的结构)。接触特性包括接触方式和接触力学特征(刚度和强度)。介质的总体力学特征取决于粒子的这些基本特性，改变这些基本特性就意味着显著改变了介质的宏观力学特性。(2)介质的初始条件如地应力场条件会影响介质的结构特征(粒子集合体的密度)从而影响其物理力学特性，因此，地应力场条件必须作为模型特性的一个与介质基本物理力学特性相关联且不是独立的因素考虑，这与以往的数值计算方法完全不同也非常自然地描述了应力环境对介质(特别是岩体)基本物理力学特性的影响。(3)由于介质的力学特性取决于介质内部粒子的结构和接触特征，因此，计算中不需要给介质赋以某种本构关系模型。介质的本构特征将全部由介质内部粒子之间状态特征的变化体现出来，粒子间接触的破坏和发展标志着介质整体力学特性由线形向非线性转化，由弹性向弹塑性转化。(4)构建PFC模型和进行相应的运算准备工作必须使用PFC的二次开发功能，可通过自编程操纵PFC实现上述目的。颗粒流模型的基本假定颗粒流方法在模拟过程中作了如下假设:(l)颗粒单元为刚性体，本身不会破坏;(2)接触发生在很小的范围内、即点接触;(3)接触特性为柔性接触，接触处允许有一定的“重叠”量;(4)“重叠”量的大小与接触力大小有关，与颗粒大小相比，“重叠”量很小;(5)接触处可以有豁结强度; (6)所有的颗粒是圆形(PFCZD)或球体(PFC3D)，也可以用到簇逻辑机理生成任意形状的超级颗粒。每一个簇单元由一系列颗粒重叠而成为边界可以变形的刚体。
