PurposeThe objective of this paper is to quantitatively assess shear band evolution by using two-dimensional discrete element method (DEM).Design/methodology/approachThe DEM model was first calibrated by retrospectively modelling existing triaxial tests. A series of DEM analyses was then conducted with the focus on the particle rotation during loading. An approach based on particle rotation was developed to precisely identify the shear band region from the surrounding. In this approach, a threshold rotation angle ω0 was defined to distinguish the potential particles inside and outside the shear band and an index g(ω0) was introduced to assess the discrepancy between the rotation response inside and outside shear band. The most distinct shear band region can be determined by the ω0 corresponding to the peak g(ω0). By using the proposed approach, the shear band development of two computational cases with different typical localised failure patterns were successfully examined by quantitatively measuring the inclination angle and thickness of shear band, as well as the microscopic quantities.FindingsThe results show that the shear band formation is stress-dependent, transiting from conjugated double shear bands to single shear band with confining stress increasing. The shear band evolution of two typical localised failure modes exhibits opposite trends with increasing strain level, both in inclination angle and thickness. Shear band featured a larger volumetric dilatancy and a lower coordination number than the surrounding. The shear band also significantly disturbs the induced anisotropy of soil.Originality/valueThis paper proposed an approach to quantitatively assess shear band evolution based on the result of two-dimensional DEM modelling.