MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 10.1155/2015/512648 512648 Research Article Reliability Analysis of High Rockfill Dam Stability Yi Ping Liu Jun Xu Chunlei Villanueva Rafael J. School of Civil Engineering Dalian University of Technology Dalian 116024 China dlut.edu.cn 2015 1692015 2015 11 09 2014 17 04 2015 1692015 2015 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.

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 ; Janbu 1973 ). 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 ; Bureau 2003 ; Yanmaz and Beser 2005 ). Duncan (2000)  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 ; Griffiths and Fenton 2004 ). 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)  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 ). 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 ) method or its extension, the point-estimate method (Rosenblueth 1975 ), 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 ). In others, two expressions are given based on the normal distribution and the lognormal distribution (Liang et al. 1999 ). 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 ). Baecher and Christian (2003)  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 ). 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 ). 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, g ( x ) , should be defined to identify the failure state ( g ( x ) < 0 ) and safety state ( g ( x ) > 0 ) , where x = ( x 1 , x 2 , , x n ) is a random variable vector. The following formulation of the performance function (Phoon 2008 ) is widely utilized and adopted in this paper: (1) g x = g x 1 , x 2 , , x n = F s x 1 , x 2 , , x n - 1 . 0 , where F s is the safety factor and the prescribed acceptable safety factor is 1.0 (Liang et al. 1999 ). The probability of failure can be defined as (2) P f = P g x < 0 = g x 0 f x d x in which f ( x ) 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 (3) P f Φ - β ; Φ ( · ) is the standard normal cumulate distribution function.

In slope stability reliability analysis, most previous reports have utilized FOSM (Hassan and Wolff 1999 ) method or its extension, the point-estimate method (Rosenblueth 1975 ), 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 ) and Baecher and Christian (2003)  have studied such problem in geotechnical engineering. Accordingly, FORM has gradually been accepted by researchers despite its complexity (Babu and Srivastava 2010 ).

In FORM, the original random vector x is transformed to a standard Gaussian vector u firstly, expressed as x = T ( u ) (Hohenbichler and Rackwitz 1981 ) and the performance function g ( x ) = g ( T ( u ) ) = G ( u ) . Then the reliability index β is the minimum distance from the coordinate origin to the limit state surface in u -space, and its computation is formulated as the following optimization problem: (4) min u s.t. G u = 0 so that β = u , where u 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 ). Among various methods that have been evaluated for application on structural reliability, an iterative algorithm, the HL-RF algorithm (Hasofer and Lind 1974 , Rackwitz and Flessler 1978 ), 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 ). The HL-RF algorithm computes the reliability index by the following iterative formulas: (5) β k = G u k - u G u k T u k u G u k , u k + 1 = - β k u G u k u G u k . 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 ), simple Janbu method (Janbu 1968 ), or the global analysis method (Zheng 2009 , Zheng 2012 ) 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)  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 ). 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, u , 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 ; Honjo and Kuroda 1991 ). However, some researchers have shown that these two kinds of slip surfaces do not generally coincide (Liang et al. 1999 ; Hassan and Wolff 1999 ). Here, we calculate β corresponding to the critical probabilistic surface. However, when the derivative u G ( u k ) 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 ).

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 ), the Shuangjiangkou core-wall rockfill dam with a height of 314 m (Wu et al. 2009 ) and the Nuozhadu core-wall rockfill dam with a height of 261.5 m (L.-H. Chen and Z.-Y. Chen 2007 ). 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 13. 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.

Comparison of stability analyses of the three examples.

Results computed by 3DSTAB Results in the literature
F s β F s β
Multilayered soil slope 1.525 5.121 1.500 4.963
Shuangjiangkou core-wall rockfill dam 1.813 5.353 5.573
Nuozhadu core-wall rockfill dam 1.949 4.185 1.888

Comparison of the critical slip circles in the example of a multilayered soil slope.

Result computed by 3DSTAB

Result in  (Chen, 2003)

Comparison of the critical slip circles of the Shuangjiangkou core-wall rockfill dam.

Result computed by 3DSTAB

Result in  (Wu et al. 2009)

Comparison of the critical slip circles of the Nuozhadu core-wall rockfill dam.

Result computed by 3DSTAB

Result in  (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: (6) τ f = c + σ n tan Ф in which τ f is the peak shear strength, σ n is the effective normal stress, c 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 ). Thus, the nonlinearity of the Mohr-Coulomb failure criterion cannot be neglected. Duncan (1980)  proposed that the drained friction angle Φ of rockfill can be evaluated from the following logarithmic expression: (7) Φ = Φ 0 - Δ Φ log σ 3 P a in which σ 3 is the smallest principal stress, P a is the atmospheric pressure, Φ 0 is the value of Φ when σ 3 is equal to P a , and Δ Ф is the reduction of Φ corresponding to a 10-fold increase in σ 3 . This logarithmic nonlinear shear strength criterion is widely utilized in China (Bai and Cui 1994 ; Lu and Yin 2004 ), and the Chinese “Design Code for Rolled Earth-rock Fill Dams” (SL274-2001)  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) , 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 ; Wu et al. 2009 ), the correlations are not considered here, and the obtained results are thus relatively conservative.

Statistical information of earth core.

γ (kN/m3) c (kN/m3) Φ (°)
Mean value μ 20.7 23.0 30.3
Standard deviation σ 1.3 2.8 1.3
Variation coefficient δ (%) 6.3 12.2 4.3
Distribution type Normal Lognormal Lognormal

Statistical information of rockfill.

Main rockfill Subrockfill
γ (kN/m3) Φ 0 (°) Δ Φ (°) γ (kN/m3) Φ 0 (°) Δ Φ (°)
Mean value μ 21.0 51.2 8.2 20.8 48.8 9.4
Standard deviation σ 0.8 4.3 3.2 0.7 4.6 2.3
Variation coefficient δ (%) 3.8 8.4 39.0 3.3 9.4 24.5
Distribution type Normal Normal Normal Normal Normal Normal

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 R , and negative values may be given to a parameter that is positive in nature, such as Φ 0 . This would give an unreasonable lower value to the reliability index (Liang et al. 1999 ). 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 ). The upper limit is μ + 3 σ , and the lower limit is the larger one of μ - 3 σ 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 H varies from 100 m to 350 m in intervals of 50 m. The slope ratio is expressed as 1 : m, where m 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 ( m H + 20 ), 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.

Cross section of a concrete faced rockfill dam.

Cross section of an earth core rockfill dam.

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 ( H ) and slope ratio parameters ( m ) 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.

Safety factors and reliability indices of downstream slopes of concrete faced rockfill dams.

Safety factor

Reliability index

Safety factors and reliability indices of upstream slopes of earth core dams.

Safety factor

Reliability index

Safety factors and reliability indices of downstream slopes of earth core rockfill dams.

Safety factor

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 ( H ) increases, as mentioned by Cao and Chen (1991) . At the same time, both quantities decrease as the slope ratio parameter ( m ) 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 ( H ) 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 ( H ) increases and the slope factor ( m ) 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 F s 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)  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

( 1 ) 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.

( 2 ) 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.

( 3 ) 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.

Chen Z.-Y. Morgenstern N. R. Extensions to the generalized method of slices for stability analysis Canadian Geotechnical Journal 1983 20 1 104 119 10.1139/t83-010 2-s2.0-0020704846 Janbu N. Hirschfeld R. C. Poulos S. J. Slope stability computation Embankment-Dam Engineering 1973 New York, NY, USA John Wiley & Sons 47 86 Wolff T. F. Probabilistic slope stability in theory and practice Proceedings of the Uncertainty in the Geologic Environment: From Theory to Practice July-August 1996 Madision, Wis, USA ASCE 419 433 Bureau G. J. Chen W. F. Scawthorn C. Dams and appurtenant facilities Earthquake Engineering Handbook 2003 CRC Press Yanmaz A. M. Beser M. R. On the reliability-based safety analysis of Porsuk dam Turkish Journal of Engineering & Environmental Sciences 2005 29 309 320 Duncan J. M. Factors of safety and reliability in geotechnical engineering Journal of Geotechnical and Geoenvironmental Engineering 2000 126 4 307 316 10.1061/(asce)1090-0241(2000)126:4(307) 2-s2.0-0034173877 Hassan A. M. Wolff T. F. Search algorithm for minimum reliability index of earth slopes Journal of Geotechnical and Geoenvironmental Engineering 1999 125 4 301 308 10.1061/(asce)1090-0241(1999)125:4(301) 2-s2.0-0033120954 Griffiths D. V. Fenton G. A. Probabilistic slope stability analysis by finite elements Journal of Geotechnical and Geoenvironmental Engineering 2004 130 5 507 518 10.1061/(asce)1090-0241(2004)130:5(507) 2-s2.0-2342514143 Xu B. Low B. K. Probabilistic stability analyses of embankments based on finite-element method Journal of Geotechnical and Geoenvironmental Engineering 2006 132 11 1444 1454 10.1061/(asce)1090-0241(2006)132:11(1444) 2-s2.0-33751242728 Luo X. F. Li X. Zhou J. Cheng T. A Kriging-based hybrid optimization algorithm for slope reliability analysis Structural Safety 2012 34 1 401 406 10.1016/j.strusafe.2011.09.004 2-s2.0-80355131622 Rosenblueth E. Point estimates for probability moments Proceedings of the National Academy of Sciences of the United States of America 1975 72 10 3812 3814 10.1073/pnas.72.10.3812 Liang R. Y. Nusier O. K. Malkawi A. H. A reliability based approach for evaluating the slope stability of embankment dams Engineering Geology 1999 54 3-4 271 285 10.1016/s0013-7952(99)00017-4 2-s2.0-0032711283 Hasofer A. M. Lind N. C. Exact and invariant second-moment code format Journal of the Engineering Mechanics Division 1974 100 1 111 121 2-s2.0-0016026184 Baecher G. B. Christian J. T. Reliability and Statistics in Geotechnical Engineering 2003 Chichester, UK John Wiley & Sons Babu G. L. S. Srivastava A. Reliability analysis of earth dams Journal of Geotechnical and Geoenvironmental Engineering 2010 136 7 995 998 10.1061/(asce)gt.1943-5606.0000313 2-s2.0-77955619376 Cho S. E. Probabilistic stability analyses of slopes using the ANN-based response surface Computers and Geotechnics 2009 36 5 787 797 10.1016/j.compgeo.2009.01.003 2-s2.0-64849106854 Phoon K.-K. Reliability-Based Design in Geotechnical Engineering 2008 New York, NY, USA Taylor and Francis Hohenbichler M. Rackwitz R. Non-normal dependent vectors in structural safety Journal of Engineering Mechanics 1981 107 6 1227 1239 2-s2.0-0019717070 Val D. Bljuger F. Yankelevsky D. Optimization problem solution in reliability analysis of reinforced concrete structures Computers & Structures 1996 60 3 351 355 10.1016/0045-7949(95)00390-8 2-s2.0-0030567803 Rackwitz R. Flessler B. Structural reliability under combined random load sequences Computers & Structures 1978 9 5 489 494 10.1016/0045-7949(78)90046-9 2-s2.0-0018033230 Liu P.-L. der Kiureghian A. Optimization algorithms for structural reliability Structural Safety 1991 9 3 161 177 10.1016/0167-4730(91)90041-7 2-s2.0-0026107548 Bishop A. W. The use of the slip circle in the stability analysis of slopes Géotechnique 1955 5 1 7 17 10.1680/geot.1955.5.1.7 Janbu N. Slope stability computations Soil Mechanics and Foundation Engineering Report 1968 Trondheim, Norway The Technical University of Norway Zheng H. Eigenvalue problem from the stability analysis of slopes Journal of Geotechnical and Geoenvironmental Engineering 2009 135 5 647 656 10.1061/(asce)gt.1943-5606.0000071 Zheng H. A threedimensional rigorous method for stability analysis of landslides Engineering Geology 2012 145-146 30 40 10.1016/j.enggeo.2012.06.010 Duncan J. M. State of the art: limit equilibrium and finite-element analysis of slopes Journal of Geotechnical Engineering 1996 122 7 577 596 10.1061/(asce)0733-9410(1996)122:7(577) 2-s2.0-0001690750 Li L. Chi S. Lin G. Search for the critical slip surface using mixed version of particle swarm optimization algorithms Industrial Construction 2007 37 2 55 59 Calle E. O. F. Probabilistic analysis of stability of earth slopes Proceedings of the 11th International Conference on Soil Mechanics and Foundation Engineering August 1985 San Francisco, Calif, USA Balkema 809 812 Honjo Y. Kuroda K. A new look at fluctuating geotechnical data for reliability design Soils and Foundations 1991 31 1 110 120 10.3208/sandf1972.31.110 2-s2.0-0026120660 Yi P. Cheng G. Jiang L. A sequential approximate programming strategy for performance-measure-based probabilistic structural design optimization Structural Safety 2008 30 2 91 109 10.1016/j.strusafe.2006.08.003 2-s2.0-36348931585 Chen Z. Y. Soil Slope Stability Analysis: Theory, Methods and Programs 2003 Beijing, China Water Pub Press (Chinese) Wu Z. Y. Chen J. K. Xu W. L. Wen B. Zhang P. Research on methodology of reliability analysis of high rock-fill dam slope stability using nonlinear strength indexes and its engineering application Chinese Journal of Rock Mechanics and Engineering 2009 28 1 130 137 2-s2.0-60149095562 Chen L.-H. Chen Z.-Y. Effect of nonlinear strength of rockfill on slope stability of high earth-rock dam Rock and Soil Mechanics 2007 28 9 1807 1810 2-s2.0-35349008500 Indraratna B. Wijewardena L. S. S. Balasubramaniam A. S. Large-scale triaxial testing of grey wacke rockfill Geotechnique 1993 43 1 37 51 10.1680/geot.1993.43.1.37 2-s2.0-0027341378 Duncan J. M. Strength, stress-strain and bulk modulus parameters for finite element analyses of stresses and movements in soil masses 1980 Berkeley, Calif, USA The University of California Bai S. T. Cui Y. H. The mechanical properties of rockfill Journal of Hydroelectric Engineering 1994 58 3 21 30 Lu Q. F. Yin Z. Z. Effect of strength nonlinearity on slope stability of high rockfill dam Chinese Journal of Rock Mechanics and Engineering 2004 23 16 2708 2711 2-s2.0-4644335340 The Professional Standards Compilation Group of People's Republic of China Design Code for Rolled Earth-Rock Fill Dams (SL274-2001) 2002 Beijing, China China Water Power Press (Chinese) Xu C. Study on Reliability of Slope Stability of High Rockfill Dam 2010 Dalian, China Dalian University of Technology (Chinese) Most T. Knabe T. Reliability analysis of the bearing failure problem considering uncertain stochastic parameters Computers and Geotechnics 2010 37 3 299 310 10.1016/j.compgeo.2009.11.003 2-s2.0-76549120180 Cao K. M. Chen Z. W. The stability analysis of rockfill dam using non-linear strength criterion Water Resources and Hydropower Engineering 1991 8 9 37 42