ANALYSIS OF LARGE AMPLITUDE MANEUVER DATA

The usual procedure for aircraft parameter estimation is to perturb it slightly from its trim position by deflecting one or more of its control surfaces and gathering flight data from such small amplitude maneuvers. However, in practice it may not be possible to trim the aircraft at certain test points in the flight envelope. Under such conditions, large amplitude maneuvers (LAMs) are more useful in determining the aerodynamic derivatives. The approach is not valid in the region where stability and control derivatives vary rapidly. The results of the LAM are generated for the fighter aircraft FA2 for the conditions: (1) simulated longitudinal LAM data at two flight conditions and analyzed using the SMLR (stepwise multiple regression method [1]) without partitioning—here the entire data set is handled as a single LAM; (2) the data as in (1) but analyzed using the model error method (MEM [1]); and (3) the real flight data (*) analyzed using the previous two approaches [28]. The results are shown in Table 9.13.

The other method consists of partitioning the LAM that covers a large AOA range into several bins, each of which spans a smaller range of AOA [29]. These data would not be contiguous and therefore the regression method is employed to

Подпись: System Identification and Parameter Estimation 295

TABLE 9.12

Fit Errors for Lift and Drag Coefficients—Parameter Estimation

RC Maneuver

Set 1 Set2 Set3 Set 4

Estimation

Mach No. =

0.6

Mach No. =

0.7

Mach No. =

0.8

Mach No. =

1.0

Method

Q

CD

Q

Co

Q

Co

Q

Co

EUDF-NMBA

0.0216

0.1586

0.0201

0.1112

0.0190

0.1036

0.0442

0.4275

EUDF-MBA

0.2597

0.9108

0.3940

0.6660

0.4939

1.0506

0.5439

0.3659

SOEM

0.2198

0.7540

0.3747

0.5979

0.4746

1.0009

0.4973

0.2996

WUT Maneuver

Set 1

Set 2

Set 3

Set 4

Mach No. =

0.6

Mach No. =

0.7

Mach No. =

0.8

Mach No. =

1.0

EUDF-NMBA

0.0411

0.5960

0.0543

0.4930

0.0624

0.4588

0.0423

0.2261

EUDF-MBA

0.6193

0.9564

0.4524

1.1780

0.3922

0.8220

0.3995

0.7903

SOEM

0.3871

0.8862

0.4811

1.2402

0.5078

0.8579

0.2487

0.4741

Note: Altitude equal to 8 km.

 

TABLE 9.13

Longitudinal LAM Data Results

Trim

AOA 0

Ref.

SMLR

MEM

Ref.

SMLR

MEM

Ref.

SMLR

MEM

11.90

0.077

0.071

0.081

-1.173

-1.295

-1.252

-0.391

-0.375

-0.379

4.3

-0.014

-0.01

– 0.020

-1.32

-1.19

-1.16

-0.415

-0.398

– 0.411

11.55*

0.077

0.058

0.061

-1.1733

-0.992

-0.922

-0.391

-0.394

-0.399

Подпись: System Identification and Parameter Estimation 295

TABLE 9.12

Fit Errors for Lift and Drag Coefficients—Parameter Estimation

RC Maneuver

Set 1 Set2 Set3 Set 4

Estimation

Mach No. =

0.6

Mach No. =

0.7

Mach No. =

0.8

Mach No. =

1.0

Method

Q

CD

Q

Co

Q

Co

Q

Co

EUDF-NMBA

0.0216

0.1586

0.0201

0.1112

0.0190

0.1036

0.0442

0.4275

EUDF-MBA

0.2597

0.9108

0.3940

0.6660

0.4939

1.0506

0.5439

0.3659

SOEM

0.2198

0.7540

0.3747

0.5979

0.4746

1.0009

0.4973

0.2996

WUT Maneuver

Set 1

Set 2

Set 3

Set 4

Mach No. =

0.6

Mach No. =

0.7

Mach No. =

0.8

Mach No. =

1.0

EUDF-NMBA

0.0411

0.5960

0.0543

0.4930

0.0624

0.4588

0.0423

0.2261

EUDF-MBA

0.6193

0.9564

0.4524

1.1780

0.3922

0.8220

0.3995

0.7903

SOEM

0.3871

0.8862

0.4811

1.2402

0.5078

0.8579

0.2487

0.4741

Note: Altitude equal to 8 km.

 

estimate derivatives. This approach was also used for the analysis of the data mentioned earlier, and the detailed results can be found in Ref. [28].

For LD data generation from LAMs (in an RTS), the modes were excited in roll and yaw axes by giving aileron and rudder inputs superimposed on a steadily increasing AOA due to slow but steady increase in elevator deflection, thereby generating LAM data. Normal practice is to estimate linear derivative models but, if necessary, the SMLR approach can be used to determine the model structure with higher-order terms for better representation of aircraft dynamics. A pilot flew an RTS (Chapter 6) and persistently excited the roll and yaw motions at increasing AOA in pullup. The data were generated for Mach numbers 0.4 and 0.6 with AOA variation between 2° and 20° [30]. Significant variation in yaw during lateral stick inputs, possibly caused by aileron-rudder interconnect, was seen in the time – history plots. Regression and data partitioning techniques were used for estimation of the LD derivatives. For each Mach, the data were partitioned into bins, where each bin corresponds to a mean value of AOA: at Mach 0.4, the total variation in AOA is from 2° to 20°. The data from 2° to 4° are put into one group; the data from 4° to 6° are put into another group, and so on until all the data are accounted for. The data in a group would have resulted from different portions of a single maneuver or from certain portions of several other maneuvers. These data in each group were analyzed using regression. Each bin should have sufficient number of data points for successful estimation of derivatives and should have sufficient variations in the aircraft motion variables so that accurate estimation of derivatives is possible. Several criteria, e. g., the number of data points in each bin, the CRBs of the estimated derivatives, and the percentage fit error, were evaluated to check the quality of the estimates. The derivatives C„b, C/, C„v and C/ were estimated and cross-plotted with the reference derivatives. The results were very encouraging. One compound derivative (for the flight condition Mach number = 0.6 and with auto slat) C = C^ + Q tan a is shown Figure 9.6 [30]. This RTS exercise established that the approach if followed would work for the data gathered from real flight tests. Using this approach it is possible to determine the aerodynamic derivatives at flight conditions that cannot be covered during 1g trim flight.