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

    葉輪口環(huán)結(jié)構(gòu)對離心泵性能及流場的影響

    2017-07-21 05:09:08牟介剛代東順谷云慶劉劍吳登昊馬藝
    關(guān)鍵詞:結(jié)構(gòu)

    牟介剛,代東順,谷云慶,劉劍,吳登昊,馬藝

    ?

    葉輪口環(huán)結(jié)構(gòu)對離心泵性能及流場的影響

    牟介剛1,代東順1,谷云慶1,劉劍1,吳登昊2,馬藝1

    (1. 浙江工業(yè)大學(xué)機(jī)械工程學(xué)院,浙江杭州,310014;2. 浙江工業(yè)大學(xué)之江學(xué)院,浙江杭州,310024)

    為了揭示葉輪口環(huán)間隙及結(jié)構(gòu)型式對離心泵性能的影響,以IS80-65-160離心泵為模型載體,采用RNG?湍流模型,對不同葉輪口環(huán)間隙及口環(huán)結(jié)構(gòu)型式下的離心泵內(nèi)部流場進(jìn)行數(shù)值模擬,分析葉輪口環(huán)間隙對離心泵的效率、壓力場和速度場的影響規(guī)律;研究3種環(huán)狀結(jié)構(gòu)型式的葉輪口環(huán)對離心泵性能的影響機(jī)制。研究結(jié)果表明:葉輪口環(huán)間隙的變化改變了離心泵內(nèi)部流場的流動狀態(tài),影響其前、后泵腔處的壓力分布,間隙內(nèi)部的速度發(fā)生明顯變化;鋸齒形口環(huán)泄漏量最多可減少16.2%,容積效率提高,軸向力降低,鋸齒形葉輪口環(huán)結(jié)構(gòu)可提升離心泵2%的水力效率;流經(jīng)兩環(huán)狀結(jié)構(gòu)空腔的流體,過流面積增加,產(chǎn)生劇烈旋渦,并進(jìn)行能量耗散,阻止了流體壓力能的恢復(fù),分散了經(jīng)密封齒隙高速射出的流體的動能,間隙內(nèi)部流體的速度梯度和壓力梯度增大,湍動狀態(tài)加劇,流動阻力增加。

    離心泵;流體動力學(xué);葉輪口環(huán);結(jié)構(gòu)型式;數(shù)值模擬

    離心泵是一種應(yīng)用極廣、種類甚多的通用機(jī)械,廣泛應(yīng)用于石油化工、動力工業(yè)、采礦和造船工業(yè)等。作為主要的輸送設(shè)備之一,對其性能及穩(wěn)定性的要求也在逐步提高。在把機(jī)械能轉(zhuǎn)化為液體能量的過程中不可避免地伴有機(jī)械損失,沖擊、脫流、速度方向急劇變化等引起的水力損失,以及液體從泵腔經(jīng)葉輪密封環(huán)間隙向葉輪進(jìn)口逆流產(chǎn)生的容積損失[1]等,造成能量浪費(fèi)。口環(huán)間隙密封是離心泵葉輪密封最常用的一種形式,口環(huán)間隙和結(jié)構(gòu)會對離心泵的容積損失產(chǎn)生直接作用,并且改變了離心泵內(nèi)部流場的流動狀態(tài),繼而對離心泵的整機(jī)性能、軸向力、徑向力等產(chǎn)生重要的影響。間隙流動作為目前流體機(jī)械研究的熱點(diǎn)問題,離心泵口環(huán)間隙處流動狀態(tài)復(fù)雜,口環(huán)間隙小,葉輪口環(huán)處的流動損失大[2]。已有很多學(xué)者對此進(jìn)行了研究,如:馬文生等[3]運(yùn)用數(shù)值方法對轉(zhuǎn)子?軸承密封系統(tǒng)動力學(xué)模型求解,研究了不同的密封間隙、密封直徑、壓差、入口損失率和密封長度對泄漏量和臨界轉(zhuǎn)速的影響,得到了泄漏量與密封間隙、密封直徑和密封長度之間的關(guān)系。王煒哲等[4]針對迷宮密封內(nèi)作渦動旋轉(zhuǎn)的轉(zhuǎn)子承受腔室內(nèi)周向不均勻壓力,從而導(dǎo)致振動特性發(fā)生變化的問題進(jìn)行了探討。由于口環(huán)間隙尺寸的限制,計(jì)算網(wǎng)格劃分困難,在現(xiàn)有的研究中一般只對葉輪、蝸殼等流體域進(jìn)行研究而忽略口環(huán)間隙、結(jié)構(gòu)變化對離心泵內(nèi)部流場的影響,導(dǎo)致預(yù)測的效率與試驗(yàn)值之間存在較大的偏差,對計(jì)算模型的簡化是造成偏差的主要原因[5?6]。計(jì)算流體力學(xué)的發(fā)展,為分析離心泵的水力性能和改進(jìn)泵的設(shè)計(jì)提供了有力的工具并且得到了廣泛的應(yīng)用[7?10]。研究發(fā)現(xiàn),在離心泵的高效區(qū)(設(shè)計(jì)流量的0.7~1.2倍)預(yù)測值和試驗(yàn)值基本吻合[11?12]。為了進(jìn)一步研究葉輪口環(huán)間隙及結(jié)構(gòu)的變化對離心泵性能的影響,本文作者采用數(shù)值模擬的方法,建立離心泵全流體域計(jì)算模型,采用商用Fluent軟件,分析了3種口環(huán)間隙和3種不同口環(huán)結(jié)構(gòu)對離心泵泄漏量、軸向力、水力效率等特性的影響,揭示了離心泵內(nèi)部流場的流動規(guī)律。

    1 數(shù)值模擬方法

    1.1 計(jì)算模型

    以IS80-65-160離心泵為模型載體,計(jì)算模型如圖1所示。計(jì)算模型主要包括進(jìn)水管、出水管、葉輪、蝸殼、前后泵腔6個(gè)部分,其主要參數(shù)為:流量=50 m3/h,揚(yáng)程=32 m,轉(zhuǎn)速=2 900 r/min,進(jìn)口直徑1=80 mm,葉輪外徑2=168 mm,葉片數(shù)為4。葉輪口環(huán)間隙長度=15 mm、間隙寬度=0.5 mm。

    (a) 口環(huán)間隙示意圖;(b) 離心泵流體域

    1.2 求解器和控制方程

    Fluent實(shí)際上是一個(gè)基于非結(jié)構(gòu)化網(wǎng)格的通用求解器,本文采用3D單精度,基于壓力的求解方法對計(jì)算模型進(jìn)行求解,并運(yùn)用多重旋轉(zhuǎn)坐標(biāo)系和有限體積法對非定常不可壓縮N-S 方程進(jìn)行離散,按順序依次求解動量方程、壓力修正方程、能量方程及其他標(biāo)量方程,如湍流方程等。在流場求解時(shí)所用連續(xù)方程、動量方程和能量方程可以寫成以下通用形式[13?14]:

    式中:為通用因變量;S為廣義源項(xiàng);為速度矢量;Γ為與相對應(yīng)的廣義擴(kuò)散系數(shù);為流體密度;為時(shí)間。

    離心泵內(nèi)部為復(fù)雜的三維湍流流動,故還需設(shè)置湍流模型,RNG?湍流模型[15]作為RANS方法中的一種形式,通過附加的湍動能、耗散率、比耗散率等湍流量進(jìn)行控制[16],考慮了湍流漩渦的影響,對強(qiáng)旋流動、脫流、漩渦等復(fù)雜流動的計(jì)算精度較高,更適用于對離心泵流場的數(shù)值模擬[17?18]。

    1.3 網(wǎng)格生成及邊界條件設(shè)置

    離心泵葉片、蝸殼等結(jié)構(gòu)扭曲程度較高,網(wǎng)格質(zhì)量將直接影響計(jì)算結(jié)果的準(zhǔn)確性,非結(jié)構(gòu)化網(wǎng)格對處理復(fù)雜問題具有很強(qiáng)的適應(yīng)性,更適用于復(fù)雜模型的數(shù)值模擬。因此,本文運(yùn)用ANSYS ICEM軟件,采用Tetra/Mixed非結(jié)構(gòu)化混合網(wǎng)格,對離心泵流體域進(jìn)行分塊網(wǎng)格劃分,并對口環(huán)間隙處進(jìn)行局部細(xì)化,使其內(nèi)部具有足夠的節(jié)點(diǎn),以保證計(jì)算精度的要求。網(wǎng)格劃分時(shí),首先將口環(huán)間隙處的最大網(wǎng)格尺寸設(shè)為0.5,口環(huán)間隙內(nèi)部的網(wǎng)格層數(shù)為2層;然后逐漸減小最大網(wǎng)格尺寸,當(dāng)口環(huán)間隙處的最大網(wǎng)格尺寸為0.2時(shí),口環(huán)間隙內(nèi)部的網(wǎng)格層數(shù)為4層,離心泵流體域的總體計(jì)算網(wǎng)格數(shù)量在120萬~220萬范圍內(nèi)發(fā)生變化。流體域網(wǎng)格及葉輪口環(huán)間隙處的網(wǎng)格如圖2所示。

    為了檢驗(yàn)網(wǎng)格劃分及口環(huán)間隙處網(wǎng)格層數(shù)對數(shù)值計(jì)算結(jié)果的影響,分別選定6組不同的網(wǎng)格數(shù)進(jìn)行數(shù)值模擬,分析網(wǎng)格的無關(guān)性。圖3所示為標(biāo)準(zhǔn)工況=50 m3/h,口環(huán)間隙寬度=0.5 mm時(shí),不同網(wǎng)格數(shù)量對泄漏量和揚(yáng)程影響的計(jì)算結(jié)果。由圖3可知,當(dāng)口環(huán)間隙內(nèi)部網(wǎng)格層數(shù)為4,網(wǎng)格總數(shù)達(dá)到200萬以后,隨著網(wǎng)格數(shù)的增加,設(shè)計(jì)工況點(diǎn)的泄漏量和揚(yáng)程變化較小,基本趨于穩(wěn)定。因此在數(shù)值計(jì)算中,離心泵模型網(wǎng)格數(shù)量控制在200萬上下波動,口環(huán)間隙內(nèi)部網(wǎng)格層數(shù)為4層。

    (a) 離心泵流體域網(wǎng)格;(b) 葉輪口環(huán)間隙網(wǎng)格

    1—H;2—q。

    離心泵流場計(jì)算域內(nèi)的流體為不可壓縮的牛頓流體,工作過程中流體的性質(zhì)不發(fā)生變化,并且不考慮重力對流場的影響。模型進(jìn)口為速度入口,出口為自由出流;近壁面采用標(biāo)準(zhǔn)壁面函數(shù),壁面取不滲透、無滑移固體壁面邊界條件,采用具有二階精度的QUICK格式離散,壓力修正為SIMPLC算法;計(jì)算過程中收斂精度設(shè)為1×10?4,外界條件為常溫常壓,介質(zhì)為水。在流量0分別為0.8,1.0和1.2等工況下,分別對離心泵的內(nèi)部流動狀態(tài)進(jìn)行定常流動數(shù)值模擬。

    2 口環(huán)間隙寬度對離心泵特性影響

    離心泵的口環(huán)、多級泵的級間軸封、填料密封、軸向力平衡裝置的動靜部件之間都有液體的泄漏,使離心泵輸出的液體體積減小并帶走了部分功率。為了驗(yàn)證數(shù)值模擬結(jié)果的準(zhǔn)確性,首先將不同口環(huán)間隙寬度情況下計(jì)算的模型水力效率與樣機(jī)試驗(yàn)值進(jìn)行對比,分析不同口環(huán)間隙寬度情況下離心泵的特性。其中樣機(jī)口環(huán)間隙寬度為0.5 mm,數(shù)值計(jì)算模型的口環(huán)間隙寬度分別為0.2,0.5和0.8 mm。離心泵水力效率的計(jì)算公式為

    h=(in?out)/(1 000)(2)

    式中:為實(shí)際流過葉輪的流量,m3/h;為葉輪扭矩,N·m;為葉輪轉(zhuǎn)動角速度,rad/s;in和out分別為離心泵的進(jìn)、出口壓力,Pa。

    離心泵水力效率試驗(yàn)值和模擬值的結(jié)果如圖4所示。由圖4可知:口環(huán)間隙的變化改變了離心泵的水力性能,不同口環(huán)間隙寬度下的模擬值和試驗(yàn)值之間均有不同程度的偏差,口環(huán)間隙寬度=0.2 mm時(shí),不同流量下的模擬值偏高;口環(huán)間隙寬度增大為0.8 mm時(shí),模擬值則低于試驗(yàn)值;而當(dāng)=0.5 mm時(shí),兩者具有較好的一致性。在標(biāo)準(zhǔn)工況=50 m3/h下,口環(huán)間隙寬度由0.2 mm增加到0.8 mm的過程中,水力效率由76%下降到69%,降低了近7%。

    1—試驗(yàn)值;2—b=0.2 mm;3—b=0.5 mm;4—b=0.8 mm。

    2.1 不同口環(huán)間隙寬度下離心泵壓力場分析

    圖5所示為標(biāo)準(zhǔn)工況下,不同口環(huán)間隙寬度在=0平面的靜壓分布云圖。由圖5(a)可知:口環(huán)間隙寬度改變了泵腔的壓力分布,其影響主要表現(xiàn)在前、后泵腔處;隨著口環(huán)間隙寬度的增大,離心泵前泵腔的低壓區(qū)向蝸室上移、擴(kuò)大,后泵腔的高壓區(qū)范圍縮小;由于泵腔內(nèi)軸向尺寸與徑向尺寸之比很小,故壓力分布只沿徑向發(fā)生明顯變化,且間隙寬度越小,前泵腔的壓力梯度越明顯。圖5(b)所示為前泵腔的靜壓分布曲線。由圖5(b)可知:前泵腔的壓力隨著半徑的減小,呈線性降低的趨勢,且隨著口環(huán)間隙寬度的增大,其斜率不斷增加;在3種不同間隙寬度的口環(huán)間隙中,泵腔進(jìn)口處的壓力相差不大,而在葉輪口環(huán)處,間隙寬度為0.2 mm時(shí),流體的靜壓力明顯高于其他間隙寬度下的壓力。前泵腔壓力的變化改變了流體的徑向流動,間隙內(nèi)部壓力越高,流體流入葉輪口環(huán)時(shí)的阻力就越大,葉輪口環(huán)間隙高壓區(qū)阻礙了流體向葉輪進(jìn)口的泄漏。

    (a) Y=0截面;(b) 靜壓分布曲線

    2.2 不同口環(huán)間隙下離心泵速度場分析

    在計(jì)算過程中,設(shè)Δ為流體在葉輪口環(huán)流道進(jìn)口壓力0與出口壓力l之差,流動方向上的長度為,間隙寬度為,則葉輪口環(huán)間隙流體平均速度m和壓差流動的阻力系數(shù)的計(jì)算公式分別為:

    (4)

    流體受到外部葉輪口環(huán)壁面剪切力的作用時(shí)發(fā)生變形流動,其內(nèi)部產(chǎn)生對抗變形的抵抗,并以內(nèi)摩擦力的形式表現(xiàn)出來。故葉輪口環(huán)間隙內(nèi)部的流體同時(shí)具有沿周向的流動(前泵蓋與葉輪的相對運(yùn)動產(chǎn)生)及軸向的壓力流動(前泵腔與葉輪進(jìn)口的壓差產(chǎn)生)。由式(3)和式(4)可知,離心泵口環(huán)間隙流體的流動速度對間隙寬度的變化更為敏感,當(dāng)口環(huán)間隙寬度減小時(shí),葉輪口環(huán)間隙流體平均速度m減小,阻力系數(shù)增大。圖6所示為不同口環(huán)間隙寬度下,葉輪口環(huán)處的平均速度分布云圖。由圖6可知:在泵腔與葉輪進(jìn)口壓差作用下的高速流體通過葉輪口環(huán)間隙流向泵的進(jìn)口,葉輪口環(huán)間隙寬度增大時(shí),間隙內(nèi)部的平均速度增加,當(dāng)口環(huán)間隙寬度=0.8 mm,葉輪口環(huán)內(nèi)部的平均速度約為13 m/s。因此,口環(huán)間隙寬度增大最終導(dǎo)致離心泵泄漏量的增加。

    b/mm:(a) 0.2;(b) 0.5;(c) 0.8

    3 口環(huán)結(jié)構(gòu)對離心泵特性影響

    為了減小泄漏量對離心泵性能的影響,保持口環(huán)間隙長度=15 mm不變,將口環(huán)間隙寬度增加到1.5 mm,并在葉輪口環(huán)上布置圓形、矩形和鋸齒形環(huán)狀結(jié)構(gòu),環(huán)狀結(jié)構(gòu)頂端距前泵蓋的間隙寬度為0.5 mm(與原模型口環(huán)間隙寬度一致)。離心泵不同口環(huán)結(jié)構(gòu)型式如圖7所示,圖中箭頭為液體的流動方向,各類型環(huán)狀結(jié)構(gòu)的軸向間距為=3 mm,高度=1 mm。

    3.1 泄漏量和容積效率分析

    分析不同口環(huán)結(jié)構(gòu)下離心泵性能的變化,取0=0.8,1.0,1.23個(gè)流量點(diǎn),對間隙泄漏量、及容積效率進(jìn)行計(jì)算,口環(huán)間隙泄漏量的基本計(jì)算公式為

    式中:a為間隙泄漏的過流面積,m3;Δc1為液體流過密封間隙的損失水頭,m;為間隙損失系數(shù);為重力加速度。

    則容積效率表示為

    v=/(+) (6)

    圖8所示為不同結(jié)構(gòu)口環(huán)下的間隙泄漏量曲線。由圖8可知:在不同流量下,離心泵泄漏量曲線變化趨勢基本一致,均隨著的增加,呈先增加后減小的變化趨勢;但在小于設(shè)計(jì)流量點(diǎn)時(shí),泄漏量的變化趨勢并不明顯;圓形環(huán)狀葉輪口環(huán)在3種工況下的泄漏量與其他形式的葉輪口環(huán)相比偏差較大,這可能與其流道內(nèi)的流動狀態(tài)有關(guān),導(dǎo)致泄漏量增加,在標(biāo)準(zhǔn)工況下泄露量增加約30%。計(jì)算結(jié)果說明采用圓形環(huán)狀葉輪口環(huán)時(shí)并沒起到減小泄漏損失的作用;矩形和鋸齒形的葉輪口環(huán)與原模型相比則具有較好的減小泄漏的效果,鋸齒形結(jié)構(gòu)的葉輪口環(huán)泄漏量最多可減少16.2%。

    (a) 原模型;(b) 矩形環(huán)狀結(jié)構(gòu)口環(huán);(c) 圓形環(huán)狀結(jié)構(gòu)口環(huán);(d) 鋸齒形環(huán)狀結(jié)構(gòu)口環(huán)

    1—原模型;2—圓形密封齒;3—矩形密封齒;4—鋸齒形密封齒。

    圖9所示為不同葉輪口環(huán)結(jié)構(gòu)的容積效率曲線。由圖9可知:圓形葉輪口環(huán)在各工況下的容積效率均低于其他3種口環(huán)結(jié)構(gòu);鋸齒形和矩形環(huán)狀結(jié)構(gòu)葉輪口環(huán)則均能有效提升離心泵的容積效率,且鋸齒形葉輪口環(huán)的提升幅度更優(yōu)于矩形結(jié)構(gòu);在不同工況下采用鋸齒形葉輪口環(huán)時(shí),離心泵的容積效率最高可提升3%。由此可知,口環(huán)間隙寬度雖然增加到了1.5 mm,但布置在葉輪口環(huán)上的鋸齒形和矩形環(huán)狀結(jié)構(gòu)也能減小泄漏,提高容積效率。

    1—原模型;2—圓形密封齒;3—矩形密封齒;4—鋸齒形密封齒。

    3.2 軸向力分析

    圖10所示為不同口環(huán)結(jié)構(gòu)模型的軸向力曲線。由圖10可知:不同工況下,4種口環(huán)結(jié)構(gòu)下軸向力的變化趨勢基本一致;對于同一葉輪口環(huán)結(jié)構(gòu)的離心泵,軸向力均隨著流量的增加而減小;同一流量點(diǎn)下,圓形葉輪口環(huán)離心泵的軸向力最大,矩形和鋸齒形葉輪口環(huán)離心泵的軸向力均小于原模型,矩形環(huán)狀結(jié)構(gòu)口環(huán)平衡軸向力的效果最佳,軸向力最大減小10%。

    1—原模型;2—圓形密封齒;3—矩形密封齒;4—鋸齒形密封齒。

    圖11所示為不同葉輪口環(huán)結(jié)構(gòu)的離心泵在標(biāo)準(zhǔn)工況1.0下,=0截面的靜壓分布云圖。由圖11可知:鋸齒形和矩形口環(huán)結(jié)構(gòu)的葉輪前泵腔壓力分布與口環(huán)間隙寬度為0.5 mm的原模型相比,低壓區(qū)范圍擴(kuò)大且后泵腔的高壓區(qū)范圍縮小,而圓形環(huán)狀結(jié)構(gòu)的葉輪口環(huán)內(nèi)部的壓力分布并沒有明顯變化。圓形結(jié)構(gòu)的葉輪口環(huán),低壓區(qū)的范圍進(jìn)一步縮小且后泵腔的高壓區(qū)消失,故其軸向力不但沒有減小反而增加。離心泵葉輪上的軸向力與前、后泵腔壓力的變化密切相關(guān),故鋸齒形和矩形口環(huán)結(jié)構(gòu)形式的密封口環(huán)可以起到平衡離心泵葉輪的軸向力的作用;同時(shí)口環(huán)間隙內(nèi)部出現(xiàn)了較大的壓力梯度,壓力梯度的變化使流體的湍流脈動和雷諾應(yīng)力增加,鋸齒形和矩形口環(huán)結(jié)構(gòu)改變了間隙內(nèi)部流體的流動狀態(tài),促進(jìn)了流體由層流流動向湍流流動的轉(zhuǎn)變,流體的湍動能增加,實(shí)現(xiàn)了對流體的阻礙。

    3.3 速度場分析

    圖12所示為不同葉輪口環(huán)結(jié)構(gòu)下口環(huán)間隙的速度場分布云圖。由圖12(a)可知:在改變?nèi)~輪口環(huán)的結(jié)構(gòu)后,兩環(huán)狀結(jié)構(gòu)之間的空腔內(nèi)存在明顯的低速區(qū),且從前泵腔到葉輪進(jìn)口低速區(qū)的范圍逐漸增大;速度從前泵腔到葉輪進(jìn)口逐漸降低,與原模型口環(huán)相比,速度梯度的變化較為明顯;一個(gè)環(huán)狀結(jié)構(gòu)和一個(gè)膨脹空腔構(gòu)成了一級迷宮密封,當(dāng)高速流動的流體遇到環(huán)狀結(jié)構(gòu)的阻礙后進(jìn)入節(jié)流齒隙,導(dǎo)致流體的湍動狀態(tài)加劇;壓力推動流體由高壓側(cè)流向低壓側(cè),流體流過節(jié)流齒隙,由于流道變窄,流速增高,壓力降低,壓力能轉(zhuǎn)變?yōu)閯幽?。以矩形結(jié)構(gòu)葉輪口環(huán)間隙為例,流動狀態(tài)如12(b)所示,流體以高速進(jìn)入兩環(huán)狀結(jié)構(gòu)之間的空腔時(shí),過流面積突然增加,流速降低,并產(chǎn)生劇烈旋渦,有效阻止了流體壓力能的恢復(fù),最大程度地分散經(jīng)密封齒隙高速射出的流體的動能;在壓力梯度作用下,流體將會產(chǎn)生回流形成邊界層分離;從下游回流的液體受到主流的沖擊折返流向下游,在分離點(diǎn)形成漩渦和尾流區(qū),尾流區(qū)的漩渦使流體的阻力劇增,減小葉輪口環(huán)間隙的泄漏。

    (a) Y=0截面;(b) 口環(huán)間隙處

    (a) 口環(huán)間隙處速度云圖;(b) 口環(huán)間隙處速度矢量圖

    4 結(jié)論

    1) 口環(huán)間隙寬度的變化改變了離心泵前后泵腔的壓力分布,口環(huán)間隙增大,前泵腔的低壓區(qū)向蝸室上移、擴(kuò)大,后泵腔的高壓區(qū)范圍縮?。婚g隙內(nèi)部流體的速度隨間隙的增大而增加,容積效率下降。

    2) 口環(huán)間隙寬度增大到1.5 mm的情況下,鋸齒形和矩形環(huán)狀結(jié)構(gòu)葉輪口環(huán)均能明顯減小泄漏量,最多減少16.2%,不同工況下采用鋸齒形葉輪口環(huán)時(shí),離心泵的容積效率最高可提升3%。

    3) 在不同結(jié)構(gòu)的口環(huán)中,圓形葉輪口環(huán)離心泵的軸向力最大,矩形和鋸齒形葉輪口環(huán)離心泵的軸向力均小于原模型,矩形環(huán)狀結(jié)構(gòu)口環(huán)平衡軸向力的效果最佳,軸向力最大減小10%。

    4) 環(huán)狀結(jié)構(gòu)葉輪口環(huán)改變了離心泵口環(huán)間隙內(nèi)部流體的流動狀態(tài),矩形和鋸齒形狀結(jié)構(gòu)之間的空腔內(nèi)存在明顯的低速區(qū),間隙內(nèi)部速度梯度、壓力梯度增加,實(shí)現(xiàn)了增加流體流動阻力的作用,最終導(dǎo)致離心泵口環(huán)處泄漏量減小,容積效率提高。

    [1] 周凌九, 王占民, 江東智, 等. 離心泵非定常流動計(jì)算及性能預(yù)測[J]. 排灌機(jī)械工程學(xué)報(bào), 2010, 28(4): 286?290. ZHOU Lingjiu, WANG Zhanmin, JIANG Dongzhi, et al. Computation of unsteady flow and performance prediction for centrifugal pump[J]. Journal of Drainage and Irrigation Machinery Engineering, 2010, 28(4): 286?290.

    [2] 吳大轉(zhuǎn), 許斌杰, 武鵬, 等. 多級離心泵內(nèi)部間隙流動與泄漏損失[J]. 浙江大學(xué)學(xué)報(bào), 2011, 45(8): 1393?1398. WU Dazhuan, XU Binjie, WU Peng, et al. Internal clearance flow and leakage loss in multistage centrifugal pump[J]. Journal of Zhejiang University, 2011, 45(8): 1393?1398.

    [3] 馬文生, 陳照波, 焦映厚, 等. 迷宮密封結(jié)構(gòu)對泄漏量和軸系臨界轉(zhuǎn)速影響分析研究[J]. 振動工程學(xué)報(bào), 2013, 26(6): 823?830. MA Wensheng, CHEN Zhaobo, JIAO Yinghou, et al. Leakage and critical speed effect of labyrinth seal structure[J]. Journal of Vibration Engineering, 2013, 26(6): 823?830.

    [4] 王煒哲, 劉應(yīng)征, 葉春, 等. 迷宮密封-轉(zhuǎn)子系統(tǒng)動力學(xué)特性的試驗(yàn)測量和數(shù)值模擬[J]. 機(jī)械工程學(xué)報(bào), 2007, 43(3): 22?27. WANG Weizhe, LIU Yingzheng, Ye Chun, et al. Experimental measurement and numerical simulation of dynamics of labyrinth seal-rotor combination[J]. Journal of Mechanical Engineering, 2007, 43(3): 22?27.

    [5] LUCIUS A, BRENNER G. Unsteady CFD simulation of a pump in part load conditions using scale-adaptive simulation[J]. International Journal of Heat and Fluid Flow, 2010, 31(6): 1113?1118.

    [6] BENIGNI H, JABERG H, YEUNG H, et al. Numerical simulation of low specific speed american petroleum institute pumps in part-load operation and comparison with test rig results[J]. ASME Journal of Fluids Engineering, 2012, 134(2): 024501.

    [7] 李曉俊, 袁壽其, 潘中永, 等. 基于結(jié)構(gòu)化網(wǎng)格的離心泵全流場數(shù)值模擬[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2013, 44(7): 50?54. LI Xiaojun, YUAN Shouqi, PAN Zhongyong, et al. Numerical simulation of whole flow field for centrifugal pump with structured grid[J]. Transactions of the Chinese Society for Agricultural Machinery, 2013, 44(7): 50?54.

    [8] 譚磊, 曹樹良, 桂紹波, 等. 帶有前置導(dǎo)葉離心泵空化性能的試驗(yàn)及數(shù)值模擬[J]. 機(jī)械工程學(xué)報(bào), 2010, 46(18): 177?182. TAN Lei, CAO Shuliang, GUI Shaobo, et al. Experiment and numerical simulation of cavitation performance for centrifugal pump with inlet guide vane[J]. Journal of Mechanical Engineering, 2010, 46(18): 177?182.

    [9] LUCIUS A, BRENNER G. Numerical simulation and evaluation of velocity fluctuations during rotating stall of a centrifugal pump[J]. Journal of Fluids Engineering, Transactions of the ASME, 2011, 133(8): 1?8.

    [10] ZHANG Y L, Li Y, CUI B L, et al. Numerical simulation and analysis of solid-liquid two-phase flow in centrifugal pump[J]. Chinese Journal of Mechanical Engineering, 2013, 26(1): 53?60.

    [11] 周嶺, 施衛(wèi)東, 陸偉剛, 等. 深井離心泵數(shù)值模擬與試驗(yàn)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2011, 42(3): 69?73. ZHOU Ling, SHI Weidong, LU Weigang, et al. Numerical simulation and experiment on deep-well centrifugal pump[J]. Transactions of the Chinese Society for Agricultural Machinery, 2011, 42(3): 69?73.

    [12] 袁壽其, 徐宇平, 張金鳳, 等. 流固耦合作用對螺旋離心泵流場影響的數(shù)值分析[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào), 2013, 44(1): 38?42. YUAN Shouqi, XU Yuping, ZHANG Jinfeng, et al. Numerical analysis for effect of fluid-structure interaction on flow field in screw centrifugal pump[J]. Transactions of the Chinese Society for Agricultural Machinery, 2013, 44(1): 38?42.

    [13] 谷云慶, 牟介剛, 代東順, 等. 基于蚯蚓背孔射流的仿生射流表面減阻性能研究[J]. 物理學(xué)報(bào), 2015, 64(2): 310?319. GU Yunqing, MOU Jiegang, DAI Dongshun, et al. Characteristics on drag reduction of bionic jet surface based on earthworm’s back orifice jet[J]. Physica Sinica, 2015, 64(2): 310?319.

    [14] GU Yunqing, ZHAO Gang, ZHENG Jinxing, et al. Experimental and numerical investigation on drag reduction of non-smooth bionic jet surface[J]. Ocean Engineering, 2014, 81: 50?57.

    [15] MOUTAOUAKIL L E, ZRIKEM Z, ABDELBAKI A. Performance of various RANS eddy-viscosity models for turbulent natural convection in tall vertical cavities[J]. Heat and Mass Transfer, 2014, 50(8): 1103?1113.

    [16] SONG X G, Park J, KIM S G, et al. Performance comparison and erosion prediction of jet pumps by using a numerical method[J]. Mathematical and Computer Modeling, 2013, 57(1): 245?253.

    [17] DING H, VISSER F C, JIANG Y, et al. Demonstration and validation of a 3D CFD simulation tool predicting pump performance and cavitation for industrial applications[J]. ASME Journal of Fluids Engineering, 2011, 133(1): 1?14.

    [18] 林琳, 王言英. 不同湍流模型下圓柱渦激振動的計(jì)算比較[J]. 船舶力學(xué), 2013, 17(10): 1115?1125. LIN Lin, WANG Yanying. Comparison between different turbulence models on vortex induced vibration of circular cylinder[J]. Journal of Ship Mechanics, 2013, 17(10): 1115?1125.

    (編輯 趙俊)

    Influences of impeller ring structure on performance and flow field of centrifugal pump

    MOU Jiegang1, DAI Dongshun1, GU Yunqing1, LIU Jian1, WU Denghao2, MA Yi1

    (1. College of Mechanical Engineering, Zhejiang University of Technology, Hangzhou 310014, China;2. Zhijiang College, Zhejiang University of Technology, Hangzhou 310024, China)

    In order to study the impeller ring clearance and structural change effects on the performance of centrifugal pump, based on IS80-65-160 centrifugal pump, internal flow field for different impeller ring clearances and structures were numerically simulated through RNG-turbulence model. The influences on efficiency, pressure field and velocity field were analyzed, and three kinds of impeller ring structures effect on the centrifugal pump were researched. The results show that impeller ring clearances changed the internal flow field of centrifugal pump and affected the pressure distribution of pump cavity. Among the three kinds of annularity structures, rectangle and serrated impeller ring can reduce leakage, (at most reduced by 16.2%), improve volumetric efficiency, balance action of axial thrust. Hydraulic efficiency increases by 2% after the change of impeller ring structure. When the fluid flows through the cavity between two rings, volume expansion and energy dissipation violently can prevent fluid restore its pressure, maximum dispersion the kinetic energy of high-speed injection fluid sealing from sealed gear. The internal fluid velocity gradient and pressure gradient are enlarged, and turbulent kinetic energy increases, and so the leakage is reduced.

    centrifugal pump; hydrodynamics; impeller ring; type of construction; numerical simulation

    10.11817/j.issn.1672?7207.2017.06.015

    TB311

    A

    1672?7207(2017)06?1522?08

    2016?06?30;

    2016?09?07

    國家自然科學(xué)基金資助項(xiàng)目(51406183);浙江省自然科學(xué)基金資助項(xiàng)目(LQ15E050005);中國博士后科學(xué)基金資助項(xiàng)目(2016M601736)(Project(51406183) supported by National Natural Science Foundation of China; Project(LQ15E050005) supported by Zhejiang Provincial Natural Science Foundation of China; Project(2016M601736) supported by the China Postdoctoral Science Foundation)

    谷云慶,博士,碩士生導(dǎo)師,從事離心泵流場理論研究;E-mail:guyunqing@hrbeu.edu.cn

    猜你喜歡
    結(jié)構(gòu)
    DNA結(jié)構(gòu)的發(fā)現(xiàn)
    《形而上學(xué)》△卷的結(jié)構(gòu)和位置
    論結(jié)構(gòu)
    中華詩詞(2019年7期)2019-11-25 01:43:04
    新型平衡塊結(jié)構(gòu)的應(yīng)用
    模具制造(2019年3期)2019-06-06 02:10:54
    循環(huán)結(jié)構(gòu)謹(jǐn)防“死循環(huán)”
    論《日出》的結(jié)構(gòu)
    縱向結(jié)構(gòu)
    縱向結(jié)構(gòu)
    我國社會結(jié)構(gòu)的重建
    人間(2015年21期)2015-03-11 15:23:21
    創(chuàng)新治理結(jié)構(gòu)促進(jìn)中小企業(yè)持續(xù)成長
    国产精品野战在线观看| 狠狠狠狠99中文字幕| 国产淫片久久久久久久久| 亚洲精品乱码久久久v下载方式| netflix在线观看网站| 两个人的视频大全免费| 亚洲内射少妇av| 97热精品久久久久久| 精品午夜福利视频在线观看一区| 日韩,欧美,国产一区二区三区 | 亚洲av不卡在线观看| 欧美日本视频| 国产真实伦视频高清在线观看 | 国语自产精品视频在线第100页| 国产真实伦视频高清在线观看 | 在线a可以看的网站| 少妇丰满av| 嫩草影院精品99| 亚洲人与动物交配视频| 国产精品一区二区免费欧美| 三级国产精品欧美在线观看| 国产成年人精品一区二区| 精品久久久久久久末码| 欧美高清性xxxxhd video| 亚洲欧美日韩东京热| 午夜福利18| 黄色视频,在线免费观看| 51国产日韩欧美| 精品久久久久久久久久久久久| 国产精品国产高清国产av| 成人三级黄色视频| 亚洲aⅴ乱码一区二区在线播放| 蜜桃亚洲精品一区二区三区| 国产毛片a区久久久久| 国产 一区精品| 国内毛片毛片毛片毛片毛片| 亚洲av不卡在线观看| 国内精品一区二区在线观看| 国产高清激情床上av| 岛国在线免费视频观看| 欧美极品一区二区三区四区| 老熟妇乱子伦视频在线观看| 一进一出好大好爽视频| 免费观看的影片在线观看| 乱码一卡2卡4卡精品| 国产视频一区二区在线看| 少妇丰满av| 亚洲国产精品成人综合色| 精品日产1卡2卡| 日韩,欧美,国产一区二区三区 | 精品久久久噜噜| 中文字幕人妻熟人妻熟丝袜美| 成年版毛片免费区| 我的女老师完整版在线观看| 国产欧美日韩一区二区精品| 成人一区二区视频在线观看| 亚洲精品一区av在线观看| 国国产精品蜜臀av免费| 国产69精品久久久久777片| 18+在线观看网站| 男人舔奶头视频| a级毛片免费高清观看在线播放| 精品人妻熟女av久视频| 在线观看舔阴道视频| av中文乱码字幕在线| 嫩草影院精品99| 麻豆精品久久久久久蜜桃| 国产精品人妻久久久久久| 日日啪夜夜撸| 俄罗斯特黄特色一大片| 亚洲av电影不卡..在线观看| 久久香蕉精品热| 亚洲一级一片aⅴ在线观看| 午夜福利视频1000在线观看| 国内精品久久久久精免费| 日本黄大片高清| 波多野结衣高清作品| 国产欧美日韩一区二区精品| 一本精品99久久精品77| 全区人妻精品视频| 色哟哟·www| 亚洲成人久久爱视频| 小蜜桃在线观看免费完整版高清| av国产免费在线观看| 岛国在线免费视频观看| 亚洲人成网站高清观看| 国产男靠女视频免费网站| 国产精品自产拍在线观看55亚洲| 成熟少妇高潮喷水视频| 国产高清不卡午夜福利| 精品久久久久久久久av| 性欧美人与动物交配| 欧美绝顶高潮抽搐喷水| 国产私拍福利视频在线观看| 日韩 亚洲 欧美在线| 色综合婷婷激情| 99国产极品粉嫩在线观看| 亚洲狠狠婷婷综合久久图片| 欧美3d第一页| 成人国产综合亚洲| 欧美黑人欧美精品刺激| 麻豆精品久久久久久蜜桃| 欧美绝顶高潮抽搐喷水| 99热只有精品国产| 午夜激情欧美在线| 内射极品少妇av片p| www日本黄色视频网| 女人十人毛片免费观看3o分钟| 女的被弄到高潮叫床怎么办 | 别揉我奶头~嗯~啊~动态视频| 真人做人爱边吃奶动态| 日韩中字成人| 亚洲av免费高清在线观看| 成人一区二区视频在线观看| 国产三级在线视频| 99riav亚洲国产免费| 精品久久久久久久久久免费视频| 午夜免费成人在线视频| 69av精品久久久久久| av黄色大香蕉| 精品人妻偷拍中文字幕| 搡老妇女老女人老熟妇| 亚洲天堂国产精品一区在线| av中文乱码字幕在线| 狠狠狠狠99中文字幕| 尤物成人国产欧美一区二区三区| 日韩一本色道免费dvd| 极品教师在线免费播放| 亚洲美女搞黄在线观看 | 日本精品一区二区三区蜜桃| 两人在一起打扑克的视频| 人妻制服诱惑在线中文字幕| 久久这里只有精品中国| 能在线免费观看的黄片| 日韩精品中文字幕看吧| 欧美精品国产亚洲| 99久国产av精品| 在线a可以看的网站| 天堂√8在线中文| aaaaa片日本免费| 看免费成人av毛片| АⅤ资源中文在线天堂| 特大巨黑吊av在线直播| av在线蜜桃| 露出奶头的视频| 中出人妻视频一区二区| 色哟哟哟哟哟哟| 亚洲狠狠婷婷综合久久图片| 国产大屁股一区二区在线视频| 两个人的视频大全免费| 天美传媒精品一区二区| 天堂网av新在线| 99国产极品粉嫩在线观看| 欧美在线一区亚洲| 国产伦精品一区二区三区四那| 自拍偷自拍亚洲精品老妇| 国产精品久久久久久av不卡| 免费人成在线观看视频色| 1000部很黄的大片| 国产精品av视频在线免费观看| 免费观看的影片在线观看| 午夜日韩欧美国产| 国产免费av片在线观看野外av| 一本一本综合久久| 成人国产综合亚洲| 天堂√8在线中文| 免费看av在线观看网站| 熟女电影av网| 日韩欧美免费精品| 国产综合懂色| 97碰自拍视频| 中亚洲国语对白在线视频| 精品一区二区三区av网在线观看| 亚洲欧美日韩无卡精品| 一本一本综合久久| 亚洲av五月六月丁香网| 亚洲无线在线观看| 精品人妻一区二区三区麻豆 | 婷婷色综合大香蕉| 午夜爱爱视频在线播放| 欧美激情久久久久久爽电影| 少妇的逼好多水| 男人狂女人下面高潮的视频| 12—13女人毛片做爰片一| 校园人妻丝袜中文字幕| 99热精品在线国产| 99热网站在线观看| 干丝袜人妻中文字幕| 女同久久另类99精品国产91| 免费看美女性在线毛片视频| 国产大屁股一区二区在线视频| 国产高潮美女av| www日本黄色视频网| 成人永久免费在线观看视频| 亚洲成人精品中文字幕电影| 18禁在线播放成人免费| 热99re8久久精品国产| 欧美高清成人免费视频www| 欧美+亚洲+日韩+国产| 久久久精品欧美日韩精品| 成人av在线播放网站| 联通29元200g的流量卡| 不卡视频在线观看欧美| 成人国产一区最新在线观看| 亚洲一级一片aⅴ在线观看| videossex国产| 尤物成人国产欧美一区二区三区| 国产乱人伦免费视频| 国产不卡一卡二| 亚洲欧美日韩卡通动漫| 亚洲精品456在线播放app | 免费av毛片视频| 亚洲精品久久国产高清桃花| 国产精品永久免费网站| 波野结衣二区三区在线| 草草在线视频免费看| 欧美黑人巨大hd| 国产在视频线在精品| 日本黄色视频三级网站网址| 国产一区二区亚洲精品在线观看| 国内精品久久久久久久电影| 看片在线看免费视频| av中文乱码字幕在线| 午夜免费激情av| 尾随美女入室| 床上黄色一级片| 色av中文字幕| 日韩人妻高清精品专区| 亚洲五月天丁香| 人人妻人人看人人澡| 色综合亚洲欧美另类图片| 亚洲人成伊人成综合网2020| 一区二区三区激情视频| 我要搜黄色片| 午夜精品一区二区三区免费看| 日本与韩国留学比较| 俺也久久电影网| 亚洲无线观看免费| 久久中文看片网| 麻豆精品久久久久久蜜桃| 三级男女做爰猛烈吃奶摸视频| 琪琪午夜伦伦电影理论片6080| 午夜激情福利司机影院| 久久久久久久久久久丰满 | 国产欧美日韩一区二区精品| 成人一区二区视频在线观看| 日韩一本色道免费dvd| 精品免费久久久久久久清纯| 一本精品99久久精品77| 亚洲一区二区三区色噜噜| 床上黄色一级片| 国产精品一区www在线观看 | 国产精品无大码| 欧美成人a在线观看| 黄色欧美视频在线观看| 麻豆久久精品国产亚洲av| 日韩强制内射视频| 国产v大片淫在线免费观看| 国产精品亚洲一级av第二区| 免费高清视频大片| 一级av片app| 少妇高潮的动态图| 免费av毛片视频| 国产亚洲av嫩草精品影院| 日韩一区二区视频免费看| 婷婷精品国产亚洲av在线| 69人妻影院| 男人和女人高潮做爰伦理| 啦啦啦观看免费观看视频高清| 一级黄色大片毛片| avwww免费| 97热精品久久久久久| 99热这里只有是精品50| 韩国av一区二区三区四区| 亚洲一区二区三区色噜噜| 哪里可以看免费的av片| 啦啦啦啦在线视频资源| 91狼人影院| 亚洲精品一区av在线观看| 18禁在线播放成人免费| 最近在线观看免费完整版| 99热这里只有精品一区| 伦理电影大哥的女人| 午夜a级毛片| 久久久久国内视频| 精品久久久久久久久久免费视频| 麻豆久久精品国产亚洲av| 成人特级av手机在线观看| 熟女人妻精品中文字幕| 国产亚洲精品av在线| 美女免费视频网站| 狠狠狠狠99中文字幕| 精品久久久久久久久久久久久| av.在线天堂| 久久精品影院6| 国产麻豆成人av免费视频| 欧美绝顶高潮抽搐喷水| 久久精品久久久久久噜噜老黄 | 午夜福利成人在线免费观看| 久久欧美精品欧美久久欧美| 真人一进一出gif抽搐免费| 看十八女毛片水多多多| 在线观看美女被高潮喷水网站| 啦啦啦观看免费观看视频高清| 成人国产麻豆网| 国产麻豆成人av免费视频| 在线观看免费视频日本深夜| 亚洲国产精品成人综合色| 欧美高清成人免费视频www| 我的女老师完整版在线观看| 高清在线国产一区| 一级a爱片免费观看的视频| 成人性生交大片免费视频hd| 国产精品久久久久久久电影| 亚洲五月天丁香| 不卡视频在线观看欧美| a级毛片a级免费在线| 国内精品久久久久精免费| 亚洲av.av天堂| 精品99又大又爽又粗少妇毛片 | 国产久久久一区二区三区| 婷婷色综合大香蕉| 91麻豆av在线| 一a级毛片在线观看| 久久草成人影院| 日本与韩国留学比较| 99热6这里只有精品| 12—13女人毛片做爰片一| 国内精品宾馆在线| 1000部很黄的大片| 白带黄色成豆腐渣| 日本a在线网址| netflix在线观看网站| 国国产精品蜜臀av免费| 少妇的逼水好多| 九色成人免费人妻av| 国产色婷婷99| 欧美黑人欧美精品刺激| 天堂网av新在线| 亚州av有码| 老师上课跳d突然被开到最大视频| 成人国产麻豆网| 一卡2卡三卡四卡精品乱码亚洲| 亚洲av中文av极速乱 | 美女免费视频网站| 欧美极品一区二区三区四区| 精华霜和精华液先用哪个| 成人二区视频| 精品久久久久久,| 在线看三级毛片| 欧美又色又爽又黄视频| 午夜激情欧美在线| 国产三级在线视频| 亚洲精品456在线播放app | 中亚洲国语对白在线视频| 亚洲狠狠婷婷综合久久图片| 国产精品野战在线观看| 国产一区二区激情短视频| 人妻丰满熟妇av一区二区三区| 国产91精品成人一区二区三区| 精品久久久久久久久久久久久| 男女做爰动态图高潮gif福利片| 国产v大片淫在线免费观看| 久久国产乱子免费精品| 美女 人体艺术 gogo| 久久精品人妻少妇| 亚洲成a人片在线一区二区| 亚洲精品在线观看二区| 99热这里只有是精品在线观看| av中文乱码字幕在线| 亚洲成a人片在线一区二区| 午夜免费成人在线视频| 久久草成人影院| av国产免费在线观看| 精品日产1卡2卡| 色噜噜av男人的天堂激情| 欧美日韩亚洲国产一区二区在线观看| 久久精品夜夜夜夜夜久久蜜豆| 免费黄网站久久成人精品| 国产高清激情床上av| 欧美国产日韩亚洲一区| 国产精品免费一区二区三区在线| 天堂网av新在线| 久久人人爽人人爽人人片va| 天堂av国产一区二区熟女人妻| 男女边吃奶边做爰视频| 在线国产一区二区在线| 国产伦人伦偷精品视频| 网址你懂的国产日韩在线| 少妇猛男粗大的猛烈进出视频 | 国产伦在线观看视频一区| 国内毛片毛片毛片毛片毛片| 看免费成人av毛片| 国模一区二区三区四区视频| 身体一侧抽搐| 午夜激情欧美在线| 国产 一区精品| 欧美极品一区二区三区四区| 99热网站在线观看| 免费观看精品视频网站| 亚洲av日韩精品久久久久久密| 成人精品一区二区免费| 日韩欧美免费精品| 国产乱人伦免费视频| 午夜老司机福利剧场| 亚洲成人免费电影在线观看| 别揉我奶头~嗯~啊~动态视频| 男女啪啪激烈高潮av片| 国产单亲对白刺激| 久久精品国产亚洲av香蕉五月| 午夜久久久久精精品| 美女高潮喷水抽搐中文字幕| 99久久精品一区二区三区| 九九热线精品视视频播放| 韩国av一区二区三区四区| 搡老妇女老女人老熟妇| 在线a可以看的网站| 国产亚洲av嫩草精品影院| 免费在线观看影片大全网站| 在线看三级毛片| 日本色播在线视频| 欧美国产日韩亚洲一区| 久久人人爽人人爽人人片va| 在线观看舔阴道视频| 最近中文字幕高清免费大全6 | 美女xxoo啪啪120秒动态图| 亚洲va日本ⅴa欧美va伊人久久| 久久6这里有精品| 国语自产精品视频在线第100页| 亚洲va在线va天堂va国产| 国产午夜福利久久久久久| 乱系列少妇在线播放| 自拍偷自拍亚洲精品老妇| 老司机福利观看| 精品一区二区免费观看| 99精品在免费线老司机午夜| 极品教师在线免费播放| h日本视频在线播放| 日本一二三区视频观看| 日韩一本色道免费dvd| 亚洲一区高清亚洲精品| 久久99热这里只有精品18| 国产高清三级在线| 露出奶头的视频| 天天一区二区日本电影三级| 亚洲va在线va天堂va国产| 日韩欧美 国产精品| 亚洲欧美激情综合另类| 久久亚洲精品不卡| 午夜精品一区二区三区免费看| a级毛片免费高清观看在线播放| 熟妇人妻久久中文字幕3abv| videossex国产| 我的老师免费观看完整版| aaaaa片日本免费| 国产不卡一卡二| 真实男女啪啪啪动态图| 露出奶头的视频| 日韩,欧美,国产一区二区三区 | 91精品国产九色| 久久久久久国产a免费观看| 国产av一区在线观看免费| 日本 av在线| 免费黄网站久久成人精品| eeuss影院久久| 精品日产1卡2卡| 国产 一区精品| 久久精品国产99精品国产亚洲性色| 人人妻,人人澡人人爽秒播| 最近最新免费中文字幕在线| 日韩欧美在线乱码| 国产私拍福利视频在线观看| 国产成人av教育| 亚洲第一区二区三区不卡| 亚洲最大成人av| 男女那种视频在线观看| 国产精品久久久久久精品电影| x7x7x7水蜜桃| 国产一区二区三区视频了| 偷拍熟女少妇极品色| 可以在线观看毛片的网站| 一区福利在线观看| 最新在线观看一区二区三区| 麻豆久久精品国产亚洲av| 嫩草影院入口| 国产黄片美女视频| 舔av片在线| 国产亚洲精品综合一区在线观看| 亚洲精品国产成人久久av| 色综合婷婷激情| 国产精品自产拍在线观看55亚洲| 在线播放无遮挡| 最近视频中文字幕2019在线8| 国产伦精品一区二区三区四那| 国产欧美日韩精品亚洲av| 精品免费久久久久久久清纯| 日韩大尺度精品在线看网址| 亚洲中文字幕日韩| 男人的好看免费观看在线视频| 三级男女做爰猛烈吃奶摸视频| 免费人成视频x8x8入口观看| 黄色日韩在线| 国产成人影院久久av| 国产午夜福利久久久久久| 精品人妻熟女av久视频| 亚洲精品亚洲一区二区| 美女 人体艺术 gogo| 日韩欧美精品v在线| 女的被弄到高潮叫床怎么办 | 午夜免费激情av| 日韩欧美 国产精品| 少妇人妻精品综合一区二区 | 91精品国产九色| 久久精品夜夜夜夜夜久久蜜豆| 久久精品国产鲁丝片午夜精品 | 99精品在免费线老司机午夜| 丰满乱子伦码专区| 91狼人影院| 欧美zozozo另类| 欧美日韩综合久久久久久 | 国产精品亚洲美女久久久| 最近中文字幕高清免费大全6 | 国产av不卡久久| 亚洲成av人片在线播放无| 少妇裸体淫交视频免费看高清| 黄色配什么色好看| 日本 欧美在线| 999久久久精品免费观看国产| 精品欧美国产一区二区三| 成人国产一区最新在线观看| 97超级碰碰碰精品色视频在线观看| 美女大奶头视频| 亚洲精品乱码久久久v下载方式| 99久久中文字幕三级久久日本| 大又大粗又爽又黄少妇毛片口| 露出奶头的视频| 日日摸夜夜添夜夜添av毛片 | 国产探花极品一区二区| 国产精品久久久久久av不卡| 久久精品91蜜桃| 高清毛片免费观看视频网站| 国产精品1区2区在线观看.| 色精品久久人妻99蜜桃| 日韩国内少妇激情av| 国产精品久久久久久av不卡| 国产老妇女一区| 午夜福利视频1000在线观看| 久久热精品热| 99热6这里只有精品| 日本五十路高清| 成年免费大片在线观看| 日日啪夜夜撸| 国产69精品久久久久777片| 美女免费视频网站| 九九爱精品视频在线观看| 免费电影在线观看免费观看| 久久久精品大字幕| 成年女人永久免费观看视频| 色综合婷婷激情| 欧美日韩乱码在线| 看十八女毛片水多多多| 久久久久九九精品影院| 精品一区二区三区人妻视频| 国产精品美女特级片免费视频播放器| 婷婷精品国产亚洲av在线| 国产v大片淫在线免费观看| 久久婷婷人人爽人人干人人爱| 国产aⅴ精品一区二区三区波| 长腿黑丝高跟| 中文字幕av在线有码专区| 欧美日韩乱码在线| 精品福利观看| 亚洲av免费高清在线观看| 此物有八面人人有两片| 99久国产av精品| 国产91精品成人一区二区三区| 日日夜夜操网爽| 狠狠狠狠99中文字幕| 国产午夜福利久久久久久| 亚洲成a人片在线一区二区| 精品国产三级普通话版| 淫妇啪啪啪对白视频| 丰满乱子伦码专区| av天堂中文字幕网| 成年免费大片在线观看| 天堂√8在线中文| .国产精品久久| 免费看av在线观看网站| 有码 亚洲区| 国产精品亚洲一级av第二区| 免费看日本二区| 国产精品日韩av在线免费观看| 麻豆成人av在线观看| 干丝袜人妻中文字幕| 国产黄a三级三级三级人| 日本成人三级电影网站| 欧美中文日本在线观看视频| 亚洲精品粉嫩美女一区| 波野结衣二区三区在线| 欧美中文日本在线观看视频| 成年女人永久免费观看视频| 全区人妻精品视频| 校园春色视频在线观看| av在线天堂中文字幕| 麻豆久久精品国产亚洲av| 高清在线国产一区| 亚洲av第一区精品v没综合| 97碰自拍视频| 天堂网av新在线| 精品久久久久久久末码| 美女大奶头视频| 色综合婷婷激情| 国产欧美日韩一区二区精品| 丰满的人妻完整版| 日本熟妇午夜|