Development of a stable two-phase contact MPM algorithm for saturated soil-structure interaction problems
share
The simulation of soil-structure interaction problems involving two-phase materials poses significant challenges in geotechnical engineering. These challenges arise due to differences in material stiffnesses, the interaction between multiple phases, high bulk modulus of pore fluid, and low permeability. The conventional explicit time integration scheme is limited by its conditional stability, necessitating small time step sizes and resulting in pressure oscillations under rapid loading conditions. To address these issues, we propose a stable two-phase contact algorithm within the Material Point Method (MPM) framework for soil-structure interaction problems. Our algorithm models the soil as a fully saturated porous media with incompressible pore fluid. We introduce three main advancements over conventional MPM methods. We employ Chorin's projection method to solve coupled formulations and reduce numerical oscillations. By implicitly handling a diffusion term, our algorithm permits larger stable time step sizes, independent of the bulk modulus and permeability of the pore fluid. Lastly, We integrate a rigid algorithm to model solid bodies accurately and a precise contact detection algorithm. We provide detailed formulations and time increment processes of the two-phase contact MPM algorithm. Furthermore, we compare the proposed algorithm with Finite Element Method (FEM) and explicit MPM to assess its accuracy and performance in simulating coupled hydro-mechanical problems. The two-phase contact algorithm offers a more stable and efficient approach to simulate soil-structure interaction problems.
READ FULL TEXT