Multiphase laminates of extremal effective conductivity in two dimensions

This paper deals with two-dimensional composites made of several isotropic linearly conducting phases in prescribed volume fractions. The primary focus is on the three-phase case; the generalization to a larger number of phases is straightforward.A class of high- but finite-rank laminates is introduced. The laminates saturate the known inequality bounds—due to the work of Hashin and Shtrikman, Lurie and Cherkaev, Tartar, and Murat and Tartar—on the effective conductivity tensor of any composite. These bounds depend only on the constituent material properties and volume fractions and not on the placement of these materials in the composite. The bounds are known not to be optimal for all admissible choices of the conductivities and volume fractions. However, they are now known to be realizable in a much larger range of these parameters than was previously known.The range of effective properties of our multiphase laminates strictly includes those corresponding to the composites found earlier by Milton and Kohn, Lurie and Cherkaev, and Gibiansky and Sigmund. The new optimal laminates are found in a systematic fashion by satisfying sufficient conditions on the fields in each layer. This leads to a simple algorithm for generating optimal laminates.In addition a new supplementary bound for multiphase structures is also proven which must be satisfied by composites with smooth interfaces.     

Three-phase plane composites of minimal elastic stress energy: High-porosity structures

The paper establishes exact lower bound on the effective elastic energy of two-dimensional, three-material composite subjected to the homogeneous, anisotropic stress. It is assumed that the materials are mixed with given volume fractions and that one of the phases is degenerated to void, i.e., the effective composite is porous. Explicit formula for the energy bound is obtained using the translation method enhanced with additional inequality expressing certain property of stresses. Sufficient optimality conditions of the energy bound are used to set the requirements which have to be met by the stress fields in each phase of optimal effective material regardless of the complexity of its microstructural geometry. We show that these requirements are fulfilled in a special class of microgeometries, so-called laminates of a rank. Their optimality is elaborated in detail for structures with significant amount of void, also referred to as high-porosity structures. It is shown that geometrical parameters of optimal multi-rank, high-porosity laminates are different in various ranges of volume fractions and anisotropy level of external stress. Non-laminate, three-phase microstructures introduced by other authors and their optimality in high-porosity regions is also discussed by means of the sufficient conditions technique. Conjectures regarding low-porosity regions are presented, but full treatment of this issue is postponed to a separate publication. The corresponding “G-closure problem” of a three-phase isotropic composite is also addressed and exact bounds on effective isotropic properties are explicitly determined in these regions where the stress energy bound is optimal.     

Quasiconformal mappings as a tool to study certain two-dimensional G-closure problems

A recent theorem due to Astala establishes the best exponent for the area distortion of planar K-quasiconformal mappings. We use a refinement of Astala's theorem due to Eremenko and Hamilton to prove new bounds on the effective conductivity of two-dimensional composites. The bounds are valid for composites made of an arbitrary finite number n of possibly anisotropic phases in prescribed volume fractions. For n= 2 we prove the optimality of the bounds under certain additional assumptions on the G-closure parameters.     

Optimal lower bounds on the hydrostatic stress amplification inside random two-phase elastic composites

Composites made from two linear isotropic elastic materials are subjected to a uniform hydrostatic stress. It is assumed that only the volume fraction of each elastic material is known. Lower bounds on all rth moments of the hydrostatic stress field inside each phase are obtained for r?2. A lower bound on the maximum value of the hydrostatic stress field is also obtained. These bounds are given by explicit formulas depending on the volume fractions of the constituent materials and their elastic moduli. All of these bounds are shown to be the best possible as they are attained by the hydrostatic stress field inside the Hashin-Shtrikman coated sphere assemblage. The bounds provide a new opportunity for the assessment of load transfer between macroscopic and microscopic scales for statistically defined microstructures.     

Explicit cross-property correlations for anisotropic two-phase composite materials

Explicit correlations between two groups of anisotropic effective properties—conductivity and elasticity—are established for two-phase composite materials with anisotropic microstructures (non-randomly oriented inclusions of non-spherical shapes). The correlations are derived in the framework of the non-interaction approximation. The elasticity tensor is expressed in terms of the conductivity tensor in closed form. Applications to realistic microstructures, containing mixtures of diverse inclusion shapes are given. Compliance/stiffness contribution tensors of an inclusion, that characterize the inclusion's contribution to the overall elastic response, are derived in the course of analysis; these results are of interest on their own.     

Universal bounds on the electrical and elastic response of two-phase bodies and their application to bounding the volume fraction from boundary measurements

Universal bounds on the electrical and elastic response of two-phase (and multiphase) ellipsoidal or parallelopipedic bodies have been obtained by Nemat-Nasser and Hori. Here we show how their bounds can be improved and extended to bodies of arbitrary shape. Although our analysis is for two-phase bodies with isotropic phases it can easily be extended to multiphase bodies with anisotropic constituents. Our two-phase bounds can be used in an inverse fashion to bound the volume fractions occupied by the phases, and when the volume fraction is asymptotically small reduce to those of Capdeboscq and Vogelius, for electrical conductivity, and Capdeboscq and Kang, for elasticity. Other volume fraction bounds derived here utilize information obtained from thermal, magnetic, dielectric or elastic responses. One bound on the volume fraction can be obtained by simply immersing the body in a water filled cylinder with a piston at one end and measuring the change in water pressure when the piston is displaced by a known small amount. This bound may be particularly effective for estimating the volume of cavities in a body. We also obtain new bounds utilizing just one pair of (voltage, flux) electrical measurements at the boundary of the body.     

Strong-contrast expansion correlation approximations for the effective elastic moduli of multiphase composites

Following the previous approach of Pham and Torquato (J Appl Phys 94:6591–6602, 2003) and Torquato (J Mech Phys Solids 45:1421–1448, 1997; Random heterogeneous media, Springer, Berlin, 2002), we derive the strong-contrast expansions for the effective elastic moduli K _e,G _e of d-dimensional multiphase composites. The series consists of a principal reference part and a fluctuation part (perturbation about a homogeneous reference or comparison material), which contains multi-point correlation functions that characterize the microstructure of the composite. We propose a three-point correlation approximation for the fluctuation part with an objective choice of the reference phase moduli, such that the fluctuation terms vanish. That results in the approximations for the effective elastic moduli of isotropic composites, which coincide with the well-known self-consistent and Maxwell approximations for two-phase composites having respective microstructures. Applications to some two-phase materials are given.     

Basic equations for the microscopic theory of plasticity of polycrystalline materials with given texture

An expression for the yield stress of anisotropic materials is applied to the anisotropic strength of hard rolled copper foils whose crystallographic texture is known. We assume that this crystallographic texture is the only cause of the anisotropic plastic behaviour of the material. The model used for the yield stress is also used to deduce:

1. Stress-strain relations for isotropic polycrystalline materials;
2. A formula for the fully plastic strain tensor, applied to anisotropic hard rolled copper foils.

For the anisotropic copper foils considered the calculated curves of the yield stress and of the strain tensor as a function of the angle x between rolling and tensile direction agree qualitatively with the measured values. However, the theory is not complete, since the yield stress and the plastic strain tensor are both a function of a parameter Q, the fraction of the number of available crystallographic slip planes on which the maximum shear stress has reached the critical value τ_a. We assume that for “fully” plastic deformation a certain critical fraction Q _e of the total number of slip planes has to be active. The fraction Q _e is called the critical active quantity. With the parameter Q _e we adjust the calculated curves to the measured ones. The dependence of Q _e on the properties of the material (e.g. the crystallographic texture) is discussed in Appendix I.     

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.     

An upper bound to the free energy of mixing by twin-compatible lamination for n-variant martensitic phase transformations

Modeling the energetic behavior of martensitic (phase transforming) materials usually leads to non quasiconvex energy formulations. For this reason, researchers often employ quasiconvex relaxation methods to improve the character of the formulation. Unfortunately, explicit expressions for the relaxed free energy density for multi-variant martensitic materials are typically not available. Thus, some researchers have employed a Reuβ-like convex lower bound, which neglects compatibility constraints, as an estimate on the free energy of mixing. To be confident with such a technique, one needs a measure of the quality of the lower bound. In this paper, we seek such a measure by comparing the Reuβ-like lower bound to an upper bound. The upper bound is constructed upon assumptions on the type of microstructures that form in such alloys. In particular, we consider lamination type microstructures which form by temperature- or stress-induced transformation in monoclinic and orthorhombic Copper-based alloys with cubic austenitic symmetry. Our results display a striking congruence of upper and lower bounds in the most relevant cases.     

Bounds on the effective conductivity of statistically isotropic multicomponent materials and random cell polycrystals

Upper and lower bounds on the effective conductivity of statistically isotropic multicomponent materials in d dimensions (d=2 or 3) are constructed from the minimum energy principles and appropriate trial fields. The trial fields, involving harmonic potentials and free parameters to be optimized, lead to the bounds containing up to three-point correlation information about the microgeometry of a composite. The bounds are applied to give estimates for the symmetric cell materials, which are optimal over some ranges of parameters, and asymmetric multicoated spheres, which yield the exact effective conductivity in certain cases. The results also agree with many known ones. New bounds for random cell polycrystals are obtained and illustrated on a number of polycrystalline aggregates.     

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.     

Microstructures réalisant toute compressibilité ou conductivité encadrée par les bornes de Hashin et Shtrikman

The conductivity or bulk modulus of an isotropic composite consisting of two well-ordered isotropic phases is situated between the lower and upper bounds of Hashin and Shtrikman. This work is concerned with the construction of two new microstructures capable of realizing every conductivity or bulk modulus comprised between these bounds. The two microstructures obtained, appearing to be simpler than those presented in the literature to achieve the same purpose, have the property that the variation of the corresponding conductivity or bulk modulus is strictly monotone with respect to the geometric parameter characterizing them. To cite this article: Q.-C. He, H.L. Quang, C. R. Mecanique 333 (2005).     

An upper bound to the free energy of n-variant polycrystalline shape memory alloys

In the last two decades, the problem of computing the elastic energy of phase transforming materials has been studied by a variety of research groups. Due to the non-quasiconvexity of the underlying multi-well landscape, different relaxation methods have been used in order to estimate the quasiconvex envelope of the energy density, for which no explicit expression is known at present.This paper combines a recently developed lamination bound for monocrystalline shape memory alloys which relies on martensitic twinned microstructures with the work of Smyshlyaev and Willis [1998a. A ‘non-local’ variational approach to the elastic energy minimization of martensitic polycrystals. Proc. R. Soc. London A 454, 1573–1613]. As a result, a lamination upper bound for n-variant polycrystalline martensitic materials is obtained.The lamination bound is then compared with Reuß- and Taylor-type estimates. While, for given volume fractions, good agreement of lamination upper and convexification lower bounds is obtained, a comparison using energy-minimizing volume fractions computed from the various bounds yields larger differences. Finally, we also investigate the influence of the polycrystal's texture. For a strong ellipsoidal texture, we observe even better agreement of upper and lower bounds than for the case of isotropic statistics.     

Mindlin second-gradient elastic properties from dilute two-phase Cauchy-elastic composites. Part I: Closed form expression for the effective higher-order constitutive tensor

It is shown that second-order homogenization of a Cauchy-elastic dilute suspension of randomly distributed inclusions yields an equivalent second gradient (Mindlin) elastic material. This result is valid for both plane and three-dimensional problems and extends earlier findings by Bigoni and Drugan [Bigoni, D., Drugan, W.J., 2007. Analytical derivation of Cosserat moduli via homogenization of heterogeneous elastic materials. J. Appl. Mech. 74, 741–753] from several points of view: (i) the result holds for anisotropic phases with spherical or circular ellipsoid of inertia; (ii) the displacement boundary conditions considered in the homogenization procedure is independent of the characteristics of the material; (iii) a perfect energy match is found between heterogeneous and equivalent materials (instead of an optimal bound). The constitutive higher-order tensor defining the equivalent Mindlin solid is given in a surprisingly simple formula. Applications, treatment of material symmetries and positive definiteness of the effective higher-order constitutive tensor are deferred to Part II of the present article.     

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.     

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.     

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.     

Constitutive behaviour of anisotropic materials under shock loading

Thermodynamically and mathematically consistent constitutive equations suitable for shock wave propagation in an anisotropic material are presented in this paper. Two fundamental tensors α_ij and β_ij which represent anisotropic material properties are defined and can be considered as generalisations of the Kronecker delta symbol, which plays the main role in the theory of isotropic materials. Using two fundamental tensors α_ij and β_ij, the concept of total generalised “pressure” and pressure corresponding to the thermodynamic (equation of state) response are redefined. The equation of state represents mathematical and physical generalisation of the classical Mie–Grüneisen equation of state for isotropic material and reduces to the Mie–Grüneisen equation of state in the limit of isotropy. Based on the generalised decomposition of the stress tensor, the modified equation of state for anisotropic materials, and the modified Hill criteria, combined with the associated flow rule, a system of constitutive equations suitable for shock wave propagation is formulated. The behaviour of aluminium alloy 7010-T6 under shock loading conditions is considered. A comparison of numerical simulations with existing experimental data shows good agreement of the general pulse shape, Hugoniot Elastic Limits (HELs), and Hugoniot stress levels, and suggests that the constitutive equations are performing satisfactorily. The results are presented and discussed, and future studies are outlined.     

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.