Open Access

ARTICLE

Caveats on the Implementation of the Generalized Material Point Method

by O. Buzzi¹, D. M. Pedroso², A. Giacomini¹

1 Centre for Geotechnical and Materials Modelling, The University of Newcastle, Newcastle NSW 2308 Australia.
2 Golder Geomechanics Centre, School of Engineering, The University of Queensland, Brisbane QLD 4072 Australia.

Computer Modeling in Engineering & Sciences 2008, 31(2), 85-106. https://doi.org/10.3970/cmes.2008.031.085

Abstract

The material point method (MPM) is a numerical method for the solution of problems in continuum mechanics, including situations of large deformations. A generalization (GMPM) of this method was introduced by Bardenhagen and Kober (2004) in order to avoid some computational instabilities inherent to the original method (MPM). This generalization leads to a method more akin of the Petrov-Galerkin procedure. Although it is possible to find in the literature examples of the deduction and applications of the MPM/GMPM to specific problems, its detailed implementation is yet to be presented. Therefore, this paper attempts to describe all steps required for the explicit implementation of the material point method, including its generalization. Moreover, some caveats during the implementation are addressed. For example, the setting up of boundary conditions and the steps for the computation of values at nodes and material points are discussed. The influences of the time and space discretization are also verified, basing on numerical analyses. Two strategies for the update of stress, known as update stress first (USF) and update stress last (USL) are numerically investigated. It is shown that both the order for the computation of boundary conditions and the way that the grid values are extrapolated have high impact on the accuracy of the solution. The complete 3D algorithm is detailed and summarized in order to make easier the implementation of the GMPM/MPM.

---

Keywords

---