HUSCAP logo Hokkaido Univ. logo

Hokkaido University Collection of Scholarly and Academic Papers >
Graduate School of Engineering / Faculty of Engineering >
Peer-reviewed Journal Articles, etc >

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

Files in This Item:
Article-NAG-revision-NormalCopy.pdf4.15 MBPDFView/Open
Please use this identifier to cite or link to this item:

Title: Development of a GPGPU-parallelized hybrid finite-discrete element method for modeling rock fracture
Authors: Fukuda, Daisuke Browse this author →KAKEN DB
Mohammadnejad, Mojtaba Browse this author
Liu, Hongyuan Browse this author
Dehkhoda, Sevda Browse this author
Chan, Andrew Browse this author
Cho, Sang-Ho Browse this author
Min, Gyeong-Jo Browse this author
Han, Haoyu Browse this author
Kodama, Jun-ichi Browse this author →KAKEN DB
Fujii, Yoshiaki Browse this author →KAKEN DB
Keywords: CUDA C
C plus plus
fracture process analysis
hybrid FDEM
impact loading
quasi-static loading
Issue Date: Jul-2019
Publisher: John Wiley & Sons
Journal Title: International journal for numerical and analytical methods in geomechanics
Volume: 43
Issue: 10
Start Page: 1797
End Page: 1824
Publisher DOI: 10.1002/nag.2934
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.
Rights: This is the peer reviewed version of the following article: Fukuda, D, Mohammadnejad, M, Liu, HY, et al. Development of a GPGPU‐parallelized hybrid finite‐discrete element method for modeling rock fracture. Int J Numer Anal Methods Geomech. 2019; 43: 1797– 1824. , which has been published in final form at This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Use of Self-Archived Versions.
Type: article (author version)
Appears in Collections:工学院・工学研究院 (Graduate School of Engineering / Faculty of Engineering) > 雑誌発表論文等 (Peer-reviewed Journal Articles, etc)

Submitter: 福田 大祐

Export metadata:

OAI-PMH ( junii2 , jpcoar )

MathJax is now OFF:


 - Hokkaido University