為了研究環(huán)境參數(shù)及外形布局對平流層飛艇氣動阻力的影響,在驗證CFD數(shù)值模擬方法的基礎(chǔ)上,從氣動阻力包括壓差阻力與摩擦阻力的角度探討了風(fēng)速、動力粘度系數(shù)、空氣密度、Re數(shù)、長細(xì)比及尾翼對飛艇氣動阻力的影響規(guī)律及機(jī)理。結(jié)果表明:氣動阻力系數(shù)隨風(fēng)速與空氣密度的增加而減小,隨動力粘度系數(shù)的增加而增加;氣動阻力系數(shù)隨Re數(shù)減小的趨勢,取決于摩擦阻力系數(shù)隨Re數(shù)的減小趨勢;隨長細(xì)比的增加,摩擦阻力系數(shù)呈現(xiàn)增加趨勢,但氣動阻力系數(shù)呈現(xiàn)先減小后增加的趨勢;尾翼對氣動阻力系數(shù)的影響主要體現(xiàn)在壓差阻力系數(shù)的改變。 

【文章來源】：工程力學(xué). 2019,36(01)北大核心EICSCD

【文章頁數(shù)】：9 頁

【部分圖文】：

風(fēng)洞模型布置圖

250工程力學(xué)驗完全一致的幾何模型；2)將模型導(dǎo)入Gambit中，定義計算域：設(shè)定模型距離地面的高度為H，模型迎風(fēng)前緣距風(fēng)速入口5H，出口距模型后緣15H，模型側(cè)面及頂部距計算域壁面5H；3)劃分網(wǎng)格：計算域采用結(jié)構(gòu)網(wǎng)格與非結(jié)構(gòu)網(wǎng)格相結(jié)合，模型周圍的核心區(qū)域采用非結(jié)構(gòu)網(wǎng)格。單元最小尺寸為0.01m。設(shè)定邊界層網(wǎng)格劃分為5層，起始邊界層網(wǎng)格厚度為單元最小尺寸的一半，即0.005m。計算域尺寸及網(wǎng)格布置如圖4所示。5H5H15HH圖4模型網(wǎng)格劃分Fig.4Meshingofthemodel邊界條件設(shè)置：采用速度入口，均勻入口流速為10m/s。出口條件為自由流出口。模型表面設(shè)置為No-Slip非滑移壁面約束，計算域四周壁面設(shè)置為Slip滑移約束。流體介質(zhì)參數(shù)與試驗環(huán)境參數(shù)一致：密度ρ0=1.225kg/m3，動力粘性系數(shù)μ0=1.7894×105kg/(m·s)。湍流模型：對于不可壓縮流體的湍流運(yùn)動，雷諾平均納維-斯托克(RANS)方程是目前工程中常用的方法，該方法的核心是將湍流運(yùn)動分為平均運(yùn)動與脈動運(yùn)動，通過湍流模式化來考慮脈動運(yùn)動對平均運(yùn)動的影響。湍流模型分為代數(shù)模型(零方程模型)、一方程模型(Spalart-Allmaras，S-A)與二方程模型(k-ε、k-ω)。其中，一方程模型是代數(shù)模型與二方程模型的折合，可考慮較復(fù)雜的湍流流動，計算效率比較高，數(shù)值穩(wěn)定性較好，已經(jīng)航空航天器的氣動研究中得到了廣泛的應(yīng)用[18－23]。計算方法：對于RANS方程，F(xiàn)luent中包含密度求解器與壓力求解器。壓力求解器主要用于低速不可壓縮流體，而密度求解器主要用于高速可壓縮流體。平流層飛?

器主要用于低速不可壓縮流體，而密度求解器主要用于高速可壓縮流體。平流層飛艇是典型的低速航空器[23―24]，其飛行速度一般不超過100km/h，采用壓力求解器對控制微分方程進(jìn)行離散，采用有限體積法將微分方程分解成一系列關(guān)于多個變量的非線性耦合代數(shù)方程組。在運(yùn)用有限體積法時，采用二階迎風(fēng)格式離散對流項，采用中心差分格式離散擴(kuò)散項，壓力速度耦合聯(lián)合采用SIMPLEC算法與PISO算法。1.3數(shù)值方法的驗證取水平姿態(tài)的試驗結(jié)果與數(shù)值結(jié)果進(jìn)行對比分析�；赥ecplot繪制艇體風(fēng)壓分布圖，圖5與圖6分別給出對應(yīng)的俯視圖與仰視圖。觀察發(fā)現(xiàn)，數(shù)值與試驗所得艇體的風(fēng)壓分布規(guī)律一致。(a)試驗(b)數(shù)值圖5試驗與數(shù)值風(fēng)壓分布對比圖(俯視圖)Fig.5Comparisonofwindpressuredistributionbetweenexperimentandsimulation(topview)(a)試驗(b)數(shù)值圖6試驗與數(shù)值風(fēng)壓分布對比圖(仰視圖)Fig.6Comparisonofwindpressuredistributionbetweenexperimentandsimulation(bottomview)為了進(jìn)一步驗證數(shù)值方法的正確性，提取子午線φ=30°與φ=150°上測點對應(yīng)的風(fēng)壓系數(shù)進(jìn)行對比。規(guī)定艇體下垂尾對應(yīng)的子午線為φ=180°，如圖7所示。圖7子午線位置示意圖Fig.7Locationofmeridians

【參考文獻(xiàn)】：
期刊論文
[1]平流層飛艇氣動阻力的數(shù)值模擬及公式擬合[J]. 李天娥,孫曉穎,張中昱,林國昌.  工程力學(xué). 2017(08)
[2]平流層飛艇的多目標(biāo)優(yōu)化與決策[J]. 孫曉穎,李天娥,陸正爭,武岳,王長國.  工程力學(xué). 2015(06)
[3]平流層飛艇充氣柔性膜結(jié)構(gòu)的區(qū)間不確定優(yōu)化[J]. 仇翯辰,邱志平.  工程力學(xué). 2015(04)
[4]飛艇大攻角繞流氣動特性模擬及湍流模型與參數(shù)影響研究[J]. 吳小翠,王一偉,黃晨光,杜特專,于嫻嫻,廖麗涓.  工程力學(xué). 2014(08)
[5]平流層飛艇大迎角氣動特性分析[J]. 黃龍?zhí)?張超,姜琬.  航空工程進(jìn)展. 2013(04)
[6]基于響應(yīng)面法的某雙體飛艇囊體氣動外形優(yōu)化[J]. 曹鵬鈞,姜琬,張華.  航空計算技術(shù). 2012(01)
[7]飛艇空氣動力學(xué)及其相關(guān)問題[J]. 任一鵬,田中偉,吳子牛.  航空學(xué)報. 2010(03)
[8]駐點引射飛艇減阻數(shù)值模擬[J]. 謝飛,葉正寅.  工程力學(xué). 2010(02)
[9]平流層雙軸橢球體飛艇繞流場的數(shù)值分析[J]. 張丹,郭雪巖.  力學(xué)季刊. 2008(04)
[10]平流層飛艇放飛、回收過程初步分析[J]. 趙攀峰,王永林,劉傳超.  航空科學(xué)技術(shù). 2007(04)

博士論文
[1]張拉膜結(jié)構(gòu)氣彈失穩(wěn)機(jī)理研究[D]. 陳昭慶.哈爾濱工業(yè)大學(xué) 2015



本文編號：3365866