An approach for determination of lateral limit angle in kinematic planar sliding analysis for rock slopes

Planar sliding is one of the frequently observed types of failure in rock slopes.Kinematic analysis is a classic and widely used method to examine the potential failure modes in rock masses.The accuracy of planar sliding kinematic analysis is significantly influenced by the value assigned to the lateral limit angleγlim.However,the assignment ofγlim is currently used generally based on an empirical criterion.This study aims to propose an approach for determining the value ofγlim in deterministic and probabilistic kinematic planar sliding analysis.A new perspective is presented to reveal thatγlim essentially influences the probability of forming a potential planar sliding block.The procedure to calculate this probability is introduced using the block theory method.It is found that the probability is correlated with the number of discontinuity sets presented in rock masses.Thus,different values ofγlim for rock masses with different sets of discontinuities are recommended in both probabilistic and deterministic planar sliding kinematic analyses;whereas a fixed value ofγlim is commonly assigned to different types of rock masses in traditional method.Finally,an engineering case was used to compare the proposed and traditional kinematic analysis methods.The error rates of the traditional method vary from 45%to 119%,while that of the proposed method ranges between 1%and 17%.Therefore,it is likely that the proposed method is superior to the traditional one.

Development of Finite Element Limiting Analysis Method and Its Applications in Geotechnical Engineering

The Finite Element Limiting Analysis Method(LELAM) has the advantage of combining a numerical analysis method with traditional limiting equilibrium methods.It is particularly applicable to the analysis and design of geotechnical engineering.In the early 20th century,FELAM has been developed vigorously in domestic geotechnical engineering over international common finite element procedures.It has made great achievements in basic theory research and computational precision,thus broadening the application fields in practical projects.In order to gradually make innovations in geotechnical design methods,some of our research results are presented,mainly including geotechnical safety factor definitions,the principles for use of the method concerned,the overall failure criterion,the deduction and selection of the yield criterion,and the measurement to improve the computational precision,etc..The application field has been broadened from two-dimensional to three-dimensional,from soil slope to jointed rock slope and foundation,from stable seepage to non-stable seepage,from slope and foundation to tunnel.This method has also been used in search of many hidden sliding surfaces of complex landslides,conducting the structural support design considering the interaction between the soil and the structure,and computing simulation foundation bearing plates load tests,etc..

Application of dynamic analysis of strength reduction in the slope engineering under earthquake

At present,the methods of analyzing the stability of slope under earthquake are not accurate and reasonable because of some limitations. Based on the real dynamic tensile-shear failure mechanism of slope,the paper proposes dynamic analysis of strength reduction FEM (finite element method) and takes the reduction of shear strength parameters and tensile strength parameters into consideration. And it comprehensively takes the transfixion of the failure surface,the non-convergence of calculation and mutation of displacement as the criterion of dynamic instability and failure of the slope. The strength reduction factor under limit state is regarded as the dynamic safety factor of the slope under earthquake effect and its advantages are introduced. Finally,the method is applied in the seismic design of anchors supporting and anti-slide pile supporting of the slope. Calculation examples show that the application of dynamic analysis of strength reduction is feasible in the seismic design of slope engineering,which can consider dynamic interaction of supporting structure and rock-soil mass. Owing to its preciseness and great advantages,it is a new method in the seismic design of slope supporting.

Engineering Geological and Geotechnical Studies of Taprang Landslide, West-central Nepal: An Approach for Slope Stability Analysis

Detailed investigation of Taprang landslide was carried out in order tounderstand the surface, subsurface lithological information and physicalproperties of soil by using multi-disciplinary methods such as engineeringgeological, geophysical and geotechnical studies for the determinationof factor of safety for slope stability analysis. Geological study wascarried out by detail mapping of surface geology, soil condition, propertiesof bedrock and its discontinuities. The geophysical survey (ElectricalResistivity Tomography-ERT) were carried out to know the electricalresistivity of soil for identifying the groundwater table and slip surface ofthe landslide. Geotechnical analysis such as grain size analysis, liquid limitand direct shear test were carried out in order to evaluate soil classification,moisture content, cohesion and the angle of internal friction of soil forknowing the strength the soil. These soil parameters indicate the soil is verylow strength. The combination of these results were used for calculatingthe factor of safety (FoS) by Limit Equilibrium Method (LEM) proposedby Bishop and Janbu methods. The result of factor of safety in the Tapranglandslide demonstrates that the slope become stable in drained (dry)condition, remain ultimate stage in undrained (wet) condition and finallyfailure occurs if applied the seismic load in both drained and undrainedconditions.

Combined influence of nonlinearity and dilation on slope stability evaluated by upper-bound limit analysis

The combined influence of nonlinearity and dilation on slope stability was evaluated using the upper-bound limit analysis theorem.The mechanism of slope collapse was analyzed by dividing it into arbitrary discrete soil blocks with the nonlinear Mohr–Coulomb failure criterion and nonassociated flow rule.The multipoint tangent(multi-tangent) technique was used to analyze the slope stability by linearizing the nonlinear failure criterion.A general expression for the slope safety factor was derived based on the virtual work principle and the strength reduction technique,and the global slope safety factor can be obtained by the optimization method of nonlinear sequential quadratic programming.The results show better agreement with previous research result when the nonlinear failure criterion reduces to a linear failure criterion or the non-associated flow rule reduces to an associated flow rule,which demonstrates the rationality of the presented method.Slope safety factors calculated by the multi-tangent inclined-slices technique were smaller than those obtained by the traditional single-tangent inclined-slices technique.The results show that the multi-tangent inclined-slices technique is a safe and effective method of slope stability limit analysis.The combined effect of nonlinearity and dilation on slope stability was analyzed,and the parameter analysis indicates that nonlinearity and dilation have significant influence on the result of slope stability analysis.

Three-dimensional stability analysis of anisotropic and non-homogeneous slopes using limit analysis

A method of three-dimensional loaded slope stability for anisotropic and nonhomogeneous slopes was presented based on the upper-bound theorem of the limit analysis approach. The approach can be considered as a modification and extension of the solutions. The influences of friction angle, anisotropy factor, nonhomogeneous factor, slope angle, ratio of width to depth, and load on the slope crest were investigated. The results show that solutions are suitable to deal with the purely cohesive soils and frictional/cohesive soils, isotropic and anisotropic, homogeneous and nonhomogeneous, loaded and unloaded cases.

Force analysis of pile foundation in rock slope based on upper-bound theorem of limit

Based on the characteristic that the potential sliding surfaces of rock slope are commonly in the shape of either line or fold line,analysis thought of conventional pile foundation in the flat ground under complex load condition was applied and the upper-bound theorem of limit analysis was used to compute thrust of rock layers with all possible distribution shapes. The interaction of slope and pile was considered design load in terms of slope thrust,and the finite difference method was derived to calculate inner-force and displacement of bridge pile foundation in rock slope under complex load condition. The result of example shows that the distribution model of slope thrust has certain impact on displacement and inner-force of bridge pile foundation. The maximum displacement growth rate reaches 54% and the maximum moment and shear growth rates reach only 15% and 20%,respectively,but the trends of inner-force and displacement of bridge pile foundation are basically the same as those of the conventional pile foundation in the flat ground. When the piles bear the same level lateral thrust,the distribution shapes of slope thrust have different influence on inner-force of pile foundation,especially the rectangle distribution,and the triangle thrust has the smallest displacement and inner-force of pile foundation.

A Method Combining Numerical Analysis and Limit Equilibrium Theory to Determine Potential Slip Surfaces in Soil Slopes

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.

A pseudo-coupled numerical approach for stability analysis of frozen soil slopes based on finite element limit analysis method

To simplify the stability analysis of frozen soil slope, a pseudo-coupled numerical approach is developed. In this approach, the coupled heat transfer and water flow in frozen soils are simulated first, and based on the computed thermal-hydro field, the stability of frozen soil slope is evaluated. Although the shear strength for frozen soil is very complicated and is usually represented by a nonlinear MC failure criterion, a simple linear MC yield criterion is utilized. In this method, the internal friction angle is expressed as a function of volumetric ice content and the cohesion is fitted as a simple bilinear expression of Tand volumetric water content. To assess slope stability, the limit analysis is employed in conjunction with the recently developed a-section search algorithm. A frozen soil slope example is used to examine the proposed pseudo-coupled numerical approach, and numerical studies validate its effectiveness. Based on numerical results, it is seen that slope stability may be remarkably influenced by warming air （or grotmd surface） temperature. With increasing ground surface temperature, slope stability indicated by FOS may reduce to 1.0, implying that wanning air temperature could be a trigger of frozen soil slope failure.

New algorithm of mine slope reliability based on limiting state hyper-plane and its engineering application

Due to the influence of joint fissure, mining intensity, designed slope angle, underground water and rainfall, the failure process of mine slope project is extremely complicated. The current safety factor calculation method has certain limitations, and it would be difficult to obtain the reliability index when the performance function of reliability analysis is implicit or has high order terms. Therefore, with the help of the logistic equation of chaos theory, a new algorithm of mine slope reliability based on limiting state hyper-plane is proposed. It is shown that by using this new reliability algorithm the calculation of partial derivative of performance function is avoided, and it has the advantages of being simple and easy to program. The new algorithm is suitable for calculating the reliability index of complex performance function containing high order terms. Furthermore, the limiting state hyper-plane models of both simplified Bishop's and Janbu's method adaptive to slope project are obtained, and have achieved satisfactory effect in the study of mine slope stability in Dexing copper open pit.

3D limit analysis of rock slopes based on equivalent linear failure criterion with tension cut-off

Hoek-Brown (HB) failure criterion is widely used to predict the strength of intact or heavily jointed rock mass. For stability analysis of rock slopes governed by the HB failure criterion, the equivalent linearity to Mohr-Coulomb (MC) criterion is often adopted, leading to the well-known equivalent Mohr-Coulomb method (EMCM). Existing studies on EMCM analysis mainly consider the shear strength of rock material, while consideration of the tensile strength is rare. This contradicts the fact that the underlying tensile strength of rock mass has considerable impact on the rock slope stability in real world. In this regard, this paper proposes a limit analysis-based approach that can account for tension in the three-dimensional (3D) stability analysis of HB rock slope. This approach is established on the equivalent linearity of the HB criterion with consideration of tensile strength, known as the equivalent tension cut-off MC method (ETMCM), and using a horn-like 3D mechanism of limit analysis. The safety factor solutions given by the proposed approach are validated by previous studies and numerical results. Parametric studies are conducted to investigate the effect of rock tensile strength on slope stability. Results show that the consideration of tension leads to a more conservative safety factor and a sharper curvature of the failure surface, and these impacts tend to be more obvious with the increases in slope inclination and slope width. Finally, the stability of the HB rock slope under seepage conditions is studied using the proposed approach. The results indicate that the effect of tensile strength is highly remarkable in seepage circumstances.

Probabilistic approach for open pit bench slope stability analysis——A mine case study

The geotechnical slope design of an open pit wall starts at the bench scale configuration.At this scale,the rock slope stability is governed primarily by the geological discontinuities within the rock mass and as a result,structurally-controlled failures(e.g.planar,wedge or toppling)are most likely to occur.The probabilistic approach offers a major advantage over the traditional deterministic method in that it accounts for the different degrees of variability and uncertainty often encountered in rock properties.This paper presents a bench slope stability assessment for an open pit mine in Peru using a probabilistic-based approach by coupling a kinematic analysis based on stereographic projection techniques followed by a kinetic analysis by means of the limit equilibrium method.Finally,these two probabilities are combined to provide an overall measure of the probability of failure(PoF)of the bench slope system.The case study is characterized by significant scatter in the geometrical and mechanical properties of the joints.Extensive surface mapping was conducted at 36 different sites following the ISRM suggested procedures.Several direct shear tests were carried out.It is shown that by combining field and laboratory measurements and engineering judgment,the probability density functions(PDF)of the discontinuity parameters can be obtained.These are then used in a Monte Carlo simulation process to compute both kinematic and kinetic probabilities of failure.The overall probability of failure aims to provide the design engineer with a tool to critically evaluate the bench performance from a geotechnical risk perspective and to provide a basis for future bench design optimization.

Upper bound analysis of slope stability with nonlinear failure criterion based on strength reduction technique

Based on the upper bound limit analysis theorem and the shear strength reduction technique, the equation for expressing critical limit-equilibrium state was employed to define the safety factor of a given slope and its corresponding critical failure mechanism by means of the kinematical approach of limit analysis theory. The nonlinear shear strength parameters were treated as variable parameters and a kinematically admissible failure mechanism was considered for calculation schemes. The iterative optimization method was adopted to obtain the safety factors. Case study and comparative analysis show that solutions presented here agree with available predictions when nonlinear criterion reduces to linear criterion, and the validity of present method could be illuminated. From the numerical results, it can also be seen that nonlinear parameter rn, slope foot gradient ,β, height of slope H, slope top gradient a and soil bulk density γ have significant effects on the safety factor of the slope.

Application of cusp catastrophe theory to reliability analysis of slopes in open-pit mines

A method of slope reliability analysis was developed by imposing a state equation on the limit equilibrium theory, given the basis of a fixed safety factor technique. Among the many problems of reliability analysis, the most important problem is to find a performance function. We have created a new method of building a limit state equation for planar slip surfaces by applying the mathematical cusp catastrophe theory. This new technique overcomes the defects in the traditional rigid limit equilibrium theory and offers a new way for studying the reliability problem of planar slip surfaces. Consequently, we applied the technique to a case of an open-pit mine and compared our results with that of the traditional approach. From the results we conclude that both methods are essentially consistent, but the reliability index calculated by the traditional model is lower than that from the catastrophic model. The catastrophe model takes into consideration two possible situations of a slope being in the limit equilibrium condition, i.e., it may or may not slip. In the traditional method, however, a slope is definitely considered as slipping when it meets the condition of a limit equilibrium. We conclude that the catastrophe model has more actual and instructive importance compared to the traditional model.

Three-dimensional analysis of slopes reinforced with piles

Based on the upper bound of limit analysis, the plane-strain analysis of the slopes reinforced with a row of piles to the 3D case was extended. A 3D rotational failure mechanism was adopted to yield the upper bound of the factor of safety. Parametric studies were carried out to explore the end effects of the slope failures and the effects of the pile location and diameter on the safety of the reinforced slopes. The results demonstrate that the end effects nearly have no effects on the most suitable location of the installed piles but have significant influence on the safety of the slopes. For a slope constrained to a narrow width, the slope becomes more stable owing to the contribution of the end effects. When the slope is reinforced with a row of piles in small space between piles, the effects of group piles are significant for evaluating the safety of slopes. The presented method is more appropriate for assessing the stability of slopes reinforced with piles and can be also utilized in the design of plies stabilizing the unstable slopes.

Pseudo-static stability analysis of rock slopes reinforced by passive bolts using the generalized Hoek-Brown criterion

The stability analysis of passive bolt-reinforced rock slopes under seismic loads is investigated within the framework of the kinematic approach of limit analysis theory.A pseudo-static method is adopted to account for the inertial forces induced in the rock mass by seismic events.The strength properties of the rock material are described by a modified Hoek-Brown strength criterion,whereas the passive bolts are modeled as bar-like inclusions that exhibit only resistance to tensile-compressive forces.Taking advantage of the ability to compute closed-form expressions for the support functions associated with the modified Hoek-Brown strength criterion,a rotational failure mechanism is implemented to derive rigorous lower bound estimates for the amount of reinforcement strength to prevent slope failure.The approach is then applied to investigating the effects of relevant geometry,strength and loading parameters in light of a preliminary parametric study.The accuracy of the approach is assessed by comparison of the lower bound estimates with finite element limit analysis solutions,thus emphasizing the ability of the approach to properly predict the stability conditions and to capture the essential features of deformation localization pattern.Finally,the extension of the approach to account for slipping at the interface between reinforcements and surrounding rock mass is outlined.

Limit load and failure mechanisms of a vertical Hoek-Brown rock slope

The problem considered in this short note is the limit load determination of a vertical rock slope.The classical limit theorem is employed with the use of adaptive finite elements and nonlinear programming to determine upper and lower bound limit loads of a Hoek-Brown vertical rock slope.The objective function of the mathematical programming problem is such as to optimize a boundary load,which is known as the limit load,resembling the ultimate bearing capacity of a strip footing.While focusing on the vertical slope,parametric studies are carried out for several dimensionless ratios such as the dimensionless footing distance ratio,the dimensionless height ratio,and the dimensionless rock strength ratio.A comprehensive set of design charts is presented,and failure envelopes shown with the results explained in terms of three identified failure mechanisms,i.e.the face,the toe,and the Prandtl-type failures.These novel results can be used with great confidence in design practice,in particularly noting that the current industry-based design procedures for the presented problem are rarely found.

Multivariate adaptive regression splines analysis for 3D slope stability in anisotropic and heterogenous clay

Little research can be found in relation to the stability of anisotropic and heterogenous soils in three dimensions.In this paper,we propose a study on the three-dimensional(3D)undrained slopes in anisotropic and heterogenous clay using advanced upper and lower bounds finite element limit analysis(FELA).The obtained stability solutions are normalized,and presented by a stability number that is a function of three geometrical ratios and two material ratios,i.e.depth ratio,length ratio,slope angle,shear strength gradient ratio and anisotropic strength ratio.Numerical results are compared with experimental data in the literature,and charts are presented to cover a wide range of design parameters.Using the multivariate adaptive regression splines(MARS)analysis,the respective influence and sensitivity of each design parameter on the stability number and the failure mechanism are investigated.An empirical equation is also developed to effectively estimate the stability number.

Seismic stability analysis of slopes with pre-existing slip surfaces

In analyzing seismic stability of a slope with upper bound limit analysis method, the slip surface is often assumed as a log-spiral or plane slip surface. However, due to the presence of a weak layer and unfavorable geological structural surface or a bedrock interface with overlying soft strata, the preexisting slip surface of the slope may be irregular and composed of a series of planes rather than strictly logspiral or plane shape. A computational model is developed for analyzing the seismic stability of slopes with pre-existing slip surfaces. This model is based on the upper bound limit analysis method and can consider the effect of anchor bolts. The soil or rock is deemed to follow the Mohr-Coulomb yield criterion. The slope is divided into multiple block elements along the slip surface. According to the displacement compatibility and the associated flow rule, a kinematic velocity field of the slope can be obtained computationally. The proposed model allows not only calculation of the rate of external work owing to the combined effect of self-weight and seismic loading, but also that of the energy dissipation rate caused by the slip surface, interfaces of block elements and anchorage effect of the anchors. Considering a direct relationship between the rate of external work and the energy dissipation rate, the expressions of yield acceleration and permanent displacement of anchored slopes can be derived. Finally, the validity of this proposed model is illustrated by analysis on three typical slopes. The results showed that the proposed model is more easily formulated and does not need to solve complex equations or time consuming iterations compared with previous methods based on the conditions of force equilibrium.

The limit analysis in soil and rock:a mature discipline of geomechanics

The solution of a slope stability problem can be approached by its least upper-bound and maximum lower-bound with high accuracy. The limit equilibrium methods that employ vertical slices imply a lower bound of the factor of safety. It has been successfully extended to the area of active earth pressure analysis that accounts for different input of locations of earth pressure applications. Those methods that employ slices with inclined interfaces give an upper-bound approach to the stability analysis. It enjoys a sound mechanical background and is able to provide accurate solutions of soil plasticity. It has been successfully extended to the area of bearing capacity analysis in which various empirical coefficients are no longer necessary. The 3D upper- and lower-bound methods under this framework have been made possible and show great potential for solving various engineering problems.