Exponential time propagators for elastodynamics

Paavai Pari,Bikash Kanungo,Vikram Gavini
2024-05-09
Abstract:We propose a computationally efficient and systematically convergent approach for elastodynamics simulations. We recast the second-order dynamical equation of elastodynamics into an equivalent first-order system of coupled equations, so as to express the solution in the form of a Magnus expansion. With any spatial discretization, it entails computing the exponential of a matrix acting upon a vector. We employ an adaptive Krylov subspace approach to inexpensively and and accurately evaluate the action of the exponential matrix on a vector. In particular, we use an apriori error estimate to predict the optimal Kyrlov subspace size required for each time-step size. We show that the Magnus expansion truncated after its first term provides quadratic and superquadratic convergence in the time-step for nonlinear and linear elastodynamics, respectively. We demonstrate the accuracy and efficiency of the proposed method for one linear (linear cantilever beam) and three nonlinear (nonlinear cantilever beam, soft tissue elastomer, and hyperelastic rubber) benchmark systems. For a desired accuracy in energy, displacement, and velocity, our method allows for $10-100\times$ larger time-steps than conventional time-marching schemes such as Newmark-$\beta$ method. Computationally, it translates to a $\sim$$1000\times$ and $\sim$$10-100\times$ speed-up over conventional time-marching schemes for linear and nonlinear elastodynamics, respectively.
Numerical Analysis,Computational Engineering, Finance, and Science
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to achieve efficient and accurate long - time simulations in elastodynamic simulations. Specifically, the author proposes a method based on exponential - time propagators for fast and accurate long - time simulations of linear and nonlinear elastodynamic systems. ### Main Problems 1. **Limitations of Existing Methods**: - **Explicit Methods**: Although easy to implement and with low computational cost, the time step is very small due to stability limitations. - **Implicit Methods**: Although stable at larger time steps, the computational cost per step is high because it is necessary to solve linear equations. - **Highly Nonlinear Systems**: Existing time - integration methods perform poorly when dealing with highly nonlinear systems and are prone to becoming unstable during long - time calculations. 2. **Objectives**: - **Improve Accuracy**: Maintain high precision at larger time steps. - **Improve Efficiency**: Reduce computational time, especially in large - scale systems. - **Maintain System Characteristics**: Such as energy conservation and momentum conservation to ensure long - time numerical stability. ### Solutions The author proposes a method based on exponential - time propagators. The main steps are as follows: 1. **Problem Reformulation**: - Convert second - order dynamic equations into first - order coupled equations so that the solution can be expressed using the Magnus expansion. - Use the adaptive Krylov subspace method to efficiently evaluate the action of the exponential matrix on vectors. 2. **Construction of Exponential Propagators**: - Use the Magnus expansion to construct exponential propagators for nonlinear dynamic systems. - Achieve second - order convergence by truncating the first term of the Magnus expansion. 3. **Numerical Implementation**: - **Integral Evaluation**: Use the midpoint quadrature rule to approximate the integral and maintain second - order accuracy. - **Predictor - Corrector Method**: Used to evaluate the H - matrix at future time points. - **Krylov Subspace Method**: Efficiently evaluate the action of the exponential matrix on vectors and reduce computational cost. 4. **Performance Verification**: - Verify the accuracy and efficiency of the method through linear and nonlinear benchmark systems (such as linear cantilever beams, soft - tissue elastomers, and hyperelastic rubbers). - The results show that compared with the traditional Newmark - β method, this method allows a time step that is 10 - 100 times larger under the same accuracy requirements, and the computational speed is increased by 1000 times to 100 times. ### Key Formulas - **Magnus Expansion**: \[ \Omega(t, t_0) = \int_{t_0}^{t} A(\tau) d\tau - \frac{1}{2} \int_{t_0}^{t} \left[ \int_{t_0}^{\tau} A(\sigma) d\sigma, A(\tau) \right] d\tau+\cdots \] where \([X, Y]=XY - YX\) is the commutator. - **Exponential Propagator**: \[ \bar{w}(t)=\exp(\Omega(t, t_0)) \bar{w}(t_0) \] - **Krylov Subspace Method**: \[ \exp(A)w\approx \|w\| V_m \exp(H_m) e_1 \] where \(V_m\) and \(H_m\) are generated by the Lanczos or Arnoldi methods. ### Summary This paper proposes a new method based on exponential - time propagators, aiming to solve the limitations of existing time - integration methods when dealing with highly nonlinear elastodynamic systems. Through theoretical analysis and numerical experiments, it is proved that this method significantly improves computational efficiency while maintaining high accuracy.