多体系统动力学微分-代数方程广义-α投影法

高效、稳定的微分-代数方程数值求解方法是多体系统动力学领域的关键问题之一。该文针对多体系统动力学指标3微分-代数方程,对目前多体系统动力学中引入的隐式时域逐步积分方法进行了深入研究,提出了适用于一般质量矩阵的广义-α-S法,并结合约束投影方法,构造了广义-α-S投影法。该方法既能较好地保持系统总能量,又能较高程度地同时满足位移约束、速度级约束和加速度级约束,并且在步长较大时可稳定求解,计算效率较高。