• <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ī)攻角傳感器安裝定位研究
    国产一区二区三区在线臀色熟女| 国产一区二区三区视频了| 又黄又爽又免费观看的视频| 国产成人精品久久二区二区免费| 一级毛片女人18水好多| 中文字幕av在线有码专区| 午夜精品一区二区三区免费看| 国产精品电影一区二区三区| 午夜精品一区二区三区免费看| 欧美激情久久久久久爽电影| 在线a可以看的网站| 不卡av一区二区三区| 精品国内亚洲2022精品成人| 99热这里只有精品一区 | 国产精品免费一区二区三区在线| 日本精品一区二区三区蜜桃| 五月玫瑰六月丁香| 国产黄a三级三级三级人| 亚洲精品中文字幕一二三四区| 亚洲无线在线观看| 国产亚洲欧美在线一区二区| 午夜福利免费观看在线| 五月玫瑰六月丁香| 精品久久久久久久末码| 欧美日韩一级在线毛片| 精品熟女少妇八av免费久了| 国产97色在线日韩免费| 久久香蕉激情| 男人的好看免费观看在线视频 | 精品不卡国产一区二区三区| 老汉色∧v一级毛片| 小说图片视频综合网站| 又紧又爽又黄一区二区| 少妇熟女aⅴ在线视频| 中文资源天堂在线| 日日爽夜夜爽网站| 老熟妇乱子伦视频在线观看| 在线观看66精品国产| 两个人的视频大全免费| 99国产精品一区二区三区| 午夜激情av网站| 成人欧美大片| 国产亚洲精品av在线| 变态另类成人亚洲欧美熟女| 女同久久另类99精品国产91| 国产亚洲av高清不卡| 国产亚洲精品久久久久5区| 99国产综合亚洲精品| 色播亚洲综合网| 国产精华一区二区三区| 69av精品久久久久久| av超薄肉色丝袜交足视频| 夜夜爽天天搞| 国产欧美日韩一区二区精品| 久久久久国产一级毛片高清牌| 亚洲欧美日韩高清在线视频| 亚洲av中文字字幕乱码综合| 啦啦啦观看免费观看视频高清| 精品久久久久久久久久免费视频| 国产精品久久电影中文字幕| 搞女人的毛片| 久久香蕉精品热| 久久久久久九九精品二区国产 | 18美女黄网站色大片免费观看| 老鸭窝网址在线观看| tocl精华| 国产单亲对白刺激| 精品久久久久久久久久免费视频| 国产熟女xx| 国产精品亚洲美女久久久| 国产亚洲精品一区二区www| 在线十欧美十亚洲十日本专区| 白带黄色成豆腐渣| 此物有八面人人有两片| 他把我摸到了高潮在线观看| 国产精品日韩av在线免费观看| 精品久久久久久久毛片微露脸| 国产精品久久视频播放| 欧美丝袜亚洲另类 | 亚洲成人中文字幕在线播放| 久久人人精品亚洲av| 久久久国产欧美日韩av| 两个人免费观看高清视频| 听说在线观看完整版免费高清| 男女视频在线观看网站免费 | 国产精品乱码一区二三区的特点| 变态另类成人亚洲欧美熟女| www.自偷自拍.com| 国产精品av视频在线免费观看| 久久久国产精品麻豆| 天天添夜夜摸| 男人舔女人的私密视频| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩欧美在线乱码| 不卡一级毛片| 日韩欧美精品v在线| 18禁黄网站禁片免费观看直播| 亚洲精品色激情综合| 丝袜美腿诱惑在线| 免费人成视频x8x8入口观看| 久久久国产成人免费| 国产成人啪精品午夜网站| 丰满人妻熟妇乱又伦精品不卡| 国产精品久久电影中文字幕| 日本一本二区三区精品| 欧美日本视频| 老司机靠b影院| 国产精品香港三级国产av潘金莲| 国产精品九九99| 久久香蕉国产精品| 99久久无色码亚洲精品果冻| 精品免费久久久久久久清纯| 黄色片一级片一级黄色片| 久久伊人香网站| 亚洲免费av在线视频| aaaaa片日本免费| 成年女人毛片免费观看观看9| 成人亚洲精品av一区二区| 国产激情欧美一区二区| 国产成+人综合+亚洲专区| 欧美乱色亚洲激情| 脱女人内裤的视频| 免费在线观看视频国产中文字幕亚洲| 在线观看免费视频日本深夜| 午夜福利成人在线免费观看| 国产精品 国内视频| 欧美又色又爽又黄视频| 欧美日韩黄片免| 老司机在亚洲福利影院| 韩国av一区二区三区四区| 九色国产91popny在线| 99久久综合精品五月天人人| av福利片在线| 免费看十八禁软件| 免费人成视频x8x8入口观看| 国产一区二区激情短视频| 久久天堂一区二区三区四区| 少妇人妻一区二区三区视频| 国产高清视频在线播放一区| 亚洲专区国产一区二区| 757午夜福利合集在线观看| 一级毛片精品| 欧美久久黑人一区二区| 国产aⅴ精品一区二区三区波| 精品欧美一区二区三区在线| 中文字幕精品亚洲无线码一区| 亚洲五月婷婷丁香| 色精品久久人妻99蜜桃| 亚洲七黄色美女视频| 午夜亚洲福利在线播放| 嫩草影院精品99| 色综合站精品国产| 长腿黑丝高跟| 精品国产美女av久久久久小说| 日韩 欧美 亚洲 中文字幕| 欧美日本亚洲视频在线播放| 国产三级在线视频| 全区人妻精品视频| 亚洲精品一卡2卡三卡4卡5卡| 一本大道久久a久久精品| 国产av麻豆久久久久久久| xxx96com| √禁漫天堂资源中文www| 可以免费在线观看a视频的电影网站| 久久久久国产一级毛片高清牌| 老鸭窝网址在线观看| 香蕉丝袜av| 亚洲av第一区精品v没综合| 性色av乱码一区二区三区2| 91成年电影在线观看| 19禁男女啪啪无遮挡网站| 99在线人妻在线中文字幕| 99精品欧美一区二区三区四区| 美女免费视频网站| 色综合站精品国产| 国产成人影院久久av| 最近最新中文字幕大全免费视频| 黄频高清免费视频| 叶爱在线成人免费视频播放| 视频区欧美日本亚洲| 一个人免费在线观看的高清视频| 中文字幕人妻丝袜一区二区| 十八禁网站免费在线| 亚洲国产欧美人成| 久久久久久国产a免费观看| 91av网站免费观看| 免费看美女性在线毛片视频| 一个人免费在线观看的高清视频| 久久香蕉国产精品| 亚洲无线在线观看| 亚洲无线在线观看| a级毛片a级免费在线| 国产一区二区在线av高清观看| 99精品欧美一区二区三区四区| 国产精品,欧美在线| 午夜免费观看网址| 国产免费av片在线观看野外av| 2021天堂中文幕一二区在线观| 国产三级在线视频| 99国产极品粉嫩在线观看| 嫩草影视91久久| 成人国产一区最新在线观看| 欧美一区二区精品小视频在线| 伊人久久大香线蕉亚洲五| 此物有八面人人有两片| 国产三级在线视频| 国产精品久久久久久精品电影| 丰满的人妻完整版| 午夜免费观看网址| 国产私拍福利视频在线观看| 国产午夜精品久久久久久| 亚洲熟妇中文字幕五十中出| 欧美黑人欧美精品刺激| 男女床上黄色一级片免费看| 亚洲精品色激情综合| 天堂动漫精品| 久久中文字幕一级| 亚洲av第一区精品v没综合| 国产高清有码在线观看视频 | 亚洲七黄色美女视频| 亚洲精华国产精华精| 午夜免费激情av| 欧美日韩一级在线毛片| 亚洲人与动物交配视频| av片东京热男人的天堂| 成人三级做爰电影| 日韩欧美在线二视频| 级片在线观看| 男男h啪啪无遮挡| 人人妻人人看人人澡| 91九色精品人成在线观看| 国产精品乱码一区二三区的特点| 欧美精品啪啪一区二区三区| 激情在线观看视频在线高清| 在线观看66精品国产| 国产视频一区二区在线看| 国产99白浆流出| 大型黄色视频在线免费观看| 国产精品影院久久| www日本黄色视频网| 午夜福利欧美成人| 黄色女人牲交| 国产欧美日韩精品亚洲av| 国产精品日韩av在线免费观看| 91成年电影在线观看| 18禁美女被吸乳视频| 一级片免费观看大全| 亚洲五月婷婷丁香| 午夜精品久久久久久毛片777| 美女 人体艺术 gogo| 成年免费大片在线观看| 首页视频小说图片口味搜索| 蜜桃久久精品国产亚洲av| 精品一区二区三区av网在线观看| 国产99久久九九免费精品| netflix在线观看网站| 亚洲av五月六月丁香网| 久久香蕉精品热| 国产精品综合久久久久久久免费| 国产精品影院久久| avwww免费| 亚洲18禁久久av| 最近视频中文字幕2019在线8| 大型av网站在线播放| 国产av又大| 亚洲成a人片在线一区二区| 亚洲精品国产一区二区精华液| 久久人妻av系列| 一边摸一边做爽爽视频免费| 麻豆久久精品国产亚洲av| 男人舔奶头视频| 亚洲国产精品999在线| 久久亚洲精品不卡| 亚洲中文av在线| 欧美成人性av电影在线观看| 免费搜索国产男女视频| 99在线视频只有这里精品首页| 非洲黑人性xxxx精品又粗又长| 久久精品国产99精品国产亚洲性色| 美女扒开内裤让男人捅视频| 久久久精品大字幕| 亚洲国产日韩欧美精品在线观看 | 色综合亚洲欧美另类图片| 一个人免费在线观看的高清视频| 久久久久久久精品吃奶| 久久婷婷成人综合色麻豆| 波多野结衣高清无吗| 午夜a级毛片| 亚洲成人国产一区在线观看| 久久久久久大精品| 色播亚洲综合网| 国产成人aa在线观看| 狂野欧美激情性xxxx| 亚洲人成网站高清观看| 一本综合久久免费| 日韩中文字幕欧美一区二区| 亚洲乱码一区二区免费版| 日韩欧美三级三区| 久久中文字幕一级| 亚洲成av人片免费观看| 久久久久性生活片| 九色国产91popny在线| 白带黄色成豆腐渣| 日本免费一区二区三区高清不卡| 99re在线观看精品视频| 亚洲精品美女久久久久99蜜臀| 狠狠狠狠99中文字幕| 俺也久久电影网| 亚洲第一欧美日韩一区二区三区| 精品一区二区三区av网在线观看| 国产午夜精品论理片| 亚洲熟女毛片儿| 2021天堂中文幕一二区在线观| 成人av一区二区三区在线看| 欧美日韩精品网址| 亚洲欧美日韩东京热| 亚洲精品av麻豆狂野| 亚洲人成网站高清观看| 亚洲国产精品合色在线| 国产激情欧美一区二区| 亚洲av成人精品一区久久| 国产精品 欧美亚洲| 欧美3d第一页| av片东京热男人的天堂| 午夜精品久久久久久毛片777| 免费看十八禁软件| 久久精品国产亚洲av高清一级| 最近视频中文字幕2019在线8| 亚洲熟妇中文字幕五十中出| 天堂动漫精品| 国产精品 国内视频| 亚洲成a人片在线一区二区| 99热这里只有是精品50| 成人国产综合亚洲| 午夜福利在线在线| 少妇熟女aⅴ在线视频| 成人特级黄色片久久久久久久| 丁香欧美五月| 亚洲国产欧美人成| 一本久久中文字幕| 亚洲欧美精品综合久久99| 亚洲天堂国产精品一区在线| 日本熟妇午夜| 国产爱豆传媒在线观看 | 制服诱惑二区| 色av中文字幕| 亚洲成av人片在线播放无| 老汉色∧v一级毛片| 日本 av在线| 亚洲狠狠婷婷综合久久图片| 99riav亚洲国产免费| 免费看日本二区| 久久精品国产综合久久久| 亚洲精品在线美女| 啦啦啦免费观看视频1| 一个人免费在线观看的高清视频| 麻豆成人午夜福利视频| 中文在线观看免费www的网站 | 午夜免费观看网址| 两个人视频免费观看高清| 亚洲一区高清亚洲精品| 欧美日韩国产亚洲二区| 人妻夜夜爽99麻豆av| 午夜免费观看网址| 欧美黑人巨大hd| 亚洲第一欧美日韩一区二区三区| 欧美乱码精品一区二区三区| 好男人电影高清在线观看| 88av欧美| 日韩高清综合在线| 午夜成年电影在线免费观看| 在线十欧美十亚洲十日本专区| 国产av在哪里看| 亚洲av日韩精品久久久久久密| 禁无遮挡网站| 99久久综合精品五月天人人| 国产熟女xx| 18禁观看日本| 长腿黑丝高跟| 欧美日韩一级在线毛片| 亚洲午夜理论影院| 久久久久久人人人人人| 亚洲国产欧洲综合997久久,| a在线观看视频网站| 国内久久婷婷六月综合欲色啪| 欧美精品啪啪一区二区三区| 日韩 欧美 亚洲 中文字幕| 制服丝袜大香蕉在线| 制服人妻中文乱码| 国产亚洲精品综合一区在线观看 | 九色成人免费人妻av| 午夜精品在线福利| 级片在线观看| 变态另类成人亚洲欧美熟女| 成人精品一区二区免费| 三级国产精品欧美在线观看 | 国产一区二区三区视频了| 欧美黄色片欧美黄色片| 国产亚洲av嫩草精品影院| 一本精品99久久精品77| 老熟妇乱子伦视频在线观看| 国产在线观看jvid| 日本免费a在线| 97超级碰碰碰精品色视频在线观看| 久久精品国产综合久久久| 男人舔奶头视频| 91成年电影在线观看| 狂野欧美白嫩少妇大欣赏| 国产99久久九九免费精品| 91字幕亚洲| a级毛片a级免费在线| 美女 人体艺术 gogo| 好看av亚洲va欧美ⅴa在| 亚洲乱码一区二区免费版| 免费看美女性在线毛片视频| 国产区一区二久久| 99国产精品一区二区三区| 一级a爱片免费观看的视频| 狠狠狠狠99中文字幕| 亚洲欧美精品综合久久99| 亚洲一区中文字幕在线| 免费av毛片视频| 桃红色精品国产亚洲av| 中文字幕久久专区| 老熟妇仑乱视频hdxx| 长腿黑丝高跟| 免费在线观看亚洲国产| 国产成+人综合+亚洲专区| 无遮挡黄片免费观看| cao死你这个sao货| 久久草成人影院| 成人国产综合亚洲| 男女做爰动态图高潮gif福利片| 亚洲五月天丁香| 精品久久久久久久久久久久久| 成人特级黄色片久久久久久久| 亚洲熟妇熟女久久| 99热6这里只有精品| 窝窝影院91人妻| 国产精品久久久久久亚洲av鲁大| 精品久久久久久久毛片微露脸| 亚洲人与动物交配视频| 久久久国产欧美日韩av| 日日干狠狠操夜夜爽| 国产v大片淫在线免费观看| 日韩精品中文字幕看吧| 国产一区在线观看成人免费| 免费人成视频x8x8入口观看| 黑人欧美特级aaaaaa片| 我的老师免费观看完整版| 国产黄a三级三级三级人| 亚洲男人天堂网一区| 日韩大码丰满熟妇| 精品国产乱子伦一区二区三区| 大型黄色视频在线免费观看| 天天添夜夜摸| 成熟少妇高潮喷水视频| 国产精品免费视频内射| 99热只有精品国产| 午夜精品久久久久久毛片777| 搡老熟女国产l中国老女人| 2021天堂中文幕一二区在线观| 精品国产乱码久久久久久男人| 精品午夜福利视频在线观看一区| 国语自产精品视频在线第100页| 少妇被粗大的猛进出69影院| 俺也久久电影网| 国产亚洲精品av在线| 国产精品香港三级国产av潘金莲| 香蕉国产在线看| 成人精品一区二区免费| 日韩成人在线观看一区二区三区| 一个人免费在线观看电影 | 日本 欧美在线| 国产aⅴ精品一区二区三区波| 999久久久国产精品视频| 免费看美女性在线毛片视频| 国产成人啪精品午夜网站| 成人特级黄色片久久久久久久| 琪琪午夜伦伦电影理论片6080| 老司机午夜福利在线观看视频| 国产精品1区2区在线观看.| 日本精品一区二区三区蜜桃| 三级毛片av免费| 亚洲欧美精品综合一区二区三区| 国产欧美日韩精品亚洲av| 精品无人区乱码1区二区| 757午夜福利合集在线观看| 丰满的人妻完整版| 日本一区二区免费在线视频| 757午夜福利合集在线观看| 1024香蕉在线观看| 美女 人体艺术 gogo| 欧美日韩瑟瑟在线播放| 久久久久久久久久黄片| 色av中文字幕| 在线视频色国产色| 久久这里只有精品19| 国产三级在线视频| 国产单亲对白刺激| av中文乱码字幕在线| 亚洲人成伊人成综合网2020| 久久香蕉精品热| 一边摸一边抽搐一进一小说| 欧美日本视频| 国产午夜精品久久久久久| 成人一区二区视频在线观看| 亚洲精品美女久久av网站| 色综合欧美亚洲国产小说| 亚洲电影在线观看av| 淫秽高清视频在线观看| svipshipincom国产片| ponron亚洲| 一边摸一边抽搐一进一小说| 一区二区三区国产精品乱码| 啦啦啦免费观看视频1| 国产av又大| 国产成人影院久久av| 色综合婷婷激情| 国产一区二区在线观看日韩 | 啪啪无遮挡十八禁网站| 久久久精品大字幕| 日韩欧美精品v在线| 亚洲五月婷婷丁香| 国产精品亚洲美女久久久| 亚洲,欧美精品.| 久久久国产成人精品二区| 精品国产超薄肉色丝袜足j| 丝袜美腿诱惑在线| 91老司机精品| av国产免费在线观看| 国产熟女午夜一区二区三区| 成年版毛片免费区| 国产亚洲欧美在线一区二区| 一卡2卡三卡四卡精品乱码亚洲| 不卡一级毛片| 免费搜索国产男女视频| 波多野结衣巨乳人妻| 丰满的人妻完整版| 两人在一起打扑克的视频| 老汉色∧v一级毛片| 亚洲国产精品成人综合色| 国产激情欧美一区二区| 亚洲精品久久成人aⅴ小说| 三级毛片av免费| 一区二区三区高清视频在线| av视频在线观看入口| 中文字幕av在线有码专区| 少妇被粗大的猛进出69影院| 久久这里只有精品19| 中文亚洲av片在线观看爽| 香蕉丝袜av| 可以在线观看毛片的网站| 五月伊人婷婷丁香| 国产真实乱freesex| 国产精品98久久久久久宅男小说| 亚洲av电影在线进入| 欧美黄色淫秽网站| 亚洲精品粉嫩美女一区| 男人的好看免费观看在线视频 | 又粗又爽又猛毛片免费看| 日韩大码丰满熟妇| 看片在线看免费视频| 日韩精品青青久久久久久| 亚洲精品国产精品久久久不卡| 亚洲欧美日韩高清专用| 久久久精品大字幕| 久久久久久九九精品二区国产 | 给我免费播放毛片高清在线观看| 欧美国产日韩亚洲一区| 一边摸一边做爽爽视频免费| 在线免费观看的www视频| 一本久久中文字幕| 九色国产91popny在线| 动漫黄色视频在线观看| 十八禁网站免费在线| 国产熟女午夜一区二区三区| 一级毛片精品| 国产精品久久久久久亚洲av鲁大| 成在线人永久免费视频| 久久久国产成人精品二区| 久久久久国产精品人妻aⅴ院| 美女午夜性视频免费| 精品久久久久久久末码| netflix在线观看网站| 中文字幕人妻丝袜一区二区| 99热这里只有精品一区 | 桃色一区二区三区在线观看| 成人av一区二区三区在线看| 一边摸一边抽搐一进一小说| 在线十欧美十亚洲十日本专区| 一个人免费在线观看的高清视频| 亚洲av成人精品一区久久| 男女做爰动态图高潮gif福利片| 首页视频小说图片口味搜索| ponron亚洲| 久久香蕉激情| 国产一区二区三区视频了| 久久天躁狠狠躁夜夜2o2o| 国产免费男女视频| 床上黄色一级片| 最近最新中文字幕大全电影3| 午夜影院日韩av| 久久久精品国产亚洲av高清涩受| 国产精品综合久久久久久久免费| 午夜精品在线福利| 一级黄色大片毛片| 欧美zozozo另类| 老司机深夜福利视频在线观看| 视频区欧美日本亚洲| 日日爽夜夜爽网站| 欧美黑人欧美精品刺激| 级片在线观看| 久久99热这里只有精品18| 久久久久国产一级毛片高清牌|