Premium
A semi‐implicit material point method based on fractional‐step method for saturated soil
Author(s) -
Kularathna Shyamini,
Liang Weijian,
Zhao Tianchi,
Chandra Bodhinanda,
Zhao Jidong,
Soga Kenichi
Publication year - 2021
Publication title -
international journal for numerical and analytical methods in geomechanics
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.419
H-Index - 91
eISSN - 1096-9853
pISSN - 0363-9061
DOI - 10.1002/nag.3207
Subject(s) - compressibility , consolidation (business) , material point method , projection method , interpolation (computer graphics) , permeability (electromagnetism) , porosity , porous medium , mathematics , numerical analysis , plane stress , mechanics , pore water pressure , geotechnical engineering , finite element method , mathematical analysis , mathematical optimization , computer science , geology , dykstra's projection algorithm , engineering , structural engineering , physics , computer graphics (images) , membrane , business , biology , genetics , accounting , animation
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.