Analysis of Anomalous Dynamic Responses of Fiber Metal Laminates Under Pulse Loading

Composed of metallic layers and composite plies, fiber metal laminates (FMLs) combine inherent advantages of metals and composites. The phenomenon of anomalous response behavior, in which the permanent deformation is in the opposite direction to the incoming pressure wave, has been discovered in dynamic responses of monolithic plate subjected to impulsive loading. In this study, dynamic responses of FML plates of various configurations are examined and compared against the responses of a monolithic plate using finite element analysis. It is found that under the load condition in which an anomalous dynamic behavior of an aluminum plate is observed, FMLs oscillate a few cycles before resting on permanent deformation, and the laminate thickness strongly affects the final deflection, which could be in the counter-intuitive direction or near the initial position of zero deflection depending on configuration and specification of FMLs. In addition, interaction of damage accumulation with deflection responses is investigated. The findings of this study can be useful for optimal design of FMLs intended for usage under extreme loadings.


INTRODUCTION
Composite materials are composed of two or more different substances, each with its own characteristics, combined to achieve superior properties than the constituent materials (Lapczyk and Hurtado, 2007;Guo et al., 2021). Having witnessed a rapidly development, composite materials are widely used in erospace, automotive, electrical, construction, and other fields (Ameri, Moradi and Talebitooti, 2020;Wanchoo et al., 2021). Active research on the dynamic behavior of composite materials includes rate dependence of mechanical properties , constitutive models (Zhan et al., 2021), damage evolution such as shear softening upon dynamic loading (Tang et al., 2017), and etc.
Fiber-metal laminates (FMLs) are made of thin metal layers alternating with composite layers having a unidirectional, cross-ply, or fabric structure (Vasiliev and Morozov, 2018). FMLs combine advantages of fiber-reinforced composite materials and metals, possessing high specific strength, fatigue resistance, and excellent impact resistance (Andrew et al., 2019;Nicolinco et al., 2021). A particular important application of FMLs is using them in conditions that involving severe and extreme loads to absorb energies (Sasso et al., 2019;Jia et al., 2021). The dynamic responses of FMLs under various loading conditions have drawn research interest. One category is the study of projectile impacts, including FMLs under low-speed and high-speed projectile impacts (Lee et al., 2018;Li et al., 2018;Sharma et al., 2021), multiple impacts with the same total energy (Yao et al., 2019), the influence of projectile deformability (Sangsefidi et al., 2021), and etc. Impulsive loading is characterized by high intensity and short duration. While metal plates made of aluminum alloys or steel (Chen and Hao, 2014;Micallef et al., 2016;Cerik, 2017;Zhang et al., 2020), stiffened structures (Zhao et al., 2020;Fu et al., 2021;Kong et al., 2021), as well as other innovative forms of structures Zhou and Jing, 2020;Jing et al., 2021) under impulsive loads have been intensively studied, FMLs subjected to intensive short duration loading, however, are far from being well understood. The involved complexities include the wide range of changes in strain rate under blast loading Xu et al., 2020), interaction of fragments with blast waves (Dhari, 2021), modeling of interfacial debonding between adjacent plies (Soutis et al., 2011), and sensitivity to the load characteristics (Karagiozova et al., 2010).
Anomalous dynamic responses, also knowns as counterintuitive behaviors, have been found in ductile metal objects subjected to a pressure pulse. The final deflections are in the direction opposite to that of the pressure applied during the pulse. Since the discovery of this interesting dynamic behavior (Symonds and Yu, 1985), theoretical, experimental, and numerical studies have been conducted on the counterintuitive behavior of beams and plates under impact and impulse loads. It has been realized that the counter-intuitive behavior is related with the elastic-plastic non-linearity and the compressive instability (Dong et al., 2011). Under impulsive loading conditions, the negative phase due to the momentum of air that generates an overexpansion has been found to have significant influence on the occurrence of counter-intuitive behaviors (Aune et al., 2016). Attention to counter-intuitive behavior is deserved because it is an unstable configuration that a small loading less than the failure level might leads to catastrophic results (Ma et al., 2015). However, the responses of FMLs are not clear in the dynamically anomalous regime that has been observed for monolithic metal plates. Hence, unexpected but potentially disastrous consequences of counter-intuitive behavior can be a critical factor restricting extensive usage of FMLs.
In this study, dynamic responses of FML plates subjected to impulsive loads are investigated, with special focus to anomalous response behavior. First, dynamic responses of a monolithic metal plate are studied and validated. Then, same impulsive loading conditions are applied to FML plates of varying configurations. The differences between monolithic and FML plates are examined. Lastly, the damage induced during the dynamic responses is investigated, and implications on FML design and structural applications are discussed. Materials and Methods is devoted to description devotes to the materials and modeling considerations. In Results and Discussions, numerical results are discussed. Conclusions are summarized in Conclusion.

MATERIALS AND METHODS
In this study, FMLs are composed of two thin aluminum layers bonded with fiber composite layers. A coding system denoted by t 1 /t 2 [θ 1 /θ 2 ]/t 3 is adopted. For example, 0.22/0.36[0°/90°]/0.22 indicates that the bottom and top aluminum layers are of a thickness t 1 t 3 0.22 mm and the middle composite layer is of a thickness t 2 0.36 mm with two plies of unidirectional fibers laid orthogonally to one another [0°/90°] with the angle θ 0°p arallel to the x-axis.
In the present study, three types of FML panels with varying layer thickness and fiber orientation are studied as shown in Figure 1. The thickness versus weight trade-off is always an important consider in FML applications. Thus, the total thickness is set to the same value for all FML panels; the performance of different FML configurations can then be evaluated.  Dynamic responses of FML plates subjected to impulsive loads are studied via numerical simulations. The plate dimensions are 300 × 300 × 0.8 mm 3 , and edges are fully constrained in the finite element simulations to model clamped boundary conditions. The density of aluminum is 2.7 g/cm 3 , the Young's modulus is 70 GPa, the Poisson's ratio is 0.33, and the specific heat capacity is 910J/(kg · K), and the material properties of the fiber-reinforced composite are shown in Table 1.

Metal Constitutive Relation
In order to simulate large elastoplastic deformation of metal, the rate of deformation is decomposed into elastic and plastic parts as where D e is the elastic deformation rate, and D P is the plastic deformation rate. The elastic constitutive relationship is defined as where E is Young's modulus, ] is Poisson's ratio, σ ∇J is the Jaumann rate of Cauchy stress, and I is the second order identity tensor.For isotropic materials, the yield function is defined as where σ eq is the equivalent stress, ε p is the equivalent plastic strain, T is the temperature, and where σ ' is the stress deviation, and the plastic strain increment is related to the stress deviator.The yield stress function is defined as (Aune et al., 2017).
where σ 0 is the initial yield stress, Q 1 , Q 2 , C 1 , C 2 are the material hardening parameters, _ ε p 0 is the reference plastic strain rate, and the expression of T p is given as where T r is the ambient temperature, and T m is the melting temperature of the material. Assuming that the plastic potential function coincides with the yield function, the plastic deformation rate tensor can be expressed as Since materials undergo rapid plastic deformation subjected to impulsive loads, the adiabatic process is assumed. The rate of temperature rise is defined as where χ is the Taylor-Quinney coefficient, ρ is the density, and C p is the specific heat. The material parameters used in the metal constitutive relation are shown in Table 2.

Modeling of Fiber Reinforced Polymer
Because of geometrical feature of the FML plates having a thickness significantly smaller than the other dimensions, the aluminum and composites layers are modeled using shell elements. Composite layups can be composed of plies made of different materials in different orientations. To simplify the model, the shell composite layup model in ABAQUS that discretizes only the reference surfaces for each ply is used, and thus fibers and matrix are not explicitly modeled.
The failure criteria of fiber-reinforced composite materials adopt the Hashin's theory (Hashin and Rotem, 1973;Hashin, 1980) that has been implemented in the commercial finite element package ABAQUS. Four different modes of failure are considered including fiber rupture in tension, fiber buckling and kinking in compression, matrix cracking under transverse tension and shearing, and matrix crushing under transverse compression and shearing. The initiation criterion is expressed as where X T denotes the longitudinal tensile strength, X C denotes the longitudinal compressive strength, Y T denotes the transverse tensile strength, Y C denotes the transverse compressive strength, S L denotes the longitudinal shear strength, S T denotes the transverse shear strength, α is a coefficient that determines the contribution of the shear stress to the fiber tensile initiation criterion, and σ 11 , σ 22 , τ 12 are components of the effective stress tensor. A value of d f ≥ 1.0 indicates that the initiation criterion has been met. When damage begins to accumulate, each failure mode has its own damage evolution and a damage index in the range of 0-1 as the measure of damage. Using the energy-based linear damage evolution law, the fracture energies G c ft , G c fc , G c mt and G c mc associated with the energy consumed in the damage process of each mode, are set in the simulations. In addition, material models with stiffness degradation usually lead to difficulties in  (Aune et al., 2017). the convergence of analyses, so the viscosity regularization scheme is used to facilitate convergence. Viscosity coefficients, η ft , η fc , η mt and η mc for each damage mode, are summarized in Table 1.

Modeling Considerations for Adhesive Layers
The degradation of bonding adhesive layer between metal and composite layers in an FML significantly reduce the rigidity and eventually lead to interlaminar delamination. Thus, the adhesive behavior of bonded neighboring layers must be taken into consideration. In this study, the surface-based cohesive behavior model available in ABAQUS is used. The linear elastic traction-separation behavior is assumed in the surface-based cohesive model. Material properties used to define the elastic behavior of adhesive are summarized in Table 3. The debonding process consists of three phases, i.e., damage initiation, damage evolution and complete failure. This study uses the quadratic separation criterion to characterize damage initiation represented as where δ 0 n , δ 0 s , and δ 0 t represent the values of separation, when the separation is either purely along the contact normal or purely in the first or the second shear direction, respectively. δ n , δ s , and δ t represent distances of separation in different directions.
The law of damage evolution describes the degradation rate of stiffness once the initiation criterion is reached. The energy-based exponential damage evolution law is adopted, and the power law criterion is used to describe the interaction of mixed modes as (Camanho and Davila, 2002) where G C n , G C s and G C t refer to the critical fracture energies required to cause failure in the normal, the first, and the second shear directions, respectively. G n , G s and G t represent the work done by the tractions and their conjugate separations. The mixed-mode fracture energy is where where δ f m denotes the displacement at complete failure, δ 0 m denotes the displacement at the beginning of the damage, T eff and δ are the traction and separation, respectively, and G 0 is the elastic energy at damage initiation.

Modeling of Impulsive Loading
A primary goal of this study is to investigate the anomalous dynamic response that is highly sensitive to the applied loads, hence the load function is of critical importance. While many studies use simplified pulse shapes such as triangles or rectangles in their numerical simulations, an in-house code specifically developed for modeling the impulsive pressure is ported to commercial finite element software via user defined subroutines in the present study.
The pressure-time history is comprised of positive and negative phases as shown in Figure 2. The pressure profile is characterized by the ambient pressure p 0 , peak reflected overpressure p rα with the subscript α denoting the angle of incidence, negative overpressure p r,min , arrival time t a , positive duration t + d , negative duration t − d , and decay parameter b. The required input parameters for the impulsive load function are the location of source and equivalent charge weight. The scaled distance using the Hopkinson-Cranz scaling law can be used to describe the intensity of shock loading as where R is the stand-off distance and W is the equivalent weight. The positive phase of the pressure-time history can be described by the modified Friedlander equation expressed as The reflected pressure considering the angle of incidence α can be calculated by p rα p inc 1 + cos α − 2 cos 2 α + p ref cos 2 α where p inc and p ref are incidence and reflected pressures, respectively. By integrating the positive phase of pressure over time, the impulse i + r is obtained Since p inc , p ref , i rα+ are given as functions of the stand-off distance R and weight W, the decay parameter b can be implicitly solved.  (Lapczyk and Hurtado, 2007). The negative phase of the pressure is critical to the occurrence of counter-intuitive behavior; thus, it must be considered in the loading function. A bilinear approximation of the negative phase is adopted (Aune et al., 2017), and the full expression of pressure is defined as Where t − d is the time duration of negative phase, which is typically longer than that of the positive phase, and the coefficient c is introduced to control the proportion of the negative pressure rise time in the entire negative pressure duration.

RESULTS AND DISCUSSIONS Dynamic Response of a Monolithic Metal Plate
In this subsection, the dynamic response of an aluminum plate under impulsive loading is studied and compared against experimental results available in literature. The dimensions of the plate are 300 × 300 × 0.8 mm 3 with edges fully clamped. Two impulsive loads of varying intensity defined in terms of the Hopkinson-Cranz scaled distance, Z 1.46 m/kg 1/3 and Z 1.82 m/kg 1/3 , are considered.
The deflection at the plate center is shown in Figure 3, indicating radically different dynamic responses of the plate subjected to impulsive loads of varying intensity. The final deflection of the plate under the impulsive load of Z 1.46 m/kg 1/3 is around 25 mm. On the other hand, the impulsive load of Z 1.82 m/kg 1/3 yields a negative final deflection around -20 mm with the plate resting on the opposite direction to the incoming impulsive load. This phenomenon is referred to as anomalous or counter-intuitive behavior. Figure 4 shows the final configurations of the plate. Different signs in the out-of-plane (i.e., the z-direction) displacement can easily be noticed. The predicted permanent deflections are validated against the experimental results reported in (Aune et al., 2016). Good agreement between numerical predictions and experiments is reached. It must be emphasized that since the counter-intuitive behavior only occurs in a tight range of impulsive loads, the accuracy of loading function is crucial to the prediction of counter-intuitive behavior. A slight discrepancy might lead to wrong results in predictions of anomalous responses. Consistency with the experimental results in both intuitive and counter-intuitive behaviors underpins the following studies of FMLs subjected to impulsive loading.

Dynamic Responses of Fiber-Metal Laminate Plates
We then study dynamic responses of FML plates and examine the difference in comparison to the responses of the monolithic metal plate. The geometries and material  properties have been summarized in Materials and Methods. The thickness of the FML plates under consideration is identical to that of the aluminum plate studied in Dynamic Response of a Monolithic Metal Plate, so a quantitative comparison in dynamic responses can be made between monolithic and FML plates. Figure 5A plots the deflection time histories at the plate center with the impulsive load of Z 1.82 m/kg 1/3 . Once reaching the maximum positive deflection, all three FML plates rebound toward the opposite direction to the incoming impulsive load. After about two cycles of oscillation from the lowest negative deflection, FML plates progressively arrive at their final configuration. Like the monolithic aluminum plate, the FML plate with a 0.18 mm composite layer reaches a final configuration with a negative deflection at the plate center.
But the deflection is reduced by 42% from 19.16 mm to 11.20 mm. Two FML plates with a 0.22 mm composite layer approach their initial configuration of zero deflection at the plate center. What stands out in Figure 5B is that the thickness of composite layer in an FML only mildly affects the maximum positive deflection but imposes a strong effect on the final configuration. Figure 6 shows the final deformation profiles of FML plates. For the 0.31/0.18[0°]/0.31 FML plate, the deformation profiles along the x-and y-directions demonstrate a significant difference due to the directionality of composite layer. On the other hand, the FML with two orthogonally laid plies of fibers, i.e., the 0.22/0.36[0°/90°]/0.22 plate, demonstrate a much smaller difference in the deformation profiles along the x-and y-directions. More importantly, the slope of the y-direction FIGURE 4 | Final permantent deformation of metal plate subjected to impulsive loads of (A) Z 1.82 m/kg 1/3 that leads to an counter-intuitive behavior and (B) Z 1.46 m/kg 1/3 that leads to intuitive deformation. deformation profile near the edge is considerably higher in the FML plate with only unidirectional lamina. As large deflection slope might eventually lead to tearing of the plate at the boundary, special attention in structural design is thus desired.
Given the considerable difference in the slope of deflection near the plate edges, we further inspect the state of stress triaxiality defined as the ratio of the hydrostatic mean stress to the von Mizes equivalent stress, i.e., σ p σ H /σ eq in which the hydrostatic mean stress defined in terms the principal stresses as σ H (σ 1 + σ 2 + σ 3 )/3. It is generally accepted that shear fracture is the failure mechanism in the negative stress triaxiality regime and fracture due to void formation in the high stress triaxiality regime (Bao and Wierzbicki, 2004). Figure 7A plots the evolution of stress triaxiality at the mid-point of edge perpendicular to the fiber direction. Upon the arrival of pressure pulse, the stress triaxiality abruptly increases to a state of 0.6 in the vicinity of the top cohesive interface, and a negative stress triaxiality of −0.4 is reached near the bottom interface. During the dynamic responses to the impulsive loading, the stress triaxiality intensely oscillates between positive and negative regimes, indicating different failure mechanisms can be activated once a threshold value, such as critical strain, is reached. When the plate eventually rests on its final configuration, the stress triaxialities are positive on the top and negative on the bottom interfaces, respectively. Figure 7B shows the evolution of stress triaxiality at the mid-point of edge parallel to the fiber direction. Positive stress triaxialities develop near both top and bottom cohesive interfaces upon the arrival of pressure pulse. Two phases of oscillation in triaxiality are noticed, first in the range of 0.6 to −0.6 and −0.2 to −0.6 subsequently. When the plate rests on its final configuration, the stress triaxialities are negative on the top and positive on the bottom interfaces, respectively. The clear differences in the state of stress triaxiality observed in Figures 7A,B is closely related with the fiber orientation.
Next, we examine the dynamic responses of FML plates subjected to higher intensity of impulsive loading with the scaled distance Z 1.46 m/kg 1/3 . Figure 8 shows the deflection versus time at the plate center. Recall that the permanent displacement of the aluminum plate is in the intuitive direction under the impulsive load of Z 1.46 m/kg 1/3 (see Figure 3). The 0.31/0.18 [0°]/0.31 plate, however, demonstrate a counter-intuitive behavior with the permanent deflection in the opposite direction to the incident impulsive loading. Hence, it is inferred that the dynamic responses of monolithic and FML plates can be entirely different. Furthermore, unlike the 0.31/0.18 [0°]/0.31 plate, both 0.22/0.36 [0°/0°]/0.22 and 0.22/0.36 [0°/90°]/ 0.22 plates demonstrate an intuitive permanent deformation with positive deflection at the plate center, suggesting significant influences of FML configurations on dynamic responses. The presence or absence of counter-intuitive behavior is affected by many factors, such as material properties (Li et al., 2008), load magnitudes and profiles (Symonds and Yu, 1985;Aune et al., 2017), boundary conditions (Bassi et al., 2003), etc. The ultimate state is extremely sensitive to the initial conditions (Ng and Daolin, 2002). The analyses in the present study highlight additional, but potentially enormous, intricacy induced by   FMLs to the nonlinear responses that deserve attention in structural design and applications.

Damage of Fiber-Metal Laminate Plates
For all three FML plates subjected to the scaled pulse Z 1.82 m/kg 1/3 , composite layers suffer matrix tension damage as shown in Figure 9. Comparing Figures 9A,B reveals the influence of the ratio of the metal to composite layer thickness while the total plate thickness is kept the same for both cases. It is found that the 0.31/0.18[0°]/0.31 plate suffers less damage than the 0.22/0.36[0°/0°]/0.22 plate. Therefore, a thicker aluminum layer attenuates the matrix damage resulted from impulsive loading. Additionally, reducing aluminum layer thickness tends to shift damage from the plate center to the edges along the direction perpendicular to the fiber direction as shown in Figure 9B. Two orthogonally laid plies have a beneficial effect in reducing damage as shown in Figure 9C. In general, the plate center and clamping edges are of high propensity to accumulate damage. To further reduce the extent of damage, the thickness of the FML plates needs to increase.
The degradation and eventual failure of bonding between two cohesive surfaces are major concerns for applications using FMLs. Figure 10 shows the damage evolution at various locations of the 0.22/0.36 [0°/0°]/0.22 plate. Overall, the top cohesive interface, which is closer to the source of impulsive loading, suffers severer damage than the bottom cohesive interface; plate edges sustain more damage than the center, which is closer to the source though. Although damage is produced, no debonding takes place for both Z 1.82 m/kg 1/3 and Z 1.46 m/kg 1/3 .
It is worth pointing out that although the impulsive loading is applied normal to the plate surface, localized sliding can occur between layers. To illustrate this phenomenon, the displacements along the x-and y-centerlines of the plate are plotted in Figure 11. A particular noteworthy finding is that localized sliding can arise at locations that are often overlooked (see the rectangle annotations in Figure 11 The displacements along the x-and y-centerlines of the 0.22/0.36 [0°/0°]/0.22 plate under the impulsive load of Z 1.46 m/kg 1/3 .) especially if the final permanent deflection is small or even close to the zero. The localized sliding due to inconsistency in in-plane displacements might lead to fiber kinking and several other failure types (Greenhalgh, 2009).
Lastly, in order to clearly understand the damage of FML during the dynamic response, the energy dissipated by damage is plotted as shown in Figure 12. In the case of impulsive load Z 1.82 m/kg 1/3 , the energy dissipation due to damage is very close in value for the 0.31/0.

CONCLUSION
Dynamic responses of FML plates subjected to impulsive loads are investigated, with special focus to anomalous response behavior. Modeling considerations include high strain-rate ductile plasticity, failure criteria for fiber reinforced polymer, degradation of adhesive layers, and validated pressure pulses that strongly affect predictions of counter-intuitive behaviors under impulsive loadings. Different FML configurations under impulsive loadings are examined.
First, dynamic responses of a monolithic plate made of aluminum are examined. A counter-intuitive behavior with the final permanent deflection in the opposite direction to the incoming load is predicted numerically. This anomalous behavior vanishes as the loading intensity, which is characterized by the scaled distance, decreases by 20%. The numerical predictions are validated against experimental results available in literature.
Under the same loading conditions, dynamic responses of FML plates of different configurations are studied. Remarkably different responses are noticed. While the maximum positive deflections of FML plates are consistently reduced in all cases, the final permanent deformation shows a wide disparity among FML configurations and load intensities. Under the load condition in which an anomalous dynamic behavior of the aluminum plate is observed, FMLs oscillate before resting on permanent deformation, and the laminate configuration and specification strongly affect the final deflection, which might hold position in the counter-intuitive direction opposite to the incoming pulse or in the vicinity of the initial position with small permanent deflection. Interestingly, after decreasing the scaled distance such that the anomalous behavior vanishes for the aluminum plate, a counter-intuitive behavior still exists in the FML plate that has a thinner composite layer than the other FML plates.
On one hand, the anomalous behavior is related with the elastic-plastic non-linearity and the compressive instability; on the other, interaction with damage accumulation in plates made of FMLs deserves attention. It is found that in the studied loading cases, stress triaxialities in different layers as well as along different directions are noticed to be quite different, suggesting different state of deformation such as tension or shear. Lastly, although matrix damage due to tension is the main failure mechanism in the present study, localized sliding is noticed between layers.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
HZ, JZ, and WL conceived and designed the study. JZ, WL, and BM carried out the simulation work. HZ, JZ, and WL reviewed and edited the manuscript. All authors contributed to the article and approved the submitted version.

FUNDING
This work was financially supported by the Open Collaboration and Innovation Fund of Xi'an Modern Chemistry Research Institute.