Abstract:Peridynamics is a nonlocal generalization of continuum mechanics theory which adresses discontinuous problems without using partial derivatives and replacing its by an integral operator. As a consequence, it finds applications in the framework of the development and evolution of fractures and damages in elastic materials.
In this paper we consider a one-dimensional nonlinear model of peridynamics and propose a suitable two-dimensional fast-convolution spectral method based on Chebyshev polynomials to solve the model. This choice allows us to gain the same accuracy both in space and time. We show the convergence of the method and perform several simulations to study the performance of the spectral scheme.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: how to effectively solve the one - dimensional nonlinear Peridynamics model. Specifically, the paper proposes a fast convolution spectral method based on Chebyshev polynomials to obtain numerical solutions with high precision simultaneously in space and time. This method not only avoids the time - integration steps required in traditional methods, but also improves the computational efficiency by using the Fast Fourier Transform algorithm, thus being able to maintain high efficiency and high precision when dealing with non - periodic problems.
### Background and Problem Description of the Paper
Peridynamics is a non - local continuum mechanics theory that can describe the development and evolution processes of cracks and damage in materials. Different from traditional partial differential equations, Peridynamics uses integral operators to replace partial derivatives, so it can handle discontinuous problems. However, existing numerical methods such as mesh - free methods and finite - difference methods have high computational costs when dealing with non - local problems, especially when the non - local area covers a large range.
### Main Contributions of the Paper
1. **Proposing a New Spectral Method**: The paper proposes a fast convolution spectral method based on Chebyshev polynomials, which can obtain numerical solutions with high precision simultaneously in space and time.
2. **Avoiding Time Integration**: By expanding the unknown function in the two - dimensional space - time domain, this method does not require explicit time - integration steps, but is achieved by solving an algebraic system.
3. **Efficient Computation**: By using the Fast Fourier Transform algorithm, this method has an advantage in computational efficiency, especially when dealing with large - scale problems.
4. **Theoretical Analysis**: The paper proves the convergence of the proposed method and verifies its performance through numerical simulations.
### Mathematical Model
The nonlinear Peridynamics model studied in the paper can be expressed as:
\[
\frac{\partial^2 u(x,t)}{\partial t^2}=\int_{B_\delta(x)} f(x - x', u(x', t)-u(x, t), t) \, dx', \quad x\in\Omega, \quad t > 0
\]
where \( u(x,t) \) represents the displacement field, and \( f \) is a pairwise force function, which is assumed to be in the form of:
\[
f(\xi,\eta)=C(\xi)H(\eta)
\]
Here, \( C(\xi) \) is a micromodulus function, and \( H(\eta) \) is a globally Lipschitz - continuous odd function. In particular, the paper selects \( H(\eta)=\eta^3 \), and this power - nonlinear form is very useful numerically because it can utilize the properties of Chebyshev transforms and convolution products.
### Construction of the Spectral Method
The paper proposes a fully spectral method, expanding the unknown function \( u(x,t) \) in space and time in the form of Chebyshev polynomials:
\[
u_N(x,t)=\sum_{j = 0}^N\sum_{k = 0}^N u_{jk}T_j(x)T_k(t)
\]
By this method, the original problem can be transformed into an algebraic system, thus avoiding explicit time - integration steps. In addition, the choice of Chebyshev polynomials makes the method no longer dependent on periodic boundary conditions.
### Convergence and Numerical Simulation
The paper proves the convergence of the proposed method and verifies its performance through numerical simulations. Numerical results show that this method has high precision and high efficiency when dealing with non - periodic problems.
### Conclusion
The paper proposes a fast convolution spectral method based on Chebyshev polynomials for solving the one - dimensional nonlinear Peridynamics model. This method obtains numerical solutions with high precision simultaneously in space and time and has high computational efficiency. Through theoretical analysis and numerical simulations, the effectiveness and superiority of the method are verified.