Force heterogeneities in particle assemblies: from order to disorder

Bulk modulus of soft particle assemblies under compression

Using a numerical approach based on the coupling of the discrete and finite element methods, we explore the variation of the bulk modulus K of soft particle assemblies undergoing isotropic compression. As the assemblies densify under pressure-controlled boundary conditions, we show that the non-linearities of K rapidly deviate from predictions standing on a small-strain framework or the, so-called, Equivalent Medium Theory (EMT). Using the granular stress tensor and extracting the bulk properties of single representative grains under compression, we propose a model to predict the evolution of K as a function of the sample’s solid fraction and a reference state as the applied pressure P→0. The model closely reproduces the trends observed in our numerical experiments confirming the behavior scalability of soft particle assemblies from the individual particle scale. Finally, we present the effect of the interparticle friction on K’s evolution and how our model easily adapts to such a mechanical constraint.

Microstructural analysis of sheared polydisperse polyhedral grains

This article presents an analysis of the shear strength of numerical samples composed of polyhedra presenting a grain size dispersion. Previous numerical studies using, for instance, disks, polygons, and spheres, have consistently shown that microstructural properties linked to the fabric and force transmission allow granular media to exhibit a constant shear resistance although packing fraction can dramatically change as a broader grain-size distribution is considered. To have a complete picture of such behavior, we developed a set of numerical experiments in the frame of the discrete element method to test the shear strength of polydisperse samples composed of polyhedral grains. Although the contact networks and force transmission are quite more complex for such generalized grain shape, we can verify that the shear strength independence still holds up for 3D regular polyhedra. We make a particular focus upon the role of different contact types in the assemblies and their relative contributions to the granular connectivity and sample strength. The invariance of shear strength at the macroscopic scale results deeply linked to fine compensations at the microstructural level involving geometrical and force anisotropies of the assembly.

Scaling properties of force networks for compressed particulate systems

We consider, computationally and experimentally, the scaling properties of force networks in the systems of circular particles exposed to compression in two spatial dimensions. The simulations consider polydisperse and monodisperse particles, both frictional and frictionless, and in experiments we use monodisperse and bidisperse frictional particles. While for some of the considered systems we observe consistent scaling exponents describing the behavior of the force networks, we find that this behavior is not universal. In particular, we find that frictionless systems, independently of whether they partially crystallize under compression or not, show scaling properties that are significantly different compared to the frictional disordered ones. The findings of nonuniversality are confirmed by explicitly computing fractal dimension for the considered systems. The results of the physical experiments are consistent with the results obtained in simulations of frictional disordered systems.

Evolution of force networks in dense particulate media

We discuss sets of measures with the goal of describing dynamical properties of force networks in dense particulate systems. The presented approach is based on persistent homology and allows for extracting precise, quantitative measures that describe the evolution of geometric features of the interparticle forces, without necessarily considering the details related to individual contacts between particles. The networks considered emerge from discrete element simulations of two-dimensional particulate systems consisting of compressible frictional circular disks. We quantify the evolution of the networks for slowly compressed systems undergoing jamming transition. The main findings include uncovering significant but localized changes of force networks for unjammed systems, global (systemwide) changes as the systems evolve through jamming, to be followed by significantly less dramatic evolution for the jammed states. We consider both connected components, related in a loose sense to force chains, and loops and find that both measures provide a significant insight into the evolution of force networks. In addition to normal, we consider also tangential forces between the particles and find that they evolve in the consistent manner. Consideration of both frictional and frictionless systems leads us to the conclusion that friction plays a significant role in determining the dynamical properties of the considered networks. We find that the proposed approach describes the considered networks in a precise yet tractable manner, making it possible to identify features which could be difficult or impossible to describe using other approaches.

Rheology of three-dimensional packings of aggregates: microstructure and effects of nonconvexity

We use three-dimensional contact dynamics simulations to analyze the rheological properties of granular materials composed of rigid aggregates. The aggregates are made from four overlapping spheres and described by a nonconvexity parameter depending on the relative positions of the spheres. The macroscopic and microstructural properties of several sheared packings are analyzed as a function of the degree of nonconvexity of the aggregates. We find that the internal angle of friction increases with the nonconvexity. In contrast, the packing fraction first increases to a maximum value but declines as the nonconvexity increases further. At a high level of nonconvexity, the packings are looser but show a higher shear strength. At the microscopic scale, the fabric and force anisotropy, as well as the friction mobilization, are enhanced by multiple contacts between aggregates and interlocking, thus revealings the mechanical and geometrical origins of shear strength.

Persistence of force networks in compressed granular media

We utilize the tools of persistent homology to analyze features of force networks in dense granular matter, modeled as a collection of circular, inelastic frictional particles. The proposed approach describes these networks in a precise and tractable manner, allowing us to identify features that are difficult or impossible to characterize by other means. In contrast to other techniques that consider each force threshold level separately, persistent homology allows us to consider all threshold levels at once to describe the force network in a complete and insightful manner. We consider continuously compressed system of particles characterized by varied polydispersity and friction in two spatial dimensions. We find significant differences between the force networks in these systems, suggesting that their mechanical response may differ considerably as well.

Force chains and contact network topology in sheared packings of elongated particles

By means of contact dynamic simulations, we investigate the contact network topology and force chains in two-dimensional packings of elongated particles subjected to biaxial shearing. The morphology of large packings of elongated particles in quasistatic equilibrium is complex due to the combined effects of local nematic ordering of the particles and orientations of contacts between particles. The effect of elongation on shear behavior and dilatancy was investigated in detail in a previous paper [Azéma and Radjai, Phys. Rev. E 81, 051304 (2010)]. Here, we show how particle elongation affects force distributions and force-fabric anisotropy via various local structures allowed by steric exclusions and the requirement of force balance. We find that the force distributions become increasingly broader as particles become more elongated. Interestingly, the weak force network transforms from a passive stabilizing agent with respect to strong force chains to an active force-transmitting network for the whole system. The strongest force chains are carried by side-side contacts oriented along the principal stress direction.