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

    超聲速可調(diào)進(jìn)氣道內(nèi)流雙解現(xiàn)象及其節(jié)流特性

    2023-06-28 09:12:14金毅孫姝郭赟杰譚慧俊張悅
    航空學(xué)報(bào) 2023年7期
    關(guān)鍵詞:進(jìn)氣道喉道流態(tài)

    金毅,孫姝,郭赟杰,譚慧俊,張悅

    1.南京航空航天大學(xué) 能源與動(dòng)力學(xué)院 江蘇省航空動(dòng)力系統(tǒng)重點(diǎn)實(shí)驗(yàn)室,南京 210016 2.南京航空航天大學(xué) 民航學(xué)院,南京 211106

    一般地,當(dāng)進(jìn)氣道能通過調(diào)節(jié)自身幾何型面改變其進(jìn)氣量或內(nèi)流波系結(jié)構(gòu)以匹配飛行條件和發(fā)動(dòng)機(jī)工況時(shí),稱其為可調(diào)進(jìn)氣道。對(duì)于在寬?cǎi)R赫數(shù)范圍內(nèi)工作的超聲速飛行器而言,可調(diào)進(jìn)氣道能兼顧低馬赫數(shù)下的自起動(dòng)能力和高馬赫數(shù)下的氣動(dòng)性能,是工程應(yīng)用中的一個(gè)重要發(fā)展趨勢(shì)。目前可調(diào)進(jìn)氣道主要通過旋轉(zhuǎn)[1-2]、平移[3]唇罩或調(diào)節(jié)喉道[4-8]等方式實(shí)現(xiàn)。其中喉道調(diào)節(jié)作為最常見的調(diào)節(jié)方式,能輕易改變進(jìn)氣道內(nèi)通道中的激波位置和內(nèi)收縮比(ICR),使進(jìn)氣道達(dá)到最佳工作性能,得到了廣泛的應(yīng)用。

    由于可調(diào)進(jìn)氣道在喉道調(diào)節(jié)過程中往往會(huì)盡可能增大ICR以提高其臨界性能,因此極易面臨大ICR導(dǎo)致的不起動(dòng)問題[9-11]。當(dāng)不起動(dòng)現(xiàn)象出現(xiàn)后需立即減小ICR使進(jìn)氣道實(shí)現(xiàn)再起動(dòng)。然而由于喉道調(diào)節(jié)過程存在遲滯現(xiàn)象[11],導(dǎo)致進(jìn)氣道的再起動(dòng)ICR滯后于不起動(dòng)ICR,且其再起動(dòng)過程中出現(xiàn)了豐富的流態(tài)變化。進(jìn)氣道的再起動(dòng)過程實(shí)質(zhì)上是口部分離包逐步被吞入內(nèi)通道的動(dòng)態(tài)演化過程。具體地,隨著ICR減小,進(jìn)氣道依次經(jīng)歷口部大尺度分離誘導(dǎo)的不起動(dòng)流態(tài)、內(nèi)收縮段中小尺度分離誘導(dǎo)的非設(shè)計(jì)流態(tài)和再起動(dòng)流態(tài)。特別地,在進(jìn)氣道內(nèi)收縮段中小尺度分離誘導(dǎo)的非設(shè)計(jì)流態(tài)下,雖內(nèi)收縮段中存在偏離設(shè)計(jì)流態(tài)的復(fù)雜波系結(jié)構(gòu),但其并未影響進(jìn)氣道的口部流量捕獲特性。因此根據(jù)經(jīng)典教科書中的定義[12]判斷進(jìn)氣道此時(shí)處于起動(dòng)狀態(tài)。該類特殊的非設(shè)計(jì)流態(tài)往往出現(xiàn)在可調(diào)超聲速進(jìn)氣道的再起動(dòng)過程中[11,13-14],甚至在定幾何進(jìn)氣道[15]的節(jié)流變化過程中同樣存在此類流態(tài)。因此該類流態(tài)的內(nèi)流結(jié)構(gòu)及其工作性能值得關(guān)注和研究。

    在進(jìn)氣道的性能參數(shù)中,下游節(jié)流特性直接決定了進(jìn)氣道的臨界總壓恢復(fù)性能和穩(wěn)定裕度,是評(píng)價(jià)進(jìn)氣道工作性能優(yōu)劣的重要標(biāo)準(zhǔn)之一。學(xué)者們通過在出口設(shè)置堵錐[16]或擋板[17]的方式對(duì)超聲速進(jìn)氣道的下游節(jié)流特性開展了大量試驗(yàn)或仿真研究;結(jié)果表明下游節(jié)流造成的擾動(dòng)以結(jié)尾激波的形式向上游傳播[18-23]。當(dāng)來流馬赫數(shù)較低(一般小于2)時(shí)結(jié)尾激波為單道正激波[20],當(dāng)來流馬赫數(shù)較高(一般大于2.5)時(shí)結(jié)尾激波以正激波串[21]或斜激波串[16-17,19,22-23]的形式存在。此外當(dāng)結(jié)尾激波頭波位于喉道附近時(shí)進(jìn)氣道處于臨界狀態(tài)。一旦下游節(jié)流度進(jìn)一步增加,結(jié)尾激波就不能穩(wěn)定存在于進(jìn)氣道內(nèi)收縮段中,且會(huì)被迅速推出進(jìn)氣道口部,從而形成大范圍的不起動(dòng)流場(chǎng)結(jié)構(gòu)。然而上述研究均是針對(duì)設(shè)計(jì)流態(tài)下的固定幾何進(jìn)氣道開展的,其與可調(diào)進(jìn)氣道在非設(shè)計(jì)流態(tài)下的節(jié)流特性的區(qū)別仍有待研究。

    針對(duì)以上問題,本文設(shè)計(jì)并加工了一個(gè)典型的矩形喉道可調(diào)超聲速進(jìn)氣道,并在其出口設(shè)置了可移動(dòng)堵錐以模擬其下游節(jié)流過程。同時(shí)利用高速紋影觀測(cè)系統(tǒng)和動(dòng)態(tài)壓力測(cè)量系統(tǒng)在來流馬赫數(shù)為2.9的自由射流超聲速風(fēng)洞中開展了試驗(yàn)研究。筆者先對(duì)比研究超聲速可調(diào)進(jìn)氣道在設(shè)計(jì)和非設(shè)計(jì)流態(tài)下的通流內(nèi)流結(jié)構(gòu),接著對(duì)比研究設(shè)計(jì)和非設(shè)計(jì)流態(tài)下進(jìn)氣道的節(jié)流特性,最后對(duì)比研究設(shè)計(jì)和非設(shè)計(jì)流態(tài)下進(jìn)氣道節(jié)流過程及流動(dòng)機(jī)理。

    1 研究對(duì)象

    研究對(duì)象為喉道可調(diào)的矩形超聲速進(jìn)氣道,其試驗(yàn)?zāi)P腿鐖D1所示。進(jìn)氣道主體由外壓縮面、唇罩、側(cè)板、壓縮面?zhèn)群淼揽烧{(diào)機(jī)構(gòu)和泄流系統(tǒng)組成;其中進(jìn)氣道兩側(cè)為豎直側(cè)板,在試驗(yàn)中具備便于流動(dòng)光學(xué)檢測(cè)和易于模型安裝的優(yōu)點(diǎn)。此外進(jìn)氣道的喉道調(diào)節(jié)過程主要通過四連桿機(jī)構(gòu)實(shí)現(xiàn)。同時(shí)在進(jìn)氣道出口處設(shè)置可移動(dòng)堵錐以模擬其下游節(jié)流過程。進(jìn)氣道具體設(shè)計(jì)參數(shù)如下:工作馬赫數(shù)范圍為0~4、激波封口馬赫數(shù)為4、外壓縮面楔角為9°、唇罩壓縮角為7°、總長(zhǎng)為365.20 mm、捕獲高度為40.42 mm、入口高度為24.16 mm、出口高度為27 mm、矩形內(nèi)通道的寬度為40 mm??烧{(diào)喉道的高度能在8.22~19.97 mm變化,對(duì)應(yīng)的ICR變化范圍為1.21~2.94。由于研究的可調(diào)超聲速進(jìn)氣道具有工作ICR大和工作馬赫數(shù)范圍寬的特點(diǎn),其內(nèi)收縮段中的激波/邊界層干擾現(xiàn)象嚴(yán)重且內(nèi)流波系結(jié)構(gòu)變化范圍大。因此出于實(shí)際工程應(yīng)用的控制需求,內(nèi)收縮段壓縮面?zhèn)炔贾昧朔植际叫沽飨到y(tǒng)。該泄流系統(tǒng)由分布式泄流縫、獨(dú)立集氣腔和限流喉道組成,共使用了16條流向?qū)挾葹?.5 mm、展向長(zhǎng)度為32 mm的分布式泄流縫。當(dāng)通道內(nèi)的低能流經(jīng)過泄流縫后由2個(gè)獨(dú)立集氣腔收集,隨后從寬度為2 mm的限流喉道排出。此外為獲得進(jìn)氣道的動(dòng)態(tài)壓力數(shù)據(jù),在試驗(yàn)?zāi)P屯鈮嚎s面和唇罩上分別開設(shè)了C1~C13和R1~R12共25個(gè)動(dòng)態(tài)壓力測(cè)點(diǎn);同時(shí)為獲得進(jìn)氣道流場(chǎng)的試驗(yàn)紋影,在側(cè)板上開設(shè)了尺寸為145 mm×29 mm的矩形光學(xué)觀察窗。而且為保證進(jìn)氣道內(nèi)通道的氣密性,在試驗(yàn)?zāi)P瓦B接處均進(jìn)行了密封處理。

    圖1 試驗(yàn)?zāi)P虵ig.1 Test model

    在風(fēng)洞試驗(yàn)過程中,進(jìn)氣道的喉道調(diào)節(jié)機(jī)構(gòu)和可移動(dòng)堵錐分別由獨(dú)立的直線步進(jìn)電機(jī)驅(qū)動(dòng)。在喉道調(diào)節(jié)過程中ICR定義為進(jìn)氣道入口面積和可調(diào)喉道面積之比。在下游節(jié)流過程中對(duì)堵塞度(TR)的定義為

    式中:Ath,plug為進(jìn)氣道出口和堵錐之間的最小截面積;Aout為進(jìn)氣道的出口面積。此外為觀察到節(jié)流過程中的穩(wěn)定流場(chǎng)結(jié)構(gòu),試驗(yàn)中使用臺(tái)階式調(diào)節(jié)喉道和進(jìn)錐的方式,對(duì)應(yīng)的ICR和堵塞度變化示意圖如圖2所示,圖中t為時(shí)間,電機(jī)信號(hào)的上升沿表示堵錐固定于當(dāng)前ICR或TR,下降沿表示喉部四連桿驅(qū)動(dòng)機(jī)構(gòu)或堵錐以12 mm/s的速度向上游運(yùn)動(dòng)。

    圖2 試驗(yàn)過程中ICR和TR變化示意圖Fig.2 Schematic diagrams of variations of ICR and TR during test

    2 試驗(yàn)條件及測(cè)試手段

    試驗(yàn)均在南京航空航天大學(xué)的自由射流超聲速風(fēng)洞中完成。該風(fēng)洞運(yùn)行系統(tǒng)主要由入口處的蝶閥、整流段、名義馬赫數(shù)為3.0的拉瓦爾噴管、測(cè)試段、擴(kuò)壓器、閘板閥、真空罐和真空泵組成;其中方形拉瓦爾噴管出口尺寸為200 mm×200 mm,測(cè)試段兩側(cè)開設(shè)有尺寸為280 mm×205 mm的矩形光學(xué)觀察窗,真空罐的體積為400 m3。風(fēng)洞運(yùn)行過程如下:首先開啟真空泵將真空罐內(nèi)抽至真空條件,隨后開啟閘板閥使測(cè)試段內(nèi)同樣達(dá)到真空條件,緊接著開啟入口處的蝶閥,大氣環(huán)境的高壓氣流依次經(jīng)過整流段、拉瓦爾噴管、測(cè)試段、擴(kuò)壓器向真空罐內(nèi)的低壓區(qū)流動(dòng),由此在測(cè)試段內(nèi)建立了穩(wěn)定的超聲速流場(chǎng)。同時(shí)試驗(yàn)過程中記錄了風(fēng)洞測(cè)試段的來流條件,其中風(fēng)洞的實(shí)際馬赫數(shù)為2.9,來流總壓和總溫分別為102.6 kPa和290.5 K,風(fēng)洞穩(wěn)定運(yùn)行時(shí)間不少于20 s。該風(fēng)洞運(yùn)行系統(tǒng)的可靠性在前期的研究[24-26]中已經(jīng)得到了充分驗(yàn)證。

    此外為獲得進(jìn)氣道的壁面壓力數(shù)據(jù)和流場(chǎng)結(jié)構(gòu),試驗(yàn)中分別使用了動(dòng)態(tài)壓力采集系統(tǒng)和高速紋影觀測(cè)系統(tǒng)。在動(dòng)態(tài)壓力采集系統(tǒng)中使用了昆山雙橋傳感器有限公司生產(chǎn)的固有頻率為50 kHz、滿量程為100 kPa、測(cè)量精度為0.1%滿量程的CYG-503A型動(dòng)態(tài)壓力傳感器。其采集的動(dòng)態(tài)壓力信號(hào)由與之匹配的DAQ PCI-6255型數(shù)據(jù)采集卡獲得,該采集卡最多能以256 kB/s的采樣速率同時(shí)采集40個(gè)通道的壓力數(shù)據(jù),試驗(yàn)中采集卡的采樣率為20 kHz。需要注意的是C6、R7和R12等動(dòng)態(tài)壓力測(cè)點(diǎn)損壞,因此在后續(xù)分析過程中將其剔除。在高速紋影觀測(cè)系統(tǒng)中選擇由Nikon公司生產(chǎn)的MEMRECAM HX-3型高速攝影儀,試驗(yàn)過程中將其幀率設(shè)置為6 000、分辨率為1 152×336、快門時(shí)間為2 μs。同時(shí)為更清晰地顯示進(jìn)氣道的邊界層及流動(dòng)分離現(xiàn)象,紋影采集過程中均使用水平刀口。此外試驗(yàn)中使用外觸發(fā)信號(hào)觸發(fā)高速攝影儀,同時(shí)該信號(hào)被動(dòng)態(tài)壓力采集系統(tǒng)記錄以實(shí)現(xiàn)動(dòng)態(tài)壓力采集系統(tǒng)和高速紋影觀測(cè)系統(tǒng)的同步,便于后續(xù)數(shù)據(jù)分析。

    3 進(jìn)氣道內(nèi)流雙解現(xiàn)象

    在試驗(yàn)過程中進(jìn)氣道ICR的變化通過調(diào)節(jié)喉道面積而實(shí)現(xiàn)。初始狀態(tài)下進(jìn)氣道處于小ICR的通流起動(dòng)狀態(tài),隨后正向調(diào)節(jié)喉道以增加ICR,當(dāng)ICR過大導(dǎo)致進(jìn)氣道不起動(dòng)后反向調(diào)節(jié)喉道以減小ICR,直至進(jìn)氣道再起動(dòng)。在上述過程中,當(dāng)ICR=1.79~2.04時(shí)出現(xiàn)了流動(dòng)響應(yīng)遲滯現(xiàn)象。其中當(dāng)ICR=1.79時(shí)存在設(shè)計(jì)流態(tài)和非設(shè)計(jì)流態(tài)的特殊雙解流場(chǎng)。因此選擇ICR=1.79的進(jìn)氣道流態(tài)開展針對(duì)性研究。

    進(jìn)氣道設(shè)計(jì)和非設(shè)計(jì)流態(tài)下通流流場(chǎng)的試驗(yàn)紋影如圖3所示,對(duì)應(yīng)的壁面測(cè)點(diǎn)壓力分布曲線如圖4所示,圖中p為壓力,x為距離。首先進(jìn)氣道設(shè)計(jì)流態(tài)下的通流流場(chǎng)結(jié)構(gòu)如圖3(a)所示。其主要由外壓縮激波①、唇罩激波②、口部壓縮激波③和內(nèi)通道反射波系④組成。當(dāng)來流經(jīng)過口部激波系減速后其靜壓顯著增加,C3測(cè)點(diǎn)靜壓達(dá)到了約6.2倍來流靜壓(p0);隨后由于受肩部膨脹的影響,C4測(cè)點(diǎn)靜壓下降至4.8p0;接著在內(nèi)通道中反射激波的作用下,C5測(cè)點(diǎn)靜壓增至6.1p0;最后在喉道下游的擴(kuò)張段中超聲速氣流逐漸加速,其沿程靜壓呈下降趨勢(shì)。該流態(tài)下進(jìn)氣道內(nèi)流波系結(jié)構(gòu)穩(wěn)定,且內(nèi)通道中無明顯分離,是一類常見的設(shè)計(jì)流態(tài)。

    圖3 進(jìn)氣道雙解通流流場(chǎng)(ICR=1.79)Fig.3 Dual-solution inlet unthrottled flowfields(ICR=1.79)

    圖4 進(jìn)氣道通流流場(chǎng)壁面壓力分布Fig.4 Surface pressure distributions of inlet unthrottled flowfields

    其次進(jìn)氣道非設(shè)計(jì)流態(tài)下的通流流場(chǎng)結(jié)構(gòu)如圖3(b)所示。由于內(nèi)收縮段C1測(cè)點(diǎn)附近存在局部唇罩側(cè)分離包⑤,導(dǎo)致其分離激波與唇罩激波②匯聚形成了一道更強(qiáng)的入射激波。該入射激波與口部壓縮激波③相交后內(nèi)收縮段壓縮面?zhèn)刃纬闪司植狂R赫桿結(jié)構(gòu)⑥,且出現(xiàn)了小范圍的壓縮面?zhèn)确蛛x包⑦。同樣地,上述流場(chǎng)結(jié)構(gòu)在進(jìn)氣道內(nèi)通道中也形成了反射波系④。受唇罩側(cè)分離包的分離激波和再附激波的影響,C1和C2測(cè)點(diǎn)的靜壓分別高達(dá)8.9p0和9.0p0;C3測(cè)點(diǎn)靜壓則受肩部膨脹的影響逐漸下降至6.5p0;隨后C4測(cè)點(diǎn)靜壓在內(nèi)通道中反射激波的作用下增至7.6p0;最后超聲速氣流在擴(kuò)張段中逐漸加速,且沿程靜壓呈下降趨勢(shì)。該狀態(tài)下內(nèi)收縮段中存在多處激波/邊界層干擾導(dǎo)致的分離區(qū),是一類特殊的非設(shè)計(jì)流態(tài)。

    接著對(duì)該類特殊的非設(shè)計(jì)流態(tài)能穩(wěn)定存在于內(nèi)收縮段中的原因開展分析。當(dāng)口部分離包進(jìn)入內(nèi)通道時(shí),泄流系統(tǒng)能通過排除邊界層分離內(nèi)的低能流大幅減小壓縮面?zhèn)鹊姆蛛x尺寸,使內(nèi)收縮段中的分離區(qū)偏向唇罩側(cè)。此時(shí)唇罩側(cè)分離的肩部位置存在一個(gè)氣動(dòng)喉道,且該氣動(dòng)喉道的面積小于幾何喉道。但從該氣動(dòng)喉道所在位置的波系結(jié)構(gòu)可知主流仍然為超聲速流動(dòng)(并未減速至聲速),表明氣動(dòng)喉道并未導(dǎo)致主流發(fā)生流動(dòng)壅塞現(xiàn)象。此外唇罩側(cè)分離和口部波系結(jié)構(gòu)共同誘導(dǎo)出了壓縮面?zhèn)鹊木植狂R赫桿結(jié)構(gòu),并在壓縮面?zhèn)刃纬闪诉h(yuǎn)小于唇罩側(cè)分離尺寸的小范圍分離現(xiàn)象。該壓縮面?zhèn)确蛛x誘導(dǎo)的再附激波恰好入射到唇罩側(cè)分離的再附點(diǎn)附近,導(dǎo)致唇罩側(cè)分離下游形成了強(qiáng)逆壓力梯度,阻止了唇罩側(cè)分離繼續(xù)向下游運(yùn)動(dòng)。也就是說在有限的通道內(nèi),唇罩側(cè)分離誘導(dǎo)的復(fù)雜波系結(jié)構(gòu)反作用于該分離包自身,形成了穩(wěn)定“自持”的流動(dòng)結(jié)構(gòu)。

    需要注意的是非設(shè)計(jì)流態(tài)下內(nèi)收縮段中的壓力較高,其喉道靜壓(pth)和泄流量均大于設(shè)計(jì)值。此外由于內(nèi)收縮段中的復(fù)雜波系造成了額外總壓損失,非設(shè)計(jì)流態(tài)下進(jìn)氣道喉道總壓(p*th)小于設(shè)計(jì)值。結(jié)合式(2)中的理想氣體總、靜壓關(guān)系式可得非設(shè)計(jì)流態(tài)下進(jìn)氣道喉道馬赫數(shù)(Math)小于設(shè)計(jì)值。

    式中:k為理想氣體比熱比。

    圖4補(bǔ)充給出了兩類進(jìn)氣道流場(chǎng)壁面壓力測(cè)點(diǎn)的脈動(dòng)幅度對(duì)比。顯然非設(shè)計(jì)流態(tài)下的內(nèi)通道流場(chǎng)脈動(dòng)幅度遠(yuǎn)大于設(shè)計(jì)流態(tài)下的脈動(dòng)幅度。以內(nèi)收縮段中壓力脈動(dòng)幅度最大的C1測(cè)點(diǎn)為例,使用傅里葉變換和連續(xù)小波變換進(jìn)行了頻域分析(如圖5所示),圖中功率譜密度(PSD)以其均方幅值(MSA)表征,f為頻率。同時(shí)借助Morlet小波[27-28]以提高對(duì)頻率的分辨能力,其波數(shù)設(shè)置為20。從圖5(a)中可看出非設(shè)計(jì)流態(tài)下的進(jìn)氣道發(fā)生了明顯的流場(chǎng)振蕩,其呈現(xiàn)出了較寬的頻率分布范圍,對(duì)應(yīng)的主頻為72.9 Hz;而設(shè)計(jì)流態(tài)下的進(jìn)氣道并未發(fā)生明顯的流場(chǎng)振蕩現(xiàn)象。圖5(b)中非設(shè)計(jì)流態(tài)下的小波頻譜分析可進(jìn)一步證明其流場(chǎng)振蕩呈現(xiàn)出了寬頻特征,且頻率主要集中在100 Hz以下。這是因?yàn)镃1測(cè)點(diǎn)附近分離包自身具備振蕩特性,導(dǎo)致內(nèi)收縮段中出現(xiàn)了一定頻率的小幅振蕩;此外下壁面?zhèn)鹊木植狂R赫桿結(jié)構(gòu)具有不穩(wěn)定特性,導(dǎo)致振蕩呈現(xiàn)出了寬頻特征。

    圖5 通流狀態(tài)下進(jìn)氣道C1壓力測(cè)點(diǎn)的頻域分析Fig.5 Frequency domain analyses based on pressure signal slices from C1 of unthrottled inlet

    4 進(jìn)氣道節(jié)流特性

    第4節(jié)主要對(duì)第3節(jié)設(shè)計(jì)和非設(shè)計(jì)流態(tài)下進(jìn)氣道的節(jié)流特性開展對(duì)比研究。由于進(jìn)氣道構(gòu)型兩側(cè)為豎直側(cè)板,一旦內(nèi)通道中的結(jié)尾激波越過側(cè)板最前緣位置,那么其兩側(cè)將產(chǎn)生溢流,從而改變口部流量捕獲特性。因此當(dāng)監(jiān)測(cè)到外壓縮面最下游測(cè)點(diǎn)R11產(chǎn)生清晰的壓力上升信號(hào)時(shí)判斷進(jìn)氣道開始陷入不起動(dòng)狀態(tài),此時(shí)對(duì)應(yīng)的堵塞度為臨界堵塞度,以TRC表示。此外定義在不起動(dòng)狀態(tài)出現(xiàn)的上一時(shí)刻進(jìn)氣道依然能維持穩(wěn)定起動(dòng)的狀態(tài)為臨界狀態(tài);定義臨界狀態(tài)下進(jìn)氣道的臨界壓比為出口靜壓(pC,以C13測(cè)點(diǎn)標(biāo)定)與來流靜壓(p0)之比,以pC/p0表示。

    圖6給出了節(jié)流過程中的進(jìn)氣道壁面壓力分布曲線。對(duì)于設(shè)計(jì)流態(tài)下的節(jié)流過程(圖6(a)),當(dāng)TR增加到32.4%時(shí)節(jié)流造成的壓力擾動(dòng)傳遞到C7測(cè)點(diǎn)附近;隨后當(dāng)TR逐步增加到41.7%時(shí)擾動(dòng)傳遞到喉道附近;下一時(shí)刻,當(dāng)TR為42.4%時(shí)C1及其下游測(cè)點(diǎn)靜壓均顯著增加,且外壓縮面上R11測(cè)點(diǎn)的靜壓升高,表明節(jié)流造成的擾動(dòng)已經(jīng)影響到了進(jìn)氣道的口部波系結(jié)構(gòu),進(jìn)氣道開始陷入不起動(dòng)狀態(tài),該堵塞度對(duì)應(yīng)為臨界堵塞度;當(dāng)TR為45.0%時(shí)外壓縮面上大部分測(cè)點(diǎn)靜壓遠(yuǎn)高于通流狀態(tài),進(jìn)氣道不起動(dòng)程度加劇。由于試驗(yàn)中采用臺(tái)階式進(jìn)錐的方式,取TR為41.7%的穩(wěn)定臺(tái)階下進(jìn)氣道的節(jié)流狀態(tài)為臨界狀態(tài),此時(shí)對(duì)應(yīng)的臨界壓比為15.8。對(duì)于非設(shè)計(jì)流態(tài)下的節(jié)流過程(圖6(b)),雖內(nèi)通道中的初始?jí)荷^高,但其擾動(dòng)傳播過程和設(shè)計(jì)流態(tài)下相似。具體地在TR由0增加至40.0%的過程中擾動(dòng)逐漸從進(jìn)氣道出口傳遞至喉道附近;隨后當(dāng)TR為41.7%時(shí)內(nèi)通道中所有測(cè)點(diǎn)靜壓增加,且R11測(cè)點(diǎn)的靜壓升高,進(jìn)氣道開始陷入不起動(dòng)狀態(tài),該堵塞度對(duì)應(yīng)為臨界堵塞度;當(dāng)TR為45.0%時(shí)進(jìn)氣道不起動(dòng)程度加劇。取TR為40.0%的穩(wěn)定臺(tái)階下進(jìn)氣道的節(jié)流狀態(tài)為臨界狀態(tài),此時(shí)對(duì)應(yīng)的臨界壓比為16.0。

    圖6 節(jié)流過程中的進(jìn)氣道壁面壓力分布曲線Fig.6 Surface pressure distribution curves of inlet during throttling processes

    作為總結(jié),表1列出了設(shè)計(jì)和非設(shè)計(jì)流態(tài)的進(jìn)氣道節(jié)流特性對(duì)比,其中pth/p0為進(jìn)氣道通流狀態(tài)下的喉道壓比,定義為喉道靜壓(pth,以C4測(cè)點(diǎn)標(biāo)定)和來流靜壓(p0)之比。不難看出非設(shè)計(jì)流態(tài)下內(nèi)收縮段中的復(fù)雜流場(chǎng)使其喉道壓比高于設(shè)計(jì)流態(tài)。而且雖非設(shè)計(jì)流態(tài)下喉道上游的復(fù)雜波系結(jié)構(gòu)造成了額外的總壓損失,但其整體節(jié)流性能與設(shè)計(jì)流態(tài)相當(dāng)(臨界壓比分別為15.8和16.0)。

    表1 設(shè)計(jì)和非設(shè)計(jì)流態(tài)的進(jìn)氣道節(jié)流特性對(duì)比Table 1 Comparison of throttling characteristics of inlet under designed and undesigned conditions

    5 進(jìn)氣道節(jié)流過程及流動(dòng)機(jī)理

    本節(jié)分析設(shè)計(jì)和非設(shè)計(jì)流態(tài)下進(jìn)氣道的節(jié)流過程及流動(dòng)機(jī)理。首先進(jìn)氣道設(shè)計(jì)和非設(shè)計(jì)流態(tài)下節(jié)流過程中的壁面壓力變化曲線如圖7所示。選擇唇罩最上游測(cè)點(diǎn)C1、喉道稍上游測(cè)點(diǎn)C3、唇罩最下游測(cè)點(diǎn)C13和外壓縮面最下游測(cè)點(diǎn)R11等關(guān)鍵測(cè)點(diǎn)進(jìn)行對(duì)比分析,同時(shí)給出了堵錐電機(jī)的運(yùn)動(dòng)信號(hào)及其穩(wěn)定臺(tái)階對(duì)應(yīng)的堵塞度作為參考。對(duì)于設(shè)計(jì)流態(tài)下的節(jié)流過程,當(dāng)TR達(dá)16.8%時(shí)C13測(cè)點(diǎn)壓力開始呈現(xiàn)出隨TR增大而增大的趨勢(shì);至TR達(dá)到41.7%的過程中,圖7中其余測(cè)點(diǎn)壓力值均保持不變;接著當(dāng)TR為42.4%時(shí),R11測(cè)點(diǎn)感受到了明顯壓升,喉道上游C1、C3測(cè)點(diǎn)壓力產(chǎn)生了突變,表明進(jìn)氣道開始陷入不起動(dòng)狀態(tài);隨后進(jìn)一步增加TR,各測(cè)點(diǎn)振蕩幅度逐漸增加,進(jìn)氣道不起動(dòng)程度加劇。非設(shè)計(jì)流態(tài)下C1和C3測(cè)點(diǎn)的靜壓和脈動(dòng)幅值均高于設(shè)計(jì)流態(tài),但其節(jié)流過程中的各測(cè)點(diǎn)靜壓變化趨勢(shì)均與設(shè)計(jì)流態(tài)相似。當(dāng)TR>TRC時(shí)C1、C3和R11等壓力測(cè)點(diǎn)均產(chǎn)生突變,進(jìn)氣道陷入不起動(dòng)狀態(tài)。這表明兩類流態(tài)的進(jìn)氣道下游節(jié)流程度一旦超過某一臨界值將迅速陷入不起動(dòng)狀態(tài),即其均為“突變型”的節(jié)流不起動(dòng)模式。

    圖7 進(jìn)氣道節(jié)流過程中的壁面壓力變化曲線Fig.7 Surface pressure time histories during throttling processes

    兩類流態(tài)的進(jìn)氣道在節(jié)流過程中的試驗(yàn)紋影如圖8所示。對(duì)于設(shè)計(jì)流態(tài)下的節(jié)流過程,下游節(jié)流造成的擾動(dòng)以結(jié)尾激波串的形式向上游傳播。為便于分析,圖8中標(biāo)出了喉道下游的背景激波a及其反射激波b、c,同時(shí)標(biāo)出了壓縮面?zhèn)鹊呐蛎洸╠。當(dāng)TR達(dá)32.4%時(shí)紋影觀察窗內(nèi)開始觀察到結(jié)尾激波串頭波,且頭波呈“X”形,其上、下半支基本對(duì)稱分布。由于內(nèi)通道下壁面邊界層較厚,其抵抗逆壓力梯度的能力較差,當(dāng)背景激波b入射在擴(kuò)張段下壁面時(shí)會(huì)形成具有較大初始逆壓力梯度的干擾區(qū)。在TR逐漸增至36.7%的過程中激波串頭波下半支接觸到了背景激波b在擴(kuò)張段下壁面造成的干擾區(qū),隨后下壁面分離將快速向上游增長(zhǎng)。而頭波上半支受到了膨脹波d的影響,其邊界層分離發(fā)展受限。因此該堵塞度下激波串核心區(qū)向上壁面偏轉(zhuǎn)。隨著TR繼續(xù)增加,激波串頭波上半支接觸到背景激波a入射在上壁面造成的激波/邊界層干擾區(qū),上壁面邊界層分離快速增長(zhǎng)。在此過程中激波串的非對(duì)稱性減弱。當(dāng)TR達(dá)41.7%時(shí)結(jié)尾激波串頭波運(yùn)動(dòng)到喉道附近,激波串頭波形態(tài)基本呈對(duì)稱分布。進(jìn)一步增加TR,結(jié)尾激波串不能穩(wěn)定存在于內(nèi)收縮段中,其迅速向上游運(yùn)動(dòng)。當(dāng)TR達(dá)42.4%時(shí)結(jié)尾激波串頭波造成的初始分離位置剛好位于側(cè)板前緣處,內(nèi)通道完全被結(jié)尾激波造成的低能流占據(jù),此時(shí)進(jìn)氣道開始陷入不起動(dòng)狀態(tài)。當(dāng)TR為45.0%時(shí)進(jìn)氣道處于不起動(dòng)狀態(tài),此時(shí)外壓縮面上出現(xiàn)了大尺寸邊界層分離,其誘導(dǎo)的分離激波導(dǎo)致口部的唇罩側(cè)出現(xiàn)了溢流。

    圖8 進(jìn)氣道節(jié)流過程中的試驗(yàn)紋影Fig.8 Schlieren images of inlet during throttling processes

    對(duì)于非設(shè)計(jì)流態(tài)下的節(jié)流過程,下游節(jié)流造成的擾動(dòng)同樣以結(jié)尾激波串的形式向上游傳播,且一旦結(jié)尾激波串越過喉道,其將迅速向上游運(yùn)動(dòng)至進(jìn)氣道陷入不起動(dòng)狀態(tài)。需要注意的是在節(jié)流過程中設(shè)計(jì)和非設(shè)計(jì)流態(tài)下結(jié)尾激波串的分布形態(tài)存在區(qū)別。具體地對(duì)比兩類流態(tài)的進(jìn)氣道在節(jié)流過程中TR為36.7%的流場(chǎng),可看出設(shè)計(jì)流態(tài)下節(jié)流過程中的結(jié)尾激波串頭波下半支明顯前伸,然而非設(shè)計(jì)流態(tài)下節(jié)流過程中的結(jié)尾激波串頭波上、下半支基本呈對(duì)稱分布。這是因?yàn)榉窃O(shè)計(jì)流態(tài)下的進(jìn)氣道喉道馬赫數(shù)較低,一方面使喉道下游擴(kuò)張段中的反射激波/邊界層干擾強(qiáng)度減弱,另一方面也使結(jié)尾激波串上游馬赫數(shù)降低,從而減小了結(jié)尾激波串的強(qiáng)度。

    基于上述認(rèn)識(shí),可從流動(dòng)損失的角度進(jìn)一步分析在設(shè)計(jì)和非設(shè)計(jì)流態(tài)下進(jìn)氣道整體節(jié)流性能相當(dāng)?shù)臋C(jī)理。由于兩類進(jìn)氣道在臨界狀態(tài)下結(jié)尾激波串頭波均位于喉道附近,因此取該狀態(tài)進(jìn)行分析。此時(shí)整個(gè)進(jìn)氣道的流動(dòng)損失由喉道上游波系結(jié)構(gòu)和喉道下游結(jié)尾激波串共同造成。首先由于非設(shè)計(jì)流態(tài)下內(nèi)收縮段中存在復(fù)雜波系結(jié)構(gòu),因此其喉道上游波系結(jié)構(gòu)造成的流動(dòng)損失大于設(shè)計(jì)值。其次由于非設(shè)計(jì)流態(tài)下的通流流場(chǎng)喉道馬赫數(shù)小于設(shè)計(jì)值,使喉道下游結(jié)尾激波串造成的流動(dòng)損失小于設(shè)計(jì)值。因此在喉道上、下游波系結(jié)構(gòu)造成的流動(dòng)損失的共同影響下設(shè)計(jì)和非設(shè)計(jì)流態(tài)下整個(gè)進(jìn)氣道的流動(dòng)損失相近,進(jìn)而導(dǎo)致兩類進(jìn)氣道的整體節(jié)流性能相當(dāng)。

    此外由于激波串具備非定常特性[29-32],因此進(jìn)一步對(duì)設(shè)計(jì)和非設(shè)計(jì)流態(tài)下進(jìn)氣道在節(jié)流過程中的結(jié)尾激波串振蕩特性進(jìn)行分析。以進(jìn)氣道喉道下游的C7測(cè)點(diǎn)為例,對(duì)其動(dòng)態(tài)壓力信號(hào)開展頻域分析(如圖9所示)??煽闯鰧?duì)于設(shè)計(jì)流態(tài)下的節(jié)流過程,在TR為32.4%~41.7%的范圍區(qū)間內(nèi)C7測(cè)點(diǎn)能感受到結(jié)尾激波串的振蕩。其中當(dāng)TR在38.3%附近時(shí)脈動(dòng)能量最強(qiáng),此時(shí)振蕩呈現(xiàn)出寬頻特征,其頻率主要分布在600 Hz以下。對(duì)于非設(shè)計(jì)流態(tài)下的節(jié)流過程,C7測(cè)點(diǎn)僅能在TR為32.4%~34.9%時(shí)感受到結(jié)尾激波串的振蕩。此時(shí)振蕩同樣呈現(xiàn)出寬頻特征,但其頻率主要分布在400 Hz以下。也就是說在節(jié)流過程中設(shè)計(jì)和非設(shè)計(jì)流態(tài)下進(jìn)氣道結(jié)尾激波串的振蕩均具備寬頻特征,但后者的激波串的振蕩堵塞度范圍區(qū)間和頻率分布范圍區(qū)間均更窄。造成該現(xiàn)象的原因在于相比于設(shè)計(jì)流態(tài),非設(shè)計(jì)流態(tài)下進(jìn)氣道通流喉道馬赫數(shù)更低,致使喉道下游激波/邊界層干擾強(qiáng)度減弱、結(jié)尾激波串強(qiáng)度降低,進(jìn)而導(dǎo)致結(jié)尾激波串振蕩強(qiáng)度降低。

    圖9 進(jìn)氣道節(jié)流過程中C7壓力測(cè)點(diǎn)的功率譜密度Fig.9 Power spectra density based on pressure signal slices from C7 of inlet during throttling processes

    6 結(jié)論

    1) 試驗(yàn)結(jié)果表明在超聲速可調(diào)進(jìn)氣道的喉道正、反向調(diào)節(jié)過程中出現(xiàn)了設(shè)計(jì)和非設(shè)計(jì)流態(tài)的起動(dòng)雙解通流流場(chǎng)。其中設(shè)計(jì)流態(tài)下進(jìn)氣道內(nèi)流結(jié)構(gòu)正常建立,內(nèi)通道無明顯邊界層分離。但非設(shè)計(jì)流態(tài)下進(jìn)氣道內(nèi)收縮段中存在分離誘導(dǎo)的復(fù)雜波系,并形成了局部馬赫桿結(jié)構(gòu),導(dǎo)致其整體壓升相較于設(shè)計(jì)流態(tài)更高、喉道馬赫數(shù)更低,且存在100 Hz以下的寬頻、低頻小幅振蕩。

    2) 雖非設(shè)計(jì)流態(tài)下喉道上游的復(fù)雜波系結(jié)構(gòu)造成了額外的總壓損失,但其整體節(jié)流性能與設(shè)計(jì)流態(tài)相當(dāng)。設(shè)計(jì)和非設(shè)計(jì)流態(tài)的臨界堵塞度分別為42.4%和41.7%,臨界壓比分別為15.8和16.0。

    3) 在設(shè)計(jì)和非設(shè)計(jì)流態(tài)下進(jìn)氣道的下游節(jié)流過程中擾動(dòng)均以結(jié)尾斜激波串的形式向上游傳播,且當(dāng)其頭波位于喉道附近時(shí)進(jìn)氣道處于臨界狀態(tài)。此外設(shè)計(jì)流態(tài)下的結(jié)尾激波串出現(xiàn)了非對(duì)稱的分布特征,但非設(shè)計(jì)流態(tài)下的結(jié)尾激波串基本呈對(duì)稱分布。在振蕩特性方面,設(shè)計(jì)和非設(shè)計(jì)流態(tài)下進(jìn)氣道結(jié)尾激波串的振蕩均具備寬頻特征,但后者的振蕩堵塞度范圍區(qū)間和頻率分布范圍區(qū)間均更窄。

    猜你喜歡
    進(jìn)氣道喉道流態(tài)
    側(cè)邊機(jī)組故障對(duì)泵站前池流態(tài)的影響
    基于AVL-Fire的某1.5L發(fā)動(dòng)機(jī)進(jìn)氣道優(yōu)化設(shè)計(jì)
    基于輔助進(jìn)氣門的進(jìn)氣道/發(fā)動(dòng)機(jī)一體化控制
    改進(jìn)邊界條件的非恒定流模型在城市河流橡膠壩流態(tài)模擬中的應(yīng)用
    U型渠道無喉道量水槽流動(dòng)規(guī)律數(shù)值模擬
    勝利油田致密砂巖油藏微觀孔隙結(jié)構(gòu)特征
    亞聲速二喉道流場(chǎng)不對(duì)稱現(xiàn)象研究
    The coupling characteristics of supersonic dual inlets for missile①
    動(dòng)態(tài)流態(tài)冰蓄冷系統(tǒng)在千級(jí)凈化廠房的應(yīng)用
    基于TM遙感影像的河口流態(tài)信息半定量化研究
    国产片特级美女逼逼视频| 91成人精品电影| 精品视频人人做人人爽| 欧美xxⅹ黑人| 男人操女人黄网站| 青青草视频在线视频观看| 免费日韩欧美在线观看| 亚洲精品中文字幕在线视频| 亚洲av美国av| av网站免费在线观看视频| 亚洲欧洲国产日韩| 亚洲,一卡二卡三卡| 麻豆乱淫一区二区| 黄片小视频在线播放| 满18在线观看网站| 高清av免费在线| 亚洲欧美清纯卡通| 悠悠久久av| av视频免费观看在线观看| 精品国产乱码久久久久久男人| bbb黄色大片| 午夜免费鲁丝| 水蜜桃什么品种好| 国产欧美日韩一区二区三 | 男女边摸边吃奶| 一本—道久久a久久精品蜜桃钙片| 色视频在线一区二区三区| 国产精品99久久99久久久不卡| 中文字幕制服av| 男女午夜视频在线观看| 欧美日韩av久久| 91麻豆av在线| 香蕉国产在线看| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品一二三| 国产免费又黄又爽又色| av在线app专区| 美女视频免费永久观看网站| 飞空精品影院首页| 国产在线一区二区三区精| 国产免费又黄又爽又色| 国产91精品成人一区二区三区 | 人成视频在线观看免费观看| 久久精品亚洲熟妇少妇任你| 国产亚洲午夜精品一区二区久久| 亚洲av男天堂| 欧美日韩视频精品一区| 婷婷色麻豆天堂久久| 婷婷色综合大香蕉| 十分钟在线观看高清视频www| 免费av中文字幕在线| 91精品国产国语对白视频| 又黄又粗又硬又大视频| 精品人妻1区二区| 成人亚洲精品一区在线观看| 亚洲天堂av无毛| 精品欧美一区二区三区在线| 多毛熟女@视频| 天天操日日干夜夜撸| 嫩草影视91久久| 亚洲国产精品成人久久小说| 亚洲国产最新在线播放| 久久国产精品人妻蜜桃| 啦啦啦中文免费视频观看日本| 中文字幕精品免费在线观看视频| 午夜两性在线视频| 亚洲av在线观看美女高潮| 精品人妻在线不人妻| 美女高潮到喷水免费观看| 免费观看人在逋| 黄频高清免费视频| 狂野欧美激情性bbbbbb| 国产激情久久老熟女| 大香蕉久久成人网| 最近最新中文字幕大全免费视频 | 成人手机av| 国产成人一区二区三区免费视频网站 | 纵有疾风起免费观看全集完整版| 又大又黄又爽视频免费| 一本久久精品| 国产欧美日韩一区二区三区在线| 精品亚洲成国产av| 亚洲精品第二区| 男女高潮啪啪啪动态图| 美女主播在线视频| 欧美精品一区二区大全| 看十八女毛片水多多多| 国产有黄有色有爽视频| 电影成人av| 又大又黄又爽视频免费| 久久影院123| 老司机在亚洲福利影院| 亚洲成av片中文字幕在线观看| 两性夫妻黄色片| 日本av手机在线免费观看| 一级a爱视频在线免费观看| 午夜视频精品福利| 少妇 在线观看| 一本色道久久久久久精品综合| 国产精品一国产av| 这个男人来自地球电影免费观看| 天天躁夜夜躁狠狠久久av| 国产精品秋霞免费鲁丝片| 精品第一国产精品| 两个人看的免费小视频| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品久久久久久人妻精品电影 | 80岁老熟妇乱子伦牲交| 久久久亚洲精品成人影院| 免费久久久久久久精品成人欧美视频| 亚洲五月色婷婷综合| 午夜福利视频在线观看免费| 1024香蕉在线观看| 色婷婷久久久亚洲欧美| 99久久99久久久精品蜜桃| 中国国产av一级| 久久精品国产综合久久久| 国产精品一区二区在线不卡| 中国国产av一级| 国产精品久久久久久精品电影小说| 精品国产一区二区三区四区第35| 欧美在线一区亚洲| 亚洲综合色网址| 日日爽夜夜爽网站| 国产精品一区二区免费欧美 | 午夜福利乱码中文字幕| 狠狠精品人妻久久久久久综合| 亚洲精品在线美女| 成年动漫av网址| 久久精品成人免费网站| 婷婷丁香在线五月| 丁香六月欧美| 欧美乱码精品一区二区三区| 一级黄色大片毛片| 日韩 亚洲 欧美在线| 一本综合久久免费| av天堂久久9| 男女高潮啪啪啪动态图| 国产精品 欧美亚洲| 欧美日韩福利视频一区二区| 日韩伦理黄色片| 欧美大码av| 老汉色∧v一级毛片| 国产男女超爽视频在线观看| av福利片在线| 91麻豆av在线| 多毛熟女@视频| 国产熟女午夜一区二区三区| 人成视频在线观看免费观看| 男女床上黄色一级片免费看| 老汉色av国产亚洲站长工具| 久久精品国产亚洲av高清一级| 免费在线观看黄色视频的| 亚洲国产精品一区三区| 黄片播放在线免费| 国产精品秋霞免费鲁丝片| 亚洲精品自拍成人| 美女脱内裤让男人舔精品视频| 成人免费观看视频高清| 男女之事视频高清在线观看 | 国产精品国产三级国产专区5o| 久热这里只有精品99| 天堂中文最新版在线下载| 国产伦理片在线播放av一区| 一个人免费看片子| 青春草视频在线免费观看| 黄色视频不卡| 一本一本久久a久久精品综合妖精| cao死你这个sao货| 亚洲精品国产一区二区精华液| 国产精品国产av在线观看| 在线精品无人区一区二区三| 91精品三级在线观看| 色视频在线一区二区三区| 久久综合国产亚洲精品| 国产1区2区3区精品| 国产精品久久久久成人av| 一本大道久久a久久精品| 悠悠久久av| 国产伦人伦偷精品视频| 国产av国产精品国产| 丁香六月天网| 亚洲人成77777在线视频| 观看av在线不卡| 日韩一卡2卡3卡4卡2021年| 国产精品偷伦视频观看了| 国产精品一国产av| 极品人妻少妇av视频| 18禁黄网站禁片午夜丰满| 成年女人毛片免费观看观看9 | 亚洲色图 男人天堂 中文字幕| 国产精品av久久久久免费| 女性生殖器流出的白浆| 母亲3免费完整高清在线观看| 一区在线观看完整版| 国产97色在线日韩免费| 99九九在线精品视频| 啦啦啦啦在线视频资源| 亚洲精品第二区| 另类精品久久| 国产精品一二三区在线看| 亚洲av在线观看美女高潮| www.999成人在线观看| 一级,二级,三级黄色视频| 夫妻午夜视频| 高潮久久久久久久久久久不卡| 亚洲精品国产一区二区精华液| 国产精品一区二区在线不卡| 黄片播放在线免费| 欧美人与性动交α欧美软件| 免费日韩欧美在线观看| 国产高清视频在线播放一区 | 天天躁夜夜躁狠狠躁躁| av视频免费观看在线观看| 久久久精品免费免费高清| 国产精品香港三级国产av潘金莲 | 老司机午夜十八禁免费视频| 如日韩欧美国产精品一区二区三区| 黑人巨大精品欧美一区二区蜜桃| a级毛片黄视频| 免费在线观看日本一区| 777久久人妻少妇嫩草av网站| 老司机靠b影院| 两个人免费观看高清视频| 黄片小视频在线播放| 美女主播在线视频| 中文欧美无线码| 电影成人av| 中文字幕人妻丝袜制服| 丝袜人妻中文字幕| 亚洲av成人精品一二三区| 免费日韩欧美在线观看| 色网站视频免费| 黑人欧美特级aaaaaa片| 国产淫语在线视频| 亚洲免费av在线视频| 亚洲精品第二区| 一级毛片黄色毛片免费观看视频| 人人澡人人妻人| 国产精品麻豆人妻色哟哟久久| 男人操女人黄网站| 日韩一区二区三区影片| 最近最新中文字幕大全免费视频 | 久久久久精品人妻al黑| 男女国产视频网站| 精品亚洲乱码少妇综合久久| 国产欧美日韩综合在线一区二区| 亚洲欧美日韩另类电影网站| 国产亚洲精品久久久久5区| 黄片播放在线免费| 久久久欧美国产精品| 各种免费的搞黄视频| 999精品在线视频| 伊人亚洲综合成人网| 国产1区2区3区精品| 国产不卡av网站在线观看| 成人国语在线视频| 日韩制服丝袜自拍偷拍| 精品熟女少妇八av免费久了| 国产精品一二三区在线看| 久久99一区二区三区| 亚洲中文日韩欧美视频| 十八禁网站网址无遮挡| 久久午夜综合久久蜜桃| 少妇被粗大的猛进出69影院| 自拍欧美九色日韩亚洲蝌蚪91| 男人舔女人的私密视频| 免费看十八禁软件| www日本在线高清视频| 国产免费现黄频在线看| 日日爽夜夜爽网站| 欧美日本中文国产一区发布| 我的亚洲天堂| 两个人免费观看高清视频| 亚洲精品乱久久久久久| 国产高清视频在线播放一区 | 精品第一国产精品| kizo精华| 久久久久国产一级毛片高清牌| 91字幕亚洲| 免费高清在线观看视频在线观看| 欧美日韩精品网址| 国产av精品麻豆| 高潮久久久久久久久久久不卡| 人妻 亚洲 视频| 一级毛片 在线播放| 久久久久久久国产电影| 欧美日韩福利视频一区二区| 女人被躁到高潮嗷嗷叫费观| 欧美日韩成人在线一区二区| 99久久精品国产亚洲精品| 黄色 视频免费看| 国产女主播在线喷水免费视频网站| 精品高清国产在线一区| 成年人免费黄色播放视频| www.熟女人妻精品国产| 一区福利在线观看| 国产男女内射视频| 久久精品熟女亚洲av麻豆精品| 精品国产一区二区久久| 亚洲 国产 在线| 在线观看免费日韩欧美大片| 水蜜桃什么品种好| 天天躁夜夜躁狠狠久久av| 久久精品亚洲熟妇少妇任你| 波野结衣二区三区在线| 大片免费播放器 马上看| 欧美人与善性xxx| 亚洲精品国产色婷婷电影| 日韩 欧美 亚洲 中文字幕| e午夜精品久久久久久久| 欧美 日韩 精品 国产| 久久精品亚洲av国产电影网| 秋霞在线观看毛片| 丝袜人妻中文字幕| 亚洲三区欧美一区| 国产黄色免费在线视频| 97精品久久久久久久久久精品| 51午夜福利影视在线观看| 咕卡用的链子| 久久国产精品大桥未久av| 丁香六月欧美| 中文字幕人妻丝袜一区二区| 天堂8中文在线网| 午夜福利一区二区在线看| 一本—道久久a久久精品蜜桃钙片| 中文字幕制服av| 后天国语完整版免费观看| 午夜两性在线视频| 免费高清在线观看视频在线观看| 国产精品免费视频内射| 高清黄色对白视频在线免费看| 美国免费a级毛片| 波野结衣二区三区在线| 咕卡用的链子| 成人亚洲欧美一区二区av| 亚洲中文日韩欧美视频| 国产精品一区二区在线观看99| 看免费av毛片| 热99久久久久精品小说推荐| 亚洲国产精品一区三区| 在线精品无人区一区二区三| 亚洲精品久久午夜乱码| 人妻一区二区av| 欧美少妇被猛烈插入视频| 国产日韩欧美在线精品| 日韩中文字幕视频在线看片| 久久狼人影院| 久久久久国产一级毛片高清牌| 波多野结衣av一区二区av| 日本av免费视频播放| 建设人人有责人人尽责人人享有的| 久久精品成人免费网站| 少妇 在线观看| 日韩av不卡免费在线播放| 无限看片的www在线观看| 国产精品av久久久久免费| 午夜91福利影院| 天天添夜夜摸| 王馨瑶露胸无遮挡在线观看| 精品国产乱码久久久久久男人| 国语对白做爰xxxⅹ性视频网站| 欧美日韩综合久久久久久| 亚洲欧美一区二区三区久久| 少妇裸体淫交视频免费看高清 | a 毛片基地| 亚洲av欧美aⅴ国产| 色婷婷久久久亚洲欧美| 国产精品 国内视频| 国产片特级美女逼逼视频| 亚洲国产精品一区二区三区在线| 777久久人妻少妇嫩草av网站| 欧美精品一区二区大全| 秋霞在线观看毛片| 亚洲熟女精品中文字幕| 女人被躁到高潮嗷嗷叫费观| 久久久精品国产亚洲av高清涩受| 亚洲少妇的诱惑av| 精品欧美一区二区三区在线| 成人午夜精彩视频在线观看| 国产一区二区 视频在线| 黄片小视频在线播放| 青春草亚洲视频在线观看| a级片在线免费高清观看视频| 国产人伦9x9x在线观看| 亚洲,欧美,日韩| 日本av手机在线免费观看| 日韩,欧美,国产一区二区三区| 亚洲欧美日韩高清在线视频 | 蜜桃在线观看..| 伦理电影免费视频| 黄片小视频在线播放| 91老司机精品| 国产一区二区三区综合在线观看| 黄色怎么调成土黄色| 国产精品二区激情视频| 国产av国产精品国产| 亚洲av日韩精品久久久久久密 | 国产又色又爽无遮挡免| 久久综合国产亚洲精品| 欧美激情极品国产一区二区三区| 午夜日韩欧美国产| 国产成人精品久久久久久| 热99国产精品久久久久久7| 亚洲七黄色美女视频| 亚洲精品日本国产第一区| 老鸭窝网址在线观看| 一级a爱视频在线免费观看| 精品国产一区二区三区久久久樱花| 天堂8中文在线网| 久久精品国产亚洲av涩爱| 天堂8中文在线网| 啦啦啦 在线观看视频| 国产成人a∨麻豆精品| 久久久精品免费免费高清| 七月丁香在线播放| 欧美精品高潮呻吟av久久| 日韩制服骚丝袜av| 丝袜在线中文字幕| 午夜av观看不卡| 欧美日韩亚洲高清精品| 桃花免费在线播放| 青春草视频在线免费观看| svipshipincom国产片| 国产av一区二区精品久久| 久久午夜综合久久蜜桃| 黄色一级大片看看| 天天躁日日躁夜夜躁夜夜| 亚洲欧美精品自产自拍| 国产精品99久久99久久久不卡| 国产精品国产三级专区第一集| 婷婷成人精品国产| 精品久久久久久电影网| 久久狼人影院| 欧美久久黑人一区二区| 大香蕉久久网| 人妻人人澡人人爽人人| 亚洲一码二码三码区别大吗| 亚洲三区欧美一区| 国产国语露脸激情在线看| 人人妻,人人澡人人爽秒播 | 久久精品国产a三级三级三级| 一本—道久久a久久精品蜜桃钙片| 精品视频人人做人人爽| 久久久久久久久免费视频了| 久久女婷五月综合色啪小说| 欧美 亚洲 国产 日韩一| 国产女主播在线喷水免费视频网站| 侵犯人妻中文字幕一二三四区| 成年人黄色毛片网站| 久久久久精品国产欧美久久久 | 亚洲av日韩精品久久久久久密 | 性高湖久久久久久久久免费观看| 午夜福利影视在线免费观看| 国产成人精品久久二区二区91| 一级黄片播放器| 观看av在线不卡| 极品少妇高潮喷水抽搐| 男女边吃奶边做爰视频| 男女床上黄色一级片免费看| 亚洲欧美一区二区三区久久| 亚洲国产欧美一区二区综合| 国产无遮挡羞羞视频在线观看| 亚洲欧洲国产日韩| 青春草视频在线免费观看| 制服人妻中文乱码| 亚洲成色77777| www日本在线高清视频| 国产精品国产三级专区第一集| 午夜视频精品福利| 大码成人一级视频| 91麻豆精品激情在线观看国产 | 国产亚洲av高清不卡| 美女主播在线视频| 搡老岳熟女国产| 成人影院久久| 90打野战视频偷拍视频| 丝袜脚勾引网站| 香蕉国产在线看| av天堂在线播放| 国产精品三级大全| 91精品三级在线观看| 欧美日韩国产mv在线观看视频| 亚洲欧美一区二区三区黑人| 天堂俺去俺来也www色官网| 亚洲欧洲国产日韩| h视频一区二区三区| svipshipincom国产片| 亚洲精品美女久久av网站| 欧美日韩国产mv在线观看视频| 天天操日日干夜夜撸| 超碰97精品在线观看| 桃花免费在线播放| 亚洲视频免费观看视频| 国产1区2区3区精品| 大话2 男鬼变身卡| 久久久精品区二区三区| 国产极品粉嫩免费观看在线| 日韩电影二区| 一级毛片女人18水好多 | 99热全是精品| 日韩制服骚丝袜av| 国产精品久久久人人做人人爽| 丁香六月欧美| 欧美在线一区亚洲| 久久久久久人人人人人| 久久久国产一区二区| 一二三四社区在线视频社区8| 亚洲精品国产av成人精品| 这个男人来自地球电影免费观看| 老熟女久久久| 精品一区在线观看国产| 国产成人精品在线电影| 少妇 在线观看| 亚洲av成人精品一二三区| av欧美777| 啦啦啦在线免费观看视频4| 满18在线观看网站| 欧美成人精品欧美一级黄| 精品欧美一区二区三区在线| www.999成人在线观看| 婷婷色综合www| 亚洲五月婷婷丁香| 免费黄频网站在线观看国产| 青青草视频在线视频观看| 中文字幕另类日韩欧美亚洲嫩草| 久久精品久久久久久久性| 十分钟在线观看高清视频www| 欧美老熟妇乱子伦牲交| 美国免费a级毛片| 欧美日韩成人在线一区二区| 丁香六月天网| 人人澡人人妻人| 国产伦人伦偷精品视频| 婷婷成人精品国产| 又粗又硬又长又爽又黄的视频| a级片在线免费高清观看视频| av国产久精品久网站免费入址| 国产黄色免费在线视频| videos熟女内射| 欧美精品av麻豆av| 国产精品av久久久久免费| 日韩中文字幕欧美一区二区 | 国产午夜精品一二区理论片| 精品国产国语对白av| 色精品久久人妻99蜜桃| 嫩草影视91久久| 精品一区在线观看国产| 叶爱在线成人免费视频播放| 啦啦啦啦在线视频资源| xxx大片免费视频| 中文字幕人妻丝袜一区二区| 老司机影院成人| 午夜av观看不卡| 日韩 欧美 亚洲 中文字幕| 国精品久久久久久国模美| 99国产精品一区二区蜜桃av | 2018国产大陆天天弄谢| 女人久久www免费人成看片| 丝袜脚勾引网站| 亚洲av成人精品一二三区| 国产欧美日韩一区二区三 | 日韩一区二区三区影片| 亚洲精品乱久久久久久| 日韩,欧美,国产一区二区三区| 亚洲午夜精品一区,二区,三区| 女警被强在线播放| 亚洲视频免费观看视频| 欧美人与性动交α欧美软件| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美黑人精品巨大| 中文字幕高清在线视频| 超色免费av| 国产日韩欧美在线精品| 免费一级毛片在线播放高清视频 | 久久精品国产亚洲av高清一级| 日本一区二区免费在线视频| 欧美变态另类bdsm刘玥| 美国免费a级毛片| 后天国语完整版免费观看| 一区福利在线观看| netflix在线观看网站| 中文欧美无线码| 久久狼人影院| 精品国产国语对白av| 欧美精品一区二区免费开放| 美女中出高潮动态图| 国产1区2区3区精品| 国产精品国产三级专区第一集| 国产成人系列免费观看| 久久 成人 亚洲| 最新的欧美精品一区二区| 91成人精品电影| 亚洲 欧美一区二区三区| 日韩av不卡免费在线播放| 精品亚洲成国产av| 免费在线观看日本一区| 1024视频免费在线观看| 少妇精品久久久久久久| 久久 成人 亚洲| 欧美亚洲 丝袜 人妻 在线| 巨乳人妻的诱惑在线观看| 美女福利国产在线| 国产精品偷伦视频观看了| 美女中出高潮动态图| 精品国产国语对白av| 亚洲熟女精品中文字幕| 免费人妻精品一区二区三区视频| 亚洲第一av免费看| 精品少妇内射三级| 欧美日韩成人在线一区二区| 男人添女人高潮全过程视频| av国产久精品久网站免费入址| 亚洲熟女精品中文字幕| 少妇被粗大的猛进出69影院|