Abstract
3D braided composite joints are the important components in CFRP truss, which have significant influence on the reliability and lightweight of structures. To investigate the mechanical performance of 3D braided composite joints, a numerical method based on the microscopic mechanics is put forward, the modeling technologies, including the material constants selection, element type, grid size, and the boundary conditions, are discussed in detail. Secondly, a method for determination of ultimate bearing capacity is established, which can consider the strength failure. Finally, the effect of load parameters, geometric parameters, and process parameters on the ultimate bearing capacity of joints is analyzed by the global sensitivity analysis method. The results show that the main pipe diameter thickness ratio , the main pipe diameter , and the braided angle are sensitive to the ultimate bearing capacity .
1. Introduction
CFRP truss is composed by members and joints, and the composite joint is the most important component, where the force and deformation are very complex, moreover there exists phenomenon of stress concentration, which is the weak point of the loading process [1, 2]. According to statistics, there are 70% spacecraft structure damages that occurred at the connecting part [3, 4]. Therefore, in the analysis and optimization of composite structures, the key point is the prediction of the joint strength and the influence on each parameter.
According to forming process it can be divided into Molding joint, Winding joint, Layer-Molding joint, and 3D braided joint. The first three joints belong to the laminated structure, and the theoretical research and manufacturing process are relatively mature, and they are commonly used in the CFRP truss joints currently [5, 6]. But the laminated CFRP joint strength between layers is small, delamination defects occurs easily. While the 3D braided CFRP also has fiber in thickness direction, overcoming the weakness of delamination in the former three joints [7]. Therefore, the 3-D braided joints are important forms to study on in all kinds of CFRP truss joints.
Recently, researches on 3D braided joint are mainly concentrated on material constants and numerical simulation and experiment. Zheng et al., and so forth, [8] taking CFRP 3-D braided spherical joint as the research object, studied its damage mode under complex loads. C.-Y. Yang and H.-N. Yang [9] studied the bending stiffness on three connected specimens. Zheng et al. and so forth [10] analyzed the lug load capacity of 3D four-directional braided. Sun et al. [11] studied the carbon/epoxy 3D multidirectional braided tubular joints by finite element analysis. The current study provides numerical analysis methods, but most of them do not get tested, and no one do correlation analysis on the influence of the CFRP parameters in numerical simulation.
The mechanical properties of K type 3D braided joints and the influence of CFRP parameters were selected for research in this paper. Numerical analysis of K type 3D braided joint is studied firstly, verified through an existing test, and then various parameters were discussed including the effects of load parameters, geometric parameters, and material parameters on weighing the sensitivity impact of the joint ultimate bearing capacity.
2. Numerical Analysis Method Based on Micromechanics
2.1. Material Model
On the force along the fiber direction, Kelly and Davies put forward the hypothesis that all fibers have the same strength and more fragile than that of matrix (see Figure 1), if the composite has more than a minimum fiber volume content of , the composite is reaching its ultimate stress when the fiber deformation reaching its maximum strength. If the fiber strain along the fiber direction is equal to the matrix strain, the ultimate strength of the composite is In (1), is the ultimate stress of composite materials; is ultimate stress of fiber; and is the stress of ultimate matrix strain.
Because the fiber is brittle, it cannot exhibit elongation as matrix does. While the fiber damaged by longitudinal tension, composite material would be damaged, the strength of the composite is
In (2), is the fiber tensile strength; is the matrix stress which is equal to the matrix stress when fiber tensile strain reached limit; is cross-sectional area of the fiber; and is the cross-sectional area of the matrix.
For the 3D braided joints, 3-cells model which was presented by De-long and Shen [12] is accepted in engineering in the prediction of elastic constants for 3D braided composite material. The basic idea is as follows.
Calculate the elastic modulus of unidirectional composites.
Reference [12] gives the semiempirical formula which has been widely used in composite material area, as
In (3), , ; , ; , , respectively, are represented as the longitudinal elastic modulus, transverse elastic modulus; longitudinal shear modulus, transverse shear modulus; Poisson’s ratio, transverse Poisson’s ratio; , , , , , respectively, are represented as the longitudinal elastic modulus of fiber, transverse elastic modulus, longitudinal shear modulus and Transverse shear modulus, longitudinal Poisson’s ratio; , respectively, are represented as elastic modulus, shear modulus, Poisson’s ratio; is fiber volume content and is matrix volume content.
Calculate the axial flexibility matrix of unidirectional composites.
Unidirectional composite axial flexibility matrix is shown in (4). Among them, , , , , are determined, respectively, by (3). Consider:
Calculate the axial stiffness matrix of unidirectional composite
Calculate transformation matrix of unidirectional composite from the local coordinate to global coordinate. The Internal cell and fiber direction are shown in Figure 2.
In order to obtain the unified material constants, the stiffness matrix in the local coordinates ( coordinate) of unidirectional composites needs to be transformed into the global coordinate system. . Consider: , , , , , , , .
, are defined as shown in Figure 3 (for coordinates of the fiber as an example): as the fiber angle between the projection of unit cell longitudinal plane and longitudinal axis, as the fiber angle between the projection of unit cell transverse plane and horizontal axis, of the four fibers are , , , .
Calculate the unidirectional composite stiffness matrix of partial axis equivalent in the global coordinate system
Calculate the unit cell overall stiffness matrix.
Unit cell overall stiffness is obtained by the average matrix stiffness of all fiber direction, and this method is called stiffness average method. Consider:
Calculate the elastic constants of 3D braided composite materials by unit cell overall flexibility matrix.
3-cells model was used to predict the elastic constants of 3D braided composite material, which calculate stiffness matrix at each element according to the method which is shown in (3)–(8); then the whole stiffness matrix of the material could be calculated by stiffness matrix of each element according to the proportion of each unit cell by using the weighted sum method. Because of [13], unit cell stiffness matrix which occupies the largest proportion of the composite materials represents the stiffness matrix of braided material. Consider:
The elements of the flexibility matrix in whole unit cell are composed of elastic constants, as shown in (4). Therefore, the relationship between elastic constants of 3D braided composite material and the elements of unit cell overall flexibility matrix is as follows: , , are the elastic modulus of 3D braided composites; , , are the shear modulus of 3D braided composites; , , are the Poisson’s ratio of 3D braided composites; is row and column in the overall flexibility matrix .
After the MATLAB programming calculation, material elastic constants of this paper are , , , , and , . Test differs by 4.15% from values in [14]. (Unified listed in Table 1.)
2.2. Geometric Analysis Model
In the K-joint size as shown in Figure 4, the diameter of main pipe is 120 mm, 4 mm of thickness; the diameter of branch pipe is 54 mm, 3.2 mm of thickness.
The K-joint adopts four step 1 × 1 (making one transverse and longitudinal movement for each time) woven technology (as shown in Figure 5), which braided angle is 25°, fiber volume content is 59%. Properties of the matrix and fiber (is shown in Table 1).
(a) Braided indicate
(b) Initial state
(c) The first step in the state
(d) The second step in the state
(e) The third step in the state
(f) The fourth step in the state
2.3. Unit Type and Mesh Density
2.3.1. The Unit Type
SHELL93 and the SOLID92 were used to analyze the K-joint in this paper, the calculation results are shown in Table 2. It showed that the calculation results of 3D braided joint were less influenced by the type units. In order to improve the calculation efficiency, SHELL93 is used to analyze the K-joint in this paper.
2.3.2. The Mesh Dividing
Woven fibers have their own direction, branch and main pipe of the fiber direction differs, so as to overall performance. In order to simulate orthotropic character of material, the unit material properties should be corresponding to the fiber orientation. ANSYS provides the unit coordinate system to simulate the fiber direction, there are 2 ways: one is to mesh under different local coordinate; the other is to mesh together, then revise the unit coordinate to the local coordinate in unit character. The second method has been used in this paper, branch pipe divided unit coordinates is shown in Figure 6, in which the axial, tangential, and normal represent the unit coordinate three directions .
The mesh precision has influence on the calculation results. The K-joints are divided along in five ways: sparse, slightly sparse, moderate, slightly dense, and dense, quantitative indicators which are shown in Table 3.
It can be seen that the joints ultimate bearing capacity differs 2% when grid is slightly dense compared with the densest, however, the error of joints’ ultimate bearing capacity is larger when it divided sparsely. Considering the calculation of time and accuracy, the slightly dense grid is adopted in the paper.
2.4. Applying the Boundary Conditions and Load
The common used boundary conditions and loading method in analysis of K-joint is shown in Figure 7. In (a) one end of the main pipe is fixed, the other is connected with the sliding hinge, the branch pipe is hinged; in (b) the main pipe is fixed at both ends, the branch pipe is hinged; in (c) one end of the main pipe is fixed, the other is connected with the sliding hinge, the branch pipe is free; in (d) the main pipe is fixed at both ends, the branch is free. When the main pipe is stressed, it becomes two-way loading, and is one-way loading when it does not stress.
(a)
(b)
(c)
(d)
When loaded in one-way direction, the Loading-deformation curves in different boundary conditions are shown in Figure 8. Constraint on (a) and (b), (c) and (d) the results are similar, which illustrates constraints of the joints are less influenced the ultimate bearing capacity; constraints on (a) and (c), (b) and (d) differed greatly, which illustrates constraints at end of the main pipe have less influence on the joints ultimate bearing capacity.
When the branch is free, the ultimate bearing capacity of the joints is pretty less than the branch pipe which is hinged; at the same time branch pipe of K-joint deformation is larger, which has deviated from the K-plane, as shown in Figure 9. Considering the actual structure of the branch pipe ends restrained by abdominal rod, the main pipe could move to some extent, therefore the boundary conditions of constraint (a) has been chosen in this paper, considering the one-way and two-way loading modes.
2.5. Failure Criterion and Ultimate Bearing Capacity
Wang et al. [13] did statistic for usage of the composite materials failure criterion (as shown in Figure 10). The maximum stress criterion application frequency is ranked in second place, next to the maximum strain criterion. The maximum stress criterion was adopted in this paper, for this criterion is not only simple, practical, but also the test data of the maximum stress was verified in most of the reference for comparison.
is the maximum stress of Von. Mises of the joint, is the allowable material stress. is solved using Zuo [15] proposed 3D braided material strength of two order Tsai-Wu strength criterion based on prediction method and the maximum stress criterion. Prediction of the joint strength value is 610.7 MPa, the test value in the reference is 629 MPa, and the difference between the two is .
When the maximum equivalent stress is higher than the material allowable stress, joints damage, and this load can be seen as the joint bearing capacity.
2.6. Numerical Methods for Verifying the Accuracy
The lug tensile test in [10] as shown in Figure 11. Carbon fiber material for the test is T700-12K, matrix material is TDE-85 epoxy resin, the braided angle is 20°, and fiber volume content is 45%. Calculation results show that the materials elastic constants are , , , , , and . The material strength .
3D braided joint numerical analysis method was used to make numerical simulation for the lug, calculating the ultimate bearing capacity which is 74.9 kN, testing value which is 79.4 kN, which differed 5.66%, the maximum stress position and joint breaking position are basically the same (as shown in Figure 12); that is to say, the numerical simulation method achieves certain precision, which can be used to simulate the 3D braided joint.
3. Analysis on the Parameter Sensitivity of 3D Braided Joint
After the analysis on the mechanical properties of 3D braided joint, it should make parameter sensitivity analysis on the parameters affecting the bearing capacity.
3.1. The Parameters and Value of 3D Braided Joint
In this paper, the effect of load parameters, geometric parameters, and process parameters on the ultimate bearing capacity of was considered. The main parameters and their values are as shown in Table 4.
3.2. Analysis on the Global Parameter Sensitivity of 3D Braided Joint
According to whether to consider the interaction between parameters, parameter sensitivity analysis methods can be divided into local sensitivity analysis and global sensitivity analysis. Local sensitivity analysis refers to one parameter changed, and the other parameters remain unchanged, which can test the degree of influence on certain parameters variation to target results. Global sensitivity analysis refers to the influence on target results by changing one certain parameter when all the parameters changed. Because the global sensitivity analysis method considers the interaction between parameters, it is generally believed to be more accurate and scientific than the local sensitivity analysis method.
Since the parameter units of each node are not consistent, it needs to normalize the various parameters then make the parameters normalized to the range by using the Equation: [16], to get the curve of normalized parameters and target results. The Equation is as follows: —sensitivity coefficient; —variation target, taking the ultimate bearing capacity of as a target; —target reference value; —variations of parameter , and —designed domain of parameter .
Then it takes nonparametric statistical methods to make sensitivity analysis which is proposed by Marivoet and Saltelli in 1990. The method carried out regressing to analyze the parameters and the results, calculating the parameters and the standard variance and relative coefficient between the two, using the following equation to obtain sensitivity coefficients under different parameters: —coefficient parameter of parameter ,—coefficient parameter of parameter , —standard variance of parameter , and —standard variance of parameter .
This paper selected , , , , , , , and as the joints reference state.
Figure 13 shows an example of ultimate bearing capacity changed under different parameters. The slope of the curve is greater, the change of parameters has more influence on the ultimate bearing capacity, and the sensitivity is stronger.
Linear fitting for each parameter then take the average value, calculating sensitivity coefficient in Table 5.
Equation (13) defined the parameter (the relative contribution rate) as target sensitivity, that is, the ratio of either sensitivity coefficient or the sum of all the parameter sensitivity:
According to the definition above, the relative contribution of various parameters on the sensitivity of is the ratio of the sensitivity coefficient and the sum of all parameters sensitivity coefficient, which is shown in Figure 14.
4. The Conclusion
(1)By studying the ultimate bearing capacity of 3D braided composite joint with numerical analysis method from the material model, element type, mesh size, boundary conditions, loads, and failure criteria six aspects, and compared the 3D braided composites with the lug, the numerical analysis method is proved feasible in this essay.(2)Based on these results, the main results about the sensitivity analysis of parameters on the 3D braided composite joint is as follows: the main truss diameter thickness ratio fiber braiding angle load parameters competent diameter main and branch diameter ratio fiber volume content branch and main thickness ratio . The ratio of thickness between the branch and main pipe has small influence on the ultimate bearing capacity.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
The work described in this paper was jointly funded by the National Natural Science Foundation of China under Grant nos. 51378150 and 91215302.