Reference: E.J. Garboczi and J.G. Berryman, Mechanics of Materials, 33, No. 8, 455-470 (2001).

(PDF Version of Original Paper)

Go back to Part II Chapter 7 Sec.(7) Go back to Table of Contents

Elastic moduli of a material containing composite inclusions: Effective medium theory and finite element computations

E.J. Garboczi
National Institute of Standards and Technology
Building Materials Division
100 Bureau Drive, Stop 8621
Gaithersburg, Maryland 20899-8621 USA


J.G. Berryman
Lawrence Livermore National Laboratory
P. O. Box 808 L-200
Livermore, CA 94551-9900

Abstract

Concrete is a good example of a composite material in which the inclusions (rocks and sand) are surrounded by a thin shell of altered matrix material and embedded in the normal matrix material. Concrete, therefore, may be viewed as consisting of a matrix material containing composite inclusions. Assigning each of these phases different linear elastic moduli results in a complicated effective elastic moduli problem.

A new kind of differential effective medium theory (D-EMT) is presented in this paper that is intended to address this problem. The key new idea is that each inclusion particle, surrounded by a shell of another phase, is mapped onto an effective particle of uniform elastic moduli. The resulting simpler composite, with a normal matrix, is then treated in usual D-EMT. Before use, however, the accuracy of this method must be determined, as effective medium theory of any kind is an uncertain approximation. One good way to assess the accuracy of effective medium theory is to compare to exact results for known microstructures and phase moduli. Exact results, however, only exist for certain microstructures (e.g., dilute limit of inclusions) or special choices of the moduli (e.g., equal shear moduli). Recently, a special finite element method has been developed that can compute the linear elastic moduli of an arbitrary digital image in 2-D or 3-D. If a random microstructure can be represented with enough resolution by a digital image, then its elastic moduli can be readily computed. This method is used, after proper error analysis, to provide stringent tests of the new D-EMT equations, which are found to compare favorably to numerically exact finite element simulations, in both 2-D and 3-D, with varying composite inclusion particle size distributions.