Investigation of the influence of model control parameters on fracture characteristics of GPU parallel FDEM

Surrounding rock mass fracture characteristics play a significant role in the understanding of the CO 2 geological storage and utilization (CGSU) engineering practices in abandoned mines. The combined finite discrete element method (FDEM) shows advantages in simulating fracture and fragmentation of rock-like materials, however, many computational parameters and a lack of basis for accurate values affect the simulation results. For systematic explorations of the influence of the effect of model parameters with different time steps, this study conducted different loading rate specimen tests and unloading rate tests both in laboratory-scale tests and field-scale based on the CUDA-based GPU parallel FDEM program. In laboratory-scale uniaxial compressive strength (UCS) tests, a smaller loading rate ensures a quasi-static loading process and should be below 0.1 ms −1 for simulation. Then, continuous optimization and improvement of the GPU parallel FDEM in tunnel excavation were proposed, and the model parameters reflect the continuous improvement of the simulation results. In the process, the tunnel excavation simulation method with the reduction rate of the opening zone’s Young's modulus in the excavation was performed to investigate the unloading mode and rate by continuing to optimize the GPU parallel FDEM program and model parameters. Besides, the main factor in fracture mode and failure mechanism of the surrounding rock mass also was calibrated. The results indicate that the system kinetic energy of the model is maintained at a small level with the reasonable unloading mode and critical threshold set at a small value, the damping parameters, dissipation mechanism, and excavation fractures are clearer and reasonable, and the computational cost is significantly reduced with GPU parallel FDEM.