Development of a GPGPU‐parallelized hybrid finite‐discrete element method for modeling rock fracture

The hybrid finite‐discrete element method (FDEM) is widely used for engineering applications, which, however, is computationally expensive and needs further development, especially when rock fracture process is modeled. This study aims to further develop a sequential hybrid FDEM code formerly proposed by the authors and parallelize it using compute unified device architecture (CUDA) C/C++ on the basis of a general‐purpose graphics processing unit (GPGPU) for rock engineering applications. Because the contact detection algorithm in the sequential code is not suitable for GPGPU parallelization, a different contact detection algorithm is implemented in the GPGPU‐parallelized hybrid FDEM. Moreover, a number of new features are implemented in the hybrid FDEM code, including the local damping technique for efficient geostatic stress analysis, contact damping, contact friction, and the absorbing boundary. Then, a number of simulations with both quasi‐static and dynamic loading conditions are conducted using the GPGPU‐parallelized hybrid FDEM, and the obtained results are compared both quantitatively and qualitatively with those from either theoretical analysis or the literature to calibrate the implementations. Finally, the speed‐up performance of the hybrid FDEM is discussed in terms of its performance on various GPGPU accelerators and a comparison with the sequential code, which reveals that the GPGPU‐parallelized hybrid FDEM can run more than 128 times faster than the sequential code if it is run on appropriate GPGPU accelerators, such as the Quadro GP100. It is concluded that the GPGPU‐parallelized hybrid FDEM developed in this study is a valuable and powerful numerical tool for rock engineering applications.

[1]  Qianbing Zhang,et al.  A numerical study of rock scratch tests using the particle-based numerical manifold method , 2018 .

[2]  Chengzeng Yan,et al.  A 2D fully coupled hydro-mechanical finite-discrete element model with real pore seepage for simulating the deformation and fracture of porous medium driven by fluid , 2018 .

[3]  Jiansheng Xiang,et al.  A numerical investigation of mesh sensitivity for a new three-dimensional fracture model within the combined finite-discrete element method , 2016 .

[4]  Hong Zheng,et al.  FDEM-flow3D: A 3D hydro-mechanical coupled model considering the pore seepage of rock matrix for simulating three-dimensional hydraulic fracturing , 2017 .

[5]  Dwayne D. Tannant,et al.  Discrete element method contact model for dynamic simulation of inelastic rock impact , 2007, Comput. Geosci..

[6]  John-Paul Latham,et al.  On the validation of DEM and FEM/DEM models in 2D and 3D , 2009 .

[7]  Jeen-Shang Lin,et al.  On the critical failure mode transition depth for rock cutting , 2013 .

[8]  Giovanni Grasselli,et al.  Continuum–discontinuum analysis of failure mechanisms around unsupported circular excavations in anisotropic clay shales , 2014 .

[9]  Davide Elmo,et al.  An Integrated Numerical Modelling–Discrete Fracture Network Approach Applied to the Characterisation of Rock Mass Strength of Naturally Fractured Pillars , 2010 .

[10]  Ning Guo,et al.  A coupled FEM/DEM approach for hierarchical multiscale modelling of granular media , 2014 .

[11]  Giovanni Grasselli,et al.  A review of discrete modeling techniques for fracturing processes in discontinuous rock masses , 2014 .

[12]  Hongyuan Liu,et al.  Hybrid finite-discrete element modelling of dynamic fracture and resultant fragment casting and muck-piling by rock blast , 2017 .

[13]  Ning Guo,et al.  Multiscale Insights Into Borehole Instabilities in High‐Porosity Sandstones , 2018 .

[14]  O. K. Mahabadi,et al.  Three-dimensional FDEM numerical simulation of failure processes observed in Opalinus Clay laboratory samples , 2014 .

[15]  Glaucio H. Paulino,et al.  Extrinsic cohesive modelling of dynamic fracture and microbranching instability in brittle materials , 2007 .

[16]  Ning Guo,et al.  Parallel hierarchical multiscale modelling of hydro-mechanical problems for saturated granular soils , 2016 .

[17]  Jiansheng Xiang,et al.  A numerical study of fracture spacing and through-going fracture formation in layered rocks , 2017 .

[18]  O. K. Mahabadi,et al.  Y-Geo: New Combined Finite-Discrete Element Numerical Code for Geomechanical Applications , 2012 .

[19]  O. K. Mahabadi,et al.  Numerical simulation of acoustic emission in brittle rocks by two-dimensional finite-discrete element analysis , 2013 .

[20]  B. S. A. Tatone,et al.  A calibration procedure for two-dimensional laboratory-scale hybrid finite–discrete element simulations , 2015 .

[21]  Chaoshui Xu,et al.  Discrete element modelling of rock cutting: from ductile to brittle transition , 2015 .

[22]  Ning Guo,et al.  Multiscale modeling and analysis of compaction bands in high-porosity sandstones , 2018 .

[23]  Jiansheng Xiang,et al.  Implementation of an Empirical Joint Constitutive Model into Finite-Discrete Element Analysis of the Geomechanical Behaviour of Fractured Rocks , 2016, Rock Mechanics and Rock Engineering.

[24]  Weiji Liu,et al.  The analysis of ductile-brittle failure mode transition in rock cutting , 2018 .

[25]  Sang-ho Cho,et al.  Influence of the applied pressure waveform on the dynamic fracture processes in rock , 2004 .

[26]  Jiansheng Xiang,et al.  The Virtual Geoscience Workbench, VGW: Open Source tools for discontinuous systems ☆ , 2010 .

[27]  Davide Elmo,et al.  Volumetric Fracture Intensity Measurement for Improved Rock Mass Characterisation and Fragmentation Assessment in Block Caving Operations , 2015, Rock Mechanics and Rock Engineering.

[28]  W. Benz,et al.  Simulations of brittle solids using smooth particle hydrodynamics , 1995 .

[29]  Hongyuan Liu,et al.  Hybrid finite-discrete element modelling of asperity degradation and gouge grinding during direct shearing of rough rock joints , 2016 .

[30]  G. G. Schiava D’Albano,et al.  Space decomposition based parallelization solutions for the combined finite–discrete element method in 2D , 2014 .

[31]  Giovanni Grasselli,et al.  Fracture development around deep underground excavations: Insights from FDEM modelling , 2014 .

[32]  Hong Zheng,et al.  A fully coupled three-dimensional hydro-mechanical finite discrete element approach with real porous seepage for simulating 3D hydraulic fracturing , 2017 .

[33]  Steven F. Quigley,et al.  A fast scalable implementation of the two-dimensional triangular Discrete Element Method on a GPU platform , 2013, Adv. Eng. Softw..

[34]  Michael Garland,et al.  Designing efficient sorting algorithms for manycore GPUs , 2009, 2009 IEEE International Symposium on Parallel & Distributed Processing.

[35]  Hong Zheng,et al.  A coupled thermo-mechanical model based on the combined finite-discrete element method for simulating thermal cracking of rock , 2017 .

[36]  Hongyuan Liu,et al.  Hybrid finite–discrete element modeling of geomaterials fracture and fragment muck-piling , 2015 .

[37]  Mojtaba Mohammadnejad,et al.  An overview on advances in computational fracture mechanics of rock , 2018 .

[38]  Zhijun Wu,et al.  Investigation of the Rock Fragmentation Process by a Single TBM Cutter Using a Voronoi Element-Based Numerical Manifold Method , 2018, Rock Mechanics and Rock Engineering.

[39]  Jiansheng Xiang,et al.  Finite strain, finite rotation quadratic tetrahedral element for the combined finite–discrete element method , 2009 .

[40]  O. K. Mahabadi,et al.  A 2D, fully-coupled, hydro-mechanical, FDEM formulation for modelling fracturing processes in discontinuous, porous rock masses , 2017 .

[41]  A. Munjiza.  The Combined Finite-Discrete Element Method: Munjiza/Discrete Element Method , 2004 .

[42]  Jeen-Shang Lin,et al.  Modeling the ductile–brittle failure mode transition in rock cutting , 2014 .

[43]  Antonio Munjiza,et al.  A framework for grand scale parallelization of the combined finite discrete element method in 2d , 2014, CPM 2014.

[44]  Antonio Munjiza,et al.  Validation of a three-dimensional Finite-Discrete Element Method using experimental results of the Split Hopkinson Pressure Bar test , 2014 .

[45]  J. Latham,et al.  The use of discrete fracture networks for modelling coupled geomechanical and hydrological behaviour of fractured rocks , 2017 .

[46]  O. K. Mahabadi,et al.  Acceleration of a 2D/3D finite-discrete element code for geomechanical simulations using General Purpose GPU computing , 2018, Computers and Geotechnics.

[47]  Davide Elmo,et al.  Damage characterization during laboratory strength testing: A 3D-finite-discrete element approach , 2014 .