Viscoelastic response of a unidirectional composite containing two viscoelastic constituents

Viscoelastic response of a unidirectional aramid fiber-reinforced epoxy was measured. Procedures to measure all five time-dependent material properties necessary to describe behavior of a transversely isotropic continuous-fiber unidirectional lamina were implemented. The Iosipescu shear method was used to measure in-plane and interlaminar shear viscoelastic response. Applicability of the Schapery single-integral nonlinear viscoelastic constitutive model to describe time-dependent mechanical behavior of a laminated composite material containing two viscoelastic phases was explored. Linear and nonlinear viscoelastic parameters for this two-viscoelastic constituent composite were measured and data summaries are presented in the paper. The time-dependent behavior of this two-viscoelastic constituent composite material was found to be complex, but the Schapery nonlinear viscoelastic model did adequately fit the response of such a composite to uniaxial applied loads.Paper was presented at the 1988 SEM Spring Conference on Experimental Mechanics held in Portland, OR on June 5–10.     

Neo-Hookean fiber-reinforced composites in finite elasticity

The response of a transversely isotropic fiber-reinforced composite made out of two incompressible neo-Hookean phases undergoing finite deformations is considered. An expression for the effective energy-density function of the composite in terms of the properties of the phases and their spatial distribution is developed. For the out-of-plane shear and extension modes this expression is based on an exact solution for the class of composite cylinder assemblages. To account for the in-plane shear mode we incorporate an exact result that was recently obtained for a special class of transversely isotropic composites. In the limit of small deformation elasticity the expression for the effective behavior agrees with the well-known Hashin-Shtrikman bounds. The predictions of the proposed constitutive model are compared with corresponding numerical simulation of a composite with a hexagonal unit cell. It is demonstrated that the proposed model accurately captures the overall response of the periodic composite under any general loading modes.     

A composites-based hyperelastic constitutive model for soft tissue with application to the human annulus fibrosus

This paper presents a composites-based hyperelastic constitutive model for soft tissue. Well organized soft tissue is treated as a composite in which the matrix material is embedded with a single family of aligned fibers. The fiber is modeled as a generalized neo-Hookean material in which the stiffness depends on fiber stretch. The deformation gradient is decomposed multiplicatively into two parts: a uniaxial deformation along the fiber direction and a subsequent shear deformation. This permits the fiber-matrix interaction caused by inhomogeneous deformation to be estimated by using effective properties from conventional composites theory based on small strain linear elasticity and suitably generalized to the present large deformation case. A transversely isotropic hyperelastic model is proposed to describe the mechanical behavior of fiber-reinforced soft tissue. This model is then applied to the human annulus fibrosus. Because of the layered anatomical structure of the annulus fibrosus, an orthotropic hyperelastic model of the annulus fibrosus is developed. Simulations show that the model reproduces the stress-strain response of the human annulus fibrosus accurately. We also show that the expression for the fiber-matrix shear interaction energy used in a previous phenomenological model is compatible with that derived in the present paper.     

Fiber push-out study of a copper matrix composite with an engineered interface: Experiments and cohesive element simulation

The fiber push-out test is a basic method to probe the mechanical properties of the fiber/matrix interface of fiber-reinforced metal matrix composites. In order to estimate the interfacial properties, parameters should be calibrated using the measured load–displacement data and theoretical models. In the case of a soft matrix composite, the possible plastic yield of the matrix has to be considered for the calibration. Since the conventional shear lag models are based on elastic behavior, a detailed assessment of the plastic effect is needed for accurate calibration. In this paper, experimental and simulation studies are presented regarding the effect of matrix plasticity on the push-out behavior of a copper matrix composite with strong interface bonding. Microscopic images exhibited significant local plastic deformation near the fibers leading to salient nonlinear response in the global load–displacement curve. For comparison, uncoated interface with no chemical bonding was also examined where the nonlinearity was not observed. A progressive FEM simulation was conducted for a complete push-out process using the cohesive zone model and inverse fitting. Excellent coincidence was achieved with the measured push-out curve. The predicted results confirmed the experimental observations.     

形状记忆合金纤维复合材料的等效力学行为

在Aboudi提出的胞元模型以及Liu等建立的形状记忆合金的本构模型的基础上，由Legendre多项式，假设每个子胞元的位移场、应变场和应力场，再由子胞元间交界面的应力连续条件和外荷载边界条件推导出基体为弹塑性材料的形状记忆合金纤维复合材料的胞元模型；模拟了呈周期对称的形状记忆合金纤维复合材料受轴向单向拉伸、横向拉伸和横向剪切荷载作用下的等效力学行为，与有限元解进行了比较，结果基本一致。与有限元法比较起来，本文推导出的形状记忆合金纤维复合材料的胞元模型更具高效性。     

Homogenization of long fiber reinforced composites including fiber bending effects

This paper presents a homogenization method, which accounts for intrinsic size effects related to the fiber diameter in long fiber reinforced composite materials with two independent constitutive models for the matrix and fiber materials. A new choice of internal kinematic variables allows to maintain the kinematics of the two material phases independent from the assumed constitutive models, so that stress–deformation relationships, can be expressed in the framework of hyper-elasticity and hyper-elastoplasticity for the fiber and the matrix materials respectively. The bending stiffness of the reinforcing fibers is captured by higher order strain terms, resulting in an accurate representation of the micro-mechanical behavior of the composite. Numerical examples show that the accuracy of the proposed model is very close to a non-homogenized finite-element model with an explicit discretization of the matrix and the fibers.     

A multi-scale framework for layered composites with thermo-rheologically complex behaviors

A multi-scale modeling framework is proposed for generating the effective nonlinear thermo-viscoelastic responses of multi-layered and thick-section composites. The modeling framework is demonstrated for multi-layered composite systems consisting of two alternating fiber-reinforced polymeric layers: unidirectional fiber (roving) and continuous filament mat (CFM). The viscoelastic behavior of the polymeric matrix is considered as stress-dependent and thermo-rheologically complex. Two simplified 3D micromechanical models with periodic boundary conditions are formulated for the roving and CFM. In addition, a sublaminate model is developed for the 3D effective homogenized through-thickness response of the roving and CFM layers. This framework provides an effective time-dependent material response while simultaneously recognizing the corresponding deformation at the microconstituents under overall thermo-mechanical loadings at the macro-level. Stress correction algorithms are developed at each scale to enhance computational efficiency and accuracy. Short-term (30 min) creep tests on off-axis multi-layered specimens are also conducted under combined stresses and temperatures to calibrate in-situ fiber and matrix properties and verify the predictions of the multi-scale framework. A time-shifting method is applied to create long-term material behaviors from the available short-term creep data. Verification of the multi-scale material framework is also done for the overall long-term responses. Finally, long-term structural analyses under thermo-mechanical loadings are conducted by integrating the multi-scale material framework to finite element (FE) structural analyses.     

On thermodynamically consistent constitutive equations for fiber-reinforced nonlinearly viscoelastic solids with application to biomechanics

In this paper, we are interested in developing thermodynamically consistent constitutive equations for fiber-reinforced nonlinearly viscoelastic bodies, in particular for transversely isotropic nonlinearly viscoelastic solids, in isothermal processes. It follows from results in the theory of algebraic invariants that constitutive equations for such materials can be expressed in terms of functions of 18 independent invariants associated with deformation and fiber orientation: 10 of them are isotropic invariants and 8 of them are associated with the deformation and the orientation of the fiber. Among the 8 anisotropic invariants just 6 are related to the viscoelastic response. The terms in the Cauchy stress tensor associated to these 6 invariants are analyzed with respect to thermodynamical consistency, and we obtain restrictions for the corresponding constitutive coefficients. This framework is applied to viscoelastic potentials within the context of biomaterials.     

兼具基体蠕变和轻度界面蠕变的复合材料的本构特性

金属基复合材料的基体蠕变特性在一定的温度与应力下表现明显，当界面为非理想并具有粘弹性性质时，其对复合材料的整体蠕变亦有重要影响．本文利用具有非理想界面的复合材料的Ｍｏｒｉ－Ｔａｎａｋａ方法，研究了既具有基体蠕变又具有界面轻度蠕变的复合材料的本构关系，给出了各微结构参量对复合材料整体蠕变特性的影响     

兼具基体蠕变和轻度界面蠕变的复合材料的本构特性

金属基复合材料的基体蠕变特性在一定的温度与应力下表现明显，当界面为非理想并具有粘弹性性质时，其对复合材料的整体蠕变亦有重要影响．本文利用具有非理想界面的复合材料的Ｍｏｒｉ－Ｔａｎａｋａ方法，研究了既具有基体蠕变又具有界面轻度蠕变的复合材料的本构关系，给出了各微结构参量对复合材料整体蠕变特性的影响     

A general hyperelastic model for incompressible fiber-reinforced elastomers

This work presents a new constitutive model for the effective response of fiber-reinforced elastomers at finite strains. The matrix and fiber phases are assumed to be incompressible, isotropic, hyperelastic solids. Furthermore, the fibers are taken to be perfectly aligned and distributed randomly and isotropically in the transverse plane, leading to overall transversely isotropic behavior for the composite. The model is derived by means of the “second-order” homogenization theory, which makes use of suitably designed variational principles utilizing the idea of a “linear comparison composite.” Compared to other constitutive models that have been proposed thus far for this class of materials, the present model has the distinguishing feature that it allows consideration of behaviors for the constituent phases that are more general than Neo-Hookean, while still being able to account directly for the shape, orientation, and distribution of the fibers. In addition, the proposed model has the merit that it recovers a known exact solution for the special case of incompressible Neo-Hookean phases, as well as some other known exact solutions for more general constituents under special loading conditions.     

Modeling of fiber-reinforced cement composites: Discrete representation of fiber pullout

The discrete modeling of individual fibers in cement-based materials provides several advantages, including the ability to simulate the effects of fiber dispersion on pre- and post-cracking composite performance. Recent efforts in this direction have sought a balance between accurate representation of fiber behavior and computational expense. This paper describes a computationally efficient approach to representing individual fibers, and their composite behavior, within lattice models of cement-based materials. Distinguishing features of this semi-discrete approach include: (1) fibers can be positioned freely in the computational domain, irrespective of the background lattice representing the matrix phase; (2) the pre- and post-cracking actions of the fibers are simulated with little computational expense, since the number of system degrees of freedom is independent of fiber count. Simulated pullouts of single fibers are compared with theory and test results for the cases of perfectly-plastic and slip-hardening behavior of the fiber–matrix interface. To achieve objective results with respect to discretization of the matrix, pullout forces are distributed along the embedded lengths of fibers that bridge a developing crack. This is in contrast to models that lump the pullout force at the crack surfaces, which can lead to spurious break-off of matrix particles as the discretization of the matrix is refined. With respect to fracture in multi-fiber composites, the proposed model matches theoretical predictions of post-cracking strength and pullout displacement corresponding to the load-free condition. The work presented herein is a significant step toward the modeling of strain-hardening composites that exhibit multiple cracking.     

Creep-rupture lifetime simulation of unidirectional metal matrix composites with and without time-dependent fiber breakage

A numerical simulation for predicting the axial creep-rupture lifetime of continuous fiber-reinforced metal matrix composites is proposed, based on the finite element method. The simulation model is composed of line elements representing the fibers and four-node isoparametric plane elements representing the matrix. While the fibers behave as an elastic body at all times, the matrix behaves as an elasto-plastic body at the loading process and an elasto-plastic creep body at the creep process. It is further assumed in the simulation that the fibers are fractured not only in stress criterion but time-dependently with random nature. Simulation results were compared with the creep-rupture lifetime data of a boron-aluminum composite with 10% fiber volume fraction experimentally obtained. The simulated creep-rupture lifetimes agreed well with the averages of the experimental data. The proposed simulation is further carried out to predict a possibility of creep-rupture for the composite without time-dependent fiber breakage. It is finally concluded that the creep-rupture of a boron-aluminum composite is closely related with the shear stress relaxation occurring in the matrix as well as time-dependent fiber breakage.     

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

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

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

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

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.     

The averaged mechanical properties of prismatic lattice formed by ellipsoidal inclusions

Summary The objective of this paper is to evaluate the averaged elastic properties of 3-D grained composites in which identical inclusions form a prismatic network interacting with the matrix material. The inclusions are of ellipsoidal shape with transverse circular sections located at the nodes of a doubly-periodic lattice with an orthogonal elementary cell. When the arrays of inclusions are set at equal spacings in normal directions through the thickness of the matrix, the material formed is an anisotropic composite with tetragonal symmetry at planes transverse to the fiber axis. The longitudinal and transverse elastic and shear moduli as well as the longitudinal Poisson's ratios of such composites are evaluated in this paper. The averaged properties are studied in terms of the aspect ratio and volume fraction of the inclusions as well as the relative rigidity of the constituent phases. Employing the Eshelby's theory for the stress field around a single ellipsoidal inhomogeneity, which is surrounded by the effective anisotropic material, and considering the Mori-Tanaka's concept for the mutual interaction of the neighboring inclusions, we may evaluate the averaged elastic properties of grained composites with aligned ellipsoidal inclusions at finite concentrations. The results provided in a closed-form solution concern the stiffness of 3-D grained composites with parallely dispersed ellipsoidal inclusions forming a prismatic network inside the principal material. It is shown that the stiffness is affected by both the geometry of the inclusions and their concentration. The use of different composite models in the analysis shows that intense variations of stiffness occur mainly in hard composites weakened by soft ellipsoidal inclusions. These findings come in full verification with experimental or theoretical results from the literature. Received 10 February 1998; accepted for publication 27 November 1998     

Constitutive relation of discontinuous reinforced metal-matrix composites

A micromechanical model is developed to simulate the mechanical behaviors of discontinuous reinforced composites. The analysis for a representative unit cell is based on the assumption of a periodic array of aligned reinforcements. The minimum energy principle is used to determine the unknown coefficients of the displacement field of the unit cell. The constitutive behavior of composites is studied to obtain the relationship between the main variables of matrix and reinforcements. It is concluded that the flow strength of composites is strongly influenced by volume fraction, aspect ratio of reinforcement, and the strain hardening exponent of matrix. An analytical constitutive relation of composites is obtained. The predicted results are in agreement with the existing experimental and numerical results. The project supported by the National Natural Science Foundation of China (19704100) and National Science Foundation of Chinese Academy of Sciences (KJ951-1-20)     

Onset of failure in finitely strained layered composites subjected to combined normal and shear loading

A limiting factor in the design of fiber-reinforced composites is their failure under axial compression along the fiber direction. These critical axial stresses are significantly reduced in the presence of shear stresses. This investigation is motivated by the desire to study the onset of failure in fiber-reinforced composites under arbitrary multi-axial loading and in the absence of the experimentally inevitable imperfections and finite boundaries.By using a finite strain continuum mechanics formulation for the bifurcation (buckling) problem of a rate-independent, perfectly periodic (layered) solid of infinite extent, we are able to study the influence of load orientation, material properties and fiber volume fraction on the onset of instability in fiber-reinforced composites. Two applications of the general theory are presented in detail, one for a finitely strained elastic rubber composite and another for a graphite-epoxy composite, whose constitutive properties have been determined experimentally. For the latter case, extensive comparisons are made between the predictions of our general theory and the available experimental results as well as to the existing approximate structural theories. It is found that the load orientation, material properties and fiber volume fraction have substantial effects on the onset of failure stresses as well as on the type of the corresponding mode (local or global).     

纤维增强复合材料弹性性能预测的域分解方法及应用

提出了新的有限元建模方法,即域分解方法,用于预测纤维增强复合材料单向带T300/BSL914C(环氧树脂)和AS4/3501-6(环氧树脂)的弹性性能。域分解方法基于区域叠合技术,分别建立单胞的整体域与纤维域模型用于代替传统有限元建模方法中单胞的基体域与纤维域模型。整体域是真实基体体积与纤维体积的叠加,两区域网格独立划分,互不影响。采用MSC.Nastran中的多节点约束Explicit单元,在整体域与纤维域节点之间建立位移连接属性模拟单胞基体域与纤维域之间的位移约束关系,从而实现两区域的耦合计算。计算结果表明:域分解方法单胞模型纤维增强方向弹性模量Ez预测值与试验值误差在7%以内,其余弹性常数也都与试验值吻合较好。域分解方法不仅可以大大简化纤维增强复合材料的细观力学建模,而且可以准确地预测纤维增强复合材料的弹性性能。