含摩擦滑移铰及驱动约束多刚体系统数值算法

研究了具有驱动约束及非光滑滑移铰多体系统动力学方程的建模与数值计算方法．将驱动约束视为非定常约束，非光滑滑移铰视为双边定常约束，滑移铰的摩擦模型采用库仑摩擦模型；应用第一类Lagrange方程建立系统的动力学方程，应用距离函数建立滑移铰的约束方程；将线性互补方法和Baumgarte约束稳定化方法引入，以解决滑移铰法向约束力的计算以及约束方程违约问题．最后应用曲柄摇杆机构作为算例，说明该方法的有效性．     

Modeling and analysis of planar rigid multibody systems with translational clearance joints based on the non-smooth dynamics approach

The main purpose of this paper is to present and discuss a methodology for a dynamic modeling and analysis of rigid multibody systems with translational clearance joints. The methodology is based on the non-smooth dynamics approach, in which the interaction of the elements that constitute a translational clearance joint is modeled with multiple frictional unilateral constraints. In the following, the most fundamental issues of the non-smooth dynamics theory are revised. The dynamics of rigid multibody systems are stated as an equality of measures, which are formulated at the velocity-impulse level. The equations of motion are complemented with constitutive laws for the normal and tangential directions. In this work, the unilateral constraints are described by a set-valued force law of the type of Signorini’s condition, while the frictional contacts are characterized by a set-valued force law of the type of Coulomb’s law for dry friction. The resulting contact-impact problem is formulated and solved as a linear complementarity problem, which is embedded in the Moreau time-stepping method. Finally, the classical slider-crank mechanism is considered as a demonstrative application example and numerical results are presented. The results obtained show that the existence of clearance joints in the modeling of multibody systems influences their dynamics response.     

A differential approach for modeling revolute clearance joints in planar rigid multibody systems

A non-penetration approach of frictional contact analysis is presented for modeling revolute clearance joints of planar rigid multibody systems. In the revolute clearance joint, the motion modes of the journal are divided into three categories, namely, the free motion, collision, and permanent contact modes. The switch between different contact modes is identified by the state of the journal and bearing, including the gap and the normal relative velocity. When impact in the revolute clearance joint is detected, the collision process is simulated by the impulse-based differential approach, where Stronge’s improved model for restitution is employed to determine the relative velocity after impact. Instead of algebraic equations, the impact process is described by a set of ordinary differential equations (ODEs), which avoids solving complementarity problems. Moreover, in the permanent contact mode, the constraint-based approach and modified Coulomb’s friction law are adopted. The permanent contact mode maintains for most of the time and the governing ODEs are non-stiff. There is general agreement that the constraint-based approach is more efficient than the force-based method. A slider–crank mechanism with a revolute clearance joint is considered as a demonstrative application example where the comparison with the continuous contact force model is investigated.     

Dynamics of spatial flexible multibody systems with clearance and lubricated spherical joints

A computational methodology for analysis of spatial flexible multibody systems, considering the effects of the clearances and lubrication in the system spherical joints, is presented. The dry contact forces are evaluated through a Hertzian-based contact law, which includes a damping term representing the energy dissipation. The frictional forces are evaluated using a modified Coulomb’s friction law. In the case of lubricated joints, the resulting lubricant forces are derived from the corresponding Reynolds’ equation. An absolute nodal formulation is utilized in flexible body formulation. The generalized-α method is used to solve the resulting equations of motion. The effectiveness of the methodology is demonstrated by two numerical examples.     

Modeling and simulation of longitudinal dynamics coupled with clutch engagement dynamics for ground vehicles

During the engagement of the dry clutch in automotive transmissions, clutch judder may occur. Vehicle suspension and engine mounts couple the torsional and longitudinal models, leading to oscillations of the vehicle body that are perceived by the driver as poor driving quality. This paper presents an effective formulation for the modeling and simulation of longitudinal dynamics and powertrain torsional dynamics of the vehicle based on non-smooth dynamics of multibody systems. In doing so friction forces between wheels and the road surface are modeled along with friction torque in the clutch using Coulomb’s friction law. First, bilateral constraint equations of the system are derived in Cartesian coordinates and the dynamical equations of the system are developed using the Lagrange multiplier technique. Complementary formulations are proposed to determine the state transitions from stick to slip between wheels and road surface and from the clutch. An event-driven scheme is used to represent state transition problem, which is solved as a linear complementarity problem (LCP), with Baumgarte’s stabilization method applied to reduce constraint drift. Finally, the numerical results demonstrate that the modeling technique is effective in simulating the vehicle dynamics. Using this method stick-slip transitions between driving wheel and the road surface and from the clutch, as a form of clutch judder, are demonstrated to occur periodically for certain values of the parameters of input torque from engine, and static and dynamic friction characteristics of tire/ground contact patch and clutch discs.     

Frictional contact analysis of spatial prismatic joints in multibody systems

The contact analysis of spatial prismatic joints remains a hard problem due to its complex nature. In this paper, a methodology for the frictional contact analysis of rigid multibody systems with spatial prismatic joints is presented, which is free of calculating the relative motion between the slider and guide, and is particularly suitable to the case of clearances being tiny. Under the assumption of the slider and guide being rigid, we prove that all types of contacts in the joint can be converted to point-to-point contacts. At each of the candidate points, two gap functions are introduced. However, in the proposed method, not the values of these gap functions but the relations between them are essential. In view of the non-colliding contacts being predominant when clearances of joints are tiny, we formulate the contact forces in terms of resultant frictional forces in the joint, resulting in a linear complementarity problem. By the proposed method, details about the contacts including the impact instants can be obtained, although impacts are not taken into consideration explicitly, as indicated by the numerical examples in this paper.     

多体系统中典型铰的摩擦力计算模型

在铰内间隙很小的前提下,多体系统中铰仍具有运动学约束作用.但由于铰内接触形式与系统状态相关,铰内摩擦力与约束反力之间具有复杂的函数关系.本文在假设铰内接触为刚性接触的前提下,基于分布接触反力与点接触反力之间的等效关系,给出了几种典型铰内摩擦力的计算模型,并通过数值算例验证了所提模型的正确性.     

Dynamic Analysis for Planar Multibody Mechanical Systems with Lubricated Joints

The dynamic analysis of planar multibody systems with revolute clearance joints, including dry contact and lubrication effects is presented here. The clearances are always present in the kinematic joints. They are known to be the sources for impact forces, which ultimately result in wear and tear of the joints. A joint with clearance is included in the multibody system much like a revolute joint. If there is no lubricant in the joint, impacts occur in the system and the corresponding impulsive forces are transmitted throughout the multibody system. These impacts and the eventual continuous contact are described here by a force model that accounts for the geometric and material characteristics of the journal and bearing. In most of the machines and mechanisms, the joints are designed to operate with some lubricant fluid. The high pressures generated in the lubricant fluid act to keep the journal and the bearing surfaces apart. Moreover, the lubricant provides protection against wear and tear. The equations governing the dynamical behavior of the general mechanical systems incorporate the impact force due to the joint clearance without lubricant, as well as the hydrodynamic forces owing to the lubrication effect. A continuous contact model provides the intra-joint impact forces. The friction effects due to the contact in the joints are also represented. In addition, a general methodology for modeling lubricated revolute joints in multibody mechanical systems is also presented. Results for a slider-crank mechanism with a revolute clearance joint between the connecting rod and the slider are presented and used to discuss the assumptions and procedures adopted.     

Docking dynamics between two spacecrafts with APDSes

Contact/impact models are established mainly for the capturing phase of a docking system. All the interactions between them are supposed as point contacts/impacts. For contacts, instead of using compliant models, constraint equations are formulated on the thorough investigation of the contour shapes of the docking system. All the contact constraints are originated from the sufficient and necessary conditions of point contacts. The constraints and Coulomb’s frictional law are incorporated into the contact dynamic equations of the system to obtain the contact forces by Lagrangian multipliers. For impacts, multiple impacts are considered in the term of the complex configuration of the docking dynamics. Impulsive equations are deduced, and the distributing law which is also called the LZB method is applied for multiple impacts. Finally, simulation is realized and the numerical results for one working condition are presented.     

含铰链间隙的刚-柔机械臂动力学模型

根据Hertz接触定律和Coulomb摩擦定律,建立了含间隙平面旋转铰的力学模型;采用几何变形约束法和模态缩聚技术描述柔性机械臂的非线性变形;同时考虑两个旋转铰的间隙特性和柔性臂的弹性变形,最终采用Kane方程建立了含铰链间隙的刚-柔机械臂系统的动力学模型．     

Computed torque control of fully-actuated nondeterministic multibody systems

In this paper we are interested in the dynamic behavior of a slider-crank mechanism with single and two revolute clearance joints. Due to the clearance existence in the revolute joints, it is important to choose an appropriate contact force model in analyzing the dynamic response of a slider-crank mechanism with clearances. The dynamic equations are established by combining the Newton–Euler equations with modified contact force model and improved Coulomb friction force model, and the Baumgarte stabilization approach is used to improve the numerical stability. According to numerical and experimental results, the method of continuous contact can be verified to be reasonable. Comparing dynamic the response between one clearance joint and two clearance joints in a crank-slider mechanism, it is easy to find a significant mutual coupling region due to the presence of two clearance joints by simply contact figures. The dynamic response in a mechanism with two clearance joints is not a simple superposition of that in mechanism with one clearance joint. Therefore, all the joints in a multibody system should be modeled as clearance joints.     

Improved bushing models for general multibody systems and vehicle dynamics

The development and computational implementation, on a multibody dynamics environment, of a constitutive relation to model bushing elements associated with mechanical joints used in the models of road and rail vehicles is presented here. These elements are used to eliminate vibrations in vehicles, due to road irregularities, to allow small misalignment of axes, to reduce noise from the transmission, or to decrease wear of the mechanical joints. Bushings are made of a special rubber, used generally in energy dissipation, which presents a nonlinear viscoelastic relationship between the forces and moments and their corresponding displacements and rotations. In the methodology proposed here a finite element model of the bushing is developed in the framework of the finite element code ABAQUS to obtain the constitutive relations of displacement/rotation versus force/moment for different loading cases. The bushing is modeled in a multibody code as a nonlinear restrain that relates the relative displacements between the bodies connected with the joint reaction forces, and it is represented by a matrix constitutive relation. The basic ingredients of the multibody model are the same vectors and points relations used to define kinematic constraints in any multibody formulation. One particular, and relevant, characteristic of the formulation now presented is its ability to represent standard kinematic joints, clearance, and bushing joints just by defining appropriate constitutive relations. Spherical, revolution, cylindrical, and translational bushing joints are modeled, implemented, and demonstrated through the simulation of two multibody models of a road vehicle, one with perfect kinematic joints for the suspension sub-systems, and other with bushing joints. The tests conducted include an obstacle avoidance maneuver and a vehicle riding over bumps. It is shown that the bushing models for vehicle multibody models proposed here are accurate and computationally efficient so that they can be included in the vehicle models leading reliable simulations.     

Impulsive motion of multibody systems

This article uses the piecewise model and Kane’s method to present a procedure for studying impulsive motion of multibody systems. Impulsive motion occurs when the system is subject to either impulsive forces or impulsive constraints, or when subjected to both simultaneously. The Appellian classification of impulsive constraints and the corresponding equations of impulsive motion of the multibody system are discussed. The governing equations are derived based upon multibody formulation procedures developed by Huston. Constraint impulses associated with finite and impulsive constraints are incorporated into impact dynamical equations through the impulsive Lagrange multipliers. The kinetic energy change of the scleronomic multibody system due to the impact is derived. Newton’s impact law is treated as an impulsive constraint equation to study single-point frictionless collision between two multibody systems. Several examples are used to demonstrate and validate the procedure.     

Recursive formulations for multibody systems with frictional joints based on the interaction between bodies

In practice, the clearances of joints in a great number of mechanical systems are well under control. In these cases, some of the existing methods become unpractical because of the little differences in the order of magnitude between relative movements and computational errors. Assuming that the effects of impacts are negligible, we proved that both locations and forces of contacts in joints can be fully determined by parts of joint reaction forces. Based on this fact, a method particularly suited for multibody systems possessing frictional joints with tiny clearances is presented. In order to improve the efficiency of computation, recursive formulations are proposed based on the interactions between bodies. The proposed recursive formulations can improve the computation of joint reaction forces. With the methodology presented in this paper, not only the motion of bodies in a multibody system but also the details about the contacts in joints, such as forces of contacts and locations of contact points, can be obtained. Even with the assumption of impact free, the instants of possible impacts can be detected without relying upon any ambiguous parameters, as indicated by numerical examples in this paper.     

A contact force solution for non-colliding contact dynamics simulation

Rigid-body impact modeling remains an intensive area of research spurred on by new applications in robotics, biomechanics, and more generally multibody systems. By contrast, the modeling of non-colliding contact dynamics has attracted significantly less attention. The existing approaches to solve non-colliding contact problems include compliant approaches in which the contact force between objects is defined explicitly as a function of local deformation, and complementarity formulations in which unilateral constraints are employed to compute contact interactions (impulses or forces) to enforce the impenetrability of the contacting objects. In this article, the authors develop an alternative approach to solve the non-colliding contact problem for objects of arbitrary geometry in contact at multiple points. Similarly to the complementarity formulation, the solution is based on rigid-body dynamics and enforces contact kinematics constraints at the acceleration level. Differently, it leads to an explicit closed-form solution for the normal forces at the contact points. Integral to the proposed formulation is the treatment of tangential contact forces, in particular the static friction. These friction forces must be calculated as a function of microslip velocity or displacement at the contact point. Numerical results are presented for four test cases: (1) a thin rod sliding down a stationary wedge; (2) a cube pushed off a wedge by an applied force; (3) a cube rotating off the wedge under application of an external moment; and (4) the cube and the wedge both moving under application of a moment. To ascertain validity and correctness, the solutions to frictionless and frictional scenarios obtained with the new formulation are compared to those generated by using a commercial simulation tool MSC ADAMS.     

A computer-oriented method for reducing linearized multibody system equations by incorporating constraints

Consider a spatial multibody system with rigid and elastic bodies. The bodies are linked by rigid interconnections (e.g. revolute joints) causing constraints, as well as by flexible interconnections (e.g. springs) causing applied forces. Small motions of the system with respect to a given nominal configuration can be described by linearized dynamic equations and kinematic constraint equations. We present a computer-oriented procedure which allows to develop a minimum number of these equations. There are three problems. First: algorithmic selection of position coordinates; second: condensation of the dynamic equations; third: evaluation of the constraint forces. To demonstrate the procedure, a closed loop multibody system is used as an example.     

Optimization of a hyper-elastic structure with multibody contact using continuum-based shape design sensitivity analysis

In this paper, a continuum-based shape design sensitivity formulation is presented for a hyper-elastic structure with multibody frictional contact. A nearly incompressible constraint is treated using the pressure projection method that projects a hydrostatic pressure into a lower order space to avoid a volumetric locking. The variational formulation for multibody frictional contact is developed using a penalty method that regularizes the solution of the variational inequality. The material derivative of continuum mechanics is utilized to develop the continuum-based shape design sensitivity analysis for the hyper-elastic constitutive relation and penalized contact formulation. The sensitivity equation is solved at each converged load step using the same tangent stiffness of response analysis due to the path dependency of the sensitivity of the frictional contact problem. A very accurate and efficient sensitivity results are shown through shape optimization of a windshield wiper. Received August 8, 1999