多体系统动力学仿真向后微分公式算法

• 论文 •

多体系统动力学仿真向后微分公式算法

马秀腾1,2，翟彦博1，罗书强1

1.西南大学工程技术学院，重庆北碚400715;2.现代汽车零部件技术湖北省重点实验室，湖北武汉430070

收稿日期:2013-01-25 修回日期:2013-01-25 出版日期:2013-01-25 发布日期:2013-01-25

Dynamics simulation of multi-body system based on backward differentiation formulas

MA Xiu-teng1，2，ZHAI Yan-bo1，LUO Shu-qiang1

1.College of Engineering and Technology, Southwest University, Chongqing 400715, China；2.Hubei Key Laboratory of Advanced Technology for Automotive Components,Wuhan 430070,China

Received:2013-01-25 Revised:2013-01-25 Online:2013-01-25 Published:2013-01-25

摘要/Abstract

摘要： 针对多体系统动力学仿真中需要求解DAEs形式运动方程的问题，提出基于向后微分公式的DAEs算法。将已有的求解指标-3DAEs形式运动方程的NStiff算法推广到求解非完力学整系统指标-2DAEs形式的运动方程；提出了NStiff-v算法，分别求解完整约束多体系统指标-3DAEs和非完整力学系统指标-2DAEs形式的运动方程。离散过程中对约束方程进行缩放，消除了Newton-Raphson方法求解非线性方程组时由于步长较小而存在的Jacobian矩阵病态问题。通过曲柄滑块机构和Snakeboard模型两个实例，比较NStiff算法和NStiff-v算法的结果来验证算法的有效性，并比较了效率。算例说明算法求解DAEs时具有二阶精度。

关键词: 多体系统, 约束方程, 微分-代数方程, 向后微分公式, 算法

Abstract: Aiming at the problem of solving DAEs'motion equation in multi-body system dynamics simulation, new methods were proposed based on Backward Differentiation Formulas (BDF). The existing algorithm of NStiff which solved index-3DAEs form motion equation was extended to the index-2DAEs form motion equation of non-holonomic mechanics system. A new algorithm of NStiff-v was proposed to solve index-3DAEs form motion equation of holonomic constraint multi-body system and index-2DAEs form motion equation of non-holonomic mechanical system. During the discretization of motion equations, the constraint equations were scaled to eliminate ill-conditioning of Jacobian matrix caused by small time steps on solving nonlinear equations with Newton-Raphson method. Through two examples of slider-crank mechanism and Snakeboard model, the effectiveness of proposed algorithm was validated by comparing NStiff with NStiff-v. In addition, experiments illustrated that the algorithm had second-order accuracy in solving DAEs.

Key words: multi-body system, constraint equation, differential-algebraic equations, backward differentiation formulas, algorithms

中图分类号:

马秀腾1,2，翟彦博1，罗书强1. 多体系统动力学仿真向后微分公式算法[J]. .

MA Xiu-teng1，2，ZHAI Yan-bo1，LUO Shu-qiang1. Dynamics simulation of multi-body system based on backward differentiation formulas[J]. .

[1]	陈年,金涛,王建民. 基于变异反馈的临床路径优化[J]. 计算机集成制造系统, 2021, 27(9): 2556-2564.
[2]	陈志远,伍章俊,童珊珊,刘晓. 基于改进双归档进化算法的多目标动态软件项目调度[J]. 计算机集成制造系统, 2021, 27(9): 2565-2574.
[3]	原佳怡,朱锐,林雷蕾,李彤,郑明. 单触发序列突发漂移检测算法[J]. 计算机集成制造系统, 2021, 27(9): 2636-2646.
[4]	汤雅惠,李彤,朱锐,南峰涛,付会林. 基于执行者过程树的双维度遗传过程挖掘方法[J]. 计算机集成制造系统, 2021, 27(9): 2680-2690.
[5]	王沙沙,张则强,刘俊琦,陈凤. 多路径交互环形过道布置问题建模及改进蚁狮算法优化[J]. 计算机集成制造系统, 2021, 27(8): 2237-2247.
[6]	秦旋,房子涵,张赵鑫. 考虑资源约束的预制构件多目标生产调度优化[J]. 计算机集成制造系统, 2021, 27(8): 2248-2259.
[7]	尚正阳,顾寄南,王建平. 求解带能力约束车辆路径优化问题的改进模拟退火算法[J]. 计算机集成制造系统, 2021, 27(8): 2260-2269.
[8]	李国明,李军华. 带软时间窗随机需求车辆路径问题的算法研究[J]. 计算机集成制造系统, 2021, 27(8): 2270-2281.
[9]	顾九春,姜天华,朱惠琦. 多目标离散灰狼优化算法求解作业车间节能调度问题[J]. 计算机集成制造系统, 2021, 27(8): 2295-2306.
[10]	黎英杰,刘建军,陈庆新,毛宁. 多层级装配作业车间等量分批策略与调度算法[J]. 计算机集成制造系统, 2021, 27(8): 2307-2320.
[11]	冯春,张祎伟,黄成,姜文彪,武之炜. 双足机器人步态控制的深度强化学习方法[J]. 计算机集成制造系统, 2021, 27(8): 2341-2349.
[12]	冯国奇,崔东亮,代学武,俞胜平. 数据视角下神经网络增量学习支持的涡轮盘多目标优化[J]. 计算机集成制造系统, 2021, 27(8): 2393-2404.
[13]	汪逸晖,高亮. 乌鸦搜索算法的改进及其在工程约束优化问题中的应用[J]. 计算机集成制造系统, 2021, 27(7): 1871-1883.
[14]	刘景森,马义想,李煜. 改进鲸鱼算法求解工程设计优化问题[J]. 计算机集成制造系统, 2021, 27(7): 1884-1897.
[15]	赵转哲,叶国文,张宇,刘永明,张振,何康. 基于改进的二进制蚁狮算法的特征选择模型及应用[J]. 计算机集成制造系统, 2021, 27(7): 1898-1908.