Accurate estimations of pre-failure deformations and post-failure responses of geostructures require that the simulation tool possesses at least three main ingredients: 1) a constitutive model that is able to describe the macroscopic stress-strain-strength behavior of soils subjected to complex stress/strain paths over a wide range of confining pressures and densities, 2) an embedded length scale that accounts for the intricate physical phenomena that occur at the grain size scale in the soil, and 3) a computational platform that allows the analysis to be carried out beyond the development of an initially "contained" failure zone in the soil. In this paper, a two-scale micropolar plasticity model will be used to incorporate all these ingredients. The model is implemented in a finite element platform that is based on the mechanics of micropolar continua. Appropriate finite elements are developed to couple displacement, micro-rotations, and pore-water pressure in form of $u_n-{phi}_m$ and $u_n-p_m-{phi}_m$ (n > m) elements for analysis of dry and saturated soils. Performance of the model is assessed in a biaxial compression test on a slightly heterogeneous specimen of sand. The role of micropolar component of the model on capturing the post-failure response of the soil is demonstrated.