Critical State of Granular Materials Based on the Sliding-Rolling Theory

By representing the assembly by a simplified column model, a constitutive theory was recently developed for a two-dimensional assembly of rods. This theory, referred to as the sliding-rolling theory, is extended in this paper to represent the triaxial stress-strain behavior of granular materials. The sliding-rolling theory provides a dilatancy rule and an expression for the slope of the line of zero dilatancy in the stress space. These rules are then combined with triaxial observations to provide a microstructural interpretation of the critical state of granular materials. According to the theory, the slope of the critical state line in the stress space depends on the interparticle friction angle and the degree of contact normal anisotropy. To verify the basic ideas of the sliding-rolling theory, numerical experiments are conducted using the discrete-element method on three-dimensional assemblies of spheres.     

Sliding and Rolling Constitutive Theory for Granular Materials

By analyzing a microelement based on four spheres, equations governing the equilibrium of the microelement are developed. By examining these equations more closely, two primary mechanisms of failure of the microelements, one based on particle sliding and the other based on particle rolling, are identified. For each primary mechanism, two separate mechanisms, one based on collapse of the microelement in the vertical direction and the other based on collapse in the horizontal direction, are recognized. With the aid of these concepts, constitutive equations are developed for a two-dimensional assembly of granular particles. The assembly is considered to consist of four-sphere microelements. Taking the plastic strain to be a consequence of the collapse of some of the microelements, equations are developed for plastic strain. The formulation yields loading criteria and flow directions. With suitable hardening rules, it is shown that the microstructural model is capable of simulating most of the salient features of the stress-strain behavior of granular materials. In particular, the stress-dilatancy relation, taking into consideration phenomena of phase transformation, and critical state failure are simulated satisfactorily.     

Evaluation of Two Homogenization Techniques for Modeling the Elastic Behavior of Granular Materials

This paper discusses the capabilities of two homogenization techniques to accurately represent the elastic behavior of granular materials considered as assemblies of randomly distributed particles. The stress-strain relationship for the assembly is determined by integrating the behavior of the interparticle contacts in all orientations, using two different homogenization methods, namely the kinematic method and the static method. The numerical predictions obtained by these two homogenization techniques are compared to results obtained during experimental studies on different granular materials. Relations between elastic constants of the assembly, interparticle properties, and fabric parameters are discussed, as well as the capabilities of the models to take into account inherent and stress-induced anisotropy for different stress conditions.     

Elastic Model for Partially Saturated Granular Materials

This paper presents the development of an elastic model for partially saturated granular materials based on micromechanical factor consideration. A granular material is considered as an assembly of particles. The stress-strain relationship for an assembly can be determined by integrating the behavior at all interparticle contacts and by using a static hypothesis, which relates the average stress of the granular assembly to a mean field of particle contact forces. As for the nonsaturated state, capillary forces at grain contacts are added to the contact forces created by an external load. These are then calculated as a function of the degree of saturation, depending on the grain size distribution and on the void ratio of the granular assembly. Hypothesizing a Hertz-Mindlin law for the grain contacts leads to an elastic nonlinear behavior of the particulate material. The prediction of the stress-strain model is compared to experimental results obtained from several different granular materials in dry, partially saturated and fully saturated states. The numerical predictions demonstrate that the model is capable of taking into account the influence of key parameters, such as degree of saturation, void ratio, and mean stress.     

Capillary Force and Water Retention between Two Uneven-Sized Particles

Capillary force and water retention between two uneven-sized spherical particles are investigated. Previous studies have been limited to systems with even-sized particles. The appropriate definition of the boundary value problem for a water lens between two uneven-sized particles is presented under the consideration of thermodynamic free energy at the microscopic level. Capillary force and water retention under the consideration of toroidal approximation are also derived for a system with two uneven-sized particles. Comparison of normalized capillary force and water retention calculated by the free energy approach and toroidal approximation are conducted. The quantitative analysis shows that for a system with two identical particles, the behavior of water retention and normalized capillary force is very similar to some recent studies by others, confirming that the toroidal approximation provides reasonably good estimations for both capillary force and water retention. For a system with uneven-sized particles, it is shown that error in normalized capillary force could be significant as the matric suction approaches zero and the particle sizes become very different. The errors for the mean curvature of the meniscus for the toroidal approximation are significant where the matric suction is near zero. Thus for soils with varying particle sizes, it may be necessary to employ the exact solution to meniscus shape in order to accurately quantify normalized capillary force and water retention. The induced normalized capillary force increases inversely with the particle size, and is generally insensitive to the water content. For soil assembly with particle size of 0.01?mm, the normalized capillary force could reach 10?kPa, whereas for soil assembly with particle size of 1?mm, the normalized capillary force is on the order of 100?Pa.     

Modeling Liquefaction by a Multimechanism Model

An anisotropic constitutive model was recently presented for describing the stress–strain behavior of granular materials with considerations for the initial and induced anisotropy. The model was developed within the framework of a microstructural theory known as the sliding–rolling theory. The resulting model falls within the definition of multimechanism models. The model was shown to satisfactorily represent the drained and undrained behaviors under monotonic loading. The framework used in the model allows extension to describe the behavior under cyclic loading, which is the subject of the present paper. Specifically, the model is further developed for representing the undrained behavior of granular materials under one- and two-way cyclic loading, some of which cause liquefaction resulting in large strain accumulations and the others lead to limited pore pressure and strain accumulations. The validity of the model is verified using triaxial data on Nevada sand.     

Micro-Macro Quantification of the Internal Structure of Granular Materials

We have attempted a multiscale quantification of the internal structure of granular materials. The internal structure of granular materials, i.e., the geometrical information on granular particles and their spatial arrangement, was described mathematically on the particle scale using Voronoi–Delaunay tessellations. These tessellations were further modified into two cell systems: a solid cell system and a void cell system, with the internal supporting structure properly reflected. By doing so, the two cell systems were geometrically and physically significant. Taking solid/void cells as the microscopic basic elements, the behavior of granular materials was expressed as the volumetric average of the microcell behavior. Macroscopically, the internal structure could be characterized by the statistical measures from the geometry of the microcells. Our approach was used to investigate the anisotropic behavior of granular materials. A study on the void cells explains how the spatial arrangement affects the strength and dilatancy of granular materials. A new anisotropic fabric tensor was defined based on the void cell anisotropy. The correlation between the anisotropic fabric tensor and the macro behavior of granular materials was verified with numerical simulations. The results showed that the new material anisotropic tensor is a more effective definition than the existing ones based on particle orientations and contact normals.     

Anisotropic Nonlinear Elastic Model for Particulate Materials

This paper presents the development of an elastic model for particulate materials based on micromechanics considerations. A particulate material is considered as an assembly of particles. The stress–strain relationship for an assembly can be determined by integrating the behavior of the interparticle contacts in all orientations and using a static hypothesis which relates the average stress of the granular assembly to a mean field of particle contact forces. Hypothesizing a Hertz–Mindlin law for the particle contacts leads to an elastic nonlinear behavior of the particulate material, we were able to determine the elastic constants of the granular assembly based on the properties of the particle contacts. The numerical predictions, compared to the results obtained during experimental studies on different granular materials, show that the model is capable of taking into account both the influence of the inherent anisotropy and the influence of the stress-induced anisotropy for different stress conditions.     

Specimen Size Effect in Discrete Element Simulations of Granular Assemblies

The paper addresses the question of whether the number of particles in a noncemented granular assembly will affect the mechanical characteristics of the assembly: its strength and stiffness. The question is answered by applying the discrete element method to assemblies of different sizes. To isolate the effect of assembly size, apart from the scatter that usually accompanies such simulations, multiple assemblies were tested. The two-dimensional assemblies had nearly identical initial porosities and fabrics, and they were all loaded in biaxial compression. Two different boundaries were tested: periodic and wall boundaries. We find that the peak compressive strength decreases with assembly size for both types of boundaries and over a range of assembly sizes that contain 256 particles to over 66,000 particles. Stiffness is only slightly reduced and only with wall boundaries. Deformation is less uniform in the larger assemblies, with deformation concentrated in a smaller fraction of the assembly area. An analysis of deformation patterning shows that at least a few thousand particles are required for realistic microband patterning.     

Three-Dimensional Digital Representation of Granular Material Microstructure from X-Ray Tomography Imaging

This paper presents the use of x-ray tomography imaging to reconstruct a three-dimensional (3D) digital representation of individual particles in a granular system. The granular system is represented by the mass center coordinates and the morphology representation of each particle. An automated procedure using pattern recognition to identify related particle cross sections in adjacent serial images was developed. Procedures to calculate quantities needed for subsequent simulation of particle behavior including the volume and the momentum of inertia of each particle are also presented. The developments described in the paper enable modeling and simulation of the behavior, and experimental observations of the particle kinematics of real microstructures of granular materials in a true 3D platform.     

Modified Direct Shear Test for Anisotropic Strength of Sand

This paper presents a simple method to estimate the directional dependency of granular soil strength using a modified shear box and a special specimen preparation procedure. This method is used to investigate the strength anisotropy of granular materials with particle shapes varying from spherical to angular. The experimental results show that the friction angle of granular materials varies with the orientation of shear plane relative to the bedding plane, and the degree of anisotropy is affected by particle shape. Comparison of the data from direct shear tests in this study with those of plane strain and torsional simple shear tests in the literature shows that direct shear using the modified direct shear box can reasonably capture the directional dependency of the friction angle for cohesionless materials.     

Influence of Particle Shape and Surface Friction Variability on Response of Rod-Shaped Particulate Media

Discrete element methods are important tools for the investigation of the mechanics of granular materials. In two dimensions, the reliability of these numerical approaches can be explored using physical tests on rod assemblies. This work highlights the importance of representing the actual distribution of rod shapes and surface friction in numerical simulations. The sensitivity of the response of hexagonally packed rods to minor changes in particle geometry and friction is investigated using a combination of laboratory tests and discrete element simulations. Laboratory test results highlight the influence of small variations in rod geometry on the global response, with the peak friction angle decreasing significantly as the standard deviation of the rod size distribution increased. Small changes in rod shape are also seen to be important. The numerical simulations indicate that the peak friction angle decreases as the standard deviation of the distribution of particle surface friction increases. This paper illustrates the way in which laboratory tests and numerical simulations can be used in a complementary manner to better understand the micromechanics of the response of granular materials.     

Micromechanical Analysis for Interparticle and Assembly Instability of Sand

Instability of granular material may lead to catastrophic events such as the gross collapse of earth structures, and thus it is an important topic in geotechnical engineering. In this paper, we adopt the micromechanics approach for constitutive modeling, in which the soil is considered an assembly of particles, and the stress-strain relationship for the assembly is determined by integrating the behavior of the interparticle contacts in all orientations. Although analyses regarding material instability have been extensively studied for a soil element at the constitutive level, it has not been considered at the interparticle contact level. Through an eigenvalue analysis, two modes of instability are identified at the local contact level: the singularity of tangential stiffness matrix and the loss of positiveness of second-order work. The constitutive model is applied to simulate drained and undrained triaxial tests on Toyoura sand of various densities under various confining pressures. The predictions are compared with experimentally measured instability at the assembly level. The modes of stability at the interparticle contact level and their relations to the overall instability of the assembly are also analyzed.     

Model for Granular Materials with Surface Energy Forces

In light of environmental differences (such as gravitational fields, surface temperatures, atmospheric pressures, etc.), the mechanical behavior of the subsurface soil on the Moon is expected to be different from that on the Earth. Before any construction on the Moon can be envisaged, a proper understanding of soil properties and its mechanical behavior in these different environmental conditions is essential. This paper investigates the possible effect of surface-energy forces on the shear strength of lunar soil. All materials, with or without a net surface charge, exhibit surface-energy forces, which act at a very short range. Although, these forces are negligible for usual sand or silty sand on Earth, they may be important for surface activated particles under extremely low lunar atmospheric pressure. This paper describes a constitutive modeling method for granular material considering particle level interactions. Comparisons of numerical simulations and experimental results on Hostun sand show that the model can accurately reproduce the overall mechanical behavior of soils under terrestrial conditions. The model is then extended to include surface-energy forces between particles in order to describe the possible behavior of lunar soil under extremely low atmospheric pressure conditions. Under these conditions, the model shows that soil has an increase of shear strength due to the effect of surface-energy forces. The magnitude of increased shear strength is in reasonable agreement with the observations of lunar soil made on the Moon’s surface.     

Generalized Model for Geosynthetic-Reinforced Granular Fill-Soft Soil with Stone Columns

This paper pertains to the development of a mechanical model to predict the behavior of a geosynthetic-reinforced granular fill over soft soil improved with stone columns. The saturated soft soil has been idealized by Kelvin–Voight model to represent its consolidation behavior. The stone columns are idealized by stiffer springs. Pasternak shear layer and rough elastic membrane represent the granular fill and geosynthetic reinforcement layer, respectively. The nonlinear behavior of the granular fill and the soft soil is considered. Effect of consolidation of the soft soil due to inclusion of the stone columns has also been included in the model. Plane strain conditions are considered for the loading and reinforced foundation soil system. An iterative finite difference scheme is applied for obtaining the solution, and results are presented in nondimensional form. Comparison between the results from the present study and the analytical solution using theory of elasticity shows reasonable agreement. The advantage of using geosynthetic reinforcement is highlighted. Results indicate that inclusion of the geosynthetic layer effectively reduces the settlement. Nonlinearity in the behavior of the soft soil and the granular fill is reduced due to the use of geosynthetic reinforcement layer.     

Input Parameters of Discrete Element Methods

This paper presents a sensitivity analysis of the input parameters of a program based on the discrete element method (DEM). Triaxial compression simulations were conducted on an assembly of ellipsoids with two particle shapes. We examine four input parameters including shear modulus of particles, density of particles, time step, and damping. Generally, these parameters are chosen by calibrating the result with certain known behavior of granular materials. In dynamic simulations, these input parameters are bounded by their physical attributions that should not be altered. However, in static simulations, they do not have the same physical implication. Validity of results may be questionable when input parameters are used without justification. A sensitivity analysis of the input parameters should shed light on this issue. In this paper, we will study the effect of the input values within the range of (1/10)–10 times the benchmark value. The benchmark values are commonly used by the writer. The results are presented against the benchmark simulation. The unbalanced forces in the simulations are kept below a prescribed value to enforce equilibrium. The result shows that the effect of all input parameters used in this paper is negligible as long as the small unbalanced forces in the system can be achieved. The runtimes are different. However, there are two simulations (one with low damping and the other with a large time step size) that cannot maintain the required small unbalanced force. In other words, equilibrium cannot be achieved for these simulations.     

Stochastic Modeling of Load-Induced Settlement in Loose Granular Materials

The response of loose cohesionless granular material to surface applied loads is investigated from the viewpoint of probabilistic mechanics of particulate media. A model is proposed that is based on the combined propagation of intergranular forces and an excess volume of voids. In this regard, it provides a bridge between earlier theories developed independently for the diffusion of stresses and for the propagation of settlements. In its general formulation, the theory can model three-dimensional, transient effects. However, the model is believed to be limited to normally consolidated or noncompacted, fully drained or dry, granular materials that do not exhibit dilatancy effects. The derived numerical modeling of steady state deflection patterns under a rigid footing is found to be in good agreement with x-ray images of laboratory model tests using noncompacted silt. The proposed theory recognizes the discrete and inherently random nature of natural granular materials such as cohesionless soils and builds upon these fundamental characteristics to predict responses of such materials to boundary applied load. This is achieved by modeling intergranular force and excess pore volume propagation as Markovian diffusion-advection processes. This approach, which departs from traditional continuum mechanics models, seems to have potential for addressing some of the challenging aspects of granular material mechanics in lunar or Martian environments.     

Normalizing Behavior of Unsaturated Granular Pavement Materials

One of the important components of a flexible pavement structure is granular material layers. Unsaturated granular pavement materials (UGPMs) in these layers influence stresses and strains throughout the pavement structure, and have a large effect on asphalt concrete fatigue and pavement rutting (two of the primary failure mechanisms for flexible pavements). The behavior of UGPMs is dependent on water content, but this effect has been traditionally difficult to quantify using either empirical or mechanistic methods. This paper presents a practical mechanistic framework for determining the behavior of UGPMs within the range of water contents, densities, and stress states likely to be encountered under field conditions. Both soil suction and generated pore pressures are determined and compared to confinement under typical field loading conditions. The framework utilizes a simple soil suction model that has three density-independent parameters, and can be determined using conventional triaxial equipment that is available in many pavement engineering laboratories.     

Visualization and Analysis of Microstructure in Three-Dimensional Discrete Numerical Models

A computer program has been developed to allow for the virtual slicing of irregularly spaced and irregularly shaped three-dimensional image data. The program was used to virtually slice three-dimensional particle assemblies from discrete element method (DEM) simulations, allowing, for the first time, direct comparison to two-dimensional slices extracted from solidified physical specimens. Based on slices obtained from the numerical specimens, it is possible to compare quantitatively numerical microstructure directly to its physical analog, which should lead to greatly improved calibrations of granular mechanics models, and could facilitate the calibration of models across all scales of interest rather than solely at specimen boundaries. Improved confidence in the ability of the DEM to realistically simulate the microstructure of granular assemblies (through improved multiscale calibration) should result in increased confidence in microstructural parameters measurable in numerical simulations but inaccessible in the laboratory. Algorithm development within the framework of the open-source Visualization Toolkit is described and performance of the algorithm is quantified for two platforms. Results from virtual slices of a test assembly with regular particle packing are verified against known analytical solutions. A slice of a more complex assembly comprised of nearly 40,000 spheres is quantified statistically and compared to an analogous slice from a physical specimen of uniform sand.