The compaction behavior of deformable grain assemblies beyond jamming remains bewildering, and existing models that seek to find the relationship between the confining pressure P and solid fraction ϕ end up settling for empirical strategies or fitting parameters. Using a coupled discrete-finite element method, we analyze assemblies of highly deformable frictional grains under compression. We show that the solid fraction evolves nonlinearly from the jamming point and asymptotically tends to unity. Based on the micromechanical definition of the granular stress tensor, we develop a theoretical model, free from ad hoc parameters, correctly mapping the evolution of ϕ with P. Our approach unveils the fundamental features of the compaction process arising from the joint evolution of grain connectivity and the behavior of single representative grains. This theoretical framework also allows us to deduce a bulk modulus equation showing an excellent agreement with our numerical data.

This paper analyzes the compaction behavior of assemblies composed of soft (elastic) spherical particles beyond the jammed state, using three-dimensional non-smooth contact dynamic simulations. The assemblies of particles are characterized...

We analyze the isotropic compaction of mixtures composed of rigid and deformable incompressible particles by the nonsmooth contact dynamics approach. The deformable bodies are simulated using a hyperelastic neo-Hookean constitutive law by means of classical finite elements. We characterize the evolution of the packing fraction, the elastic modulus, and the connectivity as a function of the applied stresses when varying the interparticle coefficient of friction. We show first that the packing fraction increases and tends asymptotically to a maximum value φ max , which depends on both the mixture ratio and the interparticle friction. The bulk modulus is also shown to increase with the packing fraction and to diverge as it approaches φ max . From the micromechanical expression of the granular stress tensor, we develop a model to describe the compaction behavior as a function of the applied pressure, the Young modulus of the deformable particles, and the mixture ratio. A bulk equation is also derived from the compaction equation. This model lays on the characterization of a single deformable particle under compression together with a power-law relation between connectivity and packing fraction. This compaction model, set by well-defined physical quantities, results in outstanding predictions from the jamming point up to very high densities and allows us to give a direct prediction of φ max as a function of both the mixture ratio and the friction coefficient.

We analyze the isotropic compaction of assemblies composed of soft pentagons interacting through classical Coulomb friction via numerical simulations. The effect of the initial particle shape is discussed by comparing packings of pentagons with packings of soft circular particles. We characterize the evolution of the packing fraction, the elastic modulus, and the microstructure (particle rearrangement, connectivity, contact force, and particle stress distributions) as a function of the applied stresses. Both systems behave similarly: the packing fraction increases and tends asymptotically to a maximum value φ max , where the bulk modulus diverges. At the microscopic scale we show that particle rearrangements occur even beyond the jammed state, the mean coordination increases as a square root of the packing fraction, and the force and stress distributions become more homogeneous as the packing fraction increases. Soft pentagons experience larger particle rearrangements than circular particles, and such behavior decreases proportionally to the friction. Interestingly, the friction between particles also contributes to a better homogenization of the contact force network in both systems. From the expression of the granular stress tensor we develop a model that describes the compaction behavior as a function of the applied pressure, the Young modulus, and the initial shape of the particles. This model, settled on the joint evolution of the particle connectivity and the contact stress, provides outstanding predictions from the jamming point up to very high densities.