SUMMARYThe paper deals with numerical investigations of the behaviour of granular bodies during shearing. Shearing of a narrow layer of sand between two very rough boundaries under constant vertical pressure is numerically modelled with a "nite element method using a hypoplastic constitutive relation within a polar (Cosserat) continuum. The constitutive relation was obtained through an extension of a non-polar one by polar quantities, viz. rotations, curvatures, couple stresses using the mean grain diameter as a characteristic length. This relation can reproduce the essential features of granular bodies during shear localization. The material constants can be easily determined from element test results and can be estimated from granulometric properties. The attention is laid on the in#uence of the initial void ratio, pressure level, mean grain diameter and grain roughness on the thickness of shear zones. The results of shearing are also compared to solutions without the polar extensions. The FE-calculations demonstrate that polar e!ects manifested by the appearance of grain rotations and couple stresses are signi"cant in the shear zone, and its thickness is sensitive to the initial void ratio, mean grain diameter and layer height. The e!ect of the pressure level is rather low within the considered range.