%0 Journal Article %T 基于\pmb\theta_1方法的多体动力学数值算法研究 %A 马秀腾 %A 翟彦博 %A 罗书强 %J 力学学报 %P 931-938 %D 2011 %R 10.6052/0459-1879-2011-5-lxxb2010-730 %X 将结构动力学领域的\theta_1方法拓展到数值求解多体系统运动方程------微分--代数方程(DAEs),分别求解指标-3DAEs形式的运动方程和指标-2超定DAEs(ODAEs)形式的运动方程.通过数值算例验证了方法的有效性,并得到\theta_1方法中参数\theta_1的选取与数值耗散量之间的关系.数值算例还说明对于同一个多体系统,采用指标-3的DAEs描述时存在速度违约,用指标-2的ODAEs描述时,从计算机精度上讲,位置和速度约束方程同时满足,并且\theta_1方法在求解非保守系统DAEs和ODAEs形式的运动方程时都具有2阶精度.最后\theta_1方法与其他直接积分法求解DAEs和ODAEs形式运动方程的CPU时间进行了比较. %K \theta %K _1 %K 方法 %K 多体系统 %K 微分--代数方程 %K (DAEs) %K 数值耗散 %K 2阶精度 %U http://lxxb.cstam.org.cn/CN/abstract/abstract142692.shtml