# A Time-Frequency Research for Ultrasonic Guided Wave Generated from the Debonding Based on a Novel Time-Frequency Analysis Technique

Authors: Junhua Wu; Zheshu Ma; Yonghui Zhang
Journal: Shock and Vibration
 (2017)

---

## Abstract

Carbon fibre composites have a promising application in the future of the vehicle, because of their high strength and light weight. Debonding is a major defect of the carbon fibre composite. The time-frequency analysis is fundamental to identify the defect on ultrasonic nondestructive evaluation and testing. In order to obtain the instantaneous frequency and the peak time of modes of the ultrasonic guided wave, an algorithm based on the Smoothed Pseudo Wigner-Ville distribution and the peak-track algorithm is presented. In the algorithm, a masking step is proposed, which can guarantee that the peak-track algorithm can automatically exact the instantaneous frequency and the instantaneous amplitude of different modes on the Smoothed Pseudo Wigner-Ville distribution. An experiment for detecting the debonding for a type of carbon fibre composite is done. The presented algorithm is employed on the experimental signals. The processed result of experimental signals reveals that the defect can stimulate new modes, and there is a quantitative relationship between the defect size and the frequency of the new mode. The presented technique provides a valuable way to detect the presentence, calculate the size, and locate the position of the debonding defect.


---

## Body

## 1. Introduction

Carbon fibre composite is a type of composite material with a high strength. The specific strength and specific stiffness are better than other lightweight materials. Therefore, it is widely employed in aerospace domain and military products. The material can obviously reduce the weight of automobile. The material has been applied on construction elements and cover parts of transmission system of chassis and automobile bodies. Debonding defect is a major problem of the carbon fibre composites. A great number of scholars have done investigations of the nondestructive evaluation and testing (NDE/NDT) for the defect [1–4].At present, ultrasonic guided wave (GW) testing has emerged as a very prominent active scheme. It offers an effective method to estimate the location, severity, and type of defects and has been used on NDE/NDT of the carbon fibre composites [2, 5, 6]. Separating different modes of the GW and obtaining the corresponding time-frequency distributions (TFD) are vital for the inspection of the defect.A number of time-frequency analysis techniques have been explored for analysis of GW signals. Prosser et al. [7] used a pseudo Wigner-Ville distribution to process GW signals for material characterization of composites. Kercel et al. [8] used Bayesian parameter estimates to separate multiple modes in GW signals obtained from laser ultrasonic on a workpiece manufacturing assembly line. Mahadev et al. [9] employed the short time Fourier transform (STFT) to extract a suitable parameter for a tomographic image reconstruction mapping structural defects. Osegueda et al. [10] used Hilbert–Huang transform (HHT) to process GW signals in plate structures. This technique allows for separating of the GW signal into intrinsic mode functions (IMF) and a residue. Hilbert transform is done for IMFs to determine the energy time signal of each mode, enabling the easy location and characterization of the notch. The wavelet analysis is widely used [11–17]; a good many works were done about applications of wavelet transform (WT) for the GW signals. Dhital and Lee [13] used WT for detection of cracks, which was also realized. Paget et al. [14] proposed a new damage-detection technique based on WT with a new basis. Y. Y. Kim and E.-H. Kim evaluated the effectiveness of WT analysis for studying the wave dispersion [16].The instantaneous frequency (IF) and the peak time of the instantaneous amplitude (IA) are significant parameters of the GW mode. Tracking the IF change may provide way to identify the presence and the size of the defect. Peak time can evaluate the distance between the defect and the probe; on this basic, we can locate the defect. HHT can realize these operations. But the frequency resolution is a limitation. Ref. [18] reveals that when the ratio between a relatively low frequency and a relatively high frequency becomes close to 0.75 (low frequency/high frequency) two components of signals cannot be separated.Although Ensemble Empirical Mode Decomposition (EEMD) can overcome the mode mixing in certain degree, this limitation still exists in EEMD. Wigner-Ville distribution has a fine time-frequency resolution and can reach the low boundary of Heisenberg uncertainty principle. And Smoothed Pseudo Wigner-Ville distribution (SPWVD) can remove the cross-term in frequency direction and time direction of the time-frequency panel. The peak-track algorithm can extract the IF and IA of a monocomponent signal. For multicomponent signals, mashing operation should be conducted to guarantee that just one mode in time-frequency panel is processed by using the peak-track algorithm. To realize automatically mashing operation, a step based on the minimal and maximal values of a statistic parameter of the time-frequency distribution is presented to extract and separate valuable modes on the time-frequency panel. Combining with the peak-track algorithm and SPWVD, an algorithm that can compute the IF and IA of the GW signal is proposed. After that, a time-frequency analysis of debonding of a type of carbon fibre composites is conducted to discuss the feasibility of identifying the presence, size, and location of the defect.The rest of this paper is organized as follows. Section2 illustrates the principle of algorithm. Section 3 provides an illustration of the presented method and a contrast between it and several traditional time-frequency analysis techniques. Section 4 describes the detail of the experiment. Section 5 shows an analysis of the time-frequency analysis of the experimental signals employing the presented technique. Finally, Section 6 concludes the paper.


## 2. Principle of the Algorithm

The processing steps of the extension algorithm are shown in Figure1 and are as follows.Figure 1
The different steps of the algorithm presented in the paper.(1) Smoothed Pseudo Wigner-Ville distribution: Wigner-Ville distribution has a fine time-frequency resolution and can reach the low boundary of Heisenberg uncertainty principle. It is defined as follows [19]:(1)WVDst,f=∫-∞∞st+τ2st-τ2e-j2πfτdτ.SPWVD can suppress the cross-term in frequency and time directions. To get time-frequency panel with a promising time-frequency resolution, we employ SPWVD on testing signals. And the formula of SPWVD is as follows [20]:(2)SPWsn,Θ=∑k=-L+1L-1hk∑l=M+1M-1glsn+l+ks∗n+l-ke-j2kΘ,where g(l) and h(k) are smoothing window functions in time direction and frequency direction, respectively. s is an analyzed signal and n and Θ are time variable and frequency variable, respectively.(2) Extracting and separating valuable areas of the time-frequency panel. As mentioned above, significant time-frequency coefficients tend to cluster in groups or regions. To extract and separate them, we do following steps: (1) sum coefficients of each frequency-indexed group, and the calculation formula is(3)EΘ=∑n=1NSPWsn,Θ,where N denotes the number of the time variable of the time-frequency panel. (2) Calculate minimal values and maximum values of summations. Lines of the minimal values in the time-frequency panel can separate the valuable areas. And the maximum value evaluates whether the energy of a mode is strong enough to be kept. (3) Conduct threshold processing for the maximum values. And the equation is defined as follows:(4)EΘmax=keepEΘmax>T0discardEΘmax≤T0,where E(Θ)max means the maximum of E(Θ). T0 denote the threshold for sift valuable modes. A mode with a strong energy is more valuable for detecting the debonding defect, so a threshold process is introduced to improve calculation efficiency. Moreover, the noise can also be removed by this operation. After this step, the areas of the time-frequency panel corresponding to the significant modes can be extracted and separated.(3) Peak-track algorithm. Conduct the peak-track algorithm for the significant areas of SPWVD to obtain the primary IF and the primary IA of different models. The principle of the peak-track algorithm can be found in [21]. And then, the IFs and IAs of significant models can be obtained, which provides detailed information of the time-frequency distribution of the GWs for detecting defects. On the basis of these information, we can discuss the debonding of testing objects.(4) Threshold process for the primary IF and the primary IA: to suppress the influence of noise, a threshold process for the primary IA is conducted. The equation is defined as(5)A=A0A0>TA0A0≤TA,where A denote the final IAs of modes and A0 denote the primary IAs of modes. TA is a threshold to suppress the influence of noise. Correspondingly, the primary frequency is where the final IA is set to 0. And then, we can research the peak time tpeak of the final IA as follows:(6)tpeak=arg⁡max⁡IAt.In Section 3, the detail of the present algorithm will be illustrated with a sample signal.


## 3. Illustration of the Presented Algorithm

We construct a sample signal to illustrate the presented algorithm as follows:(7)st=c1s1t+c2s2t+c3s3t,s1t=0,0≤t≤0.003612503t-0.0036sin⁡75000×2πt,0.0036<t≤0.006-25009t-0.006sin⁡75000×2πt,0.006<t≤0.00960,0.0096<t≤0.012,s2t=0,0≤t≤0.002412503t-0.0024sin⁡100000×2πt,0.0024<t≤0.0048-25009t-0.0048sin⁡100000×2πt,0.0048<t≤0.00840,0.0084<t≤0.012,s3t=0,0≤t≤0.001212503t-0.0012sin⁡200000×2πt,0.0012<t≤0.0036-25009t-0.0036sin⁡200000×2πt,0.0036<t≤0.00720,0.0072<t≤0.012.The sampling frequency is 150 kHz. A noise with a strength of 0.1 is added to the signal to test the suitability of the algorithm. The sample signal contains three frequencies at 75, 100, and 200 kHz. In (7),   c1 and c2 are 1, and c3 is 0.4. The curve of the sample signal in time domain is shown in Figure 2.Figure 2
The curve of the sample signal in time domain.(  1) Smoothed Pseudo Wigner-Ville Distribution. We employ SPWVD for the sample signal to obtain the corresponding time-frequency panel, which is shown in Figure 3. A contrast among SPWVD, STFT, and EEMD is done. STFTs representation with long and short windows of the sample signal is shown in Figures 4 and 5. The STFT representation can be obtained by [22](8)STFTst,f=∫-∞∞sτh∗τ-te-j2πfτdτ,where h∗(τ-t) is the conjugate function of window function h(τ-t). It is different to keep a high frequency resolution and a high time resolution at the same time for time-frequency representation techniques, which is illustrated in Figures 4 and 5. When the window is long, the frequency resolution is good at the cost of the time resolution, as shown in Figure 4. By contrast, as can be seen in Figure 5, when the window is short, the time resolution is good, but the modes at 75 and 100 kHz cannot be separated.Figure 3
The SPWVD representation of the sample signal: (a) two-dimensional figure and (b) three-dimensional figure.
(a) (b)Figure 4
The STFT representation with a long window of the sample signal: (a) two-dimensional figure and (b) three-dimensional figure.
(a) (b)Figure 5
The STFT representation with a short window of the sample signal: (a) two-dimensional figure and (b) three-dimensional figure.
(a) (b)Wigner-Ville distribution has a fine time-frequency resolution and can reach the low boundary of Heisenberg uncertainty principle. However, for multicomponent signals, it suffers from an inevitable cross-term interference. SPWVD can remove it in frequency direction and time direction of the time-frequency panel. As can be seen in Figure3, the time resolution and frequency resolution of the time-frequency are promising. Moreover, no interference is in the representation.In EEMD, signals are decomposed as(9)st=∑i=1ncit+rnt,where ci(t) is the ith IMF and rn(t) is the residual signal which represents the slowly varying or constant trend of the signal. According to [23], an IMF with a high coefficient of correlation between it and the analyzed signal means that this IMF is an important component. Figure 6 presents the coefficients of correlation between different IMFs and the sample signal. We can learn that IMFs 1–3 are vital components of the signal, as the coefficients of correlation are relatively greater. So these IMFs are shown in Figure 7. It is visible that the mode mixing occurs in EEMD. Ref. [18] reveals that, when the ratio between a relatively low frequency and a relatively high frequency becomes close to 0.75 (low frequency/high frequency), two components of signals cannot be separated. Therefore, the poor frequency resolution is a limitation of EEMD for its application.Figure 6
Correlation coefficients between each IMF and the sample signal.Figure 7
The IMFs 1–3 in the EEMD of the sample signal.(  2) Extracting and Separating Valuable Areas of the Time-Frequency Panel. As mentioned above, peak-track algorithm can obtain the peak values and the corresponding coordinates in different time points of the time-frequency panel. However, when more than one mode exist on the panel, some errors may occur. To resolve it, we should extract and separate different modes from the time-frequency representation. At first, we sum coefficients of each frequency-indexed group by using (3). And the different frequency-group summations of the time-frequency panel of the sample signal are illustrated in Figure 8.Figure 8
The different frequency-group summations of the time-frequency panel of the sample signal.Afterward, we calculate minimal values and maximum values of summations. In the calculation process of the sample signal, 139 minimal values are obtained, which means that 138 modes will be extracted. According to (7), the sample signal includes 3 modes. The error root is from the noise. To remove the effectiveness of noises, we employ (4) to sift the valuable mode and set the threshold to 100 in (4). Three groups of minimal values are kept, meaning three modes are kept.Finally, we filter the SPWVD representation with the kept groups of minimal values. Figure9 shows the filtering result. As shown in Figure 9, three modes are extracted and separated.Figure 9
The filtering result of the SPWVD of the sample signal.
(a) (b) (c)(  3) Peak-Track Algorithm. We employ the peak-track algorithm on the filtering result of the SPWVD of the sample signal shown in Figure 9. The primary IAs and IFs are obtained and are shown in Figure 10. Because of noise, calculation errors occur in IFs.Figure 10
The primary IFs and IAs obtained by using the peak-track algorithm: (a) IF and (b) IA.
(a) (b)(  4) Threshold Process for the Primary IF. To suppress the influence of the noise, (5) is applied on the IA. The threshold TA is set to 0.1. The final IFs are got, as shown in Figure 11. As can be seen in Figure 11, the IFs of the three mode are 75, 100, and 200 kHz, respectively. The calculations of IFs are in a high accuracy. Moreover, the time resolution is also in a promising accuracy. For example, the peak times of the different modes are 0.004201 s, 0.005402 s, and 0.00006597 s, and the calculation errors are 0.00001 s, 0.00002 s, and 0.00003 s. These results prove the validity of the presented algorithm.Figure 11
The final IFs of the sample signal by using the presented algorithm.


## 4. Detail of the Experimental

The material of the specimen is T300/QY8911 composite materials. The size is 500 mm × 300 mm × 2 mm and contains 10 layers. The physical parameters of the material are as follows:E1=135 GPa, E2=8.8 GPa, E3=8.8 GPa, G12=4.47 GPa, G23=3.45 GPa, G13=4.47 GPa, ν12=0.3, ν23=0.34, ν13=0.3, and ρ=1560 kg/m3. The photograph of the experimental specimen is shown in Figure 12. The corresponding size diagram is presented in Figure 13. Two debonding defects are in the specimen. The side lengths of the defects are, respectively, 30 mm and 40 mm.Figure 12
The photograph of the specimen.Figure 13
The size diagram of the specimen.Figure14 is the diagram of the testing principle. In testing process, the exciting probe sends an exciting wave, and, then, the receiving probe will receive excited lamb waves. The excitation wave is a square wave with a frequency of 500 kHz. The sampling frequency is 25 MHz. Three groups of signal are collected in the experiment, responding to three situations, that is, no defect, 30 mm defect, and 40 mm defect, respectively. The corresponding GWs collected in the experiment are presented in Figure 15.Figure 14
The diagram of the testing principle.Figure 15
The GWs collected in the experiment: (a) no defect; (b) 30 mm; and (c) 40 mm.
(a) (b) (c)


## 5. Result and Discussion

Figures16–18 illustrate the final IFs and IAs of three groups of experimental signals by using the presented algorithm. First, two modes exist in the experimental signal of no defect and are at 50 kHz and 25 kHz. Modes at 50.8 kHz, 32.4 kHz, and 24.2 kHz are in the signal of 30 mm defect. Modes at 48.9 kHz, 35.6 kHz, and 24.3 kHz are in the signal of 40 mm defect. Considering that 50 kHz, 50.8 kHz, and 48.9 kHz pose little differences and the same for 25 kHz, 24.2 kHz, and 24.3 kHz, it seems that defects stimulate new modes, and 32.4 kHz corresponds to 30 mm defect, and 35.6 kHz corresponds to 40 mm defect. On the basis of this phenomenon, we can try to detect the defect. Moreover, the frequency of the new mode becomes greater along with the defect size (32.4 kHz for 30 mm, and 35.6 kHz for 40 mm). Maybe, we can do a try to evaluate the size of the defect, according to this relationship. The peak time of the mode becomes less along with the mode frequency increasing, which brings into correspondence with the conclusion that a higher frequency poses a higher propagation velocity. This phenomenon can provide a way to locate the defect. The discussions above reveal that, on the basis of the algorithm presented in this paper, we may detect the defect, calculate the size of the defect, and locate the position of the defect.Figure 16
The final IFs and IAs of no defect by using the presented algorithm: (a) IF and (b) IA.
(a) (b)Figure 17
The final IFs and IAs of 30 mm defect by using the presented algorithm: (a) IF and (b) IA.
(a) (b)Figure 18
The final IFs and IAs of 40 mm defect by using the presented algorithm: (a) IF and (b) IA.
(a) (b)


## 6. Conclusion

This paper presents a time-frequency algorithm aiming to analyze the characteristics of ultrasonic GWs generated in a NDT for the debonding in a type of composite material. The time resolution and frequency resolution of SPWVD are promising, which is advantageous to obtain the frequencies and peak times of different GW modes. On the basis of the presented algorithm, the characteristics of the experimental signals were investigated. The result shows that new modes can be stimulated by the debonding defect, and the frequencies of new modes also change along with the size of defect. Based on the presented algorithm, we can detect the presentence, calculate size, and locate the debonding defect. The technique also can be applied in analogue NDTs and NDEs on the basis of the ultrasonic GWs. Further research will be done to validate the feasibility for locating defects by the algorithm.



---

*Source: 5686984-2017-02-09.xml*