含摩擦柱铰链平面多体系统动力学的建模和数值方法-Engineering.pdf

29 卷 5 期 Vol.29 No.5 工 程 力 学 2012 年 5 月 May 2012 ENGINEERING MECHANICS 193 文章编号：1000-4750(2012)05-0193-07 含摩擦柱铰链平面多体系统动力学的 建模和数值方法 庄方方，王 琪 (北京航空航天大学动力学与控制系，北京 100191) 摘 要：以含摩擦柱铰链平面多体系统为研究对象，建立其动力学方程并给出相应的数值计算方法。首先，建立 了含摩擦转动柱铰链的力学模型。在此基础上，应用 一类 Lagrange 方程给出了该类系统的动力学方程，将 Lagrange 乘子与柱铰链的法向约束力建立了对应关系，并给出了柱铰链摩擦力的广义力。由于摩擦力的存在，使 得该方程是关于 Lagrange 乘子的分段连续的非线性代数方程组，该文对此采用混合算法：对于连续段(物体相对 转动的角速度不为零时) ，采用拟牛顿算法和龙格-库塔法求解方程；在不连续点(物体相对转动的角速度为零时) ， 通过粒子群算法(PSO)、试算法和龙格-库塔法求解方程，克服了方程在不连续处Lagrange 乘子(法向约束力) 的初 值不易选取的困难。最后，通过算例说明了该算法的有效性和可行性。 关键词：多体系统；转动铰链；Lagrange 乘子；摩擦；粒子群算法 中图分类号：O313.7 文献标志码：A MODELING AND NUMERICAL ALGORITHM FOR PLANAR MULTIBODY SYSTEM WITH FRICTION ON REVOLUTE JOINTS ZHUANG Fang-fang , WANG Qi (Department of Dynamics and Control, Beijing University of Aeronautics and Astronautics, Beijing 100191, China) Abstract: The modeling and numerical algorithm for a planar multi-body system with friction on revolute j oints are developed. Based on the model of a revolute joint, the dynamic equations of the system are derived by the first kind of Lagrange s equations, then the relationship between Lagrange multipliers and the normal constraint force acting on the joint is obtained and the generalized force of friction is given. For the equations are piecewise continuous nonlinear algebraic equations about Lagrange multipliers, Quasi-Newton and Runge-Kutta algorithms are used to simulate the dynamical systems during the continuous interval (the relative angular velocity of the body does not equal to zero), and Particle Swarm Optimization, the trial and error method and Runge-Kutta algor