摘要

采用了微分-代数方程求解的雷同策略,只在离散区段内满足无约束的积分,而在离散节点处则用约束冲量(Lagrange乘子向量)让轨道发生转折,给出了非完整约束动力系统的离散积分格式.数值例题的效果表明,本文给出的积分算法严格满足节点的非完整约束,具有良好的长时间积分性能.