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.     

Stationary variational estimates for the effective response and field fluctuations in nonlinear composites

This paper presents a variational method for estimating the effective constitutive response of composite materials with nonlinear constitutive behavior. The method is based on a stationary variational principle for the macroscopic potential in terms of the corresponding potential of a linear comparison composite (LCC) whose properties are the trial fields in the variational principle. When used in combination with estimates for the LCC that are exact to second order in the heterogeneity contrast, the resulting estimates for the nonlinear composite are also guaranteed to be exact to second-order in the contrast. In addition, the new method allows full optimization with respect to the properties of the LCC, leading to estimates that are fully stationary and exhibit no duality gaps. As a result, the effective response and field statistics of the nonlinear composite can be estimated directly from the appropriately optimized linear comparison composite. By way of illustration, the method is applied to a porous, isotropic, power-law material, and the results are found to compare favorably with earlier bounds and estimates. However, the basic ideas of the method are expected to work for broad classes of composites materials, whose effective response can be given appropriate variational representations, including more general elasto-plastic and soft hyperelastic composites and polycrystals.     

Estimation of very narrow bounds to the behavior of nonlinear incompressible elastic composites

Variational bounds for the effective behavior of nonlinear composites are improved by incorporating more-detailed morphological information. Such bounds, which are obtained from the generalized Hashin–Shtrikman variational principles, make use of a reference material with the same microstructure as the nonlinear composite. The geometrical information is contained in the effective properties of the reference material, which are explicitly present in the analytical formulae of the nonlinear bounds. In this paper, the variational approach is combined with estimates for the effective properties of the reference composite via the asymptotic homogenization method (AHM), and applied to a hexagonally periodic fiber-reinforced incompressible nonlinear elastic composite, significantly improving some recent results.     

Unified formulae of variational bounds for multiphase anisotropic elastic composites

Using the spherical and deviator decomposition of the polarization and strain tensors, we present a general algorithm for the calculation of variational bounds of dimension d for any type of anisotropic linear elastic composite as a function of the properties of the comparison body. This procedure is applied in order to obtain analytical expressions of bounds for multiphase, linear elastic composites with cubic symmetry where the geometric shapes of the inclusions are arbitrary. For the validation, it can be proved that for the isotropic particular case, the bounds coincide with those recently reported by Gibiansky and Sigmund. On the other hand, based on this general procedure some, classical bounds reported by Hashin for transversely isotropic composites, are reproduced. Numerical calculations and some comparisons with other models and experimental data are shown.     

An exact solution for micromechanical analysis of periodically layered composites

This paper presents an exact solution for effective properties and local fields of periodic layered composites, obtained based on variational asymptotic method for unit cell homogenization, a recently developed micromechanics modeling framework. The layered composites could be made of multiple general anisotropic layers. This solution reproduces some published results when specialized to layered composites made of two orthotropic or isotropic layers and is located within the correctly interpreted bounds of the rules of mixtures according to Voigt and Reuss. It is also emphasized in this paper that the rules of mixtures according to Voigt and Reuss are absolute upper and lower bounds for the effective properties of layered composites if the bounds are interpreted in terms of the stiffness matrix or the fourth-order elasticity tensor. The confusion in a few recent papers regarding effective Young's modulus of layered isotropic composites could exceed the Voigt and Reuss estimates is clarified.     

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.     

A method of plasticity for general aligned spheroidal void or fiber-reinforced composites

Based on the general concept of the secant moduli method, together with a new way of evaluating the average matrix effective stress originally proposed by Qiu and Weng (“A Theory of Plasticity for Porous Materials and Particle-Reinforced Composites”, ASME J. Appl. Mech. (1992), 59, 261.), a method for nonlinear effective properties of general aligned fiber or void composites is proposed. The method is capable of predicting composite (especially for porous materials) yielding under a hydrostatic load. Compared to the Tandon and Weng (“A Theory of Particle-Reinforced Plasticity,” ASME J. Appl. Mech. (1988), 55, 126.), model the proposed method always gives softer prediction in the uniaxial tension. The proposed method will predict the same nonlinear stress and strain relation as the Ponte Castaneda (“The Effective Mechanical Properties of Nonlinear Isotropic Composite,” J. Mech. Phys. Solids (1991), 39, 45.) variational model if the same estimates or bounds for the linear comparison composite are adopted.     

New bounds and estimates for porous media with rigid perfectly plastic matrixNouvelles bornes et estimations pour les milieux poreux à matrice rigide parfaitement plastique

We derive new rigorous bounds and self-consistent estimates for the effective yield surface of porous media with a rigid perfectly plastic matrix and a microstructure similar to Hashin's composite spheres assemblage. These results arise from a homogenisation technique that combines a pattern-based modelling for linear composite materials and a variational formulation for nonlinear media. To cite this article: N. Bilger et al., C. R. Mecanique 330 (2002) 127–132.     

A new variational estimate for the effective response of hyperelastic composites

A variational method for predicting the effective properties of hyperelastic composites in terms of available estimates for “hyperelastic comparison composites” is proposed. In some cases this estimate can produce a lower bound on the effective energy-density function. This nonlinear-comparison variational procedure is specialized to classes of fiber and statistically isotropic composites with the aid of appropriate choices of comparison composites with neo-Hookean phases. The end results are given in terms of closed-form expressions for the effective strain energy-density functions, from which the stress-strain relations can be extracted analytically. Explicit analytical estimates for the overall responses of composites whose phases behaviors are governed by the Gent model are obtained. The results for the fiber composites are compared with corresponding finite element simulations of periodic fiber composites as well as with other available estimates. A fine agreement between the predictions obtained via the various estimates is revealed even in the limit of infinite contrast between the properties of the phases.     

Homogenization estimates for multi-scale nonlinear composites

This work is concerned with the generalization of the “variational linear comparison” method of Ponte Castañeda (J. Mech. Phys. Solids 39 (1991) 45)) to multi-scale, random, heterogeneous material systems with nonlinear isotropic constituents. This method has the distinguishing feature of allowing the conversion of bounds or estimates that might be available for linear systems into corresponding bounds or estimates for the nonlinear composites of interest. Furthermore, the method is fairly simple to implement and quite general. General estimates are developed for two-scale systems and applied to various model composites with “particulate” and “granular” micro- and meso-structures, and compared with the corresponding results for their single-scale counterparts. It is found that the way that the material heterogeneity is distributed at the two separate scales can in most cases have a significant effect on the macroscopic behavior of the composite system.     

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.     

Variational bounds for anisotropic elastic multiphase composites with different shapes of inclusions

In the present work, unified formulae for the overall elastic bounds for multiphase transversely isotropic composites with different geometrical types of inclusions embedded in a matrix are calculated, including the spherical and long or short continuous cylindrical fiber cases. The influence of the different geometrical configurations of the inclusions on the composites is studied. The transversely isotropic effective bounds are obtained by applying the variational formulation for anisotropic composites developed by Willis, which relies on expressions for the static transversely isotropic Green’s function. Some numerical calculations and comparisons with the effective coefficients derived from the self-consistent approach, asymptotic homogenization method, and finite element method (FEM) are shown for different aspect ratio values, exhibiting good agreement.     

Unconventional Hamilton-type variational principles for nonlinear elastodynamics of orthogonal cable-net structures

According to the basic idea of classical yin-yang complementarity and modem dual-complementarity,in a simple and unified new way proposed by Luo,the unconven- tional Hamilton-type variational principles for geometrically nonlinear elastodynamics of orthogonal cable-net structures are established systematically,which can fully charac- terize the initial-boundary-value problem of this kind of dynamics.An important in- tegral relation is made,which can be considered as the generalized principle of virtual work for geometrically nonlinear dynamics of orthogonal cable-net structures in mechan- ics.Based on such relationship,it is possible not only to obtain the principle of virtual work for geometrically nonlinear dynamics of orthogonal cable-net structures,but also to derive systematically the complementary functionais for five-field,four-field,three-field and two-field unconventional Hamilton-type variational principles,and the functional for the unconventional Hamilton-type variational principle in phase space and the poten- tial energy functional for one-field unconventional Hamilton-type variational principle for geometrically nonlinear elastodynamics of orthogonal cable-net structures by the general- ized Legendre transformation given in this paper.Furthermore,the intrinsic relationship among various principles can be explained clearly with this approach.     

On the effect of small fluctuations in the volume fraction of constituents on the effective properties of composites

This study deals with three-scale composite materials comprised of nonlinear constituents. At the meso scale the composite can be considered as locally homogeneous with a macroscopic spatial variation of the constituents volume fraction. When these variations about a mean value are small, a Taylor expansion to second-order of the effective properties of the composite with respect to the fluctuations is given. This expansion can be used to discuss the beneficial or deleterious effects of clusters of inhomogeneities. It can also be used to derive new upper and lower bounds for the effective properties of nonlinear composites from dilute results. To cite this article: P. Suquet, C. R. Mecanique 333 (2005).     

Macroscopic strain potentials in nonlinear porous materials

By taking a hollow sphere as a representative volume element (RVE), the macroscopic strain potentials of porous materials with power-law incompressible matrix are studied in this paper. According to the principles of the minimum potential energy in nonlinear elasticity and the variational procedure, static admissible stress fields and kinematic admissible displacement fields are constructed, and hence the upper and the lower bounds of the macroscopic strain potential are obtained. The bounds given in the present paper differ so slightly that they both provide perfect approximations of the exact strain potential of the studied porous materials. It is also found that the upper bound proposed by previous authors is much higher than the present one, and the lower bounds given by Cocks is much lower. Moreover, the present calculation is also compared with the variational lower bound of Ponte Castañeda for statistically isotropic porous materials. Finally, the validity of the hollow spherical RVE for the studied nonlinear porous material is discussed by the difference between the present numerical results and the Cocks bound.     

Improved rigorous bounds on the effective elastic moduli of a composite material

A new method for deriving rigorous bounds on the effective elastic constants of a composite material is presented and used to derive a number of known as well as some new bounds. The new approach is based on a presentation of those constants as a sum of simple poles. The locations and strengths of the poles are treated as variational parameters, while different kinds of available information are translated into constraints on these parameters. Our new results include an extension of the range of validity of the Hashin-Shtrikman bounds to the case of composites made of isotropic materials but with an arbitrary microgeometry. We also use information on the effective elastic constants of one composite in order to obtain improved bounds on the effective elastic constants of another composite with the same or a similar microgeometry.     

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.     

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.     

Bounds for the non-local effective properties of random media—II

The paper deals with a random medium subjected to a static field with inhomogeneous mean values. Then, effective linear material parameters show dispersion, i.e. they depend on the ‘wave vector’ k of the mean field. Starting from a variational method previously developed by the authors, upper and lower bounds for k-dependent scalar effective parameters are derived in terms of two-point and three-point correlation functions of the stochastic material parameters. Taking into consideration the three-point correlation function gives a substantial improvement of the generalized Hashin-Shtrikman bounds obtained previously. In particular, composites with cell structure and arbitrary binary systems are considered. In order to illustrate the general results, numerical evaluations are carried out for effective permittivity of a binary cell material composed of nearly spherical grains of equal size.