颗粒离散元法及其在球床高温堆球流研究中应用的可能性分析.pdfVIP

http :// 颗粒离散元法及其在球床高温堆球流研究中应用的 可能性分析 胡文平，杨星团，姜胜耀 清华大学核能与新能源技术研究院，北京 ( 100084 ) Emai：hwp04@ 摘 要：简述了离散元法（DEM ）的基本原理及其在国内外的发展现状，着重介绍了离散 元法在工程上诸多应用中的有关颗粒元介质的研究情况。分析了离散元法应用于球床式高温 气冷堆内球流运动分析的可行性，指出离散元法应用于研究球床式高温气冷堆内球流运动规 律是可行的，但若用于工程实践仍需解决一些重要的技术问题。 关键词：离散元法；颗粒元；高温气冷堆；球流运动 1. 引言 离散体在自然界和工程中普遍存在，如泥石流、制药工程中的粉体、料仓卸料问题中的 颗粒体及岩土工程中的块体等。它们以力学特性上的颗粒物性、形状的分散性及运动特性上 的瞬态、波动、聚团和破裂等特性，而区别于连续介质。工程上研究连续介质的方法很多且 较成熟，离散体的上述特性使得难以采用连续介质力学方法分析工程上的离散体问题。 随着计算力学及计算机技术的发展，出现了许多用于离散体力学计算的数值方法，较早 [1]主要用于分析微观的原子、分子世界的非连续 的分子动力学模型（Molecular Dynamics ） 介质数值方法；离散元法（Discrete/Distinct Element Method, DEM ）[2] 、不连续变形分析 [3]方法、光滑粒子法 （Smoothed Particle Hydrodynamics, （Deformed Discrete Analysis, DDA ） [4-5] [6] [7] SPH） 、改进的光滑粒子法（CSPH ） 以及扩展单元法（Diffuse Element Method ） 等 许多数值方法都可用于宏观的非连续体介质的数值研究。本文介绍其中的离散元法，试图探 讨其用于球床式高温气冷堆内球流运动研究的可能性。 离散元法是美国学者Cundall于上世纪七十年代首先提出的，称之块体离散元法，最初 的研究对象是岩石块体的渐进运动过程。随着离散元法30 多年的发展，其在岩石工程、粉 体工程、矿山工业和农业机械等许多涉及离散介质的领域内都发挥了其它数值算法无法取代 的作用。但是，从总体上看对离散元法的研究主要集中在利用其计算、模拟工程问题，对方 法本身的理论和算法的研究相对较少。继提出块体离散元法之后，Cundall和Strack[8-10]又提 出了更适合于颗粒体介质的颗粒离散元法，并开发了二维圆盘（BALL ）和三维圆球 （TRUBALL ）程序；英国Aston大学的Thornton研究组对离散元法的发展做出了重要贡献。 他抛弃Cundall 的弹簧— 阻尼器模型，发展出球体接触力学模型，并对TRUBALL程序进行全 面改造，形成了GRANULE程序。 在离散元商用软件的开发方面，最有名的要属 Cundall 加盟的国际工程咨询公司 (ITASCA) ，该公司开发了一系列离散元商业软件。其开发的FLAC(Fast Lagrangian Analysis of Continua)和 FLAC3D 离散元软件主要用于模拟二维和三维岩土或其他材料的非线性力学行 为，解决大变形大应变、非线性及非稳定系统如完全塌方等问题；二维PFC2D 和PFC3D 软 件分别为基于二维圆盘和三维圆球单元的离散元软件，主要用于研究颗粒集合体的破裂和破 裂发展问题以及颗粒的流动 （大位移）和材料的混合等问题。Thornton 的GRANULE 程序， 可进行不同形状的干、湿颗粒的流动研究，离散本构关系的细观力学分析，结块碰撞及料仓 -1-