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.

Three Dimensional Absolute Nodal Coordinate Formulation for Beam Elements: Theory

The purpose of this paper is to present formulations for beam elements based on the absolute nodal co-ordinate formulation that can be effectively and efficiently used in the case of thin structural applications. The numerically stiff behaviour resulting from shear terms in existing absolute nodal co-ordinate formulation beam elements that employ the continuum mechanics approach to formulate the elastic forces and the resulting locking phenomenon make these elements less attractive for slender stiff structures. In this investigation, additional shape functions are introduced for an existing spatial absolute nodal co-ordinate formulation beam element in order to obtain higher accuracy when the continuum mechanics approach is used to formulate the elastic forces. For thin structures where bending stiffness can be important in some applications, a lower order cable element is introduced and the performance of this cable element is evaluated by comparing it with existing formulations using several examples. Cables that experience low tension or catenary systems where bending stiffness has an effect on the wave propagation are examples in which the low order cable element can be used. The cable element, which does not have torsional stiffness, can be effectively used in many problems such as in the formulation of the sliding joints in applications such as the spatial pantograph/catenary systems. The numerical study presented in this paper shows that the use of existing implicit time integration methods enables the simulation of multibody systems with a moderate number of thin and stiff finite elements in reasonable CPU time.

Analysis of Thin Beams and Cables Using the Absolute Nodal Co-ordinate Formulation

This part of these two companion papers demonstrates the computer implementation of the absolute nodal coordinate formulation for three-dimensional beam elements. Two beam elements that relax the assumptions of Euler-Bernoulli and Timoshenko beam theories are developed. These two elements take into account the effect of rotary inertia, shear deformation and torsion, and yet they lead to a constant mass matrix. As a consequence, the Coriolis and centrifugal forces are identically equal to zero. Both beam elements use the same interpolating polynomials and have the same number of nodal coordinates. However, one of the elements has two nodes, while the other has four nodes. The results obtained using the two elements are compared with the results obtained using existing incremental methods. Unlike existing large rotation vector formulations, the results of this paper show that no special numerical integration methods need to be used in order to satisfy the principle of work and energy when the absolute nodal coordinate formulation is used. These results show that this formulation can be used in manufacturing applications such as high speed forming and extrusion problems in which the element cross section dimensions significantly change.

Three Dimensional Absolute Nodal Coordinate Formulation for Beam Elements: Implementation and Applications

1. Introduction 2. Kinematics 3. Forces and stresses 4. Constitutive equations 5. Plasticity formulations 6. Finite element formulation: large deformation, large rotation problem 7. Finite element formulation: small deformation, large rotation problem.

/pdf/computational-continuum-mechanics-4a3n1xh7xx.pdf

Computational Continuum Mechanics

https://repositorium.sdum.uminho.pt/bitstream/1822/23521/1/6.7.15%202009.pdf

Dynamics of spatial flexible multibody systems with clearance and lubricated spherical joints

In a previous publication, procedures that can be used with the absolute nodal coordinate formulation to solve the dynamic problems of flexible multibody systems were proposed. One of these procedures is based on the Cholesky decomposition. By utilizing the fact that the absolute nodal coordinate formulation leads to a constant mass matrix, a Cholesky decomposition is used to obtain a constant velocity transformation matrix. This velocity transformation is used to express the absolute nodal coordinates in terms of the generalized Cholesky coordinates. 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 motion. The implementation of a computer procedure based on the absolute nodal coordinate formulation and Cholesky coordinates is discussed in this paper. Numerical examples are presented in order to demonstrate the use of Cholesky coordinates in the simulation of the large deformations in flexible multibody applications.

Use of Cholesky coordinates and the absolute nodal coordinate formulation in the computer simulation of flexible multibody systems

: The description of a beam element by only the displacement of its center line 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 multi body 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 effects, 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.

http://www.dtic.mil/cgi-bin/GetTRDoc?AD=ADA374867&Location=U2&doc=GetTRDoc.pdf

R. Y. Yakoub

Papers

Three Dimensional Absolute Nodal Coordinate Formulation for Beam Elements: Theory

Three Dimensional Absolute Nodal Coordinate Formulation for Beam Elements: Implementation and Applications

Use of Cholesky coordinates and the absolute nodal coordinate formulation in the computer simulation of flexible multibody systems

An Isoparametric Three Dimensional Beam Element Using the Absolute Nodal Coordinate Formulation

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