Bounds on the non-local effective elastic properties of composites

The paper deals with the effective linear elastic behaviour of random media subjected to inhomogeneous mean fields. The effective constitutive laws are known to be non-local. Therefore, the effective elastic moduli show dispersion, i.e¹ they depend on the “wave vector” k of the mean field. In this paper the well-known Hashin-Shtrikman bounds (1962) for the Lamé parameters of isotropic multi-phase mixtures are generalized to inhomogeneous mean fields k ≠ 0. The bounds involve two-point correlations of random elastic moduli. In the limit k → ∞ the bounds converge to the exact result. The interest is focussed on composites with cell structures and on binary mixtures. To illustrate the results, numerical evaluations are carried out for a binary cell material composed of nearly spherical grains of equal size.     

Rigorous bounds on the effective moduli of composites and inhomogeneous bodies with negative-stiffness phases

We review the theoretical bounds on the effective properties of linear elastic inhomogeneous solids (including composite materials) in the presence of constituents having non-positive-definite elastic moduli (so-called negative-stiffness phases). Using arguments of Hill and Koiter, we show that for statically stable bodies the classical displacement-based variational principles for Dirichlet and Neumann boundary problems hold but that the dual variational principle for traction boundary problems does not apply. We illustrate our findings by the example of a coated spherical inclusion whose stability conditions are obtained from the variational principles. We further show that the classical Voigt upper bound on the linear elastic moduli in multi-phase inhomogeneous bodies and composites applies and that it imposes a stability condition: overall stability requires that the effective moduli do not surpass the Voigt upper bound. This particularly implies that, while the geometric constraints among constituents in a composite can stabilize negative-stiffness phases, the stabilization is insufficient to allow for extreme overall static elastic moduli (exceeding those of the constituents). Stronger bounds on the effective elastic moduli of isotropic composites can be obtained from the Hashin–Shtrikman variational inequalities, which are also shown to hold in the presence of negative stiffness.     

Improved bounds for the elastic moduli of perfectly-random composites

New upper and lower bounds are constructed for the elastic moduli of a class of isotropic composites with perfectly-random microgeometries ([1–3]), which improve upon the bounds on the elastic shear modulus given in [1].     

Bounds for the effective properties of heterogeneous plates

This paper presents new bounds for heterogeneous plates which are similar to the well-known Hashin–Shtrikman bounds, but take into account plate boundary conditions. The Hashin–Shtrikman variational principle is used with a self-adjoint Green-operator with traction-free boundary conditions proposed by the authors. This variational formulation enables to derive lower and upper bounds for the effective in-plane and out-of-plane elastic properties of the plate. Two applications of the general theory are considered: first, in-plane invariant polarization fields are used to recover the “first-order” bounds proposed by Kolpakov [Kolpakov, A.G., 1999. Variational principles for stiffnesses of a non-homogeneous plate. J. Meth. Phys. Solids 47, 2075–2092] for general heterogeneous plates; next, “second-order bounds” for n-phase plates whose constituents are statistically homogeneous in the in-plane directions are obtained. The results related to a two-phase material made of elastic isotropic materials are shown. The “second-order” bounds for the plate elastic properties are compared with the plate properties of homogeneous plates made of materials having an elasticity tensor computed from “second-order” Hashin–Shtrikman bounds in an infinite domain.     

The boundaries of the effective elastic moduli for inhomogeneous solids

The calculation of the effective elastic moduli of inhomogeneous solids, which connect the stresses and strains averaged for the material, is accompanied by certain mathematical difficulties owing to correlation relationships of arbitrary orders. Neglect of correlation relationships leads to average elastic moduli, where averaging according to Voigt and Reuss establishes boundaries containing the effective elastic moduli [1]. Approximate values of the latter can be found by taking into account the correlation relationships of the second order in both calculation schemes [2, 3]. Another method of evaluating the true moduli consists of narrowing the boundaries of Voigt and Reuss on the basis of model representations [4-6]. The approximate effective elastic moduli for a series of polycrystals with various common-angle values are presented in [7]. An analysis of the effect of the correlation relationships between the grains of a mechanical mixture of isotropic components on the effective elastic moduli is carried out in [8], although in all the papers just mentioned the use of correlative corrections to narrow the range of elastic moduli is not investigated. Below it is shown that the calculation of the correlation corrections in the second approximation allows the range for the effective moduli to be narrowed.     

Estimates for the elastic moduli of random trigonal polycrystals and their macroscopic uncertainty

Explicit expressions of the upper and lower estimates on the macroscopic elastic moduli of random trigonal polycrystals are derived and calculated for a number of aggregates, which correct our earlier results given in Pham [Pham, D.C., 2003. Asymptotic estimates on uncertainty of the elastic moduli of completely random trigonal polycrystals. Int. J. Solids Struct. 40, 4911–4924]. The estimates are expected to predict the scatter ranges for the elastic moduli of the polycrystalline materials. The concept of effective moduli is reconsidered regarding the macroscopic uncertainty of the moduli of randomly inhomogeneous materials.     

Bounds on the overall properties of composites with ellipsoidal inclusions

A new model is put forward to bound the effective elastic moduli of composites with ellipsoidal inclusions. In the present paper, transition layer for each ellipsoidal inclusion is introduced to make the trial displacement field for the upper bound and the trial stress field for the lower bound satisfy the continuous interface conditions which are absolutely necessary for the application of variational principles. According to the principles of minimum potential energy and minimum complementary energy, the upper and lower bounds on the effective elastic moduli of composites with ellipsoidal inclusions are rigorously derived. The effects of the distribution and geometric parameters of ellipsoidal inclusions on the bounds of the effective elastic moduli are analyzed in details. The present upper and lower bounds are still finite when the bulk and shear moduli of ellipsoidal inclusions tend to infinity and zero, respectively. It should be mentioned that the present method is simple and needs not calculate the complex integrals of multi-point correlation functions. Meanwhile, the present paper provides an entirely different way to bound the effective elastic moduli of composites with ellipsoidal inclusions, which can be developed to obtain a series of bounds by taking different trial displacement and stress fields.     

A micromechanics-based nonlocal constitutive equation incorporating three-point statistics for random linear elastic composite materials

A variational formulation employing the minimum potential and complementary energy principles is used to derive a micromechanics-based nonlocal constitutive equation for random linear elastic composite materials, relating ensemble averages of stress and strain in the most general situation when mean fields vary spatially. All information contained in the energy principles is retained; we employ stress polarization trial fields utilizing one-point statistics so that the resulting nonlocal constitutive equation incorporates up through three-point statistics. The variational structure is developed first for arbitrary heterogeneous linear elastic materials, then for randomly inhomogeneous materials, then for general n-phase composite materials, and finally for two-phase composite materials, in which case explicit variational upper and lower bounds on the nonlocal effective modulus tensor operator are derived. For statistically uniform infinite-body composites, these bounds are determined even more explicitly in Fourier transform space. We evaluate these in detail in an example case: longitudinal shear of an aligned fiber or void composite. We determine the full permissible ranges of the terms involving two- and three-point statistics in these bounds, and thereby exhibit explicit results that encompass arbitrary isotropic in-plane phase distributions; we also develop a nonlocal “Milton parameter”, the variation of whose eigenvalues throughout the interval [0, 1] describes the full permissible range of the three-point term. Example plots of the new bounds show them to provide substantial improvement over the (two-point) Hashin–Shtrikman bounds on the nonlocal operator tensor, for all permissible values of the two- and three-point parameters. We next discuss further applications of the general nonlocal operator bounds: to any three-dimensional scalar transport problem e.g. conductivity, for which explicit results are given encompassing the full permissible ranges of the two- and three-point statistics terms for arbitrary three-dimensional isotropic phase distributions; and to general three-dimensional composites, where explicit results require future research. Finally, we show how the work just summarized, treating elastostatics, can be generalized to elastodynamics, first in general, then explicitly for the longitudinal shear example.     

Bounds on the elastic and transport properties of two-component composites

Weshow in detail how the McCoy bounds on the effective shear modulus of a statistically isotropic composite, can be simplified and expressed in terms of the volume fraction, f₁, and two geometric parameters. ζ₁ and η₁. We simplify Silnutzer's bounds on the effective elastic moduli of fibre-reinforced composites and find they can be expressed in terms f₁and two geometric parameters, ζ'₁ and η'₁. The parameter ζ'₁ also determines bounds on the transport and optical constants of such composites. Also, the Elsayed-McCoy bounds on the transport properties of fibre-reinforced, symmetric-cell materials are shown to depend on three geometric parameters.     

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.     

Iteration method in linear elasticity of random structure composites containing heterogeneities of noncanonical shape

We consider a linearly elastic composite medium, which consists of a homogeneous matrix containing a statistically inhomogeneous random set of heterogeneities of arbitrary shape. The general integral equations connecting the stress and strain fields in the point being considered with the stress and strain fields in the surrounding points are obtained for the random fields of heterogeneities. The method is based on a recently developed centering procedure where the notion of a perturbator is introduced and statistical averages are obtained without any auxiliary assumptions such as, e.g., effective field hypothesis implicitly exploited in the known centering methods. Effective elastic moduli and the first statistical moments of stresses in the heterogeneities are estimated for statistically homogeneous composites with the general case of both the shape and inhomogeneity of the heterogeneities moduli. The explicit new representations of the effective moduli and stress concentration factors are built by the iteration method in the framework of the quasicristallite approximation but without basic hypotheses of classical micromechanics such as both the EFH and “ellipsoidal symmetry” assumption. Numerical results are obtained for some model statistically homogeneous composites reinforced by aligned identical homogeneous heterogeneities of noncanonical shape. Some new effects are detected that are impossible in the framework of a classical background of micromechanics.     

Elastic fields of 2D and 3D misfit particles in an infinite medium

In micromechanics, accurate quantification of the elastic field (stress, strain, and displacement) caused by the presence of an inclusion in an infinite body is desired for both the particle and matrix materials. Ideally, the solution should be applicable to any particle geometry or shape and for any distribution of misfit along the interface (i.e. misfit profile). This work presents a dislocation-based numerical method, that is an extension to earlier work in this journal [Lerma, J.D., Khraishi, T., Shen, Y.L., Wirth, B.D., 2003. The elastic fields of misfit cylindrical particles: a dislocation-based numerical approach. Mech. Res. Commun. 30, 325–334], for determining the elastic fields of volume misfit particles with arbitrary misfit distribution or particle shape.     

Stability of an elastoplastic sphere under external pressure

The stability of ideally plastic, elastoplastic, and reinforced elastoviscoplastic bodies subjected to large subcritical strains was investigated in [1–4], The problems solved in these papers were related to the stability of systems in which homogeneous stress and strain fields arise in the initial state. The stability of an elastic thick-walled spherical shell subjected to external pressure leading to large subcritical strains was investigated in [5]. The stability of an axisymmetric sphere of elastoplastic material subjected to large plastic strains is examined below.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 5, pp. 155–159, September–October, 1977.     

Bounds for effective elastic moduli of disordered materials

Recently P.H. Dederichs and R. Zeller (1973) have developed a formal theory of the bounds of odd order n for the effective elastic moduli of linearly elastic, disordered materials. The bounds are established by use of statistical information given in terms of correlation functions up to order n (= 1, 3, 5,…). This theory is extended to include the bounds of even order n. It is indicated how these bounds can be made optimum under the given statistical information. The results for bounds of even and odd order are obtained in forms which resemble Neumann series, containing multiple integrals up to order (n?1). These integrals can be calculated for certain materials which are classified in terms of a gradual statistical homogeneity, isotropy and disorder. Materials which possess these properties up to the correlation functions of nth order are called overall grade n materials. The optimum bounds for overall grade 2 and grade 3 materials are given explicitly. Optimum bounds for materials which are of grade ∞ in homogeneity and isotropy (i.e. (statistically) perfectly homogeneous and isotropic) and, at the same time, disordered of grade 2 or 3 are also derived. Those for grade 2 in disorder are the Z. Hashin and S. Shtrikman's (1963) bounds. Those for grade 3 are the narrowest, explicit bounds so far derived for random elastic materials. They contain within themselves the so-called self-consistent elastic moduli.     

Bounds on the elastic moduli of statistically isotropic multicomponent materials and random cell polycrystals

Minimum energy and complementary energy principles are used to derive the upper and lower bounds on the effective elastic moduli of statistically isotropic multicomponent materials in d ($d=2$ or 3) dimensions. The trial fields, involving harmonic and biharmonic potentials, and free parameters to be optimized, lead to the bounds containing, in addition to the properties and volume proportions of the material components, the three-point correlation information about the microgeometries of the composites. The relations and restrictions among the three-point correlation parameters are explored. The upper and lower bounds are specialized to symmetric cell materials and asymmetric multi-coated spheres, which are optimal or even converge in certain cases. New bounds for random cell polycrystals are constructed with particular results for random aggregates of cubic crystals.     

Effective pressure-sensitive elastoplastic behavior of particle-reinforced composites and porous media under isotropic loading

The composite under investigation consists of an elastoplastic matrix reinforced by elastic particles or weakened by pores. The material forming the matrix is pressure-sensitive. The Drucker–Prager yield criterion and a one-parameter non-associated flow rule are employed to formulate the yield behavior of the matrix. The objective of this work is to estimate the effective elastoplastic behavior of the composite under isotropic tensile and compressive loadings. To achieve this objective, the composite sphere assemblage model of Hashin [Z. Hashin, The elastic moduli of heterogeneous materials, ASME J. Appl. Mech. 29 (1962) 143–150] is used. Exact solutions are thus derived as estimations for the effective secant and tangent bulk moduli of the composite. The effects of the loading modes and phase properties on the effective elastoplastic behavior of the composite are analytically and numerically evaluated.     

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

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

Exact Results for the Homogenization of Elastic Fiber-Reinforced Solids at Finite Strain

This work is concerned with the homogenization of solids reinforced by aligned parallel continuous fibers or weakened by aligned parallel cylindrical pores and undergoing large deformations. By alternatively exploiting the nominal and material formulations of the corresponding homogenization problem and by applying the implicit function theorem, it is shown that locally homogeneous deformations can be produced in such inhomogeneous materials and form a differentiable manifold. For every macroscopic strain associated to a locally homogeneous deformation field, the effective nominal or material stress–strain relation is exactly determined and connections are also exactly established between the effective nominal and material elastic tangent moduli. These results are microstructure-independent in the sense that they hold irrespectively of the transverse geometry and distribution of the fibers or pores. A porous medium consisting of a compressible Mooney–Rivlin material with cylindrical pores is studied in detail to illustrate the general results. This work was the first time presented at the Euromech Colloqium 464 on “Fiber-reinforced Solids: Constitutive Laws and Instabilities”, September 28–October 1, 2004, Cantabria, Spain.     

Revised Bounds on the Elastic Moduli of Two-Dimensional Random Polycrystals

Our earlier derived bounds on the elastic moduli of two-dimensional random polycrystals [1, 2] involve a geometric restriction through an assumption on the form of an isotropic eight-rank tensor. The general form of the tensor is used in this study to reconstruct the bounds, which are expected to approach the scatter range for the moduli of the irregular aggregate.     

Dynamic stability analysis of an elastic composite material having a negative-stiffness phase

The rigorous classical bounds of elastic composite materials theory provide limits on the achievable composite stiffnesses in terms of the properties and arrangements of the composite's constituents. These bounds result from the assumption, presumably made for stability reasons, that each constituent material must have positive-definite elastic moduli. If this assumption is relaxed, recently published elasticity analyses and experimental measurements show these bounds can be greatly exceeded, resulting in new materials of enormous potential.The key question is whether a composite material having a non-positive-definite constituent can be stable overall in the practically useful situation of applied traction boundary conditions. Drugan [2007. Elastic composite materials having a negative-stiffness phase can be stable. Phys. Rev. Lett. 98 (5), article no. 055502] first proved the answer is yes, by applying the energy criterion of elastic stability to the basic two- and three-dimensional composites consisting of a cylinder or sphere having non-positive-definite (but strongly elliptic) moduli with a thin positive-definite coating and proving overall stability provided the coating is sufficiently stiff.Here, we perform a complete and direct dynamic stability analysis of the plane strain fundamental elastic composite consisting of a circular cylinder of non-positive-definite material firmly bonded to a positive-definite concentric coating, for the full range of coating thicknesses (i.e., volume fractions). We determine quantitatively the full permissible range of inclusion and coating moduli, as a function of coating thickness, for which the overall composite is stable under dead traction boundary conditions. Among the results, we show that in the thin-coating case, the present dynamic stability analysis leads to precisely the same analytical stability requirements as those derived via the energy criterion by Drugan [2007. Elastic composite materials having a negative-stiffness phase can be stable. Phys. Rev. Lett. 98 (5), article no. 055502], and we derive new analytical stability requirements that are valid for a wider range of coating thickness. At the other extreme, we show that in the case of very thick coatings (corresponding to the dilute case of a matrix-inclusion composite), even an inclusion with merely strongly elliptic moduli can be stabilized by a positive-definite matrix satisfying weak requirements, for which we derive analytical expressions. Overall, our results show that surprisingly weak restrictions on the moduli and thickness of the positive-definite coating are sufficient to stabilize a non-positive-definite inclusion, even one whose moduli are merely strongly elliptic. These results legitimize expanding the search for novel materials with extreme properties to those incorporating a non-positive-definite constituent, and they provide quantitative restrictions on the constituent materials’ moduli and volume fractions, for the geometry examined here, that ensure overall stability of such composite materials.