多体系统指标2运动方程HHT方法违约校正

作     者：马秀腾 翟彦博 谢守勇 Ma Xiuteng;Zhai Yanbo;Xie Shouyong

作者机构：西南大学工程技术学院重庆400715 现代汽车零部件技术湖北省重点实验室武汉430070 

出 版 物：《力学学报》 (Chinese Journal of Theoretical and Applied Mechanics)

年 卷 期：2017年第49卷第1期

页      面：175-181页

核心收录：

学科分类：08[工学] 080101[工学-一般力学与力学基础] 0801[工学-力学（可授工学、理学学位）] 

基　　金：国家自然科学基金(51605391) 现代汽车零部件技术湖北省重点实验室开放基金(2012-04)资助项目 

主　　题：运动方程 指标 HHT方法 Moore--Penrose广义逆 违约校正 

摘      要：采用Cartesian绝对坐标建模方法,完整约束多体系统运动方程是指标3的微分-代数方程(differentialalgebraic equations,DAEs),数值求解指标3的DAEs属于高指标问题,通过对位置约束方程求导,可使运动方程的指标降为2.位置约束方程求导得到的是速度约束方程.直接求解指标3的运动方程,速度约束方程得不到满足,而且高指标DAEs的数值求解存在一些问题.论文首先采用HHT(Hilber--Hughes--Taylor)直接积分方法求解降指标得到的指标2运动方程,此时速度约束方程参与离散计算,从机器精度上讲速度约束自然得到满足,而位置约束方程没有参与计算,存在 违约 .针对违约问题,采用基于Moore--Penrose广义逆理论的违约校正方法,消除位置约束方程的违约.指标2运动方程HHT方法违约校正,将HHT方法和违约校正方法很好地结合,在数值求解指标2运动方程的过程中,位置约束方程和速度约束方程都不存在违约问题,而且新方法没有引入新的未知数向量,离散得到的非线性方程组的方程数量与原指标2运动方程的方程数量相同,求解规模没有扩大.新方法的实用和有效性通过算例的数值实验得到验证,数值实验也说明新方法保持了HHT方法本身具有的数值阻尼可以控制和二阶精度的特性.最后从非线性方程组的求解规模和计算速度上与其他方法进行了比较分析,说明新方法的优势所在.