NASA/TP--2000-209908 DOT/FAA/AR-99/85

U.S. Department of Transportation Federal Aviation Administration

NASA/FAA Icing Program: Flight Test Report

Thomas P. Ratvasky Glenn Research Center, Cleveland, Ohio

Judith Foss Van Zante Dynacs Engineering Company, Inc., Brook Park, Ohio

Alex Sim Dryden Flight Research Center Edwards Air Force Base, California

March 2000 The NASA STI Program Office... in Profile

Since its founding, NASA has been dedicated to CONFERENCE PUBLICATION. Collected the advancement of aeronautics and space papers from scientific and technical science. The NASA Scientific and Technical conferences, symposia, seminars, or other Information (STI) Program Office plays a key part meetings sponsored or cosponsored by in helping NASA maintain this important role. NASA.

The NASA STI Program Office is operated by SPECIAL PUBLICATION. Scientific, Langley Research Center, the Lead Center for technical, or historical information from NASNs scientific and technical information. The NASA programs, projects, and missions, NASA STI Program Office provides access to the often concerned with subjects having NASA STI Database, the largest collection of substantial public interest. aeronautical and space science STI in the world. The Program Office is also NASA's institutional TECHNICAL TRANSLATION. English- mechanism for disseminating the results of its language translations of foreign scientific research and development activities. These results and technical material pertinent to NASA's are published by NASA in the NASA STI Report mission. Series, which includes the following report types: Specialized services that complement the STI TECHNICAL PUBLICATION. Reports of Program Office's diverse offerings include completed research or a major significant creating custom thesauri, building customized phase of research that present the results of data bases, organizing and publishing research NASA programs and include extensive data results.., even providing videos. or theoretical analysis. Includes compilations of significant scientific and technical data and For more information about the NASA STI information deemed to be of continuing Program Office, see the following: reference value. NASA's counterpart of peer- reviewed formal professional papers but • Access the NASA STI Program Home Page has less stringent limitations on manuscript at http:/lwww.sti.nasa.gov length and extent of graphic presentations. • E-mail your question via the Internet to TECHNICAL MEMORANDUM. Scientific [email protected] and technical findings that are preliminary or of specialized interest, e.g., quick release • Fax your question to the NASA Access reports, working papers, and bibliographies Help Desk at (301) 621-0134 that contain minimal annotation. Does not contain extensive analysis. • Telephone the NASA Access Help Desk at (301) 621-0390 CONTRACTOR REPORT. Scientific and technical findings by NASA-sponsored Write to: contractors and grantees. NASA Access Help Desk NASA Center for AeroSpace Information 7121 Standard Drive Hanover, MD 21076 NASA/TP--2000-209908 DOT/FAA/AR-99/85

U.S. Department of Transportation Federal Aviation Administration

NASA/FAA Tailplane Icing Program: Flight Test Report

Thomas P. Ratvasky Glenn Research Center, Cleveland, Ohio

Judith Foss Van Zante Dynacs Engineering Company, Inc., Brook Park, Ohio

Alex Sim Dryden Flight Research Center Edwards Air Force Base, California

National Aeronautics and Space Administration

Glenn Research Center

March 2000 Acknowledgments

The NASA/FAA Tailplane Icing Program flight testing phases were accomplished only through the dedicated effort of the outstanding flight operations and technical staff at the NASA Glenn Research Center. Special recognition is made to Richard Ranaudo, who served as chief pilot and lead pilot for the TaJlplane Icing Program flight tests. His unparalleled knowledge of icing effects on aircraft performance and stability and control provided an increased safety margin in this risky flight test activity. The contributions of pilots Bill Rieke, Bob McKnight and Dennis O'Donoghue, and operations engineers Ed Emery, Chris Hegedus, Mike Ernst and Jon Yencho were key to the productivity of the flight tests. Special acknowledgment to crew chief Dan Gorman and fabrication specialist Steve Hayes for their accomplishments of preparing the aircraft for this unprecedented type of flight testing. Regarding the dynamic data analysis, the authors would like to thank Tom Nettleton of DeHavilland for suggesting the cross-plot method, and also Mike Hinson of Learjet and Dana Herring of Raytheon for their analysis suggestions. The authors would also like to recognize Jim Riley of the FAA William J. Hughes Technical Center for his continuous encouragement throughout this four-year program and careful review of this document. He provided much guidance and motivation throughout the entire program. Finally, we would like to express our appreciation to the sponsors of this research effort: the NASA Aerospace Operations Systems Program and the FAA William J. Hughes Technical Center.

Trade names or manufacturers' names are used in this report for identification only. This usage does not constitute an official endorsement, either expressed or implied, by the National Aeronautics and Space Administration.

Available from

NASA Center for Aerospace Information National Technical Information Service 7121 Standard Drive 5285 Port Royal Road Hanover, MD 21076 Springfield, VA 22100 Price Code: A08 Price Code: A08 Table of Contents

Executive Summary ...... viii

List of Symbols and Abbreviations ...... ix

1.0 Introduction ...... 1

2.0 Research Aircraft and Ice Shapes ...... 3

3.0 Instrumentation System ...... 5

3.1 Aircraft Dynamics System ...... 5 3.2 Tail Aeroperformance System ...... 5 3.3 Video Systems ...... 6 3.4 Calibration Procedures ...... 6

4.0 Data Handling Process ...... 7

5.0 Steady State Analysis ...... 9

5.1 Flight Test Procedures ...... 9 5. i. 1 Steady-Wings-Level (SWL): ...... 9 5.1.2 Steady-Heading Sideslips (SHSS): ...... 9 5.2 Method of Analysis for Steady State Data ...... 9 5.2.1 General ...... 9 5.2.2 Data Consistency Check ...... 9 5.3 Results ...... 10 5.3.1 Basic Aircraft Performance Description ...... 10 5.3.2 Tailplane Angle of Attack Description ...... 12 5.3.3 Tailplane Lift Description ...... 14 5.3.4 Hinge Moment Description ...... 17 5.3.5 Icing Effects on Tailplane Lift ...... 19 5.3.6 Icing Effects on Elevator Hinge Moment ...... 24 5.3.7 Center of Gravity Effects on Tailplane Lift ...... 27 5.3.8 Thrust Effects on Tailplane Flow Field and Performance ...... 29 5.3.9 Spanwise Variation of Flow at the Horizontal Tail ...... 31 5.3.10 Effect of Steady Heading Sideslips on Tailplane Flow Field and Performance ...... 33

6.0 Dynamic Maneuver Analyses ...... 37

6.1 Maneuver Description and Analysis ...... 38 6.1.1 Data Compatibility ...... 38 6.1.2 Flap Transition ...... 39 6.1.3 Speed Transition ...... 42 6.1.4 Thrust Transition ...... 45 Propwash Effect ...... 47 6.1.5 Pushover ...... 49 Shaped Maneuver ...... 53 Stall ...... 53 Advanced Analysis Methods ...... 54 6.1.6 Elevator Doublet ...... 58 6.1.7 Balked Landing ....,...... 63

NASA/TP--2000-209908 iii 6.1.8 A Case of Turbulence ...... 64 6.2 Comparison Between the Pushover and Elevator Doublet ...... 65 6.3 Repeatability Analysis for the Pushover ...... 68 6.3.1 Description of the Analysis ...... 68 6.3.2 Description of the Results ...... 69 6.4 Dynamic Maneuvering Conclusions ...... 72 6.4.1 Transition Maneuvers ...... 72 6.4.2 Dynamic Maneuvers ...... 72

7.0 Parameter Estimation Analysis ...... 75

7.1 Flight Test Procedures ...... 75 7.2 Method of Analysis for Parameter Estimation ...... 76 7.3 Results ...... 76 7.4 Parameter Estimation Conclusions ...... 84

8. Guest Pilot Workshop ...... 85

8.1 Attendees ...... 85 8.2 Guest Pilot Workshop Flight Test Procedures and Maneuvers ...... 85 8.3 Handling Qualities Assessment ...... 86

9. Concluding Remarks ...... 89

9.1 Data Analysis Techniques ...... 89 9.2 Maneuvers ...... 90 9.3 Aircraft Configuration ...... 90 9.4 Closing Commentary ...... 91

10. References ...... 93

Appendix A. Flight Test Summaries, Research Aircraft, Ice Shapes and Instrumentation ...... 95

Appendix B. Select Steady State Data Plots ...... 101

Appendix C. Time Histories of Select Dynamic Maneuvers ...... 105

Appendix D. Pressure Coefficient Distributions of the Full Tail Stall ...... 147

Appendix E. Recommendations to the Pilot ...... 153

NASA/TP---2000-209908 iv List of Figures

Figure 1• Total lift coefficient vs. aircraft angle of attack for various flap settings ...... 10 Figure 2. Typical lift curve during landing transition ...... 10 Figure 3. Force and moment diagram ...... 11 Figure 4. Tailplane angle of attack vs. aircraft angle of attack ...... 12 Figure 5. Tailplane angle of attack vs. airspeed ...... 13 Figure 6. Typical tail AOA during landing transition ...... 13 Figure 7. Horizontal tail diagram ...... 14 Figure 8. Tail lift coefficient from wind tunnel data ...... 14 Figure 9. Cp distribution for Baseline tail ...... 15 Figure 10. Ctrait vs. o_aitfrom flight data ...... 15 Figure 11. Ct_raitvs. VTAS from flight data ...... 15 Figure 12. Tail lift coefficient from wind tunnel data and flight data ...... 16 Figure 13. Horizontal stabilizer, elevator and ...... 17 Figure 14. Elevator hinge moment coefficient vs. tail angle of attack ...... 18 Figure 15. Elevator hinge moment coefficient vs. airspeed ...... 18 Figure 16. Tail lift coefficients for _F = 0 O ...... 19 Figure 17. Tail lift coefficients for SF = 10°...... 19 Figure 18. Tail lift coefficients for SF = 20 °...... i...... 20 Figure 19. Tail lift coefficients for SF = 30 °...... 20 Figure 20. Tail lift coefficients for SF = 40 °...... 20 Figure 21. Elevator deflection for SF = 0°...... _...... 21 Figure 22. Elevator deflection for 61: = 10 °...... 21 Figure 23. Elevator deflection for 6F = 20 °...... 22 Figure 24. Elevator deflection for 61== 30 °...... 22 Figure 25. Elevator deflection for 6F = 40 °...... 22 Figure 26. Elevator hinge moment coefficient for 8F = 0°...... 24 Figure 27. Elevator hinge moment coefficient for _F = 10°...... 24 Figure 28. Elevator hinge moment coefficient for _F = 20 °...... 25 Figure 29. Elevator hinge moment coefficient for d_F= 30"...... 25 Figure 30. Elevator hinge moment coefficient for cSF= 40 °...... 25

Figure 31. CG effect: Cl_Tail VS. O_ail ...... 27 Figure 32. CG effect: Ct_rai t vs. VTAS ...... 27 Figure 33. CG effects on C_r,,t with respect to Cz,vc ...... 28

Figure 34. Thrust effect on O_ail...... 29 Figure 35. Thrust effect on Vl,_il...... 30 Figure 36. Thrust effect on Cl_r_il...... 30 Figure 37. Span variation of c_,t ...... 31 Figure 38. Span variation of titbit...... 32 Figure 39. Span variation of Vigil...... 32 Figure 40. Sideslip effect on o_it ...... 33 Figure 41. Sideslip effect on a,vc...... 34 Figure 42. Sideslip effect on CiTail ...... 35 Figure 43. Data compatibility time traces ...... 38 Figure 44. Flap transition; clean tail, VIAS = 85kts ...... 40 Figure 45. Flap transition; Failed Boot, VIAS = 80kts (VTAS = 91kts) ...... 41 Figure 46. Speed transition time histories; Failed Boot at 6F = 20 ° and Cr = 0.1 ...... 43 Figure 47. Detail of speed transition: (a) ai time history with two highlighted data points, and (b) CLr,, vs. al ...... 44 Figure 48. Thrust transition time histories; Failed Boot at d_F= 40 ° and V= 1.5Vs ...... 46 Figure 49. Propwash effect for a clean tail at different flap settings and constant speed, V = 1.25Vs ...... 48 Figure 50. Non-critical pushover time histories: clean tail, SF = 0 ° and V= 1.5Vs = 100kts ...... 51 Figure 5 I. Critical pushover time histories: Failed Boot ice shape, _F = 20 ° and V = 1.06Vs = 55kts ...... 52

NASA/TP--2000-209908 v Figure 52. Example of elevator shaping during pushover maneuver ...... 53 Figure 53. Example of a trailing edge tail stall. Baseline, _F = 40 °, V = 0.94Vs = 50kts, Cr = 0.23 ...... 54 Figure 54. Co-plots of 8E and FYE for pushover. Failed Boot, VTAS = 80kts, increasing flap deflection ...... 55 Figure 55. Co-plots of 8E and q for pushover. Same data base as in Figure 54 ...... 55 Figure 56. Cross-plots for Baseline case at 8F = 0° and V = 1.5Vs ...... 56 Figure 57. Cross-plots for same database as Figure 54. Top row is cSF= 0% middle 10 ° and bottom 20 °...... 57 Figure 58. Cross-plots corresponding to data in Figure 53, the trailing edge stall pushover case ...... 57 Figure 59. Elevator doublet for a clean tail, SF = 0 °, V = 1.5Vs = 100kts ...... 59 Figure 60. Elevator doublet for the Failed Boot, 8F = 20 °, V = 1.1Vs = 55kts ...... 60 Figure 61. Pitch response and damping characteristics for extreme elevator doublet cases ...... 61 Figure 62. Cross-plots of q vs. d_Efor elevator doublet data from Figure 61 ...... 62 Figure 63. Cross-plots of N_ vs. FYE for elevator doublet data from Figure 61 ...... 62 Figure 64. Balked landing tracking task ...... 63 Figure 65. Comparison between quiescent and turbulent environment. Failed Boot, _SF= 0° flaps ...... 64 Figure 66. Three maneuvers (thin lines) and their average (thick line) vs. scaled time for an Inter-Cycle ice, _SF= 0°, V=I .5Vs = 100kts test point. Error bars are indicated on select traces ...... 70 Figure 67. Additional elevator required to trim ...... 77 Figure 68. Icing effects on pitching moment derivatives, SF = 0°...... 79 Figure 69. Icing effectson pitching moment derivatives, 8F = 10°...... 80 Figure 70. Icing effects on pitching moment derivatives, 8F = 20 °...... 81 Figure 71. Icing effects on pitching moment derivatives, cSF= 30 °...... 82 Figure 72. Icing effects on pitching moment derivatives, _F = 40 °...... 83 Figure 73. Cooper-Harper handling qualities rating flow chart ...... 86 Figure 74. Cooper-Harper ratings for approach and go-around MTE ...... 87

Figure A- 1. NASA Glenn Research Center icing Research Aircraft ...... 97 Figure A-2. Inter-Cycle Ice on tailplane ...... 98 Figure A-3. Failed Boot Ice on tailplane ...... 98 Figure A-4. S&C Ice on tailplane ...... 99 Figure A-5. Tail flow probes and pressure belt ...... 99 Figure A-6.5-hole tail flow probe head ...... 100 Figure A-7. Tall camera position ...... 100

Figure B-1. 1995 data consistency results ...... 101 Figure B-2. 1997 data consistency results ...... 102 Figure B-3. Tail flow field during SHSS maneuvers ...... 103

Figure C- 1. Pushover Baseline, 5F = 0°, VIAS = 1.50Vs = 100kts ...... 106 Figure C-2. Pushover Baseline, 8F = 40 °, VIAS = 1.06Vs = 55kts ...... 108 Figure C-3. Pushover Inter-Cycle Ice, _F = 0°, VIAS = 1.50V s = 100kts ...... 110 Figure C-4. Pushover Inter-Cycle Ice, fiF = 40 °, VIAS = 1.06Vs = 55kts ...... 112 Figure C-5. Pushover Failed Boot, 8F = 20 °, VIAS = 1.33Vs = 75kts ...... 114 Figure C-6. Pushover Failed Boot, 5F = 20 °, VIAS = 0.97Vs = 55kts ...... 116 Figure C-7. Pushover S&C, _F = 0% VIAS = 1.50Vs = 100kts ...... 118 Figure C-8. Pushover S&C, _F = 20 °, VIAS = 0.97Vs = 55kts ...... 120 Figure C-9. Elevator doublet, Baseline, _F = 0°, VIAS = 1.50Vs = 100kts ...... 122 Figure C-10. Elevator doublet Baseline, 15F= 40 °, VIAS = l_06Vs = 55kts ...... ' ...... 124 Figure C- 11. Elevator doublet Inter-Cycle Ice, 8F = 30 °, VIAS = 1.59Vs = 85kts ...... 126 Figure C- 12. Elevator doublet Inter-Cycle Ice, _iF = 40 °, VIAS = 1.63Vs = 85kts ...... 128 Figure C- 13. Elevator doublet Failed Boot, _SF= 30 °, VIAS = 1.59Vs = 85kts ...... 130 Figure C- 14. Elevator doublet Failed Boot, 8F = 30 °, VIAS = 1.03Vs = 55kts ...... 132 Figure C- 15. Elevator doublet S&C, 8F = 20 °, VIAS = 0.97Vs = 55kts ...... 134 Figure C- 16. Elevator doublet S&C, _SF= 30 °, VIAS = 1.40Vs = 75kts ...... 136

NASA]TP--2000-209908 vi FigureC-17.Thrusttransition,Baseline,8F=0°, VIAS = 1.20Vs = 80kts ...... 138 Figure C-18. Thrust transition, Baseline, 5F = 40 °, VIAS = 1.63Vs = 85kts ...... 140 Figure C-19. Thrust transition, Failed Boot, 8F = 0°, VIAS = 1.20Vs = 80kts ...... 142 Figure C-20. Thrust transition, Failed Boot, _SF= 40 °, VIAS = 1.63Vs = 85kts ...... 144

Figure D-1. Select time histories of the stall event. Circles indicate times presented in Figure D-3 ...... 147 Figure D-2. Location of pressure taps on clean airfoil ...... 148 Figure D-3. Plots of Cp for full tail stall event at various times. Solid, blue line for lower surface ...... 149

List of Tables

Table 1. Dynamic maneuvers flown and number of test points for each ...... 37 Table 2. Pushover - elevator doublet comparison ...... 65 Table 3. Pass/Fail maps for both the control force and pitch response ratings during a pushover maneuver ...... 67 Table 4. Repeatability analysis - precision levels for several flight conditions to minumum Nz ...... 71 Table 5. Repeatability analysis - accuracy of achieving the target values for non-critical and critical cases ...... 71 Table 6. Comparison of the advantages and disadvantages for both the pushover and elevator doublet ...... 74 Table 7. Conditions for the parameter estimation maneuvers ...... 75 Table 8. Guest pilot attendees ...... 85 Table 9. Critical and non-critical aircraft configurations for tailplane stall ...... 90

Table A-1. 1995 flight test summary ...... 95 Table A-2. 1997 flight test summary ...... 95 Table A-3. Physical characteristics of research aircraft ...... 96 Table A-4. Instrumentation specifications ...... 96

NASA/TP--2000-209908 vii ExecutiveSummary

The purpose of this report is to present the results from research flights that explored the characteristics of an ice- contaminated tailplane. Using the NASA Twin Otter Icing Research Aircraft, fifty-one test flights were conducted with various simulated ice shapes attached to the of the horizontal tailplane. A clean leading edge pro- vided the baseline case, then three ice shapes were flown in order of increasing severity. Flight tests with each ice shape included both steady state and dynamic maneuvers. The steady state points were 1G wings level and steady heading sideslips. The primary dynamic maneuvers were pushovers to approximately 0.5G, 0.25G, and 0G; elevator doublets; thrust transitions and simulated balked landings. These maneuvers were conducted for a full range of flap positions and aircraft angle of attack where possible. In some iced cases, full flaps could not be deployed without leading to a full tail stall.

The analysis of this data set has clearly demonstrated the detrimental effects of ice contamination on aircraft stability and controllability. Paths to tailplane stall were revealed through parameter isolation and transition studies. These paths are (1) increasing ice shape severity, (2) increasing flap deflection, (3) high or low speeds, depending on whether the aircraft is in a steady state (high speed) or pushover maneuver (low speed), and (4) increasing thrust. It is suspected that the adverse response to thrust is configuration specific, i.e., it occurs when the thrust line is positioned above the aircraft's center of gravity.

The two dynamic maneuvers, the pushover and elevator doublet, were studied in depth. Analysis of the pushover provided an understanding of its effect on tailplane angle of attack and the imminent tail stall characteristics. A pos- sible alternate maneuver, the elevator doublet, was compared to the pushover. From this comparison, it became ap- parent that the extreme tailplane angles of attack attained by the pushover are not a necessary precursor to tail stall. Rather, it is the combination of tail angle of attack and elevator deflection; i.e., the more extreme the elevator deflec- tion trailing edge up, the lower the stall margin. Finally, the elevator doublet was also used for a stability and control parameter estimation process. These results indicated that the elevator control effectiveness derivative, C,,_E, was severely degraded by the ice contamination.

Although it was not the intention of this research program, a full tail stall was experienced during one of the thrust transition test points. The data acquired during this event provided great insight into the recovery procedure. It is interesting to note from this one experience that even though corrective actions - yoke back, power to idle and flaps up - were applied within 0.4 seconds, the aircraft continued to lose 300 ft of altitude over 8 seconds before the air- craft recovered.

The concluding phase of the Tailplane Icing Program was a Guest Pilot workshop. It was conducted to quickly transfer "lessons learned" on ice-contaminated tailplane stall to the FAA and Transport Canada, airframe manufac- turers, and the general public through the aviation media. The Guest Pilot Workshop served as an excellent forum to transfer this information and to obtain handling quality ratings of the test aircraft from the guest pilots. The handling quality results unanimously demonstrate the degradation of the aircraft with increased flap angle. In addition, an edu- cational video for pilots entitled Tailplane Icing was produced and disseminated internationally. Of crucial impor- tance to the pilot is the point that although the difference in tactile cues between an impending wing and tail stall are subtle, the recovery procedures are diametric.

The flight research effort was very comprehensive, but limited to one aircraft, and therefore did not examine effects of tailplane design and location, or other aircraft geometry configuration effects. However, this effort provided the role of some of the parameters in promoting tailplane stall. The lessons learned from this effort will provide guidance to regulatory agencies, aircraft manufacturers, and operators on ice-contaminated tailplane stall in the effort to in- crease aviation safety and reduce the fatal accident rate.

NASMTP--2000-209908 viii List of Symbols and Abbreviations

A Aspect Ratio, wing VIAS Indicated airspeed, kts AOA angle of attack, deg Vre flap extension speed alt altitude, ft Vs stall speed C14e elevator hinge moment coefficient V. Vprobe tailplane velocity, kts

Cl_Tail tail section coefficient of lift C.. A/C pitching moment coefficient Cp pressure coefficient Greek: CT thrust coefficient _, alpha angle of attack, deg CG aircraft center of gravity t, beta angle of sideslip, deg FAA Federal Aviation Administration SE, deIE elevator deflection angle, deg FYE yoke force, lbs 54, deIA deflection angle, deg G acceleration due to gravity SR, deIR deflection angle, deg IRT NASA Icing Research Tunnel 6F, deIF flap deflection angle, deg MAC mean aerodynamic aelT trim tab deflection angle, deg MTE mission task element E downwash, deg NASA National Aeronautics and Space O, theta pitch angle, deg Administration N_ normal acceleration Subscripts: OSU Ohio State University a/c aircraft Pbeta tailplane sideslip angle L,I total lift, section lift q pitch rate, deg/s tail, t horizontal tailplane ROD rate of descent, ft/s i refers to one of three tailplane S&C stability and control 5-hole probes: 1 = outboard, SHSS steady heading sideslip 2 = midspan, 3 = inboard TAOA tailplane angle of attack, deg TED trailing edge down TEU trailing edge up

Aircraft Axis System and Sign Convention:

Stabilizer and Elevator Sign Convention:

+ Lift

"...... nunzontal _taD,llZt_r i_ -_

+ 8t + 8E

- Lift

NASA/TP--2000-209908 ix

1.0 Introduction

Aircraft accident analyses have revealed ice contamination on horizontal as the primary cause of 16 acci- dents resulting in 139 fatalities. Ice can lead to a premature tail stall that causes the aircraft to pitch nose-down, which at low altitude, may not be recoverable prior to impact with the ground. Three International Tailplane Icing Workshops were convened to appraise the collective experience on ice-contaminated tailplane stall (ICTS) from airframe manufacturers, operators, aviation regulators, and other interested parties. Workshop attendees provided recommendations to reduce the number of accidents attributed to ICTS. In response to some of these recommenda- tions, the Federal Aviation Administration (FAA) requested the National Aeronautics and Space Administration (NASA) to conduct research into the characteristics of ice-contaminated tailplane stall and to develop techniques and methodologies to minimize the hazard.

NASA responded to the FAA request by developing the NASA/FAA Tailplane Icing Program (TIP). _ The TIP was a four-year research program that utilized a combination of icing experts and test facilities that included the NASA Glenn (formerly NASA Lewis) Icing Research Tunnel (IRT), The Ohio State University (OSU) Low Speed Wind Tunnel, and the NASA Glenn DeHavilland DHC-6 Twin Otter Icing Research Aircraft. These resources were used to accomplish the following program goals: 1) improve understanding of iced tailplane aeroperformance and aircraft aerodynamics, 2) develop analytical tools to help discriminate tailplane sensitivity to icing, and 3) develop training aids to expand the awareness of the ICTS aviation hazard.

The purpose of this report is to document the flight data results acquired during this four-year program. Subsequent reports will document results from the wind tunnel tests and analytical code development.

NASA conducted extensive flight tests using its Twin Otter Icing Research Aircraft with three simulated ice shapes attached to the leading edge of the horizontal tail, and a clean tail baseline. Testing the Twin Otter was advantageous because of low operating costs, adaptability for the intended experiment, and known flight characteristics from pre- vious iced stability and control experiments. 2, 3 The Twin Otter also has a known ice-contaminated tail stall suscep- tibility. Fifty-one research flights were conducted in the two-phase flight test program (Table A-1 and Table A-2). Flight test maneuvers with each ice shape included 1G wings level points; steady-heading sideslips; pushovers to 0.5G, 0.25G, and 0G; elevator doublets; and thrust transitions. These maneuvers were conducted for a full range of flap positions and aircraft angles of attack until tail stall was imminent. Over 2000 test points were collected. Spe- cific test points were selected for an in-depth analysis and are reported in this document.

This report is organized in the following sections: description of the research aircraft, instrumentation systems, and three data analysis sections: steady state, dynamic maneuvers, and parameter estimation. Within each of these sections, the flight test procedures, data analysis methods, results, and conclusions are presented. Concluding remarks on the analysis methods and maneuvers are offered.

NASA/TP--2000-209908 1

2.0 Research Aircraft and Ice Shapes

The NASA Icing Research Aircraft is a modified DeHavilland DHC-6 Twin Otter (Appendix A, Figure A-l). The horizontal tailplane has a fixed stabilizer with an elevator and trim tab. It is powered by two 550 SHP Pratt and Whitney PT6A-20A turbine engines driving three-bladed Hartzell constant speed propellers. The flight controls are mechanically operated through a system of cables and pulleys. Control surfaces consist of elevator, , rudder, and wing flaps. Physical characteristics of the aircraft are listed in Table A-3.

The Twin Otter was tested in a Baseline configuration (clean tail) and with three simulated ice shapes. One ice shape represented an Inter-Cycle ice accretion (Figure A-2), the second represented a Failed Boot ice accretion (Figure A-3), and the third represented a double-horn, glaze ice shape (Figure A-4). The Inter-Cycle and Failed Boot shapes resulted from an Icing Research Tunnel (IRT) test on a full scale Twin Otter tailplane model. These ice shapes were urethane casts made from molds of the ice accreted during the IRT test and retained the overall shape and rough texture of the actual ice shape. The third shape was used extensively in NASA's previous stability and control tests and was named the S&C ice shape. The S&C ice shape was cut from Styrofoam blocks and did not incorporate surface roughness or 3D effects. Each of the simulated ice shapes was attached to the leading edge of the horizontal stabilizers only. No other surfaces were contaminated.

NASA/TP--2000-209908 3

3.0 Instrumentation System

To gain a complete understanding of the parameters that lead to ice-contaminated tailplane stall, the aircraft was instrumented to acquire three distinct types of data: l) aircraft dynamics, 2) tail aeroperformance and 3) tailplane flow visualization and pilot visual and tactile cues. Each type of data and instrumentation used to acquire it is discussed in detail below.

3.1 Aircraft Dynamics System

The aircraft dynamics data set included: inertial data, air data, control surface deflection data, pilot forces, and engine parameter data. The inertial data consisted of three orthogonally-mounted linear accelerometers, three orthogonally- mounted rate gyros, and a vertical gyro to provide pitch and roll angle data. All sensors were aligned with the body-axes of the aircraft, and contained in a single box near the aircraft center of gravity (CG). Yaw angle was provided by the aircraft directional gyro. Instrumentation specifics are listed in Table A-4.

Air data consisted of airspeed, angle of attack, angle of sideslip, pressure altitude, and outside air temperature (OAT). All air data parameters (except OAT) were sensed by a Rosemount 858 probe head extended from the aircraft on a 9- foot noseboom.

The control surface deflections, _SE,_A, 15R,5F, were measured using linear control position transducers (CPTs) located near the control horns, which eliminated cable stretching errors. Pilot forces were measured using a custom-designed yoke wheel and rudder pedals that incorporated load cells to measure forces required to deflect the ailerons, elevator and rudder. 4

Engine data consisted of torque pressure, propeller RPM, and fuel flow from both engines. By recording fuel flow, the aircraft mass and CG were adjusted for fuel burn during the test flights. All aircraft dynamic sensor data signals (except engine torque pressure and RPM) were passed through a Precision 6000 Filters System unit for conditioning, digitized and recorded using a Science Engineering and Associates (S.E.A) data acquisition system. The engine torque pressure and RPM data were digitized and recorded directly to the S.E.A. system without the signal conditioning of the Precision 6000 system.

3.2 Tail Aeroperformance System

Tail aeroperformance data consisted of tail inflow angles and velocities as well as surface pressures. To measure the former, three 5-hole probes were mounted to the leading edge of the left-side horizontal tail (Figure A-5). One probe was mounted near the tail tip, the second was mounted mid-span, and the third was mounted near the tail root. To measure the latter, a pressure belt wrapped chordwise around the stabilizer and elevator.

Each 5-hole probe measured the local angle of attack, angle of sideslip, and airspeed at the designated spanwise location. The probes were a variation on a design used by Raytheon for the Beech Starship _ but modified by increasing the probe size to reduce vibration errors (Figure A-6). Each 5-hole probe was extensively calibrated at NASA and OSU to determine angle of attack, sideslip, and dynamic pressure over the anticipated ranges encountered in flight. These 5-hole probes related pressure sensed at holes on the probe tip to the relative incoming velocity and flow angle.

On the aircraft, the pressure on each hole (18 total) was transmitted through tubing to two Scannivalve electronic scanning differential pressure transducers. From the Scannivalve transducers, analog signals were transmitted to an S.E.A. pressure multiplexer interface unit, where the flow probe signals were digitized and recorded with the aircraft dynamic parameters.

NASA/TP--2000-209908 5 The other aeroperformance data consisted of chordwise distributed surface pressures around the horizontal tail. These data were obtained using a pressure belt made from a series of strip-a-tube tygon tubing wrapped completely around the horizontal tail in a mid-span location (Figure A-5). The end of each tube on the tail was plugged. Each tube in the belt had a single hole cut at a specific chordwise location where the surface pressure was sensed. The other end of each tube was routed to Scannivalve electronic scanning pressure transducers. In this way, the surface pressures about the tail were sensed at 50 holes on the pressure belt and measured by the Scarmivalve unit inside the aft section of the aircraft. Like the flow probe data, the Scannivadve transmitted the pressure belt data to the S.E.A. pressure multiplexer interface unit, where the signals were digitized and recorded with the other parameters. A theoretical analysis of the pressure belt and flow probe system dynamic response characteristics indicated an approximate lag of 0.2 sec for step input changes.

3.3 Video Systems

Flow visualization on the tailplane was accomplished by attaching yarn tufts to the lower surface of the horizontal tail and mounting a video camera to the bottom aft section of the (Figure A-7). The camera field-of-view could be adjusted by research personnel through a tilt/roll mounting mechanism during the flight. The video lens also incorporated a zoom capability to provide a close-up or far-fidd perspective. The video signal was annotated with real- time engineering unit values of specific parameters. This annotation improved the information content of the video by informing the viewer of information such as airspeed, thrust coefficient, elevator deflection and pilot forces. Lastly, the video signal with annotation was recorded jn SVHS format.

Another unique video system was instailed to record the pilot actions during the maneuvers and also record the view through the windscreen to obtain the pilots perspective. These two views were merged onto a single screen format by using a screen splitter so that the upper part of the screen showed the view through the windscreen, while the lower part of the screen presented an over-the-shoulder look at the pilot controlling the aircraft. This single screen presentation was annotated with engineering unit data to indicate the aircraft pitch and roll angles, pilot forces, thrust coefficient and elevator angle, and then recorded in SVHS format.

3.4 Calibration Procedures

Extensive calibrations were conducted on all components of the data systems. Prior to the flight seasons, individual calibrations were conducted for all sensors, and through-put calibrations (sensor-filter-data acquisition-recording) were performed prior to each flight. Specifically, the Scarmivalve pressure transducers and all 32 analog channels underwent a three-point check prior to each flight. The transducer coefficients obtained from these checks were used in the post-flight data reduction process.

NASA/TP--2000-209908 6 4.0 Data Handling Process

All together, the data from the flight program consisted of the three video signals and 95 data signals. The data signals were recorded at 100 Hz sampling frequency and had 16-bit resolution. These data were recorded onto 8-mm tape using a ruggedized PC compatible data acquisition system. The data on the 8-mm tapes were converted to ASCII files post- flight for further processing on ground-based computer systems. The initial ASCII files consisted of row and column data for each sensor voltage. A FORTRAN program was written to read these files, apply the calibration conversion equations for each sensor, and perform correction algorithms for sensor position errors and other relevant calculations for data analysis. The resultant files were uploaded into Excel for plotting and other data analysis to be described in the subsequent sections.

NASA/TP--2000-209908 7 ? 5.0 Steady State Analysis

Steady flight test points enabled an in-depth examination of the effects of tailplane icing on tailplane lift and hinge moments, and elevator deflections required to trim the aircraft. The parameter space examined in this study varied the tailplane configuration [clean, Inter-Cycle ice, Failed Boot ice, S&C ice], the flap deflection - 6F, the thrust set- ting - Cr, the aircraft airspeed - VIAS (hence angle of attack - _vc), and the aircraft sideslip angle - ft. A compre- hensive matrix of test points were conducted with a forward center of gravity (CG) location. In addition, a limited number of test points were conducted with an aft CG location. The following section details the flight test proce- dures, method of analysis, shows results and provides conclusions.

5.1 Flight Test Procedures

5.1.1 Steady-Wings-Level (SWL): The aircraft was configured to a specific tail condition (clean or iced) and flap deflection angle (fF), and then the thrust setting (Cr) was set at the appropriate test speed. The aircraft was trimmed at this initial test speed, so that the elevator hinge moment was hulled. With the aircraft in a steady condition, data records were taken for approximately 15 seconds. Due to the fixed thrust setting, the aircraft was sometimes in a steady climb, level flight, or in a steady descent. After the initial test point, the next test speed was reached by adjusting the yoke position and resetting the throttles to obtain the consistent Cr. The elevator was not retrimmed, so that the yoke force required to hold the ele- vator in that position could be translated into an elevator hinge moment (Cne). This procedure was repeated for six airspeeds at each thrust setting, flap deflection and tail ice configuration. All of these test points were done with minimal sideslip (13= 0) on the aircraft.

5.1.2 Steady-Heading Sidesfips (SHSS): A limited number of steady-heading sideslips (SHSS) test points were obtained. The procedure was similar to the steady-wing-level points, except that the aircraft was put into approximately a 17 ° sideslip to the right and to the left, while maintaining a specified airspeed, flap angle and thrust coefficient. The sideslip was accomplished by the pilot applying either right or left rudder and cross-controlling with left or right aileron to yaw and roll the aircraft into a steady-heading sideslip.

5.2 Method of Analysis for Steady State Data

5.2.1 General Each 15 second data file was converted and corrected as stated in the data handling process section. For each meas- ured and calculated parameter, the average and standard deviation were calculated. These averaged values were then input to a spreadsheet program for plotting and analysis.

5.2.2 Data Consistency Check

Data consistency check points were added to the test matrix after the fifth flight in 1995 ('FLT 95-17). These steady state points were performed with a consistent aircraft configuration and flight condition, i.e., 8F = 10°, VIAS = 100knots, and Cr = 0.10 or 0.05. These check points were analyzed using the steady state processing techniques described above to determine flight-date-dependent, systematic errors in various instruments. These in-flight check points enabled the discovery of problems in the pressure instrumentation for the tail flow probes and pressure belt. Selected parameters from the data consistency results are shown in Figure B-1 and Figure B-2 in Appendix B.

NASA/'TP--2000-209908 9 5.3 Results

5.3.1 Basic Aircraft Performance Description

The total lift coefficient of the test aircraft was measured for flap deflections [0, 5, 10, I5, 20, 25, 30, 35 and 40°], and is displayed in Figure 1. Clearly, the effect of increasing angle of attack within a specific flap deflection in- creases total lift coefficient, and the effect of increasing flap deflection shifts the lift curve up and to the left. This is well known and understood.

Baseline Baseline

CL m/c V8 (X a/c CL m/cV80_ a/c 3.5 3.5 T 3 3 I _ 1 i t 2.5 2.5 i '¢' .... ,Ira dF--O ) 2 ,i dF=lO dF=5 .a o (11dF=-20 dF=15 i 1.5 i i dF=30 I_ dF=25 J _¢, dF=40 o dF=35j 1

0.5 ) 0 o! -10 -5 0 5 10 15 -10 -5 0 5 10 15

(x (deg) a (dog)

Figure 1. Total lift coefficient vs. aircraft angle of attack for various flap settings.

As a typical aircraft transitions from cruise to landing configurations, the path along the CL VS. angle of attack curve looks similar to Figure 2.

Beginning at point 1, the cruise speed is reduced to flap extension speed (VFe) or slower by reducing thrust and increasing angle of attack (point 2). To maintain altitude, the aircraft CL also increased. Flap deflection drives the situation from point 2 to point 3 causing a significant reduction in aircraft angle of attack. Next, the speed is again decreased by reducing thrust and increasing angle of attack (point 3 to 4). When full flaps are selected (point 4 to 5), the minimum angle of attack is attained assuming the transition took place at VFE. As the aircraft airspeed is further reduced for landing, the lift condition and angle of attack increases as shown in point 5 to 6.

This scenario depicts some of the benefits of wing flaps. Flaps allow an aircraft to land at reduced airspeeds, which in turn reduces the runway length required for landing. Flaps also reduce the angle of attack required for a (eeg) high lift situations, which provides better forward visibil- ity to the pilot during landing. Figure 2. Typical lift curve during landing transition

NASA/TP--2000-209908 10 However,withthebenefitsofflaps,therearedrawbacksaswell.Sinceflapsincreasethewinglift byincreasingthe circulation,thereisanincreaseinthedownwashangleintothehorizontaltail.Flapextensionalsoincreasesthenose- downpitchingmomentbecauseof both the increased wing lift vector and the increased moment arm due to increased chord length. Figure 3 indicates how the pitching moment increases with flap deflection, and how the tail down lift must be increased to balance this moment. The double lines indicate vectors after flap deflection. Note the magni- tudes of the vectors in Figure 3 are for illustration purposes only.

WING LIFT

TAIL LIFT

WEIGHT

Figure 3. Force and moment diagram.

The function of the horizontal tail is to stabilize the aircraft in pitch. Therefore, when flaps are deployed, the nose- down pitching moment due to increased wing lift and aft movement of the wing center of lift requires the tail to pro- vide an increased down lift to balance the pitching moment. The increased downwash angle from the wing assists the horizontal tail in performing its function by increasing the tail angle of attack. This increase in tail angle of attack may be too much or too little for the particular trim airspeed, in which case the elevator angle is adjusted by the pilot to either decrease or increase the tail down lift.

NASA/TP--2000-209908 11 5.3.2 Tailplane Angle of Attack Description

The test aircraft was used to measure the effects of wing flap and aircraft angle of attack on the tail angle of attack. The effects of aircraft angle of attack and downwash can be seen in Figure 4 for the clean tail configuration. Note that the lines in the figures are drawn to highlight trends; they are not meant for quantitative analysis.

Baseline

tail2 VS(7. Mc

5_ ......

i oi I

I O dF=20 _ ! -10

-15 t--- j I I I -20 t 1 -10 -5 0 5 10 15

c((deg)

Figure 4. Tailplane angle of attack vs. aircraft angle of attack.

For a given flap deflection, it is clearly seen that as the aircraft angle of attack was decreased, the tailplane angle of attack also decreased. For the test aircraft in a given flap setting, this change in aircraft angle of attack was greater than the change in tailplane angle of attack such that the slope Aat, it / Aot,g_ < 0.6 for 6F = I0 °, and much less than 0.6 for the other flap settings. Likewise, aS_the flap deflection was increased, there was a bias shift in the tailplane angle of attack in the negative direction. This bias shift was due to the increase in wing downwash. The largest in- cremental downwash contribution due to flap deflection was observed when flaps were extended from 0-10 °. This caused approximately a -5 ° shift in tailplane angle of attack. As flaps were further deflected, the incremental change in tailplane angle of attack decreased such that when flaps deflected from 30-40 °, the shift in tailplane angle of at- tack was less than -1 °

The relationship between tailplane angle of attack, aircraft angle of attack and downwash is described in Ref. 6 as the following equation: o_ = a,e,.- e- i, where _ = tailplane angle of attack a_/_ = aircraft angle of attack t; = downwash angle it = tail incidence angle (fixed at O°for test aircraft)

We also know from Ref. 7 that wing downwash is a function of lift coefficient (CL). In theory, the downwash at the i i4.6C L tail is approximated by (A is the wing aspect ratio) e ° = zrA

NASA/TP--2000-209908 12 Astheaircraftangleof attack(aojc)decreases,thereis alsoa decreasein downwashangle(e),suchthatif the changesinthesetwoparameterswerethesame,therewouldbenonetchangein thetailalpha.Butweseein this datathatastheaircraftalphaandcorrespondingdownwashdecreasewithinagivenflapsetting,thatthetailalpha alsodecreases.Thuswithinagivenflapsetting,theaircraftalphadominatesthetailalphaequation. Analternativewayoflookingatthevariationin tailalphawiththeflapconfigurationiswithrespecttoaircraftair- speed(Figure5).Forallflapsettings,anincreaseinairspeedresultsinamorenegativetailangleofattack.

Baseline (_ tail 2 vs VTAS e/c

• dF=O

-5 z, dF=lO , • dF=20

q -10 i _, dF=30 I ' i .odF=40 -15 1 -20 I 2O 40 60 80 I00 120 140 160

VTAS (kt)

Figure 5. Tailplane angle of attack vs. airspeed.

Tying this information together with the typical transi- tion from cruise to landing scenario developed above (see Figure 2), we see in Figure 6 that the tail angle of attack is greatly affected by both the flap transition and aircraft angle of attack changes. In the same way as before, the cruise aircraft (point 1) is slowed to VFE o (point 2). Then by selecting the first flap setting of 10°, the tail AOA is reduced from 0 ° to -8 ° (point 3). The -5 tail AOA increases as the aircraft is slowed further to -10 the next flap VFe (point 4), and then reduced from -6 ° to -13 ° as the flaps are fully deflected (point 5). As the aircraft flares for landing, the tail AOA again increases -15 for landing (point 6). -2O The relevance of this information is not in the exact -10 -5 0 5 10 15 values of tailplane AOA, as that will vary from aircraft a (deg) to aircraft. The relevance is in the relationship between the aircraft angle of attack, flap setting and tailplane Figure 6. Typical tail AOA during landing transition. angle of attack. This information is positive evidence that as the aircraft angle of attack is reduced, or as the aircraft is at or near VFE , the tailplane AOA is at the most negative value for that flap setting. Therefore in non- maneuvering flight, maximum flap deflection and lowest aircraft AOA or VrE (point 5) clearly subjects the tailplane to the most negative angle of attack. This condition is where the tailplane stall margin is a minimum. This is exactly opposite of the wing stall scenario where slow speeds and high aircraft angle of attack put the wing near its stalling

NASA/TP--2000-209908 13 angle. Therefore, in steady-state (non-maneuvering) flight, the pilot has two stalling limits. A low speed limit for the wing, and a high speed limit for the tail.

5.3.3 Tailplane Lift Description

As stated in Section 5.3.1, the function of the horizontal tailplane is to provide stability and control in the pitch axis. Conventional horizontal tails incorporate a fixed stabilizer and hinged elevator with trim tab on the elevator (Figure 7). Pitch stability and control is accomplished by modifying the amount of lift generated at horizontal tail by moving the elevator. The elevator cambers the tail to modify the amount and direction of the lift. Conventionally, a positive elevator deflection moves the trailing edge down (TED) and reduces the down lift. A negative elevator deflection moves the trailing edge up (TEU) and increases the down lift. Previous wind tunnel tests documented in Refs. 8 and 9 show how the elevator effects the magnitude of lift and stalling angle of attack. Data from Ref. 9 are shown in Figure 8 where the lift coefficient becomes more negative and the stall angle of attack becomes less negative (moves to the right) as SE decreases. Also, for a given elevator deflection (t_E), the tail lift becomes more negative as the tail AOA becomes more negative. This behavior is similar to a wing, but inverted.

+ Lift

t (C"_ _...... Horizontal Stabilizer levator

+ St + _E + CHe

- Lift

Figure 7. Horizontal tail diagram.

Baseline Tail Lift Coefficient

os 4

J _0_

-LE I t

-20 -15 -10 -5 0

Angle Of Attack (degrmm)

Figure 8. Tail Hft coefficient from wind tunnel data.

NASA/TP--2000-209908 14 SimilartothewindtunneltestinginRef.9,thetestaircraftwasequippedwithapressurebeltontheleftsideofthe horizontaltailasdescribedin Section3.2.Thepressurebeltenabledtheresearcherstoobserveandmeasurethesur- facepressures(Figure9)andcalculatethesectionlift coefficientbyintegratingtheCp distribution.

Baseline, dF = 20, Ct = 0.1 xcAifoil CP L_ ! ,,, .Lower erSuSurface ace

-0.2 0 0.2 04 0.6 O.B 1 xJc 1.2

Figure 9. Cp distribution for Baseline tail.

Figure 10 shows the tail lift coefficients from the flight data plotted versus tailplane angle of attack measured at the mid-span (probe 2). All test points in Figure 10 were with a forward CG location (22% MAC) and moderate thrust setting (Cr = 0.10).

Baseline C I tail V8 (1 tail 2

o • dF=0 -o.1 ...... -o.2 I & dF=10 -0,3 • dF=20 m -o,4 d -05 II dF=30 --- I / I i -o.6 o dF=40 -o.7 _1 / I -0.8 - I I ! -14 -12 -10 -8 -6 -4 -2 0 2 _x 1=*1(deg)

Figure 10. Ct_roit vs. c_il from flight data.

Another way of representing these data are with respect to aircraft airspeed. Figure I 1 shows the variation in tail download required to maintain specific airspeeds at various flap settings.

Baseline C_tail Required for Speed c.g.=22% MAC; C --0.10 0 -0.1 IF-----, • '. -0.2 dF=l 0 -0.3 -0.4 dF--2.0 (5 dF=O -0.5 : I dF=30 -0.6 t -0.7 t Lo dF---40 -0.8 ! 60 70 B0 90 100 110 120 130 140 150 160 vTAS (knots)

Figure 11. CLT=Uvs. VTAS from flight data.

NASA/TP--2000-209908 15 OneofthekeyobservationsfromthesedatapresentedinFigureI0andFigure11istheincreaseddownloadrequired forsteadyflightastheflapswereextended.Forthetestaircraft,therewasasubstantialincreasein themagnitudeof tailsectionlift coefficientfromthenoflaptofullflapextension.Withinagivenflapsetting,thetaillift coefficient variedwiththeflightcondition.Especiallyforthelargerflapdeflections,theCLraiI required to stabilize the aircraft was most negative at the higher angle of attack or slow speed flight.

Another key observation was that at the larger flap deflection cases, the CLr_il magnitude decreased as tail AOA be- came more negative. This seems counterintuitive to the wind tunnel results in Ref. 9, but the difference is because the wind tunnel data were with constant elevator deflection angles (BE), whereas the flight Cl_r_ildata presented in Figure 10 had various elevator deflections based on the tail download requirement. Since the aircraft was held in a steady wings level condition, the elevator was adjusted to hold the aircraft at the specified airspeed. In the higher speed conditions, the elevator was held in a TED position, thereby unloading the horizontal tail. In order to hold the low speed points, a greater down load is required from the tail. For that condition, the elevator was held in the TEU posi- tion to create a higher down load and slow the aircraft. The dashed line in Figure 10 indicates the approximate loca- tion of _E = 0°. The arrow pointing to the upper left indicates the elevator position with the trailing edge down, and the opposing arrow indicates the elevator in a trailing edge up position.

Figure 12 compares the wind tunnel to the flight test results. Flight data from the _F = 40 ° case are superimposed on the wind tunnel data with the elevator deflection angle noted. Also, a dashed line has been added to the plot to indi- cate the stalling angles of attack for each elevator angle. Note that the wind tunnel data are from a 2-D test and have not been corrected for 3-D effects of flight. For these reasons, the actual flight and wind tunnel CL values for a given elevator deflection do not agree well.

Tail Lift Coefficient

1 -*-dE:l,2 I I os - _/--=-dE=-10 t_...... j _.____J

.... Flight, dF= 40 -_o- dE= -20 F __k 0 I - L

-0 5 - _a_ _'_:-_-_ ± ° a _'_ _.*J _ jp/I

-_0 _15 -10 -$ o

Angle Of Attack (degrees)

Figure 12. Tail lift coefficient from wind tunnel data and flight data.

Figure 12 does illustrate, however, relative stall margins. Even though slow speed flight test points required a greater tail down load (Ct_r_il = -0.7), the tail AOA was well within the linear range of Cl_r_i_ for that particular elevator de- flection (tE - -lC). On the other hand, the high speed test point required less tail down load (CLr_i I = -0.55), but subjected the tailplane to an angle of attackcloserto the stalling angle for the respective elevator deflection (SE = 7). One can see that the high speed point is closer to the stall angle line than the low speed point.

Although this comparison of windtunnel data to flight data is not fully developed, it is used to illustrate the decrease in tail stall margin as flaps are deployed andthe aircraft is flown at speeds near the flap extension speed.

NASA/TP--2000-209908 16 5.3.4 Elevator Hinge Moment Description

Recall from Section 5.1.1 that the steady state points were initiated with the elevator trimmed using the pilot- controlled elevator trim tab for the specified speed (Figure 13). This initial test speed was typically near 1.1Vs, so the elevator was positioned TEU. In order to hold the elevator in that position, the trim tab was adjusted TED. The aero- dynamic forces over the trim tab provided a hinge moment to counteract the hinge moment of the whole elevator. In this way, the elevator force felt by the pilot was zeroed for the initial test point.

4 T HorizontalStabilizer] E)e_ator ] ITrim Tab

Figure 13. Horizontal stabilizer, elevator and trim tab.

Test points following the initial test point were obtained by the pilot positioning the yoke more forward to reach and maintain the next higher test speed. The force applied by the pilot was measured and transformed into the elevator hinge moment coefficient through the relation below:

CH, = g qtail Se Ce

where: g = elevator gearing ratio = 0.073 lb/in-lb for test aircraft (obtained experimentally) F, = pilot stick force to move elevator q t=t= dynamic pressure at the tailplane S, = area of elevator aft of hinge line ce = elevator chord aft of hinge Iine

As the pilot set the yoke more forward, the elevator moved to a more TED position. Typically, the pilot provided a push force to move the yoke forward which resulted in a negative elevator hinge moment. The negative hinge mo- ment on an elevator that is being moved TED is a restoring moment. The restoring moment would return the elevator to the trimmed position and the aircraft would return to the trim speed if the pilot were to release the controls.

The useful way of expressing elevator hinge moment from Ref 11 is below. The terms in this model will be used to help describe the icing effects on tailplane performance in later sections.

Crte = b0 + bl0q + b2_ + b3(_'r

where: oq is the change in tail angle of attack is the change in elevator deflection angle is the chance in trim tab angle bo is the value of the hinge moment when txt= _ = _= 0 b_ is referred to as the elevator "floating" parameter and is usually negative, but can positive by design. If negative, the elevator floats with the wind when the angle of attack is changed. b2 is referred to as the elevator 'qaeaviness" parameter since it relates the amount of pilot force required to deflect the elevator. A preferred control characteristic is to have the hinge moment acting to oppose the mo- tion of the control surface, so the term b2 should always be negative. If b2 is positive, the control surface would be considered "overbalanced." b3 is the trim tab effectiveness parameter and is typically positive so that the b38-rwill oppose the b28t_ term.

NASA/TP--2000-209908 17 Figure14andFigure15showthetotalelevatorhingemomentCue for the Baseline aircraft with respect to tail AOA and aircraft airspeed.

The test aircraft in the baseline configuration experienced the restoring moment for all flap settings. Note that in the _F = 40 ° case, the aircraft was not able to be trimmed at the slow speed test point and the pilot had to pull slightly (3 lb.) to maintain the 55 knot indicated airspeed. This suggests that the pilot-controlled elevator trim tab was unable to produce enough hinge moment (magnitude of b3_ was too small) to counter the effective elevator hinge moment caused by the elevator deflection angle (b2_ term).

Baseline CHe V80_ tail 2 0.02 ...... i ...... i"...... I ...... " o ..... t I a • dF=-O -0.01 -0.02 ,t, dF=lO

t_ -0.03 • dF--20 -0.04 -o.o5 dF--30 -0.08 -0.07 k o dF--40 -0.08 l I ...... t -14 -I2 -10 -8 -6 -4 -2 0 2 4

OC=. (deg)

Figure 14. Elevator hinge moment coefficient vs. tail angle of attack.

Baseline CHe vsTrue Airspeed 002 0.01 0 _e dF=O -0.01 -0.02 -0.03 -0.04 dF=20 -0.05 i dF=lOdF--30 -0.06 -0.07 i o dF--40 -0.08 6O 70 80 90 100 110 120 130 140 150 160

VTAS (k nots)

Figure 15. Elevator hinge moment coefficient vs. airspeed.

NASA/TP--2000-209908 18 5.3.5 Icing Effects on Tailplane Lift

Using one of the techniques in Section 5.3.3, the effects of various levels of icing on tailplane lift are examined in this section. Since the effects of the various levels of icing are of interest, the data series on each plot consists of the clean baseline, Inter-Cycle Ice, Failed Boot Ice, and the S&C Ice. The tailplane section lift coefficient versus aircraft true airspeed data are presented in Figure 16 through Figure 20. Each figure is for a specified flap deflection.

Tail Lift Coefficient Required for Speed 8F=0; c.g.--22%MAC C_--0.10 0 -0.1 • Baseline -0.2 [] Intercycle -0.3 1 O- -0.4 Failed Boot -0.5 ,L S&C -0.6 -0.7 ! 60 70 80 90 100 110 120 130 140 150 60

VTAS (knots)

Figure 16. Tail lift coefficients for _F = 0 °.

Tail Lift Coefficient Required for Speed 5F=10; c.g.--22% MAC C_---0.10 0 -0.1 • Baseline -0.2 -0.3 [] Intercycle o- -0.4 Failed Boot -0.5 A S&C -0.5 -0.7 60 70 80 90 1O0 110 120 130 140 150 160

VTAS (knots)

Figure 17. Tail lift coefficients for 8F = 10%

NASA/TP--2000-209908 19 Tail Lift Coefficient Required for Speed _F=20; c.g.=22%MAC CT=0.10 0 -0.1 • Baseline -0.2 -o.3 ill Intercycle o- -0.4 _ Failed Boot -0.5 -0.6 1 S&C -0.7 11 60 70 80 90 1O0 110 120 130 140 150 160

VTAS (knots)

Figure 18. Tail lift coefficients for o_F= 20 °.

Tail Lift Coefficient Required for Speed 8F=30; c.g.=22%MAC C_-0.10

-ol __I 1 • Baseline -0.2 1 -0.3 ...... m _tercycle o- -0.4 _Failed Boot S&C -0.7 60 70 80 90 100 110 120 130 140 150 160

VTAS (knots)

Figure 19. Tail lift coefficients for 8F = 30 °.

Tail Lift Coefficient Required for Speed 8F--40; c.g.--22% MAC C_--0.10 0 -0.1 I I -0.2 J I .... I • Baseline -0.3 i I I 1 _ _- -0.4 ill Intercycle 0 -0.5 I 1 -0.6 Failed Boot -0.7 I -0.8 ! i 60 70 80 90 100 110 I20 130 140 150 16O

VTAS (knots)

Figure 20. Tail lift coefficients for 8F = 40 °.

NASA/TP--2000-209908 20 Althoughthereweresomevariationsinthelift coefficientsdevelopedwiththeicecontamination,thegeneralamount of downloadgeneratedforthegivenflapdeflectionaresimilar.Thisisnotsurprisingsincetheaircraftrequiresa specificamountofdownloadfromthetailfortheflightconditionandaircraftconfiguration.NotethattheS&Cice shapewasnottestedin thetSF = 40 ° configuration because the tailplane could not stabilize or control the aircraft adequately,

Having similar tail down loads did not mean the tailplane was not affected by the ice contamination. The lift loss due to the ice was accounted for by deflecting the elevator and cambering the tail more than in the clean baseline cases. Figure 21 through Figure 25 illustrate how much elevator deflection was required for the iced cases and various flap configurations. In the SF = 0° case, there is virtually no difference between the iced cases. As the flaps extend to 10°, the effects of the ice can be seen in the 1-2 ° decrease in 6E from the baseline case.

Elevator Deflection Required for Speed 8F=0; c.g.=22%MAC CT=0.10 68 ...... I I I I • Baseline ,_, 4 m Intercycle "o .... O_-IlI_L-q _l-iir: _; 4_Failed'Boot ILl -2 '0 -4 -6 I -8 I I .... 60 70 80 90 100 110 120 130 140 150 160

VTAS (knots)

Figure 21. Elevator deflection for 6F = 0%

Elevator Deflection Required for Speed aF=10; c.g.=22%MAC C_--0.10

4 ...... i...... -e I,TI...... t...... Ii...... • Baseline II Intercycle Failed Boot Ill A S&C

60 70 80 90 100 110 120 130 140 150 160

VTAS (knots)

Figure 22. Elevator deflection for 8F = 10 °.

NASA/TP--2000-209908 21 Elevator Deflection Required for Speed 8 F=20; c.g.=22% MAC CT=0.10 8

• Baseline r [] lntercycle _ Failed Boot __!, -1'5 " =o -4 i I ,_ S&C -6 1 , i , I -8 i i i p ..t I 60 70 80 90 1O0 110 120 130 140 150 160

VTAS (knots)

Figure 23. Elevator deflection for b'F = 20%

Elevator Deflection Required for Speed 5F--30; c.g,--22% MAC C_-0.10

e Baseline ti htercycle e 1 uJ -2 .... i ...... _=' -4 -6 t t l _ S&CFailed B°°t -8 ! .! 60 70 80 90 100 110 120 130 140 150 160

VTAS (knots)

Figure 24. Elevator deflection for &F = 30 °.

Elevator Deflection Required for Speed 5F_10; ¢.g.=22% MAC CT--0,10 8

Ir'---C73: ' I Intercycle ] 1 I! FailedBaselineBoot -6 J -8 i I 60 70 80 90 100 110 120 130 140 150 160

VTAS (knots)

Figure 25. Elevator deflection for b'F = 40 °.

NASA/TP--2000-209908 22 Thetrendcontinuesasflapsarefurtherextended.With_F= 30°,theS&Cicedtailrequiredabout5°moreTEUele- vatortostabilizetheairplaneat80knots.Duetothelargeamountofflowseparationandcontrolforcereversal,the testpointswiththeS&Cicewerediscontinued'beforecompletingtheSF= 30° speedpoints.NotethatfortheS&C icecasein boththe_F= 20° and 30 °, the difference in SE requirement from the baseline increased as the speed in- creased. This again signifies that as the tailplane angle of attack became more negative, the flow separation and lift loss increased, and reduced the elevator effectiveness.

Another interesting observation was the similarity in the Inter-Cycle and Failed Boot cases in these plots. This sug- gests that the level of lift loss in both of these iced cases was accounted for by the same amount of elevator deflec- tion.

NASA/TP--2000-209908 23 5.3.6 Icing Effects on Elevator Hinge Moment

Similar to the last section, the techniques of Section 5.3.4 will be used to demonstrate the effects of various levels of icing on elevator hinge moment. The elevator hinge moment coefficient is plotted with respect to the aircraft true airspeed for each-ice shape case. Figure 26 represents the _F = 0° case; one can see that there is little difference be- tween the ice cases. The Inter-Cycle is slightly offset, but that is due to the test run starting without fully trimming the elevator. Likewise, Figure 27 shows the data from the 6F = l0 ° case. Again an offset is seen but in the baseline case instead of the Inter-Cycle ice. The offset in hinge moment resulted from trimming the baseline case at a higher airspeed. For these flap settings, the ice did not have a strong effect on the hinge moment sensed by the pilot.

1 i Elevator Hinge Moment Coefficient Required for Speed ] 5F=-0; c.g.--22% MAC CT=0.10 0.05 0.04 m_ i iit .... J i ] • Baseline 0.03 I l 0.02 Ill Intercycle J 0.01_ Failed Boot -0.01 -0.02 S&C -0.03 , ...... j ! !::IL ,,_]i 60 70 80 90 100 110 120 130 140 150 160

VTAS (knots)

Figure 26. Elevator hinge moment coefficient for 8F = 0°.

Elevator Hinge Moment Coefficient Required for Speed 5F=-10; c.g.--22%MAC C_--0.10 0.05 0.04 0.03 • Baseline 0.02 i Intercycle J 0.01 4, Failed Boot -0.01 -0.02 _. S&C -0.03 6O 70 80 90 100 110 120 130 140 150 160

VTAS (knots)

Figure 27. Elevator hinge moment coefficient for b'F = 10 °.

NASA/TP--2000-209908 24 Elevator Hinge Moment Coefficient Required for Speed 6F=-20; c.g.--22%MAC C_--0.10

• Baseline 0.02 J lit Intercycle O.Ol ,-- Failed Boot S&C -0.03 - _ i 60 70 80 90 100 110 120 130 140 150 160

VTAS (knots)

Figure 28. Elevator hinge moment coefficient for _F = 20 °.

Elevator Hinge Moment Coefficient Required for Speed 8F=-30; c.g.=22%MAC CT=0.10 0.05 ...... i...... i...... q 0.04 1 0.03 • Baseline / 0.02 i t IR Intercycle -b i J 0.0(_ _ Failed Boot -001 -0.02 _, S&C I' -0.03 ! I i_ -_ 60 70 80 90 100 110 120 130 140 150 160

VTAS (knots)

Figure 29. Elevator hinge moment coefficient for 8F = 30 °.

Elevator Hinge Moment Coefficient Required for Speed _F=-40; c.g.--22%MAC C_-0.10 0.05 0.04 ]----_ 0.03 0.02 • Baseline i J 0.01 II lntercycle -0.01 I ,_, Failed Boot -0.02

-0.03 I 60 70 80 90 100 110 120 130 140 150 160

VTAS (knots)

Figure 30. Elevator hinge moment coefficient for 6F = 40 °.

NASA/TP--2000-209908 25 In Figure28witht_F = 20 °, the effects of the S&C ice became significant. The slope of Cne to VTAS reversed signs causing the final Cm value to be nearly zero. This switch in slope is indicative of a control force lightening. Note this is only lightening instead of a CFR because the Cne did not cross the zero line (initial trim point).

When the flaps were defected to 30 ° (Figure 29), three observations can be made. One observation is that the yoke forces could not be trimmed in all iced configurations. The pilot had to pull the yoke to establish the low speed test point. This indicates that the pilot-actuated trim tab could not produce enough of a moment to counter the hinge mo- ment required to hold the elevator in position for that airspeed. The trim tab creates a counter-moment through the aerodynamic force created over its surface. When the airflow upstream of the trim tab is perturbed, the trim tab is not as effective (the term b3 from Section 5.3.4 is reduced). The second observation concerns the control force reversal caused by the S&C ice. The slope of Cn, to VTAS is positive and increases to a more positive value than the initial test point. In other words, as the pilot moved the yoke forward to increase airspeed, he had to pull back with an in- creasing amount of force to hold the yoke in that forward position. The reason for this is because the yoke was transmitting the force caused by the separation over the lower surface of the elevator. In terms of the math model in section 5.3.4, the b_ and bz terms became positive so that the total hinge moment, Cn, became positive. The pressure distribution about the tailplane caused the elevator to float against the Wind (bt > 0) and to behave as if overbalanced (b2 > 0). The third observation is to note the beginning of control force lightening in the Failed Boot ice case.

The 61== 40 ° configuration is depicted in Figure 30. As in the previous configuration, neither ice nor the baseline cases were trimmable. The Failed Boot ice shape clearly causes a control force lightening and nearly a reversal be- cause the final Cne was very close to the initial Cne.

Comparing the data sets in this way indicates the seriousness of the control issues caused by tailplane icing. In the worst case, the control forces recorded to hold the aircraft steady were on the order of 35 lbs. This was not excessive, but the forces were in the opposite direction from the norm and represent poor control "float" and "heaviness" pa- rameters. Also, the degradation caused by the various ice shapes was clearly observed such that the S&C ice shape was by far the worst, followed by the Failed Boot. As expected, the Inter-Cycle ice had the least aerodynamic effect on hinge moment, but still caused a reduction in trim tab effectiveness.

NASA/TP--2000-209908 26 5.3.7 Center of Gravity Effects on Tailplane Lift

The effect of aircraft center of gravity (CG) on Ct_rait was also examined and the results are shown below. The for- ward CG test points have a solid trend line associated with them, and all the symbols are filled in except ¢5F=40° case. The aft CG test points have similar symbols to designate flap settings, but the symbols are not filled except the _SF = 40 ° case. The forward CG locations were approximately 22% MAC whereas the aft CG locations were about 33% MAC.

The key observation from Figure 31 is the decreased magnitude of C_r_it required with an aft center of gravity. This result was expected; as the moment arm from the wing center of lift to the CG location was reduced, so was the tail

Center of Gravity Effects on Cl tai_ vs ct tail

Baseline, CT=0.I 1, cg=22%MAC or 33%MAC

0.1

v 0.0 [ 0 • dF=O A _ <> ,5 O -0.1 i -02 4- [_ dF=lO -0.3 dF=20 Ill -0.4 dF=30 5 -0.5 0 dF=40

-0.6 0 dF=O -0.7

-0.8 _, dF=10dF=20

-0.9 _ dF=30 -14 -12 -10 -8 -6 -4 -2 2 4 • dF=40 = t,. (deg)

Figure 31. CG effect: CZ_ToiZVS. O_.a.

Center of Gravity Effects on C I tail Required for Speed

Baseline, CT=0.11, cg=22% MAC or 33% MAC 0.1 <>. _ i 0.0 b - i • dF=O -0.1 ,_ dF=lO LA, A • i -0.2

-0.3 • A_ ! ,t_ lh, dF=30 -0.4 .... • Ill dF=30 0 -0.5 O dF=40 -0.6 i 0 dF=O -0.7

-0.8

-0.9 dF=20

50 75 100 125 150 1 dF=30

VTAS (kt$) • dF=40

Figure 32. CG effect: G_r,u vs. VTAS.

NASA/TP--2000-209908 27 download required to hold the aircraft steady. However, for the aft CG case and within a given flap position, the tail down load increased as a_ai_ became more negative. This trend is opposite the one with the forward CG, so further investigation was warranted. The C L r, it data were also plotted versus aircraft speed to gain further insight (Figure 32).

Again the decreased C__rait is observed with the aft CG, and again the divergence in tail lift requirement is observed. Examination of the tail surface pressures verified the trends. These data were consistent with the results shown in Ref. 12. Since the CG is aft of the wing aerodynamic center, the tail down lift needed to be reduced as the wing lift increased to balance the pitching moment. This effect is further demonstrated in Figure 33. For the forward CG test points, all the trend lines (except tSF = 0°), have a negative slope, which indicates more tail down load required as the aircraft lift increases. In the aft CG cases, the trend lines have a positive slope, indicating less tail down load as the aircraft lift increases. The aircraft is still statically stable since the CG is forward of the neutral point; however, there were some flight conditions (SF = 0°, CL,gc > 1.0) where the tail was required to have a positive lift.

Center of Gravity Effects on C_ tai! Required for Trim Baseline, CT=0.11, cg=21% MAC or 33% MAC

0 dF=O

A dF=l 0

• dF=20

-_,dF=30 ,J I_ dF=30

° dF=40

° dF=O

,', dF=lO

© dF=20 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0 13 dF=30 C ! ta II 41,dF=40

Figure 33. CG effects on Ct_roit with respect to CL _c.

The effect of CG on tailplane lift was demonstrated above. Greater down loads were required for the forward CG position than the aft CG position. As the aircraft lift was increased, the tail lift was increased for the forward position and decreased for the aft position. The greatest tail down load required for steady wings level flight was with flaps fully deflected, a forward CG position, and the aircraft in a high lift flight condition (near Vs).

NASA/TP--2000-209908 28 5.3.8 Thrust Effects on Tailplane Flow Field and Performance

The steady-state effect of propeller thrust (Cr) on the flow field at the tail and on the Ct_ra_twere examined and the results are shown in Figure 34 through Figure 36.

As the thrust was increased, and the aircraft velocity held constant (V = 1.3Vs), the tail AOA at the midspan nomi- nally remained constant for a given flap setting (Figure 34). With 6F = 0°, the trend indicates a small increase in negative t;_,_ as thrust increased, but considering the scatter in the data, this was not a significant result. Other con- stant velocity cases were examined and similar results were observed for the midspan flow probe. From the spanwise analysis (next section), there may have been a more substantial thrust effect on o_a, observed at the tail tip. This analysis was not performed on the flow probe near the tail tip (Probe 1) due to initial problems caused by a blocked pressure port.

Thrust Effect on Tail AOA Baseline; V~1.3Vs, c.g.=20%MAC

5

• dF---0 A 0 D'J dF=l 0

-5 , dF=20 J 4 _t A A ,_..=A dF--30

O dF--40 -10 _ -15 , , 0.00 0.05 0.10 0.15 0.20 0.25

CT

Figure 34. Thrust effect on c_u.

Another observation of the thrust effects on the flow field is the impingement of propeller slipstream on areas of the horizontal tail. With 6F = 0° and 10°, there is a clear increase in flow velocity at the midspan probe position (Figure 35). This increase in velocity was attributed to the propeller slipstream. As the flaps were extended greater than 10°, there were only minor velocity increases observed at the midspan location. This indicated that the propeller slip- stream moved below the horizontal tail for flap positions greater than 10°.

NASA/TP--2000-209908 29 Thrust Effect on Tail Velocity Baseline; V~1.3Vs, c.g.=20%MAC

130 120 • dF=O

dF=10 i _ 100110 I = 90 o dF=20 [ > 80 dF=30

70 0 0 dF=40 60 0.00 0.05 0.10 0.15 0.20 0.25

CT

Figure 35. Thrust effect on V=a.

The effect of thrust on the tailplane section lift coefficient was of greater interest. For all the forward CG cases ex- amined, the tail lift coefficient increased as the thrust was increased (Figure 36). In the 5F = 0° and 10° cases, there were only small increases in the negative Ct_ra_Zrequired as the thrust increased. However, for the &F = 30 ° and 40 ° cases, there was a substantial increase in the Cl_r_t required as the thrust increased. This behavior was consistent with the requirement for greater tail down load to balance the nose-down pitching moment caused by the vertical dis- placement of the thrust vector to the CG position. Since the thrust line on the Twin Otter is above the CG, the in- creased thrust caused a nose-down pitching moment. With flaps deflected, the thrust was also deflected downwards, causing an even greater nose-down moment that required counter-balancing by the tail down load.

Thrust Effect on Tail Section Lift Coefficient Baseline; V~1.3Vs, c.g.--20%MAC

!e dF--O I l _, dF=l 0 I • dF=20

It dF--30

o dF=401

p q 0.05 0.10 0.15 0.20 0.25

CT

Figure 36. Thrust effect on CLrou.

NASAFFP--2000-209908 30 : ._ 5.3.9 Spanwise Variation of Flow at the Horizontal Tail

The variation of the flow field across the semi-span of the tailplane was examined through use of the three 5-hole probes. For each probe, the _it, fl,_, and Vt_l are plotted (Figure 37 through Figure 39) for 6[" = [0 °, 10°, 20 °, 30 °, 40 °] to display the complex flow field at the tailplane. The primary cause of the spanwise variation is likely the com- bination of the propeller slipstream and flap deflection. Secondary contributors to the variation are the fuselage and vertical tail, and other aircraft features like engine nacelles and wing may also play a role in the complex na- ture of the tailplane flow field. It should also be noted that the data presented here are for one set of conditions, and if those conditions are varied, the flow field can be altered significantly.

For a given flap setting, the _l typically varied 1°-2 ° across the semi-span (Figure 37). For all flap settings except 6/" = 0°, the outboard position measured the most negative _. With 6/' = 0°, the outboard position measured the most positive angle, indicating the presence of some flow disturbance ahead of this probe location for that flap posi- tion. Such a flow disturbance may be caused by the propeller wash, engine exhaust, nacelle, or wing . With the low level of variation in a_,_t, particularly with the larger flaps settings, the leading edge flow separation caused by ice should be fairly uniform across the span. From observation of yarn tufts during the flight, this was commonly found to be true during fl = 0° test points.

Span Effect on Flow Angle Va:c=1.3V,, I],/o=0,C.F--0.11

4.0 2.0 0.0 [+ dF=O -2.0 •---_.-.-dF=lO ,,,."° -4.0 _ ...... _,__-- -6.0 + dF=20 -8.o ---M, dF=30 --e-- dF=401 -12.0 _ o .e -- -14.0

,,D £ £ £ (3.. (3,. 0_ Tail Probe No

Figure 37. Span variation of o;,u.

Looking at the ]3t,it (Figure 38), the typical variation for a given flap setting was 5° across the semi-span, with z_[_tail = 7° for the 6/" = 10° case. Typically the trend was an outboard side-wash near the tail root, a small inboard side-wash at the midspan, and a larger inboard side-wash near the tail tip. The one exception to this trend was again the SF = 0" case, where the largest inboard side-wash was measured at the midspan position. The outboard side-wash near the root was most likely due to the flow around the vertical tail. With the flow probe on the left side of the tail, the flow circulating around the vertical tail caused a positive _ measurement on the tail probe. The inboard side-wash near the tip was most likely caused by the propeller wash or tip vortex flow. The change in side-wash direction near the mid- span location indicates that there was a spanwise location of confluence between the outboard wash caused by the vertical tail and the inboard wash caused by the propeller wash. At that location, there would be a small acceleration in the flow speed. The values of these _ angles were fairly small as expected since the horizontal tail had a 0° sweep angle, and would not play a significant factor in the ice-contaminated tail stall event.

NASA/TP--2000-209908 31 Span Effect on Flow Angle V_c -1.3V,, [3_=0, CT=0.11 6.0 ,0 t + dF--O 2"0 -..t, ... dF=lO --.6,.----.dF=20 - -_---- dF=30 -2.0 --.-e..- dF=40

-4.0 A" v -6.0[

2 _o o

Tail Probe No

Figure 38. Span variation of/_,_m The most significant variation in the flow field was in the velocities. The propeller wash effects are clearly observed in the Vt,,itplots (Figure 39)• With t_F= 0, the mid-span probe measured a significant increase (20 knots) in the flow velocity (V_/c= 92 knots). As the flap angle increased to 10°, the propeller wash was more concentrated at the out- board position (probe 1), where a 20 knot increase was observed above the aircraft airspeed• As flaps were extended between 20° to 40°, there was a small velocity increase above the aircraft airspeed observed at the midspan probe, but the effect was much smaller• These observations further support the idea that the propeller slipstream moves be- low the tailplane as flaps extend beyond 10°. Also, the atypical behavior in c_,il and fltaitfor t_F= 0° and 10° is better established as a propeller wash effect. Finally, it should be noted that for the &F = 0° and 10° cases, the dynamic pressure at the tail (q,,,it)can be higher than the dynamic pressure of the aircraft (qa/_).

Span Effect on Flow Speed V,zc ---1.3V,, J],vc--O,CT=0.11

120.0

110.0 -- _dF=O m 100.0 --._-.---dF=l 0 "-" 90.0 dF=20 ...... ,m- -dF--30 > 80.0 --,e-- d F=40 70.0 60.0 O,J 03

2 2 2 Q. 0. Q_ Tail Probe No

Figure 39. Span variation of Vt,m As statedabove, the flow field at the tail can be very complex with mu|tip]e causeand effect relationships.The variation across the tail span was demonstrated in this section and possible explanations were provided. For the cases examined, the variation does not appear to be significant enough to cause a spanwise variation of ice-induced sepa- ration.

NASA/TP--2000-209908 32 5.3.10 Effect of Steady Heading Sideslips on Tailplane Flow Field and Performance

Steady heading sideslip (SHSS) maneuvers were included as part of the test matrix since an SHSS may help to de- termine susceptibility to ice-contaminated tailplane stall for some aircraft _3.The SHSS maneuvers conducted as part of this investigation revealed more about the complex flow field presented to the horizontal tailplane. Figure 40 shows the effects of aircraft sideslip on tailplane angle of attack measured by the mid-span flow probe. Note that these data are for the no-ice baseline and all iced tail configurations with the test point velocity approximately at 1.3Vs speed. As noted in the previous section, there is a strong influence of the propeller wash for the _F = 0° case at the mid-span probe. However, as the flaps deflected beyond 10°, the trends due to aircraft sideslip were evident. The positive SHSS caused the tail angle of attack to become more positive (approximately +4°). In comparison, the negative SHSS caused only a small positive increase in the tail angle of attack. Similar trends were observed for the probe near the tail root (probe 3). The outboard probe (probe l) generally had similar trends, but the positive in- crease with negative SHSS was comparable with the increase due to positive SHSS (approximately +4°).

Sideslip Effects on tan All Configurations, V-l.3V_ ,o...... 5 dF=0

A dF=10

• dF=-20

_-10 U dF=30

-15 <>dF=40

-20 -20 -15 -10 -5 0 5 10 15 20

(deg)

Figure 40. Sideslip effect on c_oa. The significance of the positive increase in tail angle of attack with sideslip is that the tail stall margin is reduced with the SHSS maneuver. The cause of this positive increase in a,_i;was most likely due to the increase in aircraft angle of attack (Figure 41) to maintain constant lift for holding altitude or vertical speed. In performing the SHSS maneuver, a bank angle and a crab angle were held relatively constant. With the wings banked (up to 14° bank an- gle), the lift component normal to the aircraft was increased by the increase in angle of attack to maintain a relatively constant vertical lift component. This increase in aircraft angle of attack effectively made the tail angle of attack less negative and reduced the tail stall margin.

NASA/TP--2000-209908 33 Sideslip Effects on cx a/c

All Configurations, V~1.3V_

• dF=-0

_ dF=l 0

o dF=20

I_l_iidF=30

<>dF--40

-20 -15 -10 -5 0 5_ 10 15 20

13 (deg)

Figure 41. Sideslip effect on ¢Y_/c.

Observations of the yarn tufts during the SHSS maneuvers confirmed this finding. At the larger flap settings with the Failed Boot or S&C ice shapes, a significant chordwise separation bubble existed across the entire span of the tail- plane with the aircraft at ]3 = 0°. As the aircraft was sideslipped in either direction, the size of the separation bubble was reduced and moved away from the direction of the sideslip. Also, there was significant pulsing in the elevator caused by flow separation and reattachment for the _ = 0° cases described. This pulsing was effectively eliminated when the aircraft was placed in a sideslip angle greater than approximately 12°.

The results from the pressure belt data for the SHSS maneuvers are shown in Figure 42. The trends suggest that less down lift is needed with the positive sideslip and more down lift with the negative sideslip. However, this is an arti- fact of the test setup and the flow field presented to the pressure belt. The primary cause of these trends was the flow around the vertical stabilizer and rudder. The positive sideslips were accomplished with a rudder deflection to the left side horizontal tall - where the pressure belt was located. The leading edge of the left side of the vertical tail had a large suction peak with a pressure rise over the left side of the rudder. The flow field around the vertical stabilizer and rudder affected the pressures near the pressure belt such that the integrated values showed a more positive lift. On th e other hand, when the rudder was moved to the right side for the negative sideslips, the flow field_arou_ nd the vertical/rudder combination caused the integrated pressures about the horizontal tail to become more negative than the zero sideslip case. Since this behavior is symmetrical about the zero sideslip case, it is reasonable toassume that the increased section lift coefficient on one side is equally decreased on the other side during the sideslip maneuver to accomplish a constant tail lift during the SHSS.

NASA/TP--2000-209908 34 Sideslip Effect on Tail Section Lift Coefficient All Configurations, V~1.3Vs, C,T=0.10

o.o t ...... • dF=O dF=l 0 -0.2 • dF=20 O -0.4 El dF=30

-0.6 0 dF=40

-0.8 -20 -15 -10 -5 0 5 10 15 20

13(deg)

Figure 42. Sideslip effect on Cl_Taa.

The SHSS maneuver was performed for wide range of test conditions and ice configurations on the NASA Twin Otter. None of the SHSS maneuvers led to a stalled condition. To the contrary, the SHSS maneuvers significantly reduced flow separation caused by leading edge contamination. It should be noted no ice contamination was tested on the vertical stabilizer during these tests. If ice contamination had been present on the vertical stabilizer, it would have affected the pressure field about the vertical and horizontal tail junction. It is unknown whether such a distur- bance would have led to a greater sensitivity to sideslip. In addition, since there were some aircraft known to have a sensitivity to tail stall preceded by sideslip, the results shown above should not be considered universal. The reduc- tion of flow separation experienced in these tests may have been a configuration specific phenomenon, and as such, the SHSS is not a recommended procedure for recovering from an impending tail stall.

NASA/TP--2000-209908 35 5.4 Steady-State Analysis Conclusions

Steady state data points with the test aircraft's tail configured with various levels of ice contamination were acquired and analyzed as described in the previous sections, For steady state flight, the following conclusions were made:

1. For any ice configuration, the most negative tailplane AOA and smallest tailplane stall margin occur at full flap deflection while flying at flap extension speed, Vre. 2. The difference in elevator deflection angle between a clean and iced tail is a good indicator of reduced tailplane performance. The larger the difference, the greater the performance degradation. This assumes that aircraft con- figuration (flaps, , CG, etc) and flight condition (airspeed, thrust setting) are held constant for the clean and iced cases considered. This method clearly discriminated between the iced and clean cases, but it did not discriminate between the Failed Boot ice and the Inter-Cycle ice. 3. The difference in elevator hinge moment is also a good indicator of reduced tailplane performance. It may be better than the elevator deflection method since it can discriminate between smaller levels of degradation as in the Failed Boot and Inter-Cycle ice cases. The elevator hinge moment method requires consistency in aircraft configuration and flight condition between the baseline and iced cases, and requires the pilot to fly the aircraft in an untrimmed condition to acquire the appropriate forces. This method also provided data on the progressive hinge moment lightening and reversal as the tailplane stall margin was reduced. 4. As the tailplane approached stall, the trim tab effectiveness was reduced such that the tab could not effectively balance the elevator hinge moment. This provides another good indicator of approaching ice-contaminated tail- plane stall. In other words, the inability or difficulty to trim the iced aircraft with flaps deflected could indicate an imminent tail stall. 5. A forward CG location demanded a higher performance from the horizontal tail to provide stability. It was shown that the most negative tail section lift coefficient was required with a forward CG, the flaps fully de- flected, and at a speed near stall speed, V = 1.0Vs. 6. As thrust increased, a more negative tail section lift coefficient was required of the horizontal tail to maintain stability. This was most likely due to the vertical displacement between the CG and the high thrust line of the test aircraft causing a nose-down pitching moment due to thrust increase. It is reasonable to assume that high thrust line aircraft would have similar results. Thrust did not have a substantial effect on c_,_ or Vt,_t when flaps were deflected beyond _SF= 10 °. 7. The flow field presented to the tail was quite complex especially with the flaps retracted or at 8F = 10°. The complex nature at these flap settings was most likely due to the propeller wash impinging on the horizontal tail. However, at these flap settings, the tail stall margin was large, so the complexity in flow field did not provide in- sight into the tail stall phenomenon for the test aircraft.

8. For flaps greater than _SF = 10°, there was relatively little spanwise variation in c_,_, _tail or Vtail SO that this tail would stall relatively symmetrically across the span. This result was verified through use of the yarn tuft visual- izing the separation bubble on the lower surface of the tail. 9. Steady heading sideslips did not precipitate an ICTS condition for the test aircraft. To the contrary, SHSS alle- viated flow separation caused by leading edge ice contamination and eliminated control pulsing when the side- slip angle reached approximately t3 = ±12 °. The reduction in flow separation is likely due to an increase in air- craft AOA and subsequent positive increase in tail AOA, causing an increase in tail stall margin. However, since no ice contamination was present on the vertical tail, and the interaction of flow separation on the vertical may play a role in flow separation on the horizontal tail, this conclusion is limited to the clean vertical tail cases.

NASA/TP--2000-209908 36 6.0 Dynamic Maneuver Analyses

The primary purpose of the dynamic maneuvers was to measure the aerodynamic behavior of the horizontal tailplane while driving it toward stall. This was achieved by either intentionally manipulating the elevator or changing one of several parameters. For convenience, the parameter space is reiterated below: • ice shape severity • flap deflection, tSF • speed, VIAS (aircraft angle of attack, trwc) • thrust, Cr Increasing any one of these parameters will put the aircraft on a path toward tail stall. The "paths to stall" maneuvers flown include the flap transition, speed transition and thrust transition. Another category of maneuver involves the controlled manipulation of the elevator. These include the elevator doublet and pushover maneuvers. Finally, a closed-loop tracking task was performed: the balked landing maneuver, which simulated an ILS approach and go- around. In contrast to the other maneuvers, the closed-loop task offered a clear criterion for whether or not the pilot met the performance objective. Table 1 lists these dynamic maneuvers, along with the number of test points flown and the parameter space. It should be noted that the location of the aircraft's center of gravity was also briefly stud- ied. As expected, the more forward the CG, the greater the loading on the tailplane. This is discussed in Section 5.3.7. The other parameter mentioned in the steady state section, the sideslip angle, was not an element of the dy- namic maneuvers.

Maneuver Test Points Comments Data Compatibility 22 All Shapes, tSF = 0 °, V = 1.8Vs, Thrust for level flight Flap Transition 6 Clean Tail and Failed Boot, Moderate Thrust Speed Transition 2 Failed Boot, Full Flaps, Moderate Thrust Thrust Transition 98 All Shapes, Flaps & Speeds, Mod Thrust Elevator Doublet 103 All Shapes, Flaps & Speeds, Mod & Hi Thrust Pushover 133 All Shapes, Flaps & Speeds, Mod & Hi Thrust Balked Landing Failed Boot, Flaps, V= 1.3Vs, Thrust as Required

Table 1. Dynamic maneuvers flown and number of test points for each.

Airfoils stall due to either a leading or trailing edge separation. For this series of flight tests, evidence of stall was documented by both the tailplane tuft video and the pressure distributions from the static pressure belt around the tailplane. Preliminary observations suggest agreement with the Ohio State University wind tunnel tests 9. That is, the clean leading edge was found to undergo only trailing edge separation at sufficiently negative tail angles of attack and elevator deflections (more cambered). For the ice shapes tested, on the other hand, separation occurred from the leading edge. From the tuft video, three distinct zones of activity along the underside of the tail could be identified. There were several cases where both the leading and trailing edge zones exhibited flow unsteadiness while the tufts in the middle zone were less active.

The dynamic maneuvers are introduced in Section 6.1. Herein, the pushover is also examined with some advanced analysis techniques. In Section 6.2, a direct comparison between the pushover and elevator doublet is presented. This section also contains a comparison between two relevant metrics regarding the pass/fall criteria: the control force reversal versus the undamped pitch response. In Section 6.3, the precision and accuracy with which one pilot flew the pushovers is investigated. The main findings are summarized in Section 6.4. Appendix C contains time histories of 21 quantities from select test points within the parameter space. Finally, Appendix D presents a time history of Cp distributions around the tailplane for one of the thrust transition test points.

NASA/TP--2000-209908 37 6.1 Maneuver Description and Analysis

6.1.1 Data Compatibility This maneuver was the first test point of each flight. The purpose was to check for instrument errors by performing this data consistency check. Beyond noting that the data are within range, further analysis had not been done. Description: As seen from the time histories in Figure 43, the maneuver began with a sharp, stick back input (about 1" column movement) so that the short period was excited. The column was held there for a short time (3 sec). Then, a sharp, forward stick movement (about 2") was applied and the column frozen there for a slightly longer time (5 see). Next, the column was sharply returned to trim and frozen. As soon as the elevator reached trim, sharp rudder step inputs on the order of 1-2" pedal deflection were applied and timed with the lateral response of the aircraft. The control column was kept frozen and the ailerons in the neutral position. Six alternating rudder kicks were applied in all, After the last kick, the pedals were returned to neutral and all controls were held frozen for a few more seconds. The elevator step inputs mark the pitch response and damping characteristics (q) of the Twin Otter, while the rudder kicks highlight the dutch-roll and damping characteristics (r). Configuration: This maneuver was flown with all ice shapes at 6F = 0 °, V= 1 .SVs = 120kts, Cr as required for level flight, and propeller RPM = 90%.

6 , , - 8

...... __ i- ,- 4 -, ...... delR 4 (9 "1o v iii "6 ...... : _ : -4 -8 t :..... "" , .B

A

"O 1o 0 v , _" i _ ?_,d. _ I _alpha -3 ...... \s -V _d -v ..... ___..._,. "_ JD -6 '_ ' " -I 0

2o -- 275 A 10 ...... _- .., - ...... _- 265 1D v i' / i /'";" ...... , - t "1o J_ 0 ...... "F"'"" I 255

o. -10 a,i ...... '--- _ ...... : I--Ph I -i-24_ r- -20 : .... ; I-:- P'_ i 2a8 10 7 3O A I __ ..... :...... Y-" - r,-: ..... :- 15 m == 0 "o os, b'=''_J --; " =-11i ...... !'...... _'_ ' _ "' ' v 0. .o4...... i" :t=! -15 ,Z i V : : ,[ --P -10 J -30

0 5 10 15 20 25 30 time (e) 35

Figure 43. Data compatibility time traces.

NASA/TP--2000-209908 38 6.1.2 Flap Transition One of the key parameters, of course, was the flap deflection angle. All of the test points flown included this setting as a factor. To more efficiently isolate the effect of the flap deflection, the flap transition was constructed and flown six times. Description: The flaps were deployed from zero to full deflection, 00-37.5 ° (hereinafter rounded to 40 °) while the remainder of the parameter space remained unchanged. Configuration: This maneuver was flown for only the Baseline and Failed Boot cases at a moderate thrust setting, Cr = 0.10, and a constant airspeed, 85 KIAS, which corresponded to 1.3-1.6Vs as tSF increased steadily from 0° to 40 °. Analysis: The flap transition profile for the aircraft is depicted in the sketch to the right. Selected time histories are shown in Figure 44 for the __ t_F = 40 Baseline case. As the flap angle increased, the aircraft angle of attack was C La/c decreased to maintain velocity. Since neither the speed nor the aircraft lift (= weight) changed, CL,vc was also constant. Furthermore, a constant CL,V_ also implies that the downwash is constant (see Section 5.3.2). The de- crease in u.,vc drove a_ more negative. The more negative c_ caused the / / Ct_rait to be more negative than required; thus the SE was increased (de- flected TED) to maintain the speed condition. Had _E been maintained, V c_ the aircraft would have pitched up and slowed down.

The nominally same flight conditions with the Failed Boot ice shape are presented in Figure 45. Also plotted on Figure 45 are the equivalent trim cases (o, *). These points were added because the speed was not as tightly con- trolled during this test point (the flight was for demonstration purposes). In general, the trends are very similar to those of the Baseline. Some key differences between the clean and iced tail, however, appear in the traces of elevator deflection (delE) and control force (FYE). For the Baseline case, a push force was required throughout the entire maneuver, and the elevator deflection changed only slightly for tSF > 20 °. For the Failed Boot case, on the other hand, the required force transitioned from a push to a pull force at SF = 20 °. Beyond 6F > 35 °, oscillations in the elevator appeared both in the data and are dramatically apparent from a video shot by a chase plane. The video also captured that the tufts across the whole underside of the tailplane had lifted off the surface and were moving fore to aft.

While maintaining a high airspeed, increasing the flap deflection with an ice contamination drove the tailplane to- ward stall by driving o_ more negative.

NASA/TP--2000-209908 39 40

30

'10 2O LL 10

0

10 I ...... 60

= 6 ...... _ - - E(d_ - 30 %-

,,, (lbs) ' LLI >- of ...... _ .... _ , _ , , , _ , ',o LL

-5 _" "_'_" .... - ...... - -30

47 - 0 0 __.. , , , , f--TAOAI(deo) -- -0.2 .o-_'!iiiiii_i iii--°-°_ -16 " -1

8___

_ -4 ......

< -8 ..... i -12

1,5 7 CLac _- 10 _dwnwsh (deg) _ -_ I _ -----8

0.5 ..... 6 • c

7600 _ 115 I --ALT (ft) }

7300"°°i_...... 12...... _ ...... 85°'_

7200 I I t I t ! I I 75 0 5 I0 15 20 25 30 35 time 40

Figure 44. Flap transition; clean tail, VIAS = 85kts.

NASA/TP--2000-209908 40 40- , , ©

30 - " ',...... i" - - ---_ - _detF(deo) _b "13 LL .O ioo iiiiiiiiiilli, .!e,..

10 - , , _ 60

W W

"" _', _" _ I © delEtrim / u. - : '_,%,_,,_K'''_' : [--..-L-.-FYE(Ibs) I -5 - -30

4 - ...... ,...... -...... --TAOAI (deg) I 0 A 0 _ ...... i""'_,_, i I I O TAOA1 trim _-0.2 m _ ...... %.,--'-;, ----CLIail /

< -12 - -_ -0.8 p- -16 - -1

131 @

=-

-12

1.5 - • , -- _ 12

10 , O : C'.'. .O. - _ - ..... I 8 - ;_ ', _,,,._A.,_,,,_.-_-_'-_._--_.,_.AL._,,,_.,._ / 0.5 - , _._ .... - ...... _CLac _ 4 ,_. _ , _.._..._,_._ . 0 CLac trim I ...... _'; ...... dw nwsh (deg) 0 t i ! i I t X_ dwnwshlrim .-_ 0

7800 .... 125

7700 ..... " ...... ' ._._ I _A L_ t 115 "_ 7600 - ' - - _,_"_'-<,.- ...... I ..... VTAS (kts) I . --_ 105

7500 ' • ' • ....._.i,- -_ _ • ...... - / 95 > 7400 I I I I I t I I | 85 0 5 10 15 20 25 30 35 tim • 40

Figure 45. Flap transition; Failed Boot, VIAS = 80kts (VTAS = 9lkts).

NASA/TP--2000-209908 41 6.1.3 Speed Transition The speed transition was a singular test point until the Guest Pilot Work- 4O shop. The original intent for this test point was to fly steady wings level 1 Ct_/c with the Failed Boot ice shape at full flaps, moderate thrust CT = 0.10, and V = 85kts = 1.5Vs, the highest target speed. The experience from the pre- vious test point at V = 75kts had indicated this might be difficult. The pilot / started from a lower speed (65kts = 1.2Vs) and gradually brought the nose down to increase to the target speed. The graph to the right depicts the Ii, transition of aircraft lift coefficient during for this maneuver. tXa/c

The data record was initiated at 73kts = 1.4Vs because the tufts were indicating highly unsteady flow, see Figure 46. As the aircraft angle of attack decreased from -5.5 ° to -7.5 ° and the elevator deflected TED (gross trend from 2° to 4°), the control yoke started to oscillate back and forth. As noted from the tuft video, this was associated with flow separation and reattachment beyond the hinge point. With increasing speed, the oscillation of the control forces grew rapidly. By the target speed, the force cycled between 0-50 lbf at a nominally 0.3 Hz cycle; the traces also show smaller oscillations on the order of 2 Hz. Note that the mean _ line remained unchanged throughout the entire ma- neuver. It is interesting to note that the c_ features lead the aojc features by about 1.4 sec. This suggests that the events at the tailplane are dominating the flight mechanics.

The events at the tailplane presented in Figure 47 demonstrate that indeed, the tailplane was going in and out of stall. Following the discussion outlined Section 5.3.3, two data points representing the extrema of one of these cycles were selected. These are depicted as circles in Figure 47(a); they are the local maximum at t = 26.07s, (c_ = -11.4 °, tSE = 0.5 °) and the local minimum at t = 28.41s (_ = -13.3 °, tSE = 6.3°). In Figure 47(b), these two points from the speed transition are graphed against the Cl_rail vs. _ plot from the wind tunnel. The correction to the wind tunnel data re- fers to the standard 3-D correction,

NASA/TP--2000-209908 42 9O _VIAS (kts)! 85 l

8O

75

7O I , I I I

8 I 6 • ' ' ' ' _delE(deg)

4

2

0 I

-2

-11 (deg) I

-12 . ', , TA OA 1

-13

-14

-4 • --Alpha (deg) I -5 _ '/_ ...... ; /_ ...... dw nw sh (deg)_l_ -6

-7

-8

-9 p _ t 4

4

, • - T , .J 2

0

-2

-4

20 - - -

0 • ..: i-- Eibs,l

-20

0.9

0.8 E I i 0 5 10 15 20 25 30 35 40 45 time 50

Figure 46. Speed transition time histories; Failed Boot at _F = 20° and Cr = 0,1.

NASA/TP--2000-209908 43 -11

O -13

2 J

-14 0 5 10 15 20 25 30 35 40 45 50 tim e (s) (.) Wind Tunnel Comparison

0 , , . . . ,_ , .

_1 -0.4 _-%_"_/+_/_-'--.- _-+,_-_--"_-___+Raw Wind Tunneh deIE-O O _'_ ,,,, 1 _ _ ..... +----Raw Wind Tunnel: delE=10 i _,_ _" .. . -0.6 _ _ .- .r f "" _ - + Speed Transition Flight Data ..__ _'_---._ _ _.ff" _ _ Corrected Wind TunneldelE=0

i 2 _---- _'_._ __ _ Corrected W,nd Tunnel delE=lO

-16 -14 -12 -10 -8 -6 -4 -2 0

c_ Tail

(b)

Figure 47. Detail of speed transition: (a) o_ time history with two highlighted data points, and (b) Ct_r, it vs. o_.

NASA/TP--2000-209908 44 6.1.4 Thrust Transition The thrust transition was added in mid-program after the Twin Otter's sensitivity to thrust effects became apparent. Descriotion: The thrust was increased from idle to maximum continuous while the remainder of the parameter space was maintained. To maintain airspeed while increasing thrust, the pilot had to pitch the nose up. The thrust line for the Twin Otter is above the CG. Configuration: This maneuver was flown with all ice shapes, flap deflections and speeds; except that the S&C shape was not flown beyond tSF = 30 °. Analysis: Because the thrust line is positioned above the CG, this maneuver challenged the horizontal stabilizer by increasing the download requirement at the tail. This maneuver also drives the elevator toward the common landing position - TEU. While the intent of the program was only to fly toward the stall angle, a full tail stall did occur dur- ing one execution of this maneuver. This stall event occurred for the Failed Boot ice shape at full flap deflection and maximum speed, 85kts = 1.SVs. Figure 48 presents selected time histories as well as the probable stall point and the time at which thrust was reduced. Corrective actions - thrust to idle, yoke back and flaps up - were applied within 0.4 sec of stall, yet the aircraft lost a total of 300 ft in 8 sec and reached a (nose down) pitch angle of -37 °. The "kink" in both the a and o_ curves at 32.09 sec was identified to be the probable stall point. From the tuft video, there was a point where the last row of tufts (x/c = 0.95) all pointed upstream. Prior to stall, the _ signal led the tr._c signal by about 0.7 sec.

During the recovery effort, note that even with considerable and increasing pull force, 75-100 lbs, the elevator con- tinued to move TED. With the thrust reduced to idle and the flaps rising, the pilot was able to continue to pull, 100- 172 lbs, to bring the elevator TEU. Eventually, according to the video, the flow reattached around 35 sec and the pitch over ended. Note that even though the pilot had regained control of the aircraft, the greatest altitude loss (200 ft) occurred between t = 35-37 sec because of the increased airspeed.

Select static pressure distributions are shown in Appendix D for various times. From these figures, the pressure dis- tribution on the lower surface flat lined - indicating full separation - after 32.34 sec. Recall from Section 3.2 that the system response lag was estimated to be about 0.2 sec; this corresponds well with the stall point as identified by the kink in the ct and at traces. The distribution remained fiat until after 34.86 sec.

NASA/TP--2000-209908 45 115 0,15 . --VIAS 105 95 85 r..9 < 75 > 65 55 ' _,{ Powe_ Reduce_] 0 12 .,.

UJ 0 i

"o -4 - - -

-8

._. -8 -12

-16 --

-4 i st'"Po'n'!

0 -t2 < I.- -16 ,: 200 , :

.a v 100 l.IJ

u_ 50

0

2.5 ': ,;

" o.s

0 -0.5 -- -._

10 ': S 0 ' i i ...... i } i : ] ' " _ ' _ -_o ..... ,..... , -' .... : .... ' -:i--'- '......

-20 ° -30 iiiiiiiiii' - ' '

-40 -- ,o 7400 - ,: 45 ,,._..___ , ! ' _ t--ALT

¢D

,_ 88oo ...... ,_ ..... 7a>.L I ...... 3o U- 6600 ...... ------ii _._- _ .... 25 1D 8400 T .... ! , , 20 0 5 10 15 20 25 30 35 40 45 tim • 50

Figure 48. Thrust transition time histories; Failed Boot at 8F = 40 ° and V = 1.5Vs.

NASAJTP--2000-209908 46 Propwash Effect The thrust transition maneuver is also an ideal method to investigate the effect of propwash on the horizontal stabi- lizer. The assumption is that if or when the flow from the propeller impacts the tailplane, increasing the thrust will identify the propwash impact. The output of the three 5-hole probes on the tailplane was compared to the output of the noseboom and the coefficient of thrust profile.

These comparisons were made primarily for the clean tail; an ice shape on the leading edge of the tailplane should have minimal effect on the propwash effect. However, a few of the observed trends were also checked against some of the ice shape results. The results show the same trend.

Figure A-1 details the aircraft probe locations with respect to the propeller. The outboard probe, probe #1, is aligned directly behind the hub of the propeller. The midspan probe, probe #2, is aligned behind the propeller tip, and the inboard probe, probe #3 is located upstream of the vortex generators on the horizontal stabilizer. The relative loca- tion between the horizontal stabilizer and wing is also identified.

There is a propwash effect in the span; it is a function of flap setting and airspeed. This effect was determined by plotting the three tail probe velocities in time along with the aircraft velocity and also noting the thrust increase. For the Twin Otter, the greatest effect occurred for t_F = 10°. Figure 49 presents a limited set of these results: the clean tail at different flap settings for nominally 1.2Vs. This set is representative of all the speeds flown with clean tail. This figure demonstrates that for tSF = 10°, the true airspeeds from the noseboom (VTAS =- Va/c) and two inner probes (Vprobe2 & 3) are flat and parallel to each other. The outboard true airspeed (Vprobel), however, climbs like the thrust (Cr) profile. The calculations indicate a 30 knot increase in the outboard velocity over that of the nose- boom: Vprobel = 120kts, Va/c = 90kts. Note that the probes were calibrated to 120kts, and values were extrapolated above 100kts. Even though the numerical values should not be accepted with certainty, the trend is indisputable. Also, there is an acknowledged difference between VTAS and Vprobe values beyond _F = 20 °. The Vprobe values are believed to be nominally 10% greater than actual. This was found to be a result of calibration methodology.

Now the issue concerning the mechanism that drove the Twin Otter tailplane toward stall (see Figure 48) at large flap deflections and moderately high thrust settings can be resolved. The two competing hypotheses were these: (1) turbulence from the propeller impacted the horizontal tail, destabilized the boundary layer and caused separa- tion, or (2) because the thrust line is above the CG, thrust addition increases the nose down pitching moment. To maintain airspeed, the tailplane had to provide more downward lift, i.e., increase it's camber by deflecting the elevator TEU. Knowing that the propwash at t_F = 40 ° slips well below the tail eliminates the first explanation. The leading hy- pothesis, therefore, is the high thrust line. Ultimately, this should be verified with other aircraft. Until this issue is resolved fully, pilots of aircraft with thrust lines above the CG should be aware that adding thrust in a potential tail stall situation might only exacerbate the situation.

NASA/TP--2000-209908 47 120 _VTAS (kts) Ba,a,_' = o,vms_'== z.2o ,_,_"------_,%_,%,_., ...... Vpr obel (kts) ...... Vprobe2 (kts) Vpr obe3 (k_) 2

70 ...... _.., I

130 _VTAS (kts) Ba, d.F = ZO,V_Kg/V's = 132 .... Vprobel (kts) ]20 ...... V_obe2 (kls) --.. V pr obe3 (kts) 110

]00 _J 2 i

90

8O

II0 ...... _VTAS (k_) Ba,41F = 2O,V3_._N's = :I.,24 ...... VproDel I00 V_o_3 (kt_) 90

80 3 70

60

]10 B a,d_' - 40,Vl_gN., = 1.25 qO0 _VTAS (kts) -- V pr ot::,eq (k_) 90 ...... Vprobe2 (k_) Vproloe3 (kts) 8O ___ .... - _-: _ _i__'_ 1,3

70

6O | 0 I0 20 30 ,40 50 60

0.3

0.2 _ _ - "_ "_ _*_-'_--'--'_

J--'- , ' __-- ' -- 40 0 k 1 f 0 I0 20 30 40 50 60 Time (sec)

Figure 49. Propwash effect for a clean tail at different flap settings and constant speed, V = 1.25Vs.

NASA/TP--2000-209908 48 6.1.5 Pushover Description: The pushover was flown in a manner similar to that described in current aircraft certification programs. The aircraft was configured at the proper thrust and flap deflection, then trimmed at the target airspeed, V,arge, The pitch angle, Otarge,,was noted. The task was to hit a specified pitch rate, qtarget, (or G) at Vtarget and as the nose of the aircraft tracked through the horizon, Otarse,.(The pitch rate and vertical acceleration are coupled.) In order to achieve this from trim, the pilot first pushed forward on the column to gain speed, then pulled back on the column to raise the pitch angle prior to the pushover at V_get. Just prior to Vtarse, the control column was pushed forward to hit qtarget at 0_se,. One such maneuver is called a parabola. For the Twin Otter, the parabolas typically lasted up to 20 sec with a maximum of five seconds at the minimum gravity. For this flight program, three parabolas were flown in a row dur- ing one test point. Configuration: This maneuver was flown for all ice shapes and the full range of flap deflections and speeds with Cr = 0.10. However, because a control force reversal (CFR) was experienced at tSF = 20 ° for both the Failed Boot and S&C ice shapes, no pushovers were flown beyond that flap setting. Anal sy__: The two extrema of a "pass" and a "fail" case will be discussed here. More examples can be found in Ap- pendix C, which contains time histories of eight test points. In addition, and an example of a "shaped" maneuver and a trailing edge stall on a clean tailplane will be presented. Finally, the more advanced analysis techniques of co- and cross-plots will be introduced.

The pushover maneuver was designed to achieve high angles of attack at the tailplane dynamically. The tailplane angle of attack is given by Ref. 6.

a,=a,o+Aa,=[aw-e-i,]+ Aot_,, 1- + , where a_ois the trim (steady state) angle of attack and Aft is the dynamic component. The change in wing angle of attack is denoted Ate, and It is the distance from the CG to the tail center of lift. The (v3d_ga) values vary from 0.4- 0.7 as a function of SF for the Twin Otter. Therefore, the dynamic tailplane angle of attack may be increased primarily by increasing the pitch rate or reducing the speed. Moreover, high flap settings and high speed increase C_o.

For ease of discussion, it is prudent to define some terminology. Let a pushover maneuver where the values of o_ are far from O_:aU, be considered a "non-critical" pushover. Likewise, denote a case where the values of _ are close to __statt as "critical". For a target pitch rate, the most non-critical pushover therefore occurs at zero flap deflection and high speed. Conversely, the most critical is with full flap deflection and at a low speed.

Figure 50 illustrates a limited set of the time histories for the most non-critical case with no ice. Note that this was the very first set of pushover maneuvers during this flight test series. The three pushovers were flown to increasing lower G; specifically, the task was to increase the target q: qt_se, = -5 deg/s, -10 deg/s and required q to achieve 0G. The pitch rate mimicked the vertical acceleration, albeit not perfecdy. The target speed was VIAS = 1.5Vs = 100kts. In all cases, the pilot accurately hit qt_sa (or Nz) at the Or_rg_= -1 ° but under flew the speed." The most negative o_ was -6 °, far from stall, and the most negative Ct_r_i_ was -0.3.

Note that in Figure 50, both SE & FYE and tSE & q are co-plotted on the same graph. Co-plots of certain parameters are perhaps the best and most immediate way to interpret dynamic data. For example, a co-plot of tSE & FYE imme- diately indicates whether or not a CFR has occurred. There are two inputs to the force measured at the yoke: (1) the pilot input and (2) the pressure-field around the elevator. If the elevator is fixed TED, yet the control force crosses the trim point, then it must be due solely to changes in the pressure field, and a CFR has been experienced. For this non-critical case, the push force remained fairly flat during the maneuver. Similarly, from a t_E & q co-plot one can see whether the pitch rate response is damped or undamped. In this case, the response is damped; the most negative q peaked shortly after the elevator is deflected and then self-restores.

"In subsequent test points, the information to the pilot was improved, and the target velocity was better met.

NASA/TP--2000-209908 49 AcriticalpushoverisdepictedinFigure51.WiththeFailedBooticeshapeatt_F = 20 °, the target was to hit N z = 0G at V = 55kts = 1.06Vs. While the pilot hit V,a,ge_at 0t_rgez--=--5° fairly well, the minimum N z occurred only for a much steeper pitch angle, 0 = -30 °. For this configuration, the most negative a; was -14 °. The most negative C_r_il was -0.9; the instability evidenced between t = 51-53 sec suggests stall. This is confirmed by the tuft video. Also observe that the yoke force (FYE) started to lighten as soon as the elevator was pushed TED. Indeed, during this maneuver, full control force reversals (CFR) were experienced. In fact, the control force crossed neutral at t = 17.6, 34.0 and 47.5 sec while the elevator was fixed forward (tSE = 10°). The tSE & q co-plot also clearly indicated an undamped response; the decrease in pitch rate only ceased when the elevator was returned TEU.

NASA/TP--2000-209908 50 2.5 _Nz (G)

1.5 1 0,5 0 -0.5

160

_VIAS (kts) I

120

100

80 ,

40 , _ . ! _ theta (deg) I

-20

-4O

10 --delE (deg) FY E (Ibs) j _ 80

-5'-...... /< i I '°40 -10 -80

S i 4 ...... '...... _ --'raoA1 (deg)

0 -

-4

-8

0 i _ CLtal] -0,1 .... ' -- J_"

-0,3 ......

-0.4

8

(deg) I

-4 ......

-8

0 10 20 30 40 50 60 70 tim • 80

Figure 50. Non.critical pushover time histories: clean tail, 8F = 0° and V = 1.5Vs = 100kts.

NASA/TP--2000-209908 51 Nz (g) 1.25 1 0.75 0.5 0.25 0 ! I

95 , - -- (kts)

75

65

55

45

30 ..... ',...... theta 20 [deg] 10 0 -10

-20 -30

14 ...... delE (deg) FY E (Ibs) ] ...... 50

7

0

-7

-14

20

10

0

-10

-20

0 1 I I I I I 1 . __

-4

-8

-12

-16

0

-0.2

-0.4

-0.6

-0.8

-1

8

4

0 a = I

-4 -8 I v -12 0 10 20 30 40 50 60 time (,,) 70

Figure 51. Critical pushover time histories: Failed Boot ice shape, ¢_F = 20 ° and V = 1.06Vs = 55kts.

NASA/TP--2000-209908 52 Shaped Maneuver Figure 52 illustrates how a maneuver may be "shaped" to pass the no CFR criteria. At 18 sec, an adjustment of the elevator is sufficient to reattach the flow and cause the force to remain on the push side. The other two maneuvers show definite lightening and indicate a CFR.

--ddE 90

!\ ...... FYE 60

U.I " -30 _ -60

0 10 20 30 40 50 time (s) 60

Figure 52. Example of elevator shaping during pushover maneuver.

Trailing Edge Stall The pushover maneuver has an additional challenge. During the pull up, the elevator deflects TEU suddenly. Cam- bering the tail increases the download capacity, but toward a reduced stall margin (see inset figure). If the tailplane were close to the stall margin, the rapid cambering might push a_ through the stall point. Evidence of this is pre- sented from the 1995 flight campaign. An extreme Baseline case of full flaps, tSF = 40 °, low speed, VIAS = 50kts = 0.94Vs and TED, 8E=const I high thrust, Cr = 0.23, is depicted in Figure 53. Note that this test point was flown near maximum thrust, not an insignificant factor for the Twin Otter. The trim condition is established for TEU, 8E=const2 the first 12 sec. It should be noted that the pilot was actually pulling 30 lbs at this trim point because the trim tab had run out of authority (see Section 5.3.4); in the plot, this offset has been

subtracted (FYE- FYEO). The pilot then pushed the yoke for- Cl_Tail ward to dive and gain speed. The pitch over portion was initi- ated at t = 27.5 sec with the elevator push TED. The elevator is returned TEU at t = 30 sec. The control force lightened, but clearly did not reverse. [The FYE data were clipped above a certain value, unfortunately.] Inspire of the pull-up, however, N_ continued to decrease and reached its mini- mum at 31 sec. Moreover, the aircraft continued to gain speed and pitched over to 0 = -28 °. Examination of the pressure belt data indicates that immediately after the elevator was pulled TEU, the tailplane experienced a full stall starting from the trailing edge. What is important to note here is that this pushover would have passed the CFR criterion, yet the tailplane had stalled. There are more advanced techniques, however, that can identify tailplane stall. They are introduced next.

NASA/TP--2000-209908 53 0.5 j 0

85 •

75 - - - '- -

45 I , :

20 _ " _ ] theta (de9) -,o 171- i°

-30 '

25 10 I _deIE(deg) © dEdef FYE(Ibs} ,#,_'_

0

-25 , --", _ _ ),,_ -50

-I0

-t5 OA2 (deg) -20

-25

,30 L

-0.4 , , ,

-0.6 - - - ...... t

-0,8

-1,2 - 0 S 10 I5 20 25 30 35 40 45 tim • 50

Figure 53. Example of a trailing edge tail stall. Baseline, 6F = 40 °, V = 0.94Vs = 50kts, CT = 0.23.

Advanced Analysis Methods To examine some of the more advanced analysis techniques, the example of increasing flap deflection will be used. Three pushover test points, all with the Failed Boot ice shape and achieved conditions of N_ = 0G and VTAS = 80kts were flown with three flap deflections, tSF = 0°, 10° and 20 °. This advanced analysis will revisit some co-plots and introduce some cross-plots.

Figure 54 depicts the tSE & FYE time histories for the three flap deflections. Note in particular the push and hold portions where the elevator went TED (&E increased). This occurred at the following times: for tSF = 0°; 2-5, 22-25, and 42-45 sec; for 6F = 10°; 14-17, 25-28 and 37-40 sec, and for SF = 20°; 22-24, 43-45 and 67-69 see. One can see that as the elevator was pushed forward, the control force remained a nominally constant push (-40 lbs) for 8F = 0% lightened to a point for t_F = 10°, and reversed, or crossed the neutral axis, for tSF = 20 °. [The circles appearing in Figure 54 are for reference in Figure 57.] The 6E & q time histories appear in Figure 55. One can see the pitch response is strongly damped for &F = 0 °, weakly damped for 10°, and undamped for 20 °.

NASA/TP--2000-209908 54 10 ...... r...... ooE 40

0 - 0

_u .5 -20

-15 ..... 60 0 5 10 15 20 25 30 35 40 45 50 55

15 60

uJ

-5 -20 -_o ...... i...... i...... ,o

0 5 I 0 I5 20 25 30 35 40 45 50 55 F ...... :...... :...... l ''°

o _ .._ _ o

-15 -EO 20 25 30 35 40 45 50 55 60 65 70 75

tim • (=)

Figure 54. Co-plots of ¢_E and FYE for pushover. Failed Boot, VTAS = 80kts, increasing flap deflection.

=_ -8 -I0 =r

-16 ...... L...... • ...... J- -20

0 5 10 15 20 25 30 35 40 45 50 55

16 20

o o

"o -8 -10 o-

-16 ...... 20

0 5 10 15 20 25 30 35 40 45 50 55

_ 8 10 _ 0 0 _w _ -8 -10 er

-16 -20 20 25 30 35 40 45 50 55 60 65 70 75 tim • (=)

Figure 55. Co-plots of _E and q for pushover. Same data base as in Figure 54.

NASA/TP--2000-209908 55 Perhaps a more elegant presentation of these data is via the cross-plot. The question posed is whether or not there is a distinct footprint in these cross-plots which would indicate problems. Before continuing with the effect of flap de- flection on a contaminated leading edge, a well-behaved non-critical pushover is examined first. A set of cross-plots for the Baseline case at 6F = 0° and V = 1.5Vs = lOOkts is presented in Figure 56. Time histories of these data are depicted in Figure 50; the particular pushover examined occurs between t = 53-68 sec. Note, for this well-behaved case, the character of these cross-plots, both the q vs. rE and FYE vs. _E, exhibit nearly flat straight lines.

50 15 5o

25 10 _ _ . ° . 2sl J 5 0 0 -25 ...... ir -5 It. ,,¢--2s i ..... -50 -10 -50i ...... -75 -15 -75 :

-0.4 0 0.4 0.8 1.2 -10 -5 0 5 lO -10 -5 0 5 10_

Nz (g) dE (deg) dE (deg)

Figure 56. Cross-plots for Baseline case at 8F = 0 ° and V = 1.5Vs.

Now let us return to Failed Boot data presented in Figure 54 through Figure 55. Figure 57 displays the three corre- sponding cross-plots of a single pushover - those identified by the circles in Figure 54. The SF = 0° case appears in the top row, 8F = 10° in the middle and SF = 20 ° in the bottom row. Like any phase plot, the temporal information is lost. However, the circles, which correspond to the circles found on the 8E & FYE co-plots from Figure 54, give some indication of events. The first circle marks the beginning of the maneuver - the maximum N z. After this point, the elevator is pushed TED and the 0G portion of the pushover begins. The second circle marks the point where the elevator is returned or pulled TEU. The arrows indicate the path from the first to the second circle. [Note that for the FYE vs. N_ cross-plot, only the push portion is indicated; only the second circle is included in these plots.]

Some of the interesting features of each cross-plot are discussed next. FYE vs. N_ : Data in this format have been presented by Ref. 14, and is used within the community fairly extensively. One key feature of this graph is to note on which side of the axis the elevator is returned. A circle in the upper half of the graph indicates a CFR. By noting the trend after the elevator is returned, one can also check for a stalled condi- tion. If, after crossing the trim force, the data trend continues to the left, then the aircraft has not responded to the elevator TEU command. A plot of FYE vs. q yields similar information. Q vs. SE : This graph highlights the pitch response to the elevator deflection. Of most interest is the lower right hand quadrant where the elevator is returned. For 8F = 0 °, the short period is evidenced by the overshoot in q while 8E is still moving TED. When the elevator is returned, q responds immediately. For 8F = 10°, the short period response is still evident, albeit much diminished. For 8F = 20 °, on the other hand, the elevator traveled over 10° before any ap- preciable decrease in the magnitude of q was noted. FYE vs. 6E : This cross-plot is an easy check of CFR. If, as the elevator returns across the trim force, the slope of the line is positive, then a CFR is indicated.

NASA/TP--2000-209908 56 75 _---- 15

50 " - - 50 ...... i 105 ...... i goo ._ 25 ...... i I i _.2 -s UJ >. 0 i >_ 0 LL b iii-25 -25 i " " -_5 t - - - ...... Z;...... -50 _ ...... i...... J -20 _ ...... -50 *0 4 04 {38 12 -10 -5 5 ]0 ]5 -15 -t0 -5 0 5 ]O 15 Nz (g) de1 E (deg) del E (deg/

75 _...... l 5 ,- ...... _ ...... = ...... 75

10 50 - - 50

_ 25 - - oi = 25 W "i ,/li -5 _ .... W >. o i >. 0 LL LL -_o ._ .... -25 -25 -15 ! ....

-20 -50

0,4 08 I 2 -I0 -5 5 I0 15 -15 -I0 -5 0 5 I0 15 Nz [g] del E (deg) der E (deg)

75 75 r-- ...... " ...... 15 _ ...... _ ...... , ......

I0 _ .... 50 - - 5O 5 _ .... I

25 - • _" 25

W W >_ o i--- >.. 0 u_ i LL -26 _ - o" -10 -_ -25 -15 _, ....

-50 _ ...... -2(1 ...... -5O

-0 4 0 4 08 1 2 -t -5 0 5 I0 15 -_5 -}o -5 o 5 }o _s Nz [g) del E (deg) del E (deg)

Figure 57. Cross-plots for same database as Figure 54. Top row is 6F = 0 °, middle 10° and bottom 20 °.

Recall the earlier discussion that the CFR criterion was inadequate for warning against a trailing edge stall. Figure 58 i]Iustrates the FYE vs. Nz and q vs. 6E cross-plots corresponding to the data set presented in Figure 53. From the FYE vs. Nz cross-plot, the elevator was clearly returned TEU prior to CFR; however, the subsequent negative slope of the line indicates that the aircraft did not respond immediately. The lower right-hand quadrant (indicated with the dashed box) of the q vs. 6E cross-plot immediately indicates a lack of pitch response.

25 5

r 0 0

n I ...I , _ -5 -25 III -10

II -50 -15

-75 -20

0 0.25 0.5 0.75 1 1.25 -15 -10 -5 0 5 10 15

Nz (g} dE (deg)

Figure 58. Cross-plots corresponding to data in Figure 53, the trailing edge stall pushover case.

NASA/TP--2000-209908 57 6.1.6 Elevator Doublet Elevator doublets demonstrate the effect of tailplane ice on the longitudinal stability and control. This maneuver is typically associated with a parameter estimation study. The discussion contained herein introduces the maneuver and discusses some of the aspects and features as related to pitch rate characteristics. Section 7 discusses the parameter estimation analysis. Descripti.on.: This maneuver consisted of a series of elevator inputs initiated from straight and level flight. The air- craft was trimmed for the target flap deflection, speed and thrust. One series of negative and positive deflections was made immediately followed by a second series whose amplitude was heightened and period shortened. Configuration: All ice shapes, flap deflections, and speeds were flown. For the majority of cases, the moderate thrust setting of Cr = 0.10 (0.05 per engine) was chosen; however, a higher setting of Cr = 0.18 was also flown for a very limited number of test points. For the Failed Boot and S&C ice shapes, elevator doublets were not flown beyond tSF = 30 °. Recall the corresponding pushovers were limited to t_F = 20 °. Anal sy__: The elevator doublet is directly compared to the pushovers previously shown. Figure 59 depicts an eleva- tor doublet for the same aircraft configuration and flight condition (Baseline, _F = 0 °, V = 1.5Vs) as the pushover in Figure 50. Likewise, the elevator doublet presented in Figure 60 is similar to that of the pushover in Figure 51, Failed Boot shape at tSF = 20 ° and V = 1.06Vs. Perhaps the most direct comparison is between the second push and elevator TED portion of the elevator doublet to the push and elevator TED portion of the pushover, Finally, the most critical elevator doublet is examined in Figure 61. Recall that for steady level flight, the most negative _ occurs for a high flap setting and high speed. Appendix C contains additional elevator doublet data.

As the t_E vs. FYE is the key relationship for the pushover, o°E vs. q is key for the elevator doublet. The primary question concerns whether the pitch response is damped or if it continues to increase after the elevator has been fro- zen, As an illustration of a damped response, examine q and t_E for the second "push and hold" portion of the eleva- tor doublet depicted between t = 7-9 sec in Figure 59. While the elevator was fixed at t_E = 0.5 °, the q-magnitude increased (q decreased), peaked and receded. This peak occurred when the short period was excited. The q- magnitude then lessened because the pitch response was damped. An example of an undamped response is provided by all four of the responses in Figure 60. Again, focussing on the last push between 7-8 sec, one sees that the q- magnitude increases continuously the whole time the elevator is fixed TED at tSE = -1 °. Only when the pilot returned the elevator to neutral did the q magnitude decrease.

For the non-critical case of Figure 59, the control force was flat and the pitch response damped. The minimum o_ was -1.8 °, and the minimum Ct_railwas -0.4. Similar to the corresponding pushover, this elevator doublet clearly passed all tests. For the more critical case in Figure 60, the control force lightened and came very close to reversing on the last elevator push (t = 7.75 sec). In addition, the pitch response is clearly undamped. The most negative _ = -7.5 and Cl..rail = -0.7.

NASA/TP---2000-209908 58 1.5 1 ' ' | _Nz l

ol

110 "

100

95 1 90 _ _ ,

10 - , • - ) 1 8 _ .' .... ' _ _ _theta (deg) I

o!

8 4O delE (deg) ..... FY E (Ibs) 4 ...... r:-':_...... J_--,.,...... 2O

0

-20

-8 -40

5 ' [- - , ..... J ? ...... ] - / _ J __ ', ...... Cl (deg/s) 1 -°i::::s - - uJ ./'f'iLJ\ . 5..... -10 - '- "" - ......

0.2

0

-0.2

-0.4

8 [deg) ! 6 J

4

2

0 2 4 6 8 10 12 time 14

¢

Figure 59. Elevator doublet for a clean tail, 8F = 0 °, V = 1.SVs = 100kts.

NASA/TP--2000-209908 59 Nz (9) 1.5 ...... , ...... --

t 0.5 ......

0 i

65 --VIAS (kts) 60 ......

55 ...... '-- _ ......

50 ...... , ......

45 .... T -

12

.4 I

14 _delE (deg) ...... FYE (Ibs) H T 40 ' - - - f_2 ...... itS-< ; .... 20

-7 - = - -" .... 20

-I4 -40

14

, --delE (deg) o7 ...... L ._,..-_. '<, _/ t-\_. ,,...... : f. -...... q (degls) ! -7 - = "'" -'--" _ " _ ...... I -14 I

-2 ' TAOA 1 (dag]

-6

-8

0.2

Cltail

-0.2

-0.4 -0.6 1 -0.8

12 I (deg)

106_8 ; i . _ . Alpha

4

2 4 6 8 10 12 14 tim a 16

Figure 60. Elevator doublet for the Failed Boot, 8F = 20 °, V = 1.1Vs = 55kts.

NASA/TP--2000-209908 60 Thepreviousdiscussionallowedforadirectcomparisonof plotsbetweenthepushoverandelevatordoubletgiven anaircraftconfiguration.Thiscomparisonisdiscussedin Section6.2.Nowtheq vs. ¢5Eco-plots and some corre- sponding cross-plots will be compared for the extremes of the elevator doublet parameter space. The time histories of these are depicted in Figure 61 for a clean and an iced shape. The Baseline configuration is compared to the worst-case ice shape, S&C, at flap deflections of zero and maximum flap deflection flown for the S&C ice shape, 30 °. The speeds were nominally 1.6Vs. Several items stand out from these plots. One is the relative similarity be- tween three of the comers and the drastic difference for the worst case - S&C, _F = 30 °. Note for both Baseline cases and the S&C _F = 0° cases that the pitch response is damped. In all cases the pitch response relaxes more on the nose down portion (SE moving positive or TED). For the S&C ice shape at _F = 30 °, however, the pitch re- sponse is undamped. The aircraft continues to increase its pitch rate magnitude until the elevator position is changed. The fluctuations of the SE plot indicate the difficulty in holding the elevator still. This is due to the unsteadiness of the separation bubble washing back and forth over the hinge point. Another thing to note is the elevator deflection at trim for the different ice shapes. For SF = 0 °, the elevator deflection is nominally the same, 8E = -2.5 °, regardless of ice shape. For _F = 30 °, however, _E rises from 7 ° for the Baseline to -1 ° for the S&C ice shape. The ice shape has caused an 8° loss of elevator authority. .

Cross-plots are presented in the next two figures. As with the pushover data, the circles in Figure 61 denote the por- tion of time being cross-plotted. Note that for the elevator doublet, the first circle occurs with the elevator deflecting TED, not with the maximum Nz, as was the case for the pushover maneuver. Both Figure 62 and Figure 63 bring out some features that are not as easily picked up from the time history co-plots. For example, both the Baseline at SF = 30 ° and S&C ice shape at ¢5F= 0° cross-plots indicate that the final push portion is only weakly damped.

_q Ideg/s) Ba=eline, dF = O, V = 1.5W= S & C, dF = O, V = 1.5Vs _tl (deg/s)

...... clelE (Oeg)

____j-: oo:';:oo' L O dEdef is 1 10 ...... 10 ...... ,

0 "_-......

-'15

0 5 I0 15 0 5 10 15

time'(s) time (s)

Baseline, dF = 30, V = 1.S V$ _ q (tleg/s) S & C, dF : 30, V = 1 .SVs _cl (deg/s)

_, delE (deg) I ..... delE (deg)

15 T...... Z ...... _......

0

s I i',11......

0 5 10 15 0 5 10 15 time (s) time (s)

Figure 61. Pitch response and damping characteristics for extreme elevator doublet cases.

NASA/TP--2000-209908 61 Baseline, dF = O, V = 1.5 Vs S & C, dF= 0, V = 1.5Vs

15 15 1...... _ ...... ,...... ] 10 ...... _ .... r .... 10 jo A

0

o" O" -5[ -10 L ...... I -15 -8 -4 0 4 8 1'_ 0 4 B 12

dalE (deg) dale (dag)

Baseline, dF = 30, V = 1.6 Vs S&C, dF=30, V= 1.6Vs

15 ...... • ...... t ] . , 10

A

o : - i ', o_iiiiii:,I ...... i _:'_:!:!il;i ii cr

_:o!...... : ...... -10

-t5 ......

-8 -4 0 4 8 12 -8 -4 0 4 8 12

dalE (dag) dalE (deg)

Figure 62. Cross-plots of q vs. _E for elevator doublet data from Figure 61.

Baseline, dF = 0, V = 1.5 Vs S & C, dF = 0, V = 1.5 V=

4O 40 ...... ,...... , ......

2oI ...... '...... :...... 20

0 &]

tlJ uJ >- -20 :i- Lk u.

-40 -40 _ ...... J -60 L -60 0 0.5 1 1.5 2 0.5 1 1.5 2

Nz (G) Nz (G)

Baseline, dF = 30, V = 1.6 Vs S & C, dF = 30, V = 1.6Vs

40 • 2oI .o60 I...... iJ,...... 1 A A == 20I ..... '...... _.--.:,...... J

W uJ >- >- U. -20 ...... '.... - - - u...... 20| ' \ \ / - -40 ...... 1 ...... _ E °I -40 l ...... J -60 ......

0 0.5 1 1.5 2 0,5 1 1.5 2

Nz (G) Nz (G)

Figure 63. Cross-plots of N=vs. FYE for elevator doublet data from Figure 61.

NASA/TP--2000-209908 62 6.1.7 Balked Landing This maneuver was born from the debriefing discussion during the first sequence of the Guest Pilot Workshop with the FAA and Transport Canada pilots (see Section 8). The question was raised concerning the distinction between stability and controllability, and its relationship to a typical piloting task. Thus, a closed-loop task was constructed to determine safe flight characteristics and controllability margins with an ice-contaminated tailplane. Description: The balked landing maneuver consisted of a simulated instrument approach and go-around. This ma- neuver was flown at altitude and with "heads down." To simulate course and glide slope corrections, commands to change both rate of decent (ROD) and heading were made every 20 sec all the while maintaining a constant 1.3Vs velocity. Specifically, heading changes of +5 ° off a reference heading, and RODs of-500 + 500 fpm were ordered. The RODs required the pilot to adjust the thrust requirements. The idealized flight path appears as the dashed line in Figure 64. The 20 sec intervals required the pilot to make fairly aggressive control and thrust inputs. At the conclu- sion of the approach portion of the task, a go-around was commanded. Unlike the maneuvers discussed previously (open loop), the ability to determine whether or not the pilot actually met the targets closed the task loop. Configuration: This maneuver was only flown with the Failed Boot ice shape with flaps set to _SF = 10 °, 20 ° and 30 °. Test points were conducted with V= 1.3Vs and thrust set as required (Cr = 0-0.24). Anal sy__: In addition to the idealized task requirements, Figure 64 illustrates flight data from two pilots. The data set depicted in the left column was flown by a NASA pilot. One can see that this pilot met the task criteria within the 20 sec allotted. The data set in the right column was flown by a 600-hour guest pilot. This guest pilot was able to stay on task for the relatively low thrust requirements but encountered difficulty when near maximum thrust had to be added to attain level flight. Instead of leveling off, the aircraft lost 280 ft of altitude when the flow at the tailplane separated (the safety pilot recovered the aircraft). It is felt that this comparison highlights the fabled "dark and stormy night" scenario where there is a considerable pilot workload in addition to some level of fatigue at the end of a flight. Should the tailplane start to stall, would the pilot be able to recognize the situation, and respond quickly and appropriately?

I500 1500

,_lllm Sire ul.ted _ _ ,_ll :ipmp_lol::: 1000 1000

E 500 Approach P'i I_ _F_ 500 | u |

_ 0 0 , l : ! , i o -5oo -500 1 Go Around m -1000 -1000

-1500 i - -1500

10 10

v 5 5

c 0 0 I I I

-5

I -10 -10

-15 -15

100 100 I i 90 l 90 l < 80 80

70 70

20 40 60 80 100 120 140 160 0 20 40 60 80 "100 120 140 'T60

tim • (8) tim • (o)

Figure 64. Balked landing tracking task.

NASA/TP--2000-209908 63 6.1.8 A Case of Turbulence The opportunity to document a "light to moderate" turbulence encounter arose during the guest pilot workshop. During a flight home, 16 sec of data were taken with the Failed Boot ice shape at zero flaps and 140kts. In Figure 65, it is compared to 13 sec of a steady, wings level flight for the same nominal conditions. Two key configuration dif- ferences are (1) a much higher thrust setting for the turbulent case, Cr = 0.15, as opposed to 0.09 for the quiescent case, and (2) significantly offset FYE values. For the quiescent case, the pilot was not allowed to change the trim from the lowest speed configuration. The corresponding difference in yoke force was -63 Ibs (push). What is plotted, then, is the difference (FYE - FYEo), where FYEo is the value at the start of the maneuver.

One can see the turbulence tossed the aircraft +0.2G, with a corresponding pitch rate oscillation of approximately +l°/s. Also illustrated is the effect of the turbulence on the elevator and tailplane angle of attack. The standard de- viation of the signal for the turbulence vs. quiescent case is 0.125 ° vs. 0.040 ° for d_E and 0.25 ° vs. 0.09 ° for _.

15

I 25

I

0 75 i ...... , .... ,.... O5

2

I

Q o" -1 _,,,_.y...... v_..:..,,_.:....

2

! 0 m -1 " ...... :.... :.... .2

2

-6

1,15

140 > 135

130

8 g" ,

-2

15

-0 5 I

15 w 10

5 w 0

-5 -10 2 4 e 8 !0 12 14 IE time 1_

Figure 65. Comparison between quiescent and turbulent environment. Failed Boot, &F = 0° flaps.

NASA/TP--2000-209908 64 6.2 Comparison Between the Pushover and Elevator Doublet The NASA/FAA Tailplane Icing Program contains the data to directly compare the pushover (PO) to the elevator doublet (ED). Six configurations were compared; two representative configurations are listed in Table 2. These con- figurations are the previously shown (1) Baseline, 6F = 0 °, V = 1.5Vs and (2) Failed Boot, 6/7 = 20 °, V = Vs. Recall that for the latter configuration, the Twin Otter failed both tests; it experienced a CFR during the pushover and an undamped pitch rate response during the elevator doublet. The excursion range, minimum values and a difference are compared. The range is the difference between the maximum and minimum values recorded during the test point, the minimum is simply the minimum value recorded, and the difference is that between the trim and minimum values. Note that under the Comparison of Minimums header, elevator deflection (delE) data appear. The values recorded are those of the elevator at its most TED position, i.e., the maximum value.

For the range data, ED/PO represents the relative excursion percentage of the elevator doublet compared to the pushover. In all cases, except for Cm, the elevator doublet range excursions are less than those of the pushover. Fur- thermore, the pitch, speed and altitude excursions are significantly lower, less than 20%. Neither the N z, q nor FYE ranges exceed those of the pushover. The Cm range is larger for the elevator doublet because the nature of the ma- neuver requires a step input to the elevator.

For the minimum and difference data, PO-ED represents the difference between the two maneuvers. While the tail- plane angles of attack are not as negative during the elevator doublet as for the pushover and their difference from trim is not as great, the tail loading factors, minimum CLr, ir and ACLTa,t, are comparable. Note from the delE data that during the pushover, the tailplane is much less cambered than for the elevator doublet. Recall the less the cam- ber, the greater the stall margin (Figure 12). Therefore, it is the combination of g and SE that determines the air- craft's longitudinal stability, not large o_ alone. This combination is manifest in the CLr_il value.

Baseline, 5F = O, V = 1.5Vs = lOOkts Failed Boot_ 5F = 20, V = Vs = 55kts

Comparison of Ranges: Comparison of Ranges: PO :ED ED/PO ED!PO theta (deg) 70 7 0.10 theta (deg) 50 10 0.20 FYE (Ibs) 94 58 0.62 FYE (Ibs) 80 50 0.63 VlAS (kts) 60 4 0.07 VlAS (kts) 30 5 0.17 alt (ft) 925 38 0.04 air (ft) 480 60 0.13 Nz (g) 2.5 1 0.40 Nz (g) 1.5 0.3 0.20 q (deg/s) 23 17 0.74 q (deg/s) 28 14 0.50 delE (deg) 13 9 0.69 delE (deg) 23 15 0.65 Cm (about CG) o.11 0.4 4.00 Cm (about CG) 0.5 0.7 1.40

Comparison of Minimums: Comparison of Minimums: ...... PO ED PQ-ED ::PO ,ED P_ED alpha (deg) -4 0.7 -4.7 ialpha (deg) -9 4 -13.0 TAOA (deg) -8 -3 -5.(3 TAOA (deg) -14 -7 -7.0 Cl_Tail -0.3 -0.4 0.1 Cl_Tail -0.9 -0.7 -0.2 delE (max) 6 0.5 delE (max) 10 -1

Comparison of Differences: Comparison of Differences: PO ED PO-ED ..:;. PO ED PO, ED .a, ,,, A(alpha) deg 6.6 2.7 3.9 A(alpha) deg 12.2 ' 2.6 9.6 A(TAOA) de£ 6.6 2.5 4.1 IA(TAOA) de 9 7 2.5 4.5 &(Cl_Tail) 0.2 0.2 (3 _,(Cl_Tail) 0.5 0.4 0.1

Table 2. Pushover - elevator doublet comparison.

NASA/TP--2000-209908 65 Nowthatsomequantitativedifferencesbetweenthesetwomaneuvershavebeenestablished,thenextquestioncon- cernsapass/failcriterion.Toattemptto answerthisquestion,aPass/FailMapwasconstructedforboththecontrol forceandpitchresponsecriteria.Thesejudgementsweremadeononesetofpushovermaneuvers.Caution:thismap attemptstopresenta4-Ddatasetin a2-Dplane.Thethreeinputsaret_F, minimum VTAS/Vs and minimum G-load for each maneuver. The output has three levels: pass, marginal and fail for each criterion. In color copies of this map, note that blue italic text corresponds to a pass, black to a marginal response and bold red to a fail. Careful inspection of this map will reveal a nearly one-to-one correspondence between the two criteria; i.e., they are nearly interchange- able.

The control force rating examined the relationship between FYE and t_E. The pass level is denoted by P; this means the control force remained nominally flat while the elevator was deflected TED. The marginal level is denoted CFL; this means the control force lightened while the elevator was moving TED. This level was further subdivided into two categories: a weak CFL (wCFL), which means that the force lightened, but only to some point shy of neutral, and a strong CFL (sCFL), which means that the force crossed the axis within a second of the elevator returning TEU. Finally for the fail level, a CFR means that the force crossed the axis while the elevator was fixed TED. This level is also presented in bold face in Table 3. For the pitch response rating, q and tSE were examined. For this rating, the three categories were damped D for pass, weakly damped wD for marginal and undamped U for fail.

Table 3 presents a tSF vs. VTAS/Vs chart for each of the ice shapes (the Baseline case passed all configurations). The map on this plane consists of the above rating with the G-load next to it in parentheses. For example, for the Failed Boot ice shape at _F = 20 ° and VTAS/Vs = 1.4, the entry is CF(L.5,R.2). This means that for a 0.5G pushover, a CFL was experienced, but when it went to 0.2G, the response was a CFR. For the pitch response of the exact same data set, the entry is wD(.5), U(.2). This means that the response was weakly damped during the 0.5G pushover, and undamped for the 0.2G maneuver.

As can be seen from Table 3, these rating criteria yield nearly identical results. The pitch response criteria will yield the same information as the control force response.

NASA/TP--2000-209908 66 Control Force Ra:in.nq 1.8 P (OG) P (0) 1.7 1.6 1.5 CFR (.4) CFR (.5) 1.4 sCFL (.2) CFR (.S) 1.3 CFR(,t5) CFR (.4) t.2 P(O) P(O) CFL(.5w,0s) 1.1 P (.6) CFR (.6) 1.0 sCFL (.2) CFR(.4) 0 10 20 30 40

Pitch Rate Ratihq 1.8 O (OG) 0 (0) ...... 1.7 1.6 1.5 u {.4) u (.s) 1.4 u {°2) u (.5) 1.3 wD (.3, .15) u (.4) 1.2 D (0) 0 (0) wD(.5), U(0) 1.1 0(.6) u (,e) 1.0 wD (.2) u (,4) 0 10 20 30 40 dF

_! Force Ratinq ...... Control Force Ratinq ..... 1.8 1.8 CFR(0) 1.7 P(OG) 1.7 CFR (0) 1.6 P(.2) 1.6 P (OG) CF(L.4,R.2) 1.5 1.5 CF(L.4,R0) CFR (0) 1.4 CF(L.5,R.2) 1.4 CF(L.S,R0) 1.3 CFL (0) CFL (.5) 1.3 CFL (.4,.2) CFL (.6) 1.2 CFR (.2,0) 1.2 CFa (0) 1.1 P(o) P(.6), sCFL(.4) 1.1 P (0) CFR 64,0) 1.0 CFR (.2) 1.0 CFR 0 10 20 0 10 i5 20

pitch Rate Rating 1.8 1.8 wD (0) 1.7 D (OG) 1.7 u (o) 1.6 D (.2) 1.6 D (OG) D(.4),wD(.2) 1.5 1.5 U(.4,0) u (o) 1.4 wD(.5), U(.2) 1.4 u(.&o) 1.3 wD (0) wD (.5) 1.3 u (.4,2) u (.6) 1.2 U (.2,0) 1.2 u (o) 1.1 D (0) D(.6), wD(.4) 1.1 D (0) U (,4,0) 1.0 u (.2) 1.0 U (0) 0 10 20 0 10 15 2O dF dF

Table 3. Pass/Fail maps for both the control force and pitch response ratings during a pushover maneuver.

NASA/TP--2000-209908 67 6.3 Repeatability Analysis for the Pushover

This section introduces an analysis of the pushover data introduced in Section 6.1.5 and includes a precision and accuracy study.

An important question regarding the pushover is how precisely it must be flown. If, for example, a point is flown five knots below the lowest target speed and a control force reversal is experienced, should the aircraft be certified if it would have passed at the target speed? What if the point is flown five knots too fast? In general, what magnitude of error is tolerable, and what is not? It is beyond the provision of this report to provide comment on such guidelines. What can be discussed herein, however, are the variations experienced by one NASA pilot flying the Twin Otter aircraft.

The nature in which most of these pushovers were flown allows for such a repeatability analysis. Within each test point three parabolas were flown, each targeting the same conditions of pitch rate (and correspondingly, vertical ac- celeration), speed and pitch angle. If the three parabolas can be superimposed on each other and synchronized to a distinguishing feature, then a point-by-point repeatability analysis can be conducted.

The distinguishing feature was chosen to be the vertical acceleration Nz. This choice of an output parameter was based on the fact that the criteria for the pushover maneuver is written to this parameter, and the fact that it has obvi- ous demarcations between maneuvers. The choice of the input parameter, elevator deflection, was not chosen be- cause it was not "clean" in that the elevator was moved differently for different configurations. This choice would have required more filters or decisions to properly define the maneuver.

Reported herein are precision and accuracy analyses for a range of critical and non-critical configurations. This sec- tion starts with a detailed description of the analysis.

6.3.1 Description of the Analysis The goal is to superimpose the three maneuvers of each test point on one axis. It was decided that a maneuver shall begin and end with the pull up, i.e., from maximum Nz to maximum N z. This choice also puts the dynamics of interest in the middle of the maneuver. Once the three maneuvers are properly defined, each is scaled to run from 0 to 1. Typically, maneuvers lasted 9-15 sec; within a given configuration, the variation was usually within 2 sec. The three individual maneuvers, plus their average, are depicted in Figure 66 for a non-critical configuration that demonstrated no CFR, and therefore easily "passed" the certification criteria: Inter-Cycle Ice with tSF = 0° and V = 1.5Vs.

It became apparent that simply identifying the end points was insufficient; the minimum scaled Nz appeared any- where from 50% to 80% of the total time; therefore, it needed to be fixed as well. The optimal place to fix the scaled N__rainis the average N__,_in location of the three maneuvers under consideration. For the example in Figure 66, the choice of 50% was made. The full procedure is then: 1) Find the actual times of the two Nz_,,_, to and tl, to determine the period T. 2) Find the actual time of Nz_,,_, train. 3) To fix the scaled time of N__rainat 0.5, subtract 0,5 Tfrom trainto find the new start time, to'. Add Tto the new start time to find the new end time, t_'. 4) With the bounds of the maneuver determined, scale the actual time, t, in this manner: tscaled = (t -- to')/(tl'- to'). And, of course, tscaled (N__ra,,) = 0.5 = (train -- to')/(tl'- to').

Once the traces are scaled and synchronized, statistical analyses may be performed. To make the task of directly comparing all three scaled time traces more manageable, the time axis was subdivided into equidistant bins. For this exercise, the number of maneuvers, N,_n = 3, and the number of bins, Nb,n = 20, which left the number of points per bin, Ns_z = 70, Let Q(i,j,k) represent the i'h point of the fh maneuver in the k'h bin for a quantity Q. The scaled data within a bin were averaged, and the average value was placed at the center of that bin. This was done for each of the three maneuvers.

NASA/TP--2000-209908 68 N ¸.

siz i=l

Within a bin, the values of the three averages were themselves averaged to determine the overall average for that time step.

N : , Q(k )= Nma-----_ j=l

As seen in Figure 66, the overall average for each quantity is represented with the thick blue line. This line is simply a straight-line fit through the Q . The error bars within each bin, s(k), come from taking the standard deviation of the three differences between each trace average and the overall average.

sCk)= 2.,j=, !--,Nma.r--(N,,_n - 1 2],,2

The final error, _, is the average of all of the bins.

To isolate the target portion of the maneuver and eliminate the entry and exit variations, the values centered around

the Nz_mi n bin can be averaged. The quantity _ ta, is defined similarly to s, but only averaged over the immediate target neighborhood. Let ko define the bin which contains Nz_,,_,, and n define the number of bins in the neighbor- hood.

ko+-". - 1 z Star =- Zs(k) n

k=ko- 2

6.3.2 Description of the Results

For the Inter-Cycle ice case presented in Figure 66, the s tar value was averaged from tscolea= 0.375-0.625. This neighborhood is depicted with the vertical dashed lines.

In Table 4, both the _ t_, and s- errors are listed for the above and several other configurations. These are divided into cases where no CFR was experienced (non-critical) and cases where CFR was experienced (critical). Focussing on the s tar values, Table 4 suggests that N, was very repeatable regardless of the maneuver criticality. On the other hand, the repeatability of q, V and to some extent 0 was effected by criticality. Note that the speed was normalized by its target, V/Vtars, t.

With the average values established, the accuracy of achieving the targets may now be addressed. The accuracy shall be defined as the difference (Dif) between the average of the ko bin (Actual) and the target quantity (Target), IQ (ko) - Qt_rml. Again, the speed was normalized by its target, %Dif = I(V (ko) - Vtarm)/Vt_gat* 100. These values are pre- sented in Table 5. The results suggest the pilot was able to achieve the target N, to within 0.1 G, regardless of diffi- culty performing the maneuver. The pilot averaged 5% accuracy in hitting V_rS_t. Theta was clearly the most difficult target to meet. The average of the values reported is -18 ° off target. Typically, to achieve the target N z, the Twin Otter had to be in a dive, not at the horizon.

NASA/TP--2000-209908 69 2,5 2 ::,"...... I ...... I --',...... "...... /:

0.5 "-...... 7 F ...... 43-56s | 0 r , - , . • • • j,,,,,,_I.,,--AVG Nz r' -o.5 .... | --| ' - J _6o I- ! "; 14o...... t ...... I- - -"...... :.:

u_ 120 ...... JL ...... ___._ " "_ ...... 28-39s ....."" = ..... _ m , ...... 43-56s > 100 ...... | - -', ..... _AVG VIAS 80 ' I ' -....7--- I " 40 I- I _14.28, _...... _ ...... _i ...... 28-3gs s 2o-_...... I ...... I , _...... 43-s8, 0 _ "_'*_L ,,,=,_===AVG theta -- _".... '..... '...... ,'- :l . -20 - -

-40 -- 8 I I | ...... _ ...... " ...... _ ...... ! : --14-288 o_ 4 ...... " 28 39s

0 . -==lk,==AVG del E,

-8-- J-- ! so ..... ! I

"z _---,,_..._: J t : __ ...... _

0 , . .- ,_, ' | ' ' ' | " _4-28s uJ . . _ ,_ ...... , f_#//\ __ __ _ J...... 28-39s _ " _...... ,._...7:._.._-...... ,_-, I,--41,,--AVGFYE -Ioo-- -I- -I ' ' 12 I I ' ...... ======,: ® 0 "' t i ' _ ¢ _14-28s- -4 - _\ ...... : ...... "__._.,,_ff_../...... 28-3_s cr - " ,__ -- _,,,_L-,-_ -...... 43-56s -8 ..... _,., ...... _ " "' ...... ,,,,,,,,I,--AVG Q -12-- l I 8 I I "o 4J ...... '_" - I ' 4 28_ L, < ....._- I...... _3-s8=l _- -8_---- ...... I" :- I--,--T., AOAi 8-- I !

m _AVG alpha[

-8 -- 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 08 0.g 1 t_l ¢=led

Figure 66. Three maneuvers (thin lines) and their average (thick line) vs. scaled time for an Inter-Cycle ice, _F = 0°, V = 1.5Vs = 100kts test point. Error bars are indicated on select traces.

NASA/TP--2000-209908 70 Configuration Nz(G) q(dug/s) V/Vtarget (%) theta (dug) Ice Shape, _SF,V/Vs s t_ [ s ] s _ [ s- ] s"tar [ S ] S"_ [ S ] NO CFR

IC, 0°, 1.5 0.12 [0.17] 1.2 [1.5] 4.2 [3.t] 2.8 [3.1]

FB, 0°, 1.5 0.09 [0.24] 0.9 [2.0] 5.0 [7.2] 2.5 [4.0]

SC, 0°, 1.5 0.05 [0.21] 0.5 [1.9] 3.1 [3.6] 1.4 [2.6]

Ba, 20 °, 1.0 0.06 [0.07] 1.0 [1.1] 2.0 [4.9] 1.5 [1.5]

IC, 20 °, 1.0 0.06 [0.05] 1.2 [1.0] 3.5 [2.4] 1.7 [1.3]

IC, 40 °, 1.0 0.07 [0.09] 1.9 [1.7] 9.1 [7.5] 2.6 [3.2]

FB, 20 °, 1.15 0.07 [0.12] 2.4 [2.3] 2.6 [5.2] 1.8 [3.2]

FB, 20 °, 1.0 0.08 [0.19] 2.9 [3.4] 5.1 [8.2] 8.3 [5.6]

SC, 15 °, 1.0 0.24 [0.24] 3.5 [3.4] 6.7 [6.9] 2.5 [4.2]

SC, 20 °, 1.0 0.10 [0.12] 2.3 [1.7] 8.5 [9.3] 3.5 [3.5]

Table 4. Repeatability analysis - precision levels for several flight conditions to minumum Nz,

Configuration Nz (G) VIAS (kts) theta (dug) Ice Shape, 8F, V/Vs Actual [Target Dif Actual Target [ %Dif Actual ]Target [ Die No CFR ....

IC, 0°, 1.5 0.07 0 0,07 108 100 8.0 -4.2 3 -7.2 , ,,,/ FB, 0°, 1.5 -0.07 0 -0,07 102 100 2.0 -12.5 3 -15.5

SC, 0°, 1.5 0.05 0 0.05 105 100 5.0 5.5 3 2.5 Ba, 20 °, 1.0 0.22 0.17 0105 57 55 .... 3.6 -25.9 3 -2819

IC, 20 °, 1.0 0.31 0.25 0,06 60 55 9.1 -20.7 3 _23.7

...... CUR ...... ¸¸ ¸¸¸ ¸¸ ¸¸¸¸ i

IC, 40 °, 1.0 0.58 0.50 0108 55 55 0.0 -15.7 -4 _11.7

FB, 20 °, 1.15 0.15 0.17 66 65 1.8 -18.4 -0.5 -17.9

FB, 20 °, 1.0 0.24 0.17 0.07 56 55 1.3 -18.3 3 -21.3 sc, 15o, 1.0 0.01 0 010i 58 55 5.4 -23.3 5 -28.3

SC, 20 °, 1.0 0.34 0.25 0.09 62 55 12.7 -23.9 3 -26.9

Table 5. Repeatability analysis - accuracy of achieving the target values for non-critical and critical cases.

NASA/TP--2000-209908 71 6.4 Dynamic Maneuvering Conclusions

Several types of dynamic maneuvers were flown during the NASA/FAA Tailplane Icing Program. The primary inter- est was to study the pushover maneuver. A potential alternate maneuver, the elevator doublet, was also studied. In addition to these two maneuvers which intentionally manipulated the elevator, several other conceptually simplified maneuvers were flown. These maneuvers were three types of transitions - increasing flaps, speed and thrust - flown to isolate their effect of the tailplane. Each of these transitions, by itself, drove an ice-contaminated Twin Otter hori- zontal tailplanc toward stall. It should be noted that a thrust increase might not affect other aircraft as adversely; the phenomenon observed for the Twin Otter is believed to be caused by the thrust line positioned above the CG.

A tail stall occurs, of course, when the tailplane angle of attack, o_, exceeds the stall angle, o_ s,an. The elevator de- flection, BE, also plays a role in determining the stall angle. For example, the stall point of a more cambered airfoil " occurs with more downward lift but at a reduced stall angle. A high flap deflection always corresponds to a more negative tail angle of attack. However, it was also noted that there are opposite trends in the speed criteria for ex- treme tail angles depending on whether the aircraft is in a steady state or dynamic maneuver. For steady, wings level flight, high speed will drive O_omore negative. During maneuvering, on the other hand, the dynamic component be- comes dominant, Aat, This corresponds to low speeds and high pitch rates. With these thoughts in mind, the single variable paths to stall will be discussed first, then the pushover and elevator doublet maneuvers.

6.4.1 Transition Maneuvers The flap transition is clearly a path toward stall. Flap deflection drives the tailplane toward stall by making a, uc more negative. The downwash remains nominally constant because CL_,ucremains nominally constant. To maintain speed, the elevator moves TED which increases the stall margin slightly, but the dominant factor is the substantial decrease in ¢x,uc, which drives the o_ more negative.

The mechanism driving the speed transition toward stall, on the other hand, is less obvious. A microscopic view of the data is necessary since looking at only the gross trends of a nominally constant t_ and a slight 6E deflection TED suggest that the stall margin should increase. A careful inspection of the oscillations in the o_ data coupled with the 6E data, revealed that the most negative extremes of the _ trace drive the tail toward stall, whereas the least negative extremes increase the stall margin. For this extreme case, the tailplane is going into and out of stall.

For the full stall of the thrust transition, the features are again subtle. The ty-,u_decreased about 2° with a slight de- crease (-11 ° to -12 °) in the _. To maintain speed, the elevator was required to move TEU. With this, the stall mar- gin decreased and the plane plunged into a full tail stall. ' .....

6.4.2 Dynamic Maneuvers The main focus of the study was on the pushover. To this end, nearly 400 pushover maneuvers were flown in a man- ner consistent with the FAA ANM-100 memorandum E'T°e Bookmrknot d_n,,_d..Note that in the ANM-100 description, the pushover maneuver begins with the elevator TED deflection and ends with the minimum load factor. The pass/fail criteria are that (1) a push control force is required throughout the maneuver (to minimum-G) and (2) "the airplane should demonstrate suitable controllability and ma_neuverability throughout the maneuver with no force re- versal and no tendency to diverge in pitch". For this series of flight tests, it was interpreted that the "maneuver" ended when the elevator was returned TEU.

Time history analyses immediately indicate how well the pilot met the target conditions of the maneuver (N_, V, 0) and whether or not a control force reversal (CFR) was experienced. This was determined by inspecting the co-plot of _E & FYE; a CFR occurred when the control force crossed the neutral axis (or trim point) while the elevator was fixed TED, prior to the elevator return TEU. As expected, the incidence of CFRs for the pushover increased with:

NASA/TP--2000-209908 72 (1) increasingiceshapeseverity, (2) increasingflapdeflection, (3) decreasingnormalloadfactor,and (4) decreasingspeed. Datawerepresenteddemonstratingthatmerelylookingfor aCFRis agoodbutnotfullyinclusiveindicatorof whetherornotthetailplanestalledorhadundampedcharacteristics.In fact,thesedataconfirmedtheANM-100 guidancetonotethe"tendencytodivergeinpitch".Thatis,inspectingq vs. t for an undamped response seemed to be a slightly more robust indicator of a destabilized tailplane than the CFR. More conclusive results are obtained by the following more advanced analyses: (1) q vs. t_E - Examine the pitch rate response after the elevator is returned TEU. If stalled, the elevator could travel quite a distance before the pitch rate starts to return to zero. This would also appear as a reduced ele- vator effectiveness, C,._. (2) N z vs. FYE - Check the slope of the line after the elevator is returned TEU and after it crosses the trim point. If stalled, the aircraft will continue toward 0G regardless of the force applied by the pilot to return toward 1G flight.

Also presented were two examples of special cases. One was of a "shaped" maneuver that erroneously passed the CFR criteria because the elevator was deflected slightly but sufficiently to maintain flow attachment. The other ex- ample passed the CFR criteria but stalled after the elevator deflection TEU. The flow remained attached until the elevator was pulled TEU and the tailplane became more cambered. It should be noted that this example failed the pitch rate response criterion. Analysis of the q vs. tSE and Nz vs. FYE cross-plots helped resolve this issue. More ex- amples of pushovers can be found in Appendix C.

Since the pushover was a difficult maneuver to fly accurately, a repeatability analysis was synthesized to determine the variations in how well the NASA pilot was able repeat the maneuver and how accurately he simultaneously hit the three targets. The severity of the ice shape did not seem to play a major role in how accurately he met the targets, although it did, to some extent, affect the repeatability of the velocity and pitch attitude. The pilot met the Nz ,,i, ob- jective within +0.1G. The simultaneous VIAS varied but was typically within about 5% of V,_rg_,.The most difficult target to meet simultaneously with N_ ,,i, was/9. For the Twin Otter, this was typically -20 ° beyond the trim 0.

Part of the request from the FAA was to investigate potential alternate maneuvers to the pushover. To this end, ele- vator doublets were flown. Of course, the elevator doublet is typically flown for parameter estimation studies, but it seems to have applicability for tailplane stall certification as well; this maneuver could also destabilize the tailplane. Appendix C also contains a number of elevator doublet time histories. One issue that needs to be addressed for the elevator doublet, or any alternate maneuver, is that of "equivalency" to the pushover. That is, what is the relevant parameter for comparison?

For the elevator doublet, as with the pushover, examining either the t_E & FYE or 6E & q co-plots yielded essentially the same conclusion as to whether or not the tailplane destabilized; again, the pitch rate response, q, seems to be a more robust indicator. Based on an undamped response, it was found, as expected, that the elevator doublet was most likely to fail with: (1) increasing ice shape severity, (2) increasing flap deflection, and (3) increasing speed.

The pushover and elevator doublet maneuvers were compared directly for the same aircraft configurations and flight conditions. A test point that destabilized the tailplane for the pushover also destabilized the tailplane for the elevator doublet, even though the pushover generated much more negative tailplane angles of attack. What became apparent is that the extreme o_ values alone are not an adequate indicator of impending tail stall, but rather the combination of o_ and _E that drive the tailplane toward stall. For the Twin Otter, this seems to translate to an equivalency parameter of Ct_T_,t.This should be checked with other aircraft, however.

NASA/TP--2000-209908 73 A furtherinvestigationofthedifference,if any,betweenthecontrolforceandpitchresponsecriteriawasconducted. Thetwocriteriawerecompareddirectlyusingthepushoverdata.A nearlyone-to-onecorrespondencewasobserved foralliceshapes,flapdeflectionsandspeeds.Inotherwords,foragivendataset,if aCFRisobserved,thensoisan undampedpitchresponse.FortheTwinOtter,thecontrol force and pitch rate response criteria yield nearly the same results. Further analysis of any discrepancies indicated that the pitch response was the more robust meas- ure.

In summary, it was found that for the Twin Otter, the same certification pass/fail decisions could be made using the elevator doublet rather than the pushover. In general, the elevator doublet is an easier and safer maneuver to fly. Summarized in Table 6 are the advantages and disadvantages of the pushover compared to the elevator doublet. It would be of great interest to investigate the response of other aircraft to both the pushover and elevator dou- blet.

...... _jl_ Ilrmrnn _[_!1_ I_'rl _

....' ...... _ightio 0G does n0irequire insirumenta-Excursions arei _+015G ...... tion to detect. ' • Much less risky than Pushover; Advantages • Control Force Reversal criteria straight- much easier to recover from a stall. forward to detect PROVIDED the maneu- • Easier to fly; less pilot variability; ver is flown properly. not susceptible to shaping. • Excursions are 1 + 1.1G • Hardware required to measure • Flight to less than 0G requires instrumen- pitch rate and elevator deflection. tation for documentation. • Difficult to hit all three targets (Nz, V & 0) as aircraft tracks through horizon. Disadvantages • Possible to manipulate the elevator to keep flow attached. • If tail stalls, aircraft will be in a difficult position for recovery (0 -_ -40 °, near VFE). • Concern for fluid system operations during 0G flight (e.g., hydraulics, fuel, oil).

Table 6. Comparison of the advantages and disadvantages for both the pushover and elevator doublet.

NASA/TP--2000-209908 74 7.0 Parameter Estimation Analysis

Elevator doublet maneuver test points were added to the test matrix in Flight 97-39 as a potential alternate maneuver for discerning icing effects on tailplane performance. The benefit of the maneuver was twofold. Prior experience with the elevator doublet suggested it could also be used for pilot evaluation of the aircraft response. In addition, the elevator doublet maneuver can be used in parameter estimation analyses to derive mathematical aerodynamic models and stability and control derivatives of the aircraft. The focus of this section is to provide the results of the parameter estimation analysis on the Twin Otter with the various levels of contamination.

7.1 Flight Test Procedures

The elevator doublet maneuver was a series of elevator inputs initiated from straight and level flight. The aircraft was trimmed for the target flap deflection, speed and thrust. The elevator input was a near perfect classic square wave consisting of two doublet cycles followed by 6 to 10 sec of no ccyntrol input. Elevator doublets were conducted with all ice shapes and covered the speed range for the flap settings tested. For the majority of cases, the moderate thrust setting of Cr = 0.10 (0.05 per engine) was chosen; however, a higher setting of Cr = 0.18 was also flown for a very limited number of test points. For the Failed Boot and S&C ice shapes, elevator doublets were not flown beyond 6F = 30 °.

Since the elevator doublet was also used for pilot evaluation, the maneuvers were of relatively large amplitude. The incremental angle of attack range typically varied +3 ° to +5°; while, the incremental normal acceleration often varied +0.5G from 1G flight. The maneuvers are also discussed in Section 6.1.6, and were very similar to those previously performed for and documented in Ref. 3.

Table 7 shows the tailplane ice and aircraft configuration where data maneuvers were performed. A complete matrix of data was obtained for the 0°, 10°, 20 °, and 30 ° flap settings. Data were also obtained for the Baseline and Inter- Cycle ice configurations with the flaps at 40 °. For the S&C ice configuration, additional data were obtained at the 15° and 25 ° flap settings. Except for the Baseline data with the higher thrust coefficient of 0.18, an "X" in the table represents 3 to 4 maneuvers over the angle of attack range. At the higher thrust setting, the "x" represents 1 to 2 ma- neuvers at the high angle of attack for each flap setting. Because of the limited quantity of the data with higher thrust, it will not be presented with the results. The remaining data were all obtained with the same lower thrust setting. For the aft CG baseline data, a linear aerodynamic correcT;ion was applied to the static stability parameter and to the ele- vator control parameter to correct them to the forward CG for data presentation.

Data will be presented at the 0°, 10°, 20 °, 30 °, and 40 ° flap settings as the effects of tailplane icing are more promi- nent when the various ice configurations are viewed at a common flap setting.

Flap Deflection (deg) Tail Configuration CG CT 0 10 15 20 25 30 40 Clean AFT 0.10 X X X X X Clean FWD 0.18 x x x x x Clean FWD 0.10 X X X X X Inter-Cycle ice FWD 0.10 X X X X X Failed Boot Ice FWD 0.10 X X X X S&C Ice FWD 0.10 X X x X x X I

Table 7. Conditions for the parameter estimation maneuvers.

NASA/TP--2000-209908 75 7.2 Method of Analysis for Parameter Estimation

A parameter estimation analysis of the flight data was performed using version 2.3.4 of pEst, a parameter estimation program similar to an earlier version in Ref. 15. The pEst program is an interactive, nonlinear program for the analy- sis of dynamic systems. This analysis used the standard pEst user routines that define the aircraft six-degree-of- freedom nonlinear set of differential equations. However, full coupling of the equations was not needed as the ma- neuvers were performed in a decoupled manner. This allowed the aerodynamic part of the equations to be split into a longitudinal subset, while the inertial part of the equations maintained their full nonlinearity. For the longitudinal subset, the elevator time history measurement was the only input control variable. The output control variables con- sisted of state variables and additional response variables. The state variables used were angle of attack, pitch angle, and pitch rate. The additional time responses used were normal acceleration, axial acceleration, and pitch angular acceleration. The measured lateral-directional response variables were also used to supplement the inertial part of the equations. Dynamic pressure and velocity varied with time throughout the analysis.

The parameters used to define the longitudinal aerodynamic model contained the standard longitudinal partial de- rivatives that are often referred to as stability and control derivatives. In addition, the angle of attack squared pa- rameters for pitching moment and axial force were used to better model the large size of the maneuvers. The results of the pEst analysis were sets of longitudinal parameters for each maneuver. The resulting parameters are a best value averaged over the oscillation range of the maneuver. For example, for maneuvers where angle of attack oscil- lated +5 °, the resulting parameters are the averaged values over the 10° angle of attack range. For plotting presenta- tion, the parameters are then plotted verses the average angle of attack for each maneuver.

It should be noted that the analysis process works best when the model used is consistent with the vehicle character- istics. The mostly linear aerodynamic model used for the analysis is very adequate to describe the piecewise linear behavior exhibited by a conventional configuration below its stall angles of attack. However, when the wing flaps are deployed much more than about 20 °, the resulting separated flow induces a modeling error that shows up as scatter in the parameter results.

The analysis process consisted of an initial analysis of all the maneuvers followed by a final analysis. The initial analysis used starting values from the prior maneuver in sequence and yielded a complete set of the longitudinal de- rivatives. The data were then plotted as a function of angle of attack for each flap position and tail ice condition. Trends from the initial analysis were used as starting values and predicted reference values for the final analysis. A low weighting was placed on the predicted values to lessen the scatter in the final results. Care was taken to limit the cost due to the weighted predictions to approximately 2% of the total cost function. Data from a small number of the maneuvers were not used because the analysis was not able to obtain an acceptable fit between the measured and the calculated responses.

7.3 Results

The stability and control characteristics of the Twin Otter have been previously documented in Ref. 3. This analysis expands on the prior results by presenting the effects of tailplane icing using three ice shapes and with wing flap ex- tensions as high t_F = 40 °.

A very basic review of the effects of tailplane icing on the flight data will enhance interpretation of the parameter results. Moving the CG forward requires the horizontal tail to produce more down load to balance (trim) the aircraft. Deploying the flaps creates a negative pitching moment that also requires the horizontal tail to produce more down load to trim the aircraft. Thus for a given CG and flap position, the elevator or engine thrust can be modulated to attain the desired angle of attack or velocity. As most of the parameter estimation data were at the same CG, thrust was at a constant setting, and the flaps were at discrete settings, the elevator was the only remaining longitudinal control. Thus, as the various ice shapes were placed on the horizontal tall, a rough measure of the degradation in the tail's ability to generate a down force was the elevator position (tSE). However, the elevator required for trim varies with flap deflection because the flaps also change the flow angle into the horizontal tail. While this negates using the absolute value of t_E as metric, the incremental value between the ice free tail and the various ice shapes can still be a

NASA/TP--2000-209908 76 validmetricwithinthelinearloadrangeof thehorizontaltail.Incremental8Eisnotavalidmetricwhenthehori- zontaltailisclosetostall.FortheforwardCGandthrustcoefficientCT--- 0.10, the additional _E required to main- tain angle of attack as a function of flap position is shown in Figure 67 for the various ice shapes. Note that the Baseline configuration is the zero elevator axis. Given the qualifications of this plot, it is noted that the Inter-Cycle ice had no effect (within the plot resolution) on the horizontal tail's ability to produce down load until the flaps were deflected past SF = 20 °. Likewise the Failed Boot and S&C ice shapes had only a limited effect until the flaps were deflected beyond SF = 10°. While the figure is only valid prior to tailplane stall, the maneuvers for parameter esti- mation were also conducted without encountering tailplane stall. Assuming the parameter variation with ice shapes follows the trends of Figure 67, only slight changes in the parameters will result at the lo.wer flap settings where the ice shapes have only a limited effect on the function of the horizontal tail.

-5 -} -_"" Intercycle __ _,__ | -6 | -._-,,,,Failed Boot "_ |

0 10 20 30 40

Flap Position (deg)

Figure 67. Additional elevator required to trim.

While the tailplane ice caused a significant degradation of the general flying qualities in the areas of buffet and the associated control wheel shaking, its effects were explicit in a relatively small number of parameters. The elevator control effectiveness, C,,aE, was the principal parameter. The effects were also apparent in the longitudinal static sta- bility parameter, C_a_ and the pitch damping parameter, C,,,q. While lesser effects could be inferred from the axial and normal force parameters, the effects were generally of limited significance to the flying qualities and, thus, are not be presented. Angle of attack typically oscillated over a 6° to 10° range; while the elevator control inputs oscil- lated over a 10 ° to 14° range. The analysis returns the best value of the parameter averaged over the oscillation range. Thus as tailplane stall was approached, the parameters were only marginally degraded as only the nose down limits of the elevator doublets were near tailplane stall. While smaller doublets would have been more desirable for the parameter analysis, the doublets analyzed for this study needed to be large for pilot evaluation.

The resulting parameters obtained near SF = 0° setting are shown in Figure 68. Of significance is the approximately

8% reduction in C,_sE with either the Failed Boot or S&C ice shapes. The effects on C_a and Cmq due to the ice shapes are within the scatter of the data. The reduction in Cma at angles of attack greater than 7° is related to the wing characteristics.

Figure 69 presents the data obtained with &F = 10°. The data again show an approximate 8% reduction in C,,6_ with either the Failed Boot or S&C ice shapes, and there is also about a 4% reduction in C,,8_ due to the Inter-Cycle ice shape. Slight reductions in C,_ and Cmq are also starting to become apparent.

NASA/TP--2000-209908 77 WiththetSF = 20 ° flap setting, (Figure 70) the parameter variations with ice shape become more dramatic. The re- ductions in the elevator control effectiveness parameter, C,,_z, are about a 10% for the Inter-Cycle ice, 17% for the Failed Boot ice, and 27% for the S&C ice. The static stability parameter, C,,c,, becomes unstable (positive values) at the higher angles of attack (lower speeds) for both the Failed Boot and the S&C ice shapes. C,,q shows reduced pitch damping for the S&C ice shapes at the lower angles of attack.

At the tSF = 30 ° flap setting (Figure 71), the effects of ice on the horizontal tail significantly degrade the flying qualities. C,._e is reduced about 10% for the Inter-Cycle ice shape, 22% for the Failed Boot ice shape, and 33% for the S&C ice shape. The effects on C,,c, are not explicit for the Inter-Cycle or Failed Boot ice shapes. For the S&C ice shape, C,,,_ becomes statically unstable above about -2 ° angle of attack. The analysis of C,,q is strongly effected by the separated flow from the flaps and contains enough scatter to mask any trends in the data.

With 6F = 40 ° (Figure 72), data were obtained for only the Baseline configuration and with the Inter-Cycle ice shape. With the Inter-Cycle ice shape, C,,SE shows a strong reduction as angle of attack is decreased (increasing air- speeds). This reduction is about 30% at -4 ° angle of attack. The trends in C,,a and C,,,q are not apparent.

NASA/TP--2000-209908 78 0.04

0.02 ...... , ...... J ......

L_ -0.02

-0.04

-0.06

3 4 5 6 7 8 9 10 11 12

-0.005

-0.01

-0.015

-0.02 ...... J ...... _ ...... rOE

-0.025

-0.03

-0.035 -._,',,,,- ...... • ,_ "-__"':'"CE_.:,_=-,,_ "_....

-0.04 3 4 5 6 7 8 9 10 11 12

0 Baseline _ Intercycle Ice :_, Failed Boot & S&C

-10 ...... ' ...... !

-20

L_

-30

-4O

J -5O

3 4 5 6 Aircraft Angle of Attack(deg) 9 10 11 12

Figure 68. Icing effects on pitching moment derivatives, ¢_F = 0°.

NASA/TP--2000-209908 79 0.04

0.02 ......

-0.02 . _- ...... -_-_n .... ' ..... i o : ..... :'- : - - o- -,., ......

-0.04 i ...... :..... i ...... i..... i ..... i i J -0.06 : -2 -1 0 1 2 3 4 5 6 7

I J -0.005 - ) )

-0.01 ...... I

-0.015 ...... )

J -0.02 ......

-0.025 ......

-o.o3 ---,_ ...... :.... i ...... ',..... _ - -o.o35 _ _-:-:_--'_-- - '..... _ .....

0 , -0.04 ...... -2 -1 0 1 2 3 4 5 6 7

/ o Baseline {3 Intercycle Ice .: Failed Boot A S&C I

-10

-20 ...... J

-30 ...... 0 - - o ....Q,,_ o D 0 -40_ ......

.50 _ -2 -1 0 1 Aircraft Angle of Attack (deg) 4 5 6

Figure 69. Icing effects on pitching moment derivatives, 6F = 10 °.

NASAfrP--2000-209908 80 0.04

0.02 ...... ,...... : --_ ...... , ..... _ L.>I=,---

0 i :

-0.02 i e _._o_-"'_' o o

-0.04

-0.06 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8

-0.005 ......

-0.01 ......

-0.015 ......

-0.02 ...... -o.o25--,-..:..:_:_=.=..=..:;: ...... -oo..... -0.035 _ ...... ,...... O ...... '--- O_ G-C) =- U -0.04 ...... -5 -4 -3 -2 -1 0 I 2 3 4 5 6 7 B

o Baseline _ Intercycle Ice o Failed Boot z_ S&C]

-10

-20 J -3O

-4O

-5O -5 -4 -3 -2 -1 Aircraft Angle of Attack (deg) 4 5 6 7

Figure 70, Icing effects on pitching moment derivatives, b'F = 20 °.

NASA/TP--2000-209908 81 o.o2i°°4F-...... _...... i"...... t ...... _.i_' ......

d .... -- o ...... i

-0.06 a...... _'...... '-...... i -6 -5 -4 -3 -2 -1 " 0 1 2 3 4 5

0

-O.0O5

-0.01

-0.015

o_ -o.o2 ...... ; ...... A ......

-0.025

-0.03

-0.035 ,----'...... '...... 6--- l- o ...... 0 -_ -0.04 .J -6 -5 -4 -3 -2 -1 0 1 2 3 4 5

t o Baseline _ Intercycle Ice _ Failed Boot A S&C

-10

-2O

J -3O ._.Q.o A ;_._ _' A o

-4O _ _ -, ......

-5O -6 -5 -4 -3 -2 Aircraft Angle of Attack (deg) 2 3 4 5

Figure 71. Icing effects on pitching moment derivatives, _/F = 30 °.

NASA/TP--2000-209908 82 0.04 ...... I...... ,...... _......

0.02

J

-0.02 ...... '...... _ ...... t l/ ......

0 D o 0 0 0 0 ...... 0 ..... _ 0 -0.04

-0.06 _ ...... -8 -7 -6 -5 -4 -3 -2 -1 0 2

-0.005 ......

-0.01

-0.015

a_ o' -0.02 -0.025 ...... @......

-0.03

-0.035

-0.04 -8 -7 -6 -5 -4 -3 -2 -1 0 1 2

0 ' I o Baseline _ Intercycle Ice

-10 ......

-20 i ...... J i O -30 _ _ _ 0 ..... O O ' -_ ...... --O ol

-40

-50 -B -7 -6 -5 Aircraft Angle of Attack (deg) -1 0 1 2

Figure 72. Icing effects on pitching moment derivatives, 5F = 40 °.

NASA/TP--2000-209908 83 7.4 Parameter Estimation Conclusions

A parameter analysis of the flight data was used to quantify the degradation resulting from Inter-Cycle, Failed Boot, and the S&C ice shapes attached to the leading edge of the Twin Otter tailplane. The analysis was conducted as a function of angle of attack with flap settings as high as t_F = 40 °.

The Inter-Cycle ice shape was shown to slightly degrade the baseline airplane. Up through tSF = 10°, the effects of the Inter-Cycle ice were minimal. With SF = 20 ° and tSF = 30 ° the elevator control effectiveness was reduced about 10%, while changes in static stability and pitch damping were small. With _F = 40 °, the Inter-Cycle ice shape began to significantly decrease the elevator control effectiveness.

The Failed Boot ice shape showed a larger degradation. The ice shape effects were minimal up through t_F = 10° flap setting with the primary effect being about an 8% decrease in elevator control effectiveness. The elevator control effectiveness droped by 17% at tSF = 20 ° and by 27% at tSF = 30 °. Significant reductions in static stability are also apparent,

The S&C ice shape yielded the largest degradation shown in the parameters. At 8F = 0° and tSF = 10°, the primary effect was about an 8% decrease in elevator control effectiveness. The elevator control effectiveness was reduced by 27% at 8F = 20 ° and by 33% at 8F = 30 °. At both 6F = 20 ° and 6F = 30°s, the vehicle became strongly statically unstable with increasing angle of attack. Trends in pitch damping were less certain due.to data scatter.

NASA/TP--2000-209908 84 8. Guest Pilot Workshop

At the conclusion of the flight tests, a Guest Pilot Workshop was conducted to provide a forum to rapidly dissemi- nate valuable lessons learned in the Tailplane Icing Program to the user community and to gather flight data using non-TIP pilots. An international group representing various facets of the aviation industry - aviation regulatory agencies, aircraft manufacturers and aviation media pilots/reporters - were invited. In total, 15 guest pilots and engi- neers had the opportunity to fly the NASA Twin Otter with the Failed Boot ice shape on the tail and experience the unique flying qualities of an aircraft with an ice-contaminated tailplane. The purpose of this section is to provide an overview of the workshop and show results from a handling qualities assessment of the simulated approach and balked landing maneuver under specified conditions.

8.1 Attendees A total of 15 guest pilots and engineers attended the workshop. Table 8 lists the affiliation of each guest pilot that attended the workshop. Due to the relatively short duration time of each flight, the workshop was divided into three sequences to reduce the wait time for each guest pilot. The attendees for the first sequence were the FAA, Transport Canada, and non-TIP NASA pilots. Attendees for the second sequence covered the airframe manufacturers. The final sequence was devoted to the aviation media pilots.

Sequence 1 ] Sequence 2 Sequence 3 FAA, Test Pilot [ Bombardier, Test Pilot "Airline Pilot", Pilot FAA, Test Pilot [ Bombardier, Design Engineer/Pilot "Airline Pilot", Pilot/Reporter Transport Canada, Test Pilot Cessna, Design Engineer/Pilot "Aviation Week", Pilot/Reporter NASA, Test Pilot Raytheon, Design Engineer/Pilot "Commercial & Business Aviation", Pilot/Reporter Raytheon, Test Pilot "Flying Magazine", Pilot/Reporter "Professional Pilot", Pilot/Reporter.

Table 8. Guest pilot attendees.

8.2 Guest Pilot Workshop Flight Test Procedures and Maneuvers

Guest pilots performed each maneuver listed below from the left seat with a NASA safety pilot occupying the right seat. The data acquisition system was running during all maneuvers. The video cameras recorded the horizon, pilot control inputs as well as comments, and tailplane tuft activity while performing test maneuvers.

In order to provide a useful experience for multiple pilots in one flight session, a truncated test matrix was developed from the full test matrix to demonstrate the findings of the TIP. The maneuvers were: • Data Compatibility • Flap Transitions and Elevator Doublets • Steady Heading Sideslips • Wind up turns • Constant Airspeed Thrust Transitions • Pushover Maneuvers • Simulated Instrument Approach and Balked Landing (For more detail on maneuver description, see Section 6).

At the conclusion of the demonstration flight, debrieflngs were held and recorded with the intent of capturing imme- diate pilot impressions. NASA received written reports on the workshop from most participants that contained the following information on: • unusual or unexpected flight characteristics observed during the course of the demonstration program • suggestions to improve the technical content of NASA' s tailplane icing research program

NASA/TP--2000-209908 85 • impressionsofthepushovermaneuvervs.thesimulatedapproachtaskorrepeatdoubletmaneuversforassessing safeflightcharacteristics

8.3 Handling Qualities Assessment

The simulated approach and balked landing maneuver was developed after the first sequence of guest pilots. The discussion in those debriefings suggested that there needed to be a task-oriented maneuver to assess the reduced tail- plane performance. The simulated instrument approach and balked landing maneuver was provided two mission task elements (MTE) that could be used for handling qualities assessment. The maneuver was flown at altitude and with "heads down." To simulate course and glide slope corrections, commands to change both rate of descent (ROD) and heading were made every 20 seconds all the while maintaining a constant 1.3Vs velocity. Specifically, heading changes of +5 ° off a reference heading, and RODs of 0, 500, or 1000 ft/min were ordered. The expected tolerance in ROD was +I00 ft/min. The RODs required the pilot to adjust the throttle setting. See Figure 64 for the idealized flight path (dashed line). The 20 sec intervals required the pilot to make fairly aggressive control and thrust inputs. At the conclusion of the simulated approach MTE, a go-around was commanded. The go-around was accomplished by increasing thrust and raising the nose without changing the flap configuration. After the pilot accomplished a steady positive rate of ascent, the flaps were raised. The pilot immediately ranked the Twin Otter's flying qualities based on the Cooper-Harper handling qualities rating scale (Figure 73) for the approach MTE and go-around MTE. This maneuver was only flown with the Failed Boot ice shape and with flaps set to t_F = 20 ° and 30 °,

Aircraft Demands on the pilot in 8elected task Charectarl stlcs or required operation RATING !

Good: negligible IPilot compensationno_ a factor for desired Is it satisfactory deficiencies without improvement? DesirableEXCellent:HighlyperformancePil°tcompensationnota factorfordesired [ Deficlendes IFair: soma mildly IMinimsl pilotcompensation required for the 111 warrant lunpleasant deflcendes .l.pIdesired.e._ormancperformancee improvement

YES defieencies compensa_on , 4 l Moderately objectionahie Adequate performancerequlros considerabls Is deficendes p ot compensation IVer,/objectionable Adequate performance requtres extensk,e pilot attainable with a tolerable 3ilot workload? Deficiencies I Ideticen_es compensation require I

YES

Major deficiencies IIntense lYd_ compensationIs required to retain [Maio rdeficiencies icontrbl Is it controllable? Improvement _ Control w_llbe lost during soma portion of the [ __1 mandatory ItMaj°rdeficiencieslrequred operation I '01

PILOT DECISIONS

Figure 73. Cooper-Harper handling qualities rating flow chart.

The handling qualities ratings (HQR) for five of the guest pilots and the average value are shown in Figure 74. Sepa- rate ratings were given for the simulated instrument approach and go-around MTE's with the aircraft configured with flaps at tSF = 20 ° and at t_F = 30 °. By doing so, it is clearly seen that the handling qualities of aircraft with tSF = 20 °

NASA/TP--2000-209908 86 aremuchbetterthanthehandlingqualitiesoftheaircraftwithattSF = 30 °. Likewise, in both flap cases examined, it can be seen that the go-around MTE required more pilot compensation than the approach MTE.

HQR for Tracking Task

9

6 r- |m m

r_-- ...... G ...... I- •

0 Sim Appr Go Around Sim Appr Go Around

• 1 2 3 3 7

& 2 2 3 6 8

4, 3 3 6 8 10 o .m 6 8 a. O 4 2 3

[] 5 2 4 4 8

-41_ AVG 2.2 3.8 5.4 8.2

5F = 20 ° 5F = 30 °

Figure 74. Cooper-Harper ratings for approach and go-around MTE.

The significance of these HQRs is the agreement of the pass/fail criteria set up in the elevator doublet analysis. With the Failed Boot ice shape on the tail, CFR's were experienced during elevator doublets with tSF = 30 °, but not with tSF = 20 °. Likewise, the HQR for the 6F = 30 ° suggested that there were major deficiencies in the aircraft character- istics that required considerable pilot compensation to control the aircraft during the go-around task. However, the HQR for the tSF = 20 ° case, suggested that there were minor but annoying deficiencies in the aircraft characteristics.

NASA/TP--2000-209908 87

9. Concluding Remarks

The purpose of this section is to provide general comments and conclusions on the entire project and draw connec- tions across the previous sections.

9.1 Data Analysis Techniques

The flight data acquired in the TIP included steady-state 1G test points and dynamic maneuvering test points about the pitch axis with the intention of increasing the steady-state tailplane angle of attack. Data from the steady-state test points were analyzed through visual inspection of trends of a variety of parameters. Likewise, data from the dynamic maneuvering test points were also analyzed through visual inspection of time history plots, co-plots, and cross-plots of selected variables. In addition to these visual inspection analysis techniques, a computational analysis method was used with a subset of the dynamic maneuvering test points to estimate stability and control coefficients for a mathe- matical model that was representative of the test aircraft. Lastly, a different subset of the dynamic maneuvering test points was analyzed using visual inspections of pilot rankings of the aircraft flying qualities.

Each of these methods proved useful in some capacity by depicting the degradation in iced tailplane performance and the resultant degradation in flight dynamics and handling qualities. • For the steady state data, two key methods were developed that indicated degradation in tailplane performance and control anomalies. First, the difference in the elevator deflection (¢5E) required for speed for the clean and iced case for consistent aircraft configuration (CG, _F, Cr, landing gear, etc) clearly indicated significant tail- plane performance losses. Second, the difference in elevator hinge moment (Cne) required for speed for a clean and iced case for a consistent aircraft configuration, indicated tailplane performance losses and provided insight as to impending controllability problems. • For all of the dynamic maneuvering data, visual inspection of time history co-plots of stick force and elevator deflection (FYE & SE), and pitch rate and elevator deflection (q & SE), were key in the assessment of control force anomalies and the reduced pitch damping characteristics of an iced tailplane. This technique was improved by examining cross-plots of stick force versus normal acceleration (FYE vs. N_) and pitch rate versus elevator deflection (q vs. _SE). These techniques clearly indicated control force lightening and reversal, reduced pitch damping when the elevator was in the TED position, and reduced elevator effectiveness. • For the elevator doublet maneuvers, the estimation of stability and control parameters proved useful by provid- ing numerical values reflecting the reduced stability and control coefficients due to tailplane icing. The results clearly showed a reduced elevator effectiveness (C,,_) and static longitudinal stability (Cm_,) as ice shape sever- ity increased and tailplane AOA increased due to flap deflection. This analysis method did not reveal significant changes in the pitch damping derivative (Cmq), which were revealed by the visual techniques described above. This is likely a result of the time periods examined for each analysis method. The technique used by pEst esti- mated stability and control coefficients that best represented the entire doublet maneuver - both pitch up and pitch down. The visual technique focused on the part of the maneuver when the elevator was fixed TED and the tailplane AOA was most negative. Since the pest technique fitted the whole maneuver, the values may not be representative of events that took place within only part of the maneuver. Therefore, the Cmq results from pEst may not accurately reflect the damping characteristics at the most critical condition for the tailplane. To improve these results, maneuvers with smaller amplitudes could be developed specifically to explore this flight regime. • For the simulated approach and balked landing, the handling qualities results provided numerical indicators of the reduction of the pilot's ability to accomplish a specified task with an iced tailplane. These results are very significant since they provide insight into the ultimate question with an iced tailplane - how well can a pilot ac- complish a task with the aircraft in the potentially critical condition? Inherently, the ratings provide feedback on whether the aircraft has deficiencies that are tolerable or require improvement. Another significant aspect of this analysis method is that it employs a maneuver that is part of the normal flight operation.

NASA/TP--2000-209908 89 9.2 Maneuvers

The maneuvers chosen in this flight project ranged from steady state to dynamic maneuvers. Each type had value as revealed through the analysis techniques described above, The dynamic maneuvers provided a means to push the beyond the steady state values c_0 by some increment, Ac_. The increment depended on a_ns, V and q and was largest in magnitude for the minimum-G pushover maneuver at low speeds. However, for a given ice shape, the tailplane stall margin is not only a function of the _ but also _E. With the pushover maneuver, the elevator is TED, which by itself increases the tail stall margin. But since the Ac_ is large in the pushover maneuver at low speeds, the combined change in c_ and ¢5Eindicated ice-contaminated tail stall characteristics typically at lower flap settings than the other dynamic and steady state maneuvers. This may appear as an added level of conservatism, but it should be clearly noted that the tSE during the pushover is not in a position likely to be observed in normal flight operations. The ele- vator doublet provided information on tallplane performance degradation equivalent to the pushover at similar con- figurations by analyzing the pitch response characteristics. During the push part of the elevator doublet, 8E was less trailing edge down, and exhibited stalling tendencies at a less negative o_ than the pushover. Thrust transitions and the balked landing maneuvers were performed with a more appropriate 6E and also revealed deficiencies due to the tailplane icing, but the aircraft was configured with larger flap settings than the pushover. In summary, each of these maneuvers can provide an indication of reduced tailplane performance. The difficulty is not in the choice of the ma- neuver, but in determining the threshold of acceptability that is universally applicable and widely acceptable. Deter- mining thresholds was never the intent of this project and the nature of these results are far too limited to suggest those thresholds.

9.3 Aircraft Configuration

Results presented in this report are for one specific aircraft - a modified DHC-6 Twin Otter, and as such, the abso- lute values in parameters such as _ are limited to this aircraft. Likewise, it is not possible without further study to comment on aspects of tailplane design or other aircraft features such as wing and tail positions, wing and tail sweep and taper, airfoil sections, or propulsion system type. Generalizations on design aspects such as those would require a different means than the type of experimental testing performed in this project.

However, there are some aspects of aircraft configuration to which comments can be made. Since the data were col- lected over a wide range of configurations on the Twin Otter, some general insights were gained on critical and non- critical aircraft configurations for tailplane stall with respect to a few of the maneuvers flown. These generalizations are listed in Table 9. In all cases, the forward CG and maximum flap deflection are critical to achieving the greatest demand on the horizontal tallplane. Likewise, the aft CG and minimum flap deflection were non-critical in all cases. Airspeeds varied depending on the maneuver used for testing. Regarding thrust settings, it is logical that the critical- ity of high thrust was caused by the vertical displacement of the thrust line above the CG in the case of the Twin Ot- ter. So the Cr listed in the table are likely to be limited to similar high thrust line aircraft. This table may serve to reduce tests for tailpiane icing Sensitivity at the non-critical aircraft configuration.

Maneuver I Critical Non-Critical Steady wings level 1G CG=fwd, cSF=max, V=Vre, Cr=max CG=aft, _F=0, V=Vs, Cr=min Pushover CG=fwd, _SF=max, V=Vs, Cr=max, N_=min CG=aft, _SF=0, V=V,,_is,, Cr=min, Nz>min Elevator doublet CG=fwd, SF=max, V=Vre, Cr=max CG=aft, _SF=0, V=Vs, Cr=min Thrust transition [ CG=fwd, SF=max, V=Vre, Cr=max CG=aft, _SF=0, V=Vs, Cr=rnin

Table 9. Critical and non-crltical aircraft c0nfignraii0_fo_ ta-]lplanestall.

NASA/TP--2000-209908 90 9.4 Closing Commentary

The body of knowledge on ice-contaminated tailplane stall has been significantly advanced through the NASA/FAA Tailplane Icing Program by gaining a better understanding of the basic aeroperformance issues of a tailplane and the aerodynamic characteristics of the pushover and several alternate maneuvers. Although the progress made in this program was noteworthy, it did not provide answers to some of the key questions that remain about ice-contaminated tailplane stall. Some of these questions are: • What aircraft design aspects yield a greater or lesser susceptibility to ICTS? • What experimental techniques or computational tools can be used in the design phases to screen susceptibility to ICTS? • What aspects of the ice contamination cause critical degradation of tailplane performance? • What are universally applicable and widely acceptable pass/fail criteria for flight test maneuvers that screen sus- ceptibility to ICTS?

Although these are very difficult questions to provide absolute answers, efforts are in progress at NASA, the FAA and other organizations to improve the state of our knowledge and perhaps give guidance to answering some of these difficult questions.

NASAFFP--2000-209908 91

10. References

1 Ratvasky, T.P., NASA/FAA/OSU Tailplane Icing Program Requirements Document, May 1994.

2 Ranaudo, R.J., Baterson, J.G., Reehorst, A.L., Bond, T.H. & O'Mara, T.M., Determination of Longitudinal Aerodynamic Derivatives Using Flight Data From and Icing Research Aircraft, AIAA Paper 89-0754, NASA TM 101427, 1989.

3 Ratvasky, T.P. & Ranaudo, R.J., Icing Effects on Aircraft Stability and Control Determined from Flight Data, Preliminary Results. AIAA Paper 93-0398, NASA TM 105977, 1993.

4 Jordan, J.L., Platz, S.J. & Schinstock, W.C., Flight Test Report of the NASA Icing Research Airplane, KSR 86- 01, Kohlman Systems Research, Inc., Lawrence KS, NASA CR-179515, 1986.

5 Pfeiffer, N.J., Wake Rake Studies Behind a Swept Surface, CanardAircraft, Beechcraft Report, Sept. 1993.

6 Etkin, Bernard, Dynamics of Fliaht - Stability and Control, John Wiley & Sons, 1982.

7 Perkins, C.D. & Hage, R.E., Airplane Performance Stability and Control, John Wiley & Sons, 1949.

8 Ingelman:Sundberg, M. & Trunov, O.K., Wind Tunnel Investigation of the Hazardous Tail Stall Due to Icing, Swedish-Soviet Working Group Report No. JR-2, 1979.

9 Hiltner, D., McKee, M. & La No6, K., DHC-6 Twin Otter Tailplane Airfoil Section Testing In The Ohio State University 7'xlO' Wind Tunnel, NASA/CR--2000-209921/VOL1, 2000.

10 Gregorek, G.M., Dreese, J.J., _ La No6, K., Additional Testing of the DHC-6 Twin Otter Iced Airfoil Section at the Ohio State University 7'x10' Low Speed Wind Tunnel, NASA/CR--2000-209921/VOL2, 2000.

I 1 Dickson, B., Aircraft Stability and Control for Pilots and Engineers, Sir Isaac Pitman & Sons LTD, 1968.

12 Obert, E., Low-Speed Stability and Control Characteristics of Transport Aircraft With Particular Reference to Tailplane Design, AGARD Take-off and Landing, Jan 01, 1975.

13 Wojnar, Ronald T., Recommended Method of Identification, Susceptibility to Ice Contaminated Tailplane Stall, FAA Memorandum ANM-100, Apr 29, 1994.

14 Trunov, O.K. & Ingelman-Sundberg, M., On the Problem of Horizontal Tail Stall Due to Ice, Sw_ish-Soviet Working Group Report No. JR-3, 1985.

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

NASA/TP--2000-209908 93

Appendix A. Flight Test Summaries, Research Aircraft, Ice Shapes and Instrumentation

Flight Date Flight No. Tail Config ..... CG Comment 14-Sep-95 FLT 95-12 Baseline fwd A/C systems check in 0G 27-Sep-95 FLT 95-13 Baseline fwd SWL, steady pull up 28-Sep-95 FLT 95-14 Baseline fwd ISWL 29-Sep-95 FLT 95-15 Baseline fwd SWL 02-Oct-95 FLT 95-16 Baseline fwd SWL, DC added 03-Oct-95 FLT 95-17 Baseline fw_d_ SWL_ CHK added 10-Oct-95 FLT 95-18 Baseline fwd CHK points only, called back 10-Oct-95 FLT 95-19 Baseline fwd SWL, CTL=0 11 -Oct-95 FLT 95-20 Baseline fwd SHSS 11-Oct-95 FLT 95-21 Baseline fwd SHSS 12-Oct-95 FLT 95-22 Baseline fwd SHSS 17-Oct-95 FLT 95-23 Baseline fwd PO 18-Oct-95 FLT 95-24 Baseline fwd PO 19-Oct-95 FLT 95-25 Baseline fwd PO 20-Oct-95 FLT 95-26 Baseline fwd PO 23-Oct-95 FLT 95-27 Baseline fwd PO CTL=0 23-Oct-95 FLT 95-28 Baseline fwd PO, TPR, WT 25-Oct-95 FLT 95-29 Baseline fvyd SWL_ SHSS 26-Oct-95 FLT 95-30 Baseline fwd AIt maneuvers, SR, go-around_ PO max thrust 26-Oct-95 FLT 95-31 Baseline fwd AIt maneuvers, SR, go-around, PO max thrust 30-Oct-95 FLT 95-32 Baseline fwd air-air video document of maneuvers

Table A-1. 1995 flight test summary.

Flight Date Flight No:, Tell Conflg CG Comment 10-J ul-97 FLT 97-31 Baseline fwd Trailing Cone Fit 23-Jul-97 FLT 97-32 Baseline fwd bad ZOCs 23-Jui-97 FLT 97-33 Baseline fwd badZOCs 24-Jul-97 FLT 97-34 Baseline fwd bad ZOCs 25-Jul-97 FLT 97-35 Baseline fwd bad ZOCs 19-Aug-97 FLT 97-36 Baseline fwd ZOC's repaired 20-Aug-97 FLT 97-37 Baseline fwd 21-Aug-97 FLT 97-38 Baseline aft 21-Aug-97 FLT 97-39 Baseline aft 26-Aug-97 FLT 97-40 Inter-Cycle ice fwd 02-Sep-97 FLT 97-41 Inter-Cycle ice fwd thrust transition developed 05-Sep-97 FLT 97-42 Inter-Cycle ice fwd 05-Sep-97 FLT 97-43 Inter-Cycle ice fwd 13-Sep-97 FLT 97-44 Failed Boot Ice fwd 17-Sep-97 FLT 97-45 Failed Boot Ice fwd full tail stall 18-Sep-97 FLT 97-46 Failed Boot Ice fwd 18-Sep-97 FLT 97-47 Failed Boot Ice fwd air-air video document of maneuvers 22-Sep-97 FLT 97-48 S&C Ice fwd 22-Sep-97 FLT 97-49 S&C Ice fwd 24-Sep-97 FLT 97-50 S&C Ice fwd 24-Sep-97 FLT 97-51 S&C Ice fwd 26-Sep__7 FLT 97-52 Baseline fwd limited repeat w/good ZOCs 06-Oct-97 FLT 97-53 Failed Boot Ice fwd Guest Pilot Workshop I 07-Oct-97 FLT 97-54 Failed Boot Ice fwd Guest Pilot Workshop I 15-Oct-97 FLT 97-55 Failed Boot Ice fwd develo p balked landing 21-Oct-97 FLT 97-56 Failed Boot Ice fwd Guest Pilot Workshop II 23-Oct-97 FLT 97-57 Failed Boot Ice fwd Guest Pilot Workshop II 28-Oct-97 FLT 97-58 Failed Boot Ice fwd G u_estPilot Workshop III 28-Oct-97 FLT 97-59 Failed Boot Ice fwd Guest Pilot Workshop III 30-Oct-97 FLT 97-60 Failed Boot Ice fwd Guest Pilot Workshop III

Table A-2. 1997 flight test summary.

NASA/TP--2000-209908 95 characteristic Low High Mass, kg 4,5'i0' 4,970 Inertia IX, kg-m 2 26,190 26,660 IY, kg-m z 33,460 34,650 IZ kg-m 2 47,920 51,650 IXZ, kg-m 2 1,490 1,560

Area, m (ft) 39.02 (422.5) Aspect Ratio 10.06 Span, m (ft) 19.81 (65.0) Mean Geometric Chord, m (ft) 1.98 (6.5) Airfoil Section "DeHavilland High Lift" 17% thickness Horizontal Tail: Area, m (ft) 9.10 (98.18) Aspect ratio 4.35 Span, m (ft) 6.30 (20.67) Mean geometric chord, m (fi) 1.45 (4.75) Airfoil section NACA 63A213 Tail Volume, Vn 0.91

Table A-3. Physical characteristics of research aircraft.

Parameter Sensor Range R_lution (16 bit) AX & AY Sundstrand QA-700 + 1G 0.00003G AZ Sundstrand QA-700 +3G, -1G 0.00006G p, q, r Humphrey RG02-2324-I + 60 °Is 0.0018 °/s Humphrey VG24-0636-I + 90° 0.0027 ° 0 Humphrey VG24-0636-1 + 60° 0.0018 ° otaircraft Rosemount 858 probe + 15° NA [Baircraft Rosemount 858 probe _+15° NA V aircraft, airspeed Rosemount 542K2 0 to 190 knots 0.01 knots Altitude Rosemoum 542K2 0 to 15,000 fi 8.2 ft OAT Rosemount 102AU1P -20 ° to 30° C 0.041 ° C

6AL & _AR SAC series 160 +19°, -16 ° 0.0091 ° _SE SAC series 160 +14°, -260 0,0128 ° /SR SAC series 160 +_t6° 0.0080 ° /SF SAC series 160 0 to 40° 0.0091 ° Tail Flow Probe Scannivalve __0.72 psi 4.4 e'5 psi ZOC 14IPTCU/8DPx ...... i Tail Pressure Belt Scannivalve 0-0.72 psi 2.2e-5 psi ZOC 14IPTCU/32Px

Table A-4. Instrumentation specifications.

NASA/TP--2000-209908 96 L_,,_,-- _...... : • - "T

! • , J

Figure A-1. NASA Glenn Research Center Icing Research Aircraft.

NASA/TP--2000-209908 97 .__------inter-cycle IRT Shape 0 • V=-135 kts, alpha=-2.9 _ • LWC=O.§g/m 3, _D=20_ _'-"---...... To=-4° C ...... --. • time=15 min, with boot ------cycle every 3 minutes

Figure A-2. Inter-Cycle Ice on tailplane.

_-::-:--:--::-:_

._"- Failed Boot IRT Shape ,_ • V=13fi kts, alpha=-2_q °

_;i_._...___ • LWC=0.Sg/rr,3, MVD=20_n \ • T0=-4 ° C, time=22 min

Figure A-3. Failed Boot Ice on tailplane.

NASA/TP--2000-209908 98 derived from in-flight photos _ and ADS-4 "_;_ • used in previous stability l ...... control flight tests

Figure A-4. S&C Ice on tailplane.

Figure A-5. Tail flow probes and pressure belt.

NASA/TP--2000-209908 99 L_.P "7

i_ plo Tin., 14_lp ,.--'.,_ _.-' //_//.2_,,, ,//I_,,El...... ,,,,_ _-

A_

/ ,X-

NOTES ,,..,_,B- B

1. l_r_-_ _ 0_1_1_ I_ '

_,_,,_.,_1,__, -/ ",,.-_ IlllfII If I ",,=- B i B a_.__U..a

Figure A-6. 5-hole tail flow probe head.

i/ .

\\

Lo'_Aer T allIm_gir_ Syslern

Figure A-7. Tail camera position.

NASA/TP--2000-209908 100 Appendix B. Select Steady State Data Plots

_II_.AIC =lphi _ Inflight Conllilten©y Check

Tel[ lip hi. 1 [ (ch01 =CT 14, c h02=CT C7)

.-._-.. Tel elphL 2 ...... -'j-'sHe'p'=3["--'"[°--"-F-'_--]'_".... 'T. [ _1_ _ . [ _ _I_ i

'°°_ - - -_ ...... ----':-''"-: "_'T _ ....:T-: =".' iiiiiiiii!iii!l!l!l!I!l!l!l!ill

Flight Num bet, till point

I_:L°.:,. '°",'o'to'L:::'.Fo"E_,?,;'°k

5O".--41--TllbeteO--'_-TLir betl 32 _'--_ ...... _ .... I" .... "I _ "'T" _ _ _T.__!

00 :021...... L_LI-_ -LL_-_LL_LJ ===,======,=;_[}li[l,,o,, ,°

Flloht Num bet, lest point

-_,_ Ap'C VTAS Inflifiht Con=i=tency Check

[ --_r--Prbl VTAS I (ch(_I=CT 14= ch02=CT 07} +_'"'t-IT-_lr_ Prb2 VTAS T-I _ [1...... 1_;

iii!iiiii!lll!i!l!!!i!l!l!l!ill

Flight Num bet, feet point

---_-.. Cltsil I Inflight Conllilltency Chick ---_- Xc g I (oh0 =CT 14 ch02=CT 07)

O

02 _- ...... +...... +.--_:F:::-.,..--,--.--.--..-,....,-:_---_2,,-::¢_-+ --4,-_

01

CO

-01

-02

-03

-04

-06 -i!iiiiiii liliIllillllllliiii:._ _ . _ . _ . _ . _ . =

--_-.-FIp ._gle Inflight Conoistency Check

I _/_:/El_ph, {ch0I=CT ..... 02=CT 07)

' t I *---*-_*``-_*--*-_---1_-_-_-_-*-_---*_-*___--*_*-_--_`_ 11- ...... - £ " , ...... --_...... 2 I :11__ - ..... :" iilillIlll liililltllililttl _ Figure B-1. 1995 data consistency results.

NASA/TP--2000-209908 101 --_-,A/C alphl i Inflight ConkiClency Check --/-- Tii; alpha t f {Oh01 =CT 10 ch02=CT 04)

.,or-'-*_,_ \...... :...... I .I [ ...... T--T- I-- ..... !I,;-}-_.

-80 .,oo I_.__.I._._L o t...... L_L__J___ ...... L.J i

Flight Numbir, teal point

.. A/C beta I Inflight Conliltency Check _--ToII bitll 1 1 {ch01=CT 10 ch02=CT (}41 --i---Tei! bitl 2

'o_ '

!___ii '_°___!!__,___I E

Flag ht Num b,r, tilt palnt

....@--A/CVTAS t Inflight Coriblitency Check Prbl VTAS --IP-- P'b2 VTA$ t (ch01=CT 10 ch02=CT 04}

--4--- _'_b3 VTAS ...... , : ',

i ii!i!i ii i!i

Flight Numblr, lilt point

,4t--.C.I till I lnflillht Coneiltcncy Check

.... _-,--- CTI l o.oJ . 711 _f_ ; ;-T--CTTT_-/]/--;TT/]T-7 7-TT[-i

-o 100 t, _ ...... _ . - - i

Flight Hum bi r, tail point

.... }'' F_ il_all ] tnfliliht Cllrllilti!il_¢y Check ---i-- dll EIIv _ I T 4 --_iA/CIIphl I (ch01=v, 0, ch02=C 0 ) --I--- Pitch A nglI ! ......

i_ - -= iim;t .l- -i.-: .]=- -t:.{ r i-..i._ - .-i.--_-i-t:_-4_ i i..+i.,, t;.4ici..i...#.. _ _- -i- ii.i-_ ; :::::::: --::::::-:::::::::: :::::::::: : ::;::::i:::::: ; -_"-'-:._-__ .__*_._.__ -_: ° _L _ , - , . _.._ ..._4. .-- _ ._..-_..,,._- --_ ' -L.-,,. _ ' "4

-4

- B .... "......

Figure B-2. 1997 data consistency results.

NASA/TP--2000-209908 102 o

"o 'o "o

= i _>" . o

c _ c'7 J 0 > - m

uJ 3 - _, _8 .... 1 _J _--" = _ i: I=n._"__ _ i____,/"' " ° "_ !11 jit

i

(Gap) _ P (_,p)"" _1

NASAfI'P--2000-209908 103

Appendix C. Time Histories of Select Dynamic Maneuvers

This appendix contains the time histories of 20 quantities for select test points. Two test points for each of the ice contamination cases are presented for both the pushover and elevator doublet maneuvers. For the thrust transition, only the extreme Baseline and Failed Boot cases are shown. These points were chosen to allow comparisons across several different parameters, including head-to-head comparisons of the pushover and elevator doublet. Note that there are two pages per figure. Also note that the velocities from the three tailplane probes (Vprobe) are true air speeds measured in knots, regardless of how indicated.

Also provided is a parameter space finder box for each ice contamination and maneuver case. As depicted in the figure on the right, the axes of this space 1.6 critical are the flap deflection and the speed. For each maneuver, the most critical region corner of this space is identified with a gray box. The dots mark where the V/Vs Most two test points reside. As the flyable parameter space shrinks, so does the 1.0 parameter box. The particular example depicted is of the Baseliae pushover. 0 [iF 40 The most critical region is at higher flap deflections and lower speed. The two test points indicated occurred at ([iF, V/Vs) = (0, 1.5) and (40, 1.0).

List of Figures Figure C-1. Pushover, Baseline, [iF = 0°, VIAS = 1.50Vs = 100kts. R Figure C-2. Pushover, Baseline, [iF -- 40 °, VIAS = 1.06Vs = 55kts. Figure C-3. Pushover, Inter-Cycle Ice, [iF = 0 °, VIAS = 1.50Vs = 100kts. Figure C-4. Pushover, Inter-Cycle Ice, [iF 40 °, VIAS = 1.06Vs = 55kts. U Figure C-5. Pushover, Failed Boot, [iF = 20 °, VIAS = 1.33Vs = 75kts. ["j Figure C-6. Pushover, Failed Boot, [iF = 20 °, VIAS = 0.97Vs = 55kts. Figure C-7. Pushover, S&C, [iF = 0 °, VIAS = 1.5Vs = 100kts. k-'] Figure C-8. Pushover, S&C, [iF = 20 °, VIAS = 0.97Vs = 55kts. U

Figure C-09. Elevator doublet Baseline, [iF = 0°, VIAS = 1.50Vs = 100kts. Figure C- 10. Elevator doublet Baseline, [iF = 40 °, VIAS = 1.06Vs = 55kts. Figure C- 11. Elevator doublet Inter-Cycle Ice, [iF = 30 °, VIAS = 1.59Vs = 85kts. Figure C- 12. Elevator doublet Inter-Cycle Ice, [iF = 40 °, VIAS = 1.63Vs = 85kts. Figure C-13. Elevator doublet Failed Boot, [iF = 30 °, VIAS = 1.59Vs = 85kts. 2 Figure C-14. Elevator doublet Failed Boot, [iF = 30 °, VIAS = 1.03Vs = 55kts. Figure C-15. Elevator doublet S&C, [iF = 20 °, VIAS = 0.97Vs = 55kts. Figure C-16. Elevator doublet S&C, [iF = 30 °, VIAS = 1.40Vs = 75kts.

Figure C-17. Thrust transition, Baseline, [iF = 0 °, VIAS = 1.20Vs = 80kts. Figure C-18. Thrust transition, Baseline, [iF = 40 °, VIAS = 1.63Vs = 85kts. Figure C-19. Thrust transition, Failed Boot, [iF = 0°, VIAS = 1.20Vs = 80kts. Figure C-20. Thrust transition, Failed Boot, [iF = 40 °, VIAS = 1.63Vs = 85kts. D

NASA/TP--2000-209908 105 1o

5 j ...... : . _ _ ' .--_°'E_°°_li

0

-5

-I0

2.5 2 ...... i_ i ._ rl__N, cGII 1.5 I 0,5 0 -0,5 r

t i (deg/s)j

B 7 41 ...... i...... ?

.2 i i_- _ -4 (deg) ! -61[ .... -!----_---I_L_------'---I---;-- , .... I oA_ I -8

8

4 2 0

-2

-6 _

4O 1 ' _ ' _theta (deg)]

-20o - _ -40 --

160 --VIAS (kts) I I

120

1,o100 :i

80 '

3000

280O (ft)] 2600

2400

2200

2000 0 10 20 30 40 50 60 70 time 80

Figure C-la. Pushover, Baseline, 6F = 0% VIAS = 1.50Vs = 100kts.

NASA/TP--2000-209908 106 40 1 (Ibs) i 0

-40

-80

t %,_ '_ ._d" •...... TAOA2 (deg) -10 t ...... _'"_ _ , ...... TAOA3 (deg) !j

0.03 ...... ,- -

-0.03 o - - - _c. -0.06

1 _L ...... -- CLac o.e ...... '...... '--- I 0.6 0.4 0.2°l -0,2 - ' ....

-0.1

-0.2

-0.3

0.1

0 05

0

-0.05

-0.1

160

t40 • . ' !_, • 120 --Vprobel-- " - _ .... _? .... _ :' " '- _ (m/s)[ 100 ...... _.' _ _ ...... i'_ ,_ __ ...... _-_ _V.._" Vprobe2 (mis) f

8O

, , _--Pbetal ()].deg., 4 _ ...... , , / .... Pbeta2 (deg)| . ._._o _ _ .....,¢_. _ L_I_,s_ _ . _ ---- (deg)j 0 t '_ "_"_'_ _ ' ' " " 1

-6 ......

0 10 20 30 40 50 60 70 lime 80

Figure C-lb. Pushover, Baseline, _F = 0 °, VIAS = 1.50Vs = 100kts.

NASA/TP--2000-209908 107 157

lool5I_ _,E(deQ) -5 4 ...... -10 J

2

0.75 0.5 .... 0.25 ...... ' o i

15 10

-5 -10 -15 (deg/s) -20

-6 , I , , , -B -10

-14 0 1 (deg) -16

B

0 ......

-4 -6 -8 - - -10

40 20

0

-20 lh!--,,_ta(deg) -40 -60

B5 75

65 ...... - ' S 55 (kts) 45 35 25

530O

4800

4300

3800 10 20 30 40 50 60 70 50 time 90

Figure C-2a. Pushover, Baseline, _F = 40 °, VIAS = 1.06Vs = 55kts.

NASA/'rP--2000-209908 108 -5

-10

-15 ...... " ...... I ...... TAOA2 (deg) -20 L TAOA3 (deg)

0.06 o.o,i ...... :_V--c_o 0.02o -0,02

-0.04

-1.5 J -

0.4

Cm

-0.2 - -, -0.4 '

120 --Vprobel (m/s) / Vprobe2 (m/s) 100 -[ ...... Vprobe3 (m/s)

.o] .... -.,i --_j" "_.,,_:...... ".,.,j......

204oI ...... •...... -- ......

I__ Pbetal (deg) 15 1 I ...... Pbeta2 (deg) I0 t ...... I Pbets3 (deg)

-5o__l._SJ...... /, -10 10 20 30 40 50 60 70 80 time g0

Figure C-2b. Pushover, Baseline, _F = 40 °, VIAS = 1.06Vs = 55kts.

NASA/TP--2000-209908 109 lO (deg) 5 I ...... _ delE

0

-5

-10

2,5 2 Nz 1.5 1 0.5 o -0.5

15 10 (deg/s) 5 0 -5 -10 -15

6 ...... - _ _ A I--TAOA 1 4 (deg) 2 0 -2 -4 -6

6 6 ...... _ ...... _._ .... _ ...... j--Alpha (deg) 4 2 0 -2 -4 -6

40 _heta (deg) 20

0

-20

-4o

15o

130

110

90 V!AS (kts)

7O

6500 (ft) 5500111111160001...... : ...... I___ALT 5000 ! ......

4500 " 0 10 20 30 40 50 60 time 70

Figure C-3a. Pushover, Inter-Cycle Ice, _F = 0°, VIAS = 1.50Vs = 100kts.

NASA/TP--2000-209908 110 5O

-50o -100 --

lo,'t ::2: :-_ ....._J ...._ (°_°"d°°"°°°' -10

0.02

0 L ,__ //_ -_ _k'_ Che

-0.02

-0.04

1.5

1

0.5

0

-0.5

0,1

-0.1

-0.2

-0.3

0.1

-0.05

-0.1

6

j _'W _ _ _ ..... T _ _'] _ _ Y _ _ _ ...... Pbeta2 -2 j , / Pbeta3 {d:_] -4 0 10 20 30 40 50 60 time 70

Figure C-3b. Pushover, Inter-Cycle Ice, _F = 0°, VIAS = 1.50Vs = 100kts.

NASA/TP--2000-209908 l I 1 10 _--

-15

1.5 -- Nz

1 _ ..... 0.5

10 5 I_ ...... A ...... {_q (d,g/,)

-5 -10 -15

-6

-8

-I0

-12

-I4

4 (deg) ! J

-2 -4 -6 -8

I0 -- theta (deg) 0 -10

-20

-30 55.575 _v.s.., 45

5700

5600 -I ...... [--ALT (ft) 5500

5300 5200 5100 0 10 20 30 40 50 60 time 70

Figure C-4a. Pushover, [nter-Cyc|e Ice, _F = 40 °, _v_[AS = 1.06Vs = 55kts.

NASA/TP--2000-209908 112 -4 ] _--] _% _TAOA1 (dug) -6 / ...... _..L--J " "_._. .... I ...... TAOA2 (dug) / ...... _ __ _1_ :: _ TAOA3 (dug)

-14

0.1 Che 0,05

0

-0.05

p Cltal " -0.5

-1

-I ,5

S5 ....

65 _ "+-_'-'>'- ...... "_ .. -'_._._,;,_- -_> - m--_ Vprobel (m/s) / _"_: ...... Vprobe2 (m/s) 55 U ...... "...... Vprobe3 (m/s) 45 |

10

5

0 " ,-.. '-_-- , .-. -_ "..... W_ _ _'_ -_ _< %_d _

-5 J ..... __/ ...... Pbeta.2 (dug)

-10 I ...... Pbeta3 (dug) 0 10 20 30 40 50 60 time 70

Figure C-4b. Pushover, Inter-Cycle Ice, 5F = 40 °, VIAS = 1.06Vs = 55kts.

NASA/FP--2000-209908 113 15 t0 ...... ' delE (deg)

0 -5 -10 -15

1.75 , 1,5 Nz 1.25 1 0.75 0,5 0.25 0

2O

10 q (deg/s)

0

-10

-20

-4 (deg)

-8 -12 !!o!1 -16

8

4 I ha (deg) 0

-4

-6 - -

-12

20

-10 -20 -30

95 " VIAS (kts) 85

75

65

55

7300 7200 (ft) 7100

7000

6900 6800 10 20 30 40 50 60 70 80 time g0

Figure C-Sa. Pushover, Failed Boot, _F = 20 °, VIAS = 1.33Vs = 75kts.

NASAfTP--2000-209908 114 ...... 20 _1o -40 -10

0

-1o 2"_,...... "" _ - - _- _ - "_-_ ' _ - __._ _. _-_- - - _ _ ...... l %_ Y %_,,_", I--TAOA1

0.15 .... •

0,05

0

-0.05 ......

2.5 2 Clac ] 1.5

1 0.5 iiiiii 0

-0.2 ...... _ ...... " -0.4 o,-0.6 iii;L°'"' -1 ---

0.3

Om 0.1 0 -0.1 -02 4 ...... I -o.3 J

105 ] _ _ ...... A ' Vprobel (m/s)

95 ] ..'_.,;'_ "_ ..... _ff';._.:'_'_ _ _ - " _, " - - - ...... Vprobe3 (m/s) 85 ......

6s | Z

lO ...... _1 -- Pbetal ((::leg) , / _" Pbeta2 (deg) 5 ._ - - _"_'. _._ _ _ ...... Pbeta3 (deg) __ _'_ _'_- _ ,_ ,_._ ,1_, _ _,_ -- 0 -5 ..... - ---V 121 -10 0 10 20 30 40 50 60 70 80 lime 90

Figure C-Sb. Pushover, Fai]ed Boot, 5F = 20% VIAS = 1.33Vs = ?Skts.

NASA/TP--2000-209908 115 15 10 5 0 -5 -10 L- ...... _...... :jiiiiiiiill -15

1.75 1.5 (g) 1.25 1 0.75 0,5 0.25 0

2O

10 _ q (deg/s) I ...... ! ...... 0

-10

-20

0 1 (deg) -4

-8

-12

-16

, i ..... : .... (deg]

-4 -8 -12

20 10 (deg) 0

-10 -20 -30

85 (kts) 75

65

55

45

6700 6600 6500

6400

6300 6200 0 I0 20 30 40 50 60 time 70

Figure C-6a. Pushover, Failed Boot, 5F = 20°, VIAS = 0.97Vs = 55kts.

NASA/TP--2000-209908 116 50 ] ill ^ -_ I _FYE(Ibs) 14 _^ .... ![1_I ...... ,o,E _ o1,.,_=" o / -25 i - , ._ ...... 7 -50 ...._ V . 14

5 I'_ i I--TAOA1 (deg)t 0 • _-- -L._J ' I...... TAOA2 (deg)'

-5

-10

-15

0.1

0.05 CHe ]

0

-0.05 1

-0.1 J ...... _....

2.5

1.5 '_1 iiii °"°_J 0.5 _ ...... ol

0

-O.2 ; __ CLIaItj -0.4

-0.6 -0.8

0,3 -

1 -o.3 '

95 ] ...... ' _Vprobel (kts)] / _ .,A_=_,. ,_.),_i ' "-_ I ...... Vprobe2 (kts) ,s7 -- _: ..... ff _ -- -.1%_ : Cm_l ...... v,,o_o,(,,,_]

55

lO _ , _Pbetal (deg) . , I _. Pbeta2 (deg) 5 .... ._,*_-_": - _;< ..._'_. _- - ._ "_'- ..... __,_,,_._ : ..... _._.... "_';. " " - -_,,,/ i ...... Pbeta3 (deg) 0 ,_,_#*_ " ____t % __ _ _,"t:,_ ..... _ I __

o5

-10 0 10 20 30 40 50 60 time 70

Figure C-6b. Pushover, Failed Boot, _F = 20 °, VIAS = 0.97Vs = 55kts.

NASA/TP--2000-209908 l 17 lO _delE (deg) 5

0

-5

-10

25 2 1.5 1 0,5 0 -0.5

15 10 (deg/s) 5 o -5 -10 -15 -20

(cieg)

6 (deg)

0 -2 -4 i'i -6

30 {deg}

-10 -20 -30

135

125

115

106

95

6400 6200

6000

5800

5600 5400 10 20 30 40 50 tim • 60

Figure C-7a. Pushover, S&C, _F = 0 °, VIAS = 1.50Vs = 100kts.

NASA/TP--2000-209908 118 40 I _ FYE (Ibs)

0

-40

-80

5 -I ...... _ ...... ]-- TAOA2

0 ..... , .

-10

0.02 --Che

-0.02

-0.04 _

1,5 i ...... ' r _ Clac

0.5

0

-0.5

0.1 0 -0.1

-0.2 -0.3 -0.4

0,1

Cm

-0.05o.o5o -0.1

150 _ ._ _ I_ _ f_ --Vprob-el (m/s) I40 I ...... _-"% - ;. - - ,;_,,_'%," = " 3t" [ ,__ -1.-£ - /'1 _" ; ...... Vprobe2 (m/s) J - _-/- - _ ,._.._1 _ _',_J"" _#_,P_l_i IV-I,, L._,, •...... vprobe3 (m/s) 13oI ___.._,'_,.,,.,"% ' I- _,_,",'_,,_ t _ t T E,_.... £ _ .... "_. I _ _._"__ T_.!...... <7.... '_'.. ..___"

6 .... _--- Pbetal {deg) _. / ...... Pbeta2 (deg) 4 ...... _, ...... , ...... 1 _- - - _, - - ,_ _ = _, _ _e_"_,_ , _ I ...... Pbeta3 (deg) 2 _.- ..... _-.=_-_- - - -f_#,_ =-_._ -!.._::,i_i_I_I _- i-; ;_,_ I_ 0 _, . .,," _ . _ _:Wr'?: _. : __ -2 -4 10 20 30 40 50 tim • 60

Figure C-'Tb. Pushover, $&C, _F = 0°, ¥]_$ = l.$0Vs = 100kts.

NASA/TP--2000-209908 119 lO

5 --detE (deg) o

-5 -IO -I5

2 1.75 1.5 1.25 1 0.75 0.5 0.25 0

15 (deg/s)

-10 -15

-2 -4 (deg) -6 -8 -10 -12 -14

6 4 2 0 -2 -4 -6 ...... i -6 -10

20 lO (deg) o -10 -20 -30 -40

95 65 ...... ' ...... _ S (kts)

65 5575 [VIA 45

560O (ft) l I

500051005200530054005500 ...... ' ...... i i ! ! ! ALT

0 10 20 30 40 50 50 time 70

Figure C-8a. Pushover, S&C, 5F = 20 °, VIAS = 0.97Vs = 55kts.

NASA/TP--2000-209908 120 75 18 5O i _FYE (Ibs) ...... delE (deg) 25 6

0 -120 -25 -6 -50 J ' _ _':_ -12

0 -2 L .. _ '_"_ r_ I--TAOA1 (deg)l _, _'_ - - - _'%_;-':...... ,_ _...... '_.b_ _ I ...... TAOA2 (deg) -4 "_,- ...... _,,.:_,_ - - _.',,_,_ - _-_._. "---'- i_::_£_ - •...... TAOA3 (deg) -6 -8 -10 -12 -14 ..... , k

0.1 Che] 0.05

o

-0.05

2.5 k r 2 1.5

1 0.5 iiiii -cac 0

0 -0.2 .... : ...... : .... C,tei, -0.4 -0.6 -0.8 -1

0.2

-OA

-0.2

95 _ ..... ,,,_ , _,_ [_Vprobel (m/s)? 85 _ .... _ _ _._*°"_. _ , . _,s'-.,, - /[_--Vpr°be2 (m/s)l! : - _ • ": ----- Vprobe3

55 ] ...... ', _"...... I_1 45 ' '

lO ...... _...... i I - PbJal(deg) I 5 _:_%_' ...... _ ...... _ ...... _ ...... _ . ,_%_.. J ...... Pbeta2 (deg)

0

-5

-10 0 10 20 30 40 50 60 time 70

Figure C-Sb. Pushover, S&C, 5F = 20 °, VIAS = 0.97Vs = 55kts.

NASA/TP--2000-209908 121 5 --delE (deg)

-5°I -I0

2

0,5 ...... ! oT

(deg/s) 5 ...... I

0

°6 _ i .- ......

-10

024__TAOA1 (deg)

-2 ] ...... -4 _

8

6

4

2

0

10 6 (deg) 6

4 2 0

110

105

100

95

90

5700 5680 5660

5640

5620 5600 0 2 4 6 6 10 I2 time 14

Figure C-9a. Elevator doublet, Baseline, 5F = 0°, VIAS = 1.50Vs = 100kts.

NASA/TP--2000-209908 122 4O I -- FY E (tbs)

-20 ] . - -4(?

_ --TAOAI (deg) J _ _ _ _ I...... TAOA2 (deg)

-2 -" ...... _ *'_._.-_.,jf .... t-'_- " r<..<*'_.>'r` '_.... _ ..... I z I -4

0.02

0,01

0

-0.01

-0.02

1.5 -- Clac 1

0.5

0

0.2

o

-o.2

-0,4

-o.6

02

o

-o.2 ...... ,cm

-0.4

140 --Vprobel (m/s) .>._-_'_. _ _<_ ...... Vprobe2 (m/s) 120 %.._=. _ 7%. _<_+._'<'-'_"_ ...... Vprobe3 (m/s)

100

80

2 , J--Pbetal (deg) 1 _ , I ...... Pbeta2 (deg) 0 4e____.+_._>_,_. it _ __ ...... Pbeta3 (deg) t ;_;_l_ =l,l. ' _,_._._# .,._ -tw .... c# _ _ __

-6 0 2 4 8 8 10 12 time 14

Fi_Llre C-gb. Elevator doublet_ ]]asellrie_ _F = 0°_ ViAS = 1.50Vs = 100kts.

NASA/TP--2000-209908 123 5

0

-5

-10

-15

2

1'51 ...... _- - - - - i i i ...... i \ i ...... [ --Nz

0

10

(deg/s)

0

-5

-t0

-8

-10

-12

-14

-16

8

6 ...... i ...... [-- Alpha (deg)

0 . -2

-4 ......

5 theta (deg) 0

-10

60 --

5658 VIAS (kts) 62 60

6740

6720 ...... _[_A LT (ft) 6700 6680

6660 ...... 6640 0 2 4 6 8 10 12 tim e 14

Figure C-10a. Elevator doublet, Baseline, _F = 40 °, VTAS = 1.06Vs = 55kts.

NASA/TP--2000-209908 124 ' I--TAOA

-4 1 .... _; - - _- ...... ; j';', ...... I...... TAOA3 (deg)(deg)

-16

0.1

0.05

0

-0.05

-0.1

...... :i , mCaoi

0,4 -- Cltail

-0.4

-0.B

-1.2

0.6 0.4 0.2 0 -0.2 -0.4 -0.6

9O 1 --Vprobel (m/s) I / ...... Vprobe2 (m/s) B0 1 ...... -...... Vprobe3 (m/s) _oE__'_'*/_-"7.'__"_'_'_"*_"_.-._y'-_¢*'''...._"-'_-" ....

60 ...... :"%"","_ -'_-...... _oI

8 Pbetal [deg) _._,. Pbeta2 [deg) 6 4 ...... _ ...... _...... Pbeta3 2 0 -2 -4 0 2 4 6 B 10 12 time 14

Figure C-10b. Elevator doublet, Baseline, 8F = 40°, VIAS = 1.06Vs = 55kts.

NASA/TP--2000-209908 125 (deg)

_Nz

0.5 ......

0

(deg/s) ! I

° T -441J ...... -8 !

-6

-8 ' ' t_TAOA1 (deg)

-t0

-12

-14

o t -2 ' ' Alpha (deg) -4

-6 -8 -10

-4 I _ , .... i ...... --theta (deg) ] -8

-12

-16

95

7 I | , , , VIAS (kts) i

85

80 4 ...... l 75 +1 , : ,

7400 7380 ...... ',_ i ..... _[___Atr (ft} I

7360 ii

7320 _, _ 2 4 6 8 10 12 time 14

Figure C-11a. Elevator doublet, Inter-Cycle Ice, 5F = 30 °, VIAS = 1.59Vs = 85kts.

NASA/TP--2000-209908 126 40 ] '_,..--.... ir,,.,___..,_ , • I0

-20 -- -5

I .'_. .._ ,' /.,-,:_'...... _ ...... J...... TAOA3 (deg)[

0.03 ...... • ......

0.02 I Oh _

oo,0 -O.Ol -0.02

I -- Clac ] I

0.5 ' 11 " - 0 o.°.°'°"o.1 ]- _c., -0,9

0.2

0.1 om]

0

-0,1

-0,2

I00 98 96

94 92 90

• , , . _ _ Pbetal (deg) !

86 _ _.,.,_ _. - ' ..... ' ,.=_"_-_I_l _'...... Pbeta2Pbeta3 (deg)[

0 2 4 6 8 I0 12 time 14

Figure C-11b. Elevator doublet, Inter-Cycle Ice, 5F = 30 °, VIAS = 1.59Vs = 85kts.

NASA/TP--2000-209908 127 lO 8t i...... _ ...... _oe,E(deg)

0

2

0.5 ± ...... ' ...... o{ I

12 1 ,_' , ] --q (deg/,) 6 J .... ,_="_._ .... -._,_ - - ( _;..... _ - - - : - - - ] ...... delE(deg)

-4 I

-6

(deg) I -I0 . , __!AoAI -12

-14

-2

-4 ' _Alpha (deg) ,° iiii I -10

o_ _theta (deg) ]

-12

-16

90

I _VIAS (kts) J!

80 '!iiiiiii ii - -

75 _ :

7040

7000

6980

6960 ......

6940 , : : 2 4 6 6 I0 12 tim e 14

Figure C-12a. Elevator doublet, Inter-Cycle Ice, _F = 40°, VIAS = 1.63Vs = 85kts.

NASA/TP--2000-209908 128 60 12

-15 - -- -3

-6 I ' _TAOAI (deg)[ -8 ' - ...... _o_ ...... ' ...... TAOA2 ...._ _'_\ ' ' I...... TAOA3 (deg) -10 (deg) -12 -14

-16

0.06 ...... Che ]

-0.02 ......

2 = ' _ Clac' ! 1.5

1 ___ - J

0.5

0 ; i

Cltail -0.2 ......

-0.8

0,2 . _Cm

-0.1

-0.2

105 _ ...... _ --- _Vprobel (kts)[ | ' ' ------Vprobe2 (kts) 100 _=.._.','_,_:_._.w_. _ ...... _ w_.^_^,_ (kts) J

85

_. / -- Pbetal (deg) 6 _ _._ , _ .... Pbeta2 (deg) ....

0 2 4 6 8 10 12 time I4

Figure C-12b. Elevator doublet, Inter-Cycle Ice, 5F = 40 °, VIAS = 1.63Vs = 85kts.

NASA/TP--2000-209908 129 lO

5

o _ _ _ _ _delE(deg)

-5 L_

2

15

1

0,5

0

12 8 (deg/s) 4

0 -4

-6

-6 1 (deg) -6

-10

-12

-14

0

-2

-4

-6

-8

0

-2

-4

-6 _ _thela (deg)

-8

9O

85 •

8O 75 v,;s kts) 70

6050

......

6010

5990

59706030 t ALT (ft) 5950 0 2 4 6 8 10 12 time 14

Figure C-13a. Elevator doublet, Failed Boot, _F = 30 °, VIAS = 1.59Vs = 85kts.

NASA/TP--2000-209908 130 60i A : i - m-_-_o_-,)-- 12 40 i :_'--_-;-_-,I IN_ ('_.._-...."_ ...... delE (deg) 8 2o,-:-----_.....F',. _ 1 "_-i--- -...... :"--J"-...... ' ° i -20 - ' _ -- -4

-6 _ , , --TAOA 1 (deg) _ -8 .... _, ...... _._ - ...... TAOA2 (deg)i

_,._,_ ,.,_ ...... TAOA3 (deg)! i -10

-12

-14

0.06

0.04

0.02

0

-0,02

2 Clac 1.5

1

0.5

0

0

-0.2 Cltait]

-0.4

-0.6

-0,8

0.2

0.1

0

-0,1

-0.2

95 93 _,,..= , , , _ ...... --_el (kts)] _r_¢_l_._?'_l:,_ ...... , ...... ; ...... Vprobe2 (kts) = 91 " _'_'_:'_ ...... Vprobe3 (kts) ij 89 87 85

8 " - _ _ Pbetal (deg) 6 [ _ - - - . _._ _._._._ _.,_ _, ;,,._,. -._._.._;._.._ - J ...... Pbeta2 (deg)', 4 . :_,'_._'._._".'_".'_._'"=';_'_:'_'_"___'_"_...... _'_ _'_'_k._ ' ,,_,__ ' ...... _ | ...... Pbeta3(deg) i 2 0 -2 -4 0 2 4 6 8 10 12 lime 14

_igure C-13b. Elevator doublet, Failed Boot, 6F = 30 °, ¥]AS = 1.59_s = 85kts.

NASA/TP--2000-209908 131 5 (de9) 0

-5

-10

-15

14

1,2 mNz

I

0.8

0,6

8

4 (deg/s) 0 -4 -8 -12

-2

-4 '. - - i ...... '_ _ _ _TAOAI (deg) -6

-8 -I0 -12

12

8

4

0

-4

10

5

0

-5

-10

85L60 I ...... ! : : _s,k,s,

4550| ...... :::i

5500 (ft) 5450

5400

5350 0 2 4 6 8 10 12 14 time 16

Figure C-14a. Elevator doublet, Failed Boot, _F = 30 °, VIAS = 1.03Vs = 55kts.

NASA/TP--2000-209908 132 , --14 40 _ ¢1_ _Jl_, " , -- FYE (;bs) 20 T ...... #lf_,_ |V _ ...... I_.#_,...... delE (deg) - 7

0 _...... 0

-40 "J" " ' " _ -- -14

I0 , r

5 /_,:-';_- ...... _ _'4_- _,...... ,...... _._-. TAOA2TAOAI (deg)I(deg) 0 . .'_'; L .... "*"" _- I ._...... TAOA3 (deg) _ 5

-15 ] 2 k.l

0 -- , 0.05o., __ ...... :-!iiiii -0.05 ....

-0,1 -L . , ,

4 - = 1

3 ...... -- Clac

2

1 ...... ' ...... ' ......

0 ; ;

0,2 _ - _ 0 I i i !

-1

0.4 0.2 0

-0.2 -0.4 -0.6

75 L " ' _ ' _.... --Vprobel (kts) I -., _ ._,_*."4r'___"__ Vprobe2 (kts) 65

55 1 _'_.... , ; 45

15 ] ...... i--Pbetal (deg) I. / e< ' ' -- ...... Pbeta2 (deg)I 10 _ .... ; _-._-_...... _ ...... ,...... Pbeta3 (deg) J 5 ._,;h'_*_ _'__.,.C_ _"_-_ = ._, _&_,_,_,_+,_,-_'_-_'-_ _ ._-'_.,_,_*;_._,_. 13"-F o _ ! q.L_l

0 2 4 6 8 10 12 I4 lime 16

Figure C-14b. Elevator doublet, Failed Boot, 5F = 30 °, VIAS -- 1.03Vs = 55kts.

NASAgFP--2000-209908 133 2 _

1.5 ...... I 1 .... _.- .... i !__Nzla_ 0.5 ...... ' ...... !

lO

(deg/s) 5

0

-5

-10

-2 • . (deg)

-10 -12 -14

(deg)

4 0 ,

is] --theta (deg) 10

5

0

-5

65

60

55

50

45

646O 6440 6420

6400

6380 6360 0 2 4 6 8 10 12 time 14

Figure C-15a. Elevator doublet, S&C, 8F = 20 °, VIAS = 0.97Vs = 55kts.

NASA/TP--2000-209908 134 40 _ A _. i _FYE(Ibs) ----,delE(deg) 16

-4oi \"....." " ', '. -16

lO

5

0 ...... •,_ _ _ i , ! .... TAOA3 (deg)j -5 ._,_,:_ '_ , 12l -10 -15 U

0.1

0,05

0

-o,05

-o.1

4 ' ' -- Clac 3

2

1

0

0.2 0 Cltail -0.2 -0.4 -0.6 -0.8 -1 -1.2

0.4 0.2 0

-0,2 -0.4 -0.6

_ i" _,# , _ vp_uuu, (kts) ] 50 ...... _,_ .... ._,4 ..... _,y...... : ...... Vprobe2 (kts) 40 I _!_I; _'_ ', ', Vprobe3 (kts)

10 T ;_ .i --= _ _--Pbetal (deg) l I ..._._& _,_, ,_%_ ' ',_._,.,...... I.... Pbeta2 (deg) 5 -L ...... _ J" .'_;_: ';: t .__,_..... _.t __% . _,e_, _,__ ,_,_ . _'_L-_=_:_ _:z,_ ...... Pbeta3 1 o

-10 0 2 4 6 8 10 12 time 14

Figure C-15"0. Elevator doublet, $&C, _F = 20 °, YL_S = 0.97_s = 55kts.

NASA/TP--2000-209908 135 15 1 10 (deg) J 5 0 -6 -10 -15

2

1.5

1

0.5

0

15 I0 i 5 ;_, ....."\ ...... IE (deg) T I ' ' '-_ i i Imq (deg,'s) ' 0 -5 ..... t_ _ _ _ _"___, ...... -10

-15 J

-6 -8 -10

-12 -14 -16

4 .... --Alpha (deg] 0

-4

-6

-12

2O --theta (deg) 10

0

-10

-20

85 8O 1 : : ' ' VIAS (kts) ] 75

7O 65 6O

7050

7000

6950

6900

6850 0 2 4 6 8 10 12 14 16 18 time 20

Figure C-16a. Elevator doublet, S&C, _F = 30 °, VIAS = 1.40Vs = 75kts.

NASA/TP--2000-209908 136 9O 12 ; _ _kp\ L i ', [_FYE(Ibs) ...... delE(deg) 45 _..,

-90 , , v" , , -12

0 t .... _ [_TAOA1 (deg) I / .... f__:,_ I_1 ]..... TAOA2 (deg) .5 t ...... _;- - ..... , -_ ...... ,..... ,-/ .... "_ I I -I ...... TAOA3(deg)

0.15 , • • m 0.1 ...... ' _ ' ..... '..... '_ _ I _ChoP 0.05 0

-0.05 -OA

2,5 ...... _ • . • _ .....

1.5

1 0.5 0

I L--Cltail' ' _ _ _ 1

-1 ......

0.4

0.2

0

-0.2

-0.4 "] . . : IT] (k,,)] , . .._.,._._-,_.-_-_,_. , ," ,, _..__'_:: _Vprobe2 (kts) 85 V probe3 (kts)

65 _ .... _"

10 I .... ---I_Pbetal (deg) I ' ' _ ' ' J ...... Pbeta2 (deg) 5 t _,- % ; ..*.R. • ..... _-_,_._., I ...... beta3 (deg)

0 2 4 6 8 10 12 14 16 18 tim_ 20

Figure C-16b. Elevator doublet, S&C, 5F = 30 °, VIAS = 1.40Vs = 75kts.

NASA/TP--2000-209908 137 -4

-5

-6

-7

-8

2 1 1.5

1

0.5

0

t,...... _ ...... -- i

o

-I

-2

9 _

--Alpha (deg) I J

6

20

15 J

10

5

0 ,'-

0.3

g0 _ ...... [ --VIAS(kts) Ct

75 0

5000 4g00 4800

4700

4600 4500 10 20 30 40 50 tlme 60

Figure C-17a. Thrust transition, Baseline, _F = 0 °, VIAS = 1.20Vs = 80kts.

NASA/TP--2000-209908 138 5

0 (Ibs)I

-5

-10

-15

3 ...... _ - I--TAOA I (deg)] 2

1 _:' ' :" _':_: '_;_ :;' : . " " _ " ': " ,_;_:_::-:*_J_) ...... TA OA 3 [deg) 0 -1 -2

0.005

o

-0.005 - -

-o.01

2 , -- is t __c,aoI

0.5

0

0 , .

-0.2 Cltaill I -o,_] ',

-0.8"0'6t ...... -I

0.02 _Cm ] 0,01 ...... j oo: ...... !......

°0102 ....

120 ...... _ ......

rT-I ...... Vprobe2 (kts) 100 i ...... Vprobe3 (kts) 90 110 " _' -3_-_Vpr°bel (kts)i 80 7O

10 ] ...... _ ...... _I (deg) p ...... Pbeta2 (deg) 5_ ,---_ .... _ .... lOl_ ...... Pbeta3 (deg) _j_._.:_._ _ ....._._._,_,......

10 20 30 40 50 time 60

Figure C-]To. Thrust transition, Baseline, SF = 0°, VIAS = |.20Vs = 80kts.

NASA/TP--2000-209908 139 " ' --delE (deg) I

1.51 _ ......

0.5 , ...... ' .... : ...... t 0

2 1 1 j - . ,r_ ,,_ - - _ /_ ..... i ...... _------q(deg/s)j

0

-1

-2

-10

-11 ' ' I--TAOA1 (deg) !i -12

-13 -14

-15

-6

-7

-8 ...... !i

-9

0

-5

-10 -15 .....i -20

95 O.3

______,:_.,,,n-,,.,J",_-,_,-_,--,,_....,-.,--,,_h,,_,,,.,.,n.:_l.r_ 9O ...... ,__ ._,_.._,,-..__ .... , ...... _ ...... 0.2 _VIAS (kts) ...... Ct , . 85

, :_ ...... 0.1 8O 0 8600---_-16'°° ...... :.... i ...... I--ALT""]

°zoo6300 _iii 6200 10 20 30 40 50 60 time 70

Figure C-18a. Thrust transition, Basefine, _F = 40 °, VIAS = 1.63Vs = 85kts.

NASA/TP--2000-209908 140 (lbs) ] J

-lO I ...... '...... i...... [--TAOA_(_ogll -11 | " , , | ..... TAOA2 (deg) -12 _ , , [ ...... TAOA3 (deg)

-14 cI ..... "" _ _ , -15

0.01

0.005 0 ....t.oheE

-o.005

1.5 I 1 iiiiiii °'"' 0.5 I

-0,2 ...... , ...... ' ...... Cttail I -0.4

-0.6

-0.8 ...... -1

-0.01

-0.02

110 W _, ..... :._-,,J,,_, _ - -, - - -..=_,_,_,,_,.,_ -...... Vprobe2 (kts)]

70

IO I _ - - - / _ Pbetal (deg)

l ...... Pbeta2 (deg)

-5 0 10 20 30 40 50 60 time 70

Figure C-18b. Thrust transition, Basefine, _F = 40°, VIAS = 1.63Vs = 85kts.

NASA/TP--2000- 209908 141 -5 [deg) -6

-7

-8

-9

2

Nz "1.5

1

0.5

0

2 _q (deg/s) 1

0

-t

-2

_ , !AOA1 (deg)

0

10 AIpha 9 '-°'7 8

7

15

10

5

0

-5

90 0,3

85 t _VIAS (kts) ...... Ct 80

75

64O0 ,

6300 +I ...... ' ' [--ALT (ft) li 6200 _ ......

5900 4 ; ; , 0 5 10 15 20 25 30 35 time 40

Figure C-19a. Thrust transition, Failed Boot, 5F = 0_, VIAS 1.20Vs = 80kts.

NASA/TP--2000-209908 142 10 ......

-15 _ , m

6

4 _ -"...... --_...... ", .... •...... TAOA3 3 ,'_' _:_:" _,_':_._'...... _._. _._._!.,:_,_. _:,_:._...... (deg)

0.01

0.005 . . Pche i 0

-0.005

-0.01

2

1.5

1

0.5

0

0.2 Clta il

-0.2 ...... '...... ' ......

-0.4 ...... -0.6 ...... '...... ' ...... -0.8 .... _ • m

0.02 °°_o_,:,:,-,._, :_._, _ _/_V:_:L_:_i_,_,. _J °mi oo,J,"_ V:V_ _ V WITY: VV_ °_T_' -0.02 _ _

120 1 . _ --vp,ob_ (kts)[ 110 i ...... , _ - -.._ , - " --.-.--. Vprobe2 (kts) | _ " , ...... Vprobe3 (kts)

80 _,.,_.,,,_ _',_ _ - ..... _,_,_,_' __'_ ...... __-'_ ...... _"

i F3

...... I ...... Pbeta3 (deg)l

0 5 10 15 20 25 30 35 lime 40

Figure C-19b. Thrust transition, Failed Boot, _F = 0°, VIAS = 1.20Vs = 80kts.

NASA/TP--2000-209908 143 IO - ', : . ! .[_ delE (deg) ! 5 I .... _. _ J 0

-5

-10

2

1.5

1

0.5

0 -' '. : : ...._- ,

20 j ...... : ...... :. I -t_; ..... [_q (deg/s) ] 1°ii122211_:2_222:22 122222'1211 ' 12 --Fi.p_u_2,_ 6_. . . _ __' ,-,.,. ,,,,,_.. _ _ _ v-,.,_,., " -5 -lo° I'...... :'._....._ _'Z_:: ...... ""_ _ _ __"__ _ :.....:_"'_- :', _ i ...... -_- _ -15

-6 ;

-10 -12 I -14 1 (deg) j -16

2

-6 -8 -10 -12 deg) -14 -16

-20 -30o theta (deg) 1 -40

9585 ; -_ ,i _: _ _, _ _ ...... 0,20.3

75

65 0

7400 7 , . 50 7200 ...... ;...... ;...... ', _ .. _._...... 40 7000 =. 2d 6800

6600 ] - _ALT(ft) ..delF(deg) _ _ _ : .... ._ : . - _ - 10 6400 . ; ; 0 0 5 10 I5 20 25 30 36 40 45 time 50

Figure C-20a. Thrust transition, Failed Boot, 5F = 40 °, VIAS = 1.63Vs = 85kts.

NASA/TP--2000-209908 144 20O (rbs)! 150 I 100

50

0

-4

-8

-12 I .... '_ ,d_ --TAOA1 (deg)] -16 ...... •\_'_ .... J.... TAOA2 (dag)] -20 J : : : ; "'_ D _ 1--" rAoA3(_,g_]

0.15 , , , i __c.o 0.1 0.05 i- 0

2

-- CI c ] 1 m 5 J

0.5 J 0 q 2

__i ...... ', ..... ', ...... '..... :_ .. L:_ [k_,,_"" _ -- Cltail ]', -o.4 ...... _ ...... -- -0.6 4_ -_ - - v'_-'_t-.,___ I_.J ------'- - - -0.8 I ,, ,,':-y .1t i

.... _ --cm I o.1 I...... _J_ ..' ^.it^ _ _A_',,. A^'_ ,!

-0.1 [ _ v_j_ _' • -o.2/ ..... :..... : ...... :..... : i v

12o] : ] [ : !r,--I --vp,ob_1 (kt_)' 110 J - - - , - :L.J -_-_,-, Vprobe2 (kts) tf 100 ' ' ' ' _,_ -- Vprobe3 (kts) j

80 70

10 ...... ,,. -- Pbetal (deg) I • , ...... Pbeta2 (deg) I '. '. ' ' .-" ,.,, I---. Pbeta3 (deg)

o r...... _-_'_ - _ _ !, _/! X_ I_,°°*1' -5J...... :___ ,...... ;...... :...... "..._..."" "" ' ' 0 5 10 15 20 25 30 35 40 45 time 50

Figure C-20b. Thrust transition, Failed Boot, 5F = 40°, VIAS = 1.63Vs = 85kts.

NASA/TP--2000-209908 145

Appendix D. PressureCoefficient Distributions of the Full Tail Stall

This appendix contains a detail of select time histories in Figure D 1. The events before the stall to after the recovery are detailed. Indicated on the time histories are the times where the pressure distributions are presented. The stall point is also noted with the vertical line.

Figure D2 depicts the pressure tap locations. In Figure D3, the information presented around the Cp plots includes o_, time stamp, CLr_i1and &E. The suction surface pressures are denoted with the solid, blue line.

10 -4

5 '-S m ,,c LIJ e,l m 0 ? -12 IB ¢D 1D -- del E avg ,2 del E avg -5 I'-- ...... Stall - - -16

¢, TAOA1 -10 -20

27 28 29 30 31 32 33 34 35 36 37

-0.2

-0.3

-0.4 t'_ I'-- -0.5 I -0.6

-0.7

-0.8

27 28 29 30 31 32 33 34 35 36 37

200

I60 m t,n 100 I,LI >- II 50

0

27 26 29 30 31 32 33 34 35 36 37

0.16 40

0.12 37.5

U,. I" 0.08 35 ID "0

0.04 32.5

0 3O

27 28 29 30 31 32 33 34 35 36 37 Tim e (s)

Figure D-1. Select time histories of the stall event. Circles indicate times presented in Figure D-3.

NASA/TP--2000-209908 147 0.2

0.1

0.0

-0.1

-0.2 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2

Figure D-2. Location of pressure taps on clean airfoil.

NASA/TP--2000-209908 148 alpha tail= -11.56 ° alpha tail= -11.91 °

2 2 15 1 I[ G o 0

-I -i

-2 -15 -2 20 AO 60 0 10 20 30 4C C _- -0.4860 tim e=o x (inches) x (inches) Cl= -0.63 3O dE = 7.53 dE= 5,06

alpha tail= -11,71a delF = 39 alpha tail = -11.82 °

1 11

o I ..... !'/ -1 t -2 __t -_/ 0 20 40 60 0 i0 20 3_ _0 50 60 x (inches) CI= -0.56 Cl= -0.57 21.38 30.21 x (inches) dE = 4.94 dE 4.63

alpha tail = -12.59 ° atpha tail = -11.65 °

2 _ ......

I t 1 4,

-1

-1 -_-_7 ,,,_L_

0 20 40 60 o 1o 20 3c (o o 50 x (inches) CI= -0.60 x (inches) CI-_- -056 27.96 30.53 dE= 6.11 dE= 3.78

alpha tail = -12.62" alpha tail= -12.92 = 2

I .....

o

-2 [ ......

o 10 20 30 40 50 60 0 20 20 30 ¢0 50 50 CI= -0.65 X (inches) Ct= -0.62 28.27 31.95 x (inches) dE = 5,98 dE= 6,28

Figure D-3a. Plots of Cp for full tail stall event at various times. Solid, blue line for lower surface.

NASA/TP--2000-209908 149 [ alpha tail = -15.99" alpha tail = - 12.69 ° i 2 ...... 2

1 L

o O-

o I0 20 30 40 50 60 I0 20 30 40 50 60 X _',,-=1"no"--" CI-- -0.60 32.18 x (inches) CI= -0.45 33.43 dE= 4.57 dE= 8.08

alpha tail= -16.07 ° alpha tail = -13,09 °

2 ----_T

2 ...... l ......

: ,%

2+- r I

-2 2 O 10 20 30 40 50 60

x (inches) Ct= -060 0 I0 20 30 40 O _=50 "041 60 32.34 33.72 x (inches)

dE= 5.88 dE= 7_73

alpha tail = -13.27 ° alpha tail = -16.25 ° 2 ......

1 -

b 0 :_. "" -h_E__-_ o _....___.. 1_ I ' lJl ;

o 1o 20 30 40 ciSO 6,_ x (inches) = -0 6 0 20 20 30 40 C=} SO -0.4860 32.49 33.95 x (inches) dE= 5.73 dE= 734

alpha tail= -1573 = alpha tail = -14.89 °

2 I , I

o -_- "---- _,-_ •

-L - . J i -2 ...... J -2 ...... 0 I0 30 40 50 60 0 lo 20 30 4o Cl= '° - 0% 20 x (inches) 33 34.24 x(inches) Cl= -0.47 dE= 8.10 dE= 6.29

Figure D-3b. Plots of Cp for full tag stall event at various times. Solid, blue line for lower surface.

NASA/TP--2000-209908 150 alpha tail= -12.850

2-

1-

0-

-1

-2 20 20 30 _;0 50 60 x(inches) CI= -0.62 34.73 dE= -2.20

alpha tail = -11.79 °

2-] ......

;,... ,:\

-2

20 20 30 _0 __50 60 34.86 x (inches) CI -0.60 dE= -3,99

alpha tail = -9.32 ° iI...... i...... /......

L,/

o :o _o 30 _o Ci5O= 6o x (inches) -0.69 35.33 dE= -256

alpha tail = -5.86"

1 ...... "......

l/ I,

o "-o _o 30 _o 015o= 60 x (inches) -0.45 35,67 dE= -2.85

Figure D-3c. Plots of Cp for full tail stall event at various times. Solid, blue line for lower surface.

NASA/TP--2000-209908 151

Appendix E. Recommendations to the Pilot

Be aware that the tailplane is a more efficient collector of ice than a wing. If ice is detected on the windshield wiper or other objects with sharp leading edges, then it is probably also on the tailplane. The first recommendation is to activate the de-ice/anti-ice protection system. Another recommendation is to manually fly the plane. The cues for an impending horizontal tailplane stall are subtle and much more difficult to detect through the autopilot. Symptoms of contamination include changes in the trim point, difficulty in trimming, onset of PIO, and buffeting in the control yoke as the flow separates and reattaches aft of the hinge point. Also be aware that wing flap extension aggravates a tailplane stall situation. While the differences in tactile cues between a wing stall and a tailplane stall are subtle, the recovery procedure is diametric. When the tail stalls, recovery is achieved be reattaching the flow on the lower sur- face. This is opposite a wing stall, which requires that the flow reattach on the upper surface. To recover the aircraft, pull the nose back, raise flaps and be judicious with thrust. The reduction of thrust played a significant role in the recovery of the (modified) DHC-6. This is believed to be configuration specific; namely, for aircraft with thrust lines above the CG, adding thrust increases the nose-down pitching moment.

In other words:

IF ICE ACCUMULATION ISNOTED ON OBJECTSWITH SHARP LEADLNGEDGES (LE., WINq3SHIELDWIPER) THEN

Recommendations • Activate Ice Protection • Hand Fly Aircraft • Extend Flaps with Extreme Caution

Symptoms of Ice Accumulation on Horizontal Stabilizer (Listed in order of increasing seve_rit)_ If flaps have been fully or partially extended, look for: • Change in Trim Setting from Non-icing Conditions • Difficulty in Trimming • Onset of PIO • Buffet in Control Yoke (not airframe)

Recovery Procedure • Yoke Back • Flaps Up • Thrust Reduce (may be configuration specific)

NASA/TP--2000-209908 153 REPORT DOCUMENTATION PAGE Fo_AVprove_ OMB No. 0704-0188

Public reporting burden for this collection of information is estimated to average t hour per response, including fhe t'ime'ior reviewing instruCtions, searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the coIlection of information Sen0 comments r_garding this burden estimate or any other aspect of this CoI|ection of information, inciuding suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503 1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED March 2000 Technical Paper 4. TITLE AND SUBTITLE 5. FUNDING NUMBERS

NASA/FAA Tailplane Icing Program: Flight Test Program

WU-548-21-23--00 6. AUTHOR(S) Interagency Agreement DTFA03-95-90001 Thomas P. Ratvasky, Judith Foss Van Zante, and Alex Sire

i]7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION REPORT NUMBER National Aeronautics and Space Administration John H. Glenn Research Center at Lewis Field E-12126 Cleveland, Ohio 44135-3191

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

National Aeronautics and Space Administration Washington, DC 20546- 0001 NASA TP--2000-209908 DOT FAA AR-99-85

11. SUPPLEMENTARY NOTES

Thomas P. Ratvasky, NASA Glenn Research Center; Judith Foss Van Zante, Dynacs Engineering Company, Inc., 2001 Aerospace Parkway, Brook Park, Ohio 44142; and Alex Sim, NASA Dryden Flight Research Center, Edwards Air Force Base, California. Responsible person, Thomas P. Ratvasky, organization code 5840, (216) 433-3905.

12a. DISTRIBUTION/AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE

Unclassified - Unlimited Subject Categories: 08, 05, and 02 Distribution: Standard

This publication is available from the NASA Center for AeroSpace Information, (301) 621-0390. 13. ABSTRACT (Maximum 200 words) This report presents results from research flights that explored the characteristics of an ice-contaminated tailplane using various simulated ice shapes attached to the leading edge of the horizontal tailplane. A clean leading edge provided the baseline case, then three ice shapes were flown in order of increasing severity. Flight tests included both steady state and dynamic maneuvers. The steady state points were 1G wings level and steady heading sideslips. The primary dynamic maneuvers were pushovers to various G-levels; elevator doublets; and thrust transitions. These maneuvers were conducted for a full range of flap positions and aircraft angle of attack where possible. The analysis of this data set has clearly demonstrated the detrimental effects of ice contamination on aircraft stability and controllability. Paths to tailplane stall were revealed through parameter isolation and transition studies. These paths are (1) increasing ice shape severity, (2) increasing flap deflection, (3) high or low speeds, depending on whether the aircraft is in a steady state (high speed) or pushover maneuver (low speed), and (4) increasing thrust. The flight research effort was very comprehensive, but did not examine effects of tailplane design and location, or other aircraft geometry configuration effects. However, this effort provided the role of some of the parameters in promoting tailplane stall. The lessons learned will provide guidance to regulatory agencies, aircraft manufacturers, and operators on ice-contaminated tailplane stall in the effort to increase aviation safety and reduce the fatal accident rate. 14. SUBJECT TERMS 15. NUMBER OF PAGES 166 Aircraft icing; Tailplane icing; Stability and control; Aerodynamics; Aeroperformance 16. PRICE CODE A08 17. SECURITY CLASSIFICATION 18. SECU RITY CLASSIFICATION 19. SECURITY CLASSIFICATION 20. LIMITATION OF ABSTRACT OF REPORT OF THIS PAGE OF ABSTRACT Unclassified Unclassified Unclassified

NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) Prescribed by ANSI Std. Z39-18 298-102