含夹杂复合材料宏观性能研究

本文综述并评价了有关含夹杂复合材料的有效弹性模量研究的代表性工作，包括自洽理论，微分法，Ｅｓｈｅｌｂｙ－Ｍｏｒｉ－Ｔａｎａｋａ法，Ｈａｓｈｉｎ和Ｓｈｔｒｉｋｍａｎ的变分法等。指出上述理论由于没有充分考虑复合材料内部的微结构特征，如夹杂的形状、几何尺寸、分布和夹杂间的相互影响，在夹杂的体积份数较大，如大于０．３时已不能有效地预报复合材料的有效弹性模量，随后介绍了近来才发展起来的一种新方法─—相关函数积分法，理论与实验的结果的比较表明，该方法在夹杂体积份数较大时仍然有效。     

基于等效特征应变原理的复合材料有效模量细观力学分析

基于等效特征应变原理,提出了一种新的复合材料有效模量细观力学分析方法。首先,在等效特征应变原理基础上提出平均等效特征应变原理,它可用于解决有限体下任意形状(无论是凸或凹形)的单个夹杂或多个夹杂的弹性变形问题。其次,将平均等效特征应变原理与细观力学直接均匀法相结合,来分析确定复合材料的有效模量。最后利用复合材料纤维与基体的力学性能参数及纤维的体分比,借助MATLAB编程方法,预测其有效模量。通过将理论预测值与已有的的试验值、其它理论预测值进行对比,验证了新分析方法的合理性和分析精度。     

含正交排列夹杂和缺陷材料的等效弹性模量和损伤

研究含正交排列夹杂和缺陷材料的等效弹性模量和损伤,推导了以Eshelby－Mori－Tanaka方法求解多相各向异性复合材料等效弹性模量的简便计算公式,针对含三相正交椭球状夹杂的正交各向异性材料,得到了由细观参量（夹杂的形状、方位和体积分数）表示的等效弹性模量的解析表达式．在此基础上,提出了一个宏细观结合的正交各向异性损伤模型,从而建立了以细观量为参量的含损伤材料的应力应变关系．最后,对影响材料损伤的细观结构参数进行了分析．     

Prediction of the overall moduli of a cylindrical short-fiber reinforced composite

With respect to obtaining the effective elastic moduli of the composite, the present theory differs from both Eshelby's equivalent inclusion method and Hill's self-consistent one, both of which only consider the mechanical properties of the matrix and inclusions (fibers). In fact, the inclusion-inclusion interaction is more pronounced when the volume fraction of inclusions of the composite increases. Hence, in this paper the effective elastic moduli of the composite are derived by taking into account the shapes, sizes and distribution of inclusions, and the interactions between inclusions. In addition, it is more convincing to assume short-fibers as cylindrical inclusions as in the present paper than as ellipsoidal ones as in others^[7,8]. Finally, numerical results are given.     

考虑孔隙影响的混凝土有效弹性模量预测分析研究

本文通过纳米压痕实验技术得到混凝土材料细观各相参数，基于渐进均匀化理论，采用蒙特卡洛方法和双向游走方法建立了含孔隙混凝土的胞元模型．分析了孔隙在冻融循环次数增加情况下对混凝土有效弹性模量的影响，同时与有限元模拟分析进行了比较．结果表明：随着冻融循环次数增加，孔隙体积分数增大，界面与砂浆压痕模量相对降低，但对骨料影响较小，导致混凝土宏观弹性模量随之降低；理论分析预测的混凝土有效弹性模量与有限元模拟结果吻合良好．应用含孔隙混凝土胞元模型能有效地预测混凝土宏观弹性模量，进而也为其在冻融作用下老化演变机理的研究评估提供了基础．     

On the stiffness of materials containing a disordered array of microscopic holes or hard inclusions

We study the macroscopic mechanical behavior of materials with microscopic holes or hard inclusions. Specifically, we deal with the effective elastic moduli of composites whose microgeometry consists of either soft or hard isolated inclusions surrounded by an elastic matrix. We approach this problem by taking the stiffness of the inclusion phase to be a complex variable, which we eventually evaluate at the soft or hard limits. Our main result states that there is a certain class of non-physical, negative-definite values of the elastic moduli of the inclusion phase for which the effective tensor does not have infinities or become otherwise singular.We present applications of this result to the estimation of effective moduli and to homogenization theorems. The first application involves using complexanalytic methods to obtain rigorous and accurate bounds on the effective moduli of the high-contrast composites under consideration. We also discuss the variational estimates of Rubenfeld & Keller, which yield a complementary set of bounds on these moduli. The best bounds are given by a combination of the analytical and variational results. As a second application, we show that certain known theorems of homogenization for materials with holes are simple consequences of our main result, and in this connection we establish corresponding new theorems for materials with hard inclusions. While our rederivation of the homogenization theorems for materials with holes can be closely related to other known constructions, it appears that certain elements provided by our main result are essential in the proof of homogenization for the hard-inclusion case.     

含有随机夹杂非均匀体的有效弹性模量

在对含有随机夹杂的非均匀体求有效弹性模量时,一般多根据Eshelby的等效夹杂法,但由于该方法没有充分考虑非均匀体内部的微结构,所以其理论具有一定的局限性。本文认为Kunin的微结构理论与Eshelby的等效夹杂法相比更具一般性,因而本文采用了文[9]中一些合理的思想,摒弃了其中不合理的假设,并且建立了一种新的理论模型.最后,本文针对球夹杂的情况给出了非均匀体有效弹性模量依赖于夹杂体积份数的关系,并将该结果与文[10]中的结果进行了比较.     

碳纳米管增强Nb-Si基复合材料界面应力传递的研究

基于均匀化理论,建立了碳纳米管增强Nb-Si基复合材料的代表体积元模型,并采用剪切滞后模型对碳纳米管增强Nb-Si基复合材料界面上的应力分布和传递机制进行了研究,探讨了分子间作用力、杨氏模量比β、长径比α、体积分数?等对其应力分布和传递机制的影响。结果表明,复合材料界面应力分布的变化主要集中在碳纳米管的两端,最大的应力都是分布在加载端或拔出端,然后向另一端传递;分子间作用力、杨氏模量比、长径比、体积分数等参数对界面应力的传递均有一定的影响,其中长径比和体积分数的影响最明显,体积分数为0.02时拔出端的界面剪切应力值相对于体积分数为0.0025时增大幅度达到近7倍,而长径比从200减小到50时,其应力传递距离增大了近8倍。     

Dynamic bulk and shear moduli due to grain-scale local fluid flow in fluid-saturated cracked poroelastic rocks: Theoretical model

Grain-scale local fluid flow is an important loss mechanism for attenuating waves in cracked fluid-saturated poroelastic rocks. In this study, a dynamic elastic modulus model is developed to quantify local flow effect on wave attenuation and velocity dispersion in porous isotropic rocks. The Eshelby transform technique, inclusion-based effective medium model (the Mori–Tanaka scheme), fluid dynamics and mass conservation principle are combined to analyze pore-fluid pressure relaxation and its influences on overall elastic properties. The derivation gives fully analytic, frequency-dependent effective bulk and shear moduli of a fluid-saturated porous rock. It is shown that the derived bulk and shear moduli rigorously satisfy the Biot-Gassmann relationship of poroelasticity in the low-frequency limit, while they are consistent with isolated-pore effective medium theory in the high-frequency limit. In particular, a simplified model is proposed to quantify the squirt-flow dispersion for frequencies lower than stiff-pore relaxation frequency. The main advantage of the proposed model over previous models is its ability to predict the dispersion due to squirt flow between pores and cracks with distributed aspect ratio instead of flow in a simply conceptual double-porosity structure. Independent input parameters include pore aspect ratio distribution, fluid bulk modulus and viscosity, and bulk and shear moduli of the solid grain. Physical assumptions made in this model include (1) pores are inter-connected and (2) crack thickness is smaller than the viscous skin depth. This study is restricted to linear elastic, well-consolidated granular rocks.     

松质骨弹性模量计算的均匀化方法

对松质骨建立了六种单胞微观结构模型，采用均匀化方法和有限元方法计算松质骨的宏观等效弹性模量。给出了六种单胞模型的松质骨弹性模量与材料密度(体分比）的关系，与实验数据进行了对比，分析了不同微观结构模型在不同骨骼中的应用。结果表明，本文方法及六种单胞模型可以对松质骨微观结构和材料性能进行有效的模拟计算。同时本文又着重对松质骨的宏观等效弹性模量与体分比的指数关系进行了探讨。     

Magneto-electro-elastic coated inclusion problem and its application to magnetic-piezoelectric composite materials

In this work, a micromechanical model for the estimate of the magneto-electro-elastic behavior of the magnetic-piezoelectric composites with coated reinforcements is proposed. The coating is considered as a thin layer with properties different from those of the inclusion and the matrix. The micromechanical approach based on the Green’s functions techniques and on the interfacial operators is designed for solving the magneto-electro-elastic inhomogeneous coated inclusion problem. The effective magneto-electro-elastic properties of the composite containing thinly coated inclusions are obtained through the Mori–Tanaka’s model. Numerical investigations into magneto-electro-elastic moduli responsible for the magneto-electric coupling are presented as functions of the volume fraction and characteristics of the coated inclusions. Comparisons with existing models are presented for various shape and orientation of the coated inclusions.     

Stress intensity factors for ring-shaped crack surrounded by spherical inclusions

Interaction of a ring-shaped crack with inhomogeneities such as inclusions is analyzed for the resulting three-dimensional stress field. Considered for the composite solid with a given volume fraction of inclusions are the two cases of (a) spherical voids and (b) spherical inclusions with elastic moduli different from the matrix. A ring-shaped crack is initiated at the equator of one of the voids or inclusions. A three-phase model is used to examine the interaction between the crack and surrounding inhomogeneities. Finite element method is then applied to calculate the stress intensity factor for different configurations. The effects of volume fraction of inhomogeneities, relative size of crack to inclusions, and material constants on crack behavior are discussed.     

Asymptotic determination of effective elastic properties of composite materials with fibrous square-shaped inclusions

We propose an asymptotic approach for evaluating effective elastic properties of two-components periodic composite materials with fibrous inclusions. We start with a nontrivial expansion of the input elastic boundary value problem by ratios of elastic constants. This allows to simplify the governing equations to forms analogous to the transport problem. Then we apply an asymptotic homogenization method, coming from the original problem on a multi-connected domain to a so called cell problem, defined on a characterizing unit cell of the composite. If the inclusions' volume fraction tends to zero, the cell problem is solved by means of a boundary perturbation approach. When on the contrary the inclusions tend to touch each other we use an asymptotic expansion by non-dimensional distance between two neighbouring inclusions. Finally, the obtained “limiting” solutions are matched via two-point Padé approximants. As the results, we derive uniform analytical representations for effective elastic properties. Also local distributions of physical fields may be calculated. In some partial cases the proposed approach gives a possibility to establish a direct analogy between evaluations of effective elastic moduli and transport coefficients. As illustrative examples we consider transversally-orthotropic composite materials with fibres of square cross section and with square checkerboard structure. The obtained results are in good agreement with data of other authors.     

Differential schemes for the elastic characterisation of dispersions of randomly oriented ellipsoids

The paper deals with the elastic characterisation of dispersions of randomly oriented ellipsoids: we start from the theory of strongly diluted mixtures and successively we generalise it with a differential scheme. The micro-mechanical averaging inside the composite material is carried out by means of explicit results which allows us to obtain closed-form expressions for the macroscopic or equivalent elastic moduli of the overall composite materials. This micromechanical technique has been explicitely developed for describing embeddings of randomly oriented not spherical objects. In particular, this study has been applied to characterise media with different shapes of the inclusions (spheres, cylinders and planar inhomogeneities) and for special media involved in the mixture definition (voids or rigid particles): an accurate analysis of all these cases has been studied yielding a set of relations describing several composite materials of great technological interest. The differential effective medium scheme (developed for generally shaped ellipsoids) extends such results to higher values of the volume fraction of the inhomogeneities embedded in the mixture. For instance, the analytical study of the differential scheme for porous materials (with ellipsoidal zero stiffness voids) reveals a universal behaviour of the effective Poisson ratio for high values of the porosity. This means that Poisson ratio at high porosity assumes characteristic values depending only on the shape of the inclusions and not on the elastic response of the matrix.     

多相材料微结构多目标拓扑优化设计

在采用多尺度均匀化方法求解微结构等效特性的基础上，提出了多相材料 微结构的多目标优化设计模型. 以组分材料用量为约束，采用周长控制消除棋盘格，结合有 限元方法和对偶凸规划求解技术，对两相和三相材料微结构多项等效模量的组合进行了优化 设计. 研究比较了微结构网格粗细、材料组分以及三相材料微结构优化中的两相实体材料弹 性模量相对比例不同对优化结果的影响. 数值算例验证了优化模型和优化算法的有效性，表 明了相关因素对优化结果的影响.     

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     

Homogenization methods to approximate the effective response of random fibre-reinforced Composites

In this article a fibre-reinforced composite material is modelled via an approach employing a representative volume element with periodic boundary conditions. The effective elastic moduli of the material are thus derived. In particular, the method of asymptotic homogenization is used where a finite number of fibres are randomly distributed within the representative periodic cell. The study focuses on the efficacy of such an approach in representing a macroscopically random (hence transversely isotropic) material. Of particular importance is the sensitivity of the method to cell shape, and how this choice affects the resulting (configurationally averaged) elastic moduli. The averaging method is shown to yield results that lie within the Hashin–Shtrikman variational bounds for fibre-reinforced media and compares well with the multiple scattering and (classical) self-consistent approximations with a deviation from the latter in the larger volume fraction cases. Results also compare favourably with well-known experimental data from the literature.     

Plasticity of a two-phase composite with partially debonded inclusions

The effective elastoplastic behavior of a two-phase composite consisting of partially debonded elastic inclusions and a ductile matrix is investigated by a homogenization method. The method drew information from a recent study by the authors on the effective elastic moduli of the said composite and from an energy approach suggested by Qui and Weng, J. Appl. Mech., 59, 261 [1992] to address the homogenized plastic state of the heterogeneously deformed ductile matrix. Two types of partial debonding configuration are considered; the first is on the top and bottom of the aligned oblate inclusions and the other is on the lateral surface of the prolate ones, with special reference to spherical inclusions for both types of debonding. The transversely isotropic elastoplastic properties of the partially debonded composite are found to be highly dependent upon the debonding mode and the volume concentration and shape of inclusions. A damage mechanics based on Weibull's statistical function is also proposed to study the progressive partial debonding of the initially bonded composite under pure tension and under biaxial tension, respectively, for these two types of partial debonding. It is found that the interfacial strength, particle concentration, inclusion shape and debonding mode all play significant role in the overall response of the heterogeneous system during the progressive debonding process.     

多相材料有效性能预测的高精度方法

有效介质方法是常用的细观力学方法之一.其可用于计算多相材料的有效性能,并建立材料微细观结构和宏观性能的定量关系;有助于指导新材料设计,减少试验工作量等.然而,当夹杂含量升高时,传统有效介质方法的计算精度下降.本文以两相材料为研究对象,提出一种新的参考介质,即:为更合理考虑不同夹杂颗粒间的相互作用,假定参考介质的应变是基体相平均应变和某一修正张量的双点积.在此基础上,推导了新参考介质下两相材料的有效模量表达式,并给出该修正张量的近似计算方法;通过反复更新参考介质,采用多层次均匀化思路,将本文方法进一步用于多相材料性能的预测.为验证方法的有效性,将预测结果与已有模型结果和试验数据进行对比.结果表明本文方法较已有方法更为合理、有效.当夹杂含量升高时,本文方法较传统有效介质方法的计算精度有所提升.     

Three-dimensional analysis of transient thermal stresses in functionally graded plates

An analytical solution is presented for three-dimensional thermomechanical deformations of a simply supported functionally graded (FG) rectangular plate subjected to time-dependent thermal loads on its top and/or bottom surfaces. Material properties are taken to be analytical functions of the thickness coordinate. The uncoupled quasi-static linear thermoelasticity theory is adopted in which the change in temperature, if any, due to deformations is neglected. A temperature function that identically satisfies thermal boundary conditions at the edges and the Laplace transformation technique are used to reduce equations governing the transient heat conduction to an ordinary differential equation (ODE) in the thickness coordinate which is solved by the power series method. Next, the elasticity problem for the simply supported plate for each instantaneous temperature distribution is analyzed by using displacement functions that identically satisfy boundary conditions at the edges. The resulting coupled ODEs with variable coefficients are also solved by the power series method. The analytical solution is applicable to a plate of arbitrary thickness. Results are given for two-constituent metal-ceramic FG rectangular plates with a power-law through-the-thickness variation of the volume fraction of the constituents. The effective elastic moduli at a point are determined by either the Mori–Tanaka or the self-consistent scheme. The transient temperature, displacements, and thermal stresses at several critical locations are presented for plates subjected to either time-dependent temperature or heat flux prescribed on the top surface. Results are also given for various volume fractions of the two constituents, volume fraction profiles and the two homogenization schemes.