Interface Failure of Heated GLARE TM Fiber–Metal Laminates under Bird Strike

: Many high-strength composite materials have been developed for aircraft structures. GLAss ﬁber REinforced aluminum (GLARE) is one of the high-performance composites. The review of articles, however, yielded no study on the impact damage of heated GLARE laminates. This study, therefore, aimed at developing a numerical model that can delineate the continuum damage of GLARE 5A-3 / 2-0.3 laminates at elevated temperatures. In the ﬁrst stage, the inter-laminar interface failure of heated GLARE laminate had been investigated at room temperature and 80 ◦ C. The numerical analysis employed a three-dimensional GLARE 5A-3 / 2-0.3 model that accommodated volumetric cohesive interfaces between mating material layers. Lagrangian smoothed particles populated the projectile. The model considered the degradation of tensile and shear modulus of glass ﬁber reinforced epoxy (GF / EP) at 80 ◦ C, while incorporated temperature-dependent critical strain energy release rate of cohesive interfaces. When coupled with the material particulars, an 82 m / s bird impact at room temperature exhibited delamination ﬁrst in the GF / EP 90 ◦ / 0 ◦ interface farthest from the impacted side. Keeping the impact velocity, interface failure propagated at a slower rate at 80 ◦ C than that at room temperature, which was in agreement with the impact damage determined in the experiments. The outcomes of this study will help optimize a GLARE laminate based on the anti-icing temperature of aircraft.


Impact Damage
Laminated composites are prone to damage when they encounter low-energy impacts, e.g., bird strikes and hailstones [1][2][3]. The frequent damage has undermined to some extent the potential of composite materials, although, their high stiffness, strength-to-weight ratio, good fatigue resistance, corrosion resistance, and formability are appreciable. The impact damage of composite materials appears in dissimilar modes, to name a few, delamination, fiber breakage, and matrix cracks. In most low velocity impact events, the nonvisible damage of composite laminates can be a concern [4], as it can grow depending on the future load spectrum.

GLARE
By contrast, the impact damage of GLAss fiber REinforced aluminium (GLARE) cannot go undetected attributed to the visible dents in the ductile outer aluminum (Al) skin. GLARE is a laminated fiber/metal composite and comprises alternating S2-glass fiber reinforced epoxy-FM94 composite laminates and Al 2024-T3 alloy sheets. The advantages of GLARE manifest in its excellent fatigue tolerance, high impact resistance, and lightening behavior [5]. Each Airbus A380 has about 380 m 2 GLARE in the upper fuselage [6,7]. Strides have been taken to elaborate the thermal functions Aerospace 2020, 7, 28 2 of 25 of GLARE, in particular, to assist aircraft de/anti-icing. At the anti-icing temperature, e.g., 80 • C [8,9], the strength and the stiffness of GF/EP degrade [10], though the ductility of epoxy inclines [11]. Low velocity bird impacts, as a result, are believed to alter the impact resistance of GLARE at an elevated temperature. Yet, the ductile epoxy should delay the inter-ply failure by plasticizing the inter-laminar interfaces. The merits of these conjecture need to be interrogated before drawing a conclusion.

Existing Knowledge
Numerous studies have shown that the extent of composite delamination reduces at an elevated temperature. For example, the spectra/epoxy composites delaminate to a lower extent at 100 • C [12]. Because, the delamination energy augments twofold when the temperature reaches +100 • C from −50 • C [12]. Levin et al. [13] found a lower extent of delamination in carbon fiber reinforced epoxy (CF/EP) composites, when impacted at a test temperature of 70 • C. Hot-wet exposures combined with a frequent temperature change (−55/+80 • C) exacerbate the damage in composites [14][15][16], though, reinforce the delamination resistance by increasing the fracture toughness of resin [11,[17][18][19].

Critical Strain Energy Release Rate
Efforts have been made to characterize the properties of inter-laminar interfaces of various composite laminates at dissimilar temperatures. Quasi-static experiments at 80 • C have demonstrated an increase of critical strain energy release rate (CSERR = G c ) of S2-glass/FM94-epoxy and E-glass/araldite laminate interfaces [10,30]. End notch flexure (ENF) tests on CF/EP (HTA/6376C) composites resulted in an increase of G IIc at 100 • C [31]. In conformity, a test temperature of 90 • C enhanced the G Ic of a symmetric [0 • ] 16/s CF/EP (AS4/3501-6) laminate attributed to the fiber bridging and ductility of epoxy [32], even though, degraded the G Ic of XN1244 adhesive [33]. CF/EP (IM7/8552) laminates, in a four-point quasi-static bending test at 80 • C, showed an increase of G c by 39% compared to that at room temperature [34]. Similarly, G Ic of unreinforced epoxy-3501-6 increased by 27% when the temperature inclined from −50 to +100 • C [35], due to the fact that the fracture evolved to be stable and ductile, and the failure locus shifted to the adhesive interface [36]. The overall outcomes depicted the increase of CSERR of epoxy interfaces at an elevated test temperature.
Noting the temperature as well as the strain rate, the fiber direction of mating composite plies, and the crack length impact the inter-laminar interface properties, a higher strain rate raises the least threshold of G Ic of CF/EP (AS4/3501-6) composites [37], while a 0 • /0 • laminate interface is tougher than a 0 • /90 • interface [38]. As a crack grows, G c exceeds its initial magnitude to reach a plateau that characterizes a steady crack propagation [39][40][41].
deformation of the projectile. The model could approximate the breach of GLARE inter-laminar interfaces by taking the fiber orientation of adjacent plies, temperature-dependent material properties and strain-rate of material in consideration. In companion, the low velocity impact experiments exhibited continuum damage in correspondence to the predicted ones, which demonstrated the reliability of the numerical model.

Geometry
3D CAD software CATIA V5 furnished the sketch of GLARE geometry. The commercial mesh generator HYPERMESH 11.0 was implemented to discretize the geometry. The explicit nonlinear finite element (FE) code LS-DYNA accommodated the FE model of GLARE 5A-3/2-0.3 grade. This cross-ply GLARE configuration was preferred, because it is optimized for structures likely to experience impact loads [43,44]. The stacking sequence [Al/0°/90°/90°/0°/Al/0°/90°/90°/0°/Al] incorporated material layers of in-plane dimension of 2000 × 2000 mm. Respectively, 0.3, 0.125, and 0.01 mm thick Al-layers, GF/EP-plies, and cohesive interfaces comprised the GLARE thickness of 1.9 mm. The cohesive elements were placed between consecutive volume elements of the material layers with nodal coincidence. Figure 1 demonstrates the GLARE laminate and its cohesive interfaces

Discretization
A unique element set discretized each Al-layer and GF/EP-ply. Total 395,032 full integrated Lagrangian solid elements discretized the GLARE laminate. Keeping the interfacial strength, the finer CZEs fail faster [45]. The tendency is opposite when the CZEs are coarser [45]. This non-physical effect of discretization tailbacks to the parameters govern the damage initiation of a cohesive element.
Trial and error convinced the author to implement an element size of 15 × 15 mm, which offered a smooth stress gradient and an intermediate failure rate of cohesive interfaces. It is worth mentioning that this element size is relevant for the geometry, the material, the boundary and the loading conditions specified in this study. Through the thickness, one solid element per GF/EP ply and per interface allowed a layerwise discretization. The solid elements of Al-layers and GF/EP-plies had

Discretization
A unique element set discretized each Al-layer and GF/EP-ply. Total 395,032 full integrated Lagrangian solid elements discretized the GLARE laminate. Keeping the interfacial strength, the finer CZEs fail faster [45]. The tendency is opposite when the CZEs are coarser [45]. This non-physical effect of discretization tailbacks to the parameters govern the damage initiation of a cohesive element. Trial and error convinced the author to implement an element size of 15 × 15 mm, which offered a smooth stress gradient and an intermediate failure rate of cohesive interfaces. It is worth mentioning that this element size is relevant for the geometry, the material, the boundary and the loading conditions specified in this study. Through the thickness, one solid element per GF/EP ply and per interface allowed a layerwise discretization. The solid elements of Al-layers and GF/EP-plies had eight nodes Aerospace 2020, 7, 28 4 of 25 and one integration point, while eight nodes and four integration points charcterized the cohesive elements. Density per unit volume delineated the mass of cohesive elements. 7.86 × 10 5 Lagrangian discrete particles populated the flat faced cylindrical projectile of dimension 63 mm (radius (R b )) × 152 mm (length (L b )) ( Figure 2). The cylindrical projectile emulated a simplified bird. Numerous studies have validated the cylindrical projectile as a bird-analog for a bird strike analysis [46][47][48].
Aerospace 2020, 7, x FOR PEER REVIEW 4 of 25 eight nodes and one integration point, while eight nodes and four integration points charcterized the cohesive elements. Density per unit volume delineated the mass of cohesive elements. 7.86 × 10 5 Lagrangian discrete particles populated the flat faced cylindrical projectile of dimension 63 mm (radius (Rb)) × 152 mm (length (Lb)) ( Figure 2). The cylindrical projectile emulated a simplified bird. Numerous studies have validated the cylindrical projectile as a bird-analog for a bird strike analysis [46][47][48].

GLARE Material
Tables 1 and 2 include the mechanical properties of a unidirectional GF/EP-ply and Al-layer. The composite damage material card_022 accommodated the mechanical properties of S2-glass/FM94epoxy composite plies. The piecewise linear plasticity material card_024 incorporated the strain-rate dependent plasticity of Al. Figures 3a and 3b illustrate the plastic response of Al layer to a tensile stress at different strain rates at room temperature and 80 °C. Tensile stress regime was chosen, since tension promotes crack initiation. The tensile strength of Al was circa 5 MPa lower at 80 °C than that at room temperature for the same strain-rate. Figures 3a,b, therefore, look nearly alike. The preassigned plastic response of Al at different strain rates allowed the numerical code to extrapolate the plastic behavior of Al layers for a beyond threshold strain-rate. It is to notice that the dissipated impact energy in GLARE differs by less than 1.5% when residual curing stresses exist in GLARE [49]. This analysis, in accordance, neglected the minor influence of residual thermal stress on the impact resistance of GLARE. The material properties outlined in Tables 1-3 are taken from an in-house technical report, which has not been made available for public release.

GLARE Material
Tables 1 and 2 include the mechanical properties of a unidirectional GF/EP-ply and Al-layer. The composite damage material card_022 accommodated the mechanical properties of S2-glass/FM94-epoxy composite plies. The piecewise linear plasticity material card_024 incorporated the strain-rate dependent plasticity of Al. Figure 3a,b illustrate the plastic response of Al layer to a tensile stress at different strain rates at room temperature and 80 • C. Tensile stress regime was chosen, since tension promotes crack initiation. The tensile strength of Al was circa 5 MPa lower at 80 • C than that at room temperature for the same strain-rate. Figure 3a,b, therefore, look nearly alike. The pre-assigned plastic response of Al at different strain rates allowed the numerical code to extrapolate the plastic behavior of Al layers for a beyond threshold strain-rate. It is to notice that the dissipated impact energy in GLARE differs by less than 1.5% when residual curing stresses exist in GLARE [49]. This analysis, in accordance, neglected the minor influence of residual thermal stress on the impact resistance of GLARE. The material properties outlined in Tables 1-3 are taken from an in-house technical report, which has not been made available for public release.

Cohesive Zone Material
The cohesive mixed-mode material card_138 incorporated the material properties of FM94epoxy of the CZEs [51]. In line with the literature [52], the stiffness and the strength of FM94-epoxy were decreased with the temperature increase, while the CSERR of FM94-epoxy was increased. Table  3 shows the properties of CZEs at room temperature and 80 °C.

Projectile Material
The null material card_009 accommodated the material properties of the projectile. The projectile volume 1.9 × 10 6 mm 3 and the projectile density 9.38 × 10 −7 kg/mm 3 resulted in the projectile mass (mp) of circa 1.8 kg. The projectile material emulated gelatin, which had a dynamic viscosity coefficient of 4 × 10 −9 . The constitutive law of homogeneous fluidic gelatin reads [53] ij ij ij e P  ργ δ σ where the stress tensor, the pressure, and the Kronecker delta symbol are σij, P, and δij, respectively; the dynamic mass density, the strain, and the strain rate tensor are ρ, γ, and ėij, respectively. The material properties of projectile evolved with the dynamic fluid pressure. The thermodynamic equation of state (EOS) concatenated the fluid pressure (P) with the natural logarithm of volumetric strain (εv) and internal energy (E) of projectile [48,53] (see Equation 2). Here,

Cohesive Zone Material
The cohesive mixed-mode material card_138 incorporated the material properties of FM94-epoxy of the CZEs [51]. In line with the literature [52], the stiffness and the strength of FM94-epoxy were decreased with the temperature increase, while the CSERR of FM94-epoxy was increased. Table 3 shows the properties of CZEs at room temperature and 80 • C.

Projectile Material
The null material card_009 accommodated the material properties of the projectile. The projectile volume 1.9 × 10 6 mm 3 and the projectile density 9.38 × 10 −7 kg/mm 3 resulted in the projectile mass (m p ) of circa 1.8 kg. The projectile material emulated gelatin, which had a dynamic viscosity coefficient of 4 × 10 −9 . The constitutive law of homogeneous fluidic gelatin reads [53] where the stress tensor, the pressure, and the Kronecker delta symbol are σ ij , P, and δ ij , respectively; the dynamic mass density, the strain, and the strain rate tensor are ρ, γ, andė ij , respectively.
Aerospace 2020, 7, 28 6 of 25 The material properties of projectile evolved with the dynamic fluid pressure. The thermodynamic equation of state (EOS) concatenated the fluid pressure (P) with the natural logarithm of volumetric strain (ε v ) and internal energy (E) of projectile [48,53] (see Equation (2)). Here, C is the compression modulus of gelatin. Since the temperature (T) change of gelatin was unlikely in low velocity impact events, two state-variables (C and ε v ) defined the EOS (see Equation (3)).
The compression modulus followed the relative volume of the projectile, as assigned by the tabular EOS in Table 4. C increased at the phase of projectile fragmentation. Table 4. Compressibility and volumetric strain of the gelatin projectile.

Projectile Hydrodynamics
The projectile acted as a fluid domain of moving particles. The SPH kernel evenly dispersed the projectile mass among the constituent particles. Each particle replicated an interpolation point, where the hydrodynamic and the thermodynamic fluid properties were updated in each computation cycle. The analysis scheme interpolated the fluid properties from one particle to the neighbor ones located within a spatial distance, i.e., smoothing length. The smoothing length reads [54] where h is the smoothing length; σ and d are the volumetric stress and the three dimensions of the projectile; m i and ρ i are the mass and the density of a particle denoted as i. The computation cycles updated the density of particles as a function of the compression modulus of gelatin, making the smoothing length time-dependent. A smoothing kernel collected the updated smoothing length in each time step and implemented a cubic spline interpolation function to interpolate the particle properties within the smoothing length. The smoothing kernel (W) was defined as [55] where θ is the incidence angle of the projectile relative to the normal of the GLARE front surface; x i and x j are the spatial positions of neighboring particles within the smoothing length. The constituent hydrodynamic particles of the projectile were described by [55] ( where x i (t) and m i (t) are the time-dependent spatial position and mass of a SPH particle; P in the subscript is the set of moving particles. Temporal evolution of the particle mass in the impact phase and the rebound phase was computed from [54] where, V is the velocity of the ith particle.

GF/EP Composites
The analysis presumed GF/EP-laminates linear elastic in mechanical behavior, according to the findings in [56,57]. For the sake of brevity, the numerical framework did not contemplate the composite damage mechanism, since only inter-laminar interface failures were in the focus of this study. To account for the through-thickness deformation, the constitutive equation of a GF/EP ply was orthotropic, for which the stress-strain relationship in the incremental form was [58,59] where, x, y = 1, 2, 3 are the material directions; ∆σ xy and ∆ε xy are the stress increment and the strain increment, respectively; C xy is the stiffness coefficient. The X-axis (fiber direction) and the Y-axis (transverse-fiber) of a unidirectional GF/EP ply were oriented, respectively, in the 11-direction and the 22-direction. The global Z-axis followed the 33-direction, i.e., the through-thickness direction.

Cohesive Zone
A bilinear traction-separation law governed the failure of cohesive zones both in the peel mode and the shear mode [60], see Figure 4. Damage initiated in the cohesive element given that the total mixed-mode relative displacement (δ 0 m ) was satisfied [61,62] δ 0 m = δ 0 δ 0 I = T/E N and δ 0 II = S/E T denote the displacements to onset the interface failure in mode-I (normal only mode) and mode-II (shear only mode). The ratio of mode mixity (β) is expressed by β was assumed as 0.2, which meant 16.66 % contribution of G IIc to the total G c . This assumption is appropriate for an impact load on laminated composites [63,64], however, is load-case specific.  When the interface failure initiated, two different approaches were implemented to determine the linear softening and the critical (maximum) mixed-mode displacement, at which the cohesive element entirely failed. The power law determined the critical mixed-mode displacement ( max when XMU > 0; where XMU resembles the fracture toughness of the resin interface. If XMU < 0, the Benzeggagh-Kenane law computed max m δ from [62,65] ( ) (13) where EN and ET are the normal and the tangential stiffness of epoxy resin, respectively; T and S are the peak traction in the normal direction and the shear direction, respectively; GIc and GIIc are the CSERR in the corresponding directions.

Translational Constraints
Clamps at the edges prevented the rigid body motion of the GLARE panel. A stand-off distance of 8 mm defined the initial gap between the projectile and the GLARE panel. The projectile was launched with an impact velocity (VI) of 82 m/s at an incidence obliquity of 0° related to the normal of the GLARE surface. The 0° angle of attack (AoA) made sure the most transfer of projectile momentum to the GLARE panel. The general landing speed of passenger aircraft is 82 m/s [66], which was employed as the projectile velocity to reconstruct a bird strike event. When the interface failure initiated, two different approaches were implemented to determine the linear softening and the critical (maximum) mixed-mode displacement, at which the cohesive element entirely failed. The power law determined the critical mixed-mode displacement (δ max m )

Projectile/Target Interaction
when XMU > 0; where XMU resembles the fracture toughness of the resin interface. If XMU < 0, the Benzeggagh-Kenane law computed δ max m from [62,65] δ max where E N and E T are the normal and the tangential stiffness of epoxy resin, respectively; T and S are the peak traction in the normal direction and the shear direction, respectively; G Ic and G IIc are the CSERR in the corresponding directions.

Translational Constraints
Clamps at the edges prevented the rigid body motion of the GLARE panel. A stand-off distance of 8 mm defined the initial gap between the projectile and the GLARE panel. The projectile was launched with an impact velocity (V I ) of 82 m/s at an incidence obliquity of 0 • related to the normal of the GLARE surface. The 0 • angle of attack (AoA) made sure the most transfer of projectile momentum to the GLARE panel. The general landing speed of passenger aircraft is 82 m/s [66], which was employed as the projectile velocity to reconstruct a bird strike event.

Projectile/Target Interaction
An automatic node-to-surface contact defined the interaction between the projectile and the GLARE panel. The contact setting SOFT 1 prevented the penetration of SPH particles in the GLARE surface by updating the contact stiffness based on the global time step size and the nodal mass. The projectile/target interaction was presumed frictionless as contact friction is not discernible at a high strain rate of gelatin projectile [25].

Numerical Implementation
The reference frame of numerical analysis coupled Lagrangian and SPH methods, because they did not involve many derivatives in iteration [67]. For the sake of simplicity, the numerical model did not include the fiber breakage and the matrix cracks of GF/EP plies, thus, eliminated the computational complexity, which arises due to distorted Lagrangian elements. A time step size of 0.01 ms satisfied the convergence criterion of Courant-Friedrichs-Lewy (CFL ≤ 1). Hourglass control type-6 was chosen to prevent the non-physical zero-energy degrees of freedom of adherends after the failure of cohesive elements [68].

Elastoplastic Response of GLARE
Low velocity impact on a GLARE laminate induces multiple delamination through the laminate thickness [44,64,69]. Posterior to impact loading, the delamination area, the axisymmetry and the thickness of generated sub-laminates decide the post-damage buckling behavior of the disintegrated laminate under compression [70]. Thicker sub-laminates possess a higher bending stiffness. Similarly, the impact analysis of GLARE at room temperature and 80 • C predicted debonding/delamination-induced sub-laminates of dissimilar thickness, which will be discussed in the upcoming discourse. An automatic node-to-surface contact defined the interaction between the projectile and the GLARE panel. The contact setting SOFT 1 prevented the penetration of SPH particles in the GLARE surface by updating the contact stiffness based on the global time step size and the nodal mass. The projectile/target interaction was presumed frictionless as contact friction is not discernible at a high strain rate of gelatin projectile [25].

Numerical Implementation
The reference frame of numerical analysis coupled Lagrangian and SPH methods, because they did not involve many derivatives in iteration [67]. For the sake of simplicity, the numerical model did not include the fiber breakage and the matrix cracks of GF/EP plies, thus, eliminated the computational complexity, which arises due to distorted Lagrangian elements. A time step size of 0.01 ms satisfied the convergence criterion of Courant-Friedrichs-Lewy (CFL ≤ 1). Hourglass control type-6 was chosen to prevent the non-physical zero-energy degrees of freedom of adherends after the failure of cohesive elements [68].

Elastoplastic Response of GLARE
Low velocity impact on a GLARE laminate induces multiple delamination through the laminate thickness [44,64,69]. Posterior to impact loading, the delamination area, the axisymmetry and the thickness of generated sub-laminates decide the post-damage buckling behavior of the disintegrated laminate under compression [70]. Thicker sub-laminates possess a higher bending stiffness. Similarly, the impact analysis of GLARE at room temperature and 80 °C predicted debonding/delaminationinduced sub-laminates of dissimilar thickness, which will be discussed in the upcoming discourse.   The stress wave traversed the GLARE panel at a velocity of 4.754 km/s ascribed to the GLARE bulk modulus of 32.97 GPa, which was derived from the elastic modulus at room temperature. The stress wave reverberated at the rear surface (non-impacted side). Since the stress wave got released at the inter-laminar interfaces and the critical stresses of the Al/GF/EP-0 • and the GF/EP 90 • /0 • interface were reached, the entire GLARE laminate disintegrated into sub-laminates within the impact zone ( Figure 5b). Besides, a lower G c at room temperature, compared to that at 80 • C, allowed a faster propagation of interface failure towards the GLARE periphery (Figure 5d-e). As seen in Figure 5f, multiple interface failures enveloped the impact zone, because the mismatch of bending stiffness of sub-laminates reinforced the mode-I (normal mode) interface loading.

Model Verification
To corroborate the outcomes of bird impact analysis, the delineated predictions at room temperature were assessed qualitatively using comparable experimental results published in open literature. Yaghoubi et al. [64] investigated the low velocity impact damage of GLARE 5-3/2-0.3 laminate at room temperature. The 101.6 × 101.6 mm specimens were manufactured from S2-glass fiber reinforced epoxy and Al 2024-T3 alloy. The square specimens were clamped in a fixture, leaving an open circular area of 76.2 mm diameter. Next, a hemispherical steel projectile of 16 mm diameter and 12.91 kg mass was dislodged to exert 30 J impacts on the clamped specimens in a drop-weight impact tower.
In a complementary study, Chen et al. [71] conducted low velocity impact experiments on GLARE 5-3/2-0.25 specimens at room temperature. The GLARE laminate was manufactured from S4-glass fiber reinforced epoxy plies and Al 2024-T3 alloy sheets. The 150 × 100 mm rectangular specimens were clamped in a steel fixture with a cutout of 125 × 75 mm. Next, a solid spherical head of 16 mm diameter and 5.36 kg mass was dislodged in a guided drop-weight tower to impart 70 J energy to the clamped specimens. These two independent experimental campaigns exhibited delamination of the GF/EP-90 • /0 • interface farthest from the impacted side, which agreed with the prediction of the numerical analysis of this study (see Figure 6). Figure 5 exemplifies the response of GLARE laminate to the impact energy of 6.051 kJ. The impact energy resulted from the mp = 1.8 kg and the VI = 82 m/s of the gelatin projectile. The impact response of GLARE corresponded to the properties of material and cohesive interface at room temperature.
The stress wave traversed the GLARE panel at a velocity of 4.754 km/s ascribed to the GLARE bulk modulus of 32.97 GPa, which was derived from the elastic modulus at room temperature. The stress wave reverberated at the rear surface (non-impacted side). Since the stress wave got released at the inter-laminar interfaces and the critical stresses of the Al/GF/EP-0° and the GF/EP 90°/0° interface were reached, the entire GLARE laminate disintegrated into sub-laminates within the impact zone (Figure 5b). Besides, a lower Gc at room temperature, compared to that at 80°C, allowed a faster propagation of interface failure towards the GLARE periphery (Figures 5d-e). As seen in Figure 5f, multiple interface failures enveloped the impact zone, because the mismatch of bending stiffness of sub-laminates reinforced the mode-I (normal mode) interface loading.

Model Verification
To corroborate the outcomes of bird impact analysis, the delineated predictions at room temperature were assessed qualitatively using comparable experimental results published in open literature. Yaghoubi et al. [64] investigated the low velocity impact damage of GLARE 5-3/2-0.3 laminate at room temperature. The 101.6 × 101.6 mm specimens were manufactured from S2-glass fiber reinforced epoxy and Al 2024-T3 alloy. The square specimens were clamped in a fixture, leaving an open circular area of 76.2 mm diameter. Next, a hemispherical steel projectile of 16 mm diameter and 12.91 kg mass was dislodged to exert 30 J impacts on the clamped specimens in a drop-weight impact tower.
In a complementary study, Chen et al. [71] conducted low velocity impact experiments on GLARE 5-3/2-0.25 specimens at room temperature. The GLARE laminate was manufactured from S4glass fiber reinforced epoxy plies and Al 2024-T3 alloy sheets. The 150 × 100 mm rectangular specimens were clamped in a steel fixture with a cutout of 125 × 75 mm. Next, a solid spherical head of 16 mm diameter and 5.36 kg mass was dislodged in a guided drop-weight tower to impart 70 J energy to the clamped specimens. These two independent experimental campaigns exhibited delamination of the GF/EP-90°/0° interface farthest from the impacted side, which agreed with the prediction of the numerical analysis of this study (see Figure 6).  [71], and (c) numerical analysis of this study; 1 and 2 stand for the impacted side and the nonimpacted side, respectively; 3 stands for the GF/EP-90°/0° interface delamination at room temperature.  [71], and (c) numerical analysis of this study; 1 and 2 stand for the impacted side and the non-impacted side, respectively; 3 stands for the GF/EP-90 • /0 • interface delamination at room temperature.
Looking at the differences, the projectile of numerical analysis was fluidic gelatin by contrast to the rigid metal projectiles used in the selected experimental studies. On the other hand, the impact energy in numerical analysis was two orders of magnitude higher than that of the experiments. The rigid projectiles of experiments displayed no deformation upon impact. However, the fragmentation of gelatin dissipated a large amount of impact energy, which otherwise would have been allocated to the impact damage of GLARE. The numerical prediction, therefore, was in fair agreement with the experimental results.
3.1.3. 80 • C (G c80 = G c0 + (39%*G c0 )) Keeping the impact energy at 6.051 kJ, Figure 7 exhibits the impact response of GLARE laminate for the material and interface properties at 80 • C (see Tables 1-3). The analysis considered a 39% increase of G c of cohesive interfaces at 80 • C compared to G c0 at room temperature, as recommended in [34].
Looking at the differences, the projectile of numerical analysis was fluidic gelatin by contrast to the rigid metal projectiles used in the selected experimental studies. On the other hand, the impact energy in numerical analysis was two orders of magnitude higher than that of the experiments. The rigid projectiles of experiments displayed no deformation upon impact. However, the fragmentation of gelatin dissipated a large amount of impact energy, which otherwise would have been allocated to the impact damage of GLARE. The numerical prediction, therefore, was in fair agreement with the experimental results.

80 °C (Gc80 = Gc0 + (39%*Gc0))
Keeping the impact energy at 6.051 kJ, Figure 7 exhibits the impact response of GLARE laminate for the material and interface properties at 80 °C (see Tables 1-3). The analysis considered a 39% increase of Gc of cohesive interfaces at 80 °C compared to Gc0 at room temperature, as recommended in [34]. At 80 °C, the velocity of stress wave in GLARE declined to 4.388 km/s, since the bulk modulus of GLARE degraded to 28.08 GPa at 80 °C. Upon impact, the peel and the shear strength of cohesive interfaces were not exhausted, as a result (Figure 7b). Next, the progressive bending deformation of GLARE reinforced the inter-laminar stresses beyond the failure limit, leading to the onset of GF/EP 90°/0° and 0°/90° interface delamination (Figure 7c). In the successive time steps, the delamination grew at a slower rate in the lateral direction (Figures 7d-f) ascribed to the higher Gc of cohesive interfaces at 80 °C. Because of no mismatch angle between mating plies, the shear stress was the least in the GF/EP 90°/90° interface, which accumulated namely the lowest extent of failure. At 0.14 ms, these two inter-ply delaminations divided the GLARE laminate into three sub-laminates.

80 °C (Gc80 = Gc0 + (45%*Gc0))
Compared to the findings in [34], the author determined through DCB and ENF experiments a 45% increase of Gc of cohesive interfaces at 80°C relative to that at room temperature. To assess the influence of tougher inter-laminar interfaces on the integrity of GLARE laminate, a third analysis was executed by implementing the updated GLARE constituent properties. Applying the identical impact energy, the GLARE deformation and the interface delamination were akin to that exhibited in Figure   Figure 7. Sub-laminates through the GLARE thickness attributed to delamination at (a) 0 ms (initial state), (b) 0.099 ms, (c) 0.11 ms, (d) 0.12 ms, (e) 0.13 ms, (f) 0.14 ms; 1 indicates the impacted side; 2 stands for the GF/EP 90 • /0 • interface delamination; material and interface properties correspond to 80 • C, considering a 39% increase of G c of cohesive interfaces; impact velocity 82 m/s; the fragmented projectile is not included in the imagery to facilitate the pictorial analysis.
At 80 • C, the velocity of stress wave in GLARE declined to 4.388 km/s, since the bulk modulus of GLARE degraded to 28.08 GPa at 80 • C. Upon impact, the peel and the shear strength of cohesive interfaces were not exhausted, as a result (Figure 7b). Next, the progressive bending deformation of GLARE reinforced the inter-laminar stresses beyond the failure limit, leading to the onset of GF/EP 90 • /0 • and 0 • /90 • interface delamination (Figure 7c). In the successive time steps, the delamination grew at a slower rate in the lateral direction (Figure 7d-f) ascribed to the higher G c of cohesive interfaces at 80 • C. Because of no mismatch angle between mating plies, the shear stress was the least in the GF/EP 90 • /90 • interface, which accumulated namely the lowest extent of failure. At 0.14 ms, these two inter-ply delaminations divided the GLARE laminate into three sub-laminates.
Compared to the findings in [34], the author determined through DCB and ENF experiments a 45% increase of G c of cohesive interfaces at 80 • C relative to that at room temperature. To assess the influence of tougher inter-laminar interfaces on the integrity of GLARE laminate, a third analysis was executed by implementing the updated GLARE constituent properties. Applying the identical impact energy, the GLARE deformation and the interface delamination were akin to that exhibited in Figure 7. The onset of delamination surfaced at the same instance (0.12 ms). The delamination growth was very much alike. Furthermore, the sub-laminates were comparable to the sub-laminates of the preceding analysis in terms of thickness and stacking sequence, despite the tougher cross-ply interfaces (as were in this case) appeared to be more associative with the accommodating plies. Due to the slim difference of GLARE deformation for the explicated tougher interface, the pictorial images are not repeated here. Figures 8-10 show the temporal evolution of failure in each GLARE inter-laminar interface. Different colors distinguish the inter-ply failure of distinct interfaces. The interface failure was simulated in the cohesive elements by joining nodes of upper and lower volume elements of mating material layers. The time instance of interface failure was akin to the time step chosen to capture the elastoplastic deformation of GLARE, as aforementioned in Figures 5 and 7. At the end of each time step, all interface failure areas were demarcated to create a 'numerical C-scan'. The interface failure was displayed over an area of 255 × 255 mm, which circumscribed the impact zone.

Delamination of Cohesive Interfaces
The interface failure started, when the numerical scheme eroded the cohesive elements as soon as linear softening exhausted the energy of CZEs. In the incremental time steps, kinetic energy of the projectile transformed in the internal energy of GLARE laminate, making more energy available to reinforce the interface failure in the lateral direction. Noticeably, the interface failure area surfaced in dissimilar shapes, if the material and interface properties were temperature-dependent and the fiber orientation of mating GF/EP plies was dissimilar. The difference in transverse stiffness of Al-layers and GF/EP-plies promoted the mismatch of interface failure enclosures.
More interfaces delaminated with the elastoplastic deformation of GLARE (see , as the imbalance between compression at the impacted side and tension at the non-impacted side augmented the inter-laminar shear stresses. Yet at 80 • C, the analysis predicted significant delamination only of the GF/EP 0 • /90 • interfaces (see the interface failure until 0.14 ms). The GF/EP 90 • /0 • cohesive interface, at quarter normalized GLARE-thickness downstream from the impacted side, incurred maximum failure at room temperature and 80 • C (Figures 8-10). This outcome collapsed on the same line of theoretical approximation, which identified the most delaminated interfaces also at the quarter normalized thickness of a carbon fiber-reinforced composite [70]. Though, cross-ply GF/EP laminates were incorporated in the stacking sequence, the insertion of isotropic Al-layers resulted in a quasi-isotropic behavior of the GLARE laminate. The delamination shape of each inter-ply interface, therefore, was not necessarily oriented in the fiber direction of the lower ply, which is frequently seen in monolithic composite laminates [45].

Manufacturing of GLARE Specimens
The numerical predictions were assessed by conducting low velocity impact experiments on square GLARE 5A-3/2-0.3 specimens of dimension 110 mm (length) × 110 mm (width) × 1.9 mm (thickness). The areal density of the GLARE specimen was 5.6 mg/mm 2 ± 1%. The GLARE panels were manufactured from Al 2024-T3 alloy sheets and S2-glass/FM94-epoxy pre-pregs. The material layers were stacked manually (see Figure 11) and cured in an autoclave at 120 • C temperature and 6 bar pressure for 90 minutes. Square specimens with dimensions of circa 110 × 110 mm were cut using a diamond saw afterwards. Ultrasonic C-scan of the specimens a priori of impact experiments elucidated no manufacturing-induced flaws or cutting-induced interface failure (delaminations or debonding) at the specimen edges.

Manufacturing of GLARE Specimens
The numerical predictions were assessed by conducting low velocity impact experiments on square GLARE 5A-3/2-0.3 specimens of dimension 110 mm (length) × 110 mm (width) × 1.9 mm (thickness). The areal density of the GLARE specimen was 5.6 mg/mm 2 ± 1%. The GLARE panels were manufactured from Al 2024-T3 alloy sheets and S2-glass/FM94-epoxy pre-pregs. The material layers were stacked manually (see Figure 11) and cured in an autoclave at 120 °C temperature and 6 bar pressure for 90 minutes. Square specimens with dimensions of circa 110 × 110 mm were cut using a diamond saw afterwards. Ultrasonic C-scan of the specimens a priori of impact experiments elucidated no manufacturing-induced flaws or cutting-induced interface failure (delaminations or debonding) at the specimen edges.

Impact Test Method
Due to the absence of a bird strike test set-up, the experiments were conducted in a drop-weight impact tower fitted with an automated data acquisition system. This data system registered the impact force from a piezoelectric force sensor placed inside the projectile. An infrared timing gate (see Figure 12d) triggered the data acquisition system and, meanwhile, tracked the travel distance of the projectile to measure the projectile velocity.

Test Set-Up
Penetrating and non-penetrating impacts were performed by clamping the specimens between two rigid plates. The plates had a 100 × 100 mm square aperture. The clamp covered circa 3 mm wide areas of the specimen edges and provided a fixed boundary condition at four ends of the specimen (see Figure 12c). Two clamp reinforcements prevented the dynamic motion of the fixture on the testbed and exerted more contact pressure to decrease slippage at the periphery of specimen (see Figure  12c). A blunt nose conical-shaped steel projectile impacted the GLARE specimen at the center. A rebounce brake prevented multiple impacts on the specimen. The projectile had an approximate dimension of 4 cm (height) × 2.54 cm (diameter) and weighed 0.09 kg. The dead weight attached to the projectile resulted in a total projectile mass of 1.09 kg albeit.
To perform the impact tests at 80 °C, first the specimens were conditioned in an air-circulated environmental chamber (see Figure 12a). The specimens were placed inside the chamber after the temperature and the relative humidity in the chamber reached circa 94.5 °C and 70% RH, respectively. The conditioning time was 45 minutes in line with the early experience that concluded: a longer exposure time is imperative to homogenize the temperature through the thickness of the GLARE specimen. A circa 10 °C higher conditioning temperature was preferred to reach the planned specimen temperature in the test rig. To make sure the specimen temperature was near 80 °C, the

Impact Test Method
Due to the absence of a bird strike test set-up, the experiments were conducted in a drop-weight impact tower fitted with an automated data acquisition system. This data system registered the impact force from a piezoelectric force sensor placed inside the projectile. An infrared timing gate (see Figure 12d) triggered the data acquisition system and, meanwhile, tracked the travel distance of the projectile to measure the projectile velocity.

Test Set-Up
Penetrating and non-penetrating impacts were performed by clamping the specimens between two rigid plates. The plates had a 100 × 100 mm square aperture. The clamp covered circa 3 mm wide areas of the specimen edges and provided a fixed boundary condition at four ends of the specimen (see Figure 12c). Two clamp reinforcements prevented the dynamic motion of the fixture on the test-bed and exerted more contact pressure to decrease slippage at the periphery of specimen (see Figure 12c). A blunt nose conical-shaped steel projectile impacted the GLARE specimen at the center. A re-bounce brake prevented multiple impacts on the specimen. The projectile had an approximate dimension of 4 cm (height) × 2.54 cm (diameter) and weighed 0.09 kg. The dead weight attached to the projectile resulted in a total projectile mass of 1.09 kg albeit.
To perform the impact tests at 80 • C, first the specimens were conditioned in an air-circulated environmental chamber (see Figure 12a). The specimens were placed inside the chamber after the temperature and the relative humidity in the chamber reached circa 94.5 • C and 70% RH, respectively. The conditioning time was 45 minutes in line with the early experience that concluded: a longer exposure time is imperative to homogenize the temperature through the thickness of the GLARE specimen. A circa 10 • C higher conditioning temperature was preferred to reach the planned specimen temperature in the test rig. To make sure the specimen temperature was near 80 • C, the specimens were quickly transferred from the chamber to the test rig. Despite these circumspections, thermal energy of the GLARE specimen partially dissipated in the thermal conductive steel fixture. A visual infrared thermometer showed that the specimen cooled to around 75 • C (see Figure 12b) over the brief time of transferring and clamping.
Aerospace 2020, 7, x FOR PEER REVIEW 17 of 25 specimens were quickly transferred from the chamber to the test rig. Despite these circumspections, thermal energy of the GLARE specimen partially dissipated in the thermal conductive steel fixture. A visual infrared thermometer showed that the specimen cooled to around 75 °C (see Figure 12b) over the brief time of transferring and clamping.  Table 5 demonstrates the test conditions applied to demarcate the severity of damage of GLARE at room temperature and circa 75 °C. Keeping the projectile mass, the drop height (DH) was changed to furnish an impact velocity within 4.7-8.08 m/s. Ascribed to the limited number of specimens, two specimens were tested in each test condition.

Non-Destructive Inspection
A customized water-jet ultrasonic C-scan system checked the impacted specimens. The C-scan system comprised a flat-faced transducer and a receiver with a nozzle diameter of 8 mm and frequency of 10 MHz. The C-scan images elucidated the in-plane damage areas, however, disclosed nothing about the depth of impact damage or the delamination distribution through the GLARE  Table 5 demonstrates the test conditions applied to demarcate the severity of damage of GLARE at room temperature and circa 75 • C. Keeping the projectile mass, the drop height (DH) was changed to furnish an impact velocity within 4.7-8.08 m/s. Ascribed to the limited number of specimens, two specimens were tested in each test condition.

Non-Destructive Inspection
A customized water-jet ultrasonic C-scan system checked the impacted specimens. The C-scan system comprised a flat-faced transducer and a receiver with a nozzle diameter of 8 mm and frequency of 10 MHz. The C-scan images elucidated the in-plane damage areas, however, disclosed nothing about the depth of impact damage or the delamination distribution through the GLARE thickness. This pitfall of the C-scan system is the scattered reflection of sound waves at the interfaces of differentially dense GLARE constituents. Figure 13 exemplifies the scanned damage area of impacted GLARE specimens. As seen in the ultrasonic C-scan images, the blunt-nose projectile inscribed in most of the impact events a pyramid shape damage in the GLARE specimens both at room temperature and 75 • C for impact energies 12 and 20 J. This impact energies corresponded to 2 m and 3 m DH, respectively. For these non-penetration impact events, the projectile kinetic energy dissipated in part in the material damage. The elastic deformation of GLARE conserved rest of the projectile kinetic energy and later, released the stored energy in the bounce-back of projectile. GLARE specimen impacted at room temperature ( Figure 13a) accumulated a 76.89 mm 2 larger damage area (see Table 6) than that inscribed in the GLARE plate at 75 • C (Figure 13b) for an impact energy of 12 J. Similar were the outcomes for the 20 J impacts, in which the damage area of GLARE at room temperature was 55.78 mm 2 larger than that at 77.4 • C (compare Figure 13c to Figure 13d). The difference in the damage area was anticipated due to the plasticizing of epoxy at circa 75 • C that restrained the later damage growth. It is worth noting that, in 20 J impacts, the damage area was relatively smaller compared to that in the 12 J impacts for a comparable test temperature, ascribed to the local concentration of the impact energy at the expense of higher strain rate of material. The augmented material strain rate resulted in shear failure of matrix in the deep dent. By contrast, the 12 J impact energy dissipated via flexural deformation of GLARE and punched a shallow dent in a larger so-called plasticity-induced damage area. In the non-penetration impact events, the damage area enveloped the debonding of rear Al-skin at the non-impacted side.

Experimental Outcomes
When the impact energy sufficed to penetrate the GLARE laminate (as seen for 35.5 J impact energy related to the 4 m DH), the circular cusp of the rear Al-skin cracked parallel to the rolling direction (see Figure 13e), whereas crack in the front Al-skin (at the impacted side) followed the periphery of the indentation zone. Besides, 35.5 J impacts caused fiber breakage and fiber pull out, bending, stretching and thinning of the front and the rear Al-skins. The penetrated rear Al-skin exhibited no debonding. The material damage was visible from outside. The effect of temperature on GLARE damage area nearly disappeared in the penetration events, however (compare Figure 13e to Figure 13f and see Table 6). Of notes, the measured impact force was higher at the elevated test temperature if contrasted to that at room temperature. The relation of impact force to the test temperature persisted for all the inspected impact energies. That meant that the GLARE plate offered a higher impact resistance at the elevated test temperature, which attenuated the damage efficiency of the projectile to a large extent.  Table 5; rows (i), (ii), (iii), and (iv) stand for rear view, front view, side view, and c-scan image, respectively; the legend indicates the attenuation of ultrasonic waves in C-scan.   Table 5; rows (i), (ii), (iii), and (iv) stand for rear view, front view, side view, and c-scan image, respectively; the legend indicates the attenuation of ultrasonic waves in C-scan.

Discussion
High temperature conditioning released the curing-induced residual stress in GF/EP. The attenuation of residual stress suggested softening and re-solidification of fiber-epoxy interfaces and Aerospace 2020, 7, 28 20 of 25 further cross-linking of the epoxy resin, which reinforced the bending stiffness of GLARE to an extent. Due to the increased bending stiffness, the GLARE specimens endured a lower bending deformation and exerted a higher resistant force for a comparable impact energy at 75 • C test temperature. Moreover, the S2-glass fibers were 83 GPa stiffer than the epoxy resin. The out-of-plane deformation disintegrated the fiber/resin interfaces through shear-induced stresses in the indentation zone, as a result. The projectile stretched the detached fibers to break them under tension, thus, less energy was allocated to the interface failure. As a result, the damage area of GLARE surfaced smaller at elevated test temperatures.
High temperature, besides, mobilized the segments of epoxy molecular chain and improved the resin ductility. The enhanced ductility of resin allowed resin plasticization, which implied that the GLARE inter-laminar interfaces behaved tougher. The resin plasticization delayed the interface failure accumulation. Perhaps the progression of interface failure was stable by comparison to an abrupt interface failure of unconditioned GLARE specimens at room temperature [11,17].
Since debonding and delamination of GLARE started earlier at room temperature, the sub-laminates took advantage of membrane stretching. The mismatch in bending stiffness of sub-laminates promoted the mode-I interface loading, inducing a large interface failure. By comparison, the pre-conditioned GLARE specimens, when subjected to a comparable impact energy at 75 • C, were at first intact at the impact center. Because, the mixed modes governed the interface failure and demanded more energy to disintegrate the interfaces.
It is to mention that the intra-laminar transverse matrix cracks can leverage the onset of inter-laminar interface failure. Besides, while cracks are present in the Al layers of GLARE, the glass fibers act as a secondary load path over the crack. This fiber bridging mechanism introduces shear stresses at the Al-GF/EP interfaces, which promote the interface debonding. Excluding the composite damage mechanism did make the model computationally less expensive. The proposed model, however, could have underestimated the extent of debonding/delamination, since the contribution of material damage in the interface failure had not been accounted for.
Another disparity to identify is the gelatin projectile of the numerical analysis compared to the solid rigid steel projectile of the experiment. The numerical impact event, with a gelatin projectile, allowed the dissipation of the projectile kinetic energy to a large extent in the projectile fragmentation. In such a soft-body impact event, the material damage-e.g., fiber failure of GF/EP-would have been negligible. By contrast, the steel projectile of the experiment did not undergo any deformation and transferred the impact energy in large to the material damage. A one-to-one comparison of a soft-body impact to a hard body one would not be ideal, therefore. Yet, the effect of service temperature and fiber orientation on the GLARE interface integrity should be comparable, since these two particulars are independent of projectile type and related strain rate of materials.
Looking at the outcomes, it can be concluded that the numerical predictions conformed to the experimental outcomes, at least in terms of the severity of material damage at different test temperatures. Of notes, this study addressed only the cross-ply GLARE 5A-3/2-0.3 configuration because of its superior impact resistance [63]. A companion study on other GLARE grades of quasi-isotropic or angular fiber orientation will ascertain the extent to which the findings of this study can be extrapolated. Besides, the metal volume fraction (MVF) of GLARE influences the ILSS and the energy dissipation. As found, the permanent dissipated impact energy decreases for increasing the MVF [72], while a single Al sheet bonded with adhesive offers the best ILSS [73]. It would be wise, therefore, to evaluate the combined influence of MVF and operating temperature on the impact resistance of GLARE.

Concluding Remarks
This study highlighted the coupled influence of temperature, fiber orientation, and high strain rate on the interface integrity and plastic deformation of GLARE 5A-3/2-0.3. The strain rate corresponded to low velocity impacts. The investigation on the impact response of heated GLARE outlined some interesting findings, summarized as:

1.
In the non-penetration impact events at room temperature, first the GF/EP 90 • /0 • cohesive interface farthest from the impacted side of the GLARE model disintegrated. Because stress waves reverberated at the rear face, which augmented the cross-ply induced shear stress beyond the limit threshold upon impact. At room temperature, the brittle behavior of cohesive elements accelerated the inter-laminar interface degradation.

2.
At 80 • C, the GLARE model endured interface failure at the periphery of impact zone, though the panel remained coherent at the impact center. In analogy to the interface degradation at room temperature, the GF/EP 0 • /90 • interfaces suffered severe failures at 80 • C. Nonetheless, the delamination growth was delayed at 80 • C. 3.
As explicated in the experiments, FM94-epoxy behaved ductile at 75 • C. As a result, the critical strain energy release rate of FM94-epoxy was higher compared to that at room temperature. It is believed that the higher energy release rate to cohesive failure changed the stick-slip nature of interface degradation to a stable one at 75 • C. The numerical analysis, when implemented a 39% higher G c of cohesive interfaces by contrast to that at room temperature, corroborated the stable failure of the inter-laminar interface.

4.
Thermal conditioning at circa 85 • C for 45 minutes reinforced the bending stiffness of GLARE. As a result, the GLARE plate deflected only locally and offered a higher impact resistance, which restricted the impact damage with the support of resin plasticization.

5.
For the sake of brevity, stress-based intra-ply fiber failure and matrix failure criteria were not included in the composite material model, although intra-ply matrix cracks can be the embryos of inter-laminar interface failure [74]. Scopes are there to complement the proposed model. However, if composite plies fail and are eroded prior to the failure of cohesive elements, convergence difficulties are inevitable. 6.
This study illustrates that a lower scale of interface failure at an elevated temperature may ascertain the residual strength of GLARE at the acceptable level. Therefore, the risk of a total collapse of a thin-walled GLARE airfoil at an anti/de-icing temperature under a soft body impact is not unexpectedly high. 7.
The literature review gives hundreds if not thousands of scientific articles interrogated the impact resistance and the related interface failure of GLARE and other fiber-metal laminates. Most of the studies accommodated experiments at room temperature. Thus, the question arises of whether the delineated energy partition between the different damage modes of GLARE holds also for higher operating temperatures. This study clearly shows the necessity to re-evaluate the impact performance of different GLARE grades based on the operating temperature, than just focusing on the cross-ply GLARE variant. 8.
The outcomes of this study allocate new fields of application for GLARE, e.g., the energy storage of aircraft or shielding against moderate temperature. The thermal energy stored in GLARE could be converted in electric power to execute the secondary operations of an aircraft. In such cases, the landing impact should not impair the interface integrity of GLARE, if the findings of this study are extrapolated in a conservative manner.