Abstract The Discrete Element Method (DEM) is increasingly used to study the failure behavior of rock. Despite DEM's intrinsic capability to capture the mechanical behavior of discontinua, there remains several open questions that include the numerical modelling of the meso-fracture evolution and behavior of brittle rock bodies. A Cohesive Fracture Model (CFM) designed explicitly for polyhedral shaped DEM particles is proposed for the simulation of the fracture and behavior of brittle rock bodies. A rock body is discretized into a series of rigid polyhedral blocks which are bonded along the boundary faces along the normal and shear directions. A cohesive criterion dictates the normal and shear break strengths of the bonds for the CFM. The computational efficiency of the CFM for polyhedral shaped DEM simulations is enhanced by parallelization over multiple graphical processing units (GPUs) in the Blaze-DEM framework. A series of numerical and laboratory tests are conducted. For marble, these include three-point bending and uniaxial tests to investigate the relationship between the meso-mechanical parameters and macro-mechanical behavior. Following a sensitivity analysis of the meso-mechanical parameters, an inversion procedure is established to estimate the numerical meso-mechanical parameters of the CFM model. Two different failure modes for different rocks are proposed.

中文翻译：

---

基于GPU加速离散元法的岩石破裂过程模拟

摘要 离散元法（DEM）越来越多地用于研究岩石的破坏行为。尽管 DEM 具有捕捉不连续体力学行为的内在能力，但仍有几个悬而未决的问题，包括细观裂缝演化和脆性岩体行为的数值建模。提出了一种专为多面体形 DEM 颗粒设计的内聚断裂模型 (CFM)，用于模拟脆性岩体的断裂和行为。岩石体被离散成一系列刚性多面体块，这些块体沿着边界面沿法向和剪切方向粘合。内聚标准规定了 CFM 键的法向和剪切断裂强度。通过在 Blaze-DEM 框架中的多个图形处理单元 (GPU) 上进行并行化，提高了 CFM 对多面体形状 DEM 模拟的计算效率。进行了一系列数值和实验室测试。对于大理石，这些包括三点弯曲和单轴试验，以研究细观力学参数与宏观力学行为之间的关系。在对细观力学参数进行敏感性分析之后，建立反演程序以估计 CFM 模型的数值细观力学参数。针对不同的岩石提出了两种不同的破坏模式。其中包括三点弯曲和单轴试验，以研究细观力学参数与宏观力学行为之间的关系。在对细观力学参数进行敏感性分析之后，建立反演程序以估计 CFM 模型的数值细观力学参数。针对不同的岩石提出了两种不同的破坏模式。其中包括三点弯曲和单轴试验，以研究细观力学参数与宏观力学行为之间的关系。在对细观力学参数进行敏感性分析之后，建立反演程序以估计 CFM 模型的数值细观力学参数。针对不同的岩石提出了两种不同的破坏模式。