Summary
Within the framework of the displacement‐based virtual element method (VEM), namely, for plane elasticity, an important topic is the development of optimal techniques for the evaluation of the stress field. In fact, in the classical VEM formulation, the same projection operator used to approximate the strain field (and then evaluate the stiffness matrix) is employed to recover, via constitutive law, the stress field. Considering a first‐order formulation, strains are locally mapped onto constant functions, and stresses are piecewise constant. However, the virtual displacements might engender more complex strain fields for polygons, which are not triangles. This leads to an undesirable loss of information with respect to the underlying virtual stress field. The recovery by compatibility in patches, originally proposed for finite element schemes, is here extended to VEM, aiming at mitigating such an effect. Stresses are recovered by minimizing the complementary energy of patches of elements over an assumed set of equilibrated stress modes. The procedure is simple, efficient, and can be readily implemented in existing codes. Numerical tests confirm the good performance of the proposed technique in terms of accuracy and indicate an increase of convergence rate with respect to the classical approach in many cases.