Numerical study of one-dimensional compression of granular materials. I. Stress-strain behavior, microstructure, and irreversibility

The behavior of a model granular material, made of slightly polydisperse beads with Hertz-Mindlin elastic-frictional contacts, in oedometric compression (i.e., compression along one axis, with no lateral strain) is studied by grain-level numerical simulations. We systematically investigate the influence of the (idealized) packing process on the microstructure and stresses in the initial, weakly confined equilibrium state, and prepare both isotropic and anisotropic configurations differing in solid fraction Phi and coordination number z. Phi (ranging from maximally dense to moderately loose), z (which might vary independently of Phi in dense systems), fabric and force anisotropy parameters, and the ratio K-0 of lateral stresses sigma(2) = sigma(3) to stress sigma(1) in the compression direction are monitored in oedometric compression in which sigma(1) varies by more than three orders of magnitude. K-0 reflects the anisotropy of the assembling process and may remain nearly constant in further loading if the material is already oedometrically compressed (as a granular gas) in the preparation stage. Otherwise, it tends to decrease steadily over the investigated stress range. It is related to force and fabric anisotropy parameters by a simple formula. Elastic moduli, separately computed with an appropriate matrix method, may express the response to very small stress increments about the transversely isotropic well-equilibrated states along the loading path, although oedometric compression proves an essentially anelastic process, mainly due to friction mobilization, with large irreversible effects apparent upon unloading. While the evolution of axial strain epsilon(1) and solid fraction Phi (or of the void ratio e = -1 + 1/ Phi) with axial stress sigma(1) is very nearly reversible, especially in dense samples, z is observed to decrease (as previously observed in isotropic compression) after a compression cycle if its initial value was high. K-0 relates to the evolution of internal variables and may exceed 1 in unloading. The considerably greater irreversibility of oedometric compression reported in sands, compared to our model systems, should signal contact plasticity or damage.