EDEM软件介绍

基于离散元方法的EDEM软件介绍2012年09月1离散元方法简介传统的力学研究都是建立在连续性介质假设的基础上的，即认为研究对象是由相互连接没有间隙的大量微团构成。然而，这种假设在有些领域并不适用，如：岩土力学。1971年，CUNDALL提出的一种处理非连续介质问题的数值模拟方法，离散元方法（DiscreteElementMethod，简称DEM），理论基础是结合不同本构关系（应力-应变关系）的牛顿第二定律。随后，这种方法被越来越广泛的应用于涉及颗粒系统地各个领域。通过求解系统中每个颗粒的运动学和动力学方程（碰撞力及场力），不断地更新位置和速度信息，从而描述颗粒系统行为。EDEM软件介绍EDEM主要由三部分组成：Creator、Simulator和Analyst。Creator是前处理工具，完成几何结构导入和颗粒模型建立等工作；Simulator是求解器，用于模拟颗粒体系的运动过程；Analyst是后处理工具，对计算结果进行各种处理。图1.1EDEM结构框架及功能Creator——EDEM的前处理工具EDEM的前处理工具Creator主要完成建模工作，包括：材料参数设置，确定颗粒形状、颗粒产生方法、几何设备导入及运动特性描述等。2Creator的颗粒几何形状建模现实世界中，颗粒状物质形状各异、千差万别，而形状对颗粒体系的运动情况又有着重要的影响。EDEM的前处理工具可以精确描述颗粒的几何外形，Creator通过球面填充技术，将颗粒的表面用若干球面的组合表征，不仅能体现颗粒的非球形特征，又可以使颗粒的接触满足球面接触的物理模型。图1.2颗粒建模界面图1.3采用球面填充方法表征颗粒形状3图1.4各种形状的颗粒颗粒工厂技术EDEM特有的颗粒工厂技术（ParticleFactoryTM），可以根据用户需要，设置颗粒的初始位置、生成速率、颗粒种类、粒径分布等。图1.5按正态分布生成的颗粒4图1.6指定颗粒生成的位置（红色区域）EDEM的材料数据库EDEM的材料数据库允许客户将所关注领域内的各种材料整理成库，在每次建模仿真时，直接从库里导出，不仅减少了用户建模时查找数据的繁琐工作，实现了相关数据的管理和积累。此外，EDEM也内置了一些材料的特性参数供用户参考。图1.7Creator中的材料库丰富的接触模型接触模型是离散单元法的重要基础，其实质就是拟静态下颗粒固体的接触力学弹塑性分析结果。接触模型的分析计算直接决定了粒子所受的力和力矩的大小。需要指出，尽管接触关系是非线性的，仍近似采用叠加原理。离散元法的接触模型有多种，接触力的计算方法也各不相同，但是整体计算的原理都是相同的。EDEM中内置的接触模型列举如下：Hertz-Mindlin(noslip)——基本接触模型，计算颗粒接触时的基本作用力；5Hertz-MindlinwithBonding——在Hertz-Mindlin(noslip)的基础上，考虑了颗粒内聚力的影响；Hertz-MindlinwithHeatConduction——在Hertz-Mindlin(noslip)的基础上，计算颗粒接触后的热量传递情况；LinearCohesion——描述颗粒结块，粘结的物理模型；LinearSpring——基本接触模型，计算颗粒接触是的基本作用力；MovingPlane——颗粒与运动平面间接触过程的物理模型；Tribocharging——带电颗粒相互作用的物理模型；激励的施加EDEM可以施加各种不同的复杂激励，如：1）连续运动（平动、转动及复合运动）2）间歇性运动激励（振动）EDEM的几何文件接口EDEM支持多种格式的几何文件。通过几何文件，用户可以将其他CAD软件建立的机械设备几何模型导入EDEM，减少了重复建模和对模型的修整，提高建模效率。目前，EDEM所支持的几何文件格式如下表所示。图1.8导入几何模型文件扩展名类型应用软件源.igs.iges几何文件IGES6.msh几何文件FluentMesh.stl几何文件STL.stp.step几何文件STEP表1.1EDEM支持的标准文件类型EDEM先进的求解技术利用DEM求解器Simulator进行模拟，可以快速、有效地监测离散颗粒间的碰撞；能够选用动态时间步长；软件既可以在单个处理器上运行，也支持多处理器并行计算；可以通过模型参数的可视化图表来分析模拟结果，从而快速地识别趋向和修正结果。图1.9Simulator界面EDEM利用离散单元法进行计算，其基本思想是把介质看作由一系列离散的独立运动的单元（粒子）所组成，单元的尺寸是细观的，利用牛顿第二定律建立每个单元的运动方程，并用显示中心差分法求解，整个介质的变形和演化由各单元的运动和相互位置来描述。在解决连续介质力学问题时，除了边界条件以外，还有3个方程必须满足，即本构方程、平衡方程和变形协调方程。进行离散元数值计算时，往往通过循环计算的方式，跟踪计算材料颗粒的移动状况。其内部计算关系如图1.10所示。7图1.10计算流程示意图每一次循环包括两个主要的计算步骤：(1)由作用力、反作用力原理和相邻颗粒间的接触本构关系确定颗粒间的接触作用力和相对位移；(2)由牛顿第二定律确定由离散单元法及其在EDEM上的实践相对位移在相邻颗粒间产生的新的不平衡力，直至要求的循环次数或颗粒移动趋于稳定或颗粒受力趋于平衡。并且计算过程按照时步迭代遍历整个颗粒体，计算时间的长度可以根据需要自行设定。EDEM并行计算能力对于某些特定的问题，DEM模拟可能会所需要很大的计算量，单处理器计算往往难于满足仿真的要求。EDEM的并行计算模块是分析该类问题的有力工具。图1.11给出Processor数目对计算速度的影响，能够知道随着处理器数目的增加，计算速度近似线性增长。该测试是在2xQuad-core3.2GHzIntelXeon、8GB内存和WindowsXPProfessionalx64Edition(SP2)操作系统上进行的。8图1.11Processor个数变化对计算速度影响Analyst——EDEM强大的后处理工具后处理模块（Analyst）提供了对仿真结果进行分析和判断的非常丰富的工具。可以用动画方式显示仿真过程，用图表显示仿真结果，创建动画文件并可输出结果数据到其他处理程序中。图1.123D动画显示图1.13剖面显示9图1.14各种类型的统计图表（a）原型显示粒子（b）向量显示粒子（c）流方式显示粒子（d）球锥方式显示粒子图1.15不同类型的粒子显示方式EDEM软件的LICENSE管理及运行平台EDEM具有先进license管理和软件配置机制，可以最大限度利用现有的license资源，避免不必要的license资源浪费。FLUENT软件的license管理机制和运行平台具有以下优势：1)EDEM的license可以在支持TCP/IP协议的局域网环境浮动，license服务可以安装在WINDOWS、LINUX平台。2)所支持的操作系统：WindowsXPServicePack3–32bit(微软XP32位操作系统)10WindowsXPServicePack2–64bit(微软XP64位操作系统)RedhatEnterpriseLinuxWorkstation5–64bit(红帽企业版Linux第5版)11流体－颗粒系统数值模拟的FLUENT－EDEM解决方案北京海基嘉盛科技有限公司2012年10月08日12一、概述绝大多数固态物质的个体是以颗粒状的外形存在的，即：有特定的尺寸和形状，与外界有有限的边界。自然界中的矿石，种子，沙粒，工业产品中的药片、糖果等都是典型的颗粒。通常，无论是在自然界，还是人类生产实践中，都会涉及到了流体与颗粒相互作用（包括：质量交换、动量交换和能量交换等）。如：沙尘暴，水土流失，农作物的干燥，工业上使用的各种流化床，旋流分离器以及气力输运设备等。研究这种相互作用，对人们的生产生活有着重要意义：不仅为提高生产力，更能为改善人类的生存环境提供指导依据。我们将涉及流体流动换热和颗粒运动的体系称为“流体-颗粒系统”。该类系统的研究难点在于：1）流体本身就具有形态不固定，变化无常，难于观察和测量的特点；2）大量颗粒进行相互碰撞：不同时刻和位置，每个颗粒的运动、受力情况都有所不同；3）流体与颗粒相互影响，形成强烈的耦合作用，更加大了系统的复杂度。在以往的研究中，实验研究占很大的比重，主要通过测定或统计的方法来获取系统的宏观指标。另一些则是通过模型简化，进行机理性的研究。随着计算机技术和数值算法的发展，越来越多的科学家和研究人员投入到数值仿真的研究中来，FEA（有限元分析）方法和CFD（计算流体力学）技术成为应用力学中发展最为迅速、活跃的分支。针对流体-颗粒系统的数值模拟研究，主要采用基于CFD方法的多相流技术和CFD-DEM耦合方法。二、气固（液固）两相流技术发展状况在研究初期，由于没有很好的描述颗粒系统的计算模型，人们更倾向于以研究流体为切入点（研究该类系统的科学家和研究人员通常是流体力学专业出身），将系统中大量的颗粒假设为一种准流体——颗粒流，从而产生了气固（液固）两相流技术。气固两相物质所组成的流动系统称为气固两相流系，其中气相通常以连续相形式出现，固相以颗粒或团块的形式处于气相中。气固两相流的流动形态有多种。除了同单相流动那样区分为层流和湍流外，还可以依据两相相对含量（常称为相比）、相界面的分布特性、运动速度、流场几何条件划分流动形态。因此气固两相流研究的首要课题是判断流动形态，并进而分析分散相在连续相中的运动规律及其影响，尤其是对设备的性能影响，包括摩擦阻力、振动和稳定性等。气固两相流的理论分析比单相流困难得多，描述两相流的通用微分方程组至今尚未建立。大量理论工作主要采用以下两类简化模型：均相模型。将两相介质看成是一种混合得非常均匀的混合物，假定处理单相流动的概念和方法仍然适用于两相流，但需对它的物理性质及传递性质作合理的假定；13分相模型。认为单相流的概念和方法可分别用于两相系统的各个相，同时考虑两相之间的相互作用。在气固两相流动研究中，两相间的相互作用是控制流体流动的主要因素，该相互作用主要体现为两相间受力的关系，诸如气固曳力、Magnus升力、Saffman力、Basset力等等。因此，对这些受力的准确描述才能够有效地提高数值预报的准确度。目前，对这些作用力的研究主要采用实验研究与理论分析的形式进行。但实验方法测量颗粒的受力需要有较高分辨率的测量仪器才可能测准颗粒所受到的较小的作用力。随着计算力学的发展，数值计算方法可以在一定范围内替代现有的实验测量，以颗粒受力问题为例，通过计算不同工况条件下的流体绕颗粒流动问题，就可以通过积分颗粒表面的粘性力与压力获得颗粒的受力情况。目前，主流的商业CFD软件都包含气固（液固）两相流模型，用以分析流体—颗粒系统。如：Fluent，CFX等。三、CFD-DEM耦合方法气固两相流模型本质上是将流体—颗粒系统假设为相互掺混的两种流体组成的系统，通过每种流体的浓度变化及分布特性来描述系统，揭示规律。这在研究初期，人们对该类系统知之甚少的情况下，起到了一定的积极作用。但由于该模型本身的局限性，它不能表征颗粒级别上的各种属性，如：颗粒的形状，粒径分布，相互碰撞，运动轨迹等。而且，当颗粒浓度大到某一值时，（固体体积分数为60%，多相流中称为“密相”），基于多相流模型的计算结果也不甚理想。随着研究的深入，人们需要一种新的计算模型，对流体—颗粒系统进行更全面的模拟，CFD-DEM耦合方法应运而生。3.1DEM简介传统的力学研究都是建立在连续性介质假设的基础上的，即认为研究对象是由相互连接没有间隙的大量微团构成。然而，这种假设在有些领域并不适用，如：岩土力学。1971年，CUNDALL提出的一种处理非连续介质问题的数值模拟方法，离散元方法（DiscreteElementMethod,简称DEM），理论基础是结合不同本构关系（应力-应变关系）的牛顿第二定律。随后，这种方法被越来越广泛的应用于涉及颗粒系统地各个领域。通过求解系统中每个颗粒的