返回
Exponential time propagators for elastodynamics
DOI:10.1016/j.jmps.2024.105871.png)
摘要
En 中文
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 accurately evaluate the action of the exponential matrix on a vector. In particular, we use an apriori error estimate to predict the optimal Krylov 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-100x larger time-steps than conventional time-marching schemes such as Newmark-fl method. Computationally, it translates to a similar to 1000x and similar to 10-100x speed-up over conventional time-marching schemes for linear and nonlinear elastodynamics, respectively.
Keyword:
Dynamics
Numerical algorithms
Exponential integrators
期刊
IF:
6
论文数:
5.2K
被引数:
3.0W
机构
引用论文
Magnetic study of large Apollo samples: Possible evidence for an ancient centered dipolar field on the Moon大型阿波罗样品的磁性研究: 月球上古代中心偶极场的可能证据
Seasonal Variation of Picocyanobacterial Community Composition in Seawaters Around Dokdo, Korea
The Sea
IF0
Real time time-dependent density functional theory using higher order finite-element methods
PHYSICAL REVIEW B
IF3.7

