Evaluation on the Stability of Slope at Faculty of Engineering and Built Environment ( FKAB ) using Slope / w

Slope stability topics have become one of the main issues due to the nature of topography and weather conditions in Malaysia. This study examines the stability of the slope at the Faculty of Engineering and Built Environment (FKAB) and eventually to identify the best approach for the analysis. The three main objectives of this study are to evaluate the slope of FKAB using the Limit Equilibrium method, using Slope/w to evaluate the critical soil parameters obtain through back analysis, and eventually to obtain the best method to be use for further analysis on the stability of slope. The steps involve in solving the model of the slope were drawn, defined and assigned accordingly before the final analysis is being carried out in accordance to the methodology of the Slope/w 2007 software. The minimum Factor of Safety (FOS) obtained from the analysis using the unsaturated soil model were found to be 0.999, 1.001 and 1.002 from the Janbu Method, Ordinary Method, Bishop Method and Sarma Method respectively. These Factor of Safety (FOS) were obtained by using back analysis of a two layered model. Layer 1 with soil parameters of 18 kN/m3 of unit weight, 90 kPa of cohesion and 20° of angle of friction were assigned. Correspondingly for layer 2, the values of 15 kN/m3 unit weight above the groundwater table and 18 kN/m3 below groundwater table, 3 kPa of cohesion and 8° of angle of friction were assigned. In this study, the Sarma Method was found to be the best method as it satisfies both force and moment equilibrium and is recommended for further analysis to be taken at the slope.


INTRODUCTION
Slope failure has been acknowledged as one of the most persistent natural catastrophe that can lead to immense loss in property and life.Hence, the analysis of slope stability is very important in order to protect the slopes from failing and minimizing the likelihood of slope failure.
The purpose of this study is to evaluate the slope at FKAB from the Limit Equilibrium method, using Slope/w to evaluate the critical soil parameters obtain through back analysis and eventually to obtain the best method to be use for further analysis of the slope stability.This study has been carried out at one of the slope failure at Faculty of Engineering and Built Environment (FKAB), Universiti Kebangsaan Malaysia located in Bangi, Selangor, Malaysia which had occurred along the slope behind the new faculty building.THEORY Limit Equilibirum (LE) methods are based on certain assumptions for the interslice normal (E) and interslice shear (T) forces, and the basic difference among the methods is how these forces are determined or assumed.
In a typical slice of a slope stability analysis, the forces that act on the typical slice are shown in Figure 1.W is the effective weight of the slice.The forces N' and S' are the normal and tangential components of the reaction W, respectively.E 1 and E 2 are the normal forces that act on the sides of the slice.Similarly, the shearing forces that act on the sides of slice are T 1 and T 2 .
All limit equilibrium methods involves the Mohr-Coulomb expression to determine the shear strength (τ f ) along the sliding surface.The shear stress at which a soil fails in shear is defined as the shear strength of the soil.According to Janbu (1954), a state of limit equilibrium exists when the generalize shear stress (τ) is expressed as a fraction of the shear strength.At the moment of failure, the shear strength is fully generalized along the failure surface when the critical state conditions are reached.FIGURE 1. Forces on typical slice slip surface of arbitrary shape Several limit equilibrium (LE) methods have been developed for slope stability analyses.Fellenius (1936) introduced the first method, referred to as the Ordinary or the Swedish method, for a circular slip surface.Bishop (1955) advanced the first method introducing a new relationship for the base normal force.The equation for the FOS hence became non-linear.At the same time, Janbu developed a simplified method for non-circular failure surfaces, dividing a potential sliding mass into several vertical slices.The generalized procedure of slices (GPS) was developed at the same time as a further development of the simplified method.
Later, Morgenstern-Price (1965), Spencer (1967), Sarma (1973) and several others made further contributions with different assumptions for the interslice forces.A procedure of General Limit Equilibrium (GLE) was developed as an extension of the Spencer (SM) and Morgenstern-Price methods (M-PM), satisfying both moment and force equilibrium conditions (Abramson 2002).These developments are reviewed in the following section, which aims to find out the key differences in the various approaches for FOS determination.A summary of selected LE methods and their assumptions are presented in Tab.Slope failure has been acknowledged as one of the most persistent natural catastrophe that can lead to immense loss in property and life.Hence, the analysis of slope stability is very important in order to protect the slopes from failing and minimizing the likelihood of slope failure.
The purpose of this study is to evaluate the slope at FKAB from the Limit Equilibrium method, using Slope/w to evaluate the critical soil parameters obtain through back analysis and eventually to obtain the best method to be use for further analysis of the slope stability.This study has been carried out at one of the slope failure at Faculty of Engineering and Built Environment (FKAB), Universiti Kebangsaan Malaysia located in Bangi, Selangor, Malaysia which had occurred along the slope behind the new faculty building.THEORY Limit Equilibirum (LE) methods are based on certain assumptions for the interslice normal (E) and interslice shear (T) forces, and the basic difference among the methods is how these forces are determined or assumed.
In a typical slice of a slope stability analysis, the forces that act on the typical slice are shown in Figure 1.W is the effective weight of the slice.The forces N' and S' are the normal and tangential components of the reaction W, respectively.E 1 and E 2 are the normal forces that act on the sides of the slice.Similarly, the shearing forces that act on the sides of slice are T 1 and T 2 .FIGURE 1. Forces on typical slice slip surface of arbitrary shape All limit equilibrium methods involves the Mohr-Coulomb expression to determine the shear strength (τ f ) along the sliding surface.The shear stress at which a soil fails in shear is defined as the shear strength of the soil.According to Janbu (1954), a state of limit equilibrium exists when the generalize shear stress (τ) is expressed as a fraction of the shear strength.At the moment of failure, the shear strength is fully generalized along the failure surface when the critical state conditions are reached.
Several limit equilibrium (LE) methods have been developed for slope stability analyses.Fellenius (1936) introduced the first method, referred to as the Ordinary or the Swedish method, for a circular slip surface.Bishop (1955)  where, the interslice slope, tan θ = λ.f(x) is assigned to determine the interslice forces (Krahn 2004, Abramson 2002).This figure is good to compare the most common methods in a FOS versus λ diagram as shown in Figure 2. The most likely inclinations of force equilibrium FOS (F f ) and moment equilibrium FOS (F m ) has been indicated particularly for circular SS analysis.When λ = 0, the FOS is obtained for BSM and JSM, as indicated in Figure 2. Similarly, the intersection point gives the FOS for SM or M-PM.According to Fredlund and Krahn (1977), Janbu's corrected method (JCM) and Janbu's generalised methods (GPS) are close to the intersecting point (see Fig. 2).
Interslice forces depend on a number of factors, including stress-strain and deformation characteristics of the materials.Their evaluation, however, becomes complicated in the LE methods.Therefore, simplified assumptions are made in most methods either to neglect both or to one of them.Nevertheless, the most advanced methods consider these forces in the LE analyses.
The Ordinary Method (OM) is limited to hand calculations and demonstration purposes only, whereas Bishop Simplified (BS) and Janbu Simplified Method (JSM) have been widely used for stability analyses in many years.These methods are common because the Factor of Safety (FOS) in most cases can be calculated with adequate accuracy.However, these methods have limitations in satisfying both force and moment equilibrium.
Interslice force function is defined as the ratio of tangential to normal forces acting on a slice.This function varies according to the assumptions mode.The fundamental difference in the assumptions has been summarized in Table 2. Morgenstern-Price Method assumes an arbitrary function (f(x)) with a scale factor (λ). Similarly, Janbu's GPS method relates the function to the inclination of the line of thrust.Likewise, Spencer assumed a constant function thought out the sliding surface.However, Sarma (1973) relates the relationship similar to the Mohr-Coulomb expression for shear strength.Moreover, the function in Lowe-Karafiath (1960) method is assumed as a tangent to the average inclination of the slope and sliding plane.However, the Corps of Engineers method considers the function as the average angle of the entry and exit points of the sliding surface.
The strength parameters c and φ can be total strength parameters or effective strength parameters.Slope/w makes no distinction between these two sets of parameters.Which set is appropriate for a particular analysis is project-specific, and is something we as the software user, need to decide.The software cannot do this for us.From a slope stability analysis point of view, effective strength parameters give the most realistic solution, particularly with respect to the position of the critical slip surface.
An analysis of slope stability begins with the hypothesis that the stability of a slope is the result of downward or motivating forces (i.e., gravitational) and resisting (or upward) forces.The resisting forces must be greater than the motivating forces in order for a slope to be stable.The relative stability of a slope (or how stable it is at any given time) is typically conveyed by geotechnical engineers through a factor of safety, FS defined as equation [1]; Equation [1] states that the factor of safety (FS) is the ratio between the forces/moments resisting (R) movement and the forces/moments motivating (M) movement.Sarma developed a method for a non-vertical slice or for general blocks.This method satisfies both equilibrium condition.In addition, the interslice force relationship is assumed as a linear Mohr-Coulomb expression.The interslice forces are adjusted until the Factor of Safety for force and moment equilibrium is satisfied.In summary, Sarma's method considers both interslice normal and shear forces and it satisfies both moment and force equilibrium.
Many practical problems involve assessing the shear strength of unsaturated soils.Fredlund and Morgenstern (1977) showed that the shear strength of unsaturated soils can be described by any two of stress state variables, namely, (σ n -u a ) and (u a -u w ), where u a is the pore-air pressure.Fredlund (1987) proposed the following equation [2] for the strength of unsaturated soil; Lowe-Karafiath T = E tanθ, where θ = ½(α + β) Coprs of Engrs.T = E tanθ, where θ = ½(α 1 + α 2 ) Where: ϕ b = angle indicating the rate of increase in shear strength relative to a change in matric suction (u au w ), when using (σ n -u a ) and (u a -u w ) as the two state of variables, and ϕ ' = angle indicating the rate of increase in shear strength with respect to net normal stress, (σ n -u w ) when using (σ n -u w ) and (u a -u w ) as the two state variables.

METHODOLOGY
The process in the Slope/w Analysis is briefly explained and summarized as shown in Figure 3.
A limit equilibrium analysis was carried out using the Slope/w software for the slope stability of the natural slope.The geometry was created in .gszformat and imported into the Slope/w program.The analysis type is then selected and Figure 4 shows the circular slip surface with the overall minimum factor of safety of 1.002 from the Sarma Method in Slope/w.Since the minimum value is inside the Grid, this is often use as a it is assigned with that failure of the slip circle to follow from left to right path.The Sarma analysis and half-sine function were selected and the software also gives the result of the factor of safety for the Ordinary, Bishop and Janbu method of analysis.After the material input and pore pressure were assigned, the probable slip surface for the circular failure using the grid and radius method was selected for the analysis.When the slip surface has been specified, then Slope/w runs several checks to verify the input data using the verify/optimize data command in the Tools menu.
When the verification is completed and there are no errors, then Slope/w computes the factor of safety using the method of slice selected.The minimum factor of safety is obtained for that particular analysis and its corresponding critical slip surface is displayed.
One of the problems with this method is how to locate the critical slip surface.This difficulty is solved by attempting multiple trials until we have obtained a satisfied value of overall minimum Factor of Safety of about unity.Figure 4 shows the circular slip surface with the overall minimum factor of safety of 1.002 from the Sarma Method in Slope/w.Since the minimum value is inside the Grid, this is often use as a guide to indicate that the minimum factor of safety has been found.
In order to identify the critical slip surface, the smallest Factor of Safety is to be taken from several try and error from the analysis.Table 3 summarizes all the FOS that have been obtained with 1.001 for Ordinary method, 1.002 for Bishop method, 0.999 for Janbu method and 1.002 for Sarma method.If each method gives very nearly the same critical slip surface, then this is a good indication that the true critical circular slip surface has been located.The  From the preliminary report with the borelog profile given of the slope based on Figure 5, evaluation has been made to produce the Volumetric Water Content Versus Matric Suction Curve applied in the slope stability model shown in Figure 6.
parameters of 18 kN/m³ of unit weight, 90 kPa of cohesion and 20° of angle of friction were assigned.Correspondingly for layer 2, the values of 15 kN/m³ unit weight above the groundwater table and 18 kN/m³ below groundwater table, 3 kPa of cohesion and 8° of angle of friction were assigned.TABLE 5. Comparison of residual soil strength parameters at various sites From the preliminary report with the borelog profile given of the slope based on Figure 5, evaluation has been made to produce the Volumetric Water Content Versus Matric Suction Curve applied in the slope stability model shown in Figure 6.
From the result of the slope stability analysis of Slice 1 as shown in Figure 7, the base shear mobile force to be 0.43989 kPa that has been obtained in comparison to the value of 0.52633 kPa as computed using conventional formula of equation [2].The difference that occurred in the study can be attributed to the many combinations of inputs that were used when running the software.Nevertheless, the difference obtained is relatively small.
Figure2incorporates all the assumptions and development made by the latest LE methods.In fact, this is an extension of Spencer and Morgenstern-Price methods where, the interslice slope, tan θ = λ.f(x) is assigned to determine the interslice forces(Krahn 2004, Abramson 2002).This figure is good to compare the most common methods in a FOS versus λ diagram as shown in Figure2.The most likely inclinations of force equilibrium FOS (F f ) and moment equilibrium FOS (F m ) has been indicated particularly for circular SS analysis.When λ = 0, the FOS is obtained for BSM and JSM, as indicated in Figure2.Similarly, the intersection point gives the FOS for SM or M-PM.According toFredlund and Krahn (1977), Janbu's corrected method (JCM) and Janbu's generalised methods (GPS) are close to the intersecting point (see Fig.2).Interslice forces depend on a number of factors, including stress-strain and deformation characteristics of the materials.Their evaluation, however, becomes complicated in the LE methods.Therefore, simplified assumptions are made in most methods either to neglect both or to one of them.Nevertheless, the most advanced methods consider these forces in the LE analyses.The Ordinary Method (OM) is limited to hand calculations and demonstration purposes only, whereas Bishop Simplified (BS) and Janbu Simplified Method (JSM) have been widely used for stability analyses in many years.These methods are common because the Factor of Safety (FOS) in most cases can be calculated with adequate accuracy.However, these methods have limitations in satisfying both force and moment equilibrium.Interslice force function is defined as the ratio of tangential to normal forces acting on a slice.This function varies according to the assumptions mode.The fundamental difference in the assumptions has been summarized in Table2.Morgenstern-Price Method assumes an arbitrary function (f(x)) with a scale factor (λ). Similarly, Janbu's GPS method relates the function to the inclination of the line of thrust.Likewise, Spencer assumed a constant function thought out the sliding surface.

FIGURE 2 .
FIGURE 2. Effects on FOS on force and moment equilibrium in relation to lamda (λ)(Fredlund and Krahn 1977)

TABLE 2 .
(Aryal 2006)orces and relationships(Aryal 2006) back analysis parameters used to calculate the critical slip surface for the soil layers above and below the water table in terms of unit weight, angle of friction and cohesion are summarized for the purpose of evaluation as shown in Table4.Table5compares the result of various studies from different sites on residual soil with the results obtained from this study using back analysis of a two layered model.Layer 1 with soil parameters of 18 kN/m³ of unit weight, 90 kPa of cohesion and 20° of angle of friction were assigned.Correspondingly for layer 2, the values of 15 kN/m³ unit weight above the groundwater table and 18 kN/m³ below groundwater table, 3 kPa of cohesion and 8° of angle of friction were assigned.

TABLE 4 .
Critical Soil Parameters used from back analysis

TABLE 5 .
Comparison of residual soil strength parameters at various sites