• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    柵格翼氣動特性及其應(yīng)用研究綜述

    2015-04-22 05:42:32張為華
    固體火箭技術(shù) 2015年4期
    關(guān)鍵詞:超聲速柵格氣動

    彭 科,胡 凡,張為華,周 張

    (1.國防科學(xué)技術(shù)大學(xué) 航天科學(xué)與工程學(xué)院,長沙 410073;2.江南工業(yè)集團(tuán)有限公司,長沙 410200)

    ?

    柵格翼氣動特性及其應(yīng)用研究綜述

    彭 科1,胡 凡1,張為華1,周 張2

    (1.國防科學(xué)技術(shù)大學(xué) 航天科學(xué)與工程學(xué)院,長沙 410073;2.江南工業(yè)集團(tuán)有限公司,長沙 410200)

    從實(shí)驗(yàn)研究、工程估算、數(shù)值模擬三方面綜述了柵格翼氣動特性研究方法及研究現(xiàn)狀;概述了柵格翼主要?dú)鈩犹匦?;總結(jié)了弧形、整體后(前)掠、局部后掠等格翼新穎幾何構(gòu)型研究進(jìn)展;全面概述了柵格翼在國內(nèi)外各類型飛行器上應(yīng)用情況及主要用途。針對格翼設(shè)計(jì)應(yīng)用過程的瓶頸問題,分析展望了柵格翼相關(guān)領(lǐng)域需要關(guān)注和解決的關(guān)鍵問題。

    柵格翼;氣動特性;幾何構(gòu)型;應(yīng)用動態(tài)

    0 引言

    柵格翼是由外部框架和內(nèi)部柵格布置形成的空間多升力面系統(tǒng)[1],其基本構(gòu)型分為框架式與蜂窩式,蜂窩式又分為正置和斜置。柵格翼起源于飛行器設(shè)計(jì)初期所采用的由多翼面組成的升力系統(tǒng),但由于多翼面升力系統(tǒng)在當(dāng)時的技術(shù)條件下存在結(jié)構(gòu)、重量、阻力、工藝等方面的缺陷,長期以來飛行器設(shè)計(jì)中單翼升力面占主導(dǎo)地位。隨著技術(shù)的進(jìn)步,人們對起源于多翼面升力系統(tǒng)的柵格翼的性能特點(diǎn)認(rèn)識不斷深入,用其作為各類飛行器氣動穩(wěn)定面與控制面的優(yōu)勢不斷凸顯。俄羅斯、美國、德國和中國等都加強(qiáng)了柵格翼相關(guān)研究工作,取得了一系列理論與應(yīng)用成果。

    柵格翼幾何結(jié)構(gòu)復(fù)雜,氣動問題是其設(shè)計(jì)過程中的核心關(guān)鍵問題,對柵格翼氣動特性國內(nèi)外研究進(jìn)展進(jìn)行綜述,把握其研究方法、性能特點(diǎn)、應(yīng)用前景,具有十分重要的意義。本文綜述柵格翼氣動特性研究方法、研究成果與國內(nèi)外應(yīng)用動態(tài),展望柵格翼相關(guān)領(lǐng)域需要關(guān)注和解決的關(guān)鍵問題,為相關(guān)領(lǐng)域研究人員提供參考。

    1 柵格翼氣動特性研究方法進(jìn)展

    柵格翼氣動特性研究方法包含實(shí)驗(yàn)、工程估算以及數(shù)值計(jì)算三方面。

    1.1 實(shí)驗(yàn)研究

    柵格翼外形復(fù)雜,其外流場結(jié)構(gòu)、氣動特征、氣動分析方法與平面翼差異較大,為掌握柵格翼外流場特征及基本氣動特性,建立其氣動特性分析方法。國內(nèi)外從風(fēng)洞實(shí)驗(yàn)、水洞實(shí)驗(yàn)、彈道靶實(shí)驗(yàn)與飛行實(shí)驗(yàn)等方面廣泛開展了柵格翼氣動特性研究。

    美國從1985年至上世紀(jì)末密集開展了大量不同構(gòu)形柵格翼亞、跨、超聲速風(fēng)洞實(shí)驗(yàn),研究了柵格翼基本構(gòu)型氣動特性、弧形柵格翼氣動特性、跨聲速壅塞現(xiàn)象、減阻構(gòu)型、幾何參數(shù)對氣動特性影響、后掠效應(yīng)、平面翼與柵格翼性能異同等[2];1996年,美國進(jìn)行了柵格翼彈頭的火箭助推飛行試驗(yàn)與空投試驗(yàn),驗(yàn)證了柵格翼作為氣動穩(wěn)定部件的良好性能[2];近年來,美國進(jìn)行了“獵戶座”逃逸飛行器縮比模型亞、跨、超聲速氣動特性風(fēng)洞實(shí)驗(yàn)研究,驗(yàn)證了柵格翼布局逃逸飛行器靜穩(wěn)定性,對比了常規(guī)與整體后掠?xùn)鸥褚須鈩有阅躘3]。英國采用風(fēng)洞實(shí)驗(yàn)研究對比了柵格翼與平面翼縱橫向靜穩(wěn)定特性、控制特性、鉸鏈力矩特性與阻力特性等方面的異同。國外還開展了柵格翼彈頭彈道靶實(shí)驗(yàn)研究,測得了柵格翼彈頭滾轉(zhuǎn)力矩系數(shù)與滾轉(zhuǎn)阻尼系數(shù)[4-6],得到了流場結(jié)構(gòu)陰影圖像。

    國內(nèi),南京理工大學(xué)經(jīng)過多次實(shí)驗(yàn),研究了亞、跨、超聲速條件下格寬翼弦比對柵格翼氣動特性的影響;探索了減小柵格翼超聲速阻力的方法和途徑;分別研究了采用無尾翼、“T”型尾翼、柵格尾翼時鴨式布局遠(yuǎn)程彈在跨/超聲速段的氣動特性,得到了不同馬赫數(shù)及攻角狀態(tài)下不同尾翼構(gòu)形對遠(yuǎn)程彈氣動特性尤其是鴨舵滾轉(zhuǎn)控制特性的影響[7-11]。北京航空航天大學(xué)陸中榮教授等采用氫氣泡法和絲線法在水洞和低速風(fēng)洞中進(jìn)行了柵格翼繞流特性實(shí)驗(yàn)[12],根據(jù)實(shí)驗(yàn)結(jié)果分析了柵格翼自由渦系與各網(wǎng)孔內(nèi)流態(tài),分析了流動機(jī)理。中國空氣動力研究與發(fā)展中心建立了國內(nèi)首個全尺寸柵格舵試驗(yàn)平臺[13],實(shí)現(xiàn)了全尺寸柵格舵模型氣動外形、舵控系統(tǒng)、結(jié)構(gòu)的一體化試驗(yàn)驗(yàn)證;此外,該中心還成功研制出了動態(tài)響應(yīng)能力良好、滿足小量鉸鏈力矩測量要求的天平[14]。中國航天空氣動力技術(shù)研究院研究了柵格翼展開過程氣動力測量試驗(yàn)方法[15],實(shí)現(xiàn)了全尺寸柵格翼超聲速條件下自動連續(xù)展開。

    實(shí)驗(yàn)手段加深了人們對柵格翼繞流特性與氣動特性的理解,研究成果促進(jìn)了工程估算與數(shù)值計(jì)算的不斷發(fā)展,為柵格翼設(shè)計(jì)提供有力支撐。

    1.2 工程估算

    以大量實(shí)驗(yàn)現(xiàn)象、實(shí)驗(yàn)數(shù)據(jù)為基礎(chǔ),國內(nèi)外學(xué)者結(jié)合不同馬赫數(shù)狀態(tài)柵格翼流場結(jié)構(gòu),根據(jù)部件組拆法,考慮各部件間干擾建立了柵格翼翼身組合體氣動特性工程估算模型。

    1.2.1 柵格翼流場結(jié)構(gòu)

    柵格翼幾何結(jié)構(gòu)復(fù)雜,其繞流結(jié)構(gòu)隨馬赫數(shù)變化而發(fā)生演變,對其氣動特性影響顯著,需根據(jù)不同馬赫數(shù)條件下流場結(jié)構(gòu)特點(diǎn)確定相應(yīng)氣動特性工程估算方法[2]。引入臨界馬赫數(shù)表征柵格翼流場結(jié)構(gòu),不同馬赫數(shù)狀態(tài)柵格翼流場結(jié)構(gòu)與流動特征如表1所示。

    表1 柵格翼流場結(jié)構(gòu)與流動特征

    Table 1 Grid fin flow field Sstructure and characteristics

    各臨界馬赫數(shù)的計(jì)算是確定特定馬赫數(shù)下柵格翼流場結(jié)構(gòu)的關(guān)鍵。文獻(xiàn)[2]將柵格翼內(nèi)部流動等效為二維噴管流動建立了柵格翼第一臨界馬赫數(shù)計(jì)算方法。第二臨界馬赫數(shù)Macr2與格片前緣半徑等參數(shù)相關(guān),計(jì)算較為復(fù)雜,目前多依靠實(shí)驗(yàn)數(shù)據(jù)比擬得到。第三臨界馬赫數(shù)可根據(jù)格片鍥角采用理論公式計(jì)算,當(dāng)格片為鈍前緣時采用等效鍥方法計(jì)算,精度要求較高時可采用數(shù)值方法計(jì)算[16]。

    1.2.2 單獨(dú)柵格翼氣動特性工程估算模型

    對于單獨(dú)柵格翼,其氣動力包含表面法向壓力與切向摩擦力兩部分貢獻(xiàn)。摩擦力貢獻(xiàn)目前多采用經(jīng)驗(yàn)公式計(jì)算[17-18];表面法向壓力貢獻(xiàn)采用基于無粘假設(shè)的工程模型計(jì)算,本節(jié)對該領(lǐng)域國內(nèi)外研究方法進(jìn)行綜述。

    亞聲速時,根據(jù)實(shí)驗(yàn)結(jié)果,由柵格翼格片脫出的自由渦系形成互相平行的旋渦壁[12,19],該類流動的氣動力計(jì)算適于采用基于升力面理論的渦格法。文獻(xiàn)[18,20-23]采用該方法計(jì)算了不同類型柵格翼亞聲速氣動特性,文獻(xiàn)[18]對馬赫數(shù)影響作了修正。此外,臧勇基于第二代低階面元法VSAERO,建立了亞聲速、小攻角狀態(tài)飛行器流場模型與氣動載荷計(jì)算模型,進(jìn)而計(jì)算了柵格翼亞聲速氣動特性[24]。

    柵格翼跨聲速段格孔內(nèi)部流動復(fù)雜,氣動特性工程估算難度較大,目前國內(nèi)外相關(guān)研究較少,計(jì)算理論尚不成熟。Macr1≤Ma≤1情況下多采用考慮邊界層修正的一維流理論計(jì)算[25],在1≤Ma≤Macr2根據(jù)理論公式計(jì)算脫體激波前后氣流參數(shù),近似認(rèn)為波后流動方向不變,根據(jù)波后參數(shù)采用亞聲速狀態(tài)方法計(jì)算柵格舵所受氣動力。

    關(guān)于柵格翼超聲速氣動特性工程估算,沈遐齡等不考慮格片間激波干擾,根據(jù)各格片理論升力線斜率求得正置柵格翼超聲速氣動特性[20]。文獻(xiàn)[26-27]給出了斜置蜂窩式柵格翼對正置蜂窩式柵格翼的等效條件。文獻(xiàn)[25]根據(jù)平面翼超聲速氣動特性計(jì)算模型,考慮端板效應(yīng),求得柵格翼各格片表面壓力分布,進(jìn)而得到其超聲速氣動特性。Theerthamalai P等考慮格片間激波干擾,基于激波-膨脹波理論發(fā)展了較為完整的柵格翼超聲速氣動特性工程估算方法[17]。文獻(xiàn)[28]基于面元法建立柵格翼翼身組合體高超聲速氣動特性計(jì)算模型。

    1.2.3 柵格翼翼身組合體部件干擾模型

    對于柵格翼翼身組合體,需考慮彈身對柵格翼干擾、柵格翼對彈身干擾、柵格翼各格片間干擾。彈身對柵格翼的氣動干擾計(jì)算方法常用的方法有:偶極子模擬彈身法、上洗流場模型法、平面渦格法,其中偶極子模擬彈身法所做的簡化最少,精度最高,是用來計(jì)算柵格翼與彈身間相互氣動干擾的最佳方法;上洗流場模型法、平面渦格法計(jì)算較簡單,在對計(jì)算精度要求不高的情況下可采用[18]。對大攻角狀態(tài)彈體脫體渦影響,目前國內(nèi)外多采用經(jīng)驗(yàn)公式估算[18,22]。柵格翼對彈身干擾主要基于鏡像法與經(jīng)驗(yàn)公式計(jì)算[17]。柵格翼各格片間干擾在亞聲速狀態(tài)需做考慮,文獻(xiàn)[25]采用1.2.2節(jié)方法將各格片同時劃分面元聯(lián)合求解,所得結(jié)果即反應(yīng)了柵格舵舵片間干擾,跨、超聲速段小型運(yùn)載火箭各舵片間干擾可忽略不計(jì)。

    現(xiàn)有工程估算方法效率較高,但精度不能滿足柵格翼外形設(shè)計(jì)與氣動特性仿真等工程需要,相關(guān)理論方法尚待進(jìn)一步研究。

    1.3 數(shù)值模擬

    基于大量假設(shè)與簡化的工程估算方法存在精度瓶頸,數(shù)值模擬方法被引入到柵格翼氣動特性研究。該方法能用較少假設(shè)給出詳盡、準(zhǔn)確、廣泛的信息,在柵格翼氣動問題研究過程中發(fā)揮了不可替代的作用,拓展了柵格翼氣動特性研究方式。

    關(guān)于柵格翼氣動特性數(shù)值模擬,國內(nèi)外公開發(fā)表的文獻(xiàn)所取數(shù)學(xué)模型、求解方法多為常規(guī)CFD方法。其幾何結(jié)構(gòu)復(fù)雜,數(shù)值模擬難度主要集中于網(wǎng)格劃分。目前國內(nèi)外柵格翼翼身組合體氣動特性數(shù)值模擬網(wǎng)格劃分方法主要包括:多塊對接結(jié)構(gòu)網(wǎng)格[29]、四面體非結(jié)構(gòu)網(wǎng)格[30]、直角非結(jié)構(gòu)網(wǎng)格[31]、結(jié)構(gòu)/非結(jié)構(gòu)混合網(wǎng)格[32-33]重疊結(jié)構(gòu)網(wǎng)格[34]等,后3種方法最為常用。直角非結(jié)構(gòu)網(wǎng)格劃分方法簡單、省時、自動化程度較高,便于網(wǎng)格自適應(yīng)處理;結(jié)構(gòu)/非結(jié)構(gòu)混合網(wǎng)格方法結(jié)合了結(jié)構(gòu)網(wǎng)格和非結(jié)構(gòu)網(wǎng)格各自的優(yōu)點(diǎn),減小網(wǎng)格生成難度,提高網(wǎng)格質(zhì)量,計(jì)算精度較高;重疊結(jié)構(gòu)網(wǎng)格劃分方法可有效提高計(jì)算精度。

    數(shù)值模擬方法對于柵格翼的氣動特性預(yù)示精度可滿足工程需求,但計(jì)算效率較低。

    2 柵格翼氣動特性研究成果綜述

    2.1 柵格翼主要?dú)鈩犹匦?/p>

    柵格翼作為新型氣動部件,在升力特性、鉸鏈力矩特性、折疊特性、強(qiáng)度特性等方面較傳統(tǒng)平面翼具有不可替代的優(yōu)點(diǎn),主要性能特點(diǎn)為[38]:

    (1)升力特性好,大攻角狀態(tài)仍能保持良好線性,不易失速;

    (2)壓心漂移小,可有效減小鉸鏈力矩,降低對舵機(jī)的要求,減輕舵機(jī)質(zhì)量;

    (3)弦向尺寸小,可緊貼彈體折疊安裝,減小尺寸,利于存儲,折疊狀態(tài)可在自身空氣動力矩作用下自動打開;

    (4)強(qiáng)度質(zhì)量比大,有利于降低結(jié)構(gòu)質(zhì)量;

    (5)阻力較大,可用作減速部件,合理設(shè)計(jì)降低阻力后可拓寬使用范圍。

    此外,柵格翼還有以下特點(diǎn):

    (1)作為空間多升力面系統(tǒng),柵格尾翼縱向格片對升力仍有貢獻(xiàn),比平板尾翼利于提高全箭靜穩(wěn)定性,可提供更大法向過載,提高機(jī)動性能。值得關(guān)注的是,背風(fēng)區(qū)柵格翼大攻角時其根部出現(xiàn)流動分離,導(dǎo)致部分升力面當(dāng)?shù)毓ソ菫樨?fù)值,對升力的貢獻(xiàn)低于其他翼面[30];

    (2)對于采用鴨舵控制滾轉(zhuǎn)的火箭/導(dǎo)彈,在攻角非零情況下,受脫體渦影響,尾翼為平面翼時易出現(xiàn)滾轉(zhuǎn)控制耦合問題,甚至產(chǎn)生控制反效,采用柵格翼時滾轉(zhuǎn)控制更為有效,為解決鴨式布局滾轉(zhuǎn)耦合問題提供了思路[39];

    (3)在大攻角狀態(tài)柵格尾翼仍可保持偏航方向穩(wěn)定性,也具備大攻角狀態(tài)橫向機(jī)動的能力;

    (4)柵格翼跨聲速階段產(chǎn)生氣流壅塞現(xiàn)象[40-41],對其氣動特性影響顯著,控制效率降低[38],阻力增加,且第一臨界馬赫數(shù)附近法向力系數(shù)、俯仰力矩系數(shù)、壓心系數(shù)變化皆較大,對飛行器靜穩(wěn)定性產(chǎn)生較大影響[42-43]。

    2.2 柵格翼幾何構(gòu)型研究進(jìn)展

    對于柵格翼基本構(gòu)形,影響其氣動性能的主要因素有邊框厚度、格片厚度、柵格長寬比、格寬翼弦比與翼型等。陳少松等采用風(fēng)洞實(shí)驗(yàn)方法對不同形狀柵格翼的減阻特性進(jìn)行了分析,得到了柵格翼邊框尺寸、邊框截面形狀、格柵形式、格寬翼弦比等對柵格翼阻力的影響,結(jié)果表明,(1)柵格翼的邊框剖面形狀對柵格翼的阻力影響較大,遠(yuǎn)大于格片厚度帶來的影響;(2)減少柵格數(shù)目的情況下,柵格翼阻力降低,同時帶來升力的損失;(3)柵格翼格寬翼弦比通過影響內(nèi)部流態(tài)影響其升阻力特性,在一定馬赫數(shù)、攻角狀態(tài)下,存在最佳格寬翼弦比使得柵格翼升阻比最大[8-10]。吳曉軍等運(yùn)用數(shù)值模擬方法研究了柵格長寬比對柵格氣動特性的影響,結(jié)果表明,在所選擇的長寬比范圍內(nèi),長寬比對阻力的影響較小,而對法向力影響明顯,可根據(jù)需要優(yōu)化長寬比[44]。

    隨著研究的深入,出現(xiàn)了弧形柵格翼、整體后掠?xùn)鸥褚砼c前緣局部后掠?xùn)鸥褚淼确浅R?guī)構(gòu)形柵格翼,提升了性能。

    2.2.1 弧形柵格翼

    文獻(xiàn)[45]介紹了弧形柵格翼概念,弧形柵格翼折疊安裝可與彈體完全貼合,最大限度減小尺寸?;⌒螙鸥褚砜稍O(shè)計(jì)為凸面迎風(fēng)與凹面迎風(fēng)兩種形式。兩種形式的弧形柵格翼在跨聲速段升阻比與平面柵格翼差異不大,亞、超聲速段皆高于平面柵格翼[46];迎風(fēng)面為凸面時其升力系數(shù)與升阻比皆略大于迎風(fēng)面為凹面的情況;迎風(fēng)面為凹面時弧形柵格翼展開所需機(jī)構(gòu)力與其所受氣動力方向一致,便于展開。目前弧形柵格翼設(shè)計(jì)多以凹面為迎風(fēng)面(圖1)。

    圖1 弧形柵格翼Fig.1 Curved grid fin

    2.2.2 整體后(前)掠?xùn)鸥褚?/p>

    文獻(xiàn)[45]引出了整體后掠與整體前掠?xùn)鸥褚砀拍?,各自又分為“斜”與“平”2種類型(圖2、圖3)。

    “斜型”柵格翼格片平面與彈體軸線夾角δ等于翼面后掠角Λ,來流攻角為零時格片當(dāng)?shù)毓ソ遣粸榱?。對于“斜型”柵格翼,?dāng)后(前)掠角一致時,后掠與前掠方式升、阻力特性基本一致,在前/后掠角在約±30°以內(nèi)時,法向力基本不隨后(前)掠角變化,30°后(前)掠角狀態(tài)阻力約為零后(前)掠角狀態(tài)的2~3倍,“斜型”后掠?xùn)鸥褚砜稍诒3诛w行器靜穩(wěn)定性的同時用作減速部件[45]。“平型”后(前)掠?xùn)鸥褚砜捎行p少阻力,后掠狀態(tài)較前掠狀態(tài)減阻效果更好?!矮C戶座”逃逸飛行器采用“平型”整體后掠?xùn)鸥褚硖岣哽o穩(wěn)定性能[3]。

    (a)斜型 (b)平型

    (a)斜型 (b)平型

    2.2.3 前緣局部后掠?xùn)鸥褚?/p>

    受平面翼后掠減阻效應(yīng)啟發(fā),Guyot等提出了局部后掠?xùn)鸥褚砀拍頪47](圖4),其零升阻力較常規(guī)構(gòu)型柵格翼約減少30%~40%,比整體后掠?xùn)鸥褚頊p阻效果更好,并可有效減緩跨聲速段壅塞現(xiàn)象,減小阻力[48-49]。局部后掠?xùn)鸥褚矸譃闁鸥窠唤狱c(diǎn)后掠(模型GP)與格片中心點(diǎn)后掠兩類(模型GV),前者減阻更為有效[11]。

    前緣局部后掠?xùn)鸥褚硪话阋硐逸^長,一定程度上影響了其折疊性能,另外加工更為復(fù)雜,目前尚未見前緣局部后掠?xùn)鸥褚淼南嚓P(guān)工程應(yīng)用。

    (a)常規(guī)構(gòu)型

    (b)后掠?xùn)鸥褚?/p>

    3 柵格翼國內(nèi)外應(yīng)用動態(tài)

    柵格翼因其獨(dú)特優(yōu)勢引起了世界各國的重視,前蘇聯(lián)/俄羅斯、美國、德國以及中國等都加強(qiáng)了柵格翼的理論與試驗(yàn)研究工作,將其成功用于各類飛行器[50]。

    3.1 前蘇聯(lián)/俄羅斯

    前蘇聯(lián)于20世紀(jì)50年代初就對柵格翼進(jìn)行了研究,進(jìn)行了大量氣動計(jì)算和實(shí)驗(yàn),開展了柵格翼結(jié)構(gòu)、強(qiáng)度、生產(chǎn)工藝的研究,形成了柵格翼設(shè)計(jì)方法,并成功將柵格翼用作聯(lián)盟號宇宙飛船救生逃逸系統(tǒng)的穩(wěn)定、減速部件,1963年實(shí)現(xiàn)了首次發(fā)射。N-1探月火箭一級采用柵格翼作為穩(wěn)定面,在轉(zhuǎn)運(yùn)過程中柵格翼折疊,1969年進(jìn)行了首次發(fā)射。

    70年代前后,前蘇聯(lián)科學(xué)家開始在各類導(dǎo)彈設(shè)計(jì)中廣泛應(yīng)用柵格翼。戰(zhàn)略彈道導(dǎo)彈中,野人(Savage,北約代號SS-13)I型與II型一級使用4個柵格翼作為穩(wěn)定面;罪人(Sinner,北約代號SS-16)、佩刀(Saber,北約代號SS-20)I型與II型、鐮刀(Sickle,北約代號SS-25)、鐮刀 B(Sickle B,北約代號SS-27)一級尾段有8個柵格翼,其中4個為穩(wěn)定面,另外4個為控制面。

    戰(zhàn)術(shù)彈道導(dǎo)彈中,薄板(Scaleboard,北約代號SS-12,兩級固體導(dǎo)彈)二級尾部安裝有4個弧形柵格翼,一級飛行時柵格翼折疊,以減小一級飛行靜不穩(wěn)定度,改善一級操縱性,二級飛行時柵格翼展開,并作為控制面與燃?xì)舛媛?lián)動;圣甲蟲(Scarab,北約代號SS-21)、蜘蛛(Spider,北約代號SS-23)采用可折疊柵格翼,作為控制面與燃?xì)舛媛?lián)動。

    空空導(dǎo)彈中蝰蛇(Adder,北約代號AA-12,俄羅斯稱R-77)及其各改進(jìn)型采用可折疊柵格翼作為控制面,R-77導(dǎo)彈是世界上首型采用柵格翼作為控制面的空空導(dǎo)彈,充分利用了柵格翼各方面優(yōu)點(diǎn),是目前柵格翼最成功的應(yīng)用案例[50]。對于空射導(dǎo)彈,隨著現(xiàn)代戰(zhàn)爭對隱身性能的需求不斷增強(qiáng),機(jī)載武器、傳感器等掛載實(shí)現(xiàn)內(nèi)埋是該類武器系統(tǒng)發(fā)展的必然趨勢,使用可折疊的柵格翼是減小空射導(dǎo)彈體積、實(shí)現(xiàn)內(nèi)埋的有效手段[51]。

    俄羅斯最新的潛射巡航導(dǎo)彈熾熱(Sizzler,北約代號SS-N-27,俄羅斯稱3M-54 Klub)助推級采用折疊柵格翼,待發(fā)狀態(tài)緊貼彈體折疊,發(fā)射后柵格翼打開,可有效解決導(dǎo)彈水中潛行靜穩(wěn)定問題。

    3.2 美國

    美國柵格翼相關(guān)研究始于20世紀(jì)80年代中后期,1991年產(chǎn)生了第一個關(guān)于柵格翼技術(shù)的專利[51]。美國在MOAB(Massive Ordnance Air Blast bomb)炸彈上實(shí)現(xiàn)了柵格翼的首次應(yīng)用,將柵格翼用作穩(wěn)定面與控制面,MOAB炸彈于2003年3月11日開展了首次試驗(yàn)。此外,美國在MOP(Massive Ordnance Penetrator)侵徹彈中同樣使用了柵格翼。美國空軍干擾彈MALD(Miniature Air-Launched Decoy)采用復(fù)雜構(gòu)形可折疊柵格翼,可從B-52轟炸機(jī)、F-16戰(zhàn)斗機(jī)等平臺發(fā)射以干擾對方防空力量[52]。為提高反潛作戰(zhàn)能力,解決聲吶浮標(biāo)的高空準(zhǔn)確投放問題,美國于2005年開展了相關(guān)技術(shù)研究,設(shè)計(jì)了空投飛行器(ADV,Air Delivery Vehicle),該飛行器在高空從飛機(jī)上投放后,可在柵格翼控制下機(jī)動到指定位置并釋放聲吶浮標(biāo),該技術(shù)也可用于其他作戰(zhàn)環(huán)境下的載荷空投。此外,美國在“獵戶座”逃逸飛行器[3]及無逃逸塔的逃逸飛行器方案[53]中同樣使用了柵格翼。

    3.3 其他

    德國也開展了柵格翼的理論研究和試驗(yàn)研制工作。2002年2月,德國試飛了其首款使用柵格翼進(jìn)行姿態(tài)控制的高超聲速導(dǎo)彈HFK-E0,其飛行馬赫數(shù)達(dá)6.66。

    國內(nèi)自20世紀(jì)90年代初開始了對柵格翼的研究,“長征二號F”型運(yùn)載火箭逃逸飛行器使用柵格翼作為氣動穩(wěn)定部件,快舟小型固體運(yùn)載火箭在國內(nèi)實(shí)現(xiàn)了柵格舵作為氣動控制部件的首次應(yīng)用,標(biāo)志著我國柵格翼設(shè)計(jì)應(yīng)用水平邁上了新的臺階。

    4 關(guān)鍵技術(shù)問題分析與展望

    在柵格翼獨(dú)特性能優(yōu)勢的驅(qū)動下,國內(nèi)外關(guān)于其氣動特性研究與應(yīng)用實(shí)踐方面都取得一系列的進(jìn)展與突破,體現(xiàn)在以下方面:

    (1)氣動特性研究手段不斷發(fā)展。實(shí)驗(yàn)研究投入不斷增大,為準(zhǔn)確把握其流場結(jié)構(gòu)與性能特點(diǎn)發(fā)揮了不可替代的作用;初步建立了柵格翼氣動特性工程估算方法;采用數(shù)值求解方法在其性能與新穎構(gòu)型研究方面取得諸多成果。

    (2)對于其氣動特性的認(rèn)識不斷深入,研究提出了弧形、全局后掠、局部后掠等新穎柵格翼構(gòu)型,為進(jìn)一步提高柵格翼性能提供了新的思路。

    (3)在多個國家、不同類型飛行器、多種型號上得以成功應(yīng)用。

    為進(jìn)一步發(fā)揮柵格翼作為氣動穩(wěn)定面與控制面的性能優(yōu)勢,克服其缺點(diǎn),以下方面值得進(jìn)一步關(guān)注與研究:

    (1)柵格翼氣動特性工程計(jì)算方法研究。工程計(jì)算方法效率較高,是柵格翼初步設(shè)計(jì)過程不可或缺的仿真手段?;诖罅繉?shí)驗(yàn)與數(shù)值仿真結(jié)果修正現(xiàn)有工程計(jì)算模型,建立其氣動特性計(jì)算經(jīng)驗(yàn)?zāi)P?,是亟待開展的工作。

    (2)柵格翼外形設(shè)計(jì)準(zhǔn)則研究。目前公開發(fā)表的文獻(xiàn)鮮有直接面向柵格翼設(shè)計(jì)的成果,亟需建立柵格翼格數(shù)、格柵形式、格寬翼弦比等基本參數(shù)設(shè)計(jì)準(zhǔn)則,加強(qiáng)柵格翼設(shè)計(jì)模型、流程、方法研究。

    (3)跨聲速段動態(tài)不穩(wěn)定性研究。柵格翼在跨聲速段特定馬赫數(shù)條件下俯仰力矩距離變化,動態(tài)不穩(wěn)定問題突出,采用非定常數(shù)值模擬等手段對其產(chǎn)生過程、產(chǎn)生機(jī)理進(jìn)行研究,對降低柵格翼應(yīng)用成本與應(yīng)用風(fēng)險具有現(xiàn)實(shí)意義。

    (4)隱身性能研究。柵格翼隱身性能較差,開展其隱身特性定量研究對其應(yīng)用具有指導(dǎo)意義。

    (5)高超聲速狀態(tài)氣動熱特性研究、氣動力/熱耦合計(jì)算研究、防熱結(jié)構(gòu)設(shè)計(jì)。柵格翼優(yōu)良的高超聲速性能必將帶來其高超聲速領(lǐng)域應(yīng)用的不斷拓展,但高超聲速條件下其復(fù)雜繞流結(jié)構(gòu)導(dǎo)致的氣動熱問題同樣突出,開展高超聲速狀態(tài)氣動熱特性研究、氣動力/熱耦合計(jì)算研究與防熱結(jié)構(gòu)設(shè)計(jì)研究重要理論與實(shí)踐意義。

    (6)結(jié)構(gòu)與加工工藝研究。柵格翼加工工藝較為復(fù)雜,高超聲速條件下防熱結(jié)構(gòu)的引入使得該問題更為突出,在當(dāng)前最新工藝成果的基礎(chǔ)上開展其加工工藝研究對推動?xùn)鸥褚響?yīng)用領(lǐng)域的拓展具有重要意義。

    5 結(jié)束語

    隨著相關(guān)領(lǐng)域技術(shù)的不斷發(fā)展,柵格翼作為氣動穩(wěn)定面與控制面的獨(dú)特性能優(yōu)勢必將使其應(yīng)用領(lǐng)域不斷拓展。結(jié)合運(yùn)載火箭、彈道導(dǎo)彈、巡航導(dǎo)彈、智能彈藥、再入飛行器、水下航行體等多種不同對象,開展柵格翼應(yīng)用研究必將產(chǎn)生巨大社會與軍事效益。

    [1] 吳甲生,雷娟棉.制導(dǎo)兵器氣動布局與氣動特性[M].北京:國防工業(yè)出版社,2008.

    [2] Wm David Washington,Mark S Miller.Experimental investigations of grid fin aerodynamics:a synopsis of nine wind tunnel and three flight tests[C]//RTO AVT Symposium on Missile Aerodynamics.Sorrento,Italy:1998.

    [3] Daniel A Pruzan,Michael R Mendenha,William C Rose,David M Schuster.Grid fin stabilization of the orion launch abort vehicle[C]//29th AIAA Applied Aerodynamics Conference.Honolulu,Hawaii:2011.

    [4] Dupuis A,Berner C.Aerodynamic aspects of a grid finned projectile at subsonic and supersonic velocities[C]//19th International Symposium of Ballistics.Interlaken,Switzerland:2001.

    [5] Berner C,Dupuis A.Wind tunnel tests of a grid finned projectile configuration[C]//39th Aerospace Sciences Meeting & Exhibit.Reno,Nevada:2001.

    [6] Fournier E Y.Wind tunnel investigation of grid fin and conventional planar control surfaces[C]//39th Aerospace Sciences Meeting & Exhibit.Reno,Nevada:2001.

    [7] 鄧帆,陳少松,譚獻(xiàn)忠,等.不同尾翼鴨式布局遠(yuǎn)程彈在跨/超聲速的氣動特性實(shí)驗(yàn)研究[J].實(shí)驗(yàn)流體力學(xué),2010,24(2):46-50.

    [8] 陳少松,張建葉.格寬翼弦比對柵格翼氣動特性的影響[J].南京理工大學(xué)學(xué)報(bào),2005,29(6):693-696.

    [9] 陳少松,徐琴,王福華,等.格柵翼減阻特性研究[J].流體力學(xué)實(shí)驗(yàn)與測量, 2001,15(4):7-11.

    [10] 陳少松,徐琴,王福華,等.格柵翼組合體的超音速氣動特性研究[J].彈道學(xué)報(bào), 2000,12(2):50-54.

    [11] 鄧帆,陳少松.柵格翼外形特征對減阻影響的研究[J].實(shí)驗(yàn)流體力學(xué),2011,25(3):10-15.

    [12] 陸中榮,王海文.柵格翼繞流特性的實(shí)驗(yàn)研究[J].流體力學(xué)實(shí)驗(yàn)與測量,2002,16(1):21-26.

    [13] 陳建中,趙忠良,涂正光,等.柵格舵氣動與操縱特性高速風(fēng)洞試驗(yàn)技術(shù)研究[J].實(shí)驗(yàn)流體力學(xué),2011,25(6):82-87.

    [14] 陳建中,易國慶,彭超,等.全尺寸柵格舵鉸鏈力矩天平研制[J].實(shí)驗(yàn)流體力學(xué), 2012,26(1):79-82.

    [15] 吳凱,張尚彬,蔣坤,等.柵格翼展開過程測力試驗(yàn)方法[C]//中國空氣動力學(xué)會測控技術(shù)專委會第六屆四次學(xué)術(shù)交流會文集.三亞:2013.

    [16] 張亮,王淑華,姜貴慶.鈍化前緣對柵格翼激波干擾與熱流分布的影響[J].宇航學(xué)報(bào),2010,31(2):324-328.

    [17] Theerthamalai P,Nagarathinam M.Aerodynamic analysis of grid-fin configurations at supersonic speeds[J].Journal of Spacecraft and Rockets,2006,43(4):750-756.

    [18] Theerthamalai P.Aerodynamic characterization of grid fins at subsonic speeds[J].Journal of Aircraft,2007,44(2):694-697.

    [19] LU Zhong-rong,SU Wen-han.Hairpin vortices in the wake of grid fin[J].Chinese Journal of Aeronautics,2001,14(4):200-204.

    [20] 沈遐齡,王麗麗.柵格翼和機(jī)身組合體的氣動特性計(jì)算分析[J].北京航空航天大學(xué)學(xué)報(bào),2001,27(1):54-56.

    [21] 余永剛.亞跨音速柵格翼氣動特性研究[D].南京:南京理工大學(xué),2007.

    [22] John E Burkhalter,Roy J Hartfield,Todd M Leleuxt.Nonlinear aerodynamic analysis of grid fin configurations[J].Journal of Aircraft,1995,32(3):547-554.

    [23] Ross A Brooks,John E Burkhaltert.Experimental and analytical analysis of grid fin configurations[J].Journal of Aircraft,1989,26(9):885-887.

    [24] 臧勇.柵格舵/彈身組合體亞音速氣動特性研究[D].北京:北京理工大學(xué),2000.

    [25] Richard W Kretzschmar,John E Burkhalter.Aerodynamic prediction methodology for grid fins[C]//RTO AVT Symposium on Missile Aerodynamics.Sorrento,Italy:1998.

    [26] 任懷宇,蔡昱,孫洪森.高超聲速柵格翼氣動特性研究[J].宇航學(xué)報(bào),2010,31(1):44-49.

    [27] 陸中榮,沈遐齡,童自力,等.柵格翼空氣動力特性的計(jì)算與分析[J].北京航空航天大學(xué)學(xué)報(bào),1996,22(5):575-580.

    [28] 周張,胡凡,于勇.基于面元法的柵格翼翼身組合體高超聲速氣動特性研究[J].固體火箭技術(shù),2012,35(4):442-446.

    [29] Sun Y,Khalid M.A Cfd investigation of grid fin missiles[C]//34th AIAAlASME/SAE/ASEE Joint Propulsion Conference & Exhibit.Cleveland,OH:1998.

    [30] James Despirito,Jr Milton E Vaughn,W David Washington.Cfd investigation of canard-controlled missile with planar and grid fins in supersonic flow[C]//AIAA Atmospheric Flight Mechanics Conference and Exhibit.Monterey,California:2002.

    [31] James E Kless,Michael J Aftosmis.Analysis of grid fins for launch abort vehicle using a cartesian euler solver[C]//29th AIAA Applied Aerodynamics Conference.Honolulu,Hawaii:2011.

    [32] Deng You-qi,Ma Ming-sheng,Zheng Ming.Navier-stokes computation of grid fin missile using hybrid structured-unstructured grids[J].Chinese Journal of Aeronautics,2006,19(4):304-308.

    [33] 劉剛,肖中云,江雄,等.混合網(wǎng)格方法在柵格翼數(shù)值模擬中的應(yīng)用研究[J].空氣動力學(xué)學(xué)報(bào),2007,25(2):230-277.

    [34] 王學(xué)占,王立強(qiáng).基于結(jié)構(gòu)重疊網(wǎng)格的柵格翼導(dǎo)彈流場數(shù)值仿真[J].航空兵器,2012(1):44-47.

    [35] Montgomery C Hughson,Eric L Blades,Edward A Luke,Gregg L Abate.Analysis of lattice grid tailfin missiles in high-speed flow[C]//25th AIAA Applied Aerodynamics Conference.Miami,FL:2007.

    [36] Philippe Reynier,Uwe Reisch,José-Maria Longo,Rolf Radespiel.Flow predictions around a missile with lattice wings using the actuator disc concept[J].Aerospace Science and Technology,2004(8):377-389.

    [37] 李錦,江定武,毛枚良,等.BGK-NS格式在復(fù)雜流動中的應(yīng)用研究[J].空氣動力學(xué)報(bào),2013,31(4):449-453.

    [38] Simpson G M,Sadler A J.Controls:A comparison with conventional,planar fins[C]//RTO AVT Symposium.Sorrento,Italy:1998.

    [39] James Despirito,Jr Milton E Vaughn,W David Washington.Numerical investigation of canard-controlled missile with planar and grid fins[J].Journal of Spacecraft and Rockets,2003,40(3):363-370.

    [40] Montgomery C Hughson,Eric L Blades,Gregg L Abate.Transonic aerodynamic analysis of lattice grid tail fin missiles[C]//24th Applied Aerodynamics Conference.San Francisco,California:2006.

    [41] Karl S Orthner.Aerodynamic analysis of lattice grid fins in transonic flow[D].Ohio:Air Force Institute Of Technology,2004.

    [42] Gregg Abate,Gerald Winchenbach,Wayne Hathaway.Transonic aerodynamic and scaling issues for lattice fin projectiles tested in a ballistics range[C]//19th International Symposium of Ballistics.Interlaken,Switzerland:2001.

    [43] Gregg Abate,Ralf Duckerschein,Wayne Hathaway.Subsonic/transonic free-flight tests of a generic missile with grid fins[C]//38th AIAA Aerospace Sciences Meeting and Exhibit.Reno,NV:2000.

    [44] 吳曉軍,陳紅全,鄧有奇,等.柵格繞流數(shù)值模擬研究[J].空氣動力學(xué)學(xué)報(bào),2009,27(1):78-82.

    [45] Wm David Washington,Pamela F Booth,Mark S Miller.Curvature and leading edge sweep back effects on grid fin aerodynamic characteristics[J].AIAA-93-3480-Cp,1993:659-667.

    [46] 譚獻(xiàn)忠,鄧帆,陳少松.翼面氣動外形對柵格翼減阻的影響[J].實(shí)驗(yàn)力學(xué),2013,28(2):255-260.

    [47] Daniel Guyot,Erich Schulein.Novel locally swept lattice wings for missile control at high speeds[C]//45th AIAA Aerospace Sciences Meeting and Exhibit.Reno,Nevada:2007.

    [48] Marco Debiasi,Zeng Yan,Chng Tat Loon.Swept-back grid fins for transonic drag reduction[C]//28th AIAA Applied Aerodynamics Conference.Chicago,Illinois:2010.

    [49] Jinsheng Cai.Numerical study on choked flow over grid-fin configurations[J].Journal of Spacecraft and Rockets,2009,46(5):949-956.

    [50] 黎漢華,石玉紅.柵格翼國內(nèi)外研究現(xiàn)狀及發(fā)展趨勢[J].導(dǎo)彈與航天運(yùn)載技術(shù),2008(6):27-30.

    [51] Jerome P Fanucci,Michael J King,James J Gorman.Method for producing lattice fin for missiles or other fluid-born bodies[P].US006928715B2.Woburn,MA,US:2005.08.16.

    [52] Stephen R Perillo,Donald J Atkins.Challenges and emerging trends in store separation engineering-an air force seek eagle office perspective[C]//47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition.Orlando,Florida:2009.

    [53] Paul V Tartabini,Michael G Gilbert,James R Beaty.Flight performance feasibility studies for the max launch abort system[C]//AIAA Atmospheric Flight Mechanics (AFM)Conference.Boston,MA:2013.

    (編輯:呂耀輝)

    Review of aerodynamic characteristics and application of grid fin

    PENG Ke1,HU Fan1,ZHANG Wei-hua1,ZHOU Zhang2

    (1.College of Aerospace Science and Engineering, National University of Defense Technology, Changsha 410073, China;2.Jiangnan Industries Group CO.,LTD,Changsha 410200,China)

    Researd methods on aerodynamic characteristics of grid fin,were summarized in fields of experimental research,engineering prediction,and numerical simulation.Major aerodynamic characteristics of grid fin,latest geometrical configurations research progresses about curved grid fin and global&local swept grid fin,as well as application trends at home and abroad were discussed.Critical concerns about grid fin were analyzed in view of the bottleneck problems of design and application.

    grid fin;aerodynamic characteristics;geometrical configurations;application trends

    2014-06-20;

    :2014-10-20。

    國家自然科學(xué)基金(51105368);國防科技大學(xué)優(yōu)秀研究生創(chuàng)新資助項(xiàng)目(S130105)。

    彭科(1989—),男,博士生,研究方向?yàn)轱w行器總體設(shè)計(jì)與氣動外形設(shè)計(jì)優(yōu)化。E-mail:pengke_pk@163.com

    V412

    A

    1006-2793(2015)04-0458-07

    10.7673/j.issn.1006-2793.2015.04.002

    猜你喜歡
    超聲速柵格氣動
    中寰氣動執(zhí)行機(jī)構(gòu)
    高超聲速出版工程
    高超聲速飛行器
    基于鄰域柵格篩選的點(diǎn)云邊緣點(diǎn)提取方法*
    基于NACA0030的波紋狀翼型氣動特性探索
    基于反饋線性化的RLV氣動控制一體化設(shè)計(jì)
    超聲速旅行
    不同剖面形狀的柵格壁對柵格翼氣動特性的影響
    高超聲速大博弈
    太空探索(2014年5期)2014-07-12 09:53:28
    基于CVT排布的非周期柵格密度加權(quán)陣設(shè)計(jì)
    99re6热这里在线精品视频| 中文乱码字字幕精品一区二区三区| 午夜精品一区二区三区免费看| av在线老鸭窝| 亚洲av成人精品一区久久| 超碰97精品在线观看| 麻豆精品久久久久久蜜桃| 啦啦啦啦在线视频资源| 成人亚洲精品av一区二区| 精品人妻熟女av久视频| 免费在线观看成人毛片| 亚洲综合精品二区| av女优亚洲男人天堂| 日韩强制内射视频| 嫩草影院精品99| 美女国产视频在线观看| 国产av不卡久久| 国产男女内射视频| 日韩成人av中文字幕在线观看| 精品午夜福利在线看| av一本久久久久| 亚洲av中文av极速乱| 成人国产麻豆网| 2018国产大陆天天弄谢| 一区二区av电影网| 亚洲精品色激情综合| 日韩一本色道免费dvd| 日韩精品有码人妻一区| 日韩不卡一区二区三区视频在线| 全区人妻精品视频| 天天躁夜夜躁狠狠久久av| 国产精品偷伦视频观看了| 全区人妻精品视频| 2021天堂中文幕一二区在线观| 亚洲精品aⅴ在线观看| 国产精品成人在线| 婷婷色综合大香蕉| 99久久中文字幕三级久久日本| 激情五月婷婷亚洲| 色播亚洲综合网| 国产 精品1| 亚洲国产成人一精品久久久| 真实男女啪啪啪动态图| 欧美高清成人免费视频www| 国产精品一区二区在线观看99| 亚洲精品自拍成人| 精品久久久久久电影网| 自拍欧美九色日韩亚洲蝌蚪91 | 成人美女网站在线观看视频| 伊人久久国产一区二区| 国产亚洲av片在线观看秒播厂| 亚洲成色77777| 欧美日韩综合久久久久久| 午夜日本视频在线| 黄片无遮挡物在线观看| 亚洲av男天堂| 欧美少妇被猛烈插入视频| 啦啦啦中文免费视频观看日本| 国产亚洲精品久久久com| 韩国高清视频一区二区三区| 新久久久久国产一级毛片| 国产成人一区二区在线| 黄色配什么色好看| 人妻夜夜爽99麻豆av| 欧美日韩精品成人综合77777| xxx大片免费视频| 在线观看美女被高潮喷水网站| 亚洲精品国产av成人精品| 婷婷色av中文字幕| 国产成人精品一,二区| 高清视频免费观看一区二区| 欧美日韩一区二区视频在线观看视频在线 | 久久久久久久久久成人| 人体艺术视频欧美日本| 2022亚洲国产成人精品| 亚洲一级一片aⅴ在线观看| 深爱激情五月婷婷| 好男人视频免费观看在线| 亚洲欧美一区二区三区国产| 蜜桃久久精品国产亚洲av| 免费看不卡的av| 人体艺术视频欧美日本| 亚洲av国产av综合av卡| 欧美日韩精品成人综合77777| 国产av不卡久久| 少妇猛男粗大的猛烈进出视频 | av福利片在线观看| 亚洲精品乱码久久久久久按摩| 九九久久精品国产亚洲av麻豆| 欧美少妇被猛烈插入视频| 国语对白做爰xxxⅹ性视频网站| 色视频www国产| 99热这里只有是精品50| 不卡视频在线观看欧美| 成人毛片a级毛片在线播放| 99热网站在线观看| 深爱激情五月婷婷| 日本欧美国产在线视频| 国产一级毛片在线| 黄色欧美视频在线观看| 中国美白少妇内射xxxbb| 永久网站在线| 又大又黄又爽视频免费| 91久久精品国产一区二区三区| 日韩欧美精品免费久久| 亚洲精品中文字幕在线视频 | 国产一区二区亚洲精品在线观看| 亚洲aⅴ乱码一区二区在线播放| 日韩一区二区三区影片| 国产色爽女视频免费观看| 亚洲国产日韩一区二区| 国产黄片美女视频| 肉色欧美久久久久久久蜜桃 | 熟妇人妻不卡中文字幕| 中国美白少妇内射xxxbb| 1000部很黄的大片| 五月伊人婷婷丁香| 激情 狠狠 欧美| 麻豆成人午夜福利视频| 中文字幕亚洲精品专区| 日韩视频在线欧美| 在线a可以看的网站| 亚洲成人中文字幕在线播放| 久久精品综合一区二区三区| 美女cb高潮喷水在线观看| 国产爽快片一区二区三区| 亚洲,欧美,日韩| 久久热精品热| 我的女老师完整版在线观看| 精品少妇久久久久久888优播| 亚洲精品第二区| 观看免费一级毛片| 国产欧美亚洲国产| 在线观看av片永久免费下载| 亚洲欧美成人综合另类久久久| 18禁在线播放成人免费| 好男人在线观看高清免费视频| 亚洲人成网站在线播| 人人妻人人看人人澡| 欧美xxⅹ黑人| 久久久久精品久久久久真实原创| 欧美另类一区| 一本一本综合久久| 精品亚洲乱码少妇综合久久| 国产欧美亚洲国产| 欧美日韩在线观看h| 中文精品一卡2卡3卡4更新| 国语对白做爰xxxⅹ性视频网站| 国产91av在线免费观看| 91狼人影院| 亚洲真实伦在线观看| 七月丁香在线播放| 少妇猛男粗大的猛烈进出视频 | 亚洲天堂av无毛| 久久人人爽av亚洲精品天堂 | 直男gayav资源| 少妇高潮的动态图| 国产精品成人在线| av在线天堂中文字幕| 简卡轻食公司| 精品久久久噜噜| 在线亚洲精品国产二区图片欧美 | 国产高潮美女av| 日韩大片免费观看网站| 大片免费播放器 马上看| 又爽又黄无遮挡网站| 久热这里只有精品99| 国产熟女欧美一区二区| 国产成人91sexporn| 91久久精品国产一区二区成人| 美女视频免费永久观看网站| 97超碰精品成人国产| 久久久亚洲精品成人影院| 大片电影免费在线观看免费| 91久久精品国产一区二区成人| 免费看av在线观看网站| 尾随美女入室| 国产黄色免费在线视频| 亚洲精品国产成人久久av| 精品一区二区三区视频在线| 涩涩av久久男人的天堂| 欧美最新免费一区二区三区| 男女下面进入的视频免费午夜| 国产精品女同一区二区软件| 久久精品国产亚洲网站| av专区在线播放| 网址你懂的国产日韩在线| 久久久久精品久久久久真实原创| 日本爱情动作片www.在线观看| 超碰97精品在线观看| 国产成人freesex在线| 久久国产乱子免费精品| 亚洲av一区综合| 啦啦啦啦在线视频资源| 麻豆成人av视频| 日韩人妻高清精品专区| 大片电影免费在线观看免费| 下体分泌物呈黄色| 黄色欧美视频在线观看| 91久久精品电影网| 国产av不卡久久| 下体分泌物呈黄色| 日韩欧美一区视频在线观看 | 丝袜喷水一区| 精品国产露脸久久av麻豆| www.色视频.com| 偷拍熟女少妇极品色| 三级国产精品片| 国产成人一区二区在线| 久久国产乱子免费精品| 日本av手机在线免费观看| 欧美极品一区二区三区四区| 人妻制服诱惑在线中文字幕| 3wmmmm亚洲av在线观看| 一区二区三区四区激情视频| 极品教师在线视频| 日本熟妇午夜| 少妇丰满av| 欧美成人a在线观看| 伊人久久国产一区二区| 人体艺术视频欧美日本| 亚洲自偷自拍三级| 亚洲av免费高清在线观看| 亚洲人成网站高清观看| 搞女人的毛片| 午夜日本视频在线| 亚洲av电影在线观看一区二区三区 | 久久久欧美国产精品| 亚洲精品久久久久久婷婷小说| 一级二级三级毛片免费看| 国产v大片淫在线免费观看| 国产精品.久久久| 国产精品久久久久久精品电影| 插阴视频在线观看视频| 成人美女网站在线观看视频| 国内少妇人妻偷人精品xxx网站| av一本久久久久| 美女视频免费永久观看网站| 网址你懂的国产日韩在线| 日韩大片免费观看网站| 亚洲国产精品成人久久小说| 联通29元200g的流量卡| 中文字幕av成人在线电影| 国产精品精品国产色婷婷| 99久久人妻综合| 肉色欧美久久久久久久蜜桃 | 人妻 亚洲 视频| 国产 一区精品| 精品久久久精品久久久| 成人亚洲精品av一区二区| 国产乱人偷精品视频| 赤兔流量卡办理| 国产白丝娇喘喷水9色精品| 一本色道久久久久久精品综合| 亚洲综合色惰| a级毛片免费高清观看在线播放| 一级黄片播放器| 尾随美女入室| 亚洲精品成人久久久久久| 少妇裸体淫交视频免费看高清| 麻豆成人午夜福利视频| 波多野结衣巨乳人妻| 国产 一区精品| 欧美日韩国产mv在线观看视频 | 日韩电影二区| 亚洲成人精品中文字幕电影| 久久99热这里只频精品6学生| 丝瓜视频免费看黄片| 日本-黄色视频高清免费观看| 精品一区在线观看国产| 亚洲成人一二三区av| 亚洲经典国产精华液单| 国产成年人精品一区二区| 99热全是精品| 国产黄片美女视频| 2018国产大陆天天弄谢| 热re99久久精品国产66热6| 久久女婷五月综合色啪小说 | 少妇人妻精品综合一区二区| 丝袜脚勾引网站| 在线亚洲精品国产二区图片欧美 | 国产精品久久久久久精品古装| 少妇人妻 视频| 联通29元200g的流量卡| 久久精品国产鲁丝片午夜精品| 国产欧美亚洲国产| 国产欧美亚洲国产| 黄片无遮挡物在线观看| 男男h啪啪无遮挡| 久久久久久久亚洲中文字幕| 91午夜精品亚洲一区二区三区| 超碰av人人做人人爽久久| 亚洲精品国产av蜜桃| 九九在线视频观看精品| 一级片'在线观看视频| 国产视频内射| 女人久久www免费人成看片| 搡女人真爽免费视频火全软件| 国产精品嫩草影院av在线观看| 亚洲人成网站在线观看播放| 黄片无遮挡物在线观看| 亚洲精品国产av成人精品| 直男gayav资源| 日韩av不卡免费在线播放| 听说在线观看完整版免费高清| 精品久久久久久久末码| 人妻少妇偷人精品九色| 亚洲无线观看免费| 国产精品.久久久| 少妇人妻久久综合中文| 精品人妻熟女av久视频| 99久久人妻综合| av在线蜜桃| 国产男人的电影天堂91| 日本wwww免费看| 中文乱码字字幕精品一区二区三区| 亚洲无线观看免费| 亚洲av.av天堂| 爱豆传媒免费全集在线观看| 黄色视频在线播放观看不卡| 亚洲无线观看免费| 国产亚洲5aaaaa淫片| 毛片一级片免费看久久久久| 韩国av在线不卡| 欧美潮喷喷水| 亚洲欧洲国产日韩| 亚洲最大成人av| videos熟女内射| 亚洲经典国产精华液单| 欧美最新免费一区二区三区| av一本久久久久| 亚洲久久久久久中文字幕| 熟女av电影| 又大又黄又爽视频免费| 99久久精品国产国产毛片| 精品亚洲乱码少妇综合久久| 在线看a的网站| 色视频在线一区二区三区| 亚洲国产成人一精品久久久| 午夜福利在线观看免费完整高清在| 免费黄频网站在线观看国产| 国产成人freesex在线| 国产男人的电影天堂91| 美女视频免费永久观看网站| 亚洲av日韩在线播放| 日韩人妻高清精品专区| 18禁在线无遮挡免费观看视频| 日日啪夜夜爽| 国产精品.久久久| 99热国产这里只有精品6| 亚洲无线观看免费| 国产永久视频网站| 久久人人爽人人片av| 成人亚洲精品av一区二区| 久久精品国产a三级三级三级| 网址你懂的国产日韩在线| 男女国产视频网站| 极品少妇高潮喷水抽搐| 久久99精品国语久久久| 成人无遮挡网站| 欧美高清成人免费视频www| 日韩制服骚丝袜av| 天堂俺去俺来也www色官网| 肉色欧美久久久久久久蜜桃 | 免费在线观看成人毛片| 久久久久久久久久人人人人人人| 国产精品久久久久久精品电影| 99久久精品国产国产毛片| 欧美激情国产日韩精品一区| 国产成人免费观看mmmm| 欧美 日韩 精品 国产| h日本视频在线播放| 在线观看美女被高潮喷水网站| 嘟嘟电影网在线观看| 男的添女的下面高潮视频| 王馨瑶露胸无遮挡在线观看| 高清在线视频一区二区三区| 国产精品国产三级专区第一集| 亚洲欧美清纯卡通| 精品酒店卫生间| 纵有疾风起免费观看全集完整版| 激情 狠狠 欧美| 丝袜喷水一区| 街头女战士在线观看网站| 亚洲最大成人手机在线| 国产成人午夜福利电影在线观看| 天天一区二区日本电影三级| 日本猛色少妇xxxxx猛交久久| 亚洲图色成人| 亚洲真实伦在线观看| 国产午夜福利久久久久久| 最近最新中文字幕免费大全7| 久久久久久久国产电影| 欧美精品一区二区大全| 国产探花极品一区二区| 国产男女内射视频| 在线亚洲精品国产二区图片欧美 | 国产成人freesex在线| 一边亲一边摸免费视频| 国产精品久久久久久久久免| 三级国产精品欧美在线观看| 国产男女内射视频| 中文在线观看免费www的网站| 日韩人妻高清精品专区| 五月伊人婷婷丁香| 午夜爱爱视频在线播放| 国产高清不卡午夜福利| 七月丁香在线播放| 亚洲av不卡在线观看| 亚洲av欧美aⅴ国产| 成人无遮挡网站| 午夜福利视频1000在线观看| 看十八女毛片水多多多| 成人特级av手机在线观看| av女优亚洲男人天堂| 国产午夜精品一二区理论片| 国产成人精品婷婷| 欧美xxxx性猛交bbbb| 成人高潮视频无遮挡免费网站| 亚洲天堂av无毛| 一级毛片电影观看| 国产高清三级在线| 一区二区av电影网| 国产精品熟女久久久久浪| 日韩电影二区| 亚洲自拍偷在线| 国产中年淑女户外野战色| 国产91av在线免费观看| av播播在线观看一区| 国产精品一区www在线观看| 亚洲第一区二区三区不卡| 高清av免费在线| 别揉我奶头 嗯啊视频| 汤姆久久久久久久影院中文字幕| 日韩 亚洲 欧美在线| 伊人久久国产一区二区| 日韩中字成人| 亚洲成人久久爱视频| 伦理电影大哥的女人| 亚洲精品国产成人久久av| 在线观看人妻少妇| 国产毛片在线视频| a级毛片免费高清观看在线播放| 亚洲在线观看片| 男人舔奶头视频| 亚洲国产欧美人成| 国产精品麻豆人妻色哟哟久久| 国产精品偷伦视频观看了| 王馨瑶露胸无遮挡在线观看| 久久久精品免费免费高清| 久久人人爽人人爽人人片va| 国产一区二区三区av在线| 亚洲精品成人av观看孕妇| 国产免费一级a男人的天堂| 人妻 亚洲 视频| 亚洲欧美日韩无卡精品| 午夜激情福利司机影院| 国产日韩欧美在线精品| 成人国产av品久久久| 干丝袜人妻中文字幕| 精品国产一区二区三区久久久樱花 | 极品少妇高潮喷水抽搐| 成人国产麻豆网| 中文字幕亚洲精品专区| 亚洲成人中文字幕在线播放| 一级毛片久久久久久久久女| 另类亚洲欧美激情| 一级片'在线观看视频| 观看免费一级毛片| 国产成人免费观看mmmm| 99热这里只有精品一区| 免费播放大片免费观看视频在线观看| 国产成人精品婷婷| 国产色婷婷99| 国产欧美另类精品又又久久亚洲欧美| 观看免费一级毛片| 三级经典国产精品| 少妇猛男粗大的猛烈进出视频 | 麻豆国产97在线/欧美| 在线观看一区二区三区激情| 国产一区亚洲一区在线观看| 久热久热在线精品观看| 免费av不卡在线播放| 欧美成人a在线观看| 欧美老熟妇乱子伦牲交| 欧美亚洲 丝袜 人妻 在线| 日韩亚洲欧美综合| 麻豆久久精品国产亚洲av| 国产精品国产三级国产av玫瑰| 老女人水多毛片| 亚洲精品久久久久久婷婷小说| 亚洲精品第二区| 最近中文字幕2019免费版| 色哟哟·www| 男女无遮挡免费网站观看| 国产av不卡久久| 日韩电影二区| 久久久午夜欧美精品| 久久久久性生活片| 亚洲精品亚洲一区二区| 国内精品宾馆在线| 身体一侧抽搐| av国产精品久久久久影院| 久久久久国产网址| 交换朋友夫妻互换小说| 最近2019中文字幕mv第一页| 精品国产露脸久久av麻豆| 春色校园在线视频观看| videos熟女内射| 黑人高潮一二区| 久久精品熟女亚洲av麻豆精品| 男人舔奶头视频| 天堂俺去俺来也www色官网| 久久久久久久久久久免费av| 王馨瑶露胸无遮挡在线观看| 午夜精品一区二区三区免费看| 高清av免费在线| 成人美女网站在线观看视频| 大码成人一级视频| 欧美日韩亚洲高清精品| 中文欧美无线码| 久久久久久久大尺度免费视频| av在线蜜桃| 狂野欧美白嫩少妇大欣赏| 国产黄色免费在线视频| 麻豆乱淫一区二区| 国产精品福利在线免费观看| 国产淫片久久久久久久久| 日韩人妻高清精品专区| 日韩免费高清中文字幕av| 精品国产露脸久久av麻豆| 深爱激情五月婷婷| 成人高潮视频无遮挡免费网站| 久久韩国三级中文字幕| 91狼人影院| 国产免费视频播放在线视频| 欧美+日韩+精品| 免费观看无遮挡的男女| 国产成人精品一,二区| 色哟哟·www| 97超碰精品成人国产| 久久久久久久精品精品| 三级国产精品片| 少妇人妻一区二区三区视频| 欧美极品一区二区三区四区| 免费av不卡在线播放| 禁无遮挡网站| 亚洲成人av在线免费| 一本一本综合久久| 天堂中文最新版在线下载 | 成年女人在线观看亚洲视频 | 欧美高清性xxxxhd video| 我的老师免费观看完整版| 韩国av在线不卡| 日韩欧美精品v在线| 午夜精品国产一区二区电影 | av在线播放精品| 国产毛片在线视频| 婷婷色综合大香蕉| 97在线人人人人妻| 日韩精品有码人妻一区| 天美传媒精品一区二区| 一个人看的www免费观看视频| 男插女下体视频免费在线播放| 久久久成人免费电影| 亚洲在线观看片| av一本久久久久| 免费av观看视频| av免费观看日本| 精品亚洲乱码少妇综合久久| 日韩三级伦理在线观看| 国产日韩欧美在线精品| 久久久精品94久久精品| 国产精品久久久久久精品电影小说 | 水蜜桃什么品种好| 制服丝袜香蕉在线| 韩国av在线不卡| 日日啪夜夜撸| 国产精品久久久久久精品电影小说 | 一个人看视频在线观看www免费| 国产成人福利小说| 国产精品99久久久久久久久| av卡一久久| 伦精品一区二区三区| 少妇猛男粗大的猛烈进出视频 | 亚洲婷婷狠狠爱综合网| av在线老鸭窝| 精品久久国产蜜桃| 国产探花极品一区二区| 一级毛片我不卡| 男女啪啪激烈高潮av片| 极品教师在线视频| 国产亚洲5aaaaa淫片| 亚洲精品国产av蜜桃| 五月玫瑰六月丁香| 国产伦理片在线播放av一区| 久久久成人免费电影| av专区在线播放| 在线观看国产h片| 成年版毛片免费区| 在线 av 中文字幕| 欧美日韩一区二区视频在线观看视频在线 | 下体分泌物呈黄色| 久久国产乱子免费精品| 亚洲aⅴ乱码一区二区在线播放| 天堂中文最新版在线下载 | 久久精品国产a三级三级三级| 国产av不卡久久| 精品一区二区三卡| 天美传媒精品一区二区| 色网站视频免费| 久久韩国三级中文字幕| 国产女主播在线喷水免费视频网站| 午夜爱爱视频在线播放| 免费观看a级毛片全部|