Singularity-Free Lie Group Integration and Geometrically Consistent Evaluation of Multibody System Models described in terms of Standard Absolute Coordinates

Author(s):  
Andreas Mueller

Abstract A classical approach to the MBS modeling is to use absolute coordinates, i.e. a set of (possibly redundant) coordinates that describe the absolute position and orientation of the individual bodies w.r.t. to an inertial frame (IFR). A well-known problem for the time integration of the equations of motion (EOM) is the lack of a singularity-free parameterization of spatial motions, which is usually tackled by using unit quaternions. Lie group integration methods were proposed as alternative approach to the singularity-free time integration. Lie group integration methods, operating directly on the configuration space Lie group, are incompatible with standard formulations of the EOM, and cannot be implemented in existing MBS simulation codes without a major restructuring. A framework for interfacing Lie group integrators to standard EOM formulations is presented in this paper. It allows describing MBS in terms of various absolute coordinates and at the same using Lie group integration schemes. The direct product group SO(3)xR3; and the semidirect product group SE(3) are use for representing rigid body motions. The key element of this method is the local-global transitions (LGT) transition map, which facilitates the update of (global) absolute coordinates in terms of the (local) coordinates on the Lie group. This LGT map is specific to the absolute coordinates, the local coordinates on the Lie group, and the Lie group used to represent rigid body configurations. This embedding of Lie group integration methods allows for interfacing with standard vector space integration methods.

2021 ◽  
Author(s):  
Andreas Müller

Abstract A classical approach to the modeling of multibody systems (MBS) is to use absolute coordinates, i.e. a set of (redundant) coordinates that describe the absolute position and orientation of the individual bodies w.r.t. to an inertial frame (IFR). A well-known problem for the time integration of the equations of motion (EOM) is the lack of a singularity-free parameterization of spatial motions, which is usually tackled by using unit quaternions. Lie group integration methods were proposed as alternative approach to the singularity-free time integration. Lie group methods are inherently coordinate free and thus incompatible with any absolute coordinate description. In this paper, an integration scheme is proposed that allows describing MBS in terms of arbitrary absolute coordinates and at the same using Lie group integration schemes, which allows for singularity-free time integration. Moreover, the direct product group SO (3) × ℝ3 as well as the semidirect product group SE (3) can be use for representing rigid body motions, which is beneficial for constraint satisfaction. The crucial step of this method, which renders the underlying Lie group integration scheme applicable to EOM in absolute coordinates, is the update of the (global) absolute coordinates in terms of the (local) coordinates on the Lie group by means of a local-global transitions (LGT) transition map. This LGT map depends on the used absolute coordinates and the local coordinates on the Lie group, but also on the Lie group itself used to represent rigid body configurations (respectively the deformation field of flexible bodies), i.e. the geometry of spatial frame motions. The Lie group formulation is thus embedded, which allows interfacing with standard vector space integration methods.


Author(s):  
Andreas Müller

Geometric methods have become increasingly accepted in computational multibody system (MBS) dynamics. This includes the kinematic and dynamic modeling as well as the time integration of the equations of motion. In particular, the observation that rigid body motions form a Lie group motivated the application of Lie group integration schemes, such as the Munthe-Kaas method. Also established vector space integration schemes tailored for structural and MBS dynamics were adopted to the Lie group setting, such as the generalized α integration method. Common to all is the use of coordinate mappings on the Lie group SE(3) of Euclidean motions. In terms of canonical coordinates (screw coordinates), this is the exponential mapping. Rigid body velocities (twists) are determined by its right-trivialized differential, denoted dexp. These concepts have, however, not yet been discussed in compact and concise form, which is the contribution of this paper with particular focus on the computational aspects. Rigid body motions can also be represented by dual quaternions, that form the Lie group Sp̂(1), and the corresponding dynamics formulations have recently found a renewed attention. The relevant coordinate mappings for dual quaternions are presented and related to the SE(3) representation. This relation gives rise to a novel closed form of the dexp mapping on SE(3). In addition to the canonical parameterization via the exponential mapping, the noncanonical parameterization via the Cayley mapping is presented.


Author(s):  
Andreas Müller ◽  
Zdravko Terze

Recently various numerical integration schemes have been proposed for numerically simulating the dynamics of constrained multibody systems (MBS) operating. These integration schemes operate directly on the MBS configuration space considered as a Lie group. For discrete spatial mechanical systems there are two Lie group that can be used as configuration space: SE(3) and SO(3) × ℝ3. Since the performance of the numerical integration scheme clearly depends on the underlying configuration space it is important to analyze the effect of using either variant. For constrained MBS a crucial aspect is the constraint satisfaction. In this paper the constraint violation observed for the two variants are investigated. It is concluded that the SE(3) formulation outperforms the SO(3) × ℝ3 formulation if the absolute motions of the rigid bodies, as part of a constrained MBS, belong to a motion subgroup. In all other cases both formulations are equivalent. In the latter cases the SO(3) × ℝ3 formulation should be used since the SE(3) formulation is numerically more complex, however.


Author(s):  
Zdravko Terze ◽  
Andreas Müller ◽  
Dario Zlatar

The paper presents two novel second order conservative Lie-group geometric methods for integration of rigid body rotational dynamics. First proposed algorithm is a fully explicit scheme that exactly conserves spatial angular momentum of a free spinning body. The method is inspired by the Störmer–Verlet integration algorithm for solving ordinary differential equations (ODEs), which is also momentum conservative when dealing with ODEs in linear spaces but loses its conservative properties in a nonlinear regime, such as nonlinear SO(3) rotational group. Then, we proposed an algorithm that is an implicit integration scheme with a direct update in SO(3). The method is algorithmically designed to conserve exactly both of the two “main” motion integrals of a rotational rigid body, i.e., spatial angular momentum of a torque-free body as well as its kinetic energy. As it is shown in the paper, both methods also preserve Lagrangian top integrals of motion in a very good manner, and generally better than some of the most successful conservative schemes to which the proposed methods were compared within the presented numerical examples. The proposed schemes can be easily applied within the integration algorithms of the dynamics of general rigid body systems.


Sign in / Sign up

Export Citation Format

Share Document