Abstract: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.

A trigonometrically fitted explicit Numerov-type method for second-order initial value problems with oscillating solutions

Trigonometrically Fitted Explicit Numerov-type Method for Periodic IVPs with Two Frequencies

A trigonometrically fitted explicit hybrid method for the numerical integration of orbital problems

New embedded pairs of explicit Runge–Kutta methods with FSAL properties adapted to the numerical integration of oscillatory problems

A Parametric Six-Step Method for Second-Order IVPs with Oscillating Solutions

A Phase-Fitted Runge-Kutta-Nyström method for the Numerical Solution of Initial Value Problems with Oscillating Solutions

A New Pair of Explicit ARKN Methods for the Numerical Integration of General Perturbed Oscillators

Extended Phase Properties and Stability Analysis of RKN-type Integrators for Solving General Oscillatory Second-Order Initial Value Problems

A class of two-step explicit methods for periodic IVPs

Efficient Multistep Algorithms for First-Order IVPs with Oscillating Solutions: II Implicit and Predictor–Corrector Algorithms

A Multistep Method for Integration of Perturbed and Damped Second-Order ODE Systems

A New Methodology for the Development of Efficient Multistep Methods for First-Order IVPs with Oscillating Solutions

A Filon-type asymptotic approach to solving highly oscillatory second-order initial value problems.

A Numerical Method For Solving The One-Dimensional Vlasov-Poisson Equation In Phase Space

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

Reproducing Kernel-Based Piecewise Methods for Efficiently Solving Oscillatory Systems of Second-Order Initial Value Problems

Functionally-fitted energy-preserving methods for solving oscillatory nonlinear Hamiltonian systems

A new high efficient and high accurate Obrechkoff four-step method for the periodic nonlinear undamped Duffing's equation

Efficient second-order semi-implicit finite element method for fourth-order nonlinear diffusion equations

Numerical methods for highly oscillatory Volterra integral equations with general oscillators

High asymptotic order numerical methods for highly oscillatory ODEs with large initial data