Abstract:We investigate numerically a quasi-static elasticity system of Kachanov-type. To do so we propose an Euler time discretization combined with a suitable finite elements scheme (FEM) to handle the discretization is space. We use ODE-type arguments to prove the consistency of the scheme as well as its convergence rate. We rely on the computational platform FEniCS to perform the FEM discretizations in space needed to compute the model output. The simulation results show a good agreement with both the physics of the problem and with our previous qualitative mathematical analysis results obtained for precisely the same problem setting. Furthermore, our implementation recovers nicely the theoretically expected convergence rate. This is a preliminary study preparing the framework for the rigorous numerical identification of the damage process in Kachanov-type models.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to study the damage evolution process in quasi - static elastic systems through numerical simulation, especially the Kachanov - type damage model. Specifically, the author hopes to verify and understand the behavior of this model within the physically meaningful parameter range through numerical methods, and lay the foundation for the subsequent more rigorous numerical identification framework.
### Main Objectives of the Paper
1. **Numerical Simulation**: Study the characteristics of Kachanov - type damage evolution in quasi - static elastic systems through numerical simulation.
2. **Discretization Scheme**: Propose and verify a method combining Euler time discretization and finite element space discretization (FEM) to handle the numerical simulation of damage evolution.
3. **Consistency and Convergence**: Use ordinary differential equation (ODE) - type arguments to prove the consistency and convergence rate of the proposed discretization scheme.
4. **Physical Verification**: Ensure that the numerical results are consistent with the physical properties of the problem and the previous mathematical analysis results, and achieve the theoretically expected convergence rate.
5. **Framework Preparation**: Provide a preliminary research framework for the strict numerical identification of the damage process in the Kachanov - type model.
### Mathematical Model and Numerical Method
The model equations of the Kachanov - type system described in the paper are as follows:
\[
\sigma(t, x)=(1 - d(t, x))(\lambda(x)\text{tr}(\varepsilon(u(t, x)))I + 2\mu(x)\varepsilon(u(t, x)))
\]
\[
-\text{div}(\sigma(t, x)) = f(t, x)
\]
\[
d'(t, x)=(1 - d(t, x))^{-\alpha}g(t, x,\nabla_\mu u(t, x))
\]
The boundary conditions and initial conditions include:
- \(u(t, x) = 0\) on \(S\times\Gamma_0\)
- \(\sigma(t, x)\nu(t, x) = 0\) on \(S\times\Gamma_1\)
- \(\sigma(t, x)\nu(t, x)=\tau(t, x)\) on \(S\times\Gamma_2\)
- \(u(0, x)=u_0(x)\) in \(\Omega\)
- \(d(0, x)=d_0(x)\) in \(\Omega\)
### Research Contributions
1. **Verification of Numerical Method**: The effectiveness of the proposed discretization scheme has been verified through numerical simulation, especially its performance within the physically reasonable parameter range.
2. **Combination of Theory and Practice**: It shows a good consistency between the numerical results and the theoretical expectations, verifying the previous mathematical analysis results on this problem.
3. **Framework Construction**: It provides a solid foundation for the subsequent further research, especially for the numerical identification of nonlinear coupling problems.
In conclusion, this paper aims to deeply understand the behavior of the Kachanov - type damage model through numerical simulation and provide support for future numerical identification research.