基于流形覆盖的岩体宏细观破裂的颗粒离散元法-物理学报.PDF

物理学报 Acta Phys. Sin. Vol. 63, No. 5 (2014) 050202 基于流形覆盖的岩体宏细观破裂的颗粒离散元法∗ 李树忱 平洋 李术才 寇强 马腾飞 冯丙阳 (山东大学岩土与结构工程研究中心, 济南 250061) ( 2013 年10 月24 日收到; 2013 年11 月12 日收到修改稿) 本文基于流形覆盖的思想, 结合颗粒接触模型, 构造一种新的颗粒离散元数值方法即MPD (Manifold Particle Discrete), 该方法可以从宏细观角度分析岩体破裂的全过程. 通过分析流形覆盖和球颗粒模型, 构造 三维流形覆盖的四面体单元, 得到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 . 由于有限元方法的局限性和离 散法的诸多优点, 使得离散元法被不断的继承和发 岩体是由多种矿物颗粒、胶结物及孔隙缺陷等 12 展, 经过近30 年无数专家学者 、研究机构的努 组成的复杂混合体, 一般都随机地赋存大量地质缺 力, 离散元法已经逐步发展成熟. 岩石、混凝土类 陷, 如裂纹、断裂面、节理、孔洞、充填物等, 并且随 脆性材料的破坏过程是一个由连续体向不连续体 机地分布于岩石中. 在外载荷作用下, 岩石内部微 过渡的过程, 采用一般的有限元方法对其进行细观 缺陷的成核、扩展以及相互之间的作用决定了岩石 力学数值模拟时, 存在一些不足和困难之处; 而离 的宏观变形、破裂特性1 . 多年来, 各种理论一直在 散元法等用于散体分析的数值方法在处理结构破 均匀介质假设基础上研究岩体局部破坏, 很少涉及 碎后的动力分析时有很大的优势, 所以很多研究 岩石自身非均匀性引起的应力分布非均匀性而诱 者采用离散元法与其他方法结合的算法来模拟材 发的局部破坏15 . 料破裂的动力学行为68 . Owen, Munjiza 和唐春 对于非均匀材料问题, 国内外学术界进行了长 安等做了大量的研究914 . Klerck 和Cottrell 构建 期探索, 特别是近20 年来取得了很大的进展, 并逐 了准脆性材料在拉压应力场中破裂过程的力学模 渐形成了一个新型学科分支2 . 从研究方法来看, 型, 并发现随着材料的破裂同时产生了肉眼可见 比较流行的主要有两种类型: 一种是建立在各种各 的裂纹10 . 唐春安等1314 运用岩石破裂过程分析 样自洽模型基础上的有限元与边界元方法, 另一种 REPA2D 系统, 通过对岩石试样中预置的倾斜裂纹 是建立在统计技术上的格形化有限元方法24 . 这 扩展过程的数值模拟, 研究了材料非均匀性对岩石 两种方法已被成功地应用于各种非均匀材料问题, 介质中裂纹扩展模式的影响. 尽管有关非均匀性材 特别在复合材料力学性能研究和脆性材料破坏过