Two-Level preconditioning method for solving saddle point systems in contact computation

Xiaoyu Duan,Hengbin An

2024-09-24

Abstract:In contact mechanics computation, the constraint conditions on the contact surfaces are typically enforced by the Lagrange multiplier method, resulting in a saddle point system. Given that the saddle point matrix is indefinite, solving these systems presents significant challenges. For a two-dimensional tied contact problem, an efficient two-level preconditioning method is developed. This method utilizes physical quantities for coarsening, introducing two types of interpolation operators and corresponding smoothing algorithms. Additionally, the constructed coarse grid operator exhibits symmetry and positive definiteness, adequately reflecting the contact constraints. Numerical results show the effectiveness of the method.

Numerical Analysis

What problem does this paper attempt to address?

This paper attempts to solve the difficult problem of solving saddle - point systems in contact mechanics calculations due to the use of the Lagrange multiplier method to impose contact surface constraint conditions. Specifically:

1. **Problem Background**:
- In contact mechanics calculations, the Lagrange multiplier method is usually used to enforce the constraint conditions on the contact surface, which results in a saddle - point system.
- The saddle - point matrix is indefinite, so there are significant challenges in solving these systems.

2. **Limitations of Existing Methods**:
- Direct methods are impractical for large - scale problems because they require a large amount of computational resources and memory.
- Although iterative methods are effective, for saddle - point systems, the classical Algebraic Multi - Grid (AMG) method often exhibits slow convergence or failure to converge due to the indefiniteness and coupling characteristics of the saddle - point matrix.

3. **The Method Proposed in the Paper**:
- The paper develops an efficient two - level preconditioning method, which is especially suitable for two - dimensional binding contact problems.
- This method uses physical quantities for coarsening and introduces two interpolation operators and corresponding smoothing algorithms.
- The constructed coarse - grid operator has symmetric positive definiteness and can fully reflect the contact constraint conditions.

4. **Advantages of the Method**:
- Numerical results show the effectiveness of this method.
- Compared with commonly used preconditioning methods, the two - level method shows advantages in various contact models and greatly promotes the solution process of saddle - point equations.

### Formula Summary

- **Matrix Form of Saddle - Point System**:
\[
\begin{bmatrix}
K & G^T \\
G & 0
\end{bmatrix}
\begin{bmatrix}
d \\
\lambda
\end{bmatrix}
=
\begin{bmatrix}
f \\
0
\end{bmatrix}
\]
where \( K \) is the stiffness matrix, \( G \) is the matrix representing the contact constraint, \( d \) is the displacement vector, \( \lambda \) is the Lagrange multiplier, and \( f \) is the external force vector.

- **Ideal Interpolation Operator**:
\[
\hat{P} =
\begin{bmatrix}
I_{NN} & 0 & 0 \\
0 & I_{MM} & 0 \\
0 & P & -D^{-\top}K_{SN} \\
0 & 0 & -D^{-\top}K_{SS}P
\end{bmatrix}
\]

- **Coarse - Grid Operator**:
\[
\hat{A}_H = \hat{R} A \hat{P} = S
\]
where \( S \) is the Schur complement matrix of the original saddle - point system.

Through this method, the paper effectively solves the difficult problem of solving saddle - point systems in contact mechanics calculations and shows its superior performance in numerical experiments.