A modified damage index probability imaging algorithm based on delayandsum imaging for synthesizing timereversed Lamb waves
Qingwei Xia^{1} , Yanyan Liu^{2} , Yu Lu^{3} , Shuhao Cao^{4} , Hanfei Zhang^{5} , Shiwei Ma^{6}
^{1, 2, 3, 4, 5, 6}Shanghai University, Shanghai, China
^{6}Corresponding author
Journal of Vibroengineering, Vol. 21, Issue 8, 2019, p. 21402147.
https://doi.org/10.21595/jve.2019.20568
Received 1 February 2019; received in revised form 2 June 2019; accepted 12 June 2019; published 31 December 2019
JVE Conferences
Imaging for damage in plate structure by Lamb waves is one of the most effective methods in the field of structural health monitoring. In order to improve the accuracy of damage localization, a novel method is proposed to modify damage exponent probability imaging algorithm based on delayandsum imaging by using time reversal Lamb waves. A new probability distribution function is introduced to improve the damage index probability method and is combined with delayandsum method for damage localization. Experimental results on aluminum plate show that the hybrid algorithm achieves better accuracy of damage location and imaging quality than the conventional delayandsum method.
Keywords: damage probability, improved distribution function, time reversal, delayandsum.
1. Introduction
Structural health monitoring (SHM) has been an important research field in material integrity nondestructive testing (NDT). As a special ultrasonic wave, Lamb wave has been widely used in the health monitoring of thin plate or shell structure in recent years because of its long propagation distance and low energy attenuation. When the Lamb wave propagates through the defect of the structure, the scattering phenomenon occurs. The damage can be regarded as a secondary source, and the scattered signal propagates in all directions, which pass through the defect again and form a larger echo [13].
Time reversal based on the principle of acoustic reciprocity can focus energy in space and time, through which the position of sound source can be obtained, and the signal of sound source can be reconstructed. Ing and Fink [4] compensate the dispersion of Lamb wave automatically by time reversal method. The experiments show the spatial and temporal focusing characteristics of Lamb wave. Wang et al. [5] synthesizes time reversal and synthetic aperture techniques for damage image in structures.
Damage Index (DI) is proposed to characterize changes caused by structural damage and identify the existence of the damage. Wang D et al. [6] proposes a probabilistic imaging algorithm for damage index based on sensor network. Wang et al. [5] proposes a delayandsum (DAS) algorithm which effectively images the received signals but results in artifacts for multiple defects. Sohn et al. [2] studies the relationship between the damage condition and the consistency of shape of the source and the reconstructed signal and uses the damage index to characterize the damage degree. Hay et al. [7] proposes a reconstruction algorithm for probabilistic inspection of damage (RAPID), which is a probabilitybased damage imaging method.
In this paper, we propose a modified damage index probability imaging algorithm based on delayandsum imaging algorithm by synthesizing timereversed Lamb waves. A new probability distribution is investigated to better describe the probability of defect occurrence. The validity of the new algorithm is verified by experiments and the damage location is accurately displayed.
2. Principle of time reversal method
Time reversal means that the receiver receives a signal transmitted by the sound source in time domain, and then transmits the signal from the corresponding receiving sensor, that is, first to the second, and then to the first. By focusing energy in space and time, the position of sound source can be obtained and the reconstruction of signal from sound source can be realized.
We refine the above process. First, the electrical signal is converted to mechanical strain through the transmitter, where other factors are ignored, and only piezoelectric conversion is considered [8]:
where $\omega $ is angular frequency, ${k}_{a}\left(\omega \right)$ is electromechanical efficiency coefficient, $V\left(\omega \right)$ is driving voltage, and $\epsilon \left(\omega \right)$ is mechanical strain. The Lamb waves generated by the mechanical strain propagates in the medium, and the signals received at the receiver are as follows:
where ${k}_{s}\uff08\omega \uff09$ is mechanicalelectro efficiency coefficient and $G\left(\omega \right)$ is a transfer function. The time reversal of the signal ${V}_{B}\left(\omega \right)$ is equivalent to the complex conjugation ${V}_{B}^{*}\left(\omega \right)$ in frequency domain. Then the signal that reverses the time is transmitted from the receiver, at which point the signal received at the transmitter is the refactoring signal:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=\mathrm{}{k}_{s}\uff08\omega \uff09\cdot G\left(\omega \right)\cdot {k}_{a}\uff08\omega \uff09\cdot {k}_{s}^{*}\uff08\omega \uff09\cdot {G}^{*}\left(\omega \right)\cdot {k}_{a}^{*}\uff08\omega \uff09\cdot {V}^{*}\left(\omega \right)$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=T\left(\omega \right)\cdot K\left(\omega \right)\cdot {K}^{*}\left(\omega \right)\cdot {V}^{*}\left(\omega \right),$
where $G\left(\omega \right)$ is the same transfer function as in Eq. (2) based on the reciprocity of elastodynamics [9]. For simplification, let $\left(\omega \right)=G\left(\omega \right)\cdot {G}^{*}\left(\omega \right)$ and $K\left(\omega \right)={k}_{a}\uff08\omega \uff09{\cdot k}_{s}\uff08\omega \uff09$.
For thin plate without defects, $G\left(\omega \right)$ is a constant, so the reconstructed signal is similar to the original signal. For the case of thin plate with defect,$\mathrm{}G\left(\omega \right)$ is nonlinear, and the frequency dispersion occurs, so the reconstructed signal is distorted compared with the excitation signal after the final time reversal [10], as shown in Fig. 1. In this paper, the reconstructed signals are normalized and enveloped.
Fig. 1. Reconstruction signal dispersion
3. Proposed method
3.1. Delayandsum imaging
The delay and sum imaging algorithm is a combination of elliptic imaging algorithm and time reversal theory. The position of damage is the intersection of multiple ellipses as shown in Fig. 2. In detection region with damage, the time of the point $X\left(x,y\right)$ can be calculated by Eq. (4):
where ${S}_{0}$ is the transmitting sensor, ${S}_{i}$ is the $i$th receiving sensor, and ${c}_{g}$ is Lamb wave group velocity.$\mathrm{}\Vert {S}_{0}X\Vert $ is the distance from ${S}_{0}$ to $X\left(x,y\right)$.$\mathrm{}\Vert {S}_{i}X\Vert $ is the distance from ${S}_{i}$ to $X\left(x,y\right)$. The pixel value of point $X\left(x,y\right)$ can be calculated by Eq. (5) for a known time reversal of the reconstruction signal:
where $I(x,y)$ is the damage imaging value of the position $X\left(x,y\right)$, $N$ is the number of the signals, ${A}_{i}$ is a coefficient for balancing the transducers’ output,$\mathrm{}t(x,y)$ is the time from point $X\left(x,y\right)$ to a set of transmitterreceive pairs, ${f}_{i}$ is the reconstruction signal, and ${t}_{s}$ is delay time.
Fig. 2. Delayandsum imaging method
3.2. Improved damage index probability imaging
According to time reversal method, the damage index of a sensing path can be defined as the similarity between the reconstructed signal from this transmitting  receiving pair and the original excitation signal. There are several methods to calculate the damage index of a sensing path, and in this paper the classical method [11] is used. The damage index reflects the influence of defect on the sensing path, as follows:
where $i$ denotes the $i$th damage sensing path, ${t}_{0}$ and ${t}_{1}$ are the signal start time and the end time, respectively. The probability of defect occurrence at a certain point can be reconstructed from a ${DI}_{i}$ mapped probability distribution function. For $N$ sensing paths in detection region, the damage probability of each pixel $X(x,y)$ can be expressed as Eq. (7):
where $P(x,y)$ is the defect probability at position $\left(x,y\right)$ within detection region,$\mathrm{}{p}_{i}\left(x,y\right)$ is the estimation from the$i$th sensing path, and ${W}_{i}\left[{R}_{i}\left(x,y\right)\right]$ is the probability distribution function. In this paper, the distribution function is improved on the basis of the study of Sheen et al. [11]. The relative distance between the transmitterreceiver pair ${R}_{i}\left(x,y\right)$ is defined as:
where ${d}_{Ai}$ is the distance from $X\left(x,y\right)$ to the transmitter in the $i$th sensing path. Similarly, ${d}_{Ri}$ is the distance to the receiver. ${d}_{i}$ is the distance from the transmitter to receiver. $m$ is an index of distribution decline rate which determines the decreasing rate of damage index in elliptical distribution. By controlling the descent rate of the distribution function, the corresponding distribution boundary is controlled by Eq. (9):
where $\beta $ is a shape parameter of a certain sensing path, which determines the distribution boundary. The best distribution function in DAS weighted damage probability index algorithm is explored for the optimum combination of $m$ and $\beta $.
3.3. Improved DAS weighted damage probability imaging
In DAS method, damage is located at where the amplitude of ellipse is maximum. Therefore, the distribution function of damage probability is further improved. According to Eq. (5) and Eq. (9), the improved distribution function is expressed as Eq. (10):
So the overall image by the proposed algorithm is as follows:
4. Experiment and results
4.1. Experimental setup
Experimental instrument consists of a host computer, an arbitrary waveform signal generator (Agilent 33520A), an independently designed active signal amplifier (based on the TI op amp chip OPA657), and a 4channel oscilloscope (Tektronix DPO5034B). The system is shown in Fig. 3.
Fig. 3. Schematic illustration of experimental setup
a) Experimental system
b) Experiment setup
Two pieces of 1000×1000×1.5 mm^{3} aluminum plate with a single circular hole defect and double circular hole defects, respectively, are used in the experiment. The diameter of all the defects is 15 mm. Taking the lower left corner of aluminum plate as the origin, the coordinate of single circular hole defect is (550, 550), the coordinates of double circular hole defects are (550, 550) and (450, 500), and the unit is $\mathrm{m}\mathrm{m}$, as shown in Figs. 4(a) and 4(b).
Fig. 4. Diagram of sample plates
a) Single flaw sample
b) Double flaws sample
c) Sensing paths layout
At the same time, the strategy of limited detection region is used to physically avoid the influence of edge reflection. The detection region is divided into 400×400 grids with an area of 2.5×2.5 mm^{2} for each grid. The arrangement of eight sensors $\left({S}_{1}~{S}_{8}\right)\mathrm{}$is shown in Fig. 4(c), with coordinates of (mm):$\mathrm{}{S}_{1}\left(\mathrm{250,250}\right)$, ${S}_{2}\left(\mathrm{500,250}\right)$, ${S}_{3}\left(\mathrm{750,250}\right)$, ${S}_{4}\left(\mathrm{500,750}\right)$, ${S}_{5}\left(\mathrm{750,750}\right)$, ${S}_{6}\left(\mathrm{750,500}\right)$, ${S}_{7}\left(\mathrm{750,250}\right)$, and ${S}_{8}\left(\mathrm{500,250}\right)$. The sensor (Bat brand by Changzhou Ultrasonic Electronics Co. Ltd.) with an oblique probe with a dip angle of 30° is used in the experiment.
The attenuation of Lamb wave is inevitable in the transmission process. In order to reduce the attenuation on the reconstructed signal, a 5cycle Hanning windowmodulated sine wave is used as the excitation signal, and the frequency is 500 kHz.
The time reversal reconstruction signal of all the paths as shown in Fig. 4(c) is obtained. All signals are normalized and enveloped. The parameter $m$ is set to 1/2, 1 and 2, and the optimum imaging combination is found by the experimental regulation of parameter $\beta $.
4.2. Data processing and experimental results
According to Eq. (6), the damage index along different sensing paths for single flaw is shown in Table 1, and for double flaws is shown in Table 2. Six sensing paths along which the amplitudes of reconstructed signals are the first six largest are selected for imaging by Eq. (11).
Table 1. Damage index of single flaw
Paths

${S}_{4}{S}_{7}$

${S}_{1}{S}_{5}$

${S}_{1}{S}_{6}$

${S}_{3}{S}_{6}$

${S}_{3}{S}_{7}$

DI

0.1402

0.1336

0.1330

0.1287

0.1242

Paths

${S}_{3}{S}_{8}$

${S}_{2}{S}_{5}$

${S}_{1}{S}_{4}$

${S}_{2}{S}_{7}$

${S}_{5}{S}_{8}$

DI

0.1212

0.1202

0.1189

0.1160

0.1010

Table 2. Damage index of double flaws
Paths

${S}_{4}{S}_{7}$

${S}_{1}{S}_{5}$

${S}_{3}{S}_{7}$

${S}_{1}{S}_{4}$

${S}_{3}{S}_{6}$

DI

0.1695

0.1450

0.1429

0.1386

0.1322

Paths

${S}_{3}{S}_{8}$

${S}_{1}{S}_{6}$

${S}_{2}{S}_{7}$

${S}_{2}{S}_{5}$

${S}_{5}{S}_{8}$

DI

0.1299

0.1228

0.1192

0.1167

0.1129

For singleflaw sample, with distribution drop rate parameter $m$ = 1/2, 1and 2, the best shape parameter is calculated as $\beta =$ 0.974, 0.948 and 0.895, respectively. The corresponding images are shown in Figs. 5(a), 5(b) and 5(c), where the same black circle is the actual damage position, and the white circle is the damage location estimated by the proposed algorithm. It is seen that the three sets of parameter combination of $m$ and $\beta $ show a good localization. The estimated damage position is (557.5, 562.5). In order to facilitate observation, A 90 % threshold value is used for Fig. 5(c) for better illustration, and the enlarged image is as shown in Fig. 5(d), where the black circle is the actual damage position, and white circle is the estimated damage position with the largest pixel value.
For doubleflaw sample, the parameter combination of $m$ and $\beta $ are the same as for singleflaw. The images for doubleflaw are shown in Fig. 6(a), 6(b) and 6(c). For $m$ = 1/2 and $\beta $ = 0.974, the estimated positions of the two defects are (560, 550) and (435, 512.5), respectively. For $m=$ 1 and $\beta $ = 0.948, the estimated position of damage 1 is (560, 550); there are two possible positions for damage 2, which are (435, 512.5) or (440, 505). For $m$ = 2 and $\beta $ = 0.895, the image of double flaws shows the best localization. The estimated positions of double flaws are (560, 550) and (440, 505), respectively. An enlarged image by applying a 90 % threshold value to Fig. 6(c) is shown in Fig. 6(d), where the black circle is the actual damage position, and the white circle is the estimated damage position with the largest pixel value.
The proposed method is compared with DAS imaging using the same setting. The images of DAS method are shown in Fig. 7(a) and 7(b).
Fig. 5. Imaging results of singleflaw sample
a)$m$ = 1/2,$\beta $ = 0.974
b)$m$ = 1,$\beta $ = 0.948
c)$m$ = 2,$\beta $ = 0.895
Fig. 6. Imaging results of doubleflaw sample
a)$m$ = 1/2,$\beta $ = 0.974
b)$m$ = 1,$\beta $ = 0.948
c)$m$ = 2,$\beta $ = 0.895
Fig. 7. Imaging results of DAS
a) Singleflaw sample
b) Doubleflaw sample
The errors of our proposed method and DAS method for single flaw and double flaws sample are listed in Table 3. Generally, our proposed method has a better damage localization than DAS method for both singleflaw and double flaw cases. For parameter combination of $m=$ 2 and $\beta =$ 0.895, our proposed method shows the most accurate damage position.
Table 3. Error of our proposed algorithm and DAS method
Algorithm

Single defect

Double defects


D1

D1

D2


Our proposed method ($\mathrm{m}\mathrm{}$= 2,$\mathrm{}\mathrm{\beta}\mathrm{}$= 0.895)

14.58 mm

10.00 mm

11.18 mm

Delayandsum

19.04 mm

29.15 mm

40.31 mm

5. Conclusions
Based on the delayandsum weighted damage index probability method, the proposed method improves the distribution function by introducing a distribution drop rate parameter $m$. The time reversal lamb wave is used to focus the scattering energy, and the damage index is obtained by calculating the reconstruction deviation from the original signal. Experimental result shows that the method has better location accuracy compared with the delayandsum weighted damage index probability method, which validates the efficiency of the proposed method, and the best parameter combination are distribution drop rate parameter $m$ = 2 and the shape parameter $\beta $ = 0.895. In future studies, more precise positioning can be achieved through more accurate group velocity measurements.
Acknowledgements
This work was supported by the National Natural Science Foundation of China with Grant Nos. 61671285.
References
 Xu B., Giurgiutiu V. Single mode tuning effects on lamb wave time reversal with piezoelectric wafer active sensors for structural health monitoring. Journal of Nondestructive Evaluation, Vol. 26, Issues 24, 2007, p. 123134. [Publisher]
 Park H. W., Kim S. B., Sohn H. Understanding a time reversal process in Lamb wave propagation. Wave Motion, Vol. 46, Issue 7, 2009, p. 451467. [Publisher]
 Liqiang G., Weiguang Z., Yifeng L. Research on hybrid techniques of timereversal ellipse location and tomographic imaging of lamb wave. Journal of Nanjing University (Natural Science), 2019. [CrossRef]
 Ing R. K., Fink M. Timereversed Lamb waves. IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency, Control, Vol. 45, Issue 4, 1998, p. 10321043. [Publisher]
 Wang C. H., Rose J. T., Chang F. K. A synthetic timereversal imaging method for structural health monitoring. Smart Materials and Structures, Vol. 13, Issue 2, 2004, p. 415423. [Publisher]
 Wang D., Ye L., Lu Y., et al. Probability of the presence of damage estimated from an active sensor network in a composite panel of multiple stiffeners. Composites Science and Technology, Vol. 69, Issue 13, 2009, p. 20542063. [Publisher]
 Hay T. R., Royer R. L., Gao H., et al. A comparison of embedded sensor Lamb wave ultrasonic tomography approaches for material loss detection. Smart Materials and Structures, Vol. 15, Issue 4, 2006, p. 946. [Publisher]
 Su Z., Ye L. Identification of damage using lamb waves: from fundamentals to applications. Lecture Notes in Applied and Computational Mechanics, SpringerVerlag London, 2009. [Publisher]
 Crawley E. F., De Luis J. Use of piezoelectric actuators as elements of intelligent structures. AIAA Journal, Vol. 25, 1987, p. 13731385. [Publisher]
 Sohn H., Park H. W., Law K. H., et al. Damage detection in composite plates by using an enhanced time reversal method. Journal of Aerospace Engineering, Vol. 20, Issue 3, 2007, p. 141151. [Publisher]
 Seyedpoor S. M., Montazer M. A damage identification method for truss structures using a flexibilitybased damage probability index and differential evolution algorithm. Inverse Problems in Engineering, Vol. 8, 2015, p. 13031322. [CrossRef]