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

    XTER內(nèi)收縮組合進氣道設計理念及氣動特性

    2022-03-16 05:31:38蔡澤君胡占倉余聯(lián)郴洪唐寶朱呈祥尤延鋮
    空氣動力學學報 2022年1期
    關鍵詞:型面進氣道馬赫數(shù)

    蔡澤君,胡占倉,余聯(lián)郴,洪唐寶,朱呈祥,尤延鋮

    (廈門大學 航空航天學院,廈門 361000)

    0 引 言

    隨著吸氣式推進技術的深入發(fā)展,基于國防以及經(jīng)濟等各方面的戰(zhàn)略需求,可重復使用寬速域飛行器成為各國的必爭高地[1-3]。從推進系統(tǒng)角度出發(fā),單一模態(tài)的動力系統(tǒng)已不再滿足寬速域飛行需求,因此近些年國內(nèi)外多種組合動力形式應運而生,如渦輪基組合循環(huán)發(fā)動機(Turbine Based Combined Cycle engine,TBCC)[4-6]、火 箭 基 組 合 循 環(huán) 發(fā) 動 機(Rocket Based Combined Cycle engine,RBCC)[7],以及火箭增強的組合循環(huán)發(fā)動機:TRRE(Turbo-aided Rocket-augmented Ramjet combined cycle Engine)[8-10]、預冷空氣渦輪火箭PATR(Pre-cooling Air Turbo Rocket)[11]、SERJ(Super Charged Ejector Ramjet)[12]等。

    組合循環(huán)發(fā)動機的發(fā)展不斷推動著進氣道等部件的設計技術進步,而進氣道作為吸氣式推進系統(tǒng)的關鍵部件之一,其性能也決定了推進系統(tǒng)的工作狀態(tài)。為了匹配組合循環(huán)發(fā)動機寬速域的工作需求,進氣道也必須從單一模態(tài)轉(zhuǎn)向多模態(tài)工作,譬如“SR-71”采用的通過移動激波錐匹配不同工況的串聯(lián)型組合進氣道[13]。針對馬赫數(shù)4及以上的寬域工作需求,目前國際上已經(jīng)形成了以二元為主的各類組合進氣道方案,典型的包括NASA的X43-B進氣道[14]、洛馬的外并聯(lián)型組合進氣道[15-16]、南京航空航天大學的內(nèi)并聯(lián)型組合進氣道[17]等。

    由于強壓縮、高流量捕獲、迎風面積小等優(yōu)勢,三維內(nèi)收縮組合進氣道近年來正在被高度關注,新近披露的SR-72[18]、HAWC[19]、波音高超聲速飛機[20]等均不約而同采用了三維內(nèi)收縮組合進氣道方案。事實上,Kothari等在20世紀90年代開發(fā)高超聲速飛行器設計代碼時,就發(fā)現(xiàn)內(nèi)收縮進氣道在寬域內(nèi)表現(xiàn)出更高的性能[21]。美國Aerojet公司在“Trijet”三通道組合動力系統(tǒng)的研究中,專門對其三維內(nèi)收縮進氣道開展了大量風洞試驗,獲得了可觀的該類進氣道的穩(wěn)態(tài)/動態(tài)性能數(shù)據(jù)[22]。閔浩等也根據(jù)“Trijet”的理念設計了一種三通道內(nèi)并聯(lián)可調(diào)內(nèi)收縮進氣道,并對其通道間的干擾特性展開了研究[23]。田方超等則研究了內(nèi)收縮進氣道進出口形狀對其氣動特性的影響,并通過控制泄流槽位置以及泄流量改善了內(nèi)收縮進氣道的低馬赫數(shù)起動問題[24]。張航等從工程角度研究了不同前緣鈍化半徑對內(nèi)轉(zhuǎn)進氣道性能特性的影響,發(fā)現(xiàn)鈍化半徑在3 mm以內(nèi)時,對進氣道性能及流場結構影響較小[25]。鄔鳳林等針對進氣道的調(diào)節(jié)機構展開了研究,通過調(diào)節(jié)二元斜板,拓寬了單通道可調(diào)內(nèi)轉(zhuǎn)進氣道的工作范圍[26]。朱呈祥等基于Xiamen Turbine Ejector Ramjet(XTER)組合發(fā)動機的總體需求,設計了工作速域馬赫數(shù)0~6的內(nèi)收縮多通道組合進氣道,并針對低速通道開展了風洞試驗,得到了渦輪通道的出口氣流特性[27]。胡占倉等[28]進一步針對XTER進氣道開展了模態(tài)轉(zhuǎn)換過程的特性分析,發(fā)現(xiàn)在區(qū)間馬赫數(shù)完成模態(tài)轉(zhuǎn)換將具有更優(yōu)的氣動性能。

    本文將詳細梳理并介紹XTER內(nèi)收縮組合進氣道的設計理念,重點分析XTER內(nèi)收縮組合進氣道的復雜多通道流動結構、全速域流量特性及反壓規(guī)律,旨在進一步探究寬速域內(nèi)收縮組合進氣道在模態(tài)轉(zhuǎn)換過程中的性能變化規(guī)律,為組合進氣道的穩(wěn)定工作及協(xié)同控制提供依據(jù)。

    1 XTER內(nèi)收縮組合進氣道設計理念

    1.1 內(nèi)收縮組合進氣道設計要素分析

    寬速域組合進氣道的設計難點在于兼顧進氣道的高、低馬赫數(shù)性能,并同時在受限空間約束下滿足不同發(fā)動機的流量需求。因此,組合進氣道設計過程應當充分考慮全速域的綜合性能。圖1給出了XTER內(nèi)收縮組合進氣道的設計理念,其中橢圓框均為進氣道設計的輸入條件,圓角矩形框均為設計約束。

    內(nèi)收縮組合進氣道的設計輸入為基本流場、入口形狀,內(nèi)收縮組合進氣道的約束條件為幾何尺寸約束、全速域性能需求以及驅(qū)動功率限制。由圖1可以看出,內(nèi)收縮組合進氣道是否滿足約束條件與進氣道的設計輸入條件直接相關,但輸入條件與約束條件之間存在復雜的信息傳遞,因此難以直接通過入口形狀及基本流場判斷內(nèi)收縮組合進氣道是否能夠滿足各類約束要求。而一旦在進氣道設計的中間過程發(fā)現(xiàn)輸入?yún)?shù)無法滿足約束要求,則需修改相應輸入?yún)?shù)并進行進氣道的重新設計,這將極大延長設計周期。因此,首先有必要梳理各設計要素對內(nèi)收縮組合進氣道尺寸約束、性能需求以及驅(qū)動功率限制的影響機制。

    圖1 XTER內(nèi)收縮組合進氣道的設計理念圖Fig. 1 Design concept of XTER inlet

    內(nèi)收縮組合進氣道的尺寸約束主要包括長度約束及周向輪廓約束。圖2(a)為基本流場及入口形狀對內(nèi)收縮組合進氣道幾何特性的影響機制,圖2(b)為入口形狀在基本流場中進行流線追蹤得到進氣道壓縮型面的示意圖。進氣道壓縮型面長度主要受入口面積及基本流場的影響。基本流場的設計參數(shù)包括設計點條件、流場構型及中心體高度。其中,設計點條件(包含來流馬赫數(shù)、來流靜壓、來流靜溫)決定了進氣道的入射激波角度,入射激波與中心體的交點則決定了基本流場的唇口點,唇口點及反射激波的角度最終決定了基本流場的出口高度及位置。由此,基本流場的無量綱長度得以確定,再根據(jù)入口面積及捕獲流量需求確定縮放因子,即可計算出進氣道的壓縮型面長度。

    內(nèi)收縮組合進氣道周向輪廓約束主要對入口形狀進行了限制。內(nèi)收縮組合進氣道入口形狀僅決定組合進氣道的最小周向輪廓,組合進氣道的最大周向輪廓需結合調(diào)節(jié)機構引起的型面變化考慮。

    入口形狀影響了調(diào)節(jié)機構的周向形狀設計,基本流場影響了調(diào)節(jié)機構的流向形狀設計,換言之,基本流場也受到了周向輪廓約束的間接限制。若將傳統(tǒng)的ICFC、ICFD等基本流場應用于組合進氣道的設計,在型面上任意位置設計調(diào)節(jié)機構后,即使小幅偏轉(zhuǎn)也將出現(xiàn)明顯的氣流膨脹角,并且彎曲壁面也增加了大幅偏轉(zhuǎn)的實現(xiàn)難度。而雙波入射基本流場則表現(xiàn)出明顯的優(yōu)勢,雙波入射基本流場是采用兩道入射激波和一道反射激波構成壓縮波系的軸對稱基本流場,其構型如圖2(b)所示。該構型的二級壓縮型面適合作為繞二級壓縮角處轉(zhuǎn)軸旋轉(zhuǎn)的調(diào)節(jié)機構,其型面僅在二級壓縮型面偏轉(zhuǎn)量大于二級壓縮角度時才會導致氣流膨脹,因此雙波入射基本流場既具有內(nèi)收縮進氣道高壓縮效率的優(yōu)點,又具有滿足組合進氣道型面調(diào)節(jié)需要的優(yōu)秀可調(diào)能力。內(nèi)收縮組合進氣道的高可調(diào)性對入口形狀提出的特殊要求主要體現(xiàn)在調(diào)節(jié)機構對流量的分配上,即組合進氣道的分流面積比變化,該部分內(nèi)容將在下文結合調(diào)節(jié)機構的具體構型展開介紹。

    圖2 幾何約束要素分析圖Fig. 2 Analysis of geometric constraints

    在性能約束方面,本文著重考慮組合進氣道最關鍵的流量分配問題。由于無量綱型面與實際尺寸模型的縮放因子是根據(jù)設計點的流量需求確定的,因此設計點通常不會出現(xiàn)流量不足的問題,以下將著重梳理非設計點的情況。圖3給出了基本流場及入口形狀對通道間流量分配的影響機制。組合進氣道寬速域流量特性與入口形狀的關系主要通過調(diào)節(jié)機構體現(xiàn),而與基本流場的關系則主要體現(xiàn)在非設計點的捕獲流量。在總捕獲流量不變的情況下,調(diào)節(jié)機構的流量分配作用在組合進氣道中顯得格外重要。以下將重點梳理組合進氣道的調(diào)節(jié)機構對流量分配的影響機制。

    圖3 流量特性要素分析圖Fig. 3 Analysis of flow characteristics

    雖然不同構型的分流調(diào)節(jié)機構在調(diào)節(jié)方式、調(diào)節(jié)規(guī)律上存在較大差異,但其影響機制是大同小異的,因此本文以圖4(a)所示的雙分流板構型分流調(diào)節(jié)機構為例進行影響機制分析。該調(diào)節(jié)機構由布置于進氣道上壁面的兩塊旋轉(zhuǎn)分流板構成,其中前分流板繞軸1旋轉(zhuǎn),在控制通道開閉的同時也起著調(diào)節(jié)喉道面積的作用,當上通道完全關閉時,前分流板構成了進氣道的內(nèi)收縮型面。后分流板為前分流板的輔助機構,同樣采用繞軸旋轉(zhuǎn)的運動方式,其主要作用在于使流過前分流板的氣流經(jīng)可控擴張比進入上通道。兩個通道在分流截面前為共用流道,其流通面積隨前分流板的運動而相應改變。

    圖4 分流調(diào)節(jié)機構要素分析圖Fig. 4 Mass distribution mechanism

    通道入口截面積在分流截面的占比是計算通道流量的重要參數(shù)。為建立組合進氣道分流板偏轉(zhuǎn)量與截面積的關系,根據(jù)式(1)定義分流調(diào)節(jié)機構的調(diào)節(jié)效率:

    式中:η為調(diào)節(jié)效率,Amax為分流截面最大面積,其中分流截面為經(jīng)過分流點的豎直截面,Amin為分流截面最小面積,Ai為分流調(diào)節(jié)機構在任意位置時的分流截面積 ,θi為分流調(diào)節(jié)機構從當前位置偏轉(zhuǎn)至最小面積經(jīng)過的角度。顯然,前分流板的調(diào)節(jié)效率與其長度、寬度以及偏轉(zhuǎn)角度有關。但上述參數(shù)均為幾何參數(shù),要確定不同通道的流量還需結合分流截面的氣動參數(shù)。

    分流調(diào)節(jié)機構對進氣道流量特性的影響機制如圖4(b)所示。通過偏轉(zhuǎn)前分流板雖然能夠增大上通道在分流截面的面積占比,從而增加上通道分得的流量,但型面變化導致進氣道收縮比減小的同時也導致共用流道內(nèi)的流場結構發(fā)生變化。因此,隨著前分流板偏轉(zhuǎn)量的增大,分流機構的流量調(diào)節(jié)效率逐漸減小。調(diào)節(jié)效率的減小將使調(diào)節(jié)機構消耗更多的能量。

    調(diào)整分流板的轉(zhuǎn)軸位置,以及分流板的形狀同樣會對調(diào)節(jié)機構的功率需求和能量需求產(chǎn)生影響。將分流板轉(zhuǎn)軸向上游移動,將使分流板的長度增加,能夠使分流板在小幅度的偏轉(zhuǎn)內(nèi)實現(xiàn)更大范圍的流量控制,但會使分流板受到的力矩增大,增加了驅(qū)動機構的負擔。反之,將分流板轉(zhuǎn)軸向下游移動,將降低分流板的流量調(diào)節(jié)效率,但也使分流板所需的驅(qū)動功率減小。換言之,當分流調(diào)節(jié)機構的功率需求高于功率限制時,可通過該方法調(diào)節(jié)分流調(diào)節(jié)機構方案。類似的,通過改變分流板的形狀,增加分流板的寬度將提高調(diào)節(jié)機構的調(diào)節(jié)效率,但會增加分流板彎扭的風險。

    在完成調(diào)節(jié)機構的設計后,隨之而來的問題是組合進氣道的通道設計以及調(diào)節(jié)機構的調(diào)節(jié)規(guī)律。根據(jù)上文可知,內(nèi)收縮組合進氣道通常以高馬赫數(shù)為設計點,因此其高速通道設計與傳統(tǒng)的單通道進氣道類似。以下將重點分析組合進氣道的調(diào)節(jié)規(guī)律及低速通道的設計。

    首先需要確定的是組合進氣道不同非設計點的分流截面積占比,從而確定調(diào)節(jié)機構的控制規(guī)律。由于調(diào)節(jié)機構引起的型面變化既改變了各通道分流截面積的占比,又改變了組合進氣道的波系結構,因此無法簡單根據(jù)流量需求確定機構的調(diào)節(jié)量,但可以根據(jù)改進的MFQE方法[29]進行快速評估。該方法可以通過流線確定調(diào)節(jié)量與流量分配的定量關系,從而明確組合進氣道的調(diào)節(jié)規(guī)律。

    其次是低速通道的型面設計。在確定的調(diào)節(jié)規(guī)律下,低速通道的入口最大截面形狀也是確定的,低速通道是由該截面形狀光順過渡至出口截面的變截面通道。為了保證組合進氣道的模態(tài)轉(zhuǎn)換性能,低速通道的擴張比必須可調(diào)且可控,如圖4所示,通過后分流板隨前分流板的同步偏轉(zhuǎn),通道內(nèi)的擴張比變化連續(xù)且平緩。

    基于以上分析,內(nèi)收縮組合進氣道要實現(xiàn)兼顧設計點及非設計點的各類約束要求,需對各設計參數(shù)進行全面的統(tǒng)籌,其中入口形狀、基本流場、分流調(diào)節(jié)機構尤為重要。入口形狀、基本流場作為設計輸入,在設計過程中演化出進氣道的型面,并使分流調(diào)節(jié)機構的設計受到限制。不同構型的基本流場具有不同的型面調(diào)節(jié)能力,雙波入射基本流場是一種具有良好可調(diào)性并且適用于組合進氣道設計的基本流場。入口形狀的設計將影響型面調(diào)節(jié)對應的組合進氣道分流截面積變化量,從而影響調(diào)節(jié)機構的流量分配效果。分流調(diào)節(jié)機構是進氣道設計過程中的一個重要節(jié)點,需要同時滿足幾何要求、性能要求、驅(qū)動功率限制的要求,導致其設計難度較大。分流調(diào)節(jié)機構設計難度大的另一個原因是分流調(diào)節(jié)機構是非設計點型面的設計的直接輸入,對非設計點的型面設計影響極大。分流調(diào)節(jié)機構的構型在設計點型面設計過程中已經(jīng)確定,但其工作特性是否滿足要求則需非設計點型面設計完成后才能全部確定。這不僅使分流調(diào)節(jié)機構從設計到工作特性確定的中間過程變得更為復雜,也使內(nèi)收縮組合進氣道設計參數(shù)間的相互制約、相互影響問題更為突出。此外,調(diào)節(jié)機構也是低速通道設計的直接輸入,據(jù)此設計擴張比可控可調(diào)的低速通道能夠提高模態(tài)轉(zhuǎn)換過程中的氣流品質(zhì)。

    1.2 XTER內(nèi)收縮組合進氣道方案

    根據(jù)XTER組合動力系統(tǒng)的幾何及性能要求,其內(nèi)收縮組合進氣道需滿足兩側(cè)雙渦輪發(fā)動機、上側(cè)引射亞燃發(fā)動機及下側(cè)超燃發(fā)動機的工作需求?;诘?.1節(jié)設計要素分析,首先采用具有兩道入射激波的雙波入射內(nèi)收縮基本流場開展馬赫數(shù)0~6的進氣道氣動型面設計。為保證分流調(diào)節(jié)機構的可實現(xiàn)性,進氣道入口形狀設計為類矩形;為了兼顧分流調(diào)節(jié)機構的流量控制能力、提高擴張通道的抗反壓能力,采用雙旋轉(zhuǎn)分流板的分流調(diào)節(jié)機構。通過流線追蹤、附面層修正等技術[24],生成高速通道壓縮型面,隔離段則設置設計為類矩形轉(zhuǎn)橢圓形的微擴通道。最終得到的XTER內(nèi)收縮組合進氣道模型如圖5所示。

    圖5 三維內(nèi)轉(zhuǎn)組合進氣道氣動模型示意圖Fig. 5 Schematic diagram of inward-turning combined cycle inlet

    XTER內(nèi)收縮組合進氣道采用雙波入射基本流場進行壓縮型面設計,根據(jù)通道出口的畸變及性能要求,對兩側(cè)渦輪通道及引射亞燃通道采用基于中心線方程、面積變化規(guī)律及圓角規(guī)律的CST(Class Shape function Transformation)[30]理論方法進行快速構型,同步實現(xiàn)了“S”彎的通道畸變抑制及出口性能提升。XTER內(nèi)收縮進氣道的分流機構為位于渦輪通道兩側(cè)及引射亞燃通道上側(cè)的雙旋轉(zhuǎn)分流板[29],兩次模態(tài)轉(zhuǎn)換的分流板旋轉(zhuǎn)方式如圖6所示。到目前為止,已對XTER內(nèi)收縮組合進氣道的全速域工況進行了仿真,并對其低速特性開展了試驗驗證。下文將著重介紹XTER內(nèi)收縮組合進氣道與下游發(fā)動機的流量及反壓匹配特性。

    圖6 模態(tài)轉(zhuǎn)換過程示意圖Fig. 6 Schematic diagram of mode transition

    2 數(shù)值方法及驗證

    2.1 數(shù)值方法及邊界條件

    本文所涉及的數(shù)值仿真均采用ANSYS Fluent軟件進行模擬,求解基于密度的三維雷諾平均Navier-Stokes方程,湍流模型選用k-ωSST模型,無黏對流通量采用二階迎風Roe-FDS差分格式求解。計算過程以連續(xù)方程、動量方程、能量方程及k-ω方程殘差下降至少三個數(shù)量級且進氣道出口截面流量穩(wěn)定作為收斂判據(jù)。假定空氣為理想氣體,定壓比熱cp選用Piecewise-Polynomial擬合,黏度采用Sutherland公式求解。定義來流迎風面為壓力遠場條件,各通道出口截面為壓力出口條件,壁面均不考慮傳熱采用絕熱無滑移邊界條件。

    為了滿足不同馬赫數(shù)的計算域尺度需求且兼顧組合進氣道的通道調(diào)節(jié)工況,本文共涉及11套計算網(wǎng)格。其中,渦輪模態(tài)、引射亞燃模態(tài)各2套,超燃模態(tài)1套,兩個模態(tài)轉(zhuǎn)換過程各3套。每套網(wǎng)格均對近壁面進行加密,且保證y+不大于30。表1匯總了本文涉及的計算工況的來流條件、工作模態(tài)及網(wǎng)格量。表中Ma為來流馬赫數(shù),H為飛行高度,Mesh為網(wǎng)格量,Mode為工作模態(tài),其中TUR為渦輪模態(tài),ERJ為引射亞燃模態(tài),SCR為超燃模態(tài),MT-TUR-ERJ為渦輪-引射亞燃模態(tài)轉(zhuǎn)換,MT-ERJ-SCR為引射亞燃-超燃模態(tài)轉(zhuǎn)換,百分數(shù)代表模態(tài)轉(zhuǎn)換進程。

    表1 本文計算狀態(tài)匯總表Table 1 Summary of numerical cases

    2.2 數(shù)值驗證

    為了驗證計算結果的可靠性,首先開展了馬赫數(shù)4.5典型工況的網(wǎng)格無關性驗證。通過圖7所示粗網(wǎng)格(200萬)、中等網(wǎng)格(400萬)和細網(wǎng)格(600萬)結果的對比,發(fā)現(xiàn)采用中等網(wǎng)格可以較好地兼顧收斂性與計算效率。其他計算狀態(tài)網(wǎng)格量則參考該工況,并根據(jù)計算域及通道數(shù)進行適度增減調(diào)整。

    圖7 不同網(wǎng)格量計算結果圖Fig. 7 Calculation results of different mesh quality

    采用2.1節(jié)數(shù)值方法對文獻[31]中的試驗模型進行了仿真計算,模型截面如圖8(a)所示,來流條件參考文獻[31]。圖8(b)為試驗與仿真的上下壁面壓力分布結果對比,可以發(fā)現(xiàn),無論是數(shù)值還是曲線變化趨勢,仿真結果均與試驗結果吻合良好。以上驗證表明了所用計算方法可以支撐本文氣動性能的仿真研究。

    圖8 沿程壁面壓力分布的對比Fig. 8 Streamwise variation of wall-pressure

    3 結果分析與討論

    3.1 XTER內(nèi)收縮組合進氣道典型流動結構

    根據(jù)上文介紹的XTER組合進氣道的設計理念及氣動型面,可以發(fā)現(xiàn)該進氣道無論是壓縮型面還是各個通道的型面,都具有突出的三維特征。在馬赫數(shù)0~6的全速域范圍內(nèi),進氣道幾何外形隨工作模態(tài)不斷變化,內(nèi)部三維流動表現(xiàn)出差異性,但最具代表性的要數(shù)各通道全部打開(即渦輪模態(tài))的流動結構。

    如前所述,在渦輪發(fā)動機可能的工作速域范圍(馬赫數(shù)0~2.5),渦輪通道分流板的開合角度取決于總體推力需求及各子動力的工作特性。為了體現(xiàn)XTER內(nèi)收縮進氣道的典型多通道流動結構,選取渦輪模態(tài)最高工作馬赫數(shù)2.5且各通道均打開的通流狀態(tài)進行剖析。圖9為該狀態(tài)下的進氣道流向切片馬赫數(shù)云圖。圖中顯示進氣道入口產(chǎn)生一道入射激波,此時由于來流馬赫數(shù)相對較低,入射激波的激波角較大,激波并未封口。氣流進入壓縮段,由于分流板打開分流,在渦輪、引射亞燃通道分流板轉(zhuǎn)軸處產(chǎn)生膨脹波使氣流加速,另外引射亞燃通道分流段側(cè)面與渦輪通道分流段上壁形成的“臺階”結構,使氣流在分流段受到橫向及縱向的壓縮,因此在分流段形成壓縮波并逐漸匯聚成一道激波,且越靠近下壁面,激波強度越大。這是由于唇口之后為擴張通道,氣流進入唇口加速,靠近下壁面氣流的馬赫數(shù)上升,進而提高了激波的強度。經(jīng)過分流,氣流過渡到各個通道,從圖中可以看出,渦輪通道上壁面以及下壁面均存在分離區(qū),增加了通道內(nèi)氣流的不均勻性。結合圖10的密度梯度云圖發(fā)現(xiàn),渦輪通道上側(cè)的分離是由于通道入口激波與附面層發(fā)生干擾導致的,而下壁面的流動分離主要是由于S彎型面引起的局部逆壓梯度。此外,在引射亞燃通道出口處可以看到上壁面的附面層較厚,這是因為引射亞燃通道入口上壁面的局部分離逐漸向下游發(fā)展所致,引射亞燃通道與超燃通道型面相對平緩過渡,通道內(nèi)并無明顯分離,且通道出口氣流相對均勻。

    圖9 馬赫數(shù)2.5進氣道流向切片馬赫數(shù)云圖Fig. 9 Mach number distributions of streamwise slices in the inlet at Mach 2.5

    在圖10所示密度梯度云圖中可以清晰地分辨通道內(nèi)波系結構的形狀,包括入射激波、唇口處的反射激波及分流板轉(zhuǎn)軸處的膨脹波等。同時,在對稱面以及渦輪通道切片圖中,可以清晰地看到上文所述分流段內(nèi)壓縮波匯聚形成的激波,該道激波呈現(xiàn)明顯的三維彎曲特征,并逐漸發(fā)展交匯于超燃通道下壁面,與下壁面附面層干擾導致局部分離,并誘導產(chǎn)生分離激波。此外,各個通道的入口均存在一道激波,是由于在分流面位置,氣流要過渡到各個通道,通過該道激波使超聲速氣流的流動方向改變。在引射亞燃通道的上側(cè)壁面,由于通道入口的誘導激波與上側(cè)壁面附面層干擾,導致了局部的小分離,使附面層變厚并逐漸發(fā)展到下游出口。

    圖10 馬赫數(shù)2.5進氣道密度梯度的切片圖Fig. 10 Density gradient in the inlet at Mach 2.5

    通過以上分析,XTER內(nèi)收縮組合進氣道內(nèi)部存在多波系結構,通道內(nèi)有明顯的激波相交以及激波-邊界層干擾現(xiàn)象,并出現(xiàn)不同尺寸的分離,渦輪模態(tài)下S彎通道內(nèi)的上述特征尤為明顯。但需要指出的是,在實際XTER組合動力工作中,渦輪通道受出口壓力影響會產(chǎn)生結尾激波,且結尾激波的位置波動會帶來通道內(nèi)流動結構的變化,并促使S彎通道展現(xiàn)出氣流摻混降畸變的優(yōu)勢。只要通道出口壓力不高于許用壓力,結尾激波便不會對壓縮段及通道入口處的波系結構產(chǎn)生影響。

    當XTER組合發(fā)動機工作于引射亞燃模態(tài)時,內(nèi)收縮組合進氣道僅剩上下兩個通道。根據(jù)前文介紹的XTER工作模式,引射亞燃模態(tài)下通道的分流板也會根據(jù)總體工作需求進行適度調(diào)節(jié)。圖11給出了引射亞燃模態(tài)、馬赫數(shù)4、通流,且引射亞燃通道分流板旋轉(zhuǎn)至與上游壓縮面光順過渡狀態(tài)的流向馬赫數(shù)切片云圖,此時可以明顯觀察到入射激波的三維彎曲特征。與馬赫數(shù)2.5狀態(tài)下的入射激波相比,此時入射激波距離唇口更近,進氣道具有較高的流量捕獲能力。經(jīng)過壓縮段,氣流經(jīng)過分流面過渡到兩個通道,在引射亞燃通道的上側(cè)壁面產(chǎn)生局部分離。結合圖12對稱面內(nèi)的馬赫數(shù)及密度梯度云圖,發(fā)現(xiàn)這時是由于引射亞燃通道入口的激波與后分流板壁面附面層干擾較強,局部邊界層顯著增厚,且在角區(qū)逐漸卷起低能低速橫向渦向通道出口上壁面發(fā)展。同樣在超燃通道內(nèi),由于型面從類矩形過渡到橢圓形出口,在角區(qū)過渡位置形成局部渦并向下游發(fā)展,在圖中超燃通道出口上側(cè)壁面可以看到。

    圖11 馬赫數(shù)4進氣道流向切片馬赫數(shù)云圖Fig. 11 Mach number contours at several streamwise slices in the inlet at Mach 4

    圖12 馬赫數(shù)4對稱面馬赫數(shù)及密度梯度云圖Fig. 12 Mach number and density gradient on symmetry plane at Mach 4

    圖12為該狀態(tài)下進氣道對稱面馬赫數(shù)和密度梯度云圖。此時入射激波及唇口處的反射激波在馬赫數(shù)和密度梯度云圖中均可觀測到。由于此時馬赫數(shù)較高,黏性效應增強,壓縮型面上的附面層沿流向明顯變厚,采用抽吸方式進行近壁流動控制,保證喉道處的氣動性能。在密度梯度云圖中,唇口上方存在一道壓縮波,這是由于側(cè)面的二級壓縮型面向內(nèi)收縮產(chǎn)生的。在喉道分流位置,兩側(cè)通道入口均存在一道誘導激波使氣流過渡到兩個通道,引射亞燃通道入口處的激波經(jīng)反射相交于通道下壁面,由于通道型面擴張,激波強度也逐漸增強。此外,由于引射亞燃通道后分流板的旋轉(zhuǎn),其與通道間存在明顯的轉(zhuǎn)角過渡,因此產(chǎn)生了較強且長的彎曲激波結構,該激波與上述反射激波均交匯于引射亞燃通道下壁面,誘導附面層發(fā)生了大尺度分離。在圖12的馬赫數(shù)云圖中可以明顯觀察到引射亞燃通道下壁面的低速低能區(qū),而從密度梯度云圖可以清晰地發(fā)現(xiàn)激波/附面層干擾及誘導分離。

    以上對XTER內(nèi)收縮組合進氣道渦輪模態(tài)與引射亞燃模態(tài)的兩個典型狀態(tài)進行了流動結構剖析。在單通道的超燃模態(tài)下,進氣道接近設計狀態(tài),入射激波幾乎封口,反射激波在隔離段內(nèi)反射,波系相對常規(guī),在此不再贅述。

    模態(tài)轉(zhuǎn)換對組合進氣道的全速域工作至關重要,XTER內(nèi)收縮組合進氣道在工作過程中將完成渦輪-引射亞燃模態(tài)轉(zhuǎn)換及引射亞燃-超燃模態(tài)轉(zhuǎn)換。由于渦輪-引射亞燃模態(tài)轉(zhuǎn)換過程中涉及的通道數(shù)多,流動也更具代表性,因此下面重點分析該模態(tài)轉(zhuǎn)換過程中的流動特征。

    圖13為馬赫數(shù)2.5渦輪-引射亞燃模態(tài)轉(zhuǎn)換的密度梯度云圖。百分數(shù)代表模態(tài)轉(zhuǎn)換的進程,渦輪模態(tài)為0%,引射亞燃模態(tài)為100%。隨著模態(tài)轉(zhuǎn)換的進行,入射激波未發(fā)生變化,前分流板轉(zhuǎn)軸處由分流板偏轉(zhuǎn)產(chǎn)生的膨脹波逐漸減弱。由于分流板的偏轉(zhuǎn),進氣道的內(nèi)收縮比逐漸增大,渦輪通道入口截面馬赫數(shù)降低,入射激波角度逐漸增大,并在后分流板上壁面形成誘導分離。在模態(tài)轉(zhuǎn)換75%進程時,該激波未與后分流板相交,而是由兩分流板間的泄流口排出,但由于渦輪通道此時的擴張比較大,氣流在通道內(nèi)膨脹加速,在后分流板轉(zhuǎn)軸處發(fā)生偏轉(zhuǎn),形成強度較大的激波。

    圖13(b)為組合進氣道唇口點處的切片密度梯度云圖。前文所述角區(qū)匯聚形成的三維激波隨著模態(tài)轉(zhuǎn)換過程的進行,激波角度逐漸增大,因此切片圖中激波在該截面的截曲線逐漸向?qū)ΨQ面靠近。可以看出75%進程時的激波截曲線與0%~50%進程的方向相反,該截曲線是截取了由內(nèi)收縮進氣道左右半模的角區(qū)激波在該截面的上游相交并形成的透射波。結合圖13(a)可以看出,在0%~50%進程中,該透射波進入超燃沖壓通道,未對渦輪通道產(chǎn)生影響,而75%進程時,該透射波進入渦輪通道。

    圖13 渦輪-引射亞燃模態(tài)轉(zhuǎn)換的密度梯度云圖Fig. 13 Density gradient at lip section during the mode transition of turbine to ejector

    本節(jié)對兩個典型工況通流狀態(tài)下的進氣道流動結構以及馬赫數(shù)2.5渦輪-引射亞燃模態(tài)轉(zhuǎn)換過程流動特性進行了介紹,相關剖析有利于進一步揭示XTER內(nèi)收縮組合進氣道的設計理念,也為后續(xù)型面優(yōu)化指明了方向。在對XTER內(nèi)收縮組合進氣道流動結構形成清晰認識的基礎上,下文將重點介紹該進氣道在全速域范圍內(nèi),尤其是模態(tài)轉(zhuǎn)換過程中的性能變化規(guī)律。

    3.2 全速域特性

    XTER內(nèi)收縮組合進氣道涉及的流道多,工作狀態(tài)復雜,需結合發(fā)動機的總體需求及各子動力工作特性確定實際工況。但就進氣道部件本身而言,必須對其全速域的各種可能工況進行分析,以便為XTER發(fā)動機的集成工作提供更為全面的數(shù)據(jù)支撐。

    圖14給出了XTER進氣道在馬赫數(shù)0~6全速域范圍內(nèi)的流量系數(shù)及總壓恢復特性。這里涉及的渦輪模態(tài)均為四個通道全開的通流狀態(tài);馬赫數(shù)3時引射亞燃通道全開,馬赫數(shù)4時引射亞燃通道半開;超燃模態(tài)引射亞燃通道開小角度泄流。圖14中柱狀圖代表XTER進氣道的總流量系數(shù),空心點為各通道出口的流量系數(shù)φ,而實心點為各通道出口的總壓恢復系數(shù)σ。由圖可以發(fā)現(xiàn),在整個飛行包線范圍內(nèi),XTER內(nèi)收縮組合進氣道的總流量系數(shù)變化比較緩和,數(shù)值一直在0.75以上,表明該進氣道優(yōu)秀的捕獲能力。同時,隨著來流馬赫數(shù)的增加,進氣道的總流量系數(shù)逐步增大至1,能夠確保XTER組合進氣道具有好的高速特性,且全馬赫數(shù)范圍內(nèi)具備較小的溢流阻力。另一反面,隨著來流馬赫數(shù)的增加,進氣道內(nèi)的流動損失也逐漸升高,例如馬赫數(shù)0.5時各通道出口接近1的總壓恢復系數(shù)在馬赫數(shù)6超燃通道出口降至不到0.4。

    圖14 XTER進氣道全速域工作特性圖Fig. 14 Operating characteristics of XTER inlet in full speed range

    對于組合進氣道而言,模態(tài)轉(zhuǎn)換過程中的特性規(guī)律十分關鍵。為此,本文針對該多通道進氣道在渦輪-引射亞燃模態(tài)轉(zhuǎn)換及引射亞燃-超燃模態(tài)轉(zhuǎn)換過程中的性能規(guī)律進行了研究。

    考慮到XTER發(fā)動機的工作模態(tài),分析了馬赫數(shù)2及馬赫數(shù)2.5時進行渦輪-引射亞燃模態(tài)轉(zhuǎn)換的兩種情況。圖15為該進氣道在馬赫數(shù)2進行渦輪-引射亞燃模態(tài)轉(zhuǎn)換時的工作特性圖。其中,橫坐標為分流板關閉渦輪通道面積的百分比,左側(cè)縱坐標代表實際出口流量與理論捕獲流量之比的流量系數(shù),右側(cè)縱坐標代表出口總壓與來流總壓之比的總壓恢復系數(shù)。圖中實心點均表示流量系數(shù)變化,空心點均表示總壓恢復系數(shù)變化;紅色、藍色與綠色分別代表引射亞燃通道參數(shù)、超燃通道參數(shù)和渦輪通道參數(shù);黑色實線為進氣道的總流量系數(shù)??梢园l(fā)現(xiàn),在模態(tài)轉(zhuǎn)換的前50%進程中,進氣道總流量系數(shù)受分流板泄流量增大的影響由0.75小幅下降至0.74。具體到各通道而言,隨著渦輪通道分流板的關閉,渦輪通道流量系數(shù)由0.32下降至0.18;引射亞燃通道流量系數(shù)由0.26增加至0.32;超燃通道流量系數(shù)由0.16增加至0.23。相應地,渦輪通道總壓恢復系數(shù)由0.63下降至0.38;引射亞燃通道總壓恢復系數(shù)由0.84增加至0.87;超燃通道總壓恢復系數(shù)由0.78增至0.92。但是當模態(tài)轉(zhuǎn)換進一步進行至75%時,可以發(fā)現(xiàn)進氣道總流量系數(shù)出現(xiàn)了0.15的急劇下降,且各通道總壓損失均在增加,這表明此時進氣道已經(jīng)進入流量擁塞的不起動狀態(tài)。換言之,該組合進氣道在馬赫數(shù)2進行渦輪-引射亞燃模態(tài)轉(zhuǎn)換時將面臨困難,在設計XTER發(fā)動機控制律時需重點關注。

    圖15 馬赫數(shù)2渦輪-引射亞燃模態(tài)轉(zhuǎn)換特性Fig. 15 Mode transition from turbine to ejector at Mach 2

    鑒于馬赫數(shù)2要進行渦輪-引射亞燃的完整模態(tài)轉(zhuǎn)換存在困難,因此進一步考慮馬赫數(shù)2.5的可行性。圖16為馬赫數(shù)2.5時渦輪-引射亞燃模態(tài)轉(zhuǎn)換過程中的工作特性。隨著渦輪通道的逐漸關閉,進氣道壓縮激波角度略有增加,進氣道總流量系數(shù)穩(wěn)中有降。具體到各通道出口參數(shù)可以發(fā)現(xiàn),隨著渦輪通道流通面積逐漸減小,渦輪通道流量系數(shù)顯著降低,由0.32下降至0;引射亞燃通道及超燃通道流量系數(shù)均線性增加且增長幅度相當??倝夯謴头矫妫瑴u輪通道逐漸關閉時總壓恢復系數(shù)相應降低,結合圖13可以看出,在模態(tài)轉(zhuǎn)換50%進程時,后分流板與通道入口氣流形成明顯的膨脹角度,通道入口激波在上壁面形成分離并沿膨脹型面迅速向下游發(fā)展,導致渦輪通道出口總壓恢復性能顯著下降,由0.58下降至0.38,在關閉75%時渦輪通道內(nèi)已出現(xiàn)大范圍低能區(qū),總壓恢復系數(shù)僅為0.16,此時渦輪發(fā)動機的推力貢獻將十分有限。另一方面,渦輪通道的逐漸關閉使得引射亞燃通道及超燃通道總壓恢復系數(shù)均穩(wěn)步提升,至模態(tài)轉(zhuǎn)換結束時,引射亞燃通道與超燃通道的出口總壓恢復系數(shù)均接近或超過0.85。當然,以上分析僅針對XTER工作模態(tài)要求中的兩個典型馬赫數(shù)進行了渦輪-引射亞燃模態(tài)轉(zhuǎn)換特性的分析,實際進氣道調(diào)節(jié)規(guī)律及馬赫數(shù)轉(zhuǎn)換區(qū)間還需根據(jù)XTER發(fā)動機的總體規(guī)劃進行優(yōu)選。

    圖16 馬赫數(shù)2.5渦輪-引射亞燃模態(tài)轉(zhuǎn)換特性Fig. 16 Mode transition from turbine to ejector at Mach 2.5

    針對引射亞燃-超燃的模態(tài)轉(zhuǎn)換,考慮了另一種寬馬赫數(shù)區(qū)間轉(zhuǎn)級的可能性。XTER進氣道在馬赫數(shù)3時仍保持引射亞燃通道全開狀態(tài),模態(tài)轉(zhuǎn)換的主要過程在馬赫數(shù)4進行,由0%進行至62.5%,并在62.5%之后逐漸加速至馬赫數(shù)4.5再幾乎完全關閉,如圖17所示。圖中曲線同樣表示引射亞燃通道及超燃通道出口流量及總壓恢復性能的變化??梢园l(fā)現(xiàn),從馬赫數(shù)3至馬赫數(shù)4,來流馬赫數(shù)的增加帶來入射激波角度的減小,因此進氣道的總流量系數(shù)由0.81提升至0.90。在隨后的模態(tài)轉(zhuǎn)換過程中,總流量系數(shù)均維持在0.9以上,這也表明了XTER進氣道的高流量捕獲能力。具體到各通道而言,隨著引射亞燃-超燃模態(tài)轉(zhuǎn)換過程的進行,引射亞燃通道及超燃通道的流量分配發(fā)生了顯著變化,引射亞燃通道流量系數(shù)由0.48降至0,超燃通道流量則相應增加。在總壓恢復方面,引射亞燃通道在整個模態(tài)轉(zhuǎn)換過程中總壓恢復系數(shù)是逐漸下降的,而且在前62.5%過程中總壓恢復由0.81下降至0.58,在后37.5%的過程中才完全降至0。這也表明引射亞燃通道在前半程的模態(tài)轉(zhuǎn)換過程中仍然具有較強的工作能力,XTER發(fā)動機應結合通道特性綜合考慮其推力分配。超燃通道在引射亞燃-超燃模態(tài)轉(zhuǎn)換的全過程中總壓恢復系數(shù)由0.83降至0.64,該性能并不特別理想,后續(xù)可考慮在馬赫數(shù)4.5時進一步排移壓縮面的近壁低能流。

    圖17 引射亞燃-超燃模態(tài)轉(zhuǎn)換特性Fig. 17 Mode transition from ejector to scramjet

    基于以上分析,XTER內(nèi)收縮組合進氣道可以實現(xiàn)寬速域范圍內(nèi)的連續(xù)正常工作,具體的模態(tài)轉(zhuǎn)換規(guī)律應結合XTER發(fā)動機的總體規(guī)劃與部件匹配進行綜合考量。

    4 總 結

    本文對XTER組合動力系統(tǒng)的三維內(nèi)收縮進氣道進行了梳理與分析,形成以下主要結論:

    1)XTER內(nèi)收縮組合進氣道眾多設計要素的耦合關系增大了設計難度,尤其是分流調(diào)節(jié)機構的設計。分流調(diào)節(jié)機構受到幾何尺寸、性能需求、驅(qū)動功率等約束,同時其作為非設計點型面設計的直接輸入,還需滿足非設計點的工作需要。因此,分流調(diào)節(jié)機構設計難度大,也使設計參數(shù)間的相互制約問題更為突出。

    2)XTER內(nèi)收縮進氣道的內(nèi)部流動具有明顯的三維特征。在馬赫數(shù)2.5的渦輪模態(tài)下,三維特征體現(xiàn)在分流段內(nèi)壓縮波匯聚形成的三維彎曲激波。該彎曲激波在馬赫數(shù)2.5渦輪-引射亞燃模態(tài)轉(zhuǎn)換過程中激波角逐漸增大,透射波在0%~50%進程中進入超燃沖壓通道,在75%進程時進入渦輪通道。在馬赫數(shù)4的引射亞燃模態(tài)下,三維特征體現(xiàn)在入射激波的彎曲,以及對稱面唇口上方出現(xiàn)一道強度較大的壓縮波,該壓縮波是由二級壓縮型面的側(cè)壁面的內(nèi)收縮引起的。

    3)XTER內(nèi)收縮組合進氣道的總流量系數(shù)在整個飛行包線范圍內(nèi)變化比較緩和,數(shù)值一直在0.75以上。在模態(tài)轉(zhuǎn)換過程中,通過控制分流調(diào)節(jié)機構可實現(xiàn)流量的平穩(wěn)過渡,并保證了通道出口的性能參數(shù)。

    4)馬赫數(shù)2.5渦輪-引射亞燃模態(tài)轉(zhuǎn)換過程中,總流量系數(shù)維持在0.75以上,且在模態(tài)轉(zhuǎn)換完成時引射亞燃通道與超燃通道的出口總壓恢復系數(shù)均接近或超過0.85。在馬赫數(shù)3~4.5進行的引射亞燃-超燃模態(tài)轉(zhuǎn)換過程中,引射亞燃通道在前62.5%總壓恢復系數(shù)下降有限,由0.81下降至0.58,表明該通道在前半程的模態(tài)轉(zhuǎn)換過程中仍具有較強的工作能力。

    致謝:感謝1912項目組對本文工作的支持。

    猜你喜歡
    型面進氣道馬赫數(shù)
    一維非等熵可壓縮微極流體的低馬赫數(shù)極限
    基于AVL-Fire的某1.5L發(fā)動機進氣道優(yōu)化設計
    基于輔助進氣門的進氣道/發(fā)動機一體化控制
    載荷分布對可控擴散葉型性能的影響
    基于數(shù)值分析的汽車A柱加強板模具型面改進
    模具型面數(shù)控加工自動化編程系統(tǒng)開發(fā)
    基于鋁擠壓模具倒扣型面的三軸加工應用
    The coupling characteristics of supersonic dual inlets for missile①
    某柴油機進氣道數(shù)值模擬及試驗研究
    汽車零部件(2014年2期)2014-03-11 17:46:30
    NF-6連續(xù)式跨聲速風洞馬赫數(shù)控制方式比較與研究
    久久久a久久爽久久v久久| 91aial.com中文字幕在线观看| 又粗又硬又长又爽又黄的视频| 亚洲国产最新在线播放| 亚洲人成网站在线播| 少妇的逼水好多| 欧美日韩一区二区视频在线观看视频在线 | 青春草亚洲视频在线观看| 男人的好看免费观看在线视频| 国产单亲对白刺激| 久久久久免费精品人妻一区二区| 国产爱豆传媒在线观看| 国产成人福利小说| 国产黄色视频一区二区在线观看 | 青春草亚洲视频在线观看| 成人av在线播放网站| 国产色爽女视频免费观看| 国产高清不卡午夜福利| 少妇被粗大猛烈的视频| 少妇猛男粗大的猛烈进出视频 | 少妇熟女欧美另类| 亚洲性久久影院| 男人舔奶头视频| 亚洲一级一片aⅴ在线观看| 91精品伊人久久大香线蕉| 女的被弄到高潮叫床怎么办| 亚洲色图av天堂| 精华霜和精华液先用哪个| 深爱激情五月婷婷| 国产v大片淫在线免费观看| 日日摸夜夜添夜夜添av毛片| 精品人妻视频免费看| 欧美性猛交╳xxx乱大交人| 午夜激情福利司机影院| 欧美人与善性xxx| 国产亚洲5aaaaa淫片| a级毛片免费高清观看在线播放| 亚洲国产精品成人久久小说| 国产伦精品一区二区三区四那| 亚洲第一区二区三区不卡| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成人毛片a级毛片在线播放| 国产精品国产三级专区第一集| 日韩精品青青久久久久久| 精品不卡国产一区二区三区| 国产女主播在线喷水免费视频网站 | 亚洲性久久影院| 在线天堂最新版资源| 大话2 男鬼变身卡| 少妇人妻一区二区三区视频| 七月丁香在线播放| 精品久久国产蜜桃| 91在线精品国自产拍蜜月| 久久久色成人| 亚洲国产成人一精品久久久| 国产精品野战在线观看| 汤姆久久久久久久影院中文字幕 | 亚洲精品影视一区二区三区av| 午夜精品国产一区二区电影 | 在线播放国产精品三级| 中文亚洲av片在线观看爽| 免费av毛片视频| 六月丁香七月| 亚洲在线自拍视频| 寂寞人妻少妇视频99o| 五月玫瑰六月丁香| 日韩亚洲欧美综合| 国产激情偷乱视频一区二区| 亚洲成人精品中文字幕电影| 六月丁香七月| 国产淫语在线视频| 亚洲国产精品sss在线观看| 视频中文字幕在线观看| 亚洲图色成人| 人妻夜夜爽99麻豆av| 在线天堂最新版资源| 七月丁香在线播放| 最近2019中文字幕mv第一页| 中文字幕av在线有码专区| 神马国产精品三级电影在线观看| 日本av手机在线免费观看| 嘟嘟电影网在线观看| 日韩,欧美,国产一区二区三区 | 欧美一级a爱片免费观看看| 精品一区二区三区人妻视频| 免费av毛片视频| 午夜福利视频1000在线观看| 少妇的逼水好多| 人人妻人人澡欧美一区二区| 一个人观看的视频www高清免费观看| 亚洲图色成人| 在线观看av片永久免费下载| 国产黄片视频在线免费观看| 国产 一区 欧美 日韩| 日本-黄色视频高清免费观看| 亚洲18禁久久av| 久99久视频精品免费| 小蜜桃在线观看免费完整版高清| 色噜噜av男人的天堂激情| 亚洲精品影视一区二区三区av| 午夜久久久久精精品| 日本免费在线观看一区| 一边亲一边摸免费视频| 国产精品一二三区在线看| 国产精品99久久久久久久久| 免费电影在线观看免费观看| 久久人人爽人人片av| 欧美+日韩+精品| 3wmmmm亚洲av在线观看| 精品欧美国产一区二区三| a级一级毛片免费在线观看| 汤姆久久久久久久影院中文字幕 | 国产亚洲精品久久久com| 啦啦啦韩国在线观看视频| 久久久亚洲精品成人影院| 嫩草影院精品99| 又爽又黄无遮挡网站| 国产精品国产高清国产av| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲经典国产精华液单| 国产黄a三级三级三级人| 欧美xxxx性猛交bbbb| 尾随美女入室| 日韩大片免费观看网站 | 在线a可以看的网站| 亚洲自偷自拍三级| 热99re8久久精品国产| 色5月婷婷丁香| 男人舔女人下体高潮全视频| 精品人妻熟女av久视频| 免费无遮挡裸体视频| 99久久中文字幕三级久久日本| 好男人在线观看高清免费视频| 国产又色又爽无遮挡免| 搡女人真爽免费视频火全软件| 插阴视频在线观看视频| kizo精华| 九九久久精品国产亚洲av麻豆| 免费观看的影片在线观看| 美女被艹到高潮喷水动态| 亚洲精品,欧美精品| 一夜夜www| av线在线观看网站| 久久久亚洲精品成人影院| 观看美女的网站| 高清毛片免费看| 国产成人精品久久久久久| 久久精品综合一区二区三区| 亚洲美女搞黄在线观看| 99视频精品全部免费 在线| 天堂√8在线中文| 久久99热6这里只有精品| h日本视频在线播放| av又黄又爽大尺度在线免费看 | av.在线天堂| 97在线视频观看| 成人二区视频| 日本黄大片高清| 国产精品精品国产色婷婷| 如何舔出高潮| 在线观看av片永久免费下载| av国产久精品久网站免费入址| 精品国内亚洲2022精品成人| 久久亚洲精品不卡| 欧美xxxx性猛交bbbb| 国产毛片a区久久久久| 最近最新中文字幕免费大全7| 简卡轻食公司| 日韩,欧美,国产一区二区三区 | 91狼人影院| 亚洲欧美成人综合另类久久久 | 美女大奶头视频| 亚洲精品影视一区二区三区av| 建设人人有责人人尽责人人享有的 | 日本免费一区二区三区高清不卡| 亚洲欧洲国产日韩| 国产黄色小视频在线观看| a级毛色黄片| 波多野结衣巨乳人妻| 日韩欧美 国产精品| 国产成年人精品一区二区| 亚洲国产最新在线播放| 高清毛片免费看| 国产白丝娇喘喷水9色精品| 欧美日韩精品成人综合77777| 欧美性感艳星| 午夜激情福利司机影院| 免费av观看视频| 亚洲av中文av极速乱| 丰满乱子伦码专区| 一级黄片播放器| 黄色配什么色好看| 一级毛片电影观看 | 自拍偷自拍亚洲精品老妇| 欧美xxxx黑人xx丫x性爽| 国产精品野战在线观看| 国产精品,欧美在线| 午夜激情福利司机影院| 国产一级毛片七仙女欲春2| 在线观看美女被高潮喷水网站| 国产av一区在线观看免费| 婷婷色麻豆天堂久久 | 欧美一区二区精品小视频在线| 91在线精品国自产拍蜜月| 国产色婷婷99| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲国产高清在线一区二区三| 狂野欧美激情性xxxx在线观看| 九九爱精品视频在线观看| 国产高潮美女av| 青春草亚洲视频在线观看| 免费看av在线观看网站| 国内少妇人妻偷人精品xxx网站| 亚洲国产精品成人综合色| 国产精品电影一区二区三区| 九九久久精品国产亚洲av麻豆| 三级毛片av免费| 久久久久久久久久久免费av| 免费大片18禁| 日韩一区二区视频免费看| 亚洲欧洲国产日韩| 一级av片app| 色网站视频免费| 嫩草影院精品99| 中文在线观看免费www的网站| 亚洲精品久久久久久婷婷小说 | 黄色一级大片看看| 嫩草影院精品99| 国产色爽女视频免费观看| 深夜a级毛片| 精品一区二区三区视频在线| 全区人妻精品视频| 一二三四中文在线观看免费高清| 91精品一卡2卡3卡4卡| 99久久中文字幕三级久久日本| 免费看av在线观看网站| 欧美又色又爽又黄视频| 国产亚洲av嫩草精品影院| 国产探花极品一区二区| 中文字幕熟女人妻在线| 亚洲三级黄色毛片| 色综合亚洲欧美另类图片| 三级国产精品片| 国产精品久久久久久久电影| 国产免费视频播放在线视频 | 国产av在哪里看| 好男人视频免费观看在线| 国产精品人妻久久久久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 岛国在线免费视频观看| 日韩在线高清观看一区二区三区| 亚洲国产欧洲综合997久久,| 在线播放国产精品三级| 日韩精品有码人妻一区| 国产v大片淫在线免费观看| 精品人妻偷拍中文字幕| 国产精品熟女久久久久浪| 岛国毛片在线播放| 国产精品女同一区二区软件| 波多野结衣巨乳人妻| 国产精品一及| 国产精品久久久久久久久免| 国产亚洲5aaaaa淫片| 日韩亚洲欧美综合| 中文亚洲av片在线观看爽| 91在线精品国自产拍蜜月| 麻豆一二三区av精品| 最近最新中文字幕免费大全7| 亚洲欧美精品综合久久99| 国内少妇人妻偷人精品xxx网站| 精品久久久久久久末码| 色综合色国产| 久久精品久久精品一区二区三区| 国产精品99久久久久久久久| 日韩强制内射视频| 欧美潮喷喷水| 99在线人妻在线中文字幕| 国产av一区在线观看免费| 美女大奶头视频| 免费黄色在线免费观看| 国内少妇人妻偷人精品xxx网站| 性色avwww在线观看| 天天一区二区日本电影三级| 免费看美女性在线毛片视频| 日韩欧美在线乱码| 欧美性猛交黑人性爽| 男女那种视频在线观看| 美女黄网站色视频| 18禁动态无遮挡网站| 亚洲精品久久久久久婷婷小说 | 精品国产露脸久久av麻豆 | 日韩亚洲欧美综合| 国产极品精品免费视频能看的| 欧美zozozo另类| 亚洲国产日韩欧美精品在线观看| 日韩强制内射视频| 久久久欧美国产精品| 久久精品影院6| 国产综合懂色| 国产精品一及| 乱码一卡2卡4卡精品| 国内揄拍国产精品人妻在线| 精品久久久久久成人av| 久久久久精品久久久久真实原创| 国产精品久久久久久久久免| 精华霜和精华液先用哪个| a级毛色黄片| 中文字幕av成人在线电影| 国产精华一区二区三区| 精品熟女少妇av免费看| 禁无遮挡网站| 高清午夜精品一区二区三区| 日本一二三区视频观看| 欧美三级亚洲精品| 免费观看的影片在线观看| av福利片在线观看| 欧美成人一区二区免费高清观看| 国产高清不卡午夜福利| 午夜老司机福利剧场| 久久久久久久久久久免费av| 嫩草影院入口| 成年版毛片免费区| 一二三四中文在线观看免费高清| 色综合色国产| av卡一久久| 国产伦理片在线播放av一区| 欧美性猛交╳xxx乱大交人| 亚洲av中文字字幕乱码综合| 亚州av有码| 久久精品国产鲁丝片午夜精品| 国产午夜精品久久久久久一区二区三区| av国产久精品久网站免费入址| 麻豆一二三区av精品| 男人舔奶头视频| 偷拍熟女少妇极品色| 亚洲不卡免费看| 国产午夜福利久久久久久| 国产视频首页在线观看| 久久精品国产亚洲av天美| 欧美一区二区精品小视频在线| 一区二区三区四区激情视频| 久久久久久久亚洲中文字幕| 蜜臀久久99精品久久宅男| 毛片一级片免费看久久久久| 秋霞伦理黄片| 日韩精品有码人妻一区| av女优亚洲男人天堂| 欧美人与善性xxx| 精品午夜福利在线看| 欧美3d第一页| av在线蜜桃| 一区二区三区四区激情视频| 亚洲丝袜综合中文字幕| 亚洲精品成人久久久久久| 久久久久网色| 边亲边吃奶的免费视频| 国产人妻一区二区三区在| 纵有疾风起免费观看全集完整版 | 麻豆乱淫一区二区| 亚洲乱码一区二区免费版| 国产精品国产三级国产专区5o | 久久久久久久久久久丰满| 午夜a级毛片| 麻豆av噜噜一区二区三区| av在线观看视频网站免费| a级毛色黄片| 免费看a级黄色片| 国产精品一二三区在线看| 免费看a级黄色片| 久久99热6这里只有精品| 自拍偷自拍亚洲精品老妇| 桃色一区二区三区在线观看| 日韩欧美精品v在线| 干丝袜人妻中文字幕| 深夜a级毛片| 国产精品伦人一区二区| 3wmmmm亚洲av在线观看| 直男gayav资源| 亚洲综合色惰| 国产免费又黄又爽又色| 国产精品日韩av在线免费观看| 欧美潮喷喷水| 我要看日韩黄色一级片| 长腿黑丝高跟| 赤兔流量卡办理| or卡值多少钱| 国产探花极品一区二区| 久久韩国三级中文字幕| 人体艺术视频欧美日本| 国产精品嫩草影院av在线观看| 大话2 男鬼变身卡| 久久精品夜色国产| 91aial.com中文字幕在线观看| 嘟嘟电影网在线观看| 欧美bdsm另类| 级片在线观看| АⅤ资源中文在线天堂| 97超碰精品成人国产| 我要看日韩黄色一级片| 中文资源天堂在线| 搡老妇女老女人老熟妇| 女人久久www免费人成看片 | 国产精品国产三级专区第一集| 国产精品一二三区在线看| 深夜a级毛片| 亚洲av电影不卡..在线观看| 亚洲av.av天堂| 国产成人a区在线观看| 精品午夜福利在线看| 欧美潮喷喷水| 久久综合国产亚洲精品| 精品酒店卫生间| 亚洲五月天丁香| 高清午夜精品一区二区三区| 久99久视频精品免费| 日日啪夜夜撸| 欧美成人精品欧美一级黄| 男人舔女人下体高潮全视频| 丝袜喷水一区| 国产一区二区亚洲精品在线观看| kizo精华| 亚洲国产精品成人综合色| 小说图片视频综合网站| 1024手机看黄色片| 波多野结衣高清无吗| 免费看av在线观看网站| 国产精品av视频在线免费观看| 乱码一卡2卡4卡精品| 国产亚洲91精品色在线| 亚洲国产日韩欧美精品在线观看| 女人十人毛片免费观看3o分钟| 久久久久国产网址| 成人无遮挡网站| 亚洲欧美精品自产自拍| 99热6这里只有精品| 少妇人妻一区二区三区视频| 国产乱人偷精品视频| 日韩欧美三级三区| 国产午夜精品一二区理论片| 91在线精品国自产拍蜜月| 麻豆乱淫一区二区| 毛片女人毛片| 看十八女毛片水多多多| 天天一区二区日本电影三级| 成人午夜高清在线视频| 亚洲国产高清在线一区二区三| 亚洲欧美中文字幕日韩二区| 久久久亚洲精品成人影院| 美女高潮的动态| 最新中文字幕久久久久| videos熟女内射| 变态另类丝袜制服| 麻豆成人av视频| av卡一久久| 国产成人精品久久久久久| 最近的中文字幕免费完整| 亚洲精品,欧美精品| 亚洲人成网站在线播| 久久久久久久久久黄片| 97超碰精品成人国产| 两个人视频免费观看高清| 国产精品1区2区在线观看.| 国产精品熟女久久久久浪| 亚洲精品日韩在线中文字幕| 中文字幕免费在线视频6| 免费电影在线观看免费观看| 又爽又黄无遮挡网站| eeuss影院久久| 国产乱人偷精品视频| 变态另类丝袜制服| 精品少妇黑人巨大在线播放 | 中文亚洲av片在线观看爽| 精品人妻熟女av久视频| www日本黄色视频网| 精品人妻一区二区三区麻豆| 一级毛片电影观看 | 我要看日韩黄色一级片| 一个人看的www免费观看视频| 午夜老司机福利剧场| 久久久久久九九精品二区国产| 欧美3d第一页| 国产私拍福利视频在线观看| 日韩欧美在线乱码| 国产精品久久久久久av不卡| 两个人的视频大全免费| 91精品国产九色| 丰满人妻一区二区三区视频av| 波多野结衣巨乳人妻| 一区二区三区四区激情视频| 寂寞人妻少妇视频99o| 人妻少妇偷人精品九色| eeuss影院久久| 成人美女网站在线观看视频| 中国国产av一级| 国产黄片美女视频| 国产精品一区www在线观看| 久久国产乱子免费精品| 99热这里只有是精品50| 我要搜黄色片| 高清午夜精品一区二区三区| 99热这里只有是精品在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 国产成年人精品一区二区| 91精品一卡2卡3卡4卡| 亚洲高清免费不卡视频| 国产爱豆传媒在线观看| 中文字幕免费在线视频6| 长腿黑丝高跟| 国内精品一区二区在线观看| 中文乱码字字幕精品一区二区三区 | 日本午夜av视频| 国产大屁股一区二区在线视频| 男女国产视频网站| 国产精品一区二区在线观看99 | 免费黄网站久久成人精品| 亚洲国产欧美在线一区| 不卡视频在线观看欧美| 精品欧美国产一区二区三| 深夜a级毛片| 亚洲人成网站高清观看| 观看美女的网站| 亚洲电影在线观看av| 久久久久性生活片| 国产成人91sexporn| 亚洲av中文字字幕乱码综合| 性插视频无遮挡在线免费观看| 成人毛片a级毛片在线播放| 女的被弄到高潮叫床怎么办| 九色成人免费人妻av| 特级一级黄色大片| 国产亚洲av嫩草精品影院| 国产午夜精品久久久久久一区二区三区| 乱码一卡2卡4卡精品| 欧美性猛交╳xxx乱大交人| 人体艺术视频欧美日本| 成年免费大片在线观看| 国产一区亚洲一区在线观看| 久久精品国产亚洲av涩爱| 午夜福利视频1000在线观看| 少妇人妻一区二区三区视频| 欧美一级a爱片免费观看看| 最近中文字幕高清免费大全6| 亚洲欧美精品自产自拍| 欧美成人免费av一区二区三区| 国产亚洲91精品色在线| 大香蕉久久网| 人妻少妇偷人精品九色| 国产伦精品一区二区三区视频9| 麻豆一二三区av精品| 久久午夜福利片| 日韩精品青青久久久久久| 精品人妻一区二区三区麻豆| 久久综合国产亚洲精品| 国产免费福利视频在线观看| 晚上一个人看的免费电影| 一边亲一边摸免费视频| 三级经典国产精品| 中文字幕精品亚洲无线码一区| 菩萨蛮人人尽说江南好唐韦庄 | 色哟哟·www| 亚洲成色77777| 麻豆乱淫一区二区| 日韩三级伦理在线观看| 亚洲国产精品合色在线| 51国产日韩欧美| 久久婷婷人人爽人人干人人爱| 97超碰精品成人国产| 美女cb高潮喷水在线观看| 亚洲自拍偷在线| 中文字幕久久专区| 男女视频在线观看网站免费| 亚洲精品乱码久久久久久按摩| av黄色大香蕉| 天美传媒精品一区二区| 国产爱豆传媒在线观看| 国产乱人偷精品视频| 久久综合国产亚洲精品| 国产精品一区二区性色av| 亚洲人成网站在线播| 亚洲最大成人av| 成年女人看的毛片在线观看| 中文字幕亚洲精品专区| 淫秽高清视频在线观看| 国产精品一区二区在线观看99 | videos熟女内射| 搡老妇女老女人老熟妇| 天天躁夜夜躁狠狠久久av| 精品久久久久久久久久久久久| 久久综合国产亚洲精品| 一区二区三区高清视频在线| 国产av一区在线观看免费| 在线观看美女被高潮喷水网站| 欧美精品一区二区大全| 高清在线视频一区二区三区 | 天堂av国产一区二区熟女人妻| 丰满少妇做爰视频| 91aial.com中文字幕在线观看| 99热精品在线国产| 欧美xxxx黑人xx丫x性爽| 日韩欧美精品v在线| 久久这里有精品视频免费| 一区二区三区免费毛片| 精品人妻视频免费看| a级一级毛片免费在线观看| 国产在视频线精品| 美女大奶头视频| 波多野结衣巨乳人妻| 免费观看的影片在线观看| 午夜a级毛片| 国产午夜精品论理片| 亚洲欧美日韩高清专用| 菩萨蛮人人尽说江南好唐韦庄 | or卡值多少钱| 三级国产精品欧美在线观看| 国产在线一区二区三区精 |