The virtual element method for a contact problem with wear and unilateral constraint

This paper is dedicated to the numerical solution of a mathematical model that describes frictional quasistatic contact between an elastic body and a moving foundation, with the wear effect on the contact interface of the moving foundation due to friction. The mathematical problem is a system consisting of a time-dependent quasi-variational inequality and an integral equation. The numerical method is based on the use of the virtual element method (VEM) for the spatial discretization of the variational inequality and a variable step-size left rectangle integration formula for the integral equation. The existence and uniqueness of a numerical solution are shown, and optimal order error estimates are derived for both the displacement and the wear function for the lowest order VEM. Numerical results are presented to demonstrate the efficiency of the method and to illustrate the numerical convergence orders.