一种基于Hamel形式的无条件稳定动力学积分算法

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.