Skip to main content

Stability and Control Estimation Flight Test Results for the SR-71 Aircraft With Externally Mounted Experiments

NASA/TP-2002-210718 · NASA (NTRS) · 2002

Public domain · NASA (NTRS)Technical Reports

Overview

A maximum-likelihood output-error parameter estimation technique is used to obtain stability and control derivatives for the NASA Dryden Flight Research Center SR-71A airplane and for configurations that include experiments externally mounted to the top of the fuselage. This research is being done…

Publisher
NASA (NTRS)
Document
NASA/TP-2002-210718
Year
2002
Pages
96

Document

NASA/TP-2002-210718

Stability and Control Estimation Flight Test

Results for the SR-71 Aircraft With

Externally Mounted Experiments

Timothy R. Moes and Kenneth lliff NASA D_den Flight Research Center Edwards, California

June 2002

The NASA STI Program Office...in Profile

Since its founding, NASA has been dedicated CONFERENCE PUBLICATION.

to the advancement of aeronautics and space

Collected papers from scientific and

science. The NASA Scientific and Technical

technical conferences, symposia, seminars,

Information (STI) Program Office plays a key

or other meetings sponsored or cosponsored

part in helping NASA maintain this

by NASA.

important role.

SPECIAL PUBLICATION. Scientific,

technical, or historical information from

The NASA STI Program Office is operated by

NASA programs, projects, and mission,

Langley Research Center, the lead center for

often concerned with subjects having

NASA's scientific and technical information.

substantial public interest.

The NASA STI Program Office provides access

to the NASA STI Database, the largest collection

TECHNICAL TRANSLATION. English-

of aeronautical and space science STI in the

language translations of foreign scientific

world. The Program Office is also NASA's

and technical material pertinent to

institutional mechanism for disseminating the

NASA's mission.

results of its research and development activities.

These results are published by NASA in the

Specialized services that complement the STI

NASA STI Report Series, which includes the

Program Office's diverse offerings include

following report types:

creating custom thesauri, building customized

databases, organizing and publishing research

TECHNICAL PUBLICATION. Reports of

results...even providing videos.

completed research or a major significant

phase of research that present the results of

For more information about the NASA STI

NASA programs and include extensive data

Program Office, see the following:

or theoretical analysis. Includes compilations

of significant scientific and technical data

• Access the NASA STI Program Home Page

and information deemed to be of continuing

at http ://www.sti.nasa.gov

reference value. NASA's counterpart of

peer-reviewed formal professional papers but

• E-mail your question via the Internet to

has less stringent limitations on manuscript

help@sti.nasa.gov

length and extent of graphic presentations.

• Fax your question to the NASA Access Help

TECHNICAL MEMORANDUM. Scientific

Desk at (301) 621-0134

and technical findings that are preliminary or

of specialized interest, e.g., quick release

• Telephone the NASA Access Help Desk at

reports, working papers, and bibliographies

(301) 621-0390

that contain minimal annotation. Does not

contain extensive analysis.

Write to:

NASA Access Help Desk

CONTRACTOR REPORT. Scientific and

NASA Center for AeroSpace Information

technical findings by NASA-sponsored 7121 Standard Drive

contractors and grantees. Hanover, MD 21076-1320

NASA/TP-2002-210718

Stability and Control Estimation Flight Test

Results for the SR-71 Aircraft With

Externally Mounted Experiments

Timothy R. Moes and Kenneth lliff

NASA Dryden Flight Research Center

Edwards, California

National Aeronautics and

Space Administration

Dryden Flight Research Center

Edwards, California 93523-0273

June 2002

NOTICE

Use of trade names or names of manufacturers in this document does not constitute an official endorsement of such products or manufacturers, either expressed or implied, by the National Aeronautics and Space Administration.

Available from the following:

NASA Center for AeroSpace Information (CASI) National Technical Information Service (NTIS)

7121 Standard Drive

5285 Port Royal Road

Hanover, MD 21076-1320

Springfield, VA 22161-2171

(301) 621-0390 (703) 487-4650

CONTENTS

iii

TABLES

2. Zero-fuel weight, CG, and mass moment of inertia information for

FIGURES

5. Time history of an acceleration from subsonic to supersonic flight

6. Time history of an acceleration from subsonic to supersonic flight

7. Body-axis mass moments of inertia for the baseline configuration using the

16. Predicted and flight-determined angle-of-attack derivatives

iv

29. Flight conditions for comparison of configuration longitudinal stability

36. Predicted and flight-determined angle-of-sideslip derivatives

55. Flight conditionsfor comparisonof configurationlateral-directionalstability

vi

ABSTRACT

A maximum-likelihood output-error parameter estimation technique is used to obtain stability and control derivatives for the NASA Dryden Flight Research Center SR-71A airplane and for configurations that include experiments externally mounted to the top of the fuselage. This research is being done as part of the envelope clearance for the new experiment configurations. Flight data are obtained at speeds ranging from Mach 0.4 to Mach 3.0, with an extensive amount of test points at approximately Mach 1.0.

Pilot-input pitch and yaw-roll doublets are used to obtain the data. This report defines the parameter estimation technique used, presents stability and control derivative results, and compares the derivatives for the three configurations tested. The experimental configurations studied generally show acceptable stability, control, trim, and handling qualities throughout the Mach regimes tested. The reduction of directional stability for the experimental configurations is the most significant aerodynamic effect measured and identified as a design constraint for future experimental configurations. This report also shows the significant effects of aircraft flexibility on the stability and control derivatives.

NOMENCLATURE

normal acceleration (positive up), ft/sec 2

a n

lateral acceleration (positive toward the right), ft/sec 2

ay

b

reference span, 56.7 f'_

BL

butt line, in.

mean aerodynamic chord, 37.7 ft

CG

center of gravity, percent c

coefficient of rolling moment

Cl

rolling moment bias, ]inear coefficient estimate for 13 = 0 °

CI b

derivative of rolling moment due to nondimensional roll rate, 3CI/3(pb/2VR), rad -1

Clp

derivative of rolling moment due to nondimensional yaw rate, _Ci/3(rb/2VR), rad -1

CI_

derivative of rolling moment due to sideslip, OCI/_[5, deg -1

Cl_

derivative of rolling moment due to aileron, OCl/O_ a , deg -l

CI 8

tl

derivative of rolling moment due to rudder, OCl/O_) r , deg -l

Cl 8

r

coefficient of pitching moment

Cm

pitching moment bias, linear coefficient estimate for ot = 0 °

Cm b

C derivative of pitching moment due to nondimensional pitch rate, OCm/O(qc/2 VR), rad -1

mq

derivative of pitching moment due to angle of attack, OCm/OO_, deg -l

Cm(t

C

derivative of pitching moment due to elevon, OCm/t)_) e , deg -1

m 6 e

Ct l coefficient of yawing moment

yawing moment bias, linear coefficient estimate for 13 = 0 °

Cn b

C n derivative of yawing moment due to nondimensional roll rate, bCn/_(pb/2 VR), rad -l

P

C

derivative of yawing moment due to nondimensional yaw rate, bCn/b(rb/2VR ), rad -1

?1r

C

derivative of yawing moment due to sideslip, _Cn/_13, deg -l

nf_

derivative of yawing moment due to aileron, bCn/_8 a , deg -l

Cna

a

derivative of yawing moment due to rudder, _Cn/c_r, deg -1

Cna

r

coefficient of normal force

CN

normal force bias, linear coefficient estimate for _ = 0 °

CN b

derivative of normal force due to nondimensional pitch rate, OCN/_)(qc/2 VR), rad -1

CNq

derivative of normal force due to angle of attack, OCN/OO_, deg -1

C N_t

derivative of normal force due to elevon, OCy/O_e, deg -1

CN a

e

coefficient of side force

Cr

side force bias, linear coefficient estimate for 13 = 0 °

CY b

derivative of side force due to nondimensional roll rate, OCy/O(pb/2VR), rad -1

Cyp

C

derivative of side force due to nondimensional yaw rate, OCy/O(rb/2 VR), rad -1

Yr

derivative of side force due to sideslip, 0Cy/_)13, deg -1

Cy_

derivative of side force due to aileron, OCy/_8 a , deg -1

Cya

a

C

derivative of side force due to rudder, _)Cy/_)6 r , deg -t

Ya F

F state derivative function

FS

fuselage station, in.

acceleration of gravity, ft/sec 2

g

G

response function

rolling moment of inertia, slug-ft 2

cross product of inertia, slug-ft 2

pitching moment of inertia, slug-ft 2

/y

yawing moment of inertia, slug-ft 2

cost function

J(¢)

KEAS

equivalent airspeed, knots

LASRE

Linear Aerospike SR-71 Experiment

aircraft mass, slug

m

number of time history points used

t/t

number of response variables

//z

roll rate, deg/sec

P

roll acceleration, deg/sec 2

P

PID

parameter identification

pitch rate, deg/sec

q

pitch acceleration, deg/sec 2

dynamic pressure, lbffft'

?-

yaw rate, deg/sec

yaw acceleration, deg/sec 2

R conversion factor, 57.2958 deg/rad

SR-71 reference area, 1605 ft 2

S

SAS

stability augmentation system

t time, sec

discrete time point at ith data point

ti

U measured control input vector

V

true airspeed, ft/sec

W

weight, lb

W

response weighting matrix (used in the cost function)

WL water line, in.

state vector

X

time derivative of the state vector

y_

normal accelerometer location, ft aft of the CG

Xa n

lateral accelerometer location, ft aft of the CG

Xay

angle-of-attack measurement location, ft aft of the CG

X(I

angle-of-sideslip measurement location, ft aft of the CG

xf_

normal accelerometer location, ft to the right of the CG

YO n

lateral accelerometer location, ft to the right of the CG

Yay Z

response vector (measurement vector)

normal accelerometer location, ft above the CG

Z-a n

lateral accelerometer location, ft above the CG

Zay

angle-of-sideslip measurement location, ft above the CG

z_

wing reference plane angle of attack, deg

time rate of change of angle of attack, deg/sec

angle of sideslip, deg

time rate of change of angle of sideslip, deg/sec

control-surface deflection, deg

aileron deflection, deg

elevon deflection, deg

8r rudder deflection, deg

lateral stick position, in.

rla

longitudinal stick position, in.

tie

pitch angle, deg

time rate of change of pitch angle, deg/sec

stability and control derivative parameter vector

roll angle, deg

time rate of change of roll angle, deg/sec

transpose

estimated response parameter

INTRODUCTION

A Mach 3.2--capable SR-71 airplane has completed a series of flight tests at the NASA Dryden Flight Research Center (Edwards, California). The series was performed to determine stability and control characteristics of the baseline corffiguration (fig. 1) and two configurations with experiments mounted on top of the fuselage. NASA Dryden previously modified the internal structure of one of its SR-71 aircraft to accommodate experiments weighing a maximum of 14,500 lb for high-speed flight research of new and unique concepts. The Linear Aerospike SR-71 Experiment (LASRE) is one example of such flight research (ref. 1) and consisted of an approximately 14,140-1b payload weight that was mounted to the SR-71 upper fuselage (fig. 2). The LASRE configuration obtained flight data at speeds to a maximum of Mach 1.75.

After the termination of the LASRE program, a four-flight test bed configuration flight program (ref. 2) was conducted to a maximum speed of Mach 3.0. The test bed configuration consists of the LASRE pod without the half-span lifting-body model (fig. 3). Because of the large size of the LASRE and test bed configurations, an incremental stability and control flight envelope expansion became necessary to ensure safe flying characteristics. The envelope expansion includes pilot-input doublet maneuvers for parameter identification (PID) of stability and control derivatives at increasing Mach

numbers. Both pitch doublets and yaw-roll doublets have been performed at each Mach number

condition. Data also have been obtained at low and high equivalent airspeeds to determine the effect of aircraft flexibility on the stability and control derivatives. A maximum-likelihood output-error program (ref. 3) has been used postfligh! to estimate stability and control derivatives from the flight data. This report presents flight-determined stability and control results for the SR-71 baseline, LASRE, and test bed configurations. Stability ar, d control derivative predictions from the SR-71 aerodynamic model (ref. 4) also are presented for the baseline configuration.

VEHICLE DESCRIPTION

An SR-71A airplane (Lockheed Martin Corporation, Palmdale, California) is used as the carrier vehicle for the LASRE and test bed configurations. Stability and control data have been obtained for the baseline, LASRE, and test bed configurations that are described in this section. Because the SR-71 aircraft is fairly flexible, aeroelastic effects on the stability and control derivatives also have been obtained. Control surfaces are not modified for the LASRE and test bed configurations.

Baseline Configuration

The SR-71A aircraft is a two-place, twin-engine aircraft capable of cruising at speeds to a maximum of Mach 3.2 and altitudes to a maximum of 85,000 ft. The aircraft is powered by two 34,000-1bf thrust-class J58 (Pratt & Whitney, West Palm Beach, Florida) afterburning turbojet engines.

Approximately 5 percent of thrust enhancement is obtained on both engines by increasing the turbine exhaust gas temperature and rotor speed (ref. 2). The engine is aligned with the wing reference plane, which has a 1.2-deg nosedown incidence compared to the fuselage centerline reference plane. The angle of attackusedhereinis referencedto the wing reference plane.The engineinlet is cantedslightly down andinwardto obtainlow local flow angles.

Control-surfaceactuatorsarepoweredusing two independenthydraulic systems.Twin all-moving, tetrahedrallyshaped,vertical fins mounted on top of the enginenacellesprovide directional control; inboardandoutboardelevonsprovidelongitudinal andlateralcontrol. The inboardandoutboardsurfaces simultaneously move; however,the outboardelevonsareriggedwith 3° moretrailing-edge-upincidence thanthe inboardelevons.The maximumcontrol-surface position limits, however,are the samefor both inboardandoutboardsurfaces. Table 1 lists the SR-71maximumcontrol-surfaceposition andratelimits.

Table 1. SR-71 aerodynamic control-surface position and rate limits.

Surface Position limit, Rate limit,

deg deg/sec

Inboard elevons:

Trailing edge down 20 30

Trailing edge up 35 30

Outboard elevons:

Trailing edge down 35 30

Trailing edge up 35 30

Rudders:

Trailing edge left 20 33

Trailing edge right 20 33

All control-surface positions have been measured except that of the right outboard elevon. Because the right outboard elevon position is not instrumented, only the inboard surface positions have been used to define the elevon and aileron deflections (6e and 8 a , respectively) for the PID analysis. The additional 3 ° of trailing-edge-up outboard elevon position has been verified using the left wing inboard and outboard position measurements. The control-surface deflections are defined as follows:

8 a = (Left inboard t5 - Right inboard t5)/2

_e = (Left inboard 8 + Right inboard 8)/2

_r = (Right 8 + Left _5)/2

The flight control system includes three-axis auto pilot and stability augmentation systems (SASes).

To minimize trim drag, the baseline longitudinal open-loop static stability is designed to be slightly positive at Mach 3.2 with the center of gravity (CG) at 0.25 mean aerodynamic chord (c) (which is the operationalaft CG limit). A redundant pitch SAS is used to provide good closed-loop handling qualities at all Mach numbers. The pitch SAS uses high-passed pitch rate to augment damping and lagged pitch rate to slightly augment stability (ref. 5). The SR-71 baseline configuration is also designed to have a minimum, but still positive, derivative of yawing moment due to sideslip _ " _e(C'q3/ at Mach 3.2. A yaw SAS is required to provide acceptable handling qualities at high Mach numbers and to prevent extreme sideslip transients caused by potential inlet "unstarts." The yaw SAS uses yaw-rate feedback for damping and uses lateral acceleration feedback to augment stability. The effective closed-loop directional stability provided by the yaw SAS can be computed using the following equation:

Closed-loop C% = Open-loop C% + 30C% _-Cyf_ (1)

r A roll SAS is used to provide roll damping through roll-rate feedback.

The empty weight of the SR-71 baseline configuration is approximately 60,700 lb. The SR-71 aircraft has a maximum fuel capacity of 80,000 lb. For these baseline configuration tests, fuel loads of a maximum of 62,000 lb were used.

Linear Aerospike SR-71 Experiment Configuration

The LASRE configuration was developed to obtain in-flight performance data on an aerospike rocket engine. The LASRE components mounted to the top of the SR-71 airplane are referred to as the "canoe," "kayak," "reflection plane," and "model" (fig. 4). Collectively, these structural components are referred to as the LASRE "pod." The canoe was installed on the SR-71 fuselage and was designed to contain the gaseous hydrogen fuel and liquid water needed for cooling. The kayak, located beneath the reflection plane and on top of the canoe, set the model incidence angle to 2 ° nosedown to align the lower part of the model with the expected local flow over the top of the SR-71 airplane. The reflection plane was mounted on top of the kayak to help promote uniform flow in the region of the model. The model was designed to approximate a half-span lifting body with a 70-deg swept cylinder leading edge and spherical nose, Liquid oxygen and igniter materials required to operate the rocket engine were stored in the model. The model was vertically mounted so that the angle of sideslip of the SR-71 airplane imparted angle of attack on the model.

With a full load of expendables, the pod weighed approximately 14,140 lb. The SR-71 fuel

distribution system was adjusted to accommodate a maximum of 67,000 lb of fuel for the LASRE flights.

This adjustment was made to prevent overloading the vehicle when carrying the added weight of the LASRE pod. However, actual fuel loads of a maximum of only 62,000 lb were used during the LASRE flight tests. To compensate for CG shifts caused by the pod weight, 5000 lb of available fuel in the forward tank was considered unusable during the flight.

The high transonic drag of the LASRE configuration showcased the ability of the SR-71 airplane to achieve stabilized data at speeds approximating Mach 1.0. Apart from the main purpose of this report, this unique capability of the SR-71 airplane to sustain nearly Mach-I test conditions deserves emphasis.

The SR-71 physical attributes (inertia, fineness ratio, control systems, and the relative characteristics of the transonic drag and propulsive forces) all combine to provide a unique platform for exposing experimentalshapes to selectedstabilizedtransonicflow conditionsin a real flight environment.These conditionscan be maintainedfor severalminutes;andbecause of the relatively large sizeof the aircraft, the candidatemodelscanbeof a respectable scaleandcanincludesignificantdetail.

Figure 5 showsatime history of a nearly level-altitudeacceleration from Mach0.9 to Mach 1.1.This figure showsthe fairly smoothtransition from subsonicto supersonic flight. Much of the acceleration wasachievedusingfull afterburnerthrust.If desired, the pilot could stabilizeat anytransonicspeed using throttle control (with the exceptionof stabilizing at free-stream MachnumbersbetweenMach 1.010and Mach 1.025,which is wherethe airdataMachjump occurs).This capability makesthe SR-71airplanea unique and versatiletransonicresearchfacility that is currently availableto the flight test community.

This researchcapability of the SR-71 airplane represents a valuable complementto its well-known capabilityfor flight research at high supersonic speeds to a maximumof Mach 3.2.

Test Bed Configuration

Four flights were flown with the model removed from the LASRE pod. This configuration became known as the test bed configuration because it can accommodate new model shapes for flight testing.

The weight of the remaining canoe, kayak, and reflection plane is approximately 9400 lb. Fuel loads of a maximum of 66,000 lb were used for these tests.

The test bed configuration also demonstrated its excellent capability for transonic flight research.

Figure 6 shows a time history of a nearly level-altitude acceleration from Mach 0.9 to Mach 1.1. As in figure 5, a smooth transition from subsonic to supersonic flight is seen. The angle-of-attack time history shows two longitudinal PID doublets performed during the acceleration. Throttle control can be used to stabilize at any transonic speed (except at the airdata Mach jump).

Mass Properties

Accurate estimates of weight, CG, and mass moments of inertia were required for each PID

maneuver. Each fuel tank is instrumented to obtain fuel quantity. The total weight is simply the sum of the zero-fuel weight and the total fuel weight recorded by the six tank sensors. Because of the symmetry of the left and right sides of the aircraft, only the rolling, pitching, and yawing moment and cross product of inertias (I x, ly, I z, and Ixz, respectively) were required. Table 2 shows a summary of zero-fuel weight, zero-fuel weight CG, and zero-fuel inertias for the three flight configurations. Figure 7 shows the body-axis mass moments of inertia for the baseline configuration (using the standard fuel burn schedule) as a function of total vehicle weight.

Fuel quantity measurements from the six fuselage fuel tanks are used to compute the CG. Each fuel tank CG is a function of both measured fuel quantity and aircraft pitch attitude. For the LASRE and test bed configurations, pod component CGs are also used to obtain the total configuration CG. Flight data were obtained at CGs ranging from 0.173 c to 0.258 c. Individual fuel quantity measurements and pod component mass distribution information are used to compute the inertias at each test condition.

Table 2. Zero-fuel weight, CG, and mass moment of inertia information for the three SR-71

configurations.

SR-71 zero-fuel weight CG

Fuselage Mean aerodynamic SR-71 zero-fuel inertias,

Configuration Flight SR-71 zero-fuel

station, in. chord, percent slug-ft 2

numbers weight, lb

Ix = 220,660

877.9 20.1 ly = 954,850

Baseline 37-44 60,728

I z = 1,172,039

Ixz = 19,200

Ix = 230,880

911.0 27.4 ly = 1,035,140

LASRE 45-48 74,032

I. = 1,252,330

Ixz = 44,640

I x = 230,880

910.6 27.3 ly = 1,035,520

LASRE 49-51 75,349

/z= 1,252,710

Ix= = 44,710

I x = 224,670

892.3 23.3 ly - 992,280

Test bed 52-55 70,158

I. = 1,209,470

Ix: . = 28,390

METHODS OF ANALYSIS

This section describes the formulation of the output-error parameter estimation technique used to analyze the flight data. The nonlinear equations of motion used in the analysis also are defined.

Parameter Identification Formulation

The primary objective of this research is to estimate from flight test the stability and control derivatives for each of these SR-71 configurations. The actual vehicle system is described by a vector set of dynamic equations of motion that are defined in the next section. The form of these equations is assumed to be known, but the time-invariant aerodynamic stability and control parameters in these equations are unknown. The PID flight test maneuvers are designed to record the response of the aircraft system to measured control inputs. The parameter estimation program known as pEst (ref. 3) is used in postflight analysis to adjust the unknown parameter values in the model until the estimated aircraft response agrees with the measured response.

The pEst program defines a cost function that can be used to quantitatively measure the agreement between the computed response and the actual measured response of the model. The pEst program searches for the unknown parameter values to minimize the cost function.

To obtain the cost function, the pEst program must solve a vector set of time-varying ordinary differential equations of motion. The equations of motion are separated into a continuous-time state equation and a discrete-time response equation:

x(t) = F[x(t), u(t), _] (2)

(3)

Z(ti) = G[x(ti), u(ti), _] where F is the state derivative function, G is the response function, x is the state vector, x is the time derivative of the state vector, z is the response or measurement vector, u is the measured control input vector, _ is the stability and control derivative parameter vector, and t is time. For this application of stability and control derivative estimation, state noise is assumed to not exist.

The output-error cost function, J(_), used by the pEst program is as follows: n t

1 *

(4)

J(_) - 2nzn t E [z(ti)-z(ti)] W[z(ti)-z(ti)] i=1 where n t is the number of time history points used, n z is the number of response variables, _ is the estimated response vector, and W is the response weighting matrix. The superscript * denotes transpose.

For each possible estimate of the unknown parameters, a probability that the aircraft response time histories attain values approximating the observed values can be defined. The maximum-likelihood estimates are defined as those estimates that maximize this probability. Minimizing the cost function gives the maximum-likelihood estimate of the stability and control parameters.

Figure 8 shows the maximum-likelihood parameter estimation process. The measured response is compared with the estimated response, and the difference between these, called the response error, is included in the cost function. The minimization algorithm is used to find the coefficient values that minimize the cost function. Each iteration of this algorithm provides a new estimate of the unknown coefficients on the basis of the response error. These new estimates are then used to update values of the coefficients of the mathematical model, providing a new estimated response and, therefore, a new response error. Updating the mathematical model iteratively continues until a convergence criterion is satisfied (in this case, the ratio of the change in total cost to the total cost, AJ(_)/J(_), must be less than 0.000001). The estimates resulting from this procedure are the maximum-likelihood estimates.

The estimator also provides a measure of the reliability of each estimate based on the information obtained from each dynamic maneuver. This measure of reliability is called the Cram6r-Rao bound (ref. 6). The Cram6r-Rao bound is a measure of relative, not absolute, accuracy. A large Cram6r-Rao bound indicates poor information content in the maneuver for the derivative estimate.

Equations of Motion

The aircraft equations of motion used in the PID analysis are derived from a general system of nine coupled, nonlinear differential equations that describe the aircraft motion (ref. 4). These equations assumea rigid vehicle anda flal, nonrotatingEarth. The time rate of change of mass and inertia is assumed negligible. The SR-71 configurations studied herein, like most aircraft, are basically symmetric about the vertical-centerline plane. This symmetry is used, along with small angle approximations, to separate the equations of motion into two largely independent sets describing the longitudinal and lateral-directional motions of the aircraft. The equations of motion are written in body axes referenced to the CG and include both state and response equations. The applicable equations of motion are as follows for the longitudinal and lateral-directional axes: Longitudinal state equations:

_R

+ (5)

sin0sinc_)

a = -. _ , C_,cosc_ + q -tanl3(pcosc_ + rsinc_)+ V-_2-_ [3 (c°s0c°s0c°sc_L.u_

mvcosp '"

(6)

Oily = ctSCCmR + [rp(l z- I x ) + ( r2- p2)Ixz]/R

(7)

0 = qcosO-rsinO

Longitudinal response equations:

& = ot + xa q (8)

Cl = q + qbias (9)

- 0 (10)

an= m _IS CNg g_[Xa,(l+Ya,P]-g-_Za, (q2+ P2)+anbi"._ (11)

where qbias and anbia s are estimates of instrumentation biases and R is a conversion factor between degrees and radians.

Lateral-directional state equations:

?ISR C gR

= _ y+psino_-rcoso_+-Q-[sinOcosOcos_-sin_(cosOcosOsinct-sinOcosot)] (12)

(13)

pI x - i'Ixz = _ISbCIR + [qr(Iy - Iz) + pqlxz]/R

(14)

i-I z - Pl xz = ?tSbCnR + [pq(l x - ly) - qrlxz]/R

(15)

= p +qtanOsinO+ rtanOcosO

Lateral-directionalresponse equations:

P -- X_ + _bias (16) :13+ z13

E rl

[9 = p + Pbias (17)

= r + rbias (18)

= ¢ (19)

(20)

mg - -_[-x a i" + ZayP] -- + r 2)

av = 77Scv + 1 g_Yay(p2 - y

where _bias' Pbias" and rbias are estimates of instrumentation biases.

coefficients. The

Equations (5)-(20) contain locally linear approximations of the aerodynamic

longitudinal aerodynamic coefficients are expanded as follows: ¢

C N = CNb + CNa(X + -_--V--RCNqq + CNSe_) e (21)

c

C m = Crab + CmaCt + _--V_Cmqq + Cmse_ e (22)

The coefficients are based on a reference area of 1605 ft 2 and a mean aerodynamic chord of 37.7 ft.

The coefficient with the subscript "b" is a linear extrapolation of the angle-of-attack derivative from the average angle of attack of the maneuver to 0 ° angle of attack (ref. 6). Axial force coefficients are not used in this analysis because the engine performance model is not well known, but this is not a concern because the axial force derivatives do not significantly affect flying qualities. All the longitudinal derivatives in equations (21)-(22) are estimated in the analysis.

The lateral-directional aerodynamic coefficients are expanded as follows:

b

(23)

Cy = CYb + Cyf3_ + "_-R( Cy pp + CYrr ) + CYsa_ a + CY_r8 r

b

(24)

C l = Clb + Cl_ + _-Q'-_(CIpP + Clrr) + ClSa_ a + ClSr_ r

b

(25)

C n = Crib + CnB_ + _-V-_(CnpP + Cnrr) + Cnsa6 a + Cnsr8 r

The reference span, b, is 56.7 ft. The coefficient with the subscript "b" is a linear extrapolation of the angle-of-sideslip derivative from the average angle of sideslip of the maneuver to 0 ° angle of sideslip. All the lateral-directional derivatives in equations (23)-(25) were estimated in the analysis.

INSTRUMENTATION AND DATA ACQUISITION

The SR-71 airplane is equipped with a complete set of research airdata and inertial instrumentation.

Free-stream pitot-static airdata are obtained from a calibrated noseboom. Angle-of-attack and -sideslip data are obtained from a four-hole hemispherical probe doglegged to the noseboom. The angle of attack is referenced to the wing reference plane, which is 1.2-deg nosedown in incidence compared to the fuselage centerline reference plane. Angle-of-attack and -sideslip measurements are lagged on the order of 0.2 to 0.4 sec because of the pneumatic plumbing. These lags are accounted for by time skews in the data analysis. Pitch and roll attitude data are obtained from the SR-71 inertial navigation system.

Three-axis angular rate and linear accelerations are measured using strapdown sensors installed at fuselage station (FS) 683.0, butt line (BL) 32.5, and water line (WL) 86.2.

Signal conditioning on the angular rate and acceleration measurements includes a first-order passive antialiasing filter with a 40-Hz rolloff frequency. This filter imparts a 45-deg phase lag that is equivalent to a 3-msec time lag. This time lag is less than the sample time interval of 5 msec. Although measured at sample rates that are higher, the flight data are thinned to 20 samples/sec for the PID analysis. Angles of attack and sideslip and linear accelerations are corrected in the mathematical model of the pEst program to the CG using angular rate measurements and sensor position information. Vehicle weight and longitudinal CG location are obtained using fuel tank measurements. Laterally and vertically, the CG is assumed to be located at BL 0 anti WL 100, respectively. All control-surface positions are measured with the exception of the right outboard elevon. For the PID analysis, only the inboard surface positions are used to define the elevon and aileron deflections.

In the PID analysis, angle of attack and normal acceleration are the primary aircraft responses used to obtain normal force coefficient estimations. Similarly, angle of sideslip and lateral acceleration are used to obtain side force coefficient estimates. Because of uncertainties in the calibration and pneumatic lags associated with the flow angle measurements, larger weights are assigned to the acceleration response measurements in the PID analysis than to the flow angle response measurements. The response weights are constant for all configurations tested, with the exception of a few cases where temporary instrumentation problems existed.

FLIGHT TEST APPROACH

The objective of this research is to obtain baseline SR-71 stability and control derivatives from PID flight data. When the baseline derivatives are known, the aerodynamic effects of the LASRE and test bed configurations can then be obtained from further PID flight testing. An "'envelope expansion" approach has been used to safely "clear" the configuration to the desired maximum Mach number. The envelope was expanded by incrementally increasing Mach number and performing PID doublet maneuvers. The pilots evaluated the airplane handling qualities in real time, and the PID maneuvers were analyzed postflight using the pEst program to obtain stability and control derivative estimates. Therefore, multiple flights were required to clear the envelope for safe operations.

Each PID maneuver consisted of pilot-performed doublet inputs. Figure 9 shows a typical

longitudinal maneuver.As the figure shows,the pilot performsthe pitch doublet by pushingthe stick forward, pulling the stick aft, and then returning the stick to the neutral position. Figure 10 showsa typical lateral-directionalmaneuver.The maneuverconsistsof the pilot performing a yaw doublet with the ruddersimmediatelyfollowed by a roll doubletwith lateralstick movements. The rudderpedalinputs were not measuredor neededfor the PID analysis.The requiredrudder surfacedeflections,however, weremeasured (fig. 10).

For PID analysis,the aircraft response measurement signal-to-noise ratio needsto be largeenoughto obtain good identifiability. Typically, signal-to-noise ratios greaterthan 5 are considered adequate. For SR-71operations,structuredloadsconstraintslimited the sizeof the doubletmaneuvers. For the LASRE and test bed configurations,the constraintswere defined as angular accelerationlimits and a normal acceleration limit. The angularacceleration limits were 43.0deg/sec 2 for roll, 4.5 deg/sec 2 for yaw, and 8.0 deg/sec 2 for pitch (note that angular accelerationswere not measuredon the airplane, but were computedin real time in the control room from angularrate measurements). The normal acceleration constraintrequired the maneuverto be performed between 0.6 and 1.4 g. Because of the initially unknown aerodynamics of the LASRE configuration, two sizes of doublet maneuvers were used for the LASRE configuration testing. At each Mach condition, an initial micro-sized doublet was performed. If the pilots and the control room personnel concurred that a larger doublet would not exceed the constraints, then the pilot would execute a larger, but still small-sized, doublet.

Ideally, the doublets would be performed at a stabilized 1-g flight condition. However, for flight test efficiency reasons, many of the maneuvers were performed during accelerated flight. Typically, PID data for supersonic conditions were obtained during a 450-knots equivalent airspeed (KEAS) climbing acceleration to Mach 2.6 that was followed by a continuous "KEAS bleed" to 390 KEAS as the airplane climbed and accelerated to approximately Mach 3.2. The supersonic deceleration typically consisted of a 365-KEAS descent during which PID data were also obtained. Because of the SR-71 aeroelastic characteristics, obtaining PID data at consistent equivalent airspeeds was desirable. For the highly flexible SR-71 aircraft, the stability and control derivatives are a secondary function of dynamic pressure in many cases. This data set emphasized obtaining supersonic data at two dynamic pressure conditions to show the influence of aircraft flexibility on the derivatives. The 450-KEAS acceleration equated to a dynamic pressure of 686 lbf/ft 2 and the 365-KEAS deceleration equated to a lbf/ft . In some cases, data also were obtained for less than 365 KEAS for dynamic pressure of 451 2 additional envelope clearance.

After each flight test, the PID maneuvers were analyzed to obtain the flight-determined stability and control derivatives. Pilot simulations before the subsequent flight included updated aerodynamic model information from the analysis of the previous flights PID maneuvers.

RESULTS AND DISCUSSION

Extensive flight testing has been completed to assess the changes in aircraft stability and control caused by adding experiments to the top of the SR-71 airplane. Flight data to a maximum speed of Mach

3.00 for the SR-71 baseline configuration, to a maximum speed of Mach 1.75 for the LASRE

configuration, and to a maximum speed of Mach 3.00 for the test bed configuration are presented herein.

To obtain these results, 283 doublet maneuvers were flown and analyzed.

Simulation predictions of the baseline derivatives have been compared with flight data. The

simulation predictions come from a workstation-basedbatch simulation. The aerodynamicmodel

incorporatedinto the simulator c:ame from the SR-71 baselineaerodynamicmodel (ref. 4). Predicted stability andcontrol derivativeswereobtainedby linearizing the aerodynamicmodel at the flight test Mach number,altitude, and massproperty conditions. Simulationpredictions will be plotted only to comparewith the baselineconfigurationflight results.Predictedincrementsto the stability andcontrol derivativescausedby the LASRE pod installation were obtainedin wind-tunnel testsand publishedin reference 7. No wind-tunnelpredictions,however,wereobtainedfor the testbedconfiguration.

The PID analysiswasusedto estimatethe open-loopstability andcontrol derivativesfrom pilot-input doublet maneuvers.Stability augmentation systemswere usedat all times in all axesto increasethe closed-loopstability anddamping.With the SASesremainingon, autopilotswereturnedoff for the axes of interestduring the PID maneuvers. Someof the initial flight test resultsfor the testbedconfiguration havebeenpublishedin reference 2 andfor the baselineandLASRE configurationsin reference 8.

For the LASRE configuration, both micro-sized and larger, but still small-sized,doublets were performed.In many cases,the PID analysisshowsgood fits of the measuredand estimatedresponse parameters for the micro-sizeddoublets.However, the Cram6r-Raoboundswere usually larger for the micro-sizeddoubletsthan for the small-sizeddoubletsbecause of smallersignal-to-noiseratios; andthe parameter estimatesometimes differedfrom multiple estimatesusingthe micro-sizeddoublet.Therefore, only small-sizeddoubletresultsarepresented in this report.

Figures11-13 showthe flight envelopeavailablefor this testing andthe Mach numberand altitude flight conditions usedfor the baseline,LASRE, and test bed configurations,respectively.The flight conditionsreferredto as"Iow-KEAS'"test pointsare shownin the shadedregionsof figures 11-13. The remainingtest points are considered "high-KEAS" test points. The distinction betweenIow-KEAS and high-KEAStest pointsis importantbecause, in somecases, aircraft flexibility affectedthe stability and controlderivatives.Theflexibilily effectswereincludedin the baselineaerodynamic model (ref. 2).

Flight data were obtained at CG values ranging between 0.173 and 0.258 c. All moment derivatives were estimated about the flight CG using the pest program. For presentation in this report, the pitching and yawing moment derivatives were corrected to the moment reference using flight-estimated normal and side force derivatives, respectively (ref. 9). The moment reference is located at 0.25 c (FS 900).

Scatter in derivative estimates could be caused by maneuvers being performed at different weights, angles of attack, trim elevon positions, and bank angles. Slight variations in maneuver sizes, flexibility effects, and not accounting for engine gyroscopic effects (which are assumed negligible) could also result in data scatter. As stated previously, the Cram6r-Rao bounds (ref. 6) are used as a measure of relative, but not absolute, accuracy. Large Cram6r-Rao bounds indicate poor information content in the maneuver for the derivative estimate. The Cram6r-Rao bounds plotted in this report have been multiplied by a factor of five to increase clarity.

Longitudinal Derivatives

Longitudinal stability and control derivatives were determined independently from lateral-directional derivatives. This section presents results from the SR-71 baseline, LASRE, and test bed configurations obtained using longitudinal PID pitch doublet maneuvers. A comparison of the stability and control derivatives obtained from the three configurations also will be shown.

Baseline Configuration

Figure 14 shows time histories from typical subsonic, transonic, and supersonic test points. These time histories include pilot stick inputs, elevon control-surface positions, and aircraft responses for angle of attack, pitch rate, pitch attitude, and normal acceleration. For the response parameters, the solid lines represent measured aircraft responses and the dashed lines represent the responses obtained by integrating the equations of motion using the pEst estimates of the stability and control derivatives. As figure 14 shows, the angle-of-attack response shows the worst fit between measured and pEst-estimated responses. This result is not surprising because in the pEst program, the angle-of-attack measurement is weighted less than the normal acceleration measurement because of high confidence in the normal acceleration measurement as explained in the "Instrumentation and Data Acquisition" section. Figures 15-18 show the baseline longitudinal stability and control derivatives.

Figure 15 shows coefficient of normal force and pitching moment bias (C N and C m ) estimates.

b b The circle symbols represent high-KEAS test points and the cross symbols represent low-KEAS test points. The solid and dashed lines are fairings of the high- and low-KEAS results, respectively. These fairings are based on the authors' interpretation of the trends in the flight estimates. The vertical bars on the plots represent the scaled Cram6r-Rao bounds. The CNb do not show significant differences caused by flexibility. The C m show reduced values at supersonic speeds for the Iow-KEAS test points,

b

especially at approximately Mach 1.8. Note that these bias values are not the traditional normal force and pitching moment coefficients at 0 ° angle of attack with no surface deflections. The parameter estimation program, pEst, uses a maximum-likelihood technique to obtain a linear fit of the data around the trim point. The bias is simply the extrapolation of the linear fit to 0 ° angle of attack. For these test points, the trimmed wing reference plane angle of attack varied between 2.5 ° and 6 ° . The baseline aerodynamic model (ref. 4) contains a reasonably linear normal force coefficient over this angle-of-attack range, but the pitching moment coefficient typically was nonlinear. Also, reference 4 shows that the nonlinear effect of aircraft flexibility on pitching moment becomes increasingly significant at supersonic Mach numbers, which is consistent with the observed change in pitching moment coefficient bias estimates at supersonic speeds caused by aircraft flexibility.

Figure 16 shows the angle-of-attack derivatives. The thin solid lines are the simulation-predicted values for the high-KEAS test points, and the thin dashed lines are the simulation-predicted values for the low-KEAS test points. The flight data show reduced values for derivatives of normal force due to angle of attack, CN, compared to the simulation values. The flight data also show the simulation-predicted flexibility effects (fig. 16(a)).

Figure 16(b) shows the deriwttive of pitching moment due to angle of attack, Cm, _ . The flight data show slightly reduced static stability (that is, less negative) compared with the simulator values, and the flexibility effects at supersonic Mach numbers are not as pronounced in the flight data. Note that C)n,_ is computed about the 0.25-c moment reference point. The positive values at subsonic Mach numbers do not indicate that the airplane was ever flown with negative static margins. For those test points, the CG was significantly forward of the 0.25-c reference location. As figure 16(b) also shows, the static stability tends toward zero as Mach number is increased. This tendency was expected because the SR-71 aircraft is designed to have minimum open-loop static stability at the design cruise Mach number of 3.2 to reduce trim drag (ref. 4). The design Math 3,2 cruise CG location for acceptable stability and trim drag is 0,25 c.

Figure 17 shows the dynamic derivatives. The baseline aerodynamic model (ref. 4) assumes a zero value for the derivative of normal force due to nondimensional pitch rate, . The flight data show CNq CNq generally decreases as Mach number increases (figure 17(a)). The derivative of pitching moment due to nondimensional pitch rate, Cmq, measured in flight was larger (that is, a more negative derivative) than predicted, with the exception of the data recorded at Mach 2.5 (fig. 17(b)). The reason for the variation in the derivative estimates at Mach 2.5 is suspected to be that the two test points with positive values were flown at an altitude of 65,000 ft and the two test points with negative values were flown at an altitude of 60,000 ft. Although the simulator predicted no difference in the Cmq for these two conditions, open-loop damping is known to decrease as altitude increases (ref. 4). Open-loop damping is also expected to be reduced as the Mach number increases toward the Mach-3.2 design condition. The test point at Mach 3.0 was flown at an altitude of 80,000 ft and showed a negative Cmq value that causes the values recorded at Mach 2.5 and an altitude of 65,000 ft to be suspect (although only one test point was obtained at Mach 3.0).

Figure 18 shows the elevator control derivatives. The flight-determined normal force derivative estimates are less than predicted (fig. 18(a)). For supersonic conditions, the flight-determined normal force derivative estimates were almost zero. The flexibility effects predicted by the simulation are evident in the flight data. The flight-estimated near-zero elevator contribution to normal force seemed anomalous to the authors; however, further scrutiny of the analysis consistently confirmed the result that deflecting the elevators had little effect on normal force. An equation-error PID technique (ref. 10) also was used to analyze these data and shows the same result. Figure 18(b) shows the pitching moment effectiveness of the elevator. The flight data and simulator predictions agree well and the flexibility effects are clearly seen, especially at supersonic-condition Mach numbers where a clear distinction exists between low- and high-KEAS test points. This aeroelastic effect was expected because control surfaces on flexible aircraft typically become less effective as dynamic pressure increases.

Linear Aerospike SR-71 Experiment Configuration

Figure 19 shows time histories from typical test points for the LASRE configuration at subsonic, transonic, and supersonic conditions. Good fits were obtained between the measured and pEst-estimated responses; angle-of-attack fit is the most noticeably off. The measured pitch rate (fig. 19(b)) shows a 2-Hz response that is caused by the fuselage first-bending mode. Some maneuvers not shown in figure 19 also show a 2-Hz response in the normal accelerometer output. The pEst implementation used in this research assumes rigid body motion and was expected to identify aircraft flexibility effects that result from changes in dynamic pressure. The pEst program was not expected to identify the high-frequency structural modes because no structural equations of motion are included in the formulation. Therefore, the pEst-estimated responses are not expected to match this observed 2-Hz motion.

Figures 20-23 show the LASRE configuration longitudinal stability and control derivatives. Only a few Iow-KEAS test points were obtained during the LASRE program; these points were at Mach numbers approximating 0.90, 1.20, and 1.45 (fig. 9).

Figures 20(a) and 20(b) show the CNb and Cmb, respectively. A negative CNb is seen across the Mach range tested. The pitching moment coefficient bias is negative (nosedown) for subsonic conditions, shows a large positive (noseup) value for transonic conditions, and then shows a constant positive value for supersonic conditions. Flexibility effects are seen in the CNb and Crab during flight at Mach 1.45.

Figure 21 shows angle-of-attack derivatives. No clear flexibility effects are evident from the flight data.

Figure 22 shows the dynamic derivatives, and figure 23 shows the elevator effectiveness derivatives. As with the baseline configuration, the elevators are consistently more effective at the Iow-KEAS test points (fig. 23(b)).

Test Bed Configuration

Figure 24 shows time histories from typical test points for the test bed configuration at subsonic, transonic, and supersonic conditions. Good fits were obtained between the measured and pEst-estimated responses; angle-of-attack fit again is the most noticeably off. Figure 24(c) shows that the airplane responses at Mach 3.02 were very small, although the stick inputs and elevon deflections were approximately the same magnitude as at low-speed conditions. These small responses resulted in some of the high-speed derivative estimates having larger Cram6r-Rao bounds (because of smaller response measurement signal-to-noise ratios) than low-speed test points.

Figures 25-28 show the test bed configuration longitudinal stability and control derivatives. Figure 13 shows the flight test conditions. As figure 13 shows, the majority of the low-KEAS points are at 365 KEAS and the majority of the high-KEAS points are at 450 KEAS.

Figure 25 shows the CNb and Cmb. Some flexibility effects have been identified. Figure 26 shows the angle-of-attack derivatives. High-KEAS test points at Mach 2.10 and Mach 2.23 did show a significantly reduced Cma for the test bed configuration (fig. 26(b)). These results were obtained on the second test bed flight; the results were not repeated on the third test bed flight when doublets were flown at similar Mach numbers.An aerodynamic explanationfor thesetwo data points likely exists because they wereflown ata slightly highernormalacceleration (approximately1.2 g), which resulted in the trim angle of attack being approximately 1o higher and the elevon trim being 1° more noseup. At these slightly different trim conditions, the shock structure is suspected to be different enough to manifest into a reduced stability. Note that these two test points also show reduced Crab values compared to other high-KEAS test points (fig. 25(b)).

Flexibility effects in the angle-of-attack derivatives generally were small but in the same direction as the baseline configuration. Figure 27 shows the dynamic derivatives, and figure 28 shows the elevator effectiveness derivatives. As with the baseline and LASRE configurations, the elevators are consistently more effective at the low-KEAS test points (fig. 28(b)), demonstrating an effect of aircraft flexibility.

Configuration Comparisons

To better assess the effects of the experimental configurations on the longitudinal aerodynamics, the derivatives from each configuration are plotted together (without Cram6r-Rao bounds for clarity).

Figure 29 shows the test points used for this comparison. The test points used for plotting were chosen so that the derivative plots would not be confusing because of flexibility effects (although a low-KEAS test point is included for the baseline configuration because no high-KEAS data were available at approximately Mach 3.0).

Figure 30 shows the CNb and Crab values. A significant amount of scatter exists in the bias coefficients at speeds less than Mach 0.6. At subsonic conditions, the LASRE configuration generally shows more negative normal force coefficient bias values (fig. 30(a)), which would indicate a down force caused by the LASRE pod. At supersonic conditions, no difference exists in normal force coefficient bias values.

Figure 30(b) shows the C,n b . The LASRE configuration has a more negative (that is, more

nosedown) Crab at subsonic speeds, and a larger Crab at the transonic Mach numbers, compared to the baseline configuration. The test bed configuration does not have as much nosedown bias at subsonic conditions, but has a similar transonic peak, as the LASRE data. The test bed Crab at Mach 1.2-1.5 is significantly less than the biases for the baseline and LASRE configurations. These trends in Crab are consistent with results documented in references 2 and 8 for the flight-measured pitching moment coefficient at 0 ° angle of attack and zero control surface deflections.

At subsonic conditions, the LASRE model is suspected to cause a reduction in the surface pressures on the upper surface of the SR-71 airplane at the aft model location, which results in a more negative Crab. At Mach numbers greater than 1.2, the shock from the LASRE model increased the surface pressures on the aft fuselage, and the shock from the leading edge of the canoe increased the pressures on the fuselage forward of the CG. This increase resulted in no change in Crab at the low, supersonic-condition Mach numbers compared with the baseline values. However, removing the model for the test bed configuration resulted in a significant nosedown Crab increment at Mach 1.2-1.6 because the model was no longer pressurizing the upper fuselage. This effect of the canoe and model on pitching moment was identified in the preflight wind-tunnel tests (ref. 7) and quantified in flight. As discussed in references 2 and 8, these transonic pitching moment changes caused by the LASRE and test bed experiments did not result in violating any safety-of-flight margins in trimming the airplane.

Figure 31 shows the angle-of-attack derivative comparison for the three configurations. As the figure shows, the experiments mounted on the SR-71 airplane caused only minor changes in the CNa and Cmc _.

Although not a flight safety concern, the test bed configuration was less stable in the Mach 1.3-1.6 region. What caused this reduced stability is unknown.

Figure 32 shows the pitch-rate derivative comparison. As seen previously, the Cmq contains a lot of scatter but shows a general reduction in damping (that is, less negative values) as Mach number increases. Also, the pitch damping for the experimental configurations was reduced from the baseline configuration values for Mach 0.9-2.5. The elevon derivative estimates (fig. 33) show excellent agreement for the three configurations with only some scatter at the low subsonic speeds. This agreement demonstrates that no significant effect exists of the upper fuselage-mounted experiments on the pitch control effectiveness of the elevons.

The data in figures 30-33 indicate that the experiments mounted on the back of the SR-71 airplane primarily affect the normal force and pitching moment coefficient bias values. No flight safety or handling quality concerns have been identified because of experiment effects on longitudinal stability, pitch damping, or elevon control effectiveness. Some reduced longitudinal stability was seen for the test bed configuration at Mach 1.3-1.6, and the pitch damping for the LASRE and test bed configurations was reduced at Mach 0.9-2.5.

Lateral-Directional Derivatives

Lateral-directional stability and control derivatives were determined independently from longitudinal derivatives. This section presents results from the SR-71 baseline, LASRE, and test bed configurations obtained using lateral-directional PID doublet maneuvers. A comparison of the stability and control derivatives obtained from the three configurations is also shown.

Baseline Configuration

Figure 34 shows time histories from typical subsonic, transonic, and supersonic test points for the SR-71 baseline configuration. These time histories include pilot lateral stick inputs, rudder and aileron control-surface positions, and airplane responses for angle of sideslip, roll rate, yaw rate, roll attitude, and

2O

lateral acceleration. For the response parameters, the solid lines representmeasuredairplaneresponses, and the dashedlines represent the responses obtainedby integratingthe equationsof motion using the pEst estimates of the stability and control derivatives. Good fits between the measured and pEst-estimated responseswere obtainedfor all of the responseparameters.Figures 35--40 show the lateral-directionalstabilityandcomrolderivativesfor the baseline configuration.

Figure 35 showsthe coefficientsof sideforce, rolling moment,and yawing momentbiases, C r , Clb, and Cnb. The SR-71 aircraft is basically symmetrical about the vertical-centerline plane and at]l maneuvers were trimmed at almost 0 ° angle of sideslip. Because the PID maneuvers were done with only small variations in angle of sideslip about 0 ° and because the aerodynamics are almost linear at these small sideslip angles, the bias values were expected to be approximately zero. For the most part, the flight data did show bias values at approximately zero (fig. 35).

Figure 36 shows the angle-of-sideslip derivatives. The estimated derivative of side force due to sideslip, Cy_ , was typically 25-percent less than predicted (figure 36(a)). The trends with Mach number and the lack of flexibility effect:_ agree with the simulation. The derivative of rolling moment due to sideslip, , agrees well with predictions (figure 36(b)). The is more negative at the low-KEAS test

C lf_ C I[_

points compared to the high-KEAS test points, as was predicted by the simulation and as was expected for a flexible aircraft. The Cnf _ estimated in flight agrees very well with predictions (fig. 36(c)), with the only difference being an improvement in stability (that is, more positive value) at Mach 2. The SR-71 aircraft was designed to have minimum directional static stability at the design cruise speed of Mach 3.2 (ref. 4) to save tail weight and drag.

Figure 37 shows the dynamic derivatives due to roll rate. The derivative of rolling moment due to nondimensional roll rate, is less than predicted and shows an aeroelastic effect (figure 37(b)).

Clp, Figure 38 shows the dynamic derivatives due to yaw rate. The derivative of yawing moment due to nondimensional yaw rate, Cnr, was higher (that is, a more negative value) than predicted for Mach numbers less than 2.5 (fig. 38(,,')). As Mach number increases beyond 2.5, the flight data suggests negative open-loop yaw damping (that is, positive values of Cnr). Reduced yaw damping is seen for some of the low-KEAS test points at supersonic conditions. The combination of low aerodynamic yaw damping and low directional stability (fig. 36(c)) at high-supersonic conditions resulted in the requirement for a yaw-axis SAS for adequate closed-loop handling qualities (refs. 4 and 5).

Figure 39 shows the rudder derivatives. Generally, slightly lower side force and higher rolling moment effectiveness values were seen in flight compared to predictions. The flight data agree well with the predicted yawing moment effectiveness values of the rudders, with the exception of lower flight-determined values at subsonic conditions.

Figure 40 shows the aileron derivatives. The simulation uses zero for the derivative of side force due to aileron. The flight data estimated a small number slightly greater than zero (fig. 40(a)). The flight data agree well with the predicted derivative of rolling moment due to aileron, Cl_ , except at approximately t/ Mach 0.9 where the flight-determined effectiveness is higher than predicted (fig. 40(b)). As predicted and expected for a flexible aircraft, the flight data show increased rolling moment effectiveness for the Iow-KEAS test points. The flight data show significantly lower derivatives of yawing moment due to aileron than predicted (fig. 40(c)).

Linear Aerospike SR-71 Experiment Configuration

Figure 41 shows time histories from typical test points for the LASRE configuration at subsonic, transonic, and supersonic conditions. Good fits between the measured and pEst-estimated responses were obtained for all of the response parameters.

Figures 42-47 show the LASRE configuration lateral-directional stability and control derivatives.

Only three low-KEAS test conditions were available for the LASRE configuration data. This lack of data makes it difficult to assess flexibility effects for this configuration.

Figure 42 shows Cyb, Clb, and Cnb. These biases were expected to be almost zero because of the symmetry of the LASRE pod. The flight data show these to be almost zero except at Mach 0.9 for the high-KEAS test point and the Clb values at subsonic conditions. The LASRE configuration would normally trim at approximately 0 ° angle of sideslip. At the high-KEAS Mach-0.9 test points, the configuration trimmed at -0.5 ° angle of sideslip, which was a result of using the rudders to trim out the large positive Crib shown in figure 42(c). What nonsymmetry in the geometry or aerodynamic flow caused this large Cnb value is uncertain.

The Iow-KEAS, Mach-0.9 test point did not show this large Crib value and trimmed at approximately 0 ° angle of sideslip. In this case, the low- and high-KEAS test points were only separated by 26 kn. The high-KEAS test point was flown at an altitude of 27,300 ft (349 KEAS), and the Iow-KEAS test point was only 3,800 ft higher at an altitude of 31,100 ft (323 KEAS). That such a large change in Crib and trim angle of sideslip would occur over such a small change in equivalent airspeeds is interesting. The Clb value at Mach 0.4 is also of interest (figure 42(b)). For this test point, the trim angle of sideslip was 0.5 °, which was the result of trimming out the positive Clb.

Figure 43 shows the angle-of-sideslip derivatives. The Clf _ becomes less negative as the Mach number increases throughout the subsonic range, and then becomes fairly constant for supersonic conditions (fig. 43(b)). Figure 43(c) shows that Cn_ decreases slightly with Mach number and remains stable (that is, positive). Figure 43(c) also shows an increase in Cnf _ for the high-KEAS test points at Mach 0.9, at which the vehicle trimmed at -0.5 ° angle of sideslip as stated above.

Figures 44 and 45 show the dynamic derivatives due to roll rate and yaw rate, respectively. The CI is constant throughout the Mach range studied (fig. 44(b)). The Cnr starts at approximately zero and the_ increases (that is, becomes more negative) as Mach number increases (fig. 45(c)). A significant increase in yaw dampingoccursas the airplanetransitionsfrom subsonicto supersonicflight. Also, aeroelastic effectscaused a reductionin yaw dampingfor the low-KEAStestpointsat supersonic conditions.

Figure 46 shows the rudder derivatives. The most variation in the derivatives occurs between

Mach0.9 andMach 1.0.The ya_ing momenteffectiveness of the ruddersremainsconstantthroughout

the Mach rangeexceptat approximatelyMach 0.9, at which the effectivenessis increased(fig. 46(c)).

Figure47 showsthe aileronderivatives.The rolling momenteffectiveness of the aileronsdecreases with increasingMach number (figure 47(b)). As expectedfor a flexible aircraft, the Iow-KEAS test points indicatehigher rolling momenteffectivenessvaluesthan existedfor the high-KEAStest points. Figure 47(c) showsproverseyaw throughoutthe Mach rangewith the curiousexceptionat Mach 1.05.

Test Bed Configuration

Figure 48 shows time histories from typical test points for the test bed configuration at subsonic, transonic, and supersonic conditions. Good fits between the measured and pEst-estimated responses were

obtained for all of the response parameters. Figures 49-54 show the test bed configuration

lateral-directional stability and control derivatives.

Figure 49 shows Crb, Clb , and Cnb. The values in general are approximately zero for all three derivatives, with the exception of slightly positive Clb values at subsonic conditions. Figure 50 shows the angle-of-sideslip derivatives. The Clf _ decreases toward zero between Mach 2.0 and 3.0 as seen in figure 50(b). The C% values (fig. 50(c)) show a sharp decrease beginning at Mach 2.2, level off at slightly positive values at Mach 2.5, and remain constant to Mach 3.0. Aeroelastic effects are seen in the Clf _ and Cnl 3 data at speeds less than Mach 2.

Figures 51 and 52 show the dynamic derivatives due to roll rate and yaw rate, respectively. The Cip values are negative and fairly constant to a maximum of Mach 2.5, at which point they becomes less negative and then slightly positive at Mach 3.0 (fig. 51(b)). The Cnr values have moderate scatter but remain negative throughout the Mach range (fig. 52(c)).

Figure 53 shows the rudder derivatives. The yawing effectiveness of the rudders reaches its maximum at almost Mach 1, then decreases as Mach number increases to Mach 3.0. Some aeroelastic effects are shown in the rolling and yawing moment rudder derivatives. Figure 54 shows the aileron derivatives. The Cl_ decreases with increasing Mach number (fig. 54(b)). As with the baseline and • , a LASRE configurations, the Iow-KEAS test points indicate higher aileron rolling moment effectiveness values than those for the high-KEAS test points. Figure 54(c) shows mostly proverse yaw at the low Mach numbers; slight adverse yaw begins at Mach 2.3.

Configuration Comparisons

To assess the effects of the experimental configurations on the lateral-directional aerodynamics, the derivatives from each configuration are plotted together (with Cram6r-Rao bounds not included for clarity). Figure 55 shows the test points used for this comparison. The test points used for plotting were chosen so that the derivative plots would not be confusing because of flexibility effects (although a Iow-KEAS test point is included for the baseline configuration because no high-KEAS data were available at approximately Mach 3.0).

Figure 56 shows a comparison of the lateral-directional coefficient bias values. No significant differences in the biases were expected because the LASRE and test bed configurations are symmetric about the vertical-centerline plane. Differences did occur, however, at subsonic-condition Mach numbers and particularly at Mach 0.9 for the Clb and Crib (figs. 56 (b) and (c)). These differences are caused by some unknown geometric or aerodynamic asymmetry.

Figure 57 shows the angle-of-sideslip derivatives. The CIf 3 is reduced for the LASRE and test bed configurations compared to the baseline configuration (fig. 57(b)). Why Clf _ would decrease for configurations with the experiments on top is not intuitive; especially with the LASRE configuration, which essentially includes an additional vertical surface. Wind-tunnel tests (ref. 8) did predict the reduction in Clb for the LASRE configuration. Because the test bed configuration saw a similar reduction in Clb, the canoe is suspected to be the significant contributor to the reduction. The test bed dihedral effect values approach zero as Mach number increases.

The C was also significantly affected by the LASRE and test bed configurations. At transonic nl_ conditions, both experimental configurations show a reduction in Cnl 3 compared to the baseline configuration. This reduction again was not intuitive; especially for the LASRE configuration, which includes the "vertical-like" surface. The test bed configuration directional stability decreased significantly beginning at Mach 2.2 and leveled out at Mach 2.5 with a slightly positive value (fig. 57(c)).

This concern was significant during the envelope expansion phase of the test bed program.

Piloted simulations were done to determine the effect of reduced open-loop stability derivatives on handling qualities and aircraft responses caused by failures of a single engine (for example, inlet unstarts). With reduced directional stability, an engine failure could result in significant sideslip leading to aircraft structural failure. Simulations were done to demonstrate the open-loop directional stability level at which point an engine failure would result in exceeding the aircraft sideslip limit. With the additional closed-loop stability provided by the yaw SAS (eq. (1)), the determination was made that flight tests could safely proceed to Mach 3.2. A new mission rule was created, however, requiring the CG to be forward of 23 percent for flight at speeds faster than Mach 2.5 (ref. 2).

Figure 58 showsa comparisonof the roll-rate dynamic derivatives. For the most part, significant differences in the data scatter are difficult to see. Some evidence exists that the Ctp decreases (that is, the derivative becomes less negative) at speeds faster than Mach 2.5 for the test bed configuration, and even shows positive values at speeds faster than Mach 2.9 (fig. 58(b)). Figure 59 shows a comparison of the yaw-rate dynamic derivatives. The Cnr (fig. 59(c)) tends to indicate that yaw damping for the test bed configuration increases at speeds faster than Mach 2.5 compared to the baseline configuration.

Figure 60 shows the rudder control derivatives. The LASRE configuration shows reduced side

force, rolling moment, and yawing moment derivatives at high subsonic-condition and low

supersonic-condition Mach numbers compared to the baseline configuration. This reduction is

especially true at Mach 0.9, a! which Crib and Ctb were known to be affected by the LASRE

experiment. At speeds greater than Mach 1.5, no difference exists in rudder control effectiveness for the three configurations.

Figure 61 shows the aileron control derivatives. The aileron derivative of primary interest is the Cl, _ (fig. 61(b)). As can be seen, the LASRE and test bed configurations had little to no effect on the U rolling moment effectiveness of the ailerons. The LASRE model, however, did have large effects on the side force and yawing moment derivatives that show significantly more proverse yaw at speeds less than Mach 1.8, except at approximately Mach 1.05.

The data in figures 56--61 indicate that the experiments mounted on the back of the SR-71 airplane did affect lateral-directional stability and control. The primary concern for the project was the reduced static stability and reduced dihedral effects at high Mach numbers for the test bed configuration. The LASRE configuration also demonstrated reduced directional control effectiveness at transonic-condition Mach numbers, but this was not a concern to the project because the SR-71 tails were sized for the Mach-3.2 design cruise condition and thus were very effective at the transonic speeds.

CONCLUDING REMARKS

Longitudinal and lateral-directional stability and control derivatives were obtained from flight test of the SR-71 baseline aircraft and fi'om SR-71 configurations with the Linear Aerospike SR-71 Experiment

(LASRE) and test bed experiment hardware mounted on the SR-71 upper fuselage. A

maximum-likelihood output-error parameter estimation program was used to analyze a total of 283 pilot-input doublet maneuvers. This work has resulted in a significant database of parameter identification maneuvers for a large, flexible aircraft over a range from Mach 0.4 to Mach 3.0. Many test points were obtained at transonic speeds because of the ability of the SR-71 aircraft to stabilize at flight conditions at almost Mach 1.O Also, the aeroelastic characteristics of the aircraft were studied by obtaining flight data at different dynamic pressures. Longitudinal and lateral-directional stability and control flight tests of the three SR-71 configurations showed the results detailed in the following sections.

Longitudinal Derivatives

The SR-71 baseline flight estimates of the longitudinal stability and control derivatives were compared with the simulation predictions. The flight-determined normal force derivative caused by angle of attack was less than predicted, but did show the predicted flexibility effects. The flight data showed slightly reduced derivatives of pitching moment due to angle of attack, Cmc _, compared to predictions and did not show significant variation in static stability caused by flexibility. The flight data showed larger derivative values of pitching moment due to nondimensional pitch rate, Cmo, than predicted. The derivative of normal force due to elevon was significantly less than predicted, but showed predicted flexibility effects. The derivative of pitching moment due to elevon, Cm_ , agreed e well with predictions and showed increased control effectiveness at low dynamic pressures, as expected for a flexible wing.

The LASRE and test bed configuration flight data were compared with the baseline configuration data to assess the effects of the experimental hardware on the longitudinal stability and control derivatives. The results obtained are as follows:

• The LASRE experiment induced a more nosedown pitching moment bias at subsonic-condition

Mach numbers compared to the baseline and test bed configurations. Both the LASRE and test

bed configurations showed larger noseup values at approximately Mach 1. The test bed

configuration showed a large nosedown increment at low supersonic Mach numbers.

• Mounting of the LASRE and test bed hardware did not significantly affect the open-loop Cma.

Some reduced stability was seen for the test bed configuration at Mach 1.3-1.6.

• The Cmq for the LASRE and test bed configurations was reduced at Mach 0.9-2.5.

• The Cm_ showed excellent agreement for the three configurations and had only some scatter at e low subsonic speeds.

In general, the experimental hardware of the LASRE and test bed configurations resulted in no significant stability, control, trim, handling quality, or safety-of-flight concerns in the longitudinal axis.

Lateral-Directional Derivatives

The SR-71 baseline flight estimates of the lateral-directional stability and control derivatives were compared with the simulation predictions. The flight-determined side force derivatives caused by angle of sideslip were approximately 25-percent lower than predicted at all Mach numbers. The derivative of rolling moment due to sideslip, Clf _, agreed well with predictions, including the prediction of increased dihedral effect at low dynamic pressures. The derivative of yawing moment due to sideslip, C,_f, agreed well with predictions except at Mach 2, at which flight data showed higher stability than predicted. As predicted in the simulation, no change in stability caused by flexibility effects existed. The flight data showed less roll damping than predicted and showed some aeroelastic effects. The yaw damping was higher than predicted for speeds less than Mach 2.5. As Mach number increased beyond Mach 2.5, the flight data showed negative open-loop yaw damping.

The flight data agreed well with predictions of the yawing moment due to rudder except at subsonic speeds, at which the flight data showed less effectiveness than predicted. The flight data agreed well with the predicted derivative of rolling moment due to aileron, CI_ , except at approximately Mach 0.9, at which the flight-determined effectiveness was higher than _predicted. As predicted and expected for a flexible aircraft, the flight data showed increased rolling moment effectiveness for the low-knots-equivalent-airspeed test points. The flight data also showed significantly lower yawing moment due to aileron than predicted.

The LASRE and test bed configuration flight data were compared with the baseline configuration data to assess the effects of the experimental hardware on the lateral-directional stability and control derivatives. The results are as follows:

The LASRE rolling moment and yawing moment coefficient bias values were significantly

different than the test bed and baseline configuration values at Mach 0.9.

The LASRE and test bed configurations demonstrated a reduced Clf _ compared to the baseline configuration values.

The LASRE and test bed Cnl _ showed reductions at transonic flight conditions compared to the baseline configuration. A1 approximately Mach 2.2, the test bed open-loop directional stability

began decreasing rapidly toward zero. At Mach 2.5, the stability leveled off and remained

relatively constant to Mach 3.0. This reduced stability was of great concern during the envelope expansion phase; howeven simulations showed enough margins to fly safely at these high Mach numbers with the test bed configuration and with stability augmentation systems fully operational.

The LASRE model had a significant effect on the rudder derivatives at transonic speeds. The side

force, rolling moment, and yawing moment derivatives caused by rudder deflection were all

smaller than the test bed and baseline derivatives at transonic speeds. In addition, at Mach 0.9, the LASRE rolling and yawing moment rudder derivatives showed significant negative peaks caused by some unknown transonic interaction. None of these rudder derivative variations resulted in safety-of-flight concerns.

The Ct_ was unaffected by the LASRE and test bed configurations.

a

The LASRE model did have a large effect on the sideforce and yawing moment aileron

derivatives. The LASRE data showed significantly more proverse yaw than the baseline and test bed configurations, except at approximately Mach 1.05.

The most significant lateral-directionaleffect of the experimentalhardwarewas the reduction of directional stability of the test bedconfiguration at speeds greaterthan Mach 2.2. Flight teststo Mach 3.0 and simulation of emergencysituationsshowedthat even with the reducedstability, the test bed configuration could be safely flown. In general,the LASRE and test bed experimentalconfigurations showedacceptable stability,control, trim, andhandlingqualitiesthroughoutthe Machregimestested.

As a concluding remark, future potential SR-71 "piggyback" flight test experimentscan use the results presentedin this report to "size" and locateexperimentalhardwareon the back of the SR-71

airplane to ensure adequate stability and control capability. For example, any reasonably

aerodynamicallyshapedhardwaremountedat the aft end of the canoeshould improve the directional stability abovethatseenfor the testbedconfiguration.

REFERENCES

, Corda, Stephen, Bradford A. Neal, Timothy R. Moes, Timothy H. Cox, Richard C. Monaghan, Leonard S. Voelker, Griffin P. Corpening, Richard R. Larson, and Bruce G. Powers, Flight Testing the Linear Aerospike SR-71 Experiment (LASRE), NASA TM-1998-206567, 1998.

, Corda, Stephen, Timothy R. Moes, Masashi Mizukami, Neal E. Hass, Daniel Jones, Richard C.

Monaghan, Ronald J. Ray, Michele L. Jarvis, and Nathan Palumbo, The SR-71 Test Bed Aircraft: A Facility for High-Speed Flight Research, NASA TP-2000-209023, 2000.

3. Murray, James E. and Richard E. Maine, pEst Version 2.1 User's Manual, NASA TM-88280, 1987.

4. Lockheed Aircraft Corporation, Handling Qualities of the SR-71, SP-508, Oct. 1964.

5. DeGrey, R. P., R. L. Nelson, and J. E. Meyer, "SR-71 Digital Automatic Flight and Inlet Control System," SAE-851977, Oct. 1985.

6. Maine, Richard E. and Kenneth W. Iliff, Application of Parameter Estimation to Aircraft Stability and Control: The Output-Error Approach, NASA RP-1168, 1986.

.

Moes, Timothy R., Brent R. Cobleigh, Timothy R. Conners, Timothy H. Cox, Stephen C. Smith, and

Norm Shirakata, Wind-Tunnel Development of an SR-71 Aerospike Rocket Flight Test

Configuration, NASA TM-4749, 1996 (also published as AIAA-96-2409).

.

Moes, Timothy R., Brent R. Cobleigh, Timothy H. Cox, Timothy R. Conners, Kenneth W. Iliff, and Bruce G. Powers, Flight Stability and Control and Performance Results from the Linear Aerospike SR-71 Experiment (LASRE), NASA TM-1998-206565, 1998 (also published as AIAA-98-4340).

.

Gainer, Thomas G. and Sherwood Hoffman, Summary of Transformation Equations and Equations

of Motion Used in Free Flight and Wind Tunnel Data Reduction and Analysis, NASA SP-3070, 1972.

10. Morelli, Eugene A., "Real-Time Parameter Estimation in the Frequency Domain," AIAA-99-4043, Aug. 1999.

EC91 520-7 Figure 1.SR-71baselineconfigurationin flight.

!i_i_ ¸.... ! if!¸ i i

EC98 44509-7 Figure2. LASRE configurationin flight.

EC99-45065-6 Figure 3.Testbedconfigurationin flight.

1023.0 FS 1196.5 FS 1041.9 1230.9 F WL 131.0 Kayak

f

Moment reference FS 1244.6 _S 102.0 FS 1355.0 _BL 340.2 BL165 _BL45.0 II --''-_._-%3BL;6.5 Canoe ........

_'/_ 010252 Figure 4. Side and planform views of the LASRE configuration.

Mach number Pressure altitude, ff KEAS !

I I I Angle of attack, deg 0 100 200 300 400 500 600 700 Time, sec 010253

Figure 5. Time history of an acceleration from subsonic to supersonic flight (LASRE

configuration).

---Airdata Mach jump 1.1 I I Mach 1,0 ....................................

number f ..............................

I I I I 2.56 x 104 I Pressure 2.54 altitude, ft 2.52 I 2.50 I I 45O I KEAS I 400 _'-_- " " _ I i I i Angle of attack, 5 ........................ _--................................. ! ..............................

deg 4 I 0 20 40 60 80 100 120 140 Time, sec 010254

Figure 6. Time history of an acceleration from subsonic to supersonic flight (test bed

configuration).

2,2 x 10 6 2.0 1.8 1.6 1.4 Mass 1.2 moment of inertia, slug.ft 2 1.0 .8 .6 .4 .2 .7 .8 .9 1.0 1.I 1.2 1.3 1.4 x 10 5 .6 Total welght, Ib 010255 Figure 7. Body-axis mass moments of inertia for the baseline configuration using the standard fuel burn schedule.

Con,_, . . _ ...ur.m.°t

'n_"___ T., l__.J__ no,.e

I l a'_rs"l - _J _-.su_

I " _ response L__.I Mathematical mode I _+

- I of aircraft _

_1 (a_te es.mator) _st,mat_-- )

response Minimization likelihood [_ Maximum- functional

method cost I--.ee_r.S e

estimate of aircraft _ Maximum-Ilkellho_ _ parameters ooo277

Figure 8. Maximum-likelihood parameter

estimation process.

_r, deg -5 11 a, in.

-2 -4 _e, in.

_a,

2 I

deg -2 --6 _e, d_l -2 deg -4 deg deg -10 ql P_ deg/sec 0 deg/sec -2 -10 -4 rl O_ deg/sec deg .1 ay, 1.0 g an'g "" 1.5 [

I I

--.I .5 0 2 4 6 8 10 2 4 6 Time, sec OlO2S8 Time, sec 010257

Figure 10. Typical yaw-roll doublet time

Figure 9. Typical pitch doublet time history

history (Mach 2.79 and an altitude of

(Mach 1.06 and an altitude of :27,900 ft).

79,800 ft).

90,000 Airspeed, KEAS 80,000 350

t

60,000 / Low-KEAS //"//J 50,000 Pressure 50,000 altitude, ft 40,000 30,000 20,000

L I w/// 0 Pitch doublets

10,000 | / _//I I-1 Yaw-roll doublets

/I _,';',/,' I I J I I I

0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010259

Figure

11. Flight conditions for stability and control test points

(baseline configuration).

90,000 -- Airspeed, KEAS 80,000 70,000 60,000 -- _ 350 -- Low-KEAS S test points _ _//, Pressure 50,000 altitude, ft 40,000

-

30,000 20,000 //pJ_l,_//_//// 0 Pitch doublets 10,000 [] Yaw-roll doublets

l _Vl l/l

I J/;/,'t I I I i I

.5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010260

Figure 12. Flight conditions for stability and control test points

(LASRE configuration).

Airspeed, KEAS 90,000 35O 80,000 .400 o 450 70,000 60,000 Low-KEAS Pressure 50,000 test altitude, ft 40,000 30,000 / 20,000 0 Pitch doublets I 0,000 I-I Yaw-roll doublets

I I I I I

2.5 3.0 3.5 0 .5 1.0 1.5 2.0 Mach number 010261

Figure 13. Flight conditions for stability and control test points (test

bed configuralion).

Measured value -- Measured value -- Measured value pest estimate .... pEst estimate .... pEst estimate 11e, 11e, _e, in. In.

in.

-3 -3 _e' (_e, _e, -2 deg

deg

d_ -2 -4 -2 --4 -6 1% CCj deg deg d_i 1 1 q9 q9 deg/sec deg/sec -2 -3 -4 -3 0_ deg deg deg 2.0 1.5 1.5 1.5 a n , aN , an , g g g 1.0 1.0

.5 .5 I I

.5 0 2 4 6 0 2 4 6 0 2 4 6 8 Time, sec Time, sec Time, sec 010311 010312 010262

(a) Mach 0.6 and an (b) Mach 1.06 and an

(c) Mach 2.99 and an

altitude of 6000 ft.

altitude of 27,900 ft.

altitude of 80,100 ft.

Figure 14. Baseline configuration longitudinal maneuver time histories.

-.10 Test points High-KEAS -- -x--- Low-KEAS -.05 -- CN b -.05 -- -.10

I

I I [ I I I

-.15 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010313

(a) C

N o " .020 Test points ----(P-- High-KEAS -- .-N--- Low-KEAS ,015 .010 -- .005 -- Cm b 0-- -.005 --

-.o,o ] ] 1 I I I I

0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010263 (b) Crab.

Figure 15. Flight-determined longitudinal coefficient biases (baseline configuration).

.06 Test points High-KEAS -,-x--. Low-KEAS .05 High-KEAS (predicted) Low-KEAS (predicted) .04 CNc(' .03 deg -1 .02 .01

I I I I I I I

.5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010314

(a) C

Ne_" 4 x 10 -3 Test points High-KEAS 3-- ---'*-- = Low-KEAS 2-- High-KEAS (predicted) .... Low-KEAS (predicted) I -- 0 -- Cm(_' -1 -- --2 -- --3 -- --4 --

I I I I I I I

-6 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 01O264 (b) Cma (moment reference at 0.25 c).

Figure 16. Predicted and flight-determined angle-of-attack derivatives

(baseline configuration).

-- Test points •---e,--- High-KEAS -5 Low-KEAS High-KEAS (predicted) Low-KEAS (predicted) -10 -15 CNq ' -20 rad-1 -25 -30 .L "" -35

I I I I ] "'"--/ I

-40 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010315 (a) C N q .

1.0 -- .5 _ 0 m ".5 -1.0 -- Cmq' -1.5 -- rad -1 -2.0 -- $ -2.5 -- -it Test points High-KEAS -3.0 -- " -- _- Low-KEAS High-KEAS (predicted) -3.5 Low-KEAS (predicted)

( 1 I I I I I

-4.0 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010265 (b) Cmq (moment reference at 0.25 c).

Figure 17. Predicted and flight-determined pitch-rate derivatives

(baseline configuration).

.015 Test points High-KEAS -- "*=- - Low-KEAS High-KEAS (predicted) .010 Low-KEAS (predicted)

, \

.005

Ii

CNSe' deg -1 i i i i _

I I

I I I I I

-.010 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010316

(a) CNa

¢ 0 x 10-3 -1 -- -2 -- -3 -- Cm_e' deg -1 -4 -- Test points -5-- High-KEAS ..-m-. Low-KEAS -6-- High-KEAS (predicted) ..... Low-KEAS (predicted)

I I I I I r I

-7 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010266 (b) Cma (moment reference at 0.25 c).

e

Figure 18. Predicted and flight-determined elevator derivatives

(baseline configuration).

Measured value Measured value Measured value pEst estimate .... pEst estimate .... pEst estimate _e, 11 e , 11 e, 0 - : in.

in.

in.

-2 -2 _e, (_e, (_e, deg deg deg -2 -4 (_, deg deg deg q_ q_ q) deg/sec deg/sec deg/sec -3 -2 -3 6.5 0_ , O, deg deg deg 5.5 6.0 5.0 1.4 1.4 1,4 a n , a n , a n , 1.0 g g g 1.0 1.0 .6 .6 .6 0 2 4 6 8 2 4 2 4 Time, see Time, sec Time, sec 010317 010318 010267

(a) Mach 0.5 and an

(b) Mach 1.03 and an

(c) Mach 1.73 and an

altitude of 12,000 ft.

altitude of 26,800 ft.

altitude of 43,600 ft.

Figure 19. LASRE configuration longitudinal maneuver time histories.

-.10 Test points High-KEAS ----_--- Low-KEAS -.05 CN b -.10 --

I

I I I I I I I I I

-.15 0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 Mach number 010319 (a) CNb.

.020 Test points High-KEAS ----_--- Low-KEAS .015 .010 -- .005 -- Cm b 0-- -.005 --

I I

-.010 I I I I I I I I

0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 Mach number 010268 (b) Cmb.

Figure 20. Flight-determined longitudinal coefficient biases (LASRE configuration).

.06 Test points _ High-KEAS .05 -- _,, --.-K--- Low-KEAS .04 -- CNc_' .03 -- deg-1 .02 -- .01 --

I r,,. I I I I I I I I

.2 .4 .6 ,8 1.0 1.2 1.4 1.6 1.8 2.0 Mach number 010320 (a) CNe L.

4x 10 -3 Test points -----0--- High-KEAS 3-- "_ --..x--- Low-KEAS 2-- 1 -- 0-- Cm_' -1 -- deg -1 -2 -- -3 -- -4 -- -5 --

I I I 1 I I 1 I I I

-6 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 Mach number 010269 (moment reference at 0.25 c).

(b) Cma

Figure 21. Flight-determined angle-of-attack derivatives (LASRE configuration).

Test points High-KEAS -5 ---M--- Low-KEAS -10 -15 CNq' -20 rad-1 -25 -30 -35

I I I I I I I I I I

-4O 0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 Mach number 010321 (a) CNq.

1.0 B Test points ---e-- Hlgh-KEAS °5 -- ---_--- Low-KEAS 0 m -,5 -1.0 m

xl

Cmq' -1.5 -- rad-1 -2.0 -- -2.5 -- -3.0-- -3.5 --

I I I I I I I I I I

-4.0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 Mach number 010270 (b) Cmq (moment reference at 0.25 c).

Figure 22. Flight-determined pitch rate derivatives (LASRE configuration).

.015 Test points High-KEAS -- --_--- Low-KEAS .010 .005

I i I I I I 1 I 1 I

-,010 0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 Mach number 010322

(a) CN_

e 0:10 .-3 -- Test points High-KEAS __ ---it--- Low-KEAS -1 -2 -3 / / / / I -4 / / / -6 --

I I I I I I I I I I

--7 0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1,8 2.0 Math number 010271 (b) Cm_ (moment reference at 0.25 c).

e Figure 23. Flight-determined elevator derivatives (LASRE configuration).

Measured value _- Measured value Measured value .... pEst estimate .... pEst estimate pEst estimate l_et 1"1 e, 1"1 e , in.

in.

in.

-1 -1 -2 _e, _e, _e, deO deg deg 1 -2 I I -2 I I -4 4.5 4.0 C_, 3.5 (_' 3.5 deg deg deg 4.0 3.0 I I 2.5 3.0 q, q_ deg/sec deoJsec -2 -2 5.5 2.0 O_ 0_ O_ 4.5 deg deg deg 5.0 4.0 3.5 1.5 1.1

'21 1

1.0 a n , a n , a n , .9 g g g 1.0 .8 .7 .9 0 1 2 3 0 2 4 0 1 2 3 Time, sec Time, sec 010323 Time, sec 010324 010272 (a) Mach 0.67 and an (b) Mach 0.97 and an (c) Mach 3.02 and an altitude of 7,100 ft. altitude of 25,100 ft.

altitude of 67,900 ft.

Figure 24. Test bed configuration longitudinal maneuver time histories.

.10 Test points High-KEAS ----_--- Low-KEAS ,05 CN b m -.05 -.10

i

I I I I I I

-.15 1.0 1.5 2.0 2.5 3.0 3.5 0 .5 Mach number 010325 (a) CNb.

.020 -- _ Test points .015 .010 Cm b .005 --.005

I I I I I I I

-.010 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010273

(b) C

m b Figure 25. Flight-determined longitudinal coefficient biases (test bed configuration).

.06 l / Test points

T T ¢ _ H,g. KE,S

.05 .O4

iI"

CNc(' .03 deg -1 .02 .01

] I I I I I I

.5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010326 (a) CN .

Test points HIgh-KEAS --.-x-- - Low-KEAS I

I I I

1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010274 (b) Cma (moment reference at 0.25 c).

Figure 26. Flight-determined angle-of-attack derivatives (test bed configuration).

0 -- Test points ---e--- High-KEAS -5 -- ..x--- Low-KEAS -10 -- -15 -- CNq' -20 -- _d-1 -25 --

I I I I I I I

-4O 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010327

(a) C

Nq" 1.0 _ Test points ---e--- HIgh-KEAS • 5 _ ----K--- Low.KEAS 0-- -1.0 -- Cmq' -1.5 -- rad-1 -2.0-- -2.5 -- -3.0 -- -3.5 --

I It- I I I I I

-4.0 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010_75 (b) Cmq (moment reference at 0.25 c).

Figure 27. Flight-determined pitch-rate derivatives (test bed configuration).

.015 "T" Test points _ High-KEAS .010 _L.I_T - .-_-- Low-KEAS .005 qu

t

I I I I ( I I

-.010 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010328

(a) CN_

e Cm(_e' f deg-1 I Test points Hlgh-KEAS ---K--- Low-KEAS

I ] I I I I

1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010276 (b) Cm_ (moment reference at 0.25 c).

e Figure 28. Flight-determined elevator derivatives (test bed configuration).

90,000 -- Airspeed, KEAS 80,000 45O 70,000 60,000 Pressure 50,000 altitude, ft 40,000 30,000 I ,,,/_ configuration 20,000 I aJ_ /J / O Baseline I.__E]_/IZ' X LASRE 10,000

i .___./_ SR-71 I._,'g/ ,_ +,,,bed

I-_';/,'1 i I t i I

0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010277

Figure 29. [:light conditions for comparison of configuration

longitudinal stability and control derivatives.

.10 -- O .05 CN b

_.._N _• "gO , x

-.05

, ¢

X x x tO SR-71 configuration -.10 o _ Baseline --m-- LASRE - - O- - - Test bed

-.15 I I I J iI I

0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010329 (a) CNb.

.O20 SR-71 ,_ configuration _' _111 + Baseline .015 -- :l , I --_-- LASRE

, _.,, --_-- Testl_,_

.010 -- .005 -- Cm b

a !_ I= _ ..... --_'_"" D a oa

GI _ '.-a--"_- -_O,: ......... -J#,a-a - a .... '_z'-fa o a ._Qm, / 13 0-- x t O 0 x • X_.y__._

_° m

-.005 --

I I I I I I I

--.010 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010278 (b) Crab.

Figure 30. Flight-determined longitudinal coefficient biases (all three configurations).

.06 SR-71 configuration Baseline LASRE .05 Test bed .04 CNc_' .03 d_ -1 .02 .01

[ f f f f f I

.5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010330 (a) CNe L .

4 x tO -3 SR-71 configuration Baseline LASRE 2 ][ --0--- Test bed ,A \'_ a a ""'o \_ o a ota,a ,._a .....

o, 't _o o _._'"'_"

Cm(_' -1 deg -1

o ;__ °o ..°'-'_'_

-2

• _'_.i" --_-7 _'-

-3

o_

--4-- ---5--

... I I I I I I I

0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number olo279 (b) Cmc t (moment reference at 0.25 c).

Figure 31. Flight-determined angle-of-attack derivatives (all three configurations).

SR-71 configuration Baseline -5 LASRE Test bed i.O --0--- -10 -15 CNq, '_113 0 "''° 13 0 _-]_ 0 Ck'"-._ rad-1 -20

x "_.o_f<_

-25 -30 -35

o

I I I I 1 g I

-40 .5 1,0 1.5 2.0 2.5 3.0 3.5 Mach number 010331 (a) CN q.

1.0 D SR-71 configuration • 5 m ---0--- Baseline LASRE Test bed --.5 m:°.

[]g

0 oO -1.0

a x a ..... _.- _ _)/---_.__.._.. -

Cmq' -1.5

o o [] . P'_.'-'_:;,'" oY [] []

rad-1 ..... _.,_ x x,; ,s-" o a,7: -2.0 -2.5 -- [] 0

o o

-3.0 -- -3.5 --

I I I I I I I

-4.0 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010280 (moment reference at 0.25 c).

(b) Cmq

Figure 32. Flight-determined pitch-rate derivatives (all three configurations).

SR-71 .015 configuration Baseline LASRE Test bed -- = ,_1. - ° .010 o "'% 0 .OO5 WxxW- "

I I I I I I I

-.010 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010332

(a) CN_

e 0 x 10-3 --1 -- O --2 --

x

-3-- CmSe' deg -1 -4 -- x x SR-71 -5 -- _ _ _ _'" configuration x ..... _ Baseline O - - _ -" LASRE -5 -- 0 Test bed O 0 -- _1---

I I I I I I I

-7 .5 1 .O 1.5 2.0 2.5 3.0 3.5 Mach number 010281 (b) Cm5 (moment reference at 0.25 c).

e Figure 33. Flight-determined elevator derivatives (all three configurations).

Measured value Measured value Measured value pEst estimate .... pEst estimate .... pEst estimate 3 2 T_ e, l_e, 11 e, 0 0 in. in. in.

-2

I I I

-3 -3 --4 t 5 5

51 .. /

& (_, deg

deg deg

I

-5 -5 -5 _t 0 0 deg deg deg -3 -3 -3 20 20 P_ P, 0 0 deg/sec deg/sec -20 -20 -20 4 4 rl r, r_ 0 0 degieeo deg/sec deg/sec -2 -2 -2 --4 -4 -4 10 10 _)_ -10 -10 -10 deg deg deg

-20 ] -20 -20

-30 -30 -30 .2 .2 .2 ay, ay, ay, g g g

I I I

--.2 --°2 too2 0 2 4 6 8 10 12 0 2 4 6 8 0 2 4 6 8 10 12 Time, sec Time, sec 010333 Time, eec 010334 010282

(c) Mach 2.97 and an

(a) Mach 0.61 and an (b) Mach 1.07 and an

altitude of 6,100 ft. altitude of 27,900 ft. altitude of 79,800 ft.

Figure 34. Baseline configuration lateral-directional maneuver time histories, .015 -- Test points ----e---- High-KEAS .010 -- ---.x--- Low-KEAS .005 -- 0 -- CY b -.005 -- -.010 --

I I

I I I I I

-.015 0 3.0 3.5 .5 1.0 1.5 2.0 2.5 Mach number 010335 (a) Crb.

1.5 x 10-3 Test points ---e--- High-KEAS 1.0 -- ----_-- Low-KEAS ,5 Cl b 0 -- "-.5 --

I I I I I I I

-1,0 0 .5 1.0 1.5 2.0 2,5 3.0 3.5 Mach number 010336 (b) Clb.

4 x 10 -3 Test points ------0-- High-KEAS 3 -- ----_--- Low-KEAS 2 -- Cnb 0--

I I I I J

I

-1 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number olo2s3 (C) C1%.

Figure 35. Flight-determined lateral-directional coefficient biases (baseline

configuration).

6!

n Test points High-KEAS -.002 -- "'--- Low-KEAS High-KEAS (predicted) -.004 Low-KEAS (predicted) CYI_' -.006 deg-1 -.008 -- --.010 --

I I I I I I

I

-.012 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010337

(a) C

YI3"

.5 x 10-3 0 -- .5

,,_ Jr''P

m -1.0 ... - _. =.,. ___,-.- --.. .....

T ,---O-,- Hi KEAS D g-h -1.5 •...m-. Low-KEAS High-KEAS (predicted) -2.0 Low-KEAS (predicted)

I

I I I I I I

-2.5 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010338 (b) Clf _.

2.5 x 10-3 Test points High-KEAS 2.0-- --'_ - Low-KEAS High-KEAS (predicted) 1.5-- 1.0 -- ow-KEAS (predicted) • 5 --

I

I I I I I I

.5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010284 (c) Cnf _ (moment reference at 0.25 c).

Figure 36. Predicted and flight-determined angle-of-sideslip derivatives (baseline

configuration).

1.5 -- Test points --,,e,-- High-KEAS •--m-, Low-KEAS High-KEAS (predicted) 1°0 -- .... Low-KEAS (predicted) Cyp, • 5 -- red -1 0 --

I I I I I I I

-.5 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010339 (a) Cyp.

.10 Test points Test points High-KEAS -- High-KEAS (predicted) .05 (predicted) -- T -- _- Low-KEAS .... Low-KEAS -.05 I I .....

Clp' -.10 rad-1 -.15 -.20 -.25

I

] I I I I I

--.30 3.5 .5 1.0 1.5 2.0 2.5 3.0 Mach number 010340 (b) Clp.

.20 -- Test points -,-0--- High-KEAS .15 -- •-,m,., Low-KEAS -- High-KEAS (predicted) .10 -- .... Low-KEAS (predicted) •05 -- Cnp, 0-- rad-1 -°05 -- -.10 -- -.15 --

I 1 I I I I

-.20 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 01o285 (c) Cnp (moment reference at 0.25 c).

Figure 37. Predicted and flight-determined roll-rate derivatives (baseline

configuration).

Cy r, -1 0 _ ....

rad-1 -2 -'3 _ High-KEAS i "_",D,.

-4 --..m.. Low.KE_A S -" _,,_1.

High-KEAS (predicted) _',,,, --5 Low-KEAS (predicted) _",,,,_=

-6 I I I I I _fl I

0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010341

(a) C

Yr" B 1.2 Test points High-KEAS 1.0i m t T ---'*--- LOw-KEAS High-KEAS (predicted) _ High-KEAE .8 .L .... Low-KEAS (predicted) .6 CI r, rad-1 D .4 .2 m

I I I I I I I

--.2 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010342 (b) Cir.

.5 Cn r, -.5 rad-1 Test points -,,-0.--- Hlgh-KEAS ,,, ,K-,. Low-KEAS -1.0 High-KEAS (predicted) .... Low-KEAS (predicted)

I I I I I I I

-1.5 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010286 (c) Cnr (moment reference at 0.25 c).

Figure 38. Predicted and flight-determined yaw-rate derivatives (baseline

configuration).

5 10-3 Test points -,-e-- High-KEAS •- ,_,,, = Low-KEAS High-KEAS (predicted) .... w-KEAS (predicted) CYSr' deg -1

I I I I i I I

.5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010343

(a) C

Y6 " r 4x10-4 3 -- Cl_r' _'_" -_ Test points deg-1 -1 High-KEAS •- _x_ = Low-KEAS -2 High-KEAS (predicted) -3 .... Low-KEAS (predicted)

I I I i i l I

-4 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number olo344

(b) Cta

r Ox 10 -3 i -.5 CnSr' -1.0 / ___LE''f-w- '" - " deg -1 _L]_ --r Test points -1.5 ,.=m--= LOW KEAS High-KEAS (predicted) Low-KEAS (predicted)

I I I i I I I

-2.0 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number Ol0287 (c) C% (moment reference at 0.25 c).

r

Figure 39. Predicted and flight-determined rudder derivatives (baseline

configuration).

3 x 10 -3 Test points High-KEAS 2-- •- _ - Low-KEAS High-KEAS (predicted) Low-KEAS (predicted) CY(_a' 0- deg -1 mmmm_

I I I I I I I

-3 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010345 (a) C Y8 " a 2.5 x 10 -3 Test points High-KEAS -, ,,_-- Low-KEAS 2.0-- High-KEAS (predicted) Low-KEAS (predicted) Cl(_a , 1.5 -- deg-1 1.0 • 5 --

I

1 I I I I I

0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010346

(b) Ct_

t_ 6 x 10-4 Test points _ High-KEAS High-KEAS (predicted) - Low-KEAS (predicted) Cn(_a, 2 ...... ......

deg-1 0 -- m m m m m m ml_ --2

I I I I I I I

-4 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010288 (c) C (moment reference at 0.25 c).

n 8 a

Figure 40. Predicted and flight-determined aileron derivatives (baseline

configuration).

Measured value Measured value Measured value .... pEst estimate .... pEst estimate .... pEst estimate T_ a, 11 a, in.

in.

In.

rla' 0 -2 -2 -2

31 F_r /

/--Sr /ItT-- 8a deg deg d_'_ 0 -3 -3 -3 deg deg d_

13, o_

-2 -2 P, P' 0 deg/sec deg/sec deg/sec -10

-lO I I I I -10

rl deg/sec deg/sec 0 degJsec 0

r, _

-1 -2

-2 I I I

(_)._ deg deg deg

_' o_

-5 -5 .O6 .O6 ay, ay, ay, g g g M°O6 MoO6 -.06 0 2 4 6 6 0 2 4 6 8 10 12 0 2 4 6 6 Time, see 010289 Time, sec Time, sec 010348 010347

(c) Mach 1.74 and an

(b) Mach 0.95 and an

(a) Mach 0.60 and an

altitude of 44,000 ft.

altitude of 27,500 ft.

altitude of 15,100 ft

Figure 41. LASRE configuration lateral-directional maneuver time histories.

.015 Test points High-KEAS .010 ---x--- Low-KEAS .005 CY b -.005 -.010

I I I I I I I I

-.015 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 Mach number 010349 (a) CYb.

1.5 x 10-3 Test points High-KEAS 1.0 -- ----K-- - Low-KEAS • 5 -- Cl b 0 -- mS --

I

1 t ] I ] I I I I

-1.0 0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 Mach number 010350 (b) Clb.

4 x 10-3 Test points _ High-KEAS 3-- -- --- Low-KEAS Cnb

I 1 I I I I I I I J

-1 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 Mach number 010290

(c) C

n b

Figure 42. Flight-determined lateral-directional coefficient biases (LASRE

configuration).

0 D Test points High-KEAS -.002 --.-K--- Low-KEAS -.004 -- Cy_, -.006 -- deg -1 -.008 -- -.010 --

I

I I I I I I I I

I

-.012 2.0 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 0 .2 Mach number 010351

(a) C

Y_" .5 x 10-3 -- Test points HIgh-KEAS 0-- -- ...K--- Low-KEAS -.5 -- -1.0 -- -1.5 -- -2.0 --

I

-2.s I t I I I I I I I

2.0 0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 Mach number 010352

(b) C

113"

2.5 x 10-3 Test points High-KEAS ----K--- Low-KEAS 2.0-- 1.5-- 1.0-- I • 5 --

I

I I I I I I I I I

2.0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 Mach number 010291 (c) Cnf _ (moment reference at 0.25 c).

Figure 43. Flight-determined angle-of-sideslip derivatives (LASRE configuration).

1.5 Test points High-KEAS ----K-- - Low-KEAS 1.0 Cyp, .5 rad-1 0 --

I

I I I I I I I I I

-.5 0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 Mach number 010353 (a) Cy p.

B .10 Test points .05 High-KEAS -- --_-- - Low-KEAS m -.05 Clp, D -.10 rod-1 -.15 m m.20 --,25

I I I I I I I I I t

-.30 0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 Mach number 010354 (b) C/p.

.20 -- Test points .15 -- High-KEAS --..x--- Low-KEAS .10 -- • 05 -- Cnp, 0 -- rad-1 _,05 -- -.10 -- -.15

-.20 J I I I I I I I I I

0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 Mach number 010292 (moment reference at 0.25 c).

(c) C.p

Figure 44. Flight-determined roll-rate derivatives (LASRE configuration).

Test points ----e--- High-KEAS m -I CY r, -2 -- Low-KEAS rsd-1 -3-- -4_ -5--

I I I I I I I I 1 I

-6 0 .2 ,4 .6 ,8 1.0 1.2 1.4 1.6 1.8 2.0 Mach number o_o35s

(a) C

Yr" 1.2 -- Test points ----e--- High-KEAS 1.0 D ---_--- Low-KEAS CI r, rad-1 °4 -- • 2 -- 0--

I I I I I I I I I I

--,2 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 Msch number 010356 (b) Cl .

.5-- Cn r, -.5 rad-1 Test points -1.0 ---e--- HIgh-KEAS ---N--- Low-KEAS

I I I

-1.s I I I I I I I

0 .2 .4 .6 .8 1.0 1.2 1.4 1,6 1.8 2.0 Mach number 010293 (moment reference at 0.25 c).

(c) C,,_

Figure 45. Flight-determined yaw-rate derivatives (LASRE configuration).

5 x 10 -3 Test points High-KEAS 4-- -- .-K--- Low-KEAS CY_ir, 3 -- deg-1 2 -- 1 --

I

I t I I I I I I I

0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 Mach number 010357

(a) C

Y8 r 4 x 10-4 3 -- 2 -- 1 -- Cl_r' 0 -- deg-1 Test points --2 -- High-KEAS ----K--- Low-KEAS --'3 --

-4 I I I I t I I I

0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 Mach number 010358 (b) Ct8 .

r Ox 10 -3 Test points Hlgh-KEAS ----_--- Low-KEAS -.5 -- CnSr' -I.0 -- deg-1 -1.5 --

-2.0 I I I I I I I I I I

0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 Mach number 010294 (moment reference at 0.25 c).

Figure 46. Flight-determined rudder derivatives (LASRE configuration).

3x 10-3 Test points High-KEAS 2-- ---.x-- - Low-KEAS 1 -- CYSa' 0-- deg -1 -1 -- -2 --

I I I I I I I I I I

-3 0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 Mach number 010359

(a) C

Y6 a 2.5x10 _ Test points 2.0-- --e-- Hlgh-KEAS i -- --x--- Low-KEAS 1.5-- C1(5 a' deg -1 1.0-- •5 --

I I I I I I I I I I

.2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 Mach number 010360

(b) Cl8

G 6x10 _ J[ Test points -(_ ---e--- High-KEAS 4-- --'J" -- --- ow-KEAS 2-- Cn(_a' deg -1 0-- -2 --

I I I I I I I I I I

0 • 2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 Mach number olo295 (c) C (moment reference at 0.25 c).

n8 Q Figure 47. Flight..determined aileron derivatives (LASRE configuration).

Measured value --- Measured value Measured value .... pEst estimate .... pEst estimate .... pest estimate 11 a , 118' 11 a, in.

in. in.

-2 _y 0 0

deg °l:--_J, ;. "--

deg deg

I _a_/ Ij "

-3 -3 _' 0 deg deg -1 -I -1 15 10 10 p, 5 P' 0 P' 0 deg/sec 0 deg/sec deoJsec -5 -10 -10 -10 2 2 r, r, r, deg/sec 0 deg/sec 0 deg/sec 0 -2 -2 -2 -10 0 -15 _)_ ¢, deg deg deg -20 -25 -10 .O2 .02 ay, ay, ay, g g -.02 g -.04 -.04 _., , , , , 4W --.06 --.08 0 2 4 6 8 10 12 0 2 4 6 8 0 2 4 6 8 Time, sec Time, sec 010361 Time, sec 010296 010362

(b) Mach 1.06 and an

(a) Mach 0.71 and an (c) Mach 3.01 and an

altitude of 25,500 ft. altitude of 20,000 ft. altitude of 68,300 ft.

Figure 48. Test bed configuration lateral-directional maneuver time histories.

-.015 -- Test points High-KEAS -.010 -- -.005 -- 0 -- CY b -.005 -- -.010 -- -.015 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010363 (a) Cy b .

1.5 x 10 -3 Test points 1.0 -- ,5 CI b 0 -- --.S --

I I I I I I I

-1.0 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010364 (b) Clb.

4 x 10-3 Test points ---e---- Hlgh-KEAS 3 -- -- -.x--- Low-KEAS 2. -- Cn b _--=---_-_a:_-_ ...... =------- _ o _-.

0--

,r _-_'_" _ - - _ "_----_'-_''-

I I I I I I I

-1 0 .5 1.0 1.5 2.0 2.5 3,0 3.5 Mach number 010297

(c) C

n b

Figure 49. Flight-determined lateral-directional coefficient biases (test bed

configuration).

Test points High-KEAS -.002 i --"_--- Low-KEAS -.004 Cy_, -.006 d¢_ -1 -.11118 -- -.010

I I 1 I I I I

-.012 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010365 (a) Cr .

.5x10 _ 0-- _.5 -- -1.0 -- J / Test points -1.5 -- t._ "r Hlgh-KEAS -2.0 -- --.-u--- Low-KEAS

t

1 I I I

-2.5 i 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010366 (b) fir _ .

2.5x10 _ Test points High-KEAS t "1- 2.0 --.-K--- Low-KEAS 1.5 Cn_, deg -1 1.0 .5 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010298 (c) C,q_ (moment reference at 0.25 c).

Figure 50. Flight-determined angle-of-sideslip derivatives (test bed configuration).

/ _ High-KEAS T 1.5 F Test points t rad-1

I

-.5 I 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010367

(a) C

Yp" .10 Test points __ .05 Hig_-KEAS JV_ I -.05_ -- -l¢--- Low-KEAS "_- _/_ Clp' -.10 rad-1 -.15

-.2ol

-.25

I I I I I i I

-.30 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010_8 (b) Cip.

.20 Test points B 1 5 ----e--- High-KEAS -- ..x--- Low-KEAS I- .10

3_

Cnp' t

i

rad -1 -.15 -.10

I

_.2o I I I I I 7-

0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number O10299 (c) Cnp (moment reference at 0.25 c).

Figure 5 I. Flight-determined roll-rate derivatives (test bed configuration).

3 m 2-- High-KEAS t Test points -1 CY r, I --_--- Low-KEAS rad -I -3-- -4-- -5--

I I I I I I

-6 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number OLO369 (a) CYr.

1.0 .8 Cir, .6 tad -1 .4 .2

I

.5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010370 (b) Cir.

.5 m Cn r, rad-1

I

3.5 010300 (C) Cnr (moment reference at 0.25 c).

Figure 52. Flight-determined yaw-rate derivatives (test bed configuration).

5 x 10-3 Test points T_- _ High-KEAS -- --x--- Low-KEAS 4-- 3-- CY(_r' deg -1 2--

1 I I I I I I

.5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010371

(a) C

yfi " r x 10-4 _]_ __ _.

CISr' 0 deg-1 -1 Test points -2 -- _ High-KEAS __ -- --_-- - Low-KEAS -3

I I I I I I I

-4 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010372 (b) Ct_ r .

0xl0 _ .5 Cn_r' -1.0 deg -1 1- Test -1.5 ----e--- High-KEAS

I I I I I I I

-2.0 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number oto3ol (c) C%,. (moment reference at 0.25 c).

Figure 53. Flight-determined rudder derivatives (test bed configuration).

3 x 10 -3 Test points t _ High-KEAS :]--I-.-_'- -- --x-- - Low-KEAS CYSa' deg -1 -1 -2

o F

I t I I I I I

--3 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010373

(a) C

y_ • Test points High-KEAS -- -.x--- Low-KEAS

I I I I I I

1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010374

(b) C

a 6 10-4 Test points High-KEAS 4-- "_" -- ..K-- - Low-KEAS 2-- Cn(_a' deg -1 0--

t

-2 --

I I I J I I I

-4 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010302 (c) Cn_ (moment reference at 0.25 c).

a Figure 54. Flight-determined aileron derivatives (test bed configuration).

90,000 Airspeed, KEAS 80,000 f 400 ,450 70,000 4" 60,000 f 50,000 Pressure altitude, I1 40,000 30,000

/

SR-71 20,000 / configuration O Baseline X LASRE 10,000 O Test bed

I I I

I I

2.5 3.0 3.5 0 .5 1.0 1.5 2.0 Mach number 010303

Figure 55. Flight conditions for comparison of configuration lateral-directional

stability and control derivatives.

.015 -- SR-71 configuration m .010 Baseline LASRE .005 Test bed

.... a B a -'_'"

0_-,,¢ ........ x o " oo,., o

CY b -.005 n m -,010

I I

I I I I I

-.015 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010375

(a) C

Yb" 1.5 x 10-3 SR-71 configuration X Baseline o ---e-- 1.0 -- LASRE Test bed

-.',,_ n --_--

.5 D CI b 0-- ,,,, - _ - l!

w X B5

I I I I I I I

-1.0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010376 (b) Clb.

4 x 10-3 SR-71 X configuration Baseline ----O--- 3-- LASRE Test bed Ii --O--- ! I 2-- I I I I Cnb 1 --

• _..._. ,,_, -

' _,_, r - 0 -- 7,; R1

k

I ol I I I I

-1 3.5 0 .5 1.0 1.5 2.0 2.5 3.0 Mach number 010304

Figure 56. Flight-determined lateral-directional coefficient biases (all three

configurations).

-- SR-71 configuration -.002 -- _ Baseline --_-- LASRE -.004 -- o- o - Test bed O

Cy _, -.006 - x _

deg -1 -.008 -- -.010 --

I I I I I I I

-.012 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number Ol0377 (a) Cyf.

.5 10-3 0-- x -- B... -J_" T'Zll H_ oJ_ x_ ._o_ -.5 CIl_' -1.0 -- deg-1 _,l_/6u I., SR-?I O/a" configuration -1.5 ,O _ Baseline - -_ - - LASRE -2.0 -- _ - - Test bed

"l I I I I I I

-2.5 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010378 (b) Ctf _ .

2.5 x 10-3 SR-71 configuration 2.0 - -_ - - LASRE Test bed _'t_. _ Baseline

Cn w 1.s --

deg-1 1.0 --

a b

.5

"_.._= o °=

I I 1 1 °r"=e-"'l I

.5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010305 (c) C% (moment reference at 0.25 c).

Figure 57. Flight-determined angle-of-sideslip derivatives (all three

configurations).

1.5 SR-71 configuration Baseline 1.0 ---)_m. LASRE o Test bed Cyp, .5 rad-1 0 -- O O x a o Oo o o

I ] " I I I I

--,5 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010379

(a) C

Yp" .10 B .05

D /

B 0 a °,/ 0 -.05 O 00 s" Clp' -.10 O 0 rad-1 0 r_ . n ° __0 m 0 SR-71 B -.15 ---; ___b'_'_0o go 0 - configuration . _x x x ---e--- Baseline --.20 - N - " O --m-- LASRE -.25 - - _3- - - Test bed

I I I I I l I

-.30 .5 1.0 1.5 2.0 2.5 3.0 3.5 Msch number 010380 (b) Clp.

.20 -- SR-71 configuration .15 ---0--- Baseline .10 -- - -_ -- LASRE 0 .05 - - Q-- - Test bed o o O 0 Cnp ' 0

a _.1_=._-. --- - ?_ °o ' "

rad-1 --.05 D x _ -.10 Q D -.15

I I I I I I I

--.20 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010306 (C) Cnp (moment reference at 0.25 c).

Figure 58. Flight-determined roll-rate derivatives (all three configurations).

3 -- SR-71 configuration 2-- Baseline 1 -- --_-- LASRE 0-- _ _, ^ O L _ O a - - "13- - - Test bed CY r, --2 red -1 --3 -- [] -_'_X [] [] O -4 -- [] -5 --

I

I I I I I

-6 3.5 .5 1.0 1.5 2.0 2.5 3.0 Mach number 010381

(a) C

Yr" 1.2 SR-71 configuration 1.0 Baseline p--e_ LASRE [] Test bed o O .6 0 oo CI r, O [] 0 0 .0 -'B''''''D'" rad-1 p .4 • 2 --

v= x -_'--#-_i3 TM - a " 0

o 0P

I

I I I I I [

--,2 3.5 .5 1.0 1.5 2.0 2.5 3.0 Mach number 010382 (b) CI .

.5

_P

Cn r, -.5 SR-71 o [] rad-1 configuration Baseline -1.0 LASRE a B _ raN,_.

Test bed

I I 1 I I I

-1.5 3.5 0 .5 1.0 1.5 2.0 2.5 3.0 Mach number 010307 (c) Cnr (moment reference at 0.25 c).

Figure 59. Flight-determined yaw-rate derivatives (all three configurations).

5 x 10-3 SR-71 configuration 4-- O - -)* - - LASRE [] - - "B- - - Test bed 3-- X I _, .,,x. _ Baseline tqL.X 2-- 1 --

I I I I I I I

.5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010383

(a) C

Y6" r 4x10 _ 3 -- 2 --

° *"

1 --

- , 4J a

ClSr' O -- m n \ _ _ SR-71 deg-1 -1 -- x _ a_ configuration --e--- Baseline I !

-2 -- I v --_-- LASRE - - "!3" - - Test bed X

I I I I I I I

.5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010384 0xlO _ -- SR-71 configuration ---e--- Baseline -.5 -- --_-- LASRE __"_-I3 --_-- Test bed j_r,_'_ "o -I Cn(_r' -1.0 deg -1 -1.5

-2.0 I I I I I I

0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010308 (c) Cn8 (moment reference at 0.25 c).

t" Figure 60. Flight-determined rudder derivatives (all three configurations).

3 x 10-3 SR-71 configuration 2-- --e--- Baseline LASRE Test bed 1 -- O m ° '_1" °" CY(_a' 0-- deg -1 --2 --

I I [ 1 I I I

-3 0 .5 1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010385

(a) C

Y8

SR-71 configuration Baseline LASRE Test bed Cl(_a' deg-1

I J I I I I

1.0 1.5 2.0 2.5 3.0 3.5 Mach number 010386

(b) Ct_

a SR-71 configuration x x -'x"_" x- _ "-'0-- Baseline 4 -- X k LASRE X l ----14'--- Test bed

" _ JL -_"

2-- CnSa' .... u.. _--_ -; - o i_'__'.. ,,,.

deg -1 n 0 ",I_.

0-- llg rl -2--

"I I " [ I I I

-4 0 .5 1.0 1.5 2.0 2.5 3.0 Mach number (c) C% (moment reference at 0.25 c).

a Figure 61. Flight-determined aileron derivatives (all three configurations).

REPORT DOCUMENTATION PAGE Fo,m Approved

OMB No. 0704-0188 PuL_)C reporting burden for this co;leclion of information is eslimated to average t hour per response, including the time for reviewing instructions, searching existing data sources, gathenng and maintaining the dala needed, and completing and reviewing the collection of information Send comments regarding this burden estimate or any other aspect of this collection ol information, including suggestions for reducing this burden, to Washington Headquarters Services, Oirectorale lor information Operations and Reports. 1215 Jefferson Daws Highway Suite 1204. Arlington.

VA 22202-4302, and to tt_e Office of Management and Bu0get, Paperwork Re0uc_ion Pro)ect (0704-0188), Washint_on, DC 20503 1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE

June 2002 Technical Publication

3. REPORTTYPE AND DATES COVERED 4.TITLE AND SUBTITLE 5. FUNDING NUMBERS

for the SR-71

Stability and Control Estimation Flight Test Results

Aircraft With Externally Mounted Experiments

WU 710-35-14-E8-OM-00-844

6. AUTHOR(S)

Timothy R. Moes and Kenneth Iliff

8. PERFORMING ORGANIZATION 7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) REPORT NUMBER

NASA Dryden Flight Research Center

EO. Box 273

H-2465

Edwards, California 93523-0273

10. SPONSORING/MONITORING 9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) AGENCY REPORT NUMBER

National Aeronautics and Space Administration

NASA/TP-2002-210718

Washington, DC 20546-0001

11.SUPPLEMENTARY NOTES 12a. DISTRIBUTION/AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE

Unclassified--Unlimited

Subject Category 08

This report is available at http://www.dfrc.nasa.gov/DTRS/

13. ABSTRACT (Maximum ZOO words) A maximum-likelihood output-error parameter estimation technique is used to obtain stability and

control derivatives for the NASA Dryden Flight Research Center SR-71A airplane and for

configurations that include experiments externally mounted to the top of the fuselage. This research is being done as part of the envelope clearance for the new experiment configurations. Flight data are obtained at speeds ranging from Mach 0.4 to Mach 3.0, with an extensive amount of test points at approximately Mach 1.0. Pilot-input pitch and yaw-roll doublets are used to obtain the data. This report defines the parameter estimation technique used, presents stability and control derivative results, and compares the derivatives for the three configurations tested. The experimental configurations studied generally show acceptable stability, control, trim, and handling qualities throughout the Mach regimes tested. The reduction of directional stability for the experimental configurations is the most significant aerodynamic effect measured and identified as a design constraint for future experimental configurations. This report also shows the significant effects of aircraft flexibility on the stability and control derivatives.

14. SUBJECT TERMS 15. NUMBER OF PAGES Control derivatives, Maximum likelihood estimates, Parameter identification, 16. PRICE CODE

SR-71, Stability derivatives

A03

18. SECURITY CLASSIRCATION 19, SECURITY CLASSIFICATION 20. LIMITATION OF ABSTRACT 17. SECURITY CLASSIFICATION OF THIS PAGE OF ABSTRACT OF REPORT

Unclassified Unclassified Unlimited Unclassified

NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) Preecobed by ANSI Std Z39_1B 298-102

Source & rights

Source: ntrs.nasa.gov. Public-domain U.S. Government work (17 USC §105) — freely reproducible.

Permanent URL — we don’t break links.

Report a problem or request removal

Document details

Doc number
NASA/TP-2002-210718
Publisher
NASA (NTRS)
Year
2002
Pages
96
File size
3.4 MB