1. Introduction
Systematic flight testing of remotely piloted aircraft (RPA) is usually done in two different situations: One is the development of the platform itself for its use as a final product, as in the case of unmanned aerial systems (UASs). The second is during the early investigation and development of a particular configuration, feature or technology that is intended to be part of a larger, often manned, system. In this second case, typically referred to as scaled or subscale flight testing (SFT), RPAs are carefully designed to reproduce certain characteristics of the full-scale system [
1,
2]. During the early stages of design or technology development, SFT can be used to increase the available knowledge in a relevant environment without the high cost and risk associated with full-scale flight testing [
2,
3,
4]. An example of an early use of this technique is shown in
Figure 1.
One of the potential areas of interest for SFT is the investigation of flight dynamics. The complexity of modern aircraft design requires both the static and dynamic flight-mechanical characteristics of the aircraft to get a full picture of the expected performance. Historically, wind-tunnel testing has been the main source of aerodynamic knowledge for aircraft design [
5]. Currently, the static aircraft characteristics can be evaluated more easily and earlier in the design thanks to the continuous increase of computational power and the advances in computational fluid dynamic (CFD) methods. However, reliable dynamic characteristics are still difficult to obtain and evaluate or they become only available at later stages of detailed design. This makes SFT an interesting and affordable method to reduce uncertainty: the stability and control characteristics of the full-scale aircraft can be assessed with a dynamically scaled model, as described in [
1,
2,
6,
7].
Flight testing a remotely piloted aircraft can in some cases require a costly infrastructure with robust and reliable ground support for the flight crew and the test team, as described in [
8]. Fortunately, modern developments in consumer radio-controlled (RC) systems and UAS equipment have enabled smaller organizations to perform relatively advanced flight tests with less resources and in shorter times [
9]. Due to the current regulations in most countries and unless a segregated airspace and the corresponding clearances can be afforded, the flight test is likely to take place at low altitude and within the visual line-of-sight (VLOS). While this type of operation is more economical, it presents two main challenges represented in
Figure 2: First, subscale demonstrators often have a high wing loading and fly faster than typical RC aircraft, leaving little time for test manoeuvres on each pass. Second, the aircraft is exposed to higher levels of turbulence near the ground. An illustrative example of the consequences of VLOS limitations for subscale experiment design is the recent FLEXOP project: even the test-bed aircraft was specifically designed to cope with the fast transitions needed to reach the desired flight condition within the designated test perimeter [
10].
The methods presented in this paper aim to improve both the productivity of the flight test and the quality of the modelling under these demanding circumstances. The focus is mainly on the initial flight envelope expansion, where usually, linear flight-mechanical characteristics are studied and primary flying qualities are compared to the existing models. At this stage, test manoeuvres typically involve relatively small disturbances from trimmed states. Single, double (doublets) or multiple pulses are some of the control inputs traditionally used for excitation in both open- and closed-loop systems.
Here, the application of multisine input signals for system excitation is proposed. Compared to classical pulses, multisine input signals enable simultaneous excitations of different control surfaces, and hence, they have the potential to maximise the information obtained during a very short time. Although uncommon in full-scale aircraft testing, this type of input signal was tested on NASA’s T-2 subscale jet transport aircraft in 2012 [
11]. While NASA’s T-2 is comparable in terms of size, flight speed and instrumentation to the aircraft tested here, the testing environment is relatively more demanding here due to VLOS operation.
The means for dealing with air turbulence, especially noticeable near the ground, also need to be considered. There are aircraft parameter estimation algorithms for system identification that take process and measurement noise into account, as described in [
12]. However, the problem formulation in this paper is of an errors-in-variable character. An instrumental variable (IV) approach is a common way of dealing with this type of problem. An algorithm using this approach was given in [
13]. This approach is used to improve the algorithm performance.
An additional issue when using SFT at an early stage of design is that only basic knowledge of the aircraft characteristics may be available by the start of the campaign. It can therefore be helpful to use online tools for monitoring the information acquired during flight in real time. One such interesting method for the estimation of the flight mechanical parameters was given by Morelli [
14]. This method is based on a frequency domain approach, which is fast and suitable for real-time applications. It can also be used to make decisions about the continuation or interruption of the testing if anything deviates from the expected conditions.
The methods proposed here are described in
Section 2. These were also tested on a real flight experiment using a jet-powered subscale aircraft denominated the Generic Future Fighter (GFF) demonstrator [
15]. The multisine input signals were designed by using a simulation model based on basic data from earlier flights [
16]. The results of the input design are presented in
Section 3, and the details of the experiment and its results are presented in
Section 4,
Section 5 and
Section 6. Conclusions are presented in
Section 7.
3. Input Design Using Simulated Data
This section describes the design of input signals applied during the practical experiment on the Generic Future Fighter (GFF) subscale aircraft. First, a linear state-space simulation model was created based on the data available from previous flights of this aircraft [
16]. This model takes the form:
where the considered variables are defined in
Figure 3.
The states are the angle-of-attack and the pitch rate, and the input consists of the elevator and canard deflections. All states are measured.
Making an estimation of the model parameters
Z and
M in (
16) using an ordinary least squares (OLS) time domain method [
17] on the existing data shown in
Figure 4 leads to the following continuous-time model:
The OLS in the time domain is used here in batch mode since the whole dataset is available, and it is a quick and easy way to get a model to use. The model (
17) was validated on other datasets. One of these is shown in
Figure 5. Even though the match was not perfect, the model was deemed accurate enough to be used for the experimental design.
Two types of signals were evaluated (see
Figure 6). The first one is a double pulse, commonly used for estimating aerodynamic derivatives from flight test data. The second one is a multisine signal:
where
and
are the amplitudes and phases of the signal components, respectively. One reason for using a multisine is that it is possible to achieve weakly correlated, simultaneously excited input signals by using different frequencies for the different control surfaces, a method described in [
11]. This will save time and reduce cost since it reduces the time spent in the air for the same amount of information gathered.
Nine different crest factor multisine input signals were created with a sample frequency of
Hz and a period of
s. The toolbox described in [
18] was used to generate the crest factor optimized input signals. An example of the amplitudes and the phases of the separated input signals is shown in
Figure 7, where the phases corresponding to low-amplitude frequencies are omitted.
Three periods of excitation were used to get good input signals: one to take care of the transient in the beginning and two steady-state periods. This type of signal was described in [
19], where it was used to detect non-linearities.
Note that the frequency span of 1–10 Hz was chosen to excite the short period flight-mechanical model for the GFF subscale demonstrator. This is higher than the 0.1–2 Hz typically used for full-sized aircraft as in [
14]. It is a well known fact that small-scale aircraft have this short period characteristic. This phenomenon can be seen in the Bode plot of
Figure 8, where the model of a full-scale Saab JAS 39 Gripen used in [
13] is compared to the subscale GFF model from (
17).
To investigate the possible benefits of the multisine signals, the inputs for the elevator and the canard will be run both separately and in parallel. They will also be compared to the double pulse input when separated and with a small time delay to mimic the parallel input without causing too much correlation between the different control surfaces. Examples of input signals from simulations, with a process noise of about ten percent of the input amplitude, are shown in
Figure 9. The total excitation time here is
s, which is the expected time window available in the real flight test.
Making simulations with the estimated models for nine realisations of the input signals and comparing the double pulse (DP) and the crest factor optimized multisine (MS) results to the true model using the relative error:
as a measure give the results in
Table 1. It is clear that the crest factor optimized input gives more accurate estimates for this type of testing.
4. Experimental Setup
Due to its size and its instrumentation, the GFF subscale demonstrator was deemed a suitable platform to evaluate the usefulness of these methods. The GFF was a research study carried out during the NFFP4 Swedish National Aeronautical Research Program, and although it originated from a conceptual study at Saab, it is a generic design with no couplings to any Saab products. The GFF subscale demonstrator, shown in
Figure 10, was originally built by Linköping University during that project and made its first flight in 2009 [
15]. It is a radio-controlled jet-powered fighter with a close-coupled delta-wing canard configuration that reproduces the geometry of the full-scale GFF design at a
scale. Dynamical scaling, using the Froude number as the similarity parameter [
1,
2], is achieved by adding extra weights at specific locations. However, the demonstrator is currently flown without extra ballast at a lower take-off weight in order to simplify its certification and operation in experiments that do not require extrapolation of results to the full-scale aircraft. Such was the case here, where the experiment was aimed at evaluating the method rather than to scale-up the obtained coefficients. The main characteristics of the subscale demonstrator at the moment of testing are given in
Table 2.
The GFF demonstrator is one of the platforms used in the MSDEMO and MESTA research programs for investigating the potential of subscale flight testing within aircraft conceptual design [
9]. It is equipped with a data acquisition system and several sensors: Two inertial measurement units (IMU) provide accelerations and angular velocities. The IMUs together with redundant GPS receivers, magnetometers and air-pressure transducers are used to estimate the position, orientation and altitude. A custom-made nose-boom provides the airspeed, angle-of-attack and angle-of-sideslip using a pitot tube and contactless flow-angle transducers. The system also registers the incoming pilot commands and measures the real deflection of each control surface. These signals are logged on-board at an average sampling frequency of 100 Hz and also sent to a ground station in real time via a separate telemetry link.
Figure 11 shows the layout of the data acquisition system.
While the GFF design represents a modern fighter aircraft with an augmented flight control system, the GFF subscale demonstrator was flown with open-loop control in a stable configuration to reduce risk and simplify its certification. After all, the methods studied here are appropriate for both stable and unstable systems and no significant information was omitted using this simplification. Furthermore, this experiment was limited to the evaluation of linear characteristics expected during the initial expansion of the flight envelope. The study of non-linear characteristics, typically found in these fighter configurations at the edges of the flight envelope, is usually a tedious and very specialized task that falls beyond the scope of this experiment.
Flight Test Procedure
The flight test campaign was carried out at a closed airfield located inside the controlled airspace of a nearby airport. Hence, all operations were carried out within VLOS and in contact with the corresponding air traffic control service. This type of operation severely constrains the available flight area, as shown in
Figure 2; however, its main advantage is that it does not require a segregated airspace, a large infrastructure or a large team.
Figure 12 describes the test procedure including the three main roles: one test conductor, responsible for controlling the test execution; one pilot, flying the aircraft from the ground control; and one test monitor, evaluating the test results between the flights and supervising the safety area. Test cards were used to predefine the manoeuvres to be performed in each flight. The corresponding commands and excitation signals were pre-programmed in advance and automatically executed from the radio-control transmitter by means of a novel flight test software program developed at Linköping University and described in [
9]. This program augments the pilot control without requiring a closed-loop flight control system on-board. On the field, the easy adjustment and reconfiguration of the test sequences between flights also contributed to the agility and efficiency of this campaign.
The flight test was done during the evening to minimize the exposure to air turbulence; see
Figure 13. Three flights of approximately 10 min of duration were completed. The first flight was dedicated to system checks and instrument calibration. The other two flights were dedicated to collect identification data. The first two of the nine multisine input signals, mentioned earlier in
Section 3, were used both for parallel and separate excitation of the elevator and canard. The test program can be seen in
Table 3. The test setup was used to get repetitions of excitation for different fuel weights. The fuel weight is about
of the take-off weight, and it can, therefore, have an effect on the results.
5. Results
Some examples of estimation results from the four different types of manoeuvres (multisine 1 separated (MSs), multisine parallel (MSp), double pulse parallel (DPp), double pulse separated (DPs)) defined in
Table 3 are shown in
Figure 14,
Figure 15,
Figure 16 and
Figure 17. The test with the separated double pulse (DPs) had to be repeated twice during flight FT.
since the first one was aborted due to ground proximity. This gave two realisations, DPs1 (aborted) and DPs2 (repeated). The repeated double pulse was used for the estimations. The figures show the flight test data in red and the model (
17), based on earlier data, in a black dashed line. The estimated parameters are given as blue squares for the original method and as black circles for the new method. The estimated uncertainty measures of two standard deviations are shown as vertical bars. It should be noted that the model (
17) is not the true system, but rather a reference that can be used to say if the flight test results show control effectiveness or stability characteristics that are different from the model. If these differences were large, they could motivate the decision to abort the test until the causes are investigated and understood.
When comparing the results in
Figure 14 and
Figure 15, the methods give quite similar estimates in the end for the separated and parallel multisine inputs. This can be seen in
Table 4. The difference is during the test manoeuvres where the estimated canard effectivenesses
and
are quite uncertain before the excitation for the separated manoeuvres, as should be expected.
The differences are bigger when comparing the parallel multisine and double pulse in
Figure 15 and
Figure 16, predominantly in the pitching moment. The biggest difference is in the canard effectiveness
where the estimated effectiveness is much lower for the double pulse input. A similar result can be seen when comparing the two double pulse inputs in
Figure 16 and
Figure 17. This can also be seen in
Table 5.
The results from the repeated separated double pulse (DPs2) are shown in
Figure 17. It gives not exact, but similar results to the separated multisine in
Figure 14 when using the new method. It can also be seen that the difference between the estimates of the original method and the new method are largest for this manoeuvre.
In
Table 6, the absolute values of correlation between the real part of the regressors are given. The multisine input signals have a much lower correlation than the double pulses. The parallel double pulse, DPp, has the highest correlation of
. This might reduce the resulting estimation accuracy.
The validation of the estimated models is an important part of system identification [
17,
20]. The flight manoeuvre FT.
, DPs1, which was aborted due to a ground proximity warning, was used as the test dataset to validate the different estimated models. A common measure for judging how well a model,
m, can explain the information in the validation dataset,
Z, is the model fit criterion defined as:
where
indicates the two-norm of a vector.
This was used to compare the original and new method. The results are shown in
Table 7 and
Figure 18. Two things can be noted. First, the new method gives in general a better model fit than the original method. This means that the addition of instrumental variables contributes in a positive way. Second, the double pulse excitation gives a lower model fit than the multisine input with one exception. This may indicate that it is beneficial to use multisine inputs for this type of subscale flight testing.
The results shown in
Table 7 show an unexpected poor fit for manoeuvre MSs2 in Flight FT.
. The flight test data reveal that the pilot commanded a sudden roll correction during the execution of this manoeuvre, as shown in
Figure 19. The disturbance caused by this correction may explain the poor fit associated with this manoeuvre.
6. Post Flight Analysis
In post flight analysis, all test manoeuvres are available and can be used together. Doing this will improve the estimation result since the information content available is much larger than for each individual manoeuvre. The aborted manoeuvre, DPs1 in Flight FT.
, again was used as the test data for validation while all 14 remaining manoeuvres given in
Table 3, including DPs2 in Flight FT.
, were used for model construction.
A data fusion approach, like the one given in [
21], can be used to combine models from different manoeuvres into one single model. The model fusion is a weighted average of estimates
, where the weights are the information matrices
from different time sequences. The model fusion is therefore given as:
For the original method, the resulting model becomes:
while for the new method, the resulting model becomes:
The validation of these models using the DPs1 dataset produces the results shown in
Figure 20, where the initial model (
17) is also included for reference. The resulting models for both the original and the new method are similar, and they give a satisfactory fit of
and
, respectively. The initial, pre-flight model (
17) presents a less damped behaviour and gives a poorer fit of
. The similarity of the results for the original and the new method could be motivated by an averaging effect, since the estimates are based on the entire dataset.
7. Conclusions
This paper describes a method to improve the flight testing of remotely piloted aircraft for system identification of linear flight mechanical characteristics, especially when flying near the ground within a limited distance. The focus is on increasing the efficiency of flight envelope expansion tests by using real-time identification and dealing with two major issues: the short time available for each manoeuvre and the large exposure to air turbulence.
The authors make use of a relatively simple software solution to command tailored multisine input signals from a conventional radio transmitter without the need for complex on-board systems. These signals, designed according to the basic data available from previous flights, present a low correlation between the different inputs and could be applied simultaneously to different control surfaces. Results from a flight test with a subscale jet aircraft support the idea that the estimated parameters from the simultaneous multisine inputs are similar to those obtained with sequential ones, while the excitation time could be reduced by almost 50%. In addition, both multisine approaches seem to give better results than more traditional pulse inputs.
Two different methods are used for the analysis of flight test data: The first uses an existing real-time frequency domain approach, denoted here as the original method. The second incorporates instrumental variables in the original method. Results from the subscale flight test suggest that this new method improves the accuracy of the estimations for data affected by turbulent conditions.
All methods, procedures and results described in this paper contribute to reducing the gap between relatively simple and cost-effective subscale flight tests and more ambitious, resource-intensive flight test campaigns. This supports both the development of UASs, as well as the use of SFT as an economical tool for risk reduction during the early stages of full-scale aircraft design.