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.     

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.     

New variational principles in plasticity and their application to composite materials

I , a variational method for bounding the effective properties of nonlinear composites with isotropic phases, proposed recently by ponte castañeda (J. Mech. Phys. Solids 39, 45, 1991), is given full variational principle status. Two dual versions of the new variational principle are presented and their equivalence to each other, and to the classical variational principles, is demonstrated. The variational principles are used to determine bounds and estimates for the effective energy functions of nonlinear composites with prescribed volume fractions in the context of the deformation theory of plasticity. The classical bounds of Voigt and Reuss for completely anisotropic composites are recovered from the new variational principles and are given alternative, simpler forms. Also, use of a novel identity allows the determination of simpler forms for nonlinear Hashin-Shtrikman bounds, and estimates, for isotropic, particle-reinforced composites, as well as for transversely isotropic, fiber-reinforced composites. Additionally, third-order bounds of the Beran type are determined for the first time for nonlinear composites. The question of the optimality of these bounds is discussed briefly.     

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.     

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.     

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.     

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.     

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.     

Second-order homogenization estimates for nonlinear composites incorporating field fluctuations: I—theory

This paper is concerned with the development of an improved second-order homogenization method incorporating field fluctuations for nonlinear composite materials. The idea is to combine the desirable features of two different, earlier methods making use of “linear comparison composites”, the properties of which are chosen optimally from suitably designed variational principles. The first method (Ponte Castañeda, J. Mech. Phys. Solids 39 (1991) 45) makes use of the “secant” moduli of the phases, evaluated at the second moments of the strain field over the phases, and delivers bounds, but these bounds are only exact to first-order in the heterogeneity contrast. The second method (Ponte Castañeda, J. Mech. Phys. Solids 44 (1996) 827) makes use of the “tangent” moduli, evaluated at the phase averages (or first moments) of the strain field, and yields estimates that are exact to second-order in the contrast, but that can violate the bounds in some special cases. These special cases turn out to correspond to situations, such as percolation phenomena, where field fluctuations, which are captured less accurately by the second-order method than by the bounds, become important. The new method delivers estimates that are exact to second-order in the contrast, making use of generalized secant moduli incorporating both first- and second-moment information, in such a way that the bounds are never violated. Some simple applications of the new theory are given in Part II of this work.     

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).