GPU-accelerated vector-form particle-element method for 3D elastoplastic contact of structures

目的 结构的三维弹塑性接触问题通常包含强非线性, 且计算比较耗时。为解决这类问题, 本文提出基于图形处理器加速的有限质点法。 创新点 1. 发展基于有限质点法的六面体缩减积分单元; 2. 提出结构的三维并行接触算法。 方法 1. 发展基于有限质点法的六面体缩减积分单元, 并采用沙漏控制技术, 用于模拟结构的弹塑性行为; 2. 提出结构的三维并行接触算法, 将包含接触面的三维空间分解为立方体单元格, 仅在相邻单元格之间进行接触搜索, 并使用链式数据结构存储接触质点; 3. 通过基于图形处理器的并行计算技术对算法进行加速。 结论 1. 本文方法与有限元软件Abaqus/Explicit相比, 在总计算时间和接触计算时间上分别提升效率约80倍和340倍; 2. 本文方法的有效性和计算效率都得到了验证。 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.