Next Article in Journal
Effect of Bending Deformation on the Lateral Force of Spinning Projectiles with Large Aspect Ratio
Next Article in Special Issue
Study on the Aerothermoelastic Characteristics of a Body Flap Considering the Nozzle–Jet Interference
Previous Article in Journal
A Turbulent Inflow Generation Method for the LES of High Re Flow by Scaling Low Re Flow Data
Previous Article in Special Issue
Aeroelastic Response of Spinning Projectiles with Large Slenderness Ratio at Supersonic Speed
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Aeroelasticity of a Cantilever Wing Equipped with the Spanwise Morphing Trailing Edge Concept

1
Department of Aerospace Engineering, Khalifa University, Abu Dhabi P.O. Box 127788, United Arab Emirates
2
Advanced Research and Innovation Center (ARIC), Department of Aerospace Engineering, Khalifa University, Abu Dhabi P.O. Box 127788, United Arab Emirates
3
Faculty of Engineering, University of Nottingham, Nottingham NG7 2RD, UK
*
Author to whom correspondence should be addressed.
Aerospace 2023, 10(9), 809; https://doi.org/10.3390/aerospace10090809
Submission received: 1 June 2023 / Revised: 6 September 2023 / Accepted: 8 September 2023 / Published: 15 September 2023
(This article belongs to the Special Issue Applied Aeroelasticity and Fluid-Structure Interaction)

Abstract

:
This paper studies the aeroelastic behavior of a rectangular, cantilever wing equipped with the spanwise morphing trailing edge (SMTE) concept. The SMTE consists of multiple trailing edge flaps that allow controlling the spanwise camber distribution of a wing. The flaps are attached at the wing’s trailing edge using torsional springs. The Rayleigh–Ritz method is used to develop the equations of motion of the wing-flap system. The use of shape functions allows for representing the wing as an equivalent 2D airfoil with generalized coordinates that are defined at the wingtip. Strip theory, based on Theodorsen’s unsteady aerodynamic model, is used to compute the aerodynamic loads acting on the wing. A representative Padé approximation for Theodorsen’s function is utilized to model the aerodynamic behaviors in a state-space form allowing time-domain simulation and analysis. The model is validated using a rectangular cantilever wing and the data are available in the literature. A comprehensive parametric comparison study is conducted to assess the impact of flap stiffness on the aeroelastic boundary. In addition, the potential of the SMTE to provide load alleviation and flutter suppression is assessed for a wide range of flight conditions, using a discrete (1-cosine) gust. Finally, the implementation and validation of a controller for a wing with SMTE for gust load alleviation are studied and controller parameters are tuned for a specific gust model.

1. Introduction

In recent years, extensive research has been carried out to develop gust load alleviation (GLA) systems. The aerodynamic forces generated by GLA systems modify the overall forces in such a way as to alleviate the structural turbulence response. Moreover, minimizing the impact of gusts by deploying an active GLA system, typically through the use of conventional control surfaces, is an integral part of modern aircraft design [1]. Bernhammer et al. [2] presented an experimental aeroservoelastic investigation of a novel load alleviation concept using trailing edge flaps. In their model, the flaps were free-floating and mass underbalanced, such that they may become unstable at operational velocities unless suppressed by their control system (trailing edge tabs). They found that limit cycle oscillation could be reached either through structural limiters or by control actions of the trailing edge tabs. In the latter case, the amplitude of the limit cycle oscillation could be adjusted to the required energy output. An energy balance between harvested power and power consumption for actuators and sensing systems was made showing that the vibration energy of limit cycle oscillations could be used to keep the amplitude of the limit cycle constant, while the electric batteries, powering up the load alleviation system, were being charged. Wildschek et al. [3] investigated the gust load response of a large 750-passenger Blended Wing Body (BWB) airliner for the identification of sizing cases for the aircraft structure. They designed a GLA system and evaluated it by numerical simulations to allow for structural weight saving. They concluded that the structural weight saving was mainly limited by the allowed load factors considered in such optimization, the finite control authority of the actuators, and the achievable reliability of the GLA system itself. Huvelin et al. [4] conducted a numerical simulation for the gust response along with experimental results at ONERA to present a gust load alleviation study. Numerical comparisons were performed using various techniques of gust modeling and finally, an application example for gust load alleviation was presented. Bekemeyer et al. [5] investigated the gust load alleviation for an airfoil and a large aircraft configuration using computational fluid dynamics. The results of the study were then used in designing a simple gust controller aimed at negating gust-induced loads via deploying conventional control surfaces. Pusch [6] studied the effectiveness of using distributed flaps for wing gust load alleviation. It was concluded that it is possible to simplify the controller design and tuning and achieve the desired performance by using the proposed allocation method. Dias and Hubbard [7] presented a novel approach to tackle the gust load alleviation problem by controlling the shape of the lift distribution profile along the span. The unsteady aerodynamics of a finite wing featuring multiple trailing edge flaps were modeled using the unsteady vortex lattice method (UVLM), yielding a linear, time-variant, high-order-state-space model. By using an aerodynamic mode shape and associated eigenvalues, this distributed approach allows control of the loads at all the points. Pusch et al. [8] discussed the optimization of control surface layout for gust load alleviation by using a nonlinear model of a large-scale flexible aircraft with unsteady aerodynamics. Compared to a GLA, the system using the original aileron configuration exhibited a 9% performance improvement. Furthermore, a trade-off study was carried out which enabled a target-oriented balancing between individual load channels. The significance of aileron size and position on overall GLA performance was demonstrated, and consideration was proposed for the preliminary aircraft design process.
Recently, spanwise morphing trailing edge (SMTE) concepts have attracted a lot of interest due to their ability to alter the spanwise camber distribution. SMTE can be considered a morphing aircraft technology that allows altering the spanwise camber distribution for a variety of purposes including control authority, load alleviation, drag reduction, and flutter suppression. For more details on camber morphing and its applications, the reader is advised to check Barbarino et al. [9] and Ajaj et al. [10]. Several studies have assessed the performance benefits of SMTE. For instance, Pankonien and Inman [11,12] proposed a modular SMTE concept that locally varied the trailing edge camber of a wing functioning as a modular replacement for conventional control surfaces without altering the spar box. Utilizing alternating active sections of Macro Fiber Composites (MFCs) driving internal compliant mechanisms and inactive sections of elastomeric honeycombs, the SMTE concept eliminated geometric discontinuities associated with shape change; increasing aerodynamic performance. The modular morphing wing is shown in Figure 1. Force and deformation analysis of a morphing wing was performed with a smooth surface and independent spanwise-varying control surface. Wind tunnel tests at a flow speed of 10 m/s were also carried out using a hardware demonstrator. Results indicated that the developed spanwise morphing trailing edge led to excellent aerodynamic and structural performance.
Some recent analytical and experimental studies show the potential and feasibility of an active aeroelastic control system. The control surfaces are operated according to a control law that relates the motion of the controls to some measurements taken on the aircraft. Haghighat et al. [13] investigated the design optimization of an active load alleviation control system which was integrated with the design optimization of the aerodynamic shape and structural sizing of a UAV. Khalil and Fezans [1] discussed a methodology for designing combined feedback/feedforward GLA systems. The methodology applied to large aircraft and the simulation results showed the effectiveness of achieving the desired objectives while ensuring both design flexibility and control robustness and optimality. Capello et al. [14] discussed a comprehensive robust adaptive controller for gust load alleviation. They implemented and validated the proposed approach on subsonic aircraft for different mass flight conditions. Magar et al. [15] explored the use of origami to achieve camber morphing for vibration suppression and gust load alleviation of a typical wing section. Bourchak et al. [16] presented an optimum design of a PID controller for the adaptive torsion wing (ATW) using a genetic algorithm optimizer. In many papers, see (Kuznetsov et al. [17], Mozaffari-Jovin et al. [18], Vindigni et al. [19] and Munoz et al. [20]), various active wing-flutter suppression and control methods have been investigated. A morphing aircraft is usually associated with significant changes in the aerodynamic loads, structural/elastic properties, inertial properties, aeroelastic behavior, and flight dynamics and stability characteristics. This necessitates effective and robust control strategies to ensure certain stability and performance criteria are met during the morphing process. In addition, wing-shape changes require effective controllers to provide suitable actuation under various flight conditions and mission profiles. In summary, morphing can lead to complex time-varying nonlinear dynamical models with internal and external uncertainties [21,22]. These uncertainties and time-varying characteristics demand sophisticated control systems to confirm the stability and performance of the morphing wings.
It should be noted that several studies investigated the aeroelasticity of an airfoil with a flap. For example, Irani et al. [23] investigated limit cycle oscillations as well as nonlinear aeroelastic analysis of three degrees of freedom (dof) aeroelastic airfoil motion with cubic restoring moments in the pitch degree of freedom. The majority of studies in the literature have focused on developing SMTE concepts and assessing their performance characteristics; however, very little has been done on studying the aeroelasticity of SMTE [11,12]. SMTE concepts can have significant effects on the aeroelastic boundaries (divergence and flutter) and can be actively used to suppress some of the critical aeroelastic phenomena. This paper aims to fill this gap and conduct a comprehensive study on the aeroelasticity of SMTE. To achieve this, a rectangular, cantilever wing is equipped with an SMTE concept (consisting of three flaps). It is common for design studies to begin with low to medium-fidelity tools and move to higher-fidelity tools at later stages [24]. A low-fidelity aeroelastic model is developed using the Rayleigh–Ritz method coupled with strip theory aerodynamics. A parametric study is conducted to assess the impact of flap stiffness on the aeroelastic boundaries of the wing. Then, the ability of the SMTE to provide load alleviation is assessed and studied. Furthermore, a GLA controller is designed for SMTE and its effectiveness is assessed.

2. Aeroelastic Modeling

To simplify the analysis, a rectangular, cantilever wing equipped with three discrete flaps attached at the trailing edge is considered. Each flap is a rectangular segment connected to the wing using a torsional spring and the mechanical and geometric properties of each flap are uniform but might differ from one flap to another. The properties associated with each flap segment are listed in Table 1. The Rayleigh–Ritz method is used to develop the equation of motion for the wing with SMTE. The shape functions used here correspond to the uncoupled first bending and first torsion modes of a uniform cantilever beam. In addition, a shape function for each flap is used. This allows the wing with SMTE to be modeled as an equivalent airfoil with generalized coordinates that are defined at the wingtip. Figure 2 and Figure 3 show the layout of the wing with SMTE and the corresponding nomenclature used. Using standard notation, the mass of wing-flap per unit span is denoted by m , the plunging deflection is denoted by h, positive in the upward direction; α is the pitch angle about the elastic axis, positive when nose-up; and β i is the control surface angle of ith flap, positive when the trailing edge (TE) flap is moved down. The elastic axis is located at a distance a b from the mid-chord, where b is semi-chord, while the wing mass center is located at a distance x α b from the elastic axis. The axis of rotation for the ith control surface is located at a distance c i from the mid-chord for the ith wing segment, while the center of mass of the ith flap is located at a distance x β i from the hinge. All distances are positive when measured towards the TE of the airfoil. In Figure 2, K h , K α and K β i are the stiffness in plunge, pitch, and the ith flap, respectively.
The continuous, multi-degree of freedom, wing structure is modeled as two degrees of freedom system via the Rayleigh–Ritz method using shape functions. These shape functions correspond to the uncoupled first bending and first torsional modes of a uniform cantilever beam. The first bending shape function, f y , is given as:
f y = 0.5 ( cosh B n y cos ( B n y ) ) σ n ( sinh ( B n y ) sin ( B n y )
where y is the spanwise position measured from the wing root and
σ n = ( cosh B n l + cos ( B n l ) ) sinh ( B n l + sin ( B n l ) )
where l is the wing semi-span and
B n l = 1.875
the torsion shape function, ϕ ( y ) , is given as
ϕ ( y ) = sin π 2 l y
Similarly, to obtain the torsion shape functions of the flap and its boundary is given as,
Ψ i y = 1 ,   Ψ 1 y = 1 , Ψ 2 y = 0 , Ψ 3 y = 0 0 y l 1 Ψ 1 y = 0 , Ψ 2 y = 1 , Ψ 3 y = 0 l 1 < y l 1 + l 2 Ψ 1 y = 0 , Ψ 2 y = 0 , Ψ 3 y = 1 l 1 + l 2 < y l 1 + l 2 + l 3
Using the variable separation approach and shape functions, the plunge displacement, speed, and acceleration at any spanwise location (y) can be related to those of the wingtip (generalized coordinates) as:
h ( t , y ) = h t t f ( y ) h ˙ ( t , y ) = h t ˙ t f ( y ) h ¨ ( t , y ) = h t ¨ t f ( y )
The pitch displacement, speed, and acceleration at any spanwise location ( y ) and time instant can now be related to those of the wingtip (generalized coordinates) as:
α ( t , y ) = α t t ϕ ( y ) α ˙ ( t , y ) = α t ˙ t ϕ ( y ) α ¨ ( t , y ) = α t ¨ t ϕ ( y )
Similarly, the flap displacement, speed, and acceleration at any spanwise location ( y ) and time instant can now be related to those of the wingtip (generalized coordinates) as:
β i ( t , y ) = β t i t Ψ i ( y ) β i ˙ ( t , y ) = β t i ˙ t Ψ i ( y ) β i ¨ ( t , y ) = β t i ¨ t Ψ i ( y )
where h t t , α t t , and β t i t represent the generalized coordinates coinciding with the wingtip and subscript ‘i’ in Equation (3) represents the flap number.

2.1. Equations of Motion

The kinetic energy ( T ) and potential energy ( U ) of the system can be expressed as:
T = 1 2   m h ˙ t 2 0 l f 2 d y + 1 2   I α α ˙ t 2 0 l ϕ 2 d y S α h ˙ t α ˙ t 0 l f ϕ d y + i = 1 n { 1 2   I β i β ˙ t i 2 0 l Ψ i 2 d y S β i h ˙ t β ˙ t i 0 l f Ψ i d y + α ˙ t β ˙ t i I β i + b c i a S β i 0 l ϕ Ψ i d y }
U = 1 2 G J α t 2 0 l d d y 2 d y + 1 2 E I h t 2 0 l d 2 f d y 2 2 d y + i = 1 n 1 2 K β i β t i 2 0 l Ψ i 2 d y
where I α , and I β i are the mass moment of inertia of the wing and flap per unit span of the wing, respectively, S α is the static moments of the wing, S β i is the static moments of the flap, and n is the total number of flaps. The equations of motion of the system using Lagrangian mechanics can be obtained as:
d d t ( T U ) h ˙ t ( T U ) h t = L
d d t ( T U ) α ˙ t ( T U ) α t = M α
d d t ( T U ) β ˙ t i ( T U ) β t i = M β i
where L and M α are the generalized lift and pitching moment of the wing, respectively, and M β i is the moment about flap hinge points for the ith flap. They can be expressed as:
L = 0 l L h y d y ,   M α = 0 l M α ϕ y d y     a n d     M β i = 0 l M β i Ψ i y d y
where L , M α , and M β i are the generalized lift, pitching moment, and moment about flap hinge points per unit length, respectively.

2.2. Aerodynamic Loads

To compute the unsteady aerodynamic loads, the strip theory, based on Theodorsen’s unsteady aerodynamic model, is used. Theodorsen’s unsteady aerodynamics model has a circulatory component to account for the effect of the wake on the airfoil and it contains the main damping and stiffness terms and a non-circulatory component to account for the acceleration of the fluid surrounding the airfoil [25]. L is the unsteady lift per unit span of the wing,   M α is the pitching moment around the elastic axis per unit span of the wing, and M β i   is the moment about flap segments hinge point of the i t h   flap per unit span of the wing. The expressions for L , M α , and M β i   are taken from NACA Report No.496 [25]. It should be noted that the expression of M β i   varies from one flap to another depending on the properties of the flap.
L = π ρ b 2 h ¨ + v α ˙ a b α ¨ v T 4 i π β i ˙ b T 1 i π β i ¨ + 2 π ρ v b C ( k )   h ˙ + v α + b 1 2 a α ˙ + v T 10 i π β i + b T 11 i 2 π β i ˙  
M α = π ρ b 2 ( a b h ¨ + a 1 2 b v α ˙ b 2 1 8 + a 2 α ¨ T 4 i + T 10 i β i T 1 i T 8 i c a T 4 i + 1 2 T 11 i v b π β i ˙ + ( T 7 i + c a T 1 i ) b 2 π β i ¨ ) + 2 π ρ v b 2 a + 1 2 C ( k )   h ˙ + v α + b 1 2 a α ˙ + v T 10 i π β i + b T 11 i 2 π β i ¨  
M β i = π ρ b 2 ( T 1 i b π h ¨ 2 T 9 i T 1 i T 4 i a 1 2 v b π α ˙ 2 T 13 i b 2 π α ¨ T 5 i T 4 i T 10 i v π 2 β i + T 4 i T 11 i v b 2 π 2 β i ˙ + T 3 i b π 2 β i ¨ ) ρ v b 2 T 12 i C ( k )   h ˙ + v α + b 1 2 a α ˙ + v T 10 i π β i + b T 11 i 2 π β i ¨  
Theodorsen’s unsteady aerodynamics model has a frequency-dependent term, Theodorsen’s transfer function which accounts for attenuation of lift amplitude and phase lag in the lift response due to sinusoidal motion. In this paper, unsteady lift per unit span and pitching moments per unit span are expressed in the time domain. Therefore, a Padé approximation developed by Brunton and Rowley [26] for Theodorsen’s transfer function is used. The approximate transfer function C s in the Laplace domain becomes
C s 0.5177 a ^ 2 s 2 + 0.2752 a ^ s + 0.01576 a ^ 2 s 2 + 0.3414 a ^ s + 0.01582
where
a ^ = c 2 V
The equivalent lift force becomes,
L = 0.0985 a ^ u ˙ + 0.0076 a ^ 2 u + 0.5177 B ^ h ˙ t 0 l f 2 d y A ^ h ¨ t 0 l f 2 d y + 0.5177 B ^ V α t 0 l f ϕ d y + A ^ V + 0.5177 B ^ V b ( 1 2 a ) α t ˙ 0 l f ϕ d y A ^ a b α ¨ t 0 l f ϕ d y + β t i 0 l 0.5177 B ^ V π   T 10 i f Ψ i d y + β t i ˙ 0 l 0.5177 B ^ b 2 π T 11 i A ^ V π T 4 i f Ψ i d y β t i ¨ 0 l A ^ b π T 1 i h f d y
and the pitching moment is given by
M α = 0.0985 a ^ b 1 2 + a v ˙ + 0.0076 a ^ 2 b 1 2 + a v 0.5177 B ^ b 1 2 + a h ˙ t 0 l f ϕ d y A ^ b a h ¨ t 0 l f ϕ d y + 0.5177 b 1 2 + a B ^ V α t 0 l ϕ 2 d y b a 1 2 A ^ V + 0.5177 B ^ V b 2 1 2 + a b 1 2 a α t ˙ 0 l ϕ 2 d y A ^ b 2 ( 1 8 + a 2 ) α ¨ t 0 l ϕ 2 d y + β t i 0 l 0.5177 b ( 1 2 + a ) B ^ V π   T 10 i A ^ V 2 π ( T 4 i + T 10 i ) ϕ Ψ i d y + β t i ˙ 0 l 0.5177 B ^ b 2 ( 1 2 + a ) 2 π T 11 i A ^ V b π ( T 1 i T 8 i c a T 4 i + T 11 i 2 ) ϕ Ψ i d y β t i ¨ 0 l A ^ b 2 π ( T 7 i + c a T 1 i ) ϕ Ψ i d y
Similarly, the moment on the flap is given by
M β i = 0.0985 a ^ b T 12 i 2 π w i ˙ + 0.0076 a ^ 2 b T 12 i 2 π v 0.5177 B ^ b h ˙ t 0 l T 12 i 2 π f Ψ i d y A ^ b h ¨ t π 0 l T 1 i f Ψ i d y + 0.5177 b B ^ V α t 0 l T 12 i 2 π ϕ Ψ i d y + α t ˙ 0 l 0.5177 B ^ V b 2 1 2 a T 12 i 2 π A ^ V b π 2 T 9 i T 1 i + a 1 2 T 4 i ϕ Ψ i d y 2 A ^ b 2 π α ¨ t 0 l T 13 i ϕ Ψ i d y + β t i 0 l 0.5177 b T 12 i 2 π B ^ V π   T 10 i A ^ V 2 π 2 ( T 5 i T 4 i T 10 i ) Ψ i 2 d y + β t i ˙ 0 l 0.5177 B ^ b 2 T 12 i 2 π 2 π T 11 i A ^ V b 2 π 2 T 4 i T 11 i Ψ i 2 d y + β t i ¨ 0 l A ^ b 2 π 2 T 3 i Ψ i 2 d y
u ¨ = 0.3414 a ^ u ˙ 0.0158 a ^ 2 u B ^ h ˙ t 0 l f 2 d y + B ^ V α t 0 l f ϕ d y + B ^ b ( 1 2 a ) α t ˙ 0 l f ϕ d y + β t i 0 l B ^ V π   T 10 i f Ψ i d y + β t i ˙ 0 l B ^ b 2 π T 11 i f Ψ i d y
v ¨ = 0.3414 a ^ v ˙ 0.0158 a ^ 2 v B ^ h ˙ t 0 l f ϕ d y + B ^ V α t 0 l ϕ 2 d y + B ^ b ( 1 2 a ) α t ˙ 0 l ϕ 2 d y + β t i 0 l B ^ V π   T 10 i ϕ Ψ i d y + β t i ˙ 0 l B ^ b 2 π T 11 i ϕ Ψ i d y
w i ¨ = 0.3414 a ^ w i ˙ 0.0158 a ^ 2 w i B ^ h ˙ t 0 l f Ψ i d y + B ^ V α t 0 l ϕ Ψ i d y + B ^ b ( 1 2 a ) α t ˙ 0 l ϕ Ψ i d y + β t i 0 l B ^ V π   T 10 i Ψ i 2 d y + β t i ˙ 0 l B ^ b 2 π T 11 i Ψ i 2 d y
A ^ = π ρ b 2
B ^ = 2 π ρ V b

2.3. Aeroelastic Equations of Motion

The governing equation for aeroelastic analysis can be expressed as:
M X ¨ + D X ˙ + K X = F a + F g
where M, D, and K represent the structural mass, structural damping (set to zero here), and structural stiffness matrices of the wing-flap system, respectively; X represents the displacement vector in generalized coordinates and the elements of these matrices are represented in Appendix A; F a represents the unsteady aerodynamic force vector; and F g gust load vector and it is included only for the gust response and alleviation analysis.

2.4. The Gust Model

To evaluate the gust load alleviation capability of SMTE, it is essential to determine the variations in root bending moment and shear force (shown in Figure 4) when the wing encounters gusts. It is noted that here only the gust velocity based on a discrete, 1-cosine gust profile is considered. The 1-cosine gust profile is defined, according to FAR Part 25, Section 25.341, as:
w g ( t ) = U d s 2 1 cos π V t H
where H, the gust gradient, is the distance parallel to the airplane’s flight path for the gust to reach its peak velocity, and it varies from 9.144 to 106.7 m. U d s is the design gust velocity, and it can be expressed as:
U d s = U r e f F g H 106 1 6
where U r e f , the reference gust velocity, has a magnitude of 17.07 m/s at sea level and reduces linearly from 17.07 to 13.4 m/s EAS at 15,000 feet. F g is the flight profile alleviation factor, and it is set to one. The airspeed is set to 30 m/s and the angle of attack is kept at zero.

2.5. Validation: Flutter Analysis

Due to the lack of available aeroelastic data on rectangular wings equipped with flap(s), the aeroelastic model developed here is validated using the Goland wing and the mechanical and geometric properties are listed in Table 2. For the sake of comparison, the stiffness of the three flaps is assumed to be very high (104 Nm/rad). This effectively limits the dynamics of the wing with SMTE to the bending-twisting of the cantilever wing. Figure 4 shows the variations in the first bending and the first torsion modes for the Goland wing with SMTE and the clean Goland wing (without SMTE). It can be seen that for the assumed high values of flap stiffness, the flutter mode for both scenarios is the 1st torsion mode. The flutter speed and frequency are 139.3 m/s and 68.87 rad/s for Goland with SMTE and 137.11 m/s and 69.9 rad/s for Goland without SMTE. This confirms the accuracy of the developed time-domain aeroelastic model. It should be noted that the flutter boundaries are estimated using the PK (frequency matching) method for validation purposes.

3. Parametric Study

The Goland wing, equipped with the SMTE concept which consists of three flaps, is used here for parametric aeroelastic analysis. The purpose of this study is to identify the influence of various design parameters on the aeroelastic boundary. The loads due to gust are set to zero in this section.

Aeroelastic Boundaries: Flutter and Divergence

The flutter boundary of the cantilever wing with SMTE is obtained by calculating the critical flutter speed, flutter frequency, and divergence speed for different flap stiffness values. Figure 5a,b shows the variation in flutter speed and frequency for different flap stiffness of a cantilever wing with a single flap. Similarly, Figure 5c,d shows flutter speed and frequency for varying flap stiffness of SMTE (where the flap stiffnesses were changed simultaneously by the same amount).
Figure 5 shows that the wing with a single flap and the SMTE wing with three-flap configurations are similar in the variation in flap stiffness. It can be seen that the flutter speed is very sensitive to flap stiffnesses for a range between 10 3 and 10 5   N m / r a d . Further, an increase in stiffness from 10 5   N m / r a d does not affect flutter speed. The flutter frequency also shows the same trend as in the case of speed as after 10 5   N m / r a d stiffness, it remains constant. It should be noted that regardless of stiffness, the flutter mode does not change and it remains torsion mode for both wings. For flap stiffnesses above 10 5 N m / r a d , the flutter speed for both configurations are increased from 90 m/s to 137 m/s and the frequency reduced from 85 rad/s to 69 rad/s. The flutter is determined to be due to the interaction of the first torsional and first bending mode, which excite each other and cause the torsional mode to become unstable. The low stiffness of the flap affects the globalized stiffness matrix and brings the first bending and first torsion frequencies closer to each other. This causes the coalescence of the modes to be brought about earlier in the presence of aerodynamic forces and causes the flutter to occur at a lower velocity.
Figure 6 shows the variation in flutter speed and frequency of a cantilever wing with SMTE (three-flap configuration) with the stiffness of each flap varied alone while keeping the other two stiffnesses at a value of baseline wing ( 10 8   N m / r a d ). It can be seen that from Figure 6, the flutter boundary is not very sensitive to flap stiffnesses, but for a range between 10 3 and 10 5   N m / r a d shows a slight change and then it remains constant even when the stiffness increses further and the variations are identical in both velocity and frequency. It should also be noted that for flap stiffnesses above 10 5   N m / r a d , the flutter speed and frequency for the three-flap configuration are very similar to those for the single-flap configuration for all three configurations.
Figure 7 shows the variation in divergence speed for a range of flap stiffness varied from 10 N m / r a d to 1010 N m / r a d . Four different configurations were considered; varying the stiffness of all three flaps of SMTE together ( K β ; blue), for varying inboard flap stiffness ( K β 1 ; red) alone whilst keeping the other two stiffness at 10 8   N m / r a d , for varying midboard flap stiffness ( K β 2 ; yellow) alone whilst keeping the other two stiffness at 10 8   N m / r a d , and for varying outboard flap stiffness ( K β 3 ; magenta) alone whilst keeping the other two stiffness at 10 8   N m / r a d . The divergence behavior of the wing can be attributed to the presence of higher torsional loads being subjected to the structure in case of low flap stiffness. For a wing at low angle of attack and significantly low flap stiffness, the trailing edge of the flap will tend to move upwards, creating a negative lift force. Since this force is applied aft of the shear center, it will create nose-up torsional moments at the wing leading edge, and add to the possibility of divergence instability. It can be seen that from Figure 7, the divergence boundary is different for all the four configurations for a range of flap stiffnesses ( 10   N m / r a d and 10 6   N m / r a d ) and above 10 6   N m / r a d it become a single line having a divergence speed of 252.28 m/s which is same as goland wing [27]. The first configuration, K β has the most deviated from the goland wing divergence speed and the inboard flap ( K β 1 ) is the least deviated from the goland wing divergence speed. It is evident from aeroelastic boundaries that the SMTE configuration is behaving very similarly to single flap configuration for stiffness values of above 106 N m / r a d , so in the following sections the analysis is carried out for stiffness value of 10 8   N m / r a d for all the flaps.

4. Response to Discrete Gusts

The response of the wing with SMTE (three-flap configuration) is assessed for upward and downward, discrete gust with the minimum (H = 9.07 m) and maximum (H = 106.7 m) gust gradients. The airspeed is set to 30 m/s, the angle of attack is fixed to 0.2 radian, and the flap stiffnesses are set to 10 8   N m / r a d . In this paper, gust velocity based on a discrete, 1-cosine gust profile is considered and it is assumed that the flaps are deflected before the gust hits the wing and the loads are used where the gust velocity is assumed to result in an instantaneous change in angle of attack. Only the maximum load point is taken into consideration for each gust gradient. The dynamic analysis is out of the scope of this paper because the flaps are deflected before the gust hits the wing. This is based on an assumption that a LIDAR or a sensor is located at the nose of the aircraft to detect gust and based on the signal from the LIDAR/sensor, the wing will be already prepared in optimum shape when a gust arrives at it. To evaluate the load alleviation capability of SMTE, the root shear force (RSF) and root bending moment (RBM) variations for different flap angles are presented using carpet plots (Figure 8, Figure 9, Figure 10 and Figure 11).
Upward-gust carpet plots are shown in Figure 8 and Figure 9 and downward-gust carpet plots are shown in Figure 10 and Figure 11. The carpet plots show an overview of the alleviation properties of SMTE for different combinations of flap deflections. For the upward gust, the upward flap deflections give alleviation in RSF and RBM while for the downward gust, the downward flap deflection gives alleviation in RSF and RBM. A more detailed gust load alleviation analysis and comparison with a single flap result are presented in the following sections.

4.1. Gust Load Alleviation Response Analysis

To identify the effectiveness of the location of different flaps in providing load alleviation capability of SMTE, each flap where deflected alone while keeping the other two flaps kept at zero flap angle in this subsection. Upward and downward gusts (1-cosine) are considered with an airspeed of 30 m/s, the angle of attack is fixed to 0.2 radian, the flap stiffnesses are set to 10 8   N m / r a d , and the gust gradient to H = 59 m. It should be noted that in this analysis, the flaps were deflected before the gust hits the wing and the loads are used where the gust velocity is assumed to result in an instantaneous change in the angle of attack.

4.1.1. Case 1: Inboard Flap Only

The inboard flap angle ( β 1 ) is deflected in a step of 10° angle up and down for upward and downward gusts, respectively. The resulting gust response for different gust gradients is shown in Figure 12 and Figure 13, respectively. It should be stressed that the deflections of the midboard and outboard flaps are set to zero.
Figure 12 (downward gust), when comparing with the baseline (blue) configuration, shows that RSF (8.14 kN) and RBM (33.24 kNm) are reduced by 33.63% and 11.32%, respectively, at a flap angle of 30° downward and H = 59 m. Increasing the flap deflection increases the alleviation percentage and at lower gust gradients a small deflection is enough to achieve more alleviation in RSF. From Figure 13 (upward gust), RSF (23.76 kN) and RBM (81.01 kNm) are reduced by 14.78% and 4.98%, respectively, for a flap angle of 30° upward and H = 59 m. Increasing the flap deflection increases the alleviation percentage.

4.1.2. Case 2: Midboard Flap Only

Similarly, the midboard flap angle ( β 2 ) is deflected in a step of 10° angle up and down for upward and downward gusts, respectively. The resulting gust response for different gust gradients for downward and upward gusts are shown in Figure 14 and Figure 15, respectively. It should be stressed that the deflections of the midboard and outboard flaps are set to zero.
Figure 14 (downward gust), when comparing with the baseline (blue) configuration, shows that RSF (8.10 kN) and RBM (24.77 kNm) are reduced by 33.91% and 33.92%, respectively, for a flap angle of 30° downward and H = 59 m. Increasing the flap deflection increases the alleviation percentage and a small deflection is enough to achieve more alleviation in both RSF and RBM. From Figure 15 (upward gust), RSF (23.73 kN) and RBM (72.55 kNm) are reduced by 14.90% and 14.91%, respectively, for a flap angle of 30° upward and H = 59 m. Increasing the flap deflection increases the alleviation percentage and both RSF and RBM vary gradually with flap deflection for all gust gradients.

4.1.3. Case 3: Outboard Flap Only

Only the outboard flap angle ( β 3 ) is deflected in a step of 10° angle up and down for upward and downward gusts, respectively. The resulting gust responses for different gust gradients for downward and upward gusts are shown in Figure 16 and Figure 17, respectively.
Figure 16 (downward gust), when compared with baseline (blue), shows that RSF (8.08 kN) and RBM (16.35 kNm) are reduced by 34.07% and 56.37%, respectively, for a flap angle of 30° downward and H = 59 m. Increasing the flap deflection increases the alleviation percentage and a small deflection is enough to achieve more alleviation in both RSF and RBM. From Figure 17 (upward gust), RSF (23.71 kN) and RBM (64.13 kNm) are reduced by 14.97% and 24.91%, respectively, for a flap angle of 30° upward and H = 59 m. Increasing the flap deflection increases the alleviation percentage and both RSF and RBM vary gradually with flap deflection for all gust gradients.
From all the above cases, for downward gust, the RBM alleviation percentages are increasing from inboard to outboard but the RSF alleviation percentages are not showing any significant variation. For upward gust, both RSF and RBM alleviation percentages increase from inboard to outboard. It should be noted that the alleviation percentage is more in a downward gust as compared to an upward gust and this is because, for the downward gust, the downgoing flap will produce more aerodynamic opposing force than the upgoing flap as in the case of an upward gust. A more detailed summary of the response analysis is tabulated in Table 3, Table 4, Table 5 and Table 6 for single-flap and multi-flap operations of SMTE.

4.2. Optimal Flap Deflections for Load Alleviation

In this work, optimal flap deflections are obtained for different load alleviation objectives and constraints. The inboard, midboard, and outboard flap angles ( β 1 ,   β 2   a n d   β 3 ) are deflected in a step of 10° angle up and down for a downward gust. The gust (1-cosine) gradient is set at H = 9.07 m, while the airspeed is fixed to 30 m/s, the angle of attack is set to 0° and the flap stiffnesses are set to 10 8   N m / r a d . A set of 343 flap configurations are considered. For each configuration, the RSF, the RBM, and the percentage alleviation are computed. The results of the two different optimizations are reported as:
  • Objective 1: Minimize RBM subject to a 30% reduction in RSF: The objective is to find the optimal configuration that minimizes RBM and reduces the RSF by 30% relative to the baseline configuration (10,198 N) (flaps angles are zero).
  • Objective 2: Minimize RSF subject to a 30% reduction in RBM: The objective is to find the optimal configuration that minimizes RSF and reduces the RBM by 30% relative to the baseline configuration (31,084 N) (flaps angles are zero).
Table 7 presents the optimal configuration for each objective. In the case of Objective 1, the inboard and midboard flaps deflect to 30° while the outboard flap deflects to −28° giving a 66% reduction in RBM whilst reducing the RSF by 30%. However, for Objective 2, the midboard and outboard flap deflects to 30° while the inboard deflects to −4° resulting in a 52% reduction in RSF and a 30% reduction in RBM. To minimize RBM, the inboard flap deflects to the maximum angle while the midboard and outboard flaps deflect downward to the maximum angle. On the other hand, to minimize RSF, the inboard and midboard flaps deflect downward to the maximum angle whereas the outboard flap deflects slightly upward. It is worth noting that the ability to vary the spanwise camber distribution gives wider control on the RBM and RSF which is not the case with a single flap running from root to tip.

4.3. Comparison between a Single Flap and an SMTE (with Three Flaps)

The gust response of a single flap is obtained and presented in Table 8. The comparison is carried out in such a way that a single flap angle at 20° downward deflection can provide a 56% reduction in RSF and a 56% reduction in RBM. However, a 3-flap SMTE is capable of providing a 56% reduction in RBM with a 63% reduction in RSF. In addition, the SMTE is capable of providing a 56% reduction in RSF with a 78% reduction in RBM. Since SMTE allows spanwise camber variation, it gives the aircraft designer a wider range of options for load alleviation when compared to a single flap concept.

5. PID Controller for the SMTE Wing

Due to its simplicity and robustness, a PID controller is designed to meet the required loads at the wing root as shown in Figure 18. The PID controller delivers actuation force on SMTE to move from one position to another while maintaining the required load alleviation. The design of the PID controller did not account for actuator dynamics. To account for actuator dynamics, the type of actuation system must be determined and this requires trade-off studies that are beyond the scope of this paper which mainly aims to conduct a parametric aeroelastic study. The closed loop ensures that the achieved tip deflection is very close to the desired one. The controller parameters are proportional gain Kp, derivative gain Kd, and integral gain Ki. Percentage load alleviation was used to estimate control gains of the control system. There are many ways to tune a PID controller and there are also special methods for direct tuning based on simple process experiments. The PID tuning has been carried out based on the gust load alleviation values from the optimization study. A trial-and-error strategy was used to obtain the optimal control gain parameters.
This parametric study aims to have an initial investigation of a controller for SMTE. The gust alleviation property of SMTE was used to design controller parameters. If SMTE is used as a gust load alleviation device, it must be fast and must settle to the targeted wing root loads in the shortest period and with minimum overshoot. The selection of PID control gains is performed using the gust load alleviation values from the optimization study in Section 4.2. A trial-and-error strategy was used to reach control gain parameters. For a flight condition of 30 m/s, an angle of attack of zero radians, and at sea level, with a gust gradient of H = 9.07 m SMTE with three-flap configurations the controller parameters have been found out. Controllers with three different proportional gains, integral and differential gains have arrived at different combinations of flap angles. The simulations of the SMTE with three PID controllers were performed for different combinations of flap angles. For a given actuation time the integral gain is independent and a −1 (constant value) for differential gain gives minimum overshoot and stabilizes the system. Controller parameters for 30 percent load alleviation in root bending moment for some combinations of flap angles are tabulated in Table 9.

6. Conclusions

This paper presented an investigation of the dynamic aeroelastic characteristics of a cantilever wing equipped with a spanwise morphing trailing edge (SMTE) concept. The structure of the wing was represented using the Euler–Bernoulli beam, the Rayleigh–Ritz method was used to derive the generalized equation of motion, and Theodorsen’s unsteady aerodynamic theory was used to estimate the aerodynamic loads. The governing aeroelastic equations of the wing with a three-flap configuration were derived in the time domain. A flutter analysis was conducted to compute the flutter boundaries. The flutter boundaries of the wing with a three-flap configuration are compared with the wing with a single flap configuration. Furthermore, several scenarios were studied to assess the feasibility of SMTE as a load alleviation device when subject to discrete gusts. Finally, the implementation and validation of a controller for gust load alleviation were studied and controller parameters are tuned for a specific gust model. The simulation results showed the potential load alleviation capabilities of the SMTE when compared to a single flap configuration.

Author Contributions

Conceptualization, J.S.P. and R.M.A.; methodology, J.S.P. and R.M.A.; validation, J.S.P. and R.M.A.; formal analysis, J.S.P.; investigation, J.S.P. and R.M.A.; resources, R.M.A.; data curation, J.S.P.; writing—original draft preparation, J.S.P.; writing—review and editing, J.S.P., R.M.A., Z.H. and M.A.; supervision, R.M.A.; project administration, J.S.P. and R.M.A.; funding acquisition, R.M.A. All authors have read and agreed to the published version of the manuscript.

Funding

The work presented herein has been funded by the Abu Dhabi Education Council Award for Research Excellence Program (AARE 2019) through grant number AARE19-213.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

anon-dimensional distance from airfoil mid-chord to an elastic axis
bairfoil semi-chord
C(s) Theodorsen transfer function in the Laplace domain
c i non-dimensional distance from the airfoil mid-chord to the ith flap hinge line
LLift
M α moment of wing-flap about an elastic axis
M β i moment about the hinge axis of the ith flap
mmass of wing-flap (per unit span)
yspanwise location measured relative to the wing root
h plunge displacement at the elastic axis
α pitch angle
β i flap angle of the ith flap
V true airspeed
f ( y ) bending shape function
ϕ ( y ) torsion shape function
Ψ ( y ) flap shape function
Ttotal kinetic energy
Utotal potential energy
ρ air density
S Laplace variable
Acronyms
Dofdegree of freedom
SMTEspanwise morphing trailing edge
GLA gust load alleviation
BWBblended wing boy
PIDproportional integral derivative

Appendix A

It should be noted that the aeroelastic equation is of the general form M X ¨ + D X ˙ + K X , where X is a vector of the system variables with an overdot representing time derivatives, the structural mass (M), and structural stiffness (K) matrices can be expressed as:
M = m S α S β i S α I α I β i + b c i a S β i S β i I β i + b c i a S β i I β i
K = K h 0 0 0 K α 0 0 0 K β i
And
X = h α β i

References

  1. Khalil, A.; Fezans, N. Performance Enhancement of Gust Load Alleviation Systems for Flexible Aircraft using H∞Optimal Control with Preview. In Proceedings of the AIAA Scitech 2019 Forum, San Diego, CA, USA, 7–11 January 2019. [Google Scholar] [CrossRef]
  2. Bernhammer, L.O.; Navalkar, S.T.; Sodja, J.; De Breuker, R.; Karpel, M. Experimental investigation of an autonomous flap for load alleviation. In Proceedings of the 56th AIAA/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference, Kissimmee, FL, USA, 5–9 January 2015. [Google Scholar] [CrossRef]
  3. Wildschek, A.; Stroscher, F.; Klimmek, T.; Šika, Z.; Vampola, T.; Valášek, M.; Gangsaas, D.; Aversa, N.; Berard, A.; Raumfahrt, L.U. Gust load alleviation on a large blended wing body airliner. In Proceedings of the 27th Congress of the International Council of the Aeronautical Sciences 2010, ICAS 2010, Nice, France, 19–24 September 2010. [Google Scholar]
  4. Huvelin, F.; Dequand, S.; Lepage, A.; Liauzun, C. On the validation and use of high-fidelity numerical simulations for gust response analysis. AerospaceLab J. 2018. [Google Scholar] [CrossRef]
  5. Bekemeyer, P.; Thormann, R.; Timme, S. Investigation into gust load alleviation using computational fluid dynamics. In Proceedings of the International Forum on Aeroelasticity and Structural Dynamics 2019, IFASD 2019, Savannah, GA, USA, 10–13 June 2019. [Google Scholar]
  6. Pusch, M. Allocation of distributed flaps for gust load alleviation. In Proceedings of the 1st Annual IEEE Conference on Control Technology and Applications, CCTA 2017, Kohala Coast, HI, USA, 27–30 August 2017; pp. 2120–2125. [Google Scholar] [CrossRef]
  7. Dias, J.N.; Hubbard, J.J.E. Active Spanwise Lift Control: A Distributed Parameter Approach; American Institute of Aeronautics and Astronautics, Inc.: Reston, VA, USA, 2020. [Google Scholar] [CrossRef]
  8. Pusch, M.; Knoblach, A.; Kier, T. Integrated optimization of control surface layout for gust load alleviation. CEAS Aeronaut. J. 2019, 10, 1059–1069. [Google Scholar] [CrossRef]
  9. Barbarino, S.; Bilgen, O.; Ajaj, R.M.; Friswell, M.I.; Inman, D.J. A Review of Morphing Aircraft. J. Intell. Mater. Syst. Struct. 2011, 22, 823–877. [Google Scholar] [CrossRef]
  10. Ajaj, R.M.; Parancheerivilakkathil, M.S.; Amoozgar, M.; Friswell, M.I.; Cantwell, W.J. Recent developments in the aeroelasticity of morphing aircraft. Prog. Aerosp. Sci. 2021, 120, 100682. [Google Scholar] [CrossRef]
  11. Pankonien, A.M.; Inman, D.J. Spanwise morphing trailing edge on a finite wing. In Active and Passive Smart Structures and Integrated Systems; SPIE: San Diego, CA, USA, 2015; p. 94310T. [Google Scholar] [CrossRef]
  12. Pankonien, A.M.; Inman, D.J. Aerodynamic performance of a spanwise morphing trailing edge concept. In Proceedings of the ICAST 2014—25th International Conference on Adaptive Structures and Technologies, The Hague, The Netherlands, 6–8 October 2014. [Google Scholar]
  13. Haghighat, S.; Martins, J.R.R.A.; Liu, H.H.T. Aeroservoelastic Design Optimization of a Flexible Wing. J. Aircr. 2012, 49, 432–443. [Google Scholar] [CrossRef]
  14. Capello, E.; Guglieri, G.; Quagliotti, F. A Comprehensive Robust Adaptive Controller for Gust Load Alleviation. Sci. World J. 2014, 2014, 1–12. [Google Scholar] [CrossRef] [PubMed]
  15. Magar, K.T.; Fuchi, K.; Pankonien, A.M.; Beblo, R.V.; Reich, G.W. Active Control of Origami Inspired Camber Morphing Airfoil for Gust Load Alleviation. In Proceedings of the ASME 2018 Conference on Smart Materials, Adaptive Structures and Intelligent Systems, SMASIS, San Antonio, TX, USA, 10–12 September 2018. [Google Scholar] [CrossRef]
  16. Bourchak, M.; Ajaj, R.M.; Flores, E.I.S.; Khalid, M.; Juhany, K.A. Optimum design of a PID controller for the adaptive torsion wing. Aeronaut. J. 2015, 119, 871–889. [Google Scholar] [CrossRef]
  17. Kuznetsov, N.V.; Andrievsky, B.; Zaitceva, I.; Kudryashova, E.V.; Kuznetsova, O.A. Adaptive Suppression of Wing Flutter Under Actuator Saturation and Time Sampling. IFAC-PapersOnLine 2022, 55, 689–694. [Google Scholar] [CrossRef]
  18. Mozaffari-Jovin, S.; Firouz-Abadi, R.; Roshanian, J. Flutter of wings involving a locally distributed flexible control surface. J. Sound Vib. 2015, 357, 377–408. [Google Scholar] [CrossRef]
  19. Vindigni, C.R.; Orlando, C. Simple adaptive V-stack piezoelectric based airfoil flutter suppression system. J. Vib. Control. 2022, 29, 2802–2816. [Google Scholar] [CrossRef]
  20. Muñoz, Á.; García-Fogeda, P. Active Flutter Suppression of a Wing Section in a Compressible Flow. Aerospace 2022, 9, 804. [Google Scholar] [CrossRef]
  21. Wu, Z.; Lu, J.; Rajput, J.; Shi, J.; Ma, W. Adaptive Neural Control Based on High Order Integral Chained Differentiator for Morphing Aircraft. Math. Probl. Eng. 2015, 2015, 1–12. [Google Scholar] [CrossRef]
  22. Seigler, T.M.; Neal, D.A.; Bae, J.-S.; Inman, D.J. Modeling and Flight Control of Large-Scale Morphing Aircraft. J. Aircr. 2007, 44, 1077–1087. [Google Scholar] [CrossRef]
  23. Irani, S.; Sarrafzadeh, H.; Amoozgar, M.R. Bifurcation in a 3-DOF Airfoil with Cubic Structural Nonlinearity. Chin. J. Aeronaut. 2011, 24, 265–278. [Google Scholar] [CrossRef]
  24. Chauhan, S.S.; Martins, J.R.R.A. Low-Fidelity Aerostructural Optimization of Aircraft Wings with a Simplified Wingbox Model Using OpenAeroStruct. In Proceedings of the EngOpt 2018 Proceedings of the 6th International Conference on Engineering Optimization, Libson, Portugal, 14 September 2018. [Google Scholar] [CrossRef]
  25. Theodorsen, T. General Theory of Aerodynamic Instability and the Mechanism of Flutter; NACA Technical Report no.496; NASA: Washington, DC, USA, 1 January 1949; Volume 496. [Google Scholar]
  26. Brunton, S.; Rowley, C. Low-Dimensional State-Space Representations for Classical Unsteady Aerodynamic Models. In Proceedings of the 49th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition, Orlando, FL, USA, 4–7 January 2011. [Google Scholar] [CrossRef]
  27. Ajaj, R.M.; I Friswell, M. Aeroelasticity of compliant span morphing wings. Smart Mater. Struct. 2018, 27, 105052. [Google Scholar] [CrossRef]
Figure 1. A finite wing with the SMTE concept [12].
Figure 1. A finite wing with the SMTE concept [12].
Aerospace 10 00809 g001
Figure 2. Schematic of an airfoil section with a control surface.
Figure 2. Schematic of an airfoil section with a control surface.
Aerospace 10 00809 g002
Figure 3. A top view of the SMTE.
Figure 3. A top view of the SMTE.
Aerospace 10 00809 g003
Figure 4. Flutter modes of Goland wing; blue for torsion and red for bending (with SMTE, K β set to very high); green for torsion and black for bending (without SMTE).
Figure 4. Flutter modes of Goland wing; blue for torsion and red for bending (with SMTE, K β set to very high); green for torsion and black for bending (without SMTE).
Aerospace 10 00809 g004
Figure 5. Variation in flutter speed and frequency with flap stiffness (a,b) for a single flap, and (c,d) for the SMTE wing.
Figure 5. Variation in flutter speed and frequency with flap stiffness (a,b) for a single flap, and (c,d) for the SMTE wing.
Aerospace 10 00809 g005
Figure 6. Flutter speed and frequency of SMTE wing; (a,b) for varying inboard flap stiffness ( K β 1 ) alone; (c,d) for varying midboard flap stiffness ( K β 2 ) alone; (e,f) for varying outboard flap stiffness ( K β 3 ) alone.
Figure 6. Flutter speed and frequency of SMTE wing; (a,b) for varying inboard flap stiffness ( K β 1 ) alone; (c,d) for varying midboard flap stiffness ( K β 2 ) alone; (e,f) for varying outboard flap stiffness ( K β 3 ) alone.
Aerospace 10 00809 g006
Figure 7. Divergence speed vs. stiffness of SMTE wing; for varying the stiffness of all flaps together ( K β ); for varying inboard flap stiffness ( K β 1 ) alone; for varying midboard flap stiffness ( K β 2 ) alone; for varying outboard flap stiffness ( K β 3 ) alone.
Figure 7. Divergence speed vs. stiffness of SMTE wing; for varying the stiffness of all flaps together ( K β ); for varying inboard flap stiffness ( K β 1 ) alone; for varying midboard flap stiffness ( K β 2 ) alone; for varying outboard flap stiffness ( K β 3 ) alone.
Aerospace 10 00809 g007
Figure 8. Carpet plot for upward gust for H = 9.07 m.
Figure 8. Carpet plot for upward gust for H = 9.07 m.
Aerospace 10 00809 g008
Figure 9. Carpet plot for upward gust for H = 106.7 m.
Figure 9. Carpet plot for upward gust for H = 106.7 m.
Aerospace 10 00809 g009
Figure 10. Carpet plot for downward gust for H = 9.07 m.
Figure 10. Carpet plot for downward gust for H = 9.07 m.
Aerospace 10 00809 g010
Figure 11. Carpet plot for downward gust for H = 106.7 m.
Figure 11. Carpet plot for downward gust for H = 106.7 m.
Aerospace 10 00809 g011
Figure 12. RSF and RBM at different gust gradients for downward gust.
Figure 12. RSF and RBM at different gust gradients for downward gust.
Aerospace 10 00809 g012
Figure 13. RSF and RBM at different gust gradients for upward gust.
Figure 13. RSF and RBM at different gust gradients for upward gust.
Aerospace 10 00809 g013
Figure 14. RSF and RBM at different gust gradients for downward gust.
Figure 14. RSF and RBM at different gust gradients for downward gust.
Aerospace 10 00809 g014
Figure 15. RSF and RBM at different gust gradients for upward gust.
Figure 15. RSF and RBM at different gust gradients for upward gust.
Aerospace 10 00809 g015
Figure 16. RSF and RBM at different gust gradients for downward gust.
Figure 16. RSF and RBM at different gust gradients for downward gust.
Aerospace 10 00809 g016
Figure 17. RSF and RBM at different gust gradients for upward gust.
Figure 17. RSF and RBM at different gust gradients for upward gust.
Aerospace 10 00809 g017
Figure 18. Block diagram of the PID controller.
Figure 18. Block diagram of the PID controller.
Aerospace 10 00809 g018
Table 1. Geometric and Mechanical Properties of the Multi-flap Wing.
Table 1. Geometric and Mechanical Properties of the Multi-flap Wing.
ParameterFlap 1Flap 2Flap 3
Torsional stiffness K β 1 K β 2 K β 3
Length l 1 l 2 l 3
Mass per unit span m f 1 m f 2 m f 3
Mass moment of inertia per unit span I β 1 I β 2 I β 3
Static moments S β 1 S β 2 S β 3
Distance of flap mass center from c x β 1 x β 2 x β 3
Spanwise position 0 y l 1 l 1 y l 2 + l 1 l 1 + l 2 < y l 1 + l 2 + l 3
Table 2. Properties of Goland Wing [27].
Table 2. Properties of Goland Wing [27].
SpecificationsGoland Wing
Half span (m)6.096
Chord (m)1.8288
Mass per unit length (kg/m)35.71
Moment of inertia (50% chord) (kgm)8.64
Spanwise elastic axis (from LE)33%
Center of gravity (from LE)43%
Spanwise bending rigidity (Nm2)9.77 × 106
Torsional rigidity (Nm2)0.987 × 106
Density of air kg/m31.225
Table 3. Summary of the percentage alleviation response of downward gust for different gust gradients.
Table 3. Summary of the percentage alleviation response of downward gust for different gust gradients.
Downward GustH = 9 mH = 59 mH = 106.7 m
% Alleviation of RSF% Alleviation of RBM% Alleviation of RSF% Alleviation of RBM% Alleviation of RSF% Alleviation of RBM
Inboard Flap only ( β 1 = 30 ° )602034112910
Midboard Flap only ( β 2 = 30 ° )606034342929
Outboard Flap only ( β 3 = 30 ° )609934562948
Table 4. Summary of the percentage alleviation response of upward gust for different gust gradients.
Table 4. Summary of the percentage alleviation response of upward gust for different gust gradients.
Upward GustH = 9 mH = 59 mH = 106.7 m
% Alleviation of RSF% Alleviation of RBM% Alleviation of RSF % Alleviation of RBM% Alleviation of RSF % Alleviation of RBM
Inboard Flap only ( β 1 = 30 ° )186154.9144.6
Midboard Flap only ( β 2 = 30 ° )181815151414
Outboard Flap only ( β 3 = 30 ° )183115251423
Table 5. Summary of the percentage alleviation response of upward gust for a combination of flap deflections.
Table 5. Summary of the percentage alleviation response of upward gust for a combination of flap deflections.
Upward GustH = 9 mH = 59 mH = 106.7 m
% Alleviation of RSF % Alleviation of RBM% Alleviation of RSF % Alleviation of RBM% Alleviation of RSF % Alleviation of RBM
Inboard and Midboard Flap ( β 1 = 30 ° )372530202818
Inboard and Outboard Flap
( β 2 = 30 ° )
373730302828
Midboard and Outboard Flap
( β 3 = 30 ° )
374915401437
Table 6. Comparison of percentage alleviation response for downward and upward gust.
Table 6. Comparison of percentage alleviation response for downward and upward gust.
Downward Gust (H = 59 m)Upward Gust (H = 59 m)
% Alleviation of RSF% Alleviation of RBM% Alleviation of RSF% Alleviation of RBM
Inboard Flap only
( β 1 = 30 ° )
3411155
Midboard Flap only ( β 2 = 30 ° )34341515
Outboard Flap only ( β 3 = 30 ° )34561525
Table 7. Parameters of the selected configuration.
Table 7. Parameters of the selected configuration.
Flap Angle (in Degrees)RSF (in N)RBM (in Nm)% Alleviation of RSF% Alleviation of RBM
β 1 β 2 β 3
Objective 1−28303010,19815,1493066
Objective 23030−4704731,0845230
Table 8. Comparison between a single flap and an SMTE (with three flaps).
Table 8. Comparison between a single flap and an SMTE (with three flaps).
Flap Angle (in Degrees)RSF
(in N)
RBM
(in Nm)
% Alleviation of RSF% Alleviation of RBM
Single Flap β
20643919,6985656
SMTE β 1 β 2 β 3
For the same RSF10303043998,3445678
For the same RBM30201834319,6986356
Table 9. Controller parameters.
Table 9. Controller parameters.
Flap Angle (in Degrees)Proportional Gain, Kp1Proportional Gain, Kp2Proportional Gain, Kp3Differential Gain, Kd
β 1 β 2 β 3
10018−5.890−10.54−1
101012−5.67−5.67−6.72−1
10206−5.45−10.90−3.18−1
10300−5.25−15.760−1
10016−11.410−9.05−1
101010−10.98−5.49−5.41−1
10202−15.40−10.27−0.94−1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Pilakkadan, J.S.; Ajaj, R.M.; Haider, Z.; Amoozgar, M. On the Aeroelasticity of a Cantilever Wing Equipped with the Spanwise Morphing Trailing Edge Concept. Aerospace 2023, 10, 809. https://doi.org/10.3390/aerospace10090809

AMA Style

Pilakkadan JS, Ajaj RM, Haider Z, Amoozgar M. On the Aeroelasticity of a Cantilever Wing Equipped with the Spanwise Morphing Trailing Edge Concept. Aerospace. 2023; 10(9):809. https://doi.org/10.3390/aerospace10090809

Chicago/Turabian Style

Pilakkadan, Jafar S., Rafic M. Ajaj, Zawar Haider, and Mohammadreza Amoozgar. 2023. "On the Aeroelasticity of a Cantilever Wing Equipped with the Spanwise Morphing Trailing Edge Concept" Aerospace 10, no. 9: 809. https://doi.org/10.3390/aerospace10090809

APA Style

Pilakkadan, J. S., Ajaj, R. M., Haider, Z., & Amoozgar, M. (2023). On the Aeroelasticity of a Cantilever Wing Equipped with the Spanwise Morphing Trailing Edge Concept. Aerospace, 10(9), 809. https://doi.org/10.3390/aerospace10090809

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop