scholarly journals Is There an Optimal Choice of Configuration Space for Lie Group Integration Schemes Applied to Constrained MBS?

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 ◽  
Milan Vrdoljak ◽  
Dario Zlatar

Dynamic simulation procedures of flight vehicle maneuvers need robust and efficient integration methods in order to allow for reliable simulation missions. Derivation of such integration schemes in Lie-group settings is especially efficient since the coordinate-free Lie-group dynamical models operate directly on SO(3) rotational matrices and angular velocities, avoiding local rotation parameters and artificial algebraic constraints as well as kinematical differential equations. In the adopted modeling approach, a state-space of the flight vehicle (modeled as a multi-body system comprising rigid bodies) is modeled as a Lie-group. The numerical algorithm is demonstrated and tested within the framework of the characteristic case study of the aircraft 3D maneuver.


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

The exponential and Cayley maps on SE(3) are the prevailing coordinate maps used in Lie group integration schemes for rigid body and flexible body systems. Such geometric integrators are the Munthe–Kaas and generalized- α schemes, which involve the differential and its directional derivative of the respective coordinate map. Relevant closed form expressions, which were reported over the last two decades, are scattered in the literature, and some are reported without proof. This paper provides a reference summarizing all relevant closed-form relations along with the relevant proofs, including the right-trivialized differential of the exponential and Cayley map and their directional derivatives (resembling the Hessian). The latter gives rise to an implicit generalized- α scheme for rigid/flexible multibody systems in terms of the Cayley map with improved computational efficiency.


2008 ◽  
Vol 2008 ◽  
pp. 1-5 ◽  
Author(s):  
Fu Jingli ◽  
Chen Benyong ◽  
Xie Fengping

This letter focuses on studying a new energy-work relationship numerical integration scheme of nonconservative Hamiltonian systems. The signal-stage, multistage, and parallel composition numerical integration schemes are presented for this system. The high-order energy-work relation scheme of the system is constructed by a parallel connection of n multistage scheme of order 2 which its order of accuracy is 2n. The connection, which is discrete analog of usual case, between the change of energy and work of nonconservative force is obtained for nonconservative Hamiltonian systems.This letter also shows that the more the stages of the schemes are, the less the error rate of the scheme is for nonconservative Hamiltonian systems. Finally, an applied example is discussed to illustrate these results.


Author(s):  
Tsu-te Wu

This paper presents an improved methodology for evaluating the dynamic responses of shipping casks subjected to the sequential HAC impact loads. The methodology utilizes the import technique of the finite-element mesh and the analytical results form one dynamic analysis using explicit numerical integration scheme into another dynamic analysis also using explicit numerical integration scheme. The new methodology presented herein has several advantages over conventional methods. An example problem is analyzed to illustrate the application of the present methodology in evaluating the structural responses of a shipping cask to the sequential HAC loading.


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.


Sign in / Sign up

Export Citation Format

Share Document