Abstract:The critical-state response of granular assemblies composed of elastic spheres under generalised three-dimensional loading conditions was investigated using the discrete element method (DEM). Simulations were performed with a simplified Hertz–Mindlin contact model using a modified version of the LAMMPS code. Initially isotropic samples were subjected to three-dimensional stress paths controlled by the intermediate stress ratio, b=[(σ2′-σ3′)/\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$b=[(\sigma '_{2}-\sigma '_{3})/$$\end{document}(σ1′-σ3′)]\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$(\sigma '_{1}-\sigma '_{3})]$$\end{document}. Three types of simulation were performed: drained (with b\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$b$$\end{document}-value specified), constant volume and constant mean effective stress. In contrast to previous DEM observations, the position of the critical state line is shown to depend on b\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$b$$\end{document}. The data also show that, upon shearing, the dilatancy post-peak increases with increasing b\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$b$$\end{document}, so that at a given mean effective stress, the void ratio at the critical state increases systematically with b\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$b$$\end{document}. Four commonly-used three-dimensional failure criteria are shown to give a better match to the simulation data at the critical state than at the peak state. While the void ratio at critical state is shown to vary with b\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$b$$\end{document}, the coordination number showed no dependency on b\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$b$$\end{document}. The variation in critical state void ratios at the same p′\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$p'$$\end{document} value is apparently related to the directional fabric anisotropy which is clearly sensitive to b\documentclass[12pt]{minimal}\usepackage{amsmath}\usepackage{wasysym}\usepackage{amsfonts}\usepackage{amssymb}\usepackage{amsbsy}\usepackage{mathrsfs}\usepackage{upgreek}\setlength{\oddsidemargin}{-69pt}\begin{document}$$b$$\end{document}.

Comparison of Continuum Stresses in Granular Material Computed by Volume Average Approach and Boundary Average Approach Under Static and Quasi-Static Conditions

Numerical Simulations of Granular Material Behavior under Rotation of Principal Stresses: Micromechanical Observation and Energy Consideration

The Influence of Particles' Aspect Ratio on the Shear Behaviour of Granular Materials

The Influence of Particle-Size Distribution on Critical State Behavior of Spherical and Non-Spherical Particle Assemblies

A New Definition on Critical State of Granular Media Accounting for Fabric Anisotropy

Contact Force Distribution and Induced Anisotropy in Granular Materials under Biaxial Test

Incremental behavior of granular soils: a strain response envelope perspective

A Comparison of Two-Dimensional and Three-Dimensional Micromechanical Discrete Element Modeling of the Splitting Tests for Asphalt Mixtures

On the Influence of Inter-Particle Friction and Dilatancy in Granular Materials: a Numerical Analysis

DEM analysis of the influence of the intermediate stress ratio on the critical-state behaviour of granular materials

Comparison of Granular Material Behaviour under Drained Triaxial and Plane Strain Conditions Using 3D DEM Simulations

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

Comparison Between Priebe's Method and Elastic Analytical Solution of Unit Cell for Composite Foundations

A virtual experiment technique on the elementary behaviour of granular materials with discrete element method

Macro-micro Performances of Granular Materials Considering the Influences of Density and Stress Path under True Triaxial Conditions: A DEM Investigation

Discrete Element Simulation of Macro and Micro Mechanical Properties of Round Gravel Material under Triaxial Stress

Numerical Study on Energy Transformation in Granular Matter under Biaxial Compression

Discrete Averaging Relations for Micro to Macro Transition

Mix-Mode Elastic Finite Element Formulation for Bonded Granular Material Considering Rotation-Of Particles

DEM Simulation of Micro-Macro Mechanical Behaviour of Granular Materials

3D DEM-based particle-scale analysis of drained and undrained triaxial behaviors of granular materials