The flutter of a two-dimensional lifting surfaces with cubic structural and aerodynamic nonlinearities in a supersonic flow field is considered. By using the maple program, the normal form of a general system of motion of a rigid airfoil with flight Mach number, nonlinear stiffness coefficient and isentropic gas coefficient is computed. The bifurcation theory is used to obtain the universal unfolding of the normal form, and all the dynamic behavior of the universal unfolding is discussed. In addition, the flutter boundary and its character are derived. This effective analysis method is applied to two numerical examples, and the influence of the structural and aerodynamic parameters on the character of flutter instability is investigated. Finally, the numerical simulations obtained by using fourth-order Runge-Kutta method will illustrate the quality of this analysis method.
In recent years, many problems concerning the 2-degree-of-freedom (2-DOF) aeroelatstic system with cubic nonlinearity have been successfully studied using the center manifold theory and the normal form method [1-3]. There exist many potential sources of nonlinearities, which can have significant effect on an aeroelastic response. The nonlinearities of the aeroelastic system can be structural [4-5] or aerodynamic .
Breitbach  presented a detailed review of some possible structural nonlinearities and their effect on aeroelastically induced vibrations. Raghothama and Narayanan  investigated periodic oscillations and bifurcations of a two-dimensional airfoil in plunge and pitch motions with cubic pitch stiffness in incompressible flow. Shahrzad and Mahzoon  studied the limit cycle oscillation (LCO) flutter of a two-dimensional wing with nonlinear pitch stiffness. Liu, Wong and Lee  adopted a point transformation technique to study the nonlinear behavior of a two-dimensional aeroelastic system with freeplay models. Zhao and Zhang  investigated the LCO flutter and chaotic motion of a 2-DOF airfoil system with a combination of freeplay and cubic nonlinearities. The Krylov-Bogoliubov-Mitropolsky (KBM) method was used by Yang  to study the LCOs of an airfoil with the freeplay nonlinearity. Kim and Lee  analyzed the dynamic behavior of a flexible airfoil with the freeplay nonlinearity. Price et al.  discovered the LCOs at flow speeds below the linear flutter boundary in a system with a hard pitch spring. Abbas et al.  investigated the flutter and post-flutter of a two-dimensional double-wedge lifting surface. The behavior of bifurcation point and the influence of the system parameters on the stability and amplitudes of the LCO flutter were analyzed in references [15-16]. The relation between frequency and airfoil parameters was obtained by means of the harmonic balance method in reference .
In this paper, an effective analysis method to the problem of flutter of a two-dimensional airfoil with cubic structural and aerodynamic nonlinearities is proposed. This analysis method enables one to make a parametric study. Firstly, the normal form and universal unfolding are deduced by means of the maple program and the bifurcation theory. Secondly, super- and sub-critical Hopf bifurcations are analyzed. Finally, numerical examples and simulations are given to show the efficiency of this analysis method. Moreover, the influence of the structural and aerodynamic parameters, representations , , and , on the airfoil safety is studied.
Fig. 1Airfoil section
2. Nonlinear model of a two-dimensional lifting surface
Consider a two-dimensional lifting surface featuring plunging and twisting degrees of freedom, elastically constrained by a linear translational spring and nonlinear torsional spring, as shown in Fig. 1, exposed to a supersonic flow field. The equations of motion of the lifting surface are :
where is the plunging displacement (positive downward), the twist angle about the pitch axis (positive nose up), the primes denote differentiation with respect to time , the structural mass per unit span, the static unbalance about the elastic axis, the mass moment of inertia about the elastic axis of the airfoil, and are linear plunging and pitching viscous damping coefficients, while is the plunging stiffness coefficient. And represents the overall restoring moment that is connected to the pitch angle by:
where and are the linear and nonlinear pitching stiffness coefficients, respectively. The tracer that identifies this type of nonlinearities can take the value 1 or 0 depending on whether the nonlinearity is accounted for or discarded, respectively. In addition, the expressions of the aerodynamic lift and moment are:
herein, is the half-chord length of the airfoil, is streamwise position of the pitch axis measured from the leading edge, , and are the air speed, the air density, and the flight Mach number of the disturbed flow, respectively. and are correction aerodynamic factor and the isentropic gas coefficient, respectively.
Upon denoting , , and the dimensionless time , Eq. (1) can be written as:
where the coefficients can be found in Appendix, and the dots denote .
3. The stability and bifurcation analysis of system (4)
3.1. Stability of the initial equilibrium point
Obviously, is an equilibrium point of Eq. (4). The Jacobian of the linearization system of system Eq. (4) evaluated at is:
The characteristic equation of the matrix (5) is:
According to the Routh-Hruwitz criterion, the initial equilibrium solution is stable if the following conditions are satisfied:
When conditions (8) are not satisfied, the initial equilibrium solution is unstable, and bifurcations may occur. First we give the following results on the eigenvalues of matrix (5).
Proposition 1. Zero is not an eigenvalue of matrix (5).
Proposition 2. There exist no two pairs of purely imaginary eigenvalues for matrix (5).
In fact, and due to all the parameters in the expressions of , are positive and . However suppose zero is an eigenvalue of matrix (5), then it follows, on the other hand, if there exist two pairs of purely imaginary eigenvalues for matrix (5), then it follows . So Propositions 1 and 2 are obvious.
For the aeroelastic stability problems, the condition should be satisfied. According to the above analysis, it is meaningful to consider the only one critical case: two eigenvalues with negative real parts and a pair of purely imaginary eigenvalues (this case corresponds to the limit cycle flutter). In this case, the following conditions are fulfilled:
3.2. Normal form, universal unfolding and Hopf bifurcation
Choosing the following dimensionless values of parameters as : , , , , , and , Eq. (4) can be rewritten as:
Substituting the above dimensionless values of parameters and the Eq. (7) into the first formula of (9), the speed on the flutter instability boundary is:
The relation of the flutter speed and flight Mach number is depicted in Fig. 2. The figure reveals that the increase of flight Mach number can result in the increase of the flutter safety.
Fig. 2The flutter speed VF vs. flight Mach number M
In the following, we assume is a bifurcation point, i.e. the conditions (9) hold. Substituting into Eq. (10), we get:
For Eq. (12), there is a coordinate transformation so that:
where , and are the eigenvalues of the matrix , and .
Using the Maple program in reference , the normal form of Eq. (13) is calculated as:
The coordinate transformation is omitted here (for details, see the Maple program in reference ).
If , the universal unfolding of Eq. (14) can be written as :
Because Eq. (16) is the universal unfolding of Eq. (14), Eq. (16) bears all the dynamic behavior of any small perturbed system of Eq. (14) in the vicinity of . Furthermore, using the transformations and , we can obtain all the dynamic behavior of any small perturbed system of Eq. (12) in the vicinity of .
Let and , then the polar coordinate form of Eq. (16) is:
Now we derive the expression of the unfolding parameter in Eq. (17).
Consider one of the perturbed systems of Eq. (12) as follows. Without loss of generality, using the parameter transformations and the state variable transformation , one may transform Eq. (12) into a new system as follows:
Let . According to reference , in Eq. (17) the unfolding parameter . We can obtain the conclusions for Eqs. (16)-(18) as follows:
Case 1: .
If , (0,0) is a stable focus of Eq. (16), so (0,0,0,0) is also a stable focus of Eq. (18); if , because the solutions of Eq. (17), , are attracted to zero when , (0,0,0,0) is an asymptotically stable, but non-hyperbolic equilibrium point of Eq. (18); if , for Eq. (16), (0,0) is an unstable focus and there is a stable LCO , , so for Eq. (18) (0,0,0,0) is an unstable focus and there is a stable LCO as well. Supercritical Hopf bifurcation takes place in this case. In practice, because the amplitude of the LCO builds up from zero gradually as increases, it doesn’t bring immediately failure of the structure. So the flutter instability in this case is benign.
Case 2: .
For Eq. (18), if , (0,0,0,0) is a local stable focus, and there is an unstable LCO; if , (0,0,0,0) is an unstable non-hyperbolic equilibrium point; if , (0,0,0,0) is an unstable focus. Subcritical Hopf bifurcation occurs in this case. In practice, when , the response with any small initial conditions will move away from the trivial point immediately and the flutter is violent. So the flutter instability in this case is catastrophic.
4. The effect of the structural and aerodynamic parameters on character of the flutter instability
One concludes that the sign of decides the type of Hopf bifurcation or the character of flutter instability, i.e. benign or catastrophic. vs. as a function of the normalized nonlinear stiffness coefficient when are shown in Fig. 3 (curves are drawn pointwisely by using the first formula of Eq. (15)).
Fig. 3The relations of δ1 and flight Mach number M as a function of the structural nonlinearity B
The Fig. 3 indicates that the aerodynamic nonlinearity, representation , and the structural nonlinearity, representation , can affect the sign of .
For , and , the whole curves lie above the axis. So the flutter takes place always as a subcritical Hopf bifurcation for any value of , and the flutter instability is catastrophic. For and 50, the curves intersect the axis at the critical values. We define the critical value as , at where the sign of changes from negative to positive. This implies that with the increase of the flight Mach number , the flutter instability changes from benign to catastrophic at last. In other respect, the increase of hard structural nonlinearity () results in the increase of and the flight safety, and for hard structural nonlinearity the decrease of flight Mach number can increase the flight safety as well.
5. Examples, numerical simulations, and the effect of the parameters on character of the flutter instability
In the following, two cases with the flight Mach number and are considered and the flutter speeds are solved from Eq. (11) as for and for , respectively.
Take for example. In this case, the parameters in Eq. (16) were determined as:
Fig. 4The relation of δ1, B and γ for M=4
Fig. 5Supercritical Hopf bifurcation for γ=1.4, B=100, M=4: a) trajectory starting from x0T=(0.5,0.28,0,0) converges to O for V=14, b) the transient time history of h/b for V=14.3 and x0T=(0.0001,0.0001,0,0), c) the transient time history of α for V=14.3 and x0T=(0.0001,0.0001,0,0), d) trajectory starting from x0T=(0.5,0.28,0,0) converges to a LCO for V=14.3
The relation of , and is depicted in Fig. 4. The character of flutter instability is benign above the slope line, but it is catastrophic below the slope line. The decrease of isentropic gas coefficient results in the increase of the flight safety, and the increase of hard structural nonlinearity () also does.
Fig. 6Subcritical Hopf bifurcation for γ=1.4, B=2.5, M=4: a) trajectory starting from x0T=(0.04,0.0001,0,0) converges to a LCO immediately for V=14, b) trajectory starting from x0T=(0.0002,0.00015,0,0) converges to a LCO for V=14.3
Fig. 7Supercritical Hopf bifurcation for γ=1.4, B=21, M=10: a) trajectory starting from x0T=(0.8,0.5,0,0) converges to O for V=22.1, b) the transient time history of h/b for V=22.5 and x0T=(0.001,0.005,0,0), c) the transient time history of α for V=22.5 and x0T=(0.8,0.7,0,0)
To verify these theoretical results, numerical simulations using the Runge-Kutta algorithm to the Eq. (10) for were carried out. Notice that , . When and , . Supercritical Hopf bifurcation occurs, as shown in Fig. 5. Because the amplitude of the LCO is small no matter the initial distance is small or large, the flutter instability is benign. When and , . Subcritical Hopf bifurcation occurs, as shown in Fig. 6. The amplitudes of the two LCOs are large, so the flutter instability yields catastrophic failure of the structure.
For the case with the flight Mach number , the flutter speed , , and . The relation of , and is similar to what is depicted in Fig. 4. Numerical simulations for the Eq. (10) when were carried out. When and , , Fig. 7 indicates that the bifurcation is supercritical and the flutter instability is benign. Whilst when and Fig. 8 indicates that the subcritical Hopf bifurcation takes place and the flutter instability is catastrophic.
Fig. 8Subcritical Hopf bifurcation for γ=1.4, B=-10, M= 10: a) trajectory starting from x0T=(0.001,0,0,0) converges to O for V=22.1, b) trajectory starting from x0T=(0.00002,0.00001,0,0) moves away from O immediately for V=22.5
This paper deals with the flutter of a two-dimensional lifting surfaces in a supersonic flow field. This Eq. (4) undergoes only Hopf bifurcations in the neighborhood of the trivial equilibrium point if this system bifurcates at this point. The normal form and universal unfolding are calculated by using the maple program and the bifurcation theory. The supercritical Hopf bifurcation is benign, because it brings less flutter to the airfoil. Comparably the subcritical Hopf bifurcation is catastrophic, because it renders a harmful violent flutter of the airfoil.
In addition, we study the influence of the structural and aerodynamic parameters, representations and on the character of flutter instability. The influence of these parameters can be described as follows:
1) The soft structural nonlinearities () yield failure of the structure (i.e. catastrophic flutter), and yet the hard structural nonlinearities () perhaps not.
2) The increase of hard structural nonlinearity results in the increase of the airfoil safety, and for hard structural nonlinearity the decrease of the flight Mach number and isentropic gas coefficient also do.
Zhang Q. C., Liu H. Y., Ren A. D. The study of limit cycle flutter for airfoil with nonlinearity. Acta Aerodynamica Sinica, Vol. 22, Issue 3, 2004, p. 332-336, (in Chinese).
Yang X. L., Sethna P. R. Local and global bifurcations in parametrically excited vibrations of nearly square plates. International Journal of Non-Linear Mechanics, Vol. 26, Issue 2, 1991, p. 199-220.
Zhang Q. C., Liu H. Y., Ren A. D. Local bifurcation for airfoil with cubic nonlinearities. Journal of Tianjin University, Vol. 37, Issue 11, 2004, p. 970-974, (in Chinese).
Breitbach E. Effect of structural nonlinearities on aircraft vibration and flutter. The 45th Structures and Materials AGARD Panel Meeting, Agard report 665, 1977.
Librescu L. Elastostatics and Kinetics of Anisotropic and Heterogeneous Shell-Type Structures, Aeroelastic Stability of Anisotropic Multilayered Thin Panels. Mechanics of Elastic Stability. Noordhoff International Publishing, Leyden, 1975, p. 53-63, p. 106-158, 543-550.
Dowell E. H., Ilgamov M. Studies in Non-Linear Aeroelasticity. Springer-Verlag, New York, 1988.
Raghothama A., Narayanan S. Non-linear dynamics of a two-dimensional airfoil by incremental harmonic balance method. Journal of Sound and Vibration, Vol. 226, Issue 3, 1999, p. 493-517.
Shahrzad P., Mahzoon M. Limit cycle flutter of airfoils in steady and unsteady flows. Journal of Sound and Vibration, Vol. 256, Issue 2, 2002, p. 213-225.
Liu L., Wong Y. S., Lee B. H. K. Non-linear aeroelastic analysis using the point transformation method, Part 1: freeplay model. Journal of Sound and Vibration, Vol. 253, Issue 2, 2002, p. 447-469.
Zhao D. M., Zhang Q. C. Bifurcation and chaos analysis for aeroelastic airfoil with freeplay structural nonlinearity in pitch. Chinese Physics B, Vol. 19, Issue 3, 2010, p. 1-10.
Yang Y. R. KBM method of analyzing limit cycle flutter of a wing with an external store and comparison with a wind-tunnel test. Journal of Sound and Vibration, Vol. 187, Issue 2, 1995, p. 271-280.
Kim S. H., Lee I. Aeroelastic analysis of a flexible airfoil with a freeplay non-linearity. Journal of Sound and Vibration, Vol. 193, Issue 4, 1996, p. 823-846.
Price S. J., Alighanbari H., Lee B. H. K. The aeroelastic response of a two-dimensional airfoil with bilinear and cubic structural nonlinearities. Journal of Fluids and Structures, Vol. 9, Issue 2, 1995, p. 175-193.
Abbas L. K., Chen Q., Donnell K., Valentine D., Marzocca P. Numerical studies of a non-linear aeroelastic system with plunging and pitching freeplays in supersonic/hypersonic regimes. Aerospace Science and Technology, Vol. 11, Issue 5, 2007, p. 405-418.
Liu L., Wong Y. S., Lee B. H. K. Application of the center manifold theory in non-linear aeroelasticity. Journal of Sound and Vibration, Vol. 234, Issue 4, 2000, p. 641-659.
Ding Q., Wang D. L. The flutter of an airfoil with cubic structural and aerodynamic non-linearities. Aerospace Science and Technology, Vol. 10, Issue 5, 2006, p. 427-434.
Lee B. H. K., Liu L., Chung K. W. Airfoil motion in subsonic flow with strong cubic nonlinear restoring forces. Journal of Sound and Vibration, Vol. 281, Issue 3, 2005, p. 699-717.
Librescu L., Chiocchia G., Marzocca P. Implications of cubic physical/aerodynamic non-linearities on the character of the flutter instability boundary. International Journal of Non-Linear Mechanics, Vol. 38, Issue 2, 2003, p. 173-199.
Bi Q. S., Yu P. Symbolic computation of normal forms for semi-simple cases. Journal of Computational and Applied Mathematics, Vol. 102, Issue 2, 1999, p. 195-220.
Lu Q. S. Advanced Series in Nonlinear Science Bifurcation and Singularity. Shanghai Scientific and Technological Education Publishing House, Shanghai, 1995, (in Chinese).
Yu P. Computation of normal forms via a perturbation technique. Journal of Sound and Vibration, Vol. 211, Issue 1, 1998, p. 19-38.
About this article
This project was supported by the National Natural Science Foundation of China (Nos. 11172125 and 11202095) and the National Research Foundation for the Doctoral Program of Higher Education of China (20133218110025).