Flutter Analysis of Control Surface Free-play by Using a 3-DOF Airfoil Model 107 where the real parts of the eigenvalues are examined to identify the flutter speed. Depending on the value of eigenvalue λ, behavior of the system may exhibit oscillatory, dissipative, or dispersive behavior. If λ<0, the system is dynamically stable; if λ >0, it becomes dynamically unstable, leading to divergent oscillations. The point where λ = 0 defines the ‘stability boundary.’ In aeroelastic systems, this instability is known as ’flutter’, and the corresponding boundary is called the ‘flutter boundary’. In flutter cases, the solution typically shows oscillatory or dispersive behavior, both indicating instability, where even a small disturbance can lead to large structural displacements. Hence, the flutter speed is determined as the point where the real part of the eigenvalue transitions from negative to positive. In state space form, equation of motion of the system can be written as follows: A˙q+Bq=0 (12a) A= 0 M M C , B= −M 0 0 K , q= ˙x x (12b) Here, q assuming harmonic form for x, qcan be represented as follows q=qeλt, (13a) Substituting Eq. (13a) into Eq. (12a) the following eigenvalue problem is obtained λAq+Bq=0. (13b) Eq. (13b) is for the linear problem. For the nonlinear problem, DFM is used to find the nonlinear stiffness matrix. This method, as demonstrated by Iannelli et al. [13], allows for stability assessment without requiring detailed frequencydomain analysis. In their study, Iannelli et al. applied DFM to aeroelastic systems with free-play nonlinearities, addressing uncertainties using structured singular values and integral quadratic constraints. Their work effectively captured the dynamic behavior of aeroelastic systems, providing an efficient way to model and analyze nonlinearities such as free-play in control surfaces. In the presence of free-play nonlinearity, the nonlinear restoring torque acting on the control surface can be written as follows by using the describing function theory TN =ν • bβ, (14) where, ν is the describing function which is defined as ν = i πAZ 2π 0 TN(Asin(ψ))e−iψdψ. (15) TN is the nonlinear restoring torque, Ais the amplitude of oscillations, ψ =ωt and bβ is the complex representation of the flap displacement. It should be noted that, in general, ν is complex quantity; however, for free play nonlinearity, it is a real quantity, and it can be obtained as follows ν(A)= −2k2 π sin−1 δ A + δ Aq1− δ A 2 +k2 A>δ 0 A≤δ , (16) Here, Ais the flap displacement amplitude, k2 is the torsional stiffness of the hinge mechanism, and δ defines the gap in free-play at one side of the equilibrium resulting total gap of 2δ. Using DFM, the stiffness matrix of the system can be written as follows K=K−Ka+ 0 0 0 0 0 0 0 0 ν(β) . (17) Results The flutter characteristics of the system are influenced by several structural and geometric parameters, which affect both flutter speed and frequency. The interaction between these parameters and the aeroelastic response is critical for understanding the stability of the system. The results demonstrate how variations in stiffness, semi-chord length, and hinge line position shape susceptibility to flutter, providing insights into optimizing aeroelastic performance.
RkJQdWJsaXNoZXIy MTMzNzEzMQ==