首先利用Hamilton原理对耦合结构进行建模,然后利用有限元方法将空间连续模型离散化,得到有限元模型,然后将模型导入到Hamilton系统中,获得Hamilton正则方程。在建模的基础上,采用半隐式辛Runge-Kutta（SRK）算法对Hamilton系统下的模型进行计算,并与传统Runge-Kutta方法进行了比较。数值仿真结果表明,采用二阶的半隐式SRK算法能够长时间保持系统的能量守恒,是一种保结构算法,而传统的Runge-Kutta方法是一种耗散算法,在求解初期具有高的精度,但是不能保持系统解的长期稳定性。从仿真过程中还可以得到一个非常有意义的结果,二阶的半隐式SRK算法对步长的要求低于传统的四阶Runge-Kutta方法。

多体系统动力学的微分/代数方程求解一般是所谓的指标-3问题,是十分困难的,可以说,目前还没有使人非常满意的关于它的数值积分方法.多体系统动力学的微分/代数方程的辛算法,是近几年出的新的数值方法,一般它具有精度高、数值稳定性等优点.笔者建立了约束多体系统动力学的微分/代数形式的约束正则方程形式,利用Runge-Kutta法合成辛算法对约束多体系统的约束哈密顿形式的方程进行仿真研究取得了较好的结果.