Application of the Absolute Nodal Coordinate Formulation to Large Rotation and Large Deformation Problems

Abstract

There are three basic finite element formulations which are used in multibody dynamics. These are the floating frame of reference approach, the incremental method and the large rotation vector approach. In the floating frame of reference and incremental formulations, the slopes are assumed small in order to define infinitesimal rotations that can be treated and transformed as vectors. This description, however, limits the use of some important elements such as beams and plates in a wide range of large displacement applications. As demonstrated in some recent publications, if infinitesimal rotations are used as nodal coordinates, the use of the finite element incremental formulation in the large reference displacement analysis does not lead to exact modeling of the rigid body inertia when the structures rotate as rigid bodies. In this paper, a simple non-incremental finite element procedure that employs the mathematical definition of the slope and uses it to define the element coordinates instead of the infinitesimal and finite rotations is developed for large rotation and deformation problems. By using this description and by defining the element coordinates in the global system, not only the need for performing coordinate transformation is avoided, but also a simple expression for the inertia forces is obtained. The resulting mass matrix is constant and it is the same matrix that appears in linear structural dynamics. It is demonstrated in this paper that this coordinate description leads to exact modeling of the rigid body inertia when the structures rotate as rigid bodies. Nonetheless, the stiffness matrix becomes nonlinear function even in the case of small displacements. The method presented in this paper differs from previous large rotation vector formulations in the sense that the inertia forces, the kinetic energy, and the strain energy are not expressed in terms of any orientation coordinates, and therefore, the method does not require interpolation of finite rotations. While the use of the formulation is demonstrated using a simple planar beam element, the generalization of the method to other element types and to the three dimensional case is straightforward. Using the finite element procedure presented in this paper, beams and plates can be treated as isoparametric elements.

Dynamic modeling and control of flexible manipulators: a review

“…In 1996, Shabana [14] proposed the absolute node coordinate method based on finite element and continuum mechanics theory. The modeling process does not need small deformation and local coordinate system assumptions, it can truly reflect the dynamic behavior of large deformation [2].…”

Earthquake dynamic response of large flexible multibody systems

“…The inventive methods for earthquake structure response analysis and design are based on precise dynamic modeling and simulation. The up-to-date methods for multibody system dynamics simulation taking into account nonlinear behavior of the large flexible structures could be classify in three basic groups: floating frame of references (Shabana, 2005), finite elements in relative coordinates (Zahariev, 2006); and absolute nodal coordinate formulation (Shabana, 1996;Shabana and Mikkola, 2003). The large structures and mechanical devices like skyscrapers, antenna, wind power generators (WPG), etc., are typical multibody systems compiled of rigid and flexible parts with nonlinear behavior.…”

A Variable-Length Beam Element Incorporating the Effect of Spinning

“…In the late 1990s, Shabana (1996) proposed a new method for modeling the dynamics of flexible multi-body systems. Known as the absolute nodal coordinate formulation (ANCF), this method is based on the finite element method (Bonet and Wood, 1997;Khan et al, 2016) and continuum mechanics (Shabana, 2012;Asemi et al, 2014).…”