• International Journal of Aeronautical and Space Sciences
• /
• 제13권1호
• /
• pp.43-57
• /
• 2012

In a modern aircraft, there are many variations in its mass, stiffness, and aerodynamic characteristics. Recently, an analytical approach was proposed, and this approach uses the idea of uncertainty to find out the most critical flight flutter boundary due to the variations in such aerodynamic characteristics. An analytical method that has been suggested to predict robust stability is the mu method. We previously analyzed the robust flutter boundary by using the mu method, and in that study, aerodynamic variations in the Mach number, atmospheric density, and flight speed were taken into consideration. The authors' previous attempt and the results are currently quoted as varying Mach number mu analysis. In the author's previous method, when the initial flight conditions were located far from the nominal flutter boundary, conservative predictions were obtained. However, relationships among those aerodynamic parameters were not applied. Thus, the varying Mach number mu analysis results required validation. Using an optimization approach, the varying Mach number mu analysis was found out to be capable of capturing a reasonable robust flutter boundary, i.e., with a low percentage difference from boundaries that were obtained by optimization. Regarding the optimization approach, a discrete nominal flutter boundary is to be obtained in advance, and based on that boundary, an interpolated function was established. Thus, the optimization approach required more computational effort for a larger number of uncertainty variables. And, this produced results similar to those from the mu method which had lower computational complexity. Thus, during the estimation of robust aeroelastic stability, the mu method was regarded as more efficient than the optimization method was. The mu method predicts reasonable results when an initial condition is located near the nominal flutter boundary, but it does not consider the relationships that are among the aerodynamic parameters, and its predictions are not very accurate when the initial condition is located far from the nominal flutter boundary. In order to provide predictions that are more accurate, the relationships among the uncertainties should also be included in the mu method.

• 한국소음진동공학회논문집
• /
• 제22권4호
• /
• pp.393-400
• /
• 2012

This paper deals with a development of 2-DOF flutter experiment equipment which represents a 2-DOF typical section model. For a conventional 2-DOF flutter experiment equipment, it is hard to observe flutter boundary clearly due to the complexity of the experiment equipment. To refine our flutter experiment equipment system, a compliant mechanism based torsional spring is used. Well-designed extruded aluminum pipe works as a torsional spring. SolidWorks and ANSYS are used for modeling, analysis and design of the torsional spring. With this designed torsional spring, the 2-DOF flutter experiment equipment is developed and wind tunnel tests are performed. Clear flutter boundary which is estimated by classical flutter analysis is observed in the experiments.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2012년도 춘계학술대회 논문집
• /
• pp.429-434
• /
• 2012

This paper deals with a development of 2-DOF flutter experiment equipment which represents a 2-DOF typical section model. For a conventional 2-DOF flutter experiment equipment, it is hard to observe flutter boundary clearly due to the complexity of the experiment equipment. To refine our flutter experiment equipment system, a compliant mechanism based torsional spring is used. Well-designed extruded aluminum pipe works as a torsional spring. SolidWorks and ANSYS are used for modeling, analysis and design of the torsional spring. With this designed torsional spring, the 2-DOF flutter experiment equipment is developed and wind tunnel tests are performed. Clear flutter boundary which is estimated by classical flutter analysis is observed in the experiments.

• Advances in aircraft and spacecraft science
• /
• 제8권4호
• /
• pp.345-372
• /
• 2021

The analysis of nonlinear vibrations, buckling, post-buckling, flutter boundary determination and post-flutter behavior of a homogeneous curved plate assuming cylindrical bending is conducted in this article. Other assumptions include simply-supported boundary conditions, supersonic aerodynamic flow at the top of the plate, constant pressure conditions below the plate, non-viscous flow model (using first- and third-order piston theory), nonlinear structural model with large deformations, and application of mechanical and thermal loads on the curved plate. The analysis is performed with constant environmental indicators (flow density, heat, Reynolds number and Mach number). The material properties (i.e., coefficient of thermal expansion and modulus of elasticity) are temperature-dependent. The equations are derived using the principle of virtual displacement. Furthermore, based on the definitions of virtual work, the potential and kinetic energy of the final relations in the integral form, and the governing nonlinear differential equations are obtained after fractional integration. This problem is solved using two approaches. The frequency analysis and flutter are studied in the first approach by transferring the handle of ordinary differential equations to the state space, calculating the system Jacobin matrix and analyzing the eigenvalue to determine the instability conditions. The second approach discusses the nonlinear frequency analysis and nonlinear flutter using the semi-analytical solution of governing differential equations based on the weighted residual method. The partial differential equations are converted to ordinary differential equations, after which they are solved based on the Runge-Kutta fourth- and fifth-order methods. The comparison between the results of frequency and flutter analysis of curved plate is linearly and nonlinearly performed for the first time. The results show that the plate curvature has a profound impact on the instability boundary of the plate under supersonic aerodynamic loading. The flutter boundary decreases with growing thermal load and increases with growing curvature.

• International Journal of Aeronautical and Space Sciences
• /
• 제12권4호
• /
• pp.305-317
• /
• 2011

This paper presents an overview on flutter boundary prediction in tests which is principally based on a system stability measure, named Jury's stability criterion, defined in the discrete-time domain, accompanied with the use of autoregressive moving-average (AR-MA) representation of a sampled sequence of wing responses excited by continuous air turbulences. Stability parameters applicable to two-, three- and multi-mode systems, that is, the flutter margin for discrete-time systems derived from Jury's criterion are also described. Actual applications of these measures to flutter tests performed in subsonic, transonic and supersonic wind tunnels, not only stationary flutter tests but also a nonstationary one in which the dynamic pressure increased in a fixed rate, are presented. An extension of the concept of nonstationary process approach to an analysis of flutter prediction of a morphing wing for which the instability takes place during the process of structural morphing will also be mentioned. Another extension of analytical approach to a multi-mode aeroelastic system is presented, too. Comparisons between the prediction based on the digital techniques mentioned above and the traditional damping method are given. A future possible application of the system stability approach to flight test will be finally discussed.

• Wind and Structures
• /
• 제1권2호
• /
• pp.175-181
• /
• 1998

Torsional flutter occurs at 2D rectangular cylinders with side ratios B/D smaller than about 8 or 10. On the other hand, slender cylinders indicate the occurrence of coupled flutter, which means the coupled derivatives of slender cylinders have more significant role for flutter instability than that of bluffer ones. In this paper, based upon so called "Step-by-step analysis", it is clarified the coupled derivatives stabilize torsional flutter instability of bluffer cylinders (e.x. B/D=5), while they destabilize torsional flutter or coupled flutter instabilities of mores slender cylinders. The boundary of them exists between B/D=5 and 8.

• 한국군사과학기술학회지
• /
• 제12권6호
• /
• pp.704-710
• /
• 2009

In this study, a comparison study of flutter analysis for the AGARD 445.6 wing with wind turnnel test data has been conducted in the subsonic, transonic and supersonic flow regions. Nonlinear aeroelastic using FSIPRO3D which is a generalized user-friendly fluid-structure analyses have been conducted for a 3D wing configuration considering shockwave and turbulent viscosity effects. The developed fluid-structure coupled analysis system is applied for aeroelastic computations combining computational structure dynamics(CSD), finite element method(FEM) and computations fluid dynamics(CFD) in the time domain. MSC/NASTRAN is used for the vibration analysis of a wing model, and then the result is applied to the FSIPRO3D module. the results for dynamic aeroelastic response using different turbulent models are presented for several Mach numbers. Calculated flutter boundary are compared with the wind-tunnel experimental and the results show very good agreements.

• 소음진동
• /
• 제7권1호
• /
• pp.91-97
• /
• 1997

An analysis is presented on the stability of an elastic cantilever column having the elastic restraints at its free end, carrying an added tip mass, and subjected to uniformly distributed follower forces. The elastic restraints are formed by both a translational spring and a rotatory spring. For this purpose, the governing equations and boundary conditions are derived by using Hamilton's principle, and the critical flutter loads and frequencies are obtained from the numerical evaluation of the eigenvalue functions of this elastic system. The added tip mass increases as a whole the critical flutter load of the elastic cantilever column, but the presence of its moment of inertia of mass has a destabilizing effect. The existence of the translational and rotatory springs at the free end increases the critical flutter load of the elastic cantilever column. Nevertheless, their effects on the critical flutter load are not uniform because of their coupling. The translational spring restraining the free end of the cantilever column decreases the critical flutter load by coupling with a large value of tip mass, while by coupling with the moment of inertia of tip pass its effect on the critical flutter load is contrary. The rotatory spring restraining the free end of the cantilever column increases the critical flutter load by coupling with the tip mass, but decreases it by coupling with the moment of inertia of the tip mass.

• 한국항공우주학회지
• /
• 제34권1호
• /
• pp.18-23
• /
• 2006

본 연구에서는 Gurney 플랩이 달린 NACA 0012 익형에 대한 플러터 해석을 시간 영역에서 수행하였다. 2차원 비정상 압축성 Navier-Stokes 방정식과 Lagrange 방정식으로 부터 유도한 2계 자유도 plunge & pitch 모델을 지배방정식으로 하여 연성 결합 기법을 통해 플러터 해석을 수행하였다. 계산 결과 Gurney 플랩을 장착할 경우 NACA 0012에 비해 플러터가 발생하는 속도가 낮아졌고, 마하수가 0.85보다 작은 영역에서는 Gurney 플랩의 플러터 경계 곡선은 안전 여유를 상회하는 영역에 위치하였다. 그러나, 마하수가 0.85에서 0.9사이일 경우에는 안전 여유에 근접하게 되므로 이러한 운용 영역에서는 Gurney 플랩의 사용에 주의를 요한다.

• Structural Monitoring and Maintenance
• /
• 제4권2호
• /
• pp.115-131
• /
• 2017

Classical flutter of wind turbine blades indicates a type of aeroelastic instability with fully attached boundary layer where a torsional blade mode couples to a flapwise bending mode, resulting in a mutual rapid growth of the amplitudes. In this paper the monitoring problem of onset of flutter is investigated from a detection point of view. The criterion is stated in terms of the exceeding of a defined envelope process of a specific maximum torsional vibration threshold. At a certain instant of time, a limited part of the previously measured torsional vibration signal at the tip of blade is decomposed through the Empirical Mode Decomposition (EMD) method, and the 1st Intrinsic Mode Function (IMF) is assumed to represent the response in the flutter mode. Next, an envelope time series of the indicated modal response is obtained in terms of a Hilbert transform. Finally, a flutter onset criterion is proposed, based on the indicated envelope process. The proposed online flutter monitoring method provided a practical and direct way to detect onset of flutter during operation. The algorithm has been illustrated by a 907-DOFs aeroelastic model for wind turbines, where the tower and the drive train is modelled by 7 DOFs, and each blade by means of 50 3-D Bernoulli-Euler beam elements.