多体系统动力学方程违约修正的数值计算方法

摘要：

多体系统动力学方程为微分代数方程,一般将其转化成常微分方程组进行数值计算,在数值积分的过程中约束方程的违约会逐渐增大.本文对具有完整、定常约束的多体系统,在修改的带乘子Lagrange正则形式的方程的基础上,根据Baumgarte提出的违约修正的方法,给出了一种多体系统微分代数方程违约修正法和系统的动力学方程的矩阵表达式.通过对曲柄-滑块机构的数值仿真,计算结果表明本文给出的方法在计算精度和计算效率上好于Baumgarte提出的两种违约修正的方法.

关键词：多体系统、完整定常约束、微分代数方程、违约修正

所属期刊栏目：24

分类号：O312(理论力学（一般力学）)

资助基金：国家自然科学基金10272008

在线出版日期：2007-04-09（万方平台首次上网日期，不代表论文的发表时间）

页数：共6页

页码：44-49