基于Ansys-CFX的飛機(jī)機(jī)翼瞬態(tài)流場(chǎng)分析

孫巍 張博 郭鵬星

(渤海大學(xué)工學(xué)院,遼寧錦州 121000)

基于Ansys-CFX的飛機(jī)機(jī)翼瞬態(tài)流場(chǎng)分析

孫巍 張博 郭鵬星

(渤海大學(xué)工學(xué)院,遼寧錦州 121000)

機(jī)翼是飛機(jī)飛行必不可少的組成部分,了解其在飛行過程中周圍流場(chǎng)的狀態(tài)對(duì)設(shè)計(jì)者在設(shè)計(jì)機(jī)翼結(jié)構(gòu)時(shí)起到非常關(guān)鍵的作用。利用ANSYS-CFX軟件建立飛機(jī)機(jī)翼有限元模型,通過建模和設(shè)置邊界條件分別對(duì)機(jī)翼低音速,超音速和臨界音速三個(gè)約束條件下機(jī)翼的狀態(tài)進(jìn)行分析,得到相應(yīng)的數(shù)據(jù)及云圖。分析結(jié)果得出在飛機(jī)飛行過程中,機(jī)翼周圍的流速與壓力大小成反比,并且機(jī)翼最前端所受的壓力最大,同時(shí)由壓力云圖可以得出,飛機(jī)在接近音速條件下其周圍所受的壓力最大,在飛機(jī)突破音障之后,機(jī)翼周圍的壓力反而減小。同時(shí)由仿真結(jié)果得到一些重要數(shù)據(jù),為高空飛行時(shí)機(jī)翼的設(shè)計(jì)和改進(jìn)提供了數(shù)據(jù)依據(jù)。

機(jī)翼設(shè)計(jì) 繞流 音障

這些年來，世界各國越來越重視高空長(zhǎng)時(shí)間飛行的飛機(jī)設(shè)計(jì)與優(yōu)化。根據(jù)長(zhǎng)時(shí)間飛行飛機(jī)的時(shí)間性能的要求，對(duì)飛機(jī)各部件的要求也越來越嚴(yán)格。機(jī)翼是飛機(jī)的重要組成部件之一。它的最主要的作用是飛行過程中提供起飛所需要的升力，同時(shí)使飛機(jī)具有操作性和橫側(cè)的穩(wěn)定性[1]。因此，機(jī)翼結(jié)構(gòu)力學(xué)穩(wěn)定性對(duì)保證整個(gè)飛機(jī)在飛行過程中的安全起著非常重要的作用[2][3]。目前很多文獻(xiàn)研究非對(duì)稱翼，對(duì)非對(duì)稱翼簡(jiǎn)單來說上表面比較凸，因此流過機(jī)翼上表面的流速較快；而機(jī)翼下表面比較平因此流過機(jī)翼下表面的氣流較慢，與上表面正好相反。根據(jù)流體力學(xué)的基本原理，流動(dòng)慢的大氣壓強(qiáng)較大，而流動(dòng)快的大氣壓強(qiáng)較小，這樣機(jī)翼下表面的壓強(qiáng)就比上表面的壓強(qiáng)高，大氣施加與機(jī)翼下表面的壓力(方向向上)比施加于機(jī)翼上表面的壓力(方向向下)大，二者的壓力差便形成了飛機(jī)的升力[4][5]。本文利用流體力學(xué)仿真軟件ANSYS-CFX對(duì)一種經(jīng)典的對(duì)稱型固定翼模型NACA0012的機(jī)翼外部氣體流動(dòng)狀態(tài)進(jìn)行模擬仿真，通過CFX模擬計(jì)算并觀察在三種不同范圍的速度下，機(jī)翼周圍的壓力大小，氣流分布及氣體軌跡狀態(tài)。

圖1 CFX 軟件結(jié)構(gòu)圖

圖2 亞音速時(shí)的速度云圖

圖3 亞音速時(shí)的壓力云圖

圖4 跨音速時(shí)的速度云圖

圖5 跨音速時(shí)的壓力云圖

1 ANSYS CFX介紹

ANSYS-CFX流體力學(xué)軟件是全球第一個(gè)通過ISO9001質(zhì)量認(rèn)證的于流體計(jì)算的大型商業(yè)CFD軟件，可以用于亞音速、跨音速和超音速流體的數(shù)值模擬。其功能的實(shí)現(xiàn)主要由5部分組成[6]:

(1)幾何建模:創(chuàng)建或?qū)霂缀斡颍豢刂品匠淘谄渲袑⒈磺蠼猓瑫r(shí)獲得在這個(gè)區(qū)域中的結(jié)果。

(2)網(wǎng)格劃分:創(chuàng)建有限體積或是單元的過程；網(wǎng)格可以在另一個(gè)模塊中創(chuàng)建或是從其他的網(wǎng)格劃分器中導(dǎo)入。

(3)前處理:求解問題描述，定義物理模型，材料屬性和邊界條件。

(4)求解求解控制方程:設(shè)置求解器選項(xiàng)；將迭代求解控制方程過程作為一個(gè)批處理的過程，求解收斂。(5)后處理分析和顯示求解結(jié)果。其結(jié)構(gòu)及過程如圖1所示。

該軟件采用有限體積法進(jìn)行離散，采用K-標(biāo)準(zhǔn)，所用控制方程如下。

其中: μeff為有效粘度；˙為修正壓力； Cε1， Cε2， σε，kσ都是常數(shù)[7][8]。

2 模型的設(shè)定及條件的設(shè)立

前人的研究表明對(duì)稱翼型在所有翼型中的阻力是最小的。為了研究對(duì)稱翼飛機(jī)在飛行過程中，其翼型阻力分布情況及其影響阻力大小的關(guān)系，因此采用典型的對(duì)稱翼型NACA0012來進(jìn)行研究。由于影響飛行中機(jī)翼的因素有很多，比如，氣體，溫度，氣流攻角，氣流速度等等。為了方便研究，在此設(shè)定飛機(jī)水平飛行即氣體攻角為0度，氣體溫度采用靜態(tài)溫度T=300K，為了進(jìn)行對(duì)比，文中選擇三個(gè)不同的氣流速度，分別是亞音速，跨音速，超音速三個(gè)范圍中的某一速度值。由于在模擬超音速時(shí)，必須保證域內(nèi)物質(zhì)設(shè)置為理想氣體，熱量傳遞模型為全熱模型，為了保證各個(gè)狀態(tài)下所處環(huán)境能相同，因此

圖6 超音速時(shí)的速度云圖

圖7 超音速時(shí)的壓力云圖

在三個(gè)不同速度的模擬中，所處環(huán)境都選擇理想氣體。

在流體域設(shè)定中，選擇理想氣體，參考?jí)毫x擇1atm，由于全熱模型適合高速流體及可壓縮流體的熱量傳輸計(jì)算，因此熱量傳輸選擇全熱傳輸，同時(shí)為了使模型適合高精度邊界層的模擬，本文中湍流模型選擇剪切壓力傳輸(SST) k-ω模型。

在邊界設(shè)定中，選擇對(duì)稱邊界及無滑移壁面邊界。在入口邊界設(shè)置中，設(shè)置湍流密度為0．01，湍流長(zhǎng)度為0．02m，靜態(tài)壓強(qiáng)為0Pa，靜態(tài)溫度為300K，速度分別為100m/s，340m/s，600m/s。

本模擬利用殘差值為0．00001的RMS，進(jìn)行100次迭代實(shí)驗(yàn)得到在不同速度下的三組仿真圖。

3 仿真結(jié)果及分析

3.1 仿真云圖(圖2-圖7)

3.2 結(jié)論及分析

(1)從3組速度云圖與壓力云圖可以看出，對(duì)于對(duì)稱翼來說，機(jī)翼上下受力情況完全相同，在機(jī)翼周圍，流速大的地方壓力小，流速小的地方壓力大，符合伯努利原理。仿真結(jié)果正確，可靠。

(2)從文中的3組圖中都可以直觀的看出，在機(jī)翼前緣處，其流速最小，并且低于飛機(jī)本來的飛行速度，同時(shí)其所受的壓力最大，因此基于這一結(jié)論，在飛機(jī)的制造中，要特別重視機(jī)翼前緣材料的抗壓能力，才能保證機(jī)翼的安全性能，并且可以得出，在飛機(jī)速度為100m/s時(shí)，其前緣的壓強(qiáng)大小為4733Pa，在飛機(jī)速度為340m/s時(shí)，其前緣的壓強(qiáng)大小為 1.023×106P a 。在飛機(jī)速度為600m/s時(shí)，其前緣的壓強(qiáng)大小為 1.987×105P a ，這些數(shù)據(jù)為機(jī)翼設(shè)計(jì)，制作提供重要的參考依據(jù)。

(3)由3組壓力云圖可以看出，當(dāng)飛機(jī)速度處于340m/s時(shí)，機(jī)翼周圍的壓力處于最大狀態(tài)，當(dāng)飛行速度超過音速時(shí)(340m/s)隨著速度的增加其機(jī)翼周圍的壓力反而減小，說明此時(shí)已經(jīng)產(chǎn)生音障現(xiàn)象，即物體的速度快要接近音速時(shí)，周邊的空氣受到聲波疊合而呈現(xiàn)非常高壓的狀態(tài)，因此一旦物體穿越音障后，周圍壓力將會(huì)陡降。音障現(xiàn)象也是曾經(jīng)阻礙飛機(jī)從亞音速到超音速突破時(shí)的最主要障礙，在飛行器設(shè)計(jì)中，不應(yīng)忽視其作用。

(4)由壓力云圖可以看出，機(jī)翼上下表面受力情況完全相同，因此如果只是以0攻角飛行的話，無法形成向上的升力，飛機(jī)是無法起飛的，因此要想起飛則需要一定的氣流攻角，讓機(jī)翼下表面的氣流速度小于機(jī)翼上表面的氣流速度，才能產(chǎn)生一定的壓力差，使之能夠起飛。

[1] 葛訊.鴨式旋翼/機(jī)翼飛行器研究[D].南京航空航天大學(xué),2012.

[2] 秦可偉,馬貴春,席園,姚光生.基于ANSYS的機(jī)翼動(dòng)力學(xué)分析[J].航空計(jì)算技術(shù),2014.2:106-109.

[3] 劉德民,劉小兵,李娟.鈦合金機(jī)翼的振動(dòng)與模態(tài)分析[J].機(jī)械設(shè)計(jì)與制造,2008.3:76-77.

[4] 杜加友.水下滑翔機(jī)本體及調(diào)節(jié)機(jī)構(gòu)研究[D].浙江大學(xué),2006.

[5] 劉鵬.翼型失速的電磁控制[D].南京理工大學(xué),2007.

[6] 高飛,李昕編著.ANSYS CFX 14.0超級(jí)學(xué)習(xí)手冊(cè)[M].人民郵電出版社,2013.

[7] 陳景秋,胡韓飛,張永祥.Star-CD對(duì)汽車外流場(chǎng)的三維數(shù)值模擬[J].重慶大學(xué)學(xué)報(bào):自然科學(xué)版,2005,28(4):99-101.

[8] 張英朝,邵書鑫.基于Star-CCM+的跑車減租設(shè)計(jì)[C]//CDAJ-China中國用戶論文集,2011.

孫巍(1992—),女,漢族,遼寧遼陽人,碩士研究生,單位:渤海大學(xué)工學(xué)院,研究方向:模式識(shí)別與智能系統(tǒng)；郭鵬星(1992—),男,漢族,河南安陽人,本科,單位:渤海大學(xué)工學(xué)院。

張博(1977—),男,漢族,副教授,單位:渤海大學(xué)工學(xué)院,研究方向:計(jì)算機(jī)網(wǎng)絡(luò),兵器科學(xué)與技術(shù)。jzhzhb@sian.com。