Internal friction and absence of dilatancy of packings of frictionless polygons

Dilatancy, shear jamming, and a generalized jamming phase diagram of frictionless sphere packings

Granular packings display the remarkable phenomenon of dilatancy, wherein their volume increases upon shear deformation. Conventional wisdom and previous results suggest that dilatancy, also being the related phenomenon of shear-induced jamming, requires frictional interactions. Here, we show that the occurrence of isotropic jamming densities φ_j above the minimal density (or the J-point density) φ_J leads both to the emergence of shear-induced jamming and dilatancy in frictionless packings. Under constant pressure shear, the system evolves into a steady-state at sufficiently large strains, whose density only depends on the pressure and is insensitive to the initial jamming density φ_j. In the limit of vanishing pressure, the steady-state exhibits critical behavior at φ_J. While packings with different φ_j values display equivalent scaling properties under compression, they exhibit striking differences in rheological behaviour under shear. The yield stress under constant volume shear increases discontinuously with density when φ_j > φ_J, contrary to the continuous behaviour in generic packings that jam at φ_J. Our results thus lead to a more coherent, generalised picture of jamming in frictionless packings, which also have important implications on how dilatancy is understood in the context of frictional granular matter.

Transition from steady shear to oscillatory shear rheology of dense suspensions

Recent studies have highlighted that oscillatory and time-dependent shear flows might help increase the flowability of dense suspensions. While most focus has been on cross-flows we here study a simple two-dimensional suspensions where we apply simultaneously oscillatory and stationary shear along the same direction. We first show that the dissipative viscosities in this set-up significantly decrease with an increasing shear-rate magnitude of the oscillations and given that the oscillatory strain is small, in a similar fashion as found previously for cross-flow oscillations. As for cross-flow oscillations, the decrease can be attributed to the large decrease in the number of contacts and an altered microstructure as one transitions from a steady shear to an oscillatory shear dominated rheology. As subresults we find both an extension to the μ(J) rheology, a constitutive relationship between the shear stresses and the shear rate, valid for oscillatory shear flows and that shear-jamming of frictional particles at oscillatory shear dominated flows occurs at higher packing fractions compared to steady shear dominated flows.

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.

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.

The physics of jamming for granular materials: a review

Granular materials consist of macroscopic grains, interacting via contact forces, and unaffected by thermal fluctuations. They are one of a class systems that undergo jamming, i.e. a transition between fluid-like and disordered solid-like states. Roughly twenty years ago, proposals by Cates et al for the shear response of colloidal systems and by Liu and Nagel, for a universal jamming diagram in a parameter space of packing fraction, ϕ, shear stress, τ, and temperature, T raised key questions. Contemporaneously, experiments by Howell et al and numerical simulations by Radjai et al and by Luding et al helped provide a starting point to explore key insights into jamming for dry, cohesionless, granular materials. A recent experimental observation by Bi et al is that frictional granular materials have a a re-entrant region in their jamming diagram. In a range of ϕ, applying shear strain, γ, from an initially force/stress free state leads to fragile (in the sense of Cates et al), then anisotropic shear jammed states. Shear jamming at fixed ϕ is presumably conjugate to Reynolds dilatancy, involving dilation under shear against deformable boundaries. Numerical studies by Radjai and Roux showed that Reynolds dilatancy does not occur for frictionless systems. Recent numerical studies by several groups show that shear jamming occurs for finite, but not infinite, systems of frictionless grains. Shear jamming does not lead to known ordering in position space, but Sarkar et al showed that ordering occurs in a space of force tiles. Experimental studies seeking to understand random loose and random close packings (rlp and rcp) and dating back to Bernal have probed granular packings and their response to shear and intruder motion. These studies suggest that rlp's are anisotropic and shear-jammed-like, whereas rcp's are likely isotropically jammed states. Jammed states are inherently static, but the jamming diagram may provide a context for understanding rheology, i.e. dynamic shear in a variety of systems that include granular materials and suspensions.

Microscopic Origins of Shear Jamming for 2D Frictional Grains

Shear jamming (SJ) occurs for frictional granular materials with packing fractions ϕ in ϕ_{S}<ϕ<ϕ_{J}^{0}, when the material is subject to shear strain γ starting from a force-free state. Here, ϕ_{J}^{μ} is the isotropic jamming point for particles with a friction coefficient μ. SJ states have mechanically stable anisotropic force networks, e.g., force chains. Here, we investigate the origins of SJ by considering small-scale structures-trimers and branches-whose response to shear leads to SJ. Trimers are any three grains where the two outer grains contact a center one. Branches occur where three or more quasilinear force chain segments intersect. Certain trimers respond to shear by compressing and bending; bending is a nonlinear symmetry-breaking process that can push particles in the dilation direction faster than the affine dilation. We identify these structures in physical experiments on systems of two-dimensional frictional discs, and verify their role in SJ. Trimer bending and branch creation both increase Z above Z_{iso}≃3 needed for jamming 2D frictional grains, and grow the strong force network, leading to SJ.

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.

Microstructure as a function of the grain size distribution for packings of frictionless disks: Effects of the size span and the shape of the distribution

This article presents a numerical study of the effects of grain size distribution (GSD) on the microstructure of two-dimensional packings of frictionless disks. The GSD is described by a power law with two parameters controlling the size span and the shape of the distribution. First, several samples are built for each combination of these parameters. Then, by means of contact dynamics simulations, the samples are densified in oedometric conditions and sheared in a simple shear configuration. The microstructure is analyzed in terms of packing fraction, local ordering, connectivity, and force transmission properties. It is shown that the microstructure is notoriously affected by both the size span and the shape of the GSD. These findings confirm recent observations regarding the size span of the GSD and extend previous works by describing the effects of the GSD shape. Specifically, we find that if the GSD shape is varied by increasing the proportion of small grains by a certain amount, it is possible to increase the packing fraction, increase coordination, and decrease the proportion of floating particles. Thus, by carefully controlling the GSD shape, it is possible to obtain systems that are denser and better connected, probably increasing the system's robustness and optimizing important strength properties such as stiffness, cohesion, and fragmentation susceptibility.

Effects of grain size distribution on the packing fraction and shear strength of frictionless disk packings

Using discrete element methods, the effects of the grain size distribution on the density and the shear strength of frictionless disk packings are analyzed. Specifically, two recent findings on the relationship between the system's grain size distribution and its rheology are revisited, and their validity is tested across a broader range of distributions than what has been used in previous studies. First, the effects of the distribution on the solid fraction are explored. It is found that the distribution that produces the densest packing is not the uniform distribution by volume fractions as suggested in a recent publication. In fact, the maximal packing fraction is obtained when the grading curve follows a power law with an exponent close to 0.5 as suggested by Fuller and Thompson in 1907 and 1919 [Trans Am. Soc. Civ. Eng. 59, 1 (1907) and A Treatise on Concrete, Plain and Reinforced (1919), respectively] while studying mixtures of cement and stone aggregates. Second, the effects of the distribution on the shear strength are analyzed. It is confirmed that these systems exhibit a small shear strength, even if composed of frictionless particles as has been shown recently in several works. It is also found that this shear strength is independent of the grain size distribution. This counterintuitive result has previously been shown for the uniform distribution by volume fractions. In this paper, it is shown that this observation keeps true for different shapes of the grain size distribution.

Flow of wet granular materials: A numerical study

We simulate dense assemblies of frictional spherical grains in steady shear flow under controlled normal stress P in the presence of a small amount of an interstitial liquid, which gives rise to capillary menisci, assumed isolated (pendular regime), and attractive forces, which are hysteretic: Menisci form at contact, but do not break until grains are separated by a finite rupture distance. The system behavior depends on two dimensionless control parameters, inertial number I and reduced pressure P*=aP/(πΓ), comparing confining forces ∼a2P to meniscus tensile strength F0=πΓa, for grains of diameter a joined by menisci with surface tension Γ. We pay special attention to the quasistatic limit of slow flow and observe systematic, enduring strain localization in some of the cohesion-dominated (P*∼0.1) systems. Homogeneous steady flows are characterized by the dependence of internal friction coefficient μ* and solid fraction Φ on I and P*. We also record normal stress differences, fairly small but not negligible and increasing for decreasing P*. The system rheology is moderately sensitive to saturation within the pendular regime, but would be different in the absence of capillary hysteresis. Capillary forces have a significant effect on the macroscopic behavior of the system, up to P* values of several units, especially for longer force ranges associated with larger menisci. The concept of effective pressure may be used to predict an order of magnitude for the strong increase of μ* as P* decreases but such a crude approach is unable to account for the complex structural changes induced by capillary cohesion, with a significant decrease of Φ and different agglomeration states and anisotropic fabric. Likewise, the Mohr-Coulomb criterion for pressure-dependent critical states is, at best, an approximation valid within a restricted range of pressures, with P*≥1. At small enough P*, large clusters of interacting grains form in slow flows, in which liquid bonds survive shear strains of several units. This affects the anisotropies associated with different interactions and the shape of function μ*(I), which departs more slowly from its quasistatic limit than in cohesionless systems (possibly explaining the shear banding tendency).