Next Article in Journal
Silicon Photonic Biosensors Using Label-Free Detection
Next Article in Special Issue
A 40-MHz Ultrasound Transducer with an Angled Aperture for Guiding Percutaneous Revascularization of Chronic Total Occlusion: A Feasibility Study
Previous Article in Journal
A Fuzzy Tuned and Second Estimator of the Optimal Quaternion Complementary Filter for Human Motion Measurement with Inertial and Magnetic Sensors
Previous Article in Special Issue
Fabrication and Characterization of High-Frequency Ultrasound Transducers Based on Lead-Free BNT-BT Tape-Casting Thick Film
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Approach to Guided Wave Ray Tomography for Temperature-Robust Damage Detection Using Piezoelectric Sensors

1
Research Center of Smart Networks and Systems, The Department of Electronic Engineering, Fudan University, Shanghai 200433, China
2
Institute of Acoustics, School of Physics Science and Engineering, Tongji University, Shanghai 200092, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(10), 3518; https://doi.org/10.3390/s18103518
Submission received: 5 September 2018 / Revised: 2 October 2018 / Accepted: 17 October 2018 / Published: 18 October 2018
(This article belongs to the Special Issue Ultrasound Transducers)

Abstract

:
In this paper, a new approach to guided wave ray tomography for temperature-robust damage detection with time-of-flight (TOF) temperature compensation is developed. Based on the linear relationship between the TOF of a guided wave and temperature, analyses show that the TOF of the baseline signal can be compensated by the temperature measurement of the inspected materials without estimating the temperature compensation parameters. The inversion is based on the optimization of the TOF misfit function between the inspected and compensated baseline TOFs of the guided waves, and is applied by the elastic net penalty approach to perform thickness change mapping in a structural health monitoring (SHM) application. Experiments that are conducted in isotropic plates by piezoelectric sensors demonstrate the effectiveness of the proposed method. According to the results, our approach not only eliminates the artefacts that are caused by a temperature variation from 25 °C to 70 °C but also provides more accurate and clearer imaging of damage than conventional ray tomography methods.

1. Introduction

Damage of industrial infrastructure, aircraft, bridges, pipes, and rails is a significant problem that involves various safety issues [1]. Ultrasonic guided waves, which are used as a reliable tool in structural health monitoring (SHM) for detecting and monitoring damage, have received considerable attention because these waves can travel over long distances with less signal loss and thus can cover large monitoring areas [2,3,4,5]. In an SHM system that uses guided waves for damage detection, ultrasonic transducers can be permanently attached to plate-like or pipe-like structures. Each transducer, in turn, generates guided waves, which propagate through the structure and interact with damage, and the remaining transducers record the receiving wave motion. These recorded signals can be processed to generate an image with the locations of damage via many algorithms such as travel-time tomography [6,7,8,9], diffraction tomography [10,11], the hybrid algorithm for robust breast ultrasound tomography (HARBUT) algorithms [12,13] and the full waveform inversion (FWI) method [14,15,16,17,18]. Damage in inspected materials can cause changes in the guided wave waveforms. However, temperature variations of materials can also greatly affect the amplitudes and phases of guided waves, especially those of baseline subtraction signals [19,20]. In other words, changes in temperature can result in changes in the guided wave waveforms. The tomography algorithms that are reported in [6,7,8,9,10,11,12,13,14,15,16,17], which image damage mainly based on the differences between the inspected and baseline signals, cannot distinguish whether the changes in the waveforms are caused by temperature variations or damage in the detected material.
To mitigate the effect of temperature on guided waves, a considerable amount of research has been undertaken on temperature compensation of guided waves in the last decades. The baseline signal stretch (BSS) method [21,22], optimal baseline selection (OBS) method [23,24], and hybrid algorithms that combine the BSS and OBS methods are the most commonly used approaches [25,26,27]. The BSS method compensates for the effect of temperature by modifying the waveform of the baseline signal to make it match the inspected signal in terms of the mean square deviation. However, it is only suitable for a temperature variation range of a few degrees [20,21]. To accommodate a larger temperature variation range, the OBS method uses criteria that are similar to those that are used in the BBS procedure to select a best-matched waveform from the baseline signal database, which contains signals that were measured at various temperatures. This method requires many baseline signals that have sufficiently low post-subtraction noise levels for various temperatures. In various situations, it may be impossible to acquire a baseline signal database with the required temperature resolution. The hybrid algorithms, which combine BBS and OBS, can achieve temperature compensation with a reduced number of baseline signals. This method selects a best-matched waveform from a baseline signal database that has low temperature resolution. The best-matched waveform is then used to compensate the inspected signal via the BBS method. In addition to these algorithms, other temperature compensation methods have been proposed in the literature, such as combinations of the OBS method and the adaptive filter algorithm [28] and of the Hilbert transform and the orthogonal matching pursuit algorithm [20]. As discussed above, most temperature compensation methods model the differences between the baseline and inspected signal waveforms that are caused by the temperature variations as an amplitude scale factor and a phase shift. Once the amplitude scale factor and the phase shift have been estimated, the guided wave waveform can be obtained via temperature compensation. These methods are limited by the estimation accuracies of the amplitude scale factor and the phase shift, which depend on the length of the time window that is used in the temperature compensation parameter estimation. In addition to altering the waveforms of guided waves, temperature can affect the piezoelectric transducer [29]. However, the effect of temperature on transducer performance is significantly less than the effect of temperature on guided wave propagation [30].
In this paper, we propose a new approach to guided wave ray tomography for temperature-robust damage detection. Using the linear relationship between the TOFs of guided waves and the temperature variations, it is shown that the TOFs of the baseline signal can be compensated by the temperature measurement of the inspected material without estimating the temperature compensation parameters. After that, an optimization problem that is based on the difference between the inspected and compensated baseline TOFs of the guided waves is formulated via an elastic net penalty approach, and is resolved to image the damage in an SHM application. To evaluate the effectiveness of our approach, experiments were conducted in isotropic plates using piezoelectric sensors. The experimental results demonstrate that our approach can eliminate the artefacts that are caused by a temperature variation from 25 °C to 70 °C. They also show that more accurate and clearer imaging of damage than that from conventional ray tomography methods can be achieved.
The remainder of this paper is organized as follows: in Section 2, a temperature compensation method that is based on the TOF of guided waves is presented, and a TOF difference model for damage monitoring based on the differences between the inspected and baseline TOFs of the guided waves is proposed. Using our model, an elastic net penalty approach is used to perform damage imaging via ray tomography. In Section 3, the experimental setup for monitoring the health of an aluminum plate is introduced. The experimental results, which are used to demonstrate the correctness of the proposed method, are reported in Section 4 and discussed in Section 5. Finally, conclusions are presented in Section 6.

2. Method

2.1. Temperature Compensation Method for TOF

Guided waves refer to waves of strain and stress that propagate along a plate or a pipe with the traction-free boundary condition. Classified according to the displacement distribution in the propagation direction, there are two categories of guided waves: symmetrical and anti-symmetrical modes [31,32]. The A0 guided wave mode is a fundamental anti-symmetric mode of guided waves. Prior to temperature compensation, the TOFs of guided waves are assumed to have been estimated based on the sampled data. The envelope of the A0 guided wave mode is extracted using the Hilbert transform of the received signal. Then, the TOFs of guided waves are determined based on the maxima of the envelope.
It has been demonstrated in [28] that the TOF of a guided wave is linearly related to the temperature, which means that the TOF difference of two guided wave signals that are obtained at different temperatures is proportional to the temperature difference. As a result, one can express it as:
t 2 t 0 t 1 t 0 = T e m 2 T e m 0 T e m 1 T e m 0
where Tem0, Tem1 and Tem2 denote three temperature measurements of the inspected materials and t0, t1 and t2 represent the TOFs of the guided waves at temperatures Tem0, Tem1 and Tem2, respectively. The TOF shift of guided waves that is caused by temperature can be compensated as follows: According to Equation (1), the TOFs of guided waves at two temperatures and the corresponding temperatures should be measured prior to temperature compensation. Assume that the baseline signal is obtained in an undamaged inspected material at temperature Temb and the TOF that is extracted from the baseline signal is denoted as tb. The reference signal is obtained in the same undamaged inspected material at temperature Temr and the TOF that is extracted from the reference signal is denoted as tr. When the temperature varies, the baseline TOF is also changed. The baseline signal’s TOF, namely, tb, is compensated according to the temperature measurement Temc of an inspected material based on Equation (1); this can be expressed as:
t b c = T e m c T e m b T e m r T e m b ( t r t b ) + t b
where tbc is the baseline TOF at temperature Temc. If there is no defect in the sensing path, tbc will be equal to tc, which represents the TOF that was extracted from the inspected signal that was obtained at temperature Temc. In contrast, when the inspected signal has passed through a defect, tbc is not equal to tc and the TOF difference between tc and tbc is not 0. Therefore, when the inspection of plate-like material damage begins, the inspection temperature Temc is measured first. The baseline TOF at temperature Temb is compensated by the inspection temperature Temc based on Equation (2). After that, the TOF difference between the compensated baseline TOF and inspected TOF is used for damage imaging. In this way, the effect of temperature on damage imaging can be eliminated and Equation (2) is viewed as a temperature compensation method for TOF.

2.2. Temperature Compensation Model for Ray Tomography

Based on the proposed temperature compensation method for TOF, a temperature compensation model for ray tomography is developed. Assume that there are two transducer arrays with M-element transmitters and receivers that are mounted on a detected plate-like structure. The transducers in the transmitter array are individually used to transmit the guided waves, while all the receivers are exploited to receive the excited guided waves. The TOFs that are extracted from all the recorded data are used to perform thickness mapping.
The monitored area is divided into N × N grids (MN) and each propagation path between the transmitter and receiver is divided into small pieces by these grids. According to ray theory, the TOFs of guided waves in the M × M propagation paths can be expressed as the following linear Equation [9]:
T = L S
where T = ( t i ) M 2 × 1 is the column vector of the TOF measurements that are obtained from the received guided wave signals, S = ( s j ) N 2 × 1 is the column vector of slowness, and L = ( l i j ) M 2 × N 2 is the length matrix of the propagation paths, which are divided into pieces by the grids. Equation (3) is a conventional linear model for ray tomography. By solving Equation (3), one can reconstruct a velocity map of a guided wave and convert it to a thickness map. However, temperature variation is not considered in this model.
Suppose that there is no damage in the monitored plate-like material before its structure is investigated. The column vectors of TOF measurements that are obtained from the received guided wave signals are denoted as the baseline TOF measurements T b at temperature Temb and the reference TOF measurements T r at temperature Temr. Once T b and T r have been obtained, the system can begin its monitoring work. By inserting T b into Equation (3), one obtains:
T b = L S b
where S b is the column vector of slowness at temperature Temb. Then, the inspection temperature Temc is measured, and based on Equation (2), one obtains:
T b c = L S b c
where T b c is the column vector of the baseline TOF measurements at temperature Temc and S b c is the column vector of the baseline slowness at temperature Temc. Meanwhile, one may continuously or intermittently obtain the inspected TOF measurements T c that are described by Equation (3). Subtracting Equation (5) from Equation (3), one obtains:
T c T b c = L ( S c S b c )
Let Δ T = T c T b c and Δ S = S c S b c denote the TOF and slowness differences, respectively. One can rewrite Equation (6) as:
Δ T = L Δ S
Compared to Equation (3), Equation (7) contains a temperature parameter, which can be used to perform temperature compensation. Hence, Equation (7) is viewed as a temperature compensation model for ray tomography.
We assume that the plate damage gradually grows from a small partial thickness. Once a small partially monitored plate has been damaged, the slowness difference of the small partial grids changes. If K grids have been corroded, only the values of K elements in Δ S are non-zero. Since K is much smaller than N2, Δ S is sparse. In addition, the small partial grids with changed slowness difference are concentrated in the corroded area; hence, the non-zero elements in Δ S will appear continuously. Figure 1 presents such an example, where the monitored area is divided into 4 × 4 grids and four grids have been corroded. The solution of our sparse model from Equation (7) for this example is:
Δ S = ( 0   0   0   0 4 Δ s 1 Δ s 2 Δ s 3 Δ s 4 4 0   0   0   0 4 0   0   0   0 4 )
Four non-zero elements in Δ S appear continuously. If the monitored area is divided into 8 × 8 grids, sixteen non-zero elements in Δ S appear continuously. The grouping effect strengthens as the number of grids increases. Therefore, the solution of Equation (7) has a grouping effect except for the sparsity.

2.3. Elastic Net Penalty Approach

Since Equation (7) is also a linear equation, the simultaneous iterative reconstruction technique (SIRT) method [9], which is widely used in ray tomography, can be used to perform damage imaging. However, the solution that is obtained via SIRT is not sparse. Regularization methods can also be used to solve Equation (7) [33]. The solution that is obtained by the regularization methods is sparse; however, the grouping effect of the solutions, which is illustrated in Figure 1, is not taken into consideration.
In [34], a new regularization and variable selection method, namely, the elastic net penalty technique, is proposed. It often outperforms the abovementioned regularization methods, while enjoying similar sparsity of the representation. In addition, it is demonstrated in [34] that the solutions that are obtained via the elastic net penalty technique can cope with the grouping effect. For a nonnegative parameter α and a parameter β that is strictly between 0 and 1, the elastic net penalty technique solves Equation (7) as follows:
m i n Δ T L Δ S 2 2 + α ( 1 β 2 Δ S 2 + β Δ S 1 )
Function 1 β 2 Δ S 2 + β Δ S 1 is called the elastic net penalty. The elastic net penalty is the same as the Lasso penalty when β = 1. As β shrinks towards 0, the elastic net penalty approaches ridge regression. Therefore, the elastic net penalty is viewed as a convex combination of the Lasso and ridge penalties. For β ( 0 , 1 ] , it is strictly convex [34,35]. Since Equation (9) is a convex problem, if α and β are specified, the solutions of Equation (9) can be obtained using the CVX toolbox [36,37].
Parameters α and β are the two key factors in our method. The imaging quality is strongly related to them, especially regarding artefacts. Since α is a regularization parameter, the number of non-zero elements in Δ S decreases when α increases. In our applications, α and β are determined in turn. If β = 1, α ranges between 0.0001 and 0.01. Meanwhile, the optimal value of α can be calculated from the L-curve via Equation (9) in [38]. In our experiments, the optimal value of α that is calculated via this method is 0.001. Once α has been determined as discussed above, β can be calculated from the L-curve method when the function 1 β 2 Δ S 2 + β Δ S 1 is viewed as a whole. Via this approach, the calculated value of β is 0.3.

3. Experimental Setup

The experimental setup is shown in Figure 2a. Experimental measurements are performed on a 300 × 300 × 1.5 m3 aluminum plate. Piezoelectric transducers (PZTs) are employed as the ultrasound transducers in our experimental system, where the diameter of each PZT is 5 mm. The central frequency of the transducer is 250 kHz. The transmitter and receiver arrays each have 20 transducers. The transducer spacing is 10 mm for both the transmitter and receiver arrays. The distance between the two arrays is 200 mm. The aluminum plate with transducers is placed into a temperature testing box (custom made by Shan zhi Technologies, Wuhan, China), and its temperature control accuracy is 0.1 °C. The largest aluminum plate that is placed in the temperature testing box is 30 × 30 mm2.
During the system operation, the transducers in the transmitter array are used one by one to transmit a 3.5-cycle Hanning-windowed sinusoidal signal, while the receiver transducers are exploited to receive the excited guided waves. For example, when transmitter transducer 1 in Figure 2a is used at a specified time to excite the guided waves, the signal that is received by receiver transducer 1 is recorded by the oscilloscope. At another time, transmitter transducer 1 excites the guided waves again, and the received signal of receiver transducer 2 is sampled by the oscilloscope. This process is repeated until each receiver transducer has operated at least once. To ensure that each receiver transducer measures at the same time, a synchronous signal is outputted each time transmitter transducer 1 excites the guided waves. The oscilloscope sampling of the received signal of each receiver transducer is triggered by the synchronous signal. Then, the remaining transmitter transducers are used one by one to excite the guided waves; the above process is carried out repeatedly. The excitation signal is generated with an initial amplitude of 5 V using an arbitrary waveform generator (33220 A, which was manufactured by Agilent Technologies, Santa Clara, CA, USA) and subsequently amplified to 50 V using a power amplifier (AG1021, which was manufactured by T&C Power Conversion, Inc., Rochester, NY, USA). The received signal is sampled by an oscilloscope (54622A, which was manufactured by Agilent Technologies). The sampling rate of the acquisitions is 250 MS/s. A computer is used to read the received data from the oscilloscope. Figure 2b shows a typical original received signal, where both the transmitter and receiver transducers are the 10th ones. The received signal contains A0 and S0 guided wave modes. The amplitude of the A0 guided wave mode is greater than that of the S0 mode. To ensure imaging quality, a time gating function similar to that in [6] is applied to remove unwanted components, and the A0 guided wave mode is reserved for extracting the TOF. The envelopes of the excited signal and the A0 guided wave mode are extracted using the Hilbert transform. Then, the TOFs of the guided waves are determined based on the maxima of the envelope, as shown in Figure 2c. Next, the TOFs of the guided waves are employed for the inversion.
The flow of the experiment, which is illustrated in Figure 2d, is as follows: The defect-free aluminum plate is placed into the temperature testing box. The temperature testing box is set to 25 °C. The received signals in the defect-free aluminum plate are recorded, and their TOFs are extracted as the baseline TOFs. The temperature testing box is set to 65 °C. The received signals in the defect-free aluminum plate are recorded, and their TOFs are extracted as the reference TOFs. Then, defects are introduced into the aluminum plate, and the plate is placed into the temperature testing box. The temperature testing box is set to a temperature that differs from the baseline and reference temperatures. The received signals in the defective aluminum plate are recorded, and their TOFs are extracted as the inspected TOFs. Our temperature compensation model is used to perform temperature compensation. After that, the defects are imaged via Equation (9).
In our experiment, the wavelength of the A0 guided wave mode is approximately 7 mm. The size of the grid is selected as half of the wavelength. The monitored area is 200 × 200 mm2. Therefore, the monitored area is divided into 60 × 60 grids.

4. Experimental Results

4.1. TOF Versus Temperature Variation Results

To verify the relationship between the TOF and temperature, a defect-free aluminum plate is placed into the temperature testing box, where the temperature is gradually increased from 25 °C to 70 °C at a heating rate of 1 °C/min and for every 5 °C step, the temperature is held constant for 40 min. Then, the transmitter transducers are excited one by one with a 3.5-cycle Hanning-windowed sinusoidal signal with a central frequency of 250 kHz, and the receiver transducers record all the received signals using the oscilloscope. Take the 10th transmitter and receiver transducers as an example. The received signals at various temperatures in the range of 25 °C to 70 °C are illustrated in Figure 3.
The TOFs of the guided waves are influenced by the temperature changes. The TOF increases when the temperature is increased. The least-squares method was employed to fit the TOF as a function of temperature; the result is shown in Figure 4. The linearity of the TOF as a function of temperature with 95% confidence bounds is 0.9923. The result demonstrates that TOF has a strong linear relationship with temperature.

4.2. Single Regular Defect

To verify the correctness of the proposed method, an experiment with a single regular defect was conducted. The TOFs that were extracted from the received signals in the defect-free aluminum plate at a temperature of 25 °C were selected as the baseline TOFs and the selected reference temperature was 65 °C. A 50-mm-diameter circular flat-bottomed hole in the center of the same aluminum plate was made such that its thickness loss was 50%, as illustrated in Figure 5.
The aluminum plate was placed in the temperature testing box, and the temperature was set to 60 °C. This temperature is viewed as the inspection temperature. The tomographic images that were reconstructed using SIRT and our Equation (9) are shown in Figure 6. According to Figure 6a,b the spatial resolution of the tomographic reconstruction via SIRT is relatively low. In contrast, according to Figure 6c, the image of the tomographic reconstruction via our Equation (9) is relatively clear, even though temperature compensation was not performed. After temperature compensation, the visible artefacts that were caused by the temperature variations in the image are not observed, as shown in Figure 6d.

4.3. Two Regular Defects

To further verify the performance of the proposed method, an experiment with two regular defects was conducted. Another aluminum plate was used. The baseline temperature in this plate was selected as 25 °C, and the reference temperature was selected as 65 °C. Two 40-mm-diameter circular flat-bottomed holes were made in the same aluminum plate, as shown in Figure 7. The center coordinates of one the holes were x = 250 mm and y = 280 mm, and the thickness loss was 50%. The center coordinates of the remaining hole were x = 250 mm and y = 220 mm, and the thickness loss was 30%. This aluminum plate was placed in the temperature testing box and the temperature was set to 55 °C. That temperature is viewed as the inspection temperature. The tomographic images that were reconstructed using SIRT and our Equation (9) are shown in Figure 8. According to Figure 8a,b the tomographic reconstruction results that were obtained via SIRT confused the two defects with artefacts. In contrast, Figure 8c, which was obtained using the method that is proposed in this paper, clearly shows both thickness losses on the plate, even though temperature compensation was not performed. After temperature compensation, the visible artefacts in the image are no longer observed, as shown in Figure 8d. Moreover, the images of damage that were obtained via our method are more accurate and clearer than those via the SIRT method. The cross-sections of Figure 7 with the damage depths are shown in Figure 9. Our method outperforms SIRT in terms of the depth estimation results. After temperature compensation, the depth estimation accuracy was further improved using our method, as shown in Figure 9b.
The inspection temperatures were set to 25 °C, 55 °C and 70 °C, and the tomographic images that were reconstructed via Equation (9) are shown in Figure 10. According to Figure 10a–c, the number of artefacts significantly increases with the inspection temperature; this occurs because the baseline temperature is 25 °C, and the TOF is changed significantly when the temperature is increased. After temperature compensation, according to Figure 10d–f, the damage image is not influenced by the temperature variation.

5. Discussion

The results that are presented in this paper demonstrate that a guided wave that passes through plate-like materials with temperature changes is influenced by the temperature variation. Our approach is not only robust to the temperature variations but also yields more accurate and clearer thickness mappings of plate-like materials compared to SIRT. The visible artefacts in the tomographic reconstruction images are eliminated by our method. The statistical results from Figure 6 and Figure 8 are listed in Table 1. According to Table 1, the number of grids calculated by our method with temperature compensation is almost the same as the area of that number covered by the damage. Meanwhile, since there are many artefacts in the damage image that is reconstructed by the SIRT method, according to Figure 6 and Figure 8, the temperature effect on the SIRT method is not significant. However, from Table 1, the number of artefacts is significantly decreased using the SIRT method after temperature compensation. Moreover, the effectiveness of our temperature compensation method is demonstrated. In addition, the materials that can be inspected by our method are not limited to aluminum plates and the temperature compensation range is not limited to 25 °C to 70 °C. Our method is applicable when the TOF of guided waves that pass through the inspected material is a linear function of the temperature variation. The preliminary experiments are also performed on carbon fiber composites under a temperature variation from 25 °C to 70 °C, and the results demonstrate that the linearity of TOF as a function of temperature with 95% confidence bounds is 0.8563. Hence, the TOF of carbon fiber composites does not have a strong linear relationship with temperature. Therefore, ray tomography for the damage detection of composite materials under temperature variation should be further investigated.
In addition, the computation time of our approach was compared with that of SIRT. The computation times for SIRT and our method in the two experiments are listed in Table 2. The computer that was used in our experiment had an Intel core i7-4770 3.4 GHz quad-core desktop processor with 16 GB of DDR3 memory. According to Table 1, the computational time of our approach is slightly longer than that of SIRT.
The resolution of our method is also limited by the width of the first Fresnel zone, which is expressed as D λ , where λ is the wavelength of a guided wave and D is the distance between the transmitter and receiver transducers. In our experiment, the wavelength of the A0 guided wave mode was approximately 7 mm and the width of the first Fresnel zone was approximately 37 mm. To satisfy ray theory, the sizes of the manufactured defects were much larger than the wavelength, and the distance between the defects exceeded the width of the first Fresnel zone. Guided wave tomography methods have been available for more than 15 years [9]. Compared to a recent paper [18] that was published in Sensors, which also investigates damage monitoring using guided wave tomography, the imaging results of our method have poor resolution. However, our method is robust to temperature variations and has low computational cost. Our method can be used for real-time monitoring with temperature variations. In addition, only isotropic defects of regular shape have been considered in this paper. In the future, we will consider irregular defects and further analyze the accuracy of our method.

6. Conclusions

In this paper, we propose a new guided wave ray tomography method that is based on TOF temperature compensation for temperature-robust damage detection. Our approach does not require the estimation of temperature compensation parameters to perform temperature compensation on the guided wave waveform. Our approach only requires the temperature measurement of an inspected material to perform temperature compensation based on the linear relationship between the TOFs of guided waves and the temperature variations. Then, we formulate an optimization problem of the difference between the inspected and baseline TOFs of the guided waves that pass through plate-like materials with temperature variations via an elastic net penalty approach to perform damage mapping of an SHM application via ray tomography. The correctness of our method has been demonstrated experimentally. Additionally, our approach is capable of reconstructing damage mappings of materials without artefacts when the temperature varies from 25 °C to 70 °C. Moreover, the proposed method achieves more accurate and clearer imaging of damage than the SIRT method.

Author Contributions

conceptualization, D.L. and J.Z.; methodology, D.L. and D.T.; software, M.S.; validation, M.S., F.X. and D.L.; writing—original draft preparation, D.L.; writing—review and editing, D.L., C.L., J.Z. and D.T.

Funding

This research was funded by the National Natural Science Foundation of China, grants number 11604055, 61571131, and 11604054.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Groysman, A. Corrosion for Everybody; Springer: New York, NY, USA, 2010. [Google Scholar]
  2. Hai-Yan, Z.; Xian-Hua, C.; Ya-Ping, C.; Jian-Bo, Y. Focusing of Time Reversal Lamb Waves and Its Applications in Structural Health Monitoring. Chin. Phys. Lett. 2010, 27, 138–141. [Google Scholar] [CrossRef]
  3. Li, W.B.; Cho, Y. Combination of nonlinear ultrasonics and guided wave tomography for imaging the micro-defects. Ultrasonics 2016, 65, 87–95. [Google Scholar] [CrossRef] [PubMed]
  4. Lin, J.; Gao, F.; Luo, Z.; Zeng, L. High-Resolution Lamb Wave Inspection in Viscoelastic Composite Laminates. IEEE Trans. Ind. Electron. 2016, 63, 6989–6998. [Google Scholar] [CrossRef]
  5. Qiu, L.; Yuan, S.; Mei, H.; Fang, F. An Improved Gaussian Mixture Model for Damage Propagation Monitoring of an Aircraft Wing Spar under Changing Structural Boundary Conditions. Sensors 2016, 16, 291. [Google Scholar] [CrossRef] [PubMed]
  6. Hinders, M.K.; Hou, J.; Leonard, K.R. Multi-Mode Lamb Wave Arrival Time Extraction for Improved Tomographic Reconstruction. In Proceedings of the AIP Conference, Salt Lake City, UT, USA, 10−11 August 2005; pp. 736–743. [Google Scholar]
  7. Volker, A.; Vos, H. Experimental results of guided wave travel time tomography. In Proceedings of the AIP Conference Proceedings, Jyväskylä, Finland, 3–7 September 2012; pp. 1968–1975. [Google Scholar]
  8. Willey, C.L.; Simonetti, F. A two-dimensional analysis of the sensitivity of a pulse first break to wave speed contrast on a scale below the resolution length of ray tomography. J. Acoust. Soc. Am. 2016, 139, 3145–3158. [Google Scholar] [CrossRef] [PubMed]
  9. Malyarenko, E.V.; Hinders, M.K. Fan beam and double crosshole Lamb wave tomography for mapping flaws in aging aircraft structures. J. Acoust. Soc. Am. 2000, 108, 1631–1639. [Google Scholar] [CrossRef] [PubMed]
  10. Belanger, P.; Cawley, P.; Simonetti, F. Guided wave diffraction tomography within the born approximation. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2010, 57, 1405–1418. [Google Scholar] [CrossRef] [PubMed]
  11. Malyarenko, E.V.; Hinders, M.K. Ultrasonic Lamb wave diffraction tomography. Ultrasonics 2001, 39, 269. [Google Scholar] [CrossRef]
  12. Huthwaite, P.; Simonetti, F. High-resolution imaging without iteration: A fast and robust method for breast ultrasound tomography. J. Acoust. Soc. Am. 2011, 130, 1721–1734. [Google Scholar] [CrossRef] [PubMed]
  13. Huthwaite, P.; Simonetti, F. High resolution guided wave tomography. Wave Motion 2013, 50, 979–993. [Google Scholar] [CrossRef]
  14. Rao, J.; Ratassepp, M.; Fan, Z. Limited-view ultrasonic guided wave tomography using an adaptive regularization method. J. Appl. Phys. 2016, 120, 113–127. [Google Scholar] [CrossRef]
  15. Rao, J.; Ratassepp, M.; Fan, Z. Guided Wave Tomography Based on Full-Waveform Inversion. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2016, 63, 737–745. [Google Scholar] [CrossRef] [PubMed]
  16. Rao, J.; Ratassepp, M.; Fan, Z. Investigation of the reconstruction accuracy of guided wave tomography using full waveform inversion. J. Sound Vib. 2017, 400, 317–328. [Google Scholar] [CrossRef]
  17. Rao, J.; Ratassepp, M.; Fan, Z. Quantification of thickness loss in a liquid-loaded plate using ultrasonic guided wave tomography. Smart Mater. Struct. 2017, 26, 125017. [Google Scholar] [CrossRef]
  18. Rao, J.; Ratassepp, M.; Lisevych, D.; Hamzah Caffoor, M.; Fan, Z. On-Line Corrosion Monitoring of Plate Structures Based on Guided Wave Tomography Using Piezoelectric Sensors. Sensors 2017, 17, 2882. [Google Scholar] [CrossRef] [PubMed]
  19. Cawley, P.; Cegla, F.; Galvagni, A. Guided waves for NDT and permanently-installed monitoring. Insight Non-Destr. Test. Condition Monit. 2012, 54, 594–601. [Google Scholar] [CrossRef]
  20. Liu, G.; Xiao, Y.; Zhang, H.; Ren, G. Baseline Signal Reconstruction for Temperature Compensation in Lamb Wave-Based Damage Detection. Sensors 2016, 16, 1273. [Google Scholar] [CrossRef] [PubMed]
  21. Michaels, J.E.; Michaels, T.E. Detection of structural damage from the local temporal coherence of diffuse ultrasonic signals. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2005, 52, 1769–1782. [Google Scholar] [CrossRef] [PubMed]
  22. Harley, J.B.; Moura, J.M. Scale transform signal processing for optimal ultrasonic temperature compensation. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2012, 59, 2226. [Google Scholar] [CrossRef] [PubMed]
  23. Lu, Y.; Michaels, J.E. A methodology for structural health monitoring with diffuse ultrasonic waves in the presence of temperature variations. Ultrasonics 2005, 43, 717. [Google Scholar] [CrossRef] [PubMed]
  24. Konstantinidis, G.; Drinkwater, B.W.; Wilcox, P.D. The temperature stability of guided wave structural health monitoring systems. Smart Mater. Struct 2006, 15, 967. [Google Scholar] [CrossRef]
  25. Clarke, T.; Cawley, P.; Wilcox, P.D.; Croxford, A.J. Evaluation of the damage detection capability of a sparse-array guided-wave SHM system applied to a complex structure under varying thermal conditions. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2009, 56, 2666. [Google Scholar] [CrossRef] [PubMed]
  26. Clarke, T.; Simonetti, F.; Cawley, P. Guided wave health monitoring of complex structures by sparse array systems: Influence of temperature changes on performance. J. Sound Vib 2010, 329, 2306–2322. [Google Scholar] [CrossRef]
  27. Croxford, A.J.; Moll, J.; Wilcox, P.D.; Michaels, J.E. Efficient temperature compensation strategies for guided wave structural health monitoring. Ultrasonics 2010, 50, 517–528. [Google Scholar] [CrossRef] [PubMed]
  28. Wang, Y.; Gao, L.; Yuan, S.; Qiu, L.; Qing, X. An adaptive filter-based temperature compensation technique for structural health monitoring. J. Intell. Mater. Syst. Struct. 2014, 25, 2187–2198. [Google Scholar] [CrossRef]
  29. Lim, Y.Y.; Kwong, K.Z.; Liew, W.Y.H.; Padilla, R.V.; Soh, C.K. Parametric study and modeling of PZT based wave propagation technique related to practical issues in monitoring of concrete curing. Constr. Build. Mater. 2018, 176, 519–530. [Google Scholar] [CrossRef]
  30. Konstantinidis, G.; Wilcox, P.D.; Drinkwater, B.W. An Investigation into the Temperature Stability of a Guided Wave Structural Health Monitoring System Using Permanently Attached Sensors. IEEE Sens. J. 2007, 7, 905–912. [Google Scholar] [CrossRef]
  31. Xu, K.; Ta, D.; Su, Z.; Wang, W. Transmission analysis of ultrasonic Lamb mode conversion in a plate with partial-thickness notch. Ultrasonics 2014, 54, 395–401. [Google Scholar] [CrossRef] [PubMed]
  32. Xu, K.; Laugier, P.; Minonzio, J.G. Dispersive Radon Transform. J. Acoust. Soc. Am. 2018, 143, 2729–2743. [Google Scholar] [CrossRef] [PubMed]
  33. Hastie, T.; Tibshirani, R.; Wainwright, M. Statistical Learning with Sparsity: The Lasso and Generalizations; CRC Press: Boca Raton, FL, USA, 2015; pp. 156–157. [Google Scholar]
  34. Zou, H.; Hastie, T. Regularization and Variable Selection via the Elastic Net. J. R. Stat. Soc. 2005, 67, 768. [Google Scholar] [CrossRef]
  35. Boyd, S.; Vandenberghe, L. Convex Optimization; Cambridge University Press: New York, NY, USA, 2004; p. 1859. [Google Scholar]
  36. Grant, M.C.; Boyd, S.P. Graph Implementations for Nonsmooth Convex Programs. In Recent Advances in Learning and Control; Springer: London, UK, 2008; pp. 95–110. [Google Scholar]
  37. Grant, M.C.; Boyd, S.P. CVX: Matlab Software for Disciplined Convex Programming. Available online: http://cvxr.com/cvx/ (accessed on 18 October 2018).
  38. Oraintara, S.; Karl, W.C.; Castanon, D.A.; Nguyen, T.Q. A method for choosing the regularization parameter in generalized Tikhonov regularized linear inverse problems. In Proceedings of the International Conference on Image Processing, Vancouver, BC, Canada, 10–13 September 2000; pp. 93–96. [Google Scholar]
Figure 1. Example of the monitored area, which is divided into 4 × 4 grids.
Figure 1. Example of the monitored area, which is divided into 4 × 4 grids.
Sensors 18 03518 g001
Figure 2. (a) SHM system diagram. (b) A typical time trace from the experiment within the chosen time window. (c) An example of TOF extraction from the transmitter and the receiver. (d) A flowchart of the experiment.
Figure 2. (a) SHM system diagram. (b) A typical time trace from the experiment within the chosen time window. (c) An example of TOF extraction from the transmitter and the receiver. (d) A flowchart of the experiment.
Sensors 18 03518 g002aSensors 18 03518 g002b
Figure 3. Received signals of the 10th transmitter and receiver transducers at various temperatures in the range from 25 °C to 70 °C.
Figure 3. Received signals of the 10th transmitter and receiver transducers at various temperatures in the range from 25 °C to 70 °C.
Sensors 18 03518 g003
Figure 4. Fit of TOF versus temperature from 25 °C to 70 °C.
Figure 4. Fit of TOF versus temperature from 25 °C to 70 °C.
Sensors 18 03518 g004
Figure 5. Aluminum plate with a thickness of 1.5 mm and a 50-mm-diameter circular flat-bottomed hole. The thickness loss is 50% within the flaw.
Figure 5. Aluminum plate with a thickness of 1.5 mm and a 50-mm-diameter circular flat-bottomed hole. The thickness loss is 50% within the flaw.
Sensors 18 03518 g005
Figure 6. Tomographic reconstruction images for a 1.5-mm-thick aluminum plate with a 50-mm-diameter circular flat-bottomed hole. The thickness loss is 50% within the flaw. (a) The tomographic reconstruction image that is obtained via SIRT without temperature compensation, (b) the tomographic reconstruction image that is obtained via SIRT with temperature compensation, (c) the tomographic reconstruction image that is obtained via our method without temperature compensation, and (d) the tomographic reconstruction image that is obtained via our method with temperature compensation.
Figure 6. Tomographic reconstruction images for a 1.5-mm-thick aluminum plate with a 50-mm-diameter circular flat-bottomed hole. The thickness loss is 50% within the flaw. (a) The tomographic reconstruction image that is obtained via SIRT without temperature compensation, (b) the tomographic reconstruction image that is obtained via SIRT with temperature compensation, (c) the tomographic reconstruction image that is obtained via our method without temperature compensation, and (d) the tomographic reconstruction image that is obtained via our method with temperature compensation.
Sensors 18 03518 g006
Figure 7. Aluminum plate with a thickness of 1.5 mm and two 40-mm-diameter circular flat-bottomed holes. Their thickness losses are 50% and 30% within the flaws.
Figure 7. Aluminum plate with a thickness of 1.5 mm and two 40-mm-diameter circular flat-bottomed holes. Their thickness losses are 50% and 30% within the flaws.
Sensors 18 03518 g007
Figure 8. Tomographic reconstruction images for a 1.5-mm-thick aluminum plate with 40-mm-diameter circular flat-bottomed holes. Their thickness losses are 50% and 30% within the flaws. (a) The tomographic reconstruction image that is obtained via SIRT without temperature compensation, (b) the tomographic reconstruction image that is obtained via SIRT with temperature compensation, (c) the tomographic reconstruction image that is obtained via our method without temperature compensation, and (d) the tomographic reconstruction image that is obtained via our method with temperature compensation.
Figure 8. Tomographic reconstruction images for a 1.5-mm-thick aluminum plate with 40-mm-diameter circular flat-bottomed holes. Their thickness losses are 50% and 30% within the flaws. (a) The tomographic reconstruction image that is obtained via SIRT without temperature compensation, (b) the tomographic reconstruction image that is obtained via SIRT with temperature compensation, (c) the tomographic reconstruction image that is obtained via our method without temperature compensation, and (d) the tomographic reconstruction image that is obtained via our method with temperature compensation.
Sensors 18 03518 g008
Figure 9. Cross-sections of the reconstructions of two circular flat-bottomed holes along the line that is marked in Figure 7. (a) SIRT and our method without temperature compensation and (b) SIRT and our method with temperature compensation.
Figure 9. Cross-sections of the reconstructions of two circular flat-bottomed holes along the line that is marked in Figure 7. (a) SIRT and our method without temperature compensation and (b) SIRT and our method with temperature compensation.
Sensors 18 03518 g009
Figure 10. Tomographic reconstruction images for a 1.5-mm-thick aluminum plate with 40-mm-diameter circular flat-bottomed holes. Their thickness losses are 50% and 30% within the flaws. (a) The inspection temperature is 25 °C, and temperature compensation is not performed for the tomographic reconstruction image; (b) the inspection temperature is 55 °C, and temperature compensation is not performed for the tomographic reconstruction image; (c) the inspection temperature is 70 °C, and temperature compensation is not performed for the tomographic reconstruction image; (d) the inspection temperature is 25 °C, and temperature compensation is performed for the tomographic reconstruction image; (e) the inspection temperature is 55 °C, and temperature compensation is performed for the tomographic reconstruction image; and (f) the inspection temperature is 70 °C, and temperature compensation is performed for the tomographic reconstruction image.
Figure 10. Tomographic reconstruction images for a 1.5-mm-thick aluminum plate with 40-mm-diameter circular flat-bottomed holes. Their thickness losses are 50% and 30% within the flaws. (a) The inspection temperature is 25 °C, and temperature compensation is not performed for the tomographic reconstruction image; (b) the inspection temperature is 55 °C, and temperature compensation is not performed for the tomographic reconstruction image; (c) the inspection temperature is 70 °C, and temperature compensation is not performed for the tomographic reconstruction image; (d) the inspection temperature is 25 °C, and temperature compensation is performed for the tomographic reconstruction image; (e) the inspection temperature is 55 °C, and temperature compensation is performed for the tomographic reconstruction image; and (f) the inspection temperature is 70 °C, and temperature compensation is performed for the tomographic reconstruction image.
Sensors 18 03518 g010
Table 1. Number of grids that are covered by the damage area compared to the numbers that are calculated by the SIRT and proposed methods.
Table 1. Number of grids that are covered by the damage area compared to the numbers that are calculated by the SIRT and proposed methods.
DefectActual Damage AreaSIRT without T. C. *SIRT with T. C.Proposed Method without T. C.Proposed Method with T. C.
Single defect177961693236186
Two defects224992732332248
* T. C. represents temperature compensation.
Table 2. Computation times of SIRT and our method (UNITS: s).
Table 2. Computation times of SIRT and our method (UNITS: s).
DefectSIRT without T. C. *SIRT with T. C.Proposed Method without T. C.Proposed Method with T. C.
Single defect0.010.010.890.90
Two defects0.020.020.970.99

Share and Cite

MDPI and ACS Style

Li, D.; Shi, M.; Xu, F.; Liu, C.; Zhang, J.; Ta, D. A New Approach to Guided Wave Ray Tomography for Temperature-Robust Damage Detection Using Piezoelectric Sensors. Sensors 2018, 18, 3518. https://doi.org/10.3390/s18103518

AMA Style

Li D, Shi M, Xu F, Liu C, Zhang J, Ta D. A New Approach to Guided Wave Ray Tomography for Temperature-Robust Damage Detection Using Piezoelectric Sensors. Sensors. 2018; 18(10):3518. https://doi.org/10.3390/s18103518

Chicago/Turabian Style

Li, Dan, Ming Shi, Feng Xu, Chengcheng Liu, Jianqiu Zhang, and Dean Ta. 2018. "A New Approach to Guided Wave Ray Tomography for Temperature-Robust Damage Detection Using Piezoelectric Sensors" Sensors 18, no. 10: 3518. https://doi.org/10.3390/s18103518

APA Style

Li, D., Shi, M., Xu, F., Liu, C., Zhang, J., & Ta, D. (2018). A New Approach to Guided Wave Ray Tomography for Temperature-Robust Damage Detection Using Piezoelectric Sensors. Sensors, 18(10), 3518. https://doi.org/10.3390/s18103518

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