Shared Memory OpenMP Parallelization of Explicit MPM and Its Application to Hypervelocity Impact
Shared Memory OpenMP Parallelization of Explicit MPM and Its Application to Hypervelocity Impact
复制标题
DOI：
10.3970/cmes.2008.038.119
复制
发表时间：
2008-12
影响因子：
2.4
通讯作者：
P. Huang;Xiong Zhang;S. Ma;H. Wang
中科院分区：
文献类型：
--
作者：
P. Huang;Xiong Zhang;S. Ma;H. Wang
材料点法是质点网格法在固体力学中的推广。本文采用FORTRAN 95和OpenMP语言开发了一个主要用于求解冲击动力学问题的并行MPM程序。为了避免节点更新阶段的数据竞争，提出了两种并行算法：数组扩展法和区域分解法。在阵列展开法中，为节点变量创建二维辅助阵列。在更新所有线程中的网格节点后，组装辅助阵列以建立全局节点阵列。在区域分解方法中，背景网格被分解成一些均匀的块，每个线程处理一个块。通过共享变量交换相邻块的信息。在更新所有面片中的节点之后，将它们的节点变量组装起来以建立全局节点变量。数值试验表明，区域分解方法比阵列展开方法具有更好的并行可扩展性和更高的并行效率。因此，一个并行计算机程序，MPM3的开发基于区域分解方法。最后，利用MPM3300对13，542，030个粒子进行了大尺度模拟，得到了超高速碰撞碎片云的高分辨率结果。
The material point method (MPM) is an extension of particle-in-cell method to solid mechanics. A parallel MPM code is developed using FORTRAN 95 and OpenMP in this study, which is designed primarily for solving impact dynamic problems. Two parallel methods, the array expansion method and the domain decomposition method, are presented to avoid data races in the nodal update stage. In the array expansion method, two-dimensional auxiliary arrays are created for nodal variables. After updating grid nodes in all threads, the auxiliary arrays are assembled to establish the global nodal array. In the domain decomposition method, the background grid is decomposed into some uniform patches, and each thread deals with a patch. The information of neighbor patches is exchanged through shared variables. After updating nodes in all patches, their nodal variables are assembled to establish the global nodal variables. The numerical tests show that the domain decomposition method has much better parallel scalability and higher parallel efficiency than the array expansion method. Therefore, a parallel computer code, MPM3DMP, is developed based on the domain decomposition method. Finally, MPM3DMP is applied to a large-scale simulation with 13,542,030 particles for obtaining the high-resolution results of debris cloud in hypervelocity impact.