Finite element analyses of slope stability problems using non-associated plasticity

In recent years, finite element analyses have increasingly been utilized for slope stability problems. In comparison to limit equilibrium methods, numerical analyses do not require any definition of the failure mechanism a priori and enable the determination of the safety level more accurately. The paper compares the performances of strength reduction finite element analysis(SRFEA) with finite element limit analysis(FELA), whereby the focus is related to non-associated plasticity. Displacement-based finite element analyses using a strength reduction technique suffer from numerical instabilities when using non-associated plasticity, especially when dealing with high friction angles but moderate dilatancy angles. The FELA on the other hand provides rigorous upper and lower bounds of the factor of safety(FoS) but is restricted to associated flow rules. Suggestions to overcome this problem, proposed by Davis(1968), lead to conservative FoSs; therefore, an enhanced procedure has been investigated. When using the modified approach, both the SRFEA and the FELA provide very similar results. Further studies highlight the advantages of using an adaptive mesh refinement to determine FoSs. Additionally, it is shown that the initial stress field does not affect the FoS when using a Mohr-Coulomb failure criterion.

A new double reduction method for slope stability analysis

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.

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.

Finite element analysis of slope stability by expanding the mobilized principal stress Mohr's circles-Development, encoding and validation

In recent years,finite element analysis is increasingly being proposed in slope stability problems as a competitive method to traditional limit equilibrium methods(LEMs)which are known for their inherent deficiencies.However,the application of finite element method(FEM)to slope stability as a strength reduction method(SRM)or as finite element limit analysis(FELA)is not always a success for the drawbacks that characterize both methods.To increase the performance of finite element analysis in this problem,a new approach is proposed in this paper.It consists in gradually expanding the mobilized stress Mohr’s circles until the soil failure occurs according to a prescribed non-convergence criterion.The present approach called stress deviator increasing method(SDIM)is considered rigorous for three main reasons.Firstly,it preserves the definition of the factor of safety(FOS)as the ratio of soil shear strength to the mobilized shear stress.Secondly,it maintains the progressive development of shear stress resulting from the increase in the principal stress deviator on the same plane,on which the shear strength takes place.Thirdly,by introducing the concept of equivalent stress loading,the resulting trial stresses are checked against the violation of the actual yield criterion formed with the real strength parameters rather than those reduced by a trial factor.The new numerical procedure was encoded in a Fortran computer code called S^(4)DINA and verified by several examples.Comparisons with other numerical methods such as the SRM,gravity increasing method(GIM)or even FELA by assessing both the FOS and contours of equivalent plastic strains showed promising results.

A Calculation Method of Double Strength Reduction for Layered Slope Based on the Reduction of Water Content Intensity

The calculation of the factor of safety(FOS)is an important means of slope evaluation.This paper proposed an improved double strength reductionmethod(DRM)to analyze the safety of layered slopes.The physical properties of different soil layers of the slopes are different,so the single coefficient strength reduction method(SRM)is not enough to reflect the actual critical state of the slopes.Considering that the water content of the soil in the natural state is the main factor for the strength of the soil,the attenuation law of shear strength of clayey soil changing with water content is fitted.This paper also establishes the functional relationship between different reduction coefficients.Then,a USDFLD subroutine is programmed using the secondary development function of finite element software.Controlling the relationship between field variables and calculation time realizes double strength reduction applicable to the layered slope.Finally,by comparing the calculation results of different examples,it is proved that the stress and displacement distribution of the critical slope state obtained by the improved method is more realistic,and the calculated safety factor is more reliable.The newly proposedmethod considers the difference of intensity attenuation between different soil layers under natural conditions and avoids the disadvantage of the strength reduction method with uniform parameters,which provides a new idea and method for stability analysis of layered and complex slopes.

Study on high and steep slope stability of surface mine based on RFPA-SRM

The instability and failure mechanism of high and steep slopes in surface mines, and the basis for some reasonable landslide prevention measures were provided using the RFPA-SRM. Based on the actual progress of the Pingzhuang Western Surface Mine and based on strength reduction method, the dynamic instability processes of the top high and steep slope was simulated. Also, the landslide mode was determined, the characteristics of the displacement distribution, the deformation, failure, and the stress distribution in the slope were demonstrated, and the stability was calculated. Conclusions can be drawn as follows： the landslide or failure of high and steep slopes on a surface mine is a gradual process, in which the slope undergoes the generation, expansion, and connection of the fractures and the displacement increases until landslide occurs; a small portion of the upper rocks fail due to the tension and the lower rocks fail due to the shear; the potential sliding surface is combined and the essential cause of the landslide is the shear stress concentration.

Rock slope stability evaluation in static and seismic conditions for left bank of Jinsha River Bridge along Lijiang-Xamgyi'nyilha railway, China

Jinsha River Bridge is located along the Lijiang-Xamgyi＇nyilha railway on the southeastern Tibet plateau; it is an area with a high prevalence of earthquakes. The bridge abutments were designed to be constructed in river bank slopes, where rocks are controlled by two sets of joint planes that significantly influence the stability of the left bank slope. According to the engineering-geological conditions and the characteristics of discontinuities, strength properties of the rock mass were obtained based on Barton model and direct shear test. Numerical analyses were performed using FLAC3D software to examine the slope＇s response to seismic loading. Then in order to evaluate the damage trends of the rock mass under the different loading conditions, a calculation model based on the geological parameters and slope stability was simulated and analyzed using the discrete element numerical simulation program UDEC （Universal Dis- tinct Element Code）, and the effect of degradation of discontinuities on the slope stability was investigated. The results show that the destruction of rock mass under the gravity, bridge foundation, and seismic load are mainly concentrated within 30 m depth of slope, and the slope under loading may slide along joint planes. In addition, the dynamic analysis by amplification of the input loading indicates that instability occurs to the bank slope at a height of about 200 m, and rock blocks will fail under seismic load. Therefore, to prevent the slope from deformation under the engineering loading and strong earthquakes, the bridge foundation should be strengthened.

Stability assessment of landslide-prone road cut rock slopes in Himalayan terrain:A finite element method based approach

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.

Slope Stability Evaluations by Limit Equilibrium and Finite Element Methods Applied to a Railway in the Moroccan Rif

Since 1930, the analysis of slope stability is done according to the limit equilibrium approach. Several methods were developed of which certain remain applicable because of their simplicity. However, major disadvantages of these methods are (1) they do not take into account the soil behavior and (2) the complex cases cannot be studied with precision. The use of the finite elements in calculations of stability has to overcome the weakness of the traditional methods. An analysis of stability was applied to a slope, of complex geometry, composed of alternating sandstone and marls using finite elements and limit equilibrium methods. The calculation of the safety factors did not note any significant difference between the two approaches. Various calculations carried out illustrate perfectly benefits that can be gained from modeling the behavior by the finite elements method. In the finite elements analysis, the shape of deformations localization in the slope is nearly circular and confirms the shape of the failure line which constitutes the basic assumption of the analytical methods. The integration of the constitutive laws of soils and the use of field’s results tests in finite elements models predict the failure mode, to better approach the real behavior of slope soil formations and to optimize its reinforcement.

APPLICATION OF STOCHASTIC FINITE ELEMENT METHOD TO STRENGTH AND STABILITY ANALYSIS OF EARTH DAMS

This paper applies the stochastic finite element method to analyse the statistics of stresses in earth dams and assess the safety and reliability of the dams. Formulations of the stochastic finite element method are briefly reviewed and the procedure for assessing dam's strength and stability is described. As an example, a detailed analysis for an actual dam Nululin dam is performed. A practical method for studying built-dams based on the prototype observation data is described.

Strength reduction and step-loading finite element approaches in geotechnical engineering

The finite element limit analysis method has the advantages of both numerical and traditional limit equilibrium techniques and it is particularly useful to geotechnical engineering.This method has been developed in China,following well-accepted international procedures,to enhance understanding of stability issues in a number of geotechnical settings.Great advancements have been made in basic theory,the improvement of computational precision,and the broadening of practical applications.This paper presents the results of research on（1） the efficient design of embedded anti-slide piles,（2） the stability analysis of reservoir slopes with strength reduction theory,and（3） the determination of the ultimate bearing capacity of foundations using step-loading FEM（overloading）.These three applications are evidence of the design improvements and benefits made possible in geotechnical engineering by finite element modeling.

Finite element reliability analysis of slope stability

The method of nonlinear finite element reliability analysis (FERA) of slope stability using the technique of slip surface stress analysis (SSA) is studied. The limit state function that can consider the direction of slip surface is given, and the formula-tions of FERA based on incremental tangent stiffness method and modified Aitken accelerating algorithm are developed. The limited step length iteration method (LSLIM) is adopted to calculate the reliability index. The nonlinear FERA code using the SSA technique is developed and the main flow chart is illustrated. Numerical examples are used to demonstrate the efficiency and robustness of this method. It is found that the accelerating convergence algorithm proposed in this study proves to be very efficient for it can reduce the iteration number greatly, and LSLIM is also efficient for it can assure the convergence of the iteration of the reliability index.

Slope analysis based on local strength reduction method and variable-modulus elasto-plastic model

Employing an ideal elasto-plastic model,the typically used strength reduction method reduced the strength of all soil elements of a slope.Therefore,this method was called the global strength reduction method(GSRM).However,the deformation field obtained by GSRM could not reflect the real deformation of a slope when the slope became unstable.For most slopes,failure occurs once the strength of some regional soil is sufficiently weakened; thus,the local strength reduction method(LSRM)was proposed to analyze slope stability.In contrast with GSRM,LSRM only reduces the strength of local soil,while the strength of other soil remains unchanged.Therefore,deformation by LSRM is more reasonable than that by GSRM.In addition,the accuracy of the slope's deformation depends on the constitutive model to a large degree,and the variable-modulus elasto-plastic model was thus adopted.This constitutive model was an improvement of the Duncan–Chang model,which modified soil's deformation modulus according to stress level,and it thus better reflected the plastic feature of soil.Most importantly,the parameters of the variable-modulus elasto-plastic model could be determined through in-situ tests,and parameters determination by plate loading test and pressuremeter test were introduced.Therefore,it is easy to put this model into practice.Finally,LSRM and the variable-modulus elasto-plastic model were used to analyze Egongdai ancient landslide.Safety factor,deformation field,and optimal reinforcement measures for Egongdai ancient landslide were obtained based on the proposed method.

Evaluation of potential surface instability using finite element method in Kharsali Village,Yamuna Valley,Northwest Himalaya

Kharsali village, located in the Northwest Himalaya near the confluence of the Yamuna River and Unta Gad, is situated on a thick(>150 m) paleolandslide deposit. The village is continuously being eroded at its base by the two rivers. Cracks are noted in most houses while the ancient Shani Temple lying to the south of the village has tilted ~5° towards the northeast. Three slope sections(S-1, S-2, S-3) were modelled and analysed to determine the displacement and shear strain patterns of the slopes. Based on surface failure conditions, potential slope instability of the Kharsali village was evaluated from 2D Finite Element Method(FEM) using Shear Strain Reduction(SSR) analysis in the Phase2 software. Results indicate a critical Stress Reduction Factor(SRF) of 1.5 for the southern edge of the village(S-1) housing the Shani Temple. The development of failure surfaces at its lower portion signifies the propagating, progressive nature of the slope. The S-2 slope section is most vulnerable to slope failure, with a critical SRF of 1.08. This has been inferred by the formation of failure surfaces with displacements of 0.05-0.08 m. The S-3 section in the northern part of the Kharsali shows highest critical SRF of 2.76. The un-metalled road section in the north of the village near S-3 hasdeveloped a failure surface with displacement of 0.003-0.004 m, and a zone of subsidence. The S-3 section is relatively stable, whereas the S-2 section is the most vulnerable portion of the village.

A New Criterion for Defining Inhomogeneous Slope Failure Using the Strength Reduction Method

A new variational method treating the system as a whole with rigorous mathematical and physical derivation was presented in this paper.Combined with classical and engineering examples,variational energy expressions of slopes were derived.In addition,the calculation programs were written in the FISH language set in FLAC3D(fast Lagrangian analysis of continua in three dimensions)software.Factors of safety(FOSs)of the models were determined by the variational method based on the strength reduction method(SRM)and then compared with other criteria or methods.The result showed that the variational method reflected the process of slope plasticity and failure uniformly and was feasible to analyze the stability of inhomogeneous slopes.The method was applicable to both two-dimensional and three-dimensional heterogeneous slopes.The small error with other criteria or methods also showed the accuracy of this method.This method unified other criteria,avoided the artificial error of other criteria,and provided a logical derivation for the instability of heterogeneous slope.This method considered the system as a whole and avoided the shortcomings of the general method of one-sided instability analysis.The proposed method is of great significance as it considers the coupling effect of stress and strain of materials and gives the mechanical basis for the instability of complex slopes.

Stress and Displacement Distribution of Soft Clay Slope with 2D and 3D Elastoplastic Finite Element Method

Based on elastoplastic model, 2D and 3D finite element method （FEM） are used to calculate the stress and displacement distribution in the soft clay slope under gravity and uniform load at the slope top. Stability analyses indicate that 3D boundary effect varies with the stress level of the slope. When the slope is stable, end effect of 3D space is not remarkable. When the stability decreases, end effect occurs; when the slope is at limit state, end effect reaches maximum. The energy causing slope failure spreads preferentially along y-z section, and when the failure resistance capability reaches the limit state, the energy can extend along x-axis direction. The 3D effect of the slope under uniform load on the top is related to the ratio of load influence width to slope height, and the effect is remarkable with the decrease of the ratio.

Optimal mining sequence for coal faces under a bedding slope:insight from landslide prevention

Repetitive mining beneath bedding slopes is identified as a critical factor in geomorphic disturbances, especially landslides and surface subsidence. Prior research has largely concentrated on surface deformation in plains due to multi-seam coal mining and the instability of natural bedding slopes, yet the cumulative impact of different mining sequences on bedding slopes has been less explored. This study combines drone surveys and geological data to construct a comprehensive three-dimensional model of bedding slopes. Utilizing FLAC3D and PFC2D models, derived from laboratory experiments, it simulates stress, deformation, and failure dynamics of slopes under various mining sequences. Incorporating fractal dimension analysis, the research evaluates the stability of slopes in relation to different mining sequences. The findings reveal that mining in an upslope direction minimizes disruption to overlying strata. Initiating extraction from lower segments increases tensile-shear stress in coal pillar overburdens, resulting in greater creep deformation towards the downslope than when starting from upper segments, potentially leading to localized landslides and widespread creep deformation in mined-out areas. The downslope upward mining sequence exhibits the least fractal dimensions, indicating minimal disturbance to both strata and surface. While all five mining scenarios maintain good slope stability under normal conditions, recalibrated stability assessments based on fractal dimensions suggest that downslope upward mining offers the highest stability under rainfall, contrasting with the lower stability and potential instability risks of upslope downward mining. These insights are pivotal for mining operations and geological hazard mitigation in multi-seam coal exploitation on bedding slopes.

A strength reduction method based on double reduction parameters and its application

In the traditional strength reduction method,the cohesion and the friction angle adopt the same reduction parameter,resulting in equivalent proportional reduction.This method does not consider the different effects of the cohesion and friction angle on the stability of the same slope and is defective to some extent.Regarding this defect,a strength reduction method based on double reduction parameters,which adopts different reduction parameters,is proposed.The core of the double-parameter reduction method is the matching reduction principle of the slope with different angles.This principle is represented by the ratio of the reduction parameter of the cohesion to that of the friction angle,described as η.With the increase in the slopeangle,ηincreases; in particular,when the slope angle is 45°,tηis 1.0.Through the matching reduction principle,different safety margin factors can be calculated for the cohesion and friction angle.In combination with these two safety margin factors,a formula for calculating the overall safety factor of the slope is proposed,reflecting the different contributions of the cohesion and friction angle to the slope stability.Finally,it is shown that the strength reduction method based on double reduction parameters acquires a larger safety factor than the classic limit equilibrium method,but the calculation results are very close to those obtained by the limit equilibrium method.

Slope stability analysis under seismic load by vector sum analysis method

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.

Two-dimensional finite element modeling of seepage-erosion coupled process within unsaturated soil slopes

Fines migration along with rainfall infiltration is a possible cause of failures of slopes composed of loose deposits.To investigate the intrinsic mechanisms,a rigid mathematical model which can fully capture the multi-phasic and multiphysical process is necessary.In this research,the macro and micro physical phenomena of fines migration process within deposited soil slopes under rainfall infiltration were summarized.Based on the mixture theory,a seepage-erosion model for unsaturated erodible soils capable to capture these phenomena mathematically was built based on a rigid theoretical framework.The model was used to simulate a set of rainfall flume tests involving fines migration phenomena with the finite element method.Two distinct slope failure modes observed experimentally,which were induced by the soil erosion-deposition properties,can be well reproduced by our numerical model.The seepage-erosion coupled process during the rainfall infiltration,as well as the intrinsic mechanism responsible for the slope failures,was illustrated in detail based on the numerical results.It was shown that the fines migration process can affect the hydro-mechanical response within unsaturated slopes significantly,and therefore special attention should be paid to those soil slopes susceptible to internal erosion.