AN UNCONDITIONALLY STABLE DYNAMICAL INTEGRATION ALGORITHM BASED ON HAMEL’S FORMALISM

被引：0
|
作者：
Gu W. [1 ]
Liu C. [2 ,3 ]
An Z. [4 ]
Shi D. [5 ,6 ]
机构：
[1] School of Information and Electronics, Beijing Institute of Technology, Beijing
[2] CAEP Software Center for High Performance Numerical Simulation, Beijing
[3] Institute of Applied Physics and Computational Mathematics, Beijing
[4] School of Aerospace Engineering, Beijing Institute of Technology, Beijing
[5] School of Mathematics and Statistics, Beijing Institute of Technology, Beijing
[6] Beijing Key Laboratory on MCAACI, Beijing
来源：
Lixue Xuebao/Chinese Journal of Theoretical and Applied Mechanics | 2022年 / 54卷 / 09期
关键词：
generalized; α; method; Hamel-generalized; Hamel’s field variational integrators; Lie group; moving frame;
D O I：
10.6052/0459-1879-22-131
中图分类号：
学科分类号：
摘要：
Time integration algorithm is a key issue in solving dynamical system. An unconditionally stable Hamel generalized α method is proposed to solve the instability issue arising in the time integration of dynamic equations and to eliminate the pseudo high order harmonics incurred by the spatial discretization of finite element simultaneously. Therefore, the development of numerical integration algorithm to solve the above-mentioned problems has important theoretical and application value. The algorithm proposed in this paper is developed based on the moving frame method and Hamel’s field variational integrators along with the strategy to construct an unconditionally stable Hamel generalized α method. It is shown that a new numerical formalism with higher accuracy can be derived under the same framework of the unconditional stable algorithm established through a special variational formalism and variational integrators. The above-mentioned formalism can be extended from general linear space to Lie group by utilizing the moving frame method and the Lie group formalism of the Hamel generalized α method has been obtained. Both the convergence and stability of the algorithm are discussed, and some numerical examples are presented to verify the conclusion. It is demonstrated by the theoretical analysis that the Hamel generalized α method proposed in the paper is unconditionally stable, second-order accurate and can quickly filter out pseudo high-frequency harmonics. Both conventional and proposed methods have been applied to numerical examples respectively. Comparisons between results of numerical examples show that the aforementioned advantages of the proposed method in terms of accuracy, dissipation and stability are tested and verified. At the same time, it can be developed that new numerical integration algorithms with even higher order accuracy. The scheme can also be proposed, which is suitable for both general linear space and Lie group space. A new way for constructing variational integrators is also obtained in this paper. © 2022 Chinese Journal of Theoretical and Applied Mechanics Press. All rights reserved.
引用
收藏
页码：2577 / 2587
页数：10
相关论文
共 41 条
  • [1] Newmark NM., A method of computation for structural dynamics, Journal of the Engineering Mechanics Division, 85, 3, pp. 67-94, (1959)
  • [2] Chung J, Hulbert G., A time integration algorithm for structural dynamics with improved numerical dissipation: The generalized-α method, ASME Journal of Applied Mechanics, 60, pp. 371-375, (1993)
  • [3] Wood WL, Bossak M, Zienkiewicz OC., An alpha modification of Newmark's method, International Journal for Numerical Methods in Engineering, 15, 10, pp. 1562-1566, (1980)
  • [4] Hilber HM, Hughes TJR, Taylor RL., Improved numerical dissipation for time integration algorithms in structural dynamics, Earthquake Engineering and Structural Dynamics, 5, 3, pp. 283-292, (1977)
  • [5] Bazilevs Y, Calo VM, Cottrell JA, Et al., Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows, Computer Methods in Applied Mechanics and Engineering, 197, 1-4, pp. 173-201, (2007)
  • [6] Behnoudfar P, Calo VM, Deng Q, Et al., A variationally separable splitting for the generalized-α method for parabolic equations, International Journal for Numerical Methods in Engineering, 121, 5, pp. 828-841, (2020)
  • [7] Gomez H, Hughes TJR, Nogueira X, Et al., Isogeometric analysis of the isothermal Navier–Stokes–Korteweg equations, Computer Methods in Applied Mechanics and Engineering, 199, 25-28, pp. 1828-1840, (2010)
  • [8] Behnoudfar P, Deng Q, Calo VM., Higher-order generalized-α methods for hyperbolic problems, Computer Methods in Applied Mechanics and Engineering, 378, (2021)
  • [9] Yen J, Petzold L, Raha S., A time integration algorithm for flexible mechanism dynamics: The DAE α-method, Computer Methods in Applied Mechanics and Engineering, 158, 3-4, pp. 341-355, (1998)
  • [10] Arnold M, Bruls O., Convergence of the generalized-α scheme for constrained mechanical systems, Multibody System Dynamics, 18, 2, pp. 185-202, (2007)
12345