|
|||||
|
|
Development of a GPGPU-parallelized hybrid finite-discrete element method for modeling rock fracture |
|
| Authors: | Daisuke Fukuda Mojtaba Mohammadnejad Hongyuan Liu Sevda Dehkhoda Andrew Chan Sang-Ho Cho Gyeong-Jo Min Haoyu Han Jun-ichi Kodama Yoshiaki Fujii |
| Institution: | 1. Faculty of Engineering, Hokkaido University, Hokkaido, 060-8628 Japan
School of Engineering, College of Sciences and Engineering, University of Tasmania, Hobart, Tasmania, 7001 Australia;2. School of Engineering, College of Sciences and Engineering, University of Tasmania, Hobart, Tasmania, 7001 Australia CSIRO Minerals Resources Business Unit, Queensland Centre for Advanced, Technologies, Brisbane, Queensland, 4069 Australia;3. School of Engineering, College of Sciences and Engineering, University of Tasmania, Hobart, Tasmania, 7001 Australia;4. CSIRO Minerals Resources Business Unit, Queensland Centre for Advanced, Technologies, Brisbane, Queensland, 4069 Australia;5. Department of Mineral Resources and Energy Engineering, Chonbuk National, University, Jeollabuk-do, 561-756 South Korea;6. Faculty of Engineering, Hokkaido University, Hokkaido, 060-8628 Japan |
| Abstract: | 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. |
| Keywords: | CUDA C/C++ fracture process analysis GPGPU hybrid FDEM impact loading quasi-static loading rocks |