Finite-volume direct averaging micromechanics of heterogeneous materials with elastic–plastic phases

In this communication, we extend the recently re-constructed micromechanics model called high-fidelity generalized method of cells (Bansal, Y., Pindera, M.-J., 2005. A second look at the higher-order theory for periodic multiphase materials. J. Appl. Mech. 72 (2), 177–195.) by incorporating inelastic response capability for the individual phases. The re-construction, based on the local/global stiffness matrix approach, has simplified the model’s theoretical framework and substantially increased its computational efficiency as well as implementability, enabling analysis of unit cells with realistic multiphase microstructures previously unattainable in the original formulation developed by Aboudi et al. (Aboudi, J., Pindera, M-J., Arnold, S.M., 2001. Linear thermoelastic higher-order theory for periodic multiphase materials. J. Appl. Mech. 68 (5), 697–707; Aboudi, J., Pindera, M-J, Arnold, S.M., 2003. Higher-order theory for periodic multiphase materials with inelastic phases. Int. J. Plasticity 19 (6), 805–847.) with an accuracy approaching finite-element solutions. Just as importantly, the re-construction has revealed the model to be based on a finite-volume, direct averaging approach with clearly discernible similarities to, and differences with, the finite-element method and the finite-volume technique used in computational fluid mechanics. Herein, easily programmable closed-form expressions have been derived for the thermo-inelastic contributions to the local stiffness matrix equations that facilitate incorporation of different inelastic constitutive theories for the phase response. The re-constructed model is then employed to investigate orientational and architectural effects in unidirectional metal matrix composites characterized by multi-inclusion unit cells. Classical incremental plasticity theory with isotropic hardening is employed for the matrix response for consistency and comparison with previously reported results by Aboudi et al. (2003). Unit cells representative of a square array of fibers rotated by an angle about the fiber axis, which lack planes of material symmetry in the rotated coordinate system in which the micromechanical analysis is performed, belong in the first category. New results are presented for such rotated unidirectional porous composites which suggest guidelines for optimizing stiffness and ductility of this class of light-weight materials relative to dominant loading directions. Strengthening effects due to fiber clustering, which require highly discretized multi-inclusion unit cells, fall in the second category. It is demonstrated that the previously documented results for particulate composites, explained by the clustering-induced alteration of stress invariants which govern plastic strain evolution, are recovered for unidirectional composites as well.     

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.     

A time dependent micro-mechanical fiber composite model for inelastic zone growth in viscoelastic matrices

In this work, a fiber composite model is developed to predict the time dependent stress transfer behavior due to fiber fractures, as driven by the viscoelastic behavior of the polymer matrix, and the initiation and propagation of inelastic zones. We validate this model using in situ, room temperature, micro-Raman spectroscopy fiber strain measurements. Multifiber composites were placed under constant load creep tests and the fiber strains were evaluated with time after one fiber break occurred. These composite specimens ranged in fiber volume fraction and strain level. Comparison between prediction and MRS measurements allows us to characterize key in situ material parameters, the critical matrix shear strain for inelastic zones and interfacial frictional slip shear stress. We find that the inelastic zone is predominately either shear yielding or interfacial slipping, and the type depends on the local fiber spacing.     

STATISTIC MODELING OF THE CREEP BEHAVIOR OF METAL MATRIX COMPOSITES BASED ON FINITE ELEMENT ANALYSIS

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｃｒｅｅｐｂｅｈａｖｉｏｒｏｆｓｈｏｒｔｆｉｂｅｒｒｅｉｎｆｏｒｃｅＭｅｔａｌＭａｔｒｉｘＣｏｍｐｏｓｉｔｅｓ (ＭＭＣｓ)ｄｅｐｅｎｄｓｏｎｔｈｅｆｏｌｌｏｗｉｎｇｆａｃｔｏｒｓ,ｓｕｃｈａｓｔｈｅｃｒｅｅｐｐｒｏｐｅｒｔｙｏｆｔｈｅｍａｔｒｉｘ ,ｅｌａｓｔｉｃａｎｄｆｒａｃｔｕｒｅｓｐｒｏｐｅｒｔｉｅｓｏｆｔｈｅｆｉｂｅｒ,ｇｅｏｍｅｔｒｉｃｐａｒａｍｅｔｅｒｓｏｆｔｈｅｆｉｂｅｒｓ,ａｒｒａｎｇｅｍｅｎｔｏｆｔｈｅｆｉｂｅｒｓａｎｄｔｈｅｐｒｏｐｅｒｔｙｏｆｔｈｅｆ…     

Elasto-plastic constitutive relations of unidirectional fiber composites for cyclic loadings

The clasto-plastic constitutive behaviors of continuous fiber reinforced composites under cyclic loadings are studied by the micromechanics method in which the equal-strain model is used in the fiber direction, the equal-stress model in the other directions. It is supposed that fiber is linearly elastic and matrix is clastic-viscoplastic. The constitutive equations of the matrix are described by Bodner-Partom's unified constitutive theory. Boron/Aluminum composite, as an example, is investigated in detail for an understanding of the stress-strain relations and initial yield behaviors of metal matrix composites. Present results are compared with the experimental data.The project was supported by the Chinese National Natural Science Foundation.     

Analysis of elastic crack bridging in ceramic matrix composites

A micromechanics analytical model based on the consistent shear lag theory is developed for predicting the failure modes in fiber reinforced unidirectional stiff matrix composites. The model accounts for a relatively large matrix stiffness and hence its load carrying capacity. The fiber and matrix stresses are established as functions of the applied stress, crack geometry, and the microstructural properties of the constituents. From the predicted stresses, the mode of failure is established based on a point stress failure criterion. The role of the microstructural parameters of the constituents on the failure modes such as self-similar continuous cracking, crack bridging and debonding parallel to the fibers is assessed.     

A multi-axial constitutive model for metal matrix composites

Metal matrix composites (MMCs) generally do not follow the classical plasticity theory, even though the matrix metals do deform plastically. A tension-compression yield asymmetry is typically observed in MMCs. For particulate-reinforced MMCs, this non-classical response is mainly due to the variation of damage evolution with loading modes. In this paper, a viscoplastic multi-axial constitutive model for plastic deformation of MMCs is constructed using the Mises-Schleicher yield criterion. The subsequent plastic flow is characterized by an associated and decomposed flow rule considering effects from both deviatoric and hydrostatic stresses. This model is capable of describing the multi-axial yield and flow behavior of MMCs by using simulated or measured asymmetric tensile and compressive stress-strain responses as input. As an example, the influence of damage evolution in terms of interfacial debonding in MMCs (obtained from FEM simulations) is incorporated through the different tensile and compressive stress-strain behaviors. Applying this model to predict the torsion and the pressure-dependant tensile responses of some commonly used MMCs provides good agreement with experimental data.     

一种预测颗粒增强复合材料界面力学性能的新方法

界面在颗粒增强复合材料中起到传递载荷的关键作用,界面性能对复合材料整体力学行为产生重要影响.然而由于复合材料内部结构较为复杂,颗粒与基体间的界面强度和界面断裂韧性难以确定,尤其是法向与切向界面强度的分别预测缺乏有效方法.本文以氧化锆颗粒增强聚二甲基硅氧烷(PDMS)复合材料为研究对象,提出一种预测颗粒增强复合材料界面力...     

纤维增强金属基复合材料整体蠕变性能的细观力学分析

高温下金属基复合材料的蠕变主要由基体蠕变和界面扩散蠕变两部分构成,以往的研究中常常只考虑其中一种蠕变机理,从而导致得到的规律具有较大的局限性．本文提出了一种可预测金属基复合材料整体蠕变性能的细观力学方法,同时考虑了基体蠕变和界面扩散蠕变两种蠕变机理,导出了具有张量形式并满足不可压缩性的界面扩散蠕变应变表达式．采用Mori-Tanaka法和自洽法二者结果的平均以便更准确地计算纤维中的应力,揭示了两种蠕变机理相互影响的竞争关系．研究了恒定双轴荷载下的总体蠕变和固定位移约束下的应力松弛这两种常见蠕变问题,探究了基体蠕变与界面扩散蠕变两种蠕变机理在总蠕变中发挥的作用,考察了不同加载条件和不同纤维体积分数对复合材料整体蠕变行为的影响．     

Stress relaxation in broken fibers in unidirectional composites: modeling and application to creep rupture analysis

This paper describes a model of stress relaxation in broken fibers in unidirectional metal matrix composites reinforced with long brittle fibers. A cylindrical cell with a broken fiber embedded in a power law creeping matrix is employed, and the broken fiber is assumed to have a bilinear distribution of axial stress. Then, on the basis of energy balance in the cell under constant overall strain, a relaxation equation of interfacial shear stress acting on stress recovery segments is derived in a simple form. The relaxation equation is approximated rationally and integrated to obtain an analytical solution, which is shown to agree fairly well with the numerical analysis of Du and McMeeking. (Du, Z.-Z., McMeeking, R.M., 1995. Creep models for metal matrix composites with long brittle fibres. J. Mech. Phys. Solids 43, 701–726.) Moreover, the relaxation equation is combined with Curtin's model (Curtin, W.A., 1991. Theory of mechanical properties of ceramic-matrix composites. J. Am. Ceram. Soc. 74, 2837–2845.), so that rupture time in long term creep is evaluated analytically and explicitly on the assumption of global load sharing. It is shown that the resulting relation represents well the dependence of creep rupture time on applied stress observed experimentally on a unidirectional metal matrix composite.     

高温环境下纤维复合材料蠕变损伤的细观机理研究

首先利用复合材料纤维断裂单胞模型，编制蠕变损伤子程序，对单胞模型进行蠕变损伤分析。分析了纤维／基体弹性模量比对蠕变变形、蠕变损伤以及应力场的影响。从计算结果发现，蠕变损伤首先在纤维断裂尖端起始，然后沿着一定的角度向基体外围延伸，直至完全损伤，而且纤维／基体模量比对高温环境下的复合材料蠕变损伤产生很大的影响；纤维与基体的模量相差越大，复合材料越容易变形，抵抗蠕变变形的能力就越小，蠕变损伤越严重。经过对不同韧性的基体材料进行研究，发现基体韧性低的复合材料蠕变损伤明显高于高韧性基体复合材料，表明低韧性基体复合材料抵抗蠕变破坏的能力较低。     

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

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

Higher-order theory for periodic multiphase materials with inelastic phases

An extension of a recently-developed linear thermoelastic theory for multiphase periodic materials is presented which admits inelastic behavior of the constituent phases. The extended theory is capable of accurately estimating both the effective inelastic response of a periodic multiphase composite and the local stress and strain fields in the individual phases. The model is presently limited to materials characterized by constituent phases that are continuous in one direction, but arbitrarily distributed within the repeating unit cell which characterizes the material's periodic microstructure. The model's analytical framework is based on the homogenization technique for periodic media, but the method of solution for the local displacement and stress fields borrows concepts previously employed by the authors in constructing the higher-order theory for functionally graded materials, in contrast with the standard finite-element solution method typically used in conjunction with the homogenization technique. The present approach produces a closed-form macroscopic constitutive equation for a periodic multiphase material valid for both uniaxial and multiaxial loading. The model's predictive accuracy in generating both the effective inelastic stress-strain response and the local stress and inelastic strain fields is demonstrated by comparison with the results of an analytical inelastic solution for the axisymmetric and axial shear response of a unidirectional composite based on the concentric cylinder model and with finite-element results for transverse loading.     

二维编织陶瓷基复合材料应力-应变行为

对二维编织陶瓷基复合材料拉伸应力-应变行为进行了试验研究和理论模拟. 将二维 编织结构简化为：正交铺层结构和纤维束波动结构. 基于基体随机开裂、纤维随机断裂分布 理论,得到正交铺层结构的应力-应变关系;基于体积平均方法,将纤维束波动部分进行分割, 引入强度分析模型,得到纤维束波动部分的应力-应变关系. 结合正交铺层部分和纤维束波动 部分的应力-应变关系,得到二维编织结构的应力-应变行为,理论与试验吻合较好.     

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.     

2D编织陶瓷基复合材料应力-应变行为的试验研究和模拟

本文对2D编织陶瓷基复合材料拉伸应力-应变行为进行了试验研究和理论模拟。将2D编织结构简化为：正交铺层结构和纤维束波动结构。基于基体随机开裂、纤维随机断裂的统计分布理论,得到正交铺层结构的应力-应变关系;基于体积平均方法,将纤维束波动部分分割为若干子单元;由于纤维束的波动使各子单元材料方向与加载方向不一致,因此考虑了各子单元的线性行为和非线性行为对材料响应的影响,同时引入强度分析模型,得到纤维束波动部分的应力-应变关系。结合正交铺层部分和纤维束波动部分的应力-应变关系,得到2D编织结构的应力-应变行为,理论与试验吻合较好。     

热环境中旋转运动功能梯度圆板的强非线性固有振动

研究热环境中旋转运动功能梯度圆板的非线性固有振动问题．针对金属－陶瓷功能梯度圆板,考虑几何非线性、材料物理属性参数随温度变化以及材料组分沿厚度方向按幂律分布的情况,应用哈密顿原理推得热环境中旋转运动功能梯度圆板的非线性振动微分方程．考虑周边夹支边界条件,利用伽辽金法得到了横向非线性固有振动方程,并确定了静载荷引起的静挠度．用改进的多尺度法求解强非线性方程,得出非线性固有频率表达式．通过算例,分析了旋转运动功能梯度圆板固有频率随转速、温度等参量的变化情况．结果表明,非线性固有频率随金属含量的增加而降低;随转速和圆板厚度的增大而升高;随功能梯度圆板表面温度的升高而降低．     

THERMOMECHANICAL RESPONSE OF SMA FIBER COMPOSITES WITH NON-UNIFORM TEMPERATURE DISTRIBUTION

A micromechanics method based on the High-Order-Theory developed by Aboudi et al. is used to predict the thermomechanical response of composites reinforced by shape memory alloy (SMA) fibers, and the non-uniform thermal distribution in composite arising from the process of heating or cooling is considered. The numerical development based on this model was coded to predict the thermomechanical response of shape memory alloy fiber/elastomer matrix composite subjected to thermal cycle loading. When the composite is heated, two heating ways, thermal gradients and heat source by passing an electric current through the SMA fibers are imposed on the composite respectively. Upon cooling, the first thermal boundary condition and the second thermal boundary condition are subjected to the composite respectively. A series of stress distributions and temperature distributions for different instants are calculated to reveal the interaction between the SMA material and matrix. It is useful to analyze and design the SMA actuator driven by heat source or the surface temperature.     

Uncertainties of unidirectional composite strength under tensile loading and variation of environmental condition

A probabilistic strength model is developed for unidirectional composites with fibers in hexagonal arrays. The model assumes that, a central core of broken fibers surrounded by unbroken fibers which are subjected to unidirectional tensile loading. The proposed approach consists in using a modified shear lag model to calculate the ineffective lengths and stress concentrations around fiber breaks. The main feature in the model lies in incorporating the variation of composite properties due to temperature and moisture, in order to predict degradation of fibers and matrix characteristics. The strength degradation is often seen as a result of changes in ineffective lengths at fiber breaks, leading to stress concentrations in intact neighboring fibers. As fiber breaks are intrinsically random, the variability of input data allows us to describe the probabilistic model by using the Monte-Carlo method. The sensitivities of the mechanical response are evaluated regarding the uncertainties in design variables such as Young’s modulus of fibers and matrix, fiber reference strength, shear yield stress, fiber volume fraction and shear parameter defining the shear stress in the inelastic region.     

用广义胞元法研究弱界面复合材料偏轴拉伸强度

本文用广义胞元法结合应力集中系数模型,从细观、宏观力学结合的角度,预测了弱界面复合材料偏轴拉伸强度值.用广义胞元法/高精度广义胞元法计算复合材料开裂前和开裂后的应力场,引入基体应力集中系数以得到基体真实应力.在计算真实应力时根据宏观试验现象考量是否对界面开裂后的复合材料进行刚度衰减,最终形成4种方案计算出复合材料的偏轴拉伸强度.通过对比芳纶纤维和亚麻纤维两种弱界面复合材料的偏轴拉伸强度试验值,找到了最可靠的预报方案并具有良好的预报精度.