编织复合材料的一种数值模型

提出编织复合材料的一种均匀化的数值模型,为工程应用提供了简便而有较高精度的数值计算方法。首先,对于编织复合材料的代表性单胞建立了精细有限元的胞元模型;其次,将单胞有限元分析结果引入基于三类变量广义变分原理的三维体罚单元,从而将一类有细观编织结构的复合材料等效为可用于宏观尺度计算的复合材料单元。通过数值算例将此方法的计算结果与实验和其他数值计算结果进行了比较,验证了该方法的有效性与实用性。     

Meshfree modeling and homogenization of 3D orthogonal woven composites

In this paper, evaluation of 3D orthogonal woven fabric composite elastic moduli is achieved by applying meshfree methods on the micromechanical model of the woven composites. A new, realistic and smooth fabric unit cell model of 3D orthogonal woven composite is presented. As an alternative to finite element method, meshfree methods show a notable advantage, which is the simplicity in meshing while modeling the matrix and different yarns. Radial basis function and moving kriging interpolation are used for the shape function constructions. The Galerkin method is employed in formulating the discretized system equations. The numerical results are compared with the finite element and the experimental results.     

Three-dimensional (3D) modeling of the thermoelastic behavior of woven glass fiber-reinforced resin matrix composites

The thermoelastic behavior of glass fiber-reinforced resin matrix composites is very important in several applications such as electronic packaging. Simulation of the composite behavior is complicated because of the complex nature of woven fiber architecture. In this study, we have conducted a numerical simulation of elastic and thermal expansion behavior of woven glass fiber-reinforced resin matrix composite. The simulations were compared to experimental data, showing excellent agreement with elastic properties and fairly good results for the thermal expansion coefficient of the composite.     

Simulation of Degraded Properties of 2D plain Woven C/SiC Composites under Preloading Oxidation Atmosphere

In this paper, a numerical model which incorporates the oxidation damage model and the finite element model of 2D plain woven composites is presented for simulation of the oxidation behaviors of 2D plain woven C/SiC composite under preloading oxidation atmosphere. The equal proportional reduction method is firstly proposed to calculate the residual moduli and strength of unidirectional C/SiC composite. The multi-scale method is developed to simulate the residual elastic moduli and strength of 2D plain woven C/SiC composite. The multi-scale method is able to accurately predict the residual elastic modulus and strength of the composite. Besides, the simulated residual elastic moduli and strength of 2D plain woven C/SiC composites under preloading oxidation atmosphere show good agreements with experimental results. Furthermore, the preload, oxidation time, temperature and fiber volume fractions of the composite are investigated to show their influences upon the residual elastic modulus and strength of 2D plain woven C/SiC composites.     

Compressive Behavior of 3D Woven Composite Stiffened Panels: Experimental and Numerical Study

The structural behavior and damage propagation of 3D woven composite stiffened panels with different woven patterns under axial-compression are here investigated. The panel is 2.5D interlock woven composites (2.5DIWC), while the straight-stiffeners are 3D woven orthogonal composites (3DWOC). They are coupled together with the Z-fibers from the stiffener passing straight thought the thickness of the panel. A “T-shape” model, in which the fiber bundle structure and resin matrix are drawn out to simulate the real situation of the connection area, is established to predict elastic constants and strength of the connection region. Based on Hashin failure criterion, a progressive damage model is carried out to simulate the compressive behavior of the stiffened panel. The 3D woven composite stiffened panels are manufactured using RTM process and then tested. A good agreement between experimental results and numerical predicted values for the compressive failure load is obtained. From initial damage to final collapse, the panel and stiffeners will not separate each other in the connection region. The main failure mode of 3D woven composite stiffened panels is compressive failure of fiber near the loading end corner.     

基于结构参数的平纹机织复合材料等效弹性性能预测

经向纤维束与纬向纤维束纵横交错引起的纤维弯曲(也称为波纹)是平纹机织复合材料固有特征。首先,提出了一种精确描述平纹机织复合材料单胞3D结构特征的数学表达式。其次,基于经典层合板理论和等应力假设,考虑平纹机织复合材料厚度方向非对称引起的弯曲-拉伸耦合效应及单胞结构特征,建立了含结构参数的平纹机织复合材料等效弹性性能多参数解析模型。通过数个典型算例验证了建立的多参数解析模型,结果表明：该多参数解析模型预测值与相关文献中有限元模型预测值、解析模型预测值、实验值等均吻合较好;该多参数解析模型预测值尤其是Z向弹性性能预测值,比文献中解析模型预测值更接近于实验值。在此基础上,进一步探讨了纤维束波纹比(包括纤维束波动方向波纹比与纤维束横截面波纹比)、经向与纬向纤维束构成的预成形体厚度、纤维束中弯曲部分的长度、相邻纤维束之间间距等结构参数对平纹机织复合材料弹性性能影响。该多参数解析模型建模方法为研究纺织复合材料力学性能提供了参考。     

Intra/inter-ply shear behaviors of continuous fiber reinforced thermoplastic composites in thermoforming processes

The thermoforming of continuous fiber reinforced thermoplastic (CFRTP) composite panels generally involves significant in-plane shear deformation. In the present work, the in-plane shear behavior of woven thermoplastic composites (Carbon/Polyphenylene Sulfide) over a range of processing temperatures is studied by bias-test experiments at different velocities. The experimental data of force versus displacement and force versus shear strain are presented for different extension velocities and temperatures. A thermo-visco-elastic model for numerical simulations of woven thermoplastic composite forming is proposed considering the influences of temperature and of strain rate. We applied a large displacement three-dimensional cohesive element with eight nodes which has been used for crack analysis in fracture mechanics by other authors, to investigate the inter-ply shear mechanism of woven thermoplastic composites. Applying three-dimensional cohesive elements, multi-plies forming simulations are performed to show inter-ply slip behaviors at different temperatures. The proposed models can be useful to predict from the properties of reinforcement and resin the intra/inter-ply shear behaviors of woven thermoplastic composites at high temperatures if experimental characterization of composite laminate behaviors is difficult to conduct.     

Spring stiffnesses of composite circular springs with extended flat contact surfaces under unidirectional line-loading and surface-loading configurations

Theoretical expressions based on the principle of minimum potential energy are presented which describe the stiffnesses of mid-surface symmetric, woven composite circular springs with extended flat contact surfaces subject to unidirectional line and surface-loading configurations. Three-dimensional finite element analysis has been employed to study the transverse shearing effects. Eighteen woven fibre/epoxy composite circular springs with extended flat contact surfaces in different radii and thicknesses were fabricated and tested. Comparison studies of the results obtained from both the analytical and numerical models are made with experimental data, and the results are found to be satisfactory. The semi-included angle of the flat contact surface is vital parameter to spring stiffnesses of the composite spring.     

A multi-scale finite element approach for modelling damage progression in woven composite structures

A significant challenge in the numerical modelling of composite structures with a multi-axis fibre architecture is the reproducibility of the textile mechanics [1]. A numerical analysis procedure for woven composite structures using a multi-scale finite element approach has been developed, and is presented in this paper. The approach is demonstrated for a flat two-dimensional woven glass/epoxy laminate. Digital microscopy is used to estimate tow cross-section and path, and quantify the amount of variation of these parameters. This data is used to generate both a meso-scale model of a single unit cell as well as a macro-scale model of the complete structure. Numerical results from the proposed approach are compared to experimental stress-strain data, which show good agreement in the lower strain range.     

Dynamic characterization of the thickness-tapered laminated plant fiber-reinforced polymer composite plates

Evolution of the laminated woven natural fiber fabric-reinforced polymer composite structures makes a way to the development of the non-uniform laminated composite structures in order to achieve the stiffness variation throughout the structure. An attempt is made in this work to carry out the experimental and numerical investigations on the dynamic characteristics of the thickness-tapered laminated woven jute/epoxy and woven aloe/epoxy composite plates. The governing differential equations of motion for the thickness-tapered laminated composite plate are developed using the h-p version FEM based on higher order shear deformation theory. The validation of the present finite element formulation is carried out by comparing the natural frequencies obtained using the finite element formulation with those natural frequencies determined experimentally. The developed model is further validated with the available literature works on tapered composite plate to confirm the efficiency of h-p version FEM. This work also explores the study of the vibrational characteristics of composite plates under the influence of plant fiber’s transverse isotropic material characteristics and porosity associated with plant fiber composites through the elastic constants evaluated in the author’s previous work. Also the influences of aspect ratios, ply orientations, and taper angles under various end conditions on the natural frequencies of the woven jute/epoxy composite plate are studied using the present finite element formulation. The forced vibration response of the thickness-tapered laminated woven jute/epoxy composite plate under the harmonic force excitation is carried out considering CFCF and CFFF end conditions.     

基于ANSYS环境的平面编织层合板拉伸破坏数值仿真

以ANSYS为平台编制了具有可移植性的APDL程序, 建立了损伤累积模型, 对平面编织层合板的损伤破坏行为进行了数值仿真。该模型对适合于单向铺层的Hanshin判据和Reddy刚度衰减方法进行了相应的修正。为验证模型的有效性, 对G803/5224平面编织光滑板、 孔板进行了相应的试验研究。结果表明, 该模型仿真结果与试验结果吻合, 并且比较简单直观, 为平面编织层合板的损伤扩展与破坏的研究提供了便于工程应用的数字化手段。      

碳纤维织物力学性能和冲压成形实验研究

在实际成形过程中,碳纤维复合材料往往处于复杂的应力状态,开展近于真实载荷环境下的力学试验分析,能够更准确地认识实际应用中材料的成形性能和变形机理.为获得碳纤维织物的基本力学特性,设计了平纹碳纤维织物拉伸试样及成形试样,进行了单轴拉伸、双轴拉伸、镜框剪切试验和方盒冲压成形实验研究,对比了不同双拉比及纱线取向对力学性能及成形性能的影响.研究结果表明:碳纤维织物具有高度的非线性、各向异性和双拉耦合特性,即经纬向纤维的力学性能会相互影响;剪切变形是成形过程中的主要变形模式,当剪切角达到临界锁死角时,织物发生起皱现象;同种织物不同纱线取向试样表现出不同的成形性能,因此可以根据零件几何形状选择合适纤维取向的织物,从而减少缺陷,优化成形零件的力学性能.研究结果为后续建立碳纤维织物本构模型和成形仿真奠定了基础.     

圆管状立体机织复合材料的多尺度分析

采用多尺度耦合的数值模型研究了圆管状立体机织复合材料的力学性能。建立了反映纤维束中纤维/基体二相材料的微观尺度单胞和反映周期性编织结构的细观尺度扇形单胞,并重点讨论了扇形单胞的周期性边界条件。通过逐级计算微观单胞、细观单胞的平均弹性常数,得到了圆管状立体机织复合材料的刚度参数,实现了由组分材料性能及编织参数预测圆管的宏观弹性性能,模型预测刚度与试验结果吻合。另一方面,研究了从大到小各尺度耦合的应力分析,对于圆管环向应力非周期分布的情况,建立了嵌入细观单胞的环状模型,进行了复杂荷载下从宏观圆管结构、到细观纤维束尺度、再到微观纤维尺度之间的逐级应力分析。     

考虑孔隙和微裂纹缺陷的C/C-SiC编织复合材料等效模量计算

通过观察C/C-SiC复合材料组元分布的扫描电子显微镜(SEM)照片 , 获得了C/C-SiC复合材料化学气相渗透(CVI)制备过程中产生孔隙和微裂纹的几何信息。在此基础上 , 建立了包含孔隙和微裂纹的C/C-SiC微结构有限元模型 , 并利用均匀化等效计算方法预测了平纹编织C/C-SiC复合材料的模量。针对CVI沉积方式制备的2组不同的C/C-SiC复合材料 , 实验测试与等效计算结果表明 : 基于 SEM照片建立的C/C-SiC纤维束和复合材料微结构有限元模型 , 能够反映CVI工艺制备C/C-SiC中孔隙和微裂纹的分布状况; 计算结果与实验数据有良好的一致性 , 数值计算可有效预测C/C-SiC编织复合材料的模量。      

Beyond plain weave fabrics – I. Geometrical model

In response to the large variety of weaving styles offered by the textile industry, a new general approach for the geometrical modeling of 2D biaxial orthogonal woven fabric reinforcements for composite materials is proposed here. New geometrical parameters are introduced in order to describe general families of twill and satin woven patterns, and a new classification of woven fabrics is proposed based on these parameters. Generation of the 3D internal geometry of the woven fabric families is achieved based on new geometrical functions that consider the actual configuration of the composite material in all its complexity. The proposed geometrical model is intended as the foundation for further analytical or numerical modeling of the mechanical properties of the composite materials reinforced with these fabrics.     

A unit cell approach of finite element calculation of ballistic impact damage of 3-D orthogonal woven composite

This paper presents ballistic impact damages of 3-D orthogonal woven composite in finite element analysis (FEA) and experimental. A unit-cell model of the 3-D woven composite was developed to define the material behavior and failure evolution. A user-defined subroutine VUAMT was compiled and connected with commercial available FEA code ABAQUS/Explicit to calculate the ballistic penetration. Ballistic impact tests were conducted to investigate impact damage of 3-D kevlar/glass hybrid woven composite. Residual velocities of conically-cylindrical steel projectiles (Type 56 in China Military Standard) and impact damage of the composite targets after ballistic perforation were compared both in theoretical and experimental. The reasonable agreements between FEA results and experimental results prove the validity of the unit-cell model in ballistic limit prediction of the 3-D woven composite. We believe such an effort could be extended to bulletproof armor design with the 3-D woven composite.     

Semi-continuous approach for the modeling of thin woven composite panels applied to oblique impacts on helicopter blades

In aeronautics, passenger safety and reliability of structures are essential aspects. In the specific case of helicopters, blades are subjected to impact loadings. Modeling these phenomena continue to be difficult and experimental tests often replace the prediction. The following work will focus on the experimental and numerical study of an oblique impact on the skin of the blade. It is equivalent in a first approach to an impact on a sandwich panel made up of a foam core and a thin woven composite skin. The objectives are to identify the mechanisms of damage in the skin for this kind of loading and to develop a representative modeling of the chronology of damage adapted to the modeling of the complete structure. Thus, a semi-continuous F.E. explicit modeling has been developed. It relies on the development of a specific damageable element at the bundles scale. Satisfactory numerical results are obtained. They allow the identification of the damage mechanism of the woven skin.     

The influence of distinct fiber arrangement on the thermo-mechanical behavior of steel fiber-reinforced thermoplastic matrix composites

The effect of different fiber arrangements on mechanical behavior was investigated by using both experimental study and finite elements analyses. In particular, this study examined resultant residual stresses and plastic strains of steel-fiber reinforced thermoplastic composite discs under constant convective air cooling conditions. Three composite discs were manufactured with an identical concentration of woven, circular and radial arrays. The thermal and mechanical properties of the composite discs were measured. The numerical and experimental cooling curves were converged to correctly describe the convective cooling condition of the finite element analyses. After the cooling, the residual stresses and plastic strains in each disc were compared with one another and the results were analyzed. No thermal residual stress or plastic strain was observed for the woven fiber array. Residual stress and plastic strain found in the circular fiber array was twice as high as those in the radial fiber array. It is concluded that the reinforcement fiber array of thermoplastic composites is an effective parameter to describe their thermo-mechanical properties for the formation of thermal residual stresses and plastic deformation.     

Shape and position effects of double holes on lateral buckling of cantilever composite beams

One of the most important mechanical behavior of composite beams subjected to certain external loads and boundary conditions is lateral buckling. The effects of hole dimension, shape and position, and beam thickness on the lateral buckling behavior of woven fabric laminated composite cantilever beams, having two square or two circular holes, were investigated. Firstly, the theoretical, experimental and numerical critical buckling loads of the beams without holes were found and compared with each other. It was shown that there is a good agreement among the theoretical, experimental and numerical results. ANSYS finite element (FEM) program was used for the numerical analyses. Therefore, the numerical analysis of some models with various hole dimensions, shapes (square or circular) and fiber directions were done by changing distance between the holes. It is concluded that the circular holes are advantageous compared to the square ones in terms of lateral buckling behavior.