Determining the Optimal Placement of Sensors on a Concrete Arch Dam Using a Quantum Genetic Algorithm

1State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Hohai University, Nanjing 210098, China 2National Engineering Research Center of Water Resources Efficient Utilization and Engineering Safety, Hohai University, Nanjing 210098, China 3College of Water-Conservancy and Hydropower, Hohai University, Nanjing 210098, China 4Zhejiang Institute of Hydraulics and Estuary, Hangzhou 310020, China 5Engineering Safety and Disaster Prevention Department, Changjiang River Scientific Research Institute, Wuhan 430010, China


Introduction
Under the effects of outburst accidents (such as earthquakes and wind) and operating loads, engineering structures can accumulate damage, which may lead to destructive accidents.In fact, internal damage of engineering structures inevitably leads to changes in structural dynamic parameters, such as the natural frequency, damping, and mode shapes.Therefore, structural modal identification has become the core technology of modern dynamic testing and on-line monitoring of complex engineering structures.
Modal identification and damage diagnosis based on structural vibrations have been widely adopted in many fields of civil engineering, but the application of this technology to hydraulic structures is still in the initial stage.In the 1960s, the construction bureau of California firstly developed the mechanical vibration machine and conducted an operational dynamic test on the Montst arch dam.The first four natural frequencies of the dam were successfully measured.In the 1990s, Houqun et al. [1] conducted prototype dynamic tests on the Dongjiang arch dam in Hunan Province and on the Longyangxia arch dam in Qinghai Province using a blasting method and extracted the corresponding modal parameters from the test results.Yifeng and Ming [2] established a weighted rubber model of a double-curvature arch dam.Jijian et al. [3] built a large-scale hydroelastic model of the Laxiwa arch dam to conduct the experimental modal analysis in order to identify the structural modal parameters.Mridha and Maity [4] investigated the nonlinear response of a concrete gravity dam-reservoir system using laboratory experiments on a small-scale model of the Koyna dam.Altunis ¸ik et al. [5] studied the variations of the modal parameters of a damaged arch dam before and after retrofitting using laboratory model experiments.Darbre et al. [6] researched the relationship between the dam natural frequency and reservoir water level using a plaster model of an arch dam.Loh and Wu [7] identified the modal parameters of the Fei-Tsui arch dam from monitoring data during a strong earthquake and studied the influence of the reservoir water level on the structural modal parameters and nonuniform input on the dynamic response of the arch dam.Mau and Wang [8] performed a system identification of an arch dam using vibrational test data.Sevim et al. [9] identified the modal parameters of the Berke arch dam using the frequency domain method based on the environmental excitations and then calibrated a three-dimensional FEM model with the results.
During modal experiments, determining the optimal number of sensors and corresponding configuration on structures has received increasing attention because inappropriate sensor placement can reduce the accuracy of the identified modal parameters.Meanwhile, modern concrete arch dams with complex structures and giant volumes may contain considerable degrees of freedom (DOFs); thus, the arrangement of sensors would particularly influence the measurement accuracy of the modal parameters [10,11].During traditional dynamic tests of dams, the sensors are arranged based on engineering experience, which does not guarantee optimal sensor placement [12][13][14].Therefore, research of the optimal sensor placement on concrete arch dams has become an important subject in recent years.
Many researchers have studied optimization criteria of sensor placement for structural modal identification in the past few years.In 1990, Kammer [15] proposed the effective independence (EFI) method to obtain the greatest spatial resolution of the targeted modes.The Fisher information matrix (FIM) is established to guarantee linear independence of interesting modal vectors and reserve the DOFs that contribute most to the independence of the targeted modes.Liu and Tasker [16] proposed the Multiple-Reference Ibrahim Time Domain method for sensor placement and developed a relationship between the sensor locations and variance of identification.By combining the Eigensystem Realization algorithm, Lim [17] arranged sensors to minimize the condition number of the Hankel matrix to maintain the independence of the targeted modes.Rafajłowicz [18] revealed a relationship between the information matrix and density of the input spectrum and measurement positions and studied the problem of sensor placement for parameter identification in the frequency domain.Xing and Bainum [19] researched the problem of optimal sensor placement based on the degree of controllability and observability of the discrete system.Reynier and Abou-Kandil [20] placed sensors by maximizing the minimum eigenvalue of the Gramian matrix.Shih et al. [21] defined the degree of controllability and observability as a second-order ordinary differential equation and deduced a relationship between the equation and frequency-response function.The sensors are arranged based on the contribution of each DOF to the index.Salama et al. [22] claimed that the sensors should be placed on the locations with great energy because it is good for modal identification and optimized the sensor locations by maximizing the modal strain energy (MSE).Baruh and Choe [23] used spline interpolation of the response of the measured points to obtain information about the unmeasured points.The sensors were optimized by minimizing the fitting error of interpolation.Breitfeld [24] arranged sensors by minimizing the off-diagonal elements of the Modal Assurance Criterion (MAC) matrix.Cruz et al. [25] established the fitness function by maximizing the natural frequency identification effectiveness and the mode shape independence and adopted the custom genetic algorithm to optimize the sensor configuration.Debnath et al. [26] evaluated the modal participation at individual degree of freedom (DOF) for the target modes and proposed the modal contribution in output energy (MCOE) as the optimization criterion.Papadimitriou [27] adopted the theory of information entropy to measure the uncertainty in the system parameters and proposed the nominal structural model to optimize the sensor configuration on the truss structure.In recent years, various optimization methods have been used in the process of sensor placement, such as the serial method [28,29], particle swarm optimization [30], simulated annealing algorithm [31], and genetic algorithm [32].
This paper studied two traditional methods for maximizing the modal information, the effective independence method based on the maximal determinant of the FIM and the Kinetic Energy method based on the maximum modal strain energy.One drawback of the EFI method is that the locations with low modal stress energy may be selected, which would result in the loss of modal information.One drawback of the Kinetic Energy method is that the optimization process is highly dependent on the partitioning of the finite element mesh.Based on the above analysis, a novel optimization method combining the EFI method with the MSE method is proposed that could effectively resolve the drawbacks of the two methods.Considering the enormous number of DOFs of the hydraulic structure, a quantum genetic algorithm (QGA) is adopted to increase the computational efficiency and accuracy in this paper.Finally, a specific concrete arch dam is used as an example, and the sensors on the upstream surface are optimized with the proposed method.A comparison of the traditional and proposed methods shows that the proposed method has a higher convergence speed and better optimization accuracy, which has theoretical and practical application values.

The Basic Principle of the Optimal Sensor Placement
2.1.The Effective Independence Method.In 1991, Kammer presented the effective independence method and introduced the Fisher information matrix, which is based on the displacement modal matrix.To make the concerned modal vectors linearly independent, the measured points are ranked based on an effective independence value of the modal matrix and the measurement point with the minimum contribution is deleted in succession.For modal experiments that only identify the structural mode shapes, the generalized coordinates  of the mode shapes can represent the identified parameters.The output value of the sensor can be expressed as follows: where   represents the output value of the sensor; Φ  is the modal matrix; and  is the vector of the modal coordinates.
The least square estimation of  can be written as follows: By introducing measurement noise, (3) is defined as where  is the Gaussian white noise of a uniform distribution with the variance of  2 .By minimizing the covariance matrix of the modal coordinates, the effective unbiased estimate of  can be obtained as follows: where  = Φ   Φ  is the Fisher information matrix (FIM).Therefore, an optimal estimation of the modal coordinates  can be obtained when the determinant of the FIM is maximized.Then, the covariance matrix  is minimized and the targeted modal vectors are linearly independent, which is good for structural modal identification.
Then, the corresponding fitness value  1 of the EFI method can be represented as follows:

The Modal Strain Energy Method. Modal strain energy (MSE
) is an index that is sensitive to the variations of the structural parameters and the ratio between the element MSE and total structural kinetic energy and is one of the system eigenvalues.Therefore, the MSE index is usually applied in structural damage identification.The element MSE is defined as follows: where {  } is the th normalized mode shape and [  ] is the element stiffness matrix of the th element.
Similarly, the th structural modal strain energy can be expressed as follows: where [] is the global stiffness matrix of the structure.This paper defines the fitness value of the optimal sensor placement based on the MSE as follows: where   represents the stiffness coefficient between the th DOF and th DOF;   is the deformation of the th element in the th mode;   is the deformation of the th element in the th mode; and  and  are the numbers of modes and measurement points, respectively.

The Combined Optimization Algorithm.
Although the EFI method can make the selected modal vectors approximately linearly independent, the measured MSE of the selected points may be low, which would result in the loss of modal information.Because of this effect, this paper combines the EFI method with the MSE method and defines the fitness value  3 as follows: where  is the adjustment parameter that scales the fitness value into the appropriate range.

The Basic Theory of the Quantum Genetic Algorithm
The quantum genetic algorithm (QGA) is a recently developed probability evolutionary algorithm that combines quantum computing with a genetic algorithm [33,34].Based on quantum theory, the quantum probability vector is used to encode chromosomes, and the population is updated and optimized by adopting quantum-rotating doors to search for a globally optimal solution.The smallest unit of information in the QGA is called a quantum bit and the state of a quantum bit can be expressed as follows: where ,  represent the probability amplitudes of the quantum bit and meet the following normalized condition: Thus, the state of a quantum bit can also be written as follows: In the quantum genetic algorithm, quantum information is encoded by pairs of complex numbers.The quantum chromosome composed of  sets of quantum bits can be written as follows: where This method can represent the random linear superposition of quantum states.For example, a chromosome with 3 quantum bits can be expressed as follows: An update of a quantum gate can be represented as where (  ) = [ cos(  ) − sin(  ) sin(  ) cos(  ) ] is the quantum-rotating gate, among which the variables are denoted as follows: where  is the adaptive variable;  is the evolutionary population; and iter max is a constant that depends on the complexity of the optimization problem.
Furthermore, to prevent the optimization from converging at local extrema, the algorithm introduces a mutation operation based on a certain probability.For instance, a quantum bit |0⟩ + |1⟩ can be transferred into |1⟩ + |0⟩ through the operation.In the actual operation, the mutation probability is generally between 0.1 and 0.01, which can both maintain the diversity of the population and prevent the algorithm from converging at local extremes.
The basic steps of the algorithm are shown in Figure 1.

The Evaluation Criterion of Optimal Sensor Placement
A good sensor network should allow the measured modal parameters to be sensitive to variations in the structural parameters and provide comprehensive information for structural damage diagnosis.In this paper, three criteria are selected to evaluate the ability of a sensor network to capture the structural dynamic response.

MAC.
Based on orthogonality theory of modal vectors, the orthogonality of modal vectors directly influences the results of modal identification.Therefore, the MAC is selected as the comparison criterion of the optimal sensor placement as follows: where MAC  is the cosine value of the intersection angle between the th modal vector Φ  and the th modal vector Φ  .The orthogonality of the modal vectors is better when the maximum value of the off-diagonal elements of the MAC matrix is smaller.Therefore, the maximum value of the offdiagonal elements of MAC matrix can be used as an index to evaluate the quality of the sensor placement scheme.

The Criterion of Antinoise Performance.
The identification of modal parameters is greatly influenced by ambient noise perturbations; under extreme conditions, the noise signal and structural vibration signal are not distinguishable.In this paper, the FIM determinant after the elimination of the DOFs of the candidate sensors is selected as the index for evaluating the modal identification ability of the sensor network under the effects of ambient noise perturbation.

The Criterion of the Condition Number.
The condition number is a criterion for evaluating the morbidity of a matrix.Larger numbers mean that the modal matrix is more morbid and the identified parameters will be severely distorted.A modal matrix with a small condition number can perfectly resist ambient distortions.Therefore, the condition number of the modal matrix is used as a criterion for evaluating each optimizing method as is shown as follows: A three-dimensional finite element model is developed to conduct vibration modal analysis of the dam-foundation system, as shown in Figure 3.The total number of elements and nodes of the model are 201,689 and 207,714, respectively.The model coordinates follow the right-hand screw rule.The direction of each axis is as follows: the -axis is perpendicular to the flow direction and the left-bank direction is positive; the -axis is parallel to the flow direction and the downstream direction is positive; the -axis is parallel to the height direction and the downward direction is positive.

Modal Analysis.
For large-scale civil structures, usually the displacements of low frequency modes could give sufficient information to describe the dynamic behavior of a structural system with sufficient accuracy to allow its health state modifications to be determined effectively [35].It has been illustrated by [36] that the first ten mode shapes are able to reflect most of the dynamic characteristics of huge engineering structures.In recent years, many scholars have adopted no more than first ten mode shapes to study the structural dynamic behavior of the large-scale concrete dams, bridges and truss structures, and so forth [26,[37][38][39].Therefore, in this paper, the first ten mode shapes of the dam are selected as the targeted modes.ABAQUS finite element software is used to conduct a dynamic modal analysis of the dam.Added-mass elements [40] are used to simulate the function of the upstream reservoir water.The first ten mode  shapes and corresponding natural frequencies are shown in Figure 4 and Table 2, respectively.

Optimal Sensor Placement.
There are 1,718 nodes on the upstream surface of the finite element model of the arch dam.In this paper, two DOFs ( direction and  direction) of each node are considered.After eliminating 57 nodes that are unsuitable for sensor installation, there are 1,681 candidate nodes and 3,362 corresponding DOFs on the upstream surface.The quantum genetic algorithm is adopted to arrange 30 sets of sensors on the upstream surface based on the three criteria described previously.
The parameters of the quantum genetic algorithm are defined as follows: the initial population size and quantum bit are set as 3,000 and 15, respectively; the initial search is superposed based on an equal probability; and the corresponding probability amplitude is defined as ± √ 2/2.To verify the computation efficiency and quality of the quantum genetic algorithm, a traditional genetic algorithm (GA) is also used.The parameters of GA are listed as follows: the initial population size is 3000; the crossover rate and the mutation rate are defined as 0.75 and 0.02, respectively.The maximum operating generations of the two algorithms are defined as 2000.Evolution lines of the two methods for different optimization criteria are shown in Figure 5.
It is clear that the QGA can always converge to a higher fitness value in fewer generations than the GA.The premature phenomenon of the GA can be effectively avoided by adopting QGA.The statistic of multiple runs of the QGA and GA based on three criteria is shown in Tables 3-5.As shown in the tables, the average fitness value and the best fitness value of the QGA out of different runs are all greater than those of the GA.The QGA possesses higher computation quality with the same operating generations comparing with the GA.
The total elapsed time of the two algorithms is shown in Figure 6.As is shown in figure, the total operating time of QGA based on three criteria is always shorter than the GA.Therefore, it can be concluded that the QGA possesses higher Bar graphs of the MAC matrix for the three criteria obtained by the quantum genetic algorithm are shown in Figure 7 and the relevant statistics of the off-diagonal elements in MAC matrixes are shown in Table 6.We can see that the average value and the maximum value of the offdiagonal elements obtained by the  3 criterion are much less than those with the other two criteria.The off-diagonal element MAC  represents the crossing angle of the two corresponding modal vectors.According to [25], the smaller the off-diagonal elements of the MAC matrix are, the more independent the different modal vectors remain.Therefore, it can be concluded that the sensor configuration obtained by the  3 criterion can improve the spatial intersection angle among the modal vectors and thus better represent the modal characteristics of the dam structure.The determinant of the FIM after removing the candidate sensors and condition number of the modal matrix are shown in Table 7 and the following conclusions could be drawn: (1) The FIM determinant after eliminating the DOFs of the candidate sensors can be seen as the measure of the uniformity between the measured mode and the results of the FEM methods.It is the indicator which evaluates the influence of the ambient noise on the measured mode [26].As is shown in Table 7, the obtained FIM determinant based on the  3 criterion is far greater than those of the other two criteria.Therefore, it can be concluded that the sensor network obtained by the proposed method possesses a better antinoise property.
(2) According to [41], the matrix condition number represents the sensibility of matrix calculation on computation error.The larger the condition number, the more ill-conditioned the modal matrix, whereas a matrix with condition number close to 1 would make the identified modal characteristics more accurate [42].As is shown in Table 7, the condition number obtained by the  3 criterion is 11, which is much less than those of the other two criteria (19 and 27).Therefore, it can be concluded that the proposed  3 criterion can effectively resist the ambient environmental perturbation and raise the accuracy of the identified modal parameters.
In conclusion, the combined optimization criterion based on a quantum genetic algorithm is shown to be better than  traditional methods.The final optimal sensor network is shown in Figure 8.

Conclusion
This paper proposes a novel optimization criterion for sensor placement by combining the effective independent method with the modal strain energy method to improve the modal identification ability.Moreover, to increase the computing efficiency and accuracy, a quantum genetic algorithm is adopted.By comparing the traditional methods, the following conclusions could be drawn: (1) By combining the quantum computing principle with a genetic algorithm, the QGA adopts the vector of a quantum state to represent the genetic encoding and updates the chromosomes with a quantum-rotating gate to improve the efficiency of the optimization.In comparison with a traditional genetic algorithm, it is clear that the algorithm possesses the property of global convergence and always converges to a higher fitness value in fewer generations than the GA.
(2) The off-diagonal elements of the MAC matrix obtained using the combined optimization criterion are much smaller than those obtained using traditional methods.It is verified that the sensor network obtained using the proposed method can maximize the spatial intersection angle among different modal vectors, which could increase the modal identification accuracy.
(3) The determinant of the FIM obtained with the proposed method is obviously greater than that obtained using traditional methods, while the condition number of the modal matrix is much lower than those obtained using the traditional methods.The obtained sensor network can effectively resist ambient perturbations, which also increase the modal identification accuracy.

Figure 1 :
Figure 1: Flow chart of the quantum genetic algorithm.

5. 1 .
Project Background.A double-curvature concrete arch dam is located midstream along the Lancang River in Yunnan Province.The elevation of the dam crest is 1,245 m, with a maximum dam height of 292 m and a dam crest length of 922.74 m.The top and bottom width of the crown cantilever are 13 m and 69.49m, respectively.The normal storage water level of the reservoir is 1,240 m, with a total storage of 15.12 billion m 3 and an effective storage of 9.895 billion m 3 , which makes it a pluriennial-regulation reservoir.The downstream elevation of the arch dam is shown in Figure 2.

Figure 2 :
Figure 2: Downstream elevation of the arch dam.

Figure 3 :
Figure 3: Three-dimensional finite element model of the damfoundation system.

Figure 4 :Figure 5 :
Figure 4: Typical vibration modes of the arch dam.

Figure 6 :
Figure 6: The comparison of total elapsed time between QGA and GA.

Figure 7 :
Figure 7: The bar graph of the optimal MAC values based on different criteria.

Figure 8 :
Figure 8: Final sensor network on the arch dam.

Table 2 :
The natural frequency of the arch dam.

Table 3 :
The statistic of the fitness value after multiple runs based on  1 criterion.

Table 4 :
The statistic of the fitness value after multiple runs based on  2 criterion.

Table 5 :
The statistic of the fitness value after multiple runs based on  3 criterion.

Table 6 :
The statistics of the off-diagonal elements of the MAC matrix based on three criteria.

Table 7 :
FIM determinant and the condition number of the vibration mode based on three criteria.