Lamb Wave Local Wavenumber Approach for Characterizing Flat Bottom Defects in an Isotropic Thin Plate

This paper aims to use the Lamb wave local wavenumber approach to characterize flat bottom defects (including circular flat bottom holes and a rectangular groove) in an isotropic thin plate. An air-coupled transducer (ACT) with a special incidence angle is used to actuate the fundamental anti-symmetric mode (A0). A laser Doppler vibrometer (LDV) is employed to measure the out-of-plane velocity over a target area. These signals are processed by the wavenumber domain filtering technique in order to remove any modes other than the A0 mode. The filtered signals are transformed back into the time-space domain. The space-frequency-wavenumber spectrum is then obtained by using three-dimensional fast Fourier transform (3D FFT) and a short space transform, which can retain the spatial information and reduce the magnitude of side lobes in the wavenumber domain. The average wavenumber is calculated, as a real signal usually contains a certain bandwidth instead of the singular frequency component. Both simulation results and experimental results demonstrate that the average wavenumber can be used not only to identify shape, location, and size of the damage, but also quantify the depth of the damage. In addition, the direction of an inclined rectangular groove is obtained by calculating the image moments under grayscale. This hybrid and non-contact system based on the local wavenumber approach can be provided with a high resolution.


Introduction
Lamb waves have shown great potential for structural health monitoring (SHM) in plate-like structures [1]. Their attractive features include sensitivity to a variety of damage types and the capability of traveling relatively long distances [2]. Damage imaging methods based on Lamb waves have been widely studied by many researchers, such as delay-and-sum (DAS) imaging [3,4], time reversal focusing imaging [5,6], diffraction tomography imaging [7], and ultrasonic phased array [8].
Lamb waves are dispersive and multimodal [9]. Moreover, the propagating Lamb waves may include incident, reflected, and converted waves when they encounter a sudden thickness variation [10,11], such as circular flat bottom holes and rectangular grooves. Obviously, various wave modes make the interpretation of Lamb waves very difficult [12].
In previous literature, the mode identification and separation of Lamb waves have been widely investigated in the time domain, time-frequency domain, and frequency-wavenumber domain [13,14]. Short-time Fourier transform (STFT) with a sliding window can make the variant distribution of the energy spectrum act as a function of time, which is capable of identifying Lamb wave mode [15]. on Snell's law. Then, the fundamental anti-symmetric mode (A0) can be actuated by an ACT with this special incidence angle. A scanning grid is predefined before using a non-contact LDV to measure the velocities of the desired mode. Obviously, the velocity U(x, y, t) is a function of both time and space. Figure 1 shows the theoretical dispersion curves obtained by a commercial software package called Disperse (version 2.0.20a, Imperial College London, London, UK) [28]. In this paper, the incidence angle is set at 12.75 • , which is suitable for an ACT to actuate the A0 mode in an aluminum plate with the frequency-thickness product of 0.3 MHz·mm. Figure 2 shows a schematic diagram, which is used for describing the process of sound field transmission and reception.

Identification and Separation of Lamb Wave Modes in the Wavenumber Domain
Only the fundamental anti-symmetric mode (A0) can be actuated by ACT with the special incidence angle; however, mode conversion still exists due to the interaction between the Lamb wave and the flat bottom defects. The received Lamb wave signals usually comprise the incident, reflected, and converted wave modes. By using three-dimensional fast Fourier transform (3D FFT) [29], the full wavefield data can be transformed into the wavenumber domain, where various wave modes can be distinguished.
U(k x , k y , f ) = F 3D (U(x, y, t)) (1) The entire wavefield signals are filtered by the wavenumber domain filter for the purpose of removing any modes other than the A0 mode.
These filtered signals in the wavenumber domain should be transformed back into the time domain by using three-dimensional inverse fast Fourier transform (3D IFFT). U(x, y, t) is the filtered wavefield.

Acquisition of the Average Wavenumber
A spatial window is used for reducing the magnitude of the side lobes. It should be noted that the spatial window size is at least twice the wavelength of the A0 mode in this paper. The centered coordinate of this spatial window is located at (x m , y n ). Under this circumstance, the process of obtaining the average wavenumber is clearly illustrated as below.
First, a short space transform is conducted in order to retain the spatial information and reduce the magnitude of side lobes in the wavenumber domain. Specifically, the wavefield data is multiplied by a spatial window to produce the windowed dataset U mn (x, y, t). The spatial window is non-zero for only a short period in space. The Hanning window W xy is usually adopted [30].
where r and D x can be expressed as: Here, λ denotes the wavelength of the Lamb wave. D x is the window length. These signals processed with the spatial window are transformed into the spatial-frequency-wavenumber domain by using 3D FFT.
The weight of k must be considered before calculating the local wavenumber k mn , as presented in the following equation: where k is the Euclidean norm of the 2D wavenumber range [31]. The average wavenumber k mn over the selected frequency band is calculated because a real signal usually contains a certain bandwidth instead of the singular frequency component. The average wavenumber of each point can be obtained by changing the centered coordinate of this spatial window and repeating the above steps. The flow chart of the local wavenumber approach is shown in Figure 3.

Depth Reconstruction
The average wavenumber k mn at each point can be obtained by the local wavenumber approach. Hence, the phase velocity C p (x m , y n ) can be calculated as below: The relationship between the phase velocity and the plate thickness can be gained from the phase velocity dispersion curve shown in Figure 1a. Consequently, the plate thickness d(x m , y n ) of the target area can be reconstructed according to this relationship. Then, the depth distribution can be written as: where D is the thickness of the pristine plate. In order to evaluate the depth, the maximum depth error E depth is defined as below: Here, h max denotes the maximum value of the reconstructed depth and h 0 is the actual value.

Finite Element Model
A commercial software package called PZFlex was used to create the finite element model, which is shown in Figure 4. PZFlex is the time domain finite element method (FEM). The complex acoustic wave propagation in a three-dimensional model can be simulated by using this method [32]. The upper and bottom surfaces of the model are free, while the remaining four sides are absorbing. The size of the 6061-T6 aluminum plate is 400 mm × 400 mm × 1.5 mm, and the other material properties are listed in Table 1. The incident point of the ACT is (100, 200), and its incidence angle is set at 12.75 • , which is suitable to actuate the A0 mode on aluminum with the frequency-thickness product of 0.3 MHz·mm. The damage located at the middle of the plate is a circular flat bottom hole with a diameter of 10 mm and a depth of 1 mm. This damage is surrounded by a rectangular scanning area (60 mm × 60 mm), which is divided into 121 × 121 points with a spatial interval of 0.5 mm. The spatial interval should not be more than half of the wavelength according to the sampling theorem. In our test, the phase velocity of the A0 mode is about 1551 m/s, and the central frequency of the received signal is 200 kHz. Therefore, the wavelength of the A0 mode is 7.775 mm (λ = 1511/200 = 7.775 mm). The sampling interval is only 0.5 mm (λ/15). Apparently, the sampling interval satisfies the sampling theorem and can provide a good spatial resolution. The out-of-plane velocities over the scanning area are recorded by the LDV.
In this paper, a five-cycle sinusoidal tone burst modulated by a Hanning window with a central frequency of 200 kHz is applied as the excitation signal, which is shown in Figure 5. Eleven specimens listed in Table 2 are used to evaluate the performance of the local wavenumber approach from different aspects. Case 8 and Case 9 are used to investigate whether the impact which the big hole exerts on the small one is connected with their positions relative to the excitation source when they are in the same line.
Case 10 and Case 11 are used to identify the direction of the rectangular groove.

Mode Separation
The interactions between Lamb waves and the damage have been analyzed by many researchers. Bhuiyan et al. adopted the wave damage interaction coefficient (WDIC) to investigate the scattering of Lamb waves and demonstrated that the damage in the plate acted as a non-axisymmetric secondary source [33]. In another report, Bhuiyan et al. described the 3D interaction (frequency, incident direction, and azimuth direction) of Lamb waves with the damage [34]. Testoni et al. adopted a wavenumber filter to filter out the main wavefield generated by the acoustic source. Thus, the wave scattered by the delamination can be clearly observed in the filtered signal [35]. In this paper, Case 1 is taken as an example to describe the interaction between Lamb waves and the damage; moreover, the mode separation is introduced in detail.
The snapshots of the Lamb waves are shown in Figure 6. It can be observed that only the A0 mode can be actuated by an ACT with a special incidence angle. However, the S0 mode will be generated when the A0 mode encounters a flat bottom defect. Following this, the incident A0, reflected A0, transmitted A0, and S0 mode will coexist for a long time. It is worth noting that the amplitude of the S0 mode is very low. The Lamb wave propagation was analyzed in detail. The next section will be focused on mode separation. The frequency-thickness product is 0.3 MHz·mm at the undamaged area. Under this situation, the theoretical wavenumbers of the A0 mode and S0 mode are 0.8102 rad/mm and 0.2312 rad/mm, respectively. The frequency-thickness product is only 0.1 MHz·mm at the damage. Hence, the theoretical wavenumbers of the A0 mode and S0 mode are 1.3082 rad/mm and 0.2311 rad/mm, respectively. Apparently, the wavenumber of the A0 mode changes more obviously than that of the S0 mode when the thickness changes. Figure 7 shows the mode separation by using 3D FFT and the wavenumber domain filtering technique. The wavenumber spectrum is obtained at a frequency of 200 kHz. The images are given in dB scale so that the wavenumber spectrum can be clearly visualized. The wavenumber spectrum for the pristine plate is shown in Figure 7a. The wavenumber spectrum for the damaged plate is shown in Figure 7b. The high-pass filter in wavenumber domain is shown in Figure 7c. The filtered wavenumber spectrum is shown in Figure 7d. By comparing the wavenumber spectrum in Figure 7a,b, we can confirm the fact that the A0 mode has interacted with the defect. In addition, the S0 mode generated by mode conversion is very weak, as the amplitude of the S0 mode is much lower than that of the A0 mode in the time domain.
Generally, a band-pass filter should be applied to reserve the A0 mode according to the theory of the local wavenumber approach. However, the wavenumber distribution of the A0 mode is wide in our test, while the wavenumber distribution of the S0 mode is very narrow. Moreover, the intensity of the S0 mode is weak. Therefore, it is easy to design a high-pass filter rather than a band-pass filter.
The high-pass filter should ensure that the wavenumbers of the A0 mode are included in the pass band.

FEM Results of a Circular Flat Bottom Hole
In this paper, the wavelength of the A0 mode is 7.775 mm. Therefore, the window radius (D x /2) should be larger than 7.775 mm. It is assumed that the window radius is 8 mm. The average wavenumber can be obtained using Equation (12). Figure 8 shows the acquisition of the average wavenumber for Case 1. The frequency interval is 15 kHz. It can be seen that the average wavenumber contributes to improve the test quality at the scanning boundary by comparing Figure 8d,h. Moreover, lower image noise is seen in the image when the average wavenumber is adopted. In fact, the best test result is not always obtained under the condition in which the window radius is equal to the wavelength. It is meaningful to analyze the test results under various window radii in order to find the optimum size. The test results shown in Figure 9 indicate that the damage region characterized by the proposed approach will expand as the window radius increases. Meanwhile, the maximum wavenumber attained under the various radii will decrease as the window radius increases, which is shown in Figure 10. It is only the window radius of 8 mm that can make the identified damage region match well with its actual size. Therefore, 8 mm is the optimum size. All of the test results were obtained under this condition.   Table 2. The simulation results of Case 1 to Case 3 are shown in Figure 11. The size and shape of the defect can be clearly visualized. The depth estimations are shown in Table 3. The depth error of the deep hole is 1.13%, while it reaches 3.26% for the shallow hole. Apparently, the depth estimation of the deep hole is more accurate than that of the shallow one. The depth error of the big hole is 1.13%, while it reaches 8.17% for the small one. Hence, the depth estimation of the big hole is more accurate than that of the small one.  Case 4 to Case 7 contain two circular flat bottom holes with the same radius and depth, but the excitation source location and the gap between two defects are different. The specimen parameters are listed in Table 2. The simulation results of Case 4 to Case 7 are shown in Figure 12. It is observed that the shape and location of the defects are recognizable. The two holes exhibit a symmetrical distribution when the excitation source is located at (100, 200). In this condition, the better result appears under a big gap between the defects. The two holes and the excitation source are in the same line when the excitation source is located at (200, 100). In this condition, it can be noticed that the hole far away from the excitation source seems smaller than the near one. Moreover, the closer the defects are, the more obvious the effect is. Table 4 shows that the depth errors of Case 4 to Case 7 are all within 2%.  Two circular flat bottom holes with different radii and the same depth are used in Case 8 and Case 9. It is important to note that these defects and the excitation source are in the same line, which means that the defects have the different positions relative to the excitation source. The specimen parameters are listed in Table 2. Simulation results are shown in Figure 13. It can be seen that the shape and location of two holes can be identified, which demonstrates that the local wavenumber approach has a high resolution even in distinguishing adjacent defects. Table 5 shows the depth evaluations of Case 8 and Case 9. The depth error of the big hole is within 1%, while the depth error of the small one is in the range of 4% to 9%. Thus, the depth estimation of the big hole is more accurate than that of the small one. Besides, the relative position is also taken into consideration. The depth error of the small hole is 8.17% when it is in front of the big one. Apparently, this level is equivalent to that of Case 3. However, the depth error of the small hole is reduced to 4.76% when it is behind the big one. This indicates that the impact which the big hole exerts on the small one is related to their positions when they and the excitation source are in the same line.  Case 10 is a rectangular groove with an inclined angle of 90 • . Case 11 is a rectangular groove with an inclined angle of 135 • . The specimen parameters are listed in Table 2. The simulation results of Case 10 and Case 11 are shown in Figure 14. The damage shape can be clearly visualized. It should be noted that the length of this defect is a little shorter than the actual value, while its width is obviously larger than the actual value. In addition, the orientation of the rectangular groove can be easily obtained and the computational process is introduced as below. In image processing, computer vision, and related fields, an image moment is a certain particular weighted average (moment) of the image pixels' intensities. Image moments and moment invariants are useful to describe global features in image processing [36]. Image moments are calculated under grayscale in order to distinguish the direction of the rectangular groove. G(x, y) denotes the pixel intensity of the gray value, where x, y are the row and column of the pixel matrix, respectively. The gray image must be processed by a threshold as we focus attention on the area where the rectangular groove is likely to exist rather than the whole area. I(x, y) is the pixel intensity obtained by setting a threshold at the level of −3 dB. The zero moment is shown in Equation (16). The first moment is shown in Equation (17). The second moment is shown in Equation (18). These image moments can be used to extract image features [37].
M 00 , M 10 , and M 01 are used to obtain the centroid coordinates of the feature image [38]. The calculation formulas are written as below.
The spindle orientation of the feature image can be obtained by calculating an angle, which is shown in the following equation.
where µ ′ 11 , µ ′ 20 , and µ ′ 02 can be given by the following equations: Figure 15 shows the orientation of an inclined rectangular groove. The depth and angle evaluations are shown in Table 6. It is easy to find that the angle errors of Case 10 and Case 11 only reach 0.00% and 0.12%, respectively, while the depth errors of Case 10 and Case 11 rise to 12.95% and 13.60%, respectively. One reasonable explanation is clarified in the discussion section.

Experimental Setup
The experimental system is composed of an ACT, an LDV, an ultrasonic power device, a signal amplifier, a digital phosphor oscilloscope, and a data acquisition interface. The ACT is a piezoelectric ceramic-based transducer with a 200 kHz center frequency, and a bandwidth of ±25% of the center frequency. The incidence angle is set at 12.75 • , which is suitable for the ACT to actuate the A0 mode. The LDV is employed to measure the out-of-plane velocity on the surface of the plate. The sampling frequency is set at 500 MHz. The signal should be executed at a suitable time when the waveform displayed on the oscilloscope tends to be stable. Each measurement is repeated 10 times in order to reduce the measurement error. The received signals need be amplified by a wideband amplifier for the purpose of improving the signal-to-noise ratio. The gain of the signal amplifier is 60 dB. Moreover, it can be adjusted according the central frequency of the excitation signal. The data acquisition subsystem provides detailed information of the recording signals, such as the sample interval, trigger point, vertical scale, horizontal scale, etc. The experimental system is shown in Figure 16. The experimental conditions are in agreement with the simulation conditions.

Experimental Results
For the purpose of removing the noise, the received signals are filtered by a Chebyshev band-pass filter [39] with a bandwidth of 140 kHz to 260 kHz. Figure 17 shows the experimental signal in the time domain and frequency domain. Figures 18-22 are the experimental results.     The experimental results have good agreement with the simulation results. Figure 18 shows that the depth estimation of the deep hole is more accurate than that of the shallow one. In addition, the depth estimation of the big hole is more accurate than that of the small one. Figure 19 reveals that the circular flat bottom hole far away from the excitation source seems smaller than the near one when the excitation source and two identical defects are in the same line. Figure 20 shows that the proposed approach is capable of identifying the adjacent defects with different sizes. Meanwhile, it indicates that the impact that the big hole exerts on the small one is connected to their positions relative to the excitation source when they are in the same line. Figure 21 shows that the width of the rectangular groove is larger than its actual value. Figure 22 shows the direction identification for an inclined rectangular groove. Tables 7 and 8 indicate that the evaluation results still need to be improved due to the existence of multi-reflection between the adjacent defects and the spatial window effect, which will be analyzed in the discussion section.

Discussion
In this paper, the size of a spatial window is 16 mm, while the diameter of a circular flat bottom hole is no more than 10 mm. Hence, this spatial window will not only contain the propagating waves located at the damage, but also contain those located at the undamaged area. According to the dispersion curves, it should be noted that the wavenumber of the A0 mode will decrease with the increase in the frequency thickness product. More specifically, the wavenumber located at the damage is larger than that located at the undamaged area. Moreover, the deeper the circular flat bottom hole is, the more prominent the difference is. Therefore, this window effect cannot be ignored.
(1) The depth estimation of the deep hole is more accurate than that of the shallow one, as revealed by comparing the results of Case 1 and Case 2.
The frequency-thickness product of the shallow hole is much larger than that of the deep hole, which means that the wavenumber of the former is smaller than that of the latter. Additionally, there is no significant difference between the shallow hole and the undamaged area in terms of the wavenumber. According to Equation (11), it is easy to understand that the shallow hole will be easily affected by the window effect relative to the deep hole.
(2) The depth estimation for the big hole is more accurate than that of the small one. This conclusion was drawn by comparing the results of Case 1 and Case 3.
One reasonable explanation for this is that the local wavenumber is affected by the difference in size between the defect and the spatial window. The small diameter of a circular flat bottom hole means that the spatial window will contain more propagating waves located at the undamaged area where the wavenumbers are small. According to Equation (11), it is easy to understand that the local wavenumber will decrease as the lower wavenumbers are present in the spatial window. The smaller the diameter of a circular flat bottom hole is, the smaller the local wavenumber and the reconstructed depth will be.
(3) The results of Case 4 and Case 5 indicate that the better test result appears under the big gap between defects when two circular flat bottom holes are symmetrical to the excitation source.
Apparently, multi-reflection between two defects is low when they are far apart from each other. Furthermore, they are affected at the same level as they are symmetrical to the excitation source.
(4) The circular flat bottom hole far away from the excitation source seems smaller than the near one when the excitation source and two identical defects are in the same line. Moreover, the closer the two defects are, the more obvious the effect is. This conclusion was obtained by comparing the results of Case 6 and Case 7.
In order to determine the explanation for the above phenomenon, several numerical simulations were conducted under different conditions, which are listed in Table 9. The wavenumber distributions of the spatial window are shown in Figure 23. There is only one defect existing in Figure 23a-c, while there are two defects coexisting in Figure 23d-f. The wavenumber distributions of Figure 23a,d are obtained under the spatial window whose center is at (200,210). Great differences between Figure 23a,d can be found. The wavenumber distributions of Figure 23b,e are obtained under the spatial window whose center is at (200, 215). There are some differences between Figure 23b,e. The wavenumber distributions of Figure 23c,f are obtained under the spatial window whose center is at (200, 220). Only a few differences between Figure 23c,f can be found. The fundamental reason for this is that the defect near to the excitation source can change the incident wave in direction and amplitude, which dominates the distribution of the wavenumber. Furthermore, it can be clearly found that the impact which the near defect exerts on the far one is gradually reduced with the increase of the distance between them when the excitation source and two identical defects are in the same line.  Table 9. (5) The results of Case 8 and Case 9 indicate that the impact which the big hole exerts on the small one is related to their positions when the excitation source and multiple defects are in the same line.
The direction and amplitude of the incident wave can be changed obviously when the big hole is in front of the small one. Hence, the wavenumber distributions located at the small hole also change. The impact that the big hole exerts on the small one comes mainly from the reflection wave when the big hole is behind the small one. However, it is the incident wave that plays a crucial role in changing the wavenumber distributions, because the incident wave has a larger amplitude than the reflection wave. Therefore, the test result for the small hole will only change a little when the big hole is behind it.
(6) The depth errors for Case 10 and Case 11 are more than 10%. The length of the rectangular groove is larger than that of the spatial window. However, the width of the former is still far smaller than that of the latter. Hence, the local wavenumber is seriously affected by the window effect.
(7) Several measures are adopted to improve the test quality when the proposed approach is applied for an unknown damage.
In practice, both the damage size and its position are usually unknown. A large amount of data will be acquired if the entire structure is inspected. Thus, a preliminary step is desired. More specifically, it is easy to acquire the approximate location of a defect by using reliable methods, such as DAS or the compact phased array (CPA) method [40]. Then a small scanning region where the defect is likely to exist can be obtained. Under the circumstance, the local wavenumber approach can be conducted without consuming much time. The spatial window size is at least twice the wavelength of the excited mode. It should be noted that the best test result is not always obtained under the minimum size of the spatial window. In fact, the best test result is obtained under the optimum size due to the fact that the optimum size contributes to diminishing the spatial window effect. The optimum size is also fixed when the central frequency of the excitation signal and the excited mode are fixed. The optimum size can be attained by adjusting the size of the spatial window until the maximum wavenumber of the wavenumber map is acquired. Therefore, patience is important during the period of seeking the optimum size. The process of improving the test quality is shown as below.
The excitation signal with a higher central frequency is adopted as the wavelength of the excited mode is small. Under this condition, the minimum size of the spatial window may be smaller than the optimum size. The wavelength of the A0 mode is about 5 mm when the central frequency of the excitation signal is 400 kHz. Hence, the minimum radius of the spatial window is 5 mm. The radius will be increased at an interval of 1 mm. The test results obtained under various window radii are shown in Figure 24. The maximum wavenumber obtained under various window radii is shown in Figure 25. Figure 24 shows that the identified damage region will expand with the increase of the window radius. It was found that the test results are poor when the radius is small. The identified damage region is almost equal to the actual size under a special window radius (8 mm). Figure 25 shows that the maximum wavenumber increases with the size of the spatial window and gradually approaches to the theoretical value at the early stage. The maximum wavenumber is very close to the theoretical value under a special window radius (8 mm). However, the maximum wavenumber will decrease when the size of the spatial window continues to be increased.
It is worth that the wavelength of the excited mode is large when the central frequency of the excitation signal is low. Hence, it may lead to the minimum size of the spatial window being larger than the optimum size. Under this circumstance, a good test result is obtained at the minimum size of a spatial window. The test results are shown in Figures 9 and 10.  The delamination of composite laminated plates was taken as an example to illustrate the proposed approach. The interactions between Lamb waves and the delamination are complex. It is difficult to identify and isolate a specific mode for the reconstruction analysis. Therefore, the average wavenumber is less accurate. In addition, the position of the delamination cannot be quantified precisely.
In summary, the test quality can be affected by the spatial window effect and the reflection between multiple defects. It is necessary to avoid making the excitation source and multiple defects in the same line in order to reduce the interference between adjacent defects. The optimum size of a spatial window need be sought based on the test results. Increasing the central frequency of the excitation signal is a way to decrease the wavelength of the excited mode and helps to improve the test quality. However, the high frequency can easily produce complex modes. Therefore, it is important to maintain a balance between increasing the central frequency of the excitation signal and reducing the mode complexity.

Future Work
Given that the local wavenumber approach has many advantages in structural health monitoring (SHM), we will apply this approach for damage detection in adhesive-bonded concrete structures. The main work will be focused on identifying the local debonding and the weakened adhesion. A non-contact system will be composed of two low frequency air coupled transducers. One will be located at a fixed place for actuating a desired guided wave mode by setting an appropriate incidence angle. The other will be used for receiving the signals with the same inclined angle and for scanning over the target area with a spatial interval. The work will be full of challenges as separating the desired mode from multi-modes is not easy and the attenuation existing in the concrete structure may lead to a low signal-to-noise ratio. The detailed solutions will be introduced in the future work.

Conclusions
The Lamb wave local wavenumber approach for characterizing flat bottom defects (including two circular flat bottom holes and a rectangular groove) in an isotropic thin plate was presented in this paper. The hybrid and non-contact system is composed of an air-coupled transducer (ACT) and a laser Doppler vibrometer (LDV). On the one hand, the A0 mode is easily actuated by the ACT with a special incidence angle. On the other hand, the out-of-plane velocity over the scanning area can be recorded by an LDV.
Both the simulation results and the experimental results demonstrated that the damage can be characterized in shape, location, and size. Compared with other methods, the proposed approach can be used to reconstruct the thickness distribution over the scanning area when the wavenumber is combined with the phase velocity dispersion curve. Consequently, the depth of the flat bottom defect also can be quantified. In addition, adjacent defects are clearly characterized and the direction of the inclined rectangular groove can be obtained by calculating the image moments under grayscale. It must be noted that an appropriate window size contributes to improving the test quality.