On geotechnical problems, many localized deformations are observed before the ultimate state. In such cases, it is important to predict how such localization is initiated and developed. In this study, the localization of deformation is simulated numerically using a constitutive model with a non-coaxial plasticity and a double-slip mechanism. Results show that material properties such as dilatancy and boundary conditions affect the local development of localization, and that the scale effect on the bearing capacity can be quantitatively predicted. As another example, this constitutive model is applied to fault problems.