EI、Scopus 收录
中文核心期刊
张润森, 王琪. 浮放物体平面多刚体动力学建模与算法研究[J]. 力学学报, 2017, 49(6): 1370-1379. DOI: 10.6052/0459-1879-17-235
引用本文: 张润森, 王琪. 浮放物体平面多刚体动力学建模与算法研究[J]. 力学学报, 2017, 49(6): 1370-1379. DOI: 10.6052/0459-1879-17-235
Zhang Runsen, Wang Qi. RESEARCH ON MODELING AND NUMERICAL METHOD OF FREE STANDING BODY ON PLANAR RIGID MULTIBODY DYNAMICS[J]. Chinese Journal of Theoretical and Applied Mechanics, 2017, 49(6): 1370-1379. DOI: 10.6052/0459-1879-17-235
Citation: Zhang Runsen, Wang Qi. RESEARCH ON MODELING AND NUMERICAL METHOD OF FREE STANDING BODY ON PLANAR RIGID MULTIBODY DYNAMICS[J]. Chinese Journal of Theoretical and Applied Mechanics, 2017, 49(6): 1370-1379. DOI: 10.6052/0459-1879-17-235

浮放物体平面多刚体动力学建模与算法研究

RESEARCH ON MODELING AND NUMERICAL METHOD OF FREE STANDING BODY ON PLANAR RIGID MULTIBODY DYNAMICS

  • 摘要: 采用非光滑多体系统动力学的方法研究浮放物体与基础平台组成的多体系统,建立其非光滑接触的动力学方程与数值算法.浮放物体由主体部分和支撑腿组成,其间通过含黏弹性阻力偶的转动铰连接.支撑腿与基础平台间的接触力简化为接触点的法向接触力和摩擦力,采用扩展的赫兹接触力模型描述接触点的法向接触力,采用库伦干摩擦模型描述其摩擦力.采用笛卡尔坐标系下的位形坐标作为系统的广义坐标.首先,将基础平台运动看作非定常约束,用第一类拉格朗日方程建立系统的动力学方程,并采用鲍姆加藤约束稳定化的方法解决违约问题.随后给出基于事件驱动法和线性互补方法的数值算法.当相对切向速度为零时,构造静滑动摩擦力的正负余量和正、负向加速度的互补关系,从而将接触点黏滞——滑移切换的判断以及静滑动摩擦力的计算转化为线性互补问题进行求解,并采用Lemke算法求解线性互补问题.最后,通过数值仿真选择合适的步长;通过仿真结果说明浮放物体运动中存在的黏滞-滑移切换现象以及基础平台运动、质心位置对浮放物体运动的影响.

     

    Abstract: A multibody system composed of free standing body and the basic platform is investigated by non-smooth dynamics of multibody system. Dynamic equations and numerical method of the system with non-smooth contacts are proposed. The free standing body consists of main body and supporting legs, which are connected by revolute joints with viscoelastic moments. The contact forces between free standing body and the basic platform are simplified as normal forces and frictional forces of contact points. Moreover, the modified Hertz contact model and Coulomb's law for dry friction are employed to describe normal forces and frictional forces, respectively. And the configuration coordinates of Cartesian coordinate system are used as the generalized coordinates. Firstly, the system's dynamic equations are established by Lagrange's equations of the first kind and the motion of basic platform is regarded as a rheonomous constraint. The problem of constraints violations is solved by Baumgarte stabilization method. Secondly, the numerical method of the multibody system are proposed, which is based on the event-driven schemes and linear complementarity formulations. The complementary formulations of friction saturations and the relative accelerations in the tangential are given, while the relative velocities in the tangential of the contact points are equal to zero. Therefore the judgements of stick-slip transitions for contact points and the solutions of frictional forces in stick situation could be solved as a linear complementarity problem. And the linear complementarity problem is solved by Lemke's algorithm. Finally, an appropriate step is chosen by the simulation. Then the numerical simulations denote the stick-slip phenomenon and the influence of basic platform as well as mass centre's position.

     

/

返回文章
返回