Modeling and Simulating of Multi-Rigid-Body System with Friction-Affected Translational Joints

The objective of this paper is to implement and test the theory presented in a companion paper for the non-smooth dynamics exhibited in a bouncing dimer. Our approach revolves around the use of rigid body dynamics theory combined with constraint equations from the Coulomb's frictional law and the complementarity condition to identify the contact status of each contacting point. A set of impulsive differential equations based on Darboux–Keller shock dynamics is established that can deal with the complex behaviours involved in multiple collisions, such as the frictional effects, the local dissipation of energy at each contact point, and the dispersion of energy among various contact points. The paper will revisit the experimental phenomena found in Dorbolo et al . ( Dorbolo et al . 2005 Phys. Rev. Lett. 95 , 044101), and then present a qualitative analysis based on the theory proposed in part I. The value of the static coefficient of friction between the plate and the dimer is successfully estimated, and found to be responsible for the formation of the drift motion of the bouncing dimer. Plenty of numerical simulations are carried out, and precise agreements are obtained by the comparisons with the experimental results.

Frictional force is a vitally important factor in the application of engineering either in macroscopic contacts or in micro/nanoscale contacts. Understanding of the influencing factors about frictional force is essential for the design of miniaturized devices and the use of minimal friction force. In the paper, dimensional analysis is used to analysis factors relative to frictional force. We show that the frictional force scales with where A is the contact area and N is the normal contact force. An experiment is carried out to verify the new friction law.

The simple formulation of “Coulomb's law of friction” which can be found in most textbooks—the force of friction is proportional to the normal force and does not depend on the contact area and velocity—has little to do with the real work of Coulomb. On the contrary, Coulomb found that Amontons' law, as well as the independence of the coefficient of friction on velocity, normal force, contact area and roughness are only a first, very rough approximation. He differentiated between material couples (e.g. metal‐metal), where Amontons' law is a good approximation, and other (wood on metal or wood on wood), where there are significant deviations from Amontons' law. In all cases, however, the dependencies are relatively weak. In contemporary language, we would say they are of logarithmic character: the geometric and loading parameters have to be changed by several orders of magnitude to achieve a change in the coefficient of friction by a factor of two. Coulomb also gives simple two‐term relations, which empirically summarize these experimental findings.

The main purpose of this paper is to present a new planar dynamic model of the biped robot-walking with the supporting foot slippage and contact-impact forces. In contrast to McGeer’s passive dynamic models, normal forces and frictional forces acting on the feet and ground have been taken into account in the proposed model. The equations of motion and of the passive dynamic biped are obtained by using the Standard Lagrange multiplier method. The dynamics equations are obtained by forward dynamics. The normal forces acting on the feet of the passive biped are described based on the viscoelastic contact model. In turn, the frictional forces and the slippage is solve by the equations of Bengisu law for dry friction. Finally, results are validated with experimental results of the literature.

This paper seeks to illustrate the development of a device for tribotesting and analysing of friction and normal forces in orthogonal cutting. In order to investigate the friction and normal forces, the objective of the device will be focused on the analyzing of the behavior of the tool-chip (rake surface) or tool-workpiece interface (flank surface) under dry metal cutting conditions. The device will simulate an orthogonal cutting process conditions that can be afterward modified for another cutting process basically with different geometries. The testing setup allows the friction coefficient to be estimated using the well-known Coulomb’s friction law. Experimental procedures were conducted to validate the new device. The experiments were carried out using a tungsten carbide P20 as insert and C45 as the workpiece material, and the results showed that the coefficient of friction is dependent on the area of contact between the tool and workpiece.

To proceed on dynamic analysis of folding boom system in its range changing case with high precision and convenience, conditionally decoupled the nonlinear flexible system dynamic equations, obtained motion equation of rigid system and flexible element kinematic equation considering rigid flexible coupling effect, an efficient dynamic analysis method for planer link system was given. Based on Lagrange equation, basic formulations of flexible and rigid multi-body dynamic for open link mechanism were established. By introducing equivalent rigid virtue element, the unified pattern of multi-sectioned rigid body system dynamic equation was obtained. Took the rigid body dynamic analysis as precondition and borrowed the thought of KED method, a high-efficiency elastokinetics method was given for planer link system. Took the range changing case of certain overhead working truck arm as example, analyzed by the elastokinetics method, obtains the dynamic response considering rigid-flexible coupling effect. Results show that the elastokinetics method for planer link system is of high efficiency and accuracy.

A theoretical and computational model has been developed for the nonlinear motion of an inextensible beam undergoing large deflections for cantilevered and free–free boundary conditions. The inextensibility condition was enforced through a Lagrange multiplier which acted as a constraint force. The Rayleigh–Ritz method was used by expanding the deflections and the constraint force in modal series. Lagrange's equations were used to derive the equations of motion of the system, and a fourth-order Runge–Kutta solver was used to solve them. Comparisons for the cantilevered beam were drawn to experimental and computational results previously published and show good agreement for responses to both static and dynamic point forces. Some physical insights into the cantilevered beam response at the first and second resonant modes were obtained. The free–free beam condition was investigated at the first and third resonant modes, and the nonlinearity (primarily inertia) was shown to shift the resonant frequency significantly from the linear natural frequency and lead to hysteresis in both modes.

Rack feeders for the automated operation of high bay rackings are of high practical importance. They are characterized by a horizontally movable carriage supporting a tall and flexible vertical beam structure, on which a cage containing the payload can be positioned in vertical direction. To shorten the transport times by using trajectories with increased maximum acceleration and jerk values, accompanying control measures can be introduced counteracting or avoiding undesired vibrations of the flexible structure. In this contribution, both the control-oriented modeling for an experimental setup of such a flexible rack feeder and the model-based design of a gain-scheduled feedforward and feedback control structure are presented. Whereas, a kinematical model is sufficient for the vertical axis, the horizontal motion of the rack feeder is modeled as a planar elastic multibody system with the cage position as scheduling parameter. For the mathematical description of the bending deflections, a one-dimensional Ritz ansatz is introduced. The tracking control design is performed separately for both the horizontal and the vertical axes using decentralized state-space representations. Remaining model uncertainties are estimated by a disturbance observer. The resulting tracking accuracy of the proposed control concept is shown by measurement results from the experimental setup. Furthermore, these results are compared to those obtained with an alternative control concept from previous work.

A new computer haptics algorithm to be used in general interactive manipulations of deformable virtual objects is presented. In multimodal interactive simulations, haptic feedback computation often comes from contact forces. Subsequently, the fidelity of haptic rendering depends significantly on contact space modeling. Contact and friction laws between deformable models are often simplified in up to date methods. They do not allow a "realistic" rendering of the subtleties of contact space physical phenomena (such as slip and stick effects due to friction or mechanical coupling between contacts). In this paper, we use Signorini's contact law and Coulomb's friction law as a computer haptics basis. Real-time performance is made possible thanks to a linearization of the behavior in the contact space, formulated as the so-called Delassus operator, and iteratively solved by a Gauss-Seidel type algorithm. Dynamic deformation uses corotational global formulation to obtain the Delassus operator in which the mass and stiffness ratio are dissociated from the simulation time step. This last point is crucial to keep stable haptic feedback. This global approach has been packaged, implemented, and tested. Stable and realistic 6D haptic feedback is demonstrated through a clipping task experiment.

In mechanics, one often describes microscopic processes such as those leading to friction between relative interfaces using macroscopic variables (relative velocity, temperature, etc.) in order to avoid models of intangible complexity. As a consequence, such macroscopic models are frequently nonsmooth, a prominent example being the Coulomb law of friction. In many cases, these models are perfectly adequate for engineering purposes. Formally, however, since the Fundamental Theorem of Existence and Uniqueness does not apply to these situations, one generally expects that these models possess forward nonuniqueness of solutions. Consequently, numerical computations of such systems might possibly unknowingly discard certain solutions. In this paper, we try to shed further light on this issue by studying solutions of a simple friction oscillator subject to stiction friction. The stiction law is a simple nonsmooth model of friction that is a modification of Coulomb based on the fundamental observation that the dynamic friction force, when the mass is in motion, is smaller than the static friction force during stick. The resulting piecewise smooth vector field of this discontinuous model does not follow the classical Filippov convention, and the concept of a Filippov solution cannot be used. Furthermore, some Carathéodory solutions, i.e., absolutely continuous solutions satisfying the differential equation in a weaker sense, are nonphysical. Therefore, we introduce the concept of stiction solutions. These are the Carathéodory solutions that are physically relevant, i.e., the ones that follow the stiction law. However, we find that some of the stiction solutions are forward nonunique in subregions of the slip onset. We call these solutions singular, in contrast to the regular stiction solutions that are forward unique. In order to further understanding of the nonunique dynamics, we then introduce a general regularization of the model. This gives a singularly perturbed problem that captures the main features of the original discontinuous problem. Using geometric singular perturbation theory, we identify a repelling slow manifold that separates the forward slipping from the forward sticking solutions, leading to high sensitivity to the initial conditions. On this slow manifold we find canard trajectories that have the physical interpretation of delaying the slip onset. Most interestingly, we find that these new solutions do not correspond to stiction solutions in the piecewise-smooth limit, and are therefore seemingly nonphysical, yet they are robust and appear generically in the class of regularizations we consider. Finally, we show that the regularized problem has a family of periodic orbits interacting with the canards. We observe that this family has a saddle stability and that it connects, in the rigid body limit, the two regular, slip-stick branches of the discontinuous problem, which are otherwise disconnected.

In order to examine the effect of the blank holding force on the friction behavior in the sheet metal forming, a sheet metal friction tester is designed and manufactured, which can measure friction forces in various forming conditions such as lubrication, die roughness, drawing speed, radius of die corner, blank holding force, etc., and the friction tests are performed, in which friction coefficients in various blank holding forces and pulling speeds are calculated using Coulomb's friction law. The friction test reveals that friction coefficient decreases as the blank holding force, the drawing speed and lubricant viscosity increase together or individually.

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.

ABSTRACTThe dynamic responses of cable-driven parallel sinking platform are investigated considering the longitudinal–torsional coupled characteristics of the cable. The dynamic equations of the system are derived through the Lagrange equations, and the geometric matching conditions and the constraint forces at the interfaces between the cables and the platform have also been considered in the equations through Lagrangian multipliers. By neglecting the slack state of the cable, the numerical results of the mathematical model are obtained through Newmark-beta method after transforming the derived differential algebraic equations (DAEs) to a system of ordinary differential equations (ODEs). An ADAMS simulation is carried out, and the ADAMS simulation results agree well with the results obtained through Newmark-beta method. The results have shown that the lateral characteristics of the cables have a great effect on the torsional vibration of the platform, and by arranging the cables correctly will help to reduce the rotation of the platform.

In early design phases of ball typed constant velocity joints (CVJ), which are widely used for the uniform transmission of rotational movement in driveshafts, the knowledge of the interactions between the components is essential. Therefore, a simulation model for a multi-body analysis is developed to consider general, transient operating conditions of the joint in the development process. The model is derived using the example of a plunging joint of the rear prop shaft, which has a straight track geometry and two track angles. On this basis, a two-dimensional, analytical contact determination between balls and tracks can be carried out, which reduces the required computation times compared to a general contact determination. Although rigid bodies are used for numerical analysis, local elastic deformations between the contact bodies are possible on the basis of Hertzian contact theory. For a higher level of detail, frictional forces are also taken into account via the Coulomb friction law, which is adapted for numerical use. To validate the simulation model, the influences of the numerical parameters are discussed and the numerical results are compared with the analytical results of the statically undeflected joint condition. Finally, the created simulation model is utilised to present the numerical joint analysis using two examples of joint kinematics analysis.

Couplings in machines and mechanisms exhibiting backlash and friction phenomena can be modeled as multibody systems with unilateral constraints and Coulomb friction. The structure of the differential-algebraic equations describing the system depends on the state of the constraints. The contact forces occurring at active constraints are taken into account in the equations of motion as Lagrange multipliers. Additionally, the kinematic conditions of all active constraints are formulated on the acceleration level. Contact and friction laws are sufficient conditions for state transitions of active constraints, and are represented by nonsmooth characteristics. Several formulations, like the linear complementarity problem, and two different nonlinear systems of equations are presented together with their solution method. The theory is applied to a mechanical system containing three-dimensional and coupled unilateral constraints with friction.