|
|
力学学报 2011
基于\pmb\theta_1方法的多体动力学数值算法研究DOI: 10.6052/0459-1879-2011-5-lxxb2010-730, PP. 931-938 Keywords: \theta,_1,方法,多体系统,微分--代数方程,(DAEs),数值耗散,2阶精度 Abstract: 将结构动力学领域的\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时间进行了比较.
|
