International Journal of Distributed Sensor Networks

Volume 2012 (2012), Article ID 454530, 7 pages

doi:10.1155/2012/454530

Research Article

Research on Optimal Sensor Placement Based on Reverberation Matrix for Structural Health Monitoring

Hai-feng Yang,¹ Zi-yan Wu,¹ Shu-kui Liu,1,2 and Hong-bin Sun¹

¹School of Mechanics, Civil Engineering and Architecture, Northwestern Polytechnical University, Xi’an 710072, China

²Cockrell School of Engineering, The University of Texas at Austin, Austin, TX 78712, USA

Received 28 June 2012; Revised 28 September 2012; Accepted 29 September 2012

Academic Editor: Hong-Nan Li

Copyright © 2012 Hai-feng Yang 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.

Adequate sensor placement plays a key role in such fields as system identification, structural control, damage detection, and structural health monitoring (SHM) of large-scale civil infrastructures. Many optimal sensor placement (OSP) methods have been developed for general optimized solution searches. Due to the limitations of equipment facilities and cost, the number of sensors to be installed in a structure is relatively few. It is very important to determine the necessary number of sensors to be installed and where to deploy these sensors. Taking into account energy attenuate during the signal propagation, combined with classic reverberation matrix method, a two-step method is proposed to determine the sensors arrangement using the scattering matrix in this paper. First, calculate the utmost distance of wave propagation on a special structure by the principle of elastic wave propagation and determine the preliminary number of sensors; second, the utmost distance and number of sensors are applied to a sensor optimization algorithm named Effective Independence Driving-Point Residue method. In the bridge benchmark model case study, it shows the validity of proposed method under a special detection system.

1. Introduction

A typical SHM system includes three major components: a sensor system, a data processing system (including the data acquisition, transmission, and storage), and a health evaluation system [1–3]. The most common used methods of data acquisition and analysis are based on vibration test (active excitation or ambient excitation). The number and location of sensors directly affect the accuracy of the identified modal parameters. Thus, the sensors must be judiciously placed in order to provide adequate information for the identification of the structural behavior. On the other hand, using too many sensors as well as accessorial facilities will also cause exorbitant cost.

In the past years, many methods have been developed for general optimized solution searches such as effective independence method [4], kinetic energy method [5], and some modern-type methods [6–10]. These methods solve the sensor optimization position through, as much as possible, information getting from the data acquisition system based on the analysis of finite element model (FEM). Another approaches range from applying the constraints on the objective function (evaluation criteria) to streamline the optimization process to applying the advanced artificial analysis techniques such as the genetic algorithms and simulated annealing.

As illustrated in [11], damage detection techniques based on modal parameters represent the “traditional” approach to the problem. This approach could appear convenient when a complete characterization of the structure is strictly required, but, in general, it may present some critical issues. For example, the modal-based identification procedures are almost never completely and often not universal, but rather specifically related to the facility under examination. Moreover these techniques are not efficient in tracing the micro damages that could early occur in the process of damage, as illustrated by Peeters and Roeck [12]: modal parameters variations can be caused both by the presence of damage, whether as a result of changes in environmental conditions. Therefore making certain optimal number of sensors and placing sensors on appropriate locations are especially important.

But up to now, a few works were studied about the optimal number of needed sensors especially for data acquisition validity and micro damage detection. In fact, the effectiveness and the amount of information are equally important to the subsequent data analysis. The energy of wave (the strength of signal) will attenuate during the propagation. On the other hand, there are additional frequency component of wave at the damage position, the higher the frequency is, the rapid the wave attenuates. So even if the sensor got the signal, it cannot be judged if there is really damage when the energy of wave does not match the requirement of data acquisition system. The results of measurement are unreliable under this condition.

Reverberation matrix method is put forward in recent years. It is an effective way to analyze the transient dynamic response of structure based on wave propagation and classic Timoshenko beam theory [13, 14]. It has been proved that there is an obvious effect of inertia due to the shear deformation at high frequencies on the dynamic properties of structure [15]. Based on reverberation matrix method, a two-step method is hereby proposed in this paper to determine the sensors arrangement using the scattering matrix in this paper. First, calculate the utmost distance of wave propagation on a special structure by the principle of elastic wave propagation and determine the preliminary number of sensors; second, the utmost distance and number of sensors are applied to a sensor optimization algorithm named Effective Independence Driving-Point Residue method. In this way, the number and location of sensors can be optimized simultaneously.

2. Basic Theory

2.1. Analysis of Antinoise Capacity of Data Acquisition System

In actual engineering, different data acquisition system required different precision of the sensor based on the system’s antinoise capacity. It is needed to analyze the system’s antinoise capacity in order to ensure the validity and reliability of the measurement. The noise level that the system can tolerate must be estimated and peak signal to noise ratio (PSNR) must be calculated.

It is assumed that the given instrument’s measurement range is , precision is , so the absolute deviation threshold is . And the value of PSNR is

represents acceleration amplitude of signal, represents acceleration amplitude of noise.

According to (1), the needed signal minimum amplitude of acceleration can be calculated:

2.2. The Utmost Distance of Sensor Arrangement

The lateral motion equation of modification Timoshenko beam taking account of the shear displacement is:

where is section shear coefficient, is the area of section, is shear module, is density, is the moment of Inertia. , are lateral displacement caused by shear and moment, respectively, and the lateral displacement .

After Fourier transformed, the solutions of (3) are

the “ ” indicates a variant in the frequency domain, , , , and are undetermined constants depend on arrival wave, , , , and are undetermined constants depend on departure wave, is the Laplace transform complex variable, is 2.718 28, represents the number of wave, is the speed of wave, , , is frequency, is 0 or 1, and .

The moment and shear can be expressed as follows in the frequency domain:

So the solutions of (5) are

In reverberation matrix method, node scattering matrix is defined as , the departing wave amplitude matrix can be calculated by arriving wave amplitude matrix at node :

where and are local scattering matrix and local source matrix. The and are local arriving wave amplitude matrix and local departing wave amplitude matrix separately. For the whole structure, it is

where is global scattering matrix, is global source matrix, and and are global arriving wave amplitude matrix and global departing wave amplitude matrix. The global scattering matrix is a constant matrix depends on special structure. Equation (8) reflects the scattering relationship between departing wave amplitude and arriving wave amplitude under external force.

The in (6) represents the distance of wave propagation. The data acquisition system can correctly identify the signal only if the peak value of wave signal is greater than . The attenuation degree of peak value of response signal can be determined by matrix . When it satisfies , becomes the utmost distance of wave propagation, that is, the utmost distance of sensor arrangement.

2.3. Effective Independence Driving-Point Residue (EFI-DPR) Method

As mentioned above, previous works addressing the issue of optimally locating a given number of sensors in a structure have been carried out by several investigators. And the comparisons of different sensor placement methods have been done [2]. Because this paper mainly discusses the number of needed sensors taking into account data acquisition validity and micro damage detection sensitivity, so we just select a common method to optimize the location of sensors. In fact, a method named Effective Independence Driving-Point Residue (EFI-DPR) method is adopted in the aspect of location optimizing process.

The effective independence method (EFI method) was proposed by Kammer at 1991 [4]. The main idea of the EFI method is to gradually eliminate the degree of freedom which contributes the smallest to linearly independent of target mode vector from all measured points and to achieve the goal of that as much as possible modal parameters information is collected with limited number of sensors. The elimination process is implemented by the computing of effective independence distribution vector :

where represents the modal shape matrix of reduction finite element model (FEM) of the structure to be detected based on all candidate measurement points.

The other optimal sensor placement method is Kinetic energy method (KEM). The basic idea of KEM is that the node’s response is relatively large with larger modal strain energy, and then the sensors are installed on these nodes according to modal strain energy distribution map of every modal shape [16].

The step to optimize sensors location of KEM is similar to effective independence method (EFI), but the KEM method is highly dependent on the mesh level of finite element model. The sensor will be distributed more far away if the mesh is divided more coarsely, vice versa.

The definition of Kinetic energy matrix is

A residue coefficient of Effective Independence Driving-Point Residue (EFI-DPR) method is defined based on the theory of Kinetic energy:

where represents stiffness matrix of structure, represents unit stiffness modal kinetic energy, represents mass matrix, is frequency, and is the total number of nodes of structure.

To multiply (9) by , the independence distribution vector of EFI-DPR is shown in

The optimal process of EFI-DPR method is similar to the classic effective independence method. The performance of EFI-DPR becomes a little bit excellent through combining the advantages of the two traditional methods (EFI and KEM) [17]. (For more information, please see the reference [17].)

3. Process of Optimal Sensor Placement

3.1. The Optimal Minimum Number of Needed Sensors

First, utmost distance of sensor arrangement is calculated based on the theory of elastic wave propagation. Then the minimum number of sensors needed is estimated by , where is the longest size of the structure ( is not the size only in one dimension, it means the longest path of wave propagation through the whole structure.) and it is assumed that the sensors are arranged uniformly at first.

3.2. The Optimal Position of Sensors

First, the optimal location of sensors is obtained using Effective Independence Driving-Point Residue (EFI-DPR) method according to the . Then the distance between every two sensors will be checked after the sensors are placed by EFI-DPR. The number of sensors will be added if there is any distance longer than the ultimate distance mentioned above. This process will continue until all the distances are satisfied. The final optimal sensor placement plan will be get after the iteration is finished.

4. Case Study

4.1. Model Analysis and Utmost Distance Calculation

This paper considers the optimal sensor placement for a bridge benchmark model structure as shown in Figure 1. The length, width, and the height of model are 12 m, 2 m, and 1 m separately. All beams of the structure are joist steel, the elastic modulus is 2.0811 N/m², mass density is 7800 kg/m³, and Poisson’s ratio is 0.31.The parameters are shown in Table 1. All the joints are connected by bolts. The nine columns are fixed on the ground. The data acquisition system is Intelligent Signal Acquisition Analyzer INV306U-5160 produced by Beijing Orient Institute of Noise & Vibration. Its measurement range is ±10 m/s², precision is 0.001. According to (2), the needed minimum signal acceleration peak value is that  m/s² (the is −5 db according to the statistics result).

Table 1: Material parameters of structure.

2.0811
7800

3.8747
0.06

Figure 1: Benchmark model of bridge structure and its finite element model (FEM).

In fact, there are various frequency components of wave at the damage position due to different damage type. According to the results of simulation of the damage due to stiffness decrease, a wave with 2000 Hertz and 0.13 m/s² acceleration was taken as an example data. The attenuation curve of wave is shown in Figure 2.

Figure 2: Wave acceleration peak value attenuation.

The utmost distance of sensor arrangement was calculated using the proposed method. The results indicate that the distance of wave propagation is 5.86 meters when the peak value of acceleration attenuates to 4.3% of the origin (from 0.13 m/s² to 0.0056 m/s²). So the utmost distance of sensor arrangement is 5.86 meters.

4.2. Final Placement Plan

Based on the utmost distance  m, the optimal sensor placement of the bridge Benchmark model was accomplished using EFI-DPR method with the given data acquisition system. The result is shown in Figure 3. (The arrows in Figure 3 represent the position and direction of the sensors.)

Figure 3: Optimal sensor placement of Benchmark model.

4.3. Results Analysis

In order to check the validity of the proposed method, the values of mean square deviation (MSD) and modal assurance criteria (MAC) between FEM and the measured data from the installed sensors are calculated. Results are shown in Tables 2 and 3 and Figures 4 and 5.

Table 2: Mean square deviation (MSD) value between FEM modal shape and curve fitting of test.

Modal number
MSD

1

2

3

4

5

6

7

8

9

10

Table 3: Modal assurance criteria (MAC) value between FEM modal shape and curve fitting of test.

Modal number
MAC

1
0.9922

2
0.9468

3
0.9863

4
0.8633

5
0.8426

6
0.9771

7
0.9680

8
0.8544

9
0.8629

10
0.9534

Figure 4: Values of mean square deviation (MSD).

Figure 5: Values of modal assurance criteria (MAC).

Table 2 shows that the mean square deviation (MSD) value is very small at the first several modes. Table 3 shows the value of modal assurance criteria (MAC) is close to 1 at the first three modes, six-order mode, and seven-order mode. MAC is a very good common tool for evaluating the orthogonality of two vectors, and in general, the data of first several modes is enough for ordinary test.

At the same time, the mode shape comparison has been made between the simulation of FEM and the data of curve fitting based on the measured data from the installed sensors (the dashed line in Figure 6). Figure 6 shows that there is no significant difference between the mode shape of FEM and the one by curve fitting.

Figure 6: Modal shapes between finite element model (FEM) and cubic spline.

All the results show that the information measured from the installed sensors was consistent with the finite element model (FEM) at the first 1 to 6 modes, that is, the optimal sensor placement with proposed method is number enough and effective under the given data acquisition system.

5. Conclusions

This paper presents an alternative way to make effective optimal sensor placement especially for data acquisition validity and micro damage detection in structure health monitoring. Taking into account energy attenuate during the signal propagation, combined with classic reverberation matrix method, a two-step method is proposed to determine the sensors arrangement. The optimal number and location of sensors are calculated simultaneously through applying the utmost distance of sensors to a sensor optimization algorithm. The bridge benchmark model test is representative and the results show the validity of proposed method.

Although the study provides an alternative way for optimal sensor placement in structure damage detection and structure health monitoring, and this may make the system more economical and effective, there is still much more work need to do. The establishment of scattering matrix is more cumbersome. There are different for bar system and plate system. A little bit variation of exists even in the similar structure. On the other hand, more damage types must be considered for the selection of excitation wave in large structures. And how to determine the order of the monitoring modes still needs further research.

Acknowledgments

This research work was jointly supported by the National High Technology Research and Development Program (“863” Program) of China (Grant no. 2006AA04Z437) and the National Natural Science Foundation of China (Grant nos. 51278420 and 50878184). The authors would also like to thank the anonymous reviewers for their valuable suggestions.

References

A. Deraemaeker, E. Reynders, G. De Roeck, and J. Kullaa, “Vibration-based structural health monitoring using output-only measurements under changing environment,” Mechanical Systems and Signal Processing, vol. 22, no. 1, pp. 34–56, 2008.

T. H. Yi and H. N. Li, “Methodology developments in sensor placement for health monitoring of civil infrastructures,” International Journal of Distributed Sensor Networks, vol. 2012, Article ID 612726, 11 pages, 2012.

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.

D. C. Kammer, “Sensor placement for on-orbit modal identification and correlation of large space structures,” Journal of Guidance, Control, and Dynamics, vol. 14, no. 2, pp. 251–259, 1991.

K. Worden and A. P. Burrows, “Optimal sensor placement for fault detection,” Engineering Structures, vol. 23, no. 8, pp. 885–901, 2001.

T. H. Yi, H. N. Li, and M. Gu, “A new method for optimal selection of sensor location on a high-rise building using simplified finite element model,” Structural Engineering and Mechanics, vol. 37, no. 6, pp. 671–684, 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 Structure, vol. 21, no. 10, Article ID 05033, pp. 1–9, 2012.

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.

W. Liu, W. C. Gao, Y. Sun, and M. J. Xu, “Optimal sensor placement for spatial lattice structure based on genetic algorithms,” Journal of Sound and Vibration, vol. 317, no. 1-2, pp. 175–189, 2008.

A. Cruz, W. Vélez, and P. Thomson, “Optimal sensor placement for modal identification of structures using genetic algorithms-a case study: the olympic stadium in Cali, Colombia,” Annals of Operations Research, vol. 181, no. 1, pp. 769–781, 2010.

Y. J. Yan, L. Cheng, Z. Y. Wu, and L. H. Yam, “Development in vibration-based structural damage detection technique,” Mechanical Systems and Signal Processing, vol. 21, no. 5, pp. 2198–2211, 2007.

B. Peeters and G. D. Roeck, “One-year monitoring of the Z24-bridge: environmental effects versus damage events,” Earthquake Engineering and Structural Dynamics, vol. 30, no. 2, pp. 149–171, 2001.

Y. H. Pao, W. Q. Chen, and X. Y. Su, “The reverberation-ray matrix and transfer matrix analyses of unidirectional wave motion,” Wave Motion, vol. 44, no. 6, pp. 419–438, 2007.

Y. H. Pao, D. C. Keh, and S. M. Howard, “Dynamic response and wave propagation in plane trusses and frames,” AIAA Journal, vol. 37, no. 5, pp. 594–603, 1999.

R. Chen, C. F. Wan, and S. T. Xue, “Modification of motion equation of Timoshenko beam and its effect,” Journal of Tongji University, vol. 33, no. 6, pp. 711–715, 2005 (Chinese).

C. B. Larson, D. C. Zimmerman, and E. L. Marek, “A comparison of modal test planning techniques: excitation and sensor placement using the NASA 8-bay truss,” in Proceedings of the 12th International Modal Analysis Conference, vol. 2251 of Proceedings of SPIE, pp. 41–46, 1994.

Z. Wu, F. Dai, J. Song, and D. Wu, “More efficient optimal sensor placement method for structure damage detection,” Journal of Northwestern Polytechnical University, vol. 25, no. 4, pp. 503–507, 2007 (Chinese).