This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is e...This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is employed to analyze the stability of 3D anisotropic soil slopes.The accuracy of the proposed method is first verified against the data in the literature.We then simulate the 3D soil slope with a straight slope surface and the convex and concave slope surfaces with a 90turning corner to study the 3D effect on slope stability and the failure mechanism under anisotropy conditions.Based on our numerical results,the end effect significantly impacts the failure mechanism and safety factor.Anisotropy degree notably affects the safety factor,with higher degrees leading to deeper landslides.For concave slopes,they can be approximated by straight slopes with suitable boundary conditions to assess their stability.Furthermore,a case study of the Saint-Alban test embankment A in Quebec,Canada,is provided to demonstrate the applicability of the proposed FE model.展开更多
Sudden and unforeseen seismic failures of coal mine overburden(OB)dump slopes interrupt mining operations,cause loss of lives and delay the production of coal.Consideration of the spatial heterogeneity of OB dump mate...Sudden and unforeseen seismic failures of coal mine overburden(OB)dump slopes interrupt mining operations,cause loss of lives and delay the production of coal.Consideration of the spatial heterogeneity of OB dump materials is imperative for an adequate evaluation of the seismic stability of OB dump slopes.In this study,pseudo-static seismic stability analyses are carried out for an OB dump slope by considering the material parameters obtained from an insitu field investigation.Spatial heterogeneity is simulated through use of the random finite element method(RFEM)and the random limit equilibrium method(RLEM)and a comparative study is presented.Combinations of horizontal and vertical spatial correlation lengths were considered for simulating isotropic and anisotropic random fields within the OB dump slope.Seismic performances of the slope have been reported through the probability of failure and reliability index.It was observed that the RLEM approach overestimates failure probability(P_(f))by considering seismic stability with spatial heterogeneity.The P_(f)was observed to increase with an increase in the coefficient of variation of friction angle of the dump materials.Further,it was inferred that the RLEM approach may not be adequately applicable for assessing the seismic stability of an OB dump slope for a horizontal seismic coefficient that is more than or equal to 0.1.展开更多
Volcanic terrains exhibit a complex structure of pyroclastic deposits interspersed with sedimentary processes,resulting in irregular lithological sequences that lack lateral continuity and distinct stratigraphic patte...Volcanic terrains exhibit a complex structure of pyroclastic deposits interspersed with sedimentary processes,resulting in irregular lithological sequences that lack lateral continuity and distinct stratigraphic patterns.This complexity poses significant challenges for slope stability analysis,requiring the development of specialized techniques to address these issues.This research presents a numerical methodology that incorporates spatial variability,nonlinear material characterization,and probabilistic analysis using a Monte Carlo framework to address this issue.The heterogeneous structure is represented by randomly assigning different lithotypes across the slope,while maintaining predefined global proportions.This contrasts with the more common approach of applying probabilistic variability to mechanical parameters within a homogeneous slope model.The material behavior is defined using complex nonlinear failure criteria,such as the Hoek-Brown model and a parabolic model with collapse,both implemented through linearization techniques.The Discontinuity Layout Optimization(DLO)method,a novel numerical approach based on limit analysis,is employed to efficiently incorporate these advances and compute the factor of safety of the slope.Within this framework,the Monte Carlo procedure is used to assess slope stability by conducting a large number of simulations,each with a different lithotype distribution.Based on the results,a hybrid method is proposed that combines probabilistic modeling with deterministic design principles for the slope stability assessment.As a case study,the methodology is applied to a 20-m-high vertical slope composed of three lithotypes(altered scoria,welded scoria,and basalt)randomly distributed in proportions of 15%,60%,and 25%,respectively.The results show convergence of mean values after approximately 400 simulations and highlight the significant influence of spatial heterogeneity,with variations of the factor of safety between 5 and 12 in 85%of cases.They also reveal non-circular and mid-slope failure wedges not captured by traditional stability methods.Finally,an equivalent normal probability distribution is proposed as a reliable approximation of the factor of safety for use in risk analysis and engineering decision-making.展开更多
A method combining finite difference method(FDM)and k-means clustering algorithm which can determine the threshold of rock bridge generation is proposed.Jointed slope models with different joint coalescence coefficien...A method combining finite difference method(FDM)and k-means clustering algorithm which can determine the threshold of rock bridge generation is proposed.Jointed slope models with different joint coalescence coefficients(k)are constructed based on FDM.The rock bridge area was divided through k-means algorithm and the optimal number of clusters was determined by sum of squared errors(SSE)and elbow method.The influence of maximum principal stress and stress change rate as clustering indexes on the clustering results of rock bridges was compared by using Euclidean distance.The results show that using stress change rate as clustering index is more effective.When the joint coalescence coefficient is less than 0.6,there is no significant stress concentration in the middle area of adjacent joints,that is,no generation of rock bridge.In addition,the range of rock bridge is affected by the coalescence coefficient(k),the relative position of joints and the parameters of weak interlayer.展开更多
Block-flexure toppling constitutes the predominant form of toppling failure in rock slopes.Although it has been extensively studied,the current theoretical models are often oversimplified by treating rock layers as ri...Block-flexure toppling constitutes the predominant form of toppling failure in rock slopes.Although it has been extensively studied,the current theoretical models are often oversimplified by treating rock layers as rigid bodies that diverge from actual conditions.The proposed Equivalent Deformation Compatibility Method(EDCM)offers a fresh approach to assess the stability of rock slopes prone to block-flexure toppling.EDCM posits that blocky rock layers,with their inability to withstand significant bending and role in merely transferring forces,can be modeled as intact layers with a reduced modulus.The method simplifies the complex issue of analyzing discrete and continuous rock layers to the study of layered soft and hard rock,establishing deformation compatibility equations subsequently.Validation of the EDCM was achieved through numerical models,physical model testing,and application to an actual slope.The factor of safety(FS)for slopes corresponds with the results from both models and the actual slope,demonstrating the method's applicability for evaluating susceptibility to block-flexure toppling.When applying the EDCM,it is advised to set the elastic modulus reduction coefficient for blocky layers at a value below 0.1.展开更多
Steep bedding slopes are widely distributed in Southwestern China’s mountainous regions and have complex seismic responses and instability risks,causing casualties and property losses.Considering the high-seismic-int...Steep bedding slopes are widely distributed in Southwestern China’s mountainous regions and have complex seismic responses and instability risks,causing casualties and property losses.Considering the high-seismic-intensity environment,the dynamic failure evolution and instability mechanism of high-steep bedding slopes are simulated via the discrete element method and shaking table test.The dynamic response characteristics and cumulative failure effects of slopes subjected to continuous ground motion are investigated.The results show that the dynamic response characteristics of slopes under continuous earthquakes are influenced by geological and topographic conditions.Elevation has a distinct impact on both the slope interior and surface,with amplification effects more pronounced on the surface.The weak interlayers have different influences on the dynamic amplification effect of slopes.Weak interlayers have dynamic magnification effects on the slope surface at relative elevations of 0-0.33 and 0.82-1.0 but have weakening effects between 0.33 and 0.82.Moreover,the weak interlayers also have controlling effects on the dynamic instability mode of slopes.The characteristics of intergranular contact failure,fracture propagation,and displacement distribution are analyzed to reveal the dynamic failure evolution and instability mechanism through the discrete-element model.The dynamic instability process of slopes includes three stages:fracture initiation(0-0.2g),fracture expansion(0.2g-0.3g),and sliding instability(0.3g-0.6g).This work can provide a valuable reference for the seismic stability and reinforcement of complex slopes.展开更多
The stability of the tunnel portal slope is crucial for ensuring safe tunnel construction.Thus,a sound stability evaluation is of significance.Given the unique geological characteristics of tunnel portal slopes,it is ...The stability of the tunnel portal slope is crucial for ensuring safe tunnel construction.Thus,a sound stability evaluation is of significance.Given the unique geological characteristics of tunnel portal slopes,it is necessary to establish a specific evaluation indicator system that differs from those used for ordinary slopes.Based on the unascertained measure method,uncertainties in the indicator are addressed by introducing the left and right half cloud asymmetric cloud model to optimize the linear membership function.The subjectivity of confidence criterion level identification is also improved by using the Euclidean distance method.Thus,a stability evaluation model for the tunnel portal slope is established based on the improved unascertained measure method.Finally,using the collected tunnel portal slope data,the results of four evaluation methods are compared with the safety factor levels.The evaluation methods include the traditional unascertained measure method,the method improved by using the left and right half cloud asymmetric cloud model,the method improved by using the Euclidean distance method,and the method improved by using both the left and right half cloud asymmetric cloud model and the Euclidean distance method.The results show that the accuracy rates of these four methods are 50%,55%,85%,and 90%,respectively.Among them,the joint improvement method has the slightest deviation,with only one level,while the other three methods had deviations of two levels.This result verifies the stability and effectiveness of the joint improvement method,providing a reference for tunnel portal slope stability evaluation.展开更多
This study proposes an alternative calculation mode for stresses on the slip surface(SS).The calculation of the normal stress(NS)on the SS involves examining its composition and expanding its unknown using the Taylor ...This study proposes an alternative calculation mode for stresses on the slip surface(SS).The calculation of the normal stress(NS)on the SS involves examining its composition and expanding its unknown using the Taylor series.This expansion enables the reasonable construction of a function describing the NS on the SS.Additionally,by directly incorporating the nonlinear Generalized Hoke-Brown(GHB)strength criterion and utilizing the slope factor of safety(FOS)definition,a function of the shear stress on the SS is derived.This function considers the mutual feedback mechanism between the NS and strength parameters of the SS.The stress constraints conditions are then introduced at both ends of the SS based on the spatial stress relation of one point.Determining the slope FOS and stress solution for the SS involves considering the mechanical equilibrium conditions and the stress constraint conditions satisfied by the sliding body.The proposed approach successfully simulates the tension-shear stress zone near the slope top and provides an intuitive description of the concentration effect of compression-shear stress of the SS near the slope toe.Furthermore,compared to other methods,the present method demonstrates superior processing capabilities for the embedded nonlinear GHB strength criterion.展开更多
The vibration characteristics and dynamic responses of rock and soil under seismic load can be estimated with dynamic finite element method （DFEM）. Combining with the DFEM, the vector sum analysis method （VSAM） is...The vibration characteristics and dynamic responses of rock and soil under seismic load can be estimated with dynamic finite element method （DFEM）. Combining with the DFEM, the vector sum analysis method （VSAM） is employed in seismic stability analysis of a slope in this paper. Different from other conventional methods, the VSAM is proposed based on the vector characteristic of force and current stress state of the slope. The dynamic stress state of the slope at any moment under seismic load can he obtained by the DFEM, thus the factor of safety of the slope at any moment during earthquake can be easily obtained with the VSAM in consideration of the DFEM. Then, the global stability of the slope can be estimated on the basis of time-history curve of factor of safety and reliability theory. The VSAM is applied to a homogeneous slope under seismic load. The factor of safety of the slope is 1.30 under gravity only and the dynamic factor of safety under seismic load is 1.21. The calculating results show that the dynamic characteristics and stability state of the slope with input ground motion can be actually analyzed. It is believed that the VSAM is a feasible and practical approach to estimate the dynamic stability of slopes under seismic load.展开更多
The attenuation factor or quality factor(Q-factor or Q) has been used to measure the energy attenuation of seismic waves propagating in underground media. Many methods are used to estimate the Q-factor. We propose a m...The attenuation factor or quality factor(Q-factor or Q) has been used to measure the energy attenuation of seismic waves propagating in underground media. Many methods are used to estimate the Q-factor. We propose a method to calculate the Q-factor based on the prestack Q-factor inversion and the generalized S-transform. The proposed method specifies a standard primary wavelet and calculates the cumulative Q-factors; then, it finds the interlaminar Q-factors using the relation between Q and offset(QVO) and the Dix formula. The proposed method is alternative to methods that calculate interlaminar Q-factors after horizon picking. Because the frequency spectrum of each horizon can be extracted continuously on a 2D time–frequency spectrum, the method is called the continuous spectral ratio slope(CSRS) method. Compared with the other Q-inversion methods, the method offers nearly effortless computations and stability, and has mathematical and physical significance. We use numerical modeling to verify the feasibility of the method and apply it to real data from an oilfield in Ahdeb, Iraq. The results suggest that the resolution and spatial stability of the Q-profile are optimal and contain abundant interlaminar information that is extremely helpful in making lithology and fluid predictions.展开更多
The core of strength reduction method(SRM) involves finding a critical strength curve that happens to make the slope globally fail and a definition of factor of safety(FOS). A new double reduction method, including a ...The core of strength reduction method(SRM) involves finding a critical strength curve that happens to make the slope globally fail and a definition of factor of safety(FOS). A new double reduction method, including a detailed calculation procedure and a definition of FOS for slope stability was developed based on the understanding of SRM. When constructing the new definition of FOS, efforts were made to make sure that it has concise physical meanings and fully reflects the shear strength of the slope. Two examples, slopes A and B with the slope angles of 63° and 34° respectively, were given to verify the method presented. It is found that, for these two slopes, the FOSs from original strength reduction method are respectively 1.5% and 38% higher than those from double reduction method. It is also found that the double reduction method predicts a deeper potential slide line and a larger slide mass. These results show that on one hand, the double reduction method is comparative to the traditional methods and is reasonable, and on the other hand, the original strength reduction method may overestimate the safety of a slope. The method presented is advised to be considered as an additional option in the practical slope stability evaluations although more useful experience is required.展开更多
A new method （kinetic vector method, KVM） is presented for analyzing the dynamic stability of wedge in rock slope. The dynamic analysis is carried out based on three dimensional distinct element code （3DEC）, and t...A new method （kinetic vector method, KVM） is presented for analyzing the dynamic stability of wedge in rock slope. The dynamic analysis is carried out based on three dimensional distinct element code （3DEC）, and the kinetic inertial force of the wedge under seismic loading can be obtained via calculating the net vectorial nodal force of the finite difference grid. Then, the factor of safety （FOS） of the wedge can be calculated based on limit equilibrium method （LEM） at each dynamic analysis step, therefore time series of the FOS for whole earthquake process can be obtained. For the purpose of evaluating the entire dynamic stability of the wedge, dynamic factor of safety （DFOS） is proposed and defined as a numerical value corresponding with a given rate of probability guarantee based on reliability theory. Consequently, the KVM inherits the merits of the LEM and also has fully nonlinear dynamic analysis capabilities, and the feasibility and correctness of the KVM are tested by an example given by Hoek and Bray （1981）. Finally, a rock slope case in Wenchuan Earthquake regions of China is presented to verify the engineering practicability of the KVM, and the results matched the actual situation well.展开更多
Overhanging rock slopes(steeper than 90°) are typically avoided in rock engineering design, particularly where the scale of the slope exceeds the scale of fracturing present in the rock mass. This paper highlight...Overhanging rock slopes(steeper than 90°) are typically avoided in rock engineering design, particularly where the scale of the slope exceeds the scale of fracturing present in the rock mass. This paper highlights an integrated approach of designing overhanging rock slopes where the relative dimensions of the slope exceed the scale of fracturing and the rock mass failure needs to be considered rather than kinematic release of individual blocks. The key to the method is a simplified limit equilibrium(LE) tool that was used for the support design and analysis of a multi-faceted overhanging rock slope. The overhanging slopes required complex geometries with constantly changing orientations. The overhanging rock varied in height from 30 m to 66 m. Geomechanical modelling combined with discrete fracture network(DFN)representation of the rock mass was used to validate the rock mass strength assumptions and the failure mechanism assumed in the LE model. The advantage of the simplified LE method is that buttress and support design iterations(along with sensitivity analysis of design parameters) can be completed for various cross-sections along the proposed overhanging rock sections in an efficient manner, compared to the more time-intensive, sophisticated methods that were used for the initial validation. The method described presents the development of this design tool and assumptions made for a specific overhanging rock slope design. Other locations will have different geological conditions that can control the potential behaviour of rock slopes, however, the approach presented can be applied as a general guiding design principle for overhanging rock cut slope.展开更多
Large-scale slope destabilization could be aggravated due to swift urbanization and ever-rising demands of geoengineering projects such as dams,tunnels,bridges and widening roads.National Highway-58 connects Delhi to ...Large-scale slope destabilization could be aggravated due to swift urbanization and ever-rising demands of geoengineering projects such as dams,tunnels,bridges and widening roads.National Highway-58 connects Delhi to Badrinath in India,which passes through complex geomorphological and geological terrain and often encounters cut slopes susceptible to slope failures.In the present investigation,a detailed geotechnical appraisal is conducted along the road cut slopes from Rishikesh to Devprayag in the Himalayas.Twenty vulnerable road cut slopes were demarcated for detailed slope stability analysis using Phase2D finite element modeling simulator.Nonlinear generalized Hoek-Brown(GHB)criterion was adopted for stability analyses.Out of 20 slopes,five slopes(S6,S7,S18,S19 and S20)are unstable with factor of safety(FoS)less than or equal to 1,and thus needs immediate attention.The FoS values of four slopes(S2,S9,S13 and S17)lie between 1 and 1.3,i.e.marginally stable,and slopes S1,S3,S4,S5,S8,S10,Sll,S12,S14,S15 and S16 are stable.Mohr-Coulomb(MC)criterion was also adopted to compare the slope stability analysis with GHB criterion.The FoS calculated from GHB criterion is close to that using MC criterion for lower values of FoS whereas for higher values,the difference is marked.For the jointed rock in the Himalayan region,the nonlinear GHB criterion gives better results as compared to MC criterion and matches with the prevailing field conditions.Accordingly,some suggestions are proposed to strengthen the stability of cut slopes.展开更多
The elliptic mild slope equation is used to simulate linear wave propagation over variable sea bed topography with mild slopes. The governing equation is discretized by the finite difference method. Based on the BI-CG...The elliptic mild slope equation is used to simulate linear wave propagation over variable sea bed topography with mild slopes. The governing equation is discretized by the finite difference method. Based on the BI-CGSTAB technique, an attractive variant bf BI-Conjugate Gradients (BI-CG) method, the obtained linear algebraic system of equations is solved. Numerical experiments show that the BI-CGSTAB method is efficient for solving the elliptic mild slope equation. The results obtained by the BI-CGSTAB-Based method are much the same as those obtained by other authors with different solution methods, but the convergence rate is much faster than that of other methods.展开更多
Although slice methods are simple and effective slope stability analysis approaches,they are statically indeterminate.Several modifications of the slice method,such as the Spencer,MorgensternPrice,and Chen-Morgenstern...Although slice methods are simple and effective slope stability analysis approaches,they are statically indeterminate.Several modifications of the slice method,such as the Spencer,MorgensternPrice,and Chen-Morgenstern methods,are statically determinate and solvable as they assume the inter-slice force inclination angle;however,there is a small gap between the assumptions and actual landslide stability analysis.Through reasonable theoretical analysis,the Su slice method provides a reliable approach for determining the inter-slice force inclination angle that can be used in slice analysis to accurately analyse,calculate,and evaluate the stability of landslides.However,the Su slice method requires further research and analysis,especially in terms of the parameter values sinλbiandρ.In this study,we investigated more accurate methods for calculating the parameters sinλbiandρ.In addition,an adjustment coefficient(μ)was introduced to improve the solution method for the inter-slice force inclination angle.The inter-slice force inclination and safety factors of three landslides with arc-shaped slip surfaces and one landslide with a polyline-shaped slip surface were analysed and compared using the different slice methods.The improved inter-slice force inclination not only satisfies the calculation of static force equilibrium condition but also satisfies the calculation of both the force and moment equilibrium conditions.The improved method for calculating inter-slice force inclination presented the best correlation.The safety factors calculated using the improved Su slice method were close to those obtained using numerical simulations and the Morgenstern-Price method.Despite negligible differences among the safety factors calculated using the Su slice,improved Su slice,and M-P methods,the accuracy of the improved Su slice method was better than the M-P method in terms of inter-slice force inclination angles which can be useful to improve protection engineering design.展开更多
One of the critical aspects in mine design is slope stability analysis and the determination of stable slopes. In the Chador- Malu iron ore mine, one of the most important iron ore mines in central Iran, it was consid...One of the critical aspects in mine design is slope stability analysis and the determination of stable slopes. In the Chador- Malu iron ore mine, one of the most important iron ore mines in central Iran, it was considered vital to perform a comprehensive slope stability analysis. At first, we divided the existing rock hosting pit into six zones and a geotechnical map was prepared. Then, the value of MRMR （Mining Rock Mass Rating） was determined for each zone. Owing to the fact that the Chador-Malu iron ore mine is located in a highly tectonic area and the rock mass completely crushed, the Hoek-Brown failure criterion was found suitable to estimate geo-mechanical parameters. After that, the value of cohesion （c） and friction angle （tp） were calculated for different geotechnical zones and relative graphs and equations were derived as a function of slope height. The stability analyses using numerical and limit equilibrium methods showed that some instability problems might occur by increasing the slope height. Therefore, stable slopes for each geotechnical zone and prepared sections were calculated and presented as a function of slope height.展开更多
This paper describes a precise method combining numerical analysis and limit equilibrium theory to determine potential slip surfaces in soil slopes. In this method, the direction of the critical slip surface at any po...This paper describes a precise method combining numerical analysis and limit equilibrium theory to determine potential slip surfaces in soil slopes. In this method, the direction of the critical slip surface at any point in a slope is determined using the Coulomb’s strength principle and the extremum principle based on the ratio of the shear strength to the shear stress at that point. The ratio, which is considered as an analysis index, can be computed once the stress field of the soil slope is obtained. The critical slip direction at any point in the slope must be the tangential direction of a potential slip surface passing through the point. Therefore, starting from a point on the top of the slope surface or on the horizontal segment outside the slope toe, the increment with a small distance into the slope is used to choose another point and the corresponding slip direction at the point is computed. Connecting all the points used in the computation forms a potential slip surface exiting at the starting point. Then the factor of safety for any potential slip surface can be computed using limit equilibrium method like Spencer method. After factors of safety for all the potential slip surfaces are obtained, the minimum one is the factor of safety for the slope and the corresponding potential slip surface is the critical slip surface of the slope. The proposed method does not need to pre-assume the shape of potential slip surfaces. Thus it is suitable for any shape of slip surfaces. Moreover the method is very simple to be applied. Examples are presented in this paper to illustrate the feasibility of the proposed method programmed in ANSYS software by macro commands.展开更多
In order to expand the application of strength reduction methods with the ubiquitous-joint criterion, the corresponding program is compiled using FLAC3D software. A procedure for strength reduction in the ubiquitous-j...In order to expand the application of strength reduction methods with the ubiquitous-joint criterion, the corresponding program is compiled using FLAC3D software. A procedure for strength reduction in the ubiquitous-joint criterion is proposed to study the safety factor of slopes as well as the relationships of the bedding plane inclination angle β and the safety factor F. The results show that： 1） for the bedding rock slope, the various failure modes cause different variations of the safety factor F; 2） a bed- ding rock slope can be divided into two types by the angle between the bedding plane inclination and slope surface inclination a; when a_〈45~, the bedding slope can be strictly defined as the subsequent bedding rock slope; when 45°〈α〈90°, the bedding slope is defined as an oblique bedding slope; 3） for bedding rock slopes, the safety factor increases with an increase in a; for inverse bed- ding slopes, when the bedding plane inclination angle fl is small, the safety factor F of the slope increases at first, then decreases with an increase in a; when β is large, the safety factor F increases with an increase in α.展开更多
Based on the limit analysis upper bound method, a new mechanism of soil slope failure has been proposed which was consisted of plastic shear zone and rigid block zone. The different zones interface were regarded as di...Based on the limit analysis upper bound method, a new mechanism of soil slope failure has been proposed which was consisted of plastic shear zone and rigid block zone. The different zones interface were regarded as discontinuity lines. Two sliding blocks of the slope were also incorporated horizontal seismic force and vertical gravity force. The velocities and forces were analyzed in two blocks, and the expression of velocity discontinuities was derived according to the principle of incompressibility. The external force done work for the blocks and the internal energy dissipated of the plastic shear zone and the velocity discontinuous were solved.The stability ratios were derived for the height of two-level slope with different rates to involve seismic and no seismic. The present stability ratios were compared to the previous study, which showed the superiority of the mechanism and the rationality of the analysis. The critical height of the slope can provide a theoretical basis for slope support and design.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.51890912,51979025 and 52011530189).
文摘This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is employed to analyze the stability of 3D anisotropic soil slopes.The accuracy of the proposed method is first verified against the data in the literature.We then simulate the 3D soil slope with a straight slope surface and the convex and concave slope surfaces with a 90turning corner to study the 3D effect on slope stability and the failure mechanism under anisotropy conditions.Based on our numerical results,the end effect significantly impacts the failure mechanism and safety factor.Anisotropy degree notably affects the safety factor,with higher degrees leading to deeper landslides.For concave slopes,they can be approximated by straight slopes with suitable boundary conditions to assess their stability.Furthermore,a case study of the Saint-Alban test embankment A in Quebec,Canada,is provided to demonstrate the applicability of the proposed FE model.
基金the financial support provided by MHRD,Govt.of IndiaCoal India Limited for providing financial assistance for the research(Project No.CIL/R&D/01/73/2021)the partial financial support provided by the Ministry of Education,Government of India,under SPARC project(Project No.P1207)。
文摘Sudden and unforeseen seismic failures of coal mine overburden(OB)dump slopes interrupt mining operations,cause loss of lives and delay the production of coal.Consideration of the spatial heterogeneity of OB dump materials is imperative for an adequate evaluation of the seismic stability of OB dump slopes.In this study,pseudo-static seismic stability analyses are carried out for an OB dump slope by considering the material parameters obtained from an insitu field investigation.Spatial heterogeneity is simulated through use of the random finite element method(RFEM)and the random limit equilibrium method(RLEM)and a comparative study is presented.Combinations of horizontal and vertical spatial correlation lengths were considered for simulating isotropic and anisotropic random fields within the OB dump slope.Seismic performances of the slope have been reported through the probability of failure and reliability index.It was observed that the RLEM approach overestimates failure probability(P_(f))by considering seismic stability with spatial heterogeneity.The P_(f)was observed to increase with an increase in the coefficient of variation of friction angle of the dump materials.Further,it was inferred that the RLEM approach may not be adequately applicable for assessing the seismic stability of an OB dump slope for a horizontal seismic coefficient that is more than or equal to 0.1.
基金the project PID2022-139202OB-I00Neural Networks and Optimization Techniques for the Design and Safe Maintenance of Transportation Infrastructures:Volcanic Rock Geotechnics and Slope Stability(IA-Pyroslope),funded by the Spanish State Research Agency of the Ministry of Science,Innovation and Universities of Spain and the European Regional Development Fund,MCIN/AEI/10.13039/501100011033/FEDER,EU。
文摘Volcanic terrains exhibit a complex structure of pyroclastic deposits interspersed with sedimentary processes,resulting in irregular lithological sequences that lack lateral continuity and distinct stratigraphic patterns.This complexity poses significant challenges for slope stability analysis,requiring the development of specialized techniques to address these issues.This research presents a numerical methodology that incorporates spatial variability,nonlinear material characterization,and probabilistic analysis using a Monte Carlo framework to address this issue.The heterogeneous structure is represented by randomly assigning different lithotypes across the slope,while maintaining predefined global proportions.This contrasts with the more common approach of applying probabilistic variability to mechanical parameters within a homogeneous slope model.The material behavior is defined using complex nonlinear failure criteria,such as the Hoek-Brown model and a parabolic model with collapse,both implemented through linearization techniques.The Discontinuity Layout Optimization(DLO)method,a novel numerical approach based on limit analysis,is employed to efficiently incorporate these advances and compute the factor of safety of the slope.Within this framework,the Monte Carlo procedure is used to assess slope stability by conducting a large number of simulations,each with a different lithotype distribution.Based on the results,a hybrid method is proposed that combines probabilistic modeling with deterministic design principles for the slope stability assessment.As a case study,the methodology is applied to a 20-m-high vertical slope composed of three lithotypes(altered scoria,welded scoria,and basalt)randomly distributed in proportions of 15%,60%,and 25%,respectively.The results show convergence of mean values after approximately 400 simulations and highlight the significant influence of spatial heterogeneity,with variations of the factor of safety between 5 and 12 in 85%of cases.They also reveal non-circular and mid-slope failure wedges not captured by traditional stability methods.Finally,an equivalent normal probability distribution is proposed as a reliable approximation of the factor of safety for use in risk analysis and engineering decision-making.
基金supported by the National Natural Science Foundation of China(No.42277175)Guangxi Emergency Management Department 2024 Innovation and Technology Research Project,China(No.2024GXYJ006)+2 种基金Hunan Provincial Department of Natural Resources Geological Exploration Project,China(No.2023ZRBSHZ056)The First National Natural Disaster Comprehensive Risk Survey in Hunan Province,China(No.2022-70)Guizhou Provincial Major Scientific and Technological Program,China(No.2023-425).
文摘A method combining finite difference method(FDM)and k-means clustering algorithm which can determine the threshold of rock bridge generation is proposed.Jointed slope models with different joint coalescence coefficients(k)are constructed based on FDM.The rock bridge area was divided through k-means algorithm and the optimal number of clusters was determined by sum of squared errors(SSE)and elbow method.The influence of maximum principal stress and stress change rate as clustering indexes on the clustering results of rock bridges was compared by using Euclidean distance.The results show that using stress change rate as clustering index is more effective.When the joint coalescence coefficient is less than 0.6,there is no significant stress concentration in the middle area of adjacent joints,that is,no generation of rock bridge.In addition,the range of rock bridge is affected by the coalescence coefficient(k),the relative position of joints and the parameters of weak interlayer.
基金financially supported by Youth Innovation Promotion Association,CAS(Grant No.2022333)Knowledge Innovation Program of Wuhan e Basic Research(Grant No.2022010801010161)Natural Science Foundation of Hubei Province,China(Grant No.2023AFD219).
文摘Block-flexure toppling constitutes the predominant form of toppling failure in rock slopes.Although it has been extensively studied,the current theoretical models are often oversimplified by treating rock layers as rigid bodies that diverge from actual conditions.The proposed Equivalent Deformation Compatibility Method(EDCM)offers a fresh approach to assess the stability of rock slopes prone to block-flexure toppling.EDCM posits that blocky rock layers,with their inability to withstand significant bending and role in merely transferring forces,can be modeled as intact layers with a reduced modulus.The method simplifies the complex issue of analyzing discrete and continuous rock layers to the study of layered soft and hard rock,establishing deformation compatibility equations subsequently.Validation of the EDCM was achieved through numerical models,physical model testing,and application to an actual slope.The factor of safety(FS)for slopes corresponds with the results from both models and the actual slope,demonstrating the method's applicability for evaluating susceptibility to block-flexure toppling.When applying the EDCM,it is advised to set the elastic modulus reduction coefficient for blocky layers at a value below 0.1.
基金Project(52108361)supported by the National Natural Science Foundation of ChinaProjects(BK20231217,BK20220265)supported by the Basic Research Program of Jiangsu Province,China+5 种基金Project(sklhse-KF-2025-D-02)supported by the Open Research Fund Program of the State Key Laboratory of Hydroscience and Engineering,ChinaProject(2023ZB15)supported by the Independent Research Project of the State Key Laboratory of Subtropical Building and Urban Science,ChinaProject(SKLGME023001)supported by the Key Laboratory of Geomechanics and Geotechnical Engineering Safety,the Chinese Academy of SciencesProject(2025A04J3992)supported by the Basic and Applied Basic Research Project of the Guangzhou Science and Technology Bureau,ChinaProject(SKLGP2022Z015)supported by the State Key Laboratory of Geohazard Prevention and Geoenvironment Protection Independent Research Project,ChinaProjects(2023YFS0436,2024NSFSC1715)supported by the Science and Technology Department of Sichuan Province,China。
文摘Steep bedding slopes are widely distributed in Southwestern China’s mountainous regions and have complex seismic responses and instability risks,causing casualties and property losses.Considering the high-seismic-intensity environment,the dynamic failure evolution and instability mechanism of high-steep bedding slopes are simulated via the discrete element method and shaking table test.The dynamic response characteristics and cumulative failure effects of slopes subjected to continuous ground motion are investigated.The results show that the dynamic response characteristics of slopes under continuous earthquakes are influenced by geological and topographic conditions.Elevation has a distinct impact on both the slope interior and surface,with amplification effects more pronounced on the surface.The weak interlayers have different influences on the dynamic amplification effect of slopes.Weak interlayers have dynamic magnification effects on the slope surface at relative elevations of 0-0.33 and 0.82-1.0 but have weakening effects between 0.33 and 0.82.Moreover,the weak interlayers also have controlling effects on the dynamic instability mode of slopes.The characteristics of intergranular contact failure,fracture propagation,and displacement distribution are analyzed to reveal the dynamic failure evolution and instability mechanism through the discrete-element model.The dynamic instability process of slopes includes three stages:fracture initiation(0-0.2g),fracture expansion(0.2g-0.3g),and sliding instability(0.3g-0.6g).This work can provide a valuable reference for the seismic stability and reinforcement of complex slopes.
基金supported by the National Natural Science Foundation of China(Grant No.42377191,42072300)“The 14th Five Year Plan”Hubei Provincial advantaged characteristic disciplines(groups)project of Wuhan University of Science and Technology(Grant No.2023A0303).
文摘The stability of the tunnel portal slope is crucial for ensuring safe tunnel construction.Thus,a sound stability evaluation is of significance.Given the unique geological characteristics of tunnel portal slopes,it is necessary to establish a specific evaluation indicator system that differs from those used for ordinary slopes.Based on the unascertained measure method,uncertainties in the indicator are addressed by introducing the left and right half cloud asymmetric cloud model to optimize the linear membership function.The subjectivity of confidence criterion level identification is also improved by using the Euclidean distance method.Thus,a stability evaluation model for the tunnel portal slope is established based on the improved unascertained measure method.Finally,using the collected tunnel portal slope data,the results of four evaluation methods are compared with the safety factor levels.The evaluation methods include the traditional unascertained measure method,the method improved by using the left and right half cloud asymmetric cloud model,the method improved by using the Euclidean distance method,and the method improved by using both the left and right half cloud asymmetric cloud model and the Euclidean distance method.The results show that the accuracy rates of these four methods are 50%,55%,85%,and 90%,respectively.Among them,the joint improvement method has the slightest deviation,with only one level,while the other three methods had deviations of two levels.This result verifies the stability and effectiveness of the joint improvement method,providing a reference for tunnel portal slope stability evaluation.
基金Project(52278380)supported by the National Natural Science Foundation of ChinaProject(2023JJ30670)supported by the National Science Foundation of and Technology Major Project of Hunan Province,China。
文摘This study proposes an alternative calculation mode for stresses on the slip surface(SS).The calculation of the normal stress(NS)on the SS involves examining its composition and expanding its unknown using the Taylor series.This expansion enables the reasonable construction of a function describing the NS on the SS.Additionally,by directly incorporating the nonlinear Generalized Hoke-Brown(GHB)strength criterion and utilizing the slope factor of safety(FOS)definition,a function of the shear stress on the SS is derived.This function considers the mutual feedback mechanism between the NS and strength parameters of the SS.The stress constraints conditions are then introduced at both ends of the SS based on the spatial stress relation of one point.Determining the slope FOS and stress solution for the SS involves considering the mechanical equilibrium conditions and the stress constraint conditions satisfied by the sliding body.The proposed approach successfully simulates the tension-shear stress zone near the slope top and provides an intuitive description of the concentration effect of compression-shear stress of the SS near the slope toe.Furthermore,compared to other methods,the present method demonstrates superior processing capabilities for the embedded nonlinear GHB strength criterion.
基金Supported by the Program of Yunnan Provincial Institute of Communications Planning,Design and Research (2011(D)11-b)
文摘The vibration characteristics and dynamic responses of rock and soil under seismic load can be estimated with dynamic finite element method （DFEM）. Combining with the DFEM, the vector sum analysis method （VSAM） is employed in seismic stability analysis of a slope in this paper. Different from other conventional methods, the VSAM is proposed based on the vector characteristic of force and current stress state of the slope. The dynamic stress state of the slope at any moment under seismic load can he obtained by the DFEM, thus the factor of safety of the slope at any moment during earthquake can be easily obtained with the VSAM in consideration of the DFEM. Then, the global stability of the slope can be estimated on the basis of time-history curve of factor of safety and reliability theory. The VSAM is applied to a homogeneous slope under seismic load. The factor of safety of the slope is 1.30 under gravity only and the dynamic factor of safety under seismic load is 1.21. The calculating results show that the dynamic characteristics and stability state of the slope with input ground motion can be actually analyzed. It is believed that the VSAM is a feasible and practical approach to estimate the dynamic stability of slopes under seismic load.
基金supported by The National Key Research and Development Program Plane(No.2017YFC0601505)National Natural Science Foundation(No.41672325)Science&Technology Department of Sichuan Province Technology Project(No.2017GZ0393)
文摘The attenuation factor or quality factor(Q-factor or Q) has been used to measure the energy attenuation of seismic waves propagating in underground media. Many methods are used to estimate the Q-factor. We propose a method to calculate the Q-factor based on the prestack Q-factor inversion and the generalized S-transform. The proposed method specifies a standard primary wavelet and calculates the cumulative Q-factors; then, it finds the interlaminar Q-factors using the relation between Q and offset(QVO) and the Dix formula. The proposed method is alternative to methods that calculate interlaminar Q-factors after horizon picking. Because the frequency spectrum of each horizon can be extracted continuously on a 2D time–frequency spectrum, the method is called the continuous spectral ratio slope(CSRS) method. Compared with the other Q-inversion methods, the method offers nearly effortless computations and stability, and has mathematical and physical significance. We use numerical modeling to verify the feasibility of the method and apply it to real data from an oilfield in Ahdeb, Iraq. The results suggest that the resolution and spatial stability of the Q-profile are optimal and contain abundant interlaminar information that is extremely helpful in making lithology and fluid predictions.
基金Project(11102218) supported by the National Natural Science Foundation of China
文摘The core of strength reduction method(SRM) involves finding a critical strength curve that happens to make the slope globally fail and a definition of factor of safety(FOS). A new double reduction method, including a detailed calculation procedure and a definition of FOS for slope stability was developed based on the understanding of SRM. When constructing the new definition of FOS, efforts were made to make sure that it has concise physical meanings and fully reflects the shear strength of the slope. Two examples, slopes A and B with the slope angles of 63° and 34° respectively, were given to verify the method presented. It is found that, for these two slopes, the FOSs from original strength reduction method are respectively 1.5% and 38% higher than those from double reduction method. It is also found that the double reduction method predicts a deeper potential slide line and a larger slide mass. These results show that on one hand, the double reduction method is comparative to the traditional methods and is reasonable, and on the other hand, the original strength reduction method may overestimate the safety of a slope. The method presented is advised to be considered as an additional option in the practical slope stability evaluations although more useful experience is required.
基金support of the National Basic Research Program of China (No. 2011CB710606)the Geological Survey Program of the China Geological Survey (No. 1212010914036)the National Natural Science Foundation of China (No. 41102195)
文摘A new method （kinetic vector method, KVM） is presented for analyzing the dynamic stability of wedge in rock slope. The dynamic analysis is carried out based on three dimensional distinct element code （3DEC）, and the kinetic inertial force of the wedge under seismic loading can be obtained via calculating the net vectorial nodal force of the finite difference grid. Then, the factor of safety （FOS） of the wedge can be calculated based on limit equilibrium method （LEM） at each dynamic analysis step, therefore time series of the FOS for whole earthquake process can be obtained. For the purpose of evaluating the entire dynamic stability of the wedge, dynamic factor of safety （DFOS） is proposed and defined as a numerical value corresponding with a given rate of probability guarantee based on reliability theory. Consequently, the KVM inherits the merits of the LEM and also has fully nonlinear dynamic analysis capabilities, and the feasibility and correctness of the KVM are tested by an example given by Hoek and Bray （1981）. Finally, a rock slope case in Wenchuan Earthquake regions of China is presented to verify the engineering practicability of the KVM, and the results matched the actual situation well.
文摘Overhanging rock slopes(steeper than 90°) are typically avoided in rock engineering design, particularly where the scale of the slope exceeds the scale of fracturing present in the rock mass. This paper highlights an integrated approach of designing overhanging rock slopes where the relative dimensions of the slope exceed the scale of fracturing and the rock mass failure needs to be considered rather than kinematic release of individual blocks. The key to the method is a simplified limit equilibrium(LE) tool that was used for the support design and analysis of a multi-faceted overhanging rock slope. The overhanging slopes required complex geometries with constantly changing orientations. The overhanging rock varied in height from 30 m to 66 m. Geomechanical modelling combined with discrete fracture network(DFN)representation of the rock mass was used to validate the rock mass strength assumptions and the failure mechanism assumed in the LE model. The advantage of the simplified LE method is that buttress and support design iterations(along with sensitivity analysis of design parameters) can be completed for various cross-sections along the proposed overhanging rock sections in an efficient manner, compared to the more time-intensive, sophisticated methods that were used for the initial validation. The method described presents the development of this design tool and assumptions made for a specific overhanging rock slope design. Other locations will have different geological conditions that can control the potential behaviour of rock slopes, however, the approach presented can be applied as a general guiding design principle for overhanging rock cut slope.
基金NRDMS Division,Department of Science and Technology,Government of India for providing financial assistance for field investigations.
文摘Large-scale slope destabilization could be aggravated due to swift urbanization and ever-rising demands of geoengineering projects such as dams,tunnels,bridges and widening roads.National Highway-58 connects Delhi to Badrinath in India,which passes through complex geomorphological and geological terrain and often encounters cut slopes susceptible to slope failures.In the present investigation,a detailed geotechnical appraisal is conducted along the road cut slopes from Rishikesh to Devprayag in the Himalayas.Twenty vulnerable road cut slopes were demarcated for detailed slope stability analysis using Phase2D finite element modeling simulator.Nonlinear generalized Hoek-Brown(GHB)criterion was adopted for stability analyses.Out of 20 slopes,five slopes(S6,S7,S18,S19 and S20)are unstable with factor of safety(FoS)less than or equal to 1,and thus needs immediate attention.The FoS values of four slopes(S2,S9,S13 and S17)lie between 1 and 1.3,i.e.marginally stable,and slopes S1,S3,S4,S5,S8,S10,Sll,S12,S14,S15 and S16 are stable.Mohr-Coulomb(MC)criterion was also adopted to compare the slope stability analysis with GHB criterion.The FoS calculated from GHB criterion is close to that using MC criterion for lower values of FoS whereas for higher values,the difference is marked.For the jointed rock in the Himalayan region,the nonlinear GHB criterion gives better results as compared to MC criterion and matches with the prevailing field conditions.Accordingly,some suggestions are proposed to strengthen the stability of cut slopes.
基金National Natural Science Foundation of China under grant No.59839330China Postdoctoral Science Foundation
文摘The elliptic mild slope equation is used to simulate linear wave propagation over variable sea bed topography with mild slopes. The governing equation is discretized by the finite difference method. Based on the BI-CGSTAB technique, an attractive variant bf BI-Conjugate Gradients (BI-CG) method, the obtained linear algebraic system of equations is solved. Numerical experiments show that the BI-CGSTAB method is efficient for solving the elliptic mild slope equation. The results obtained by the BI-CGSTAB-Based method are much the same as those obtained by other authors with different solution methods, but the convergence rate is much faster than that of other methods.
基金evolution mechanism and prevention countermeasures of the Outang landslide in the Three Gorges Reservoir Area(No.20C0023)research projectthe geological safety risk investigation,evaluation and control of key resettlement towns in the Three Gorges Reservoir Area(No.HBHDZFCG2021025)+2 种基金the National Natural Science Foundation of China(No.42077268)the Chongqing Geological Disaster Prevention and Control Center of China(No.20C0023)the open fund of state key laboratory of geohazard prevention and geoenvironment protection(No.SKLGP2020K015)。
文摘Although slice methods are simple and effective slope stability analysis approaches,they are statically indeterminate.Several modifications of the slice method,such as the Spencer,MorgensternPrice,and Chen-Morgenstern methods,are statically determinate and solvable as they assume the inter-slice force inclination angle;however,there is a small gap between the assumptions and actual landslide stability analysis.Through reasonable theoretical analysis,the Su slice method provides a reliable approach for determining the inter-slice force inclination angle that can be used in slice analysis to accurately analyse,calculate,and evaluate the stability of landslides.However,the Su slice method requires further research and analysis,especially in terms of the parameter values sinλbiandρ.In this study,we investigated more accurate methods for calculating the parameters sinλbiandρ.In addition,an adjustment coefficient(μ)was introduced to improve the solution method for the inter-slice force inclination angle.The inter-slice force inclination and safety factors of three landslides with arc-shaped slip surfaces and one landslide with a polyline-shaped slip surface were analysed and compared using the different slice methods.The improved inter-slice force inclination not only satisfies the calculation of static force equilibrium condition but also satisfies the calculation of both the force and moment equilibrium conditions.The improved method for calculating inter-slice force inclination presented the best correlation.The safety factors calculated using the improved Su slice method were close to those obtained using numerical simulations and the Morgenstern-Price method.Despite negligible differences among the safety factors calculated using the Su slice,improved Su slice,and M-P methods,the accuracy of the improved Su slice method was better than the M-P method in terms of inter-slice force inclination angles which can be useful to improve protection engineering design.
文摘One of the critical aspects in mine design is slope stability analysis and the determination of stable slopes. In the Chador- Malu iron ore mine, one of the most important iron ore mines in central Iran, it was considered vital to perform a comprehensive slope stability analysis. At first, we divided the existing rock hosting pit into six zones and a geotechnical map was prepared. Then, the value of MRMR （Mining Rock Mass Rating） was determined for each zone. Owing to the fact that the Chador-Malu iron ore mine is located in a highly tectonic area and the rock mass completely crushed, the Hoek-Brown failure criterion was found suitable to estimate geo-mechanical parameters. After that, the value of cohesion （c） and friction angle （tp） were calculated for different geotechnical zones and relative graphs and equations were derived as a function of slope height. The stability analyses using numerical and limit equilibrium methods showed that some instability problems might occur by increasing the slope height. Therefore, stable slopes for each geotechnical zone and prepared sections were calculated and presented as a function of slope height.
文摘This paper describes a precise method combining numerical analysis and limit equilibrium theory to determine potential slip surfaces in soil slopes. In this method, the direction of the critical slip surface at any point in a slope is determined using the Coulomb’s strength principle and the extremum principle based on the ratio of the shear strength to the shear stress at that point. The ratio, which is considered as an analysis index, can be computed once the stress field of the soil slope is obtained. The critical slip direction at any point in the slope must be the tangential direction of a potential slip surface passing through the point. Therefore, starting from a point on the top of the slope surface or on the horizontal segment outside the slope toe, the increment with a small distance into the slope is used to choose another point and the corresponding slip direction at the point is computed. Connecting all the points used in the computation forms a potential slip surface exiting at the starting point. Then the factor of safety for any potential slip surface can be computed using limit equilibrium method like Spencer method. After factors of safety for all the potential slip surfaces are obtained, the minimum one is the factor of safety for the slope and the corresponding potential slip surface is the critical slip surface of the slope. The proposed method does not need to pre-assume the shape of potential slip surfaces. Thus it is suitable for any shape of slip surfaces. Moreover the method is very simple to be applied. Examples are presented in this paper to illustrate the feasibility of the proposed method programmed in ANSYS software by macro commands.
基金Project 5049027 supported by the National Natural Science Foundation of China
文摘In order to expand the application of strength reduction methods with the ubiquitous-joint criterion, the corresponding program is compiled using FLAC3D software. A procedure for strength reduction in the ubiquitous-joint criterion is proposed to study the safety factor of slopes as well as the relationships of the bedding plane inclination angle β and the safety factor F. The results show that： 1） for the bedding rock slope, the various failure modes cause different variations of the safety factor F; 2） a bed- ding rock slope can be divided into two types by the angle between the bedding plane inclination and slope surface inclination a; when a_〈45~, the bedding slope can be strictly defined as the subsequent bedding rock slope; when 45°〈α〈90°, the bedding slope is defined as an oblique bedding slope; 3） for bedding rock slopes, the safety factor increases with an increase in a; for inverse bed- ding slopes, when the bedding plane inclination angle fl is small, the safety factor F of the slope increases at first, then decreases with an increase in a; when β is large, the safety factor F increases with an increase in α.
基金the National Natural Science Foundation of China(No.51478444,51574115 and 51774121)
文摘Based on the limit analysis upper bound method, a new mechanism of soil slope failure has been proposed which was consisted of plastic shear zone and rigid block zone. The different zones interface were regarded as discontinuity lines. Two sliding blocks of the slope were also incorporated horizontal seismic force and vertical gravity force. The velocities and forces were analyzed in two blocks, and the expression of velocity discontinuities was derived according to the principle of incompressibility. The external force done work for the blocks and the internal energy dissipated of the plastic shear zone and the velocity discontinuous were solved.The stability ratios were derived for the height of two-level slope with different rates to involve seismic and no seismic. The present stability ratios were compared to the previous study, which showed the superiority of the mechanism and the rationality of the analysis. The critical height of the slope can provide a theoretical basis for slope support and design.
