Microscopic study of granular material behaviours under general stress paths

The granular material behaviour is determined by the local contact behaviour between particles and the spatial arrangement of particles. Investigation of particle-scale mechanism provides fundamental insights into global granular material behaviour. A multi-scale investigation has been carried out to study granular material behaviour under general stress paths using discrete element method (DEM). The commercial software Particle Flow Code in Three Dimensions (PFC3D) is employed for numerical simulations and the linear contact model is used to describe local contact behaviour. General loading paths were achieved by implementing a boundary control programme with independent control of both the magnitudes of three principal stresses and their principal directions. The intermediate principal stress ratio , where are the major, intermediate and minor principal stresses, and material anisotropy both had significant effect on granular material strength. The true triaxial simulation results indicated that the peak stress ratio was mainly contributed by the micro-scale contact force anisotropy. A smaller stress ratio was observed at greater a b value due to smaller degree of contact force anisotropy. Fabric anisotropy was another contributor to the material stress state. A lower peak stress ratio was obtained at a larger tilting major principal stress direction from the vertical deposition direction since smaller fabric anisotropy degree developed at larger . However, the material initial anisotropy had negligible effect on the critical stress ratio owing to the same contact force anisotropy and fabric anisotropy achieved. In true triaxial simulations, the intermediate strain increment rate ratio was generally larger than the stress ratio b since the particle-scale tangential force ratio was observed to be smaller than b value. The non-coaxial deformation observed in monotonic loading with various loading direction can be explained due to the non-coincidence between the principal fabric direction and the principal stress direction. And the degree of non-coaxiality decreased against shearing as the principal fabric direction approached loading direction gradually. The granular material response to rotational shear showed significant volumetric contraction and deformation non-coaxiality. The material internal structure rotated continuously along the principal stress rotation. The principal fabric direction did not exactly follow the rotation of principal stress direction. The fabric reorganisation mechanism accompanied by irrecoverable plastic deformation, leading to non-coaxial deformation behaviour. During rotational shear, the ultimate void ratio was determined by the stress ratio and b value but independent of initial void ratios. Under otherwise identical conditions, the greater internal structure anisotropy was observed at the higher stress ratio and at a greater b value, resulting in smaller ultimate void ratio (larger volumetric contraction). The general degree of deformation non-coaxiality decreased with increasing stress ratio and b value for rotational shear. The difference between the major principal stress direction and the major principal fabric direction was smaller at higher stress ratio and greater b value. It was interesting to note that the sample could fail during rotational shear, resulting in significant deviatoric strain developed in the first few cycles. The sample failed at a stress ratio , which was lower than the peak stress ratio obtained in monotonic loading but higher than the critical stress ratio . This indicated importance of considering stress rotation in geotechnical design and the material strength should be chosen based on the critical stress ratio rather than the peak value. The multi-scale investigation of granular material explains the strength characteristics from the micromechanical point of view. Observations on the fabric evolution have been made under various loading conditions. This may be useful information for the development of an advanced constitutive model