A numerical method for the study of shear band propagation in soft rocks