二维二轴1×1编织复合材料细观结构模型及力学性能有限元分析

考虑纤维束相互挤压及横截面形状变化, 采用纤维束截面六边形假设, 建立了二维二轴1×1编织复合材料的参数化单胞结构模型。通过引入周期性位移边界条件, 基于细观有限元方法, 对编织材料的弹性性能进行预测, 讨论了编织角及纤维体积含量对面内弹性常数的影响, 并分析了典型载荷下单胞细观应力场分布。研究表明: 单胞结构模型有效反映了纤维束的空间构型和交织特征, 实现了不同编织工艺参数下模型的快速建立; 基于单胞有限元模型的弹性性能预测结果与试验结果较为吻合; 模型给出了单胞合理的应力场分布, 为二维编织复合材料的结构优化和损伤预测奠定基础。      

四步法三维矩形编织复合材料的细观结构模型

基于现有实验研究和编织工艺中携纱器的运动规律, 重点分析了材料内部区域纤维束的空间构型, 建立了材料的三维实体细观结构模型。该模型不仅体现了内部纱线因打紧工序而形成的紧密接触和截面变形, 而且考虑了内部和表面区域纱线因挤紧状态不同所造成的纱线填充因子变化。基于一种单胞取向平行于材料横截面边界方向的新划分方案, 解决了45°单胞划分方案的不足, 建立了便于力学性能分析的单胞几何模型, 并指出了编织工艺参数和模型宏细观结构参数的关系。模型数值结果与试件实测数据吻合, 表明了该模型的合理有效性,为材料后续力学性能分析奠定了基础。      

三维机织复合材料的一种梁单元细观力学模型

根据三维机织复合材料中纤维束排列和变形的周期性特点,推导了一种细观梁单元模型。该模型考虑了纤维束的拉 (压) 弯耦合效应和纤维束之间的相互作用,可以描述纤维束和基体中的细观应力分布,并得到宏观的力学性能。针对一种典型的三维机织复合材料,首先根据编织参数,确定其细观几何结构,取最小周期的一段纤维束作为分析胞元,用上述细观梁单元分析了该段纤维束在面内拉伸荷载下的细观应力分布,计算出平均模量, 并用材料试验和细观实体有限元对本模型进行了检验,结果与本文的预测吻合良好。研究表明,拉、弯耦合效应引起的纤维束中的细观弯曲应力同平均轴向应力相比,不可忽略。      

含纤维束截面形状变化的三维编织复合材料细观模型及刚度预报

基于实验观察和理论研究, 重点分析了材料内部区域纤维束的空间构型, 建立了一个新的三维实体细观结构模型, 并指出了编织工艺参数和模型细观结构参数之间的关系。该模型较真实地反映了纤维束之间的相互挤压变形方式, 纤维束横截面积沿纤维束轴向不断变化, 更符合三维四向编织复合材料的实际结构。基于刚度体积平均及柔度体积平均混合思想, 建立了相应的刚度预报模型。用该模型计算编织复合材料几何特性及工程弹性常数的数值结果与试件实测数据吻合, 表明了该模型的合理有效性, 为进一步研究三维编织复合材料的拉伸强度及破坏机制提供了基础。     

三维四向编织复合材料有效性能的预报

根据三维四向编织复合材料中纤维束的空间几何结构特征, 建立了比较合理的三胞模型。模型中考虑了3种单胞各自纤维束的空间结构和弯曲, 同时引入纤维束填充因子来描述各类单胞中纤维束的不同截面形状对材料弹性常数的影响。基于刚度体平均方法, 建立了相应的刚度预报模型, 得到了三维四向编织复合材料的工程弹性常数。用细观力学方法分析了工艺参数和尺寸效应对材料有效性能的影响规律。不同尺寸试件的数值预报结果和实验结果吻合较好。      

三维四向编织复合材料界面损伤机理数值分析

考虑界面脱粘表面压应力下摩擦力对材料界面力学性能的影响,建立损伤-摩擦相结合的界面本构模型,编写用户材料子程序VUMAT,实现其在有限元软件ABAQUS中的嵌入。基于周期性胞元分析思想,在单胞模型中纤维束/基体、纤维束/纤维束分界面引入界面单元,结合损伤-摩擦相结合的界面本构模型,建立含界面相三维四向编织复合材料的细观有限元模型。模拟典型载荷下界面损伤的起始和扩展过程,分析界面应力传递和界面破坏机理,研究界面性能对复合材料宏细观力学性能的影响规律,为实现三维四向编织复合材料界面性能优化设计和控制提供参考。      

2.5D机织复合材料的渐进损伤与失效模拟

为了研究典型2.5D机织复合材料的压缩性能,开展了复合材料单胞结构的经向和纬向压缩实验,并通过对材料编织结构的细观表征,建立了细观尺度的单胞有限元模型来模拟压缩载荷下单胞内部的变形及渐进失效过程。结果表明,2.5D机织复合材料在受压时表现出明显的非线性力学响应,材料沿经向的压缩模量和强度均高于纬向;经向压缩时材料的主要破坏模式有经纱的横向开裂、纤维束间的界面分层破坏、纬纱的压溃及基体的开裂,纬向压缩时出现的主要破坏模式是纬纱的压溃、纬纱纤维束的断裂及基体开裂;通过对比试验与有限元结果,认为所建立的细观有限元模型能够准确预测材料单胞在压缩载荷下的应力-应变响应,并且能够模拟编织结构中的损伤起始和演化过程。      

二维二轴编织复合材料几何模型及弹性性能预测

提出了二维二轴1×1和2×2编织复合材料的几何模型,模型考虑了纤维束的相互挤压及横截面的变化。基于细观分析和体积平均法,建立了预测二维二轴编织复合材料弹性性能的理论分析方法。数值结果与试验结果吻合,表明该方法行之有效,且具有运算快、精度高、适合工程分析等优点。分析了编织角、纤维体积含量和纤维束横截面形状对材料弹性常数的影响。研究表明,编织角对弹性常数的影响具有互补性,材料弹性模量与纤维体积含量成正比,纤维束截面形状变化对材料弹性常数影响不大。     

三维四步编织复合材料刚度和强度的数值模拟

为了研究三维四步法编织复合材料的力学性能,利用ANSYS有限元软件对材料的细观体胞模型进行数值模拟,计算三维编织复合材料的宏观弹性常数,讨论了纤维编织角和体积比对弹性常数的影响。采用不同的强度准则分别对纤维束和基体材料进行强度校核,从而得到材料发生破坏时失效单元的体积百分比。根据失效单元的分布情况分析材料的破坏机理,进而预报材料的拉伸强度。模拟计算结果与实验值吻合较好。     

考虑纤维束间粘接层的三维四向编织复合材料弹性性能数值预测

考虑了相邻纤维束之间的界面粘接效应, 建立了考虑纤维束间粘接层的三维四向编织复合材料单胞有限元模型, 较为真实地模拟了该材料的细观结构, 讨论了相应的边界条件和约束条件, 并采用有限元方法计算了该材料的等效弹性性能参数, 计算结果与实验值符合较好。研究了等效弹性性能参数随不同编织角及体积分数的变化关系, 得到了体胞的细观应力场, 为强度分析提供数据。      

An inverse method for material parameters determination of titanium samples under tensile loading

Different approaches used for the simulation of woven reinforcement forming are investigated. Especially several methods based on finite element approximation are presented. Some are based on continuous modelling, while others, called discrete or mesoscopic approaches, model the components of the fabric. A semi discrete finite element made of woven unit cells under biaxial tension and in-plane shear is detailed. In continuous approaches, the difficulty lies in the necessity to take the strong specificity of the fibrous material into account. The yarn directions must be strictly followed during the large strains of the fabric. This is the main goal of the non-orthogonal model and of the hypoelastic constitutive model based on the yarn rotation presented in this paper. In the case of discrete and semi-discrete approaches the directions of the yarns are “naturally” followed because the yarns are modeled. Explicitly, however, modeling each component at the mesoscopic scale can lead to high numerical cost.     

The effect of yarn distortion on the mechanical properties of 3D four-directional braided composites

A representative volume cell (RVC) is chosen to analyze the mechanical properties of 3D (3 dimensions) four-directional braided composites. Owing to braid yarns (an assembly of fibers) squeezing against each other in the braided composites, the braid yarns are distorted. Based on geometrical characteristics of the braided composites, cross-section of each braid yarn is supposed to be an octagon and divided into seven regions in the RVC. The distortion characteristics of yarns are considered in each region. Elastic properties of each region obtained by stochastic function theory are introduced into finite element model to calculate the mechanical properties of the RVC. The influences of yarn distortion on the stiffness and strength of the braided composites are obtained and discussed.     

A homogenization procedure for the numerical analysis of woven fabric composites

The paper presents a procedure for the numerical evaluation of the mechanical properties of woven fabric laminates. Woven fabrics usually present orthogonal interlaced yarns (warp and weft) and distribution of the fibers in the yarns and of the yarns in the composite may be considered regular. This allows us to apply the homogenization theory for periodic media both to the yarn and to the fabric. Three-dimensional finite element models are used in two steps to predict both the stiffness and the strength of woven fabric laminates. The model includes all the important parameters that influence the mechanical behavior: the lamina thickness, the yarn orientation, the fiber volume fraction and the mechanical characteristics of the components. The capabilities of the numerical model were verified studying the elastic behavior of a woven fabric laminate available in the literature and the ultimate strength of a glass fabric laminate experimentally investigated. The procedure, that can be implemented into commercial finite element codes, appears to be an efficient tool for the design of textile composites.     

Effect of Frictions on the Ballistic Performance of a 3D Warp Interlock Fabric: Numerical Analysis

3D interlock woven fabrics are promising materials to replace the 2D structures in the field of ballistic protection. The structural complexity of this material caused many difficulties in numerical modeling. This paper presents a new tool that permits to generate a geometry model of any woven fabric, then, mesh this model in shell or solid elements, and apply the mechanical properties of yarns to them. The tool shows many advantages over existing software. It is very handy in use with an organization of the functions in menu and using a graphic interface. It can describe correctly the geometry of all textile woven fabrics. With this tool, the orientation of the local axes of finite elements following the yarn direction facilitates defining the yarn mechanical properties in a numerical model. This tool can be largely applied because it is compatible with popular finite element codes such as Abaqus, Ansys, Radioss etc. Thanks to this tool, a finite element model was carried out to describe a ballistic impact on a 3D warp interlock Kevlar KM2? fabric. This work focuses on studying the effect of friction onto the ballistic impact behavior of this textile interlock structure. Results showed that the friction among yarns affects considerably on the impact behavior of this fabric. The effect of the friction between projectile and yarn is less important. The friction plays an important role in keeping the fabric structural stability during the impact event. This phenomenon explained why the projectile is easier to penetrate this 3D warp interlock fabric in the no-friction case. This result also indicates that the ballistic performance of the interlock woven fabrics can be improved by using fibers with great friction coefficients.     

三维机织复合材料残余应力/应变多尺度分析及工艺参数优化

为预测三维机织复合材料工艺引入的残余应力/应变,提出工艺制度优化方案,建立了一种工艺过程分析的多尺度模型.通过建立纤维尺度及纱线尺度代表体元(RVE),计算了成型过程中纤维纱线及三维机织复合材料的模量演化历程.考虑固化过程中树脂的化学收缩效应,在纱线尺度上开展热-化学-力学耦合分析,预测了细观残余应力-应变及其演化规律...     

Tensile Properties and Failure Mechanism of a New 3D Nonorthogonal Woven Composite Material

Tensile properties and failure mechanism of a newly developed three-dimensional (3D) woven composite material named 3D nonorthogonal woven composite are investigated in this paper. The microstructure of the composite is studied and the tensile properties are obtained by quasi-static tensile tests. The failure mechanism of specimen is discussed based on observation of the fracture surfaces via electron microscope. It is found that the specimens always split along the oblique yarns and produce typical v-shaped fracture surfaces. The representative volume cell (RVC) is established based on the microstructure. A finite element analysis is conducted with periodical boundary conditions. The finite element simulation results agree well with the experimental data. By analyzing deformation and stress distribution under different loading conditions, it is demonstrated that finite element model based on RVC is valid in predicting tensile properties of 3D nonorthogonal woven composites. Stress distribution shows that the oblique yarns and warp yarns oriented along the x direction carry primary load under x tension and that warp yarns bear primary load under y tension.     

A continuum shell finite element model for impact simulation of woven fabrics

A new computational approach is developed to predict the impact behaviour of fabric panels based on the detailed response of the smallest repeating unit (unit cell) in the fabric. The unit cell is constructed and calibrated using measured geometrical (weave architecture, crimp, voids, etc.) and mechanical properties of the fabric. A pre-processor is developed to create a 3D finite element mesh of the unit cell using the measured fabric cross-sectional micro-images. To render an efficient method for simulation of multi-layer packs, these unit cells are replaced with orthotropic shell elements that have similar macroscopic (smeared) mechanical properties as the unit cell. The aim is to capture the essence of the response of a unit cell in a single representative shell element, which would replace the more complicated and numerically costly 3D solid model of the yarns in a crossover. The 3D finite element analysis of the unit cell is used to provide a baseline mechanical response for calibrating the constitutive model in the equivalent shell representation. This shell element takes advantage of a simple physics-based analytical relationship to predict the behaviour of the fabric's warp and weft yarns under general applied displacements in these directions. The analytical model is implemented in the commercial explicit finite element code, LS-DYNA, as a user material routine (UMAT) for shell elements. Layers of fabric constructed from these specialized elements are stacked together to create fabric targets that are then analysed under projectile impact. This approach provides an efficient numerical model for the dynamic analysis of multi-layer fabric structures while taking into account several geometrical and material attributes of the yarns and the fabric.     

Multiscale modeling of the impact of textile fabrics based on hybrid element analysis

In this work, a multi scale modeling approach has been developed to simulate the impact of woven fabrics using a finite element (FE) analysis. A yarn level of resolution is used in the model. This approach, referred to as the hybrid element analysis (HEA) is based on decreasing the complexity of the finite element model with distance away from the impact zone based on the multiscale nature of the fabric architecture and the physics of the impact event. Solid elements are used to discretize the yarns around the impact region, which transition to shell elements in the surrounding region. A new method for modeling the shell yarns is incorporated that more accurately represents the contours of the yarn cross section. Impedances have been matched across the solid–shell interface to prevent interfacial reflections of the longitudinal strain wave. The HEA method is validated by first applying it to the FE model of a single yarn for which an analytical solution is known. The HEA method is then applied to a woven fabric model and validated by comparing it against a baseline model consisting of yarns discretized using only solid elements.     

2.5D机织复合材料压缩性能实验与数值模拟

为了研究2.5D机织复合材料的压缩损伤和失效机制,验证双尺度渐进损伤有限元数值模拟方法的有效性,对这类复合材料分别沿经纱方向和纬纱方向进行了准静态压缩实验,获得了其相应的应力-应变曲线,并测定了材料的初始弹性模量和极限强度。在此基础上,利用双尺度渐进损伤有限元数值方法模拟分析了材料的压缩应力-应变响应和损伤演化行为,取得了与实验吻合较好的模拟结果。结果表明:2.5D机织复合材料在纬向压缩下的主要失效模式是纬纱的轴向压溃与断裂,可获得相对较高的压缩强度;但在经向压缩下,经纱因弯曲会承受附加弯矩作用,从而对周围基体造成挤压,故在经纱轴向断裂之前容易出现经纱之间基体的压溃和纱线之间的分层开裂,使强度降低,不利于发挥纤维的承载优势。