论文部分内容阅读
在分析多体动力学仿真计算中广泛使用的增广法基础上,提出了一种基于稀疏矩阵技术的改进算法。算法首先修改了拉格朗日乘子的求解方程,并利用系数矩阵的稀疏性进行排序和符号分解,获得最佳计算数据结构,从而快速求解拉格朗日乘子;针对违约现象,建立直接违约修正方程,利用上述数据结构快速求解违约修正值。最后,通过一个算例验证了该算法的有效性,并与传统增广法进行比较,显示了该算法快速和高精度的特点。