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

Micromechanical Origin of Plasticity and Hysteresis in Nestlike Packings

Disordered packings of unbonded, semiflexible fibers represent a class of materials spanning contexts and scales. From twig-based bird nests to unwoven textiles, bulk mechanics of disparate systems emerge from the bending of constituent slender elements about impermanent contacts. In experimental and computational packings of wooden sticks, we identify prominent features of their response to cyclic oedometric compression: nonlinear stiffness, transient plasticity, and eventually repeatable velocity-independent hysteresis. We trace these features to their micromechanic origins, identified in characteristic appearance, disappearance, and displacement of internal contacts.

Investigation on the Influencing Factors of K0 of Granular Materials Using Discrete Element Modelling

Earth pressure coefficient at rest K0 is commonly estimated by empirical equations, which to date has had insufficient accuracy and universality. For better prediction, the investigation on the factors influencing K0 is required. A series of discrete element method (DEM) simulations of oedometer tests are conducted to verify the key factors influencing K0 of granular materials. The influences of initial fabric anisotropy, particle shape, initial void ratio, inter-particle friction angle is investigated. The evolution of microstructure is monitored during the tests to reveal the relationship between the microstructure evolution and K0 values. The results show that the effect of fabric anisotropy exists but is limited. Particle shape, initial void ratio, and inter-particle friction angle all significantly affect the K0 values alone. According to the DEM results, an attempt is made to propose a more reasonable empirical equation in which K0 is a function of relative density, critical state friction angle, and “shape factor”. This new empirical equation has higher accuracy and can consider the effect of particle shape, inspiring the determination of K0 values in practical engineering.

Understanding slow compression and decompression of frictionless soft granular matter by network analysis

We consider dense granular systems in three spatial dimensions exposed to slow compression and decompression, below, during, above and well above jamming. The evolution of granular systems under slow deformation is non-trivial and involves smooth, continuous, reversible (de)compression periods, interrupted by fast, discontinuous, irreversible transition events. These events are often, but not always, associated with rearrangements of particles and of the contact network. How many particles are involved in these transitions between two states can range from few to almost all in the system. An analysis of the force network that is built on top of the contact network is carried out using the tools of persistent homology. Results involve the observation that kinetic energy is correlated with the intensity of rearrangements, while the evolution of global mechanical measures, such as pressure, is strongly correlated with the evolution of the topological measures quantifying loops in the force network. Surprisingly, some transitions are clearly detected by persistent homology even though motion/rearrangement of particles is much weaker, i.e., much harder to detect or, in some cases, not observed at all.

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.

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.

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.

Stress driven fractionalization of vacancies in regular packings of elastic particles

Elucidating the interplay of defects and stress at the microscopic level is a fundamental physical problem that has a strong connection with materials science. Here, based on the two-dimensional crystal model, we show that the instability mode of vacancies with varying size and morphology conforms to a common scenario. A vacancy under compression is fissioned into a pair of dislocations that glide and vanish at the boundary. This neat process is triggered by the local shear stress around the vacancy. The remarkable fractionalization of vacancies creates rich modes of interaction between vacancies and other topological defects, and provides a new dimension for mechanical engineering of defects in extensive crystalline structures.

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.

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.

Stress distribution in two-dimensional silos

Simulations of a polydispersed two-dimensional silo were performed using molecular dynamics, with different numbers of grains reaching up to 64 000, verifying numerically the model derived by Janssen and also the main assumption that the walls carry part of the weight due to the static friction between grains with themselves and those with the silo's walls. We vary the friction coefficient, the radii dispersity, the silo width, and the size of grains. We find that the Janssen's model becomes less relevant as the the silo width increases since the behavior of the stresses becomes more hydrostatic. Likewise, we get the normal and tangential stress distribution on the walls evidencing the existence of points of maximum stress. We also obtained the stress matrix with which we observe zones of concentration of load, located always at a height around two thirds of the granular columns. Finally, we observe that the size of the grains affects the distribution of stresses, increasing the weight on the bottom and reducing the normal stress on the walls, as the grains are made smaller (for the same total mass of the granulate), giving again a more hydrostatic and therefore less Janssen-type behavior for the weight of the column.

Scaling laws for the mechanics of loose and cohesive granular materials based on Baxter's sticky hard spheres

We have conducted discrete element simulations (pfc3d) of very loose, cohesive, granular assemblies with initial configurations which are drawn from Baxter's sticky hard sphere (SHS) ensemble. The SHS model is employed as a promising auxiliary means to independently control the coordination number z_{c} of cohesive contacts and particle volume fraction ϕ of the initial states. We focus on discerning the role of z_{c} and ϕ for the elastic modulus, failure strength, and the plastic consolidation line under quasistatic, uniaxial compression. We find scaling behavior of the modulus and the strength, which both scale with the cohesive contact density ν_{c}=z_{c}ϕ of the initial state according to a power law. In contrast, the behavior of the plastic consolidation curve is shown to be independent of the initial conditions. Our results show the primary control of the initial contact density on the mechanics of cohesive granular materials for small deformations, which can be conveniently, but not exclusively explored within the SHS-based assembling procedure.

Effect of the microstructure on the propagation velocity of ultrasound in magnetic powders

We analyze experimentally and theoretically the sound propagation velocity of P-waves in granular media made of micrometer-size magnetite particles under an external magnetic field. The sound velocity is measured in a coherent (long-wavelength) regime of propagation after a controlled sample preparation consisting of a fluidization and the application of a magnetic field. Several different procedures are applied and result in different but reproducible particle arrangements and preferential contact orientations affecting the measured sound velocity. Interestingly, we find that the sound velocity increases when the magnetic field is applied parallel to the sound propagation direction and decreases when the magnetic field is applied perpendicular to the sound propagation direction. The observed qualitative relationship between the changes in the particle arrangement and the sound velocity is analyzed theoretically based on an effective medium theory adapted to account for the effect of the magnetic field in the preparation procedure and its influence on the medium contact fabric.

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. 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.