In this paper numerical approximation of history-dependent hemivariational inequalities with constraint is considered, and corresponding Céa’s type inequality is derived for error estimate. For a viscoelastic contact problem with normal penetration, an optimal order error estimate is obtained for the linear element method. A numerical experiment for the contact problem is reported which provides numerical evidence of the convergence order predicted by the theoretical analysis.