An improved point-to-surface contact algorithm with penalty method for peridynamics

Haoran Zhang,Lisheng Liu,Xin Lai,Jun Li

2024-08-13

Abstract:It is significantly challenging to obtain accurate contact forces in peridynamics (PD) simulations due to the difficulty of surface particles identification, particularly for complex geometries. Here, an improved point-to-surface contact model is proposed for PD with high accuracy. First, the outer surface is identified using the eigenvalue method and then we construct a Verlet list to identify potential contact particle pairs efficiently. Subsequently, a point-to-surface contact search algorithm is utilized to determine precise contact locations with the penalty function method calculating the contact force. Finally, the accuracy of this point-to-surface contact model is validated through several representative contact examples. The results demonstrate that the point-to-surface contact model model can predict contact forces and deformations with high accuracy, aligning well with the classical Hertz contact theory solutions. This work presents a contact model for PD that automatically recognizes external surface particles and accurately calculates the contact force, which provides guidance for the study of multi-body contact as well as complex contact situations.

Computational Engineering, Finance, and Science

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is the challenge of accurately obtaining contact forces in Peridynamics (PD) simulations based on non - local theory, especially when dealing with complex geometries. Specifically, the paper proposes an improved point - to - surface contact algorithm, which solves the following key problems by identifying outer - surface particles and calculating contact forces using the penalty function method:

1. **Surface Particle Identification**: In PD simulations, identifying the outer - surface particles of an object is a difficult problem, especially when dealing with complex geometries. Traditional PD methods are difficult to pre - determine the contact surface, resulting in an inability to accurately predict contact forces.

2. **Contact Force Calculation**: In order to ensure the accuracy of contact forces, an efficient method is required to identify potential contact particle pairs and accurately calculate the contact forces between them. The point - to - surface contact algorithm proposed in the paper combines the Verlet adjacency list and the penalty function method to improve the accuracy of contact force calculations.

3. **Comparison with Classical Hertzian Theory**: The paper verifies the accuracy of the proposed model through multiple representative contact examples, and the results show that the model can well predict contact forces and deformations and is highly consistent with the solutions of the classical Hertzian contact theory.

### Main Contributions of the Paper

- Proposes a new point - to - surface contact algorithm applicable to PD simulations, which can automatically identify outer - surface particles and accurately calculate contact forces.
- Uses the eigenvalue method to identify outer - surface particles, solving the problem of difficult pre - determination of the contact surface in traditional PD methods.
- Improves the search efficiency of potential contact particle pairs by constructing the Verlet adjacency list.
- Uses the penalty function method to calculate contact forces, ensuring high - precision contact force calculations.
- Verifies the effectiveness of the model through multiple contact examples, and the results are highly consistent with the classical Hertzian contact theory.

### Formula Summary

- **Eigenvalue Method for Identifying Outer - Surface Particles**:
\[
B(x_i) = \left( \sum \nabla W_{ij} \otimes \mathbf{r}_{ij} \Delta V_j \right)^{-1}
\]
\[
\mathbf{r}_{ij} = x_i - x_j, \quad \Delta V_j \text{ is the volume of the } j \text{ - th nearest - neighbor particle}
\]
\[
W_{ij} = W\left(\frac{x_i - x_j}{h}\right)
\]

- **Contact Force Calculation**:
\[
F_N = K \delta_p^n
\]
where $\delta_p$ represents the relative penetration depth, $K$ is the contact stiffness parameter determined by material and geometric properties, and $n$ is the nonlinear power exponent.

- **Contact Stiffness**:
\[
K = \frac{4}{3} \left( \frac{1}{m_i} + \frac{1}{m_j} \right) R_i
\]
\[
m_l = \frac{1 - v_l^2}{E_l}, \quad (l = i, j)
\]
where $v_l$ and $E_l$ are the Poisson's ratio and elastic modulus of the two contacting objects respectively.

Through these improvements, the paper provides an effective and accurate method for dealing with contact problems in PD, providing important guidance for subsequent research.