MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 10.1155/2014/601418 601418 Research Article Evaluation Standard for Safety Coefficient of Roller Compacted Concrete Dam Based on Finite Element Method Li Bo 1 Zhang Zhijun 2 Liu Yan 3 Yang Shengmei 1 Liu Yan-Jun 1 Engineering Safety and Disaster Prevention Department Changjiang River Scientific Research Institute, Wuhan 430010 China crsri.cn 2 Construction Design Department Changjiang Institute of Survey, Planning, Design and Research, Wuhan 430010 China 3 Center of Construction Management & Quality & Safety Supervision MWR, Beijing 100038 China 2014 672014 2014 16 03 2014 27 05 2014 10 06 2014 6 7 2014 2014 Copyright © 2014 Bo Li 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.

The lack of evaluation standard for safety coefficient based on finite element method (FEM) limits the wide application of FEM in roller compacted concrete dam (RCCD). In this paper, the strength reserve factor (SRF) method is adopted to simulate gradual failure and possible unstable modes of RCCD system. The entropy theory and catastrophe theory are used to obtain the ultimate bearing resistance and failure criterion of the RCCD. The most dangerous sliding plane for RCCD failure is found using the Latin hypercube sampling (LHS) and auxiliary analysis of partial least squares regression (PLSR). Finally a method for determining the evaluation standard of RCCD safety coefficient based on FEM is put forward using least squares support vector machines (LSSVM) and particle swarm optimization (PSO). The proposed method is applied to safety coefficient analysis of the Longtan RCCD in China. The calculation shows that RCCD failure is closely related to RCCD interface strength, and the Longtan RCCD is safe in the design condition. Considering RCCD failure characteristic and combining the advantages of several excellent algorithms, the proposed method determines the evaluation standard for safety coefficient of RCCD based on FEM for the first time and can be popularized to any RCCD.

1. Introduction

Roller compacted concrete dams (RCCDs) are now constructed with great frequency in the world. Over the past two decades, many high RCCDs have been built, such as the 130 m high Baise dam and the 200.5 m high Guangzhao dam in China. Currently the highest dam of this type is Longtan Dam at 216 m in China, with Diamer-Bhasha Dam planned at 272 m in Pakistan .

The RCCD is mainly composed of the body and interface. The interface plays an important role as the weak point of the stability of the RCCD. The stability and safety of a RCCD to guard against sliding along the dam interface are key factors that must be addressed in the design stage of the RCCD. Generally speaking, analysis methods for the RCCD stability should include the rigid-body limiting equilibrium method (LEM) and finite element method (FEM). LEM is normally adopted before considering FEM. Ministry of Water Resources  presents evaluation standards for the LEM. Compared with the LEM, the FEM used in the calculation of antisliding analysis does not require any assumption of the sliding plane, but safety coefficient based on FEM does not have corresponding evaluation standard, which limits the usage of FEM in RCCD. In order to popularize FEM to RCCD more widely than ever, it is necessary to study the evaluation standard for RCCD safety coefficient based on FEM.

When the dam is in critical state, small perturbation will lead the dam away from the equilibrium state into the instability quickly. According to the dissipative structure theory, in the course of dam failure, some abnormal deformations appear certainly with the reduction of system entropy . Thereby, in this paper, RCCD failure criterion is put forward combining the entropy theory and catastrophe theory.

RCCD failure is formed under a variety of combined factors. The sliding mechanism and failure mode are determined by the mechanical parameters (deformation and strength parameters) of RCCD interface. The sensitivity analysis of effect factors for the stability of RCCD interface can determine the major control plane of the RCCD, namely, the most dangerous sliding plane which has great significance for RCCD antisliding stability analysis. In this paper, sensitivity analysis for RCCD interface stability based on the RCCD failure criterion with a hybrid algorithm is proposed. Latin hypercube sampling (LHS) was introduced by McKay et al.  as a means of generating a distribution of plausible collections of parameter values from a multidimensional distribution and has been further developed for different purposes by several researchers [5, 6]. Partial least squares regression (PLSR) is a multivariate statistical data analysis method, which can establish effectively the model for the data involving multicollinearity or small sample size [7, 8]. In order to combine the advantages of LHS and PLSR, a hybrid algorithm which combines LHS with auxiliary analysis of PLSR is proposed for the sensitivity analysis. Through the sensitivity analysis of effect factors, the most dangerous sliding plane of RCCD failure is determined. On this basis, the evaluation standard for RCCD safety coefficient based on FEM is proposed with a hybrid optimization algorithm. Least squares support vector machines (LSSVM) are least squares versions of support vector machines (SVM), which are a set of related supervised learning methods that analyze data and recognize patterns and are used for classification and regression analysis . Particle swarm optimization (PSO) is originally attributed to Kennedy and Eberhart who were illumed by the behavior of bird swarms in 1995 . It has been demonstrated that PSO outperforms genetic algorithm in multivariable function optimization . A hybrid optimization algorithm which combines LHS, LSSVM, and PSO is proposed for evaluation standard of RCCD safety coefficient based on FEM.

This paper reports the analyses of evaluation standard for safety coefficient of the Longtan RCCD based on FEM. The result shows that the Longtan RCCD is safe in the design condition.

2. RCCD Failure Criterion 2.1. RCC Yield Criterion

The transversely isotropic and elastoplastic model is used for the constitutive relationship of RCCD material. The isotropic and elastoplastic model is used for the constitutive relationship of the rock material of the dam foundation. The most typical Drucker-Prager (D-P) criterion [12, 13] in geotechnical engineering is used for the yield criterion of the RCCD and rock of the dam foundation: (1)F=αI1+J2-K=0, where I1 and J2 are the first invariant of stress tensor and the second invariant of deviatoric stress tensor, respectively. Both  α  and K are strength constants, which are related to cohesion c and internal friction angle  φ. Different yield criterions can be realized in finite element calculations through changing the expressions of  α and K. In this paper, the D-P criterion of equivalent-area cones is used.

2.2. Strength Reserve Factor Method

Zienkiewicz and Pande  firstly introduced strength reserve factor (SRF) method to elastoplastic FEM analysis. The dam and dam foundation have strength reserve ability. When shear strength parameters (cohesion c and internal friction coefficient f) of the dam decrease, the dam will lose strength and become instable. Basic principle of SRF method is that the cohesion c and internal friction coefficient f are divided by the reduction coefficient KR: (2)cR=cKR,fR=fKR.

With the gradual increasing of the reduction coefficient KR, the process of RCCD failure can be obtained, and the dam will reach a limit state.

2.3. Entropy Catastrophe Criterion of RCCD Failure

The process of RCCD failure is the evolution process from the “disorder” to the “order.” Because of the correspondence between the system entropy and the system disorder degree, the safety state of RCCD can be easily judged by calculating the entropy value of RCCD in this process .

The information entropy is a measure of the disorder associated with a random variable. The concept of information entropy was introduced by Shannon , which is computed with the algorithm (3)S=-ϕinpilog2pi, where S is the information entropy and pi is the probability of a discrete random variable.

In the process of RCCD failure, the dam interface is main control plane. Therefore, the interface entropy is mainly researched here. The total displacement U of the interfaces can be expressed as follows: (4)U=inui,(i=1,2,,n),λi=uiU,(i=1,2,,n), where ui is the displacement of the element (FEM element) i and λi expresses the distribution of the interface displacement. Obviously, the physical variable λi has nonnegative characteristic, which is the same as the probability mass function. According to (3), the information entropy S of the interface displacements is defined as follows: (5)S=-ϕinλilnλi.

With the reduction of the strength parameters by SRF method, the degree of order increases, and the entropy decreases. When the entropy decreases to certain value, the instable interface will be completely damaged.

Because the entropy S changes with the reduction coefficient KR, it can be shown as follows: (6)S(1)(x)=a0+a1x+a2x2+a3x3++anxn, where a0,a1,,an are undetermined coefficients and x is the reduction coefficient KR. Using multiple regression analysis, the coefficients a0,a1,,a5 can be obtained.

The derivation of (6) gets the following equation: (7)V=a1+2a2x+3a3x2+4a4x3+5a5x4, where V is the potential function; when b0=a1, b1=2a2, b2=3a3, b3=4a4, and b4=5a5, V can be expressed as (8)V=b0+b1x+b2x2+b3x3+b4x4.

When p=b3/4b4, x=y-p, the cube is eliminated: (9)V=c4y4+c2y2+c1y+c0, where c4=b4, c2=6p2b4-3pb3+b2, c1=-4p3b4+3p2b3-2pb2+b1, and c0=p4b4-p3b3+p2b2-pb1+b0. When y=(Z/4c44)(c4>0) or y=(Z/-4c44)(c4<0), (9) can be written as follows: (10)V=14Z4+12uZ2+vZ, or (11)V=-14Z4+12uZ2+vZ, where u=c2/|c4|, v=c1/4|c4|4.

Equation (10) is the standard potential function of regular cusp catastrophe. Equation (11) is the standard potential function of dual cusp catastrophe [16, 17].

When V/Z=0, (10) and (11) can be expressed, respectively, as follows: (12)Z3+uZ+v=0,(Regular  cusp  catastrophe),(13)-Z3+uZ+v=0,(Dual  cusp  catastrophe).

When 2V/Z2=0, (10) and (11) can be expressed, respectively, as follows: (14)3Z2+u=0,(Regular  cusp  catastrophe),(15)-3Z2+u=0,(Dual  cusp  catastrophe).

Through (12) and (14) or (13) and (15), bifurcation set equation is shown as (16)Δ=4u3+27v2=0.

The RCCD failure criterion is defined as follows.

If Δ=4u3+27v2>0, the RCCD is in the stable state.

If Δ=4u3+27v2=0, the RCCD is in the critical state between stability and failure.

If Δ=4u3+27v2<0, the RCCD is in the failure state.

3. Determination of the Most Dangerous Sliding Plane of RCCD Failure

In this paper, the LHS method is adopted to generate the combinations of various mechanical parameters, and the elastoplastic finite element analysis is carried out based on the entropy catastrophe criterion; in the end, the most dangerous sliding plane is obtained by the sensitivity analysis of the mechanical parameters using the auxiliary analysis of PLSR.

3.1. Latin Hypercube Sampling

The LHS technique is a constrained sampling technique whereby the input parameter range is divided into equiprobable nonoverlapping intervals. Let N denote the required number of realizations and K the number of random variables. An N×K matrix P, in which each of the K columns is a random permutation of 1,,N, and an N×K matrix R of independent random numbers from the uniform (0,1) distribution are established. These matrices form the basic sampling plan represented by the matrix S as (17)S=P-RN.

Each element of S, sij, is then mapped according to its target marginal distribution as (18)x^ij=Fxj-1(sij), where Fxj-1 represents the inverse of the target cumulative distribution function for variable j. A vector X^i=[x^i1,x^i2,,x^ik] now contains input data for one deterministic computation.

3.2. Partial Least Squares Regression and Auxiliary Analysis 3.2.1. Partial Least Squares Regression

The basic steps of PLSR are given as follows .

Step 1.

E 0 and F0 are standardized matrix of independent variable X and dependent variable y, respectively.

Step 2.

The first principal component t1 is extracted from E0: t1=E0w1, w1=E0TF0/E0TF0, and w1=1. Two regression equations are solved as follows: (19)E0=t1p1T+E1,F0=t1r1T+F1, where E1 and F1 are residual matrix, and regression coefficient can be obtained as follows: (20)p1=E0Tt1t12,r1=F0Tt1t12.

Step 3.

E 0 and F0 are replaced by E1 and F1, respectively, and the same method is used to calculate w2 and t2.

It is not necessary to use all principal components to establish PLSR model, but a truncated way is adopted, which uses frontal principal components, t1,t2,,tm, to establish good prediction model. m can usually be obtained by cross validation (CV) .

3.2.2. Auxiliary Analysis

The principal component th(h=1,2,,m) is obtained using the PLSR method. On the one hand, the th as much as possible represents the information of the independent variable X={x1,,xp}; on the other hand, it explains the dependent variable y to the greatest extent. In fact, the explanation ability of the th is equivalent to the complex correlation coefficient R2 of linear regression equation on th. The explanation abilities of the th for the xj(j=1,2,,p), X, and y are expressed, respectively, as follows: (21)Rd(xj;th)=R2(xj,th),Rd(X;th)=1pi=1pRd(xj;th),Rd(y;th)=R2(y,th).

The ability that the factor xj explains the variable y is expressed by the variable projection importance index VIPj(j=1,2,,p): (22)VIPj=h=1mRd(y;th)whj2h=1mRd(y;th), where whj is the jth component of the variable wh.

3.3. The Steps of Determining the Most Dangerous Sliding Plane of the RCCD

The steps of determining the most dangerous sliding plane of the RCCD are shown as follows.

Step 1.

The deformation and strength parameters of the RCCD are taken as the effect factors, and the ranges of the factors are determined.

Step 2.

The combinations of the effect factors are established by the LHS.

Step 3.

The elastoplastic finite element model is established, and the safety coefficients of the dam failure corresponding to the combinations in Step 2 are calculated using RCCD failure criterion.

Step 4.

The sensitivities of mechanical parameters to the whole stability of the dam are determined using the auxiliary analysis of PLSR to analyze the combinations of the effect factors and corresponding safety coefficients of the dam failure. The interface of the most sensitive parameter is the most dangerous sliding plane.

4. Evaluation Standard Analysis of RCCD Safety Coefficient 4.1. Failure Process Analysis of RCCD Interface

Compared with conventional concrete dam, RCCD has obvious rolled characteristic and relatively weak construction interface. A lot of studies show that at the time of the RCCD failure, an interface is yielded completely, and the sliding channel appears . Figure 1 shows typical failure process of RCCD interface. Y-axis is strength reserve coefficient Ks and X-axis is the interface yield rate PL.

Typical failure process of RCCD interface.

( 1 ) Stability Stage (before Point A). When Ks=1.0, in addition to local small yield in the dam, the dam is basically in linear elastic state. With the increase of the Ks, yield area has a little expansion. Only when the Ks increases to certain value (corresponding point A), the yield area begins to expand. The point A is called quasielastic state.

( 2 ) Yield Deformation Stage (between Point A and Point C). With the increase of the Ks from the point A, the yield area increases rapidly to certain scope (about 30%~50%). Near the point B, the extension temporarily becomes slow until Point C.

( 3 ) Failure Stage (between Point C and Point D). The yield area expands rapidly from Point C. The upstream and downstream yield areas of the interface sharply joint together. Near Point D, sliding channel is formed and whole sliding failure of the dam leads to the loss of the ability of dam operation.

4.2. The Determination of the Evaluation Standard for RCCD Safety Coefficient

In the process of the RCCD failure, the quasielastic state is an important critical point. Before this point, the dam is mainly in elastic state and it is stable as a whole. After this point, the yield area rapidly expands until the dam failure. Therefore, the evaluation standard for RCCD safety coefficient based on the SRF method is given as (23)Ke[Ke], where Ke is the safety coefficient in the quasielastic state and [Ke] is the allowable value of the Ke.

If the allowable value [Ke]=1.0, the dam is very safe in the design condition. But in the process of dam operation, water level overload, work environment worsening, and the variation of mechanical parameters and other adverse factors will appear. Therefore, these factors must be considered into the impact of the [Ke]. Firstly, according to the results of the test and statistics, the ranges of the factors are obtained. Secondly, when the above factors change in the ranges at the same time, the [Ke] is obtained in the most adverse circumstance where the minimum of the quasielastic strength reserve coefficient appears. Specific expression is given as follows: (24)[Ke]=1.0×(1+|Kd-Km|Kd), where Kd is the quasielastic strength reserve coefficient in the design condition and Km is the minimum of the quasielastic strength reserve coefficient with the variation of various factors.

Before determining [Ke], the Km must be determined firstly. The stability of RCCD is affected by many factors, which mainly include the strength parameters (internal friction coefficient and cohesion), deformation parameters (elastic modulus), water level, and seepage condition. Among them, the deteriorations of strength parameters and seepage condition and the high water level will make the quasielastic strength reserve coefficient of RCCD become smaller. The effect of the elastic modulus is uncertain to the quasielastic strength reserve coefficient. In this paper, the LHS, LSSVM, and PSO are introduced to solve jointly the [Ke].

4.2.1. Regression Algorithm of Least Squares Support Vector Machine

Least squares support vector machine (LSSVM) is a modified version of standard SVM, in which analytical solutions can be obtained solving linear equations instead of a quadratic programming (QP) problem. Given a training data set {xi,yi}i=1N, the objective function of optimization problem is defined as [19, 21] (25)minw,b,e12w2+γ12i=1Nei2, where w is the vector of weights, γ is a regularization parameter, and ei is the random error.

The SVM inequality constraint is translated into equality constraint as follows: (26)yi=wφ(xi)+b+ei,i=1,2,,N, where w is a weight vector, b is a bias term, and φ is a nonlinear function that map the input pattern xi into a higher-dimensional feature space.

The solution of LSSVM regression will be obtained by constructing Lagrangian function as follows: (27)L(w,b,e,a)=12w2+γ12i=1Nei2-i=1Nai{wφ(xi)+b+ei-yi}, where ai is Lagrangian multiplier. The conditions for optimality are given by (28)Lw=0,Lb=0,Lei=0,Lai=0.

The solution of (28) is obtained as follows: (29)w=i=1Naiφ(xi),i=1Nai=0,ai=γei,w·φ(xi)+b+ei-yi=0.

The Elimination of w and e will yield a linear system instead of a quadratic programming problem: (30)[01NT1NTZZT+γ-1IN][ba]=[0y], where y=[y1,y2,,yN]T, 1N=[1,,1]T, a=[a1,a2,,aN]T, Z=[φ(x1),φ(x2),,φ(xN)]T, and IN is N×N identity matrix.

LSSVM function can be expressed as follows: (31)f(x)=i=1Nai(φ(x)φ(xi))+b.

In order to avoid “dimension disaster” in high-dimensional feature space, kernel k(x,xi) of input space is applied to the inner product of high-dimensional space, which solves cleverly the problem of high-dimensional calculation. Nonlinear function can be finally obtained as follows: (32)f(x)=i=1laik(x,xi)+b.

The kernel function has the following choices: polynomial function, radial basis function (RBF), and sigmoid function. In this paper, RBF is selected as follows: (33)k(x,xj)=exp{-x-xj22σ2}, where σ is a constant, which determines the scaling of the inputs in the RBF.

Main parameters of LSSVM using RBF are regularization parameter γ and kernel function width σ. To a great extent, the two parameters determine learning and generalization ability of LSSVM. Grid search algorithm is usually applied to determine these two parameters. In order to overcome the computational complexity of this method, improved grid search algorithm is applied to determine the parameters in this paper. That is, firstly rough search is applied to determine the best combination using a big step size, and then more precise search is used around this combination by small step size. For each grid point, the mean square error (MSE) from K-fold cross validation is determined, and minimum MSE is detected. In K-fold cross validation, a single subsample is used as the testing data for the LSSVM model, and the remaining subsamples are used as training data. The cross validation process is run repeatedly, with each subsample used exactly once as the testing data.

Specific steps are shown as follows.

Step 1.

Determine the ranges of the regularization parameter γ and kernel function width σ: γ[2-10,2-18,,210], σ[2-10,2-18,,28].

Step 2.

Normalize training set and use the cross validation of 1-fold to obtain the the best combination of the parameters without affecting model accuracy.

Step 3.

Enlarge the ranges of the parameters obtained in Step 2 and use finer grid search to obtain the highest precision combination of parameters by step size 20.5.

Compared to general grid search method, improved grid search method not only can reduce training volume, but also can improve forecast accuracy of the model.

4.2.2. Particle Swarm Optimization

PSO is a heuristic global optimization algorithm and has been broadly applied in optimization problems. PSO is developed on a very simple theoretical framework and can be implemented easily with only primitive mathematical operators . In PSO, a group of particles are composed by m particles in D dimension space, where the position of the particle i is Xi=(xi1,xi2,,xiD) and the speed is Vi=(vi1,vi2,,viD). The speed and position of each particle are changed in accordance with the following equation: (34)vidj+1=wividj+c1r1(pidj-xidj)+c2r2(pgdj-xidj),xidj+1=xidj+vidj+1, where i=1,2,,m, d=1,2,,D, and m is particle size. pidj is the d dimension component of the pbest that is individual optimal location of the particle i in the jth iteration. pgdj is the d dimension component of the gbest that is optimal position of all particles in the jth iteration. wi is the inertia weight coefficient. c1 and c2 are learning factors, which are commonly chosen to be 2 and will also be used in this study. r1 and r2 are random numbers in the range [0,1].

The inertia weight w is critical for the performance of PSO, which balances global and local exploitation abilities of the swarm. The most well-known algorithm for controlling inertia weight is linearly decreasing inertia weight method . The strategy of linearly decreasing weight is most commonly used, and it can improve the performance of PSO to some extent, but there are still some problems. Therefore, some adaptive algorithms for tuning inertia weight were presented [24, 25]. In this work, we propose nonlinearly decreasing inertia weight method of tuning the value of w for further performance improvement as follows: (35)w=wmax-(wmax-wmin)×(t-1)2(tmax-1)2, where wmax and    wmin are the maximum and minimum of w, respectively, t is the current iteration number, and tmax is the maximum iteration number.

4.2.3. Basic Steps

The basic steps of determining the [Ke] are shown as follows.

Step 1.

The most dangerous sliding plane of the RCCD is determined.

Step 2.

The elastoplastic model is established, and the quasielastic strength reserve coefficient Kd in the design condition is obtained using the SRF method.

Step 3.

According to the test result, the ranges of all factors are determined. On this basis, as for the factors which have an unambiguous effect on the dam stability, their values are the maximums or minimums. As for ambiguous factors, LHS is used to construct the combinations of various factors.

Step 4.

The quasielastic strength reserve coefficients corresponding to the combinations constructed in Step 3 are calculated using the SRF method.

Step 5.

The combinations of various factors and the corresponding quasielastic strength reserve coefficients in Step 4 are trained using LSSVM.

Step 6.

The LSSVM model established in Step 5 is optimized using PSO, and the corresponding minimum of quasielastic strength reserve coefficient Km is obtained.

Step 7.

According to (24), the [Ke] is obtained.

5. Example Analysis 5.1. Establishment of Finite Element Model

The Longtan dam, the highest RCCD in the world currently, was built on the Hongshui River in China. Longtan RCCD was constructed in two phases. In the first phase, normal water level is 375 m, and the elevation of dam top is 382 m. In the second phase, normal water level is 400 m, and the elevation of dam top is 406.5 m. In this paper, 11# dam section of right bank is selected in the first phase. The material partitions of Longtan 11# dam section are shown in Figure 2. According to the material partitions, RCC1, RCC2, RCC3, and foundation interfaces are looked at as main interfaces. The FEM grid of Longtan 11# dam section is shown in Figure 3.

Material partitions of Longtan 11# dam section.

The FEM grid of Longtan 11# dam section.

5.2. Determination of the Most Dangerous Sliding Plane of the Dam Failure

In this paper, RCC is regarded as transversely isotropic material. According to the material partitions shown in Figure 2, there are eight groups of parameters that need to be analyzed. The horizontal and vertical elastic modulus of RCC1, RCC2, and RCC3 interfaces are Eli and Evi  (i=1,2,3), respectively, the internal friction coefficient and cohesion of RCC1, RCC2, and RCC3 interfaces are fi and ci(i=1,2,3), respectively, the horizontal and vertical elastic modulus of foundation interface are Elj and Evj(j=4), respectively, and the internal friction coefficient and cohesion of foundation interface are fj and cj(j=4), respectively. According to the test result, the ranges of the mechanical parameters of the interfaces are shown in Tables 1 and 2. The design values of main mechanical parameters are shown in Tables 3 and 4, where Elk and Evk(k=5,6,7) are the horizontal and vertical elastic modulus of RCC1, RCC2, and RCC3 bodies, respectively, fk and ck(k=5,6,7) are the internal friction coefficient and cohesion of RCC1, RCC2, and RCC3 bodies, respectively, E8 and E9 are the elastic modulus of the conventional concrete and dam foundation rock, respectively, and f8 and c8 are the internal friction coefficient and cohesion of the conventional concrete and dam foundation rock, respectively.

Ranges of interface elastic modulus.

Elastic modulus (GPa) E l 1 E v 1 E l 2 E v 2 E l 3 E v 3 E l 4 E v 4
Minimum 12.19 9.75 10.56 8.45 9.75 7.80 12.20 9.75
Maximum 25.50 20.40 22.10 17.68 20.40 16.32 25.50 20.40

Ranges of interface strength parameters.

Strength parameter f 1 c 1 (MPa) f 2 c 2 (MPa) f 3 c 3 (MPa) f 4 c 4 (MPa)
Minimum 0.82 1.17 0.73 1.04 0.70 0.66 0.86 0.83
Maximum 1.26 2.04 1.12 1.80 1.08 1.14 1.32 1.44

Sample and calculation results of mechanical parameters.

Number E l 1 (GPa) E v 1 (GPa) E l 2 (GPa) E v 2 (GPa) E l 3 (GPa) E v 3 (GPa) E l 4 (GPa) E v 4 (GPa) f 1 c 1 (MPa) f 2 c 2 (MPa) f 3 c 3 (MPa) f 4 c 4 (MPa) K u
1 12.41 10.28 12.29 10.14 11.70 9.36 17.30 14.90 1.08 1.71 0.98 1.56 0.97 1.08 1.30 1.41 1.82
2 12.86 10.64 14.6 11.99 14.54 12.20 23.95 20.22 0.84 1.33 0.84 1.33 0.9 0.94 1.07 1.21 1.57
3 13.30 11.70 16.14 13.83 17.03 14.47 15.53 14.54 1.11 1.36 0.88 1.08 0.74 0.96 0.94 1.19 1.55
4 13.74 12.41 18.06 15.68 19.51 7.94 21.73 19.87 0.9 1.56 0.96 1.58 1.02 1.02 0.88 1.33 1.72
5 14.63 14.54 11.91 11.68 15.61 14.76 25.28 13.83 1.21 1.71 0.93 1.61 0.95 1.04 1.14 1.37 1.85
6 19.07 9.93 17.68 9.22 16.67 8.23 21.73 17.38 1.16 1.88 1.04 1.71 0.86 1.12 1.04 1.41 1.89
7 19.95 10.28 19.22 11.07 18.80 11.07 13.75 11.70 0.96 1.53 0.81 1.43 0.71 0.89 1.07 1.10 1.61
8 21.73 13.48 14.6 8.60 17.74 11.35 23.06 11.35 1.05 1.94 0.96 1.20 1.05 0.92 1.14 1.12 1.72
9 22.62 20.22 20.75 16.30 18.45 14.19 20.85 15.25 1.02 1.5 0.85 1.28 0.80 0.81 0.88 1.02 1.65
10 23.66 17.03 12.29 8.60 19.16 13.91 24.39 15.96 0.93 2.03 1.04 1.51 0.87 0.96 1.05 1.15 1.73
11 23.50 18.45 16.14 12.30 13.12 9.65 22.18 15.61 0.97 1.29 0.74 1.76 1.04 0.86 0.99 1.06 1.79
12 25.28 19.51 18.45 14.14 15.61 11.92 14.20 9.93 1.24 1.85 1.01 1.53 0.92 0.78 0.94 0.86 1.76

The ranges of the elastic modulus of the RCC bodies, conventional concrete, and dam foundation rock.

Elastic modulus (GPa) E l 5 E v 5 E l 6 E v 6 E l 7 E v 7 E 8 E 9
minimum 24.38 19.50 23.44 18.75 22.50 18.00 21.38 11.20
maximum 40.95 32.76 39.38 31.50 37.80 30.24 35.91 20.80

Using LHS to generate the combinations of the interface mechanical parameters and considering the horizontal elastic modulus bigger than the vertical, the sample results of the mechanical parameters are obtained in Table 3. Other mechanical parameters adopt design values.

Outer loads are mainly the weight load, water load, and uplift pressure. In this example, the upstream water level is the normal water level 375.00 m, and the downstream water level is 225.00 m. The strength reserve coefficients Ku of the dam failure corresponding to the combinations of the mechanical parameters are obtained using the entropy catastrophe criterion. The final results are shown in Table 3.

On this basis, using the auxiliary analysis of PLSR, the important indexes of the mechanical parameters are obtained, and the sensibilities of the mechanical parameters to the Ku are shown in Figure 4. In Figure 4, the X-axis is the mechanical parameter, and the Y-axis is the important index. The X-axis factors 1–16 correspond to the mechanical parameters in Table 3, respectively, where factors 1–8 are elastic modulus and factors 9–16 are strength parameters. Through Figure 4, it can be seen that the strength parameters have larger impact on the Ku than the elastic modulus. The 12th parameter, RCC2 interface cohesion, has the greatest impact on the Ku, which indicates that the RCC2 interface is the most dangerous sliding plane of the dam failure.

The sensibilities of the mechanical parameters to the strength reserve coefficient of the dam failure.

5.3. Analysis of Evaluation Standard for RCCD Safety Coefficient

Through the above analysis, it is shown that the RCC2 interface is the most dangerous sliding plane of the RCCD. The quasielastic strength reserve coefficient of the RCC2 interface Kd=1.76   is obtained using SRF method in the design condition.

According to the previous analysis, the deterioration of strength parameters and seepage condition and the high water level will make the quasielastic strength reserve coefficient become smaller. So the strength parameters are chosen to be the minimums, the seepage condition is without drainage in dam foundation, the uplift pressure reduction factor α is 0.3, and the upstream water level is 1.07 times higher than the design value. According to the test result, the ranges of the elastic modulus of the RCC bodies, conventional concrete, and dam foundation rock are shown in Table 4.

Adopting the LHS and SRF method, the quasielastic strength reserve coefficients corresponding to the combinations of the above deformation parameters are calculated, and the final results are shown in Table 5.

Sample and calculation results of mechanical parameters.

Number E l 1 (GPa) E v 1 (GPa) E l 2 (GPa) E v 2 (GPa) E l 3 (GPa) E v 3 (GPa) E l 4 (GPa) E v 4 (GPa) E l 5 (GPa) E v 5 (GPa) E l 6 (GPa) E v 6 (GPa) E l 7 (GPa) E v 7 (GPa) E 8 (GPa) E 9 (GPa) [Ke]
1 23.50 10.71 14.37 9.10 19.23 10.44 23.50 16.46 31.84 22.02 25.19 23.98 30.61 18.86 30.53 14.56 1.22
2 21.91 11.77 18.98 10.39 20.08 11.97 20.58 12.20 27.86 24.14 33.16 29.08 28.77 19.10 35.76 20.13 1.25
3 18.98 14.12 13.91 8.91 15.18 13.34 19.51 16.03 38.80 24.41 29.66 28.31 23.27 18.61 22.11 18.78 1.24
4 17.11 11.13 21.06 11.50 19.87 10.61 24.83 17.31 35.81 31.83 37.63 28.57 35.20 24.73 23.27 12.83 1.26
5 14.72 10.07 12.06 11.31 18.59 13.51 25.10 13.48 35.15 23.88 26.79 23.47 27.24 24.00 26.47 19.94 1.26
6 24.30 14.97 18.75 13.53 17.95 10.78 25.37 15.39 34.82 30.24 39.22 22.19 30.30 30.12 27.34 13.98 1.27
7 20.04 17.10 11.60 10.20 11.77 8.57 17.12 16.67 39.46 30.77 36.67 27.04 30.92 28.89 28.79 19.36 1.25
8 16.58 13.05 17.14 10.57 15.61 14.02 18.18 17.10 39.13 21.22 34.12 22.45 33.67 26.69 35.18 18.21 1.25
9 13.65 11.35 18.52 16.66 18.38 13.85 13.93 12.63 37.80 23.35 29.34 21.17 36.12 27.18 30.82 11.30 1.28
10 16.85 12.41 20.14 14.45 10.71 7.89 22.71 14.97 25.21 29.18 28.38 21.94 34.28 21.79 29.66 14.18 1.27
11 20.58 14.54 11.83 11.13 16.89 12.83 16.06 11.56 32.17 29.45 28.06 20.66 35.81 21.55 27.92 12.64 1.24
12 23.77 13.90 15.52 13.71 14.54 10.10 14.99 14.54 40.45 24.94 31.57 18.88 33.98 19.59 21.53 20.70 1.28

On this basis, the nonlinear relation between the deformation parameters and the quasielastic strength reverse coefficients is determined using the LSSVM. For the LSSVM, the ranges of parameters  γ and σ are [2−2, 260] and [2−2, 240], respectively. LSSVM uses improved grid search algorithm to search the global best fitness value. The optimal solutions of parameters γ and σ obtained by rough search are 254 and 238, respectively. Then, a finer search is run in the region [252, 256] and [236, 240] for γ and σ, respectively. The optimal solutions of the parameters  γ and σ of the LSSVM using improved grid search algorithm are 253 and 239, respectively.

The stepwise regression (SWR) and LSSVM performances for training samples are shown in Table 6 and Figure 5. Table 6 shows the mean squares error (MSE) from PLSR and LSSVM for the training samples. Training result of the algorithms can be clearly seen in Figure 5. Table 6 and Figure 5 explain that the training performance of LSSVM is better than SWR.

Comparison between target values and estimated values from various algorithms for the training samples.

Criterion MSE
SWR 3.3 5 e - 6
LSSVM 3.3 3 e - 15

Comparison figure of training result of the algorithms.

PSO and Genetic Algorithm (GA) are applied to optimize the established LSSVM model. For the PSO, c1 and c2 all are chosen to be 2, and the inertia weight w decreases nonlinearly from 0.9 to 0.4. For the GA, crossover probability and mutation probability are chosen to be 0.9 and 0.1, respectively. Each algorithm uses 20 particles (or populations) to search the global best fitness value with each run executed for 200 iterations.

Figure 6 depicts the great success of optimization process using PSO compared with Genetic Algorithm (GA) for the quasielastic strength reserve coefficients Km. It confirms the superiority of PSO in terms of convergence speed without the premature convergence problem. The minimums of the quasielastic strength reserve coefficients Km=1.22   are obtained using two algorithms.

Iteration process.

According to (24), the allowable value of the quasielastic strength reserve coefficients [Ke]=1.31 can be obtained. Due to Kd=1.76>[Ke]=1.31, it can be considered that the dam safety coefficient based on SRF method meets the requirement and the dam is safe in the design condition.

6. Conclusions

Evaluation standard for RCCD safety coefficient based on FEM is a problem requiring urgent solution. In this paper, the D-P criterion, SRF method, information entropy, and catastrophe theory are applied to RCCD failure criterion. On this basis, the LHS and auxiliary analysis of PLSR are used to determine the most dangerous sliding plane of RCCD failure. Finally, LHS, LSSVM, and PSO are adopted to establish evaluation standard for RCCD safety coefficient based on SRF method. The proposed method combines the advantages of several excellent algorithms and can be popularized to any RCCD. The analysis of evaluation standard for Longtan RCCD safety coefficient based on SRF method is carried out. The calculation shows that the ultimate failure is controlled by the strength of RCC2 interface. In the design condition, the quasielastic strength reserve coefficient is 1.76 and the allowable value of the quasielastic strength reserve coefficient is 1.31. Due to 1.76>1.31, the Longtan RCCD is safe in the design condition.

Conflict of Interests

The authors do not have any conflict of interests regarding the content of the paper.

Acknowledgment

The present work is supported by the Special Fund Projects supported by Basic Scientific Research Operating Expenses of Central-Level Public Academies and Institutes (no. CKSF2015022/GC, CKSF2014037/GC).

Brian F. RCC-New developments and innovations Proceedings of the Brazilian International Roller Compacted Concrete (RCC) Symposium 2008 Salvador, Brazil 1 18 Ministry of Water Resources of the People's Republic of China Design Specification for Concrete Gravity Dams (SL-319-2005) 2005 Beijing, China China Water Power Press Qin S. Q. Primary discussion on formation mechanism of dissipative structure in instability process of rock mass Chinese Journal of Rock Mechanics and Engineering 2000 19 3 265 269 2-s2.0-0033879909 McKay M. D. Beckman R. J. Conover W. J. A comparison of three methods for selecting values of input variables in the analysis of output from a computer code Technometrics 1979 21 2 239 245 10.2307/1268522 MR533252 2-s2.0-0018468345 Olsson A. M. J. Sandberg G. E. Latin hypercube sampling for stochastic finite element analysis Journal of Engineering Mechanics 2002 128 1 121 125 10.1061/(ASCE)0733-9399(2002)128:1(121) 2-s2.0-0036172940 Olsson A. Sandberg G. Dahlblom O. On Latin hypercube sampling for structural reliability analysis Structural Safety 2003 25 1 47 68 10.1016/S0167-4730(02)00039-5 2-s2.0-0037212955 Lindgren F. Geladi P. Wold S. The kernel algorithm for PLS Journal of Chemometrics 1993 7 1 45 59 Tenenhaus M. Vinzi V. E. Chatelin Y. Lauro C. PLS path modeling Computational Statistics and Data Analysis 2005 48 1 159 205 10.1016/j.csda.2004.03.005 2-s2.0-10144243373 Suykens J. A. K. Vandewalle J. Least squares support vector machine classifiers Neural Processing Letters 1999 9 3 293 300 10.1023/A:1018628609742 2-s2.0-0032638628 Kennedy J. Eberhart R. Particle swarm optimization Proceedings of the IEEE International Conference on Neural Networks December 1995 Perth, Australia 1942 1948 2-s2.0-0029535737 Hasni A. Taibi R. Draoui B. Boulard T. Optimization of greenhouse climate model parameters using particle swarm optimization and genetic algorithms Energy Procedia 2011 6 371 380 10.1016/j.egypro.2011.05.043 Drucker D. C. Prager W. Soil mechanics and plastic analysis or limit design Quarterly of Applied Mathematics 1952 10 157 165 MR0048291 ZBLl0047.43202 Wei Z. Xiaolin C. Chuangbing Z. Xinghong L. Failure analysis of high-concrete gravity dam based on strength reserve factor method Computers and Geotechnics 2008 35 4 627 636 10.1016/j.compgeo.2007.10.005 2-s2.0-46749143318 Zienkiewicz O. C. Pande G. N. Time-dependent multilaminate model of rocks-a numerical study of deformation and failure of rock masses International Journal for Numerical and Analytical Methods in Geomechanics 1977 1 3 219 247 10.1002/nag.1610010302 2-s2.0-0017515476 Shannon C. E. A mathematical theory of communication The Bell System Technical Journal 1948 27 379 423 10.1002/j.1538-7305.1948.tb01338.x MR0026286 Zhou C. Zhang Y. Research on entropy catastrophic regularity and failure criterion in the deformation and failure process of rocks Rock and Soil Mechanics 2007 28 12 2506 2510 2-s2.0-37749013467 Fu C. Chen S. Study on instability criteria of surrounding rock of underground engineering cavern based on catastrophe theory Rock and Soil Mechanics 2008 29 1 167 172 2-s2.0-39749130362 Gu C. S. Li B. O. Xu G. L. Yu H. Back analysis of mechanical parameters of roller compacted concrete dam Science China Technological Sciences 2010 53 3 848 853 10.1007/s11431-010-0053-0 2-s2.0-77950823585 Li B. Gu C. Li Z. Zhang Z. Monitoring model for dam seepage based on partial least-squares regression and partial least square support vector machine Shuili Xuebao 2008 39 12 1390 1400 2-s2.0-58249139969 Li B. Gu C. S. Wu J. K. Gradual change law of elastic mechanical parameters of roller compacted concrete dam Shuili Xuebao 2013 44 12 1488 1497 Ren X. H. Xia S. Y. Zhang S. Q. Preliminarily study on failure mechanism of RCC gravity dams Journal of Hohai University 1997 25 1 73 80 Cheng Y. M. Li L. Chi S.-C. Wei W. B. Particle swarm optimization algorithm for the location of the critical non-circular failure surface in two-dimensional slope stability analysis Computers and Geotechnics 2007 34 2 92 103 10.1016/j.compgeo.2006.10.012 2-s2.0-33947148210 Liu X. Shao C. Ma H. Liu R. Optimal earth pressure balance control for shield tunneling based on LS-SVM and PSO Automation in Construction 2011 20 4 321 327 10.1016/j.autcon.2010.11.002 2-s2.0-79954600691 Wu X. Huang Q. Zhu X. J. Thermal modeling of a solid oxide fuel cell and micro gas turbine hybrid power system based on modified LS-SVM International Journal of Hydrogen Energy 2011 36 1 885 892 10.1016/j.ijhydene.2010.08.022 2-s2.0-79251623718 Alfi A. Modares H. System identification and control using adaptive particle swarm optimization Applied Mathematical Modelling 2011 35 3 1210 1221 10.1016/j.apm.2010.08.008 MR2740862 2-s2.0-78549243263