Modeling Cross Anisotropy in Granular Materials

A constitutive model has been developed to capture the behavior of cross-anisotropic frictional materials. The elastoplastic, single hardening model for isotropic materials serves as the basic framework. Based on the experimental results of cross-anisotropic sands in isotropic compression tests, the principal stress coordinate system is rotated such that the model operates isotropically within the rotated framework. Experimental plastic work contours on the octahedral plane are plotted for a series of true triaxial tests on dense Santa Monica Beach sand to study the effects of cross anisotropy on the evolution of yield surfaces. The amount of rotation of the yield and plastic potential surfaces decreases to zero (isotropic state) with loading. The model is constructed for cases where the principal stress and material symmetry axes are collinear and no significant rotation of principal stresses occur. The model incorporates fourteen parameters that can be determined from simple experiments, such as isotropic compression, drained triaxial compression, and triaxial extension tests. A series of true triaxial and isotropic compression tests on dense Santa Monica Beach sand are used as a basis for verification of the capabilities of the proposed model.     

Modeling Permanent Deformation of Unbound Granular Materials under Repeated Loads

Finite-element analysis on a pavement structure under traffic loads has been a viable option for researchers and designers in highway pavement design and analysis. Most of the constitutive drivers used were nonlinear elastic models defined by empirical resilient modulus equations. Few isotropic/kinematic hardening elastoplastic models were used but applying thousands of repeated load cycles became computationally expensive. In this paper, a cyclic plasticity model based on fuzzy plasticity theory is presented to model the long-term behavior of unbound granular materials under repeated loads. The discussion focuses on the model parameters that control long-term behavior such as elastic shakedown. The performance of the constitutive model is presented by comparing modeled and measured permanent strain at various numbers of load cycles. Calculated resilient modulus from the complete stress-strain curve is also discussed.     

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.     

Instability in Granular Materials: Experimental Evidence of Diffuse Mode of Failure for Loose Sands

The influence of three loading paths on the collapse of loose sand is analyzed with a particular attention paid to the onset of collapse and the mode of failure exhibited. Experimental results on conventional undrained triaxial compression tests, constant shear drained tests, as well as quasi-constant shear undrained path are presented, compared, and analyzed. It is now recognized that some collapses can occur before the Mohr-Coulomb plastic limit criterion is reached, and our recent results obtained with the new arrangement built up highlight that these collapses occur under a diffuse mode of failure. An extensive experimental series of tests shows that the first negative value of the second-order work computed using experimental data corresponds to the loss of controllability. Moreover, it is shown that the stress ratios at collapse and the corresponding mobilized angles of friction are very close for all types of tests. For similar void ratios, the onset of collapse is thus largely independent of the loading path under drained and undrained conditions but depends on a stress state to bring the material inside the unstable domain and also on the current direction of the stress increment. Indeed, it appears that the orientations of the stress increments at collapse for all tests are the same, what explains, according to the second-order work criterion, that collapse occurs at the same stress ratio. A potentially unstable domain, depending on the stress increment direction, can thus be defined.     

Elastoplastic Model for the Long-Term Behavior Modeling of Unbound Granular Materials in Flexible Pavements

The constitutive modeling of cyclic plasticity of soils has made great progress, especially in the area of sands liquefaction modeling. Nowadays, the problem of rutting of flexible pavements linked to permanent deformations occurring in the unbound layers is taken into account only by empirical formulas. This paper presents an elastoplastic model with both isotropic and kinematic hardening. The yield surface, plastic potential, and isotropic hardening are based on a model for sands, which takes into account the influence of the initial void ratio and of the mean stress on the mechanical behavior. A kinematic hardening has been added in order to take into account the mechanical behavior of the material for large cycle numbers. A complete model is then developed, simulations are presented, and comparisons with repeated load triaxial tests carried out on a subgrade soil (clayey sand), have been made. These comparisons underline the capabilities of the model to take into account the monotonic, cyclic, and ratchetting behavior of unbound materials for roads.     

Microscopic Evaluation of Strain Distribution in Granular Materials during Shear

The evolution of local strains during shear of particles of a granular material is presented in this paper. A cylindrical specimen composed of 6.5-mm spherical plastic particles was loaded under an axisymmetric triaxial loading condition. Computed tomography (CT) was used to acquire three-dimensional images of the specimen at three shearing stages. The high-resolution CT images were used to identify the 3D coordinates of 400 particles. Nine strain components (normal, shear, and rotation), rotation angles, and local dilatancy angles for particle groups were calculated, and their frequency distribution histograms are presented and discussed. It was found that there is no preferred shear direction, and the standard deviation values for shear strain components (εxy, εxz, and εyz) were almost equal for the specific test shearing stage. Shear strains as high as 25.6% were recorded for some particle groups. Furthermore, granular particle groups rotated in the 3D space with almost equal amounts of rotation strains when loaded under axisymmetric triaxial condition. Rotation strain values are very close to the corresponding shear strains. Compared to particle sliding, rotation plays a major role in the shearing resistance of granular materials. The cumulative vertical rotation angles can be as high as 38° and the horizontal rotation angles have values as high as 60°. The statistical distributions of the local dilatancy angle (ψ1) of particle groups were calculated and found to be increasing as shearing continues. The “global” dilatancy angle value is very close to the mean local ψ1 during the first stage of shearing (i.e, when global εz = ?7.3%)     

Alternative Definition of Particle Rolling in a Granular Assembly

The paper presents a definition of rolling between a pair of two-dimensional or three-dimensional particles with a compliant contact. The definition of rolling movement is based upon the shapes of the objects’ surfaces as described with differential geometry. A pseudoinverse of the surface curvatures is used for producing a rolling vector that is tangent to the two objects at their contact. Matrix expressions are presented for the efficient computation of the vector. The rolling vector is objective and is independent of the reference points that are used to track the particle motions. The definition of rolling is applied in a discrete element method simulation of the triaxial compression of three large, dense cubic assemblies: one packing of spherical particles and two packings of nonspherical particles. At small strains, particle rolling was slightly less with the nonspherical particles, but the packing with the greatest coordination number had much less rolling.     

Visualization of Crushing Evolution in Granular Materials under Compression Using DEM

Granular materials forming part of natural slopes, embankments, subgrades of foundations, and pavement structures are subjected to both static and dynamic loads during their engineering lives. As a result of these loads, particle crushing may occur. The present study demonstrates that the discrete-element method (DEM) can be used to visualize the evolution of this breakage process. In particular, the evolution of crushing in a simulated granular material subjected to uniaxial compression is presented. Even though DEM does not normally consider particle breakage, it is possible to simulate crushing by replacing one particle that has failed in tension with a combination of many particles of different sizes. The results from the simulation indicate that crushing does not develop uniformly throughout the sample, but rather concentrates in certain regions. These observations agree with experimental results of uniaxial tests conducted on sand. Other results from the simulation satisfactorily agree with experimental results previously reported by other researchers. In this way, by using a simplified failure criterion, DEM can be used to visualize and understand the evolution of granular crushing. This is something that is very difficult to do with laboratory tests alone.     

Wave Velocities in Granular Materials under Microgravity

Velocities of primary (P) and shear (S) waves in granular materials are highly dependent on confining stress. These wave velocities are related to mechanical properties of the materials such as stiffness, density, and stress history. Measurements of the wave velocities using piezoelectric sensors provide scientists and engineers a technique for nonintrusive characterization of those mechanical properties. For aerospace engineering, measuring the wave velocities under microgravity, which simulates low loading and stress conditions, has a number of potential applications. It can help the understanding of the soil mechanics and the development of appropriate materials handling technologies in extraterrestrial environments, which will be crucial to meeting NASA’s future space exploration goals. This paper presents the technique and results of experiments conducted at NASA Glenn Research Center using the 2.2?s drop tower. Velocities of P and S waves in three sizes of glass beads and one size of alumina beads were measured under initially dense or loose compaction states. It was found that under microgravity, the wave signals were significantly weaker and the velocities were much slower. The material that makes up the beads has a strong influence on the wave velocities as well. The initial compaction state also has some influence on the wave velocities.     

Fluidization in an Anaerobic EGSB Reactor: Analysis of Primary Wakes and Modeling of Sludge Blanket

Fluidization of biogranules in an anaerobic expanded granular sludge bed (EGSB) reactor is stochastic in nature and it is a function of the size distribution and the frequency of generation of flow-through gas bubbles in the reactor. Other factors that contribute to the distribution of granules along the height of the reactor are the settling characteristics of granules and the fluid velocity. A simulation was conducted in a test column to obtain a relationship between the flow-through gas and granules at different heights along the column. This relationship was combined with the pattern of gas flow through an identical EGSB reactor to create a model to predict the concentration of granules at different heights along the reactor. The model can predict well the stochastic nature of the axial distribution of granules but underestimates the number of granules at different heights. The reasons for such deviations are explained. The pattern of granule shedding from the primary wake associated with spherical cap bubbles and terminal velocities of bubbles have also been studied and modeled to estimate the maximum height of ascent of granules under isolated spherical cap bubbles. The results of this model agreed well with the experimental observations.     

Quantification and Simulation of Particle Kinematics and Local Strains in Granular Materials Using X-Ray Tomography Imaging and Discrete-Element Method

Microfeatures of granular materials have significant effects on their macrobehaviors. Unfortunately, three-dimensional (3D) quantitative measurements of microfeatures are rare in literature because of the limitations of conventional techniques in obtaining microquantities such as microdisplacements and local strains. This paper presents a new method for quantifying the particle kinematics and local strains for a soft confined compression test using X-ray computed tomography and compares the experimental measurements with the simulated results using the discrete-element method (DEM). The experimental method can identify and recognize 3D individual particles automatically, which is essential for quantifying particle kinematics and local strains. 3D DEM simulations of the soft confined compression test were performed by using spherical particles and irregular particles. The simulated global deformations and particle translations that were based on irregular particles showed better agreement with the experimental measurements than those that were based on spherical particles. The simulated movements of spherical particles were more erratic, and the material composed of spherical particles showed larger vertical contraction and radial dilation.     

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.     

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.     

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.     

Simple Nonlinear Model for Elastic Response of Cohesionless Granular Materials

A constitutive model based on hyperelasticity is proposed to capture the resilient (elastic) behavior of granular materials. Resilient behavior is a widely accepted idealization of the response of unbound granular layers of pavements, following shakedown. The coupling property of the proposed model accounts for shear dilatancy and pressure-dependent behavior of the granular materials. The model is calibrated using triaxial resilient test data obtained from the literature. A statistical comparison is made between the predictions of the proposed model and a few of the prominent models of resilient response. The proposed coupled hyperelastic model yields a significantly better fit to the experimental data. It also offers a computational efficiency when implemented in a classical nonlinear finite elemental framework.     

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.     

Observations of Stresses and Strains in a Granular Material

The use of glass ballotini as a granular material provides the opportunity to simultaneously study internal stress fields and internal fields of deformation as a sample is submitted to boundary perturbations. Digital image correlation makes use of the visible fabric of the material to deduce a field of displacements from one digital photographic image to the next. If the glass granules are immersed in a fluid having the same refractive index, then observation with polarized light exploits the photoelastic properties of the glass to reveal information about the stresses. Again, comparison of digital photographs enables changes in stress conditions from one image to the next to be discovered. Tests performed in a simple loading device which forces rotation of principal axes in parts of the granular mass are presented to demonstrate the unique potential of this dual experimental configuration.     

Probabilistic Settlement Analysis by Stochastic and Random Finite-Element Methods

The paper discusses finite element models for predicting the elastic settlement of a strip footing on a variable soil. The paper then goes on to compare results obtained in a probabilistic settlement analysis using a stochastic finite element method based on first order second moment approximations, with the random finite element method based on generation of random fields combined with Monte Carlo simulations. The paper highlights the deficiencies of probabilistic methods that are unable to properly account for spatial correlation.     

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.