Fully implicit, stabilised MPM simulation of large-deformation problems in two-phase elastoplastic geomaterials

The Material Point Method (MPM) has been gaining increasing popularity as an appropriate approach to the solution of coupled hydro-mechanical problems involving large deformations. This study extends the implicit GIMP-patch method for coupled poroelastic problems recently proposed by Zheng et al. (2021b) to tackle large-deformation problems in (nearly) isochoric elastoplastic geomaterials, particularly by remedying the numerical inaccuracies caused by volumetric locking, such as spurious stress oscillations and an excessively stiff overall response of the system at hand. To overcome these difficulties in two-phase coupled analyses, the B¯ approach of Hughes (1980) is incorporated into an existing version of the implicit GIMP-patch method. Details regarding the formulation and implementation of the proposed method are provided, while several benchmark problems are numerically analysed to evaluate its performance in the presence of elastoplastic behaviour. Particular emphasis is placed on (i) mitigating effective stress oscillations and (ii) solving several two-phase, coupled, large deformation geotechnical problems. The numerical results confirm the suitability of the implicit B¯ GIMP-patch method for the solution of geotechnical problems spanning weak to strong hydro-mechanical coupling and small to large deformations.