Research Article  Open Access
Yue Wu, Bo Nan, Liang Chen, "Mechanical Performance and Parameter Sensitivity Analysis of 3D Braided Composites Joints", The Scientific World Journal, vol. 2014, Article ID 476262, 9 pages, 2014. https://doi.org/10.1155/2014/476262
Mechanical Performance and Parameter Sensitivity Analysis of 3D Braided Composites Joints
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, LayerMolding 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 3D 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 3D 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 fourdirectional 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 crosssectional area of the fiber; and is the crosssectional area of the matrix.
For the 3D braided joints, 3cells model which was presented by Delong 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.
3cells 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 Kjoint 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 Kjoint 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 Kjoint 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 Kjoint 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 Kjoints 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 Kjoint 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 twoway loading, and is oneway loading when it does not stress.
(a)
(b)
(c)
(d)
When loaded in oneway direction, the Loadingdeformation 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 Kjoint deformation is larger, which has deviated from the Kplane, 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 oneway and twoway 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 TsaiWu 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 T70012K, matrix material is TDE85 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.
 
: The ratio of pipe force and material strength. : The ratio of main pipe diameter and thickness. : The ratio of main pipe diameter and branch pipe diameter. : Branch pipe axis and main pipe axis angle (°). : The ratio of main pipe wall thickness and branch pipe wall thickness. : Fiber braiding angle (°). : Fiber volume content (%). : The main pipe diameter (mm). 
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.
References
 J. W. H. Yap, M. L. Scott, R. S. Thomson, and D. Hachenberg, “The analysis of skintostiffener debonding in composite aerospace structures,” Composite Structures, vol. 57, no. 14, pp. 425–435, 2002. View at: Publisher Site  Google Scholar
 A. Mills, “Automation of carbon fibre preform manufacture for affordable aerospace applications,” Composites A: Applied Science and Manufacturing, vol. 32, no. 7, pp. 955–962, 2001. View at: Publisher Site  Google Scholar
 S. Ju, J. Zeng, and D. Jiang, “Study progress in composite trussjoint,” Materials Review, vol. 20, no. 12, pp. 28–32, 2006 (Chinese). View at: Google Scholar
 Q. Guo, G. Zhang, and J. Li, “Process parameters design of a threedimensional and fivedirectional braided composite joint based on finite element analysis,” Materials & Design, vol. 46, pp. 291–300, 2013. View at: Publisher Site  Google Scholar
 J. M. Yang, C. L. Ma, and T. W. Chou, “Fiber inclination model of threedimensional textile structural composites,” Journal of Composite Materials, vol. 20, no. 5, pp. 472–484, 1986. View at: Google Scholar
 L. D. Zhou and Z. Zhuang, “Strength analysis of threedimensional braided Tshaped composite structures,” Composite Structures, vol. 104, pp. 162–168, 2013. View at: Publisher Site  Google Scholar
 D.L. Wu and Z.P. Hao, “5D braided structural composites,” Journal of Astronautics, vol. 3, no. 7, pp. 40–50, 1993 (Chinese). View at: Google Scholar
 B. Zheng, S. Zhang, P. He, J. Ju, and Y. Sun, “Experiments and analyses for carbon fiber reinforced composite multijoints of the satellite antenna truss,” Acta Materiae Compositae Sinica, vol. 22, no. 6, pp. 172–177, 2005 (Chinese). View at: Google Scholar
 C.Y. Yang and H.N. Yang, “Bending rigidity of a satellite antenna truss joint made of 3D woven composites,” Materials Science & Technology, vol. 16, no. 6, pp. 810–813, 2008. View at: Google Scholar
 X. Zheng, Y. Guo, Q. Sun, Y. Chai, and X. Li, “Strength predication for loadbearing joints of threedimensional braided composites,” Journal of Mechanical Strength, vol. 28, no. 6, pp. 923–926, 2006 (Chinese). View at: Publisher Site  Google Scholar
 Y. Sun, L. Chen, and J.L. Li, “Finite element analysis on loading properties of intersecting tubular joint with carbon/epoxy 3D multidirectional braided composite,” Journal of Solid Rocket Technology, vol. 31, no. 3, pp. 266–274, 2008 (Chinese). View at: Google Scholar
 W. U. Delong and H. R. Shen, “On the study of mechanical behavior of textile structural composites,” Advances In Mechanics, vol. 31, no. 4, pp. 583–591, 2011 (Chinese). View at: Google Scholar
 B. Wang, S. Wu, and J. Liang, “Strength theory and failure analysis of composite materials,” Failure Analysis and Prevention, vol. 1, no. 2, pp. 13–19, 2006. View at: Google Scholar
 B. Wang, Study on the Mechanical Behavior of Three Dimensional Braided Composites, Northwestern Polytechnical University, Xi'an, China, 2003 (Chinese).
 W.W. Zuo, Research on Mechanical Properties and Engineering Application of 3D Braided Composites, Huazhong University of Science & Technology, Wuhan, China, 2006 (Chinese).
 M.H. Wei and D.P. Yu, “Landslide susceptibility factor analysis and stability evaluation,” China Water Transport, vol. 5, no. 8, pp. 97–98, 2007 (Chinese). View at: Google Scholar
Copyright
Copyright © 2014 Yue Wu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.