Construction of Higher-Order Accurate Time-Step Integration Algorithms by Equal-Order Polynomial Projection

Abstract

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.

Keywords

Unconditionally stable algorithms higher-order differential equations generalized Padé approximations time-step integration

Get full access to this article

View all access options for this article.

References

Adjerid, S. , Devine, K. D. , Flaherty, J. E. , and Krivodonova, L. , 2002, “A posteriori error estimation for discontinuous Galerkin solutions of hyperbolic problems,” Computer Methods in Applied Mechanics and Engineering 191, 1097-1112 .

Aharoni, D. and Bar-Yoseph, P. , 1992, “Mixed finite element formulations in the time domain for solution of dynamic problems,” Computational Mechanics 9, 359-374 .

Argyris, J. H. and Scharpf, D. W. , 1969, “Finite elements in time and space,” Journal of the Royal Aeronautical Society 73, 1041-1044 .

Argyris, J. H. , Vaz, L. E. , and Willam, K. J. , 1977, “Higher-order methods for transient diffusion analysis,” Computational Methods in Applied Mechanics and Engineering 12, 243-278 .

Bailey, C.D. , 1975, “Application of Hamilton’s law of varying action,” AIAA Journal 13, 1154-1157 .

Bailey, C. D. , 1976, “The method of Ritz applied to the equation of Hamilton,” Computer Methods in Applied Mechanics and Engineering 7, 235-247 .

Borri, M. and Bottasso, C. , 1993,“A general framework for interpreting time finite element formulations,” Computational Mechanics 13, 133-142 .

Borri, M. , Ghiringhelli, G. L. , Lanz, M. , Mantegazza, P. , and Merlini, T. , 1985, “Dynamic response of mechanical systems by a weak Hamilton formulation,” Computers and Structures 20, 495-508 .

Cannarozzi, M. and Mancuso, M. , 1995, “Formulation and analysis of variational methods for time integration of linear elastodynamics,” Computational Methods in Applied Mechanics and Engineering 127, 241-257 .

10.

Fried, I. , 1969, “Finite-element analysis of time-dependent phenomena,” AIAA Journal 7, 1170-1173 .

11.

Fung, T. C. , 1996, “Unconditionally stable higher-order accurate Hermitian time finite elements,” International Journal for Numerical Methods in Engineering 39, 3475-3495 .

12.

Fung, T. C. , 1998, “Complex-time-step Newmark methods with controllable numerical dissipation,” International Journal for Numerical Methods in Engineering 41, 65-93 .

13.

Fung, T. C. , 1999a, “Weighting parameters for unconditionally stable higher-order accurate time-step integration algorithms - Part 1. First-order equations,” International Journal for Numerical Methods in Engineering 45, 941-970 .

14.

Fung, T. C. , 1999b, “Weighting parameters for unconditionally stable higher-order accurate time-step integration algorithms - Part 2. Second-order equations,” International Journal for Numerical Methods in Engineering 45, 971-1006 .

15.

Fung, T. C. , 1999c, “Complex-time-step methods for transient analysis,” International Journal for Numerical Methods in Engineering 46, 1253-1271 .

16.

Fung, T. C. , 2000, “Weighting parameters for time-step integration algorithms with predetermined coefficients,” International Journal for Numerical Methods in Engineering 49, 1397-1430 .

17.

Fung, T. C. and Leung, A. Y. T., 1996, “On the accuracy of discontinuous Galerkin methods in the time domain,” Journal of Vibration and Control 2, 193-217 .

18.

Gellert, M. , 1978, “A new algorithm for integration of dynamic systems,” Computers and Structures 9, 401-408 .

19.

Hairer, E. , Norsett, S. P. , and Wanner, G. , 1987, Solving Ordinary Differential Equations I, Springer-Verlag, Berlin .

20.

Hairer, E. and Wanner, G. , 1991, Solving Ordinary Differential Equations II, Springer-Verlag, Berlin .

21.

Hodges, D. H. and Hou, L. J. , 1991, “Shape functions for mixed p-version finite elements in the time domain,” Journal of Sound and Vibration 145, 167-178 .

22.

Hughes, T. J. R. , 1987, The Finite Element Method: Linear Static and Dynamic Finite Element Analysis, Prentice-Hall, Englewood Cliffs, NJ .

23.

Hulbert, G. M. , 1992, “Time finite element methods for structural dynamics,” International Journal for Numerical Methods in Engineering 33, 307-331 .

24.

Hulbert, G. M. , 1994, “A unified set of single-step asymptotic annihilation algorithms for structural dynamics,” Computer Methods in Applied Mechanics and Engineering 113, 1-9 .

25.

Kim, S. J. , Cho, J. Y. , and Kim, W. D. , 1997, “From the trapezoidal rule to higher-order accurate and unconditionally stable time-integration method for structural dynamics,” Computer Methods in Applied Mechanics and Engineering 149, 73-88 .

26.

Kujawski, J. and Gallagher, R. H. , 1989, “A generalized least-squares family of algorithms for transient dynamic analysis,” Earthquake Engineering and Structural Dynamics 18, 539-550 .

27.

Kujawski, J. and Wiberg, N. E. , 1987, “Least-squares schemes for time integration of thermal problems,” International Journal for Numerical Methods in Engineering 24, 159-175 .

28.

Lambert, J. D. , 1973, Computational Methods in Ordinary Differential Equations, Wiley, NewYork .

29.

Lewis, R. W. and Bruch, I. C. Jr , 1974, “An application of least squares to one dimensional transient problems,” International Journal for Numerical Methods in Engineering 8, 633-647 .

30.

Möller, P. W. , 1993, “High-order hierarchical A- and L-stable integration methods,” International Journal for Numerical Methods in Engineering 36, 2607-2624 .

31.

Oden, J. T. , 1969, “A general theory of finite element: 2. Application,” International Journal for Numerical Methods in Engineering 1, 247-259 .

32.

Penry, S. N. and Wood, W. L. , 1985, “Comparison of some single-step methods for the numerical solution of the structural dynamics equation,” International Journal for Numerical Methods in Engineering 21, 1941-1955 .

33.

Peters, D. A. and Izadpanah, A. P. , 1988, “hp-version finite elements for the space-time domain,” Computational Mechanics 3, 73-88 .

34.

Riff, R. and Baruch, M. , 1984a, “Stability of time finite elements,” AIAA Journal 22, 1171-1173 .

35.

Riff, R. and Baruch, M. , 1984b, “Time finite element discretization of Hamilton’s law of varying action,” AIAA Journal 22, 1310-1318 .

36.

Segerlind, L. J. , 1989, “Weighted residual solutions in the time domain,” International Journal for Numerical Methods in Engineering 28, 679-685 .

37.

Simkins, T. E. , 1978, “Unconstrained variational statements for initial and boundary value problems,” AIAA Journal 16, 559-563 .

38.

van der Houwen, P. J. , 1977, Construction of Integration Formulas for Initial Value Problems, North-Holland, Amsterdam .

39.

Varga, R. S. , 1962, Matrix Iterative Analysis, Prentice-Hall, Englewood Cliffs, NJ .

40.

Wood, W.L. , 1990, Practical Time-Stepping Schemes, Clarendon Press, Oxford .

41.

Zienkiewicz, O. C. and Taylor, R. L. , 1991, The Finite Element Method, 4th edition, McGraw-Hill, New York .