Explicit Time Integration of Multibody Systems Modelled With Three Rotation Parameters

2020 ◽  
Author(s):  
Stefan Holzinger ◽  
Johannes Gerstmayr
Keyword(s):  
Lie Group ◽  
Multibody Systems ◽  
Degrees Of Freedom ◽  
Time Integration ◽  
Rigid Bodies ◽  
Rotation Vector ◽  
Body Motion ◽  
Euler Parameters ◽  
Time Step ◽  
Rotation Parameters

Abstract Rigid bodies are an essential part of multibody systems. As there are six degrees of freedom in rigid bodies, it is natural but also precarious to use three parameters for the displacement and three parameters for the rotation parameters — since there is no singularity-free description of spatial rotations based on three rotation parameters. Standard formulations based on three rotation parameters avoid singularities, e.g. by applying reparameterization strategies during the time integration of the rotational kinematic equations. Alternatively, Euler parameters are commonly used to avoid singularities. State of the art approaches use Lie group methods, specifically integrators, to model rigid body motion without the need for the above mentioned solutions. However, the methods so far have been based on additional information, e.g., the rotation matrix, which has to been computed in each step. The latter procedure is thus difficult to be implemented in existing codes that are based on three rotation parameters. In this paper, we use the rotation vector to model large rotations. Whereby Lie group integration methods are used to compute consistent updates for the rotation vector in every time step. The resulting rotation vector update is finite, while the derivative of the rotation vector in the singularity becomes unbounded. The advantages of this method are shown in an example of a gyro. Additionally, the method is applied to a multibody system and the effects of crossing singularities are presented.

scholarly journals Time integration of rigid bodies modelled with three rotation parameters

2021 ◽  
Author(s):  
Stefan Holzinger ◽  
Johannes Gerstmayr
Keyword(s):  
Lie Group ◽  
Attitude Determination ◽  
Time Integration ◽  
Rigid Bodies ◽  
Singular Points ◽  
Rotation Vector ◽  
Euler Angles ◽  
Euler Parameters ◽  
Time Step ◽  
Rotation Parameters

AbstractThree rotation parameters are commonly used in multibody dynamics or in spacecraft attitude determination to represent large spatial rotations. It is well known, however, that the direct time integration of kinematic equations with three rotation parameters is not possible in singular points. In standard formulations based on three rotation parameters, singular points are avoided, for example, by applying reparametrization strategies during the time integration of the kinematic equations. As an alternative, Euler parameters are commonly used to avoid singular points. State-of-the-art approaches use Lie group methods, specifically integrators, to model large rigid body rotations. However, the former methods are based on additional information, e.g. the rotation matrix, which must be computed in each time step. Thus, the latter method is difficult to incorporate into existing codes that are based on three rotation parameters. In this contribution, a novel approach for solving rotational kinematics in terms of three rotation parameters is presented. The proposed approach is illustrated by the example of the rotation vector and the Euler angles. In the proposed approach, Lie group time integration methods are used to compute consistent updates for the rotation vector or the Euler angles in each time step and therefore singular points can be surmounted and the accuracy is higher as compared to the direct time integration of rotation parameters. The proposed update formulas can be easily integrated into existing codes that use either the rotation vector or Euler angles. The advantages of the proposed approach are demonstrated with two numerical examples.

Two Lie Group Formulations for Dynamic Multibody Systems With Large Rotations

2011 ◽  
Author(s):  
Olivier Bru¨ls ◽  
Martin Arnold ◽  
Alberto Cardona
Keyword(s):  
Numerical Solution ◽  
Lie Group ◽  
Multibody Systems ◽  
Equations Of Motion ◽  
Time Integration ◽  
Computational Cost ◽  
Differential Algebraic Equations ◽  
Rigid Body Dynamics ◽  
Time Step ◽  
Iteration Matrix

This paper studies the formulation of the dynamics of multibody systems with large rotation variables and kinematic constraints as differential-algebraic equations on a matrix Lie group. Those equations can then be solved using a Lie group time integration method proposed in a previous work. The general structure of the equations of motion are derived from Hamilton principle in a general and unifying framework. Then, in the case of rigid body dynamics, two particular formulations are developed and compared from the viewpoint of the structure of the equations of motion, of the accuracy of the numerical solution obtained by time integration, and of the computational cost of the iteration matrix involved in the Newton iterations at each time step. In the first formulation, the equations of motion are described on a Lie group defined as the Cartesian product of the group of translations R3 (the Euclidean space) and the group of rotations SO(3) (the special group of 3 by 3 proper orthogonal transformations). In the second formulation, the equations of motion are described on the group of Euclidean transformations SE(3) (the group of 4 by 4 homogeneous transformations). Both formulations lead to a second-order accurate numerical solution. For an academic example, we show that the formulation on SE(3) offers the advantage of an almost constant iteration matrix.

Dynamics of Serial Multibody Systems Using the Decoupled Natural Orthogonal Complement Matrices

1999 ◽  
Vol 66 (4) ◽  
pp. 986-996 ◽  
Author(s):  
S. K. Saha
Keyword(s):  
Multibody Systems ◽  
Degrees Of Freedom ◽  
Inverse Dynamics ◽  
Equations Of Motion ◽  
Kinematic Chain ◽  
Rigid Bodies ◽  
Dynamic Equations ◽  
Orthogonal Complement ◽  
Forward Dynamics ◽  
Velocity Constraints

Constrained dynamic equations of motion of serial multibody systems consisting of rigid bodies in a serial kinematic chain are derived in this paper. First, the Newton-Euler equations of motion of the decoupled rigid bodies of the system at hand are written. Then, with the aid of the decoupled natural orthogonal complement (DeNOC) matrices associated with the velocity constraints of the connecting bodies, the Euler-Lagrange independent equations of motion are derived. The De NOC is essentially the decoupled form of the natural orthogonal complement (NOC) matrix, introduced elsewhere. Whereas the use of the latter provides recursive order n—n being the degrees-of-freedom of the system at hand—inverse dynamics and order n3 forward dynamics algorithms, respectively, the former leads to recursive order n algorithms for both the cases. The order n algorithms are desirable not only for their computational efficiency but also for their numerical stability, particularly, in forward dynamics and simulation, where the system’s accelerations are solved from the dynamic equations of motion and subsequently integrated numerically. The algorithms are illustrated with a three-link three-degrees-of-freedom planar manipulator and a six-degrees-of-freedom Stanford arm.

Modeling Spatial Rigid Multibody Systems Using an Explicit-Time Integration Finite Element Solver and a Penalty Formulation

2004 ◽  
Author(s):  
Tamer Wasfy
Keyword(s):  
Finite Element ◽  
Degrees Of Freedom ◽  
Equations Of Motion ◽  
Time Integration ◽  
Rigid Bodies ◽  
Rotation Matrix ◽  
Explicit Time Integration ◽  
Penalty Formulation ◽  
A New Technique ◽  
Rigid Multibody

A new technique for modeling rigid bodies undergoing spatial motion using an explicit time-integration finite element code is presented. The key elements of the technique are: (a) use of the total rotation matrix relative to the inertial frame to measure the rotation of the rigid bodies; (b) time-integration of the rotational equations of motion in a body fixed (material) frame, with the resulting incremental rotations added to the total rotation matrix; (c) penalty formulation for creating connection points (virtual nodes which do not add extra degrees of freedom) on the rigid-body where joints can be placed. The use of the rotation matrix along with incremental rotation updates circumvents the problem of singularities associated with other types of three and four parameter rotation measures. Benchmark rigid multibody dynamics problems are solved to demonstrate the accuracy of the present technique.

Flexible Multibody Systems With Large Deformations Using Absolute Nodal Coordinates for Isoparametric Solid Brick Elements

2003 ◽  
Author(s):  
Lars Ku¨bler ◽  
Peter Eberhard ◽  
Johannes Geisler
Keyword(s):  
Multibody Systems ◽  
Degrees Of Freedom ◽  
Large Deformations ◽  
Flexible Multibody Systems ◽  
Body Motion ◽  
Absolute Nodal Coordinates ◽  
Shell Elements ◽  
Flexible Multibody ◽  
Absolute Coordinates ◽  
Solid Bodies

In this paper a formulation for flexible Multibody Systems (MBS) is proposed where flexible bodies are described using absolute coordinates for isoparametric brick elements. The use of absolute coordinates allows for large deformations and provides an accurate description of rigid body motion and inertia in the case of large rotations without additional considerations. Further, constant mass matrices are obtained for isoparametric elements. Brick elements are important, e. g. if general solid bodies with low stiffness, i. e. not negligible large deformations, are part of the MBS and cannot be modeled using beam, plate, or shell elements. Since only nodal translational degrees of freedom are used for brick elements additional questions arise. For example, imposing joint constraints for relative rotations between two bodies requires a nodal reference frame at connection points. An approach is proposed to define such a reference system utilizing displacement information of three finite element nodes.

Developing a Numerical Simulation Software for 3D Multibody Systems Based on a Unified Computational Modeling Technique

2009 ◽  
Author(s):  
Shahram Shokouhfar ◽  
Sayyid Mahdi Khorsandijou
Keyword(s):  
Numerical Simulation ◽  
Multibody Systems ◽  
Dynamic Models ◽  
Initial Conditions ◽  
Rigid Bodies ◽  
Simulation Software ◽  
Computational Technique ◽  
Main Task ◽  
Algebraic Equations ◽  
Time Step

This article represents the features and capabilities of a newly developed application namely MASS (Mechanisms Analysis and Simulation Software) and the formulation and techniques therein. MASS is a general C++ application program whose main task is to construct and solve the governing algebraic differential motion equations of 3D multibody systems automatically in matrix forms complying with the computational algorithms required for numerical simulation. Newton-Raphson and SVD methods have been used for kinematical assembling and producing consistent initial conditions. Adaptive time-step Runge-Kutta-Fehlberg numerical integration methods might be used for forward dynamics problems. The governing equations perfectly describe the kinematics and dynamics of multibody systems within which 3D kinematical joints and collisions between rigid bodies might be taken into consideration. The unified computational technique for mathematical modeling of kinematical joints is the most important concept on top of which MASS has been implemented. It has occurred due to the existence of thirteen basic kinematical constraint equations. Each kinematical joint might be defined by a set of algebraic equations being selected from the mentioned basic equations. The unified dynamic models for collisions and impulsive loads have been also achieved using the mentioned technique. Simulation results obtained from MASS have been compared with that of the corresponding software of Working Model ver. 6 and a discussion about the coincidences and differences has been exposed.

Dynamics of Multibody Systems With Known Configuration Changes

1991 ◽  
Vol 58 (1) ◽  
pp. 215-221 ◽  
Author(s):  
J. J. McPhee ◽  
R. N. Dubey
Keyword(s):  
Numerical Simulation ◽  
Computer Program ◽  
Multibody Systems ◽  
Equations Of Motion ◽  
Rigid Bodies ◽  
Rotation Vector ◽  
Dynamic Effects ◽  
Physical Systems ◽  
Dynamics Of Multibody Systems ◽  
Configuration Changes

The equations of motion are derived for a system with inertial properties that are varying in time as a result of known relative motions between the rigid bodies comprising the system. This vector-dyadic formulation has been encoded into a computer program, making use of the conformal rotation vector for the representation of rotations. The numerical simulation of two different physical systems is presented in order to illustrate the dynamic effects of the changing inertial properties, and the usefulness of the encoded formulation for modeling these effects.

Sign in / Sign up

Export Citation Format

Share Document