Directly self-starting higher-order implicit integration algorithms with flexible dissipation control for structural dynamics

2022 ◽  
Vol 389 ◽  
pp. 114274
Author(s):  
Jinze Li ◽  
Rui Zhao ◽  
Kaiping Yu ◽  
Xiangyang Li
2017 ◽  
Vol 84 (7) ◽  
Author(s):  
Wooram Kim ◽  
J. N. Reddy

For the development of a new family of implicit higher-order time integration algorithms, mixed formulations that include three time-dependent variables (i.e., the displacement, velocity, and acceleration vectors) are developed. Equal degree Lagrange type interpolation functions in time are used to approximate the dependent variables in the mixed formulations, and the time finite element method and the modified weighted-residual method are applied to the velocity–displacement and velocity–acceleration relations of the mixed formulations. Weight parameters are introduced and optimized to achieve preferable attributes of the time integration algorithms. Specific problems of structural dynamics are used in the numerical examples to discuss some fundamental limitations of the well-known second-order accurate algorithms as well as to demonstrate advantages of using the developed higher-order algorithms.


2017 ◽  
Vol 84 (7) ◽  
Author(s):  
Wooram Kim ◽  
J. N. Reddy

For the development of a new family of higher-order time integration algorithms for structural dynamics problems, the displacement vector is approximated over a typical time interval using the pth-degree Hermite interpolation functions in time. The residual vector is defined by substituting the approximated displacement vector into the equation of structural dynamics. The modified weighted-residual method is applied to the residual vector. The weight parameters are used to restate the integral forms of the weighted-residual statements in algebraic forms, and then, these parameters are optimized by using the single-degree-of-freedom problem and its exact solution to achieve improved accuracy and unconditional stability. As a result of the pth-degree Hermite approximation of the displacement vector, pth-order (for dissipative cases) and (p + 1)st-order (for the nondissipative case) accurate algorithms with dissipation control capabilities are obtained. Numerical examples are used to illustrate performances of the newly developed algorithms.


2019 ◽  
Vol 17 (08) ◽  
pp. 1950048 ◽  
Author(s):  
Wooram Kim ◽  
Jin Ho Lee

Two families of higher-order accurate time integration algorithms are numerically tested by using various nonlinear problems of structural dynamics, and the numerical results obtained from them are compared. To be specific, the higher-order algorithms of Kim and Reddy and the higher-order algorithms of Fung are used for this study. In linear analyses, these two different families of higher-order algorithms do not present noticeable differences. However, performances of these algorithms are quite different when they are applied to various nonlinear dynamic problems. For the numerical tests, well-known nonlinear problems are selected from the past studies. For the completeness, the two families of algorithms are briefly reviewed, and their advantageous computational structures are also explained.


2005 ◽  
Vol 11 (1) ◽  
pp. 19-49 ◽  
Author(s):  
T. C. Fung

This paper presents a new framework to construct higher-order accurate time-step integration algorithms based on weakly enforcing the differential/integral relation. The dependent variable and its time derivatives are assumed to be polynomials of equal order. A differential equation is then transformed into an algebraic equation directly. The main issue is how to approximate the integral of a polynomial by another polynomial of the same degree. Various methods to determine the optimal representation (or projection) are considered. It is shown that to reproduce numerical results equivalent to the Padé or generalized Padé approximations, the coefficients of the optimal polynomial representation are related to the weighting parameters derived previously for time-step integration algorithms with predetermined coefficients. A special feature for the present formulation is that the same procedure can be used to solve first-, second-, and even higher-order non-homogeneous initial value problems in a unified manner. The resultant algorithms are higher-order accurate and unconditionally A-stable with controllable numerical dissipation. It is also shown that for the numerical results to maintain higher-order accuracy at the end of a time interval, the higher-order terms in the excitation have to be projected as polynomials of lower degree within the present framework as well. Numerical examples are given to illustrate the validity of the present formulations.


Sign in / Sign up

Export Citation Format

Share Document