论文部分内容阅读
以数值积分求解约束力学系统指标-2的微分-代数方程(DAEs)为例,介绍了直接积分方法中的Newmark方法和HHT-α方法在约束力学系统仿真中的应用。主要从算法有无数值耗散、耗散量的多少以及如何控制耗散量的角度研究了这两种方法,并与BDF方法、Radau5方法以及ode45的数值耗散量进行了比较,得到的结果对于这些算法的使用,对于商业动力学仿真软件"求解器"中算法的选取有重要的指导意义,同时数值实验也表明辛算法也可以应用到非保守系统,虽然此时辛结构不再保持,但算法"无"数值耗散,可以更好的反映系统