Improvement on Evaluating Axial Elastic Force in Bernoulli-Euler Beam Based on the Absolute Nodal Coordinate Formulation by Accurate Mean Axial Strain Measure

Author(s):  
Tsubasa Wago ◽  
Nobuyuki Kobayashi ◽  
Yoshiki Sugawara

This paper presents an improved formulation of axial elastic force in three-dimensional Bernoulli-Euler beam element based on the absolute nodal coordinate formulation. An accurate measure of mean axial strain for evaluating the axial elastic forces characterizes the presented formulation. The presented formulation evaluates the mean axial strain accurately by calculating the length of deformed beam element along its neutral axis. A comparison of the conventional formulations of the axial elastic force and the presented formulation is performed in some numerical examples which contain large bending deformation of flexible beam. As a result, it is verified that the presented formulation can express large deformation accurately with smaller number of elements than the conventional formulation which calculates the mean axial strain with straight-line distance between both element nodes. Moreover, it is also verified that the presented formulation can avoid excessive increase in computing time to simulate the dynamic behavior of flexible beam.

Author(s):  
Riki Iwai ◽  
Nobuyuki Kobayashi

This paper establishes a new type component mode synthesis method for a flexible beam element based on the absolute nodal coordinate formulation. The deformation of the beam element is defined as the sum of the global shape function and the analytical clamped-clamped beam modes. This formulation leads to a constant and symmetric mass matrix as the conventional absolute nodal coordinate formulation, and makes it possible to reduce the system coordinates of the beam structure which undergoes large rotations and large deformations. Numerical examples show that the excellent agreements are examined between the presented formulation and the conventional absolute nodal coordinate formulation. These results demonstrate that the presented formulation has high accuracy in the sense that the presented solutions are similar to the conventional ones with the less system coordinates and high efficiency in computation.


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):  
Yoshitaka Takahashi ◽  
Nobuyuki Shimizu

Abstract There are three basic finite element formulations which are used in the dynamics analysis of flexible beams undergoing large rotation and deformation. These are the floating frame of reference approach, the finite segment method and the large rotation vector approach. Recently, the absolute nodal coordinate formulation was proposed by A.A.Shabana et al. In this formulation, elastic forces are lead by approximating the slope of the beam at an arbitrary point on the neutral axis of the beam in terms of the slope of the simple support axis of the beam. In this paper, we propose the mean axis of the planar Bernoulli-Euler beam for the absolute nodal coordinate approach. The origin and the orientation of this axis are selected so as to minimize the total deformation of the concerned beam. And the selected axis can be simply described by the nodal coordinates of the beam element. Using the mean axis instead of the simple support axis, the elastic forces of the beam element may be more precisely calculated. Finally, we show numerical examples to demonstrate effectiveness of this approach.


Author(s):  
Hiroyuki Sugiyama ◽  
Hirohisa Koyama ◽  
Hiroki Yamashita

In this investigation, a gradient deficient beam element of the absolute nodal coordinate formulation is generalized to a curved beam for the analysis of multibody systems and the performance of the proposed element is discussed by comparing with the fully parameterized curved beam element and the classical large displacement beam element with incremental solution procedures. Strain components are defined with respect to the initially curved configuration and described by the arc-length coordinate. The Green strain is used for the longitudinal stretch, while the material measure of curvature is used for bending. It is shown that strains of the curved beam can be expressed with respect to those defined in the element coordinate system using the gradient transformation and the effect of strains at the initially curved configuration is eliminated using one-dimensional Almansi strain. This property can be effectively used with non-incremental solution procedure employed for the absolute nodal coordinate formulation. Several numerical examples are presented in order to demonstrate the performance of the gradient deficient curved beam element developed in this investigation. It is shown that the use of the proposed element leads to better element convergence as compared to that of the fully parameterized element and the classical large displacement beam element with incremental solution procedures.


Author(s):  
Aki Mikkola ◽  
Oleg Dmitrochenko ◽  
Marko Matikainen

In this study, a procedure to account for transverse shear deformation in the absolute nodal coordinate formulation is presented. In the absolute nodal coordinate formulation, shear deformation is usually defined by employing the slope vectors in the element transverse direction. This leads to the description of deformation modes that are, in practical problems, associated with high frequencies. These high frequencies, in turn, complicate the time integration procedure burdening numerical performance. In this study, the description of transverse shear deformation is accounted for in a two-dimensional beam element based on the absolute nodal coordinate formulation without the use of transverse slope vectors. In the introduced shear deformable beam element, slope vectors are replaced by vectors that describe the rotation of the beam cross-section. This procedure represents a simple enhancement that does not decrease the accuracy or numerical performance of elements based on the absolute nodal coordinate formulation. Numerical results are presented in order to demonstrate the accuracy of the introduced element in static and dynamic cases. The numerical results obtained using the introduced element agree with the results obtained using previously proposed shear deformable beam elements.


Author(s):  
Oleg Dmitrochenko ◽  
Aki Mikkola

This study is an extension of a newly introduced approach to account transverse shear deformation in absolute nodal coordinate formulation. In the formulation, shear deformation is usually defined by employing slope vectors in the element transverse direction. This leads to the description of deformation modes that, in practical problems, may be associated with high frequencies. These high frequencies, in turn, could complicate the time integration procedure, burdening numerical performance of shear deformable elements. In a recent study of this paper’s authors, the description of transverse shear deformation is accounted for in a two-dimensional beam element, based on the absolute nodal coordinate formulation without the use of transverse slope vectors. In the introduced shear deformable beam element, slope vectors are replaced by vectors that describe the rotation of the beam cross-section. This procedure represents a simple enhancement that does not decrease the accuracy or numerical performance of elements based on the absolute nodal coordinate formulation. In this study, the approach to account for shear deformation without using transverse slopes is implemented for a thin rectangular plate element. In fact, two new plate elements are introduced: one within conventional finite element and another using the absolute nodal coordinates. Numerical results are presented in order to demonstrate the accuracy of the introduced plate element. The numerical results obtained using the introduced element agree with the results obtained using previously proposed shear deformable plate elements.


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.


Sign in / Sign up

Export Citation Format

Share Document