Rheology of three-dimensional packings of aggregates: microstructure and effects of nonconvexity

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.

Force transmission during repose of flexible granular chains

We study the mechanics of standing columns formed during the repose of flexible granular chains. It is one of the many intriguing behaviours exhibited by granular materials when links capable of transmitting tension exist between particles. We develop and calibrate a discrete element method contact model to simulate the mechanics of the macroscopic flexible granular chains and conduct simulations of the angle of repose experiments of these chains by extracting a chain-filled cylinder and allowing the material to flow out under gravity and repose. We evaluate various micro-mechanical, topological and macroscopic parameters to elucidate the mechanics of the repose behaviour of chain ensembles. It is the ability of the links connecting the individual particles to transmit tensile forces along the chain backbone that provides lateral stability to the column, enabling them to stand. In particular, the contact force rearrangement inside the columns generates a self-confining radial stress near the base of the columns, which provides an important stabilizing stress.

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.

Structural analysis of disordered dimer packings

Jammed disordered packings of non-spherical particles show significant variation in the packing density as a function of particle shape for a given packing protocol. Rotationally symmetric elongated shapes such as ellipsoids, spherocylinders, and dimers, e.g., pack significantly denser than spheres over a narrow range of aspect ratios, exhibiting a characteristic peak at aspect ratios of α_max ≈ 1.4-1.5. However, the structural features that underlie this non-monotonic behaviour in the packing density are unknown. Here, we study disordered packings of frictionless dimers in three dimensions generated by a gravitational pouring protocol in LAMMPS. Focusing on the characteristics of contacts as well as orientational and translational order metrics, we identify a number of structural features that accompany the formation of maximally dense packings as the dimer aspect ratio α is varied from the spherical limit. Our results highlight that dimer packings undergo significant structural changes as α increases up to α_max manifest in the reorganisation of the contact configurations between neighbouring dimers, increasing nematic order, and decreasing local translational order. Remarkably, for α > α_max our metrics remain largely unchanged, indicating that the peak in the packing density is related to the interplay of structural rearrangements for α < α_max and subsequent excluded volume effects with unchanged structure for α > α_max.

Shape or friction? Which of these characteristics drives the shear strength in granular systems?

The shape of the particles and local friction, separately, are known to strongly affect the macroscopic properties of an assembly of grains. But the combined effects of these two parameters still remain poorly described. By means of extensive two dimensional contact dynamics simulations, we perform a systematic analysis of the interplay between friction and shape on strength properties of granular systems. The shape of the particles is varied from disks to triangles, while the friction is varied from 0 to 0.7. We find that the macroscopic friction first increases with angularity, but it may decline (for low friction values), saturate (for intermediates friction values), or continue to increase (for large friction values) for the most angular shapes. In other words, the effect of the particle’s angularity on the shear strength depends on the level of sliding friction. In contrast, the effect of local friction on the shear strength does not depend on the specific properties of shape. The results presented here highlight the subtle coupling existing between shape and friction effects.

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.

Yielding, rigidity, and tensile stress in sheared columns of hexapod granules

Granular packings of nonconvex or elongated particles can form freestanding structures like walls or arches. For some particle shapes, such as staples, the rigidity arises from interlocking of pairs of particles, but the origins of rigidity for noninterlocking particles remains unclear. We report on experiments and numerical simulations of sheared columns of "hexapods," particles consisting of three mutually orthogonal sphero-cylinders whose centers coincide. We vary the length-to-diameter aspect ratio, α, of the sphero-cylinders and subject the packings to quasistatic direct shear. For small α, we observe a finite yield stress. For large α, however, the column becomes rigid when sheared, supporting stresses that increase sharply with increasing strain. Analysis of x-ray microcomputed tomography (micro-CT) data collected during the shear reveals that the stiffening is associated with a tilted, oblate cluster of hexapods near the nominal shear plane in which particle deformation and average contact number both increase. Simulation results show that the particles are collectively under tension along one direction, even though they do not interlock pairwise. These tensions comes from contact forces carrying large torques, and they are perpendicular to the compressive stresses in the packing. They counteract the tendency to dilate, thus stabilizing the particle cluster.

Universality of internal structure characteristics in granular media under shear

We examine the signatures of internal structure emerged from quasistatic shear responses of granular materials based on three-dimensional discrete element simulations. Granular assemblies consisting of spheres or nonspherical particles of different polydispersity are sheared from different initial densities and under different loading conditions (drained or undrained) steadily to reach the critical state (a state featured by constant stress and constant volume). The radial distribution function used to measure the packing structure is found to remain almost unchanged during the shearing process, regardless of the initial states or loading conditions of an assembly. Its specific form, however, varies with polydispersities in both grain size and grain shape. Set Voronoi tessellation is employed to examine the characteristics of local volume and anisotropy, and deformation. The local inverse solid fraction and anisotropy index are found following inverse Weibull and log-normal distributions, respectively, which are unique at the critical states. With further normalization, an invariant distribution for local volume and anisotropy is observed whose function can be determined by the polydispersities in both particle size and grain shape but bears no relevance to initial densities or loading conditions (or paths). An invariant Gaussian distribution is found for the local deformation for spherical packings, but no invariant distribution can be found for nonspherical packings where the distribution of normalized local volumetric strain is dependent on initial states.

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.

Active particles with desired orientation flowing through a bottleneck

We report extensive numerical simulations of the flow of anisotropic self-propelled particles through a constriction. In particular, we explore the role of the particles’ desired orientation with respect to the moving direction on the system flowability. We observe that when particles propel along the direction of their long axis (longitudinal orientation) the flow-rate notably reduces compared with the case of propulsion along the short axis (transversal orientation). And this is so even when the effective section (measured as the number of particles that are necessary to span the whole outlet) is larger for the case of longitudinal propulsion. This counterintuitive result is explained in terms of the formation of clogging structures at the outlet, which are revealed to have higher stability when the particles align along the long axis. This generic result might be applied to many different systems flowing through bottlenecks such as microbial populations or different kind of cells. Indeed, it has already a straightforward connection with recent results of pedestrian (which self-propel transversally oriented) and mice or sheep (which self-propel longitudinally oriented).

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.

Compaction of noncohesive and cohesive granular materials under vibrations: Experiments and stochastic model

We study the time evolution of the compaction of a noncohesive or cohesive granular material submitted to shaking through experiments and a stochastic model. Beyond well-known empirical expressions, we show that the characteristic time scales depend on the number of objects in the assembly. For a noncohesive granular material, the compaction time scale is governed by the number of individual grains in the system. In the case of a cohesive granular material, a two-scale model (individual particles and clusters) allows one to mimic the time evolution of the compaction of an actual cohesive powder driven by horizontal vibrations. In this case, the two time scales are associated with the numbers of clusters and grains, respectively.

Simulating competitive egress of noncircular pedestrians

We present a numerical framework to simulate pedestrian dynamics in highly competitive conditions by means of a force-based model implemented with spherocylindrical particles instead of the traditional, symmetric disks. This modification of the individuals' shape allows one to naturally reproduce recent experimental findings of room evacuations through narrow doors in situations where the contact pressure among the pedestrians was rather large. In particular, we obtain a power-law tail distribution of the time lapses between the passage of consecutive individuals. In addition, we show that this improvement leads to new features where the particles' rotation acquires great significance.

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.

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.

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.

Evolution of shear zones in granular materials

The evolution of wide shear zones or shear bands was investigated experimentally and numerically for quasistatic dry granular flows in split bottom shear cells. We compare the behavior of materials consisting of beads, irregular grains, such as sand, and elongated particles. Shearing an initially random sample, the zone width was found to significantly decrease in the first stage of the process. The characteristic shear strain associated with this decrease is about unity and it is systematically increasing with shape anisotropy, i.e., when the grain shape changes from spherical to irregular (e.g., sand) and becomes elongated (pegs). The strongly decreasing tendency of the zone width is followed by a slight increase which is more pronounced for rodlike particles than for grains with smaller shape anisotropy (beads or irregular particles). The evolution of the zone width is connected to shear-induced packing density change and for nonspherical particles it also involves grain reorientation effects. The final zone width is significantly smaller for irregular grains than for spherical beads.

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.

Particle shape effects on the stress response of granular packings

We present measurements of the stress response of packings formed from a wide range of particle shapes. Besides spheres these include convex shapes such as the Platonic solids, truncated tetrahedra, and triangular bipyramids, as well as more complex, non-convex geometries such as hexapods with various arm lengths, dolos, and tetrahedral frames. All particles were 3D-printed in hard resin. Well-defined initial packing states were established through preconditioning by cyclic loading under given confinement pressure. Starting from such initial states, stress-strain relationships for axial compression were obtained at four different confining pressures for each particle type. While confining pressure has the largest overall effect on the mechanical response, we find that particle shape controls the details of the stress-strain curves and can be used to tune packing stiffness and yielding. By correlating the experimentally measured values for the effective Young's modulus under compression, yield stress and energy loss during cyclic loading, we identify trends among the various shapes that allow for designing a packing's aggregate behavior.