In the paper, the plastic deformation of heterogeneous materials is analyzed by direct numerical simulation based on the theory of an elastic-plastic orthotropic Cosserat continuum, with the plasticity condition taking into account both the shear and rotational mode of irreversible deformation. With the assumption of a block structure of a material with elastic blocks interacting through compliant plastic interlayers, this condition imposes constraints on the shear components of the asymmetric stress tensor, which characterize shear, and on the couple stresses, which irreversibly change the curvature characteristics of the deformed state of the continuum upon reaching critical values. The equations of translational and rotational motion together with the governing equations of the model are formulated as a variational inequality, which correctly describes both the state of elastic-plastic deformation under applied load and the state of elastic unloading. The numerical implementation of the mathematical model is performed using a parallel computing algorithm and an original software for cluster multiprocessor systems. The developed approach is applied to solve the problem of compressing a rectangular brick-patterned blocky rock mass by a rough nondeformable plate rotating with constant acceleration. The effect of the yield stress of the compliant interlayers on the stress-strain state of the rock mass in shear and bending is studied. The field of plastic energy dissipation in the rock mass is analyzed along with the fields of displacements, stresses, couple stresses, and rotation angle of structural elements. The obtained results can help to validate the hypothesis about the predominant effect of curvature on plastic strain localization at the mesolevel in microstructural materials.