Next Article in Journal
Additively Manufactured Continuous Cell-Size Gradient Porous Scaffolds: Pore Characteristics, Mechanical Properties and Biological Responses In Vitro
Previous Article in Journal
Magnetic Nanomaterials as Contrast Agents for MRI
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on Properties Prediction and Braiding Optimization of Axial Braided Carbon/Carbon Composite

1
College of Electrical and Control Engineering, Shaanxi University of Science & Technology, Xi’an 710021, China
2
College of Architecture and Civil Engineering, Beijing University of Technology, Beijing 100124, China
3
The Fourth Academy of CASC, Xi’an 710025, China
4
College of Architecture and Civil Engineering, Xi’an University of Science and Technology, Xi’an 710054, China
5
Science and Technology on Reliability and Environment Engineering Laboratory, Beijing Institute of Structure and Environment Engineering, Beijing 100076, China
*
Author to whom correspondence should be addressed.
Materials 2020, 13(11), 2588; https://doi.org/10.3390/ma13112588
Submission received: 29 February 2020 / Revised: 18 May 2020 / Accepted: 29 May 2020 / Published: 5 June 2020

Abstract

:
It is well established that the microstructure has significant effects on the properties of axial braided C/C composites. In this study, a method coupling the homogenization method and finite element method (FEM) was proposed to predict the relationship between the microstructure characteristics and macroscopic properties. Based on the representative volume element (RVE) model, the periodic displacement boundary condition was introduced to predict the equivalent elastic properties of the RVE and component of C/C composite material, and the coefficient of thermal expansion (CTE) of the material was predicted by the energy prediction method. The predicted results were in good agreement with experimental results. By predicting the thermal and mechanical properties of the materials with different braiding spacing and fiber rod diameter, the variation of the properties with braiding spacing and fiber rod diameter was obtained. The research methods and results in this paper could provide important references for the optimization and rational application of composite materials.

1. Introduction

Axial braided C/C composite is a special carbon fiber reinforced carbon matrix composite. Both the reinforcing phase and matrix phase are composed of pure carbon [1,2,3,4,5]. It exhibits many superior properties, such as high specific strength, good thermal stability, wear resistance and ablative resistance. Specially, different from common materials, its mechanical properties do not decrease but rather increase with the increase of temperature, which makes it widely used in aviation and aerospace fields [6,7]. Due to its complex internal microstructures, such as porosity, cracks, pore mean and pore shape, it is difficult to quantify the effects of microstructure characteristics on its macroscopic properties based on experiments. Therefore, numerical and theoretical methods calibrated and verified by experiments are required to evaluate the relationship between its microstructure characteristics and macro properties [8,9,10,11,12].
Many methods have been introduced to predict the mechanical properties of braided C/C composites [13,14,15,16]. Among them, the stiffness performance has been well investigated. With the development of finite element method (FEM), the method for predicting the properties based on the homogenization theory provides a new idea for predicting the properties of C/C composite materials. Because the FEM is used in this method, the microstructure characteristics and distribution in the fiber rod (bundle) can be considered. Because the homogenization theory is combined, the interaction between components can be taken into account in the model, so that the stiffness prediction accuracy of C/C composite materials is higher [17,18,19]. The energy method-based prediction method of CTE of composite materials derived by Zhang [20,21,22,23] is more convenient to calculate than the homogenization method, and has been gradually popularized in the industry.
In this paper, the braided structure characteristics of axial braided C/C composites were considered, the homogenization method was combined with FEM by ABAQUS. Besides, the energy method was introduced to establish the prediction method of macroscopic properties of materials. A Python program was used to generate the matrix, interface and fiber rod (bundle) geometric models that meet the distribution law of microstructure characteristics and predict their stiffness. The predicted results and the measured results were taken as the input of FEM. Based on the RVE model of materials, periodic displacement boundary conditions were introduced to predict the equivalent elastic properties and the CTE of axial braided C/C composites [17,18,19,20]. By predicting the macroscopic properties of the material with different braiding spacing and fiber rod diameter, the variation of material properties was obtained, which provided a reference for the optimization and rational application of the material.

2. A Method for Predicting the Macroscopic Properties of Axial Braided C/C Composite

2.1. Homogenization Process and Stiffness Prediction Method

The homogenization theory has been widely recognized to predict the stiffness of composites, and its concept is shown in Figure 1. In order to solve the difficulty of large-scale macro-micro joint computation, the concept of microelement is introduced. A microelement represents a point in the composite structure, and its average stress and strain relation is regarded as the effective constitutive relation of composite. In fact, the non-uniform element is replaced by a representative volume element (RVE) with the above average stress and strain relationship. The above concepts are applied to every point of the composite structure. Based on such way, the original problem of analyzing the non-uniform composite structure was converted to the problem of analyzing the properties of the macro structure after homogenization, which greatly simplified the original problem [24,25].
Meso-level composites often exhibit a high degree of heterogeneity, and even in the neighborhood ε with a very small macroscopic position x , the field variables will vary greatly, as shown in Figure 2.
Considering the periodicity of the composite structure, these field variables can be expressed as follows:
Φ ε ( x ) = Φ ( x , y + k Y ) , y = x / ε , ϕ ε x   =   ϕ x + 1 ε ϕ y
where, x is the macroscopic scale, y is the mesoscopic scale, ε is the ratio of the two scales, Y is the period of the periodic function, and k is the integer.
The displacement is expressed as a function of double-scale coordinates and expanded according to small parameters:
u ε ( x ) = ( x , y ) + ε 1 u 1 ( x , y ) + ε 2 u 2 ( x , y ) + ,   y = x / ε
u ( i ) ( x , y ) is a periodic function of Y . According to constitutive equation and geometric equation in elastic mechanics, the above asymptotic expansion was substituted into the equilibrium equation, and the left side of the equation was arranged as an asymptotic expansion with respect to ε . We set the coefficient of ε i ( i = 1 , 0 , 1 , ) as zero to obtain a series of governing equations, in which:
σ i j ( 1 ) ( x , y ) y j   = 0
σ i j ( 0 ) ( x , y ) y j + σ i j ( 1 ) ( x , y ) x j = 0
In Equation (4), σ i j ( x , y ) is a function of stress field. By simplifying the operation, the following equation was obtained:
σ i j ( 0 ) = ( C i j k l C i j m n χ m k l / y n ) ε x k l ( u ( 0 ) )
[ C i j k l ( χ ( k , y l ) m n + I k l m n ) ] y j = 0 on   Y
χ i m n ( y ) = χ i m n ( y   + Y ) on   Y
χ i m n ( y ) = 0 on   Y vert
Among them:
I k l m n = ( δ m k δ n l + δ n k δ m l ) / 2
χ ( k , y l ) m n = 1 2 ( χ k m n y l + χ l m n y k )
The periodic function Y can be expressed as the region of a single cell in the axial braided carbon/carbon composite. Y is the boundary of the cell, and Y vert is the vertex of the cell on the boundary. The homogenization coefficient of a single cell can be expressed as:
C i j k l H = 1 | Y | Y C i j k l ( χ ( k , y l ) m n + I k l m n ) d V
σ ( 0 ) ( x , y ) is the second term of the asymptotic expansion of σ ε ( x , y ) , which is called the mesoscopic stress field. The mesoscopic stress field includes the macroscopic parameters x and the mesoscopic parameters y , which reflect the fluctuation of stress in the scale of RVE.
Gauss theorem is used to treat the governing equation of the homogenization coefficient of Equation (11), and the equation is converted into a “weak” solution equivalent to the original problem as follows:
Y C i j k l χ ( k , y l ) m n   1 2 ( δ v i y j + δ v j y i ) d V = p Y p Y p Y C i j m n n j δ v i d S
From Equation (12) can be found in the original control equation is equivalent to the three dimensional linear elastic problems in the form of the "weak" solution of equation, the difference is that the material within a plane distribution related to the material properties and boundary shape, the plane distribution of force on the internal boundary of different material in the RVE ( Y p Y p Y ), pointing to the internal material area, the direction of the i component for C i j m n n j .
The two terms in the integral sign of Equation (12) were separated to obtain the deformation of the homogenization coefficient:
C i j k l H = 1 | Y | Y C i j k l χ ( k , y l ) m n d V + 1 | Y | Y C i j m n d V
From Equation (13), it can be seen that in the boundary force method, the homogenization coefficient is equal to the volume average of the stress of the three-dimensional linear elastic problem on the RVE plus the volume average of the elastic tensor on the RVE [16,17,18]. The homogenization coefficient is the effective elastic property of macroscopic structure.

2.2. A Method for Predicting CTE of Composites Based on Energy Method

The CTE prediction method based on the energy method is more convenient to calculate than the simple homogenization method. Therefore, the following method combined with FEM was used to calculate the CTE of materials in this paper [19,22].
The basic idea of the energy method is that based on the relationship between microstructure and homogeneous equivalent body, the energy expression of the equivalent property of composite materials can be obtained by deducing the relationship between the equivalent property of composite materials and the deformation energy of microstructure [16,17].
For the microstructure shown in Figure 1, when the temperature rise value is Δ T , the average stress and average strain of the microstructure have the following relationship:
σ ¯ = C H ( ε ¯ α H Δ T )
where C H , α H = [ α 11 H , α 22 H , α 33 H ] are the stiffness matrix and CTE matrix of axial braided C/C composites.
The deformation energy of the microstructure of the material is shown in Equation (15):
C = Ω ε T D ε d Ω
For working condition 1 ( u 2 = u 3 = 0 ) as shown in Figure 3a, the microstructure deformation energy caused by Δ T temperature change is the same as that of working condition 4 without temperature load as shown in Figure 3b. The displacement load u 1 in working condition 4 is equal to the displacement load under working condition 1. The force f is the interface force caused by the difference of thermal stress between the two loads. The expressions of the two are Equations (16) and (17).
u 1 = ( α 11 H Δ T + ( D 1122 H / D 1111 H ) α 22 H Δ T + ( D 1133 H / D 1111 H ) α 33 H ) l
f = s [ D ( E 1 ) α ( α 1 ) D ( E 2 ) α ( α 2 ) ] Δ T d s
For the microstructure of working condition 1, the variable properties of its homogeneous equivalent are equal to the difference between the variable properties of working condition 1 and that of working condition 6, as shown in Equation (18):
C e m ( 1 ) = C m ( 1 ) C m ( 6 ) = D 1111 H ( α 11 H Δ T + ( D 1122 H / D 1111 H ) α 22 H Δ T + ( D 1133 H / D 1111 H ) α 33 H Δ T V ) 2
Similar to Equation (18), the variable properties of the homogeneous equivalent of condition 2 ( u 1 = u 3 = 0 ) and condition 3 ( u 1 = u 2 = 0 ) can be expressed as Equations (19) and (20):
C e m ( 2 ) = C m ( 2 ) C m ( 6 ) = D 2222 H ( α 22 H Δ T + ( D 1122 H / D 2222 H ) α 11 H Δ T + ( D 2233 H / D 2222 H ) α 33 H Δ T ) 2 C e m ( 3 ) = C m ( 3 ) C m ( 6 ) = D 3333 H ( α 33 H Δ T + ( D 1133 H / D 3333 H ) α 11 H Δ T + ( D 2233 H / D 3333 H ) α 22 H Δ T ) 2
The FEM equation of the microstructure under operating conditions 4, 5 and 6 can be expressed as Equation (20).
[ k 11 k 12 k 21 k 22 ] [ U 1 ( i ) U 2 ( i ) ] = [ F 1 ( i ) F 2 ( i ) ]   i = 4 , 5 , 6
where, U 1 ( i ) is the unknown node displacement vector under condition i , and U 2 ( i ) is the known node displacement vector. F 1 ( i ) is the known node force vector under condition i , and F 2 ( i ) is the unknown node force vector. To solve Equation (20), the following Equation (21) can be obtained.
U 1 ( i ) = K 11 1 F 1 ( i ) K 11 1 K 12 U 2 ( i ) F 2 ( i ) = K 21 K 11 1 F 1 ( i ) + ( K 22 K 21 K 11 1 K 12 ) U 2 i
The microstructural variable performance of working condition i can be expressed in the following form:
C m ( i ) = Ω ε ( i ) T D ε ( i ) d Ω = F 1 ( i ) T U 1 ( i ) + F 2 ( i ) T U 2 ( i )
The calculation method of CTE of microstructure can be obtained by solving Equation (22) as shown in Equation (23):
α 11 H = ( D 1122 H D 1111 H D 2233 H D 2222 H D 1133 H D 1111 H ) ( D 2233 H D 3333 H Q R ) ( D 2233 H D 3333 H D 2233 H D 2222 H 1 ) ( D 1122 H D 1111 H Q P ) ( D 1122 H D 1111 H D 2233 H D 2222 H D 1133 H D 1111 H ) ( D 2233 H D 3333 H D 1122 H D 2222 H D 1133 H D 3333 H ) ( D 2233 H D 3333 H D 2233 H D 2222 H 1 ) ( D 1122 H D 1111 H D 1122 H D 2222 H 1 ) α 22 H = ( D 1122 H D 2222 H D 1133 H D 1111 H D 2233 H D 2222 H ) ( D 1133 H D 3333 H P R ) ( D 1133 H D 3333 H D 1133 H D 1111 H 1 ) ( D 1122 H D 2222 H P Q ) ( D 1122 H D 2222 H D 1133 H D 1111 H D 2233 H D 2222 H ) ( D 1133 H D 3333 H D 1122 H D 1111 H D 2233 H D 3333 H ) ( D 1133 H D 3333 H D 1133 H D 1111 H 1 ) ( D 1122 H D 2222 H D 1122 H D 1111 H 1 ) α 33 H = ( D 1133 H D 3333 H D 1122 H D 1111 H D 2233 H D 3333 H ) ( D 2233 H D 3333 H P Q ) ( D 1122 H D 2222 H D 1122 H D 1111 H ) 1 ) ( D 1133 H D 3333 H P R ) ( D 1133 H D 3333 H D 1122 H D 1111 H D 2233 H D 3333 H ) ( D 1122 H D 2222 H D 1133 H D 1111 H D 2233 H D 2222 H ) ( D 1122 H D 2222 H D 1122 H D 1111 H ) 1 ) ( D 1133 H D 3333 H D 1133 H D 1111 H 1 )
In the equation P = C m ( 1 ) C m ( 6 ) D 1111 H ( Δ T ) 2 V , Q = C m ( 2 ) C m ( 6 ) D 2222 H ( Δ T ) 2 V , R = C m ( 3 ) C m ( 6 ) D 3333 H ( Δ T ) 2 V .

2.3. The Boundary Conditions

According to the braided characteristics of the axial braided C/C composites studied in this paper, each surface of the RVE model must be set as a periodic symmetric boundary condition. The displacement of the RVE should also meet the periodic boundary conditions [17]. For the periodic RVE structure with spatial translational symmetry, the periodic condition of its displacement can be expressed as follows:
u 2 ( y 1 0 , y 2 , y 3 ) = u 1 ( y 1 0 + Y 1 , y 2 , y 3 ) u 2 ( y 1 , y 2 0 , y 3 ) = u 1 ( y 1 , y 2 0 + Y 2 , y 3 ) u 2 ( y 1 , y 2 , y 3 0 ) = u 1 ( y 1 , y 2 , y 3 0 + Y 3 )

3. The Introduction of the Experiment Carried out in the Early Stage

3.1. Properties Test of Fiber Rod

Under the premise of minimum damage, the z-direction fiber rod (XY direction is transverse) was stripped, and the tensile and compression properties of the fiber rod were tested to obtain its stiffness properties. A special fixture was designed for its test. In the test, the unconstrained connection mechanism was adopted to eliminate the error caused by eccentricity when the sample was loaded. The accurate strain of fiber rods was obtained by non-contact optical strain testing system (ARAMIS, 2.3M, GOM, Brunswick, Germany). The obtained stiffness properties of the fiber rod are shown in Table 1.

3.2. Macroscopic Properties Test

The macroscopic mechanical properties of axial braided C/C composites mainly include uniaxial tensile test and uniaxial compression test. Samples under different loading conditions were designed according to the braiding mode and product size of the axial braided C/C composite. Samples with different directions including fiber rod direction (axial direction) and fiber bundle direction (radial direction) were tested and more than 10 samples were prepared for each group. The exact strain of specimens was obtained by non-contact optical strain testing system (ARAMIS). The stiffness properties of axial braided C/C composites are shown in Table 2.

4. RVE Establishment and Properties of Component

4.1. RVE Establishment

The RVE of axial braided C/C composite material is a rectangular structure composed of fiber rod, fiber bundle, matrix and interface. The braided structure and phase scale of each component are shown in Figure 4. Thin carbon fiber rigid rod formed by pultrusion forms the axial reinforcing network, and the soft carbon fiber bundle is woven into the pre-texture. The fiber rods are arranged in an equilateral triangle in the axial direction, and the fiber bundles are successively increased through the 0°, 60° and 120° channels of cambium formed by the fiber rods, and so on until the required size of pre-texture is formed. High density 4D C/C composites were prepared by asphalt impregnation, carbonization, densification and high temperature treatment. The minimum unit of this pre-texture is axially symmetric, and its braided thickness accumulates axially, so it is called axial braided C/C composite.
The generated finite element model of axial braided C/C composite material composed from the RVE model is shown in Figure 5. The division of the grid is highly related to the application of periodic boundary conditions, so it is necessary to strictly guarantee the uniformity of the grid nodes on the corresponding surface of RVE. In this paper, C3D8 elements and C3D6 elements were used to mesh the fiber rod, fiber bundle and matrix. The fiber rod/matrix interface and fiber bundle/matrix interface were all meshed by C3D8 elements [24,25,26,27,28].

4.2. Study on the Properties of Components

The premise of RVE model performance calculation with finite element technique is that the properties of each meso components (including the fiber rod and bundle, matrix and interface) are known. Due to the poor stripping and testability of the meso components, some of the necessary properties need to be obtained in combination with numerical calculations.

4.2.1. Mechanical Properties of Matrix and Interface

The pores and cracks in the matrix and interface of the axial braided C/C composite were regarded as inclusions with different shapes and sizes and introduced into the pure carbon matrix. In this paper, the random geometric models meeting objective distribution law of microstructure characteristics were generated by a subroutine coded by python and its mechanical properties were calculated by FEM, which can consider the interaction between various microstructures [28,29,30,31,32].
In order not only to simplify the calculation workload but also to contain sufficient pore information, the side length of the calculation model of the matrix was selected as 0.4 mm, and the side length of the interface was selected as 0.02 mm, and the generated model was shown in Figure 6. In the model, the pores are mainly ellipsoid, and the ratio of long axis and short axis of ellipsoid at some locations is large, resulting in cracks. In the early stage, the non-pore matrix modulus was 12.34 GPa and Poisson’s ratio was 0.23 by using nano-indentation experiment [16,17], which was used as the material input for finite element calculation.
Figure 7a shows the effect of porosity on the elastic properties of the matrix. It can be seen that the elastic properties of the matrix vary linearly with porosity. The increase of porosity leads to the decrease of both elastic modulus and Poisson’s ratio.
Figure 7b shows the influence of porosity on the interface elastic properties. It can be found that the influence trend of porosity on the interface elastic properties is the same as that of porosity on the matrix.
Based on the statistical study in [2], the porosity of the matrix was 0.036, and the modulus and Poisson’s ratio were 11.5 GPa/0.228 respectively, obtained from Figure 7a. Based on the statistical study in [2], the porosity of the interface was 0.27, and the modulus and Poisson’s ratio were 3.91 GPa/0.12, respectively, obtained from Figure 7b. The CTE of the matrix was obtained from the preliminary test results, −1.3 × 10−6/K.

4.2.2. Study on Mechanical Properties of Fiber Rod and Fiber Bundle

In this section, based on the microstructure characteristics and random distribution of fiber rod (bundle), FEM and Python language programming were used to generate the fiber rod (bundle) geometric model satisfying the distribution law of microstructure characteristics, and the homogenization method was used to predict the stiffness properties. Figure 8 shows the microscopic observation image and geometric model of the fiber rod. The model consists of fiber monofilament, matrix and pores.
Based on relevant experiments and combined with literature [16,17], the input parameters of fiber monofilament and matrix were determined as shown in Table 3. According to the literatures, these data were directly measured by the experiment.
The stiffness properties of the reinforcing phase were calculated by FEM. The results show that the reinforcing phase was transversally isotropic even if there are a few pores. Table 4 lists the calculation results of the enhancement phase and some experimental results in the early stage. It can be seen from the table that the prediction error of fiber rod’s tensile modulus, compression modulus and shear modulus is 13.4%, 13.7% and 10.4% respectively, which verifies the reliability of the calculation method for the reinforcement phase. There are two reasons for the error: the error of the test system and the test method; The inaccuracy of porosity observation leads to the inaccuracy of FEM. In the subsequent calculations of the macroscopic properties of the material, the longitudinal properties of the rod were given according to the experimental values, while the shear properties and other untested parameters were given according to the predictive values.

5. Calculation Results and Analysis

5.1. The Calculation Results of Elastic Properties

Periodic boundary conditions were applied to the RVE model of the axial C/C composite material to obtain its effective properties as shown in Table 5. By comparing the prediction results with the experimental results under different conditions, it can be found that the elastic modulus of the numerical prediction results is basically consistent with the experimental results, thus verifying the effectiveness of the prediction method.

5.2. The Calculation Result of CTE

Based on the material parameters listed in Table 3 and Table 4, the energy analysis method was used to predict the CTE of the RVE model. According to previous research experience, the effect of porosity on CTE is very weak, so this paper did not study the effect of porosity on CTE. The comparison between the predicted results and the experimental results is shown in Table 6. The experimental data of thermal expansion in the table were tested by DIL 402 PC (NETZSCH, Lanzhou, China), and the sample size was φ 6 mm × 50 mm.
The comparison between the results of prediction and experimental results shows that the prediction method based on energy method has certain prediction accuracy. However, the predicted value in both directions is larger than the experimental value, which may be caused by the fact that the CTE of the matrix is larger than the real value.

5.3. Effect of Braiding Spacing on Stiffness and CTE

The stiffness and CTE of the material can be predicted by changing the center spacing of the fiber rod. The influence of braiding spacing on the tensile modulus and CTE is shown in Figure 9. With the increase of braiding spacing, the tensile modulus and CTE in axial and radial directions decrease gradually. The main reason is that the volume content of axial fiber rod and radial fiber bundle in RVE decreases while the volume fraction of matrix increases with the increase of braiding distance. Since the modulus and CTE of the matrix are lower than that of the fiber rod and fiber bundle, the increase of the volume fraction will inevitably lead to the overall decrease of the modulus and CTE [28,29,30,31,32,33].
From the perspective of the use environment of C/C composite materials, it is necessary to reduce the modulus and CTE to reduce the thermal stress of the structural parts. Therefore, the increase of braiding spacing of materials is more in line with the requirements. However, studies have shown that the fiber volume fraction is closely related to the ablative resistance of the material. When the fiber volume fraction decreases, the ablative rate of the material increases, which is unfavorable for engineering application. Therefore, the optimization of braiding spacing of C/C composites should be combined with the ablative properties of the materials. The influence of braiding spacing on ablation rate will be studied in the future.

5.4. Effect of Fiber Rod Diameter on Stiffness and CTE

The influence of fiber rod diameter on tensile modulus and the CTE is shown in Figure 10a. As can be seen, with the increase of fiber rod diameter, the axial tensile modulus increases significantly, while the CTE decreases significantly. This is mainly because, under the condition of constant braiding spacing, the diameter of fiber rod increases, which significantly increases the content of axial fiber in RVE. In addition, due to the higher modulus of axial fiber and the lower CTE, the axial modulus of the material increases and the CTE decreases.
Radial tensile modulus and CTE are both incremental functions of fiber rod diameter, but the increment is small, as shown in Figure 10b.The reason is that for the radial tension modulus, when the braiding spacing is constant, increasing the diameter of the fiber rod will not cause the change of the volume fraction of the radial fiber bundle without changing the fiber bundle’s own size, but will cause the decrease of the volume fraction of the matrix. Because the modulus of the matrix is smaller than the transverse modulus of the fiber rod, the radial tensile modulus increases. For the CTE, since the CTE of the matrix is smaller than transverse CTE of the fiber rod, when the volume fraction of the fiber rod increases and the volume fraction of the matrix decreases, the increase of transverse CTE of the fiber rod to the radial CTE of the material is greater than the decrease of radial CTE of the matrix. Figure 10b shows that it is appropriate for the diameter of the axial braided C/C composite fiber rod to be around 1.26 mm.

6. Conclusions

In this paper, the stiffness prediction method based on homogenization method and the CTE prediction method based on energy method were established. In the prediction of component stiffness properties, Python language was used to generate a component material model that meets the characteristics of pore microstructure randomly, so that the prediction of component material properties was more accurate. The following conclusions are obtained:
(1)
The FEM was utilized for the calculation of the stiffness properties of the matrix and interface, and the results showed that the porosity was important to determine the elastic properties of the matrix and interface. It can be seen that the elastic properties of the matrix and interface change linearly with the porosity, which increases, leading to decreases of both the elastic modulus and Poisson’s ratio.
(2)
The stiffness performance of the reinforcing phase was calculated by FEM. The calculation results showed that the reinforcing phase manifested transversely isotropic characteristic even if there was a few of pores. It was found that the prediction precision of tensile modulus, compression modulus and shear modulus was reliable, being 13.4%, 13.7% and 10.4%, respectively.
(3)
Based on the homogenization theory and energy method combined with the RVE of axial braided C/C composite material, the macroscopic stiffness and CTE of the material were also predicted. By comparing the predicted results with the experimental results in different states, it can be concluded that the numerical predicted results are basically consistent with the experimental results, thus proving the effectiveness of the prediction method.
(4)
By calculating the macroscopic stiffness performance and CTE under the conditions with different braiding spacing and fiber rod diameter, the relationship between the microstructure and macroscopic properties was evaluated. The results show that the best diameter of the axial braided C/C composite fiber rod studied in this paper is 1.26 mm, and the optimization of braiding spacing needs to be carried out in combination with the ablative property.
(5)
The advantages of the prediction method in this paper are as follows: by combining the homogenization method with FEM, the pore characteristics of the fiber rod (bundle), matrix and interface can be considered, so as to obtain the properties of mesoscopic materials accurately. As the input parameters, the stiffness prediction of the material has higher precision. The disadvantage is that the pore characteristics cannot be added into the finite element model of RVE for the joint calculation of macroscopy-mesoscopic and microscopic, which is also the focus of the follow-up research of this project.

Author Contributions

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

Funding

The National Natural Science Foundation of China (Grant No. 51769028), The Natural Science Foundation of Qinghai Province in China (Grant No. 2017-ZJ-933Q), Beijing Institute of Structure and Environment Engineering Fund (Grant No. BQ2019001).

Acknowledgments

Thanks for the financial support from Shaanxi University of Science & Technology.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hui, W.-H.; Bao, F.-T.; Wei, X.-G.; Liu, Y. Ablation performance of a 4D-braided C/C composite in a parameter-variable channel of a Laval nozzle in a solid rocket motor. New Carbon Mater. 2017, 32, 365–373. [Google Scholar] [CrossRef]
  2. Wang, C.; Tang, M.; Liu, W.; Zhu, T. Study on microstructure characteristics of axially braided carbon/carbon composites based on SEM and Micro-CT. Materials 2020, 13, 1414. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Gao, B.; Tang, M.; Yang, Y.; Shi, H. Mechanical experiment of 4D in-plain C/C composites. Acta Mater. Compos. Sin. 2011, 28, 245–250. [Google Scholar]
  4. Rajneesh, S.; Ramesh, B.A.; Puneet, M. Finite element analysis for mechanical characterization of 4D inplane carbon/carbon composite with imperfect microstructure. Lat. Am. J. Solids Struct. 2014, 11, 170–184. [Google Scholar]
  5. Li, D.S.; Yao, Q.Q.; Jiang, N.; Jiang, L. Bend properties and failure mechanism of a carbon/carbon composite with a 3D needle-punched preform at room and high temperatures. New Carbon Mater. 2016, 31, 437–444. [Google Scholar] [CrossRef]
  6. Ren, J.; Zhang, Y.; Fu, Y.; Zhang, P.; Tian, S.; Zhang, L. Effects of the second phase on the microstructure and ablation resistance of HfC coating on C/C composites. Surf. Coat. Technol. 2018, 344, 250–258. [Google Scholar] [CrossRef]
  7. Shan, Z.; Kang, H.; Zang, Y.; Liu, F. Microscopic structure model and finite element simulation of mechanical properties of fiber-bar composites reinforced by three dimensional weaving. Acta Mater. Compos. Sin. 2015, 32, 138–149. [Google Scholar]
  8. 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]
  9. Castaneda, N.; Wisner, B.; Cuadra, J.; Amini, S.; Kontsos, A. Investigation of the Z-binder role in progressive damage of 3D woven composites. Compos. Part A Appl. Sci. Manuf. 2017, 98, 76–89. [Google Scholar] [CrossRef]
  10. Gerlach, R.; Siviour, R.C.; Wiegand, J.; Petrinic, N. In-plane and through-thickness properties, failure modes, damage and delamination in 3D woven carbon fibre composites subjected to impact loading. Compos. Sci. Technol. 2011, 72, 397–411. [Google Scholar] [CrossRef]
  11. Lu, Z.; Xia, B.; Yang, Z. Investigation on the tensile properties of three-dimensional full five-directional braided composites. Comput. Mater. Sci. USA 2013, 77, 445–455. [Google Scholar] [CrossRef]
  12. Zhang, P. Microstructure Modeling and Mechanical Properties Prediction of 3d Needled C/C Composites; Northwestern Polytechnical University: Xi’an, China, 2016. [Google Scholar]
  13. Yan, S.; Zhang, J.; Sun, B.; Guo, B. In situ measurement of strains at different locations in 3-D braided composites with FBG sensors. Compos. Struct. 2019, 230, 111527. [Google Scholar] [CrossRef]
  14. Dang, H.; Zhao, Z.; Liu, P.; Zhang, C.; Tong, L.; Li, Y. A new analytical method for progressive failure analysis of two-dimensional triaxially braided composites. Compos. Sci. Technol. 2020, 186, 107936. [Google Scholar] [CrossRef]
  15. Hao, W.; Huang, Z.; Zhang, L.; Zhao, G.; Luo, Y. Study on the torsion behavior of 3-D braided composite shafts. Compos. Struct. 2019, 229, 111384. [Google Scholar] [CrossRef]
  16. Xu, C.-H.; Xu, D.-S.; Song, L.-Y.; Xu, K. Microstructure characterization from X-ray micro-tomography and tensile failure mechanism of 4D in-plane carbon/carbon composites. J. Solid Rocket Technol. 2013, 36, 811–815. [Google Scholar]
  17. Tang, M.; Gao, B.; Yang, Y.-C.; Shi, H.-B. Mechanical properties prediction of 4D-in plane carbon/carbon composites based on homogenization method. J. Solid Rocket Technol. 2011, 34, 109–112. [Google Scholar]
  18. Zhu, Z.; Qiang, H. Forecast research on micromechanical properties of 4D in-plane C/C composites basedon homogenization theory. J. Solid Rocket Technol. 2019, 42, 92–96. [Google Scholar]
  19. Shokrieh, M.M.; Mazloomi, M.S. A new analytical model for calculation of stiffness of three-dimensional four-directional braided composites. Compos. Struct. 2012, 94, 1005–1015. [Google Scholar] [CrossRef]
  20. Zhang, W.; Wang, F.; Dai, G. Topology optimal design of material microstructures using strain energy-based method. Chin. J. Aeronaut. 2007, 20, 321–328. [Google Scholar] [CrossRef] [Green Version]
  21. Yevtushenko, A.A.; Grzes, P.; Adamowicz, A. The temperature mode of the carbon-carbon multi-disc brake in the view of the interrelations of its operating characteristics. Materials 2020, 13, 1878. [Google Scholar] [CrossRef]
  22. Vignoles, G.L.; Aspa, Y.; Quintard, M. Modelling of carbon-carbon composite ablation in rocket nozzles. Compos. Sci. Technol. 2010, 70, 1303–1311. [Google Scholar] [CrossRef] [Green Version]
  23. Liu, F.; Guan, Z.; Bian, T. Damage model for predicting shear strength of carbon/carbon composite fastener based on post-failure behavior. Compos. Struct. 2019, 221, 110864.1–110864.17. [Google Scholar] [CrossRef]
  24. Mohammed, A.S.K.; Sehitoglu, H.; Rateick, R. Interface graphitization of carbon-carbon composites by nanoindentation. Carbon 2019, 150, 425–435. [Google Scholar] [CrossRef]
  25. Pimazzoni, A.; Cescon, E.; Dalla, P.M. Thermo-mechanical analysis of unidirectional carbon-carbon composite for thermal imaging diagnostic of a particle beam. Fusion Eng. Des. 2019, 146, 2457–2461. [Google Scholar] [CrossRef]
  26. Li, L.; Aliabadi, M.H. Elastic property prediction and damage mechanics analysis of 3D braided composite. Theor. Appl. Fract. Mech. 2019, 104, 102338. [Google Scholar] [CrossRef]
  27. Li, W.; Fang, G.; Li, W.; Liang, J.; Li, M. Role of mesoscopic features on thermochemical ablative behavior of 3D C/C braided composites. Int. J. Heat Mass Transf. 2019, 144, 118602. [Google Scholar] [CrossRef]
  28. Lachaud, J.; Aspa, Y.; Vignoles, L.G. Analytical modeling of the steady state ablation of a 3D C/C composite. Int. J. Heat Mass Transf. 2008, 51, 2614–2627. [Google Scholar] [CrossRef] [Green Version]
  29. Wang, C.G.; Xu, G.; Gong, J.L. Buckling instability of flexible joint under high pressure in solid rocket motor. Int. J. Aerosp. Eng. 2020, 2020, 10. [Google Scholar]
  30. Chinnaraj, R.K.; Hong, S.M.; Kim, H.S.; Oh, P.Y. Ablation experiments of ultra-high-temperature ceramic coating on carbon-carbon composite using ICP plasma wind tunnel. Int. J. Aeronaut. Space Sci. USA 2020, in press. [Google Scholar] [CrossRef]
  31. Wang, C.G.; Xu, G. The Analytical method to compute the strain on the soft PSD in double-pulse SRM. Int. J. Aerosp. Eng. 2019, 2019, 7. [Google Scholar] [CrossRef] [Green Version]
  32. Wang, C.; Chen, C.; Ren, Q.; Chen, K. Investigation on metal diaphragm of pulse separation device in dual-pulse solid rocket motor. Proc. Inst. Mech. Eng. 2017, 231, 444–451. [Google Scholar]
  33. Chowdhury, P.; Sehitoglu, H.; Rateick, R. Damage tolerance of carbon-carbon composites in aerospace application. Carbon 2018, 126, 382–393. [Google Scholar] [CrossRef]
Figure 1. Composite structure, representative volume element and equivalent structure. (a) original structure (b) micro element (c) homogenized element (d) homogenized structure.
Figure 1. Composite structure, representative volume element and equivalent structure. (a) original structure (b) micro element (c) homogenized element (d) homogenized structure.
Materials 13 02588 g001
Figure 2. Composite material with mesoscopic structure. (a) Macro structure (b) Mesoscopic structure and its periodicity (c) RVE.
Figure 2. Composite material with mesoscopic structure. (a) Macro structure (b) Mesoscopic structure and its periodicity (c) RVE.
Materials 13 02588 g002
Figure 3. Microstructure forms under different operating conditions. Among of them, working condition 2 and working condition 3 are similar to working condition 1, which are not completely shown in the figure.
Figure 3. Microstructure forms under different operating conditions. Among of them, working condition 2 and working condition 3 are similar to working condition 1, which are not completely shown in the figure.
Materials 13 02588 g003
Figure 4. RVE of axial braided C/C composite.
Figure 4. RVE of axial braided C/C composite.
Materials 13 02588 g004
Figure 5. Finite element model of axial braided C/C composite material.
Figure 5. Finite element model of axial braided C/C composite material.
Materials 13 02588 g005
Figure 6. Computational model of matrix and interface.
Figure 6. Computational model of matrix and interface.
Materials 13 02588 g006
Figure 7. Effect of porosity on the elastic properties of matrix and interface.
Figure 7. Effect of porosity on the elastic properties of matrix and interface.
Materials 13 02588 g007
Figure 8. SEM morphology and geometric model of fiber rod (bundle).
Figure 8. SEM morphology and geometric model of fiber rod (bundle).
Materials 13 02588 g008
Figure 9. Effect of center spacing of fiber rod on axial and radial tensile modulus and CTE.
Figure 9. Effect of center spacing of fiber rod on axial and radial tensile modulus and CTE.
Materials 13 02588 g009
Figure 10. Effect of fiber rod diameter on radial tensile modulus and CTE of axial braided C/C composites.
Figure 10. Effect of fiber rod diameter on radial tensile modulus and CTE of axial braided C/C composites.
Materials 13 02588 g010
Table 1. Stiffness characteristics of fiber rod.
Table 1. Stiffness characteristics of fiber rod.
ItemEx, GPa ν x y = ν x z
Mean ValueStandard DeviationMean ValueStandard Deviation
Tension211.110.30.060.02
Compression99.48.470.150.03
Note: Angle mark x represents the fiber rod axial direction, Angle mark y and z represent the transverse direction.
Table 2. Stiffness characteristics of axial braided C/C composites.
Table 2. Stiffness characteristics of axial braided C/C composites.
Item E X , G P a ( G X Y , GPa ) ν X Y E Z , GPa ( G X Z , GPa ) ν X Z
Mean ValueStandard DeviationMean ValueStandard DeviationMean ValueStandard DeviationMean ValueStandard Deviation
Tension38.32.30.280.4140.54.20.090.05
Compression20.91.90.280.0422.12.10.090.05
Shear12.61.3 3.40.5
Note: Angle mark Z represents the fiber rod axial direction, Angle mark X and Y represent the transverse direction.
Table 3. Properties of fiber and matrix in reinforcing phase.
Table 3. Properties of fiber and matrix in reinforcing phase.
Item E x , GPa E y , GPa G x y , GPa ν x y ν y z αx, 10−6/Kαy = αz, 10−6/K
T300 fiber345.1☆120.0★19.026.70.20.35−0.788.14
Matrix12.34--0.23-−1.3−1.3
Note: ☆ represents the tensile modulus. ★ represents the compression modulus. Angle mark x represents the fiber filaments axial direction, Angle mark y and z represent the transverse direction. The parameters of the matrix are the properties without pores.
Table 4. Stiffness properties of the enhanced phase.
Table 4. Stiffness properties of the enhanced phase.
Item E x , GPa E y , GPa G x y , GPa ν x y ν y z αx, 10−6/Kαy = αz, 10−6/K
The predictive value of the fiber rod239.5☆85.8★17.3720.250.110.45−0.945.1
Test value of fiber rod211.1☆99.4★ 22.60.06---
The predictive value of fiber bundle220.8☆79.1★17.1118.960.120.40−0.994.52
Note: ☆ represents the tensile modulus. ★ represents the compression modulus. Angle mark x represents the fiber rod (bundle) axial direction, Angle mark y and z represent the transverse direction.
Table 5. Mechanical properties of axial braided C/C composites.
Table 5. Mechanical properties of axial braided C/C composites.
ItemElastic Modulus GPaPoisson’s RatioShear Modulus GPa
E X E Y E Z ν X Y ν X Z ν Y Z G X Y G X Z G Y Z
Predicted value (tension)37.033.644.60.370.090.1112.23.32.5
Experimental value (tension)38.3-40.50.280.09-12.63.4-
Predicted value (compression)22.320.521.10.380.120.15---
Experimental value (compression)20.9-22.10.280.09----
Note: Angle mark Z represents the fiber rod axial direction, Angle mark X and Y represent the transverse direction. Notice the difference between the direction sign here and the direction sign of the meso component.
Table 6. Coefficient of thermal expansion at room temperature.
Table 6. Coefficient of thermal expansion at room temperature.
ItemAxial Coefficient of Thermal Expansion, 10−6/KRadial Coefficient of Thermal Expansion, 10−6/K
Predicted value0.147−0.477
Experimental value0.145−0.50

Share and Cite

MDPI and ACS Style

Wang, C.; Cao, P.; Tang, M.; Tian, W.; Liu, K.; Liu, B. Study on Properties Prediction and Braiding Optimization of Axial Braided Carbon/Carbon Composite. Materials 2020, 13, 2588. https://doi.org/10.3390/ma13112588

AMA Style

Wang C, Cao P, Tang M, Tian W, Liu K, Liu B. Study on Properties Prediction and Braiding Optimization of Axial Braided Carbon/Carbon Composite. Materials. 2020; 13(11):2588. https://doi.org/10.3390/ma13112588

Chicago/Turabian Style

Wang, Chunguang, Peng Cao, Min Tang, Weiping Tian, Ketong Liu, and Baorui Liu. 2020. "Study on Properties Prediction and Braiding Optimization of Axial Braided Carbon/Carbon Composite" Materials 13, no. 11: 2588. https://doi.org/10.3390/ma13112588

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