Evaluation of Complementarity Techniques for Minimal Coordinate Contact Dynamics

2014 ◽  
Author(s):  
Harshavardhan Mylapilli ◽  
Abhinandan Jain
Keyword(s):  
Complementarity Problem ◽  
Optimization Problems ◽  
Plane Surface ◽  
Equations Of Motion ◽  
Nonlinear Complementarity Problem ◽  
Contact Dynamics ◽  
Test Case ◽  
Coulomb’S Friction ◽  
Space Formulation ◽  
Multi Body

In this article, the non-smooth contact dynamics of multi-body systems is formulated as a complementarity problem. Minimal coordinates operational space formulation is used to derive the dynamics equations of motion. Depending on the approach used for modeling Coulomb’s friction, the complementarity problem can be posed either as a linear or a nonlinear problem. Both formulations are studied in this paper. An exact modeling of the friction cone leads to a nonlinear complementarity problem (NCP) formulation whereas a polyhedral approximation of the friction cone results in a linear complementarity problem (LCP) formulation. These complementarity problems are further recast as non-smooth unconstrained optimization problems, which are solved by employing a class of Levenberg-Marquardt algorithms. The necessary theory detailing these techniques is discussed and five schemes are implemented to solve contact dynamics problems. A simple test case of a sphere moving on a plane surface is used to validate these schemes, while a twelve-link pendulum example is chosen to compare the speed and accuracy of the schemes presented in this paper.

Complementarity Techniques for Minimal Coordinate Contact Dynamics

2016 ◽  
Vol 12 (2) ◽  
Author(s):  
Harshavardhan Mylapilli ◽  
Abhinandan Jain
Keyword(s):  
Complementarity Problem ◽  
Optimization Problems ◽  
Nonlinear Complementarity Problem ◽  
Computational Cost ◽  
Contact Dynamics ◽  
Time Step ◽  
Frictional Contacts ◽  
Computational Costs ◽  
Simulation Results ◽  
Dynamics Problems

In this paper, nonsmooth contact dynamics of articulated rigid multibody systems is formulated as a complementarity problem. Minimal coordinate (MC) formulation is used to derive the dynamic equations of motion as it provides significant computational cost benefits and leads to a smaller-sized complementarity problem when compared with the frequently used redundant coordinate (RC) formulation. Additionally, an operational space (OS) formulation is employed to take advantage of the low-order structure-based recursive algorithms that do not require mass matrix inversion, leading to a further reduction in these computational costs. Based on the accuracy with which Coulomb's friction cone is modeled, the complementarity problem can be posed either as a linear complementarity problem (LCP), where the friction cone is approximated using a polygon, or as a nonlinear complementarity problem (NCP), where the friction cone is modeled exactly. Both formulations are studied in this paper. These complementarity problems are further recast as nonsmooth unconstrained optimization problems, which are solved by employing a class of Levenberg–Marquardt (LM) algorithms. The necessary theory detailing these techniques is discussed and five solvers are implemented to solve contact dynamics problems. A simple test case of a sphere moving on a plane surface is used to validate these solvers for a single contact, whereas a 12-link complex pendulum example is chosen to compare the accuracy of the solvers for the case of multiple simultaneous contacts. The simulation results validate the MC-based NCP formulations developed in this paper. Moreover, we observe that the LCP solvers deliver accuracy comparable to that of the NCP solvers when the friction cone direction vectors in the contact tangent plane are aligned with the sliding contact velocity at each time step. The theory and simulation results show that the NCP approach can be seamlessly recast into an MC OS formulation, thus allowing for accurate modeling of frictional contacts, while at the same time reducing overall computational costs associated with contact and collision dynamics problems in articulated rigid body systems.

Modeling and Simulation of Motorcycle Ride Comfort Based on Bump Road

2010 ◽  
Vol 139-141 ◽  
pp. 2643-2647 ◽  
Author(s):  
Dong Mei Yuan ◽  
Xiao Mei Zheng ◽  
Ying Yang
Keyword(s):  
Degrees Of Freedom ◽  
Ride Comfort ◽  
Lagrange Equation ◽  
Equations Of Motion ◽  
Vertical Displacement ◽  
Dynamics Model ◽  
Body Dynamics ◽  
Space Formulation ◽  
Simulation Results ◽  
Multi Body

Through analyzing the motion when motorcycle runs on the bump road, the 5-DOF multi-body dynamics model of motorcycle is developed, the degrees of freedom include vertical displacement of sprung mass, rotation of sprung mass, vertical displacement of driver, and vertical displacement of front and rear suspension under sprung mass. According to Lagrange Equation, the differential equations of motion and state-space formulation are derived. Then bump road is simulated by triangle bump, and input displacement is programmed by MATLAB. With the input of bump road, motorcycle ride comfort is simulated, and the simulation results are verified by experiment results combined with two channels tire-coupling road simulator. It indicates that the simulation results and experiment results match well; the 5-DOF model has guidance for development of motorcycle ride comfort.

An Overview of Several Formulations for Dry and Lubricated Revolute Joint Clearances in Planar Rigid-Multi-Body Mechanical Systems

2014 ◽  
Author(s):  
Paulo Flores ◽  
Hamid M. Lankarani
Keyword(s):  
Numerical Models ◽  
Equations Of Motion ◽  
Mechanical Systems ◽  
Contact Forces ◽  
Hertzian Contact ◽  
Revolute Joint ◽  
Coulomb’S Friction ◽  
Revolute Joints ◽  
Joint Clearances ◽  
Multi Body

The influence of the revolute joint model on the dynamic response of planar multibody mechanical systems is studied in this work. In the sequel of this process, under the framework of the multibody formalisms, a general methodology for modeling the main kinematic aspects of dry revolute joint clearances is revisited. The numerical models for normal and tangential contact forces developed at the clearance joints are also discussed, which are based on the Hertzian contact theory and dry Coulomb’s friction law, respectively. The fundamental kinematic and dynamic issues of the modeling lubricated revolute joints are presented in this work in order to compare them with the dry revolute joint approach. In a simple manner, the lubrication forces are obtained by integrating the pressure distribution evaluated with the aid of Reynolds’ equation corresponding to the dynamic regime. The intra-joint forces developed for both dry and lubricated cases are evaluated based on the state of variable of the system and subsequently included into the dynamic equations of motion of the multibody system as external generalized forces. The main assumptions and procedures adopted throughout this work are demonstrated through simulations of a planar slider-crank mechanism, which includes dry and lubricated revolute joint with clearance. Finally, some experimental data is also presented and analyzed.

scholarly journals The formation of vortices from a surface of discontinuity

1931 ◽  
Vol 134 (823) ◽  
pp. 170-192 ◽  
Keyword(s):  
Steady State ◽  
Plane Surface ◽  
Equations Of Motion ◽  
Steady Motion ◽  
Small Oscillations ◽  
Approximate Form ◽  
History Of ◽  
State Increase ◽  
Liquid Surfaces

Helmholtz was the first to remark on the instability of those “liquid surfaces” which separate portions of fluid moving with different velocities, and Kelvin, in investigating the influence of wind on waves in water, supposed frictionless, has discussed the conditions under which a plane surface of water becomes unstable. Adopting Kelvin’s method, Rayleigh investigated the instability of a surface of discontinuity. A clear and easily accessible rendering of the discussion is given by Lamb. The above investigations are conducted upon the well-known principle of “small oscillations”—there is a basic steady motion, upon which is superposed a flow, the squares of whose components of velocity can be neglected. This method has the advantage of making the equations of motion linear. If by this method the flow is found to be stable, the equations of motion give the subsequent history of the system, for the small oscillations about the steady state always remain “small.” If, however, the method indicates that the system is unstable, that is, if the deviations from the steady state increase exponentially with the time, the assumption of small motions cannot, after an appropriate interval of time, be applied to the case under consideration, and the equations of motion, in their approximate form, no longer give a picture of the flow. For this reason, which is well known, the investigations of Rayleigh only prove the existence of instability during the initial stages of the motion. It is the object of this note to investigate the form assumed by the surface of discontinuity when the displacements and velocities are no longer small.

On an application of the complex nonlinear complementarity problem

1976 ◽  
Vol 15 (1) ◽  
pp. 141-148 ◽  
Author(s):  
J. Parida ◽  
B. Sahoo
Keyword(s):  
Minimization Problem ◽  
Complementarity Problem ◽  
Complex Space ◽  
Nonlinear Complementarity Problem ◽  
Polyhedral Cone ◽  
Convex Minimization ◽  
Convex Minimization Problem ◽  
Polyhedral Cones ◽  
Existence Of A Solution ◽  
Nonlinear Complementarity

A theorem on the existence of a solution under feasibility assumptions to a convex minimization problem over polyhedral cones in complex space is given by using the fact that the problem of solving a convex minimization program naturally leads to the consideration of the following nonlinear complementarity problem: given g: Cn → Cn, find z such that g(z) ∈ S*, z ∈ S, and Re〈g(z), z〉 = 0, where S is a polyhedral cone and S* its polar.

Sign in / Sign up

Export Citation Format

Share Document