[[alternative]]Aircraft aerodynamic analysis of ATR-72 under ice accretion condition

碩士[[abstract]]2002年12月21日，復興航空GE-791貨機，因遭遇積冰問題而失事墜毀 於澎湖外海。本文針對此台灣境內首例因積冰失事之案件，參考飛行安全委 員會之失事報告，並以復興航空ATR-72為例，模擬二維外形下主翼與水平尾 翼遭遇積冰狀態之升阻力係數與力矩係數相關變化，更進一步的利用飛行力 學分析軟體探究在三維外形下積冰現象對飛機性能參數之影響。 本研究以CFD方法模擬積冰外形生成，利用對翼外形產生非結構網格， 將Navier-Stokes方程式簡化為Euler方程式並以有限體積法離散求解之，最 後以水滴軌跡方程式生成積冰之外形。考慮ATR-72機型因特有T型尾翼，水 平尾翼不受主翼之尾流場之影響，針對當時大氣環境分別模擬出光滑冰 (Glaze Ice)與透明冰(Clear Ice)兩種不同積冰於主翼及水平尾翼上作堆 積，並利用FLUENT軟體求其流場解分析其氣動力之影響。並利用飛行力學軟 體DATCOM模擬在三維外形下遭遇積冰現象下飛行性能參數。 根據計算之結果發現，在二維狀態下於翼面上積冰外形越趨不規則，其 升力係數將越趨遞減，阻力係數隨之提高；在水平尾翼上之積冰將減少其負 向升力，導致負向力矩之增加，水平尾翼氣流提前分離造成失速狀態，將使 機首產生一突然向下之俯仰力矩，對飛行安全上造成極大之威脅。在三維外 型外型處理下，因DATCOM為飛機設計之初步設計之工具，對於積冰後外型處 理有較大誤差，因此未來在這一部分有較大的改進空間，以更精確的模擬出 真實的飛行姿態。[[abstract]]The Trans Asia Airway flight GE-791, an ATR-72 airplane encountered ice accretion to crash at Peng Hu Island, Taiwan, December 21, 2002, was the first ice accretion accident case ever in Taiwan. This thesis is based on Aviation Safety Council accident report, and taking GE791 as the case for numerically simulating aerodynamic coefficients of main wing and horizontal tail under several two-dimensional icing conditions. This study uses existing CFD method to predict the ice accretion on airfoil. First step is the generation of unstructured grid on the wing and tail configurations. Secondly, solve Navier-Stokes equation with finite volume method. Finally, the equation of motion for supercooled water droplets is implemented for ice accumulation production on the airfoil leading edge locations. Considering the ATR-72 has distinctive T-tail shape, it is assumed that the tail will not influenced by the main wing trailing vortices, thus the glaze ice and clear ice can be simulated under proper atmospheric conditions. The two-dimensional iced airfoil shapes are achieved, and the results are extended to three-dimensional and then IV use the flight mechanics software DATCOM to analyze the aircraft performance behavior. According to the above research, in two-dimensional the airfoil ice shape is more anomalous; their lift coefficients will hastily decreasing and drag coefficient will increase. In addition, the ice accretion phenomenon will cause the horizontal tail to stall ahead of time. When tail stalls, this downward force is lessened or removed, and nose of the airplane can severely pitch down, this leads to threaten of flight safety and is the major cause of ATR-72 accident. On the other hand the DATCOM is an aircraft preliminary design tool, so it has error in analyzing complicated ice shapes. There are still major rooms for improvement in prediction of actual airplane flight performance through this software.[[tableofcontents]]目錄 第一章 緒論 ............................................. 1 1.1 引言 ............................................... 1 1.2 積冰造成之飛安事件 .................................. 1 1.3 文獻回顧 ............................................ 6 第二章 積冰現象 ........................................ 10 2.1 積冰生成之條件 .................................... 10 2.2 積冰種類 .......................................... 12 2.2.1 透明冰(CLEAR ICE 或 GLAZE ICE) ...................... 12 2.2.2 霜狀冰(RIME ICE) .................................. 13 2.2.3 混合冰(MIXED ICE) ................................. 14 2.3 飛行中積冰現象 .................................... 14 2.4 降落時積冰 ........................................ 16 2.5 預防積冰的產生 .................................... 16 2.5.1 防冰系統(ANTI-ICING) .............................. 17 2.5.2 除冰系統(DE-ICING) ................................ 17 第三章 研究方法 ......................................... 19 3.1 CFD 模擬 .......................................... 19 3.1.1 積冰模型建立 .................................... 20 3.1.2 流場解 .......................................... 26 3.2 飛行力學軟體分析 ................................... 27 第四章 結果與討論 ...................................... 29 4.1 CFD 計算結果 ....................................... 29 4.1.1 積冰外形驗證 .................................... 29 4.1.2 流場解驗證 ...................................... 31 4.1.3 網格點品質分析 ................................... 32 4.1.4 主翼積冰外形建立 ................................. 35 4.1.5 水平尾翼積冰外形生成 ............................. 46 4.2 飛行力學模擬 ...................................... 58 4.3 結果與討論 ........................................ 63 第五章 結論與未來展望 .................................. 66 5.1 結論 .............................................. 66 5.2 未來展望 .......................................... 66 參考文獻 ................................................ 69 圖目錄 圖2-1 CLEAR ICE[23] .................................... 13 圖2-2 GLAZE ICE [23] ................................... 13 圖2-3 RIME ICE [23] .................................... 14 圖2-4 MIXED ICE [23] ................................... 14 圖2-5 飛機力矩關係圖[24]................................ 15 圖2-6 飛機力矩變化關係[24] .............................. 16 圖2-7 氣壓式除冰套除冰方式[25] .......................... 18 圖2-8 氣壓式除冰套[25] ................................. 18 圖3-1 ASPECT RATIO 定義圖 ............................... 21 圖3-2 積冰模型建構流程示意圖 ............................ 26 圖4-1 積冰16.7 分鐘後[23] ............................... 30 圖4-2 各種積冰外形比較[23] .............................. 30 圖4-3 NACA63012 結構性網格全圖 .......................... 31 圖4-4 NACA63012 結構性網格局部放大圖 .................... 32 圖4-5 實驗值與數值解之升力係數關係圖 .................... 32 圖4-6 升力係數與攻角關係圖 .............................. 34 圖4-7 阻力係數與攻角關係圖 .............................. 34 圖4-8 ATR-72 主翼翼剖面外形 [UIUC AIRFOIL DATABASE] ..... 37 圖4-9 乾淨外形初始網格點生成 ............................ 38 圖4-10 乾淨外形局部放大................................. 38 圖4-11 積冰時間12 分鐘，所產生外形 ...................... 39 圖4-12 積冰產生後最終網格圖 ............................. 39 圖4-13 積冰產生後局部放大圖 ............................. 40 圖4-14 間隔4 分鐘之積冰外形生成 ......................... 40 圖4-15 積冰產生後最終網格圖 ............................. 43 圖4-16 積冰產生後局部放大圖 ............................. 44 圖4-17 主翼CLEAR ICE 間隔4 分鐘之積冰外形生成 ........... 44 圖4-18 主翼升力係數與攻角關係圖 ......................... 45 圖4-19 主翼阻力係數對攻角關係圖 ......................... 45 圖4-20 主翼力矩係數對攻角關係圖 ......................... 46 圖4-21 ATR-72 飛機水平尾翼 .............................. 47 圖4-22 NACA63012 [UIUC AIRFOIL DATABASE] ................ 47 圖4-23 乾淨外形水平尾翼翼剖面全圖 ....................... 48 圖4-24 乾淨外形水平尾翼翼剖面局部放大圖 ................. 48 圖4-25 水平尾翼CLEAR ICE 積冰外形 ....................... 49 圖4-26 ICE1：12 分鐘水平尾翼積冰局部放大圖 .............. 49 圖4-27 ICE2：積冰外形上下顛倒局部放大圖 ................. 50 圖4-28 二相流架構下蒐集係數與翼端前緣關係圖[18] ......... 50 圖4-29 ICE3：COLLECTION EFFICIENCY 作用下積冰局部放大圖 51 圖4-30 水平尾翼升力係數與攻角關係圖 ..................... 56 圖4-31 水平尾翼阻力係數與攻角關係圖 ..................... 56 圖4-32 水平尾翼力矩係數與攻角關係圖 ..................... 57 圖4-33 水平尾翼對主翼之力矩係數與攻角關係圖 ............. 57 圖4-34 ATR-72 未積冰三視圖 .............................. 59 圖4-35 ATR-72 未積冰三維外形結構圖 ...................... 59 圖4-36 積冰在主翼剖面生成圖 ............................. 60 圖4-37 三維外形下升力係數與攻角關係圖 ................... 62 圖4-38 三維外形下阻力係數與攻角關係圖 ................... 62 表目錄 表1 初始條件： ......................................... 29 表2 結構性網格與非結構性網格之比較： ................... 33 表4 主翼乾淨外型與角狀積冰比較關係： .................... 37 表6 主翼乾淨外型與平滑積冰比較關係： .................... 43 表7 乾淨水平尾翼和犄角冰第一組升阻係數比較： ............ 52 表8 乾淨水平尾翼和犄角冰第二組升阻係數比較： ............ 53 表9 乾淨水平尾翼和犄角冰第三組升阻係數比較： ............ 53 表10 乾淨水平尾翼和平滑冰升阻係數比較： ................. 54 表11 水平尾翼之力矩係數比較： ........................... 54 表12 水平尾翼力對主翼之力矩係數比較： ................... 55 表13 水平尾翼力矩係數誤差比較 ........................... 55 表14 不同模擬條件下ATR-72 之升阻比與攻角關係： .......... 60 表15 DATCOM 計算與無積冰模擬機升阻比與攻角關係： ........ 60 表16 DATCOM 計算乾淨機翼與積冰外形比較： ................ 61 表17 各種積冰狀態模組之升阻比與攻角關係： ............... 61 表18 DATCOM 積冰計算與各積冰模組數值比較： .............. 61[[note]]學號: 695430057, 學年度: 9