# A New Optimal Sensor Placement Strategy Based on Modified Modal Assurance Criterion and Improved Adaptive Genetic Algorithm for Structural Health Monitoring

Authors: Can He; Jianchun Xing; Juelong Li; Qiliang Yang; Ronghao Wang; Xun Zhang
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

Optimal sensor placement (OSP) is an important part in the structural health monitoring. Due to the ability of ensuring the linear independence of the tested modal vectors, the minimum modal assurance criterion (minMAC) is considered as an effective method and is used widely. However, some defects are present in this method, such as the low modal energy and the long computation time. A new OSP method named IAGA-MMAC is presented in this study to settle the issue. First, a modified modal assurance criterion (MMAC) is proposed to improve the modal energy of the selected locations. Then, an improved adaptive genetic algorithm (IAGA), which uses the root mean square of off-diagonal elements in the MMAC matrix as the fitness function, is proposed to enhance computation efficiency. A case study of sensor placement on a numerically simulated wharf structure is provided to verify the effectiveness of the IAGA-MMAC strategy, and two different methods are used as contrast experiments. A comparison of these strategies shows that the optimal results obtained by the IAGA-MMAC method have a high modal strain energy, a quick computational speed, and small off-diagonal elements in the MMAC matrix.


---

## Body

## 1. Introduction

With the development of constructive technology, more and more large-scale structures, such as suspension bridges, television towers, sea wharfs, and high-rise buildings, have been built all over the world. However, due to harsher environment, these civil infrastructures may be damaged by strong wind force, torrential rain, severe earthquake, explosion, and other abnormal loads [1]. The traditional routine visual inspection usually ignores the tiny failures and the flaws inside the structures. Therefore, structural health monitoring (SHM) [2–4] has become an important research topic in the engineering protection field. By monitoring the force condition of the structures real-timely, SHM system can detect the anomalies in response and find the possible deterioration at early stage to ensure the safety of the structures. A typical SHM system usually includes three subsystems: sensor subsystem, data processing subsystem, and condition evaluation subsystem [5]. Sensor subsystem is a fundamental part of SHM system. Modal parameter identification, damage detection, and structural condition evaluation are all based on the data acquired from the sensors. Generally, the more sensors are placed, the more information could be obtained. Nevertheless, the number of sensors is strictly constrained by the high cost of purchase and maintenance for the sensors. For instance, there are only 23, 65, and 72 accelerometers instrumented in the Tsing Ma Bridge, Tingkau Bridge, and Jiangyin Bridge, respectively [6]. The sensors number is far less than the available positions. Therefore, how to deploy the limited sensors becomes a challenging task.The traditional expert experience method is effective for the simple structures. Nevertheless, for a large-scale structure that has a large volume and complex geometry, it is very difficult to determine the optimal sensor locations only by experience. In the past few years, many technologies have been developed to achieve optimal sensor placement (OSP), such as modal kinetic energy method [7, 8], effective independence method [9, 10], QR decomposition method [11, 12], and information entropy method [13, 14]. Carne and Dohrmann [15] considered that distinguishing one modal vector from another is essential to realize modal parameter identification and proposed a famous OSP method named minMAC. The minMAC method includes two parts. The first section is modal assurance criterion (MAC). In the criterion, the off-diagonal element of the MAC matrix is considered to be an index for evaluating the angle of the corresponding two modal vectors. The smaller element indicates bigger angle and less correlation between two modal vectors. The second section is to achieve sensor placement using MAC. First, some locations are selected as an initial set based on experience and structural topology, and the remaining locations are considered as a candidate set. Second, the candidate locations are added one by one to the initial set and the combination that minimizes the maximum off-diagonal element of the MAC matrix is selected as the new initial set. Third, the second step is repeated till the fixed sensors number is gained.Although minMAC is considered as an effective OSP method and used widely, some defects still exist. Firstly, the core idea of MAC is to evaluate the correlation of the modal vectors for different placement scheme, but the modal energy of the selected location cannot be guaranteed. It means that some degree of freedoms (DOFs) with low signal to noise ratio may be selected as the optimal sensor positions; this defect could decrease the precision of the acquired data and the accuracy of modal parameters identification. Secondly, the minMAC method has a large computational complexity. To determine one location, the minMAC method needs to compute the MAC matrices of the variable combinations until all candidate locations have been selected. In order to improve the computational efficiency of the minMAC method, computational intelligence technology has been used, such as simulated annealing algorithm [16], ant colony optimization algorithm [17], particle swarm optimization algorithm [18], and monkey algorithm [19]. Due to easy coding method and quick evolution velocity, genetic algorithm (GA) has become a hot research direction. Liu et al. [20] presented a modified GA using two-dimensional array coding method instead of binary coding method; generalized genetic algorithm (GGA) was proposed by Yi et al. [21] to solve the OSP problem for high-rise buildings. Javadi et al. [22] adopted a hybrid intelligent method which was based on a combination of neural network and GA. Although these modifications make great progress, there are still some defects existing in GA [23]. The main drawback is that the crossover and mutation factors are invariant during the whole iteration cycle; this defect drops the searching ability of GA.The objective of this study is to present a modified minMAC method that has a big modal energy index and a quick computational speed. A new method termed IAGA-MMAC is proposed to achieve this goal. First, a modified modal assurance criterion (MMAC) is proposed to improve the modal energy of the selected locations by constructing the new modal shape matrix. And then, an improved adaptive genetic algorithm (IAGA), which uses root mean square of off-diagonal elements in the MMAC matrix as the fitness function, is adopted to determine the optimal sensor locations. The excellent calculation efficiency and global optimization ability of IAGA assist to enhance the computational speed of IAGA-MMAC.The remaining parts of the paper are organized as follows. Section2 introduces the traditional MAC and presents MMAC. The IAGA-MMAC strategy is described in Section 3. Section 4 provides a case study of sensor placement on a numerically simulated wharf structure to verify the effectiveness of the IAGA-MMAC method. Section 5 is the conclusion and the future work.


## 2. The Modified Modal Assurance Criterion

### 2.1. Introduction of the Traditional MAC

A basic requirement to distinguish the measured modes is that the measured modal vectors must be as linearly independent as possible. MAC provides a useful criterion to evaluate the correlation of modal vectors. MAC matrix is defined as(1)MACij=ΦiTΦj2ΦiTΦiΦjTΦj,
where Φi and Φj are the ith and jth column vectors in the modal shape matrix Φ.The off-diagonal elements in the MAC matrix express the correlation between two modal vectors. Identifying theith mode and the jth mode is easy if the value of MACij is small.


### 2.2. Description of MMAC

The main idea of MAC is to guarantee the linear independence of the modal vectors. However, some DOFs with low energy may be selected as sensor locations using MAC; this defect could decrease the signal to noise ratio of the acquired data and make troubles in modal shape identification and structural damage detection.As shown in1, the element of the MAC matrix is computed using the column vectors of the modal shape matrix Φ. In general, modal shape matrix should be composed by all orders of modes. But for the sake of reducing computation complexity, only a few modal orders are selected to compose the modal shape matrix. To improve the modal energy of the selected locations, the modal orders having large dynamical response should be selected. Therefore, how to evaluate the dynamical feature of different modal orders becomes an important work. In tradition, the first several modes are considered to have large dynamical response. Nevertheless, the modal energy of a structure is not always concentrated on the low-order modes. To some complex structures, such as high-rise building and suspension bridge, some high modal orders also have large dynamical response [24]. Therefore, it is not suitable that only the first several modes are selected to compose the modal shape matrix. In order to improve the modal energy of sensor locations, a modified MAC is presented to settle the issue in this study. In the new criterion, modal participation factor (MPF) is presented as the criterion to evaluate the dynamical response of different modal orders. And then, the new modal shape matrix Φ- is composed by the modal orders with big MPF. Finally, the modified MAC matrix is computed using Φ-.


#### 2.2.1. Modal Order Selection Using MPF

Different choices of modal order result in the various placement results. According to the requirement of OSP, the modal order having a strong dynamic response should be selected. Nevertheless, modal order is selected based on experience in the traditional methods. In order to solve this problem, MPF is presented as the evaluation criterion for the modal order selection.The equation of motion for then DOF dynamic system is represented as
(2)Mu¨+Cu˙+Ku=-MDu¨g,
where M is the structural mass matrix, C is the structural damping matrix, K is the structural stiffness matrix, D is the vector describing the excitation direction, u is the modal displacement vector, and u¨g is the acceleration generated by input force.The displacement vectoru can be transformed by
(3)u=Φ·q,
where Φ is the modal shape matrix and q is the modal coordinate.IfM, C, and K meet the orthogonalization of the normalization modal shape, 2 can be written as
(4)q¨i+2ζiwiq˙i+wi2qi=-ΦiTMDΦiTMΦiu¨gt,
where qi is the ith modal coordinate, wi is the natural frequency of the ith mode, and ζi is the fraction of critical damping of the ith mode.The MPF is defined as follows in [25]
(5)φir=ΦiTMDrΦiTMΦi,
where φir is the MPF of the ith mode in r direction.The mass normalization for the modal shape matrixΦ is determined by
(6)ΦiTMΦi=1.Based on5 and 6, the MPF can be expressed as
(7)φir=ΦiTMDr,
where Dr is a vector describing the excitation in r direction, and it can be computed as
(8)DrT=T·Er,
where Er is a unit vector in r direction and T is a constant matrix that indicates the relationship between the excitation and reference points:
(9)T=1000Z-Z0Y0-Y010Z0-Z0X-X0001Y-Y0X0-X0000100000010000001,
where X0, Y0, and Z0 are the global Cartesian coordinates of the reference point and X, Y, and Z are the global Cartesian coordinates of the excitation point.MPF reflects the amplitude of dynamical response for all the modes. Therefore, the modal order that has big MPF should be selected priorly to construct the new modal shape.Although all of the modes can be sorted based on their MPFs, there is still a question that how many modes should be selected. The modal participating mass ratio is presented as a criterion to determine the number of the modes. Modal participating mass is defined as(10)Mpi=ΦiTMD2ΦiTMΦi,
where Mpi is the modal participating mass of the ith mode.When mass normalization has been conducted for the modal shape matrixΦ, 10 can be transformed by
(11)Mpi=ΦiTMD2.The modal participating mass ratio of theith mode can be expressed as
(12)ri=MpiDTMD.Wilson [25] suggested that the summation of modal participating mass ratio of the selected modes should be over 90%. This criterion is employed in this study to determine the number of modes.


#### 2.2.2. Computation of the MMAC Matrix

The dynamical response of different modal orders can be sorted by comparing their MPFs; the new modal shape matrixΦ- is constructed by the modal orders with big MPFs. Therefore, the modified MAC (MMAC) matrix can be computed as
(13)MMACij=Φ¯iTΦ¯j2Φ¯iTΦ¯iΦ¯jTΦ¯j.The computation process of the modified MAC matrix is shown as Figure1.Figure 1
Computation process of the modified MAC matrix.


## 2.1. Introduction of the Traditional MAC

A basic requirement to distinguish the measured modes is that the measured modal vectors must be as linearly independent as possible. MAC provides a useful criterion to evaluate the correlation of modal vectors. MAC matrix is defined as(1)MACij=ΦiTΦj2ΦiTΦiΦjTΦj,
where Φi and Φj are the ith and jth column vectors in the modal shape matrix Φ.The off-diagonal elements in the MAC matrix express the correlation between two modal vectors. Identifying theith mode and the jth mode is easy if the value of MACij is small.


## 2.2. Description of MMAC

The main idea of MAC is to guarantee the linear independence of the modal vectors. However, some DOFs with low energy may be selected as sensor locations using MAC; this defect could decrease the signal to noise ratio of the acquired data and make troubles in modal shape identification and structural damage detection.As shown in1, the element of the MAC matrix is computed using the column vectors of the modal shape matrix Φ. In general, modal shape matrix should be composed by all orders of modes. But for the sake of reducing computation complexity, only a few modal orders are selected to compose the modal shape matrix. To improve the modal energy of the selected locations, the modal orders having large dynamical response should be selected. Therefore, how to evaluate the dynamical feature of different modal orders becomes an important work. In tradition, the first several modes are considered to have large dynamical response. Nevertheless, the modal energy of a structure is not always concentrated on the low-order modes. To some complex structures, such as high-rise building and suspension bridge, some high modal orders also have large dynamical response [24]. Therefore, it is not suitable that only the first several modes are selected to compose the modal shape matrix. In order to improve the modal energy of sensor locations, a modified MAC is presented to settle the issue in this study. In the new criterion, modal participation factor (MPF) is presented as the criterion to evaluate the dynamical response of different modal orders. And then, the new modal shape matrix Φ- is composed by the modal orders with big MPF. Finally, the modified MAC matrix is computed using Φ-.


### 2.2.1. Modal Order Selection Using MPF

Different choices of modal order result in the various placement results. According to the requirement of OSP, the modal order having a strong dynamic response should be selected. Nevertheless, modal order is selected based on experience in the traditional methods. In order to solve this problem, MPF is presented as the evaluation criterion for the modal order selection.The equation of motion for then DOF dynamic system is represented as
(2)Mu¨+Cu˙+Ku=-MDu¨g,
where M is the structural mass matrix, C is the structural damping matrix, K is the structural stiffness matrix, D is the vector describing the excitation direction, u is the modal displacement vector, and u¨g is the acceleration generated by input force.The displacement vectoru can be transformed by
(3)u=Φ·q,
where Φ is the modal shape matrix and q is the modal coordinate.IfM, C, and K meet the orthogonalization of the normalization modal shape, 2 can be written as
(4)q¨i+2ζiwiq˙i+wi2qi=-ΦiTMDΦiTMΦiu¨gt,
where qi is the ith modal coordinate, wi is the natural frequency of the ith mode, and ζi is the fraction of critical damping of the ith mode.The MPF is defined as follows in [25]
(5)φir=ΦiTMDrΦiTMΦi,
where φir is the MPF of the ith mode in r direction.The mass normalization for the modal shape matrixΦ is determined by
(6)ΦiTMΦi=1.Based on5 and 6, the MPF can be expressed as
(7)φir=ΦiTMDr,
where Dr is a vector describing the excitation in r direction, and it can be computed as
(8)DrT=T·Er,
where Er is a unit vector in r direction and T is a constant matrix that indicates the relationship between the excitation and reference points:
(9)T=1000Z-Z0Y0-Y010Z0-Z0X-X0001Y-Y0X0-X0000100000010000001,
where X0, Y0, and Z0 are the global Cartesian coordinates of the reference point and X, Y, and Z are the global Cartesian coordinates of the excitation point.MPF reflects the amplitude of dynamical response for all the modes. Therefore, the modal order that has big MPF should be selected priorly to construct the new modal shape.Although all of the modes can be sorted based on their MPFs, there is still a question that how many modes should be selected. The modal participating mass ratio is presented as a criterion to determine the number of the modes. Modal participating mass is defined as(10)Mpi=ΦiTMD2ΦiTMΦi,
where Mpi is the modal participating mass of the ith mode.When mass normalization has been conducted for the modal shape matrixΦ, 10 can be transformed by
(11)Mpi=ΦiTMD2.The modal participating mass ratio of theith mode can be expressed as
(12)ri=MpiDTMD.Wilson [25] suggested that the summation of modal participating mass ratio of the selected modes should be over 90%. This criterion is employed in this study to determine the number of modes.


### 2.2.2. Computation of the MMAC Matrix

The dynamical response of different modal orders can be sorted by comparing their MPFs; the new modal shape matrixΦ- is constructed by the modal orders with big MPFs. Therefore, the modified MAC (MMAC) matrix can be computed as
(13)MMACij=Φ¯iTΦ¯j2Φ¯iTΦ¯iΦ¯jTΦ¯j.The computation process of the modified MAC matrix is shown as Figure1.Figure 1
Computation process of the modified MAC matrix.


## 2.2.1. Modal Order Selection Using MPF

Different choices of modal order result in the various placement results. According to the requirement of OSP, the modal order having a strong dynamic response should be selected. Nevertheless, modal order is selected based on experience in the traditional methods. In order to solve this problem, MPF is presented as the evaluation criterion for the modal order selection.The equation of motion for then DOF dynamic system is represented as
(2)Mu¨+Cu˙+Ku=-MDu¨g,
where M is the structural mass matrix, C is the structural damping matrix, K is the structural stiffness matrix, D is the vector describing the excitation direction, u is the modal displacement vector, and u¨g is the acceleration generated by input force.The displacement vectoru can be transformed by
(3)u=Φ·q,
where Φ is the modal shape matrix and q is the modal coordinate.IfM, C, and K meet the orthogonalization of the normalization modal shape, 2 can be written as
(4)q¨i+2ζiwiq˙i+wi2qi=-ΦiTMDΦiTMΦiu¨gt,
where qi is the ith modal coordinate, wi is the natural frequency of the ith mode, and ζi is the fraction of critical damping of the ith mode.The MPF is defined as follows in [25]
(5)φir=ΦiTMDrΦiTMΦi,
where φir is the MPF of the ith mode in r direction.The mass normalization for the modal shape matrixΦ is determined by
(6)ΦiTMΦi=1.Based on5 and 6, the MPF can be expressed as
(7)φir=ΦiTMDr,
where Dr is a vector describing the excitation in r direction, and it can be computed as
(8)DrT=T·Er,
where Er is a unit vector in r direction and T is a constant matrix that indicates the relationship between the excitation and reference points:
(9)T=1000Z-Z0Y0-Y010Z0-Z0X-X0001Y-Y0X0-X0000100000010000001,
where X0, Y0, and Z0 are the global Cartesian coordinates of the reference point and X, Y, and Z are the global Cartesian coordinates of the excitation point.MPF reflects the amplitude of dynamical response for all the modes. Therefore, the modal order that has big MPF should be selected priorly to construct the new modal shape.Although all of the modes can be sorted based on their MPFs, there is still a question that how many modes should be selected. The modal participating mass ratio is presented as a criterion to determine the number of the modes. Modal participating mass is defined as(10)Mpi=ΦiTMD2ΦiTMΦi,
where Mpi is the modal participating mass of the ith mode.When mass normalization has been conducted for the modal shape matrixΦ, 10 can be transformed by
(11)Mpi=ΦiTMD2.The modal participating mass ratio of theith mode can be expressed as
(12)ri=MpiDTMD.Wilson [25] suggested that the summation of modal participating mass ratio of the selected modes should be over 90%. This criterion is employed in this study to determine the number of modes.


## 2.2.2. Computation of the MMAC Matrix

The dynamical response of different modal orders can be sorted by comparing their MPFs; the new modal shape matrixΦ- is constructed by the modal orders with big MPFs. Therefore, the modified MAC (MMAC) matrix can be computed as
(13)MMACij=Φ¯iTΦ¯j2Φ¯iTΦ¯iΦ¯jTΦ¯j.The computation process of the modified MAC matrix is shown as Figure1.Figure 1
Computation process of the modified MAC matrix.


## 3. Optimal Sensor Placement Based on the IAGA-MMAC Method

The traditional GA has a shortcoming that its crossover and mutation factors are stationary during the whole iteration cycle; this defect could drop the evolution velocity and searching ability of the algorithm. Therefore, a new OSP method named IAGA-MMAC is presented to settle the issue. The flowchart of the IAGA-MMAC is shown in Figure2.Figure 2
Flowchart of the IAGA-MMAC method.


### 3.1. Coding Method

In the traditional coding method, the gene number of value 1 could be changed in the processes of crossover and mutation, which may result in two or more sensors placed on one location. Therefore, the dual-structure coding method is used in this study. The chromosome of individual is composed of two rows, where the upper rows(i) named append code indicates all of the locations, and the lower row Xs(i) named variable code represents whether the location s(i) places a sensor with the value 0 or 1. The dual-structure coding method is shown as Table 1.Table 1
Dual-structure coding method.
Location code
s
(
1
)
s
(
2
)
⋯
s
(
n
)
Variable code
X
s
(
1
)
X
s
(
2
)
⋯
X
s
(
n
)


### 3.2. Selection Scheme

The nonlinear ranking selection scheme is used to choose the good individual in one generation. Firstly, the individuals are sorted by comparing their fitness values. And then, the nonlinear function14 is adopted to determine the selection probability of different individuals. Compared with the common roulette wheel selection scheme, this method does not depend on the concrete values of individual fitness; therefore, it can avoid the scale conversion of fitness values. This advantage can inhibit the algorithm premature convergence:
(14)Pxn=q1-qi-1,
where q is the selection probability of best individual and i is the sort number of the fitness value of xn.


### 3.3. Crossover and Mutation

The crossover and mutation factors can make great influence on the computation efficiency of GA. In the traditional GA, the factors keep invariant during the whole iteration cycle. This defect drops the evolution velocity and optimization ability of the algorithm. Therefore, the adaptive mechanism is introduced into this algorithm [26]. The adaptive adjustment processes of the two factors are defined as
(15)Pc=k1(fc-favg)fmax⁡-favg,fc≥favg,k2,fc<favg,Pm=k3(fm-favg)fmax⁡-favg,fm≥favg,k4,fm<favg,
where fmax⁡ is the maximum of fitness value, favg is the average of population fitness value, fc represents the bigger fitness value between the two crossover individuals, fm is the biggest fitness value among the mutation individuals, k1 and k2 are the adjustment coefficients of crossover factor, and k3 and k4 are the adjustment coefficients of mutation factor.As shown in15, when fc and fm are close or equal to fmax⁡, the Pc and Pm will decline to zero. This problem leads to the good individuals update very slowly at the beginning of the evolution process. In order to address this issue, the improved adaptive crossover and mutation factors are presented as
(16)Pc=p1-p1-p2fc-favgfmax⁡-favg,fc≥favg,p1,fc<favg,Pm=p3-(p3-p4)(fm-favg)fmax⁡-favg,fm≥favg,p3,fm<favg,
where p1 and p2 represent the maximum and minimum values of the crossover probability, respectively, and p3 and p4 are the maximum and minimum values of the mutation probability, respectively.Whenfc and fm are close to fmax⁡, the new adaptive mechanism ensures that the crossover and mutation factors of good individual increased to p2 and p4. This advantage can help the good individual update quickly in the whole evolution cycle.


### 3.4. Fitness Function

The off-diagonal element of the MMAC matrix expresses the correlation between the corresponding two modal vectors. The elements are valued between 0 and 1, and it is easy to identify theith mode and jth mode if the value of MMACij is small. When MMACij=0, the two modal vectors are orthogonal. Therefore, the off-diagonal elements of MMAC matrix can be considered as an optimization objective for the OSP problem: the smaller element value represents the better sensor placement result.The fitness function in IAGA is asked to be a maximization problem. Therefore, the fitness function is constructed as(17)f=1-1N∑i=1NXi2,
where Xi is the off-diagonal elements of the MMAC matrix.


## 3.1. Coding Method

In the traditional coding method, the gene number of value 1 could be changed in the processes of crossover and mutation, which may result in two or more sensors placed on one location. Therefore, the dual-structure coding method is used in this study. The chromosome of individual is composed of two rows, where the upper rows(i) named append code indicates all of the locations, and the lower row Xs(i) named variable code represents whether the location s(i) places a sensor with the value 0 or 1. The dual-structure coding method is shown as Table 1.Table 1
Dual-structure coding method.
Location code
s
(
1
)
s
(
2
)
⋯
s
(
n
)
Variable code
X
s
(
1
)
X
s
(
2
)
⋯
X
s
(
n
)


## 3.2. Selection Scheme

The nonlinear ranking selection scheme is used to choose the good individual in one generation. Firstly, the individuals are sorted by comparing their fitness values. And then, the nonlinear function14 is adopted to determine the selection probability of different individuals. Compared with the common roulette wheel selection scheme, this method does not depend on the concrete values of individual fitness; therefore, it can avoid the scale conversion of fitness values. This advantage can inhibit the algorithm premature convergence:
(14)Pxn=q1-qi-1,
where q is the selection probability of best individual and i is the sort number of the fitness value of xn.


## 3.3. Crossover and Mutation

The crossover and mutation factors can make great influence on the computation efficiency of GA. In the traditional GA, the factors keep invariant during the whole iteration cycle. This defect drops the evolution velocity and optimization ability of the algorithm. Therefore, the adaptive mechanism is introduced into this algorithm [26]. The adaptive adjustment processes of the two factors are defined as
(15)Pc=k1(fc-favg)fmax⁡-favg,fc≥favg,k2,fc<favg,Pm=k3(fm-favg)fmax⁡-favg,fm≥favg,k4,fm<favg,
where fmax⁡ is the maximum of fitness value, favg is the average of population fitness value, fc represents the bigger fitness value between the two crossover individuals, fm is the biggest fitness value among the mutation individuals, k1 and k2 are the adjustment coefficients of crossover factor, and k3 and k4 are the adjustment coefficients of mutation factor.As shown in15, when fc and fm are close or equal to fmax⁡, the Pc and Pm will decline to zero. This problem leads to the good individuals update very slowly at the beginning of the evolution process. In order to address this issue, the improved adaptive crossover and mutation factors are presented as
(16)Pc=p1-p1-p2fc-favgfmax⁡-favg,fc≥favg,p1,fc<favg,Pm=p3-(p3-p4)(fm-favg)fmax⁡-favg,fm≥favg,p3,fm<favg,
where p1 and p2 represent the maximum and minimum values of the crossover probability, respectively, and p3 and p4 are the maximum and minimum values of the mutation probability, respectively.Whenfc and fm are close to fmax⁡, the new adaptive mechanism ensures that the crossover and mutation factors of good individual increased to p2 and p4. This advantage can help the good individual update quickly in the whole evolution cycle.


## 3.4. Fitness Function

The off-diagonal element of the MMAC matrix expresses the correlation between the corresponding two modal vectors. The elements are valued between 0 and 1, and it is easy to identify theith mode and jth mode if the value of MMACij is small. When MMACij=0, the two modal vectors are orthogonal. Therefore, the off-diagonal elements of MMAC matrix can be considered as an optimization objective for the OSP problem: the smaller element value represents the better sensor placement result.The fitness function in IAGA is asked to be a maximization problem. Therefore, the fitness function is constructed as(17)f=1-1N∑i=1NXi2,
where Xi is the off-diagonal elements of the MMAC matrix.


## 4. Demonstration Case

In order to verify the effectiveness of the IAGA-MMAC method, a case study of sensor placement on a numerically simulated wharf structure is provided below.


### 4.1. Modeling and Modal Analysis

A high-piled wharf is selected as a case study. The total length of wharf is 55 m, and the width of anterior bearing platform is 30 m. The upper structure of the platform is composition of cross beams, longitudinal beams, and reinforced concrete decks. The platform is supported by six rows of vertical piles, each row has five piles, and the distance between adjacent piles is 5.50 m. The finite element (FE) model of the wharf is established using commercial software Ansys, and the FE model is shown in Figure3.Figure 3
FE model of the wharf.The cross beams and longitudinal beams are modeled by beam4, vertical piles are simulated by pipe16, and reinforced concrete decks are built by shell63. 15147 nodes are generated in total. The natural frequency and modal shape of the structure are received as Table2 and Figure 4 by modal analysis.Table 2
Natural frequency of the first 10 modes.
Modal order
1
2
3
4
5
Frequency (Hz)
2.592
2.620
2.818
6.902
7.339
Mode order
6
7
8
9
10
Frequency (Hz)
7.501
8.011
8.261
9.035
9.948Figure 4
The 1st–4th modes shape.
(a)
1st mode
(b)
2nd mode
(c)
3rd mode
(d)
4th mode


### 4.2. Computation of the MMAC Matrix

The MPFs of the modes are computed using7 to 9. The results of the first 50 modes are shown in Figures 5, 6, and 7.Figure 5
MPF ofx direction.Figure 6
MPF ofy direction.Figure 7
MPF ofz direction.As shown in Figures5 to 7, the 1st, 3rd, 5th, and 11th modal orders have big MPFs in x direction. The 2nd, 4th, and 6th modal orders have large MPFs in y direction. MPFs of the 4th, 6th, and 8th modal orders are big in z direction. Therefore the 1st, 2nd, 3rd, 4th, 5th, 6th, 8th, and 11th modes are selected as the optimal combination. The modal participating mass ratio of the combination in each direction is computed using 10 to 12, and the results are presented as
(18)rx=95.2%;ry=95.1%;rz=98.9%,
where rx is the modal participating mass ratio in x direction.The results show thatrx, ry, and rz are all over 90%. Therefore, the selected modal order combination is reasonable.The new modal shape matrixΦ- is composed as
(19)Φ-=Φ1,Φ2,Φ3,Φ4,Φ5,Φ6,Φ8,Φ11,
where Φi is the ith column vector of the original modal shape matrix Φ.The MMAC matrix is computed usingΦ- according to 13.


### 4.3. Optimal Sensor Placement Using IAGA-MMAC

IAGA, which uses the off-diagonal elements of the MMAC matrix as fitness function, is adopted to achieve the sensor placement. According to the advice of [26], the parameters of IAGA are set as in Table 3.Table 3
Parameter setting.
Population
Generation
P
1
P
2
P
3
P
4
100
500
0.9
0.6
0.1
0.01The optimal process of IAGA-MMAC is executed according to Figure2. The algorithm has been run for 50 times with a certain sensor number. The best fitness function values with different sensors numbers are shown in Figure 8.Figure 8
Fitness function values with different sensor numbers.As shown in Figure8, when the sensor number is set between 5 and 11, the fitness increases with sensor added. The fitness function reaches its maximum when the sensor number is set at 11. But if the sensor number is more than 11, the fitness does not increase with sensor added. The reason for the contradiction is that the row vector determined at the newly included sensor location has a strong linear relationship with the entire previous sensor set. Therefore, eleven is determined as the optimal sensor number in consideration of placement cost and efficiency. The final optimal result is in Table 4.Table 4
Optimal sensor locations.
Sensor number
1
2
3
4
5
6
Node
387
597
1445
2544
2769
3384
Direction
x
z
y
z
x
z
Sensor number
7
8
9
10
11
Node
4272
5423
7392
11060
14028
Direction
z
x
z
y
z


### 4.4. Contrast Experiment

Three different methods are used to deploy the sensors for the wharf, including AGA with the traditional MAC as fitness function (AGA-MAC) [27], IAGA with the traditional MAC (IAGA-MAC), and IAGA-MMAC. The parameters of AGA are the same as IAGA. The optimal placement results of three methods are provided in Figures 9, 10, and 11.Figure 9
Sensor placement result of AGA-MAC.Figure 10
Sensor placement result of IAGA-MAC.Figure 11
Sensor placement result of IAGA-MMAC.The modal energy index and MAC index are selected to evaluate the placement effect of three methods. The summation of modal strain energy (MSE) for all the selected locations is shown in Figure12, and the energy of single sensor location is provided by Figure 13.Figure 12
MSE summation of all the sensor locations.Figure 13
MSE of single sensor location.The MAC matrices of three methods are shown in Figures14, 15, and 16. The RMS and maximum of off-diagonal elements in the MAC matrix are provided in Table 5.Table 5
MAC index of three methods.
Method
RMS of off-diagonal elements
Maximum of off-diagonal elements
AGA-MAC
0.0349
0.1228
IAGA-MAC
0.0202
0.0631
IAGA-MMAC
0.0162
0.0371Figure 14
MAC matrix of AGA-MAC.Figure 15
MAC matrix of IAGA-MAC.Figure 16
MMAC matrix of IAGA-MMAC.The evolution processes of three methods are shown in Figure17.Figure 17
Evolution curve of three methods.


### 4.5. Result Analysis

As shown in Figures12 and 13, the MSE summation of all the selected locations by IAGA-MMAC is bigger than that by AGA-MAC and IAGA-MAC obviously, and the MSE of single sensor location by IAGA-MMAC is also much larger. These results imply that IAGA-MMAC is a better method in terms of modal energy index. The reason for the results is that the MMAC matrix is computed using the new modal shape matrix Φ-, which is composed by the modal vectors with large dynamical response.The data in Table5 show that IAGA-MMAC gets the smaller RMS and maximum of off-diagonal elements than the other two methods. A similar conclusion can be achieved from Figures 14 to 16. The result in Figure 17 shows that the IAGA-MMAC method has a faster evolution velocity and better fitness function value than AGA-MAC. The reason for the result is that the improved adaptive crossover and mutation factors enhance the searching efficiency and global optimization ability of the algorithm.


### 4.6. Discussion

The demonstration case is based on a numerically simulated wharf structure. The FE model of the wharf structure is built in strict accordance with the actual structure. Therefore, the numerical simulation can verify the effectiveness of the IAGA-MMAC method to a certain extent. Of course, the numerical FE model maintains some limitations inevitably. Firstly, noise and other environmental factors, such as temperature and wind force, may impact the dynamic feature of the real wharf structure. However, these factors are not considered in the simulation experiment. Secondly, in the real structure, there are some locations that are not available to deploy the sensors. But these nodes may be selected as the optimal sensor locations in the numerical analysis.


## 4.1. Modeling and Modal Analysis

A high-piled wharf is selected as a case study. The total length of wharf is 55 m, and the width of anterior bearing platform is 30 m. The upper structure of the platform is composition of cross beams, longitudinal beams, and reinforced concrete decks. The platform is supported by six rows of vertical piles, each row has five piles, and the distance between adjacent piles is 5.50 m. The finite element (FE) model of the wharf is established using commercial software Ansys, and the FE model is shown in Figure3.Figure 3
FE model of the wharf.The cross beams and longitudinal beams are modeled by beam4, vertical piles are simulated by pipe16, and reinforced concrete decks are built by shell63. 15147 nodes are generated in total. The natural frequency and modal shape of the structure are received as Table2 and Figure 4 by modal analysis.Table 2
Natural frequency of the first 10 modes.
Modal order
1
2
3
4
5
Frequency (Hz)
2.592
2.620
2.818
6.902
7.339
Mode order
6
7
8
9
10
Frequency (Hz)
7.501
8.011
8.261
9.035
9.948Figure 4
The 1st–4th modes shape.
(a)
1st mode
(b)
2nd mode
(c)
3rd mode
(d)
4th mode


## 4.2. Computation of the MMAC Matrix

The MPFs of the modes are computed using7 to 9. The results of the first 50 modes are shown in Figures 5, 6, and 7.Figure 5
MPF ofx direction.Figure 6
MPF ofy direction.Figure 7
MPF ofz direction.As shown in Figures5 to 7, the 1st, 3rd, 5th, and 11th modal orders have big MPFs in x direction. The 2nd, 4th, and 6th modal orders have large MPFs in y direction. MPFs of the 4th, 6th, and 8th modal orders are big in z direction. Therefore the 1st, 2nd, 3rd, 4th, 5th, 6th, 8th, and 11th modes are selected as the optimal combination. The modal participating mass ratio of the combination in each direction is computed using 10 to 12, and the results are presented as
(18)rx=95.2%;ry=95.1%;rz=98.9%,
where rx is the modal participating mass ratio in x direction.The results show thatrx, ry, and rz are all over 90%. Therefore, the selected modal order combination is reasonable.The new modal shape matrixΦ- is composed as
(19)Φ-=Φ1,Φ2,Φ3,Φ4,Φ5,Φ6,Φ8,Φ11,
where Φi is the ith column vector of the original modal shape matrix Φ.The MMAC matrix is computed usingΦ- according to 13.


## 4.3. Optimal Sensor Placement Using IAGA-MMAC

IAGA, which uses the off-diagonal elements of the MMAC matrix as fitness function, is adopted to achieve the sensor placement. According to the advice of [26], the parameters of IAGA are set as in Table 3.Table 3
Parameter setting.
Population
Generation
P
1
P
2
P
3
P
4
100
500
0.9
0.6
0.1
0.01The optimal process of IAGA-MMAC is executed according to Figure2. The algorithm has been run for 50 times with a certain sensor number. The best fitness function values with different sensors numbers are shown in Figure 8.Figure 8
Fitness function values with different sensor numbers.As shown in Figure8, when the sensor number is set between 5 and 11, the fitness increases with sensor added. The fitness function reaches its maximum when the sensor number is set at 11. But if the sensor number is more than 11, the fitness does not increase with sensor added. The reason for the contradiction is that the row vector determined at the newly included sensor location has a strong linear relationship with the entire previous sensor set. Therefore, eleven is determined as the optimal sensor number in consideration of placement cost and efficiency. The final optimal result is in Table 4.Table 4
Optimal sensor locations.
Sensor number
1
2
3
4
5
6
Node
387
597
1445
2544
2769
3384
Direction
x
z
y
z
x
z
Sensor number
7
8
9
10
11
Node
4272
5423
7392
11060
14028
Direction
z
x
z
y
z


## 4.4. Contrast Experiment

Three different methods are used to deploy the sensors for the wharf, including AGA with the traditional MAC as fitness function (AGA-MAC) [27], IAGA with the traditional MAC (IAGA-MAC), and IAGA-MMAC. The parameters of AGA are the same as IAGA. The optimal placement results of three methods are provided in Figures 9, 10, and 11.Figure 9
Sensor placement result of AGA-MAC.Figure 10
Sensor placement result of IAGA-MAC.Figure 11
Sensor placement result of IAGA-MMAC.The modal energy index and MAC index are selected to evaluate the placement effect of three methods. The summation of modal strain energy (MSE) for all the selected locations is shown in Figure12, and the energy of single sensor location is provided by Figure 13.Figure 12
MSE summation of all the sensor locations.Figure 13
MSE of single sensor location.The MAC matrices of three methods are shown in Figures14, 15, and 16. The RMS and maximum of off-diagonal elements in the MAC matrix are provided in Table 5.Table 5
MAC index of three methods.
Method
RMS of off-diagonal elements
Maximum of off-diagonal elements
AGA-MAC
0.0349
0.1228
IAGA-MAC
0.0202
0.0631
IAGA-MMAC
0.0162
0.0371Figure 14
MAC matrix of AGA-MAC.Figure 15
MAC matrix of IAGA-MAC.Figure 16
MMAC matrix of IAGA-MMAC.The evolution processes of three methods are shown in Figure17.Figure 17
Evolution curve of three methods.


## 4.5. Result Analysis

As shown in Figures12 and 13, the MSE summation of all the selected locations by IAGA-MMAC is bigger than that by AGA-MAC and IAGA-MAC obviously, and the MSE of single sensor location by IAGA-MMAC is also much larger. These results imply that IAGA-MMAC is a better method in terms of modal energy index. The reason for the results is that the MMAC matrix is computed using the new modal shape matrix Φ-, which is composed by the modal vectors with large dynamical response.The data in Table5 show that IAGA-MMAC gets the smaller RMS and maximum of off-diagonal elements than the other two methods. A similar conclusion can be achieved from Figures 14 to 16. The result in Figure 17 shows that the IAGA-MMAC method has a faster evolution velocity and better fitness function value than AGA-MAC. The reason for the result is that the improved adaptive crossover and mutation factors enhance the searching efficiency and global optimization ability of the algorithm.


## 4.6. Discussion

The demonstration case is based on a numerically simulated wharf structure. The FE model of the wharf structure is built in strict accordance with the actual structure. Therefore, the numerical simulation can verify the effectiveness of the IAGA-MMAC method to a certain extent. Of course, the numerical FE model maintains some limitations inevitably. Firstly, noise and other environmental factors, such as temperature and wind force, may impact the dynamic feature of the real wharf structure. However, these factors are not considered in the simulation experiment. Secondly, in the real structure, there are some locations that are not available to deploy the sensors. But these nodes may be selected as the optimal sensor locations in the numerical analysis.


## 5. Conclusions and Future Work

There are some shortcomings existing in the traditional minMAC method, such as a low modal energy and a large computational complexity. In order to solve these defects, a new method termed IAGA-MMAC is presented in this study. In the strategy, the traditional MAC is modified by constructing the new modal shape matrix to improve modal strain energy of the selected sensor locations. And then, IAGA, using the MMAC matrix as optimization objective function, is adopted to enhance the computation efficiency. With the case analysis, some conclusions are summarized as follows.(1)
The traditional MAC only concerns the correlation of the measured modal vectors. Therefore, it is necessary for the criterion to be modified to improve the modal strain energy of the selected locations. In the modified criterion, MPF is adopted to evaluate the dynamical feature of different modal orders firstly. And then, the new modal shape matrixΦ- is constructed by the modal orders which have big vibration response. Finally, the MMAC matrix is computed using the new modal shape matrix. As experiment results showed, MMAC can enhance the MSE of the selected locations.(2)
Due to the improved adaptive crossover and mutation factors, IAGA has a good searching ability and computation efficiency. Therefore, IAGA using the MMAC matrix as fitness function is presented in this study to achieve OSP. Experiment results show that IAGA-MMAC has a faster computational speed and a better optimization ability than AGA-MAC and IAGA-MAC.The IAGA-MMAC method is proved to be an effective OSP method as the numerical analysis. However, there are still some works to do. An experiment of sensor placement on a real wharf structure will be conducted by our team in the further research.



---

*Source: 626342-2015-03-11.xml*