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

    燃料電池氣體循環(huán)泵多葉轉(zhuǎn)子腔內(nèi)部氣動性能數(shù)值研究

    2021-03-08 05:52:00黎義斌李龍劉建峰楊悅民馮德瑋
    西安交通大學學報 2021年3期
    關(guān)鍵詞:葉數(shù)循環(huán)泵子葉

    黎義斌,李龍,劉建峰,楊悅民,馮德瑋

    (1.蘭州理工大學能源與動力工程學院,730050,蘭州;2.北京衛(wèi)星制造廠有限公司,100094,北京)

    凸輪式氣體循環(huán)泵是一種非接觸容積式泵,又常被稱為羅茨泵,通過泵內(nèi)相互嚙合的兩個轉(zhuǎn)子進行同步反向轉(zhuǎn)動,將空氣、氫氣、氧氣、氮氣、甲烷等氣體介質(zhì)由泵入口處吸入并充滿腔室,并將氣體排出。泵在工作過程中,兩轉(zhuǎn)子之間及轉(zhuǎn)子與腔壁之間存在一定間隙,避免產(chǎn)生磨損而造成損壞,延長泵的使用壽命,同時可以使泵達到較高的轉(zhuǎn)速。凸輪式氣體循環(huán)泵在輸送氣體時具有較高的真空度及耐磨、耐腐蝕性,并且可靠性高、結(jié)構(gòu)簡單、調(diào)節(jié)范圍廣,因此在航空航天、新能源汽車、船舶等領(lǐng)域具有廣泛的應用前景。

    近年來,國內(nèi)外學者對凸輪式氣體循環(huán)泵的研究不斷深入,其中主要集中于對凸輪轉(zhuǎn)子型線的設計及優(yōu)化、內(nèi)部流場的分析等。文獻[1-4]對傳統(tǒng)漸開線轉(zhuǎn)子型線進行改進,提高了轉(zhuǎn)子傳動的平穩(wěn)性和轉(zhuǎn)子傳動的承載能力。Wang等提出了一種由橢圓弧及其共軛曲線組成的新型橢圓羅茨泵轉(zhuǎn)子輪廓[5]。劉振超等將NURYS曲線運用到羅茨泵轉(zhuǎn)子型線的設計中,使轉(zhuǎn)子型線的面積利用系數(shù)提高了1%[6]。Li等提出一種新型漸變間隙的轉(zhuǎn)子腔結(jié)構(gòu),并分析了其對壓力脈動和徑向激勵力的影響[7]。Hsieh提出一種可變擺線比的擺線型線,同時基于橢圓軸比參數(shù)建立了新型橢圓形曲線,并驗證了兩種曲線的優(yōu)越性[8-9]。Yao針對羅茨鼓風機提出了一種3葉螺旋轉(zhuǎn)子,該轉(zhuǎn)子能夠產(chǎn)生更多氣流并降低運行噪聲[10]。

    隨著計算機仿真技術(shù)的不斷發(fā)展,越來越多的研究人員利用計算流體力學(CFD)對泵內(nèi)部的流動情況進行仿真。Li等通過瞬態(tài)數(shù)值模擬研究了黏度、間隙結(jié)構(gòu)和速度對流速和壓力的影響,并通過實驗進行了驗證,同時闡明了螺旋角、徑長比等幾何參數(shù)對凸輪泵轉(zhuǎn)子腔內(nèi)部流量特性的影響規(guī)律[11-13]。盧陽等對等螺距和變螺距2種螺桿真空泵分別進行了數(shù)值分析和實驗,得出變螺距真空泵性能優(yōu)于等螺距真空泵[14]。Guo等引入Zwart-Gerber-Belamri空化模型,對凸輪轉(zhuǎn)子泵內(nèi)部空化流進行數(shù)值分析[15]。翟云飛等詳細討論了3葉羅茨泵內(nèi)吸氣、輸運、反沖、排氣4個階段的特點和機理[16]。劉瑞青等對圓弧型2葉、3葉直葉和扭葉羅茨真空泵這3種轉(zhuǎn)子進行數(shù)值研究其性能[17]。Sun等使用動態(tài)網(wǎng)格方法來模擬羅茨泵內(nèi)流動,研制羅茨泵測試壓力分布和質(zhì)量流率,從而驗證CFD模型[18]。Guo等研究了以CO2和水為工作流體的氣液混合羅茨泵[19]。

    轉(zhuǎn)子葉數(shù)、轉(zhuǎn)子間隙等對凸輪式氣體循環(huán)泵性能影響較大,為了研究轉(zhuǎn)子葉數(shù)、轉(zhuǎn)子間隙對氣體循環(huán)泵綜合性能影響,本文通過建立圓弧-漸開線-圓弧型轉(zhuǎn)子型線方程,以3~6葉4種不同轉(zhuǎn)子葉數(shù)的凸輪式氣體循環(huán)泵為研究對象,采用動網(wǎng)格技術(shù)和RNGk-ε湍流模型,對不同葉數(shù)的轉(zhuǎn)子進行三維非定常數(shù)值分析,揭示葉數(shù)、轉(zhuǎn)子間隙對凸輪式氣體循環(huán)泵性能參數(shù)的影響規(guī)律,提高對其內(nèi)部流動特性的認識。

    1 數(shù)學模型

    1.1 轉(zhuǎn)子型線設計

    凸輪式氣體循環(huán)泵多葉轉(zhuǎn)子采用圓弧-漸開線-圓弧型線,以5葉轉(zhuǎn)子為例,如圖1所示。

    IG—葉根圓弧段;GK—漸開線段;KJ—葉峰圓弧段;R—漸開線GK的基圓半徑。圖1 圓弧-漸開線-圓弧型多葉轉(zhuǎn)子型線

    由于型線是對稱的,只需推導出∠BO1P內(nèi)的型線方程,即可繪制轉(zhuǎn)子型線。改變轉(zhuǎn)子葉數(shù)Z、中心距A、葉頂圓半徑Rm,通過推導即可得到不同葉數(shù)的轉(zhuǎn)子型線方程。

    (1)分別以O1、O2為圓心,轉(zhuǎn)子中心距的一半A/2為節(jié)圓半徑,繪制兩轉(zhuǎn)子的節(jié)圓。兩節(jié)圓相切于P點,過P點作兩節(jié)圓的公切線Y,再過P點做直線N1N2。N1N2分別與兩節(jié)圓相切,其與Y軸的夾角為節(jié)圓壓力角α,且

    (1)

    以R為半徑畫兩轉(zhuǎn)子的基圓,且

    (2)

    以O1為圓心畫一條射線,與節(jié)圓交于B點,且該射線與y軸正向的夾角為型線剩余角γ,γ與轉(zhuǎn)子葉數(shù)有關(guān),再過B點作基圓的切線,切點為M,且BM切線與漸開線段的交點為G,其中

    (3)

    (4)

    則GK漸開線段型線方程為

    (5)

    其中GK漸開線段展開起始角φmin和終止角φmax為

    (6)

    (2)以B為圓心、BG為半徑,繪制葉根圓弧段IG。設B點坐標為(XB,YB),有

    (7)

    則IG葉根圓弧段型線方程為

    (8)

    其中IG圓弧段展開起始角φmin和終止角φmax為

    (9)

    (3)以P點為圓心、PK為半徑繪制葉峰圓弧段KJ,設P點坐標為(XP,YP),XP=0,YP=A/2,則KJ葉峰圓弧段型線方程為

    (10)

    其中KJ圓弧段展開起始角φmin和終止角φmax為

    (11)

    1.2 幾何模型

    以圓弧-漸開線-圓弧型凸輪式氣體循環(huán)泵為研究對象,建立3葉、4葉、5葉、6葉轉(zhuǎn)子的氣體循環(huán)泵計算模型,并分別選取0.1、0.15、0.2、0.25、0.3 mm 5種不同的間隙,分析不同轉(zhuǎn)子葉數(shù)、不同間隙對轉(zhuǎn)子腔氣動性能的影響。轉(zhuǎn)子的額定參數(shù)為Rm=45 mm、轉(zhuǎn)子長度L=60 mm、進出口直徑d=40 mm、Rm/L=0.75、轉(zhuǎn)速n=3 000 r·min-1,轉(zhuǎn)子的其他設計參數(shù)如表1所示,多葉轉(zhuǎn)子三維模型如圖2所示。

    表1 多葉轉(zhuǎn)子設計參數(shù)

    (a)3葉(b)4葉

    (c)5葉(d)6葉圖2 多葉轉(zhuǎn)子三維模型

    由于轉(zhuǎn)子與轉(zhuǎn)子之間、轉(zhuǎn)子與泵殼之間、轉(zhuǎn)子端面與泵蓋之間均存在微小間隙,在泵運轉(zhuǎn)過程中都會造成氣體泄漏,影響泵的性能。轉(zhuǎn)壁間隙兩端分別與低壓腔/高壓腔和過渡腔相連接,或兩端均與過渡腔連接,而轉(zhuǎn)子間間隙兩端分別與高壓腔和低壓腔直接相連,使高壓腔與低壓腔之間形成一個直連泄漏通道,并且隨著轉(zhuǎn)子轉(zhuǎn)動泄漏通道形狀不斷變化。相比轉(zhuǎn)壁間隙,轉(zhuǎn)子間隙兩端壓差更大,更容易造成泄漏,因此轉(zhuǎn)子間隙形成的泄漏是影響凸輪式氣體循環(huán)泵性能的主要原因。本文忽略端面間隙產(chǎn)生的泄漏,只對轉(zhuǎn)子間間隙進行模擬分析。

    2 數(shù)值模擬

    2.1 網(wǎng)格劃分及邊界條件設置

    以5葉轉(zhuǎn)子為例,如圖3所示,最大網(wǎng)格單元設置為0.005,最小網(wǎng)格單元設置為0.000 2,網(wǎng)格區(qū)域分為入口、出口和旋轉(zhuǎn)流場部分,對兩個相鄰的不同域之間建立interface交互面。邊界條件設置為入口壓力202.65 kPa、出口壓力243.05 kPa、壓升40.4 kPa,輸送介質(zhì)為空氣。為了盡量減小網(wǎng)格數(shù)對數(shù)值計算結(jié)果的影響,同時參考文獻[20]中網(wǎng)格驗證方法,以5葉轉(zhuǎn)子為例,通過改變轉(zhuǎn)子的網(wǎng)格劃分尺寸,分別生成了7種網(wǎng)格數(shù)的尺寸。

    圖3 氣體循環(huán)泵計算模型的網(wǎng)格劃分

    如圖4所示,當網(wǎng)格數(shù)由41.47萬逐漸增加到71.18萬時,網(wǎng)格數(shù)對出口流量及壓力特性的計算結(jié)果影響很小。圖5a所示為泵內(nèi)高壓側(cè)的速度監(jiān)測點設置,通過速度監(jiān)測,并對41.47萬、55.18萬、71.18萬網(wǎng)格數(shù)進行計算,得出速度監(jiān)測曲線如圖5b所示。結(jié)果表明,41.47萬網(wǎng)格數(shù)與55.18萬、71.18萬網(wǎng)格數(shù)計算結(jié)果的最大誤差分別為2.32%、2.99%,即網(wǎng)格數(shù)對速度分布不敏感。這說明當轉(zhuǎn)子腔計算域網(wǎng)格數(shù)在41.47萬左右時,計算結(jié)果趨于穩(wěn)定,因此最終確定5葉轉(zhuǎn)子的網(wǎng)格數(shù)為41.47萬。按照相同的網(wǎng)格劃分設置方法,分別對3葉、4葉、6葉轉(zhuǎn)子腔模型進行網(wǎng)格劃分,最終得到3葉、4葉、6葉轉(zhuǎn)子的網(wǎng)格數(shù)分別為38.99萬、40.80萬、42.11萬。

    圖4 出口流量、壓力特性驗證

    (a)腔內(nèi)速度監(jiān)測點設置

    (b)腔內(nèi)速度驗證圖5 轉(zhuǎn)子腔內(nèi)速度監(jiān)測點設置及速度驗證

    2.2 數(shù)值計算方法

    基于動網(wǎng)格技術(shù)和RNGk-ε湍流模型,采用SIMPLES算法求解基本方程,固壁面采用無滑移壁面,附近壁面采用標準壁面函數(shù),并采用一階迎風離散格式。該模型的湍動能k和耗散率ε的輸運方程可表示為

    Gk+Gb-ρε-YM+Sk

    (12)

    (13)

    式中:Gk是由于平均速度梯度引起的湍動能k的生成項;Gb是由于浮力引起的湍動能k的生成項;YM為可壓縮湍流中脈動擴張項;αk和αε為與湍動能k和湍流耗散率ε對應的Prandtl數(shù)。

    3 實驗結(jié)果分析

    為了驗證數(shù)值方法的合理性,選取5葉轉(zhuǎn)子(轉(zhuǎn)子間隙0.15 mm)的凸輪式氣體循環(huán)泵進行實驗驗證,對4種不同進出口壓差條件下的泵腔內(nèi)流量進行測試,圖6所示為實驗系統(tǒng)圖。

    1—氣罐;2—手閥;3—熱式流量計;4—減壓閥;5—單向閥;6—溫度計;7—加濕器;8—壓力測量計;9—氣體循環(huán)泵;10—電動調(diào)節(jié)閥;11—空壓機;12—壓力表。圖6 氣體循環(huán)泵實驗系統(tǒng)

    (a)出口流量-進出口壓差曲線

    (b)容積效率-進出口壓差曲線圖7 氣體循環(huán)泵數(shù)值預測和實驗結(jié)果對比

    圖7為氣體循環(huán)泵數(shù)值預測與氣動性能實驗結(jié)果的對比。轉(zhuǎn)速為3 000 r·min-1時,通過調(diào)節(jié)球閥開度,選取4種進出口壓差,測量泵出口的實際流量。結(jié)果表明,隨進出口壓差的增大,氣體循環(huán)泵出口流量和容積效率均逐漸減小,數(shù)值預測和實驗結(jié)果誤差在3%以內(nèi)。產(chǎn)生誤差的主要原因是數(shù)值預測中忽略了端面間隙產(chǎn)生的軸向泄漏及機械損失,其相對誤差在合理范圍,能較為準確地預測泵內(nèi)氣動性能。

    4 CFD數(shù)值計算結(jié)果分析

    4.1 轉(zhuǎn)子間隙對凸輪式氣體循環(huán)泵多葉轉(zhuǎn)子出口流量特性的影響

    圖8a為出口平均流量與轉(zhuǎn)子間隙、轉(zhuǎn)子葉數(shù)的關(guān)系,圖8b為容積效率與轉(zhuǎn)子間隙、轉(zhuǎn)子葉數(shù)的關(guān)系。容積效率為

    (14)

    式中:Qop為單位時間內(nèi)實際泵出口平均流量;Qt為單位時間內(nèi)理論泵出口平均流量;q為轉(zhuǎn)子每轉(zhuǎn)一圈泵理論排除介質(zhì)的流量,即排量。

    由圖8可以看出,轉(zhuǎn)子間隙、轉(zhuǎn)子葉數(shù)對氣體循環(huán)泵出口流量特性、容積效率有顯著影響,隨著轉(zhuǎn)子間隙由0.1 mm增大到0.3 mm,不同葉數(shù)的泵出口平均流量和容積效率均呈下降趨勢。這是由于轉(zhuǎn)子葉數(shù)增多使轉(zhuǎn)子容積占比增加,泵腔有效容積減小,進而導致容積效率下降,出口平均流量也隨之下降。3葉、4葉、5葉、6葉轉(zhuǎn)子的泵出口平均流量分別下降了0.009 4、0.007 9、0.006 4和0.006 2 kg·s-1,同時容積效率分別下降了28.601%、16.963%、14.674%和15.155%,且3葉轉(zhuǎn)子的出口平均流量和容積效率下降幅度最大,5葉轉(zhuǎn)子的容積效率下降幅度最小。這是由于隨著轉(zhuǎn)子間隙的增大,由高壓腔至低壓腔的直連泄漏通道變寬,單位時間內(nèi)氣體泄漏量逐漸增多;相反,轉(zhuǎn)子間隙越小,氣體泄漏越少,但容易對轉(zhuǎn)子造成磨損,降低泵的壽命,同時轉(zhuǎn)子加工難度增大。

    (a)出口平均流量與轉(zhuǎn)子間隙、轉(zhuǎn)子葉數(shù)的關(guān)系

    (b)容積效率與轉(zhuǎn)子間隙、轉(zhuǎn)子葉數(shù)的關(guān)系圖8 出口平均流量、容積效率與轉(zhuǎn)子間隙、轉(zhuǎn)子 葉數(shù)的關(guān)系

    圖9為轉(zhuǎn)子間隙、轉(zhuǎn)子葉數(shù)和流量脈動系數(shù)的關(guān)系。流量脈動系數(shù)通常用來表示流量脈動強度,定義流量脈動系數(shù)為

    (15)

    式中:Qmax為出口最大瞬時流量;Qmin為出口最小瞬時流量;Qave為出口平均流量。

    圖9 轉(zhuǎn)子間隙、轉(zhuǎn)子葉數(shù)和流量脈動系數(shù)的關(guān)系

    由圖9可以看出:轉(zhuǎn)子間隙、轉(zhuǎn)子葉數(shù)對流量脈動強度有顯著影響,隨著轉(zhuǎn)子間隙由0.1 mm增加到0.3 mm,不同葉數(shù)的流量脈動系數(shù)均呈現(xiàn)先減小后增大的趨勢,并且均存在一個最小值,3葉、4葉轉(zhuǎn)子在0.15 mm間隙時流量脈動系數(shù)最小,分別為3.016 8、1.749 1;5葉、6葉轉(zhuǎn)子在0.25 mm間隙時流量脈動系數(shù)最小,分別為1.435 9、1.931 7。結(jié)合圖8可以得到,3葉轉(zhuǎn)子的出口平均流量和流量脈動強度均大于其他葉數(shù)。這是由于葉數(shù)較少,轉(zhuǎn)子容積占比較小,就會使泵腔容積效率和出口平均流量增加,但同時轉(zhuǎn)子葉數(shù)過少會導致入口端與出口端之間的逐級串聯(lián)過渡腔室較少,進而相鄰腔室之間的脈動壓差增大,造成流量脈動強度增大。由于6葉轉(zhuǎn)子的平均流量小于其他葉數(shù),即Qave較小,導致6葉轉(zhuǎn)子的流量脈動系數(shù)較大。同時,可以看出5葉轉(zhuǎn)子的流量脈動系數(shù)要小于4葉、6葉,說明在介質(zhì)為空氣的條件下,5葉轉(zhuǎn)子能較好地抑制流量脈動強度。

    綜合以上分析,不同間隙會導致泵的出口平均流量、容積效率和脈動強度產(chǎn)生變化,在保證泵出口平均流量和容積效率的前提下,應使泵出口流量脈動強度穩(wěn)定在較低的數(shù)值,3葉、4葉轉(zhuǎn)子在0.15 mm間隙時流量脈動強度最小,且容積效率分別保持在75.8%、70.9%以上;5葉、6葉轉(zhuǎn)子在0.25 mm間隙時流量脈動強度最小,且容積效率分別保持在60.4%、54.9%以上。

    4.2 葉數(shù)對轉(zhuǎn)子腔流動及徑向激勵力的影響

    4.2.1 葉數(shù)對轉(zhuǎn)子腔內(nèi)流動的影響 通過選取相對脈動強度較低的0.15 mm間隙時3葉、4葉轉(zhuǎn)子,0.25 mm間隙時5葉、6葉轉(zhuǎn)子,進行對比分析葉數(shù)對凸輪式氣體循環(huán)泵腔內(nèi)流動及性能的影響。

    圖10 轉(zhuǎn)子葉數(shù)和泵出口流量的關(guān)系

    圖10所示為不同葉數(shù)的氣體循環(huán)泵一個周期內(nèi)的出口流量曲線,結(jié)果表明,不同葉數(shù)的氣體循環(huán)泵的出口瞬時流量均存在周期性變化,在一個周期內(nèi),泵出口流量脈動的次數(shù)是轉(zhuǎn)子葉數(shù)的2倍。同時,3葉、4葉、5葉轉(zhuǎn)子的氣體循環(huán)泵出口均存在不同程度的次級流量脈動,隨著葉數(shù)的增加,出口流量脈動強度和次級流量脈動強度逐漸衰減,當葉數(shù)達到6時,泵出口的次級流量脈動已完全消失,因此葉數(shù)的增加有效抑制了泵出口流量脈動強度和次級流量脈動強度。隨著葉數(shù)由3依次增加到6,泵出口平均流量分別下降了0.004 8、0.011和0.015 kg·s-1,并且葉數(shù)越多,相鄰葉數(shù)之間的泵出口平均流量下降幅度越小。同時,3葉轉(zhuǎn)子出現(xiàn)回流現(xiàn)象,由于3葉轉(zhuǎn)子的流量脈動強度較高,腔內(nèi)壓力上升,在壓力的作用下產(chǎn)生回流。由此說明,多葉轉(zhuǎn)子能較好地抑制流量脈動強度,改善內(nèi)部流動,但葉數(shù)越多,流量衰減越嚴重,因此應適當增加轉(zhuǎn)子的葉數(shù)來改善氣體循環(huán)泵出口的流量特性。

    圖11所示為不同轉(zhuǎn)子葉數(shù)的氣體循環(huán)泵旋轉(zhuǎn)一周的出口總壓曲線,結(jié)果表明:在一個周期內(nèi),泵出口瞬時總壓的脈動次數(shù)是轉(zhuǎn)子葉數(shù)的2倍,且3葉、4葉、5葉轉(zhuǎn)子的氣體循環(huán)泵出口均存在不同程度的次級壓力脈動,隨著葉數(shù)的增加,泵出口次級壓力脈動強度和出口瞬時壓力脈動強度都逐漸減小,當葉數(shù)達到6時,泵出口的次級壓力脈動已完全消失,并且波動較為平穩(wěn);隨著葉數(shù)由3依次增加到6,泵出口平均總壓分別下降了0.142、0.068和0.013 kPa,并且葉數(shù)越多,相鄰葉數(shù)之間的平均總壓下降幅度越小。由此說明,葉數(shù)增加可以減小出口平均總壓的下降趨勢,削弱泵出口壓力脈動,但同時會造成泵出口壓力降低,因此應適當增加葉數(shù)來改善泵出口壓力特性。

    圖11 轉(zhuǎn)子葉數(shù)和泵出口瞬時總壓的關(guān)系

    (a)3葉轉(zhuǎn)子腔內(nèi)壓力分布

    (b)4葉轉(zhuǎn)子腔內(nèi)壓力分布

    (c)5葉轉(zhuǎn)子腔內(nèi)壓力分布

    (d)6葉轉(zhuǎn)子腔內(nèi)壓力分布圖12 不同葉數(shù)轉(zhuǎn)子腔內(nèi)壓力分布

    圖12為不同葉數(shù)的氣體循環(huán)泵轉(zhuǎn)子腔內(nèi)壓力分布,結(jié)果表明,隨著葉數(shù)的增加,轉(zhuǎn)子腔內(nèi)的逐級串聯(lián)過渡腔室也逐漸增多,由于轉(zhuǎn)子與轉(zhuǎn)壁之間存在微小間隙,和出口處相連的高壓腔氣體在壓力作用下通過轉(zhuǎn)壁間隙泄漏到過渡腔室,但由于泄漏量一定,所以過渡腔室中的壓力范圍一定。通過對比圖12a~12d可以看出,轉(zhuǎn)子葉數(shù)增多,獨立腔室的數(shù)量也隨之增多,進而導致氣體由高壓側(cè)經(jīng)轉(zhuǎn)壁間隙泄漏到低壓側(cè)時,高壓側(cè)與低壓側(cè)之間的壓差和泄漏呈逐級遞減的趨勢,形成多個緩沖區(qū),因此多葉轉(zhuǎn)子能夠有效抑制泵在運轉(zhuǎn)過程中由于壓差產(chǎn)生的沖擊力,提高泵內(nèi)部流場的穩(wěn)定性,延長泵的使用壽命。

    圖13為沿模型外緣方向,自低壓側(cè)到高壓側(cè)依次設置19個壓力監(jiān)測點,由于不同葉數(shù)轉(zhuǎn)子的外圓直徑均為90 mm,因此不同葉數(shù)轉(zhuǎn)子外緣上監(jiān)測點的位置相同。

    圖13 壓力監(jiān)測點設置

    (a)t=0.085 s

    (b)t=0.09 s

    (c)t=0.095 s

    圖14為不同葉數(shù)在一個周期內(nèi)4個時刻的腔內(nèi)壓力梯度變化規(guī)律,結(jié)果表明,在同一時刻,不同葉數(shù)的轉(zhuǎn)子腔內(nèi)壓力均呈階梯性變化,隨著轉(zhuǎn)子葉數(shù)增多,轉(zhuǎn)子腔內(nèi)部增壓效果更快,壓力變化更為平緩。圖14a、14b中3葉轉(zhuǎn)子腔內(nèi)共出現(xiàn)兩次增壓,且第二次增壓更為急劇,因此會造成較大的壓力脈動。圖14c、14d中3葉轉(zhuǎn)子只出現(xiàn)一次增壓,這是因為轉(zhuǎn)子轉(zhuǎn)動的角度不同,只形成一個獨立腔室;當葉數(shù)增加到6時,轉(zhuǎn)子腔內(nèi)共出現(xiàn)4次較為明顯的增壓,且每次增壓的幅度較小,相鄰兩個獨立腔室之間壓差減小,因此6葉轉(zhuǎn)子的壓力脈動最小。由此說明,隨轉(zhuǎn)子葉數(shù)的增多,轉(zhuǎn)子腔內(nèi)增壓過程更為穩(wěn)定,即能量轉(zhuǎn)換的穩(wěn)定性更好。

    (d)t=0.1 s圖14 不同葉數(shù)轉(zhuǎn)子腔內(nèi)壓力梯度變化規(guī)律

    4.2.2 葉數(shù)對轉(zhuǎn)子徑向激勵力的影響 由于泵腔內(nèi)兩轉(zhuǎn)子呈對稱分布,且進行同步反向轉(zhuǎn)動,因此兩轉(zhuǎn)子所受的徑向激勵力大小相等,方向相反,本文取左轉(zhuǎn)子為研究對象進行分析。

    圖15 葉數(shù)對左轉(zhuǎn)子徑向激勵力的影響

    圖15為不同葉數(shù)的左轉(zhuǎn)子所受徑向激勵力矢量的分布規(guī)律,可以看出,隨著轉(zhuǎn)子葉數(shù)的逐漸增大,左轉(zhuǎn)子所受徑向激勵力逐漸減小,分布逐漸集中,且徑向激勵力在y軸的分量均為正方向,在x軸的分量大部分為負方向。3葉轉(zhuǎn)子各點的散亂程度明顯大于其他葉數(shù),且部分點出現(xiàn)在x軸正方向,即Fx的方向發(fā)生改變。當葉數(shù)增大到6時,各點的散亂程度降為最小,因此多葉轉(zhuǎn)子能明顯減小轉(zhuǎn)子所受徑向激勵力的大小,保證泵運行更加平穩(wěn)。另外,圖15中有少量徑向激勵力矢量離散點,這是由于轉(zhuǎn)子型線圓弧段與漸開線段并非完全光滑過渡,產(chǎn)生了徑向激勵力突變。

    圖16a所示為不同葉數(shù)的凸輪式氣體循環(huán)泵在一個周期內(nèi)轉(zhuǎn)子徑向激勵力x軸分量Fx的變化規(guī)律,結(jié)果表明,不同葉數(shù)轉(zhuǎn)子的Fx均呈周期性變化,隨著葉數(shù)的增加,轉(zhuǎn)子腔內(nèi)氣體對Fx的脈動幅值呈逐漸減小的趨勢。在一個周期內(nèi),波峰與波谷的數(shù)量與轉(zhuǎn)子葉數(shù)相等,說明其主要受到轉(zhuǎn)子與腔內(nèi)氣體相互作用的影響。與圖15所示一致,3葉轉(zhuǎn)子所受的力一部分為正值,而其他葉數(shù)的轉(zhuǎn)子所受的力均為負值。這是由于3葉轉(zhuǎn)子過渡腔室較少,過渡腔與出口處高壓腔之間的壓差較大,氣體出現(xiàn)回流現(xiàn)象,導致氣體對轉(zhuǎn)子出現(xiàn)相反的作用力。3葉轉(zhuǎn)子Fx的最大波峰值為76.07 N;6葉轉(zhuǎn)子Fx的最小波峰值為10.37 N,方向指向y軸負方向,比3葉轉(zhuǎn)子Fx的最大值減小86.37%。同時,隨著轉(zhuǎn)子葉數(shù)增多,Fx的二次脈動幅值也逐漸下降,脈動幅度逐漸平穩(wěn)。

    (a)Fx的脈動幅值

    (b)Fx的脈動頻域圖圖16 多葉轉(zhuǎn)子徑向激勵力x軸分量Fx的脈動特性

    圖16b所示為不同葉數(shù)的凸輪式氣體循環(huán)泵在一個周期內(nèi)轉(zhuǎn)子徑向激勵力x軸分量Fx的脈動頻域圖,結(jié)果表明,隨著葉數(shù)的增多,徑向激勵力脈動主頻出現(xiàn)的位置逐漸后移,主頻最大振幅出現(xiàn)在Z=3時,且隨著轉(zhuǎn)子葉數(shù)的增多主頻最大振幅不斷減小。小葉數(shù)會出現(xiàn)相對明顯的高頻脈動,這是由于轉(zhuǎn)子葉數(shù)較少,腔內(nèi)有效容積增大,更容易出現(xiàn)流動分離和二次流動現(xiàn)象。由此可見,轉(zhuǎn)子葉數(shù)增多能有效抑制Fx的脈動強度。

    圖17a所示為不同葉數(shù)的凸輪式氣體循環(huán)泵在一個周期內(nèi)轉(zhuǎn)子徑向激勵力y軸分量Fy的變化規(guī)律,結(jié)果表明,不同葉數(shù)轉(zhuǎn)子的Fy均呈周期性變化,在一個周期內(nèi),波峰與波谷的數(shù)量與轉(zhuǎn)子葉數(shù)相等。3葉轉(zhuǎn)子Fy的最大波峰值為212.11 N;6葉轉(zhuǎn)子Fy的最小波峰值為191.55 N,比3葉轉(zhuǎn)子Fy的最大值減小9.69%。除5葉轉(zhuǎn)子外,隨著轉(zhuǎn)子葉數(shù)的增加,轉(zhuǎn)子腔內(nèi)高壓側(cè)氣體對Fy的脈動幅值呈逐漸減小的趨勢,而5葉轉(zhuǎn)子Fy的最大波峰值為207.41 N,僅次于3葉轉(zhuǎn)子,但各個葉數(shù)轉(zhuǎn)子Fy的峰值差較小,3葉、6葉相差最大為9.69%,3葉、5葉相差最小為2.27%。當Z=3~6時,Z對Fy影響較小。

    (a)Fy的脈動幅值

    圖17b所示為不同葉數(shù)的凸輪式氣體循環(huán)泵在一個周期內(nèi)轉(zhuǎn)子徑向激勵力y軸分量Fy的脈動頻域圖,可以看出,Fy的脈動主頻出現(xiàn)在250 Hz附近,且主頻最大振幅出現(xiàn)在Z=3時,隨著轉(zhuǎn)子葉數(shù)增多,主頻最大振幅逐漸減小,但脈動強度的變化較小,與圖16b相比,轉(zhuǎn)子葉數(shù)對Fy脈動強度的抑制不顯著。

    (b)Fy的脈動頻域圖圖17 多葉轉(zhuǎn)子徑向激勵力y軸分量Fy的脈動特性

    5 結(jié) 論

    采用數(shù)值模擬方法,分析了轉(zhuǎn)子不同間隙、不同葉數(shù)條件下凸輪式氣體循環(huán)泵的氣動性能和徑向激勵力的變化規(guī)律,得出以下結(jié)論。

    (1)隨著轉(zhuǎn)子間隙由0.1 mm逐漸增大到0.3 mm,不同葉數(shù)的泵出口平均流量和容積效率均逐漸下降,且隨著轉(zhuǎn)子葉數(shù)增多,下降幅度逐漸減小,同時轉(zhuǎn)子腔獨立腔室增多,相鄰兩腔室之間的壓差減小,當間隙增大時,多葉轉(zhuǎn)子相比3葉轉(zhuǎn)子的泄漏幅度降低;隨著轉(zhuǎn)子間隙增大,不同葉數(shù)轉(zhuǎn)子的流量脈動系數(shù)均呈先減小后增大的趨勢,同時5葉轉(zhuǎn)子的流量脈動系數(shù)最小,說明介質(zhì)為空氣時,5葉轉(zhuǎn)子能較好地抑制流量脈動強度。

    (2)轉(zhuǎn)子葉數(shù)增多會降低泵出口平均流量和出口總壓,但能夠有效削弱泵出口流量脈動強度和壓力脈動強度。隨著轉(zhuǎn)子葉數(shù)增多,逐級串聯(lián)過渡腔室增多,壓差逐級遞減,形成多個緩沖區(qū),降低了泵出口流量脈動強度和壓力脈動強度;同時,葉數(shù)增多,轉(zhuǎn)子腔內(nèi)部增壓效果更平穩(wěn),即能量轉(zhuǎn)換的穩(wěn)定性更好。因此,在保證泵的流量指標和壓力指標的前提下,應適當增多轉(zhuǎn)子葉數(shù)來抑制流量脈動和壓力脈動,保證泵的穩(wěn)定運行。

    (3)隨著轉(zhuǎn)子葉數(shù)增多,轉(zhuǎn)子徑向激勵力矢量逐漸減小,分布逐漸集中,同時Fx的大小和脈動強度均顯著降低,而轉(zhuǎn)子葉數(shù)對Fy的影響不顯著,各個葉數(shù)所受Fy的峰值差較小。

    猜你喜歡
    葉數(shù)循環(huán)泵子葉
    皖南煙區(qū)烤煙新品種中川208的適宜留葉數(shù)研究
    脫硫系統(tǒng)用漿液循環(huán)泵汽蝕分析
    煙草株高和葉數(shù)性狀QTL定位及候選基因預測
    脫硫裝置水洗循環(huán)泵機封沖洗案例分析
    摘葉對不同株型木薯品種產(chǎn)量和氮磷鉀素積累分配的影響
    脫硫裝置貧液循環(huán)泵優(yōu)化改造研究
    懸鈴木幼苗的初生維管系統(tǒng)演化結(jié)構(gòu)研究
    熱源循環(huán)泵流量對熱電聯(lián)供系統(tǒng)優(yōu)化運行影響
    煤氣與熱力(2021年2期)2021-03-19 08:55:50
    黑山羊胎盤子葉性狀結(jié)構(gòu)與繁殖性能的相關(guān)性
    留葉數(shù)和成熟度影響雪茄煙葉品質(zhì)
    亚洲三级黄色毛片| 国产一级毛片在线| 久久狼人影院| 我要看黄色一级片免费的| 国产日韩一区二区三区精品不卡| 精品一区二区免费观看| 观看av在线不卡| freevideosex欧美| 亚洲精品自拍成人| 日产精品乱码卡一卡2卡三| 男女免费视频国产| 久久免费观看电影| 久久99一区二区三区| 亚洲av日韩在线播放| 久久久久精品人妻al黑| 尾随美女入室| 久久久国产欧美日韩av| 看非洲黑人一级黄片| 纵有疾风起免费观看全集完整版| 99国产精品免费福利视频| 男女午夜视频在线观看 | 97在线视频观看| 99久久中文字幕三级久久日本| 一区二区三区精品91| 国产白丝娇喘喷水9色精品| 国产xxxxx性猛交| 女人被躁到高潮嗷嗷叫费观| 五月天丁香电影| 黄色毛片三级朝国网站| 欧美另类一区| 极品人妻少妇av视频| 一级,二级,三级黄色视频| 亚洲四区av| 免费黄网站久久成人精品| 最近手机中文字幕大全| 中文字幕最新亚洲高清| 我的女老师完整版在线观看| 人体艺术视频欧美日本| 日韩在线高清观看一区二区三区| av黄色大香蕉| 黄色 视频免费看| 欧美日韩av久久| 大片免费播放器 马上看| 日本猛色少妇xxxxx猛交久久| 18禁裸乳无遮挡动漫免费视频| 日韩制服骚丝袜av| 日本色播在线视频| 久久99蜜桃精品久久| 香蕉丝袜av| 黄色一级大片看看| 18禁国产床啪视频网站| 久久毛片免费看一区二区三区| 丁香六月天网| 丝袜在线中文字幕| 最近2019中文字幕mv第一页| 2021少妇久久久久久久久久久| 黄色毛片三级朝国网站| 国产av国产精品国产| 国产亚洲欧美精品永久| 99久久中文字幕三级久久日本| 免费播放大片免费观看视频在线观看| 草草在线视频免费看| 黑人高潮一二区| 日韩成人av中文字幕在线观看| 精品少妇内射三级| 欧美国产精品va在线观看不卡| 国产片特级美女逼逼视频| 9热在线视频观看99| 秋霞在线观看毛片| av黄色大香蕉| 国产一区有黄有色的免费视频| 天天影视国产精品| 18+在线观看网站| 在线天堂中文资源库| 国产有黄有色有爽视频| 亚洲精品aⅴ在线观看| 看十八女毛片水多多多| 国产在线视频一区二区| 国产成人一区二区在线| 人妻 亚洲 视频| 欧美97在线视频| 国产男人的电影天堂91| 美女大奶头黄色视频| 男人舔女人的私密视频| 亚洲精品成人av观看孕妇| 最近2019中文字幕mv第一页| 国产精品嫩草影院av在线观看| 亚洲欧美清纯卡通| 在线观看免费高清a一片| 国产日韩一区二区三区精品不卡| 伦理电影免费视频| xxx大片免费视频| 99久久中文字幕三级久久日本| av播播在线观看一区| 精品福利永久在线观看| 午夜福利在线观看免费完整高清在| 国产精品久久久久久精品电影小说| 日韩 亚洲 欧美在线| a 毛片基地| 国产日韩欧美视频二区| av在线app专区| 中文字幕人妻丝袜制服| 国产在线视频一区二区| 国产又色又爽无遮挡免| 久久这里有精品视频免费| 免费少妇av软件| 男女无遮挡免费网站观看| 久热久热在线精品观看| 亚洲国产精品成人久久小说| 午夜福利影视在线免费观看| 日本欧美国产在线视频| 老司机影院成人| 国产精品久久久久久精品电影小说| 欧美xxxx性猛交bbbb| 麻豆乱淫一区二区| 免费黄色在线免费观看| 丝袜脚勾引网站| 97精品久久久久久久久久精品| 亚洲av免费高清在线观看| 亚洲,一卡二卡三卡| 国产免费一区二区三区四区乱码| 午夜老司机福利剧场| 日本欧美国产在线视频| 在现免费观看毛片| 精品少妇久久久久久888优播| 国产成人免费无遮挡视频| 久久久精品94久久精品| 人人妻人人澡人人看| 久久久久精品久久久久真实原创| 少妇人妻 视频| 女性生殖器流出的白浆| 黄色一级大片看看| 伊人久久国产一区二区| 晚上一个人看的免费电影| 久久久久久人人人人人| 午夜精品国产一区二区电影| 黄色 视频免费看| 韩国精品一区二区三区 | 久久久精品94久久精品| 成人18禁高潮啪啪吃奶动态图| 丰满迷人的少妇在线观看| 中文字幕最新亚洲高清| 美女xxoo啪啪120秒动态图| 最近中文字幕高清免费大全6| av国产精品久久久久影院| 高清av免费在线| 国产高清国产精品国产三级| 成人18禁高潮啪啪吃奶动态图| 黑人欧美特级aaaaaa片| 9色porny在线观看| av在线播放精品| 美女内射精品一级片tv| 热re99久久国产66热| 亚洲精品国产av成人精品| 国内精品宾馆在线| 午夜激情av网站| 亚洲五月色婷婷综合| 国产精品久久久久久久久免| 高清毛片免费看| 国产亚洲一区二区精品| 三级国产精品片| 中文字幕制服av| 男人舔女人的私密视频| 18禁观看日本| 在线观看人妻少妇| 免费在线观看完整版高清| 成人毛片60女人毛片免费| 亚洲精品乱久久久久久| 中文天堂在线官网| 久久精品国产综合久久久 | 日韩不卡一区二区三区视频在线| 狂野欧美激情性xxxx在线观看| 一本大道久久a久久精品| 人妻少妇偷人精品九色| 精品亚洲乱码少妇综合久久| 熟女av电影| 国产福利在线免费观看视频| 国产 一区精品| 午夜激情久久久久久久| 亚洲成人手机| 国产精品久久久久久久久免| 精品一品国产午夜福利视频| 成人免费观看视频高清| 亚洲精品自拍成人| 一本—道久久a久久精品蜜桃钙片| 午夜免费鲁丝| 美女国产高潮福利片在线看| 韩国精品一区二区三区 | 少妇的丰满在线观看| 观看美女的网站| 最近的中文字幕免费完整| 高清在线视频一区二区三区| 伦精品一区二区三区| 在线看a的网站| 日日撸夜夜添| 美女内射精品一级片tv| 18禁观看日本| 蜜桃国产av成人99| 亚洲欧美清纯卡通| 久久久久久人妻| 午夜福利乱码中文字幕| 国产成人精品久久久久久| 婷婷色麻豆天堂久久| 18在线观看网站| 日韩一本色道免费dvd| 男女啪啪激烈高潮av片| 国产黄频视频在线观看| 岛国毛片在线播放| 国产一区二区激情短视频 | 日本免费在线观看一区| 两个人免费观看高清视频| 国产色爽女视频免费观看| 男人舔女人的私密视频| 亚洲成人一二三区av| 香蕉精品网在线| 激情视频va一区二区三区| 亚洲五月色婷婷综合| 最近2019中文字幕mv第一页| 久久国产精品男人的天堂亚洲 | 美女福利国产在线| 国产成人欧美| 午夜福利,免费看| 18在线观看网站| 午夜久久久在线观看| 国产精品熟女久久久久浪| 男女午夜视频在线观看 | 91久久精品国产一区二区三区| 老司机影院成人| 国产免费又黄又爽又色| 最近中文字幕2019免费版| 亚洲欧美成人综合另类久久久| 国产综合精华液| 亚洲精品456在线播放app| 日本黄色日本黄色录像| 国产极品天堂在线| 亚洲欧洲精品一区二区精品久久久 | 午夜激情久久久久久久| 五月天丁香电影| 久久国产精品男人的天堂亚洲 | 夜夜骑夜夜射夜夜干| 中文字幕制服av| 大香蕉久久网| 99香蕉大伊视频| 国产成人av激情在线播放| 黄色 视频免费看| 国产成人一区二区在线| 日本色播在线视频| 免费大片18禁| 久久久精品区二区三区| 十八禁网站网址无遮挡| av福利片在线| 一本大道久久a久久精品| 国产亚洲av片在线观看秒播厂| 国产乱人偷精品视频| 欧美亚洲日本最大视频资源| 免费黄频网站在线观看国产| 日本黄大片高清| 亚洲欧洲精品一区二区精品久久久 | 精品久久蜜臀av无| 久久国产亚洲av麻豆专区| 美女国产高潮福利片在线看| 久久久精品94久久精品| 美女国产高潮福利片在线看| 考比视频在线观看| 亚洲激情五月婷婷啪啪| 丰满乱子伦码专区| a级毛色黄片| 少妇高潮的动态图| 制服诱惑二区| 搡老乐熟女国产| 毛片一级片免费看久久久久| videos熟女内射| 成年女人在线观看亚洲视频| 大香蕉久久网| 亚洲第一区二区三区不卡| 大码成人一级视频| 97人妻天天添夜夜摸| kizo精华| 亚洲熟女精品中文字幕| 一区二区三区四区激情视频| 久久久久久人妻| 亚洲成人av在线免费| 亚洲综合色网址| 一级黄片播放器| 欧美日本中文国产一区发布| 99热国产这里只有精品6| 夜夜骑夜夜射夜夜干| 亚洲精品一二三| 丝袜喷水一区| 香蕉丝袜av| 精品久久蜜臀av无| 人妻少妇偷人精品九色| 免费看不卡的av| 观看av在线不卡| 夜夜骑夜夜射夜夜干| 国产色婷婷99| a级毛片在线看网站| 黄片播放在线免费| 日韩制服丝袜自拍偷拍| 久久人妻熟女aⅴ| 国产片内射在线| 久久久久久久国产电影| 五月开心婷婷网| 十分钟在线观看高清视频www| 熟女人妻精品中文字幕| 亚洲人成网站在线观看播放| 精品国产乱码久久久久久小说| 18禁在线无遮挡免费观看视频| 久久婷婷青草| 22中文网久久字幕| 中国国产av一级| 久久久久久久大尺度免费视频| 国产高清不卡午夜福利| av免费在线看不卡| 免费日韩欧美在线观看| 熟女av电影| 亚洲欧美成人综合另类久久久| 国产片特级美女逼逼视频| 久久av网站| 少妇人妻 视频| 女性被躁到高潮视频| 国产精品不卡视频一区二区| av福利片在线| 欧美国产精品va在线观看不卡| 国产又色又爽无遮挡免| 天堂8中文在线网| 日本-黄色视频高清免费观看| 18禁裸乳无遮挡动漫免费视频| 女性被躁到高潮视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 97在线人人人人妻| 母亲3免费完整高清在线观看 | 久久久久视频综合| 一边摸一边做爽爽视频免费| 99re6热这里在线精品视频| 日本wwww免费看| a 毛片基地| 岛国毛片在线播放| 精品久久蜜臀av无| 三级国产精品片| 最近中文字幕2019免费版| 亚洲欧洲日产国产| 飞空精品影院首页| 最近中文字幕2019免费版| 丰满乱子伦码专区| 免费黄网站久久成人精品| 天天操日日干夜夜撸| 欧美日韩视频高清一区二区三区二| 丁香六月天网| av网站免费在线观看视频| 一区二区三区四区激情视频| 一个人免费看片子| 中文天堂在线官网| 91在线精品国自产拍蜜月| 免费黄频网站在线观看国产| 视频在线观看一区二区三区| 美女福利国产在线| 午夜免费男女啪啪视频观看| 久久久久久伊人网av| 水蜜桃什么品种好| 黄色一级大片看看| 免费不卡的大黄色大毛片视频在线观看| 水蜜桃什么品种好| 国产在线免费精品| 在线观看三级黄色| 性色av一级| 久久国产精品大桥未久av| 一边亲一边摸免费视频| 99久国产av精品国产电影| 久久这里只有精品19| 国产精品久久久久成人av| 国产成人午夜福利电影在线观看| 日韩中字成人| 国产精品嫩草影院av在线观看| 亚洲性久久影院| 我的女老师完整版在线观看| 日韩 亚洲 欧美在线| 亚洲,欧美,日韩| 看非洲黑人一级黄片| 在现免费观看毛片| 国产在线免费精品| 国产视频首页在线观看| 天天躁夜夜躁狠狠躁躁| 男女下面插进去视频免费观看 | av黄色大香蕉| 久久99蜜桃精品久久| 一本久久精品| 久久av网站| 啦啦啦在线观看免费高清www| 成人二区视频| 乱码一卡2卡4卡精品| 国产免费福利视频在线观看| 9色porny在线观看| 在线亚洲精品国产二区图片欧美| 免费观看在线日韩| 久久影院123| 国产亚洲精品第一综合不卡 | 亚洲欧美一区二区三区国产| 国产成人精品福利久久| 一级片'在线观看视频| 18禁动态无遮挡网站| 交换朋友夫妻互换小说| 高清不卡的av网站| 日韩精品免费视频一区二区三区 | 免费观看无遮挡的男女| av在线老鸭窝| 国产精品国产三级专区第一集| 五月天丁香电影| 国产探花极品一区二区| 高清av免费在线| 熟妇人妻不卡中文字幕| 日韩av不卡免费在线播放| 久久午夜福利片| 尾随美女入室| 欧美精品国产亚洲| 美女视频免费永久观看网站| 国产福利在线免费观看视频| 91精品国产国语对白视频| 夜夜爽夜夜爽视频| 国产熟女午夜一区二区三区| 亚洲人与动物交配视频| 国产精品.久久久| 欧美少妇被猛烈插入视频| 99久久人妻综合| 欧美日韩成人在线一区二区| 我要看黄色一级片免费的| 国产免费现黄频在线看| 成人影院久久| av在线app专区| 五月伊人婷婷丁香| 亚洲欧美中文字幕日韩二区| 九九在线视频观看精品| 亚洲精品久久午夜乱码| 国产午夜精品一二区理论片| 日韩大片免费观看网站| 日韩制服丝袜自拍偷拍| 超色免费av| 在线 av 中文字幕| 色网站视频免费| 国产av国产精品国产| kizo精华| 在线观看三级黄色| 亚洲伊人色综图| 欧美亚洲 丝袜 人妻 在线| www日本在线高清视频| 性高湖久久久久久久久免费观看| 国产精品免费大片| 老熟女久久久| 国产精品国产三级专区第一集| 热99国产精品久久久久久7| 久久鲁丝午夜福利片| 亚洲一级一片aⅴ在线观看| 成年人午夜在线观看视频| 欧美激情国产日韩精品一区| 男女高潮啪啪啪动态图| 国产成人精品一,二区| 一级毛片 在线播放| 亚洲伊人色综图| 欧美bdsm另类| 欧美3d第一页| 视频在线观看一区二区三区| 日本欧美国产在线视频| 亚洲成色77777| 亚洲色图综合在线观看| 蜜桃在线观看..| 免费久久久久久久精品成人欧美视频 | 国产成人免费无遮挡视频| 91精品三级在线观看| 男女高潮啪啪啪动态图| 成年女人在线观看亚洲视频| 中国三级夫妇交换| 亚洲精品国产av蜜桃| 免费观看a级毛片全部| 国产女主播在线喷水免费视频网站| 欧美丝袜亚洲另类| 黄色怎么调成土黄色| 久久婷婷青草| 亚洲综合色惰| 婷婷成人精品国产| 日韩av在线免费看完整版不卡| 国内精品宾馆在线| 99精国产麻豆久久婷婷| av在线播放精品| www.av在线官网国产| 视频中文字幕在线观看| 亚洲,欧美精品.| 91久久精品国产一区二区三区| 久热这里只有精品99| 精品午夜福利在线看| 久久久国产欧美日韩av| 亚洲天堂av无毛| 国产精品久久久久久久久免| 边亲边吃奶的免费视频| 欧美人与性动交α欧美软件 | √禁漫天堂资源中文www| 日本91视频免费播放| 一级毛片 在线播放| 亚洲伊人色综图| 国产麻豆69| 香蕉精品网在线| 免费观看a级毛片全部| 亚洲精品456在线播放app| 国产精品无大码| 欧美3d第一页| 国产精品一国产av| 国产精品国产三级国产专区5o| 1024视频免费在线观看| 97精品久久久久久久久久精品| 久久久久久久久久久久大奶| 亚洲国产色片| 欧美精品人与动牲交sv欧美| av女优亚洲男人天堂| 高清黄色对白视频在线免费看| 2021少妇久久久久久久久久久| 国产有黄有色有爽视频| 欧美人与性动交α欧美软件 | 午夜老司机福利剧场| 99精国产麻豆久久婷婷| 夜夜骑夜夜射夜夜干| 国产成人精品无人区| 97人妻天天添夜夜摸| 欧美精品一区二区大全| 51国产日韩欧美| 久久精品国产自在天天线| 99热全是精品| 国产精品一区二区在线观看99| 啦啦啦视频在线资源免费观看| 精品酒店卫生间| 亚洲精品久久午夜乱码| 亚洲精品一二三| 搡老乐熟女国产| 欧美激情 高清一区二区三区| 黑人欧美特级aaaaaa片| 免费看光身美女| 国产成人免费无遮挡视频| 日韩熟女老妇一区二区性免费视频| 啦啦啦视频在线资源免费观看| 国产成人精品福利久久| 亚洲精品久久午夜乱码| 国产在线免费精品| 王馨瑶露胸无遮挡在线观看| 最近中文字幕2019免费版| 一本久久精品| 成年美女黄网站色视频大全免费| 久久久久网色| 伊人亚洲综合成人网| 久久久久精品性色| 一级,二级,三级黄色视频| a级毛色黄片| 国产精品无大码| 中文字幕制服av| 日韩制服骚丝袜av| 又黄又爽又刺激的免费视频.| 国产国拍精品亚洲av在线观看| 秋霞伦理黄片| 国产精品女同一区二区软件| 十八禁网站网址无遮挡| 91精品伊人久久大香线蕉| 99精国产麻豆久久婷婷| 国产精品嫩草影院av在线观看| 22中文网久久字幕| 韩国高清视频一区二区三区| 看非洲黑人一级黄片| 精品人妻一区二区三区麻豆| 啦啦啦在线观看免费高清www| 人妻 亚洲 视频| 中国国产av一级| 十八禁网站网址无遮挡| 久久国产亚洲av麻豆专区| 一级毛片黄色毛片免费观看视频| 男男h啪啪无遮挡| 人妻少妇偷人精品九色| 这个男人来自地球电影免费观看 | 交换朋友夫妻互换小说| 国产乱来视频区| 蜜臀久久99精品久久宅男| 欧美成人精品欧美一级黄| 日本色播在线视频| 亚洲图色成人| 亚洲精品色激情综合| 成年动漫av网址| 国产精品免费大片| 国产又色又爽无遮挡免| 一边亲一边摸免费视频| 十八禁高潮呻吟视频| 国产精品国产三级国产av玫瑰| 欧美变态另类bdsm刘玥| 成人18禁高潮啪啪吃奶动态图| 亚洲精品一二三| 日韩成人av中文字幕在线观看| 人妻 亚洲 视频| 五月玫瑰六月丁香| 成年动漫av网址| 9色porny在线观看| 久热这里只有精品99| 亚洲综合色网址| kizo精华| 卡戴珊不雅视频在线播放| 久久精品人人爽人人爽视色| 国产精品熟女久久久久浪| 精品福利永久在线观看| 国产乱来视频区| 观看av在线不卡| 老女人水多毛片| 亚洲欧洲精品一区二区精品久久久 | 日本免费在线观看一区| a 毛片基地| av播播在线观看一区| 免费黄色在线免费观看| 最黄视频免费看| 涩涩av久久男人的天堂| 有码 亚洲区| 国产探花极品一区二区| 国产在视频线精品| 视频在线观看一区二区三区|