基于HHT-α方法的非完整系统运动方程数值算法研究.PDFVIP

( ) 第 卷第 期 西 南 大 学 学报 自然科学版 年 月 33 5 2011 5 ( ) Vol.33 No.5 JournalofSouthwestUniversit NaturalScienceEdition Ma 2011 y y 文章编号: ( ) 1673 9868201105 0151 04 基于HHT-α方法的非完整 系统运动方程数值算法研究① 马秀腾, 翟彦博, 罗书强 , 西南大学 工程技术学院 重庆 400715 : ( ) , 摘要 非完整系统运动方程是微分 代数方程 DAEs.基于直接时间积分中的 HHT-α方法 将其拓展应用到数值求 , , 解 形式的非完整系统运动方程 针对一个 简化模型 通过与 算法的求解结果进行比较 验 DAEs . Snakeboard Radau5 证了此方法的正确性 数值试验也说明 方法求解 形式的非完整系统运动方程时具有二阶精度 . HHT-α DAEs . : ; ; ( ); 关 键 词 非完整系统 运动方程 微分 代数方程 方法 DAEs HHT-α 中图分类号: ; 文献标志码: O313.7 TB115 A , , 预先给定的制约自由系统状态的条件称为约束 约束可分为完整约束和非完整约束 完整约束是指几 , 何约束和可积分的微分约束 非完整约束是指不可积分的微分约束 约束系统运动方程一般可表示成微分 . 代数方程( )的形式 DAEs . [] 1 只有完整约束的力学系统运动方程是指标 的 ,此 的数值算法一直是研究的热点问 3 DAEs DAEs [ ] [ ] 2-3