Integrated Design of a Morphing Winglet for Active Load Control and Alleviation of Turboprop Regional Aircraft

: Aircraft winglets are well-established devices that improve aircraft fuel efﬁciency by enabling a higher lift over drag ratios and lower induced drag. Retroﬁtting winglets to existing aircraft also increases aircraft payload/range by the same order of the fuel burn savings, although the additional loads and moments imparted to the wing may impact structural interfaces, adding more weight to the wing. Winglet installation on aircraft wing inﬂuences numerous design parameters and requires a proper balance between aerodynamics and weight efﬁciency. Advanced dynamic aeroelastic analyses of the wing/winglet structure are also crucial for this assessment. Within the scope of the Clean Sky 2 REG IADP Airgreen 2 project, targeting novel technologies for next-generation regional aircraft, this paper deals with the integrated design of a full-scale morphing winglet for the purpose of improving aircraft aerodynamic efﬁciency in off-design ﬂight conditions, lowering wing-bending moments due to maneuvers and increasing aircraft ﬂight stability through morphing technology. A fault-tolerant morphing winglet architecture, based on two independent and asynchronous control surfaces with variable camber and differential settings, is presented. The system is designed to face different ﬂight situations by a proper action on the movable control tabs. The potential for reducing wing and winglet loads by means of the winglet control surfaces is numerically assessed, along with the expected aerodynamic performance and the actuation systems’ integration in the winglet surface geometry. Such a device was designed by CIRA for regional aircraft installation, whereas the aerodynamic beneﬁts and performance were estimated by ONERA on the natural laminar ﬂow wing. An active load controller was developed by PoliMI and UniNA performed aeroelastic trade-offs and ﬂutter calculations due to the coupling of winglet movable harmonics and aircraft wing bending and torsion.


Introduction
Morphing systems have been entered into aircraft design since the very first years of modern aviation [1,2]. Wright's Flyer [3] or the first concept of aircraft [4], had flexible wings to accomplish the different needs during flight (maneuver, climbing, dive, take-off, landing, etc.). However, the sudden operational speed increase, and the correlated load level rise, moved attention towards more rigid structural elements. Off-cruise operation solved. Indeed, as well as giving the achieved results, demonstrating a sustained actual flight with those novel structural systems, the aim of the study is to describe the entire design process, ranging from both low-fidelity and high-fidelity aerodynamic simulations and predicting the achieved performance of the winglet, to active load controller design and alleviation performance simulations, aeroelastic trade-off assessments and structural layout design and verification. As aforementioned, in fact, morphing systems should be effectively looked at as a new way to interpret the aircraft system design, overcoming the many trappings that have been formed over almost one century of aircraft engineering. Such an approach promises to be an established base for future developments in the field.

Reference Aircraft
The reference aircraft considered in this work is the one designed by Leonardo Company in the framework of the Clean Sky Green Regional Aircraft program. Such a regional aircraft is a 90 passenger turboprop aircraft, referred to as TP90 ( Figure 1) throughout this article. In the course of the Clean Sky 2 AG2 project, the TP90 A/C wing shape was redesigned by Onera in order to have extended natural laminar flow (NLF) abilities at cruise. The design considered multi-point optimizations of the airfoils for cruise, climb and low-speed flow conditions (see [14,15] for details about the NLF design). The resulting aircraft configuration is then referred to as AG2-NLF.

Numerical Methods Used for Aerodynamic Performance Evaluation
For the evaluation of the aerodynamic aircraft performances, the elsA CFD software [22] (ONERA-Airbus-Safran property) was used. This code solves the compressible three-dimensional RANS equations by using a cell-centered finite-volume spatial discretization on structured multi-block meshes.
For the spatial scheme, the one proposed by Jameson [23], is used for the conservative variables. A fourth-order linear dissipation κ 4 is generally used, with added second-order dissipation terms κ 2 for the treatment of flow discontinuities. The turbulence model used is the one-equation Spalart-Allmaras model with the QCR modification [24]. Multi-grid computations have been used for convergence acceleration.
An adaptation of the computational grid is necessary for the evaluation of the morphed winglet shapes. Among the different possible methods, the one used in in the SARISTU project [25] and described in [26], has been retained. This first considers a deformation of the surface grid of the wing, which then propagates into the 3D volumic grid002E.
Finally, it is necessary to compute the location of the transition line on the wing for an accurate estimation of the performance of the NLF wing. The elsA software has the ability to compute laminar flow regions and to determine the transition location by using the AHD compressible criterion for Tollmien-Schlichting instabilities and the C1 criterion for crossflow instabilities within the iterative convergence process of the RANS computation, as described in [27].

AG2-NLF Wing Equipped with Reference Winglets
In a second design step, the AG2-NLF wing was equipped with reference winglets, designed by the Leonardo Company, and new Karman fairings. The different design integrations made on the reference TP90 aircraft in order to create the reference AG2-NLF aircraft did not significantly affect the NLF characteristics found on the wing-alone computations. In a general overview, a laminar flow extension on about 50% of both the upper and lower surfaces is found on a large range of incidence, which satisfies the aerodynamic design constraints. On the upper surface (Figure 2), NLF flow is lost for α higher than 6.5 • roughly (C L higher than 1.0). On the lower surface ( Figure 3), a turbulent flow behavior is found at the wing tip at α = 0 • , but the laminar flow extends to nearly mid-chord for the other wing sections.   Figure 4 shows the computed aerodynamic performance of this reference aircraft wing, including winglets, at cruise conditions, for both turbulent and NLF flow conditions. In terms of performance, the NLF technology allows a global drag reduction of about 40d.c. at the flight C L (0.50), or an improvement of about 41d.c. However, note that the computations do not consider the nacelle and the propeller. Therefore, these values have to Appl. Sci. 2021, 11, 2439 5 of 27 be considered as the upper limit of the performance improvement, which will be observed on the total aircraft configuration.  The structural model of the complete aircraft was developed using NeoCASS [28]. Such a first guess included the preliminary mass and stiffness distribution on the basis of a simplified structural description of the aircraft, where the main structural components were modelled using beam elements. This model provides a useful starting point for the preliminary evaluation of the main aeroelastic properties of the aircraft and the effect of the morphing winglet on both the aircraft stability and aerodynamic loads. In addition, it enables the initial design of an active load alleviation system based on the availability of the morphing mobile surfaces on the winglet.
The structural sizing module of NeoCASS generates a model where the main structural components are represented by beam element, as shown in Figure 5. The model included both the baseline aircraft and the winglets; the winglet's mass and stiffness were included in the model to properly assess the related impact on airframe loads and aeroelastic stability, and no structural model was considered for the control surfaces, which were represented only through their estimated mass properties and as aerodynamic elements.
The structural mass and stiffness matrices are assembled by the aeroelastic analysis module of NeoCASS, which employs a Finite Volume formulation for the beam elements [29]. A modal analysis is then performed on the structural matrices to obtain a set of natural modes able to represent the dynamic response of the aircraft in free flight, such as where q is the set of natural modes and M and K are, respectively, the associated modal mass and stiffness matrices. The vector of generalized external forces, f a , contains the unsteady aerodynamic forces that will be discussed in the following paragraph.

Aerodynamic Model
To perform aeroelastic analyses it is also necessary to have a description of linear unsteady aerodynamic forces coupled with the structural response. This is obtained by the definition of an aerodynamic mesh, shown in Figure 6, which is used by the Doublet Lattice Method (DLM) [30] implemented in NeoCASS to obtain a set of frequency-dependent matrices. These matrices represent the generalized aerodynamic forces associated with the computed set of structural natural modes.  The final formulation of the generalized aerodynamic forces is where q ∞ is the dynamic pressure and v g is the gust velocity, defined as an additional forcing input for the system, along with the commanded deflection of the control surfaces δ c . The aerodynamic matrices H am (jω), H ag (jω) and H aδ (jω) are tabulated for discrete values of the frequency ω and can only be directly used in frequency domain analyses.

Morphing Winglet Concept
Besides significant aerodynamic benefits allowing reduced fuel consumption, morphing winglets have attracted growing attention in aviation because of their adaptive ability to lower wing-bending moments and increase aircraft flight stability in response to changing flight conditions. Several morphing winglet concepts have been patented [31,32], or are being developed [33,34], to alleviate gust loads and control the wing lift distribution over the wingspan through adaptive geometries. For the purpose of this research, the morphing winglet developed by CIRA is investigated for regional aircraft application. It consists of two "finger-like" mechanisms, shown in Figure 8, controlling the movement of movable surfaces (namely, upper and lower tabs), whose deflections are driven by dedicated electromechanical actuators [17]. By independently rotating the upper and lower surfaces, such an adaptive device is capable to both increase off-design aerodynamic efficiency by reducing induced drag and alleviate gust and maneuver loads by moving the surfaces, either synchronously or independently, to different angles. However, compared to more conventional plain flaps located at the trailing edge of the airfoils, the rear portion of the two morphing surfaces rotates upwards and downwards on two hinge lines of the finger-like mechanism. Such a morphing mechanism was already successfully validated on both full-scale morphing wing trailing edge [35][36][37] and aileron demonstrators [38][39][40].

Low-Fidelity Simulations
A low-order homebuilt solver (CIRA property) was used to optimize the winglet aerodynamic shape and to carry out the preliminary morphing winglet performance assessment. The Xavl code [41] is a 2.5D methodology which couples an inviscid 3D VLM solution carried out by the AVL public domain solver developed at MIT [42] with viscous 2D databases computed in a series of defined wingspanwise sections. The coupling between 3D inviscid solution and 2D viscous data is obtained by using the equivalent mean-line approach. The use of this low-fidelity code makes it possible to execute a full aerodynamic optimization with acceptable computational cost. The aerodynamic design was accomplished using the optimization chain described in Figure 9. The chain consists of an optimization tool (GAW), based on the Pareto dominance [43,44], the aerodynamic solver Xavl and a post-processor. Both the planform and the airfoil shapes of the winglet were concurrently optimized. The former aimed to maximize the aerodynamic efficiency in three different design points, including cruise, climb and climb, in one-engine-out conditions. The morphing airfoils' shape optimization had the additional goal of reducing the bending root moment in the most severe loading conditions. The winglet geometry was parametrized on a five-design section, shown in Figure 10. In each station, the GAW optimization tool is allowed to modify the sweep angle, the twist angle, the chord extension and the cant angle. The overall winglet height was not a direct design variable, but the GAW tool could modify the spanwise distance between each of the five-design section stations.

High Fidelity CFD Results and Performance Comparison at Design Points
In Figure 11, the resulting winglet aeroshape is compared with the reference one. It can be seen that the wingspan is slightly increased, and that its planform is more tapered than the reference one. Note that the same grid topology was used for the CFD computations and that the same reference wing area was considered for the aerodynamic coefficient. In a first step, it is necessary to evaluate the performance of the new winglet with no morphing function enabled. When considering the pressure field on the winglet, it can be seen ( Figure 12) that the load on the upper surface of the new winglet is more important than for the reference one at the nominal cruise conditions. Then, performances between the two winglets are compared at cruise conditions in Figure 13 for both fully turbulent and transition-free conditions. Black (turbulent) and red (laminar) curves are for the reference winglet, green (turbulent) and blue (laminar) curves for the new winglet. Performances in climb conditions have also been evaluated and similar trends have been observed. The C L (α) curve is nearly unaffected by the winglet change, and there is a slight improvement in the C D (or LoD) observed for the morphing winglet (see Table 1 for numerical values). A slight increase in pitching moment is also observed. As a conclusion, the new winglet was selected as the wing tip device for the AG2-NLF plane, as it exhibits some aerodynamic improvements when compared to the reference one.

Study of Morphing Winglet Configurations
Then, the use of morphing flaps adapted to this winglet has been assessed in term of possible aerodynamic performances. The adopted convention considers a positive deflection when it moves outward. For the present analysis, the same deflection applied at the same hinge location was considered for the upper and lower flaps. Only one "failure" case, with positive deflection on the upper flap and negative on the lower flap, has been considered. Figure 14 presents the different configurations considered by CFD evaluations for positive flap deflections. The angles considered were 3 • , 5 • and 10 • , applied either on hinge 1 (upper line) or hinge 2 (lower line). same hinge location was considered for the upper and lower flaps. Only one "failure" case, with positive deflection on the upper flap and negative on the lower flap, has been considered. Figure 14 presents the different configurations considered by CFD evaluations for positive flap deflections. The angles considered were 3°, 5° and 10°, applied either on hinge 1 (upper line) or hinge 2 (lower line). Considering the negative deflections ( Figure 15), only one case (−5°) was considered. Some problems occurred for larger negative values in the grid deformation tool, with local grid inversion leading to "negative cells" in the field, due to the concave topology used in the baseline grid. Finally, a kind of "failure" case was considered, with +10° deflection on the upper flap, and −5° for the lower flap, both applied on hinge 1. Computations were carried out around the CL value, corresponding to the design cases. Several flight conditions have been considered and numerical results have been sent Considering the negative deflections ( Figure 15), only one case (−5 • ) was considered. Some problems occurred for larger negative values in the grid deformation tool, with local grid inversion leading to "negative cells" in the field, due to the concave topology used in the baseline grid. Finally, a kind of "failure" case was considered, with +10 • deflection on the upper flap, and −5 • for the lower flap, both applied on hinge 1. Computations were carried out around the C L value, corresponding to the design cases. Several flight conditions have been considered and numerical results have been sent to partners for local force analysis on the different elements for sizing cases. In this paper, only the reference cruise (M = 0.52 @ 20,000 ft) and climb (M = 0.36 @ 15,000 ft) are presented. All the results presented in this document are at "aircraft" level. The force and moment coefficients considered the entire aircraft surfaces for pressure and friction integration. Figure 16 presents the C L (α) curves computed at cruise or climb conditions for the different morphed winglet configurations considered. The results obtained for the Reference initial winglet (black curve) and the new winglet (CIRA winglet) with no morphing (red curve) are given for comparison. These results show that the morphed winglet device can be used as a surface control on the aircraft for load alleviation. As a basic result, it can be seen that the larger the deflection or the flap length, the larger the ∆C L . It can also be noted that some combinations are nearly equivalent in terms of C L (α). For instance, applying a global deflection of 3 • on hinge 1 or 5 • on hinge 2 leads to the same ∆C L . In terms of improvements in the aerodynamic performance around the design flight point, the use of a morphed winglet seems very limited. Figure 17 presents the LoD obtained for the different configurations. At cruise conditions, there is no gain observed around the flight C L . Small positive deflections are nearly equivalent to the δ = 0 • case. For climb conditions, there are some very small gains observed in LoD found (about +0.05) for most of the positive deflections. Note that this approach underestimated the actual system benefits, as morphing effect was simulated by plain rotations at the first and second hinge lines of the tabs.  Figure 18 compares the "Moment diagrams" at cruise and climb conditions. In these diagrams, improvements are present when the absolute value of the moment coefficient decreases. For instance, this corresponds to a trend from the new winglet to the reference one. As expected, the use of negative deflections decreases the absolute values of the different moment coefficients. For positive deflections, an increase is observed, but with a small amplitude for +3 • (hinge 1 or 2) and +5 • (hinge 2) cases. Note that for a possible aero-structural application, the case with −5 • on hinge 2 could be used as a starting point.
There are indeed small losses in aerodynamic efficiency, but these are balanced with a global decrease in the absolute values of the moment coefficients.  Figure 18 compares the "Moment diagrams" at cruise and climb conditions. In these diagrams, improvements are present when the absolute value of the moment coefficient decreases. For instance, this corresponds to a trend from the new winglet to the reference one. As expected, the use of negative deflections decreases the absolute values of the different moment coefficients. For positive deflections, an increase is observed, but with a small amplitude for +3° (hinge 1 or 2) and +5° (hinge 2) cases. Note that for a possible aerostructural application, the case with −5° on hinge 2 could be used as a starting point. There are indeed small losses in aerodynamic efficiency, but these are balanced with a global decrease in the absolute values of the moment coefficients.

Active Load Alleviation Performance
Besides aerodynamic benefits, the morphing winglet design can exploit the movable control surfaces to carry out load alleviation tasks. Active load alleviation capability is obtained in this work by moving the winglet surfaces during the manoeuvre condition to reduce the resulting aerodynamic loads [45]. The main reason supporting this choice is the slow response of the winglet actuators which are designed to adjust the winglet deflection in slowly varying steady flight conditions, thus requiring a limited bandwidth. Higher actuators' bandwidth would have led to an unnecessary weight and power consumption increase.

Static Aeroelastic Analysis
A static aeroelastic analysis is initially used to evaluate the best-case system performance, that is the one that can be obtained neglecting the slow response of the winglet actuator.

Model Definition
Steady structural response is assumed in the formulation of the dynamic equations for static aeroelastic analysis [46], this means that only the degrees of freedom associated with the rigid body response will appear in dynamic equations, while the elasticity will provide a correction of the rigid body matrices, leading to a set of motion equations written as: where M bb and C bb are the rigid-body mass and damping matrices corrected with elastic and aerodynamic effects, v b is the vector containing the rigid body degrees of freedom, u c is a vector containing the control surface rotations and F c is the matrix defining the rigid body forces associated with the control surface rotation. F g is the weight force. No relationship between v b and . v b is assumed in the solution of the dynamic equations. The system of equation described above is underdetermined, and a set of additional relationships need to be introduced to fully define the manoeuvre condition, such as by imposing the angle of attack, the sideslip, the load factor or other quantities. All the constraints are expressed in the form As stated in [47], when there are more control surfaces available than unconstrained degrees of freedom it is possible to reformulate the manoeuvre problem as a constrained optimization problem, where a quadratic cost function of the form is minimized. The dynamic equations is then used as a constraint equation in the optimization, along with the other constraint equations defining the manoeuvre defined above and additional constraints that limits the control surface rotations, expressed as: where δ lower and δ upper represents the lower and upper limits for the rotation of the control surface.

Controller Definition
The optimization system defined in the previous section was used to get the controller gains. The cost function was based on the minimization of the wing root bending moment and the trim solution then contained presented the combination of main control surface rotations and winglet rotation leading to the minimization of the bending moment at wing root, while satisfying the manoeuvre definition.

Simulation Results
A steady level flight condition is considered at sea level with Mach number M ∞ = 0.41. A static pull up manoeuvre with load factor N z = 2.5 is used to evaluate the controller performance. The results in terms of structural loads are presented in Figure 19 considering the out-of-plane bending moment and normalizing the value with respect to the value at wing root. A reduction in the structural loads was obtained, with a limited reduction at wing root, and a very significant reduction at wing tip. However, since retrofitting winglets to existing aircraft leads to increased structural loads at the outboard section of the wing, a controller capable to reduce structural loads at the outer wing sections is of particular interest, even if it provides only a limited alleviation in the inboard wing sections.

Aeroelastic State-Space Model of the A/C
Dynamic effects can be important in the definition of the structural loads in manoeuvres conditions, for this reason a full evaluation of the MLA performances requires also a dynamic evaluation of the manoeuvres.
In addition to dynamic aeroelastic effects, also the dynamic response of the actuators and the digital implementation of the controller may affect the results. The slow dynamic response of the actuators is reducing the effectiveness of the controller, since it delays the deflection of the winglet surfaces with respect to that of the main control surfaces, and a nonlinear dynamical model is used to reproduce this behaviour.
To consider all the aforementioned effects, a nonlinear time simulation of each manoeuvre is performed. The frequency-domain based aerodynamic forces described above cannot be used directly in a time simulation, but an approximated system is required. The identification method described in [48] is then used to express the system in statespace form .
The aeroelastic system is then expressed in state-space form as .
where the output vector y ae contains the parameters describing the rigid motion of the aircraft and the internal loads along the wing span.

Winglet Actuation Model
The aeroelastic system is coupled with a dynamic model of the winglet servoactuators. The servoactuator system includes the definition of both the physical actuator and the control logic, its input is the commanded winglet rotation while the output is the actual winglet rotation δ c . The actual winglet rotation is also affected by the aerodynamic and inertial hinge moment, which need to be provided to the actuator model as additional output of the aeroelastic system.
Each of the two movable portions of the winglet are moved by an Electromechanical Actuator (EMA). A numerical model of the EMA is connected to the system, accepting as an input the tension command from the servocontroller, the actuator shaft loads and the actuator position. The servocontroller is then used to drive the EMA to the commanded position.

Time Simulation Results
The manoeuvres analysed with a static aeroelastic analysis were reproduced also using a dynamic simulation. A step deflection of the elevator is used to get the target load factor of 2.5g. The time history of the command is shown in Figure 20 and includes both a negative and positive command, to have both pull-up and pull-down manoeuvres. The dynamics of the elevator actuation system is considered by introducing a 10 Hz bandwidth and a saturation for the deflection rate. The slow winglet response prevents it to reach the saturation limits, and the deflection is not identical in the two sections due to the different sizes that lead to different aerodynamic hinge moments and then to different actuator responses. The larger surface of the outboard winglet section, in particular, leads to a lower rotation. The time history of wing internal loads shows that the MLA system is able to reduce structural loads even if the deflection values predicted in the static aeroelastic analysis are not reached in dynamic analyses. This is shown in Figure 21 where the bending moment both at wing root and at wing tip are shown. In the time histories the time has been nondimensionalised considering as a reference the period of the first wing structural mode.

Preliminary Structural Layout
The structural layout of the morphing winglet consists of a passive and an active part. The former is made of laminate skin panels and a torsion box consisting of spars and ribs. The two morphing tabs incorporate a "finger-like mechanism", shown in Figure 22 Winglet morphing is enabled by the relative rotations of three adjacent blocks, which are free to rotate around the hinges on the camber line, thus physically turning the camber line into an articulated chain of consecutive segments. Linking rods, hinged on non-adjacent blocks, force the camber line segments to rotate according to a specific gear ratio, which complies with the shapes that need to be achieved. Each tab is, therefore, a single-degree-of-freedom (SDOF) system; if rotation of any of the blocks is prevented, no change in shape can be obtained. On the contrary, if an actuator moves any of the blocks, all the other blocks follow the movement accordingly.

Aeroelastic Trade-Off
As for any safety-critical device, a preliminary assessment of morphing winglet on aircraft aeroelastic stability was judged to be extremely relevant to avoid the consolidation of an unsafe, and therefore unflyable, structural layout.
Flutter analyses were carried out in two consecutive phases; at first, trade-off investigations were addressed to establish the maximum allowable weight of the winglet structure.
To this aim, the entire device was modelled as a lumped mass placed at the wing tip; starting from the typical mass value of a conventional arrangement, the winglet mass was progressively increased until bending-torsion flutter of the wing was induced within the flight certification envelope (i.e., at a flight speed lower than 1.25, the aircraft dive speed VD). Bending/torsion flutter speed was estimated by referring to the Molyneux equation [49] and to the wing's natural frequencies resulting from the model described in Section 2.4.1; both symmetric and anti-symmetric coupling mechanisms were considered.
It was found that only the flutter due to the coalescence of anti-symmetric wing bending and torsion was sensibly affected by the winglet mass increase ( Figure 23); in order to assure aircraft flutter clearance, the overall mass of the imposed morphing winglet did not exceed the value of 90 Kg.
The preliminary structural design of the morphing winglet (Section 5.1) was carried out under this relevant inertial constraint and, once completed, a new loop of analyses was performed to verify the clearance from any other type of flutter induced by unstable couplings of winglet and wing modes.
The finite element model of the morphing winglet described in Section 5.3 was condensed into a dynamically equivalent set of beams, which were integrated into the overall aircraft aeroelastic model to replace those used for the simulation of the conventional winglet (Section 2.4.1). Winglet movables were also modelled through equivalent beams, and the rotation around their longitudinal axis was constrained by means of grounded springs reproducing the stiffness of the actuators in power-on mode. All flutter calculations and the elaboration of Vg plots were made in SANDY software environment [51].
A flutter due to the coupling of winglet movable harmonics and winglet bending/torsion was found at a critical speed of nearly 200 m/s ( Figure 24). The flutter dynamics were proven to be induced by a typical underbalanced behaviour of the movables. To increase the margin of safety with respect to flutter, the degree of mass-balancing of each movable was then increased up to 0.7, which successfully led to a first critical speed of nearly 210 m/s ( Figure 25). In such a configuration, a further flutter speed was also detected; nevertheless, since it was much higher than 200 m/s, it was considered not critical for the safety of the reference aircraft.

FE Model
An FE model was developed to fully represent the 3D CAD of the morphing winglet device. Here, the focus is on the passive structure of the winglet, designed to withstand the overall structural loads. The main substructures are: CFRP upper and lower panels and internal CFRP spars and ribs. Such components were FE modelled by 10,593 2D quad elements, shown in Figure 26, derived from the middle planes of the components. The upper and lower skins are composed of fabric layers. UD-plies in x direction are used for reinforcement. The spars and ribs are composed of fabric layers. The thickness of the laminates is between 2.7 and 3.6 mm. The weight was well below the preliminary target of 50 kg.
A number of load cases were considered for the structural design of the morphing winglet, in compliance with the CS25 (manoeuvres, gust, landing, engine failure) airworthiness regulation. Such loads were intended to be limit loads, i.e., the loads that the structure must be able to support without detrimental permanent deformation.
For the analysis of the laminates, a maximum/minimum strain criterion in 11-and 22-direction was used. Additionally, stress was also evaluated in 11-, 22-and 12-directions, but only to check. The ply stresses and strains for the failure analysis were read from the FEM results. The reserve factors for each ply of an element were then calculated. The enveloped results were then elaborated to allow for computation of the minimum reserve factors. An example of maximum principal strain computed on the skin for one of the sizing cases is shown in Figure 27.

Actuation Layout
A further step of the integrated design of a morphing winglet involves the preliminary assessment of the actuation system, including both the installation and performance requirements. It is worth noting that both influence the effective reproduction of the given aeroshapes under aerodynamic loads.
In light of the limited space available within the winglet box to accommodate both the winglet actuators and actuation chains, several installation trade-off studies were carried out to achieve the best configuration. A linear electromechanical actuation was considered for each mechanical drive of the two independent surfaces. Such an approach assumes that, in principle, a flight-worthy actuator of adequate size, weight, and power will be developed to withstand the hinge moments due to the aerodynamic loads.
The actuation chains addressed in this paper are controlled by two separate actuators, one for each surface, positioned at different winglet heights [17] (see Figure 28). Such a direct actuation layout, based on a crankshaft design, has the advantages of relative installation simplicity, lower weight, and higher stiffness, with additional benefits in terms of assembly time and maintenance.
Parameter studies and motion optimization were studied through a multibody model developed in MSC/Adams. The upper morphing surface deployment and its comparison with the given aeroshape is shown in Figure 29. This allowed for a comparison of the actual rotation of the tab with the nominal (expected) values enabling the load control function. A static actuation law was then implemented. As a result, the actuator rod stroke was correlated with the resulting tab rotation. Both functions are reported in Figures 30  and 31, respectively. For the given geometrical parameters, the actuators' performance requirements were then derived (see Table 2). Operating and maximum static axial loads were about 3.5 and 5 kN, respectively. The maximum speed foreseen for the load control function was in the order of 5 mm/s.

Conclusions
A morphing winglet has the potential to combine aerodynamic and structural benefits more efficiently than a conventional fixed wingtip surface. A fault-tolerant concept based on two individual (asynchronous) morphing surfaces (upper and lower) was investigated in this work in order to enhance aerodynamic efficiency and reduce maneuver loads on a regional transport aircraft.
A low-fidelity design tool was used to preliminarily design the winglet aeroshape and predict the related aerodynamic benefits without moving the adaptive tabs. Then, high-fidelity CFD analyses were proven to provide aircraft with reduced induced drag compared to the reference counterpart.
Different morphed aeroshapes were considered. It was found that such a device can be used as control surface on the aircraft for load alleviation, as a deflection of the tabs acts directly on the CL(α) curve. Concerning the load control, limited improvements were proven. However, it was found that some configurations may be used for a global aero-elastic optimization as a decrease in aerodynamic performance could be balanced by an improvement in moment coefficients at aircraft level, and, therefore, wing structure mass benefits could be obtained. The combined use of morphing devices adapted to the wing configuration equipped with the morphing winglet leads to about 3% drag reduction at climb conditions with respect to the reference winglet shape with no morphing applied.
Attention was also paid to the design and performance of a load-alleviation controller driving the hinged tabs to reduce wing loads in maneuvering conditions. A feed-forward architecture based on a direct measurement of the deflection of the main control surfaces was considered. The results demonstrate that the proposed mechanical system is able to reduce the wing loads in different flight conditions (including the sizing ones). Additionally, although the actuator dynamics reduce the performances that can be obtained with the MLA controller, the system was shown to reduce the wing structural loads.
Finally, structural analyses were carried out to confirm the actual feasibility of the device and to preliminarily estimate its weight, which was well below the preliminary target of 50 kg. For structural verification, the worst load cases within minimum and maximum internal forces (bending and torque) were considered. An FE model, with its geometry and material properties, was processed to size the laminate skin panels and the winglet torsion box made of spars and ribs. Trade-off aero-elastic assessments were carried out to predict the flutter speed with different winglet mass values. Actuation kinematics were studied through a multibody model by correlating the actuator rod strokes with the resulting tab rotations.
Further studies are planned to assess the behavior of the structure by considering the actuation chain's contribution to the overall stiffness of the morphing device. A structural optimization process is also envisaged with respect to a set of failure criteria. Because of the importance of aeroelastic phenomena, further development will focus on aeroelastic stability margins and the safety-related implications of actuator failures.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
Airgreen 2 project AG2-NLF Regional aircraft of the AIRGREEN2 project equipped with the NLF wing CFD Computational Fluid Dynamics CS2 Clean Sky 2 C L Lift coefficient C D Drag coefficient EMA Electromechanical actuators LCA Load Control and Alleviation LoD Lift over Drag ratio MLA Maneuver Load Alleviation NLF Natural Laminar Flow TP90 90-passenger Regional Turboprop Aircraft