Abstract
Structural strength and integrity of composites can be considerably affected by the low-velocity impact damage due to the unique characteristics of composites, such as layering bonded by adhesive and the weakness to impact. For such damage, there is an urgent need to develop advanced nondestructive testing approaches. Despite the fact that the second harmonics could provide information sensitive to the structural health condition, the diminutive amplitude of the measured second-order harmonic guided wave still limits the applications of the second-harmonic generation–based nonlinear guided wave approach. Herein, laminated composites suffered from low-velocity impact are characterized by use of nonlinear guided waves. An enhancement in the signal-to-noise ratio for the measure of second harmonics is achieved by a phase-reversal method. Results obtained indicate a monotonic correlation between the impact-induced damage in composites and the relative acoustic nonlinear indicator of guided waves. The experimental finding in this study shows that the measure of second-order harmonic guided waves with a phase-reversal method can be a promising indicator to impact damage rendering in an improved and reliable manner.
Keywords
Introduction
Prevalence of composites is found in the industrial field of aircraft, transportation, and marine owing to their superior material properties. 1 However, since the composite laminates often lack the through-thickness reinforcement that can subject to impact in the transverse direction, such structure has low transverse damage tolerance. 2 For low-velocity impacts in the composite, laminated structures absorb the impact energy without plastic deformation.3–5 Defects existing in the inner layers of composite laminates caused by impact are complex and cannot easily be detected effectively. Recently, transverse stresses and plasticity propagation in layered structures were analyzed by Liu and Jeffers.6,7 Invisible interior defects caused by low-velocity impacts will degrade the structural performance of composite laminates.8,9 With the impact-induced damage affecting significantly the composite materials, failures to evaluate and characterize the low-velocity impacts will result in unpredictable catastrophes of the composite structure. Therefore, nondestructive testing (NDT) approaches need to be developed for the damage assessment of low-velocity impacts.
Numerous NDT techniques were used for the assessment of impact damage in composite laminates. 10 Impact-induced delamination in such materials was imaged by digital shearography, 11 which demands sensitive and expensive instruments as well as highly skilled inspectors for the use of this technique. In addition, the lack of clarity of defects deep beneath the surface also limits the application. Recently, a combination of the tomographic and optical techniques was used to evaluate the impact damage in composites by Sfarra and colleagues.12,13 Radiographic testing is another approach for damage inspection in composites. 14 However, the radiographic approach is only available in the case of delamination with specific orientation as regards X-ray beam. Ultrasonic testing is an effective method for impact damage detection in composites. However, most earlier investigations are about macro-damage detection by use of linear ultrasonic waves. 15 The nonlinear acoustic approach has been taken as a promising means for the assessment of micro-damage.16–18 Another investigation on impact damage detection was conducted in a sandwich plate by Meo and Zumpano. 19 Aymerich and Staszewski 20 investigated the relevance of impact energy and sidebands of ultrasonic amplitude by frequency modulation of ultrasonic waves in composite laminates. Progressive damage in polymer-based composites was monitored by nonlinear dynamics and the acoustic emission technique. 21
Recently, considering the great advantages of the guided wave technique for material characterization, many applications have been realized to detect and evaluate delamination, stress fatigue, thermal fatigue, and impact damage in composite materials by guided waves.22–26 Favorable characteristics of guided wave techniques comprise remarkable cost-effectiveness and low energy consumption.22,27 Recently, increasing attention has been paid to NDT in composites by nonlinear guided waves. Generation of higher harmonic waves is the typical ultrasonic nonlinear phenomenon. Second-harmonic generation (SHG) of guided waves in metallic materials was reported in earlier publications.28–31 An experimental scheme for characterizing the thermal damage was developed by SHG of Lamb waves in composites. 32 Second-harmonic guided waves were also adopted to investigate the potential of fatigue assessment in composites. 33 To the authors’ knowledge, SHG of guided waves is a sensitive indictor for composite structural health condition at an early stage.
Nevertheless, SHG of guided waves in composites can be readily omitted for the reason of a very small amplitude of measured second-harmonic signals. Accordingly, it is practically difficult to implement the nonlinear guided wave technique for the reason of low generation efficiency of second harmonics. In view of this difficulty, a phase-reversal approach is employed to intensify the signal clarity of SHG. Second-harmonic signals of a phase-matched guided wave are measured in this study. In addition, we extracted the second harmonics after the superposition of two propagating guided waves with the opposite phase. The method proposed is used for the damage assessment in the tested specimens. Variations of the acoustic nonlinear parameter corresponding to different artificial impact energies are obtained in the composite structures. Variations of wave velocity and acoustic nonlinearity with respect to impact energies are provided as well.
General considerations
The nonlinear feature of guided waves
The nonlinearity of materials is related to various imperfections in the material microstructure or discontinuities of interfaces.34,35 Distortion of ultrasonic waveforms occurs in tested materials with damage, and double-frequency second-harmonic waves are subsequently generated. The damage-induced acoustic nonlinearity (expressed as
where
As shown above,
Principle of the phase-reversal approach
The phase-reversal method is an effective means to counteract odd harmonics while highlighting the contribution of the even ones.36,37 The fundamental of the phase-reversal method presented in this study expands to only the second-order nonlinear term. The total wavefield regarding guided wave propagation in specimens with nonlinearity is
where
With the phase of the fundamental ultrasonic wave reversed 180°, while retaining the remaining conditions, the total wavefield in the waveguide is
By summing equations (3) and (4), one can see that only the second-harmonic field is left as shown below
It is shown that the phase-reversal technique doubles the second-harmonic amplitude while counteracting that of the primary ones. Apparently, the efficiency of SHG of primary guided wave propagation can be effectively enhanced by employing the phase-reversal method proposed here in nonlinear media.
Mode selection and signal excitation
Guided wave propagation in the composite structure has a few unique features due to anisotropic properties, such as direction-dependent wave speeds. In addition, many factors (e.g. thickness, density, and elastic properties) can affect the guided waves propagating in composites. Phase and group velocity dispersion curves play a critical role in mode selection and frequency tuning of guided wave NDT. In composite laminates, guided wave propagation can be described by satisfying Navier’s displacement equations within each layer. A transfer matrix, which relates the displacement and stress of each layer, can be constructed to represent the relations. In this study, the classical transfer matrix method is employed to calculate the dispersion curves of guided waves in the tested composite samples (shown in Figure 1).

Dispersion curves of (a) phase velocity and (b) group velocity. Wave propagation is along the fiber direction of the unidirectional composites.
As one can see from Figure 1, the selected phase-matched modes comprise the fundamental S1 mode waves (2.5 MHz) and the double-frequency S2 mode (5.0 MHz). The phase velocity (9.70 km/s) of the primary mode is equal to that of the second-order one, so does the group velocity (9.50 km/s).
In general, the dispersive nature of guided waves negatively affects the SHG. Thus, cumulative SHG versus propagation distance is crucial to measure the nonlinear response of guided waves. 30 The second harmonics of guided waves synchronize with the primary waves, which results in the cumulative phenomenon of SHG. More detailed discussions on this issue can be found in our earlier work. 32 In this investigation, S1 mode meets the condition of internal resonance (i.e. its phase and group velocity equal those of the double-frequency second-harmonic wave, as indicated in Figure 1), which ensures the cumulative effect of SHG.
In this study, excitation and reception of guided wave signals are realized by piezoelectric (PZT) transducers. The wave incident angle
where
Experiments
Specimens
The experimental tests were carried out in carbon fiber lay-up architectures, of which the size is 1.0 mm thick, 200 mm long, and 200 mm wide. The specimens used in this study are provided by the same supplier (Beijing Institute of Aeronautical Materials). Mass density and elastic constants of the tested specimen are provided in Table 1. The stiffness matrix is obtained by referencing the GB/T 1477 standard. One of the specimens is undamaged for the reference experiment. The low-velocity impact damage was induced by an MTS free-fall drop test machine (CEAST 9340; Instron), conforming to the ASTM D7136 standard. The impactor weighs 2.50 kg with a smooth hemispherical striker, of which the diameter is 16 mm. In this research, composite laminates were placed in the middle of the tested machine. The impact energies are calculated according to the falling height of the striker from the top of the crosshead at resting position on each sample. The impact energies of 5.0 and 15 J were derived from the gravitational potential energy equations referencing the ASTM D7136 standard.
Mass density
The impact-induced damage in the specimens was also inspected by the infrared thermography technique, as illustrated in Figure 2. It is shown that the impact damage can hardly be found in the specimens even by infrared thermography.

Infrared thermography of the impact damage in (a) virgin specimen, (b) 5 J impact loading specimen, and (c) 15 J impact loading specimen.
The system setup
An experimental setup is provided in Figure 3. A sinusoidal tone-burst electrical signal (2.50 MHz) is excited by a RITEC system, which then passes through a 14-dB attenuator and a low pass filter. A narrow band transmitter Tx (central frequency at 2.25 MHz) is employed for the generation of the primary wave. The other broad-band PZT transducer (central frequency at 5.0 MHz) is taken as a receiver Rx to measure the primary waves and the generated second harmonics concurrently. The transducers were clamped to the acrylic wedges. A light lubrication oil (type: Ultragel II; Sonotech) was used to facilitate the transmission of sound energy between different media. The transducers were mounted on the specimens with the couplant. The relative acoustic nonlinear parameters

Experimental setup.
Results and discussion
As analyzed in section “Mode selection and signal excitation,” the phase-matched mode pair selected is crucial for nonlinear ultrasonic guided wave testing. Thus, the group velocity of the propagating signal was measured and compared with the theoretical value. Besides, the frequency of the measured signal was checked to ensure that it is the expected guided wave mode. Practically, the A1 mode is very likely to be generated when the S1 mode is to be excited owing to their adjacent positions as shown in Figure 1(a). But it is also noticed that their group velocities differ from each other, so the wave package of the double-frequency harmonic signal (S1 mode) will eventually separate from the A1 mode wave package due to the noticeable group velocity difference. 32 The typical measured signals were modulated by Hanning window and occupy a length of 16 cycles, as illustrated in Figure 4(a). Figure 4(b) shows the frequency spectrum obtained by the FFT of measured signals. As one can see, group velocity (expressed as cg) of the received propagating signal is 9.43 km/s, and its central frequency is 2.50 MHz. Thus, the expected guided wave mode can be confirmed.

Confirmation of (a) group velocity and (b) frequency spectrum of the measured signal.
Figure 5(a) clearly presents the frequency spectra of primary guided waves and second harmonics. In addition, as indicated in equation (2), the relative nonlinear parameter (the ratio of measured amplitudes) versus wave propagation distance can be constructed to demonstrate nonlinearity. The cumulative effect mentioned previously was confirmed by Figure 5(b). In the experiment, the amplitude ratio of second harmonics to the second power of the primary one, which could generally be expressed as the relative acoustic nonlinearity as indicated in equation (2), is measured three times at each marked distance. The experimental average values with error bars are presented in Figure 5(b), which suggested a monotonic relationship with the propagation distance. Generally, SHG of guided wave propagation can be affected by coupling media, instruments, and material damage. It is noted that the effects of coupling media and instruments can be minimized by representing acoustic nonlinearity as the slope ratio of the measured

(a) Frequency spectra of the measured signal at a specific position and (b) relative nonlinear parameters measured at different propagation distances.
The effect of the phase-reversal approach was also examined. Here we show a specific implementation process. For the undamaged specimen, a sine pulse wave at the frequency of 2.50 MHz was applied to the incident transducer. All the settings are kept unchanged except for a phase-reversal step of the primary waves excited from the incident transducer and then the same measurement is repeated. Obviously, after superposition of these two phase-reversed signals, the signals with opposite phases can be counteracted, while signals in phase with each other will be enhanced as illustrated in section “Principle of the phase-reversal approach.”Figure 6 shows the spectra of measured signals in the frequency domain before and after the phase-reversal signal processing. The very small amplitude of the measured second harmonics is scarcely noticed in contrast to that of the primary wave under normal circumstances. However, the second-harmonic amplitude increases significantly with the phase-reversal technique adopted, while the two fundamental waves counteract each other. In detail, the second-harmonic signal (i.e. the double-frequency component) becomes twice as much as it is before. The remained primary wave amplitude does not exactly equal zero (see Figure 6(b)), which can be attributed to the inconsistent coupling conditions during the period of repeated measurements. Despite that, this technique is still able to effectively build up the second harmonics, which is favorable for the measurement of acoustic nonlinearity.

Frequency spectra of the measured guided wave signal by (a) the normal operation and (b) the phase-reversal method.
In the conventional approach, we calculate
As indicated in equation (5), the two fundamental waves offset each other while the amplitude of the double-frequency harmonic wave has doubled by the phase-reversal method. In this method, the two time-domain signals with the opposite phase are obtained simultaneously using excitations of two out-of-phase primary waves for a given propagation distance. Then, the primary wave amplitudes of the two signals in the frequency domain can be obtained, which are denoted by A1
a
and A1b, respectively. As regards the determination of
The experimental results indicate that the phase-reversal technique is a feasible method for the enhancement in the second-harmonic signal and thus increases the resolution for NDT evaluation of micro-damages in materials (due to higher signal-to-noise ratio (SNR) of the second-harmonic signals measured). SHG of guided waves is applied in detecting the low-velocity impact damage in composite laminates with the phase-reversal technique. Figure 7 provides us with the correlation of impact energies and the relative nonlinear parameter measured by the above experiment. As the impact energies increase, the relative acoustic nonlinear parameter increases significantly.

Relative nonlinear parameters measured in composites as a function of impact energies.
With regard to the sensitivity analysis of

Variations of
For low-velocity impacts in the composite, laminated structures absorb the impact energy without plastic deformation. Defects existing in the inner layers of composite laminates caused by impact are complex and cannot easily be detected effectively. Thus, matrix cracks, delamination, fiber breakage, as well as penetration can be initiated inside the composite structures. 38 The low-velocity impact-induced internal damage can be the source of material nonlinearity. Distortion of the ultrasonic waveform occurs in the tested material with damage, and double-frequency second-harmonic waves are subsequently generated. Thus, the measure of SHG is a convincing means to qualitatively assess the damage induced by low-velocity impact.
Conclusion
Low-velocity impact damage in composites is detected by the SHG of ultrasonic guided waves. The second harmonics with cumulative effect versus propagation distance are generated by the phase-matched guided wave mode pair. The applicability of the chosen mode is verified experimentally. An enhancement in the SNR for the measure of second harmonics is achieved by a phase-reversal method. The proposed method is applied to assess the different levels of low-velocity impact energies in the composites. An increase in the acoustic nonlinear parameter with respect to the impact energies is observed experimentally, while velocity variations can be neglected. This work suggests that the measure of SHG combined with the phase-reversal method is a potential technique for the assessment of low-velocity impact damage in composites.
Footnotes
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Natural Science Foundation of China (Grant Nos 11774295, 11834008, and 11474361).
