Recent Progress on Aeroelasticity of High-Performance Morphing UAVs

Binbin Lv,Jun Zha,Kaichun Zeng,Hongtao Guo,Li Yu and Peng Zhang

Institute of High Speed Aerodynamics,China Aerodynamics Research and Development Center,Mianyang,621000,China

ABSTRACT The high-performance morphing aircrafthas become a research focus all over the world.The morphing aircraft,unlike regular unmanned aerial vehicles (UAVs), has more complicated aerodynamic characteristics, making itmore difficultto conduct its design,model analysis,and experimentation.This paper reviews the recent process and the current status of aeroelastic issues, numerical simulations, and wind tunnel test of morphing aircrafts.The evaluation of aerodynamic characteristics,mechanism,and relevant unsteady dynamic aerodynamic modeling throughout the morphing process are the primary technological bottlenecks formorphing aircrafts.The unstable aerodynamic forces have a significant impact on the aircrafthandling characteristics, control law design, and flight safety.In the past, the structural analysis of morphing aircrafts, flight dynamics modeling, computational mesh morphing technology,and aerodynamic calculation were performed in promoting the development of next generation UAVs,with nonlinear dynamic challenges includingtransonic aeroelastic problems and high angle of attack aeroelastic problems.At present,many facets of these difficulties,together with the accompanying numerical simulation studies,remain under-explored.In addition,wind tunnel experiments face significant challenges in the dynamic morphing process.Finally,dynamic unsteady aerodynamic characteristics in the continuous morphing process still need to be verified by more related experiments.

KEYWORDS Morphing aircraft;aeroelasticity;numerical simulation;wind tunnel test

1 Introduction

1.1 Morphing Aircraft Technology

Traditional fixed-wing aircrafts were designed first for the primary flight environment and tasks(such as cruise,reconnaissance,and combat,etc.), while also accounting for critical flight processes like take-off and landing, and then compromised optimized design [1–3].As a result, aircrafts that perform different functions and travel at different speeds would have extremely distinct designs,such as subsonic long-range drones represented by the Global Hawk,supersonic fighters represented by the F-16,and hypersonic vehicles represented by the HTV-2 and X-51A.Flying creatures provide excellent illumination.In cruising,hovering,diving,and attacking,they will adopt a proper shape to provide the best aerodynamic force and flight performance.In view of this,scientists and engineers recommended designing a new concept aircraft called the morphing aircraft, which can accomplish performance optimization at all phases by changing the layout shape while completing different missions under varied flight conditions[4–6].

At present, research on morphing aircrafts is focused on the following topics: (1) what type of morphing the aircraft needs to generate;(2)how to realize morphing;(3)the safety during dynamic morphing; (4) how to control; (5) the cost.Aerodynamics, flight control, materials, structure, drive,morphing control,and other basic ideas have been explored[7,8].The term“aircraft morphing”refers to a wide range of spatial scale and time scale(local,dispersed,and overall).It is divided into three categories based on the morphing scale and realized function:(1)local morphing(small morphing);(2) distributed morphing (medium-scale morphing); (3) integral morphing (large-scale morphing).Meanwhile, it can also be divided into two types based on the realization method: (1) mechanical morphing;(2)smart materialsand structures based morphing[9].

The research on morphing aircrafts began in 1903 and developed from the Wright brothers’initial aircraft,which depended on flexible deformable wings,to the variable-sweep wing technology exemplified by the F-14 in the 1960s and 1970s.In recent years,the continuous variable trailing edge camber technology used on the Gulfstream III aircraft has been a hot topic in the aerospace industry.It may triggertechnological advancements in aircrafts and have a disruptive influence.It has attracted wide attention all over the world[10].From the perspective of basic scientific problems and bottleneck technologies,the main problems include two aspects:(1)Aerodynamics,flight mechanics,and flight control of morphing aircrafts;(2)morphing structure,drive,and morphing control[6,11].

Since the 1980s,a lot of manpower and financial resources have been invested to develop morphing technologies to improve the overall performance of aircrafts[12].Some research projects onthe aircraft intelligent morphing control technology were carried out,such as the Mission-Adaptive Wing(MAW),the Active Flexible Wing (AFW), the Smart Wing Project (SWP), the Active Compliant Trailing Edge flaps(ACTE),and the Variable Continuously Camber Trailing-Edge Flaps(VCCTEF)[13–15].Later, some mature wing morphing technologies in Europe and the United States have passed wind tunnel tests and flight tests, and are getting increasingly close to practical engineering applications.However, the research work on morphing aircrafts still needs new conceptual design and is still in the exploratory research stage [16].Nevertheless, some scholars have researched the wing morphing technology.For example, Huang [17] explored the application of a new type of zero Poisson’s ratio honeycomb structure in variable camber wings.Li [18] designed a multi-joint adaptive wing trailing edge based on memory alloy differential drive.Liu et al.[19]proposed a design scheme for the front and rear edges of the variable camber wing based on a distributed flexible structure, a closed-chain skeletal mechanism, and a rigid-flexible coupling mechanism.Although specific technologies have been studied in the existing literature,there isno systematic and in-depth review on key aerodynamics and related technical issues.It is highly complex and difficult due to the lack of strong support from basic technologies such as smart materials,flexible skins,and memory alloys,as well as variant technologies involving aerodynamics,structure,control,aeroelasticity,and interdisciplinarity[20,21].

As one filiation of mechanics, aeroelasticity mainly focuses on mechanical characteristics, the law of gas flow,and concomitant physical and chemical changes of aircrafts flying in the air.Along with the development of aircraft industry and jet propulsion technology, it gradually becomes a subject.Usually,the main research contents of aeroelasticity are velocity,temperature,pressure and density of gas,lift force and drag force of aircrafts and change laws,etc.when aircrafts are flying in various conditions.Based on different flight environments,morphing UAVs can autonomously change their aerodynamic configuration aiming to efficiently finish corresponding flight mission.Therefore,aeroelasticity is the key issue formorphing UAVs.It mainly involves dynamic problems, aeroelastic numerical simulation,aeroelastic wind tunnel test,etc.,which will be reviewed in the following texts.

1.2 Numerical Simulation of Aeroelasticity

In the aerospace field,the process of aeroelastic analysis methods is similar to that of many other disciplines.Before the application of electronic computers, it mainly focuses on experimental and analytical methods[22,23].However,there are many difficulties in aeroelastic problems,especially in the high-speed wind tunnel aeroelastic experiment.Therefore,the calculation method of aeroelasticity has become an important research and development direction[24,25].

The aeroelastic problems of early aircrafts mainly lie in the wings, while the early wings were straight with relatively large chords.Therefore,the aeroelastic analysis methods can only be developed around the aerodynamic characteristics of two-dimensional(2D)wings[13].After the 1950s,with the development of high-speed aircrafts, the aeroelastic analysis method of high-speed wings based on piston theory was developed.However,the calculation method that can solve the aeroelastic problem of actual aircrafts depends on the development of computer technology[26].In the 1970s,the unsteady lift surface method based on the potential streamline theory was widely used and developed rapidly.It is still the main method in engineering nowadays and has been embedded in the famous analysis software NASTRAN.In the late 1980s, as the computer processing speed increased, aeroelastic analysis methods based on the full-speed potential equation and Euler equations were developed.In particular, Euler equations helped solve the long-term trouble (transonic aeroelastic problem), but they are under the assumption of inviscid flow, thus failing to consider the significant influence of viscous effects.In recent years, aeroelastic analysis methods based on the N-S equation have been developed.However, these methods entail huge amount of numerical calculation of the unsteady flow.The computer technology is still incapable of realizing the aeroelastic numerical simulation of the viscous flow problem [27].The development direction of aeroelasticity is thus clear.With the further increase in computer speed,the development of parallel computer speed,and the advancement of numerical calculation methods, numerical simulation to solve complex aeroelastic problems will gradually become a reality[28,29].

1.3 Aeroelastic Wind Tunnel Test

Aeroelastic wind tunnel test has the advantages of high reliability (compared to numerical calculation)and low cost(compared to flight test).The scaled-down similar model is tested in a wind tunnel,which can simulate the aeroelastic characteristics of the aircraft more accurately.Therefore,the wind tunnel test has always been one of the important means to analyze and evaluate the aeroelastic characteristics of aircrafts.In the design stage, the aeroelastic characteristics of the aircraft must be analyzed through wind tunnel tests to verify whether its strength designed and performance indicators meet the requirements,and corresponding optimization and improvement should be made.The wind tunnel flutter test is usually conducted before the flight test,so that its results can be used to provide a basis for and reduce the risk of the flight test.In addition, the aeroelastic wind tunnel test also provides the technical support for the verification of advanced aerodynamics, structural coupling design,and analysis methods.The transonic dynamics tunnel(TDT)of the Langley Center of NASA was completed in the early 1960s.By 2016,more than 600 wind tunnel tests had been carried out.Most of them are aeroelastic tests,including fixed-wing aircrafts,rotary-wing aircrafts,launch vehicles,and other configurations,which have greatly facilitated the aeroelastic design and analysis of the research models,the engineering application of new layouts and new technologies,and numerical calculation methods.Therefore,the aeroelastic wind tunnel test technology is a solid foundation for the research and development of aerospace vehicles[30,31].

This paper takes high-performance morphing aircrafts as the research object,aiming at reviewing the aeroelastic problems in the morphing process, comprehensively explaining the key problems in morphing aircrafts and the research progress, and respectively analyzing the aeroelastic numerical simulation methods.Finally,aeroelastic wind tunnel test technologies are elaborated.

2 Aeroelastic Problems of High-Performance UAVs

2.1 Development History of Aeroelasticity

The aeroelastic problems of aerospace vehicles have accompanied the development of aircrafts.As early as 1903,Langley’s monoplane suffered a broken wing during its first powered flight.After the Wright brothers succeeded in achieving a powered flight that was heavier than air,the development of aircrafts accelerated rapidly,and the aeroelastic problem became increasingly prominent.In the 1920s,the research on wing flutter analysis had achieved rapid progress[32].Representative research included Theodorsen’s unsteady aerodynamic analysis method of the simple harmonic vibration airfoil, and studies on the theory of the wing flutter and other content [33].After entering the 1950s, due to the emergence of high-speed aircrafts, the problem of thermo-aeroelasticity attracted the attention of researchers.Because the hypersonic aircraft was mainly focused on spacecraft at that time, the problem of thermo-aeroelasticity was not systematically solved[34,35].After the 1980s,the aircrafts focused on the transonic flight,and the study of nonlinear aeroelastic problems related to transonic speeds obtained systematic and in-depth investigations[36,37].In recent years,with the development of micro-aircrafts, the aeroelasticity of flexible deformable wings has also attracted wide scholarly attention[38,39].In short,the research of aerospace aeroelasticity is closely related to the development of aircrafts and mechanical analysis methods.Tables 1 and 2 respectively summarize the development history of aeroelasticity and unsteady aerodynamic calculation methods.

Table 1: The development history of aeroelasticity[40]

Table 1 (continued)

Table 2: The development history of unsteady aerodynamic calculation methods[40]

2.2 Classification of Aeroelastic Problems

2.2.1 Statics Problems

(1)Wing torsional divergence

Fig.1 shows the schematic diagram of a two-dimensional wing.The initial angle(initial angle of attack)is α0,and the wind speed is V.The angle of attack of the wing will increase under the action of the aerodynamic force (expressed as the aerodynamic torque to the rigid center), and the wing is subjected to a spring restoring torque at the same time.When V is not too large,the wing will reach the equilibrium at a new angle of attack α =α0+θ under the combined action of the aerodynamic torque and the spring restoring torque.Obviously,the additional angle of attack is the torsional deformation caused by the elastic support of the wing,and it is the torsional deformation that embodies the elastic body effect of the wing structure in the airflow(i.e.,the aeroelastic effect).If the spring stiffness is very large or the wind speed of the wind tunnel is very low, the torsion angle will be small.If the spring stiffness is very small or the wind speed of the wind tunnel is high,the torsion angle will be very large and even the torsion spring will be twisted to exceed it.The elastic limit causes damage.When the wind speed reaches a certain high value,the torsion angle of the wing will suddenly increase and even turn the wing over.This phenomenon is commonly referred to as“torsion divergence”,also known as“static divergence”or“divergence”[29].

Figure 1:Two-dimensional wing redrawn based on reference[29]

(2)Load redistribution

When the dynamic pressure is less than the torsional divergent dynamic pressure,the torsion angle and the aerodynamic lift are definite finite values and change with the change of dynamic pressure.This change will inevitably affect the flight performance of the aircraft,which is the redistribution of aerodynamic load.In the design of aircrafts,the possibility of the system static divergence instability should be avoided,but the redistribution of the load is inevitable.Engineers should design the initial airfoil that happens to be the cruise profile after loading.

(3)Reversal of the control surface

Fig.2 shows a two-dimensional wing with an aileron.The initial angle of attack is 0.Manipulate the aileron to make it deflect downward at an angle β.At this time,due to the deflection of the aileron,the aerodynamic force on the wing will change accordingly.The downward deflection of the aileron will increase the curvature of the wing;therefore,there is an increase in lift Lβ= (∂CL/∂β)βqS and its action point is generally located behind the aerodynamic center.If it is moved to the aerodynamic center,an additional torque Mβ= (∂Cm0/∂β)βqSc to the aerodynamic center will be generated,and usually Mβis a torque to lower the wing head.This means that ∂Cm0/∂β is a negative value.Since the wing is elastically supported,the torque will cause the wing to produce an elastic torsional deformation in the direction of reducing the angle of attack,which is to twist an angle θ in the direction of reducing the angle of attack.It will cause the wing to produce an additional downward aerodynamic force Lθ= (∂CL/∂α)θqS, and thus the total lift increment (usually referred to as lift increase) is ΔL = Lβ- Lθ.This effect of increasing and decreasing the actual lift caused by deflecting the aileron due to the elastic deformation of the wing is the so-called control surface efficiency (aileron efficiency)[40].

Figure 2:The statics model of a binary wing with control surfaces[40]

When the speed is lower than the reaction critical speed, the reduction in the aileron control efficiency can be expressed by the control efficiency.Assuming that the wing is rigidly supported,it is obvious that the lift increment generated after controlling the aileron deflection angle is only related to the angle.The lift coefficient is shown as Eq.(1):

The aileron control efficiency of a two-dimensional wing is defined as the ratio of the lift coefficient of the elastically supported wing to the lift coefficient of the rigidly supported wing,as Eq.(2):

where the aileron control efficiency η varies with the ratio q/qRof the air velocity pressure to the reaction critical velocity pressure,and its parameter is the ratio qR/qDof the reaction critical velocity pressure to the torsional divergence critical velocity pressure.It is worth noting that for a wing with determined parameters,and qR/qDis a fixed value:

(a)If qR＜qD,the control surface reaction occurs first,and then the torsional divergence occurs.

(b)If qR＞qD,the torsional divergence occurs first.

(c)IfqR=qD,then η ≡1 and the best handling efficiency is obtained.

2.2.2 Dynamic Problems

(1)Flutter

Flutter is a self-excited vibration with the nonattenuation amplitude due to the coupling of the aerodynamic force, the elastic force, and the inertial force when the aircraft exceeds a certain critical flight speed.Its essence is the dynamic instability of the aeroelastic system,and it often causes catastrophic damage[41–43].For example,in 1997,the first stealth fighter bomber F-117 crashed due to the flutter of the control surface.Currently,the wide application of lightweight flexible composites makes the structure of a new generation of aircrafts increasingly flexible.The structural plasticity of both the airframe and the aerodynamic rudder used for attitude control is greatly increased,which provides a sufficient internal cause for the occurrence of flutter.With the emergence of new technologies,the difficulty and breadth of aircraft flutter problems are indirectly exacerbated.Aircraft structure design and flutter suppression design will face great challenges.

(2)Gust response

Because the atmosphere is unstable,the aircraft cannot fly in the static air,and it is often affected by atmospheric turbulence.At this time, the aircraft may bear an excessively high load and cause a flight accident or seriously affect the flight performance and the safety and the comfort of the crew.There is no shortage of flight accidents caused by wind gusts in the history of aviation.Therefore,the aircraft strength and rigidity specifications stipulate that gust load calculations must be performed to verify whether the aircraft has sufficient strength and rigidity to ensure the safety and meet the operational requirements.At present, the flight experiment technology faces the challenges of great difficulty,long period,and extremely high cost.Various engineering algorithms and numerical analysis methods are mainly used to study the gust response.Gusts can be divided into vertical gusts and horizontal gusts.The disturbance motion caused by the vertical gust of the aircraft is limited to the longitudinal motion,and the disturbance motion caused by the horizontal gust of the aircraft is limited to the lateral movement,which can greatly simplify the gust response analysis.

(3)Buffeting

Wing buffeting is the dynamic response generated by the random aerodynamic excitation acting on the wing surface induced by the separation of airflow on the wing surface.The buffeting excitation comes from the generation of unsteady separated flows on the wings.It mainly includes the buffeting excitation at low speed and high angle of attack and the transonic shock wave buffeting excitation.There are three types of tail wing buffeting.The first type is the buffeting caused by the tail wing in the high-energy turbulence.Such buffeting can be caused when the tail is in the wake flow of the upstream fuselage, wing, or in the edge vortex rupture flow of the edge wing layout.The second type is the buffeting caused by the tail itself,which is similar to the buffeting of the wing at high angles of attack.At high angles of attack or sideslip,the airflow separation on the tail surface or the pressure pulsation caused by the shock-boundary layer interference when the flying Mach number is large can cause this kind of buffeting.The third type of buffeting belongs to coupled buffeting.When the frequency of other parts of the aircraft is close to the natural frequency of the tail wing,it is possible to cause the tail wing buffeting through structural coupling.The first type is the most important form of buffeting,which has the largest impact on the tail structure.

(4)Servo flutter

For the fly-by-wire aircrafts,the introduction of fly-by-wire flight control systems will cause a new flutter problem,namely servo flutter.The reason is that the modern aircrafts generally broaden the frequency band of the flight control system,especially for the aircrafts with active control functions.The frequency band of the controller often covers the frequency band of the low-order structural modal distribution,which will cause the high-frequency vibration response to be included in the sensor.Among the signals, it cannot be effectively attenuated when passing through the control loop, and even forms positive feedback,and then through the control surface and other actuating mechanisms to further excite structural vibration,thereby forming structural damage similar to flutter.

Since the actual aircraft is an elastic body, the sensor not only receives signals from the flight status but also receives vibration signals from wings and fuselage.These signals refed back to the rudder surface through the controller to form the high-frequency vibration of the rudder surface.This high-frequency vibration mode will be coupled with the original structural vibration mode,which may reduce the stability of the original aeroelastic system,resulting in servo flutter.

2.3 Bottlenecks in the Research of Aerodynamic Characteristics of Morphing Aircrafts

At this stage, there are three main types of flight dynamics models for morphing aircrafts: (1)the use of linearized small disturbance hypothetical flight dynamic model combined with the partial freezing method of the quasi-steady aerodynamic model; (2) the use of linearized small disturbance flight dynamic model for small-scale rapid morphing; (3) the flight dynamic model for large-scale slow morphing.The perturbed flight dynamic model is combined with the unsteady aerodynamics correction, which is suitable for large-scale and slow-speed morphing, considering the time change rate of inertia and the flight mechanics model of quasi-steady aerodynamics.For the case of the large morphing scale and the high speed,there is little related research.It is thus necessary to adopt a flight dynamic model that considers the time change rate of inertia combined with unsteady aerodynamic forces.The difficulty lies in the modeling of unsteady aerodynamic forces during the morphing process[44–46].The ultimate goal of the research on aerodynamic and flight mechanics characteristics during the morphing process is to achieve reasonable and effective flight control during the morphing process and ensure the safety of the flight.On the basis of the effective control of the flight mode,it is necessary to ensure the safe and effective control during the dynamic conversion process.At present,this part of the research work is still in the exploratory stage for modern morphing aircrafts and especially for the flight control during the morphing process[47,48].

To achieve effective control during the morphing process, the first prerequisite is to accurately predict and deeply understand the unsteady dynamic aerodynamic characteristics of the aircraft during the morphing process,and to obtain an accurate and effective aerodynamic model[49].This is very important for morphing aircrafts that exhibit strong unsteady nonlinear interference aerodynamic characteristics during complex morphing processes.However, this is also a great challenge for the current aerodynamic prediction technology[50].Although many research institutions and universities have carried out wind tunnel experiments, numerical simulations, and theoretical analyses on the aerodynamic characteristics of morphing aircrafts, there is still a long way to go before achieving a mature technology state[51,52].

3 Current Status of Aeroelastic Numerical Simulation of High-Performance UAVs

In the morphing process,the unsteady dynamic aerodynamic features lack acknowledged analysis and support from relevant aerodynamic basic theories [53,54].In terms of numerical simulation methodologies, the modern numerical simulation technology is capable of simulating the moving boundary’s unstable dynamic aerodynamic properties[55,56].However,there is no effective theory or verification and confirmation of wind tunnel experimental results,and their rationality and correctness lack solid support[57].

3.1 Morphing Structure Analysis

The foundation of aeroelastic numerical simulation is morphing structural analysis[58,59].There are two approaches for structural analysis: the analytical homogenization and the finite element method.The homogenization method[60]is based on partial differential equations with fast-changing factors.Its essence is to indefinitely magnify a specific spot in a heterogeneous material, revealing a periodic unit cell stacking structure at a meso-scale (the hypothesis of periodic structure is also the premise of the homogenization theory).The finite element method [61], often called finite element analysis(FEA)in practical applications,is a numerical method for solving partial differential equations.It is a high-efficiency and commonly used calculation method.It discretizes the continuum into a collection of several finite-size unit bodies to solve the continuum mechanics problem[62–64].

In order to obtain better performance of the deformed skin based on 2D thin beam elements,Ermakova et al.[65]proposed the shape optimization of corrugations(as shown in Fig.3).Considering that the ratio of the Young’s modulus of the elastomer coating to the Young’s modulus of the composite corrugated core is small, the elastomer coating is ignored in the contact area.In Abaqus 6.12, S4R elements are used for finite element analysis,and the homogenization method[60]is used to evaluate the accuracy of the corrugated plate buckling mode,as shown in Fig.4.Based on Castigliano’s second theorem, Dayyani et al.[66] proposed two analytical solutions for calculating the equivalent tensile and bending properties of the coated corrugated skin(as shown in Fig.5).Since the elastomer coating and the composite corrugated core are well bonded together and the two have the same displacement,the strain energy term of the contact between the elastomer and the glass fiber can be ignored.

Figure 3:Three typical shapes of corrugation[65]

Figure 4:Homogenization method to evaluate the buckling mode of corrugated plate[60]:(a)bending of coated corrugated skin based on thin beam elements; (b) out-of-plane global buckling on S4R elements

Figure 5:The coated composite corrugated core[66]

Dayyani et al.[67] derived a general-purpose super unit for a corrugated core unit with elastic coating that can be applied to morphing aircraft structures.Through the direct stiffness method and Castigliano’s second theorem,together with appropriate boundary conditions,the stiffness matrix of an arbitrary corrugation curve(as shown in Fig.6)is obtained.The matrix ignores the strain energy term of the elastic body in contact with the composite material and can be applied to small morphing of a 2D thin bending beam with the variable curvature.For beams with large deflection, geometric nonlinearity plays an important role.By introducing von Kármán strain,geometric nonlinear effects can be included.Assuming that there is no change in the displacement,strain and stress in the width direction of the beam,combined with the Ritz method,the principle of minimum total potential energy is employed to obtain its approximate solution.For plates such as folding wings[68],it is assumed that Kirchhoff’s plate theory is applicable.Compared with the von Kármán hypothesis,the plate rotation is negligible,and the structural equation of motion can be established by the principle of virtual work.

Figure 6:Schematic diagram of the corrugation unit cell and its curve function[67]

Structural analysis is combined with other system-level analyses such as aerodynamic analysis and control analysis[69,70].The structural equation of motion can be obtained by Newton’s principle,and the stability can be determined by using the Strut diagram based on Mathieu’s equation.Assuming that the structure under consideration is in a state of plane stress,the Lagrangian equation[71]can be used to derive the system equations,and structural nonlinearities can also be considered at the system level,such as piecewise linear stiffness and cubic stiffness.Commercial multibody code LMS Virtual.Lab Motion can be used to perform nonlinear dynamic simulation of flexible multi-body systems.

3.2 Aircraft Dynamics Modeling

The current dynamic modeling of morphing aircrafts can be divided into two categories:(1)the multi-rigid body dynamics model (such as Newton-Euler method [72], Lagrange method [73], and Kane method[74]);(2)the flexible-body dynamics model(such as floating frame of formulation[75]and absolute nodal coordinate formulation [76,77]).Compared with the multi-rigid body dynamics model, the flexible-body dynamics model is closer to the actual engineering conditions and can be applied to the aeroelastic analysis of folding wings,so it has received more attention.

In the calculation and evaluation of the aerodynamic characteristics of morphing aircrafts,Wickenheiser et al.[78]extended Prandtl’s lift line theory to a wing with arbitrary camber and chord length distribution and a non-ideal airfoil.Bowman et al.[79]found that at low speeds,the increase in the camber of the airfoil causes the point of least drag on the lift coefficient-drag coefficient curve to move to the lower right.In contrast, at high speeds, the increasein type camber causes the point of least resistance on the lift coefficient-drag coefficient curve to move to the upper right.Samareh et al.[80] developed an analysis tool using an efficient parameterized model formula, which can automatically generate models of aerodynamic parameters,geometric parameters,and shape change parameters for systems that undergo large deformations.Cho et al.[81] used the Panel method to solve the bionic double-wing problem and compared it with the CFD method.Gao[82]theoretically derived the unsteady aerodynamic characteristics of the two-dimensional airfoil during the morphing process through the analytical solution of the Zhukovsky.

Different from the aeroelastic analysis of conventional wings,the aeroelastic analysis of folding wings requires re-establishing the structural dynamic model and the aerodynamic model at each angle,and coupling them to obtain the aeroelastic model.For example,Snyder et al.[83]studied the flutter characteristics of the Goland folding wingand found that the hinge stiffness of the folding wing has a great influence on the flutter frequency and flutter speed of the morphing wing.Selitrennik et al.[84] proposed an efficient aeroelastic modeling method for rapidly morphing folding wing systems.This method uses virtual mass method to establish the structural dynamic equation of the system, and establishes the aeroelastic equation of the foldable morphing wing by coupling the dynamic equation and CFD technology.Lee et al.[85]proposed a numerical calculation method that can efficiently predict the subcritical limit cycle oscillation of a folding wing.Although the abovementioned aeroelastic dynamics modeling methods have been successfully applied to the aeroelastic modeling of the foldable variant wing,it is still necessary to carry out repeated as well as non-repetitive structural dynamics modeling for different folding angles in the aeroelastic analysis process.Zhao et al.[86] and Ni et al.[87] used the substructure synthesis and the dipole grid method to construct a parametric aeroelastic model of the folding wing (as shown in Fig.7), and analyzed the influence of modal damping, folding angle, and hinge stiffness on the dynamics of the folding wing.Huang et al.[88] proposed a parametric aero-servo-elastic modeling method, which efficiently obtained the aero-servo-elastic mathematical model of wings under different folding angles, and realized the design of control law and closed-loop aero-servo for active flutter suppression.Zhan et al.[89],based on the manifold tangent space interpolation method, established a folding wing parameter variable structure dynamic model with the folding angle as a parameter,and coupled the parametric unsteady aerodynamic model based on the dipole grid method to establish the parametric aeroelastic model of folding wings.

Figure 7:Parameterized aeroelastic modeling and flutter analysis for a folding wing[86]

3.3 Aerodynamic Calculation Method

Regarding the current mainstream CFD technology, whether it is solving steady or unsteady flow problems,its control equations are all implemented in the unsteady form through a certain time iterative algorithm[90–92].The main difference is that in the calculation of the steady flow field,it is the virtual time variable that has no practical significance, while the variable in the evolution of the unsteady flow field is the real time variable that has the exact physical meaning[93,94].

3.3.1 Steady Flow Field Solution Technology

The steady flow solution technology has a longer development process.It has an inherited relationship with the unsteady flow solution technology[95].Many steady flow field solution techniques are the basis of unsteady flow solutions.

The linear method is based on the theory of potential flow and is most suitable for thin lift surfaces with small angles of attack (such as cruise flight).The panel potential flow method is widely used to simulate the steady aerodynamic behavior in morphing wing design and optimization[96,97].This method uses singularity panels to replace wing boundaries,such as source panels,bimodal panels,and vortex panels.Lafountain et al.[98]used XFOIL to perform fast and relatively accurate 2D steadystate flow simulations for different morphing configurations,and used arc-controlled morphing wings to maneuver.Woods et al.[99]compared the linear method XFOIL with the open source CFD solver OpenFOAM,and found that the two have very similar aerodynamic performance prediction methods,but the computational cost is only a small part of the CFD solver.Fincham et al.[100]used XFOIL as the source of aerodynamic predictions,and found that the performance of a morphing airfoil(as shown in Fig.8)is almost close to an ideal airfoil whose shape can be changed arbitrarily.

Figure 8:Example of a NACA 0012 morphed to different cambers through the RBF method[100]

The extended lift line theory decomposes a 3D wing into a series of 2D wing shapes,which is a common selection method.Wickenheiser et al.[78]proposed an extension of the Weissinger method,which extended Prandtl’s lift line theory to wings with arbitrary curvature and chord distribution and non-ideal airfoil cross-sections.Ajaj et al.[101]presented the Zigzag wingbox concept(as shown in Fig.9) that allows the wing span to be varied by 44% (22% extension and 22% retraction), and used the Tornado code to predict the aerodynamics of morphing wings along the span.Li et al.[102]took the seamless leading edge and trailing edge control surfaces as the research object, and used the Tornado code to calculate the lift coefficient and pitching moment coefficient of the light flexible deformable wing.

Figure 9:Top-view of the Zigzag wingbox concept[101]

Static aeroelastic analysis methods include structural analysis and aerodynamic analysis[103,104].The structural analysis method usually takes less time than the aerodynamic analysis method.Due to cost considerations,the linear aerodynamic method is preferred when optimizing morphing wings for static aeroelastic analysis.Bilgen et al.[105]proposed a static aeroelastic modeling and optimization method for a variable camber wing.Two cascaded bimorph actuators were used on the top and bottom surfaces of the wing, and they were fixed behind the wing.Then based on XFOIL and ANSYS, a semi-solid variable camber airfoil coupling system was established.The optimization problem of 3D morphing wings is often dealt with by the 3D panel method or the extended lift line theory.Molinari et al.[106]optimized a tailless morphing aircraft driven by MFC based on the 3D panel method.The results showed that the rolling torque was increased by 54%,the pitch moment at 30 m/s was increased by 43%,and the pitch moment at 15 m/s was increased by 378%.

The CFD method is more suitable when the angle of attack is large and the flow shows a strong viscous effect.Woods et al.[99] compared the Open FOAM solver with the XFOIL solver for the morphing structure and found that the calculation results of the two methods werethe most consistent under the 8-degree angle of attack.One of the many difficulties in CFD analysis is the prediction of the change from laminar flow to turbulent flow,which may be related to the flow conditions of the designed morphing wing.Gabor et al.[107]compared the CFD calculation results of the morphing wing(as shown in Fig.10)with the wind tunnel test results,and found that numerical results of the pressure distribution were consistent with the wind tunnel test results, but they failed to accurately explain the change mechanism of drag in experiments.Although there are some computational difficulties, the CFD is often better than linear methods when performing aero bomb prediction analysis based on a limited number of grids.Spadoni et al.[108] studied various configurations in airfoil design based on static aeroelastic analysis methods, and used Euler equations for flow field calculations.Ismail et al.[109] studied the torsional bending deformation of perimeter reinforced(PR)-compliant wing (as shown in Fig.11).By solving the 3D incompressible Reynolds Average Navier-Stokes (RANS) shear stress transfer equation, the aerodynamic force was calculated and its aerodynamic characteristics were identified.Studies have shown that the PR-compatible wing can produce a higher nonlinear lift distribution than the PR or rigid wing type.Lyu et al.[110]used the CFD solver sumb to optimize the adaptive deformable trailing edge wing.

Figure 10:Layout of the morphing wing[107]

Figure 11:PR-compliant wing[109]

3.3.2 Unsteady Flow Field Solution Technology

The definition is based on the description of the Euler’s field.If the flow field parameters change with time, which means the derivative of the physical quantity with respect to time is not zero, then this flow can be regarded as unsteady [103,111].For the flow control, the unsteady characteristics are mainly reflected by the time derivative term.Therefore, in the calculation of the unsteady flow, although the numerical dispersion of the time derivative term is similar to that of the steady flow, the essence of time advancement is different.The sources of unsteady effects can be divided into three categories.The first one is that caused by the instability of the flow itself[112].Typical phenomena include the shedding and evolution of separation vortices,shock-boundary layer interference oscillations,etc.The second one is that caused by the pulsation of incoming flow,which mainly includes random gust response, pulsating jet flow control, change of incoming flow velocity and direction with time, etc.The third category is the unsteady effect caused by the morphing of the solid wall boundary.Typical phenomena include the turbulence during the flexible vibration of the wing, the multi-body separation turbulence, and the control surface maneuvering deflection turbulence.In fact, a specific unsteady flow is often the result of the combined effects of the above three types of sources.From the perspective of the complexity of the numerical method,the numerical simulation of the first type of problems is easiest to implement,and it is also the basis for the second and third types[113].Table 3 summarizes the aerodynamic methods applied to morphing wings.

Table 3: Aerodynamic methods applied to morphing wings[114]

Peters et al.[120] proposed a state-space aerodynamic theory of flexible airfoils, which can be combined with other dynamic codes for time advancement, eigenvalue, or control system design.This theory has been applied to the aeroelastic analysis of flexible airfoils [121].Murua et al.[122]numerically investigated the effect of chord-wiseflexibility on the dynamic stability of compliant airfoils.A classical two-dimensional aeroelastic model was expanded with an additional degree of freedom to capture time-varying camber deformations, defined by a parabolic bending profile of the mean aerodynamic chord (as shown in Fig.12).Berci et al.[123] explored the combined aeroelastic behavior and gust response of a flexible airfoil theoretically(as shown in Fig.13)through simulating the aerodynamic load by Peters’state space model.This method can also be applied to morphing airfoils[124].In addition,there are some other aerodynamic theories that take into account the flexibility of the wing [125,126].For example, Johnson et al.[125] summarized the systematic application of von Kármán and Sears’extended unsteady thin airfoil theory to general morphing airfoils.The theory proposed by Huang et al.[127]is limited to rigid body motion and can be extended to unstable aerodynamics of thin morphing airfoils.

Figure 12:Airfoil camber line with two spring-supported rigid-body(plunge and pitch)and one elastic(camber deformations)degrees of freedom[122]

Figure 13: Flexible airfoil: Cross-sectional schematic diagram showing key parameters and load distribution[123]

The unsteady vortex lattice method(UVLM)is often used to simulate the morphing process in the time domain(as shown in Fig.14[116]).It is assumed that the flow field is inviscid,incompressible,and non-rotational [128,129].In the case of morphing wings, the influence coefficient matrix will change, and must be modified and solved at each time step.Therefore, the UVLM can easily adapt to different wing geometries and morphing processes[115].Ghommem et al.[116]considered active shape morphing in optimizing the flying performance of flapping wings based on the UVLM.The results showed that the spline-based morphing,which requires more design variables to be specified,produces a significant improvement in propulsion efficiency.Ren et al.[117,127] combined the reduced-order UVLM and the Euler-Bernoulli beam theory with transient boundary conditions to develop a first-order state-space aeroelastic model.The results indicated that the critical flutter speed of a variable-span wing is very sensitive to the span length because it has a significant impact on the structural stiffness and aerodynamic characteristics.

Figure 14:Wing discretization and an example of a typical vortex element[116]

The double lattice method (DLM) is used to simulate aerodynamics in the frequency domain.The DLM is usually used as the first step in aeroelastic analysis[118,130,131].Castrichini et al.[132]studied a flexible wing folding device for load reduction.The commercial multi-body code LMS Virtual Lab Motion was used to build the structural model,and the aerodynamic influence coefficient matrix was evaluated by the DLM.Shrestha et al.[133] studied the influence of the change of the folding angle on the flutter speed and frequency of awing.The structural and aerodynamic models were created by using MSC.PATRAN Flight Loads.Via changing the folding angle of the outer wing independently of that of the inner wing, a proposal was made to enhance the flutter characteristics of the aluminum folding wing.Pecora et al.[134] were dedicated to studying the aeroelastic stable configuration of the morphing wing trailing edge driven by traditional electromechanical actuators.The unsteady aerodynamic influence coefficient matrix was evaluated by the DLM,and the structural model was established based on the finite element method.The results obtained showed that for the most practical combination of the trailing edge stiffness and the inertia distribution,flutter could be avoided if the actuation chain provides sufficient stiffness.The reduced-order method(ROM)based on the CFD method is very popular in dynamic aeroelasticity.Liu et al.[119] conducted a ROMbased aeroelastic continuous dynamic simulation of a morphing wing.The results showed that the LCO amplitude predicted by the ROM method is very consistent with the direct CFD method, but the phase shift would occur with the development of LCO.

4 Current Status of Aeroelastic Wind Tunnel Test on Morphing Aircrafts

Due to the unsteady aerodynamic characteristics caused by the shape change of the variant structure,great challenges are presented to the prediction technology of the aerodynamic performance and load environment of aircrafts.Compared with the traditional wind tunnel test, the morphing model design needs to consider the consistency of the morphing scale and the response time in the real flight environment.Considering that the characteristic size of the aircraft changes greatly after the morphing,the revision of the similarity theory to the test results requires the establishment of a new test simulation method,and the measurement and evaluation of unsteady aerodynamic characteristics under dynamic morphing conditions have also become a major difficulty.

4.1 Overview of Wind Tunnel Test

The wind tunnel test of the aircraft is to fix the actual aircraft or its reduced model in the ground wind tunnel test environment and blow the airflow at a certain pressure and wind speed,and to use the relative movement of the airflow and the test object to simulate the actual flight situation in the air [135].The purpose is to form the interaction between the airflow and the model, and then to study the aerodynamic characteristics of the actual aircraft.According to the test object used, it is divided into the physical test and the model test.The former will not be distorted, but it needs more test consumables,sets higher requirements for the test environment,and requires more stringent conditions.The latter is based on similarity theory, performing an experiment with a model that is geometrically similar to the real thing.At the same time, the two flow fields are related to similar criteria,such as Reynolds number,Mach number,and Prandtl number.It should be ensured that the flow field similar to the real flow field is used for testing.The dynamic performance of the aircraft is simulated and tested.Wind tunnels can be classified according to the airflow velocity.Those with an airflow velocity ranging from 0.8 Ma to 1.4 Ma are called transonic wind tunnels;those with an air velocity of about 0.4 Ma to 0.8 Ma are called subsonic wind tunnels[136,137].

In the design stage,the aeroelastic characteristics of the aircraft must be analyzed through wind tunnel tests to verify whether the strength design and performance indicators meet the requirements,and corresponding improvements should be made.The wind tunnel flutter test is carried out before the flight test, so that its results can be used to provide a basis for and reduce the risk of the flight test.In addition,the aeroelastic wind tunnel test also provides support for the verification of advanced aerodynamic/structural coupling design and analysis methods.Taking the development of advanced aircrafts in the United States as an example,the transonic dynamics wind tunnel(TDT)of the Langley Center of NASA was completed in the early 1960s.By 2016,more than 600 wind tunnel tests had been carried out.Most of them are aeroelastic tests, including fixed-wing aircrafts, rotary wing aircrafts,launch vehicles,and other configurations,which have contributed to the aeroelastic design and analysis of the research models, the engineering application of new layouts and new technologies, and the development of numerical calculation methods.Therefore,the aeroelastic wind tunnel test technology is a solid foundation for the research and development of aerospace vehicles.

The main vibration forms of the wind tunnel model can be decomposed into pitch vibration,the roll vibration, and the yaw vibration.The air pulsation (which should be referred to as velocity pulsation) is one of the vibration factors of the wind tunnel test and comes from two aspects.The main frequency of the air pulsation is generally low, and it is easy to be coupled with the low-order natural frequency of the support system.However,the original air pulsation is difficult to reduce.On the other hand,the airflow pulsation caused by the airflow separation phenomenon generated by the model during the test(such as the separation of the model airflow in the high angle of attack test)is unavoidable.

The vibration of the wind tunnel test is mainly improved from two perspectives.One is the improvement on the wind tunnel aerodynamics, and the other is the improvement on the dynamic characteristics of the model support system.It is very difficult to improve the excitation source from the aerodynamic aspect, especially in such cases as modifying the shape of the aircraft model or changing the flow field of the wind tunnel.Therefore,damping vibration by improving the dynamic characteristics of the model support system will be an effective way to solve this problem.According to the principle,the vibration suppression method of the wind tunnel tail support model can be divided into the passive control method based on a damper and the active control method based on the inverse piezoelectric effect of a piezoelectric element.

4.2 Similarity Law Theory

Similar phenomena refer to physical phenomena that follow the same physical laws and equations.Single-valued conditions are conditions distinguished from the same type of phenomena(conditions that determine the uniqueness of solutions to physical equations), and the similarity criterion is the relationship between certain physical quantities in similar phenomena.The similarity criteria of Reynolds number, Mach number, mass ratio, reduction frequency, etc., are derived from geometric conditions, physical conditions, boundary conditions, and initial conditions.Given the physical equations, we can use the equation analysis method (similar transformation method and integral analogy method)to derive the similarity criterion.If the physical equation is not established,we can then use the dimensional analysis method(Buckinghamπtheorem)to derive the similarity criterion.

The similarity theory requires that all similar parameters of the model in the model experiment must be completely consistent with the prototype, but it is almost impossible to completely meet the requirements of the similarity theory in the aeroelastic model experiment.Therefore, the local simulation method is adopted to simulate the similar parameters that control the main physical phenomena and processes,and the requirements for other secondary similar parameters are relaxed.The low-speed flutter test generally simulates Reynolds number,mass ratio,and stiffness ratio,whereas the transonic flutter test generally simulates Mach number,mass ratio,and stiffness ratio.

For the high-speed wind tunnel flutter test model,the first three basic scales are:length,density,and dynamic pressure.The length ratio is determined according to the experimental section of the wind tunnel and the actual size,the density ratio is determined according to the air density in the wind tunnel and in the actual flight state,and the dynamic pressure ratio is estimated according to the reachable limit of the stagnant point pressure in front of the wind tunnel and the actual object.The high-speed wind tunnel flutter test is much more difficult than the low-speed wind tunnel flutter test.The reason is that the scale of the high-speed wind tunnel test section is much smaller than that of the low-speed wind tunnel,the test time is shorter,and the strength requirements of the test model are stricter.

4.3 Design of Wind Tunnel Experimental Model

In recent years, in the field of structural optimization design, the respective advantages of comprehensive sensitivity methods and genetic algorithms have complemented each other, forming a sensitivity/genetic hybrid optimization design method with efficient global search capabilities.In the design process of the flutter model,the method can be applied not only to the final modification of the structure,but also to multiple links in the entire design process.Neal et al.[137,138]designed the wind tunnel experimental models of morphing aircrafts to achieve five goals: establishing a quasisteady model of large morphing aircrafts,optimizing research on‘efficient flight layout’,establishing transient models for fast-morphing aircrafts, conducting evaluation research on ‘morphing as a maneuver control device’,and simulating the flight control of morphing aircrafts.There are two types of morphing that can occur in the final experimental model.The first type is to adapt the aircraft to large deformations for different tasks,including changes in extension,sweep angle,and tail layout.The second type is to enable the aircraft to maneuver.For the controlled morphing of the wing including wing torsion,the combination of the two types of morphing can significantly enhance the morphing ability of the experimental model.Henry et al.[139]used wind tunnel experiments to study the effect of variable aspect ratio on the stability of UAVs.The calculated stability coefficient verified that the variable aspect ratio wing can maintain stability when the span change symmetrically, and the roll mode becomes more stable as the span increases.In addition,a wind tunnel experiment verified that asymmetrical elongation can generate a rolling moment that controls the rolling maneuver.Hall et al.[140] studied the possibility of improving the flight performance of micro-aircrafts by changing the sweep angle.High aerodynamic efficiency and the ability to launch quickly and easily were required.To achieve these goals,a variable swept-back wing was investigated.

An exceptionally fast morphing rate will have a great impact on the unsteady aerodynamic characteristics, stability characteristics, and flight control of the aircraft during the morphing process.Rapid morphing also sets higher requirements for deformable materials, structures, actuation systems, and morphing control.Conversely, a too slow morphing rate will bring additional costs and disadvantages to the deformable materials, structures, actuation systems, and flight control,and cannot meet the requirements of the flight mission and the environment when the overall flight performance is required to respond quickly.Therefore, it is necessary to comprehensively consider factors such as aerodynamics,flight control,materials,structure,and action,and carry out comprehensive calculations and experimental analysis,so as to obtain a realistic and feasible morphing rate.For example,Nextgen’s MFX-1 disclosed that the time required to deform the wing is 15 s,the F-111 variable sweep speed can reach 3.8°/s, and the F-14 variable sweep speed can reach 7.5°/s.It is still unclear whether the designer chooses such a morphing rate out of the consideration of the structural realization factors and flight mission requirements,or out of the intention to avoid causing unsteady aerodynamics and related flight control problems.However, in the preliminary analysis of the morphing of birds and variable-sweeping aircrafts, the safe morphing rate should have a certain relationship with the incoming flow velocity and the characteristic scale of the aircraft.

5 Conclusion

This paper has reviewed the development process and the current status of aeroelastic issues,numerical simulations,and wind tunnel test of morphing aircrafts.It aims at discussing the potential aeroelastic common problems in the morphing process, and comprehensively elaborates the key problems faced by the morphing aircrafts and the progress.In addition,the wind tunnel test technology has been analyzed in detail.The concept,origin,and development of aircraft morphing technology,aeroelastic numerical simulation, and aeroelastic wind tunnel test have been summarized, which provide background introduction and pavement for the follow-up review.

The aeroelastic problems and their solutions in the development of high-performance UAVs have been explained from three aspects: the development history of aeroelasticity, the classification of aeroelastic problems, and the bottlenecks of the aerodynamic characteristics of morphing aircrafts.The technical bottlenecks of the aerodynamic characteristics of morphing aircrafts are mainly embodied in the evaluation of aerodynamic characteristics, mechanism, and applicable unsteady dynamic aerodynamic modeling during the morphing process.Unsteady aerodynamic forces have an important impact on the handling characteristics of the aircraft,control law design,and flight safety.Due to this reason,it is necessary to study the unsteady aerodynamic characteristics,dynamic characteristics,and generation mechanism of the typical morphing methods during the morphing process, to establish the unsteady dynamic aerodynamic model,and to develop the analysis and evaluation method of the stability characteristics.

The current research status of aeroelastic numerical simulation of high-performance UAVs has been described from four aspects: structural analysis of morphing, modeling of flight dynamics,computational mesh deformation technology, and aerodynamic calculation methods.Traditional aeroelastic problems are divided into static aeroelastic problems and dynamic aeroelastic problems.They mainly focus on the stability of elastic structures under the combined action of the inertial force, the elastic force, and the aerodynamic force.With the application of the fly-by-wire control technology in aircrafts,the servo aeroelasticity has gradually emerged,and the research on hypersonic aircrafts and the thermo aeroelasticity has also been paid attention to.For high-speed aircrafts, the transonic problem is the speed range that all will face.Therefore,the transonic aeroelasticity is a matter of widespread concern.The military aircraft, micro-aircraft, and fluid machinery also involve the aeroelastic problem under high angle of attack separated flow.Both transonic aeroelastic problems and high angle of attack aeroelastic problems belong to nonlinear dynamic problems.The corresponding numerical simulations still have many aspects(such as more exact computing methods,less computing time, multi-physics simulations)to be studied and explored.In the future,these problems should be deeply investigated to promote the development of morphing UAVs.

The research status of the aeroelastic wind tunnel test of morphing aircrafts has been described from three aspects: the wind tunnel test, the similarity law theory, and the wind tunnel test model design.Due to the unsteady aerodynamic characteristics caused by the shape change of the morphing structure,great challenges are presented to the prediction technology of the aerodynamic performance and load environment of the aircrafts.Compared with the traditional wind tunnel test,the morphing model design needs to consider the consistency of the morphing scale and the response time in the real flight conditions.Considering that the characteristic size of the aircraft changes greatly after the morphing,similar theories have little effect on the test results.The correction requires the establishment of a suitable test simulation method,and the measurement and evaluation of unsteady aerodynamic characteristics under dynamic morphing conditions have also become a major difficulty.In the process of dynamic morphing, there are great difficulties in wind tunnel experiments.At present, there are more tests to examine steady or quasi-steady aerodynamic characteristics under different morphing states,but less to examine dynamic unsteady aerodynamic characteristics in the continuous morphing process.Therefore,testing dynamic unsteady aerodynamic characteristics is of interest and value for promoting the development of morphing UAVs.

Funding Statement:The authors received no specific funding for this study.

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.