| Title: | Development of a 3D Hybrid Finite-Discrete Element Simulator Based on GPGPU-Parallelized Computation for Modelling Rock Fracturing Under Quasi-Static and Dynamic Loading Conditions |
| Authors: | Fukuda, Daisuke Browse this author →KAKEN DB |
| Mohammadnejad, Mojtaba Browse this author |
| Liu, Hongyuan Browse this author |
| Zhang, Qianbing Browse this author |
| Zhao, Jian Browse this author |
| Dehkhoda, Sevda Browse this author |
| Chan, Andrew Browse this author |
| Kodama, Jun-ichi Browse this author →KAKEN DB |
| Fujii, Yoshiaki Browse this author →KAKEN DB |
| Keywords: | Rocks |
| 3D fracture process analysis |
| FDEM |
| Quasi-static loading |
| Dynamic loading |
| Parallel computation |
| GPGPU |
| CUDA C |
| C plus plus |
| Issue Date: | Mar-2020 |
| Publisher: | Springer |
| Journal Title: | Rock Mechanics and Rock Engineering |
| Volume: | 53 |
| Issue: | 3 |
| Start Page: | 1079 |
| End Page: | 1112 |
| Publisher DOI: | 10.1007/s00603-019-01960-z |
| Abstract: | As a state-of-the-art computational method for simulating rock fracturing and fragmentation, the combined finite-discrete element method (FDEM) has become widely accepted since Munjiza (2004) published his comprehensive book of FDEM. This study developed a general-purpose graphic-processing-unit (GPGPU)-parallelized FDEM using the compute unified device architecture C/C ++ based on the authors' former sequential two-dimensional (2D) and three-dimensional (3D) Y-HFDEM IDE (integrated development environment) code. The theory and algorithm of the GPGPU-parallelized 3D Y-HFDEM IDE code are first introduced by focusing on the implementation of the contact detection algorithm, which is different from that in the sequential code, contact damping and contact friction. 3D modelling of the failure process of limestone under quasi-static loading conditions in uniaxial compressive strength (UCS) tests and Brazilian tensile strength (BTS) tests are then conducted using the GPGPU-parallelized 3D Y-HFDEM IDE code. The 3D FDEM modelling results show that mixed-mode I-II failures are the dominant failure mechanisms along the shear and splitting failure planes in the UCS and BTS models, respectively, with unstructured meshes. Pure mode I splitting failure planes and pure mode II shear failure planes are only possible in the UCS and BTS models, respectively, with structured meshes. Subsequently, 3D modelling of the dynamic fracturing of marble in dynamic Brazilian tests with a split Hopkinson pressure bar (SHPB) apparatus is conducted using the GPGPU-parallelized 3D HFDEM IDE code considering the entire SHPB testing system. The modelled failure process, final fracture pattern and time histories of the dynamic compressive wave, reflective tensile wave and transmitted compressive wave are compared quantitatively and qualitatively with those from experiments, and good agreements are achieved between them. The computing performance analysis shows the GPGPU-parallelized 3D HFDEM IDE code is 284 times faster than its sequential version and can achieve the computational complexity of O(N). The results demonstrate that the GPGPU-parallelized 3D Y-HFDEM IDE code is a valuable and powerful numerical tool for investigating rock fracturing under quasi-static and dynamic loading conditions in rock engineering applications although very fine elements with maximum element size no bigger than the length of the fracture process zone must be used in the area where fracturing process is modelled. |
| Rights: | This is a post-peer-review, pre-copyedit version of an article published in Rock Mechanics and Rock Engineering. The final authenticated version is available online at: http://dx.doi.org/10.1007/s00603-019-01960-z |
| Type: | article (author version) |
| URI: | http://hdl.handle.net/2115/80506 |
| Appears in Collections: | 工学院・工学研究院 (Graduate School of Engineering / Faculty of Engineering) > 雑誌発表論文等 (Peer-reviewed Journal Articles, etc)
|
