This article presents a model of grain fragmentation to be implemented in discrete element methods: the Split-Cell Method (SCM). In this method, the particles are of polygonal shape, and they split into polygo- nal cells once a certain failure criterion, depending on the forces exerted at the contacts, the size and shape of the grain, and the tensile strength of the material, is satisfied. The SCM is an improvement compared to other methods currently available in the literature, given that it does not restrict the shape of the grains or their fragments, mass is conserved throughout the fragmentation events, and it does not introduce artificial length scales in the system. In order to validate the proposed method, an experiment using plaster particles was conducted and its results were compared to those of a numerical simulation of the same system, finding a good match between both the experiment and the simulation.

A single circular particle subjected to diametral compression (Brazilian test)

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.

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.

Polydisperse in size and Shape

We present a detailed analysis of the morphology of granular systems composed of frictionless pentagonal particles by varying systematically both the size span and particle shape irregularity, which represent two polydispersity parameters of the system. The microstructure is characterized in terms of various statistical descriptors such as global and local packing fractions, radial distribution functions, coordination number, and fraction of floating particles. We find that the packing fraction increases with the two parameters of polydispersity, but the effect of shape polydispersity for all the investigated structural properties is significant only at low size polydispersity where the positional and/or orientational ordering of the particles prevail. We focus in more detail on the class of side/side contacts, which is the interesting feature of our system as compared to a packing of disks. We show that the proportion of such contacts has weak dependence on the polydispersity parameters. The side- side contacts do not percolate but they define clusters of increasing size as a function of size polydispersity and decreasing size as a function of shape polydispersity. The clusters have anisotropic shapes but with a decreasing aspect ratio as polydispersity increases. This feature is argued to be a consequence of strong force chains (forces above the mean), which are mainly captured by side-side contacts. Finally, the force transmission is intrinsically multiscale, with a mean force increasing linearly with particle size.

Snapshot of radial forces for S5. Line thickness is proportional to the radial force. ee contacts are in green, f v contacts in red, f e contacts in blue, and ff contacts in yellow.

We present a systematic numerical investigation of the shear strength and structure of granular packings composed of irregular polyhedral particles. The angularity of the particles is varied by increasing the number of faces from 8 (octahedronlike shape) to 596. We find that the shear strength increases with angularity up to a maximum value and saturates as the particles become more angular (below 46 faces). At the same time, the packing fraction increases to a peak value but declines for more angular particles. We analyze the connectivity and anisotropy of the microstructure by considering both the contacts and branch vectors joining particle centers. The increase of the shear strength with angularity is shown to be due to a net increase of the fabric and force anisotropies but at higher particle angularity a rapid falloff of the fabric anisotropy is compensated by an increase of force anisotropy, leading thus to the saturation of shear strength.

triaxial compaction

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.

Shear test with platy particles

This is the first paper of a series devoted to the micro-mechanical modeling of clayey soils, by means of discrete element simulations. We specifically focus here on the effect of the platy shape of particles by reducing the interactions between particles to mechanical contact forces (i.e., neither electrostatic repulsion nor van der Waals forces are taken into account). The particles are three-dimensional square plates, approximated as spheropolyhedra. Several samples composed of particles of different levels of platyness (related to the ratio of length to thickness) were numerically prepared and sheared up to large deformations. We analyzed the shear strength, packing fraction, orientation of the particles, connectivity, fabric of the interactions network, and interaction forces as functions of the platyness. We find that both the mechanical behavior and microstructure are strongly dependent on the degree of platyness. The principal underlying phenomenon is the alignment of particle faces along a particular direction. This ordering phenomenon, which emerges even for shapes that deviate only slightly from that of a sphere, enhances the ability of the packing to develop an anisotropic structure leading to large shear strength, especially as a consequence of the fabric and mobilization of friction forces. Moreover, the connectivity of the packings and their packing fraction also evolve with the platyness. In particular, the packing fraction evolves in a nonmonotonic fashion, as observed in other granular materials composed of elongated or angular particles.

The side-vertex contacts are in red and side-side contacts are in green. The line thickness is proportional to the normal force.

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.