Natural Coordinates in the Optimal Control of Multibody Systems

Author(s):  
Peter Betsch ◽  
Ralf Siebert ◽  
Nicolas Sänger

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.

Author(s):  
Peter Betsch ◽  
Ralf Siebert ◽  
Nicolas Sa¨nger

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.


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.


1999 ◽  
Vol 121 (4) ◽  
pp. 594-598 ◽  
Author(s):  
V. Radisavljevic ◽  
H. Baruh

A feedback control law is developed for dynamical systems described by constrained generalized coordinates. For certain complex dynamical systems, it is more desirable to develop the mathematical model using more general coordinates then degrees of freedom which leads to differential-algebraic equations of motion. Research in the last few decades has led to several advances in the treatment and in obtaining the solution of differential-algebraic equations. We take advantage of these advances and introduce the differential-algebraic equations and dependent generalized coordinate formulation to control. A tracking feedback control law is designed based on a pointwise-optimal formulation. The stability of pointwise optimal control law is examined.


Author(s):  
Olivier A. Bauchau ◽  
André Laulusa

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.


1999 ◽  
Vol 123 (2) ◽  
pp. 272-281 ◽  
Author(s):  
B. Fox ◽  
L. S. Jennings ◽  
A. Y. Zomaya

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.


2019 ◽  
Vol 87 (1) ◽  
Author(s):  
X. M. Xu ◽  
J. H. Luo ◽  
Z. G. Wu

Abstract The inertia representation of a constrained system includes the formulation of the kinetic energy and its corresponding mass matrix, given the coordinates of the system. The way to find a proper inertia representation achieving better numerical performance is largely unexplored. This paper extends the modified inertia representation (MIR) to the constrained rigid multibody systems. By using the orthogonal projection, we show the possibility to derive the MIR for many types of non-minimal coordinates. We present examples of the derivation of the MIR for both planar and spatial rigid body systems. Error estimation shows that the MIR is different from the traditional inertia representation (TIR) in that its parameter γ can be used to reduce the kinetic energy error. With preconditioned γ, numerical results show that the MIR consistently presents significantly higher numerical accuracy and faster convergence speed than the TIR for the given variational integrator. The idea of using different inertia representations to improve the numerical performance may go beyond constrained rigid multibody systems to other systems governed by differential algebraic equations.


Author(s):  
Olivier A. Bauchau ◽  
Andre´ Laulusa

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. While the dynamic behavior of constrained systems is well understood, the numerical solution of the resulting equations, potentially of differential-algebraic nature, remains problematic. Many different approaches have been proposed over the years, all presenting advantages and drawbacks: the sheer number and variety of methods that have been proposed indicate the difficulty of the problem. A cursory survey of the literature reveals that the various methods fall within broad categories sharing common theoretical foundations. This paper summarizes the theoretical foundations to the enforcement in constraints in multibody dynamics problems. Next, methods based on the use of Lagrange’s equation of the first kind, which are index-3 differential algebraic equations are reviewed. Methods leading to a minimum set of equations are discussed; in view of the numerical difficulties associated with index-3 approaches, reduction to a minimum set is often performed, leading to a number of practical algorithms using methods developed for ordinary differential equations. Finally, alternative 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 has 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.


1996 ◽  
Vol 2 (2) ◽  
pp. 161-192 ◽  
Author(s):  
M. Wösle ◽  
F. Pfeiffer

In couplings of machines and mechanisms, backlash and friction phenomena are always occurring. Whether stick-slip phenomena take place depends on the structure of such couplings. These processes can be modeled as multibody systems with a time-varying topology. Making use of Lagrange multiplier methods with a mathematical formulation of the contact problem is very efficient for large systems with many constraints. The differential-algebraic equations of a system are transformed into a resolvable mathematical form by means of contact laws. In the following, rigid multibody systems with dependent constraints and planar friction will be considered. For the evaluation of such problems, an iterative algorithm is presented. This method is based on transformations of the kinematic secondary conditions in the form of inequalities to equations. In mathematical sense, these transformations are projections of the constraint forces on convex sets. Ultimately, we have a solvable nonlinear system of equations consisting of the differential equations of motion, the constraint equations and the projections of the constraint forces.


Author(s):  
J. P. Meijaard ◽  
V. van der Wijk

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.


Sign in / Sign up

Export Citation Format

Share Document