• <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ī)攻角傳感器安裝定位研究
    久热爱精品视频在线9| 亚洲,一卡二卡三卡| 久久精品久久精品一区二区三区| 国产精品99久久99久久久不卡 | 久久人妻熟女aⅴ| 一区二区三区精品91| 十分钟在线观看高清视频www| 日本爱情动作片www.在线观看| 在现免费观看毛片| 久久久久久人人人人人| 国产成人系列免费观看| 人体艺术视频欧美日本| 看非洲黑人一级黄片| 青春草亚洲视频在线观看| 精品久久久精品久久久| 精品人妻一区二区三区麻豆| 免费高清在线观看日韩| 七月丁香在线播放| 大片免费播放器 马上看| 久久久精品区二区三区| 国产一卡二卡三卡精品 | 久久午夜综合久久蜜桃| 九色亚洲精品在线播放| 国产精品香港三级国产av潘金莲 | 欧美日韩一级在线毛片| 国精品久久久久久国模美| 亚洲欧美日韩另类电影网站| 欧美激情 高清一区二区三区| 最近中文字幕2019免费版| 男人添女人高潮全过程视频| 国产福利在线免费观看视频| 在线观看三级黄色| 又黄又粗又硬又大视频| 少妇精品久久久久久久| 看十八女毛片水多多多| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品久久午夜乱码| 在线观看www视频免费| 少妇精品久久久久久久| 国产成人精品在线电影| 如日韩欧美国产精品一区二区三区| 人体艺术视频欧美日本| 免费高清在线观看视频在线观看| 丰满迷人的少妇在线观看| 一级毛片我不卡| 九九爱精品视频在线观看| 99国产精品免费福利视频| 亚洲国产欧美一区二区综合| 国产亚洲午夜精品一区二区久久| 成年人午夜在线观看视频| 久久久国产欧美日韩av| 欧美日韩综合久久久久久| 如日韩欧美国产精品一区二区三区| 最新在线观看一区二区三区 | 在线观看人妻少妇| 1024香蕉在线观看| 日韩av不卡免费在线播放| 777米奇影视久久| 毛片一级片免费看久久久久| 中文字幕制服av| 亚洲久久久国产精品| 欧美亚洲 丝袜 人妻 在线| 国产成人啪精品午夜网站| 免费在线观看视频国产中文字幕亚洲 | 日韩人妻精品一区2区三区| 亚洲国产欧美在线一区| 国产精品久久久久久精品电影小说| 一级爰片在线观看| 欧美日韩亚洲综合一区二区三区_| 日本wwww免费看| 中文字幕高清在线视频| 亚洲熟女毛片儿| 赤兔流量卡办理| 最近最新中文字幕大全免费视频 | 王馨瑶露胸无遮挡在线观看| 女人被躁到高潮嗷嗷叫费观| 日本黄色日本黄色录像| 少妇人妻 视频| 国产免费视频播放在线视频| 国产欧美日韩一区二区三区在线| 天天添夜夜摸| 国产极品粉嫩免费观看在线| 中文字幕高清在线视频| 精品亚洲成国产av| 欧美最新免费一区二区三区| 各种免费的搞黄视频| 日韩欧美精品免费久久| 在线观看一区二区三区激情| 男的添女的下面高潮视频| 熟女少妇亚洲综合色aaa.| 女性被躁到高潮视频| 搡老岳熟女国产| 中文字幕另类日韩欧美亚洲嫩草| 国产成人啪精品午夜网站| 国产精品一区二区精品视频观看| 黄片小视频在线播放| 欧美日韩成人在线一区二区| 黄片播放在线免费| 亚洲专区中文字幕在线 | 老司机亚洲免费影院| 国产一区有黄有色的免费视频| 一级毛片黄色毛片免费观看视频| 欧美日韩亚洲国产一区二区在线观看 | av一本久久久久| 一级毛片黄色毛片免费观看视频| 性高湖久久久久久久久免费观看| av网站免费在线观看视频| 亚洲精品视频女| 人人妻人人澡人人看| 国产男人的电影天堂91| 制服人妻中文乱码| 99久久综合免费| 下体分泌物呈黄色| 激情视频va一区二区三区| av在线app专区| 国产色婷婷99| av有码第一页| 国产精品久久久久久人妻精品电影 | 岛国毛片在线播放| av.在线天堂| 如日韩欧美国产精品一区二区三区| 日韩不卡一区二区三区视频在线| 中文字幕最新亚洲高清| 亚洲国产最新在线播放| 久久精品aⅴ一区二区三区四区| 欧美精品人与动牲交sv欧美| 国产高清国产精品国产三级| 18禁观看日本| 人体艺术视频欧美日本| 国产亚洲最大av| 久久久久久免费高清国产稀缺| 日韩,欧美,国产一区二区三区| 操出白浆在线播放| www日本在线高清视频| 精品国产一区二区三区四区第35| 秋霞伦理黄片| 欧美日韩视频精品一区| 欧美日韩国产mv在线观看视频| 精品亚洲成a人片在线观看| 黄频高清免费视频| 久久影院123| 久久精品久久久久久久性| 不卡视频在线观看欧美| 亚洲免费av在线视频| 在现免费观看毛片| 中文字幕人妻丝袜一区二区 | 国产成人系列免费观看| av在线观看视频网站免费| 国产成人av激情在线播放| 国产成人欧美| videosex国产| 国产男女内射视频| 国产一级毛片在线| 国产av码专区亚洲av| 精品久久蜜臀av无| 超碰97精品在线观看| 看免费成人av毛片| av又黄又爽大尺度在线免费看| 欧美精品亚洲一区二区| 老汉色av国产亚洲站长工具| 亚洲精品乱久久久久久| 日韩一区二区三区影片| 国产成人精品在线电影| 亚洲欧洲精品一区二区精品久久久 | 考比视频在线观看| 99精品久久久久人妻精品| 亚洲国产最新在线播放| 另类亚洲欧美激情| av福利片在线| 欧美精品亚洲一区二区| 久久女婷五月综合色啪小说| 午夜激情久久久久久久| 国产成人欧美| 中文天堂在线官网| 亚洲 欧美一区二区三区| 午夜精品国产一区二区电影| 天天操日日干夜夜撸| 久久久久人妻精品一区果冻| 伊人久久国产一区二区| 久久免费观看电影| 男男h啪啪无遮挡| 国产精品久久久av美女十八| 国产黄频视频在线观看| 日韩免费高清中文字幕av| h视频一区二区三区| 欧美日韩成人在线一区二区| 91精品国产国语对白视频| 十八禁高潮呻吟视频| 满18在线观看网站| 男男h啪啪无遮挡| 免费黄频网站在线观看国产| 久久精品亚洲av国产电影网| 国产精品久久久久成人av| 亚洲国产看品久久| 亚洲欧美成人综合另类久久久| 中文字幕精品免费在线观看视频| 亚洲精品美女久久av网站| 亚洲成国产人片在线观看| 侵犯人妻中文字幕一二三四区| 捣出白浆h1v1| 欧美少妇被猛烈插入视频| 久久久久久久大尺度免费视频| 成人亚洲欧美一区二区av| 可以免费在线观看a视频的电影网站 | 午夜福利在线免费观看网站| 亚洲成国产人片在线观看| 在线免费观看不下载黄p国产| 宅男免费午夜| 国产成人欧美在线观看 | 久久精品aⅴ一区二区三区四区| 1024视频免费在线观看| 日本av手机在线免费观看| 国产亚洲av片在线观看秒播厂| 一级,二级,三级黄色视频| 高清不卡的av网站| av视频免费观看在线观看| 少妇被粗大猛烈的视频| 久热爱精品视频在线9| 亚洲中文av在线| 777米奇影视久久| 黄网站色视频无遮挡免费观看| 黄色一级大片看看| 国产亚洲午夜精品一区二区久久| 涩涩av久久男人的天堂| 亚洲一卡2卡3卡4卡5卡精品中文| 我要看黄色一级片免费的| 国产成人一区二区在线| 久久久久视频综合| 国产精品久久久久久久久免| 一级片免费观看大全| 国产99久久九九免费精品| 亚洲国产毛片av蜜桃av| 在线观看免费午夜福利视频| 丝袜脚勾引网站| 人妻一区二区av| 日韩一卡2卡3卡4卡2021年| 亚洲欧美一区二区三区国产| 免费观看人在逋| 久久综合国产亚洲精品| 日韩欧美一区视频在线观看| 亚洲国产av新网站| 国产亚洲一区二区精品| 久久人妻熟女aⅴ| 久久99热这里只频精品6学生| 国产又爽黄色视频| 国产一区二区在线观看av| 久久人人爽人人片av| 国产一区二区激情短视频 | 亚洲伊人久久精品综合| 老司机靠b影院| 中文天堂在线官网| 欧美 亚洲 国产 日韩一| 高清视频免费观看一区二区| 多毛熟女@视频| 国产在线一区二区三区精| 最新的欧美精品一区二区| 人妻人人澡人人爽人人| 深夜精品福利| 免费av中文字幕在线| 韩国精品一区二区三区| 久久精品亚洲av国产电影网| 一二三四在线观看免费中文在| 欧美日韩综合久久久久久| 黄片播放在线免费| 人妻 亚洲 视频| 精品国产乱码久久久久久小说| 欧美日韩av久久| 成人免费观看视频高清| 国产精品一区二区在线不卡| 国产亚洲精品第一综合不卡| 天天躁夜夜躁狠狠久久av| 国产精品女同一区二区软件| 成人午夜精彩视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 麻豆乱淫一区二区| 久久天躁狠狠躁夜夜2o2o | 99香蕉大伊视频| 欧美日韩视频高清一区二区三区二| 国产色婷婷99| 亚洲,欧美,日韩| 国产一区二区在线观看av| 黄网站色视频无遮挡免费观看| 亚洲欧美精品自产自拍| 亚洲国产最新在线播放| 看免费成人av毛片| 91老司机精品| 国产一卡二卡三卡精品 | 自线自在国产av| 国产精品偷伦视频观看了| 美女国产高潮福利片在线看| 免费观看人在逋| 黄色毛片三级朝国网站| 精品少妇黑人巨大在线播放| 亚洲精品国产色婷婷电影| 久久久国产欧美日韩av| 精品一区二区免费观看| 天天躁夜夜躁狠狠躁躁| 一区在线观看完整版| 免费黄色在线免费观看| 亚洲av电影在线进入| 亚洲欧美成人综合另类久久久| 亚洲欧美日韩另类电影网站| 男女边摸边吃奶| 国产精品嫩草影院av在线观看| 看十八女毛片水多多多| 久久国产亚洲av麻豆专区| 亚洲人成网站在线观看播放| 黑人猛操日本美女一级片| 搡老岳熟女国产| 欧美日韩一级在线毛片| 女人被躁到高潮嗷嗷叫费观| 这个男人来自地球电影免费观看 | 青草久久国产| 美女国产高潮福利片在线看| 国产成人精品久久二区二区91 | 亚洲欧美成人精品一区二区| 欧美少妇被猛烈插入视频| 国产又爽黄色视频| 18禁裸乳无遮挡动漫免费视频| 九九爱精品视频在线观看| 大香蕉久久成人网| 国产日韩欧美视频二区| a级毛片在线看网站| 五月天丁香电影| 久久女婷五月综合色啪小说| 美女福利国产在线| 国产一区二区激情短视频 | √禁漫天堂资源中文www| 丁香六月天网| 日韩大片免费观看网站| 国产野战对白在线观看| 亚洲伊人久久精品综合| 国产日韩欧美亚洲二区| 日本vs欧美在线观看视频| 男人舔女人的私密视频| 日韩伦理黄色片| 超碰97精品在线观看| 最新的欧美精品一区二区| 9热在线视频观看99| 丝袜美腿诱惑在线| 亚洲婷婷狠狠爱综合网| 亚洲成人av在线免费| 久久青草综合色| 久久精品aⅴ一区二区三区四区| xxx大片免费视频| 女的被弄到高潮叫床怎么办| 国产精品国产三级专区第一集| 欧美日韩av久久| 亚洲七黄色美女视频| 日韩免费高清中文字幕av| 亚洲av电影在线进入| 一区二区三区四区激情视频| 久久鲁丝午夜福利片| 久久久久久久久免费视频了| 精品少妇黑人巨大在线播放| 精品少妇久久久久久888优播| 欧美激情 高清一区二区三区| √禁漫天堂资源中文www| 国产淫语在线视频| 国产男人的电影天堂91| h视频一区二区三区| 亚洲国产欧美在线一区| 精品一品国产午夜福利视频| 中文精品一卡2卡3卡4更新| 亚洲少妇的诱惑av| 七月丁香在线播放| 国产日韩欧美在线精品| 国产男女内射视频| 两个人免费观看高清视频| 又粗又硬又长又爽又黄的视频| 新久久久久国产一级毛片| 少妇人妻久久综合中文| 久久精品国产a三级三级三级| 亚洲欧洲国产日韩| 熟妇人妻不卡中文字幕| 性少妇av在线| 女人高潮潮喷娇喘18禁视频| 少妇人妻久久综合中文| 精品亚洲成国产av| 国语对白做爰xxxⅹ性视频网站| 丁香六月天网| 99精国产麻豆久久婷婷| 国产成人系列免费观看| 国产国语露脸激情在线看| 亚洲三区欧美一区| 国产免费现黄频在线看| 搡老乐熟女国产| 久久婷婷青草| 免费不卡黄色视频| 韩国高清视频一区二区三区| 色94色欧美一区二区| 亚洲精品国产色婷婷电影| 伊人久久国产一区二区| 国产精品偷伦视频观看了| 国产精品熟女久久久久浪| 高清视频免费观看一区二区| xxx大片免费视频| 另类精品久久| 肉色欧美久久久久久久蜜桃| 欧美精品av麻豆av| 中文字幕亚洲精品专区| 亚洲激情五月婷婷啪啪| 一二三四中文在线观看免费高清| 两个人免费观看高清视频| 亚洲av成人精品一二三区| 水蜜桃什么品种好| 老司机影院毛片| 青春草亚洲视频在线观看| 天天影视国产精品| 亚洲一卡2卡3卡4卡5卡精品中文| 9热在线视频观看99| 久久精品人人爽人人爽视色| av一本久久久久| 国产精品久久久av美女十八| 国产熟女午夜一区二区三区| 水蜜桃什么品种好| 肉色欧美久久久久久久蜜桃| 久久 成人 亚洲| 久久久欧美国产精品| 日本一区二区免费在线视频| 极品少妇高潮喷水抽搐| 国产精品亚洲av一区麻豆 | 十分钟在线观看高清视频www| 91成人精品电影| 国产精品亚洲av一区麻豆 | 日韩成人av中文字幕在线观看| 国产精品香港三级国产av潘金莲 | 欧美亚洲日本最大视频资源| 久久av网站| 亚洲国产av新网站| 欧美亚洲日本最大视频资源| 男女下面插进去视频免费观看| 亚洲色图综合在线观看| 亚洲第一av免费看| 街头女战士在线观看网站| av女优亚洲男人天堂| 1024香蕉在线观看| av片东京热男人的天堂| 麻豆精品久久久久久蜜桃| 宅男免费午夜| 精品久久蜜臀av无| 人人妻,人人澡人人爽秒播 | 自线自在国产av| 色精品久久人妻99蜜桃| 日韩制服骚丝袜av| 韩国av在线不卡| 亚洲欧美成人精品一区二区| 色婷婷av一区二区三区视频| 丝袜在线中文字幕| 久久99精品国语久久久| 国产一级毛片在线| 日韩制服骚丝袜av| 国产xxxxx性猛交| 免费av中文字幕在线| 亚洲成人国产一区在线观看 | 飞空精品影院首页| 久久久精品国产亚洲av高清涩受| 精品亚洲成国产av| 国产亚洲午夜精品一区二区久久| 午夜福利,免费看| 久久久久久久精品精品| 大香蕉久久网| 亚洲国产最新在线播放| 欧美精品一区二区大全| 国产精品国产三级国产专区5o| 午夜av观看不卡| 中文字幕亚洲精品专区| 免费不卡黄色视频| 美女午夜性视频免费| 丰满乱子伦码专区| 亚洲精品美女久久久久99蜜臀 | 久久久久久人妻| 日本wwww免费看| 欧美av亚洲av综合av国产av | 午夜福利一区二区在线看| 国产一区亚洲一区在线观看| 精品人妻一区二区三区麻豆| 久久久精品国产亚洲av高清涩受| 丝袜在线中文字幕| 久久99精品国语久久久| 日韩 亚洲 欧美在线| 国产男女超爽视频在线观看| 最近中文字幕高清免费大全6| 80岁老熟妇乱子伦牲交| 狂野欧美激情性xxxx| 中文字幕亚洲精品专区| 亚洲熟女精品中文字幕| 日韩 亚洲 欧美在线| 香蕉丝袜av| 亚洲精品国产av成人精品| 90打野战视频偷拍视频| 综合色丁香网| kizo精华| 777米奇影视久久| 日本91视频免费播放| 国产精品秋霞免费鲁丝片| 国产无遮挡羞羞视频在线观看| 色精品久久人妻99蜜桃| 亚洲欧美精品综合一区二区三区| 免费观看人在逋| 久久精品国产综合久久久| 人人妻人人爽人人添夜夜欢视频| 男女边吃奶边做爰视频| 日日啪夜夜爽| 日韩伦理黄色片| 免费看av在线观看网站| 久久性视频一级片| 天堂俺去俺来也www色官网| 十八禁人妻一区二区| 午夜日韩欧美国产| 欧美黑人欧美精品刺激| 久久午夜综合久久蜜桃| 精品人妻在线不人妻| 99热网站在线观看| 久久久久久久国产电影| 免费高清在线观看视频在线观看| 美女国产高潮福利片在线看| 精品国产一区二区三区久久久樱花| 亚洲精品av麻豆狂野| 欧美黑人欧美精品刺激| xxxhd国产人妻xxx| 亚洲成av片中文字幕在线观看| www.熟女人妻精品国产| 成年动漫av网址| 精品一区在线观看国产| 久久久精品免费免费高清| √禁漫天堂资源中文www| 精品一品国产午夜福利视频| 日日撸夜夜添| 久久久国产欧美日韩av| 日本vs欧美在线观看视频| 亚洲少妇的诱惑av| 亚洲国产av影院在线观看| 伊人久久国产一区二区| 亚洲欧洲国产日韩| 人人妻,人人澡人人爽秒播 | 午夜福利免费观看在线| 美女福利国产在线| 久久久久久久久久久免费av| 久久精品国产综合久久久| 成人亚洲精品一区在线观看| 一级毛片我不卡| 午夜福利免费观看在线| 咕卡用的链子| 国产成人精品在线电影| 国产在视频线精品| 精品免费久久久久久久清纯 | 免费久久久久久久精品成人欧美视频| 久久精品aⅴ一区二区三区四区| 免费av中文字幕在线| 久久亚洲国产成人精品v| 观看美女的网站| av有码第一页| 9191精品国产免费久久| 你懂的网址亚洲精品在线观看| av电影中文网址| 久久99一区二区三区| 女性被躁到高潮视频| 麻豆av在线久日| 成人影院久久| 午夜精品国产一区二区电影| 久久国产精品大桥未久av| 久久热在线av| 午夜影院在线不卡| 国产极品粉嫩免费观看在线| 国产日韩欧美视频二区| 中文字幕人妻丝袜一区二区 | 一边摸一边抽搐一进一出视频| kizo精华| 伦理电影免费视频| 巨乳人妻的诱惑在线观看| 久久性视频一级片| 久久精品亚洲熟妇少妇任你| 欧美日韩亚洲综合一区二区三区_| 国产成人免费无遮挡视频| 青春草视频在线免费观看| 国产精品 国内视频| 在现免费观看毛片| 看免费成人av毛片| www.熟女人妻精品国产| 啦啦啦啦在线视频资源| 老司机深夜福利视频在线观看 | videosex国产| 欧美精品av麻豆av| 久久狼人影院| 国产成人免费观看mmmm| 欧美亚洲日本最大视频资源| 国产精品 国内视频| 久久精品国产综合久久久| 免费在线观看黄色视频的| 免费黄色在线免费观看| 91aial.com中文字幕在线观看| 亚洲成人一二三区av| av网站免费在线观看视频| 大片电影免费在线观看免费| 1024香蕉在线观看| 成人漫画全彩无遮挡| 19禁男女啪啪无遮挡网站| 亚洲欧美清纯卡通| 日韩av不卡免费在线播放| 女的被弄到高潮叫床怎么办| 亚洲欧洲国产日韩| 久久婷婷青草| 黄色毛片三级朝国网站| av.在线天堂| 一边摸一边做爽爽视频免费| 亚洲天堂av无毛| 国产深夜福利视频在线观看| 亚洲国产欧美在线一区| 国产毛片在线视频| 天天躁日日躁夜夜躁夜夜| 久久天躁狠狠躁夜夜2o2o | 欧美成人精品欧美一级黄| 久久久久视频综合| kizo精华|