离散元的基本原理、方法和应用.pdfVIP

7 ( 1一镪 第 13卷第 2期 河 ．海 科 技 进 展 1993年 6月 离散单元法的基本原理、方法及应用 卓家寿 赵 宁 己 ， — — — 一 — — — — — — ～ (河海大学工程力学系 南京 210024) 摘要 本文介绍了离散单元法的基本原理及其基本解法。该法可用于研究边坡失稳、国岩稳定 等不连续介质的大位移问题．采用动卷松驰法求解 ．离散元能够模拟出结构的失稳、凿道爆破的 1 引 言 离散单元法 (discreteeleraentmethod，简称DEM)是近年来解决非连续介质问题的一种颇 具特色、富有发展前景的数值方法。它已在采矿工程、岩土工程 (如边坡及围岩的稳定)等 领域中得到成功的应用，并愈来愈引起工程界和学术界的重视。 DEM 是美国Cundall，P．A于 1971年提出的0 】，最初是从研究具有裂隙节理的岩体开 始的。它把岩体视为被节理切割而成的若干个块体的组合体 l基于岩体的变形主要依赖于软 弱结构面 (如裂隙、节理及层面等)的客观事实，提出了岩块为刚性的假定，以刚性元及其 周界的几何，运动和本构方程为基础 ，采用动态松驰迭代格式 ，建立了求解节理岩块非连续 介质大变形的差分方程。 1978年，Main[，T．和 Cundall，P．A等进一步改进 了原来的剐块离散元模型 ，考虑 了岩 块 自身的变形，提出了可变形块体模型的通用程序UDEC(Universaldicsreteclementcode)】，并 将其推广至模拟岩块破碎和爆炸的运动过程一。同一时期内，Cundall与Straek等还开发了二 维圆形块体的Ball程序r ，用于研究颗粒介质的力学行为，所得结果与Drescher等人用光弹 技术所获得的实验结果相当吻合，这使离散元方法名声大振 ，并为研究颗粒散体介质的本构 关系开辟了一条新途径 。离散单元法可以与边界单元法或有限单元法藕合应用口 ]，解决了远 场岩体为连续介质，近场为不连续介质的问题，大大拓宽 了应用范围，此外，离散单元法还 用于研究散体动力学和边坡的动力稳定性 。 WiBame，J．R，Hocking，G在 l985年发表了 “离散元的理论基础”一文 ；日本的Kawai 在离散元基础上提出了固体力学数值分析的一种新方法——刚体一弹簧元 (rigidelement—spring 收碥 ¨JIl，『l992-Il09 method) l2-，高效且成功地解决了裂隙附近弹塑性应力场等非线性问题的分析 ．由此 ，人们 认识到离散元的深刻的学术价值 。 离散元的理论依据是牛顿第二定律，正如有限元的理论依据是能量原理 ，边界元的理论 依据是 Bctti功互等定理一样，是经典的固体力学原理 ，其基础是牢靠的。从块体单元的几何 形状上 ，离散元分为任意多边形，圆形和椭圆离散元 }从块体单元的变形特征上 ，又有刚性 块离散元和可变形块离散元。 本文主要介绍任意多边形刚性块离散元及其应用 。 2 离散元的基本假定与计算模型 2．1 刚性离散元的基本假定 ‘离 ＼、 圈 I 块体接触 。 图3 块体的集合体及作用于 I块体上的力 由牛顿第二定律 +描述 I块形心的运动方程为 n(1)+ ·m(￡)=∑ (+ ) (1) r 式 中 m为 I块 的惯性