Numerical Analysis of Guided Waves to Improve Damage Detection and Localization in Multilayered CFRP Panel

Multilayered carbon fiber-reinforced polymers (CFRP) are increasingly used in aircraft components because of their superior mechanical properties. However, composite materials are vulnerable to impact loads, resulting in delamination-type damage which, if unnoticed, could lead to catastrophic structural failure. The objective of this research was to investigate possibilities to improve damage detection and the localization using signal processing methods. Numerical modeling using the semi-analytical finite element (SAFE) method was performed to obtain guided wave dispersion curves and to perform modal analysis. From the modal analysis, A0 mode for inspection of the composite with delamination type defects was selected. From the numerical simulation, A0 mode interaction with delamination along the longitudinal direction was analyzed and the location of the defect was estimated by measuring the time of flight (ToF) of the signal using Hilbert transform (HT) and continuous wavelet transform (CWT). The CWT has shown better results in estimating the delamination location compared with HT. The depth of delamination was characterized in the frequency domain by comparing the amplitude of the A0 mode. Inverse fast Fourier transform (IFFT) is recommended to reconstruct the reflected and transmitted modes for better damage detection and to reduce the complexity of signal interpretation.


Introduction
Composite materials such as carbon fiber-reinforced polymers (CFRP) and glass fiberreinforced polymers (GFRP) are widely used in aerospace components, wind turbine blades, and other engineering structures [1][2][3]. The use of multilayered CFRP materials has increased significantly in the aerospace sector due to the high strength-to-weight ratio and because these materials do not corrode. Ultimately, aircraft weight can be reduced, and in return, the reduction of airframe weight allows fuel economy and a reduction of flight operation costs.
However, multilayered CFRP structures are vulnerable to delamination and fiber breakages. Delamination-type defects due to impact damage can occur during assembly and exploitation. Key issues in aircraft wing structures are delamination and disbond-type damage. Such types of defects are difficult to detect during inspection, and this can lead to catastrophic structural failure.
In aerospace applications, structural health monitoring (SHM) could increase the safety, reliability, and quality control of the structure [2]. However, there is no universal SHM technique that can effectively solve all problems. Generally, each SHM technique has its advantages and disadvantages depending on the problem. Among potential SHM techniques such as vibration-based and acoustic emission, ultrasonic-guided wave techniques are widely investigated due to their better damage detection and localization [4]. Different

Modeling Setup
The CFRP plate (USN150B) of 200 mm in length and 3.5 mm in thickness was selected for analysis. The CFRP panel was constructed from 18 plies, and the laminate orientation was (-45 • , 0 • , +45 • , 0 • , 90 • , 0 • , -45 • , 0 • , and +45 • ) 2 . The normal force was applied on an excitation zone of 8.4 mm to excite A 0 mode by setting the transmitter-receiver configuration. The specific length of the excitation zone obtained a maximum amplitude of A 0 mode. In this study, delamination-type damages were analyzed, and the length of the delamination was 20 mm. The excitation zone was at the leading edge of the plate, while three receivers were positioned at 70, 80, and 90 mm from the leading edge (X r ). Two cases of delamination positions were considered: the delamination position (X d ) at 40 mm, i.e., between the excitation zone and the receivers, and at 120 mm, i.e., further than the receivers. Cases studied, positions, and distances are presented in Table 1. The defect-free CFRP plate is presented in Figure 1a. In case 1, the delamination is positioned in between the excitation zone and the receiver, as shown in Figure 1b. In case 2, the delamination is positioned after the receiver, as shown in Figure 1c. The mechanical properties of the CFRP (USN150B) used in the numerical study are presented in Table 3. Here, E is the Young's modulus, G is the shear modulus, is the Poisson's ratio, is the density of CFRP, and is the ply thickness. In this 2D model, the directions X and Y are defined as the X-axis along the fibers and the Y-axis along the thickness of the plate [21]. Table 3. The mechanical properties of the USN150B carbon fiber-reinforced polymer (CFRP) material [22]. E1 (GPa) 131 E2 = E3 (GPa) 8 In another investigation, the delamination position in the through-thickness direction was changed. The delamination position (X d ) was fixed at 100 mm, and the delamination location in the through-thickness direction (X h ) was modified: 1.1, 1.7, and 2.3 mm. The cases investigated are presented in Table 2. The mechanical properties of the CFRP (USN150B) used in the numerical study are presented in Table 3. Here, E is the Young's modulus, G is the shear modulus, ν is the Poisson's ratio, ρ c is the density of CFRP, and t ply is the ply thickness. In this 2D model, the directions X and Y are defined as the X-axis along the fibers and the Y-axis along the thickness of the plate [21].  Table 3. The mechanical properties of the USN150B carbon fiber-reinforced polymer (CFRP) material [22].

Guided Wave Dispersion Curves
It is essential to obtain dispersion curves to select suitable modes and the central frequency for the inspection [23][24][25][26]. The semi-analytical finite element (SAFE) method is one of the most common approaches to calculate eigenvalues of guided wave modes in arbitrary cross-sectional waveguides [27]. In this method, the one-dimensional crosssection of the plate is discretized, and an analytical solution is applied along the wave propagation direction. The governing equation for wave propagation in an elastic medium can be expressed by the following equation: where λ and µ are the Lamé constants, ρ is the density, and t is the time. The u x , u y , and u z are in-plane, shear horizontal, and out-of-plane displacement vectors. In this analysis, the length and width of the plate were assumed to be infinitely long to consider the guided wave mode as symmetrically distributed in the throughthickness direction. Eighteen elements across the plate were assigned, and the thickness of one element replicates one ply of the multilayered plate. The thickness of each ply was 0.194 mm and the total thickness of the 18 plies of the CFRP plate was 3.5 mm.
The SAFE method was derived to obtain the dispersion curves of guided wave phase velocity, c ph ( f ), and group velocity, c g ( f ), dispersion curves. The phase velocity versus frequency and group velocity versus frequency dispersion curves are presented in Figure 2a,b. From the dispersion curves, it follows that at a 100 kHz frequency, the phase and group velocities of A 0 mode were 1369 and 1625 m/s, respectively. The selected frequency range does not obtain higher-order modes and enables only fundamental modes. The maximum amplitude of each mode varies depending on the frequency. In this analysis, the guided wave mode at higher frequency had a smaller wavelength, but in return, attenuation increased, and superposition of different modes led to complex signal interpretation due to dispersion. 2a,b. From the dispersion curves, it follows that at a 100 kHz frequency, the phase and group velocities of A mode were 1369 and 1625 m/s, respectively. The selected frequency range does not obtain higher-order modes and enables only fundamental modes. The maximum amplitude of each mode varies depending on the frequency. In this analysis, the guided wave mode at higher frequency had a smaller wavelength, but in return, attenuation increased, and superposition of different modes led to complex signal interpretation due to dispersion.

Guided Wave Modal Analysis
The purpose of the SAFE method is not only to calculate dispersion curves, but in addition, to analyze displacement distribution in the through-thickness of the plate and to select a particular mode, which would be sensitive to delamination-type defects. In guided waves, out-of-plane displacements are dominant in asymmetrical modes, and inplane displacements are dominant in symmetrical modes.
The in-plane and out-of-plane displacements of A mode in a 3.5 mm-thick CFRP plate are presented in Figure 3. The maximum amplitude of in-plane displacements of A mode at 100 kHz was relatively small compared to out-of-plane displacements of A mode. Therefore, the maximum amplitude of out-of-plane displacements of the A mode was dominant and relatively higher compared to in-plane displacements.

Guided Wave Modal Analysis
The purpose of the SAFE method is not only to calculate dispersion curves, but in addition, to analyze displacement distribution in the through-thickness of the plate and to select a particular mode, which would be sensitive to delamination-type defects. In guided waves, out-of-plane displacements are dominant in asymmetrical modes, and in-plane displacements are dominant in symmetrical modes.
The in-plane and out-of-plane displacements of A 0 mode in a 3.5 mm-thick CFRP plate are presented in Figure 3. The maximum amplitude of in-plane displacements of A 0 mode at 100 kHz was relatively small compared to out-of-plane displacements of A 0 mode. Therefore, the maximum amplitude of out-of-plane displacements of the A 0 mode was dominant and relatively higher compared to in-plane displacements.  The in-plane and out-of-plane displacements of the S mode in a 3.5 mm-thick CFRP plate are presented in Figure 4. In this analysis, the maximum amplitude of in-plane and out-of-plane displacements of the S mode at 100 kHz was relatively small compared to A mode.  The in-plane and out-of-plane displacements of the S 0 mode in a 3.5 mm-thick CFRP plate are presented in Figure 4. In this analysis, the maximum amplitude of in-plane and out-of-plane displacements of the S 0 mode at 100 kHz was relatively small compared to A 0 mode. The in-plane and out-of-plane displacements of the S mode in a 3.5 mm-thick CFRP plate are presented in Figure 4. In this analysis, the maximum amplitude of in-plane and out-of-plane displacements of the S mode at 100 kHz was relatively small compared to A mode. The mode was selected for the inspection of the component with delaminationtype defects, as the maximum amplitude of particle displacements of A mode was significantly higher compared to the S mode. In addition, the asymmetric mode was sensitive due to its displacement distribution in the thickness of the plate and A mode velocity changes over delamination. The out-of-plane displacements of the A mode were uniformly distributed throughout the thickness of the plate. The FEA of the guided wave The A 0 mode was selected for the inspection of the component with delamination-type defects, as the maximum amplitude of particle displacements of A 0 mode was significantly higher compared to the S 0 mode. In addition, the asymmetric mode was sensitive due to its displacement distribution in the thickness of the plate and A 0 mode velocity changes over delamination. The out-of-plane displacements of the A 0 mode were uniformly distributed throughout the thickness of the plate. The FEA of the guided wave propagation in the CFRP plate was performed to analyze guided wave interactions with delamination at different positions and depths, as presented in the following sections.

Numerical Simulation of Guided Waves
Experimental inspection of every guided wave propagation scenario for different damages, loading, and environmental conditions requires many mock-up trials. Therefore, physics-based numerical modeling is an efficient way to understand the mechanism of guided wave propagation and interactions with delamination in composite structures.
The guided wave interaction with delamination was analyzed with the transmitterreceiver configuration, as shown in Figure 1. A numerical simulation was performed to analyze the sensitive mode to estimate the position of the delamination in the longitudinal and through-thickness directions. The simulation was performed in the 2D approach, where the plate is infinite along the z-axis. The multilayered CFRP panels were anisotropic and the GW propagation velocity was angle-dependent. The 3D modeling could give propagation in all directions, which would be closer to a real investigation. However, 3D modeling of ultrasonic wave propagation, taking into account 20 elements per wavelength criterion, requires huge computational and time resources. In this research, numerical modeling was performed to investigate the interaction of guided waves with delaminationtype defects at different depths. Therefore, modeling of guided wave propagation using the 2D plane strain approach was most suitable for the given task. The simulation was carried out with 2100 time-steps and the duration of each step was 0.1 µs. The total modeled duration of the A 0 mode propagation was 210 µs.
In FEA, mesh shape and size are important factors. The CPE4R-type elements were selected for the analysis in ABAQUS explicit. These types of elements are plain-strain quadrilateral and allow less integration time, and better hourglass control. In FEA, hourglass control attempts to minimize severe mesh distortion and provides better mesh control. On the other hand, a smaller mesh size improves the accuracy of the analysis but increases the computation time. In guided wave simulation, the mesh size is significantly dependent on the number of elements per wavelength (N/λ). To maintain a good balance between accuracy and computation time, 20 elements per wavelength (λ) were selected for the mesh calculation. Therefore, the mesh size was 0.6 mm.
The fundamental asymmetrical (A 0 ) mode was excited by applying normal force on the excitation zone. The input signal was a five-cycle sine wave with a Gaussian envelope at a driving frequency of 100 kHz.
Guided wave analysis at 100 kHz is an appropriate way to achieve a good balance between dispersion, mode complexity, and attenuation. In order to properly identify A 0 mode, the signals on the surface of the plate were recorded and corresponding B-scan data were obtained. Using 2D Fourier transform, B-scan data were converted into phase velocity versus frequency domain, as presented in Figure 5. The dispersion characteristics of the guided waves propagating in the CFRP plate from B-Scan data were obtained using the two-dimensional Fourier transform. The propagated guided wave along the plate is expressed as a function of the coordinates and time u(x,t), and can be transformed into the wavenumber, k, and the frequency space using 2D FFT [28]: where x is the coordinate, t is the time, ω is the angular frequency, and k is the wavenumber. The dispersion characteristics of the guided waves propagating in the CFRP plate from B-Scan data were obtained using the two-dimensional Fourier transform. The propagated guided wave along the plate is expressed as a function of the coordinates and time, u(x,t), and can be transformed into the wavenumber, k, and the frequency space using 2D FFT [28]:

Guided Wave Interaction with Delamination
where x is the coordinate, t is the time, ω is the angular frequency, and k is the wavenumber.

Guided Wave Interaction with Delamination
When the guided wave travels along the sample and encounters a delamination, part of excited A 0e mode is reflected at the leading edge of the delamination and part is converted to S 0c mode, as shown in Figure 6. A vertical, out-of-plane displacement component of the propagated wave was recorded in the study, as it is the dominant component of the A 0 mode. Meanwhile, the vertical displacement component of the S 0 mode is practically close to zero (see Figure 4). While traveling over the delamination, the converted S 0c mode was converted back to A 0t mode at the trailing edge of the delamination. Finally, strong reflected (A 0r ) and transmitted (A 0t ) modes were obtained after interacting with the delamination. Therefore, only reflected (A 0r ), and transmitted (A 0t ) modes can be observed in the B-scan images.

Localization of Delamination Location Using Hilbert Transform (HT)
The guided wave interaction with delamination at 40 and 120 mm was analyzed. T analyze the delamination position (Xd), the signals were acquired at three positions of th receivers, as shown in Figure 1. In the case when delamination was in between the excit tion zone and the receiver, i.e., Xd was 40 mm, the guided wave signals were received (X at 70, 80, and 90 mm with respect to the leading edge of the CFRP plate. The delaminatio

Localization of Delamination Location Using Hilbert Transform (HT)
The guided wave interaction with delamination at 40 and 120 mm was analyzed. To analyze the delamination position (Xd), the signals were acquired at three positions of the receivers, as shown in Figure 1. In the case when delamination was in between the excitation zone and the receiver, i.e., Xd was 40 mm, the guided wave signals were received (Xr) at 70, 80, and 90 mm with respect to the leading edge of the CFRP plate. The delamination location can be estimated by analyzing the reflected wave from the delamination.
In this case, the A mode propagates to the receiver, which is considered as direct A mode (t1), and the reflection from the delamination must be considered as t2. The delamination position can be estimated considering the time of flight (ToF) between t1 and

Localization of Delamination Location Using Hilbert Transform (HT)
The guided wave interaction with delamination at 40 and 120 mm was analyzed. To analyze the delamination position (X d ), the signals were acquired at three positions of the receivers, as shown in Figure 1. In the case when delamination was in between the excitation zone and the receiver, i.e., X d was 40 mm, the guided wave signals were received (X r ) at 70, 80, and 90 mm with respect to the leading edge of the CFRP plate. The delamination location can be estimated by analyzing the reflected wave from the delamination.
In this case, the A 0 mode propagates to the receiver, which is considered as direct A 0 mode (t 1 ), and the reflection from the delamination must be considered as t 2 . The delamination position can be estimated considering the time of flight (ToF) between t 1 and t 2 . The Hilbert envelope of the guided wave signals received at different receiver positions are presented in Figure 8. The envelope of the guided wave signal obtained using Hilbert transform, H(t), can be expressed using the following equation [29][30][31][32][33]: The distance between the receiver and delamination can be estimated by calculating ToF between and [34]: ∆ = − is the time difference between direct and reflected A mode, and is the group velocity.
The location of delamination was estimated by measuring the time difference between direct and reflected A mode from delamination. If the delamination is located between the excitation zone and the receiver, the excited A mode is reflected from the delamination leading edge, then from the leading edge of the sample, and only then does it propagate to the receiver (Figure 1b). Therefore, guided wave mode approaches delamination multiple times until it reaches the receiver. The estimation of the delamination location is presented in Table 4. In this case, the estimated distance is the delamination position from the leading edge. The relative error between the actual and the measured value can be calculated by Equation (5): where is the relative error, and Da and Dm are the actual and measured values, respectively.  The distance between the receiver and delamination can be estimated by calculating ToF between t 1 and t 2 [34]: ∆t = t 2 − t 1 is the time difference between direct and reflected A 0 mode, and c g is the group velocity. The location of delamination was estimated by measuring the time difference between direct and reflected A 0 mode from delamination. If the delamination is located between the excitation zone and the receiver, the excited A 0 mode is reflected from the delamination leading edge, then from the leading edge of the sample, and only then does it propagate to the receiver (Figure 1b). Therefore, guided wave mode approaches delamination multiple times until it reaches the receiver. The estimation of the delamination location is presented in Table 4. In this case, the estimated distance is the delamination position from the leading edge. The relative error between the actual and the measured value can be calculated by Equation (5): where ε R is the relative error, and D a and D m are the actual and measured values, respectively. In the case of X d = 120 mm, the guided wave signals were received at 70, 80, and 90 mm with respect to the leading edge of the CFRP plate. The distance between delamination and the receiver (X dr ) was 50, 40, and 30 mm. In this approach, if the delamination is located after the receiver position, the excited A 0 mode is received at the receiver as the first wave packet (t 1 ) and the first reflection from the leading edge of the delamination is received as the second wave packet (t 2 ) (Figure 1c), as shown in Figure 9.
In this case, the estimated distance is the distance between delamination and the receiver. The delamination position (X d ) can be estimated by adding the distance between the edge of the sample and the receiver and the estimated distance between the receiver and delamination (X e ). The localization of delamination and the relative error of this approach are presented in Table 5. first wave packet (t1) and the first reflection from the leading edge of the delamination is received as the second wave packet (t2) (Figure 1c), as shown in Figure 9. In this case, the estimated distance is the distance between delamination and the receiver. The delamination position (Xd) can be estimated by adding the distance between the edge of the sample and the receiver and the estimated distance between the receiver and delamination (Xe). The localization of delamination and the relative error of this approach are presented in Table 5.

Localization of Delamination Using Continuous Wavelet Transform (CWT)
The wavelet analysis is one of the well-known signal processing approaches for nonstationary signals. In the CWT, a wavelet, , is the analyzing function. This method compares the signal with the shift and dilation of a wavelet [35][36][37][38][39][40]. Here, dilation refers to stretched or compressed versions of the wavelet. It provides a linear time-frequency representation by pre-selecting the mother wavelet (Morse), . The Morse wavelet has the advantage of determining both time and frequency components at once due to the better symmetry and time-bandwidth product of the mother wavelet. The CWT of the analytic signal to represent in the time-frequency domain can be defined by the following equation [10,41,42]: where CWT denotes continuous wavelet transform, t is the time, f is the frequency, and

Localization of Delamination Using Continuous Wavelet Transform (CWT)
The wavelet analysis is one of the well-known signal processing approaches for nonstationary signals. In the CWT, a wavelet, ψ, is the analyzing function. This method compares the signal with the shift and dilation of a wavelet [35][36][37][38][39][40]. Here, dilation refers to stretched or compressed versions of the wavelet. It provides a linear time-frequency representation by pre-selecting the mother wavelet (Morse), ψ. The Morse wavelet has the advantage of determining both time and frequency components at once due to the better symmetry and time-bandwidth product of the mother wavelet. The CWT of the analytic signal to represent in the time-frequency domain can be defined by the following equation [10,41,42]: where CWT denotes continuous wavelet transform, t is the time, f is the frequency, and f s is the sampling frequency. The inspection setup to estimate the delamination location is presented in Figure 1b,c. The ToF of reflected A 0 mode from different delamination locations in the time-frequency domain is presented in Figure 10a The guided wave mode response obtained from the delamination position (Xd) at 40 mm is post-processed and the estimated location using CWT is presented in Table 6. The guided wave mode response obtained from the delamination position (Xd) at 120 mm and the signals received at different receiver positions to estimate the delamination location using CWT are presented in Table 7. The location of the delamination can be estimated using Equation (3). The maximum amplitude of the wave packet was used for the analysis. Therefore, the group velocity of The guided wave mode response obtained from the delamination position (X d ) at 40 mm is post-processed and the estimated location using CWT is presented in Table 6. The guided wave mode response obtained from the delamination position (X d ) at 120 mm and the signals received at different receiver positions to estimate the delamination location using CWT are presented in Table 7.
The location of the delamination can be estimated using Equation (3). The maximum amplitude of the wave packet was used for the analysis. Therefore, the group velocity of A 0 mode was used to estimate the location of the delamination. The accuracy of the measured delamination location improved using CWT compared with HT, the comparison is presented in Table 8.
The defect location measurement using CWT provided more accurate results compared with HT. Furthermore, the time-frequency distribution of the guided wave signal along with the amplitude were obtained simultaneously. The proposed method has the potential to measure the ToF of a signal by considering the excitation frequency as a reference even in the presence of other frequency content in the signal.

Characterization of Delamination at Different Depths
In this section, the delamination position in the through-thickness direction was examined. The delamination position in the longitudinal direction was 100 mm, and the positions at different depths of 1.1, 1.7, and 2.3 mm were analyzed in the CFRP plate.

Mode Separation
The amplitude of the transmitted and reflected guided wave mode varies depending on the position and size of the defect. Therefore, it is necessary to examine transmitted and reflected modes separately. In order to separate superimposed modes, the signal processing method was applied on the received guided wave signal. The flow chart of the

Mode Separation
The amplitude of the transmitted and reflected guided wave mode varies depending on the position and size of the defect. Therefore, it is necessary to examine transmitted and reflected modes separately. In order to separate superimposed modes, the signal processing method was applied on the received guided wave signal. The flow chart of the signal processing method to separate reflected and transmitted modes is presented in Figure 12. In this analysis, the inverse fast Fourier transform (IFFT) method was used to characterize mode conversion, reflection, and transmission at the leading and trailing edges of delamination.

Mode Separation
The amplitude of the transmitted and reflected guided wave mode varies depending on the position and size of the defect. Therefore, it is necessary to examine transmitted and reflected modes separately. In order to separate superimposed modes, the signal processing method was applied on the received guided wave signal. The flow chart of the signal processing method to separate reflected and transmitted modes is presented in Figure 12. In this analysis, the inverse fast Fourier transform (IFFT) method was used to characterize mode conversion, reflection, and transmission at the leading and trailing edges of delamination. Initially, a two-dimensional fast Fourier transform (2D FFT) was applied to the Bscan data to obtain the wavenumber of the reflected and transmitted modes. The wavenumber of the signal using 2D fast Fourier transform can be expressed by the following equation [28]: where t is the guided wave propagation time and x is the propagation distance. The bandpass filter with a cosine-tapered window function was applied to the frequency-wavenumber domain to separate the modes. The wavenumber of reflected A mode was isolated and the filtered mode was reconstructed in the B-scan representation using IFFT. The same procedure was followed to reconstruct the transmitted mode as well. The reflected and transmitted modes are separated and reconstructed using IFFT [28]: where B(f, k) is the 2D bandpass filter of the selected mode. The results obtained using IFFT to separate the guided wave modes are presented in Figure 13a-e. In the cases of delamination at 1.1 and 2.3 mm, the A mode interacted with delamination and the top/bottom surface of the plate simultaneously. Therefore, an additional A mode appeared in the transmitted mode, as shown in Figure 13b,f. In the case of delamination at 1.7 mm, there was a normal interaction with the top/bottom surface of the plate. Therefore, only the transmitted mode was observed. Thus, the signal processing method has shown a significant improvement to differentiate the additional A mode obtained from the 1.1 and 2.3 mm delamination in the 3.5 mm CFRP plate. Initially, a two-dimensional fast Fourier transform (2D FFT) was applied to the B-scan data to obtain the wavenumber of the reflected and transmitted modes. The wavenumber of the signal using 2D fast Fourier transform can be expressed by the following equation [28]: where t is the guided wave propagation time and x is the propagation distance. The bandpass filter with a cosine-tapered window function was applied to the frequencywavenumber domain to separate the modes. The wavenumber of reflected A 0 mode was isolated and the filtered mode was reconstructed in the B-scan representation using IFFT. The same procedure was followed to reconstruct the transmitted mode as well. The reflected and transmitted modes are separated and reconstructed using IFFT [28]: where B(f, k) is the 2D bandpass filter of the selected mode. The results obtained using IFFT to separate the guided wave modes are presented in Figure 13a-e. In the cases of delamination at 1.1 and 2.3 mm, the A 0 mode interacted with delamination and the top/bottom surface of the plate simultaneously. Therefore, an additional A 0 mode appeared in the transmitted mode, as shown in Figure 13b,f. In the case of delamination at 1.7 mm, there was a normal interaction with the top/bottom surface of the plate. Therefore, only the transmitted mode was observed. Thus, the signal processing method has shown a significant improvement to differentiate the additional A 0 mode obtained from the 1.1 and 2.3 mm delamination in the 3.5 mm CFRP plate.

Distinguishing Near-Surface Delamination
Guided waves are frequency-dependent, and the amplitude of each mode varies o different frequency ranges [43]. In this analysis, the reflected and transmitted A mo response was presented in the frequency spectrum to estimate the amplitude. The quency spectrum of reflected A mode interacting with delamination at different dep was compared with the defect-free signal.
The amplitude of A mode versus frequency can be determined from the filte two-dimensional spectrum presented in the following equation [28]: The frequency spectrum of reflected A mode compared with the defect-free sig is presented in Figure 14. The reflected signal amplitude was reduced by 31% and 3 compared to the defect-free signal when the delamination was present near the top a bottom surfaces (at 1.1 mm depth), respectively. It is an indication of strong reflection, a

Distinguishing Near-Surface Delamination
Guided waves are frequency-dependent, and the amplitude of each mode varies over different frequency ranges [43]. In this analysis, the reflected and transmitted A 0 mode response was presented in the frequency spectrum to estimate the amplitude. The frequency spectrum of reflected A 0 mode interacting with delamination at different depths was compared with the defect-free signal.
The amplitude of A 0 mode versus frequency can be determined from the filtered two-dimensional spectrum presented in the following equation [28]: The frequency spectrum of reflected A 0 mode compared with the defect-free signal is presented in Figure 14. The reflected signal amplitude was reduced by 31% and 32% compared to the defect-free signal when the delamination was present near the top and bottom surfaces (at 1.1 mm depth), respectively. It is an indication of strong reflection, and is due to the delamination edge and surface of the CFRP plate. It is important to note that the delamination position at 1.1 and 2.3 mm has the same distance from the top/bottom of the plate. Therefore, the amplitude of the reflected A 0 mode is similar. In the case where delamination presented at the center (1.7 mm) of the plate, the signal amplitude was reduced by 48%. When A 0 mode propagated over the delamination without the interference of the boundaries, a weak reflection occurred. delamination presented at the center (1.7 mm) of the plate, the signal amplitude was reduced by 48%. When A mode propagated over the delamination without the interference of the boundaries, a weak reflection occurred. Considering the transmitted mode, the signal amplitude was reduced by 60% and 50% when the delamination was present at the top and bottom surfaces (1.1 mm depth), respectively. As mentioned above, the delamination at 1.1 and 2.3 mm was at the same distance from the top/bottom of the plate. Hence, the transmitted A mode amplitude was similar for 1.1 and 2.3 mm, as shown in Figure 15. In the case where the delamination presented in the center (1.7 mm) of the plate, the signal amplitude was reduced by 18.5%. A mode propagated over the delamination with a normal interaction with the top/bottom surface of the plate.

Discussion
The guided wave dispersion curves and modal analysis were performed using the one-dimensional SAFE method. In this analysis, the guided wave was bounded by two surfaces and propagated in a single direction. However, guided wave modes have displacements in three directions and are dominant in a single direction. Therefore, from the analysis, A mode had a dominant displacement in the y-direction and had a higher amplitude compared to S mode.
The plain-strain numerical simulation was performed on a multilayered CFRP panel to analyze the guided wave mode response with delamination. In this analysis, two dif- Considering the transmitted mode, the signal amplitude was reduced by 60% and 50% when the delamination was present at the top and bottom surfaces (1.1 mm depth), respectively. As mentioned above, the delamination at 1.1 and 2.3 mm was at the same distance from the top/bottom of the plate. Hence, the transmitted A 0 mode amplitude was similar for 1.1 and 2.3 mm, as shown in Figure 15. In the case where the delamination presented in the center (1.7 mm) of the plate, the signal amplitude was reduced by 18.5%. A 0 mode propagated over the delamination with a normal interaction with the top/bottom surface of the plate. delamination presented at the center (1.7 mm) of the plate, the signal amplitude was reduced by 48%. When A mode propagated over the delamination without the interference of the boundaries, a weak reflection occurred. Considering the transmitted mode, the signal amplitude was reduced by 60% and 50% when the delamination was present at the top and bottom surfaces (1.1 mm depth), respectively. As mentioned above, the delamination at 1.1 and 2.3 mm was at the same distance from the top/bottom of the plate. Hence, the transmitted A mode amplitude was similar for 1.1 and 2.3 mm, as shown in Figure 15. In the case where the delamination presented in the center (1.7 mm) of the plate, the signal amplitude was reduced by 18.5%. A mode propagated over the delamination with a normal interaction with the top/bottom surface of the plate.

Discussion
The guided wave dispersion curves and modal analysis were performed using the one-dimensional SAFE method. In this analysis, the guided wave was bounded by two surfaces and propagated in a single direction. However, guided wave modes have displacements in three directions and are dominant in a single direction. Therefore, from the analysis, A mode had a dominant displacement in the y-direction and had a higher amplitude compared to S mode.
The plain-strain numerical simulation was performed on a multilayered CFRP panel to analyze the guided wave mode response with delamination. In this analysis, two different inspection setups were used to locate the delamination. In one case, the delamina-

Discussion
The guided wave dispersion curves and modal analysis were performed using the onedimensional SAFE method. In this analysis, the guided wave was bounded by two surfaces and propagated in a single direction. However, guided wave modes have displacements in three directions and are dominant in a single direction. Therefore, from the analysis, A 0 mode had a dominant displacement in the y-direction and had a higher amplitude compared to S 0 mode.
The plain-strain numerical simulation was performed on a multilayered CFRP panel to analyze the guided wave mode response with delamination. In this analysis, two different inspection setups were used to locate the delamination. In one case, the delamination was located between the excitation zone and the receiver. In the other case, the delamination was located after the receiver. In the conventional way, the defect location can be estimated by measuring the ToF between the converted and transmitted modes. In this analysis, the converted S 0 mode had very low amplitude. Therefore, it was not possible to locate the delamination by evaluating converted and transmitted modes.
In guided wave inspection, the reflected mode could be useful for delamination localization when the converted mode does not provide any information due to the low amplitude. The location of the delamination can be estimated using ToF between direct A 0 mode obtained from the excitation zone and reflected A 0 mode from delamination.
In this approach, the HT and CWT were used to calculate the ToF between the direct (t 1 ) and the reflected (t 2 ) A 0 mode. The HT obtains the envelope of the signal, and the ToF of the signal can be estimated by measuring the peak-to-peak amplitude of the t 1 and t 2 . The CWT obtains the time-frequency and maximum amplitude of the signal simultaneously. The ToF between the maximum amplitudes of the two wave packets allows to estimate the defect location. In addition, the CWT has the potential to achieve accurate local transients in the signal.
ToF was similar in cases when the delamination was positioned between the excitation zone and the receiver, or after the receiver. Using three receivers at different positions enabled to determine the delamination location. In the case when the delamination was in between the excitation zone and the receivers, the estimated distances for all the receiver positions were similar and showed the delamination position from the leading edge. Meanwhile, in the case when the delamination was further than the receiver, the estimated distances for the three receiver positions were different and showed the delamination position from the respective receivers.
When analyzing near-surface delamination at 1.1 and 2.3 mm in the through-thickness direction, strong reflections were obtained due to the delamination edge and the top/bottom surface of the CFRP plate. In the case of delamination at a depth of 1.7 mm, weak reflection was obtained.
When considering the transmitted mode, a weak transmitted mode was obtained due to the delamination edge and the top/bottom surface of the CFRP plate. The presented method can be used to determine the characterization of near-surface delamination, but it has limitations when it comes to measuring the exact delamination depth.

Conclusions
The A 0 mode was selected for the delamination inspection because the dominant and maximum amplitude of particle displacements of A 0 mode is significantly higher compared to S 0 mode. In addition, the asymmetrical mode is sensitive to delamination-type damages due to its displacement distribution through the thickness of the plate and guided wave mode velocity shift as it approaches delamination.
In this investigation, the inspection setup was proposed to detect and localize the delamination. This method offers to locate delamination when the converted modes do not have a significant amplitude. The reflected mode provides useful information to precisely locate the delamination. The localization of the delamination was estimated using six different scenarios to assess the accuracy of the proposed inspection setup. All inspection scenarios have shown consistent results and obtained an accurate estimate of the location of the delamination, with a relative error of less than 7.5%. Using three receivers at different positions allowed to determine where delamination was located, either in between the excitation zone and the receivers, or further than the receivers.
In this article, signal processing methods to improve damage detection and localization using guided waves have been proposed. The localization of delamination using HT and CWT was compared, where the relative error using HT and CWT was 2.5% and 1.2%, respectively. CWT obtained more accurate results to estimate the location of the defect compared to HT. In addition, CWT provided better linear time-frequency representation, higher resolution of time-frequency components, a finer spectrum, and it is superior for localizing the transients of the guided wave signal.
While characterizing near-surface delamination at 1.1 and 2.3 mm in the throughthickness direction, a strong reflection was obtained from the delamination edge and the top/bottom surface of the CFRP plate. Such a mode response occurred as a result of the variation in propagation of mode speed along the upper and lower surfaces of the delamination. In the case of 1.7 mm, the distance between the upper and lower surfaces of the plate was identical. Therefore, the mode propagation speed was identical along the upper and lower surfaces of the delamination.
Damage detection was improved using 2DFFT. When characterizing near-surface delamination at 1.1 and 2.3 mm, consecutive A 0 modes were obtained due to the boundaries and delamination. Since the consecutive A 0 modes were superimposed, it is difficult to differentiate the additional A 0 mode in the B-scan image from the reflected and transmitted modes. The additional A 0 mode was revealed after separating the reflected and transmitted modes and reconstructed using the IFFT. Thus, the complexity of signal interpretation was reduced, and damage detection was improved.
Future work will focus on the experimental verification of the proposed approach, which will be performed as the next step of this investigation.