EQUATIONS OF MOTION OF COMPLEX MULTIBODY SYSTEMS USING KINEMATICAL DIFFERENTIALS

1989 ◽  
Vol 13 (4) ◽  
pp. 113-121 ◽  
Author(s):  
M. HILLER ◽  
A. KECSKEMETHY

In complex multibody systems the motion of the bodies may depend on only a few degrees of freedom. For these systems, the equations of motion of minimal order, although more difficult to obtain, give a very efficient formulation. The present paper describes an approach for the automatic generation of these equations, which avoids the use of LAGRANGE-multipliers. By a particular concept, designated “kinematical differentials”, the problem of determining the partial derivatives required to state the equations of motion is reduced to a simple re-evaluation of the kinematics. These cover the solution of the global position, velocity and acceleration problems, i.e. the motion of all bodies is determined for given generalized (independent) coordinates. For their formulation and solution, the multibody system is mapped to a network of nonlinear transformation elements which are connected by linear equations. Each transformation element, designated “kinematical transformer”, corresponds to an independent multibody loop. This mapping of the constraint equations makes it possible to find closed-form solutions to the kinematics for a wide variety of technical applications, and (via kinematical differentials) leads also to an efficient formulation of the dynamics. The equations are derived for holonomic, scleronomic systems, but can also be extended to general nonholonomic systems.

Author(s):  
Daniel Stadlmayr ◽  
Wolfgang Witteveen ◽  
Wolfgang Steiner

Commercial multibody system simulation (MBS) tools commonly use a redundant coordinate formulation as part of their modeling strategy. Such multibody systems subject to holonomic constraints result in second-order d-index three differential algebraic equation (DAE) systems. Due to the redundant formulation and a priori estimation of possible flexible body coordinates, the model size increases rapidly with the number of bodies. Typically, a considerable number of constraint equations (and physical degrees-of-freedom (DOF)) are not necessary for the structure's motion but are necessary for its stability like out-of-plane constraints (and DOFs) in case of pure in-plane motion. We suggest a combination of both, physical DOF and constraint DOF reduction, based on proper orthogonal decomposition (POD) using DOF-type sensitive velocity snapshot matrices. After a brief introduction to the redundant multibody system, a modified flat Galerkin projection and its application to index-reduced systems in combination with POD are presented. The POD basis is then used as an identification tool pointing out reducible constraint equations. The methods are applied to one academic and one high-dimensional practical example. Finally, it can be reported that for the numerical examples provided in this work, more than 90% of the physical DOFs and up to 60% of the constraint equations can be omitted. Detailed results of the numerical examples and a critical discussion conclude the paper.


Author(s):  
A. Avello ◽  
E. Bayo

Abstract When a multibody system reaches a singular position, one or more degrees of freedom appear instantaneously and the jacobian matrix of the constraint equations becomes rank-deficient. The classical kinematic formulation is based on the factorization of the jacobian and, therefore, fails in singular positions. In this paper we develop an efficient method, which uses a penalty and an augmented Lagrangian formulation, and successfully handles singular positions. This formulation automatically copes with redundant incompatible constraints and guarantees the stability of the constraints during numerical integration. Critical numerical examples are shown which corroborate these findings.


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):  
Martin Schulze ◽  
Stefan Dietz ◽  
Bernhard Burgermeister ◽  
Andrey Tuganov ◽  
Holger Lang ◽  
...  

Current challenges in industrial multibody system simulation are often beyond the classical range of application of existing industrial simulation tools. The present paper describes an extension of a recursive order-n multibody system (MBS) formulation to nonlinear models of flexible deformation that are of particular interest in the dynamical simulation of wind turbines. The floating frame of reference representation of flexible bodies is generalized to nonlinear structural models by a straightforward transformation of the equations of motion (EoM). The approach is discussed in detail for the integration of a recently developed discrete Cosserat rod model representing beamlike flexible structures into a general purpose MBS software package. For an efficient static and dynamic simulation, the solvers of the MBS software are adapted to the resulting class of MBS models that are characterized by a large number of degrees of freedom, stiffness, and high frequency components. As a practical example, the run-up of a simplified three-bladed wind turbine is studied where the dynamic deformations of the three blades are calculated by the Cosserat rod model.


1998 ◽  
Vol 65 (3) ◽  
pp. 719-726 ◽  
Author(s):  
S. Djerassi

This paper is the third in a trilogy dealing with simple, nonholonomic systems which, while in motion, change their number of degrees-of-freedom (defined as the number of independent generalized speeds required to describe the motion in question). The first of the trilogy introduced the theory underlying the dynamical equations of motion of such systems. The second dealt with the evaluation of noncontributing forces and of noncontributing impulses during such motion. This paper deals with the linear momentum, angular momentum, and mechanical energy of these systems. Specifically, expressions for changes in these quantities during imposition and removal of constraints are formulated in terms of the associated changes in the generalized speeds.


Author(s):  
Junjie Gu ◽  
Xiaoting Rui ◽  
Jianshu Zhang ◽  
Gangli Chen

Riccati transfer matrix method for multibody systems (RMSTMM) has lower matrix order and better numerical stability than transfer matrix method for multibody systems (MSTMM). In order to make technicians more convenient to apply RMSTMM in practical engineering to improve the computational efficiency of dynamics, in this paper, a linear RMSTMM solver is developed based on the linear RMSTMM theory. A solver input document with good compatibility and extensibility is designed based on extensible markup language (XML); The data structure of multibody system is designed based on object-oriented programming method. The technique of auto selecting the cut hinges of closed-loops of the multibody system is established by introducing the correlation matrix and the dynamic connectivity matrix which depict the connecting state of elements. The automatic generation of the derived tree system by cutting off the closed-loops in the multibody system is realized based on the technique. The automatic regularly numbering of dynamics elements of multibody systems is realized based on the depth first recursive traversal algorithm; Finally, the Riccati transfer matrix recursive technique is implemented based on the regular numbers of dynamics elements of the multibody system. An example is given to verify the effectiveness of the solver which provides a powerful tool for extending the application of RMSTMM in practical engineering.


2014 ◽  
Vol 611 ◽  
pp. 40-45
Author(s):  
Darina Hroncová ◽  
Jozef Filas

The paper describes an algorithm for automatic compilation of equations of motion. Lagrange equations of the second kind and the transformation matrices of basic movements are used by this algorithm. This approach is useful for computer simulation of open kinematic chains with any number of degrees of freedom as well as any combination of bonds.


2013 ◽  
Vol 321-324 ◽  
pp. 1725-1729 ◽  
Author(s):  
Zhong Shuang Wang ◽  
Yang Yang Tao ◽  
Quan Yi Wen

In order to increase the reliability and efficiency of the kineto-static analysis of complex multibody systems, the corresponding vector bond graph procedure is proposed. By the kinematic constraint condition, spatial multibody systems can be modeled by vector bond graph. For the algebraic difficulties brought by differential causality in system automatic kineto-static analysis, the effective decoupling method is proposed, thus the differential causalities in system vector bond graph model can be eliminated. In the case of considering EJS, the unified formulae of driving moment and constraint forces at joints are derived based on vector bond graph, which are easily derived on a computer in a complete form and very suitable for spatial multibody systems. As a result, the automatic kineto-static analysis of spatial multibody system on a computer is realized, its validity is illustrated by the spatial multibody system with three degrees of freedom.


2011 ◽  
Vol 346 ◽  
pp. 324-331
Author(s):  
Wei Jiang ◽  
Xin Luo ◽  
Wen Chuan Jia ◽  
Yuan Tai Hu ◽  
Hong Ping Hu

A new algorithm is presented to calculate the degrees of freedom (DOFs) of spatial complex mechanisms by using the coefficient matrix of the linear constraint equations. A joint constraint matrix is firstly put forward for each kind of joint to formulate linear constraint equations in terms of spatial fine displacements of joint acting point with respect to joint frame. Two kinds of transformation are then proposed to rewrite all the constraint equations in terms of a set of fine displacements of all bodies and it leads to a set of homogeneous linear equations. The rank of the resulting coefficient matrix stands for the number of effective constraints and therefore the DOFs of the mechanism can be easily figured out. The proposed method can be widely used to solve the problem of DOFs for many spatial complex mechanisms, which may not be correctly solved with traditional approaches. Besides, the proposed method is very easy for implementation.


Author(s):  
Gilbert Gede ◽  
Dale L. Peterson ◽  
Angadh S. Nanjangud ◽  
Jason K. Moore ◽  
Mont Hubbard

Symbolic equations of motion (EOMs) for multibody systems are desirable for simulation, stability analyses, control system design, and parameter studies. Despite this, the majority of engineering software designed to analyze multibody systems are numeric in nature (or present a purely numeric user interface). To our knowledge, none of the existing software packages are 1) fully symbolic, 2) open source, and 3) implemented in a popular, general, purpose high level programming language. In response, we extended SymPy (an existing computer algebra system implemented in Python) with functionality for derivation of symbolic EOMs for constrained multibody systems with many degrees of freedom. We present the design and implementation of the software and cover the basic usage and workflow for solving and analyzing problems. The intended audience is the academic research community, graduate and advanced undergraduate students, and those in industry analyzing multibody systems. We demonstrate the software by deriving the EOMs of a N-link pendulum, show its capabilities for LATEX output, and how it integrates with other Python scientific libraries — allowing for numerical simulation, publication quality plotting, animation, and online notebooks designed for sharing results. This software fills a unique role in dynamics and is attractive to academics and industry because of its BSD open source license which permits open source or commercial use of the code.


Sign in / Sign up

Export Citation Format

Share Document