---

Existence of discontinuities causes higher deformability and lower strength in rock masses. Thus joints can change the rock mass behaviour due to the applied loads. For this reason properties and orientation of the joint sets have a great effect on the stability of rock slopes. In this paper, after introducing some numerical methods for evaluating the factor of safety for the stability of slopes, stability of jointed rock slopes in the plane strain condition is investigated with the strength reduction technique this method is modified and applied in the multilaminate framework. First of all, stability of one homogeneous rock slope is investigated and compared with the limit equilibrium method. Then stability of a layered rock slope is analyzed with some modifications in the strength reduction technique. Effects of orientation, tensile strength and dilation of layered joint sets on the factor of safety and location of the sliding block are explained.

---

Presented is a method of three-dimensional stability analysis of convex slopes in plan view based on the Lower-bound theorem of
the limit analysis approach. The method’s aim is to determine the factor of safety of such slopes using numerical linear finite
element and lower bound limit analysis method to produce some stability charts for three dimensional (3D) homogeneous convex
slopes. Although the conventional two and three dimension limit equilibrium method (LEM) is used more often in practice for
evaluating slope stability, the accuracy of the method is often questioned due to the underlying assumptions that it makes. The
rigorous limit analysis results in this paper together with results of other researchers were found to bracket the slope stability
number to within ±10% or better and therefore can be used to benchmark for solutions from other methods. It was found that using
a two dimensional (2D) analysis to analyze a 3D problem will leads to a significant difference in the factors of safety depending
on the slope geometries. Numerical 3D results of proposed algorithm are presented in the form of some dimensionless graphs which
can be a convenient tool to be used by practicing engineers to estimate the initial stability for excavated or man-made slopes

---

Seismic stability of slopes is typically evaluated by conventional methods under the assumption that the slope is subjected to an

earthquake just for one time. In general, time histories of loadings on slopes are unknown and loads are of variable repeated

nature. Shakedown phenomenon can be considered as a safe state for slopes subjected to variable repeated loadings. In this study,

lower bound dynamic shakedown theorem is employed for the seismic stability of slopes as a comprehensive verification. A

numerical method applied previously to evaluate roads under the traffic loads was modified to make it appropriate for dynamic

shakedown analysis in the present study. The numerical method is based on the combination of finite element and linear

programming methods. Critical PGA is employed as a comparative parameter to compare shakedown and pseudostatic methods.

Results show that, unlike pseudostaic method, shakedown approach is able to consider dynamic properties of load and slope.

Also, it is indicated that contrary to pseudostaic approach, shakedown solutions are different for slopes and embankments.

Shakedown and pseudostaic critical PGA versus dynamic properties of load and slope creates four distinct zones. It is shown that

the forgoing zones can be used as appropriate tools for seismic zonation of slopes based on their short term and long term safety

---

The purpose of this study is to compare the results of different probabilistic methods such as the perturbation method, Stochastic Finite Element Method (SFEM) and Monte Carlo Method. These methods were used to study the convergence of direct approach for slope stability analysis and are developed for a linear soil behavior. In this study, two dimensional random fields are used and both the First Order Reliability Method (FORM) and Limited Step Length Iteration Method (LSLIM) have been adopted to evaluate the reliability index. The study found that the perturbation method of the second order is easy to apply using the field’s theory because accuracy is reached even with different coefficients of variation of input variables, while the spectral finite element method yields accurate results only for high levels of solution development.

---

The paper present the results of various experimental and numerical studies on slopes, small scale physical modeling of slope under surcharge load were performed in loose sand environment. Digital images were captured during the incremental loading from side of model simultaneously. The Particle Image Velocimetry (PIV) and 3D numerical model was applied to investigate the slope under surcharge loading and some of the other important factors that control the performance of piles are investigated. The factors of safety and location of critical failure surfaces of reinforced and unreinforced slopes obtained and compared for various slopes. For homogenous slope, it is found for stabilized slope with piles, the 3D failure surface shape in front of piles is triangle, unlike its conical shape in laterally loaded piles. The paper also studies numerically the effect of soft bound effect, soil properties, pile spacing, pile position and surcharge distance effects on stabilized and none stabilized slopes. The results indicate that the depth of soft soil layer from ground surface and its angles with horizontal direction has importance effect on optimum location of stabilized pile. Studies on pile space effects shows with increasing pile space, arching phenomenon didn’t developed and partial pressure of supported portion transferred to un supported soil portion and the potential failure volume of the slope becomes large.