A Stabilised Semi-Implicit Double-Point Material Point Method for Soil-Water Coupled Problems

Mian Xie,Pedro Navas,Susana Lopez-Querol

2024-01-22

Abstract:A semi-implicit two-phase double-point Material Point Method (MPM) formulation, based on the incremental fractional-step method to model large deformation geotechnical problems has been derived. The semi-implicit formulation has two advantages compared with the explicit approach: the time step is independent of the water phase, and the pore pressure field is more stable. The semi-implicit MPM models based on the incremental fractional-step method available in the literature consist of modelling the soil and water mixture using a single set of material points only, in order to save computational time. In this study, we further derive this formulation with two sets of material points to represent the soil and water phases separately. The stress oscillations that are frequently found in the water and soil phases are stabilised with this approach. A new stabilisation method is developed based on the modified F-bar method. The proposed method is validated with two numerical examples under small and large deformations, respectively. After that, Nor-Sand constitutive soil model is used to simulate landslides. Numerical examples show an excellent performance of the proposed coupled MPM and the stabilisation method. The formulation with two sets of material points yields significantly different but more reliable results in the landslides analysis, compared with the single-point approach. Additionally, this research shows that the additional computational cost caused by the additional water material points is acceptable. Therefore, it is recommended to use two sets of material points for certain large deformation geotechnical problems.

Numerical Analysis

What problem does this paper attempt to address?

This paper aims to address several key challenges in the soil - water coupling problem. Specifically, the paper attempts to improve the performance of existing methods in simulating large - deformation geotechnical problems by proposing a stabilized semi - implicit two - point Material Point Method (MPM). This method is based on the incremental fractional - step method and mainly solves the following problems:

1. **Time step independent of the water phase**: Compared with explicit methods, the proposed semi - implicit method makes the time step no longer affected by the water phase, so that a larger time step can be used for simulation, improving computational efficiency.

2. **Stability of the pore pressure field**: The semi - implicit method can provide a more stable pore pressure field and reduce the stress oscillation problems common in numerical simulations.

3. **Introduction of the two - point method**: Different from the single - point method, the method proposed in the paper uses two sets of material points to represent the soil and water phases respectively. This not only improves the accuracy of the model, but also can better simulate the relative motion between the soil and water phases, such as seepage, erosion and submarine landslides.

4. **Improvement of numerical stability**: In order to further improve the stability of the model, the paper develops a new stabilization technique based on the modified F - bar method, which effectively eliminates the volumetric locking problem, especially in the case of very large deformations.

5. **Application of advanced soil constitutive models**: The paper also implements the Nor - Sand constitutive model to show the different results of the two - point method and the single - point method in analyzing complex geological problems such as landslides, and emphasizes the importance of using the two - point MPM.

Through these improvements, the method proposed in the paper performs well in dealing with large - deformation geotechnical problems. Especially in simulating phenomena such as landslides, compared with the single - point method, the two - point method provides more reliable results, and the additional computational cost is acceptable. Therefore, for some specific large - deformation geotechnical problems, the two - point Material Point Method is recommended.