Design and flight test of a linear parameter varying flight controller

Future aircraft generations require improved performance and efficiency to enable a reduced environmental footprint. To acquire this goal, for example new material and wing concepts are perused at the moment by the aircraft industry. These developments, which include aspects such as over-actuation and lowly damped flexible modes, give rise to more complex, multi-objective control problems. One candidate method, which delivers a solution to these problems for the whole flight envelope, is linear parameter varying (LPV) control. It naturally incorporates the controller scheduling in the synthesis process, guarantees stability and robustness over the entire parameter envelope, and enables intuitive multi-objective, multiple-input multiple-output (MIMO) controller designs. This paper proves the concept of LPV control in practice: The paper presents and discusses the LPV controller design process, simulation results, motion simulator test and finally, the in-flight validation of the control system on a Cessna Citation II aircraft. The developed inner loop controller structures are inspired by classical flight controllers used on state-of-the-art fly-by-wire airliners. The longitudinal aircraft motion is augmented with load-factor command and the lateral motion controller features a roll rate command with attitude hold behavior. The control laws are validated in flight by automated and actual pilot inputs with respect to functionality, flying and handling qualities. Test results are encouraging with the provided key findings and lessons learned aiming to provide a simplification for future LPV flight controller development and testing campaigns.


Introduction
The design of flight control laws for transport aircraft is still mostly relying on classical control methods, such as proportional-integral-derivative controllers in combination with straightforward gain scheduling. This is due to the stringent certification requirements but also due to the good performance characteristics of these classical control approaches. Future aircraft controller designs will, however, need to incorporate multi-objective design problems, as well as tackle their multiple-input multiple-output (MIMO) nature.
One modern control design technique, which showed already promising results in flight is linear parameter varying (LPV) control [16,17]. In addition to robustness against model uncertainties, the LPV technique results in a controller which is already scheduled over flight envelope parameters [1]. Moreover, the LPV design can be seen as an extension to a previously conducted H ∞ control test campaign, which was performed using an unmanned aerial vehicle [25].
For this work, a CS 25 certified Cessna Citation II is chosen as testbed. On this aircraft, extensive experience of testing different fly-by-wire control laws has been gained during previous test campaigns for (Incremental) non-Linear Dynamic Inversion control [8]. Thus, the aircraft is a perfect test bed to develop new control algorithms on real aircraft.
The contribution of this paper is to extensively present the whole development process of an LPV controller, from its initial model based controller synthesis [26] to the inflight validation [27]. Therefore, LPV control for both axes is assessed within this work. LPV controller synthesis is employed to derive robust and scheduled control laws for both longitudinal and lateral augmentation of the test aircraft. The 1 3 LPV technique is an advancement of classical robust control approaches and therefore allows achieving satisfactory robustness against model uncertainties via the closed loop shaping technique. Sufficient robustness with respect to model uncertainties and variations in the weight and balance configuration is considered in the control design. Moreover, the additional performance objectives defined herein are satisfactory handling qualities for augmented flight. Specifications for satisfactory flying qualities defined in [24] are considered. A strategy on how these design requirements are transformed into frequency-domain weighting filters used by the LPV closedloop shaping synthesis is discussed. For testing and evaluation of the controller, the DLR Robotic Motion Simulator within the Institute of System Dynamics and Control [3] is used for software in the loop (SIL) evaluations before transferring the control laws to the test aircraft for Hardware in the Loop (HIL) testing and flight experiments.
Firstly, an introduction of the used means of LPV controller synthesis is given in Sect. 2. For the design of the control laws, an LPV representation of the nonlinear aircraft model of the Citation II is obtained within Sect. 3. Following, an LPV controller is developed based on the gridded LPV model. Sect. 4 gives insight in the process of control law validation leading towards the flight testing as well as test execution and a summary of results.

LPV control synthesis
In this section, the LPV approach used for control design is described. As specifications on the longitudinal and lateral flying qualities are given in [24] already in the frequency domain, this is a promising starting point to derive a closed loop shaping process, which results in good handling qualities. H ∞ control [23] is a useful method as it allows to 'shape' the closed loop frequency response. As extension to this, LPV control synthesis permits to use a similar scheme for LPV models [18,28]. A considerable advantage of the resulting controller is the ability to be directly scheduled with the varying parameter vector while guaranteeing stability and performance for all parameter trajectories within the bounds of its defined values and rate of change.
Firstly, LPV models are introduced in Eq. (1) similar to Linear Time Invariant (LTI) models with the addition that all four matrices are dependent on a scheduling parameter vector (t) . Each entry of the parameter vector (t) is a function of time with allowable trajectories, defined by limits on each parameter as well its derivative.
To explain the idea of LPV control, the induced L 2 norm of an LPV system G from input d to output z over all allowed trajectories in a set P, assuming zero initial conditions, is defined in [28][Def. 3.

3.1] as
The induced L 2 norm expresses the maximum input-output gain of G , i. e., the largest amplification of L 2 input signals over all frequencies, input/output directions, and admissible parameter trajectories. For the design of linear parameter varying control laws, a closed loop interconnection including a weighting structure as seen in Fig. 1 is introduced to shape the closed loop response of the system. In the generalized plant of Fig. 1  (3) In Eq. (3) K y, is the feedback part and K r, the feed-forward part of the controller, G is the transfer matrix from command inputs to plant outputs. Further, S = (I + G K y, ) −1 is the sensitivity function, and S i, = (I + K y, G ) −1 and T i, = K y, G (I + K G ) −1 are the input sensitivity and input complementary sensitivity functions, respectively T r, = K r, G S complimentary sensitivity on reference signal. W e (s) and W u (s) are the diagonal output weighting matrices, and W d i and W d o are the diagonal input weighting matrices to be chosen in the design process. Note that the aircraft models G are scaled as proposed in [23] by input-output scaling. As indicated in Fig. 1, W e (s) , W h (s) and W u (s) are selected as dynamic filters, while all other filters are constant. W h (s) denotes a handling quality filter in the form of Eq. 5 and allows to set a design response in terms frequency h and damping h .
The entries for the shaping filter W u are selected as first-order transfer functions with unit gain up to the available bandwidths a,i for i = 1, … , n of each control input and approximately differentiating behavior beyond that frequency. In W e only the first element is dynamic to weight the tracking behavior, while all other channels are set to constant values to weight the outputs. As in this paper integral behavior is demanded in the tracking channel, the choice for the first entry is a first-order transfer function with approximately integral behavior up to the desired bandwidth b , which reduces the sensitivity up to this defined bandwidth. The weighting filters are selected in such manner that the general design goals for pitch and roll control, which are taken from the level 1 flying quality specifications for piloted aircraft [24], will be met. The demands for the flying qualities are defined in frequency domain for both longitudinal and lateral motion. For example, for design of the weighting filters of the pitch controller the control anticipation parameter (CAP) is used as a design parameter and sets a desired short period (SP) frequency. A more detailed description can be found in Sect. 3.2. Further, the controllers shall meet common requirements for gain and phase margins (6 dB, 45 • ).
The resulting parameter dependent controller is synthesized solving a convex optimization problem [29]. The obtained controller therefore minimizes the upper bound on the induced L 2 -norm of the interconnection in Fig. 1 Assumptions necessary to solve this state feedback and observer problem (no full state feedback available) are described in [21].
As proven in [21,28], an output feedback controller K that stabilizes the closed loop interconnection G cl, = F L (P, K) and guaranteed ‖G cl, ‖ < can be found if there exist symmetric positive definite matrix functions X ∶ ℙ → ℝ n×n and where ⋆ denotes symmetric completion and As Eq. 9 also depends on the rate of change in the parameter ̇ , for the search of the functional relations X( ) and Y( ) mth order polynomial base functions of the form are selected [29]. The choice of the base function parameters is up to the designer. It is useful to choose the base functions as simple as possible, as the number of unknowns in these functions correlates to the unknowns in the Linear Matrix Inequalities (LMIs) in Eq. 8 and therefore increases problem complexity. The obtained solutions for X( ) and Y( ) allow to specify a controller K according to Eq. 6. A detailed calculation of the four controller matrices based on X( ) and Y( ) can be found in [29,Sect. 4]. Design examples for LPV controllers are described in literature, for example in [19,29]. The LPVTools MATLAB toolbox [9] implements the solution to the LPV control problem in a similar way as described above.

Modeling and control design
For the model-based control design strategy a mathematical model of the considered test-aircraft is derived using an available non-linear simulation model. The aircraft, the process of deriving the LPV design model and finally, the designed controllers, including a first design-verification, are discussed in this section.

Aircraft model
The Cessna Citation II depicted in Fig. 2 is a research aircraft operated by the Delft University of Technology. The basic parameters of the aircraft are listed in Table 1.
A simulation model of the investigated aircraft was developed using the Delft University Aircraft Simulation Model and Analysis Tool (DASMAT) [12] and allows simulation and control design within the MATLAB/Simulink environment. The used DASMAT model was originally generated for a Cessna Model 500 Citation; nonetheless the modified simulation model shows good compliance with the Citation II [10]. The used aircraft model features standard six Degree of Freedom (DOF) nonlinear equations with aerodynamic lookup table data based on [12]. Note that the aircraft features relatively low aspect ratio paired with high stiffness due to a metal structure. Thus, the rigid body approximation holds especially in the envelope flown during controller assessment (medium altitudes and airspeeds) and effects like aeroelastic modes can be neglected. The actuators were modeled as second-order transfer functions according to [5]. The sensors are modeled as a combination of time delay, bias and noise. The basic flight control system is included and features a servo controller, which is able to set the actuator positions commanded by the flight control laws.
To obtain an LPV model of the Citation II aircraft, the nonlinear model is trimmed and linearized at a grid of operating points. Thus, a model according to Eq. 1 is obtained. As scheduling variable, dynamic pressure is chosen, where is air density and V true airspeed. For longitudinal and lateral dynamics, the LPV model of the full plant is fractioned into two models, as for controller design it is appropriate to consider the two axes as decoupled [6].
In the scheduling parameter q is considered in a range of q ∈ [ 2.5, 6.5] kPa, which corresponds to a range of true airspeeds from 76 m s −1 to 129 m s −1 at an altitude of 4500 m (Flight Level (FL) 150)). The parameter rate is bounded to be within ±500 Pa s −1 . This corresponds to a rate of change in airspeed of approximately 7 m s −1 when at FL 450 with initial airspeed of 90 m s −1 . For the controller development, an equally spaced parameter grid with 18 points was selected and verified against a denser grid of more than 50 points. Figure 3 depicts the poles of the open loop plant 4 × 4 Linear Time Invariant (LTI) models of the aircraft at four equidistantly chosen grid points, with dynamic pressure values of [2.5, 3.5, 5.0, 6.5] kPa. From the figure, the conjugate complex poles of SP mode can be read with frequencies ranging from 2 to 3 rad s −1 and a damping ratio of approximately 0.45. The Dutch Roll (DR) motion is located at a similar frequency range, but with a damping ratio of slightly less than 0.2, this mode is only poorly damped. The roll time constant t r can be read to a range from 0.33 s at highest dynamic pressure to 0.5 s at the lowest dynamic pressure value. The phugoid and spiral poles have large time constants and are located near (0,0) in Fig. 3.
The following control design is carried out with the LPV models for longitudinal and lateral dynamics obtained here.

Longitudinal control laws
For the longitudinal control augmentation, load factor command n z,cmd is selected as tracking reference for control design. The controller commands the maneuver load factor of the aircraft set by the pilot's pitch control input.  As the n z,b sensor measurement includes gravity, the feedback needs to be corrected by the aircraft's pitch and roll attitudes ( , ): This corrected value n z,man results in a value of 1 g, independent of the aircraft's attitude and flight path. This design approach has the advantage that a zero command input results in an attitude hold behavior, which is convenient for augmented flight. Therefore, a neutral stick input is related to a 1 g flight with flight path angle hold. As the test aircraft is not equipped with an auto-throttle system, manual adjustment of the throttle is mandatory to keep the airspeed and dynamic pressure inside the test envelope bounds. To obtain a stable and fully controllable plant, control design only considers elevator input together with the SP model, with angle of attack and pitch rate as states and load factor n z as controlled plant output. Due to the missing auto-throttle, the 4th order longitudinal model is not fully controllable (attitude and path variables cannot be decoupled). For this reason, instead the 2nd order SP model is used for controller synthesis. When integrating the resulting SP controller into the full model, the result is a phugoid motion with two real poles of which one may be in the unstable region, but with a time constant larger than 30 s. This leads to the requirement (13) n z,man = n z,cg cos( ) cos( ) of manual adjustment of thrust settings accordingly to avoid stall or overspeed regions of the envelope. The limits of the controlled output are ± 0.3 g for the commanded load factor. The design process for the maneuver load factor requires r = n z,man , as the tracking reference variable. The measured outputs are the load factor in body z-axis and the pitch rate of the aircraft. Design goals for the longitudinal motion are defined in [24] via the control anticipation parameter (CAP), as well as a minimum SP damping of = 0.3 and an optimal SP damping value of = 0.7.
The CAP was originally defined in [4] and relates the immediate pitch acceleration to a control input and the steady state load factor which will be the result of a longitudinal control input: where SP is the SP frequency and n z, the (stationary) load factor resulting from a change in angle of attack. The value of n z, is characteristic for each aircraft, thus only SP can be influenced. The CAP has a goal value of one, which means that the pilot is able to determine the amount of stick input needed for a certain load factor easily from the pitch acceleration ̇q 0 , which is a direct result of an elevator deflection . The load factor n z in turn takes time to build up as it has two integral stages in between: where q is the pitch rate and the angle of attack. In case of a CAP smaller than one, the pitch acceleration after a control input is low and the pilot will have the tendency to underestimate the stationary load factor, thus giving higher pitch input as needed and overshooting the desired load factor. The handling quality level 1 boundaries and results of the chosen controller design for load factor control can be seen in Fig. 4. The resulting CAP is in the level 1 region close to the optimal value for the whole design envelope, which was achieved in this case by a constant handling quality filter with a frequency of 1.75 rad s −1 and damping ratio of 0.8. The load factor longitudinal controller features 12 states and shows integrating behavior up to a frequency of 2 rad s −1 . Figure 5 depicts the sensitivity and complementary sensitivity functions of the load factor controller. The plots indicate a maximum frequency for reference tracking of 0.4 rad s −1 . Figure 6 depicts singular value plots of the synthesized controller K at four equally spaced grid points. It can be seen that until 2 rad s −1 the controller features integrating behavior, whereas constant/slightly increasing gain beyond until about 20 rad s −1 ensures robustness of the closed loop.
The time domain results of the linear n z command design exhibit rise times smaller than 4 s for low and smaller than (15) →̇q ∫ → q →̇∫ → → n z , 2 s for high dynamic pressures. In terms of robustness, the multivariate multiple-input multiple-output (MIMO) margin has a minimum at the upper end of the parameter envelope with values of 6.4 dB and 41 deg respectively, which can be considered as sufficient. The margins listed for the load factor control in Table 2 are a comparison of standard gain-, phase-, and delay-margins of loop cuts at the sensor or actuator signal(s) with more significant disk margins and robustness margins used in [22]. The disk margin is computed according to [2] and values for the phase of the disk margin vs. the scheduling parameter of the controller are depicted in Fig. 7. The drop of the margins for the high dynamic pressure region calls for a careful flight testing at the end of the envelope. As the maneuverability of the experimental flight control system (FCS) is limited by a relatively moderate power of the pitch servo, most maneuvers are flown in a low dynamic pressure region. When looking at frequencies, it can be noted that the frequency of the minimum disk margins in Fig. 7 is located several magnitudes above the ) and level 1 flying qualities region boundaries ( ) according to [24] are depicted

Lateral control laws
For the lateral inner loop the roll rate p and angle of side slip are chosen as control variables. In the chosen setup, the experimental sidestick on the right hand pilot's seat allows the control of a roll rate up to a bank angle of 27 • , beyond constant stick input is needed for higher bank angles. The maximum bank angle of 35 • shall never be exceeded and for a neutral stick the roll attitude shall return to 27 • . The aircraft is not equipped with a yaw input on the experimental FCS, thus the yaw axis can be assessed via differential thrust settings or slip inputs triggered by the flight test crew.
For the inner loop control laws designed in this section, flying quality specifications for lateral axis [24] request a maximum time of 2 s for a 30 deg bank angle change. Thus, a roll rate limit of p max = 15 deg s −1 is desirable, but due to limitations in fly-by-wire (FBW) actuation of the test aircraft only p max = 10 deg s −1 can be achieved.
The roll time constant t r , which defines rise time of roll rate, shall be less than one second, which allows a sufficiently fast buildup of the roll rate following the pilot's lateral input. For the Dutch roll (DR) mode, specifications request a minimum damping of = 0.08 and a minimum frequency of = 0.4 rad s −1 as well as the product of damping and frequency to be at least = 0.15 rad s −1 . The handling quality filter for the roll rate is implemented as a second order filter in the form of Eq. 5 with a frequency of SP = 5 rad s −1 and a damping of = 0.9, to provide tracking of a commanded roll rate. The same form of filter, with a lower corner frequency, is chosen for the angle of sideslip channel. For both filters the performance weight on the tracking error is represented by a transfer function with integrating behavior up to a frequency of 12 rad/s, which is sufficiently lower than the actuator bandwidth of 30 rad/s. The actuator weights have differentiating behavior above the actuator bandwidth. All other weights are selected as constants. The resulting sensitivity and complementary sensitivity functions are depicted in Fig. 8. For both reference variables, nearly parameter independent behavior has been achieved, which results in similar aircraft responses for the roll motion over the whole envelope. Roll rate is tracked up to a frequency of 8 rad s −1 and angle of sideslip up to 1 rad s −1 , above this frequency tracking accuracy diminishes due to higher prioritized yaw damping. The lateral controller has a number of 20 internal states. Figure 9 shows the closed loop poles of the controlled LPV system versus the poles of the uncontrolled LPV plant. According to the previously explained tracking behavior, it can be clearly seen that a unique roll time constant for the full envelope is generated. This ensures similar times for the build-up of roll rate independent of the dynamic pressure, thus the pilot can expect similar behavior of the aircraft's roll motion without considering the actual airspeed and altitude. Moreover, an increased damping of the DR motion is achieved when comparing the closed loop with open loop poles.
Furthermore, Fig. 10 depicts the singular values over frequency of the resulting LPV lateral inner loop controller. Up to the DR frequency all curves show similar, integrating behavior whereas above this frequency, until about 10 rad/s the dependency on the value of the scheduling parameter can be seen in the controller gain. Figure 11 depicts the phase of the minimal disk margin of the lateral controller is depicted. For the comparison of classical gain / phase margins with the robustness margin computations introduced in Sect. 3.2, the values shown in Table 3 are calculated. With a minimum of 4.2 dB, sufficient robustness even against simultaneous variations in all channels is ensured.

Controller validation
This section gives an overview of the means of control law validation, staring with SIL simulations, then moving to integration of the algorithms onboard the Cessna Citation research aircraft leading to flight testing.

Software in the loop verification
For verification of the designed control laws, an extensive simulation campaign has been performed, including nonlinear simulation on the DLR Robotic Motion Simulator (RMS). For this first controller verification step, a set of different maneuvers, e.g., step-like commands and augmented piloted flights, have been carried out. Additionally, changes in the operating point have been used to verify smooth transitions of the controlled aircraft between the test points. The Institute of System Dynamics and Control's RMS [3] has been used as a platform for flight test preparation via non-linear simulations. The RMS is a six degree of freedom industrial robot with an additional linear axis, as seen in Fig. 12. The simulator seats a single pilot and can be equipped with aircraft controls and virtual reality glasses in order to provide the user with a realtime response in terms of accelerations, body rates and attitudes. The pilot input, given via sidestick and thrust levers, is processed in the main simulation model and its 6 DOF outputs are delivered as motion command inputs to the RMS. Filtering algorithms ensure the gravity  Fig. 12. The SIL setup features a distributed computing environment, which consists of the above described RMS, pilot controls, the aircraft model described in Sect. 3 and detailed actuator models, which include the behavior of the aircraft's auto-pilot servos. The experimental FCS behavior, sensor and sensor bus environment with similar characteristics (transport delays, signal update rates) as the test aircraft were included. For convenience a visualization system was added. As a computing infrastructure, i7 CPUs running windows OS were used. The communication between the subsystems was implemented via network protocols and, in contrast to the flight test middleware onboard the test aircraft, no error handling of real-time violations was implemented.
The experimental flight control computer on the test aircraft is able to execute the control laws in continuous time, thus no explicit discretization of the LPV controller has been performed. All non-linear closed loop simulations as well as the SIL experiments on the RMS, however, were performed with a communication frequency of 100 Hz, the same frequency as used on the test aircraft. The SIL setup has been mainly used for purposes of • closed control loop simulations with pilot and automated control commands; • investigation of transient behavior during controller (de-) activation and operating point changes; • procedure training for flight test crew.
After successful completion of the SIL testing campaign, the control software has been considered ready to be integrated into the DUECA distributed real-time computation framework [15] which is used on the test aircraft. After this step, the DUECA framework manages communication between sensor, actuation signals and the FCS.

Test aircraft and control system
The Cessna Citation PH-LAB in Fig. 2, jointly operated by TU Delft and the Dutch Aerospace Center (NLR), serves as a multi-functional research platform. The aircraft is certified according to CS25 specifications [7] for large airplanes and equipped with a conventional, fully reversible flight control system providing a fix-geared link between the pilot's controls and the control surfaces of the aircraft. Additionally, an autopilot which has authority over the primary flight controls (elevator, aileron and rudder)   Fig. 12 Robotic motion simulator at the DLR Institute of system dynamics and control [3] is available. The test aircraft is equipped with an experimental FBW system [30], which uses the autopilot servos as control actuators. This system also has been explicitly certified under CS 25 specifications [7]. In addition, a flight test instrumentation system [14] including further sensors (angle of attack, sideslip etc.) is available for data acquisition and logging. The hardware setup is described in more detail in Refs. [8,14,30]. Figure 13 shows a block diagram of the flight control software, including the LPV baseline controller, implemented on the test aircraft. As shown in Fig. 13, the command module translates the stick signal to the reference signals needed for the LPV inner loop. The described control system setup has been tested in previous experiments on the Citation aircraft [8,11,20] and has been reused in the shown flight tests. Additionally, an attitude hold option via an outer loop can be activated to the lateral LPV inner loop control system. It additionally corrects for disturbances in roll attitude. As a result, it is possible to test a typical Airbus-like system behavior with attitude hold, or rate-command-only mode, which corresponds to a control philosophy in which the pilot needs to correct errors in attitude due to external disturbances.
The output signals of the LPV control feed the FBW servo control system [5] explained above. The servo system internally uses position feedback from the control surfaces to set the control surface deflections demanded by the LPV inner loop control laws.

Test cases and test execution
For in-flight validation a number of test cases were defined in order to verify all control design objectives. The flight test cases used for the test flights are described in Table 4.   Fig. 14). Both flights followed a similar scheme: after successful activation of the controller followed by careful maneuvers by the experimental pilot, the commands were slowly expanded to the full range limited by maximum control authority, including sweeps through the full range of the reference signal. An additional analysis of the flight test data is provided in [27].
To evaluate the tracking behavior and stability (rise time, settling time) the pilot was asked to command step inputs with small amplitudes followed by the same maneuver being repeated with automated input signals. Since the tested LPV method should provide robustness against model uncertainties, configuration changes were performed during flight to investigate stability during transition as well as responses after the configuration change.

Flight test data
A selection of test data is presented and discussed herein. Firstly, the responsiveness of longitudinal and lateral control channels is discussed, followed by the analysis of a simulated engine failure case as an example for robustness.
Longitudinal flight test results Due to the autopilot pitch servo actuator having limited authority, the load factor reference controller is validated during a bank maneuver. The load factor command response is displayed in Fig. 15. This setup allows load factor commands of multiple seconds to be fed without changing the aircraft's operating point in terms of airspeed and altitude. In Fig. 15, the pilot first banks to the left performing two steps and after approximately ten seconds of steady turning, steps in bank attitude to the right are performed. The left upper plot of Fig. 15 shows the associated load factor demand and resulting load factor for the maneuver. It can be seen, that the commanded load factor during the bank maneuver is tracked nicely by the longitudinal controller, keeping the flight path angle (left lower plot) at around zero degrees.
Lateral flight test results For verifying the lateral tracking performance, small step inputs are used. Figure 16 depicts an automatic step input given to the system with the respective aircraft and as control surface responses. The left diagram in Fig. 16 depicts roll and yaw rate as well as the reference signal for roll rate. The roll rate ( ̇ ) follows the step in expected PT-1 like behavior with a time constant of approximately 1 s as designed for the closed loop (see Fig. 9). Yaw rate ( ̇ ) (blue line) reacts without oscillations. Thus, the DR mode is sufficiently stabilized by the lateral MIMO control law. The right diagram shows the aileron and rudder deflections. As the lateral control law is designed for the multi-input multi-output problem with aileron and rudder as outputs, both control surfaces are deflected simultaneously. In a conventional control law, the roll input would develop an angle of sideslip first, which would then be corrected by the yaw controller. In this case, rudder deflection follows immediately after initiating the roll maneuver. Following, the developed sideslip has a maximum of 0.3 deg after 4 s of the roll maneuver.
Simulated single engine failure Figure 17 shows the response of a (simulated) single engine failure. In this scenario the pilot manually decreased the power setting of one engine and at the same time increased the power of the remaining one to keep up constant level flying with constant indicated airspeed. The bottom left plot shows the actual engine power setting N1 whilst the top plots depict the response in angle of sideslip and roll attitude. The desired aircraft behavior is zero angle of sideslip and a steady roll attitude. Both goals are achieved by simultaneous deflections of rudder and aileron, which is the expected reaction of the lateral MIMO controller. With no angle of sideslip, the aircraft continues straight and level flying despite a differential thrust setting of 40 % engine power on the left and 95 % power on the right engine.
In addition to testing the different maneuvers listed in Table 4, the augmented inner loop control law was also enabled in the initial approach phase where aircraft configuration (flaps extended, gear down), deviated from the clean configuration used for the LPV controller synthesis. As the pilot was capable of maneuvering without difficulties even in landing configuration, a qualitative evidence for robustness against changes in the aerodynamic characteristics of the aircraft has been provided.  Figure 18 shows the comparison between flight test data presented in Fig. 15 and simulation data produced with recorded flight test commands as inputs. For re-simulation, the aircraft model was trimmed at the airspeed and altitude gathered from test execution. The inputs to the inner loop LPV commands in pitch and roll axis were recorded during flight and these values are now fed to the closed loop simulation model. The resulting outputs are depicted in Fig. 18 where dashed lines represent simulation data and solid lines the flight test data. For the simulation case, only roll rate command recorded from the flight test is fed to the inner loop, the attitude hold option is not activated. This causes a slight drift of simulated roll angle as seen in the lower left plot, but with this approach it is ensured that only the relevant LPV inner loop is compared. The disturbance reduction provided by the outer loop thus is only activated in flight test case. Due to the increased roll attitude (starting at 20 s), a higher load factor is needed for a coordinated turn. As a result, load factor and elevator deflection in the simulation data are slightly increased compared to the flight test data. The pitch attitude and path angle show similar values and behavior. This indicates a good compliance of simulated and flight test behavior of the closed loop system.

Conclusions
In this work a Linear Parameter Varying based flight control system has been developed for a CS 25 certified aircraft. The obtained flight controller has been validated in an extensive flight test campaign, by means of a selection of pitch and roll maneuvers. The tested flight control algorithms have been proven to provide very good handling qualities according to the test pilot's feedback. This observation is also backed by the evaluation of test data. The gain scheduling feature of the controller, which is an intrinsic property due to the design process, has been verified in flight through maneuvering at different airspeeds. Finally, configuration changes were used for demonstration of satisfactory robustness of the control laws against unknown model changes.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.

Fig. 18
Simulation and flight test data of a load factor response during a roll maneuver at 95 kts VTAS. Command ( ) and response values ( ) are color coded, flight test data are solid and simulation data dashed curves