Vol.129, No.2, 2021, pp.595-626, doi:10.32604/cmes.2021.017321
A GPU-Based Parallel Algorithm for 2D Large Deformation Contact Problems Using the Finite Particle Method
  • Wei Wang1,2, Yanfeng Zheng1,3, Jingzhe Tang1, Chao Yang1, Yaozhi Luo1,*
1 College of Civil Engineering and Architecture, Zhejiang University, Hangzhou, 310058, China
2 Center for Balance Architecture, Zhejiang University, Hangzhou, 310028, China
3 Architectural Design and Research Institute of Zhejiang University Co., Ltd., Hangzhou, 310028, China
* Corresponding Author: Yaozhi Luo. Email:
Received 01 May 2021; Accepted 26 July 2021; Issue published 08 October 2021
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.
Finite particle method; graphics processing unit (GPU); parallel computing; contact algorithm; large deformation
Cite This Article
Wang, W., Zheng, Y., Tang, J., Yang, C., Luo, Y. (2021). A GPU-Based Parallel Algorithm for 2D Large Deformation Contact Problems Using the Finite Particle Method. CMES-Computer Modeling in Engineering & Sciences, 129(2), 595–626.
This work is licensed under a Creative Commons Attribution 4.0 International License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.