A Lagrange–Eulerian formulation of an axially moving beam based on the absolute nodal coordinate formulation

In the scope of this paper, a finite-element formulation for an axially moving beam is presented. The beam element is based on the absolute nodal coordinate formulation, where position and slope vectors are used as degrees of freedom instead of rotational parameters. The equations of motion for an axially moving beam are derived from generalized Lagrange equations in a Lagrange–Eulerian sense. This procedure yields equations which can be implemented as a straightforward augmentation to the standard equations of motion for a Bernoulli–Euler beam. Moreover, a contact model for frictional contact between an axially moving strip and rotating rolls is presented. To show the efficiency of the method, simulations of a belt drive are presented.     

Experimental validation of rigid-flexible coupling dynamic formulation for hub–beam system

The aim of this paper is to compare the accuracy of the absolute nodal coordinate formulation and the floating frame of reference formulation for the rigid-flexible coupling dynamics of a three-dimensional Euler–Bernoulli beam by numerical and experimental validation. In the absolute nodal coordinate formulation, based on geometrically exact beam theory and considering the torsion effect, the material curvature of the beam is derived, and then variational equations of motion of a three-dimensional beam are obtained, which consist of three position coordinates, two slope coordinates, and one rotational coordinate. In the floating frame of reference formulation, the displacement of an arbitrary point on the beam is described by the rigid-body motion and a small superimposed deformation displacement. Based on linear elastic theory, the quadratic terms of the axial strain are neglected, and the curvatures are simplified to the first order. Considering both the linear damping and the quadratic air resistance damping, the equations of motion of the multibody system composed of air-bearing test bed and a cantilevered three-dimensional beam are derived based on the principle of virtual work. In order to verify the results of the computer simulation, two experiments are carried out: an experiment of hub–beam system with large deformation and a dynamic stiffening experiment. The comparison of the simulation and experiment results shows that in case of large deformation, the frequency result obtained by the floating frame of reference formulation is lower than that obtained by the experiment. On the contrary, the result obtained by the absolute nodal coordinate formulation agrees well with that obtained by the experiment. It is also shown that the floating frame of reference formulation based on linear elastic theory cannot reveal the dynamic stiffening effect. Finally, the applicability of the floating frame of reference formulation is clarified.     

A geometrically exact beam element based on the absolute nodal coordinate formulation

In this study, Reissner’s classical nonlinear rod formulation, as implemented by Simo and Vu-Quoc by means of the large rotation vector approach, is implemented into the framework of the absolute nodal coordinate formulation. The implementation is accomplished in the planar case accounting for coupled axial, bending, and shear deformation. By employing the virtual work of elastic forces similarly to Simo and Vu-Quoc in the absolute nodal coordinate formulation, the numerical results of the formulation are identical to those of the large rotation vector formulation. It is noteworthy, however, that the material definition in the absolute nodal coordinate formulation can differ from the material definition used in Reissner’s beam formulation. Based on an analytical eigenvalue analysis, it turns out that the high frequencies of cross section deformation modes in the absolute nodal coordinate formulation are only slightly higher than frequencies of common shear modes, which are present in the classical large rotation vector formulation of Simo and Vu-Quoc, as well. Thus, previous claims that the absolute nodal coordinate formulation is inefficient or would lead to ill-conditioned finite element matrices, as compared to classical approaches, could be refuted. In the introduced beam element, locking is prevented by means of reduced integration of certain parts of the elastic forces. Several classical large deformation static and dynamic examples as well as an eigenvalue analysis document the equivalence of classical nonlinear rod theories and the absolute nodal coordinate formulation for the case of appropriate material definitions. The results also agree highly with those computed in commercial finite element codes.     

Analysis of Large Flexible Body Deformation in Multibody Systems Using Absolute Coordinates

To consider large deformation problems in multibody system simulations afinite element approach, called absolute nodal coordinate.formulation,has been proposed. In this formulation absolute nodal coordinates andtheir material derivatives are applied to represent both deformation andrigid body motion. The choice of nodal variables allows a fullynonlinear representation of rigid body motion and can provide the exactrigid body inertia in the case of large rotations. The methodology isespecially suited for but not limited to modeling of beams, cables andshells in multibody dynamics.This paper summarizes the absolute nodal coordinate formulation for a 3D Euler–Bernoulli beam model, in particular the definition of nodal variables, corresponding generalized elastic and inertia forces and equations of motion. The element stiffness matrix is a nonlinear function of the nodal variables even in the case of linearized strain/displacement relations. Nonlinear strain/displacement relations can be calculated from the global displacements using quadrature formulae.Computational examples are given which demonstrate the capabilities of the applied methodology. Consequences of the choice of shape.functions on the representation of internal forces are discussed. Linearized strain/displacement modeling is compared to the nonlinear approach and significant advantages of the latter, when using the absolute nodal coordinate formulation, are outlined.     

Effect of the centrifugal forces on the finite element eigenvalue solution of a rotating blade: a comparative study

In this study, the effect of the centrifugal forces on the eigenvalue solution obtained using two different nonlinear finite element formulations is examined. Both formulations can correctly describe arbitrary rigid body displacements and can be used in the large deformation analysis. The first formulation is based on the geometrically exact beam theory, which assumes that the cross section does not deform in its own plane and remains plane after deformation. The second formulation, the absolute nodal coordinate formulation (ANCF), relaxes this assumption and introduces modes that couple the deformation of the cross section and the axial and bending deformations. In the absolute nodal coordinate formulation, four different models are developed; a beam model based on a general continuum mechanics approach, a beam model based on an elastic line approach, a beam model based on an elastic line approach combined with the Hellinger–Reissner principle, and a plate model based on a general continuum mechanics approach. The use of the general continuum mechanics approach leads to a model that includes the ANCF coupled deformation modes. Because of these modes, the continuum mechanics model differs from the models based on the elastic line approach. In both the geometrically exact beam and the absolute nodal coordinate formulations, the centrifugal forces are formulated in terms of the element nodal coordinates. The effect of the centrifugal forces on the flap and lag modes of the rotating beam is examined, and the results obtained using the two formulations are compared for different values of the beam angular velocity. The numerical comparative study presented in this investigation shows that when the effect of some ANCF coupled deformation modes is neglected, the eigenvalue solutions obtained using the geometrically exact beam and the absolute nodal coordinate formulations are in a good agreement. The results also show that as the effect of the centrifugal forces, which tend to increase the beam stiffness, increases, the effect of the ANCF coupled deformation modes on the computed eigenvalues becomes less significant. It is shown in this paper that when the effect of the Poisson ration is neglected, the eigenvalue solution obtained using the absolute nodal coordinate formulation based on a general continuum mechanics approach is in a good agreement with the solution obtained using the geometrically exact beam model.     

A corotational procedure that handles large rotations of spatial beam structures

A practical motion process of the three dimensional beam element is presented to remove the restriction of small rotations between two successive increments for large displacement and large rotation analysis of space frames using incremental-iterative methods. In order to improve convergence properties of the equilibrium iteration, an n-cycle iteration scheme is introduced.

The nonlinear formulation is based on the corotational formulation. The transformation of the element coordinate system is assumed to be accomplished by a translation and two successive rigid body rotations: a transverse rotation followed by an axial rotation. The element formulation is derived based on the small deflection beam theory with the inclusion of the effect of axial force in the element coordinate system. The membrane strain along the deformed beam axis obtained from the elongation of the arc length of the beam element is assumed to be constant. The element internal nodal forces are calculated using the total deformational nodal rotations. Two methods, referred to as direct method and incremental method, are proposed in this paper to calculate the total deformational rotations.

An incremental-iterative method based on the Newton-Raphson method combined with arc length control is adopted. Numerical studies are presented to demonstrate the accuracy and efficiency of the present method.     

Definition of the Elastic Forces in the Finite-Element Absolute Nodal Coordinate Formulation and the Floating Frame of Reference Formulation

The equivalence of the finite-element formulations used inflexible multibody dynamics is the focus of this investigation. Thisequivalence will be used to address several fundamental issues related tothe deformations, flexible body coordinate systems, and the geometriccentrifugal stiffening effect. Two conceptually different finite-elementformulations that lead to exact modeling of the rigid body dynamics will beused. The first one is the absolute nodal coordinateformulation in which beams and plates can be treated as isoparametricelements. This formulation leads to a constant and symmetric mass matrix andhighly nonlinear elastic forces. In this study, it is demonstrated thatdifferent element coordinate systems which are used for the convenience ofdescribing the element deformations lead to similar results as the elementsize is reduced. In particular, two element frames are used;the pinned and the tangent frames. The pinned frame has one ofits axes passing through two nodes of the element, while the tangent frame isrigidly attached to one of the ends of the element. Numerical resultsobtained using these two different frames are found tobe in good agreement as the element size decreases. The relationshipbetween the coordinates used in the absolute nodal coordinate formulationand the floating frame of reference formulation is presented. Thisrelationship can be used to obtain the highly nonlinear expression of thestrain energy used in the absolute nodal coordinate formulation from thesimple energy expression used in the floating frame of referenceformulation. It is also shown that the source of the nonlinearityis due to the finite rotation of the element. The result of the analysispresented clearly demonstrates that the instability observedin high-speed rotor analytical models due to the neglect of the geometriccentrifugal stiffening is not a problem inherent to a particular finite-element formulation. Such a problem can only be avoided by considering the known linear effect of the geometric centrifugal stiffening or by using a nonlinear elastic model as recently demonstrated. 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. Another method forformulating the elastic forces in the absolute nodal coordinate formulationbased on a continuum mechanics approach is also presented.     

A formal procedure and invariants of a transition from conventional finite elements to the absolute nodal coordinate formulation

In this paper, finite elements based on the absolute nodal coordinate formulation (ANCF) are studied. The formulation has been developed by various authors for the dynamical simulation of large-displacement and large-rotation problems in flexible multibody dynamics. This study introduces a procedure to track the general geometrical properties of ANCF elements back to their prototypes in the conventional finite-element method (FEM), which deals with small-displacement problems. In this study, it is shown that each known ANCF element can be derived from a conventional FEM using a universal transform. Moreover, some important static and dynamic properties of the elements in small-displacement problems are automatically preserved. In the past, the authors of each newly proposed ANCF element have made unnecessary efforts to show the consistency of the above mentioned properties.     

Poisson modes and general nonlinear constitutive models in the large displacement analysis of beams

Most existing formulations for structural elements such as beams, plates and shells do not allow for the use of general nonlinear constitutive models in a straightforward manner. Furthermore, such structural element models, due to the nature of the generalized coordinates used, do not capture some Poisson modes such as the ones that couple the deformation of the cross section of the structural element and stretch and bending. In this paper, beam models that employ general nonlinear constitutive equations are presented using finite elements based on the nonlinear absolute nodal coordinate formulation. This formulation relaxes the assumptions of the Euler–Bernoulli and Timoshenko beam theories, and allows for the use of general nonlinear constitutive models. The finite elements based on the absolute nodal coordinate formulation also allow for the rotation as well as the deformation of the cross section, thereby capturing Poisson modes which can not be captured using other beam models. In this investigation, three different nonlinear constitutive models based on the hyper-elasticity theory are considered. These three models are based on the Neo–Hookean constitutive law for compressible materials, the Neo–Hookean constitutive law for incompressible materials, and the Mooney–Rivlin constitutive law in which the material is assumed to be incompressible. These models, which allow capturing Poisson modes, are suitable for many materials and applications, including rubber-like materials and biological tissues which are governed by nonlinear elastic behavior. Numerical examples that demonstrate the implementation of these nonlinear constitutive models in the absolute nodal coordinate formulation are presented. The results obtained using the nonlinear and linear constitutive models are compared in this study. These results show that the use of nonlinear constitutive models can significantly enhance the performance and improve the computational efficiency of the finite element models based on the absolute nodal coordinate formulation. The results also show that when linear constitutive models are used in the large deformation analysis, singular configurations are encountered and basic formulas such as Nanson’s formula are no longer valid. These singular deformation configurations are not encountered when the nonlinear constitutive models are used.     

Rigid-flexible coupling dynamics of three-dimensional hub-beams system

In the previous research of the coupling dynamics of a hub-beam system, coupling between the rotational motion of hub and the torsion deformation of beam is not taken into account since the system undergoes planar motion. Due to the small longitudinal deformation, coupling between the rotational motion of hub and the longitudinal deformation of beam is also neglected. In this paper, rigid-flexible coupling dynamics is extended to a hub-beams system with three-dimensional large overall motion. Not only coupling between the large overall motion and the bending deformation, but also coupling between the large overall motion and the torsional deformation are taken into account. In case of temperature increase, the longitudinal deformation caused by the thermal expansion is significant, such that coupling between the large overall motion and the longitudinal deformation is also investigated. Combining the characteristics of the hybrid coordinate formulation and the absolute nodal coordinate formulation, the system generalized coordinates include the relative nodal displacement and the slope of each beam element with respect to the body-fixed frame of the hub, and the variables related to the spatial large overall motion of the hub and beams. Based on precise strain-displacement relation, the geometric stiffening effect is taken into account, and the rigid-flexible coupling dynamic equations are derived using velocity variational principle. Finite element method is employed for discretization. Simulation of a hub-beams system is used to show the coupling effect between the large overall motion and the torsional deformation as well as the longitudinal deformation. Furthermore, conservation of energy in case of free motion is shown to verify the formulation.     

Nonlinear finite element analysis of elastic frames

A very simple and effective formulation and numerical procedure to remove the restriction of small rotations between two successive increments for the geometrically nonlinear finite element analysis of in-plane frames is presented. A co-rotational formulation combined with small deflection beam theory with the inclusion of the effect of axial force is adopted. A body attached coordinate is used to distinguish between rigid body and deformational rotations. The deformational nodal rotational angles are assumed to be small, and the membrane strain along the deformed beam axis obtained from the elongation of the arc length of the deformed beam element is assumed to be constant. The element internal nodal forces are calculated using the total deformational nodal rotations in the body attached coordinate. The element stiffness matrix is obtained by superimposing the bending and the geometric stiffness matrices of the elementary beam element and the stiffness matrix of the linear bar element. An incremental iterative method based on the Newton-Raphson method combined with a constant arc length control method is employed for the solution of the nonlinear equilibrium equations. In order to improve convergence properties of the equilibrium iteration, a two-cycle iteration scheme is introduced. Numerical examples are presented to demonstrate the accuracy and efficiency of the proposed method.     

Definition of the Slopes and the Finite Element Absolute Nodal Coordinate Formulation

In the finite element literature, beams and platesare not treated as isoparametric elements becausethe displacement gradients at the nodes are used todescribe infinitesimal rotations. Such alinearization of the slopes leads to incorrect rigidbody equations of motion when the beams and platesrotate as rigid bodies. A simple and efficient absolute nodal coordinate formulation, in which noinfinitesimal or finite rotations are used as nodalcoordinates, can be systematically developed andefficiently used in many large deformationapplications as well as in the analysis of curvedstructures. This formulation leads to a constantmass matrix, and as a consequence, an efficientprocedure can be used for solving for the nodalaccelerations. In the absolute nodal coordinateformulation, the displacement gradients at the nodesare determined in the undeformed referenceconfiguration using simple rigid body kinematics.The problems that arise from the linearization ofthe slopes and the use of the finite rotations asnodal coordinates are discussed in this paper.     

The development of a sliding joint for very flexible multibody dynamics using absolute nodal coordinate formulation

In this paper, a formulation for a spatial sliding joint is derived using absolute nodal coordinates and non-generalized coordinate and it allows a general multibody move along a very flexible cable. The large deformable motion of a spatial cable is presented using absolute nodal coordinate formulation, which is based on the finite element procedures and the general continuum mechanics theory to represent the elastic forces. And the nongeneralized coordinate, which is related to neither the inertia forces nor the external forces, is used to describe an arbitrary position along the centerline of a very flexible cable. Hereby, the non-generalized coordinate represents the arc-length parameter. The constraint equations for the sliding joint are expressed in terms of generalized coordinate and nongeneralized coordinate. In the constraint equations for the sliding joint, one constraint equation can be systematically eliminated. There are two independent Lagrange multipliers in the final system equations of motion associated with the sliding joint. The development of this sliding joint is important to analyze many mechanical systems such as pulley systems and pantograph-catenary systems for high speed-trains.     

A new locking-free formulation for planar,shear deformable,linear and quadratic beam finite elements based on the absolute nodal coordinate formulation

Many widely used beam finite element formulations are based either on Reissner’s classical nonlinear rod theory or the absolute nodal coordinate formulation (ANCF). Advantages of the second method have been pointed out by several authors; among the benefits are the constant mass matrix of ANCF elements, the isoparametric approach and the existence of a consistent displacement field along the whole cross section. Consistency of the displacement field allows simpler, alternative formulations for contact problems or inelastic materials. Despite conceptional differences of the two formulations, the two models are unified in the present paper.     

Three- and four-noded planar elements using absolute nodal coordinate formulation

This paper investigates two new types of planar finite elements containing three and four nodes. These elements are the reduced forms of the spatial plate elements employing the absolute nodal coordinate approach. Elements of the first type use translations of nodes and global slopes as nodal coordinates and have 18 and 24 degrees of freedom. The slopes facilitate the prevention of the shear locking effect in bending problems. Furthermore, the slopes accurately describe the deformed shape of the elements. Triangular and quadrilateral elements of the second type use translational degrees of freedom only and, therefore, can be utilized successfully in problems without bending. These simple elements with 6 and 8 degrees of freedom are identical to the elements used in conventional formulation of the finite element method from the kinematical point of view. Similarly to the famous problem called “flying spaghetti” which is used often as a benchmark for beam elements, a kind of “flying lasagna” is simulated for the planar elements. Numerical results of simulations are presented.