论文部分内容阅读
多体系统动力学微分／代数混合方程组又称为Ｅｕｌｅｒ－Ｌａｇｒａｎｇｅ方程。其数值积分的困难之一是由违约引起的数值不稳定。基于对约束方程左部的Ｔｙｌｏｒ展开，根据积分步长提出了一种能对约束误差自动修正的小扰动我稳定方法。该方法大大改善了传统违修正法的数值性态，并具有简单、实用、高效的特点。最后对该方法与传统增广方法及其违约的修正方法进行了数值比较。