Abstract
In compression tests of geomaterials, many patterns of localized deformation are observed depending on initial imperfections and boundary conditions. In this paper, especially the developments of such localized deformation are traced numerically using a constitutive model with non-coaxial plasticity and double sliding effect. As a typical result, the final patterns of localization with initial imperfections lose their symmetry and the body shows global softening even though the local constitutive law holds hardening. Furthermore many phenomena like mode-switching are also observed along their loading paths.