Impact Damage Detection in Laminate and Honeycomb CFRPs using Lamb Wave Ultrasonic Sensing

The paper presents the results on application of Lamb waves based technique for impact damage detection and severity identification. The PZT network operates in the round-robin mode changing the actuator and sensor roles of the transducers in order to detect the response of the system in the presence of damage. The monitoring is performed via the analysis of three parameters: change of the amplitude (dA), change of the energy (dP) and cross-correlation (NCC) of the signals in baseline and damaged state. Testing of laminate CFRPs shows that the damage location is estimated with an error of 5–15 mm, while the computed Damage index is linearly dependent on the applied impact energy. For honeycomb CFRPs the NCC parameter do not provide accurate results, however, other parameters allow identification within the 5–20 mm error and reflect accurate data on the severity of the damage.


INTRODUCTION
Non-destructive testing is a field of science and engineering that has great relevance for safety in many industries. But ensuring safety by NDT is often associated with a stop of operation and brings economic losses due to the unit downtime. So in industries demanding high quality inspection the Structural Health monitoring (SHM) concept becomes of great importance [1][2][3][4]. The development of SHM in aerospace is a serious engineering task that got a boost with the design of the aircrafts with wide CFRP usage (Boeing 787 and Airbus A350XWB) in main load-bearing structures: wings, fuselage, empennage, etc. The operation of such aircrafts has shown that the aircraft suffers large amount of external impacts of different origin: debris on the runaway, thunderstorms, bird strikes, low-quality or improper maintenance, etc. Such events lead to the formation of damages the vast majority of which are classified as barely visible impact damages (BVID) and do not lead to the urgent operation halt. The BVID of carbon fiber reinforced polymers are complex internal structure defects combining delamination, fiber breakage and pull out, and matrix debonding. Thus due to operational loading and environmental effect the BVIDs can evolve to large-scale weakened areas which are more difficult to repair or even they can become a threat to load bearing capacity of the structure. Thereby the monitoring systems that can provide timely data on the emergence of damages are of great importance in aerospace.
The scientists and engineers examine various physical principles to be utilized in SHM: vibration [5][6][7], ultrasonic waves [8], thermography [9,10], eddy currents [11,12], electromagnetism [13], etc. In the case of aerospace the SHM development task is more complicated due to hard limitations of the aircraft weight: the installed sensors, wires, and electronics can nullify the advantages of SHM application [14,15]. Thus opposite requirements for reduced weight and high safety promote using of the most high-tech solutions. For example, the aerospace considers acoustic emission [16][17][18][19][20], fiber-optic strain gauging [21][22][23], and other methods to be used in monitoring systems. Artificial neural networks are studied extensively to process raw data and provide efficient decision making. Recently, an increasing number of studies has been devoted to the investigation of ultrasonic Lamb waves in SHM [24,25]. Due to low attenuation these waves can ACOUSTIC METHODS cover extensive areas, that is extremely important in aircraft industry, where the area of CFRP skins can be very large. Lamb wave ultrasonic methods are widely studied for the detection of impact damages in composites [26,27]; the optimal networks of transducers are investigated [28][29][30]. However, the physical mechanism of Lamb waves propagation is quite complex [31,32] requiring consideration of many factors such as the stability of the adhesive layer and the degradation of its properties [33,34], temperature deviations [35,36], sensor failure [37], etc. as it can significantly affect the quality of obtained data and reduce the monitoring efficiency.
The technique for monitoring of impact damages studied in the paper was developed earlier and applied on the honeycomb specimen [38,39]. The algorithm was established using the artificial damages and subsequent multiple impact damaging of the specimen.
The present paper is a result of application of this technique for monitoring of different materials using modified testing procedure. The aim of the paper is investigation of the algorithms of detection of impact damages with different severity in laminate and honeycomb CFRPs using Lamb waves.

MATERIALS AND METHODS
There were two types of specimens used for testing of ultrasonic technique for the detection of barely visible impact damages. The first type is a laminate carbon fiber reinforced polymer with the following layup: [90 G /(0/-45/90/+45/0) 3 ] S . The specimens were produced with prepreg unidirectional carbon fiber fabrics using vacuum bagging and autoclave pressing (on the outer surfaces of the CFRP the glass fiber protective layer marked with "G" letter was placed). The total thickness of the 32-ply CFRP is 5.2 mm. The size of the specimen is 300 × 300 mm 2 . This specimen will be further denoted with "L" letter.
The second specimen is a sandwich with aluminum honeycomb filler and CFRP skins with the following layup: [90 G /0/-45/90/+45/0] S . Each 12-ply skin is about 1.7 mm thick and the total thickness of the specimen is about 18.2 mm. The size is the same as for the laminate specimens -300 × 300 mm 2 . They are denoted as "H" specimens below.
Each specimen was instrumented with two networks (on the front surface to be impacted and on the rear surface) of PZT transducers allowing both actuating and sensing. The network had a size of 4 × 4 with a 45 mm offset from the specimen boundary and a step of 70 mm. The sketch depicting the PZT network is shown in Fig. 1. The PZTs transducers used are AW1E12G-190EFL1Z by Audiowell with the diameter of 9 mm and thickness of ceramic layer of 0.19 mm bonded to the steel substrate. The use of two similar PZT networks on both surfaces allow us to estimate the possibility to detect damage of the outer surface of structure from the sensing system installed inside and protected from the environmental conditions. This is quite important for honeycomb CFRP which has higher attenuation of the ultrasonic waves traveling from one skin to another by aluminum honeycomb core. 3M DP490 epoxy adhesive was used to bond the transducers on the surface of CFRP specimens.

Offset
Instrumented with PZT Step Y X The waveform used is 5 cycles of sine modulated by Hanning window. The frequency distribution of the waveform results in moderate dispersion during wave propagation compared with sine or pulse signal. This makes the analysis of the obtained results easier because the reduced dispersion allows classifying the components of the received waveform, e.g. direct passed wave and edge reflections. The actuating and sensing hardware are arbitrary waveform generator AWG-4105 and 4-channel USB oscilloscope Handyscope HS4 correspondingly. AWG-4105 allows generating waveforms with amplitude of 10 V which was enough for the experimental testing of the ultrasound technique for 300 × 300 mm 2 CFRP specimens. However the magnitude of sensed signal for honeycomb specimens was insufficient, thus the high-frequency amplifier was used for oscilloscope.
During the preliminary tests of "actuator-sensor" pairs bonded on the surface of CFRP specimens the frequencies of 100 and 300 kHz were chosen. According to the literature data of Lamb waves dispersion curves and simulation of ultrasonic waves propagation in LAMSS-Composites software [40] it was established that 100 kHz frequency corresponds to the A0 mode with a wavelength of 18 mm while the S0 mode is negligible. At 300 kHz S0 is generated much better (9 mm diameter transducer is comparable to the half of wavelength) instead. However A0 is also presented with much lower amplitude. The A0 mode wavelength at 100 kHz is quite large compared to the size of the BVIDs but it was expected that the waves would interact with such BVIDs with increased attenuation. This phenomenon is due to decrease of stiffness and emergence of defects like delamination, matrix cracking, fiber breakage, etc.
It should be noted that literature review [41] and the tests conducted earlier by the authors show that the difference in wave propagation velocity in 0 and 90 fiber direction differs by 1-3% thus the anisotropy of the specimens with quasi-isotropic layup can be neglected.
In order to increase the S/N ratio the received signals were averaged 100 times. The signal acquisition by HS4 oscilloscope was triggered by the reference signal from the generator using the first channel of the oscilloscope. Three remaining channels were used for data capturing. Each position of the generator (one of the 12 transducers) requires the data from 11 sensors to be captured thus the relay box was designed for sequential registering the data from all sensors. The relay box allows changing the generator position simplifying the experimental procedure. Summary there are 132 signals for the 4 × 4 pattern of PZTs to be captured. Using of two frequencies and testing of two sides resulted in the set of 528 waveforms for ultrasonic testing of one specimen.
The proposed technique was evaluated for possibility of detection of damage position and severity of impact of varying energies. The specimens were impacted according to the ASTM D7136 standard using a drop weight technique. The 2 kg impactor having a hemispherical striker tip with a diameter of 16 mm was used. The impact energies were 10, 15 and 20 Joules for laminates while honeycomb CFRPs were impacted with 2 and 4 Joules.
Before the impact testing the ultrasonic monitoring was performed creating the baseline set of signals. After the impacting of the specimens, the same procedure was performed resulting in the set of signals for the damaged state. The mathematical comparison of these sets of waveforms allows to detect the presence of a defect and its location as well as to estimate roughly the severity of the damage.
Numerical evaluation of the ultrasonic signals in order to detect the changes occurring due to impact damaging of the CFRP was performed by calculation of the set of informative parameters (between the signals for two states-damaged and baseline): difference of maximums of envelopes, normalized correlation coefficient and energy difference.
The calculation of the signal envelope of the received signal is performed using Hilbert transform in the frequency domain. On the first step the n-point Discrete Fourier Transform (DFT) of the N-length signal is calculated and spectrum values X(m) corresponding to the negative frequency are zeroed. The values X(0) and X(N/2) are divided by 2. Then the reverse N-point DFT is performed. The signal has real and imaginary parts. The values of envelopes are determined as the absolute values of the obtained signals. The difference of envelopes maximums dA is calculated as follows: where А bas is an envelope maximum for the baseline state signal and A dam is an envelope maximum for the damaged state signal.
Normalized correlation coefficient if obtained by the following equation: where x bas and x dam are the signal magnitudes of the baseline and damaged states correspondingly.
In order to obtain the power spectrum the DFT was performed and the values of the spectrum were calculated by the expression: (3) where X are values of the Fourier spectrum. Then for each signal the energy within the range of 0.4f was calculated: (4) where f is a point corresponding to the central frequency of the signal. Then for two signals the difference dP is calculated: where P bas is an envelope maximum of the baseline state signal and P dam is an envelope maximum of the damaged state signal.
After the testing the database for each specimen contains the following data: locations (coordinates) of PZT transducers, sets of registered signals for all actuator-sensor pairs, informative parameters calculated by comparison of baseline and damaged state data. Then using developed software the data is processed and analyzed in order to obtain the location and severity of impact damage. The algorithm of processing is described below: (1) the wave paths are sorted according to the ascending of absolute relative change of dA, dP or NCC value in assumption of higher attenuation due to presence of damage; (2) in the sorted list the first N paths are chosen (in the present paper 10 and 20 paths were used for comparison); (3) for each pair of wave paths the coordinates of intersection points are obtained; (4) then according to the equations (6) and (7) the coordinates of the position of the damage and its severity are calculated.
The coordinates of impact damage are calculated using equation: where r i are coordinates of i-th intersection point of wave paths, w i are values of chosen informative parameters. The damage severity is assessed using dimensionless damage index (DI): where w i are values of chosen informative parameters, n is a number of wave paths. The higher the DI the more severe the damage is. Both coordinates of the damage and DI can be obtained using each informative parameter (dA, dP, or NCC).

RESULTS AND DISCUSSION
The section presents the experimental results obtained after investigation. Table 1 presents the data on studied specimens. The laminates and honeycombs are designated as L and H correspondingly.
In order to describe operation of the experimental technique the detailed results for the specimen L2/20 are provided. Figure 2 shows the 4 × 4 sensor network with a generator defined in lower corner: possible wavepaths and the differences of the maximums of envelopes between baseline and damaged states are shown. This example corresponds to the dA equal to 0.41, i.e. the drop of amplitude by ~41% for "actuatorsensor" pair 10→5 of the L2/20 specimen.
It can be seen that the neighboring sensors 4 and 6 show drops of 12 and 23% correspondingly. The change of the signals in presence of damage for other sensors is quite low. The example of initial waveform with a corresponding envelope for L2/20 specimen is shown in the upper plot in Fig. 3. Lower plot describes the envelopes for both baseline and damaged state and the difference between maximums is graphically shown.  (6). This dot array is calculated for each combination of testing frequency, side of the CFRP specimen, amount of calculation paths and informative parameter (red for NCC; green for dA; and blue for dP). Thus there are 2 3 = 8 damage location prediction points for each informative parameter. These points are used to obtain the resulting X and Y coordinates which are plotted as crossed circles of the corresponding color.
Quantitative results for all laminate specimens are summarized in Table 2. It can be seen that the technique allows detecting the location damage with an errors of 5-15 mm. Monitoring of the impacted side and back surface provides the same results about impact damage prognosis. The scatter of Δ avg for any parameter is nearly the same allowing to conclude that all of them can be utilized in damage location procedure. The last column in the table shows the Damage Indexes calculated by the equation (7) and averaged for each informative parameter. These DI have different ranges of variation. DI dA is the simplest to analyze varying from 0 to 1, where 0 stands for an absolutely non-damaged state resulted from the comparison of baseline and damaged ultrasonic data with no amplitudes changes while 1 is obtained if the all damaged signals will be nullified. Each number between 0 and 1 describes the average drop of amplitude of signals used for damage location. The same idea is for DI dP except that the energy of the signals is calculated and compared instead of amplitude. Both these parameters are aimed to detect higher attenuation   of the Lamb waves in presence of damage. DI NCC varies in the range of 0 to 1 as well, but demonstrates different behavior because the correlation can change due to different reasons, like phase shift. The 0 value stands for absolutely non-damaged state while 1 for damaged but the dependence within this range is highly nonlinear. Thus DI NCC for the damages obtained due to 15-20 Joules impacts reach 0.8-0.9 but the BVIDs obtained are quite small. Thus it can be concluded that NCC will have very low classification ability for larger damages. Figure 6 presents DI dA for all laminate specimens are presented. The data are fitted using linear function resulting in good linear dependence of the damage index versus impact energy.
Then Lamb wave ultrasonic testing of all honeycomb specimens has been performed. The preparation and testing procedure for honeycomb specimens was absolutely the same as one utilized for laminates but the results showed some differences. Firstly it has been found that damage location procedure using NCC parameter provides unreliable results for any H specimen. The calculation of the damage prediction points for every combination of frequency, side of the specimen and number of wavepaths provides high scattered set of points which after averaging resulted in X and Y which are very close to the center of the specimen (150; 150). Above said can be seen in Figs. 7 and 8. For both specimens H1/2 and H2/4, as well as for remaining H3-H6, the scatter of red dots is high and the averaged damage prediction is nearly in the center.
Secondly the testing of honeycomb specimens showed that damage prediction points obtained by the network on the rear side (opposing the side of the impact) demonstrate high scatter for dA and dP informative parameters. This can be seen for nearly all honeycomb specimens and it should be concluded that the detection of the damage on the side opposite to external impacted surface is not possible when using such technique. The detailed results of high scatter of rear side network location are shown in Table 3 while for laminates such data has not been provided due to similarity of the numbers for front and rear PZT networks.   Thus for the calculation of averaged X, Y coordinates for dA and dP parameters only the front PZT networks has been used. All the results are summarized in Table 4. It can be seen that the damage location prediction error for dA varies within the range of 10-20 mm while dP has higher scatter of ∆ avg : 4-35 mm.
Damage indexes for various parameters have been calculated showing the sensitivity and classification ability as well as for laminate specimens (Fig. 9). The scatter of the DI dA is higher than for laminates but the direct dependence is clearly distinguishable.

CONCLUSIONS
The experimental technique utilizing Lamb waves for BVID detection of CFRP has been studied. Piezotransducers adhesively bonded to the surface were used to actuate and detect the Lamb waves. The round robin testing of laminate specimens shows that all informative parameters allow detection of damage location prediction within the range of 5-15 mm. The damage index shoes linear dependence on impact energy. Finally both front and back surfaces can be used for PZT system having a little effect on accuracy. NCC can be used to detect the damage but correct evaluation of severity degree using this parameter is complicated.
Testing of honeycomb specimens showed sufficient prognosis accuracy in the range of 5-20 mm in spite of the larger scatter than for laminates. The monitoring of damages by the PZT on the side that is opposite to the impacted one is more complicated and the accuracy is much less due to the presence of aluminum honeycombs and their large thickness. Thus the Lamb waves are insensitive to damages on the opposite side. Damage Index allows to classify the damage severity. NCC parameter cannot be used for honeycomb CFRP because of the influence of phase shift on the NCC value.

FUNDING
The work was performed according to the government research assignment for ISPMS SB RAS, project FWRW-2021-0010.
OPEN ACCESS This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.