Development of constitutive laws for thermomechanical behaviors of composites containing multi-type ellipsoidal reinforcements

Composite materials are widely used in industrial applications because of their excellent properties and behaviors. While a composite material is defined as a mixture of two or more different materials, many research works in the literature dealt with composites of only two constituents, which are matrix and one type of particles. On the other hand, the theoretical research works that dealt with more than two constituents are rare. Using some additives affects the sintering behavior, the tribological behavior and the fracture mechanics behavior of composites. For example, a suitable amount of additives as sintering aids, in the sintering process, could lower the sintering temperature, enhance phase wettability and bonding strength and improve the interlaminar fracture resistance of a composite. Therefore, it is worthwhile to develop the constitutive laws that describe the behavior of such composite materials. Accordingly, the aim of this paper is to modify the previous paper, Shabana (2003) [Shabana, Y.M., 2003. Incremental constitutive equation for discontinuously reinforced composites considering reinforcement damage and thermoelastoplasticity. Computational Materials Science 28, 31–40], in order to propose constitutive laws that predict the thermomechanical behavior of composites containing multi-type ellipsoidal reinforcements. This includes reinforcements with different materials and/or different shapes that are represented by aspect ratios. These constitutive laws not only predict the macroscopic and microscopic thermoelastoplastic behaviors of composites containing multi-type ellipsoidal reinforcements, but also characterize their different overall effective properties such as modulus of elasticity, Poison’s ratio and thermal expansion coefficient in different directions. Beside this, they are applicable for porous materials and composites with multiple reinforcements and porosities of different shapes and distributions. In the present numerical analyses, composites with two, three and four constituents considering different materials and aspect ratios as well as reinforcement damage are discussed.     

Plastic constitutive behavior of short-fiber/particle reinforced composites

A cylindrical cell model based on continuum theory for plastic constitutive behavior of short-fiber/particle reinforced composites is proposed. The composite is idealized as uniformly distributed periodic arrays of aligned cells, and each cell consists of a cylindrical inclusion surrounded by a plastically deforming matrix. In the analysis, the non-uniform deformation field of the cell is decomposed into the sum of the first order approximate field and the trial additional deformation field. The precise deformation field are determined based on the minimum strain energy principle. Systematic calculation results are presented for the influence of reinforcement volume fraction and shape on the overall mechanical behavior of the composites. The results are in good agreement with the existing finite element analyses and the experimental results. This paper attempts to stimulate the work to get the analytical constitutive relation of short-fiber/particle reinforced composites.     

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)     

On the influence of material non-linearities in geometric modeling of kink band instabilities in unidirectional fiber composites

The axial compressive failure of aligned fiber composites triggered by kink band instabilities is the topic of investigation herein. Particular emphasis is put on the accurate prediction of the post-failure regime, where fiber composites are known to exhibit substantial post-failure strength. In this regard, a previous analytical model, based on geometric relationships and energy principles, is enhanced by consistently taking into account material non-linearities. Therefore, a non-linear constitutive law is introduced herein based on a newly developed exponential formulation. This non-linear constitutive law is subsequently implemented into the stress–strain response in interlaminar shearing as well as the compression response. The model enhancements are validated against published experimental data yielding excellent comparisons. Furthermore, the relevance of modeling non-linear material behavior in interlaminar dilation and bending is assessed using a bilinear constitutive law. However, implementing non-linear material behavior does not yield any improvements and can therefore be neglected.     

Nonlinear Fractional Order Viscoelasticity at Large Strains

In this paper, we formulate a fractional order viscoelastic model for large deformations and develop an algorithm for the integration of the constitutive response. The model is based on the multiplicative split of the deformation gradient into elastic and viscous parts. Further, the stress response is considered to be composed of a nonequilibrium part and an equilibrium part. The viscous part of the deformation gradient (here regarded as an internal variable) is governed by a nonlinear rate equation of fractional order. To solve the rate equation the finite element method in time is used in combination with Newton iterations. The method can handle nonuniform time meshes and uses sparse quadrature for the calculations of the fractional order integral. Moreover, the proposed model is compared to another large deformation viscoelastic model with a linear rate equation of fractional order. This is done by computing constitutive responses as well as structural dynamic responses of fictitious rubber materials.     

Nonlinear Fractional Order Viscoelasticity at Large Strains

In this paper, we formulate a fractional order viscoelastic model for large deformations and develop an algorithm for the integration of the constitutive response. The model is based on the multiplicative split of the deformation gradient into elastic and viscous parts. Further, the stress response is considered to be composed of a nonequilibrium part and an equilibrium part. The viscous part of the deformation gradient (here regarded as an internal variable) is governed by a nonlinear rate equation of fractional order. To solve the rate equation the finite element method in time is used in combination with Newton iterations. The method can handle nonuniform time meshes and uses sparse quadrature for the calculations of the fractional order integral. Moreover, the proposed model is compared to another large deformation viscoelastic model with a linear rate equation of fractional order. This is done by computing constitutive responses as well as structural dynamic responses of fictitious rubber materials.     

缝合式C/SiC复合材料非线性本构关系及断裂行为研究

C/SiC复合材料具有高比强度、高比模量和优良的热稳定性能等一系列优点, 广泛应用于航空航天领域中. 裂纹扩展进而引起的脆性断裂是其主要失效形式之一, 因而材料的断裂性能分析对材料的结构设计和应用有重要的指导意义. 本文开展了缝合式C/SiC复合材料简单力学试验和断裂试验, 研究了材料在不同载荷下的力学响应及断裂特征. 基于缝合式C/SiC复合材料简单力学试验, 建立了材料宏观非线性损伤本构方程, 并模拟了缝合式C/SiC复合材料单边切口梁和双悬臂梁的断裂行为. 本构方程采用简单函数描述了材料在复杂应力状态下的非线性应力-应变曲线, 并考虑了反向加载过程中造成的裂纹闭合. 基于商业有限元软件ABAQUS, 通过编写UMAT子程序实现非线性损伤本构方程, 采用单个单元验证了建立的本构方程的有效性. 在此基础上, 采用线弹性损伤本构和非线性损伤本构分别模拟了缝合式C/SiC复合材料单边切口梁和双悬臂梁的断裂行为. 采用非线性损伤本构方程模拟的力-位移曲线结果与试验结果更为吻合, 非线性损伤本构预测的失效载荷与试验失效载荷更为接近, 验证了所建立的非线性损伤本构方程的准确性, 为C/SiC复合材料断裂行为的研究提供了借鉴, 为缝合式C/SiC复合材料结构的设计和应用提供了理论基础.      

STUDY ON NONLINEAR CONSTITUTIVE RELATIONSHIP AND FRACTURE BEHAVIOR OF STITCHED C/SiC COMPOSITES$^{\bf 1)}$

C/SiC复合材料具有高比强度、高比模量和优良的热稳定性能等一系列优点, 广泛应用于航空航天领域中. 裂纹扩展进而引起的脆性断裂是其主要失效形式之一, 因而材料的断裂性能分析对材料的结构设计和应用有重要的指导意义. 本文开展了缝合式C/SiC复合材料简单力学试验和断裂试验, 研究了材料在不同载荷下的力学响应及断裂特征. 基于缝合式C/SiC复合材料简单力学试验, 建立了材料宏观非线性损伤本构方程, 并模拟了缝合式C/SiC复合材料单边切口梁和双悬臂梁的断裂行为. 本构方程采用简单函数描述了材料在复杂应力状态下的非线性应力-应变曲线, 并考虑了反向加载过程中造成的裂纹闭合. 基于商业有限元软件ABAQUS, 通过编写UMAT子程序实现非线性损伤本构方程, 采用单个单元验证了建立的本构方程的有效性. 在此基础上, 采用线弹性损伤本构和非线性损伤本构分别模拟了缝合式C/SiC复合材料单边切口梁和双悬臂梁的断裂行为. 采用非线性损伤本构方程模拟的力-位移曲线结果与试验结果更为吻合, 非线性损伤本构预测的失效载荷与试验失效载荷更为接近, 验证了所建立的非线性损伤本构方程的准确性, 为C/SiC复合材料断裂行为的研究提供了借鉴, 为缝合式C/SiC复合材料结构的设计和应用提供了理论基础.     

A macroscopic model for plastic flow in metal-matrix composites

A micromechanically based continuum model is developed to analyze the enhancement of plastic properties of particulate-reinforced metal-matrix composites over matrix materials. The composite is idealized as uniformly distributed periodic arrays of unit cells. Each unit cell consists of a rigid inclusion surrounded by a plastically deforming material. An energy method is adopted to obtain the overall constitutive relation for the composite on the basis of the local nonuniform deformation fields. Effects of particle volume fractions and shapes (e.g. whiskers, discs, etc.) as well as the matrix properties on the flow properties of the composite are obtained. The results are in good agreement with experimental observations and finite element analyses found in the literature. An explicit expression is also proposed, providing a means for evaluating various factors affecting the strength of composites.     

Micromechanical analysis of the finite elastic–viscoplastic response of multiphase composites

Nonlinear thermoelastic–viscoplastic constitutive equations for large deformations with isotropic and directional hardening, are incorporated into a micromechanical finite strain analysis. As a result of this analysis, which is based on the homogenization technique for periodic microstructures, a global thermoinelastic constitutive law is established that governs the overall response of multiphase materials under finite deformations. This constitutive law is expressed in terms of the instantaneous effective mechanical and thermal stress tangent tensors together with the instantaneous global inelastic stress tensor that represents the viscoplastic effects. Results for a thermoinelastic matrix reinforced by a hyperelastic compressible material are given that illustrate the response of fibrous and particulate composites to various types of loading.     

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.     

纺织结构复合材料冲击拉伸研究进展

纺织结构复合材料是以纺织结构作为增强体的一类复合材料, 其在一系列实际应用时往往要承受着高速冲击拉伸、冲击压缩等冲击载荷(冲击加载) 的作用. 因为纺织结构的整体性能, 纺织结构复合材料具有优异的抗冲击、抗分层与高损伤容限性能. 研究复合材料的冲击性能对于纺织结构复合材料的设计与应用具有指导作用. 本文详细介绍了纺织结构复合材料的发展, 纺织结构的种类及纺织结构复合材料的冲击拉伸性能实验与有限元分析研究情况, 同时也分析了纺织结构复合材料冲击拉伸破坏下的破坏机理研究进展, 并对纺织结构复合材料冲击拉伸性能研究的发展进行了展望.     

一种脆性材料动态本构损伤参数的计算反求方法

提出一种利用平板撞击的瞬态速度剖面响应来反求陶瓷脆性材料的动态本构损伤参数的方法。通过数值模拟再现平板撞击试验中的正平面波在飞片、样品和窗口材料中的传播情况,分析了“样品/窗口”界面瞬态速度剖面响应与材料参数之间的映射关系,进而建立了用响应反求本构参数的反问题,且通过一些数值试验验证该反问题是适定的。基于平板撞击速度剖面响应反求的方法,为快速地获得高应变率下一些传统方法难以测定出的本构参数提供新途径。     

Micro–macro analysis of shape memory alloy composites

The aim of the paper is to develop a micro–macro approach for the analysis of the mechanical behavior of composites obtained embedding long fibers of Shape Memory Alloys (SMA) into an elastic matrix. In order to determine the overall constitutive response of the SMA composites, two homogenization techniques are proposed: one is based on the self-consistent method while the other on the analysis of a periodic composite. The overall response of the SMA composites is strongly influenced by the pseudo-elastic and shape memory effects occurring in the SMA material. In particular, it is assumed that the phase transformations in the SMA are governed by the wire temperature and by the average stress tensor acting in the fiber. A possible prestrain of the fibers is taken into account in the model.Numerical applications are developed in order to analyze the thermo-mechanical behavior of the SMA composite. The results obtained by the proposed procedures are compared with the ones determined through a micromechanical analysis of a periodic composite performed using suitable finite elements.Then, in order to study the macromechanical response of structural elements made of SMA composites, a three-dimensional finite element is developed implementing at each Gauss point the overall constitutive laws of the SMA composite obtained by the proposed homogenization procedures. Some numerical applications are developed in order to assess the efficiency of the proposed micro–macro model.     

Nonlinear viscoelastic-degradation model for polymeric based materials

This study presents a phenomenological constitutive model for describing response of solid-like viscoelastic polymers undergoing degradation. The model is expressed in terms of recoverable and irrecoverable time-dependent parts. We use a time-integral function with a nonlinear integrand for the recoverable part and another time-integral function is used for the irrecoverable part, which is associated with the degradation evolution in the materials. Here, the degradation is attributed to the secondary and tertiary creep stages. An ‘internal clock’ concept in viscoelastic materials is used to incorporate the accelerated failure in the materials at high stress levels. We ignore the effect of heat generation due to the dissipation of energy and possible healing in predicting the long-term and failure response of the polymeric materials. Experimental data on polymer composites reported by Drozdov (2011) were used to characterize the material parameters and validate the constitutive model. The model is shown capable of predicting response of the polymer composites under various loading histories: creep, relaxation, ramp loading with a constant rate, and cyclic loadings. We observed that the failure time and number of cycles to failure during cyclic loadings are correlated to the duration of loading and magnitude of the prescribed mechanical loads. A scalar degradation variable is also introduced in order to determine the severity of the degradation in the materials, which is useful to predict the lifetime of the structures subject to various loading histories during the structural design stage.     

Micromechanical simulation on strength and ductility of two kinds of Al-based nanostructural materials

The nanostructured Al-based composites possess the combination of high yield strength and good ductility. In this paper, a micromechanical model is presented to simulate the mechanical response of bimodal nanostructured Al and the particle-reinforced aluminum matrix composite(PAMC). The constitutive relations for different phases are addressed in the model, as well as the contribution of microcracks. Numerical results show that the model can successfully describe the enhanced strength and ductility of the bimodal nanostructured Al, and the predictions of the PAMC are in good agreement with the experimental data. It is worth noting that the strength and ductility are sensitive to the volume fraction of constituents and the distribution of microcracks in both bimodal nanostructured Al and PAMC. Therefore, the present theoretical results can be used to optimize the microstructure for improving the mechanical properties of nanostructured Al-based composites.     

基于J-C本构模型的Comp.B炸药落锤冲击数值模拟

为研究Comp.B炸药在惯性冲击下的力学响应特性, 基于Johnson-Cook本构模型, 对 点火前的Comp.B炸药大落锤(400kg)冲击实验进行了数值模拟. 在不同应力加载条件下模 拟计算得到了用状态方程和不用状态方程的$\sigma$-$t$曲线, 对比应力峰值、应 力上升时间等主要力学参数, 结果表明考虑了状态方程的数值模 拟结果与实验吻合最好, 可以很好地模拟炸药惯性冲击下的力学响应.     

Investigation of stress-strain behavior of single walled carbon nanotube/rubber composites by a multi-scale finite element method

The excellent properties of carbon nanotubes have generated technological interests in the development of nanotube/rubber composites. This paper describes a finite element formulation that is appropriate for the numerical prediction of the mechanical behavior of rubber-like materials which are reinforced with single walled carbon nanotubes. The considered composite material consists of continuous aligned single walled carbon nanotubes which are uniformly distributed within the rubber material. It is assumed that the carbon nanotubes are imperfectly bonded with the matrix. Based on the micromechanical theory, the mechanical behavior of the composite may be predicted by utilizing a representative volume element. Within the representative volume element, the reinforcement is modeled according to its atomistic microstructure. Therefore, non-linear spring-based line elements are employed to simulate the discrete geometrical structure and behavior of the single-walled carbon nanotube. On the other hand, the matrix is modeled as a continuum medium by utilizing solid elements. In order to describe its behavior an appropriate constitutive material model is adopted. Finally, the interfacial region is simulated via the use of special joint elements of variable stiffness which interconnect the two materials in a discrete manner. Using the proposed multi-scale model, the stress-strain behavior for various values of reinforcement volume fraction and interfacial stiffness is extracted. The influence of the single walled carbon nanotube addition within the rubber is clearly illustrated and discussed.     

帘线/橡胶复合材料各向异性黏-超弹性本构模型

帘线/橡胶复合材料广泛应用于轮胎等重要工程领域,为了描述其在服役条件下的大变形、非线性、各向异性和高应变率等材料力学行为,基于纤维增强复合材料连续介质力学理论,提出了一种考虑应变率效应的帘线/橡胶复合材料各向异性黏-超弹性本构模型. 该模型中单位体积的应变能被解耦为便于参数识别的基体等容变形能、帘线拉伸变形能、剪切应变能和黏性应变能四部分. 给出了模型参数的确定方法,并通过拟合文献中单轴拉伸、偏轴拉伸实验数据,得到了模型参数. 利用该模型预测了不同加载和变形条件下的力学行为,并将预测结果与实验结果对比分析. 结果表明, 考虑黏性模型和不考虑黏性模型对不同应变率变形条件下的预测结果相差很大,且考虑黏性模型的预测结果与实验结果吻合很好. 因此,与不考虑黏性模型相比,所提出的各向异性黏-超弹性本构模型能更好地表征帘线/橡胶复合材料在大变形、高应变率条件下的力学特性.