The effect of slope height and angle on the safety factor and modes of failure of 3D slopes analysis using limit equilibrium method

It is a well-known fact that the safety of slopes majorly depends on several factors such as geometry, soil properties. The safety factor might change significantly depending on the soil type and the slope’s shape. The knowledge of the nature of the change in the safety factor due to the change in the slope’s height and angle is essential for implementing an effective strategy of increasing the safety factor for any slope stability problem. The influence of geometric shapes on the stability of the slope needs to be properly investigated through three-dimensional slope stability analysis, as the three-dimensional analysis is suitable for all slopes, even those which invalidate the plane-strain conditions. To calculate the three-dimensional safety factor, multiple analyses of three homogenous soil slopes with different soil properties were conducted by varying slope height, angle, and combinations. Each slope's height and angle were recorded to identify the types of slope failure. The analysis’s findings showed that while a decrease in height raises the safety factor nonlinearly, a decrease in slope angle increases the safety factor almost linearly. Base failure is the most likely failure for slopes with a height less than 4.0 m and an angle of inclination less than 18°. On clay and sandy clay soils, toe slide is the most common type of slope failure. The expected failure type will be either toe or face failure when the slope's height and base angle exceeds 5.0 m and 22°, respectively. This study also found that the three-dimensional safety factor for soil slope is generally 10–20% higher than the two-dimensional factor of slope safety. The slope’s nature depended on the soil type and slope form, but the safety factor increased as the slope angle and height decreased. To determine the most efficient method for slope stabilization, it is necessary to do an extensive study on slope height and angle reduction techniques. It should be ensured that the sliding mass of soil does not rise, resulting in a potential slope failure. The present study will help identify the correlation between the height and base inclination of the slope with the expected nature of slope failure. The present study helps to investigate the variation of the safety factor of a three-dimensional homogenous soil slope subjected to self-weight only. The study can be further extended to observe the variation of the factor of safety for a 3D slope subjected to pore water pressure and seismic loading also.

financial losses.Roads, embankment cuts, fills, and dams must be investigated with due care to ensure the project's safety.Slopes should be adequately inspected and fixed so that they do not fall apart in a way that causes a disaster.Failures of slopes rely on soil type, geometry, soil stratification, groundwater, and infiltration [2,3].Failure of a slope can be caused by translational slide, rotational slide, or flow [4,5].Fine-grained and homogenous soils are prone to rotational slides, as shown in Fig. 1.In the analysis of soil slopes, the limit equilibrium method (LEM), limit analyses (LA), and finite element method (FEM) are the most often used methods [6].In slope analysis, limit equilibrium approaches are widely used to establish the slope factor of safety (FS) against failure [7][8][9][10][11][12].Initially, the stability checks of a slope were formulated in two dimensions (2D) under the assumption that plane-strain conditions prevailed.But the plane-strain assumption is often wrong when the section changes along the slope's longitudinal direction.Singh et al. [13] conducted a comprehensive analysis of conventional and soft computing techniques for slope stability analysis.In these cases, a three-dimensional (3D) slope stability study is needed to determine how the slope will fail [14].A few scholars have developed a limit equilibrium theory-based 3D slope stability study [15][16][17][18][19].The collapsing mass is symmetrical and divided into vertical columns.Researchers reported 3D method-of-columns techniques based on a symmetrical plane for failure mass [16,[20][21][22][23][24][25].Recently Rao et al. [18] studied the Box Search technique for 3D Slope Stability Analysis Using the LEM based on the symmetric plane.Some formulations based on limit equilibrium techniques did not use a symmetrical plane [26][27][28].Since [29], other scholars [30][31][32] have used LEM to solve the FS of asymmetrical slopes [33,34].presented a comprehensive literature study involving slope stability analysis.
With limit analysis, slope analysis is looked at regarding energy balance, and the results are accurate.The upper bound LA looks for the slope failure mechanism by making an acceptable velocity field from a kinetic point of view [35].The upper bound theorem of LA is used to evaluate slope stability because no assumptions regarding interaction forces and specified failure surfaces are necessary [36][37][38][39].Limit finite element analysis (LFEA) was also utilized to investigate the topic of slope stability.Numerous studies have accounted for slopes' complex geometry and constitutive links [40,41].The lower bound theorem is attractive because it provides a reliable estimate of the load capacity of a structure by assuming a stiff plastic material model with an associated flow rule.Most numerical implementations of the lower bound theorem are based on discretizing the continuum using finite elements.Using the LA method's upper bound, many researchers have done 3D slope stability analyses [42,43].The displacement finite element method (FEM) and finite difference analysis (FDA), which are preferable for deformation investigations, have also been applied to slope stability evaluation based on the strength reduction method [44,45].
The search procedure given by [29] for locating the critical failure surface (CFS) and corresponding minimum FS is compatible with other well-known column techniques, as demonstrated by [22].Previously, the 3D simplified Janbu approach [46] was utilized to develop an efficient search strategy for locating important slip surfaces.Horizontal and vertical force equilibrium conditions were satisfied by Ugai's 3D Janbu method.Consequently, it was discovered that critical slip surfaces reported using the 3D Janbu method were significantly deeper than those found using variational formulations [47] in a number of situations, particularly for slopes in cohesive soil with low friction.In addition, the 3D simplified Janbu technique produced minimum safety factor values that were more conservative than variational solutions.
Sustainability concepts are in Policy 418 of the American Society of Civil Engineers."Do the Right Project; Do the Project Right; Perform Life Cycle Assessment from Planning to Reuse; Use Resources Wisely; Plan for Resiliency and Validate Application of Principles" are the policies for civil engineers [48].This motto highlights the need for project safety analysis.Adding support structures and modifying slope geometry can increase slope stability.This research aims to study the height, angle parameters and mode of failure surface for a 3D slope using LEM.Slope safety can be increased by decreasing slope height/angle; however, the safety depends on the soil type.The strategy of reducing the slope angle is more efficient for some soils than decreasing slope height, This research aims to study the height, angle parameters, and mode of failure surface for a 3D slope using LEM.
The safety factor might change significantly depending on the soil type and the slope's shape.The knowledge of the nature of the change in the safety factor due to the change in the slope's height and angle is essential for implementing an effective strategy of increasing the safety factor for any slope stability problem.

Design
This study uses 3D slope stability analyses using the limit equilibrium method.For 3D slope stability, investigations were carried out using the Scoops3D source program.Scoops3D allows Ordinary and Bishop simplified methods to compute the FS, but in this research work, Bishop's simplified method is used for chosen problem.

Setting of the study
This study investigates the critical failure surface and safety factor for different heights and base inclination angle combinations for three types of soils: clayey, sandy clay, and sandy soil.To produce the 3D geometric profile of the slope, a digital elevation model (DEM) input file is prepared.The horizontal resolution of the DEM grid (m) is taken at 0.50 m intervals.An initial height variation is performed in this study to determine the mode of failure surface and safety factor.The second analytical approach involves gradually changing the slope angle while maintaining the slope height and soil characteristics constant.
In a third analysis, the author investigated different slope height and angle combinations while keeping the soil parameters the same.

Geometric design of 3D soil slope
Scoops3D is a software from USGS (United States Geological Society) that can perform three-dimensional slope stability analysis.Scoops3D analyses slope stability by extending the traditional 2D limit equilibrium formulation to three dimensions.It computes the stability of a stiff soil mass covered by the spherical trial surfaces (potential sliding surfaces).For three-dimensional analysis, the entire slope domain is divided into several columns, as shown in Fig. 2. The partial columns for which at least two sides fall inside the failure surface are considered while estimating the failure mass.If partial columns are considered while calculating the failure mass, the final output results' accuracy is better [49].The column width, also known as DEM cell size, is given as input by the user.The DEM data includes the slope surface's top elevation (i.e., z coordinate).The geometry data of any intermediate layer inside the slope are also specified in the DEM file.The average of the four surrounding DEM cell heights is used to calculate column corner elevations.The piezometric surface profile, if any, also needs to be defined using a DEM input file.

3D formulation of bishop simplified analysis
In this research, the FS of a 3D slope is determined using a Scoops3D-based computer program.Scoops3D utilizes Ordinary and Bishop simplified methods (BSM) to compute FS, but in this research work, Bishop's simplified method is used for the chosen problem. Figure 3 represents the free body diagram of the j,k column as no external force was acting on the column subjected to all possible combinations of forces.
Where W is the weight of the column; E x j,k , E y j,k = x and y directions inter-column normal force;H x j,k , H y j,k = Hori- zontal shear force in y-z plane; X x j,k , X y j,k = Inter-column shear force in x-z plane;N j,k , U j,k = effective normal force and base pore water force; S j,k = Mobilized shear force acting on the base; α j,k = Slide angle relative to the x-y plane; α x , α y = base inclination in x-z and y-z planes at the middle of each column.
The final safety factor expression for Bishop's simplified method with no groundwater condition is expressed in Eq. (1).Here, m α j,k = cos ε j,k + tan ϕ ′ d m z , and m z = sin α j,k ,c j,k = the effective cohesion;ϕ j,k = the effective internal friction angle; R j,k is the distance from the j,k column's trial slip region to its axis of rotation; A j,k is the trial surface area of the column; W j,k is the weight of the column; and ε j,k is the angle between the inclined surface at the bottom of the slice and the horizontal x-axis.For 3D formulation, the summation of normal and shear forces acting along the sides of the columns are assumed to be equal to zero for both the x and y directions.

Search for critical failure surface and associated FS of 3D slope profile
In Scoops3D, the three-dimensional slope profile is created using the DEM technique.The user inputs the column width, also known as the DEM cell size.DEM input files contain DEM cell surface elevation data.In 3D slope stability analyses, a search method known as Box Search Method is used to find out the critical failure surface and corresponding minimum FS. Figure 4 shows a 3D search lattice of a DEM profile.During the search phase, Scoops3D keeps track of the DEM cells for which the smallest safety factor was determined among all trial surfaces containing that cell.A sphere with a rotatable center point above the DEM and a predetermined radius must contain each trial surface to facilitate the search process.The trial surface with the lowest safety factor is called the critical surface for that DEM cell. (

Correlation coefficient of slope height and slope angle with safety factor (FS)
The term correlation coefficient (r) can be defined as the degree of relationship between two variables.Pearson established the Pearson correlation, often known as the correlation coefficient, based on the work of others, including Galton [50], who first found the concept of correlation [51].The value of the r is always between -1 and + 1. Pearson's correlation coefficient (r) is the covari- ance ratio of two variables to the product of their standard deviations.A positive correlation coefficient indicates if one variable increases, the other value also increases, while a negative correlation coefficient indicates if one variable increases, the other value decreases.A zero-correlation coefficient indicates no relationship between the two variables.The expression of (r) is expressed below: where cov is the covariance, σ x is the standard deviation of variable x, and σ y is the standard deviation of variable y.
The correlation coefficient can be determined using the formula if x and y are the variables under consideration.
where n is the number of variables.

Results
Several factors, including soil type and slope geometry, determine the various modes of slope failure.This study investigates the critical failure surface and safety factor for different heights and base inclination angle combinations for three types of soils: clayey, sandy clay, and sandy soil.The soil properties of these soils are shown in Table 1, taken by Shiferaw [3].To produce the 3D geometric profile of the slope, a DEM input file is prepared.2 and 3 summarize slope failure modes as slope height and angle change.

3D Safety factor variation with a slope height
An initial height variation is performed in this study to determine the mode of failure surface and safety factor.The height variation is taken from 12.0 m to 1.0 m at a 1.0 m interval while other parameters have been kept constant.For the three types of soil, the safety factor increases with a decrease in the height of the slope, as shown in Table 4. Shiferaw [3] performed two-dimensional slope stability analyses using PLAXIS-2D for three soil types: clayey, sandy clay, and sandy soils.The 3D slope stability analyses performed for a similar problem reveal that the 3D FS values are 10-20% higher than their 2D values.The differences between 2 and 3D are frequently in the range of 10% to 20% [21,25,33].In this study, the correlation coefficient (r) is carried out Unit weight (kN/m 3 ) 18 16 18 between the height of the slope and the 3D safety factor for all three types of soil, and a negative value is obtained, which indicates a good correlation between the height of the slope and the safety factor.

3D Safety factor variation with slope angle
The second analytical approach involves gradually changing the slope angle while maintaining the slope height and soil characteristics constant.For the three types of soil, the safety factor increases with a decrease in the slope angle, as shown in Table 5.In this case, a strong opposite correlation coefficient was found between the slope angle and 3D safety factor for three soil types.

3D Safety factor variation with slope height and angle
In a third analysis, the author investigated different slope height and angle combinations while keeping the soil parameters the same.This analysis finds slope heights of 12 m, 9 m, 6 m, and 3 m.The slope angle considered for each slope height is 30.96• ,25.64 • , 21.8 • , and 18.77 • , respectively.In this analysis, sandy clay soil is the only type of soil considered for analysis.The calculated safety factor value for the third analysis is shown in Table 6.The reported 3D safety factor values for sandy clay with slope and height variation are compared to the 2D safety factor [3]. So, the increase in safety factor values is observed at 10-20% compared to 2D slope analysis.

Discussion
Shiferaw [3] performed two-dimensional slope stability analyses using PLAXIS-2D for three soil types: clayey, sandy clay, and sandy soils.The 3D slope stability analyses performed for a similar problem reveal that the 3D FS values are 10-20% higher than their 2D values.

The effect of slope angle and height on the failure mode of slope
Slope failure can occur in three ways: slope/face failure, toe slide, and base slide.Both the kind of soil and the slope's geometry influence the failure mode.Tables 2 and  3 show that the slopes fail differently depending on their height and angle.Toe slide is the most common slope  failure on clay and sandy soils.The most common type of failure in sandy soil is a slope/face slide.In this study, all three soil types: clayey, sandy clay, and sandy soil, fail by base sliding at heights less than 4 m.Slope/face slide failure is typical for sandy clay soils on steep slopes, while toe slide failure dominates on mild slopes.Toe sliding is the predominant mode of failure for clayey soils.The base slide happens when the slope angle is less than 18 0 .In sandy soil conditions, slope/face slides are the most typical form of failure.

The effect of slope height on safety factor
When the slope height increases in the later stage, there is a general trend toward a slower rate of rise in the safety factor.Figure 8 shows that the safety factor increases quickly for height changes of less than 3 m. Figure 8 depicts a reasonably steady slope of the line for a range of slope heights from 12 m down to 4 m.The correlation examination reveals a significant negative linear relationship between slope height and FS, with r = − 0.96159.Observe the abrupt shift in the slope of the graph between height and safety factor at 3 m.Correlation analysis between 3 and 1 m heights yields r = − 0.95437.The overall r is − 0.7540, which means the safety factor decreases slower as height goes from 12 to 3 m.After a slope is 3 m high, the safety factor increases faster.

Slope angle's effect on the safety factor
As the slope gets steeper, the safety factor drops almost in a straight line.Figure 9 demonstrates that the correlation between the slope angle and the safety factor is almost straight while the angle gets steeper.The correlation coefficient between the slope angle and the safety factor is − 0.9793.According to the correlation coefficient, there is a significant and inverse association between slope angle and FS for the three types of soil.

The effect of slope height and slope on safety factor
The safety factor can be maximized by optimizing slope height and angle.Calculating the FS shown in Table 6 allows researchers to examine the effects of decreasing slope height and angle.As shown in Fig. 10, a slope height of 3 m or less increases the FS at a higher rate than  a height of 3 m or more.So, decreasing the slope angle on slopes with less height increases the FS at a higher rate.Figure 11 illustrates the difference in the FS with respect to height at various slope angles.When slope height increases, the safety factor decreases.Further, it is seen that the safety factor increases when the slope angle decreases.
The variation of FS with slope's height and base angle has been earlier studied by Shiferaw [3].In two dimensions, in 3D, the authors could not find any similar study.Therefore, it is difficult to compare these results with previously published results.

Conclusions
The present study is useful for investigating the variation of the safety factor of a three-dimensional homogenous soil slope subjected to self-weight only.The results indicate that the slope's height and angle influence how the three soil types fail.Toe slip is the predominant slope collapse for clay and sandy clay soils.It is noticed that for sandy soil, slope failure is the dominant mode of failure.All three soil types: clayey, sandy clay, and sandy soil, fail due to base sliding at heights below 4.0 m.Sandy clay soils collapse by slope slide on steep slopes and toe slide on mild slopes.Clayey soils fail mostly by toe sliding.Slopes under 18 ° cause base slides.Most sandy soil failures are slope/face slides.Slope safety generally increases linearly as the slope angle decreases, although the safety factor rises at variable rates as the slope height decreases.Slope heights under 3.0 m raise safety factors faster than those over 3.0 m.Thus, decreasing the slope angle on lower slope height enhances the safety factor.The expected failure type will be either toe or face failure when the slope's height and base angle exceeds 5.0 m and 22 °, respectively.This study also found that the three-dimensional safety factor for soil slope is generally 10-20% higher than the two-dimensional factor of slope safety.The correlation coefficient shows a strong inverse correlation between the safety factor and slope angle.Knowing the failure mode and the influence of geometric change on the slope safety factor can help researchers choose the best method for improving slope stability.The study can be further extended to observe the variation of the FS for a 3D slope subjected to pore water pressure and seismic loading.

Fig. 1
Fig. 1 Different forms of rotating slope failure

Fig. 2
Fig. 2 Planar view of the potential sliding mass of divided columns for a 3D slope

x 2 n y 2 − y 2 Fig. 3 Fig. 4
Fig. 3 Free body diagram of j,k column

Fig. 11
Fig.11 Slope height vs. safety factor for slope angle

Table 1
Geotechnical properties of soil

Table 3
Critical failure mode based on slope angle at slope height H = 12 m

Table 4
3D Safety factor based on slope height at slope angle •Slope height (m) Safety factor

Table 5
3D Safety factor based on slope angle at slope height H = 12 m

Table 6
3D safety factor based on Height and Angle of slope