連續(xù)旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)隔離段流場(chǎng)數(shù)值模擬研究.docx

1、彈道學(xué)報(bào)JournalofBallisticsDOLIO.12115j,is$n.l004-499X（20H）l）2-007連續(xù)旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)隔離段流場(chǎng)數(shù)值模擬研究郭凱欣廳翁春生廳武郁文（南京理工大學(xué)瞬態(tài)物理國(guó)家重點(diǎn)實(shí)驗(yàn)室祥江蘇南京210094）摘要:為了研究連續(xù)旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)（CRDE）隔離段的反壓特性府建立了等直隔離段以及帶擴(kuò)張段的隔離段模型由用FLURNT對(duì)這2種構(gòu)型的隔離段進(jìn)行三維數(shù)值模擬由研究了旋轉(zhuǎn)爆轟波在隔離段內(nèi)向上游傳播的特性由分析了不同結(jié)構(gòu)對(duì)流場(chǎng)結(jié)構(gòu)'激波串前傳位置的影響小對(duì)同一結(jié)構(gòu)府分析了反壓大小、旋轉(zhuǎn)速度對(duì)其的影響用研究結(jié)果表明府反壓在隔離段內(nèi)與附面層相互作用府誘導(dǎo)

2、出復(fù)雜激波和膨脹結(jié)構(gòu)府形成了激波串用反壓擾動(dòng)在隔離段中前傳的距離與隔離段的構(gòu)型、反壓旋轉(zhuǎn)速度和大小有關(guān):等直隔離段的抗壓能力相對(duì)更強(qiáng)R反壓旋轉(zhuǎn)速度或壓力增加標(biāo)都會(huì)使得擾動(dòng)區(qū)域增加M數(shù)值模擬結(jié)果對(duì)CRDE隔離段反壓特性研究以及結(jié)構(gòu)設(shè)計(jì)具有一定參考價(jià)值JH關(guān)鍵詞:爆轟發(fā)動(dòng)機(jī)0旋轉(zhuǎn)爆轟反壓0隔離段0流場(chǎng)結(jié)構(gòu)中圖分類號(hào):V231.3文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：10。廠499X（2019）0廠0035-。7NumericalInvestigationofFlowFieldinContinuousRotatingDetonationEngineIsolatorGUOKa:iinWENGChunshenWUYu

3、wen（NationaKeyLaboriloryofTramien:PhyiicjiilNaniingUnirefjiiyofScience&TechnolonNiniin?210094»IChini）Abstract：Toinvestisatethecharacteristicsofthebackpressureinacontinuousrotatingdetonationengine（CRDE）ijolaloriii'themodelsofslraightisolalorandisolatorwithanexpansionseclionwereestablisbe

4、d.Thethreedimensionalnumericalemulationontheisolatoroitheiwoconiigurationswascarriedou:byFLURNT.Thecharacteristicsoftherotatingdetonationwavepropagatingupdreamintheisolatorwereanalyzed.TheinfluencesofdifferentitructuresonlhefhwfieldHruclureandtheforwardpositionoftheshocktrainwereanalyzed.Forthesames

5、tructureittheiofIuencesofthebackpressureandtherotationa1speedonitwereanalyzed.TheresultsshowthatlhebackpressureinleractingwiththeboundarylayerinlheitolalorcaninducecomplexshockandexpansionJiruciureslitformingashocklia:n.Thedistanceforwardedbythebackprenureintheisolatorisrela:edIotheconfirmationofthe

6、isolatorHiandlhespeedandmagnitudeoflhebackpressurerotation.TherelativelyenhancementollhepressureresistanceofthestraightisoIatoritbeincreasingofthebackpressurespeedorpressurecanincreasethedisturbancearea.Thenumericalsimulationresultshavecertainreiereneevalueforthebackpressurecharaclerislicsandstructu

7、raldesignoftheisolatedrotatingdetonationengine.Keywords：detonationengineRotatingdelonalionbackpieisureJPisolatorjPflowfieldstructure收稿日期：20)8-10-25基金項(xiàng)目：國(guó)家自然科學(xué)基金項(xiàng)目(702143)JI國(guó)防預(yù)研基金項(xiàng)目(6142604D201162604002)作者簡(jiǎn)介：郭凱欣(1994-女撲頁(yè)士研究生R研究方向?yàn)檫B續(xù)旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)IfE_mi'：461708539l#通信作者：翁春生(1964-)if男it教授it博士生導(dǎo)師il研兗方向?yàn)橥七M(jìn)技

8、術(shù)NE_miil;wenicj,連戟轉(zhuǎn)像轟發(fā)動(dòng)機(jī)(continuousrotatingdetonationenginetfCRDE)是一種新型發(fā)動(dòng)機(jī)仲現(xiàn)有的航空航天動(dòng)力系統(tǒng)主要基于常規(guī)的等壓燃燒方式且此類發(fā)動(dòng)機(jī)的發(fā)展己經(jīng)相對(duì)成熟R性能難以再取得重大突破沖與傳統(tǒng)的發(fā)動(dòng)機(jī)相tWCRDE基于爆轟燃燒£爆轟波近似于等容燃燒貴理論上具有比等壓燃燒更高的熱循環(huán)效率和熱釋放速率村專播速度可達(dá)千米每秒量級(jí)R具有結(jié)構(gòu)緊湊、燃燒速率快'熱效率高、高比沖等特點(diǎn)用此外R爆轟波沿與流向垂直的圓周方向的動(dòng)態(tài)傳播也增加了燃燒的穩(wěn)定性1，_23I但是在工作過(guò)程中V燃燒室入口處的壓力在空間上非均勻分布f且每個(gè)

9、位置的壓力都隨時(shí)間變化麻其反壓環(huán)境會(huì)對(duì)上游來(lái)流產(chǎn)生較大影響R因此在燃燒室上游需設(shè)置隔離段月隔離段位于進(jìn)氣道與燃燒室之間#主要作用是承受燃燒室的反壓環(huán)境R同時(shí)防止燃燒室產(chǎn)生的擾動(dòng)前傳而影響進(jìn)氣道的流動(dòng)狀態(tài)苜溫玉芬5通過(guò)數(shù)值模擬與實(shí)驗(yàn)相結(jié)合的方法府對(duì)脈沖爆轟發(fā)動(dòng)機(jī)進(jìn)氣道的流動(dòng)特性進(jìn)行了研究府分析了周期性爆轟壓力的擾動(dòng)對(duì)進(jìn)氣道流場(chǎng)結(jié)構(gòu)的干擾作用玳以及不同爆轟波特性、不同進(jìn)氣道布局形式對(duì)進(jìn)氣道流動(dòng)特性的影響就王超等e開展了吸氣式連續(xù)旋轉(zhuǎn)爆轟試驗(yàn)竹對(duì)隔離段和燃燒室壓力進(jìn)行測(cè)量針?lè)治隽诉B續(xù)旋轉(zhuǎn)爆轟對(duì)空氣來(lái)流的作用竹并考慮了燃燒室尺寸的影響用試驗(yàn)發(fā)現(xiàn)了吸氣式連續(xù)旋轉(zhuǎn)爆轟與來(lái)流的3種相互作用：來(lái)流總壓不變府且

10、不存在高頻壓力震蕩0來(lái)流總壓不變玳但存在高頻壓力振蕩府且與連續(xù)旋轉(zhuǎn)爆轟波傳播頻率相同川來(lái)流存在與連續(xù)旋轉(zhuǎn)爆轟波傳播頻率相同的高頻壓力振蕩浦來(lái)流總壓升高礦且隨著燃燒室流通面積的減小仃連續(xù)旋轉(zhuǎn)爆轟波引起的隔離段中的高頻壓力振蕩峰值增加巾頻率上升*對(duì)來(lái)流的影響增強(qiáng)劇Su等也通過(guò)數(shù)值模擬研究了正弦脈動(dòng)反壓下超燃沖床發(fā)動(dòng)機(jī)隔離段壁面的壓力振蕩特性用上述文獻(xiàn)中的研究沒(méi)有涉及到連續(xù)旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)的隔離段的構(gòu)型對(duì)其反壓特性的影響用隔離段內(nèi)的流動(dòng)特性不僅受入口與出口條件的影響面還與其幾何構(gòu)型有著密切的關(guān)系用常見(jiàn)的隔離段構(gòu)型為等直管道叫旦超燃沖壓發(fā)動(dòng)機(jī)燃燒室里普遍采用了擴(kuò)張角和后臺(tái)階的結(jié)構(gòu)面此類結(jié)構(gòu)一般也承擔(dān)著一

11、定的抗反壓的作用用因此本文用FLUENT軟件對(duì)等直及漸擴(kuò)隔離段進(jìn)行了三維數(shù)值模擬府研究了旋轉(zhuǎn)爆轟波在隔離段內(nèi)向上游傳播的特性面分析了隔離段結(jié)構(gòu)對(duì)流場(chǎng)結(jié)構(gòu)和激波串前傳位擔(dān)的影響仲對(duì)同一結(jié)構(gòu)R分析了反壓大小和旋轉(zhuǎn)速度對(duì)其的影響玳為后續(xù)的研究奠定了基礎(chǔ)用1計(jì)算模型與數(shù)值方法1.1計(jì)算模型連續(xù)旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)的燃燒室通常為環(huán)形結(jié)構(gòu)就因此選擇環(huán)形結(jié)構(gòu)的隔離段用隔離段包括等直構(gòu)型、帶擴(kuò)張段構(gòu)型竹如圖1所示據(jù)其中等直隔離段總長(zhǎng)度為90mm#內(nèi)徑為40mmHT外徑為60mmji1帶擴(kuò)張段構(gòu)型的隔離段總長(zhǎng)度為170mm由其中擴(kuò)張段長(zhǎng)度為100mm#入口處外徑為44mm#擴(kuò)張角度為9.1°I采用網(wǎng)格前處理

12、軟件1CEM對(duì)模型進(jìn)行網(wǎng)格劃分宿吏用六面體結(jié)構(gòu)化網(wǎng)格府為了保證計(jì)算精度®同時(shí)控制網(wǎng)格數(shù)量面在隔離段主流方向、半徑方向和圓周方向分別布置了180、36、200個(gè)網(wǎng)格訴網(wǎng)格總數(shù)約為128萬(wàn)用為了更好地反映近壁面的附面層流動(dòng)面還對(duì)內(nèi)外壁網(wǎng)格進(jìn)行了局部加密沖壁面第一層網(wǎng)格高度為0.Olmmttt保證壁面處y'UORy為第一層網(wǎng)格中心到壁面的無(wú)量綱距離映圖1計(jì)算模型1.2計(jì)算方法本文利用FLUENT求解器府基于理想氣體假設(shè)if采用UDF模擬旋轉(zhuǎn)爆轟反壓府不考慮化學(xué)反應(yīng)祥考慮粘性作用砰采用密度基顯式算法求解N-S方程明對(duì)流項(xiàng)選用對(duì)激波具有較高捕捉精度的三階MUSCL格式進(jìn)行離散命物理通量

13、利用AUSM格式進(jìn)行分解#時(shí)間項(xiàng)采用Runge'Kutla法進(jìn)行處理山中為了能夠較好地捕捉旋轉(zhuǎn)爆轟反壓與邊界層的相互作用的時(shí)間步長(zhǎng)取0.01uslit以保證計(jì)算精度用1.3邊界條件與初始條件本文計(jì)算區(qū)域的邊界條件有如下3種網(wǎng) 壁面邊界條件用隔離段內(nèi)外壁面均為無(wú)滑移、無(wú)穿透、絕熱固壁邊界用 壓力入口邊界條件描入口馬赫數(shù)Ma=1.5HT總壓p.=0.4MPaiJ靜壓Pi=0.107MP河總溫為T=606KW 壓力出口邊界條件用初始時(shí)刻面對(duì)隔離段內(nèi)流場(chǎng)以入口條件初始化評(píng)為避免出口反壓在超聲速條件下壓力外推府對(duì)出口處網(wǎng)格分區(qū)初始化府并對(duì)其賦予亞聲速初始條件用本文考察了不同反壓大小、爆轟波傳播速

14、度對(duì)流場(chǎng)的影響由反壓大小由恢復(fù)區(qū)壓力和壓力振幅組成R分別對(duì)應(yīng)圖2(a)中的ppB值府相應(yīng)的算例初始參數(shù)見(jiàn)表1M表中ifv為爆轟波傳播速度N為了模擬由于爆轟波在環(huán)形燃燒室內(nèi)持續(xù)旋轉(zhuǎn)傳播產(chǎn)生的旋轉(zhuǎn)反壓e*本文用指數(shù)函數(shù)擬合了反壓振型if反壓沿圓周方向的分布如圖2(a)所示巾并通過(guò)UDF導(dǎo)入算例作為隔離段出口處的壓力邊界條件巾圖2(D中祀為隔離段環(huán)面的極角坐標(biāo)用由函數(shù)擬合得到的隔離段出口邊界壓力云圖與文獻(xiàn)8中旋轉(zhuǎn)爆轟燃燒室入口處的壓力云圖如圖2(b).圖2(c)所示#由圖可知W沿周向分布的壓力變化趨勢(shì)是一致的用表1算例及初始條件參數(shù)算例隔離段構(gòu)型v(m.'1>p./MPapi/MPa1

15、等直180(10.150.302等直180(10.300.603擴(kuò)張18000.150,304擴(kuò)張16000.150.305擴(kuò)張1OOO0.150.306擴(kuò)張18000.130.267擴(kuò)張18000.100.20圖2隔離段出口壓力分布2計(jì)算結(jié)果及其討論2.1等直隔離段內(nèi)流場(chǎng)結(jié)構(gòu)分析圖3和圖4給出了算例1和算例2外壁面、隔離段截面上的壓力與馬赫數(shù)分布云圖用由圖3可知面當(dāng)爆轟反壓較低時(shí)仃旋轉(zhuǎn)反壓上升至一定高度以后的不再繼續(xù)向上游傳播祐甚至被抑制在隔離段出口附近用由圖4可知帥增大反壓后的反壓不斷向上游傳播用拖曳出一道斜激波面的沿著軸線方向螺旋上升祥受擾動(dòng)的區(qū)域逐漸向上游擴(kuò)展用從壓力及馬赫數(shù)分布云圖上

16、可以看出價(jià)沿氣流流動(dòng)的方向的附面層逐漸加厚命靜壓逐漸升高的馬赫數(shù)逐漸降低用隨著反壓的擾動(dòng)沿附面層向上游傳播仃在附面層附近產(chǎn)生壓縮波面從而形成斜激波面在斜激波后的壓力升高很多用斜激波在中心線相交以后府被互相壓縮府又產(chǎn)生了斜激波系麻斜激波系反射向壁面府并和壁面撞擊井壁面處的壓強(qiáng)急劇升高府從而使得附面層分離用氣流分離區(qū)與核心區(qū)壓強(qiáng)的差異導(dǎo)致激波在壁面附近發(fā)生普朗特-邁耶反射討形成稀疏膨脹波討膨脹波和斜激波重復(fù)出現(xiàn)5K形成了激波串成這與文獻(xiàn)1。-12中在圓截面等直隔離段中觀察到的激波串結(jié)構(gòu)是相似的由沿著氣流的方向浦隔離段中線處的靜壓交替上升、下降質(zhì)馬赫數(shù)的變化趨勢(shì)則剛好相反明可見(jiàn)廳反壓的增加會(huì)使激波串

17、出現(xiàn)并向上游移動(dòng)用算例2中反壓不斷前傳礦直至將首道斜激波推出隔離段R整個(gè)流場(chǎng)變成亞聲速進(jìn)氣條件府因此在該邊界條件下砰隔離段滿足不了抑制旋轉(zhuǎn)反壓的要求出圖3算例1流場(chǎng)分布圖4算例2流場(chǎng)分布2.2帶擴(kuò)張段的隔離段內(nèi)流場(chǎng)結(jié)構(gòu)分析通常認(rèn)為入口馬赫數(shù)越大R隔離段的抗反壓性能越強(qiáng)府但是引入擴(kuò)張角將會(huì)逐漸改變隔離段的高度、邊界層的厚度、馬赫數(shù)等條件R實(shí)際上沿程不僅馬赫數(shù)增大了府壓力也減小了質(zhì)而壓力的減小又對(duì)抗反壓能力不利竹綜合考慮情況較為復(fù)雜用因此府考慮在原本的等直隔離段構(gòu)型前面加一段擴(kuò)張段廳觀察其流動(dòng)特性用圖5為帶擴(kuò)張段的隔離段在不加反壓時(shí)的純流動(dòng)流場(chǎng)M從圖中可以看出府在擴(kuò)張段內(nèi)面來(lái)流不斷加速就馬赫數(shù)達(dá)到

18、3以上府壓力降低目經(jīng)過(guò)等直段時(shí)面壓力和速度基本保持不變用圖6為算例3的流場(chǎng)分布圖用圖5帶擴(kuò)張段隔離段在無(wú)反壓時(shí)的純流動(dòng)流場(chǎng)圖6算例3流場(chǎng)分布從圖6中可以看出府由于擴(kuò)張段的存在祁付面層的分離更加嚴(yán)重由而且厚度迅速增加帥激波與附面層的干擾更強(qiáng)府反壓螺旋著向上游傳播面相比等直隔離段更容易沿著附面層向前傳面且反壓傳至靠近擴(kuò)張段入口處的某一位置后面不再繼續(xù)向上游發(fā)展府該位置即代表爆轟波所能影響的最上游的位置由但是本文算例中該位置并不是一個(gè)平面府暫且將其最上游的位置稱為“極限位置”用以算例3作為基準(zhǔn)算例府控制變量府分別改變爆轟波的旋轉(zhuǎn)速度和壓力峰值由與算例4、算例5和算例6、算例7做比較用2.2.1反壓旋

19、轉(zhuǎn)速度的影響圖7、圖8中反壓的旋轉(zhuǎn)速度分別為1600m/s#1000m/sf云圖形狀基本與算例3類似W由于擴(kuò)張段單側(cè)擴(kuò)張#內(nèi)壁面的附面層分離現(xiàn)象要遠(yuǎn)小于外壁面面在經(jīng)過(guò)極限位置的一道斜激波以后B氣流速度降為亞聲速寸在擴(kuò)張段內(nèi)逐漸減速#到等直段時(shí)面馬赫數(shù)己經(jīng)非常小了不利于進(jìn)氣*將算例3、算例4、算例5的云圖對(duì)比可以發(fā)現(xiàn)戶隨著反壓旋轉(zhuǎn)速度的下降R壓力云圖上這條螺旋的高壓線的“螺距”增加了由反壓能到達(dá)的“極限位置”降低用2.2.2反壓大小的影響圖9和圖10分別給出了算例6、算例7外壁面、軸截而上的壓力與馬赫數(shù)分布云圖用算例6、算例7中反壓的大小分別為0.26MPa和0.20MPaf隨著出口反壓的降低府

20、受擾動(dòng)的區(qū)域減少W首道斜激波的位置下降面即“極限位置”向后推移了甘齡圖9算例6流場(chǎng)分布圖7和圖8分別給出了算例4、算例5外壁面、軸截面上的壓力與馬赫數(shù)分布云圖做潮融蹣花圖10算例7流場(chǎng)分布做潮融蹣花圖10算例7流場(chǎng)分布2.3旋轉(zhuǎn)爆轟波壓力及傳播速度對(duì)隔離段流動(dòng)特性的影響2.3.1對(duì)等直隔離段的影響圖11和圖12為不同反壓下等直隔離段沿圓環(huán)中心軸線及內(nèi)外壁面的沿程壓力和馬赫數(shù)分布用1.加電Kgfit微叔國(guó)魏微線、值鋼圖1】算例1中心面及內(nèi)外壁面沿程壓力和馬赫數(shù)分布圖12算例2中心面及內(nèi)外壁面沿程壓力和馬赫數(shù)分布圖11右端為隔離段入口用算例1沿程壓力一直緩慢增加廳直至距離隔離段出口0.1m處竹中心

21、線壓力突升府內(nèi)外壁面相對(duì)變化較為緩慢針且外壁面的壓力波動(dòng)略小于內(nèi)壁面用馬赫數(shù)首次突然下降的位置對(duì)應(yīng)于靜壓上升的位置竹且可以明顯看出府中心線的馬赫數(shù)最高討外壁面次之針內(nèi)壁面最低用算例2的旋轉(zhuǎn)反壓高于算例1府所以等直隔離段內(nèi)反壓向上游傳播的距離遠(yuǎn)大于算例14且壓力峰值逐漸降低并趨于平緩討隔離段內(nèi)變?yōu)閬喡曀倭鲃?dòng)用2.3.2對(duì)漸擴(kuò)隔離段的影響圖13為算例3中3條貫穿隔離段的軸線的沿程壓力和馬赫數(shù)分布映y=0.20m為內(nèi)壁面上的一條軸線（Jy=0.21m為擴(kuò)張段入口處的中心軸線ijy=0.22m則相當(dāng)于擴(kuò)張段入口處的“外壁面，就由壓力圖可以看出面3條線兒乎重合從馬赫數(shù)圖可見(jiàn),y=0.21m處的馬赫數(shù)值高

22、于另外2條線面且尸=0.20m處#即內(nèi)壁面靠近擴(kuò)張段出口處有一段馬赫數(shù)趨勢(shì)相反仲圖14馬赫數(shù)圖的峰值所對(duì)應(yīng)的z坐標(biāo)則代表反壓前傳的“極限位置啊算例3、算例4、算例5的“極限位置”依次遠(yuǎn)離入口*即隨著反壓旋轉(zhuǎn)速度的降低府反壓前傳的距離減少目算例3、算例6、算例7的“極限位置”也依次遠(yuǎn)離入口府即隨著反壓的降低巾反壓前傳的距離也減少用但是府壓力對(duì)于極限位置的影響大于旋轉(zhuǎn)速度的影響用圖13算例3中心面及內(nèi)外壁面沿程壓力和馬赫數(shù)分布I_n_I_i】Rg岫蚣iLW心咖故株密侃伐iiw四您1潮m圖14漸擴(kuò)隔離段中心面及內(nèi)外壁面的沿程壓力和馬赫數(shù)分布在等直隔離段中竹激波基本沿著隔離段中心環(huán)面與內(nèi)外壁面對(duì)稱府且

23、在等直隔離段環(huán)面中心形成了一道正激波竹隨著反壓的增大激波與附面層的干擾作用增強(qiáng)麻附面層厚度增加用帶擴(kuò)張段的隔離段中的激波串則呈非對(duì)稱分布狀態(tài)祥并且在隔離段內(nèi)未形成正激波行但是在首道斜激波之后的與內(nèi)外壁面發(fā)生了不同程度的附面層分離府且附面層分離較為嚴(yán)重用3結(jié)論本文采用2種隔離段構(gòu)型研究了旋轉(zhuǎn)爆轟反壓在隔離段內(nèi)的流動(dòng)特性廳得出以下結(jié)論： 在等直隔離段中由隨著旋轉(zhuǎn)爆轟反壓的增加府旋轉(zhuǎn)爆轟波向上游傳播府在隔離段內(nèi)形成斜激波系小在首道激波的下游面附面層逐漸加厚府靜壓逐漸升高府而馬赫數(shù)逐漸降低用由于激波與附面層的相互作用府發(fā)生普朗特-邁耶反射府形成了激波串用 同樣在0.30MPa的反壓條件下®在

24、帶擴(kuò)張段的隔離段中反壓前傳的距離比等直隔離段要更靠近入口面因此等直隔離段抗擾動(dòng)的能力比帶擴(kuò)張段的隔離段更強(qiáng)用 在帶擴(kuò)張段的隔離段中竹反壓一直前傳至擴(kuò)張段內(nèi)R并基本穩(wěn)定在某一位置針且該位置與反壓的旋轉(zhuǎn)速度及大小有關(guān)用隨著反壓旋轉(zhuǎn)速度的減小祐反壓向上游的擾動(dòng)也減小0隨著反壓的增大府反壓最終穩(wěn)定的位置也越靠近入口處用改變壓力比改變速度引起的變化更明顯用參考文獻(xiàn)劉世杰林志勇&覃蔥B等.連續(xù)旋轉(zhuǎn)爆震波發(fā)動(dòng)機(jī)研究進(jìn)展J.E航導(dǎo)彈»i2010ii2:7D-75.LIUSbijieilLINZhiioD；TANHuillelal.Renircbprogiessofconlinuouirola

25、iydetonationwaveengineI.ManeuverableMis"sileif20101(2：70-75.(inChinese)王健平*1周蕊R武丹.連續(xù)旋轉(zhuǎn)爆轟發(fā)動(dòng)機(jī)的研究進(jìn)展】J.實(shí)驗(yàn)流體力學(xué)il2015i|29(4)：2-25.WANGJibnpin?ifZHOURu宿WUDid.Profreuolcontinuouslyrotatingdetontlionengineiesear<hJ.lournt!ofExperiment!inFluidMechinicsil2H5i|29(4)：1225.(inChinese)溫玉芬.吸氣式脈沖爆震發(fā)動(dòng)機(jī)進(jìn)氣道的流動(dòng)特

26、性研究D.南京:南京航空航天大學(xué)2011.WENYifen.Reseaichoniht(lowcharsclerislicsofIheinleliorairbreathngpulsedetonationengineD.Ninjins：NanjingUniiefjiiyofAeronauticstndAiiroaauiitiOOi1.(nChinese)4王超劉衛(wèi)東劉世杰*等.吸氣式連續(xù)旋轉(zhuǎn)爆震與來(lái)流相互作用Jj.航空學(xué)報(bào)料。16明7(5)：411-1418.WANGChaoitLIUWeidongilLIUShijt.Inleratiionofairbruihin|tonnuousro'

27、;.a.'ngde.onationwhnowI,Acl8Aeronau:cieAstroQau:caSinic»il201637(5):Hll-Hl8.(nChinese)5SUWeijiilJIYonanitCHENYun.Effecuofdymmicbickpreuureonpseudoshockoscillationinnrimjeline:isoilorI.JoninilofPfopullionindPowe甫2。1632(2：516-529.(下轉(zhuǎn)第59頁(yè))點(diǎn)法的攔截幾何竹且導(dǎo)彈過(guò)裁指令變化平緩砰對(duì)三點(diǎn)法全彈道數(shù)值仿真具有較強(qiáng)的實(shí)用價(jià)值明參考文獻(xiàn)1 宙虎民.導(dǎo)彈制導(dǎo)

28、與控制原理M.北京：國(guó)防工業(yè)出版社2006.LEIHumin.M$：tsuidinteandcontrol$y$lemM.Beijing：NalionilDefenseIndustryPreislf2OI6.(inChinese)朱衛(wèi)兵N張耀良.戰(zhàn)術(shù)導(dǎo)彈三點(diǎn)法遙控制導(dǎo)彈道方程的解析解J.彈箭與制導(dǎo)學(xué)報(bào)f2006if26(1)=245-250.ZHUWeibingilZHANGYao.inj.AnakLcsoulionol:heballiiiiceqiHliontotheremotecontrolledticlicilminilebjusingihrefpoin:im«:hodJ.Jo

29、urnalolProjectilesRockeliiiMissilesandGudancelt20l)626(1):245-250.(ioChinese)周慈鐘il李忠應(yīng).有翼導(dǎo)彈飛行動(dòng)力學(xué)M.北京:北京航空航天大學(xué)出版社R1987.ZHOUHuizhonji)LIZhongyng,Fgh:aerodinimcjo'wingedmisjileM.Beijing：BeijingUnifersilyofAeronautics&Ailroniniic!PiessHl987.(nChinese)張大元。趙玉芹.地空導(dǎo)彈三點(diǎn)法三維運(yùn)動(dòng)學(xué)彈道建模與仿真飛行力學(xué)il2H5il3l(l)：57

30、-62.ZHANGDtiutnitZHAOYuqin.ModelingandsimulkingonkinemjticIrijecoryofgroundtoa：rmin.leconlro.edbyihret-pointme:hodJ.FlightDjnamittil201Sit30(1)：57-62.(inChnese)5RATN00A.Threepointgu:dinceforinterceptingweavingtJ.JournalofGuidanceifControlindDymmiciit2016if39(S)：1879-I884.(inChinae)6朱戰(zhàn)霞韓沛打陳鵬.基于非線性Tem

31、ina滑模的動(dòng)能攔截器末制導(dǎo)律設(shè)計(jì)11.西北工業(yè)大學(xué)學(xué)報(bào)»l2013i31(2):233-238.ZHUZhimiiilHANPeWCHENPeng.Designoinonlineir'.ermiaalSMGL!orKKVJ.JournalofNorthwejiernPoly.echnica!Un:reisityil?(ll3il3<2)：233-238.(inChinese)7KUMARSRifRAOSifGH0SED.SlidingmodeguidanceandcontrolforaJ'aspulimeicepiorwiihteimini!i嘰kconstr

32、iinuJj.JoarnilolGii4inceiilControltndDynsmic(il2D12it35(4)f1230-1246.8 孫勝K周荻.有限時(shí)間收斂變結(jié)構(gòu)導(dǎo)引律門，宇航學(xué)報(bào)也。8。29(4)：1258-|262.SUNShen詢ZHOUD：.Afinite1metonver;emvaiiibiestrntture;udance1jwJ.lournilolAstronaoticsil2008il29(4)：1258_1262.(inChinese)螢飛春.導(dǎo)彈魯棒高階滑模制導(dǎo)控制一體化方法研究】.航空學(xué)報(bào)|20131|34(9)J212-2218.DONGFeiyao.Rtua

33、rtholn:egratedrobusthighorderslidingmodeguidanceandconirolformissilesJ.AtitAeronjulici&Astion8uticii|2(ll3b4(9>：22)2-22i8.(iiChinese)(上接第41頁(yè))6馬虎.壓力條件對(duì)旋轉(zhuǎn)爆震發(fā)動(dòng)機(jī)的影響J.彈道學(xué)報(bào)H2012»l24(4)：94-9S.MAHu.ESfectofpreure(ondi:ononroiaiin;deiontlionengineJ.JournalolBiLisiitii2l)l2i|24(4)：94-98,(inChinete)T李寶星X翁春生.進(jìn)氣總壓對(duì)連續(xù)旋轉(zhuǎn)爆表發(fā)動(dòng)機(jī)爆爰影響的二維數(shù)值模擬J.固體火箭技術(shù)20169(5):612-618.LIBaoiinHWENGChuntheng.Two4imensionilnumericalsimulationof:heinletstagnationpressureinfluenceon:hecoDlioDoailohtingdetooelionmineJ.JournalofSolidRocketTechnolofyii2«16«i39(5)：612-61S.(inChine”>邱玉杰R