The Scientific World Journal

Volume 2014 (2014), Article ID 483243, 14 pages

http://dx.doi.org/10.1155/2014/483243

Research Article

Ant Colony Optimization Analysis on Overall Stability of High Arch Dam Basis of Field Monitoring

Peng Lin,¹ Xiaoli Liu,¹ Hong-Xin Chen,² and Jinxie Kim¹

¹State Key Laboratory of Hydroscience and Engineering, Tsinghua University, Beijing 100084, China

²Department of Civil and Environmental Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Hong Kong

Received 30 December 2013; Accepted 16 May 2014; Published 17 June 2014

Academic Editor: Fei Kang

Copyright © 2014 Peng Lin 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 dam ant colony optimization (D-ACO) analysis of the overall stability of high arch dams on complicated foundations is presented in this paper. A modified ant colony optimization (ACO) model is proposed for obtaining dam concrete and rock mechanical parameters. A typical dam parameter feedback problem is proposed for nonlinear back-analysis numerical model based on field monitoring deformation and ACO. The basic principle of the proposed model is the establishment of the objective function of optimizing real concrete and rock mechanical parameter. The feedback analysis is then implemented with a modified ant colony algorithm. The algorithm performance is satisfactory, and the accuracy is verified. The groups of feedback parameters, used to run a nonlinear FEM code, and the displacement and stress distribution are discussed. A feedback analysis of the deformation of the Lijiaxia arch dam and based on the modified ant colony optimization method is also conducted. By considering various material parameters obtained using different analysis methods, comparative analyses were conducted on dam displacements, stress distribution characteristics, and overall dam stability. The comparison results show that the proposal model can effectively solve for feedback multiple parameters of dam concrete and rock material and basically satisfy assessment requirements for geotechnical structural engineering discipline.

1. Introduction

Scale models [1, 2], intelligent cracking control [3–5], numerical simulation [6–8], and prototype monitoring structures [9–11] have been used to investigate the behavior of concrete arch dams and foundation for the last 80 years. The purpose of these investigations is to understand the behavior of a high arch dam and their foundations by studying the external and internal actions during operation period. Failures of large dams, such as at Malpasset, Vaiont, Teton, and Kolnbrein [12–15], can be highly catastrophic, posing a great danger to people and properties. Hence, much attention should be paid to long-term safe operation of dams based on prototype surveillance and inverse analysis.

Parameters such as elastic modulus, unit weight, Poisson’s ratio, friction coefficient, and cohesion are parameters in structural analysis intrinsic to the determination of stress distributions and displacements, especially when the design of the structure is based on elasticity considerations. In a dam-foundation system, these parameters, for mass concrete, are hard to determine directly from tests due to the necessity for large specimens and large testing machines. The parameters for rock are also hard to determine because of the complicated nature of most geological situations. Currently, the inversed parameter is focused on Young's modulus of concrete and rock material, which may limit the variety of materials in the inversion Young's modulus. If we can find an algorithm with capability of inversion Young's modulus of more material, and even a variety of mechanical parameters of various materials, this will be a powerful tool for determining the mechanical parameters of dam-foundation systems. Through inverse analysis, the exact parameter values can be determined, and a precise evaluation of dam cracking mechanism, the overall stability of dams, and underground excavations can be made [16–20].

In recent years, inverse analysis has been mainly based on the two approaches of neural networks [21, 22] and optimization algorithms [23–27]. Three types of optimization algorithm have been used in feedback analysis. The first is the gradient-based direct search algorithm, such as the Levenberg-Marquardt method, conjugate gradient method, and thrust region method. The second is the relatively simple direct search method such as the simplex search method. The last type is the intelligent global search algorithm, such as genetic algorithms [24, 25], monkey algorithms [26, 27], differential evolution, particle swarm optimization, and ant colony optimization (ACO) [28, 29]. The first and the second types of algorithms both have the advantage of estimating solutions in relatively short computational times, but the results are affected by the initial values chosen and premature convergence is likely to occur. As an alternative to direct search algorithms, intelligent global search algorithms are being widely used for reverse analysis, but they are time-consuming. Ant colony optimization is one of the most widely used optimization algorithms, since it can avoid problems of premature convergence and is not influenced by the initial values [30, 31]. Ant colony optimization (ACO) is a metaheuristic algorithm for combinatorial optimization problems. The ACO algorithm was first introduced in the early 1990s [28, 32, 33] and was successfully applied to some academic problems and to real-word applications [34–36]. The ACO algorithm has been used for inverse analysis by some researchers investigating the material parameters of dams and embankment [29]. Several extensions and improvements of the original ant system (AS) algorithms have been introduced over the years. These successful ACO variants include elitist AS (EAS), rank-based AS (RAS), MAX-MIN ant system (MMAS), ant colony system (ACS), and hypercube framework (HCF). However, successful application of the ACO algorithm in the case of real large arch dams has been rather limited. The objective of this study was to develop a modified ACO algorithm for determining the mechanical material parameters of a large arch dam based on the monitored deformation data.

The paper is organized as follows. Firstly, the modified ACO algorithm for inverse analysis is proposed. Secondly, based on actual operational conditions and the monitored Lijiaxia arch dam deformation data collected over decades [37, 38], inverse analysis of the dam was carried out by employing the FEM code [39]. The optimized mechanical parameters of concrete and rock of the dam and abutments were found to be obtainable through this process. Finally, 3D numerical analysis of the entire dam was carried out using these optimized parameters to evaluate the effects of reinforcement in the dam and to investigate the crack initiation mechanism, how cracks propagated in the downstream face of the dam, and the overall safety of the dam foundation.

2. Nonlinear Feedback Numerical Model

2.1. Numerical Inverse Analysis of a Large Arch Dam

Structural health monitoring of large concrete arch dams is based on the acquisition of displacement measurements. These displacements are interpreted to identify significant deviations from what could be considered as the normal response based on statistical or deterministic models of dam behavior.

The finite element analysis method is adopted for solving the dam-foundation system. The analytical model in finite element formulation is

where is the structural stiffness matrix, is the displacement vector, and is the load vector. Given and , can be obtained. Generally, is known beforehand; depends on parameters such as Young’s modulus, cohesion, and Poisson’s ratio.

Actual deformations of arch dams can be obtained through monitoring. Given a group of concrete-rock mechanical parameters, displacements can be computed by (1). The parameters should be optimized to match the real scenario by minimizing the objective function, which is expressed as the sum of the squares of the errors between the computed displacements and the field monitored displacements. The mathematical model for inverse analysis of arch dams can be expressed as

where is the mechanical parameters vector; is the total number of observation points; is the monitored displacement at the th observed point; and is the corresponding numerical analysis displacement at the th monitoring point.

2.2. Parameter Inverse Analysis Model

This study is similar to the application of ACO in the case of travelling salesman problem (TSP) [30, 40]. Inverse analysis in the search for mechanical dam foundation parameters can be defined as follows.

There are parameters and ants. Since ACO is good at solving discretized optimization problems [36], the parameters are discretized first. Each parameter has a range and can be divided into some segments. Then the center of each segment is represented by a point. For example, the elastic modulus, , typically ranges from 1 to 50 GPa and 100 points can be selected within this range. Similarly, 100 points for each of , , and , which are Poisson’s ratio, unit weight, friction coefficient, and cohesion coefficient, can be selected. More points mean higher accuracy but with more computation time. The goal is to find a group of parameters that contains every parameter and minimizes the value of . If there is more than one kind of material involved, the number of parameters increases.

is the vector of the mechanical parameters to be analyzed, the dimension of which is . is the search space of , which is defined as

where is the th mechanical parameter of . The difference between and defines the range of the th mechanical parameter.

An artificial ant is an agent which moves between parameter points. It chooses the next point by using a probabilistic function determined by both the pheromone value, , and the heuristic value, , which is the standard deviation.

is defined as

where is a constant which ranges from 0 to 1 and is the number of parameters. Every ant has a tour. Tour corresponds to . If tour includes edge (, ), it can be described as . represents the number of ants that choose edge (, ). is the average value of all the relating to ants choosing edge (, ). It is given as

where reflects the sensitivity of the edge (, ) and must be updated. Faster convergence is prevented if more ants have a higher probability of choosing those more sensitive edges.

Without , if the number of ants choosing edge (, ), is zero, the search space is smaller and smaller. Thus, the move probability from site- to site- for ant- at time can be described by

where is a constant, which ensures that edges with zero transition probability also have a poor chance of being selected by the ants.

Additionally, can be updated as

where is the pheromone evaporation rate which is a constant parameter and ranges between zero and one. Consider

where is a constant parameter and is the tour length of ant-. The combination of , , , and should be determined to be the specific case:

where is the value corresponding to ant- and is a constant to avoid the overflow of when is too small.

2.3. Nonlinear Constitutive Model

The characteristics of high arch dams and dam foundations lead to 3D complex mesh model, highly nonhomogeneous material distributions, and very high loadings. All of these have detrimental impacts on the convergence of elastic-plastic analyses. The convergence of FEM is a principal characteristic of a stable geotechnical structure, for example, the specific and widely used strength reduction method. In this study, the back-analysis adopted a nonlinear constitutive model based on the Drucker-Prager (D-P) criterion [41, 42]. A brief introduction is as follows.

In this study, robustness of iteration and an integration policy based on the D-P criterion increased the stability of calculation. This method improved computational convergence and ensured that computation converged to the correct solution. The yielding condition for the ideal elastic-plastic model adopted the D-P criterion:

where

where , , and are major principal stresses, respectively. is the first invariant stress, is the second invariant stress, and can be obtained by the fitting Mohr-Coulomb (MC) yielding criteria. On the plane, if the D-P criterion is the circle circumscribing the Coulomb hexagon, then

where and are friction angle and cohesion coefficient of the material, respectively. If the D-P criterion is the inscribed circle to the Coulomb hexagon, then

In this study

Concrete and rock masses are materials low in tensile strength. Conditions of tension are

where is the uniaxial tension strength of the materia. If the tension strength is not given, it is evaluated according to (16) to determine its tensile strength. The program flow, first of all, estimates whether the tension condition is satisfied. If not, stresses are adjusted until the condition is satisfied. The stresses are again estimated after the adjustment.

The literature [23, 39] shows that the analytic solution for transferred stress based on the D-P criterion is equivalent to a linear prediction-radical adjustment algorithm as far the stress adjustment process is concerned. It is equivalent to the so-called closest projection algorithm from the aspect of constitutive relationships integration policy. The closest projection algorithm is of first-order accuracy and is stable unconditionally. As a special case of the generalized midpoint method, it can achieve high accuracy over large strain increments.

2.4. Modified ACO Inverse Analysis for Large Arch Dam

To verify the accuracy of the method, a group of parameters are chosen to obtain computed displacements employing the FEM code [39]. These computed displacements can be viewed as the measured ones. The inverse displacements are then compared with prototype monitored values.

In this modified ACO model, let artificial ants search such that each ant will have a tour which contain a group of parameters. An ant cycle has two half-motions. Particularly, in first half-motion, the artificial ants start at random points between and . Once the artificial ants arrive at , they then choose points to move to in tracing a path towards , which is the second half-motion, completing the ant cycle. The probability of ant choosing the next point is determined by and . The initial values of pheromone between the points are the same. The same is true of . Thus, in the first half-motion, parameter values are unlike those applying to the prototype of the ACO algorithm. Obviously, once the first half-motion is accomplished, then groups of parameters will be obtained which are used to run the FEM code to obtain numerical displacements. Then, and will be updated. Importantly, the second half-motion is not random. Only the first half-motion of the first cycle is random. The procedure above is repeated until a satisfactory result is obtained. The procedure of ACO process is presented in Figure 1. The key concept and parameters of the modified ACO model correspond with the arch dam and are descripted in Table 1.

Table 1: The key concept and parameters of the ACO model corresponding with the arch dam.

Concept
Define corresponding with the arch dam

Parameter discrete
Each group of concrete or rock mechanical parameters corresponds to a group of displacements and these displacements are defined as computed data. As the value of each parameter is continuous, some discrete points are chosen within the range.

Ant
An ant is an agent which moves between parameter discrete points. Each ant can choose a group of parameters during each search.

Objective function is expressed as the sum of the squares of the errors between the computed displacements and the field monitored displacements. is the vector of the mechanical parameters to be analyzed.

Pheromone update
The object of pheromone updating is to study the influence of material parameters on the modified ACO algorithm.

Inversion parameter
To find a group of parameters that contains every parameter and minimizes the value of .

Edge
Relationship of the discrete points between two adjacent parameters, for example, each of 50 discrete points of two adjacent parameters, there are 2,500 edges between parameters.

One cycle
Including two searches: “moveahead” and “moveback,” and output is two times parameters

One tour
A group of parameters, the tour is optimized depending on , average value of all the relating to ants choosing edge (, )

Figure 1: Flow chart of D-ACO inverse analysis for deformation of large arch dam.

3. Case Study

3.1. Brief of Lijiaxia Power Station

Lijiaxia hydropower station is located on the Yellow River, border between the villages of Jianzha and Hualong in Qinghai province, about 100 km southeast of Xining (Figure 2(a)). The dam is a concrete double curvature arch dam with maximum height, 155 m. The elevation level (EL) of the dam crest is at 2185 m and the maximum crest width is 45 m (Figure 2(b)). The total installed capacity is of 2,000 MW. The Lijiaxia arch dam has multiple functions including the provision of hydropower, flood control, and irrigation. A gravity block was installed at the left abutment, and a discharge structure, a station diversion system and a power house structure of units in two rows were constructed at the dam downstream, as shown in Figure 2(b). The dam is EL 2185 m, as stated above, and the normal reservoir water surface is EL 2180.0 m. The storage volume is 16.5 hundred million m³. Project construction began in April 1988 and the reservoir began impounding on December 26, 1996. The dam has been operating safely after about twenty years.

(a)

(b)

Figure 2: Location map and snapshot of the Lijiaxia arch dam: (a) schematic map of the location; (b) a snapshot of Lijiaxia arch dam.

The dam lies across the middle of the Lijiaxia Valley, about 5 km long, a deep, narrow, and V-shaped gorge. The valley is also symmetrical with a slope angle for two side abutments of 45~50°. The dam foundation is very complicated, and the bedrock consists of simian black mica schist and chlorite schist in continuous bands interlocked with granite rock. The joints are comparatively well developed. The metamorphic rock at the dam-heel is influenced by multistage tectonic activity, and faults and joints are well developed; see Figure 3.

Figure 3: Faults distribution in dam zone of Lijiaxia arch dam [38].

3.2. Numerical Model and Analysis Cases

On the basis of the actual Lijiaxia project sites (Figure 2), a 3D finite element model was constructed as in Figure 4(a). The 3D finite element model includes the whole arch dam (Figure 4(b)) and the foundation, which includes the main faults shown in Figure 3. The size of the 3D numerical model represents 480 × 700 × 305 m (length × width × height), which is much bigger than the dam size itself, since a vast area of the abutments is included in the model. All of the rock masses together with the related geological structures are represented in the 3D numerical model. The rock masses are of types A2, A4, A1, B2, B4, C3, and D around the dam, and under the dam are of types A2, A4, A1, B2, and B4. The geological structures in these ranges include faults f33, f20, f35, F32, F26, F27, F20, F50, and F201, as shown in Figures 3 and 4(b). The main elastic modulus parameters of the rock masses, of the geological structures (i.e., weak zones), and of the arch dam material are listed in Table 2 based on the monitoring feedback analysis. The total number of hexahedral elements was 25376 and 29529 of these were used to model the dam and foundation in Figure 4(a).

Table 2: The elastic modulus inversion results using various feedback methods (unit: GPa).

Number
Material
Design value
DPIP-ACO
Generalized least squares
Neural networks

1
Dam concrete
20
27.2
30
28

2
A2
20
23.3
28
26 (A)

3
A4
12
11.2
17

4
B2
8
10
11.6
22 (B)

5
B4
12
16.4
18

6
C3
5
7
7
13 (C)

7
A1
15
16
22.5
26 (A)

8
F20
0.62
—
—
—

9
F20-1
0.31
—
—
—

10
F26
0.52
—
—
—

11
F27
0.57
—
—
—

12
F32
0.243
—
—
—

13
F50
0.21
—
—
—

14
D
2.75
3.4
3.85
—

15
f20
0.94
—
—
—

16
f33
1.05
—
—
—

17
f35
0.8
—
—
—

18
Gravity pier
20
29.6
30
—

19
Foundation reinforcement zone
10
27.2
15
—

(a)

(b)

Figure 4: 3D FEM mesh model of Lijiaxia arch dam: (a) 3D overall mesh model; (b) main faults distribution.

During analysis, the node displacements of the overall model are applied as the boundary conditions. The upstream/downstream surfaces of foundation are employed displacement constraints along the river direction () and constraint of two abutment side surfaces, using a transverse river direction to the (). The bottom surface of foundation is taken as the vertical displacement constraints ( direction).

The numerical analysis especially takes reinforcement parameters into account. Analysis and evaluation mainly concerns dam displacements, stresses, safety failure locations in the dam model and foundation, stability evaluation of the abutments, and the riverbed interface. In this study, the dam self-weight, water and silt loadings, and temperature loadings were taken into account in the analyses. The upstream water level is EL 2178 m and the downstream level is EL 2050 m. The main analysis cases are listed below.

 Analysis case 1: dam self-weight + normal water load. Analysis case 2: dam self-weight + normal water load + silt load + temperature dropped loading (the temperature loading determined by the average March temperature). Analysis case 3: dam self-weight + normal water load + silt load + temperature increased load (the temperature loading determined by the average September temperature).

In order to prove the ant colony optimization analysis is effective in relation to the monitored feedback as far as concrete and rock mechanical parameters are concerned, a comparison analysis was carried out employing various material parameters obtained from design values, generalized least squares, and neural networks methods.

3.3. Feedback Analysis on Material Parameter

3.3.1. D-ACO Initial Value Determination

According to the proposed inverse analysis model in Section 2.2, each group of parameters corresponds to a group of displacements, and these displacements are defined as computed data. As the value of each parameter is continuous, some discrete points are chosen within the range. In the inverse analysis program, the number of discrete points for each parameter is 110. For the inverse analysis, it is better to set 1.5 as the ratio of site to ant number, and the ant number is chosen as 10. In addition, it was found that the control values are very important in the D-ACO algorithm, and it is very essential to determine them by test computing. In this case, the values of control parameters were set as , , , and .

3.3.2. D-ACO Pheromone Update

This is the updating stage of the pheromone value, , and the heuristic value, , following the first ant cycle. The object of pheromone updating is to study the influence of material parameters on the D-ACO algorithm. The heuristic algorithm parameters are , , and . According to the proposed feedback model in Section 2, three times the discrete space is constructed for solving material parameter feedback of Lijiaxia arch dam.

3.3.3. Established of the Objective Function

Based on discussion of Section 2.2, comparing site monitoring and numerical displacements, the objective function is established as shown in (2) in Section 2.1 above. Then, the minimized objective function is then determined. Table 3 shows the field monitoring and numerical displacements used for solving the objective function .

Table 3: Field monitored displacement values and numerical results for solving the objective function .

Monitoring point
EL (m)
Dam monolith
Displacement of case  3 (mm)
Displacement of case  2 (mm)

Monitoring
Numerical simulation
Monitoring
Numerical simulation

2
2150
Number 6
19.52
17.75
24.33
25.1

3
2114
Number 6
15.64
15.05
17.71
18.6

4
2087
Number 6
7.29
10.75
8.45
12.2

5
2150
Number 11
25.53
23.7
33.42
38

6
2114
Number 11
26.84
24
29.74
32

7
2087
Number 11
22.22
20
22.87
25

8
2059
Number 11
13.65
14
15.30
16

9
2035
Number 11
8.18
7.8
8.74
8.8

10
2150
Number 16
12.91
19.1
14.05
25.8

11
2114
Number 16
12.20
14.2
15.32
17.1

12
2087
Number 16
8.08
8.3
9.05
9.3

3.3.4. Parameter Feedback Results

In this study, the elastic modulus, , for thirteen material types was inversed. The feedback analysis results for the various material parameters are shown in Table 2.

The feedback results show (1) an elastic modulus of the dam concrete which is about 36% greater; (2) an elastic modulus of abutment and riverbed rock 4~85% which were not included feedback in the comparisons. Based on feedback results shown in Table 2, the D-ACO feedback parameter results are better than for the other two methods.

3.4. Feedback Analysis Results Discussion

Based on the material parameters obtained from the various methods (original design value, generalized least squares), comparative analyses were conducted on dam displacements, stress characteristics, and overall stability of dam-foundation system.

3.4.1. Comparison on Dam Displacement

Table 4 shows the main numerical displacement results under case 2. Table 5 shows the displacement differences in the various numerical results using various feedback material parameters. Based on the Tables 4 and 5 results, the largest dam displacement along the direction of the river decreased under analysis case 2. Comparing with that simulated using original design values, the difference is 9.5 mm (Table 5). The abutment displacements in the direction along the river also show decreasing values of the left and right hand abutments of 1.68 mm and 2.27 mm, respectively. These values are basically equivalent.

Table 4: Dam displacements in the direction along/perpendicular the river (analysis case  2, unit: mm).

Analysis case
EL (m)
Design value
Generalized least squares
D-ACO

Right arch side
Arch crown
Left arch side
Right arch side
Arch crown
Left arch side
Right arch side
Arch crown
Left arch side

Case  2
2185
1.5
0.7
48.8
1.3
3.4
1.4
0.9
−0.4
35.0
0.9
2.2
0.9
0.9
−0.7
39.6
0.7
2.6
1.7

2148
4.5
−2.6
48.6
1.2
5.3
2.4
3.0
−1.7
34.1
0.8
3.6
1.7
2.6
−1.8
37.9
0.6
3.9
2.5

2119
7.3
−3.3
42.8
1.0
8.0
3.0
5.0
−2.2
30.1
0.7
5.6
2.2
4.5
−2.3
33.2
0.5
5.9
2.8

2100
8.7
−3.5
37.3
0.7
9.3
3.2
6.1
−2.4
26.2
0.5
6.6
2.4
5.6
−2.6
28.7
0.3
7.0
3.0

2075
10.4
−3.3
27.7
0.5
8.5
2.6
7.3
−2.3
19.4
0.4
6.1
2.0
7.4
−3.0
21.0
0.2
6.8
2.9

2050
9.4
−2.0
17.9
0.3
8.4
2.5
6.7
−1.4
12.6
0.2
6.0
1.9
6.6
−2.0
13.4
0.1
6.5
2.6

2035
7.7
−0.9
11.4
0.2
7.3
1.8
5.6
−0.6
8.2
0.2
5.3
1.4
5.5
−1.1
8.6
0.1
5.6
1.8

2030
6.9
−0.4
9.4
0.2
6.8
1.5
5.0
−0.3
6.8
0.2
4.9
1.1
5.2
−0.9
7.2
0.1
5.3
1.4

Max. value
10.4
−3.5
49.1
1.4
9.3
3.2
7.3
−2.4
35.0
0.9
6.6
2.4
7.4
−3.0
39.6
0.7
7.0
3.0

represents displacements in the direction perpendicular the river; positive value means deformation to left abutment; represents displacements in the direction along the river; and positive value means deformation to downstream direction.

Table 5: The difference of downstream displacement in the direction perpendicular, along river obtained from various feedback analysis methods (analysis case  2, unit: mm).

Analysis case
EL (m)
Between generalized least squares and design
Between D-ACO and design

Right arch side
Arch crown
Left arch side
Right arch side
Arch crown
Left arch side

Case  2
2185
0.6
1.1
13.8
0.4
1.2
0.5
0.6
1.4
9.2
0.6
0.8
−0.3

2148
1.5
−0.9
14.5
0.4
1.7
0.7
1.9
−0.8
10.7
0.6
1.4
−0.1

2119
2.3
−1.1
12.7
0.3
2.4
0.8
2.8
−1
9.6
0.5
2.1
0.2

2100
2.6
−1.1
11.1
0.2
2.7
0.8
3.1
−0.9
8.6
0.4
2.3
0.2

2075
3.1
−1
8.3
0.1
2.4
0.6
3
−0.3
6.7
0.3
1.7
−0.3

2050
2.7
−0.6
5.3
0.1
2.4
0.6
2.8
0
4.5
0.2
1.9
−0.1

2035
2.1
−0.3
3.2
0
2
0.4
2.2
0.2
2.8
0.1
1.7
0

2030
1.9
−0.1
2.6
0
1.9
0.4
1.7
0.5
2.2
0.1
1.5
0.1

Max. value
3.1
−1.1
14.1
0.5
2.7
0.8
3
−0.5
9.5
0.7
2.3
0.2

Under the same analysis case, at the upper EL 2050 m level, in the direction perpendicular to the river, the difference in displacements between design material parameters and optimized material parameters (D-ACO) is greater on the right side than on the left abutment. Below EL 2050 m, the difference in displacement of right side is greater than that of a shift toward the left side. The results are consistent with field monitoring results. The displacements monitoring of each dam monolith showed a shift to the right bank before 1998, while in 1998, and especially after 2000, except for the foundation displacement which continued to shift toward the right bank, the displacement of the upper part of the arch crown shifted towards the left bank. The shifted displacements were 1.53 mm at EL 2150 m and 2.45 mm at EL 2185 m.

Figure 5 shows numerical and monitoring displacements of the crown cantilever, at the left and right hand span of the dam. The comparison results are consistent.

(a) Crown cantilever

(b) Right 1/4-span

(c) Left 1/4-span

Figure 5: Numerical and field survey results of dam displacement in the direction along river.

Contrasting curves of temperature rise and fall at different elevation levels for each dam monolith are shown in Figure 5. The numerical results show the following.

(1)Field monitoring results and feedback analysis results fit particularly well in the arch crown beam and right hand arch (approximately 1/4 arch, number 6 dam monolith). Monitored results and feedback analysis for the temperature drop condition (case 2) are basically the same.(2)Errors exist between field monitoring results and feedback evaluation at EL 2150 m in the left hand arch (approximately 1/4 arch results, number 16 dam monolith).(3)Generally after ten years dam prototype observations, monitoring and reinforcement numerical analysis of the dam downstream are basically the same.

3.4.2. Comparison of Dam Stresses

For analysis case , stress distributions (Figure 6) on both upstream and downstream dam faces satisfy the usual stress regularity of a double curvature arch dam. The largest compressive stresses are between 3.3~3.6 MPa on the downstream face and 1.2~2.6 MPa tensile on the upstream face. Numerical results show that the tensile stress on the upstream surface is greater than that on the downstream face. The stress distribution characteristics with various calculated material parameters are shown in Table 6.

Table 6: Characteristic value of dam stresses under analysis case  2.

Position
Stress type
Design (MPa)
Generalized least squares (MPa)
D-ACO (MPa)

Upstream surface
Maximum tensile stress
1.24 (left side)
1.35 (left dam abutment)
2.56 (left side)

Maximum compressive stress
−3.6 (EL 2100 m)
−3.67 (EL 2100 m)
−3.33 (EL 2093 m)

Downstream surface
Maximum tensile stress
0.61 (EL 2150 m)
0.69 (EL2141 m)
0.7 (EL 2151 m)

Maximum compressive stress
−7.52 (EL 2060 m)
−7.73 (EL 2060 m)
−8.12 (EL 2060 m)

Interface
Maximum tensile stress
0.61 (dam heel)
0.38 (dam heal)
—

Maximum compressive stress
−1.87 (EL 2030 m)
−1.88 (EL 2030 m)
—

(a)

(b)

Figure 6: Principal stress vector contour of upstream/downstream surface under analysis case 2: (a) upstream surface; (b) downstream surface.

The comparison analysis is as follows.

(1)Under analysis case 2, the dam downstream surface is predominantly in a compressive stress state adopting the optimal parameters evaluation. The maximum compressive stress when using D-ACO parameters is −7.22 MPa, greater than the cases of design parameters and generalized least squares.(2)Analyzing the two cases of temperature drop and temperature rise, there is a tensile zone in the same direction as the beam at the dam downstream surface. The stress level is below 1 MPa in both cases.(3)The dam is in a stable stress state. Comparing the results obtained with different numerical material parameters, the characteristic of dam stress is in agreement with that obtained by D-ACO.(4)The comparison results show that the improved ant colony algorithm (D-ACO) can effectively determine the material mechanical parameters and reflects well the actual dam deformation and dam stress distribution. The algorithm accuracy satisfies project safety evaluation requirements.

3.4.3. Comparison on Overall Safety Factors

By employing the material parameter obtained from D-ACO and generalized least squares, the overall dam safety factors are shown in Table 7. Based on comparison results, the following conclusions are illustrated.

(1)Under analysis case 2, dam foundation safety factors are basically symmetric. After reinforcement, point safety factor at the elevation EL 2150 m at the left bank is a little lower than that at the right. No overall yielding appears.(2)The riverbed always has the lowest safety factors, and cracks first occur at the dam-foundation interface. Once yielding zones and cracks occur, unbalanced thrust forces will transfer to those zones with comparatively higher safety factors at both abutments, where high bearing capacity levels develop. As the height-width ratio of Lijiaxia arch dam is high, and there is a large pedestal of thickness 30 m, the safety factor in the valley upstream is at least 2.0 in any case under normal load, and there is no cracking.(3)During dam overloading process, unbalanced forces in the arch dam will transfer to the two banks, under the water load of 3 times normal loads, the depth of the crack upstream is about 1/4 of the dam thickness, and the point safety factor is about 1.0~2.0. The carrying capacity of the riverbed will reduce and transfers to areas with higher safety factors at both banks. Under a water load of 3 times normal loads, the dam can still work and without any yielding zone at the downstream interface. Under a water load of 5 times, there is yielding at the downstream interface.(4)Based on Table 7, the overall safety factor if using D-ACO parameters is slightly lower than with the case of generalized least squares. Abutment safety factors are reduced by approximately 5% to 15%.

Table 7: Point safety factor of abutments.

EL (m)
Generalized least squares
D-ACO

Left abutment
Right abutment
Left abutment
Right abutment

2185
3.0 (gravity pier), 1.5~2.0 (out of gravity pier)
2
2.0~3.0 (gravity pier); 1.1~1.5 (out of gravity pier)
2.0

2148
1.5~3, 1.2~1.5 (fault)
2~3
1.5~3.0; 1.0~1.5 (fault)
1.5~2.0

2100
2~3
2, 1.2~1.5 (fault)
1.5~3.0
1.5~2.0; 1.2~1.5 (fault)

2050
2~3, 1.2~1.5 (fault)
2
2.0~3.0; 1.1~1.5 (fault)
1.5~2.0

2030

2~3
1.5~3.0
1.5~5.0

Above all, modified ACO showed effectively distributed computing capabilities, strong robustness, and easiness to combine with other algorithms, or FEM numerical code, and can well avoid premature convergence phenomenon. The proposal model can effectively solve for feedback multiple parameters of dam concrete and rock material. Through inverse analysis, the exact parameter values can be determined, and a precise evaluation of dam cracking and deformation mechanism and the overall stability of dam-foundation can be made.

4. Conclusions

Resulting from this study, a modified dam ant colony optimization (D-ACO) model is proposed for obtaining concrete and rock mechanical parameters of a large arch dam. Based on field monitored deformations and the ant colony optimization technique, a typical dam parameters feedback problem was solved using a nonlinear back-analysis numerical model. The basic D-ACO principle is introduction of distributed computing, D-ACO initial value determination, and D-ACO pheromone update, which is used to establish the objective function. The numerical analysis is implemented through proposed dam ant colony algorithm. The mechanical parameters are determined using this algorithm, and construct solutions combined with nonlinear constitutive relations.

By employing the proposed back analysis model, calculated deformations and a stability evaluation of the Lijiaxia arch dam were compared with the monitoring results taken over 10-year supervision period of the dam prototype. The results demonstrated that the proposed model can effectively solve for feedback multiple parameters of dam concrete and rock material. Through inverse analysis, the exact parameter values can be determined, and a precise evaluation of dam cracking and deformation mechanism and the overall stability of dam-foundation can be made.

Conflict of Interests

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

Acknowledgments

This research work was supported by the National Natural Science Foundation of China (Nos. 11272178 and 51339003), National Basic Research Program of China (973 Program) Grants Nos. 2011CB013503 and 2013CB035902, and Tsinghua University Initiative Scientific Research Program. The authors are very grateful to Professor W. Y. Zhou for supporting this study, and to Guest Editors of this special issue and three reviewers for their critical recommendations which helped the authors to improve this paper significantly.

References

P. Lin, W. Zhou, and H. Liu, “Experimental study on cracking, reinforcement, and overall stability of the Xiaowan super-high arch dam,” Rock Mechanics and Rock Engineering, 2014.

J. M. Raphael and R. H. Carlson, Measurement of Structure Action in Dams, James J.Gillick and Co, Berkeley, Calif, USA, 1956.

P. Lin, Q. B. Li, and H. Hu, “A flexible network structure for temperature monitoring of a super high arch dam,” International Journal of Distributed Sensor Networks, vol. 2012, Article ID 917849, 10 pages, 2012.

P. Lin, Q. B. Li, S. W. Zhou, and Y. Hu, “Intelligent cooling control method and system for mass concrete,” Journal of Hydraulic Engineering, vol. 44, no. 8, pp. 950–957, 2013.

P. Lin, Q. B. Li, and P. Y. Jia, “A real-time temperature data transmission approach for intelligent cooling control of mass concrete,” Mathematical Problems in Engineering, vol. 2014, Article ID 514606, 10 pages, 2014.

A. R. Ingraffea, “Case studies of simulation of fracture in concrete dams,” Engineering Fracture Mechanics, vol. 35, no. 1–3, pp. 553–564, 1990.

Z. Z. Liang, C. A. Tang, H. X. Li, T. Xu, and Y. B. Zhang, “Numerical simulation of 3-d failure process in heterogeneous rocks,” International Journal of Rock Mechanics and Mining Sciences, vol. 41, no. 1, pp. 2–6, 2004.

L. C. Li, C. A. Tang, W. C. Zhu, and Z. Z. Liang, “Numerical analysis of slope stability based on the gravity increase method,” Computers and Geotechnics, vol. 36, no. 7, pp. 1246–1258, 2009.

ICOD, Dam Monitoring: General Considerations, Bulletin 60, International Commission on Large Dams, Paris, France, 1988.

Z. Z. Liang, N. W. Xu, K. Ma, S. B. Tang, and C. A. Tang, “Microseismic monitoring and numerical simulation of rock slope failure,” International Journal of Distributed Sensor Networks, vol. 2013, Article ID 845191, 10 pages, 2013.

P. Novak, A. I. B. Moffat, and C. Nalluri, Hydraulic Structures, Unwin Hyman, London, UK, 1990.

Foundations for dams, Asilomar Conference Grounds Pacific Grove, American Society of Civil Engineers, New York, NY, USA, 1974.

G. A. Kiersch, S. Brooks, and M. Anderson, Eds., “Vaiont reservoir disaster,” Civil Engineering, vol. 34, pp. 32–39, 1964.

D. N. Petley and R. J. Allison, “The mechanics of deep-seated landslides,” Earth Surface Processes and Landforms, vol. 22, no. 8, pp. 747–758, 1997.

K. Baustadter and R. Widmann, “The behaviour of the Kolnbrein arch dam,” in Fifteenth International Congress on Large Dams, (Lausanne, Switzerland: June 1985), vol. 2, pp. 633–651, 1985.

P. Lin, R. Wang, S. Kang, H. Zhang, and W. Zhou, “Key problems of foundation failure, reinforcement and stability for superhigh arch dams,” Chinese Journal of Rock Mechanics and Engineering, vol. 30, no. 10, pp. 1945–1958, 2011.

P. Lin, Y. Zhou, H. Liu, and C. Wang, “Reinforcement design and stability analysis for large-span tailrace bifurcated tunnels with irregular geometry,” Tunnelling and Underground Space Technology, vol. 38, pp. 189–204, 2013.

L. C. Li, T. H. Yang, Z. Z. Liang, and C. A. Tang, “Numerical investigation of groundwater outbursts near faults in underground coal mines,” International Journal of Coal Geology, vol. 85, no. 3-4, pp. 276–288, 2011.

P. Lin, X. Liu, Y. Hu, W. Xu, and Q. Li, “Deformation stability analysis of Xiluodu arch dam under stress-seepage coupling condition,” Chinese Journal of Rock Mechanics and Engineering, vol. 32, no. 6, pp. 1145–1156, 2013.

H. Y. Liu, J. C. Small, and J. P. Carter, “Full 3D modelling for effects of tunnelling on existing support systems in the Sydney region,” Tunnelling and Underground Space Technology, vol. 23, no. 4, pp. 399–420, 2008.

E. M. R. Fairbairn, N. F. F. Ebecken, C. N. M. Paz, and F.-J. Ulm, “Determination of probabilistic parameters of concrete: Solving the inverse problem by using artificial neural networks,” Computers and Structures, vol. 78, no. 1, pp. 497–503, 2000.

S. J. Li and Y. X. Liu, “Parameter identification of concrete dam with neural networks,” Journal of Dalian University of Technology, vol. 40, no. 5, pp. 531–535, 2005.

P. Huang and Y. G. Meng, Optimal Theories and Methods, Tsinghua University Press, 2009.

T.-H. Yi, H.-N. Li, and M. Gu, “Optimal sensor placement for structural health monitoring based on multiple optimization strategies,” Structural Design of Tall and Special Buildings, vol. 20, no. 7, pp. 881–900, 2011.

T.-H. Yi, H.-N. Li, and M. Gu, “Optimal sensor placement for health monitoring of high-rise structure based on genetic algorithm,” Mathematical Problems in Engineering, vol. 2011, Article ID 395101, 11 pages, 2011.

T.-H. Yi, H.-N. Li, and X.-D. Zhang, “A modified monkey algorithm for optimal sensor placement in structural health monitoring,” Smart Materials and Structures, vol. 21, no. 10, Article ID 105033, 9 pages, 2012.

T.-H. Yi, H.-N. Li, and X.-D. Zhang, “Sensor placement on Canton Tower for health monitoring using asynchronous-climb monkey algorithm,” Smart Materials and Structures, vol. 21, no. 12, Article ID 125023, 12 pages, 2012.

M. Dorigo and L. M. Gambardella, “Ant colonies for the travelling salesman problem,” BioSystems, vol. 43, no. 2, pp. 73–81, 1997.

M. Dorigo and S. Thomas, Ant Colony Optimization, Massachusetts Institute of Technology, 2004.

G. Cai and E. Dong, “Comparison and analysis of generation algorithm and ant colony optimization on TSP,” Computer Engineering and Application, vol. 43, no. 10, pp. 96–98, 2007.

C. Blum and M. Dorigo, “The hyper-cube framework for ant colony optimization,” IEEE Transactions on Systems, Man, and Cybernetics, Part B, vol. 34, no. 2, pp. 1161–1172, 2004.

B. Bullnheimer, R. F. Hartl, and C. Strauß, “A new rank-based version of the ant system: a computational study,” Central European Journal for Operations Research and Economics, vol. 7, no. 1, pp. 25–38, 1999.

C. Blum, “Ant colony optimization: Introduction and recent trends,” Physics of Life Reviews, vol. 2, no. 4, pp. 353–373, 2005.

F. Kang, J. Li, and Q. Xu, “Structural inverse analysis by hybrid simplex artificial bee colony algorithms,” Computers and Structures, vol. 87, no. 13-14, pp. 861–870, 2009.

M. R. Jalali, A. Afshar, and M. A. Mariño, “Multi-colony ant algorithm for continuous multi-reservoir operation optimization problem,” Water Resources Management, vol. 21, no. 9, pp. 1429–1447, 2007.

K. Bongik and K. Junghyun, “Discrete optimization of plane frame structures using genetic algorithms,” Journal of the Korea Society for Marine Environmental Engineering, vol. 16, no. 4, pp. 25–31, 2002.

Z. Q. Jian, J. G. Bai, and Y. F. Chen, “Outlines on the optimal design of Lijiaxia Hydropower Station,” Journal of Hydro Electric Power, vol. 18, no. 3, pp. 4–12, 2002.

J. Bai, P. Lin, P. Li, and W. Zhou, “Reinforcement effects of Lijiaxia arch dam on complicated rock foundations and its back analysis,” Chinese Journal of Rock Mechanics and Engineering, vol. 27, no. 5, pp. 902–912, 2008.

W. Y. Zhou and R. Q. Yang, Nonlinear Analysis on Concrete Arch Dam By 3D FEM (TFINE), Department of Hydraulic Engineering, Tsinghua University, 1988.

P. Zou, Z. Zhou, G.-L. Chen, and J. Gu, “Multilevel reduction algorithm to TSP,” Journal of Software, vol. 14, no. 1, pp. 35–42, 2003.

J. J. López Cela, “Material identification procedure for elastoplastic Drucker-Prager model,” Journal of Engineering Mechanics, vol. 128, no. 5, pp. 586–591, 2002.

R. Wang, P. Lin, and W. Zhou, “Study on cracking and stability problems of high arch dams on complicated foundations,” Chinese Journal of Rock Mechanics and Engineering, vol. 26, no. 10, pp. 1951–1958, 2007.