Next Article in Journal
A Secure LEACH-PRO Protocol Based on Blockchain
Previous Article in Journal
Matching-Based Resource Allocation for Satellite–Ground Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical and Experimental Research on Non-Reference Damage Localization Based on the Improved Two-Arrival-Time Difference Method

1
College of Aeronautical Engineering, Civil Aviation University of China, Tianjin 300300, China
2
College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Sensors 2022, 22(21), 8432; https://doi.org/10.3390/s22218432
Submission received: 23 September 2022 / Revised: 27 October 2022 / Accepted: 27 October 2022 / Published: 2 November 2022
(This article belongs to the Section Physical Sensors)

Abstract

:
The identification of damage based on Lamb waves can hardly avoid obtaining the reference signal under healthy conditions. A non-reference damage localization (NRDL) method is proposed in this paper. The NRDL method is established by the improved two-arrival-time difference method (2/ATDM) and BFGS method. The layout principles of the piezoelectric ceramic lead zirconate titanate (PZT) transducer array in the specimen are proposed. In contrast to existing methods, the damage outside the array in the specimen is identified based on the NRDL. The full-area damage location in the specimen is realized. Furthermore, the optimization of the excitation signal center frequency and transducer array layout is carried out. The damage location accuracy is greatly improved.

1. Introduction

With the rapid development of industry, large plate or plate-shaped metal engineering structures are widely used in the field of aviation, in components such as the fuselage and hatch. In the complex and changeable service environment and working conditions, these structures have sufficient safety and reliability. However, the accumulation of damage and the degradation of function will occur under the long-term influence of human factors and material aging. If the damage is not found in time, it will not only affect the normal use of the structure, but also lead to catastrophic accidents and heavy casualties in extreme cases. Damages cannot be directly identified by visual inspection due to the diversity of forms. With the development of non-destructive testing (NDT) technology (such as eddy current testing and ultrasonic testing), this problem has been alleviated to a certain extent. NDT technology is a form of passive monitoring. It can not immediately obtain the source and severity of damage. With the emergence of intelligent materials, structure health monitoring (SHM) technology has been rapidly developed. Active monitoring makes it possible to effectively estimate the remaining life of the structure. Among these approaches, an active detection technology consisting of Lamb waves based on Piezoelectric ceramic lead zirconate titanate (PZT) transducers is widely used by scholars because of its wide application range and high monitoring accuracy.
Two methods are usually used to analyze the arrival time of Lamb waves: the time domain method and time-frequency domain method. Kehlenbach et al. [1] proposed a method for damage identification based on an elliptic algorithm. It used the time of flight (ToF) to determine the damage localization in the structure. Baptista et al. [2] investigated the performance of time reversal for the localization of impact sources on fiber-reinforced plastic composite structures with embedded piezoelectric sensors. Furthermore, they then studied the feasibility of a piezoelectric diaphragm instead of piezoelectric ceramic and compared the effect of temperature variation on both sensors [3]. The variations in the amplitude, energy and cross-correlation of signals in health and damage states was researched by Burkov et al. [4]. Furthermore, the severity of the damage was obtained. Zhou et al. [5] used the rectangular actuator to excite axisymmetric and non-axisymmetric Lamb waves. They proposed an analytic expression to predict the amplitude of each circular ridge Lamb mode. The sensitivity is improved by selecting the appropriate rectangular actuator for damage detection. Perelli et al. [6] proposed a flat-plate acoustic emission localization strategy with dispersion and reverberation. The method extracted the differences in the distances traversed by stress-guided waves. The difficulties associated with the arrival time detection based on the classical thresholding procedure were overcome. Zhu et al. [7,8] evaluated fatigue damage by using the group velocities of non-linear Lamb waves. An effective mode was selected and a parameter was proposed to quantify the efficiency of the cumulative second harmonic generation (SHG) of Lamb waves. The results show that the S3-S6 modes were sensitive to the evolution of fatigue damage. The integrated amplitude of the second harmonic increased by nearly 300 percent with increasing fatigue period. The non-linear Lamb waves with group velocity mismatch were validated as candidates for effective fatigue damage assessment. Zhang et al. [9] proposed a method for detecting corrosion depth at the edges of holes by using the A0 mode of Lamb waves. Zhan et al. [10] focused on the non-linear response of two different frequencies of Lamb waves in a plate structure with micro-cracks. It was demonstrated that the mixed-frequency Lamb waves technique was sensitive to micro-cracks. The length of the cracks had a relationship with the harmonic amplitude. Rai and Mitra [11,12] proposed a deep learning structure based on a one-dimensional convolutional neural network, which could directly process the recorded original discrete time-domain signals of Lamb waves. The damage detection results were improved by a hybrid physics-assisted multi-layer feed-forward neural network model. Zonzini et al. [13] proposed a multi-layer convolutional neural network and a capsule structure with dedicated time retrieval logic to estimate the time of arrival of acoustic emission signals. It had the particular advantage of providing consistent results in the presence of significantly low signal-to-noise ratios. Ma and Yu [14] proposed quantifying the complex damage of plate structures by using spatiotemporal Lamb waves. The outline information of the damage was clearly obtained. Giurgiutiu et al. [15] applied guided Lamb waves to detect cracks in aluminum specimens. Three wave propagation methods were applied to detect three types of damage. Borate et al. [16] proposed a data-driven structural health monitoring method for the in situ assessment of the localization and extent of damage. Wang et al. [17] used an active sensing method based on Lamb waves and an improved sensor network to detect corrosion expansion at the hole edge. Some scholars have noticed the damage localization method without a baseline. Hameed et al. [18,19] proposed a damage estimation method based on the continuous wavelet transform of the normalized Lamb wave signal. The proposed method can estimate the localization and severity of multiple types of plate damage without baseline signals from undamaged plates. Zhou et al. [20] proposed a reference-free crack-detection method based on Lamb wave mode conversion. The Lamb wave mode conversion technique was adopted and the cracks could be quickly detected. Deng et al. [21] proposed a focused multiple signal classification (MUSIC) algorithm based on the virtual time reversal technique. The baseline free Lamb wave damage localization of isotropic materials was achieved. However, these methods cannot realize damage identification outside the transducer array. De Fenza et al. [22] proposed an astroid algorithm with the two-arrival-time difference method (2/ATDM). In our previous article, a damage shape recognition algorithm was proposed based on the 2/ATDM [23]. The layout and number of sensors have been studied by several scholars. Chen et al. [24] proposed an optimal load-dependent sensor arrangement scheme considering multi-source uncertainty. Furthermore, they [25] proposed an algorithm for structural health monitoring sensor layout based on a subclustering strategy. The algorithm is able to reduce redundancy and improve the performance of sensor configurations.
At present, damage identification based on Lamb waves mainly focuses on damage localization. There are also some studies focused on the imaging of different damage types or multiple damages. None of these efforts can avoid the existence of a problem. It is necessary to obtain the reference Lamb wave signal under healthy conditions. Healthy reference signals are not available in all monitoring processes. Moreover, the layouts of the PZT transducer array in the structure are usually required. Under the same monitoring requirements, the fewer transducers and more concise layout in the array can improve the ability to identify damage. Therefore, a non-reference damage localization (NRDL) method for aluminum alloy plate is presented in this paper. The actual aviation structure is complex, and not all sensors are suitable for deployment. If the monitoring range of the existing transducer array can be expanded, the application of SHM will be improved. Through the further study of the proposed NRDL method, the damage outside the array is also well-monitored in the specimen. The following parts of this article are organized as follows: In Section 2, the proposed NRDL method is presented with the detailed process of combining the improved 2/ATDM and BFGS methods. In Section 3, the basic theory of Lamb waves is introduced. The processes of the experiment are presented and the steps of the simulation are presented in as much detail as possible. The convergence of mesh size and the time increment step are considered. In Section 4, the analysis process of the NRDL method is proposed in detail. Based on NRDL, the numerical results are in good agreement with those of the experiment. Furthermore, the layout of the transducers and the center frequency of the excitation signal are optimized, and better damage localization results are obtained. Finally, some conclusions and limitations are drawn in Section 5.

2. Non-Reference Damage Localization Method

An NRDL method is proposed by improving the governing equation of 2/ATDM. The schematic diagram of the improved 2/ATDM is shown in Figure 1. PZT0 is used as the actuator. PZT1, PZT2 and PZT3 are taken as the sensors. When there is damage in the specimen, the signals received by PZT1, PZT2 and PZT3 contain the direct wave signal from the PZT0 and the reflected signal after the damage. When the damage is close to the boundary, the signal reflected by the boundary will also be received.
Let the time taken by the excitation signal to reach the sensor through damage be t. The distance of the arrived direct wave signal is Z. The distance of the damage reflection signal is D, and the distance of the possible boundary reflection signal is R; then, there is
Z n + D n + R n = c g t n
where n is the sensor number n = 1 , 2 , 3 . c g is the group velocity of the Lamb waves. As shown in Figure 1, the distances from PZT0 to PZT1, PZT2 and PZT3 are the same, and Z 1 = Z 2 = Z 3 . Depending on the layout of the actuator and sensor, the arrival time of the boundary reflection is consistent. Then, there exists R 1 = R 2 = R 3 . Let the received signals of the three groups of sensors differ from each other. The governing equation of the improved 2/ATDM is given by
D 1 D 2 = c g ( t 1 t 2 ) D 2 D 3 = c g ( t 2 t 3 ) D 3 D 1 = c g ( t 3 t 1 )
where D 1 = D ( P 0 d ) + D ( D P 1 ) , D 2 = D ( P 0 d ) + D ( D P 2 ) , D 3 = D ( P 0 d ) + D ( D P 3 ) . D ( P 0 d ) denotes the distance from the actuator to the sensor. D ( D P n ) is the distance from the damage to each sensor. Therefore, the above governing equation can be rewritten as
D ( D P 1 ) D ( D P 2 ) = c g ( t 1 t 2 ) D ( D P 2 ) D ( D P 3 ) = c g ( t 2 t 3 ) D ( D P 3 ) D ( D P 1 ) = c g * ( t 3 t 1 )
According to the Equation (3), any two sets of equations differ from each other. Three sensors can obtain three groups of hyperbolas. These hyperbolas focus on the position of any two sensors. Every two groups of hyperbolas will intersect at one point, and these intersections may be the localization where the damage occurs. Then, the accuracy of t n in the governing equation will seriously affect the final results of localization. The arrival time is easily affected by crosstalk and other factors. Let t n * = t n + t , where t is the peak time of the excitation signal. Equation (3) is rewritten as
D ( D P 1 ) D ( D P 2 ) = c g ( t 1 * t 2 * ) D ( D P 2 ) D ( D P 3 ) = c g ( t 2 * t 3 * ) D ( D P 3 ) D ( D P 1 ) = c g ( t 3 * t 1 * )
The continuous wavelet transform (CWT) is introduced to obtain the accurate t n * [26,27]. The governing equation of CWT is given by
W ( u , s ) = 1 s f ( t ) φ * ( t u s ) d t
where f ( t ) denotes the input signal. s denotes the number of decimal places. u is the position. W ( u , s ) is a soft function of wavelet coefficients f ( t ) , and φ * ( t u s ) is a wavelet function in the complex domain.
Although multiple intersections are found, the localization of damage is still uncertain. A quasi-Newton optimization algorithm, the BFGS method, is introduced. The BFGS method was jointly studied by C.G. Broyden, R. Fletcher, D. Goldfarb, and D.F. Shanno in 1970. It is a method of symmetric positive definite iterative matrix [28]. The BFGS method has excellent numerical stability. It is necessary to introduce the quasi-Newton method before describing the BFGS method. For an objective function f ( x ) existing in the plane, make a vertical line through any point x 1 on the x-axis. The point f ( x 1 ) is obtained, which is also a point on the objective function. Make an intersecting line to f ( x ) passing through f ( x 1 ) . The intersection point x 2 is obtained. The above steps are iterated repeatedly. The obtained f ( x n ) becomes closer to the root of the objective function. However, the computing time also increases. It is necessary to determine a threshold such that when n = k , x k + 1 subtracted from x k is less than this threshold. x k is approximately considered the root of f ( x ) . If f ( x ) is a differentiable function, then there is
x k + 1 = x k f ( x k ) f ( x k )
Based on the second-order Taylor expansion, the Equation (6) is rewritten as
x k + 1 = x k f ( x k ) f ( x k )
The optimized objective function f ( x ) is a multivariate function. The Hessian matrix is introduced and the equivalent matrix of Equation (7) is given by
x k + 1 = x k H k 1 g k , k = 0 , 1
where g k is f ( x ) in Equation (7) and H k is the second derivative f ( x ) . The BFGS method is an algorithm to find H k 1 by iteration. The governing equation of the BFGS method is given by
D k + 1 = ( I s k y k T y K T s k ) D k ( I y k s k T y k T s k ) + s k s k T y k T s k
where D k is H k 1 , s k = x k + 1 x k , y k = g k + 1 g k , and g k is the derivative function of the original function. When k = 1 , g k is the unit matrix.
The objective function of the NRDL method is constructed. j intersections can be obtained through the improved 2/ATDM. The intersection is considered as valid in the specimen. The invalid intersection can be removed artificially. ( x i , y i ) ( i = 1 , 2 , 3 j ) is the i-th coordinate. ( x * , y * ) is the coordinate of the optimal point. The sum of the distance P from the optimal point to each intersection can be expressed as
P = i = 1 j ( x i x * ) 2 + ( y i y * ) 2
The objective function of the NRDL method can be defined as
F ( x , y ) = m i n { q } = m i n { i = 1 j ( x i x * ) 2 + ( y i y * ) 2 }

3. Experiment and Simulation

3.1. Lamb Waves

In elastic solids, Lamb waves are excited in the form of elastic waves. The dispersion equation of Lamb waves is derived from elasticity [23].
t a n ( β d ) t a n ( α d ) = 4 ε 2 α β ( ε 2 β 2 ) 2 t a n ( β d ) t a n ( α d ) = ( ε 2 β 2 ) 2 4 ε 2 α β
where d denotes the thickness of the plate. In the variables α = ω 2 c l 2 k 2 and β = ω 2 c t 2 k 2 , ε denotes the wave number, and ω is the angular frequency of the resonant wave. E, υ , and ρ denote the Young’s modulus, Poisson’s ratio, and density of the plate, respectively. c l = E ( 1 ν ) ρ ( 1 + ν ) ( 1 2 ν ) and c t = E 2 ρ ( 1 + ν ) are the propagation velocity of the plane stress longitudinal wave.
As can be seen from Equation (12), the Lamb wave has the characteristics of multiple modes and dispersion. When the thickness of the plate is 1 mm, the phase velocity dispersion curves of Lamb wave are shown in Figure 2a. The relationship between group velocity and phase velocity is given by
c g = c p [ 1 ω c p c p ω ] 1
The group velocity dispersion curves are shown in Figure 2b. When the frequency is less than 100 kHz, only the S0 and A0 mode are generated. Although the wavelength of the A0 mode is larger than that of the S0 mode, the propagation speed of the S0 mode is greater, and can reduce the superposition of boundary reflection signals and damage signals. Therefore, the S0 mode is adopted.

3.2. Details of Experiment

To validate the NRDL method, an experiment with the S0 mode is conducted on a 600 mm × 600 mm × 1 mm aluminum alloy plate. The material properties of the plate specimen are given in Table 1. Lamb waves are excited by a Keysight 33500b function generator. A Tektronix 3024 four-channel digital storage oscilloscope is used to receive the response signals. In experiment, the plate specimen is placed on a soft towel to simulate the free vibration. The experimental setup is shown in Figure 3.
The narrow-band excitation signal is a five-cycle sine burst centered at 350 kHz, as shown in Figure 4. The amplitude of the excitation signal is 10 V. Five PZT transducers (PZT0, PZT1, PZT2, PZT3 and PZT4) are arranged on the specimen according to the NRDL method.
Taking the lower left corner of the specimen as the coordinate origin, the coordinates of those transducers are (300,300), (200,200), (400,200), (400, 400) and (200, 400), respectively. These PZT transducers are selected by the D33 instrument and WK6500B low-frequency impedance analyzer, and their properties are almost the same. The diameter of these PZT transducers is 10 mm and the thickness is 0.5 mm. The properties are shown in Table 1. A 6 mm diameter hole damage is prefabricated. Generally, the PZT transducer array divides the specimen into three regions: A, B and C. The three regions are represented by green, yellow and blue, respectively. A schematic diagram is shown in Figure 5.
The prefabricated hole damage centers in the three regions are D A ( 220 , 360 ) , D B ( 320 , 440 ) and D C ( 460 , 440 ) , respectively. In order to distinguish the damages in the three regions, these damages are called Case I, Case II and Case III, as shown in Figure 6.

3.3. Numerical Simulation

The NRDL method is numerically studied in this paper. Only the finite element analysis model of Case I is given, as shown in Figure 7. The in-plane mesh size of the three-dimensional model is 1 mm × 1 mm and that in the thickness direction is 0.5 mm. The hole damage and the surroundings of the piezoelectric sheet are encrypted. The SOLID5, SOLID185 and SOLID45 elements are employed to simulate the plate, PZT transducer and bonding material, respectively. The total numbers of nodes and elements of the proposed model are 715,488 and 1,082,965, respectively.
In our previous research, the requirements for the convergence of mesh size and time increment steps were given [23]. The mesh size is affected by the group velocity of the S0 mode. The time increment step is subject to the center frequency of the excitation signal. The convergence analysis results are shown in Figure 8.
Finally, as shown in Figure 8, the mesh size is 1 mm, and the time increment step is 1.5 × 10 7 s.

4. Results and Discussion

4.1. Validations

The responses of the Lamb waves in Case I, Case II and Case III are experimentally and numerically compared. The results are shown in Figure 9, Figure 10 and Figure 11, where (a–d) indicate the excitation signals received by PZT1, PZT2, PZT3 and PZT4 as sensors, respectively.
As shown in Figure 9, Figure 10 and Figure 11, the results of the simulation are in good agreement with those of the experiment in each case. Taking Case I as an example, the difference between the six groups of received signals is calculated. The six groups of difference signals are processed by CWT. The CWT coefficients of the numerical results are shown in Figure 12. In order to further effectively extract the information of the difference signal, three principles are followed,
(1) The same properties of the transducers make the direct arrival signals of S0 mode and A0 mode the same. The actuator and sensors are symmetrically laid in the specimen. The signals reflected from the boundary are almost the same. When damage occurs, the only difference between the signals received by each sensor is the reflected signal of the damage. In order to avoid the superimposed signal interfering with the extraction of the damage signal, the S0 mode with faster propagation speed is adopted as the signal for identifying the damage.
(2) Taking Case I (damage occurs in the PZT transducers array) as an example, there are two extremes in the localization of damage: (a) the damage completely overlaps with the position of the actuator PZT0; and (b) the damage overlaps with the position of any sensor (PZT1, PZT2, PZT3 and PZT4). For these two extremes, there are also two possibilities for the signal reflected by the damage: (a) when the damage completely overlaps with the position of the actuator PZT0, the signals received by all sensors are consistent. The arrival time of the excitation signal is 2.6 × 10 5 s. The t i * in Equation (4) is 3.33 × 10 5 s; (b) when the damage overlaps with any sensor position, the longest arrival time of the excitation signal is 7.83 × 10 5 s. The t i * in Equation (4) is 8.55 × 10 5 s. Therefore, the CWT coefficient of the difference signal is valid only when 3.33 × 10 5 s < t i * < 8.55 × 10 5 s.
(3) The following conditions shall be met in order to prevent the differential signal from being interfered with:
D ( D P i ) D ( D P j ) c g T ( i j )
where the subscripts i and j denote the numbers of the two different piezoelectric sensors; c g denotes the group velocity of S0; and D ( D P i ) and D ( D P j ) are the distance between the damage and the different sensor, respectively. T denotes one period of the excitation signal. This can effectively ensure that the signals reflected by the damage are not superimposed. There is no interference between the acquired difference signals.
According to the principles of the proposed NRDL method, the times of the two different signal peaks reflected by the damage can be obtained from Figure 12a,c,e,f, respectively. Each set of time is brought into Equation (4). Multiple groups of hyperbolas are calculated and drawn. The number of sensors corresponding to the time of the damage signal can be clearly marked in Figure 12, because the hole damage is prefabricated. However, the number of sensors in actual monitoring cannot be determined. The proposed method requires all signals to be different from one another. The calculation result of Equation (4) is not affected. It is only necessary to calibrate the serial numbers of the sensor in a certain order. As shown in Figure 12b, two peaks are not found in the effective time range. This group of signals is discarded. As shown in Figure 12d, although the time of two peaks can be found within the set time range, the times are too close. The signal superposition may occur and the real time cannot be determined. Therefore, this group of signals is also discarded according to principle (3). Multiple groups of hyperbolas can be obtained on the basis of complying with the above principles, as shown in Figure 13.
As shown in Figure 13, seven intersections can be obtained. The coordinates of the intersections are shown in Table 2.
Similarly, the results of the experiment are also processed. The CWT coefficients are shown in Figure 14.
By extracting the time domain information of the effective CWT coefficient, the hyperbolas of the experimental results are obtained, as shown in Figure 15. Seven intersections are also obtained. The coordinates of the intersections are shown in Table 3.
The results in Table 2 and Table 3 are brought into the objective function F ( x , y ) . The coordinates of the predicted damage localization center are shown in Table 4. The coordinates of the actual damage center are (220, 360).
As shown in Table 4, the localization of damage is accurately obtained in experiments and numerical simulations. This shows the effectiveness of the proposed method. However, there are still some errors in the experimental results from the perspective of absolute distance, for several reasons: (a) there may be some errors in the properties of PZT transducers in the experiment; (b) the existence of crosstalk makes some small disturbances unable to be eliminated by subtraction; and (c) the boundary of the prefabricated hole damage is not smooth. Usually, Lamb waves are reflected at the edge of the damage. A rough boundary will cause error, but this error may be positive. According to the above process, the damages in Case II and Case III are predicted. The coordinates of the damage center in Case II and Case III are (320, 440) and (460, 440), respectively. In these two cases, the predictions of damage localization are shown in Table 5.
As shown in Table 4 and Table 5, there is a large error in damage localization in Case III. The excitation signal takes a long time to reach the sensors, and is easily superimposed with the reflected signal of the boundary. The proposed method needs to be optimized. The intersections of curves determined by the NRDL method are the possible locations of the damage. If more curves can be obtained, more possible positions of the damage can be obtained. More accurate damage prediction can be obtained by introducing the BFGS method. According to the principle of the proposed method, the relationship between the difference signal and the period of the excitation signal needs to be considered. Two optimization methods are proposed: (a) frequency optimization, that is, the reason for choosing 350 kHz as the center frequency; and (b) the layout of eight PZT transducers is adopted.

4.2. Frequency Optimization

Taking Case I as an example, four groups of sinusoidal modulated signals with center frequencies of 200 kHz, 250 kHz, 300 kHz and 350 kHz are selected as excitation signals, respectively. According to the four center frequencies, the corresponding T in Equation (14) is calculated as 2.5 × 10 5 s, 2 × 10 5 s, 1.67 × 10 5 s and 1.428 × 10 5 s respectively. The numerical results of received signals with different center frequencies (200 kHz, 250 kHz and 300 kHz) are compared with the experimental results, as shown in Figure 16, Figure 17 and Figure 18.
As shown in Figure 16, Figure 17 and Figure 18, the results of the simulations are in good agreement with those of the experiment. The signal processed by CWT is compared with the actual damage localization. When the center frequency is 200 kHz, only two groups of difference signals (PZT2-PZT4 and PZT3-PZT4) can be obtained. When the center frequency is 250 kHz and 300 kHz, three groups of difference signals (PZT1-PZT4, PZT2-PZT4 and PZT3-PZT4) can be obtained. Based on the above time-domain information, the predictions of damage center coordinates are obtained, as shown in Table 6. The absolute distance between the predicted damage center and the actual damage center is also given, as shown in Table 7.
According to Table 6 and Table 7, when the center frequency of the excitation signal is 200 kHz, only two curves can be obtained. Only one intersection cannot be optimized by the BFGS algorithm. The intersection of these two curves will be considered as the predicted damage center. This obviously shows that the error is large in comparison with the results under other central frequencies. In particular, the damage location of the experiment is far beyond the actual damage area. The absolute distance is 17.5 times that of the hole radius. When the center frequency of the excitation signal is 250 kHz and 300 kHz, three hyperbolas can be obtained. Three hyperbolas intersect at one point in pairs. The prediction of damage localization can be obtained by the BFGS method. Since the number of intersections is the same, the two groups of results have the same number of differential signals. The prediction error is basically the same. The absolute distances in the two groups of experiments are about 4.5 times that of the hole radius. Therefore, the excitation signal with a central frequency of 350 kHz is finally adopted. Four groups of differential signals can be obtained. The prediction of the damage localization is more accurate. A large center frequency is not necessarily better. When the center frequency increases, the sampling points used for the excitation function need to be increased significantly. The received signal is prone to distortion. A valid ToF may not be obtained, which puts forward high requirements for the sampling frequency of the experiment. The amount of calculation increases exponentially for the simulation. Therefore, the most reasonable center frequency shall be determined.

4.3. Layout Optimization

According to the governing equation of the proposed NRDL method, the subtraction between the received signals can filter out the direct waves of the excitation signal. PZT transducers with the same distance from the actuator in the specimen are arranged. More response results for differential signals can be obtained. Therefore, based on the previous layout, a circle is drawn with the original actuator as the center and the distance between the actuator and the sensor as the radius. New PZT transducers are added in the vertical and horizontal directions on the circle, respectively. Taking Case I as an example, the new layout of the sensor is shown in Figure 19a.
The transmission of Lamb waves is affected by the form of meshes in the numerical model. The mesh of the central radiation is adopted to ensure that it is the same from the actuator to each sensor, as shown in Figure 19b. The numbers of elements and nodes of the new model are 630,650 and 958,955, respectively. The calculation scale is reduced compared with the previous model. The number of sensors is increased to eight. The original red PZT transducers (PZT1, PZT3, PZT5 and PZT7) are taken as one group, and the new blue PZT transducers (PZT2, PZT4, PZT6 and PZT8) are taken as another group. A new set of hyperbolas can be drawn by another four sensors. The hyperbolas of the experiment and simulation are shown in Figure 20.
As shown in Figure 20, a total of 32 intersections are obtained from the experimental hyperbola and 34 intersections are obtained from the numerical hyperbola. The coordinate values of the corresponding intersections are shown in Table 8 and Table 9, respectively. These intersections can be used to obtain a new prediction of damage localization, as shown in Table 10. After the optimization of the layout, the absolute distance between the center of the predicted damage localization and the actual center is significantly reduced.
The number of sensors was increased to eight. More hyperbolas can be obtained by using the received signals of all sensors. The prediction of damage localization through the intersection obtained by eight sensors is shown in Table 11.
The prediction of a group of eight sensors is more accurate than that of a group of four sensors. However, it is no better than the prediction of two groups of four sensors. The global error increases due to some intersections with large errors. The amount of calculation using the BFGS method also increases. Hence, the layout optimization of two groups of four sensors is used to predict the damage in Case II and Case III. The optimized prediction is shown in Table 12.
Compared with Table 4 and Table 5, the identification of two groups of four sensors has been greatly improved. In Case I, it is obvious that the absolute distances are greatly reduced in the experiment and simulation, and are far less than the radius of hole damage. In particular, the experimental results show that more intersections are obtained near the actual center of the damage. The local error caused by individual points is reduced. This also occurs in the numerical results of Case III. Similarly, the experimental results of Case III are optimized to some extent. However, the influence of the experimental environment and the insufficient number of intersections lead to the limitation of optimization, which also puts forward new challenges for further optimization.

5. Conclusions

A numerical simulation model of an aluminum alloy plate specimen with hole damage is proposed based on the Lamb waves technique, which can be used to determine the localization of the damage without healthy reference signals. The NRDL method is proposed based on the improved 2/ATDM and BFGS methods. The numerical results of the NRDL method are verified by those of the experiment. Two optimization models based on frequency improvement and layout improvement are established. From the conducted numerical and experimental investigations, the following main conclusions emerge:
  • The damage location without a reference is realized based on the proposed method. The error of the damage location is effectively reduced by the optimization of the center frequency and array layout. This is of great significance for the application of practical engineering.
  • The full-area damage localization is realized in the specimen. The damage identification range of the PZT transducer array is improved. The accuracy of damage identification is higher in Region B outside the array than that in Region C. Furthermore, damage shape identification should use more points in the array and Region B.
  • The principles for extracting good time-frequency domain information are established. These principles can be used to effectively obtain the characteristics of the response signal. The influence of the signal superposition is reduced.
Some limitations of the proposed method are as follows:
  • Small differences in the properties of the PZT transducers will lead to errors in the results during the experiment, which puts forward higher requirements for the manufacture and selection of PZT transducers.
  • The proposed method requires that the distance between each sensor and the actuator is consistent. The accuracy of mesh discretization is high in numerical simulation.
  • If the prefabricated damage is not smooth during the experiment, the damage reflection signal and A0 signal may be superimposed, so the damage signal cannot be effectively separated.

Author Contributions

Conceptualization, Y.W.; Data curation, K.L.; Formal analysis, D.L.; Funding acquisition, X.S.; Methodology, P.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Science Foundations of China number 12072364, Natural Science Foundations of Tianjin (21JCQNJC00880, 19JCQNJC02500).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study is available on request from the corresponding author.

Acknowledgments

The authors would like to thank to all anonymous reviewers and editors for their helpful suggestions for the improvement of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kehlenbach, M.; Hanselka, H. Automated structural integrity monitoring based on broadband lamb wave excitation and matched filtering. In Proceedings of the 44th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference, Norfolk, VA, USA, 7–10 April 2003; p. 1563. [Google Scholar]
  2. Coles, A.; Castro, B.; Andreades, C.; Baptista, F.G.; Ciampa, F. Impact Localization in Composites Using Time Reversal, Embedded PZT Transducers, and Topological Algorithms. Front. Built Environ. 2020, 6, 27. [Google Scholar] [CrossRef] [Green Version]
  3. Campeiro, L.M.; Budoya, D.E.; Baptista, F.G. Lamb wave inspection using piezoelectric diaphragms: An initial feasibility study. Sens. Actuators A Phys. 2021, 331, 112859. [Google Scholar] [CrossRef]
  4. Burkov, M.; Eremin, A.; Byakov, A.; Lyubutin, P.; Panin, S. Impact Damage Detection in Laminate and Honeycomb CFRPs using Lamb Wave Ultrasonic Sensing. Russ. J. Nondestruct. Test. 2021, 57, 114–124. [Google Scholar] [CrossRef]
  5. Zhou, K.; Zhou, Z.; Yang, Z.; Xu, X.; Wu, Z. Axisymmetric and non-axisymmetric Lamb wave excitation using rectangular actuators. Smart Mater. Struct. 2019, 28, 115024. [Google Scholar] [CrossRef]
  6. Perelli, A.; De Marchi, L.; Marzani, A.; Speciale, N. Acoustic emission localization in plates with dispersion and reverberations using sparse PZT sensors in passive mode. Smart Mater. Struct. 2012, 21, 025010. [Google Scholar] [CrossRef]
  7. Zhu, W.; Xiang, Y.; Liu, C.J.; Deng, M.; Xuan, F.Z. A feasibility study on fatigue damage evaluation using nonlinear Lamb waves with group-velocity mismatching. Ultrasonics 2018, 90, 18–22. [Google Scholar] [CrossRef]
  8. Gao, G.; Liu, C.; Hu, N.; Deng, M.; Chen, H.; Xiang, Y. Response of second-harmonic generation of Lamb wave propagation to microdamage thickness in a solid plate. Wave Motion 2020, 96, 102557. [Google Scholar] [CrossRef]
  9. Zhang, W.; Xu, D.; Wang, X. Corrosion Depth Monitoring of Hole-Edge Based on Lamb wave. J. Phys. Conf. Ser. 2021, 1930, 012014. [Google Scholar] [CrossRef]
  10. Zhan, Z.; Wang, F.; Zheng, Y.; Wang, Z.; Huang, S.; Li, X.; Wang, S. Nonlinear Responds of Lamb Waves in Plate Structure with Micro-Crack Using Frequency-Mixing Technique. J. Phys. Conf. Ser. 2021, 1877, 012014. [Google Scholar] [CrossRef]
  11. Rai, A.; Mitra, M. A hybrid physics-assisted machine-learning-based damage detection using Lamb wave. Sādhanā 2021, 46, 1–11. [Google Scholar] [CrossRef]
  12. Rai, A.; Mitra, M. Lamb wave based damage detection in metallic plates using multi-headed 1-dimensional convolutional neural network. Smart Mater. Struct. 2021, 30, 035010. [Google Scholar] [CrossRef]
  13. Zonzini, F.; Bogomolov, D.; Dhamija, T.; Testoni, N.; De Marchi, L.; Marzani, A. Deep Learning Approaches for Robust Time of Arrival Estimation in Acoustic Emission Monitoring. Sensors 2022, 22, 1091. [Google Scholar] [CrossRef]
  14. Ma, Z.; Yu, L. Lamb wave imaging with actuator network for damage quantification in aluminum plate structures. J. Intell. Mater. Syst. Struct. 2021, 32, 182–195. [Google Scholar] [CrossRef]
  15. Giurgiutiu, V.; Cuc, A. Embedded non-destructive evaluation for structural health monitoring, damage detection, and failure prevention. Shock Vib. Dig. 2005, 37, 83. [Google Scholar] [CrossRef]
  16. Borate, P.; Wang, G.; Wang, Y. Data-driven structural health monitoring approach using guided Lamb wave responses. J. Aerosp. Eng. 2020, 33, 04020033. [Google Scholar] [CrossRef]
  17. Wang, X.; Dai, W.; Xu, D.; Zhang, W.; Ran, Y.; Wang, R. Hole-edge corrosion expansion monitoring based on lamb wave. Metals 2020, 10, 1469. [Google Scholar] [CrossRef]
  18. Hameed, M.S.; Li, Z.; Zheng, K. Damage Detection Method Based on Continuous Wavelet Transformation of Lamb Wave Signals. Appl. Sci. 2020, 10, 8610. [Google Scholar] [CrossRef]
  19. Chen, J.; Li, Z.; Gong, K. Nondestructive testing method based on lamb waves for localization and extent of damage. Acta Mech. Solida Sin. 2017, 30, 65–74. [Google Scholar] [CrossRef]
  20. Zhou, L.; He, Z.; Sun, H. Lamb wave mode conversion-based crack detection for plate-like structures without baseline information. J. Vibroeng. 2013, 15, 647–657. [Google Scholar]
  21. Xu, C.; Wang, J.; Yin, S.; Deng, M. A focusing MUSIC algorithm for baseline-free Lamb wave damage localization. Mech. Syst. Signal Process. 2022, 164, 108242. [Google Scholar] [CrossRef]
  22. De Fenza, A.; Sorrentino, A.; Vitiello, P. Application of Artificial Neural Networks and Probability Ellipse methods for damage detection using Lamb waves. Compos. Struct. 2015, 133, 390–403. [Google Scholar] [CrossRef]
  23. Wu, Y.; Shen, X.; Li, D. Numerical and experimental research on damage shape recognition of aluminum alloy plate based on Lamb wave. J. Intell. Mater. Syst. Struct. 2021, 32, 2273–2288. [Google Scholar] [CrossRef]
  24. Yang, C.; Ouyang, H. A novel load-dependent sensor placement method for model updating based on time-dependent reliability optimization considering multi-source uncertainties. Mech. Syst. Signal Process. 2022, 165, 108386. [Google Scholar] [CrossRef]
  25. Yang, C.; Liang, K.; Zhang, X.; Geng, X. Sensor placement algorithm for structural health monitoring with redundancy elimination model based on sub-clustering strategy. Mech. Syst. Signal Process. 2019, 124, 369–387. [Google Scholar] [CrossRef]
  26. Rioul, O.; Vetterli, M. Wavelets and signal processing. IEEE Signal Process. Mag. 1991, 8, 14–38. [Google Scholar] [CrossRef] [Green Version]
  27. Sinha, S.; Routh, P.S.; Anno, P.D.; Castagna, J.P. Spectral decomposition of seismic data with continuous-wavelet transform. Geophysics 2005, 70, 19–25. [Google Scholar] [CrossRef]
  28. Li, D.H.; Fukushima, M. A modified BFGS method and its global convergence in nonconvex minimization. J. Comput. Appl. Math. 2001, 129, 15–35. [Google Scholar] [CrossRef]
Figure 1. The schematic of the improved 2/ATDM.
Figure 1. The schematic of the improved 2/ATDM.
Sensors 22 08432 g001
Figure 2. Dispersion curve of aluminum alloy plate.
Figure 2. Dispersion curve of aluminum alloy plate.
Sensors 22 08432 g002
Figure 3. Experimental setup.
Figure 3. Experimental setup.
Sensors 22 08432 g003
Figure 4. Narrow-band excitation signal with a central frequency of 350 kHz.
Figure 4. Narrow-band excitation signal with a central frequency of 350 kHz.
Sensors 22 08432 g004
Figure 5. The schematic of specimen area division.
Figure 5. The schematic of specimen area division.
Sensors 22 08432 g005
Figure 6. (a) Details of specimen, and (b) schematic.
Figure 6. (a) Details of specimen, and (b) schematic.
Sensors 22 08432 g006
Figure 7. Details of numerical model.
Figure 7. Details of numerical model.
Sensors 22 08432 g007
Figure 8. Results of convergence analysis: (a) mesh size and (b) time increment step.
Figure 8. Results of convergence analysis: (a) mesh size and (b) time increment step.
Sensors 22 08432 g008
Figure 9. Comparison of experimental and simulation results under Case I.
Figure 9. Comparison of experimental and simulation results under Case I.
Sensors 22 08432 g009
Figure 10. Comparison of experimental and simulation results under Case II.
Figure 10. Comparison of experimental and simulation results under Case II.
Sensors 22 08432 g010
Figure 11. Comparison of experimental and simulation results under Case III.
Figure 11. Comparison of experimental and simulation results under Case III.
Sensors 22 08432 g011
Figure 12. CWT coefficient of the difference signal of any two sensors under the simulation of Case I.
Figure 12. CWT coefficient of the difference signal of any two sensors under the simulation of Case I.
Sensors 22 08432 g012
Figure 13. Intersection of simulations under Case I.
Figure 13. Intersection of simulations under Case I.
Sensors 22 08432 g013
Figure 14. CWT coefficient of the difference signal of any two sensors in the experiment under Case I.
Figure 14. CWT coefficient of the difference signal of any two sensors in the experiment under Case I.
Sensors 22 08432 g014
Figure 15. Intersection of experiment for Case I.
Figure 15. Intersection of experiment for Case I.
Sensors 22 08432 g015
Figure 16. Comparison of experimental and simulation results under the center frequency of 200 kHz.
Figure 16. Comparison of experimental and simulation results under the center frequency of 200 kHz.
Sensors 22 08432 g016
Figure 17. Comparison of experimental and simulation results under the center frequency of 250 kHz.
Figure 17. Comparison of experimental and simulation results under the center frequency of 250 kHz.
Sensors 22 08432 g017
Figure 18. Comparison of experimental and simulation results under the center frequency of 300 kHz.
Figure 18. Comparison of experimental and simulation results under the center frequency of 300 kHz.
Sensors 22 08432 g018
Figure 19. (a) The schematic of layout optimization, and (b) the numerical model.
Figure 19. (a) The schematic of layout optimization, and (b) the numerical model.
Sensors 22 08432 g019
Figure 20. Hyperbolas of two groups of 4 sensors under Case I: (a) experiment and (b) simulation.
Figure 20. Hyperbolas of two groups of 4 sensors under Case I: (a) experiment and (b) simulation.
Sensors 22 08432 g020
Table 1. Material properties of plate specimen and PZT transducers.
Table 1. Material properties of plate specimen and PZT transducers.
Properties of the Plate Specimen
Young’s modulus ( GPa ) 67
Poisson’s ratio0.34
Mass density ( kg / m 3 ) 2730
Properties of the PZT transducers
Young’s modulus ( GPa ) 12.8
Piezoelectric constant ( c / m 2 ) −5.3
Relative permittivity ( ε 33 ) 1000
Mass density ( kg / m 3 ) 7500
Table 2. Coordinates of the intersections under simulation.
Table 2. Coordinates of the intersections under simulation.
Point 1Point 2Point 3Point 4Point 5Point 6Point 7
x 219.6936219.7961221.7283220.275224.5164226.1663126.9352
y 357.5862357.3172352.2234357.6784358.427362.453581.902
Table 3. Coordinates of the intersections under experiment.
Table 3. Coordinates of the intersections under experiment.
Point 1Point 2Point 3Point 4Point 5Point 6Point 7
x 213.2933221.7281222.1135222.6515225.2257241.8125168.6522
y 350.3837369.0463367.9631369.2788359.1532375.9530509.9387
Table 4. Comparison of experiment and simulation results under Case I.
Table 4. Comparison of experiment and simulation results under Case I.
Experimental ResultsNumerical Results
CoordinateErrorAbsolute DistanceCoordinateErrorAbsolute Distance
Case I x 222.24261.02%8.8625 x 220.27510.13%2.3377
y 368.57412.38% y 357.67850.64%
Table 5. Comparison of experiment and simulation results under Case II and Case III.
Table 5. Comparison of experiment and simulation results under Case II and Case III.
Experimental ResultsNumerical Results
CoordinateErrorAbsolute DistanceCoordinateErrorAbsolute Distance
Case II x 318.90590.256%2.2472 x 319.39250.189%3.5082
y 438.03710.117% y 436.54470.785%
Case III x 449.16472.35%10.8504 x 447.45722.73%12.8374
y 439.42730.130% y 437.26560.621%
Table 6. Comparison of experimental and simulation results under different center frequencies.
Table 6. Comparison of experimental and simulation results under different center frequencies.
Center FrequencyExperimental ResultsNumerical Results
CoordinateErrorCoordinateError
200 kHz x 159.730027.40% x 223.39141.54%
y 323.401610.17% y 361.60550.45%
250 kHz x 235.68267.13% x 222.35101.07%
y 371.27153.13% y 359.04240.27%
300 kHz x 218.86260.52% x 222.28831.04%
y 341.88125.03% y 358.67380.37%
Table 7. Absolute distance of experiment and simulation under different center frequencies.
Table 7. Absolute distance of experiment and simulation under different center frequencies.
200 kHz250 kHz300 kHz350 kHz
Numerical results3.75222.53852.64482.3377
Experimental results70.511819.313018.15458.8625
Table 8. Coordinates of the experimental intersections of two groups of 4 sensors under Case I.
Table 8. Coordinates of the experimental intersections of two groups of 4 sensors under Case I.
PointCoordinatesPointCoordinates
x(mm)y(mm)x(mm)y(mm)
1221.4553364.546317223.3166364.9534
2224.0102357.471518222.7508357.7103
3223.2137359.685519222.6048357.7382
4223.6227358.549520224.0965357.4553
5224.0238357.433821223.1458357.6351
6221.2460365.122622220.2911358.1851
7181.7809468.554023222.9364358.4981
8222.2292362.411324221.1110356.5792
9156.6495531.742125223.0379358.0812
10201.0410362.255026218.5730354.6916
11140.8030378.051227223.2491357.2055
12224.4695365.221228219.2470356.0378
13233.8489367.822329222.0560362.0188
14350.6063434.999930143.4028547.3921
15221.4017364.535131297.7728294.1370
16230.9606398.105032297.6905294.3483
Table 9. Coordinates of the simulation intersections of two groups of 4 sensors under Case I.
Table 9. Coordinates of the simulation intersections of two groups of 4 sensors under Case I.
PointCoordinatesPointCoordinates
x(mm)y(mm)x(mm)y(mm)
1220.5253360.302218229.7778357.0491
2221.1313358.676519227.2649357.5072
3210.3383387.105820218.4681359.2059
4220.3237360.842321206.4816351.6225
5222.5655354.812122206.2354349.9729
6232.3564327.632723211.8637410.1102
7219.8613362.078924200.2313348.2598
8152.9875529.18832593.8966321.1719
9216.3509359.635826225.4863365.1087
10105.3407389.832727218.6479359.5679
11207.4560358.706728227.5217356.3419
12219.3574360.101329214.9374352.5160
13337.1930414.848130223.8542371.5136
14226.3602361.468631171.8906507.0888
15218.8712360.021032213.5070450.0110
16207.8005361.449933203.3735333.9378
17218.2258359.254734143.9420260.5652
Table 10. Experimental and simulation results of two groups of 4 sensors under Case I.
Table 10. Experimental and simulation results of two groups of 4 sensors under Case I.
Experimental ResultsNumerical Results
CoordinateErrorAbsolute DistanceCoordinateErrorAbsolute Distance
Case I x 218.88060.51%1.1768 x 222.38151.080%2.7567
y 360.36310.10% y 358.61150.385%
Table 11. Experimental and simulation results of 8 sensors under Case I.
Table 11. Experimental and simulation results of 8 sensors under Case I.
Experimental ResultsNumerical Results
CoordinateErrorAbsolute DistanceCoordinateErrorAbsolute Distance
Case I x 218.90520.50%1.1648 x 222.88931.31%2.5595
y 359.60220.111% y 358.42270.438%
Table 12. Experimental and simulation results of two groups of 4 sensors under Case II and Case III.
Table 12. Experimental and simulation results of two groups of 4 sensors under Case II and Case III.
Experiment ResultsNumerical Results
CoordinateErrorAbsolute DistanceCoordinateErrorAbsolute Distance
Case II x 320.36610.475%2.1066 x 319.00080.3123%2.9088
y 437.92540.080% y 437.26820.621%
Case III x 467.06551.53%7.7949 x 460.44870.097%2.9586
y 443.29240.748% y 437.07560.664%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wu, Y.; Liu, K.; Li, D.; Shen, X.; Lu, P. Numerical and Experimental Research on Non-Reference Damage Localization Based on the Improved Two-Arrival-Time Difference Method. Sensors 2022, 22, 8432. https://doi.org/10.3390/s22218432

AMA Style

Wu Y, Liu K, Li D, Shen X, Lu P. Numerical and Experimental Research on Non-Reference Damage Localization Based on the Improved Two-Arrival-Time Difference Method. Sensors. 2022; 22(21):8432. https://doi.org/10.3390/s22218432

Chicago/Turabian Style

Wu, Yaogang, Kangwei Liu, Dinghe Li, Xing Shen, and Pengcheng Lu. 2022. "Numerical and Experimental Research on Non-Reference Damage Localization Based on the Improved Two-Arrival-Time Difference Method" Sensors 22, no. 21: 8432. https://doi.org/10.3390/s22218432

APA Style

Wu, Y., Liu, K., Li, D., Shen, X., & Lu, P. (2022). Numerical and Experimental Research on Non-Reference Damage Localization Based on the Improved Two-Arrival-Time Difference Method. Sensors, 22(21), 8432. https://doi.org/10.3390/s22218432

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