scholarly journals A concise nodal-based derivation of the floating frame of reference formulation for displacement-based solid finite elements

2019 ◽  
Vol 49 (3) ◽  
pp. 291-313 ◽  
Author(s):  
Andreas Zwölfer ◽  
Johannes Gerstmayr

AbstractThe Floating Frame of Reference Formulation (FFRF) is one of the most widely used methods to analyze flexible multibody systems subjected to large rigid-body motion but small strains and deformations. The FFRF is conventionally derived via a continuum mechanics approach. This tedious and circuitous approach, which still attracts attention among researchers, yields so-called inertia shape integrals. These unhandy volume integrals, arising in the FFRF mass matrix and quadratic velocity vector, depend not only on the degrees of freedom, but also on the finite element shape functions. That is why conventional computer implementations of the FFRF are laborious and error prone; they require access to the algorithmic level of the underlying finite element code or are restricted to a lumped mass approximation. This contribution presents a nodal-based treatment of the FFRF to bypass these integrals. Each flexible body is considered in its spatially discretized state ab initio, wherefore the integrals are replaced by multiplications by a constant finite element mass matrix. Besides that, this approach leads to a simpler and concise but rigorous derivation of the equations of motion. The steps to obtain the inertia-integral-free equations of motion (in 2D and 3D spaces) are presented in a clear and comprehensive way; the final result provides ready-to-implement equations of motion without a lumped mass approximation, in contrast to the conventional formulation.

2020 ◽  
Author(s):  
Andreas Zwölfer ◽  
Johannes Gerstmayr

AbstractIn a recent paper of the authors, a novel nodal-based floating frame of reference formulation (FFRF) for solid finite elements has been proposed. The nodal-based approach bypasses the unhandy inertia shape integrals ab initio, i.e. they neither arise in the derivation nor in the final equations of motion, leading to a surprisingly simple derivation and computer implementation without a lumped mass approximation, which is conventionally employed within commercial multibody codes. However, the nodal-based FFRF has so far been presented without modal reduction, which is usually required for efficient simulations. Hence, the aim of this follow-up paper is to bring the nodal-based FFRF into a suitable form, which allows the incorporation of modal reduction techniques to reduce the overall system size down to the number of modes included in the reduction basis, which further reduces the computational complexity significantly. Moreover, this exhibits a way to calculate the so-called FFRF invariants, which are constant “ingredients” required to set up the FFRF mass matrix and quadratic velocity vector, without integrals and without a lumped mass approximation.


Author(s):  
Andreas Zwölfer ◽  
Johannes Gerstmayr

Abstract The conventional continuum-mechanics-based floating frame of reference formulation involves unhandy so-called inertia-shape-integrals in the equations of motion, which is why, commercial multibody software codes resort to a lumped mass approximation to avoid the evaluation of these integrals in their computer implementations. This paper recaps the conventional continuum mechanics floating frame of reference formulation and addresses its drawbacks by summarizing recent developments of the so-called nodal-based floating frame of reference formulation, which avoids inertia shape integrals ab initio, does not rely on a lumped mass approximation, and exhibits a way to calculate the so-called invariants, which are constant “ingredients” required to set up the equations of motion, in a consistent way.


Author(s):  
D García-Vallejo ◽  
H Sugiyama ◽  
A A Shabana

The fact that incorrect unstable solutions are obtained for linearly elastic models motivates the analytical study presented in this paper. The increase in the number of finite elements only leads to an increase in the critical speed. Crucial in the analysis presented in this paper is the fact that the mass matrix and the form of the elastic forces obtained using the absolute nodal coordinate formulation remain the same under orthogonal coordinate transformation. The absolute nodal coordinate formulation, in contrast to conventional finite element formulations, does account for the effect of the coupling between bending and extension. Based on the analytical results obtained using the absolute nodal coordinate formulation, a new correction is proposed for the finite element floating frame of reference formulation in order to introduce coupling between the axial and bending displacements. In this two-part paper, two- and three-dimensional finite element models are used to study the problem of rotating beams. The models are developed using the absolute nodal coordinate formulation that allows for accurate representation of the axial strain, thereby avoiding the ill-conditioning problem that arises when classical displacement-based finite element formulations are used. In the first part of the paper, the case of linear elasticity is considered and assumptions used in the finite element floating frame of reference formulation are investigated. In the second part of the paper, non-linear elasticity is considered. A rotating helicopter blade is simulated, and the complexity of the motion suggests the inclusion of rotary inertia, shear deformation, and non-linear elastic forces in order to obtain an accurate solution that does not suffer from the instability problem regardless of the number of finite elements used.


Author(s):  
H. Irschik ◽  
M. Nader ◽  
M. Stangl ◽  
H.-G. v. Garssen

Formulations in rotordynamics are usually based on the assumption that the displacements of the bearings of the rotor are small, such that, besides the axial rotation, no large rigid-body motions have to be taken into account. This results in linear equations of motion with gyroscopic terms. When the axial angular speed of a rotor is increased, however, as well as for rapidly changing transient conditions, a non-linear coupling between the large axial rotation and the small rigid body motion induced by the compliance of the bearings and the small elastic deformation of the rotor body itself is to be expected. It is the scope of the present contribution to present a rational strategy for dealing with this situation. First, we present a problem-oriented version of the floating-frame-of-reference formulation (FFRF). We use a co-rotating rigid rotor as reference configuration, which allows using linear modes of the non-rotating elastic rotor as Ritz approximations. The position vector of the origin of a body-fixed coordinate system and three suitable Bryant angles are used as rigid body coordinates, and free elastic modes of the rotor are considered as elastic Ritz approximations. The properties of the latter and their consequences upon simplifying the necessary spatial integrals in the FFRF are addressed in some detail. The free modes are obtained from a Finite Elements pre-processing of the elastic rotor body. The non-linear equations of motion of the rotor are obtained afterwards by means of symbolic computation This formulation leads to a set of relations, in which the rigid-body degrees of freedom need not to be small, and which is integrated using an implicit scheme. Results for a rotor with unbalance forces, accelerated by external forces and having linear visco-elastic bearings are successfully compared to a commercial multi-body dynamics code.


Author(s):  
Markus Dibold ◽  
Johannes Gerstmayr ◽  
Hans Irschik

In the present paper, a comparison of the absolute nodal coordinate formulation (ANCF) and the floating frame of reference formulation (FFRF) is performed for standard static and dynamic problems, both in the small and large deformation regime. Special emphasis is laid on the converged solutions and a comparison to analytical and numerical solutions from the literature. In addition to the work of previous authors, the computational performance of both formulations is studied for the dynamic case, where detailed information is provided concerning the different effects influencing the single parts of the computation time. In case of the ANCF finite element, a planar formulation based on the Bernoulli-Euler theory is utilized, consisting of two position and two slope coordinates in each node only. In the FFRF beam finite element, the displacements are described by the rigid body motion and a small superimposed transverse deflection. The latter is described by means of two static modes for the rotation at the boundary and a user-defined number of eigenmodes of the clamped-clamped beam. In numerical studies, the accuracy and computational costs of the two formulations are compared for a cantilever beam, a pendulum and a slider-crank mechanism. It turns out that both formulations have comparable performance and that the choice of the optimal formulation depends on the problem configuration. Recent claims in the literature that the ANCF would have deficiencies compared to the FFRF thus can be refuted.


Author(s):  
Markus Dibold ◽  
Johannes Gerstmayr ◽  
Hans Irschik

In extension to a former work, a detailed comparison of the absolute nodal coordinate formulation (ANCF) and the floating frame of reference formulation (FFRF) is performed for standard static and dynamic problems, both in the small and large deformation regimes. Special emphasis is laid on converged solutions and on a comparison to analytical and numerical solutions from the literature. In addition to the work of previous authors, the computational performance of both formulations is studied for the dynamic case, where detailed information is provided, concerning the different effects influencing the single parts of the computation time. In case of the ANCF finite element, a planar formulation based on the Bernoulli–Euler theory is utilized, consisting of two position and two slope coordinates in each node only. In the FFRF beam finite element, the displacements are described by the rigid body motion and a small superimposed transverse deflection. The latter is described by means of two static modes for the rotation at the boundary and a user-defined number of eigenmodes of the clamped-clamped beam. In numerical studies, the accuracy and computational costs of the two formulations are compared for a cantilever beam, a pendulum, and a slider-crank mechanism. It turns out that both formulations have comparable performance and that the choice of the optimal formulation depends on the problem configuration. Recent claims in literature that the ANCF would have deficiencies compared with the FFRF thus can be refuted.


2017 ◽  
Vol 88 (3) ◽  
pp. 1953-1968 ◽  
Author(s):  
Grzegorz Orzechowski ◽  
Marko K. Matikainen ◽  
Aki M. Mikkola

Sign in / Sign up

Export Citation Format

Share Document