The interplay between anisotropy and strain localisation in granular soils: a multiscale insight

Abstract

This paper presents a multiscale investigation on the interplay among inherent anisotropy, fabric evolution and strain localisation in granular soils, based on a hierarchical multiscale framework with rigorous coupling of the finite-element method (FEM) and discrete-element method (DEM). DEM assemblies with elongated particles are generated to simulate inherent anisotropy and are embedded to the Gauss points of the FEM mesh to derive the required constitutive relation. Specimens prepared with different bedding plane angles are subjected to biaxial shear under either smooth or rough loading platens. Key factors and physical mechanisms contributing towards the occurrence and development of strain localisation are examined. The competing evolutions of two sources of anisotropy, one related to particle orientations and the other related to contact normals, are found to underpin the development of the shear band. A single band pattern is observed under smooth boundary conditions, and its orientation relative to the bedding plane depends critically on the relative dominance between the two anisotropies. Under rough boundary conditions, the non-coaxial material response and the boundary constraint jointly lead to cross-shaped double shear bands. The multiscale simulations indicate that the DEM assemblies inside the shear band(s) undergo extensive shearing, fabric evolution and particle rotation, and may reach the critical state, while those located outside the shear band(s) experience mild loading followed by unloading. The particle-orientation-based fabric anisotropy needs significantly larger shear and dilation for mobilisation than the contact-normal based one. The asynchrony in evolution of the two fabric anisotropies can cause non-coaxial responses for initially coaxial packings, which directly triggers strain localisation.