其它数值方法.pptVIP

四、UDEC（Universal Discrete Element Code）应用简介 Itasca 公司的主要产品 two dimensional continuum, with joints three dimensional continuum, with joints three dimensional DEM* spheres clumps two dimensional DEM* disks clumps two dimensional DEM* polygonal bodies three dimensional DEM* polyhedral bodies UDEC 非连续介质被作为是凸块或凹块的集合体，块体可以是刚体或变形体。 不连续面被看作是块体之间的边界条件。 沿着不连续面的运动是由切向和法向的线性及非线性力-位移关系控制的。 地质体或类似材料可以由很多内置的块体和节理本构模型来描述；还有可供开发的自定义模型。 是一种非常适合于用来模拟具有非线性力学行为的非连续材料（包含多个相互切割的非连续面）。尤其是，它具有如下特征： 5. 有平面-应变、平面-应力及轴对称三种问题模型。 6. 具有用来描述岩石-结构相互作用的结构单元模型，如锚杆、桩、梁、衬砌，喷射混凝土等。 7. 能够进行静力和动力问题分析 8. 沿着节理的暂态和稳态流. 9. 粘弹性和粘塑性（蠕变）模型. 10. 热分析能力，流固耦合分析能力 Overview of UDEC operation (1) - 工程模拟通常由很长的一系列操作组成 - UDEC数据文件可以很容易地利用文本编辑器进 行修改，而且多个文件还可以连接在一起 这种输入文件的方式提供了一个很好的可以用来 保存分析稳定记录的方法. 命令流驱动结构允许通过操作程序的输入与输出 来开发其前后处理程序 UDEC 是命令流驱动程序 * 资源与地球科学学院 * 第六章 其它数值方法 第一节 离散单元法（详细） 第二节 边界单元法（自学） 第三节 数值方法的耦合（自学） 主要内容 第一节 离散单元法 一、概述 二、基本原理 四、UDEC（Universal Discrete Element Code）应用简介 三、运动方程的解法 推荐参考书：王泳嘉.离散单元法及其在岩土力学中的应用. 东北工学院出版社，1991 一、概述 有限元和边界元法作为一种数值方法，在研究工程地质问题中是以连续介质力学为基础的。对于不连续面的模拟，有限元法一般采用节理单元，边界元法使用位移间断法，有时采用节理化岩体模型。但是这些方法在模拟岩质边坡的失稳运动、矿山开采的岩层移动等这类块体可能沿弱面作大位移或大角度旋转的运动则显得无能为力。 离散元法是专门用来解决不连续介质问题的数值模拟方法。该方法把节理岩体视为由离散的岩块和岩块间的节理面所组成,允许岩块平移、转动和变形,而节理面可被压缩、分离或滑动。因此,岩体被看作一种不连续的离散介质。其内部可存在大位移、旋转和滑动乃至块体的分离,从而可以较真实地模拟节理岩体中的非线性大变形特征。 离散元法的一般求解过程为:将求解空间离散为离散元单元阵,并根据实际问题用合理的连接元件将相邻两单元连接起来;单元间相对位移是基本变量,由力与相对位移的关系可得到两单元间法向和切向的作用力;对单元在各个方向上与其它单元间的作用力以及其它物理场对单元作用所引起的外力求合力和合力矩,根据牛顿运动第二定律可以求得单元的加速度;对其进行时间积分,进而得到单元的速度和位移。 产生与发展：该方法是20世纪70年代发展起来的，目前已在数值模拟理论与工程应用方面取得了很大进展。该方法的基本特征在于允许各离散块体发生平动和转动，甚至发生分离，弥补了有限元法或边界元法的介质连续和小变形的限制，因而特别适合块状及裂隙介质的大变形及破坏问题的分析。 该方法首先由P.A.Cundall首先提出，于20世纪80年代中期由王泳嘉等人引入到我国。目前，离散元法发展迅速，美国有二维和三维的离散元程序UDEC和3DEC。我国有2D-BLOCK和3D-BLOCK。 应用领域：边坡、巷道与采场、地下开采、地震、爆炸、核废料储存、散体介质运动、断裂、地下水渗流、热传导等。 离散单元法所用的求解方法有静态松弛法和动态松弛法。 静态松弛法直接寻找块体失去平衡后达到再平衡时力和位移的关系，通过求解联立平衡方程组，并不断迭代使块体的残余余力和力矩趋于零，从而获得结果。这一方法不考虑粘性阻尼和计算时步，避免了动态松弛法确定这些参数时所遇到的困难。但该方法不足之处是在求