Ahmed A. Shabana

Three Dimensional Absolute Nodal Coordinate Formulation for Beam Elements: Theory

The description of a beam element by only the displacement of its centerline leads to some difficulties in the representation of the torsion and shear effects. For instance such a representation does not capture the rotation of the beam as a rigid body about its own axis. This problem was circumvented in the literature by using a local coordinate system in the incremental finite element method or by using the multibody floating frame of reference formulation. The use of such a local element coordinate system leads to a highly nonlinear expression for the inertia forces as the result of the large element rotation. In this investigation, an absolute nodal coordinate formulation is presented for the large rotation and deformation analysis of three dimensional beam elements. This formulation leads to a constant mass matrix, and as a result, the vectors of the centrifugal and Coriolis forces are identically equal to zero. The formulation presented in this paper takes into account the effect of rotary inertia, torsion and shear, and ensures continuity of the slopes as well as the rotation of the beam cross section at the nodal points. Using the proposed formulation curved beams can be systematically modeled.

A non-incremental finite element procedure for the analysis of large deformation of plates and shells in mechanical system applications

In this investigation, a non-incremental solution procedure for the finite rotationand large deformation analysis of plates is presented. The method, whichis based on the absolute nodal coordinate formulation, leads to plateelements capable of representing exact rigid body motion. In thismethod, continuity conditions on all the displacement gradients areimposed. Therefore, non-smoothness of the plate mid-surface at the nodalpoints is avoided. Unlike other existing finite element formulationsthat lead to a highly nonlinear inertial forces for three-dimensionalelements, the proposed formulation leads to a constant mass matrix, andas a result, the centrifugal and Coriolis inertia forces are identicallyequal to zero. Furthermore, the method relaxes some of the assumptionsused in the classical and Mindlin plate models and automatically satisfiesthe objectivity requirements. By using a generalcontinuum mechanics approach, a relatively simple expression for theelastic forces is obtained. Generalization of the formulation to thecase of shell elements is discussed. As examples of the implementationof the proposed method, two different plate elements are presented; oneplate element does not guarantee the continuity of the displacementgradients between the nodal points, while the other plate elementguarantees this continuity. Numerical results are presented in order todemonstrate the use of the proposed method in the large rotation anddeformation analysis of plates and shells. The numerical results, whichare compared with the results obtained using existing incrementalprocedures, show that the solution obtained using the proposed methodsatisfies the principle of work and energy. These results are obtainedusing explicit numerical integration method. Potential applications ofthe proposed method include high-speed metal forming, vehiclecrashworthiness, rotor blades, and tires.

Analysis of Thin Plate Structures Using the Absolute Nodal Coordinate Formulation

Abstract 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.

Use of the finite element absolute nodal coordinate formulation in modeling slope discontinuity

A large rigid body rotation of a finite element can be described by rotating the axes of the element coordinate system or by keeping the axes unchanged and change the slopes or the position vector gradients. In the first method, the definition of the local element parameters (spatial coordinates) changes with respect to a body or a global coordinate system. The use of this method will always lead to a nonlinear mass matrix and non-zero centrifugal and Coriolis forces. The second method, in which the axes of the element coordinate system do not rotate with respect to the body or the global coordinate system, leads to a constant mass matrix and zero centrifugal and Coriolis forces when the absolute nodal coordinate formulation is used. This important property remains in effect even in the case of flexible bodies with slope discontinuities. The concept employed to accomplish this goal resembles the concept of the intermediate element coordinate system previously adopted in the finite element floating frame of reference formulation. It is shown in this paper that the absolute nodal coordinate formulation that leads to exact representation of the rigid body dynamics can be effecitively used in the analysis of complex structures with slope discontinuities. The analysis presented in this paper also demonstrates that objectivity is not an issue when the absolute nodal coordinate formulation is used due to the fact that this formulation automatically accounts for the proper coordinate transformations.

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

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 J 2 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.

General Method for Modeling Slope Discontinuities and T-Sections Using ANCF Gradient Deficient Finite Elements

Slope discontinuities and T-sections can be modeled in a straight forward manner using fully parameterized absolute nodal coordinate formulation (ANCF) finite elements that have a complete set of gradient vectors. Linear transformations that define the element connectivity can always be obtained and used to preserve ANCF desirable features that include constant mass matrix and zero Coriolis and centrifugal forces in the case of spinning structures. The objective of this paper is to develop a general method that allows for modeling slope discontinuities and T-sections using gradient deficient ANCF finite elements that do not have a complete set of coordinate lines and gradient vectors. Linear connectivity conditions that preserve all the ANCF desirable features including the constant mass matrix are developed at the nodes of slope discontinuities. At these nodes of discontinuity, one can always define a complete set of independent coordinate lines that lie on the structure. These coordinate lines can be used to define a complete set of independent gradient vectors at these nodes. Since the proposed method is based on linear coordinate transformations, the method can be implemented in a preprocessor computer program. The application of the proposed general method is demonstrated using ANCF gradient deficient beam element example.