Modelling of dynamic rock fracture process using the finite-discrete element method with a novel and efficient contact activation scheme
Fukuda, D and Liu, H and Zhang, Q and Zhao, J and Kodama, J and Fujii, Y and Chan, AHC, Modelling of dynamic rock fracture process using the finite-discrete element method with a novel and efficient contact activation scheme, International Journal of Rock Mechanics and Mining Sciences, 138 Article 104645. ISSN 1365-1609 (2021) [Refereed Article]
Within the framework of combined finite-discrete element method (FDEM), we propose a novel and efficient semi-adaptive contact activation approach (semi-ACAA). The semi-ACAA adaptively activates contact calculations for continuum solid elements around the cohesive element which has just been subjected to shear softening while its softening function has just satisfied a prescribed threshold. The semi-ACAA is implemented in a three-dimensional (3D) FDEM code parallelized based on general-purpose graphic processing units (GPGPU) to model dynamic fracture processes of marbles in dynamic Brazilian tensile strength (BTS) and uniaxial compressive strength (UCS) tests. The semi-ACAA not only overcomes spurious fracturing mode associated with FDEM simulations using ACAA but also is more physically sound and computationally efficient compared with brute-force contact activation approach (BCAA), which has been prevalent in FDEMs with intrinsic cohesive zone model (ICZM)-based cohesive elements. Furthermore, it is proven that the spurious fracturing mode is caused by unphysical movements of mesh due to missing contacts during the shear softening of cohesive elements instead of the sudden activation of contact force calculations as indicated in some literature. Finally, a series of simulations of the dynamic BTS and UCS tests conducted by split Hopkinson pressure bar (SHPB) facility reveals that, if the threshold is set to around unity, not only the spurious fracturing mode can be overcome but also the modelling precision of stress waves in intact rocks can be improved. The speed-up times of semi-ACAA compared with BCAA are between 1.23 and 13.5 depending on the intensity of shear cracking in the simulations including mixed-mode cracking. Since the proposed semi-ACAA is simple and can be easily implemented into any FDEM codes with ICZM-based cohesive elements, it is expected to accelerate future developments and applications of FDEM with ICZM-based cohesive elements in rock dynamics such as drilling and blasting.
3D FDEM, GPGPU, SHPB, dynamic fracture, novel activation of contact calculation