Global aero-structural design optimization of composite wings with active manoeuvre load alleviation

In the scope of the DLR project VicToria (Virtual Aircraft Technology Integration Platform), an integrated process for aero-structural wing optimization based on high fidelity simulation methods is continuously developed and applied. Based upon a parametric geometry, flight performance under transonic flight conditions and manoeuvre loads are computed by solving the Reynolds-averaged Navier–Stokes equations. Structural mass and elastic characteristics of the wing are determined from structural sizing of the composite wing box for essential manoeuvre load cases using computational structural mechanics. Static aeroelastic effects are considered in all flight conditions and active manoeuvre load alleviation is integrated in the process. Global aero-structural wing optimizations are successfully performed for wings with and without active manoeuvre load alleviation. The active manoeuvre load alleviation is introduced with a simplified modelling of control surface deflections using a mesh deformation technique. The minimization of the fuel consumption for three typical flight missions represents the objective function. Wing optimizations are performed for variable and constant wing planform parameters as well as for wings with conventional composite wing box structure and for more flexible wings. The latter is accomplished by introducing modifications of the structural concept and the strain allowable. A significant mass reduction of the optimized wing box is obtained for wings with active manoeuvre load alleviation, resulting in a drop in fuel consumption of about 3%. For wing optimizations with the more flexible wing concept, the active manoeuvre load alleviation shows an additional reduction of the fuel consumption in the order of 2%. The wings with active manoeuvre load alleviation results in optimized wing geometries with increased aspect ratio and reduced taper ratio.


Introduction
The environmental impact of commercial aviation increases with the rapid growth of air travel and the CO 2 share of aviation will increase due to the increase of renewable energies in other transport sectors and in industry in general. For environmental protection and conservation of resources, the main goal of the aeronautical research in Europe and in the United States of America is a strong reduction of the CO 2 emissions per passenger kilometre [1][2][3][4].
To achieve this challenging goal, the development time scales for new technologies have to be reduced significantly. In this context, the methodologies and processes for physics based aircraft design and optimization have to be improved. Furthermore, an assessment of new technologies with consideration of all relevant disciplines and their interactions on overall aircraft level will be essential in the future.
The efficiency of commercial aircraft is determined by aerodynamic performance in terms of lift to drag ratio, aircraft empty mass and thrust specific fuel consumption of the engine. For the accurate drag prediction under cruise flight conditions the flow physics of transonic and turbulent flow can be taken into account using RANS-based computational fluid dynamics (CFD). To reduce the structural mass composite materials like carbon fibre reinforced polymers (CFRP) have been introduced in aircraft manufacturing. The corresponding structural concepts and sizing criteria have to be considered in the structural analysis and sizing process using structural mechanics solvers based on the finite element method (FEM).
Within the aero-structural wing optimization the optimum trade-off between the aerodynamic performance and the wing mass is achieved through combining high fidelity methods for numerical flow simulation of the aircraft outer shape and structural sizing of the wing box with an appropriate optimization algorithm. Thereby, the interaction of aerodynamic forces and wing deformations have to be considered for accurate flight performance and static manoeuvre loads prediction using fluid-structure coupling.
Wing design and optimization is a multidisciplinary task with a lot of practical constraints. For example, the size of the tail has to fulfil all stability and control constraints and thus depends on center of gravity range and wing geometry. Furthermore, the landing gear integration and the space allocation for the control surfaces including their actuators have to be considered. Neglecting the landing gear integration leads to unrealistic optimization results as shown in [5].
The technology of active manoeuvre load alleviation (MLA) reduces the aerodynamic loads using trailing edge control surface deflections to adapt the lift distribution under manoeuvre flight conditions. This technology has been published by White [6] for example and successfully applied to the Lockheed C-5A [7] and the Lockheed L-1011 [8]. In modern airliners, the manoeuvre load alleviation functions are an integral part of the flight control system. To reduce the loads due to atmospheric disturbances in flight, active gust load alleviation systems have been developed. An overview of applications of active control technologies for gust load alleviation has been given by Regan and Jutte [9]. The potential of manoeuvre and gust load alleviation to reduce the fuel burn and the direct operating costs has been investigated by Xu and Kroo [10] on aircraft conceptual design level. The active and passive load alleviation technologies have to be integrated into the sizing process of the aircraft structure and result in longer maintenance intervals and mass reductions. Passive load alleviation technologies use specifically designed wing geometries and stiffness characteristics to reduce the loads due to aeroelastic deformations. The anisotropic material properties of carbon fibre reinforced plastics allow improving the aeroelastic behaviour. With the industrial utilisation of automated processes for fibre placements, new technologies in the field of unconventional composites featuring spatially varying tow orientation become possible [11]. Therefore, a physical modelling of active and passive manoeuvre load alleviation technology is pursued in preliminary aircraft design.
With increasing knowledge of composite materials further mass reduction potentials can be exploited by better adaptation of fibre direction to internal loads, introduction of advanced structural concepts, and new manufacturing processes. The more flexible wing concept is a result of mass reduction due to new structural concepts with increased strain allowable and applied to the current generation of aircraft from Boeing (Boeing 787 and Boeing 777-8/9). In addition, the passive load alleviation due to static aeroelastic effects leads to further wing mass reduction. With the increased wing deflections of more flexible wings, the geometric non-linearities affect the internal loads [12] with increased aspect ratio. These geometric non-linearities begin to play a role in current wing design.
The enhancement of the aeroelastic design by introducing the interactions with the active flight controls is known as aeroservoelastic design. The fundamental aspects of aeroservoelastic analysis, design and optimization has been summarized by Livne [13]. Applications of aeroservoelastic optimization including active flutter suppression of metallic transport aircraft wing box structures have been shown by Stanford [14,15] for example. In the work published by Binder et al. [16], the interactions of manoeuvre load alleviation, gust load alleviation and aeroelastic tailoring of the composite wing box structure have been investigated on the basis of different aeroservoelastic optimizations. Thereby, the outer wing shape and the topology of the wing structure have not been optimized and the deflections of the spoilers and ailerons are used for active manoeuvre and gust load alleviation. The results show that 95% of the maximum achievable mass reduction of a generic long range transport aircraft configuration can be achieved with the combination of manoeuvre load alleviation and aeroelastic tailoring of the wing box structure.
Improvements in automation and coupling of accurate simulation methods in combination with advances in numerical optimization strategies lead to the emergence of multidisciplinary design optimization (MDO) based on high fidelity methods. Examples for the development and application of modern MDO frameworks are the publications of Sgueglia et al. [17] and Papageorgiou et al. [18]. In addition to the performance indicators based on flight physics, the economic aspects have been introduced in the context of MDO [19,20].
The challenge in using MDO based on high fidelity methods is the large number of design parameters and constraints and the increased computing effort. To overcome this issue, the adjoint method enables the efficient calculation of the flow variable gradients as a function of the design parameters for gradient-based optimization [21,22]. Up-to-date applications of the adjoint approach for multidisciplinary wing optimization have been shown in the publications of Kenway and Martins [5]; Liem, Kenway and Martins [23]; Keye et al. [24] and Abu-Zurayk et al. [25]. These publications show that the gradientbased optimization using the adjoint approach is an adequate method for multidisciplinary wing optimization with high fidelity simulation programs and a large number of design parameters.
In this work, an alternative MDO approach is introduced for cases in which gradients cannot be computed efficiently for all relevant disciplines. This applies particularly to cases which involve CAD modelling and structural sizing of composite structures using proprietary codes. Additionally, the gradient-based approaches have to be reconsidered for flows near the maximum lift including flow separations and the usage of active control surface deflections. Furthermore, a certain degree of flexibility in the process architecture and optimization strategy is desired. Especially the option to use optimization strategies seeking for the global optimum is important.
The focus of this publication is the investigation and assessment of active manoeuvre load alleviation technology in the context of multidisciplinary wing design. This includes the application of active manoeuvre alleviation to more flexible wings.

Methods
In the DLR project VicToria [25][26][27][28][29][30], an integrated process for aero-structural wing optimization based on high fidelity simulation methods is continuously developed and improved. A detailed description of the original process chain and their successful application is published by Wunderlich et al. [31,32]. The improvements relate to the introduction of grid deformation techniques for large geometry changes and simplified control surface deflections. Further extensions include a landing gear integration, a tail sizing based on handbook methods and a trim drag estimation functionality.
The MDO architecture of the integrated process chain falls in the category of MDF optimizations (Multi-Disciplinary Feasible) and can be described as ASO (Asymmetric Subspace Optimization) according to Martins and Lambe [33].
A detailed description of the process chain is outlined in the publication of Wunderlich et al. [34] and only the top level is presented here again.

Process chain for aero-structural wing optimization
The process chain applied is illustrated in terms of XDSMdiagrams (Extended Design Structure Matrix) [35]. Each component in the diagram receives input data in vertical direction and provides output data in horizontal direction. Input and output data are marked by parallelograms. Thick gray lines show the data flow and thin black arrows indicate the process flow. The numbering system defines the order in which the components are executed.
The flow chart of the process chain for aero-structural wing optimization is shown in Fig. 1. In every optimization step, the geometrical aircraft description of the baseline configuration is recalculated and updated in accordance with the current values of the design parameters. The resulting aircraft geometry is transferred to the subsequent simulation programs by using the Common Parametric Aircraft Configuration Schema (CPACS) [36,37]. The recalculation of the aircraft geometry includes the wing positioning relative to the fuselage, the integration of the main landing gear, the sizing of the tail and the check of geometrical constraints.
In the next step, the parametric CAD model is updated, the aerodynamic volume mesh is deformed and the structural model is generated. The parametric CAD model has been built in the commercial software CATIA ® V5, which enables accurate surface representation and robust and time efficient geometry changes.
Within the CFD volume mesh deformation process, the mesh representing the baseline configuration is deformed in two stages for all flight conditions in parallel. In the first stage, the geometrical changes between the baseline and the current geometry are computed based on the corresponding structured multi-block (SMB) surface meshes. The latter result from the automatic surface mesh generation and have an identical mesh topology with the same number of points. For automatic surface mesh generation, the commercial software Pointwise ® is applied. In the second stage, the control surface deflections are taken into account. Corresponding to the control surface deflection to be produced, the surface mesh displacement field is computed for each flight condition. It is propagated to the CFD volume mesh using the Elasticity Analogy (EA) mesh deformation method [38] available in the FlowSimulator [39][40][41] environment.
For the generation of the structural model the DLR inhouse tool DELiS (Design Environment for thin-walled Lightweight Structures) [42] is used. Based on the central data format CPACS, DELiS automatically generates  [43]. The finite element model is made up of shells elements enriched with physical properties of the wing spars, ribs and skin cells and finally exported for the commercial FE solver MSC Nastran ™ .
The fluid-structure coupling loop is marked with a rounded yellow box and the values of the design mission lift-to-drag ratio, the wing mass and the objective function value are evaluated for the convergence examination. The fluid-structure interaction belongs to the category of loosely coupled analysis [44,45], with the main difference of replacing the structural analysis of a sized wing structure by a combined structural wing analysis and sizing process. The integration of the structural sizing process into the fluid-structure coupling loop reduces the number of iterations by introducing a damper like behaviour.
For all flight conditions the aerodynamic forces and coefficients are computed using RANS-based CFD simulations. The flow simulations are performed using the DLR TAU-Code [46,47] which is integrated in the HPC framework FlowSimulator [39]. The solver's capabilities with respect to accurate flow predictions, also in near off-design regions, have been demonstrated in numerous publications, including those of the AIAA Drag Prediction Workshop Series [48]. The approach ensures that flight performance under cruise flight conditions and selected manoeuvre loads with consideration of flow separations in the presence of control surface deflections are analysed accurately and efficiently.
Based on the aerodynamic loads computed for the flight conditions considered, the wing-box structure is sized. Within the structural analysis and sizing process the objective is to fulfil the structural constraints in terms of failure criteria and converge the margins of safety (MoS) and wing mass. Hence, the structural analysis and sizing process represents a subspace optimization, which is described in detail in the publication of Wunderlich et al. [34]. Different design criteria are applied to ensure a valid structural design. As proposed by Dähne et al. [49] for stiffened panels, the criteria for strength, maximum strain and local and global buckling are used for skin and all stringer components. The main results of this process are the wing mass and the deformed wing shapes for the flight conditions considered. The structural analysis and sizing process uses the linear analysis of the commercial software MSC Nastran ™ for computing the internal loads and stresses. The commercial software HyperSizer ® is applied for sizing the composite wing box.
The structural deformations form the input for the CFD volume mesh deformation. The mesh deformation method based on radial basis functions (RBF) [50] available in the FlowSimulator is used. Afterwards, the objective function is evaluated and the convergence criteria of the static aeroelastic analysis are examined. In Table 1 a list of all considered physical quantities and their corresponding convergence criteria is shown.
The selected values represent an appropriate trade-off between accuracy and computing time for the static aeroelastic analyses. Once convergence of the fluid-structure coupling loop is reached, the objective function value is given to the global optimizer.
After the optimization run has finished the optimized vector of design parameters represents the main result of the process chain for the corresponding optimization problem.

Global optimization strategy
The evaluation of the objective function with the introduced process chain requires a relatively high computational effort. Derivative information is not available. A survey of global optimization methods for such problems is given by Vu et al. [51]. In this publication optimization methods are presented, which construct surrogate models or meta models of the objective function and constraints as a function of the design parameters. These surrogate models are then used to find the global optimum. The sampling strategy for the design of experiments (DoE) and the infill sampling criterion influences the efficiency of the optimizer [52]. For high dimensional constrained problems, the computing time for building the surrogate models increases and alternative approaches for building the surrogate models have been published by Bouhlel et al. [53] for example.
For the wing optimizations in this work an in-house surrogate-based optimization (SBO) method implemeted by Wilke [54] has been selected. This global optimization strategy represents an adequate compromise between exploring the design space and locating the optimum. The selected optimization method is an implementation of the optimization method EGO (Efficient Global Optimizer), which has been introduced by Jones et al. [55] and is discussed in Forrester et al. [56].
At the beginning of the optimization a design of experiments (DoE) for a selected number of samples is performed. In this work, the central Voronoi tessellated (CVT) Latin hypercube [57] has been selected as primary DoE technique. For the calculated objective function value and for each selected constraint, a surrogate model based on kriging [58] is built. These surrogate models are able to model the nonlinear behaviour of the objective and constraints. Additionally, a statistical error estimation is included. Based on the surrogate models of the objective function and constraints, a hybrid optimization strategy is used to find the optimum in terms of expected improvement (EI), which combines the predictions of objective function value and model error. The hybrid optimization strategy starts with a global optimization method and the localization of the optimum is improved by the application of a local optimization method. For the global optimization the differential evolutionary (DE) algorithm published by Storn and Price [59] is used. The simplex pattern search method from Nelder and Mead [60] has been selected for the local search in the surrogate models. For the resulting global optimum in terms of expected improvement a recalculation with the physical model is performed. The result of this recalculation is then used to improve the surrogate models for the objective function value and constraints. The described optimization procedure is iterated until convergence is reached.

Active manoeuvre load alleviation
The active manoeuvre load alleviation (MLA) reduces the aerodynamic loads using trailing edge control surface deflections to adapt the lift distribution under manoeuvre flight conditions. In a pull up manoeuvre for example, an inboard load shift can be achieved by increasing the lift in the inboard region with positive control surface deflections and decreasing the lift in the outboard wing region with negative control surface deflections. The result of the inboard load shift are reduced aerodynamic loads in terms of wing bending moment. Thereby, the inboard load shift is related to a forward shifting of the center of pressure in the case of backward swept wings and the aircraft trimming redistributes the lift between the wing and the horizontal tail. As a result, the lift of the horizontal tail increases and simultaneously the lift of the wing decreases.
The aerodynamic limits of the lift distribution adaptation are the minimum and maximum local lift coefficients with deflected control surfaces. These limits can be explained with flow separations, which occur in a viscous fluid with adverse pressure gradient. Furthermore, the geometrical extension of the control surfaces are limited by the rear spar position of the wing box and the required actuator size, mass and actuation power.
In this work, the control surface deflections are modelled using a mesh deformation approach. This approach allows the consideration of viscous flow effects including flow separation in the aerodynamic loads computation, but neglects the complex flow physics around the edges of the deflected control surfaces. The control surface deflection angles are limited to values between −20 • and +20 • due to the robustness of the used volume mesh deformation method.
In Fig. 2 the surface solution for a pull up manoeuvre with control surface deflections for active manoeuvre load alleviation is shown. The corresponding lift and lift coefficient distribution with the center of lift (black dot) are presented in Fig. 3. Furthermore, the elliptical lift distribution (dashed line) with its center of lift (gray square) is shown as reference. With the consideration of active manoeuvre load alleviation further wing box mass reduction and fuel consumption reduction will be expected.
It should be noted that no constraints for stability and control are considered in this work. The presented results show the physical limits of active manoeuvre load alleviation with the prescribed layout of control surfaces and neglect the control surface effectiveness for roll and yaw control. Furthermore, in the used process no structural sizing of the control surfaces is taken into account and no design method of the kinematics and actuators is included.

Structural concept of the more flexible wing
For the more flexible wing, the structural concept and the maximum strain allowable have been changed. The structural concept of the conventional composite wing   [61]. Through a detailed consideration of stringer constraints and stiffness, the evaluation of a more flexible wing becomes possible, while relevant structural constraints are considered. For stiffened composite panels, the more flexible wing concept has been investigated by Bach and Hühne [62].
In this work, the more flexible wing has been modelled with a stringer dominant structural concept of the upper cover. This includes a selected upper skin percentage ply share of (10/80/10) and the usage of I-stringers. Based on the modified structural concept a value of 5000 μm∕m has been selected for the strain allowable of the more flexible wing. The percentage ply share of the lower skin, spars and ribs has been optimized in a preceding aero-structural wing optimization [34] and the values has been transferred to the presented optimizations with more flexible wings. This transfer has been done on the basis of predefined wing box regions as shown in Fig. 4. For each wing box region, the percentage ply share is equal to the prescribed value. In Table 2 the differences between the structural concepts of the conventional composite wing and the more flexible wing have been summarized.

Results
In the DLR project VicToria, multi-mission aero-structural wing optimizations have been successfully applied to optimize wing planform, twist and thickness distribution of the Airbus XRF1 research configuration. To investigate the impact of active manoeuvre load alleviation, the optimizations have been performed for wings with and without the consideration of active manoeuvre load alleviation. Thereby, the XRF1 is an Airbus provided industrial standard multi-disciplinary research test case representing a typical configuration for a long range wide body aircraft.
In the first step, the wing optimizations without consideration of active manoeuvre load alleviation have been performed. The result of the twist and thickness distribution optimization represents the baseline. With the optimization of wing planform, twist and thickness distribution, the design space has been further extended and shows the full potential of multidisciplinary wing optimization with the introduced optimization approach. Additionally, a result for the optimization of wing planform, twist and thickness distribution for a more flexible wing by changing the structural concept and the maximum strain allowable is presented.
In the next step, the active manoeuvre load alleviation has been introduced using active control surface deflections in all manoeuvre load cases. Thereby, the control surface deflection angles have been used as additional design parameters. In these wing optimizations, the wing planform, twist and thickness distribution have been optimized for wings with conventional composite wing structure and for a more flexible wing, respectively.

Design task
The design task describes the objective function, the design space and the constraints. In this work, the wing design for a long range commercial aircraft configuration has been selected.

Objective function, flight missions and load cases
The objective function of the multi-mission aero-structural wing optimizations is the combined fuel consumption of three selected flight missions. In this work, the fuel consumption is defined in terms of fuel burn per range and payload. Hence, the combined fuel consumption is the weighted sum of the corresponding mission fuel consumption as given in Eq. (1).
In Table 3 an overview of the selected flight missions and weighting factors is shown. With the selected weighting factors the expected relative frequency of the missions in operation has been considered. For the study and design mission, the design Mach number of the Airbus XRF1 has been selected. The design mission range is set to 6500nm and the corresponding payload is a result of the aero-structural wing analysis. The selection of range and payload for the study mission is based on a typical long range mission with a passenger load factor of 0.85 and represents the mission for which the aircraft will be optimized primarily. The difference between high speed and the study mission is the increased cruise Mach number to consider off-design conditions in the wing optimization.
For the structural sizing of the wing box, the manoeuvre load cases with the maximum loads have to be defined. These manoeuvre load cases have been derived from the flight envelope limits and the limits of the manoeuvring load factor resulting from the certification regulations CS-25/ FAR 25. In Table 3 an overview of the selected manoeuvre load cases is given.
To compute the fuel consumption of each flight mission, a modelling from conceptual design [63,64] has been used. Thereby, the flight mission has been divided into five segments and the corresponding aircraft mass fractions have been transferred from typical values given in the textbook published by Jenkinson [65] to the Airbus XRF1 reference aircraft configuration. The flight mission segments are summarized in Table 4.
For the cruise segment of the flight mission, a constant Mach number and constant lift-to-drag ratio have been assumed. Furthermore, the thrust-specific fuel consumption has been modelled by a formula published by Mattingly [66]. This formula describes the dependency of the thrust-specific fuel consumption from the flight conditions for a given engine and has been adapted to a typical engine map in the Rolls-Royce Trent 1000 class. The formula with the adopted parameters are given in Eq. (2).
(2) TSFC = C 1 + C 2 Ma g √ SL with C 1 = 0.245h −1 and C 2 = 0.415h −1 .  Landing and taxi m 5 ∕m 4 0.997 [65] The aircraft mass fraction for the cruise segment is calculated with Eq. (3), which has been derived from the wellknown Breguet range equation and the thrust specific fuel consumption of Eq. (2).
For each flight mission, the corresponding lift-to-drag ratio for the cruise segment is a result of the aerodynamic coefficients of the flow simulation for the wing body configuration, the estimated aerodynamic coefficients of the tail, and the given residual drag coefficient as shown in Eq. (4).
The lift coefficient of the horizontal tail is a result of aircraft trimming for the prescribed center of gravity position. This aircraft trimming loop based on the equilibrium of forces and moments around the center of gravity is not described in detail here. For the drag coefficient prediction of the tail a simplified approach from conceptual design [67] based on Prandtl's lifting-line theory and flat plate analogy has been used. The take-off mass of the aircraft is the sum of the residual mass m Res (structural mass without the wing and tail including the operating items mass), the wing mass m W , the tail mass, the payload and the fuel masses (mission and reserve fuel) as shown in Eq. (5).
The wing mass is a result of the structural sizing of the wing box and the tail mass is estimated by scaling the tail mass of the reference aircraft with the tail surface ratio after tail sizing. Thereby, the tail sizing based on conceptual design methods using constant tail volume coefficients [64]. The fuel mass follows directly from the aircraft mass difference for the complete flight mission.
For the calculation of the fuel consumption, the required equations are listed in Table 5. Thereby, the fuel mass ratio m F ∕m TO is computed from the aircraft mass fractions with the given range R and the lift-todrag ratio L/D for the cruise segment. For the study mission and the high speed mission, the payload is specified and the take-off mass has to be calculated. In the design mission the take-off mass equals the maximum take-off mass and the payload is resulting from the cruise flight performance and wing mass. For both cases, the corresponding equations are evaluated in terms of the payload ratio m P ∕m TO . With the fuel mass ratio and the payload ratio the fuel consumption per range and payload follows directly from the last equation in Table 5.

Design parameters and constraints
The outer shape of the wing has been parameterized with the design parameters shown in Fig. 5.
Thereby, the wing planform is calculated from wing area, aspect ratio, leading edge sweep angle and the taper ratios of the inboard, mid wing and outboard wing region. Furthermore, the twist and relative thickness distribution are defined in the corresponding wing sections. In the wing sections between these sections the values of the twist and relative thickness are interpolated linearly. The fuselage shape has been held constant during the wing optimizations. For the belly fairing an adaptation to the root section of the wing has been considered by scaling the middle section of the belly fairing. The positioning of the wing in x-direction has been performed by maintaining the x-position of the aerodynamic center.
The wing box regions are shown in Fig. 4. Thereby, the wing box has been separated into four regions for which the percentage ply share of the lower skin, spars and ribs has been optimized for the more flexible wing concept in a preceding aero-structural wing optimization as mentioned before. The definition of the spars is based on the relative positioning of spar points, which are given in relative span and relative chord coordinates. For the ribs, a constant rib spacing has been considered. In the center wing and inboard region the ribs are oriented in flight direction and in the middle and outboard wing region the ribs are positioned normal to the front spar. These rib orientations are typical for Airbus aircraft. In Table 6 all selected global design parameters are summarized with their corresponding lower and upper bound.
The definition of the control surfaces and the fuel tanks are shown in Fig. 6. For the consideration of the active manoeuvre load alleviation, the deflections of the inboard During the wing optimization, the required fuel tank volume is calculated for all selected flight missions and compared with the useable fuel tank volume. The fuel tank volume constraint has been considered in all wing optimizations.
With the consideration of geometry constraints for the integration of a landing gear and the control surfaces, a better comparability of the optimization results with the baseline aircraft configuration is achieved. Figure 7 gives an overview of the geometrical constraints, which have to be fulfilled for each wing design. This includes the positioning of the main gear wheel on the ground with a given relative x-position while maintaining the minimal allowed distances between the main landing gear, the control surfaces and the wing box.
In Fig. 8 further geometrical constraints for the landing gear integration are shown. For each aircraft category, the outer main gear wheel span has to be within the given limits. Furthermore, the geometrical constraints for nose down engine clearance h NDEC , touch down tail clearance h TC and engine and wing clearances h EC , h WC for a bank angle of ML = 5 • have to be fulfilled. The selected values are given in Table 6.
The introduced landing gear integration consists of a design loop for automatic main gear positioning. Thereby, the main gear wing attachment point is shifted from inboard to outboard position and from front to rear position for each span location. All geometrical constraints are checked for each prescribed position, until a feasible design is found.
In Table 6 the introduced design parameters and constraints are summarized. The design parameters include wing planform and wing section parameters. In addition the inboard rear spar position has been used as design parameter. The control surface deflection angles have been introduced for the wing optimizations with active manoeuvre load alleviation.
The constraints consist of mass constraints, propulsion constraints, geometrical constraints for airport conformity, landing gear and control surface integration constraints, flight mission constraints and structural sizing constraints. In Table 6 the used values and there limits are given.

Wing optimization results without MLA
In this section, the wing optimization results without the consideration of active manoeuvre load alleviation are presented. The selected global design parameters of the wing optimizations are summarized in Table 7. Based on the Airbus XRF1 reference aircraft geometry, the twist and thickness distribution have been optimized with the result representing the baseline. In the next step, the wing planform including the twist and thickness distribution have been optimized. The last optimization additionally considers the structural concept of a more flexible wing by changing the structural concept and the maximum strain allowable.
The resulting wing geometries are presented in Fig. 9 and the corresponding twist and thickness distributions are shown in Fig. 10. Thereby, the wing planform, the wing box geometry with the spars and ribs, the landing gear including the support beam and the control surface geometries are presented. The twist distribution for the rigid "jig-shape" and the elastic "flight-shape" are shown for the optimized wings. The shown "flight-shape" results from the wing deformations at the beginning of the "study  Table 3). For the structural interpretation of the results, the absolute wing thickness distribution is given. In addition, the relative thickness distribution for the aerodynamic interpretation is presented.
The result of the twist and thickness optimization shows a wing with a significant thin inboard wing section. With the aero-structural wing optimization an optimal trade-off between cruise flight performance and wing mass in terms  Fig. 10 and Fig. 16 Relative thickness distribution (t∕c) 1 , (t∕c) 5 , (t∕c) 9 , (t∕c) 14 -See Fig. 10  of combined fuel consumption has been achieved. In Table 8 an overview of all relevant values is given. A value lower than one of the used fuel tank volume for the design mission indicates that the fuel volume constraint is fulfilled. The results of the wing planform optimizations show increased wingspan and reduced taper ratio. This leads to increased cruise flight performance without drawbacks due to wing mass changes. Based on the significant mass reduction of the more flexible wing concept, the wing optimization lead to an optimized wing design with further increased wingspan and further decreased relative aerofoil thickness in the inboard wing region. The corresponding cruise flight performance in terms of lift-to-drag ratios is increased and the wing mass is reduced. In comparison to the wing with conventional wing structure the more flexible wing shows a fuel consumption reduction in the order of 3%.
In Figs. 11 and 12 an overview of the aerodynamic results is given for the optimized wings. For each lift distribution in Fig. 12  This leads in combination with the relative aerofoil thickness to higher values for the isentropic Mach number in the inboard wing region (see Fig. 11). In the pull up manoeuvre, a significant inboard load shift occurs due to static aeroelastic effects of the backward swept wing [69,70]. This inboard load shift is stronger for the more flexible wing in comparison to the optimized wings with conventional wing structure. Thereby, the inboard load shift is related to a forward shifting of the center of pressure in the case of backward swept wings and the aircraft trimming redistributes the lift between the wing and the horizontal tail. As a result, the lift of the horizontal tail increases and simultaneously the lift of the wing decreases.
The resulting wing deformations are presented in Fig. 13 for the cruise flight condition and the 2.5g symmetric pull up manoeuvre in comparison to the rigid jig-shape. The more flexible wing show higher deflections due to the modified structural concept and the increased strain allowable.
In addition an overview of the structural results is presented in Fig. 14. The element thickness (sum of skin thickness and stringer height) distribution is shown for the spars, ribs and covers of the sized wing box.
All optimized wings with the conventional wing box structure show a similar thickness distribution. The upper and lower covers have the highest thickness of all components. Therefore, the covers have the biggest mass contribution. All optimization results show a thickness peak at the engine position and the kink.   In the case of the more flexible wing with the increased maximum strain allowable the mass of the upper and lower covers decreases significantly, because the maximum strain criteria is the dominant criteria for the covers on the conventional wing structure. Furthermore, the spars of the more flexible wing are getting thinner due to the reduced aerofoil thickness and the corresponding reduction of buckling areas, where stability is still the sizing criteria for the spars. While optimizing the aerofoil thickness the thickness of the covers increase slightly with decreasing aerofoil thickness.

Wing optimization results with MLA
In this section, the wing optimization results with the consideration of active manoeuvre load alleviation are presented. The selected global design parameters of the wing optimizations are summarized in Table 9. Based on the Airbus XRF1 reference aircraft geometry, the twist and thickness distribution and the inboard and outboard control surface deflections have been optimized and the result is representing the baseline with active manoeuvre load alleviation. In the next optimization, the wing planform including the twist and thickness distribution have been optimized. The last optimization additionally considers the structural concept   Fig. 16 for the wing optimizations with active manoeuvre load alleviation. As a comparison to the results of the wing optimizations without consideration of active manoeuvre load alleviation, the outer shape and inner wing structure are presented for the wings with active manoeuvre load alleviation.
The results show very thin inboard wing sections for the wing with optimized twist and thickness distribution and the optimized more flexible wing. Based on the significant mass reduction due to the active manoeuvre load alleviation system, the aero-structural wing optimizations lead to wings with decreased relative aerofoil thickness in the inboard wing region. In Table 10 an overview of all relevant values is given. The wing mass of the baseline without manoeuvre load alleviation has been reduced in the order of 2.2t in comparison to the twist and thickness optimized wing with active load alleviation.
The results of the wing planform optimizations show increased wingspan and reduced taper ratio here again. This leads to increased cruise flight performance without drawbacks due to wing mass changes. In comparison to the wings without manoeuvre load alleviation the optimization results with consideration of with active manoeuvre load alleviation show wing geometries with increased wingspan and aspect ratio, decreased taper ratio and improved cruise flight performance in terms of lift-to-drag ratios. The corresponding fuel consumptions have been reduced due to improved cruise flight performance without drawbacks due to operational empty mass changes.
Furthermore, the wing optimizations with active manoeuvre load alleviation result in more inboard center of pressure positions for manoeuvre flight and similar center of pressure positions under cruise flight conditions. This larger inboard load shift is the result of the control surface deflections given in Table 10 and the reduced wing stiffness. All optimized wings with active manoeuvre load alleviation fulfil the fuel volume constraint.
An overview of the aerodynamic results is given in Figs. 17 and 18 for the optimized wings with active manoeuvre load alleviation in comparison to the optimization results of the wings without the consideration of active manoeuvre load alleviation. The resulting wing deformations are presented in Fig. 19 for the cruise flight condition and the 2.5g symmetric pull up manoeuvre in comparison to the rigid jig-shape. The wings with active manoeuvre load alleviation show reduced deflections for the pull up manoeuvre and higher deflections under cruise flight conditions due to reduced wing stiffness. Furthermore, the more flexible wing show higher deflections due to the modified structural concept and the increased strain allowable.
The results of the wing planform optimizations show cruise flight lift distributions with increased inboard loading in comparison to the elliptical lift distribution as shown in Fig. 18. In the pull up manoeuvre, a significant inboard load shift occurs due to active control surface deflections and static aeroelastic effects of the backward swept wing. This inboard load shift of the optimized wings with active manoeuvre load alleviation is larger in comparison to the optimized wings without the active load alleviation system. As expected, the inboard load shift increases further for the more flexible wing. The wing planform optimization of the more flexible wing with active manoeuvre load alleviation shows a strong reduction of wing bending in the pull up manoeuvre. Thereby, the inboard load shift is related to a forward shifting of the center of pressure in the case of backward swept wings and the aircraft trimming redistributes the lift between the wing and the horizontal tail. As a result, the lift of the horizontal tail increases and simultaneously the lift of the wing decreases.  In the wing planform optimization of the more flexible wing, the control surface deflections of the outboard flap and the aileron reach the lower bound of the design space and maximize the inboard load shift. The impact of the corresponding torsional loads is considered in the structural sizing of the wing box structure and results in adequate torsional stiffness of the wing to preserve control surface effectiveness for active load alleviation. Hence, it can be assumed that the aileron effectiveness of the more flexible wing is still acceptable. It should be noted in this context that modern airliners use a combination of inboard and outboard trailing edge devices and spoilers for yaw and roll control with actuations depending on aircraft mass and flight condition. a Values of 3.15 kg CO2 ∕kg Fuel for the fuel burn in a turbofan engine [68] and 105 kg for the passenger mass with baggage are assumed    In analogy to the results without manoeuvre load alleviation, an overview of the structural results for the wings with an active manoeuvre load alleviation system is presented in Fig. 20. With the reduced level of aerodynamic loads the mass of the upper and lower covers decreases significantly. In the case of the wing planform, twist and thickness optimization a similar load level is obtained for an increased wingspan by introducing the active manoeuvre load alleviation system. Furthermore, the spars are getting thinner due to the reduced aerofoil thickness and the corresponding reduction of buckling areas, where stability is still the sizing criteria for the spars. While optimizing the aerofoil thickness the thickness of the covers increase slightly with decreasing aerofoil thickness. By introducing the more flexible wing concept, the increased maximum strain allowable leads again to reduced thickness of the wing box structure.

Result overview and assessment of active manoeuvre load alleviation
The main results of the aero-structural wing optimizations are summarized in Table 11. The percentage of reduction in combined fuel consumption is presented in the last row and shows the potential of active manoeuvre load alleviation to reduce the CO 2 emissions per passenger kilometre. With the introduced active manoeuvre load alleviation a reduction of the combined fuel consumption between 2% and 3% has been achieved for the optimizations with conventional composite wing box structure. The active manoeuvre load alleviation shows for wing optimizations with the more flexible wing concept an additional reduction of the combined fuel consumption in the order of 2% due to snowball effects. In comparison to the twist optimized reference aircraft Airbus XRF1, the introduction of the active manoeuvre load alleviation technology and of the more flexible wing concept, leads to reductions of the combined fuel consumption in the order of 13% after aero-structural wing optimization. The results for the reference aircraft Airbus XRF1 are not presented here, because of restricted publication rules.
In Fig. 21 an overview for the results of the aero-structural wing optimizations is shown. For all aero-structural wing analysis the cruise flight performance, the wing mass ratio and the corresponding combined fuel consumption are summarized. Due to the reusability of many results for constructing the surrogate models of different wing optimizations, the overall sum of aero-structural wing analysis in this work is in the order of 750.
The combined fuel consumption depends on cruise flight performance and wing mass. The global search in the design space shows different combinations of these two counterparts with the same objective function value. However, the global search has been converged with a acceptable accuracy for practical wing design and optimization.

Conclusion and Outlook
In this work, the assessment of the active manoeuvre load alleviation has been successfully demonstrated using an integrated process for aero-structural wing optimization based on high fidelity simulation methods. The comparison of optimization results with the same objective function, design parameters and constraints allows a proper technology assessment. To find the optimum trade-off between aerodynamic performance and wing mass, the twist and thickness distribution and the wing planform design parameters have been involved in the wing optimization. The results of this optimization show the expected reduction of the combined fuel consumption due to decreased wing mass based on active adaptation of the load distribution in manoeuvre flight and an increased aerodynamic performance under cruise flight conditions. This increase in aerodynamic performance  has been achieved with higher aspect ratio and reduced taper ratio of the wing. With the application of active manoeuvre load alleviation, the significance of the aerodynamic limits of control surface deflections due to flow separations for the prediction of manoeuvre loads has been shown. The structural wing sizing based on a reduced level of aerodynamic   Another aspect of the wing optimization results obtained, is the resulting thin inboard wing, which contradicts conceptual design estimations based on handbook methods and is significantly lower than the values of current aircraft designs. One reason for this is the requested flexibility of airlines for extended range missions, which requires additional fuel tank volume.
In the future, the processes for wing optimization have to be extended with engine design and integration, gust loads prediction and take-off and landing constraints. The assumptions for the usage of linear structural analysis have to be verified for the increased deformations of more flexible wings and compared with non-linear structural analysis results. The corresponding impact on the optimization results and the computing times have to be investigated. Furthermore, the aerofoil design has to be integrated with the usage of an aerofoil catalogue or by gradient-based shape optimization.
In addition a design method of the actuator system for the active manoeuvre load alleviation functions should be implemented. The structural analysis and sizing process have to be extended for the sizing of the control surfaces itself and for a more local load introduction of the control surface loads to the wing box structure.
With increasing computational resources and progress in numerical processes, based in high fidelity methods, it