Extended B-spline-based implicit material point method for saturated porous media

Abstract

The large deformation and fluidization process of a solid–fluid mixture includes significant changes to the temporal scale of the phenomena and the shape and properties of the mixed material. This paper presents an extended B-spline (EBS)-based implicit material point method (EBS-MPM) for the coupled hydromechanical analysis of saturated porous media to enhance the overall versatility of MPM in addressing such diverse phenomena. The proposed method accurately represents phenomena such as high-speed motion in both the quasi-static and dynamic states by employing a full formulation of coupled hydromechanical modeling. The weak imposition of boundary conditions based on Nitsche's method allows representing the boundary conditions independent of the relative position of the particles and computational grid. In addition, it enables dynamic changes in the boundary domain based on the deformation. The robustness of this boundary representation is reinforced using EBS basis functions, which prevent the degradation of the condition number of the system matrices regardless of the position of the boundary domain with respect to the computational grid. Furthermore, a stabilization method based on a variational multiscale method (VMS) approach is employed to provide the flexibility in choosing arbitrary basis functions for spatial discretization, facilitating the effective construction of EBS. Numerical examples including comparisons between a full formulation and a simplified formulation are presented to demonstrate the performance of the developed method under various boundary conditions and loading states across different time scales.