Research Article  Open Access
I. P. G. Sopan Rahtika, I. N. G. Wardana, A. A. Sonief, E. Siswanto, "Experimental Investigation on Flutter Similitude of ThinFlat Plates", Advances in Acoustics and Vibration, vol. 2017, Article ID 7091425, 8 pages, 2017. https://doi.org/10.1155/2017/7091425
Experimental Investigation on Flutter Similitude of ThinFlat Plates
Abstract
This paper shows the experimental results of the flutter speed of thinflat plates with free leading edge in axial flow as a function of plates’ geometry, fluid densities, and viscosities, as well as natural frequencies of the plates. The experiment was developed based on similitude theory using dimensional analysis and Buckingham Pi Theorem. Dimensional analysis generates four dimensionless numbers. Experiment was conducted by placing the thinflat plates in a laminar flow wind tunnel in order to obtain the relationship among those dimensionless numbers. The flutter speed was measured by varying the flow velocity until the instability occurred. The dimensional analysis gives a map of the flutter Reynolds number as a function of a new type of dimensionless number that is hereby called flutter fluid structure interaction number, thicknesstolength, and aspect ratios as the correcting factors. This map is a very useful tool for predicting the flutter speed of thinflat plates in general. This investigation found that the flutter Reynolds number is very high at the region of high flutter fluid structure and thicknesstolength ratios numbers; however, it is very sensitive to the change of those two dimensionless numbers. The sensitivity is higher at lower aspect ratio.
1. Introduction
Flutter is a potentially damaging dynamic aeroelastic phenomenon where aerodynamic forces with the natural modes of vibration cause a periodic motion of a structure going unstable. In a certain fluid structure interaction, the aerodynamic forces serve as input energy to the structural vibration. When the system is not damped by the aerodynamic damping, the vibration amplitude will increase and eventually will lead to structural failure. Flutter can occur in a variety of structures, such as in aircraft wings and turbine blades or even on the bridge. Most previous studies on flutter were focused on the prediction of the flutter speed using numerical methods. Very few literatures are available that discuss flutter using experimental similitude method as a tool to predict flutter phenomenon, despite the fact that similitudes offer cost savings in the investigation of the fluid flow phenomena. This scarcity was one that motivated this study.
The flutter phenomenon observed in this research focused on the flutter of thinflat plates with free leading edge in axial flow. There were numerous studies which have been done on the flutter of the thinflat plates due to the canonical characteristics of the problem, such as researches by Chad Gibbs et al. [1], Tang et al. [2], Tang and Païdoussis [3], Tang and Dowell [4], Howell et al. [5], and Zhao et al. [6].
Chad Gibbs et al. [1] performed experimental and theoretical work on the flutter of a flat plate with a fixed leading edge using the threedimensional vortexlattice method. Gibbs’ report comprehensively described the characteristics of plate flutter as a function of the mass ratio and the aspect ratio. Tang and Païdoussis [3] discussed the flutter of two flat plates which were positioned in parallel with the axial direction of the fluid flow.
The work of Tang and Dowell [4] was about nonlinear flutter and Limit Cycle Oscillation (LCO) of twodimensional panels in low subsonic flow and later was extended to threedimensional panels by Tang et al. [2]. The dynamics of the system was built to produce nonlinear models. Zhao et al. [6] discussed both theory and experiment of the flutter of a flat plate. The discussion emphasizes nonlinear analysis to produce Poincaré maps. The nonlinear analysis results were compared to experimental results. Howell et al. [5] used fluid flow interactions to discuss the nature of a cantilevered plate with ideal flow.
The aeroelastic instability of a flexible plate has been investigated using weakly nonlinear analyses by Eloy et al. [7]. Later, a deeper investigation was focused on the origin of the instability hysteresis [8].
Despite the detrimental effect of flutter, there are new research trends to utilize flutter for wind harvesting. The utilization of flutter phenomena for energy harvesting has been explored by Doaré and Michelin [9], Makihara and Shimose [10], and Dunmon et al. [11]. The other works on energy harvesting using a slender structure in the wake of a bluff body were also conducted by Allen and Smits [12]. As a matter of fact, this research is oriented as a theoretical base to explore this application further.
This research chose to explore the free leading edge instead of a fixed leading edge because the chosen configuration experiences flutter at a lower wind speed. A free leading edge plate will experience flutter in its first mode shape as shown later on the experimental results, while the fixed leading edge will experience flutter at the second mode [1, 5, 7]. This advantage of free leading edge is the reason of the selection of this configuration in this research.
Furthermore, the utilization of thinflat plates for wind harvester requires a map of flutter speed as a function of plates’ geometry for design optimization. One significant part of this research is to generate a flutter speed map for wind harvester design optimization. The authors currently have an ongoing development of a micropower generator utilizing flutter of a free leading edge configuration. This is new expansion of flutter energy conversion utilization into the field of microelectromechanical systems. This paper puts a benchmark for this research.
Flutter phenomena of plates continuously have been observed in the last several years. More recent studies on plate or panel flutter were also done by Fernandes and Mirzaeisefat [13], CunhaFilho et al. [14], Peng and DeSmidt [15], and Yaman [16].
Similitude theory has been well developed and widely used in the field of fluid dynamic. The use of similitude theory continues to grow in the field of structural vibration, for example, Torkamani et al. [17], and acoustic, for example, De Rosa et al. [18]. A work has also been done on structural similitude for flutter of composite plates by Yazdi and Rezaeepazhand [19].
In this study, a testing method is developed using dimensional similitude analysis based on the Buckingham Pi Theorem. Dimensional analysis generated four dimensionless numbers. Experiment was conducted on thinflat plate which is placed in a wind tunnel. The problem in this research was to formulate relationships between the flutter speed and the affecting parameters. Flutter speed is defined as the velocity of fluid flow at which flutter started happening. The plates’ parameters that affect the flutter speed that were taken into account in this experiment were their length, width, thickness, and natural frequencies, while the fluid’s parameters (in this case, fluid is air) that were taken into account were its density and viscosity. By using the similitude principle, the experimental results were used to obtain the relationship between the dimensionless numbers , , , and .
Later on, it has been observed during the analysis that the first dimensionless number is the Reynolds number measured during flutter and then called flutter Reynolds number. The second dimensionless number is called fluid structure interaction number since it contains the interacting forces. The third and fourth numbers are the thicknesstolength and aspect ratios, consecutively. This study finally discovered the relationship of the flutter Reynolds number as a function of the fluidstructure interaction number and the geometric ratios.
2. Theory of Flutter and Similitude
Theoretically, this research was about implementing the similitude theory to observe the flutter phenomenon of thinflat plates which then was used to generate a map for predicting the flutter speed. The conceptual theory of this research was developed from the theory of flutter and the similitude theory.
2.1. Flutter
Theoretical work on flutter has been recognized as early as 1878 by Rayleigh [20]. However, the practical work on flutter was later on reported by Theodorsen in 1934. Theodorsen explained the flutter mechanism theoretically and experimentally of the aircraft wing and also the combination of wingailerontab [21–23].
Flutter is an unstable fluid and structure interaction. The dynamics of a thinflat plate’s structure is stable by itself. However, when it is placed in a moving air, the aerodynamic forces shift the stability of the system. The system will become unstable when the air speed reaches a certain speed. This speed is called the flutter speed.
The dynamics of a thinflat plate can be modeled as a matrix equation of motion where the plate’s structural dynamics is subjected to aerodynamic force . In this case, are the element nodal displacement, velocity, and acceleration vectors, consecutively. , , and are the mass, damping, and stiffness matrices.
The aerodynamic force is nonlinear in nature. For the purpose of predicting the flutter speed, is often linearized to beThe denotation “” on and represents the aerodynamic contribution to the damping and stiffness matrices.
Substituting the linearized into the full plate’s aeroelastic equation of motion yields the plate’s linearized equation of motion:Equation (3) shows how the flutter can occur. The flutter will occur if (3) is unstable. The eigenvalue analysis of this linearized equation can give the value of the flutter speed of the plate.
2.2. Similitude Requirements for Modeling in Fluid Mechanics
The similitude requirements for the fluid dynamic problems are already well developed. Similitude is usually used in analyzing fluid dynamic problems such as lift and drag forces. In this study, the similitude method is used for analyzing the fluid structure interaction.
Similitude deals with the similarity of the actual system with its labscaled model or prototype. The similitude requirements for fluid dynamic problems have been well defined by Wolowicz et al. [24]. Similarity in geometric configuration is a fundamental requirement. Prototype and actual objects have to be geometrically congruent.
Another requirement is kinematic similarity. Two flows are kinematically similar if the associated velocities at the same point are related to a constant factor in the direction and magnitude [24]. This means that two streams are equal in their kinematics when streamline pattern associated with a constant factor.
Further requirements that must be met are the dynamics similarity. Two flows are dynamically similar when the associated forces at the same point are related to a constant factor in the direction and magnitude.
For the purpose of this research, an extra similitude requirement should be considered. Flutter problem is a moving boundary problem. The shape of the boundary varies with time. Ideally, in a steady case, the change is periodic or constant in frequency spectrum. Flutter will occur in a certain shape that is associated with the structural natural mode. Due to this reason, it is required to add a fourth similarity requirement when similitude theory is implemented to flutter problem. This fourth requirement is that the mode shape of the actual and the prototype must be congruent.
2.3. Buckingham Pi Theorem
The development of a similitude method relays very much on the Buckingham Pi Theorem. In this section, the Beckingham Pi Theorem is recalled as the base for the dimensional analysis.
For every physical phenomenon dependent parameter which can be expressed by a function of independent parameters, we can express the relationship between the parameters to formwhich is dependent parameters and is the independent parameters. In mathematics the above functional relationship can be expressed by an equivalent functionwhere is an unspecified function, different from .
Buckingham Pi Theorem states that [25] if there are parameters in the functionthen parameters can be grouped into different dimensionless ratio, or parameters, which can be expressed in the form of the functionor
Theorem does not predict the form of functional s or . The form of the functional relationship between dimensionless independent parameters must be determined experimentally.
2.4. Group for Flutter Similitude ThinFlat Plate
In this experimental study, the dependent variable is the flutter speed . Measurements were performed in SI units; thus, unit is m/s. Variable definitions for the plate’s dimensions can be referred to Figure 1. Independent parameters that are expected to affect the value of are = plate length (m), = plate thickness, = plate width, μ = fluid viscosity (N·s/m^{2} = kg/(m·s)), = fluid mass density (kg/m^{3}), = natural frequency of the thinflat plate (rad/s).
Relationship between the variables studied can be expressed by the symbolic function as follows:
Based on this function then . Primary dimensional dimension used in this analysis is the mass, length, and time. The selected repeating parameters are , , . So = 3. Analysis results obtained two s that are
The formulation of the dimensional analysis on (12) will work for geometrically similar plates. For the cases of thinflat plates, the thicknesses of the plates are relatively much smaller than the other geometric parameters and the variations in the thickness will not affect the geometric similarity of the flow significantly. Based on this reasoning—also approved by the experimental data—it is reasonable to do similitude analysis of the plates with various thicknesses but having the same aspect ratios. Allowing the variations in the plates’ thickness requires the thickness () to be included in the dimensional analysis. An extra dimensionless number is introduced to the analysis with the addition of thickness () to the dimensional analysis with the value:
Hence, the experiment was used to obtain the relationship of as a function of and .or
Consider the control volume of the fluidstructure interaction as shown in Figure 2. It is shown that the system is geometrically congruent from the top view projection. Then, it is possible to include the plates with different aspect ratio to analysis by considering the width to dimensional analysis. Inclusion of all plates with different aspect ratio to the analysis then introduces a fourth dimensionless number:
The form of first dimensionless number can be viewed as a form of comparison between inertial forces to the viscous forces of the fluid. It is a kind of Reynolds number. Since it is calculated during the flutter speed, it can be called flutter Reynolds number. Dimensionless number contains the inertia forces of both fluid and structure as well as the structural elastic forces and fluid viscous forces. Hence, represents the fluidstructural interaction of the system.
Since the experiments allow the variation in the thickness and aspect ratios free from the requirement of structural similarity, the thickness factor and aspect ratio have to be included in the empirical equation (18).
2.5. Structural Analysis
The experiment shows the flutter of the plates in this free leading edge configuration occurring in their first natural frequencies. The calculation of the second dimensionless number requires the value of the plate’s natural frequency. The mode shapes and the natural frequencies of the plates can be approximated using EulerBernoulli beam model as in (19). Using the beam approximation for a plate, the displacement of the plate is restricted to the transverse direction with its value as a function of the distance from the end of the beam and time : with
The natural frequency of a clampedfree beam [26] can be derived from partial differential equation (19) to become the following formula:
The calculated natural frequency from (21) is used to calculate the dimensionless number . In (20) and (21), is plate’s Young’s modulus of elasticity; is plate’s mass density; is plate’s crosssectional area; is second moment of plate’s crosssectional area. The value of constant depends on the boundary conditions and the length of the beam. For a cantilevered beam, the value of is 1.8751 for the first natural frequency.
3. Experimental Setup
The purpose of the experiment was to determine the function as in (18) that describes the relationship of the three dimensionless numbers. Therefore, the test data should be taken with a plate that has a variety of different length and different natural frequencies. Form of the function will be obtained through multivariable nonlinear regression analysis. All plates tested were copper plates with Young’s modulus of elasticity of 110 GPa and mass density of 8960 kg/m^{3}.
Experimental setup can be seen in Figure 2. The flutter speed of every single thinflat plate was measured one by one. The plate was placed on the wind tunnel that has maximum speed of 40 m/s. Each plate was positioned in vertical direction in the width direction. The trailing edge of the plate was clamped and the leading edge was set free. Air was blown by a blower into the test section. Air velocity in the wind tunnel was set by regulating the air intake into the blower. The flow velocity was measured using a pitot tube.
For each plate, the testing is done by varying the air velocity from zero to flutter speed. The image of the plate motion was captured by high speed camera that was set on 210 framepersecond speed as shown in Figure 3. The flutter speed determined from that image is the air speed when the plate was stopped at the infinite deflection. The results of measurements were presented in the form of the flutter speed as a function of plates’ dimensions.
4. Result and Discussion
The experimental results are tabulated in two ways: the first way is the tabulation without dimensional analysis as in Table 1 and the second one is the tabulation with dimensional analysis as in Table 2. Table 1 is divided into three subtables according to the plates’ aspect ratio. The list of flutter speed for various dimensions was obtained from the experiment and tabulated on Table 1. The natural frequency of each plate was also calculated using (21). The data from Table 1 is plotted on Figure 4 in the form of flutter speeds verses the lengths of the plates. In this graph, the data are also classified according to the plates’ aspect ratio: green dots for aspect ratio of 1 : 6; yellow ones for that of 1 : 5; and red for that of 1 : 3.
(a) Aspect ratio 1 : 6  
 
(b) Aspect ratio 1 : 5  
 
(c) Aspect ratio 1 : 3  

(a) = 1 : 6  
 
(b) = 1 : 5  
 
(c) = 1 : 3  

Not all physical trends can be observed in Figure 4. One trend can be observed from Figure 4: that is, a longer plate has a lower flutter speed. The length of the plate relates to the area of the interaction surface between structure and the fluid. A longer plate has a larger area of interaction surface so that it is easier to transfer energy from fluid; thus, the flutter will occur at lower speed.
Similitude analysis on flutter speed offers a comprehensive view of all data trends. The values of dimensionless variable are calculated using the data from Table 1 and tabulated on Table 2 using (10), (11), (13), and (16). Table 2 is also divided into three subtables according to the value of . The calculated dimensionless numbers from Table 2 are plotted on 2dimensional chart of versus as in Figure 5. In Figure 5, the data are also classified into three different aspect ratios: black dots are for = 1 : 6; blue dots are for = 1 : 5; and red dots are for = 1 : 3.
Some physical trends can be observed from Figure 5. The dimensionless number represents flutter speed and 2 represents the plate’s natural frequency. The relationship of the flutter speed with the plate’s natural frequency can be described using energy scheme. The higher the frequency of a vibrating plate, the higher its vibration energy. A plate with higher natural frequency will need a higher energy from the fluid to reach flutter.
Also from Figure 5, a plate with relative higher thickness with respect to its length will tend to have a higher flutter Reynolds number, because thicker plates tend to be stiffer and have a higher value of structural stiffness matrices. By thinking in the static stability scheme, to invoke the instability requires a higher speed with the stiffer plate.
The trends of the data are not clearly identified if the data are plotted on 2dimensional graph as in Figure 5 without the classification of their geometric ratios. This implies that the geometric ratios should be included in the analysis as the correcting factors as in (17). The data were apparently nonlinear and approximating the power trends. Then, a multivariable nonlinear regression analysis is conducted to obtain the form of function in (17). The result of the regression analysis is shown on (22) and (23). Substituting the values of as a function of , , and into (23) yields
The illustration of trends body of (22) requires a 4dimensional plot. The 4dimensional plot is projected as isoaspectratio surfaces on 3dimensional space as shown in Figure 6. The isoaspectratio surfaces are the relationship of as a function of and with constant value of .
Equation (23) is 4dimensional map of the tendency of the flutter speed for various plates’ dimensions. The projection of this map to 3dimensional space on Figure 6 shows the characteristic of the flutter Reynolds numbers as a function of flutter fluid structure interaction numbers , with the geometric ratios as the correcting factors. This map shows that flutter Reynolds number is very high at the region of high flutter fluid structure interaction number and thicknesstolength ratios number . However, it is very sensitive to the change of those two dimensionless numbers at the region. The sensitivity is higher at lower aspect ratio following the power trends. A lower aspect ratio means the plate is narrow and long (see (16)) so it has smaller interaction surface. The sensitivity of values of flutter Reynold number increases to the change of the values of fluid structure interaction number for the plate with lower aspect ratio. This happens as a result of decreasing contact surface for lower aspect ratio plates. Thus, the change of natural frequencies, which is also identical to the change of flutter fluid interaction number, has more dominant effect on flutter Reynolds number. On the contrary, at higher aspect ratio the contact surface is broader and it has dominant effect. As a result flutter Reynolds number is less sensitive to the change in natural frequency .
The map generated in this research is a very useful tool in optimizing the design of wind harvester given its range of operating wind speed. The map illustrates the flutter speed trends for various plates’ dimensions. A specific plate’s dimension will be optimum for a specific range of wind speed.
An experimental setup is often limited to a certain flutter Reynolds number due to the limited speed of the wind tunnel. Numerical investigations should have more freedom to exceed this limitation. The dimensional analysis and similitude method formulated in this paper can be used in a future research as method of presenting a result numerical investigation to generate a broader map, of course with some minor drawbacks of the numerical method limitations in representing the real case.
The similitude analysis presented in this paper is focused on the free leading edge plates in axial flow with application orientation for wind harvesting. However, there are also possibilities of using this method for other applications such as for the flutter speed analysis of airfoil of an aircraft wings or turbine blades.
5. Conclusions
Similitude method has been developed to generate an empirical map to predict the flutter speed of a thinflat plate. The experiment has showed that the trend of flutter speed has been traced in various geometric parameters. The relationship between the flutter speeds of thin flatplates in a variety of sizes that are geometrically congruent and the parameters that affect them such the plates’ width, thickness, length, fluid density, fluid viscosity, and their natural frequencies has been described.
This research also has found an empirical relationship between the flutter Reynolds number and the fluidstructure interaction number with the geometric ratios as correcting factors.
Departing from the testing of thinflat plate, this study provides a formulation in testing the similitude theory for a future study on flutter analysis of an airfoil that is used in the aircraft or turbomachinery industries. In addition, flutter similitude method can also be used to test the possibility of a detrimental effect of flutter on narrow and long bridges.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
References
 S. Chad Gibbs, I. Wang, and E. Dowell, “Theory and experiment for flutter of a rectangular plate with a fixed leading edge in threedimensional axial flow,” Journal of Fluids and Structures, vol. 34, pp. 68–83, 2012. View at: Publisher Site  Google Scholar
 D. M. Tang, H. Yamamoto, and E. H. Dowell, “Flutter and limit cycle oscillations of twodimensional panels in threedimensional axial flow,” Journal of Fluids and Structures, vol. 17, no. 2, pp. 225–242, 2003. View at: Publisher Site  Google Scholar
 L. Tang and M. P. Païdoussis, “The coupled dynamics of two cantilevered flexible plates in axial flow,” Journal of Sound and Vibration, vol. 323, no. 3–5, pp. 790–801, 2009. View at: Publisher Site  Google Scholar
 D. Tang and E. H. Dowell, “Limit cycle oscillations of twodimensional panels in low subsonic flow,” International Journal of NonLinear Mechanics, vol. 37, no. 7, pp. 1199–1209, 2002. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 R. M. Howell, A. D. Lucey, P. W. Carpenter, and M. W. Pitman, “Interaction between a cantileveredfree flexible plate and ideal flow,” Journal of Fluids and Structures, vol. 25, no. 3, pp. 544–566, 2009. View at: Publisher Site  Google Scholar
 W. Zhao, M. P. Païdoussis, L. Tang, M. Liu, and J. Jiang, “Theoretical and experimental investigations of the dynamics of cantilevered flexible plates subjected to axial flow,” Journal of Sound and Vibration, vol. 331, no. 3, pp. 575–587, 2012. View at: Publisher Site  Google Scholar
 C. Eloy, R. Lagrange, C. Souilliez, and L. Schouveiler, “Aeroelastic instability of cantilevered flexible plates in uniform flow,” Journal of Fluid Mechanics, vol. 611, pp. 97–106, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 C. Eloy, N. Kofman, and L. Schouveiler, “The origin of hysteresis in the flag instability,” Journal of Fluid Mechanics, vol. 691, pp. 583–593, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 O. Doaré and S. Michelin, “Piezoelectric coupling in energyharvesting fluttering flexible plates: linear stability analysis and conversion efficiency,” Journal of Fluids and Structures, vol. 27, no. 8, pp. 1357–1375, 2011. View at: Publisher Site  Google Scholar
 K. Makihara and S. Shimose, “Supersonic flutter utilization for effective energyharvesting based on piezoelectric switching control,” Smart Materials Research, vol. 2012, Article ID 181645, 10 pages, 2012. View at: Publisher Site  Google Scholar
 J. A. Dunnmon, S. C. Stanton, B. P. Mann, and E. H. Dowell, “Power extraction from aeroelastic limit cycle oscillations,” Journal of Fluids and Structures, vol. 27, no. 8, pp. 1182–1198, 2011. View at: Publisher Site  Google Scholar
 J. J. Allen and A. J. Smits, “Energy harvesting eel,” Journal of Fluids and Structures, vol. 15, no. 34, pp. 629–640, 2001. View at: Publisher Site  Google Scholar
 A. C. Fernandes and S. Mirzaeisefat, “Flow induced fluttering of a hinged vertical flat plate,” Ocean Engineering, vol. 95, pp. 134–142, 2015. View at: Publisher Site  Google Scholar
 A. G. CunhaFilho, A. M. G. de Lima, M. V. Donadon, and L. S. Leão, “Flutter suppression of plates using passive constrained viscoelastic layers,” Mechanical Systems and Signal Processing, vol. 79, pp. 99–111, 2016. View at: Publisher Site  Google Scholar
 M. Peng and H. A. DeSmidt, “Stability analysis of a flutter panel with axial excitations,” Advances in Acoustics and Vibration, vol. 2016, Article ID 7194764, 7 pages, 2016. View at: Publisher Site  Google Scholar
 K. Yaman, “Subsonic flutter of cantilever rectangular PC plate structure,” International Journal of Aerospace Engineering, vol. 2016, Article ID 9212364, 10 pages, 2016. View at: Publisher Site  Google Scholar
 S. Torkamani, H. M. Navazi, A. A. Jafari, and M. Bagheri, “Structural similitude in free vibration of orthogonally stiffened cylindrical shells,” ThinWalled Structures, vol. 47, no. 11, pp. 1316–1330, 2009. View at: Publisher Site  Google Scholar
 S. De Rosa, F. Franco, X. Li, and T. Polito, “A similitude for structural acoustic enclosures,” Mechanical Systems and Signal Processing, vol. 30, pp. 330–342, 2012. View at: Publisher Site  Google Scholar
 A. A. Yazdi and J. Rezaeepazhand, “Structural similitude for flutter of delaminated composite beamplates,” Composite Structures, vol. 93, no. 7, pp. 1918–1922, 2011. View at: Publisher Site  Google Scholar
 L. Rayleigh, “On the instability of jets,” Proceedings of the London Mathematical Society, vol. s110, no. 1, pp. 4–13, 1878. View at: Publisher Site  Google Scholar  MathSciNet
 T. Theodorsen, “General theory of aerodynamic instability and the mechanism of flutter,” NACA Technical Report 469, 1934. View at: Google Scholar
 T. Theodorsen and L. E. Garrick, “Mechanism of flutter: a theoretical and experimental investigation of the flutter problem,” NACA Technical Report 685, 1940. View at: Google Scholar
 T. Theodorsen and L. E. Garrick, “Nonstationary flow about awingaelerontab combination including aerodynamic balance,” NACA Technical Report 736, 1941. View at: Google Scholar
 C. H. Wolowicz, J. S. Bowman Jr., and W. P. Gilbert, “Similitude requirement and scaling relationships as applied to model testing,” NASA Technical Paper, 1979. View at: Google Scholar
 R. W. Fox and A. T. McDonald, Introduction to Fluid Mechanics, John Wiley & Son, 1994.
 D. J. Inman, Engineering Vibration, PrenticeHall, Upper Saddle River, NJ, USA, 1994.
Copyright
Copyright © 2017 I. P. G. Sopan Rahtika et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.