#### Abstract

A program 3DSTAB combining slope stability analysis and reliability analysis is developed and validated. In this program, the limit equilibrium method is utilized to calculate safety factors of critical slip surfaces. The first-order reliability method is used to compute reliability indexes corresponding to critical probabilistic surfaces. When derivatives of the performance function are calculated by finite difference method, the previous iteration’s critical slip surface is saved and used. This sequential approximation strategy notably improves efficiency. Using this program, the stability reliability analyses of concrete faced rockfill dams and earth core rockfill dams with different heights and different slope ratios are performed. The results show that both safety factors and reliability indexes decrease as the dam’s slope increases at a constant height and as the dam’s height increases at a constant slope. They decrease dramatically as the dam height increases from 100 m to 200 m while they decrease slowly once the dam height exceeds 250 m, which deserves attention. Additionally, both safety factors and reliability indexes of the upstream slope of earth core rockfill dams are higher than that of the downstream slope. Thus, the downstream slope stability is the key failure mode for earth core rockfill dams.

#### 1. Introduction

Rockfill dams are commonly used geotechnical infrastructures for water management. One of the critical aspects of rockfill dam design is stability analysis, that is, the computation of safety factors (Chen and Morgenstern 1983 [1]; Janbu 1973 [2]). A deterministic approach is traditionally utilized for this analysis. However, as natural materials, rockfill and earth core exhibit large uncertainties in their shear strength parameters, dry bulk density, and other pertinent properties, which cannot be handled in the traditional deterministic methods. Therefore, the application of probabilistic reliability concepts to the stability analysis of dams has drawn increasing attention over the past two decades (Wolff 1996 [3]; Bureau 2003 [4]; Yanmaz and Beser 2005 [5]). Duncan (2000) [6] pointed out that the reliability analysis offers a useful supplement to conventional stability analyses, because the resultant reliability index contains more information than the deterministic safety factor.

The limit equilibrium method (LEM) and the strength reduction method (SRM) on the basis of the finite element method (FEM) or finite difference method (FDM) are currently popular methods among engineers for slope stability analysis (Hassan and Wolff 1999 [7]; Griffiths and Fenton 2004 [8]). Whether LEM or SRM is used, the stability analysis of rockfill dams is a time-consuming process and the calculated safety factor is an implicit function of basic variables, such as materials’ shear strength parameters and dry bulk density. Then in the reliability analysis, the performance function is implicit and many iterations of calculating the safety factor have to be performed to obtain the reliability index. Therefore, despite its potential value, reliability theory has not been widely adopted in geotechnical engineering because of huge computational cost.

The response surface method (RSM) has been developed to deal with implicit performance functions and Xu and Low (2006) [9] utilized RSM combined with FEM to calculate the reliability index of slopes. However, the surrogate performance function from the RSM may cause a deviation from the exact model (Luo et al. 2012 [10]). When no such surrogate performance function is used, most previous reports have utilized the mean value first-order second-moment (FOSM, Hassan and Wolff 1999 [7]) method or its extension, the point-estimate method (Rosenblueth 1975 [11]), to calculate the approximate reliability index of slope stability. In some papers, the distributions of the random variables are not mentioned at all, and the reliability index is calculated by the performance function’s mean value divided by its standard deviation (Duncan 2000 [6]). In others, two expressions are given based on the normal distribution and the lognormal distribution (Liang et al. 1999 [12]). However, it is well known in reliability analysis that different values of the reliability index might be obtained for different mathematical forms of the same limit state function when FOSM is used (Hasofer and Lind 1974 [13]). Baecher and Christian (2003) [14] have studied such problem in geotechnical engineering. Accordingly, the first-order reliability method (FORM) is recommended and has gradually been accepted by researchers despite its complexity (Babu and Srivastava 2010 [15]). In this paper, FORM is utilized.

Although there are some commercial software programs, for example, PLAXIS (FEM) and FLAC (FDM), for the slope stability analysis, software package combining the slope stability and reliability analysis is rare (Cho 2009 [16]). In this paper, a self-developed program 3DSTAB, which integrates slope stability analysis and reliability analysis, is introduced. LEM is employed in this program to calculate safety factors and FORM, more precisely, the HL-RF iterative algorithm, is utilized to compute reliability indices. Three examples reported in recent literature are studied and the comparisons prove the validity and accuracy of 3DSTAB. Then two typical kinds of rockfill dams, that is, concrete faced rockfill dams and earth core rockfill dams, are considered. Rockfill dams with different heights over 100 m and different slope ratios are analyzed to study the relationships of the safety factors and reliability indices to the dam height and dam slope ratio parameter.

#### 2. Reliability Analysis

To perform the reliability analysis of the slope stability of dams, a performance function or limit state function, , should be defined to identify the failure state and safety state , where is a random variable vector. The following formulation of the performance function (Phoon 2008 [17]) is widely utilized and adopted in this paper:where is the safety factor and the prescribed acceptable safety factor is 1.0 (Liang et al. 1999 [12]). The probability of failure can be defined asin which is the joint probability density function of** x**. Because the multidimensional integral in (2) can be very difficult and nearly impossible, the reliability index is generally calculated in engineering, and the failure probability is estimated by is the standard normal cumulate distribution function.

In slope stability reliability analysis, most previous reports have utilized FOSM (Hassan and Wolff 1999 [7]) method or its extension, the point-estimate method (Rosenblueth 1975 [11]), to calculate the reliability index to avoid heavy computational burden. However when FOSM is used, different values of the reliability index might be obtained for different mathematical forms of the same limit state function (Hasofer and Lind 1974 [13]) and Baecher and Christian (2003) [14] have studied such problem in geotechnical engineering. Accordingly, FORM has gradually been accepted by researchers despite its complexity (Babu and Srivastava 2010 [15]).

In FORM, the original random vector is transformed to a standard Gaussian vector firstly, expressed as (Hohenbichler and Rackwitz 1981 [18]) and the performance function . Then the reliability index is the minimum distance from the coordinate origin to the limit state surface in -space, and its computation is formulated as the following optimization problem:where is the design point on the limit state surface in** u**-space. The design point can be located by various optimization algorithms, such as the gradient projection (GP) method, the augmented Lagrangian (AL) method, and the sequential quadratic programming (SQP) method (Val et al. 1996 [19]). Among various methods that have been evaluated for application on structural reliability, an iterative algorithm, the HL-RF algorithm (Hasofer and Lind 1974 [13], Rackwitz and Flessler 1978 [20]), is found to be very efficient as it requires the least amount of storage and computation in each step in comparison with other methods (Liu and der Kiureghian 1991 [21]). The HL-RF algorithm computes the reliability index by the following iterative formulas:The iteration is completed when the absolute change of the design points between two subsequent iterations is smaller than a prescribed small quantity, for example, 0.001. The HL-RF algorithm can generally obtain the results with enough accuracy after some iterations and is utilized in this paper.

#### 3. 3DSTAB Program

The slope stability analysis and reliability analysis are combined in the self-developed program 3DSTAB. LEM, which has received wide acceptance because of its simplicity, is utilized to calculate the slope stability safety factor. FORM, more precisely, the HL-RF iterative algorithm, is utilized to compute the reliability index.

The Bishop method (Bishop 1955 [22]), simple Janbu method (Janbu 1968 [23]), or the global analysis method (Zheng 2009 [24], Zheng 2012 [25]) can be chosen under static conditions. Although the program can perform 2D and 3D slope stability analyses, just 2D analysis is performed in this paper because Duncan (1996) [26] pointed out that the safety factors resulting from 3D analyses are normally greater than those from 2D analyses. To calculate the safety factor of the slope stability, mean values are assigned to the variables and the particle swarm optimization (PSO) method is applied to search the critical slip surface, which is referred to as the critical deterministic surface (Li et al. 2007 [27]). In reliability analysis, the critical slip surface is searched, and the corresponding safety factor is determined for every realization of the random variables in every iteration. Several iterations are needed to calculate the reliability index through the HL-RF algorithm. The critical slip surface for the final realization, , is referred to as the critical probabilistic slip surface. Some previous papers only calculated the reliability index corresponding to the critical deterministic surface (Calle 1985 [28]; Honjo and Kuroda 1991 [29]). However, some researchers have shown that these two kinds of slip surfaces do not generally coincide (Liang et al. 1999 [12]; Hassan and Wolff 1999 [7]). Here, we calculate corresponding to the critical probabilistic surface. However, when the derivative is computed by the finite difference method in every iteration, we do not search the critical slip surface anymore, and the previous critical slip surface is saved and utilized. This kind of sequential approximation strategy is widely utilized in structural optimization and reliability-based optimization (Yi et al. 2008 [30]).

To prove the accuracy and validity of 3DSTAB, three examples that have recently appeared in literature are restudied and compared. The examples include an example of multilayered soil (Chen 2003 [31]), the Shuangjiangkou core-wall rockfill dam with a height of 314 m (Wu et al. 2009 [32]) and the Nuozhadu core-wall rockfill dam with a height of 261.5 m (L.-H. Chen and Z.-Y. Chen 2007 [33]). The computational models and the statistical information of the random parameters are the same as those in the corresponding literature. The comparisons of the critical slip circles are shown in Figures 1–3. The safety factors and the reliability indices are listed in Table 1. The figures and the table show that the results computed by 3DSTAB program are similar to the results reported in other publications and are accurate and reliable.

**(a) Result computed by 3DSTAB**

**(b) Result in [31] (Chen, 2003)**

**(a) Result computed by 3DSTAB**

**(b) Result in [32] (Wu et al. 2009)**

**(a) Result computed by 3DSTAB**

**(b) Result in [33] (L.-H. Chen and Z.-Y. Chen 2007)**

#### 4. Stability Reliability Analysis of High Rockfill Dams

As natural materials, rockfill and earth core exhibit large uncertainties. In this paper, the dry bulk densities and shear strength parameters, including the linear strength parameters of earth core material and the nonlinear strength parameters of rockfill material, are taken as random variables. Based on the statistical properties of the random variables reported in literature, stability reliability analyses of two typical kinds of rockfill dams with different heights over 100 m and different slope ratios are performed using the self-developed program 3DSTAB. The effects of dam height and slope ratio on the safety factor and the reliability index of concrete faced rockfill dams and earth core rockfill dams are discussed.

##### 4.1. Statistical Properties of Random Variables

The linear Mohr-Coulomb failure criterion is conventionally utilized to perform dam stability analyses as follows:in which is the peak shear strength, is the effective normal stress, is the cohesion or snap-in force, and is the internal frictional angle. However, the results of many large-scale triaxial tests indicate that Mohr failure envelopes for rockfill exhibit significant nonlinearity, particularly under low and medium normal stresses. And the slope failure of a rockfill dam is just generally associated with relatively low normal stresses (Indraratna et al. 1993 [34]). Thus, the nonlinearity of the Mohr-Coulomb failure criterion cannot be neglected. Duncan (1980) [35] proposed that the drained friction angle of rockfill can be evaluated from the following logarithmic expression:in which is the smallest principal stress, is the atmospheric pressure, is the value of when is equal to , and Δ is the reduction of corresponding to a 10-fold increase in . This logarithmic nonlinear shear strength criterion is widely utilized in China (Bai and Cui 1994 [36]; Lu and Yin 2004 [37]), and the Chinese “Design Code for Rolled Earth-rock Fill Dams” (SL274-2001) [38] recommends that the logarithmic expression (7) should be utilized to describe the nonlinear shear strength of rough grain materials, such as rockfill.

In this paper, we consider earth core material to have linear shear strength expressed in (6) and the main and subrockfill material to have the logarithmic nonlinear shear strength expressed in (7). The shear strength parameters in these expressions and the dry bulk density, , are taken as random parameters and their statistical information, which is reported by Xu (2010) [39], is listed in Tables 2 and 3. Because considering correlations between the parameters would only marginally increase the slope stability and its reliability (Babu and Srivastava 2010 [15]; Wu et al. 2009 [32]), the correlations are not considered here, and the obtained results are thus relatively conservative.

From Tables 2 and 3, we can see that most of the stochastic variables obey normal distribution. It has been pointed out that a normally distributed variable is defined on the real number space , and negative values may be given to a parameter that is positive in nature, such as . This would give an unreasonable lower value to the reliability index (Liang et al. 1999 [12]). The situation will be worse when the variation coefficient of the parameter is not small. So a truncated normal distribution is adopted in this study for all normally distributed parameters to avoid sampling negative values (Most and Knabe 2010 [40]). The upper limit is , and the lower limit is the larger one of and zero, in which is the mean value and is the standard deviation.

##### 4.2. Stability Reliability Analysis

Two kinds of rockfill dams, that is, concrete faced rockfill dams and earth core rockfill dams, are analyzed and their typical cross-sectional shapes are shown in Figures 4 and 5. The dam height varies from 100 m to 350 m in intervals of 50 m. The slope ratio is expressed as 1 :* m,* where can be 1.3, 1.4, 1.5, or 1.6 for concrete faced rockfill dams and 1.8, 1.9, or 2.0 for earth core rockfill dams, respectively. Cushion zone and transition zone are ignored due to little effect on slope stability. The Bishop method is utilized to calculate safety factors under static conditions. The number of soil strips is (), and the soil strip width is less than 1.0 m. The reliability index is computed through the HL-RF iterative algorithm with a maximum of 20 iterations.

The water level is set to the maximum pool level. The stability of the downstream slope is considered for concrete faced rockfill dams. The safety factors and reliability indices at different heights () and slope ratio parameters () are shown in Figure 6. For earth core rockfill dams, the stabilities of both the upstream and downstream slopes under steady-state seepage conditions are calculated. The results are shown in Figures 7 and 8.

**(a) Safety factor**

**(b) Reliability index**

**(a) Safety factor**

**(b) Reliability index**

**(a) Safety factor**

**(b) Reliability index**

Figure 6 shows that the safety factors of the downstream slopes of concrete faced rockfill dams range from 1.34 to 1.85 and the reliability indices range from 1.60 to 5.67. Both quantities decrease as the dam height () increases, as mentioned by Cao and Chen (1991) [41]. At the same time, both quantities decrease as the slope ratio parameter () decreases; that is, the safety factors and the reliability indices decrease as the dam slope becomes steeper. This corresponds with the fact that the higher and the steeper dams are, the more dangerous they become. It should be noted in Figure 6 that when the dam height is relatively low (i.e., from 100 m to 200 m), the safety factors and the reliability indices decrease dramatically as the dam height increases. Strangely, when the dam height () exceeds 250 m, the safety factors and the reliability indices decrease very slowly with further increases of the dam height. The stability reliability analyses of both the upstream and downstream slopes of earth core rockfill dams present the same behavior, and this deserves extra attention.

Figures 7 and 8 show the safety factors and reliability indices of the upstream and downstream slopes of earth core rockfill dams, respectively. These quantities also decrease as the dam height () increases and the slope factor () decreases. The safety factors of the upstream slopes of earth core rockfill dams range from 1.80 to 2.35, and the reliability indices range from 4.03 to 11.52. The ranges are from 1.70 to 2.23 for and from 3.20 to 8.37 for of downstream slopes. Both the safety factors and reliability indices of the upstream slopes of earth core rockfill dams are higher than those of the downstream slopes. Thus, the stability of the downstream slope is the key failure mode of earth core rockfill dams.

Figures 6 to 8 also show that the safety factors and reliability indices of earth core rockfill dams are higher than those of concrete faced rockfill dams. This is mainly because earth core rockfill dams considered here have relatively gentle slopes compared to concrete faced rockfill dams.

It should be pointed out that the reliability indices in this section are a bit small, especially for dams higher than 250 m. This happens because the statistical parameters reported by Xu (2010) [39] are the statistical results from 44 dams in China, and their divergences are relatively larger than those of a single dam. However, the behavior of the quantities is more important than their specific values.

#### 5. Conclusions

A self-developed program 3DSTAB combining the slope stability analysis and reliability analysis is introduced. In this program, LEM is utilized to calculate the safety factor of the critical slip surface under static conditions. FORM, more precisely, the HL-RF iterative algorithm, is utilized to compute the reliability index corresponding to the critical probabilistic surface. When the derivatives of the performance function are calculated by the finite difference method in every iteration, a sequential approximation strategy is utilized. In other words, we do not search the critical slip surface under this circumstance, and the previous critical slip surface is saved and used, which improves the efficiency notably. Three examples reported in recent literature are studied and the comparisons prove the accuracy and validity of 3DSTAB.

The dry bulk densities and shear strength parameters, including the linear strength parameters of earth core material and the nonlinear strength parameters of rockfill material, are taken as random variables. The truncated normal distribution is used for all normally distributed parameters to prevent negative sampling values. Because considering correlations between the parameters would only marginally increase the slope stability and its reliability, the correlations are not considered in this paper.

Using the self-developed program 3DSTAB, the stability reliability analyses of concrete faced rockfill dams and earth core rockfill dams with different heights over 100 m and different slope ratios are performed. The results show that both the safety factors and reliability indices decrease as the dam’s slope increases at a constant dam height and as the dam height increases at a constant slope ratio. More precisely, they decrease dramatically as the dam height increases from 100 m to 200 m, while decrease slowly once the dam height exceeds 250 m, which deserves extra attention.

Both the safety factors and reliability indices of the upstream slopes of earth core rockfill dams are higher than those of the downstream slopes. Thus, the stability of the downstream slope is the key failure mode of earth core rockfill dams.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This work is funded by the Key Project of Chinese National Programs for Fundamental Research and Development (2015CB057703), the National Natural Scientific Foundation of China (51138001), and the Fundamental Research Funds for the Central Universities (DUT14ZD206), which are gratefully acknowledged.