Numerical analysis of history-dependent quasivariational inequalities with applications in contact mechanics

¹ Department of Mathematics, University of Wisconsin Oshkosh, Oshkosh, WI 54901, USA
kazmis@uwosh.edu
² Laboratoire de Mathématiques et Physique, University of Perpignan, 52 Avenue Paul Alduy, 66860 Perpignan, France
barboteu@univ-perp.fr; sofonea@univ-perp.fr
³ Department of Mathematics, University of Iowa, Iowa City, IA 52242, USA
weimin-han@uiowa.edu

Received: 25 July 2012
Revised: 10 May 2013

Abstract

A new class of history-dependent quasivariational inequalities was recently studied in [M. Sofonea and A. Matei, History-dependent quasivariational inequalities arising in contact mechanics. Eur. J. Appl. Math. 22 (2011) 471–491]. Existence, uniqueness and regularity results were proved and used in the study of several mathematical models which describe the contact between a deformable body and an obstacle. The aim of this paper is to provide numerical analysis of the quasivariational inequalities introduced in the aforementioned paper. To this end we introduce temporally semi-discrete and fully discrete schemes for the numerical approximation of the inequalities, show their unique solvability, and derive error estimates. We then apply these results to a quasistatic frictional contact problem in which the material’s behavior is modeled with a viscoelastic constitutive law, the contact is bilateral, and friction is described with a slip-rate version of Coulomb’s law. We discuss implementation of the corresponding fully-discrete scheme and present numerical simulation results on a two-dimensional example.