Rheology of granular materials composed of nonconvex particles

Rheology of nonconvex granular flows based on particle rotational characteristics

Particle shape has a profound impact on the flow behaviors of granular materials, yet effectively incorporating the role of particle shape into granular rheology remains challenging. In this study, we employ three representative types of nonconvex particles generated through the multisphere approach and identify a consistent one-to-one relationship between the rescaled friction coefficient and the inertial number I across both inertial and quasistatic flow regimes. However, variations in particle shape cause notable deviations in rheological data compared to their spherical counterparts. Based on the observed dependence of rheological data on I for various nonconvex particles and their convergence at high I, we propose an inertial number I_{s} to effectively capture the impact of particle shape on flow states. The model parameters defining I_{s} are shown to be nearly independent of flow states and configurations, with physical interpretations related to particle rotational characteristics during shear. For practical application, we propose an empirical formula to capture the dependence of model parameters on particle geometrical shapes. The robustness of the proposed model is validated by predicting flow in an inclined flow configuration and applying it to additional nonconvex particles with more irregular and asymmetric features. This establishes a crucial foundation for extending the application of this generalized rheological model to other complex granular flows.

Dynamic compaction of cohesive granular materials: scaling behavior and bonding structures

The compaction of cohesive granular materials is a common operation in powder-based manufacture of many products. However, the influence of particle-scale parameters such as bond strength on the packing structure and the general scaling of the compaction process are still poorly understood. We use particle dynamics simulations to analyze jammed configurations obtained by dynamic compaction of sticky particles under a fixed compressive pressure for a broad range of system parameter values. We show that relative porosity, representing the relative importance of porosity with respect to its minimum and maximum values, is a unique function of a modified cohesion number that combines adhesion force, confining pressure, and particle size, as well as contact stiffness, which is often assumed to be ineffective but is shown here to play an essential role in compaction. An asymmetric sigmoidal form based on two power laws provides an excellent fit to the data. The statistical properties of the bond network reveal self-balanced force structures and an exponential fall-off of the number of both tensile and compressive forces. Remarkably, the properties of the bond network depend on the cohesion number rather than the modified cohesion number, implying that similar bond network characteristics are compatible with a broad range of porosities mainly due to the effect of contact stiffness. We also discuss the origins of data points escaping the general scaling of porosity and show that they reflect either finite system size or rigid confining walls.

Contact networks and force transmission in aggregates of hexapod-shaped particles

Hexapods, consisting of three mutually orthogonal arms, have been utilized as a representative nonconvex shape to demonstrate the impact of interlocking on the strength properties of granular materials. Nevertheless, the microstructural characteristics of hexapod packings, which underlie their strength, have remained insufficiently characterized. We use particle dynamics simulations to build isotropically-packed aggregates of hexapods and we analyze the effects of aspect ratio and interparticle friction on the microstructure and force transmission. We find that the packing fraction is an unmonotonic function of aspect ratio due to competition between steric exclusions and interlocking. Interestingly, the contact coordination number declines considerably with friction coefficient, showing the stronger effect of friction on the stability of hexapod packings as compared with sphere packings. The pair distribution functions show that local ordering due to steric exclusions disappears beyond the aspect ratio 3 and the hexapods touch their second neighbors. Remarkably, hexapods of aspect ratio 3 tend to align with their neighbors and form locally ordered structures, implying a contact coordination number which is highly sensitive to the confining pressure. We also show that the probability density function of forces between hexapods is similar to that of sphere packings but with broadening exponential fall-off of strong forces as aspect ratio increases. Finally, the elastic bulk modulus of the aggregates is found to increase considerably with aspect ratio as a consequence of the rapid increase of contact density and the number of contacts with second neighbors.

Elasticity and rheology of auxetic granular metamaterials

The flowing, jamming, and avalanche behavior of granular materials is satisfyingly universal and vexingly hard to tune: A granular flow is typically intermittent and will irremediably jam if too confined. Here, we show that granular metamaterials made from particles with a negative Poisson's ratio yield more easily and flow more smoothly than ordinary granular materials. We first create a collection of auxetic grains based on a re-entrant mechanism and show that each grain exhibits a negative Poisson's ratio regardless of the direction of compression. Interestingly, we find that the elastic and yielding properties are governed by the high compressibility of granular metamaterials: At a given confinement, they exhibit lower shear modulus, lower yield stress, and more frequent, smaller avalanches than materials made from ordinary grains. We further demonstrate that granular metamaterials promote flow in more complex confined geometries, such as intruder and hopper geometries, even when the packing contains only a fraction of auxetic grains. Moreover, auxetic granular metamaterials exhibit enhanced impact absorption. Our findings blur the boundary between complex fluids and metamaterials and could help in scenarios that involve process, transport, and reconfiguration of granular materials.

Secondary flow in ensembles of nonconvex granular particles under shear

Studies of granular materials, both theoretical and experimental, are often restricted to convex grain shapes. We demonstrate that a nonconvex grain shape can lead to a qualitatively novel macroscopic dynamics. Spatial crosses (hexapods) are continuously sheared in a split-bottom container. Thereby, they develop a secondary flow profile that is completely opposite to that of rod-shaped or lentil-shaped convex grains in the same geometry. The crosses at the surface migrate towards the rotation center and sink there mimicking a "reverse Weissenberg effect." The observed surface flow field suggests the existence of a radial outward flow in the depth of the granular bed, thus, forming a convection cell. This flow field is connected with a dimple formed in the rotation center. The effect is strongly dependent on the particle geometry and the height of the granular bed.

The effect of grain shape and material on the nonlocal rheology of dense granular flows

Nonlocal rheologies allow for the modeling of granular flows from the creeping to intermediate flow regimes, using a small number of parameters. In this paper, we report on experiments testing how particle properties affect the model parameters used in the Kamrin & Koval cooperative nonlocal model, using particles of three different shapes (circles, ellipses, and pentagons) and three different materials, including one which allows for the measurement of stresses via photoelasticity. Our experiments are performed on a quasi-2D annular shear cell with a rotating inner wall and a fixed outer wall. Each type of particle is found to exhibit flows which are well-fit by nonlocal rheology, with each particle having a distinct triad of the local, nonlocal, and frictional parameters. While the local parameter b is always approximately unity, the nonlocal parameter A depends sensitively on both the particle shape and material. The critical stress ratio μ_s, above which Coulomb failure occurs, varies for particles with the same material but different shape, indicating that geometric friction can dominate over material friction.

The effect of particle shape on discharge and clogging

Granular flow is common across different fields from energy resource recovery and mineral processing to grain transport and traffic flow. Migrating particles may jam and form arches that span constrictions and hinder particle flow. Most studies have investigated the migration and clogging of spherical particles, however, natural particles are rarely spherical, but exhibit eccentricity, angularity and roughness. New experiments explore the discharge of cubes, 2D crosses, 3D crosses and spheres under dry conditions and during particle-laden fluid flow. Variables include orifice-to-particle size ratio and solidity. Cubes and 3D crosses are the most prone to clogging because of their ability to interlock or the development of face-to-face contacts that can resist torque and enhance bridging. Spheres arriving to the orifice must be correctly positioned to create stable bridges, while flat 2D crosses orient their longest axes in the direction of flowlines across the orifice and favor flow. Intermittent clogging causes kinetic retardation in particle-laden flow even in the absence of inertial effects; the gradual increase in the local particle solidity above the constriction enhances particle interactions and the probability of clogging. The discharge volume before clogging is a Poisson process for small orifice-to-particle size ratio; however, the clogging probability becomes history-dependent for non-spherical particles at large orifice-to-particle size ratio and high solidities, i.e., when particle–particle interactions and interlocking gain significance.

Link between packing morphology and the distribution of contact forces and stresses in packings of highly nonconvex particles

An external load on a particle packing is distributed internally through a heterogeneous network of particle contacts. This contact force distribution determines the stability of the particle packing and the resulting structure. Here, we investigate the homogeneity of the contact force distribution in packings of highly nonconvex particles both in two-dimensional (2D) and three-dimensional (3D) packings. A recently developed discrete element method is used to model packings of nonconvex particles of varying sphericity. Our results establish that in 3D packings the distribution of the contact forces in the normal direction becomes increasingly heterogeneous with decreasing particle sphericity. However, in 2D packings the contact force distribution is independent of particle sphericity, indicating that results obtained in 2D packings cannot be extrapolated readily to 3D packings. Radial distribution functions show that the crystallinity in 3D packings decreases with decreasing particle sphericity. We link the decreasing homogeneity of the contact force distributions to the decreasing crystallinity of 3D packings. These findings are complementary to the previously observed link between the heterogeneity of the contact force distribution and a decreasing packing crystallinity due to an increasing polydispersity of spherical particles.

Representativity of 2D Shape Parameters for Mineral Particles in Quantitative Petrography

This paper introduces an assessment of the representation of shape parameter measurements on theoretical particles. The aim of the study was to establish a numerical method for estimating sphericity, roundness, and roughness on artificially designed particles and to evaluate their interdependence. The parameters studied included a fractal dimension (FD), solidity (So), Wadell’s roundness (Rw), a perimeter-area normalized ratio (¥), and sphericity (S). The methods of the work included: (a) the design of theoretical particles with different shapes, (b) the definition of optimal analysis conditions for automated measurements, (c) the quantification of particle parameters by computer vision-based image processing, and (d) the evaluation of interdependence between the parameters. The study established the minimum sizes required for analysis of the particle shape. These varied depending on the method used (150 pixels or 50 pixels). Evaluating the relationships between the parameters showed that FD and So are independent of S. Nevertheless, Rw and ¥ are clearly dependent on S and, thus, must be numerically corrected to Rwc and ¥c. FD, So, Rwc, and ¥c were used to establish, mathematically, a new regularity parameter (RBC) that reflects the degree of roundness of a particle. The process was applied to a case study and the evaluation of all parameters corroborated previous petrographic characterizations.

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.

Shear-driven flow of athermal, frictionless, spherocylinder suspensions in two dimensions: Stress, jamming, and contacts

We use numerical simulations to study the flow of a bidisperse mixture of athermal, frictionless, soft-core two-dimensional spherocylinders driven by a uniform steady-state shear strain applied at a fixed finite rate. Energy dissipation occurs via a viscous drag with respect to a uniformly sheared host fluid, giving a simple model for flow in a non-Brownian suspension and resulting in a Newtonian rheology. We study the resulting pressure p and deviatoric shear stress σ of the interacting spherocylinders as a function of packing fraction ϕ, strain rate γ[over ̇], and a parameter α that measures the asphericity of the particles; α is varied to consider the range from nearly circular disks to elongated rods. We consider the direction of anisotropy of the stress tensor, the macroscopic friction μ=σ/p, and the divergence of the transport coefficient η_{p}=p/γ[over ̇] as ϕ is increased to the jamming transition ϕ_{J}. From a phenomenological analysis of Herschel-Bulkley rheology above jamming, we estimate ϕ_{J} as a function of asphericity α and show that the variation of ϕ_{J} with α is the main cause for differences in rheology as α is varied; when plotted as ϕ/ϕ_{J}, rheological curves for different α qualitatively agree. However, a detailed scaling analysis of the divergence of η_{p} for our most elongated particles suggests that the jamming transition of spherocylinders may be in a different universality class than that of circular disks. We also compute the number of contacts per particle Z in the system and show that the value at jamming Z_{J} is a nonmonotonic function of α that is always smaller than the isostatic value. We measure the probability distribution of contacts per unit surface length P(ϑ) at polar angle ϑ with respect to the spherocylinder spine and find that as α→0 this distribution seems to diverge at ϑ=π/2, giving a finite limiting probability for contacts on the vanishingly small flat sides of the spherocylinder. Finally, we consider the variation of the average contact force as a function of location on the particle surface.

Reduction of the bulk modulus with polydispersity in noncohesive granular solids

We study the effect of grain polydispersity on the bulk modulus in noncohesive two-dimensional granular solids. Molecular dynamics simulations in two dimensions are used to describe polydisperse samples that reach a stationary limit after a number of hysteresis cycles. For stationary samples, we obtain that the packing with the highest polydispersity has the lowest bulk modulus. We compute the correlation between normal and tangential forces with grain size using the concept of branch vector or contact length. Classifying the contact lengths and forces by their size compared to the average length and average force, respectively, we find that strong normal and tangential forces are carried by large contact lengths, generally composed of at least one large grain. This behavior is more dominant as polydispersity increases, making force networks more anisotropic and removing the support, from small grains, in the loading direction thus reducing the bulk modulus of the granular pack. Our results for two dimensions describe qualitatively the results of three-dimensional experiments.

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.

Optimizing packing fraction in granular media composed of overlapping spheres

What particle shape will generate the highest packing fraction when randomly poured into a container? In order to explore and navigate the enormous search space efficiently, we pair molecular dynamics simulations with artificial evolution. Arbitrary particle shape is represented by a set of overlapping spheres of varying diameter, enabling us to approximate smooth surfaces with a resolution proportional to the number of spheres included. We discover a family of planar triangular particles, whose packing fraction of ϕ ∼ 0.73 is among the highest experimental results for disordered packings of frictionless particles. We investigate how ϕ depends on the arrangement of spheres comprising an individual particle and on the smoothness of the surface. We validate the simulations with experiments using 3D-printed copies of the simplest member of the family, a planar particle consisting of three overlapping spheres with identical radius. Direct experimental comparison with 3D-printed aspherical ellipsoids demonstrates that the triangular particles pack exceedingly well not only in the limit of large system size but also when confined to small containers.

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.

Celebrating Soft Matter's 10th Anniversary: toward jamming by design

In materials science, high performance is typically associated with regularity and order, while disorder and the presence of defects are assumed to lead to sub-optimal outcomes. This holds for traditional solids such as crystals as well as for many types of nanoscale devices. However, there are circumstances where disorder can be harnessed to achieve performance not possible with approaches based on regularity. Recent research has shown opportunities specifically for soft matter. There, the phenomenon of jamming leads to unique emergent behavior that enables disordered, amorphous systems to switch reversibly between solid-like rigidity and fluid-like plasticity. This makes it possible to envision materials that can change stiffness or even shape adaptively. We review some of the progress in this direction, discussing examples where jamming has been explored from micro to macro scales in colloidal systems, suspensions, granular-materials-enabled soft robotics, and architecture. We focus in particular on how the jammed aggregate state can be tailored by controlling particle level properties and discuss very recent ideas that provide an important first step toward actual design of specifically targeted jamming behavior.

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.

Effect of size polydispersity versus particle shape in dense granular media

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.

Local fluctuations and spatial correlations in granular flows under constant-volume quasistatic shear

We investigate the local fluctuations in dense granular media subjected to athermal, quasistatic shearing, based on three-dimensional discrete element method simulations. By shearing granular assemblies of different polydispersities under constant-volume constraint, we quantify the characteristics of local structures (in terms of local volume and local anisotropy) and local deformation (using local shear strain and nonaffine displacement). The distribution of the local volume in a granular medium is found unchanged during the entire shearing process, which indicates a constant temperaturelike compactivity for the material. The compactivity is not, however, equilibrated among different particle groups in a polydisperse assembly. The local structures of a disordered granular assembly are inherently anisotropic. The fluctuations in local anisotropy can be well captured by a gamma or mixed-gamma distribution function, which is also unchanged during the shear. The local anisotropic orientation evolves towards the coaxial direction of the stress anisotropy with shear. The deformation characteristics of a jammed granular medium have their origins in the structural amorphousness. The local shear strain field depicts clear shear transformation zones which act as plasticity carriers. The spatial correlation of the local shear strains exhibits a fourfold pattern which is stronger in the stress deviatoric planes than in the stress isotropic plane. The fluctuations of nonaffine displacement suggest an isotropic granular temperature and an isotropic spatial correlation independent of the stress state. Both the local strain and the nonaffine displacement exhibit a power-law decayed distribution with a long-range correlation. We further modify the shear-transformation-zone theory to predict the pressure-dependent constitutive behavior of a sheared granular material and compare its prediction with our simulation data.

Packings of irregular polyhedral particles: strength, structure, and effects of angularity

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.

Cohesive granular materials composed of nonconvex particles

The macroscopic cohesion of granular materials made up of sticky particles depends on the particle shapes. We address this issue by performing contact dynamics simulations of 2D packings of nonconvex aggregates. We find that the macroscopic cohesion is strongly dependent on the strain and stress inhomogeneities developing inside the material. The largest cohesion is obtained for nearly homogeneous deformation at the beginning of unconfined axial compression and it evolves linearly with nonconvexity. Interestingly, the aggregates in a sheared packing tend to form more contacts with fewer neighboring aggregates as the degree of nonconvexity increases. We also find that shearing leads either to an isotropic distribution of tensile contacts or to the same privileged direction as that of compressive contacts.

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.

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.

Discrete simulation of dense flows of polyhedral grains down a rough inclined plane

The influence of grain angularity on the properties of dense flows down a rough inclined plane are investigated. Three-dimensional numerical simulations using the nonsmooth contact dynamics method are carried out with both spherical (rounded) and polyhedral (angular) grain assemblies. Both sphere and polyhedra assemblies abide by the flow start and stop laws, although much higher tilt angle values are required to trigger polyhedral grain flow. In the dense permanent flow regime, both systems show similarities in the bulk of the material (away from the top free surface and the substrate), such as uniform values of the solid fraction, inertial number and coordination number, or linear dependency of the solid fraction and effective friction coefficient with the inertial number. However, discrepancies are also observed between spherical and polyhedral particle flows. A dead (or nearly arrested) zone appears in polyhedral grain flows close to the rough bottom surface, reflected by locally concave velocity profiles, locally larger coordination number and solid fraction values, smaller inertial number values. This dead zone disappears for smooth bottom surfaces. In addition, unlike sphere assemblies, polyhedral grain assemblies exhibit significant normal stress differences, which increase close to the substrate.

Stability, deformation, and variability of granular fills composed of polyhedral particles

By means of extensive contact dynamics simulations, we investigate the mechanical equilibrium and deformation of a granular material composed of irregular polyhedral particles confined between two horizontal frictional planes. We show that, as a consequence of mobilized wall-particle friction forces at the top and bottom boundaries, the transient deformation induced by a constant vertical load is controlled by the aspect ratio (thickness over width) of the packing as well as the stress ratio. The transient deformation declines considerably for increasingly smaller aspect ratios and grows with the stress ratio. From the simulation data for a large number of independent configurations, we find that sample-to-sample fluctuations of the deformation have a broad distribution and they scale with the average deformation. We also analyze the evolution of particle connectivity during settlement and with the applied force. The face-face and edge-face contacts between polyhedral particles concentrate strong force chains with a growing proportion as a function of the applied force.

Strong interlocking of nonconvex particles in random packings

We present a numerical study of random packings made of nonconvex grains. These particles are built by the agglomeration of overlapping spheres in order to control their sphericity φ. The contact number C is found to be much larger than the coordination number Z, providing a significant difference with convex grains. The packing properties are found to be highly dependent on the morphological parameters of the grains : packing fractions as low as 0.3 have been reached. More importantly, the way nonconvex grains develop multiple contacts, i.e., interlocking, is found to be a relevant effect in such packings. Interlocking provides more stability to loose packings.

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.