In this paper, a computational framework based on the Smoothed Particle Finite Element Method is developed to study the coupled seepage-deformation process in unsaturated porous media. Governing equations are derived from the balance laws of solid and fluid phases considering partial saturation effects in porous media. Moreover, an hourglass control method is implemented to avoid the rank-deficiency issue in SPFEM and the moving least squares approximation technique (MLS) is implemented to eliminate the pore pressure oscillations when the low-order triangle element is used. The proposed coupled SPFEM formulation is validated through four elastic examples and one elasto-plastic example. Good agreement with the numerical or analytical results reported in the literature is obtained. Further, the rainfallinduced slope failure is studied, in which a suction-dependent elasto-plastic Mohr-Coulomb model is adopted to take account of the suction effect in unsaturated soil. The evolution of the suction and soil deformation during the rainfall period and the whole slope failure process are obtained. It is demonstrated that the proposed method is a promising tool in numerical investigations of both the triggering mechanisms and post-failure behavior of the rainfall-induced slope failure.

A landslide is a common natural disaster that causes environmental damage, casualties and economic losses, which seriously affects the sustainable development of society. In geomechanics, it is one of the largest deformation problems. Herein, the GPU-accelerated explicit smoothed particle finite element method (eSPFEM) for large deformation analysis in geomechanics was developed on the CUDA platform based on high-performance computing using a self-designed eSPFEM program code. The eSPFEM combines the strain smoothing nodal integration techniques found in the particle finite element method (PFEM) framework, which allows for the use of low-order triangular elements without volume locking and avoids frequent information transfer and mapping errors between Gaussian points and particles in PFEM. A numerical simulation of slope instability using the eSPFEM and based on a strength reduction technique was conducted using various examples, including a cohesive homogeneous slope, a non-cohesive homogeneous slope, a non-homogeneous slope and a slope with a thin soft band. The calculation results show that the eSPFEM can be applied to slope stability analysis under different working conditions, simulating the entire process of slope instability initiation, sliding and reaccumulation, and obtaining reliable FOS values. A numerical simulation was conducted to analyse a landslide that occurred in the Zhangjiazhuang tunnel on the Lanzhou-Xinjiang high-speed railway line on 18 January 2016. A natural unsaturated soil slope, a soil slope with a high moisture content and a soil slope with a high moisture content subjected to an earthquake were analysed. The findings of this study are in good agreement with the actual slope failure conditions. The primary triggers identified for the landslide were heavy rainfall and earthquakes. The verification results indicate that the eSPFEM can effectively simulate an actual landslide case, showcasing high accuracy and applicability in simulating the large deformation behaviour of landslides.