A GPU-Based Parallel Algorithm for 2D Large Deformation Contact Problems Using the Finite Particle Method

• Published in: Computer modeling in engineering & sciences Vol. 129; no. 2; pp. 595 - 626
• Main Authors: Wang, Wei, Zheng, Yanfeng, Tang, Jingzhe, Yang, Chao, Luo, Yaozhi
• Format: Journal Article
• Language: English
• Published: Henderson Tech Science Press 2021
• Subjects: Algorithms; Central processing units; Computing time; Contact force; CPUs; Data structures; Deformation effects; Equations of motion; Finite element method; Graphics processing units; Internal forces; Mathematical analysis; Microprocessors; Parallel processing; Quadrilaterals; Search algorithms
• ISSN: 1526-1506; 1526-1492; 1526-1506
• DOI: 10.32604/cmes.2021.017321

Large deformation contact problems generally involve highly nonlinear behaviors, which are very time-consuming and may lead to convergence issues. The finite particle method (FPM) effectively separates pure deformation from total motion in large deformation problems. In addition, the decoupled procedures of the FPM make it suitable for parallel computing, which may provide an approach to solve time-consuming issues. In this study, a graphics processing unit (GPU)-based parallel algorithm is proposed for two-dimensional large deformation contact problems. The fundamentals of the FPM for planar solids are first briefly introduced, including the equations of motion of particles and the internal forces of quadrilateral elements. Subsequently, a linked-list data structure suitable for parallel processing is built, and parallel global and local search algorithms are presented for contact detection. The contact forces are then derived and directly exerted on particles. The proposed method is implemented with main solution procedures executed in parallel on a GPU. Two verification problems comprising large deformation frictional contacts are presented, and the accuracy of the proposed algorithm is validated. Furthermore, the algorithm's performance is investigated via a large-scale contact problem, and the maximum speedups of total computational time and contact calculation reach 28.5 and 77.4, respectively, relative to commercial finite element software Abaqus/Explicit running on a single-core central processing unit (CPU). The contact calculation time percentage of the total calculation time is only 18% with the FPM, much smaller than that (50%) with Abaqus/Explicit, demonstrating the efficiency of the proposed method.