scholarly journals The significance of the configuration space Lie group for the constraint satisfaction in numerical time integration of multibody systems

2014 ◽  
Vol 82 ◽  
pp. 173-202 ◽  
Author(s):  
Andreas Müller ◽  
Zdravko Terze
Keyword(s):  
Configuration Space ◽  
Lie Group ◽  
Constraint Satisfaction ◽  
Multibody Systems ◽  
Time Integration ◽  
Numerical Time Integration

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.

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.

Direct Sensitivity Analysis of Multibody Systems: A Vehicle Dynamics Benchmark

2019 ◽  
Vol 14 (2) ◽  
Author(s):  
Alfonso Callejo ◽  
Daniel Dopico
Keyword(s):  
Sensitivity Analysis ◽  
Multibody Systems ◽  
Automatic Differentiation ◽  
Time Integration ◽  
General Purpose ◽  
Contact Forces ◽  
Design Parameters ◽  
Computational Techniques ◽  
Computationally Efficient ◽  
Integration Algorithm

Algorithms for the sensitivity analysis of multibody systems are quickly maturing as computational and software resources grow. Indeed, the area has made substantial progress since the first academic methods and examples were developed. Today, sensitivity analysis tools aimed at gradient-based design optimization are required to be as computationally efficient and scalable as possible. This paper presents extensive verification of one of the most popular sensitivity analysis techniques, namely the direct differentiation method (DDM). Usage of such method is recommended when the number of design parameters relative to the number of outputs is small and when the time integration algorithm is sensitive to accumulation errors. Verification is hereby accomplished through two radically different computational techniques, namely manual differentiation and automatic differentiation, which are used to compute the necessary partial derivatives. Experiments are conducted on an 18-degree-of-freedom, 366-dependent-coordinate bus model with realistic geometry and tire contact forces, which constitutes an unusually large system within general-purpose sensitivity analysis of multibody systems. The results are in good agreement; the manual technique provides shorter runtimes, whereas the automatic differentiation technique is easier to implement. The presented results highlight the potential of manual and automatic differentiation approaches within general-purpose simulation packages, and the importance of formulation benchmarking.

DQFEM Analyses of Static and Dynamic Nonlinear Elastic-Plastic Problems Using a GSR-Based Accelerated Constant Stiffness Equilibrium Iteration Technique

2001 ◽  
Vol 123 (3) ◽  
pp. 310-317 ◽  
Author(s):  
Chang-New Chen
Keyword(s):  
Finite Element ◽  
Numerical Technique ◽  
Time Integration ◽  
Differential Quadrature ◽  
Element Discretization ◽  
Constant Stiffness ◽  
Structural Problems ◽  
Nonlinear Structural ◽  
Numerical Time Integration ◽  
Nonlinear Elastic

An integrated numerical technique for static and dynamic nonlinear structural problems adopting the equilibrium iteration is proposed. The differential quadrature finite element method (DQFEM), which uses the differential quadrature (DQ) techniques to the finite element discretization, is used to analyze the static and dynamic nonlinear structural mechanics problems. Numerical time integration in conjunction with the use of equilibrium iteration is used to update the response history. The equilibrium iteration can be carried out by the accelerated iteration schemes. The global secant relaxation-based accelerated constant stiffness and diagonal stiffness-based predictor-corrector equilibrium iterations which are efficient and reliable are used for the numerical computations. Sample problems are analyzed. Numerical results demonstrate the algorithm.

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

2021 ◽  
Author(s):  
Andreas Mueller
Keyword(s):  
Rigid Body ◽  
Lie Group ◽  
Time Integration ◽  
Integration Methods ◽  
Product Group ◽  
Local Coordinates ◽  
Group Integration ◽  
The Absolute ◽  
Lie Group Integration ◽  
Absolute Coordinates

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.

Energy-Momentum Conserving Time Integration of Modally Reduced Flexible Multibody Systems

2013 ◽  
Author(s):  
Alexander Humer ◽  
Johannes Gerstmayr
Keyword(s):  
Angular Momentum ◽  
Multibody Systems ◽  
Mass Matrix ◽  
Time Integration ◽  
Flexible Multibody Dynamics ◽  
Integration Scheme ◽  
Integration Schemes ◽  
Flexible Multibody ◽  
Floating Frame ◽  
Energy And Momentum

Many conventional time integration schemes frequently adopted in flexible multibody dynamics fail to retain the fundamental conservation laws of energy and momentum of the continuous time domain. Lack of conservation, however, in particular of angular momentum, may give rise to unexpected, unphysical results. To avoid such problems, a scheme for the consistent integration of modally reduced multibody systems subjected to holonomic constraints is developed in the present paper. As opposed to the conventional approach, in which the floating frame of reference formulation is combined with component mode synthesis for approximating the flexible deformation, an alternative, recently proposed formulation based on absolute coordinates is adopted in the analysis. Owing to the linear relationship between the generalized coordinates and the absolute displacement, the inertia terms in the equations of motion attain a very simple structure. The mass matrix remains independent of the current state of deformation and the velocity dependent term known from the floating frame approach vanishes due to the absence of relative coordinates. These advantageous properties facilitate the construction of an energy and momentum consistent integration scheme. By the mid-point rule, algorithmic conservation of both linear and angular momentum is achieved. In order to consistently integrate the total energy of the system, the discrete derivative needs to be adopted when evaluating the strain energy gradient and the derivative of the algebraic constraint equations.

Sign in / Sign up

Export Citation Format

Share Document