Next Article in Journal
Carboxymethyl Cellulose/Gelatin Hydrogel Films Loaded with Zinc Oxide Nanoparticles for Sustainable Food Packaging Applications
Next Article in Special Issue
The Injected Foaming Study of Polypropylene/Multiwall Carbon Nanotube Composite with In Situ Fibrillation Reinforcement
Previous Article in Journal
Cellulose Nanocrystals (CNCs) and Its Modified Form from Durian Rind as Dexamethasone Carrier
Previous Article in Special Issue
Experimental and Numerical Investigation on Fatigue Properties of Carbon Fiber Cross-Ply Laminates in Hygrothermal Environments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Parametric Study on Low-Velocity Impact (LVI) Damage and Compression after Impact (CAI) Strength of Composite Laminates

1
AVIC Manufacturing Technology Institute, Beijing 101300, China
2
School of Transportation Science and Engineering, Beihang University, Beijing 100191, China
3
Jingdezhen Research Institute of Beihang University, Jingdezhen 333000, China
*
Author to whom correspondence should be addressed.
Polymers 2022, 14(23), 5200; https://doi.org/10.3390/polym14235200
Submission received: 11 October 2022 / Revised: 12 November 2022 / Accepted: 26 November 2022 / Published: 29 November 2022

Abstract

:
A full-scale model for predicting low-velocity impact (LVI) damage and compression after impact (CAI) strength was established based on a subroutine of the material constitutive relationship and the cohesive elements. The dynamic responses of the laminate under impact load and damage propagation under a compressive load were presented. The influences of impact energy and ply thickness on the impact damage and the CAI strength were predicted. The predicted results were compared with the experimental ones. It is shown that the predicted value of the CAI strength is in good agreement with the experimental result. As the impact energy reaches a certain value, the CAI strength no longer decreases with the increase in the impact energy. Decreasing the ply thickness can effectively improve the damage resistance and CAI strength.

1. Introduction

One of the advantages of composite laminates is that different fiber-reinforced materials and design angles can be selected according to the load requirements of various structures, which provides greater freedom for the structural design and ply scheme optimization. However, the obvious anisotropy and low inter-laminar strength characteristics of the composites can also lead to multiple failure modes of structures, including fiber breakage, resin cracking [1,2] and delamination [3,4,5].
The impact damage and residual strength of composite laminates have attracted lots of attention. A lot of related works have been published, and experimental testing is one of the most direct and effective research methods. Using an experimental method, Sergii [6] investigated the impact and compression after impact (CAI) behavior of thin carbon fiber epoxy plates with two different discontinuous fiber lengths. Compared with the behavior of the continuous fiber composite laminate, the results show that longer fibers provide better damage tolerance. Gliszczynski [7] presented the global failure modes of channel section profiles subjected to CAI during experimental investigations, and the laminates with different angle ply arrangements and different impact angles were discussed. Elamin [8] studied the CAI behavior of carbon-fiber foam-core sandwich composites in low temperature conditions, and the experimental tests were conducted across temperature ranging from −70 to 23 °C. The results showed that low temperature lead to a significant reduction in the CAI strength by 60%. Abdulhamid [9] experimentally studied the CAI values of tapered composite laminate with a thickness from 4 to 6 mm. It was found that the compressive behavior of the specimens is mainly governed by the discontinuity of the neutral axis in the tapered region. Wang [10] analyzed the effect of stacking sequences on delamination morphology and the compression damage propagation mode of thin laminates. The CAI tests were performed with a non-standard device to suppress the global buckling of the plates. It was concluded that the position of the initially propagating delamination is located on one side or both sides of the 0° layer. Liu [11] investigated the CAI behavior of composite panels with foam-filled hat stiffener. The experimental results show that material compression damage is caused by impact damage, and the impact location has an obvious influence on the CAI strength. Zhang [12] compared the response of tubular composite structures with/without a honeycomb core under low-velocity impact (LVI) and CAI, and the results illustrate that the core material has a greater energy absorption and a higher CAI strength.
With the development and improvement of the finite element method, numerical simulation has become more and more popular. Tian [13] proposed an equivalent damage model to predict the CAI strength of composite laminates and stiffened panels, and the numerical model achieved a high computational efficiency by utilizing the soft inclusion method. Shah [14] studied the damage tolerance and CAI values of resin-infused thermoplastic three-dimensional (3D) woven composites which were compared with the conventional resin-infused thermoset composite, and the results show that thermoplastic matrix has a better damage tolerance. Combining the experimental tests and the numerical simulation, Tuo [15] investigated the damage and failure mechanism of thin composite laminates under the CAI. Zhang [16] analyzed the influence of the off-axis angles on the LVI and CAI damage mechanism of three-dimensional woven carbon/epoxy composites, and the results indicate that the off-axis samples exhibited higher damage volumes and out-of-plane deformation. Applying the finite element method, Bull [17] studied the effect of the parameters observed from the previous studies on the CAI. It was found that toughness plays a more significant role than permanent out-of-plane deformation does, and they evidenced and quantified the role of the undamaged cone. Lu [18] numerically investigated the effect of geometries on the CAI behavior of carbon reinforced composites. The results show that width-to-thickness ratio has a significant effect on the failure load, whereas the effect of plate aspect ratio has a smaller effect. To effectively predicting the CAI strength of large-scale composite structures, Johannes [19] developed a model using stacked sub-laminates and the explicit dynamics approach. It was found that the CAI strength in [45/0/90/−45]4s sub-laminate scaled laminates can be predicted accurately, and the maximum error is 10%. Yang [20] analyzed the effect of the impactor diameter on the impact damage and CAI strength, and it was found that the CAI strength increased by approximately 34~53% with the decrease in the impactor diameter.
To improve the impact resistance of the composites, researchers have developed ultra-thin prepregs in recent years, which can delay the initiation and propagation of inter-laminar and intra-laminar damage, such as delamination, matrix cracking and fiber breakage [21,22]. Cameron [23] studied the effect of using thin plies to increase the bearing strength of composite laminates, and it was found that shifting from a 100% conventional ply laminate to a 100% thin ply laminate gives an increase of 47% in the strength at the initiation of the damage. Using different formulations of thin ply composites, Cugnoni [24] evaluated the effect of the fiber and matrix constituents on the initiation of the damage and their strength. The results show that the delamination threshold load of thin laminates was increased by 15% compared to that of the conventional laminates. Saito [25] compared the impact resistance of carbon fiber reinforced polymer laminates with two different ply thicknesses, and the results show that the impact resistance of thin ply laminates is 23% higher than that of the laminates with a conventional ply thickness. However, by investigating the LVI resistance of thin ply in comparison with conventional aluminium-carbon laminates, Drodziel [26] stated a different conclusion, which was that the thin ply carbon fiber reinforced polymer neither increases the LVI resistance nor significantly changes the low impact response.
In order to evaluate the impact damage resistance and residual strength of composite materials, it is necessary to use numerical methods to simulate the progressive damage and predict the bearing capacity of the composite structure that is damaged and to study the influence of material parameters on its impact response. In this study, a Finite Element Model (FEM) for predicting the impact damage and the CAI strength of the composite laminates was established. The whole analysis process of impact damage and CAI strength was conducted using the FEM. The stress state of each element was controlled through a subroutine to simulate the initiation and propagation of in-plane damage and delamination between the plies. The impact damage and the CAI strength of the carbon fiber composite laminates were simulated, and the influence of the impact energy and ply thickness on the impact damage area and CAI strength were studied.

2. Experimental Tests

Standard specimens were prepared by CCF800/AC531 unidirectional reinforced prepreg with 0.14 mm thickness, and the basic mechanical properties are presented in Table 1. Where E11 and E22 are the modulus in the fiber and transverse direction, respectively, G12 is the in-plane shear modulus, v12 is Poisson’s ratio, Xt, Xc, Yt and Yc are the tensile strength in the fiber direction, the compressive strength in the fiber direction, the tensile and compressive strength in the transverse direction, respectively, and S is the in-plane shear strength. In this study, the fiber was considered to be uniformly distributed, and the composite laminate was simplified to a uniform orthotropic material. According to the ASTM standards [27,28], the specimens were cured in the autoclave. The laminates were heated to 180 °C, and then, they were cured under the curing pressure of 0.6 MPa for 150 min. The length of specimens was 150 mm, and the width was 100 mm. Thirty-two plies constitute the layup of the thick section with a stacking sequence of [45/0/−45/90]4s.
Impact tests were performed using a drop weight testing device. The impactor weight was 5.5 kg, and it has a hemispherical shape puncher with a 16 mm diameter. Impact energy was 6.67 J/mm according to the ASTM D7136 standard, and impact energy of the specimens was 30 J in this paper. During the impact process, the specimens were calmed by a fixture with rectangular opening of 125 × 75 mm2. Figure 1a,b shows the impact test device and compression test device, respectively.
Figure 2 illustrates the damage of specimens after impact. From Figure 2, it can be seen that a pit appeared on the upper surface of the laminate, and fiber breakage and matrix cracking damage appeared on the lower surface. The shape of the pit on the upper surface was approximately circular, the matrix cracking was along the fiber direction (45° direction) on the lower surface. The delamination, matrix cracking and fiber breakage damage caused by impact in the thickness direction were obtained using a C-scan technique. The damage area obtained from the C-scan is shown in Figure 2c. The damage shape was approximately elliptical, and the damage ranges were different for the different specimens. The contact forces of 5 specimens during impact are illustrated in Figure 3a, and the maximum contact force varies from 10 to 12 kN, and it can be seen that the numerical predict curve with the 3D criterion agrees well with experimental ones, and the result with 2D criterion is much larger. The damage areas and CAI values of the specimens are listed in Table 2. It can be seen that the impact-induced damage area is negatively correlated with the maximum contact force. The relationship between the damage area and the contact force is listed in Figure 3b, and the correlation coefficient is about −0.88. It should be pointed out that there is dispersion among the different samples, but the average properties of the material are basically stable, and we verified the numerical model by comparing them with the experimental average values.
The impact tests were followed by the CAI tests, and the CAI device was set up according to the ASTM D7137 standard. Supports in the length and width directions of the specimens were provided on both sides to prevent bend failure due to the bending deformation, and the compressive displacement rate was 1.0 mm/min. Figure 4 shows the comparison of the typical specimens after compression failure between the experimental test and the numerical simulation. It clear that the failure modes between the experimental test and the numerical simulation are in good agreement. The failure mode is a lateral failure, and the failure area runs through the width, and the failure location appears to be in the middle of the specimen’s length. This failure mode is one of the common acceptable failure modes of the CAI, as mentioned in the ASTM D7137 standard. Figure 5a illustrates the compressive displacement versus the load curves, and the six curves are well correlated. Table 2 lists the failure compressive displacements and the CAI values of the specimens. It can be concluded that the impact damage area is negatively correlated with the CAI strength, and the relationship between the damage area and the contact force is listed in Figure 3b; the correlation coefficient is about −0.96. Additionally, the maximum difference of the CAI strengths is about 50 MPa for the five specimens.

3. Numerical Model

3.1. Failure Criterion

3.1.1. Intra-Laminar Damage Model

There are a lot of investigations on the failure criteria of composite laminates. Many researchers have proposed different failure criteria, which can be divided into two types: a single mode failure criterion and a multi-mode failure criterion. The multi-mode criterion considers the contribution of different types of loads to structural failure. In this study, the 3D Hashin criterion [29] was used, and the composite damages in this criterion were classified into four types: a fiber tensile failure, a fiber compression failure, a transverse tensile failure and a transverse compression failure. The equations are:
( σ 11 X c ) 2 + ( τ 12 S 12 ) 2 + ( τ 13 S 13 ) 2 1         σ 11 0
( σ 11 X c ) 2 1 σ 11 < 0
( σ 22   +   σ 33 Y t ) 2 + ( τ 23 2     σ 22 σ 33 ) S 23 2 + ( τ 12 S 12 ) 2 + ( τ 13 S 13 ) 2 1 ( σ 22 + σ 33 ) 0
( σ 22   +   σ 33 2 S 12 ) 2 + σ 22   +   σ 33 Y c ( ( Y c 2 S 12 ) 2 1 ) + ( τ 23 2     σ 22 σ 33 ) S 23 2 + ( τ 12 S 12 ) 2 + ( τ 13 S 13 ) 2 1 ( σ 22 + σ 33 ) < 0
where σ i i (i = 1,2,3) are the normal stresses in the i direction, τ i j (I,j = 1,2,3) are the shear stresses in the ij plane, Xt, Xc, Yt and Yc are the tensile strength in the fiber direction, the compressive strength in the fiber direction and the tensile and compressive strength in the transverse direction, respectively, and Sij is the in-plane shear strength. The mechanical properties of the CCF800/AC531 composite ply are listed in Table 1.
When a certain failure mode appears, the corresponding material properties should be degraded, and there will be no other failure modes. The material properties after degradation can be expressed by state variables, and the Camanho model [30] was used to characterize the material state when an in-plane failure occurred.

3.1.2. Inter-Laminar Damage Model

The interface between two adjacent plies was simulated by the cohesive elements. When the interlaminar stresses exceed the strength, the failure process of the composite material is characterized by the continuous reduction of the material stiffness. The exponential constitutive model and the bilinear constitutive model are the two most popular damage constitutive models for interface damage [31,32]. The bilinear constitutive model was used in this paper, and the linear elastic stress versus strain relationship of the interface without delamination can be written as:
[ σ n τ s τ t ] = [ K n 0 0 0 K s 0 0 0 K t ] [ δ n δ s δ t ]
where σ n , τ s and τ t are the interlaminar normal stress and the interlaminar shear stresses, respectively, K n , K s and K t are the corresponding stiffness, and δ n , δ s and δ t are the displacements in three directions. The delamination initiation and propagation of the interface are determined by the second nominal stress criterion and BK criterion [33], as follows:
( σ n T ) 2 + ( τ s S ) 2 + ( τ t S ) 2 = 1
G C = G n C + ( G s C G n C ) ( G s C + G t C G s C + G n C ) η
where T and S are the normal tensile strength and shear strength, respectively, G C is the energy release rate for mixed delamination propagation, G n C , G s C and G t C are the energy release rate in the normal and shear directions, and η is the parameter related to interface property, which was selected to be 1.0 in this paper. The property parameters of interface are listed in Table 3.

3.2. Finite Element Model

Composite laminates are made by bonding multiple fiber-reinforced composite plies. The lay-up of the laminate was [45/0/−45/90]4s, and the ply thickness was 0.14 mm. According to the size of composite laminate specimen specified in the test standard, an FEM of the composite laminates was established. Since the stiffness of the impactor is much greater than that of the composite laminate, the deformation of the impactor can be ignored, and the impactor was set as the rigid body in the FEM. The FEM of the composite laminate is shown in Figure 6. The composite plies were modeled by solid elements, and the interfaces were modeled by cohesive elements. Since the angles of the adjacent plies of laminates are not the same in this paper, a total of 31 interface plies were utilized. The interface plies were numbered from 1 to 31 from the impact surface to the back surface.

4. Impact Damage and Compressive Strength

4.1. Impact-Induced Damage

An impact damage prediction model of the composite laminate was established to calculate the response under the impact load. According to the ASTM D 7136 standard, the kinetic energy of the impactor was implemented using the initial velocity, and the impact energy was set at 30 J. Figure 7 shows intra-laminar fiber failure and the matrix cracking states at different impact times. It can be seen that matrix cracking appeared earlier than it did in the fiber due to its low strength, and the final damage zone of the fiber and matrix caused by the impact is roughly the same, which is approximately circular.
The damage states of typical interfaces are shown in Figure 8. It is shown that the delamination shapes of these interfaces are relatively close, which are similar to the number “8”. The propagation directions of the delamination are not the same, and the delamination tends to propagate along the direction that is perpendicular to the fiber. Figure 3 shows the time history of contact force between the impactor and the laminate during the impact process. The calculated values are in good agreement with the experimental ones, which validates the established FEM.

4.2. Simulation of Compression Failure

The laminate state at the end of impact was used as the initial condition of the compression simulation. The failure state of the fiber, matrix and delamination at the beginning of compression was defined by field variables. Considering the storage space and computing efficiency of the computer, it was impossible to simulate the compression loading rate under real experimental conditions. The loading rate needs to be increased to ensure that the calculation time is within an acceptable range, and the enforced displacement loading rate used in this research was 50 mm/s.
Figure 9 illustrate comparison between the fiber failure propagation and matrix failure propagation of typical plies. It can be seen that the failure range of fiber and matrix has a certain propagation under the compression load. The matrix damage area is much larger than the fiber damage area is due to the lower matrix strength compared with the fiber strength.
Figure 5 compares the reaction forces predicted by the numerical model and the experiments, and the predicted curves are in agreement with experimental ones. From Figure 5, it shows that the reaction force is linearly proportional to the compression displacement in the initial stage. The reaction force begins to gradually decrease when the compression displacement reaches about 0.81 mm, which indicates that the damage began to propagate in a large area. When the displacement reaches 0.85 mm, the reaction force drops sharply, and the whole structure fails. The result indicates that when the compression displacement is less than 0.81 mm, the damage does not propagate significantly. While the compression displacement increases from 0.81 mm to 0.85 mm, the damage rapidly propagates and the structure fails, as shown in Figure 10. Table 2 shows the comparison between the CAI prediction values and the experimental average values. It illustrates that the predicted results are in good agreement with the experimental ones. The error between the predicted value using the 3D criterion and the experiment is about 4.7% for the CAI strength, which validates the presented model.

4.3. Influences of Damage Criterion

The selection of the damage criterion affects the prediction of the in-plane properties of the composite plies, and then, it affects the prediction of structural failure process and final damage state. The influences of the damage criteria on the prediction results of the impact damage and the CAI strength were studied through the FEM that was established in this paper. The 3D Hashin criterion is shown in Section 3.1.1, and the 2D Hashin criterion can be obtained by setting the parameters related to three directions (i.e., the thickness direction) in Equations (1)–(4) to 0. The damage characteristics and CAI strengths using the 3D and 2D Hashin criteria were calculated and compared.
Figure 11 shows the comparison of the matrix cracking obtained by using two in-plane failure criteria with an impact energy of 30 J. It can be seen that the matrix cracking area predicted using the 3D criterion is larger than that for the 2D criterion, and the contact force versus time curves and the reaction force versus displacement curves are shown in Figure 12. The contact time between the impactor and the laminate is shorter when the 2D criterion was used, while the maximum contact force increases by about 50%. The reason for this phenomenon is that when 2D criterion was used, the energy absorbed by the laminate during the impact process is less than it was while we were using the 3D criterion, as shown in Figure 13. The reduction of the energy absorption makes the damage area smaller, thus, the remaining stiffness of the laminate is relatively larger, resulting in a short duration of the impact process and a larger impact contact force. The impact damage area obtained using the 2D criterion is relatively small, and the CAI strength is relatively larger, which illustrates that the 3D criterion is better than the 2D criterion for predicting the CAI behavior. The damage caused by out-of-plane stress should be considered in the FEM. In the following calculation, the 3D Hashin criterion was employed to simulate the composite damage of each ply.

5. Parameter Influences

5.1. Influence of Impact Energy on Compression after Impact

The composite structures may be subjected to different impact energies, which will cause the bearing capacity of the composite structures to decline by different degrees. To study the influence of impact energy on CAI strength of the composite materials, the impact damage and CAI strengths under impact energies of 10, 20, 40 and 50 J were compared with the result of the test using 30 J. The impact damage of the composite laminate was calculated, and Figure 14 shows the comparison of delamination damage of typical interfaces caused by different impact energy conditions.
According to Figure 8 and Figure 14, it can be seen that the delamination area gradually increases when the impact energy increases from 10 to 30 J. Compared with impact energy of 30 J, increasing the impact energy to 40 and 50 J has little effect on the delamination damage area. With the increase in impact energy, the delamination damage area will not continue to increase after the impact energy exceeds 30 J. Figure 15a shows the variation of the contact force versus contact time curves under different impact energy conditions, and it can be seen that the contact force increases with the increase in the impact energy. From Figure 15b, it is clear that the reaction force gradually decreases when the impact energy increases from 10 to 30 J. After the impact energy exceeds 30 J, the maximum reaction force remains basically unchanged, illustrating that CAI strength does not decrease continually after the impact energy exceeds 30 J, as shown in Figure 15c.
To further explain the above phenomenon, Figure 16 shows the changes to the impactor energy and the energy absorption of the composite laminate under different impact energy conditions. It can be seen that the contact time between the impactor and the composite laminate becomes longer with the increase in the impact energy, and the energy absorption of the composite laminate also gradually increases. The in-plane damage of the typical laminates under three impact energy conditions, i.e., 30, 40 and 50 J, is shown in Figure 17, and the in-plane damage area gradually increases with the increase in the impact energy. The location of the fiber damage and matrix damage is the same as that in the delamination test, which is near the impact point. The results show that when the impact energy exceeds 30 J, the in-plane damage area of some of the plies gradually increases, which means that more impactor energy is absorbed. However, the in-plane damage and the inter-laminar delamination projection area of all of the plies do not increase significantly, indicating that the effect of impact energy on the CAI strength is not obvious after it exceeds 30 J.

5.2. Influence of Ply Thickness on Compression after Impact

Raman [34] shows that reducing the ply thickness may reduce the damage range of the composite structure. To study the influence of the ply thickness on compression after the impact, in this section, the impact response and CAI strength of laminates with different ply thicknesses were calculated. The results were compared with the one using conventional ply thickness of 0.14 mm, and the cured ply thicknesses and corresponding layups are shown in Table 4.
Figure 18 illustrates the compressive load versus displacement curves and the CAI strengths with different ply thicknesses. The reaction force versus displacement curves with different single ply thicknesses are shown in Figure 18a. It can be seen that the final reaction force increases with the decrease in single ply thickness of the composite laminate. According to Figure 18b, the CAI strength increases gradually with the decrease in the single ply thickness. Compared with the result of the 0.14 mm single ply thickness, the CAI strength increases by 75 MPa when the single ply thickness decreases to 0.04 mm, namely, the CAI strength increases by 22.4%. The results indicate that decreasing the single ply thickness can effectively improve the damage resistance and CAI strength of the composite laminates. This is consistent with Wagih [35], Huang [36] and Saito [25] who studied the effect of ply thickness on damage resistance and CAI through experiments, and they found that the CAI strength of thinner ply was 23% higher than that of a traditional ply thickness.
To further analyze the mechanism of the effect of the single ply thickness on the CAI strengths, Figure 19 shows the comparison of the delamination, fiber breakage and matrix cracking areas caused by the impact. It can be seen that the single ply thickness has no obvious effect on the fiber breakage area, but the matrix cracking and delamination areas decrease with the decrease in the single ply thickness. Compared with the other conditions, the matrix cracking area with ply thicknesses of 0.06 and 0.04 mm decreases, obviously. It shows that decreasing the single ply thickness can restrain the propagation of matrix cracking caused by impact, and thus, the residual strength of the laminate in compression stage (i.e., CAI strength) is improved.
The comparison of contact forces and energy absorption with different single ply thicknesses are shown in Figure 20. From Figure 20a, it can be seen that the contact force increases with the decrease in the single ply thickness, indicating that the laminate stiffness becomes larger, that is, the damage area of laminates with a thinner ply is relatively small. Figure 20b shows that the energy absorbed by the laminate decreases gradually with the decrease in the single ply thickness, which is one of the reasons for the smaller damage propagation range during the impact process.

6. Conclusions

An FEM for predicting the impact damage and CAI strength of composite laminates was established. The whole analysis process of the impact damage and CAI strength was conducted using the FEM. The predicted CAI strength of the composite laminate was compared with the experimental one. The influences of impact energy, damage criterion and single ply thickness on the impact damage area and CAI strength were investigated. Four important results emerging from the research are as follows:
  • The damage zone of the fiber breakage and matrix cracking caused by the impact is approximately circular.
  • The predicted CAI value is in good agreement with the experimental one, and the error is 4.7%.
  • The CAI strength decreases as impact energy increases from 10 to 30 J. After the impact energy exceeds 30 J, the CAI strength remains basically unchanged.
  • Decreasing the single ply thickness can effectively improve the damage resistance and CAI strength of the composite laminates. When the single ply thickness decreases to 0.04 mm, the CAI strength increases by 22.4%.

Author Contributions

Conceptualization, S.G. and X.L.; methodology, S.G. and T.L.; validation, G.B.; writing—original draft preparation, S.G. and X.L.; writing—review and editing, J.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Defense Basic Research Program of China (Grant No. JCKY2019205C002) and the National Natural Science Foundation of China (Grant No. 51875026 and Grant No. 52275231).

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. González, E.; Maimí, P.; Camanho, P.; Lopes, C.; 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]
  2. Zobeiry, N.; Vaziri, R.; Poursartip, A. Characterization of strain-softening behavior and failure mechanisms of composites under tension and compression. Compos. Part A Appl. Sci. Manuf. 2015, 68, 29–41. [Google Scholar] [CrossRef]
  3. Abrate, S. Impact on Composite Structures; Cambridge University Press: Cambridge, UK, 2005. [Google Scholar]
  4. Gong, W.; Chen, J.; Patterson, E.A. An experimental study of the behaviour of delaminations in composite panels subjected to bending. Compos. Struct. 2015, 123, 9–18. [Google Scholar] [CrossRef]
  5. Cantwell, W.; Morton, J. The impact resistance of composite materials—A review. Composites 1991, 22, 347–362. [Google Scholar] [CrossRef]
  6. Sergii, G.K.; Chris, V.; Oleksandr, G.K. An experimental investigation on low-velocity impact response and compression after impact of a stochastic, discontinuous prepreg tape composite. Compos. Part A 2021, 149, 106524. [Google Scholar]
  7. Gliszczynski, A.; Bogenfeld, R.; Degenhardt, R.; Kubiak, T. Corner impact and Compression After Impact (CAI) of thin-walled compo-site profile-an experimental study. Compos. Struct. 2020, 248, 112502. [Google Scholar] [CrossRef]
  8. Elamin, M.; Li, B.; Tan, K.T. Compression after impact performance of carbon-fiber foam-core sandwich composites in low tem-perature arctic conditions. Compos. Struct. 2021, 261, 113568. [Google Scholar] [CrossRef]
  9. Abdulhamid, H.; Bouvet, C.; Michel, L.; Aboissière, J.; Minot, C. Experimental study of compression after impact of asymmetrically tapered compo-site laminate. Compos. Struct. 2016, 149, 292–303. [Google Scholar] [CrossRef] [Green Version]
  10. Wang, X.; Huang, J.; Tan, R.; Su, Y.; Guan, Z.; Guo, X. Experimental investigation on damage mechanisms and buckling behaviors of thin composite laminates in compression after impact. Compos. Struct. 2020, 256, 113122. [Google Scholar] [CrossRef]
  11. Liu, D.; Bai, R.; Lei, Z.; Guo, J.; Zou, J.; Wu, W.; Yan, C. Experimental and numerical study on compression-after-impact behavior of composite panels with foam-filled hat-stiffener. Ocean Eng. 2020, 198, 106991. [Google Scholar] [CrossRef]
  12. Zhang, C.; Tan, K. Low-velocity impact response and compression after impact behavior of tubular composite sandwich structures. Compos. Part B: Eng. 2020, 193, 108026. [Google Scholar] [CrossRef]
  13. Ouyang, T.; Bao, R.; Sun, W.; Guan, Z.; Tan, R. A fast and efficient numerical prediction of compression after impact (CAI) strength of composite laminates and structures. Thin-Walled Struct. 2020, 148, 106588. [Google Scholar] [CrossRef]
  14. Shah, S.; Megat-Yusoff, P.; Karuppanan, S.; Choudhry, R.; Din, I.U.; Othman, A.; Sharp, K.; Gerard, P. Compression and buckling after impact response of resininfused thermoplastic and thermoset 3D woven composites. Compos. Part B Eng. 2021, 207, 108592. [Google Scholar] [CrossRef]
  15. Tuo, H.; Lu, Z.; Ma, X.; Xing, J.; Zhang, C. Damage and failure mechanism of thin composite laminates under low-velocity impact and com-pression-after-impact loading conditions. Compos. Part B 2019, 163, 642–654. [Google Scholar] [CrossRef]
  16. Zhang, D.; Gu, Y.; Zhang, Z.; Jia, M.; Yue, S.; Li, G. Effect of off-axis angle on low-velocity impact and compression after impact damage mechanisms of 3D woven composites. Mater. Des. 2020, 192, 108672. [Google Scholar] [CrossRef]
  17. Bull, D.J.; Spearing, S.M.; Sinclair, I. Image-enhanced modelling of residual compressive after impact strength in laminated composites. Compos. Struct. 2018, 192, 20–27. [Google Scholar] [CrossRef] [Green Version]
  18. Lu, T.; Chen, X.; Wang, H. Predicting compression-after-impact behavior of thermoplastic composite laminates by an experiment-based approach. Compos. Sci. Technol. 2021, 213, 108952. [Google Scholar] [CrossRef]
  19. Johannes, R.; Navid, Z.; Reza, V. Efficient finite element simulation of compression after impact behaviour in quasi-isotropic composite laminates. Compos. Commun. 2021, 28, 100967. [Google Scholar]
  20. Yang, B.; Chen, Y.; Lee, J.; Fu, K.; Li, Y. In-plane compression response of woven CFRP composite after low-velocity impact: Modelling and experiment. Thin-Walled Struct. 2020, 158, 107186. [Google Scholar] [CrossRef]
  21. Mencattelli, L.; Pinho, S. Realising bio-inspired impact damage-tolerant thin-ply CFRP Bouligand structures via promoting dif-fused sub-critical helicoidal damage. Compos. Sci. Technol. 2019, 182, 107684. [Google Scholar] [CrossRef]
  22. Sebaey, T.; Mahdi, E. Using thin-plies to improve the damage resistance and tolerance of aeronautical CFRP composites. Compos. Part A Appl. Sci. Manuf. 2016, 86, 31–38. [Google Scholar] [CrossRef]
  23. Cameron, C.J.; Larsson, J.; Loukil, M.S.; Murtagh, T.; Wennhage, P. Bearing strength performance of mixed thin/thick-ply, quasi-isotropic composite laminates. Compos. Struct. 2020, 261, 113312. [Google Scholar] [CrossRef]
  24. Cugnoni, J.; Amacher, R.; Kohler, S.; Brunner, J.; Kramer, E.; Dransfeld, C.; Smith, W.; Scobbie, K.; Sorensen, L.; Botsis, J. Towards aerospace grade thin-ply composites: Effect of ply thickness, fibre, matrix and interlayer toughening on strength and damage tolerance. Compos. Sci. Technol. 2018, 168, 467–477. [Google Scholar] [CrossRef]
  25. Saito, H.; Morita, M.; Kawabe, K.; Kanesaki, M.; Takeuchi, H.; Tanaka, M.; Kimpara, I. Effect of ply-thickness on impact damage morphology in CFRP laminates. J. Reinf. Plast. Compos. 2011, 30, 1097–1106. [Google Scholar] [CrossRef]
  26. Drodziel, M.; Jakubczak, P.; Bienia, J. Low-velocity impact resistance of thin-ply in comparison with conventional alumini-um-carbon laminates. Compos. Struct. 2021, 256, 113083. [Google Scholar] [CrossRef]
  27. D 7136/D 7136M-07[S] ASTM Standards; Standard Test Method for Measuring the Damage Resistance of a Fiber-Reinforced Polymer Matrix Composite to a Drop-Weight Impact Event. ASTM International: West Conshohocken, PA, USA, 2007.
  28. D 7137/D 7137M-07[S] ASTM Standards; Standard Test Method for Compressive Residual Strength Properties of Damaged Polymer Matrix Compo-site Plates. ASTM International: West Conshohocken, PA, USA, 2007.
  29. Hashin, Z. Failure Criteria for Unidirectional Fiber Composites. J. Appl. Mech. 1980, 47, 329–334. [Google Scholar] [CrossRef]
  30. Camanho, P.P.; Matthews, F.L. A Progressive Damage Model for Mechanically Fastened Joints in Composite Laminates. J. Compos. Mater. 1999, 33, 2248–2280. [Google Scholar] [CrossRef]
  31. Zhang, J.; Zhang, X. Simulating low-velocity impact induced delamination in composites by a quasi-static load model with surface-based cohesive contact. Compos. Struct. 2015, 125, 51–57. [Google Scholar] [CrossRef]
  32. Guo, S.X.; Li, X.Q. The influences of cohesive element compressive stress and its thickness on delamination prediction of compo-site. Fiber Reinf. Plast. Compos. 2019, 2, 20–25. (In Chinese) [Google Scholar]
  33. Borg, R.; Nilsson, L.; Simonsson, K. Simulation of delamination in fiber composites with a discrete cohesive failure model. Compos. Sci. Technol. 2001, 61, 667–677. [Google Scholar] [CrossRef]
  34. Raman, V.; Drissi-Habti, M.; Guillaumat, L.; Khadhour, A. Numerical simulation analysis as a tool to identify areas of weakness in a turbine wind-blade and solutions for their reinforcement. Compos. Part B Eng. 2016, 103, 23–29. [Google Scholar] [CrossRef] [Green Version]
  35. Wagih, A.; Maimí, P.; González, E.; Blanco, N.; de Aja, J.S.; de la Escalera, F.; Olsson, R.; Alvarez, E. Damage sequence in thin-ply composite laminates under out-of-plane loading. Compos. Part A Appl. Sci. Manuf. 2016, 87, 66–77. [Google Scholar] [CrossRef]
  36. Huang, C.; He, M.; He, Y.; Xiao, J.; Zhang, J.; Ju, S.; Jiang, D. Exploration relation between interlaminar shear properties of thin-ply laminates under short-beam bending and meso-structures. J. Compos. Mater. 2017, 52, 2375–2386. [Google Scholar] [CrossRef]
Figure 1. Devices for the CAI test. (a) Impact test device. (b) Compression test device.
Figure 1. Devices for the CAI test. (a) Impact test device. (b) Compression test device.
Polymers 14 05200 g001
Figure 2. Damage of specimens after impact. (a) Upper surface. (b) Lower surface. (c) Image of C-scan.
Figure 2. Damage of specimens after impact. (a) Upper surface. (b) Lower surface. (c) Image of C-scan.
Polymers 14 05200 g002
Figure 3. Contact forces and damage areas of specimens. (a) Comparison of the contact force during impact. (b) Relationship between damage area and contact force.
Figure 3. Contact forces and damage areas of specimens. (a) Comparison of the contact force during impact. (b) Relationship between damage area and contact force.
Polymers 14 05200 g003
Figure 4. Images of specimens after compression failure. (a) Experimental compression failure. (b) Predicted compression failure.
Figure 4. Images of specimens after compression failure. (a) Experimental compression failure. (b) Predicted compression failure.
Polymers 14 05200 g004
Figure 5. Reaction forces and CAI strengths of specimens. (a) Reaction forces versus displacement curves. (b) Relationship between damage area and CAI strength.
Figure 5. Reaction forces and CAI strengths of specimens. (a) Reaction forces versus displacement curves. (b) Relationship between damage area and CAI strength.
Polymers 14 05200 g005
Figure 6. FEM of the composite laminate.
Figure 6. FEM of the composite laminate.
Polymers 14 05200 g006
Figure 7. Intra-laminar failure of the laminate. (a) Fiber failure zone (red zone). (b) Matrix failure zone (red zone).
Figure 7. Intra-laminar failure of the laminate. (a) Fiber failure zone (red zone). (b) Matrix failure zone (red zone).
Polymers 14 05200 g007aPolymers 14 05200 g007b
Figure 8. Delamination of typical interfaces. (a) The 15th interface. (b) The 18th interface. (c) The 20th interface. (d) The 26th interface.
Figure 8. Delamination of typical interfaces. (a) The 15th interface. (b) The 18th interface. (c) The 20th interface. (d) The 26th interface.
Polymers 14 05200 g008
Figure 9. Intra-laminar failure propagation. (a) Intra-laminar failure of the 7th ply. (b) Intra-laminar failure of the 18th ply.
Figure 9. Intra-laminar failure propagation. (a) Intra-laminar failure of the 7th ply. (b) Intra-laminar failure of the 18th ply.
Polymers 14 05200 g009
Figure 10. Delamination propagation with different compressive displacements.
Figure 10. Delamination propagation with different compressive displacements.
Polymers 14 05200 g010
Figure 11. Comparison of the matrix damage. (a) Damage with 3D Hashin criterion. (b) Damage with 2D Hashin criterion.
Figure 11. Comparison of the matrix damage. (a) Damage with 3D Hashin criterion. (b) Damage with 2D Hashin criterion.
Polymers 14 05200 g011
Figure 12. Results with two damage criteria. (a) The contact force. (b) The reaction force.
Figure 12. Results with two damage criteria. (a) The contact force. (b) The reaction force.
Polymers 14 05200 g012
Figure 13. Kinetic energy variation of the impactor.
Figure 13. Kinetic energy variation of the impactor.
Polymers 14 05200 g013
Figure 14. Delamination under different impact energy conditions. (a) Impact energy 10 J. (b) Impact energy 20 J. (c) Impact energy 40 J. (d) Impact energy 50 J.
Figure 14. Delamination under different impact energy conditions. (a) Impact energy 10 J. (b) Impact energy 20 J. (c) Impact energy 40 J. (d) Impact energy 50 J.
Polymers 14 05200 g014
Figure 15. The impact and compression response under different impact energy conditions. (a) The contact force curves. (b) The compression load versus displacement curves. (c) CAI strength with different impact energy.
Figure 15. The impact and compression response under different impact energy conditions. (a) The contact force curves. (b) The compression load versus displacement curves. (c) CAI strength with different impact energy.
Polymers 14 05200 g015
Figure 16. Variation of impactor energy and energy absorption. (a) Impactor energy curve. (b) Energy absorbed by composite laminate.
Figure 16. Variation of impactor energy and energy absorption. (a) Impactor energy curve. (b) Energy absorbed by composite laminate.
Polymers 14 05200 g016
Figure 17. The intra-laminar failure under different impact energy conditions. (a) Impact energy 30 J. (b) Impact energy 40 J. (c) Impact energy 50 J.
Figure 17. The intra-laminar failure under different impact energy conditions. (a) Impact energy 30 J. (b) Impact energy 40 J. (c) Impact energy 50 J.
Polymers 14 05200 g017
Figure 18. Compressive responses with different ply thicknesses. (a) Reaction force curve. (b) CAI strength with different ply thicknesses.
Figure 18. Compressive responses with different ply thicknesses. (a) Reaction force curve. (b) CAI strength with different ply thicknesses.
Polymers 14 05200 g018
Figure 19. Comparison of impact damage with different single ply thicknesses.
Figure 19. Comparison of impact damage with different single ply thicknesses.
Polymers 14 05200 g019
Figure 20. Comparison of contact forces and energy absorption with different single ply thicknesses. (a) Contact force curves. (b) Energy absorption.
Figure 20. Comparison of contact forces and energy absorption with different single ply thicknesses. (a) Contact force curves. (b) Energy absorption.
Polymers 14 05200 g020
Table 1. Properties of CCF800/AC531 composite ply.
Table 1. Properties of CCF800/AC531 composite ply.
ParameterValue
E11/GPa171.0
E22/GPa8.76
G12/GPa4.2
v120.25
Xt/MPa2889.0
Xc/MPa1562.0
Yt/MPa84.1
Yc/MPa221.0
S/MPa163.0
Table 2. Damage areas and CAI specimens of the specimens.
Table 2. Damage areas and CAI specimens of the specimens.
SpecimenImpact Damage Area (mm2)Maximum Contact Force (kN)CAI Strength (MPa)Failure Displacement (mm)
143812.03240.86
245410.93150.83
337812.43450.89
443611.33160.83
547310.22950.77
Average43611.43200.84
Numerical41911.83340.81
Error (%)3.92.64.73.6
Table 3. Properties of the interface.
Table 3. Properties of the interface.
ParameterValue
GⅠC/(N/m)329.0
GⅡC/(N/m)2215.0
Tensile strength/MPa60.0
Shear strength/MPa112.0
Table 4. Layups for different nominal cured ply thicknesses.
Table 4. Layups for different nominal cured ply thicknesses.
ItemCured Ply Thickness/mmLayup
10.12[45/0/−45/90]5s
20.10[45/0/−45/90]6s
30.08[45/0/−45/90]7s
40.06[45/0/−45/90]9s
50.04[45/0/−45/90]13s
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Guo, S.; Li, X.; Liu, T.; Bu, G.; Bai, J. Parametric Study on Low-Velocity Impact (LVI) Damage and Compression after Impact (CAI) Strength of Composite Laminates. Polymers 2022, 14, 5200. https://doi.org/10.3390/polym14235200

AMA Style

Guo S, Li X, Liu T, Bu G, Bai J. Parametric Study on Low-Velocity Impact (LVI) Damage and Compression after Impact (CAI) Strength of Composite Laminates. Polymers. 2022; 14(23):5200. https://doi.org/10.3390/polym14235200

Chicago/Turabian Style

Guo, Shuangxi, Xueqin Li, Tianwei Liu, Guangyu Bu, and Jiangbo Bai. 2022. "Parametric Study on Low-Velocity Impact (LVI) Damage and Compression after Impact (CAI) Strength of Composite Laminates" Polymers 14, no. 23: 5200. https://doi.org/10.3390/polym14235200

APA Style

Guo, S., Li, X., Liu, T., Bu, G., & Bai, J. (2022). Parametric Study on Low-Velocity Impact (LVI) Damage and Compression after Impact (CAI) Strength of Composite Laminates. Polymers, 14(23), 5200. https://doi.org/10.3390/polym14235200

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