1. Introduction
The traditional nondestructive testing (NDT) methods are useful, but they are time consuming and expensive as they require the structures to be off-service and disassembled in components. Alternatively, on-line structural health monitoring (SHM) provides a continuous assessment of structures and predicts the remaining lifespan using permanently embedded sensing systems [
1]. The SHM techniques based on Lamb waves (a type of guided wave) have been extensively used for thin plate-like structures because they can travel a long distance with low attenuation, and their propagation is sensitive to the presence of small damage. The piezoelectric (lead zirconate titanate, PZT) transducers can both excite and sense the Lamb wave signals and can be permanently attached to the host structure in the form of a network of sensors for fast multipoint measurement in situ SHM [
2]. These PZT transducers are small, lightweight, consume low power, and produce response in a wide frequency region [
3]. Furthermore, PZT transducers can be used in harsh temperature and radiation environments for SHM of various types of damage in metallic and composite plates [
4].
There are several damage estimation methods based on the sparse array of PZT transducers, which have shown a great potential for the identification of local damage. However, the Lamb wave signals generated by PZT transducers are dispersive and multimodal. An accurate interpretation of transducer data is important for the extraction of damage features such as location and severity from complicated Lamb wave signals. Therefore, sophisticated signal processing techniques such as fast Fourier transform [
5], short-time Fourier transform [
6], Hilbert transform [
7], warped frequency transform [
8], and signal-difference-based correlation coefficient [
9,
10] have been used to accurately analyze the dynamic response signals scattered by damage. Among these methods, the continuous wavelet transformation (CWT) [
11,
12] has recently gained attention and popularity due to its powerful time–frequency feature extraction from complicated signals, its ability to detect the signal singularity for physical damage inspection, and to depress the dispersive characteristics of the Lamb waves [
13,
14].
In recent years, Lamb-wave-based imaging methods such as phased array [
15,
16] and delay-and-sum [
17,
18,
19] algorithms have been widely used for damage localization in plate-like structures. The damage imaging methods have been further empowered by probabilistic approaches such as Probability-Based Diagnostic Imaging (PDI) [
20,
21,
22], Reconstruction Algorithm for Probabilistic Inspection of Damage (RAPID) [
23,
24], Algebraic Reconstruction Techniques (ART), and Lamb Wave Tomography (LWT) [
25,
26]. These methods have used different parameters to define the probability regions for individual sensing paths, which were then used to reconstruct the damage images. For example, the PDI methods have used the weight distribution functions [
20,
21,
27,
28,
29,
30] to define the probability region between actuator–sensor pairs. These distribution functions depended on a scale/shape factor, β, which controlled the effective size of the elliptical probability area for each sensing path. Similarly, the RAPID and LWT methods used the difference between damaged and undamaged signals to define the signal difference coefficient (SDC) [
23,
24,
31,
32,
33,
34]. The SDC also depended on β, which controlled the effective size of the probability region. In most of these methods, the probability of finding the damage depended on the intersection of direct paths between actuator–sensors pairs [
20,
21,
23,
24,
25,
30,
33,
34,
35,
36]. The probability of accurately locating the damage increases if more direct path intersections fall around the damage region. Consequently, the use of more transducers increases the path intersections and obtains more accurate damage imaging. That is one of the reasons why these methods faced difficulties locating the multiple damage with fewer transducers. Furthermore, these methods were based on the signal difference obtained by the comparison of the damaged signal with the baseline signal, which required additional data to be recorded and processed [
27,
29,
30,
31,
33,
34,
35,
36], and the baseline data may not always be available for practical applications.
In the current research, instead of using any probability distribution function, the possibilities of damage presence are directly defined by CWT of the normalized damage-scattered signals. No separate parameter is used to define the size of the possibility region, it is directly defined by the damage-scattered peak in the normalized CWT signals. The damage possibilities from individual sensing pairs are fused using conjunctive and compromised fusion schemes to image the damage. The method was tested for pitch–catch and pulse–echo arrangements of transducers. The damage localization was improved by a time compensation in all the original output signals. First, the single and multiple damage cases were tested using the difference signals, and then the results were obtained for the same cases directly from the damaged signals, without the baseline.
The organization of the article is as follows:
Section 2 introduces the background and basic theory of the proposed method,
Section 3 shows the implementation of the proposed method using difference signals, the baseline-free implementation is shown in
Section 4, the influence of the environmental noise is presented in
Section 5, experimental validation of the proposed method is presented in
Section 6, a comprehensive discussion of results is presented in
Section 7, and concluding remarks are summarized in
Section 8.
2. Damage Estimation Method Based on CWT
The proposed damage estimation method describes a color-mapped image for each actuator–sensor pair, which indicates the possibility of the presence of damage in the discretized monitoring region. The region of maximum possibilities in the aggregate image of all the pairs represents the damage features. In time-of-flight (ToF)-based triangulation of a point damage
, when actuator
excites the Lamb waves, sensor
receives two signals: a direct signal from
to
, which travels a distance,
, and a damage-scattered signal, which travels a longer distance,
, as shown in
Figure 1.
The ToF from
to
for the damaged-scattered signal,
, can be written as,
where,
is group velocity of the incident Lamb wave and,
The theoretical solution of Equation (1) represents an elliptical path locus with focal points at
and
, indicating the possible damage locations on it, as shown in
Figure 1. When transducers are arranged as an active sensor network, the actual damage location falls at the intersection of multiple elliptical path loci from different actuator–sensor pairs. For damage localization using Equation (1), the accurate extraction of ToF is important therefore, the CWT is introduced to process the signals due to its high time–frequency resolution.
2.1. Normalized Coefficient of CWT
The coefficient of CWT for square-integrable function
in the time domain can be written as,
where
and
are scale and time-shift parameters of wavelet function
, respectively. During the transform, a dynamic signal is represented using these dual parameters. In this research, the scale is set as the ratio between sampling frequency and central frequency as
where
is the sampling frequency calculated using the sampling period as
, and
is the central frequency of the excitation signal. The bar accent represents complex conjugate of the wavelet function
. The selection of the mother wavelet function
is crucial for CWT applications. The Gabor wavelet function was chosen in this research for its higher time resolution to capture the ToF accurately. It has an accurate time–frequency feature extraction capability in ToF-based damage localization problems [
11]. The Gabor wavelet function can be expressed as
where
and
are the real constants.
, so that the value
becomes equal to frequency
, and
was set to satisfy the following admissibility condition:
For the monitoring region discretized into a number of grid points, the possibility of damage existence at each grid point is directly defined by the coefficient of CWT for each actuator–sensor pair. The application of the coefficient of CWT, based on the Gabor wavelet, on damage-scattered signals constructs a CWT envelope for each pair. The possibilities of damage existence are expressed within the range from 0 to 1 with the maximum possibility at a grid point with value 1. Therefore, all damage-scattered signals,
were normalized by dividing with the maximum absolute value of amplitude of the corresponding signal as,
The Gabor wavelet coefficient modulus of the normalized damaged-scattered signal then extracts the damage information for further processing of damage image. It is performed as
2.2. Time Compensation
The value of
in Equation (1) represents the total ToF of excitation signal from actuator to damage and then from damage to sensor. This is the time in the dynamic response signal when the sensor starts receiving the damage-scattered signal [
21]. Consider a normalized damage-scattered signal presented in
Figure 2 as an example, the traveling time from actuator to sensor,
should be from the beginning to the point when damage-scattered signal starts to become non-zero as shown in
Figure 2a. However, this non-zero point is very difficult to identify in practice. The envelope of the coefficient of CWT of this signal, based on Equation (7), is presented with red line in
Figure 2b. The excitation signal used throughout this research was a 5.5-cycle sinusoidal tone burst modulated by Hann window at central frequency 383 kHz, as shown in
Figure 3. For the central frequency of 383 kHz and sampling time of 0.02
, the scale parameter
was calculated as 130.5 using Equation (3). If
is the total acting time of the excitation signal as shown in
Figure 3a, a good time–frequency analysis feature of the CWT can depress the dispersion of Lamb wave effectively to make the signal contain the same duration of time during wave propagation.
The maximum value 1 in the CWT envelope in
Figure 2b appears at the center of the normalized damage-scattered peak, which is half of the time
ahead of the time
when sensor actually received the non-zero point of the damage-scattered signal. Therefore, a time lag was introduced in all the CWT signals by subtracting
as time-compensation
, shown with blue line in
Figure 2b. If
is the total number of cycles and
is the central frequency of the excitation signal,
can be expressed as
2.3. Imaging Algorithm
The monitoring region containing the actuator–sensor pairs is divided into a uniform mesh of
grid points. The possibility of damage existence at a grid point
for an actuator–sensor pair is represented by the damage index for a sensing pair (
) where,
,
and
is the total number of actuator–sensor paths. The solution of Equation (1) on the mesh grid for any actuator–sensor pair provides an image with an infinite number of path loci. The grid points that have the same value of
fall on one path locus. The time-compensated CWT envelope of the normalized damage-scattered signal
defined by Equation (7) provides
) at different time values,
. The value of
can be expressed as
The image defined by provides possible damage locations with the maximum possibility of damage existence at the time when the peak value 1 appears in the CWT envelope. An imaging algorithm was designed in MATLAB to implement Equation (9) in three steps: In the first step, the envelopes provided by the time-compensated coefficient of CWT of the normalized damage-scattered signals were plotted based on Equation (7). In the second step, the infinite number of elliptical path loci provided by Equation (1) were plotted on the discretized domain, which provided a color-mapped image for each sensing pair. Since, Equation (1) represents the ToF of the damage-scattered signals, in step 3, the field values plotted in step 2 were replaced by the values from the CWT envelopes at the corresponding times. The image obtained in step 3 represented possible damage locations in the discretized domain for individual sensing pairs, with the maximum possibility corresponding to the same time and location as in CWT envelopes.
A damage index,
then defines the aggregate image by fusion of
images for each actuator–sensor pair. Two different image fusion schemes were tested for the proposed method: conjunctive fusion scheme based on the multiplication of
values from all pairs at a grid point and compromised fusion scheme based on addition of the values. The damage indices for compromised fusion
and conjunctive fusion
schemes can be written as
4. Implementation Using the Damaged Signals
The baseline signals from the undamaged plate may not always be available to calculate the difference signals in the practical environment. It also increases the computational cost by calculating two additional signals for each sensing path, i.e., undamaged and difference signals. In addition, a tiny incorrect baseline can make the whole estimation wrong. Therefore, the proposed method was implemented using signals directly from the current state of the plate. The baseline-free implementation presented here is for the multiple damage scenario shown in
Figure 16. The damaged signal and coefficient of CWT for sensing path A8B8 are shown in
Figure 20a,b, respectively. The damage-scattered fluctuations in these signals are too small compared to the excitation fluctuation; therefore, the region of these plots containing the damage-scattered fluctuations is shown separately in
Figure 20c,d, respectively. The damage-scattered peaks can be clearly seen in
Figure 20d.
The damage possibilities for sensing path A8B8 using baseline-free estimations are presented in
Figure 21 along with actual D1 and D2.
Figure 21a contains the possibilities without time compensation and
Figure 21b with time compensation. Since, the damage-scattered peaks are too small, the damage possibility values for D1 is hardly visible, whereas for D2 it is so small that it is not even visible. It is because the possibility values for the excitation peak is very high. Similar to
Figure 20c,d, the damage-focused regions of
Figure 21a,b are shown separately in
Figure 22a,b, respectively. The highest damage possibility region near the actuator–sensor locations has been left out of these images, and the distributions around D1 and D2 are more visible now compared to surrounding regions.
Similarly, the damage possibilities from the other seven sensing paths were plotted and fused for damage imaging. The high-possibility regions for each sensing pair were near the actuator–sensor pairs, just as for the case of A8B8. The aggregate images without and with time compensation are shown in
Figure 23a,b, respectively. The time-compensated signals could accurately image both the damage, whereas the localization of damage without compensation did not work for this case either. It should be noticed that the high-possibility regions near each actuator–sensor pair disappeared in the aggregate image. It is because the possibility values from other sensing pairs at these locations were almost 0, which reduced the aggregate value to 0 using multiplication-based
. The individual damage possibility values from all 8 sensing pairs and the aggregate value at damage locations based on
are presented in
Table 1 for the pulse–echo configuration. Non-zero damage possibility values from all sensing pairs existed at the damage locations, which made the non-zero aggregate value appear only at the damage locations.
5. Influence of Environmental Noise
It is essential to test the noise tolerance of the proposed method for the real-world signals, which may be noisy and can influence the results when damage-scattered peaks are too small. The white noise signals in Gaussian distribution were artificially added to the original damaged signals. A high value of noise, i.e., signal-to-noise ratio (SNR) of 30 dB, was added to all the damaged signals, and envelopes formed by the coefficient of CWT were reevaluated. The signals shown in
Figure 20a,b are presented again with SNR 30 dB in
Figure 24a,b, respectively. A high influence of the noise in the signals can be observed in damaged-focused region of these plots separately shown in
Figure 24c,d, respectively.
The possibilities for damage existence for all sensing paths in the pulse–echo configuration were evaluated with noisy signals and then fused to get the aggregate damage image. The results for multiple damage scenario with noisy baseline-free signals are shown in
Figure 25a,b, respectively, without and with the time compensation. Although the noise has a little influence on the resolution of damage imaging, the proposed method could still accurately locate both the damage, and also estimate the severity of one damage with respect to the other damage.
6. Experimental Evaluation
The accuracy of the proposed method was further tested through experimental investigation performed on an 800 × 800 × 1.5 mm
3 aluminum plate with PZT wafers arranged in the form of two square cells of 300 × 300 mm
2 and 240 × 240 mm
2, as shown in
Figure 26a. The minimum size of the damage that can be estimated using the Lamb wave should be equal to or greater than half of its wavelength. For 383 kHz frequency excitation and
m/s group velocity, the wavelength was about 14 mm, therefore, a damage was formed as a punch hole of 7 mm diameter located at (−90, 30) mm. The size of each PZT wafer used in experiments was 5.4 mm diameter, and they were made up of PSN-33 with a density of 7.70 × 10
3 kg/m
3. The PZT wafers had a resonance frequency of 383 kHz, and they could excite/sense the Lamb wave signals in radially outward directions. The real-time laboratory setup of the experiment is shown in
Figure 26b. The signal generator (RIGOL, DG1022) was programmed to generate the excitation signal shown in
Figure 3a. The signal was amplified to 4.2 V using a power amplifier (KROHN-HITE, 7602M). The dynamic response signals were recorded at each receiving transducer using a four-channel oscilloscope (LeCroy, LC574AL). The response signals were captured over a sampling time of 200
at the sampling rate of 1 GHz. In order to reduce the noise, each signal was acquired for an average of 500 times.
When A1 excited the Lamb wave signal, the damaged signal recorded at B1 is presented in
Figure 27a, which contains three fluctuations: excitation, damage-scattered, and edge-scattered. Similarly, when A2 was excited, the recorded signal at B2 is presented in
Figure 27b, which contains two fluctuations: one is excitation fluctuation and another is damage-scattered overlapped with edge-scattered. The envelopes of the coefficient of CWT for signals at A1B1 and A2B2 are presented in
Figure 27c,d, respectively. The time compensation was already made in the experimental signals during the data recording as the oscilloscope was set to trigger or start recording when peak of the voltage signal appears. The damage possibilities for sensing paths A1B1 and A2B2 are presented in
Figure 28a,b, respectively, along with actual size and location of damage D3. The damage-scattered fluctuation in
Figure 28b is overlapped with edge-scattered fluctuation, therefore it does not represent the correct possibilities of damage existence. In the current experimental evaluation, the baseline signals from the undamaged plate were not utilized, therefore the damage possibilities from sensing pair A2B2 were not used in the aggregate damage imaging. Similarly, the damage possibilities for all useful sensing paths in the pulse–echo configuration were evaluated and fused to get the aggregate damage image, presented in
Figure 29.
The results validate that the proposed method can locate the damage with relatively good accuracy also through experimental data. The method proposes the arrangement of PZT transducers as a network of square detection cells on a large structure. Since, the wave velocity and location of the plate edges are already known, the total time for wave propagation can be calculated accordingly to avoid the edge scatterings for all the central cells. For detection cells near the physical boundary, there may be both the damage-scattered and edge-scattered fluctuations, which can be identified separately as shown in
Figure 27a for sensing pair A1B1.
In a rare case, when distances between damage location and receiver and plate’s edge and receiver are same for one of the sensing pairs, the damage-scattered and edge-scattered fluctuation may overlap, as shown in
Figure 27b for sensing pair A2B2. The possibilities of damage existence from such sensing paths where the damage scatterings are weak and overlapped with incident waves or edge scatterings can be ignored in the aggregate damage imaging. Alternatively, the baseline can be used for cells near the physical boundary of the plate in practical applications.
7. Discussion of Results
Two different excitation configurations, pitch–catch and pulse–echo, and damage index based on two different image fusion schemes, compromised fusion
and conjunctive fusion
, were tested in the proposed method. The results presented in
Figure 10 and
Figure 11 verified that the pulse–echo configuration with the conjunctive fusion scheme was most suitable for the proposed method. The pulse–echo configuration was also computationally inexpensive as it provided better accuracy with only 8 excitations compared to the pitch–catch configuration with 28 excitations.
The comparison of signals presented for pulse–echo configuration in
Figure 8,
Figure 13,
Figure 17, and
Figure 20 with pitch–catch configuration signals presented in
Figure 7 showed that the damage scatterings were much stronger than the edge scatterings in case of pulse–echo configuration, whereas the edge scatterings were stronger in case of pitch–catch configuration, for the same size of damage in both cases. This made the high-possibilities region appear at the edge location instead of the damage location during pitch–catch estimation, as shown in
Figure 9a. This is one of the reasons why inaccuracies appeared in damage estimation based on the pitch–catch configuration.
The damage possibilities envelopes defined by the coefficient of CWT were normalized in order to range the possibility values from minimum 0 to maximum 1. The damage index
produced better estimation of the damage region because it is based on the multiplication of the possibility values at grid points from each sensing path. When the damage possibilities for an individual path were too high at a non-damaged location, it still reduced to zero in the aggregate image when multiplied with the values from other sensing paths at the same region. The non-zero values of possibility in the aggregate image appeared only at the damage location because the values from all sensing paths were non-zero at that location. For the same reason,
worked well even when the damage-scattered fluctuations were too small compared to the excitation fluctuation S
0 for baseline-free signals in
Figure 20. Although, the damage-scattered peaks and damage possibilities were too small for an individual sensing path, all the higher-value regions disappeared in the aggregate image because values from other paths at those locations were 0, leaving behind only the damage region having a non-zero contribution from all sensing paths, as evident from
Figure 23b.
The results presented in
Figure 19b and
Figure 23b indicated that multiple damage could be accurately estimated with and without the baseline signals. The method could not only accurately locate the damage but also provided the severity of each damage with respect to the other. In these imaging results, the damage D2 of 12 mm diameter appeared to have a higher possibility values than the smaller damage D1 of 10 mm diameter. It made the damage image for D2 more severe in intensity and extent than that of D1.
The damage index was based on CWT signals, which have a better signal-to-noise ratio and can extract the time–frequency characteristics exactly from the noisy real-world signals. That is why CWT could process the Lamb signals directly from the damaged plate without a baseline signal from an undamaged plate. In the proposed pulse–echo configuration, the damage-scattered peaks in CWT envelopes were much more dominant than other peaks, therefore any additional noise in real-world problems would still ensure that the damage-scattered peaks have the maximum value of 1 in the CWT envelopes for accurate damage estimation.
Generally, the Lamb waves are dispersive in nature, which causes an error in determining its characteristics. The time compensation equal to the half period of the excitation signal is valid when the Lamb waves are not dispersive. In case of a large dispersion, the wave signal will have more packets than the excitation after a certain propagation distance, and the time compensation will be larger than the half period of the excitation signal. CWT is also useful as it can suppress the dispersive nature of the Lamb waves and can also accurately process noisy signals. Furthermore, a lower frequency of the excitation signal was chosen in this method, and in the lower-frequency range of dispersion curves shown in
Figure 5, the Lamb waves were less dispersive and the S
0 mode was well separated from the A
0 mode.
8. Conclusions
A damage estimation method based on coefficient of CWT of the normalized damage-scattered signals is proposed in this research. The possibilities of damage existence are defined by the envelopes of coefficient of CWT without the need of any additional parameter, which makes it simple and easy to implement. In order to avoid the mismatch between actual damage location and maximum possibility peak in CWT envelopes, a time compensation equal to the half period of the original excitation signal is required in all the CWT signals. The proposed method can be used to estimate the multiple damage scenario with difference signals and also with damaged signals without the baseline. The pulse–echo excitation configuration with the damage index based on the conjunctive fusion scheme provides the best estimation with fewer excitations. Due to being baseline-free and accurate with fewer excitations compared to existing methods, the propose method is computationally inexpensive for on-line monitoring, which can reduce the evaluation cost and time during in situ SHM. Many large modern structures have PZT transducers already embedded in them for a continuous SHM during their service life. In the proposed method, the PZT transducers can be arranged as a network of square detection cells on a large structure, and multiple damage inside each of the cells can be estimated separately. Each of the actuator–sensor pairs should be in a close vicinity, and the excitation frequency should in the lower range of the dispersion curves. There may also be additional fluctuations in the dynamic response signals caused by the geometric discontinuities. However, the location of the geometric discontinuities, such as a stiffener or plate edges, are already known, therefore it is easy to identify the fluctuations caused by the scattering from those discontinuities in the time domain. The proposed method can be applied for the estimation of corrosion damage in large metallic plate-like structures.