Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 218026, 10 pages

http://dx.doi.org/10.1155/2013/218026

Research Article

A Three-Dimensional Model of the Effective Electromechanical Impedance for an Embedded PZT Transducer

Chunyuan Zuo, Xin Feng, and Jing Zhou

Faculty of Infrastructure Engineering, Dalian University of Technology, Dalian 116024, China

Received 15 October 2013; Accepted 4 December 2013

Academic Editor: Ting-Hua Yi

Copyright © 2013 Chunyuan Zuo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A three-dimensional model of the effective electromechanical impedance for an embedded PZT transducer is proposed by considering the interaction between a PZT patch and a host structure. By introducing an effective mechanical impedance, the coupled electromechanical admittance formulations are derived using the piezoelectric constitutive equations. Then, a modified methodology for monitoring structure changes using an electromechanical impedance (EMI) technique is proposed. In the proposed method, the changes in the host structure are monitored by using the “active” part associated with the structural mechanical impedance, which is extracted from the measured raw admittance signatures. The strength gain of a concrete beam with embedded PZT transducers during the curing age was monitored with the proposed methodology. The experimental results demonstrate that the use of the “active” part is more sensitive as opposed to the raw admittance signatures for structural health monitoring (SHM).

1. Introduction

In general, a typical structural health monitoring (SHM) system includes three major components: a sensor system, a data processing system (including data acquisition, transmission, and storage), and a health evaluation system. The sensor system is the key component in the SHM system [1, 2]. The advent of commercially available smart piezoelectric materials has provided new opportunities for SHM. The key characteristic of PZT materials is the utilization of the converse piezoelectric effect to actuate the structure in addition to the direct piezoelectric effect to sense structural deformation [3, 4]. By the direct and converse piezoelectric effects of PZT materials, the electromechanical impedance (EMI) technique simultaneously employs piezoelectric ceramic (PZT) patches as both sensors and actuators.

In recent years, the EMI technique has gained wide attention and emerged as a favorable tool for structural damage detection. This technique can diagnose several types of structures, including bolted connections [5, 6], cable-stayed bridges [7], aerospace structures [8], reinforced concrete structures [9, 10], and pipeline structures [11]. The basic concept of this approach is to use high-frequency excitations to monitor the local area of the host structure for changes in the structural mechanical impedance caused by damage [12]. In the EMI technique, a small-size PZT patch is bonded to or embedded into the monitored structure and excited by an alternating voltage sweep signal. The PZT patch transfers the vibrations to the host structure, which is deformed and produces a local dynamic response. Simultaneously, the response of the host structure is transferred back to the PZT patch. Because of the electromechanical coupling property of PZT materials, the structural mechanical impedance of the host structure can be determined by the electrical impedance of the PZT patch. By monitoring the variations of the measured electrical impedance, the change in the structural mechanical impedance induced by structural damage can be detected. The advantages of EMI-based approaches include but are not limited to a lightweight and inexpensive sensor, high sensitivity to incipient damage for high-frequency actuating, and immunity to ambient or far-field loads.

The interaction between the PZT patch and the host structure is the key issue for sensor design and data interpretation in EMI-based methods. During the last two decades, several modeling approaches have been proposed for describing the PZT-structure interaction. Among them, two well-known approaches are the static approach and the impedance approach. The former one developed by Crawley and de Luis [16] uses a statically determined equivalent force and moment as the amplitude of the forcing function. It is assumed that the mechanical interaction between the bonded actuator and the host structure occurs at the ends of the actuator in the form of concentrated forces. This approach is simple and easy to use but it often leads to misleading results. More importantly, the static approach does not correctly capture the physical essence of the PZT-structure interaction. In order to improve the accuracy, Liang et al. [13] presented the impedance approach instead of the static approach. This approach modeled the PZT patch connected to the host structure represented by the mechanical impedance , as shown in Figure 1. The multidimensional effects associated with PZT vibrations are ignored in Liang’s model. Zhou et al. [14] extended the one-dimensional (1D) model to a generic two-dimensional (2D) PZT-structure model by considering the coupled vibration in the plane. This model replaced the structural impedance with a matrix consisting of the direct impedances and and the cross-impedances and , as shown in Figure 2. Experimental difficulties prevent the direct application of Zhou’s model for extraction of the structural mechanical impedance. To overcome this shortcoming, a new concept of “effective impedance” was introduced in the effective impedance model of Bhalla and Soh [15]. In this model, the PZT-structure mechanical interaction extends all over the finite-sized PZT patch, as shown in Figure 3. This model condenses the two-directional physical coupling between the PZT patch and the host structure into a single impedance term, which can be determined from the measured admittance signatures. Hence, the model bridges the gap between Liang’s 1D model and Zhou’s 2D model.

Figure 1: Liang’s 1D impedance model [13].

Figure 2: Zhou’s 2D impedance model [14].

Figure 3: 2D effective impedance model of Bhalla and Soh [15].

Conventionally, the impedance models are usually considered as 1D or 2D, and the coupling between the extensional and longitudinal vibrations is neglected. To satisfy the requirements of the 1D model, the thickness of the PZT patch must be much smaller than its radius (length and width). However, the PZT patch has finite dimensions of thickness and radius (length and width) in practical applications, and both the extensional and longitudinal vibrations are important [17]. Hence, there is a need for the development of a three-dimensional (3D) PZT-structure interaction model, which considers the coupling between the extensional and longitudinal vibrations. The study presented here aims to model the PZT-structure interaction considering the 3D coupled vibration effect. In this paper, a 3D effective electromechanical impedance model for an embedded PZT plate is proposed by adopting the “effective impedance” concept. The new model describes the extensional actuations along the radial direction and longitudinal actuation along the thickness direction of the PZT patch. The advantage of the 3D effective impedance model is that it can only capture the “active part” associated with the structural mechanical impedance by filtering the PZT contribution. Then, the inherent properties of the monitored structure can be detected with a higher accuracy than with the raw admittance signal only. To validate the effectiveness of the proposed method, a concrete beam with embedded PZT transducers during the curing age was measured to monitor the strength gain.

2. Theoretical Background

In a linear vibration system, the steady-state response , that is, displacement, velocity, or acceleration, is caused by the harmonic force . and have the same frequency, and lags behind by a phase angle . The structural mechanical impedance is defined as the ratio of the harmonic excitation to its steady-state response. Conventionally, the velocity impedance is regarded as the mechanical impedance. Therefore, the structural mechanical impedance is expressed as

where is the amplitude of the excitation force, ω is the circular frequency of the excitation force, is the amplitude of the velocity, and is the phase angle of the velocity response.

The existing models are based on the definition that the structural mechanical impedance refers to a particular point, and the point considered is the PZT end point. However, the PZT-structure interaction is not restricted at the PZT end point, and it extends all over the PZT patch. The force transmission between the PZT patch and the host structure occurs along the entire boundary of the PZT patch [15]. Ignoring the effect of the PZT vibrations in the thickness direction and assuming the frequency range of interest to be much lower than the dominant modes of thickness vibration, a new concept of “effective mechanical impedance” was defined by Bhalla and Soh [15]. However, the embedded PZT patch has both the extensional and longitudinal vibrations in practical applications. Hence, the formulations of the “effective mechanical impedance” of 3D model need to be established.

We consider a circular PZT patch embedded into an unknown host structure, as shown in Figure 4, subjected to an electric field. The interaction of the PZT patch with the host structure is represented in the form of the boundary traction , and the PZT-structure interaction force is assumed along the entire boundary of the PZT patch. Further, the PZT patch is assumed to be extremely small as compared to the host structure, and it has negligible mass and stiffness. Hence, the PZT patch can be regarded as an effective drive point. In addition, the impedance of the effective drive point is defined as the effective mechanical impedance.

Figure 4: 3D effective impedance model.

The “effective mechanical impedance” of the PZT patch is defined as

where is the effective force; is the unit vector normal to the boundary; is defined as effective displacement of the PZT patch, where is the change in volume of the PZT patch; and is the surface area of the PZT patch. Differentiation of the effective displacement with respect to time yields the effective velocity . The “effective mechanical impedance” of the host structure can be defined based on applying the force on the surface of the host structure, along the boundary of the location of the embedded PZT patch.

3. Three-Dimensional Electromechanical Impedance Models

The existing impedance models are usually considered to be 1D or 2D models, and they neglect the longitudinal vibration of the PZT patch. Hence, there is a need for the development of a 3D PZT-structure interaction model that considers both the extensional and longitudinal vibrations of the PZT patch. By applying the concept of “effective mechanical impedance,” a 3D effective impedance model is proposed. The features of the proposed impedance model are that there are no restrictions of the PZT dimensions, and the model more accurately represents the actual situation.

3.1. 3D Electromechanical Admittance Formulation by the Effective Impedance

A circular PZT patch of radius and thickness is shown in Figure 5. The polarization direction is along the axis, and the external exciting electric field is also applied along the -direction and parallel to the polarization direction. Then, the vibration of the PZT patch can be regarded as a symmetrical vibration, and the shearing effect can be ignored. Owing to the symmetrical vibration, it is sufficient to use of the PZT patch and host structure. Because of the mechanical and piezoelectric anisotropy of the PZT patch, the formulations , , , and are obtained. In cylindrical coordinates, the constitutive relations of the PZT patch are expressed as [18]

where , , and are, respectively, the radial, tangential, and axial stress; , , and are, respectively, the radial, tangential, and axial strain; is the elastic compliance constant in the short-circuited condition; is the complex Young’s modulus of elasticity at a constant electric field; , , and are Poisson’s ratios; and are the piezoelectric strain coefficients; is the electric field; is the electric displacement; and is the complex electric permittivity at a constant stress.

(a)

(b)

Figure 5: (a) Geometry of the circular PZT patch and (b) stresses and displacements on the PZT patch.

We let . is defined as the mechanical coupling coefficient between different vibrational modes [17]. By algebraic manipulation, we obtain

where . Then, the effective force of the PZT-structure interaction is given by

In the short-circuited condition, the effective force is expressed as

As derived in [18, 19], the displacements of the PZT patch in cylindrical coordinates are provided by

where the wave number of the radial direction and the wave number of the longitudinal direction . and are constants determined by the boundary conditions. is the Bessel function of the first kind and first order, and . By differentiating (7) with respect to time, the velocities are obtained as follows:

The strains are obtained by differentiating (7) with respect to the coordinate axes:

Under the multimode vibration and free boundary conditions, the frequency equations for the PZT patch in a coupled vibration are obtained as follows:

Then, the following relations are obtained by solving (8) and (9):

where is the root for (9), and it is a function of . By substituting , , and and solving, is obtained as follows:

where .

Based on the definition of the effective displacement in (1), the displacement at the boundaries of of the PZT patch is expressed as

Differentiating with respect to time, the effective velocity is obtained:

To simplify the calculation, we assume that

Substituting (7), the following formulations are derived:

Substituting (6), (8), (9), and (14) into (2), the effective impedance of the PZT patch is expressed as

The effective force is related to the effective impedance of the host structure by

By substituting (8), (9), and (18) into (2), with , the constants and are obtained:

The instantaneous electric current is expressed as

Substituting (4), (9), and (19) into (3), the electromechanical admittance of the 3D effective electromechanical impedance model is derived as

where and ; is the instantaneous voltage across the PZT patch.

This formulation contains and the effective structural mechanical impedance .   is constant for certain material parameters. The variations in caused by a change in the inherent properties of the host structure directly influence the electromechanical impedance. Hence, the variations in the inherent properties can be monitored by the measured electromechanical impedance.

3.2. Decomposition of the Coupled Electromechanical Admittance Signatures

In Section 3.1, the complex electromechanical admittance of the proposed impedance model has been derived. By rearranging the various terms in (21), it can be separated into two parts as

where . Part I only relates to the parameters of the PZT patch and is independent of the host structure. Part II represents the coupling between the PZT patch and the host structure. Hence, (22) is written as

where describes the PZT contribution that is insensitive to any structural changes; it is the “passive” part; represents the contribution of the PZT-structure interaction and is sensitive to any structural changes near the PZT patch; it is the “active” part [20]. By substituting   and , is expressed as

where

The measured electrical admittance signals contain the “passive” and “active” parts. The existence of the “passive” part influences the sensitivity of the raw electrical admittance signals for SHM. Hence, if the PZT parameters are predicted, the “passive” part can be filtered out, and the “active” part can be derived from the raw electrical admittance signals:

Substituting (25) and (26) into (22), the complex form of the “active” part is expressed as

4. Experimental Investigations

In this study, a concrete beam with an embedded PZT patch during the curing age was monitored to measure the strength gain. First, an Agilent4294A precision impedance analyzer was used for collecting the electric admittance signals. By a series of computational procedures, the “active” part associated with the mechanical impedance was derived from the raw electrical admittance signal. Finally, the “active” part was utilized to monitor the strength gain of the concrete beam during the curing age. The objective of this experiment was to demonstrate the feasibility and sensitivity of utilizing the “active” part by filtering out the PZT contribution for SHM.

4.1. Experimental Setup

In the experiment, a concrete beam with dimensions of was used for the strength gain monitoring of concrete in the early age by the proposed approach. To measure the EMI signatures during the curing age, a PZT patch having a thin epoxy adhesive that can be regarded as a waterproof layer was embedded into the concrete structure. The diameter of the PZT patch was 20 mm, and the thickness was approximately 2 mm. The properties of the PZT patch are listed in Table 1. The experimental setup for collecting the electrical admittance signals is shown in Figure 6. The patch was wired to the Agilent4294A precision impedance analyzer. Then, the Agilent4294A precision impedance analyzer was connected to a computer via a GPIB interface.

Table 1: Properties of the PZT patch.

Physical parameter
Value

Density, (kg/m³)
7600

Young’s modulus, (N/m²)

Piezoelectric strain coefficient, (m/V)

Piezoelectric strain coefficient, (m/V)

Figure 6: Experimental setup.

By applying a high-frequency excitation voltage to the PZT patch with the Agilent4294A analyzer, the electromechanical admittance signature produced by the coupling vibration between the PZT patch and the host structure was measured in the frequency range of 20–200 kHz, including the real part (conductance ) and the imaginary part (susceptance ). In the monitored frequency range, 801 data points were taken with a 1-V sweeping sine wave for the excitation. However, the measured signals are often contaminated by noise originating from the measurement system [21, 22]. To minimize the influence of noise, eight averages of the signatures were calculated per frequency point. The impedance signature at the age of 1 day was measured and defined as the baseline, and then subsequent tests were taken at the ages of 2, 3, 4, and 7 days. By analyzing the measured electrical admittance signatures, the strength gain of the concrete beam during the curing age can be monitored by the proposed approach.

4.2. Determination of the PZT Parameters

Owing to the influence of the production process and the high temperature in the welding process between the PZT and the wires, the properties of the PZT patches provided by the manufacturer are unreliable. In the EMI technique, the “active” part of the electrical admittance signals is predicted based on the PZT parameters. Hence, more accurate PZT parameters need to be derived before extracting the “active” part from the measured raw admittance signatures.

For a PZT patch in the free-free condition prior to its embedding into the host structure, the complex electromechanical admittance is derived as ()

By substituting ,   ,  , , , and , we obtain

where

Further, for low frequencies (typically <1/5 of the first resonance frequency of the PZT), and (i.e., , , , and ); therefore, , [13], and the quasistatic sensor approximation is expressed as [23]

Hence, the electrical constants and can be determined by the measured electric admittance signals and for low frequencies (typically <10 kHz). From this analysis, and were determined to be and 0.0096, respectively.

By using (29) with these PZT parameters, the electrical admittance signals of the PZT patch in the free-free condition in the frequency range of 20–220 kHz were obtained and compared with the measured free PZT signals () in Figure 7. Although the experimental and analytical free PZT signals have the good agreement in the figures, there are fundamental differences that warrant closer examination. For instance, the analytical resonance frequency is slightly higher than the experimental frequency.

(a) Conductance

(b) Susceptance

Figure 7: Comparison between the experimental and theoretical signals of the PZT patch in the free condition.

In order to achieve much better agreement between the analytical and experimental signatures, the correction factors and are introduced in the terms and , and these terms may be replaced by and , respectively. By trial and error, the values of the correction factor were and . Hence, very good agreement is observed between the analytical and experimental signatures. Using curve-fitting techniques, it was found that and . From these PZT parameters, the “active” part associated with the structural mechanical impedance by filtering the PZT contribution can be extracted from the measured raw admittance signatures.

4.3. Results and Analysis

As mentioned above, the PZT parameters can be accurately predicted by the measured admittance signatures of the PZT patches in the “free-free” condition. Hence, the “passive” part (PZT contribution) is filtered out, and the “active” part associated with the structural mechanical impedance can be extracted from the measured raw admittance signatures. The spectra are shown in Figure 7. With an increase in the curing time for the concrete specimen, the amplitudes of the and curves gradually diminished, and the curves shifted to the right. This indicates that the “active” part is sensitive to changes in the inherent properties of a concrete beam with embedded PZT transducers during the curing age.

Figure 8 depicts the variations of the “active” part curve, but the extent of the changes in the inherent properties cannot be directly determined. In other words, the spectra of the “active” part only qualitatively represent the strength gain of the concrete beam during the curing age. Although the spectra are qualitatively useful, they lack any quantitative accuracy. Traditionally, statistical pattern recognition techniques have been employed to quantify the changes in the admittance signatures, such as relative deviation, root-mean-square deviation, and mean-absolute-percentage deviation [24]. In order to quantitatively analyze the strength gain for the concrete beam during the curing age, two statistical indices are defined as follows.

(a) Conductance

(b) Susceptance

Figure 8: “Active” part extracted from the measured raw admittance signatures.

Root-mean-square deviation (RMSD)

Correlation coefficient (CC)

where and separately denote the vector and its transpose for the “active” part of the electrical admittance ( and ) for curing times of 2, 3, 4, 5, and 7 days; and and separately denote the vector and its transpose for the “active” part of the electrical admittance for a curing time of 1 day.

Figure 9 shows the plots of the RMSD and CC indices for the “active” part and the measured raw admittance signatures, which are calculated at various stages during the curing age with respect to the baseline signature. From Figure 8, we find that the RMSD index increased and the CC index decreased during the curing age for the concrete beam. These indexes qualitatively represent the variations in the inherent properties of the interior of the concrete beam.

(a) and

(b) and

Figure 9: Comparison between the “active” part and the measured raw admittance signatures.

The existing EMI techniques only depend on the real part of the measured electrical signatures, but the information in the imaginary part is lost. By eliminating the “passive” part, the admittance signatures are more sensitive to any structural changes than the raw signatures. By comparing the “active” part and the measured raw admittance signatures, the indices of the real part are similar, whereas the indices of the imaginary part have a large difference. Further, the indices of the imaginary component are similar to the indices of the real component in the “active” admittance signatures. This indicates that utilization of the imaginary component for damage identification is feasible. The experimental results demonstrate that the “active” part of the admittance signatures is more sensitive than the raw admittance signatures for SHM. Furthermore, the imaginary component can also be used for SHM.

5. Conclusions

In this paper, a three-dimensional effective electromechanical impedance model for an embedded PZT transducer is proposed. Based on the concept of effective impedance, the 3D coupled electromechanical admittance equation is derived. By rearranging the various terms, the admittance equation is separated into an “active” part and a “passive” part. Because the “active” part is sensitive to any structural changes near the PZT patch, it is feasible to use the “active” part for SHM. Then, the proposed method is used to monitor the strength gain of a concrete beam with embedded PZT transducers during the curing age. The results show that the “active” part of the admittance signatures is more sensitive to any structural damage than the raw admittance signatures, which demonstrates the validity of the proposed method.

Acknowledgments

This work was supported by the Sate Key Development Program for Basic Research of China with Grant no. 2013CB035906 and Grants from the National Natural Science Foundation of China (51121005, 51078060). The authors gratefully acknowledge the support from the Education Department of Liaoning Province with Grant no. L2012010.

References

T. Yi, H. 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, 12 pages, 2011.

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

I. Chopra, “Review of state of art of smart structures and integrated systems,” AIAA Journal, vol. 20, no. 11, pp. 2145–2187, 2002.

S. Yan, N. Z. Zhao, Q. J. Wang, and G. B. Song, “Dynamic mechanical model of surface-bonded PZT actuator: theory and experiment,” Applied Mechanics and Materials, vol. 303–306, pp. 1732–1735, 2013.

K. Nguyen, J. Park, and J. Kim, “Imote2-based multi-channel wireless impedance sensor nodes for local SHM of structural connections,” in Sensors and Smart Structures Technologies for Civil, Mechanical, and Aerospace Systems, vol. 7981 of Proceedings of SPIE, pp. 79811P–79811P-11, San Diego, Calif, USA, March 2011.

J. T. Kim, K. D. Nguyen, and J. H. Park, “Wireless impedance sensor node and interface washer for damage monitoring in structural connections,” Advances in Structural Engineering, vol. 15, no. 6, pp. 871–885, 2012.

K. D. Nguyen, P. Y. Lee, and J. T. Kim, “Smart PZT-interface for SHM in tendon-anchorage of cable-stayed bridge,” in Sensors and Smart Structures Technologies for Civil, Mechanical and Aerospace System, vol. 8345 of Proceedings of SPIE, pp. 834519–834519-11, San Diego, Calif, USA, March 2012.

Z. A. Chaudhry, T. Joseph, F. P. Sun, and C. A. Rogers, “Local-area health monitoring of aircraft via piezoelectric actuator/sensor patches,” in Smart Structures and Materials 1995: Smart Structures and Integrated Systems, vol. 2443 of Proceedings of SPIE, pp. 268–276, San Diego, Calif, USA, March 1995.

J. W. Ayres, F. Lalande, Z. Chaudhry, and C. A. Rogers, “Qualitative impedance-based health monitoring of civil infrastructures,” Smart Materials and Structures, vol. 7, no. 5, pp. 599–605, 1998.

G. Song, H. Gu, Y. L. Mo, T. T. C. Hsu, and H. Dhonde, “Concrete structural health monitoring using embedded piezoceramic transducers,” Smart Materials and Structures, vol. 16, no. 4, pp. 959–968, 2007.

G. Park, H. H. Cudney, and D. J. Inman, “Feasibility of using impedance-based damage assessment for pipeline structures,” Earthquake Engineering and Structural Dynamics, vol. 30, no. 10, pp. 1463–1474, 2001.

G. Park, H. Sohn, C. R. Farrar, and D. J. Inman, “Overview of piezoelectric impedance-based health monitoring and path forward,” Shock and Vibration Digest, vol. 35, no. 6, pp. 451–463, 2003.

C. Liang, F. P. Sun, and C. A. Rogers, “Coupled electro-mechanical analysis of adaptive material systems—determination of the actuator power consumption and system energy transfer,” Journal of Intelligent Material Systems and Structures, vol. 5, no. 1, pp. 12–20, 1994.

S. Zhou, C. Liang, and C. A. Rogers, “An impedance-based system modeling approach for induced strain actuator-driven structures,” Journal of Vibration and Acoustics, vol. 118, no. 3, pp. 323–331, 1996.

S. Bhalla and C. K. Soh, “Structural health monitoring by piezo-impedance transducers. I: modeling,” Journal of Aerospace Engineering, vol. 17, no. 4, pp. 154–165, 2004.

E. F. Crawley and J. de Luis, “Use of piezoelectric actuators as elements of intelligent structures,” AIAA Journal, vol. 25, no. 10, pp. 1373–1385, 1987.

S. Y. Lin, “Coupled vibration analysis of piezoelectric ceramic disk resonators,” Journal of Sound and Vibration, vol. 218, no. 2, pp. 205–217, 1998.

IEEE, “IEEE standard on piezoelectricity,” IEEE 176-1987, 1988.

A. H. Meitzler, H. M. O'Bryan, and H. F. Tiersten, “Definition and measurement of radial mode coupling factors in piezoelectric ceramic materials with large variations in Poisson's ratio,” IEEE Transactions on Sonics and Ultrasonics, vol. 20, no. 3, pp. 233–239, 1973.

S. Bhalla and C. K. Soh, “Structural health monitoring by piezo-impedance transducers. II: applications,” Journal of Aerospace Engineering, vol. 17, no. 4, pp. 166–175, 2004.

T. H. Yi, H. N. Li, and X. Y. Zhao, “Noise smoothing for structural vibration test signals using an improved wavelet thresholding technique,” Sensors, vol. 12, no. 8, pp. 11205–11220, 2012.

T. H. Yi, H. N. Li, and H. M. Sun, “Multi-stage structural damage diagnosis method based on “energy-damage” theory,” Smart Structures and Systems, vol. 12, no. 3-4, pp. 345–361, 2013.

V. Giurgiutiu and A. N. Zagrai, “Embedded self-sensing piezoelectric active sensors for on-line structural identification,” Journal of Vibration and Acoustics, vol. 124, no. 1, pp. 116–125, 2001.

K. K.-H. Tseng and A. S. K. Naidu, “Non-parametric damage detection and characterization using smart piezoceramic material,” Smart Materials and Structures, vol. 11, no. 3, pp. 317–329, 2002.