Finite element analysis of the geometric stiffening effect. Part 1: A correction in the floating frame of reference formulation

Abstract

The fact that incorrect unstable solutions are obtained for linearly elastic models motivates the analytical study presented in this paper. The increase in the number of finite elements only leads to an increase in the critical speed. Crucial in the analysis presented in this paper is the fact that the mass matrix and the form of the elastic forces obtained using the absolute nodal coordinate formulation remain the same under orthogonal coordinate transformation. The absolute nodal coordinate formulation, in contrast to conventional finite element formulations, does account for the effect of the coupling between bending and extension. Based on the analytical results obtained using the absolute nodal coordinate formulation, a new correction is proposed for the finite element floating frame of reference formulation in order to introduce coupling between the axial and bending displacements. In this two-part paper, two- and three-dimensional finite element models are used to study the problem of rotating beams. The models are developed using the absolute nodal coordinate formulation that allows for accurate representation of the axial strain, thereby avoiding the ill-conditioning problem that arises when classical displacement-based finite element formulations are used. In the first part of the paper, the case of linear elasticity is considered and assumptions used in the finite element floating frame of reference formulation are investigated. In the second part of the paper, non-linear elasticity is considered. A rotating helicopter blade is simulated, and the complexity of the motion suggests the inclusion of rotary inertia, shear deformation, and non-linear elastic forces in order to obtain an accurate solution that does not suffer from the instability problem regardless of the number of finite elements used.