Slope Stability Analysis (slope + stability_analysis)

Distribution by Scientific Domains


Selected Abstracts


Slope stability analysis based on elasto-plastic finite element method

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 14 2005
H. Zheng
Abstract The paper deals with two essential and related closely processes involved in the finite element slope stability analysis in two-dimensional problems, i.e. computation of the factors of safety (FOS) and location of the critical slide surfaces. A so-called ,,v inequality, sin ,,1 , 2v is proved for any elasto-plastic material satisfying Mohr,Coulomb's yield criterion. In order to obtain an FOS of high precision with less calculation and a proper distribution of plastic zones in the critical equilibrium state, it is stated that the Poisson's ratio v should be adjusted according to the principle that the ,,v inequality always holds as reducing the strength parameters c and ,. While locating the critical slide surface represented by the critical slide line (CSL) under the plane strain condition, an initial value problem of a system of ordinary differential equations defining the CSL is formulated. A robust numerical solution for the initial value problem based on the predictor,corrector method is given in conjunction with the necessary and sufficient condition ensuring the convergence of solution. A simple example, the kinematic solution of which is available, and a challenging example from a hydraulic project in construction are analysed to demonstrate the effectiveness of the proposed procedures. Copyright © 2005 John Wiley & Sons, Ltd. [source]


Geophysical investigation and dynamic modelling of unstable slopes: case-study of Kainama (Kyrgyzstan)

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2008
G. Danneels
SUMMARY The presence of massive Quaternary loess units at the eastern border of the Fergana Basin (Kyrgyzstan, Central Asia) makes this area particularly prone to the development of catastrophic loess earthflows, causing damages and injuries almost every year. Efficient disaster management requires a good understanding of the main causes of these mass movements, that is, increased groundwater pressure and seismic shaking. This paper focuses on the Kainama earthflow, mainly composed of loess, which occurred in 2004 April. Its high velocity and the long run-out zone caused the destruction of 12 houses and the death of 33 people. In summer 2005, a field survey consisting of geophysical and seismological measurements was carried out along the adjacent slope. By combination and geostatistical analysis of these data, a reliable 3-D model of the geometry and properties of the subsurface layers, as shown in the first part of the paper, was created. The analysis of the seismological data allowed us to point out a correlation between the thickness of the loess cover and the measured resonance frequencies and associated amplification potential. The second part of this paper is focused on the study of the seismic response of the slope by numerical simulations, using a 2-D finite difference code named FLAC. Modelling of the seismic amplification potential along the slope confirmed the results obtained from the seismological survey,strong amplifications at the crest and bottom of the slope where there is a thick loess cover and almost no amplification in the middle part of the slope. Furthermore, dynamic slope stability analyses were conducted to assess the influence of local amplifications and increased groundwater pressures on the slope failure. The results of the dynamic modelling, although preliminary, show that a combination of seismic and hydrologic origin (pore pressure build-up during the seismic shaking) is the most probable scenario responsible for the 2004 failure. [source]


Unsaturated slope stability analysis with steady infiltration or evaporation using elasto-plastic finite elements

INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 3 2005
D. V. Griffiths
Abstract The paper presents results of unsaturated slope stability analyses using elasto-plastic finite elements in conjunction with a novel analytical formulation for the suction stress above the water table. The suction stress formula requires four parameters, three for the soil type and one for the steady infiltration (or evaporation) due to environmental effects. The suction stress approach enables the analysis to proceed in the context of classical effective stress, while maintaining the advantages of a general non-linear finite element approach in which no advance assumptions need to be made about the shape or location of the critical failure surface. The results show the extent to which suctions above the water table can increase the factor of safety of a slope for a variety of different soil types and infiltration rates. All stability analyses that include the effects of suction stresses are contrasted with more traditional approaches in which water pressures above the water table are ignored. Copyright © 2005 John Wiley & Sons, Ltd. [source]


On the boundary conditions in slope stability analysis

INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 11 2003
Ashok K. Chugh
Abstract Boundary conditions can affect computed factor of safety results in two- and three-dimensional stability analyses of slopes. Commonly used boundary conditions in two- and three-dimensional slope stability analyses via limit-equilibrium and continuum-mechanics based solution procedures are described. A sample problem is included to illustrate the importance of boundary conditions in slope stability analyses. The sample problem is solved using two- and three-dimensional numerical models commonly used in engineering practice. Copyright © 2003 John Wiley & Sons, Ltd. [source]


A low-dimensional physically based model of hydrologic control of shallow landsliding on complex hillslopes

EARTH SURFACE PROCESSES AND LANDFORMS, Issue 13 2008
Ali Talebi
Abstract Hillslopes have complex three-dimensional shapes that are characterized by their plan shape, profile curvature of surface and bedrock, and soil depth. To investigate the stability of complex hillslopes (with different slope curvatures and plan shapes), we combine the hillslope-storage Boussinesq (HSB) model with the infinite slope stability method. The HSB model is based on the continuity and Darcy equations expressed in terms of storage along the hillslope. Solutions of the HSB equation account explicitly for plan shape by introducing the hillslope width function and for profile curvature through the bedrock slope angle and the hillslope soil depth function. The presented model is composed of three parts: a topography model conceptualizing three-dimensional soil mantled landscapes, a dynamic hydrology model for shallow subsurface flow and water table depth (HSB model) and an infinite slope stability method based on the Mohr,Coulomb failure law. The resulting hillslope-storage Boussinesq stability model (HSB-SM) is able to simulate rain-induced shallow landsliding on hillslopes with non-constant bedrock slope and non-parallel plan shape. We apply the model to nine characteristic hillslope types with three different profile curvatures (concave, straight, convex) and three different plan shapes (convergent, parallel, divergent). In the presented model, the unsaturated storage has been calculated based on the unit head gradient assumption. To relax this assumption and to investigate the effect of neglecting the variations of unsaturated storage on the assessment of slope stability in the transient case, we also combine a coupled model of saturated and unsaturated storage and the infinite slope stability method. The results show that the variations of the unsaturated zone storage do not play a critical role in hillslope stability. Therefore, it can be concluded that the presented dynamic slope stability model (HSB-SM) can be used safely for slope stability analysis on complex hillslopes. Our results show that after a certain period of rainfall the convergent hillslopes with concave and straight profiles become unstable more quickly than others, whilst divergent convex hillslopes remain stable (even after intense rainfall). In addition, the relation between subsurface flow and hillslope stability has been investigated. Our analyses show that the minimum safety factor (FS) occurs when the rate of subsurface flow is a maximum. In fact, by increasing the subsurface flow, stability decreases for all hillslope shapes. Copyright © 2008 John Wiley & Sons, Ltd. [source]


Unsaturated slope stability analysis with steady infiltration or evaporation using elasto-plastic finite elements

INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 3 2005
D. V. Griffiths
Abstract The paper presents results of unsaturated slope stability analyses using elasto-plastic finite elements in conjunction with a novel analytical formulation for the suction stress above the water table. The suction stress formula requires four parameters, three for the soil type and one for the steady infiltration (or evaporation) due to environmental effects. The suction stress approach enables the analysis to proceed in the context of classical effective stress, while maintaining the advantages of a general non-linear finite element approach in which no advance assumptions need to be made about the shape or location of the critical failure surface. The results show the extent to which suctions above the water table can increase the factor of safety of a slope for a variety of different soil types and infiltration rates. All stability analyses that include the effects of suction stresses are contrasted with more traditional approaches in which water pressures above the water table are ignored. Copyright © 2005 John Wiley & Sons, Ltd. [source]


On the boundary conditions in slope stability analysis

INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 11 2003
Ashok K. Chugh
Abstract Boundary conditions can affect computed factor of safety results in two- and three-dimensional stability analyses of slopes. Commonly used boundary conditions in two- and three-dimensional slope stability analyses via limit-equilibrium and continuum-mechanics based solution procedures are described. A sample problem is included to illustrate the importance of boundary conditions in slope stability analyses. The sample problem is solved using two- and three-dimensional numerical models commonly used in engineering practice. Copyright © 2003 John Wiley & Sons, Ltd. [source]


Inter-relations between experimental and computational aspects of slope stability analysis

INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 5 2003
R. Baker
Abstract Most conventional slope stability calculations are based on the linear Mohr,Coulomb failure criterion. However, a substantial amount of experimental evidence suggests that failure criteria of many soils are not linear particularly in the range of small normal stresses. This departure from linearity is significant for slope stability calculations since for a wide range of practical stability problems, critical slip surfaces are shallow and normal stresses acting on such surfaces are small. There exists a technical difficulty in performing strength measurements in the range of small normal stresses relevant to such slope stability problems. As a result, in many practical situations strength measurements are performed at much larger normal stresses then those relevant for the stability problem under consideration. When this is the case, use of the Mohr,Coulomb criterion amounts to a linear extrapolation of experimental information (obtained at large normal stresses), into the range of small normal stresses, which is relevant to the problem. This extrapolation results with very significant overestimation of calculated safety factors in cases when there is large mismatch between experimental and relevant ranges of normal stresses. The present work delineates the extent of this problem and suggests a practical way to overcome it. Copyright © 2003 John Wiley & Sons, Ltd. [source]


The use of an SQP algorithm in slope stability analysis

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 1 2005
Jian Chen
Abstract In the upper bound approach to limit analysis of slope stability based on the rigid finite element method, the search for the minimum factor of safety can be formulated as a non-linear programming problem with equality constraints only based on a yield criterion, a flow rule, boundary conditions, and an energy-work balance equation. Because of the non-linear property of the resulting optimization problems, a non-linear mathematical programming algorithm has to be employed. In this paper, the relations between the numbers of nodes, elements, interfaces, and subsequent unknowns and constraints in the approach have been derived. It can be shown that in the large-scale problems, the unknowns are subject to a highly sparse set of equality constraints. Because of the existence of non-linear equalities in the approach, this paper applies first time a special sequential quadratic programming (SQP) algorithm, feasible SQP (FSQP), to obtain solutions for such non-linear optimization problems. In FSQP algorithm, the non-linear equality constraints are turned into inequality constraints and the objective function is replaced by an exact penalty function which penalizes non-linear equality constraint violations only. Three numerical examples are presented to illustrate the potentialities and efficiencies of the FSQP algorithm in the slope stability analysis. Copyright © 2004 John Wiley & Sons, Ltd. [source]


Slope stability analysis based on elasto-plastic finite element method

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 14 2005
H. Zheng
Abstract The paper deals with two essential and related closely processes involved in the finite element slope stability analysis in two-dimensional problems, i.e. computation of the factors of safety (FOS) and location of the critical slide surfaces. A so-called ,,v inequality, sin ,,1 , 2v is proved for any elasto-plastic material satisfying Mohr,Coulomb's yield criterion. In order to obtain an FOS of high precision with less calculation and a proper distribution of plastic zones in the critical equilibrium state, it is stated that the Poisson's ratio v should be adjusted according to the principle that the ,,v inequality always holds as reducing the strength parameters c and ,. While locating the critical slide surface represented by the critical slide line (CSL) under the plane strain condition, an initial value problem of a system of ordinary differential equations defining the CSL is formulated. A robust numerical solution for the initial value problem based on the predictor,corrector method is given in conjunction with the necessary and sufficient condition ensuring the convergence of solution. A simple example, the kinematic solution of which is available, and a challenging example from a hydraulic project in construction are analysed to demonstrate the effectiveness of the proposed procedures. Copyright © 2005 John Wiley & Sons, Ltd. [source]