Application of LCP Contact Formulation on Mechanical Systems With Nonrigid Coupled Unilateral Contact Points Illustrated With a Simple Model of a Windscreen Wiper

Author(s):  
S. Tatzko ◽  
P. Gro¨nefeld ◽  
M. Wangenheim

Technical applications with unilateral contact are often modeled as multiple coupled discrete mass points. The transition between contact states is often described by transforming the contact formulation into a linear complementary problem (LCP). In case of compliant materials such as rubber, the LCP can be simplified so that no algorithm is needed to solve the equation system. Thus, computational effort can be reduced considerably. In this paper a windscreen wiper lip is modeled as a simple mechanical system with unilateral contact points. The system consists of masses which are coupled by linear springs and dampers. The masses can come into contact with a rigid surface. The equations of motion are derived and transformed into an LCP. The modeling of the coupled, compliant system leads to a simplification of the equation system. Therefore, it can be solved line by line as single independent scalar LCP’s. Also at the transition from separation to contact, when an impact occurs, the contact points can be considered individually. It will be shown, that the coupling can be neglected during the infinitesimal small time of impact. The LCP formulation in combination with simple models of compliant structures therefore yields an effective method for treating multibody systems or discretized continua with several unilateral contact points.

Author(s):  
Andreas Müller ◽  
Shivesh Kumar

AbstractDerivatives of equations of motion (EOM) describing the dynamics of rigid body systems are becoming increasingly relevant for the robotics community and find many applications in design and control of robotic systems. Controlling robots, and multibody systems comprising elastic components in particular, not only requires smooth trajectories but also the time derivatives of the control forces/torques, hence of the EOM. This paper presents the time derivatives of the EOM in closed form up to second-order as an alternative formulation to the existing recursive algorithms for this purpose, which provides a direct insight into the structure of the derivatives. The Lie group formulation for rigid body systems is used giving rise to very compact and easily parameterized equations.


1995 ◽  
Vol 62 (3) ◽  
pp. 685-691 ◽  
Author(s):  
F. Ma ◽  
T. K. Caughey

The coefficients of a linear nonconservative system are arbitrary matrices lacking the usual properties of symmetry and definiteness. Classical modal analysis is extended in this paper so as to apply to systems with nonsymmetric coefficients. The extension utilizes equivalence transformations and does not require conversion of the equations of motion to first-order forms. Compared with the state-space approach, the generalized modal analysis can offer substantial reduction in computational effort and ample physical insight.


Author(s):  
Jesús Cardenal ◽  
Javier Cuadrado ◽  
Eduardo Bayo

Abstract This paper presents a multi-index variable time step method for the integration of the equations of motion of constrained multibody systems in descriptor form. The basis of the method is the augmented Lagrangian formulation with projections in index-3 and index-1. The method takes advantage of the better performance of the index-3 formulation for large time steps and of the stability of the index-1 for low time steps, and automatically switches from one method to the other depending on the required accuracy and values of the time step. The variable time stepping is accomplished through the use of an integral of motion, which in the case of conservative systems becomes the total energy. The error introduced by the numerical integrator in the integral of motion during consecutive time steps provides a good measure of the local integration error, and permits a simple and reliable strategy for varying the time step. Overall, the method is efficient and powerful; it is suitable for stiff and non-stiff systems, robust for all time step sizes, and it works for singular configurations, redundant constraints and topology changes. Also, the constraints in positions, velocities and accelerations are satisfied during the simulation process. The method is robust in the sense that becomes more accurate as the time step size decreases.


Author(s):  
Kai Sorge ◽  
Friedrich Pfeiffer

Abstract A preprocessor providing the inertia data of elastic bodies for multibody algorithms is considered. The desired capability of including stress stiffening terms leads to a high computational effort. The method presented allows a straight forward implementation of the basic formulas by starting with the development of a special purpose symbolic calculator. A rotating elastic ring serves as an example.


Author(s):  
Werner Schiehlen ◽  
Makoto Iwamura

In this paper, we consider the problem to minimize the energy consumption for controlled multibody systems utilizing passive elastic elements for energy storage useful for robot systems in manufacturing. Firstly, based on the linearized equations of motion, we analyze the relationship between the consumed energy and the operating time, and the optimal trajectory using optimal control theory. Then, we verify the analytical solution by comparing with the numerical one computed considering the full nonlinear dynamics. After that we derive a condition for the operating time to be optimal, and propose the optimal design method for springs. Finally, we show the effectiveness of the design method by applying it to a 2DOF manipulator.


Author(s):  
Stefan Reichl ◽  
Wolfgang Steiner

This work presents three different approaches in inverse dynamics for the solution of trajectory tracking problems in underactuated multibody systems. Such systems are characterized by less control inputs than degrees of freedom. The first approach uses an extension of the equations of motion by geometric and control constraints. This results in index-five differential-algebraic equations. A projection method is used to reduce the systems index and the resulting equations are solved numerically. The second method is a flatness-based feedforward control design. Input and state variables can be parameterized by the flat outputs and their time derivatives up to a certain order. The third approach uses an optimal control algorithm which is based on the minimization of a cost functional including system outputs and desired trajectory. It has to be distinguished between direct and indirect methods. These specific methods are applied to an underactuated planar crane and a three-dimensional rotary crane.


Author(s):  
Maren Eriksen Eia ◽  
Elise Mari Vigre ◽  
Thorstein Ravneberg Rykkje

Abstract A Knuckle Boom Crane is a pedestal-mounted, slew-bearing crane with a joint in the middle of the distal arm; i.e. boom. This distal boom articulates at the ‘knuckle (i.e.: joint)’ and that allows it to fold back like a finger. This is an ideal configuration for a crane on a ship where storage space is a premium. This project researches the motion and control of a ship mounted knuckle boom crane to minimize the pendulum motion of a hanging load. To do this, the project leverages the Moving Frame Method (MFM). The MFM draws upon Lie group theory — SO(3) and SE(3) — and Cartan’s Moving Frames. This, together with a compact notation from geometrical physics, makes it possible to extract the equations of motion, expeditiously. The work reported here accounts for the masses and geometry of all components, interactive motor couples and prepares for buoyancy forces and added mass on the ship. The equations of motion are solved numerically using a 4th order Runge Kutta (RK4), while solving for the rotation matrix for the ship using the Cayley-Hamilton theorem and Rodriguez’s formula for each timestep. This work displays the motion on 3D web pages, viewable on mobile devices.


2021 ◽  
Author(s):  
Friedrich Pfeiffer

Abstract Constraints in multibody systems are usually treated by a Lagrange I - method resulting in equations of motion together with the constraint forces. Going from non-minimal coordinates to minimal ones opens the possibility to project the original equations directly to the minimal ones, thus eliminating the constraint forces. The necessary procedure is described, a general example of combined machine-process dynamics discussed and a specific example given. For a n-link robot tracking a path the equations of motion are projected onto this path resulting in quadratic form linear differential equations. They define the space of allowed motion, which is generated by a polygon-system.


Author(s):  
E. Bayo ◽  
J. M. Jimenez

Abstract We investigate in this paper the different approaches that can be derived from the use of the Hamiltonian or canonical equations of motion for constrained mechanical systems with the intention of responding to the question of whether the use of these equations leads to more efficient and stable numerical algorithms than those coming from acceleration based formalisms. In this process, we propose a new penalty based canonical description of the equations of motion of constrained mechanical systems. This technique leads to a reduced set of first order ordinary differential equations in terms of the canonical variables with no Lagrange’s multipliers involved in the equations. This method shows a clear advantage over the previously proposed acceleration based formulation, in terms of numerical efficiency. In addition, we examine the use of the canonical equations based on independent coordinates, and conclude that in this second case the use of the acceleration based formulation is more advantageous than the canonical counterpart.


Sign in / Sign up

Export Citation Format

Share Document