Abstract
Considering the disadvantages of the slice method commonly employed in reliability analysis of slopes, a novel method (Spoke model) was proposed for reliability analysis and safety factor calculation of slopes in this work based on geometrical relationship among slices. The safety factor and the coefficients of limit state function of slopes could be achieved with the Gaussian integral method. The minimum safety factor and the minimum reliability index, as well as their corresponding coordinates on the slip surface, can be calculated with the improved JC method and the searching method. A novel and practical method for reliability analysis of slopes has been achieved. With this method the slice process could be avoided, which helps to eliminate some calculation errors caused by oversimplified assumption. Moreover, the explicit expression of safety factor in this method requires no repeated iterative solution, which is employed in traditional slice methods, as well as can be developed into a limit state function required by calculation of the reliability index. It is demonstrated that this method works efficiently and succinctly in evaluation of reliability index and safety factor for soil slopes.
1. Introduction
Each slope (including stable and unstable slopes) possesses an inherent tendency to degrade to a more stable form and eventually make the slope tend to horizontal. Mass movement (slope rupture, landslide, etc.) may take place as a result of the shear failure. The Spoke model proposed in this work is mainly focused on the stability analysis and reliability analysis of simple slopes, where the concept of simple slope discussed in this work is of the slope morphology as shown in Figure 1, and the concept includes the following aspects: the soil is homogeneous, and considering its geotechnical parameters it only includes weight density “”, cohesion “”, and internal friction angle “”; the underground water are not taken into consideration.
2. Literature Review
There are many methods to determine a slope slip surface. The center coordinate and the radius of a sliding circle were employed to determine the slip surface in one method [1, 2], while the center coordinate and the height of a sliding circle were introduced in another method [3]. A variety of methods have been applied to analyze the reliability of slopes. Based on traditional slice methods (such as Bishop method [3], Spencer model, Sarma model, and Janbu model [4, 5]), the reliability index can be achieved with the calculated safety factor and the limit state function . In other cases the minimum reliability index was obtained with the limit state function on the slip surface where the safety factor reaches its minimum value [6]. These approaches are approximate methods for calculation of the reliability index, because all the expressions of safety factor in the slide methods are implicit, which require complicated iterative process when solving. Moreover, the slip surface corresponding to the minimum reliability index may not be the slip surface corresponding to the minimum safety factor [7].
The parallel slide method was employed to analyze the stability of slopes in Bishop method. The safety factor is expressed as follows:
Obviously, the equation is implicit where the parameter contains as well; therefore, it can be solved only with the iterative method. Since the expression is implicit, the corresponding limit state function cannot be achieved when analyzing the stability of slopes with traditional slice methods.
3. Spoke Model
To achieve the limit state equation, a novel slide model (the Spoke model) is proposed in this work. With this model we can get the limit state equation from which the reliability index of slopes can be obtained.
Assuming that a slope is homogeneous, the height, slope angle, slope weight, angle of internal friction, and cohesion are , , , , and , respectively. The slip surface is supposed to be a circular arc surface (Figure 2). Point is the center of the slope line. Points , , and are on the same line. The angle between line and the slope line is . Since and , the arc is the slip line. Translate left the center , the arc , segments, and intersections with a distance of . The coordinates correspond to all possible slip surfaces, respectively.
The slip soil can be divided into three parts , , and . Taking the second part , for example, the center angle has crossing points AFDG with the slip line and slope line. The area of the th AFDG is . The radius of slip circle is which can be expressed in Considering is a small value, Assume that In triangle , according to Sine theorem, it can be achieved that Combining (5) and (7), Similarly, in triangle , Combining (6), (8), (9), and (4), Combining (3) and (10),
Assuming that the area of the th slice is , bulk density and weight of the soil are and , respectively, as follows:
In the second part, is determined by (11).
Similarly, the expression of the first and third parts can be achieved, respectively, as follows:
Shear forces and pass through center , perpendicular to the slip surface. The sideward normal forces are and , respectively. According to the balance condition along the direction of each radial line,
Since is a definitely small value,
Combining (16), (17), and (18),
In the limit equilibrium, the sum of all the force moments on the center should be zero; therefore, where and are the inertia force of earthquake and its tension arm, respectively. is the distance between center of the circle and the weight center of a slice. The shear stress at the bottom of the slice can be expressed as
Combining (19), (20) and (21), the safety factor can be achieved. where the sum of all the force moments of the slide friction forces on the center of the circle resulting from shear stresses should be zero. Therefore,
Compared with the expression in Bishop method, (23) possesses three obvious advantages as follows: (i) the shear stress between the soil strips is not involved in (23); (ii) it is an absolute expression without in its right part, which requires no repeated iterative solution; (iii) (23) has a clear physical meaning. is the ratio of the total resisting moment to the total sliding moment.
When the slope is sliced infinitely, as well as the interstitial pressure and the seismic inertial force are not taken into consideration, the safety factor can be expressed in the integration form. Combining (13), (16), and (23), where .
Therefore, where () in (24) and (25) is the integrated area element, () is the distance between the slice center and the center of the circle, and () represents the quadrature. Combining (11)–(15), they could be expressed as the functions of the parameters , , , , and as follows: where
In (24) and (25), the upper/lower limit of integral represents the central angle between the line and another segment whose two endpoints are the center of circle and the endpoint of the slide slope segment in the circle arc, which can be expressed as the functions of the parameters , , , , and as follows:
In (24) and (25), can be achieved by solving differential equations. The differential equation in the tangential direction of the th slice slope can be expressed as follows:
To simplify the calculation, let , , and (32) can be developed into a differential equation: where is the area element, which can be expressed as follows with (26): Let . The general solution of the linear differential (33) is where is an arbitrary constant. Considering the complexity of , it is definitely complicated to solve the indefinite integral of the transcendental function. To simplify the calculation, is replaced with a proper constant . Within the interval in (26), is replaced with . Equation (32) can be developed into
With (36) can be developed into Assuming when ,
When , we can obtain the maximum value of broadside normal stress and its corresponding angle . The normal stress would reduce to zero gradually from the cut-off point toward both directions. It can be analyzed that is in the interval . Combining and (24), (25), the safety factor can be calculated, and the limit state function can be determined.
Subsequently, combine (24)–(30). As for a certain slope (i.e., and are fixed.), each group of would result in its corresponding safety factor and limit state function. Namely, safety factor and limit state function can be expressed as a function containing variables , , and , respectively as follows:
Considering the random variability of soil parameters, we could obtain the reliability index of each group of . By varying the numerical values of , the minimum values of reliability index and safety factor can both be obtained. The calculation of structure reliability index has been discussed in many reported works [9–13]. In a typical failure mode of a slope (Figure 3), the slip surface passes through both the top plane and bottom plane. Based on geometrical condition and the constraint relationship among , , and , inequalities (41) can be obtained theoretically as follows:
Pan Jia-zheng has stated that the center of critical slip surface is in the area formed by two circulars with radius of and , respectively, central normal line and perpendicular bisector of the slope surface (Figure 3) [14]. By varying the values of , , and within this area, we can achieve the minimum values of reliability index , and safety factor , and their corresponding .
3.1. Example Analysis
Example (I) (see [15]). There is a slope whose height and slope angle are 4.5 m and 20°, respectively, with a sand bed beneath which helps to drain water away adequately. The natural unit weight of soil is kNm−3. The average efficacious shear strength and variation coefficient are kPa, , , and , respectively. Find out the failure probability of the given slid surface. The safety factor was calculated approximately to be with a vertical slice method in [15]. The failure probability was calculated to be () with statistics formulas.
With the novel method (Spoke model) proposed in this work, the safety factor is calculated to be . If and are regarded as normal random variables, and they are not interrelated, the corresponding reliability index is calculated to be 2.883 with the Spoke model.
Example (II) (see [8]). There is a homogeneous simple slope whose height () is 50 m, slope proportion () is 3.25, bulk density of soil () is 19.62 kNm−3, coefficient of friction is , and adhesion stress () is 58.66 kPa. Work out the location of the most dangerous slide arc and minimum value of safety factor.
The calculated results with the Spoke model proposed in this work and the corresponding results in [8] are listed in Table 1.
Example (III) (see [7]). There is a clay slope whose height is 10.0 m, and the proportion of slope is 1 : 2. The average value and standard deviation of and are 20 kPa, 0.1817 and 4.0 kPa, 0.05, respectively, and all of them are considered to be of normal distribution and not interrelated. The unit weight of soil is 16.0 kNm−3, considered as a fixed value. Find out the minimum value of safety factor and reliability index of the slope.
With the Spoke model in this work, we take °, , and as the step length to search solution. The results are listed in Table 2.
According to the analysis of Examples (I) and (II), the calculated results with Spoke model and traditional slice methods are subequal. Example (III) indicates that minimum safety factor and minimum reliability index are not achieved on the same slip surface. In addition, the variability of soil parameters affects the minimum reliability index and the position of its corresponding slip surface but possesses little influence on the minimum safety factor and the position of its corresponding slip surface.
It is worth noting that the expression of safety factor has a clear physical meaning in this work. With the calculation program built, it is convenient for data input, and the calculation process is quick. With this method for reliability analysis of slopes, the limit state function can be achieved which directly results in the reliability index.
4. Conclusion
The explicit expression of safety factor in Spoke model proposed in this work requires no repeated iterative solution, as well as can be developed into a limit state function required by calculation of the reliability index. The Spoke model avoids complicated repeated iterative solution as required in traditional slice methods as well as eliminates the errors arising from oversimplified assumption.
According to the example analysis, the slip surface with minimum reliability index is not in coincidence with, but similar to, the one with minimum safety factor. We can find out the minimum safety factor and its corresponding coordinates first, then search the minimum reliability index nearby.
The minimum reliability index and location of its corresponding slip surface are related with soil parameters such as weight density “,” cohesion “,” and internal friction angle “.” The minimum safety factor and location of its corresponding slip surface have no relation with variability of soil parameters.
This work only discusses the calculation analysis of reliability of simple slopes. The failure surface is supposed as a circular arc, and some simplifications are utilized on the resistance calculation of slip surface caused by the stress between slices. The influence of simplified treatment in this work on calculation results could be a direction of further research.