基于流形覆盖的岩体宏细观破裂的颗粒离散元法.pdfVIP

物 理 学报 ActaPhys．Sin．Vo1．63，No．5(2014) 050202 基于流形覆盖的岩体宏细观破裂的颗粒离散元法冰 李树忱 平洋十 李术才 寇强 马腾飞 冯丙阳 (山东大学岩土与结构工程研究中心，济南 250061) (2013年 10月24日收到；2013年 l1月12日收到修改稿) 本文基于流形覆盖的思想，结合颗粒接触模型，构造一种新的颗粒离散元数值方法即MPD (Mani~ld ParticleDiscrete)，该方法可以从宏细观角度分析岩体破裂的全过程．通过分析流形覆盖和球颗粒模型，构造 三维流形覆盖的四面体单元，得到MPD模型的平衡方程和相容方程．通过数值算例，验证了所提出的数值方 法的正确性和分析岩石破裂的可行性． 关键词：岩体破裂，MPD模型，应力一应变关系，流形覆盖 PACS：02．40．Sf，45．10．一b，61．43．Bn，75．40．Mg DOI：10．7498／aps．63．050202 与运动规律的数值计算方法，是上世纪末 由Peter 1 引 言 Cundall提出的[6】-由于有限元方法的局限性和离 散法的诸多优点，使得离散元法被不断的继承和发 岩体是由多种矿物颗粒、胶结物及孔隙缺陷等 展，经过近30年无数专家学者 [112】、研究机构的努 组成的复杂混合体，一般都随机地赋存大量地质缺 力，离散元法已经逐步发展成熟．岩石、混凝土类 陷，如裂纹、断裂面、节理、孔洞、充填物等，并且随 脆性材料的破坏过程是一个 由连续体向不连续体 机地分布于岩石 中．在外载荷作用下，岩石 内部微 过渡的过程，采用一般的有限元方法对其进行细观 缺陷的成核、扩展 以及相互之间的作用决定了岩石 力学数值模拟时，存在一些不足和困难之处；而离 的宏观变形、破裂特性 1【】．多年来，各种理论一直在 散元法等用于散体分析的数值方法在处理结构破 均匀介质假设基础上研究岩体局部破坏，很少涉及 碎后的动力分析时有很大的优势，所 以很多研究 岩石 自身非均匀性引起的应力分布非均匀性而诱 者采用离散元法与其他方法结合的算法来模拟材 发的局部破坏 [1-5】． 料破裂的动力学行为 6[-s】．Owen，Munjiza和唐春 对于非均匀材料问题，国内外学术界进行了长 安等做了大量的研究[9--14]．Klerck和Cottrell构建 期探索，特别是近20年来取得 了很大的进展，并逐 了准脆性材料在拉压应力场 中破裂过程 的力学模 渐形成了一个新型学科分支 2【]．从研究方法来看， 型，并发现随着材料的破裂 同时产生了肉眼可见 比较流行的主要有两种类型：一种是建立在各种各 的裂纹 [101．唐春安等 [1a,141运用岩石破裂过程分析 样 自洽模型基础上的有限元与边界元方法，另一种 REPA2D系统，通过对岩石试样中预置的倾斜裂纹 是建立在统计技术上的格形化有限元方法 [2-4】．这 扩展过程的数值模拟，研究了材料非均匀性对岩石 两种方法已被成功地应用于各种非均匀材料问题， 介质中裂纹扩展模式的影响．尽管有关非均匀性材 特别在复合材料力学性能研究和脆性材料破坏过 料的三维破裂过程分析的研究受到人们的关注，但 程 []的演化方面取得了令人瞩 目的成果． 对于岩体类脆性材料破裂全过程的研究进展却极 离散元法是一种用于分析散粒物料力学行为 其缓慢． 国家重点基础