A Large Deformation Planar Finite Element for Pipes Conveying Fluid Based on the Absolute Nodal Coordinate Formulation

2009 ◽  
Vol 4 (3) ◽  
Author(s):  
Michael Stangl ◽  
Johannes Gerstmayr ◽  
Hans Irschik
Keyword(s):  
Finite Element ◽  
Large Deformation ◽  
Absolute Nodal Coordinate Formulation ◽  
Equations Of Motion ◽  
Directional Derivatives ◽  
Absolute Nodal Coordinate ◽  
Lagrange’S Equations ◽  
The Absolute ◽  
Lagrange's Equations ◽  
Conveying Fluid

A novel planar pipe finite element conveying fluid with steady flow, suitable for modeling large deformations in the framework of the Bernoulli–Euler beam theory, is presented. The element is based on a third order planar beam finite element, introduced by Berzeri and Shabana (2000, “Development of Simple Models for the Elastic Forces in the Absolute Nodal Co-Ordinate Formulation,” J. Sound Vib., 235(4), pp. 539–565), applying the absolute nodal coordinate formulation. The equations of motion of the pipe finite element are derived using an extended version of Lagrange’s equations of the second kind taking into account the flow of fluid; in contrast, most derivations in the literature are based on Hamilton’s principle or the Newtonian approaches. The advantage of this element in comparison to classical large deformation beam elements, which are based on rotations, is the direct interpolation of position and directional derivatives, which simplifies the equations of motion considerably. As an advantage, Lagrange’s equations of the second kind offer a convenient connection for introducing fluids into multibody dynamic systems. Standard numerical examples show the convergence of the deformation for increasing number of elements. For a cantilever pipe, the critical flow velocities for increasing number of pipe elements are compared with existing works, based on Euler elastica beams and moving discrete masses. The results show good agreement with the reference solutions applying only a small number of pipe finite elements.

A Large Deformation Finite Element for Pipes Conveying Fluid Based on the Absolute Nodal Coordinate Formulation

2007 ◽  
Author(s):  
Michael Stangl ◽  
Johannes Gerstmayr ◽  
Hans Irschik
Keyword(s):  
Finite Element ◽  
Large Deformation ◽  
Absolute Nodal Coordinate Formulation ◽  
Equations Of Motion ◽  
Directional Derivatives ◽  
Absolute Nodal Coordinate ◽  
Lagrange’S Equations ◽  
The Absolute ◽  
Lagrange's Equations ◽  
Conveying Fluid

A novel pipe finite element conveying fluid, suitable for modeling large deformations in the framework of Bernoulli Euler beam theory, is presented. The element is based on a third order planar beam finite element, introduced by Berzeri and Shabana, on basis of the absolute nodal coordinate formulation. The equations of motion for the pipe-element are derived using an extended version of Lagrange’s equations of the second kind for taking into account the flow of fluids, in contrast to the literature, where most derivations are based on Hamilton’s Principle or Newtonian approaches. The advantage of this element in comparison to classical large deformation beam elements, which are based on rotations, is the direct interpolation of position and directional derivatives, which simplifies the equations of motion considerably. As an advantage Lagrange’s equations of the second kind offer a convenient connection for introducing fluids into multibody dynamic systems. Standard numerical examples show the convergence of the deformation for increasing number of elements. For a cantilever pipe, the critical flow velocities for increasing number of pipe elements are compared to existing works, based on Euler elastica beams and moving discrete masses. The results show good agreements with the reference solutions applying only a small number of pipe finite elements.

Three-Dimensional Solid Brick Element Using Slopes in the Absolute Nodal Coordinate Formulation

2013 ◽  
Vol 9 (2) ◽  
Author(s):  
Alexander Olshevskiy ◽  
Oleg Dmitrochenko ◽  
Chang-Wan Kim
Keyword(s):  
Finite Element ◽  
Finite Elements ◽  
Absolute Nodal Coordinate Formulation ◽  
Mass Matrix ◽  
Equations Of Motion ◽  
Test Problems ◽  
Absolute Nodal Coordinate ◽  
Highly Nonlinear ◽  
Brick Element ◽  
The Absolute

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.

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]

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.

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.

Analysis of Thin Plate Structures Using the Absolute Nodal Coordinate Formulation

2005 ◽  
Vol 219 (4) ◽  
pp. 345-355 ◽  
Author(s):  
K Dufva ◽  
A A Shabana
Keyword(s):  
Finite Element ◽  
Thin Plate ◽  
Absolute Nodal Coordinate Formulation ◽  
Mass Matrix ◽  
Plate Element ◽  
Shell Elements ◽  
Absolute Nodal Coordinate ◽  
Plate And Shell ◽  
Spatial Parameters ◽  
The Absolute

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.

Study on Elastic Forces of the Absolute Nodal Coordinate Formulation for Deformable Beams

1999 ◽  
Author(s):  
Yoshitaka Takahashi ◽  
Nobuyuki Shimizu
Keyword(s):  
Absolute Nodal Coordinate Formulation ◽  
Mass Matrix ◽  
Equations Of Motion ◽  
Nonlinear Function ◽  
Absolute Nodal Coordinates ◽  
Large Rotation ◽  
Inertia Effects ◽  
Absolute Nodal Coordinate ◽  
Elastic Forces ◽  
The Absolute

Abstract There are three basic finite element formulations which are used in the dynamics of flexible beams. These are the floating frame of reference approach, the finite segment method and the large rotation vector approach. Recently, the absolute nodal coordinate formulation was proposed by A.A.Shabana et al. In this procedure, there is no need to transform the element matrices since the equations of motion are defined in terms of absolute nodal coordinates. The mass matrix becomes constant, whereas the stiffness matrix becomes nonlinear function of time, even in case of linear elastic problems. One possible method to avoid such cumbersome of the absolute nodal coordinate formulation in calculating clastic forces is to assume the infinitesimal deformation theory against beams undergoing large rotation. In this paper, a new formulation to calculate the elastic forces and add the rotary inertia effects in the expression of the inertia forces. This formulation is based on the assumption that the deformations within each element remain very small. The expression of the resulting clastic force is simple, and the need for performing coordinate transformation is avoided. As the method assumes that the deformation of the beam from a selected beam axis is very small, a large number of finite elements is required for large deformation problems. However, the formulation has been found to be efficient for large rotation and medium deformation problems. Numerical examples are demonstrated for this formulation by using planar flexible pendulum problems.

A new higher-order locking-free beam element based on the absolute nodal coordinate formulation

2017 ◽  
Vol 232 (19) ◽  
pp. 3410-3423 ◽  
Author(s):  
Haidong Yu ◽  
Chunzhang Zhao ◽  
Bin Zheng ◽  
Hao Wang
Keyword(s):  
Large Deformation ◽  
Cross Section ◽  
Absolute Nodal Coordinate Formulation ◽  
Beam Element ◽  
Higher Order ◽  
Continuity Condition ◽  
Large Rotation ◽  
Absolute Nodal Coordinate ◽  
The Absolute ◽  
The Cross

The beam elements based on the absolute nodal coordinate formulation are widely used in large deformation and large rotation problems. Some of them lead to shear and Poisson locking problems when the continuum mechanics method is employed to deduce the generalized elastic force of the element. To circumvent these locking problems, a new higher-order beam element is proposed that may capture the warping and non-uniform stretching distribution of the cross-section by introducing the trapezoidal cross-section deformation mode and increasing the order of interpolation polynomials in transverse direction. The curvature vectors are chosen as the nodal coordinates of the new element that improve the continuity condition at the element interface. Static and dynamic analyses are conducted to investigate the performance of the new element. Poisson locking phenomena may be eliminated effectively for the new element even when Poisson’s ratio is greater than zero. Meanwhile, the distortion deformation of the cross-section may be described directly. The new element has a better convergence performance compared with the spatial absolute nodal coordinate formulation beam element for that shear locking issue is eliminated. The results also show that the new element fulfills energy conservation and may be applied to the dynamics of both straight and initial curved structures with large deformation.

scholarly journals Finite Element Method-Based Elastic Analysis of Multibody Systems. A Review

Mathematics ◽  
2022 ◽  
Vol 10 (2) ◽  
pp. 257
Author(s):  
Sorin Vlase ◽  
Marin Marin ◽  
Negrean Iuliu
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Mechanical System ◽  
Multibody Systems ◽  
Equations Of Motion ◽  
Alternative Methods ◽  
Engineering Applications ◽  
Lagrange’S Equations ◽  
Lagrange's Equations ◽  
Element Method

This paper presents the main analytical methods, in the context of current developments in the study of complex multibody systems, to obtain evolution equations for a multibody system with deformable elements. The method used for analysis is the finite element method. To write the equations of motion, the most used methods are presented, namely the Lagrange equations method, the Gibbs–Appell equations, Maggi’s formalism and Hamilton’s equations. While the method of Lagrange’s equations is well documented, other methods have only begun to show their potential in recent times, when complex technical applications have revealed some of their advantages. This paper aims to present, in parallel, all these methods, which are more often used together with some of their engineering applications. The main advantages and disadvantages are comparatively presented. For a mechanical system that has certain peculiarities, it is possible that the alternative methods offered by analytical mechanics such as Lagrange’s equations have some advantages. These advantages can lead to computer time savings for concrete engineering applications. All these methods are alternative ways to obtain the equations of motion and response time of the studied systems. The difference between them consists only in the way of describing the systems and the application of the fundamental theorems of mechanics. However, this difference can be used to save time in modeling and analyzing systems, which is important in designing current engineering complex systems. The specifics of the analyzed mechanical system can guide us to use one of the methods presented in order to benefit from the advantages offered.

Sign in / Sign up

Export Citation Format

Share Document