Open Access Repository

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

Fukuda, D, Mohammadnejad, M ORCID: 0000-0001-8570-5074, Liu, H ORCID: 0000-0002-5437-4695, Dehkhoda, S, Chan, A ORCID: 0000-0003-0042-8448, Cho, S-H, Min, G-J, Han, H, Kodama, J-i and Fujii, Y 2019 , 'Development of a GPGPU‐parallelized hybrid finite-discrete element method for modeling rock fracture' , International Journal for Numerical and Analytical Methods in Geomechanics, vol. 43, no. 10 , pp. 1797-1824 , doi: 10.1002/nag.2934.

Full text not available from this repository.


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.

Item Type: Article
Authors/Creators:Fukuda, D and Mohammadnejad, M and Liu, H and Dehkhoda, S and Chan, A and Cho, S-H and Min, G-J and Han, H and Kodama, J-i and Fujii, Y
Keywords: CUDA C/C++, fracture process analysis, GPGPU, hybrid FDEM, impact loading, quasi-static loading, rocks
Journal or Publication Title: International Journal for Numerical and Analytical Methods in Geomechanics
Publisher: John Wiley & Sons Ltd
ISSN: 1096-9853
DOI / ID Number: 10.1002/nag.2934
Copyright Information:

© 2019 John Wiley & Sons, Ltd.

Related URLs:
Item Statistics: View statistics for this item

Actions (login required)

Item Control Page Item Control Page