Abstract
1. Introduction
Despite tremendous technological headway in commercial aviation over the past 60 years, the increase in the amount of air traffic requires further reductions in engine and noise emissions to render commercial aviation sustainable in the long-term perspective. In response, Advisory Council for Aeronautics Research in Europe (ACARE) put forth a set of ambitious goals: a 50% reduction in carbon dioxide, an 80% reduction in nitrogen oxides and a 50% reduction in perceived noise relative to the emission and noise levels reported in the year 2000 (ACARE, 2010). A similar initiative was also started in the United States by NASA’s Environmentally Responsible Aviation (ERA) project. As a result, many research projects investigating, among others, novel airframe configurations, propulsion and laminar flow technologies were funded by both the EU and the United States (Graham et al., 2014; Suder, 2012). According to Barbarino et al. (2011) and Weisshaar (2013), airfoil morphing can contribute to aircraft performance in several ways. It can improve laminar and basic turbulent flow conditions. It can be utilised for continuous lift-to-drag ratio optimisation or adaptation to specific flight phase requirements such as take-off, cruise and landing. It can also provide for more efficient flight control mechanisms. Hence, it is no surprise that morphing wings, airfoils and high-lift devices were the focus of research in several research projects, such as, among others, SADE, SARISTU and NOVEMOR (Kintscher et al., 2011, 2014; Vasista et al., 2016).
In particular, morphing high-lift devices such as the morphing leading edge investigated in this article bear great potential for improving the aerodynamic efficiency of future aircraft as well as reducing their acoustic intensity by allowing for flow laminarisation (Morgan, 1986) and by eliminating major noise sources associated with conventional high-lift devices (Casalino et al., 2008; Dobrzynski et al., 2001). These technologies could, therefore, prove essential in achieving ambitious environmental goals posed by ACARE and NASA.
A lot of research effort has already been invested into various aspects concerning morphing droop nose. As reported by Monner et al. (2009), the ‘active rib’ concept initially proposed by Dornier (Patent DE2907912-A1, Dornier company, 1979) was investigated and further developed within the SADE project. A dedicated integrated workflow was established to first optimise the droop nose shape in combination with the flap settings, derive the morphing skin properties and finally develop the kinematic mechanism. Kintscher et al. (2011) used the developed workflow to design a demonstrator for the wind tunnel test envisioned within the project. The demonstrator was ground tested, and Monner et al. (2012) reported good overall agreement with the numerical model concerning morphing accuracy. Kirn and Storm (2014) further improved the ‘active rib’ concept by reducing its mass and kinematic complexity within the SARISTU project.
Opposed to conventional morphing mechanisms investigated in SADE, two morphing droop nose concepts based on compliant mechanism were investigated in the NOVEMOR project (Vasista et al., 2016). The first concept was developed by DLR following a similar workflow previously used in the development of the ‘active rib’ concept (Vasista et al., 2015). First, the morphing skin was designed to optimally match the required target shape, followed by the topology optimisation of the pertinent compliant mechanism and the actuator support. A demonstrator was built and tested on a test bench focusing on the performance of the compliant mechanism in combination with different morphing skins of various thicknesses. However, other loads, such as consolidated aerodynamic loads, were not applied during the test. The second concept was developed by De Gaspari and Ricci (2011, 2014) using a two-level approach with a significant difference being that in the first level the droop nose shape is optimised for maximum
Just recently Rudenko et al. (2018) presented an extremely deformable leading edge concept. The design of the inner morphing mechanism of the concept is based on the methodology developed by DLR within the SADE project. Nevertheless, the main novelty of the concept is the introduction of the ‘curvature morphing skin’ which is characterised by highly anisotropic behaviour therewith allowing the concept to achieve a high degree of droop. A section of the full-scale leading edge was built and tested on the test bench. Moderate agreement between finite element simulations and experiment was observed. Sensitivity to manufacturing inaccuracies along with simplified numerical modelling of skin-stringer connections was identified as the primary origin of discrepancies.
The morphing leading edge demonstrator presented here was developed during the Leading Edge Topology Optimised Design and Demonstrator (LeaTop) project as part of the EU CleanSky Green Regional Aircraft (GRA) initiative. Similar to previously presented projects, one of the objectives was to develop, manufacture and test a functional morphing leading edge demonstrator that can effectively transition from its original configuration, the cruise shape, to its target configuration, the high-lift shape, while subjected to full aerodynamic loads. In addition, an important objective of the project was to address also more practical considerations which limit the integration of such morphing devices in commercial aviation at a larger scale and demonstrate a feasible solution. Such considerations are whether an accurate cruise shape can be maintained during flight; how to transfer aerodynamic loads from the morphing skin of the leading edge into the wing primary structure without causing excessive stress concentrations at the skin-spar connection; and whether conventional materials, such as aluminium, can be used for the morphing skin, therewith mitigating bird-strike and de-icing issues commonly associated with composite structures.
The demonstrator was designed by Mauchle et al. (2012) using an integrated aeroelastic-morphing framework developed by Thuwis et al. (2010, 2011) at Delft University of Technology. The framework allows for concurrent optimisation of the morphing skin properties and the morphing mechanism in order to meet the required target shape as accurately as possible while, at the same time, also minimising strains in the morphing skin and corresponding reaction and actuation forces in the morphing mechanism. The design is performed in two steps. First, topology optimisation is used to optimise the skin properties and obtain the initial mechanism layout. In the second step, a detailed design is performed to account for practicalities such as mechanism connection to the skin and pertinent offsets due to the finite dimension of mechanical components in the load introduction points.
This work presents the experimental evaluation of the morphing leading edge concept with respect to the considerations addressed in the LeaTop project. The experiments were therefore divided into three parts: strain evaluation in the airfoil skin, morphing mechanism locking and load transfer capability, and shape measurements of the morphed skin. In addition, a finite element model (FEM) of the demonstrator was developed to compare with the experimental results and to obtain better insight into the behaviour of the demonstrator.
The article is structured as follows. First, the experimental setup and the employed measurement techniques are presented. Next, the FEM is presented followed by the discussion of the obtained results. Finally, all the findings and conclusions are summed up in the ‘Conclusion’ section.
2. Experimental setup
The experimental work was divided into three parts: evaluation of strain in the airfoil skin, locking capabilities of the morphing mechanism and shape measurements of the morphed skin. Each part is designed to test the demonstrator for important requirements that were imposed during the conceptual design process. Strain measurements validate the requirement that the airfoil skin undergoes primarily bending deformation during the morphing process to keep the induced strain as low as possible. Evaluation of the locking capability of the morphing mechanism validates the ability of the demonstrator to preserve the shape of the airfoil under aerodynamic loads at cruise conditions. Within this scope, the ability of the morphing mechanism to transfer the applied loads efficiently from the airfoil skin into the main spar is also tested. Finally, shape measurements are used to assess to what extent the morphing mechanism can match the target shape which was prescribed during the design of the demonstrator.
2.1. Morphing leading edge demonstrator
The morphing leading edge demonstrator was designed as a section of high-lift system of the Clean Sky GRA Integrated Technology Demonstrator (ITD) (Mauchle et al., 2012). The selected section was positioned

LeaTop morphing leading edge demonstrator.
Overview of the experimental setup.
A Linak LA 23 series linear actuator was used for actuation of the demonstrator. The selected actuator is capable of developing 2500 N of actuation force and has a rated maximum stroke of 49 mm. The actuator is equipped with a hall potentiometer sensor in order to monitor its extension and a simple on/off electronic control switch to extend or contract its pushrod (Linak, 2012). Furthermore, by monitoring the electric current consumption of the actuator, the actuation force could be estimated. Therefore, a LabView virtual instrument was developed to drive and monitor the electric actuation signal provided to the actuator. This way, data about the actuation behaviour, namely actuator’s extension, actuating speed and required actuation force, were obtained.
2.2. Strain measurements
The airfoil skin was instrumented with a strain measurement system. A Rayleigh backscattering fibre optic distributed sensing system (ODISI B from LUNA Technologies; Luna Innovations Inc. 2012) was employed in order to acquire as much information as possible on the strain distribution along the perimeter of the leading edge skin (Güemes et al., 2010). The integration of the fibre optic strain sensor is depicted in Figure 1.
A fibre optic sensor was attached to the airfoil skin using epoxy resin. The Rayleigh Backscattering system allows for a strain measurement every 1.25 mm along the length of the optical fibre up to a total length of 10 m with a maximum data acquisition rate of 23.8 Hz. Thus, it was decided to glue the fibre along the entire perimeter of the leading edge both on the inner and outer surface of the skin. Placement of the sensing fibre along the entire perimeter of the leading edge enabled a detailed analysis of bending and tensional deformation of the airfoil skin during morphing. On the other hand, placing the fibre symmetrically on both sides of the demonstrator allowed for detection of potential misalignment or asymmetry in the morphed demonstrator. A schematic path of the sensing fibre is depicted in Figure 1. The fibre started at location 1 and then followed the inner contour of the skin up to location 2 where it was guided to the outer side of the skin. At the stiffener locations, the fibre had to be detached from the skin, guided over the stiffener and then re-attached to the skin. The measurements over these segments of the fibre were discarded since they did not correspond to the actual strain in the skin. The fibre placed on the outside of the skin then followed the contour of the skin up to location 3 where it was guided to the far edge of the skin and then placed in the reversed order.
2.3. Shape measurements
The deformed shape of the airfoil skin was measured using photogrammetry. A digital camera was used to photograph the morphing demonstrator from a fixed position. It was ensured that the camera was far enough from the demonstrator such that the perspective distortions were negligible and that the same scaling factor could be used all around the perimeter of the skin. In order to enhance the visibility of the edge of the skin, the edge was painted in orange and masked with grey foam as can be observed in Figure 2(a). The captured images were then converted to black and white images. Since the skin edge was coloured orange, it was decided to substitute the red tones with the white colour during the conversion of the images. This way a good contrast was achieved between the orange skin edge and the grey masking foam. An image recognition algorithm, developed in MATLAB, was then used in order to process the contour of the black and white images and thus extract the contours of the morphed skin. An example of a processed image is shown in Figure 2(b) where the result of the edge recognition process is plotted by a red line. One can observe that the skin edge along with the stringers was correctly identified along the entire perimeter of the edge.

Edge detection of the morphing demonstrator: (a) original picture and (b) processed picture.
2.4. Morphing mechanism locking
Morphing mechanism features several 4 bar sub-mechanisms which were realised by adding links
It was important to demonstrate that the morphing mechanism was capable of transferring the applied loads to the main spar effectively. Due to its design, the mechanism is effectively locked in position if the actuator is not moving, which makes the mechanism a statically determinate system. Therefore, it is sufficient to measure a displacement of one node in the mechanism in order to assess how the mechanism copes with the applied load. A small or negligible displacement of the measured node under applied load indicates that the mechanism is effectively locked. The displacement of the selected member was measured using a linear voltage displacement transducer (LVDT). The position of the LVDT and the monitored mechanism member is shown in Figure 3(a). Representative distributed aerodynamic forces were reduced to a set of four discrete load vectors

Setup for testing of mechanism locking: (a) mechanism displacement measurement using an LVDT and (b) test rig for application of condensed aerodynamic loads.
Maximum discrete applied loads.
Total maximum applied load in comparison to aerodynamic load.
In addition to the displacement of the mechanism, the deformation of the airfoil skin in terms of induced strains was also measured at a number of critical locations using strain gauges. In this case, strain gauges were used instead of the optical fibre since it was easier to place them close to the stiffeners. The selected locations were close to the load introduction points, at the top and bottom skin-spar connection and at the tip of the leading edge where the maximum strain was expected according to the finite element analysis. The exact location of the placed strain gauges is also shown in Figure 4. Strain gauges were placed in pairs on both the inner and outer surface of the skin to be able to differentiate between axial and bending strains.

Strain gauge locations and load introduction points.
3. Finite Element Model
For the sake of comparison, the morphing leading edge was modelled in a commercial finite element (FE) analysis tool ANSYS Mechanical®, Academic Research, Release 15.0. The geometry of the FEM was defined using the 3D CATIA® model previously used for the manufacturing of the morphing leading edge demonstrator. To ensure that the CATIA® model and the actual demonstrator were indeed similar in dimensions, the actuation mechanism of the demonstrator was measured and compared to the CAD model.
The FEM was simplified by cleaning-up and simplifying the underlying CAD model where possible. For instance, in the original CAD model, the stringers and the load introduction joints were drawn as separate components and bolted to the skin. In the FEM, these components were merged and attached to the skin using the bonded contact provided by ANSYS Mechanical®. Remaining attachment holes in the skin and stringers were eliminated from the FEM, allowing for a quality quadrilateral mesh.
Members of the actuation mechanism, load introduction joints and stringers are much stiffer than the airfoil skin; hence, they were modelled as rigid. Moreover, the skin thickness is an order of magnitude smaller than the in-plane dimensions. Consequently, it was modelled as a shell. The skin was assigned isotropic material properties. Young’s and shear modulus were set to
SHELL181 elements were used to model the skin since large deflections are expected to occur during morphing. SHELL181 element is a quadrilateral four-node element with 6 degrees of freedom per node and it is well suited for applications dealing with geometric nonlinearities due to large deflections (ANSYS, 2015a). The rigid morphing mechanism was modelled using the TARGE170 elements connected together with multi-point constraint type of elements, MPC184. MPC184 elements were used to model mechanism hinges as revolute joints with unconstrained rotational degree of freedom about the hinge axis while the remaining degrees of freedom, the three translations and the remaining two rotations were constrained. In order to prevent the rigid body motion of the model, the main spar was connected to the ground using the fixed option of the MPC184 element with all 6 degrees of freedom being fixed. Contacts between the airfoil skin and the lateral stringers were modelled using pairs of CONTA174 and TARGE170 elements which are well suited for general rigid-flexible contact analysis (ANSYS, 2015a). The assembled FEM is shown in Figure 5. The cutout in Figure 5 shows the size of the shell elements used to mesh the airfoil skin. In addition, the finite elements used to model the actuation mechanism are overlaid over the mechanism members.

FEM model.
The mesh of the FEM was checked for convergence using the recommendations by Shah (2002). Shah’s recommendation is that the global error norm for the entire FEM, Criterion 1A, is less than 15%, that the local error norm in a high stress area, Criterion 1B, is less than 10%, and that the local stress coefficient of variance in a high stress area, Criterion 1C, is less than 7%. To evaluate the local convergence criteria, the high stress area was chosen between stiffeners 2 and 3 whose location is indicated in Figure 4. The Mesh metric of the FEM and the convergence criteria are summarised in Table 4. According to the listed results, the element size of 3 mm is sufficiently small to obtain a mesh-converged solution.
FEM summary.
Optimum aspect ratio: 1 (ANSYS, 2015b).
Recommended maximum aspect ratio: ≤20 (ANSYS, 2015b).
According to Shah (2002).
4. Results and discussion
Experimental and numerical results obtained during testing are compared and discussed in this section. Shape measurements are presented first, followed by the discussion of the strain results. The morphing mechanism locking test is discussed finally.
4.1. Shape measurements
A preliminary study of the measurement system was performed first to confirm that the demonstrator was actuated consistently and accurately. The study was also used to investigate the consistency and reliability of the photogrammetric method used in the experiment. The leading edge was morphed from zero to a prescribed actuator extension three times. Each time, a picture of the morphed skin was taken and analysed. The results are shown in Figure 6. It can be observed in Figure 6(a) that the edge recognition algorithm was accurate enough to detect both the outer and inner perimeter of the skin. Consequently, an additional averaging of the acquired measurements was required to obtain the final skin shape. The average shape of the skin is distorted near the lateral stiffeners due to the averaging technique used in the process which also accounted for the points on the stiffeners. The distortions can be observed in the cutout in Figure 6(a). The magnitude of the distortion is greatly dependent on the portion of the lateral stringer that is recognised from the pictures. These distortions represent the only significant source of uncertainty in the edge recognition process, which can be observed in the cutout in Figure 6(b). Consequently, it was decided to physically mask out the lateral stiffeners on the demonstrator before continuing the measurements. The measurement error along the perimeter of the leading edge is calculated as a standard deviation of the measurements at each arc length location. The resulting average measurement error with its standard deviation is

Shape measurement: (a) processed data obtained from the captured pictures of the morphed demonstrator and (b) consistency of the shape acquisition system at different actuator extensions.
In the following step, the demonstrator was compared to the prescribed initial and target shape of the droop nose which were used during the design stage (Figure 5 in Mauchle et al., 2012). The comparison is shown in Figure 7. Geometric comparison of the prescribed original and target airfoil and the morphed demonstrator is shown in Figure 7(a). The chord-normalised minimum Euclidean distance,

Comparison of the morphed leading edge to the original and target airfoil: (a) experimental results against original and target shape and (b) relative geometric difference.

Curvature comparison between the prescribed and measured leading edge shape: (a) actuator extensions of 0 mm and (b) actuator extensions of 39.2 mm.
For a fully retracted actuator, at 0 mm, the demonstrator shape corresponds very well to the prescribed original airfoil everywhere except at the vicinity of the top and bottom stiffener, which is also clearly observable by the two peaks at arc length of 340 and 740 mm in Figure 7(b). At the time when the demonstrator was manufactured, the only suitable actuator available in the market was the LINAK LA 23. This actuator was, however, 7 mm too long to allow for the demonstrator to reach its original state even when the actuator is fully contracted.
For a fully extended actuator, at 39.2 mm, the demonstrator follows the contour of the prescribed droop nose very well, especially over the top surface of the airfoil, beyond the arc length of 650 mm. However, there are two notable differences between the morphed demonstrator and the prescribed shape. First, the morphed demonstrator does not maintain the chord length relative to the prescribed target airfoil which leads to the gap of 0.6%
According to Mauchle et al. (2012), the skin was constrained as clamped at the top and bottom edge, at
Although the demonstrator was not investigated for aerodynamic performance, the comparison of the curvatures derived from the prescribed original and target leading edge shape and from the measured shapes is shown in Figure 8. Curvature importantly affects the pressure distribution over the wing leading edge region (Thuwis et al., 2010); hence, it is considered an important facet of shape matching. Three remarks can be made about the presented results. First, at the tip of the leading edge at arc length of 520 mm, the curvature of the demonstrator follows the curvature of the prescribed airfoils very well for both actuator extensions of 0.0 and 39.2 mm. Second, oscillations at arc lengths less than 100 mm and more than 900 mm are considered numerical artefacts due to numerical calculation of the second derivative of the leading edge shape with respect to arc length, which is needed for the calculation of the curvature. Numerical differentiation is very sensitive to slight variations in the measured shape if derivatives are close to zero. Finally, one can observe significant deviations of the measured curvatures from the prescribed curvatures in the vicinity of stiffeners on the top surface at arc lengths between 550 and 800 mm. These are attributed to a sudden stiffness change due to the presence of the stiffeners and the fact that the mechanism connection to the skin can only transmit forces and no moments. Such oscillations in the curvature of the leading edge can have detrimental effects on the aerodynamic performance of the proposed morphing concept and should be investigated in more detail.
In the second part of the shape morphing study, the demonstrator was compared to the FEM. The results for the maximum actuator deflection of 39.2 mm are presented in Figure 9. Figure 9(a) compares the morphed shape of the demonstrator to the initial and modified FEM. The pertinent chord-normalised Euclidean distance between the compared results along the perimeter of the leading edge is shown in Figure 9(b). In general, one can observe good agreement between the demonstrator and both the initial and the modified FEM. The average distance shown in Figure 9(b) is less than 0.15%

Comparison of the morphed leading edge to FEM at actuator extension of 39.2 mm: (a) experimental results against FEM and (b) relative geometric difference.
The main difference between the initial FEM and the demonstrator is present in the part of the skin up to arc length of 590 mm with the maximum deviation of 0.41%
The agreement between the modified FE model and the demonstrator is significantly improved with the maximum deviation being reduced to just 0.14%
During every morphing action, the actuator extension and the actuation force were monitored. Four different extensions and their respective actuation forces over time are shown in Figure 10(a) and (b), respectively. According to the measurements, there is an initial force build-up period, up to 1 s, at the beginning of every morphing action. During this period the actuation force is rapidly rising while the actuator is still not moving. The actuator must first reach the break-away force to overcome the static friction in the actuation mechanism. Once the mechanism starts moving, a sharp drop in the actuation force,

Actuation in time: (a) actuator extension and (b) actuation force.
The continuous decrease in the required actuation force is attributed to the kinematic properties of the morphing mechanism. While the actuator is being extended, the morphing loads are increasing. However, the leverage between individual mechanism members is also changing due to their relative rotations. Initially, the change in the mechanism leverage alleviates the loads on the actuator faster than the morphing loads increase.
Actuation force is also compared to the FEM results. Even though the material properties of the airfoil skin in the FEM simulations were selected to resemble aluminium, the overall agreement is not good. For most of the actuator extension, the FEM forces significantly underpredict the measured actuation force, which can be attributed, first, to the fact that friction in the mechanism joints was not taken into account in the FEM, and second, the actuation force was measured indirectly by measuring the current drawn by the actuator, which also includes all the internal losses in the actuator itself. However, at actuator extension of 39.2 mm, the FEM actuation force is almost equal to the measured actuation force.
Both the actuator extension and the actuation force were also analysed for the measurement error. The measured actuator extension was compared to the requested actuator extension. The average error with its standard deviation of
The assessment of the measurement error in the case of the actuation force was split in two parts. Over the force build-up period, up to
Next, the actuation force is integrated with respect to the actuator extension in order to obtain the overall actuation energy required to execute a morphing action. The actuation energy as a function of actuator extension is compared to the strain energy of the morphing skin obtained by the FE analysis in Figure 11. The strain-to-actuation-energy ratio indicates that, on average, only 35% of the actuation energy is used for morphing the airfoil skin, while the rest is lost in the actuation mechanism.

Actuation energy.
There are two further comments required regarding the actuation and strain energy shown in Figure 11. First, the actuation force over the force build-up period does not contribute to the calculated actuation energy since the actuator displacement remains zero. In reality, some electric energy is spent on powering the actuator to overcome the initial static friction of the mechanism. Hence, the overall energy spent on morphing will be higher than the predicted actuation energy. Second, the strain energy of the morphing skin was reduced by 1.1 J to account for the initial deformation of the morphing skin due to the actuator size as explained in the preceding section on shape measurements. The resulting strain energy therefore corresponds to the strain energy due to the actuator extension only.
4.2. Strain
Geometrical analysis of the morphing leading edge demonstrator showed good agreement with both prescribed original and target shape as well as with the FEM. Another critical aspect to investigate is how and to what extent is the airfoil skin deformed during the morphing process by analysing the induced strains. For the sake of clarity, only the results pertinent to the modified FEM are shown.
The consistency of the measured strains is investigated first by comparing the measurements on the inner and outer surface of the airfoil skin at two spanwise locations and by assessing the measurement error. The results for two different actuator extensions, 20 mm and 39.2, are shown in Figure 12. The strain is given as a function of arc length along the airfoil skin starting at the bottom of the airfoil next to the main spar. As mentioned, strains were measured at two spanwise locations. The corresponding results are hence labelled 1 and 2. Moreover, strains measured on the outer and inner surface of the airfoil skin are labelled out. skin and in. skin, respectively. The strain measured on the inner surface is multiplied with a factor of −1 to ease the comparison with the measurements on the outer surface of the airfoil skin. A factor of −1 is selected since the respective strains are expected to be of the opposite sign due to the premise of the actuation mechanism design that morphing should be mostly driven by bending deformation of the airfoil skin. Finally, one can observe that zero strains were measured on the inner surface along some sections of the leading edge perimeter. These measurements coincide with the locations where the fibre optic sensor is not attached to the structure as it has to traverse the lateral stiffeners.

Symmetry analysis: (a) actuator extension of 20 mm and (b) actuator extension of 39.2 mm.
The comparison of the measured strains pertinent to locations 1 and 2 shows that the strains match remarkably well for both outer and inner surface of the skin and both actuation extensions. The measured strains are almost identical in this regard. This leads to a conclusion that the morphing mechanism is rigid enough to provide exact and uniform skin deformation across the entire width of the leading edge without any unwanted skewing or twisting. Moreover, the magnitude of the strain on the inner surface of the airfoil skin represents an almost perfect reflection of the strains measured on the outer surface, which suggests that the induced morphing is indeed predominantly driven by bending.
At the maximum actuator extension of 39.2 mm, three strain measurements were acquired. The measurement error along the fibre optic was then estimated as the standard deviation of the measured strains across the three measurement sets. The average error with its standard deviation was found to be 10 micro-strains ± 113 micro-strains. The high standard deviation of the estimated measurement error is mainly due to the high measurement error in the regions of high strain rates, between arc lengths of 550 and 620 mm. However, if the maximum 1% of the measurement errors is omitted from the error analysis, the pertinent standard deviation drops to ±1 micro-strain.
Strain data obtained by the fibre optic sensor on the outer surface of the airfoil are compared with the results from the FEM. Raw strains at two different actuator displacements are shown in Figure 13. The measured and the numeric strains are, in general, in good agreement except near the stiffeners and close to the bottom connection to the main spar, at the arc length of 50 mm. In these areas, the numeric results predict higher strain values. Strain is also concentrated closer to the stiffeners. Such behaviour is attributed to the stiffeners being modelled as rigid, forcing the skin to strain more. Otherwise, good agreement was expected since the measured and the calculated morphed geometries of the leading edge presented in Figure 9(b) also show good overall agreement.

Comparison of measured and FEM strains: (a) actuator extension of 20 mm and (b) actuator extension of 39.2 mm.
Strains corresponding to the actuator extension of 39.2 mm are also mapped to the morphed airfoil in Figure 14 highlighting regions of high deformations. Strain values along the outward normal of the airfoil contour are considered positive. It shows that the highest strains, both positive, tensile strains as well as negative, compressive strains, are encountered at the very tip of the leading edge between the first and the second stiffener. High compressive strains are also encountered between the second and the third stiffener. It is noteworthy that Figures 13 and 14 present only the strains on the outer surface, meaning that the nature of the underlying deformation of the airfoil skin, such as stretching or bending, cannot be deduced based on these data only.

Measured and FEM strains mapped to the morphed airfoil skin.
Further analysis of the available strain data was performed to obtain better insight into the amount of bending and axial strain present in the induced deformation. Typical strain behaviour in the airfoil skin during morphing is shown in Figure 15. Since the skin was relatively thin and there were only two strain measurements available through the skin thickness, namely, on the inner and outer surface of the skin, a linear through-thickness behaviour of the strains was assumed. Consequently, axial strain,

Strain decomposition.

Comparison of bending strains: (a) actuator extension of 20 mm and (b) actuator extension of 39.2 mm.

Comparison of axial strains: (a) actuator extension of 20 mm and (b) actuator extension of 39.2 mm.
Experimental and numeric bending strains exhibit very good agreement, except in the very vicinity of the stiffeners, where the FEM tends to over-predict the experimental results. According to the FEM, the bending deformation is more localised near the stiffeners, which is reflected by narrower and sharper peaks in the FEM strain plots.
Agreement in the axial strains is worse. The numeric strains are close to zero showing little variation along the perimeter of the airfoil skin, except at the stiffener locations, where the numerical strains show sharp peaks which are explained by sudden stiffness change due to rigid stiffeners. As a result, there is no strain induced on the inner surface where the skin attaches to the stiffeners. Therefore, following the procedure outlined in Figure 15 leads to high axial strains which equal to half of the strain obtained on the outer surface of the airfoil skin. Experimental strains show more variation along the perimeter of the leading edge especially for the higher actuator displacement of 39.2 mm. Relatively high compressive strains are observed at the arc length of about 500 mm, between the first and the second stiffener. In combination with the stiffness and geometrical properties of the skin panels in such regions, high compressive strains might indicate a buckling critical area which would have to be accounted for in the design process. Nevertheless, no buckling was detected during the current set of experiments. Finally, one can observe some oscillations in the experimental strain, such as those at arc lengths of 88 and 476 mm. These oscillations appear in the areas of very high strain rates as can be observed in Figure 12. Hence, a slight misalignment of the measurement points along the arc length between the outer and the inner surface can lead to a significant error in the acquired axial strain. These oscillations are therefore largely attributed to the misalignment error in positioning the inner and outer strains with respect to each other.
Comparing the magnitude of the axial and bending strains, the axial strains appear to be an order of magnitude smaller than the corresponding bending strains which confirms one of the main design premises that drooping of the leading edge should be driven predominantly by bending deformation in order to be able to achieve large deformations while keeping the actuation forces as low as possible.
It is important to point out that no aerodynamic or equivalent distributed loads were applied on the airfoil skin in the current experiments. Therefore, it remains to be investigated what the effect of these loads in combination with morphing deflection on the overall strain values induced in the morphing skin is.
Although fatigue life assessment of the presented demonstrator is beyond the scope of this article, given the strain levels observed at a maximum actuator displacement of 39.2 mm, a remark about possible fatigue issues is necessary. For reference, the 7075T6 aluminium comprising the airfoil skin yields at 8860 micro-strain (U.S. Department of Defense, 1998). The maximum axial strain of 681 micro-strains represents less than 10% of the yield strain. It is therefore not a concern. On the other hand, the maximum bending strain of 3403 micro-strains is of comparable magnitude as the yield strain. Even though the flexural fatigue is less severe than axial fatigue due to the lesser volume fraction of the material being exposed to high strain/stress levels, the airfoil skin could be susceptible to low-cycle fatigue (Manson and Muralidharan, 1987). However, in the current case, the airfoil skin was manufactured from a constant thickness isotropic material. Therefore, the high bending strains could be resolved, if necessary, by allowing for a variable thickness and stiffness tailoring of the morphing skin along the perimeter of the leading edge.
4.3. Locking of the morphing mechanism
Experimental results in terms of the nodal displacement and skin strain due to external load are presented in this section to evaluate the locking of the morphing mechanism. During this experiment, a total load of 2500 N was gradually applied via a system of pulleys and weights as shown in Figure 3. The experiment was repeated four times to be able to assess the consistency of the mechanism behaviour.
The displacement of the selected node is presented in Figure 18. The results show two interesting features. First, at the maximum load, the mechanism undergoes a maximum displacement of 0.4 mm. Second, the displacement measurements are very consistent over the performed set of experiments except for Experiment 1. Experiment 1 deviates from the others, which is attributed to the presence of some initial free-play in the mechanism before applying the load for the first time. The measurement error as a function of applied load was estimated as standard deviation of the measured nodal displacement across the Experiments 2, 3 and 4. Experiment 1 was excluded from the error analysis due to being an outlier. The average error with its standard deviation was found to be

Deformation of the morphing mechanism.
From the collected data it is not possible to discern whether the 0.4 mm displacement is due to the mechanism deformation or due to the remaining free-play in the joints of the mechanism. Such mechanism displacements can have a detrimental effect on the ability of the mechanism to maintain an accurate aerodynamic shape of the leading edge over the range of operating conditions. In addition, changing values of mechanism free-play can also lead to higher actuation forces due to the increased friction. The effect of mechanism stiffness and production tolerances therefore remains an important aspect to include in the future design studies.
During the application of the load, strains were also measured using strain gauges at the selected locations on the morphing skin as shown in Figure 4. The measured axial and bending strains are depicted in Figure 19(a) and (b), respectively. Strain gauge 2 and 8 were damaged during the experiment, hence their values are not shown. For the remaining strain gauges, it can be observed that the maximum magnitude of the bending and axial strain at the maximum applied load reaches roughly 2000 and 350 micro-strain, respectively. Both peaks are reached at the location of strain gauge 4, which is located near the most forward load introduction point, where also most of the consolidated aerodynamic load was applied. Moreover, the bending strain measured by strain gauge 1, close to the top skin-spar connection, is only 100 micro-strain which is an order of magnitude less than the maximum bending strain observed by the strain gauge 4. Based on the obtained results, one can conclude that most of the applied load was taken by the morphing mechanism while the airfoil skin remained relatively unstrained. This demonstrates the capability of the mechanism to provide structural support to the morphing skin. As a result, such mechanism design would allow for a more compliant morphing skin and would prevent the formation of stress/strain concentrations in the vicinity of the skin-spar connections due to external loads.

Measured strains at applied consolidated load of 255.2 kg: (a) axial strain and (b) bending strain.
Finally, the measurement errors for both axial and bending strains were evaluated similar to the nodal displacement of the mechanism. The average errors in axial and bending strain were found to be 8 micro-strains ± 11 micro-strains and 54 micro-strains ± 62 micro-strains, respectively.
4.4. Discussion and future work
The performed experimental investigation improves the understanding of the structural and morphing performance of the tested morphing leading edge demonstrator and to a certain extent validates the design process employed to design the demonstrator. Nevertheless, a number of important issues with the current demonstrator itself as well as with the design process in general have been identified which still have to be addressed in order to be able to improve the design of such morphing devices and increase their maturity.
One of the most important aspects that remains to be investigated is aerodynamic performance of the designed morphing leading edge. Such an investigation should comprise an aeroelastic analysis, either numeric or experimental, of the leading edge shape in its cruise as well as high-lift configurations in order to provide at least three insights. First, how does the morphing skin deformation under the aerodynamic loads affect the aerodynamic performance of the leading edge in both cruise and high-lift configurations. Second, how do the aerodynamic loads influence the overall strain levels in the airfoil skin during morphing. Third, a quality metrics for the achieved morphed shape such that one can gauge the quality of the actuation mechanism design in a more objective way.
From structural and morphing standpoint, an important concept for future investigations is stiffness and thickness tailoring of the morphing skin. As shown by Thuwis et al. (2010), tailoring can be used for better matching the prescribed target shape as well as for mitigating the high strain levels due to morphing. Both aspects, however, remain to be validated.
Moreover, the comparison of the curvatures between the measured and target airfoil shape revealed that the stiffeners significantly distort the curvature of the morphing skin in their vicinity which could be critical from aerodynamic point of view. Therefore, it is important to investigate whether it is possible to achieve a more blended skin-stiffener connection using skin tailoring in order to suppress this kind of curvature distortions.
In terms of actuation mechanism design, several aspects such as free-play, friction and scalability remain to be investigated. In this study, some free-play was detected in the mechanism joints. It is not yet clear, first, what is the effect of such free-play on the accuracy of the leading edge shape in both cruise and high-lift configurations, and second, what amount of free-play is still acceptable. Moreover, the effect of free-play in combination with friction on the required actuation force once the leading edge is subjected to aerodynamic loads must be also investigated. Finally, the demonstrator was built as a two-dimensional (2D) section of a wing. Hence, it remains to perform a scale-up study to investigate aspects such as the required pitch of the actuation mechanism along the wing span and three-dimensional (3D) effects such as skin warping when deformed.
Finally, to be able to advance the technology readiness level of such a concept, one needs to consider also operational and certification requirements. A fatigue and life cycle analysis is required in order to ensure sufficient longevity of the morphing skin and actuation mechanism. In addition to stress/strain effects, such an assessment should also address effects related to material ageing and abrasion due to exposure to environment. In terms of certification requirements, bird-strike and de-icing, among others, have to be considered. These issues have been already addressed to a certain extent within the SARISTU project (Wölcken and Papadopoulos, 2016). Nevertheless, it remains to be investigated whether any of the proposed solutions can be adapted to the currently investigated concept.
5. Conclusion
Experimental evaluation of the morphing leading edge concept is presented in this article. In addition, the experimental results are compared to the FE results. The experimental trials were designed to test the morphing demonstrator with respect to the most important requirements imposed during the design stage of the LeaTop project. The primary design requirements were as follows: (1) capability of matching the prescribed target airfoil, (2) ensure that morphing is dominated by bending deformation and (3) ensure that the transfer of aerodynamic loads from the skin to the main spar is achieved via the morphing mechanism. The acquired results show that the morphing leading edge demonstrator successfully meets all of the design requirements.
Airfoil shape measurements reveal a good match with the prescribed target shape of the demonstrator as well as with the morphed shape obtained using the FEM. The only significant difference with respect to the prescribed target shape was detected near the bottom skin-spar connection due to the spar design protruding in the leading edge. A noticeable difference between the experiment and the initial FEM was observed in this area as well which was attributed to the airfoil skin peeling off from the skin-spar connection. However, once the FEM was updated to account for this effect, the match between the experiments and the FE results was very good along the entire perimeter of the leading edge for all actuator displacements.
Strain measurements show that morphing is dominated by bending the airfoil skin rather than stretching it. The axial strains are at least an order of magnitude smaller than the bending strains over most of the airfoil skin. Moreover, comparison with the FE results shows very good agreement in general for both bending and axial strains at all levels of actuator displacements. Substantial differences are only observed close to mechanism attachment points which are attributed to the stiffeners being modelled as rigid allowing for zero strain on the contact surface.
The measured actuation force shows that a significant break-away force is required to overcome static friction in the actuation mechanism. Comparison of the measurements to the numerical results reveals that the actuation forces significantly disagree even though the morphed shape of the leading edge and the corresponding strains agree very well. There are two main reasons for the observed differences between the numeric and experimental actuation force. First, the actuation force was measured indirectly by measuring electric power consumed by the actuator which in addition to morphing loads and friction forces accounts also for all the electrical losses in the actuator. Second, the FEM did not account for friction in the actuation mechanism.
Finally, in the case of the morphing mechanism locking capabilities, the results show that the mechanism is effectively locked under consolidated aerodynamic loads such that only minimum displacement in the measured mechanism node is observed. Strain measurements show that the mechanism can support the skin by absorbing and transferring most of the applied load to the main spar as indicated by the low strain values. Consequently, a more compliant morphing skin can be employed and the stress/strain concentrations close to the skin-spar connections can be reduced.
