Development of a quasi-linear parameter varying model for a tiltrotor aircraft

The research presented in this paper focuses on the development of a quasi-Linear Parameter Varying (qLPV) model for the XV-15 tiltrotor aircraft. The specific category of qLPV modeling technique, known as the model stitching technique, is employed to model the time-varying dynamics of XV-15 tiltrotor aircraft over the entire flight envelope. In this modeling approach, discrete linear state-space models are interpolated through lookup tables as function of scheduling parameters with the implementation of nonlinear equations of motion. The XV-15 qLPV model is configured with four scheduling parameters: altitude, nacelle incidence angle, wing flap angle and velocity. Additionally, a computational complexity analysis is presented. In particular, computational sensitivity of qLPV models configured with lookup tables to number of states and number of scheduling parameters is demonstrated. This is done to show the feasibility of real-time implementation of qLPV models with increasing fidelity (number of states) and expanding dynamic flight envelope (number of scheduling parameters).


Introduction
Tiltrotor aircraft represent a promising solution to future civil transportation requirements [1,2] due to their broader flight envelope compared to conventional and compound helicopters. Vertical Takeoff and Landing (VTOL) capabilities of a helicopter combined with flight at relatively high cruise speeds and range like a fixed wing airplane, make tiltrotors an effective point-to-point fast means of transportation. Their multi-role capability not only makes these rotorcraft ideal for urban air mobility, but also for energy and medical services, and search and rescue missions.

3
To improve the design of future tiltrotor aircraft, it is essential that the handling qualities are assessed at the initial design phase through real-time piloted simulations and improved through the development of Stability and Control Augmentation Systems (SCAS). Furthermore, aeroelastic loads must also be predicted and alleviated through the design of Structural Load Alleviation (SLA) control. To perform such piloted simulations and for the purpose of Flight Control System (FCS) design, it is essential to develop a full flight-envelope high-fidelity flight dynamics model. The dynamics of a tiltrotor aircraft do not only change with the flight conditions but also with the aircraft configurations, and hence it is particularly challenging to develop a continuous flight dynamics model for a tiltrotor aircraft due to three distinct flight modes, namely helicopter, conversion, and airplane. A number of studies have presented various approaches to model the dynamics of tiltrotor aircraft. These models include detailed mathematical simulation models [3][4][5] for the purpose of real-time flight simulations, aeroelastic analytical [6][7][8] and numerical models [9,10], Linear Time Varying (LTV) model [11], and frequency domain models identified from flight test data [12][13][14]. The flight dynamics models listed in the foregoing sentence are either low-fidelity models or are only valid for discrete flight conditions and aircraft configurations.
A better approach to develop a continuous model for the varying dynamics of tiltrotor aircraft is to employ the model stitching technique [15] that comes under the category of quasi-Linear Parameter Varying (qLPV) modeling technique [16]. In this technique, linear state-space models at discrete flight conditions, or aircraft configurations, are stitched together to form a continuous and time-varying flight dynamics model that encompasses the entire flight envelope. The model is quasi-nonlinear in that the linear stability and control derivatives and the corresponding trim data are interpolated through the implementation of lookup tables as function of time-varying scheduling parameters, but nonlinear equations for rigid body motion and gravitational force equations are implemented.
The qLPV modeling technique has already been used for a wide array of aircraft. These include fixed wing airplane such as Cessna CJ1 [15], Boeing 747 [16], and Learjet-25 [17]; and helicopters such as Sikorsky UH-60 Black Hawk [15], Bell 206 [18], and Eurocopter EC 135 ACT-FHS [19]. With reference to modeling a tiltrotor aircraft using qLPV modeling technique, only a limited number of studies have been conducted. For the purpose of handling qualities analyses, a qLPV simulation model including 13 states (six-degree-of-freedom rigid body states and rotor flapping states) was developed for NASA' LCTR2 (Large Civil Tiltrotor, 2nd generation) using only two scheduling parameters: velocity V and nacelle incidence angle i , at hover through low speed by Lawrence et al. [20]. Similarly, a qLPV model containing 51 states (six-degrees-of-freedom rigid body states, rotor states, inflow states and nacelle dynamics) scheduled with three parameters: velocity V, nacelle incidence angle i , and altitude, was developed for a generic tiltrotor configuration by Berger et al. [21]. In both studies, engine-governor dynamics were missing, limited number of rotor elastic states were modeled and wing aeroelastic states were neglected, which play an important role in Rotorcraft Pilot Coupling (RPC) events [8,10].
The qLPV model can be utilized for the purpose of fullmission real-time piloted simulations and Flight Control System (FCS) design. The flight envelope of such simulations can be extended by increasing the range and the number of scheduling parameters in a qLPV model. Additionally, the number of states of linear state-space point models can be increased to improve the fidelity of a qLPV model. However, increasing the number of scheduling parameters will result in multivariate interpolation and the increase in number of states will lead to increase in the number of interpolations. In either case, the computational complexity will increase, limiting the real-time simulation capabilities.
The contribution and focus of the current paper is to develop a high-order aeroelastic qLPV model for tiltrotor aircraft representative of NASA's XV-15 with Advanced Technology Blades (ATB) [7], with four scheduling parameters: altitude h, nacelle incidence angle i , wing flap angle f , and velocity V. The high-order model includes engine-governor dynamics, rotor elastic states, and wing elastic states. The aeroelastic qLPV model for the tiltrotor aircraft is developed such that it can be used not only for the purpose of real-time piloted simulations and Flight Control System (FCS) design to meet the handling qualities requirement, but also for the prediction and subsequently, alleviation of aeroelastic loads at early design stage. Furthermore, a comprehensive computational complexity analysis of the qLPV modeling technique assembled with lookup tables is presented with the aim to investigate the effect of increasing number of scheduling parameters and states on the computational efficiency.
The paper is arranged into sections as follows: Theoretical background on LPV and qLPV modeling approaches is discussed in Sect. 2. Section 3 presents the local state-space models for XV-15 at discrete trim points. The qLPV model for XV-15 is detailed in Sect. 4. A thorough computational complexity analysis of the qLPV modeling approach is presented in Sect. 5. Last, a brief conclusion is presented in Sect. 6.

Theory
Linear state-space models that depend on time-varying scheduling parameters (t) are called Linear Parameter Varying (LPV). Formally, a continuous time LPV model is defined as follows [16,22]: and Δy(t) = y(t) − y trim ( (t)) are the state, control, and output perturbation vectors, respectively. The LPV model represented by Eq. 1 can be obtained through analytical methods (such as Jacobian linearization, state transformation and function substitution) or experimental methods or combination of both [16,22].
A particular case of LPV system is a quasi-Linear Parameter-Varying (qLPV) system, where a subset of scheduling parameter vector is also state of the system. Consider a state vector x(t) that can be decomposed into scheduling state vector z(t) ⊂ (t) and non-scheduling state vector w(t) ⊄ (t) ; then the state equation of a qLPV system can be written as follows: The scheduling parameter vector can be composed of both scheduling states z(t) and exogenous scheduling variables can be rewritten in terms of total states and total control inputs as follows: In the above equation, x trim ( (t)) = z T trim ( (t)) w T trim ( (t)) T and u trim ( (t)) are the vector of trim states and trim control inputs, respectively. Since z(t) is scheduling parameter and also state of the system, z(t) − z trim ( (t)) = 0 is true at all times. Therefore, all the stability derivatives with respect to z(t) are set to zero. However, the effect of scheduling states z(t) are preserved implicitly in the model by variation of the trim states and trim control inputs, such that the dynamic response of the model at anchor points is preserved [15,17,19,20].
(1) x(t) = A( (t))Δx(t) + B( (t))Δu(t) Δy(t) = C( (t))Δx(t) + D( (t))Δu(t) A further extension to the qLPV model, referred to as "model stitching technique", is proposed by Tobias and Tischler [15], where a set of linear state-space models and trim data at discrete trim points (also known as anchor points) are combined or stitched together to form a continuous full flightenvelope simulation model. The state-space entries of linear systems and the corresponding trim data, obtained at discrete anchor points, are interpolated through the implementation of lookup tables as function of instantaneous values of the scheduling parameters. Furthermore, the qLPV model is augmented with the rigid body nonlinear equations of motion and the nonlinear gravitational equations. This gives an advantage to the stitched model of being essentially nonlinear with linear time varying aerodynamics. In the current paper, this technique is employed to develop the flight dynamics model for XV-15. Figure 1 shows the model structure of high-order qLPV model implemented through lookup tables augmented with the nonlinear equations of motion, the so-called model stitching technique. First, the total number of states are divided into six rigid body states x 6 = u v w p q r T and n higher- Next, the rigid body mass matrix M composed of aircraft mass m and inertia tensor J is introduced as follows: Mass matrix M is multiplied with the rigid body accelerations ẋ 6 , obtained through Eq. 4, to obtain the perturbed aerodynamic forces ΔF aero and perturbed aerodynamic moments ΔM aero , since the linearized flight mechanics equation reads as follows: The perturbed aerodynamic forces are summed with the trim aerodynamic forces [23] of Eq. 7, to yield the following total aerodynamic forces: F aero = ΔF aero + F trim and total aerodynamic moments: M aero = ΔM aero .
Subsequently, the total aerodynamic forces are added to the nonlinear gravitational forces [24], shown in Eq. 8, to obtain the following total forces: F total = F aero + F grav and following total moments: M total = M aero .
Notice that, in order to determine the aerodynamic trim forces in Eq. 7, the trim Euler angles are required that are obtained through lookup table interpolation of the trim data. On the other hand, for the computation of gravitational forces in Eq. 8, no lookup of trim Euler angles is performed, but rather the current instantaneous values of the aircraft Euler angles are used.
In order to obtain the rigid body state derivatives ẋ 9 (composed of six rigid body accelerations ẋ 6 and three Euler angle rates), the total forces and moments are passed through the six-degree-of-freedom body-axes nonlinear equations of motion [24]. As a last step in the qLPV model structure, the rigid body state derivatives ẋ 9 along with the higher-order state derivatives ẋ H , obtained from Eq. 4, are integrated to obtain the current states.
It should be noted that, in this particular technique, i.e. qLPV model augmented with the nonlinear equations of motion, the state matrix A 66 contains only the stability X grav = −mg sin Y grav = mg cos sin Z grav = mg cos cos derivatives; it does not contain the gravity, the Coriolis and the kinematic terms, and it does not include the Euler angle states [ ] . The effect of gravity is later incorporated in the nonlinear gravitational force equations, and the Coriolis effect and the kinematics are added within the nonlinear equations of motion, represented by the yellow blocks in Fig. 1.

XV-15 state-space point models
Linear state-space models at discrete trim points serve as the building blocks for developing a continuous qLPV model. In the current work, simulation tool MASST (Modern Aeroservoelastic State Space Tools), developed at Politecnico di Milano [25,26], is employed to generate the local statespace models of XV-15 using the data published in Refs. [4,7]. MASST implements the rotor aeroelastic models obtained from CAMRAD/JA [27] and the flexible airframe is implemented from aeroelastic NASTRAN model.
Local state-space models obtained through MASST contain 91 states and 11 control inputs, and are listed in Table 1. The model includes the flexible rotor blades with three coupled flap/lag bending modes, two torsional modes (one rigid pitch and one elastic torsion), and one gimbal mode. Furthermore, the model also includes the symmetrical wing bending mode that is the main source of wing-pilot vertical bounce RPC event [8,10]. Figure 2 shows a zoomed-out and zoomed-in view of the eigenvalues of the full-order linear state-space models in helicopter mode at 0 kts (hover) and in airplane mode at 180 kts at sea level. For two rotors, there are two of each rotor modes which are labeled by their dominant motion. All the rotor modes, except for regressive gimbal and regressive lag, are placed above 10 Hz, far away from RPC frequency range. However, these modes are important to accurately model and perhaps alleviate the critical aeroelastic loads [29]. Symmetrical wing bending, labeled as SWB in the zoomed-in view, is placed at −0.67 ± 15.26j in hover. Lastly, the rigid body modes are listed in Table 2. Similar to a conventional helicopter, the longitudinal phugoid and dutch roll modes are unstable in helicopter mode in hover. As the aircraft transitions from helicopter to airplane mode and the speed of the aircraft increases, these modes become stable. Furthermore, the pitch and heave subsidences in helicopter mode are placed on the real axis at −0.64 and −0.26 , respectively. As the speed increases, these modes combine to form shortperiod mode.

Choice of scheduling parameters
In order to develop a qLPV simulation model, a minimal set of scheduling parameters need to be defined. It is essential to choose only the key scheduling parameters that affect the system dynamics significantly, because the computational complexity of implementing a qLPV model, as will be shown later in Sect. 5, and the numerical complexity of subsequent LPV control synthesis increases rapidly with the increasing number of scheduling parameters [30]. In the current study, qLPV model for XV-15 is developed using four scheduling parameters: altitude h, nacelle incidence angle i , wing flap angle f , and velocity V. Velocity and nacelle incidence angle are chosen as scheduling parameters for qLPV model of tiltrotor aircraft to accurately capture the changing dynamics, such that the entire conversion corridor of tiltrotor aircraft is encompassed.
Wing flap serves a dual purpose, i.e., functioning as a high lift device as well as reducing the wing download at hover and low speeds. Wing flap position not only affects the lower limit of conversion corridor of tiltrotor aircraft by delaying the wing stall [31], but also affects the dynamic behavior of the aircraft. The effect of wing flap setting on the dynamics of XV-15 is shown in Fig. 3, where the stability derivative X u is plotted as function of velocity at four wing flap settings in helicopter mode, Fig. 3a, and in airplane mode, Fig. 3b. It is evident from the figure that the stability derivative in both configurations have different magnitudes at different wing flap settings. Therefore, for the sake of completeness and in order to explore the edges of conversion corridor, it is essential to add wing flap angle as a scheduling parameter. Furthermore, this will allow for simulating the failure scenarios, where the wing flap is stuck in one position.
Lastly, altitude is added as scheduling parameter to ensure continuous and smooth change in altitude during simulations. It is suggested in Ref. [17] that the variation in trim and dynamics with altitude can be accounted for by applying air density ratio scaling, without explicitly adding altitude as a scheduling parameter. This is true for fixed wing aircraft, however, for rotorcraft not all the derivatives scale correctly with air density. Therefore, altitude is also added as an additional scheduling parameter.

Anchor points
A collection of linear state-space models and trim data is generated on a four-dimensional grid of scheduling parameters. Models and the corresponding trim data are obtained by trimming the aircraft in symmetric flight ( v = p = r = = = 0 ) at discrete anchor points spanning the entire conversion corridor. Furthermore, each model is obtained at two altitudes h = [0 10000] ft and four wing flap settings f = [0 20 40 75] deg. Figure 4 shows the 2-D grid (velocity-nacelle incidence angle) of linear state-space models computed for each wing flap setting and altitude, along with the conversion corridor of XV-15 [32]. The figure also shows the grid points that are outside the conversion corridor where the models and the corresponding trim data are not available. The models at these grid points are obtained by clipping and keeping the edge models of either lowest or highest speed at a given nacelle incidence angle. This is done to construct a regular rectangular grid, which is a prerequisite to develop a model stitching simulation architecture [15].
To obtain the state-space models and trim data, grid spacing of 5 knots between 0 and 50 knots, 10 knots between 60 and 170 knots and 20 knots between 180 and 280 knots was deemed satisfactory. Finer grid spacing of 5 knots was chosen near hover to accurately capture the dynamics of the aircraft for the simulation of low-speed flight [15]. Similarly, lookup table implementation of the trim data and stability and control derivatives with models at only two altitudes, sea level and 10,000 ft (FL100), is sufficient to yield acceptable results. As an example, in Fig. 5, trim pitch attitude is shown as function of velocity at varying altitudes. The interpolated and extrapolated trim pitch attitudes at 6000 ft (FL60) and 15,000 ft (FL150), respectively, show good agreement with the original trim pitch angle obtained by MASST. The agreement between original pitch attitude and interpolated/extrapolated pitch attitude is consistent in all configurations, i.e., helicopter mode, Fig. 5a, conversion mode, Fig. 5b, c, and airplane mode, Fig. 5d.
Trim data and the elements of state and control matrices are interpolated offline by means of spline interpolation to obtain a finer grid, with a grid spacing of ΔV = 5 knots and Δ i = 5 • . This data processing resulted in 8664 linear state-space models and corresponding trim data, spread over a four-dimensional grid of scheduling parameter vec-

Trim data
Trim surfaces as function of nacelle incidence angle and velocity at sea level and wing flap angle f = 20 • are presented in Fig. 6, featuring pitch attitude of aircraft, gimbal pitch for right rotor, elevator deflection, and collective pitch for right rotor. Values at trim points, where the models are obtained through MASST, are also depicted in the figure by black points. Aircraft trim pitch attitude, Fig. 6a, decreases with the increasing speed at all nacelle incidence angles. At a given speed, the trim pitch attitude increases with the decreasing nacelle incidence angle because as the aircraft transitions from helicopter to airplane mode, more lift is generated by the wings compared to the rotors. The trim rotor gimbal pitch (right rotor) is depicted in Fig. 6b; it decreases (flaps backward) with nacelle incidence angle from helicopter mode through to conversion mode. However, in airplane mode the trim gimbal pitch increases and grows further (flaps forward) with the increasing speed to compensate for the decreasing pitch attitude of the aircraft. The trim elevator deflection, Fig. 6c, keeps on increasing with the increasing speed (trailing edge down) to trim the aircraft at increasing pitch down attitude. Lastly, trim collective pitch (right rotor) is presented in Fig. 6d. As expected, the trim rotor collective pitch at high nacelle incidence angles ( i ≥ 75 • ) follows the trend similar to that of a conventional helicopter. The trim rotor collective pitch decreases with the increasing speed until the minimum power required, and subsequently it starts to grow to overcome the drag. At lower nacelle incidence angles ( i ≤ 60 • ), the trim rotor collective pitch keeps the increasing speed of the aircraft constant by contributing towards the generation of thrust.

Quasi-linear parameter varying model for XV-15
Referring to Fig. 1  should be pointed out that, in the current study only the symmetric flight conditions are considered, therefore, altitude h and velocity V only depend on the longitudinal states. Lookup tables are implemented to interpolate the trim data and the state-space matrices as function of scheduling parameters to ensure smooth transition of the system dynamics. It should be noted that, in the case of statespace matrices, interpolation is based on low-pass filtered velocity, i.e., the scheduling parameter vector for stability and control derivatives is where V filtered is defined as follows: with a cutoff frequency of c = 0.2 rad/s [15]. This is done to ensure constant derivatives for short term motion, such that the accurate dynamic response at a local trim point is retained. Lookup tables of trim states, trim controls and state-space matrices are implemented and nonlinear gravitational force equations and nonlinear rigid body equations of motion are added to develop the complete qLPV simulation model for XV-15, as explained in Sect. 2 and depicted in Fig. 1.
It is important to note that the control derivatives associated with the wing flap angle f are set to zero in control matrix B , because f is one of the scheduling parameters. The effect of change in wing flap angle is preserved implicitly in the model, as mentioned previously in Sect. 2. Furthermore, unlike the f derivative, the u and w derivatives are not zeroed out because the body x−axis velocity u and the body z−axis velocity w are not considered as the scheduling parameters explicitly, but rather as a combination in the form of total velocity V.

Additional simulation elements
Additional simulation elements are incorporated within the qLPV model for XV-15. Details of these elements are discussed in the following sub-sections.

Rotor speed governor
The main purpose of rotor speed governor is to maintain a constant rotor speed. Figure 7 shows the block diagram of rotor speed governor that is based on the blade pitch governing scheme, the so called beta-governor [33]. The governor is a proportional-integral (PI) controller that operates on the rotor speed error feedback and outputs the desired changes in the blade collective pitch. The PI gains of the rotor speed governor are shown in Table 3. In helicopter mode, the governor pitch angle gov is added to the blade collective pitch, commanded by the pilot through the collective stick X col through the gearing ratio K 0 i . As the aircraft transitions from helicopter mode to airplane mode, the blade collective pitch from the collective stick is gradually phased out [3,4]. In airplane mode, the collective stick only controls the throttle t through the gearing ratio K t and the blade collective pitch is exclusively regulated by the governor. The reference rotor speed is Ω ref = 62.93 rad/s until 160 knots and at higher speeds the rotor speed is maintained at Ω ref = 50.35 rad/s.

Actuator dynamics
The qLPV model is incorporated with a first-order actuator dynamics model, Eq. 10. Table 4 lists the time constants and saturation limits [4,34,35] for control inputs. Figure 8 shows a frequency response comparison of the XV-15 roll rate, pitch rate, and yaw rate to lateral stick, longitudinal stick and pedal input, respectively, in helicopter mode in hover. The frequency response comparison is shown between the high-order linear state-space point model and the qLPV model developed in the current study. The response is measured in deg/s, while the input is in inches of stick deflection. The qLPV model is incorporated with a first-order actuator dynamics as presented in Sect. 4.1.2, therefore, the same actuator dynamics is also added to the linear state-space point model. The qLPV model and the linear state-space point model compare extremely well with each other in all axes, showing that the grid spacing of trim data is implemented correctly in the qLPV model. The roll rate and the pitch rate response are dominated by the lowfrequency phugoid pole.

Frequency response analysis
The comparison of the primary on-axis frequency response in airplane mode at 170 kts is shown in Fig. 9. The frequency response of the roll rate, pitch rate, and yaw rate in deg/s is shown to aileron, elevator, and rudder deflection

Time response analysis
The qLPV simulation model developed for XV-15 using the model stitching architecture is validated by comparing the SCAS off time responses with the data presented in the existing literature. Figures 10 and 11 present the time history correlation with NASA's Generic Tilt Rotor Simulation (GTRS) model [5] in helicopter mode and airplane mode, respectively. In Fig. 10, pitch and yaw time histories in helicopter mode are shown. Figure 11 shows the comparison of pitch, roll, and yaw time histories between qLPV model and the GTRS model in airplane mode at 175 kts. Lastly, Fig. 12 shows the comparison with the XV-15 simulation model developed in FLIGHTLAB ® [36] in conversion mode at nacelle incidence angle of i = 60 • and 120 kts. Table 5  Except for yaw rate r in conversion mode, the RMSE is less than 2 deg/s (or deg for Euler angles) for all the three axes and aircraft modes, which reflects an acceptable level of accuracy, as suggested in Ref. [23]. The yaw rate response of qLPV model is slightly slower than the FLIGHTLAB Â® model, as can be seen in Fig. 12c. This shift in the time response is the reason for a slightly higher yaw rate RMSE in conversion mode.

Computational complexity
In order to use qLPV models for the purpose of fullmission real-time piloted simulations, it is necessary to identify performance bottlenecks by examining the computational complexity of such models and provide recommendations for improvement. The theory of computational complexity is concerned with the required computational resources to solve a given task in terms of space and time [37]. In the current paper, focus of the complexity analysis is on the computational efficiency of the algorithm in terms of time complexity. The efficiency of an algorithm is quantified by the number of basic operations, or the primitive operations, as function of input size [37]. The time to execute these basic operations (e.g., addition, multiplication etc.) is constant and does not depend on the input size. Therefore, total time of an algorithm is proportional to the number of primitive operations [38].

Complexity analysis of the qLPV modeling technique
Implementation of qLPV model that is augmented with the nonlinear equations of motion and is configured with the  (c) Response to pedal input Therefore, the total number of primitive operations for implementing lookup tables for n x trim states and n u trim control inputs is N 1 n , n x , n u = n x + n u � n ∑ i=1 log 2 k i + 3(2 n − 1) � .
itive operations, respectively. 3. The number of primitive operations to compute state and control perturbations is N 4 n x , n u = n x + n u . 4. Considering the mathematical definition of matrix multiplication, the modified state equation, Eq. 4, require N 5 n x , n u = 2n x 2 − 15n x − 6n u + 2n x n u + 31 primitive operations. 5. The number of primitive operations for nonlinear gravitational force equations and nonlinear equations of motion is constant and is denoted by N eom . 6. The number of primitive operations for integrating the state derivatives N int depend on the integration scheme. Integrating all the state derivatives require N 6 n x = n x N int primitive operations.
The number of time steps N t depend on the type of the solver, either variable-step or fixed-step. Furthermore, in case of continuous time implementation, the solver might subdivide the time step into minor steps. Following the discussion above, the total number of primitive operations required to implement the qLPV modeling technique configured with lookup tables is as follows: The number of primitive operations, and hence the computational efficiency, of qLPV model clearly depend on the number of states and the number of scheduling parameters. It should be noted that, the increase in number of states means improvement in fidelity and that the additional degrees of freedom can be incorporated into the model. On the other hand, increase in the number of scheduling parameters means that a broader dynamic flight envelope can be explored through the implementation of qLPV model. Figure 13 shows the total number of primitive operations N p n , n x , n u required to implement a qLPV model in semilog scale as function of number of scheduling parameters n with different number of states n x . Number of primitive operations N p n , n x , n u are approximated by assuming n u = 10 , N eom = 50 , N int = 4 , and N t = 1773 . Similarly, Fig. 14 shows the simulation time in semi-log scale to run a 10-second simulation of XV-15 qLPV model for three different numbers of scheduling parameters with three different number of states. Scheduling parameter vector for n = 2 is: (t) = i V , n = 3 is: (t) = h i V and n = 4 (11) N p n , n x , n u = N 1 n , n x , n u + N 2 n , n x +N 3 n , n x , n u + N 4 n x , n u +N 5 n x , n u + N 6 n x + N eom N t is: The simulation is implemented in Matlab Â® /S imulink Â® R2019a on an Intel Â® Xeon Â® , 3.70 GHz processor running Windows Â® 10 with simulation parameters listed in Table 6. The simulation step size is chosen to be sufficiently small enough to accurately capture the dynamics of the aircraft, while at the same time large enough to ensure real-time implementation.
It can be seen from Fig. 14 that the simulation time is less than 10 s, guaranteeing real-time implementation; however, further increase in the number of scheduling parameters may not allow the execution of XV-15 qLPV model in real time. Figures 13 and 14 show that the simulation time is proportional to the number of primitive operations. Furthermore, these figures show that the increase in number of scheduling parameters to extend the flight envelope or increase in number of states to improve the fidelity of model will result in increase in the number of primitive operations and hence, the simulation time. Following the discussion above, it can be concluded that the computational complexity of a qLPV model implemented through lookup tables is O N t ⋅ n x 2 ⋅ 2 n .

Conclusion
A high-order full flight-envelope quasi-Linear Parameter Varying (qLPV) simulation model for NASA's XV-15 with Advanced Technology Blades (ATB) tiltrotor aircraft using the model stitching technique is developed and presented in this paper. The varying dynamics of a tiltrotor aircraft with flight conditions as well as aircraft configurations make it particularly hard to develop a continuous full flight-envelope flight dynamics model for the purpose of real-time piloted simulations and control synthesis. Therefore, a quasi-Linear Parameter Varying (qLPV) modeling approach is employed. The qLPV model for XV-15 is developed by stitching together the discrete linear state-space models and adding six-degree-of-freedom body-axes nonlinear equations of motion and nonlinear gravitational force equations to build a continuous and quasi-nonlinear simulation model. The statespace models are obtained on a grid of scheduling parameters: altitude, nacelle incidence angle, wing flap deflection, and velocity. It is shown that the linear state-space models at only two altitudes, sea level and 10000 ft, are sufficient to yield acceptable results through lookup table implementation with linear interpolation and extrapolation.
Furthermore, additional simulation modules, such as rotor speed governor and actuator dynamics, are added to the qLPV model of XV-15. Rotor speed governor is a proportional-integral (PI) controller that regulates the blade collective pitch to maintain a constant rotor speed.
The computational complexity of qLPV model is quantified in terms of total number of primitive operations. It is shown that the computational complexity of qLPV model mainly depends on the number of parameters and number of states. Number of parameters can be increased to expand the dynamic flight envelope of a qLPV model; however, this means that the lookup tables with multivariate interpolation are required. On the other hand, fidelity can be improved by increasing the number of states, but this will lead to increase in the number of interpolations. Overall, it was concluded that the computational complexity of a qLPV model is O N t ⋅ n x 2 ⋅ 2 n . In order to ensure the real-time implementation of qLPV model, a compromise has to be made between the number of states (degrees of freedom & fidelity) and number of parameters (flight envelope coverage).
In the literature thus far and the current paper, trim data and linear state-space models are linearly interpolated using lookup tables. However, the fundamental limitation of lookup tables is the lack of continuity and adaptability. Updating such models is a complex and time consuming task. Therefore, the future work will extend to developing qLPV models based on multivariate simplex spline models. Multivariate simplex splines are function approximators that are capable of fitting highly nonlinear datasets and naturally capable of fitting scattered datasets on non-rectangular domains. They have the advantage of being continuous analytical functions that can allow local model modification without affecting the global model structure. Furthermore, qLPV models based on multivariate simplex splines will assist in the development of fault-tolerant Flight Control Systems (FCS), especially in the presence of local and global nonlinearities.
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/.