Abstract As an extension of the finite element method, virtual method (VEM) can handle very general polygonal meshes, making it suitable for non-matching meshes. In (Wriggers et al. in Comput. Mech. 58:1039–1050, 2016), lowest-order was applied to solve contact problem two elastic bodies on The numerical experiments showed robustness and accuracy scheme. this paper, we establish a priori error ...