Studies on the Stiffness Properties of the Absolute Nodal Coordinate Formulation for Three-Dimensional Beams

Author(s):  
Jussi T. Sopanen ◽  
Aki M. Mikkola

The objective of this study is to investigate the accuracy of elastic force models that can be used in the absolute nodal coordinate finite element formulation for the analysis of threedimensional beams. The elastic forces of the absolute nodal coordinate formulation can be derived using a continuum mechanics approach. This study investigates the accuracy and usability of such an approach for the three-dimensional absolute nodal coordinate beam element. This study also presents an improvement proposal for the use of a continuum mechanics approach in deriving the expression of the elastic forces of the beam element. The improvement proposal is verified using several numerical examples. Numerical examples show that the proposed elastic force model of the beam element agrees with analytical results as well as with solutions obtained using existing finite element formulation. The results also imply that the beam element does not suffer from the phenomenon called shear locking. In the beam element under investigation, global displacements and slopes are used as the nodal coordinates, which resulted in a large number of nodal degrees of freedom. This study provides a physical interpretation of the nodal coordinates used in the absolute nodal coordinate beam element. It is shown that the beam element based on the absolute nodal coordinate formulation relaxes the assumption of the rigid cross-section and is capable of representing a distortional deformation of the cross-section.

Author(s):  
Jussi Sopanen ◽  
Marko Matikainen ◽  
Aki Mikkola

In the absolute nodal coordinate formulation, the transverse shear deformation can be accounted for by using a fully parametrized element or, alternatively, by replacing longitudinal slope coordinates by a vector that describes the orientation of the cross-section. The use of a fully parametrized element allows the description of cross-section deformations in the case of beams and, correspondingly, fiber deformations in the case of plates and shells. It is noteworthy, however, that cross-section or fiber deformations are usually associated with high natural frequencies complicating the time integration of a fully parametrized element. A procedure to replace longitudinal slope coordinates by the vector that describes cross-section orientation was recently applied to a two-dimensional beam element based on the absolute nodal coordinate formulation. In this study, the procedure to account for shear deformation using the vector that describes cross-section orientation is extended to account for the nonlinear strain-displacement relationship in the definition of the elastic forces of the beam element. To accomplish this, the exact displacement field is used in the description of element kinematical and strain measures. This makes it straightforward to implement the non-linear strain-displacement relationship in the description of the elastic forces. Numerical results demonstrate that the enhanced beam element yields accurate results in eigenfrequency analysis. Results obtained in large deformation cases are in line with previously proposed elements based on the absolute nodal coordinate formulation.


Author(s):  
Haidong Yu ◽  
Chunzhang Zhao ◽  
Bin Zheng ◽  
Hao Wang

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.


Author(s):  
Marko K. Matikainen ◽  
Aki M. Mikkola

In this study, the improved description of elastic forces for the absolute nodal coordinate based plate element is introduced. The absolute nodal coordinate formulation, which utilizes global displacements and slope coordinates as nodal variables, can be used in large rotation and deformation dynamic analysis of beam and plate structures. The formulation avoids difficulties that arise when a rotation is interpolated in three-dimensional applications. In the absolute nodal coordinate formulation, a continuum mechanics approach has become the dominating procedure when elastic forces are defined. It has recently been perceived, however, that the continuum mechanics based absolute nodal coordinate elements suffer from serious shortcomings, including Poisson’s locking and poor convergence rate. These problems can be circumvented by modifying the displacement field of a finite element in the definition of elastic forces. This allows the use of the mixed type interpolation technique, leading to accurate and efficient finite element formulations. This approach has been previously applied to two- and three-dimensional absolute nodal coordinate based finite elements. In this study, the improved approach for elastic forces is extended to the absolute nodal coordinate plate element. The introduced plate element is compared in static examples to the continuum mechanics based absolute nodal coordinate plate element, as well as to commercial finite element software.


Author(s):  
Abdel-Nasser A. Mohamed ◽  
Jeff Liu

In this investigation, a three dimensional gradient deficient beam element (BEAM9) using the absolute nodal coordinate formulation (ANCF) is introduced. This element has nine coordinates per node, this includes the position vector and the two gradient vectors rx and ry. Like most of the ANCF elements, this element has constant mass matrix and zero centrifugal and Coriolis inertia forces. The plane strain elastic force model and the elastic line approach are two elastic force models presented in this paper in order to simulate the element internal resistance. Both models support resistance to the general bending and twist moments. The possibilities of employing nonlinear material models will be discussed in future work. Furthermore, the proposed element has the advantage of easy integration over general cross section area that is not easy to perform using the fully parameterized ANCF beam element (BEAM12). Comparing to the ANCF cable element (BEAM6), the proposed element can resist general bending and twist loads. Moreover, shear deformations in the xy plane due to shear force and in the yz plane due to twist moment are considered with the gradient deficient beam element proposed in this work. However, no shear deformations are considered with the ANCF cable element. Comparing to the fully parameterized ANCF beam element, the gradient deficient beam element (BEAM9) avoids some locking issues, shows better computational efficiency and offers better convergency characteristics. Numerical examples are presented in order to validate the proposed gradient deficient beam element and to compare with other ANCF beam elements.


Author(s):  
Zhenxing Shen ◽  
Xiaofeng Xing ◽  
Boyu Li

A novel modelling approach to beams with thin cross-sections is proposed in the absolute nodal coordinate formulation (ANCF), where the Lagrange interpolating and curve fitting techniques of numerical analysis are utilized for construction of the thin beam cross-section contour. Although the slope vector with respect to the coordinate line on cross-section contour is not considered in nodal coordinates, the cross-section distortion could be adequately captured through selecting an appropriate degree of polynomial. The main advantages of the present ANCF thin beam element are that the computational costs are more inexpensive than the ANCF shell element due to less generalized coordinates, there is very small amount of input data because slope vectors of the cross-section are eliminated, and the cross-sectional stress distribution may always be continuous on account of the fact that the cross-section is not discretized into a set of finite elements. Moreover, the formulations of elastic forces and Jacobian of thin laminated composite beam are also derived based on the continuum mechanics. Finally, several examples including both static and dynamic problems are performed to verify the new element and meanwhile demonstrate its general characteristics.


2005 ◽  
Vol 1 (2) ◽  
pp. 103-108 ◽  
Author(s):  
Aki M. Mikkola ◽  
Marko K. Matikainen

Dynamic analysis of large rotation and deformation can be carried out using the absolute nodal coordinate formulation. This formulation, which utilizes global displacements and slope coordinates as nodal variables, make it possible to avoid the difficulties that arise when a rotation is interpolated in three-dimensional applications. In the absolute nodal coordinate formulation, a continuum mechanics approach has become the dominating procedure when elastic forces are defined. It has recently been perceived, however, that the continuum mechanics based absolute nodal coordinate elements suffer from serious shortcomings, including Poisson’s locking and poor convergence rate. These problems can be circumvented by modifying the displacement field of a finite element in the definition of elastic forces. This allows the use of the mixed type interpolation technique, leading to accurate and efficient finite element formulations. This approach has been previously applied to two- and three-dimensional absolute nodal coordinate based finite elements. In this study, the improved approach for elastic forces is extended to the absolute nodal coordinate plate element. The introduced plate element is compared in static examples to the continuum mechanics based absolute nodal coordinate plate element, as well as to commercial finite element software. A simple dynamic analysis is performed using the introduced element in order to demonstrate the capability of the element to conserve energy.


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):  
R. Y. Yakoub ◽  
A. A. Shabana

Abstract By utilizing the fact that the absolute nodal coordinate formulation leads to a constant mass matrix, a Cholesky decomposition of the mass matrix can be used to obtain a constant velocity transformation matrix. This velocity transformation can be used to express the absolute nodal coordinates in terms of the generalized Cholesky coordinates. In this case, 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 motions. The implementation of a computer procedure based on the absolute nodal coordinate formulation and Cholesky coordinates is discussed in this paper. A flexible four-bar linkage is presented in this paper in order to demonstrate the use of Cholesky coordinates in the simulation of the small and large deformations in flexible multibody applications. The results obtained from the absolute nodal coordinate formulation are compared to those obtained from the floating frame of reference formulation.


Sign in / Sign up

Export Citation Format

Share Document