Structure of the constitutive relations for hereditarily elastic materials reinforced by hard fibers

The problem of simplifying the nonlinear hereditary elasticity relations is considered for strongly anisotropic materials such as fiber-reinforced composites. This is done using their property that the material stiffness is high along the reinforcement and is low in the cross-sectional direction. The material is considered transversally isotropic. The simplification is performed by analyzing asymptotic representations of creep relations. Relations of various degrees of accuracy for various types of composites and stress states are obtained. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 3, pp. 120–127, May–June, 2005.     

宏微观渐进展开损伤本构模型及其在碳纤维增强复合材料中的应用

基体开裂、纤维拔出、界面剥离等是碳纤维增强复合材料常出现的局部各向异性损伤现象,这些损伤逐渐扩展,削弱了材料的强度和刚度,影响材料的承载能力.对此利用宏微观摄动理论对位移进行双范围渐进展开,在微观位移中引入损伤应变,通过计算损伤应变集中因子,得到了含损伤的均质化损伤弹性常数(宏观有效刚度矩阵),用平均法和混合法检验了无...     

Mesoscale modeling of dynamic compression of boron carbide polycrystals

An anisotropic nonlinear elastic model is advanced for crystals belonging to either of two polytypes of boron carbide ceramic. Crystals undergo transformation to an isotropic, amorphous phase upon attainment of a local state-based criterion associated with a loss of intrinsic stability. The model is implemented using the dynamic finite element method, and is demonstrated on a representative volume consisting of fifty polyhedral grains subjected to uniaxial strain at a uniform high strain rate and shock compression at axial pressures ranging from 10 to 50 GPa. Predicted stress–strain behavior is in close agreement with experimental data. For polycrystals consisting of both polytypes, amorphization initiates at stress levels slightly below the experimental Hugoniot elastic limit, and occurs more readily than observed in experiment. For polycrystals consisting only of the CBC (polar) polytype, amorphization initiates at impact pressures similar to those suggested by experiment. In either case, transformation is promoted by dynamic stress interactions and elastic coefficient mismatch among anisotropic crystals. Results support a previous conjecture that amorphization is related to shear instability and cross-linking of the CBC chain in the polar polytype.     

弱界面复合材料中的基体裂纹

利用二维弹性力学模型研究了纤维增强复合材料中基体裂纹与弱界面的相互作用机理．文中首先导出各向异性弹性多层介质中刃型位错的基本解，然后运用这些基本解建立了弱界面复合材料中典型的Ｈ型缺陷的奇异积分方程组，通过求解这些方程得到外载荷的大小、弱界面的结合强度、界面的残余压力和摩擦系数、纤维与基体的弹性模量比等微结构参量与基体裂纹附近的应力场的关系     

A failure criterion for fiber reinforced polymer composites under combined compression–torsion loading

A fracture mechanics based failure criterion for unidirectional composites under combined loading has been developed. The predictions from this criterion have been compared with experimental data obtained from combined compression–torsion loading of glass and carbon fiber reinforced polymer composites of 50% fiber volume fraction. The specimens were loaded under rotation control and displacement control in a proportional manner. Comparison of the Budiansky–Fleck kinking model, specialized to a solid circular cylinder, and the new failure model against experimental data suggests that the Budiansky–Fleck model predictions do not capture the variation of compressive strength as a function of shear stress for glass fiber composites. This is because these composites fail predominantly by compressive splitting. The Budiansky–Fleck model predictions are appropriate for composites that fail by compressive kinking. The new model predictions capture the experimental results for glass composites where the compression strength is initially unaffected by shear stress but undergoes a drastic reduction when a critical value of shear stress is reached.     

PP/CF增强珊瑚砂水泥基复合材料冲击压缩力学性能研究

在人工海水制备珊瑚砂水泥基复合材料中混杂加入碳纤维和聚丙烯纤维,得到4种不同纤维掺量的碳-聚丙烯混杂纤维增强珊瑚砂水泥基复合材料。采用直径100 mm的分离式Hopkinson压杆,对材料进行5种应变率下的冲击压缩试验,采用LS-DYNA进行相应的冲击压缩数值模拟。结果表明:(1) 试验应变率临界值为200 s?1,当试验应变率大于200 s?1时,混杂碳纤维和聚丙烯纤维所形成的纤维网络对试块的增韧效果加强;(2) 碳-聚丙烯混杂纤维增强珊瑚砂水泥基复合材料峰值应力具有明显的应变率效应,且动态增强因子对应变率的敏感程度较高;(3) 使用珊瑚砂细骨料导致试块内微裂纹和微空洞等缺陷较多,在珊瑚砂水泥基复合材料内混杂掺加碳纤维和聚丙烯纤维后,试块冲击抗压强度的提升有限,但珊瑚砂水泥基复合材料的抗冲击韧性显著提升;(4) 通过试验数据和参数调试确定了HJC模型的参数,试块峰值应力的模拟结果与试验结果的误差在5.97 %以内。     

Estimation of the effective thermoelastic moduli of fibrous nanocomposites with cylindrically anisotropic phases

Recent developments in nanotechnology make it possible to fabricate nanofibers and identify their mechanical fibers. In particular, nanofibers are used as reinforcement in composites. The present work concerns unidirectional nanofibrous composites with cylindrically anisotropic phases and aims to analytically estimate their effective thermoelastic moduli. This objective is achieved by extending the classical generalized self-consistent model to the setting of thermoelasticity, to the case of cylindrically anisotropic phases, and to the incorporation of interface stress effect. Analytical closed-form estimations are derived for all the effective thermoelastic moduli, showing that these moduli depend on the fiber cross-section size. Numerical examples are provided to illustrate this size-dependent effect.     

纤维增强陶瓷基复合材料的D失效判据

经典唯象强度理论适用于正交各向异性线弹性体。对于非线性纤维增强复合材料,通过加卸载试验和损伤力学的分析方法,可以得到一种虚拟的线性化应力-应变关系;依据损伤等效假设,针对线性损伤和非线性损伤,对基于应力的经典二次失效准则进行变换,建立了一种基于损伤的强度理论,即“D失效判据”,这一强度理论可以作为经典判据的补充和扩展。针对平纹编织C/SiC复合材料的拉/剪组合试验,进行了实例计算,结果表明：利用D失效判据预测的失效包络线比蔡-希尔准则的预测曲线低,而且,失效曲线的形式与材料的损伤演化规律相关。     

Macroscopic behavior and field fluctuations in viscoplastic composites: Second-order estimates versus full-field simulations

This work presents a combined numerical and theoretical study of the effective behavior and statistics of the local fields in random viscoplastic composites. The full-field numerical simulations are based on the fast Fourier transform (FFT) algorithm [Moulinec, H., Suquet, P., 1994. A fast numerical method for computing the linear and nonlinear properties of composites. C. R. Acad. Sci. Paris II 318, 1417-1423], while the theoretical estimates follow from the so-called “second-order” procedure [Ponte Castañeda, P., 2002a. Second-order homogenization estimates for nonlinear composites incorporating field fluctuations: I—Theory. J. Mech. Phys. Solids 50, 737-757]. Two-phase fiber composites with power-law phases are considered in detail, for two different heterogeneity contrasts corresponding to fiber-reinforced and fiber-weakened composites. Both the FFT simulations and the corresponding “second-order” estimates show that the strain-rate fluctuations in these systems increase significantly, becoming progressively more anisotropic, with increasing nonlinearity. In fact, the strain-rate fluctuations tend to become unbounded in the limiting case of ideally plastic composites. This phenomenon is shown to correspond to the localization of the strain field into bands running through the composite along certain preferred orientations determined by the loading conditions. The bands tend to avoid the fibers when they are stronger than the matrix, and to pass through the fibers when they are weaker than the matrix. In general, the “second-order” estimates are found to be in good agreement with the FFT simulations, even for high nonlinearities, and they improve, often in qualitative terms, on earlier nonlinear homogenization estimates. Thus, it is demonstrated that the “second-order” method can be used to extract accurate information not only for the macroscopic behavior, but also for the anisotropic distribution of the local fields in nonlinear composites.     

A periodic unit-cell simulation of fiber arrangement dependence on the transverse tensile failure in unidirectional carbon fiber reinforced composites

The effect of fiber arrangement on transverse tensile failure in unidirectional carbon fiber reinforced composites with a strong fiber-matrix interface was studied using a unit-cell model that includes a continuum damage mechanics model. The simulated results indicated that tensile strength is lower when neighboring fibers are arrayed parallel to the loading direction than with other fiber arrangements. A shear band occurs between neighboring fibers, and the damage in the matrix propagates around the shear band when the interfacial normal stress (INS) is sufficiently high. Moreover, based on the observation of Hobbiebrunken et al., we reproduced the damage process in actual composites with a nonuniform fiber arrangement. The simulated results clarified that the region where neighboring fibers are arrayed parallel to the loading direction becomes the origin of the transverse failure in the composites. The cracking sites observed in the simulation are consistent with experimental results. Therefore, the matrix damage in the region where the fiber is arrayed parallel to the loading direction is a key factor in understanding transverse failure in unidirectional carbon fiber reinforced composites with a strong fiber/matrix interface.     

界面脱粘对陶瓷基复合材料疲劳迟滞回线的影响

采用细观力学方法对脆性纤维增强的陶瓷基复合材料拉-拉疲劳载荷下应力-应变迟滞回线进行了研究,将拉梅公式与库仑摩擦法则相结合分析了界面脱粘区以及粘结区复合材料细观应力场.根据卸载与重新加载时纤维相对基体滑移机制,分析了加卸载纤维轴向应力分布,结合断裂力学界面脱粘准则确定了初始加载界面脱粘长度ls、卸载界面反向滑移长度y以及重新加载界面滑移长度z',讨论了界面脱粘能和界面摩擦系数对初始界面脱粘、卸载界面反向滑移、重新加载界面滑移以及加卸载迟滞回线的影响.并与Pryce-Smith模型和试验数据进行对比表明:该文模型与试验曲线吻合的较好.     

SMA短纤维增强复合材料在热载下的相交特性

在形状记忆合金（ＳＭＡ）复合材料研究中,相变特性的研究是一个主要的工作．基于Ｅｓｈｅｌｂｙ的等效夹杂模型和Ｍｏｒｉ和Ｔａｎａｋａ的场平均法,考虑到ＳＭＡ材料的强物理非线性,发展了增量型的等效夹杂模型（ＩｎｃｒｅｍｅｎｔａｌＥｑｕｉｖａｌｅｎｔＩｎｃｌｕｓｉｏｎＭｏｄｅｌ）．考虑在某一温度循环条件下讨论形状记忆合金短纤维增强的铝基复合材料在热载下的相变行为．特别研究了ＳＭＡ短纤维复合材料在变温过程中纤维几何尺寸、体积分数等参数对ＳＭＡ复合材料的相变行为和ＳＭＡ内残余应力等的影响．这些工作对于指导材料设计和了解ＳＭＡ复合材料热机械特性是颇有意义的．     

采用瓣叶复合材料的生物心脏瓣膜的力学分析

采用微观组织结构分析及宏观复合材料分析结合的方法,分析了猪主动脉瓣的非线性复合材料性质,提出了一种适用于猪主动脉瓣的非线性复合材料本构模型,用提出的非线性复合材料本构模型,对闭合承载状态下的等厚度与变厚度几何模型的猪主动脉瓣的应力分布及变形进行了有限元数值模拟,发现：与各向同性瓣叶相比,单向增强复合材料的瓣叶不但具有较强的承载能力,而且具有较大的柔软性。     

界面特性对短纤维金属基复合材料蠕变行为的影响

基于短纤维增强金属基复合材料(MMC)的单纤维三维模型(三相)，利用粘弹性有限元分析方法对影响金属基复合材料的蠕变行为的因素进行了较为系统的分析。研究中主要讨论了界面特性和纤维取向角对金属基复合材料的蠕变性能的影响。研究结果发现，界面特性诸如厚度、模量和应力指数都对纤维最大轴应力和稳定蠕变率产生影响：稳态蠕变率随界面模量的增大而逐渐减小，当高于基体模量时基本保持不变；纤维轴应力的变化与蠕变率正好相反。稳态蠕变率随界面厚度、应力指数的增加而增大；而轴应力则随之减小。同时不同的纤维取向也影响金属基复合材料蠕变时的轴应力分布和稳态蠕变率。     

Piezoresistive fiber-reinforced composites: A coupled nonlinear micromechanical–microelectrical modeling approach

Piezoresistive composites are materials that exhibit spatial and effective electrical resistivity changes as a result of local mechanical deformations in their constituents. These materials have a wide array of applications from non-destructive evaluation to sensor technology. We propose a new coupled nonlinear micromechanical-microelectrical modeling framework for periodic heterogeneous media. These proposed micro-models enable the prediction of the effective piezoresistive properties along with the corresponding spatial distributions of local mechanical–electrical fields, such as stress, strain, current densities, and electrical potentials. To this end, the high fidelity generalized method of cells (HFGMC), originally developed for micromechanical analysis of composites, is extended for the micro-electrical modeling in order to predict their spatial field distributions and effective electrical properties. In both cases, the local displacement vector and electrical potential are expanded using quadratic polynomials in each subvolume (subcell). The equilibrium and charge conservations are satisfied in an average volumetric fashion. In addition, the continuity and periodicity of the displacements, tractions, electrical potential, and current are satisfied at the subcell interfaces on an average basis. Next, a one way coupling is established between the nonlinear mechanical and electrical effects, whereby the mechanical deformations affect the electrical conductivity in the fiber and/or matrix constituents. Incremental and total formulations are used to arrive at the proper nonlinear solution of the governing equations. The micro-electrical HFGMC is first verified by comparing the stand-alone electrical solution predictions with the finite element method for different doubly periodic composites. Next, the coupled HFGMC is calibrated and experimentally verified in order to examine the effective piezoresistivity of different composites. These include conductive polymeric matrices doped with carbon nano-tubes or particles. One advantage of the proposed nonlinear coupled micro-models is its ability to predict the local and effective electro-mechanical behaviors of multi-phase periodic composites with different conductive phases.     

Modeling the anisotropic finite-deformation viscoelastic behavior of soft fiber-reinforced composites

This paper presents constitutive models for the anisotropic, finite-deformation viscoelastic behavior of soft fiber-reinforced composites. An essential assumption of the models is that both the fiber reinforcements and matrix can exhibit distinct time-dependent behavior. As such, the constitutive formulation attributes a different viscous stretch measure and free energy density to the matrix and fiber phases. Separate flow rules are specified for the matrix and the individual fiber families. The flow rules for the fiber families then are combined to give an anisotropic flow rule for the fiber phase. This is in contrast to many current inelastic models for soft fiber-reinforced composites which specify evolution equations directly at the composite level. The approach presented here allows key model parameters of the composite to be related to the properties of the matrix and fiber constituents and to the fiber arrangement. An efficient algorithm is developed for the implementation of the constitutive models in a finite-element framework, and examples are presented examining the effects of the viscoelastic behavior of the matrix and fiber phases on the time-dependent response of the composite.     

Microstructure and mechanical properties of three-dimensional five-directional braided composites

Three-dimensional (3D) five-directional braided composites are significant structural materials in the fields of astronauts and aeronautics. On the basis of the 3D five-directional braiding process, three types of microstructural unit cell models are established with respect to the interior, surface and corner regions. The mathematical relationships among the structural parameters, such as fiber orientation, fiber volume fraction, the yarn packing factor, are derived. By using these three unit cell models, a micromechanical prediction procedure is described to simulate the stiffness and strength properties of 3D five-directional braided composites. Only the in situ constituent fiber and matrix properties of the composites and the fiber volume proportion are required in the simulation. The stress states generated in the constituent fiber and matrix materials are explicitly correlated with the overall applied load on the composites. The predictive stiffness and strength are in good agreement with available experimental data, which demonstrates the applicability of the present analytical model.     

碳纤维增强复合材料微观烧蚀行为数值模拟

碳纤维增强复合材料已广泛应用于烧蚀热防护系统中,其微观结构直接影响材料的烧蚀行为,因此材料微观烧蚀机制分析对材料设计和制备具有重要意义.本文采用有限体积法并引入固体相体积分数与界面重构技术,建立了碳纤维增强复合材料微观烧蚀数值模型.利用该微观烧蚀数值模型对碳基体包裹单根碳纤维进行烧蚀模拟,将数值模拟的烧蚀形貌与解析解结果进行对比,验证了数值求解方法的正确性.对单向碳纤维增强复合材料在不同碳纤维倾斜角下的微观烧蚀行为进行了模拟分析,得到了碳纤维倾斜角对微观烧蚀行为的影响规律.研究发现:对于碳纤维的抗氧化性能比基体强的情况,当氧扩散速率远远大于碳氧反应速率时,碳纤维将出现"笋尖"状的烧蚀形貌;当碳氧反应速率远远大于氧扩散速率时,纤维和基体将以相同速率烧蚀.当碳纤维的抗氧化性能比基体强时,纤维倾斜角会对材料微观烧蚀行为产生较大影响;相反,则影响不显著.      

A variational method for unidirectional fiber-reinforced composites with matrix creep

A variational method is developed for analyzing the matrix creep induced time-dependent change in fiber stress profiles in unidirectional composites. A functional of admissible profiles of fiber stress rate is presented by supposing a fiber broken in matrix as well as a fiber pulled out from matrix. The functional is shown to have the stationary function satisfying an incremental differential equation based on the shear lag assumption. Then, the stationary function is approximately determined by assuming bilinear profiles of fiber stress and a power law of matrix creep, leading to analytical solutions for the time-dependent change in fiber stress profiles. The solutions are verified on the basis of an energy balance equation and a finite difference computation. Moreover, it is shown that the solution for the fiber pull-out model agrees well with an experiment on a single carbon fiber/acrylic model composite if the initial slip at fiber/matrix interface is taken into account. In addition, the solution for the fiber breakage model is used for evaluating the characteristic time in long-term creep rupture of unidirectional composite.     

改进的Jones—Nelson物理非线性材料模型与橡胶复合材料层板分析

Ｊｏｎｅｓ－Ｎｅｌｓｏｎ模型是复合材料物理非线性应力应变关系的一种描述方法，它通过建立材料刚度与应变能密度的关系以及考虑纤维增强复合材料的物理非线性问题，非线性的材料矩阵只为应变能的子数，使得材料模型可以方便地应用地复杂应力状态下，通过 在大变形的应变应力变量和变泊松比概念，对Ｊｏｎｅｓ－Ｎｅｌｓｏｎ模型进行了改进，解决了材料特性的扩展问题和收敛问题，同时考虑了大变形中纤维铺设角度的重新取向，使其