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.     

Tailored heterogeneity increases overall stability regime of composites having a negative-stiffness inclusion

Recent theoretical and experimental results have shown the possibility of enormous increases in composite material overall elastic stiffness, damping, thermal expansion, piezoelectricity, etc., when the composite contains a tuned non-positive-definite (i.e., negative stiffness) constituent. For such composite materials to have practical utility, they must be stable. Recent research has shown they can be, for a limited range of constituent negative stiffness. This research has treated linear elastic composite materials with homogeneous phases, via the energy method and full dynamic stability analyses.In the present work, we first show how to analyze the composites previously treated by the comprehensive but simpler static stability approach, obtaining closed-form results. We then employ this approach to show that permitting heterogeneity of the positive-definite phase can substantially increase the range of constituent negative stiffness while maintaining overall composite stability. We first treat the positive-definite phase heterogeneity as piecewise homogeneous, and then treat it as continuously-varying. In the continuously-varying heterogeneity case, we seek the radially optimal distribution of the elastic moduli in the coatings, under constant coating average moduli constraint, to permit the most negative possible inclusion stiffness while maintaining overall composite stability. This is accomplished for three coating cases: constant bulk modulus but arbitrarily radially-varying shear modulus; constant shear modulus but arbitrarily radially-varying bulk modulus; and both moduli arbitrarily radially varying. We find the optimal coatings to be: a heterogeneous one with shear modulus being a specific continuously decreasing function of radius for the first case; a homogeneous one for the second case; and a heterogeneous one with both moduli being either Dirac-delta or Heaviside-step decreasing functions of radius for the last case (if the coating moduli are unrestricted in magnitude or have upper limits, respectively). The results show a substantial increase in the permissible inclusion negative stiffness range is provided by coating heterogeneity, while maintaining overall composite stability. Such an increased range of constituent negative stiffness provides an enlarged tuning parameter range for the development of novel, high-performance composite materials.     

Material instability-induced extreme damping in composites: A computational study

We investigate the effective viscoelastic performance of particle-reinforced composite materials whose particulate phase undergoes a material instability resulting in temporarily non-positive-definite elastic moduli. Recent experiments have shown that phase transitions in geometrically-constrained composite phases (such as in particles embedded in a stiff matrix) can lead to stable non-positive-definite elastic moduli, and they hinted at strong damping increases that can be achieved from such metastable composite phases. All previous theoretical efforts to explain such phenomena have used simplistic one-dimensional models or they were based on composite bounds and specific two-phase solids. Here, we study particle–matrix composites with periodic randomized particle dispersion. A finite element discretization is used in combination with a sophisticated nonlinear solver in order to perform the numerous calculations in a feasible amount of computing time. Our computational analysis shows that stable non-positive-definite inclusion moduli can indeed lead to extreme damping increases (i.e. greatly exceeding the intrinsic damping of each composite phase) and that such extreme damping arises from a shift in microstructural mechanisms.     

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.     

Dramatically stiffer elastic composite materials due to a negative stiffness phase?

Composite materials of extremely high stiffness can be produced by employing one phase of negative stiffness. Negative stiffness entails a reversal of the usual codirectional relationship between force and displacement in deformed objects. Negative stiffness structures and materials are possible, but unstable by themselves. We argue here that composites made with a small volume fraction of negative stiffness inclusions can be stable and can have overall stiffness far higher than that of either constituent. This high composite stiffness is demonstrated via several exact solutions within linearized and also fully nonlinear elasticity, and via the overall modulus tensor estimate of a variational principle valid in this case. We provide an initial discussion of stability, and adduce experimental results which show extreme composite behavior in selected viscoelastic systems under sub-resonant sinusoidal load. Viscoelasticity is known to expand the space of stability in some cases. We have not yet proved that purely elastic composite materials of the types proposed and analyzed in this paper will be stable under static load. The concept of negative stiffness inclusions is buttressed by recent experimental studies illustrating related phenomena within the elasticity and viscoelasticity contexts.     

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.     

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.     

Size-dependent effective thermoelastic properties of nanocomposites with spherically anisotropic phases

Composites made of semi-crystalline polymers and nanoparticles have a spherulitic microstructure which can be reasonably represented by a spherically anisotropic volume element. Due to the high surface-to-volume ratio of a nanoparticle, the particle-matrix interface stress, usually neglected in determining the effective elastic moduli of particle-reinforced composites, may have a non-negligible effect. To account for the latter in estimating the effective thermoelastic properties of a composite consisting of nanoparticles embedded in a semi-crystalline polymeric matrix, this work adopts a coherent interface model for the nanoparticle-matrix interface and proposes an extended version of the classical generalized-self consistent method. In particular, Eshelby's formulae widely used to calculate the elastic energy change of a homogeneous medium due to the introduction of an inhomogeneity are extended to the thermoelastic case. The nanoparticle size effect on the effective thermoelastic moduli of the composite are theoretically shown and numerically illustrated.     

Two exact micromechanics-based nonlocal constitutive equations for random linear elastic composite materials

A Hashin-Shtrikman-Willis variational principle is employed to derive two exact micromechanics-based nonlocal constitutive equations relating ensemble averages of stress and strain for two-phase, and also many types of multi-phase, random linear elastic composite materials. By exact is meant that the constitutive equations employ the complete spatially-varying ensemble-average strain field, not gradient approximations to it as were employed in the previous, related work of Drugan and Willis (J. Mech. Phys. Solids 44 (1996) 497) and Drugan (J. Mech. Phys. Solids 48 (2000) 1359) (and in other, more phenomenological works). Thus, the nonlocal constitutive equations obtained here are valid for arbitrary ensemble-average strain fields, not restricted to slowly-varying ones as is the case for gradient-approximate nonlocal constitutive equations. One approach presented shows how to solve the integral equations arising from the variational principle directly and exactly, for a special, physically reasonable choice of the homogeneous comparison material. The resulting nonlocal constitutive equation is applicable to composites of arbitrary anisotropy, and arbitrary phase contrast and volume fraction. One exact nonlocal constitutive equation derived using this approach is valid for two-phase composites having any statistically uniform distribution of phases, accounting for up through two-point statistics and arbitrary phase shape. It is also shown that the same approach can be used to derive exact nonlocal constitutive equations for a large class of composites comprised of more than two phases, still permitting arbitrary elastic anisotropy. The second approach presented employs three-dimensional Fourier transforms, resulting in a nonlocal constitutive equation valid for arbitrary choices of the comparison modulus for isotropic composites. This approach is based on use of the general representation of an isotropic fourth-rank tensor function of a vector variable, and its inverse. The exact nonlocal constitutive equations derived from these two approaches are applied to some example cases, directly rationalizing some recently-obtained numerical simulation results and assessing the accuracy of previous results based on gradient-approximate nonlocal constitutive equations.     

The total creep of viscoelastic composites under hydrostatic or antiplane loading

The problem of bounding the total creep (or total stress relaxation) of a composite made of two linear viscoelastic materials and subjected to a constant hydrostatic or antiplane loading is considered. It is done by coupling the immediate and the relaxed responses of the composite, which are pure elastic. The coupled bounds provide the possible range of the total deformation at infinite time as a function of the initial deformation of the composite. For antiplane shear existing bounds for coupled two-dimensional conductivity yield the required coupled bounds, and these are attained by doubly coated cylinder assemblages. The translation method is used to couple the effective bulk moduli of a viscoelastic composite at zero and infinite time. A number of microgeometries are found to attain the bulk modulus bounds. It is shown that the Hashin's composite sphere assemblage does not necessarily correspond to the maximum or minimum overall creep, although it necessarily attains the bounds for effective bulk moduli. For instance, there are cases when the doubly coated sphere microstructure or some special polycrystal arrangements attain the bounds on the total creep.     

Mechanical properties of single-walled carbon nanotube bundles as bulk materials

Single-walled carbon nanotubes (SWNTs) in crystalline bundles may exhibit a transition in which the cross-sections of tubes turn from perfectly circular to hexagonal, depending upon the tube diameter and externally applied pressure, and this structural instability leads to an abrupt change in the bulk elastic properties of SWNT bundles. This paper presents a hybrid atom/continuum model to study the bulk elastic properties of SWNT bundles, and the predicted characteristics of this structural instability agree well with the experimental observations available in the literature. Linearized bulk elastic properties of SWNT bundles with respect to a stable configuration are transversely isotropic and hence can be characterized by five independent elastic moduli. A complete set of these five moduli is predicted for the first time. It is found that the deformability of tube cross-sections play a dominant role in characterizing the transverse moduli.     

Micromechanics-based elastic damage modeling of particulate composites with weakened interfaces

A micromechanical framework is proposed to predict the effective elastic behavior and weakened interface evolution of particulate composites. The Eshelby’s tensor for an ellipsoidal inclusion with slightly weakened interface [Qu, J., 1993a. Eshelby tensor for an elastic inclusion with slightly weakened interfaces. Journal of Applied Mechanics 60 (4), 1048–1050; Qu, J., 1993b. The effect of slightly weakened interfaces on the overall elastic properties of composite materials. Mechanics of Materials 14, 269–281] is adopted to model spherical particles having imperfect interfaces in the composites and is incorporated into the micromechanical framework. Based on the Eshelby’s micromechanics, the effective elastic moduli of three-phase particulate composites are derived. A damage model is subsequently considered in accordance with the Weibull’s probabilistic function to characterize the varying probability of evolution of weakened interface between the inclusion and the matrix. The proposed micromechanical elastic damage model is applied to the uniaxial, biaxial and triaxial tensile loadings to predict the various stress–strain responses. Comparisons between the present predictions with other numerical and analytical predictions and available experimental data are conducted to assess the potential of the present framework.     

含柔性涂层的颗粒增强复合材料弹性模量估计

采用线弹簧型弱界面模型来模拟柔性涂层,研究柔性涂层对复合材料宏观弹性模量的影响。首先利用Mori-Tanaka方法和弱界面球形夹杂问题的弹性解,获得单夹杂内部的平均应力和平均应变,进而求得具有柔性涂层的复合材料的宏观弹性模量,并研究界面柔度对复合材料弹性模量的影响。     

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.     

New exact results for the effective electric, elastic, piezoelectric and other properties of composite ellipsoid assemblages

In this paper we present a unified treatment of composite ellipsoid assemblages in the setting of uncoupled phenomena like conductivity and elasticity and coupled phenomena like thermoelectricity and piezomagnetoelectricity. The building block of this microgeometry is a confocal ellipsoidal particle consisting of a (possibly void) core and a coating. All space is filled up with such units which have different sizes but possess the same aspect ratios. The confocal ellipsoids may have the same orientation in space or may be randomly oriented. The resulting microgeometry simulates two-phase composites in which the reinforcing components are short fibers or elongated particles. Our main interest is in obtaining information of an exact nature on the effective moduli of this microgeometry whose effective tensor symmetry structure depends on the packing mode of the coated ellipsoids. This information will sometimes be complete like the full effective thermoelectric tensor of an assemblage which contains aligned ellipsoids in which the coating is isotropic and the core is arbitrarily anisotropic. In the majority of the cases however the maximum achievable exact information will be only partial and will appear in the form of certain exact relations between the effective moduli of the microgeometry. These exact relations are obtained from exact solutions for the fields in the microstructure for a certain set of loading conditions. In all the considered cases an isotropic coating can be combined with a fully arbitrary core. This covers the most important physical case of anisotropic fibers in an isotropic matrix. Allowing anisotropy in the coating requires the fulfillment of certain constraint conditions between its moduli. Even though in this case the presence of such constraint conditions may render the anisotropic coating material hypothetical, the value of the derived solutions remains since they still provide benchmark comparisons for approximate and numerical treatments. The remarkable feature of the general analysis which covers all treated uncoupled and coupled phenomena is that it is developed solely on the basis of potential solutions of the conduction problem in the same microgeometry.     

Determination of elastic and plastic characteristics of TiC–TiNi alloys by the ultrasonic resonance method

Elastic characteristics and propagation velocities of ultrasonic waves in a TiC–TiNi composite material are determined by the ultrasonic resonance method. The values of the elastic moduli of the solid composite obtained are used to estimate its plastic properties. The effect of various additives on the elastic and plastic properties of the composite is studied.     

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.     

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.     

Differential geometrical method in elastic composite with imperfect interfaces

I.IntroductionWhethertheinterfacesofcompositematerialsareperfectornotwillaffectitsmacromechanicaloreffectivepropertiesimportantly.Butsofar,almostallofthestudiesontheeffectivepropertiesofcompositematerialsarebasedontheassumptionthattheinterfacesareperfectl"2].Infact,thisisnotappropriateforallinterfaces[31.Thusthestudiesonmechanicalpropertyofcompositematerialswithimperfaceintert'acehavebeenconsideredrecentlyinsomeliteratures.Hashin16]hasextendedtheelasticextremumprinciplesofminimumpotentialandm…