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

    類(lèi)X-43A 高超聲速飛行器機(jī)體/推進(jìn)一體化氣動(dòng)設(shè)計(jì)分析和地面試驗(yàn)問(wèn)題評(píng)述

    2021-09-01 06:46:42秦開(kāi)宇
    載人航天 2021年4期
    關(guān)鍵詞:進(jìn)氣道風(fēng)洞試驗(yàn)超聲速

    梁 捷, 秦開(kāi)宇, 陳 力

    (1.電子科技大學(xué)航空航天學(xué)院, 成都 611731; 2.中國(guó)空氣動(dòng)力研究與發(fā)展中心, 綿陽(yáng) 621000;3.福州大學(xué)機(jī)械工程及自動(dòng)化學(xué)院, 福州 350108)

    1 引言

    20 世紀(jì)中期以來(lái),航空航天技術(shù)快速發(fā)展并日趨成熟,提升進(jìn)入空間與開(kāi)發(fā)利用空間的能力已成為各航天大國(guó)的國(guó)家戰(zhàn)略[1-3]。 與之契合,高超聲速飛行器作為航空航天領(lǐng)域的標(biāo)志性創(chuàng)新成果受到各航天強(qiáng)國(guó)的廣泛重視,特別在2004年,以超燃沖壓發(fā)動(dòng)機(jī)(Supersonic Combustion Ramjet,Scramjet)為動(dòng)力的X-43A 飛行試驗(yàn)取得成功,吸氣式高超聲速飛行器技術(shù)受到國(guó)內(nèi)外學(xué)者的極大關(guān)注[4-6]。

    吸氣式推進(jìn)對(duì)于在大氣層中持續(xù)飛行的巡航飛行器是不可或缺的重要組成部分,Scramjet 發(fā)動(dòng)機(jī)是吸氣推進(jìn)高超聲速飛行器的支撐技術(shù)。 近幾十年來(lái)雖然對(duì)Scramjet 發(fā)動(dòng)機(jī)開(kāi)展了廣泛研究并取得重大進(jìn)展,但在設(shè)計(jì)方法、實(shí)驗(yàn)數(shù)據(jù)庫(kù)、試驗(yàn)設(shè)備及發(fā)動(dòng)機(jī)性能驗(yàn)證等方面都存在局限性。

    美國(guó)在成功研制洲際彈道導(dǎo)彈、載人飛船和航天飛機(jī)這些以火箭為動(dòng)力的高超聲速飛行器之后,于20 世紀(jì)80 年代中期實(shí)施了以吸氣式超燃沖壓發(fā)動(dòng)機(jī)為動(dòng)力的國(guó)家空天飛機(jī)計(jì)劃(National Aero-Space Plane Program, NASP),推進(jìn)/機(jī)體一體 化(Propulsion/Airframe Integration, PAI) 是NASP 計(jì)劃的興趣所在,并提出了如圖1 所示的高超聲速飛行器PAI 系統(tǒng)一體化設(shè)計(jì)概念[7-9]。 吸氣式高超聲速飛行器設(shè)計(jì)很大程度上受到高效PAI 需求驅(qū)動(dòng),帶動(dòng)力條件下高超聲速飛行器氣動(dòng)分析和試驗(yàn)更是PAI 的一個(gè)重要組成部分。 但“一體化”帶來(lái)的機(jī)體和推進(jìn)系統(tǒng)強(qiáng)耦合效應(yīng),給飛行器設(shè)計(jì)和驗(yàn)證提出了極高的要求。 NASP 計(jì)劃試圖建造一種單級(jí)入軌(SSTO)試驗(yàn)飛機(jī)X-30,用于飛行試驗(yàn),驗(yàn)證雙模態(tài)Ramjet-Scramjet 發(fā)動(dòng)機(jī)技術(shù)[2]。 由于設(shè)計(jì)要求超越了當(dāng)時(shí)科技發(fā)展水平,該計(jì)劃雖幾經(jīng)調(diào)整,最終于1995 年1 月終止。

    圖1 高超聲速飛行器機(jī)體和推進(jìn)系統(tǒng)“一體化”概念[10]Fig.1 PAI concept of hypersonic vehicle[10]

    NASP 計(jì)劃下馬后,美國(guó)于20 世紀(jì)末制定了新的高超聲速技術(shù)發(fā)展戰(zhàn)略,提出了“從無(wú)人到有人、從一次性使用到重復(fù)使用”循序漸進(jìn)的指導(dǎo)思想和“以高超聲速巡航導(dǎo)彈為突破口,再發(fā)展軍用飛機(jī)和天地往返系統(tǒng)”的三步走發(fā)展戰(zhàn)略。 1996 年,在NASP 基礎(chǔ)上,NASA 啟動(dòng)了研究?jī)?nèi)容更加廣泛的高超聲速試驗(yàn)計(jì)劃,即Hyper-X。Hyper-X 計(jì)劃先從設(shè)計(jì)建造小型低成本驗(yàn)證飛行器入手開(kāi)展飛行試驗(yàn),以便驗(yàn)證高超聲速飛行器相關(guān)技術(shù)和設(shè)計(jì)方法論。 NASA 指出Hyper-X 計(jì)劃的研究核心是X-43 系列驗(yàn)證機(jī),該系列驗(yàn)證機(jī)布局實(shí)質(zhì)上是根據(jù)NASP 研究的概念飛行器縮尺獲得的。 這為Hyper-X 計(jì)劃提供了很多便利,可使用NASP 數(shù)據(jù)庫(kù)、經(jīng)驗(yàn)及改進(jìn)型的任務(wù)研究計(jì)劃[11]。 因此,Hyper-X 計(jì)劃發(fā)展過(guò)程可視為NASP 發(fā)展計(jì)劃的“逆向執(zhí)行”。

    2004 年11 月,以Scramjet 發(fā)動(dòng)機(jī)為推進(jìn)的X-43A 高超聲速飛行器創(chuàng)造了10Ma速度持續(xù)飛行10 s 的飛行紀(jì)錄。 X-43A 飛行器外部、內(nèi)部輪廓如圖2 所示。 該飛行器大約12 ft(3.66 m)長(zhǎng),5 ft(1.52 m)寬,2 ft(0.61 m)高,重約3000 lb(13345 N)[12]。 X-43A 驗(yàn)證了Scramjet 發(fā)動(dòng)機(jī)在既定試驗(yàn)條件下的氣動(dòng)力穩(wěn)定性、控制性能、飛行性能及飛行器機(jī)體結(jié)構(gòu)一體化和系統(tǒng)設(shè)計(jì);X-43A 飛行數(shù)據(jù)驗(yàn)證了基于風(fēng)洞試驗(yàn)和CFD 增量預(yù)測(cè)的氣動(dòng)特性數(shù)據(jù)。

    圖2 X-43A 飛行器外形[13]Fig.2 The X-43A aircraft configuration[13]

    2 吸氣式高超聲速飛行器氣動(dòng)分析

    Hyper-X 研究表明,要使吸氣式高超聲速飛行器發(fā)展成為有實(shí)用意義的飛行器,還需解決許多技術(shù)難題,包括帶Scramjet 動(dòng)力的推進(jìn)系統(tǒng)和PAI 設(shè)計(jì)、性能預(yù)測(cè)、試驗(yàn)驗(yàn)證和計(jì)算方法等。

    2.1 Scramjet 發(fā)動(dòng)機(jī)

    在一體化布局下,升力體/乘波體整個(gè)下腹部都是推進(jìn)系統(tǒng)的一部分,前體下表面作為發(fā)動(dòng)機(jī)入流的預(yù)壓縮面,為進(jìn)氣道入流提供經(jīng)過(guò)預(yù)壓縮的高質(zhì)量氣流,后體下表面為噴管出口,供更大的氣流膨脹面積(圖1)。 在設(shè)計(jì)馬赫數(shù)下飛行時(shí),要求前體發(fā)出的弓形激波搭在進(jìn)氣道唇緣上,不產(chǎn)生進(jìn)氣道溢流損失,具有最大推力[14]。

    長(zhǎng)期以來(lái),航天運(yùn)載器使用火箭發(fā)動(dòng)機(jī)作為動(dòng)力,積累了非常成熟的經(jīng)驗(yàn)[15]。 然而,現(xiàn)有的火箭動(dòng)力系統(tǒng),在運(yùn)行過(guò)程中必須自帶氧化劑,致使起飛總重大,有效載荷低。 研究表明,火箭比沖Isp已經(jīng)接近極限,進(jìn)一步提高的可能性十分有限[16-17]。與火箭動(dòng)力相比,吸氣推進(jìn)系統(tǒng)無(wú)需攜帶氧化劑,大大降低了起飛重量。 在同樣有效載荷下,二者起飛重量之比約為Wrocket/Wairplane=4[18]。 可見(jiàn),在同樣有效載荷下,吸氣式高超聲速飛行器的起飛總重大大降低了,而且吸氣式系統(tǒng)提供了更大安全性和靈活性并降低了飛行成本。

    2.2 升阻比

    升阻比是飛行器氣動(dòng)布局設(shè)計(jì)的重要參數(shù)。傳統(tǒng)的飛行器薄機(jī)翼布局設(shè)計(jì),在執(zhí)行高超聲速飛行條件下,升阻比極低,不能滿足設(shè)計(jì)要求。 目前許多研究者都把注意力轉(zhuǎn)移到PAI 布局的升力體和乘波體布局設(shè)計(jì)上[19-21]。

    圖3 中,實(shí)線表示一般飛行器在高超聲速條件下,難以打破的L/D 阻礙,表達(dá)式為式(1)[22]:

    式中,L為升力,D為阻力,Ma為馬赫數(shù);從上式可以看出,當(dāng)Ma無(wú)窮大時(shí), (L/D)max只有4,這就是飛行器難以跨越的高超聲速升阻比屏障(圖3 中,空心符號(hào)表示一般飛行器的升阻比)。 而采用“乘波體”布局設(shè)計(jì),可使飛行器在高超聲速下升阻比值達(dá)到5~8,其(L/D)max表示為式(2)[22]:

    圖3 中,虛線表示采用乘波體飛行器的最大升阻比線,圓形黑實(shí)心符號(hào)和方形黑實(shí)心符號(hào)分別表示半指數(shù)流乘波體和錐形流乘波體在高超聲速飛行時(shí)的升阻比[22]。 文獻(xiàn)研究結(jié)果表明,即使在非設(shè)計(jì)條件下,乘波體飛行器仍然保持很高的升阻比值[22]。 而且前緣鈍度對(duì)升阻比影響很??;乘波體布局對(duì)進(jìn)氣道提供算法均勻的流場(chǎng),并不產(chǎn)生橫向流動(dòng),這一點(diǎn)對(duì)發(fā)動(dòng)機(jī)有效燃燒非常重要。

    圖3 傳統(tǒng)飛行器和兩種乘波體飛行器的升阻比[22]Fig.3 Lift-to-drag ratio of conventional aircraft and two types of waverider[22]

    高超聲速概念飛行器分為巡航飛行器和加速飛行器。 “全球到達(dá)”的巡航飛行器設(shè)計(jì)目標(biāo)是達(dá)到最大航程。 根據(jù)式(3)所示Breguet 航程公式[14]:

    式中,R為航程;V0為來(lái)流速度;Isp為比沖;WI為巡航起始重量;WF為巡航結(jié)束重量。 從式(3)可以看出,在高超聲速飛行器概念方案設(shè)計(jì)階段,為了達(dá)到最大航程,必須對(duì)(L/D)×Isp實(shí)行優(yōu)化,單獨(dú)采用L/D或Isp來(lái)評(píng)估飛行器氣動(dòng)外形或發(fā)動(dòng)機(jī)性能是不合理的,因?yàn)轱w行器氣動(dòng)界面與推進(jìn)系統(tǒng)劃分方式不同會(huì)導(dǎo)致升阻比不同。乘波體具有較大升阻比,是全球到達(dá)的巡航飛行器理想的布局形式。

    2.3 吸氣式高超聲速飛行器PAI 設(shè)計(jì)分析

    基于PAI 概念設(shè)計(jì)的X-43A 吸氣式高超聲速飛行器(如圖4 所示),在氣動(dòng)布局上的一個(gè)顯著特點(diǎn)是機(jī)體和推進(jìn)系統(tǒng)高度綜合成一體,產(chǎn)生升力的部件、推進(jìn)系統(tǒng)和機(jī)身的功能不再相互獨(dú)立,而是以高度集成的升力體出現(xiàn),這是目前高超聲速飛行器的一種優(yōu)選的布局型式。

    圖4 X-43A 高超聲速飛行器PAI 布局[23]Fig.4 Airframe/propulsion integrated aerodynamic configuration of X-43A hypersonic vehicle[23]

    前體與進(jìn)氣道一體化設(shè)計(jì)不僅僅考慮發(fā)動(dòng)機(jī)推力和流場(chǎng)要求,還需綜合考慮氣動(dòng)力、結(jié)構(gòu)和內(nèi)部容積等其他限制條件。

    機(jī)體/推進(jìn)一體化的另一特點(diǎn)是機(jī)體的后體被用作發(fā)動(dòng)機(jī)的延伸部分。 后體噴管的作用使發(fā)動(dòng)機(jī)出口氣流得到進(jìn)一步膨脹產(chǎn)生的推力貢獻(xiàn),這是由于發(fā)動(dòng)機(jī)循環(huán)性能改善引起的,其設(shè)計(jì)性能關(guān)鍵參數(shù)包括:

    1)捕獲面積。 前體對(duì)進(jìn)入發(fā)動(dòng)機(jī)氣流提供預(yù)壓縮,提高了發(fā)動(dòng)機(jī)的推進(jìn)效率。 進(jìn)氣道特性可用兩個(gè)參數(shù)衡量,一是發(fā)動(dòng)機(jī)捕獲的質(zhì)量流量比,即A1/A0要大,二是進(jìn)氣道的壓力恢復(fù)要高。高超聲速飛行要求Scramjet 發(fā)動(dòng)機(jī)產(chǎn)生足夠大推力,因此,前體進(jìn)氣道一體化設(shè)計(jì)中,必須考慮到增大進(jìn)氣道氣流的捕獲面積。 高度集成的高超聲速飛行器設(shè)計(jì)概念吸引力在于它可以提供最大的進(jìn)氣道捕獲面積和最大的噴管膨脹面積。

    如式(5)所示,吸氣式發(fā)動(dòng)機(jī)的推力F與前體捕獲的氣流質(zhì)量流量直接相關(guān)[24]:

    式中,ρ0V0是飛行器前方未受擾動(dòng)自由流每單位面積質(zhì)量流速率,A1是發(fā)動(dòng)機(jī)捕獲面積。

    發(fā)動(dòng)機(jī)進(jìn)氣道置于飛行器下表面,位于前體頭部弓形激波產(chǎn)生的預(yù)壓縮面內(nèi),進(jìn)一步提高了發(fā)動(dòng)機(jī)進(jìn)氣道入口處每單位面積質(zhì)量流速率ρ0V0,它對(duì)全機(jī)推力貢獻(xiàn)約占總推力的35%[24]。圖5 為X-43A 的三級(jí)前體斜板和進(jìn)氣道捕獲面積示意圖。

    圖5 三級(jí)前體斜板和進(jìn)氣道捕獲面積示意圖[24]Fig.5 Schematic diagram of X-43A three-stage prea cu[2r4s]o r inclined plate and inlet captured are-

    從式(5)可以看出,捕獲面積增加,可增加捕獲流量進(jìn)而增加發(fā)動(dòng)機(jī)推力;但捕獲流量增加必增大前體面積,導(dǎo)致氣動(dòng)阻力、抬頭力矩及靜不穩(wěn)定性均增加。 因此,必須從總體性能上綜合考慮,當(dāng)捕獲面積在一定范圍內(nèi)增加時(shí),發(fā)動(dòng)機(jī)推力增加更明顯,推阻裕量變大;若捕獲面積超過(guò)一定值,氣動(dòng)阻力增加更明顯,推阻裕量變小。 捕獲面積的最優(yōu)值需結(jié)合總體性能來(lái)確定。

    2)前體最大寬度。 根據(jù)發(fā)動(dòng)機(jī)運(yùn)行要求,進(jìn)入發(fā)動(dòng)機(jī)的入流品質(zhì)應(yīng)最佳,沿展向氣流所有參數(shù)變化梯度小,即畸變度小。 飛行器的前體寬度和發(fā)動(dòng)機(jī)的相對(duì)寬度對(duì)進(jìn)氣道入口的均勻性影響很大,如果發(fā)動(dòng)機(jī)質(zhì)量流量梯度很大,就不可能獲得發(fā)動(dòng)機(jī)的最佳性能。 發(fā)動(dòng)機(jī)尺寸和流場(chǎng)要求,對(duì)于良好的前體設(shè)計(jì)不是唯一必須考慮的因素,對(duì)于最佳的前體設(shè)計(jì)必需綜合考慮氣動(dòng)力、結(jié)構(gòu)等多方面的要求。

    3)前體長(zhǎng)度。 X-43A 只能在十分狹小的吸氣走廊范圍(見(jiàn)圖6 所示的X-43A 速度-高度曲線)沿恒值動(dòng)壓q0的飛行軌跡飛行,其相應(yīng)動(dòng)壓范圍[500, 2000](psf),其下邊界受最大允許升力面面積限制,上邊界受結(jié)構(gòu)極限限制。 與單位面積質(zhì)量流速率、當(dāng)?shù)芈曀賏0、來(lái)流速度V0、來(lái)流馬赫數(shù)M0有式(7)所示關(guān)系:

    圖6 X-43A 速度-高度曲線[25]Fig.6 Velocity-altitude curve for X-43A[25]

    從式(7)中可以看出,隨M0增大,每單位面積質(zhì)量流速率ρ0V0迅速下降(X-43A 速度-質(zhì)量流速率曲線如圖7 所示),雖然這與大氣密度減少有關(guān),然而這樣的變化趨勢(shì)與發(fā)動(dòng)機(jī)工作需求是相違背的。 為維持質(zhì)量流量和推力不變,就必須在設(shè)計(jì)上采取相應(yīng)措施,譬如,高超聲速飛行器使用很長(zhǎng)的前體壓縮面,其最重要作用就是增大自由流捕獲面積;或在特定M0下,增大動(dòng)壓q0,也可維持推力不變。

    圖7 X-43A 速度-質(zhì)量流速率曲線[25]Fig.7 Velocity-mass flow rate curve for X-43A[25]

    4)前體下表面斜板傾斜角。 前體是產(chǎn)生飛行器升力的主要部件,利用機(jī)體下表面三級(jí)前體斜板和進(jìn)氣道產(chǎn)生的壓力升力和有限的前體阻力設(shè)計(jì)可有效提高飛行器升阻比。 借助前體弓形激波和前體設(shè)計(jì)可為升力體/乘波體飛行器提供足夠升力,因此,前體下表面斜板傾斜角的合理設(shè)計(jì),是前體進(jìn)氣道一體化設(shè)計(jì)的另一個(gè)重要參數(shù),其對(duì)發(fā)動(dòng)機(jī)流量系數(shù)、全機(jī)升阻比特性及設(shè)計(jì)條件下全機(jī)配平特性都有重要影響。

    為協(xié)調(diào)這些不同要求,前體下表面斜板從機(jī)頭到進(jìn)氣道口可分成若干段,稱(chēng)多級(jí)斜板,每級(jí)斜板相對(duì)來(lái)流都是小角度偏轉(zhuǎn)。 如圖5 所示,X-43A 用三級(jí)前體斜板,其偏角δ1=5°,δ2=6°,δ3=3° 。 既協(xié)調(diào)了不同需求,也滿足了在進(jìn)氣道入口處所要求的減速、增壓氣流條件。 相比而言,如果為滿足全機(jī)升力特性要求必須使用δ1=5°的前體斜板,而不用多級(jí)斜板設(shè)計(jì);為滿足在巡航設(shè)計(jì)條件下,弓形激波搭在進(jìn)氣道唇緣上,而不產(chǎn)生進(jìn)氣道溢流阻力,前體長(zhǎng)度需進(jìn)一步向前延伸,盡管流量系數(shù)增加了,但會(huì)導(dǎo)致飛行器結(jié)構(gòu)重量增加、柔性變形大、前體壓力中心前移,可能導(dǎo)致配平損失增加、飛行器性能惡化。 前體流場(chǎng)品質(zhì)對(duì)發(fā)動(dòng)機(jī)流場(chǎng)品質(zhì)也有重要影響,進(jìn)入燃燒室的氣流流場(chǎng)均勻性將直接影響發(fā)動(dòng)機(jī)推力[26]。

    5)后體/噴管一體化設(shè)計(jì)。 PAI 另一個(gè)特點(diǎn)是機(jī)體后部部分被用作發(fā)動(dòng)機(jī)的延伸部分(或稱(chēng)外噴管)。 如圖4 所示,把進(jìn)氣道集成在前體下表面是為發(fā)動(dòng)機(jī)獲得最大質(zhì)量流量,后體-噴管設(shè)計(jì)也是為從膨脹的后體排氣流中獲得最大推力。 機(jī)體后部下表面作為單斜壁噴管的上邊界,其下邊界由周?chē)諝獾淖杂闪骷羟袑訕?gòu)成,剪切層位置由排氣羽流靜壓與自由靜壓相平衡決定。擴(kuò)張的外噴管使發(fā)動(dòng)機(jī)出口的高溫、高壓氣體進(jìn)一步膨脹,如果噴管長(zhǎng)度太短,不足以使儲(chǔ)存的位勢(shì)能得到充分恢復(fù),稱(chēng)“欠膨脹”,其結(jié)果是降低了推力。 一旦后體過(guò)度膨脹,致使后體壓力降低,從而使前、后體正向壓差增大,則使全機(jī)阻力增加。 NASP 研究表明,后體產(chǎn)生的推力貢獻(xiàn)約占總推力的25~35%[27-28]。 這種貢獻(xiàn)是由于改善了發(fā)動(dòng)機(jī)性能增大比沖引起的。

    后體噴管設(shè)計(jì)主要受推力和穩(wěn)定性支配。 為獲得全機(jī)的推力減阻力的差額盡可能大,則后體噴管的推力貢獻(xiàn)將起到關(guān)鍵作用。 飛行器后體作為噴管的膨脹面,作用著由發(fā)動(dòng)機(jī)排出的高壓排氣流,對(duì)后體產(chǎn)生的升力和俯仰力矩也有顯著影響。 因此,發(fā)動(dòng)機(jī)位置、推力矢量方向和所產(chǎn)生的配平損失都必須在整個(gè)飛行包線范圍內(nèi)進(jìn)行檢查。

    后體噴管提供的附加推力,基本上是在不增加阻力情況下獲得的“純凈推力”。 對(duì)高超聲速飛行器,全機(jī)“凈推力”只是全機(jī)總推力和全機(jī)總阻力兩個(gè)大量中很小的差量,所以任何阻力增加都可能顯著影響飛行器性能。 這種獨(dú)特的后體構(gòu)型設(shè)計(jì)對(duì)于全機(jī)凈推力和升阻比都很重要。

    高超聲速飛行器為了獲得大推力,要求很大的膨脹比,必然導(dǎo)致噴管出口面積很大。 為減輕重量,在實(shí)際設(shè)計(jì)中,去除噴管側(cè)壁和罩的實(shí)體部分,長(zhǎng)度也被截短,噴管的長(zhǎng)度必須在性能要求和綜合約束之間尋求折衷[29-30]。 另外,在后體下表面受發(fā)動(dòng)機(jī)高壓排氣流施加的力和力矩(圖5)可能很大,在布局設(shè)計(jì)中必須合理配置,如果設(shè)計(jì)不當(dāng)可能導(dǎo)致配平損失很大。

    3 關(guān)于機(jī)體/推進(jìn)系統(tǒng)一體化與機(jī)體氣動(dòng)力干擾評(píng)述

    Hyper-X 等吸氣式高超聲速飛行器利用Scramjet 進(jìn)行推進(jìn),發(fā)動(dòng)機(jī)位于機(jī)體下。 對(duì)于這類(lèi)推進(jìn)系統(tǒng)和機(jī)體高度一體化的飛行器的布局設(shè)計(jì)和使用,一個(gè)關(guān)鍵問(wèn)題是推進(jìn)系統(tǒng)和飛機(jī)機(jī)體空氣動(dòng)力學(xué)之間的相互干擾作用,該問(wèn)題主要與PAI 有關(guān)。 通常前體作為Scramjet 進(jìn)氣道,后體作為Scramjet 噴管。

    高超聲速飛行要求發(fā)動(dòng)機(jī)能夠產(chǎn)生足夠大的推力,而且與飛行器重心位置間的力臂很大,導(dǎo)致發(fā)動(dòng)機(jī)推力對(duì)俯仰力矩影響。 同時(shí),推進(jìn)系統(tǒng)對(duì)飛行器橫向操穩(wěn)特性也有重大影響[31-32]。

    吸氣式高超聲速飛行器的氣動(dòng)力也會(huì)影響推進(jìn)系統(tǒng)。 飛行器前體下表面弓形激波性質(zhì)確定了進(jìn)入進(jìn)氣道的壓縮氣流和前體進(jìn)氣道捕獲的空氣質(zhì)量流量,而弓形激波性質(zhì)是由氣流的迎角、動(dòng)壓和自由流特性確定的[33]。 氣流的迎角和動(dòng)壓也影響燃燒動(dòng)力學(xué)和排氣流/自由流的剪切層。

    吸氣式高超聲速飛行器機(jī)體和推進(jìn)系統(tǒng)動(dòng)態(tài)耦合來(lái)源于結(jié)構(gòu)的氣動(dòng)彈性模態(tài)、前體和后體及其在整個(gè)機(jī)體上的傳播,影響了通過(guò)進(jìn)氣道和排氣噴管的流量,因此也影響了氣動(dòng)力性能。 精確確定結(jié)構(gòu)的氣動(dòng)彈性模態(tài)對(duì)于飛行控制,尤其是精確控制迎角是很苛刻的要求。

    4 高超聲速飛行器地面試驗(yàn)要求與設(shè)備能力限制

    在高超聲速飛行器設(shè)計(jì)中,由于使用非傳統(tǒng)、高集成度布局,導(dǎo)致各部件之間相互影響靈敏度增強(qiáng),加之飛行速度范圍寬(Ma=3~25),飛行環(huán)境嚴(yán)酷。

    圖8 給出了幾類(lèi)高超聲速地面試驗(yàn)設(shè)備在X-43A 飛行走廊區(qū)域內(nèi)的模擬能力。 這些地面設(shè)備包括下吹式風(fēng)洞、暫沖式風(fēng)洞和脈沖風(fēng)洞。 暫沖式風(fēng)洞設(shè)備運(yùn)行持續(xù)時(shí)間較長(zhǎng),但焓值增加受到限制,最大馬赫數(shù)不超過(guò)8Ma,實(shí)際可達(dá)到的Re數(shù)以及溫度都受到很大限制,風(fēng)洞總溫和試驗(yàn)Re數(shù)都低于飛行條件;脈沖加熱設(shè)備,例如激波風(fēng)洞,可以短時(shí)間模擬更高焓值,最高馬赫數(shù)可以達(dá)到10Ma或更高,但為減少高能氣體損耗并防護(hù)噴管燒蝕,其運(yùn)行時(shí)間極短,以毫秒計(jì)算,這就限制了激波風(fēng)洞性能,無(wú)法模擬更多飛行狀態(tài)下熱力學(xué)參數(shù),無(wú)法測(cè)熱流及火焰穩(wěn)定性。 對(duì)于運(yùn)行時(shí)間較長(zhǎng)的設(shè)備,速度為8Ma或小于8Ma,風(fēng)洞設(shè)備質(zhì)量流量能力限制了試驗(yàn)段尺寸,也就限制了試驗(yàn)?zāi)P统叽纭?對(duì)于激波加熱設(shè)備,速度為10Ma或更高,試驗(yàn)持續(xù)時(shí)間限制了試驗(yàn)件長(zhǎng)度[23]。

    圖8 X-43A 飛行走廊及目前地面設(shè)備模擬能力圖[23]Fig. 8 Simulation capability diagram of flight corridor and current ground equipment for X-43A[23]

    根據(jù)目前美國(guó)地面設(shè)備的尺寸能力和設(shè)備持續(xù)運(yùn)行時(shí)間可以推斷:現(xiàn)有的地面設(shè)備只適合于小型高超聲速飛行器試驗(yàn),不能滿足全尺寸飛行器性能評(píng)估,也不適用于“全球到達(dá)”高超聲速飛行器的試驗(yàn)工作。 對(duì)于高超聲速飛行器PAI 的風(fēng)洞試驗(yàn)驗(yàn)證,由于受現(xiàn)有地面設(shè)備尺寸限制,試驗(yàn)?zāi)P统叽缧。荒苓M(jìn)行帶Scramjet 發(fā)動(dòng)機(jī)運(yùn)行的、完整的一體化布局試驗(yàn),因此,無(wú)法獲得對(duì)“一體化”飛行器性能有重大影響的機(jī)體和發(fā)動(dòng)機(jī)系統(tǒng)之間全部干擾作用。

    當(dāng)前高超聲速飛行器大量風(fēng)洞試驗(yàn)是采用縮尺模型進(jìn)行的,但高超聲速飛行器縮尺模型試驗(yàn)存在兩個(gè)問(wèn)題值得特別注意:①縮尺模型性能和適應(yīng)性不同于全尺寸試驗(yàn)件,這里所指的“性能”是指用氣流流量、燃油流量和推力表示的總推力性能。 這說(shuō)明用風(fēng)洞縮尺模型試驗(yàn)評(píng)估得到的推力性能和全尺寸試驗(yàn)得來(lái)的總推力性能是不同的[34];②受風(fēng)洞總溫、總壓等條件限制,縮尺模型和全尺寸試驗(yàn)件的自然現(xiàn)象和流態(tài)會(huì)有明顯差異,沿著前體下表面流路將導(dǎo)致邊界層厚度、熵層厚度、粘性層厚度、層流到湍流的轉(zhuǎn)捩位置、激波和邊界層干擾以及激波和激波干擾等流動(dòng)特性發(fā)生變化,特別是當(dāng)邊界層流入發(fā)動(dòng)機(jī)內(nèi)部,由于湍流摻混和化學(xué)反應(yīng)相互作用,也將影響燃燒過(guò)程和化學(xué)組分。

    地面推進(jìn)系統(tǒng)試驗(yàn)只能利用控制體積和參考平面將把發(fā)動(dòng)機(jī)系統(tǒng)從飛行器上分割出來(lái),分解成無(wú)連系的“可試驗(yàn)單元”,這種試驗(yàn)方法可能涉及兩個(gè)問(wèn)題:①分界面上的邊界條件必須用高保真度方法仿真計(jì)算,特別在中、高的高超聲速條件下(即Ma≥10)下,燃燒室或噴管的進(jìn)口條件依靠現(xiàn)有地面設(shè)備進(jìn)行仿真試驗(yàn)是很難實(shí)現(xiàn)的;②凈推力很難在地面試驗(yàn)中測(cè)定。 因?yàn)閮敉屏κ强偼屏涂傋枇蓚€(gè)大量中的小差量,由每個(gè)可試驗(yàn)單元估算整個(gè)系統(tǒng)總性能的過(guò)程有太多不確定性,如果估算誤差為1%~2%,都將使凈推力產(chǎn)生很大變化[35],因此,凈推力只能由飛行試驗(yàn)測(cè)定。

    與大氣中飛行相比,所有地面設(shè)備氣體成分都有污染,這種污染會(huì)影響發(fā)動(dòng)機(jī)燃燒過(guò)程,包括點(diǎn)火特性、火焰穩(wěn)定性和燃料混合特性等,其終將對(duì)發(fā)動(dòng)機(jī)推力貢獻(xiàn)產(chǎn)生影響。 X-43A 在速度為7Ma飛行時(shí),為降低飛行風(fēng)險(xiǎn),在不同設(shè)備中進(jìn)行了重復(fù)性試驗(yàn),辨識(shí)出了數(shù)據(jù)反常是由于設(shè)備污染所致[35]。

    地面試驗(yàn)的另外一個(gè)限制體現(xiàn)在研究?jī)?nèi)容上,Scramjet 發(fā)動(dòng)機(jī)的試驗(yàn)研究幾乎都集中在發(fā)動(dòng)機(jī)推力、部件效率和燃燒效率等方面。 實(shí)際上,帶動(dòng)力條件下,高超聲速吸氣飛行器的氣動(dòng)力分析和試驗(yàn)是PAI 的一個(gè)重要方面。 但設(shè)備尺寸限制,且燃燒過(guò)程不是幾何可縮比,排除了帶Scramjet 運(yùn)行的可能性。

    實(shí)現(xiàn)Scramjet 發(fā)動(dòng)機(jī)縮尺試驗(yàn),必須處理兩個(gè)問(wèn)題:①反應(yīng)時(shí)間(或長(zhǎng)度),②混合長(zhǎng)度。 由于燃料注射、點(diǎn)火和燃燒機(jī)理很復(fù)雜,而且湍流度模擬也是很不確定的,因此需進(jìn)行全尺寸的Scramjet 發(fā)動(dòng)機(jī)地面試驗(yàn)。

    5 增量方法論發(fā)展

    5.1 增量方法論

    高超聲速飛行器PAI 布局是一種非傳統(tǒng)布局型式,機(jī)體和推進(jìn)系統(tǒng)之間強(qiáng)耦合,必須評(píng)估推進(jìn)系統(tǒng)和機(jī)體之間相互干擾作用及其對(duì)飛行器性能影響。 為了滿足這類(lèi)飛行器研制需要,必須綜合利用現(xiàn)有地面設(shè)備、飛行試驗(yàn)以及大型計(jì)算機(jī)開(kāi)發(fā)出適用的新設(shè)計(jì)方法。 Boudreau 等提出“增量方法論”,即依靠地面試驗(yàn)和計(jì)算的增量,再與基準(zhǔn)的流動(dòng)相加[36]。 這種一體化方法的關(guān)鍵是在增量形式下,采用CFD 方法把風(fēng)洞試驗(yàn)環(huán)境下的氣動(dòng)力數(shù)據(jù)外推到飛行環(huán)境的數(shù)據(jù),可表示為式(8):

    式中,下標(biāo)飛行指飛行條件,基準(zhǔn)指風(fēng)洞試驗(yàn)基準(zhǔn)值,Test 指風(fēng)洞試驗(yàn)值,CFD 指CFD 計(jì)算值。

    這種方法論既立足于現(xiàn)有風(fēng)洞試驗(yàn)值并能計(jì)及由于地面試驗(yàn)環(huán)境與飛行環(huán)境差異,地面模擬不足產(chǎn)生的增量差。 文中所指“增量”意義是明確的,不是單獨(dú)依靠CFD 代碼絕對(duì)值外插能力,因此不受當(dāng)前CFD 在絕對(duì)意義上的模擬分離流、湍流度和化學(xué)作用的弱點(diǎn)限制。

    飛行器基準(zhǔn)數(shù)據(jù)應(yīng)在高質(zhì)量、完全氣體高超聲速風(fēng)洞中試驗(yàn)確定,并用增量設(shè)計(jì)代碼修正與設(shè)備有關(guān)試驗(yàn)條件引起的偏差,如氣流非均勻性、支撐系統(tǒng)影響等。 此外,在開(kāi)發(fā)增量方法中,第一步工作就是核實(shí)CFD 工具,要求其包括足夠物理和化學(xué)能力,能精確預(yù)測(cè)增量效應(yīng),計(jì)及地面試驗(yàn)和飛行之間差異。

    NASA 的Hyper-X 計(jì)劃,亦稱(chēng)高超聲速增量技術(shù)計(jì)劃,目的是驗(yàn)證和確認(rèn)PAI 飛行器研制所需的試驗(yàn)技術(shù)、計(jì)算方法和設(shè)計(jì)工具。 為達(dá)此目的,NASA 還開(kāi)展了帶Scramjet 動(dòng)力的飛行試驗(yàn),利用飛行試驗(yàn)數(shù)據(jù)驗(yàn)證設(shè)計(jì)預(yù)測(cè)值和設(shè)計(jì)方法論。 依托Hyper-X 計(jì)劃,高超聲速飛行器研制歷史上第一次建立了PAI 布局分析、CFD 數(shù)值計(jì)算、地面風(fēng)洞試驗(yàn)以及飛行試驗(yàn)數(shù)據(jù)的相關(guān)性研究,更重要的是X-43A 飛行數(shù)據(jù)驗(yàn)證了基于風(fēng)洞試驗(yàn)和計(jì)算流體力學(xué)給出的預(yù)測(cè)數(shù)據(jù),也驗(yàn)證了高超聲速飛行器設(shè)計(jì)方法,使得高超聲速飛行獲得成功。

    5.2 數(shù)據(jù)庫(kù)開(kāi)發(fā)

    NASA 為分析X-43A 的氣動(dòng)特性,在飛行器設(shè)計(jì)初期方案設(shè)計(jì)階段使用了分析或數(shù)值方法評(píng)估,最終性能評(píng)估則是利用基于風(fēng)洞試驗(yàn)數(shù)據(jù)建立的數(shù)據(jù)庫(kù)方法。 該數(shù)據(jù)庫(kù)發(fā)展流程如圖9 所示,圖中數(shù)據(jù)庫(kù)由風(fēng)洞試驗(yàn)數(shù)據(jù)、CFD 數(shù)值分析及飛行數(shù)據(jù)組成。

    圖9 Hyper-X 氣動(dòng)力數(shù)據(jù)庫(kù)發(fā)展過(guò)程Fig. 9 Development process of Hyper-X aerodynamic database

    基于風(fēng)洞數(shù)據(jù)進(jìn)行計(jì)算機(jī)分析的目的在于:①對(duì)非飛行試驗(yàn)件的地面風(fēng)洞試驗(yàn)數(shù)據(jù)進(jìn)行修正以計(jì)入環(huán)境、試驗(yàn)件和試驗(yàn)持續(xù)時(shí)間等影響因素;②預(yù)測(cè)飛行包線范圍內(nèi)馬赫數(shù)變化影響,包括內(nèi)插獲得非試驗(yàn)點(diǎn)性能。

    由于風(fēng)洞模型尺寸相對(duì)較小,不可能進(jìn)行進(jìn)氣道打開(kāi)試驗(yàn)(包括無(wú)動(dòng)力或帶動(dòng)力的),因此,現(xiàn)有風(fēng)洞試驗(yàn)氣動(dòng)力試驗(yàn)數(shù)據(jù)庫(kù)都是用關(guān)閉進(jìn)氣口的幾何外形獲得的,然而,缺少帶動(dòng)力風(fēng)洞試驗(yàn)數(shù)據(jù)是不足以建構(gòu)一個(gè)完整氣動(dòng)力數(shù)據(jù)庫(kù)的。 對(duì)于進(jìn)氣道打開(kāi),發(fā)動(dòng)機(jī)點(diǎn)火情況下氣動(dòng)力增量,利用CFD 技術(shù),把計(jì)算結(jié)果以相對(duì)于數(shù)據(jù)增量形式引入數(shù)據(jù)庫(kù),可見(jiàn),這樣的氣動(dòng)力數(shù)據(jù)庫(kù)基礎(chǔ)仍然是風(fēng)洞試驗(yàn)數(shù)據(jù)。 NASA Langley 8 英尺高溫風(fēng)洞全尺寸推進(jìn)流道試驗(yàn)證實(shí)了該方法所得到預(yù)測(cè)結(jié)果的正確性[37]。

    此外,利用CFD 數(shù)值分析方法研究還可去除風(fēng)洞試驗(yàn)影響,例如Re數(shù)的尺度效應(yīng)、Scramjet推進(jìn)流場(chǎng)誘導(dǎo)效應(yīng)(這與進(jìn)氣口打開(kāi)引起的力和力矩變化有關(guān))以及風(fēng)洞支撐干擾影響等。 飛行數(shù)據(jù)也會(huì)被用來(lái)修正和改善Hyper-X 飛行器的氣動(dòng)力數(shù)據(jù)庫(kù)。

    6 結(jié)論

    通過(guò)分析評(píng)述類(lèi)X-43A 吸氣式高超聲速飛行器PAI 相關(guān)設(shè)計(jì)、“一體化”構(gòu)型地面設(shè)備要求、風(fēng)洞能力和限制以及增量方法論發(fā)展,可得出如下結(jié)論:

    1)對(duì)于X-43A 這類(lèi)PAI 布局的高超聲速飛行器,傳統(tǒng)飛機(jī)設(shè)計(jì)方法、試驗(yàn)方法和理論分析方法都有其限制性。 為了滿足這類(lèi)飛行器研制需要,必須綜合利用現(xiàn)有地面設(shè)備、飛行試驗(yàn)以及大型計(jì)算機(jī),依靠地面試驗(yàn)和計(jì)算的增量,再與基準(zhǔn)的流動(dòng)相加,這就是“增量方法論”。 這種“一體化”設(shè)計(jì)方法論的關(guān)鍵點(diǎn)是:在增量形式下,采用CFD 方法把風(fēng)洞試驗(yàn)環(huán)境下的氣動(dòng)力數(shù)據(jù)外推到飛行環(huán)境的數(shù)據(jù)。

    2)由于飛行器氣動(dòng)界面與推進(jìn)系統(tǒng)劃分方式不同會(huì)導(dǎo)致升阻比不同,因此,高超聲速飛行器概念方案設(shè)計(jì)階段,單獨(dú)采用L/D或Isp來(lái)評(píng)估飛行器氣動(dòng)外形或發(fā)動(dòng)機(jī)性能是不合理的,而必須對(duì)(L/D)×Isp實(shí)行優(yōu)化。

    3)高超聲速飛行器的前體是產(chǎn)生升力的主要部件,利用壓力升力和有限阻力的前體設(shè)計(jì)是提高升阻比的主要途徑;后體-噴管一體化設(shè)計(jì)對(duì)于提高全機(jī)升阻比,增大凈推力具有重要意義;同時(shí),后體-噴管一體化設(shè)計(jì)還必須考慮配平影響,減輕配平阻力損失;

    4)由于受設(shè)備尺寸能力和持續(xù)運(yùn)行時(shí)間短的限制,現(xiàn)有地面風(fēng)洞設(shè)備只能實(shí)現(xiàn)小型、低高超聲速馬赫數(shù)試驗(yàn),全尺寸或高馬赫數(shù)實(shí)驗(yàn)則超出了現(xiàn)有設(shè)備能力。

    5)美國(guó)現(xiàn)有設(shè)備不能實(shí)現(xiàn)PAI 布局帶動(dòng)力運(yùn)行的完整一體化布局試驗(yàn)。 這意味著現(xiàn)有設(shè)備不能滿足具有“一體化”布局特點(diǎn)的所有試驗(yàn)要求,包括機(jī)體與推進(jìn)系統(tǒng)的相互干擾、發(fā)動(dòng)機(jī)位置影響以及后體設(shè)計(jì)影響等。

    6)現(xiàn)有設(shè)備對(duì)于推進(jìn)系統(tǒng)試驗(yàn)有更多限制,無(wú)法測(cè)定前、后流路對(duì)發(fā)動(dòng)機(jī)影響,無(wú)法確定發(fā)動(dòng)機(jī)“凈推力”,無(wú)法確定試驗(yàn)氣體成分污染影響。

    猜你喜歡
    進(jìn)氣道風(fēng)洞試驗(yàn)超聲速
    高超聲速出版工程
    高超聲速飛行器
    基于AVL-Fire的某1.5L發(fā)動(dòng)機(jī)進(jìn)氣道優(yōu)化設(shè)計(jì)
    基于輔助進(jìn)氣門(mén)的進(jìn)氣道/發(fā)動(dòng)機(jī)一體化控制
    超聲速旅行
    低風(fēng)壓架空導(dǎo)線的風(fēng)洞試驗(yàn)
    電線電纜(2017年5期)2017-10-18 00:52:03
    滾轉(zhuǎn)機(jī)動(dòng)載荷減緩風(fēng)洞試驗(yàn)
    The coupling characteristics of supersonic dual inlets for missile①
    遮擋條件下超高層建筑風(fēng)洞試驗(yàn)研究
    重慶建筑(2014年12期)2014-07-24 14:00:32
    高超聲速大博弈
    太空探索(2014年5期)2014-07-12 09:53:28
    水蜜桃什么品种好| 国产精品一区二区免费欧美 | 国产黄频视频在线观看| 人人妻人人爽人人添夜夜欢视频| 国产视频一区二区在线看| 亚洲欧洲日产国产| 狠狠婷婷综合久久久久久88av| 久久久久久久精品精品| 99久久人妻综合| 别揉我奶头~嗯~啊~动态视频 | 在线观看www视频免费| 欧美97在线视频| 老汉色av国产亚洲站长工具| 999精品在线视频| 国产成人欧美在线观看 | 女人高潮潮喷娇喘18禁视频| 国产老妇伦熟女老妇高清| 日韩中文字幕欧美一区二区 | 男女国产视频网站| 狂野欧美激情性bbbbbb| 波多野结衣一区麻豆| 精品国产一区二区三区久久久樱花| 亚洲av男天堂| 国产片特级美女逼逼视频| 国产免费福利视频在线观看| 999精品在线视频| 国产精品一区二区免费欧美 | 男女无遮挡免费网站观看| 青草久久国产| 在线观看www视频免费| 国产成人av教育| 国产欧美日韩一区二区三 | 欧美性长视频在线观看| 天堂8中文在线网| 丝袜人妻中文字幕| 国产在线视频一区二区| 亚洲,欧美,日韩| 国产亚洲欧美在线一区二区| 欧美在线黄色| 日韩视频在线欧美| 久久亚洲精品不卡| 在线观看人妻少妇| 午夜免费鲁丝| 亚洲国产欧美在线一区| 国产麻豆69| 国产亚洲欧美精品永久| 我的亚洲天堂| 夜夜骑夜夜射夜夜干| 中文字幕精品免费在线观看视频| 性色av一级| 精品人妻1区二区| 女警被强在线播放| 免费在线观看黄色视频的| 久9热在线精品视频| 亚洲色图 男人天堂 中文字幕| 大香蕉久久网| 国产一区有黄有色的免费视频| 丰满饥渴人妻一区二区三| 亚洲精品久久午夜乱码| 国产亚洲欧美精品永久| 如日韩欧美国产精品一区二区三区| 在线观看免费日韩欧美大片| 成年美女黄网站色视频大全免费| 99国产精品一区二区蜜桃av | 水蜜桃什么品种好| 少妇猛男粗大的猛烈进出视频| 一区二区三区四区激情视频| 久久综合国产亚洲精品| 久久久精品94久久精品| 在线av久久热| 国产免费视频播放在线视频| 精品人妻一区二区三区麻豆| 国产成人精品无人区| 亚洲av男天堂| 婷婷丁香在线五月| 国产成人啪精品午夜网站| 18禁黄网站禁片午夜丰满| 久久天躁狠狠躁夜夜2o2o | e午夜精品久久久久久久| 亚洲人成电影免费在线| 日本一区二区免费在线视频| 精品久久久久久电影网| 色综合欧美亚洲国产小说| 不卡av一区二区三区| 丁香六月天网| 19禁男女啪啪无遮挡网站| 满18在线观看网站| 国产成人一区二区在线| 91麻豆精品激情在线观看国产 | 性高湖久久久久久久久免费观看| 亚洲久久久国产精品| 免费高清在线观看日韩| 国产男女内射视频| 久热爱精品视频在线9| 1024香蕉在线观看| 天堂中文最新版在线下载| 国产极品粉嫩免费观看在线| 老鸭窝网址在线观看| 国产亚洲欧美在线一区二区| 精品人妻一区二区三区麻豆| 久久久久久人人人人人| 啦啦啦啦在线视频资源| 午夜视频精品福利| 成年动漫av网址| 欧美在线一区亚洲| 国产成人一区二区三区免费视频网站 | 人体艺术视频欧美日本| 人妻一区二区av| 免费看不卡的av| 久久毛片免费看一区二区三区| 日本欧美视频一区| 亚洲视频免费观看视频| 久久精品国产亚洲av涩爱| 黑丝袜美女国产一区| 香蕉丝袜av| 精品免费久久久久久久清纯 | 亚洲av成人不卡在线观看播放网 | 男的添女的下面高潮视频| 成人亚洲精品一区在线观看| av有码第一页| 热re99久久精品国产66热6| 黄网站色视频无遮挡免费观看| 丝袜脚勾引网站| 一区二区三区精品91| 精品免费久久久久久久清纯 | 亚洲av在线观看美女高潮| 日韩中文字幕欧美一区二区 | 母亲3免费完整高清在线观看| 51午夜福利影视在线观看| 亚洲精品久久久久久婷婷小说| 日本vs欧美在线观看视频| 看免费av毛片| 国产精品国产三级国产专区5o| 日韩av免费高清视频| 亚洲 欧美一区二区三区| 啦啦啦啦在线视频资源| 国产又爽黄色视频| 人体艺术视频欧美日本| 精品久久久久久电影网| 天天躁日日躁夜夜躁夜夜| tube8黄色片| 天天添夜夜摸| 69精品国产乱码久久久| 久久综合国产亚洲精品| 亚洲精品一卡2卡三卡4卡5卡 | 王馨瑶露胸无遮挡在线观看| 好男人电影高清在线观看| 亚洲av电影在线进入| 精品国产乱码久久久久久小说| 脱女人内裤的视频| 久久精品久久久久久久性| 日韩精品免费视频一区二区三区| 色网站视频免费| 日韩伦理黄色片| 久久综合国产亚洲精品| 欧美日韩黄片免| 老熟女久久久| 夜夜骑夜夜射夜夜干| 三上悠亚av全集在线观看| 久久亚洲精品不卡| 免费女性裸体啪啪无遮挡网站| 久久99精品国语久久久| 国产精品亚洲av一区麻豆| 老司机亚洲免费影院| 97人妻天天添夜夜摸| 青春草视频在线免费观看| 只有这里有精品99| 我要看黄色一级片免费的| 亚洲人成网站在线观看播放| 极品人妻少妇av视频| 亚洲国产精品成人久久小说| 亚洲人成77777在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 欧美成人午夜精品| 国产亚洲精品久久久久5区| 久久ye,这里只有精品| 中文字幕亚洲精品专区| a级毛片在线看网站| 啦啦啦中文免费视频观看日本| 国产成人系列免费观看| 精品国产乱码久久久久久小说| 欧美国产精品va在线观看不卡| 91麻豆精品激情在线观看国产 | 免费av中文字幕在线| 黄色 视频免费看| 久久中文字幕一级| 黄片播放在线免费| 日韩 欧美 亚洲 中文字幕| 1024视频免费在线观看| 天天躁日日躁夜夜躁夜夜| 国产精品国产av在线观看| 少妇人妻 视频| 777久久人妻少妇嫩草av网站| 精品人妻熟女毛片av久久网站| 日本wwww免费看| 久久精品国产综合久久久| 下体分泌物呈黄色| 国产高清videossex| 国产精品一区二区精品视频观看| 黄片小视频在线播放| 久久天躁狠狠躁夜夜2o2o | 亚洲欧美色中文字幕在线| 男女之事视频高清在线观看 | 黄色视频在线播放观看不卡| 老汉色∧v一级毛片| 女性生殖器流出的白浆| 十分钟在线观看高清视频www| 国产老妇伦熟女老妇高清| a级片在线免费高清观看视频| 欧美亚洲 丝袜 人妻 在线| 啦啦啦在线免费观看视频4| 亚洲国产精品国产精品| 成人免费观看视频高清| 亚洲av日韩在线播放| 菩萨蛮人人尽说江南好唐韦庄| 免费在线观看视频国产中文字幕亚洲 | 老司机深夜福利视频在线观看 | www.999成人在线观看| 欧美人与性动交α欧美精品济南到| 一级毛片电影观看| 飞空精品影院首页| 少妇被粗大的猛进出69影院| 国产精品成人在线| 91精品国产国语对白视频| 国产精品久久久久久精品电影小说| 熟女av电影| 亚洲欧美成人综合另类久久久| 黄色毛片三级朝国网站| 久久久久久久国产电影| 国产亚洲精品第一综合不卡| 久久久久精品人妻al黑| 国产在线免费精品| 这个男人来自地球电影免费观看| 人人澡人人妻人| 久久国产精品影院| 蜜桃在线观看..| 中文欧美无线码| 欧美另类一区| 99国产精品99久久久久| 日日摸夜夜添夜夜爱| 老汉色∧v一级毛片| 日韩精品免费视频一区二区三区| 精品一区二区三区四区五区乱码 | 国产午夜精品一二区理论片| 老司机午夜十八禁免费视频| 99热全是精品| 国产成人精品久久二区二区91| 欧美精品av麻豆av| 老司机在亚洲福利影院| 国产极品粉嫩免费观看在线| 777米奇影视久久| 国产成人av激情在线播放| 丝袜在线中文字幕| 久久av网站| 纵有疾风起免费观看全集完整版| 观看av在线不卡| 亚洲色图综合在线观看| 天天影视国产精品| 久久中文字幕一级| 色婷婷av一区二区三区视频| 一本久久精品| 亚洲视频免费观看视频| 狂野欧美激情性xxxx| 另类精品久久| 国产精品久久久久久人妻精品电影 | 国产一级毛片在线| cao死你这个sao货| 亚洲欧美清纯卡通| av有码第一页| 国产视频首页在线观看| 搡老岳熟女国产| 最新的欧美精品一区二区| 国产成人a∨麻豆精品| av视频免费观看在线观看| 视频区图区小说| 欧美精品亚洲一区二区| 亚洲七黄色美女视频| 色视频在线一区二区三区| 搡老岳熟女国产| 免费观看a级毛片全部| 国产精品99久久99久久久不卡| 成人亚洲欧美一区二区av| 精品少妇久久久久久888优播| 久久久精品免费免费高清| 一本—道久久a久久精品蜜桃钙片| 免费在线观看影片大全网站 | 2021少妇久久久久久久久久久| 日韩,欧美,国产一区二区三区| 免费av中文字幕在线| 国产在线观看jvid| 色婷婷av一区二区三区视频| 成人国产一区最新在线观看 | 波多野结衣av一区二区av| 青草久久国产| 中文字幕精品免费在线观看视频| 欧美黑人精品巨大| 麻豆av在线久日| 一级毛片电影观看| 热re99久久精品国产66热6| 纵有疾风起免费观看全集完整版| 婷婷色综合www| 国产一区二区激情短视频 | 天堂俺去俺来也www色官网| 一级毛片黄色毛片免费观看视频| 欧美另类一区| 丝袜在线中文字幕| 看十八女毛片水多多多| 老司机午夜十八禁免费视频| 亚洲一区中文字幕在线| 久久久精品免费免费高清| 这个男人来自地球电影免费观看| 精品少妇久久久久久888优播| 久久久久久久大尺度免费视频| 久久精品成人免费网站| 97精品久久久久久久久久精品| 一本综合久久免费| 欧美黑人精品巨大| 操出白浆在线播放| 日本猛色少妇xxxxx猛交久久| 高清欧美精品videossex| 国产又色又爽无遮挡免| 麻豆国产av国片精品| 色婷婷久久久亚洲欧美| 尾随美女入室| e午夜精品久久久久久久| 午夜福利影视在线免费观看| 亚洲欧洲精品一区二区精品久久久| 日韩中文字幕欧美一区二区 | 国产欧美日韩综合在线一区二区| 777久久人妻少妇嫩草av网站| 老汉色av国产亚洲站长工具| 亚洲综合色网址| 999久久久国产精品视频| 免费观看av网站的网址| 欧美大码av| 999久久久国产精品视频| 欧美精品人与动牲交sv欧美| 国产精品免费大片| 国产精品久久久久久精品古装| 久久女婷五月综合色啪小说| 成年女人毛片免费观看观看9 | 国产精品一二三区在线看| 国产精品二区激情视频| 国产熟女午夜一区二区三区| 老司机亚洲免费影院| 捣出白浆h1v1| 婷婷丁香在线五月| 国产成人欧美| 国产亚洲欧美在线一区二区| 精品欧美一区二区三区在线| 日韩 欧美 亚洲 中文字幕| 爱豆传媒免费全集在线观看| 精品亚洲成a人片在线观看| 午夜老司机福利片| 97在线人人人人妻| 少妇的丰满在线观看| 精品亚洲成国产av| 亚洲精品乱久久久久久| 男女无遮挡免费网站观看| 成年女人毛片免费观看观看9 | 免费观看a级毛片全部| 亚洲精品国产区一区二| 日本色播在线视频| 每晚都被弄得嗷嗷叫到高潮| 国产欧美日韩综合在线一区二区| 丝袜在线中文字幕| bbb黄色大片| 日韩制服丝袜自拍偷拍| 精品国产一区二区久久| 午夜免费成人在线视频| 中文字幕色久视频| h视频一区二区三区| 高清av免费在线| 国产真人三级小视频在线观看| xxx大片免费视频| 日韩视频在线欧美| 久久久精品94久久精品| 国产91精品成人一区二区三区 | 中文精品一卡2卡3卡4更新| 国产一区有黄有色的免费视频| 少妇的丰满在线观看| 免费少妇av软件| 成年女人毛片免费观看观看9 | 少妇人妻久久综合中文| videos熟女内射| 亚洲成国产人片在线观看| 无遮挡黄片免费观看| 爱豆传媒免费全集在线观看| 亚洲第一av免费看| 久久精品久久久久久噜噜老黄| 黄网站色视频无遮挡免费观看| 亚洲欧美精品自产自拍| 丁香六月欧美| 国产男女超爽视频在线观看| 91九色精品人成在线观看| 中文字幕人妻丝袜制服| 欧美日韩精品网址| 丝袜人妻中文字幕| 老司机靠b影院| 免费在线观看黄色视频的| 成人国产一区最新在线观看 | 国产一区有黄有色的免费视频| 又大又黄又爽视频免费| 午夜激情av网站| 女人久久www免费人成看片| 免费av中文字幕在线| 天天躁狠狠躁夜夜躁狠狠躁| 国产一区二区三区av在线| 亚洲欧美中文字幕日韩二区| 亚洲人成77777在线视频| 国产高清国产精品国产三级| 在线亚洲精品国产二区图片欧美| 欧美在线一区亚洲| 亚洲国产最新在线播放| 欧美激情 高清一区二区三区| 精品久久久久久电影网| 黄色 视频免费看| 日韩人妻精品一区2区三区| 人妻人人澡人人爽人人| 另类亚洲欧美激情| 在线观看免费午夜福利视频| 亚洲五月色婷婷综合| 捣出白浆h1v1| 看十八女毛片水多多多| 校园人妻丝袜中文字幕| 80岁老熟妇乱子伦牲交| 男女边吃奶边做爰视频| 欧美黑人精品巨大| 欧美日韩亚洲国产一区二区在线观看 | 日日夜夜操网爽| 欧美人与善性xxx| 亚洲国产精品一区三区| 亚洲国产看品久久| 亚洲免费av在线视频| 亚洲国产欧美日韩在线播放| 男的添女的下面高潮视频| 黄色毛片三级朝国网站| 久久久欧美国产精品| 可以免费在线观看a视频的电影网站| 男女无遮挡免费网站观看| 久久ye,这里只有精品| 久久久久精品人妻al黑| 美国免费a级毛片| 国产福利在线免费观看视频| 亚洲熟女毛片儿| 免费观看a级毛片全部| 亚洲伊人色综图| 婷婷丁香在线五月| 欧美激情 高清一区二区三区| 日日夜夜操网爽| 国产精品av久久久久免费| 老司机亚洲免费影院| 国产黄色免费在线视频| 久久久久久久国产电影| av又黄又爽大尺度在线免费看| av视频免费观看在线观看| 久久久久久久精品精品| 香蕉丝袜av| 最近手机中文字幕大全| 色94色欧美一区二区| 啦啦啦中文免费视频观看日本| 亚洲精品一二三| 亚洲精品国产一区二区精华液| 2018国产大陆天天弄谢| 自线自在国产av| 91国产中文字幕| 一级黄片播放器| 亚洲av国产av综合av卡| 亚洲精品成人av观看孕妇| 国产成人影院久久av| 一级a爱视频在线免费观看| 久久天躁狠狠躁夜夜2o2o | 国产伦理片在线播放av一区| 欧美少妇被猛烈插入视频| 青春草视频在线免费观看| 亚洲免费av在线视频| av线在线观看网站| 免费不卡黄色视频| www.自偷自拍.com| 日本午夜av视频| 婷婷色av中文字幕| 男人操女人黄网站| 亚洲五月色婷婷综合| 日韩一本色道免费dvd| 亚洲成国产人片在线观看| 国产高清视频在线播放一区 | 亚洲国产av影院在线观看| 在线亚洲精品国产二区图片欧美| 别揉我奶头~嗯~啊~动态视频 | 欧美日韩亚洲高清精品| 嫁个100分男人电影在线观看 | 欧美精品一区二区大全| 99国产精品一区二区三区| 成在线人永久免费视频| 日日摸夜夜添夜夜爱| 日本欧美视频一区| 国产成人91sexporn| 成年av动漫网址| 久久久国产一区二区| 国产亚洲av片在线观看秒播厂| 在线观看免费午夜福利视频| 亚洲欧美日韩另类电影网站| 一级片免费观看大全| 91九色精品人成在线观看| 国产福利在线免费观看视频| 日韩中文字幕欧美一区二区 | 美女午夜性视频免费| 天堂8中文在线网| 最近中文字幕2019免费版| 国产黄色视频一区二区在线观看| 国产国语露脸激情在线看| 久久国产精品男人的天堂亚洲| 亚洲av电影在线观看一区二区三区| 午夜福利免费观看在线| 一二三四在线观看免费中文在| 欧美日韩一级在线毛片| 美女中出高潮动态图| 脱女人内裤的视频| 国产成人精品无人区| 欧美亚洲日本最大视频资源| 精品久久蜜臀av无| 一二三四社区在线视频社区8| 精品卡一卡二卡四卡免费| 纵有疾风起免费观看全集完整版| 国产成人欧美| 久久久久国产一级毛片高清牌| 国产伦理片在线播放av一区| 一区二区三区乱码不卡18| 麻豆av在线久日| 自线自在国产av| 极品人妻少妇av视频| 精品人妻在线不人妻| 久久 成人 亚洲| 天天操日日干夜夜撸| 国产一级毛片在线| 午夜激情久久久久久久| 在线天堂中文资源库| 亚洲午夜精品一区,二区,三区| 亚洲熟女精品中文字幕| 国产真人三级小视频在线观看| xxxhd国产人妻xxx| 天堂中文最新版在线下载| 日韩一区二区三区影片| 操美女的视频在线观看| 成年动漫av网址| 久久 成人 亚洲| 欧美精品人与动牲交sv欧美| 丝瓜视频免费看黄片| 男女国产视频网站| 国产高清不卡午夜福利| 建设人人有责人人尽责人人享有的| 免费不卡黄色视频| 亚洲,欧美精品.| 咕卡用的链子| 久久九九热精品免费| 波多野结衣av一区二区av| 久久热在线av| 日本a在线网址| 欧美 日韩 精品 国产| 欧美日韩亚洲国产一区二区在线观看 | 天天躁夜夜躁狠狠久久av| 男人舔女人的私密视频| 亚洲国产中文字幕在线视频| 国产亚洲欧美在线一区二区| 亚洲精品第二区| 七月丁香在线播放| 99国产精品免费福利视频| 久久久久网色| 精品熟女少妇八av免费久了| av视频免费观看在线观看| 咕卡用的链子| 亚洲美女黄色视频免费看| 精品国产一区二区三区四区第35| 久久 成人 亚洲| 久久精品成人免费网站| 波多野结衣av一区二区av| 久久精品久久精品一区二区三区| 伊人亚洲综合成人网| 亚洲伊人色综图| 人人澡人人妻人| 大型av网站在线播放| 嫩草影视91久久| 亚洲国产精品一区三区| 国产免费福利视频在线观看| 成年av动漫网址| 制服人妻中文乱码| 天堂中文最新版在线下载| 999久久久国产精品视频| 国产视频首页在线观看| 国产亚洲精品久久久久5区| 午夜视频精品福利| 伦理电影免费视频| 女性生殖器流出的白浆| 成年人午夜在线观看视频| 久久这里只有精品19| 午夜福利,免费看| 国产老妇伦熟女老妇高清| 国产麻豆69| 久久亚洲国产成人精品v| 在线 av 中文字幕| 18禁黄网站禁片午夜丰满| 国产91精品成人一区二区三区 | 在线观看免费高清a一片| 婷婷色麻豆天堂久久| 免费少妇av软件| 国产成人精品久久二区二区免费| 久久性视频一级片| 日韩电影二区| 91九色精品人成在线观看| 69精品国产乱码久久久| av视频免费观看在线观看| 亚洲自偷自拍图片 自拍| 亚洲欧洲国产日韩| 人妻一区二区av| 下体分泌物呈黄色|