Explicit time-discretisation of elastodynamics with some inelastic processes at small strains

Tomas Roubicek,Christos Panagiotopoulos,Chrysoula Tsogka

DOI: https://doi.org/10.48550/arXiv.1903.11654

2019-03-30

Abstract:The 2-step staggered (also called leap-frog) time discretisation of linear 2nd-order Hamiltonian systems (typically linear elastodynamics in a stress-velocity form) is extended for a 3-step staggered discretisation applicable for systems involving some internal variables subjected to a dissipative evolution. After spatial discretisation, a-priori estimates and convergence is proved under the usual CFL-condition. Applications to specific problems in continuum mechanics of solids at small stains are considered, in particular linearized plasticity, diffusion in poroelastic media, damage, or adhesive contact. Numerical implementation and some computational 2-dimensional simulation of waves emitted by a rupture (delamination) of an adhesive contact illustrate the abstract theory and efficiency of the explicit method.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop a time - discretization method for small - strain elastodynamics problems applicable to dissipative processes containing certain internal variables. Specifically, the author proposes a three - step staggered explicit time - discretization scheme to handle the evolution of internal variables in linear second - order Hamiltonian systems (such as linear elastodynamics in stress - velocity form). This method not only maintains the high efficiency of explicit methods but can also handle nonlinear dissipative processes.

### Main problems and objectives

1. **Limitations of existing methods**:
- Although implicit time - discretization methods can ensure energy conservation, they have high computational costs in high - dimensional problems.
- Although the simple explicit central - difference scheme is efficient, it has difficulties when introducing internal variables and nonlinear dissipative processes.

2. **Objectives of the paper**:
- Develop a new explicit staggered time - discretization method that can maintain the high efficiency of explicit methods and handle the complex evolution of internal variables.
- Prove the numerical stability and convergence of this method under specific conditions.
- Apply this method to specific continuum mechanics problems, such as plasticity, diffusion, damage, and adhesive contact.

### Solutions

The author proposes a three - step staggered time - discretization scheme. The specific steps are as follows:

1. **Calculate the generalized stress \(\Sigma^{k + 1}_{\tau h}\)**:
\[
\frac{\Sigma^{k + 1}_{\tau h}-\Sigma^k_{\tau h}}{\tau}=C E_h v^k_{\tau h}+D^k_{\tau h}
\]

2. **Calculate the internal variable \(z^{k + 1}_{\tau h}\)**:
\[
\partial\Psi\left(\frac{z^{k + 1}_{\tau h}-z^k_{\tau h}}{\tau}\right)+\Phi'_z\left(\Sigma^{k + 1}_{\tau h},\frac{z^{k + 1}_{\tau h}+z^k_{\tau h}}{2}\right)\ni0
\]

3. **Calculate the velocity \(v^{k + 1}_{\tau h}\)**:
\[
T'\frac{v^{k + 1}_{\tau h}-v^k_{\tau h}}{\tau}+E^*_h S^{k + 1}_{\tau h}=F^{k + 1}_{\tau h},\quad\text{where}\ S^{k + 1}_{\tau h}=C^*\Phi'_\Sigma(\Sigma^{k + 1}_{\tau h},z^{k + 1}_{\tau h})
\]
And update the displacement \(u^{k + 1}_{\tau h}\):
\[
u^{k + 1}_{\tau h}=u^k_{\tau h}+\tau v^{k + 1}_{\tau h}
\]

### Numerical stability and convergence

To ensure numerical stability, the author introduces an improved CFL condition:
\[
\langle [W']^{-1}\Sigma_h,\Sigma_h\rangle\geq\frac{\tau^2}{4}\langle E^*_h S_h,(T')^{-1}E^*_h S_h\rangle
\]
where \(S_h = C^*\Phi'_\Sigma(\Sigma_h,\tilde{z}_h)\).

Through strict mathematical derivations, the author proves that under the above CFL condition, this method has good numerical stability and convergence and can be applied to a variety of practical problems.

### Application examples

Finally, the paper verifies the effectiveness and efficiency of this method through several two - dimensional wave - propagation simulation experiments.