Numerical analysis of an evolutionary variational–hemivariational inequality with application in contact mechanics☆