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

    載荷分布對可控擴散葉型性能的影響

    2020-12-04 05:06:30李夢雪周正貴
    科學技術(shù)與工程 2020年30期
    關(guān)鍵詞:附面層葉型馬赫數(shù)

    李夢雪,周正貴

    (南京航空航天大學能源與動力學院,南京 210016)

    可控擴散葉型(CDA)自20世紀70年代末發(fā)展至今,由于能很好地可控吸力面氣流擴散、消除或減弱激波強度,避免或推遲附面層分離,因而損失更小、低損失工作范圍更寬,在壓氣機上得到廣泛應(yīng)用[1-6]。

    中外許多學者對CDA葉型的正、反問題設(shè)計方法以及優(yōu)化設(shè)計方法進行了研究[7-12]。王琦等[13]采用正問題設(shè)計方法,在葉型參數(shù)化方面使用雙圓弧中弧線、多段圓弧厚度分布以及橢圓前緣,設(shè)計出了更不易發(fā)生流動分離、具有更高壓比及效率的CDA葉型。張小玉等[14]發(fā)展了基于目標壓力分布的壓氣機葉型黏性反問題設(shè)計方法,實現(xiàn)了對某CDA動葉包含前、尾緣的吸力面和壓力面整體葉型的反問題修正。Dunker等[15]在勢、流函數(shù)法基礎(chǔ)上開發(fā)了多葉高截面可控擴散葉型的設(shè)計方法,用該方法設(shè)計得到的壓氣機級在滿足非設(shè)計點性能要求的情況下葉片數(shù)比原始的NACA 65靜子級減少了30%。Sanger[16]于1982年首次將優(yōu)化設(shè)計方法運用于CDA設(shè)計中。趙鵬程等[17]對某大折轉(zhuǎn)角靜子葉柵在設(shè)計狀態(tài)下進行優(yōu)化設(shè)計得到的葉型吸力面馬赫數(shù)(Ma)在峰值后平緩擴散至尾緣,局部激波消除,轉(zhuǎn)捩位置由初始的40%弦長位置前移至20%弦長位置,實驗結(jié)果顯示優(yōu)化葉型設(shè)計點損失大幅下降且具有更寬廣的低損失攻角范圍。劉波等[18-19]采用反問題方法設(shè)計得到一個符合可控擴散規(guī)律的初始葉型,隨后以葉柵總壓損失系數(shù)為目標進行設(shè)計狀態(tài)下的優(yōu)化計算,并采用正問題方法對非設(shè)計點進行氣動驗算,通過人工調(diào)整葉型參數(shù)達到全局優(yōu)化。對于高亞音壓氣機葉柵流動,以設(shè)計點性能為目標進行優(yōu)化設(shè)計通??色@得符合可控擴散規(guī)律葉型,但非設(shè)計點性能不一定滿足設(shè)計要求。因此,此后研究人員發(fā)展了考慮非設(shè)計點性能或全工況優(yōu)化方法。

    K?ller等[20]兼顧設(shè)計點與非設(shè)計點性能進行葉型優(yōu)化設(shè)計,與滿足相同設(shè)計要求的CDA葉型進行對比,優(yōu)化設(shè)計葉型最大彎度更靠前、葉型中部彎度減小,葉型載荷更為靠前,吸力面峰值馬赫數(shù)略靠前;全攻角范圍損失減小、喘振裕度至少增大1.5°。徐全勇等[21]對CDA葉型進行全攻角優(yōu)化,優(yōu)化葉型吸力面峰值馬赫數(shù)位置由40%弦長附近前移至20%弦長處,優(yōu)化結(jié)果發(fā)現(xiàn)載荷前移使得擴壓更為平緩、附面層抗分離能力增強。

    吸力面峰值馬赫數(shù)小于1.3、峰值馬赫數(shù)位置靠前是可控擴散葉型重要載荷特征。但未見研究人員針對此重要特征進行詳細研究,分析峰值點馬赫數(shù)及其位置的影響規(guī)律,為可控擴散葉型設(shè)計提供有益參考。為此,針對高亞音轉(zhuǎn)子和靜子葉型,以給定吸力面峰值點馬赫數(shù)及其位置且載荷分布符合可控擴散規(guī)律為目標函數(shù),構(gòu)建自動優(yōu)化葉型設(shè)計方法;獲得不同載荷分布規(guī)律的可控擴散葉型,探究吸力面峰值馬赫數(shù)及其位置對葉柵全工況氣動性能的影響。

    1 優(yōu)化設(shè)計方法

    葉型設(shè)計采用基于遺傳算法和正問題流場計算結(jié)合的自動優(yōu)化方法,優(yōu)化軟件為自主研發(fā),包括遺傳算法優(yōu)化模塊、葉型參數(shù)化模塊、正問題流場數(shù)值計算模塊、目標函數(shù)計算模塊。以下介紹目標函數(shù)設(shè)置和葉型參數(shù)化方法,遺傳算法和正問題流場數(shù)值計算方法可見參考文獻[22]。二維回轉(zhuǎn)面葉柵流場計算采用空間中心差分和四步Runge-Kutta法時間推進進行定常求解,紊流模型采用Baldwin-Lomax雙層代數(shù)模型。為了檢驗流場數(shù)值計算方法的可靠性,進行進口馬赫數(shù)為0.6的葉柵流場計算,與實驗值進行對比[23]。從圖1可以看出,計算與實測設(shè)計點葉片表面等熵馬赫數(shù)分布及損失氣流轉(zhuǎn)角隨攻角變化特性吻合較好。

    圖1 數(shù)值計算與實驗結(jié)果對比Fig.1 Comparison of computational results and experimental results

    1.1 目標函數(shù)

    可控擴散葉型的等熵馬赫數(shù)分布具有吸力面氣流連續(xù)加速至峰值馬赫數(shù)隨后無激波連續(xù)減速至葉型尾緣、減速曲線稍向上凹,壓力面等熵馬赫數(shù)近于不變的流動特點。

    如圖2所示,為了采用優(yōu)化方法實現(xiàn)上述葉片表面馬赫數(shù)分布規(guī)律,吸力面等熵馬赫數(shù)分布用兩段曲線構(gòu)成:第一段加速段由給定葉柵進口馬赫數(shù)(Main)、吸力面峰值馬赫數(shù)(Mamax)及其位置(x1)確定的直線;第二段減速段由吸力面峰值馬赫數(shù)及其位置、出口馬赫數(shù)(Maout)及出口處斜率為0確定的二次曲線;壓力面在0.1、0.4、0.8倍相對軸向弦長位置等熵馬赫數(shù)差值為極小,以實現(xiàn)壓力面等熵馬赫數(shù)近于不變。由等熵馬赫數(shù)及進口總壓可獲得對應(yīng)的葉片表面靜壓分布。由此可得吸力面目標壓力分布關(guān)系Pi_obj。

    圖2 等熵馬赫數(shù)分布目標示意圖Fig.2 Schematic diagram of isentropic Mach number distribution objective function

    根據(jù)以上思路,優(yōu)化目標設(shè)置為

    (1)

    式(1)中:F為目標函數(shù)值;?為總壓損失系數(shù);Pi、Pi_obj為葉片吸力面第i個網(wǎng)格節(jié)點的計算和目標靜壓;Δβ、Δβobj分別為計算及目標氣流轉(zhuǎn)角;BLx為軸向相對弦長,Δx=xi+1-xi;Ppj為壓力面壓力,其中j=1、2、3,分別對應(yīng)軸向相對弦長為0.1、0.4、0.8位置處的壓力面壓力。權(quán)重系數(shù)初始值設(shè)置為C1=100,C2=10,C3=20,C4=5。再根據(jù)優(yōu)化結(jié)果,考慮是否調(diào)整權(quán)重系數(shù)重新優(yōu)化。如優(yōu)化葉型的損失過大,則增加C4數(shù)值。以上目標函數(shù)是要達到4個目標:①達到給定吸力面壓力(等熵馬赫數(shù))分布;②壓力面壓力(等熵馬赫數(shù))近于不變;③達到給定氣流轉(zhuǎn)角;④流動損失較小,該目標僅為引導優(yōu)化向葉型合理方向進展,因此權(quán)重系數(shù)較小。

    1.2 優(yōu)化設(shè)計變量

    采用基于修改量參數(shù)化方法對葉型進行參數(shù)化,就是在原始葉型型面上疊加修改量,修改量采用貝塞爾曲線描述,并利用葉型當?shù)睾穸冗M行無量綱化。修改位置為無量綱弦線方向,區(qū)間為[0,1],前緣為0,尾緣為1。給出若干個修改位置及對應(yīng)的修改量,可確定出無量綱修改量關(guān)于無量綱弦長的貝塞爾曲線,將修改量疊加到原始葉型上即可得到修改葉型。表1為吸力面和壓力面修改量參數(shù)化設(shè)置。此外葉型弦長和安裝角也作為設(shè)計變量,其變化范圍分別為[-0.1,0.1](與原始葉型弦長比值)和[-7°,7°]。設(shè)計變量范圍根據(jù)實驗室經(jīng)驗數(shù)值確定,所有算例的修改量均未達到極限值、設(shè)計變量修改范圍設(shè)置合理。

    表1 葉型型面修改范圍Table 1 Range variation of blade airfoils

    2 載荷分布規(guī)律對可控擴散靜子葉型性能影響

    與轉(zhuǎn)子葉型比較,靜子葉片安裝角較小、氣流轉(zhuǎn)角較大,兩者流動特征有明顯差別。因此分別對高亞音轉(zhuǎn)子和靜子葉型進行載荷分布規(guī)律研究。

    為對比不同載荷分布規(guī)律對可控擴散靜子葉型性能的影響,針對進口馬赫數(shù)0.7、進氣氣流角42°、氣流轉(zhuǎn)角42°(即軸向出氣)葉柵,采用自動優(yōu)化方法設(shè)計出吸力面峰值馬赫數(shù)1.0、峰值位于0.07、0.1、0.2、0.3倍相對軸向弦長;吸力面峰值位于0.1倍相對軸向弦長、峰值馬赫數(shù)為0.95、1.0、1.08、1.14可控擴散靜子葉型。

    2.1 設(shè)計工況下流場分析

    圖3為優(yōu)化所得不同載荷分布規(guī)律葉型在軸向弦長放縮至1后得到的葉型對比圖。圖3(a)為相同峰值馬赫數(shù)、不同峰值馬赫數(shù)位置的4個葉型,0.07-1.0、0.1-1.0、0.2-1.0、0.3-1.0表示峰值馬赫數(shù)位置為0.07、0.1、0.2、0.3倍相對軸向弦長處、峰值馬赫數(shù)為1.0的葉型;圖3(b)為相同峰值馬赫數(shù)位置、不同峰值馬赫數(shù)的4個葉型,0.1-0.95、0.1-1.08、0.1-1.14表示峰值馬赫數(shù)位置為0.1倍相對軸向弦長處、峰值馬赫數(shù)為0.95、1.08和1.14的葉型。由圖3結(jié)合表2可以看出:峰值馬赫數(shù)不變時,隨著載荷后移,葉型厚度增加;峰值位置不變時,隨著峰值馬赫數(shù)增大,葉型前部凸起程度增加。各葉型最大相對厚度為0.4~0.5、最大厚度相對位置為0.3~0.4;最大相對撓度為0.1~0.14、最大撓度相對位置為0.3~0.45。

    表2 可控擴散靜子葉型幾何參數(shù)Table 2 Geometry parameters of controlled diffusion stator airfoils

    表3為設(shè)計進氣條件下不同載荷分布葉型氣動性能參數(shù),總壓損失系數(shù)計算公式如式(2)所示,總壓損失系數(shù)取值位置隨式(2)進行說明。由表3可知,采用上述優(yōu)化方法所得葉型可保持氣流轉(zhuǎn)角和進口馬赫數(shù)近于不變;吸力面峰值馬赫數(shù)不變,隨著載荷前移(直至0.1倍軸向相對弦長位置),設(shè)計工況損失減??;載荷繼續(xù)前移時損失迅速增大,最低損失大約位于0.1倍軸向相對弦長處;峰值馬赫數(shù)位于0.1倍軸向相對弦長處,存在最低損失對應(yīng)的峰值馬赫數(shù)(約為1.0)。

    表3 設(shè)計進氣條件下氣動參數(shù)對比Table 3 Comparison of aerodynamic parameters under design inlet condition

    圖4為設(shè)計進氣條件下不同載荷分布靜子葉型表面等熵馬赫數(shù)分布,圖5、圖6為設(shè)計進氣條件下不同載荷分布靜子葉型馬赫數(shù)云圖。從圖4~圖6可以看出,采用上述優(yōu)化方法可實現(xiàn)可控擴散葉型設(shè)計、以及吸力面峰值點馬赫數(shù)位置和數(shù)值的可控。由圖4(a),保持吸力面峰值馬赫數(shù)不變、峰值位置后移,壓力面前緣區(qū)域馬赫數(shù)增加,葉型前部載荷變小;由圖4(b),保持吸力面峰值馬赫數(shù)位置不變、增大峰值馬赫數(shù),壓力面馬赫數(shù)近于不變,葉型前部載荷增加、后部載荷減小。從圖4可以看出,由于進氣條件不變、氣流轉(zhuǎn)角不變,各葉型總載荷(吸力面與壓力面等熵馬赫數(shù)包圍的面積)近于不變。

    圖4 設(shè)計進氣條件下等熵馬赫數(shù)分布對比Fig.4 Comparison of isentropic Mach number distribution under design inlet condition

    圖5 設(shè)計進氣條件下不同峰值馬赫數(shù)位置葉型馬赫數(shù)云圖對比Fig.5 Comparison of Mach number contours of different suction side peak Mach number location airfoils

    圖6 設(shè)計進氣條件下不同峰值馬赫數(shù)葉型馬赫數(shù)云圖對比Fig.6 Comparison of Mach number contours of different suction side peak Mach number airfoils

    結(jié)合表4、圖4、圖5可知,吸力面峰值馬赫數(shù)不變,隨著峰值位置前移(直至0.1倍相對軸向弦長處),吸力面減速區(qū)更長、逆壓力梯度減小,附面層增長越緩慢、附面層損失減小;吸力面峰值馬赫數(shù)位置不變、隨著峰值馬赫數(shù)減小(直至馬赫數(shù)1.0),峰值后逆壓力梯度減小、吸力面附面層損失越小。

    為進一步探究不同載荷分布規(guī)律葉型流動損失構(gòu)成,取出各部分流動損失。葉柵通道內(nèi)各部分損失劃分如圖7所示,計算公式為

    圖7 葉柵通道內(nèi)損失劃分Fig.7 Partition of loss in cascade

    (2)

    (3)

    (4)

    (5)

    ?wake=?t-?s-?p-?shock

    (6)

    由表4可知,在設(shè)計進氣條件下,吸力面損失大于壓力面和尾跡損失之和;尾跡損失略大于壓力面附面層損失;當載荷很靠前時(0.07-1.0葉型),尾跡損失較大,由圖3(a),該葉型安裝角較小、尾部彎度較大,尾跡范圍較大[圖5(a)];當峰值馬赫數(shù)過小(0.1-0.95葉型),吸力面附面層損失較大,由圖2(b),該葉型安裝角較大且最大撓度位置靠近尾部,尾部載荷較大[圖4(b)]。尾部載荷大,導致吸力面或壓力面附面層快速增長;同時會使吸力面和壓力面附面層在尾緣交匯時方向差別增大,增加尾跡區(qū)摻混損失[圖6(a)];最終造成尾跡損失增大。載荷沿弦向逐漸減小,有利于減小流動損失。

    表4 設(shè)計進氣條件下不同葉型損失構(gòu)成Table 4 Loss component for different airfoils under design inlet condition

    2.2 全工況性能分析

    圖8~圖10為在進氣馬赫數(shù)分別為0.6、0.7(設(shè)計進氣馬赫數(shù))、0.8下的葉柵特性圖。由圖8~圖10可知,除去峰值馬赫數(shù)最大(0.1-1.14)葉型,不同進氣馬赫數(shù)下最低損失氣流角都是接近設(shè)計進氣角(42°),因此葉型表面等熵馬赫數(shù)符合可控擴散規(guī)律,則設(shè)計進氣角對應(yīng)于最低損失進氣角。由圖9(a)可以看出,在設(shè)計進氣馬赫數(shù)下,峰值馬赫數(shù)相同時,峰值位置前移,設(shè)計進氣角損失下降、低損失攻角范圍增大;但峰值位置過于靠前(位于0.07相對弦長處),最低損失對應(yīng)的進氣角減小、損失曲線整體左移。進氣馬赫數(shù)為0.6時,損失隨攻角變化規(guī)律與設(shè)計進氣馬赫數(shù)相同[圖8(a)]。進氣馬赫數(shù)為0.8時,相對位置在0.1~0.3,損失與攻角關(guān)系近于不隨相對位置變化;當相對位置為0.07時,損失曲線整體左移、且損失明顯增大。由圖8(b)、圖9(b)、圖10(b)可知,峰值馬赫數(shù)位置不變數(shù)值增加,最低損失對應(yīng)的進氣角增大;并且峰值馬赫數(shù)對性能影響較大,峰值馬赫數(shù)為1.0性能較好。綜合以上結(jié)果,吸力面峰值馬赫數(shù)及其位置直接影響設(shè)計工況性能;進而影響非設(shè)計工況性能。設(shè)計工況損失較低、對應(yīng)的低損失攻角范圍通常也相應(yīng)較大。

    圖8 0.6Ma下特性對比Fig.8 Comparison of performance at 0.6Ma

    圖9 0.7Ma下特性對比Fig.9 Comparison of performance at 0.7Ma

    圖10 0.8Ma下特性對比Fig.10 Comparison of performance at 0.8Ma

    圖11表明,進口馬赫數(shù)降至0.6時,在設(shè)計進氣角下,仍可保持與設(shè)計馬赫數(shù)近于相同的葉型表面等熵馬赫數(shù)分布規(guī)律;圖12表明,當進氣馬赫數(shù)增至0.8時,等熵馬赫數(shù)雖然仍保持與進氣馬赫數(shù)0.7近于相同的規(guī)律,但隨峰值點位置前移,峰值馬赫數(shù)下降[圖12(a)];隨峰值馬赫數(shù)增加,峰值點位置前移[圖12 (b)]。

    圖11 0.6Ma,設(shè)計進氣角下等熵馬赫數(shù)分布Fig.11 Isentropic Mach number distribution under design inlet flow angle at 0.6Ma

    圖12 0.8Ma,設(shè)計進氣角下等熵馬赫數(shù)分布Fig.12 Isentropic Mach number distribution under design inlet flow angle at 0.8Ma

    損失由吸力面、壓力面附面層損失、激波損失及尾跡損失構(gòu)成。由于壓力面馬赫數(shù)近于不變,損失也近于不變、并且占比較??;而以上設(shè)計例大部分葉型吸力面峰值馬赫數(shù)不超音或低超音,激波損失也較小。因此吸力面附面層損失和尾跡損失占比較大且變化較大,對總損失影響較大。在氣流轉(zhuǎn)角、進氣馬赫數(shù)一定情況下,葉型總載荷不變。改變吸力面峰值馬赫數(shù)和位置,則改變了載荷沿弦向分布。吸力面峰值點馬赫數(shù)位置前移、數(shù)值減小,吸力面由峰值減速流程較長、峰值后逆壓力梯度下降,有益于降低吸力面附面層發(fā)展、減小損失。但峰值馬赫數(shù)位置越靠前,則吸力面靠前緣凸起程度越大,有可能造成低損失攻角范圍減小、正攻角損失增加(如0.07-1.0葉型)。峰值馬赫數(shù)大于1.0,隨著峰值馬赫數(shù)增加,局部激波增強,不僅激波自身產(chǎn)生損失、其與吸力面附面層干擾造成附面層增厚會使附面層損失增加。因此,對于給定工況,吸力面馬赫數(shù)峰值及其位置存在最佳值。

    3 載荷分布規(guī)律對可控擴散轉(zhuǎn)子葉型性能影響

    選取進氣馬赫數(shù)為0.7、進氣角為56°、氣流轉(zhuǎn)角36°的轉(zhuǎn)子葉柵進行載荷分布規(guī)律研究。以吸力面峰值馬赫數(shù)0.85、峰值馬赫數(shù)位置為0.1、0.15、0.2倍相對軸向弦長;峰值位置為0.1倍相對軸向弦長、峰值馬赫數(shù)為0.85、0.9、0.95為目標,進行可控擴散轉(zhuǎn)子葉型優(yōu)化設(shè)計。

    圖13為優(yōu)化所得轉(zhuǎn)子葉型,圖13(a)中0.1-0.85、0.15-0.85、0.2-0.85表示峰值位置為0.1、0.15、0.2倍軸向相對弦長處、峰值馬赫數(shù)為0.85的葉型;圖13(b)中0.1-0.85、0.1-0.9、0.1-0.95表示峰值位置為0.1倍軸向弦長處、峰值馬赫數(shù)為0.85、0.9、0.95的葉型。由圖13、表5可知,峰值馬赫數(shù)不變時,隨著載荷后移,葉型厚度增加(趨勢與靜子一致);峰值位置不變時,隨著峰值馬赫數(shù)增大,葉型前部凸起程度增加(趨勢與靜子一致)。各葉型最大相對厚度為0.6~0.7(明顯高于靜子)、最大厚度相對位置為0.3~0.45(與靜子相近);最大相對撓度在0.9~0.10(與靜子相近)、最大厚度相對位置在0.32~0.36(與靜子相近)。

    圖13 不同載荷分布轉(zhuǎn)子葉型對比Fig.13 Comparison of rotor airfoils with different load distribution

    表5 可控擴散轉(zhuǎn)子葉型幾何參數(shù)對比Table 5 Comparison of geometry parameters of controlled diffusion airfoils

    3.1 設(shè)計工況流動分析

    由表6可以看出,采用與靜子相同優(yōu)化方法所得轉(zhuǎn)子葉型氣流轉(zhuǎn)角及進口馬赫數(shù)也近于不變。由圖14~圖16,設(shè)計進氣條件下,當峰值馬赫數(shù)不變時,隨著載荷前移,吸力面減速區(qū)變長、擴壓梯度減小,損失減?。环逯雕R赫數(shù)位于0.1倍軸向弦長處,隨著峰值馬赫數(shù)增大,峰值后逆壓梯度增大,損失增大。

    表 6 設(shè)計進氣條件下氣動參數(shù)對比Table 6 Comparison of aerodynamic parameters under design inlet condition

    圖14 設(shè)計進氣條件下等熵馬赫數(shù)分布對比Fig.14 Comparison of isentropic Mach number distribution under design inlet condition

    圖15 設(shè)計進氣條件下不同峰值馬赫數(shù)位置葉型馬赫數(shù)云圖對比Fig.15 Comparison of Mach number contours of different suction side peak Mach number location airfoils

    圖16 設(shè)計進氣條件下不同峰值馬赫數(shù)葉型馬赫數(shù)云圖對比Fig.16 Comparison of Mach number contours of different suction side peak Mach number airfoils

    圖14為設(shè)計進氣條件下不同載荷分布轉(zhuǎn)子葉型表面等熵馬赫數(shù)分布。由圖14(a)可知,保持吸力面峰值馬赫數(shù)不變,峰值位置后移,壓力面前緣區(qū)域馬赫數(shù)增大,葉型前部載荷減??;由圖14(b)可知,保持吸力面峰值馬赫數(shù)位置不變、峰值馬赫數(shù)增大,葉型前部載荷增加、后部載荷減小。由圖14可以看出,由于進氣條件不變、氣流轉(zhuǎn)角不變,各葉型總載荷近于不變。

    3.2 全工況流動分析

    圖17~圖19為進氣馬赫數(shù)為0.6、0.7、0.8下的葉柵特性圖。由圖17(a)、圖18(a)和圖19(a)可以看出,峰值馬赫數(shù)相同時,隨著峰值位置前移,設(shè)計進氣條件損失增大、且最低損失對應(yīng)的進氣角增加。由圖18(b)可知,相同峰值位置,隨峰值馬赫數(shù)增加,大攻角損失近于不變、小攻角損失增加明顯。由圖17(b)可知,進氣馬赫數(shù)為0.6時,也近似有以上規(guī)律;圖19(b)表明,進氣馬赫數(shù)為0.8時,峰值馬赫數(shù)為0.85和0.90損失曲線近于相同,峰值馬赫數(shù)為0.95時,損失大幅度增加。

    圖17 0.6Ma下特性對比Fig.17 Comparison of performance at 0.6Ma

    圖18 0.7Ma下特性對比Fig.18 Comparison of performance at 0.7Ma

    圖19 0.8Ma下特性對比Fig.19 Comparison of performance at 0.8Ma

    取出各部分流動損失,進一步對不同載荷分布轉(zhuǎn)子葉型流動損失構(gòu)成進行分析。由表7可以看出,在設(shè)計進氣條件下,吸力面附面層損失大于壓力面附面層損失與尾跡損失之和;壓力面附面層損失略大于尾跡損失;載荷越靠后,尾跡損失增長越明顯,以上規(guī)律與靜子相同。由圖15、表7可知,隨著載荷后移,轉(zhuǎn)子葉型厚度增大,使得附面層損失增加、尾跡摻混劇烈;由圖16、表7可知,由于設(shè)計進氣條件下峰值馬赫數(shù)越大的葉型前緣角越大,0.1-0.95葉型壓力面流動不能較好地附著在葉型表面、壓力面附面層損失較大。

    表7 設(shè)計進氣條件下不同葉型損失組成Table 7 Loss component for different airfoils under design inlet condition

    從圖20可以看出,當進氣馬赫數(shù)降至0.6時,設(shè)計進氣角下葉型表面等熵馬赫數(shù)分布與設(shè)計進氣馬赫數(shù)下規(guī)律性一致;由圖21可以看出,當進氣馬赫數(shù)增至0.8時,葉型表面等熵馬赫數(shù)分布與設(shè)計進氣馬赫數(shù)規(guī)律近于一致,但各葉型峰值馬赫數(shù)隨峰值位置前移而增大[圖21(a)]。

    圖20 0.6Ma,設(shè)計進氣角下等熵馬赫數(shù)分布Fig.20 Isentropic Mach number distribution under design inlet flow angle,at 0.6Ma

    圖21 0.8Ma,設(shè)計進氣角下等熵馬赫數(shù)分布Fig.21 Isentropic Mach number distribution under design inlet flow angle,at 0.8Ma

    4 結(jié)論

    建立了可控擴散葉型自動優(yōu)化設(shè)計方法,并對葉柵進口馬赫數(shù)為0.7的靜子和轉(zhuǎn)子葉型進行多個載荷分布設(shè)計,分析載荷分布規(guī)律對葉柵性能的影響。得出以下主要結(jié)論。

    (1)給定葉柵進出口馬赫數(shù)、吸力面峰值點等熵馬赫數(shù)及其位置,構(gòu)造直線和二次曲線組合的吸力面目標等熵馬赫數(shù)分布;以壓力面在0.1、0.4、0.8倍相對軸向弦長位置等熵馬赫數(shù)差值極小為目標,以實現(xiàn)壓力面等熵馬赫數(shù)近于不變。采用自動優(yōu)化方法可實現(xiàn)可控擴散葉型設(shè)計、以及吸力面峰值點馬赫數(shù)位置和數(shù)值的控制。

    (2)對于靜子和轉(zhuǎn)子葉型,葉型表面等熵馬赫數(shù)符合控制擴散規(guī)律,則設(shè)計進氣角對應(yīng)于最低損失進氣角。

    (3)對于可控擴散靜子葉型,吸力面峰值馬赫數(shù)位于0.1~0.3相對軸向弦長范圍內(nèi),峰值馬赫數(shù)相同峰值位置前移,則峰值后擴壓段變長、吸力面附面層逆壓力梯度減小,總損失減小且低損失攻角范圍增大;峰值馬赫數(shù)位置不變數(shù)值增加,吸力面靠前緣凸起程度越大,小攻角損失增加、最低損失對應(yīng)的進氣角增大。

    (4)對于可控擴散轉(zhuǎn)子葉型,在設(shè)計進氣馬赫數(shù)時,設(shè)計進氣角接近低損失左邊界,因此在轉(zhuǎn)子葉型設(shè)計時可考慮在大于設(shè)計進氣角下進行可控擴散葉型設(shè)計。峰值馬赫數(shù)相同時,隨著峰值位置后移,設(shè)計進氣角損失增大、且最低損失對應(yīng)的進氣角增加;相同峰值位置,隨峰值馬赫數(shù)增加,大攻角損失近于不變、小攻角損失增加明顯。

    猜你喜歡
    附面層葉型馬赫數(shù)
    高馬赫數(shù)激波作用下單模界面的Richtmyer-Meshkov不穩(wěn)定性數(shù)值模擬
    爆炸與沖擊(2024年7期)2024-11-01 00:00:00
    葉片前緣對吸力面邊界層3維流動影響分析
    先進動葉平面葉柵試驗和數(shù)值分析
    東方汽輪機(2022年1期)2022-04-13 02:25:00
    楔形葉片旋轉(zhuǎn)空化器葉型改進數(shù)值模擬研究
    基于網(wǎng)格框架的非結(jié)構(gòu)附面層網(wǎng)格生成技術(shù)
    一維非等熵可壓縮微極流體的低馬赫數(shù)極限
    基于數(shù)值模擬的流場附面層邊緣識別方法
    超聲壓氣機葉柵流場的數(shù)值模擬與試驗驗證
    汽輪機葉型幾何特性及機械特性研究
    熱力透平(2014年1期)2014-12-03 05:36:30
    NF-6連續(xù)式跨聲速風洞馬赫數(shù)控制方式比較與研究
    亚洲欧美精品综合久久99| 欧美日韩一级在线毛片| 国产真实乱freesex| aaaaa片日本免费| 色精品久久人妻99蜜桃| 国产精品 国内视频| 黄色片一级片一级黄色片| 免费在线观看亚洲国产| 99热这里只有精品一区 | 亚洲国产看品久久| 91成年电影在线观看| 欧美日韩亚洲综合一区二区三区_| 十八禁人妻一区二区| 桃色一区二区三区在线观看| 禁无遮挡网站| 国产成人啪精品午夜网站| 精品国内亚洲2022精品成人| 中文亚洲av片在线观看爽| 2021天堂中文幕一二区在线观 | 美女扒开内裤让男人捅视频| 国产精品电影一区二区三区| 伊人久久大香线蕉亚洲五| 国产精品久久久人人做人人爽| 国产精品精品国产色婷婷| 99riav亚洲国产免费| 国内精品久久久久久久电影| 日本三级黄在线观看| 精品卡一卡二卡四卡免费| 在线观看免费视频日本深夜| 久久国产精品男人的天堂亚洲| 国产又爽黄色视频| 香蕉丝袜av| 制服诱惑二区| 久久伊人香网站| 人人澡人人妻人| 啦啦啦韩国在线观看视频| 国产激情久久老熟女| www.精华液| 久久久久久亚洲精品国产蜜桃av| 日韩一卡2卡3卡4卡2021年| 国产片内射在线| 国产精品亚洲美女久久久| 老汉色∧v一级毛片| 午夜a级毛片| 淫秽高清视频在线观看| 一进一出抽搐动态| xxx96com| 国产精品亚洲一级av第二区| 国产高清激情床上av| 午夜福利免费观看在线| 精品熟女少妇八av免费久了| 91大片在线观看| 91在线观看av| 又紧又爽又黄一区二区| 久久精品aⅴ一区二区三区四区| 国产精品精品国产色婷婷| 后天国语完整版免费观看| 久久久久国产一级毛片高清牌| 99国产精品99久久久久| 每晚都被弄得嗷嗷叫到高潮| 一二三四社区在线视频社区8| 人成视频在线观看免费观看| 精品久久蜜臀av无| 99国产精品一区二区三区| 国产成人欧美| 99热只有精品国产| 亚洲五月婷婷丁香| 人人澡人人妻人| 白带黄色成豆腐渣| 18禁黄网站禁片免费观看直播| 亚洲精品美女久久久久99蜜臀| 免费电影在线观看免费观看| 成人亚洲精品一区在线观看| 夜夜夜夜夜久久久久| 怎么达到女性高潮| 精品国产亚洲在线| 岛国视频午夜一区免费看| 欧美zozozo另类| 男人操女人黄网站| 十分钟在线观看高清视频www| 日本在线视频免费播放| 久热爱精品视频在线9| 国产精品久久久av美女十八| 狠狠狠狠99中文字幕| 老司机午夜十八禁免费视频| 色播在线永久视频| av中文乱码字幕在线| 久久久水蜜桃国产精品网| 1024视频免费在线观看| 午夜成年电影在线免费观看| 亚洲av五月六月丁香网| 久久精品国产亚洲av香蕉五月| 黑丝袜美女国产一区| 日本一区二区免费在线视频| 99久久99久久久精品蜜桃| 久久久久久免费高清国产稀缺| 久久狼人影院| 亚洲aⅴ乱码一区二区在线播放 | 一边摸一边抽搐一进一小说| av片东京热男人的天堂| 欧美日本视频| 欧美av亚洲av综合av国产av| 黄色视频,在线免费观看| 一级作爱视频免费观看| 成人亚洲精品av一区二区| 久久精品亚洲精品国产色婷小说| 亚洲国产精品999在线| 国产精品免费一区二区三区在线| 日韩高清综合在线| 最近最新中文字幕大全电影3 | 亚洲精品久久成人aⅴ小说| 天天添夜夜摸| 国产av又大| 久久久久精品国产欧美久久久| 亚洲国产精品合色在线| 国产主播在线观看一区二区| 色精品久久人妻99蜜桃| 亚洲第一电影网av| 国产亚洲av高清不卡| 亚洲国产精品sss在线观看| 一进一出好大好爽视频| 18禁裸乳无遮挡免费网站照片 | 久久久久久人人人人人| 大型黄色视频在线免费观看| 丁香欧美五月| 一a级毛片在线观看| 老汉色av国产亚洲站长工具| 国产伦人伦偷精品视频| 最近最新免费中文字幕在线| 99久久无色码亚洲精品果冻| 日韩大码丰满熟妇| 国产成人影院久久av| 欧美人与性动交α欧美精品济南到| 好男人在线观看高清免费视频 | 欧美绝顶高潮抽搐喷水| 欧美 亚洲 国产 日韩一| 侵犯人妻中文字幕一二三四区| 亚洲一区二区三区色噜噜| 免费看a级黄色片| 婷婷六月久久综合丁香| 哪里可以看免费的av片| 一区二区三区精品91| 亚洲欧洲精品一区二区精品久久久| 婷婷精品国产亚洲av| 草草在线视频免费看| cao死你这个sao货| 久久人妻av系列| 亚洲五月色婷婷综合| 久久精品国产清高在天天线| 午夜精品在线福利| 成人18禁高潮啪啪吃奶动态图| 精品少妇一区二区三区视频日本电影| 亚洲熟妇中文字幕五十中出| 久久性视频一级片| 美女扒开内裤让男人捅视频| 熟女电影av网| 亚洲国产欧洲综合997久久, | 在线观看66精品国产| 国产一区在线观看成人免费| 精品久久久久久久末码| 久久午夜综合久久蜜桃| 亚洲全国av大片| 18禁国产床啪视频网站| 国产熟女午夜一区二区三区| 国产成人av激情在线播放| 久久精品91蜜桃| 欧美黄色片欧美黄色片| 操出白浆在线播放| 色综合婷婷激情| av中文乱码字幕在线| 国产亚洲精品久久久久5区| 成年女人毛片免费观看观看9| 熟妇人妻久久中文字幕3abv| 国产高清激情床上av| 一级片免费观看大全| 久久九九热精品免费| 一本精品99久久精品77| 12—13女人毛片做爰片一| 88av欧美| 亚洲欧美激情综合另类| 欧美乱码精品一区二区三区| 一区福利在线观看| 日本 欧美在线| 国产片内射在线| 国产午夜福利久久久久久| 色在线成人网| 久久久久久久精品吃奶| 亚洲国产精品久久男人天堂| 不卡一级毛片| avwww免费| 一级片免费观看大全| 国产aⅴ精品一区二区三区波| 精品不卡国产一区二区三区| 老汉色av国产亚洲站长工具| 人人妻人人澡人人看| 日本一本二区三区精品| 久久久久久人人人人人| 91av网站免费观看| 伦理电影免费视频| 色哟哟哟哟哟哟| 啦啦啦 在线观看视频| 热99re8久久精品国产| 欧美乱码精品一区二区三区| av在线播放免费不卡| 1024视频免费在线观看| 免费观看精品视频网站| 好看av亚洲va欧美ⅴa在| 在线观看66精品国产| 国产欧美日韩一区二区三| 大香蕉久久成人网| 国产成人精品久久二区二区免费| 18禁黄网站禁片免费观看直播| cao死你这个sao货| 国产成人精品久久二区二区免费| 亚洲精品色激情综合| 久久久久久亚洲精品国产蜜桃av| 老司机在亚洲福利影院| 亚洲国产看品久久| 后天国语完整版免费观看| 人人妻人人看人人澡| 波多野结衣av一区二区av| 日韩av在线大香蕉| 欧美精品亚洲一区二区| 视频区欧美日本亚洲| 校园春色视频在线观看| 午夜福利在线在线| 一本久久中文字幕| 亚洲av五月六月丁香网| 亚洲人成77777在线视频| 黄色丝袜av网址大全| 久久欧美精品欧美久久欧美| tocl精华| 黄色成人免费大全| videosex国产| 国产日本99.免费观看| 久久人妻av系列| 在线国产一区二区在线| 亚洲熟女毛片儿| 波多野结衣av一区二区av| 成人亚洲精品一区在线观看| 成人亚洲精品av一区二区| 一级黄色大片毛片| 国产一区二区三区在线臀色熟女| 精品午夜福利视频在线观看一区| 在线十欧美十亚洲十日本专区| 色综合亚洲欧美另类图片| 别揉我奶头~嗯~啊~动态视频| 在线观看舔阴道视频| av欧美777| 国产色视频综合| 国产亚洲精品一区二区www| 国产精品香港三级国产av潘金莲| 精品欧美国产一区二区三| 午夜日韩欧美国产| 久热这里只有精品99| 男人的好看免费观看在线视频 | 一本综合久久免费| 日韩三级视频一区二区三区| 久久午夜综合久久蜜桃| 久久亚洲精品不卡| 成年免费大片在线观看| 亚洲欧美精品综合一区二区三区| 18禁美女被吸乳视频| 99精品在免费线老司机午夜| 男人舔女人下体高潮全视频| 国产成人一区二区三区免费视频网站| 亚洲色图av天堂| 久久久水蜜桃国产精品网| 亚洲七黄色美女视频| 中出人妻视频一区二区| 亚洲成人免费电影在线观看| 草草在线视频免费看| 日本a在线网址| 黄片播放在线免费| 波多野结衣av一区二区av| 两性夫妻黄色片| 91在线观看av| 天堂动漫精品| 校园春色视频在线观看| 自线自在国产av| 欧美成人性av电影在线观看| 亚洲精品av麻豆狂野| 日韩 欧美 亚洲 中文字幕| 欧美色视频一区免费| 欧美乱色亚洲激情| 天堂√8在线中文| 久久草成人影院| 亚洲电影在线观看av| 国产一区二区在线av高清观看| 亚洲第一电影网av| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲一码二码三码区别大吗| 中文字幕久久专区| 免费av毛片视频| 免费在线观看成人毛片| 国产精品亚洲一级av第二区| 黄色视频,在线免费观看| 一本精品99久久精品77| 91成人精品电影| 青草久久国产| 无限看片的www在线观看| 久久久久久久久久黄片| 免费看日本二区| 久久久久久免费高清国产稀缺| 制服丝袜大香蕉在线| 91大片在线观看| 亚洲成人精品中文字幕电影| 亚洲一卡2卡3卡4卡5卡精品中文| 国产黄片美女视频| 免费女性裸体啪啪无遮挡网站| 久久热在线av| 在线观看一区二区三区| 国产97色在线日韩免费| 国产高清视频在线播放一区| 欧美日韩乱码在线| 欧美一级a爱片免费观看看 | 韩国精品一区二区三区| 中亚洲国语对白在线视频| 热re99久久国产66热| 一级毛片女人18水好多| 亚洲av片天天在线观看| 日韩精品免费视频一区二区三区| 亚洲成av人片免费观看| 亚洲专区国产一区二区| 日本一本二区三区精品| 男人舔女人的私密视频| 国产精品 国内视频| 亚洲国产日韩欧美精品在线观看 | 亚洲av熟女| 老司机福利观看| 久久婷婷成人综合色麻豆| 丝袜美腿诱惑在线| 国产片内射在线| 亚洲精品中文字幕在线视频| 亚洲国产精品成人综合色| 日日爽夜夜爽网站| 久久精品国产清高在天天线| 国产视频内射| 身体一侧抽搐| 亚洲成人精品中文字幕电影| 亚洲av电影在线进入| 99久久精品国产亚洲精品| 久久久精品国产亚洲av高清涩受| 天天添夜夜摸| 非洲黑人性xxxx精品又粗又长| 中文字幕久久专区| 日韩中文字幕欧美一区二区| 日日夜夜操网爽| 12—13女人毛片做爰片一| 亚洲欧美日韩无卡精品| 日韩欧美一区视频在线观看| 18美女黄网站色大片免费观看| 成人18禁在线播放| 国产精品98久久久久久宅男小说| 日本 av在线| 黑丝袜美女国产一区| 无人区码免费观看不卡| 女人爽到高潮嗷嗷叫在线视频| 日本五十路高清| 欧美黑人欧美精品刺激| 午夜激情福利司机影院| 国产激情偷乱视频一区二区| 日本免费一区二区三区高清不卡| 成人亚洲精品一区在线观看| 国产成人欧美在线观看| 婷婷亚洲欧美| 色综合亚洲欧美另类图片| bbb黄色大片| 窝窝影院91人妻| 啦啦啦免费观看视频1| 禁无遮挡网站| 级片在线观看| 悠悠久久av| 亚洲国产精品成人综合色| 国内精品久久久久精免费| 一级a爱片免费观看的视频| 久久精品91无色码中文字幕| 看黄色毛片网站| 国产主播在线观看一区二区| 狂野欧美激情性xxxx| 熟妇人妻久久中文字幕3abv| 亚洲男人天堂网一区| 操出白浆在线播放| 国产熟女xx| 色综合亚洲欧美另类图片| 久久精品91无色码中文字幕| 成年女人毛片免费观看观看9| 亚洲成人免费电影在线观看| 视频区欧美日本亚洲| 精品欧美一区二区三区在线| 国产乱人伦免费视频| 深夜精品福利| 可以免费在线观看a视频的电影网站| 美女国产高潮福利片在线看| 一二三四社区在线视频社区8| 老司机午夜福利在线观看视频| 国产午夜福利久久久久久| 亚洲男人天堂网一区| 欧美黄色淫秽网站| 精品国产超薄肉色丝袜足j| 最近最新免费中文字幕在线| 少妇熟女aⅴ在线视频| 99在线人妻在线中文字幕| 中出人妻视频一区二区| 日本成人三级电影网站| 亚洲av成人av| 两人在一起打扑克的视频| 黑人巨大精品欧美一区二区mp4| 久久久久国内视频| 亚洲精品国产区一区二| 午夜老司机福利片| 久久国产乱子伦精品免费另类| 十分钟在线观看高清视频www| 97人妻精品一区二区三区麻豆 | 在线国产一区二区在线| 搡老妇女老女人老熟妇| 老熟妇乱子伦视频在线观看| 桃红色精品国产亚洲av| 每晚都被弄得嗷嗷叫到高潮| 日本撒尿小便嘘嘘汇集6| 在线视频色国产色| 成人18禁在线播放| 久久精品影院6| 国产激情偷乱视频一区二区| 免费观看人在逋| 少妇熟女aⅴ在线视频| 亚洲国产欧美一区二区综合| 亚洲成人免费电影在线观看| 可以在线观看的亚洲视频| 亚洲精品在线观看二区| 欧美日韩亚洲国产一区二区在线观看| 国产精品久久久人人做人人爽| 男女那种视频在线观看| 视频区欧美日本亚洲| 国产精品亚洲美女久久久| 亚洲av熟女| 日本成人三级电影网站| 在线观看日韩欧美| 久久人人精品亚洲av| 欧美日韩亚洲国产一区二区在线观看| 欧美激情 高清一区二区三区| 午夜免费鲁丝| 欧美性猛交╳xxx乱大交人| 国产精品乱码一区二三区的特点| 亚洲五月天丁香| 精品国产亚洲在线| 很黄的视频免费| 黄片小视频在线播放| 久久午夜综合久久蜜桃| 哪里可以看免费的av片| 亚洲国产欧美网| 俄罗斯特黄特色一大片| 成人国产综合亚洲| 99riav亚洲国产免费| 少妇粗大呻吟视频| e午夜精品久久久久久久| 成人三级黄色视频| 日韩欧美在线二视频| 欧美黑人欧美精品刺激| 男男h啪啪无遮挡| 国产精品二区激情视频| 亚洲天堂国产精品一区在线| 精品国内亚洲2022精品成人| 亚洲全国av大片| 真人一进一出gif抽搐免费| av免费在线观看网站| 国产精品自产拍在线观看55亚洲| 夜夜躁狠狠躁天天躁| 黄片小视频在线播放| 麻豆一二三区av精品| 国内毛片毛片毛片毛片毛片| 国产黄片美女视频| 制服诱惑二区| 人人妻人人澡欧美一区二区| 国产成人啪精品午夜网站| 99精品欧美一区二区三区四区| 亚洲中文日韩欧美视频| 搡老妇女老女人老熟妇| 国产亚洲av高清不卡| 日本一区二区免费在线视频| 精品国产一区二区三区四区第35| 黄色a级毛片大全视频| 嫩草影院精品99| 婷婷亚洲欧美| 99在线视频只有这里精品首页| 国产三级在线视频| 午夜成年电影在线免费观看| 亚洲av成人av| 欧美 亚洲 国产 日韩一| 18禁裸乳无遮挡免费网站照片 | 成人免费观看视频高清| 老司机深夜福利视频在线观看| 丝袜人妻中文字幕| 午夜福利高清视频| 欧美一级毛片孕妇| 1024香蕉在线观看| 国产激情偷乱视频一区二区| 欧美黄色淫秽网站| 久久99热这里只有精品18| 我的亚洲天堂| 日本免费a在线| 成人国产综合亚洲| 亚洲精华国产精华精| 又黄又爽又免费观看的视频| 黄色视频,在线免费观看| 亚洲国产欧美网| 午夜福利免费观看在线| 麻豆一二三区av精品| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品美女久久久久99蜜臀| 香蕉久久夜色| 视频在线观看一区二区三区| 久久久久久人人人人人| 色综合欧美亚洲国产小说| 欧美午夜高清在线| 99国产精品99久久久久| 久久亚洲真实| 久久久久久人人人人人| 成人三级做爰电影| 久久精品国产清高在天天线| 久久精品国产综合久久久| 久久精品91无色码中文字幕| 伊人久久大香线蕉亚洲五| 国产精品久久电影中文字幕| 国产欧美日韩精品亚洲av| 女人爽到高潮嗷嗷叫在线视频| 久久久久国内视频| 人人妻,人人澡人人爽秒播| 亚洲国产精品久久男人天堂| 欧美性猛交黑人性爽| 国语自产精品视频在线第100页| 欧美日韩一级在线毛片| 在线观看一区二区三区| 成人欧美大片| 母亲3免费完整高清在线观看| 午夜视频精品福利| 狂野欧美激情性xxxx| 亚洲第一青青草原| 日日爽夜夜爽网站| 母亲3免费完整高清在线观看| 国产高清视频在线播放一区| 熟女少妇亚洲综合色aaa.| 天堂√8在线中文| 搞女人的毛片| 在线观看免费日韩欧美大片| 日韩成人在线观看一区二区三区| xxx96com| 久久精品国产亚洲av高清一级| 老汉色∧v一级毛片| 色综合欧美亚洲国产小说| 亚洲熟妇熟女久久| 老司机在亚洲福利影院| 三级毛片av免费| 黄色女人牲交| 91九色精品人成在线观看| 午夜免费激情av| 亚洲欧美日韩无卡精品| 怎么达到女性高潮| 免费人成视频x8x8入口观看| 嫩草影视91久久| 国产高清videossex| 亚洲狠狠婷婷综合久久图片| 国产精品亚洲一级av第二区| 欧美丝袜亚洲另类 | 欧美日本视频| 一边摸一边抽搐一进一小说| 最好的美女福利视频网| 亚洲天堂国产精品一区在线| 中文字幕另类日韩欧美亚洲嫩草| 国产精品免费视频内射| 亚洲成a人片在线一区二区| 1024香蕉在线观看| 久久精品亚洲精品国产色婷小说| 91国产中文字幕| 神马国产精品三级电影在线观看 | 免费看日本二区| 亚洲男人天堂网一区| 久久精品国产综合久久久| 国产成人精品久久二区二区91| 久久香蕉精品热| 最新在线观看一区二区三区| 长腿黑丝高跟| 老熟妇乱子伦视频在线观看| 很黄的视频免费| 亚洲男人天堂网一区| 又黄又爽又免费观看的视频| 国产成人一区二区三区免费视频网站| 香蕉av资源在线| 两性午夜刺激爽爽歪歪视频在线观看 | 国产成年人精品一区二区| 不卡av一区二区三区| 一本一本综合久久| avwww免费| 91在线观看av| 真人一进一出gif抽搐免费| 一区福利在线观看| 韩国av一区二区三区四区| 国产av又大| 久久精品aⅴ一区二区三区四区| 韩国av一区二区三区四区| 男人操女人黄网站| 久久精品国产亚洲av香蕉五月| 亚洲av成人一区二区三| 欧美成狂野欧美在线观看| 国内揄拍国产精品人妻在线 | 深夜精品福利| 亚洲av五月六月丁香网| 亚洲成a人片在线一区二区| 深夜精品福利| 国产成+人综合+亚洲专区| 亚洲精华国产精华精| 悠悠久久av| 免费在线观看完整版高清| 成熟少妇高潮喷水视频| 1024香蕉在线观看|