A semi‐implicit material point method based on fractional‐step method for saturated soil

• Published in: International journal for numerical and analytical methods in geomechanics Vol. 45; no. 10; pp. 1405 - 1436
• Main Authors: Kularathna, Shyamini, Liang, Weijian, Zhao, Tianchi, Chandra, Bodhinanda, Zhao, Jidong, Soga, Kenichi
• Format: Journal Article
• Language: English
• Published: Bognor Regis Wiley Subscription Services, Inc 01.07.2021
• Subjects: Algorithms; Computational fluid dynamics; Deformation; Exact solutions; Fluid flow; Formability; fractional‐step method; Hydrostatic pressure; Incompressibility; Incompressible flow; incompressible pore fluid; Interpolation; material point method; Oscillations; Permeability; Pore water pressure; Porous media; Pressure oscillations; saturated soil; Saturated soils; slope instability; Soil; Soil permeability; Splitting; Water pressure
• ISSN: 0363-9061; 1096-9853
• DOI: 10.1002/nag.3207

In this paper, a new formulation of material point method (MPM) to model coupled soil deformation and pore fluid flow problems is presented within the framework of the theory of porous media. The saturated porous medium is assumed to be consisting of incompressible pore fluid and deformable soil skeleton made up of incompressible solid grains. The main difference of the proposed MPM algorithm is the implicit treatment of pore‐water pressure which satisfies its incompressibility internal constraint. The resulting solid‐fluid coupled equations are solved by using a splitting algorithm based on the Chorin's projection method. The splitting algorithm helps to mitigate numerical instabilities at the incompressibility limit when equal‐order interpolation functions are used. The key strengths of the proposed semi‐implicit coupled MPM formulation is its capability to reduce pressure oscillations as well as to increase the time step size, which is independent of the fluid incremental strain level and the soil permeability. The proposed semi‐implicit MPM is validated by comparing the numerical results with the analytical solutions of several numerical tests, including 1D and 2D plane‐strain consolidation problems. To demonstrate the capability of the proposed method in simulating practical engineering problems involving large deformations, a hydraulic process leading to slope failure is studied, and the numerical result is validated by the monitored data.