Multi-scale modeling, stress and failure analyses of 3-D woven composites

The very complex, multi-level hierarchical construction of textile composites and their structural components commonly manifests via significant property variation even at the macro-level. The concept of a “meso-volume” (introduced by this author in early 1990s) is consistently applied in this work to 3-D stress/strain and failure analyses of 3-D woven composites at several levels of structural hierarchy. The meso-volume is defined as homogeneous, anisotropic block of composite material with effective elastic properties determined through volumetrically averaged 3-D stress and strain fields computed at a lower (“finer”) level of structural hierarchy and application of generalized Hooke’s law to the averaged fields. The meso-volume can represent a relatively large, homogenized section of a composite structural component, a lamina in laminated composite structure, a homogenized assembly of several textile composite unit cells, a single homogenized unit cell, a resin-impregnated yarn, a single carbon fiber, even a carbon nanotube assembly. When composed together, distinct meso-volumes constitute a 3-D Mosaic model at the respective hierarchy level. A multi-scale methodology presented in this paper first illustrates 3-D stress/strain analysis of the Mosaic unidirectional composite, computation of its effective elastic properties and their further use in 3-D stress/strain analysis of the Mosaic model of 3-D woven composite Unit Cell. The obtained 3-D stress/strain fields are then volumetrically averaged within the Unit Cell, and its effective elastic properties are computed. The predicted effective elastic properties of 3-D woven composite are compared with experimental data and show very good agreement. Further, those effective elastic properties are used in 3-D simulations of three-point bending tests of 3-D woven composite; theoretical predictions for central deflection show excellent agreement with experimental data. Finally, a 3-D progressive failure analysis of generic 3-D Mosaic structure is developed using ultimate strain criterion and illustrated on the 3-D woven composite Unit Cell. The predicted strength values are compared to experimental results. The presented comparisons of theoretical and experimental results validate the adequacy and accuracy of the developed material models, mathematical algorithms, and computational tools.     

Process parameters design of a three-dimensional and five-directional braided composite joint based on finite element analysis

This paper presents an analytical method for designing the configuration of composite joint with three-dimensional (3D) five-directional braided composites. Based on the analysis of 3D braided structure characteristics, the elastic properties of the 3D five-directional braided composites were determined by the volume averaging method. The effects of the braiding angle and fiber volume fraction on the elastic constants of the braided composites were also discussed. Finite element analysis on the load capacity of the 3D five-directional braided composite joint was implemented using the software ANSYS Workbench 14.0. The influence of braiding angle on the stress, strain and deformation of the composite joint under tensile loading were calculated. The results show that when the fiber volume fraction of the 3D five-directional braided preform is given, the equivalent stress of the composite joint decreases monotonically as the braiding angle increases, while the normal stress, maximum principal stress and total deformation firstly decreases and then increases. Based on the finite element analysis, we found that at the fiber volume fraction of 60%, the braiding angle within the range of 30–35° are the optimum processing parameters for the 3D five-directional braided composite joint structure that used in the tensile load 320 N condition.     

In-plane Compressive Behaviors of 3-D Textile Composites at Various Strain Rates

The in-plane compressive behaviors of 3-D textile composites, which including 3-D woven composite, multi-axial multi-layer warp knitted (MMWK) composite and 3-D braided composite, were studied at quasi-static and high strain rate compression loading. The compression behaviors at high strain rates (600∼2,500/s) were tested with split Hopkinson pressure bar (SHPB). The quasi-static compressive tests were conducted on a MTS 810.23 tester and compared with the results at high strain rates. The comparisons indicate that the compression stiffness, failure stress and failure strain for the three kinds of 3-D textile composites are sensitive to strain rate. The MMWK composite has higher failure stress than the 3-D woven composite and 3-D braided composite at the same strain rate; however, the failure strain of the 3-D braided composites is higher than that of the 3-D woven composite and 3-D knitted composite at quasi-static compression because of the quasi-isotropic structure feature in the 3-D braided composite. The compressive failure modes of the 3-D woven composite, MMWK composite and 3-D braided composite are totally different because of the different preform structure.     

A meso-scale finite element model for simulating free-edge effect in carbon/epoxy textile composite

Textile composites are well known for their excellent through thickness properties and impact resistance. In this study, a representative unit cell model of a triaxial braided composite is developed based on the composite fiber volume ratio, specimen thickness and microscopic image analysis. A meso-scale finite element (FE) mesh is generated based on the detailed unit cell dimensions and fiber bundle geometry parameters. The fiber bundles are modeled as unidirectional fiber reinforced composites. A micromechanical finite element model was developed to predict the elastic and strength material properties of each unidirectional composite by imposing correct boundary conditions that can simulate the actual deformation within the braided composite. These details are then applied in the meso-mechanical finite element model for a 0°/+60°/−60° triaxially braided T700s/E862 carbon/epoxy composite. Model correlations are conducted by comparing numerical predicted and experimental measured axial tension and transverse tension response of a straight-sided, single-layer (one ply thick) coupon. By applying a periodic boundary condition in the loading direction, the meso model captures the local damage initiation and global failure behavior, as well as the periodic free-edge warping effect. The failure mechanisms are studied using the field damage initiation contours and local stress history. The influence of free-edge effect on the failure behaviors is investigated. The numerical study results reveal that this meso model is capable of predicting free-edge effect and allows identification of its impact on the composite response.     

Effective elastic properties of solids with irregularly shaped inclusions

A unit rectangular cell is usually cut out from a medium for investigating fracture mechanism and elastic properties of the medium containing an array of irregularly shaped inclusions. It is desirable to clarify the geometrical parameters controlling the elastic properties of heterogeneous materials because they are usually embedded with randomly distributed particulate. The stress and strain relationship of the rectangular cell is obtained by an ad hoc hybrid-stress finite element method. By matching the boundary condition requirements, the effective elastic properties of composite materials are then calculated, and the effect of shape and arrangement of inclusions on the effective elastic properties is subsequently considered by the application of the ad hoc hybrid-stress finite element method through examining three types of rectangular cell models assuming rectangular arrays of rectangular or diamond inclusions. It is found that the area fraction (the ratio of the inclusion area over the rectangular cell area) is one dominant parameter controlling the effective elastic properties.     

Homogenization of braided fabric composite for reliable large deformation analysis of reinforced rubber hose

High pressure rubber hose is in the lamination structure composed of pure rubbers and braided fabric composite layers to have the sufficient strength against the excessive radial expansion and the large deformation, in which the braided fabric layer is woven with wrap and fill tows inclined to each other with the predefined helix angle in the complex periodic pattern. The consideration of detailed geometry of braided fabric layer in the numerical analysis leads to a huge number of finite elements so that the braided fabric layer has been traditionally simplified as an isotropic cylindrical one with the homogeneous isotropic material properties of braid spun tread. However, this simple model leads to the numerical prediction and design with the questionable reliability. In this context, this paper addresses the development of an in-house module, which is able to be interfaced with commercial FEM code, for the reliable large deformation analysis of the reinforced rubber hose with the element number at the level of the traditional simple model. The in-house module is able to not only automatically generate 3-D unit cell (or RVE) model of the braided fabric layer but evaluate the homogenized orthotropic material properties by automatically performing a serious of unit cell finite element analyses based on the superposition method. The validity of the in-house module and the reliability of the homogenization method are verified through the illustrative numerical experiments.     

三维四向编织复合材料力学性能的有限元分析

在已有研究的基础上,提出了一个新的三维编织复合材料单元胞体模型,该模型正确地反映了纤维束的交织方式,十分接近三维编织复合材料的真实结构,可用于三维四向编织复合材料有效模量的有限元数值预报,并合理确定复合材料内部全场应力分布。采用有限元软件对该模型进行了力学分析,得到了相关等效弹性性能参数。结果表明:有限元计算得到的三维编织复合材料的等效弹性性能与实验结果和理论预测值都吻合较好,从而验证了该模型的有效性。此外,基于新的单元胞体模型还确定了三维四向编织复合材料的应力场,为进一步的强度计算奠定了基础。      

Eigenstrain and Fourier series for evaluation of elastic local fields and effective properties of periodic composites

The elastic stress and strain fields and effective elasticity of periodic composite materials are determined by imposing a periodic eigenstrain on a homogeneous solid, which is constrained to be equivalent to the heterogeneous composite material through the imposition of a consistency condition. To this end, the variables of the problem are represented by Fourier series and the consistency condition is written in the Fourier space providing the system of equations to solve. The proposed method can be considered versatile as it allows determining stress and strain fields in micro-scale and overall properties of composites with different kinds of inclusions and defects. In the present work, the method is applied to multi-phase composites containing long fibers with circular transverse section. Numerical solutions provided by the proposed method are compared with finite element results for both unit cell containing a single fiber and unit cell with multiple fibers of different sizes.     

Effect of microscopic damage events on static and ballistic impact strength of triaxial braid composites

The reliability of impact simulations for aircraft components made with triaxial braided carbon fiber composites is currently limited by inadequate material property data and lack of validated material models for analysis. Methods to characterize the material properties used in the analytical models from a systematically obtained set of test data are also lacking. A macroscopic finite element based analytical model to analyze the impact response of these materials has been developed. The stiffness and strength properties utilized in the material model are obtained from a set of quasi-static in-plane tension, compression and shear coupon level tests. Full-field optical strain measurement techniques are applied in the testing, and the results are used to help in characterizing the model. The unit cell of the braided composite is modeled as a series of shell elements, where each element is modeled as a laminated composite. The braided architecture can thus be approximated within the analytical model. The transient dynamic finite element code LS-DYNA is utilized to conduct the finite element simulations, and an internal LS-DYNA constitutive model is utilized in the analysis. Methods to obtain the stiffness and strength properties required by the constitutive model from the available test data are developed. Simulations of quasi-static coupon tests and impact tests of a represented braided composite are conducted. Overall, the developed method shows promise, but improvements that are needed in test and analysis methods for better predictive capability are examined.     

周期性单胞复合材料有效弹性性能的边界力方法

均匀化方法是一种适应于周期性构造复合材料有效性能预测的有效方法。然而均匀化方程数学表达形式复杂, 均匀化方法很难直接应用通用有限元软件进行计算, 因此本文中提出一种便于求解均匀化方程的边界力方法, 利用高斯定理将原均匀化问题转化为普通的三维应力问题, 给出了单胞中不同材料交界面上作用的面分布力形式。运用有限元软件求解了均匀化系数, 预报了单向复合材料和三维四向编织复合材料的有效性能, 计算结果与实验吻合很好。