Investigations of passive wing technologies for load reduction

Wings of modern aircraft have to be designed to give optimal response with respect to loads, comfort and performance. An essential part of the wing development is thus a design process which can take all these aspects into consideration. In the “Adaptive Wing” work package of the CleanSky “Smart Fixed Wing Aircraft” project, a multi-fidelity wing design method using aeroelastic tailoring has been developed. In the article, the process is presented in detail. The approach is based on a parametric wing design approach. Both beam models and shell models are derived and optimized in separate optimization environments. Investigations of the use of unbalanced laminates in aeroelastic tailoring are presented, employing the optimization of lamination parameters. The applications are demonstrated on two aircraft configurations, a long range and a short range transport aircraft. Further developments presented in the article include the introduction of CFD-based aerodynamics in the tailoring process, and a process extension to assess the influence of aeroelastic tailoring on fatigue.


Background
In the framework of the CleanSky "Smart Fixed Wing Aircraft" (SFWA) project, one field of activity has been the investigation of technologies for loads reduction for transport aircraft. Aim of the work package "Adaptive Wing" in SFWA was to create design processes and solutions for aircraft wings giving optimal response with respect to loads, comfort and performance. This article presents the work in the field of aeroelastic tailoring performed in the project by the Institute of Aeroelasticity of the German Aerospace Center (DLR), the Aerospace Structures and Materials Group of the Delft University of Technology, and the Fraunhofer Institute for Structural Durability and System Reliability.
The structure of the article mirrors the approach taken in the project. In the first phase, the focus was on the development of key technologies of all partners and the establishment and improvement of cooperation among the partners and in the work package. In this phase, the basic multi-fidelity aeroelastic tailoring process was developed, making use of DLR's parametric modelling process, see Sect. 2. In the second phase of the project, existing design technologies were matured on a common numerical aircraft platform, the so-called XRF1 model, a long range aircraft configuration provided by Airbus [1,2], see Sect. 3. The focus of the investigations was on the assessment of the use of unconventional laminates for aeroelastic tailoring and its influence of the resulting wing mass. In the third phase, the developed technologies were applied on a second aircraft, the so-called "new short range" (NSR) concept, see Sect. 4, which was used for extending the aeroelastic tailoring process by the introduction of CFD-based aerodynamics into the process. Finally, an assessment of the influence of This paper will be part of a Special Issue on the EU project Smart Fixed Wing Aircraft.

3
passive load alleviation schemes on fatigue-relevant loads was performed, see Sect. 5.
Various forms of load control have been employed in aircraft over the years. Regan [3] lists a number of examples for active manoeuvre and gust control on series aircraft; however, for commercial aircraft it is difficult to obtain technical information concerning the exact implementation of the systems. At DLR, research on active load control was performed on the ATTAS aircraft since the 1990s [4]. A wellknown set of investigations is also the Active Aeroelastic Wing (AAW) flight research program, a joint US Air Force/ NASA program [5]. Recent studies include a gust load alleviation approach using a folding wingtip device [6] and aeroservoelastic investigations on the so-called SUGAR Truss Braced Wing configuration [7]. Livne gives an extensive literature overview on load control in his recent article [8].
The optimization of composite materials used to form the load carrying structure of a wing is a research topic that has been discussed for decades; see for example Starnes Jr et al. [9], where the authors apply an approximation-based optimization to minimize the weight of a wing structure subjected to different combinations of buckling, strength, displacement and twist constraints. Hollowell et al. [10] demonstrate the considerable effects induced by the coupling terms in the bending stiffness matrix by analysing the effect of bending-twist coupling on the aeroelastic behaviour of composites' plates with various symmetric, unbalanced stacking sequences. A valuable survey on aeroelastic tailoring effects as a result of the directional stiffness in orthotropic composites is provided in Shirk et al. [11]. The authors investigate the possible influence of rotated fibre angles on minimum weight, twist, normal modes, flutter and aerodynamic performance for various configurations, reaching from fighter aircraft to forward swept configurations. A general overview on optimization technologies reaching from optimizations on panel level to the aeroelastic optimization of composites in aircraft wings is provided in Vanderplaats et al. [12]. In a more recent work, Leon et al. [13] maximize the flutter eigenfrequency of a composite plate wing with ply angles as design variables, considering symmetric and non-symmetric stackings.
Aeroelastic tailoring studies based on lamination parameters have been attempted previously, though not to the same extent as the previously depicted optimizations comprising discrete stacking sequences. Kameyama et al. [14] use a composite plate wing to demonstrate the influence of lamination parameters on the flutter and divergence characteristics. Minimization of compliance of a variable stiffness slender wing that is represented as a beam is demonstrated in Abdalla et al. [15]. The cross-sectional properties of the beam are parametrized using lamination parameters that define the membrane stiffnesses of the box cross section. Results show that variable stiffness layups can significantly outperform structures comprising constant stiffness. Unbalanced laminates in particular, comprising bending-torsion coupling, lead to the best performance in terms of compliance. A two-level optimization strategy has been proposed by Liu [16]. Allowing for symmetric and unbalanced laminates, a lamination parameter-based weight minimization is performed before applying genetic algorithms for the derivation of stacking sequences. The objective function in the second optimization step comprises a minimization of the square distance of lamination parameters. Dillinger [17] presents an optimization approach for composite aircraft wings based on lamination parameters. The introduction of that publication also gives an extensive overview over the development in the field. Next to those sources mentioned, the already cited article by Livne [5] also contains a section on aeroelastic tailoring.
Both the processes developed on the work, i.e. the multifidelity aeroelastic tailoring process, and the investigations performed, i.e. the assessment of the use of unconventional laminates for aeroelastic tailoring, are considered very significant contributions to the field of aeroelastic tailoring, especially as they have been demonstrated not only on academic applications, but on realistic aircraft configurations based on industrial models. Furthermore, a technology readiness level (TRL [18]) of three and higher could be demonstrated for the wing design process by DLR and TU Delft.
2 Multi-fidelity aeroelastic tailoring process for passive load alleviation

Approach
Passive load alleviation can be achieved by using the directional stiffness properties of composite materials to tailor the aeroelastic response of the wing. This results in a large design freedom for the designer, making it a challenge to explore the aeroelastic design space efficiently. Therefore, a multi-fidelity multidisciplinary approach to the optimization of a composite wing structure has been developed in the "Adaptive Wing" work package, which employs a two-step optimization procedure. The first step is the optimization of the wing structure using a low-fidelity nonlinear aeroelastic beam model. The result of this optimization is then used for a more detailed optimization of the wing structure using a shell model coupled to doublet-lattice (DLM) aerodynamics implemented in NASTRAN. A comparison between the beam model and the shell model clearly shows the validity of this approach, thus making it suitable for the optimization of aeroelastically tailored wing box structures. The work was performed in cooperation between the DLR Institute of Aeroelasticity and Delft University of Technology [19,20].
The underlying analysis model is defined on the basis of a parameterized geometry. That way, a convenient process for the investigation of parameter spaces can be set up. The parametric structural design tool employed in the process is ModGen [21], developed at DLR. The tool is capable of generating a parametrized structural representation on the basis of the geometry, which is then processed to a NASTRAN sizing optimization.
Two types of models can be derived. First, the initial ModGen shell model is reduced to a computationally less expensive beam model, on which a pre-sizing optimization is called, see the upper branch of Fig. 1. The optimized skin and spar stiffness is calculated in the form of ABD-matrices, representing stacking sequences. Second, the initial ModGen shell model is transferred into a finite element model, on which a structural optimization with aeroelastic constraints is performed, see the lower branch of Fig. 1. Here, too, the stiffness distribution is optimized in the form of ABD-matrices.
Thus, the designer has several choices, depending on the complexity of the models or the design and analysis task-the multi-fidelity process provides representative beam models for wing structures which can be used for the allocation of meaningful start values for the optimization of complex shell models. However, in many applications, aeroelastic analyses on high aspect ratio wings can be performed with good accuracy on beam models directly. For more complex geometries, wing structures can be designed and optimized based on the more detailed shell models.
The following sections will explain the process in more detail and give application examples.

Parametric wing definition
For the generation of the structural wing model, the tool ModGen was used. ModGen has been specifically developed for the generation of parametric representations of structural and aeroelastic models to be used in aeroelastic analysis and optimization tasks; see for example the application on the so-called FERMAT, a structural model for the CRM configuration [22], or on a forward-swept wing aircraft [23].
Based on a definition of the wing surface, ModGen is capable of generating a finite element shell model within this contour, involving wing skins, spars, ribs and stringers, the latter one being represented by beam elements. The wing surface can either be provided as surfaces in the well-known Initial Graphics Exchange Specification (IGES) format or by sets of airfoil contours and their positions in space. ModGen also provides a DLM model for aerodynamics and a splining model between the structure and DLM grid; see Fig. 2.
Furthermore, mass models considering fuel distributions can be generated. Finally, ModGen provides the required input data for the sizing as an optimization problem, i.e. the design and optimization variables. Optimization can be performed in NASTRAN [24], using SOL 200 for optimization based on the native NASTRAN elements. However, in the aeroelastic tailoring process developed in SFWA and presented here, structural optimization is performed by an optimizer provided by TU Delft.

Wing design optimization using fast design methods
TU Delft has focussed on the passive wing optimization on a beam level. The work was carried out in close collaboration with the DLR Institute of Aeroelasticity. The passive The intricate problem of composite wing skin and spar optimization, which inherently is associated with a large number of design variables and constraints, was handled by setting up a multi-fidelity design loop, making up the "upper" branch in Fig. 1. TU Delft developed a low-fidelity aeroelastic optimization routine which can quickly design a spanwise and chordwise stacking sequence distribution; see [25]. This low-fidelity design can then serve as an initial guess for the DLR medium-fidelity optimization.
The TU Delft low-fidelity aeroelastic optimization strategy consists of the monolithic coupling between a continuous-time unsteady vortex lattice aerodynamic code and a geometrically nonlinear Timoshenko beam code. The structural dynamics are linearized around the static nonlinearly deflected shape. The composite laminates can be varied in chordwise and spanwise direction over the wing in various design zones. The discrete composite stacking sequence discretization is parameterized in a continuous fashion using lamination parameters. This way, gradient-based optimization is realizable, which is an essential ingredient when dealing with an optimization problem of multiple hundreds or thousands of design variables. Using this approach, the wing can be sized for a variety of manoeuvre load cases including constraints such as buckling, strength, divergence, aileron efficiency and flutter. More details of the approach are given in [20].
The methodology was applied to the new-short range (NSR) aircraft, see Sect. 4, and it was shown that anisotropic laminates can reduce the wing structural weight by 30% as compared to quasi-isotropic composite laminates. A typical result skin thickness for the wing is shown in Fig. 3, left, together with the optimized stiffness distribution, Fig. 3, right. Here, the dark curves indicate the main stiffness direction in a design field, a slender oval representing a unidirectional local stiffness, a rosetta or a square representing a more isotropic local behaviour.

Aeroelastic design and tailoring process for dynamic global finite element models
For the application of the aeroelastic design process based on the shell model, according to the lower branch of Fig. 1, it is assumed that a basic aircraft design exists. Figure 4 provides an overview of the design loop. The aircraft configuration and global aircraft data form an input to the procedure. For the description of the wing, the basic wing layout and the wing topology should be available. A given basic structural design, material properties, loads envelope, formulation of constraints, non-structural masses and fuel can be included if available, otherwise suitable data are created as part of the wing design process (Fig. 4, "input" step). Based on the input, a wing model is created, with a focus on the parametric structural model (Fig. 4, "model generation" step). The model consists of the structure ("FEM"), a DLM model ("Aero"), and the definition of the design variables ("DVar") for the subsequent optimization step. A preliminary sizing of the wing structure, to get a good starting point for the optimization, is part of the model generation. In the next step ( Fig. 4, "sizing and optimization" step), the optimization problem is formulated, and the sizing of the wing structure is performed using optimization of the stiffness distribution according to the requirements; such requirements are usually minimum weight, a prescribed flight shape, minimum static and dynamic loads, as well as required control surface efficiency. The output of the design process can be a separate wing analysis model; however, models for all the components, main wing, empennage and fuselage, can be generated. In that case, a full aircraft model will be assembled from the components ("analysis model"), which is then available for further design studies (Fig. 4, "analysis" step), for example the design of control laws, aeroelastic analysis or more detailed load analyses, e.g. with CFD-based calculations. The results of the latter can be included in the sizing loop, if desired. A wing optimization in the presence of aeroelastic constraints is often called "aeroelastic tailoring", usually associated with the use of composites, i.e. fibre-reinforced materials. This process is the core of the "sizing and optimization" block of the aeroelastic design loop shown in Fig. 4. The aeroelastic tailoring process is depicted in more detail in Fig. 5. The parameterized wing model which forms the input of the tailoring step is taken from the modelling step "model generation" (Fig. 4), i.e. a NASTRAN model including wing structure and mass model ("FEM"), aerodynamic ("DLM") grid and design and optimization variables ("DVar"). The optimization problem is defined in MATLAB [26] and the tailoring loop is started. A number of load cases relevant to the optimization are defined in MSC.NASTRAN. This tool is used both for the load analysis of those cases and for the subsequent generation of the sensitivities of the objective function and responses with respect to the structural design parameters, e.g. stiffness distribution and layout of laminates. Those sensitivities ("SENS") are converted to linear and reciprocal sensitivities with respect to the stiffness matrices. The approximation model ("gradient-based approximation") replaces the actual analysis model in the search of the optimizer for a minimum of the objective function, greatly accelerating the function evaluations required during the optimization. The step "structural optimization" consists of a gradient-based optimization in the TU Delft optimizer ALDO [19], which generates feasible stiffness distributions for the wing structure, described in the form of the so-called ABD-matrices, the stiffness matrix in classical laminate theory [19], along with consideration of their feasible region. Each optimization step results in a new set of design variables ("VAR") that represent the global optimum of the convex approximation subproblem. Convergence is monitored by determining the change in the objective function for subsequent feasible iterations.
The optimization implemented in the process is capable of considering balanced and unbalanced laminates. For all subsequent analyses using the resulting model ("optimized model") described in the following sections of the article, the general formulation (stiffness, mass distribution, "ABD opt ") was used for the FEM. However, an optimized laminate stacking sequence which represents the optimal ABD-solution ("ABD opt ") can be generated in a post-processing step ("optimization of stacking sequ."). A more detailed description of the optimization process can be found in [27]. The process has been used for wings of several aircraft configurations. The applications will be described in the following sections.

Process validation: the ICW wing
The aeroelastic tailoring process was built up and verified using the academic example intermediate complexity wing, ICW, for tools and process development. The ICW is a reference composite wing, which has been examined in several optimization papers [28]. ModGen has been used to generate a finite element model with 32 upper and 32 lower skin shell elements with PSHELL/MAT2 property definition. Design variables are lamination parameters and laminate thicknesses. The design variable sensitivities for specified responses are computed with NASTRAN SOL200 and the sensitivities passed to the optimizer. Stacking sequences of the wing box are extracted, and the stiffness matrices of the box plates calculated. Figure 6, left, shows the layup for an optimization with respect to a required control layout efficiency. The lamination parameters indicate the direction of the main stiffness, V1 and V3 correspond to the A matrix (membrane stiffness) and the parameters W1 and W3 to the D matrix (bending stiffness). The coupling matrix B is zero because a symmetric laminate is chosen. The parameters V2, V4, W2, W4 are also zero because a balanced laminate was regarded. The resulting global stiffness properties were compared to literature results, and a good accordance is found.

Use of balanced and unbalanced laminates on the XRF1 with aeroelastic constraints
The first transport aircraft application of the aeroelastic design and tailoring process in the SFWA project was the so-called XRF1, a model of a widebody configuration, which has been distributed by Airbus in the framework of SFWA. The goal was to provide an aircraft model to research institutes which includes structure and aerodynamics data as a base to allow investigations with a model of relevant complexity and behaviour. However, the industrial model turned out to be too complex for a quick variation of stiffness properties needed for the work in the Adaptive Wing work package. It was thus decided to generate a similar, but not identical, parametric aircraft model, see Fig. 7, using the aeroelastic design process described in Sect. 2, and based on the operational parameters provided by Airbus with the XRF1. Part of the model generation was the definition of a representative flight envelope, e.g. cruise and manoeuvre speeds, and altitude. Several investigations were performed on the XRF1 model. First, a statically and dynamically equivalent shell model of the XRF1 wing was built in the parametric model generator, suitable for subsequent structural adaptations and modifications. The topology followed that of the original XEF1 wing; see Fig. 8. A finite element model with a defined set of thickness and layup parameters was generated, where the parameters were derived from the original XRF1 wing. A static validation using the application of forces and a comparison of wing deformation between the original XRF1 and the DLR variant was performed, as well as a dynamic validation including non-structural masses which were derived from the XRF1 wing. Finally, a structural optimization of the DLR wing with respect to aeroelastic objectives and constraints, especially weight and control surface efficiency, was performed. The basic DLR model was set up assuming balanced laminates.
In a second step, a study comparing the influence of using balanced and unbalanced laminates on wing mass in the presence of aeroelastic constrains was performed. Aeroelastic constraints include a required aileron efficiency as well as constraints on maximum local wing twist. Two wing designs are compared-a first model consists of standard balanced laminates, a second, otherwise identical, model consists of unbalanced laminates. Both models are optimized for minimum weight under the requirement of varying aileron efficiency and a maximum acceptable wing twist. The laminates remain symmetric, i.e. the matrix B is zero. Figure 9 shows the discretization of the wing for optimization. The wing is partitioned in design fields spanwise and chordwise, 34 fields per skin, i.e. 68 fields total. 13 NAS-TRAN design variables are defined per field, resulting in a total of 884 design variables. 486 elements are used per skin, i.e. a total of 972. The responses taken from the analysis are the weight response (i.e. mass), stress responses from the design fields (based on which strain and buckling failure are computed externally), twist responses of the wing and aileron efficiency response. The objective function is Seven load cases taken from the assumed operational envelope were selected as wing sizing load cases for the current investigation, i.e. two symmetric pull-up manoeuvres at n z = 2.5, a symmetric pushdown manoeuvre at n z = − 1, three rolling cases and one trim case at cruise speed. Figure 10 shows two typical results for the thickness distribution of the wing, derived from the optimizations result not including the aileron efficiency constraint, Fig. 10, left, and including the constraint, Fig. 10, right. It can be seen well that the optimizer increases wing thickness outboard to reduce wing twist under aileron deflection. Figure 11 shows a typical stiffness distribution of an optimized wing with unbalanced laminates. Similar to Fig. 3, right, the red curves in Fig. 11 indicate the main stiffness directions in a design field (see [17] for a more in-depth explanation of the visualization of distributed laminate stiffness). It can be seen how the main stiffness direction points backwards on the wing tip, indicating that a pronounced bending torsion coupling, opposing the geometric wash-out effect, is induced by the optimizer to increase the aileron effectiveness and thus fulfil its constraint, noting, however, that this also leads to an increase in inboard loading for highg manoeuvres or gust loads.
To verify whether a sufficiently global optimum of the structural optimization has been reached, the optimization was started from various initial conditions. It could be shown that optimization results are independent from the Significant differences can be seen for the use of balanced and unbalanced laminates under consideration of aeroelastic constraints, here especially the high speed aileron efficiency. By using balanced laminates, the aeroelastic constraints are mainly met by an increase of torsional stiffness in the outer wing, resulting effectively in a mass increase because of increasing material thickness. In case of unbalanced laminates, the bending torsion coupling on the outer wing becomes more pronounced, the material thickness remains thinner and the wing mass remains lower, overcompensating the necessary additional material needed for the increase in torsional stiffness on the outer wing.
For a wing with balanced laminates, the mass increase due to the introduction of the aeroelastic constraints is about 50%. This trend is generally well known, of course, and a common solution is the use of split ailerons on production aircraft to secure the required aileron efficiency at high speeds. However, it is interesting to see that the mass penalty for a wing with unbalanced laminates is 10% smaller than the penalty of the design with balanced laminates.
As a consequence of the work, it can be shown that aeroelastic constraints are essential for adequate and meaningful optimization of wing structures. A requirement on aileron efficiency in particular greatly affects the structural weight; the effect could be quantified for different requirements. Second, the investigations prove that the introduction of unbalanced laminates has a great potential for weight reduction and loads reduction. As the XRF1 model is based on an Airbus design, numerical values for the optimization results of the XRF1 application example used in the project cannot be given in the article. However, an identical approach was used for a generic forward-swept wing design, resulting in similar observations, where the results are published in [29].

CFD-based corrections for aeroelastic tailoring of the new short range aircraft configuration
The so-called new short range (NSR) aircraft is a configuration developed in the SFWA project by Airbus for global design studies. For the use in the Adaptive Wing work package, Airbus supplied a geometric database of a half-model consisting of wing and fuselage (Fig. 12) [30].
The half-span of the model is approximately 18 m, and the fuselage length is approximately 40 m. In the Adaptive Wing work package, a representative structure for the design is developed.
Using the aeroelastic design process, an initial wing design is defined which is then sized using the aeroelastic tailoring process described above. Figure 13 shows the layout of the wing box of the NSR wing. The engine is assumed to be located under the wing similar to the current A320 layout, rigidly attached to the wing structure as a point mass. The leading edge and trailing edge components of the wing as well as the systems and the landing gear are considered as mass points.
The wing is sized using a number of selected static load cases. The load analysis as well as the generation of sensitivities necessary for optimization is performed using NASTRAN. In Fig. 14, the colours indicate the design field distribution of the NSR wing. Properties of neighbouring shell elements (sharing the same colour) were combined as design variables to reduce the total number of variables changed by the optimizer.
In the studies following the design, the goal of the analysis was the investigation of the influence of aileron effectiveness constraint on minimized wing skin mass. Calculations were performed for a set of required aileron effectiveness constraints of η ail ≥ 0.00 and η ail ≥ 0.03. The aileron efficiency is expressed in terms of the so-called "helix angle", the angle between flight path velocity v ∞ and the wing tip velocity p•s, see Fig. 15. The wing rotates with a roll rate of p (taken from the NASTRAN analysis), and s is the semispan. The aileron effectiveness is equal to the arc tangent of the helix angle for unit aileron deflection.
The load cases which were regarded included manoeuvre load cases of n z = − 1 and 2.5 at a Mach number of 0.597 and altitudes of 0 m and 6700 m. Load cases relevant to aileron effectiveness were defined at several Mach numbers and altitudes, as well as for several mass cases.
The optimization model contained parameters representing the structural responses (e.g. mass, strain failure, buckling failure), as well as the aeroelastic responses, i.e. The approach includes the consideration of an aero load correction to improve the aerodynamic quality of the NASTRAN internal DLM by means of a higher order CFD method; see Figs. 16 and 17 below. For the CFD calculations, the DLR TAU code was used, see [31]. For the current study, only CFD calculations for symmetric manoeuvres were included, i.e. the calculation of aileron effectiveness according to Fig. 15 was determined with DLM only.   Details on the modelling and the analysis results for a forward swept, but in terms of dimensions similar wing, can be found in [32].
As a sample result it was found that increasing the aileron effectiveness constraint from η ail ≥ 0.00 (no reversal) to η ail ≥ 0.03 results in a weight increase of wing skins and spars of approximately 11%.

Assessment criteria
The present investigations focus on a methodology for the evaluation of load control concepts with respect to their potential for a reduction of structural weight. Fatigue-relevant loads are sizing for a number of structural components. A decrease in static loads or low-frequency loads should not lead to higher fatigue-relevant loads on those components. Even if the implementation of load control concepts leads to less severe loading of structural components whose mass should or cannot be reduced, the benefit of the load control concept could alternatively be formulated, e.g. in terms of longer maintenance intervals and a resulting reduction in operating costs. The effect of active or passive measures to reduce dynamic loads on fatigue has thus to be taken into consideration.

Assessment methodology
The purpose of load control concepts is to change the interaction of the aircraft with the airflow. It is therefore liable to change many aspects of aircraft loading from the overall flight mechanics to the distribution of loads or stresses among structural members. This holds true for active as well as passive load control. To assess the overall potential for weight reduction available through the implementation of load control concepts it is therefore necessary to use computational methods modelling these mutually interactive effects. Furthermore, the level of detail of the analyses needs to be adaptable to the stage of aircraft design they are used for.
Methodologies with these features have been available in SFWA, so the wing and aircraft design process described in the sections above, based on ModGen and NASTRAN. A finite element (FE) model of the aircraft is subjected to modal reduction and transferred to the MATLAB Simulink environment. Additional submodels for propulsion, gravity, airflow, etc. are added so various operating conditions and loads can be simulated based on the reduced aircraft model. Among the results are load-time histories at certain predefined monitoring stations. In preliminary development stages, submodels may be represented by lower order formulations than in later high-fidelity computations. Thus, computational cost can be offset against the accuracy of results depending on the level of detail of the current design phase. Also, additional submodels can be added in a modular fashion.
In a first step towards an evaluation of load control, baseline analyses are conducted for the relevant load cases without load control. The resulting load-time histories at the monitoring stations are then subjected to a fatigue evaluation consisting of rainflow-counting, mean-value transformation, and damage accumulation, see Fig. 18. In early design phases, generic stress-cycle (so-called SN-or Wöhler-) curves could be used to obtain a virtual damage sum.
In a second step, the load control concept is incorporated into the aircraft model. In case of passive load control (aeroelastic tailoring), the FE-model would be modified accordingly before being subjected to modal reduction. Aerodynamic simulations are then used to modify the airflow submodel as well. The analysis then proceeds as before. Active load control systems would result in an additional submodel including their actuators and control laws. Additional aerodynamic simulations can be used to incorporate the aerodynamic effect of the load control system into the airflow submodel or a separate submodel. In any case, the load-time histories at the monitoring stations will be obtained for the relevant load cases as before and then subjected to the fatigue evaluation mentioned above.
A comparison between the virtual damage sums with and without load control will then be used to evaluate the effect of the load control concept. Depending on the level of detail, this can be expressed in a comparison of preliminary weights of structural members, if possible. The approach presented here is documented in [33]. The work performed on the NSR aircraft presented in the following section was performed on the basis of the suggested approach.

Investigation of the influence of aeroelastic tailoring on fatigue loads
An investigation of the influence of aeroelastic tailoring on fatigue was performed as a cooperation study between the DLR Institute of Aeroelasticity and the Fraunhofer Institute for Structural Durability and System Reliability. Analysis of fatigue flight loads is usually based on turbulence spectra for trimmed flight conditions. The assumption is that continuous turbulence occurs as a spatially varying stochastic process. The spectra are defined as power spectral density (PSD) functions. Well-known examples for turbulence spectra are the von Kármán wind turbulence model and the Dryden wind turbulence model. When a wing structure is optimized with a heavy focus on load reduction, its natural frequencies and mode shapes might differ considerably from a standard wing design. The same is true for the use of unbalanced laminates in comparison to the use of balanced laminates. The frequency response of the wing when excited by atmospheric turbulence will thus be affected. Consequently, an activity in the SFWA work package Adaptive Wing was an investigation whether a passive wing design, tailored for minimum loads, would influence the fatigue behaviour of the resulting wing with respect to a standard layout.
In this study the new short range (NSR) configuration described in Sect. 4 is used. The two wing designs presented above, with aileron effectiveness constraints of η ail ≥ 0.00 ("wing 1") and η ail ≥ 0.03 ("wing 2"), were compared. On the wing, so-called "monitoring stations", i.e. nodes on the wing surface, are defined for which the analysis is performed. First, the transfer functions for the gust input on the response of the two wings at the monitoring station, here a vertical load, are calculated in NASTRAN. The transfer functions are exported, the further process steps are performed in a MATLAB environment. Second, the continuous turbulence input is defined as a PSD function. For the study, the input is a standard turbulence spectrum at moderate excitation for cruise flight conditions. Next, for the given transfer functions, PSD load analyses using the given turbulence spectrum are performed. The results are PSD representations of the wing response at the selected monitoring points, see for example Fig. 19, for the vertical load at the wing root of the wing with η ail ≥ 0.03, i.e. wing 2. In the plot, the dotted red lines indicate the natural frequencies of the wing.
Fatigue analyses usually work based on the evaluation of time series. Thus, an equivalent time history for the PSD Fig. 19 Results for stochastic gust analysis: PSD of the vertical load on the wing root; red lines: wing natural frequencies results was generated for each monitoring point; see Fig. 20 for loads at the wing root.
Finally, exceedance curves for load levels were created. Here, discretized load levels are plotted over the number of crossings of those load levels (i.e. exceedance). Those plots will then be the basis for a fatigue evaluation of the wing concepts. The results for two wing designs are presented in Fig. 21. For the given excitation, the wings show very similar behaviour for higher-amplitude responses with low occurrence. For a higher number of exceedance, wing  2 shows higher responses. This might be due to the fact that the wing has a higher torsional stiffness to maintain the required aileron efficiency. Summarizing the results of the study, the process showed its capability for the required investigations, however, for the given wing designs, no great difference in fatigue behaviour would probably be expected.
It could thus be shown that the aeroelastic design and tailoring process is capable of generating the necessary input for a subsequent fatigue analysis which can than evaluate whether a wing design optimized for load reduction suffers potential disadvantages when regarding fatigue. However, such a final evaluation is heavily based on the selected material, and thus outside the scope of the activities performed in the Adaptive Wing work package.

Summary and outlook
In the CleanSky "Smart Fixed Wing Aircraft" project, a multi-fidelity wing design process based on parametric modelling and aeroelastic tailoring has been developed. The process could be demonstrated on a long range and a short range aircraft configuration, with investigations on the use of unbalanced laminates in structural optimization, the influence of the aileron efficiency on wing mass, the introduction of CFD-based aerodynamics in the process, and an important step towards the assessment of the influence of passive load control on fatigue. The wide field of work could only be covered because of the fruitful cooperation of the partners whose expertise was complementary.
Next steps to be taken are the integration of the aeroelastic tailoring process in a comprehensive loads loop, in order to increase the number of relevant load cases in the optimization. This effort was started by DLR [1,34]. Furthermore, the current optimization is based is load cases without the application of active load control. Thus, a combined optimization of structure and control parameters (sometimes termed "aeroservoelastic tailoring") is an important further development. Finally, the investigations of the relevance of aeroelastic tailoring on fatigue will be pursued by the partners.