P. Huang1,2, X. Zhang1,3, S. Ma1, H.K. Wang1
CMES-Computer Modeling in Engineering & Sciences, Vol.38, No.2, pp. 119-148, 2008, DOI:10.3970/cmes.2008.038.119
Abstract 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… More >