编织复合材料弹性性能的细观力学模型

提出了编织复合材料弹性性能分析的细观力学模型，这个力学模型考虑了实际编织结构中的纬向和经向纤维束的曲屈，相邻纤维束之间的间隙和纤维束的横截面尺寸对编织复合材料弹性性能的影响，并探讨了在纤维束间纯树脂区内孔隙的含量和两种叠层结构对材料弹性性能的影响．理论计算结果与实测值的比较，表明所提出的细观力学模型是合理的．根据理论分析的结果，提出了优化单层和叠层编织结构的结构参数选择方法     

三维编织复合材料渐进损伤的非线性数值分析

基于考虑纤维束相互挤压的八边形纤维束截面单胞模型，引入周期性位移边界条件，采用 细观非线性有限元方法，建立了三维四向编织复合材料的渐进损伤拉伸强度模型. 该模型考 虑了增强体纤维束纵向非线性剪切应力-应变关系，采用Hashin型损伤失效准则定义了纤维 束的典型损伤类型，并根据纤维束和纯基体相应损伤类型所造成的材料性能退化，模拟了不 同编织角试件各类损伤产生、扩展及材料最终破坏的整个过程. 模型数值结果与实验数据吻 合较好，证明了该模型的合理有效性. 探讨了组分材料剪切非线性、损伤对材料宏观非线性 本构行为的影响，结果表明：随着编织角增大，纤维束剪切非线性效应和累积损伤对材料非 线性力学行为的影响明显增强.     

树脂基三维编织复合材料粘弹性性能的有限元分析

根据材料的细观结构,采用APDL语言分别建立了纤维束和三维编织复合材料两级单胞的参数化几何模型;推导了Prony级数表示的树脂粘弹性本构方程,对模型进行了组分材料参数设置;对纤维束单胞模型进行扫掠式网格划分,对三维编织复合材料单胞模型进行线-面-体式网格划分;对两级单胞模型均施加合理的边界条件,使单胞边界上的位移满足周期性和连续性。以有限元模型为基础,计算了三维编织复合材料的粘弹性能,并给出了材料粘弹性效应随工艺参数变化的规律。计算结果表明:三维编织复合材料编织方向的粘弹性效应随编织角的增大而增强,随纤维体积比的增大而减弱。该结果与已有实验结论一致。     

基于渐进均匀化的平纹编织复合材料低速冲击多尺度方法

针对平纹编织复合材料低速冲击响应和损伤问题,提出了一种多尺度分析方法. 首先, 建立微观尺度单胞模型,引入周期性边界条件,采用最大主应力失效准则和直接刚度退化模型表征纤维丝和基体的损伤起始与演化,预测了纤维束的弹性性能和强度性能. 其次,将这些性能参数代入介观尺度单胞模型,基于Hashin和Hou的混合失效准则以及连续介质损伤模型对介观尺度单胞进行6种边界条件下的渐进损伤模拟.然后采用渐进均匀化方法,以介观尺度单胞为媒介预测了0$^\circ$和90$^\circ$子胞的性能参数,并建立平纹编织复合材料的子胞模型,进而扩展成为材料的宏观尺度低速冲击模型. 在此基础上,研究了平纹编织复合材料低速冲击下的力学响应与损伤特征.结果表明:宏观冲击仿真和试验吻合较好, 验证了多尺度方法的正确性;最大接触力、材料吸能和分层面积均随冲击能量的增大而增大,分层损伤轮廓逐渐从椭圆形向圆形转化;基体拉伸和压缩损伤的长轴方向分别与子胞材料主方向正交和一致,损伤面积前者远大于后者.      

高温下编织复合材料热相关参数识别方法研究

为了获取高温下编织复合材料的准确弹性参数与热膨胀系数,提出一种基于均匀化理论的热相关参数识别方法. 首先,在编织复合材料单胞有限元模型基础上,基于均匀化理论和热弹性理论,施加周期性位移边界条件和温度边界条件,预测编织复合材 料的热弹性相关参数. 然后,考虑到等效过程中编织复合材料应力分布不均匀等因素引起的误差,将复合材料精细模型的热模态数据作为补 充信息,识别编织复合材料热相关参数,对预测的材料参数进行校准. 本文在二维编织结构单胞模型基础上,开展等效预测和识别方法研 究,验证所提出方法的有效性和准确性. 对比等效和识别后热模态的误差,结果表明:本文提出的基于等效预测的参数识别方法,能够 准确识别高温下编织复合材料宏观热相关参数.      

纺织复合材料和结构多尺度耦合的数值分析

研究了纺织复合材料和结构多尺度耦合的数值分析模型。建立了微、细观单胞,给出了纺织复合材料平均弹性常数的逐级分析方法,着重研究了由宏观结构、到细观纤维束、再到微观纤维三个尺度耦合的应力分析方案。对于常用的板壳状纺织复合材料结构,在面内载荷下,假设每层细观单胞的平均面内应变是一致的,在弯曲、横向剪切及扭曲等非面内载荷下,在内力等效条件下将沿厚度方向连续分布的宏观应力简化为阶梯状分布,忽略了每层细观单胞范围内宏观应力沿厚度方向的梯度变化,由此利用细观单胞模型实现宏观应力与细观应力之间的传递,再利用微观单胞可得到纤维尺度的微观应力。最后以一种三维机织复合材料为例,用上述多尺度耦合的模型逐级分析了材料的平均弹性常数,并沿相反方向,由宏观结构分析逐级计算出纤维束尺度和纤维尺度的细、微观应力的局部波动。     

倾斜内锁型三维机织陶瓷基复合材料的细观力学模型

利用平均化方法提出了倾斜内锁型三维机织陶瓷基复合材料弹性性能分析的三维细观力学模型,对材料的弹性性能进行了预测。这个力学模型考虑了倾斜内锁型三维机织陶瓷基复合材料经向纤维束的弯曲和纬向纤维束的平直,纤维束的横截面形状尺寸和相邻纤维束之间的孔洞以及材料制造过程中碳纤维性能下降对弹性性能的影响。基于层合板理论,提出两种单胞应变状态假设分别对材料的九个弹性常数进行了推导计算,结果表明两种方法理论的预测值非常接近。计算结果与实验值比较吻合,表明所提出的细观力学模型是合理的,可以为纺织陶瓷基复合材料的优化设计提供有价值的参考。     

含大量随机分布细小纤维的单向纤维增强复合材料横向弹性性能参数计算的内嵌区域模型法

利用有限元方法求取单向纤维增强复合材料的横向弹性性能参数的计算模型包括三维模型、两维平面应变模型、单胞模型等等.由于单胞模型仅仅适用于纤维规则排列情况.在纤维随机分布且纤维大小亦为随机时,单向纤维增强复合材料横向弹性性能参数必须通过对于复合材料块体的计算才能获得.同时在随机分布纤维的数量增大时,三维模型和二维平面应变模型的计算量急剧增加,模型的处理能力不强.该文提出一种利用内嵌区域模型来计算含大量随机大小、随机分布细小纤维的单向纤维增强复合材料块体的横向弹性性能参数的方法,有效降低了计算量.在较低的计算费用下,能够快速获得单向纤维增强复合材料的横向弹性性能参数.     

三维编织复合材料圆柱壳的扭转后屈曲分析

随着三维编织复合材料应用的日益广泛,且由于三维编织复合材料的优异性能是与细观结构密切相关的,深入研究其细观结构显得尤为重要.本文通过研究三维圆型编织中所形成的空间纱线交织结构的特征,给出改进的三胞结构单胞模型,通过单胞的组装和变换,得到编织预制件的整体力学性能模型.基于Reddy高阶剪切变形理论导得广义Krmn型大挠度方程,计及边界层效应,采用奇异摄动法,考虑非线性前屈曲、大挠度和初始几何缺陷的影响,给出三维四向编织复合材料圆柱壳在扭转载荷作用下的严格满足边界条件的大挠度渐近解,讨论了纤维体积含量、编织角和几何参数等因素对圆柱壳屈曲和后屈曲行为的影响.     

二维三轴编织复合材料的弹性性能分析

提出了二维三轴编织复合材料的几何模型,模型考虑了纤维束的弯曲扭转状态及空间交错特性等几何元素。基于体积平均法,建立了预测二维三轴编织复合材料弹性性能的理论分析模型;通过引入更普遍适用的周期性位移边界条件,结合二维三轴编织复合材料的细观实体结构,建立了分析其力学性能的有限元模型。两种模型预测结果均与试验结果吻合,证明了方法的合理有效性。分析了材料受载下的细观应力分布,并讨论了编织参数对材料性能的影响。研究表明,二维三轴编织复合材料轴向性能得到了增强,应力分布更均匀,编织角以及纤维体积含量对材料弹性性能影响较大。     

三维机织复合材料的弹性性能预报模型

建立了基于等效响应比拟技术的三维机织复合材料弹性性能预报模型．首先将三维机织物的结构单元分解为4个子元(经纱、纬纱、填充纱和接结纱),用几何模型去估算这些子元的体积分数．然后依据不同的外载形式,将复合材料的应力-应变关系等效地表达为3组诸子元所组成的三维弹簧网络．根据刚度系数的物理意义,采用不同的弹簧网络连接形式,并按体积平均化方法获得材料总体刚度矩阵中相应的刚度系数,进而计算得到三维机织复合材料的9个弹性系数．该模型考虑了层内交织经纱、层间交织接结纱的弯曲以及材料内部纯树脂区对三维机织复合材料弹性性能的影响．试验结果与模型的理论预测值进行比较,表明这个模型是有效的。     

Materially and geometrically non-linear woven composite micro-mechanical model with failure for finite element simulations

A computational micro-mechanical material model of woven fabric composite material is developed to simulate failure. The material model is based on repeated unit cell approach. The fiber reorientation is accounted for in the effective stiffness calculation. Material non-linearity due to the shear stresses in the impregnated yarns and the matrix material is included in the model. Micro-mechanical failure criteria determine the stiffness degradation for the constituent materials. The developed material model with failure is programmed as user-defined sub-routine in the LS-DYNA finite element code with explicit time integration. The code is used to simulate the failure behavior of woven composite structures. The results of finite element simulations are compared with available test results. The model shows good agreement with the experimental results and good computational efficiency required for finite element simulations of woven composite structures.     

Local elastodynamic stresses in the unit cell of a woven fabric composite

Summary A theoretical study of the local elastodynamic stresses of woven fabric composites under dynamic loadings is presented in this article. The analysis focuses on the unit cell of an orthogonal woven fabric composite, which is composed of two sets of mutually orthogonal yarns of either the same fiber (nonhybrid fabric) or different fibers (hybrid fabric) in a matrix material. Using the mosaic model for simplifying woven fabric composites and a shear lag approach to account for the inter-yarn deformation, a one-dimensional analysis has been developed to predict the local elastodynamic and elastostatic behavior. The initial and boundary value problems are formulated and then solved using Laplace transforms. Closed form solutions of the dynamic displacements and stresses in each yarn and the bond shearing stresses at the interfaces between adjacent yarns are obtained in the time domain for any type of in-plane impact loadings. When time tends to infinity, the dynamic solutions approach to their corresponding static solutions, which are also developed in this article. Solutions of certain special cases are identical to those reported in the literature. Lastly, the dynamic stresses and bond shearing stresses of plain weave composites subjected to step uniform impacts are presented and discussed as an example of the general analytical model. Received 3 May 1999; accepted for publication 22 September 1999     

An analytical model to predict residual thermal stress in 2D orthogonal plain weave fabric composites

This paper seeks to address a practical rectangular truss model to predict residual thermal stress in a 2 D plain weave fabric (PWF) composite. The two orthogonal yarns in a micromechanical unit cell are idealized as straight rods subjected to tensile or compression loading resulting in extension or shortening deformation. The residual thermal stresses and equivalent thermal expansion coefficients in a PWF layer are derived from the thermal constitutive equations and the deformation compatibility condition. Based on the deformation compatibility equations, the thermal constitutive relationships for PWF composites are obtained to derive the residual thermal stresses between PWF plies and pure resin. In order to validate the model, experiments have been performed to investigate the mechanical properties of two-dimensional (2D) orthogonal EW220/5284 PWF composites fabricated by resin transfer moulding (RTM). It is shown that the experimental results correlate well with predictions from the new model.     

Equivalent mechanical properties of textile monolayers from discrete asymptotic homogenization

The determination of the effective mechanical moduli of textiles from mechanical measurements is usually difficult due to their discrete architecture, which makes micromechanical analyses a relevant alternative to access those properties. Micropolar continuum models describing the effective mechanical behavior of woven fabric monolayers are constructed from the homogenization of an identified repetitive pattern of the textile within a representative unit cell. The interwoven yarns within the textile are represented as a network of trusses connected by nodes at their crossover points. These trusses have extensional and bending rigidities to allow for yarn stretching and flexion, and a transverse shear deformation is additionally considered. Interactions between yarns at the crossover points are captured by beam segments connecting the nodes. The woven fabric is modeled after homogenization as an anisotropic planar continuum with two preferred material directions in the mean plane of the textile. Based on the developed methodology, the effective mechanical properties of plain weave and twill are evaluated, including their bending moduli and characteristic flexural lengths. A satisfactory agreement is obtained between the effective moduli obtained by homogenization and numerical values obtained by finite element simulations performed over periodic unit cells.     

Equivalent properties of monolayer fabric from mesoscopic modelling strategies

A general and systematic approach for the development of mesostructurally-based continuum model of woven fabrics has been elaborated, relating the fabric behavior at the macroscopic continuum scale to the response and geometry of the fabric’s mesostructure (geometrical configuration of the weave and the yarn properties). Mesoscopic discrete models of dry fabric have been developed based on a discretization of the yarn geometry, accounting for the yarn–yarn interactions at the yarns crossing points. The yarns are modeled within a unit cell consisting of the repetitive fabric pattern as curved planar beams submitted to the reaction forces of the transverse yarns at discrete crossover points. Those reaction forces are expressed in semi-analytical form versus the yarn geometry and mechanical properties for general armour from beam theory. The equilibrium shape of the woven fabric is obtained by minimizing its total potential energy, accounting for the work of the reaction forces due to the transverse yarns. The absolute minimum of the structure’s total potential energy is achieved by a classical genetic algorithm. Simulation results show that plain weave presents a nonlinear response in the early deformation stage due to the crimp change, whereas twill shows a quasi linear response due to yarn extension being the dominant deformation mechanism. Plain weave fabric overall exhibits an orthotropic constitutive law, as biaxial simulations show. The transverse behavior of plain weave fabric is presently evaluated in terms of Poisson’s ratio, based on virtual simulations at the mesoscopic scale of analysis. Simulation results show that Poisson’s ratio first increases towards a maximum due to the rapid shrinkage of the sample in the transverse direction, and decreases thereafter when the crimp changes become limited by the reaction forces of the transverse yarns. The influence of the mechanical properties of both warp and weft on Poisson’s coefficient is assessed. The predictions of the mesoscopic models regarding the impact of yarn geometry and mechanical properties on the overall behavior provide a guideline for the design of woven fabrics.     

Convexity of the Strain-Energy Function in a Two-Scale Model of Ideal Fabrics

A two-scale model is used to generate the macro-scale constitutive response of a sheet of woven fabric from a micro-scale model of interacting yarns regarded as crossed elasticae in contact. The model furnishes a macro-scale strain-energy function for an orthotropic membrane idealized as being weak in shear compared to the extensional resistance of material curves representing the yarns. The operative Legendre–Hadamard inequality for the membrane is derived and shown to be satisfied by a suitably relaxed version of the computed strain-energy function.     

Consideration of the yarn–yarn interactions in meso/macro discrete model of fabric. Part I: Single yarn behaviour

A mesoscopic discrete model of dry fabric has been developed, based on the yarn–yarn interactions occurring at the yarns crossing points. The fabric yarns, described initially by a Fourier series development, are discretized into elastic straight bars represented by stretching springs and connected at frictionless hinges by rotational springs. The motion of each node is described by a lateral displacement and a rotation. The expression of the reaction force exerted by the transverse yarns at the contact points is assessed, from which the work of the reaction forces is established. The equilibrium shape of the yarn is obtained as the minimum of its total potential energy, accounting for the work of the reaction forces due to the transverse yarns. Simulations of a traction curve of a single yarn are performed, that evidence the effect of the yarn interactions. The two principal deformation mechanisms, the variation of undulation and the yarn stretching, are separately analysed.     

Consideration of the yarn–yarn interactions in meso/macro discrete model of fabric: Part II: Woven fabric under uniaxial and biaxial extension

A mesoscopic discrete model of fabric has been developed, accounting for the yarn–yarn interactions occurring at the yarn crossing points. The fabric yarns, described in their initial state by a Fourier series development, are discretized into elastic straight bars represented by stretching springs, and connected at frictionless hinges by rotational springs. In the first part of the paper, the behavior under uniaxial tension of a single yarn has been investigated, and the impact of the interactions of the transverse yarns has been quantitatively assessed. The consideration of the yarn interactions is extended in this second part at the scale of the whole network of interwoven yarns, under uniaxial and biaxial loading conditions. The effect of the transverse yarns properties under uniaxial tension is evidenced, as well as the impact of the biaxial loading ratio.     

NONLINEAR MICRO－MECHANICAL MODEL FOR PLAIN WOVEN FABRIC

The warp yarns and weft yarns of plain woven fabric which, being the principal axes of material of fabric, are orthogonal in the original configuration, but are obliquely crossed in the deformed configuration in general. The orthotropic constitutive model is unsuitable for fabric. In the oblique principal axes system the relations between loaded stress vectors and stress tensor are investigated, the stress fields of micro-weaving structures of fabric due to pure shear are carefully studied and, finally, a nonlinear micro-mechanical model for plain woven fabric is proposed. This model can accurately describe the nonlinear mechanical behavior of fabric observed in experiments. Under the assumption of small deformation and linearity of mechanical properties of fabric the model will degenerate into the existing linear model.