• <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ù)控制方式比較與研究
    18禁观看日本| 国产一区二区在线观看av| 亚洲av日韩精品久久久久久密 | 99香蕉大伊视频| 大香蕉久久网| 亚洲av成人不卡在线观看播放网 | 久久女婷五月综合色啪小说| 久久人人97超碰香蕉20202| 美女国产高潮福利片在线看| 一级毛片我不卡| 国产成人av激情在线播放| 一本色道久久久久久精品综合| 国产 一区精品| 高清黄色对白视频在线免费看| 亚洲第一青青草原| 热re99久久精品国产66热6| 亚洲国产最新在线播放| 2018国产大陆天天弄谢| 中文欧美无线码| 老司机影院毛片| 人人妻人人爽人人添夜夜欢视频| 国产亚洲av高清不卡| 9191精品国产免费久久| 最近的中文字幕免费完整| 可以免费在线观看a视频的电影网站 | a级片在线免费高清观看视频| 日韩精品免费视频一区二区三区| 赤兔流量卡办理| 久久狼人影院| 精品视频人人做人人爽| 狠狠精品人妻久久久久久综合| 超碰成人久久| 天堂8中文在线网| 国产极品粉嫩免费观看在线| 天天操日日干夜夜撸| av有码第一页| 精品亚洲成a人片在线观看| 亚洲国产欧美日韩在线播放| 国产精品久久久久成人av| 人人妻人人添人人爽欧美一区卜| 国产亚洲午夜精品一区二区久久| 在线观看免费日韩欧美大片| 国产福利在线免费观看视频| 国产成人精品在线电影| 国产精品久久久久成人av| 国产精品免费大片| 日韩制服骚丝袜av| 一区二区三区精品91| 亚洲一区二区三区欧美精品| 伊人亚洲综合成人网| 亚洲欧美精品自产自拍| 最近中文字幕2019免费版| 丰满迷人的少妇在线观看| 精品国产露脸久久av麻豆| 人成视频在线观看免费观看| 午夜av观看不卡| 成年美女黄网站色视频大全免费| 99国产精品免费福利视频| 亚洲av在线观看美女高潮| 国产精品一二三区在线看| 日本欧美国产在线视频| 午夜日韩欧美国产| 99久久综合免费| 极品人妻少妇av视频| 久久久久久久久久久久大奶| 国产一区二区激情短视频 | 少妇人妻久久综合中文| 国产麻豆69| 超色免费av| 这个男人来自地球电影免费观看 | 国精品久久久久久国模美| 久久午夜综合久久蜜桃| e午夜精品久久久久久久| 又大又黄又爽视频免费| 中文欧美无线码| 免费人妻精品一区二区三区视频| 美女脱内裤让男人舔精品视频| 老汉色av国产亚洲站长工具| 国产又色又爽无遮挡免| 欧美在线黄色| 欧美 日韩 精品 国产| 国产精品国产三级国产专区5o| 下体分泌物呈黄色| 国产精品久久久久成人av| 久久韩国三级中文字幕| 久久人人爽av亚洲精品天堂| 男女边吃奶边做爰视频| 日韩大片免费观看网站| 操出白浆在线播放| 国产精品国产三级专区第一集| 纵有疾风起免费观看全集完整版| 黄频高清免费视频| 久久天躁狠狠躁夜夜2o2o | 亚洲av日韩精品久久久久久密 | 欧美 日韩 精品 国产| 婷婷色av中文字幕| 免费黄色在线免费观看| 亚洲三区欧美一区| 啦啦啦啦在线视频资源| 婷婷色av中文字幕| 亚洲国产成人一精品久久久| 侵犯人妻中文字幕一二三四区| 亚洲精品国产av成人精品| 五月天丁香电影| 亚洲精品日本国产第一区| 在线 av 中文字幕| 国产精品免费视频内射| 秋霞在线观看毛片| 999精品在线视频| 亚洲自偷自拍图片 自拍| 午夜免费鲁丝| 91老司机精品| 精品一区二区三区av网在线观看 | 日韩av不卡免费在线播放| 九草在线视频观看| 美女脱内裤让男人舔精品视频| 亚洲七黄色美女视频| 中文字幕另类日韩欧美亚洲嫩草| 久久99精品国语久久久| 亚洲精品久久久久久婷婷小说| 久久久精品国产亚洲av高清涩受| 少妇人妻 视频| 亚洲成国产人片在线观看| 欧美另类一区| 国产亚洲最大av| 欧美日韩视频高清一区二区三区二| 大片免费播放器 马上看| 一本大道久久a久久精品| 最近2019中文字幕mv第一页| 亚洲国产欧美日韩在线播放| 日韩 欧美 亚洲 中文字幕| 99国产精品免费福利视频| 国产成人精品久久久久久| 狠狠婷婷综合久久久久久88av| 永久免费av网站大全| 日本一区二区免费在线视频| 成人国语在线视频| 一级毛片电影观看| 国产精品一国产av| 欧美在线一区亚洲| 大香蕉久久网| 国产伦理片在线播放av一区| 欧美激情极品国产一区二区三区| 久久鲁丝午夜福利片| 国产精品.久久久| 欧美日韩成人在线一区二区| 亚洲国产av新网站| 亚洲av电影在线观看一区二区三区| 久久精品国产综合久久久| 黄色视频在线播放观看不卡| 男男h啪啪无遮挡| 一级黄片播放器| 亚洲精品国产区一区二| 伊人久久大香线蕉亚洲五| 国产色婷婷99| 亚洲精品一二三| 在线天堂最新版资源| 国产极品天堂在线| av片东京热男人的天堂| 青春草亚洲视频在线观看| 美国免费a级毛片| 国产乱来视频区| 99久久人妻综合| 精品少妇黑人巨大在线播放| 欧美精品人与动牲交sv欧美| 最新在线观看一区二区三区 | 在线精品无人区一区二区三| 亚洲中文av在线| 男女边摸边吃奶| 桃花免费在线播放| 久久热在线av| 午夜福利影视在线免费观看| 国产成人欧美| 久久久欧美国产精品| 国产男女超爽视频在线观看| 男女下面插进去视频免费观看| 街头女战士在线观看网站| 99热网站在线观看| 久久99热这里只频精品6学生| 五月天丁香电影| 精品亚洲乱码少妇综合久久| 国产成人午夜福利电影在线观看| 51午夜福利影视在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲综合色网址| 日本av手机在线免费观看| 女性被躁到高潮视频| 极品少妇高潮喷水抽搐| 亚洲视频免费观看视频| 99国产综合亚洲精品| 国产一区二区三区av在线| 欧美另类一区| 女人爽到高潮嗷嗷叫在线视频| 嫩草影视91久久| 美女高潮到喷水免费观看| 国产极品粉嫩免费观看在线| 国产成人欧美在线观看 | 久久女婷五月综合色啪小说| 午夜日本视频在线| 久久影院123| 国产精品久久久久久久久免| 国产成人精品福利久久| 日本av手机在线免费观看| 亚洲欧美一区二区三区久久| 婷婷成人精品国产| 成年动漫av网址| 另类亚洲欧美激情| 亚洲av日韩在线播放| 国产成人91sexporn| 青春草亚洲视频在线观看| 男男h啪啪无遮挡| 青春草国产在线视频| 亚洲精品久久成人aⅴ小说| 亚洲成人国产一区在线观看 | 十八禁网站网址无遮挡| 欧美精品人与动牲交sv欧美| bbb黄色大片| 夫妻午夜视频| 不卡av一区二区三区| 999久久久国产精品视频| 亚洲精华国产精华液的使用体验| 精品免费久久久久久久清纯 | 国产在线一区二区三区精| 啦啦啦 在线观看视频| 男女边摸边吃奶| 国产极品粉嫩免费观看在线| 久久久久精品国产欧美久久久 | 午夜福利乱码中文字幕| 老司机影院毛片| 看非洲黑人一级黄片| 日日爽夜夜爽网站| 亚洲精品在线美女| 蜜桃在线观看..| 老鸭窝网址在线观看| 精品国产一区二区久久| 久久97久久精品| 丝袜在线中文字幕| 亚洲,欧美精品.| 老汉色av国产亚洲站长工具| 纵有疾风起免费观看全集完整版| 亚洲欧美成人综合另类久久久| 校园人妻丝袜中文字幕| 久热爱精品视频在线9| 国产午夜精品一二区理论片| 日韩一本色道免费dvd| 在线观看三级黄色| 免费日韩欧美在线观看| 国产一区有黄有色的免费视频| 最近中文字幕2019免费版| 99精品久久久久人妻精品| 伦理电影免费视频| 免费日韩欧美在线观看| 丁香六月天网| 久久99一区二区三区| 七月丁香在线播放| 日本爱情动作片www.在线观看| 久久久精品区二区三区| 亚洲人成网站在线观看播放| 日韩大码丰满熟妇| av又黄又爽大尺度在线免费看| 最近中文字幕高清免费大全6| 国产亚洲av高清不卡| 考比视频在线观看| 久久久亚洲精品成人影院| 亚洲综合色网址| 在现免费观看毛片| 久久国产精品男人的天堂亚洲| 久久人妻熟女aⅴ| 99热网站在线观看| 两个人看的免费小视频| 下体分泌物呈黄色| 国产不卡av网站在线观看| 色精品久久人妻99蜜桃| 在线观看三级黄色| 无限看片的www在线观看| 国产乱人偷精品视频| 欧美日本中文国产一区发布| 美女主播在线视频| 亚洲,一卡二卡三卡| 999久久久国产精品视频| 在线精品无人区一区二区三| 18在线观看网站| av天堂久久9| 性色av一级| 精品少妇黑人巨大在线播放| 一区二区三区精品91| 久久99热这里只频精品6学生| 久久久久久久国产电影| 亚洲人成77777在线视频| 免费人妻精品一区二区三区视频| 搡老乐熟女国产| 日本欧美国产在线视频| 亚洲图色成人| 亚洲国产日韩一区二区| 亚洲国产av新网站| 曰老女人黄片| avwww免费| 欧美激情高清一区二区三区 | 久久精品国产综合久久久| 晚上一个人看的免费电影| 欧美激情 高清一区二区三区| 人成视频在线观看免费观看| 最新在线观看一区二区三区 | 波多野结衣一区麻豆| 久久97久久精品| 国产乱来视频区| 人妻 亚洲 视频| 国产成人精品久久久久久| 极品人妻少妇av视频| 日日爽夜夜爽网站| 亚洲精品国产区一区二| 免费看不卡的av| 男的添女的下面高潮视频| 一级毛片 在线播放| 亚洲精品av麻豆狂野| 精品国产超薄肉色丝袜足j| 热99久久久久精品小说推荐| 国产女主播在线喷水免费视频网站| 少妇人妻久久综合中文| 嫩草影院入口| 国产精品久久久久久精品古装| 一区二区日韩欧美中文字幕| 肉色欧美久久久久久久蜜桃| 亚洲av成人精品一二三区| 久久人人爽av亚洲精品天堂| 少妇人妻精品综合一区二区| 两个人免费观看高清视频| 天堂中文最新版在线下载| 夫妻午夜视频| 亚洲精品在线美女| 精品酒店卫生间| 中国三级夫妇交换| 亚洲精品日本国产第一区| 国产精品香港三级国产av潘金莲 | 91精品三级在线观看| 日本一区二区免费在线视频| 国产成人免费无遮挡视频| 久久精品国产a三级三级三级| 色网站视频免费| 日韩中文字幕视频在线看片| a级毛片黄视频| 国产av一区二区精品久久| 又大又爽又粗| 人人妻人人添人人爽欧美一区卜| 一二三四在线观看免费中文在| 一区福利在线观看| 久久天躁狠狠躁夜夜2o2o | 精品第一国产精品| 亚洲欧美激情在线| 两个人免费观看高清视频| 亚洲av国产av综合av卡| 搡老乐熟女国产| 在线观看www视频免费| 精品少妇黑人巨大在线播放| 日本色播在线视频| 黄色毛片三级朝国网站| 波多野结衣av一区二区av| 黄频高清免费视频| 中国三级夫妇交换| 亚洲国产最新在线播放| 精品亚洲成国产av| 久久久久久久久久久久大奶| 成人漫画全彩无遮挡| 国产一级毛片在线| 黄色一级大片看看| 久久这里只有精品19| 曰老女人黄片| 捣出白浆h1v1| 国产 一区精品| 啦啦啦中文免费视频观看日本| 亚洲国产精品国产精品| 香蕉丝袜av| 在线观看三级黄色| 亚洲欧美中文字幕日韩二区| 成人午夜精彩视频在线观看| 18在线观看网站| 精品少妇久久久久久888优播| 老司机靠b影院| 国产高清不卡午夜福利| 男女国产视频网站| 丁香六月天网| 国产一区二区三区av在线| 久久热在线av| 国产精品久久久久久久久免| 亚洲中文av在线| 国产成人91sexporn| 免费在线观看视频国产中文字幕亚洲 | 午夜影院在线不卡| 亚洲一码二码三码区别大吗| 2018国产大陆天天弄谢| 国产精品一区二区精品视频观看| 国产精品久久久人人做人人爽| 国产熟女午夜一区二区三区| 热99国产精品久久久久久7| 国产午夜精品一二区理论片| 国产国语露脸激情在线看| 色婷婷久久久亚洲欧美| 亚洲精品美女久久av网站| 美女扒开内裤让男人捅视频| 亚洲精品中文字幕在线视频| av.在线天堂| 久久久久久久久免费视频了| 新久久久久国产一级毛片| 制服丝袜香蕉在线| 国产激情久久老熟女| 久久久久久免费高清国产稀缺| 国产在线免费精品| 亚洲一级一片aⅴ在线观看| 久久精品亚洲熟妇少妇任你| 成人午夜精彩视频在线观看| 欧美亚洲 丝袜 人妻 在线| 精品一区二区三区av网在线观看 | 在线亚洲精品国产二区图片欧美| 十八禁人妻一区二区| 一区二区三区乱码不卡18| kizo精华| 亚洲精品乱久久久久久| 丰满饥渴人妻一区二区三| 午夜福利在线免费观看网站| 久久ye,这里只有精品| 无遮挡黄片免费观看| 亚洲精品,欧美精品| 国产精品99久久99久久久不卡 | 一二三四中文在线观看免费高清| 青春草视频在线免费观看| 看非洲黑人一级黄片| 欧美日韩av久久| 性高湖久久久久久久久免费观看| 男女高潮啪啪啪动态图| 乱人伦中国视频| 成人手机av| 色婷婷久久久亚洲欧美| 九草在线视频观看| 久久久久精品国产欧美久久久 | 国产免费现黄频在线看| 一边亲一边摸免费视频| 亚洲精品第二区| 人妻一区二区av| 国产精品免费大片| 亚洲欧美精品综合一区二区三区| 秋霞伦理黄片| 国产精品一国产av| 亚洲色图综合在线观看| 国产色婷婷99| 老司机影院毛片| e午夜精品久久久久久久| 国产精品久久久久久精品古装| 亚洲在久久综合| 亚洲欧洲国产日韩| 午夜福利影视在线免费观看| 在线观看一区二区三区激情| 亚洲精品日韩在线中文字幕| 成人18禁高潮啪啪吃奶动态图| 亚洲成人免费av在线播放| 国产福利在线免费观看视频| svipshipincom国产片| 日韩制服丝袜自拍偷拍| 国产精品偷伦视频观看了| 国产精品免费视频内射| 亚洲欧美精品自产自拍| 热re99久久国产66热| 亚洲,欧美精品.| 2018国产大陆天天弄谢| 亚洲av成人不卡在线观看播放网 | 国产无遮挡羞羞视频在线观看| 在线看a的网站| 99久久99久久久精品蜜桃| 午夜福利,免费看| 成人三级做爰电影| 亚洲精品一区蜜桃| 日本午夜av视频| 韩国av在线不卡| 亚洲av日韩在线播放| 亚洲婷婷狠狠爱综合网| 国产毛片在线视频| 中文精品一卡2卡3卡4更新| 老熟女久久久| 观看av在线不卡| 老汉色∧v一级毛片| 国产精品成人在线| 欧美成人午夜精品| 久久精品国产亚洲av高清一级| 一级毛片我不卡| 中文字幕另类日韩欧美亚洲嫩草| 日韩免费高清中文字幕av| 五月开心婷婷网| 亚洲情色 制服丝袜| 久久亚洲国产成人精品v| 青春草亚洲视频在线观看| 欧美黄色片欧美黄色片| 天堂8中文在线网| 人妻一区二区av| 人人妻人人添人人爽欧美一区卜| 91精品三级在线观看| 久久青草综合色| 久久鲁丝午夜福利片| 日日爽夜夜爽网站| 黄色一级大片看看| 中文字幕av电影在线播放| 看免费成人av毛片| 国产99久久九九免费精品| 在线精品无人区一区二区三| 欧美日韩视频高清一区二区三区二| 国产毛片在线视频| 我的亚洲天堂| 亚洲精品美女久久av网站| 日日爽夜夜爽网站| 女人久久www免费人成看片| 久久韩国三级中文字幕| 色94色欧美一区二区| 亚洲,一卡二卡三卡| 啦啦啦 在线观看视频| 亚洲情色 制服丝袜| 久久av网站| 精品一区二区三区av网在线观看 | 99久久人妻综合| 午夜av观看不卡| 少妇 在线观看| 男女免费视频国产| 午夜91福利影院| 黄色视频不卡| 9色porny在线观看| 欧美黑人精品巨大| 亚洲精品乱久久久久久| 熟妇人妻不卡中文字幕| 黄片无遮挡物在线观看| 国产精品一区二区精品视频观看| 叶爱在线成人免费视频播放| 亚洲,一卡二卡三卡| 18禁国产床啪视频网站| 久久午夜综合久久蜜桃| 日韩一区二区视频免费看| 中文字幕色久视频| 伊人久久国产一区二区| bbb黄色大片| 欧美成人精品欧美一级黄| 这个男人来自地球电影免费观看 | 天天添夜夜摸| 中文字幕高清在线视频| 看免费成人av毛片| 久久青草综合色| 高清视频免费观看一区二区| 亚洲av电影在线观看一区二区三区| 日韩精品免费视频一区二区三区| 亚洲精品乱久久久久久| 日日啪夜夜爽| 精品一区二区免费观看| 亚洲少妇的诱惑av| 亚洲精品乱久久久久久| videos熟女内射| 一区二区av电影网| 亚洲精品国产av蜜桃| 69精品国产乱码久久久| 老汉色∧v一级毛片| 久久精品亚洲av国产电影网| 狂野欧美激情性bbbbbb| 交换朋友夫妻互换小说| 五月天丁香电影| 最近最新中文字幕大全免费视频 | 天天躁日日躁夜夜躁夜夜| 国产在线视频一区二区| 天美传媒精品一区二区| 午夜福利,免费看| 欧美久久黑人一区二区| 男女床上黄色一级片免费看| 精品少妇久久久久久888优播| 国产成人一区二区在线| 欧美久久黑人一区二区| 欧美老熟妇乱子伦牲交| av天堂久久9| 国产一卡二卡三卡精品 | av女优亚洲男人天堂| 天天操日日干夜夜撸| 国产片特级美女逼逼视频| 满18在线观看网站| 男人爽女人下面视频在线观看| 制服丝袜香蕉在线| 不卡av一区二区三区| 久久精品国产亚洲av高清一级| 中文乱码字字幕精品一区二区三区| 51午夜福利影视在线观看| 久久久久精品性色| 欧美精品一区二区免费开放| 伊人久久大香线蕉亚洲五| 精品午夜福利在线看| 久久毛片免费看一区二区三区| 日韩成人av中文字幕在线观看| 国产免费现黄频在线看| www日本在线高清视频| 国产成人精品无人区| 9热在线视频观看99| 久久久亚洲精品成人影院| 亚洲精品第二区| 18在线观看网站| 热re99久久精品国产66热6| 成人黄色视频免费在线看| 亚洲国产成人一精品久久久| 久久久久精品久久久久真实原创| 黑人巨大精品欧美一区二区蜜桃| 国产一区二区在线观看av| 国产野战对白在线观看| 狂野欧美激情性xxxx| 精品一区二区三卡| 夫妻午夜视频| 久久亚洲国产成人精品v| 9色porny在线观看| 久久鲁丝午夜福利片| 少妇猛男粗大的猛烈进出视频| 国产欧美日韩一区二区三区在线| 69精品国产乱码久久久| 成年av动漫网址| 日韩大片免费观看网站| 精品一区二区三卡| 国产人伦9x9x在线观看| 中文字幕人妻丝袜一区二区 |