Numerical analysis of variational-hemivariational inequalities with applications in contact mechanics

Abstract: Variational-hemivariational inequalities are an important mathematical framework for nonsmooth problems. The framework can be used to study application problems from physical sciences and engineering that involve non-smooth and even set-valued relations, monotone or non-monotone, among physical quantities. Since no analytic solution formulas are expected for variational-hemivariational inequalities from applications, numerical methods are needed to solve the problems. This paper focuses on numerical analysis of variational-hemivariational inequalities, reporting new results as well as surveying some recent published results in the area. A general convergence result is presented for Galerkin solutions of the inequalities under minimal solution regularity conditions available from the well-posedness theory, and C\'{e}a's inequalities are derived for error estimation of numerical solutions. The finite element method and the virtual element method are taken as examples of numerical methods, optimal order error estimates for the linear element solutions are derived when the methods are applied to solve three representative contact problems under certain solution regularity assumptions. Numerical results are presented to show the performance of both the finite element method and the virtual element method, including numerical convergence orders of the numerical solutions that match the theoretical predictions.