Geometric cohesion in two-dimensional systems composed of star-shaped particles.

Using a discrete element method, we investigate the phenomenon of geometric cohesion in granular systems composed of star-shaped particles with 3 to 13 arms. This was done by analyzing the stability of columns built with these particles and by studying the microstructure of these columns in terms of density and connectivity. We find that systems composed of star-shaped particles can exhibit geometric cohesion (i.e., a solidlike behavior, in the absence of adhesive forces between the grains), depending on the shape of the particles and the friction between them. This phenomenon is observed up to a given critical size of the system, from which a transition to a metastable behavior takes place. We also have evidence that geometric cohesion is closely linked to the systems' connectivity and especially to the capability of forming interlocked interactions (i.e., multicontact interactions that hinder the relative rotation of the grains). Our results contribute to the understanding of the interesting and potentially useful phenomenon of geometric cohesion. In addition, our work supplements an important set of experimental observations and sheds light on the complex behavior of real, three-dimensional, granular systems.

Nonlinear effect of grain elongation on the flow rate in silo discharge.

By means of two-dimensional numerical simulations based on contact dynamics, we present a systematic analysis of the joint effects of grain shape (i.e., grain elongation) and system size on silo discharge for increasing orifice sizes D. Grains are rounded-cap rectangles whose aspect ratio are varied from 1 (disks) to 7. In order to clearly isolate the effect of grain shape, the mass of the grains is keeping constant as well as the condition of the discharge by reintroducing the exiting grains at the top of the silo. In order to quantify the possible size effects, the thickness W of the silos is varied from 7 to 70 grains diameter, while keeping the silos aspect ratio always equal to 2. We find that, as long as size effects are negligible, the flow rate Q increases as a Beverloo-like function with D, also for the most elongated grains. In contrast, the effects of grain elongation on the flow rate depend on orifice size. For small normalized orifice sizes, the flow rate is nearly independent with grain elongation. For intermediate normalized orifice sizes the flow rate first increases with grain elongation up to a maximum value that depends on the normalized size of the orifice and saturates as the grains become more elongated. For larger normalized orifice size, the flow rate is an increasing function of grains' aspect ratio. Velocity profiles and packing fraction profiles close to the orifice turn out to be self-similar for all grain shapes and for the whole range of orifice and system sizes studied. Following the methodology introduced by Janda et al. [Phys. Rev. Lett. 108, 248001 (2012)PRLTAO0031-900710.1103/PhysRevLett.108.248001], we explain the nonlinear variation of Q with grain elongation, and for all orifice sizes, from compensation mechanisms between the velocity and packing fraction measured at the center of the orifice. Finally, an equation to predict the evolution of Q as a function of the aspect ratio of the grains is deduced.

Compacting an assembly of soft balls far beyond the jammed state: Insights from three-dimensional imaging.

Very soft grain assemblies have unique shape-changing capabilities that allow them to be compressed far beyond the rigid jammed state by filling void spaces more effectively. However, accurately following the formation of these systems by monitoring the creation of new contacts, monitoring the changes in grain shape, and measuring grain-scale stresses is challenging. We developed an experimental method that overcomes these challenges and connects their microscale behavior to their macroscopic response. By tracking the local strain energy during compression, we reveal a transition from granular-like to continuous-like material. Mean contact geometry is shown to vary linearly with the packing fraction, which is supported by a mean field approximation. We also validate a theoretical framework which describes the compaction from a local view. Our experimental framework provides insights into the granular micromechanisms and opens perspectives for rheological analysis of highly deformable grain assemblies in various fields ranging from biology to engineering.

Grain size distribution does not affect the residual shear strength of granular materials: An experimental proof.

Granular materials are used in several fields and in a wide variety of processes. An important feature of these materials is the diversity of grain sizes, commonly referred to as polydispersity. When granular materials are sheared, they exhibit a predominant small elastic range. Then, the material yields, with or without a peak shear strength depending on the initial density. Finally, the material reaches a stationary state, in which it deforms at a constant shear stress, which can be linked to the residual friction angle ϕ_{r}. However, the role of polydispersity on the shear strength of granular materials is still a matter of debate. In particular, a series of investigations have proved, using numerical simulations, that ϕ_{r} is independent of polydispersity. This counterintuitive observation remains elusive to experimentalists, and especially for some technical communities that use ϕ_{r} as a design parameter (e.g., the soil mechanics community). In this Letter, we studied experimentally the effects of polydispersity on ϕ_{r}. In order to do so, we built samples of ceramic beads and then sheared these samples in a triaxial apparatus. We varied polydispersity, building monodisperse, bidisperse, and polydisperse granular samples; this allowed us to study the effects of grain size, size span, and grain size distribution on ϕ_{r}. We find that ϕ_{r} is indeed independent of polydispersity, confirming the previous findings achieved through numerical simulations. Our work fairly closes the gap of knowledge between experiments and simulations.

Experimental validation of a micromechanically based compaction law for mixtures of soft and hard grains.

In this Letter, we report on an experimental study which analyzes the compressive behavior of two-dimensional bidisperse granular assemblies made of soft (hyperelastic) and hard grains in varying proportions (κ is the portion of soft grains). By means of a recently developed uniaxial compression setup [Vu and Barés, Phys. Rev. E 100, 042907 (2019)]2470-004510.1103/PhysRevE.100.042907 and using an advanced digital image correlation method, we follow, beyond the jamming point, the evolution of the main mechanical observables, from the global scale down to the strain field inside each deformable grain. First, we validate experimentally and extend to the uniaxial case a recently proposed micromechanical compaction model linking the evolution of the applied pressure P to the packing fraction ϕ [Cantor et al., Phys. Rev. Lett. 124, 208003 (2020)]0031-900710.1103/PhysRevLett.124.208003. Second, we reveal two different linear regimes depending on whether the system is above or below a crossover strain unraveling a transition from a discrete to a continuous-like system. Third, the evolution of these linear laws is found to vary linearly with κ. These results provide a comprehensive experimental and theoretical framework that can now be extended to a more general class of polydisperse soft granular systems.

Three-dimensional compaction of soft granular packings.

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 using the evolution of the packing fraction, the coordination number, and the von Misses stress distribution within the particles as the confining stress increases. The packing fraction increases and tends toward a maximum value close to 1, and the mean coordination number increases as a square root of the packing fraction. As the confining stress increases, a transition is observed from a granular-like material with exponential tails of the shear stress distributions to a continuous-like material characterized by Gaussian-like distributions of the shear stresses. We develop an equation that describes the evolution of the packing fraction as a function of the applied pressure. This equation, based on the micromechanical expression of the granular stress tensor, the limit of the Hertz contact law for small deformation, and the power-law relation between the packing fraction and the coordination of the particles, provides good predictions from the jamming point up to very high densities without the need for tuning any parameters.

Micromechanical description of the compaction of soft pentagon assemblies.

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.

Compaction of mixtures of rigid and highly deformable particles: A micromechanical model.

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.

Combined effects of contact friction and particle shape on strength properties and microstructure of sheared granular media.

We present a systematic numerical investigation concerning the combined effects of sliding friction and particle shape (i.e., angularity) parameters on the shear strength and microstructure of granular packings. Sliding friction at contacts varied from 0 (frictionless particles) to 0.7, and the particles were irregular polygons with an increasing number of sides, ranging from triangles to disks. We find that the effect of local friction on shear strength follows the same trend for all shapes. Strength first increases with local friction and then saturates at a shape-dependent value. In contrast, the effect of angularity varies, depending on the level of sliding friction. For low friction values (i.e., under 0.3), the strength first increases with angularity and then declines for the most angular shapes. For high friction values, strength systematically increases with angularity. At the microscale, we focus on the connectivity and texture of the contact and force networks. In general terms, increasing local friction causes these networks to be less connected and more anisotropic. In contrast, increasing particle angularity may change the network topology in different directions, directly affecting the macroscopic shear strength. These analyses and data constitute a first step toward understanding the joint effect of local variables such as friction and grain shape on the macroscopic rheology of granular systems.

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.

Compaction Model for Highly Deformable Particle Assemblies.

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.

Effects of particle shape mixture on strength and structure of sheared granular materials.

Using bi-dimensional discrete element simulations, the shear strength and microstructure of granular mixtures composed of particles of different shapes are systematically analyzed as a function of the proportion of grains of a given number of sides and the combination of different shapes (species) in one sample. We varied the angularity of the particles by varying the number of sides of the polygons from 3 (triangles) up to 20 (icosagons) and disks. The samples analyzed were built keeping in mind the following cases: (1) increase of angularity and species starting from disks; (2) decrease of angularity and increase of species starting from triangles; (3) random angularity and increase of species starting from disks and from polygons. The results show that the shear strength vary monotonically with increasing numbers of species (it may increase or decrease), even in the random mixtures (case 3). At the micro-scale, the variation in shear strength as a function of the number of species is due to different mechanisms depending on the cases analyzed. It may result from the increase of both the geometrical and force anisotropies, from only a decrease of frictional anisotropy, or from compensation mechanisms involving geometrical and force anisotropies.

Rheology of granular materials composed of crushable particles.

We investigate sheared granular materials composed of crushable particles by means of contact dynamics simulations and the bonded-cell model for particle breakage. Each particle is paved by irregular cells interacting via cohesive forces. In each simulation, the ratio of the internal cohesion of particles to the confining pressure, the relative cohesion, is kept constant and the packing is subjected to biaxial shearing. The particles can break into two or more fragments when the internal cohesive forces are overcome by the action of compressive force chains between particles. The particle size distribution evolves during shear as the particles continue to break. We find that the breakage process is highly inhomogeneous both in the fragment sizes and their locations inside the packing. In particular, a number of large particles never break whereas a large number of particles are fully shattered. As a result, the packing keeps the memory of its initial particle size distribution, whereas a power-law distribution is observed for particles of intermediate size due to consecutive fragmentation events whereby the memory of the initial state is lost. Due to growing polydispersity, dense shear bands are formed inside the packings and the usual dilatant behavior is reduced or cancelled. Hence, the stress-strain curve no longer passes through a peak stress, and a progressive monotonic evolution towards a pseudo-steady state is observed instead. We find that the crushing rate is controlled by the confining pressure. We also show that the shear strength of the packing is well expressed in terms of contact anisotropies and force anisotropies. The force anisotropy increases while the contact orientation anisotropy declines for increasing internal cohesion of the particles. These two effects compensate each other so that the shear strength is nearly independent of the internal cohesion of particles.

Inertial shear flow of assemblies of frictionless polygons: Rheology and microstructure.

Motivated by the understanding of shape effects in granular materials, we numerically investigate the macroscopic and microstructural properties of anisotropic dense assemblies of frictionless polydisperse rigid pentagons in shear flow, and compare them with similar systems of disks. Once subjected to large cumulative shear strains their rheology and microstructure are investigated in uniform steady states, depending on inertial number I, which ranges from the quasistatic limit ([Formula: see text]) to 0.2. In the quasistatic limit both systems are devoid of Reynolds dilatancy, i.e., flow at their random close packing density. Both macroscopic friction angle [Formula: see text], an increasing function of I , and solid fraction [Formula: see text], a decreasing function of I, are larger with pentagons than with disks at small I, but the differences decline for larger I and, remarkably, nearly vanish for [Formula: see text]. Under growing I , the depletion of contact networks is considerably slower with pentagons, in which increasingly anisotropic, but still well-connected force-transmitting structures are maintained throughout the studied range. Whereas contact anisotropy and force anisotropy contribute nearly equally to the shear strength in disk assemblies, the latter effect dominates with pentagons at small I, while the former takes over for I of the order of 10-2. The size of clusters of grains in side-to-side contact, typically comprising more than 10 pentagons in the quasistatic limit, very gradually decreases for growing I.

Shear strength and microstructure of polydisperse packings: The effect of size span and shape of particle size distribution.

By means of extensive contact dynamics simulations, we analyzed the effect of particle size distribution (PSD) on the strength and microstructure of sheared granular materials composed of frictional disks. The PSDs are built by means of a normalized ß function, which allows the systematic investigation of the effects of both, the size span (from almost monodisperse to highly polydisperse) and the shape of the PSD (from linear to pronouncedly curved). We show that the shear strength is independent of the size span, which substantiates previous results obtained for uniform distributions by packing fraction. Notably, the shear strength is also independent of the shape of the PSD, as shown previously for systems composed of frictionless disks. In contrast, the packing fraction increases with the size span, but decreases with more pronounced PSD curvature. At the microscale, we analyzed the connectivity and anisotropies of the contacts and forces networks. We show that the invariance of the shear strength with the PSD is due to a compensation mechanism which involves both geometrical sources of anisotropy. In particular, contact orientation anisotropy decreases with the size span and increases with PSD curvature, while the branch length anisotropy behaves inversely.

Binary mixtures of disks and elongated particles: Texture and mechanical properties.

We analyze the shear strength and microstructure of binary granular mixtures consisting of disks and elongated particles by varying systematically both the mixture ratio and degree of homogeneity (from homogeneous to fully segregated). The contact dynamics method is used for numerical simulations with rigid particles interacting by frictional contacts. A counterintuitive finding of this work is that the shear strength, packing fraction, and, at the microscopic scale, the fabric, force, and friction anisotropies of the contact network are all nearly independent of the degree of homogeneity. In other words, homogeneous mixtures have the same strength properties as segregated packings of the two particle shapes. In contrast, the shear strength increases with the proportion of elongated particles correlatively with the increase of the corresponding force and fabric anisotropies. By a detailed analysis of the contact network topology, we show that various contact types contribute differently to force transmission and friction mobilization.

Effects of shape and size polydispersity on strength properties of granular materials.

By means of extensive contact dynamics simulations, we analyze the combined effects of polydispersity both in particle size and in particle shape, defined as the degree of shape irregularity, on the shear strength and microstructure of sheared granular materials composed of pentagonal particles. We find that the shear strength is independent of the size span, but unexpectedly, it declines with increasing shape polydispersity. At the same time, the solid fraction is an increasing function of both the size span and the shape polydispersity. Hence, the densest and loosest packings have the same shear strength. At the scale of the particles and their contacts, we analyze the connectivity of particles, force transmission, and friction mobilization as well as their anisotropies. We show that stronger forces are carried by larger particles and propped by an increasing number of small particles. The independence of shear strength with regard to size span is shown to be a consequence of contact network self-organization, with the falloff of contact anisotropy compensated by increasing force anisotropy.

Bonded-cell model for particle fracture.

Particle degradation and fracture play an important role in natural granular flows and in many applications of granular materials. We analyze the fracture properties of two-dimensional disklike particles modeled as aggregates of rigid cells bonded along their sides by a cohesive Mohr-Coulomb law and simulated by the contact dynamics method. We show that the compressive strength scales with tensile strength between cells but depends also on the friction coefficient and a parameter describing cell shape distribution. The statistical scatter of compressive strength is well described by the Weibull distribution function with a shape parameter varying from 6 to 10 depending on cell shape distribution. We show that this distribution may be understood in terms of percolating critical intercellular contacts. We propose a random-walk model of critical contacts that leads to particle size dependence of the compressive strength in good agreement with our simulation data.

Internal friction and absence of dilatancy of packings of frictionless polygons.

By means of numerical simulations, we show that assemblies of frictionless rigid pentagons in slow shear flow possess an internal friction coefficient (equal to 0.183±0.008 with our choice of moderately polydisperse grains) but no macroscopic dilatancy. In other words, despite side-side contacts tending to hinder relative particle rotations, the solid fraction under quasistatic shear coincides with that of isotropic random close packings of pentagonal particles. Properties of polygonal grains are thus similar to those of disks in that respect. We argue that continuous reshuffling of the force-bearing network leads to frequent collapsing events at the microscale, thereby causing the macroscopic dilatancy to vanish. Despite such rearrangements, the shear flow favors an anisotropic structure that is at the origin of the ability of the system to sustain shear stress.

Particle alignment and clustering in sheared granular materials composed of platy particles.

By means of molecular dynamics simulations, we investigate the texture and local ordering in sheared packings composed of cohesionless platy particles. The morphology of large packings of platy particles in quasistatic equilibrium is complex due to the combined effects of local nematic ordering of the particles and anisotropic orientations of contacts between particles. We find that particle alignment is strongly enhanced by the degree of platyness and leads to the formation of face-connected clusters of exponentially decaying size. Interestingly, due to dynamics in continuous shearing, this ordering phenomenon emerges even in systems composed of particles of very low platyness differing only slightly from spherical shape. The number of clusters is an increasing function of platyness. However, at high platyness the proportion of face-face interactions is too low to allow for their percolation throughout the system.