Implementation of HHT algorithm for numerical integration of multibody dynamics with holonomic constraints

The Hilber–Hughes–Taylor (HHT) algorithm is implemented in this paper to solve the dynamic problem of flexible multibody system with holonomic constraints. The incremental formulas of governing equations of motion, a set of differential algebraic equations of index 3, are presented and integrated di...

Full description

Saved in:
Bibliographic Details
Published in:Nonlinear dynamics Vol. 80; no. 1-2; pp. 817 - 825
Main Authors: Wang, Jielong, Li, Zhigang
Format: Journal Article
Language:English
Published: Dordrecht Springer Netherlands 01.04.2015
Springer Nature B.V
Subjects:
Algebra
Algorithms
Automotive Engineering
Classical Mechanics
Computer simulation
Control
Damping
Differential equations
Dynamical Systems
Engineering
Equations of motion
Iterative methods
Linear equations
Linearization
Mathematical analysis
Mathematical models
Mechanical Engineering
Multibody systems
Nonlinear dynamics
Nonlinear equations
Nonlinearity
Numerical integration
Original Paper
Parameterization
Parametrization
Robustness (mathematics)
Rotation
Truncation errors
Vibration
DAE of index 3
Multibody dynamics
Matrix-free
HHT
GMRES
ISSN:0924-090X, 1573-269X
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:The Hilber–Hughes–Taylor (HHT) algorithm is implemented in this paper to solve the dynamic problem of flexible multibody system with holonomic constraints. The incremental formulas of governing equations of motion, a set of differential algebraic equations of index 3, are presented and integrated direct by current approach. The parameterization of rotation has been paid a special attention, which makes the formula of HHT algorithm much more complicated because of the nonlinearity of rotations. The associated nonlinear algebraic equations are solved with the application of simplified Newton’s iteration to simplify the linearization of finite rotation. Two options are available for solving the incremental linear equations, one is the direct method based on LU decomposition and the other GMRES-based matrix-free method. Based on the estimated truncation error, the time adaptive scheme is designed to enhance the integrating speed. Numerical simulations prove the modified HHT algorithm behavior as the second-order accuracy with numerical damping. The application of conformal rotation vector improves the robustness of the current approach which makes it possible to integral rotation of arbitrary magnitude.
Bibliography:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 14
content type line 23
ISSN:0924-090X
1573-269X
DOI:10.1007/s11071-015-1908-5