Abstract

Despite significant technical advances, pilot-induced oscillation (PIO) continues to occur in both flight tests and operational aircrafts. Such a phenomenon has led to significant research activities that aim to alleviate this problem. In this paper, the L1 adaptive controller has been introduced to suppress the PIO, which is caused by rate limiting and pure time delay. Due to its architecture, the L1 adaptive controller will achieve a desired response with fast adaptation. The analysis of PIO and its suppression by L1 adaptive controller are presented in detail in the paper. The simulation results indicate that the L1 adaptive control is efficient in solving this kind of problem.

1. Introduction

The high performance demands of modern aircrafts, especially highly maneuverable military jets, require the implementation of advanced control systems. The use of a modern electronic flight control system could provide great potential for improvement in the aircraft’s performance. In spite of all the improvements, a significant handling quality problem arose with the introduction of electronic flight control system: pilot-induced oscillation, also known as pilot in-the-loop oscillation.

Pilot-induced oscillations are described as pilot-aircraft dynamic couplings, which could lead to instability in the systems [1]. Both previous and current research has attempted to explain, predict, and avoid these oscillations. Almost every modern aircraft has experienced PIO, which is well known by the public for the catastrophic event it caused, such as the YF-22 [2] and Olympic Airways Falcon 900. The occurrence of such events has led to significant research activities that are intended to alleviate the negative effects due to PIO. Despite the research efforts made, PIOs continue to occur, and reports of PIOs on operational aircrafts are increasing.

The focus of this paper is to suppress the pilot-induced oscillation caused by both rate limiting and pure time delay by using the 𝐿1 adaptive controller. The effects of rate limiting and other system nonlinearities are considered to be the main factors that result in the occurrence of PIOs [3–5]. Some of the existing methods, to some extent, could handle the PIO caused by rate limiting, but not necessarily work with the presence of pure time delay in the actuator model. The 𝐿1 adaptive control is known for its fast adaptation and smooth control implementation due to its powerful control architecture [6, 7]. With proper design, the 𝐿1 adaptive controller will make the system respond in the desired manner.

The objective of this paper is to design the 𝐿1 adaptive controller to make the inner loop of the system respond according to a given first-order system, while suppressing the PIO phenomenon caused by both rate limiting and pure time delay in the pilot dynamic model. Section 2 gives a brief introduction of pilot-induced oscillation and pilot model. In Section 3, the design of an inner loop adaptive controller to suppress the PIO and to track desired response is given. The inner loop adaptive controller cannot handle the adversely high pilot command inputs because they are the commands to the inner closed loop system. The inner loop controlled by 𝐿1 adaptive controller solely tracks the command inputs from the pilots.

2. Pilot-Induced Oscillation

Pilot-induced oscillations can be regarded as a closed-loop instability of the pilot-aircraft loop, and it often happens when the pilot proves to be unable to adapt himself to a sudden change of the vehicle dynamics during a high demanding flight task. PIOs are complex interactions between the pilot and the aircraft dynamics.

2.1. Types of PIOs

According to the report by National Research Council (NRC) Committee in 1997 [8], PIOs can be separated into three categories.

(i) Category I PIO: Linear pilot-vehicle system oscillation. These PIOs result from linear phenomena such as excessive time delay and excessive phase loss between the pilot’s control input and the aircraft response. They are the simplest to model, understand, and prevent. Before the introduction of fly-by-wire technology, almost all PIOs were this type, and a large amount of research has been conducted on them. But they are currently the least common in operational flight. The main causes of Category I PIO are excessive lags due to time delays and various digital filter dynamics in the flight control system. These effects lead to a high frequency phase roll-off in the frequency response. There are roughly two groups of design criteria for preventing Category I PIO. One is according to the flight control stability, and the other is the handling quality requirements. In the first group of design criteria, the open loop aircraft frequency response needs to be checked as well as the phase and magnitude margin. If the phase rate of the system is too high, a small increase in frequency will result in a strong additional phase delay. This is usually related with time delay in the whole system.

(ii) Category II PIO: Quasi-linear events with some nonlinear contribution, such as rate or position limiting. For the most part, these PIOs can be modeled as linear events, with an identifiable nonlinear contribution that may be treated separately. Figure 1 shows two typical positions where rate-limiting blocks are installed in the flight control system. The block after the pilot model is installed to prevent the system from receiving a high input rate by the pilot. The other block after the controller is used to protect the actuator against overload. In the case of rate limiting, nonlinear system theory should be utilized for the analysis of the flight system, such as the Lyapunov theory and the phase plane method.

The majority of aircraft crashes due to PIOs are caused by Category II PIO through activating actuator rate limits. When the rate limiter is saturated, phase lag occurs and the aircraft dynamics change suddenly, which cause PIO and threaten the flight safety. Rate limiting adds additional phase lag, increasing the delay between the pilot input and the aircraft response. This tends to make the pilot compensate with faster responses, often worsening the situation. Rate limiting also reduces the gain, which the pilot interprets as a lack of control response and therefore makes larger command inputs, again making the situation worse.

(iii) Category III PIO: Nonlinear PIOs with transients. Such events are difficult to recognize and rarely occur but are always severe. Mode switching or rapid changes in effective vehicle characteristics could be reasons for this type of PIO. Fortunately, these events that result in nonoscillatory divergence and loss of control of the aircraft are rare.

2.2. Rate Limiting

Rate limiting of the actuator is one of main factors that lead to pilot-induced oscillation. Actuator rate limiting occurs when pilot input command error requires a higher rate than the actuator can actually provide. A simplified model of a rate-limited actuator is shown in Figure 2.

Rate limiting adds additional phase lag, increasing the delay between the pilot input and aircraft response. This tends to make the pilot compensate with faster responses, which will often worsen the situation.

2.3. Pilot Modeling

With rate limiting considered, another important factor that causes the PIO phenomenon is the pilot model, which is also regarded as the weakest point in the analysis due to its high nonlinearity and complexity. However, there are several specific pilot behavioral patterns [9] that could be used to analyze the cause of these PIOs. The pilot model is the source factor that distinguishes severe PIO problems from most aircraft feedback control design problems. The difference resides in unique human properties related to the adaptive characteristics of the human pilot.

Pilots exhibit peculiar transitions in the organizational structure of the pilot vehicle system. These transitions can involve both the pilot’s compensation and effective architecture of the pilot’s control strategy. The human pilot dynamics can be roughly separated into the following types.

(i) Compensatory behavior: essentially, the pilot has generated a lead to cancel out the lag in the aircraft model. However, the higher frequency lags of the pilot model can be approximated at the lower frequencies by the pure time delay.

(ii) Pursuit behavior: the introduction of the pursuit behavior permits an open-loop control in conjunction with the compensatory closed-loop error correcting action. The pilot model of both compensatory behavior and pursuit behavior will be superior to that where only compensatory operations are possible.

(iii) Precognitive behavior: this kind of pilot model gives a higher level of control performance. Based on the knowledge of the system dynamics, the pilot model will generate proper control signals at the right time so as to result in machine outputs that are almost as desired. This operation also appears in company with compensatory behavior as well as pursuit behavior. Most highly skilled movements will automatically fall into this category.

3. PIO Suppression by 𝐿1 Adaptive Control

3.1. System Modeling

In most PIO cases that have happened in recent years, rate- and position-limiting phenomena appear. Hence, numerous research works have been conducted to analyze rate and position limiting [10–12]. In this situation, PIOs can be explained by limit cycles occurring in a nonlinear system where the nonlinearities cause a sustained, constant amplitude oscillation. In those studies, rate limiting of the actuator is considered as one possible cause for PIOs. There are some other factors, such as relatively high pilot command gain and lag in the pilot model, contributing to the onset of PIOs. Rate limiting is commonly met in aviation practice. Hence, if an adaptive inner loop can compensate for the nonlinearities caused by rate limiting and eliminate the limit-cycle oscillations, it can be considered that the PIO is suppressed for such commonly seen situations. We apply the 𝐿1 adaptive controller on the F-14 model taken from [10] in order to test this idea. A model of the Grumman Aircraft Company (GAC) F-14 was developed for this investigation. The approximated model in [10] is compared with data of a real F-14 that has severe PIOs during flight tests. The rate-limiting phenomenon we considered is due to sudden off-nominal conditions such as hydraulic pressure failure as in the F-14 example. Therefore, it is hard to estimate the nonlinearity caused by rate limiting in advance, like most other research work has done. Adaptive control is needed for such situations.

The control structure is shown in Figure 3. It is a longitudinal dynamic model. The system states are [π›Όπœƒπ‘ž], and control input is the elevator. The 𝐿1 adaptive inner loop controller is a pitch-rate augmentation system, which provides the desired dynamics for the [π›Όπ‘ž] subsystem. The outer loop takes the πœƒ feedback signal and injects it into the pilot model. The pilot model outputs a pitch rate command for the inner loop to achieve the pitch attitude control (the πœƒ angle). The pilot model is given by𝑝(𝑠)=𝐾(𝑠+𝛽)π‘’βˆ’πœπ‘ ,(1) which describes the compensatory behavior of the pilot.

3.2. Delay Margin Detection

Consider the Pilot model described by (1); we assume the desired response of the inner loop is given by𝑔𝑐1(𝑠)=𝑇𝑝.𝑠+1(2) With the Pilot model, we could derive the transfer function of the closed-loop as𝑔(𝑠)=𝑝(𝑠)𝑔𝑐(𝑠)𝑝(𝑠)𝑔𝑐=(𝑠)+1𝐾(𝑠+𝛽)π‘’βˆ’πœπ‘ π‘‡π‘π‘ +1+𝐾(𝑠+𝛽)π‘’βˆ’πœπ‘ .(3)

In the PIO problem concerned, the pure time delay deviates between 0.2 and 0.3 second. The time delay term can be approximated by π‘’βˆ’πœπ‘ =(1βˆ’(πœπ‘ /2))/(1+(πœπ‘ /2)) according to Pade approximation. Equation (3) can be further written as=𝑔(𝑠)𝐾(𝑠+𝛽)1+πœπ‘ /2𝑇𝑝𝑒𝑠+1(1+πœπ‘ /2)+𝐾(𝑠+𝛽)(1βˆ’πœπ‘ /2)βˆ’πœπ‘ =(𝐾𝜏/2)𝑠2+(1+πœπ›½πΎ/2)𝑠+π›½πΎπ‘‡ξ€·ξ€·π‘ξ€Έξ€Έπ‘ βˆ’πΎπœ/22+𝑇𝑝𝑒+πΎβˆ’πœ/2βˆ’πœπ›½πΎ/2𝑠+1+πΎπ›½βˆ’πœπ‘ ,(4) according to Routh stability criterion, we can get the following inequalities:𝑇𝑝𝑇>𝐾,1+𝐾𝛽>0,𝑝>πΎπœπ›½/2βˆ’πœ/2βˆ’πΎ.(5)

If the parameter 𝑇𝑝 chosen is satisfied with unequal conditions in (5), the system will achieve the desired response with stability condition. From many PIO studies, the onset of PIO often was accompanied by high pilot command inputs. If some off-nominal conditions happen, such as hydraulic failure of the actuators in F-14 PIO accidents, the rate limiting comes into play. When the 𝐿1 inner loop adaptive controller is applied to this model, it can compensate for the unknown effects caused by rate limiting and time delay, no matter whether they are linear or nonlinear. The design of the 𝐿1 adaptive controller does not need a priori information for the rate limits of the actuator.

3.3. 𝐿1 Adaptive Control Design to Achieve Desired Response

Consider the 𝐿1 architecture given in Figure 4, the 𝐿1 adaptive control [13] could be divided into three parts: the adaptive law, the state predictor, and the control law.

State Predictor Design
The state predictor of the aircraft longitudinal dynamic system in Figure 3 is designed as follows: Μ‡Μ‚π‘₯(𝑑)=π΄π‘šΜ‚π‘₯(𝑑)+𝐡𝑒(𝑑)+𝜎π‘₯(𝑑),Μ‚π‘₯(0)=π‘₯0,(6) where π΄π‘šβˆˆπ‘…3Γ—3,π΅βˆˆπ‘…3Γ—1, Μ‡βŽ‘βŽ’βŽ’βŽ’βŽ’βŽ£Μ‚βŽ€βŽ₯βŽ₯βŽ₯βŽ₯βŽ¦Μ‚π‘₯(𝑑)=πœƒ(𝑑)Μ‚π‘ž(𝑑)𝛼(𝑑),𝜎π‘₯⎑⎒⎒⎒⎒⎣(𝑑)=𝜎1(𝑑)𝜎2(𝑑)𝜎3⎀βŽ₯βŽ₯βŽ₯βŽ₯⎦.(𝑑)(7) Matrix A should be Hurwitz to make sure of the stability of the model. 𝜎 can be divided into two parts, the model matched part and the model unmatched part 𝜎=π΅ξπœŽπ‘š+π΅ξπœŽπ‘’π‘š,(8) where 𝐡 is the null space of B. ξπœŽπ‘š is the matched part, and ξπœŽπ‘’π‘š indicates the unmatched part. ξπœŽπ‘š and ξπœŽπ‘’π‘š can be calculated as follows: βŽ‘βŽ’βŽ’βŽ£ξπœŽπ‘šξπœŽπ‘’π‘šβŽ€βŽ₯βŽ₯⎦=ξ‚ƒπ΅π΅ξ‚„βˆ’1𝜎.(9)

Adaptive Law Design
Assume the system dynamics is as follows: Μ‡π‘₯(𝑑)=𝐴π‘₯(𝑑)+𝐡𝑒+𝜎(𝑑).(10) Given any 𝑇>0, we have Ξ¦π‘₯ξ€œ(𝑇)=𝑇0π‘’π΄π‘š(π‘‡βˆ’πœ)π‘‘πœ.(11) Letting Μƒπ‘₯(𝑑)=Μ‚π‘₯(𝑑)βˆ’π‘₯(𝑑), the adaptive law for 𝜎π‘₯(𝑑) is given by 𝜎π‘₯(𝑑)=𝜎π‘₯[],(𝑖𝑇),π‘‘βˆˆπ‘–π‘‡,(𝑖+1)π‘‡ξπœŽπ‘₯(𝑖𝑇)=βˆ’Ξ¦π‘₯βˆ’1(𝑇)πœ‡π‘₯(𝑖𝑇),𝑖=0,1,2,…,(12) where Ξ¦(𝑇) is defined above, and πœ‡π‘₯(𝑖𝑇)=π‘’π΄π‘šπ‘‡Μƒπ‘₯(𝑖𝑇),𝑖=1,2,β€¦ξπœŽπ‘₯(𝑑)=βˆ’Ξ¦π‘₯βˆ’1(𝑇)π‘’π΄π‘šπ‘‡Μƒπ‘₯(𝑑).(13)

Control Law Design
The output estimation of the system is given by ̂𝑦=𝐢̂π‘₯.(14) The control signal can be calculated by 𝑒=π‘’π‘Ÿ+π‘’π‘š+π‘’π‘’π‘š,π‘’π‘Ÿπ‘Ÿ=βˆ’πΆπ΄π‘šβˆ’1𝐡,π‘’π‘š=βˆ’ξπœŽπ‘š,π‘’π‘’π‘š=ξ‚ƒπΆξ€·π‘†πΌβˆ’π΄π‘šξ€Έβˆ’1π΅ξ‚„βˆ’1πΆξ€·π‘†πΌβˆ’π΄π‘šξ€Έβˆ’1π΅ξπœŽπ‘’π‘š.(15)π‘’π‘Ÿ is the control signal used to track the reference signal with desired response designed aforementioned. π‘’π‘š and π‘’π‘’π‘š are control signals used to cancel the uncertainty part of the system, where π‘’π‘š is the matched control signal and π‘’π‘’π‘š is the unmatched one.

4. Simulation Study

In general, rate limiting results in an amplitude reduction and a significant added phase lag. With the introduction of pure time delay into the pilot model, the closed-loop system could even become unstable.

Based on the delay margin detection architecture shown in Figure 5, the time variable in the desired system response model can be obtained from the result derived in the previous section. There is another alternative method that could be used for the detection of the delay margin by searching for the optimal value of 𝑇𝑝 in (2). With an increase in 𝑇𝑝, the system response will change gradually in accordance. At some critical point, the oscillation will disappear and the system response will converge to the set-point value instead. If we further increase the time variable, the damping ratio of the closed-loop system will increase as well, which will affect the system’s transient response including overshoot, rise time and so forth (Figure 6).

In Figure 7, the F-14 system responses of different time constants 𝑇𝑝 are given. If the time variable 𝑇𝑝 is small as shown in Figure 7(a), the system responds rapidly, but the overshot is also very large. With increasing 𝑇𝑝, the overshoot decreases, but transient performance reduces as well. Thus, we need to build a balance between the overshoot and the rise time, which can be simply implemented by introducing an objective function asmin𝐽=πœ”1𝛿+πœ”2π‘‡π‘Ÿ,(16) where πœ”1,πœ”2 are weighting coefficients and 𝛿,π‘‡π‘Ÿ represent the overshoot and rise time, respectively.

5. Conclusion and Future Research

This paper proposes an alternative method for suppressing the pilot-induced oscillation phenomenon caused by rate limiting and time delay in the pilot dynamic model. The delay margin is detected at the beginning, and the critical part of the method is design the 𝐿1 adaptive controller to achieve the desired response of the system.

To further improve the performance of this method, we need to consider the following issues. First, certain criterion need to be built to determine the choice of 𝑇𝑝. The system response can be evaluated by stability, overshoot, and rise time, and so forth. For the robustness of this method, a time-varying pilot model should be considered.

Acknowledgment

The authors would like to thank Dr. Jiang Wang from Zona Inc. for his useful suggestions and help.