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}.

Increased OXT in the hypothalamic PVN in male mood disorder patients: A possible direct inhibition of testosterone on OXT expression

Discrete Element Modeling of Sand Behavior in a Biaxial Shear Test

A Coupled FEM/DEM Approach for Hierarchical Multiscale Modelling of Granular Media

Discrete Element Simulations of Drained Granular Material Response under Multidirectional Rotational Shear

Undrained Cyclic Behavior of Granular Materials Considering Initial Static Shear Effect: Insights from Discrete Element Modeling

A Coupled SPFEM/DEM Approach for Multiscale Modeling of Large‐deformation Geomechanical Problems

Novel Approach for Generating Homogeneous Samples for Discrete-Element-Method Studies

Discrete Element Modelling of a Soil-Rock Mixture Used in an Embankment Dam

Multiscale Simulation of the Responses of Discrete Nanostructures to Extreme Loading Conditions Based on the Material Point Method

Microscopic Mechanism Analysis of Influence of High Effective Confining Pressure on Mechanical Properties of Hydrate-Bearing Sediments

Assessing the effect of layered spatial variability on soil behavior via DEM simulation

DEM Study on the Development of the Earth Pressure of Granular Materials Subjected to Lateral Cyclic Loading

Unified discrete‐element approach applying arbitrary undrained loading paths in element testing for granular soils

Mechanical behavior of clayey methane hydrate-bearing sediment: Effects of pore size and physicochemical characteristics using a novel DEM contact model

Study on the Multiphase Fluid-Solid Interaction in Granular Materials Based on an LBM-DEM Coupled Method

DEM simulation of microscopic structure and macroscopic mechanical behavior of clay in oedometer and triaxial compression tests

Coupled Lbm-Dem Modeling Of Bingham Fluids With Suspended Particles

DEM Analyses of the Mechanical Behavior of Soil and Soil-Rock Mixture Via the 3D Direct Shear Test

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

DEM investigation on flow instability of particulate assemblies under coupling between volumetric and axial strains

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