Abstract: A graphics processing unit (GPU)‍-accelerated vector-form particle-element method, i.‍e., the finite particle method (FPM), is proposed for 3D elastoplastic contact of structures involving strong nonlinearities and computationally expensive contact calculations. A hexahedral FPM element with reduced integration and anti-hourglass is developed to model structural elastoplastic behaviors. The 3D space containing contact surfaces is decomposed into cubic cells and the contact search is performed between adjacent cells to improve search efficiency. A connected list data structure is used for storing contact particles to facilitate the parallel contact search procedure. The contact constraints are enforced by explicitly applying normal and tangential contact forces to the contact particles. The proposed method is fully accelerated by GPU-based parallel computing. After verification, the performance of the proposed method is compared with the serial finite element code Abaqus/Explicit by testing two large-scale contact examples. The maximum speedup of the proposed method over Abaqus/Explicit is approximately 80 for the overall computation and 340 for contact calculations. Therefore, the proposed method is shown to be effective and efficient.

结构三维弹塑性接触问题的图形处理器并行有限质点法求解

[1]BartezzaghiA, CremonesiM, ParoliniN, et al., 2015. An explicit dynamics GPU structural solver for thin shell finite elements. Computers & Structures, 154:29-40.

[2]BensonDJ, HallquistJO, 1990. A single surface contact algorithm for the post-buckling analysis of shell structures. Computer Methods in Applied Mechanics and Engineering, 78(2):141-163.

[3]BiabanakiSOR, KhoeiAR, WriggersP, 2014. Polygonal finite element methods for contact-impact problems on non-conformal meshes. Computer Methods in Applied Mechanics and Engineering, 269:198-221.

[4]BorjaRI, 2013. Plasticity: Modeling & Computation. Springer, Berlin, Germany.

[5]CaiY, WangGP, LiGY, et al., 2015. A high performance crashworthiness simulation system based on GPU. Advances in Engineering Software, 86:29-38.

[6]CaoXG, CaiY, CuiXY, 2020. A parallel numerical acoustic simulation on a GPU using an edge-based smoothed finite element method. Advances in Engineering Software, 148:102835.

[7]ChenH, LeiZ, ZangMY, 2014. LC-Grid: a linear global contact search algorithm for finite element analysis. Computational Mechanics, 54(5):1285-1301.

[8]ChengJ, GrossmanM, McKercherT, 2014. Professional CUDA C Programming. John Wiley & Sons, Indianapolis, USA.

[9]DongYK, WangD, RandolphMF, 2015. A GPU parallel computing strategy for the material point method. Computers and Geotechnics, 66:31-38.

[10]FlanaganDP, BelytschkoT, 1981. A uniform strain hexahedron and quadrilateral with orthogonal hourglass control. International Journal for Numerical Methods in Engineering, 17(5):679-706.

[11]GaoB, WuJ, JiaPC, et al., 2022. Experimental and numerical investigation of the polyurea-coated ultra-high-performance concrete (UHPC) column under lateral impact loading. International Journal of Structural Stability and Dynamics, 22(5):2250037.

[12]HallquistJO, GoudreauGL, BensonDJ, 1985. Sliding interfaces with contact-impact in large-scale Lagrangian computations. Computer Methods in Applied Mechanics and Engineering, 51(1-3):107-137.

[13]HeschC, BetschP, 2011. Transient 3D contact problems-NTS method: mixed methods and conserving integration. Computational Mechanics, 48(4):437-449.

[14]HouJM, ShiBS, LiangQH, et al., 2021. A graphics processing unit-based robust numerical model for solute transport driven by torrential flow condition. Journal of Zhejiang University-SCIENCE A (Applied Physics & Engineering), 22(10):835-850.

[15]KangH, KimJ, 2015. Progressive collapse of steel moment frames subjected to vehicle impact. Journal of Performance of Constructed Facilities, 29(6):04014172.

[16]LacastaA, Morales-HernándezM, MurilloJ, et al., 2014. An optimized GPU implementation of a 2D free surface simulation model on unstructured meshes. Advances in Engineering Software, 78:1-15.

[17]LuoYZ, YangC, 2014. A vector-form hybrid particle-element method for modeling and nonlinear shell analysis of thin membranes exhibiting wrinkling. Journal of Zhejiang University-SCIENCE A (Applied Physics & Engineering), 15(5):331-350.

[18]NetoDM, OliveiraMC, MenezesLF, 2017. Surface smoothing procedures in computational contact mechanics. Archives of Computational Methods in Engineering, 24(1):37-87.

[19]SeitzA, PoppA, WallWA, 2015. A semi-smooth newton method for orthotropic plasticity and frictional contact at finite strains. Computer Methods in Applied Mechanics and Engineering, 285:228-254.

[20]SharmaA, MishraR, JainS, et al., 2018. Deformation behavior of single and multi-layered materials under impact loading. Thin-Walled Structures, 126:193-204.

[21]SimoJC, HughesTJR, 1998. Computational Inelasticity. Springer, New York, USA.

[22]TangJZ, ZhengYF, YangC, et al., 2020. Parallelized implementation of the finite particle method for explicit dynamics in GPU. Computer Modeling in Engineering & Sciences, 122(1):5-31.

[23]WangSQ, ZhangQ, JiSY, 2021. GPU-based parallel algorithm for super-quadric discrete element method and its applications for non-spherical granular flows. Advances in Engineering Software, 151:102931.

[24]WangW, ZhengYF, TangJZ, et al., 2021. A GPU-based parallel algorithm for 2D large deformation contact problems using the finite particle method. Computer Modeling in Engineering & Sciences, 129(2):595-626.

[25]WriggersP, 2006. Computational Contact Mechanics. Springer, Berlin, Germany.

[26]YangC, ShenYB, LuoYZ, 2014. An efficient numerical shape analysis for light weight membrane structures. Journal of Zhejiang University-SCIENCE A (Applied Physics & Engineering), 15(4):255-271.

[27]YuY, LuoYZ, 2009. Finite particle method for kinematically indeterminate bar assemblies. Journal of Zhejiang University-SCIENCE A, 10(5):669-676.

[28]ZangMY, GaoW, LeiZ, 2011. A contact algorithm for 3D discrete and finite element contact problems based on penalty function method. Computational Mechanics, 48(5):541-550.

[29]ZhengYF, WanHP, ZhangJY, et al., 2021. Local-coordinate representation for spatial revolute clearance joints based on a vector-form particle-element method. International Journal of Structural Stability and Dynamics, 21(7):2150093.