We study a class of shear-free, homogeneous but anisotropic cosmological models with imperfect matter sources in the context of f (R) gravity. We show that the anisotropic stresses are related to the electric part of the Weyl tensor in such a way that they balance each other. We also show that within the class of orthogonal f (R) models, small perturbations of shear are damped, ...