广义自洽Mori-Tanaka模型及涂层夹杂体复合材料的有效模量

经典广义自洽模型的最大不足是需要确定相材料的位移及应变场,这一过程十分繁杂,而且最后得到的有效剪切模量无法显式表达,难以应用.该文摈弃这一经典做法,而从广义自洽模型的应变等价条件出发,在夹杂应变均匀的近似假定下,将Hill界面条件应用于整个二相体内,从而得到一种可以预报涂层夹杂体复合材料有效模量的广义自洽Mori-Tanaka模型.与已有的实验及理论结果的比较表明,该模型准确可靠,而且有效体积和剪切模量均能显式表达.同时证实,以上的Hill界面条件应用于二相体内相当于Mori-Tanaka平均场近似.     

An explicit expression of the effective moduli for composite materials filled with coated inclusions

The obvious shortcoming of the generalized self-consistent method (GSCM) is that the effective shear modulus of composite materials estimated by the method can not be expressed in an explicit form. This is inconvenient in engineering applications. In order to overcome that shortcoming of GSCM, a reformation of GSCM is made and a new micromechanical scheme is suggested in this paper. By means of this new scheme, both the effective bulk and shear moduli of an inclusion-matrix composite material can be obtained and be expressed in simple explicit forms. A comparison with the existing models and the rigorous Hashin-Shtrikman bounds demonstrates that the present scheme is accurate. By a two-step homogenization technique from the present new scheme, the effective moduli of the composite materials with coated spherical inclusions are obtained and can also be expressed in an explicit form. The comparison with the existing theoretical and experimental results shows that the present solutions are satisfactory. Moreover, a quantitative comparison of GSCM and the Mori-Tanaka method (MTM) is made based on a unified scheme. The project supported by the National Natural Science Foundation of China under the Contract NO. 19632030 and 19572008, and China Postdoctoral Science Foundation     

A generalized self-consistent polycrystal model for the yield strength of nanocrystalline materials

Inspired by recent molecular dynamic simulations of nanocrystalline solids, a generalized self-consistent polycrystal model is proposed to study the transition of yield strength of polycrystalline metals as the grain size decreases from the traditional coarse grain to the nanometer scale. These atomic simulations revealed that a significant portion of atoms resides in the grain boundaries and the plastic flow of the grain-boundary region is responsible for the unique characteristics displayed by such materials. The proposed model takes each oriented grain and its immediate grain boundary to form a pair, which in turn is embedded in the infinite effective medium with a property representing the orientational average of all these pairs. We make use of the linear comparison composite to determine the nonlinear behavior of the nanocrystalline polycrystal through the concept of secant moduli. To this end an auxiliary problem of Christensen and Lo (J. Mech. Phys. Solids 27 (1979) 315) superimposed on the eigenstrain field of Luo and Weng (Mech. Mater. 6 (1987) 347) is first considered, and then the nonlinear elastoplastic polycrystal problem is addressed. The plastic flow of each grain is calculated from its crystallographic slips, but the plastic behavior of the grain-boundary phase is modeled as that of an amorphous material. The calculated yield stress for Cu is found to follow the classic Hall-Petch relation initially, but as the gain size decreases it begins to depart from it. The yield strength eventually attains a maximum at a critical grain size and then the Hall-Petch slope turns negative in the nano-range. It is also found that, when the Hall-Petch relation is observed, the plastic behavior of the polycrystal is governed by crystallographic slips in the grains, but when the slope is negative it is governed by the grain boundaries. During the transition both grains and grain boundaries contribute competitively.     

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.     

Effective elastic moduli of inhomogeneous solids by embedded cell model

An embedded cell model is presented to obtain the effective elastic moduli for three-dimensional two-phase composites which is an exact analytic formula without any simplified approximation and can be expressed in an explicit form. For the different cells such as spherical inclusions and cracks surrounded by sphere and oblate ellipsoidal matrix, the effective elastic moduli are evaluated and the results are compared with those from various micromechanics models. These results show that the present model is direct, simple and efficient to deal with three-dimensional two-phase composites. The project supported by the National Natural Science Foundation of China (No. 19704100) and the National Natural Science Foundation of Chinese Academy of Sciences (No. KJ951-1-201)     

The effective properties of piezocomposites, part II: The effective electroelastic moduli

Since piezoelectric ceramic/polymer composites have been widely used as smart materials and smart structures, it is more and more important to obtain the closed-from solutions of the effective properties of piezocomposites with piezoelectric ellipsoidal inclusions. Based on the closed-from solutions of the electroelastic Eshelby's tensors obtained in the part I of this paper and the generalized Budiansky's energy-equivalence framework, the closed-form general relations of effective electroelastic moduli of the piezocomposites with piezoelectric ellipsoidal inclusions are given. The relations can be applicable for several micromechanics models, such as the dilute solution and the Mori-Tanaka's method. The difference among the various models is shown to be the way in which the average strain and the average electric field of the inclusion phase are evaluated. Comparison between predicted and experimental results shows that the theoretical values in this paper agree quite well with the experimental results. These expression can be readily utilized in analysis and design of piezocomposites. The project supported by the National Natural Science Foundation of China     

Micromechanics predictions of the effective moduli of magnetoelectroelastic composite materials

In this paper, the self-consistent, generalized Mori–Tanaka and dilute micromechanics theories are extended to study the coupled magnetoelectroelastic composite materials. The heterogeneous inclusion problem of magnetoelectroelastic behavior is formulated in terms of five interaction tensors related to the Green's functions for an infinite three-dimensional transversely isotropic magnetoelectroelastic solid. These tensors are then used to predict the effective moduli of the magnetoelectroelastic solid based on the self-consistent, Mori–Tanaka and the dilute approaches. Numerical results are obtained for various types of inclusions. These results are employed to study the effects of the inclusion properties, such as moduli, volume fractions, shapes, etc., on the effective moduli of magnetoelectroelastic composites, in particular, the related magnetic properties. The results obtained using the self-consistent model, the generalized Mori–Tanaka's model and the dilute approach are compared with the existing experimental and theoretical results.     

双周期圆截面纤维复合材料平面问题的解析法

结合双准周期Riemann边值问题理论与Eshelby等效夹杂原理，为双周期圆截面纤维复合材 料平面问题发展了一个实用有效的解析方法，获得了问题的全场级数解并与有限元结果进行 了比较. 该方法为非均匀材料的力学性质分析和复合材料等新材料的微结构设计提供了 一个有效的计算工具，也可用来评估有限元等数值与近似方法的精度.     

Bounds for interpolating complex effective moduli of viscoelastic materials from measured data

The complex effective moduli of viscoelastic materials can be experimentally measured by a number of techniques giving relaxation data for materials at discrete frequencies. In this paper we present a method which allows one to find bounds on the moduli of the material at the unmeasured frequencies. Our bounds can be very tight, and are, therefore, good approximations to the data. Using our bounds, experimentalists can concentrate their efforts on gathering highly accurate data at only a few frequencies and numerically generate data at other frequencies that may be of interest. Received: 10 August 1999 Accepted: 21 December 2000     

The effective properties of piezoelectric composite materials with transversely isotropic spherical inclusions

ＩｎｔｒｏｄｕｃｔｉｏｎＷｉｔｈｔｈｅｄｅｖｅｌｏｐｍｅｎｔｏｆｉｎｆｏｒｍａｔｉｏｎｉｎｄｕｓｔｒｙａｎｄｔｈｅａｐｅａｒａｎｃｅｏｆｓｍａｒｔｍａｔｅｒｉａｌｓａｎｄｓｍａｒｔｓｔｒｕｃｔｕｒｅｓ,ｉｔｂｅｃｏｍｅｓｍｏｒｅａｎｄｍｏｒｅｉｍｐｏ...     

Prediction of the effective elastic moduli of materials with irregularly-shaped pores based on the pore projected areas

Statistical modeling is used to correlate geometric parameters of pores with their contributions to the overall Young’s moduli of linearly elastic solids. The statistical model is based on individual pore contribution parameters evaluated by finite element simulations for a small pore subset selected using the design of experiments approach, so there is no need to solve the elasticity problem for all pores in the material. A polynomial relating pore geometric parameters to the contribution parameters is then fitted to the results of the simulations. We found a good correlation between normalized projected areas of the pores on three coordinate planes and their contributions to the corresponding effective Young’s moduli. The model is applied and validated for two large sets of pore geometries obtained by X-ray microcomputed tomography of a carbon/carbon and a 3D woven carbon/epoxy composite specimens.     

A unified energy approach to a class of micromechanics models for composite materials

Several micromechanics models for the determination of composite moduli are investigated in this paper, including the dilute solution, self-consistent method, generalized self-consistent method, and Mori-Tanaka's method. These micromechanical models have been developed by following quite different approaches and physical interpretations. It is shown that all the micromechanics models share a common ground, the generalized Budiansky's energy-equivalence framework. The difference among the various models is shown to be the way in which the average strain of the inclusion phase is evaluated. As a bonus of this theoretical development, the asymmetry suffered in Mori-Tanaka's method can be circumvented and the applicability of the generalized self-consistent method can be extended to materials containing microcracks, multiphase inclusions, non-spherical inclusions, or non-cylindrical inclusions. The relevance to the differential method, double-inclusion model, and Hashin-Shtrikman bounds is also discussed. The application of these micromechanics models to particulate-reinforced composites and microcracked solids is reviewed and some new results are presented.     

反平面剪切作用下双材料滑动界面的细观力学模型

非理想粘结界面对多相材料力学性能具有重要影响。对于双材料间含众多随机分布微裂纹的界面，宏观上可以等效为连续损伤的弱界面，其两侧的面力连续而位移有间断。只有切线方向的位移间断，而法线方向位移连续的弱界面称之为滑动界面。在反平面剪切的作用下，我们证明了对于含有随机分布微裂纹的弹性双材料界面在宏观上等效为线弹簧型滑动界面，并获得了滑动界面柔度的一般表达式。利用Mori—Tanaka方法和广义自洽方法，我们研究了滑动界面柔度系数和微裂纹密度的关系。对这两种方法所得的结果进行比较发现，Mori—Tanaka方法得到的界面柔度比广义自洽方法得到的界面柔度大。当裂纹密度比较小时，这两种方法求得的界面柔度很接近。两种方法的结果都表明，界面柔度随裂纹密度的增加而增加。Mori—Tanaka方法比广义自治方法求解更为简便。     

快速多极虚边界元法对含圆孔薄板有效弹性模量的模拟分析

针对虚边界元法,引入快速多极展开和广义极小残值法(GMRES)的思想,以形成快速多极虚边界元法的求解思想,并将此方法用于含圆孔薄板有效弹性模量的模拟分析.由于本文方法采用了"源点"多极展开和"场点"局部展开的组合处理方案,从而使得原问题方程组求解的计算耗时量和储存量降至与所求问题的计算自由度数成线性比例.本文工作的研究目的在于:提高虚边界元法在普通台式机上的运算能力和拓宽虚边界元法对大规模复杂问题的求解(或数值模拟).文中给出了均布圆孔的正方形薄板和之字形分布圆孔薄板二个算例,以验证该方法的可行性,计算精度和计算效率.     

适用于自洽强度方法的冲击加载-再加载实验技术

针对自洽强度方法存在的冲击加载-再加载的难题,提出了一种采用较高硬度材料为支撑制作组合飞片的简便方法。利用该方法获得了铝、锡和锆基金属玻璃较理想的冲击加载-再加载粒子速度剖面,验证了该方法的有效性。由本文获得的冲击加载-再加载粒子速度剖面,并根据自洽方法,计算得到了铝、锡和锆基金属玻璃再加载过程剪应力变化数据。进一步分析表明,在本文涉及的压力范围内,仅由冲击加载-卸载实验得到的铝、锡和锆基金属玻璃屈服强度将比实际结果降低20%~50%。因此,在采用自洽方法计算高压强度时,冲击加载-再加载数据不可或缺。     

A novel implementation algorithm of asymptotic homogenization for predicting the effective coefficient of thermal expansion of periodic composite materials

Asymptotic homogenization (AH) is a general method for predicting the effective coefficient of thermal expansion (CTE) of periodic composites. It has a rigorous mathematical foundation and can give an accurate solution if the macrostructure is large enough to comprise an infinite number of unit cells. In this paper, a novel implementation algorithm of asymptotic homogenization (NIAH) is devel-oped to calculate the effective CTE of periodic composite materials. Compared with the previous implementation of AH, there are two obvious advantages. One is its implemen-tation as simple as representative volume element (RVE). The new algorithm can be executed easily using commercial finite element analysis (FEA) software as a black box. The detailed process of the new implementation of AH has been provided. The other is that NIAH can simultaneously use more than one element type to discretize a unit cell, which can save much computational cost in predicting the CTE of a complex structure. Several examples are carried out to demonstrate the effectiveness of the new implementation. This work is expected to greatly promote the widespread use of AH in predicting the CTE of periodic composite materials.     

Effective elastic moduli of periodic granular composite with transversely isotropic phases

We find a rigorous solution describing the macroscopically uniform stress state of a periodic granular composite with transversely isotropic phases. The structure of the composite is modeled by a cube containing a finite number of arbitrarily arranged and oriented, transversely isotropic spherical inclusions. This provides the model with a flexible means of describing the microstructure. Applying periodic vector solutions and local expansion formulas reduces the initial boundary-value problem to a system of linear algebraic equations. By averaging the solution over the unit cell, we derived exact finite expressions for the components of the effective stiffness tensor. The numerical data presented help to evaluate the efficiency of the method and the limits of applicability of available approximate theories.Translated from Prikladnaya Mekhanika, Vol. 40, No. 9, pp. 123–130, September 2004.     

A numerical simulation for effective elastic moduli of plates with various distributions and sizes of cracks

A fast convergent numerical model is developed to calculate the effective moduli of plates with various distributions and sizes of cracks, in which the crack line is divided into M parts to obtain the unknown traction on the crack line. When M=1, the model reduces to Kachanov's approximation method [Int. J. Solids Struct. 23 (1987) 23]. Six types of crack distributions and three kinds of crack sizes are considered, which are four regular (equilateral triangle, equilateral hexagon, rectangle, and diamond) and two random distributions (random location and orientation, and parallel orientation and random location), and one, two and random crack sizes. Some typical examples are also analyzed using the finite element method (FEM) to validate the present model. Then, the effective moduli associated with the crack distributions and sizes are calculated in detail. The present results for the regular distributions show some very interesting phenomena that have not been revealed before. And for the two random distributions, as the effective moduli depend on samples due to the randomness, the effect of the sample size and number are analyzed first. Then, effective moduli for plates with the three sizes of cracks are calculated. It is found that the effect of crack sizes on the effective moduli is significant for high crack densities, and small for low crack densities, and the random crack size leads to the lowest effective moduli. The present numerical results are compared with several popular micromechanics models to determine which one can provide the optimum estimation of the effective moduli of cracked plates with general crack densities. Furthermore, some existing numerical results are analyzed and discussed.     

A micromechanics based damage model for composite materials

The predictive capacity of ductile fracture models when applied to composite and multiphase materials is related to the accuracy of the estimated stress/strain level in the second phases or reinforcements, which defines the condition for damage nucleation. Second phase particles contribute to the overall hardening of the composite before void nucleation, as well as to its softening after their fracture or decohesion. If the volume fraction of reinforcement is larger than a couple of percents, this softening can significantly affect the resistance to plastic localization and cannot be neglected. In order to explicitly account for the effect of second phase particles on the ductile fracture process, this study integrates a damage model based on the Gologanu–Leblond–Devaux constitutive behavior with a mean-field homogenization scheme. Even though the model is more general, the present study focuses on elastic particles dispersed in an elasto-plastic matrix. After assessing the mean-field homogenization scheme through comparison with two-dimensional axisymmetric finite element calculations, an extensive parametric study is performed using the integrated homogenization-damage model. The predictions of the integrated homogenization-damage model are also compared with experimental results on cast aluminum alloys, in terms of both the fracture strain and overall stress–strain curves. The study demonstrates the complex couplings among the load transfer to second phase particles, their resistance to fracture, the void nucleation mode, and the overall ductility.     

A pattern-based method for bounding the effective response of a nonlinear composite

This paper deals with the prediction of the effective properties of nonlinear composites. Rather than bounding the effective energy, this work aims at bounding directly the effective stress-strain response, by extending a method originally introduced by Milton and Serkov (J. Mech. Phys. Solids 48 (2000) 1295) and recently refined by Talbot and Willis (Proc. Roy. Soc. 460 (2004) 2705). In this paper, bounding the effective response is achieved by introducing a linear comparison composite with the same micro-geometry as the given nonlinear composite, as Ponte Castañeda (J. Mech. Phys. Solids 39 (1991) 45) did for the energy. It is found that any lower bound for the energy of the linear comparison composite generates a corresponding bound for the stress-strain response of the nonlinear composite. A selection of examples is presented to illustrate the method and compare the bounds obtained with existing results.