Two fundamental challenges in investigation of nonlinear behavior of cantilever beam are the reliability of developed theory in facing with the reality and selecting the proper assumptions for solving the theory-provided equation. In this study, one of the most applicable theory and assumption for analyzing the nonlinear behavior of the cantilever beam is examined analytically and experimentally. The theory is concerned with the slender inextensible cantilever beam with large deformation nonlinearity, and the assumption is using the first-mode discretization in dealing with the partial differential equation provided by the theory. In the analytical study, firstly the equation of motion is derived based on the theory of large deformable inextensible beam. Then, the partial differential equation of motion is discretized using the Galerkin method via the assumption of the first mode. An exact solution to the obtained nonlinear ordinary differential equation is developed, because the available semi analytical and approximated methods, due to their limitations, are not always sufficiently reliable. Finally, an experiment set-up is developed to measure the nonlinear frequency of oscillations of an aluminum beam within a domain of initial displacement. The results show that the proposed analytical method has excellent convergence with experimental data.
The central theme of this work was to analyze high aspect ratio structure having structural nonlinearity in low subsonic flow and to model nonlinear stiffness by finite element-modal approach. Total stiffness of high aspect ratio wing can be decomposed to linear and nonlinear stiffnesses. Linear stiffness is modeled by its eigenvalues and eigenvectors, while nonlinear stiffness is calculated by the method of combined Finite Element-Modal approach. The nonlinear modal stiffness is calculated by defining nonlinear static load cases first. The nonlinear stiffness in the present work is modeled in two ways, i.e., based on bending modes only and based on bending and torsion modes both. Doublet lattice method (DLM) is used for dynamic analysis which accounts for the dependency of aerodynamic forces and moments on the frequency content of dynamic motion. Minimum state rational fraction approximation (RFA) of the aerodynamic influence coefficient (AIC) matrix is used to formulate full aeroelastic state-space time domain equation. Time domain dynamics analyses show that structure behavior becomes exponentially growing at speed above the flutter speed when linear stiffness is considered, however, Limit Cycle Oscillations (LCO) is observed when linear stiffness along with nonlinear stiffness, modeled by FE-Modal approach is considered. The amplitude of LCO increases with the increase in the speed. This method is based on cantilevered configuration. Nonlinear static tests are generated while wing root chord is fixed in all degrees of freedom and it needs modification if one requires considering full aircraft. It uses dedicated commercial finite element package in conjunction with commercial aeroelastic package making the method very attractive for quick nonlinear aeroelastic analysis. It is the extension of M.Y. Harmin and J.E. Cooper method in which they used the same equations of motion and modeled geometrical nonlinearity in bending modes only. In the current work, geometrical nonlinearities in bending and in torsion modes have been considered.
The paper presents the author’s non-linear FEM solution of an initially stressless deformed flat frame element, in which the nodes are situated along the axis of the bar initially straight. It has been assumed that each node may sustain arbitrary displacements and rotation. The solution takes into account the effect of shear, the geometrical non-linearity with large displacements (Green-Lagrange’s strain tensor) and moderate rotations (i.e. such ones which allow a linear-elastic behaviour of the material) and alternative small rotations when the second Piola-Kirchhoff stress tensor is applied. This solution is based on , concerning beams without any initial bow imperfections. The convergence of the obtained results at different numbers of nodes and Gauss points in the element was tested basing on the example of circular arcs with a central angle of 120°÷180°. The analysis concerned elements with two, three, five, seven, nine and eleven nodes, for the same number of points of numerical integration and also with one more or less. Moreover, the effect of distributing the load on the convergence of the results was analyzed.