Psuedo static stability analysis of rock slope using patton’s shear criterion

A deterministic model for the factor of safety of an idealized rock mass for planar mode of failure is developed adopting Limit Equilibrium Method (LEM) using Patton’s shear strength criterion and considering practically occurring conditions such as the effect of tension crack, water filled up in tension crack, horizontal and vertical seismic acceleration, rock bolt stabilizing force and surcharge. In the Pseudo-static analysis horizontal seismic acceleration is taken outward from the slope and vertical seismic acceleration is considered in both the direction i.e. towards the direction of gravity (downward) and opposite to the direction of gravity (upward). An expression of normal stresses as limiting criterion has been derived in order to compare the field normal stresses along the failure surface. A detailed parametric study has been presented to investigate the influence of vertical seismic coefficient for both the direction on the stability of rock slope using developed expression. For high normal stress along the failure plane, it is observed that the factor of safety decreases with increase in the value of vertical seismic coefficient towards the direction of gravity and increases linearly with increase in the value of vertical seismic coefficient against the direction of gravity and the opposite trend has been found for lower normal stress. The vertical seismic coefficient against the direction of gravity has predominant effect on factor of safety of rock slope as the rate of increase/decrease of factor of safety with vertical seismic coefficient is more against the direction of gravity. Hence in determining the critical factor of safety, effect of vertical seismic coefficient against the direction of gravity should be considered.

slopes with well-defined discontinuities and joints e.g., layered sedimentary rock, volcanic flow rock, block jointed granite and foliated metamorphic rock [3]. Planar failure also occurs when structural sliding plane dips towards the slope face with dip angle lower angle than that of slope face and higher angle than friction angle of the discontinuity surface [3,14,27].
Over the years, the stability analyses of rock slopes, susceptible to failure have been analysed numerously [5,11,12,18]. Various methods for the rock slope stability analysis are available, namely limit equilibrium method, numerical methods and kinematic methods. Limit equilibrium method is most commonly adopted method to analyse factor of safety due to its simplicity [26,27].
Hoek and Bray [6] proposed 2D analytical limit equilibrium approach for the calculation of the factor of safety of a plane failure assuming the vertical tension crack and horizontal top surface of sliding mass. Sharma et al. [24] modified the expression considering inclined surface of sliding mass and tension crack. Shukla et al. [25] derived the expression for factor of safety of anchored rock slope considering surcharge and horizontal seismic effect against plane failure. In earlier studies, importance of vertical seismic component was ignored while many recent studies suggested considering vertical seismic component as it affects stability of slopes significantly, particularly in the case of high horizontal acceleration [9,17]. Some researchers [2,8,22] have also studied the influence of vertical seismic coefficient on slope stability adjacent to plane failure and factor of safety is found to be critical considering seismic coefficient against the direction of gravity. Das and Maheshwari [2] adopted the pseudo static approach to study the influence of vertical seismic coefficient on its stability. Ahmadi and Eslami [1] used an analytical approach in which the effect of water forces was analysed.
Many researchers proposed limit equilibrium approach based on the Mohr-Coulomb's linear shear strength criterion [5,6,17,20,23,27]. Although Coulomb's shear criterion has widely been used, it is not satisfactory to consider only peak shear strength criterion for rock material. Coulomb's shear criterion has limitations as: (a) Coulomb's shear criterion implies that a major shear fracture exists at peak strength. However, it is not proven to be general for rock failure. (b) Coulomb's shear criterion also implies a specific direction of shear failure which is not always in good agreement with field observation for rock failure (c) Coulomb's shear criterion considers only the shearing during failure; however, rock generally accounts for both shearing and sliding during its failure. As a result of these reasons, other shear strength criteria are preferred for the analysis of shear failure on discontinuity such as Patton's shear criterion and Barton Bandis shear criterion. However, major limitation of Patton's shear strength criterion is that there is a strong difficulty in identifying the normal stress that determines the rupture of the rock asperities. Due to this limitation, Patton's criterion is not so popularly used for rock slope stability. To overcome this limitation, a limiting criterion is derived for the normal stresses in order to compare the field normal stresses along the failure surface in this study and stability of an idealised rock slope is analysed using Patton's criterion. Patton's shear criterion is used for the study due to its accountability of sliding and shearing phenomena during failure instead of only shearing failure considered by Coulomb's criterion.
In this study, analytical expression for the factor of safety has been developed based on Patton's shear strength criterion to consider both sliding and shearing along the failure surface under different stress levels. To develop factor of safety model, a basic geometry of rock slope is considered from Hoek and Bray [6]. Effect of practically occurring conditions such as vertical tension crack, water filled up in tension crack, rock anchor stabilising force, surcharge, horizontal and vertical seismic forces have also been considered in developing the deterministic model for factor of safety calculation associated with a planar mode of failure. Slope geometry is analysed as 2-dimensional problem considering unit thickness throughout the slope. A limiting criterion is derived for the normal stresses in order to compare the field normal stresses along the failure surface. Then two factor of safety expressions have been developed for both low and high normal stresses using the limiting criteria of stress along the failure plane. Series of parametric studies have also been conducted to analyse the influence of vertical seismic coefficient on the factor of safety.
List of symbols used in the study Limit equilibrium model Figure 1 shows an idealised rock slope of height H inclined at an angle Ѱ f with horizontal. In Fig. 1, PQRS shows an unstable rock mass of weight W with failure plane at an angle Ѱ p with horizontal, anchored rock bolt (Rock anchor stabilizing force T) inclined at an angle θ from the normal to the failure plane, carrying a surcharge q over width B at the top of rock slope. The unstable rock mass is subjected to gravitational acceleration induced by earthquake and is bounded at its back by a tension crack (RQ) of depth z filled with water up to depth z w . PS and PQ are the slope face and toe failure surface respectively. V shows the pressure exerted by filled up water in tension crack. U shows the hydrostatic uplift pressure exerted by filled-up water along the failure plane. α h and α v are the horizontal seismic coefficient and vertical seismic coefficient respectively.
Hoek and Bray [6] developed a factor of safety model of rock slope considering planar failure using Mohr-Coulomb shear strength linear criterion. Coulomb considered the smooth surface of asperities in determining the friction angle along the rock joints. In hard rock, the surface of discontinuities is never as smooth as considered by Coulomb [6]. As already mentioned, failure along the discontinuities is due to both sliding and shearing. Patton [19] considered the saw tooth surface of rock specimen to simulate the asperities and to determine the friction angle. Patton considered two failure modes, pure sliding at low normal stresses and shearing at high normal stresses. At low normal stress, sliding on inclined tooth face results in dilation of rock mass associated with rotation of blocks. Cohesion is negligible at low normal stresses and friction angle is high as the undulation at sliding surface of rock joint is accounted. With the increase in normal stresses, dilation is suppressed and shearing of intact rock commences. It results in relatively higher cohesion and lower friction angle. Patton proposed the shear strength criterion as (1) τ = σ n tan(φ b + i) (for low normal stress) Fig. 1 An anchored rock slope section with tension crack partially filled with water and all the forces acting on the slope where, σ n is the normal stress acting on the failure plane. φ b is the basic friction angle, i is the inclination angle of asperities along failure plane, c is the apparent cohesion along asperities, φ r is the residual friction angle. Shear parameters c, φ b and φ r can be obtained by direct shear test as suggested by Patton on saw tooth specimen. Inclination angle of asperities i can be measured as angle of saw tooth face of the specimen to horizontal Fig. 2.
As already mentioned, Patton's model has a limitation that there is no clear division between low and high normal stress. In this study, by comparing Eq. 1 and 2, a limiting criterion for normal stress is developed as (σ n ) lim is the limiting value of normal stress along failure plane. Based on this limiting criterion, this study considered two cases as field normal stress along the failure plane is (1) less than limit value of normal stress using Eqs. (3) and (2) more than limit value of normal stress using Eq. (3).
Factor of safety (FS) of rock slope, FS [6], F r is the summation of forces which resists sliding and F i is the summation of forces which induces sliding.

Case 1: Normal stress along the failure plane is less than the limiting value of normal stress
Resisting force F r = τ A Using Patton's shear strength criterion for low stress, Resisting force F r = [σ n tan(φ b + i)]A From Fig. 1, the total resisting force, (2) τ = c + σ n tan(φ r ) for high normal stress Resisting force(F r ) Driving force(F i ) where, γ r is the unit weight of rock and γ w is the unit weight of water. Combining Eqs. (7) and (8a, 8b, 8c and 8d) Equation (9) presents the factor of safety expression of rock slope for the considered combination of forces including the horizontal and vertical seismic inertia forces. For the ease of parametric study, this equation of factor of safety has been converted into non-dimensional form by dividing the numerator and denominator by γ r H 2 .
where,z * = z H ,q * = q γ H , z * w = z w H , γ * = γ r γ w , q * = q γ H , T * = T γ H 2 Equation (10) presents the factor of safety expression for rock slope with the considered parameters in the non-dimensional forms for low normal stress.
Case 2: Normal stress along the failure plane is more than limit value of normal stress Using Patton's shear strength criterion for high stress, Resisting force,F r = [c + σ n tan(φ r )]A From Fig. 1, the total resisting force, Combining Eqs. (4), (6) and (11) where Combining Eqs. (8a, 8b, 8c and 8d), (12) and (13) Equation (14) presents the factor of safety expression for the rock slope for the considered combination of forces for high normal stress. Again, for the ease of parametric study this equation of factor of safety has been converted into non-dimension form by dividing the numerator and denominator by γ r H 2 Equation (15) presents the factor of safety expression for the rock slope with the considered parameters in their non-dimensional forms for high normal stress. This expression can be utilized to study the influence of any particular parameter on the stability of rock slope.

Special cases
When there is no water in the tension crack, surcharge, stabilizing force and seismic force.
When there is no water in tension crack, surcharge, stabilizing force (z w = 0, q = 0, When there is no water in tension crack, surcharge and seismic forces ( for high normal stress

Illustrative example
The main influential factors for the stability of rock slopes are geometry of rock slopes (angle of slope face, angle of failure plane, height of slope, and depth of tension crack), cohesion, peak and residual angle of internal friction, inclination of rock discontinuities, drainage condition, and external forces such as rainfall, seismicity and of course the different man-made activities over the slope [1,5,7,15,24]. Characterization of these geometrical parameters is quite simple. These can be defined by field mapping. Determination of water depth and seismic acceleration is relatively difficult in rock slope evaluation because these parameters do not have single fixed value. Different parametric analyses have thus been carried out to study the influence of vertical seismic coefficient on the stability of rock slope considering the practical possible range of parameters in their non-dimensional form as given in Table 1. These ranges of parameters considered in this study are taken from the previous studies made by Hossain et al. [8] and Shukla et.al. [25]. While performing the different parametric studies, values for the parameters which are kept constant in a particular study are taken from the listed values for case 1 and case 2 in the Table 1 and for the parameter, of which the influence on factor of safety is to be studied, the listed range in the Table 1 has been considered. Two of the parameters value q* and T* are kept different in case 1 and case 2, in order to study of low stress and high stress conditions along the failure plane respectively.
Horizontal seismic coefficient is taken as per the suggestion of various authors summarised by Kramer [15]. The horizontal seismic forces are considered to act outward of the slope i.e. in the direction of failure of slope. Horizontal seismic coefficient is always positive. While the vertical seismic coefficient can be positive and negative. Positive and negative vertical seismic coefficient represents vertical seismic forces directed towards and opposite to the direction of gravity respectively. In addition to the direction, magnitude of seismic forces is also varied.

Results and discussion
In the following sections, results are presented in the graphical form for FS vs α v with varying parameters. Studies 1-9 show the analysis of the variation of factor of safety (FS) with vertical seismic coefficient for different values of horizontal seismic coefficient (α v ), angle of inclination of slope face (Ѱ f ), angle of inclination of failure plane (Ѱ p ), depth of tension crack (z*), depth of water in tension crack (z w *), unit weight of rock ( γ*), surcharge (q*), anchor stabilising force (T*), angle of inclination of rock anchor stabilising force to the normal at the failure plane (θ), in their non dimensional form for low (case 1) and high (case 2) normal stresses, respectively. Study 10 shows the analysis of the variation of factor of safety with vertical seismic coefficient for different values of basic frictional angle (φ b ) for case 1, and for different values of residual frictional angle (φ r ) for case 2. Study 11 shows the analysis of the variation of factor of safety with vertical seismic coefficient for different values of inclination angle of asperities along failure plane (i) for case 1 and for different values of cohesion (c) for case 2.

Study 1: Effect of α v on FS for different α h
In this section, as recommended by several researchers, it is studied that how the factor of safety of the considered rock slope can vary with the vertical seismic acceleration coefficient (α v ) for different value of horizontal seismic coefficient (α h ) and whether the change is significant or not. As already mentioned, Fig. 3 Table 1.
In this study, the vertical seismic coefficient is taken as not more than the half of the horizontal seismic coefficient following IS:1893-2016 [10], where, it is mentioned that the value of vertical seismic coefficient should be taken as half or two-thirds of the horizontal seismic coefficient as usually vertical motion is weaker than the horizontal motion.
It is observed from Fig. 3a that the factor of safety decreases with increase in the values of horizontal seismic coefficient significantly. At α v = 0.0, the percentage decrease in FS with α h , from 0.05 to 0.1, is nearly 10%. It is also observed that for lower value of α h less than 0.1, FS increases with an increase in the value of α v against the direction of gravity and decreases with increase in the value of α v towards the direction of gravity although the percentage change is not so significant. For higher value of α h (greater than 0.1), FS increases with an increase in the value of α v towards the direction of gravity and decreases with increase in the value of α v against the direction of gravity. As the value of α h increases, variation of FS with α v also increases. Maximum variation of FS with α v is observed at α h = 0.3. At α h = 0.3, the percentage change of relative decrease in FS with α v , from 0.0 to 0.15 against the direction of gravity is 28.2% whereas that towards the direction of gravity is 22.07%.
From Fig. 3b, it is noted that the value of FS is greater than unity when the value of α h is smaller than 0.1. Here, the value of FS decreases almost linearly with increase in α v towards the direction of gravity and increases with increase in α v against the direction of gravity, which is just opposite to the earlier case (case 1). This is due to the consideration of the sliding phenomenon along with the discontinuities of rock mass along the failure plane at low stress. But increase of factor of safety against the direction of gravity is more than the rate of decrease of factor of safety towards the direction of gravity. Variation of FS with α v is found to be maximum at α h = 0.3. At α h = 0.3, the percentage change of relative decrease in FS with α v from 0.0 to 0.15 towards the direction of gravity is 27.27%, whereas the percentage change of relative increase in FS with α v from 0.0 to 0.15 against the direction of gravity is 35.13%.

Study 2: Effect of α v on FS for different Ѱ f
In this study, behaviour of rock slope with vertical seismic coefficient for different values of angle of slope face is of concern. As reported by many researchers [7,25], the angle of rock slope face affects the slope stability to a great extent due to the fact that factor of safety is inversely proportional to angle of slope face. Figure 4 presents the results for different value of Ѱ f as 40°, 45°, 50°, 55°, keeping other parameters as mentioned in Table 1.
From Fig. 4a, it is observed that FS reduces as Ѱ f increases. The rate of reduction in FS is more when Ѱ f from 40° to 45° and less when Ѱ f greater than 45°. For Ѱ f = 40°, it is observed that FS reduces with an increase in α v towards the direction of gravity and with the increase in α v against the direction of gravity. For greater value of Ѱ f (greater than 40°), FS reduces with increase in α v against the direction of gravity although the variation is insignificant. Maximum variation of factor of safety with α v is observed at Ѱ f = 40°. At Ѱ f = 40°, the percentage change of relative increase in FS with α v from 0.0 to 0.10 against the direction of gravity is 56.33% whereas that towards the direction of gravity is 38.22%.
From Fig. 4b, it is observed that the factor of safety is very high for Ѱ f = 40°. FS decreases as Ѱ f increases. FS decreases significantly when Ѱ f increases from 40° to 45°. For the value of Ѱ f = 40°, FS varies nonlinearly with the change in the value of α v and decreases with increase in α v towards the direction of gravity and increases with increase in α v against the direction of gravity. For Ѱ f greater than 45°, FS varies almost linearly with α v and the rate of change of FS with α v is almost insignificant. For Ѱ f greater than 50°, FS becomes less than unity making the rock slope unstable. Maximum variation of FS with α v is observed at Ѱ f = 40°. At Ѱ f = 40°, the percentage change of relative decrease in FS with α v from 0.0 to 0.1 towards the direction of gravity is 65.13% and that against the direction of gravity is 22.30%.
From this study, it can be stated as factor of safety of rock slope decreases with increase in angle of inclination of slope face with horizontal as steepness of slope face makes the slope unstable.

Study 3: Effect of α v on FS for different Ѱ p
This study focuses on the influence of vertical seismic coefficients on factor of safety for different values of angle of inclination of slope failure plane (Ѱ p ). Figure 5 presents the results for different Ѱ p as 45°, 42.5°, 40°, 37.5°, 35 o considering specific set of influential parameters in their non-dimensional form as mentioned in Table 1.
From Fig. 5a it is observed that FS decreases as angle of inclination of slope failure decreases. As Ѱ p decreases, it means that the weight of sliding rock mass increases. For the lower value of Ѱ p , FS increases as α v increases towards the direction of gravity and decreases as α v increases against the direction of gravity. However, the rate of variation of FS with α v is almost insignificant.
From Fig. 5b, it is observed that FS decreases as α v increases towards the direction of gravity and increases as α v increases against the direction of gravity. For Ѱ p = 45°, variation of FS with α v is nonlinear and rate of increase with α v against the direction of gravity is much higher than the rate of decrease with α v towards the direction of gravity. FS is greater than unity for all the values of Ѱ p greater than 40°. For Ѱ p less than 40°, the rate of change of FS with vertical seismic coefficient is again insignificant. Variation of FS is found to be maximum at Ѱ p = 45°. At Ѱ p = 45°, the percentage change of relative decrease in FS with α v from 0.0 to 0.1 towards the direction of gravity is 26.83% and that against the direction of gravity is 60.27%.
From this study it is observed that factor of safety of rock slope increases with increase in angle of slope failure plane with horizontal as the weight of rock mass along the sliding plane reduces with increase in inclination angle of failure plane to horizontal.

Study 4: Effect of α v on FS for different z*
This study presents the variation of factor of safety with vertical seismic coefficients for different value of depth of tension crack (z*). Figure 6 presents the results for different z* as 0.05, 0.10, 0.15, 0.20, 0.25 considering a specific set of influential parameters in their non-dimensional form as mentioned in Table 1 From Fig. 6a, it is observed that FS increases as the depth of tension crack increases. It is due to geometrical changes in the sliding rock mass over failure plane.
FS increases with an increase in α v towards the direction of gravity and decreases with increase in α v against the direction of gravity although the amount of variation of FS with α v is not so significant. But the rate of decrease in FS with increase in α v against the direction of gravity is more than the rate of increase in FS with increase in α v towards the direction of gravity. Variation of FS with α v is found to be maximum at z* = 0.05. At z* = 0.05, the percentage change of relative decrease in FS with α v from 0.0 to 0.15 against the direction of gravity is 16.3% whereas that towards the direction of gravity is 13.7%.
From Fig. 6b, it is observed that FS decreases with an increase in the value of α v towards the direction of gravity and increases with increase in α v against the direction of gravity almost linearly for the considered value of depth of tension crack. Here again FS increases as the depth of tension crack (z*) increases. FS is greater than unity for z* greater or equals to 0.25. Variation of FS with α v increases as z* increases. For z* = 0.25, the percentage change of relative decrease in FS with α v from 0.0 to 0.10 towards the direction of gravity is 66.63% and that against the direction of gravity is 81.88%.

Study 5: Effect of α v on FS for different z w *
This study presents the variation of factor of safety with vertical seismic coefficient (α v ) for different value of depth of water in tension crack (z w *). This study is conducted to know how different depth of water in tension crack affects in FS with different vertical acceleration coefficient. Figure 7 presents the results for different z w * as 0.0, 0.025, 0.05, 0.075, 0.10 considering a specific set of influential parameters in their non-dimensional form as mentioned in Table 1.
It is observed from Fig. 7 that, FS reduces with an increase in depth of water in tension crack. Water present in tension crack forces the surface of discontinuity apart, as hydrostatic pressure act in the direction of movement of sliding block. From Fig. 7a, it is also observed that the FS increases with increase in α v towards the direction of gravity and decreases against the direction of gravity almost linearly. Variation of FS is found to be maximum at z w = 0.1. At z w = 0.1, the percentage change of relative decrease in FS with α v from 0.0 to 0.15 against the direction of gravity is 21.3% whereas the percentage change of relative increase in FS with α v from 0.0 to 0.15 that towards the direction of gravity is 17.8%.
From Fig. 7b, it is observed that FS is higher for smaller value of z w *. FS decreases with increase in α v towards the direction of gravity and increases with increase in α v against the direction of gravity almost linearly. At z w * = 0.0, the percentage change of relative decrease in FS with α v from 0.0 to 0.10 towards the direction of gravity is 44.4% and the percentage change of relative increase in FS with α v from 0.0 to 0.10 against the direction of gravity is 53.7%.

Study 6: Effect of αvon FS for different γ*
In this study influence of unit weight of rock with varying vertical seismic coefficient on factor of safety has been presented. Figure 8 presents the result for different values of unit weight of rock ( γ*) as 2.0, 2.25, 2.5, 2.75, 3.0 considering a specific set of influential parameters in their non-dimensional form as mentioned in Table 1.
From Fig. 8a, it is observed that the variation of unit weight of rock has no significant effect on FS. From Fig. 8a it is observed FS reduces with an increase in α v against the direction of gravity and increases with α v towards the direction of gravity. But the variation of FS with α v is almost insignificant. For γ* = 2.0, the percentage change of relative increase in FS is 16.8% with increase in α v towards the direction of gravity and that is 14% against the direction of gravity. As γ * increases, the variation of FS with α v increases in either direction of gravity. Variation of FS with α v is found to be maximum at higher value of γ* = 3.0. At γ* = 3.0, the percentage change of relative decrease in FS with α v from 0.0 to 0.1 against the direction of gravity is 21.3% whereas that towards the direction of gravity is 17.88%.
From Fig. 8b Factor of safety decreases as α v increases towards the direction of gravity whereas it increases with increase in α v against the direction of gravity. Variation of FS with α v is observed to be more for higher value of γ *. At γ* = 3.0, the percentage change of relative decrease in FS with α v from 0.0 to 0.10 towards the direction of gravity is 43.2% and that against the direction of gravity is 52.4%.

Study 7: Effect of α v on FS for different q*
Study 7 presents the behaviour of rock slope in terms of factor of safety with the application of different amount of surcharge with varying vertical seismic coefficient. Figure 9 presents the result for different magnitude of surcharge (q*) considering a specific set of influential parameters in their non-dimensional form as mentioned in Table 1.  Fig. 9a, it is observed that FS decreases with an increase in surcharge q* as surcharge directly adds weight to the sliding rock mass consequently increases the shear stress along the failure plane. Variation of FS with α v is almost linear. FS increases with increase in the value of α v towards the direction of gravity whereas FS decreases with increase in the value of α v against the direction of gravity. But the rate of increase of FS with α v towards the direction of gravity is less than the rate of decrease of FS with α v against the direction of gravity. However, this variation of FS with α v is not so significant and maximum variation is observed at q* = 0.3. At q* = 0.3, the percentage change of relative decrease in FS with α v from 0 to 0.1 is only 1.89% against the direction of gravity and that towards the direction of gravity is 1.59%.
From Fig. 9b, it is observed that for lower value of q*, variation of FS is more with α v and for higher value of q*, variation of FS is less with α v in both of the direction of gravity and observed maximum at q* = 0.25. Variation of FS with α v is almost linear. For q* greater than 0.5, with the variation in vertical seismic coefficient, FS does not change significantly in both of the direction of gravity. FS reduces with increase in surcharge (q*). When q* varies from 0.0 to 0.5, the rate of decrease of FS is more than that when q* varies from 0.5 to 1.5. At q* = 0.25, the percentage change of relative decrease in FS with α v from 0.0 to 0.10 towards the direction of gravity is 66.8% and that against the direction of gravity is 82%.

Study 8: Effect of α v on FS for different T*
Study 8 presents the influence of rock anchor stabilising force (T*) with varying vertical seismic coefficient on factor of safety of rock slope. Figure 10 presents the results for different value of rock anchor stabilizing force (T*) as 0.0, 0.05, 0.1, 0.15, 0.20 considering a specific set of influential parameters in their nondimensional form as mentioned in Table 1.
From Fig. 10, it is observed that FS increases with increase in rock anchor stabilising force significantly. From Fig. 10a, it is observed that when T* greater than 0.1, FS becomes more than unity at all the considered value of α v hence slope becomes stable. For lower value of T* less than 0.1, FS increases with increase in the value of α v towards the direction of gravity and decreases with increase in the value of α v against the direction of gravity. But for higher value of T*, FS decreases with increase in α v towards the direction of gravity and increases with increase in α v against the direction of gravity and the maximum variation of FS is observed at T* = 0.15. At T* = 0.15, the percentage change of relative decrease in FS with α v , from 0.0 to 0.1 towards the direction of gravity is 75.8% whereas that against the direction of gravity is 99.9%.
From Fig. 10b, it is observed that FS decreases with the increase in α v towards the direction of gravity and increases with the increase in α v against the direction of gravity. But variation in FS with α v for lower value of T* is less, more for higher value of T* and maximum at T* = 0.15. At T* = 0.15, the percentage change of relative increase in FS with α v from 0.0 to 0.10 against the direction of gravity is 9.74% and the percentage change of relative decrease in FS with α v from 0.0 to 0.10 towards the direction of gravity is 7.58%.

Study 9: Effect of α v on FS for different θ
Study 9 presents the influence of angle of inclination of the rock anchor stabilising force to the normal at the failure plane (θ) with varying vertical seismic coefficient (α v ) on factor of safety. Figure 11 presents the results for different value of θ as 0°, 20°, 40°, 60°, 80° considering a specific set of influential parameters in their non-dimensional form as mentioned in Table 1.
It is observed from Fig. 11a that the Factor of safety increases with increase in the value of θ from 0° to 40°, then becomes nearly constant up to 60° and then decreases when θ more than 60°. FS increases with increase in the value of α v towards the direction of gravity and decreases with increases in the value of α v against the direction of gravity. But rate of decrease in FS with α v against the direction of gravity is more than the rate of increase of FS with α v towards the direction of gravity. As θ increases from 0 o to 60°, variation of FS with α v becomes lesser, but it is found highest at θ = 80°. At θ = 80°, the percentage change of relative increase in FS with α v from 0.0 to 0.1 towards the direction of gravity is 1.69% and that against the direction of gravity is 2.05%.
From Fig. 11b, it is observed that with increase in value of θ, FS increases. FS decreases with increase in α v towards the direction of gravity and increases with increase in α v against the direction of gravity. With variation of θ from 0° to 40 o the rate of increase of FS is more. Variation of FS with α v is observed to be maximum at θ = 80°. At θ = 80°, the percentage change of relative decrease in FS with α v from 0.0 to 0.10 towards the direction of gravity is 50.3% and that against the direction of gravity is 62.2%.

Study 10: Effect of α v on FS for different φ b and φ r
In this study, the stability condition of rock slope for different basic frictional angle and residual frictional angle with varying vertical seismic coefficient have been studied. Figure 12 presents variation of factor of safety (FS) of rock slope with vertical seismic coefficient (α v ) for different values of (a) basic friction angle (φ b ) for case 1 and (b) residual friction angle (φ r ) for case 2, considering a specific set of influential parameters in their non-dimensional form as mentioned in Table 1.
From Fig. 12a, it is observed that FS increases significantly with increase in the angle of friction angle. For φ b greater than 30°, FS becomes more than unity. Hence slope becomes stable. FS increases with increase in the value of α v towards the direction of gravity and decreases with increase in the value of α v against the direction of gravity. This variation of FS with α v is more for higher value of φ b . At φ b = 80°, the percentage change of relative increase in FS with α v from 0.0 to 0.1 towards the direction of gravity is 2.19% and that against the direction of gravity is 2.64%.
From Fig. 12b, it is observed that FS increases with an increase in the value of φ r . Rate of increase of FS with φ r remains nearly constant. It is because of the linear behaviour of Patton's model of shear strength for high normal stress. As φ r increases, the rate of variation of FS decreases but that is not very significant. FS decreases with increase in α v towards the direction of gravity and increases with increase in α v against the direction of gravity almost linearly. Variation of FS with α v is observed to be maximum at φ r = 5°. At φ r = 5°, the percentage change of relative decrease in FS with α v from 0.0 to 0.10 towards the direction of gravity is 44.7% and the percentage change of relative increase in FS with α v from 0.0 to 0.10against the direction of gravity is 53.6%.

Study 11: Effect of α v on FS for different i and c*
In this Study, the influence of inclination angle of asperities along the joint surface and cohesion with varying vertical seismic coefficient on rock slope stability have been observed.  Figure 13 presents the results for different value of (a) inclination angle (i) as 0°, 5°, 10°, 15°, 20°, and (b) cohesion (c*) as 0.0, 0.04, 0.08, 0.12, 0.16 considering a specific set of influential parameters in their non-dimensional form as mentioned in Table 1.
From Fig. 13a, it is observed that FS increases significantly as angle of inclination of jointed rock increases. FS increases with increase in the value of α v towards the direction of gravity and decreases with increase in the value of α v against the direction of gravity. Variation of FS with α v is linear. Rate of decrease of FS with α v against the direction of gravity is more than the rate of increase of FS with α v towards the direction of gravity. As the inclination angle of jointed rock increases, variation of FS with α v increases in both the direction of gravity. At i = 20°, the percentage change of relative increase in FS with α v from 0.0 to 0.1 towards the direction of gravity is 1.83% and that against the direction of gravity is 2.21%, thus the variation is not so significant.
From Fig. 13b, it is observed that FS increases with increase in cohesion. FS decreases with increase in α v towards the direction of gravity and increases with increase in α v against the direction of gravity almost linearly. As cohesion increases, rate of variation of FS with α v increases and observed to be maximum at c* = 0.16. At c* = 0.16, the percentage change of relative decrease in FS with α v from 0.0 to 0.10 towards the direction of gravity is 69.8% and that against the direction of gravity is 84.6%.

Conclusion
This paper presents the results of a number of pseudo-static seismic stability analyses of an idealised rock slope, considering a planar failure mechanism and adopting the Limit Equilibrium Method (LEM) and the Patton's shear strength criterion. In the analyses, some actually occurring conditions were considered, i.e. the presence of a tension crack at the slope back, hydrostatic pressures acting along both the basal failure plane and the tension crack, surcharge and stabilising rock bolt. Also, an expression of normal stresses as limiting criterion has been derived in order to compare the field normal stresses along the failure surface. Based on the limiting criterion, expressions for factor of safety for low and high normal stresses along failure plane are developed. The pseudo-static seismic stability analyses were carried out in order to evaluate the influence of the vertical seismic acceleration on the factor of safety of the slope, for various reference configurations.
Based on the studies undertaken in this paper, the following concluding remarks can be made: 1. For high normal stress along failure plane, it is observed that the factor of safety decreases with increase in the value of vertical seismic coefficient towards the direction of gravity and increases linearly with increases in the value of vertical seismic coefficient against the direction of gravity. In the case of high normal stress along the failure plane, driving mechanism is gravitational body force, since the vertical seismic force affects the normal stress, and hence the shear strength along the failure slope. when the vertical seismic acceleration increases the slope becomes potentially unstable. Opposite trend of factor of safety with vertical seismic coefficient have been found for low normal stress. It is due to occurrence of pure sliding phenomena along the discontinuities of rock mass along the failure plane at low stress. Vertical seismic acceleration has to be considered in the analysis especially when the horizontal seismic acceleration is high. Effect of vertical seismic coefficient against the direction of gravity has predominant effect on factor of safety of rock slope as the rate of increase/decrease of factor of safety with vertical seismic coefficient is more against the direction of gravity. Hence in determining the critical factor of safety, the vertical seismic coefficient against the direction of gravity should be considered. 2. Factor of safety of rock slope reduces with increase in horizontal seismic coefficients.
Reduction in factor of safety appears to be in proportional to the increase in horizontal seismic coefficient. 3. Factor of safety of rock slope increases as the depth of tension crack increases. It is due to geometrical changes in the sliding rock mass over failure plane. Results show that the presence of water in tension crack can reduce the factor of safety significantly. In practical situation, measurement of depth of tension crack and water depth present in it may be difficult. Uncertainty of depth of tension crack and water present in it may have significant effect on stability of rock slope. Probabilistic analysis might be the better way to explore the effect of depth of tension crack and presence of water in tension crack on the stability of rock slope. 4. Unit weight of rock does not have significant effect on factor of safety of rock slope.
But over the long span of life, to quantify the effect of unit weight, sensitivity analysis using probabilistic techniques would be the better approach. In the rock slope analysis, strength primarily depends on the cohesion and basic and residual frictional angles and increase in factor of safety appears to be proportional to increase in cohesion and friction angle. 5. Presence of surcharge reduces factor of safety significantly as surcharge directly adds weight to the sliding rock mass consequently increases the shear stress along the failure plane. Result shows that rate of decrement in factor of safety is higher for low value of surcharge and lesser for higher value of surcharge with the increase in seismic forces. Hence limiting the surcharge over rock slope can improve the safety of slope significantly. 6. Rock anchor stabilising force with angle of inclination of the rock anchor stabilising force to the normal (θ) to failure plane has important role in making the slope stable. Factor of safety increases with increase in rock anchor stabilising force but it increases with θ up to approximately 60° then it decreases beyond 60°. 7. Factor of safety of rock slope increases with residual friction angle and cohesion. Rate of increase of factor of safety of rock slope increases with increase in residual friction angle but it remains constant with cohesion. 8. In case of factor of safety calculation for low normal stress, FS increases with basic friction angle and angle of inclination of asperities along the basal discontinuity with horizontal. Undulations along the discontinuities has significant influence on shear behaviour.
For the stability analysis of rock slope, the Psuedo-static analysis has been adopted which does not account for the cyclic nature of loading. The effect of cyclic loading will be considered in future studies.