Performance of the Incremental and Non-Incremental Finite Element Formulations in Flexible Multibody Problems

1999 ◽  
Vol 122 (4) ◽  
pp. 498-507 ◽  
Author(s):  
Marcello Campanelli ◽  
Marcello Berzeri ◽  
Ahmed A. Shabana
Keyword(s):  
Finite Element ◽  
Multibody Systems ◽  
Absolute Nodal Coordinate Formulation ◽  
Flexible Multibody Systems ◽  
Large Displacement ◽  
Body Motion ◽  
Absolute Nodal Coordinate ◽  
Flexible Multibody ◽  
The Absolute ◽  
Finite Element Formulations

Many flexible multibody applications are characterized by high inertia forces and motion discontinuities. Because of these characteristics, problems can be encountered when large displacement finite element formulations are used in the simulation of flexible multibody systems. In this investigation, the performance of two different large displacement finite element formulations in the analysis of flexible multibody systems is investigated. These are the incremental corotational procedure proposed in an earlier article (Rankin, C. C., and Brogan, F. A., 1986, ASME J. Pressure Vessel Technol., 108, pp. 165–174) and the non-incremental absolute nodal coordinate formulation recently proposed (Shabana, A. A., 1998, Dynamics of Multibody Systems, 2nd ed., Cambridge University Press, Cambridge). It is demonstrated in this investigation that the limitation resulting from the use of the infinitesmal nodal rotations in the incremental corotational procedure can lead to simulation problems even when simple flexible multibody applications are considered. The absolute nodal coordinate formulation, on the other hand, does not employ infinitesimal or finite rotation coordinates and leads to a constant mass matrix. Despite the fact that the absolute nodal coordinate formulation leads to a non-linear expression for the elastic forces, the results presented in this study, surprisingly, demonstrate that such a formulation is efficient in static problems as compared to the incremental corotational procedure. The excellent performance of the absolute nodal coordinate formulation in static and dynamic problems can be attributed to the fact that such a formulation does not employ rotations and leads to exact representation of the rigid body motion of the finite element. [S1050-0472(00)00604-8]

Incremental Finite Element Formulations and Flexible Multibody Dynamics

1999 ◽  
Author(s):  
Marcello Berzeri ◽  
Marcello Campanelli ◽  
A. A. Shabana
Keyword(s):  
Finite Element ◽  
Multibody Dynamics ◽  
Multibody Systems ◽  
Absolute Nodal Coordinate Formulation ◽  
Flexible Multibody Dynamics ◽  
Flexible Multibody Systems ◽  
Large Displacement ◽  
Absolute Nodal Coordinate ◽  
Flexible Multibody ◽  
Finite Element Formulations

Abstract In this investigation, the performance of two different large displacement finite element formulations in the analysis of flexible multibody systems is investigated. These are the incremental corotational procedure proposed by Rankin and Brogan [8] and the non-incremental absolute nodal coordinate formulation recently proposed [9]. It is demonstrated in this investigation that the limitation resulting from the use of the nodal rotations in the incremental corotational procedure can lead to simulation problems even when very simple flexible multibody applications are considered.

Application of Plasticity Theory and Absolute Nodal Coordinate Formulation to Flexible Multibody System Dynamics

2003 ◽  
Vol 126 (3) ◽  
pp. 478-487 ◽  
Author(s):  
Hiroyuki Sugiyama ◽  
Ahmed A. Shabana
Keyword(s):  
Finite Element ◽  
Absolute Nodal Coordinate Formulation ◽  
Plasticity Theory ◽  
Position Vector ◽  
Element Formulation ◽  
Absolute Nodal Coordinate ◽  
Flexible Multibody ◽  
The Absolute ◽  
Return Algorithm ◽  
Finite Element Formulations

The objective of this investigation is to develop a nonlinear finite element formulation for the elastic-plastic analysis of flexible multibody systems. The Lagrangian plasticity theory based on J2 flow theory is used to account for the effect of plasticity in flexible multibody dynamics. It is demonstrated that the principle of objectivity that is an issue when existing finite element formulations using rate-type constitutive equations are used is automatically satisfied when the stress and strain rate are directly calculated in the Lagrangian descriptions using the absolute nodal coordinate formulation employed in this investigation. This is attributed to the fact that, in the finite element absolute nodal coordinate formulation, the position vector gradients can completely define the state of rotation and deformation within the element. As a consequence, the numerical algorithm used to determine the plastic deformations such as the radial return algorithm becomes much simpler when the absolute nodal coordinate formulation is used as compared to existing finite element formulations that employ incrementally objective algorithms. Several numerical examples are presented in order to demonstrate the use of the formulations presented in the paper.

Use of Plasticity Theory in Flexible Multibody System Dynamics

2003 ◽  
Author(s):  
Hiroyuki Sugiyama ◽  
Ahmed A. Shabana
Keyword(s):  
Finite Element ◽  
Absolute Nodal Coordinate Formulation ◽  
Flexible Multibody Dynamics ◽  
Plasticity Theory ◽  
Position Vector ◽  
Absolute Nodal Coordinate ◽  
Flexible Multibody ◽  
The Absolute ◽  
Return Algorithm ◽  
Finite Element Formulations

The objective of this investigation is to develop a general nonlinear finite deformation formulation for the elastic-plastic analysis of flexible multibody systems. The Lagrangian plasticity theory based on J2 flow theory is used to account for the effect of plasticity in flexible multibody dynamics. In addition, it is demonstrated that the principle of objectivity that is an issue when existing finite element formulations using ratetype constitutive equations are used can be fully satisfied when the stress and strain rate are directly calculated in the Lagrangian descriptions using the absolute nodal coordinate formulation employed in this investigation. This is attributed to the fact that, in the finite element absolute nodal coordinate formulation, the position vector gradients can completely define the state of rotation and deformation within the element. As a consequence, the numerical algorithm used to determine the plastic deformations such as the Radial Return Algorithm becomes much simpler when the absolute nodal coordinate formulation is used as compared to existing finite element formulations that employ incrementally objective algorithms. Several numerical examples are presented in order to demonstrate the use of the formulations presented in the paper.

Elastic Forces and Coordinate Systems in the Finite Element Formulations

1999 ◽  
Author(s):  
Marcello Berzeri ◽  
Marcello Campanelli ◽  
A. A. Shabana
Keyword(s):  
Finite Element ◽  
Absolute Nodal Coordinate Formulation ◽  
Frame Of Reference ◽  
Coordinate Systems ◽  
Absolute Nodal Coordinate ◽  
Floating Frame Of Reference ◽  
Elastic Forces ◽  
Floating Frame ◽  
The Absolute ◽  
Finite Element Formulations

Abstract The equivalence of the elastic forces of finite element formulations used in flexible multibody dynamics is the focus of this investigation. Two conceptually different finite element formulations that lead to exact modeling of the rigid body dynamics will be used. These are the floating frame of reference formulation and the absolute nodal coordinate formulation. It is demonstrated in this study that different element coordinate systems, which are used for the convenience of describing the element deformations in the absolute nodal coordinate formulation, lead to similar results as the element size is reduced. The equivalence of the elastic forces in the absolute nodal coordinate and the floating frame of reference formulations is shown. The result of this analysis clearly demonstrates that the instability observed in high speed rotor analytical models due to the neglect of the geometric centrifugal stiffening is not a problem inherent to a particular finite element formulation but only depends on the beam model that is used. Fourier analysis of the solutions obtained in this investigation also sheds new light on the fundamental problem of the choice of the deformable body coordinate system in the floating frame of reference formulation. A new method is presented and used to obtain a simple expression for the elastic forces in the absolute nodal coordinate formulation. This method, which employs a nonlinear elastic strain-displacement relationship, does not result in an unstable solution when the angular velocity is increased.

Comparison of Finite Element Solutions of Non-Rational B-Spline and ANCF Elements in the Analysis of Flexible Multibody Systems

2011 ◽  
Author(s):  
Hiroki Yamashita ◽  
Hiroyuki Sugiyama
Keyword(s):  
Finite Element ◽  
Multibody Systems ◽  
Absolute Nodal Coordinate Formulation ◽  
Beam Element ◽  
Flexible Multibody Systems ◽  
Numerical Convergence ◽  
Absolute Nodal Coordinate ◽  
B Spline ◽  
Numerical Examples ◽  
The Third

In this investigation, comparison of finite element solutions obtained using the B-spline approach and the absolute nodal coordinate formulation (ANCF) is performed. Furthermore, equivalence of the two formulations with different orders of polynomials and degrees of continuity is demonstrated by several numerical examples. The degree of continuity can be easily controlled in B-spline elements by changing knot multiplicities, while continuity conditions associated with higher order derivatives need to be imposed to achieve C2 and higher continuities in ANCF elements. In order to compare element performances of the third and quartic B-spline and ANCF elements, the three-node quartic ANCF beam element is developed. It is demonstrated in several numerical examples that use of B-spline and ANCF elements with same orders and continuities leads to identical results. Furthermore, effects of polynomial orders and continuities on the accuracy and numerical convergence are demonstrated.

Finite element analysis of the geometric stiffening effect. Part 1: A correction in the floating frame of reference formulation

2005 ◽  
Vol 219 (2) ◽  
pp. 187-202 ◽  
Author(s):  
D García-Vallejo ◽  
H Sugiyama ◽  
A A Shabana
Keyword(s):  
Finite Element ◽  
Absolute Nodal Coordinate Formulation ◽  
Frame Of Reference ◽  
Reference Formulation ◽  
Absolute Nodal Coordinate ◽  
Floating Frame Of Reference ◽  
Elastic Forces ◽  
Floating Frame ◽  
The Absolute ◽  
Finite Element Formulations

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.

Two Simple Triangular Plate Elements Based on the Absolute Nodal Coordinate Formulation

2008 ◽  
Vol 3 (4) ◽  
Author(s):  
Oleg Dmitrochenko ◽  
Aki Mikkola
Keyword(s):  
Finite Element ◽  
Numerical Solutions ◽  
Absolute Nodal Coordinate Formulation ◽  
Plate Theory ◽  
Body Motion ◽  
Absolute Nodal Coordinate ◽  
Triangular Plate ◽  
Plate Elements ◽  
Important Addition ◽  
The Absolute

In this paper, two triangular plate elements based on the absolute nodal coordinate formulation (ANCF) are introduced. Triangular elements employ the Kirchhoff plate theory and can, accordingly, be used in thin plate problems. As usual in ANCF, the introduced elements can exactly describe arbitrary rigid body motion when their mass matrices are constant. Previous plate developments in the absolute nodal coordinate formulation have focused on rectangular elements that are difficult to use when arbitrary meshes need to be described. The elements introduced in this study have overcome this problem and represent an important addition to the absolute nodal coordinate formulation. The two elements introduced are based on Specht’s and Morley’s shape functions, previously used in conventional finite element formulations. The numerical solutions of these elements are compared with previously proposed rectangular finite element and analytical results.

A Numerical Approach to Solving Flexible Multibody Systems Using the Absolute Nodal Coordinate Formulation

1999 ◽  
Author(s):  
R. Y. Yakoub ◽  
A. A. Shabana
Keyword(s):  
Absolute Nodal Coordinate Formulation ◽  
Sparse Matrix ◽  
Mass Matrix ◽  
Numerical Approach ◽  
Absolute Nodal Coordinates ◽  
Absolute Nodal Coordinate ◽  
Velocity Transformation ◽  
Flexible Multibody ◽  
Nodal Coordinates ◽  
The Absolute

Abstract By utilizing the fact that the absolute nodal coordinate formulation leads to a constant mass matrix, a Cholesky decomposition of the mass matrix can be used to obtain a constant velocity transformation matrix. This velocity transformation can be used to express the absolute nodal coordinates in terms of the generalized Cholesky coordinates. In this case, the inertia matrix associated with the Cholesky coordinates is the identity matrix, and therefore, an optimum sparse matrix structure can be obtained for the augmented multibody equations of motions. The implementation of a computer procedure based on the absolute nodal coordinate formulation and Cholesky coordinates is discussed in this paper. A flexible four-bar linkage is presented in this paper in order to demonstrate the use of Cholesky coordinates in the simulation of the small and large deformations in flexible multibody applications. The results obtained from the absolute nodal coordinate formulation are compared to those obtained from the floating frame of reference formulation.

Sign in / Sign up

Export Citation Format

Share Document