The influence of particle geometry and the intermediate stress ratio on the shear behavior of granular materials

The behavior of granular materials is very complex in nature and depends on particle shape, stress path, fabric, density, particle size distribution, amongst others. This paper presents a study of the effect of particle geometry (aspect ratio) on the mechanical behaviour of granular materials using the discrete element method (DEM). This study discusses 3D DEM simulations of conventional triaxial and true triaxial tests. The numerical experiments employ samples with different particle aspect ratios and a unique particle size distribution (PSD). Test results show that both particle aspect ratio (AR) and intermediate stress ratio \((b=({\upsigma }_{2}'-{\upsigma }_{3}')/({\upsigma }_{1}'-{\upsigma }_{3}'))\) affect the macro- and micro-scale responses. At the macro-scale, the shear strength decreases with an increase in both aspect ratio and intermediate stress ratio b values. At the micro-scale level, the fabric evolution is also affected by both AR and b. The results from DEM analyses qualitatively agree with available experimental data. The critical state behaviour and failure states are also discussed. It is observed that the position of the critical state loci in the compression \((e-p')\) space is only slightly affected by aspect ratio (AR) while the critical stress ratio is dependent on both AR and b. It is also demonstrated that the influence of the aspect ratio and the intermediate stress can be captured by micro-scale fabric evolutions that can be well understood within the framework of existing critical state theories. It is also found that for a given stress path, a unique critical state fabric norm is dependent on the particle shape but is independent of critical state void ratio.     

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.     

Micromechanical modelling of oval particulates subjected to bi-axial compression

One of the questions that still remain unanswered among researchers dealing with granular materials is how far the particle shape affects the micro-macroscopic features of granular assemblies under mechanical loading. The latest advances made with particle instrumentation allow us to capture realistic particle shapes and size distribution of powders to a fair degree of accuracy at different length scales. Industrial applications often require information on the micromechanical behaviour of granular assemblies having different particle shapes and varying surface characteristics, which still remains largely unanswered. Traditionally, simulations based on discrete element method (DEM) idealise the shape of individual particles as either circular or spherical. In the present investigation, we analyse the influence of particle shape on the shear deformation characteristics of two dimensional granular assemblies using DEM. We prepared the assemblies having nearly an identical initial packing fraction (dense), but with different basic shapes of the individual particles: (a) oval and (b) circular for comparison purposes. The granular assemblies were subjected to bi-axial compression test. We present the evolution of macroscopic strength parameters and microscopic structural/topological parameters during mechanical loading. We show that the micromechanical properties of granular systems are significantly influenced by the shape of the individual particles constituting the granular assemblies.     

DEM simulation of the shear behaviour of breakable granular materials with various angularities

Particle shape is an important factor that affects particle breakage and the mechanical behaviour of granular materials. This report explored the effect of angularity on the mechanical behaviour of breakable granular materials under triaxial tests. Various angular particles are generated using the quasi-spherical polyhedron method. The angularity α is defined as the mean exterior angle of touching faces in a particle model. A breakable particle is constructed as an aggregate composed of coplanar and glued Voronoi polyhedra. After being prepared under the densest conditions, all assemblies were subjected to triaxial compression until a critical state was reached. The macroscopic characteristics, including the shear strength and dilatancy response, were investigated. Then, particle breakage characteristics, including the extent of particle breakage, breakage pattern and correlation between the particle breakage and energy input, were evaluated. Furthermore, the microscopic characteristics, including the contact force and fabric anisotropy, were examined to probe the microscopic origins of the shear strength. As α increases, the peak shear strength increases first and then remains constant, while the critical shear strength generally increases. Assemblies with larger angularity tend to cause more serious particle breakage. The relative breakage is linearly correlated with α under shear loading. Compared with unbreakable particles, the peak shear strength and the critical volumetric strain decline, and the degree of decline linearly increases with increasing α.     

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.     

Critical state behaviour of granular materials from isotropic and rebounded paths: DEM simulations

This paper presents the results of numerical simulations using the three-dimensional discrete element method (DEM) on the critical state behaviour of isotropically compressed and rebounded assemblies of granular materials. Drained and undrained (constant volume) numerical simulations were carried out. From these numerical simulations of drained and undrained tests, it has been shown that the steady state is same as the critical state. Critical state for both isotropically compressed and rebounded assemblies form unique curved line that can be approximated by a bilinear line as proposed by Been et al. [Géotechnique 41(3): 365–381, 1991]. Further more, evolution of the internal variables such as average coordination number and induced anisotropy coefficients during shear deformation has been studied.     

Triaxial compression tests of QH-E lunar soil simulant under constant mean principal stress path using discrete element method simulations

In this paper, discrete element method (DEM) simulations are applied to investigate the triaxial compression tests of QH-E lunar soil simulant (developed by Tsinghua University, China) under constant mean principal (P) stress path. The P stress path is achieved by controlling the loading speed and direction of the axial stress and confining stress in the DEM simulations, the strain softening and dilatancy characteristics of QH-E lunar soil simulant both at low and conventional P stress levels are discussed. The results show that the deviatoric stress–strain curve is divided into hardening, softening descending, and residual strain stages, and the shear strength and residual strength increase with the increase of P stress, which are similar to the conventional triaxial compression tests. The volumetric strain versus shear strain curve shows a good linear relationship both at the linear shear dilatancy and residual shear dilatancy stages, but the slope of the linear dilatancy stage is larger than that of the residual dilatancy stage. Furthermore, it is found that the variations of shear dilatancy characteristic parameters with regard to P stress also show a good linear relationship, and the values of these parameters increase with the increase of P stress.     

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.     

Numerical analysis of compression mechanics of highly polydisperse granular mixtures with different PSD-s

The deformation of three-dimensional highly polydisperse packings of frictional spheres with continuous unimodal and discrete uniform particle size distributions (PSDs) under uniaxial compression was investigated, using the discrete element method. The stress transmission and elasticity parameters, i.e., the effective elastic modulus and the Poisson’s ratio of the granular assemblies were examined; these parameters are important in many engineering disciplines. The influence of the shape of the PSD on the porosity and the transmission of pressures in samples was observed; however, the shape of the PSD did not affect the stiffness and the Poisson’s ratio of polydisperse granular packings. The results revealed that, in some cases, knowledge of the grain size distribution is not a critical issue for modeling granular packings composed of non-uniformly sized spheres.     

Influence of particle shape on the shear behavior of superellipsoids by discrete element method in 3D

This paper aims to study the influence of particle shape on the shear strength of superellipsoidal particles by Discrete Element Method (DEM) simulations of triaxial tests in 3D. A total of forty-nine types of equiaxed superellipsoidal particles from three evolution paths have been created. The definition of effective porosity has been proposed. Our findings show that both the particle sphericity and roundness affect the shear strength of the superellipsoidal particle system. Under the mutual impact of initial porosity and particle shape, the simulation results of shear strength and volumetric strains present a trend of initially decreasing and subsequently increasing. The microstructure evolution of superellipsoidal particles during the shearing process is observed microscopically. The anisotropy of fabric reveals the mechanism of effective porosity and sphericity influencing the macroscopic shear strength at the particle scale.     

Stress-chain based micromechanics of sand with grain shape effect

The mechanical behaviors of granular media are controlled by grain properties and microstructure. The primary property of granular media is denoted by its grain shape, grain size distribution, stiffness, and interparticle friction. The grain shape itself is of particular importance. Microstructures are formed in the connection paths of contact points between grains. In this paper, the deformation of granular materials with different grain shapes was simulated using two-dimensional DEM under different stress-levels and densities. After analyzing the results, the authors investigated fabric changes. The evolution rule of stress-induced anisotropy and its limitation as well as the existence of a critical state of fabric are revealed.     

A hybrid approach for modeling of breakable granular materials using combined finite-discrete element method

It is well known that particle breakage plays a critical role in the mechanical behavior of granular materials and has been a topic subject to intensive studies. This paper presents a three dimensional fracture model in the context of combined finite-discrete element method (FDEM) to simulate the breakage of irregular shaped granular materials, e.g., sands, gravels, and rockfills. In this method, each particle is discretized into a finite element mesh. The potential fracture paths are represented by pre-inserted non-thickness cohesive interface elements with a progressive damage model. The Mohr–Coulomb model with tension cut-off is employed as the damage initiation criterion to rupture the predominant failure mode at the particle scale. The particle breakage modeling using combined FDEM is validated by the qualitative agreement between the results of simulated single particle crushing tests and those obtained from laboratory tests and prior DEM simulations. A comprehensive numerical triaxial tests are carried out on both the unbreakable and breakable particle assemblies with varied confining pressure and particle crushability. The simulated stress–strain–dilation responses of breakable granular assembly are qualitatively in good agreement with the experimental observations. The effects of particle breakage on the compressibility, shear strength, volumetric response of the fairly dense breakable granular assembly are thoroughly investigated through a variety of mechanism demonstrations and micromechanical analysis. This paper also reports the energy input and dissipation behavior and its relation to the mechanical response.     

Influence of particle shape on sheared dense granular media

We study by means of molecular dynamics simulations of periodic shear cells, the influence of particle shape on the global mechanical behavior of dense granular media. At large shear deformation samples with elongated particles, independent of their initial orientation, reach the same stationary value for both shear force and void ratio. At the micro-mechanical level the stress, the fabric and the inertia tensors of the particles are used to study the evolution of the media. In the case of isotropic particles the direction of the principal axis of the fabric tensor is aligned with the one of the principal stress, while for elongated particles the fabric orientation is strongly dependent on the orientation of the particles. The shear band width is shown to depend on the particle shape due to the tendency of elongated particles to preferential orientations and less rotation.     

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.     

Stress-induced anisotropy in sand under cyclic loading

The anisotropy of a granular material’s structure will influence its response to applied loads and deformations. Anisotropy can be either inherent (e.g. due to depositional process) or induced as a consequence of the applied stresses or strains. Discrete element simulations allow the interactions between individual particles to be explicitly simulated and the fabric can be quantified using a fabric tensor. The eigenvalues of this fabric tensor then give a measure of the anisotropy of the fabric. The coordination number is a particle scale scalar measure of the packing density of the material. The current study examines the evolution of the fabric of a granular material subject to cyclic loading, using two-dimensional discrete element method (DEM) simulations. Isotropic consolidation modifies and reduces the inherent anisotropy, but anisotropic consolidation can accentuate anisotropy. The ratio of the normal to shear spring stiffness at the particle contacts in the DEM model affects the evolution of anisotropy. Higher ratios reduce the degree of anisotropy induced by anisotropic consolidation. The anisotropy induced by cyclic loading depends on the amplitude of the loading cycles and the initial anisotropy.     

Exploring the critical state properties and major principal stress rotation of sand in direct shear test using the distinct element method

This paper provides insights into the critical state properties and major principal stress rotation of sands in direct shear tests using the distinct element method (DEM), in which a three dimensional contact model considering rolling and twisting resistances proposed by the authors was implemented. Firstly, the DEM was used to simulate a series of direct shear tests. Then the macroscopic characteristics of the DEM samples and the critical state properties within the shear band were analyzed, including the major principal stress rotation and the stress level dependency of peak strength. Finally the effect of the grain angularity was investigated. The rolling and twisting resistances were separately considered to check their individual effects. The results demonstrate the strong capability of DEM with the 3D contact model to capture the mechanical behavior of sands in direct shear tests. After a global shear strain of 15%, the shear zone reached the critical state while the whole sample did not. The major principal stress rotation angle evolved in a similar trend to that of the shear stress, and the stress level dependency of the peak strength of dense sand can be attributed to the fact that the vertical stress affects the rotation of the major principal stress. Higher shear strength, more obvious dilatancy and upward move of critical state line on the \(e\hbox {-lg}p\) plane were obtained when larger values of shape parameter (leading to greater rolling and twisting resistances) were used. And twisting resistance is much less important relatively in the presence of rolling resistance.     

Segregation due to particle shape of a granular mixture in a slowly rotating tumbler

Using DEM particle simulations we consider segregation of a binary granular particle mixture in a slowly rotating cylindrical tumbler where the particles differ only in their shape—spherical versus more cubical particles. We find that the more cubical particles segregate to the inner core of the particle bed while the spherical particles segregate to the curved walls of the tumbler. The main mechanism for this segregation is different energy dissipation rates for the different particle shape types when avalanching down along the free surface. The cubical particles, due to their sharper corners, dissipate energy much faster than the spherical particles. This results in spherical particles reaching the bottom end of the sloped, free surface which are then transported around the cylinder adjacent to the cylinder wall, as rigid body motion. In contrast to size or density segregation, the segregation due to shape is much weaker and takes longer to reach its equilibrium or steady state. In addition, the segregation occurs along the top surface rather than through the top surface (as occurs for size and density segregation). In general, in situations where two particles differ in their ease of flow (viz flowability) the more rapidly flowing particle will segregate to the base of the free surface (which in the case of the tumbler results in spherical particles near the periphery) and the more slowly flowing particle will segregate underneath.     

DEM analysis of the size effects on the behavior of crushable granular materials

Four sets of individual-particle crushing tests were carried out on sandstone grains of different size with geometric similarity. The tensile strength was analyzed using Weibull statistics, and the size-hardening law was obtained. The experimental data also validated that the Weibull modulus is independent of the grain size. Considering both the shear and tensile fracture modes of the particle, the Mohr–Coulomb model with a tension cut-off was employed as the fracture criterion of a single particle. When the particle stresses satisfied the fracture criterion, three new fragments modeled by the ‘clump’ were generated to replace the broken particle. Nine spheres with four different sizes were released from the clump and allowed to continue crushing if the fragment stresses fulfilled the criterion again. Two polydisperse assemblies with different particle sizes but same initial fabrics were prepared. DEM simulations of triaxial shear tests with different grain sizes were carried out on the crushable granular material with varied confining pressures. The simulated stress–strain–dilation responses were in agreement with the experimental observations. The macro–micro responses of the two samples, including the stress–strain–dilation behavior, the particle crushing, and the normal contact force distribution, were discussed in detail. The cause of the size effect on the shear strength and deformation was thoroughly investigated through a variety of mechanism demonstrations and micromechanical analysis.