• <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ù)控制方式比較與研究
    欧美一区二区精品小视频在线| 熟女人妻精品中文字幕| 亚洲精品一区av在线观看| 久久精品国产亚洲av香蕉五月| a级毛片a级免费在线| 久久久久国内视频| 成熟少妇高潮喷水视频| 欧美激情久久久久久爽电影| av国产免费在线观看| 亚洲av日韩精品久久久久久密| 一本一本综合久久| 99热这里只有是精品在线观看| 日韩大尺度精品在线看网址| 不卡一级毛片| 欧美精品国产亚洲| 精品欧美国产一区二区三| 国产在视频线在精品| 老熟妇乱子伦视频在线观看| 精品国产三级普通话版| 夜夜爽天天搞| 亚洲图色成人| 又紧又爽又黄一区二区| 99热网站在线观看| 男人的好看免费观看在线视频| 亚洲不卡免费看| 亚洲欧美日韩东京热| 自拍偷自拍亚洲精品老妇| 国产精品人妻久久久影院| 精品久久久久久久末码| 国内久久婷婷六月综合欲色啪| 黄色一级大片看看| 男女做爰动态图高潮gif福利片| 日本爱情动作片www.在线观看 | 简卡轻食公司| 最新在线观看一区二区三区| 中文字幕免费在线视频6| 中文在线观看免费www的网站| 真实男女啪啪啪动态图| 免费在线观看日本一区| 真人做人爱边吃奶动态| 老司机深夜福利视频在线观看| 精品福利观看| 国产一级毛片七仙女欲春2| 国产三级中文精品| 午夜福利欧美成人| 高清日韩中文字幕在线| 欧美一区二区国产精品久久精品| 观看美女的网站| 国产成人影院久久av| 中文字幕精品亚洲无线码一区| 欧美中文日本在线观看视频| 九色国产91popny在线| 午夜精品一区二区三区免费看| 禁无遮挡网站| 美女 人体艺术 gogo| 国产在线精品亚洲第一网站| av黄色大香蕉| 黄色配什么色好看| 亚洲av中文av极速乱 | 能在线免费观看的黄片| 18禁黄网站禁片午夜丰满| 一级av片app| 欧美潮喷喷水| 全区人妻精品视频| av国产免费在线观看| 国产一区二区亚洲精品在线观看| 69人妻影院| av女优亚洲男人天堂| 免费黄网站久久成人精品| av专区在线播放| 黄色配什么色好看| 两个人的视频大全免费| 日韩精品有码人妻一区| 精品久久久久久久久av| 欧美黑人巨大hd| 一本一本综合久久| 色av中文字幕| 赤兔流量卡办理| 99久久九九国产精品国产免费| 日韩中文字幕欧美一区二区| 精品久久久久久久久亚洲 | 国内久久婷婷六月综合欲色啪| 免费在线观看日本一区| 一个人看视频在线观看www免费| 免费观看在线日韩| 国产男靠女视频免费网站| 成人国产麻豆网| 国产成人福利小说| 日本a在线网址| 亚洲最大成人手机在线| 99热6这里只有精品| 日本三级黄在线观看| 国产av麻豆久久久久久久| 日韩欧美在线乱码| 国产白丝娇喘喷水9色精品| 欧美一区二区精品小视频在线| 日韩欧美 国产精品| 精品人妻视频免费看| 高清日韩中文字幕在线| 国产精品爽爽va在线观看网站| 男女之事视频高清在线观看| 色精品久久人妻99蜜桃| 日韩亚洲欧美综合| 国产中年淑女户外野战色| 亚洲精品456在线播放app | 一区福利在线观看| 久久久久久久久久黄片| 亚洲熟妇熟女久久| 偷拍熟女少妇极品色| 岛国在线免费视频观看| 欧美一区二区亚洲| 69人妻影院| 色播亚洲综合网| av视频在线观看入口| 亚洲国产高清在线一区二区三| 黄色欧美视频在线观看| 国产精品一区二区免费欧美| 午夜久久久久精精品| 中文在线观看免费www的网站| 联通29元200g的流量卡| 成人一区二区视频在线观看| 亚洲在线观看片| 十八禁国产超污无遮挡网站| 女人被狂操c到高潮| 日本撒尿小便嘘嘘汇集6| 麻豆成人午夜福利视频| 女人被狂操c到高潮| 欧美另类亚洲清纯唯美| 欧美国产日韩亚洲一区| 在线播放国产精品三级| 欧美另类亚洲清纯唯美| 熟女人妻精品中文字幕| 俺也久久电影网| 午夜爱爱视频在线播放| 成人性生交大片免费视频hd| 亚洲欧美清纯卡通| 99riav亚洲国产免费| 久久久久九九精品影院| 一边摸一边抽搐一进一小说| 亚洲精华国产精华精| 亚洲欧美日韩卡通动漫| 欧美性感艳星| 噜噜噜噜噜久久久久久91| 国产av不卡久久| 波多野结衣高清作品| 亚洲av美国av| 香蕉av资源在线| 国产午夜精品久久久久久一区二区三区 | 欧美一级a爱片免费观看看| 黄色丝袜av网址大全| 精品久久久久久成人av| 99热这里只有精品一区| 亚洲图色成人| 日本一本二区三区精品| 国产在视频线在精品| 国产 一区 欧美 日韩| 免费观看精品视频网站| 亚洲人成伊人成综合网2020| 欧美色欧美亚洲另类二区| 蜜桃久久精品国产亚洲av| 成熟少妇高潮喷水视频| 久久久久久久亚洲中文字幕| 精品久久久久久久久久免费视频| 99热这里只有是精品50| 精品99又大又爽又粗少妇毛片 | 日本-黄色视频高清免费观看| 婷婷亚洲欧美| 欧美精品啪啪一区二区三区| 日韩精品中文字幕看吧| 免费av不卡在线播放| 国产精品乱码一区二三区的特点| 国产欧美日韩精品一区二区| 在线观看舔阴道视频| 看免费成人av毛片| 内射极品少妇av片p| 亚洲电影在线观看av| 美女黄网站色视频| 精品久久国产蜜桃| 黄色视频,在线免费观看| 成人特级av手机在线观看| 欧美黑人欧美精品刺激| 高清日韩中文字幕在线| 波多野结衣高清作品| 18禁在线播放成人免费| 91在线精品国自产拍蜜月| 成人特级av手机在线观看| 亚洲av一区综合| 亚洲精华国产精华精| 成年人黄色毛片网站| 别揉我奶头~嗯~啊~动态视频| 在现免费观看毛片| 久久国内精品自在自线图片| 国产成人av教育| 婷婷丁香在线五月| 国产成人a区在线观看| 亚洲在线观看片| 美女被艹到高潮喷水动态| 在线观看舔阴道视频| 国产白丝娇喘喷水9色精品| 亚洲中文字幕日韩| 亚洲人成网站在线播放欧美日韩| 一个人看视频在线观看www免费| 身体一侧抽搐| 亚洲成人精品中文字幕电影| 丰满的人妻完整版| 国产欧美日韩精品亚洲av| 校园人妻丝袜中文字幕| а√天堂www在线а√下载| 亚洲人成伊人成综合网2020| 国内久久婷婷六月综合欲色啪| 久久国产精品人妻蜜桃| 久久久久久国产a免费观看| 久久精品国产清高在天天线| 久久久久久久久久黄片| 国产精品精品国产色婷婷| 欧美成人性av电影在线观看| 中国美女看黄片| 欧美xxxx黑人xx丫x性爽| 亚洲国产精品久久男人天堂| 久久久久久久久中文| 亚洲av二区三区四区| 国产男靠女视频免费网站| 国产男人的电影天堂91| 校园人妻丝袜中文字幕| 国产真实乱freesex| 久久国产精品人妻蜜桃| 黄色配什么色好看| 91午夜精品亚洲一区二区三区 | 国产麻豆成人av免费视频| 别揉我奶头 嗯啊视频| 国产伦在线观看视频一区| 国产极品精品免费视频能看的| 美女大奶头视频| 国产欧美日韩精品亚洲av| 一边摸一边抽搐一进一小说| 日本撒尿小便嘘嘘汇集6| 精品久久久噜噜| 国产精品国产高清国产av| 神马国产精品三级电影在线观看| av国产免费在线观看| 天堂动漫精品| 嫁个100分男人电影在线观看| 成年免费大片在线观看| 国产伦精品一区二区三区四那| 精品福利观看| 日本五十路高清| 天堂动漫精品| 禁无遮挡网站| 黄色丝袜av网址大全| 成人精品一区二区免费| 亚洲专区国产一区二区| 日本色播在线视频| 丰满乱子伦码专区| 真人一进一出gif抽搐免费| 午夜福利成人在线免费观看| 亚洲中文日韩欧美视频| 国产爱豆传媒在线观看| 波多野结衣高清作品| 国产高清有码在线观看视频| 性欧美人与动物交配| 看片在线看免费视频| 久久精品人妻少妇| 国产精品,欧美在线| 久久亚洲真实| 日本三级黄在线观看| 少妇熟女aⅴ在线视频| 神马国产精品三级电影在线观看| 国产精品,欧美在线| 久久天躁狠狠躁夜夜2o2o| 久久精品国产亚洲av天美| 日本熟妇午夜| 久久精品人妻少妇| 九色成人免费人妻av| 69av精品久久久久久| 中文字幕精品亚洲无线码一区| 天堂网av新在线| 亚洲成人免费电影在线观看| 久久这里只有精品中国| 一个人观看的视频www高清免费观看| 91麻豆av在线| 国产激情偷乱视频一区二区| 九九在线视频观看精品| eeuss影院久久| 少妇猛男粗大的猛烈进出视频 | 国产色婷婷99| 日韩中文字幕欧美一区二区| 日本精品一区二区三区蜜桃| 亚洲最大成人手机在线| 精品国内亚洲2022精品成人| 久久欧美精品欧美久久欧美| 日本色播在线视频| 琪琪午夜伦伦电影理论片6080| 久99久视频精品免费| 中文字幕高清在线视频| av在线蜜桃| 美女大奶头视频| 制服丝袜大香蕉在线| 亚洲精品456在线播放app | 国产91精品成人一区二区三区| a级毛片免费高清观看在线播放| 久久九九热精品免费| 日本三级黄在线观看| 天天一区二区日本电影三级| 波多野结衣高清无吗| 成年女人永久免费观看视频| 69av精品久久久久久| 日本黄色片子视频| 人人妻人人看人人澡| 波多野结衣高清无吗| 成年女人永久免费观看视频| 午夜激情福利司机影院| 国产伦精品一区二区三区四那| 999久久久精品免费观看国产| 国产一区二区在线观看日韩| 久久精品影院6| 久久久久九九精品影院| 国产成人a区在线观看| 欧美日韩黄片免| 丰满人妻一区二区三区视频av| 国产精品久久视频播放| 亚洲av不卡在线观看| 天天一区二区日本电影三级| 最新中文字幕久久久久| 免费搜索国产男女视频| 91精品国产九色| 日韩,欧美,国产一区二区三区 | 波多野结衣高清无吗| 精品久久久久久久久av| 在线观看av片永久免费下载| 久久精品国产清高在天天线| 欧美成人a在线观看| 蜜桃久久精品国产亚洲av| 欧美激情久久久久久爽电影| 亚洲国产精品sss在线观看| a级毛片免费高清观看在线播放| 久99久视频精品免费| 亚洲av成人精品一区久久| 中文亚洲av片在线观看爽| 黄色女人牲交| 禁无遮挡网站| 精品久久久久久久久久久久久| 国产精品98久久久久久宅男小说| 国产不卡一卡二| 欧美性猛交╳xxx乱大交人| 久久精品国产亚洲网站| 色综合婷婷激情| 国产一区二区三区av在线 | 国产激情偷乱视频一区二区| 精品一区二区三区人妻视频| 国产精品人妻久久久影院| 麻豆国产av国片精品| 人妻丰满熟妇av一区二区三区| 12—13女人毛片做爰片一| 成人特级av手机在线观看| 国产欧美日韩一区二区精品| 国产精品亚洲美女久久久| 精品人妻视频免费看| 噜噜噜噜噜久久久久久91| 91精品国产九色| 少妇熟女aⅴ在线视频| 日韩欧美一区二区三区在线观看| 日本 欧美在线| 日韩大尺度精品在线看网址| 亚洲中文字幕一区二区三区有码在线看| 色吧在线观看| 久久九九热精品免费| 精品久久久噜噜| 国产女主播在线喷水免费视频网站 | 别揉我奶头~嗯~啊~动态视频| 日韩欧美在线乱码| 亚洲成人精品中文字幕电影| 久久午夜福利片| 国产伦精品一区二区三区四那| 午夜福利在线在线| 女人十人毛片免费观看3o分钟| 亚洲av成人精品一区久久| 国产真实伦视频高清在线观看 | 精品不卡国产一区二区三区| 五月玫瑰六月丁香| 免费看av在线观看网站| 国产高清不卡午夜福利| 亚洲熟妇熟女久久| 99热6这里只有精品| 亚洲av二区三区四区| 亚洲精品国产成人久久av| 精品99又大又爽又粗少妇毛片 | 久久久久国内视频| 日日摸夜夜添夜夜添小说| 精品人妻1区二区| 在线观看午夜福利视频| 可以在线观看毛片的网站| 午夜影院日韩av| 97热精品久久久久久| 91在线观看av| 色综合亚洲欧美另类图片| x7x7x7水蜜桃| 成人高潮视频无遮挡免费网站| 久久久久久伊人网av| 久久午夜亚洲精品久久| 亚洲中文日韩欧美视频| 国产精品av视频在线免费观看| 日本欧美国产在线视频| 久久久久久久久久成人| 日韩欧美三级三区| 一级黄片播放器| 亚洲av二区三区四区| 亚洲最大成人av| 欧美另类亚洲清纯唯美| 亚洲第一区二区三区不卡| 嫩草影院精品99| 乱码一卡2卡4卡精品| 亚洲精品一区av在线观看| 国产亚洲av嫩草精品影院| 国产精品人妻久久久久久| 久久久久久伊人网av| 国产成人aa在线观看| 波野结衣二区三区在线| 久久久精品欧美日韩精品| 黄色欧美视频在线观看| 久久这里只有精品中国| 成人二区视频| 在线天堂最新版资源| 亚洲av一区综合| 精品国内亚洲2022精品成人| 一区福利在线观看| 18禁黄网站禁片午夜丰满| 久久精品91蜜桃| 男女之事视频高清在线观看| 免费人成在线观看视频色| 深夜a级毛片| 熟女人妻精品中文字幕| 在线观看免费视频日本深夜| 我要搜黄色片| 看黄色毛片网站| 精品一区二区免费观看| 午夜激情欧美在线| 国产亚洲精品久久久久久毛片| а√天堂www在线а√下载| 成人综合一区亚洲| 精品久久久久久久久亚洲 | 婷婷亚洲欧美| 欧美性猛交黑人性爽| 午夜久久久久精精品| 国产 一区精品| 国产精品美女特级片免费视频播放器| 精品福利观看| 国产精品精品国产色婷婷| 国产色婷婷99| 精品国产三级普通话版| 中文字幕精品亚洲无线码一区| 亚洲电影在线观看av| 日本免费a在线| 精品久久久噜噜| 深夜a级毛片| 日韩 亚洲 欧美在线| 国产精品久久久久久久电影| 天美传媒精品一区二区| 精品久久久久久久人妻蜜臀av| 国产精品人妻久久久影院| 婷婷丁香在线五月| 日韩精品青青久久久久久| 色精品久久人妻99蜜桃| 国产亚洲精品av在线| 女生性感内裤真人,穿戴方法视频| 亚洲avbb在线观看| 波多野结衣高清无吗| 午夜福利视频1000在线观看| 欧美激情久久久久久爽电影| 无遮挡黄片免费观看| 美女 人体艺术 gogo| 欧美区成人在线视频| 淫秽高清视频在线观看| 亚洲在线自拍视频| 午夜福利欧美成人| 别揉我奶头~嗯~啊~动态视频| 波野结衣二区三区在线| 女人十人毛片免费观看3o分钟| 日日夜夜操网爽| 男人舔女人下体高潮全视频| 99热这里只有是精品在线观看| 国产精品综合久久久久久久免费| 一本久久中文字幕| 欧美绝顶高潮抽搐喷水| 亚洲无线在线观看| 成人综合一区亚洲| 搡女人真爽免费视频火全软件 | 久久久精品欧美日韩精品| 国产色婷婷99| 精品一区二区免费观看| 亚洲精品456在线播放app | 午夜精品久久久久久毛片777| 网址你懂的国产日韩在线| 校园人妻丝袜中文字幕| av专区在线播放| 国产精品免费一区二区三区在线| 欧美极品一区二区三区四区| 日本与韩国留学比较| 三级毛片av免费| 成人av一区二区三区在线看| 天天躁日日操中文字幕| 亚洲国产精品合色在线| 色综合站精品国产| 成人美女网站在线观看视频| 免费av毛片视频| 欧美另类亚洲清纯唯美| 国产 一区精品| 麻豆av噜噜一区二区三区| 在线观看av片永久免费下载| 国产亚洲精品久久久com| 色综合色国产| 校园人妻丝袜中文字幕| 中文资源天堂在线| 久久久久久大精品| 久久久久久久亚洲中文字幕| 亚洲av.av天堂| 久久久久久久久久久丰满 | 免费搜索国产男女视频| 国产精品亚洲美女久久久| 噜噜噜噜噜久久久久久91| 亚洲av中文字字幕乱码综合| 亚洲自偷自拍三级| 午夜久久久久精精品| 久久久成人免费电影| 国产精品伦人一区二区| 欧美又色又爽又黄视频| 少妇人妻精品综合一区二区 | 国产精品自产拍在线观看55亚洲| 亚洲第一电影网av| 草草在线视频免费看| 人妻久久中文字幕网| 3wmmmm亚洲av在线观看| 久久草成人影院| 亚洲自拍偷在线| 中出人妻视频一区二区| 五月伊人婷婷丁香| 成熟少妇高潮喷水视频| 精品人妻视频免费看| 在线看三级毛片| 久久欧美精品欧美久久欧美| 国产精品美女特级片免费视频播放器| 12—13女人毛片做爰片一| 久久国内精品自在自线图片| 国产老妇女一区| 日本 欧美在线| av天堂中文字幕网| 一个人观看的视频www高清免费观看| 日韩国内少妇激情av| 搡老妇女老女人老熟妇| 久久久久久久久久成人| 亚洲av熟女| 特级一级黄色大片| 欧美性猛交黑人性爽| 亚洲中文日韩欧美视频| 欧美日韩中文字幕国产精品一区二区三区| 99精品在免费线老司机午夜| 亚洲内射少妇av| 免费人成在线观看视频色| 天堂√8在线中文| 欧美潮喷喷水| 国产一区二区在线av高清观看| 亚洲人成网站高清观看| 波多野结衣高清作品| 国产aⅴ精品一区二区三区波| 日本 av在线| 亚洲avbb在线观看| 久久这里只有精品中国| 亚洲不卡免费看| 久久亚洲精品不卡| 日本三级黄在线观看| 国产精品久久久久久精品电影| 内地一区二区视频在线| 国产主播在线观看一区二区| 99热这里只有精品一区| 少妇被粗大猛烈的视频| 精品99又大又爽又粗少妇毛片 | 两人在一起打扑克的视频| 日本与韩国留学比较| 亚洲av不卡在线观看| 日韩一区二区视频免费看| 哪里可以看免费的av片| 韩国av一区二区三区四区| 狠狠狠狠99中文字幕| 久久香蕉精品热| x7x7x7水蜜桃| 亚洲精品成人久久久久久| 欧美丝袜亚洲另类 | 亚洲欧美日韩高清专用| 人妻夜夜爽99麻豆av| 成人av在线播放网站| 日韩一本色道免费dvd| 狂野欧美白嫩少妇大欣赏| 人妻夜夜爽99麻豆av| 精品人妻熟女av久视频| 人人妻,人人澡人人爽秒播| 久久久久性生活片| 午夜影院日韩av| 中出人妻视频一区二区| 男人舔女人下体高潮全视频| 欧美国产日韩亚洲一区| 国产91精品成人一区二区三区| 十八禁国产超污无遮挡网站| 亚洲18禁久久av| 午夜视频国产福利| 美女xxoo啪啪120秒动态图| 一个人看视频在线观看www免费| 亚洲av第一区精品v没综合| 不卡视频在线观看欧美| 少妇猛男粗大的猛烈进出视频 | 国产精品一区二区性色av| 国产人妻一区二区三区在| 久久欧美精品欧美久久欧美| 麻豆av噜噜一区二区三区| 久久国产乱子免费精品| 亚洲精品乱码久久久v下载方式| 色吧在线观看|