非完整系统动力学仿真θ_1方法研究

θ1方法是一种直接时间积分算法,主要用于结构动力学仿真时运动方程的求解,方程形式为二阶常微分方程(ODEs)。对于非完整系统动力学仿真,从微分-代数方程(DAEs)的角度看,系统的运动方程是二阶DAEs。基于此,提出非完整系统仿真的θ1方法,也就是数值求解指标-2的运动方程—DAEs的新算法。通过对双轮机器人θ1方法仿真结果与DASSL和Radau5算法结果的比较,验证新算法的有效性。数值实验也说明θ1方法求解非完整系统DAEs时具有2阶精度。