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

    總溫畸變下跨聲壓氣機失速過程非定常模擬

    2023-08-31 02:36:26趙紅亮張文強邱佳慧張敏杜娟聶超群
    航空學報 2023年14期
    關(guān)鍵詞:總溫總壓周向

    趙紅亮,張文強,邱佳慧,張敏,杜娟,*,聶超群

    1.華北電力大學 能源動力與機械工程學院,北京 102206

    2.北京理工大學 機電學院,北京 100081

    3.中國科學院 工程熱物理研究所 數(shù)字孿生研究中心,北京 100190

    4.中國科學院 先進能源動力重點實驗室,北京 100190

    5.中國科學院 輕型動力創(chuàng)新研究院,北京 100190

    6.中國科學院大學 工程科學學院,北京 100049

    航空發(fā)動機實際運行(如戰(zhàn)斗機發(fā)射導(dǎo)彈、飛機編隊飛行、艦載機在航母上垂直/短距起降、客機飛過閃電區(qū)等)中都會在發(fā)動機進口產(chǎn)生溫度畸變。吸入高溫燃氣或者蒸汽后會降低發(fā)動機穩(wěn)定裕度,易于引發(fā)失穩(wěn)[1-3]。研究總溫畸變對發(fā)動機(壓氣機)的穩(wěn)定性和性能影響具有重要意義。為此,國內(nèi)外的研究機構(gòu)和研究人員們利用試驗、理論分析和數(shù)值仿真手段對發(fā)動機(壓氣機)進氣總溫畸變展開了大量研究工作并積累了豐富成果。

    試驗方面,美國航空航天局劉易斯研究中心的Rudey 和Antl[4]利用氫氣燃燒器產(chǎn)生總溫畸變,對某雙軸渦扇發(fā)動機進行了總溫畸變試驗;試驗結(jié)果表明高壓壓氣機是對溫度畸變最敏感的部件,總溫畸變強度和高溫區(qū)范圍均會造成穩(wěn)定邊界下移。中國燃氣渦輪研究院的葉培梁和劉大響[5]對某小涵道比渦扇發(fā)動機進行了周向總溫畸變下穩(wěn)定性評定試驗;試驗確定了對該發(fā)動機的穩(wěn)定性影響最嚴重的一組溫度畸變參數(shù),發(fā)現(xiàn)失速首先發(fā)生在高壓壓氣機中。

    理論分析方法是基于壓氣機的數(shù)學模型把進氣溫度畸變作為輸入?yún)?shù),利用模型的輸出參數(shù)實現(xiàn)快速分析壓氣機性能變化。目前研究進氣畸變對壓氣機穩(wěn)定性影響的模型有平行壓氣機模型[6]、激盤模型[7]、體積力模型[8-9]等。美國阿諾德工程發(fā)展中心的Hale 等[10-12]基于體積力思想,將通過流線曲率法獲得的葉輪機械源項耦合到三維歐拉求解器中,發(fā)展了三維計算程序TEACC(Turbine Engine Analysis Compressor Code),利用該程序定量評估了進氣畸變對某三級軍用風扇性能和穩(wěn)定性的影響。郭晉等[13]通過結(jié)合三維體積力和經(jīng)典二維平行壓氣機模型提出了一種混合保真度計算模型;利用該模型以較低計算成本定量評估了復(fù)雜進氣畸變對某大涵道比發(fā)動機的性能影響,并捕捉到內(nèi)部流場的三維特征。李秋實等[14]利用體積力模型將葉片排和進氣道-風扇相互作用等作為源項耦合到COMSOL-CFD 程序中,建立了進氣道-風扇耦合模擬方法;該模型較好地模擬出了進氣道-風扇系統(tǒng)的特性性能及流場性能。葉巍等[15-16]通過前人提出的畸變傳遞模型[17]研究了進氣溫度畸變對某發(fā)動機穩(wěn)定性和性能的影響;結(jié)果表明穩(wěn)定裕度的損失更多取決于畸變強度,而不是畸變的周向范圍;較高轉(zhuǎn)速下發(fā)動機抗溫度畸變能力明顯強于中低轉(zhuǎn)速;在發(fā)射的溫度畸變強度下發(fā)動機已失去氣動穩(wěn)定性,必須使用防喘系統(tǒng)。張百靈等[18]以劉易斯研究中心設(shè)計的5 級壓氣機的前3 級為研究對象建立了一個壓縮系統(tǒng)的模型,結(jié)合體積力方法,通過求解非穩(wěn)態(tài)三維歐拉方程較為準確地計算了周向總溫畸變的傳遞特性;分析發(fā)現(xiàn)周向總溫畸變經(jīng)過壓氣機后強度衰減了19%,高溫區(qū)的周向位置偏移約80°,且總溫畸變強度在轉(zhuǎn)子內(nèi)發(fā)生局部增大現(xiàn)象,第1 級轉(zhuǎn)子失速的可能性最大。

    數(shù)值仿真方面,Weston 等[19]對進氣總壓畸變條件下的某3 級風扇展開全環(huán)數(shù)值仿真,研究表明經(jīng)過風扇后總壓畸變引發(fā)了總溫畸變,越接近失速工況點,總溫畸變強度越大,這與轉(zhuǎn)子葉片周圍輸入的功有關(guān),運行在近失速工況點時轉(zhuǎn)過全環(huán)過程中第1 級轉(zhuǎn)子的單個葉片做功量變化范圍可達25%。李志平等[20]利用數(shù)值模擬手段對某單級軸流壓氣機進行了不同畸變強度的壓力-溫度組合畸變對穩(wěn)定性和性能影響的研究;通過近失速工況點轉(zhuǎn)子葉頂區(qū)二維流線的分析發(fā)現(xiàn)組合畸變下壓氣機失穩(wěn)機理仍是泄漏流從轉(zhuǎn)子葉片前緣溢出。尤延鋮等[21]總結(jié)概述了國內(nèi)外在航空發(fā)動機進氣總溫畸變方面的研究工作,包括進氣總溫畸變來源、溫度畸變發(fā)生裝置及各類總溫畸變評定標準,最后列舉分析了針對總溫畸變的發(fā)動機防喘控制技術(shù)優(yōu)缺點,指出對進氣總溫畸變產(chǎn)生機理、傳遞和影響等研究急需加強。

    通過上述研究發(fā)現(xiàn)在進氣總溫畸變試驗方面主要是針對整機展開研究,但由于整機試驗成本昂貴、周期漫長等問題,近年來公開的文獻較少;采用理論分析可實現(xiàn)進氣畸變對發(fā)動機(壓氣機)穩(wěn)定性和性能影響的快速預(yù)測,但是對內(nèi)部流場細節(jié)研究有所不足;此外,目前總溫畸變下的高精度非定常計算工作較少,主要聚焦于壓氣機進口總壓畸變引發(fā)的總溫畸變或總壓總溫復(fù)合畸變條件下,畸變區(qū)在壓氣機內(nèi)部的傳播規(guī)律及對壓氣機性能的影響。而針對進氣總溫畸變條件下壓氣機性能變化的背后物理機理和失速現(xiàn)象有待進一步深入研究。相關(guān)試驗數(shù)據(jù)[22]表明徑向溫度畸變影響較小,通常忽略不計。因此本文對Darmstadt 跨聲單級壓氣機展開了周向總溫畸變條件下全環(huán)非定常數(shù)值模擬研究,分析失速過程中壓氣機內(nèi)部流場變化,以期揭示總溫畸變下該壓氣機發(fā)生失速的動力學機理。

    1 研究對象

    研究對象為德國Darmstadt 跨聲單級壓氣機試驗臺(如圖1 所示)。該壓氣機由德國MTU 航空發(fā)動機公司和德國航空太空中心設(shè)計并優(yōu)化[23-25],轉(zhuǎn)子葉片在輪轂附近彎曲度較高且在葉頂附近很薄,整個工作范圍內(nèi)靜子葉片沒有流動分離跡象且出口角分布均勻,是渦扇高壓壓氣機前級的典型代表。壓氣機采用800 kW 直流驅(qū)動,通過變速器改變軸速度,最大轉(zhuǎn)速為20 500 r/min,最大扭矩為350 N·m。壓氣機測量截面有級進口面(ME15)、轉(zhuǎn)子進口面(ME21)、級出口面(ME30),壓氣機性能特性是根據(jù)級進口面和級出口面的面平均參數(shù)計算得到的,已公布的試驗數(shù)據(jù)包含100%和65%設(shè)計轉(zhuǎn)速結(jié)果。表1總結(jié)了試驗臺的關(guān)鍵參數(shù),更詳細的數(shù)據(jù)見文獻[26-27]。

    表1 Darmstadt 跨聲壓氣機的關(guān)鍵參數(shù)Table 1 Key parameters of Darmstadt transonic compressor

    圖1 Darmstadt 跨聲壓氣機試驗臺Fig.1 Test rig of Darmstadt transonic compressor

    2 數(shù)值模擬方法

    2.1 數(shù)值設(shè)置

    通過Ansys CFX 軟件求解雷諾平均Navier-Stokes 方程。計算設(shè)置采用了標準k-ω湍流模型[28],壁面為絕熱無滑移,流體介質(zhì)為理想氣體。對于均勻來流,采用具有周期性邊界條件的單通道定常計算,進口邊界條件設(shè)置為公布的試驗測量值[27],即進口總壓為99 kPa、進口總溫為295 K、湍流強度為4%、湍流長度尺度為0.09 m。轉(zhuǎn)靜交界面類型為Stage(Mixing-Plane)。對進氣總溫畸變采用全環(huán)非定常計算,進口總溫設(shè)置為180°周向范圍的500 K 高總溫區(qū)和295 K 的低總溫區(qū)(如圖2 所示,圖中θ為周向位置),轉(zhuǎn)靜交界面類型為Transient Rotor Stator。出口邊界條件采用固定背壓為30 kPa 的噴嘴,通過調(diào)整噴嘴喉部大小獲取不同工況點的特性,大量研究[29-31]已證實此出口邊界條件的設(shè)置可在壓氣機較小流量工況點使計算穩(wěn)定,適合壓氣機失穩(wěn)計算。求解格式均為高階求解模式,非定常計算采用隱式雙時間步,轉(zhuǎn)子葉片轉(zhuǎn)過一個葉片通道所需的物理時間步長為80 步,每個物理時間步長的虛擬時間步長為10 步。

    圖2 全環(huán)計算域示意圖Fig.2 Schematic of annulus computation domain

    2.2 網(wǎng)格無關(guān)性

    葉片網(wǎng)格是通過NUMECA/AutoGrid5 生成,采用O4H 拓撲,根據(jù)文獻[32-33],網(wǎng)格最大壁面距離y+<3。針對單通道計算域進行了粗糙、中等、精細3 種不同網(wǎng)格方案的模擬測試。3 組網(wǎng)格方案的網(wǎng)格節(jié)點在流向、周向、徑向具體分布情況見表2。將100%轉(zhuǎn)速下計算特性和試驗數(shù)據(jù)進行對比,如圖3 所示,可見粗糙網(wǎng)格雖可求解出更小流量點,但整體特性上級總壓比和等熵效率低于試驗數(shù)據(jù);隨網(wǎng)格量增大,特性線逐漸向右上角移動,折合流量、級總壓比和等熵效率的數(shù)值更大,與試驗值吻合較好;同時中等網(wǎng)格和精細網(wǎng)格求解的失速裕度及特性差距很小,沒有明顯變化。從求解精度和計算資源角度綜合考慮,最終選擇了中等網(wǎng)格量,即單通道計算域的網(wǎng)格量約為247 萬,全環(huán)計算域的網(wǎng)格量約為5 900 萬。

    表2 網(wǎng)格流線、周向、徑向分布Table 2 Mesh distributions at streamline, circumferential and radial directions

    圖3 不同網(wǎng)格下壓氣機特性與試驗數(shù)據(jù)[27]對比Fig.3 Comparison of compressor characteristics for different meshes with test data[27]

    3 均勻進氣結(jié)果與分析

    3.1 數(shù)值結(jié)果驗證

    驗證過程是將數(shù)值模擬結(jié)果和試驗數(shù)據(jù)進行對比,量化數(shù)值仿真的精度,也是進行研究的第一步。65%和100%轉(zhuǎn)速壓氣機級總壓比和等熵效率與試驗數(shù)據(jù)的對比如圖4 所示,可看出數(shù)值模擬得到的特性和試驗值吻合較好,同時能準確捕捉不同轉(zhuǎn)速下壓氣機的近失速點和堵塞點。

    圖4 數(shù)值結(jié)果與試驗數(shù)據(jù)[27]對比Fig.4 Comparison of numerical results with test data[27]

    為進一步驗證數(shù)值模擬方法的準確性,圖5 給出了100%轉(zhuǎn)速下最高效率(Peak Efficiency,PE)點和近失速工況(Near Stall,NS)點的級出口面參數(shù)徑向分布對比,圖中mC為折合流量。不同工況點級出口面的總壓比和總溫比徑向分布對比如圖5(a)和圖5(b)所示,模擬值的總壓比在葉根和葉尖處相較于試驗數(shù)據(jù)略大,而總溫比整體吻合較好。從PE 點到NS 點級總壓比徑向分布變得不均勻,NS 點靠近葉根和葉頂區(qū)域,均有較大的總壓分布;而級總溫比在60%葉高以上溫梯度大且變化范圍較大。盡管數(shù)值模擬結(jié)果與試驗數(shù)據(jù)在部分葉高有一定偏差,但結(jié)果與其他研究人員采用不同求解器的計算結(jié)果[27,33]趨勢基本一致。

    圖5 級出口面(ME30)參數(shù)與試驗數(shù)據(jù)[27]對比Fig.5 Comparison of parameters at stage exit (ME30)with test data[27]

    試驗臺的靜子輪轂空腔會存在泄漏流,根據(jù)文獻[34],考慮0.5%的泄漏流量后數(shù)值模擬結(jié)果和試驗數(shù)據(jù)在徑向分布上的對比會更好,工作中未考慮泄漏流的影響。

    最后從內(nèi)部流場角度對3 個不同的典型工況(近失速點、最高效率點和近堵塞點)的靜壓參數(shù)進行對比。數(shù)值模擬結(jié)果提取的99.9%葉高靜壓云圖與試驗轉(zhuǎn)子葉頂區(qū)的壁面靜壓云圖[27]對比如圖6 所示,圖中LE 為前緣(Leading Edge),TE 為尾緣(Trailing Edge),可見數(shù)值模擬能精確捕捉脫體激波、葉頂泄漏流的軌跡和通道激波的位置。圖3~圖6 的對比充分驗證了數(shù)值仿真設(shè)置的正確性。

    圖6 數(shù)值模擬所得內(nèi)部流場與試驗數(shù)據(jù)[27]的對比Fig.6 Comparison of internal flow fields obtained from numerical simulation with test data[27]

    3.2 均勻進氣失速特征

    文獻[35]的數(shù)值研究表明誘發(fā)該壓氣機失速的先兆波為突尖型,計算模擬捕捉到前緣溢出和尾緣回流的流動特征。在失速模擬中從失速前最后一個穩(wěn)定解開始節(jié)流。失速特征如下:在約1 rev 后壓比上升,質(zhì)量流量開始下降;在3 rev 左右達最小值后質(zhì)量流量和壓比又開始增加,在約6 rev 后達最終失速狀態(tài)。通過失速模擬中攻角變化的研究確定了該壓氣機失速擾動發(fā)生在轉(zhuǎn)子葉頂區(qū)域,且在失速先兆階段動葉周圍的流動發(fā)生前緣溢出和強烈的尾緣回流現(xiàn)象,為突尖型失速先兆。模擬過程中失速團的周向傳播速度約為1/3 轉(zhuǎn)子轉(zhuǎn)速。文獻[36]的試驗研究表明該壓氣機以突尖型先兆啟動失速,形成一個單獨的失速團并繞環(huán)旋轉(zhuǎn)。試驗中,100%轉(zhuǎn)速下從穩(wěn)定極限點關(guān)小閥門。在約5 rev 時失速團尺寸最大,對應(yīng)局部流量最小,傳播速度約為0.44 倍轉(zhuǎn)子轉(zhuǎn)速;隨后失速團尺寸開始減小,約15 rev 時失速團尺寸達最小,對應(yīng)局部流量最大,傳播速度約為0.52 倍轉(zhuǎn)子轉(zhuǎn)速。試驗觀測到隨流量減小失速團尺寸增大、傳播速度變慢的特性,失速團的尺寸和傳播速度隨壓氣機接近穩(wěn)定的失速點而震蕩。

    4 總溫畸變下壓氣機失速分析

    4.1 主要參數(shù)變化

    為分析進氣總溫畸變影響壓氣機性能的物理機理,表3 展示了PE 點(折合流量為14.58 kg/s)壓氣機進口面高溫畸變區(qū)和低溫非畸變區(qū)關(guān)鍵參數(shù)的對比,可見高溫畸變區(qū)總溫為500 K,低溫非畸變區(qū)總溫為295 K;總壓均為99 kPa;高溫區(qū)氣流密度為0.656 8 kg·m-3,低溫區(qū)氣流密度1.109 4 kg·m-3;高溫區(qū)軸向速度為141 m·s-1,而低溫區(qū)軸向速度為109 m·s-1;高溫區(qū)密流為93 kg·m-2·s-1,低溫區(qū)密流為121 kg·m-2·s-1??梢娍倻鼗儗簹鈾C的主要影響為高溫區(qū)氣流密度大幅減小,軸向速度增大。同時高溫畸變區(qū)工作在低折合轉(zhuǎn)速線上,因此總溫畸變下壓氣機工作特性點往左下角移動。

    表3 高、低溫區(qū)關(guān)鍵參數(shù)對比Table 3 Comparison of key parameters in high and low temperature regions

    4.2 不同周向位置總壓徑向分布

    為總結(jié)總溫畸變下壓氣機內(nèi)部流場隨不同工況點的變化規(guī)律,對不同周向位置處時均總壓徑向分布情況進行分析。在總溫畸變非定常計算的PE 點(折合流量為14.58 kg/s)和NS 點(折合流量為13.46 kg/s)的轉(zhuǎn)子出口面,每隔45°周向角度提取總壓徑向數(shù)據(jù)進行對比。周向角度位置定義如圖7 所示,背景為NS 點轉(zhuǎn)子出口面時均總壓云圖。

    圖7 轉(zhuǎn)子出口面周向角度位置Fig.7 Position of circumferential angle at rotor exit surface

    圖8 (a)、圖8(b)分別展示了總溫畸變下PE工況點0°~180°和180°~360°不同周向位置的總壓徑向分布情況,可見在0°、45°、90°、135°和180°的高溫區(qū)周向位置,總壓徑向分布逐漸變窄。不同周向位置的總壓徑向分布形式比較類似,葉中區(qū)域分布均勻,而葉頂葉根區(qū)域由于端壁邊界層、角區(qū)等導(dǎo)致流動損失進而使總壓存在較大梯度;在180°、225°、270°、315°和360°的周向位置,總壓徑向分布逐漸變寬。

    圖8 PE 點轉(zhuǎn)子出口面不同周向位置處總壓徑向分布Fig.8 Radial distributions of total pressure at different circumferential positions of rotor exit for PE condition

    圖9(a)、圖9(b)分別展示了總溫畸變下NS工況點0°~180°和180°~360°不同周向位置的總壓徑向分布情況,可見在0°、45°、90°、135°和180°的高溫區(qū)周向位置,總壓徑向分布呈現(xiàn)逐漸變窄的趨勢。在180°、225°、270°、315°和360°的低溫區(qū)周向位置,總壓徑向分布呈現(xiàn)逐漸變寬的趨勢且分布較為不均勻。轉(zhuǎn)子出口面總壓云圖如圖7所示,在225°周向位置處,葉根區(qū)域的總壓較小,而葉頂區(qū)域存在一個高總壓區(qū),總壓徑向分布不均勻。這是由于畸變引起氣流角沿周向分布不同,且沿徑向也發(fā)生遷移,從而加功量不同,導(dǎo)致總壓在不同周向位置的徑向分布差異。

    圖9 NS 點轉(zhuǎn)子出口面不同周向位置處總壓徑向分布Fig.9 Radial distribution of total pressure at different circumferential positions of rotor exit for NS condition

    圖10 為不同葉高處和具有相同噴嘴喉部面積的、均勻來流的、絕對周向氣流角的周向分布對比。定義逆轉(zhuǎn)子旋轉(zhuǎn)方向為絕對周向氣流角正方向。結(jié)合圖11 中的速度三角形,可見在0°~180°范圍絕對周向氣流角α由正值減小到負值,相對氣流角β逐漸減小,因此攻角減小,葉片載荷減小,總壓減小。在180°~360°范圍,絕對周向氣流角α由負值增大到正值,相對氣流角β逐漸增大,因此攻角增大,葉片載荷增大,總壓增大。并且對于該壓氣機而言,總溫畸變下周向氣流角變化范圍在葉根區(qū)域最大,從10%葉高到90%葉高逐步減小。以上研究結(jié)果表明,在轉(zhuǎn)子由低溫非畸變區(qū)轉(zhuǎn)入高溫畸變區(qū)時攻角最大,因此猜測此處可能最先發(fā)生失穩(wěn)。為驗證這一猜想,接下來進一步對失速過程進行了分析。

    圖10 轉(zhuǎn)子進口絕對周向氣流角Fig.10 Absolute circumferential flow angle at rotor inlet

    圖11 速度三角形中絕對周向氣流角Fig.11 Absolute circumferential flow angle in velocity triangle

    4.3 總溫畸變數(shù)值結(jié)果與分析

    進氣總溫畸變條件下,壓氣機的級總壓比特性曲線如圖12 所示。從圖12 中可見由一個穩(wěn)定解過渡到失速狀態(tài),總壓比和折合流量先略微升高,然后迅速下降,再回升,以此循環(huán)往復(fù)呈現(xiàn)周期性變化特征。

    圖12 總溫畸變下壓氣機特性Fig.12 Compressor performance characteristics with total temperature distortion

    進氣總溫畸變條件下失速過程中壓氣機進口物理流量隨轉(zhuǎn)數(shù)的變化情況如圖13 所示,物理流量和折合流量的換算關(guān)系為

    圖13 總溫畸變下失速過程壓氣機進口物理流量Fig.13 Physical flow rate at compressor inlet during stall process under total temperature distortion

    式中:mC為折合流量;m為物理流量為級進口總壓;為級進口總溫。

    從10 rev 開始由一個穩(wěn)定的工況解輕微減小噴嘴喉部面積。由圖13 可見從11 rev 開始物理流量迅速從10.8 kg/s 下降到7.6 kg/s,隨后升高到11.3 kg/s,往復(fù)變化,且物理流量的震蕩幅度有輕微減小,變化趨勢與圖12 中壓氣機特性變化相符合。

    4.4 失速過程軸向速度分布

    圖14 展示了不同轉(zhuǎn)數(shù)t下轉(zhuǎn)子出口面的軸向速度分布云圖,可見在11 rev 軸向速度沿周向分布較均勻,未出現(xiàn)回流;在13 rev 出現(xiàn)2 個較大失速團和較小失速團(黑虛框),且低軸向速度區(qū)域近乎占據(jù)了整個轉(zhuǎn)子出口面一半周向范圍,此時壓氣機已經(jīng)處于失速狀態(tài);轉(zhuǎn)到16 rev 和20 rev時失速團合并為1 個,且占據(jù)的周向范圍大幅縮小;轉(zhuǎn)到23 rev 和27 rev 時失速團變?yōu)?~2 個較大失速團和若干小失速團,大約占據(jù)整個轉(zhuǎn)子出口面1/3 周向范圍。結(jié)合圖13 物理流量隨轉(zhuǎn)數(shù)時間的震蕩變化情況看,當物理流量較?。ㄈ?3 rev 和27 rev)時失速團的數(shù)量較多、占據(jù)范圍較大,堵塞程度較嚴重;當物理流量較大(如16 rev 和20 rev)時失速團數(shù)量較少、占據(jù)范圍較小,堵塞程度較輕;由此計算出流量振蕩頻率約為75.0 Hz。

    圖14 失速過程的軸向速度云圖Fig.14 Axial velocity contours during stall process

    4.5 失速過程軸向速度為0 的等值面分布

    為進一步分析失速過程中失速團的體積大小和空間分布,圖15 展示了不同轉(zhuǎn)數(shù)瞬態(tài)解的轉(zhuǎn)子域中軸向速度為0 的三維等值面分布云圖。在全環(huán)域瞬態(tài)解中生成軸向速度為0 的等值面,這個等值面可看作失速團的邊界,等值面圍成的空間大小可看作失速團的體積,這種方法是由Zhang 和Vahdati[31]提出的。從圖15 中可清晰看到進氣總溫畸變條件下壓氣機失速過程中失速團的空間分布集中在葉頂區(qū)域。隨著轉(zhuǎn)子旋轉(zhuǎn)時間推進,失速團的空間位置在葉頂沿周向移動,失速團的體積也在不斷變化。

    圖15 失速過程軸向速度為0 的等值面云圖Fig.15 Iso-surface of zero axial velocity contours during stall process

    4.6 靜壓探針信號分析

    為進一步確定失速先兆最先出現(xiàn)的周向位置和先兆波類型,在數(shù)值計算設(shè)置中沿周向均勻布置了12 支絕對靜壓探針監(jiān)測轉(zhuǎn)子前緣的靜壓值變化。探針距離葉片前緣為10%轉(zhuǎn)子軸向弦長,探針的徑向位置為99%葉高,探針布局方式和序號如圖16 所示。

    圖16 靜壓探針的周向布局Fig.16 Circumferential arrangement of static pressure probes

    圖17 展示了12 支絕對靜壓探針的壓力信號隨時間的波動情況,可看出CH10 探針(即壓氣機轉(zhuǎn)子從低總溫非畸變區(qū)轉(zhuǎn)入高總溫畸變區(qū)的周向位置處)最先監(jiān)測到明顯壓力波動,隨后順轉(zhuǎn)子旋轉(zhuǎn)方向的探針壓力波動幅值增大,表明轉(zhuǎn)子葉片前緣的非定常波動增強,壓氣機由穩(wěn)定的近失速工況點轉(zhuǎn)入失速狀態(tài)。失速先兆的周向傳播速度約為88.9%轉(zhuǎn)子轉(zhuǎn)速,失速初期失速團的傳播速度約為66.0%轉(zhuǎn)子轉(zhuǎn)速。圖18(a)和圖18(b)展示了11.0 rev 和11.5 rev 時98%葉高的速度矢量和軸向速度云圖,可見在11.0 rev 時只觀測到葉頂泄漏流從轉(zhuǎn)子前緣輕微溢出;在11.5 rev 時可看到葉頂泄漏流在轉(zhuǎn)子前緣明顯溢出和尾緣回流現(xiàn)象,此時葉片通道發(fā)生堵塞,最終形成失速擾動。這符合Vo 等[37]通過單/多通道壓氣機定常/非定常數(shù)值計算提出的突尖型失速先兆的準則:葉頂泄漏流從轉(zhuǎn)子葉片前緣溢出和尾緣跨通道回流。因此總溫畸變下誘發(fā)旋轉(zhuǎn)失速的先兆波為突尖型失速先兆。

    圖17 周向靜壓探針的壓力信號Fig.17 Pressure signals of circumferential static pressure probes

    圖18 不同時刻速度矢量和軸向速度云圖Fig.18 Velocity vector and axial velocity contour at different times

    5 結(jié) 論

    對德國Darmstadt 跨聲壓氣機進行了均勻進氣條件下單通道定常數(shù)值模擬,數(shù)值結(jié)果與試驗數(shù)據(jù)吻合較好,驗證了計算設(shè)置的準確性;進行了總溫畸變條件下的壓氣機全環(huán)非定常數(shù)值模擬,對失速過程進行了研究,得到以下結(jié)論。

    1) 在周向范圍180°、高溫畸變區(qū)500 K 的周向總溫畸變條件下壓氣機總壓比-流量特性大幅下降,在失速過程中特性曲線呈周期性震蕩現(xiàn)象,當物理流量較小時對應(yīng)壓比較小,失速團數(shù)量較多、占據(jù)范圍較大,堵塞程度較重;當物理流量較大時對應(yīng)壓比較大,失速團數(shù)量較少、占據(jù)范圍較小,堵塞程度較輕。

    2) 總溫畸變下不同周向位置處的總壓徑向分布不同。順轉(zhuǎn)子葉片旋轉(zhuǎn)方向總壓在高溫區(qū)逐漸減小,在低溫區(qū)逐漸增大,且在葉根區(qū)域變化最為明顯。

    3) 總溫畸變下最先發(fā)生壓力擾動的周向位置是轉(zhuǎn)子轉(zhuǎn)入高溫畸變區(qū)處,誘發(fā)旋轉(zhuǎn)失速的先兆波為突尖型失速先兆。

    4) 總溫畸變下失速先兆的周向傳播速度約為88.9%轉(zhuǎn)子轉(zhuǎn)速,失速初期失速團的傳播速度約為66.0%轉(zhuǎn)子轉(zhuǎn)速,高于均勻進氣下失速團的傳播速度。

    猜你喜歡
    總溫總壓周向
    周向拉桿轉(zhuǎn)子瞬態(tài)應(yīng)力分析與啟動曲線優(yōu)化
    總壓探針性能結(jié)構(gòu)敏感性分析
    3m×2m結(jié)冰風洞熱流場品質(zhì)提高及評估
    大型結(jié)冰風洞熱流場符合性驗證
    可調(diào)式總壓耙設(shè)計及應(yīng)用
    亞聲速條件下總壓探針臨壁效應(yīng)的數(shù)值研究
    2 m超聲速風洞流場變速壓控制方法研究
    高超音速氣流總溫傳感器設(shè)計
    周向定位旋轉(zhuǎn)分度鉆模設(shè)計
    一種商用輕型載重汽車輪胎
    久久久久九九精品影院| av天堂在线播放| 在线观看av片永久免费下载| 日韩欧美一区二区三区在线观看| videossex国产| 久久99热6这里只有精品| 国产在线精品亚洲第一网站| 欧美日韩黄片免| 久久亚洲精品不卡| 日韩欧美免费精品| 很黄的视频免费| 亚洲在线自拍视频| 动漫黄色视频在线观看| 精品午夜福利在线看| 在线观看午夜福利视频| 精品日产1卡2卡| 真人做人爱边吃奶动态| 此物有八面人人有两片| 国产黄a三级三级三级人| 禁无遮挡网站| 日韩大尺度精品在线看网址| 成人毛片a级毛片在线播放| 99热网站在线观看| 美女cb高潮喷水在线观看| 亚洲久久久久久中文字幕| 欧美日韩国产亚洲二区| 久久国产乱子免费精品| 欧美zozozo另类| 麻豆精品久久久久久蜜桃| 日日摸夜夜添夜夜添小说| 亚洲人成网站在线播放欧美日韩| 我的女老师完整版在线观看| 一本久久中文字幕| 成人综合一区亚洲| 日本免费一区二区三区高清不卡| 欧美日韩中文字幕国产精品一区二区三区| 99精品久久久久人妻精品| 最近最新中文字幕大全电影3| 看免费成人av毛片| 久久婷婷人人爽人人干人人爱| 狠狠狠狠99中文字幕| 两人在一起打扑克的视频| 真人做人爱边吃奶动态| 国产精品无大码| 听说在线观看完整版免费高清| 一级黄片播放器| 日韩一区二区视频免费看| 色哟哟哟哟哟哟| 亚洲精品色激情综合| a在线观看视频网站| 欧美国产日韩亚洲一区| 亚洲电影在线观看av| 国产91精品成人一区二区三区| 嫩草影院入口| 免费黄网站久久成人精品| 欧美激情久久久久久爽电影| 黄片wwwwww| 日本a在线网址| 男人狂女人下面高潮的视频| 看十八女毛片水多多多| 五月玫瑰六月丁香| 别揉我奶头~嗯~啊~动态视频| 久久久精品大字幕| 日韩精品青青久久久久久| 欧美一区二区国产精品久久精品| 夜夜夜夜夜久久久久| 久久久久免费精品人妻一区二区| eeuss影院久久| 欧美成人性av电影在线观看| 黄色丝袜av网址大全| av福利片在线观看| 免费高清视频大片| 国产日本99.免费观看| 日韩欧美在线乱码| 日韩中字成人| 国产精品久久久久久av不卡| 亚洲人成网站高清观看| 男人舔奶头视频| 噜噜噜噜噜久久久久久91| 韩国av在线不卡| 国产亚洲精品av在线| av中文乱码字幕在线| 日韩欧美免费精品| 两性午夜刺激爽爽歪歪视频在线观看| 午夜福利欧美成人| 春色校园在线视频观看| 日日夜夜操网爽| 成人毛片a级毛片在线播放| 国产精品日韩av在线免费观看| 一进一出抽搐动态| 又黄又爽又免费观看的视频| 美女黄网站色视频| 热99re8久久精品国产| 99精品在免费线老司机午夜| 婷婷精品国产亚洲av在线| 我要看日韩黄色一级片| 亚洲av电影不卡..在线观看| 亚洲第一区二区三区不卡| 欧美黑人欧美精品刺激| aaaaa片日本免费| 国产午夜精品论理片| 五月伊人婷婷丁香| 春色校园在线视频观看| 国产精品野战在线观看| 色在线成人网| 男女视频在线观看网站免费| 久久久精品大字幕| 一进一出抽搐gif免费好疼| 黄色女人牲交| 欧美一区二区国产精品久久精品| 成人无遮挡网站| 免费在线观看成人毛片| 欧美成人免费av一区二区三区| 亚洲三级黄色毛片| 1000部很黄的大片| 国产精品三级大全| 精品午夜福利视频在线观看一区| 日本黄色视频三级网站网址| 少妇裸体淫交视频免费看高清| 亚洲欧美精品综合久久99| 成年版毛片免费区| 少妇的逼好多水| bbb黄色大片| av在线亚洲专区| av国产免费在线观看| 亚洲黑人精品在线| 日韩一区二区视频免费看| 国产男靠女视频免费网站| 一夜夜www| 尤物成人国产欧美一区二区三区| 国内少妇人妻偷人精品xxx网站| 亚洲无线在线观看| 直男gayav资源| 久久精品91蜜桃| 亚洲内射少妇av| 色视频www国产| 黄色一级大片看看| 亚洲av免费在线观看| 简卡轻食公司| 又黄又爽又免费观看的视频| 我的女老师完整版在线观看| 成年免费大片在线观看| 亚洲中文字幕日韩| 18禁在线播放成人免费| 男女那种视频在线观看| av在线观看视频网站免费| 日本爱情动作片www.在线观看 | 女的被弄到高潮叫床怎么办 | 人妻少妇偷人精品九色| 婷婷亚洲欧美| 天堂动漫精品| 久久草成人影院| 黄色日韩在线| 黄色视频,在线免费观看| 人妻久久中文字幕网| 久久久久久久久久成人| 久久这里只有精品中国| 无人区码免费观看不卡| 又紧又爽又黄一区二区| 我要看日韩黄色一级片| 男女那种视频在线观看| 美女免费视频网站| 亚洲av免费高清在线观看| 看片在线看免费视频| 国产一区二区三区av在线 | 欧美性感艳星| 最近最新中文字幕大全电影3| 国产爱豆传媒在线观看| 午夜福利在线观看吧| 黄色欧美视频在线观看| 婷婷色综合大香蕉| 桃红色精品国产亚洲av| 久久精品国产亚洲av涩爱 | 国产一区二区三区av在线 | 久久6这里有精品| 国产aⅴ精品一区二区三区波| 精品人妻偷拍中文字幕| 午夜福利在线观看吧| 国产激情偷乱视频一区二区| 国产高清激情床上av| 亚洲国产色片| 国产精品一区二区三区四区免费观看 | 日本-黄色视频高清免费观看| 亚洲最大成人av| 欧美另类亚洲清纯唯美| 久久精品人妻少妇| 欧美成人一区二区免费高清观看| 可以在线观看毛片的网站| 亚洲欧美日韩无卡精品| 久久久国产成人免费| 亚洲va日本ⅴa欧美va伊人久久| 香蕉av资源在线| 18禁黄网站禁片午夜丰满| 天天一区二区日本电影三级| 成人特级av手机在线观看| 欧美成人免费av一区二区三区| 成年女人看的毛片在线观看| 日日干狠狠操夜夜爽| 真实男女啪啪啪动态图| 亚洲成人中文字幕在线播放| 能在线免费观看的黄片| av福利片在线观看| 搡老妇女老女人老熟妇| 久久久久久久久中文| 久久精品久久久久久噜噜老黄 | 亚洲美女黄片视频| 久99久视频精品免费| 免费无遮挡裸体视频| 国产一区二区激情短视频| 欧美成人免费av一区二区三区| 精品一区二区三区人妻视频| 欧美黑人巨大hd| 国产av麻豆久久久久久久| 精华霜和精华液先用哪个| 国产精华一区二区三区| 成人国产综合亚洲| 舔av片在线| 国产av一区在线观看免费| 欧美又色又爽又黄视频| av女优亚洲男人天堂| 久久久久国产精品人妻aⅴ院| 日本精品一区二区三区蜜桃| 亚洲av电影不卡..在线观看| 最好的美女福利视频网| 午夜福利在线在线| 久9热在线精品视频| 成人国产综合亚洲| 亚洲av一区综合| 淫妇啪啪啪对白视频| 国产一区二区三区av在线 | 中国美女看黄片| netflix在线观看网站| 精品乱码久久久久久99久播| 一卡2卡三卡四卡精品乱码亚洲| 国产高潮美女av| 欧美黑人巨大hd| 日韩在线高清观看一区二区三区 | 婷婷丁香在线五月| 我的老师免费观看完整版| 岛国在线免费视频观看| 日韩欧美在线乱码| 麻豆国产97在线/欧美| 久久精品国产亚洲av天美| 亚洲专区中文字幕在线| 真人一进一出gif抽搐免费| 色综合色国产| 一级毛片久久久久久久久女| 日本一本二区三区精品| 婷婷六月久久综合丁香| 中出人妻视频一区二区| 自拍偷自拍亚洲精品老妇| 女人十人毛片免费观看3o分钟| 亚洲av电影不卡..在线观看| 麻豆一二三区av精品| 最好的美女福利视频网| 欧美性猛交╳xxx乱大交人| 亚洲欧美激情综合另类| 小蜜桃在线观看免费完整版高清| 搞女人的毛片| 黄片wwwwww| 91久久精品国产一区二区成人| 国产精品久久电影中文字幕| 午夜影院日韩av| 亚洲av免费高清在线观看| 欧美精品啪啪一区二区三区| 亚洲自偷自拍三级| 黄色丝袜av网址大全| 国产 一区精品| 别揉我奶头~嗯~啊~动态视频| 亚洲精品一卡2卡三卡4卡5卡| 国内毛片毛片毛片毛片毛片| 国产精品不卡视频一区二区| 国产男人的电影天堂91| 亚洲天堂国产精品一区在线| 成人国产麻豆网| 国产日本99.免费观看| 国产精品一区二区性色av| 亚洲精品影视一区二区三区av| 在线播放无遮挡| 欧美精品国产亚洲| 久久人妻av系列| www日本黄色视频网| 18禁黄网站禁片午夜丰满| 一区二区三区激情视频| 高清毛片免费观看视频网站| 看免费成人av毛片| 嫩草影院新地址| 国产在线男女| 国产一区二区亚洲精品在线观看| 亚洲欧美日韩高清在线视频| 国产91精品成人一区二区三区| 蜜桃亚洲精品一区二区三区| 精品人妻视频免费看| 能在线免费观看的黄片| 一进一出抽搐gif免费好疼| 久久久久九九精品影院| 九九在线视频观看精品| 色精品久久人妻99蜜桃| 亚洲精品一卡2卡三卡4卡5卡| 搡老妇女老女人老熟妇| 欧美激情国产日韩精品一区| 我要看日韩黄色一级片| 91在线观看av| 露出奶头的视频| 日日夜夜操网爽| 亚州av有码| 高清在线国产一区| 国产精品av视频在线免费观看| 美女高潮的动态| 国产精品久久久久久亚洲av鲁大| 乱系列少妇在线播放| 久久香蕉精品热| 亚洲欧美日韩卡通动漫| 国产精品一区二区三区四区免费观看 | 在线观看美女被高潮喷水网站| 精品国内亚洲2022精品成人| 国产高清视频在线观看网站| av黄色大香蕉| 国产精品免费一区二区三区在线| 久久国产精品人妻蜜桃| 国产主播在线观看一区二区| 欧洲精品卡2卡3卡4卡5卡区| 国产伦精品一区二区三区视频9| 如何舔出高潮| 在线观看一区二区三区| 色av中文字幕| 精品久久久久久久久亚洲 | 精品一区二区免费观看| 特大巨黑吊av在线直播| 美女免费视频网站| 成人亚洲精品av一区二区| 欧美日本亚洲视频在线播放| 亚洲一区高清亚洲精品| 国产 一区精品| 美女高潮的动态| av在线观看视频网站免费| 51国产日韩欧美| 国内少妇人妻偷人精品xxx网站| 一本一本综合久久| 成人av一区二区三区在线看| 在线免费十八禁| 国内精品美女久久久久久| 春色校园在线视频观看| 在线观看免费视频日本深夜| 男插女下体视频免费在线播放| 小说图片视频综合网站| 春色校园在线视频观看| 亚洲一区二区三区色噜噜| 午夜福利视频1000在线观看| 成年女人看的毛片在线观看| 欧美色视频一区免费| 国产精品久久久久久久电影| 亚洲欧美清纯卡通| 三级国产精品欧美在线观看| 成年人黄色毛片网站| 小说图片视频综合网站| 中国美女看黄片| .国产精品久久| 99国产精品一区二区蜜桃av| 亚洲av五月六月丁香网| 中国美女看黄片| 在线a可以看的网站| 毛片一级片免费看久久久久 | 在线观看美女被高潮喷水网站| 中国美白少妇内射xxxbb| 欧美性猛交黑人性爽| av专区在线播放| 天堂√8在线中文| 国产精品一区二区免费欧美| 国产高清视频在线观看网站| 人妻丰满熟妇av一区二区三区| 亚洲国产欧洲综合997久久,| 国产激情偷乱视频一区二区| 听说在线观看完整版免费高清| 欧美日本视频| 91在线精品国自产拍蜜月| 欧美日韩瑟瑟在线播放| 免费搜索国产男女视频| 国产精品1区2区在线观看.| 亚洲五月天丁香| 精品福利观看| 成人无遮挡网站| 熟女电影av网| 午夜精品在线福利| 不卡一级毛片| 一区二区三区免费毛片| 中文资源天堂在线| av国产免费在线观看| 久久午夜亚洲精品久久| 在线观看免费视频日本深夜| 亚洲四区av| 成年女人看的毛片在线观看| 真人一进一出gif抽搐免费| 国产女主播在线喷水免费视频网站 | aaaaa片日本免费| 欧美日本亚洲视频在线播放| 蜜桃亚洲精品一区二区三区| 国产一区二区在线观看日韩| 日本三级黄在线观看| 亚洲精品在线观看二区| aaaaa片日本免费| 国产精品无大码| 欧美成人免费av一区二区三区| 一a级毛片在线观看| 精品99又大又爽又粗少妇毛片 | 真人做人爱边吃奶动态| 一级毛片久久久久久久久女| 老熟妇乱子伦视频在线观看| 一本久久中文字幕| 午夜福利高清视频| 欧美色欧美亚洲另类二区| 婷婷精品国产亚洲av在线| 欧美丝袜亚洲另类 | 干丝袜人妻中文字幕| 日韩欧美精品v在线| 色av中文字幕| 国产精品,欧美在线| 日韩大尺度精品在线看网址| 亚洲最大成人av| 欧美人与善性xxx| 成人午夜高清在线视频| 又紧又爽又黄一区二区| 麻豆成人午夜福利视频| av在线观看视频网站免费| 成年女人毛片免费观看观看9| 男人舔女人下体高潮全视频| 久久久久精品国产欧美久久久| 国产91精品成人一区二区三区| 欧美一级a爱片免费观看看| 女生性感内裤真人,穿戴方法视频| 国产精品人妻久久久久久| 欧美日韩乱码在线| 亚洲专区国产一区二区| 毛片女人毛片| 国产爱豆传媒在线观看| 国内精品久久久久久久电影| 999久久久精品免费观看国产| 91在线精品国自产拍蜜月| 18+在线观看网站| 香蕉av资源在线| 色尼玛亚洲综合影院| 日韩中字成人| 99热6这里只有精品| 小说图片视频综合网站| а√天堂www在线а√下载| 99久久精品一区二区三区| 精品久久久久久久久亚洲 | av在线蜜桃| 日韩高清综合在线| 成年版毛片免费区| 我要搜黄色片| 波多野结衣巨乳人妻| 国产精品人妻久久久影院| 一卡2卡三卡四卡精品乱码亚洲| 久久久久久九九精品二区国产| www.www免费av| 亚洲国产精品合色在线| 欧美性猛交╳xxx乱大交人| 国产激情偷乱视频一区二区| 亚洲18禁久久av| 国产精华一区二区三区| av在线蜜桃| 日本撒尿小便嘘嘘汇集6| 美女大奶头视频| 国产69精品久久久久777片| 成人国产麻豆网| 亚洲在线自拍视频| 欧美性感艳星| 日本成人三级电影网站| 国产真实乱freesex| 国产三级在线视频| 黄色丝袜av网址大全| 哪里可以看免费的av片| 内射极品少妇av片p| 一边摸一边抽搐一进一小说| 一进一出抽搐gif免费好疼| 麻豆av噜噜一区二区三区| 国产高清激情床上av| 啦啦啦观看免费观看视频高清| 亚洲性夜色夜夜综合| 久久久久国产精品人妻aⅴ院| 桃红色精品国产亚洲av| av天堂在线播放| 免费黄网站久久成人精品| 天堂动漫精品| 久久久久免费精品人妻一区二区| 中国美女看黄片| 日韩人妻高清精品专区| 国产v大片淫在线免费观看| 色吧在线观看| 日韩在线高清观看一区二区三区 | 尾随美女入室| 亚洲精品影视一区二区三区av| 1000部很黄的大片| 中出人妻视频一区二区| 国产成年人精品一区二区| 国产精品不卡视频一区二区| 欧美成人a在线观看| 国产女主播在线喷水免费视频网站 | 女人被狂操c到高潮| 国产精品av视频在线免费观看| 99热精品在线国产| 日本a在线网址| 久久香蕉精品热| x7x7x7水蜜桃| 热99在线观看视频| 99热精品在线国产| 亚洲欧美日韩卡通动漫| 国产女主播在线喷水免费视频网站 | 国产精品一区二区性色av| 色视频www国产| 久久久久久久久大av| 国产欧美日韩精品一区二区| 亚洲欧美日韩卡通动漫| 国产女主播在线喷水免费视频网站 | 99久久成人亚洲精品观看| 少妇的逼好多水| 国产精品乱码一区二三区的特点| 亚洲黑人精品在线| 成年女人看的毛片在线观看| 性欧美人与动物交配| 免费不卡的大黄色大毛片视频在线观看 | 又黄又爽又免费观看的视频| 最近在线观看免费完整版| 午夜福利18| 天堂影院成人在线观看| 听说在线观看完整版免费高清| 又紧又爽又黄一区二区| 日本撒尿小便嘘嘘汇集6| 亚洲av.av天堂| 日本一本二区三区精品| 哪里可以看免费的av片| 中文资源天堂在线| 午夜免费激情av| 亚洲精品亚洲一区二区| 麻豆精品久久久久久蜜桃| 午夜免费成人在线视频| 精品久久久久久久久久久久久| 久久精品综合一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 狂野欧美白嫩少妇大欣赏| 久久这里只有精品中国| 日韩欧美在线乱码| 性色avwww在线观看| 最好的美女福利视频网| 国产高清三级在线| 搡老妇女老女人老熟妇| 成人av一区二区三区在线看| 亚洲人成网站在线播| 成年女人永久免费观看视频| 白带黄色成豆腐渣| 在线免费观看不下载黄p国产 | 99精品在免费线老司机午夜| 亚洲国产精品成人综合色| 国产乱人视频| 精品久久久久久久久久久久久| 我要看日韩黄色一级片| 亚洲国产高清在线一区二区三| 国产视频内射| 美女黄网站色视频| 伊人久久精品亚洲午夜| 国产欧美日韩一区二区精品| 最新中文字幕久久久久| 久久久久久久久中文| 国产成人影院久久av| 亚洲欧美日韩无卡精品| 99热精品在线国产| 日本 欧美在线| 日本五十路高清| 又黄又爽又刺激的免费视频.| 搡老妇女老女人老熟妇| ponron亚洲| 欧洲精品卡2卡3卡4卡5卡区| 久久久精品欧美日韩精品| 欧美区成人在线视频| 干丝袜人妻中文字幕| 哪里可以看免费的av片| 婷婷精品国产亚洲av| 欧美色欧美亚洲另类二区| 亚洲自偷自拍三级| 99久久九九国产精品国产免费| 午夜激情欧美在线| 亚洲av成人av| 欧美不卡视频在线免费观看| 很黄的视频免费| 88av欧美| 欧美日韩瑟瑟在线播放| 狂野欧美白嫩少妇大欣赏| 高清日韩中文字幕在线| 免费看光身美女| 日韩 亚洲 欧美在线| 国产免费av片在线观看野外av| 91狼人影院| 国产淫片久久久久久久久| 久久99热这里只有精品18| 亚洲,欧美,日韩| 别揉我奶头~嗯~啊~动态视频| 黄色女人牲交| 国产精品99久久久久久久久| 欧美色视频一区免费| 久久精品国产亚洲av香蕉五月| 色5月婷婷丁香| 午夜精品久久久久久毛片777| 国产蜜桃级精品一区二区三区| 伦精品一区二区三区| 人妻少妇偷人精品九色| 少妇丰满av| 91在线观看av| 欧美色视频一区免费| www.色视频.com| 欧美+日韩+精品| 午夜精品久久久久久毛片777| 他把我摸到了高潮在线观看| 国产探花极品一区二区| 午夜爱爱视频在线播放|