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.

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.

Finding the grain size distribution that produces the densest arrangement in frictional sphere packings: Revisiting and rediscovering the century-old Fuller and Thompson distribution.

By means of discrete-element methods, we investigate the joint effects of the grain size distribution (GSD) and contact friction on the structure of three-dimensional samples composed of spherical grains. Specifically, we compress these systems isotropically until jamming and then analyze their structure in terms of density, connectivity, coefficients of uniformity and curvature, and parameters of grading entropy. Our study focuses on power-law GSDs and particularly on the Fuller and Thompson distribution, proposed over a century ago. First, we show that, among the set of GSDs investigated, this particular distribution produces the densest and best-connected systems, falsifying a conjecture recently posed in the literature. Second, we find that the jamming packing fraction can be accurately predicted as a function of simple descriptors of the GSD, but among these descriptors the granular entropy concept proves to be the most useful. This allows for an alternative interpretation of both jamming and grading entropy concepts. Finally, we compare the Fuller and Thompson distribution with two well-known GSDs: that of the Apollonian sphere packing and that towards which granular systems evolve after intensive grain fracturing. Surprisingly, we find that these three GSDs are practically coincident in the limit of large size spans, despite having been introduced or discovered in different scientific contexts (i.e., engineering, mathematics, and earth sciences, respectively).

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.

Granular column collapse: Analysis of grain-size effects.

The column collapse experiment is a simplified version of the complex granular flows observed in both natural and industrial contexts. Due to its simple setup and rich behavior, the column collapse has been studied experimentally and numerically by several authors. The purpose of the analyses presented in this paper is to verify whether some of the results presented in these publications are affected by grain-size effects. In order to do so, we simulate two-dimensional granular columns by means of a discrete-element method, i.e., contact dynamics. Specifically, we study the influence of the grain size, as compared to the system size, on the resultant deposit geometry and collapse duration. We show that (i) some of the previously published results may be affected by grain-size effects, (ii) in order to avoid these effects, the system-size to grain-size ratio must be larger than 75 for short columns and larger than 50 for tall columns, and (iii) the quantities that are the most affected by grain-size effects are the column mobility and the collapse duration. Our findings serve as a tool for comparing results obtained by different researchers and draw guidelines on the number of grains that must be used in order to avoid grain-size effects.

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.

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.

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.

Nonlinear effects of particle shape angularity in sheared granular media.

We analyze the effects of particle shape angularity on the macroscopic shear behavior and texture of granular packings simulated by means of the contact dynamics method. The particles are regular polygons with an increasing number of sides ranging from 3 (triangles) to 60. The packings are analyzed in the steady shear state in terms of their shear strength, packing fraction, connectivity, and fabric and force anisotropies, as functions of the angularity. An interesting finding is that the shear strength increases with angularity up to a maximum value and saturates as the particles become more angular (below six sides). In contrast, the packing fraction declines towards a constant value, so that the packings of more angular particles are looser but have higher shear strength. We show that the increase of the shear strength at low angularity is due to an increase of both contact and force anisotropies and the saturation of the shear strength for higher angularities is a consequence of a rapid falloff of the contact and normal force anisotropies compensated for by an increase of the tangential force anisotropy. This transition reflects clearly the rather special geometrical properties of these highly angular shapes, implying that the stability of the packing relies strongly on the side-side contacts and the mobilization of friction forces.

Identification of rolling resistance as a shape parameter in sheared granular media.

Using contact dynamics simulations, we compare the effect of rolling resistance at the contacts in granular systems composed of disks with the effect of angularity in granular systems composed of regular polygonal particles. In simple shear conditions, we consider four aspects of the mechanical behavior of these systems in the steady state: shear strength, solid fraction, force and fabric anisotropies, and probability distribution of contact forces. Our main finding is that, based on the energy dissipation associated with relative rotation between two particles in contact, the effect of rolling resistance can explicitly be identified with that of the number of sides in a regular polygonal particle. This finding supports the use of rolling resistance as a shape parameter accounting for particle angularity and shows unambiguously that one of the main influencing factors behind the mechanical behavior of granular systems composed of noncircular particles is the partial hindrance of rotations as a result of angular particle shape.

Simulation of cemented granular materials. I. Macroscopic stress-strain response and strain localization.

This is the first of two papers investigating the mechanical response of cemented granular materials by means of contact dynamics simulations. In this paper, a two-dimensional polydisperse sample with high-void ratio is constructed and then sheared in a simple shear numerical device at different confinement levels. We study the macroscopic response of the material in terms of mean and deviatoric stresses and strains. We show that the introduction of a local force scale, i.e., the tensile strength of the cemented bonds, causes the material to behave in a rigid-plastic fashion, so that a yield surface can be easily determined. This yield surface has a concave-down shape in the mean:deviatoric stress plane and it approaches a straight line, i.e., a Coulomb strength envelope, in the limit of a very dense granular material. Beyond yielding, the cemented structure gradually degrades until the material eventually behaves as a cohesionless granular material. Strain localization is also investigated, showing that the strains concentrate in a shear band whose thickness increases with the confining stress. The void ratio inside the shear band at the steady state is shown to be a material property that depends only on contact parameters.

Simulation of cemented granular materials. II. Micromechanical description and strength mobilization at the onset of macroscopic yielding.

This is the second of two papers investigating the mechanical response of cemented granular materials by means of contact dynamics simulations. In this paper, a two-dimensional polydisperse sample with high void ratio is sheared in a load-controlled simple shear numerical device until the stress state of the sample reaches the yield stress. We first study the stress transmission properties of the granular material in terms of the fabric of different subsets of contacts characterized by the magnitude of their normal forces. This analysis highlights the existence of a peculiar force carrying structure in the cemented material, which is reminiscent of the bimodal stress transmission reported for cohesionless granular media. Then, the evolution of contact forces and torques is investigated trying to identify the micromechanical conditions that trigger macroscopic yielding. It is shown that global failure can be associated to the apparition of a group of particles whose contacts fulfill at least one of the local rupture conditions. In particular, these particles form a large region that percolates through the sample at the moment of failure, evidencing the relationship between macroscopic yielding and the emergence of large-scale correlations in the system.

Shear strength and force transmission in granular media with rolling resistance.

We investigate a class of granular materials characterized by the possibility of interlocking between the particles. The interlocking is modeled by its effect through rolling resistance depending on relative rotation and normal force at the contact points and involving a single parameter analogous to the sliding friction coefficient. The model, which is introduced in the framework of the contact dynamics method, is applied to simulate the simple shear of a large granular sample. We present a detailed analysis regarding the influence of rolling and sliding friction parameters on the macroscopic response in terms of shear strength, fabric properties, and force transmission. Interestingly, two distinct regimes can be distinguished in which the steady-state shear strength is controlled by either rolling resistance or sliding friction. The relative contributions of rolling and sliding contacts to the shear strength are consistent with the same two regimes. Interlocking strongly affects the force network by enhancing the arching effect and thus increasing the relative importance of weak contact forces and torques, which is reflected in a decreasing power-law probability distribution of the contact forces and torques below the mean. Due to the combined effect of friction and interlocking, the force-carrying backbone takes an increasingly columnar aspect involving a low fraction of particles. Our data suggest that the nature of the weak contact network is strongly affected by the formation of these columns of particles which do not need to be propped laterally. In particular, in the limit of high rolling resistance and sliding friction, the role of the weak network of contacts is no longer to prop the force chains, but, like the strong contact network, to actively sustain the deviatoric load imposed on the system.

Additive decomposition of shear strength in cohesive granular media from grain-scale interactions.

We study cemented granular media by introducing cohesive bonding (sliding or rolling friction and tensile strength) between grains in the framework of the contact dynamics method. We find that, for a wide range of bond parameters, the macroscopic angle of friction at the peak state can be split into three distinct terms of collisional, frictional and dilational origins. Remarkably, the macroscopic tensile strength depends only on the bond tensile strength, and the friction angle at the peak state is proportional to the dilatancy angle which varies linearly with sliding friction.