Non-conforming FEM for the quasi-static contact problem

Kamana Porwal,Tanvi Wadhawan
2024-01-04
Abstract:In this article, we addressed the numerical solution of a non-linear evolutionary variational inequality, which is encountered in the investigation of quasi-static contact problems. Our study encompasses both the semi-discrete and fully-discrete schemes, where we employ the backward Euler method for time discretization and utilize the lowest order Crouzeix-Raviart nonconforming finite element method for spatial discretization. By assuming appropriate regularity conditions on the solution, we establish \emph{a priori} error analysis for these schemes, achieving the optimal convergence order for linear elements. To illustrate the numerical convergence rates, we provide numerical results on a two-dimensional test problem.
Numerical Analysis
What problem does this paper attempt to address?
This paper aims to solve the numerical solution in the quasi - static contact problem and its convergence analysis. Specifically, the author focuses on how to obtain effective approximate solutions through numerical methods when dealing with such nonlinear evolution variational inequalities, and ensure that the errors of these approximate solutions reach the optimal convergence order under certain assumptions. ### Background and Problem Description of the Paper 1. **Research Background** - Contact problems play a central role in mechanical structures and are involved in many aspects of daily life, such as the interactions between brake pads and wheels, tires and roads, pistons and skirts. - Variational inequalities are effective tools for dealing with such contact problems, especially in the contact problems of elastic materials. - Quasi - static contact problems refer to the situation where the external forces on the system change gradually over time, but the acceleration can be ignored. Such problems are particularly important in the mechanics of materials. 2. **Research Objectives** - The objective of this paper is to conduct a numerical analysis of quasi - static contact problems, especially using the Tresca friction law to model the friction process. - The research includes spatially semi - discrete and fully - discrete schemes, using the finite difference method for time discretization and the non - conforming finite element method for spatial discretization respectively. ### Mathematical Model and Method 1. **Strong Form** - It describes a series of equations that the displacement vector \( \mathbf{u}: \Omega\times[0, T]\to\mathbb{R}^2 \) needs to satisfy, including equilibrium equations, boundary conditions and initial conditions. - The contact conditions are given by the Tresca friction law and are divided into the sticking zone and the sliding zone. 2. **Variational Form** - An appropriate Hilbert space \( V \) is introduced, and the variational inequality is defined in this space. - The objective is to find the displacement vector \( \mathbf{u}(t) \) that satisfies the variational inequality within a given time interval. ### Numerical Method and Analysis 1. **Semi - discrete Scheme** - The Crouzeix - Raviart non - conforming finite element method is used for spatial discretization. - By introducing the discrete bilinear form \( a_h \) and the norm \( |||\cdot||| \), the semi - discrete variational inequality is established. - A priori error analysis is carried out, and it is proved that under appropriate regularity assumptions, the error converges at the order of \( O(h) \). 2. **Fully - discrete Scheme** - Combining time discretization (backward Euler method) and spatial discretization (Crouzeix - Raviart non - conforming finite element method), the fully - discrete variational inequality is obtained. - Through detailed error estimation, it is proved that under appropriate regularity assumptions, the error converges at the order of \( O(h + k) \), where \( h \) is the spatial step size and \( k \) is the time step size. ### Conclusion Through the numerical analysis of quasi - static contact problems, this paper successfully establishes semi - discrete and fully - discrete schemes, and under appropriate regularity assumptions, proves the error convergence of these schemes. These results provide theoretical basis and numerical method support for dealing with contact problems in practical engineering.