Next Article in Journal
DFT Study on the Enhancement of Isobaric Specific Heat of GaN and InN Nanosheets for Use as Nanofluids in Solar Energy Plants
Next Article in Special Issue
Fatigue Behaviour and Its Effect on the Residual Strength of Long-Fibre-Reinforced Thermoplastic PP LGF30
Previous Article in Journal
Surface and Structural Characterization of PVTMS Films Treated by Elemental Fluorine in Liquid Perfluorodecalin
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental and Numerical Study of the Influence of Pre-Existing Impact Damage on the Low-Velocity Impact Response of CFRP Panels

1
Politecnico di Milano, Dipartimento di Meccanica, Via La Masa 1, 20156 Milano, Italy
2
PPGE3M, Federal University of Rio Grande do Sul, Porto Alegre 91501-970, Brazil
*
Author to whom correspondence should be addressed.
Materials 2023, 16(3), 914; https://doi.org/10.3390/ma16030914
Submission received: 17 December 2022 / Revised: 3 January 2023 / Accepted: 16 January 2023 / Published: 18 January 2023

Abstract

:
This paper presents an experimental and numerical investigation on the influence of pre-existing impact damage on the low-velocity impact response of Carbon Fiber Reinforced Polymer (CFRP). A continuum damage mechanics-based material model was developed by defining a user-defined material model in Abaqus/Explicit. The model employed the action plane strength of Puck for the damage initiation criterion together with a strain-based progressive damage model. Initial finite element simulations at the single-element level demonstrated the validity and capability of the damage model. More complex models were used to simulate tensile specimens, coupon specimens, and skin panels subjected to low-velocity impacts, being validated against experimental data at each stage. The effect of non-central impact location showed higher impact peak forces and bigger damage areas for impacts closer to panel boundaries. The presence of pre-existing damage close to the impact region leading to interfering delamination areas produced severe changes in the mechanical response, lowering the impact resistance on the panel for the second impact, while for non-interfering impacts, the results of the second impact were similar to the impact of a pristine specimen.

1. Introduction

The response of composite materials to out-of-plane impact loads during their service life is of concern in many applications. The low-velocity impact may cause interlaminar and intralaminar failure modes, such as matrix or fiber failure in tension and compression, crushing, and delamination. Matrix failure and delamination are particularly important since they can be invisible to the naked eye and still cause severe material degradation [1,2,3]. The low-velocity impact response of laminates has been the focus of many recent studies, especially related to a single hit [4,5,6,7,8]. Although repeated impacts are not uncommon in real-world applications of composites [9,10], there are not many studies in the literature addressing that, either numerically- or experimentally-based, as most focus on multiple hits at the same point [11]. The studies on the repeated impact behavior of composites mostly concern repeated impact at the same location by considering different impactor or target features, such as impactor geometry [12,13,14] and mass [15], reinforcement type [16,17,18,19], laminate thickness [17,20], and stacking sequence [21,22]. Atas et al. [23] investigated the effect of thickness on the repeated low-velocity impact response of E-glass/epoxy composites reporting the energy profiling diagram and the perforation limit. Liao et al. [13] performed an experimental study on damage accumulation of a repeated low-velocity impact in CFRP laminates considering different diameters of the hemispherical impactor. They presented a damage index based on stiffness reduction caused by repeated impacts which correlated well with the damage observed by a C-scan ultrasound.
Numerical simulations have proven to be an effective tool for the prediction of repeated low-velocity impacts on composites, despite the complexity of the loading condition, the constitutive model, and the multi-failure modes of composites [24]. For instance, Zhou et al. [25] developed a numerical model based on Hashin failure criteria to study CFRP laminates and observed a reduction in energy absorption for the second impact. Additionally, in a previous publication of the authors [26], the effect of repeated impacts on the initiation and propagation of interlaminar and intralaminar damage in CFRP laminates with each impact was investigated. By developing a model based on the Puck failure criterion with a constant fracture angle, the accumulation of damage was assessed. The CFRP laminate was impacted five times at the center, and the presence of damage from the first impact led to an increase in peak impact force [27]. High-fidelity macro-mechanical models have been able to predict intralaminar and interlaminar damage and permanent indentation [4], and also matrix damage and delamination [28,29]. The Puck inter-fiber failure criterion has been accurately applied to CFRP [29,30], as in the work of [31], who developed a strain-based progressive damage model with IFF of Puck to simulate laminates subjected to low-velocity impacts. Progressive delamination in a woven composite subjected to repeated impacts was simulated in [32] using a cohesive zone model, showing that delamination propagates rapidly after the first impact.
Regarding the repeated impacts at different locations, only a few experimental studies can be found [33,34], despite the practical engineering interest. Nassir et al. [35] studied multiple impacts in plain-weave glass fiber/epoxy laminates and showed that the critical force for the generation of damage does not change for a second hit at a different location. In addition, Liao et al. [36] investigated double impacts in CFRP laminates and used four different distances for the sequential impacts, concluding that the interference between the first and the second impact damage areas has a pivotal role in the response of the second impact. For impacts with interfering damaged areas, lower energy absorption was observed in the second impact, while the opposite phenomenon was observed for non-interfering damaged areas. More recently, Huang et al. [37] showed that for repeated impacts at different locations, the maximum displacement better characterizes the damage interference compared to the bending stiffness.
Considering the cited work, it is clear that further studies are required, especially numerically, to understand the complex damage behavior of composites subjected to multiple impacts. Additionally, there is a lack of studies on repeated impacts at different locations even though they are more likely to occur compared to repeated impacts at the same location. Thus, the aim of the current study is to numerically investigate the effect of the damage caused by a low-velocity impact on the structural integrity and impact performance of CFRP panels subjected to a subsequent hit at different locations, focusing on their interaction.

2. Experimental Study

The whole experimental campaign included three different levels of experimental tests on CFRP specimens obtained from unidirectional epoxy-based prepregs (CYCOM® 977-2-34%-24K IMS-196-T1, Solvay, Brussels, Belgium), namely: (i) tensile and shear tests, (ii) low-velocity impact tests on coupon specimens (150 mm × 100 mm), and (iii) single-hit and multiple-hit low-velocity impact tests on skin panels (340 mm × 210 mm). The tensile and shear tests have already been thoroughly reported in a previous reference [29] and will not be discussed here.
Low-velocity impact tests were performed on coupon level samples (dimensions: 150 mm × 100 mm) according to ASTM D7136 [38]. Laminates with 24 plies (thickness: 4.8 mm) were manufactured with [45/0/−45/90]3s layup and tested at impact energies of 8 J, 10 J, and 12 J with a drop weight equipment using a hemispherical impactor (diameter: 16 mm). The force–time histories were recorded using a load cell (Kistler 9331B, Winterthur, Switzerland), an amplifier (Kistler 5011B, Winterthur, Switzerland), and NI 9239 data acquisition system (National Instruments Corp, Austin, TX, USA). Besides, skin panels (dimensions: 400 mm × 270 mm × 3.42 mm) produced with 18 plies and a stacking sequence of [±45/0/±45/∓45/90/0]S were impacted first at the center, with an energy of 35 J or 40 J, followed by off-central impacts, with an impact energy of 25 J. Each of these panels experienced up to two low-velocity impacts at different locations according to Figure 1.

3. Numerical Simulation

3.1. Material Constitutive and Damage Model

Puck inter-fiber failure criterion [39,40] for unidirectional composites was considered for the CFRP plies. Puck failure equations (Equations (3) and (4)) calculate the stress in tension ( σ n θ > 0 ) and compression ( σ n θ < 0 ), which are used to compute the index ranging from 0 to 1, for unloaded condition to failure initiation, respectively.
F F 1 = σ 11 R t
F F 2 = σ 11 R c
I F F 1 θ = 1 R t P Ψ t R Ψ A σ n 2 + τ n t R A 2 + τ n l R 2 + P Ψ t R Ψ A σ n   f o r   σ n 0
I F F 2 θ = τ n t R A 2 + τ n l R 2 + P Ψ c R Ψ A σ n 2 + P Ψ c R Ψ A σ n   f o r   σ n < 0
where θ is the angle of rotation of the coordinate system around the x 1 axis, σ n θ , τ n t θ , τ n 1 θ are the stresses on the action plane, and R t and R c are the material strength in the fiber direction in tension and compression, respectively. R t are R are the material strengths in normal and shear directions, respectively. While:
P Ψ t R Ψ A = P t R A c o s 2 Ψ + P t R s i n 2 Ψ
P Ψ c R Ψ A = P c R A c o s 2 Ψ + P c R s i n 2 Ψ
and:
c o s 2 Ψ = 1 s i n 2 Ψ = τ n t 2 τ n t 2 + τ n l 2
R A = R c 2 1 + P c
where P t , P t , P c , and P c are the inclination parameters for the master fracture body.
As Puck’s IFF requires the determination of the fracture plane with the maximum exposure value, an algorithm should be included during the material implementation to search for that plane [29]. The Simple Parabolic Interpolation Search (SPIS), a non-iterative algorithm for the search of a fracture angle without the need for iterative methods, has been used for each element during the explicit simulations. The SPIS is described in detail in a previous publication of the group [41].
The non-linear behavior was considered here using a third-degree polynomial (Equation (9)) [4].
τ 12 = c 1 γ 12 3 + s i g n γ 12 · c 2 γ 12 2 + c 3 γ 12
where the coefficients c 1 , c 2 , and c 3 were obtained by fitting Equation (9) to a stress–strain curve from the in-plane shear test. Since the non-linear shear strain is irreversible, the total strain is separated into elastic and plastic strain [42]:
γ 12 = γ 12 e + γ 12 p
The irreversible plastic strain, γ 12 p , is calculated as:
γ 12 p = γ 12 τ 12 G 12 0
where G 12 0 is the initial shear modulus, and the unloading takes place according to G 12 0 .
After initiation of fiber or inter-fiber damage, the material is degraded according to a strain-based post-damage law. The degradation is controlled by defining several damage variables based on strain softening. The stiffness parameters are degraded selectively, based on the failure mode that occurred. The relationship between the damage variables and the stiffness parameters is shown in the damage matrix (Equation (12)):
d C = ( d f E 1 ( 1 d m v 23 v 32 ) ) / ( d f E 1 ( d m v 2 1 d m v 31 v 23 ) ) / ( d f E 1 ( v 3 1 d m v 21 v 32 ) ) / 0 0 0 ( d m E 2 ( 1 d f v 13 v 31 ) ) / ( d m E 2 ( v 3 2 d f v 12 v 31 ) ) / 0 0 0 ( E 3 ( 1 d m d f v 12 v 21 ) ) / 0 0 0 d m G 12 0 0 0 s y m m e t r i c d m G 13 0 d m G 23 for   Δ = 1 d f d m v 12 v 21 d m v 23 v 32 d f v 31 v 13 2   d f d m v 21 v 32 v 13
where d f and d m are the damage variables for the fiber and matrix, respectively, defined by Equation (13).
d f = 1 d F F 1 1 d F F 2 d m = 1 d I F F 1 1 s m c · d I F F 2
and s m c is equal to 0.6, representing a lower limit for the matrix damage variable to avoid numerical instability and element distortion. d F F 1 , d F F 2 , d I F F 1 , and d I F F 2 are the damage variables corresponding to the failure Equations (1)–(4). These variables follow the strain-softening Equation (14) and are greater than 0 when the corresponding failure equation is met.
d i = ε i , u ε i , u ε i , o n 1 ε i , o n ε i
where the subindex i represents each of the failure modes; ε i , u is the ultimate failure strain, calculated by 2 G i S i l , and G i and S i are the fracture energy and failure stress for each mode, respectively; l is the element characteristic length calculated from the cube root of the element volume at the beginning of the simulation, included in the simulations according to the crack-band model [43] to avoid mesh sensitivity problem; ε i is the equivalent strain input from simulation, being equal to ε 1 , < ε 1 > , ε n 2 + γ n t 2 + γ n l 2 , ε n 2 + γ n t 2 + γ n l 2 for d F F 1 , d F F 2 , d I F F 1 , and d I F F 2 , respectively; ε i , o n is the failure onset strain, calculated when the failure criterion is met.

3.2. Damage Model for the Interface

Here, the damage on the interface is based on the bi-linear traction-separation law which uses a quadradic stress failure criterion (Equation (15)) to predict damage initiation and the B-K mixed model fracture energy law (Equation (16)) [44] to predict delamination propagation.
t n 2 N 2 + t s 2 S 2 + t t 2 S 2 = 1
where <·> is the Macaulay bracket to only consider the contribution of tensile normal traction at the interface; and N and S are the normal and shear strengths of the interface, respectively.
G n C + G s C G n C G s G T η = G C
where G n C and G s C are the critical normal and shear fracture energies, respectively; η is the B-K material constant; and G S = G s + G t and G T = G n + G S .

3.3. Finite Element Models

The finite element models were developed in four levels of complexity. The simplest simulations were performed on single elements subjected to tensile, compression, and shear cyclic loads, designed to verify the accuracy of the progressive damage model, particularly related to damage initiation, fracture energy, and degradation of properties. Tensile specimens were simulated to verify the interaction among different failure models and to enable direct comparison with experimental results. The results for these two levels have already been thoroughly reported in a previous reference [29] and will not be discussed here.
For the low-velocity impact model on coupon specimens of the 24-ply laminate (Figure 2a), a detailed layer-by-layer description of inter-laminar and intra-laminar damage and a comparison with experimental data is provided for the impacts at 8 J, 10 J, or 12 J. This is the same model used in [26].
The model was discretized with C3D8R (8-node linear brick with reduced integration) elements with enhanced hourglass control. The element size was defined as 1 mm × 1 mm at the impactor/laminate contact region, with a coarser mesh outside this region, following a mesh sensitivity analysis which found that suitable for this material and loading condition [26]. The interface between plies in the laminate was simulated using layers of the zero-thickness cohesive element (COH3D8) whose properties are compiled in Table 1.
Since the skin panels experienced up to two impacts at different locations, see Figure 1, two different impactors were defined at each impact location moving towards the panel with predefined initial velocities. The second impact was only applied after the debouncing of the previous impactor. The mechanical properties of the CFRP composites used in the simulation are presented in Table 1, and they were obtained either from experiments or from similar materials in the literature.
For the simulation of multiple impacts, illustrated in Figure 2b, the material state of the laminate at the end of each loading step becomes the initial material state for the next loading step. This should include element stress and strain state, displacement and velocity, and damage indicators for the plies and interfaces. An initial predefined field option of Abaqus/Explicit [47] was used, as in [26]. For multiple impacts at the same location, a multi-step approach is found in the literature [25,26], i.e., applying the low-velocity impact in one simulation and then moving the deformed geometry to a new simulation and applying the subsequent impact. However, for the simulation of multiple impacts at different locations, as in the current work, the panel was hit by the impactors at the same simulation step.
This approach is considered valid if enough time passes after each impact event to mitigate unwanted oscillation. To investigate that, the force–time curve obtained for the second impact using two different simulation schemes is compared in Figure 3, that is, the multi-step simulation approach of hitting the panel once in each step, and the single-step approach with all impacts in one step. Since the response of the panel from both models is identical, the rest step was considered unnecessary for the simulation of low-velocity impacts on CFRP composite, as previously reported in [26]. Based on that, the multiple impacts were simulated in a single-step explicit FE model.

3.4. Parametric Study

Different impact locations and impact energies were simulated. These are identifiable according to the notation where the first term represents the impact event type (S for single impact and D for double impact), the second the impact location according to Figure 1 (i.e., C, S1, S2, S3, S4), and the third term the impact energy (25 J, 35 J, 40 J). For instance, S-S1-25J means a single 25 J impact at the S1 location in Figure 1. For the double impact cases, the first impact always occurs at the center of the skin panel at 40 J and the subsequent impact happens at one of the four different off-central locations with three different impact energies. So, only the second impact differs in these cases, and the notation was used to represent that second impact. A summary of all the parametric study cases is presented in Table 2.

4. Results and Discussion

In this section, the experimental and numerical results of the coupon under a low-velocity impact and skin panels subjected to single and multiple impacts at different locations are discussed.

4.1. Verification of the FE Simulations

The experimental and numerical curves force–displacement curves for the 8 J and 12 J impacts are compared in Figure 4. The numerical model provided accurate predictions of peak force, absorbed energy, peak displacement, and impact time duration.
The absorbed energies for the CFRP laminates subjected to impact at different energies are compared in Table 3. The underestimation in energy absorption can be attributed to the underestimation in damage by the numerical model. The numerically observed delamination area was compared to the projected delamination area measured using the Ultrasonic C-scan technique in Table 3. The error in the delamination area was calculated as less than −7% for all impact energies.
Figure 5 shows the layer-by-layer matrix damage predicted by the progressive damage model with the Puck failure criterion. The layer-by-layer analysis of matrix damage shows that due to the presence of dominant shear stresses at the plies in the middle of the laminate, the prediction of the damaged area did not follow the fiber orientation at each layer. The bending deformation in thin laminates results in high tensile bending stress at the rear side of the specimen, which is the location of tensile matrix damage initiation [29]. The matrix tensile damage propagated in-plane and through the thickness of the laminate until the rebounding of the impactor and no further propagation was observed during the rebounding. On the other hand, the compressive matrix damage that initiates in the contact point of the impactor and the laminate propagates towards the rear side. The dominant out-of-plane shear stresses in the middle plies result in the misalignment of the damaged area with the fiber orientation in the plies. Different damage patterns at the top, middle, and bottom layers were also reported in [48]. Matrix compression damage was only found in the layers near the impact side and the area under the impactor tip, while matrix tensile damage was greater in the bottom layers due to the bending effect. A similar observation was reported by [49,50] for different failure criteria such as Puck and Hashin.

4.2. Low-Velocity Impact on Skin Panels

Figure 6a–i compares the experimental and numerical curves of force–time, force–displacement, and energy–time for the CFRP skin panels subjected to a single impact at 40 J or 35 J at the center, and to a second impact off-center at 25 J (after the first impact at 40 J at the center). The model was capable of accurately predicting the response of the panel subjected to impacts at the center or off-center. The values of impact time duration, impact peak force, maximum displacement, and absorbed energy were consistent with the corresponding experimental data. Additionally, the numerical model accurately predicted the drops in the force–time curves at different energies, attributed to the onset and propagation of the damage in the panel, with a detrimental effect on the mechanical properties, including stiffness. More fluctuations and force drop due to the presence of fiber failure are seen in comparison to the coupon specimens subjected to 8, 10, and 12 J impacts.

4.3. Effect of the Impact Location

A total of 15 low-velocity impact cases were considered, three of them at the center and the others off-central, i.e., at S1, S2, S3, or S4 (see Figure 1). Figure 7a,b compares the force–time and force–displacement curves of the skin panels subjected to a single 40 J impact and an increase in peak force and a decrease in impact time duration is observed when the impact location shifts from the center of the plate towards the boundary. Changing the impact location from S3 to S4, which is the closest to the plate’s boundary, led to a more pronounced change in the response. A significant increase in impact bending stiffness of the panel can be observed in Figure 7b when the impact location is closer to the panel boundaries, and the maximum displacement decreased in the off-center impacts being the lowest for location S4. The energy–time curve in Figure 7c shows that shifting the impact location towards the clamped edge increased total absorbed energy. This was also observed for the 25 J and 40 J impacts. A similar finding was reported by [35] for woven glass fiber-reinforced epoxy subjected to center and off-center impacts.
Figure 8a–d summarizes the response parameters for all locations. The ratio of absorbed energy to impact energy is higher for off-center locations, which is consistent with the increase in the projected delamination area (Figure 8d). This has also been observed in the experimental findings of [35]. Furthermore, for the 35 J impact, the absorbed energy for the central impact was 9.02 J, which increased by 49% to 14.09 J for the off-center impact at S4 location. Accordingly, a 33.7% increase in the delamination area was observed (Figure 8d).

4.4. Effect of Pre-Existing Impact Damage

Figure 9 compares the force–time, force–displacement, and energy–time curves between a pristine panel and a panel with pre-existing damage subjected to impact at 40 J at off-center locations. Figure 9a compares the force–time responses for impact at S1, showing significant changes in the impact response of the pristine specimen compared to the specimen with pre-existing damage. The higher peak force and maximum impact displacement for the repeated impact in Figure 9b were also reported in [26,37]. Figure 9c shows that the interference of the damaged area from the first and second impact significantly reduces the energy absorption of the specimen [26]. Similarly, the force–time/displacement response in Figure 9d,e for impacts at S2 shows higher peak force and maximum displacement for impact on the panel with pre-existing damage; however, it can be observed that the trends are mitigated compared to the impact at S1. Lower absorbed energy for repeated impact at S2 was observed in Figure 9f, which indicated the interference of damaged areas. The impact force–time, force–displacement, and energy–time responses at S3 are shown in Figure 9g–i, respectively. At this location, the pre-existing impact damage shows a less severe influence on the response. It can be observed from Figure 9j–l that for the impact at S4, the responses are identical, while for the impact at S1 the panel with pre-existing damage, showed higher peak force and displacement, and lower absorbed energy (see Figure 9a–h), as previously reported for CFRP panels subjected to repeated impacts [25,26].
The greater maximum displacement for double impacts due to the change in bending stiffness was discussed by [36], while here the bending stiffness reduction of the first impact did not influence the impact response for far-away impact locations, see Figure 10b. The low-velocity impact in CFRP materials can thus be considered to cause local damage. However, the extent of the local damage needs to be investigated since it depends on different factors, especially the impact energy.
To find the minimum distance at which the pre-existing damage becomes relevant, the responses to the impact at 25 J and 40 J of a pristine panel and a panel with pre-existing damage are compared in Figure 10. It can be observed that their mechanical response was very similar far from the center for both impact energies. The 25 J impact shows the same result also for S2 and S3 locations, while the pre-existing damage area was the same for all cases of this figure.
Figure 10c,d compares the energy absorption ratio and total delamination area for the panels subjected to a single hit and double hits. As before, for panels with pre-existing damage, a significant difference can be observed when the impact happens in the vicinity of the pre-existing damage. For instance, for the second impact at 40 J, the presence of the pre-existing damage led to a 22.4%, 20.1%, and 3.8% decrease in the absorbed energy compared to the pristine panel for impact locations S1, S2, and S3, respectively. This is consistent with the numerical results of [25,26] for repeated impacts at the same location.
The damage accumulation pattern shows a smaller increase from the first to the second impact for the S1 location, while for impacts further away from the pre-existing damage, the delaminated area at the end of the second impact nearly doubled compared to the first impact. This is due to the change in the energy absorption mechanism of the CFRP specimen with pre-existing damage. Since the damage is already present in the vicinity of the contact area, the panel deforms more, and the peak impact force increases to absorb the initial impact energy, which is later released back to the impactor as elastic deformation energy. Therefore, the energy absorption capacity of the panel is decreased [26].
A comparison of the propagation of damage for different repeated impact scenarios on skin panels is shown in Figure 11. In this case, the projected damaged areas of D-S1-40 J and D-S4-40 J models after the first (at center) and second (off-center) impacts are compared. It can be seen that for the first impact, the damaged region is similar in shape and area since in both models, the panel was subjected to a 40 J impact at the center. The damaged area for the impact at the center is compared with the c-scan measurement to demonstrate the accuracy of the model in predicting the damage. Significant changes in the damaged area can be observed after the second impact. Two district and non-interfering damaged regions were observed for the D-S4-40 J model, while in D-S1-40 J, there was interference between damaged regions resulting in a single region of damaged material. This was consistent with the results of Figure 9, where the second impact at S1 resulted in bigger displacements and lower energy absorption. It is also worth mentioning that for the second impact at S4, the damage propagated towards the boundaries of the panel instead of the pre-existing damaged region; this behavior was similar to the single impact at the same location, meaning the pre-existing damage area did not inflict a significant change in the damage propagation of the off-center impact.

5. Conclusions

In this paper, a continuum damage mechanics-based material model was used to investigate the mechanical response and damage in CFRP composites subjected to low-velocity impacts. The progressive damage model was implemented in Abaqus/Explicit and validated in different stages, from coupon specimens to skin panels subjected to low-velocity impacts, where the constitutive model based on the Puck failure criterion led to accurate results.
The analysis showed that the impact location has an important effect on the mechanical response and damage of composite skin panels due to the significant increase in impact bending stiffness when the impact location moves towards the panel boundaries. The presence of pre-existing impact damage of 805 mm2 from a 40 J impact at the center of the skin panel resulted in a more complex impact response of the panel that depended on the impact location with respect to the pre-existing damage area. For an impact far away from the previous impact, the panel shows no difference in response compared to a similar impact on the pristine specimen. For an impact in the vicinity of a pre-existing damaged area, there was a significant change in the response of the panel compared to that of a pristine specimen. For the former, greater impact peak force and displacement and lower energy absorption capacity were observed.

Author Contributions

Conceptualization, M.R. and A.M.; methodology, M.R.; software, M.R.; validation, M.R.; investigation, M.R.; resources, A.B. and C.S.; data curation, M.R., A.B. and C.S.; writing—original draft preparation, M.R.; writing—review and editing, A.M., A.B. and S.C.A.; visualization, M.R.; supervision, A.M.; project administration, A.M.; funding acquisition, A.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sławski, S.; Szymiczek, M.; Kaczmarczyk, J.; Domin, J.; Świtoński, E. Low Velocity Impact Response and Tensile Strength of Epoxy Composites with Different Reinforcing Materials. Materials 2020, 13, 3059. [Google Scholar] [CrossRef] [PubMed]
  2. Cao, H.; Ma, M.; Jiang, M.; Sun, L.; Zhang, L.; Jia, L.; Tian, A.; Liang, J. Experimental Investigation of Impactor Diameter Effect on Low-Velocity Impact Response of CFRP Laminates in a Drop-Weight Impact Event. Materials 2020, 13, 4131. [Google Scholar] [CrossRef] [PubMed]
  3. Rezasefat, M.; Gonzalez-Jimenez, A.; Ma, D.; Vescovini, A.; Lomazzi, L.; da Silva, A.A.X.; Amico, S.C.; Manes, A. Experimental Study on the Low-Velocity Impact Response of Inter-Ply S2-Glass/Aramid Woven Fabric Hybrid Laminates. Thin Walled Struct. 2022, 177, 109458. [Google Scholar] [CrossRef]
  4. Tan, W.; Falzon, B.G.; Chiu, L.N.S.; Price, M. Predicting Low Velocity Impact Damage and Compression-After-Impact (CAI) Behaviour of Composite Laminates. Compos. Part A Appl. Sci. Manuf. 2015, 71, 212–226. [Google Scholar] [CrossRef] [Green Version]
  5. Bogenfeld, R.; Kreikemeier, J.; Wille, T. Review and Benchmark Study on the Analysis of Low-Velocity Impact on Composite Laminates. Eng. Fail. Anal. 2018, 86, 72–99. [Google Scholar] [CrossRef]
  6. Li, N.; Chen, P.H. Micro-Macro FE Modeling of Damage Evolution in Laminated Composite Plates Subjected to Low Velocity Impact. Compos. Struct. 2016, 147, 111–121. [Google Scholar] [CrossRef] [Green Version]
  7. Rezasefat, M.; Amico, S.C.; Giglio, M.; Manes, A. A Numerical Study on the Influence of Strain Rate in Finite-Discrete Element Simulation of the Perforation Behaviour of Woven Composites. Polymers 2022, 14, 4279. [Google Scholar] [CrossRef]
  8. Singh, H.; Mahajan, P. Modeling Damage Induced Plasticity for Low Velocity Impact Simulation of Three Dimensional Fiber Reinforced Composite. Compos. Struct. 2015, 131, 290–303. [Google Scholar] [CrossRef]
  9. Balasubramaniam, K.; Ziaja, D.; Jurek, M.; Fiborek, P.; Malinowski, P. Experimental and Numerical Analysis of Multiple Low-Velocity Impact Damages in a Glass Fibered Composite Structure. Materials 2021, 14, 7268. [Google Scholar] [CrossRef]
  10. Lee, S.E.; Kim, D.U.; Cho, Y.J.; Seo, H.S. Multiple Impact Damage in GLARE Laminates: Experiments and Simulations. Materials 2021, 14, 7800. [Google Scholar] [CrossRef]
  11. Sadighi, M.; Alderliesten, R. Impact Fatigue, Multiple and Repeated Low-Velocity Impacts on FRP Composites: A Review. Compos. Struct. 2022, 297, 115962. [Google Scholar] [CrossRef]
  12. Katunin, A.; Pawlak, S.; Wronkowicz-Katunin, A.; Tutajewicz, D. Damage Progression in Fibre Reinforced Polymer Composites Subjected to Low-Velocity Repeated Impact Loading. Compos. Struct. 2020, 252, 112735. [Google Scholar] [CrossRef]
  13. Liao, B.; Zhou, J.; Li, Y.; Wang, P.; Xi, L.; Gao, R.; Bo, K.; Fang, D. Damage Accumulation Mechanism of Composite Laminates Subjected to Repeated Low Velocity Impacts. Int. J. Mech. Sci. 2020, 182, 105783. [Google Scholar] [CrossRef]
  14. Dogan, A. Single and Repeated Low-Velocity Impact Response of E-Glass Fiber-Reinforced Epoxy and Polypropylene Composites for Different Impactor Shapes. J. Thermoplast. Compos. Mater. 2022, 35, 320–336. [Google Scholar] [CrossRef]
  15. Sugun, B.S.; Rao, R.M.V.G.K. Impactor Mass Effects in Glass–Epoxy Composites Subjected to Repeated Drop Tests. J. Reinf. Plast. Compos. 2016, 23, 1547–1560. [Google Scholar] [CrossRef]
  16. Sugun, B.S.; Rao, R.M.V.G.K. Low-Velocity Impact Characterization of Glass, Carbon and Kevlar Composites Using Repeated Drop Tests. J. Reinf. Plast. Compos. 2004, 23, 1583–1599. [Google Scholar] [CrossRef]
  17. de Morais, W.A.; Monteiro, S.N.; d’Almeida, J.R.M. Effect of the Laminate Thickness on the Composite Strength to Repeated Low Energy Impacts. Compos. Struct. 2005, 70, 223–228. [Google Scholar] [CrossRef]
  18. Saleh, M.N.; El-Dessouky, H.M.; Saeedifar, M.; De Freitas, S.T.; Scaife, R.J.; Zarouchas, D. Compression after Multiple Low Velocity Impacts of NCF, 2D and 3D Woven Composites. Compos. Part A Appl. Sci. Manuf. 2019, 125, 105576. [Google Scholar] [CrossRef]
  19. Donadio, F.; Ferraro, P.; Lopresto, V.; Pagliarulo, V.; Papa, I.; Rippa, M.; Russo, P. Response of Glass/Carbon Hybrid Composites Subjected to Repeated Low Velocity Impacts. J. Compos. Mater. 2022, 56, 2789–2802. [Google Scholar] [CrossRef]
  20. Boria, S.; Scattina, A.; Belingardi, G. Experimental Investigation on a Fully Thermoplastic Composite Subjected to Repeated Impacts. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2019, 233, 6985–7002. [Google Scholar] [CrossRef]
  21. Icten, B.M. Repeated Impact Behavior of Glass/Epoxy Laminates. Polym. Compos. 2009, 30, 1562–1569. [Google Scholar] [CrossRef]
  22. David-West, O.S.; Nash, D.H.; Banks, W.M. An Experimental Study of Damage Accumulation in Balanced CFRP Laminates due to Repeated Impact. Compos. Struct. 2008, 83, 247–258. [Google Scholar] [CrossRef]
  23. Atas, C.; Icten, B.M.; Küçük, M. Thickness Effect on Repeated Impact Response of Woven Fabric Composite Plates. Compos. Part B Eng. 2013, 49, 80–85. [Google Scholar] [CrossRef]
  24. Rezasefat, M.; Giglio, M.; Manes, A. Numerical Investigation of the Effect of Open Holes on the Impact Response of CFRP Laminates. Appl. Compos. Mater. 2022, 29, 1555–1578. [Google Scholar] [CrossRef]
  25. Zhou, J.; Wen, P.; Wang, S. Numerical Investigation on the Repeated Low-Velocity Impact Behavior of Composite Laminates. Compos. Part B Eng. 2020, 185, 107771. [Google Scholar] [CrossRef]
  26. Rezasefat, M.; Gonzalez-Jimenez, A.; Giglio, M.; Manes, A. Numerical Study on the Dynamic Progressive Failure due to Low-Velocity Repeated Impacts in Thin CFRP Laminated Composite Plates. Thin Walled Struct. 2021, 167, 108220. [Google Scholar] [CrossRef]
  27. Zhou, J.; Liu, B.; Wang, S. Finite Element Analysis on Impact Response and Damage Mechanism of Composite Laminates under Single and Repeated Low-Velocity Impact. Aerosp. Sci. Technol. 2022, 129, 107810. [Google Scholar] [CrossRef]
  28. Liao, B.B.; Liu, P.F. Finite Element Analysis of Dynamic Progressive Failure of Plastic Composite Laminates under Low Velocity Impact. Compos. Struct. 2017, 159, 567–578. [Google Scholar] [CrossRef]
  29. Rezasefat, M.; Gonzalez-Jimenez, A.; Giglio, M.; Manes, A. An Evaluation of Cuntze and Puck Inter Fibre Failure Criteria in Simulation of Thin CFRP Plates Subjected to Low Velocity Impact. Compos. Struct. 2021, 114654. [Google Scholar] [CrossRef]
  30. Hale, P.; Ng, E. gene Non-Linear Material Characterization of CFRP with FEM Utilizing Cohesive Surface Considerations Validated with Effective Tensile Test Fixturing. Mater. Today Commun. 2020, 23, 100872. [Google Scholar] [CrossRef]
  31. Donadon, M.V.; Iannucci, L.; Falzon, B.G.; Hodgkinson, J.M.; de Almeida, S.F.M. A Progressive Failure Model for Composite Laminates Subjected to Low Velocity Impact Damage. Comput. Struct. 2008, 86, 1232–1252. [Google Scholar] [CrossRef]
  32. Wang, C.; Chen, Z.; Silberschmidt, V.V.; Roy, A. Damage Accumulation in Braided Textiles-Reinforced Composites under Repeated Impacts: Experimental and Numerical Studies. Compos. Struct. 2018, 204, 256–267. [Google Scholar] [CrossRef] [Green Version]
  33. Hu, J.; Ji, C.; Chen, S.; Li, S.; Zou, Y.; Zhou, Z.; Wang, B. Two-Position Impact Behavior and Interference Mechanism of CFF/PEEK Thermoplastic Composites. Int. J. Mech. Sci. 2022, 232, 107644. [Google Scholar] [CrossRef]
  34. Reis, P.N.B.; Coelho, S.R.M.; Bezazi, A. Effect of Impact Position on Repaired Composite Laminates Subjected to Multi-Impacts. Materials 2022, 15, 8039. [Google Scholar] [CrossRef]
  35. Nassir, N.A.; Guan, Z.W.; Birch, R.S.; Cantwell, W.J. Damage Initiation in Composite Materials under Off-Centre Impact Loading. Polym. Test. 2018, 69, 456–461. [Google Scholar] [CrossRef]
  36. Liao, B.; Wang, P.; Zheng, J.; Cao, X.; Li, Y.; Ma, Q.; Tao, R.; Fang, D. Effect of Double Impact Positions on the Low Velocity Impact Behaviors and Damage Interference Mechanism for Composite Laminates. Compos. Part A Appl. Sci. Manuf. 2020, 136, 105964. [Google Scholar] [CrossRef]
  37. Huang, L.; Sun, J.; Zhang, D.; Zhao, J. Numerical Investigation on Composite Laminates under Double-Position Low-Velocity Impacts. J. Reinf. Plast. Compos. 2022, 073168442211474. [Google Scholar] [CrossRef]
  38. ASTM D7136/D7136M-15; Standard Test Method for Measuring the Damage Resistance of a Fiber-Reinforced Polymer Matrix Composite to a Drop-Weight. Impact Event, ASTM Int; ASTM International: West Conchhocken, PA, USA, 2015.
  39. Puck, A.; Kopp, J.; Knops, M. Guidelines for the Determination of the Parameters in Puck’s Action Plane Strength Criterion. Compos. Sci. Technol. 2002, 62, 371–378. [Google Scholar] [CrossRef]
  40. Puck, A.; Schürmann, H. Failure Analysis of FRP Laminates by Means of Physically Based Phenomenological Models. In Failure Criteria in Fibre-Reinforced-Polymer Composites; Elsevier Ltd.: Amsterdam, The Netherlands, 2004; pp. 832–876. ISBN 9780080531571. [Google Scholar]
  41. Rezasefat, M.; Torres, D.B.; Gonzalez-Jimenez, A.; Giglio, M.; Manes, A. A Fast Fracture Plane Orientation Search Algorithm for Puck’s 3D IFF Criterion for UD Composites. Mater. Today Commun. 2021, 28, 102700. [Google Scholar] [CrossRef]
  42. Falzon, B.G.; Apruzzese, P. Numerical Analysis of Intralaminar Failure Mechanisms in Composite Structures. Part I: FE Implementation. Compos. Struct. 2011, 93, 1039–1046. [Google Scholar] [CrossRef]
  43. Bažant, Z.P.; Oh, B.H. Crack Band Theory for Fracture of Concrete. Matériaux Constr. 1983, 16, 155–177. [Google Scholar] [CrossRef] [Green Version]
  44. Benzeggagh, M.L.; Kenane, M. Measurement of Mixed-Mode Delamination Fracture Toughness of Unidirectional Glass/Epoxy Composites with Mixed-Mode Bending Apparatus. Compos. Sci. Technol. 1996. [Google Scholar] [CrossRef]
  45. González, E.V.; Maimí, P.; Camanho, P.P.; Lopes, C.S.; Blanco, N. Effects of Ply Clustering in Laminated Composite Plates under Low-Velocity Impact Loading. Compos. Sci. Technol. 2011, 71, 805–817. [Google Scholar] [CrossRef] [Green Version]
  46. 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]
  47. ABAQUS 2016 Documentation; ABAQUS Theory Manual. Dassault Systèmes: Paris, France, 2017.
  48. Lin, S.; Thorsson, S.I.; Waas, A.M. Predicting the Low Velocity Impact Damage of a Quasi-Isotropic Laminate Using EST. Compos. Struct. 2020, 251, 112530. [Google Scholar] [CrossRef]
  49. Li, X.; Ma, D.; Liu, H.; Tan, W.; Gong, X.; Zhang, C.; Li, Y. Assessment of Failure Criteria and Damage Evolution Methods for Composite Laminates under Low-Velocity Impact. Compos. Struct. 2019, 207, 727–739. [Google Scholar] [CrossRef]
  50. Liu, P.F.; Liao, B.B.; Jia, L.Y.; Peng, X.Q. Finite Element Analysis of Dynamic Progressive Failure of Carbon Fiber Composite Laminates under Low Velocity Impact. Compos. Struct. 2016, 149, 408–422. [Google Scholar] [CrossRef]
Figure 1. Impact locations on skin composite laminates.
Figure 1. Impact locations on skin composite laminates.
Materials 16 00914 g001
Figure 2. The FE models: (a) Impact on coupon specimen, (b) Multiple impacts on skin panel.
Figure 2. The FE models: (a) Impact on coupon specimen, (b) Multiple impacts on skin panel.
Materials 16 00914 g002
Figure 3. Comparison of force–time curves for the panel subjected to multiple impacts from simulations with the rest step (RS) and without the rest step (nRS).
Figure 3. Comparison of force–time curves for the panel subjected to multiple impacts from simulations with the rest step (RS) and without the rest step (nRS).
Materials 16 00914 g003
Figure 4. Comparison of experimental and numerical force–displacement curves for a low-velocity impact at 8 J and 12 J on a coupon specimen [29].
Figure 4. Comparison of experimental and numerical force–displacement curves for a low-velocity impact at 8 J and 12 J on a coupon specimen [29].
Materials 16 00914 g004
Figure 5. Layer-by-layer prediction of the matrix damage area, layer 1 being the impact side and layer 24 being the rear side of the specimen.
Figure 5. Layer-by-layer prediction of the matrix damage area, layer 1 being the impact side and layer 24 being the rear side of the specimen.
Materials 16 00914 g005
Figure 6. Experimental and numerical force–time, force–displacement, and energy–time curves for (ac) Single impact at 40 J at the center, (df) Single impact at 35 J at the center, and (gi) Second impact off-center at 25 J (after the first impact at 40 J at the center).
Figure 6. Experimental and numerical force–time, force–displacement, and energy–time curves for (ac) Single impact at 40 J at the center, (df) Single impact at 35 J at the center, and (gi) Second impact off-center at 25 J (after the first impact at 40 J at the center).
Materials 16 00914 g006aMaterials 16 00914 g006b
Figure 7. Comparison of force–time (a), displacement–time (b), and energy–time (c) curves for 40 J impact at the center and off-center locations.
Figure 7. Comparison of force–time (a), displacement–time (b), and energy–time (c) curves for 40 J impact at the center and off-center locations.
Materials 16 00914 g007
Figure 8. Comparison of the center and off-center low-velocity impact responses. (a) Peak impact force, (b) Maximum displacement, (c) Ratio of absorbed energy to impact energy, and (d) Projected delamination area.
Figure 8. Comparison of the center and off-center low-velocity impact responses. (a) Peak impact force, (b) Maximum displacement, (c) Ratio of absorbed energy to impact energy, and (d) Projected delamination area.
Materials 16 00914 g008
Figure 9. Comparison of force–time, displacement–time, and energy–time curves between a single and a double 40 J low-velocity impact at: (ac) S1, (df) S2, (gi) S3, and (jl) S4.
Figure 9. Comparison of force–time, displacement–time, and energy–time curves between a single and a double 40 J low-velocity impact at: (ac) S1, (df) S2, (gi) S3, and (jl) S4.
Materials 16 00914 g009
Figure 10. Comparison impact response for pristine panel and panel with pre-existing damage: (a) Peak impact force, (b) Maximum displacement, (c) Ratio of absorbed energy to impact energy, and (d) Sum of projected delamination area.
Figure 10. Comparison impact response for pristine panel and panel with pre-existing damage: (a) Peak impact force, (b) Maximum displacement, (c) Ratio of absorbed energy to impact energy, and (d) Sum of projected delamination area.
Materials 16 00914 g010
Figure 11. Comparison of the projected damaged area of D-S1-40 J and D-S4-40 J models after the first (at center) and second (off-center) impacts.
Figure 11. Comparison of the projected damaged area of D-S1-40 J and D-S4-40 J models after the first (at center) and second (off-center) impacts.
Materials 16 00914 g011
Table 1. Mechanical properties of CFRP composite.
Table 1. Mechanical properties of CFRP composite.
Laminae PropertiesInterface Properties
E 11 [GPa]157.5 * E n = E s = E t [ GPa / mm ] 5 ***
E 22 = E 33 [GPa]9.9 *N [MPa]33.0 ***
G 12 = G 13 [GPa]4.95 *T [MPa]54.0 ***
G 23 [GPa]3.21 * G n C   [ N / mm ] 0.6 ***
c 1 ,   c 2 ,   c 3 [MPa]10.2 *, −0.5 *, 0.1 * G s C [ N / mm ] 2.1 ***
v 12 = v 13 0.24 *
v 23 0.35 **
R t [MPa]2550.0 *
R c [MPa]1350.0 **
R t [MPa]57.5 *
R c [MPa]199.8 **
R [MPa]97.0 *
G f t ,   G f c [N/mm]133.0 **, 40.0 **
G m t ,   G m c [N/mm]0.6 **, 2.1 **
* Data from test; ** Data from [4,45]; *** Data from [25,46].
Table 2. Parametric study cases for different single low-velocity impact events at different locations, and for different multiple low-velocity impacts at different locations (for the double impacts, the notations only represent the results of the second impact since the first impact is equal to 40 J for all impact cases).
Table 2. Parametric study cases for different single low-velocity impact events at different locations, and for different multiple low-velocity impacts at different locations (for the double impacts, the notations only represent the results of the second impact since the first impact is equal to 40 J for all impact cases).
Single Impacts at Different LocationsDouble Impact at Different Locations
LocationImpact Energy (J)CodeImpact Energy (J)Code
Center25S-C-25 J--
35S-C-35 J--
40S-C-40 J--
S125S-S1-25 J25D-S1-25 J
35S-S1-35 J35D-S1-35 J
40S-S1-40 J40D-S1-40 J
S225S-S2-25 J25D-S2-25 J
35S-S2-35 J35D-S2-35 J
40S-S2-40 J40D-S2-40 J
S325S-S3-25 J25D-S3-25 J
35S-S3-35 J35D-S3-35 J
40S-S3-40 J40D-S3-40 J
S425S-S4-25 J25D-S4-25 J
35S-S4-35 J35D-S4-35 J
40S-S4-40 J40D-S4-40 J
Table 3. Comparison of the experimental and numerical absorbed energies and projected delamination areas [29].
Table 3. Comparison of the experimental and numerical absorbed energies and projected delamination areas [29].
Impact EnergyAbsorbed EnergyDelamination Area
Experimental (J)NumericalExperimental (mm2)Numerical
Value (J)Error (%)Value (mm2)Error (%)
8 J3.412.98−12.6420406−3.3
10 J4.433.82−13.7605567−6.3
12 J4.884.57−6.3750739−1.5
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rezasefat, M.; Beligni, A.; Sbarufatti, C.; Amico, S.C.; Manes, A. Experimental and Numerical Study of the Influence of Pre-Existing Impact Damage on the Low-Velocity Impact Response of CFRP Panels. Materials 2023, 16, 914. https://doi.org/10.3390/ma16030914

AMA Style

Rezasefat M, Beligni A, Sbarufatti C, Amico SC, Manes A. Experimental and Numerical Study of the Influence of Pre-Existing Impact Damage on the Low-Velocity Impact Response of CFRP Panels. Materials. 2023; 16(3):914. https://doi.org/10.3390/ma16030914

Chicago/Turabian Style

Rezasefat, Mohammad, Alessio Beligni, Claudio Sbarufatti, Sandro Campos Amico, and Andrea Manes. 2023. "Experimental and Numerical Study of the Influence of Pre-Existing Impact Damage on the Low-Velocity Impact Response of CFRP Panels" Materials 16, no. 3: 914. https://doi.org/10.3390/ma16030914

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