多体系统动力学变分数值积分方法及辛算法.pdf

摘 要 多体系统动力学研究主要包括多体系统动力学建模方法、多体系统动力学数值 积分方法以及通用或者专用软件的设计。近年来，多体系统的数学建模方法和理论 得到了飞速的发展，取得了大量的研究成果。但是由于多体系统动力学本身具有的 一些问题，使得高效、稳定、高精度的数值积分方法成为多体系统动力学新的研究 的热点。常用的多体系统动力学数值积分方法均基于线性局部分析的数值稳定性概 念，且离散后的数学模型没有考虑原连续模型固有的特性，严重限制了数值积分步 长，不适宜长时间仿真分析。几何数值积分方法是指离散模型保持原连续模型固由 特性或不变量的数值方法，这些不变量包括能量、动量、辛结构和李群结构等。基 于离散力学变分原理的变分数值积分方法能够自然的同时保持动量和能量这两种不 变量，辛算法则可以保持系统的辛结构。本文对多体系统动力学变分数值积分方法 和辛算法进行了深入的研究。主要包括以下几个方面：首先，本文详细阐述了多体 系统动力学中数学建模和数值积分的基本原理和方法。然后，介绍了多体系统动力 学中传统的数值积分方法，并通过对单摆系统的实验仿真结果分析得出结论。最后， 深入研究低阶变分数值积分方法和辛算法,并使用这两种方法对动力学系统进行实 验仿真，将得到的实验结果进行对比。 关键词：多体系统动力学；变分数值积分方法；辛算法 Abstract The research of multibody system dynamics contains: the dynamical models of multibody systems, the dynamical numerical analysis methods of multibody system and general or specific software design. In recent years , the dynamic numerical analysis methods and basic theory of multibody system is rapid development and get a lot of achievement. But due to some problems of multibody system dynamics itself makes that an efficient, stable, high-precision numerical integration in multibody system dynamic become a hot issue. All these methods were based on concepts of local linear analysis for numerical stabilization, and consistencies with their counterpart continuous models have not considered, thus not applicable to long term simulations with larger time steps. Geometric integrators are the numerical methods which can inherit the intrinsic invariants of the original continuous models, these invariants are energy, momentum, sympleticity, Lie group structure, etc.. But traditional geometric integrator usually preserves a specific invariant which can be improved in stabilization by preserving multiple invariants simultaneously. This issue is in-depth study in this paper. Include the following aspects: First