A high-order implicit time integration method for linear and nonlinear dynamics with efficient computation of accelerations

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: O'Shea, Daniel, Zhang, Xiaoran, Mohammadian, Shayan, Song, Chongmin
Natura: Preprint
Pubblicazione: 2024
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866909320935899136
author O'Shea, Daniel
Zhang, Xiaoran
Mohammadian, Shayan
Song, Chongmin
author_facet O'Shea, Daniel
Zhang, Xiaoran
Mohammadian, Shayan
Song, Chongmin
contents An algorithm for a family of self-starting high-order implicit time integration schemes with controllable numerical dissipation is proposed for both linear and nonlinear transient problems. This work builds on the previous works of the authors on elastodynamics by presenting a new algorithm that eliminates the need for factorization of the mass matrix providing benefit for the solution of nonlinear problems. The improved algorithm directly obtains the acceleration at the same order of accuracy of the displacement and velocity using vector operations (without additional equation solutions). The nonlinearity is handled by numerical integration within a time step to achieve the desired order of accuracy. The new algorithm fully retains the desirable features of the previous works: 1. The order of accuracy is not affected by the presence of external forces and physical damping; 2. numerical dissipation in the algorithm is controlled by a user-specified parameter, leading to schemes ranging from perfectly nondissipative A-stable to L-stable; 3. The effective stiffness matrix is a linear combination of the mass, damping, and stiffness matrices as in the trapezoidal rule. The proposed algorithm is shown to replicate the numerical results demonstrated on linear problems in previous works. Additional numerical examples of linear and nonlinear vibration and wave propagation are presented herein. Notably, the proposed algorithms show the same convergence rates for nonlinear problems as linear problems, and very high accuracy. Second-order time integration methods commonly used in commercial software produce significantly polluted acceleration responses for a common class of wave propagation problems. The high-order time integration schemes presented here perform noticably better at suppressing spurious high-frequency oscillations and producing reliable and useable acceleration responses.
format Preprint
id arxiv_https___arxiv_org_abs_2409_13397
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle A high-order implicit time integration method for linear and nonlinear dynamics with efficient computation of accelerations
O'Shea, Daniel
Zhang, Xiaoran
Mohammadian, Shayan
Song, Chongmin
Numerical Analysis
65M22
J.2; G.1.3; G.1.8
An algorithm for a family of self-starting high-order implicit time integration schemes with controllable numerical dissipation is proposed for both linear and nonlinear transient problems. This work builds on the previous works of the authors on elastodynamics by presenting a new algorithm that eliminates the need for factorization of the mass matrix providing benefit for the solution of nonlinear problems. The improved algorithm directly obtains the acceleration at the same order of accuracy of the displacement and velocity using vector operations (without additional equation solutions). The nonlinearity is handled by numerical integration within a time step to achieve the desired order of accuracy. The new algorithm fully retains the desirable features of the previous works: 1. The order of accuracy is not affected by the presence of external forces and physical damping; 2. numerical dissipation in the algorithm is controlled by a user-specified parameter, leading to schemes ranging from perfectly nondissipative A-stable to L-stable; 3. The effective stiffness matrix is a linear combination of the mass, damping, and stiffness matrices as in the trapezoidal rule. The proposed algorithm is shown to replicate the numerical results demonstrated on linear problems in previous works. Additional numerical examples of linear and nonlinear vibration and wave propagation are presented herein. Notably, the proposed algorithms show the same convergence rates for nonlinear problems as linear problems, and very high accuracy. Second-order time integration methods commonly used in commercial software produce significantly polluted acceleration responses for a common class of wave propagation problems. The high-order time integration schemes presented here perform noticably better at suppressing spurious high-frequency oscillations and producing reliable and useable acceleration responses.
title A high-order implicit time integration method for linear and nonlinear dynamics with efficient computation of accelerations
topic Numerical Analysis
65M22
J.2; G.1.3; G.1.8
url https://arxiv.org/abs/2409.13397