This paper presents a method for analyzing slope stability in anisotropic and heterogeneous clay using a strength reduction finite element method (SRFEM) integrated with the level set method (LSM). Anisotropy refers to the inherent anisotropy in the clay's strength, while heterogeneity describes the spatial variability in strength parameters. The static LSM uses a zero level set function to model heterogeneous clay slopes. The method is validated through undrained slope stability analyses on different types of anisotropic clay and heterogeneous fields, showing its effectiveness in modeling anisotropic shear strength and capturing the characteristics of heterogeneous regions. The results indicate that the proposed method accurately predicts factors of safety and slip surfaces across various soil conditions, accounting for both anisotropic and heterogeneous characteristics.

The structural characteristics of soil-rock mixture (SRM) slopes, including the content, shape, size, and spatial distribution of rock blocks, can significantly influence their failure mechanisms and factor of safety (FOS). Defining the structural characteristics of SRM slopes for stability analysis remains challenging. This study proposes a method for establishing random models and evaluating the statistical properties of the FOS values of SRM slopes. Accordingly, the SRM slope models were constructed by considering the random properties of the shape, size, and spatial distribution of rock blocks in the slope domain. A slope failure criterion based on energy changes and the combined subroutines of USDFLD and URDFIL was implemented in the ABAQUS finite element software to determine the FOS values of the SRM slopes. Monte Carlo simulations were performed to assess the statistical properties of the FOS for random SRM slopes varying rock block properties. The results indicated that when the rock block content was greater than 30%, the stability of SRM slopes considerably increased. For a rock block content of 40%, the effect of rock block size on the SRM slope stability followed two different trends: the mean FOS value tended to decline and subsequently increased as rock block size increased. However, this trend was not observed on SRM slopes with a 30% rock block content. Besides, the dispersion of the FOS values gradually increased with increasing rock content and rock block size. Furthermore, the soil-rock interface strength affected the stability and failure mechanism of SRM slopes. These findings enhance comprehension of the SRM slope stability assessment and demonstrate improved accuracy in predicting and mitigating damage.

The determination of the factor of safety (FoS) of slopes during seismic excitation can be complex if the relevant effects of pore water pressure accumulation, nonlinear material response and variable shear strength are duly accounted for. A rational two-step approach to tackle this task based on a hydro-mechanically coupled dynamic simulation and finite element limit analyses is henceforth introduced. To ensure accurate transfer of the hydro-mechanical soil state, a mapping concept is presented, accounting for spatial distributions of stresses, excess pore water pressures, inertial forces and shear strength. The proposed approach is compared to limit equilibrium method (LEM) for the case of a large-scale water-saturated open cast mine slope subjected to seismic loading. In comparison with LEM, the new approach to assess seismic slope stability proves to be simpler in its implementation and straightforward, which could be an important asset for practitioners.

The rock-soil mass, subjected to complex and lengthy geological processes, exhibits heterogeneity which induces variations in mechanical properties, thereby affecting the overall stability of slopes. In this paper, a novel numerical model that incorporates the Weibull distribution function into the meshless numerical manifold method based on the strength reduction method (MNMM-SRM) to account for the slope soils heterogeneity and their influence on the factor of safety (Fs) and the critical sliding surface (CSS). Initially, the Weibull distribution is introduced into the MNMM-SRM model based on the complementary theory of subspace tracking, addressing the issue of multiple yield surface corners in the Mohr-Coulomb framework while simultaneously considering the heterogeneous nature of rock and soil formations. Subsequently, an intelligent method based on unsupervised learning is proposed to obtain reasonable CSS, utilizing the total displacement field at slope nodes and the equivalent plastic strain field as input variables. The results serve as criteria for terminating the strength reduction in the MNMM-SRM. The applicability of this method is verified through three typical examples, demonstrating its potential for widespread application in the assessment of heterogeneous slope stability.

Ensuring the stability of slopes is critical to the safe operation of geotechnical engineering. Evaluating slope stability to minimize geologic risks induced by destabilization is significant in reducing casualties and property damage. A conventional, single-coefficient strength reduction method is widely applied in slope stability analyses, but this method ignores the attenuation degree of different parameters in the slope destabilization. A new double-strength reduction method considering different contributions of the mechanics' parameters is proposed in this study for evaluating the stability of nonhomogeneous slope. First, the role of each mechanic's parameters in the slope destabilization was investigated theoretically and numerically using ABAQUS software 2022. The results indicate that the effect of elasticity (E), Poisson's ratio (v), and soil gravity (gamma) on the evolution of factor of safety (FOS) are insignificant and can be neglected compared with cohesive force (c), and angle of internal friction (phi). Next, an improved method was constructed to correlate the FOS with cohesive force (c) and the angle of internal friction (phi). Then, a numerical method was constructed based on the computation of the mathematical-mechanical relationship between FOS and the mechanical parameters, and the stability of slope is estimation based on the Mohr-Coulomb yield criterion. Finally, the double-strength reduction coefficient method proposed in this study, the limit equilibrium method, and the traditional finite element strength reduction coefficient method were applied to nonhomogeneous slopes and slopes containing a soft underlying layer for comparison, and the difference between them was within the range of +/- 5%. The results indicate that both the limit equilibrium method and the traditional finite element strength reduction method tend to overestimate the FOS of intricate slopes compared with the evaluated method proposed in this study. This parallel comparison serves to validate the accuracy of the double-strength reduction method proposed in the present study. Further, based on the proposed method, the relationship between slope stability and slope displacement is established, which provides a theoretical basis for the safety assessment of slope engineering.

Numerous manually excavated loess caves are present within a cultural relic protection zone in the northwestern region of China. The collapse of these caves frequently leads to the cracking, tilting, and even collapse of ancient buildings above, posing a severe threat to the safety of cultural architectural relics. Investigating the stability and characteristics of deformation and failure in loess caves is essential for effectively reinforcing and protecting cultural relics. A twodimensional model of a loess underground cavern was developed using OptumG2. The stability and modes of deformation and failure in the underground cavern were analyzed through the augmentation of soil gravity and the strength reduction method. This analysis determined the cavern's safety factor, force, deformation and damage mode, and the plastic zone's progression. Numerical simulations analyzed the force characteristics of the support structure under different stress release ratios. The findings revealed that, with the implementation of an anchor rod concrete lining support scheme, the most probable failure mode is a shear failure, initiating at the arch foot. The ground's stress release rate does not influence the safety factor of the cavern but rather the material, design, and strength of the support structure. However, the magnitude of the internal forces acting on the supporting structure by the soil in the cavern is related to the degree of ground stress release. When applied during significant stress release, support structures may experience reduced internal forces, albeit with more substantial stratum displacement; opting for an appropriate stress release when applying support structures is crucial for achieving optimal stratum displacement and lining internal forces.

Among various available methods for slope analysis, the limit equilibrium method is very popular because of its simple concepts. The limit analysis method and the finite element method (FEM) also can perform stability analysis of a slope. Increasing computing power and the easy accessibility of inexpensive numerical modeling codes have made the finite element method a very attractive tool for the practical assessment of slope stability. The present study reports the results of slope stability analysis of a few problems analyzed using a developed program utilizing FEM. This program employs a strength reduction technique based on FEM. Mohr-Coulomb strength criterion of soil is used for predicting the stress state, while the viscoplastic algorithm is used for stress redistribution. Non-convergence of the algorithm to achieve the desired equilibrium of all forces in the system is adopted as a marker of slope failure. Further, to put the proposed method to the test, a few examples from the literature are analyzed using the developed program. The example problems cover a homogenous slope with water loading, an inclined layered slope, and a staged embankment subjected to different forms of loading including earthquake forces, pore water pressure, external water pressure, etc. The results of each analysis are compared with other researchers work, and it is found that the obtained results are in good agreement. Deformed mesh, equivalent viscoplastic strain contour plots, and failure function contour plots are used for illustrating the failure state.

Fissures form the channel for rainwater infiltration, which accelerate the infiltration of rainwater into slope bodies, hence its important impact on the seepage field and stability of the slope. In this paper, taking one landslide of Liang-Wan freeway as the research object, firstly, the equivalent permeability coefficient method is used to homogenize the fissured soil. Then considering the boundary conditions of rainfall infiltration and groundwater level, a fluid-structure coupling model is established based on saturated-unsaturated seepage theory, and evolution characteristics of seepage, displacement and stress of the slope are studied. Based on these, the slope stability coefficient is determined. The results show that the rising rate of pore water pressure and volume water content of topsoil increases when multi-fissure seepage is considered, and the pore water velocity is larger in the local seepage range of fissures. With the increase of buried depth, the closer to groundwater level, the influence of multi-fissure seepage gradually weakens. The theoretical calculation results of slope displacement are more consistent with the field monitoring results. With the increase of rainfall time, the stability coefficient of slope decreases gradually, and the rate and range of decrease are greater.

There are many factors affecting the stability of open-pit mine slopes, among which slopes with soft interlayers have become an important factor inducing deformation and instability due to their poor mechanical properties. In this paper, for the destabilization of a slope with soft interlayer under rainfall infiltration in an open-pit mine in Zhejiang Province of China during the rainy season, a sudden change of displacement at the monitoring point and the formation of a continuous plastic deformation of the slope are proposed as the criteria for slope destabilization through the establishment of a rainwater seepage-stress coupling model combined with the strength reduction method, and by using COMSOL Multiphysics finite element numerical simulation software to establish a three-dimensional numerical model based on the actual mining slope conditions. The safety principle is introduced through the strength reduction method to analyze the stability of the slope with soft interlayer under the coupling effect of seepage and stress. The safety coefficients of the slope with soft interlayer are calculated by using Numerical methods under the coupling effect of seepage and stress. The influence of rainfall intensity and duration on the stability or safety factor of the slope with soft interlayers was studied through analysis of stress, saturation, displacement, and pore pressure evolution. The mechanism of rainfall affecting slope stability was investigated, and the findings were validated using an engineering case study of a slope with a soft interlayer. The findings show that rainfall intensity is the main factor affecting slope stability in the open-pit mine slope with soft interlayer. The higher the rainfall intensity, the faster the shallow soil forms a saturation zone, whereas soft interlayers speed up the process and endanger slope stability. The slope prevention and control technology of the prestressed anchor cable (rod) framework was proposed based on the slope management construction conditions of the mine, and the effectiveness of the measure was verified through on-site industrial tests. The research findings provide a reference for preventing and controlling slope with soft interlayers in open-pit mines under similar conditions.

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.