Next Article in Journal
Compartmentalized Quantitative Analysis of Concrete Sulfate-Damaged Area Based on Ultrasonic Velocity
Next Article in Special Issue
Core–Shell Rubber Nanoparticle-Modified CFRP/Steel Ambient-Cured Adhesive Joints: Curing Kinetics and Mechanical Behavior
Previous Article in Journal
Municipal Sewage Sludge as a Source for Obtaining Efficient Biosorbents: Analysis of Pyrolysis Products and Adsorption Tests
Previous Article in Special Issue
Study on Optimization of Drilling Parameters for Laminated Composite Materials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on Dynamic Mechanical Properties of Carbon Fiber-Reinforced Polymer Laminates at Ultra-Low Temperatures

1
State Key Laboratory of Explosion Science and Technology, Institute of Advanced Structure Technology, Beijing Institute of Technology, Beijing 100081, China
2
Beijing Institute of Aerospace Systems Engineering, Beijing 100076, China
*
Authors to whom correspondence should be addressed.
Materials 2023, 16(7), 2654; https://doi.org/10.3390/ma16072654
Submission received: 3 March 2023 / Revised: 17 March 2023 / Accepted: 21 March 2023 / Published: 27 March 2023
(This article belongs to the Special Issue Carbon Fiber Reinforced Polymers (2nd Edition))

Abstract

:
This study uses experimental methods, theoretical research, and numerical prediction to study the dynamic mechanical properties and damage evolution of CFRP laminates at ultra-low temperatures. Based on the Split Hopkinson Pressure Bar (SHPB) device, we set up an ultra-low temperature dynamic experimental platform with a synchronous observation function; the dynamic mechanical properties of laminates were tested, and the damage evolution process was observed. The experimental results are as follows: The compression strength and modulus increase linearly with the increase in strain rate and show a quadratic function trend of increasing and then decreasing with the decrease in temperature. The damage degree of the dynamic bending sample increases obviously with the impact velocity and decreases first and then increases with the decrease in temperature. Based on the low-temperature dynamic damage constitutive, failure criterion, and interlayer interface damage constitutive of the laminates, a numerical model was established to predict the dynamic mechanical properties and damage evolution process of CFRP laminates at ultra-low temperatures, and the finite element analysis (FEA) results are consistent with the experimental results. The results of this paper strongly support the application and safety evaluation of CFRP composites in extreme environments, such as deep space exploration.

1. Introduction

Due to their remarkable advantages, such as high specific strength, elevated impact resistance, light weight, and low-cost manufacturing, CFRP composites have been widely used in critical fields, including aerospace, the automotive industry, and rail transit in recent years [1,2,3]. With the development of deep space exploration projects such as the Lunar and Mars exploration, CFRP composites used in deep space exploration face complex work environments involving ultra-low temperatures and dynamic impact [4,5]. The mechanical properties of CFRP composites under ultra-low temperatures and dynamic impact differ from those under room temperature and static environments. There may be some problems under ultra-low temperatures, such as matrix toughness deterioration and thermal expansion mismatch at the internal interface. The high strain rate impact will also cause great changes in material properties and failure forms. Extreme environmental temperature and load conditions make the service security of CFRP composites face great challenges in deep space exploration projects. Therefore, studying the dynamic mechanical properties and damage evolution of CFRP composites at ultra-low temperatures is very important for future applications.
Under ultra-low temperatures, researchers conducted static mechanical experiments on the CFRP laminates, such as tensile, compression, and three-point bending [6,7,8,9]. They found that the interface between the fibers and matrix is strengthened at low temperatures, which improves the static mechanical properties of the laminates. As a relatively mature loading method, the SHPB was used to test the dynamic mechanical properties of CFRP composites as early as 1980 [10]. Since then, researchers [11,12,13] have studied the dynamic mechanical properties of laminates using the SHPB and found that the dynamic compression shows a noticeable strain rate effect, and the compression strength and modulus are improved, to a certain extent, at high strain rates. Gomezdelrio et al. [14] study the dynamic tensile behavior of CFRP laminates at low temperatures, and the results show that temperature and strain rate have little influence on the tensile strength along the fiber direction, but the transverse strength increases obviously at low temperatures and high strain rates. Jia et al. [9] study the dynamic bending mechanical properties of CFPR laminates at 100   ° C , and the results show that the laminates can withstand a greater impact load at low temperatures. It is worth noting that some researchers [15,16] used a high-speed camera to observe the damage evolution process of samples during dynamic compression experiments of composites. In general, there are few studies on the mechanical properties of CFRP laminates considering both temperature and strain rate effects, and there is also a lack of observation methods for the full-field deformation of samples at low temperatures.
The prediction of mechanical properties and the damage evolution of composites using constitutive models and FEA is an essential direction for future research. Wysmulski et al. [17,18,19] established a numerical model to predict the stability, load-carrying capacity, and damage evolution of laminates and their structural members, and the numerical and experimental results show satisfactory agreement. Some researchers introduce strain rate and temperature effect into the classical constitutive model of laminates to describe the response of CFRP laminates and use the cohesive zone model to simulate the interlayer interaction [20,21,22]. Ren et al. [23,24] developed a multi-scale finite element model for predicting the failure of composites at low temperatures, but this microstructure-based finite element model is expensive to calculate. Currently, the proposed constitutive model does not consider the strain rate and temperature effects simultaneously, resulting in a significant error between the FEA and experimental results. Although carbon fiber is independent of strain rate and temperature [25,26], experimental results show that laminates have noticeable temperature and strain rate effects. Therefore, the numerical study of the dynamic mechanical properties and damage evolution of CFRP laminates at low temperatures requires a constitutive model that considers strain rate and temperature effects.
In this paper, the dynamic mechanical properties and the damage evolution of CFRP laminates at ultra-low temperatures are studied using experimental methods, theoretical research, and numerical prediction. Firstly, based on the SHPB device, we set up an ultra-low temperature dynamic experimental platform with a synchronous observation function; the dynamic mechanical properties of CFRP laminates were tested, and the damage evolution process was observed. Secondly, the low-temperature dynamic constitutive model of CFRP laminates was established by introducing the strain rate and temperature effect into the classical constitutive model of laminates. Then, we also considered the 3D Hashin criterion and the interlayer interface damage constitutive based on the bilinear cohesive zone model. Finally, according to the theoretical model mentioned above, the numerical prediction model was established to predict the dynamic mechanical properties and damage evolution process of CFRP laminates at ultra-low temperatures, and the FEA results are consistent with the experimental results.

2. Experiment

2.1. Experimental Platform

Based on the SHPB device, we independently set up an ultra-low temperature dynamic experimental platform with a synchronous observation function. The schematic diagram of some critical equipment of the platform is shown in Figure 1. It mainly includes three parts: the SHPB device (high-pressure air pump, Hopkinson pressure bar, and digital signal collector), the ultra-low temperature system (cryogenic box, temperature controller, and liquid nitrogen pressure tank), and the synchronous observation system (synchronous signal generator, high-speed camera, and light source). The solenoid valve connecting the temperature controller controls the amount of liquid nitrogen supplied to the cryogenic box, and the temperature in the cryogenic box can be adjusted. When the real-time temperature of the cryogenic box reaches the set value, put the samples into the cryogenic box and ensure that each sample is kept at the set temperature for more than 30 min. The SHPB device is activated for dynamic impact, and the high-speed camera is triggered to shoot the whole impact process. The sample will be sprayed with speckles on the photographed side for deformation field analysis.

2.2. Material Preparation and Experimental Scheme

The quasi-isotropic CFRP laminates studied in this paper are prepared from carbon fiber produced by Toray Corporation of Japan and TDE-85 epoxy resin made in China using autoclave process [7]. The specific preparation process is as follows: First, cut the prepreg and place it in the mold according to the layup direction [0/45/90/−45/0/45/90/−45/0/45/90/−45]2s. Then, vacuum package the mold containing prepreg. Finally, push the mold into the autoclave and cure it according to the specific curing curve (the maximum curing temperature is 180   ° C , and the maximum pressure is 0.6   MPa in the autoclave). Some studies have pointed out that the laminate samples used in dynamic compression experiments are cylinders or cuboids, which has little influence on the final experimental results, but the length–diameter ratio is suitable in the range of 0.5 ~ 2.0 . Considering the processing difficulty, the thickness of the original laminates, and the diameter of the Hopkinson pressure bar, the dynamic compression sample is finally processed into a cuboid of 10   mm × 10   mm × 8   mm , and the surface of each sample needs to be polished to ensure it is smooth and parallel. According to the requirements of the composite material three-point bending test standard GB/T1449-2005 [27] and considering the small size of the cryogenic box and the support of the Hopkinson bending bar, the dynamic bending sample is finally processed into a cuboid of 60   mm × 10   mm × 2   mm . The shapes of the two samples are shown in Figure 2.
CFRP composites for deep space exploration need to operate in ultra-low temperatures, such as 233   ° C ~ 123   ° C on the surface of the moon and 183   ° C ~ 127   ° C on the surface of Mars, and the self-designed cryogenic box system can only achieve the lowest temperature of 180   ° C , so we chose three temperatures ( RT , 80   ° C ,   and 180   ° C ) to study the mechanical properties of CFRP laminates at low temperatures. The dynamic compression experiments of three strain rates ( 800 / s ,   1200 / s ,   and   1600 / s ) and three temperatures ( RT , 80   ° C ,   and 180   ° C ) were carried out for each group of dynamic compression samples. The dynamic bending experiments of three impact velocities (the impact velocity was selected according to three damage degrees of the sample: slight damage, obvious damage, and complete fracture) and three temperatures ( RT , 80   ° C , and   180   ° C ) were carried out for each group of dynamic bending samples. Each group of experiments was repeated three times to prevent excessive deviations or data errors and ensure the reliability of the experimental results.

3. Simulation

3.1. Low-Temperature Dynamic Constitutive Model

There is a weak coupling between temperature effects and strain rate effects of CFRP composites [9], indicating that the effects of strain rate and temperature can be calculated separately and superimposed. Research shows that the mechanical properties of resin matrix are sensitive to temperature, which is the main reason for the strength and modulus of laminates changing with temperature [28,29]. Assuming that the single-layer plate is orthotropic, the relationship between its mechanical properties in each direction and temperature can be described by the following linear expressions:
E i T = k i T + b i E i , i = 1 , 2 , 3
Y j T = k j T + b j Y j ,     j = t , c
where k and b are the material parameters that can be confirmed by experiment, E i is the elasticity modulus in different directions of the single-layer plate, and Y j is strength perpendicular to the fiber direction. t and c stand for tensile and compression, respectively. Xie et al. [30] study the dynamic compression mechanical properties of CFRP laminates under one-dimensional stress using SHPB devices. The stress–strain curves of the material at different strain rates were obtained, and the variation rule of the material strength with the strain rate was given. As our material system is the same as in the work by Xie et al., the coefficient of strain rate effect of composite is obtained
η D = 1 + A ε ˙ ε 0 ˙ n ,     i f   ε ˙ > ε 0 ˙           1         ,       o t h e r w i s e
where ε ˙ is the current strain rate of the sample, ε 0 ˙ is the reference strain rate, and A and n are the parameters fitted by the experimental results. The modulus and strength of the single-layer plate considering the strain rate effect are as follows:
E i D = η D E i E i , i = 1 , 2 , 3
Y j D = η D Y j Y j ,     j = t , c
Finally, the expressions of modulus and strength of the single-layer plate considering temperature effect and strain rate effect are
E i D T = η D E i k i T + b i E i , i = 1 , 2 , 3
Y j D T = η D Y j k j T + b j Y j ,     j = t , c
The corresponding stiffness matrix C D T can be calculated by following Equation (5):
C D T = 1 a E 1 D T 1 υ 23 υ 32 E 2 D T υ 12 υ 23 υ 32 E 3 D T υ 13 υ 23 υ 32 E 1 D T υ 21 υ 31 υ 23 E 2 D T 1 υ 13 υ 31 E 3 D T υ 23 υ 21 υ 13 E 1 D T υ 31 υ 21 υ 32 E 2 D T υ 32 υ 12 υ 31 E 3 D T 1 υ 12 υ 21 a G 12 a G 23 a G 13
where a = 1 υ 12 υ 21 υ 23 υ 32 υ 13 υ 31 2 υ 12 υ 23 υ 31 .
When the failure condition is reached, damage will occur, and the mechanical properties of the single-layer plate will decrease with the accumulation of damage. This study used different damage variables to control the stiffness degradation caused by different failure modes. Three damage variables d f , d m , and d s were introduced, where d f represents the damage along the fiber direction, d m is damage perpendicular to the fiber direction, and d s denotes the damage of shear. The damaged stiffness matrix can be calculated according to these damage variables:
C d D T = 1 a 1 d f C 11 D T 1 d s C 12 D T 1 d f C 13 D T 1 d m C 22 D T 1 d s C 23 D T C 33 D T 1 d s C 44 D T C 55 D T C 66 D T

3.2. Failure Criterion and Damage Evolution

3D Hashin criterion can consider various failure modes of different components of composites [31,32,33]. For carbon fiber tension failure   σ 11 0 :
f f t = σ 11 X t D T 2 + τ 12 S 12 2 + τ 31 S 31 2 1
For carbon fiber compression failure   ( σ 11 < 0 ) :
f f c = σ 11 X c D T 2 1
For resin matrix tension failure   σ 22 + σ 33 0 :
f m t = σ 22 + σ 33 Y t D T 2 + τ 23 2 σ 22 σ 33 S 23 2 + τ 12 S 12 2 + τ 31 S 31 2 1
For resin matrix compression failure   ( σ 22 + σ 33 < 0 ) :
f m c = σ 22 + σ 33 Y c D T Y c D T 2 S 23 2 1 + σ 22 + σ 33 2 S 23 2 + τ 23 2 σ 22 σ 33 S 23 2 + τ 12 S 12 2 + τ 31 S 31 2 1
where X t D T and X c D T , respectively, are tensile and compressive strengths along the fiber direction at different temperatures and strain rates, Y t D T and Y c D T are transverse tensile and compressive strengths at different temperatures and strain rates, and S 12 , S 23 ,   and   S 31 are shear strengths in each direction. The material parameters take into account the influence of temperature and strain rate.
The damage evolution laws degrading the stiffness matrix will be calculated when the stress fields satisfy the 3D Hashin criterion. The damage variables can be expressed as follows:
d f = 1 1 d f t 1 d f c d m = 1 1 d m t 1 d m c d s = 1 1 d f 1 d m
where d f t and d f c are the damage variables of fiber tensile and compression, respectively, and d m t and d m c are the damage variables of matrix tensile and compression, respectively.
d f i = 1 1 f f i exp 1 f f i C 11 ε 11 i L c W f i ,       i = t , c d m i = 1 1 f m i exp 1 f m i C 22 ε 22 i L c W m i ,       i = t , c
where L c is the characteristic length of the mesh to prevent mesh sensitivity problems, and W f i   and   W m i are the fracture toughness values of the fiber and matrix, respectively.

3.3. Bilinear Cohesive Zone Model for Interface

The interlayer interface of CFRP laminates is characterized by the bilinear cohesive zone model. The damaged constitutive relation is defined as the expressions below:
t i = 1 d K i δ i   ,     i = n , s , t
where d is the damage variable, t n is the normal traction, t s   and   t t are the shear tractions, δ n is the normal separation, and δ s   and   δ t are the shear separations. K n is the interfacial normal stiffness, and K s   and   K t are the interfacial shear stiffness. It should be noted that normal compression does not lead to interface damage. Here, the quadratic nominal stress criterion is applied to characterize the damage initiation as follows [31]:
σ n σ n 0 2 + τ s τ s 0 2 + τ t τ t 0 2 = 1
where σ n 0 ,   τ s 0 ,   and   τ t 0 are normal, transverse shear, and longitudinal shear stresses, with σ n = max   0 , σ n . Once the damage of the interface is activated, the stiffness will gradually reduce to zero linearly, and the damage variable d increases linearly from zero to one. Here, Benzeggagh–Kenane (BK) damage evolution based on fracture energy is employed [34]:
G C = G n C + G s C G n C 2 G s G n + 2 G s η B K
where G n and G s are the monitored fracture energies in the normal and shear directions, respectively, G n C and G s C are the corresponding critical fracture energies, and η B K is the Benzeggagh–Kenane parameter.

3.4. Finite Element Modeling

FEA is widely used to evaluate the mechanical properties and failure behavior of fiber-reinforced composites because it can provide full-field stress and strain distribution at a low cost. According to the actual size of the dynamic compression and bending samples, the commercial analysis software ABAQUS 2017 was used to establish the finite element model of laminates. As shown in Figure 3, the dynamic compression and bending samples were 48   and 12 layers, respectively, and the thickness of the single-layer plate was 0.167   mm . Each single-layer plate was assigned corresponding material properties and orientation. The interlayer interface property can be assigned using a bilinear cohesive zone model, and the quadratic nominal stress criterion was selected as the interface damage criterion. Explicit dynamic analysis was used to conduct the simulation of the dynamic compression and bending of the samples. It is worth noting that the Hopkinson pressure bar device was modeled in a 1:1 ratio to restore the actual dynamic loading environment and the boundary conditions. The material property of the Hopkinson pressure bar device (including bullet, incident bar, transmission bar, etc.) was ρ = 7.9   g / c m 3 ,   E = 210   G P a ,   υ = 0.3 ,   and   c 0 = 5155   m / s . In the explicit dynamic analysis module, the initial velocity of the bullet was given according to the real velocity measured ( 10   m / s 800 / s , 15   m / s 1200 / s , and 20   m / s 1600 / s for dynamic compression; 7   m / s , 8.5   m / s , and 10   m / s for dynamic bending), and the simulation time corresponded exactly to the real experimental time. The contact between the sample surface and the end face of the pressure bar adopted surface-to-surface contact. We considered that the temperature does not change throughout the impact process and simply set the initial temperature field in the “Predefined Field” tool of ABAQUS. These models were meshed using linear reduced integration solid elements (C3D8R) of reduced integration technique. The mesh independence verification was carried out, and it was found that the influence of the current mesh density on the calculated results can be ignored.
According to the static mechanical experimental results of the single-layer plate at different temperatures [7], the temperature-related material parameters in the above equations can be obtained: k 1 = 4.52 × 10 4 ,     k 2 = k 3 = 2.96 × 10 3 ,     k t = k c = 8.98 × 10 4 , b 1 = 0.865 , b 2 = b 3 = 1.86 ,   and   b t = b c = 0.8 ; the static material properties of single-layer plate at room temperature and the tensile properties at different temperatures are shown in Table 1 and Table 2. When the reference strain rate was defined as 0.001 / s , the strength of the single-layer plate in different directions and the corresponding material parameters are shown in Table 3 [30]. The cohesive zone model property of each composite was different. The parameters of the bilinear cohesive zone model at room temperature were determined by Ref. [35] and our experimental test data, as shown in Table 4.
Based on the proposed dynamic constitutive model introducing temperature and strain rate effects, user subroutine VUMAT was developed and implemented in FEA models containing the failure criterion, damage evolution, and interlayer interface damage constitutive. The classical step-by-step iterative method in strain incremental forms was used. After each incremental step, the stress states of the elements were updated and examined under the failure criterion. When the failure criterion was met, the damage variable was called, and the element stiffness degenerated to describe the occurrence of damage. The calculation flowchart of ABAQUS calling user subroutine VUMAT is shown in Figure 4.

4. Results and Discussion

4.1. Dynamic Compression Experimental Results

Figure 5a shows the damage evolution process captured by the high-speed camera and corresponding DIC analysis results of the dynamic compression sample of CFRP laminates during compression at a 1200 / s strain rate. The whole damage evolution process of the sample from the beginning of the compression load to delamination can be observed in Figure 5a, including the strain distribution, cracking time, and delamination position of the sample surface. At 50   μ s , a large strain appears in the lower left corner, indicating that the crack may be initiated in one layer; at 60   μ s , the crack growth between the layers accelerates; at 70   μ s , the layers on both sides of the crack are entirely separated. Similar damage processes also occur in other samples at different temperatures and strain rates, but the initiation time, growth speed, and expansion degree of the cracks are different, and the final damage morphology of the other samples is shown in Figure 5b.
According to the original data of the dynamic compression experiment, the equations provided by Ref. [36] can calculate the strain, strain rate, and stress of the dynamic compression sample. Figure 6 shows the stress–strain curves of the dynamic compression sample at different strain rates and temperatures. The stress increases approximately linearly with the increase in strain at first. When the critical strain is exceeded, the samples lose their load-carrying capacity due to damage, and the stress plummets after reaching the peak stress. In addition, we can only find that strain rates and temperatures significantly affect compression strength and modulus in Figure 6, and it is difficult to see the variation trend of vital mechanical properties (compressive strength and modulus) with temperature and strain rate. Therefore, we calculated two mechanical properties (compressive strength and modulus) according to the stress–strain curve and plotted their variation trend with temperature and strain rate in Figure 7.
As shown in Figure 7a,b, the compression strength and modulus increase linearly with the increase in strain rate, indicating that CFRP laminates exhibit a noticeable strain rate effect during dynamic compression. The mechanism of the strain rate effect can be found in Figure 5b; the failure degree of the dynamic compression samples increases with the increase in strain rate at room temperature, indicating that more impact energy can be absorbed, thus improving the compressive strength of the samples. As shown in Figure 7c,d, the compression strength and modulus show a quadratic function variation trend of first increasing and then decreasing with the increase in temperature. In order to further analyze the temperature effect and damage mechanism, the morphology at the damage location of the samples after dynamic compression was characterized using scanning electron microscopy (SEM). The characterization results of Figure 8 are as follows: At room temperature, the surface of the damage location is smooth, and only a small amount of matrix can be seen, indicating that the interlayer bonding is weak and delamination easily occurs. At 80   ° C , a single fiber and the matrix can be seen; this is because the contraction of the matrix at a low temperature enhances the interface bonding between the matrix and the fiber, which not only causes the sample to delaminate but also causes some fibers to suffer compression damage, and the interlayer properties are enhanced. At 180   ° C , there is more matrix adhesion to the fiber, and the matrix suffers compression damage together with the fiber, indicating that the interface bonding between the matrix and the fiber is enhanced more obviously. However, the decrease in temperature makes the toughness of the matrix deteriorate, and its influence gradually begins to be greater than the increase in the interface bonding strength caused by the decrease in temperature, resulting in a decrease in compression strength at a lower temperature.

4.2. Dynamic Bending Experimental Results

Figure 9 shows the damage evolution process captured by the high-speed camera of the dynamic bending sample of CFRP laminates at different temperatures and impact velocities, as well as the DIC analysis results at 1   ms , and it can be seen that the maximum bending strain of the dynamic bending sample appears at the impact position of the incident bar. The damage degree of the dynamic bending sample increases obviously with the increase in impact velocity, but the temperature has little influence on the damage degree of the sample, which roughly shows that the damage degree decreases first and then increases with the decrease in temperature. According to the morphology analysis of the samples after impact, the damage degree of the dynamic bending sample varies from complete fracture to incomplete fracture to complete fracture with the decrease in temperature under the 8.5   m / s impact velocity, and then the relationship between the damage degree of the dynamic bending sample and temperature can be proved. In Figure 10, the optical microscope image of the fracture position further explains the influence mechanism of the temperature effect. Delamination is obvious at room temperature, indicating that the interlayer interface strength is lower, and it is easier to completely fracture. However, the interfacial debonding between the fiber and the matrix decreases with the temperature decrease, indicating that the interlayer strength increases. It is found that the fracture is shorter at 180   ° C , and the main reason is that the matrix is more prone to brittle fracture with the temperature decrease, so the damage degree increases relative to that at 80   ° C .
In the dynamic bending experiment, only the samples with an impact velocity of 10   m / s are completely fractured, so we calculated the load–displacement curve applied to this group of samples using the equations mentioned in Ref. [37], as shown in Figure 11. The load increases with the increase in displacement at first, and the load drops sharply after reaching the peak load, indicating that the bending sample is damaged or fractured at this time. When the impact displacement is identical, the load applied on the dynamic bending sample first increases and then decreases with the decrease in temperature, which is caused by the combined action of the interface strengthening and the matrix toughness deterioration caused by the decrease in temperature.

4.3. Numerical Prediction Results

4.3.1. Numerical Prediction Results of Dynamic Compression

Figure 12 shows the FEA results of the dynamic compression samples of CFRP laminates under 1600 / s strain rate compression. As shown in Figure 12a, when the impact process reaches 0.21   ms , a large strain appears in the lower left corner, indicating that the crack may be initiated in one layer; the layers on both sides of the crack are entirely separated at 0.25   ms ; further impact produces a greater crack and new delamination; and buckling of the laminates is observed at 0.40   ms . Similarly, the stress nephogram of Figure 12b shows the above damage evolution law. By comparing the experimental results of the 1200 / s strain rate in Figure 5a, it can be seen that the FEA results of the damage evolution process of the dynamic compression sample are consistent with the experimental results.
Figure 6a presents the FEA and experimental results of the stress–strain curves of the dynamic compression samples at room temperature and different strain rates. It can be seen that the stress–strain curves with 800 / s and 1600 / s strain rates are in good agreement, and the errors of compression strength are within 5 % . The error of the 1200 / s strain rate is large, possibly due to the deviation between the material parameters and some experimental results. Figure 6f presents the FEA and experimental results of the stress–strain curves of the dynamic compression samples at different temperatures and a 1600 / s strain rate. It can be seen that the trend of the curve obtained by FEA is consistent with the experiment, and the error for compression strength is within 10 % .

4.3.2. Numerical Prediction Results of Dynamic Bending

Figure 13 shows the stress nephogram of the FEA results of the dynamic bending sample applied at an impact velocity of 8.5 m / s . When the impact load is applied to the sample, the single-layer plate on the right side of the bending will be subjected to greater tensile stress. Once the damaged condition is reached, the right single-layer plate is damaged, then is gradually damaged by layer until complete fractures occur. The damage process of FEA is consistent with the experimental results of Figure 9.
.
Figure 11 presents the FEA and experimental results of the load–displacement curves of the dynamic bending samples at different temperatures and a 10   m / s impact velocity. It can be seen that the trend of the load–displacement curve changing with temperature is consistent with the experimental results, and the errors for peak load are within 10 % at all temperatures.

5. Conclusions

In this study, the dynamic mechanical properties of CFRP laminates at ultra-low temperatures are systematically studied through experimental testing, a theoretical model, and simulation analysis, and the damage evolution process is discussed. There are the following important conclusions:
  • We independently designed the observable cryogenic box and improved the Hopkinson bending bar. Based on the SHPB device, we set up an ultra-low temperature dynamic experimental platform with a synchronous observation function; the dynamic mechanical properties of CFRP laminates at ultra-low temperatures were tested, and the damage evolution process was observed simultaneously.
  • The experimental results are as follows: CFRP laminates exhibit a noticeable strain rate effect during dynamic compression; the compression strength and modulus increase linearly with the increase in strain rate and show a quadratic function variation trend of first increasing and then decreasing with the increase in temperature. The damage degree of the dynamic bending sample increases obviously with the increase in impact velocity and decreases first and then increases with the decrease in temperature, and the damage degree is minimum at 80   .
  • Based on the ultra-low temperature dynamic constitutive, failure criterion, and interlayer interface damage constitutive of CFRP laminates, a numerical prediction model was established to predict the mechanical properties and damage evolution process of the dynamic compression and bending of CFRP laminates at ultra-low temperatures. The predicted results of the relationship between the dynamic mechanical properties and strain rate and temperature agree with the experimental results. The FEA results of the damage evolution process of CFRP laminates are basically consistent with the experimental observations.

Author Contributions

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

Funding

This work was supported by the National Key Research and Development of China (2018YFA0702804).

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Islam, M.S.; Melendez-Soto, E.; Castellanos, A.G. Prabhakar, Investigation of woven composites as potential cryogenic tank materials. Cryogenics 2015, 72, 82–89. [Google Scholar] [CrossRef] [Green Version]
  2. Soutis, C. Carbon fiber reinforced plastics in aircraft construction. Mater. Sci. Eng. A 2005, 412, 171–176. [Google Scholar] [CrossRef]
  3. Grodzki, W.; Łukaszewicz, A. Design and manufacture of umanned aerial vehicles (UAV) wing structure using composite materials. Mater. Werkst. 2015, 46, 269–278. [Google Scholar] [CrossRef]
  4. Murray, B.R.; Doyle, A.; Feerick, P.J.; Semprimoschnig, C.O.A.; Leen, S.B.; Brádaigh, C.M.Ó. Rotational moulding of PEEK polymer liners with carbon fibre/PEEK over tape-placement for space cryogenic fuel tanks. Mater. Des. 2017, 132, 567–581. [Google Scholar] [CrossRef] [Green Version]
  5. Tomassetti, G.; Barboni, R.; de Benedetti, M. Optimisation methodology for cryotanks. Comput. Struct. 2005, 83, 2293–2305. [Google Scholar] [CrossRef]
  6. Takeda, T.; Narita, F.; Shindo, Y.; Sanada, K. Cryogenic through-thickness tensile characterization of plain woven glass/epoxy composite laminates using cross specimens: Experimental test and finite element analysis. Compos. Part B Eng. 2015, 78, 42–49. [Google Scholar] [CrossRef]
  7. Meng, J.; Wang, Y.; Yang, H.; Wang, P.; Lei, Q.; Shi, H.; Lei, H.; Fang, D. Mechanical properties and internal microdefects evolution of carbon fiber reinforced polymer composites: Cryogenic temperature and thermocycling effects. Compos. Sci. Technol. 2020, 191, 108083. [Google Scholar] [CrossRef]
  8. Kim, M.-G.; Kang, S.-G.; Kim, C.-G.; Kong, C.-W. Tensile response of graphite/epoxy composites at low temperatures. Compos. Struct. 2007, 79, 84–89. [Google Scholar] [CrossRef]
  9. Jia, Z.; Li, T.; Chiang, F.-P.; Wang, L. An experimental investigation of the temperature effect on the mechanics of carbon fiber reinforced polymer composites. Compos. Sci. Technol. 2018, 154, 53–63. [Google Scholar] [CrossRef]
  10. Griffiths, L.J.; Parry, D.J.; Worthington, R.P. A Comparison of Optical and Strain Gauge Techniques in The Determination of the Dynamic Mechanical Behaviour of Carbon-Fibre Composite using a Split Hopkinson Pressure Bar. NASA STI/Recon Tech. Rep. A 1980, 80, 62–70. [Google Scholar]
  11. Guedes, R.M.; de Moura, M.F.S.F.; Ferreira, F.J. Failure analysis of quasi-isotropic CFRP laminates under high strain rate compression loading. Compos. Struct. 2008, 84, 362–368. [Google Scholar] [CrossRef]
  12. Naik, N.K.; Kavala, V.R. High strain rate behavior of woven fabric composites under compressive loading. Mater. Sci. Eng. A 2008, 474, 301–311. [Google Scholar] [CrossRef]
  13. Hosur, M.V.; Alexander, J.; Vaidya, U.K.; Jeelani, S. High strain rate compression response of carbon/epoxy laminate composites. Compos. Struct. 2001, 52, 405–417. [Google Scholar] [CrossRef]
  14. Gomezdelrio, T. Dynamic tensile behaviour at low temperature of CFRP using a split Hopkinson pressure bar. Compos. Sci. Technol. 2005, 65, 61–71. [Google Scholar] [CrossRef] [Green Version]
  15. Tan, H.; Huang, X.; Liu, L.; Guan, Y.; Chen, W. Dynamic compressive behavior of four-step three-dimensional braided composites along three directions. Int. J. Impact Eng. 2019, 134, 103366. [Google Scholar] [CrossRef]
  16. Li, T.; Duan, Y.; Jin, K.; Suo, T.; Yu, X.; Li, Y. Dynamic compressive fracture of C/SiC composites at different temperatures: Microstructure and mechanism. Int. J. Impact Eng. 2017, 109, 391–399. [Google Scholar] [CrossRef]
  17. Gliszczynski, A.; Wiącek, N. Experimental and numerical benchmark study of mode II interlaminar fracture toughness of unidirectional GFRP laminates under shear loading using the end-notched flexure (ENF) test. Compos. Struct. 2021, 258, 113190. [Google Scholar] [CrossRef]
  18. Wysmulski, P. Non-linear analysis of the postbuckling behaviour of eccentrically compressed composite channel-section columns. Compos. Struct. 2023, 305, 116446. [Google Scholar] [CrossRef]
  19. Debski, H.; Samborski, S.; Rozylo, P.; Wysmulski, P. Stability and Load-Carrying Capacity of Thin-Walled FRP Composite Z-Profiles under Eccentric Compression. Materials 2020, 13, 2956. [Google Scholar] [CrossRef]
  20. 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]
  21. Liu, Y.J.; Jiang, Z.; Wen, H.M. Predicting impact induced delamination of FRP laminates. Int. J. Impact Eng. 2020, 137, 103436. [Google Scholar] [CrossRef]
  22. Yan, M.; Jiao, W.; Yang, F.; Ding, G.; Zou, H.; Xu, Z.; Wang, R. Simulation and measurement of cryogenic-interfacial-properties of T700/modified epoxy for composite cryotanks. Mater. Des. 2019, 182, 108050. [Google Scholar] [CrossRef]
  23. Ren, M.-F.; Zhang, X.-W.; Huang, C.; Wang, B.; Li, T. An integrated macro/micro-scale approach for in situ evaluation of matrix cracking in the polymer matrix of cryogenic composite tanks. Compos. Struct. 2019, 216, 201–212. [Google Scholar] [CrossRef]
  24. Huang, C.; Ren, M.; Li, T.; Chang, X.; Cong, J.; Lei, Y. Trans-scale modeling framework for failure analysis of cryogenic composite tanks. Compos. Part B Eng. 2016, 85, 41–49. [Google Scholar] [CrossRef]
  25. Zhou, Y.; Wang, Y.; Xia, Y.; Jeelani, S. Tensile behavior of carbon fiber bundles at different strain rates. Mater. Lett. 2010, 64, 246–248. [Google Scholar] [CrossRef]
  26. Schutz, J.B. Properties of composite materials for cryogenic applications. Cryogenics 1998, 38, 3–12. [Google Scholar] [CrossRef]
  27. GB/T1449-2005; Fibre-Reinforced Plastic Composites-Determination of Flexural Properties. Beijing FRP Research Institute: Beijing, China, 2005.
  28. Atli-Veltin, B. Cryogenic performance of single polymer polypropylene composites. Cryogenics 2018, 90, 86–95. [Google Scholar] [CrossRef]
  29. Wei, R.; Wang, X.; Chen, C.; Zhang, X.; Xu, X.; Du, S. Effect of surface treatment on the interfacial adhesion performance of aluminum foil/CFRP laminates for cryogenic propellant tanks. Mater. Des. 2017, 116, 188–198. [Google Scholar] [CrossRef]
  30. Xie, W.; Zhang, W.; Yang, H.; Jiang, X.; Wang, L.; Gao, Y.; Cai, X. Mechanical behavior and damage kinetics of a unidirectional carbon/epoxy laminated composite under dynamic compressive loading. Polym. Compos. 2022, 43, 2909–2923. [Google Scholar] [CrossRef]
  31. He, C.; Ge, J.; Qi, D.; Gao, J.; Chen, Y.; Liang, J.; Fang, D. A multiscale elasto-plastic damage model for the nonlinear behavior of 3D braided composites. Compos. Sci. Technol. 2019, 171, 21–33. [Google Scholar] [CrossRef]
  32. Jiang, B.; Cao, L.; Zhu, F. Dynamic tensile behavior of polypropylene with temperature effect. Compos. Part B Eng. 2018, 152, 300–304. [Google Scholar] [CrossRef]
  33. Hashin, Z. Cumulative damage theory for composite materials: Residual life and residual strength methods. Compos. Sci. Technol. 1985, 23, 1–19. [Google Scholar] [CrossRef]
  34. 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]
  35. Zhang, J.; Hu, M.; Liu, S.; Wang, L.; Gu, B.; Sun, B. High strain rate compressive behaviors and adiabatic shear band localization of 3-D carbon/epoxy angle-interlock woven composites at different loading directions. Compos. Struct. 2019, 211, 502–521. [Google Scholar] [CrossRef]
  36. Kolsky, H. An Investigation of the Mechanical Properties of Materials at very High Rates of Loading. Proc. Phys. Soc. Sect. B 1949, 62, 676. [Google Scholar] [CrossRef]
  37. Wang, Z.; Suo, T.; Sheikh, M.Z.; Li, Y.; Wang, X.; Wang, Y. Quasi-static and dynamic flexural behavior of annealed and chemically strengthened aluminosilicate glass with notch defects. J. Non-Cryst. Solids 2019, 521, 119479. [Google Scholar] [CrossRef]
Figure 1. Ultra-low temperature dynamic experimental platform with synchronous observation function.
Figure 1. Ultra-low temperature dynamic experimental platform with synchronous observation function.
Materials 16 02654 g001
Figure 2. (a) Dynamic compression sample and (b) dynamic bending sample of CFRP laminates.
Figure 2. (a) Dynamic compression sample and (b) dynamic bending sample of CFRP laminates.
Materials 16 02654 g002
Figure 3. Finite element model. (a) Hopkinson pressure bar and dynamic compression sample. (b) Hopkinson bending bar and dynamic bending sample.
Figure 3. Finite element model. (a) Hopkinson pressure bar and dynamic compression sample. (b) Hopkinson bending bar and dynamic bending sample.
Materials 16 02654 g003
Figure 4. Calculation flowchart of ABAQUS calling user subroutine VUMAT.
Figure 4. Calculation flowchart of ABAQUS calling user subroutine VUMAT.
Materials 16 02654 g004
Figure 5. The dynamic compression samples of CFRP laminates: (a) the damage evolution process captured by high-speed camera and corresponding DIC analysis results during compression at a 1200 / s strain rate; (b) the final damage morphology at different temperatures and strain rates.
Figure 5. The dynamic compression samples of CFRP laminates: (a) the damage evolution process captured by high-speed camera and corresponding DIC analysis results during compression at a 1200 / s strain rate; (b) the final damage morphology at different temperatures and strain rates.
Materials 16 02654 g005
Figure 6. Stress–strain curves of the dynamic compression samples at different strain rates and temperatures: (a) room temperature; (b) 80   ° C ; (c) 180   ° C ; (d) 800 / s strain rate; (e) 1200 / s strain rate; (f) 1600 / s strain rate.
Figure 6. Stress–strain curves of the dynamic compression samples at different strain rates and temperatures: (a) room temperature; (b) 80   ° C ; (c) 180   ° C ; (d) 800 / s strain rate; (e) 1200 / s strain rate; (f) 1600 / s strain rate.
Materials 16 02654 g006
Figure 7. The relationship between mechanical properties of the dynamic compression sample and strain rate and temperature: (a) compression strength–strain rate, (b) compression modulus–strain rate, (c) compression strength–temperature, (d) compression modulus–temperature.
Figure 7. The relationship between mechanical properties of the dynamic compression sample and strain rate and temperature: (a) compression strength–strain rate, (b) compression modulus–strain rate, (c) compression strength–temperature, (d) compression modulus–temperature.
Materials 16 02654 g007
Figure 8. SEM image of the damage location after dynamic compression at different temperatures.
Figure 8. SEM image of the damage location after dynamic compression at different temperatures.
Materials 16 02654 g008
Figure 9. The damage evolution process captured by high-speed camera of the dynamic bending samples of CFRP laminates at different temperatures and impact velocities.
Figure 9. The damage evolution process captured by high-speed camera of the dynamic bending samples of CFRP laminates at different temperatures and impact velocities.
Materials 16 02654 g009
Figure 10. Optical microscope image of the fracture position after dynamic impact at different temperatures.
Figure 10. Optical microscope image of the fracture position after dynamic impact at different temperatures.
Materials 16 02654 g010
Figure 11. Load–displacement curves of the dynamic bending samples at different temperatures and 10   m / s impact velocity.
Figure 11. Load–displacement curves of the dynamic bending samples at different temperatures and 10   m / s impact velocity.
Materials 16 02654 g011
Figure 12. The FEA results of dynamic compression samples of CFRP laminates under 1600⁄s strain rate compression: (a) strain nephogram, (b) stress nephogram.
Figure 12. The FEA results of dynamic compression samples of CFRP laminates under 1600⁄s strain rate compression: (a) strain nephogram, (b) stress nephogram.
Materials 16 02654 g012
Figure 13. The stress nephogram of the FEA results of the dynamic bending sample applied at the impact velocity of 8.5   m / s .
Figure 13. The stress nephogram of the FEA results of the dynamic bending sample applied at the impact velocity of 8.5   m / s .
Materials 16 02654 g013
Table 1. Static material properties of single-layer plate at room temperature [7].
Table 1. Static material properties of single-layer plate at room temperature [7].
Modulus E 1 E 2   =   E 3 G 12   =   G 13 G 23 v 12   =   v 13 v 23
138   GPa 7   GPa 4.8   GPa 3.8   GPa 0.3 0.35
Strength X t X c Y t Y c S 12 = S 23 = S 13
2500   MPa 800   MPa 80   MPa 150   MPa 110   MPa
Fracture toughness W f t W f c W m t W m c
12.5   N / mm 12.5   N / mm 0.1   N / mm 0.1   N / mm
Table 2. Tensile properties of single-layer plate at different temperatures [7].
Table 2. Tensile properties of single-layer plate at different temperatures [7].
SamplesTemperatureModulus G P a Strength M P a
0 ° RT 150 2181
173   K 138 1867
77   K 139 1828
90 ° RT 8.8 56.6
173   K 11.3 55.6
77   K 15.4 53.2
Table 3. The strength of the single-layer plate in different directions and the corresponding material parameters at 0.001 / s reference strain rate.
Table 3. The strength of the single-layer plate in different directions and the corresponding material parameters at 0.001 / s reference strain rate.
σ 0 A n
Fiber direction 800   Mpa 2.04 0.41
Transverse direction 150   Mpa 4.65 0.22
Thickness direction 150   Mpa 0.11 0.50
Table 4. The parameters of the bilinear cohesive region model at room temperature.
Table 4. The parameters of the bilinear cohesive region model at room temperature.
K n   =   K s   =   K t σ n 0 τ s 0 = τ t 0 G n C G s C η B K
1 × 10 6   N / mm 2 45   MPa 60   MPa 1   N / mm 1.5   N / mm 1
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

Zhao, W.; Lin, S.; Wang, W.; Yang, Y.; Yan, X.; Yang, H. Study on Dynamic Mechanical Properties of Carbon Fiber-Reinforced Polymer Laminates at Ultra-Low Temperatures. Materials 2023, 16, 2654. https://doi.org/10.3390/ma16072654

AMA Style

Zhao W, Lin S, Wang W, Yang Y, Yan X, Yang H. Study on Dynamic Mechanical Properties of Carbon Fiber-Reinforced Polymer Laminates at Ultra-Low Temperatures. Materials. 2023; 16(7):2654. https://doi.org/10.3390/ma16072654

Chicago/Turabian Style

Zhao, Wenhao, Sanchun Lin, Wenfeng Wang, Yifan Yang, Xuan Yan, and Heng Yang. 2023. "Study on Dynamic Mechanical Properties of Carbon Fiber-Reinforced Polymer Laminates at Ultra-Low Temperatures" Materials 16, no. 7: 2654. https://doi.org/10.3390/ma16072654

APA Style

Zhao, W., Lin, S., Wang, W., Yang, Y., Yan, X., & Yang, H. (2023). Study on Dynamic Mechanical Properties of Carbon Fiber-Reinforced Polymer Laminates at Ultra-Low Temperatures. Materials, 16(7), 2654. https://doi.org/10.3390/ma16072654

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