Simple algebraic approximations for the effective elastic moduli of cubic arrays of spheres

Recently, Cohen and Bergman (Phys. Rev. B 68 (2003a) 24104) applied the method of elastostatic resonances to the three-dimensional problem of nonoverlapping spherical isotropic inclusions arranged in a cubic array in order to calculate the effective elastic moduli. The leading order in this systematic perturbation expansion, which is related to the Clausius-Mossotti approximation of electrostatics, was obtained in the form of simple algebraic expressions for the elastic moduli. Explicit expressions were derived for the case of a simple cubic array of spheres, and comparison was made with some accurate results. Here, we present explicit expressions for the effective elastic moduli of base-centered and face-centered cubic arrays as well, and make a comparison with other estimates and with accurate numerical results. The simple algebraic expressions provide accurate results at low volume fractions of the inclusions and are good estimates at moderate volume fractions even when the contrast is high.     

Effective properties of elastic periodic composite media with fibers

A series solution to obtain the effective properties of some elastic composites media having periodically located heterogeneities is described. The method uses the classical expansion along Neuman series of the solution of the periodic elasticity problem in Fourier space, based on the Green's tensor, and exact expressions of factors depending on the shape of the inclusions. Some properties of convergence of the solution are presented, more specifically concerning the elasticity tensor of the reference medium, showing that the convergence occurs even for empty fibers. The solution is extended for rigid inclusions. A comparison is made with previous exact solutions for a fiber composite made of cylindrical fibers with circular cross-sections and with previous estimates. Different examples are presented for new situations concerning the study of fiber composites: composites with elliptic cross-sections and multi-phase fibrous composites.     

A strong contrast homogenization formulation for multi-phase anisotropic materials

Homogenization relations, linking a material's properties at the mesoscale to those at the macroscale, are fundamental tools for design and analysis of microstructure. Recent advances in this field have successfully applied spectral techniques to Kroner-type perturbation expansions for polycrystalline and composite materials to provide efficient inverse relations for materials design. These expansions have been termed ‘weak-contrast’ expansions due to the conditionally convergent integrals, and the reliance upon only small perturbations from the reference property. In 1955, Brown suggested a different expansion for electrical conductivity that resulted in absolutely convergent integrals. Torquato subsequently applied the method to elasticity, with good results even for high-contrast materials; thus it is commonly referred to as a ‘strong contrast’ expansion. The methodology has been applied to elasticity for two phases of isotropic material, generally assuming macroscopic isotropy (with noted exceptions), thus resulting in a rather elegant form of the solution.

More recently, a multi-phase form of the solution was developed for conductivity. This paper builds upon this result to apply the method to elasticity of polycrystalline materials with both local and global anisotropy. New spectral formulations are subsequently developed for both the weak and strong contrast solutions. These form the basis for efficient microstructure analysis using these frameworks, and subsequently for inverse design applications. The process is taken through to demonstration of a property closure, which acts as the basis for materials design; the closure delineates the envelope of all physically realizable property combinations for the chosen properties, based on the particular homogenization relation being used.     

Eshelby's problem of non-elliptical inclusions

The Eshelby problem consists in determining the strain field of an infinite linearly elastic homogeneous medium due to a uniform eigenstrain prescribed over a subdomain, called inclusion, of the medium. The salient feature of Eshelby's solution for an ellipsoidal inclusion is that the strain tensor field inside the latter is uniform. This uniformity has the important consequence that the solution to the fundamental problem of determination of the strain field in an infinite linearly elastic homogeneous medium containing an embedded ellipsoidal inhomogeneity and subjected to remote uniform loading can be readily deduced from Eshelby's solution for an ellipsoidal inclusion upon imposing appropriate uniform eigenstrains. Based on this result, most of the existing micromechanics schemes dedicated to estimating the effective properties of inhomogeneous materials have been nevertheless applied to a number of materials of practical interest where inhomogeneities are in reality non-ellipsoidal. Aiming to examine the validity of the ellipsoidal approximation of inhomogeneities underlying various micromechanics schemes, we first derive a new boundary integral expression for calculating Eshelby's tensor field (ETF) in the context of two-dimensional isotropic elasticity. The simple and compact structure of the new boundary integral expression leads us to obtain the explicit expressions of ETF and its average for a wide variety of non-elliptical inclusions including arbitrary polygonal ones and those characterized by the finite Laurent series. In light of these new analytical results, we show that: (i) the elliptical approximation to the average of ETF is valid for a convex non-elliptical inclusion but becomes inacceptable for a non-convex non-elliptical inclusion; (ii) in general, the Eshelby tensor field inside a non-elliptical inclusion is quite non-uniform and cannot be replaced by its average; (iii) the substitution of the generalized Eshelby tensor involved in various micromechanics schemes by the average Eshelby tensor for non-elliptical inhomogeneities is in general inadmissible.     

Homogenization for wave propagation in periodic fibre-reinforced media with complex microstructure. I—Theory

The effective elastic properties of periodic fibre-reinforced media with complex microstructure are determined by the method of asymptotic homogenization via a novel solution to the cell problem. The solution scheme is ideally suited to materials with many fibres in the periodic cell. In this first part of the paper we discuss the theory for the most general situation—N arbitrarily anisotropic fibres within the periodic cell. For ease of exposition we then restrict attention to isotropic phases which results in a monoclinic composite material with 13 effective moduli and expressions for each of these are determined. In the second part of this paper we shall discuss results for a variety of specific microstructures.     

Scaling function, anisotropy and the size of RVE in elastic random polycrystals

This article is focused on the identification of the size of the representative volume element (RVE) in linear elastic randomly structured polycrystals made up of cubic single crystals. The RVE is approached by setting up stochastic Dirichlet and Neumann boundary value problems consistent with the Hill(-Mandel) macrohomogeneity condition. Within this framework we introduce a scaling function that relates the single crystal anisotropy to the scale of observation. We derive certain exact characteristics of the scaling function and postulate others based on detailed calculations on copper, lithium, tantalum, magnesium oxide and antimony-yttrium. In deriving the above, we make use of the fact that cubic crystals and polycrystals have a uniquely determined scale-independent bulk modulus. It turns out that the scaling function is exact in the single crystal anisotropy. A methodology to develop a material selection diagram that clearly separates the microscale from the macroscale is proposed. The proposed scaling function not only bridges the length scales but also unifies the treatment of a wide spectrum of cubic crystals. Although the scope of this article is restricted to aggregates made up of cubic-shaped and cubic-symmetry single crystals, the concept of the scaling function can be generalized to other crystal shapes and classes as well as to scaling of other elastic/inelastic properties.     

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.     

Effective elastic moduli of periodic granular composite with transversely isotropic phases

We find a rigorous solution describing the macroscopically uniform stress state of a periodic granular composite with transversely isotropic phases. The structure of the composite is modeled by a cube containing a finite number of arbitrarily arranged and oriented, transversely isotropic spherical inclusions. This provides the model with a flexible means of describing the microstructure. Applying periodic vector solutions and local expansion formulas reduces the initial boundary-value problem to a system of linear algebraic equations. By averaging the solution over the unit cell, we derived exact finite expressions for the components of the effective stiffness tensor. The numerical data presented help to evaluate the efficiency of the method and the limits of applicability of available approximate theories.Translated from Prikladnaya Mekhanika, Vol. 40, No. 9, pp. 123–130, September 2004.     

Elastic stability of inhomogeneous thin plates on an elastic foundation

We study the elastic stability of infinite inhomogeneous thin plates on an elastic foundation under in-plane compression. The elastic stiffness constants depend on the coordinate variable in the thickness direction of the plate. The elastic foundation is represented as a Winkler-type model characterized by linear and nonlinear spring constants. First we derive the Föppl–von Kármán equations by taking variations of the elastic strain energy. Next we develop the linear stability analysis of the plate under uniform in-plane compression and explicitly derive the critical loads and wave numbers for particular three cases. The effects of the material inhomogeneity, material orthotropy and loading orthotropy on the critical states are examined independently. Finally, we perform a weakly nonlinear analysis of the plate at the onset of the buckling instability. With the multiple scales method, the amplitude equations for the unstable modes that provide insight into the mode type and its amplitude are derived and then the effect of the material inhomogeneity on buckling modes are evaluated qualitatively.     

Surface instability of an elastic half space with material properties varying with depth

If a body with a stiffer surface layer is loaded in compression, a surface wrinkling instability may be developed. A bifurcation analysis is presented for determining the critical load for the onset of wrinkling and the associated wavelength for materials in which the elastic modulus is an arbitrary function of depth. The analysis leads to an eigenvalue problem involving a pair of linear ordinary differential equations with variable coefficients which are discretized and solved using the finite element method.The method is validated by comparison with classical results for a uniform layer on a dissimilar substrate. Results are then given for materials with exponential and error-function gradation of elastic modulus and for a homogeneous body with thermoelastically induced compressive stresses.     

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.     

Effective elastic properties of porous materials: Homogenization schemes vs experimental data

In this paper, we focus on the prediction of elastic moduli of isotropic porous materials made of a solid matrix having a Poisson's ratio vm of 0.2. We derive simple analytical formulae for these effective moduli based on well-known Mean-Field Eshelby-based Homogenization schemes. For each scheme, we find that the normalized bulk, shear and Young's moduli are given by the same form depending only on the porosity p. The various predictions are then confronted with experimental results for the Young's modulus of expanded polystyrene (EPS) concrete. The latter can be seen as an idealized porous material since it is made of a bulk cement matrix, with Poisson's ratio 0.2, containing spherical mono dispersed EPS beads. The Differential method predictions are found to give a very good agreement with experimental results. Thus, we conclude that when vm=0.2, the normalized effective bulk, shear and Young's modulus of isotropic porous materials can be well predicted by the simple form (1 − p)² for a large range of porosity p ranging between 0 and 0.56.     

General integral equations of thermoelasticity in micromechanics of composites with imperfectly bonded interfaces

One considers a linear thermoelastic composite medium, which consists of a homogeneous matrix containing a statistically inhomogeneous random set of heterogeneities with various interface effects and subjected to essentially inhomogeneous loading by the fields of the stresses, temperature, and body forces (e.g., for a centrifugal load). The general integral equations connecting the stress and strain fields in the point being considered and the surrounding points are obtained for the random and deterministic fields of inclusions. The method is based on a centering procedure of subtraction from both sides of a new initial integral equation their statistical averages obtained without any auxiliary assumptions such as the effective field hypothesis (EFH), which is implicitly exploited in the known centering methods. The new initial integral equation is presented in a general form of perturbations introduced by the heterogeneities and taking into account both the spring-layer model and coherent imperfect one. Some particular cases, asymptotic representations, and simplifications of proposed equations as well as a model example demonstrating the essence of two-step statistical average scheme are considered. General integral equations for the doubly and triply periodical structure composites are also obtained.     

On the thermodynamic consistency of the equivalence principle in continuum damage mechanics

We consider theories of continuum damage mechanics involving damage effect variables of different tensorial ranks. It turns out that orthotropic damage together with the use of Lemaitre's equivalence principle for the elastic part does not allow thermodynamic potentials such as the free enthalpy to exist. As the existence of these potentials is, however, a strict thermodynamic requirement, a theory employing orthotropic damage in this way is inconsistent. We show that the use of a rank-4 damage effect variable allows a consistent use of the equivalence principle.     

Fundamental solutions in antiplane elastodynamic problem for anisotropic medium under moving oscillating source

In this paper we study the antiplane problem of concentrated point force moving with constant velocity and oscillating with constant frequency in unbounded homogeneous anisotropic elastic medium.The explicit representation of the elastodynamic Green's function is obtained by using Fourier integral transform techniques for all rates of source motion as a sum of the integrals over the finite interval. The dynamic and quasistatic components of the Green's function are extracted. The stationary phase method is applied to derive an asymptotic approximation at the far wave field. The simple formulae for Poynting energy flux vectors for moving and fixed observers are presented too.It is shown that the motion brings some differences in the far field properties, such as, for example, fast and slow waves appearance under superseismic motion and modification of the wave propagation zones and their numbers.The case of isotropic medium is considered separately. For isotropic material all main formulae are obtained in explicit forms.     

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.     

The indentation modulus of elastically anisotropic materials for indenters of arbitrary shape

The contact of an indenter of arbitrary shape on an elastically anisotropic half space is considered. It is demonstrated in a theorem that the solution of the contact problem is the one that maximizes the load on the indenter for a given indentation depth. The theorem can be used to derive the best approximate solution in the Rayleigh-Ritz sense if the contact area is a priori assumed to have a certain shape. This approach is used to analyze the contact of a sphere and an axisymmetric cone on an anisotropic half space. The contact area is assumed to be elliptical, which is exact for the sphere and an approximation for the cone. It is further shown that the contact area is exactly elliptical even for conical indenters when a limited class of Green's functions is considered. If only the first term of the surface Green's function Fourier expansion is retained in the solution of the axisymmetric contact problem, a simpler solution is obtained, referred to as the equivalent isotropic solution. For most anisotropic materials, the contact stiffness determined using this approach is very close to the value obtained for both conical and spherical indenters by means of the theorem. Therefore, it is suggested that the equivalent isotropic solution provides a quick and efficient estimate for quantities such as the elastic compliance or stiffness of the contact. The “equivalent indentation modulus”, which depends on material and orientation, is computed for sapphire and diamond single crystals.     

An energy-based anisotropic yield criterion for cellular solids and validation by biaxial FE simulations

The initial yield surface of 2D lattice materials is investigated under biaxial loading using finite element analyses as well as by analytical means. The sensitivity of initial yield surface to the dominant deformation mode is explored by using both low- and high-connectivity topologies whose dominant deformation mode is either local bending or strut stretching, respectively. The effect of microstructural irregularity on the initial yield surface is also examined for both topologies. A pressure-dependent anisotropic yield criterion, which is based on total elastic strain energy density, is proposed for 2D lattice structures, which can be easily extended for application to 3D cellular solids. Proposed criterion uses elastic constants and yield strengths under uniaxial loading, and does not rely on any arbitrary parameter. The analytical framework developed allows the introduction of new scalar measures of characteristic stresses and strains that are capable of representing the elastic response of anisotropic materials with a single elastic master line under multiaxial loading.