Internal states of model isotropic granular packings. III. Elastic properties

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.

Compacting an assembly of soft balls far beyond the jammed state: Insights from three-dimensional imaging

Very soft grain assemblies have unique shape-changing capabilities that allow them to be compressed far beyond the rigid jammed state by filling void spaces more effectively. However, accurately following the formation of these systems by monitoring the creation of new contacts, monitoring the changes in grain shape, and measuring grain-scale stresses is challenging. We developed an experimental method that overcomes these challenges and connects their microscale behavior to their macroscopic response. By tracking the local strain energy during compression, we reveal a transition from granular-like to continuous-like material. Mean contact geometry is shown to vary linearly with the packing fraction, which is supported by a mean field approximation. We also validate a theoretical framework which describes the compaction from a local view. Our experimental framework provides insights into the granular micromechanisms and opens perspectives for rheological analysis of highly deformable grain assemblies in various fields ranging from biology to engineering.

Three-dimensional compaction of soft granular packings

This paper analyzes the compaction behavior of assemblies composed of soft (elastic) spherical particles beyond the jammed state, using three-dimensional non-smooth contact dynamic simulations. The assemblies of particles are characterized using the evolution of the packing fraction, the coordination number, and the von Misses stress distribution within the particles as the confining stress increases. The packing fraction increases and tends toward a maximum value close to 1, and the mean coordination number increases as a square root of the packing fraction. As the confining stress increases, a transition is observed from a granular-like material with exponential tails of the shear stress distributions to a continuous-like material characterized by Gaussian-like distributions of the shear stresses. We develop an equation that describes the evolution of the packing fraction as a function of the applied pressure. This equation, based on the micromechanical expression of the granular stress tensor, the limit of the Hertz contact law for small deformation, and the power-law relation between the packing fraction and the coordination of the particles, provides good predictions from the jamming point up to very high densities without the need for tuning any parameters.

Shear response of granular packings compressed above jamming onset

We investigate the mechanical response of jammed packings of repulsive, frictionless spherical particles undergoing isotropic compression. Prior simulations of the soft-particle model, where the repulsive interactions scale as a power law in the interparticle overlap with exponent α, have found that the ensemble-averaged shear modulus 〈G(P)〉 increases with pressure P as ∼P^{(α-3/2)/(α-1)} at large pressures. 〈G〉 has two key contributions: (1) continuous variations as a function of pressure along geometrical families, for which the interparticle contact network does not change, and (2) discontinuous jumps during compression that arise from changes in the contact network. Using numerical simulations, we show that the form of the shear modulus G^{f} for jammed packings within near-isostatic geometrical families is largely determined by the affine response G^{f}∼G_{a}^{f}, where G_{a}^{f}/G_{a0}=(P/P_{0})^{(α-2)/(α-1)}-P/P_{0}, P_{0}∼N^{-2(α-1)} is the characteristic pressure at which G_{a}^{f}=0, G_{a0} is a constant that sets the scale of the shear modulus, and N is the number of particles. For near-isostatic geometrical families that persist to large pressures, deviations from this form are caused by significant nonaffine particle motion. We further show that the ensemble-averaged shear modulus 〈G(P)〉 is not simply a sum of two power laws, but 〈G(P)〉∼(P/P_{c})^{a}, where a≈(α-2)/(α-1) in the P→0 limit and 〈G(P)〉∼(P/P_{c})^{b}, where b≳(α-3/2)/(α-1), above a characteristic pressure that scales as P_{c}∼N^{-2(α-1)}.

A micro-mechanical compaction model for granular mix of soft and rigid particles

We use bi-dimensional non-smooth contact dynamics simulations to analyze the isotropic compaction of mixtures composed of rigid and deformable incompressible particles. Deformable particles are modeled using the finite-element method and following a hyper-elastic neo-Hookean constitutive law. The evolution of the packing fraction, bulk modulus and particle connectivity, beyond the jamming point, are characterized as a function of the applied stresses for different proportion of rigid/soft particles and two values of friction coefficient. Based on the granular stress tensor, a micro-mechanical expression for the evolution of the packing fraction and the bulk modulus are proposed. This expression is based on the evolution of the particle connectivity together with the bulk behaviour of a single representative deformable particle. A constitutive compaction equation is then introduced, set by well-defined physical quantities, given a direct prediction of the maximum packing fraction φmax as a function of the proportion of rigid/soft particles.

Quasistatic response of loose cohesive granular materials

DEM-simulated model cohesive assemblies of spherical grains of diameter d, with contact tensile strength F0, once prepared in loose states, are quasistatically subjected to growing isotropic pressure P, and then to triaxial compression, maintaining lateral stresses σ2 = σ3 = P while increasing axial stress σ1 = P + q and strain є1. Reduced pressure P* = d2P/F0 varies from 0.1 (cohesion dominated case, for which systems typically equilibrate with solid fraction Ф ≃ 0.35), to large values for which the cohesionless behavior is retrieved. In triaxial compression, while the moderate strain response (є1 ~ 0.1) is influenced by initial coordination numbers and mesoscale heterogeneities, the approach to the critical state, as both q (deviator) and Ф steadily increase, gets slower for smaller P*. Critical ratio q/P strongly increases for decreasing P*, as roughly predicted in an “effective stress” scheme. Anomalously small elastic moduli are observed in the gel-like structures. While extensive geometric rearrangements take place, no shear banding is observed. Loose cohesive granular assemblies are thus capable of large quasistatic stable plastic strains and ductile rupture.

Rheology of soft granular materials: uniaxial compression

Soft granular materials are assemblies of highly deformable grains interacting via surface forces. The large grain deformations of these materials differ them from hard granular systems, in which, their behaviors are essentially governed by grain rearrangements. In this paper, we study the uniaxial compression of soft granular materials using a numerical approach based on the Material Point Method allowing for large grain deformations, coupled with the Contact Dynamics method for the treatment of unilateral frictional contacts between grains. Considering the neo-Hookean and elasto-plastic grains, the compaction of 2D soft granular packings is analyzed. We focus essentially on the evolution of the packing vertical stress as a function of the packing fraction and the predictive models are proposed.

Elastic wave velocity and attenuation in granular material

Discrete Elements Method simulations are carried out to investigate waves propagation in isotropic, frictional granular media. The focus is on the effects of confining pressure, microstructure and input frequency on both wave velocity and attenuation. The latter is described via the seismic quality factor Q and three different measurement approaches are compared, in time and frequency domain. The simulation data validate previous findings on the scaling of wave velocity with confining pressure and coordination number. The quality factor Q shows a non-monotonic behavior with input frequency.

Optimal packing in 2D and 3D granular systems: Density, connectivity, and force distributions

In this work, we explore the influence of the grain size distribution (GSD) on density, connectivity and internal forces distributions, for both 2D and 3D granular packings built mechanically. For power law GSDs, we show that there is an exponent for which density and connectivity are optimized, and this exponent is close to those that characterize other well known GSDs such as the Fuller and Thompson distribution and the Appollonian packing. In addition, we studied the distributions of normal forces, finding that these can be well described by a power-law tail, specially for the GSDs with large size span. These results highlight the role of the GSD on internal structure and suggest important consequences on macroscopic properties.

Compaction of mixtures of rigid and highly deformable particles: A micromechanical model

We analyze the isotropic compaction of mixtures composed of rigid and deformable incompressible particles by the nonsmooth contact dynamics approach. The deformable bodies are simulated using a hyperelastic neo-Hookean constitutive law by means of classical finite elements. We characterize the evolution of the packing fraction, the elastic modulus, and the connectivity as a function of the applied stresses when varying the interparticle coefficient of friction. We show first that the packing fraction increases and tends asymptotically to a maximum value ϕ_{max}, which depends on both the mixture ratio and the interparticle friction. The bulk modulus is also shown to increase with the packing fraction and to diverge as it approaches ϕ_{max}. From the micromechanical expression of the granular stress tensor, we develop a model to describe the compaction behavior as a function of the applied pressure, the Young modulus of the deformable particles, and the mixture ratio. A bulk equation is also derived from the compaction equation. This model lays on the characterization of a single deformable particle under compression together with a power-law relation between connectivity and packing fraction. This compaction model, set by well-defined physical quantities, results in outstanding predictions from the jamming point up to very high densities and allows us to give a direct prediction of ϕ_{max} as a function of both the mixture ratio and the friction coefficient.

Compaction Model for Highly Deformable Particle Assemblies

The compaction behavior of deformable grain assemblies beyond jamming remains bewildering, and existing models that seek to find the relationship between the confining pressure P and solid fraction ϕ end up settling for empirical strategies or fitting parameters. Using a coupled discrete-finite element method, we analyze assemblies of highly deformable frictional grains under compression. We show that the solid fraction evolves nonlinearly from the jamming point and asymptotically tends to unity. Based on the micromechanical definition of the granular stress tensor, we develop a theoretical model, free from ad hoc parameters, correctly mapping the evolution of ϕ with P. Our approach unveils the fundamental features of the compaction process arising from the joint evolution of grain connectivity and the behavior of single representative grains. This theoretical framework also allows us to deduce a bulk modulus equation showing an excellent agreement with our numerical data.

Oscillatory instabilities in three-dimensional frictional granular matter

The dynamics of amorphous granular matter with frictional interactions cannot be derived in general from a Hamiltonian and therefore displays oscillatory instabilities stemming from the onset of complex eigenvalues in the stability matrix. These instabilities were discovered in the context of one- and two-dimensional systems, while the three-dimensional case was never studied in detail. Here we fill this gap by deriving and demonstrating the presence of oscillatory instabilities in a three-dimensional granular packing. We study binary assemblies of spheres of two sizes interacting via classical Hertz and Mindlin force laws for the longitudinal and tangent interactions, respectively. We formulate analytically the stability matrix in three dimensions and observe that a couple of complex eigenvalues emerge at the onset of the instability as in the case of frictional disks in two dimensions. The dynamics then shows oscillatory exponential growth in the mean-square displacement, followed by a catastrophic event in which macroscopic portions of mechanical stress and energy are lost. The generality of these results for any choice of forces that break the symplectic Hamiltonian symmetry is discussed.

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.

Shear strength of wet granular materials: Macroscopic cohesion and effective stress : Discrete numerical simulations, confronted to experimental measurements

Rheometric measurements on assemblies of wet polystyrene beads, in steady uniform quasistatic shear flow, for varying liquid content within the small saturation (pendular) range of isolated liquid bridges, are supplemented with a systematic study by discrete numerical simulations. The numerical results agree quantitatively with the experimental ones provided that the intergranular friction coefficient is set to the value [Formula: see text], identified from the behaviour of the dry material. Shear resistance and solid fraction [Formula: see text] are recorded as functions of the reduced pressure [Formula: see text], which, defined as [Formula: see text], compares stress [Formula: see text], applied in the velocity gradient direction, to the tensile strength [Formula: see text] of the capillary bridges between grains of diameter a, and characterizes cohesion effects. The simplest Mohr-Coulomb relation with [Formula: see text]-independent cohesion c applies as a good approximation for large enough [Formula: see text] (typically [Formula: see text]. Numerical simulations extend to different values of μ and, compared to experiments, to a wider range of [Formula: see text]. The assumption that capillary stresses act similarly to externally applied ones onto the dry granular contact network (effective stresses) leads to very good (although not exact) predictions of the shear strength, throughout the numerically investigated range [Formula: see text] and [Formula: see text]. Thus, the internal friction coefficient [Formula: see text] of the dry material still relates the contact force contribution to stresses, [Formula: see text], while the capillary force contribution to stresses, [Formula: see text], defines a generalized Mohr-Coulomb cohesion c, depending on [Formula: see text] in general. c relates to [Formula: see text] , coordination numbers and capillary force network anisotropy. c increases with liquid content through the pendular regime interval, to a larger extent, the smaller the friction coefficient. The simple approximation ignoring capillary shear stress [Formula: see text] (referred to as the Rumpf formula) leads to correct approximations for the larger saturation range within the pendular regime, but fails to capture the decrease of cohesion for smaller liquid contents.

Dynamic induced softening in frictional granular materials investigated by discrete-element-method simulation

A granular system composed of frictional glass beads is simulated using the discrete element method. The intergrain forces are based on the Hertz contact law in the normal direction with frictional tangential force. The damping due to collision is also accounted for. Systems are loaded at various stresses and their quasistatic elastic moduli are characterized. Each system is subjected to an extensive dynamic testing protocol by measuring the resonant response to a broad range of ac drive amplitudes and frequencies via a set of diagnostic strains. The system, linear at small ac drive amplitudes, has resonance frequencies that shift downward (i.e., modulus softening) with increased ac drive amplitude. Detailed testing shows that the slipping contact ratio does not contribute significantly to this dynamic modulus softening, but the coordination number is strongly correlated to this reduction. This suggests that the softening arises from the extended structural change via break and remake of contacts during the rearrangement of bead positions driven by the ac amplitude.

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.

Numerical study of one-dimensional compression of granular materials. I. Stress-strain behavior, microstructure, and irreversibility

The behavior of a model granular material, made of slightly polydisperse beads with Hertz-Mindlin elastic-frictional contacts, in oedometric compression (i.e., compression along one axis, with no lateral strain) is studied by grain-level numerical simulations. We systematically investigate the influence of the (idealized) packing process on the microstructure and stresses in the initial, weakly confined equilibrium state, and prepare both isotropic and anisotropic configurations differing in solid fraction Φ and coordination number z. Φ (ranging from maximally dense to moderately loose), z (which might vary independently of Φ in dense systems), fabric and force anisotropy parameters, and the ratio K_{0} of lateral stresses σ_{2}=σ_{3} to stress σ_{1} in the compression direction are monitored in oedometric compression in which σ_{1} varies by more than three orders of magnitude. K_{0} reflects the anisotropy of the assembling process and may remain nearly constant in further loading if the material is already oedometrically compressed (as a granular gas) in the preparation stage. Otherwise, it tends to decrease steadily over the investigated stress range. It is related to force and fabric anisotropy parameters by a simple formula. Elastic moduli, separately computed with an appropriate matrix method, may express the response to very small stress increments about the transversely isotropic well-equilibrated states along the loading path, although oedometric compression proves an essentially anelastic process, mainly due to friction mobilization, with large irreversible effects apparent upon unloading. While the evolution of axial strain ε_{1} and solid fraction Φ (or of the void ratio e=-1+1/Φ) with axial stress σ_{1} is very nearly reversible, especially in dense samples, z is observed to decrease (as previously observed in isotropic compression) after a compression cycle if its initial value was high. K_{0} relates to the evolution of internal variables and may exceed 1 in unloading. The considerably greater irreversibility of oedometric compression reported in sands, compared to our model systems, should signal contact plasticity or damage.

Numerical study of one-dimensional compression of granular materials. II. Elastic moduli, stresses, and microstructure

The elastic moduli of a transversely isotropic model granular material, made of slightly polydisperse elastic-frictional spherical beads, in equilibrium along a one-dimensional (oedometric) compression path, as described in the companion paper [M. H. Khalili et al., Phys. Rev. E 95, 032907 (2017)]10.1103/PhysRevE.95.032907, are investigated by numerical simulations. The relations of the five independent moduli to stresses, density, coordination number, fabric and force anisotropies are studied for different internal material states along the oedometric loading path. It is observed that elastic moduli, as in isotropic packs, are primarily determined by the coordination number, with anomalously small shear moduli in poorly coordinated systems, whatever their density. Such states also exhibit faster increasing moduli in compression, and larger off-diagonal moduli and Poisson ratios. Anisotropy affects the longitudinal moduli C_{11} in the axial direction and C_{22} in the transverse directions, and the shear modulus in the transverse plane C_{44}, more than the shear modulus in a plane containing the axial direction C_{55}. The results are compared to available experiments on anisotropic bead packs, revealing, despite likely differences in internal states, a very similar range of stiffness level (linked to coordination), and semiquantitative agreement as regards the influence of anisotropy. Effective medium theory (the Voigt approach) provides quite inaccurate predictions of the moduli. It also significantly underestimates ratios C_{11}/C_{22} (varying between 1 and 2.2) and C_{55}/C_{44} (varying from 1 to 1.6), which characterize elastic anisotropy, except in relatively weakly anisotropic states. The bulk modulus for isotropic compression and the compliance corresponding to stress increments proportional to the previous stress values are the only elastic coefficients to be correctly estimated by available predictive relations. We discuss the influences of fabric and force anisotropies onto elastic anisotropy, showing in particular that the former dominates in sample series that are directly assembled in anisotropic configurations and keep a roughly constant lateral to axial stress ratio under compression.

Erosion and deposition in depth-averaged models of dense, dry, inclined, granular flows

We derive expressions for the rates of erosion and deposition at the interface between a dense, dry, inclined granular flow and an erodible bed. In obtaining these, we assume that the interface between the flowing grains and the bed moves with the speed of a pressure wave in the flow, for deposition, or with the speed of a disturbance through the contacting particles in the bed, for erosion. We employ the expressions for the rates of erosion and deposition to show that after an abrupt change in the angle of inclination of the bed the characteristic time for the motion of the interface is much shorter than the characteristic time of the flow. This eliminates the need for introducing models of erosion and deposition rate in the mass balance; and the instantaneous value of the particle flux is the same function of the instantaneous value of the flow depth as in a steady, uniform flow.

Strength of anisotropy in a granular material: Linear versus nonlinear contact model

In this paper, we deal with anisotropy in an idealized granular material made of a collection of frictional, elastic, contacting particles. We present a theoretical analysis for an aggregate of particles isotropically compressed and then sheared, in which two possible contacts laws between particles are considered: a linear contact law, where the contact stiffness is constant; and a nonlinear contact law, where the contact stiffness depends on the overlapping between particles. In the former case the anisotropy observed in the aggregate is associated with particle arrangement. In fact, although the aggregate is initially characterized by an isotropic network of contacts, during the loading, an anisotropic texture develops, which is measured by a fabric tensor. With a nonlinear contact law it is possible to develop anisotropy because contacting stiffnesses are different, depending on the orientation of the contact vectors with respect to the axis of the applied deformation. We find that before the peak load is reached, an aggregate made of particles with a linear contact law develops a much smaller anisotropy compared with that of an aggregate with a nonlinear law.

Beyond linear elasticity: jammed solids at finite shear strain and rate

The shear response of soft solids can be modeled with linear elasticity, provided the forcing is slow and weak. Both of these approximations must break down when the material loses rigidity, such as in foams and emulsions at their (un)jamming point - suggesting that the window of linear elastic response near jamming is exceedingly narrow. Yet precisely when and how this breakdown occurs remains unclear. To answer these questions, we perform computer simulations of stress relaxation and shear start-up tests in athermal soft sphere packings, the canonical model for jamming. By systematically varying the strain amplitude, strain rate, distance to jamming, and system size, we identify characteristic strain and time scales that quantify how and when the window of linear elasticity closes, and relate these scales to changes in the microscopic contact network.

Softening of stressed granular packings with resonant sound waves

We perform numerical simulations of a two-dimensional bidisperse granular packing subjected to both a static confining pressure and a sinusoidal dynamic forcing applied by a wall on one edge of the packing. We measure the response experienced by a wall on the opposite edge of the packing and obtain the resonant frequency of the packing as the static or dynamic pressures are varied. Under increasing static pressure, the resonant frequency increases, indicating a velocity increase of elastic waves propagating through the packing. In contrast, when the dynamic amplitude is increased for fixed static pressure, the resonant frequency decreases, indicating a decrease in the wave velocity. This occurs both for compressional and for shear dynamic forcing and is in agreement with experimental results. We find that the average contact number Zc at the resonant frequency decreases with increasing dynamic amplitude, indicating that the elastic softening of the packing is associated with a reduced number of grain-grain contacts through which the elastic waves can travel. We image the excitations created in the packing and show that there are localized disturbances or soft spots that become more prevalent with increasing dynamic amplitude. Our results are in agreement with experiments on glass bead packings and earth materials such as sandstone and granite and may be relevant to the decrease in elastic wave velocities that has been observed to occur near fault zones after strong earthquakes, in surficial sediments during strong ground motion, and in structures during earthquake excitation.

Ultrasonic evaluation of the morphological characteristics of metallic powders in the context of mechanical alloying

An ultrasonic method is proposed to characterize the morphological (geometrical) aspects of powders through the elastic modulus dependence of their packing on the factors of polydispersity, coordination number and particle shape. During the mechanical alloying process, the variation in geometrical characteristics of powders provides critical information. Ultrasonic parameters are shown to be sensitive not only to the average contact number per bead (i.e. the coordination number) but also to characteristics of the bead size distribution, when given the same sample preparation and confining pressure. These parameters, in turn, are sensitive to both the granular medium polydispersity and particle shapes. A non-monotonic behavior of the ultrasonic velocity (and of the derived compressional wave modulus) is observed throughout the alloying process, which thus offers possibilities for powder structure monitoring.

On granular elasticity

Mesoscopic structures form in dense granular materials due to the self-organisation of the constituent particles. These structures have internal structural degrees of freedom in addition to the translational degree of freedom. The resultant granular elasticity, which exhibits intrinsic variations and inevitable relaxation, is a key quantity that accounts for macroscopic solid- or fluid-like properties and the transitions between them. In this work, we propose a potential energy landscape (PEL) with local stable basins and low elastic energy barriers to analyse the nature of granular elasticity. A function for the elastic energy density is proposed for stable states and is further calibrated with ultrasonic measurements. Fluctuations in the elastic energy due to the evolution of internal structures are proposed to describe a so-called configuration temperature T(c) as a counterpart of the classical kinetic granular temperature T(k) that is attributed to the translational degrees of freedom. The two granular temperatures are chosen as the state variables, and a fundamental equation is established to develop non-equilibrium thermodynamics for granular materials. Due to the relatively low elastic energy barrier in the PEL, granular elasticity relaxes more under common mechanical loadings, and a simple model based on mean-field theory is developed to account for this behaviour.

Three-dimensional simulations of nanopowder compaction processes by granular dynamics method

In order to describe and to study the processes of cold compaction within the discrete element method a three-dimensional model of nanosized powder is developed. The elastic forces of repulsion, the tangential forces of "friction" (Cattaneo-Mindlin), and the dispersion forces of attraction (van der Waals-Hamaker), as well as the formation and destruction of hard bonds between the individual particles are taken into account. The monosized powders with the size of particles in the range 10-40 nm are simulated. The simulation results are compared to the experimental data of the alumina nanopowders compaction. It is shown that the model allows us to reproduce experimental data reliably and, in particular, describes the size effect in the compaction processes. A number of different external loading conditions is used in order to perform the theoretical and experimental researches. The uniaxial compaction (the closed-die compaction), the biaxial (radial) compaction, and the isotropic compaction (the cold isostatic pressing) are studied. The real and computed results are in a good agreement with each other. They reveal a weak sensitivity of the oxide nanopowders to the loading condition (compaction geometry). The application of the continuum theory of the plastically hardening porous body, which is usually used for the description of powders, is discussed.

Contact anisotropy and coordination number for a granular assembly: a comparison of distinct-element-method simulations and theory

We study an ideal granular aggregate consisting of elastic spherical particles, isotropic in stress and anisotropic in the contact network. Because of the contact anisotropy, a confining pressure applied at zero deviatoric stress, produces shear strain as well as volume strain. Our goal is to predict the coordination number k, the average number of contacts per particle, and the magnitude of the contact anisotropy ɛ, from knowledge of the elastic moduli of the aggregate. We do this through a theoretical model based upon the well known effective medium theory. However, rather than focusing on the moduli, we consider their ratios over the moduli of an equivalent isotropic state. We observe good agreement between numerical simulation and theory.

Protocol dependence of mechanical properties in granular systems

We study the protocol dependence of the mechanical properties of granular media by means of computer simulations. We control a protocol of realizing disk packings in a systematic manner. In 2D, by keeping material properties of the constituents identical, we carry out compaction with various strain rates. The disk packings exhibit the strain rate dependence of the critical packing fraction above which the pressure becomes non-zero. The observed behavior contrasts with the well-studied jamming transitions for frictionless disk packings. We also observe that the elastic moduli of the disk packings depend on the strain rate logarithmically. Our results suggest that there exists a time-dependent state variable to describe macroscopic material properties of disk packings, which depend on its protocol.

Rheology of granular materials composed of nonconvex particles

By means of contact dynamics simulations, we investigate the shear strength and internal structure of granular materials composed of two-dimensional nonconvex aggregates. We find that the packing fraction first grows as the nonconvexity is increased but declines at higher nonconvexity. This unmonotonic dependence reflects the competing effects of pore size reduction between convex borders of aggregates and gain in porosity at the nonconvex borders that are captured in a simple model fitting nicely the simulation data both in the isotropic and sheared packings. On the other hand, the internal angle of friction increases linearly with nonconvexity and saturates to a value independent of nonconvexity. We show that fabric anisotropy, force anisotropy, and friction mobilization, all enhanced by multiple contacts between aggregates, govern the observed increase of shear strength and its saturation with increasing nonconvexity. The main effect of interlocking is to dislocate frictional dissipation from the locked double and triple contacts between aggregates to the simple contacts between clusters of aggregates. This self-organization of particle motions allows the packing to keep a constant shear strength at high nonconvexity.

Elasticity of floppy and stiff random networks

We study the linear and nonlinear elastic behavior of amorphous systems using a two-dimensional random network of harmonic springs as a model system. A natural characterization of these systems arises in terms of the network coordination (average number of springs per node) relative to that of a marginally rigid network deltaz: a floppy network has deltaz<0, while a stiff network has deltaz>0. Under the influence of an externally applied load, we observe that the response of both floppy and stiff networks is controlled by the critical point corresponding to the onset of rigidity. We use numerical simulations to compute the exponents which characterize the shear modulus, the heterogeneity of the response, and the network stiffening as a function of deltaz and derive these theoretically, thus allowing us to predict aspects of the mechanical response of glasses and fibrous networks.

Solidlike behavior and anisotropy in rigid frictionless bead assemblies

We investigate the structure and mechanical behavior of assemblies of frictionless, nearly rigid equal-sized beads, in the quasistatic limit, by numerical simulation. Three different loading paths are explored: triaxial compression, triaxial extension, and simple shear. Generalizing a recent result, we show that the material, despite rather strong finite sample size effects, is able to sustain a finite deviator stress in the macroscopic limit, along all three paths, without dilatancy. The shape of the yield surface in principal stress space differs somewhat from the Mohr-Coulomb prediction, and is more adequately described by the Lade-Duncan or Matsuoka-Nakai criteria. We study geometric characteristics and force networks under varying stress levels within the supported range. Although the scalar state variables stay equal to the values observed in systems under isotropic pressure, the material, once subjected to a deviator stress, possesses some fabric and force distribution anisotropies. Each kind of anisotropy can be described, in good approximation, by a single parameter. Within the supported stress range, along each one of the three investigated stress paths, among those three quantities: deviator stress to mean stress ratio, fabric anisotropy parameter, force anisotropy parameter, any one determines the values of the two others. The pair correlation function also exhibits short range anisotropy, up to a distance between bead surfaces of the order of 10% of the diameter. The tensor of elastic moduli is shown to possess a nearly singular, uniaxial structure related to stress anisotropy. Possible stress-strain relations in monotonic loading paths are also discussed.

Computer simulation of model cohesive powders: plastic consolidation, structural changes, and elasticity under isotropic loads

The quasistatic behavior of a simple two-dimensional model of a cohesive powder under isotropic loads is investigated by discrete element simulations. We ignore contact plasticity and focus on the effect of geometry and collective rearrangements on the material behavior. The loose packing states, as assembled and characterized in a previous numerical study [Gilabert, Roux, and Castellanos, Phys. Rev. E 75, 011303 (2007)], are observed, under growing confining pressure P , to undergo important structural changes, while solid fraction Phi irreversibly increases (typically, from 0.4-0.5 to 0.75-0.8). The system state goes through three stages, with different forms of the plastic consolidation curve, i.e., Phi as a function of the growing reduced pressure P;{*}=PaF_{0} , defined with adhesion force F0 and grain diameter a . In the low-confinement regime (I), the system undergoes negligible plastic compaction, and its structure is influenced by the assembling process. In regime II the material state is independent of initial conditions, and the void ratio varies linearly with lnP [i.e., Delta(1Phi)=lambdaDelta(lnP;{*}) ], as described in the engineering literature. Plasticity index lambda is reduced in the presence of a small rolling resistance (RR). In the last stage of compaction (III), Phi approaches an asymptotic, maximum solid fraction Phi_{max} , as a power law Phi_{max}-Phi proportional, variant(P;{*});{-alpha} , with alpha approximately 1 , and properties of cohesionless granular packs are gradually retrieved. Under consolidation, while the range xi of fractal density correlations decreases, force patterns reorganize from self-balanced clusters to force chains, with correlative evolutions of force distributions, and elastic moduli increase by a large amount. Plastic deformation events correspond to very small changes in the network topology, while the denser regions tend to move like rigid bodies. Elastic properties are dominated by the bending of thin junctions in loose systems. For growing RR those tend to form particle chains, the folding of which, rather than tensile ruptures, controls plastic compaction.

Strong-weak network anisotropy switching and hysteresis in three-dimensional granular materials

We address hysteresis of three-dimensional polydisperse granular packs, comparing macro- and microscopic viewpoints, to reveal their elastic/inelastic mechanics and force network anisotropy. During the uniaxial loading-unloading cycle of an appropriately prepared pack, one can decompose the force network into weak and strong subnetworks. The first stages of loading exhibit arching, where all the fabric displays negative anisotropy. For later stages, the strong (weak) network shows positive (negative) anisotropy. On unloading, the force network progresses to a fabric wide hydrostatic point, where the anisotropies of the weak and strong subnetworks switch signs. During the loading stage, a Mohr circle analysis permits the identification of a well-defined macroscopic internal friction angle, whose value is larger than that of grain-grain interactions. To analyze unloading, a generalized local Coulomb-friction argument predicts a continuously changing friction angle, that vanishes at the hydrostatic point. A suggestive interplay between microscopic friction and fabric structure, at different loding stages, is proposed as the mechanism for the emergence of a macro internal friction angle.

Sound propagation in a constrained lattice of beads: high-frequency behavior and dispersion relation

We report on acoustic wave propagation in a regular array of nominally identical beads under isotropic static stress. The weak polydispersity of the beads makes the contact lattice random. Time-frequency analysis of the acoustic signal is performed and allows measurement of the full lattice dispersion relation. Comparison with the theoretical prediction for a perfect triangular lattice gives an indication of the level of randomness in the contact lattice. The results extend, in a consistent way, a previous study restricted to long wavelength propagation [B. Gilles and C. Coste, Phys. Rev. Lett. 90, 174302 (2003)]: The contact lattice is ordered by increasing the stress, and the smaller the wavelength, the higher the stress required to get regular lattice behavior. Measurements involving ballistic propagation of the coherent wave, whatever its frequency, evidence reversible lattice behavior under compression and/or decompression. Nevertheless, correlations of short wavelength incoherent waves are a sensitive probe of disorder, and allow us to exhibit a small irreversible evolution of the lattice.

Internal states of model isotropic granular packings. II. Compression and pressure cycles

This is the second paper of a series of three investigating, by numerical means, the geometric and mechanical properties of spherical bead packings under isotropic stresses. We study the effects of varying the applied pressure P (from 1 or 10 kPa up to 100 MPa in the case of glass beads) on several types of configurations assembled by different procedures, as reported in the preceding paper [I. Agnolin and J.-N. Roux, Phys. Rev. E 76, 061302 (2007)]. As functions of P , we monitor changes in solid fraction Phi, coordination number z, proportion of rattlers (grains carrying no force) x_(0) , the distribution of normal forces, the level of friction mobilization, and the distribution of near neighbor distances. Assuming that the contact law does not involve material plasticity or damage, Phi is found to vary very nearly reversibly with P in an isotropic compression cycle, but all other quantities, due to the frictional hysteresis of contact forces, change irreversibly. In particular, initial low P states with high coordination numbers lose many contacts in a compression cycle and end up with values of z and x_(0) close to those of the most poorly coordinated initial configurations. Proportional load variations which do not entail notable configuration changes can therefore nevertheless significantly affect contact networks of granular packings in quasistatic conditions.

Internal states of model isotropic granular packings. I. Assembling process, geometry, and contact networks

This is the first paper of a series of three, in which we report on numerical simulation studies of geometric and mechanical properties of static assemblies of spherical beads under an isotropic pressure. The influence of various assembling processes on packing microstructures is investigated. It is accurately checked that frictionless systems assemble in the unique random close packing (RCP) state in the low pressure limit if the compression process is fast enough, higher solid fractions corresponding to more ordered configurations with traces of crystallization. Specific properties directly related to isostaticity of the force-carrying structure in the rigid limit are discussed. With frictional grains, different preparation procedures result in quite different inner structures that cannot be classified by the sole density. If partly or completely lubricated they will assemble like frictionless ones, approaching the RCP solid fraction Phi_{RCP} approximately 0.639 with a high coordination number: z* approximately =6 on the force-carrying backbone. If compressed with a realistic coefficient of friction mu=0.3 packings stabilize in a loose state with Phi approximately 0.593 and z* approximately =4.5 . And, more surprisingly, an idealized "vibration" procedure, which maintains an agitated, collisional regime up to high densities results in equally small values of z* while Phi is close to the maximum value Phi_{RCP}. Low coordination packings have a large proportion (>10%) of rattlers--grains carrying no force--the effect of which should be accounted for on studying position correlations, and also contain a small proportion of localized "floppy modes" associated with divalent grains. Low-pressure states of frictional packings retain a finite level of force indeterminacy even when assembled with the slowest compression rates simulated, except in the case when the friction coefficient tends to infinity. Different microstructures are characterized in terms of near neighbor correlations on various scales, and some comparisons with available laboratory data are reported, although values of contact coordination numbers apparently remain experimentally inaccessible.