Mansor 2015
Mansor 2015
Abstract—This paper demonstrates a practical approach in For a particular flight regime, it would be possible to design
designing a longitudinal command stability augmentation system the aircraft to possess desirable flying and handling quality.
for unstable combat aircraft. The unaugmented aircraft is However, for an aircraft that will fly throughout an extended
originally unstable configuration in order to gain fast response flight envelope, the stability will be vary significantly, owing
for agility. The flight control system was designed in compliance primarily to change in the aircraft configuration (lowering of
with MIL-F-8785C and Gibson Criterion for the corresponding flaps and landing gear), or Mach and Reynolds number effects
flight case. The design case study is based on pitch rate command on the stability. The flying and handling quality also changes
control system using I controller. The design evaluation is based throughout the flight envelop. To possess desirable flying and
on time response analysis of simulated results using Matlab and handling quality for the entire operational envelop, the
Simulink. The pole-placement method was used to determine the command and stability augmentation system is usually
augmented feedback gain. Second order actuator dynamics
introduced to the aircraft. Strict requirements for military
model was introduced to see its effect on overall design. This
design case study offers a better understanding, an easier and
aircraft have been documented in the military specifications
practical approach implementation based on aircraft longitudinal MIL-F-8785C [1], which set out requirements for control
pitch rate command. system in order to achieve adequate de-coupling, stability and
maneuverability. Automatic controls have been used to
Keywords: stability augmentation system, I controller, pole- improve aircraft handling qualities [2, 3] .
placement, unstable aircraft, flight control Improving stability by changing the aircraft configuration
(such as horizontal tail volume ratio) will normally result in
decrease of aircraft maneuverability. Increase in stability
I. INTRODUCTION means increase in the difficulty of changing the initial
In recent years, the analysis and synthesis of flight control condition of the aircraft. This is certainly not favorable to
system design using state-space (multivariable modern control) combat aircraft which need high maneuverability for air-to-air
has been established. The command and stability augmentation combat, ground attack and other missions.
system (CSAS) is a closed loop system which provides The main task of this paper is to present a practical
artificial control to the aircraft. The aim is to reduce pilot’s approach in designing a longitudinal command and stability
work load of control and stabilize the aircraft. It is essential augmentation system for unstable combat aircraft. The aim is
feature for unstable combat aircraft due to two factors: 1) The to comply the aircraft characteristics to MIL-F-8785C and
combat aircrafts have an extended flight envelop, 2) They need Gibson Criterion for the corresponding level.
high maneuverability.
The inherent stability of an airplane depends upon the
aerodynamics stability derivatives. The magnitude of the II. MATHEMATICAL EQUATION OF AIRCRAFT MOTION
derivatives affect both the damping and natural frequency of
the aircraft on both longitudinal and lateral motion. The
The aircraft model used in this work is canard control
derivatives are depend on the aerodynamic and geometry
aircraft with geometry and aerodynamic data as reported in [4].
characteristic of the aircraft.
The mathematical model used in this work was derived
from a simplified generic linear aircraft model [5].
,(((
Let (u, w, q, ș) be the state vector where u is velocity, w is U=180 km/h U=270 km/h U=360 km/h U=450 km/h
normal velocity, q is pitch rate and ș is pitch angle. The
detailed assumptions and simplifications in force and moments 0.25 Imaginary axis
ଵ 0.1
ሶ ൌ ሺ ୶ ୶ ሻ െ െ Ʌ
୫ 0.05
ଵ
ሶ ൌ ሺ ሻ െ Ʌ (1) 0
୫ -4 -3 -2 -1 0 1 2 3
ͳ -0.05 Real axis
ሶ ൌ ሺ ሻ െ ሺ ୶ െ ሻ െ ሺଶ െ ଶ ሻ ୶୷ -0.1
୷
Ʌሶ ൌ -0.15
Phugoid
-0.2
subsystem of the aircraft are calculated and transformed to the Fig. 1. Open loop pole zero map for different speed
aircraft centre of gravity. The nonlinear model is linearized
about the trim condition at straight and level flight at various The longitudinal dynamics has short period and phugoid
speed conditions using analytical methods[6] .The linearized mode with poles at -1.96, 1.33 and 0.00151±0.185i,
model for longitudinal dynamics can be written in the state respectively, at 270 km/h. these poles locations show that the
space form as: poles of the short period and phugoid modes are unstable. The
pole locations for the other flight conditions are shown in Fig
ሶ ൌ 1. From Fig 1, shows that phugoid mode is marginally stable
ሺʹሻ but short period mode is always unstable and its instability
ൌ increases with increase in forward speed condition. This clearly
indicates the need for stability augmentation system (SAS).
Where, ൌ ሾǡ ǡ ǡ Ʌሿ and ൌ ሾɄୣ ሿ. The A and B are The purpose of a SAS is to provide satisfactory natural
system matrices as shown in (3) and are calculated at different frequency and damping ratio for the short period mode to
speed conditions. The linear models are derived at 180, 270, stabilize the aircraft.
360 and 450 km/hr. The actuator dynamic is second order
model and added to the linear system. The natural frequency
and damping ratio for the elevator actuator are 30 rad/s and 0.7 III. HANDLING QUALITIES AND CONTROL DESIGN
respectively. The actuator drives the input signal to the control
surface. Flying and handing qualities play a significant and
ሶ ୳ ୵ ୯ necessary role in control system design for manned aircraft
ሶ ୳ ୵ ୯ [7]. Cook describes: “The flying and handing qualities are
൦ ሶ ൪ ൌ ൦ ሾɄ ሿ
୳ ୵ ୯ ൪ ൦ ൪ ୣ ሺ͵ሻ properties that govern the ease and precision with which it
Ʌሶ Ͳ Ͳ ͳ Ͳ Ʌ Ͳ responds to pilot commands in the execution of the flight
task”[5]. In order to ensure the accomplishment of desired
mission safely and successfully with the minimum amount of
A. Longitudinal dynamics workload for the pilot, the control system needs to satisfy the
corresponding specification and standards. The automatic flight
The A, B, C and D matrices for straight and level flight control system design begins by determining the design goals
condition at 270km/hr of forward speed are given below: from the relevant handling qualities specification. The handling
qualities specifications widely applied in industry are military
flying qualities specifications MIL-SPEC-8785C [1], Military
ͲǤͲͳʹͷ ͵Ǥͷ Ͳ െͻǤͺͳ Standard MIL-STD-1787B [8]. It is customary to rate handling
ሾሿ ൌ ൦ െ͵Ǥ͵ͳ െ ͵ െͲǤͷͷͺ ͲǤͻͺͻ Ͳ qualities in terms of Cooper and Harper [9] levels.
൪
െͲǤͲͲͲ ʹǤͶͷ െͲǤͲͺͳ Ͳ
Ͳ Ͳ ͳ Ͳ
Ͳ ͳ Ͳ Ͳ Ͳ Ͳ A. Control Anticipation Parameter (CAP) Criterion
ሾሿ ൌ ൦ െǤͲͺ െ ͵ Ͳ ͳ Ͳ Ͳ Ͳ
൪ ሾሿ ൌ ൦ ൪ǡሾሿ ൌ ൦ ൪
ͳǤͷʹ Ͳ Ͳ ͳ Ͳ Ͳ
The CAP criterion was developed to predict the precision
Ͳ Ͳ Ͳ Ͳ ͳ Ͳ
that a pilot could expect in controlling an aircraft’s flight path.
Nevertheless, it is required to obtain a lower order equivalent
system to apply this criterion for highly augmented aircraft.
Although the normal acceleration response is determined by
the aerodynamic property, the CAP criterion aims to assess the
transient peak magnitude of angular pitching acceleration
which is mainly decided by the short period dynamics after the
pitch control input. Therefore, it is significant and universally used to evaluate acceptability of the short period mode feature
according to aerodynamic properties and the different transient response of the system and to cancel the poles due to
operating conditions[10]. The formal definition of CAP is the integral feedback. The command path pre-filter has been added
amount of instantaneous angular pitching acceleration per unit to reduced the over shoot of the system. The actuator model
of steady state normal acceleration. In (4) the value of CAP is used in this work represent by a second order system.
given in terms of second order like parameters and which is
currently in use. Since only an angle of attack and pitch rate are used for
feedback, a reduced order system is used. In this way the short
ሶ ሺͲሻ ɘଶ୬ୱ ɘଶ୬ୱ ǤǤ ଶ period dynamic are used and the phugoid dynamics are
ൌ ൌ ൌ ሺͶሻ eliminated from PRCAH system design process. Thus, the
ሺλሻ
reduced order system dynamics with inclusion of actuator and
controller dynamics are presented in (5).
CAP is evaluated graphically by parameters Ȧns and NĮ as
shown in Fig. 5, using reduced second order aircraft model for
the step response of short period (SPO) mode. Ƚሶ െͲǤͷͷͺ ͲǤͻͺͻ െͲǤͲͲͳ Ͳ Ͳ Ƚ
ۍሶ ۍ ې
ێɄሶ ʹ ێ ۑǤͶͷ െͲǤͲͺͳ ͳǤͷʹ Ͳ Ͳې ۍ ې
ۑ ێۑ
ۑ ێൌͲ ێ Ͳ Ͳ ͳ Ͳ ێ ۑɄ ۑ
B. Gibson’s Dropback Criterion ێሶ Ͳ ێ ۑ Ͳ െͻͲͲ െͶʹ Ͳێ ۑ ۑ
ۏɂሶ Ͳ ۏ ے ͳ Ͳ Ͳ Ͳ ۏ ےɂ ے
The Gibson’s Dropback criterion aims to design a Ͳ Ͳ
command and stability augmentation system (CSAS) which ې Ͳ ۍ ېͲۍ
ێ ۑ ۑ ێ
could give an aircraft with satisfactory handing qualities. This ۑ Ͳ ێୢ ۑ Ͳ ێሾୢ ሿ ሺͷሻ
ሾɄ ሿ
criterion is described as limiting values with pitch rate ۑͲͲͻێ ۑͲێ
overshoot ratio qm/qs versus the ratio of attitude Dropback to ے Ͳ ۏ ۏെͳے
steady pitch rate (¨șpeak/qs), which are shown in Fig. 2. Here,
criterion mappings are related to qualitative descriptions of the
response such as abruptness, sluggishness, and bobbling. (5) is the form of (6). The control law used in this flight
Negative Dropback is an indication of sluggishness, while control system is as given in (7). Substitute (7) in (6) and
large positive values of Dropback indicate abrupt and bobbling rearranging, results the closed loop state (8), where KT = [kĮ kq
tendencies[11]. kȘ kVȘ ki] and M= [km].
4
ሶ ൌ Ʉୢ ୢ ሺሻ
Continuous Ʉୢ ൌ െ ୢ ሺሻ
3.5
bobbling ሶ ൌ ሺ െ ሻ ሺ ሻୢ ሺͺሻ
3
2.5
The values of the gains vector, K, calculated using the pole
qmax/qss
satisfactory
2 Sluggish Abrupt
1.5
response bobble closed loop characteristic equation has a short period damping
tendency
ratio and natural frequency selected according to design
1 requirement to meet level 1 flying handling quality of MIL-
0.5 SPEC-8785C. The actuator dynamics are kept equal, and, and
the integrator zero is given a convenient value to be favorable
0
-0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 1.2 to the pilot. The time constant 1/6 second (f § 1Hz) is
DB/qss [sec] satisfactory.
Fig. 2. The boundaries of Dropback Criterion [11] The command path pre-filter is the last part to be designed
and it remains as the means for fine tuning the command
response characteristics for good handling, essentially by phase
adjustment, as shown in (9). Where Tș2actual is the inverse of
IV. CONTROL LAW DESIGN
actual pitch rate transfer function zero and Tș2desierd is the
replacement of Tș2actual. The value of Tș2desierd is obtained by
Flight control system (FCS) design is an important problem using the Gibson drop back criterion [11], which leads to the
for any aircraft. In this paper the design of a pitch rate requirement that to meet satisfactory region of the Gibson drop
command/attitude hold (PRCAH) flight control system is back criteria.
considered. The structure of the PRCAH system is shown in
Fig. (3). The inner loop feedback of angle of attack kĮ and pitch
rate kq is necessary to stabilize the aircraft. To ensure zero ଶ ୢୣୱ୧ୣ୰ୢ ͳ
ሺሻ ൌ ሺͻሻ
steady state error of PRCAH system a I controller has been ଶ ୟୡ୲୳ୟ୪ ͳ
added. In addition to the I controller the feed forward gain km
as suggested by Cook [3] has been added to improve the
km
q
qcom Pre-filter + I controller + + Actuator ɻ Aircraft
qd Dynamics dynamics ɲ
PF ki/s + -
-
[kɲ kq]
Fig. 3. Pitch rate command/attitude hold flight control system
described. 20
15
1. Select the design parameters Control anticipation
gains value, K
10
Parameter, CAP, short-term mode damping ratio, ȗSP 5
and Gibson Dropback criterion DB/qs=0. 0
2. Determine the gains [kĮ kq ki km] and Tș2desierd. 0 100 200 300 400 500 600
-5
3. Study the closed loop response in time domain
-10
characteristics.
-15
4. Use the MIL-SPEC-8785C and Gibson drop back Speed,U [km/hr]
1.5
A. longitudinal dynamics
q,[°/s]
classified as Category C, precisions non-aggressive, and A unit step pitch rate demand is applied and then removed
Category B, non-precision non-aggressive, mission tasks. For after 10 seconds. The corresponding pitch angle, ș and flight
each mission task, variability in flight conditions was assumed path angle, Ȗ responses are shown. These Gibson Dropback
to be due to changes in speed. Thus, the criteria were applied parameters are summarized in Table 1
for range of speed values (180 km/hr to 720 km/hr)
representative of flight conditions for each mission task. Using the definition presented in section III, the values of
pitch rate overshoot ratio, (qmax/qs), and the ratio of the attitude
Dropback to steady state pitch rate, (¨șpeak/qs) were calculated.
The result is presented in Fig.8.
• Control Anticipation Parameter (CAP) Criterion
Fig.8 shows the assessment of The Gibson Dropback
criteria for the closed loop transfer function with the Pitch rate
As the super augmentation system engaged, the specific command/attitude hold flight control system in place for all
damping and natural frequency of the short period mode are combinations of speed. Also shown Gibson Dropback
improved as Fig. 6 shows. It can be seen that the damping assessment for the for the closed loop transfer function with the
ratio, ȗSP and natural frequency, ȦnSP are satisfactory with the SAS and I controller case and one with the SAS, I and feed
MIL-F-8785C specifications constraints. The load factor (NĮ) forward for speed of 100m/s. It can be seen that the system
and short period frequency (ȦnSP ) comparison can be seen in falls within the boundaries of satisfactory performance when
Fig. 6(a). Using this criterion for the range of speed less than the super augmentation system is functional, while the system
540 km/hr appear in the Level 1 category while the range of with the SAS and I controller and system with the SAS, I
speed greater than 540 km/hr appears in the Level 2 category. controller and feed forward are falls within the Abrupt Bobble
Investigating Control Anticipation Parameter, CAP vs. tendency region. As shown in Fig. 8 for all combinations of
damping ratio ȗSP , as seen in Fig. 6(b). From these Fig it can speed, the augmented combat aircraft is assessed as satisfactory
be seen that the range of speed are set to be correlated so that response. This result was expected, since the requirement for
for the range of speed less than 540 km/hr has a predicted zero pitch attitudes Dropback was used in the flight control
handling quality requirement (HQR) of level 1, the speed of system design. Although not precisely zero, the actual pitch
720 km/hr has a predicted HQR of level 2. attitude Dropback is nevertheless very small. On the other hand
Increase of speed increase the overshoot ratio.
Fig.7 represents the application of Gibson pitch attitude
Dropback criteria on the close loop transfer function.
100
100
180
180 270 360 450 540 720
270
360
450
540
720
short period undamped natural frequancy wsp (rad/s)
10 10
1 1
Level 2
0.1 0.1
1 Na (1/rad) 10 100 Na (1/rad)
1 10 100
100 100
CAP
CAP
1 1 540
720
level 1
0.1 level 2
0.1
level3
0.01
0.01
0.1 1 10 0.1 1 10
short period mode damping ratio, ɺ short period mode damping ratio, ɺ
12
VI. CONCLUSION
10 The performance of a super augmentation system designed
for an unstable combat Aircraft in longitudinal axis is assessed
8 flight path angle,γ [ °] with respect to the CAP criterion and Gibson drop back
pitch angle, θ [°] requirements. A canard combat aircraft with conventional
q,[°/s], θ[ °], γ [°]
6
Pitch rate,q [ °/s] mechanical control is used for the simulation study. For the
Pitch rate demand,qd
4 simulation study, a linearized aircraft model at straight and
level flight condition is considered. The simulation results
2 shows the super augmentation system in longitudinal axis
meets the requirements of the level 1 handling qualities for the
0
range of speed up to 540 km/hr. Speed Greater than 540 km/hr
-2
may require and adaptive control to meet the level 1
0 5 10 15 20 requirement.
Time [sec]
Parameter 2 Sluggish
Pitch attitude response
Abrupt
Dropback, DB [ o ] 0.048 0.068 0.079 0.023 0.032 0.039 1.5 bobble
Pitch rate Peak, tendency
1
qmax [o/s] 1.203 1.216 1.255 1.247 1.382 1.720
Steady pitch rate, 0.5
qss [o/s] 1.000 1.000 1.000 1.000 1.000 1.000
0
Pitch over shoot
-0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 1.2
ratio, qmax/qss [-] 1.203 1.216 1.255 1.247 1.382 1.720 DB/qss [sec]
REFERENCES