Modeling Method and Application of Rational Finite Element Based on Absolute Nodal Coordinate Formulation

2018 ◽  
Vol 31 (2) ◽  
pp. 207-228 ◽  
Author(s):  
Chao Ma ◽  
Cheng Wei ◽  
Jing Sun ◽  
Bin Liu
1999 ◽  
Vol 122 (4) ◽  
pp. 498-507 ◽  
Author(s):  
Marcello Campanelli ◽  
Marcello Berzeri ◽  
Ahmed A. Shabana

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]


Author(s):  
Marcello Berzeri ◽  
Marcello Campanelli ◽  
A. A. Shabana

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.


Author(s):  
Alexander Olshevskiy ◽  
Oleg Dmitrochenko ◽  
Chang-Wan Kim

The present paper contributes to the field of flexible multibody systems dynamics. Two new solid finite elements employing the absolute nodal coordinate formulation are presented. In this formulation, the equations of motion contain a constant mass matrix and a vector of generalized gravity forces, but the vector of elastic forces is highly nonlinear. The proposed solid eight node brick element with 96 degrees of freedom uses translations of nodes and finite slopes as sets of nodal coordinates. The displacement field is interpolated using incomplete cubic polynomials providing the absence of shear locking effect. The use of finite slopes describes the deformed shape of the finite element more exactly and, therefore, minimizes the number of finite elements required for accurate simulations. Accuracy and convergence of the finite element is demonstrated in nonlinear test problems of statics and dynamics.


2003 ◽  
Vol 126 (3) ◽  
pp. 478-487 ◽  
Author(s):  
Hiroyuki Sugiyama ◽  
Ahmed A. Shabana

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.


Author(s):  
K Dufva ◽  
A A Shabana

The absolute nodal coordinate formulation can be used in multibody system applications where the rotation and deformation within the finite element are large and where there is a need to account for geometrical non-linearities. In this formulation, the gradients of the global positions are used as nodal coordinates and no rotations are interpolated over the finite element. For thin plate and shell elements, the plane stress conditions can be applied and only gradients obtained by differentiation with respect to the element mid-surface spatial parameters need to be defined. This automatically reduces the number of element degrees of freedoms, eliminates the high frequencies due to the oscillations of some gradient components along the element thickness, and as a result makes the plate element computationally more efficient. In this paper, the performance of a thin plate element based on the absolute nodal coordinate formulation is investigated. The lower dimension plate element used in this investigation allows for an arbitrary rigid body displacement and large deformation within the element. The element leads to a constant mass matrix and zero Coriolis and centrifugal forces. The performance of the element is compared with other plate elements previously developed using the absolute nodal coordinate formulation. It is shown that the finite element used in this investigation is much more efficient when compared with previously proposed elements in the case of thin structures. Numerical examples are presented in order to demonstrate the use of the formulation developed in this paper and the computational advantages gained from using the thin plate element. The thin plate element examined in this study can be efficiently used in many applications including modelling of paper materials, belt drives, rotor dynamics, and tyres.


Author(s):  
Marcello Berzeri ◽  
Marcello Campanelli ◽  
A. A. Shabana

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.


2006 ◽  
Vol 298 (4-5) ◽  
pp. 1129-1149 ◽  
Author(s):  
Hiroyuki Sugiyama ◽  
Johannes Gerstmayr ◽  
Ahmed A. Shabana

Author(s):  
Bassam A. Hussein ◽  
Hiroyuki Sugiyama ◽  
Ahmed A. Shabana

The finite element absolute nodal coordinate formulation (ANCF) leads to beam and plate models that relax the assumption of the classical Euler-bernoulli and Timoshenko beam and Mindlin plate theories. In these more general models, the cross section is allowed to deform and it is no longer treated as a rigid surface. The coupling between the bending and cross section deformations leads to the new ANCF-coupled deformation modes that are examined in this study. While these coupled deformation can be source of numerical and convergence problems when thin and stiff beam models are considered, the inclusion of the effect of these modes in the dynamic model is necessary in the case of very flexible structures. In order to examine the effect of these coupled deformation modes in this investigation, three different large deformation dynamic beam models are discussed. Two of these models, which differ in the way the beam elastic forces are calculated in the absolute nodal coordinate formulation, allow for systematically eliminating the coupled deformation modes, while the third allows for including these modes. The first of these models is based on a general continuum mechanics approach that leads to a model that includes the ANCF-coupled deformation modes; while the second and third methods that can be used to eliminate the coupled deformation modes are based on the elastic line approach and the Hellinger-Reissner principle. It is shown in this study that the inclusion of the ANCF-coupled deformation modes introduces geometric stiffening effects that can not be captured using other finite element models.


Sign in / Sign up

Export Citation Format

Share Document