Aggregation and caking processes of granular materials:continuum model and simulation with application to sugar

Aggregation and caking of particles are common severe problems in many operations and processing of granular materials, where granulated sugar is an important example. Prevention of aggregation and caking of granular materials requires a good understanding of moisture migration and caking mechanisms. In this paper, the modeling of solid bridge formation between particles is introduced, based on moisture migration of atmospheric moisture into containers packed with granular materials through vapor evaporation and condensation. A model for the caking process is then developed, based on the growth of liquid bridges (during condensation), and their hardening and subsequent creation of solid bridges (during evaporation). The predicted caking strengths agree well with some available experimental data on granulated sugar under storage conditions.

Shear Strength of Granular Materials, 1967

A theory was developed to allow the separate determination of the effects of the interparticle friction and interlocking of particles on the shearing resistance and deformational behavior of granular materials. The derived parameter, angle of solid friction, is independent of the type of shear test, stress history, porosity and the level of confining pressure, and depends solely upon the nature of the particle surface. The theory was tested against published data concerning the performance of plane strain, triaxial compression and extension tests on cohesionless soils. The theory also was applied to isotropically consolidated undrained triaxial tests on three crushed limestones prepared by the authors using vibratory compaction. The authors concluded that, (1) the theory allowed the determination of solid friction between particles which was found to depend solely on the nature of the particle surface, (2) the separation of frictional and volume change components of shear strength of granular materials qualitatively corroborated the postulated mechanism of deformation (sliding and rolling of groups of particles over other similar groups with resulting dilatancy of specimen), (3) the influence of void ratio, gradation confining pressure, stress history and type of shear test on shear strength is reflected in values of the omega parameter, and (4) calculation of the coefficient of solid friction allows the establishment of the lower limit of the shear strength of a granular material.

Behavior of Granular Materials Under Triaxial Compression with Pulsating Deviator Stress, HR-131, Special Report

Behavior of granular material subjected to repeated load triaxial compression tests is characterized by a model based on rate process theory. Starting with the Arrhenius equation from chemical kinetics, the relationship of temperature, shear stress, normal stress and volume change to deformation rate is developed. The proposed model equation includes these factors as a product of exponential terms. An empirical relationship between deformation and the cube root of the number of stress applications at constant temperature and normal stress is combined with the rate equation to yield an integrated relationship of temperature, deviator stress, confining pressure and number of deviator stress applications to axial strain. The experimental program consists of 64 repeated load triaxial compression tests, 52 on untreated crushed stone and 12 on the same crushed stone material treated with 4% asphalt cement. Results were analyzed with multiple linear regression techniques and show substantial agreement with the model equations. Experimental results fit the rate equation somewhat better than the integrated equation when all variable quantities are considered. The coefficient of shear temperature gives the activation enthalpy, which is about 4.7 kilocalories/mole for untreated material and 39.4 kilocalories/mole for asphalt-treated material. This indicates the activation enthalpy is about that of the pore fluid. The proportionality coefficient of deviator stress may be used to measure flow unit volume. The volumes thus determined for untreated and asphalt-treated material are not substantially different. This may be coincidental since comparison with flow unit volumes reported by others indicates flow unit volume is related to gradation of untreated material. The flow unit volume of asphalt-treated material may relate to asphalt cement content. The proposed model equations provide a more rational basis for further studies of factors affecting deformation of granular materials under stress similar to that in pavement subjected to transient traffic loads.

Micromechanical insight into the undrained instability of granular materials: deformation characteristics of geomaterials.

There is no agreement between experimental researchers whether the point where a granular material responds with a large change of stresses, strains or excess pore water pressure given a prescribed small input of some of the same variables defines a straight line or a curve in the stress space. This line, known as the instability line, may also vary in shape and position if the onset of instability is measured from drained or undrained triaxial tests. Failure of granular materials, which might be preceded by the onset of instability, is a subject that the geotechnical engineers have to deal with in the daily practice, and generally speaking it is associated to different phenomena observed not only in laboratory tests but also in the field. Examples of this are the liquefaction of loose sands subjected to undrained loading conditions and the diffuse instability under drained loading conditions. This research presents results of DEM simulations of undrained triaxial tests with the aim of studying the influence of stress history and relative density on the onset of instability in granular materials. Micro-mechanical analysis including the evolution of coordination numbers and fabric tensors is performed aiming to gain further insight on the particle-scale interactions that underlie the occurrence of this instability. In addition to provide a greater understanding, the results presented here may be useful as input for macro-scale constitutive models that enable the prediction of the onset of instability in boundary value problems.

Exploring the effects of stress history on the drained and undrained cyclic behaviour of granular materials.

This paper presents the results of 3D DEM simulations of granular materials subject to cyclic loading. While both the drained and undrained conditions are considered, the effects of depositional history and consolidation stress history on the stress-strain response are specifically evaluated. It is demonstrated that the different stress histories have a significant effect on soil response and that such effects can be attributed to differences in the initial particle arrangement (fabric).

Properties of granular analogue model materials: A community wide survey

We report the material properties of 26 granular analogue materials used in 14 analogue modelling laboratories. We determined physical characteristics such as bulk density, grain size distribution, and grain shape, and performed ring shear tests to determine friction angles and cohesion, and uniaxial compression tests to evaluate the compaction behaviour. Mean grain size of the materials varied between c. 100 and 400 μm. Analysis of grain shape factors shows that the four different classes of granular materials (14 quartz sands, 5 dyed quartz sands, 4 heavy mineral sands and 3 size fractions of glass beads) can be broadly divided into two groups consisting of 12 angular and 14 rounded materials. Grain shape has an influence on friction angles, with most angular materials having higher internal friction angles (between c. 35° and 40°) than rounded materials, whereas well-rounded glass beads have the lowest internal friction angles (between c. 25° and 30°). We interpret this as an effect of intergranular sliding versus rolling. Most angular materials have also higher basal friction angles (tested for a specific foil) than more rounded materials, suggesting that angular grains scratch and wear the foil. Most materials have an internal cohesion in the order of 20–100 Pa except for well-rounded glass beads, which show a trend towards a quasi-cohesionless (C < 20 Pa) Coulomb-type material. The uniaxial confined compression tests reveal that rounded grains generally show less compaction than angular grains. We interpret this to be related to the initial packing density after sifting, which is higher for rounded grains than for angular grains. Ring-shear test data show that angular grains undergo a longer strain-hardening phase than more rounded materials. This might explain why analogue models consisting of angular grains accommodate deformation in a more distributed manner prior to strain localisation than models consisting of rounded grains.

Application of a continuum theory to vertical vibrations of a layer of granular material

Many interesting phenomena have been observed in layers of granular materials subjected to vertical oscillations; these include the formation of a variety of standing wave patterns, and the occurrence of isolated features called oscillons, which alternately form conical heaps and craters oscillating at one-half of the forcing frequency. No continuum-based explanation of these phenomena has previously been proposed. We apply a continuum theory, termed the double-shearing theory, which has had success in analyzing various problems in the flow of granular materials, to the problem of a layer of granular material on a vertically vibrating rigid base undergoing vertical oscillations in plane strain. There exists a trivial solution in which the layer moves as a rigid body. By investigating linear perturbations of this solution, we find that at certain amplitudes and frequencies this trivial solution can bifurcate. The time dependence of the perturbed solution is governed by Mathieu’s equation, which allows stable, unstable and periodic solutions, and the observed period-doubling behaviour. Several solutions for the spatial velocity distribution are obtained; these include one in which the surface undergoes vertical velocities that have sinusoidal dependence on the horizontal space dimension, which corresponds to the formation of striped standing waves, and is one of the observed patterns. An alternative continuum theory of granular material mechanics, in which the principal axes of stress and rate-of-deformation are coincident, is shown to be incapable of giving rise to similar instabilities.

Free surface problems for static Coulomb-Mohr granular solids

This paper is concerned with some plane strain and axially symmetric free surface problems which arise in the study of static granular solids that satisfy the Coulomb-Mohr yield condition. Such problems are inherently nonlinear, and hence difficult to attack analytically. Given a Coulomb friction condition holds on a solid boundary, it is shown that the angle a free surface is allowed to attach to the boundary is dependent only on the angle of wall friction, assuming the stresses are all continuous at the attachment point, and assuming also that the coefficient of cohesion is nonzero. As a model problem, the formation of stable cohesive arches in hoppers is considered. This undesirable phenomena is an obstacle to flow, and occurs when the hopper outlet is too small. Typically, engineers are concerned with predicting the critical outlet size for a given hopper and granular solid, so that for hoppers with outlets larger than this critical value, arching cannot occur. This is a topic of considerable practical interest, with most accepted engineering methods being conservative in nature. Here, the governing equations in two limiting cases (small cohesion and high angle of internal friction) are considered directly. No information on the critical outlet size is found; however solutions for the shape of the free boundary (the arch) are presented, for both plane and axially symmetric geometries.

Classical and Cosserat plasticity and viscoplasticity models for slow granular flow

We consider models for the rheology of dense, slowly deforming granular materials based of classical and Cosserat plasticity, and their viscoplastic extensions that account for small but finite particle inertia. We determine the scale for the viscosity by expanding the stress in a dimensionless parameter that is a measure of the particle inertia. We write the constitutive relations for classical and Cosserat plasticity in stress-explicit form. The viscoplastic extensions are made by adding a rate-dependent viscous stress to the plasticity stress. We apply the models to plane Couette flow, and show that the classical plasticity and viscoplasticity models have features that depart from experimental observations; the prediction of the Cosserat viscoplasticity model is qualitatively similar to that of Cosserat plasticity, but the viscosities modulate the thickness of the shear layer.

Experimental evidence of the kinematic Cosserat effect in dense granular flows

To capture shear localization in the flow of dense granular materials in a continuum description, it has earlier been proposed that granular materials be treated as Cosserat, or micropolar, continua. Here, we provide experimental verification of the kinematic Cosserat effect, or the deviation of the particle spin from the material spin induced by the velocity gradient. Contrary to earlier belief, we find this effect to be sizable even outside the shear layers. Remarkably, the particles and material elements spin in opposite directions in flow through a hopper.

Analysis of shear bands in slow granular flows using a frictional Cosserat model

The slow flow of granular materials is often marked by the existence of narrow shear layers, adjacent to large regions that suffer little or no deformation. This behaviour, in the regime where shear stress is generated primarily by the frictional interactions between grains, has so far eluded theoretical description. In this paper, we present a rigid-plastic frictional Cosserat model that captures thin shear layers by incorporating a microscopic length scale. We treat the granular medium as a Cosserat continuum, which allows the existence of localised couple stresses and, therefore, the possibility of an asymmetric stress tensor. In addition, the local rotation is an independent field variable and is not necessarily equal to the vorticity. The angular momentum balance, which is implicitly satisfied for a classical continuum, must now be solved in conjunction with the linear momentum balances. We extend the critical state model, used in soil plasticity, for a Cosserat continuum and obtain predictions for flow in plane and cylindrical Couette devices. The velocity profile predicted by our model is in qualitative agreement with available experimental data. In addition, our model can predict scaling laws for the shear layer thickness as a function of the Couette gap, which must be verified in future experiments. Most significantly, our model can determine the velocity field in viscometric flows, which classical plasticity-based model cannot.

Experimental Studies on the Mechanics of Cohesive Frictional Granular Media

The mechanical behaviour of cohesive-frictional granular materials is a combination of the strength pervading as intergranular friction (represented as an angle of internal friction - Phi), and the cohesion (C) between these particles. Most behavioral or constitutive models of this class of granular materials comprise of a cohesion and frictional component with no regard to the length scale i.e. from the micro structural models through the continuum models. An experimental study has been made on a model granular material, viz. angular sand with different weights of binding agents (varying degrees of cohesion) at multiple length scales to physically map this phenomenon. Cylindrical specimen of various diameters - 10, 20, 38, 100, 150 mm (and with an aspect ratio of 2) are reconstituted with 2, 4 and 8% by weight of a binding agent. The magnitude of this cohesion is analyzed using uniaxial compression tests and it is assumed to correspond to the peak in the normalized stress-strain plot. Increase in the cohesive strength of the material is seen with increasing size of the specimen. A possibility of ``entanglement'' occurring in larger specimens is proposed as a possible reason for deviation from a continuum framework.

Force chains, friction, and flow: Behavior of granular media across length scales

We study the behavior of granular materials at three length scales. At the smallest length scale, the grain-scale, we study inter-particle forces and "force chains". Inter-particle forces are the natural building blocks of constitutive laws for granular materials. Force chains are a key signature of the heterogeneity of granular systems. Despite their fundamental importance for calibrating grain-scale numerical models and elucidating constitutive laws, inter-particle forces have not been fully quantified in natural granular materials. We present a numerical force inference technique for determining inter-particle forces from experimental data and apply the technique to two-dimensional and three-dimensional systems under quasi-static and dynamic load. These experiments validate the technique and provide insight into the quasi-static and dynamic behavior of granular materials.

At a larger length scale, the mesoscale, we study the emergent frictional behavior of a collection of grains. Properties of granular materials at this intermediate scale are crucial inputs for macro-scale continuum models. We derive friction laws for granular materials at the mesoscale by applying averaging techniques to grain-scale quantities. These laws portray the nature of steady-state frictional strength as a competition between steady-state dilation and grain-scale dissipation rates. The laws also directly link the rate of dilation to the non-steady-state frictional strength.

At the macro-scale, we investigate continuum modeling techniques capable of simulating the distinct solid-like, liquid-like, and gas-like behaviors exhibited by granular materials in a single computational domain. We propose a Smoothed Particle Hydrodynamics (SPH) approach for granular materials with a viscoplastic constitutive law. The constitutive law uses a rate-dependent and dilation-dependent friction law. We provide a theoretical basis for a dilation-dependent friction law using similar analysis to that performed at the mesoscale. We provide several qualitative and quantitative validations of the technique and discuss ongoing work aiming to couple the granular flow with gas and fluid flows.

Microscopic Origin of Macroscopic Strength in Granular Media: A Numerical and Analytical Approach

Constitutive modeling in granular materials has historically been based on macroscopic experimental observations that, while being usually effective at predicting the bulk behavior of these type of materials, suffer important limitations when it comes to understanding the physics behind grain-to-grain interactions that induce the material to macroscopically behave in a given way when subjected to certain boundary conditions.

The advent of the discrete element method (DEM) in the late 1970s helped scientists and engineers to gain a deeper insight into some of the most fundamental mechanisms furnishing the grain scale. However, one of the most critical limitations of classical DEM schemes has been their inability to account for complex grain morphologies. Instead, simplified geometries such as discs, spheres, and polyhedra have typically been used. Fortunately, in the last fifteen years, there has been an increasing development of new computational as well as experimental techniques, such as non-uniform rational basis splines (NURBS) and 3D X-ray Computed Tomography (3DXRCT), which are contributing to create new tools that enable the inclusion of complex grain morphologies into DEM schemes.

Yet, as the scientific community is still developing these new tools, there is still a gap in thoroughly understanding the physical relations connecting grain and continuum scales as well as in the development of discrete techniques that can predict the emergent behavior of granular materials without resorting to phenomenology, but rather can directly unravel the micro-mechanical origin of macroscopic behavior.

In order to contribute towards closing the aforementioned gap, we have developed a micro-mechanical analysis of macroscopic peak strength, critical state, and residual strength in two-dimensional non-cohesive granular media, where typical continuum constitutive quantities such as frictional strength and dilation angle are explicitly related to their corresponding grain-scale counterparts (e.g., inter-particle contact forces, fabric, particle displacements, and velocities), providing an across-the-scale basis for better understanding and modeling granular media.

In the same way, we utilize a new DEM scheme (LS-DEM) that takes advantage of a mathematical technique called level set (LS) to enable the inclusion of real grain shapes into a classical discrete element method. After calibrating LS-DEM with respect to real experimental results, we exploit part of its potential to study the dependency of critical state (CS) parameters such as the critical state line (CSL) slope, CSL intercept, and CS friction angle on the grain's morphology, i.e., sphericity, roundness, and regularity.

Finally, we introduce a first computational algorithm to ``clone'' the grain morphologies of a sample of real digital grains. This cloning algorithm allows us to generate an arbitrary number of cloned grains that satisfy the same morphological features (e.g., roundness and aspect ratio) displayed by their real parents and can be included into a DEM simulation of a given mechanical phenomenon. In turn, this will help with the development of discrete techniques that can directly predict the engineering scale behavior of granular media without resorting to phenomenology.