Paper Aiac12 Botez PDF
Paper Aiac12 Botez PDF
Paper Aiac12 Botez PDF
net/publication/260944796
CITATIONS READS
0 320
6 authors, including:
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
All content following this page was uploaded by Ruxandra Mihaela Botez on 21 March 2014.
by
1. Introduction
Stability and control derivatives calculations play a crucial role in the simulation of an
aircraft since they determine the accuracy of the flight model as they can duplicate the
actual aircraft behavior in response to random perturbations (such as gusts and
turbulence) and/or control inputs from pilots. Concerning the calculation of these
derivatives, investigations are conducted on an aircraft expressed as a linear time
invariant state space system. The main values of stability derivatives (lift, drag and their
derivatives with respect to the angle of attack) are validated on an aircraft.
These derivatives are calculated from geometrical aircraft data. In this paper, we present
the main F/A-18 aircraft geometry data. In most of cases, these data are calculated from
references and were compared, when information existed, with data given by NASA
Dryden Flight Research Laboratories.
1. The input function, which records the needed geometry for the calculations.
2. The user built functions, which find specific graphical values based on predefined
input interpolated or not.
3. The stability derivatives, which are calculated using the methods for subsonic,
transonic and supersonic speed regimes in which the user built functions and the
inputs variables, are incorporated.
In the next section, an example for the calculation of the pitching moment coefficient
derivative with respect to the angle of attack CLα is shown. All other coefficients were
calculated by methods described [2], and we will not describe here all their calculations
methods in details.
1
2. Geometrical parameters needed to calculate the stability and control
derivatives
S 2b
Cr = = (1)
(1 + λ ) b / 2 (1 + λ ) A
where S is the reference area, b is the wing span, A is the wing area, l = Ct / Cr .
b/2
2
C= ∫ C dy
2
S 0 (2)
The root chord of exposed wing surface Cre is defined with the leading edge angles (ΛLE)
and trailing edge angles (ΛTE) for any configuration:
d
Cre = Cr − [tan(Λ LE ) − tan(ΛTE )] (3)
2
where Cre is the exposed root chord, Ct is the tip chord, b is the span, d is the fuselage
diameter.
⎛ t⎞
SWe = 2 S e ⎜ 1 + 0.2 ⎟ (5)
⎝ c⎠
2
where t is the airfoil thickness, c is the airfoil chord, Se is the exposed surface of the
wing. In this paper, the report between the thickness and the chord is constant.
(b − d )
2
Ae = (6)
Se
The leading edge and the trailing edge angles are given by the following equations:
X
X tan(Λ LE ) =
b/2
Y
tan(ΛTE ) =
Cr b/2
Ct
Y
2
tan( Δ TE ) = tan( Δ LE ) + (C t − C r )
b
2(C t − C r )(C t + C r )
tan ( Δ TE ) = tan( Δ LE ) +
b (C t + C r )
2( λ − 1)(C t + C r ) (7)
tan ( Δ TE ) = tan( Δ LE ) +
b ( λ + 1)
4( λ − 1)
tan ( Δ TE ) = tan( Δ LE ) +
A ( λ + 1)
2.8 Sweepback angle of 50% and 25% chord line
We need to determine the sweep angles at 25% and 50% from the chord. We generalize
the equations (7):
4 1− λ (8)
tan( Λ n ) = tan(Λ m ) − [( n − m) ]
A 1+ λ
3
4 1− λ
tan( Λ c / 2 ) = tan(Λ LE ) − [(0.5) ]
A 1 + λ at 50% from the chord.
4 1− λ
tan( Λ c / 4 ) = tan(Λ LE ) − [(0.25) ]
A 1 + λ at 25% from the chord.
Ct (9)
λe =
Cre
⎛ Δy⊥ ⎞ Δy
δ ⊥ = arctan ⎜ ⎟ , where Δy⊥ = (10)
⎝ 5,85 ⎠ cos( Λ LE )
• Angle of attack at which section lift curve begins to deviate from linear variation
• Section zero lift angle of attack
• Angle of attack for design lift coefficient
• Section design lift coefficient
We further explain the method to calculate the lift coefficient derivative with the angle of
attack C Lα .
4
3. Lift coefficient derivative with the angle of attack C Lα
For a conventional, tail-aft airplane with subsonic Mach number, this derivative is
combination of wing W, fuselage B, and tail T contributions. The aircraft derivative can
be estimated [1]:
SH ⎛ dε ⎞
CLα = CLαWB + CLα H η H ⎜1 − ⎟ (11)
S ⎝ dα ⎠
where ε is the downwash angle at the horizontal tail, ηH is the ratio of dynamic pressure
S
at the tail to the free stream dynamic pressure and the area ratio is denoted by H . Since
S
the ratio of wing span b to fuselage diameter d is reasonably large, following
approximation may be done:
2
⎛d ⎞ ⎛d⎞
where KWB = 1 − 0.25 ⎜ ⎟ + 0.25 ⎜ ⎟ (12)
⎝b⎠ ⎝b⎠
Compressibility correction factor is written as β = 1 − M 2 . The average wing section lift
curve slope is:
π a∞ A
Clα = cos Λ LE (13)
a∞ + π A2 + 4
2π AW
CLαW = (14)
AW β 2
tan 2 Λ c / 2W
2
2+ (1 + )+4
k2 β2
2π AH
CLα H = (15)
AH 2 β 2 tan 2 Λ c / 2 H
2+ (1 + )+4
k2 β2
CLαW
where k= , the downwash ratio at the horizontal tail is
2π
dε dε CLαW
= M
where ε is the downwash angle, i.e. the mean angle at the tail
dα M dα M =0 CLαW
M =0
5
dε
plane that the airflow is deflected through by wing. The downwash derivative is a
dα
function of the wing geometry and position of the tailplane relative to the wing. The
downwash gradient at low speed,
dε 1.19
= 4.44 ⎡⎣ K A K λ K H cos Λ c / 4 ⎤⎦ where:
dα M =0
hH
1−
1 1 10 − 3λ b
KA = − , Kλ = and K H = (16)
A 1+ A 1.7
7 2lH
3
b
CLαW at Mach number equal to zero M = 0 is given by the following equation:
2π AW
CLαW =
M =0
AW 2
2+ (1 + tan 2 Λ c / 2W ) + 4
k2
Other equations are derived for the same coefficient from [2], by use of interpolation
tables and figures. Results expressed in terms of stability derivatives for four flight cases
are represented in the next section.
4. Results
Flight condition 1
6
Flight condition 2
Flight condition 3
Flight condition 4
7
5. Conclusions
Main stability derivatives (lift, drag and their derivatives with respect to the angle of
attack) were calculated for an F/A-18 aircraft for a number of four flight conditions.
These values were compared with the values given by the NASA DFRC researchers.
This work is applied and validated on the F/A-18 aircraft and is useful for future
aircraft simulation and control in aerospace industry. One advantage of this work is that
the execution time needed for the lift and drag coefficients calculations is very small
compared to the time needed to calculate these coefficients based on Computational Fluid
Dynamics (CFD) computer knowledge.
Acknowledgments
The authors would like to thank Mr. Marty Brenner at NASA Dryden Research Flight
Centre for his collaboration.
References
[1] Roskam, J., “Methods for Estimating Stability and Control Derivatives of
Conventional Subsonic Airplanes”, University of Kansas, Lawrence, Kansas.
[2] Hoak, D.E., Finck, R.D., “USAF Stability and Control Derivatives DATCOM”,
1978.