Mathematical Problems in Engineering

Volume 2015 (2015), Article ID 512648, 8 pages

http://dx.doi.org/10.1155/2015/512648

## Reliability Analysis of High Rockfill Dam Stability

School of Civil Engineering, Dalian University of Technology, Dalian 116024, China

Received 11 September 2014; Accepted 17 April 2015

Academic Editor: Rafael J. Villanueva

Copyright © 2015 Ping Yi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### 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.