A GPU parallel scheme for accelerating 2D and 3D peridynamics models

Peridynamics (PD) is prevailing in the numerical simulations of damage evolution, but with the cost of far more required computations than traditional methods. This paper proposes a massively parallel implementation scheme for PD simulations with a single-card Graphics Processing Unit (GPU) to reduce the computational cost. The GPU parallel scheme includes two-level parallel modes of particle-mapping and bond-mapping, realizing high parallelism. By reasonably setting the data structure and using the CUDA memory model, realize the efficient utilization of GPU memory resources. Three numerical experiments involving quasi-static and transient problems, 2D and 3D problems, and impact problems are performed. The results show that the proposed parallel scheme can greatly improve the computational efficiency of various PD models while ensuring accuracy. The bond-mapping unrolls the inner loop of the particle-mapping and makes full use of registers and shared memory resources for stronger performance. To further explore the fatigue behavior of rails with hole defects using the GPU parallel scheme, which is difficult to perform using the serial and OpenMP scheme. The results show that the proposed GPU parallel scheme can explore more complex structural fracture problems by greatly reducing the computational cost.