Internal states of model isotropic granular packings. III. Elastic properties

In this third and final paper of a series, elastic properties of numerically simulated isotropic packings of spherical beads assembled by different procedures, as described in the first companion paper, and then subjected to a varying confining pressure, as reported in the second companion paper, are investigated. In addition to the pressure, which determines the stiffness of contacts because of Hertz's law, elastic moduli are chiefly sensitive to the coordination number z, which should not be regarded as a function of the packing density. Comparisons of numerical and experimental results for glass beads in the 10 kPa-10 MPa pressure range reveal similar differences between dry samples prepared in a dense state by vibrations and lubricated packings, so that the greater stiffness of the latter, in spite of their lower density, can be attributed to a larger coordination number. Effective medium type approaches, or Voigt and Reuss bounds, provide good estimates of bulk modulus B, which can be accurately bracketed, but badly fail for shear modulus G, especially in low z configurations under low pressure. This is due to the different response of tenuous, fragile networks to changes in load direction, as compared to load intensity. In poorly coordinated packings, the shear modulus, normalized by the average contact stiffness, tends to vary proportionally to the degree of force indeterminacy per unit volume, even though this quantity does not vanish in the rigid limit. The elastic range extends to small strain intervals and compares well with experimental observations on sands. The origins of nonelastic response are discussed. We conclude that elastic moduli provide access to mechanically important information about coordination numbers, which escape direct measurement techniques, and indicate further perspectives.