Cross-Section Deformation in the Absolute Nodal Coordinate Formulation

Author(s):  
Hiroyuki Sugiyama ◽  
Johannes Gerstmayr ◽  
Ahmed A. Shabana

In this investigation, the deformation modes defined in the finite element absolute nodal coordinate formulation using several strain definitions are discussed. In order to accurately define strain components that can have easy physical interpretation, a material coordinate system is introduced to define the material element rotation and deformation. The results obtained in this study clearly show cross-section deformation modes eliminated when the number of the finite element nodal coordinates is systematically and consistently reduced. Using the procedure discussed in this paper, one can obtain a reduced order dynamic model, eliminate position vector gradients that introduce high frequencies to the solution of some problems, achieve the continuity of the remaining gradients at the nodal points, and obtain a formulation that automatically satisfies the principle of work and energy.

Author(s):  
Bassam A. Hussein ◽  
Hiroyuki Sugiyama ◽  
Ahmed A. Shabana

The finite element absolute nodal coordinate formulation (ANCF) leads to beam and plate models that relax the assumption of the classical Euler-bernoulli and Timoshenko beam and Mindlin plate theories. In these more general models, the cross section is allowed to deform and it is no longer treated as a rigid surface. The coupling between the bending and cross section deformations leads to the new ANCF-coupled deformation modes that are examined in this study. While these coupled deformation can be source of numerical and convergence problems when thin and stiff beam models are considered, the inclusion of the effect of these modes in the dynamic model is necessary in the case of very flexible structures. In order to examine the effect of these coupled deformation modes in this investigation, three different large deformation dynamic beam models are discussed. Two of these models, which differ in the way the beam elastic forces are calculated in the absolute nodal coordinate formulation, allow for systematically eliminating the coupled deformation modes, while the third allows for including these modes. The first of these models is based on a general continuum mechanics approach that leads to a model that includes the ANCF-coupled deformation modes; while the second and third methods that can be used to eliminate the coupled deformation modes are based on the elastic line approach and the Hellinger-Reissner principle. It is shown in this study that the inclusion of the ANCF-coupled deformation modes introduces geometric stiffening effects that can not be captured using other finite element models.


2003 ◽  
Vol 126 (3) ◽  
pp. 478-487 ◽  
Author(s):  
Hiroyuki Sugiyama ◽  
Ahmed A. Shabana

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.


2006 ◽  
Vol 298 (4-5) ◽  
pp. 1129-1149 ◽  
Author(s):  
Hiroyuki Sugiyama ◽  
Johannes Gerstmayr ◽  
Ahmed A. Shabana

Author(s):  
Carmine M. Pappalardo ◽  
Zuqing Yu ◽  
Xiaoshun Zhang ◽  
Ahmed A. Shabana

In this paper, a rational absolute nodal coordinate formulation (RANCF) thin plate element is developed and its use in the analysis of curved geometry is demonstrated. RANCF finite elements are the rational counterpart of the nonrational absolute nodal coordinate formulation (ANCF) finite elements which employ rational polynomials as basis or blending functions. RANCF finite elements can be used in the accurate geometric modeling and analysis of flexible continuum bodies with complex geometrical shapes that cannot be correctly described using nonrational finite elements. In this investigation, the weights, which enter into the formulation of the RANCF finite element and form an additional set of geometric parameters, are assumed to be nonzero constants in order to accurately represent the initial geometry and at the same time preserve the desirable ANCF features, including a constant mass matrix and zero centrifugal and Coriolis generalized inertia forces. A procedure for defining the control points and weights of a Bezier surface defined in a parametric form is used in order to be able to efficiently create RANCF/ANCF FE meshes in a straightforward manner. This procedure leads to a set of linear algebraic equations whose solution defines the RANCF coordinates and weights without the need for an iterative procedure. In order to be able to correctly describe the ANCF and RANCF gradient deficient FE geometry, a square matrix of position vector gradients is formulated and used to calculate the FE elastic forces. As discussed in this paper, the proposed finite element allows for describing exactly circular and conic sections and can be effectively used in the geometry and analysis modeling of multibody system (MBS) components including tires. The proposed RANCF finite element is compared with other nonrational ANCF plate elements. Several numerical examples are presented in order to demonstrate the use of the proposed RANCF thin plate element. In particular, the FE models of a set of rational surfaces, which include conic sections and tires, are developed.


Author(s):  
Hiroyuki Sugiyama ◽  
Ahmed A. Shabana

The objective of this investigation is to develop a general nonlinear finite deformation 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. In addition, it is demonstrated that the principle of objectivity that is an issue when existing finite element formulations using ratetype constitutive equations are used can be fully 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.


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.


1999 ◽  
Vol 122 (4) ◽  
pp. 498-507 ◽  
Author(s):  
Marcello Campanelli ◽  
Marcello Berzeri ◽  
Ahmed A. Shabana

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]


Sign in / Sign up

Export Citation Format

Share Document