Next Article in Journal
Component-Level Modeling of More Electric Auxiliary Power Units for Cooperative Control
Previous Article in Journal
Assessing the Sustainability of Liquid Hydrogen for Future Hypersonic Aerospace Flight
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Study of the Hygrothermal Effects on Low Velocity Impact Induced Indentation and Its Rebound in Composite Laminate

1
State Key Laboratory of Mechanics and Control of Mechanical Structures, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
2
Structural Damage Monitoring Laboratory, Aircraft Strength Research Institute of China, Xi’an 710065, China
*
Author to whom correspondence should be addressed.
Aerospace 2022, 9(12), 802; https://doi.org/10.3390/aerospace9120802
Submission received: 20 October 2022 / Revised: 25 November 2022 / Accepted: 28 November 2022 / Published: 6 December 2022

Abstract

:
Impact indentation is believed to be an effective indication of low-velocity impact (LVI) damage for polymer matrix composites. However, it has been discovered that an indentation can partially rebound over time. Impact indentation and its rebound behavior over a period of time are significantly affected by hygrothermal conditions, especially moisture absorption. Therefore, a good understanding of the moisture-dependent impact indentation and its rebound behavior is helpful for impact damage assessment for composites. In this paper, moisture effects are considered for both the intra-laminar transverse property model and the interlaminar interface model in the simulation of impact indentation. Then, in these two models, viscosities are introduced to represent the indentation rebound over time. In order to validate the proposed models, LVI experiments with different impact energies were conducted on dry and hygrothermal conditioned carbon fiber/epoxy matrix composite laminates. For the specimens, the initial depths of impact dents and their rebounds over time were measured. The specimens of hygrothermal conditions were found with deeper dents compared with dry ones under the same impact energy; and their rebounds were also more significant. These phenomena were explained by the fact that moisture softens epoxy in composite and meanwhile elevates its viscosity. This indentation and its rebound phenomenon were simulated in ABAQUS by considering the moisture effects and viscoelasticity with user-defined material subroutines. These experiments were simulated using the proposed models, and the numerical predictions conformed well with the experimental observations.

1. Introduction

Polymer matrix composites (PMCs), especially carbon-fiber-reinforced plastic (CFRP), have the great advantage of having a low specific weight with excellent mechanical properties. Due to this, the use of CFRP has increased widely in many industries, such as the aerospace and automobile sectors. However, CFRP has one big disadvantage too. It has poor properties in the through-thickness direction, which might degrade upon moisture absorption. Thus, making CFRP prone/vulnerable to low-velocity impacts (LVI) and when subjected to different hygrothermal conditions.
When impacted by a low-velocity object, i.e., any falling tool or impactor like a hail stone, composite laminates behave in a relatively brittle manner, resulting in surface and/or internal damage depending upon the absorption of impact energy. Moisture absorption by the composite laminates further increases the damage manifold. Damage may be in the form of material property degradation, delamination (interlaminar damage), matrix/resin cracking, fiber fracture, and matrix/fiber debonding (intralaminar damage).
A dent is a surface trace of an impact event which could relate to all kinds of damage modes.
In order to predict the interlaminar and intralaminar damage induced by LVI, several methodologies are available in the literature. These include strength-based failure methodologies [1,2,3,4], fracture-mechanics-based methodologies by using the VCCT technique [5,6] or the CZM technique by modeling cohesive (surface-based or element-based) [7,8,9] or spring elements [10] at the interface between two adjacent plies, and CDM-based methodologies [11,12,13,14].
By briefly summarizing some of these methodologies, Zubillaga presented that matrix cracks are responsible for delamination, and by considering fracture toughness and energy release rate of the interface, he proposed a failure criterion [15]. F. Aymerich, F. Dore, and P. Priolo presented a model to predict delamination induced by LVI in cross ply graphite/epoxy laminated plates using a cohesive interface in their paper [7]. C. Bouvet and S. Rivallant in their paper [16] also discussed that matrix cracking is conventionally the first sign of damage to appear upon impact, followed by delamination, which quickly occurs as damage starts to grow.
CFRP laminates exhibit similar contact force and damage magnitudes under transverse quasi-static and LVI loading [17,18,19,20]. Therefore, dynamic impact loading can also be equivalently applied as a quasi-static indentation loading. However, there are some limits where the quasi-static indentation can represent the LVI experiment [21,22].
By discussing indentation further, Karakuzu et al. [23] presented an empirical method to simulate permanent indentation using a numerical study of a glass/epoxy composite plate, which conformed well with the experimental results. In another empirical approach, a study related to indentation and penetration on CFRP laminates was carried out by Caprino et al. [24], which predicted that indentation is an empirical function of impacted energy. He et al. [25] adopted an anisotropic elasto-plasticity theory while considering delamination damage and damping effects for modeling LVI-induced permanent indentation marks on laminated composites.
Surface dents are probably the most direct form of damage observed during a visual check of the impacted site, as the depths of these dents may quantitatively relate to the degree of inside damage. These dents may also rebound depending on the viscoelastic behavior of the composite and the hygrothermal environment during the service life. This might mislead the estimations of impact energy absorbed, as well as the degree of damage during the impact. Therefore, it is necessary to develop a method for the simulation and prediction of LVI damage in terms of dent depth and the rebound of the dent with time in consideration of the environmental conditions.
One of the important properties of PMCs is viscoelasticity, mainly because of the polymer matrix/resin. This means that the dent depth at the impact site is prone to changing over time immediately after the impact. It includes three parts. The first and irreversible part is the impact damage, which causes dislocation of fibers and the matrix. The second part that recovers quickly after the load is removed is the elastic deformation of the whole impacted laminate. The third part, which is responsible for the rebound of the indentation, is the viscoelastic deformation of the resin/matrix of PMCs. It takes much longer for the dent to partially rebound due to viscoelasticity [26]. Furthermore, moisture absorption elevates the viscoelastic material response of composites. These effects are usually not considered in LVI scenarios.
Viscoelastic relations may be expressed in both integral and differential forms, as explained in many published research works. Nima Zobeiry et al. [27] presented a method for FE modeling of isotropic and transversally isotropic viscoelastic materials. Laminae, in a homogenized way, is modeled as transversally isotropic viscoelastic material.
Lin Xiao and Guanhui Wang et al. presented an experimental study related to CFRP subjected to LVI and observed its viscoelastic behavior through the rebound of the indentation dent in terms of depth in post-impact measurements [26]. The relaxation of indentation depth over time is also reported by Wagih et al. [28] in their paper. A significant relaxation of indentation depth during the first 14 days after the test was observed by them.
The material properties of PMCs degrade severely with hygrothermal conditioning, which subsequently affects their performance in LVI and post-LVI scenarios. Moisture is absorbed in composite laminates through Fickian and/or non-Fickian diffusion.
Thomas [29] showed in his paper that, due to viscoelasticity, the dent depth can decrease over time because of humidity and fatigue. Furthermore, just after the impact, this initial dent depth is three times greater than the final rebounded dent depth in some cases. The decrease in dent depth over time also depends on material properties. Wolff EG et al. [30] explained that due to this moisture-induced matrix plasticization in glass/epoxy composites, the viscoelasticity of the composite could be enhanced.
L.S. Sutherland in his review paper [31] explained that the water absorbed by the composite material may result in degradation of the material properties of the laminate. It severely affects the performance of composites in impact and post-impact scenarios. Berketis et al. [32] observed matrix dissolution and interfacial damage in GRP non-crimp glass/polyester specimens after placing these specimens in a hygrothermal condition (water baths at 65 °C) for up to 30 months. Evidence of water absorption causing interfacial or interlaminar material property degradation in CFRP laminates was also observed by Kimpara and Saito [33]. Strait et al. [34] studied the immersion effects for GRP in synthetic sea water at 60 °C and concluded that moisture-induced matrix plasticization can significantly reduce the impact resistance to failure of GRP composites.
The material properties of FRC degrade severely with hygrothermal aging, which subsequently affects the performance of FRC in LVI scenarios. Parvatareddy et al. [35] observed a 70–75% reduction in the residual strength of FRC during their study on the impact resistance and damage tolerance of FRC exposed at 150 °C. Li et al. [36] studied the hygrothermal aging effects on composites and found that under LVI the damaged area increased, which resulted in a decrease in the post-LVI performance of the composites. Mortas et al. [37] studied the effects of corrosive solutions on composites at elevated temperatures. This type of aging severely affects the residual strength of the FRC. Vieille et al. [38] investigated the effects of moisture and temperature on thermoset (TS) and thermoplastic (TP) resin-based composites.
Castaing P et al. [39] depicted that the glass transition temperature, Tg, of epoxy resin/matrix could severely reduce upon moisture absorption. An average of a 20 °C reduction in Tg could occur due to a 1% increment in the absorbed moisture, which is not desirable in a typical glassy epoxy-based polymer. Browning CE et al. [40] found that beyond Tg, the physical properties of the epoxy changed rapidly because of structural changes in the epoxy from a glassy (hard) to a rubbery (soft) state. Chateauminois A et al. [41] showed a linear reduction trend in the glass transition temperature Tg with the increase in moisture content as a result of the direct consequence of plasticization of the matrix/resin in glass/epoxy composites.
Acoustic Emission (AE), a non-destructive testing (NDT) technique used in the aerospace industry, represents the transient elastic sound waves that occur when the material undergoes stress. It is used to monitor damage in composite materials [42] by using piezoelectric sensors applied directly on the surface of samples and capturing these elastic waves. Then, an analysis used to distinguish between the damage types (matrix cracks, fiber breaks, fiber/matrix debonding, and delamination) is carried out on the basis of the collected data. This data can be further analyzed to determine different degradation mechanisms in the material during its lifetime. Life cycle analysis (LCA) of the composite materials can be considered via real time damage tracking by AE. Benzeggah and Barre [43], Chen et al. [44], Kim et al. [45], and Kotsikos et al. [46] used AE signals collected during loading to identify damage mechanisms in composite materials. Their studies proposed that the main damage mechanisms in composite materials can be ranked in ascending order based upon AE amplitude signals. While, the main damage mechanisms are matrix cracking, delamination, debonding, and fiber fracture. Moreover, these AE amplitude signals are dependent on the type of loading and material properties.
In this study, LVI testing of different impact energies (30, 40, and 50 J) is carried out on conditioned and dry specimens to investigate the hygrothermal effects on LVI-induced indentation and its rebound. Then, simulations are performed in ABAQUS with user-defined material subroutines (VUMATs and UMATs). The moisture-dependent intra-laminar viscoelastic model and the moisture-dependent viscoelastic cohesive interface model are implemented through these user-defined material subroutines to demonstrate the initial indentation upon impact and its subsequent rebound over time because of hygrothermal conditioning.

2. Moisture-Dependent Viscoelastic Constitution for Composite Laminate

2.1. Constitutive Equations for a UD Laminae Ply

Initial ideas are taken from the work carried out by Kaliske and Rothert [47] and, further, from Yousaf and Zhou [48]. The generalized Maxwell rheological model (Figure 1) is used for the derivation/formulation of viscoelastic constitutive equations.
For this model, the stress is written as
σ ( t ) = μ 0 ε ( 0 ) + i = 1 N μ i   e x p ( t τ i )   ε ( 0 ) = Γ ( t )   ε ( 0 )
Here, τ i = η μ is defined as the relaxation time for each Maxwell element chain.
σ ( t ) = σ 0 ( t ) + i = 1 N h i ( t )
Here, σ 0 ( t ) is the stress in the elastic/Hookean spring element and h i ( t ) expresses the stresses inside the Maxwell element chains, i.e.,
h i ( t ) = 0 t γ i   e x p ( t s τ i )   σ 0 ( s ) s d s
Here, γ i = μ i μ 0 .
In order to obtain a numerical solution of Equation (3), time-discretized approximation of the 2nd order was implemented in this scheme as
h i n + 1 = e x p ( Δ t τ i ) h i n + γ i ( 1 e x p ( Δ t τ i ) ) Δ t τ i [ σ 0 n + 1 σ 0 n ]
This indicates a very important result that all the values for h i ( t ) are only dependent on the previous values of h i . If these values are known, antecedent values for any given time step can be obtained with the help of the iterative formulation shown in Equation (4). Utilizing the iterative formulation, Equation (2) becomes
σ n + 1 = σ 0 n + 1 + i = 1 N h i n + 1
The three-dimensional formulation of Equation (5) is
σ M n + 1 = σ 0 M n + 1 + M K p = 1 N h M K p n + 1
Here, the elastic stress is defined as
σ 0 M n + 1 = C M K e ε K n + 1
Here, the Cauchy stress tensor is defined as
σ M = [ σ 1 σ 2 σ 3 σ 4 σ 5 σ 6 ]
with σ 11 = σ 1 , σ 22 = σ 2 , σ 33 = σ 3 , σ 12 = σ 21 = σ 4 , σ 13 = σ 31 = σ 5 , σ 23 = σ 32 = σ 6 . The strain tensor with Voigt notation can also be written in the same way as
ε K = [ ε 1 ε 2 ε 3 ε 4 ε 5 ε 6 ]
The stiffness tensor is expressed as C M K e . It defines an elastic relationship between stress and strain. This stiffness tensor is used for the transversally isotropic material as
C M K e = [ E 1 + 4 K T υ 12 2 2 K T υ 12 2 K T υ 12 0 0 0 2 K T υ 12 K T + G T K T G T 0 0 0 2 K T υ 12 K T G T K T + G T 0 0 0 0 0 0 G 12 0 0 0 0 0 0 G 12 0 0 0 0 0 0 G T ]
where K T = 1 / ( 2 ( 1 υ 23 ) E 2 ( 4 v 12 2 E 1 ) and G T = E 2 2 ( 1 + υ 23 ) are the plane strain bulk modulus and transverse shear modulus, respectively [27]. Therefore, the internal stress variables of Equation (4) can be rewritten as
h M K p n + 1 = e x p ( Δ t τ M K p ) h M K p n + γ M K p ( 1 e x p ( Δ t τ M K p ) ) Δ t τ M K p [ C M K e ε K n + 1 C M K e ε K n ]
Now the moisture concentration variables G 0 and G 1 are added in Equation (11), where G 0 indicates the increase or decrease in the instantaneous stiffness and G 1 indicates the increase or decrease in viscoelastic/transient stiffness with the change in moisture concentration. Here,
G 0 = ( C M K e α × c C 0 ) C M K e
and
G 1 = 1 β × c / C 0
where C 0 is the applied moisture concentration (ppm or gm/m3) that corresponds to the relative humidity level, α (GPa) is the slope of the graph of instantaneous modulus C M K e (GPa) versus non-dimensional moisture concentration c / C 0 with a value of c (ppm or gm/m3) which corresponds to a given material point coordinate and time, and β is a constant.
Now Equation (11) can be written as
h M K p n + 1 = e x p ( Δ t τ M K p ) G 1 h M K p n + G 1 γ M K p ( 1 e x p ( Δ t τ M K p ) ) Δ t τ M K p [ G 0 C M K e ε K n + 1 G 0 C M K e ε K n ]
After expanding Equation (6), the Cauchy stress in incremental formulation type is written as
σ M n + 1 = G 0 C M K e ε K n + 1 + M K p = 1 N exp ( Δ t τ M K p ) G 1 h M K p n + G 1 γ M K p τ M K p ( 1 exp ( Δ t τ M K p ) ) Δ t [ G 0 C M K e ε K n + 1 G 0 C M K e ε K n ]
This is an iterative formulation for the stress in the viscoelastic material and it is used in the ABAQUS user material subroutine of UMAT.
In this paper, viscoelastic UD Laminae Ply is simulated by considering only two Maxwell element chains in parallel with the Hookean springs of all the members of the stiffness tensor. Therefore, the constitutive Equation (15) takes the shape of
σ M n + 1 = G 0 C M K e ε K n + 1 + exp ( Δ t τ M K 1 ) G 1 h M K 1 n + exp ( Δ t τ M K 2 ) G 1 h M K 2 n + { 1 + G 1 γ M K 1 τ M K 1 ( 1 exp ( Δ t τ M K 1 ) ) Δ t + G 1 γ M K 2 τ M K 2 ( 1 exp ( Δ t τ M K 2 ) ) Δ t } [ G 0 C M K e Δ ε K ]
where
Δ ε K = ε K n + 1 ε K n
The internal stresses in the Maxwell element chains are updated by using the relationship defined by Equation (14) in each iteration.
Constitutive Equation (16) is further used to calculate the tangent modulus/Jacobian matrix, which is the slope for the stress–strain curve
M K n + 1 = Δ σ Δ ε = { 1 + M K p = 1 2 G 1 γ M K p ( 1 e x p ( Δ t τ M K p ) ) Δ t τ M K p } [ G 0 C M K e       n + 1 ]
where C M K e is the stiffness matrix defined in Equation (10).

2.2. Failure Law of UD Laminae Ply

The in-plane 3D Hashin Failure Criterion is used to simulate lamina-ply level failures, i.e., matrix failure and fiber failure (Table 1), to simulate impact damage during LVI testing. By adding a moisture concentration variable ( G 0 ) I m p a c t in the formulation, the instantaneous stiffness E 0 changes with moisture concentration.
( G 0 ) I m p a c t = ( E 0 α × c C 0 ) E 0
Where σ x x and σ y y are in plane stresses in fiber and transverse directions. τ xy and τ yz are shear stresses. X T and X C are fiber tension and compression strength. Y T and Y C are matrix tension and compression strength. S x y and S y z are shear strength, respectively.

2.3. Constitutive Equations for the Cohesive Interface

A similar approach was extended to formulate constitutive equations for the viscoelastic cohesive interface for initiation, evolution of damage, and rebound of deformation in an impacted laminate. To derive constitutive equations/formulations for the viscoelastic cohesive interface, let us introduce 3D formulation, vector notation adaptation in Equation (5), i.e.,
t M n + 1 = t 0 M n + 1 + M K p = 1 N h M K p n + 1
Here, the elastic stress is defined as
t 0 M n + 1 = E M K e ε K n + 1
Here, the traction stress vector is defined as
t M = [ t 1 t 2 t 3 ]
with the subscripts of 1, 2, and 3 corresponding to normal and two shear directions in the interface. The nominal strain vector can also be written as
ε K = [ ε 1 ε 2 ε 3 ]
The elasticity vector is expressed as E M K e . It defines an elastic relationship between stress and strain. This elasticity matrix is used for the isotropic material as
E M K e = [ K N ( 1 d ) 0 0 0 K S ( 1 d ) 0 0 0 K T ( 1 d ) ]
where K N , K S ,   K T ,   d   are the elastic stiffness values in normal and in two shear directions and the damage variable, respectively. Therefore, the internal stress variables of Equation (4) can be rewritten as
h M K p n + 1 = e x p ( Δ t τ M K p ) h M K p n + γ M K p ( 1 e x p ( Δ t τ M K p ) ) Δ t τ M K p [ E M K e ε K n + 1 E M K e ε K n ]
Now, the moisture concentration variables G 0 and G 1 are added in Equation (25), where G 0 indicates the increase or decrease in the instantaneous stiffness and G 1 indicates the increase or decrease in viscoelastic/transient stiffness with the change in moisture concentration. Here,
G 0 = ( E M K e α × c C 0 ) E M K e
and
G 1 = 1 β × c / C 0
where C 0 is the applied moisture concentration (ppm or gm/m3) that corresponds to the relative humidity level, α (GPa/mm) is the slope of the graph of the instantaneous modulus E M K e (GPa/mm) versus the non-dimensional moisture concentration c / C 0 with a value of c (ppm or gm/m3) that corresponds to a given material point coordinate and time, and β is a constant.
Now, Equation (25) can be written as
h M K p n + 1 = e x p ( Δ t τ M K p ) G 1 h M K p n + G 1 γ M K p ( 1 e x p ( Δ t τ M K p ) ) Δ t τ M K p [ G 0 E M K e ε K n + 1 G 0 E M K e ε K n ]
After expanding Equation (20), the traction stress vector in the incremental formulation type is written as
t M n + 1 = G 0 E M K e ε K n + 1 + M K p = 1 N exp ( Δ t τ M K p ) G 1 h M K p n + G 1 γ M K p τ M K p ( 1 exp ( Δ t τ M K p ) ) Δ t [ G 0 E M K e ε K n + 1 G 0 E M K e ε K n ]
This is an iterative formulation for the stress in a viscoelastic interface and it was also implemented in the UMAT subroutine.
In this paper, the viscoelasticity of interface was also simulated by two Maxwell element chains in parallel with the Hookean springs of all the members of the elasticity vector. Multiple values of bulk modulus and shear modulus were applied in normal and shear directions of the cohesive interface. Therefore, the constitutive Equation (29) takes the shape of
t M n + 1 = G 0 E M K e ε K n + 1 + exp ( Δ t τ M K 1 ) G 1 h M K 1 n + exp ( Δ t τ M K 2 ) G 1 h M K 2 n + { 1 + G 1 γ M K 1 τ M K 1 ( 1 exp ( Δ t τ M K 1 ) ) Δ t + G 1 γ M K 2 τ M K 2 ( 1 exp ( Δ t τ M K 2 ) ) Δ t } [ G 0 E M K e Δ ε K ]
where
Δ ε K = ε K n + 1 ε K n
The internal stresses in the Maxwell element chains are updated by using the relationship defined by Equation (28) in each iteration.
The constitutive Equation (30) is further used to calculate the tangent modulus/Jacobian matrix, which is the slope for the stress–strain curve
E M K n + 1 = Δ t Δ ε = { 1 + M K p = 1 2 G 1 γ M K p ( 1 e x p ( Δ t τ M K p ) ) Δ t τ M K p } [ G 0 E M K e n + 1 ]
where E M K e is the stiffness matrix defined in Equation (24).

2.4. Damage Law of Cohesive Interface

The stress-based quadratic traction–separation failure criterion and linear softening mixed mode B-K law are used to simulate cohesive interface failure initiation and propagation, respectively (Table 2 and Figure 2). The instantaneous elasticity K M 0 changes with the change in moisture concentration.
( G 0 ) I m p a c t = ( K M 0 α × c C 0 ) K M 0
where t n , t s , and t t are interface stresses in normal and two shear directions, respectively. T n f is the interface strength in the normal direction and T i f is the interface strength in two shear directions. δ n   or   δ 3 , δ s ,   δ t are pure mode displacement values in normal and two shear directions, respectively. For the mixed mode, δ m , δ m 0 , δ m f represent the effective displacement value, effective displacement for propagation onset, and effective displacement for total failure. Whereas d and α are variables for damage evolution in terms of delamination and maximum displacement reached historically (a state variable), respectively. K M   is the mixed mode elastic stiffness value.

3. Experiment and FEM Simulation

3.1. Experimental Methodology

Specimens were made of T700/QY9510 composite. The supplier of these composite laminates was “Xie Chuang Composite Material Co., Ltd. Dongguan, China.” The supplier manufactured these in the form of prepregs using a hot-pressing process. The fibers used in these prepregs were carbon fiber HF30 (with properties similar to T700) and the matrix/resin used had properties similar to epoxy EM817. Each specimen had dimensions of 150 mm × 100 mm × 4.97 mm (the thickness of single ply is 0.155 mm approximately) and had a ply-pattern of [45/0/-45/90]4S. A total of 33 such specimens were used.
Pre-conditioning:
All these specimens were heated in an oven at a temperature of 40 °C for 1 h in order to eliminate any moisture before being subjected to different hygrothermal conditions.
After drying, the initial weight of these specimens was measured on a digital mass balance with a precision of up to four decimal points.
Conditioning:
Hygrothermal aging on 24 specimens was completed by selecting the following conditions in an environment conditioning chamber:
I.
25 °C/RH: 85%;
II.
25 °C/RH: 100%.
To determine the moisture absorption by each specimen, the weight gain was measured on an experimental lab ledger at regular time intervals. The percentage weight gain was calculated for each specimen until saturation in accordance with the ASTM D5529/D5229M-92 [52] standard. Initially, weight measurements were undertaken in shorter time intervals of 24–48 h, and later measurement intervals were 1–2 weeks. Six specimens were used to calculate the maximum moisture gain percentage and saturation time for each hygrothermal condition.
The composite material under study here followed Fick’s Diffusion Law in moisture absorption behavior, as indicated by an initial linear increase in the moisture level, thereby steadily trending towards saturation.
Fick s   Diffusion   Law :   D = π ( h 4 M s ) 2 ( M 2 M 1 t 2 t 1 ) 2 = π ( h 4 M s ) 2 k 2
where h is the thickness of the composite, M s is the maximum moisture uptake in percentage at the saturation time, M 1 and M 2 are moisture gain values in percentage at time t 1 and t 2 , respectively. k is the slope of the linear part of the M versus t curve.
A moisture absorption curve can be numerically fitted using Equation (35). Only moisture diffusion along the thickness direction, or z-direction, was considered. Moisture concentrations were applied on both sides, i.e., at thickness t = 0 mm and at thickness t = 5.022 mm, as shown in Figure 3. Other sides/edges of the specimen were coated with sealing paint to stop moisture diffusion. Thus, reducing the case to a one-dimensional moisture diffusion problem. Figure 4 shows the maximum moisture gains of 0.32% and 0.37% at the saturation time for conditioned cases I and II, respectively. From Figure 4, using slope k, the moisture diffusion value of D z = 3.64 × 10−7 mm2/s was calculated along the thickness direction from Equation (34) at a 25 °C temperature and 1 bar pressure.
Fickian   Numerical   Fit :   M = M s [ 1 e x p { 7.3 ( D · t h 2 ) 0.75 } ]
Impact testing and dent rebound measurements:
Then, these dry (at 25 °C) and conditioned specimens were subjected to LVI testing by using an impact machine in accordance with the ASTM D7136M-15 standard [53] to determine the moisture absorption effects in terms of an LVI-induced dent and its rebound phenomenon over time.
The impactor is made of steel and is hemispherical, with a mass of 5 kg and a diameter of 16 mm.
Unidirectional (UD) laminae ply has transversally isotropic material behavior. The material properties of UD lamina ply are presented in Table 3.
The material properties of the cohesive interface are provided in Table 4.
The elasticities and strengths of the specimens presented in Table 3 and Table 4 were offered by the supplier “Xie Chuang Composite Material Co., Ltd., Dongguan, China”.
Each of the three specimens was used in 30, 40, and 50 J impact energy cases for dry and conditioned specimens. Then, averaged values were used to calculate the initial indentation depths and their recovery/rebound with a stipulated period of time in hours. Depth was measured with a depth screw gauge/micrometer.
Dent depth rebound data with time against each impact energy case for dry specimens are presented in Table 5.
The first and the last lines in Table 5 show the initial and the final indentation depths.
Table 6 and Table 7 list dent depth rebound data with time against each impact energy case for specimens at two relative humidity levels.
Similar to Table 5, the first and last lines in these two tables show the initial and final indentation depths for the specimens in two hygrothermal conditions.
The graphical representation of dent depth rebound data with time is shown in Figure 5.
Figure 6 shows a 40J dry impacted specimen from the impacted side and back side. The back side surface of the specimen remained undamaged.
By taking the data points of dent depths for the 50, 40, and 30J impact cases for dry specimens, the following two-exponential decay type equation (Equation (36)) was used in the curve fitting technique, which helped with dent depth calculation in (mm) units for a given time in (hrs) units, as shown in Figure 7:
y = A 1 × e x p ( x t 1 ) + A 2 × e x p ( x t 2 ) + y 0
Similarly, Equation (36) was also acquired in the curve fitting technique for dent depth calculation in (mm) units for a given time in (hrs) units for 50, 40, and 30J impact cases for the conditioned specimens.
The values for the variables used in Equation (36) for each impact energy case for dry and conditioned specimens are presented in Table 8.

3.2. Simulation Methodology

To simulate the indentation depth and its rebound phenomenon in a composite plate, the following methodology was adopted:
The simulation was divided into three parts, as shown in Figure 8.
Firstly, the moisture diffusion of specimen was simulated using mass diffusion analysis in ABAQUS Standard by applying different concentrations relative to different humidity levels. Additionally, a moisture concentration curve along the thickness direction (z-direction) was obtained. The degradation of instantaneous modulus and/or transient/viscoelastic modulus due to moisture concentration was incorporated in constitutive material models and implemented via user-defined material subroutines, i.e., VUMAT and UMAT.
Secondly, the impacts at different energy levels (30, 40, and 50 J) on dry and conditioned specimens were simulated in ABAQUS Explicit as a dynamic modeling case. The user-defined material subroutines (VUMAT for laminae and VUMAT for cohesive interface) were used to simulate the impact-induced damage (in terms of indentation depth) by incorporating a moisture concentration curve along the thickness direction obtained from mass diffusion analysis in the VUMAT formulation.
Thirdly, the same damage for all three cases for dry and conditioned specimens were simulated in ABAQUS Standard as a static modeling case. Moisture-dependent viscoelastic user-defined material subroutines (UMAT for laminae and UMAT for cohesive interface) were used to simulate the same damage and its rebound with time (in terms of indentation depth rebound) by incorporating a moisture concentration curve along the thickness direction obtained from mass diffusion analysis in UMAT formulation.

3.2.1. Simulation for the Moisture Diffusion Case

Each lamina ply and cohesive interface were modeled as 3D diffusivity (heat transfer) elements (designated as DC3D8 in ABAQUS). The FEM model for specimen with a dimension of 150 mm × 100 mm × 5.022 mm was modeled to simulate the moisture diffusion scenario. The moisture diffusivity of D z = 3.64 × 10−7 mm2/s was used for the composite specimen (for both laminae ply and cohesive interface). The steady state condition against each applied moisture concentration presented in Table 9 was achieved at approx. 13 months. The moisture diffusion analysis results after a conditioning time of 7200 h for both types of conditioned specimens are shown in Figure 9 and Figure 10. Once the moisture distribution curve through the thickness of the specimen was achieved, as shown in Figure 11, this curve was then incorporated into the VUMAT/UMAT formulation.
From LVI testing on a dry specimen, a maximum initial dent depth of 0.348 mm was measured against a 50 J impact energy case. So, in order to optimize the conditioning time, the normalized moisture concentration with respect to the applied moisture concentration at certain thickness values (along the z-direction), i.e., t = 0.471 mm (or t = 4.551 mm), and half of the thickness were considered. The normalized moisture concentration values of 0.921, 0.952, 0.991 and 0.726, 0.836, 0.968 were observed at a thickness of t = 0.471 mm (or t = 4.551 mm) and at one half thickness after conditioning times of 3000, 4000, and 7200 h, respectively, for specimens conditioned at 25 °C/RH:85% and 25 °C/RH:100%. The same findings are shown in Figure 12 and Figure 13.

3.2.2. Simulation for the Impact Case

The indentation depth was simulated against each impact energy level for dry and conditioned specimens.
To simulate the impact damage in Abaqus Explicit:
Aerospace 09 00802 g028
There was a much smaller deformation in the impactor relative to that in the composite laminate. Hence, it was negligible. Therefore, the impactor was modeled as a rigid body during the simulation.
The FEM model depicting the boundary conditions and meshing is shown in Figure 14.
The 30, 40, and 50J impact energies were obtained by simulating the impact with the following impactor velocities provided in Table 10.
The values of α (see Equations (19) and (33)) used to incorporate the moisture-dependent variable ( G 0 ) I m p a c t in the VUMAT formulation are presented in Table 11.
Each laminae ply was modeled as continuum solid elements (designated as C3D8 in ABAQUS), while the cohesive interface was modeled using element-based cohesive 3D elements (designated as COH3D8 in ABAQUS) with a very small thickness value of 2 × 10−6. This increases the overall thickness of the model from 4.96 mm (its original thickness) to 5.022 mm. However, this has negligible effects on the simulation results as damage from all different impact energies only involves a maximum of two cohesive interfaces. For LVI testing, the specimen with a dimension of 150 mm × 100 mm was placed at the fixture base over a hollow slot of dimension 125 mm × 75 mm and secured by 4-toggle clamps in the impact testing machine [53]. The fixture base had guiding pins to restrict displacement in x and y directions. In order to simulate the LVI testing, the specimen with a dimension of 125 mm × 75 mm was modeled only in the FEM with pinned boundary conditions (U1 = U2 = U3 = 0) at the edges. By using pinned boundary conditions at the edges, 12.5 mm from each side along the length and width was simply omitted from the FEM. Thus, simulating the specimen that lies exactly over the hollow slot. The interaction between the plate and the impactor was simulated using the general contact (explicit) regime. The smallest element size in the impact zone was 0.20 mm × 0.20 mm. A penalty friction model was included in the contact property definition and a friction coefficient of 0.3 was used. A semi-automatic mass scaling of 1 × 10−6 target time increment was used to reduce the computational time and the element deletion option was imposed.
σ 33 values (stress distribution in the z-direction) for the 50J impact energy for the specimens of three cases, i.e., dry, 85% and 100% relative humidity are shown in Figure 15.
Similarly, σ 33 values (stress distribution in the z-direction) for 40 and 30J impact energies for the specimens of three cases, i.e., dry, 85% and 100% relative humidity are shown in Figure 16 and Figure 17.

3.2.3. Simulation for the Rebound Case

In the third part of the simulation, indentation depth, which was calculated in the dynamic analysis, was achieved in ABAQUS Standard using a static general analysis. For this, the maximum pressure obtained during dynamic analysis was applied (for a certain period) to obtain the same indentation depth, and then allowed the damage (indentation depth) to evolve/rebound with time due to the viscoelastic behavior of the composite for a dry specimen and due to the moisture-dependent viscoelastic behavior of the composite for a conditioned specimen, respectively (Table 12, Table 13 and Table 14).
To simulate indentation rebound, viscoelasticity subject to moisture diffusion in Abaqus Standard was considered:
Aerospace 09 00802 g029
During each impact case for the dry specimen, the specimen absorbs corresponding impact energies. As absorbed energy is different in each case, the rebound is different for each case too. This results in different viscoelastic material behavior of specimens for each impact case. Therefore, it was deemed necessary to calibrate the viscoelastic material properties of specimens for each impact case separately. Table 15 shows the calibrated viscoelastic material properties for each impact case of the dry specimen for both UD ply and cohesive interface materials. In the 1st Maxwell Chain, elements with relatively smaller elastic stiffness values were chosen over a large period for both material types (UD laminae ply and cohesive interface). However, in the 2nd Maxwell Chain, elements with very large elastic stiffness values were chosen over a relatively smaller period for both material types. From Table 15, it can be seen that it is only necessary to slightly calibrate the 2nd Maxwell Chain for both material types as the absorbed energy in each particular impact case was different, which resulted in corresponding initial dent depths and their rebound over a period of time. Hence, the total dent depth Δ d D r y was different for different impact energies.
Viscoelastic material properties, once calibrated for the dry specimen, were used further for predicting the initial dent depth and its rebound over a period of time for the conditioned specimens.
Values of α and β (see Equations (12), (13), (26), and (27)) used to incorporate moisture-dependent variables G 0 and G 1 , respectively, in UMAT formulation are presented in Table 16.
The results for both types of simulations for different impact energy cases for dry and conditioned specimens are shown in Figure 18, Figure 19, Figure 20, Figure 21, Figure 22, Figure 23, Figure 24, Figure 25 and Figure 26.
Figure 18, Figure 19, Figure 20, Figure 21, Figure 22, Figure 23, Figure 24, Figure 25 and Figure 26 explain the simulation results obtained from the impact and rebound simulations in terms of initial indentation depth upon impact and its final indentation depth after the rebound over time, respectively. Additionally, the graphs show the dent evolution and rebound with time.

4. Results and Discussion

The simulation results agree well with the experimental observations in terms of damage (indentation), depth (initial and final), shape, and size for dry and conditioned specimens. However, the simulation results deviate from the experimental results for the case of the dent rebound path (with a maximum error of 37.90%). This decrease in the original depth of the indentation is faster with time initially after the impact, but it slows down with time and eventually stops. From Figure 27 the predicted dent rebound paths, in spite of being curved, look nearly linear due to the longer rebound time with a very small total dent rebound (in the range of 0.034–0.152 mm). However, the experimental curves are characterized as having a rapid rebound in the early period and then these curves become invariable. This phenomenon is more obvious for the 50 J impact energy cases in dry and conditioned specimens.
The rebound in dent depth is mainly due to the viscoelastic behavior of resin/matrix material in UD laminae ply as well as the viscoelastic behavior of the cohesive interface material in dry specimens. The moisture absorption by the conditioned specimens results in a significant increase in the initial dent depth and a significant decrease in its rebound over a period of time as compared to the same in dry specimens. This is mainly because of the moisture-dependent viscoelastic behavior of resin/matrix materials in UD laminae plies and the moisture-dependent viscoelastic behavior of cohesive interface material.
The graphical representation/comparison of the indentation rebound obtained through experimental data and simulations for all cases is shown in Figure 27.
The simulation results accurately predict the initial and final dent depths, when compared with experimental results, that are provided in Table 17 for the dry specimen.
Similarly, the simulation results accurately predict initial and final dent depths, when compared with experimental results, that are presented in Table 18 and Table 19 for both types of conditioned specimens.

5. Conclusions

In this study, experimental data was acquired to understand the damage incurred and its recovery (the rebound of indentation depth with time) of CFRPs subjected to the low-velocity impact of different impact energies under dry and different hygrothermal conditions with values greater than the damage threshold value, and then the same finite element models were simulated to obtain the results, and the following conclusions were obtained:
Specimens of hygrothermal conditions were found with deeper dents compared with dry ones under the same impact energy, and their rebounds were also more significant. These phenomena were explained as moisture softens epoxy in a composite and elevates its viscosity.
In general, the initial indentation depth was proportional to the impact energy and hygrothermal conditioning. Additionally, the indentation rebound was proportional to this initial indentation depth.
Our simulation results accurately predicted the initial and final dent depths for each impact case. While the decrease in the original depth of indentation remained invariable with time, the rebound path curve obtained via simulation decayed at a slower rate compared to the experimental rebound path curve, adding error to the prediction of the rebound path (with a max. error of 37.90%). This error was proportional to the total dent rebound, which was linked with an increase in the impact energy and relative humidity levels.
For the impact cases of 50, 40, and 30J on specimens at [25 °C/RH:85%] and [25 °C/RH:100%], increases of 0.86%, 1.68%, 0.92% and 1.44%, 2.52%, 1.38% in initial indentation depths were observed, respectively, as compared to the same in the dry specimen at 25 °C. Similarly, for the rebound cases of 50, 40, and 30J on specimens at [25 °C/RH:85%] and [25 °C/RH:100%], decreases of 16.54%, 17.65%, 17.68% and 22.69%, 23.53%, 22.01% in terms of final indentation depths after the rebound over time were observed, respectively, in comparison with the same in the dry specimen at 25 °C.

Author Contributions

Conceptualization, M.Y.; methodology, M.Y.; software, M.Y.; validation, M.Y.; formal analysis, M.Y.; investigation, M.Y., Y.Y. and L.W.; resources, C.Z.; data curation, M.Y.; writing—original draft preparation, M.Y.; writing—review and editing, C.Z. and Y.Y.; visualization, M.Y. and L.W.; supervision, C.Z.; project administration, C.Z.; funding acquisition, C.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (11872205), the State Key Laboratory of Mechanics and Control of Mechanical Structures (Nanjing University of Aeronautics and Astronautics, MCMS-E-0221Y02), and the Priority Academic Program Development of Jiangsu Higher Education Institutions.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

LVILow Velocity Impact
PMCPolymer Matrix Composite
B-KBenzeggagh–Kenane
UMATUser Material Subroutine for ABAQUS/Standard
VUMATUser Material Subroutine for ABAQUS/Explicit
CFRPCarbon Fiber Reinforced Polymer
CZMCohesive Zone Modeling
VCCTVirtual Crack Closure Technique
CDMContinuum Damage Mechanics
UDUnidirectional
GRPGlass Reinforced Plastic
FRCFiber Reinforced Concrete
TSThermoset
TPThermoplastic
AEAcoustic Emission
NDTNon-Destructive Testing
LCALife Cycle Analysis

References

  1. Choi, H.Y.; Chang, F.K. A model for predicting damage in graphite/epoxy laminated composites resulting from low-velocity point impact. J. Compos. Mater. 1992, 26, 2134–2169. [Google Scholar] [CrossRef]
  2. Brewer, J.C.; Lagace, P.A. Quadratic stress criterion for initiation of delamination. J. Compos. Mater. 1988, 22, 1141–1155. [Google Scholar] [CrossRef]
  3. Tita, V.; Carvalho, J.; Vandepitte, D. Failure analysis of low velocity impact on thin composite laminates: Experimental and numerical approaches. Compos. Struct. 2008, 83, 413–428. [Google Scholar] [CrossRef]
  4. Farooq, U.; Myler, P. Efficient computational modelling of carbon fibre reinforced laminated composite panels subjected to low velocity drop-weight impact. Mater. Des. 2014, 54, 43–56. [Google Scholar] [CrossRef]
  5. Fleming, D.C. Delamination Modeling of Composites for Improved Crash Analysis; NASA/CR-1999-209725; National Aeronautics and Space Administration, Langley Research Center: Hampton, VA, USA, 1999. [Google Scholar]
  6. Ronald, K. Virtual crack closure technique: History, approach, and applications. Appl. Mech. Rev. 2004, 57, 109–143. [Google Scholar]
  7. Aymerich, F.; Dore, F.; Priolo, P. Prediction of impact-induced delamination in cross-ply composite laminates using cohesive interface elements. Compos. Sci. Technol. 2008, 68, 2383–2390. [Google Scholar] [CrossRef] [Green Version]
  8. Shi, Y.; Soutis, C. A finite element analysis of impact damage in composite laminates. Aeronaut. J. 2012, 116, 1131–1147. [Google Scholar] [CrossRef]
  9. Caputo, F.; De Luca, A.; Lamanna, G.; Borrelli, R.; Mercurio, U. Numerical study for the structural analysis of composite laminates subjected to low velocity impact. Compos. Part B 2014, 67, 296–302. [Google Scholar] [CrossRef]
  10. Bouvet, C.; Castanié, B.; Bizeul, M.; Barrau, J.J. Low velocity impact modelling in laminate composite panels with discrete interface elements. Int. J. Solids Struct. 2009, 46, 2809–2821. [Google Scholar] [CrossRef] [Green Version]
  11. Donadon, M.V.; Iannucci, L.; Falzon, B.G.; Hodgkinson, J.M.; de Almeida, S.F. A progressive failure model for composite laminates subjected to low velocity impact damage. Comput. Struct. 2008, 86, 1232–1252. [Google Scholar] [CrossRef]
  12. Kim, E.H.; Rim, M.S.; Lee, I.; Hwang, T.K. Composite damage model based on continuum damage mechanics and low velocity impact analysis of composite plates. Compos. Struct. 2013, 95, 123–134. [Google Scholar] [CrossRef]
  13. Riccio, A.; De Luca, A.; Di Felice, G.; Caputo, F. Modelling the simulation of impact induced damage onset and evolution in composites. Compos. Part B 2014, 66, 340–347. [Google Scholar] [CrossRef]
  14. Batra, R.C.; Gopinath, G.; Zheng, J.Q. Damage and failure in low energy impact of fiber-reinforced polymeric composite laminates. Compos. Struct. 2012, 94, 540–547. [Google Scholar] [CrossRef]
  15. Zubillaga, L.; Turon, A.; Maimi, P.; Costa, J.; Mahdi, S.; Linde, P. An energy based failure criterion for matrix crack induced delamination in laminated composite structures. J. Compos. Struct. 2014, 112, 339–344. [Google Scholar] [CrossRef]
  16. Bouvet, C.; Rivallant, S. Damage tolerance of composite structures under low-velocity impact. In Dynamic Deformation, Damage and Fracture in Composite Materials and Structures; Woodhead Publishing: Cambridge, UK, 2016. [Google Scholar] [CrossRef]
  17. Luo, G.M.; Lee, Y.J. Quasi-static simulation of constrained layered damped laminated curvature shells subjected to low-velocity impact. Compos. Part B 2011, 42, 1233–1243. [Google Scholar] [CrossRef]
  18. Sutherland, L.S.; Guedes Soares, C. The use of quasi-static testing to obtain the low-velocity impact damage resistance of marine GRP laminates. Compos. Part B 2012, 43, 1459–1467. [Google Scholar] [CrossRef]
  19. Yokozeki, T.; Kuroda, A.; Yoshimura, A.; Ogasawara, T.; Aoki, T. Damage characterization in thin-ply composite laminates under out-of-plane transverse loadings. Compos. Struct. 2010, 93, 49–57. [Google Scholar] [CrossRef]
  20. Brindle, A.R.; Zhang, X. Predicting the compression-after-impact performance of carbon fibre composites based on impact response. In Proceedings of the 17th International Conference on Composite Materials (ICCM17), Edinburgh, UK, 27–31 July 2009. [Google Scholar]
  21. Yigit, A.S.; Christoforou, A.P. Limits of asymptotic solutions in low-velocity impact of composite plates. Compos. Struct. 2007, 81, 568–574. [Google Scholar] [CrossRef]
  22. Wagih, A.; Maimí, P.; Blanco, N.; González, E.V. Scaling effects of composite laminates under out-of-plane loading. Compos. Part A Appl. Sci. Manuf. 2019, 16, 1–12. [Google Scholar] [CrossRef]
  23. Karakuzu, R.; Erbil, E.; Aktas, M. Impact characterization of glass/epoxy composite plates: An experimental and numerical study. Compos. Part B 2010, 41, 388–395. [Google Scholar] [CrossRef]
  24. Caprino, G.; Langella, A.; Lopresto, A. Indentation and penetration of carbon fibre reinforced plastic laminates. Compos. Part B 2003, 34, 319–325. [Google Scholar] [CrossRef]
  25. He, W.; Guan, Z.; Li, X.; Liu, D. Prediction of permanent indentation due to impact on laminated composites based on an elasto-plastic model incorporating fiber failure. Compos. Struct. 2013, 96, 232–242. [Google Scholar] [CrossRef]
  26. Xiao, L.; Wang, G.; Qiu, S.; Han, Z.; Lia, X.; Zhang, D. Exploration of energy absorption and viscoelastic behavior of CFRPs subjected to low velocity impact. Compos. Part B Eng. 2019, 165, 247–254. [Google Scholar] [CrossRef]
  27. Zobeiry, N.; Malek, S.; Vaziri, R.; Poursartip, A. A differential approach to finite element modelling of isotropic and transversely isotropic viscoelastic materials. Mech. Mater. 2016, 97, 76–91. [Google Scholar] [CrossRef]
  28. Wagih, A.; Maimí, P.; Blanco, N.; Costa, J. A quasi-static indentation test to elucidate the sequence of damage events in low velocity impacts on composite laminates. Compos. Part A Appl. Sci. Manuf. 2016, 82, 180–189. [Google Scholar] [CrossRef]
  29. Thomas, M. Study of the evolution of the dent depth due to impact on carbon/epoxy laminates, consequences on impact damage visibility and on in service inspection requirements for civil aircraft composite structures. In Proceedings of the MIL-HDBK 17 Meeting, Monterey, CA, USA, 29–31 March 1994. [Google Scholar]
  30. Wolff, E.G. Moisture effects on polymer matrix composites. SAMPE J. 1993, 29, 11–19. [Google Scholar]
  31. Sutherland, L.S. A review of impact testing on marine composite materials: Part III–Damage tolerance and durability. Compos. Struct. 2018, 188, 512–518. [Google Scholar] [CrossRef]
  32. Berketis, K.; Tzetzis, D.; Hogg, P.J. The influence of long term water immersion ageing on impact damage behaviour and residual compression strength of glass fibre reinforced polymer (GFRP). Mater. Des. 2008, 29, 1300–1310. [Google Scholar] [CrossRef]
  33. Kimpara, I.; Saito, H. Post-impact fatigue behavior of woven and knitted fabric CFRP laminates for marine use. In Major Accomplishments in Composite Materials and Sandwich Structures; Springer: Dordrecht, The Netherlands, 2009; pp. 113–132. [Google Scholar]
  34. Strait, L.H.; Karasek, M.L.; Amateau, M.F. Effects of seawater immersion on the impact resistance of glass fiber reinforced epoxy composites. J. Compos. Mater. 1992, 26, 2118–2133. [Google Scholar] [CrossRef]
  35. Parvatareddy, H.; Tsang, P.W.; Dillard, D. Impact damage resistance and tolerance of high-performance polymeric composites subjected to environmental aging. Compos. Sci. Technol. 1996, 56, 1129–1140. [Google Scholar] [CrossRef]
  36. Li, G.; Pang, S.S.; Helms, J.E.; Ibekwe, S.I. Low velocity impact response of GFRP laminates subjected to cycling moistures. Polym. Compos. 2000, 21, 686–695. [Google Scholar] [CrossRef]
  37. Mortas, N.; Er, O.; Reis, P.; Ferreira, J. Effect of corrosive solutions on composites laminates subjected to low velocity impact loading. Compos. Struct. 2014, 108, 205–211. [Google Scholar] [CrossRef]
  38. Vieille, B.; Aucher, J.; Taleb, L. Comparative study on the behavior of woven-ply reinforced thermoplastic or thermosetting laminates under severe environmental conditions. Mater. Des. 2012, 35, 707–719. [Google Scholar] [CrossRef]
  39. Castaing, P.; Mallard, H. Effect of plasticization with water on the behaviour of longterm exposed filament wounded composites. In Progress in Durability Analysis of Composite Systems; Cardon, A.H., Fukuda, H., Reifsnider, K., Eds.; Taylor and Francis: Abingdon, UK, 1996; pp. 233–239. [Google Scholar]
  40. Browning, C.E.; Husman, G.; Whitney, J.M. Moisture Effects in Epoxy Matrix Composites; ASTM special technical publications: West Conshohocken, PA, USA, 1977; pp. 481–496. [Google Scholar]
  41. Chateauminois, A.; Chabert, B.; Soulier, J.P.; Vincent, L. Hygrothermal ageing effects on viscoelastic and fatigue behaviour of glass/epoxy composites. In Proceedings of the Eighth International Conference on Composite Materials (ICCM8), Honolulu, HI, USA, 15–19 July 1991. [Google Scholar]
  42. Godin, N.; Reynaud, P.; Fantozzi, G. Acoustic Emission and Durability of Composites Materials; ISTE-Wiley: London, UK, 2018; ISBN 9781786300195. [Google Scholar]
  43. Barre, S.; Benzeggagh, M. On the use of acoustic emission to investigate damage mechanisms in glass-fiber reinforced polypropylene. Compos. Sci. Technol. 1994, 52, 369–376. [Google Scholar] [CrossRef]
  44. Chen, O.; Karandikar, P.; Takeda, N.; Kishi, T. Acoustic emission characterization of a glass-matrix composite. Nondestruct. Test. Eval. 1992, 8–9, 869–878. [Google Scholar] [CrossRef]
  45. Kim, S.-T.; Lee, Y.-T. Characteristics of damage and fracture process of carbon fiber reinforced plastic under loading-unloading test by using AE method. Mater. Sci. Eng. 1997, A234–A236, 322–326. [Google Scholar] [CrossRef]
  46. Kotsikos, G.; Evans, J.T.; Gibson, A.G.; Hale, J. Use of acoustic emission to characterize corrosion fatigue damage accumulation in glass fiber reinforced polyester laminates. Polym. Compos. 1999, 20, 689–696. [Google Scholar] [CrossRef]
  47. Kaliske, M.; Rothert, H. Formulation and implementation of three-dimensional viscoelasticity at small and finite strains. Comput. Mech. 1997, 19, 228–239. [Google Scholar] [CrossRef]
  48. Yousaf, M.; Zhou, C. Numerical Study on the Rebound of Low-Velocity Impact-Induced Indentation in Composite Laminate. Aerospace 2022, 9, 651. [Google Scholar] [CrossRef]
  49. Zhang, J.; Zhang, X. An efficient approach for predicting low-velocity impact force and damage in composite laminates. Compos. Struct. 2015, 130, 85–94. [Google Scholar] [CrossRef]
  50. Benzeggagh, M.L.; Kenane, M. Measurement of mixed-mode delamination fracture toughness of unidirectional glass/epoxy composites with mixed-mode bending apparatus. J. Compos. Sci. Technol. 2003, 56, 439–449. [Google Scholar] [CrossRef]
  51. Camanho, P.P.; Davila, C.G.; De Moura, M.F. Numerical simulation of mixed-mode progressive delamination in composite materials. J. Compos. Mater. 2003, 37, 1415–1438. [Google Scholar] [CrossRef]
  52. D 5229/D 5229M-92; Standard Test Method for Moisture Absorption Properties and Equilibrium Conditioning of Polymer Matrix Composite Materials. ASTM International: West Conshohocken, PA, USA, 2004.
  53. ASTM D7136; Standard Test METHOD for Measuring the Damage Resistance of a Fiber Reinforced Polymer Matrix Composite to a Drop-Weight Impact Event. Annual Book of ASTM Standards. ASTM International: West Conshohocken, PA, USA, 2015; pp. 1–16.
Figure 1. Generalized Maxwell rheological model with N Maxwell element chains [48].
Figure 1. Generalized Maxwell rheological model with N Maxwell element chains [48].
Aerospace 09 00802 g001
Figure 2. Schematic Diagram of the Cohesive Interface Linear Softening Model [50,51].
Figure 2. Schematic Diagram of the Cohesive Interface Linear Softening Model [50,51].
Aerospace 09 00802 g002
Figure 3. Moisture Diffusion Analysis—Schematic for the Applied Concentration (along the thickness or z-direction).
Figure 3. Moisture Diffusion Analysis—Schematic for the Applied Concentration (along the thickness or z-direction).
Aerospace 09 00802 g003
Figure 4. Moisture Gain vs. Time Curve.
Figure 4. Moisture Gain vs. Time Curve.
Aerospace 09 00802 g004
Figure 5. Indentation Rebound: Depth vs. Time Graph.
Figure 5. Indentation Rebound: Depth vs. Time Graph.
Aerospace 09 00802 g005
Figure 6. Images of the 40J Impact Case (Dry Specimen at 25 °C) [48]. (a) View from the Impact Side. (b) View from the Back Side.
Figure 6. Images of the 40J Impact Case (Dry Specimen at 25 °C) [48]. (a) View from the Impact Side. (b) View from the Back Side.
Aerospace 09 00802 g006
Figure 7. Fitting Curves for 30, 40, and 50J Impact Cases—Dent Depth vs. Time Graph (Dry Specimen) [48].
Figure 7. Fitting Curves for 30, 40, and 50J Impact Cases—Dent Depth vs. Time Graph (Dry Specimen) [48].
Aerospace 09 00802 g007
Figure 8. Simulation Methodology: Moisture diffusion followed by LVI-induced indentation and the rebound with time: (a) moisture diffusion; (b) before impact; (c) impact dent; (d) initial dent depth; (e) final dent depth.
Figure 8. Simulation Methodology: Moisture diffusion followed by LVI-induced indentation and the rebound with time: (a) moisture diffusion; (b) before impact; (c) impact dent; (d) initial dent depth; (e) final dent depth.
Aerospace 09 00802 g008
Figure 9. Moisture Concentration Diagram after 7200 Hrs (Specimen at 25 °C/RH:100%).
Figure 9. Moisture Concentration Diagram after 7200 Hrs (Specimen at 25 °C/RH:100%).
Aerospace 09 00802 g009
Figure 10. Moisture Concentration Diagram after 7200 Hrs (Specimen at 25 °C/RH:85%).
Figure 10. Moisture Concentration Diagram after 7200 Hrs (Specimen at 25 °C/RH:85%).
Aerospace 09 00802 g010
Figure 11. Moisture Concentration vs. Thickness Curve.
Figure 11. Moisture Concentration vs. Thickness Curve.
Aerospace 09 00802 g011
Figure 12. Moisture Concentration vs. Time Curve (at thicknesses of t = 0.471 mm or t = 4.551 mm).
Figure 12. Moisture Concentration vs. Time Curve (at thicknesses of t = 0.471 mm or t = 4.551 mm).
Aerospace 09 00802 g012
Figure 13. Moisture Concentration vs. Time Curve (at ½ thickness).
Figure 13. Moisture Concentration vs. Time Curve (at ½ thickness).
Aerospace 09 00802 g013
Figure 14. The Finite Element Model [48].
Figure 14. The Finite Element Model [48].
Aerospace 09 00802 g014
Figure 15. σ 33 values (Stress distribution in the z-direction) for 50J Impact Energy Cases. (a) Dry Specimen [48]; (b) Specimen at RH:85%; (c) Specimen at RH:100%.
Figure 15. σ 33 values (Stress distribution in the z-direction) for 50J Impact Energy Cases. (a) Dry Specimen [48]; (b) Specimen at RH:85%; (c) Specimen at RH:100%.
Aerospace 09 00802 g015aAerospace 09 00802 g015b
Figure 16. σ 33 values (Stress distribution in the z-direction) for 40J Impact Energy Cases. (a) Dry Specimen [48]; (b) Specimen at RH:85%; (c) Specimen at RH:100%.
Figure 16. σ 33 values (Stress distribution in the z-direction) for 40J Impact Energy Cases. (a) Dry Specimen [48]; (b) Specimen at RH:85%; (c) Specimen at RH:100%.
Aerospace 09 00802 g016
Figure 17. σ 33 values (Stress distribution in the z-direction) for 30J Impact Energy Cases. (a) Dry Specimen [48]; (b) Specimen at RH:85%; (c) Specimen at RH:100%.
Figure 17. σ 33 values (Stress distribution in the z-direction) for 30J Impact Energy Cases. (a) Dry Specimen [48]; (b) Specimen at RH:85%; (c) Specimen at RH:100%.
Aerospace 09 00802 g017aAerospace 09 00802 g017b
Figure 18. Case-I: 50J (Dry Specimen at 25 °C).
Figure 18. Case-I: 50J (Dry Specimen at 25 °C).
Aerospace 09 00802 g018
Figure 19. Case-II: 40J (Dry Specimen at 25 °C).
Figure 19. Case-II: 40J (Dry Specimen at 25 °C).
Aerospace 09 00802 g019
Figure 20. Case-III: 30J (Dry Specimen at 25 °C).
Figure 20. Case-III: 30J (Dry Specimen at 25 °C).
Aerospace 09 00802 g020
Figure 21. Case-I: 50J (Specimen at 25 °C/RH:85%).
Figure 21. Case-I: 50J (Specimen at 25 °C/RH:85%).
Aerospace 09 00802 g021
Figure 22. Case-II: 40J (Specimen at 25 °C/RH:85%).
Figure 22. Case-II: 40J (Specimen at 25 °C/RH:85%).
Aerospace 09 00802 g022
Figure 23. Case-III: 30J (Specimen at 25 °C/RH:85%).
Figure 23. Case-III: 30J (Specimen at 25 °C/RH:85%).
Aerospace 09 00802 g023
Figure 24. Case-I: 50J (Specimen at 25 °C/RH:100%).
Figure 24. Case-I: 50J (Specimen at 25 °C/RH:100%).
Aerospace 09 00802 g024
Figure 25. Case-II: 40J (Specimen at 25 °C/RH:100%).
Figure 25. Case-II: 40J (Specimen at 25 °C/RH:100%).
Aerospace 09 00802 g025
Figure 26. Case-III: 30J (Specimen at 25 °C/RH:100%).
Figure 26. Case-III: 30J (Specimen at 25 °C/RH:100%).
Aerospace 09 00802 g026
Figure 27. Comparison Graph of Experimental and Simulation Indentation Rebound Results.
Figure 27. Comparison Graph of Experimental and Simulation Indentation Rebound Results.
Aerospace 09 00802 g027
Table 1. Damage Criterion for UD Laminae [49].
Table 1. Damage Criterion for UD Laminae [49].
Damage TypeFailure ModeDamage Initiation
Lamina Ply LevelMatrixTension Cracking ( σ y y Y T ) 2 + ( τ x y S x y ) 2 + ( τ y z S y z ) 2 1
Compression Cracking ( σ y y Y C ) 2 + ( τ x y S x y ) 2 + ( τ y z S y z ) 2 1
FiberTension Failure ( σ x x X T ) 2 + ( τ x y S x y ) 2 + ( τ y z S y z ) 2 1
Compression Failure ( σ x x X C ) 2 1
Table 2. Damage Criterion for the Cohesive Interface [50,51].
Table 2. Damage Criterion for the Cohesive Interface [50,51].
Damage TypeDamage InitiationDamage Propagation
Cohesive Layer Interface ( t n T n f ) 2 + ( t s T i f ) 2 + ( t t T i f ) 2 1
δ m = δ n 2 + δ s 2 + δ t 2
d = δ m f ( α δ m 0 ) α ( δ m f δ m 0 )
Linear Softening Mixed Mode B-K Law
σ = D · δ
D i j = { K M δ ¯ i j ,                                                       α δ m 0 δ ¯ i j [ ( 1 d ) K M + d K M δ 3 δ 3 δ ¯ i 3 ] ,   δ m 0 < α < δ m f δ ¯ i 3 δ ¯ 3 j δ 3 δ 3 K M ,                                         α δ m f  
Table 3. Material Properties of a single Unidirectional (UD) Ply.
Table 3. Material Properties of a single Unidirectional (UD) Ply.
Elastic Constants of a Single UD PlyStrength of a Single UD Ply
E1 = 144.62 GPa X T = 2612.24 MPa
E2 = 9.76 GPa X C = 1583.47 MPa
G12 = 5.44 GPa Y T = 58.25 MPa
G23 = 3.92 GPa Y C = 161.76 MPa
v12 = 0.31 S 12 = 126.79 MPa
v23 = 0.46 S 23 = 91.84 MPa
Table 4. Material Properties of the Cohesive Interface between Plies.
Table 4. Material Properties of the Cohesive Interface between Plies.
K n
(GPa/mm)
K i
(GPa/mm)
T n f
(MPa)
T i f
(MPa)
δ n f
(mm)
δ n f
(mm)
E n b
(GPa/mm)
δ n r
(mm)
1390.0510.065.595.50.0140.02526.50.043
Table 5. Experimental Results of Indentation Depth Rebound with time (Dry Specimen at 25 °C) [47].
Table 5. Experimental Results of Indentation Depth Rebound with time (Dry Specimen at 25 °C) [47].
Dry Specimen at 25 °C Temperature
Case-I: 50JCase-II: 40JCase-III: 30J
Time
(h)
Depth
(mm)
Time
(h)
Depth
(mm)
Time
(h)
Depth
(mm)
00.34800.23800.218(di)Dry
0.750.3380.50.2330.50.211
40.31110.2291.50.209
160.28840.22660.204
240.282140.22118.330.198
480.273240.21850.50.191
760.26742.50.21477.50.187
1000.264510.213101.330.185
1240.262720.210125.50.183
1480.261980.2081500.182
1720.2611250.2061740.182
1960.2601500.205198.330.181
1740.204
2200.2601980.204222.330.181
2220.204
2420.2602440.204244.330.181(df)Dry
Table 6. Experimental Results of the Indentation Depth Rebound with time (Specimen at 25 °C/RH:85%).
Table 6. Experimental Results of the Indentation Depth Rebound with time (Specimen at 25 °C/RH:85%).
At 25 °C Temperature and RH: 85%
Case-I: 50JCase-II: 40JCase-III: 30J
Time
(h)
Depth
(mm)
Time
(h)
Depth
(mm)
Time
(h)
Depth
(mm)
00.35100.24200.220 ( d i ) R H : 85 %
10.32910.23810.212
30.31020.23020.201
160.27940.22540.194
240.268180.214180.182
480.249240.211400.171
760.238480.201720.164
1000.232720.1921080.157
1240.226960.1851400.154
1480.2211200.1791660.151
1720.2191400.1741960.15
1960.2171720.1712200.149
2200.2171980.1682480.149
2440.2172200.168
2680.217
2800.2172540.1682600.149 ( d f ) R H : 85 %
Table 7. Experimental Results of the Indentation Depth Rebound with time (Specimen at 25 °C/RH:100%).
Table 7. Experimental Results of the Indentation Depth Rebound with time (Specimen at 25 °C/RH:100%).
At 25 °C Temperature and RH: 100%
Case-I: 50JCase-II: 40JCase-III: 30J
Time
(h)
Depth
(mm)
Time
(h)
Depth
(mm)
Time
(h)
Depth
(mm)
00.35300.24400.221 ( d i ) R H : 100 %
10.32510.24110.209
30.31220.23520.198
160.27640.22940.192
240.262160.216240.174
480.231240.20650.50.158
760.217480.188720.151
980.212720.1771080.145
1220.207960.1691320.143
148.50.2031200.1631660.142
170.50.2031480.159190.50.142
1980.2011720.1572240.141
2220.2011960.157
2480.2012200.1562480.141
2720.2012480.156
2960.2012720.156271.50.141 ( d f ) R H : 100 %
Table 8. Variables Values of the Fitting Curve Equation for 30, 40, and 50J Impact Cases (Dry and Conditioned Specimens).
Table 8. Variables Values of the Fitting Curve Equation for 30, 40, and 50J Impact Cases (Dry and Conditioned Specimens).
Hygrothermal ConditionImpact
Energy
A1t1A2t2y0Adj.
R-Square
25 °C, Dry Specimen [47]50J Case0.049183.483850.0391044.228150.259790.99989
40J Case0.011350.784060.0232054.841340.203550.99823
30J Case0.010740.744170.0261452.477420.180740.99644
25 °C,
RH: 85%
50J Case0.052122.382150.0822256.327800.215840.99888
40J Case0.017132.183990.0637792.102150.161900.99726
30J Case0.029782.376390.0434770.983190.147470.99808
25 °C,
RH: 100%
50J Case0.033470.762710.1180036.084390.201490.99932
40J Case0.009871.923370.0803155.139270.154520.99905
30J Case0.026001.444050.0548044.425920.140660.99880
Table 9. Applied Concentration Data from a typical Psychometric Chart.
Table 9. Applied Concentration Data from a typical Psychometric Chart.
Relative
Humidity
(% RH)
Temperature
(°C)
Pressure
(atm, mm of Hg)
Moisture
Content
(g/m3)
Mass
Concentration
(ppmv)
Volume
(% v)
85251, 76020274302.74
100251, 76023324273.24
Table 10. Impactor Velocities for different impact energies [48].
Table 10. Impactor Velocities for different impact energies [48].
Impact Energy (J)
504030
Velocity
(m/s)
4.474.003.46
Table 11. The values of α for ( G 0 ) I m p a c t .
Table 11. The values of α for ( G 0 ) I m p a c t .
Hygrothermal ConditionMaterial Typeα
25 °C/RH: 85%UD Laminae Ply0.1420 GPa
Cohesive Interface1.3622 GPa/mm
25 °C/RH: 100%UD Laminae Ply0.1880 GPa
Cohesive Interface1.8070 GPa/mm
Table 12. Summary of both types of simulation for each impact case (dry specimen at 25 °C) [48].
Table 12. Summary of both types of simulation for each impact case (dry specimen at 25 °C) [48].
Impact CaseABAQUS Explicit AnalysisABAQUS Standard Analysis
Loading Step Time
(ms)
Max Pressure
(Pa)
Initial Indentation Depth
(mm)
Rebound Step Time
(h)
Final Indentation Depth
(mm)
Case-I: 50J0.1554.8 × 1080.3482420.260
Case-II: 40J0.1194.4 × 1080.2382440.204
Case-III: 30J0.1264.1 × 1080.218244.330.181
Table 13. Summary of both types of simulation for each impact case (specimen at 25 °C/RH:85%).
Table 13. Summary of both types of simulation for each impact case (specimen at 25 °C/RH:85%).
Impact CaseABAQUS Explicit AnalysisABAQUS Standard Analysis
Loading Step Time
(ms)
Max Contact Pressure
(Pa)
Initial Indentation Depth
(mm)
Rebound Step Time
(h)
Final Indentation Depth
(mm)
Case-I: 50J0.1574.5 × 1080.3512800.217
Case-II: 40J0.1214.6 × 1080.2422540.168
Case-III: 30J0.12723.4 × 1080.2202600.149
Table 14. Summary of both types of simulation for each impact case (specimen at 25 °C/RH:100%).
Table 14. Summary of both types of simulation for each impact case (specimen at 25 °C/RH:100%).
Impact CaseABAQUS Explicit AnalysisABAQUS Standard Analysis
Loading Step Time
(ms)
Max Contact Pressure
(Pa)
Initial Indentation Depth
(mm)
Rebound Step Time
(h)
Final Indentation Depth
(mm)
Case-I: 50J0.1586.0 × 1080.3532960.201
Case-II: 40J0.1224.9 × 1080.2442720.156
Case-III: 30J0.12773.7 × 1080.221271.50.141
Table 15. Viscoelastic Material Properties for each Impact Case [48].
Table 15. Viscoelastic Material Properties for each Impact Case [48].
Impact CaseMaterial Type1st Maxwell Chain2nd Maxwell Chain
Time 1
(h)
Value 1
(GPa)
Time 2
(h)
Value 2
(GPa)
Case-I: 50JUD Laminae Ply240100142 × 107
Cohesive Interface240 K Value 1 G Value 114 K Value 2 G Value 2
55.2089.3311.68 × 1082.85 × 107
Case-II: 40JUD Laminae Ply240100212 × 107
Cohesive Interface24055.2089.331211.93 × 1083.27 × 107
Case-III: 30JUD Laminae Ply240100162 × 107
Cohesive Interface24055.2089.331162.15 × 1083.64 × 107
Table 16. Values of α and β for G 0 and G 1 .
Table 16. Values of α and β for G 0 and G 1 .
Hygrothermal ConditionRebound
Case
Material
Type
αβ
25 °C/RH:85%Case-I: 50JUD Laminae PlyThe same as in Table 110.335
Cohesive Interface
Case-II: 40JUD Laminae Ply0.580
Cohesive Interface
Case-III: 30JUD Laminae Ply0.535
Cohesive Interface
25 °C/RH:100%Case-I: 50JUD Laminae PlyThe same as in Table 110.400
Cohesive Interface
Case-II: 40JUD Laminae Ply0.630
Cohesive Interface
Case-III: 30JUD Laminae Ply0.580
Cohesive Interface
Table 17. Comparison between Experimental and Simulation Results (Dry Specimen at 25 °C) [48].
Table 17. Comparison between Experimental and Simulation Results (Dry Specimen at 25 °C) [48].
Dry Specimen at 25 °C Temperature
Impact Energy CaseTotal Dent
Rebound
Δ d D r y = ( ( d i ) D r y ( d f ) D r y )
(mm)
Way OutExperimental
Result
Simulation
Result
Prediction
Accuracy
Case-I: 50J0.088Initial and final dent depthsMatchedMatchedAccurately
predicted
Dent rebound pathThe curve is decaying at a faster rate and soon stops decaying before the
final point
The curve is decaying at a slower rate and never stops decaying until the
final point
Poor prediction
Max error: 19.35%
Case-II: 40J0.034Initial and final dent depthsMatchedMatchedAccurately
predicted
Dent rebound pathThe curve is decaying at a faster rate and soon stops decaying before the
final point
The curve is decaying at a slower rate and never stops decaying until the
final point
Fairly inaccurate prediction
Max error: 7.97%
Case-III: 30J0.037Initial and final dent depthsMatchedMatchedAccurately
predicted
Dent rebound pathThe curve is decaying at a faster rate and soon stops decaying before the
final point
The curve is decaying at a slower rate and never stops decaying until the
final point
Fairly inaccurate prediction
Max error: 9.88%
Table 18. Comparison between Experimental and Simulation Results (Specimen at 25 °C/RH:85%).
Table 18. Comparison between Experimental and Simulation Results (Specimen at 25 °C/RH:85%).
Specimen at 25 °C Temperature and RH: 85%
Impact CaseTotal Dent
Rebound
Δ d R H : 85 % = ( ( d i ) R H : 85 % ( d f ) R H : 85 % )
(mm)
Way OutExperimental
Result
Simulation
Result
Prediction
Accuracy
Case-I: 50J0.134Initial and final dent depthsMatchedMatchedAccurately
predicted
Dent rebound pathCurve is decaying at a faster rate and soon stops decaying before
the final point
Curve is decaying at a slower rate and never stops decaying until the
final point
Poor prediction
Max error: 28.36%
Case-II: 40J0.074Initial and final dent depthsMatchedMatchedAccurately
predicted
Dent rebound pathCurve is decaying at a faster rate and soon stops decaying before
the final point
Curve is decaying at a slower rate and never stops decaying until the
final point
Fairly inaccurate prediction
Max error: 12.29%
Case-III: 30J0.071Initial and final dent depthsMatchedMatchedAccurately
predicted
Dent rebound pathCurve is decaying at a faster rate and soon stops decaying before
the final point
Curve is decaying at a slower rate and never stops decaying until the
final point
Inaccurate
prediction
Max error: 20.36%
Table 19. Comparison between Experimental and Simulation Results [Specimen at 25 °C/RH:100%].
Table 19. Comparison between Experimental and Simulation Results [Specimen at 25 °C/RH:100%].
Specimen at 25 °C Temperature and RH: 100%
Impact CaseTotal Dent
Rebound
Δ d R H : 100 % = ( ( d i ) R H : 100 % ( d f ) R H : 100 % )
(mm)
Way OutExperimental
Result
Simulation
Result
Prediction
Accuracy
Case-I: 50J0.152Initial and final dent depthsMatchedMatchedAccurately
predicted
Dent rebound pathCurve is decaying at a faster rate and soon stops decaying before the
final point
Curve is decaying at a slower rate and never stops decaying until the
final point
Poor prediction
Max error: 37.90%
Case-II: 40J0.088Initial and final dent depthsMatchedMatchedAccurately
predicted
Dent rebound pathCurve is decaying at a faster rate and soon stops decaying before the
final point
Curve is decaying at a slower rate and never stops decaying until the
final point
Inaccurate
prediction
Max error: 21.50%
Case-III: 30J0.080Initial and final dent depthsMatchedMatchedAccurately
predicted
Dent rebound pathCurve is decaying at a faster rate and soon stops decaying before the
final point
Curve is decaying at a slower rate and never stops decaying until the
final point
Poor prediction
Max error: 28.04%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yousaf, M.; Zhou, C.; Yang, Y.; Wang, L. Numerical Study of the Hygrothermal Effects on Low Velocity Impact Induced Indentation and Its Rebound in Composite Laminate. Aerospace 2022, 9, 802. https://doi.org/10.3390/aerospace9120802

AMA Style

Yousaf M, Zhou C, Yang Y, Wang L. Numerical Study of the Hygrothermal Effects on Low Velocity Impact Induced Indentation and Its Rebound in Composite Laminate. Aerospace. 2022; 9(12):802. https://doi.org/10.3390/aerospace9120802

Chicago/Turabian Style

Yousaf, Muhammad, Chuwei Zhou, Yu Yang, and Li Wang. 2022. "Numerical Study of the Hygrothermal Effects on Low Velocity Impact Induced Indentation and Its Rebound in Composite Laminate" Aerospace 9, no. 12: 802. https://doi.org/10.3390/aerospace9120802

APA Style

Yousaf, M., Zhou, C., Yang, Y., & Wang, L. (2022). Numerical Study of the Hygrothermal Effects on Low Velocity Impact Induced Indentation and Its Rebound in Composite Laminate. Aerospace, 9(12), 802. https://doi.org/10.3390/aerospace9120802

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop