Next Article in Journal
Development of Hybrid Materials Based on Chitosan, Poly(Ethylene Glycol) and Laponite® RD: Effect of Clay Concentration
Next Article in Special Issue
Does a Compatibilizer Enhance the Properties of Carbon Fiber-Reinforced Composites?
Previous Article in Journal
Hydrogel/β-FeOOH-Coated Poly(vinylidene fluoride) Membranes with Superhydrophilicity/Underwater Superoleophobicity Facilely Fabricated via an Aqueous Approach for Multifunctional Applications
Previous Article in Special Issue
Warmforming Flow Pressing Characteristics of Continuous Fibre Reinforced Thermoplastic Composites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on Low-Velocity Impact Response of Novel Short-Fiber-Reinforced Composite Laminates

1
School of Mechanical and Power Engineering, Nanjing Tech University, Nanjing 211800, China
2
Jiangsu Bi-gold New Material Stock Co., Ltd., Zhenjiang 212499, China
3
College of Aerospace Science and Engineering, National University of Defense Technology, Changsha 410073, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Polymers 2023, 15(4), 840; https://doi.org/10.3390/polym15040840
Submission received: 2 January 2023 / Revised: 31 January 2023 / Accepted: 3 February 2023 / Published: 8 February 2023
(This article belongs to the Special Issue Fiber-Reinforced Polymer Composites)

Abstract

:
Short-fiber-reinforced polymers (SFRPs) based on unidirectionally arrayed chopped strands (UACSs) have excellent formability and outstanding mechanical response. The low-velocity impact response, such as the delamination, damage tolerance and energy absorption of UACS composites, are essential to guarantee the stability and safety of composite components in service. The current study investigates the low-velocity impact response of continuous carbon-fiber-reinforced polymer (CFRP) and UACS laminates with vertical slits under drop-weight impact with various impact energies (4, 7 and 11 J). The in-plane size of the studied samples is 100 mm × 100 mm, and the stacking sequence is [0/90]4s. The time–history curves of load and energy are examined during low-velocity impact experiments, as well as the nonvisible damages are obtained by ultrasound C-scan imaging technique. A user-defined subroutine VUMAT, including the Johnson–Cook material and failure model, which is used to simulate the elastic–plastic property of the slits filled with resin, is coded in ABAQUS/Explicit. According to C-scan inspections of the impact-damaged laminates, UACS specimens show more severe delamination as impact energy increases. The damaged area of continuous CFRP laminates under impact energy of 11 J is 311 mm2, while that of UACS laminates is 1230 mm2. The slits have a negative effect on the load-bearing capacity but increase the energy absorption of UACS laminates by approximately 80% compared to the continuous CFRP laminates at 7 J. According to the variables of different damage modes in numerical simulation, cracks appear at the slits and then expand along the direction perpendicular to the slits, leading to the fracture of fiber. Nevertheless, as the damage expands to the slits, the delamination confines the damage propagation. The existence of slits could guide the path of damage propagation.

Graphical Abstract

1. Introduction

Due to their high specific stiffness and strength, fiber composites have been used in many fields over recent years [1]. Among them, carbon-fiber-reinforced polymer composite (CFRP) materials, as an advanced composite material, have sets of benefits such as their low thermal expansion coefficient, lightweight, high specific strength and high modulus. They are extensively used in the aerospace industry, wind power industry, automobile industry, structural reinforcement engineering and other fields [2]. Traditional continuous CFRP composites have certain limitations when they are used to form components with complex shapes. In the curing process, uneven fiber and resin distribution easily occur. The existence of this kind of defect will inevitably produce a complex stress distribution inside the material, which will affect the structure’s safety performance. As a result, SFRP with exceptional fiber fluidity has been the subject of extensive research since the 1960s. For example, sheet molding compound (SMC) has been developed by leaps and bounds due to its excellent fiber fluidity and lower cost [3,4]. The molding processes can achieve proper formation; however, they are unable to regulate the microstructure of short fiber distribution and result in poor mechanical properties [5]. Taketa et al. [6] introduced UACS in which continuous CFRP prepreg is chopped into regular discontinuous distribution to further improve the equilibrium relationship between mechanical property and fiber flowability. Similarly, UACSs and SMCs are all stacked from multiple sheets. However, the primary distinction between UACSs and SMCs is the high structural regularity. When fibers are arranged in a regular pattern, UACSs exhibit excellent stiffness, which is comparable to continuous CFRPs, a strength that is much greater than SMCs, and high formability that is similar to SMCs. Hang Li et al. [7] also suggested new designs of UACSs by incorporating discontinuous angled slits, and their results showed better strength and flowability compared to the conventional UACS laminates.
Due to the fact that low-velocity impact is one of the major potential hazards for laminates, and UACS laminates are now being used in composite structures primarily for their strength and damage tolerance, an impact test on them is required to understand their impact behavior. Aircraft are concerned about low-velocity impact damage during taxiing, landing, and maintenance. Low-velocity impact damage to the wing skin can result from tool drops during maintenance or small pieces of ground debris hitting the wing skin when taxiing. The kind of instrument or object that strikes the structures will determine the type of impact damage that occurs. It could have a smooth, soft edge or a sharp edge. Due to accidents that occur during the manufacturing or maintenance processes, the components could lose more than half of their strength when confronted with low-velocity impact issues [8]. Conventional nondestructive testing methods are unable to find surface flaws or damages that are scarcely perceptible, as is typical in low-velocity impacts on composites. Numerous factors, including impact velocity and energy, projectile/impactor shape, stacking sequence, and boundary conditions, affect impact behavior and damage [9]. Damage modes for composite-fiber-reinforced laminates include matrix cracking, fiber breaking, and delamination. Under impact loading, these different damage modes are exhibited collectively [10]. Since most composites are fragile, plastic deformation is not the cause of the damage. Low-velocity impacts are characterized as incidents that can happen at speeds between 1 and 10 m/s based on the stiffness, mass and material qualities of the impactor [11]. A drop-weight impact testing device can simulate a variety of actual impact situations and gather thorough performance data [12].
The low-velocity impact response of composite laminates has been studied in several pieces of research using experimental observation [13], theoretical analysis [14,15,16,17,18,19] and numerical modeling. Numerical modeling is thought to be a useful tool to obtain precise damage propagation and is also ideal for parametrical research, considering the high cost of impact experiments. The progressive damage model (PDM), considering failure initiation and stiffness degradation, is the most common among the numerical theoretical models. When the failure criterion is reached in some area of the composite, the damage begins to occur, and visible or invisible cracks begin to expand, resulting in stiffness degradation and gradually reducing load capacity. Due to the correlation with mechanical parameters, the method based on fracture energy is widely used to predict the progressive damage behavior of composites, such as the equivalent displacement method, which has been well verified by a large number of experimental results [20,21,22]. In the main damage types that occur in laminates, delamination is viewed as the crucial one because it might drastically decrease the compression strength, stability and integrity of composites [23]. Considering the essentiality of this failure mechanism, researchers have used two basic techniques to anticipate delamination in laminates: the virtual crack closure technique (VCCT) [24] and the cohesive zone model (CZM) [25]. Nevertheless, the VCCT cannot anticipate fracture initiation and requires the use of an adaptive mesh reconstruction method in the delamination front. On the other hand, the CZM utilizes criteria based on strength and fracture energy to forecast damage beginning and propagation.
The interactive mixed mode serves as the foundation for the traction–separation law of cohesive elements for simulating delamination beginning and propagation. It makes sense when the through-thickness normal stress is under tension, but the delamination inhibition effect, which is highlighted by Hou et al. [26], is neglected when the normal stress is under compression. Through-thickness compression stress should be included in CZM because it might significantly increase the interfacial shear strength. In order to investigate the impact of compressive stress on the evolution of mode II damage, Li et al. [27] suggested a contact failure model to study the maximum compression stress on mode II damage propagation. This model used a separately derived parameter to link the rise in interlaminar strength properties and mode II key crack energy to the compression stress. In the quasi-static delamination modeling of composite structures, Zhang et al. [28,29] introduced an interfacial model to appreciate the interaction-force-generated friction impact on shear-induced mode II damage. In order to forecast low-speed impact damage in composite materials utilizing a cohesion contact based on surfaces and a quasi-static stress model, Zhang and Zhang [29] devised a numerically effective approach. By including interface friction in the shear strength directions, this model studied the impact of compressive interlaminar stress on deformation.
Based on the literature, there are few reports about research on UACS laminates under low-speed impact; hence, the goal of the current investigation is to determine its damage response and explore the influence of slits on the impact performance of UACS laminates with vertical slits. The damage mechanisms underlying this topic are of importance for the application of UACS composites. This triggered the analysis of the experimental and numerical study of the influence of vertical slits on the impact behavior of UACS laminates in the present work. A series of drop-weight impact experiments are conducted on the UACS composites under different impact energies. The impact response curves are analyzed, and the internal damage is examined with C-scan inspection. Parallelly, in order to model the progressive damage process of UACS composite structures, a 3D finite element simulation model is made in ABAQUS/Explicit, and the progressive damage model (PDM) is suggested and encoded in the user-specified subroutine (VUMAT), which incorporates the Hashin criterion and a computation technique for equivalent movement in the damage propagation model. To describe the effect of delamination occurring at the slits on the damage evolution, interface cohesive pieces with the nonlinear traction–separation law are inserted among plies. Finally, a discussion on the UACS laminates with vertical slits is performed for impact performance in damage tolerance and energy absorption.

2. Materials and Experimental Setup

2.1. Materials

The composite materials used in this research were UIN12500 (T800) (Guangwei, Weihai, China), a commercial continuous CFRP prepreg of 0.15 mm thickness lamina and 60% fiber volume fraction. A quasi-isotropic lay-up sequence of [0/90]4S was tested, and all the composite laminates were fabricated using 16 layers of the UACS laminates. The laminate was placed into a hot press machine at a pressure of 0.6 MPa and a compression temperature of 130 °C, according to the prepreg manufacturer’s supplied curing curve. The specimens were cut precisely using a water-jet cutting machine, and the vertical slits on the UACS prepreg were made by a commercial cutting tool. The final slit had nominal dimensions of 25 mm × 0.2 mm. Each of the UACS laminates had vertical slits, as illustrated in Figure 1. The specimen’s standard dimensions were kept at 100 mm in length and 100 mm in width. The laminates had an average thickness of 1.99 mm.

2.2. Experimental Procedure

The cross-ply UACS composite materials were put through low-velocity impact tests for various impact energies of 4 J, 7 J, and 11 J at room temperature. Figure 2a shows the impact test setup and the impactor that was employed in this study. The specimen was fastened by a clamping steel plate with a 75 mm diameter hole in the center, and the top surface’s edges were tightened after synthetic rubber pads were positioned in the gap [30]. Samples from the UACS laminates were evaluated under impact loading using an impactor with a hemispherical shape with a 16 mm diameter. The drop-weight mass of the hemispherical impactor was 5.265 kg. As depicted in Figure 2b, test data were collected using an acceleration sensor mounted to the impactor and a drop tower impact system (Instron CEAST 9350). The preset impact energies were achieved by adjusting the drop height of the impactor and the compression degree of the spring. At least three specimens were tested for each energy level according to the ASTM D7136/D7136M-20 standard. Because the repeatability of the drop-weight impact test is acceptable, only one sample for impact energy value was used [31]. Finally, the nondestructive inspection characterizations were assessed by ultrasound C-scanning systems (Olympus NDT 2.10R25).

3. Finite Element Modeling

In order to model the low-speed impact damage in UACS composite structures, a 3D finite element simulation model was created by using the ABAQUS/Explicit Dynamic simulation method. The researchers’ main priority was to analyze the impact damage problems of composite structures, such as the accurate experimental monitoring of the barely visible impact damage and the numerical simulation of impact response, for predicting the complex structural damage mechanisms in a fairly short time. The experiment was carried out at room temperature in an enclosed space; therefore, the numerical simulation did not consider the influence of ambient temperature. The 16 mm diameter hemispherical head impactor was simulated as an analytical rigid body. The model of the laminate was cut and given different material properties, as seen in Figure 3. The material of the slits in the model was defined as resin, and the other parts were defined as received CFRP prepreg [32]. All freedoms along the periphery of the plate were constrained to zero to simulate the experimental clamped conditions, and all freedoms of the indenter except the z-axis were constrained to zero to simulate the drop-weight low-velocity impact. The distance between the external surface of the impactor and the plate was set as 0 mm. The initial velocities of 1.23 m/s, 1.65 m/s and 2.04 m/s over the z-axis and 5.265 kg of mass were given to the indenter to generate three impact energies of 4, 7 and 11 J.
The eight-node solid elements with reduced integration (C3D8R) were used in the plate. Table 1 shows the details of the material properties of the prepreg used in this research. These parameters were provided by the manufacturer of prepreg (Weihai Guangwei) based on experimental analysis and empirical calculation. In addition, the zero-thickness eight-node three-dimension cohesive elements (COH3D8) were inserted between adjacent layers to simulate delamination. In order to improve the computational accuracy and efficiency, the mesh with an element size of 1 mm × 1 mm was used in the model. Therefore, the mesh of composite laminates consisted of 99216 C3D8R elements and 93015 COH3D8 elements. Regarding the interface cohesive element, Liu [33] selected the best parameters for the cohesive element in their work after comparing the impacts of various contact strengths on crack propagation damage. Table 2 lists the parameters of the interface cohesive elements that were employed in this investigation. All freedoms along the periphery of the plate were constrained to zero to simulate the experimental clamped conditions. A general contact algorithm in ABAQUS/Explicit was adopted to simulate contact between the impactor and the plate, as well as ply to ply in the laminate. The penalty method with a friction coefficient was used to describe the tangential behavior, while the hard contact method was employed to describe the normal behavior. Between the layers, there was an average friction coefficient of 0.5, and between the impactor and the plate, it was 0.3 [22,34]. Typically, the two categories used to describe the damage types of composite laminates under low-speed impact are intralaminar damage (matrix cracking and fiber breaking) and interlaminar damage. Damage simulations, including the damage initiation criterion and damage propagation law, can simulate the impact damage process. The equivalent displacement model selected in this paper has good accuracy for simulating damage development [22,30,35].

3.1. Johnson–Cook Elastic–Plastic Model

The failure mechanism of slits filled with resin is a key factor associated with initial damage growth. The elastic–plastic phase of resin is mimicked by the Johnson–Cook model, and then, the ductile damage criterion determines the fracture to describe nonlinear stress–strain behavior and failure mechanism. In the present study, the Johnson–Cook plasticity model is expressed by the following Equation:
σ 0 = A + B ε _ pl n
where ε _ pl is the equivalent plastic strain, and σ 0 is the static yield stress. A, B and n are material parameters. Referring to the curve fitting in a previous study [36], A = 35, B = 109.7 and n = 0.259 are applied for the failure analysis of resin, as shown in Figure 4.

3.2. Damage Initiation Criterion

One of the well-known standards for forecasting failure starts in composite structures is the Hashin criterion [37], which is constructed using formulas based on material strengths. The following are the four distinct Hashin criteria failure modes:
Fiber tensile failure: ( σ 11 ≥ 0)
F ft = σ 11 X T 2 + σ 12 S 12 2 + σ 13 S 13 2 1
Fiber compression failure: ( σ 11 < 0)
F fc = σ 11 X C 2 1
Matrix tensile failure: ( σ 22 + σ 33 ≥ 0)
F mt = σ 22 + σ 33 Y T 2 + 1 S 23 2 σ 23 2 σ 22   σ 33 + σ 12 S 12 2 + σ 13 S 13 2 1
Matrix compression failure: ( σ 22 + σ 33 < 0)
F mc = σ 22 + σ 33 2 S 23 2 + σ 22 + σ 33 Y C Y C 2 S 23 2 1 + 1 S 23 2 σ 23 2 σ 22 σ 33 + σ 12 S 12 2 + σ 13 S 13 2 1
The fiber direction tensile and compressive strengths are represented by the variables X T   and X C in the equations above; S 12 , S 13 and S 23 are the corresponding shear strengths, while Y T and Y C represent the tensile and compression strengths in the cross-sectional direction.

3.3. Damage Evolution Law

Once the damage initiation standard is attained, a damage progression law is needed for the damage development. The progressive damage model needs to be encoded in the user-specified subroutine, which incorporates the Hashin criterion and a computation technique for equivalent movement in the damage propagation model. Here, a gradual deterioration approach coupled with a modified numerical model suggested by Zhou et al. [22] is utilized to describe the progression of intralaminar damage. Each failure mode’s damage variable is written as:
G I = 1 2 σ I , eq ε I , eq l c
d I = δ I , eq f δ I , eq δ I , eq 0 δ I , eq δ I , eq f δ I , eq 0   ( d I ϵ [ 0 ,   1 ] ,   I = ft ,   fc ,   mt ,   mc )
where G I is the fracture toughness of different failure modes; l c is the element’s characteristic length. The damage corresponding displacement of the related potential failure at initial and final failure is represented by δ I , eq f and δ I , eq 0 , respectively.
The constitutive Equation with nine independent cross-ply-layered composite laminates treated as orthorhombic anisotropy materials was adopted to appreciate the mechanical response. The degraded flexibility matrix S d and the corresponding degraded stiffness matrix C d , added with different damage variables in the constitutive Equations, are expressed as:
S d = 1 d f E 11 ν 21 E 22 ν 31 E 33 ν 12 E 11 1 d m E 22 ν 32 E 33 ν 13 E 11 ν 23 E 22 1 E 33 1 d f d m G 12 1 d f d m G 23 1 d f d m G 31
C = 1 Δ d f E 11 1 d m ν 23 ν 32 d f d m E 11 ν 21 + ν 23 ν 31 d f E 11 ν 31 + d m ν 21 ν 32 d m E 22 1 d f ν 13 ν 31 d m E 22 ν 32 + d f ν 12 ν 31 E 33 1 d f d m ν 12 ν 21 Δ d f d m G 12   Δ d f d m G 23 Δ d f d m G 13
d f = 1 d ft 1 d fc d m = 1 S mt d mt 1 S mc d mc Δ = 1 d f d m ν 12 ν 21 d m ν 23 ν 32 d f ν 13 ν 31 2 d f d m ν 21 ν 32 ν 13
where d I ( I = ft ,   fc ,   mt ,   mc ) are the fiber and matrix damage variables computed using the damage evolution model. S mt and S mc are newly introduced coefficients that are adjusted to 0.96 in the current study in order to manage the shear stiffness loss [33,38].

3.4. Interlaminar Damage

Before damage initiation and evolution, the traction–separation model is in the elasticity phase. Interlaminar delamination occurs with the initiation of damage simulated by the quadratic nominal stress criterion, and corresponding reduction behavior is determined by the BK fracture criterion [39]. Equations (11) and (12) express the cohesive zone model used to predict the growth of delamination damage (12).
  t n t n o 2 + t s t s o 2 + t t t t o 2 = 1
G C = G n C + G s C G n C G S G T η G S = G s + G t ,   G T = G n + G S  
where G n , G s and G t are normal and shear directions fracture energies, and t n ,   t s , and t t stand for normal and shear tractions; t n o , t s o , and t t o stand for interface normal and shear strength. The quantities G n C , G s C and G C are critical fracture energy of the normal, the shear and total directions, respectively, and η , set as 1.45, is a material parameter.

4. Results and Discussion

As part of the impact test process at various energies (4, 7 and 11 J), the force–time, force–displacement and energy–time curves were documented. For numerical simulation, the analysis was carried out from damage modes, including fiber tensile damage, intralaminar matrix cracking and fatigue cracks at the interfaces.

4.1. Impact Response

Energy–time, force–displacement and force–time histories of both continuous CFRP and UACS laminates were compared for each impact energy in order to make an accurate comparison of the impact performance of the experimental data (see Figure 5). The energy–time curve of UACS laminates associated with penetration at 11 J will be discussed in the comparison of the numerical and experimental results. For both continuous CFRP and UACS laminates, 11 J impact energy presents the highest impact time, displacement and energy, whereas laminates impacted with 4 J display similar impact response behaviors due to no obvious damage occurring in the laminates under this impact energy. Table 3 lists the final results from the experiment with various impact energies. The peak force of UACS and CFRP laminates differs by 14% at 4 J, and subsequently, with the failure of UACS laminates, the increasing gap indicates a decrease in load capacity for UACS laminates. However, in the case of non-failure (4 J), the UACS laminates show better toughness with a maximum displacement of 4.48 mm compared to 3.61 mm of CFRP. Compared to the continuous CFRP laminates, the UACS laminate shows a larger maximum displacement by about 24%, and its peak force is about 37% less than continuous CFRP laminates under the impact energy of 7 J due to the fiber breakage of UACS laminates. The final energy absorbed by the CFRP laminates is about 1.99 J for the impact energy of 4 J, which is approximately 10% lower than that of the UACS (2.21 J). For each of the three types of impact energy, the impact energy–time histories are displayed in Figure 5c. For 11 J impact energy, the ultimate energy absorbed by the laminates in the event of damage and friction is 10.38 J and 5.34 J for UACS and CFRP, respectively. In the case of 7 J, however, the UACS value of 6.05 J is closer to the CFRP value of 3.29 J. It appears that absorbed energy increases with increasing impact energy, and it is clear that the UACS laminates absorb more energy than the CFRP laminates.
Figure 6 exhibits the correlation between experimental and numerical impact force–time graphs of UACS laminates at three impact energies. Figure 6 and Figure 7 show that the initial contact stage of the impact force–time curves always exhibits substantial oscillation, which can be attributed to the early fracture of the layers on the impact side. Meanwhile, sudden drops in the force and displacement under the energy of 7 J occur at 2 ms and 3.95 ms, respectively. In the beginning, the delamination grows quickly at the slits, which causes the initial drop in force. With damage accumulation due to the continuous absorption of kinetic energy by the laminate, the failure of the laminate leads to a steeper drop than the first drop. When the initial kinetic energy increases to 11 J, the drop in force is obviously earlier. The establishment of material structure in the numerical model is based on the ideal situation compared to the experiments. Therefore, the accurate simulation of the mechanical response of brittleness composite laminates is relatively difficult. Premature failure of impact side layers causes the penetration of composite laminates, and that is why the impact process of 11 J is prolonged. As a result of the actions of laminate vibration, damage initiation and progression with various modes, the impact force increases swiftly and fluctuates dynamically. As the value gets closer to the maximum, extremely intense oscillations occur, and the rate at which the impact force increases declines. The impact force then starts to decrease, indicating the impactor will bounce after that. In contrast to the experiment test’s limited acquisition points, a large amount of data sets may be obtained in brief time intervals, and simulation can effectively depict the intensive oscillation patterns of the impact force. The simulation’s rapid oscillation history depicts the laminates’ increasing damage due to several failure types. It does, however, differ somewhat from the actual experimental method. The impact force–displacement graphs for the three impact energies are shown in Figure 7. Due to the crack formation of laminates throughout the experiment, the maximum displacement speculated by the numerical model was consistently a little less than the results of the experiment. Figure 7a–c again demonstrate that the maximum displacement of the impactor increases as the impact energy does. Due to the bending deformation of composite laminates, the displacement of the indenter exceeds the thickness of the laminates by a significant margin.
The impact energy–time graphs for three impact energies are depicted in Figure 8. After contact, the impactor’s kinetic energy quickly starts to transfer to the composite laminates. The elastic energy of the laminates absorbs a portion of the kinetic energy, while the damage and vibration of the laminates disperse the remainder. An accurate numerical model is used to simulate the peak energy absorption. As shown in Figure 8b,c, the energy absorbed by laminates is released by means of fiber brittle fracture and delamination near the slits, and the indenter takes longer in the experimental test than in the numerical simulation to reach zero velocity and stable bouncing back velocity; as a result, the last absorbed energy determined by the numerical methods is never more than the results of the experiment. The laminates absorb more energy as the impact energy rises, indicating that the damage issue worsens.

4.2. Post Impact C-Scan Inspection

Additional information on the damage mechanisms is provided by the C-scan maps of UACS and continuous CFRP laminates (see Figure 9). The impact-damaged specimens are also compared to the C-scan maps. Each damaged UACS laminate is represented by either the hemispherical impactor (front side) or a fracture surface (back side) corresponding to the 4 J, 7 J and 11 J impact energies. The 4 J impact does not cause significant damage to the specimens, whereas the 11 J impact induces the maximum damage in comparison to the other energies. At 7 J, the hemispherical indentation on the impacted side becomes clear (this is the first energy that makes an indentation over the BVID in the specimens). Different from the hemispherical indentation, the delamination area of the UACS laminates is shaped like a square influenced by slits. Additionally, there is a generalization of fiber splitting and major development of warp/weft direction cracks on the non-impacted side. Damage evolution in UACS is more severe than that of CFRP laminates. Obviously, for the highest impact energy in UACS, the delaminated areas reach about 20% of the specimen, as listed in Table 3. It can be concluded that the UACS laminates tend to delaminate to a greater extent than CFRP laminates at increasing impact energy.

4.3. Failure Modes Analysis

Matrix cracking, delamination and fiber breaking are the primary damage mechanisms for composite structures in low-velocity impacts. Even with little impact energy, matrix cracking and delamination can easily form, but fiber breaking only happens with considerably large impact energy. These damage types can sometimes be invisible from the outside of the laminates, yet they nonetheless significantly affect the remaining mechanical performance of composite structures. High tension between laminas and shear stresses in the impact zone result in delamination, which is an essential component of impact damage.

4.3.1. Fiber Tensile Damage

In the initial stage of the impact process, the resin in the slits is the first to crack after the laminate is impacted. Then, the slits are completely fractured to release energy, which is one of the reasons why UACS laminates absorb more energy. Because the slit direction of adjacent layers is vertical, the fibers are subjected to a large amount of shear force when the slits are damaged, resulting in lower impact performance of the laminate. However, in this case, the interlaminar delamination induced by the notch failure will limit the damage extending from the impact center point to the outside.
As shown in Figure 10a,b, slight damage appears in the central region at 4 J impact energy; as a comparison, large-area fiber breakages are induced with the path of damage evolution perpendicular to the slits under 7 J impact energy, which will result in the penetration of UACS laminate at higher energy. Figure 10c illustrates fiber tensile damage before penetration at 12.8 ms of the impact process. The layers on the impact and back side present similar imaging with 4 J; namely, damage accumulates primitively in the center zone because this region bears impact load in the first stage and then expands along the fiber direction subjected to shear loading.

4.3.2. Matrix Tensile Damage and Compression Damage

Figure 11 depicts how matrix tensile damage is distributed across each layer of the numerical models. The damage is evidently concentrated near the place of impact. Moreover, it is clear that when the equivalent layer is farther away from the contact side, the anticipated matrix tensile damage area is bigger. This can be defined by the way the laminates deform and fail, showing how matrix tensile damage begins upon the last laminate and then spreads toward the top layers. It is noteworthy that the addition of slits changes the damage propagation rule, especially in Figure 11a,b, where the matrix damage offers the shape of a “vertical bar” in the vicinity of the slits as the threshold of the damage area observed by ultrasound C-scanning systems. Similarly, from Figure 11c, which shows the damage condition before the failure of the laminate, it can be seen that the matrix damage of back side layers expands from the central region to damaged slits. The predictions of matrix compression damage within every ply are shown in Figure 12. To be more specific, the damage region for every ply from the 7 J impact is greater than the damage region from the 4 J impact. According to this research, matrix compression damage starts on the contact side and spreads toward the backside, which is why the top four layers of the matrix show the most damage. For the majority of layers, the matrix compression damage region is less than the previously analyzed matrix tensile damage region. Given the direction of the fiber, matrix compression damage tends to concentrate at the site of impact and rarely extends in other directions.

4.3.3. Delamination of UACS Composite

For each interlaminar interface on the impact region, Figure 13 depicts the delamination morphology. For the overall value of the scalar damage variable (SDEG), variables above that indicate whether the BK fracture criterion is satisfied or not, a value that is less than 1.0 indicates that the criterion is not satisfied, while a value of 1.0 or higher indicates that the delamination occurred. In total, the UACS cross-ply composite laminates examined here have fifteen interfaces. As can be seen under the 4 J energy level, the majority of the delamination region in each interface is not entirely damaged but only partly fails. The composite structures could unexpectedly fail as a result of this partly failed delamination continuing to spread secretly during subsequent use. The largest critical delamination can be observed at the 11 J energy level followed by 7 J.
It is common knowledge that the stress concentration in CFRP laminates responding to low-velocity impacts is only high when close to the impact area and appears to reduce along the laminates’ in-plane and thickness directions. However, for the UACS laminates, the existence of slits makes the delamination generate along the slit direction of the upper layer. This special delamination area appears at UACS laminates under different energy levels, especially at 7 J. As shown in Figure 13b, the delamination area consists of two parts, the elliptical impact region loading by hemispherical punch and the slits delamination. Unlike the delamination area under the impact energy of 7 J, the slits delamination, as shown in Figure 13a, only occurs on several backside layers. With the rise in impact energy, the delamination area propagates from central to slits, as shown in Figure 13c, resulting in the penetration of UACS laminates under the impact energy of 11 J.

5. Conclusions

This work examined the dynamic mechanical behavior and damage progression of UACS laminates with vertical slits using a series of low-velocity impact tests and a PDM-based finite element simulation integrating intralaminar damage and delamination. From the numerical and experimental results, the following can be drawn:
  • At the same impact energy level, the UACS laminates with vertical slits show lower peak force than the continuous CFRP laminates. The continuous CFRP laminate shows a limited damage area, whereas the delamination of the UACS laminate offers the shape of a “vertical bar” in the vicinity of the slits, resulting in a larger damage area. With the increase in impact energy, the maximum displacement of UACS laminates is larger than that of continuous CFRP laminates and shows superior energy absorption capability.
  • The finite element model is developed by the Johnson–Cook elastic–plastic model and PDM, and good accordance is determined between the numerical simulation and experimental measurements on the peak force of the laminates, as well as the maximum displacement and the energy absorbed during the impact event. Moreover, the effect of the dimension of slits, the stacking sequence and the thickness of laminates on the low-velocity impact performance of UACS laminates could be investigated further in future work.

Author Contributions

Y.H.: investigation, methodology, writing—original draft. F.T.E.: investigation, writing—reviewing and editing. J.H.: conceptualization, project administration, funding acquisition, methodology, supervision, writing—original draft, writing—reviewing and editing. X.Z.: visualization, investigation. J.Z.: investigation. S.Z.: investigation. R.Q.: investigation. Z.C.: investigation. D.C.: methodology. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the China Postdoctoral Science Foundation (Grant nos. 2022M721590), the National key research and development project (Grant nos. 2020YFB1506101), and the National Natural Science Foundation of China (Grant nos. 51808288). The APC was funded by Grant nos. 2022M721590.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Supporting data is available on request.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

Nomenclature G mt Fracture toughness for tensile damage of matrix
G fc Fracture toughness for compressive damage of fiber
Symboldefinition G ft Fracture toughness for tensile damage of fiber
σ 11 Component stress acting at fiber direction F ft Hashin criteria for fiber failure in tension
σ 22 Component stress acting at normal to the fiber F fc Hashin criteria for fiber failure in compression
σ 33 Component stress acting at normal to the ply plane F mt Hashin criteria for matrix failure in tension
σ 12 Shear stress acting at the ply plane F mc Hashin criteria for matrix failure in compression
σ 13 Shear stress acting transverse to ply plane l c Element characteristic length
σ 23 Shear stress acting transverse to ply plane d ft Damage variable for fiber failure in tension
σ 0 Static yield stress d fc Damage variable for fiber failure in compression
E 11 Young modulus at fiber direction d mt Damage variable for matrix failure in tension
E 22 Young modulus at normal to fiber d mc Damage variable for matrix failure in compression
E 33 Young modulus at normal to the ply plane η Material parameter of the BK fracture criterion
G 12 Shear modulus acting at the ply plane S mc Coefficients to manage the shear stiffness
G 13 Shear modulus acting at the ply plane S mt Coefficients to manage the shear stiffness
G 23 Shear modulus acting transverse to ply plane δ I , eq f Damage corresponding displacement of failure at initial failure
ν 12 Poisson’s ratio acting at the ply plane δ I , eq 0 Damage corresponding displacement of failure at final failure
ν 13 Poisson’s ratio acting at the ply plane E n Elastic modulus of normal direction
ν 23 Poisson’s ratio acting transverse to ply plane E s Elastic modulus of shear direction
X T Fiber direction tensile strength E t Elastic modulus of shear direction
X C Fiber direction compressive strength t n o Interface normal strength
Y T Tensile strength in the cross sectional direction t s o Interface shear strength
Y C Compressive strength in the cross sectional direction t t o Interface shear strength
S 12 Shear strength acting at the ply plane G n c Critical fracture energy of the normal direction
S 13 Shear strength acting at the ply plane G s c Critical fracture energy of the shear direction
S 23 Shear strength acting transverse to ply plane G C Critical fracture energy of the total direction
G mc Fracture toughness for compressive damage of matrix ε _ pl Equivalent plastic strain

References

  1. Zhang, X. Impact on Composite Structures. In The Aeronautical Journal; Abrate, S., Ed.; Cambridge University Press: Cambridge, UK, 1998; Volume 103, p. 536. ISBN 0521473896. [Google Scholar] [CrossRef]
  2. Zhang, H.; Liu, M. Development and Applications of Carbon Fiber Reinforced Polymer. Eng. Plast. Appl. 2015, 43, 132–135. [Google Scholar]
  3. Friedrich, K.; Almajid, A.A. Manufacturing aspects of advanced polymer composites for automotive applications. Appl. Compos. Mater. 2013, 20, 107–128. [Google Scholar] [CrossRef]
  4. Palmer, J.S.L.; Ghita, O.R.; Evans, K.E. Sheet moulding compound (SMC) from carbon fibre recyclate. Compos. Part A Appl. Sci. Manuf. 2010, 41, 1232–1237. [Google Scholar] [CrossRef]
  5. Thomason, J.L. The influence of fibre length and concentration on the properties of glass fibre reinforced polypropylene: 7. Interface strength and fibre strain in injection moulded long fibre PP at high fibre content. Compos. Part A Appl. Sci. Manuf. 2007, 38, 210–216. [Google Scholar] [CrossRef]
  6. Taketa, I.; Okabe, T.; Kitano, A. A new compression-molding approach using unidirectionally arrayed chopped strands. Compos. Part A Appl. Sci. Manuf. 2008, 39, 1884–1890. [Google Scholar] [CrossRef]
  7. Li, H.; Wang, W.-X.; Takao, Y.; Matsubara, T. New designs of unidirectionally arrayed chopped strands by introducing discontinuous angled slits into prepreg. Compos. Part A Appl. Sci. Manuf. 2013, 45, 127–133. [Google Scholar] [CrossRef]
  8. Richardson, M.O.W.; Wisheart, M.J. Review of low-velocity impact properties of composite materials. Compos. Part A Appl. Sci. Manuf. 1996, 27, 1123–1131. [Google Scholar] [CrossRef]
  9. Shah, S.Z.H.; Karuppanan, S.; Megat-Yusoff, P.S.M.; Sajid, Z. Impact resistance and damage tolerance of fiber reinforced composites: A review. Compos. Struct. 2019, 217, 100–121. [Google Scholar] [CrossRef]
  10. 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]
  11. Shivakumar, K.; Elber, W.; Illg, W. Prediction of low-velocity impact damage in thin circular laminates. AIAA J. 1985, 23, 442–449. [Google Scholar] [CrossRef]
  12. Mathivanan, N.R.; Jerald, J.J. Experimental investigation of low-velocity impact characteristics of woven glass fiber epoxy matrix composite laminates of EP3 grade. Mater. Des. 2010, 31, 4553–4560. [Google Scholar] [CrossRef]
  13. Ibrahim, G.R.; Albarbar, A.; Brethee, K.F. Damage degradation modelling for transverse cracking in composite laminates under low-velocity impact. Eng. Fract. Mech. 2022, 263, 108286. [Google Scholar] [CrossRef]
  14. Mishra, A.; Naik, N.K. Failure initiation in composite structures under low-velocity impact: Analytical studies. Compos. Struct. 2010, 92, 436–444. [Google Scholar] [CrossRef]
  15. Xiao, S.; Chen, P.; Ye, Q. Prediction of damage area in laminated composite plates subjected to low velocity impact. Compos. Sci. Technol. 2014, 98, 51–56. [Google Scholar] [CrossRef]
  16. Olsson, R. Analytical prediction of damage due to large mass impact on thin ply composites. Compos. Part A Appl. Sci. Manuf. 2015, 72, 184–191. [Google Scholar] [CrossRef]
  17. Bui, T.Q.; Hirose, S.; Zhang, C.; Rabczuk, T.; Wu, C.-T.; Saitoh, T.; Lei, J. Extended isogeometric analysis for dynamic fracture in multiphase piezoelectric/piezomagnetic composites. Mech. Mater. 2016, 97, 135–163. [Google Scholar] [CrossRef]
  18. Singh, H.; Mahajan, P. Analytical modeling of low velocity large mass impact on composite plate including damage evolution. Compos. Struct. 2016, 149, 79–92. [Google Scholar] [CrossRef]
  19. Yu, Z.; Gao, S. Increase of contact radius due to deflection in low velocity impact of composite laminates and prediction of delamination threshold load. Compos. Struct. 2016, 147, 286–293. [Google Scholar] [CrossRef]
  20. Demiral, M.; Abbassi, F.; Saracyakupoglu, T.; Habibi, M. Damage analysis of a CFRP cross-ply laminate subjected to abrasive water jet cutting. Alex. Eng. J. 2022, 61, 7669–7684. [Google Scholar] [CrossRef]
  21. Demiral, M.; Kadioglu, F.; Silberschmidt, V.V. Size effect in flexural behaviour of unidirectional GFRP composites. J. Mech. Sci. Technol. 2020, 34, 5053–5061. [Google Scholar] [CrossRef]
  22. Zhou, J.; Wen, P.; Wang, S. Finite element analysis of a modified progressive damage model for composite laminates under low-velocity impact. Compos. Struct. 2019, 225, 111113. [Google Scholar] [CrossRef]
  23. Feng, D.; Aymerich, F. Finite element modelling of damage induced by low-velocity impact on composite laminates. Compos. Struct. 2014, 108, 161–171. [Google Scholar] [CrossRef]
  24. Krueger, R. Virtual crack closure technique: History, approach, and applications. Appl. Mech. Rev. 2004, 57, 109–143. [Google Scholar] [CrossRef]
  25. Camanho, P.P.; Davila, C.G.; de Moura, M.F. Numerical simulation of mixed-mode progressive delamination in composite materials. J. Compos. Mater. 2003, 37, 1415–1438. [Google Scholar]
  26. Hou, J.P.; Petrinic, N.; Ruiz, C.; Hallett, S.R. Prediction of impact damage in composite plates. Compos. Sci. Technol. 2000, 60, 273–281. [Google Scholar] [CrossRef]
  27. Li, X.; Hallett, S.R.; Wisnom, M.R. Predicting the effect of through-thickness compressive stress on delamination using interface elements. Compos. Part A Appl. Sci. Manuf. 2008, 39, 218–230. [Google Scholar] [CrossRef]
  28. Zhang, X.; Bianchi, F.; Liu, H. Predicting low-velocity impact damage in composites by a quasi-static load model with cohesive interface elements. Aeronaut. J. 2012, 116, 1367–1381. [Google Scholar] [CrossRef]
  29. 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]
  30. Shi, Y.; Pinna, C.; Soutis, C. Modelling impact damage in composite laminates: A simulation of intra- and inter-laminar cracking. Compos. Struct. 2014, 114, 10–19. [Google Scholar] [CrossRef]
  31. Falcó, O.; Lopes, C.S.; Sommer, D.E.; Thomson, D.; Ávila, R.; Tijs, B. Experimental analysis and simulation of low-velocity impact damage of composite laminates. Compos. Struct. 2022, 287, 115278. [Google Scholar] [CrossRef]
  32. Hu, J.; Zhang, X.; Chen, Z.; Guo, W.; Li, H.; Deng, X. Experimental and Numerical Research on Open-Hole Strength and Damage Mechanism of Regularly Arrayed Short Fiber Reinforced Polymer Composite. Polymers 2020, 12, 1622. [Google Scholar] [CrossRef]
  33. 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]
  34. 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]
  35. Xin, S.H.; Wen, H.M. A progressive damage model for fiber reinforced plastic composites subjected to impact loading. Int. J. Impact Eng. 2015, 75, 40–52. [Google Scholar] [CrossRef]
  36. Fuller, J.D.; Wisnom, M.R. Pseudo-ductility and damage suppression in thin ply CFRP angle-ply laminates. Compos. Part A Appl. Sci. Manuf. 2015, 69, 64–71. [Google Scholar] [CrossRef]
  37. Hashin, Z.; Rotem, A. A Fatigue Failure Criterion for Fiber Reinforced Materials. J. Compos. Mater. 1973, 7, 448–464. [Google Scholar] [CrossRef]
  38. Zhang, C.; Li, N.; Wang, W.; Binienda, W.K.; Fang, H. Progressive damage simulation of triaxially braided composite using a 3D meso-scale finite element model. Compos. Struct. 2015, 125, 104–116. [Google Scholar] [CrossRef]
  39. 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, 56, 439–449. [Google Scholar] [CrossRef]
Figure 1. Unidirectional CFRP prepreg with vertical, chopped slits.
Figure 1. Unidirectional CFRP prepreg with vertical, chopped slits.
Polymers 15 00840 g001
Figure 2. Impact setup on UACS laminate using the hemispherical impactor. (a) Schematic diagram. (b) Instron CEAST 9350 Drop Tower Impact System 3. Finite element modeling.
Figure 2. Impact setup on UACS laminate using the hemispherical impactor. (a) Schematic diagram. (b) Instron CEAST 9350 Drop Tower Impact System 3. Finite element modeling.
Polymers 15 00840 g002
Figure 3. Model of UACS composite laminates under low-velocity impact.
Figure 3. Model of UACS composite laminates under low-velocity impact.
Polymers 15 00840 g003
Figure 4. Stress-equivalent plastic strain curve of matrix resin [36].
Figure 4. Stress-equivalent plastic strain curve of matrix resin [36].
Polymers 15 00840 g004
Figure 5. Comparison of experimental result of CFRP and UACS laminates under 4 J, 7 J and 11 J impact energies. (a) Force–time curve, (b) force–displacement curve and (c) energy–time curve.
Figure 5. Comparison of experimental result of CFRP and UACS laminates under 4 J, 7 J and 11 J impact energies. (a) Force–time curve, (b) force–displacement curve and (c) energy–time curve.
Polymers 15 00840 g005
Figure 6. Comparison of numerical and experimental impact force–time curves of UACS laminates under impacts at (a) 4 J, (b) 7 J and (c) 11 J.
Figure 6. Comparison of numerical and experimental impact force–time curves of UACS laminates under impacts at (a) 4 J, (b) 7 J and (c) 11 J.
Polymers 15 00840 g006
Figure 7. Comparison of numerical and experimental force–displacement curve of UACS laminates under impacts (a) 4 J, (b) 7 J and (c) 11 J.
Figure 7. Comparison of numerical and experimental force–displacement curve of UACS laminates under impacts (a) 4 J, (b) 7 J and (c) 11 J.
Polymers 15 00840 g007
Figure 8. Comparison of numerical and experimental energy–time curves of UACS laminates under impacts (a) 4 J, (b) 7 J and (c) 11 J.
Figure 8. Comparison of numerical and experimental energy–time curves of UACS laminates under impacts (a) 4 J, (b) 7 J and (c) 11 J.
Polymers 15 00840 g008
Figure 9. Comparison of C-scan maps and damaged specimens impacted at 4 J, 7 J and 11 J. (a) Damaged UACS laminates, (b) C-scan maps of UACS laminates and (c) C-scan maps of continuous CFRP laminates.
Figure 9. Comparison of C-scan maps and damaged specimens impacted at 4 J, 7 J and 11 J. (a) Damaged UACS laminates, (b) C-scan maps of UACS laminates and (c) C-scan maps of continuous CFRP laminates.
Polymers 15 00840 g009
Figure 10. Fiber tensile damage of UACS laminates under impacts at (a) 4 J, (b) 7 J and (c) 11 J.
Figure 10. Fiber tensile damage of UACS laminates under impacts at (a) 4 J, (b) 7 J and (c) 11 J.
Polymers 15 00840 g010aPolymers 15 00840 g010b
Figure 11. Matrix tensile damage of UACS laminates under impacts at (a) 4 J, (b) 7 J and (c) 11 J.
Figure 11. Matrix tensile damage of UACS laminates under impacts at (a) 4 J, (b) 7 J and (c) 11 J.
Polymers 15 00840 g011aPolymers 15 00840 g011b
Figure 12. Matrix compression damage of UACS laminates under impacts at (a) 4 J, (b) 7 J and (c) 11 J.
Figure 12. Matrix compression damage of UACS laminates under impacts at (a) 4 J, (b) 7 J and (c) 11 J.
Polymers 15 00840 g012aPolymers 15 00840 g012b
Figure 13. Delamination area of the UACS composite laminates under impacts at (a) 4 J, (b) 7 J and (c) 11 J.
Figure 13. Delamination area of the UACS composite laminates under impacts at (a) 4 J, (b) 7 J and (c) 11 J.
Polymers 15 00840 g013aPolymers 15 00840 g013b
Table 1. Material properties of continuous CFRP laminate.
Table 1. Material properties of continuous CFRP laminate.
ParametersValue
Density1600 kg/m3
Stiffness properties E 11   = 145 Gpa, E 22 = E 33 = 9 Gpa, G 12 = G 13 = 4.4 Gpa, G 23 = 3.7 Gpa, ν 12 = ν 13 = 0.3, ν 23 = 0.4
Strength properties X T = 2600 Mpa, X C = 1050 Mpa, Y T = 62 Mpa, Y C = 192 Mpa,
S 12 = S 13 = 90 Mpa, S 23 = 52 MPa
Fracture Energy G ft = 82 N/mm, G fc = 70 N/mm, G mt = 0.22 N/mm, G mc = 1.1 N/mm
Table 2. Material properties of interface cohesive elements.
Table 2. Material properties of interface cohesive elements.
ParametersValue
Elastic modulus E n = E s = E t = 5 GPa/mm
Strength propertiesN = S = T = 30 MPa
Fracture energy G n c = 0.32 N/mm; G s c = 1 N/mm
Relevant coefficient η 1.45
Table 3. Experimental values of different specimens.
Table 3. Experimental values of different specimens.
SpecimenImpact
Energy (J)
Peak Force (kN)Maximum
Displacement (mm)
Absorbed Energy (J)Damage Area (mm2)
CFRP42.33 ± 0.083.51 ± 0.181.89 ± 0.03153 ± 14.50
73.35 ± 0.114.49 ± 0.233.19 ± 0.13250.5 ± 25.25
114.56 ± 0.215.53 ± 0.465.24 ± 0.08311 ± 22.05
UACS41.97 ± 0.094.18 ± 0.302.31 ± 0.08485.75 ± 41.75
72.13 ± 0.156.86 ± 0.586.05 ± 0.38889.75 ± 56.05
111.46 ± 0.1110.79 ± 0.7110.38 ± 0.241230 ± 87.85
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

Huang, Y.; EShun, F.T.; Hu, J.; Zhang, X.; Zhao, J.; Zhang, S.; Qian, R.; Chen, Z.; Chen, D. Research on Low-Velocity Impact Response of Novel Short-Fiber-Reinforced Composite Laminates. Polymers 2023, 15, 840. https://doi.org/10.3390/polym15040840

AMA Style

Huang Y, EShun FT, Hu J, Zhang X, Zhao J, Zhang S, Qian R, Chen Z, Chen D. Research on Low-Velocity Impact Response of Novel Short-Fiber-Reinforced Composite Laminates. Polymers. 2023; 15(4):840. https://doi.org/10.3390/polym15040840

Chicago/Turabian Style

Huang, Yinyuan, Felix Thompson EShun, Junfeng Hu, Xutong Zhang, Jianping Zhao, Siqi Zhang, Rui Qian, Zhou Chen, and Dingding Chen. 2023. "Research on Low-Velocity Impact Response of Novel Short-Fiber-Reinforced Composite Laminates" Polymers 15, no. 4: 840. https://doi.org/10.3390/polym15040840

APA Style

Huang, Y., EShun, F. T., Hu, J., Zhang, X., Zhao, J., Zhang, S., Qian, R., Chen, Z., & Chen, D. (2023). Research on Low-Velocity Impact Response of Novel Short-Fiber-Reinforced Composite Laminates. Polymers, 15(4), 840. https://doi.org/10.3390/polym15040840

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