Micromechanical behavior of granular materials with inherent anisotropy under cyclic loading using 2D DEM

This paper presents the micromechanical behavior of granular materials due to different initial inherent anisotropic conditions during cyclic loading using the discrete element method (DEM). Oval particles were used to model the samples. Three samples, with three different inherent anisotropic conditions based on the particle’s bedding direction, were prepared and subjected to biaxial cyclic loading. The differences in the inherent anisotropic conditions of the samples affect the stress–strain-dilative behavior of granular materials. The width of the stress–strain cyclic loops decreases as the preferred bedding angle changes from vertical to horizontal. Contact fabric evolution is found to be dependent on the inherent anisotropic fabric of the sample during loading and unloading. The fabric anisotropy is dominant for horizontal particle bedding at the end of loading and for vertical particle bedding at the end of unloading. A change in fabric anisotropy is observed only for the first few loading–unloading cycles for the given conditions depicted in the present study.     

Fabric,force and strength anisotropies in granular materials: a micromechanical insight

In micromechanics, the stress–force–fabric (SFF) relationship is referred to as an analytical expression linking the stress state of a granular material with microparameters on contact forces and material fabric. This paper employs the SFF relationship and discrete element modelling to investigate the micromechanics of fabric, force and strength anisotropies in two-dimensional granular materials. The development of the SFF relationship is briefly summarized while more attention is placed on the strength anisotropy and deformation non-coaxiality. Due to the presence of initial anisotropy, a granular material demonstrates a different behaviour when the loading direction relative to the direction of the material fabric varies. Specimens may go through various paths to reach the same critical state at which the fabric and force anisotropies are coaxial with the loading direction. The critical state of anisotropic granular material has been found to be independent of the initial fabric. The fabric anisotropy and the force anisotropy approach their critical magnitudes at the critical state. The particle-scale data obtained from discrete element simulations of anisotropic materials show that in monotonic loading, the principal force direction quickly becomes coaxial with the loading direction (i.e. the strain increment direction as applied). However, material fabric directions differ from the loading direction and they only tend to be coaxial at a very large shear strain. The degree of force anisotropy is in general larger than that of fabric anisotropy. In comparison with the limited variation in the degree of force anisotropy with varying loading directions, the fabric anisotropy adapts in a much slower pace and demonstrates wider disparity in the evolution in the magnitude of fabric anisotropy. The difference in the fabric anisotropy evolution has a more significant contribution to strength anisotropy than that of force anisotropy. There are two key parameters that control the degree of deformation non-coaxiality in granular materials subjected to monotonic shearing: the ratio between the degrees of fabric anisotropy and that of force anisotropy and the angle between the principal fabric direction and the applied loading direction.     

Macroscopic model with anisotropy based on micro–macro information

Physical experiments can characterize the elastic response of granular materials in terms of macroscopic state variables, namely volume (packing) fraction and stress, while the microstructure is not accessible and thus neglected. Here, by means of numerical simulations, we analyze dense, frictionless granular assemblies with the final goal to relate the elastic moduli to the fabric state, i.e., to microstructural averaged contact network features as contact number density and anisotropy. The particle samples are first isotropically compressed and then quasi-statically sheared under constant volume (undrained conditions). From various static, relaxed configurations at different shear strains, infinitesimal strain steps are applied to “measure” the effective elastic response; we quantify the strain needed so that no contact and structure rearrangements, i.e. plasticity, happen. Because of the anisotropy induced by shear, volumetric and deviatoric stresses and strains are cross-coupled via a single anisotropy modulus, which is proportional to the product of deviatoric fabric and bulk modulus (i.e., the isotropic fabric). Interestingly, the shear modulus of the material depends also on the actual deviatoric stress state, along with the contact configuration anisotropy. Finally, a constitutive model based on incremental evolution equations for stress and fabric is introduced. By using the previously measured dependence of the stiffness tensor (elastic moduli) on the microstructure, the theory is able to predict with good agreement the evolution of pressure, shear stress and deviatoric fabric (anisotropy) for an independent undrained cyclic shear test, including the response to reversal of strain.     

Fabric rates of elliptical particle assembly in monotonic and cyclic simple shear tests: a numerical study

This paper aims to investigate the evolutions of microscopic structures of elliptical particle assemblies in both monotonic and cyclic constant volume simple shear tests using the discrete element method. Microscopic structures, such as particle orientations, contact normals and contact forces, were obtained from the simulations. Elliptical particles with the same aspect ratio (1.4 and 1.7 respectively for the two specimens) were generated with random particle directions, compacted in layers, and then precompressed to a low pressure one-dimensionally to produce an inherently anisotropic specimen. The specimens were sheared in two perpendicular directions (shear mode I and II) in a strain-rate controlled way so that the effects of inherent anisotropy can be examined. The anisotropy of particle orientation increases and the principal direction of particle orientation rotates with the shearing of the specimen in the monotonic tests. The shear mode can affect the way fabric anisotropy rate of particle orientation responds to shear strain as a result of the initial anisotropy. The particle aspect ratio exhibits quantitative influence on some fabric rates, including particle orientation, contact normal and sliding contact normal. The fabric rates of contact normal, sliding contact normal, contact force, strong and weak contact forces fluctuate dramatically around zero after the shear strain exceeds 4 % in the monotonic tests and throughout the cyclic tests. Fabric rates of contact normals and forces are much larger than that of particle orientation. The particle orientation based fabric tensor is harder to evolve than the contact normal or contact force based because the reorientation of particles is more difficult than that of contacts.     

A local constitutive model with anisotropy for ratcheting under 2D axial-symmetric isobaric deformation

A local constitutive model for anisotropic granular materials is introduced and applied to isobaric (homogeneous) axial-symmetric deformation. The simplified model (in the coordinate system of the bi-axial box) involves only scalar values for hydrostatic and shear stresses, for the volumetric and shear strains as well as for the new ingredient, the anisotropy modulus. The non-linear constitutive evolution equations that relate stress and anisotropy to strain are inspired by observations from discrete element method (DEM) simulations. For the sake of simplicity, parameters like the bulk and shear modulus are set to constants, while the shear stress ratio and the anisotropy evolve with different rates to their critical state limit values when shear deformations become large. When applied to isobaric deformation in the bi-axial geometry, the model shows ratcheting under cyclic loading. Fast and slow evolution of the anisotropy modulus with strain. Lead to dilatancy and contractancy, respectively. Furthermore, anisotropy acts such that it works “against” the strain/stress, e.g., a compressive strain builds up anisotropy that creates additional stress acting against further compression.     

Investigating the micromechanical evolutions within inherently anisotropic granular materials using discrete element method

This paper investigates the mechanical behavior of inherently-anisotropic granular materials from macroscopic and microscopic points of view. The study is achieved by simulating biaxial compression tests performed on granular assemblies by using numerical discrete element method. In the same category of numerical studies found in the literature, the simulations were performed by considering elliptical/oval particles. In the present study, however, the shape of particles is considered as convex polygons, which mostly resembles real sand grains. Particle assemblies with four different bedding angles were tested. Similar to what observed in experiment, inherent anisotropy has a significant effect on macroscopic mechanical behavior of granular materials. The shear strength and dilative behavior of assemblies were found to decrease as the bedding angle increases. Evolution of the microstructure of all samples and the influence of bedding angle on the fabric and force anisotropy during loading process were investigated. It is seen that the microscopic evolutions in the fabric can justify well the macroscopic behavior of granular assemblies. It is found that the long axis of particles tend to be inclined perpendicular to the loading axis, which results in generating more stable column-like microstructures in order to transfer the applied load. Moreover, the number of contacts as well as the magnitude of forces among particles varies in different directions during the loading process and the initial anisotropy condition totally evolves due to the induced anisotropy within samples.     

Strength of non-spherical particles with anisotropic geometries under triaxial and shearing loading configurations

This paper presents a study on the macroscopic shear strength characteristics of granular assemblies with three- dimensional complex-shaped particles. Different assemblies are considered, with both isotropic and anisotropic particle geometries. The study is conducted using the discrete element method (DEM), with so-called sphero-polyhedral particles, and simulations of mechanical true triaxial tests for a range of Lode angles and confining pressures. The observed mathematical failure envelopes are investigated in the Haigh–Westergaard stress space, as well as on the deviatoric-mean pressure plane. It is verified that the DEM with non-spherical particles produces results that are qualitatively similar to experimental data and previous numerical results obtained with spherical elements. The simulations reproduce quite well the shear strength of assemblies of granular media, such as higher strength during compression than during extension. In contrast, by introducing anisotropy at the particle level, the shear strength parameters are greatly affected, and an isotropic failure criterion is no longer valid. It is observed that the strength of the anisotropic assembly depends on the direction of loading, as observed for real soils. Finally simulations on a virtual shearing test show how the velocity profile within the shear band is also affected by the grain’s shape.     

The influence of particle-size distribution on critical state behavior of spherical and non-spherical particle assemblies

This paper presents an investigation into the effects of particle-size distribution on the critical state behavior of granular materials using discrete element method (DEM) simulations on both spherical and non-spherical particle assemblies. A series of triaxial test DEM simulations examine the influence of particle-size distribution (PSD) and particle shape, which were independently assessed in the analyses presented. Samples were composed of particles with varying shapes characterized by overall regularity (OR) and different PSDs. The samples were subjected to the axial compression through different loading schemes: constant volume, constant mean effective stress, and constant lateral stress. All samples were sheared to large strains to ensure that a critical state was reached. Both the macroscopic and microscopic behaviors in these tests are discussed here within the framework of the anisotropic critical state theory. It is shown that both OR and PSD may affect the response of the granular assemblies in terms of the stress–strain relations, dilatancy, and critical state behaviors. For a given PSD, both the shear strength and fabric norm decrease with an increase in OR. The critical state angle of shearing resistance is highly dependent on particle shape. In terms of PSD, uniformly distributed assemblies mobilize higher shear strength and experience more dilative responses than specimens with a greater variation of particle sizes. The position of the critical state line in the e–p′ space is also affected by PSD. However, the effects of PSD on critical strength and evolution of fabric are negligible. These findings highlight the importance of particle shape and PSD that should be included in the development of constitutive models for granular materials.     

Micromechanical behaviour of Si-based particulate assemblies: A comparative study using DEM and atomistic simulations

In this paper, we investigate the micromechanical behaviour of Si-based particulate systems subjected to tri-axial compression loading. The investigations are based on three-dimensional discrete element modelling (DEM) and simulations. At first, we compare the variation of mean compressive stress for a silicon assembly subjected to tri-axial compression, predicted at two different scales: at the particulate scale, using the DEM simulation (mono-dispersed particulates) and at the atomistic scale using the molecular dynamics (MD) simulation results for silicon mono-crystal reported by Mylvaganam and Zhang (2003) [K. Mylvaganam, L. Zhang, Key Eng. Mater. 233–236 (2003) 615–620]. Both the simulation methods considered the silicon assembly subjected to an identical (tri-axial) loading condition. We observed a good qualitative agreement between the DEM and MD simulation results for the mean compressive stress when the assembly was subjected to small volumetric strain. However, at large volumetric strain, the mean stress of the silicon assembly predicted from MD simulation did not scale-up with the DEM results. This discrepancy could be due to that MD simulation is only valid for particle contacts, which are independent of one another and does not consider the inherent ‘discrete’ nature of particulates and the induced anisotropy prevailing at particulate scale. The micromechanical behaviour of particulate assemblies strongly depends on the inherent discrete nature of the particles, their single-particle properties and the induced anisotropy during mechanical loading. At the second stage, using DEM, we present the evolution of macroscopic compressive stress and several micromechanical features for four cases of the commonly used Si based poly-dispersed particulate assemblies (Si, SiC,Si₃N₄ and SiO₂) under tri-axial compression loading. We also present the evolution of several other phenomena occurring at particulate scale, such as the energy dissipation characteristics due to sliding contacts and the features of fabric structures developed during mechanical loading. The study shows that the single-particle properties of the Si based assemblies considered here significantly affect the micromechanical behaviour of the assemblies and DEM is a powerful tool to get insights on the internal behaviour of discrete particulates under mechanical loading.     

On micromechanical characteristics of the critical state of two-dimensional granular materials

In micromechanics of quasi-static deformation of granular materials, relationships are investigated between the macro-scale, continuum-mechanical characteristics, and the micro-scale characteristics at the particle and interparticle contact level. An important micromechanical quantity is the fabric tensor that reflects the distribution of contact orientations. It also contains information on the coordination number, i.e. the average number of contacts per particle. Here, the focus is on characteristics of the critical state in the two-dimensional case. Critical state soil mechanics is reviewed from the micromechanical viewpoint. Two-dimensional discrete element method (DEM) simulations have been performed with discs from a fairly narrow particle-size distribution. Various values for the interparticle friction coefficient and for the confining pressure have been considered to investigate the effect of these quantities on critical state characteristics (shear strength, packing fraction, coordination number and fabric anisotropy). Results from these DEM simulations show that a limiting fabric state exists at the critical state, which is geometrical in origin. The contact network tessellates the assembly into loops that are formed by contacts. For each loop, a symmetrical loop tensor is defined, based on its contact normals. This loop tensor reflects the shape of the loop. An orientation is associated with each loop, based on its loop tensor. At the critical state, the frequencies with which loops with different number of sides occur depend on the coordination number. At the critical state, these loops have, on average, the following universal characteristics, i.e. independent of the coordination number: (1) loops with the same number of sides and orientation have identical anisotropy of the loop tensor, (2) the anisotropy of the loop tensor depends linearly on the number of sides of the loop, (3) the distribution of loop orientations is identical, (4) Lewis’s law for the loop areas, which is a linear relation between the number of sides of loops and their area, is satisfied (not exclusively at the critical state) and (5) the areas of the loops do not depend on their orientation.     

A coupled FEM/DEM approach for hierarchical multiscale modelling of granular media

A hierarchical multiscale framework is proposed to model the mechanical behaviour of granular media. The framework employs a rigorous hierarchical coupling between the FEM and the discrete element method (DEM). To solve a BVP, the FEM is used to discretise the macroscopic geometric domain into an FEM mesh. A DEM assembly with memory of its loading history is embedded at each Gauss integration point of the mesh to serve as the representative volume element (RVE). The DEM assembly receives the global deformation at its Gauss point from the FEM as input boundary conditions and is solved to derive the required constitutive relation at the specific material point to advance the FEM computation. The DEM computation employs simple physically based contact laws in conjunction with Coulomb's friction for interparticle contacts to capture the loading‐history dependence and highly nonlinear dissipative response of a granular material. The hierarchical scheme helps to avoid the phenomenological assumptions on constitutive relation in conventional continuum modelling and retains the computational efficiency of FEM in solving large‐scale BVPs. The hierarchical structure also makes it ideal for distributed parallel computing to fully unleash its predictive power. Importantly, the framework offers rich information on the particle level with direct link to the macroscopic material response, which helps to shed lights on cross‐scale understanding of granular media. The developed framework is first benchmarked by a simulation of single‐element drained test and is then applied to the predictions of strain localisation for sand subject to monotonic biaxial compression, as well as the liquefaction and cyclic mobility of sand in cyclic simple shear tests. It is demonstrated that the proposed method may reproduce interesting experimental observations that are otherwise difficult to be captured by conventional FEM or pure DEM simulations, such as the inception of shear band under smooth symmetric boundary conditions, non‐coaxial granular response, large dilation and rotation at the edges of shear band and critical state reached within the shear band. Copyright © 2014 John Wiley & Sons, Ltd.     

The influence of inter-particle friction and the intermediate stress ratio on soil response under generalised stress conditions

Previous research studies have used either physical experiments or discrete element method (DEM) simulations to explore, independently, the influence of the coefficient of inter-particle friction (μ) and the intermediate stress ratio (b) on the behaviour of granular materials. DEM simulations and experiments using photoelasticity have shown that when an anisotropic stress condition is applied to a granular material, strong force chains or columns of contacting particles transmitting relatively large forces, form parallel to the major principal stress orientation. The combined effects of friction and the intermediate stress ratio upon the resistance of these strong force chains to collapse (buckling failure) are considered here using data from an extensive set of DEM simulations including triaxial and true triaxial compression tests. For all tests both?μ and b affected the macro- and micro-scale response, however the mechanisms whereby the force chain stability was improved differ. While friction clearly enhances the inherent stability of the strong force chains, the intermediate stress ratio affects the contact density and distribution of orthogonal contacts that provide lateral support to the force chains.     

Macro deformation and micro structure of 3D granular assemblies subjected to rotation of principal stress axes

This paper presents a numerical investigation on the behavior of three dimensional granular materials during continuous rotation of principal stress axes using the discrete element method. A dense specimen has been prepared as a representative element using the deposition method and subjected to stress rotation at different deviatoric stress levels. Significant plastic deformation has been observed despite that the principal stresses are kept constant. This contradicts the classical plasticity theory, but is in agreement with previous laboratory observations on sand and glass beads. Typical deformation characteristics, including volume contraction, deformation non-coaxiality, have been successfully reproduced. After a larger number of rotational cycles, the sample approaches the ultimate state with constant void ratio and follows a periodic strain path. The internal structure anisotropy has been quantified in terms of the contact-based fabric tensor. Rotation of principal stress axes densifies the packing, and leads to the increase in coordination numbers. A cyclic rotation in material anisotropy has been observed. The larger the stress ratio, the structure becomes more anisotropic. A larger fabric trajectory suggests more significant structure re-organization when rotating and explains the occurrence of more significant strain rate. The trajectory of the contact-normal based fabric is not centered in the origin, due to the anisotropy in particle orientation generated during sample generation which is persistent throughout the shearing process. The sample sheared at a lower intermediate principal stress ratio \((b=0.0)\) has been observed to approach a smaller strain trajectory as compared to the case \(b=0.5\), consistent with a smaller fabric trajectory and less significant structural re-organisation. It also experiences less volume contraction with the out-of plane strain component being dilative.     

Fabric tensor based boundary element analysis of porous solids

In this paper, the boundary element analysis of porous solids (sintered materials, foams, etc.) is studied utilizing a fabric tensor. The fabric tensor provides a measure of anisotropy in the solid, as well as information concerning the geometry and distribution of the pores. The homogenized, orthotropic elastic properties of a porous solid can then be predicted with the fabric tensor. To illustrate the analysis, the effect of porosity on a trabecular bone-titanium bimaterial is studied. The boundary element analysis uses an anisotropic, bimaterial Green's function so the interface does not require discretization. It is shown that the anisotropic Stroh variables are independent of the structural density and dependent on the eigenvalues of the fabric tensor. An example calculation is presented where the effect of porosity on the in-plane maximum shear stress in a trabecular bone-titanium bimaterial is substantial.     

A 3D fourth order fabric tensor approach of anisotropy in granular media

We propose a micromechanical approach for granular media, with a particular account of the texture-induced anisotropy and of the strain localization rule. The approach is mainly based on the consideration of a fourth order fabric tensor able to capture general anisotropy which can be induced by complex distribution of contacts. Incorporation of this fourth order fabric tensor in a suitable homogenization scheme allows to determine the corresponding macroscopic elastic properties of the granular material. For this purpose, in addition to the classical Voigt upper bound, a new kinematics-based localization rule is proposed. It generalizes the one formulated by Cambou et al. [B. Cambou, Ph. Dubujet, F. Emeriault, F. Sidoroff, Eur. J. Mech. A/Solids 14 (1995) 225–276] in the case of an isotropic contact distribution. The results of the complete model compare well to numerical simulations results when available [C.S. Chang, C.L. Liao, Appl. Mech. Rev. 47 (1 Part 2) (1994) 197–207] (case of isotropic distribution of contacts). Finally, the interest of the fourth order fabric tensor based approach combined with the proposed localization rule is shown for different distributions of contacts by comparing its predictions to those given by a second order fabric tensor approach.     

A symmetry invariant formulation of the relationship between the elasticity tensor and the fabric tensor

The fabric tensor is employed as a quantitative stereological measure of the structural anisotropy in the pore architecture of a porous medium. Earlier work showed that the fabric tensor can be used additionally to the porosity to describe the anisotropy in the elastic constants of the porous medium. This contribution presents a reformulation of the relationship between fabric tensor and anisotropic elastic constants that is approximation free and symmetry-invariant. From specific data on the elastic constants and the fabric, the parameters in the reformulated relationship can be evaluated individually and efficiently using a simplified method that works independent of the material symmetry. The well-behavedness of the parameters and the accuracy of the method was analyzed using the Mori–Tanaka model for aligned ellipsoidal inclusions and using Buckminster Fuller’s octet-truss lattice. Application of the method to a cancellous bone data set revealed that employing the fabric tensor allowed explaining 75–90% of the total variance. An implementation of the proposed methods was made publicly available.     

On the influence of inter-particle friction and dilatancy in granular materials: a numerical analysis

Mechanical behavior of granular soils is a classic research realm but still yet not completely understood as it can be influenced by a large number of factors, including confining pressure, soil density, loading conditions, and anisotropy of soil etc. Traditionally granular materials are macroscopically regarded as continua and their particulate and discrete nature has not been thoroughly considered although many researches indicate the macro mechanical behavior closely depends on the micro-scale characteristics of particles. This paper presents a DEM (discrete element method)-based micromechanical investigation of inter-particle friction effects on the behavior of granular materials. In this study, biaxial DEM simulations are carried out under both ‘drained’ and ‘undrained’ (constant volume) conditions. The numerical experiments employ samples having similar initial isotropic fabric and density, and the same confining pressure, but with different inter-particle friction coefficient. Test results show that the inter-particle friction has a substantial effect on the stress-strain curve, peak strength and dilatancy characteristics of the granular assembly. Clearly, it is noted that apart from the inter-particle friction, the shear resistance is also contributed to the dilation and the particle packing and arrangements. The corresponding microstructure evolutions and variations in contact properties in the particulate level are also elaborated, to interpret the origin of the different macro-scale response due to variations in the inter-particle friction.     

Determining particulate–solid interphase strength using shear-induced anisotropy

The resistance to shear deformation developed by a granular material layer in contact with a topographically rough natural or manufactured solid material surface is critical to the stability of a variety of composite systems. By using discrete-element method numerical simulations, we show that evolution of fabric and contact force anisotropy at the boundary between the surface and the granular media controls shear behavior. Full mobilization of granular material strength occurs when the contact force anisotropy developed at the interface is equal to the maximum contact force anisotropy of the granular media.     

Scan line void fabric anisotropy tensors of granular media

Soil fabric anisotropy tensors are related to the statistical distribution of orientation of different microstructural vector-like entities, most common being the contact normal vectors between particles, which are extremely difficult to determine for real granular materials. On the other hand, void fabric based tensors can be determined by image based quantification methods of voids (graphical approaches), which are well defined and easy to apply to both physical and numerical experiments. A promising void fabric characterization approach is based on the scan line method. Existing scan line based definitions of void fabric anisotropy tensors are shown analytically to inherit a shortcoming, since numerous small void segments in a sample have an inordinate contribution towards unwarranted isotropy. Discrete Element Method (DEM) of analysis subsequently confirms this analytical proof. The fact that such scan line void fabric tensor definitions yield acceptable results when used in conjunction with physical image-based measurements, is shown to be attributed to the natural “cut off” of smaller void segments that occurs during such measurements. This is the motivation to propose using the existing definition of void fabric tensors, with exclusion of void segments less than a “cut off” value associated with an internal length of the granular assembly. In addition, an entirely new void fabric tensor was introduced using the squared length, instead of the length of a void segment, as the weighting factor for the definition of the scan line void fabric tensor. It was found by means of DEM analysis that both alternative definitions are void of the aforementioned shortcoming and compatible with existing image quantification methods of void fabric anisotropy.     

Dependence of elastic constants of an anisotropic porous material upon porosity and fabric

The elastic properties of an anisotropic porous material can be represented as functions of the material's solid volume fraction (or porosity) and the principal diameters of the material's fabric ellipsoid. The fabric ellipsoid is a measure of the anisotropy of the microstructure of a material. The definitions and measurement techniques for fabric ellipsoids in granular materials, foams, cancellous bone, and rocks are discussed. The principal results presented in this work are algebraic expressions for the dependence of the orthotropic elastic constants upon both solid volume fraction and the fabric ellipsoid.