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

    多級軸流壓氣機(jī)二維性能預(yù)測方法

    2020-04-16 01:14:50楊曉飛蔣永松潘若癡孟德君
    航空學(xué)報 2020年3期
    關(guān)鍵詞:葉型攻角壓氣機(jī)

    楊曉飛,蔣永松,潘若癡,孟德君

    中國航發(fā)沈陽發(fā)動機(jī)研究所,沈陽 110015

    壓氣機(jī)特性計算主要有一維計算、二維計算、三維數(shù)值模擬3種方法。一維計算是在壓氣機(jī)平均流面上進(jìn)行,當(dāng)壓氣機(jī)非設(shè)計點(diǎn)的流動參數(shù)徑向發(fā)生扭曲[1],或流路收縮嚴(yán)重時,其計算精度降低。三維數(shù)值模擬在20世紀(jì)90年代商業(yè)化后迅速普及,并在壓氣機(jī)氣動性能分析中扮演重要角色,尤其是三維定常雷諾平均數(shù)值模擬已是葉輪機(jī)氣動領(lǐng)域不可或缺的計算方法。然而,湍流模型、摻混面模型等所進(jìn)行的簡化會造成計算誤差,邊界條件的不同、幾何細(xì)節(jié)差異也會對計算結(jié)果產(chǎn)生影響[2-3]。隨著壓氣機(jī)級數(shù)增多、負(fù)荷升高,三維數(shù)值模擬的精度下降,性能分析完全依賴三維計算存在巨大的風(fēng)險;除此之外,三維計算花費(fèi)時間長、數(shù)據(jù)量大,在壓氣機(jī)方案設(shè)計階段完全采用三維計算使得設(shè)計周期延長。在三維數(shù)值模擬“浪潮”漸漸回歸平靜后,S2流面正問題計算仍是壓氣機(jī)特性評估及流場分析的重要手段。

    吳仲華先生在1952年提出了兩類流面迭代計算葉輪機(jī)械內(nèi)部流動的理論,基于此理論發(fā)展的流線曲率法或通流矩陣法仍是S2流面求解的主要方法[4-6]。兩種方法都采用定常、絕熱、無黏假設(shè),黏性的影響通過引進(jìn)經(jīng)驗(yàn)或半經(jīng)驗(yàn)的模型來彌補(bǔ)。這些模型必須經(jīng)過大量試驗(yàn)數(shù)據(jù)和設(shè)計經(jīng)驗(yàn)不斷改進(jìn)、校核和驗(yàn)證。隨著壓氣機(jī)設(shè)計技術(shù)的進(jìn)步及對內(nèi)部流動機(jī)理認(rèn)識的深入,Koch和Smith[7]、Konig等[8-9]研究了葉型損失計算方法,Bloch[10]、Boyer和Brien[11]對雙激波模型計算方法做了重要論述。近幾年,Milan[12]、Kim[13]、Azamar Aguirre[14]等持續(xù)對經(jīng)驗(yàn)?zāi)P瓦M(jìn)行了研究;同時,為提升計算精度,適應(yīng)不同理念葉型的損失、落后角預(yù)估,S1流面計算在模型發(fā)展中得到應(yīng)用[15-17]。中國對于經(jīng)驗(yàn)?zāi)P偷难芯恳嘣诓粩嗤七M(jìn),劉波[18-19]、季路成[20]、滕金芳[21]等開展了有關(guān)工作。對于工程設(shè)計,中國的壓氣機(jī)技術(shù)體系是較為完備的,但更新完善卻進(jìn)展緩慢[22]。隨著中國多級高壓壓氣機(jī)的發(fā)展,穩(wěn)定、高精度的S2正問題計算程序需求日益迫切。

    無論采用何種流線曲率方程,其求解過程都是相似的,計算的精度主要受經(jīng)驗(yàn)、半經(jīng)驗(yàn)?zāi)P偷挠绊懀@些模型也正是本研究的重點(diǎn)。本文研究面向工程需要,為了提高計算精度,滿足方案評估要求,在一個流線曲率計算程序的基礎(chǔ)上,通過對比分析,優(yōu)化了攻角、落后角及損失計算方法,尤其改進(jìn)了端壁損失計算方法;同時,為克服基于傳統(tǒng)平面葉柵試驗(yàn)數(shù)據(jù)的攻角、落后角模型與新設(shè)計理念葉型之間的不匹配問題,采用S1流面計算修正S2的方法,最終發(fā)展出了一個高精度的S2正問題計算程序;利用3個不同負(fù)荷水平的、經(jīng)試驗(yàn)驗(yàn)證的多級壓氣機(jī)進(jìn)行了檢驗(yàn)計算,證實(shí)了方法的計算精度和穩(wěn)定性。目前,基于該方法的流線曲率程序已在工程設(shè)計及分析中開始應(yīng)用。

    1 流線曲率程序

    利用流線曲率法求解完全徑向平衡方程是S2流面計算的主要方法之一,在工程中應(yīng)用也最廣泛。一個常見的完全徑向平衡方程為

    (1)

    式中:ρ為密度;p為壓力;V、W分別為絕對和相對速度;r為半徑;rm為子午流線的曲率半徑;φ為流線與軸向的夾角;“m”表示子午方向,“θ”表示圓周方向。式(1)中忽略了葉片力。

    本文研究基于的平臺為中國航發(fā)沈陽發(fā)動機(jī)研究所的流線曲率程序,該程序采用基于軸對稱假設(shè)的流線曲率法求解S2流場,具有包括靜子角度調(diào)節(jié)在內(nèi)的全轉(zhuǎn)速特性計算功能。輸入?yún)?shù)主要包括流路、進(jìn)出口葉型角、稠度、最大相對厚度及位置、最大撓度位置、葉尖間隙以及轉(zhuǎn)速等。因不具備摻混計算,且模型未曾改進(jìn),所以程序?qū)Ξ?dāng)前多級壓氣機(jī)的計算精度較低。該程序較為完整的框架為本文的經(jīng)驗(yàn)、半經(jīng)驗(yàn)?zāi)P脱芯刻峁┝肆己没A(chǔ)。

    2 經(jīng)驗(yàn)?zāi)P?/h2>

    2.1 最小損失攻角、落后角模型

    目前,最小損失攻角的主流計算方法為NASA SP-36[23]中的方法,本文仍沿用這一方法。本程序的最小損失攻角在來流馬赫數(shù)Ma≤1.0時,利用式(2)計算;在Ma≥1.2時,使用葉背相切法計算,見式(3);在1.0

    Ma≤1.0

    (2)

    i*=im2Ma≥1.2

    (3)

    i*=im2+(im1-im2)f(Ma) 1.0

    (4)

    落后角計算模型為

    (5)

    (6)

    (7)

    2.2 設(shè)計點(diǎn)損失模型

    與攻角、落后角模型的匱乏相反,自流線曲率法在工程中應(yīng)用以來,損失模型不斷被提出,葉型損失、激波損失、二次流損失、端壁損失計算模型相對豐富,但每種模型都有使用范圍,不同模型之間的良好匹配是二維正問題程序計算的基礎(chǔ)。本文發(fā)展的設(shè)計點(diǎn)損失計算方法為

    ωt=ωp+ωs+ωsec+ωe

    (8)

    式中:ωt表示總損失;等號右側(cè)4項(xiàng)依次為葉型損失、激波損失、二次流損失和端壁損失。

    2.2.1 葉型損失

    葉型損失由葉片表面附面層和尾跡摻混形成。Lieblein和Roudebush根據(jù)葉柵試驗(yàn)建立了計算葉型損失的基本方法,即根據(jù)擴(kuò)散程度與尾跡動量厚度之間的關(guān)系計算葉型損失[24-25],后續(xù)研究者大都沿用此方法,只是模型的擴(kuò)散程度與尾跡動量厚度之間的曲線略有不同。本文采用Wright和Miller[26]發(fā)展的葉型損失模型,該模型利用葉型損失參數(shù)ωpara與當(dāng)量擴(kuò)散因子Deq之間的關(guān)系曲線計算損失,并考慮了馬赫數(shù)帶來的影響,如圖1所示。當(dāng)量擴(kuò)散因子和葉型損失參數(shù)計算方法分別為

    (9)

    (10)

    圖1 葉型損失參數(shù)關(guān)系曲線Fig.1 Correlation curves for profile loss parameter

    2.2.2 激波損失

    高負(fù)荷多級高壓壓氣機(jī)出于效率和可靠性的要求,一級轉(zhuǎn)子葉尖進(jìn)口相對馬赫數(shù)一般不高于1.4,因此,本文的激波損失采用經(jīng)典的Miller模型計算[27]。

    Miller模型的關(guān)鍵在于如何計算激波前吸力面馬赫數(shù),文獻(xiàn)中大都采用普朗特-邁耶函數(shù)計算,但這不可避免地會遇到如何求解吸力面氣流轉(zhuǎn)角的問題。根據(jù)計算校核,本文中激波前吸力面馬赫數(shù)MaS的計算表達(dá)式為[28]

    (11)

    式中:Δβ為氣流流經(jīng)葉片的轉(zhuǎn)角;Ma1為葉柵進(jìn)口馬赫數(shù)。

    2.2.3 二次流損失

    葉片表面附面層的徑向流動及葉片之間的通道渦會造成氣流的二次流損失,這種二次流動存在于由葉根到葉尖的整個葉片通道。本文利用Griepentrog[29]方法計算二次流損失,二次流損失與總損失之間的關(guān)系為

    (12)

    式中:fs為二次流損失因子,與葉片彎角θ′及展弦比h/c有關(guān),如圖2所示。

    圖2 二次流損失因子曲線Fig.2 Curve for secondary loss factor

    2.2.4 端壁損失

    端壁區(qū)域是壓氣機(jī)內(nèi)部流動最復(fù)雜、流動機(jī)理認(rèn)識最淺薄的區(qū)域。葉尖間隙流動、端壁附面層及其與葉片表面附面層,乃至激波的相互作用,使得此部分區(qū)域的損失很難預(yù)測。一種常用的端壁損失模型認(rèn)為根、尖流面端壁損失為葉型損失的倍數(shù),用端壁損失系數(shù)表示,由端壁附面層引起的附加損失沿葉高呈三次曲線分布[11]。此模型的不足之處在于端壁損失系數(shù)的給定只憑經(jīng)驗(yàn),而端壁區(qū)域損失在總損失的占比較大且在二維正問題計算中熵積累會對計算收斂性產(chǎn)生影響。

    為了確定端壁損失系數(shù)的值,本文利用2個已有的模型,提出了新的端壁損失計算方法。

    (13)

    (14)

    對于轉(zhuǎn)子,式(13)和式(14)中的速度采用相對速度。

    圖3 端壁損失參數(shù)隨擴(kuò)散因子和間隙的變化Fig.3 Variation of endwall loss parameter with diffusion factor and clearance

    已知根、尖流面的端壁損失后,可計算出葉根和葉尖的端壁損失系數(shù)KH和KT:

    Kj=ωe/ωpj=H,T

    (15)

    在求出KH和KT的值以后,由端壁附面層引起的附加損失仍根據(jù)Hearsey方法計算,但附面層影響范圍控制在距葉根、葉尖30%葉高區(qū)域,即端壁損失計算方法為

    (16)

    (17)

    2.3 非設(shè)計點(diǎn)損失模型

    非設(shè)計點(diǎn)損失計算采用Aungier[1]模型,即

    ω=ωt[5-4(ξ+2)]ξ<-2

    (18)

    ω=ωt(1+ξ2) -2≤ξ≤1

    (19)

    ω=2ωtξξ>1

    (20)

    式中:ξ=(i-i*)/WIDTH,WIDTH為可用攻角的范圍。

    可用攻角范圍的計算參考文獻(xiàn)[18]中的方法。同時,考慮到負(fù)攻角狀態(tài)葉柵通道阻塞的影響,對可用攻角范圍進(jìn)行修正。當(dāng)葉柵喉部阻塞時,由進(jìn)口到喉部應(yīng)用連續(xù)方程,經(jīng)變換得到:

    (21)

    式中:o為葉柵喉道寬度;s為柵距;λ為速度系數(shù);k為比熱比;βcr為阻塞狀態(tài)的相對氣流角。

    葉柵負(fù)攻角狀態(tài)的可用攻角范圍應(yīng)滿足以下要求:

    (22)

    根據(jù)計算經(jīng)驗(yàn),由葉柵阻塞引起的可用攻角范圍修正量不應(yīng)超過2.5°。一個轉(zhuǎn)子葉片修正前后的各流面負(fù)攻角狀態(tài)可用攻角范圍及與S1計算結(jié)果的對比如圖4所示,可見修正后與S1計算結(jié)果吻合得更好。

    圖4 可用攻角范圍對比Fig.4 Comparison of available angle of incidence range

    2.4 徑向摻混

    為了降低端壁區(qū)域的低能量氣流堆積,提高程序的計算穩(wěn)定性,在S2正問題程序中有必要進(jìn)行摻混計算,這也符合壓氣機(jī)內(nèi)部的流動規(guī)律。本文采用Aungier[1]的方法模擬計算內(nèi)部摻混流動。

    3 S1流面計算的應(yīng)用

    隨著壓氣機(jī)設(shè)計技術(shù)的發(fā)展,不斷涌現(xiàn)出新的葉型,比如多圓弧、控制擴(kuò)散、定制葉型等,但公開文獻(xiàn)中針對這些新設(shè)計理念葉型的最小損失攻角與落后角模型尚未涉及,不同的造型方法也給解析形式的經(jīng)驗(yàn)?zāi)P桶l(fā)展帶來困難。圖5和圖6分別為NACA-65葉型和定制葉型利用經(jīng)驗(yàn)關(guān)系式計算的最小損失攻角和落后角與S1流面計算結(jié)果的對比,其中,最小損失攻角采用式(2)計算,落后角采用式(6)計算,并利用式(7)計算密流比的影響。對于NACA-65葉型,模型計算的最小損失攻角與S1流面計算結(jié)果相當(dāng),如圖5(a)所示,落后角比S1結(jié)果偏高0.5°,相對偏高8.6%,如圖5(b)所示;而對于定制葉型,模型計算的最小損失攻角偏正,如圖6(a),落后角偏高1.8°,相對偏高37.6%左右,如圖6(b)所示。最小損失攻角、落后角計算方法對NACA-65葉型計算精度相對較高,對新理念葉型的計算精度相對較低。

    圖5 NACA-65葉型計算結(jié)果Fig.5 Calculation results of NACA-65 blade profile

    圖6 定制葉型計算結(jié)果Fig.6 Calculation results of tailored blade profile

    為了解決豐富的葉型類型與相對匱乏的最小損失攻角、落后角模型之間的不匹配問題,引入了S1流面計算。事實(shí)上,利用數(shù)值計算已成為當(dāng)前改進(jìn)經(jīng)驗(yàn)?zāi)P偷氖侄沃?,例如,Schnoes和Nicke[17]利用葉柵數(shù)值模擬重新校驗(yàn)了傳統(tǒng)的落后角模型、損失模型的系數(shù)。本文提出的利用S1流面計算修正S2的計算方法僅限于對最小損失攻角和對應(yīng)狀態(tài)落后角的修正。修正方法為

    (23)

    (24)

    (25)

    式中:δ*為最小損失狀態(tài)的落后角;下標(biāo)“M”“H”“T”分別代表S2計算中的葉中、葉根和葉尖流面;“S1”表示S1流面計算結(jié)果;C1、C2為S1流面計算與實(shí)際流動之間的修正量,即流動的三維效應(yīng)對落后角的影響。根據(jù)與三維計算、試驗(yàn)特性的比較,C1、C2取值在[0.5°,2.5°]。

    圖7為利用S1修正S2計算的一個多級壓氣機(jī)的跨聲速轉(zhuǎn)子、亞聲速轉(zhuǎn)子和一排靜子葉片的氣流角分布及其與三維計算的對比,該壓氣機(jī)采用定制葉型。對比的狀態(tài)點(diǎn)均在計算的最高效率或總壓恢復(fù)系數(shù)附近。對比表明,在10%~90%相對葉高之間的主流區(qū),S2計算的出口氣流角與三維計算結(jié)果吻合較好,尤其是亞聲速轉(zhuǎn)子和靜子葉片??缏曀俎D(zhuǎn)子和靜子的進(jìn)口氣流角與三維計算的展向分布不同,但狀態(tài)點(diǎn)在最小損失點(diǎn)附近,因此,由攻角偏離最小損失攻角引起的落后角變化量較小。在近端壁區(qū)域,三維計算考慮了黏性的影響,壁面采用無滑移邊界;而S2計算出于收斂性的考慮,會避免熵增在端壁區(qū)域的堆積,近壁面的子午速度相對主流區(qū)不會迅速降低,所以,S2的計算結(jié)果與三維計算在端壁區(qū)有所不同。

    圖7 跨聲速轉(zhuǎn)子、亞聲速轉(zhuǎn)子和靜子的氣流角對比Fig.7 Comparison of airflow angles of transonic rotor, subsonic rotor and stator

    在目前攻角模型、落后角模型匱乏的狀況下,該方法解決了傳統(tǒng)的攻角、落后角模型與種類繁多、不同設(shè)計理念葉型之間的適用性問題,既不受葉型類型的約束,又保留了S2正問題程序中馬赫數(shù)、葉尖間隙等因素對最小損失攻角的影響以及攻角變化對落后角的影響。該方法及不同模型的聯(lián)合實(shí)施有效提升了二維正問題程序的精度和可靠性。

    4 失速模型

    本研究采用最大總壓升準(zhǔn)則來預(yù)估多級壓氣機(jī)的穩(wěn)定工作邊界。對任意一條等轉(zhuǎn)速線,當(dāng)壓氣機(jī)流量減小后,計算總壓比不再升高,則認(rèn)為壓氣機(jī)發(fā)生失速。

    為了驗(yàn)證失速模型,統(tǒng)計了4臺多級壓氣機(jī)在設(shè)計轉(zhuǎn)速和0.8相對換算轉(zhuǎn)速的試驗(yàn)流量-壓比特性。4臺壓氣機(jī)依次用C1、C2、C3、C4表示,其平均級壓比如圖8所示,利用堵點(diǎn)流量和失速點(diǎn)壓比作歸一化處理后的特性如圖9所示。由圖9可知,設(shè)計轉(zhuǎn)速特性線的最大壓比點(diǎn)即為失速點(diǎn);在0.8相對換算轉(zhuǎn)速下,C1、C2、C3壓氣機(jī)的等轉(zhuǎn)速線均在最高壓比點(diǎn)發(fā)生失速,而C4壓氣機(jī)的等轉(zhuǎn)速線在經(jīng)過最高壓比點(diǎn)(A點(diǎn))后,試驗(yàn)仍錄取一個狀態(tài)點(diǎn);經(jīng)對比,由A點(diǎn)計算的喘振裕度比失速點(diǎn)計算的喘振裕度僅偏低0.5%。對于多級壓氣機(jī)的喘振裕度而言,數(shù)值計算產(chǎn)生0.5%的偏差是可接受的。

    本研究以計算的最高壓比點(diǎn)作為失速點(diǎn),下文的計算證實(shí)此種方法是有效的,預(yù)測的多級壓氣機(jī)中高轉(zhuǎn)速失速邊界與試驗(yàn)結(jié)果吻合良好。

    圖8 4臺壓氣機(jī)的平均級壓比Fig.8 Average stage pressure ratio of four compressors

    圖9 歸一化試驗(yàn)特性Fig.9 Normalized test characteristics

    5 校驗(yàn)計算

    采用3個不同負(fù)荷水平的、經(jīng)試驗(yàn)驗(yàn)證的壓氣機(jī)對發(fā)展的方法進(jìn)行校核計算。壓氣機(jī)的負(fù)荷水平用葉尖負(fù)荷系數(shù)表示,計算表達(dá)式為

    (26)

    式中:H為負(fù)荷系數(shù);Lu為壓縮功;UT為葉片葉尖旋轉(zhuǎn)速度。3個壓氣機(jī)(依次用A、B、C表示)的平均負(fù)荷系數(shù)分布如圖10所示。

    圖10 壓氣機(jī)負(fù)荷系數(shù)Fig.10 Work coefficient of compressors

    為了驗(yàn)證方法的精度及穩(wěn)定性,首先利用改進(jìn)前的程序計算了A壓氣機(jī)特性,并與改進(jìn)后程序的計算結(jié)果進(jìn)行了對比。對于B、C兩臺壓氣機(jī),因?yàn)榭倝罕?、?fù)荷提高,改進(jìn)前的程序未能計算收斂。此外,針對C壓氣機(jī)開展了級間參數(shù)對比,以驗(yàn)證級間匹配狀況。

    5.1 A壓氣機(jī)

    A壓氣機(jī)為面向燃?xì)廨啓C(jī)的低壓壓氣機(jī),壓氣機(jī)前兩級轉(zhuǎn)子存在局部超聲速區(qū),一級轉(zhuǎn)子葉尖進(jìn)口相對馬赫數(shù)為1.16。壓氣機(jī)轉(zhuǎn)、靜子葉片均采用修正圓弧中線、NACA-65葉厚分布。為適應(yīng)不同葉片的進(jìn)口馬赫數(shù)及流通能力,中弧線采用了不同程度的加彎或減彎。

    改進(jìn)前后的計算特性及其與試驗(yàn)的對比如圖11所示,其中試驗(yàn)相對換算轉(zhuǎn)速nr為1.0和0.9時未錄取到喘振邊界。對比表明,在1.0和0.9轉(zhuǎn)速程序改進(jìn)前后計算的流量-壓比、流量-效率特性基本重合;轉(zhuǎn)速降低后,改進(jìn)前程序計算的效率降低,計算的流量-壓比線在改進(jìn)后計算特性線的左下方;值得注意的是,改進(jìn)前的程序計算的穩(wěn)定邊界水平較低,本研究發(fā)展的計算方法顯著提升了程序的計算穩(wěn)定性。

    圖11 A壓氣機(jī)特性Fig.11 Overall performance of Compressor A

    改進(jìn)后程序計算的堵點(diǎn)流量比試驗(yàn)偏低4%左右,最高壓比水平與試驗(yàn)相當(dāng);各轉(zhuǎn)速計算的峰值效率與試驗(yàn)相當(dāng),從0.8和0.75轉(zhuǎn)速來看,計算穩(wěn)定邊界與試驗(yàn)較為接近。

    5.2 B壓氣機(jī)

    B壓氣機(jī)為E3十級高壓壓氣機(jī)第2輪修改后的壓氣機(jī)。E3十級高壓壓氣機(jī)在1.0換算轉(zhuǎn)速的設(shè)計流量為54.4 kg/s,總壓比為25.0,效率為0.847[30],一級轉(zhuǎn)子葉尖進(jìn)口相對馬赫數(shù)為1.353。

    壓氣機(jī)前4級轉(zhuǎn)子為跨聲速葉排,采用定制流面葉型的方法設(shè)計。5~10級轉(zhuǎn)子采用修正圓弧中線,5、6級轉(zhuǎn)子采用多圓弧厚度分布,7~10級采用NACA-65厚度分布。1~4級靜子具有修正的圓弧中線和多圓弧厚度分布,5~9級靜子采用NACA-65葉型。最后一級靜子負(fù)責(zé)將氣流轉(zhuǎn)回軸向,彎角較大,采用了特殊設(shè)計的葉型。為保證喘振裕度,前6級靜子可調(diào)。

    計算的壓氣機(jī)0.9相對換算轉(zhuǎn)速以上特性及與試驗(yàn)的對比如圖12所示。對比可得,計算的流量-壓比特性與試驗(yàn)吻合良好,計算效率比試驗(yàn)略高,計算穩(wěn)定邊界與試驗(yàn)基本一致。1.0與1.025相對換算轉(zhuǎn)速的計算效率偏高,這一方面是由于轉(zhuǎn)子葉尖進(jìn)口相對馬赫數(shù)較高,激波模型計算精度降低,另一方面,文獻(xiàn)[30]認(rèn)為進(jìn)口測量儀表及多次高速失速引起的性能退化使測試效率偏低。

    圖12 B壓氣機(jī)(E3十級壓氣機(jī))特性Fig.12 Overall performance of Compressor B (E3 ten-stage compressor)

    0.9相對轉(zhuǎn)速計算與試驗(yàn)偏差較大,這除了計算給定的間隙較小外,還與進(jìn)口導(dǎo)流葉片的性能計算有關(guān),轉(zhuǎn)速降低后進(jìn)口導(dǎo)流葉片關(guān)角度值很大,出口氣流角及總壓恢復(fù)系數(shù)計算精度降低。

    5.3 C壓氣機(jī)

    C壓氣機(jī)為高負(fù)荷多級壓氣機(jī),前3級轉(zhuǎn)子葉片存在超聲速區(qū)域,一級轉(zhuǎn)子葉尖進(jìn)口相對馬赫數(shù)為1.22。各排葉片均采用多圓弧葉型。在特性計算前,首先計算了各排葉片根、中、尖3個流面的葉型攻角特性,對最小損失攻角和落后角進(jìn)行了修正,隨后開展了全轉(zhuǎn)速特性預(yù)估。計算特性如圖13所示,計算與試驗(yàn)吻合良好。

    試驗(yàn)對前5級轉(zhuǎn)子后的總壓和總溫進(jìn)行了測量。設(shè)計轉(zhuǎn)速由堵點(diǎn)到失速點(diǎn)的級間壓比、溫升比分布的計算結(jié)果與試驗(yàn)對比如圖14~圖16所示,圖中的壓比、溫升比利用試驗(yàn)出口總壓比、溫升比進(jìn)行了無量綱處理。對比表明,堵點(diǎn)、工作點(diǎn)的計算壓比和溫升比與試驗(yàn)基本吻合。在失速點(diǎn),計算的溫升比更低,這與圖13中的總特性相應(yīng),在設(shè)計轉(zhuǎn)速近失速點(diǎn),當(dāng)計算與試驗(yàn)壓比相當(dāng)時,計算的效率更高。

    圖13 C壓氣機(jī)特性Fig.13 Overall performance of Compressor C

    圖14 近堵塞點(diǎn)的壓力和溫度對比Fig.14 Comparison of pressure and temperature near choked operating point

    圖15 近設(shè)計點(diǎn)的壓力和溫度對比Fig.15 Comparison of pressure and temperature near design operating point

    圖16 近失速點(diǎn)的壓力和溫度對比Fig.16 Comparison of pressure and temperature near stall operating point

    6 結(jié) 論

    二維計算仍是多級壓氣機(jī)設(shè)計與分析的核心,經(jīng)驗(yàn)、半經(jīng)驗(yàn)?zāi)P透桦S壓氣機(jī)設(shè)計技術(shù)的進(jìn)步而不斷更新發(fā)展。本文通過對比分析,優(yōu)化了攻角、落后角及損失計算模塊,提升了經(jīng)驗(yàn)、半經(jīng)驗(yàn)計算模塊的精度,并提出了利用S1流面修正S2正問題計算的方法,解決了攻角、落后角模型與新理念葉型之間不匹配的問題。最后,利用3個經(jīng)試驗(yàn)驗(yàn)證的多級壓氣機(jī)對發(fā)展的計算方法進(jìn)行校驗(yàn)。結(jié)果表明:

    1) 除低負(fù)荷壓氣機(jī)的計算流量偏低外,采用本研究發(fā)展的方法計算的壓氣機(jī)性能與試驗(yàn)結(jié)果吻合良好,展現(xiàn)出相對較高的計算精度。

    2) 本文發(fā)展的方法合理,較大幅度提高了程序的計算穩(wěn)定性,可用于多級軸流壓氣機(jī)性能分析,為壓氣機(jī)多程序校核提供了新手段。

    猜你喜歡
    葉型攻角壓氣機(jī)
    葉片前緣對吸力面邊界層3維流動影響分析
    先進(jìn)動葉平面葉柵試驗(yàn)和數(shù)值分析
    楔形葉片旋轉(zhuǎn)空化器葉型改進(jìn)數(shù)值模擬研究
    軸流壓氣機(jī)效率評定方法
    重型燃?xì)廨啓C(jī)壓氣機(jī)第一級轉(zhuǎn)子葉片斷裂分析
    壓氣機(jī)緊湊S形過渡段內(nèi)周向彎靜子性能數(shù)值計算
    風(fēng)標(biāo)式攻角傳感器在超聲速飛行運(yùn)載火箭中的應(yīng)用研究
    大攻角狀態(tài)壓氣機(jī)分離流及葉片動力響應(yīng)特性
    附加攻角效應(yīng)對顫振穩(wěn)定性能影響
    振動與沖擊(2015年2期)2015-05-16 05:37:34
    民用飛機(jī)攻角傳感器安裝定位研究
    欧美大码av| 亚洲第一av免费看| 国产精品免费一区二区三区在线| 啪啪无遮挡十八禁网站| 午夜福利成人在线免费观看| 女性被躁到高潮视频| 伊人久久大香线蕉亚洲五| 深夜精品福利| 一本久久中文字幕| 午夜成年电影在线免费观看| 亚洲va日本ⅴa欧美va伊人久久| 乱人伦中国视频| 国产黄a三级三级三级人| 91国产中文字幕| 高清在线国产一区| 欧美不卡视频在线免费观看 | 中文字幕久久专区| 国产一区二区三区视频了| 丰满的人妻完整版| 大型av网站在线播放| 两人在一起打扑克的视频| 国产激情欧美一区二区| 999精品在线视频| 亚洲激情在线av| 一个人观看的视频www高清免费观看 | 九色国产91popny在线| 欧美乱色亚洲激情| 日韩 欧美 亚洲 中文字幕| 90打野战视频偷拍视频| 很黄的视频免费| 国产精品野战在线观看| 亚洲精品在线观看二区| 亚洲九九香蕉| 老司机福利观看| 久久久久久久久中文| 纯流量卡能插随身wifi吗| 男女床上黄色一级片免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 黄色毛片三级朝国网站| 午夜福利高清视频| 欧美日韩中文字幕国产精品一区二区三区 | 久久午夜综合久久蜜桃| 亚洲熟女毛片儿| 一进一出好大好爽视频| 91国产中文字幕| 久久人人爽av亚洲精品天堂| 老熟妇仑乱视频hdxx| 亚洲精品美女久久久久99蜜臀| 这个男人来自地球电影免费观看| 黄片播放在线免费| 一区二区三区高清视频在线| 女生性感内裤真人,穿戴方法视频| 如日韩欧美国产精品一区二区三区| 精品第一国产精品| 国产高清视频在线播放一区| 欧美日韩中文字幕国产精品一区二区三区 | 人妻丰满熟妇av一区二区三区| 国产亚洲精品久久久久久毛片| 国产精品久久久久久精品电影 | 黄色丝袜av网址大全| 成年版毛片免费区| 国产高清激情床上av| 色播亚洲综合网| 中文字幕色久视频| 亚洲av五月六月丁香网| 欧美激情 高清一区二区三区| 中国美女看黄片| 国内久久婷婷六月综合欲色啪| 老司机在亚洲福利影院| 国产精品免费一区二区三区在线| 两性夫妻黄色片| 国产在线精品亚洲第一网站| 好男人电影高清在线观看| 国产精品九九99| 精品熟女少妇八av免费久了| 国产亚洲精品av在线| 亚洲伊人色综图| 嫩草影院精品99| e午夜精品久久久久久久| 在线观看舔阴道视频| 人人妻,人人澡人人爽秒播| 一个人免费在线观看的高清视频| 久久草成人影院| 欧美 亚洲 国产 日韩一| 亚洲国产精品999在线| 亚洲精品中文字幕一二三四区| 成年女人毛片免费观看观看9| 午夜福利影视在线免费观看| 午夜a级毛片| 啪啪无遮挡十八禁网站| 夜夜躁狠狠躁天天躁| 亚洲精品中文字幕在线视频| 一个人免费在线观看的高清视频| 成人亚洲精品av一区二区| 国产精品99久久99久久久不卡| 精品国产一区二区三区四区第35| 嫁个100分男人电影在线观看| 久久国产精品人妻蜜桃| 18禁观看日本| 狠狠狠狠99中文字幕| 大码成人一级视频| 国产乱人伦免费视频| 一区二区日韩欧美中文字幕| 欧美色视频一区免费| 日日爽夜夜爽网站| 亚洲精品在线观看二区| 电影成人av| 国产片内射在线| 一边摸一边做爽爽视频免费| 亚洲精品美女久久久久99蜜臀| 免费在线观看黄色视频的| 波多野结衣巨乳人妻| 免费看a级黄色片| av视频免费观看在线观看| 电影成人av| 伦理电影免费视频| 久久精品91蜜桃| 老汉色av国产亚洲站长工具| 欧美中文日本在线观看视频| 国产精品亚洲一级av第二区| 午夜福利免费观看在线| 亚洲国产精品久久男人天堂| 成人av一区二区三区在线看| 91av网站免费观看| 国产精品 欧美亚洲| 欧美成人午夜精品| or卡值多少钱| 国产亚洲欧美精品永久| 一二三四在线观看免费中文在| 国产欧美日韩一区二区三| 亚洲人成电影观看| 99国产精品99久久久久| 在线观看免费视频网站a站| 又大又爽又粗| 亚洲午夜理论影院| 久久精品国产亚洲av香蕉五月| 久久久久国内视频| 国产成人精品久久二区二区免费| 国产精品香港三级国产av潘金莲| 国产精品久久久久久人妻精品电影| 国产精品 欧美亚洲| 久久久久久免费高清国产稀缺| 国产精品九九99| 性欧美人与动物交配| 手机成人av网站| 91麻豆精品激情在线观看国产| 成年人黄色毛片网站| 一级毛片高清免费大全| 日本 av在线| 黄片大片在线免费观看| 在线观看舔阴道视频| 日韩精品免费视频一区二区三区| 久久久国产成人免费| 精品卡一卡二卡四卡免费| 啦啦啦免费观看视频1| 中国美女看黄片| 可以在线观看毛片的网站| 最好的美女福利视频网| 91国产中文字幕| 亚洲中文日韩欧美视频| 在线永久观看黄色视频| 日本vs欧美在线观看视频| 激情在线观看视频在线高清| 欧美激情极品国产一区二区三区| 国产精品亚洲av一区麻豆| 777久久人妻少妇嫩草av网站| 国产精品一区二区三区四区久久 | 日韩欧美在线二视频| 在线观看www视频免费| 视频区欧美日本亚洲| 成年女人毛片免费观看观看9| 亚洲五月色婷婷综合| 亚洲精品美女久久av网站| 亚洲成a人片在线一区二区| 老熟妇乱子伦视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 国产亚洲精品综合一区在线观看 | 人人妻人人澡欧美一区二区 | 麻豆久久精品国产亚洲av| 亚洲人成77777在线视频| 亚洲男人的天堂狠狠| 国产99白浆流出| 99国产综合亚洲精品| 国产精品 欧美亚洲| 少妇粗大呻吟视频| 一个人观看的视频www高清免费观看 | 成人国产一区最新在线观看| 午夜福利在线观看吧| 看免费av毛片| 男女下面进入的视频免费午夜 | 夜夜夜夜夜久久久久| 国产精品国产高清国产av| 亚洲成av片中文字幕在线观看| e午夜精品久久久久久久| 侵犯人妻中文字幕一二三四区| 在线天堂中文资源库| 很黄的视频免费| 一级a爱视频在线免费观看| 琪琪午夜伦伦电影理论片6080| 校园春色视频在线观看| 亚洲伊人色综图| 在线观看午夜福利视频| 欧美精品啪啪一区二区三区| 久久久国产精品麻豆| 亚洲一区高清亚洲精品| 国产精品综合久久久久久久免费 | 精品国产国语对白av| 一区二区日韩欧美中文字幕| 国产区一区二久久| 精品少妇一区二区三区视频日本电影| 国产一区二区激情短视频| 久久久精品国产亚洲av高清涩受| 精品第一国产精品| 国产精品美女特级片免费视频播放器 | 亚洲熟妇中文字幕五十中出| 青草久久国产| 亚洲 欧美一区二区三区| 涩涩av久久男人的天堂| 国产伦人伦偷精品视频| 99久久精品国产亚洲精品| 波多野结衣巨乳人妻| 少妇粗大呻吟视频| 一a级毛片在线观看| 久久久久久国产a免费观看| 欧美乱妇无乱码| 国产在线精品亚洲第一网站| 不卡一级毛片| 99精品欧美一区二区三区四区| 久久精品成人免费网站| 国产av一区在线观看免费| 亚洲欧美精品综合一区二区三区| 婷婷六月久久综合丁香| 国产精品久久视频播放| 中文字幕最新亚洲高清| 亚洲欧美激情在线| 88av欧美| 久9热在线精品视频| 啦啦啦观看免费观看视频高清 | 90打野战视频偷拍视频| 18禁裸乳无遮挡免费网站照片 | 俄罗斯特黄特色一大片| 日韩有码中文字幕| 免费在线观看日本一区| 一级,二级,三级黄色视频| 十分钟在线观看高清视频www| 国产精品亚洲美女久久久| 亚洲欧美日韩无卡精品| 亚洲一码二码三码区别大吗| 欧洲精品卡2卡3卡4卡5卡区| 国产视频一区二区在线看| 亚洲av美国av| 黄色片一级片一级黄色片| 757午夜福利合集在线观看| 日韩精品免费视频一区二区三区| 色精品久久人妻99蜜桃| 熟女少妇亚洲综合色aaa.| 成年版毛片免费区| 女人高潮潮喷娇喘18禁视频| 变态另类丝袜制服| 最近最新免费中文字幕在线| 国产麻豆69| 看黄色毛片网站| 免费看a级黄色片| 久久精品国产亚洲av高清一级| 婷婷丁香在线五月| 成年版毛片免费区| av欧美777| 亚洲五月天丁香| 久久人妻av系列| 久久精品91蜜桃| 午夜福利,免费看| 99精品久久久久人妻精品| 91精品国产国语对白视频| 亚洲专区中文字幕在线| 一二三四社区在线视频社区8| 日本免费a在线| 午夜精品久久久久久毛片777| 黄片播放在线免费| 国产99白浆流出| 亚洲av日韩精品久久久久久密| 国产成人av教育| 香蕉久久夜色| 美女大奶头视频| 一级,二级,三级黄色视频| 黄色a级毛片大全视频| 大型黄色视频在线免费观看| 大陆偷拍与自拍| 十分钟在线观看高清视频www| 97人妻天天添夜夜摸| 亚洲成av片中文字幕在线观看| 国产麻豆成人av免费视频| 少妇熟女aⅴ在线视频| 欧美激情 高清一区二区三区| 99久久国产精品久久久| 一本大道久久a久久精品| 国产精品久久久久久亚洲av鲁大| 久久久久久国产a免费观看| 亚洲色图 男人天堂 中文字幕| 国产片内射在线| 亚洲五月色婷婷综合| 国产午夜精品久久久久久| 99香蕉大伊视频| 亚洲国产日韩欧美精品在线观看 | 午夜日韩欧美国产| 在线观看一区二区三区| 电影成人av| 国产成人啪精品午夜网站| 亚洲欧美日韩另类电影网站| 欧美激情 高清一区二区三区| 99香蕉大伊视频| 女人爽到高潮嗷嗷叫在线视频| 国产高清激情床上av| 久9热在线精品视频| 久久精品亚洲熟妇少妇任你| 黄片小视频在线播放| 可以在线观看的亚洲视频| 又黄又爽又免费观看的视频| 久久精品国产亚洲av香蕉五月| 在线观看午夜福利视频| 欧美一级a爱片免费观看看 | videosex国产| 亚洲av电影不卡..在线观看| 亚洲三区欧美一区| 国产精品香港三级国产av潘金莲| 两性夫妻黄色片| 精品久久久精品久久久| 精品国产美女av久久久久小说| 18禁黄网站禁片午夜丰满| 岛国在线观看网站| 宅男免费午夜| 国产麻豆成人av免费视频| 老熟妇乱子伦视频在线观看| 久久精品影院6| 99久久综合精品五月天人人| 亚洲精品中文字幕在线视频| 国产亚洲精品久久久久5区| 在线播放国产精品三级| 中文亚洲av片在线观看爽| 99国产精品免费福利视频| 一区二区三区精品91| 精品久久久久久久久久免费视频| 韩国精品一区二区三区| 日本a在线网址| 久久久久久久久免费视频了| 韩国精品一区二区三区| 大型av网站在线播放| 国产高清videossex| 老熟妇仑乱视频hdxx| 一个人观看的视频www高清免费观看 | 欧美在线一区亚洲| 精品国内亚洲2022精品成人| 男女下面进入的视频免费午夜 | 日日干狠狠操夜夜爽| 国产av一区在线观看免费| 亚洲 欧美 日韩 在线 免费| 成人免费观看视频高清| 性少妇av在线| 亚洲国产精品合色在线| 亚洲熟妇中文字幕五十中出| 欧美成人免费av一区二区三区| 国产一级毛片七仙女欲春2 | 久久久久久久精品吃奶| 我的亚洲天堂| 涩涩av久久男人的天堂| 中国美女看黄片| 黑人操中国人逼视频| 欧美午夜高清在线| 亚洲人成电影免费在线| 亚洲午夜精品一区,二区,三区| 欧美日韩中文字幕国产精品一区二区三区 | 精品欧美一区二区三区在线| 成人特级黄色片久久久久久久| 国产一区二区在线av高清观看| 一进一出抽搐动态| 国产精品美女特级片免费视频播放器 | 成在线人永久免费视频| 咕卡用的链子| 两个人视频免费观看高清| 亚洲自偷自拍图片 自拍| 国产一区二区三区综合在线观看| 久久久国产成人精品二区| 午夜激情av网站| 琪琪午夜伦伦电影理论片6080| 精品欧美一区二区三区在线| 久热爱精品视频在线9| 天天添夜夜摸| 午夜福利影视在线免费观看| 大型av网站在线播放| 一边摸一边抽搐一进一出视频| 日日干狠狠操夜夜爽| 亚洲国产精品合色在线| 国产精品久久视频播放| 999久久久精品免费观看国产| 后天国语完整版免费观看| 妹子高潮喷水视频| 69精品国产乱码久久久| 精品人妻1区二区| 日韩有码中文字幕| 久久香蕉精品热| 亚洲成av片中文字幕在线观看| 男人操女人黄网站| 搡老熟女国产l中国老女人| 欧美黑人欧美精品刺激| 制服人妻中文乱码| 国产精品一区二区三区四区久久 | 成人永久免费在线观看视频| 亚洲aⅴ乱码一区二区在线播放 | 成人国产综合亚洲| 亚洲人成电影观看| 午夜视频精品福利| 国产麻豆69| 色哟哟哟哟哟哟| 性色av乱码一区二区三区2| 国内久久婷婷六月综合欲色啪| 中文亚洲av片在线观看爽| 十八禁网站免费在线| 亚洲国产精品久久男人天堂| 人人妻人人爽人人添夜夜欢视频| 精品国内亚洲2022精品成人| 亚洲性夜色夜夜综合| 欧美乱码精品一区二区三区| 亚洲国产精品成人综合色| 又紧又爽又黄一区二区| 如日韩欧美国产精品一区二区三区| 久久久久久人人人人人| 国产97色在线日韩免费| 亚洲片人在线观看| 级片在线观看| 天堂√8在线中文| 淫妇啪啪啪对白视频| 一a级毛片在线观看| 午夜精品在线福利| 国产亚洲精品一区二区www| 国产av一区二区精品久久| 一区二区三区国产精品乱码| 一区福利在线观看| 三级毛片av免费| 亚洲国产高清在线一区二区三 | avwww免费| 最新在线观看一区二区三区| 免费在线观看影片大全网站| 99riav亚洲国产免费| 少妇 在线观看| 香蕉丝袜av| 国产aⅴ精品一区二区三区波| 村上凉子中文字幕在线| 久久九九热精品免费| 搡老熟女国产l中国老女人| 亚洲男人天堂网一区| 日本三级黄在线观看| 搡老岳熟女国产| 97超级碰碰碰精品色视频在线观看| 悠悠久久av| 桃红色精品国产亚洲av| av免费在线观看网站| 亚洲精品中文字幕一二三四区| 国产亚洲精品av在线| 美女高潮到喷水免费观看| 精品国产乱子伦一区二区三区| 日韩欧美一区二区三区在线观看| 成年女人毛片免费观看观看9| 精品人妻1区二区| 18禁黄网站禁片午夜丰满| 亚洲国产欧美日韩在线播放| 这个男人来自地球电影免费观看| 国产精品98久久久久久宅男小说| 51午夜福利影视在线观看| 亚洲专区字幕在线| 久久精品人人爽人人爽视色| 人人澡人人妻人| 在线永久观看黄色视频| 国产成人欧美在线观看| 国产高清视频在线播放一区| 日本 欧美在线| 熟妇人妻久久中文字幕3abv| or卡值多少钱| 一区二区三区国产精品乱码| 在线天堂中文资源库| 亚洲欧美日韩高清在线视频| 在线免费观看的www视频| 国产精品永久免费网站| avwww免费| 精品久久久久久成人av| 久久国产精品男人的天堂亚洲| 操美女的视频在线观看| 国产精品影院久久| 非洲黑人性xxxx精品又粗又长| 黄色 视频免费看| 国产av又大| 久久久久久免费高清国产稀缺| 国产成人影院久久av| 一边摸一边做爽爽视频免费| 久久久久国内视频| 女同久久另类99精品国产91| av超薄肉色丝袜交足视频| 又黄又爽又免费观看的视频| 香蕉丝袜av| 在线观看舔阴道视频| 免费高清视频大片| √禁漫天堂资源中文www| 在线观看免费视频网站a站| 我的亚洲天堂| 性色av乱码一区二区三区2| 亚洲成人免费电影在线观看| 97超级碰碰碰精品色视频在线观看| 一区二区日韩欧美中文字幕| 一级毛片精品| 人成视频在线观看免费观看| 老熟妇仑乱视频hdxx| 91麻豆av在线| АⅤ资源中文在线天堂| 色婷婷久久久亚洲欧美| 欧美色视频一区免费| 国产成人系列免费观看| 欧美一级a爱片免费观看看 | 好看av亚洲va欧美ⅴa在| 宅男免费午夜| √禁漫天堂资源中文www| 精品国产一区二区三区四区第35| 欧美一级a爱片免费观看看 | 国产精品爽爽va在线观看网站 | 日本免费一区二区三区高清不卡 | 久久精品影院6| 不卡av一区二区三区| 91国产中文字幕| 91成人精品电影| 亚洲欧洲精品一区二区精品久久久| 成人免费观看视频高清| 日日摸夜夜添夜夜添小说| 欧美激情极品国产一区二区三区| 淫妇啪啪啪对白视频| 在线免费观看的www视频| 法律面前人人平等表现在哪些方面| 欧美在线黄色| 自拍欧美九色日韩亚洲蝌蚪91| 国产又色又爽无遮挡免费看| 免费久久久久久久精品成人欧美视频| 久久精品aⅴ一区二区三区四区| 国内久久婷婷六月综合欲色啪| 国产精品久久久久久精品电影 | 欧美国产精品va在线观看不卡| 黄频高清免费视频| 女人爽到高潮嗷嗷叫在线视频| 日韩视频一区二区在线观看| 色综合婷婷激情| 淫秽高清视频在线观看| 麻豆一二三区av精品| av有码第一页| 亚洲av片天天在线观看| 黄色视频,在线免费观看| 久久精品亚洲熟妇少妇任你| 免费人成视频x8x8入口观看| 国语自产精品视频在线第100页| 老鸭窝网址在线观看| 午夜两性在线视频| 热99re8久久精品国产| 丰满的人妻完整版| 91九色精品人成在线观看| 一区二区三区精品91| 欧美av亚洲av综合av国产av| 国产99久久九九免费精品| 国产亚洲精品久久久久5区| 9热在线视频观看99| 欧美成人一区二区免费高清观看 | 韩国精品一区二区三区| 操出白浆在线播放| 亚洲国产精品999在线| 啦啦啦观看免费观看视频高清 | 国产99白浆流出| 欧美成人免费av一区二区三区| 久久香蕉国产精品| 一级黄色大片毛片| 日本vs欧美在线观看视频| 亚洲av美国av| 欧美国产精品va在线观看不卡| 在线观看日韩欧美| 久久婷婷人人爽人人干人人爱 | 满18在线观看网站| 好男人电影高清在线观看| 级片在线观看| 国产黄a三级三级三级人| 麻豆一二三区av精品| 久久久精品欧美日韩精品| 在线观看舔阴道视频| 免费观看人在逋| 一边摸一边抽搐一进一出视频| 亚洲aⅴ乱码一区二区在线播放 | 国产成人影院久久av| 女人被躁到高潮嗷嗷叫费观| 日本撒尿小便嘘嘘汇集6| 国产真人三级小视频在线观看| 亚洲欧洲精品一区二区精品久久久| 长腿黑丝高跟| 国产免费男女视频| 亚洲精品久久国产高清桃花| 每晚都被弄得嗷嗷叫到高潮| 国产麻豆69| 中文字幕人成人乱码亚洲影| 久久久久精品国产欧美久久久| 国产主播在线观看一区二区| av天堂久久9| 色综合站精品国产| 亚洲一区二区三区不卡视频| 日本 欧美在线| 女人被躁到高潮嗷嗷叫费观| 中文字幕色久视频| 午夜免费观看网址| 一级毛片高清免费大全| av视频在线观看入口| 在线观看www视频免费| x7x7x7水蜜桃| 国产一区在线观看成人免费| 黄片播放在线免费| 母亲3免费完整高清在线观看| 女人高潮潮喷娇喘18禁视频|