# Delamination Localization in Sandwich Skin Using Lamb Waves by Finite Element Method

Authors: Salah Nissabouri; Mhammed El allami; El Hassan Boutyour; Ahmed Errkik
Journal: Advances in Acoustics and Vibration
 (2018)

---

## Abstract

In this work we model by finite element method (FEM) the Lamb waves’ propagation and their interactions with symmetric and asymmetric delamination in sandwich skin. The simulations were carried out using ABAQUS CAE by exciting the fundamental A0 Lamb mode in the frequency 300 kHz. The delamination was then estimated by analysing the signal picked up at two sensors using two technics: Two-Dimensional Fast Fourier Transform (2D-FFT) to identify the propagating and converted modes, and wavelet transform (WT) to measure the arrival times. The results showed that the mode A0 is sensible to symmetric and asymmetric delamination. Besides, based on signal changes with the delamination edges, a localization method is proposed to estimate the position and the length of the delamination. In the last section an experimental FEM verification is provided to validate the proposed method.


---

## Body

## 1. Introduction

A sandwich composite is a materiel made from two thin skins bonded to a thick core. The skin studied in this paper is an orthotropic plate[0]4 with three mutually perpendicular planes of symmetry. The sandwich materials are designed to improve mechanical proprieties of structures. They are widely used in different fields especially in aeronautic industry. Their proprieties are influenced by the proportions of the matrix and the reinforcements. There are other parameters that also affect their characteristics like size, orientation, and distribution of the fibre. However, the heterogeneity of the composites structures leads to their weakness and facilitates the appearance of internal and external damage such as fibre breakage, matrix cracking, through-thickness hole, local delamination. Among these types of damage, delamination is especially easy to appear because the transverse tensile and interlaminar shear strengths are weak compared to the in-plane strength [1]. The delamination causes wave scattering, mode conversion, and multiple reflections. To understand these mechanisms, theoretical, numerical, and experimental studies are conducted.Feng, Lopes Ribeiro, and Geirinhas Ramos [2] analysed the interaction of symmetric S0 and antisymmetric A0 mode with the delamination using finite element simulations. The Lamb waves propagation in a 4-layer [0/90]s laminate were compared with propagations obtained in 1-layer [0] and in 3-layer [90/90/0] sublaminates. Chiua, Roseb, and Nadarajaha [3] investigated the scattering of the S0 mode by a delamination in quasi-isotropic fibre-composite laminate. Guo and Cawley [4] studied by finite element analysis and by experiment the interaction of the S0 Lamb mode with delamination. Nadarajah, Vien, and Chiu [5] presented results for the scatter field for various angles of incidence and for varying defect sizes. Hayashi and Kawashima [6] studied the reflections of Lamb waves at a delamination by semianalytical finite element method. Ching-Tai and Veidta [7] investigated the scattering characteristics of A0 mode Lamb wave at a delamination in a quasi-isotropic composite laminate. Bin [8] suggested an algorithm to localize and identify the damage in Woven Glass Fibre reinforced epoxy (WGF/epoxy). Mustapha [9] characterized fundamental symmetric and antisymmetric Lamb modes in terms of their velocity and magnitude variation as they change gradually in the thickness of a composite sandwich plate with a high density foam core. Ng [10] presented a theoretical and finite element (FE) investigation of the scattering characteristics of A0 at delaminations in a quasi-isotropic composite laminate. Veidt and Ng [11] studied the influence of stacking sequence on fundamental antisymmetric Lamb wave (A0) scattering characteristics through holes in composite laminates. Luca [12] developed a valid finite element model to simulate Lamb waves’ propagation in a Carbon Fibre Reinforced Plastic (CFRP) laminate for damage detection purpose and investigated the effects of the wave interaction with respect to damage parameters such as size and orientation. Yang [13] investigated some aspects of numerical simulation of excitation and detection of Lamb waves using piezoelectric disks in plate-like composite laminates.In this paper, a numerical approach is proposed to localize a delamination in sandwich skin as presented in Figure1. The aim of this work is to study the conversion by identifying the propagating modes and also to evaluate the symmetric and asymmetric delamination.Figure 1
Sandwich structure with delaminated skin.


## 2. Lamb Waves Theory

Lamb waves are most used for many reasons; they can propagate long distances alongplates and shells so they permit quick inspection of large structures, and also they are sensitive to the small variations either in material proprieties or in structure of the plate. However, they are dispersive which means that the interpretation of received signals can be complicated. So the key is to choose one pure mode to excite and to analyse its reflection, conversion, and transmission.


### 2.1. Equation of Lamb Waves Propagation in Unidirectional Lamina

The characteristic equations of the symmetrical and antisymmetric waves are [14](1)D11D23cotankα1h2-D13D21cotankα3h2=0D11D23tan⁡kα1h2-D13D21tan⁡kα3h2=0with Dr coefficients calculated by the following equations:(2)D1r=C13+C33αrWrD2r=C55αr+WrWr=ρc2-C11-C55αr2/(C13+C55)αr. For r {1,2,3,and4}, with(3)α1=-α2=-B-B2-4AC2Aα3=-α4=-B+B2-4AC2Aα5=-α6=ρc2-C66C44coefficients A, B and C are calculated by the following equations:(4)A=C33C55B=C33C11-ρc2+C55C55-ρc2-C13+C552C=C11-ρc2C55-ρc2where h is the skin thickness, k is the wavenumber corresponding to the X direction, c is the phase velocity. C11, C13, C33, C44, C55, C66 are the elastic constants, and ρ is the skin density.


### 2.2. Dispersion Curves

The dispersion curves presented in Figure2 for orthotropic plate [0]4 are plotted using DISPERSE program [15] in terms of wavenumber versus frequency.Figure 2
Dispersion curves (wavenumber versus frequency).At low frequency, two modes, A0 and S0, can propagate. As the frequency increases, more modes are possible to propagate.


### 2.3. Lamb Waves Control Technics

Among the contact technics to generate and to receive the signal, we find two methods: pitch catch and pulse echo. These two technics are simple but they require the coupling medium which limits the transmission of acoustic power.In the pulse-echo method, the reflected signal permits evaluating the damage. In the pitch catch technic the transducer placed in the transmission side permits evaluating the defect. In this paper two methods have been used to evaluate the delamination.


## 2.1. Equation of Lamb Waves Propagation in Unidirectional Lamina

The characteristic equations of the symmetrical and antisymmetric waves are [14](1)D11D23cotankα1h2-D13D21cotankα3h2=0D11D23tan⁡kα1h2-D13D21tan⁡kα3h2=0with Dr coefficients calculated by the following equations:(2)D1r=C13+C33αrWrD2r=C55αr+WrWr=ρc2-C11-C55αr2/(C13+C55)αr. For r {1,2,3,and4}, with(3)α1=-α2=-B-B2-4AC2Aα3=-α4=-B+B2-4AC2Aα5=-α6=ρc2-C66C44coefficients A, B and C are calculated by the following equations:(4)A=C33C55B=C33C11-ρc2+C55C55-ρc2-C13+C552C=C11-ρc2C55-ρc2where h is the skin thickness, k is the wavenumber corresponding to the X direction, c is the phase velocity. C11, C13, C33, C44, C55, C66 are the elastic constants, and ρ is the skin density.


## 2.2. Dispersion Curves

The dispersion curves presented in Figure2 for orthotropic plate [0]4 are plotted using DISPERSE program [15] in terms of wavenumber versus frequency.Figure 2
Dispersion curves (wavenumber versus frequency).At low frequency, two modes, A0 and S0, can propagate. As the frequency increases, more modes are possible to propagate.


## 2.3. Lamb Waves Control Technics

Among the contact technics to generate and to receive the signal, we find two methods: pitch catch and pulse echo. These two technics are simple but they require the coupling medium which limits the transmission of acoustic power.In the pulse-echo method, the reflected signal permits evaluating the damage. In the pitch catch technic the transducer placed in the transmission side permits evaluating the defect. In this paper two methods have been used to evaluate the delamination.


## 3. Postprocess Tools

### 3.1. 2D-FFT Analysis

Lamb waves’ propagation is sinusoidal in the frequency and spatial domains. For that reason, a temporal Fourier Transform is applied to go from the time to frequency domain; after that, a spatial Fourier Transform is computed to obtain the frequency wavenumber domain; see [16]. In practice, carrying out a spatial Fourier method to data obtained experimentally or numerically requires applying a 2D-FFT, using the following.(5)Hk,f=∬-∞+∞ux,te-ikx+ωtdxdt


### 3.2. Wavelet Analysis

Wavelet transform is an important tool in the time frequency domain of transient signals. The continuous wavelet transform of signal u(t) is defined as follows.(6)cwta,b=1a∫-∞+∞utψ∗t-badtThe obtained signal is a function of a and b, the translation and scale parameters, respectively. The parameter psi (t) is the transforming function named the mother wavelet. In this paper we use wavelet transform by “gaus1” to locate the peak which permits determining the arrival time of the wave at specific frequency. Moreover, for each frequency f, the localization of maximal value of the wavelet coefficients cwt(a, b) allows identifying the arrival times t1 and t2. Knowing the group velocity and the delay (delta) Δt =t2-t1 between the two modes with the same nature, we can calculate the distance by the following equation.(7)Δx=VgΔtKnowing the delay between two different modes (symmetric and antisymmetric), we calculate the distance by the following equation:(8)Δx=ΔVgΔtwith (9)ΔVg=VSVAVS-VA.Here VS and VA are the group velocities, respectively, of S0 and A0 modes. Arrival time method was used by many authors to localize the delamination; see [17].


## 3.1. 2D-FFT Analysis

Lamb waves’ propagation is sinusoidal in the frequency and spatial domains. For that reason, a temporal Fourier Transform is applied to go from the time to frequency domain; after that, a spatial Fourier Transform is computed to obtain the frequency wavenumber domain; see [16]. In practice, carrying out a spatial Fourier method to data obtained experimentally or numerically requires applying a 2D-FFT, using the following.(5)Hk,f=∬-∞+∞ux,te-ikx+ωtdxdt


## 3.2. Wavelet Analysis

Wavelet transform is an important tool in the time frequency domain of transient signals. The continuous wavelet transform of signal u(t) is defined as follows.(6)cwta,b=1a∫-∞+∞utψ∗t-badtThe obtained signal is a function of a and b, the translation and scale parameters, respectively. The parameter psi (t) is the transforming function named the mother wavelet. In this paper we use wavelet transform by “gaus1” to locate the peak which permits determining the arrival time of the wave at specific frequency. Moreover, for each frequency f, the localization of maximal value of the wavelet coefficients cwt(a, b) allows identifying the arrival times t1 and t2. Knowing the group velocity and the delay (delta) Δt =t2-t1 between the two modes with the same nature, we can calculate the distance by the following equation.(7)Δx=VgΔtKnowing the delay between two different modes (symmetric and antisymmetric), we calculate the distance by the following equation:(8)Δx=ΔVgΔtwith (9)ΔVg=VSVAVS-VA.Here VS and VA are the group velocities, respectively, of S0 and A0 modes. Arrival time method was used by many authors to localize the delamination; see [17].


## 4. Numerical Model Using FEM

Numerical simulations were carried out using finite element software ABAQUS to predict the A0 Lamb wave propagation behaviour in a [0]4 plate. A number of modelling methods had been used by previous researchers to model the propagating wave: finite difference, boundary elements, strip elements, hybrid and finite element. The finite element method has been extensively and successfully used to model the Lamb waves and study their interactions with defects in structures [18]. Many researchers have used the numerical simulation to study Lamb waves. Chiua, Roseb, and Nadarajaha [3] built a 3D model in ABAQUS of the 8-ply laminate [45/-45/0/90]S. Guo and Cawley [4] investigated the influence of delamination position through the thickness and its interaction with the S0 mode using finite element analysis. The composite laminate modelled was the eight-layer cross-ply. Nadarajah, Vien, and Chiu [5] studied the scattering of a zero-order antisymmetric (A0) Lamb wave mode by semicircular mid-plane edge delamination using the commercial FE package ABAQUS. Hayashi and Kawashima [6] studied Lamb waves’ propagation in laminated plates with delamination using the Strip Element Method and discussed the reflection and transmission for the case S0 and A0 excitation. Ng [10] modelled the delamination as a volume split in quasi-isotropic plate. The FE results were carried out by the explicit FE code LS-DYNA. Ng and Veidta [7] used a three-dimensional FE method to simulate an eight-ply [45/-45/0/90]S quasi-isotropic composite laminate with a delamination. The simulations were computed by ANSYS software. Panda [19] presented 3D FE simulations that were carried out to visualize the wave propagation and their interaction with the defect at various depths of an 8-layered Glass Fibre Reinforced Polymer (GFRP) by finite element software ABAQUS/Explicit. Ramadas [20] investigated the interaction of the primary antisymmetric mode A0 with symmetric delamination type defects in a quasi-isotropic laminated composite using 2D model ANSYS. Soleimanpour and Ng [21] used 3D explicit finite element method to study cross-ply laminated composite beams [0/90/0/90]S; they investigated the mode conversion and scattering characteristics of guided waves at delamination. Gudimetla and Kharidi [22] showed a procedure to simulate the propagation of Lamb waves in 8-layered Carbon Reinforced Fibre Plastic (CRFP) using a 2D model in ANSYS. The authors simulate the propagating Lamb waves along the plane of the structure in the form of a time dependent force excitation. Basri and Chiu [18] investigated how Lamb waves respond to the presence of material degradation in a plate-like structure using a series of finite element analyses. The propagation of these guided waves was interpreted with the dispersion characteristics, and displacement profiles were analysed in the frequency and wavenumber domain.In this paper we try to predict the interaction of A0 mode in orthotropic plate. The 2D numerical simulations were carried out using ABAQUS CAE. We consider the Lamb waves propagation only in the skin to simplify the propagating modes identification. This assumption has been made in many papers [23–25]. They consider only the propagation in skin as its acoustic impedance is more than the acoustic impedance of the core.The skin has a length of L=400 mm and a thickness of h = 1.6 mm. The mechanical properties of each lamina are shown in Table1. Once the geometry of the plate has been achieved, it remains to mesh and to define a sampling sufficient time.Table 1
Elastic coefficientsCij and density ρ of the materiel.
ρ (kg/m3) C 11(GPa) C 22 = C 33 C 12 = C 13 = C 23 C 55 = C 66  1530 56,9 14,7 9,76 4,16


### 4.1. Meshing and Time Sampling

To satisfy an accurate solution, the model has been meshed using (10), and the time step is calculated by (11).(10)max⁡Δx,Δy<λmin10(11)Δt<0.7max⁡Δx,ΔyVL


### 4.2. Excitation

In order to identify accurately the propagating modes, the simulation is performed on the range of frequencies to prevent dispersive phenomenon. For that reason we choose 300 kHz as frequency of excitation (see Figure3).Figure 3
Excitation windowed by Hann window.


## 4.1. Meshing and Time Sampling

To satisfy an accurate solution, the model has been meshed using (10), and the time step is calculated by (11).(10)max⁡Δx,Δy<λmin10(11)Δt<0.7max⁡Δx,ΔyVL


## 4.2. Excitation

In order to identify accurately the propagating modes, the simulation is performed on the range of frequencies to prevent dispersive phenomenon. For that reason we choose 300 kHz as frequency of excitation (see Figure3).Figure 3
Excitation windowed by Hann window.


## 5. Localization Method

The finite element method with ABAQUS CAE code is used to compute the temporal displacement by 2D-FFT technic (5) so as to identify the propagating and converted modes. Once the modes are identified, the next step consists in calculating the arrival times by WT (6). So there are two cases: in the first case, where there is no modes conversion, the time delay permits calculating the delamination position by (7); in the second case, where the mode converts to another (A0 to S0 or S0 to A0), (8) is used.


## 6. Results and Discussion

### 6.1. Model:04 without Defect

The A0 mode was excited from the left edge of the plate by concentrated force at the actuator C (see Figure 4). The A0 Lamb wave is the focus in this study because of its sensitivity to the small defects. To verify the excitability of pure mode A0, 2D-FFT was performed at the sensor A. A matrix was calculated at uniformly spaced points.Figure 4
Numerical model without defect, C: actuator, A and B: sensors.Figure5 shows that there is only one mode propagating which is A0.Figure 5
Superposition of 2D-FFT over analytical dispersion curves.


### 6.2. Model: Symmetric Delamination 60 x 0.1

Figure6 is the delaminated model. Figure 7 is a snapshot of FE out-of plane displacement. The Lamb waves split into two parts and travel independently with different velocities in the upper and the lower sublaminates. Mode conversions happen at the two edges of a delamination due to the change of boundary conditions. To characterize the delamination, we analysed the signals picked up at the sensors A and B.Figure 6
Numerical model of symmetric delamination.Figure 7
Snapshot of FE simulated out-of plane displacement in the upper and the lower sublaminates separated by the delamination.


### 6.3. Modes Conversion

In Figure8, the first wave group is a forward traveling wave group, which is an excited A0 mode atx= 0. When this interacted with the entrance of the delamination, there was a low amplitude reflected wave group. This reflected wave group appears as the second wave group A0A0. Delamination separates the plate into two sublaminates: upper and lower. The transmitted wave converted to the symmetric and antisymmetric mode A0S0 and A0A0 which traveled within the sublaminates with different velocities. The first wave A0S0U—U, L mean that the mode travel, respectively, in the upper and lower sublaminate—reached first the exit of delamination; this is justified by the fact that symmetric mode is faster than the antisymmetric mode. The backward mode A0S0US0U passes through the upper and then is transmitted to the main laminate as the second wave group A0S0US0UA0. By the same manner, in the lower sublaminate the mode A0S0LS0L is transmitted to the main laminate as A0S0LS0LA0. In the case of symmetric delamination the two modes A0S0US0UA0 and A0S0LS0LA0 are reflected at the same time.Figure 8
A0 reflection at the entrance and at the end of asymmetric delamination (length 60), case A0 incident; R denotes reflected mode and T denotes transmitted mode.


### 6.4. Results: Sensor A

Figure9 is A-scan image from the numerical simulations, taken at 100 mm for a delamination length of 60 mm. As explained in Figure 8, we conclude that the first arrived wave is the incident mode A0 and the second is the first reflected wave A0A0. The delamination is symmetric so two waves arrived at the same time as one packet which is A0S0US0UA0+ A0S0LS0LA0.Figure 9
Out-of-plane displacement, taken on damaged composite[0]4 at x = 100 mm, case symmetric delamination length 60.Based on the 2D-FFT computed using Matlab program applied to the surface of delaminated skin (before the delamination), Figure10 shows that the incident and the reflected modes are antisymmetric which confirms that the first reflected mode is A0A0.Figure 10
Superposition of 2D-FFT applied before the delamination over analytical dispersion curves.


### 6.5. Characterization of Delamination from Data Acquired at Sensor A

Figure11 presents wavelet coefficients versus time at sensor A. The time delay between the peaks permits estimating the position and the length of delamination.Figure 11
Coefficients of wavelet gaus1, case A0 incident, delamination length 60, between layer 2 and layer 3.Knowing the group velocity of the modes A0 and S0: Vg (A0,300kHz)=1268 m/s, Vg (S0,300kHz)=5649 m/s, we can calculate the delamination position from the sensor A noted PD by Δt1 (Delta 1) and delamination length noted LD by Δt2 (Delta 2) as presented in Table 2.Table 2
Delamination localisation from data acquired at sensor A.
Delamination Position Length P D P D T H Error % L D Error % 1 and 2 60 103,8 100 3,8 62,0 3,3 90 103,9 100 3,9 93,3 3,7 120 103,8 100 3,8 123,4 2,9 2 and 3 60 103,2 100 3,2 62,2 3,7 90 106,6 100 6,6 94,1 4,5 120 103,1 100 3,1 120,9 0,8The sensor A permits evaluating the position and the length of symmetric and asymmetric delamination. The errors associated with the predicted symmetric delamination position range from 3.1% to 6.6%. However in the case of asymmetric delamination, the error is less than 3.8%. As far as the length prediction is concerned, the errors range from 3.6% to 6.9% in symmetric delamination and from 2.9% to 5.9% in asymmetric delamination. So we conclude that the longer the delamination, the higher the estimation.


### 6.6. Results: Sensor B

Figure12 represents different wave modes in the case where A0 mode is incident on the transmission side of the main laminate at sensor B.Figure 12
Out-of-plane displacement, taken on damaged composite[0]4 at x = 300 mm, case symmetric delamination length 60.This figure represents different wave modes in the case where A0 mode is incident on the transmission side of the main laminate at sensor B. A new mode S0 was generated when the A0 interacted with the entrance of delamination as shown in Figure 8. This new mode traveled within the upper sublaminates along the length of delamination. Then it reached the exit of delamination, interacted with the exit, and generated a new mode A0 in the forward direction. This new mode propagated in the main laminate, called the A0S0UA0 mode. In the lower sublaminate the mode is called A0S0LA0. Two modes arrived at the same time in the case of symmetric delamination.


### 6.7. Characterization of Delamination from Data Acquired at Sensor B

Figure13 presents wavelet coefficients versus time at sensor B. The Δt1 is the time delay between two transmitted symmetric modes. The ∆t2 is the time delay between the first transmitted symmetric mode and the third antisymmetric mode, so knowing the group velocity of the mode S0: Vg (S0, 300 kHz) =5649 m/s, we can calculate the delamination length noted LD1 by (7) and LD2 by (8). Table 3 resumes the results found.Table 3
Delamination localisation from data acquired at sensor B.
Delamination Position Length L D 1 Error % L D 2 Error % Min Error % 1 and 2 60 49,1 18,1 62,1 3,5 3,5 90 93,5 3,9 95,2 5,7 3,9 120 127,1 5,9 127,7 6,4 5,9 2 and 3 60 45,5 24,2 62,9 4,9 4,9 90 98,9 9,8 94,0 4,5 4,5 120 121,5 1,2 122,5 2,1 1,2Figure 13
Coefficients of wavelet gaus1, case symmetric delamination length 60.The sensor B permits evaluating only the length of symmetric and asymmetric delamination. The errors associated with the predicted symmetric delamination length range from 1.2% to 4.9%. In the case of asymmetric delamination the errors range from 3.5% to 5.9%. So we conclude that the longer the delamination, the higher the estimation. Also, the sensor B gives accurate estimation in case of symmetric delamination.


## 6.1. Model:04 without Defect

The A0 mode was excited from the left edge of the plate by concentrated force at the actuator C (see Figure 4). The A0 Lamb wave is the focus in this study because of its sensitivity to the small defects. To verify the excitability of pure mode A0, 2D-FFT was performed at the sensor A. A matrix was calculated at uniformly spaced points.Figure 4
Numerical model without defect, C: actuator, A and B: sensors.Figure5 shows that there is only one mode propagating which is A0.Figure 5
Superposition of 2D-FFT over analytical dispersion curves.


## 6.2. Model: Symmetric Delamination 60 x 0.1

Figure6 is the delaminated model. Figure 7 is a snapshot of FE out-of plane displacement. The Lamb waves split into two parts and travel independently with different velocities in the upper and the lower sublaminates. Mode conversions happen at the two edges of a delamination due to the change of boundary conditions. To characterize the delamination, we analysed the signals picked up at the sensors A and B.Figure 6
Numerical model of symmetric delamination.Figure 7
Snapshot of FE simulated out-of plane displacement in the upper and the lower sublaminates separated by the delamination.


## 6.3. Modes Conversion

In Figure8, the first wave group is a forward traveling wave group, which is an excited A0 mode atx= 0. When this interacted with the entrance of the delamination, there was a low amplitude reflected wave group. This reflected wave group appears as the second wave group A0A0. Delamination separates the plate into two sublaminates: upper and lower. The transmitted wave converted to the symmetric and antisymmetric mode A0S0 and A0A0 which traveled within the sublaminates with different velocities. The first wave A0S0U—U, L mean that the mode travel, respectively, in the upper and lower sublaminate—reached first the exit of delamination; this is justified by the fact that symmetric mode is faster than the antisymmetric mode. The backward mode A0S0US0U passes through the upper and then is transmitted to the main laminate as the second wave group A0S0US0UA0. By the same manner, in the lower sublaminate the mode A0S0LS0L is transmitted to the main laminate as A0S0LS0LA0. In the case of symmetric delamination the two modes A0S0US0UA0 and A0S0LS0LA0 are reflected at the same time.Figure 8
A0 reflection at the entrance and at the end of asymmetric delamination (length 60), case A0 incident; R denotes reflected mode and T denotes transmitted mode.


## 6.4. Results: Sensor A

Figure9 is A-scan image from the numerical simulations, taken at 100 mm for a delamination length of 60 mm. As explained in Figure 8, we conclude that the first arrived wave is the incident mode A0 and the second is the first reflected wave A0A0. The delamination is symmetric so two waves arrived at the same time as one packet which is A0S0US0UA0+ A0S0LS0LA0.Figure 9
Out-of-plane displacement, taken on damaged composite[0]4 at x = 100 mm, case symmetric delamination length 60.Based on the 2D-FFT computed using Matlab program applied to the surface of delaminated skin (before the delamination), Figure10 shows that the incident and the reflected modes are antisymmetric which confirms that the first reflected mode is A0A0.Figure 10
Superposition of 2D-FFT applied before the delamination over analytical dispersion curves.


## 6.5. Characterization of Delamination from Data Acquired at Sensor A

Figure11 presents wavelet coefficients versus time at sensor A. The time delay between the peaks permits estimating the position and the length of delamination.Figure 11
Coefficients of wavelet gaus1, case A0 incident, delamination length 60, between layer 2 and layer 3.Knowing the group velocity of the modes A0 and S0: Vg (A0,300kHz)=1268 m/s, Vg (S0,300kHz)=5649 m/s, we can calculate the delamination position from the sensor A noted PD by Δt1 (Delta 1) and delamination length noted LD by Δt2 (Delta 2) as presented in Table 2.Table 2
Delamination localisation from data acquired at sensor A.
Delamination Position Length P D P D T H Error % L D Error % 1 and 2 60 103,8 100 3,8 62,0 3,3 90 103,9 100 3,9 93,3 3,7 120 103,8 100 3,8 123,4 2,9 2 and 3 60 103,2 100 3,2 62,2 3,7 90 106,6 100 6,6 94,1 4,5 120 103,1 100 3,1 120,9 0,8The sensor A permits evaluating the position and the length of symmetric and asymmetric delamination. The errors associated with the predicted symmetric delamination position range from 3.1% to 6.6%. However in the case of asymmetric delamination, the error is less than 3.8%. As far as the length prediction is concerned, the errors range from 3.6% to 6.9% in symmetric delamination and from 2.9% to 5.9% in asymmetric delamination. So we conclude that the longer the delamination, the higher the estimation.


## 6.6. Results: Sensor B

Figure12 represents different wave modes in the case where A0 mode is incident on the transmission side of the main laminate at sensor B.Figure 12
Out-of-plane displacement, taken on damaged composite[0]4 at x = 300 mm, case symmetric delamination length 60.This figure represents different wave modes in the case where A0 mode is incident on the transmission side of the main laminate at sensor B. A new mode S0 was generated when the A0 interacted with the entrance of delamination as shown in Figure 8. This new mode traveled within the upper sublaminates along the length of delamination. Then it reached the exit of delamination, interacted with the exit, and generated a new mode A0 in the forward direction. This new mode propagated in the main laminate, called the A0S0UA0 mode. In the lower sublaminate the mode is called A0S0LA0. Two modes arrived at the same time in the case of symmetric delamination.


## 6.7. Characterization of Delamination from Data Acquired at Sensor B

Figure13 presents wavelet coefficients versus time at sensor B. The Δt1 is the time delay between two transmitted symmetric modes. The ∆t2 is the time delay between the first transmitted symmetric mode and the third antisymmetric mode, so knowing the group velocity of the mode S0: Vg (S0, 300 kHz) =5649 m/s, we can calculate the delamination length noted LD1 by (7) and LD2 by (8). Table 3 resumes the results found.Table 3
Delamination localisation from data acquired at sensor B.
Delamination Position Length L D 1 Error % L D 2 Error % Min Error % 1 and 2 60 49,1 18,1 62,1 3,5 3,5 90 93,5 3,9 95,2 5,7 3,9 120 127,1 5,9 127,7 6,4 5,9 2 and 3 60 45,5 24,2 62,9 4,9 4,9 90 98,9 9,8 94,0 4,5 4,5 120 121,5 1,2 122,5 2,1 1,2Figure 13
Coefficients of wavelet gaus1, case symmetric delamination length 60.The sensor B permits evaluating only the length of symmetric and asymmetric delamination. The errors associated with the predicted symmetric delamination length range from 1.2% to 4.9%. In the case of asymmetric delamination the errors range from 3.5% to 5.9%. So we conclude that the longer the delamination, the higher the estimation. Also, the sensor B gives accurate estimation in case of symmetric delamination.


## 7. Experimental Verification

We validate our model based on Barra results [23]; the author studied by FEM the A0 Lamb mode interaction with delamination in sandwich skin. Based on Fourier Diamond [26, 27], he has estimated the delamination length. The obtained results are validated by experimental study. These studies have been done by the propagation of Lamb waves on air by using a Laser Doppler velocimetry.To validate our model and our method we followed the same methodology, as presented in Figure14.Figure 14
Fourier Diamond.The S(x, t) is the measured time-displacement matrix at the skin surface. IFFT is the Inverse Fast Fourier Transform.


### 7.1. Model: Delamination 24 x 0.1

The B-scan illustrated in Figure16 represents the S(x, t): time-displacements matrix picked up at the surface of limited zone (see Figure 15).Figure 15
Model geometry.Figure 16
(a) Wavefields of A0 mode interaction with delamination. (b) Cut of S(x, t) at the position 185 mm.Figure16(a) indicates all the interactions of A0 Lamb mode with delamination in the skin. Figure 16(b) is the temporal evolution of out-of-plane displacement on the upper surface at x= 185 mm (blue line). To evaluate the defect length, and following the Diamond method, we have first of all separated the incident and the reflected waves. Then the separated incident waves get filtered to remove the undesirable reflections.Figure17 shows the filtered incident waves. We can easily see the mode perturbation in the studied zone. This is due to the delamination presence. To evaluate the defect length we computed the temporal FFT to the filtered signal as presented in Figure 18.Figure 17
Filtered wavefields of incident waves.Figure 18
S(x, f) temporal FFT applied to the filtered incident waves.Figure18 shows the principal lobe localized at the frequency 300 kHz which verifies the excitation frequency. To measure the delamination length, a cut of the temporal FFT at the frequency 300 kHz must be done to show clearly the peaks of the lobes.Figure19 permits measuring the delamination length. We measure it between the peaks of the lobes. The delamination length is 22 mm.Figure 19
Cut of S(x, f) at the frequency 300 kHz.To conclude, the model gives an accurate result with an error of 8.3%. The errors are due to the approximated numerical model.


## 7.1. Model: Delamination 24 x 0.1

The B-scan illustrated in Figure16 represents the S(x, t): time-displacements matrix picked up at the surface of limited zone (see Figure 15).Figure 15
Model geometry.Figure 16
(a) Wavefields of A0 mode interaction with delamination. (b) Cut of S(x, t) at the position 185 mm.Figure16(a) indicates all the interactions of A0 Lamb mode with delamination in the skin. Figure 16(b) is the temporal evolution of out-of-plane displacement on the upper surface at x= 185 mm (blue line). To evaluate the defect length, and following the Diamond method, we have first of all separated the incident and the reflected waves. Then the separated incident waves get filtered to remove the undesirable reflections.Figure17 shows the filtered incident waves. We can easily see the mode perturbation in the studied zone. This is due to the delamination presence. To evaluate the defect length we computed the temporal FFT to the filtered signal as presented in Figure 18.Figure 17
Filtered wavefields of incident waves.Figure 18
S(x, f) temporal FFT applied to the filtered incident waves.Figure18 shows the principal lobe localized at the frequency 300 kHz which verifies the excitation frequency. To measure the delamination length, a cut of the temporal FFT at the frequency 300 kHz must be done to show clearly the peaks of the lobes.Figure19 permits measuring the delamination length. We measure it between the peaks of the lobes. The delamination length is 22 mm.Figure 19
Cut of S(x, f) at the frequency 300 kHz.To conclude, the model gives an accurate result with an error of 8.3%. The errors are due to the approximated numerical model.


## 8. Conclusions

In this work, we have studied the numerical model describing the propagation of antisymmetric mode A0 within a sandwich skin. Simulations based on finite element method were carried out using ABAQUS software in order to investigate the interaction of Lamb wave mode by symmetric and asymmetric through-thickness delamination. The results showed the conversion phenomenon. To quantify the conversion and to characterize the delamination, we have computed the 2D-FFT in order to identify the propagating and the converted modes; after that, the WT was calculated at sensors A and B to measure the arrival times.Table4 compared the delamination length calculated at the two sensors A and B. It shows that sensor A gives more accurate estimation than sensor B. Moreover, sensor A permits also estimating the start of the delamination. So the optimal position to evaluate the position and the length of delamination is between the actuator and the defect.Table 4
Comparison between results at sensors A and B.
Delamination Position Length Sensor A Sensor B LD Error % L D Error % 1 and 2 60 62,0 3,3 62,1 3,5 90 93,3 3,7 93,5 3,9 120 123,4 2,9 127,1 5,9 2 and 3 60 62,2 3,7 62,9 4,9 90 94,1 4,5 94,0 4,5 120 120,9 0,8 121,5 1,2 3 and 4 60 62,0 3,3 62,1 3,5 90 93,3 3,7 93,5 3,9 120 123,4 2,9 127,1 5,9Moreover, comparing the signals obtained at two sensors: A and B, it was found that the mode conversions happen at the edges of a delamination for A0 wave. In fact, this mode is sensitive to the position and to the length of delamination; this makes it an effective tool to evaluate the delamination in orthotropic materiel. Besides, based on signal changes with the delamination edges, a localization method is proposed to estimate the position and the length of the delamination. The proposed method was validated based on previous experimental work.These are the main results:(i)
The estimation is accurate when the delamination length increases.(ii)
The optimal position to evaluate the beginning and the length of delamination is between the actuator and the defect which corresponds to the pulse-echo method.As perspective, we model the whole sandwich structure and we characterize the defect by reflection and transmission coefficients.



---

*Source: 9705407-2018-11-28.xml*