EI、Scopus 收录
中文核心期刊

多体系统摩擦接触问题的DAE-LCP方法

DAE-LCP mixed method for multibody system dynamics with frictional contacts

  • 摘要: 当多体系统的约束全部是摩擦接触时, 其动力学问题可归结为一个常微分方程(ordinarydifferential equation, ODE)与线性互补问题(linear complementarityproblem, LCP)的混合动力学问题. 如果除了摩擦接触之外还增加了光滑的双边约束, 则需要将ODE-LCP混合动力学模型推广为微分代数方程(differential algebra equation, DAE)与LCP的混合动力学模型. 该文采用DAE与LCP混合动力学方法求解不考虑碰撞但同时含有持续摩擦接触及光滑等式约束的多体系统动力学问题. 在建立系统动力学模型时,首先将含摩擦的约束从系统中移去得到基本动力学系统. 由于基本系统中带有等式约束, 所以基本系统的动力学方程为一组DAE. 结合基本系统的DAE与约束的互补条件便可以得到DAE-LCP混合动力学模型. 数值计算采用基于DAE与LCP的步进(time-stepping)算法, 将系统动力学方程及其约束离散化并转化为一个混合LCP进行求解. 该算法无需进行滞-滑状态检测, 避免了事件检测导致的繁复计算. 利用所提方法对典型机构的非光滑非线性特征进行了数值分析,验证了该文方法的正确有效性.

     

    Abstract: Presently, dynamics of nonsmooth multibody systems is a hot research topic.The usual approach in treating such systems is to derive basic system fromthe original system by removing the nonsmooth constraints firstly. TheLagrange equations of the second kind of basic system combine with thecomplementarity condition of the nonsmooth constraints to set up at eachdiscrete moment in time a Linear Complementarity Problem (LCP).This article focuses on the problem of dynamic modeling and numericalsimulating of multibody systems with friction contacts. By neglecting theclearance and the effect of impact between rigid bodies and constraints, thestate variables in the differential equations are continuous. Due to theset-value mapping characteristic of dry friction forces, the differentialequations of motion have discontinuous right-hand vector fields, thereforeto allow our system to be classified as a Filippov system. In addition tothe friction constraints, our model also incorporates frictionless bilateralones. If the simulation has only friction constraints, then the problem isan ODE-LCP model. Combining friction contacts with frictionless bilateralconstaints, the ODE-LCP model has to be extended to DAE-LCP (DAE,Differential Algebra Equation) mixed model. In order to obtain the model,the basic system is derived from the original system by removing thefriction constraints firstly. Because the equality constraints are retainedin the basic system, the dynamic model of basic system is a set of DAE. Withthe aid of constraint Jacobian matrix, the normal contact forces andtangential friction forces of nonsmooth constraints, which obey thecomplementarity contact laws, are added to the DAE of the basic system toobtain the DAE-LCP mixed model.Approaches used in the past for simulating rigid multibody dynamics withfriction contacts include piecewise DAE approaches, acceleration-force LCPapproaches, and velocity-impulse LCP-based time-stepping methods.Recognizing that the nature of the frictional constraint can inducesick-slip motion, the last approach is used in this work, which has theadvantage that it does not suffer from the detection for stick-sliptransition that could appear in the first two approaches. This framework isbased on a LCP, but it is different from acceleration-force LCP approachesthat attempt to find the accelerations of the bodies.Our approach considersimpulses and velocities as the fundamental unknowns. Acceleration-force LCPapproaches solve for accelerations from the dynamics equations and then usethe accelerations in an integration procedure. Because the complementaritylaw between acceleration and friction surplus is valid only when therelative tangential velocity is zero, zero crossing detection for velocityis required. However the complementarity law between velocity and frictionsurplus remains valid through a full-range of motion, so in contrast toacceleration-force schemes, the velocity-impulse methods need noevent-detection. In the new framework, the integration and dynamicalresolution steps are combined. The main achievement of this approach is thatit has solutions for any configuration. As the time-step tends to zero, asubsequence of the numerical solutions approaches the solution of adifferential inclusion.Our method is carried out in a numerical example, and the simulation resultsindicate that this method is effective.

     

/

返回文章
返回