Micromechanical constitutive equations are developed which allow for the broad range of interparticle interactions observed in a real deforming granular assembly: microslip contact, gross slip contact, loss of contact and an evolution in these modes of contact as the deformation proceeds. This was accomplished through a synergetic use of contact laws, which account for interparticle resistance to both sliding and rolling, together with strain-dependent anisotropies in contacts and the normal contact force. By applying the constitutive model to the bi-axial test it is demonstrated that the model can correctly predict the evolution of various anisotropies as well as the formation of a distinct shear band. Moreover, the predicted shear-band properties (e.g. thickness, prolonged localisation, void ratio) are an even better fit with experimental observations than were previously found by use of previously developed micromechanical models.