Abstract

Based on the upper bound theorem of limit analysis, this paper presents a procedure for assessment of the influence of the soil anisotropy and nonhomogeneity on the stability of fissured slopes subjected to seismic action. By means of a mathematical optimization procedure written in Matlab software codes, the stability factors NS and λ are derived with respect to the best upper bound solutions. A series of stability charts are obtained in this paper, and then the critical locations of cracks are determined for cracks of known depth. The results demonstrate a significant influence of the soil anisotropy and nonhomogeneity on the stability of the fissured slopes and the location distribution of the cracks. In addition, the procedures for getting the factor of safety are put forward. It is shown that a decrease in the nonhomogeneity coefficient n0 and an increase in the anisotropy coefficient k could lead to the fissured slopes becoming unsafe. Finally, this article also illustrates the variation in the safety factor of fissured slopes under the impact of three factors (Kh, H1/H, and λ).

1. Introduction

Cracks or fissures are often found at the crest of cohesive soil slopes because of the development of tension. The presence of cracks poses significant challenges to the assessment of slope stability because cracks introduce one discontinuity in both the static and kinematic fields [1, 2]. Taylor [3] pointed out that “……the action within the tension zone is a subject that is worthy of much study.” Therefore, stability analysis of the fissured slopes is a very important practical engineering problem. The used traditional methods to analyze the effects of cracks on the stability of slopes have been categorized as follows: finite element limit analysis (FELA) methods [4, 5] and limit equilibrium (LE) methods [610]. Baker [11] described the procedure of the limit equilibrium method with variational extremization for the evaluation of the effect of tensile cracks of the cohesive-frictional soil slope. After this, other researchers also developed the variational limit equilibrium methods to obtain the factor of safety for giving fissured slopes [12, 13]. Another approach to analyzing the stability of fissured slopes is to use the upper bound theorem of limit analysis (LA). Limit analysis provides a powerful method for analyzing stability problems [14, 15]. Based on the LA approach, Utili [1] and Utili and Abd [2] investigated the stability of homogeneous slopes manifesting vertical cracks, and the depth and locations of cracks. Michalowski [16, 17] introduced a crack-opening process within the strict framework of LA. Because the LA approach is comparatively more convenient to use than the static approach [18], adopting it to solve the stability problem has almost exclusively concentrated on the kinematic theorem [19]. Moreover, the results from the upper bound theorem are nearly identical to those from the variational limit equilibrium closed-form solutions [20]. Therefore, the LA approach is appropriate for assessing the stability of the fissured slopes.

Seismic action can trigger slope instability, especially in the case of slopes with cracks [2, 21, 22]. The conventional method for estimating the effect of a seismic force on the stability of a slope is the so-called “pseudostatic analysis method,” which was employed by Terzaghi [23]. According to this method, the inertia force is treated as an equivalent concentrated horizontal force at the center of gravity of the critical sliding mass [24]. An analytical method based on the upper bound theorem in combination with the pseudostatic method is widely presented for the stability of the fissured slopes and subjected to seismic action [2, 22, 2426]. In the previous study of fissured slope stability, soil was treated as an isotropic and homogeneous geotechnical material. However, due to the impact of inherent soil structures [27], cutting action [28], consolidation [29, 30], cementation bounds [31], and soil strength linear increase with depth [1820, 24, 3234], soil deposits can exhibit remarkably anisotropic and nonhomogeneous behavior in nature. Thus, ignoring the anisotropic and nonhomogeneous behavior can cause the geotechnical materials to be unrealistic in fissured slope stability analysis.

In this paper, the upper bound theorem of limit analysis together with the pseudostatic method is used to assess the stability of fissured slopes subjected to seismic action where the anisotropy and nonhomogeneity of the soil strength is considered. A rotational two-dimensional (2D) failure mechanism proposed by Chen et al. [19], Chen and Sawada [24], and Ukritchon and Keawsawasvong [35] is adopted. The explicit functions about the stability factors of the fissured slope and the location of the cracks are derived, which are optimized to find out the minimum stability factor and the critical crack location, respectively. Next, the procedures for getting the safety factor of fissured slopes are presented where the effects of soil anisotropy and nonhomogeneity and seismic action are both taken into account. Furthermore, the variation in the safety factor of fissured slopes under the impact of the crack depth and the horizontal and vertical seismic acceleration is also given in this article.

2. Upper Bound Theorem of Fissured Slope Stability

2.1. Basic Assumption

(1) The research problem can be studied in a plane-strain problem; (2) the Mohr–Coulomb criterion is adopted in this study; (3) when fissured slope soil is in the critical state, a failure surface passes through the toe of the slope; (4) the anisotropic and nonhomogeneous qualities of slope soil lie only in its cohesive strength; (5) the seismic action of the sliding wedge is simplified to the static load.

2.2. Anisotropy and Nonhomogeneity of Soil Strength

Figure 1 illustrates a fissured slope with soil anisotropy and nonhomogeneity and its corresponding potential failure surface. In previous studies, the slope soil was regarded as an isotropic and homogeneous material. However, plenty of experimental tests showed that the actual soil exhibits remarkable anisotropy and nonhomogeneity naturally, which should not be ignored in slope stability analysis.

Soil anisotropy indicates that the behavior of soil shear strength changes in different directions. As shown in Figures 1(a) and 1(b), the cohesion strength ci is in the direction of the maximum principal stress, which is inclined at an angle of i with the vertical direction. The cohesion strength ci formulation follows the notations presented by Lo [36] and Casagrande and Carillo [37]. A brief description of the relevant formula can be expressed aswhere ch and represent the horizontal and vertical cohesive strength, respectively. The ratio of ch to can be defined as k (k = ), which represents the anisotropy of the slope soil. As Lee and Rowe [27] discussed, k = 1 for isotropic soils, k < 1 referred to normally consolidated or slightly overconsolidated clays, and k > 1 referred to heavily overconsolidated clays.

Based on the geometry illustrated in the inset of Figure 1(a), the angle, i, can be calculated aswhere θ is any counterclockwise positive angle from the initial radius and φ is the friction angle of slope soil. The parameter δ represents the angle between the failure surface and the surface perpendicular to the direction of the principal cohesion strength. As presented by Chen et al. [19], Su et al. [38], and Pan and Dias [39], δ is taken as . In addition, δ is a constant parameter regardless of the direction of the major principal stress σ1.

For nonhomogeneous soils, nonhomogeneity is the term of soil strength and is used to describe the behavior that soil cohesion varies linearly with depth from the ground surface (Figure 1(c)). As defined by Chen et al. [19], nonhomogeneity only refers to the cohesive strength. In this study, the cohesive strength nonhomogeneous coefficient is called n0. Note that herein the coefficient parameter n0 is n0 = 1, when the soil mass shows homogeneous cohesion.

2.3. Definition and Formulation of Problem

To establish the 2D limit analysis method, Utili [1] and Michalowski [16] conducted some research on the 2D rotational mechanisms for fissured slopes considering the log-spiral failure surface. Afterwards, Utili and Abd [2] adopted the 2D failure mechanisms for toe failure to investigate the effect of seismic action on fissured slopes employing the upper bound method together with the pseudostatic approach. Hence, this study utilized the 2D upper bound method of Utili and Abd [2] to derive the crack location, stability factors, and safety factors for soil fissured slopes with seismic action.

The failure mode is shown in Figure 1(a) with the slope angle β, self-weight γ, and height H. Therein, BC represents a vertical crack, which is through the upper surface of a slope. The depth of the crack BC gives H1. The length of EF is L1 and of BF is L2. The slope is represented as a rigid body that slides along the log-spiral surface with center O and angular velocity ω. The vertical crack BC forms the sliding surface combined with the shear slip surface, namely, surface BCD.

The logarithm-spiral formula can be written aswhere the parameters r0, θ0, and θh refer to the geometrical relationship in Figure 1(a) and ψ = tan φm.

For presentation of the location of the crack, the normalized parameters L1r and L2r are adopted in this study. The functions L1r and L2r are determined bywhere the parameter θe refers to the geometry in Figure 1(a).

The rate of external work for sliding EBCD, Wr, is calculated as the work of block EFD minus the work of block BFC. Moreover, the work of blocks EFD and BFC is calculated with the same method. Thus, the external work rate of the soil weight can be written aswhere

In addition to the weight force, a horizontal pseudostatic force, KhW, and a vertical one, W, are added to account for the seismic action. The method for calculating the rate of external work due to the seismic forces is similar to that for determining the rate of external work due to the soil weight; the only difference is in the velocity direction when calculating the rate of external work due to the seismic force which is horizontal [22].

The rate of external work attributing to vertical pseudostatic force can be written aswherein . From Figure 1(a), it can be seen that the + sign indicates vertical downward acceleration, whereas the sign indicates vertical upward acceleration.

Moreover, the rate of external work attributing to horizontal pseudostatic force can be written aswhere

Note that there is no internal dissipation of energy along the crack since the crack is prior to the slope failure. Thus, the energy is just dissipated internally along the log-spiral failure surface CD. The rate of internal energy dissipation is found by multiplying the differential area r dθ/cos φ by ci times the discontinuity in , and the total internal dissipation rate of energy is calculated by

The influence of anisotropy and nonhomogeneity on the cohesive strength ci can then be obtained as follows:where

Hr is calculated by

Knowledge on the known depth of the crack, H1, brings in an additional equation about θe, which can be obtained by combining equation (13) and the geometry in Figure 1(a) so that the following constraint is found:

The parameter FS is the strength reduction factor, which is selected as the factor for evaluation of slope stability. To obtain the true factor of safety, the strength reduction factor is gradually increased until the failure of the slope [20]. When this critical value has been found, the factor of safety of the fissured slope is equal to the strength reduction factor. The shear strength parameters (c and φ) are divided by the safety factor FS, which are analytically defined as

Expressions of the rate of external work done by the soil weight , the vertical pseudostatic force , the horizontal pseudostatic force , and the rate of internal energy dissipation for the rotational mechanism can be obtained. The work-energy balance equation, which defines that the rate of external work equals the rate of energy dissipation, is as follows:

A stability number, NS = c/γHFs, is defined here as presented by Taylor [40]. From equation (16), NS is obtained by the following expression:where

The functions fd1 and fd2 appearing in the above equation are calculated asin which

For fissured slopes with soil anisotropy and nonhomogeneity subjected to seismic action, a dimensionless parameter λ = cm/γHtanφm is introduced here to calculate the safety factor FS:

Based on the mathematical optimization method, the minimum of NS = f1(θ0, θe, θh, φ, β, λ, n0, k, Kh) and λ = f2(θ0, θe, θh, φ, β, λ, n0, k, Kh) over the two geometrical variables θ0 and θh provides the best upper bound on NS and λ, respectively. In this study, a computer program is developed by using the Matlab software to implement the optimum algorithm.

3. Results and Discussion

3.1. Comparisons

In the isotropy and homogeneity condition (i.e., k = 1.0, n0 = 1.0), and in the static conditions (i.e., Kh =  = 0), the calculated solutions can be compared with upper bound results given by Utili [1], as shown in Table 1. The upper bound on the stability factor, N = γH/c, is obtained. When different φ values (i.e., φ = 10°, 20°) and different crack depth (i.e., H1/H = 0.2∼0.8) were considered in fissured slope analysis, it should be noted that the upper bound solutions presented in this study are consistent with the upper bound results of Utili [1].

In the pseudostatic condition (i.e., Kh ≠ 0,  ≠ 0), Utili and Abd [2] introduced the 2D upper bound theorem to conduct fissured slope stability, and the stability factor N was also chosen as the evaluation index for parameter analysis. The stability factor is plotted in Figure 2 for comparison with the parameters corresponding to β = 40°∼90° and Kh = 0.1∼0.4. In this case, the parameter λ is taken as 0.5, and φ = 30°. As can be seen in Figure 2, the trends of the two results are the same. In addition, the dotted lines, which represent the results of this study, are in good agreement with the solid lines, representing the results of Utili and Abd [2].

The two comparisons indicated that the solution derived from this study was reasonable for evaluating the stability of the 2D fissured slopes with soil anisotropy and nonhomogeneity subjected to seismic action.

3.2. Stability Charts

Based on comparative studies in different combinations of the related parameters, the calculation results of the stability factor N are given in Figure 3 for different coefficients of anisotropy and nonhomogeneity. We select four slopes with H1/H = 0.1, 0.2, 0.3, and 0.4 as examples. Figure 3 presents that the stability factor N as the x-coordinate is the nonhomogeneity coefficient n0. Considering different nonhomogeneity coefficients n0 (n0 = 0∼1.0) and different anisotropy coefficients k (k = 0.6∼1.4), five variation lines are presented in each figure. Furthermore, a fissured slope with the following properties was assumed: slope angle β = 60°, internal friction angle φ = 30°, seismic coefficient Kh = 0.3, and λ = 0.5. In Figure 3, the slope of the calculated line decreases with an increase in the crack depth H1/H. It is noted that the stability factor N increases linearly with an increase in the cohesion nonhomogeneity coefficient n0 and decreases as the cohesion anisotropy coefficient k increases. It implies that the impact of soil anisotropy and nonhomogeneity should be considered in the excavations of cuttings.

3.3. Determination of Crack Location

With a crack of known depth H1, through the computations, the most adverse location of the crack can be obtained, as shown in Figures 4 and 5. The dimensionless crack location L/H (L = L1L2) is plotted against the dimensionless crack depth δ/H for Kh = 0.1, 0.3 and λ = −0.5, 0.5. The nonhomogeneity coefficient n0 = 1.0 for Figure 4, and the anisotropy coefficient k ranges from 0.6 to 1.4 with increments of 0.2. In Figure 5, the anisotropy coefficient k = 1.0, and the nonhomogeneity coefficient n0 ranges from 0 to 1.0 with increments of 0.5. It can be seen in Figures 4 and 5 that the location of the crack L/H increases with increasing Kh. However, the location of the crack L/H decreases with increasing λ. It implies that, with an increase in the horizontal seismic acceleration coefficient Kh and with a decrease in the vertical seismic acceleration coefficient ( = λKh), the critical location of the crack became deeper (a smaller L2). It can also be seen that as the anisotropy coefficient k and nonhomogeneity coefficient n0 increased, the location of the crack L/H decreased. Meanwhile, it was found that the bigger the values of the known depth of the crack H1/H were, the lower the difference of the location of the crack L/H was. Especially when H1/H increased to a certain value, the location of the crack L/H with various anisotropy and nonhomogeneity coefficients reached the same value.

3.4. Solution of Safety Factors

To further study the impact of anisotropy and nonhomogeneity on the stability of given fissured slopes, the factor of safety FS can be obtained from a series of stability charts. To get the factor of safety FS, some conditions are considered in this section: A fissured slope has a height of 15 m, H1/H = 0.2, β = 60°, and γ = 18 kN/m3. The frictional cohesive soil properties φ = 30° and c = 40 kPa. The seismic coefficient Kh = 0.3 and λ = 0.5. The anisotropy coefficient k = 1.0 for case 1, and the cohesion nonhomogeneity coefficient n0 ranges from 0.2 to 1.0 with the increments of 0.2. In case 2, the cohesion nonhomogeneity coefficient n0 = 1.0, and the anisotropy coefficient k ranges from 0.6 to 1.4 with the increments of 0.2.

To explore the stability charts with the change of the anisotropy and nonhomogeneity coefficient for a fissured slope, Figures 6(a) and 6(b) present the different values of λ = c/γH tan φ by taking the parameter ψ = tan φ/FS as the x-coordinate. As expected, the values of λ were not linearly distributed with the parameter ψ. In Figure 6(a) for case 1, the change of the nonhomogeneity coefficient n0 significantly affected the values of λ and ψ. It can be concluded that the value of ψ increases as the nonhomogeneity coefficient n0 decreases, and also note that as ψ increases, the value of λ decreases. However, in Figure 6(a) for case 2, the anisotropy coefficient k has minor effects on the values of λ and ψ.

Using the presented stability charts (Figures 6(a) and 6(b)), one can easily get the factor of safety for given fissured slope and soil properties. A simple example is given here to illustrate it. Firstly, take case 1 as an example and calculate λ = c/γH tan φ = 0.257. Secondly, with the anisotropy coefficient k = 1.0 and the given nonhomogeneity coefficient n0 being 0.2, 0.6, and 1.0, it is determined from Figure 6(a) that the values of ψ are 0.688, 0.630, and 0.570, respectively. Finally, because the parameter φ is given, the factors of safety can be calculated as FS = tan φ/λ = 0.84, 0.91, and 1.01, respectively.

To completely study the influence of anisotropy and nonhomogeneity, the change of the factor of safety for a fissured slope is obtained from Figure 7, which covers the overall range of anisotropy and nonhomogeneity coefficients. Note that as the nonhomogeneity coefficient n0 increases, the factor of safety FS increases, which leads to the fissured slope becoming safe. However, the factor of safety FS decreases as the anisotropy coefficient k increases, which leads to the fissured slope becoming unsafe.

Considering different anisotropy and nonhomogeneity conditions, Figure 8 gives the values of the factor of safety FS for three different factors (Kh, H1/H, and λ). It can be seen from Figure 8(a) that as the seismic acceleration increases, the factor of safety FS decreases. Taking Figure 8(b) into consideration, an interesting phenomenon is that the factor of safety FS increases as the known depth of the crack H1/H increases. One possible reason is that the variation crack depth H1/H will correspond to different crack locations (L/H) when searching for the most dangerous sliding surface in upper bound analysis, resulting in the safety factors being uncertain. Meanwhile, a higher nonhomogeneity coefficient can lead to the factor of safety FS decrease. Figure 8(c) shows the effects of the parameter λ on the factor of safety FS. When the parameter λ is negative, the factor of safety FS obtained is larger than that when the parameter λ is positive. So it can be concluded that a fissured slope under vertical downward seismic acceleration is more dangerous than that under vertical upward seismic acceleration.

4. Conclusions

Based on the upper bound theorem of limit analysis, solutions are given to evaluate the stability of the fissured slopes subjected to seismic action with soil anisotropy and nonhomogeneity properties. A numerical optimization procedure is applied to predict the adverse location of a crack and the safety factor of a fissured slope. To better estimate the stability of the fissured slope, a series of stability charts are presented, and the procedures for getting the factor of safety are put forward. Besides, from this study, the following conclusions can be drawn:(1)The stability factor N increases as the nonhomogeneity coefficient n0 increases, but it decreases with the increasing anisotropy coefficient k. It is found from a series of stability charts that the influence of soil anisotropy and nonhomogeneity should not be ignored in the theoretical analysis of fissured slope stability.(2)The location of the crack L/H decreases as the anisotropy coefficient k and nonhomogeneity coefficient n0 increase. The seismicity has dramatic influence on the adverse location of the crack. As the horizontal seismic acceleration coefficient Kh increases or the vertical seismic acceleration coefficient decreases, the critical location of the crack becomes deeper. If the values of the known depth of the crack H1/H are bigger, the difference of the location of the crack L/H is lower.(3)The change of the nonhomogeneity coefficient n0 dramatically affected the values of λ and ψ, but the anisotropy coefficient k had minor effects on them. The factor of safety FS of the fissured slope obtained from the relevant stability charts is appropriate. An increase in the nonhomogeneity coefficient n0 and a decrease in the anisotropy coefficient k could lead to an increase in the factor of safety FS. As the seismic acceleration increases, the factor of safety FS decreases. Moreover, a fissured slope under vertical downward seismic acceleration is more dangerous than that under vertical upward seismic acceleration. It is valid under the conditions studied in this article.

Data Availability

The generated or analyzed data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest related to this paper.

Acknowledgments

This study was financially supported by the National Natural Science Foundation of China (Grant no. 51708310) and Shandong Provincial Natural Science Foundation, China (Grant no. ZR2017BEE066).