• <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ì)
    欧美 日韩 精品 国产| 高清在线国产一区| 在线观看免费视频日本深夜| 亚洲国产中文字幕在线视频| 交换朋友夫妻互换小说| 精品人妻在线不人妻| 一进一出抽搐gif免费好疼 | 18禁裸乳无遮挡动漫免费视频| 99精品欧美一区二区三区四区| 国产成+人综合+亚洲专区| 久久精品91无色码中文字幕| 男女床上黄色一级片免费看| 国产一区有黄有色的免费视频| 侵犯人妻中文字幕一二三四区| 高清欧美精品videossex| 欧洲精品卡2卡3卡4卡5卡区| 国产精品 国内视频| 亚洲人成77777在线视频| 男女午夜视频在线观看| 999久久久精品免费观看国产| 精品久久久精品久久久| 欧美精品av麻豆av| 女人精品久久久久毛片| 亚洲国产欧美一区二区综合| 日本精品一区二区三区蜜桃| 嫩草影视91久久| 香蕉久久夜色| 国产成人av激情在线播放| 人人妻人人澡人人看| 欧美黄色片欧美黄色片| 久久天堂一区二区三区四区| 大型黄色视频在线免费观看| 欧美人与性动交α欧美软件| 日韩精品免费视频一区二区三区| 成年人免费黄色播放视频| 久久草成人影院| 久久中文看片网| 精品一区二区三区视频在线观看免费 | 日韩欧美免费精品| 美女高潮喷水抽搐中文字幕| 午夜福利视频在线观看免费| 老司机在亚洲福利影院| 色综合婷婷激情| 成人三级做爰电影| 黄色片一级片一级黄色片| 51午夜福利影视在线观看| 美女高潮喷水抽搐中文字幕| 日韩精品免费视频一区二区三区| 飞空精品影院首页| 一本大道久久a久久精品| 老司机福利观看| 一级作爱视频免费观看| 亚洲aⅴ乱码一区二区在线播放 | 天天躁日日躁夜夜躁夜夜| 99热网站在线观看| 亚洲精华国产精华精| 国内毛片毛片毛片毛片毛片| 亚洲av日韩精品久久久久久密| 午夜福利免费观看在线| 国产一区有黄有色的免费视频| 久久人人97超碰香蕉20202| 黄网站色视频无遮挡免费观看| 一级毛片精品| 99精品在免费线老司机午夜| 夜夜躁狠狠躁天天躁| 久久国产精品人妻蜜桃| 免费久久久久久久精品成人欧美视频| 大码成人一级视频| 精品午夜福利视频在线观看一区| 亚洲第一欧美日韩一区二区三区| 国产色视频综合| 老鸭窝网址在线观看| 狠狠婷婷综合久久久久久88av| 老司机影院毛片| 女人久久www免费人成看片| 老司机影院毛片| 极品教师在线免费播放| 久久人妻福利社区极品人妻图片| 欧美激情极品国产一区二区三区| 捣出白浆h1v1| 人人妻人人澡人人爽人人夜夜| 90打野战视频偷拍视频| 美女高潮喷水抽搐中文字幕| 热99久久久久精品小说推荐| 在线视频色国产色| 在线观看www视频免费| 美国免费a级毛片| 国产aⅴ精品一区二区三区波| 亚洲欧美激情综合另类| 午夜老司机福利片| 男女之事视频高清在线观看| 不卡一级毛片| 国产成+人综合+亚洲专区| 亚洲欧洲精品一区二区精品久久久| 亚洲精品乱久久久久久| 亚洲第一欧美日韩一区二区三区| 亚洲国产中文字幕在线视频| 国产亚洲欧美98| 丰满人妻熟妇乱又伦精品不卡| 视频区图区小说| 91国产中文字幕| 国产精品偷伦视频观看了| 日韩人妻精品一区2区三区| 亚洲九九香蕉| 老汉色av国产亚洲站长工具| 精品国产美女av久久久久小说| 久久国产精品男人的天堂亚洲| 日韩制服丝袜自拍偷拍| 每晚都被弄得嗷嗷叫到高潮| 女人被躁到高潮嗷嗷叫费观| 日韩制服丝袜自拍偷拍| 涩涩av久久男人的天堂| videos熟女内射| 欧美av亚洲av综合av国产av| tube8黄色片| 黄色 视频免费看| 婷婷精品国产亚洲av在线 | 最近最新中文字幕大全免费视频| 老鸭窝网址在线观看| cao死你这个sao货| 国产精华一区二区三区| 黄频高清免费视频| 91成年电影在线观看| 国产免费现黄频在线看| 亚洲,欧美精品.| 色在线成人网| 18禁观看日本| 亚洲少妇的诱惑av| 免费黄频网站在线观看国产| 久久人人97超碰香蕉20202| 亚洲国产看品久久| 日本欧美视频一区| 啦啦啦免费观看视频1| av在线播放免费不卡| 亚洲片人在线观看| 亚洲片人在线观看| 亚洲精品久久成人aⅴ小说| 久久狼人影院| 一区福利在线观看| 一区二区日韩欧美中文字幕| 91麻豆精品激情在线观看国产 | 欧美不卡视频在线免费观看 | 一级,二级,三级黄色视频| 久久精品亚洲av国产电影网| 看免费av毛片| 欧美最黄视频在线播放免费 | 精品午夜福利视频在线观看一区| 欧美国产精品一级二级三级| 亚洲 欧美一区二区三区| 咕卡用的链子| www日本在线高清视频| 天堂俺去俺来也www色官网| 国产精品久久久av美女十八| 国产99白浆流出| 日本五十路高清| 午夜免费观看网址| 国产亚洲欧美在线一区二区| 午夜亚洲福利在线播放| 69精品国产乱码久久久| x7x7x7水蜜桃| 十分钟在线观看高清视频www| 午夜久久久在线观看| 精品第一国产精品| 9191精品国产免费久久| 好看av亚洲va欧美ⅴa在| 欧美日韩亚洲高清精品| 国产亚洲av高清不卡| 亚洲 国产 在线| 好男人电影高清在线观看| 久久人人97超碰香蕉20202| 极品少妇高潮喷水抽搐| 日韩 欧美 亚洲 中文字幕| 他把我摸到了高潮在线观看| 老汉色av国产亚洲站长工具| 99国产精品一区二区蜜桃av | 99久久综合精品五月天人人| 欧美日韩亚洲国产一区二区在线观看 | 天天操日日干夜夜撸| 中国美女看黄片| 青草久久国产| 久99久视频精品免费| 成人三级做爰电影| 成熟少妇高潮喷水视频| xxx96com| 亚洲中文日韩欧美视频| 美女国产高潮福利片在线看| 亚洲欧美激情在线| 国产精品影院久久| 日韩欧美一区二区三区在线观看 | 久久亚洲精品不卡| 久久精品熟女亚洲av麻豆精品| 国产精品免费一区二区三区在线 | 国产亚洲精品久久久久5区| 交换朋友夫妻互换小说| 国产精品成人在线| 啪啪无遮挡十八禁网站| 99riav亚洲国产免费| 精品无人区乱码1区二区| 成熟少妇高潮喷水视频| 精品国产乱子伦一区二区三区| 亚洲一区二区三区不卡视频| 亚洲午夜精品一区,二区,三区| 女同久久另类99精品国产91| 老熟女久久久| 亚洲片人在线观看| 亚洲精品美女久久av网站| 国产欧美亚洲国产| 欧美日韩乱码在线| 免费人成视频x8x8入口观看| 国产色视频综合| 成熟少妇高潮喷水视频| 黑人操中国人逼视频| 久久香蕉精品热| 久久国产乱子伦精品免费另类| 18禁国产床啪视频网站| 国产av精品麻豆| 一边摸一边抽搐一进一小说 | 超碰成人久久| 亚洲 国产 在线| 在线天堂中文资源库| 亚洲国产欧美一区二区综合| 操出白浆在线播放| 亚洲综合色网址| 久久久久国产一级毛片高清牌| 久久久久国产精品人妻aⅴ院 | 久久久精品免费免费高清| 欧美日韩中文字幕国产精品一区二区三区 | 老司机影院毛片| 午夜精品在线福利| 亚洲一区高清亚洲精品| 免费在线观看日本一区| 九色亚洲精品在线播放| av线在线观看网站| 久久中文看片网| 男女免费视频国产| 18禁黄网站禁片午夜丰满| 美女扒开内裤让男人捅视频| 国产激情欧美一区二区| 国产精品成人在线| 国产成人欧美在线观看 | 久热这里只有精品99| 午夜福利免费观看在线| 高清毛片免费观看视频网站 | 91av网站免费观看| 后天国语完整版免费观看| 美女扒开内裤让男人捅视频| 亚洲中文日韩欧美视频| 亚洲欧美激情在线| 一级a爱片免费观看的视频| 国产欧美日韩精品亚洲av| 国产成人精品无人区| 国产国语露脸激情在线看| 高清av免费在线| 精品久久久精品久久久| 身体一侧抽搐| 人妻 亚洲 视频| 亚洲国产精品一区二区三区在线| 99久久国产精品久久久| 精品久久久久久久久久免费视频 | 国产97色在线日韩免费| 一级毛片精品| 999久久久精品免费观看国产| av在线播放免费不卡| 国产av精品麻豆| 国产成人一区二区三区免费视频网站| 9191精品国产免费久久| 伊人久久大香线蕉亚洲五| 亚洲,欧美精品.| 日韩人妻精品一区2区三区| 精品无人区乱码1区二区| 日日夜夜操网爽| 国产成人欧美| 免费在线观看完整版高清| 18禁裸乳无遮挡免费网站照片 | 久久久久久久午夜电影 | 欧美激情高清一区二区三区| av福利片在线| av国产精品久久久久影院| 亚洲中文字幕日韩| 在线视频色国产色| 90打野战视频偷拍视频| 91国产中文字幕| 法律面前人人平等表现在哪些方面| 王馨瑶露胸无遮挡在线观看| 国产视频一区二区在线看| 精品亚洲成a人片在线观看| 中文字幕人妻丝袜一区二区| 青草久久国产| 人人妻人人添人人爽欧美一区卜| 国产亚洲av高清不卡| 午夜免费成人在线视频| 免费女性裸体啪啪无遮挡网站| av欧美777| 国产亚洲精品久久久久5区| 国产成人免费观看mmmm| 久久精品国产综合久久久| 欧美 日韩 精品 国产| a级片在线免费高清观看视频| 高清黄色对白视频在线免费看| 久久草成人影院| 精品高清国产在线一区| 亚洲精品国产区一区二| 国产不卡av网站在线观看| 男人的好看免费观看在线视频 | 大香蕉久久网| 99国产精品免费福利视频| 最近最新免费中文字幕在线| 成熟少妇高潮喷水视频| www日本在线高清视频| 午夜福利视频在线观看免费| 老司机影院毛片| 国产国语露脸激情在线看| 午夜福利影视在线免费观看| 亚洲国产毛片av蜜桃av| а√天堂www在线а√下载 | 国产视频一区二区在线看| 亚洲视频免费观看视频| 啪啪无遮挡十八禁网站| 久久香蕉激情| 久久人人爽av亚洲精品天堂| 水蜜桃什么品种好| 天天影视国产精品| 国产精品自产拍在线观看55亚洲 | 在线观看免费视频日本深夜| 在线观看免费视频网站a站| 午夜亚洲福利在线播放| 亚洲成a人片在线一区二区| 看片在线看免费视频| 老司机靠b影院| 欧美激情 高清一区二区三区| 啪啪无遮挡十八禁网站| 91麻豆av在线| 亚洲午夜理论影院| 亚洲av成人不卡在线观看播放网| 成年人免费黄色播放视频| 成人手机av| 丰满迷人的少妇在线观看| 国产精品1区2区在线观看. | 亚洲全国av大片| 天堂动漫精品| 免费观看a级毛片全部| xxxhd国产人妻xxx| 激情视频va一区二区三区| 19禁男女啪啪无遮挡网站| 男人操女人黄网站| 人妻丰满熟妇av一区二区三区 | 天堂动漫精品| 一进一出抽搐动态| 亚洲专区字幕在线| 亚洲欧美色中文字幕在线| 国产亚洲精品久久久久5区| 可以免费在线观看a视频的电影网站| 国产精品乱码一区二三区的特点 | 超色免费av| 成人18禁高潮啪啪吃奶动态图| 久久国产精品影院| 亚洲视频免费观看视频| 国产极品粉嫩免费观看在线| 精品久久久精品久久久| 精品国产一区二区三区四区第35| 免费人成视频x8x8入口观看| 成人18禁在线播放| 丰满饥渴人妻一区二区三| 久久香蕉国产精品| 免费观看a级毛片全部| 黄色怎么调成土黄色| 国产欧美日韩一区二区三区在线| 黄色视频不卡| 男人操女人黄网站| av国产精品久久久久影院| 12—13女人毛片做爰片一| 亚洲中文字幕日韩| 女人高潮潮喷娇喘18禁视频| 欧美激情极品国产一区二区三区| 性少妇av在线| 热99国产精品久久久久久7| 久热这里只有精品99| 国产精华一区二区三区| 国产精品一区二区在线不卡| 国产精品亚洲一级av第二区| 亚洲熟妇中文字幕五十中出 | 最近最新中文字幕大全免费视频| 国产精品 欧美亚洲| 欧美日韩黄片免| 成人亚洲精品一区在线观看| 色综合欧美亚洲国产小说| 一区二区三区激情视频| 男女床上黄色一级片免费看| 一进一出抽搐动态| 成年人午夜在线观看视频| 亚洲第一av免费看| 欧美亚洲 丝袜 人妻 在线| xxxhd国产人妻xxx| 搡老岳熟女国产| 99香蕉大伊视频| 久久久久国产精品人妻aⅴ院 | 一区二区三区精品91| 亚洲九九香蕉| av国产精品久久久久影院| av欧美777| 美女 人体艺术 gogo| 欧美激情久久久久久爽电影 | 免费观看精品视频网站| 村上凉子中文字幕在线| cao死你这个sao货| 中文字幕高清在线视频| 一区二区日韩欧美中文字幕| 18禁黄网站禁片午夜丰满| 九色亚洲精品在线播放| 国产成人精品在线电影| 中文字幕最新亚洲高清| 极品少妇高潮喷水抽搐| 男人的好看免费观看在线视频 | 精品国产一区二区久久| 国产乱人伦免费视频| 亚洲精品美女久久av网站| 久久久国产一区二区| 热re99久久精品国产66热6| 99精国产麻豆久久婷婷| 成人影院久久| 国产亚洲精品第一综合不卡| 欧美乱妇无乱码| 亚洲一区二区三区欧美精品| av一本久久久久| 久久久久久久国产电影| 中文字幕人妻丝袜制服| 国产成人av激情在线播放| 无遮挡黄片免费观看| 波多野结衣一区麻豆| 99国产精品一区二区三区| 天堂动漫精品| 777米奇影视久久| 一进一出抽搐gif免费好疼 | 99国产综合亚洲精品| 淫妇啪啪啪对白视频| 自线自在国产av| 最近最新免费中文字幕在线| 丝袜美足系列| 在线永久观看黄色视频| 精品久久久久久电影网| 丰满饥渴人妻一区二区三| 日韩欧美免费精品| 成人18禁高潮啪啪吃奶动态图| 极品少妇高潮喷水抽搐| 亚洲黑人精品在线| 黄色毛片三级朝国网站| 国产av一区二区精品久久| 视频区图区小说| 女同久久另类99精品国产91| 久久中文字幕人妻熟女| 黑人欧美特级aaaaaa片| 亚洲国产毛片av蜜桃av| 一级作爱视频免费观看| 老司机在亚洲福利影院| 国产午夜精品久久久久久| 国产激情久久老熟女| 久久久久久久国产电影| 满18在线观看网站| 窝窝影院91人妻| 国产精品国产av在线观看| 亚洲欧美激情综合另类| 在线观看www视频免费| 狠狠婷婷综合久久久久久88av| 丰满迷人的少妇在线观看| 亚洲成人免费av在线播放| 在线观看免费日韩欧美大片| av不卡在线播放| 国产精品.久久久| 很黄的视频免费| 黄网站色视频无遮挡免费观看| 午夜亚洲福利在线播放| 啦啦啦视频在线资源免费观看| 国产蜜桃级精品一区二区三区 | 久久香蕉精品热| 一夜夜www| 日日摸夜夜添夜夜添小说| 久久精品国产亚洲av香蕉五月 | 成人国产一区最新在线观看| 99riav亚洲国产免费| av不卡在线播放| 精品久久久久久久毛片微露脸| 免费观看精品视频网站| 在线永久观看黄色视频| 国产成人欧美在线观看 | 两性夫妻黄色片| av免费在线观看网站| 午夜亚洲福利在线播放| 一级毛片高清免费大全| 免费在线观看完整版高清| 欧美中文综合在线视频| av天堂在线播放| 一级作爱视频免费观看| 又紧又爽又黄一区二区| 男女免费视频国产| 美女福利国产在线| 欧美黑人欧美精品刺激| √禁漫天堂资源中文www| 国产在线观看jvid| aaaaa片日本免费| 午夜视频精品福利| 久久精品成人免费网站| 色婷婷久久久亚洲欧美| 成在线人永久免费视频| 免费在线观看日本一区| 宅男免费午夜| 久久久国产一区二区| 国产男靠女视频免费网站| 18禁黄网站禁片午夜丰满| 999久久久精品免费观看国产| 亚洲五月天丁香| 欧美国产精品va在线观看不卡| 国产av精品麻豆| 91精品国产国语对白视频| 天堂动漫精品| 9色porny在线观看| 精品国产美女av久久久久小说| 熟女少妇亚洲综合色aaa.| 久久国产乱子伦精品免费另类| 99热网站在线观看| 亚洲 国产 在线| 操出白浆在线播放| 黄色 视频免费看| 久久九九热精品免费| 国产精品亚洲一级av第二区| 久久亚洲真实| 亚洲午夜精品一区,二区,三区| 性色av乱码一区二区三区2| 无限看片的www在线观看| 国产精品一区二区在线观看99| 日日爽夜夜爽网站| 捣出白浆h1v1| 丁香六月欧美| 天堂√8在线中文| 99久久综合精品五月天人人| 极品少妇高潮喷水抽搐| 在线观看免费午夜福利视频| 天堂俺去俺来也www色官网| 亚洲精华国产精华精| 欧美激情久久久久久爽电影 | 国产精品电影一区二区三区 | 国产精品欧美亚洲77777| 天天躁狠狠躁夜夜躁狠狠躁| aaaaa片日本免费| www.自偷自拍.com| 99热国产这里只有精品6| 热re99久久精品国产66热6| 视频在线观看一区二区三区| 757午夜福利合集在线观看| 大陆偷拍与自拍| 真人做人爱边吃奶动态| 捣出白浆h1v1| 国产高清国产精品国产三级| 亚洲片人在线观看| 欧美成人免费av一区二区三区 | 99在线人妻在线中文字幕 | 国产精品久久久人人做人人爽| 欧美日韩亚洲国产一区二区在线观看 | 色婷婷久久久亚洲欧美| 中文字幕人妻丝袜制服| 高清欧美精品videossex| 国产亚洲精品久久久久5区| 午夜福利影视在线免费观看| 精品久久久久久久毛片微露脸| 亚洲中文日韩欧美视频| 在线天堂中文资源库| 国产99久久九九免费精品| 久久久水蜜桃国产精品网| 黄色怎么调成土黄色| 日韩欧美三级三区| 淫妇啪啪啪对白视频| 日韩欧美三级三区| 精品久久久久久久久久免费视频 | 国产精品亚洲一级av第二区| 母亲3免费完整高清在线观看| 亚洲国产精品合色在线| а√天堂www在线а√下载 | 亚洲av美国av| 人人妻人人添人人爽欧美一区卜| 他把我摸到了高潮在线观看| 在线播放国产精品三级| 他把我摸到了高潮在线观看| 中文字幕最新亚洲高清| 人妻一区二区av| 人人澡人人妻人| 日本黄色日本黄色录像| 丁香六月欧美| 露出奶头的视频| 搡老熟女国产l中国老女人| 中文字幕另类日韩欧美亚洲嫩草| 国产亚洲精品一区二区www | 久久国产亚洲av麻豆专区| 久久精品人人爽人人爽视色| 王馨瑶露胸无遮挡在线观看| 大陆偷拍与自拍| av不卡在线播放| 交换朋友夫妻互换小说| 国产精品1区2区在线观看. | 国产精品永久免费网站| 国产精品免费一区二区三区在线 | 亚洲国产欧美网| 欧美黄色片欧美黄色片| 久久ye,这里只有精品| 精品视频人人做人人爽| 亚洲欧美激情综合另类| 精品电影一区二区在线| 天堂√8在线中文| 美女福利国产在线| 极品教师在线免费播放| 99国产精品一区二区三区| 电影成人av| 亚洲精品在线美女| 99精品久久久久人妻精品| 91字幕亚洲| 首页视频小说图片口味搜索| avwww免费|