A First-order Two-scale Analysis for Contact Problems with Small Periodic Configurations

Changqing Ye,Junzhi Cui

DOI: https://doi.org/10.1137/19M1252326

2019-08-14

Abstract:This paper is devoted to studying a type of contact problems modeled by hemivariational inequalities with small periodic coefficients appearing in PDEs, and the PDEs we considered are linear, second order and uniformly elliptic. Under the assumptions, it is proved that the original problem can be homogenized, and the solution weakly converges. We derive an $O(\epsilon^{1/2})$ estimation which is pivotal in building the computational framework. We also show that Robin problems--- a special case of contact problems, it leads to an $O(\epsilon)$ estimation in $L^2$ norm. Our computational framework is based on finite element methods, and the numerical analysis is given, together with experiments to convince the estimation.

Numerical Analysis

What problem does this paper attempt to address?

### What problems does this paper attempt to solve?

This paper mainly studies a class of contact problems with small - periodic coefficients, which are modeled by hemivariational inequalities. Specifically, the author focuses on such problems in linear, second - order and uniformly elliptic partial differential equations (PDEs). The main objectives of the paper are:

1. **Homogenization analysis**: Prove that the original problem can be homogenized and the solution converges weakly. The author derives an estimate of \(O(\epsilon^{1/2})\), which is crucial for constructing the computational framework.

2. **Special treatment of Robin problems**: For Robin problems (a special case of contact problems), the author shows that the difference between the original solution and the homogenized solution in the \(L^2\) norm is \(O(\epsilon)\).

3. **Numerical methods and experimental verification**: Establish a computational framework based on the finite - element method and verify the correctness of the above estimates through numerical analysis and experiments.

#### Main contributions

- **Homogenization results**: Prove that the original problem can be simplified by homogenization theory and give the results of weak convergence.
- **Error estimates**: Provide error estimates of \(O(\epsilon^{1/2})\) and \(O(\epsilon)\) for general contact problems and Robin problems respectively.
- **Numerical framework**: Propose a computational framework based on the finite - element method and verify the effectiveness of the theoretical results through numerical experiments.

### Formula summary

- **Homogenization coefficient**:
\[
\hat{A}_{il}=\int_Q A_{il}+A_{ij}\partial_j N_l \, dy
\]
where \(N_l(y)\) is the correction function, which satisfies the PDEs with periodic boundary conditions.

- **Error estimate**:
- For general contact problems:
\[
\|\nabla w_\epsilon\|_{0,\Omega}\leq C\epsilon^{1/2}\|u_0\|_{2,\Omega}
\]
- For Robin problems:
\[
\|u_\epsilon - u_0\|_{0,\Omega}\leq C\epsilon\|u_0\|_{2,\Omega}
\]

These formulas and conclusions provide a solid theoretical basis and effective numerical tools for understanding and solving contact problems with small - periodic structures.