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.

Modeling of Joints With Clearance in Flexible Multibody Systems

2001 ◽  
Author(s):  
Olivier A. Bauchau ◽  
Jesus Rodriguez
Keyword(s):  
Multibody Systems ◽  
Time Integration ◽  
Unconditional Stability ◽  
Flexible Multibody Systems ◽  
Structural Damping ◽  
Numerical Examples ◽  
Integration Schemes ◽  
Flexible Multibody ◽  
Time Integration Schemes ◽  
Multi Body

Abstract This paper is concerned with the modeling of joints with clearance within the framework of finite element based dynamic analysis of nonlinear, flexible multi-body systems. For actual joints, clearance, lubrication and friction phenomena can significantly affect the dynamic response of the system. In this work, the effects of clearance and lubrication are studied for revolute and spherical joints. The formulation is developed within the framework of energy preserving and decaying time integration schemes that provide unconditional stability for nonlinear, flexible multibody systems. Numerical examples are presented that demonstrate the efficiency and accuracy of the proposed approach. The importance of modeling structural damping and limited driving power are discussed.

scholarly journals Investigation of Air Pocket Behavior in Pipelines Using Rigid Column Model and Contributions of Time Integration Schemes

Water ◽  
2021 ◽  
Vol 13 (6) ◽  
pp. 785
Author(s):  
Arman Rokhzadi ◽  
Musandji Fuamba
Keyword(s):  
Transient Flow ◽  
Time Integration ◽  
Driving Pressure ◽  
Implicit Time ◽  
Numerical Dissipation ◽  
Integration Scheme ◽  
Time Step ◽  
Column Model ◽  
Integration Schemes ◽  
Backward Euler

This paper studies the air pressurization problem caused by a partially pressurized transient flow in a reservoir-pipe system. The purpose of this study is to analyze the performance of the rigid column model in predicting the attenuation of the air pressure distribution. In this regard, an analytic formula for the amplitude and frequency will be derived, in which the influential parameters, particularly, the driving pressure and the air and water lengths, on the damping can be seen. The direct effect of the driving pressure and inverse effect of the product of the air and water lengths on the damping will be numerically examined. In addition, these numerical observations will be examined by solving different test cases and by comparing to available experimental data to show that the rigid column model is able to predict the damping. However, due to simplified assumptions associated with the rigid column model, the energy dissipation, as well as the damping, is underestimated. In this regard, using the backward Euler implicit time integration scheme, instead of the classical fourth order explicit Runge–Kutta scheme, will be proposed so that the numerical dissipation of the backward Euler implicit scheme represents the physical dissipation. In addition, a formula will be derived to calculate the appropriate time step size, by which the dissipation of the heat transfer can be compensated.

Flexible Multibody Dynamics Formulation by Hamilton’s Equations

2010 ◽  
Author(s):  
Martin M. Tong
Keyword(s):  
Multibody Dynamics ◽  
Multibody System ◽  
Mass Matrix ◽  
Equations Of Motion ◽  
Flexible Multibody Dynamics ◽  
Flexible Body ◽  
Flexible Multibody System ◽  
Generalized Coordinates ◽  
Flexible Multibody ◽  
The Matrix

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.

scholarly journals A Semi-Implicit Version of the MPAS-Atmosphere Dynamical Core

2015 ◽  
Vol 143 (9) ◽  
pp. 3838-3855 ◽  
Author(s):  
Steven Sandbach ◽  
John Thuburn ◽  
Danail Vassilev ◽  
Michael G. Duda
Keyword(s):  
Time Integration ◽  
Standard Test ◽  
Newton Iteration ◽  
Atmospheric Modeling ◽  
Implicit Time ◽  
Integration Scheme ◽  
Baroclinic Wave ◽  
Time Step ◽  
Integration Schemes ◽  
Helmholtz Problem

Abstract An important question for atmospheric modeling is the viability of semi-implicit time integration schemes on massively parallel computing architectures. Semi-implicit schemes can provide increased stability and accuracy. However, they require the solution of an elliptic problem at each time step, creating concerns about their parallel efficiency and scalability. Here, a semi-implicit (SI) version of the Model for Prediction Across Scales (MPAS) is developed and compared with the original model version, which uses a split Runge–Kutta (SRK3) time integration scheme. The SI scheme is based on a quasi-Newton iteration toward a Crank–Nicolson scheme. Each Newton iteration requires the solution of a Helmholtz problem; here, the Helmholtz problem is derived, and its solution using a geometric multigrid method is described. On two standard test cases, a midlatitude baroclinic wave and a small-planet nonhydrostatic gravity wave, the SI and SRK3 versions produce almost identical results. On the baroclinic wave test, the SI version can use somewhat larger time steps (about 60%) than the SRK3 version before losing stability. The SI version costs 10%–20% more per step than the SRK3 version, and the weak and strong scalability characteristics of the two versions are very similar for the processor configurations the authors have been able to test (up to 1920 processors). Because of the spatial discretization of the pressure gradient in the lowest model layer, the SI version becomes unstable in the presence of realistic orography. Some further work will be needed to demonstrate the viability of the SI scheme in this case.

Computational strategies for flexible multibody systems

2003 ◽  
Vol 56 (6) ◽  
pp. 553-613 ◽  
Author(s):  
Tamer M Wasfy ◽  
Ahmed K Noor
Keyword(s):  
Multibody Systems ◽  
Reference Frames ◽  
Control Strategies ◽  
Experimental Studies ◽  
Flexible Multibody Dynamics ◽  
Flexible Multibody Systems ◽  
Flexible Multibody ◽  
The Status ◽  
Constraint Modeling ◽  
Corotational Frame

The status and some recent developments in computational modeling of flexible multibody systems are summarized. Discussion focuses on a number of aspects of flexible multibody dynamics including: modeling of the flexible components, constraint modeling, solution techniques, control strategies, coupled problems, design, and experimental studies. The characteristics of the three types of reference frames used in modeling flexible multibody systems, namely, floating frame, corotational frame, and inertial frame, are compared. Future directions of research are identified. These include new applications such as micro- and nano-mechanical systems; techniques and strategies for increasing the fidelity and computational efficiency of the models; and tools that can improve the design process of flexible multibody systems. This review article cites 877 references.

Separate Time Integration Based on the Hilber, Hughes, Taylor Scheme for Flexible Bodies With a Large Number of Modes

2020 ◽  
Vol 15 (6) ◽  
Author(s):  
Wolfgang Witteveen ◽  
Florian Pichler
Keyword(s):  
Time Integration ◽  
Flexible Multibody Dynamics ◽  
Integration Scheme ◽  
Step Size ◽  
Body Contact ◽  
Time Integration Scheme ◽  
Core Idea ◽  
Numerical Time Integration ◽  
Taylor Scheme ◽  
Local Deformations

Abstract In the current development of flexible multibody dynamics, the efficient and accurate consideration of distributed and nonlinear forces is an active area of research. Examples are, forces due to body-body contact or due to elastohydrodynamics (EHD). This leads to many additional modes for representing the local deformations in the areas on which those forces act. Recent publications show that these can be several hundred to several thousand additional modes. A conventional, monolithic numerical time integration scheme would lead to unacceptable computing times. This paper presents a method for an efficient time integration of such systems. The core idea is to treat the equations associated with modes representing local deformations separately. Using the Newmark formulas, a fixed point iteration is proposed for these separated equations, which can always be stabilized with decreasing step size. The concluding examples underline this property, as well as the fact that the proposed method massively outperforms the conventional, monolithic time integration with increasing number of modes.

Sign in / Sign up

Export Citation Format

Share Document