Implicit-Explicit Time Integration in Multibody Dynamics

Author(s):  
Martin Arnold ◽  
Gerhard Hippmann

Robust and efficient time integration methods in multibody dynamics are tailored to the specific structure of the equations of motion. In the present paper we discuss the combination of explicit methods for non-stiff solution components with implicit methods for stiff solution components and constraints. The methods are successfully used in the dynamical simulation of large-scale multibody system models that display a clear partitioning into stiff and non-stiff subsystems.

Author(s):  
Guang Dong ◽  
Zheng-Dong Ma ◽  
Gregory Hulbert ◽  
Noboru Kikuchi

The topology optimization method is extended for the optimization of geometrically nonlinear, time-dependent multibody dynamics systems undergoing nonlinear responses. In particular, this paper focuses on sensitivity analysis methods for topology optimization of general multibody dynamics systems, which include large displacements and rotations and dynamic loading. The generalized-α method is employed to solve the multibody dynamics system equations of motion. The developed time integration incorporated sensitivity analysis method is based on a linear approximation of two consecutive time steps, such that the generalized-α method is only applied once in the time integration of the equations of motion. This approach significantly reduces the computational costs associated with sensitivity analysis. To show the effectiveness of the developed procedures, topology optimization of a ground structure embedded in a planar multibody dynamics system under dynamic loading is presented.


Author(s):  
Mohammad Rezaiee-Pajand ◽  
S. A. H. Esfehani ◽  
H. Ehsanmanesh

A new family of time integration methods is formulated. The recommended technique is useful and robust for the loads with large variations and the systems with nonlinear damping behavior. It is also applicable for the structures with lots of degrees of freedom, and can handle general nonlinear dynamic systems. By comparing the presented scheme with the fourth-order Runge–Kutta and the Newmark algorithms, it is concluded that the new strategy is more stable. The authors’ formulations have good results on amplitude decay and dispersion error analyses. Moreover, the family orders of accuracy are [Formula: see text] and [Formula: see text] for even and odd values of [Formula: see text], respectively. Findings demonstrate the superiority of the new family compared to explicit and implicit methods and dissipative and non-dissipative algorithms.


Author(s):  
Martin M. Tong

Numerical solution of the dynamics equations of a flexible multibody system as represented by Hamilton’s canonical equations requires that its generalized velocities q˙ be solved from the generalized momenta p. The relation between them is p = J(q)q˙, where J is the system mass matrix and q is the generalized coordinates. This paper presents the dynamics equations for a generic flexible multibody system as represented by p˙ and gives emphasis to a systematic way of constructing the matrix J for solving q˙. The mass matrix is shown to be separable into four submatrices Jrr, Jrf, Jfr and Jff relating the joint momenta and flexible body mementa to the joint coordinate rates and the flexible body deformation coordinate rates. Explicit formulas are given for these submatrices. The equations of motion presented here lend insight to the structure of the flexible multibody dynamics equations. They are also a versatile alternative to the acceleration-based dynamics equations for modeling mechanical systems.


2021 ◽  
Vol 2090 (1) ◽  
pp. 012145
Author(s):  
Ryuma Honda ◽  
Hiroki Suzuki ◽  
Shinsuke Mochizuki

Abstract This study presents the impact of the difference between the implicit and explicit time integration methods on a steady turbulent flow field. In contrast to the explicit time integration method, the implicit time integration method may produce significant kinetic energy conservation error because the widely used spatial difference method for discretizing the governing equations is explicit with respect to time. In this study, the second-order Crank-Nicolson method is used as the implicit time integration method, and the fourth-order Runge-Kutta, second-order Runge-Kutta and second-order Adams-Bashforth methods are used as explicit time integration methods. In the present study, both isotropic and anisotropic steady turbulent fields are analyzed with two values of the Reynolds number. The turbulent kinetic energy in the steady turbulent field is hardly affected by the kinetic energy conservation error. The rms values of static pressure fluctuation are significantly sensitive to the kinetic energy conservation error. These results are examined by varying the time increment value. These results are also discussed by visualizing the large scale turbulent vortex structure.


Author(s):  
Tamer Wasfy

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.


Sign in / Sign up

Export Citation Format

Share Document