Numerical Computation of Differential-Algebraic Equations for Non-Linear Dynamics of Multibody Systems Involving Contact Forces

1999 ◽  
Vol 123 (2) ◽  
pp. 272-281 ◽  
Author(s):  
B. Fox ◽  
L. S. Jennings ◽  
A. Y. Zomaya
Keyword(s):  
Multibody Systems ◽  
Virtual Work ◽  
Equations Of Motion ◽  
Variational Problems ◽  
Differential Algebraic Equations ◽  
Contact Forces ◽  
Algebraic Equations ◽  
Differential Algebraic Equation ◽  
Lagrange Equations ◽  
Constraint Equations

The well known Euler-Lagrange equations of motion for constrained variational problems are derived using the principle of virtual work. These equations are used in the modelling of multibody systems and result in differential-algebraic equations of high index. Here they concern an N-link pendulum, a heavy aircraft towing truck and a heavy off-highway track vehicle. The differential-algebraic equation is cast as an ordinary differential equation through differentiation of the constraint equations. The resulting system is computed using the integration routine LSODAR, the Euler and fourth order Runge-Kutta methods. The difficulty to integrate this system is revealed to be the result of many highly oscillatory forces of large magnitude acting on many bodies simultaneously. Constraint compliance is analyzed for the three different integration methods and the drift of the constraint equations for the three different systems is shown to be influenced by nonlinear contact forces.

Analytical Modelling of the Dynamics of the Four-Bar Mechanism: A Comparison of Some Methods

2018 ◽  
Author(s):  
J. P. Meijaard ◽  
V. van der Wijk
Keyword(s):  
Virtual Work ◽  
Equations Of Motion ◽  
Differential Algebraic Equations ◽  
Force Balance ◽  
Dynamic Force ◽  
Algebraic Equations ◽  
Principle Of Virtual Work ◽  
Principal Vector ◽  
Reaction Forces ◽  
Principal Points

Some thoughts about different ways of formulating the equations of motion of a four-bar mechanism are communicated. Four analytic methods to derive the equations of motion are compared. In the first method, Lagrange’s equations in the traditional form are used, and in a second method, the principle of virtual work is used, which leads to equivalent equations. In the third method, the loop is opened, principal points and a principal vector linkage are introduced, and the equations are formulated in terms of these principal vectors, which leads, with the introduced reaction forces, to a system of differential-algebraic equations. In the fourth method, equivalent masses are introduced, which leads to a simpler system of principal points and principal vectors. By considering the links as pseudorigid bodies that can have a uniform planar dilatation, a compact form of the equations of motion is obtained. The conditions for dynamic force balance become almost trivial. Also the equations for the resulting reaction moment are considered for all four methods.

Trajectory Tracking of Underactuated Multibody Systems

2011 ◽  
Author(s):  
Stefan Reichl ◽  
Wolfgang Steiner
Keyword(s):  
Trajectory Tracking ◽  
Multibody Systems ◽  
Degrees Of Freedom ◽  
Inverse Dynamics ◽  
Feedforward Control ◽  
Equations Of Motion ◽  
Three Dimensional ◽  
Differential Algebraic Equations ◽  
State Variables ◽  
Algebraic Equations

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.

A Computational Investigation of the Disk-Housing Impacts of Accelerating Rotors Supported by Hydrodynamic Bearings

2010 ◽  
Vol 78 (2) ◽  
Author(s):  
Jaroslav Zapoměl ◽  
Petr Ferfecki
Keyword(s):  
Virtual Work ◽  
Equations Of Motion ◽  
Angular Speed ◽  
Contact Forces ◽  
Computational Time ◽  
Radial Clearance ◽  
Lagrange Equations ◽  
Time Step ◽  
Hydrodynamic Bearings ◽  
Nonlinear Force

As the radial clearance between disks and the casing of rotating machines is usually very narrow, excessive lateral vibration of accelerating rotors passing critical speeds can produce impacts between the disks and the housing. The computer modeling method is an important tool for investigating such events. In the developed procedure, the shaft is flexible and the disks are absolutely rigid. The hydrodynamic bearings and the impacts are implemented in the mathematical model by means of nonlinear force couplings. Most of the publications and computer codes from the field of rotor dynamics are referred only in the case when the rotor turns at a constant angular speed and in simple cases of disk-housing impacts. Moreover, if the disks turning at variable speed are investigated, the resulting form of the equations of motion derived by different authors slightly differs and the differences depend on the method used for their derivation. Therefore, particular emphasis in this article is given to the derivation of the motion equations of a continuous rotor turning with variable revolutions to explain the mentioned differences, to develop a computer algorithm enabling the investigation of cases when impacts between an arbitrary number of disks and the stationary part take place, and to analyze the mutual interaction between the impacts and the fluid film bearings. The Hertz theory is applied to determine the contact forces. Calculation of the hydrodynamic forces acting on the bearings is based on solving the Reynolds equation and taking cavitation into account. Lagrange equations of the second kind and the principle of virtual work are used to derive equations of motion. The Runge–Kutta method with an adaptive time step is applied for their solution. The applicability of the developed procedure was tested by computer simulations. The results show that it can be used for the modeling of complex rotor systems and that the short computational time enables carrying out calculations for a number of design and operation parameters.

Comparison of Different Methods to Control Constraints Violation in Forward Multibody Dynamics

2013 ◽  
Author(s):  
Paulo Flores ◽  
Parviz E. Nikravesh
Keyword(s):  
Multibody Dynamics ◽  
Multibody Systems ◽  
Generalized Inverse ◽  
Jacobian Matrix ◽  
Equations Of Motion ◽  
Control Constraints ◽  
Algebraic Equations ◽  
Kinematic Constraint ◽  
Constraint Equations ◽  
Velocity Constraint

The dynamic equations of motion for constrained multibody systems are frequently formulated using the Newton-Euler’s approach, which is augmented with the acceleration constraint equations. This formulation results in the establishment of a mixed set of differential and algebraic equations, which are solved in order to predict the dynamic behavior of general multibody systems. It is known that the standard resolution of the equations of motion is highly prone to constraints violation because the position and velocity constraint equations are not fulfilled. In this work, a general review of the main methods commonly used to control or eliminate the violation of the constraint equations in the context of multibody dynamics formulation is presented and discussed. Furthermore, a general and comprehensive methodology to eliminate the constraints violation at the position and velocity levels is also presented. The basic idea of this approach is to add corrective terms to the position and velocity vectors with the intent to satisfy the corresponding kinematic constraint equations. These corrective terms are evaluated as function of the Moore-Penrose generalized inverse of the Jacobian matrix and of the kinematic constraint equations.

Control of Constrained Systems Described by Lagrangian DAEs

2011 ◽  
Author(s):  
Jason P. Frye ◽  
Brian C. Fabien
Keyword(s):  
Motion Tracking ◽  
Controller Design ◽  
Numerical Technique ◽  
Equations Of Motion ◽  
Differential Algebraic Equations ◽  
Constrained Systems ◽  
Algebraic Equations ◽  
Nonlinear Controller ◽  
Constrained System ◽  
Constraint Equations

In this paper, a nonlinear controller design for constrained systems described by Lagrangian differential algebraic equations (DAEs) is presented. The controller design utilizes the structure introduced by the coordinate splitting formulation, a numerical technique used for integration of DAEs. In this structure, the Lagrange multipliers associated with the constraint equations are eliminated, and the equations of motion are transformed into implicit differential equations. Making use of this, a feedback linearizing controller can be chosen for successful motion tracking of the constrained system. Numerical examples demonstrate the controller design can be successfully applied to fully actuated and underactuated systems.

scholarly journals Dynamical Analysis and Simulation Validation of Incompletely Restrained Cable-Suspended Swinging System Driven by Two Cables

2016 ◽  
Vol 2016 ◽  
pp. 1-11
Author(s):  
Naige Wang ◽  
Guohua Cao ◽  
Zhencai Zhu ◽  
Weihong Peng
Keyword(s):  
Equations Of Motion ◽  
Current Method ◽  
Differential Algebraic Equations ◽  
Dynamical Analysis ◽  
Algebraic Equations ◽  
Lagrange Equations ◽  
Simulation Validation ◽  
Geometric Matching ◽  
Adams Simulation ◽  
Suspension Cable

The flexibility of the suspension multicables and driven length difference between two cables cause the translation and rotation of the platform in the incompletely restrained cable-suspended system driven by two cables (IRCSWs2), which are theoretically investigated in this paper. The suspension cables are spatially discretized using the assumed modes method (AMM) and the equations of motion are derived from Lagrange equations of the first kind. Considering all the geometric matching conditions are approximately linear with external actuator, the differential algebraic equations (DAEs) are transformed to a system of ordinary differential equations (ODEs). Using linear boundary conditions of the suspension cable, the current method can obtain not only the accurate longitudinal displacements of cable and posture of the platform, but also the tension between the platform and cables, and the current method is verified by ADAMS simulation.

Review of Contemporary Approaches for Constraint Enforcement in Multibody Systems

2007 ◽  
Vol 3 (1) ◽  
Author(s):  
Olivier A. Bauchau ◽  
André Laulusa
Keyword(s):  
Multibody Systems ◽  
Equations Of Motion ◽  
Differential Algebraic Equations ◽  
Algebraic Equations ◽  
Constraint Violation ◽  
New Techniques ◽  
Index Reduction ◽  
Reduction Techniques ◽  
Constraint Enforcement ◽  
The Relationship

A hallmark of multibody dynamics is that most formulations involve a number of constraints. Typically, when redundant generalized coordinates are used, equations of motion are simpler to derive but constraint equations are present. Approaches to dealing with high index differential algebraic equations, based on index reduction techniques, are reviewed and discussed. Constraint violation stabilization techniques that have been developed to control constraint drift are also reviewed. These techniques are used in conjunction with algorithms that do not exactly enforce the constraints. Control theory forms the basis for a number of these methods. Penalty based techniques have also been developed, but the augmented Lagrangian formulation presents a more solid theoretical foundation. In contrast to constraint violation stabilization techniques, constraint violation elimination techniques enforce exact satisfaction of the constraints, at least to machine accuracy. Finally, as the finite element method has gained popularity for the solution of multibody systems, new techniques for the enforcement of constraints have been developed in that framework. The goal of this paper is to review the features of these methods, assess their accuracy and efficiency, underline the relationship among the methods, and recommend approaches that seem to perform better than others.

Natural Coordinates in the Optimal Control of Multibody Systems

2011 ◽  
Author(s):  
Peter Betsch ◽  
Ralf Siebert ◽  
Nicolas Sa¨nger
Keyword(s):  
Optimal Control ◽  
Multibody Systems ◽  
Mass Matrix ◽  
Equations Of Motion ◽  
Differential Algebraic Equations ◽  
Open Loop ◽  
Algebraic Equations ◽  
Natural Coordinates ◽  
Joint Forces ◽  
New Energy

The formulation of multibody dynamics in terms of natural coordinates (NCs) leads to equations of motion in the form of differential-algebraic equations (DAEs). A characteristic feature of the natural coordinates approach is a constant mass matrix. The DAEs make possible (i) the systematic assembly of open-loop and closed-loop multibody systems, (ii) the design of state-of-the-art structure-preserving integrators such as energy-momentum or symplectic-momentum schemes, and (iii) the direct link to nonlinear finite element methods. However, the use of NCs in the optimal control of multibody systems presents two major challenges. First, the consistent application of actuating joint-forces becomes an issue since conjugate joint-coordinates are not directly available. Secondly, numerical methods for optimal control with index-3 DAEs are still in their infancy. The talk will address the two aforementioned issues. In particular, a new energy-momentum consistent method for the optimal control of multibody systems in terms of NCs will be presented.

Natural Coordinates in the Optimal Control of Multibody Systems

2011 ◽  
Vol 7 (1) ◽  
Author(s):  
Peter Betsch ◽  
Ralf Siebert ◽  
Nicolas Sänger
Keyword(s):  
Optimal Control ◽  
Multibody Systems ◽  
Mass Matrix ◽  
Equations Of Motion ◽  
Differential Algebraic Equations ◽  
Open Loop ◽  
Algebraic Equations ◽  
Natural Coordinates ◽  
Joint Forces ◽  
New Energy

The formulation of multibody dynamics in terms of natural coordinates (NCs) leads to equations of motion in the form of differential-algebraic equations (DAEs). A characteristic feature of the natural coordinates approach is a constant mass matrix. The DAEs make possible (i) the systematic assembly of open-loop and closed-loop multibody systems, (ii) the design of state-of-the-art structure-preserving integrators such as energy-momentum or symplectic-momentum schemes, and (iii) the direct link to nonlinear finite element methods. However, the use of NCs in the optimal control of multibody systems presents two major challenges. First, the consistent application of actuating joint-forces becomes an issue since conjugate joint-coordinates are not directly available. Second, numerical methods for optimal control with index-3 DAEs are still in their infancy. The talk will address the two aforementioned issues. In particular, a new energy-momentum consistent method for the optimal control of multibody systems in terms of NCs will be presented.

Sign in / Sign up

Export Citation Format

Share Document