Elastic waves in particulate glass-rubber mixtures.

We investigate the propagation of waves in dense static granular packings made of soft and stiff particles subjected to hydrostatic stress. Physical experiments in a triaxial cell equipped with broadband piezoelectric wave transducers have been performed at ultrasound frequencies. The time of flight is measured in order to study the combined effect of applied stress and rubber content on the elastic properties of the mixtures. The bulk stiffness deduced from the wave speed is nonlinear and non-monotonic with the increasing percentage of rubber with a more prominent effect at higher pressures. Moreover, in the frequency domain, a spectral analysis gives insights on the transition from a glass- to a rubber-dominated regime and the influence of rubber particles on the energy dissipation. Mixtures with rubber content below 30% show enhanced damping properties, associated with slightly higher stiffness and lighter weight.

Collapse modes in simple cubic and body-centered cubic arrangements of elastic beads.

Collapse modes in compressed simple cubic (SC) and body-centered cubic (BCC) periodic arrangements of elastic frictionless beads were studied numerically using the discrete element method. Under pure hydrostatic compression, the SC arrangement tends to transform into a defective hexagonal close-packed or amorphous structure. The BCC assembly exhibits several modes of collapse, one of which, identified as cI16 structure, is consistent with the behavior of BCC metals Li and Na under high pressure. The presence of a deviatoric stress leads to the transformation of the BCC structure into face-centered cubic (FCC) one via the Bain path. The observed effects expand the knowledge on possible packings of soft elastic spheres and transformations between them, while providing an unexpected link with the mechanical behavior of certain atomic systems.

A master equation for the probability distribution functions of forces in soft particle packings.

We study the microscopic response of force-chain networks in jammed soft particles to quasi-static isotropic (de)compressions by molecular dynamics simulations. We show that not only contacts but also interparticle gaps between the nearest neighbors must be considered for the stochastic evolution of the probability distribution functions (PDFs) of forces, where the mutual exchange of contacts and interparticle gaps, i.e. opening and closing contacts, are also crucial to the incremental system behavior. By numerically determining the transition rates for all changes of contacts and gaps, we formulate a Master equation for the PDFs of forces, where the insight one gets from the transition rates is striking: the mean change of forces reflects non-affine system responses, while their fluctuations obey uncorrelated Gaussian statistics. In contrast, interparticle gaps react mostly affine in average, but imply multi-scale correlations according to a much wider stable distribution function.

Effect of cohesion on shear banding in quasistatic granular materials.

Cohesive powders have widely different bulk behavior due to their peculiar interactions. We use discrete element simulations to investigate the effect of contact cohesion on the steady state flow of dense powders in a slowly sheared split-bottom Couette cell, which imposes a wide stable shear band. The intensity of cohesive forces can be quantified by the granular Bond number (Bo), namely the ratio between maximum attractive force and average force due to external compression. We find that the shear banding phenomenon is almost independent of cohesion for Bond numbers Bo<1, however for Bo≥1 cohesive forces start to play an important role, as both width and center position of the band increase. Inside the shear band, the mean normal contact force is independent of cohesion and depends only on the confining stress. In contrast, when the behavior is analyzed focusing on the eigendirections of the local strain rate tensor, a dependence on cohesion shows up. Forces carried by contacts along the compressive and tensile directions are symmetric about the mean force (larger and smaller respectively), while the force along the third, neutral direction follows the mean force. This anisotropy of the force network increases with cohesion, just like the heterogeneity in all (compressive, tensile and neutral) directions.

Micro-macro correlations and anisotropy in granular assemblies under uniaxial loading and unloading.

The influence of contact friction on the behavior of dense, polydisperse granular assemblies under uniaxial (oedometric) loading and unloading deformation is studied using discrete element simulations. Even though the uniaxial deformation protocol is one of the "simplest" element tests possible, the evolution of the structural anisotropy necessitates its careful analysis and understanding, since it is the source of interesting and unexpected observations. On the macroscopic, homogenized, continuum scale, the deviatoric stress ratio and the deviatoric fabric, i.e., the microstructure behave in a different fashion during uniaxial loading and unloading. The maximal stress ratio and strain increase with increasing contact friction. In contrast, the deviatoric fabric reaches its maximum at a unique strain level independent of friction, with the maximal value decreasing with friction. For unloading, both stress and fabric respond to unloading strain with a friction-dependent delay but at different strains. On the micro-level, a friction-dependent non-symmetry of the proportion of weak (strong) and sliding (sticking) contacts with respect to the total contacts during loading and unloading is observed. Coupled to this, from the directional probability distribution, the "memory" and history-dependent behavior of granular systems is confirmed. Surprisingly, while a rank-2 tensor is sufficient to describe the evolution of the normal force directions, a sixth order harmonic approximation is necessary to describe the probability distribution of contacts, tangential force, and mobilized friction. We conclude that the simple uniaxial deformation activates microscopic phenomena not only in the active Cartesian directions, but also at intermediate orientations, with the tilt angle being dependent on friction, so that this microstructural features cause the interesting, nontrivial macroscopic behavior.

Contact anisotropy and coordination number for a granular assembly: a comparison of distinct-element-method simulations and theory.

We study an ideal granular aggregate consisting of elastic spherical particles, isotropic in stress and anisotropic in the contact network. Because of the contact anisotropy, a confining pressure applied at zero deviatoric stress, produces shear strain as well as volume strain. Our goal is to predict the coordination number k, the average number of contacts per particle, and the magnitude of the contact anisotropy ɛ, from knowledge of the elastic moduli of the aggregate. We do this through a theoretical model based upon the well known effective medium theory. However, rather than focusing on the moduli, we consider their ratios over the moduli of an equivalent isotropic state. We observe good agreement between numerical simulation and theory.