The determination of the elastic properties of composite materials (multiphase aggregates, polycrystals, and porous or cracked solids) from the elastic properties of the components may be approached in several ways. The problem may be treated statistically, via scattering theory, through variational principles, or by the assumption of specific geometries for the material under consideration. Each of these methods is reviewed in turn. The widely used Voigt-Reuss-Hill average can be a poor approximation for both two-phase composites and polycrystals, and its replacement by the two Hashin-Shtrikman bounds is recommended. For pore-containing or crack-containing media, specific geometry models must be considered if useful results are to be obtained. If aggregate theory is used to estimate the moduli of individual components of a composite whose bulk properties are known, the shear moduli of the component phases must be matched (within a factor of 2 or 3) for the method to be useful. Results for nonlinear composites (which allow calculation of the pressure variation of aggregate moduli) have been obtained for only a few special cases.