張傳俠，呂水燕，謝波濤，王寶林

(中國(guó)兵器工業(yè)試驗(yàn)測(cè)試研究院， 陜西 華陰 714200)

地面效應(yīng)[1-2]是一種使飛行器誘導(dǎo)阻力減小，同時(shí)能獲得比空中飛行更高升阻比的流體力學(xué)效應(yīng)，當(dāng)飛行器貼地飛行時(shí)，飛行器體上下表面壓力差增加，升力急劇增加。低速地面效應(yīng)主要是由于地面影響導(dǎo)致下表面壓力提高而升力增加，展向繞翼梢的下洗氣流受地面阻擋減小了機(jī)翼下洗；而超音速地面效應(yīng)主要由于復(fù)雜的激波反射、激波-激波干擾、激波-邊界層干擾等影響而造成氣動(dòng)力及力矩劇烈變化，特別是復(fù)雜模型的貼地飛行器地面效應(yīng)問(wèn)題更加突出。針對(duì)簡(jiǎn)單模型的地面效應(yīng)問(wèn)題，國(guó)外主要通過(guò)風(fēng)洞實(shí)驗(yàn)和地面火箭橇試驗(yàn)的手段研究，而對(duì)于超音速?gòu)?fù)雜模型的地面效應(yīng)研究主要依靠數(shù)值模擬計(jì)算方法。

火箭橇試驗(yàn)[3-4]是介于風(fēng)洞實(shí)驗(yàn)與飛行試驗(yàn)之間的一種以火箭發(fā)動(dòng)機(jī)為動(dòng)力，沿專用軌道滑行來(lái)模擬被試品或部件需要的速度、過(guò)載、力學(xué)環(huán)境等參數(shù)的貼地飛行試驗(yàn)。較飛行器高空飛行相比，有翼火箭橇貼地高速滑行且結(jié)構(gòu)外形復(fù)雜，受到地面效應(yīng)的影響更加強(qiáng)烈，側(cè)翼氣動(dòng)特性將會(huì)劇烈變化。于是，有翼火箭橇的側(cè)翼不僅要充分借助機(jī)翼升力原理，而且要合理利用地面效應(yīng)優(yōu)勢(shì)，使之能夠持續(xù)提供氣動(dòng)壓力或升力，提高主體橇在軌高速運(yùn)行的穩(wěn)定性。因此，針對(duì)有翼火箭橇的側(cè)翼有必要開展強(qiáng)地效環(huán)境下的氣動(dòng)特性研究。

基于國(guó)外超音速火箭橇復(fù)雜結(jié)構(gòu)模型地面效應(yīng)CFD研究方法[5]，本文針對(duì)強(qiáng)地效環(huán)境下超音速有翼火箭橇側(cè)翼氣動(dòng)特性變化情況，采用滑移壁面模擬地面相對(duì)運(yùn)動(dòng)邊界條件，通過(guò)數(shù)值模擬的方法對(duì)有無(wú)地面效應(yīng)、側(cè)翼攻角變化、側(cè)翼連接位置變化等工況條件側(cè)翼氣動(dòng)特性開展研究，尋找多種工況條件下側(cè)翼氣動(dòng)特性影響規(guī)律，為強(qiáng)地效環(huán)境下有翼火箭橇側(cè)翼優(yōu)化設(shè)計(jì)提供理論依據(jù)。

1 計(jì)算模型

計(jì)算模型[6](圖1)以有翼火箭橇結(jié)構(gòu)設(shè)計(jì)標(biāo)準(zhǔn)模型為研究對(duì)象，該標(biāo)準(zhǔn)模型由主體橇、側(cè)翼、發(fā)動(dòng)機(jī)及滑靴等結(jié)構(gòu)部件組成；側(cè)翼為超音速常用尖前緣三角翼型沿展向拉伸生成；翼根與主體橇連接，起到支撐主體橇作用，以提高在軌高速運(yùn)行穩(wěn)定性；翼尖與滑靴連接。

具體研究如下內(nèi)容：

1) 側(cè)翼繞坐標(biāo)原點(diǎn)(以x軸為轉(zhuǎn)軸)旋轉(zhuǎn)-6°、-3°、3°、6°，研究側(cè)翼攻角姿態(tài)變化對(duì)氣動(dòng)特性的影響情況；

2) 側(cè)翼繞翼尖前緣頂點(diǎn)(以y軸為轉(zhuǎn)軸)向上偏轉(zhuǎn)5°、向下偏轉(zhuǎn)5°，研究側(cè)翼與主體橇連接位置上下偏移對(duì)氣動(dòng)特性的影響情況；

3) 側(cè)翼翼根連接位置前移300 mm、后移300 mm，研究側(cè)翼與主體橇連接位置前后變化對(duì)氣動(dòng)特性的影響情況。

采用ICEM軟件建立數(shù)值計(jì)算遠(yuǎn)場(chǎng)邊界，并對(duì)計(jì)算區(qū)域進(jìn)行非結(jié)構(gòu)四面體網(wǎng)格生成，計(jì)算區(qū)域選擇直徑約為40倍側(cè)翼展長(zhǎng)的半圓柱體，且圓柱高度約為70倍側(cè)翼平均氣動(dòng)弦長(zhǎng)，網(wǎng)格拓?fù)涔?jié)點(diǎn)數(shù)約180萬(wàn)，網(wǎng)格單元規(guī)模約1 000萬(wàn)。為了更加詳盡的捕捉計(jì)算模型側(cè)翼附近、側(cè)翼與地面軌道系統(tǒng)之間及側(cè)翼與主體橇拐角區(qū)的流動(dòng)特征，對(duì)于側(cè)翼及拐角區(qū)附近采用局部網(wǎng)格加密處理，特別是側(cè)翼尖銳前緣與后緣線進(jìn)行線節(jié)點(diǎn)加密處理。計(jì)算模型網(wǎng)格如圖2所示。

2 計(jì)算方法[7-10]

本文采用非結(jié)構(gòu)網(wǎng)格求解可壓縮雷諾平均Navier-Stokes方程，對(duì)計(jì)算模型附近流動(dòng)進(jìn)行數(shù)值模擬計(jì)算，湍流模型選擇SST兩方程，湍流動(dòng)能采用高階迎風(fēng)格式離散。

進(jìn)出口邊界條件：流場(chǎng)入口設(shè)定速度進(jìn)口邊界條件給定速度、相對(duì)壓力與溫度，流場(chǎng)出口設(shè)定壓力出口邊界條件給定壓力。

壁面邊界條件：計(jì)算模型表面與壓力遠(yuǎn)場(chǎng)為壁面邊界，模型表面設(shè)定為固壁絕熱無(wú)滑移邊界(No Slip Wall)，遠(yuǎn)場(chǎng)設(shè)定自由滑移邊界(Free Slip Wall)。

地面軌道系統(tǒng)設(shè)定固定速度滑移邊界以模擬火箭橇計(jì)算模型與地面軌道系統(tǒng)之間的相對(duì)運(yùn)動(dòng)。

流場(chǎng)入口速度680 m/s，溫度設(shè)定為288.15 K，大氣壓強(qiáng)設(shè)定為101 325 Pa。

自由狀態(tài)(模擬標(biāo)準(zhǔn)計(jì)算模型無(wú)地面效應(yīng)工況)，地面試驗(yàn)狀態(tài)(模擬標(biāo)準(zhǔn)計(jì)算模型有地面效應(yīng)工況)，側(cè)翼攻角及連接位置變化均為地面試驗(yàn)狀態(tài)。

3 計(jì)算結(jié)果對(duì)比分析

為驗(yàn)證計(jì)算方法的可靠性[11]，本文選取超音速菱形翼型沿展向拉伸1 m的三維機(jī)翼模型為驗(yàn)證對(duì)象，在馬赫數(shù)為2.0、雷諾數(shù)250萬(wàn)(基于平均氣動(dòng)弦長(zhǎng)1 m)狀態(tài)下進(jìn)行數(shù)值模擬，對(duì)比了機(jī)翼攻角為0°時(shí)，機(jī)翼任意一截面菱形翼型的壓力系數(shù)分布。圖3給出了機(jī)翼表面采用本文計(jì)算方法所計(jì)算壓力系數(shù)分布與理論計(jì)算結(jié)果，從菱形翼型壓力系數(shù)分布的對(duì)比可以看出，數(shù)值計(jì)算結(jié)果與理論計(jì)算結(jié)果吻合較好，驗(yàn)證了本文計(jì)算方法的可靠性。

1) 有無(wú)地面效應(yīng)側(cè)翼氣動(dòng)特性研究

為了研究有無(wú)地面效應(yīng)工況條件下，有翼火箭橇側(cè)翼氣動(dòng)力的變化情況，表1給出了有翼火箭橇側(cè)翼自由狀態(tài)(無(wú)地面效應(yīng))與地面試驗(yàn)狀態(tài)(有地面效應(yīng))工況條件下氣動(dòng)升力與氣動(dòng)阻力的數(shù)據(jù)。從表中數(shù)據(jù)可以看出有無(wú)地面效應(yīng)工況條件下，側(cè)翼氣動(dòng)阻力基本沒有變化；氣動(dòng)升力變化較大，且表現(xiàn)為氣動(dòng)下壓力，較無(wú)地面效應(yīng)相比，氣動(dòng)升力增加41%。

表1 自由狀態(tài)與地面試驗(yàn)狀態(tài)側(cè)翼氣動(dòng)力

為了更加詳細(xì)的描述有無(wú)地面效應(yīng)工況條件下有翼火箭橇側(cè)翼展向截面壓力系數(shù)分布及側(cè)翼上下翼面的流動(dòng)變化情況，圖4給出了側(cè)翼表面展向50%處翼型無(wú)地面效應(yīng)狀態(tài)與有地面效應(yīng)狀態(tài)壓力系數(shù)分布曲線，圖5給出了無(wú)地面效應(yīng)狀態(tài)與有地面效應(yīng)狀態(tài)側(cè)翼上下翼面壓力云圖。從圖4、圖5可以看出有無(wú)地面效應(yīng)條件下側(cè)翼上翼面氣動(dòng)壓力基本沒有變化；下翼面變化較大，呈現(xiàn)出受強(qiáng)地面效應(yīng)影響后氣動(dòng)壓力明顯增大(負(fù)壓變?yōu)檎龎?，增幅接近80%，且翼根下翼面附近受地面激波反射及主體橇氣動(dòng)繞流影響，拐角區(qū)流動(dòng)更加復(fù)雜，氣動(dòng)壓力變化更為劇烈。

2) 側(cè)翼攻角變化氣動(dòng)特性研究

為了研究強(qiáng)地效環(huán)境下有翼火箭橇側(cè)翼攻角變化對(duì)其氣動(dòng)特性的影響情況，圖6給出了地面試驗(yàn)狀態(tài)下不同攻角姿態(tài)側(cè)翼的氣動(dòng)升力與氣動(dòng)阻力系數(shù)曲線。從圖6數(shù)據(jù)可以看出隨著側(cè)翼攻角變大，氣動(dòng)升力不斷增加(從-6°下壓力變?yōu)?°上升力)，氣動(dòng)阻力不斷變小，并且氣動(dòng)升力變化比氣動(dòng)阻力變化更為劇烈；側(cè)翼升阻比絕對(duì)值表現(xiàn)為先減小后增大。

3) 側(cè)翼連接位置變化氣動(dòng)特性研究

為了研究強(qiáng)地效環(huán)境下有翼火箭橇側(cè)翼與主體橇連接位置變化對(duì)側(cè)翼氣動(dòng)特性的影響，圖7給出了地面試驗(yàn)狀態(tài)側(cè)翼上下位置變化氣動(dòng)升力與氣動(dòng)阻力系數(shù)變化曲線，圖8給出了地面試驗(yàn)狀態(tài)側(cè)翼前后位置變化氣動(dòng)升力與氣動(dòng)阻力系數(shù)變化曲線。從圖7、圖8可以看出，側(cè)翼上下位置的變化對(duì)氣動(dòng)升力有一定的影響，表現(xiàn)為隨著側(cè)翼翼根位置的上偏氣動(dòng)升力(下壓力)逐漸減小，且位置向上變化比向下變化氣動(dòng)升力變化幅度大；側(cè)翼前后位置的變化對(duì)其氣動(dòng)力基本沒有影響，即側(cè)翼氣動(dòng)特性對(duì)前后位置的變化不敏感。

4 結(jié)論

1) 基于強(qiáng)地面效應(yīng)環(huán)境下側(cè)翼翼型氣動(dòng)特性變化劇烈；2) 通過(guò)設(shè)計(jì)合理的超音速有翼火箭橇側(cè)翼翼型及結(jié)構(gòu)布置形式，可充分利用地面效應(yīng)優(yōu)勢(shì)，有效避免強(qiáng)地效激波反射作用，為火箭橇提供足夠大小的氣動(dòng)升力或壓力，提高火箭橇在軌高速運(yùn)行過(guò)程中的穩(wěn)定性。