Next Article in Journal
Simultaneous Life Detection and Localization Using a Wideband Chaotic Signal with an Embedded Tone
Previous Article in Journal
Traffic Vehicle Counting in Jam Flow Conditions Using Low-Cost and Energy-Efficient Wireless Magnetic Sensors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Energy-Efficient Heterogeneous Wireless Sensor Deployment with Multiple Objectives for Structural Health Monitoring

1
Department of Civil and Environmental Engineering, Harbin Institute of Technology Shenzhen Graduate School, Shenzhen 518055, China
2
School of Engineering, San Francisco State University, San Francisco, CA 94132, USA
3
Department of Electronic and Information Engineering, Harbin Institute of Technology Shenzhen Graduate School, Shenzhen 518055, China
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(11), 1865; https://doi.org/10.3390/s16111865
Submission received: 12 September 2016 / Revised: 30 October 2016 / Accepted: 2 November 2016 / Published: 6 November 2016
(This article belongs to the Section Physical Sensors)

Abstract

:
Heterogeneous wireless sensor networks (HWSNs) are widely adopted in structural health monitoring systems due to their potential for implementing sophisticated algorithms by integrating a diverse set of devices and improving a network’s sensing performance. However, deploying such a HWSN is still in a challenge due to the heterogeneous nature of the data and the energy constraints of the network. To respond to these challenges, an optimal deployment framework in terms of both modal information quality and energy consumption is proposed in this study. This framework generates a multi-objective function aimed at maximizing the quality of the modal information identified from heterogeneous data while minimizing the consumption of energy within the network at the same time. Particle swarm optimization algorithm is then implemented to seek solutions to the function effectively. After laying out the proposed sensor-optimization framework, a methodology is present to determine the clustering of the sensors to further conserve energy. Finally, a numerical verification is performed on a four-span pre-stressed reinforced concrete box-girder bridge. Results show that a set of strategically positioned heterogeneous sensors can maintain a balanced trade-off between the modal information accuracy and energy consumption. It is also observed that an appropriate cluster-tree network topology can further achieve energy saving in HWSNs.

1. Introduction

Structural health monitoring (SHM) has been widely used to monitor and diagnose the health status of civil structures/infrastructures in real-time [1]. Typically, various sensors are deployed on critical locations of the structure to periodically collect different types of data relevant to the health status. For example, to monitor a bridge, engineers deploy sensors such as accelerometers, strain gauges, and displacement transducers. The collected data will be post-processed according to time/frequency domain algorithms to assess structural conditions, after which the SHM system will estimate the residual life of the structure and, possibly, send out alerts when it exceeds some pre-defined threshold.
Traditional SHM systems rely on wired sensors. The installation and maintenance cost of these sensors represents a large portion of the total cost of the system. With the advances in technology, the development of wireless sensor technology in recent years offers new opportunities for SHM applications [2,3,4,5]. Wireless communication eliminates the cost of wiring and increases the scalability of the SHM systems. Wireless sensor nodes are generally powered by batteries and mounted with multiple sensing devices placed inside and around the structure to periodically collect various types of data. Consequently, wireless sensor networks (WSNs) particularly designed for SHM purpose are evolving into heterogeneous systems [6]. Among various advantages, for example, heterogeneous wireless sensor networks (HWSNs) can potentially provide higher levels of computational power, network densities and lifetimes [7].
As a large amount of existing structures are suffering from deterioration, excessive loadings and unpredicted incidents, the development of dense, yet low-cost wireless sensor arrays to monitor the operation of structural systems has become an attractive research subject in recent years. Some notable experimental applications implemented on large-scale structures can be seen, for instance, in the Golden Gate Bridge [8] and the Jindo Bridge [9,10] as well as in several other long-span [11,12] and medium-length span bridges [13,14]. However, WSNs are not yet commonly adopted in permanent monitoring systems largely due to the challenges of providing necessary power for the monitoring system through the limited power supply of sensor nodes with relatively small size. Thus, WSNs are still considered as less popular alternatives to conventional wired SHM systems. Therefore, maintaining an optimal size of sensor network that can provide the desired modal information is deemed as one of the most critical challenges in the current deployment practice of WSNs in SHM.
In this paper, a study attempting to implement optimal sensor deployment (OSD) techniques to determine the number and locations of heterogeneous wireless sensors on a structure of interest has been performed. This implementation of such system intends to produce reliable information about the structure’s health status while maintaining reasonable energy consumption.
As an important issue in SHM, the OSD techniques have been demonstrated to be effective in helping to achieve accurate estimation of modal parameters by placing the sensors properly. Within the context of OSD, quite a few methods, such as the Effective Influence (EI) method [15], Modified Variance (MV) method [16], and Kinetic Energy (KE) method [17], have been proposed and verified by practical SHM implementations. In addition to these existing OSD techniques designed for modal parameter estimation in wired SHM systems, a handful of works have attempted to address the OSD issues for wireless sensor networks. For instance, Bhuiyan et al. [18,19] proposed a three-phase sensor placement method for SHM that addressed the quality of sensor placements, communication efficiency, and fault tolerance. Onoufriou et al. [20] presented a two-step strategy to optimize the number of sensors and their locations to satisfy both specific structural engineering requirements and energy constraint imposed by a WSN. Zhou et al. [21] formulated an energy-aware wireless sensor placement framework and developed a hybrid discrete firefly algorithm to solve complex optimization problem. Fu et al. [22] performed a study to optimize wireless sensor placement for SHM in terms of both the quality of the modal information and network energy consumption. These OSD algorithms were feasible because most WSNs in these studies were homogeneous (sensors have the same type, storage, processing, battery power, sensing, and communication capabilities). However, when applying OSD to a HWSN that supports multi-type sensor applications, the OSD should further consider the following issues:
i. Placement performance metrics issue
In a multi-type sensor network, the orthogonality of the modal vectors cannot be exploited as in single-type sensor networks. Therefore, approaches that do not rely on the orthogonality of the modal vectors are needed.
ii. Network topology issue
Network heterogeneity allows for sophisticated operations over a larger region. As most of the energy costs come from the transmission and receiving of data packets, network topology optimization is a challenging issue for energy efficient coordination when WSNs become heterogeneous.
To the best of the authors’ knowledge, none of the existing works attempt to fulfill all the above objectives simultaneously. By employing the ratio of Modal Clarity Index (MCI) [23] and Mode Shape Expansion (MSE) [24], Jalsan et al. [25] took the information quality of the measured data collected from a HWSN into account to represent placement quality of strain and acceleration sensors. However, their work did not consider the allocation of data packets along the end-to-end path from a sensor to the base station, and therefore the optimization within HWSNs could be further improved.
In this paper, a framework is proposed to overcome the insufficiency of existing OSD techniques for HWSNs with the consideration of both placement quality and clustering issues. Typical heterogeneous devices (accelerometers and strain gauges) are considered herein since they are commonly available and have proven to be effective in detecting changes of structural properties, particularly in civil structures with low vibration amplitudes. The proposed approaches are in twofold. Firstly, a framework is proposed to optimize the layout of multi-type wireless sensors at critical locations of the structure in terms of the modal information quality and network energy consumption. Secondly, a clustering algorithm is further proposed to determine the best number of clusters in the HWSN. After the framework is laid out, a numerical verification is performed on a four-span pre-stressed reinforced concrete box-girder bridge.
The remaining sections of the paper are organized as follows: Section 2 presents the formulation of the optimization problem in HWSNs. Section 3 describes the multi-objective function and its optimization. Section 4 discusses the clustering issues in HWSN in terms of communication cost. Section 5 evaluates the performance of the proposed approach via a numerical simulation of a Heterogeneous wireless SHM system deployed on a bridge. Conclusions are drawn in Section 6.

2. Optimization Problem Formulation in HWSNs

The optimization for HWSN configuration is aimed to select particular physical locations for the sensors on the monitored structure such that the resulting HWSN exhibits optimal performance. To be specific, the optimization goal is set to find node locations in order to reliably identify modal information (basis for diagnosing the health of a structure) while consuming minimum energy during data collection. The formulation of optimization goal is described in the following sections.

2.1. Heterogeneous Sensor Placement Quality

In this study, we focus on modal identification performance based on vibration information, specifically acquired through two common types of sensors in SHM—Wireless accelerometers and strain gauges. The assessment of sensor placement quality based on modal identification accuracy has been previously studied in homogeneous WSNs. For example, the Modal Assurance Criterion (MAC) which considers the orthogonality of the modal vectors was used to evaluate the sensor placement [26]. However, in the case of heterogeneous data consisting of acceleration and strain, the orthogonality of the modal vectors is no longer available. Therefore, performance metrics such as the Modal Clarity Index (MCI) and Modal Relative Error (MRE) that do not rely on the orthogonality of the modal vectors need to be employed. A few heterogeneous sensor placement indices based on these metrics have been proposed in the literature, including the weighted sum of MCI and MRE [27], and the ratio of MCI and MRE [25]. In this work, the ratio of MCI and MRE is employed as the placement quality index to quantify the quality of the sensor placement.

2.1.1. Modal Clarity Index

The objective of MCI is to determine the sensor locations that maximize the clarity between adjacent modes of response. The MCI is based on the least squares method. The best-fit amplitude matrix l is constructed using Equation (1):
λ p , q = i = 1 N a + N s α i , p β i , q i = 1 N a + N s α i , p 2 ,
where p and q are the modes being compared, a and b are the analytical and experimental modal matrices comprising of strain measurements and acceleration measurements, respectively.
The MCI can be obtained as the difference between the excited mode p and the best fit mode q using Equation (2):
M C I p , q = [ β p ( λ p , q α q ) ] T [ β p ( λ p , q α q ) ] ,
The MCI matrix is a square matrix with dimensions equal to the modes of interest. For simplification, mean value of the sum of matrix elements are calculated as the IMCI:
I M C I ( n a , n s ) = 1 m 2 p = 1 m q = 1 m I M C I p , q ,
where m is the number of measured modes, Na is the number of acceleration nodes deployed, and Ns is the number of strain nodes deployed.

2.1.2. Modal Relative Error

In general, it is impossible to identify all the vibration modes of a structure from measured data. Therefore, modal shape expansion (MSE) can be used to estimate the response of the structure at the Degrees of Freedoms (DOFs) where no sensor is equipped, from a limited number of measured DOFs [27].
The MSE method can be expressed as
y α = α [ ( β T β ) 1 β ] y β ,
in which ya is the estimated response, a is the system modal matrix, b is the modal matrix for the measured DOFs, and yb is the measured response of the structure.
The MRE can then be calculated using the estimated response and measured response via Equation (5):
I M R E = | y α y β | | y β | ,

2.1.3. Placement Quality Index

The larger the corresponding MRE, the clearer the calculated modal information would be (and vice versa) [25]. To consider both effects of MRE and MCI, the ratio of the two metrics is defined as the placement quality index to measure the modal identification accuracy [9]:
F M = I M R E I M C I ,
The ratio of IMRE and IMCI represents a IMRE gain per unit IMCI value. Minimizing Equation (6) yields the best combination of both measurement types. This formulation overcomes the orthogonality issue arising in heterogeneous sensors deployment and thus is used in this study.

2.2. Network Model

Generally, energy consumption of a typical wireless sensor node in operation mainly occurs in the phases of data acquisition, data processing, data reception, and data transmission [28]. In this study, all of these phases are considered in the energy consumption formulation for acceleration node and strain node separately. Table 1 shows the formulations of energy consumption of each part for different sensor nodes. Table 2 lists the nomenclature of parameters and their typical values used in general WSN platforms for SHM applications [29] and the following equations in this study.
In order to evaluate the energy consumption of HWSNs composed of acceleration and strain nodes, a multi-layer spherical network model with a base station located in the center as shown in Figure 1 is used. Ideal communication without errors and delays is assumed. This network model, presented in the authors’ previous work [29], is now extended with heterogeneous sensor nodes deployed. If the nodes are evenly distributed in the network (in other word, obey a uniform distribution), the node density in the spherical space for acceleration pa and strain ps can be defined as follows:
p a = N a V a , p s = N s V s ,
where Va and Vs are the volume of space covered by the acceleration nodes and strain nodes respectively. Na and Ns have the same definition to the corresponding parameters used in Equation (1).
When considering the pre-defined sensing and routing function of sensor nodes, a multi-layer HWSN contains two main portions of energy consumption. The first portion is the energy consumed by a layer of sensor nodes due to data acquisition, processing, and information transition to the upper level. The other portion is the energy consumed by the same layer of sensor nodes for receiving data from the lower layer and passing them to the upper layer.
In the first portion, the energy consumed is estimated as:
F E ( n , n ) = i = 1 N n a F n o d e ( V i ) a + i = 1 N n s F n o d e ( V i ) s = i = 1 N n a { α a K + β a ( d a ) m K + [ γ s + δ s ( d s ) m ] K } + i = 1 N n s { α s K + β s ( d i ) m K + [ γ s + δ s ( d i ) m ] K } ,
where Nna is the number of acceleration nodes in the nth layer; and Nns is the number of strain nodes in the nth layer.
Assuming the sensor nodes are uniformly deployed in the network, the number of sensor nodes in each layer can be calculated as:
N n a = ( 3 n 2 3 n + 1 ) R a ' 3 R a 3 N a ,   N n s = ( 3 n 2 3 n + 1 ) R s ' 3 R s 3 N s ,
where Ra and Rs are the distances between layers deployed with acceleration nodes and strain nodes, respectively; R’a and R’s are the coverage ranges of all acceleration nodes and strain nodes in the network, respectively.
The expectation of data transmission distance in this process is:
E ( d i ) = i = 1 N n d i p i = ( n 1 ) R n R [ x ( n 1 ) R ] 4 π x 2 4 π { ( n R ) 3 [ ( n 1 ) R ] 3 } / 3 d x = 3 2 n 2 n + 1 4 3 n 2 3 n + 1 R ,
Similarly, the energy consumed in the nth layer in the second portion is calculated by:
F E ( n , n 1 ) = i = 1 N n a ( γ a + δ a d a m + ε a ) K + i = 1 N n s ( γ s + δ s d s m + ε s ) K = ( 3 n 2 3 n + 1 ) R a ' 3 R a 3 N a [ γ a + δ a ( 3 2 n 2 n + 1 4 3 n 2 3 n + 1 R a ' ) m + ε a ] K + ( 3 n 2 3 n + 1 ) R s ' 3 R s 3 N s [ γ s + δ s ( 3 2 n 2 n + 1 4 3 n 2 3 n + 1 R s ' ) m + ε s ] K ,
Note that the data collected in the nth layer data needs (n − 1) hops to reach the base station. By summarizing these two portion of energy consumed, the HWSN energy consumption per unit time can be expressed by:
F E = n = 1 M F E ( n , n ) + n = 2 M ( n 1 ) F E ( n , n 1 ) = n = 1 M { ( 3 n 2 3 n + 1 ) R a ' 3 R a 3 N a { α a K + β a ( 3 2 n 2 n + 1 4 3 n 2 3 n + 1 R a ' ) m K + [ γ a + δ a ( 3 2 n 2 n + 1 4 3 n 2 3 n + 1 R a ' ) m ] K } + ( 3 n 2 3 n + 1 ) R s ' 3 R s 3 N s { α s K + β s ( 3 2 n 2 n + 1 4 3 n 2 3 n + 1 R s ' ) m K + [ γ s + δ s ( 3 2 n 2 n + 1 4 3 n 2 3 n + 1 R s ' ) m ] K } } + n = 2 M { ( n 1 ) ( 3 n 2 3 n + 1 ) R a ' 3 R a 3 N a [ γ a + δ a ( 3 2 n 2 n + 1 4 3 n 2 3 n + 1 R a ' ) m + ε a ] K + ( n 1 ) ( 3 n 2 3 n + 1 ) R s ' 3 R s 3 N s [ γ s + δ s ( 3 2 n 2 n + 1 4 3 n 2 3 n + 1 R s ' ) m + ε s ] K } ,
In this study, we assume the acceleration and strain nodes conform to the same deployment rules, that is:
R a ' = R s ' = R , R a = R s = M R ,
By taking Equation (13) and its typical values of the corresponding parameters into Equation (12), the energy consumption formulation for a HWSN is obtained as Equation (14):
F E = n = 1 M { 3 n 2 3 n + 1 M 3 × ( 3 2 n 2 n + 1 4 3 n 2 3 n + 1 R ) 2 × 1.36 × 10 9 × ( N a + N s ) + 3 n 2 3 n + 1 M 3 × ( 1.05 × 10 7 N a + 9 × 10 8 N s ) } K + n = 2 M { ( n 1 ) 3 n 2 3 n + 1 M 3 [ 1.8 × 10 7 + 10 11 × ( 3 2 n 2 n + 1 4 3 n 2 3 n + 1 R ) 2 ] × ( N a + N s ) } K ,

3. Objective Function and Solution

Trade-offs exist between modal identification accuracy and energy consumption in WSNs for SHM [21,30]. Therefore, the objective of the optimization problem herein is to determine the optimal number and location of heterogeneous sensor nodes in order to balance the amount of energy consumed and the modal identification quality in the HWSNs. The following sections describe the formulation of the multi-objective optimization function and its resolving algorithm.

3.1. Multi-Objective Optimization Function

A multi-objective function is formulated to describe the optimization problem as follows:
F ( F M , F E ) = ϕ M ω M F M ( N a , N s , R a , R s ) + ϕ E ω E F E ( N a , N s , R a , R s ) ,
where fM is the weight coefficient of modal identification, fE is the weight coefficient of energy consumption, wM is the adjustment coefficient of modal identification index, and wE is the adjustment coefficient of energy consumption. Note that each single objective function (FM and FE) has its own numerical dimension, thus adjustment coefficients are utilized to normalize the final calculation results of the multi-object function.
As discussed in Section 2.1, the ratio of MSE and MCI is used to formulate the objective function FM for the sensor placement quality (see Equation (6)). Meanwhile, Equation (14) in Section 2.2 is used to formulate the objective function FE for energy consumption.

3.2. Particle Swarm Optimization Algorithm

The optimization of the above multi-objective function for large SHM systems is computationally challenging due to the large number of DOFs. In recent years, multiple advanced optimization algorithms, such as the genetic algorithm (GA) [31], particle swarm optimization (PSO) [32], and the monkey algorithm [33], have been adopted to determine sensor optimal placements. Unlike other computation algorithms, PSO has no evolution operators, such as crossover and mutation, and therefore it has a faster convergence speed. In this study, the PSO algorithm is employed also due to its unique information diffusion capability and interaction mechanisms which enable PSO to solve complex optimization problem with high performance at low computational cost. Refer to Section 5.2 for the details of the utilization of PSO algorithm in this study.

4. Network Topology Optimization

HWSNs applications, such as SHM, require the equipped sensors to form a multi-hop network to collect the environmental data in real-time. Such a network typically generates a cluster-tree type topology. As most of the energy cost comes from transmitting and receiving data packets, clustering optimization become an important issue for these applications in that it benefits efficient energy coordination, thus reducing the transmission power and maximizing the network life. To this end, the topology optimization problem is considered in this section to determine an appropriate number of clusters to minimize the energy consumption after selecting the sensor number and location.
In this study, the energy consumption model for a single sensor node is assumed to conform to the first order wireless communication mode [34] as shown in Figure 2. Note that the acceleration and strain nodes have the same radio frequency module, therefore these two types of sensors share the same communication mode.
In the data transmission process, the energy consumption of transmitting K bits of data to a sensor node can be expressed as:
E T X = { K × γ + K × δ m p × d 4 K × γ + K × δ f s × d 2 ,
The energy consumption of receiving K bits of data from a sensor node can be expressed as
E R X = K × ε ,
where K is the number of data packet being sent and received; γ is the energy needed to transmit a unit bit of data; ε is the energy consumption of receiving a unit bit of data; δmp is the power amplification factor of multi-path attenuation model (d ≥ d0); δfs is the power amplification factor of the free space model (d < d0); d 0 = d f s / d m p ; d is the data transmission distance.

4.1. Two-Phase Energy Consumption Formulation

Consider a synchronized two-level cluster-tree HWSN featuring a tree-based logical topology where nodes are organized in different groups, called clusters. Each member node only interacts with its pre-defined cluster head (CH) node, and cannot be connected to other member nodes. The formulation of energy consumption in such a cluster-tree HWSN can be divided into two phases, the network initialization phase and the stable transmission phase.
Assuming that acceleration nodes, Na, and strain nodes, Ns, are randomly distributed in a L × L × L cube space, in which the number of CH nodes is C. Then each cluster has one CH node and (Na + Ns)/C − 1 member nodes on average. The energy consumption of CH nodes broadcasting message to member nodes is given by:
E c h 11 = K × γ + K × δ m p × d 4 ,
The determination of what cluster for a non-CH node to join is based on the strength of the signal received by the broadcast message, which is proportional to the distance d between the member nodes and the CH node. The energy consumption of each CH node receiving information from member nodes is given by:
E c h 12 = K × ε × ( N a + N s C 1 ) ,
After receiving the information, the CH node creates a Time Division Multiplexing Access (TDMA) [24] to be sent back to the member nodes, which also consumes energy. If the distance between the CH node and the member node is dtoch, the energy consumption in CH is given by:
E c h 13 = K × γ + K × δ f s × d t o c h 2 ,
According to Equations (18)–(20), the energy consumption of each CH node is given by:
E c h 1 = [ ε × ( N a + N s C 1 ) + 2 γ + δ f s × d t o c h 2 + δ m p × d 4 ] × K ,
Similarly, the energy consumption of a member node receiving broadcast information from the CH node is given by:
E m n 11 = K × ε ,
The energy consumption of a member node sending a message to join the cluster is given by:
E m n 12 = K × γ + K × δ f s × d t o c h 2 ,
The energy consumption of a member node receiving TDMA time table from the CH node is given by:
E m n 13 = K × ε ,
Therefore, the energy consumption of each member node in an established cluster is given by:
E m n 1 = E m n 11 + E m n 12 + E m n 13 = ( 2 ε + γ + δ f s × d t o c h 2 ) × K ,
The total energy consumption for a cluster in the network initialization phase is given by:
E 1 = C × [ E c h 1 + ( N a + N s C 1 ) × E m n 1 ] C × E c h 1 + ( N a + N s ) × E m n 1 ,
After the clusters are established, the network turns into the stable phase of data transmission. In each cluster, the member nodes send the information to the CH node, and the energy consumption of the CH node receiving this information is given by:
E c h 21 = K × ε × ( N a + N s C 1 ) ,
Then the CH node receives the information of the member nodes for information fusion. The energy consumption for such activity is given by:
E c h 22 = K × E D A × N a + N s C ,
where EDA is the energy consumed for fusing a unit bit of information.
Defining the distance between the CH node and base station as dtoBs, the energy consumption that CH node transmits the fused information to the base station is given by:
E c h 23 = K × γ + K × δ f s × d t o B s 2 ,
As a result, the energy consumption of a CH node in the stable phase is given by:
E c h 2 = E c h 21 + E c h 22 + E c h 23 = [ ε × ( N a + N s C 1 ) + γ + δ f s × d t o B s 2 + E D A × N a + N s C ] × K ,
Each member node only sends information to the CH node. Its energy consumption is:
E m n 2 = K × γ + K × δ f s × d t o c h 2 ,
Summarizing Equations (30) and (31), the energy consumption in the stable transmission phase is given by:
E 2 = C × [ E c h 2 + ( N a + N s C 1 ) × E n o c h 2 ] C × E c h 2 + ( N a + N s ) × E m n 2 ,
Consequently, the total energy consumption for a HWSN with two-level cluster-tree is shown in Equation (33):
E T o t a l = K × { ε × [ 4 ( N a + N s ) 2 C ] + [ 2 ( N a + N s ) + C ] × δ f s × d t o c h 2 + γ × [ 3 C + 2 ( N a + N s ) ] + C × δ m p × d 4 + C × δ f s × d t o B s 2 + E D A × ( N a + N s ) } ,

4.2. Optimal Clustering

In Equation (33), dtoch is an uncertain value since the number of CH nodes is unknown. Therefore, we use the mathematical expectation to estimate the value of dtoch.
Assuming that clusters composed of a CH node and multiple member nodes in a sphere volume are evenly distributed in the L × L × L cube space. In this case, for each cluster, the average volume is L3/C, the coverage radius is 3 L 3 / 4 pC 3 , and the probability density of member nodes is C/L3. The expectation of dtoch is:
E [ d t o c h 2 ] = ( x 2 + y 2 + z 2 ) ρ ( x , y , z ) d x d y d z = r 2 ρ ( r , φ , θ ) r 2 sin φ d r d φ d θ = C L 3 r 4 sin φ d r d φ d θ = 0.25 ( L 3 π C ) 2 / 3 0.116 L 2 C ,
With partial derivative on the number of CH nodes on both sides of Equation (34), it gives
E T o t a l C = 3 γ 2 ε 0.233 L 2 ( N a + N s ) δ f s C 2 + δ f s d t o B s 2 + δ m p d 2 ,
After setting E T o t a l C = 0 , it would now be possible to obtain the mathematical expression of the optimal number of CH nodes as
C = 0.233 L 2 ( N a + N s ) δ f s 3 γ 2 ε + δ f s d t o B s 2 + δ m p d 4 ,
As can be seen in Equation (36), the optimal number of CH nodes in a HWSN is determined by the total number of sensor nodes, the distance d that CH nodes broadcast information to their connected member nodes, and the distance dtoBs between the CH nodes and the base station.

5. Performance Evaluation

A numerical verification is performed on a four-span continuous reinforced concrete box-girder bridge. The optimization goal is to obtain a HWSN configuration that simultaneously minimize the network energy consumption and maximize the information quality objectives. In addition, this example also intends to show the efficiency of the proposed framework by comparing the performance of the cluster-tree topology to flat topology, which assumes that data packets are transmitted from member nodes to base station directly.

5.1. Simulation Setup

The bridge is located in Shenzhen, China, across the Pinghu railway. It is a historical bridge currently being retrofitted. A SHM system composed of wireless acceleration and strain nodes will be implemented on the bridge to monitor its long-term vibration behaviors. The bridge has span length of 42.5 m + 2 × 65 m + 42.5 m as seen from the bridge elevation view in Figure 3.
Bearing conditions of the bridge on a plan view are shown in Figure 4, where circles represent the double-column piers, DX represents an expansion bearing providing the bridge vertical constraint only, ZX represents an expansion bearing providing the bridge lateral and vertical constraints, and GX represents a fixed bearing.
To obtain modal parameters of the bridge, a detailed finite element model was built in ANSYS. The main box girder is modeled with 9854 solid elements and 3354 nodes. In this study, the first 10 mode shapes of the main girder (exhibited in lateral and vertical modes) are considered as the targeted modes and obtained through modal analysis. The first ten modes are shown in Figure 5 with their natural frequencies are listed in Table 3.

5.2. Sensor Layout Optimization Results

The proposed OSD techniques are applied to determine both the optimal number of heterogeneous sensors and their best sensing locations. In order to solve the multi-objective function established in Section 5.1, it is assumed that the accuracy of the modal parameter identification and network energy consumption for this SHM system has equal significance, meaning that the weight coefficients fM and fE in Equation (15) are both set to be 0.5. Since the two objective values are not in the same scale of magnitude, scale factors wR = 1 and wE = 106 are used for normalization. The MathWorks MATLAB is used to apply the PSO algorithm and carry out the calculations. The PSO parameters of initial population size, maximum number of iterations, and particle velocity, are set to be 300, 150, and [−3, 3], respectively. The PSO is run for 10 times to generate the network layouts. The objective function values, accelerometers, strain gauges, and coverage ranges of the layouts obtained in the ten cases are listed in Table 4.
Herein, the mean values are used as the optimal sensor layout deployed on the bridge plan. The optimal number of accelerometers, strain gauges, and coverage range are 14, 10, and 96.7 m, respectively. The sensor placement layout found from the analysis is shown in Figure 6, where black circles indicate wireless accelerometers and red rectangles represent the wireless strain gauges, respectively. Numbers 1~22 listed under the bottom of the bridge represent selected cross sections that are equipped with sensors. Numbers listed above the bridge represent longitudinal distance (cm) between sensors. The results showed non-symmetric sensor configurations in the longitudinal direction. Detailed sensor locations on each section are shown subsequently in Figure 7.

5.3. Clustering Optimization Results

The study herein is aimed to determine the best data transmission path of the HWSN based on the sensor configuration from initial optimization results, and to compare the performance of network topology effects on energy consumption. The related network parameters on the bridge monitoring system are listed in Table 5.
Since the probability density ρ(x, y, z) of the CH nodes is the reverse of the HWSN space V, the distance expectation dtoBs is calculated using Equation (35):
E [ d t o B s 2 ] = [ ( x 107.5 ) 2 + y 2 + ( z 6.875 ) 2 ] ρ ( x , y , z ) d x d y d z = 1 V [ ( x 107.5 ) 2 + y 2 + ( z 6.875 ) 2 ] d x d y d z = 7.8 × 10 3 m 2 ,
Subsequently, taking these parameter values into Equation (36), the optimal number of clusters is obtained as five. That is, the wireless sensors selected to monitor the bridge are divided into five clusters. Nevertheless, finding specific member nodes that result in best data transmission path in a large-scale HWSN with multiple clusters is a quite complex problem. In this study, for fast optimization, the nodes close to the base station are assumed to be selected as CH nodes and the member nodes in each cluster are allocated as even as possible. Figure 8 shows the cluster-tree based HWSN topology on the bridge in a 2-D configuration. Each cluster has four or five member nodes. The dotted line indicates the data transmission path.
In order to evaluate the effectiveness of clustering optimization, the simulation results using the proposed method are compared to those obtained by a flat network topology (one layer). According to Equation (14) derived in Section 2.2, the network energy consumption for both scenarios is calculated and shown in Table 6.
In the data transmission path produced by the cluster-tree HWSN topology, the coverage radius is increased by 28.3% compared to that of the flat topology. In addition, it is also found that the energy consumption of the sensors is significantly lower than that of the flat topology (33.3%). This reduction is expected since, although the coverage radius increases, the transmission distance of data packets is reduced in a cluster-tree routing path, which indicates that the energy consumption depends not only on the coverage radius of network, but also other factors such as the data transmission distance and the amount of transmitted data packets. Since there exists a linear relationship between energy consumption and the square of data transmission distance, the energy saving in a cluster-tree topology can be achieved by reducing total data transmission distance through suitable clustering mechanism, and, therefore, it is independent of the network size.
It should be emphasized that the optimization of network topology in this study only involves a fast selection of CH nodes according to the algorithm proposed. Further determination of specific member nodes in each cluster will lead to another complex problem which is out of the scope of this paper. In the expense of larger amount of computation time, the data packets transmission path shown in this study could be further optimized. However, the OSD framework proposed is applicable to other monitoring scenarios that the deployment of multi-type wireless sensors is based on the principle of modal information quality and minimum of energy consumption.

6. Conclusions

The increasing interest in employing HWSNs for SHM applications requires an efficient sensor placement methodology. In this study, first, an OSD framework for HWSNs is proposed. A multi-objective layout optimization problem is presented and resolved using the PSO algorithm to determine the trade-off between modal information quality and energy consumption with the consideration of the application requirements. Furthermore, a clustering optimization approach is proposed to conserve extra energy through selecting appropriate CH nodes in multi-hop HWSNs. A continuous-span bridge is used as an example to evaluate the performance of the proposed approach for designing a wireless SHM system comprised of acceleration nodes and strain nodes on the bridge.
The proposed optimization approaches can effectively determine the best number and location of heterogeneous wireless sensors for the purpose of SHM. Simulation results demonstrate that a set of strategically positioned heterogeneous sensors can maintain an optimal balance between the modal information accuracy and energy consumption. With the determined sensor deployment configuration, the proposed clustering optimization approach can be easily implemented to further conserve energy by selecting appropriate CH nodes in a multi-hop HWSN.

Acknowledgments

This work is financially supported by the National Science & Technology Support Program of China (No. 2014BAL05B06), Guangdong International Co-op Project (No. 2014A050503016), Shenzhen Science and Technology Projects (No. JCYJ20160525163140206, No. GJHZ20150312114346635), and the authors are grateful to the authorities for their support.

Author Contributions

Chengyin Liu proposed the optimization framework and techniques; Zhaoshuo Jiang contributed to the performance evaluation of the proposed method; Fei Wang provided constructive suggestions on potential problems about WSNs; Hui Chen analyzed the bridge in simulation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Farrar, C.R.; Worden, K. An introduction to structural health monitoring. Phil. Trans. R. Soc. A 2007, 365, 303–315. [Google Scholar] [CrossRef]
  2. Feltrin, G.; Jalsan, K.E.; Flouri, K. Vibration monitoring of a footbridge with a wireless sensor network. J. Vib. Control 2013, 19, 2285–2300. [Google Scholar] [CrossRef]
  3. Casciati, S.; Faravelli, L.; Chen, Z. Energy harvesting and power management of wireless sensors for structural control applications in civil engineering. Smart Struct. Syst. 2012, 10, 299–312. [Google Scholar] [CrossRef]
  4. Spencer, F.B.; Ruiz-Sandoval, M.; Kurata, N. Smart sensing technology: Opportunities and challenges. Struct. Control Health Monit. 2004, 11, 349–368. [Google Scholar] [CrossRef]
  5. Jian, L.; Mechitov, A.K.; Kim, R.E.; Spencer, F.B. Efficient time synchronization for structural health monitoring using wireless smart sensor networks. Struct. Control Health Monit. 2016, 23, 470–486. [Google Scholar]
  6. Chin, J.C.; Rautenberg, J.M.; Ma, Y.T.; Pujol, S.; Yau, K.Y. An Experimental Low-Cost, Low-Data-Rate Rapid Structural Assessment Network. IEEE Sens. 2009, 9, 11. [Google Scholar]
  7. Fok, C.L.; Roman, G.C.; Lu, C.Y. Servilla: A flexible service provisioning middleware for heterogeneous sensor networks. Sci. Comput. Program. 2012, 77, 663–684. [Google Scholar] [CrossRef]
  8. Kim, S.; Pakzad, S.; Culler, D.; Demmel, J.; Fenves, G.; Glaser, S.; Turon, M. Health monitoring of civil infrastructures using wireless sensor networks. In Proceedings of the Sixth International Symposium on Information Processing in Sensor Networks, Cambridge, MA, USA, 25–27 April 2007; pp. 254–263.
  9. Jo, H.; Sim, S.H.; Mechitov, K.A.; Kim, R.; Li, J.; Moinzadeh, P.; Spencer, J.B.; Park, J.W.; Cho, S.; Jung, H.J.; et al. Hybrid wireless smart sensor network for full-scale structural health monitoring of a cable-stayed bridge. SPIE Proc. 2011, 7981. [Google Scholar] [CrossRef]
  10. Shinae, J.; Jo, H.; Cho, S.; Mechitov, K.; Rice, A.R.; Sim, S.-H.; Jung, H.-J.; Yun, C.-B.; Spencer, B.F.; Agha, G. Structural health monitoring of a cable-stayed bridge using smart sensor technology: Deployment and evaluation. Smart Struct. Syst. 2010, 6, 439–459. [Google Scholar]
  11. Weng, J.H.; Loh, C.H.; Lynch, J.P.; Lu, K.C.; Lin, P.Y.; Wang, Y. Output-only modal identification of a cable-stayed bridge using wireless monitoring systems. Eng. Struct. 2008, 30, 1820–1830. [Google Scholar] [CrossRef]
  12. Chae, M.J.; Yoo, H.S.; Kim, J.Y.; Cho, M.Y. Development of a wireless sensor network system for suspension bridge health monitoring. Autom. Constr. 2012, 21, 237–252. [Google Scholar] [CrossRef]
  13. Whelan, M.J.; Gangone, M.V.; Janoyan, K.D.; Jha, R. Real-time wireless vibration monitoring for operational modal analysis of an integral abutment highway bridge. Eng. Struct. 2009, 31, 2224–2235. [Google Scholar] [CrossRef]
  14. Whelan, M.J.; Gangone, M.V.; Janoyan, K.D.; Jha, R. Operational modal analysis of a multi-span skew bridge using real-time wireless sensor networks. J. Vib. Control 2011, 17, 1952–1963. [Google Scholar] [CrossRef]
  15. Kammer, D.C. Sensor placement for on-orbit modal identification and correlation of large space structures. In Proceedings of the IEEE American Control Conference, San Diego, CA, USA, 23–25 May 1990; pp. 2984–2990.
  16. Minwoo, C.; Shamim, N.P. Optimal Sensor Placement for Modal Identification of Bridge Systems Considering Number of Sensing Nodes. J. Bridge Eng. 2014, 19, 04014019. [Google Scholar]
  17. Heo, G.; Wang, M.L.; Satpathi, D. Optimal transducer placement for health monitoring of long span bridge. Soil Dyn. Earthq. Eng. 1997, 16, 495–502. [Google Scholar] [CrossRef]
  18. Bhuiyan, M.Z.A.; Wang, G.; Cao, J.N. Sensor deployment with multiple objectives for structural health monitoring. ACM Trans. Sens. Netw. 2014, 10, 68. [Google Scholar] [CrossRef]
  19. Bhuiyan, M.Z.A.; Cao, J.N. Deploying Wireless Sensor Networks with Fault-Tolerance for Structural Health Monitoring. IEEE Trans. Comput. 2015, 64, 382–395. [Google Scholar] [CrossRef]
  20. Onoufriou, T.; Soman, R.N.; Votsis, R.; Chrysostomou, C.; Kyriakides, M. Optimization of wireless sensor locations for SHM based on application demands and networking limitations. In Proceedings of the 6th International IABMAS Conference, Stresa, Italy, 8–12 July 2012.
  21. Zhou, G.D.; Yi, T.H.; Zhang, H.; Li, H.N. Energy-aware wireless sensor placement in structural health monitoring using hybrid discrete firefly algorithm. Struct. Control Health Monit. 2015, 22, 648–666. [Google Scholar] [CrossRef]
  22. Fu, T.S.; Ghosh, A.; Johnson, E.A.; Krishnamachari, B. Energy-efficient deployment strategies in structural health monitoring using wireless sensor networks. Struct. Control Health Monit. 2013, 20, 971–986. [Google Scholar] [CrossRef]
  23. Jalsan, K.E.; Soman, R.N.; Flouri, K.; Kyriakides, M.A.; Feltrin, G.; Onoufriou, T. Layout optimization of wireless sensor networks for structural health monitoring. Smart Struct. Syst. 2014, 14, 39–54. [Google Scholar] [CrossRef]
  24. Natarajan, S.; Howells, H.; Deka, D.; Walters, D. Optimization of sensor placement to capture riser VIV response. OMAE 2006, 5, 821–849. [Google Scholar]
  25. Levine-West, M.; Milman, M.; Kissil, A. Mode shape expansion techniques for prediction-experimental evaluation. AIAA J. 1996, 34, 821–829. [Google Scholar] [CrossRef]
  26. Ewins, D.J. Modal Testing, Theory, Practice and Application; Research Studies Press: Baldock, UK, 2000. [Google Scholar]
  27. Soman, R.N.; Onoufriou, T.O.; Votsis, R.; Chrysostomou, C.Z.; Kyriakides, M.A. Multi-type, multi-sensor placement optimization for structural health monitoring of long span bridges. Smart Struct. Syst. 2014, 14, 1–9. [Google Scholar] [CrossRef]
  28. Liu, C.Y.; He, X.Y. Wireless sensor integration for bridge model health monitoring. J. Vibroeng. 2013, 15, 1028–1040. [Google Scholar]
  29. Ye, F.; Zhong, G.; Cheng, J. PEAS: A Robust Energy Conserving Protocol for Long-lived Sensor Networks. In Proceedings of the International Conference on Distributed Computing Systems (ICDCS), Providence, RI, USA, 19–22 May 2003; pp. 28–37.
  30. Liu, C.Y.; Fang, K.; Teng, J. Optimum wireless sensor deployment scheme for structural health monitoring: A simulation study. Smart Mater. Struct. 2015, 24, 115034. [Google Scholar] [CrossRef]
  31. Yao, L.; Sethares, W.A.; Kammer, D.C. Sensor placement for on orbit modal identification via a genetic algorithm. AIAA J. 1993, 31, 1922–1928. [Google Scholar] [CrossRef]
  32. Kennedy, J. Small Worlds and Mega-Minds Effects of Neighborhood Topology on Particle Swarm Performance. In Proceedings of the IEEE Congress on Evolutionary Computation, Las Vegas, NV, USA, 6–9 July 1999; pp. 1931–1938.
  33. Yi, T.H.; Li, H.N.; Zhang, X.D. A modified monkey algorithm for optimal sensor placement in structural health monitoring. Smart Mater. Struct. 2012, 21, 105033. [Google Scholar] [CrossRef]
  34. Guo, Q.G.; Zhou, D.X.; Zhou, X.F. Calculation Method of Optimal Cluster Head in LEACH Routing Protocol. Netw. Commun. 2013, 2, 61–66. [Google Scholar] [CrossRef]
Figure 1. Spherical WSN model.
Figure 1. Spherical WSN model.
Sensors 16 01865 g001
Figure 2. The first-order wireless communication mode.
Figure 2. The first-order wireless communication mode.
Sensors 16 01865 g002
Figure 3. Bridge Elevation View (unit: cm).
Figure 3. Bridge Elevation View (unit: cm).
Sensors 16 01865 g003
Figure 4. Bearing forms constraint of main bridge (unit: cm).
Figure 4. Bearing forms constraint of main bridge (unit: cm).
Sensors 16 01865 g004
Figure 5. The first ten bending mode shapes of the bridge. (a) Lateral mode of the bridge; (b) Vertical mode of the bridge.
Figure 5. The first ten bending mode shapes of the bridge. (a) Lateral mode of the bridge; (b) Vertical mode of the bridge.
Sensors 16 01865 g005
Figure 6. Sensor deployment on the bridge girder.
Figure 6. Sensor deployment on the bridge girder.
Sensors 16 01865 g006
Figure 7. Sensor deployment on bridge cross sections.
Figure 7. Sensor deployment on bridge cross sections.
Sensors 16 01865 g007aSensors 16 01865 g007b
Figure 8. Cluster-tree topology configuration.
Figure 8. Cluster-tree topology configuration.
Sensors 16 01865 g008
Table 1. The energy consumption functions of each unit about different sensor nodes.
Table 1. The energy consumption functions of each unit about different sensor nodes.
Energy ConsumptionWireless Accelerometer NodeWireless Strain Gauge Node
Data acquisitionCca = aaKCcs = asK
Data processingCpa = badamKCps = bsdsmK
Data communicationCta = (ga + dadam)KCts = (gs + dsdsm)K
Data receivingCra = eaKCrs = esK
Table 2. Nomenclature and typical values.
Table 2. Nomenclature and typical values.
ParameterDefinitionValue
αaenergy consumption related to data collecting about accelerometer60 × 10−9 J/bit
βaenergy consumption related to data processing about accelerometer45 × 10−9 J/bit
γaenergy consumption related to the transmission distance about accelerometer45 × 10−9 J/bit
δaenergy consumption related to the transmission distance about accelerometer10 × 10−12 J/bit.m2
εaenergy consumption related to data receiving about accelerometer135 × 10−9 J/bit
daaverage distance between all nodes in a network layer to the base station or cluster head node in upper network layer/
αsenergy consumption related to data collecting about strain sensor45 × 10−9 J/bit
βsenergy consumption related to data processing about strain sensor1.35 × 10−9 J/bit
γsenergy consumption related to data transmitting about strain sensor45 × 10−9 J/bit
δsenergy consumption related to the transmission distance about strain sensor10 × 10−12 J/bit.m2
εsenergy consumption related to data receiving about strain sensor135 × 10−9 J/bit
dsaverage distance between all nodes in a network layer to the base station or cluster head node in upper network layer /
Kamount of data per second transmitting200
Table 3. The first ten natural frequencies of the bridge.
Table 3. The first ten natural frequencies of the bridge.
ModeFrequency (Hz)Modal Shape
10.669lateral
20.746lateral
33.291vertical
43.541vertical
54.051vertical
64.367lateral
75.118lateral
85.370vertical
97.935vertical
108.450vertical
Table 4. Results of the PSO running 10 times.
Table 4. Results of the PSO running 10 times.
Number12345678910Mean
Function Value1541881731381571851229916263144.1
N1
(Accelerometer)
1414151415131513131414
N2
(Strain Gauge)
1010910911911111010
R(m)
(Coverage Range)
103978998991009110110010296.7
Table 5. Network parameters.
Table 5. Network parameters.
ParameterValue
HSWN Space V(0, 0, 0)~(215, −3.12, 13.75)
Accelerometer14
Strain Sensor10
Base Station Location(107.5, 0, 6.875)
γ45 × 10−9 J/bit
ε135 × 10−9 J/bit
δfs10 × 10−12 J/bit/m2
δmp0.0013 × 10−12 J/bit/m2
d100 m
Table 6. Comparison of clustering and non-clustering.
Table 6. Comparison of clustering and non-clustering.
Network TopologyEnergy Consumption (J)Coverage Radius of Network (m)
First LayerSecond LayerEntire Network
Flat0.0280105.4
Cluster-tree0.021032.7102.5135.2
Difference (%)33.328.3

Share and Cite

MDPI and ACS Style

Liu, C.; Jiang, Z.; Wang, F.; Chen, H. Energy-Efficient Heterogeneous Wireless Sensor Deployment with Multiple Objectives for Structural Health Monitoring. Sensors 2016, 16, 1865. https://doi.org/10.3390/s16111865

AMA Style

Liu C, Jiang Z, Wang F, Chen H. Energy-Efficient Heterogeneous Wireless Sensor Deployment with Multiple Objectives for Structural Health Monitoring. Sensors. 2016; 16(11):1865. https://doi.org/10.3390/s16111865

Chicago/Turabian Style

Liu, Chengyin, Zhaoshuo Jiang, Fei Wang, and Hui Chen. 2016. "Energy-Efficient Heterogeneous Wireless Sensor Deployment with Multiple Objectives for Structural Health Monitoring" Sensors 16, no. 11: 1865. https://doi.org/10.3390/s16111865

APA Style

Liu, C., Jiang, Z., Wang, F., & Chen, H. (2016). Energy-Efficient Heterogeneous Wireless Sensor Deployment with Multiple Objectives for Structural Health Monitoring. Sensors, 16(11), 1865. https://doi.org/10.3390/s16111865

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop