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

    水泵進(jìn)水池底部壓力脈動(dòng)特性試驗(yàn)

    2017-12-04 09:22:12宋希杰查智力嚴(yán)天序黃佳衛(wèi)
    關(guān)鍵詞:軸流泵漩渦水池

    宋希杰 劉 超 楊 帆 查智力 嚴(yán)天序 黃佳衛(wèi)

    (揚(yáng)州大學(xué)水利與能源動(dòng)力工程學(xué)院, 揚(yáng)州 225100)

    水泵進(jìn)水池底部壓力脈動(dòng)特性試驗(yàn)

    宋希杰 劉 超 楊 帆 查智力 嚴(yán)天序 黃佳衛(wèi)

    (揚(yáng)州大學(xué)水利與能源動(dòng)力工程學(xué)院, 揚(yáng)州 225100)

    為了研究泵站進(jìn)水池中喇叭管下方漩渦產(chǎn)生的機(jī)理,探索進(jìn)水池底部喇叭管下方水流壓力脈動(dòng)與漩渦產(chǎn)生的關(guān)系,在進(jìn)水池底部喇叭管口下方布置20個(gè)壓力脈動(dòng)監(jiān)測(cè)點(diǎn),在轉(zhuǎn)速2 200 r/min時(shí)進(jìn)行了多個(gè)不同工況點(diǎn)的壓力脈動(dòng)試驗(yàn)。通過(guò)對(duì)水流壓力脈動(dòng)的時(shí)域特性、頻域特性及進(jìn)水池底部喇叭管下方的壓力分布進(jìn)行分析,揭示了軸流泵喇叭管下方水流壓力脈動(dòng)的規(guī)律及壓力分布與漩渦之間的關(guān)系。進(jìn)水池底部壓力脈動(dòng)的主頻是葉輪轉(zhuǎn)頻的2倍,說(shuō)明軸流泵進(jìn)水池底部喇叭管下方水流的壓力脈動(dòng)仍然受到葉輪旋轉(zhuǎn)的影響;大流量工況下,喇叭管下方容易產(chǎn)生漩渦,且發(fā)生頻率高,持續(xù)時(shí)間長(zhǎng),主要集中在喇叭管下方右前方的位置;小流量工況下,喇叭管下方也能產(chǎn)生漩渦,但發(fā)生頻率低,持續(xù)時(shí)間短,主要集中在喇叭管內(nèi)部正前方的位置。在不同工況下發(fā)生漩渦的位置與進(jìn)水池底部喇叭管下方低壓區(qū)位置相對(duì)應(yīng)。

    軸流泵; 進(jìn)水池底部; 壓力傳感器; 壓力脈動(dòng); 模型試驗(yàn)

    引言

    軸流泵流量大、揚(yáng)程低,廣泛應(yīng)用于灌排泵站以及南水北調(diào)等大規(guī)??缌饔蛘{(diào)水工程。進(jìn)水池及水泵內(nèi)部復(fù)雜的運(yùn)行狀況不僅會(huì)產(chǎn)生陣發(fā)性渦帶,導(dǎo)致水泵性能的大幅下降,甚至引起機(jī)組產(chǎn)生嚴(yán)重的振動(dòng)和噪聲,嚴(yán)重影響泵站的安全穩(wěn)定運(yùn)行[1-2]。

    近年來(lái)國(guó)內(nèi)外學(xué)者針對(duì)水泵的壓力脈動(dòng)做了大量試驗(yàn)研究[3-5],通過(guò)CFD數(shù)值模擬及模型試驗(yàn),對(duì)軸流泵內(nèi)部非穩(wěn)定流場(chǎng)的壓力脈動(dòng)做了大量的研究,在水泵進(jìn)出口管道上布置壓力脈動(dòng)監(jiān)測(cè)點(diǎn),對(duì)軸流泵內(nèi)部水流的壓力脈動(dòng)特性進(jìn)行了試驗(yàn),對(duì)其有了充分的了解[6-9]。但壓力脈動(dòng)不僅發(fā)生在泵內(nèi)部,在水流進(jìn)入喇叭管之前就已經(jīng)發(fā)生,目前針對(duì)泵外部壓力脈動(dòng)對(duì)漩渦產(chǎn)生影響的研究很少。

    針對(duì)這一問(wèn)題,本文通過(guò)進(jìn)水池底部喇叭管下方壓力脈動(dòng)特性試驗(yàn)研究,揭示水泵葉輪旋轉(zhuǎn)對(duì)進(jìn)水池底部壓力脈動(dòng)的影響。

    1 試驗(yàn)裝置

    立式軸流泵裝置的能量性能試驗(yàn)和壓力脈動(dòng)試驗(yàn)均在Ф120 mm立式軸流泵裝置試驗(yàn)臺(tái)上進(jìn)行。 葉輪直徑為120 mm,葉頂間隙為0.1 mm,輪轂直徑為48 mm,葉片數(shù)為4,葉片安放角為0°,導(dǎo)葉數(shù)為7。整個(gè)試驗(yàn)臺(tái)由開(kāi)敞式進(jìn)水池、ISW150-200A型不銹鋼離心泵、PVC管道、穩(wěn)壓圓柱形水箱、D341型法蘭式不銹鋼軟密封蝶閥組成,如圖1所示。在進(jìn)水池底部喇叭管下方布置了20個(gè)壓力脈動(dòng)監(jiān)測(cè)點(diǎn),測(cè)孔直徑為5 mm,安裝20個(gè)壓力傳感器,測(cè)點(diǎn)區(qū)域及具體位置如圖2所示。

    圖1 立式軸流泵模型裝置試驗(yàn)臺(tái)Fig.1 Experiment bench of axial-flow pump model device

    圖2 壓力脈動(dòng)測(cè)試位置圖Fig.2 Sketches of pressure pulsation test position1.測(cè)點(diǎn)區(qū)域 2.喇叭管

    在能量性能試驗(yàn)中,流量測(cè)量采用上海光華儀表有限公司生產(chǎn)的LDG-SDN150型電磁流量計(jì),揚(yáng)程測(cè)量采用EJA型智能壓差傳感器,扭矩及轉(zhuǎn)速測(cè)試采用JCO型轉(zhuǎn)速轉(zhuǎn)矩傳感器,各傳感器輸出端與Power link JW-3型顯示儀相連。在壓力脈動(dòng)試驗(yàn)中,采用成都泰斯特公司研制的高精密度、高穩(wěn)定度CY302系列智能數(shù)字壓力傳感器,其表壓測(cè)試范圍為0~60 MPa,測(cè)試精度為0.1%,采樣時(shí)間間隔為0.001~1 s。CY302標(biāo)準(zhǔn)輸出為485總線(xiàn),通過(guò)該轉(zhuǎn)換器與集線(xiàn)器連接,一臺(tái)集線(xiàn)器可同時(shí)接入20支CY302傳感器,集線(xiàn)器通過(guò)網(wǎng)線(xiàn)與計(jì)算機(jī)相連,壓力測(cè)試軟件采用Smart Sensor壓力測(cè)試系統(tǒng)。

    2 試驗(yàn)結(jié)果與分析

    本試驗(yàn)的目的是研究進(jìn)水池底部喇叭管下方水流的壓力脈動(dòng)特性以及探究漩渦出現(xiàn)的位置與壓力分布之間的關(guān)系。通過(guò)觀察發(fā)現(xiàn)進(jìn)水池水深為300 mm時(shí)漩渦產(chǎn)生情況最明顯,所以試驗(yàn)方案進(jìn)水池水深選為300 mm。本試驗(yàn)方案是在轉(zhuǎn)速為2 200 r/min的情況下,進(jìn)行0.8Qd、Qd和1.2Qd(Qd為設(shè)計(jì)工況點(diǎn))3個(gè)不同工況點(diǎn)下的壓力脈動(dòng)試驗(yàn)。

    2.1 試驗(yàn)誤差分析

    采用標(biāo)準(zhǔn)誤差對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行可靠性評(píng)定。由于試驗(yàn)數(shù)據(jù)量很大,選取部分實(shí)測(cè)數(shù)據(jù)進(jìn)行試驗(yàn)誤差分析。

    (1)

    (2)

    式中m——測(cè)點(diǎn)數(shù)據(jù)組數(shù)

    σ——試驗(yàn)測(cè)量數(shù)據(jù)的標(biāo)準(zhǔn)誤差

    通過(guò)對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)誤差處理,得到σ=±0.08,試驗(yàn)誤差為±0.08。

    2.2 能量性能試驗(yàn)結(jié)果

    為了確定模型泵裝置的設(shè)計(jì)工況點(diǎn),通過(guò)試驗(yàn)得到模型泵裝置在轉(zhuǎn)速為2 200 r/min時(shí)的性能曲線(xiàn)。如圖3所示,泵裝置設(shè)計(jì)工況點(diǎn)Qd為32 L/s。

    其中,壓力脈動(dòng)測(cè)試試驗(yàn)在不同工況點(diǎn)下的性能參數(shù)數(shù)據(jù)如表1所示。

    圖3 立式軸流泵模型裝置性能曲線(xiàn)Fig.3 Performance curves of model pump

    工況點(diǎn)流量Q/(L·s-1)揚(yáng)程H/m效率η/%轉(zhuǎn)速N/(r·min-1)扭矩/(N·m)功率/kW08Qd257224677182200372086Qd32121628767220030607112Qd385605375122200248054

    2.3 壓力脈動(dòng)試驗(yàn)結(jié)果

    由于測(cè)試環(huán)境的變化,采集設(shè)備在空載時(shí)可能存在零點(diǎn)漂移值,為了能更加準(zhǔn)確地測(cè)試壓力的實(shí)際值,在測(cè)試之前采集零點(diǎn)值,在采集過(guò)程中扣除基準(zhǔn)零點(diǎn)值。因?yàn)殇鰷u一旦產(chǎn)生持續(xù)時(shí)間很短,為準(zhǔn)確地測(cè)量出漩渦發(fā)生時(shí)的壓力變化,設(shè)置采樣時(shí)長(zhǎng)t為5 s,為準(zhǔn)確地采集到進(jìn)水池底部喇叭管下方的壓力,設(shè)置采樣時(shí)間間隔ts為傳感器最小采集時(shí)間,即ts=1 ms,針對(duì)3 000 r/min以下的轉(zhuǎn)速,均在本傳感器采集范圍內(nèi),所以采樣頻率fs=1/ts=1 000 Hz能夠滿(mǎn)足2 200 r/min轉(zhuǎn)速下的壓力采集。通過(guò)參考相關(guān)關(guān)于壓力脈動(dòng)試驗(yàn)的文獻(xiàn),在一個(gè)葉輪旋轉(zhuǎn)周期內(nèi)壓力脈動(dòng)數(shù)據(jù)采集點(diǎn)不低于16個(gè)[10]。本試驗(yàn)水泵裝置,葉輪旋轉(zhuǎn)1周的時(shí)間T=1/N=0.027 s,而T/ts=27gt;16,所以能滿(mǎn)足數(shù)據(jù)采集條件。

    2.3.1壓力脈動(dòng)時(shí)域特性分析

    根據(jù)進(jìn)水池底部喇叭管下方漩渦發(fā)生時(shí)所對(duì)應(yīng)的位置,在不同工況點(diǎn)下,選取監(jiān)測(cè)點(diǎn)7、16、20進(jìn)行壓力脈動(dòng)時(shí)域特性分析。為了更好地分析壓力脈動(dòng)的時(shí)域特性,引入無(wú)量綱壓力脈動(dòng)系數(shù)

    (3)

    其中

    式中p——各測(cè)點(diǎn)的瞬時(shí)壓力,kPa

    ρ——水的密度,kg/m3

    Utip——葉頂?shù)膱A周速度,m/s

    r——葉輪半徑

    定義葉輪旋轉(zhuǎn)周期數(shù)

    式中ti——任一點(diǎn)信號(hào)的采集時(shí)間

    本文時(shí)域圖中均為5個(gè)葉輪旋轉(zhuǎn)周期的壓力脈動(dòng)數(shù)據(jù)。

    (1)設(shè)計(jì)流量工況

    設(shè)計(jì)流量工況點(diǎn)下3個(gè)監(jiān)測(cè)點(diǎn)的壓力脈動(dòng)時(shí)域特性如圖4所示,有明顯的波峰和波谷,峰峰值達(dá)到0.004 5,在一個(gè)葉輪旋轉(zhuǎn)周期內(nèi)具有2個(gè)波峰和波谷,測(cè)點(diǎn)7、16、20的壓力脈動(dòng)幅值隨著距喇叭管中心距離的增大而增大,測(cè)點(diǎn)20最先達(dá)到壓力脈動(dòng)波峰值,測(cè)點(diǎn)16和測(cè)點(diǎn)20壓力脈動(dòng)時(shí)域特性相似性很大。單個(gè)壓力脈動(dòng)周期內(nèi),3個(gè)特征測(cè)點(diǎn)的壓力脈動(dòng)系數(shù)逐漸降低。在2個(gè)壓力脈動(dòng)周期交接的地方,測(cè)點(diǎn)7會(huì)出現(xiàn)微小的波動(dòng),而測(cè)點(diǎn)16和測(cè)點(diǎn)20會(huì)出現(xiàn)較大的波動(dòng)。測(cè)點(diǎn)20壓力脈動(dòng)最大幅值為測(cè)點(diǎn)7最大幅值的1.24倍。設(shè)計(jì)流量工況下喇叭管下方的壓力脈動(dòng)整體在一定范圍內(nèi)。

    (2)大流量工況

    圖4 設(shè)計(jì)流量工況監(jiān)測(cè)點(diǎn)壓力脈動(dòng)時(shí)域圖Fig.4 Time domain diagrams of monitoring points under design flow conditions

    圖5 大流量工況監(jiān)測(cè)點(diǎn)壓力脈動(dòng)時(shí)域圖Fig.5 Time domain diagrams of monitoring points under large flow conditions

    大流量工況點(diǎn)下3個(gè)監(jiān)測(cè)點(diǎn)的壓力脈動(dòng)時(shí)域特性如圖5所示,在一個(gè)葉輪旋轉(zhuǎn)周期內(nèi)具有2個(gè)波峰和波谷;大流量工況下壓力脈動(dòng)波動(dòng)時(shí)域特性曲線(xiàn)會(huì)明顯紊亂,不同測(cè)點(diǎn)在不同葉輪旋轉(zhuǎn)周期內(nèi)的壓力脈動(dòng)幅值不同,在大流量工況下測(cè)點(diǎn)7的峰峰值最大,達(dá)到0.01,測(cè)點(diǎn)16峰峰值0.008 2,測(cè)點(diǎn)20的峰峰值達(dá)到0.008,可以看出隨著半徑的增大壓力脈動(dòng)幅值逐漸減小,同一測(cè)點(diǎn)的不同壓力脈動(dòng)周期內(nèi)壓力脈動(dòng)曲線(xiàn)變化也不同,相較于測(cè)點(diǎn)7,測(cè)點(diǎn)16和測(cè)點(diǎn)20的不同壓力脈動(dòng)周期內(nèi)的壓力脈動(dòng)相似性差,說(shuō)明大流量工況進(jìn)水池底部喇叭管下方壓力梯度大,存在劇烈的能量交換,這與此工況下泵裝置運(yùn)行狀況差及容易產(chǎn)生漩渦相符。

    (3)小流量工況

    小流量工況點(diǎn)下3個(gè)監(jiān)測(cè)點(diǎn)的壓力脈動(dòng)時(shí)域特性如圖6所示,3個(gè)時(shí)域圖均為近似正弦波形圖,壓力脈動(dòng)曲線(xiàn)很有規(guī)律,峰峰值為0.006,在一個(gè)葉輪旋轉(zhuǎn)周期內(nèi)具有2個(gè)波峰和2個(gè)波谷。小流量工況下壓力脈動(dòng)幅值隨著與喇叭管中心距離的增大壓力脈動(dòng)幅值逐漸減小,在喇叭管中心附近的測(cè)點(diǎn)7的峰峰值為測(cè)點(diǎn)20峰峰值的1.2倍,且在喇叭管附近的壓力脈動(dòng)時(shí)域圖曲線(xiàn)明顯比測(cè)點(diǎn)20的壓力脈動(dòng)時(shí)域曲線(xiàn)圖整齊,規(guī)律性強(qiáng)。

    圖6 小流量工況監(jiān)測(cè)點(diǎn)壓力脈動(dòng)時(shí)域圖Fig.6 Time domain diagrams of monitoring points under small flow conditions

    可以看出,隨著流量的增大進(jìn)水池底部的壓力脈動(dòng)幅值先減小后增大,在設(shè)計(jì)工況點(diǎn)時(shí)峰峰值最小,在大流量工況點(diǎn)時(shí)的峰峰值最大,說(shuō)明進(jìn)水池底部的壓力脈動(dòng)受到葉輪旋轉(zhuǎn)的影響,在大流量工況下進(jìn)水池內(nèi)水流能量交換強(qiáng)度大,這與大流量時(shí)水泵裝置振動(dòng)和噪聲顯著加劇的現(xiàn)象相符,這是因?yàn)楸醚b置振動(dòng)和噪聲一部分是由于壓力脈動(dòng)造成的[8]。

    2.3.2壓力脈動(dòng)頻域特性分析

    將壓力脈動(dòng)試驗(yàn)數(shù)據(jù)經(jīng)過(guò)加窗傅里葉函數(shù)變換處理[10-13],分別得到1.2Qd、Qd、0.8Qd3個(gè)不同工況點(diǎn)下的壓力脈動(dòng)頻域圖,如圖7~9所示。對(duì)不同流量工況下壓力脈動(dòng)頻域分析發(fā)現(xiàn),各測(cè)點(diǎn)壓力脈動(dòng)的主頻為水泵葉輪轉(zhuǎn)頻的2倍,并且水泵葉頻為各測(cè)點(diǎn)壓力脈動(dòng)主頻的2倍,這與葉輪內(nèi)部進(jìn)出口壓力脈動(dòng)的主頻等于葉頻有所不同[13-17],這是由于葉輪內(nèi)部的壓力脈動(dòng)直接受到葉片旋轉(zhuǎn)的作用。但是仍然說(shuō)明進(jìn)水池底部喇叭管口下方水流的壓力脈動(dòng)與葉輪旋轉(zhuǎn)有很大的關(guān)系。

    (1)設(shè)計(jì)流量工況

    如圖7所示,在設(shè)計(jì)流量工況下右側(cè)測(cè)點(diǎn)1、2、4、5、8、9、11、20的壓力脈動(dòng)幅值明顯高于左側(cè)測(cè)點(diǎn)3、6、7、10、13、14、15的壓力脈動(dòng)幅值,壓力脈動(dòng)系數(shù)最大值為0.000 56,不同測(cè)點(diǎn)壓力脈動(dòng)系數(shù)曲線(xiàn)都有規(guī)律的波動(dòng)減小,在3倍葉輪轉(zhuǎn)頻處會(huì)出現(xiàn)壓力脈動(dòng)二次頻,說(shuō)明喇叭管下方進(jìn)水池底部右側(cè)水流壓力脈動(dòng)比左側(cè)水流壓力脈動(dòng)劇烈,進(jìn)水池底部喇叭管下方在受葉輪旋轉(zhuǎn)作用的影響逐漸減弱。

    圖7 設(shè)計(jì)流量工況監(jiān)測(cè)點(diǎn)壓力脈動(dòng)頻域圖Fig.7 Frequency domain diagram of monitoring points under design flow conditions

    (2)大流量工況

    如圖8所示,在大流量工況下會(huì)發(fā)現(xiàn)壓力脈動(dòng)頻域圖中不同測(cè)點(diǎn)之間壓力脈動(dòng)主頻幅值出現(xiàn)了明顯的波動(dòng),測(cè)點(diǎn)5~13的壓力脈動(dòng)幅值明顯低于其他測(cè)點(diǎn)的壓力脈動(dòng)幅值,測(cè)點(diǎn)5~13主要集中在喇叭管下方來(lái)水側(cè),這說(shuō)明大流量工況下,整體上喇叭管下方背水側(cè)壓力脈動(dòng)變化要強(qiáng)于來(lái)流側(cè)的壓力脈動(dòng)變化,壓力脈動(dòng)幅值大說(shuō)明在此區(qū)域能量交換大,區(qū)域壓力梯度大。測(cè)點(diǎn)7、16、20壓力脈動(dòng)的主頻出現(xiàn)在葉輪3倍轉(zhuǎn)頻處,說(shuō)明在進(jìn)水池底部喇叭管下方右前側(cè)出現(xiàn)最大壓力脈動(dòng)變化要比其他測(cè)點(diǎn)出現(xiàn)最大壓力脈動(dòng)變化晚1個(gè)葉輪旋轉(zhuǎn)周期,這導(dǎo)致明顯的壓差區(qū),為漩渦的產(chǎn)生創(chuàng)造了條件。

    圖8 大流量工況監(jiān)測(cè)點(diǎn)壓力脈動(dòng)頻域圖Fig.8 Frequency domain diagram of monitoring points under large flow conditions

    (3)小流量工況

    如圖9所示,在小流量工況下,壓力脈動(dòng)頻域圖中不同測(cè)點(diǎn)之間壓力脈動(dòng)主頻幅值整體變化均勻,最大壓力脈動(dòng)幅值為0.000 69。次頻幅值小,說(shuō)明小流量工況下壓力脈動(dòng)發(fā)生后,壓力脈動(dòng)的影響消失得快。

    圖9 小流量工況監(jiān)測(cè)點(diǎn)壓力脈動(dòng)頻域圖Fig.9 Frequency domain diagram of monitoring points under small flow conditions

    由3個(gè)不同流量工況下的壓力脈動(dòng)頻域分析可以看出,隨著流量的增大,喇叭管下方進(jìn)水池底部的壓力脈動(dòng)幅度先逐漸減小后逐漸增大,在設(shè)計(jì)工況下對(duì)應(yīng)的各測(cè)點(diǎn)振動(dòng)幅值最小,且振動(dòng)幅值相對(duì)均勻整齊,由于大流量下進(jìn)水池內(nèi)水流湍動(dòng)能大,水流流場(chǎng)復(fù)雜,所以進(jìn)水池底部壓力脈動(dòng)在大流量工況下對(duì)應(yīng)的各測(cè)點(diǎn)振動(dòng)幅值最大。

    2.3.3進(jìn)水池底部壓力分布分析

    相關(guān)研究已經(jīng)采用流線(xiàn)的形式描述了流場(chǎng)中漩渦與流速之間的關(guān)系[14-18],而本文目的是探究漩渦與壓力之間的關(guān)系。

    為便于分析,根據(jù)來(lái)流方向定義左側(cè)和右側(cè),左側(cè)來(lái)流方向與葉輪旋轉(zhuǎn)方向一致稱(chēng)為順?biāo)畟?cè),同理右側(cè)稱(chēng)為逆水側(cè)[19-20]。

    在試驗(yàn)過(guò)程中當(dāng)觀察到漩渦發(fā)生時(shí)開(kāi)始采集壓力。利用壓力測(cè)試軟件Smart Sensor將采集的壓力進(jìn)行加窗傅里葉變換并導(dǎo)出。在每一個(gè)工況下的壓力數(shù)據(jù)中,分別選取3組當(dāng)漩渦發(fā)生時(shí)所對(duì)應(yīng)的1個(gè)葉輪旋轉(zhuǎn)周期內(nèi)20個(gè)測(cè)點(diǎn)的平均壓力值,利用軟件Tecplot進(jìn)行插值處理獲得進(jìn)水池底部不同時(shí)刻的壓力云圖,如圖10~13所示。對(duì)于喇叭管吸入水流的進(jìn)水結(jié)構(gòu),實(shí)際上是水流匯入的紊流混摻區(qū),極易發(fā)生漩渦。

    (1)設(shè)計(jì)流量工況

    圖10為設(shè)計(jì)流量工況下喇叭管下方進(jìn)水池底部的壓力隨時(shí)間的壓力分布云圖。可以看出,此工況下進(jìn)水池底部壓力分布均勻。左側(cè)整體的壓力略高于右側(cè)的壓力,高壓區(qū)集中在喇叭管中心處,低壓區(qū)主要分布在四周,隨著時(shí)間的改變,進(jìn)水池底部的壓力也在改變,但壓力分布形式基本不變。這是由于在設(shè)計(jì)工況點(diǎn)下喇叭管下方流場(chǎng)的能量分布均勻。

    圖10 設(shè)計(jì)流量工況下進(jìn)水池底部壓力云圖Fig.10 Pressure cloud charts at bottom of suction passage under design flow conditions

    圖11 大流量工況下進(jìn)水池底部壓力云圖Fig.11 Pressure cloud charts at bottom of suction passage under large flow conditions

    (2)大流量工況

    圖11為喇叭管下方進(jìn)水池底部在大流量工況下的壓力分布云圖,來(lái)水側(cè)的壓力要高于背水側(cè),高壓區(qū)主要集中在右后側(cè),低壓區(qū)集中在右前側(cè)的位置,如圖11所標(biāo)注處。水流從水平來(lái)流方向轉(zhuǎn)向垂直方向進(jìn)入喇叭管,在慣性作用下加快了能量的交換;通過(guò)壓力脈動(dòng)頻域分析得到進(jìn)水池底部喇叭管下方仍受到葉輪旋轉(zhuǎn)的影響,由于大流量工況下順?biāo)畟?cè)效應(yīng)和逆水側(cè)效應(yīng)[8],葉輪進(jìn)口的水流受水流入角的影響不斷擴(kuò)大,左側(cè)的入流角減小,右側(cè)的入流角增大,使水流條件異常紊亂。由圖11可以看出,伴隨著能量的交換,右后側(cè)的壓能不斷增加,右前側(cè)的壓力降低導(dǎo)致旋轉(zhuǎn)動(dòng)能積累,最終形成漩渦。這表明在大流量工況下,喇叭管下方壓力梯度大,水流能量交換變化大,水流異常紊亂形成較強(qiáng)的漩渦,進(jìn)水條件惡化,軸流泵運(yùn)行效率的下降能夠充分證明漩渦的產(chǎn)生,嚴(yán)重影響水泵機(jī)組的運(yùn)行安全,應(yīng)采取消渦措施進(jìn)行消渦處理。

    圖12 大流量工況進(jìn)水池底部漩渦位置圖Fig.12 Location map of vortex at bottom of suction passage under large flow conditions

    在大流量工況下,隨著水流流態(tài)的惡化,在喇叭管下方進(jìn)水池底板開(kāi)始形成漩渦,如圖12所示,順著水流方向迅速向上發(fā)展,延伸至導(dǎo)水錐附近。渦管形成時(shí)間在1 s左右,持續(xù)時(shí)間為3~5 s,并且不斷地在圖12標(biāo)出區(qū)域不斷移動(dòng)。在漩渦出現(xiàn)期間,會(huì)明顯感覺(jué)到泵裝置噪聲強(qiáng)度增大。根據(jù)圖11中低壓區(qū)的位置和圖12漩渦發(fā)生的位置,可以發(fā)現(xiàn)大流量工況下喇叭管下方漩渦發(fā)生的位置與進(jìn)水池底部喇叭管下方的低壓區(qū)位置一致。

    (3)小流量工況

    如圖13所示,在小流量工況下,進(jìn)水側(cè)整體壓力略高于背水側(cè)的壓力,在喇叭管中心后側(cè)存在一個(gè)低壓區(qū),在喇叭管中心的左側(cè)及右后側(cè)各存在一個(gè)高壓區(qū),在小流量工況下喇叭管入口處的水流受順?biāo)畟?cè)效應(yīng)和逆水側(cè)效應(yīng)的影響減小,喇叭管口左右兩側(cè)的能量分布均勻,導(dǎo)致小流量工況下低壓區(qū)的位置和大流量工況下低壓區(qū)位置從右前方轉(zhuǎn)移到喇叭管正前方。

    小流量工況下水流流態(tài)平穩(wěn),較大流量工況下感覺(jué)到泵裝置產(chǎn)生的噪聲明顯減小,在喇叭管中心后側(cè)出現(xiàn)一條極細(xì)的渦帶,持續(xù)時(shí)間短,且出現(xiàn)頻率很小,難以捕捉,漩渦出現(xiàn)的位置與喇叭管中心的距離增大,不過(guò),此工況下漩渦發(fā)生的位置與圖13中低壓區(qū)的位置一致。

    圖13 小流量工況進(jìn)水池底部壓力云圖Fig.13 Pressure cloud charts at bottom of suction passage under small flow conditions

    在不同工況下,運(yùn)行工況改變時(shí),旋轉(zhuǎn)動(dòng)能變化,易于積聚,導(dǎo)致壓能減少,很容易產(chǎn)生漩渦。在大流量工況下,喇叭管下方附底渦出現(xiàn)的次數(shù)最多,漩渦持續(xù)時(shí)間長(zhǎng),相鄰漩渦之間時(shí)間間隔短,附底渦產(chǎn)生的位置靠近導(dǎo)水錐,影響水泵的安全運(yùn)行,泵裝置振動(dòng)強(qiáng)度大。在設(shè)計(jì)工況下水流流動(dòng)平穩(wěn),附底渦不易產(chǎn)生。在小流量工況下,雖然會(huì)產(chǎn)生附底渦,但漩渦出現(xiàn)頻率很低,且持續(xù)時(shí)間短。漩渦的發(fā)生嚴(yán)重影響水泵機(jī)組的安全運(yùn)行,應(yīng)采取消渦措施進(jìn)行消渦處理。

    3 結(jié)論

    (1)根據(jù)壓力脈動(dòng)時(shí)域特性分析得到,隨著流量的增加進(jìn)水池底部水流的壓力脈動(dòng)幅值先減小后增大,大流量工況下水流壓力脈動(dòng)幅值最大,這與設(shè)計(jì)流量工況點(diǎn)下未產(chǎn)生漩渦和漩渦主要發(fā)生在大流量工況點(diǎn)下的現(xiàn)象一致。

    (2)進(jìn)水池底部的壓力脈動(dòng)受到水泵葉輪旋轉(zhuǎn)影響,進(jìn)水池底部水流壓力脈動(dòng)的主頻為水泵葉輪轉(zhuǎn)頻的2倍,為水泵葉頻的0.5倍。

    (3)在設(shè)計(jì)流量和小流量工況下進(jìn)水池底部喇叭管下方壓力分布變化不大;在大流量工況下,壓力分布變化大,在喇叭管下方右前側(cè)的位置存在低壓區(qū)。在大流量工況下漩渦發(fā)生在喇叭管內(nèi)右前方的位置,持續(xù)時(shí)間長(zhǎng),發(fā)生次數(shù)多。小流量工況下漩渦發(fā)生在喇叭管內(nèi)正前方的位置,持續(xù)時(shí)間短,發(fā)生次數(shù)少,可以發(fā)現(xiàn)漩渦發(fā)生的位置與對(duì)應(yīng)工況下壓力分布中的低壓區(qū)位置相對(duì)應(yīng),說(shuō)明漩渦發(fā)生的位置與進(jìn)水池底部喇叭管下方低壓區(qū)的位置相一致。

    1 劉超.軸流泵系統(tǒng)技術(shù)創(chuàng)新與發(fā)展分析[J/OL].農(nóng)業(yè)機(jī)械學(xué)報(bào),2015, 46(6): 49-59.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1amp;file_no=20150608amp;journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2015.06.008.

    LIU Chao. Researches and developments of axial-flow pump system[J/OL].Transactions of the Chinese Society for Agricultural Machinery,2015, 46(6): 49-59. (in Chinese)

    2 何耘.水泵進(jìn)水池漩渦研究的主要進(jìn)展[J].水力發(fā)電學(xué)報(bào),2004,23(5):92-96.

    HE Yun. Main advances of research on vortex in pump sumps[J].Journal of Hydroelectric Engineering, 2004,23(5):92-96. (in Chinese)

    3 施衛(wèi)東,冷洪飛,張德勝,等.軸流泵內(nèi)部流場(chǎng)壓力脈動(dòng)性能預(yù)測(cè)與試驗(yàn)[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2011,42(5):44-48.

    SHI Weidong, LENG Hongfei, ZHANG Desheng,et al. Performance prediction and experiment for pressure fluctuation of interior flow in axial-flow pump[J].Transactions of the Chinese Society for Agricultural Machinery, 2011,42(5):44-48. (in Chinese)

    4 劉超,成立,周濟(jì)人,等.水泵站開(kāi)敞進(jìn)水池三維紊流數(shù)值模擬[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2002,33(6):53-55.

    LIU Chao, CHENG Li, ZHOU Jiren,et al. Numerical simulation of three-dimensional turbulent flow for opening pump sump[J].Transactions of the Chinese Society for Agricultural Machinery, 2002,33(6):53-55. (in Chinese)

    5 叢國(guó)輝,王福軍.湍流模型在泵站進(jìn)水池漩渦模擬中的適用性研究[J].農(nóng)業(yè)工程學(xué)報(bào),2008,24(6):31-35.

    CONG Guohui, WANG Fujun. Applicability of turbulence models in numerical simulation of vortex flow in pump sump[J].Transactions of the CSAE, 2008,24(6):31-35.(in Chinese)

    6 楊帆,劉超,湯方平,等.箱涵式進(jìn)水流道的立式軸流泵裝置水動(dòng)力特性分析[J].農(nóng)業(yè)工程學(xué)報(bào),2014,30(4):62-69.

    YANG Fan, LIU Chao, TANG Fangping, et al. Analysis of hydraulic performance for vertical axial-flow pumping system with cube-type inlet passage[J]. Transactions of the CSAE,2014, 30(4): 62-69. (in Chinese)

    7 成立,劉超,周濟(jì)人.基于RNG湍流模型的雙向泵站出水流道流動(dòng)計(jì)算[J].水科學(xué)進(jìn)展, 2004,15(1):109-112.

    8 仇寶云, 黃季艷, 林海江, 等. 立式軸流泵出水流道流場(chǎng)試驗(yàn)研究[J], 機(jī)械工程學(xué)報(bào),2005,41(11):115-126.

    9 張麗萍.非均勻來(lái)流條件下軸流泵內(nèi)部壓力脈動(dòng)數(shù)值模擬研究[D].揚(yáng)州:揚(yáng)州大學(xué),2013.

    10 楊帆,金燕,劉超,等. 雙向潛水貫流泵裝置性能試驗(yàn)與數(shù)值分析[J].農(nóng)業(yè)工程學(xué)報(bào),2012,28(16):60-67.

    YANG Fan,JIN Yan,LIU Chao,et al. Numerical analysis and performance test on diving tubular pumping system with symmetric aerofoil blade[J]. Transactions of the CSAE,2012,28(16):60-67. (in Chinese)

    11 SHINJI E, HIROYUKI T, HIDETOSHI H, et al. Characteristics of flow field and pressure fluctuation in complex turbulent flow in the third elbow of a triple elbow piping with small curvature radius in three-dimensional layout[J]. International Journal of Hydrogen Energy, 2016,41(17):7139-7145.

    12 ARNDT N,ACOSTA A J,BRENNEN C E,et al. Frequency spectrum analysis of high frequency cycle square wave signal based on discrete fourier transform[J].ASME Journal of Turbomachinery,1999,120(3): 157-178.

    13 ARNDT N,ACOSTA A J,BRENNEN C E,et al. Experimental investigation of rotor-stator interaction in a centrifugal pump with several vaned diffusers[J]. ASME Journal of Turbomachinery,1990,112(1): 98-108.

    14 劉超,梁豪杰,金燕,等.立式軸流泵進(jìn)水流場(chǎng)PIV測(cè)量[J/OL].農(nóng)業(yè)機(jī)械學(xué)報(bào),2015,46(8):33-41. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1amp;file_no=20150806amp;journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2015.08.006.

    LIU Chao,LIANG Haojie, JIN Yan, et al.PIV measurements of intake flow field in axial-flow pump[J/OL] .Transactions of the Chinese Society for Agricultural Machinery, 2015,46(8):33-41. (in Chinese)

    15 ARNDT N,ACOSTA A J,BRENEN C E,et al. Rotor-stator interaction in a diffuser pump[J]. ASME Journal of Turbomachinery,1989,111(3): 213 - 221.

    16 張德勝,王海宇,施衛(wèi)東,等.軸流泵多工況壓力脈動(dòng)特性試驗(yàn)[J/OL].農(nóng)業(yè)機(jī)械學(xué)報(bào),2014,45(11):139-144. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?flag=1amp;file_no=20141122amp;journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2014.11.022.

    ZHANG Desheng, WANG Haiyu, SHI Weidong, etal. Experimental investigation of pressure fluctuation with multiple flow rates in scaled axial flow pump[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2014,45(11):139-144. (in Chinese)

    17 王福軍,張玲,張志民,等.軸流泵不穩(wěn)定流場(chǎng)的壓力脈動(dòng)特性研究[J].水利學(xué)報(bào),2007,38(8):1003-1009.

    WANG Fujun, ZHANG Ling, ZHANG Zhimin, et al. Analysis on pressure fluctuation of unsteady flow in axial-flow pump[J]. Journal of Hydraulic Engineering, 2007,38(8):1003-1009. (in Chinese)

    18 黃歡明, 高紅, 沈楓, 等.軸流泵內(nèi)流場(chǎng)的數(shù)值模擬與試驗(yàn)[J].農(nóng)業(yè)機(jī)械學(xué)報(bào), 2008, 39(8):66-69.

    HUANG Huanming, GAO Hong, SHEN Feng, et al.Numerical simulation and experimental validation of the flow field in axial flow pump[J].Transactions of the Chinese Society for Agricultural Machinery, 2008, 39(8):66-69.(in Chinese)

    19 姚捷,施衛(wèi)東,吳蘇青,等. 軸流泵壓力脈動(dòng)數(shù)值計(jì)算與試驗(yàn)[J/OL].農(nóng)業(yè)機(jī)械學(xué)報(bào),2013,44(增刊1): 119-124.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=2013s122amp;flag=1amp;journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2013.S1.022.

    YAO Jie,SHI Weidong,WU Suqing,et al. Numerical calculation and experiment on pressure fluctuation in axial flow pump[J/OL].Transactions of the Chinese Society for Agricultural Machinery,2013,44(Supp.1): 119-124.(in Chinese)

    20 張德勝,施衛(wèi)東,張華,等.不同湍流模型在軸流泵性能預(yù)測(cè)中的應(yīng)用[J].農(nóng)業(yè)工程學(xué)報(bào),2012,28(1): 66-70.

    ZHANG Desheng,SHI Weidong,ZHANG Hua,et al.Application of different turbulence model for predicting performance of axial flow pump[J].Transactions of the CSAE,2012,28(1): 66-70. (in Chinese)

    ExperimentonCharacteristicsofPressureFluctuationatBottomofPumpingSuctionPassage

    SONG Xijie LIU Chao YANG Fan ZHA Zhili YAN Tianxu HUANG Jiawei

    (SchoolofHydraulicEnergyandPowerEngineering,YangzhouUniversity,Yangzhou225100,China)

    In order to study the mechanism of bottom vortex generation in the suction passage of pumping station, and explore the relationship between pressure pulsation and vortex generation under flare pipe at bottom of suction passage, 20 pressure pulsation monitoring points were evenly arranged at the bottom of suction passage in the form of concentric circular ring. At speed of 2 200 r/min, the pressure pulsations experiment under different flow conditions was done. By analyzing the time domain and frequency domain of pressure fluctuation and pressure distribution at the bottom of suction passage, the law of pressure pulsation at the bottom of suction passage of the axial flow pump was revealed. According to the analysis of pressure fluctuation time domain, the pressure fluctuation amplitude of the bottom water flow at the bottom of inlet basin was decreased and then increased with the increase of flow rate. The maximum pressure of the flow pressure was fluctuated under the large flow condition, which was the same as the design flow. No vortex and whirlpool occurred mainly in the case of large flow conditions. The experimental results showed that the pressure pulsation at the bottom of the suction passage was twice as high as that of the impeller, and there was no dominant frequency at impeller frequency, which indicated that the pressure pulsation at the bottom of the suction passage under the bell-tube of the axial flow pump was still affected by impeller rotation. Under large flow conditions, the vortex under the bell-tube was easy to produce with high frequency of occurrence and long duration. The position of the vortex under different working conditions was corresponded to the position of the low pressure zone at the bottom of suction below the bell-tube.

    axial-flow pump; suction passage bottom; pressure sensor; pressure pulsation; model test

    10.6041/j.issn.1000-1298.2017.11.024

    TH311

    A

    1000-1298(2017)11-0196-08

    2017-03-02

    2017-04-28

    “十二五”國(guó)家科技支撐計(jì)劃項(xiàng)目(2015BAD20B01-02)、國(guó)家自然科學(xué)基金項(xiàng)目(51279173、51609210)、江蘇省自然科學(xué)基金項(xiàng)目(BK20150457)、江蘇高校優(yōu)勢(shì)學(xué)科建設(shè)工程項(xiàng)目(PAPD)和2017江蘇省研究生科研創(chuàng)新工程項(xiàng)目(KYCX17_1881)

    宋希杰(1990—),男,博士生,主要從事水利工程研究,E-mail: 1227321904@qq.com

    劉超(1950—),男,教授,主要從事泵站工程理論和應(yīng)用技術(shù)研究,E-mail: liuchao@yzu.edu.cn

    猜你喜歡
    軸流泵漩渦水池
    潛水軸流泵運(yùn)行故障分析與排除研究
    潛水軸流泵電機(jī)運(yùn)行工況的特點(diǎn)及可靠性探討
    小區(qū)的水池
    FF陷控制權(quán)爭(zhēng)奪漩渦
    魚(yú)群漩渦
    中外文摘(2017年19期)2017-10-10 08:28:41
    把住醫(yī)?;鹚亻l門(mén)
    找水池
    中醫(yī)教育陷“量升質(zhì)降”漩渦
    濃縮軸流泵干氣密封改造
    中醫(yī)教育陷“量升質(zhì)降”漩渦
    免费观看在线日韩| 国产精品秋霞免费鲁丝片| 一边亲一边摸免费视频| 99九九在线精品视频| 蜜桃国产av成人99| 1024香蕉在线观看| 亚洲第一区二区三区不卡| 色婷婷av一区二区三区视频| 国产精品久久久久久av不卡| 亚洲精品国产av蜜桃| 人妻 亚洲 视频| 成年美女黄网站色视频大全免费| 国产极品粉嫩免费观看在线| 中文字幕av电影在线播放| 久久久久国产精品人妻一区二区| 日韩精品免费视频一区二区三区| 欧美日韩一区二区视频在线观看视频在线| 久久久久视频综合| 免费日韩欧美在线观看| 精品人妻在线不人妻| 大话2 男鬼变身卡| 国产激情久久老熟女| 亚洲av免费高清在线观看| www.自偷自拍.com| 国产熟女欧美一区二区| 国产97色在线日韩免费| www日本在线高清视频| videos熟女内射| 久久女婷五月综合色啪小说| 成人国产麻豆网| 一区二区三区四区激情视频| 18禁观看日本| 国产高清不卡午夜福利| 少妇的丰满在线观看| 欧美日韩国产mv在线观看视频| 女的被弄到高潮叫床怎么办| 中文字幕人妻丝袜一区二区 | 男女免费视频国产| 免费久久久久久久精品成人欧美视频| 亚洲欧美一区二区三区国产| 亚洲欧美成人精品一区二区| 国产熟女午夜一区二区三区| 日韩成人av中文字幕在线观看| 久久精品久久久久久久性| 欧美老熟妇乱子伦牲交| 国产精品成人在线| 高清不卡的av网站| 极品少妇高潮喷水抽搐| 中文字幕人妻丝袜制服| 日韩av免费高清视频| 亚洲国产精品成人久久小说| 日韩一区二区视频免费看| 男女下面插进去视频免费观看| av网站在线播放免费| 国产老妇伦熟女老妇高清| 日韩视频在线欧美| 1024视频免费在线观看| 久久青草综合色| 一边摸一边做爽爽视频免费| 国产成人午夜福利电影在线观看| 波多野结衣一区麻豆| 校园人妻丝袜中文字幕| 精品人妻一区二区三区麻豆| 欧美精品一区二区大全| 亚洲欧洲精品一区二区精品久久久 | 成年av动漫网址| 国产在线免费精品| 七月丁香在线播放| 免费日韩欧美在线观看| 蜜桃在线观看..| 亚洲一区二区三区欧美精品| 亚洲精品久久午夜乱码| 精品亚洲乱码少妇综合久久| 边亲边吃奶的免费视频| 黑人猛操日本美女一级片| 九色亚洲精品在线播放| 高清视频免费观看一区二区| 少妇人妻 视频| 卡戴珊不雅视频在线播放| 最近2019中文字幕mv第一页| 免费看av在线观看网站| 免费观看无遮挡的男女| xxxhd国产人妻xxx| 春色校园在线视频观看| 精品国产超薄肉色丝袜足j| 1024香蕉在线观看| 18在线观看网站| 亚洲国产欧美在线一区| 日韩成人av中文字幕在线观看| 免费观看性生交大片5| 丝袜脚勾引网站| 中文字幕亚洲精品专区| 一级爰片在线观看| 欧美激情极品国产一区二区三区| 男女下面插进去视频免费观看| 色婷婷av一区二区三区视频| 国产精品三级大全| 亚洲婷婷狠狠爱综合网| 一级黄片播放器| 国产精品国产av在线观看| 日韩一本色道免费dvd| 亚洲综合色网址| av国产久精品久网站免费入址| 亚洲国产欧美日韩在线播放| 成人午夜精彩视频在线观看| 亚洲一码二码三码区别大吗| 中文字幕色久视频| av国产久精品久网站免费入址| 亚洲精品在线美女| 欧美成人精品欧美一级黄| 男人舔女人的私密视频| 在线观看一区二区三区激情| 久久精品熟女亚洲av麻豆精品| 久久精品aⅴ一区二区三区四区 | 久久精品国产亚洲av高清一级| 久久99精品国语久久久| 成年女人毛片免费观看观看9 | 女人高潮潮喷娇喘18禁视频| 国产成人精品一,二区| 热99久久久久精品小说推荐| videosex国产| 91aial.com中文字幕在线观看| 26uuu在线亚洲综合色| xxxhd国产人妻xxx| 国产黄色视频一区二区在线观看| 久久久a久久爽久久v久久| 成人午夜精彩视频在线观看| 国产白丝娇喘喷水9色精品| 免费女性裸体啪啪无遮挡网站| 制服诱惑二区| 成人18禁高潮啪啪吃奶动态图| 国产 精品1| 丝袜喷水一区| 亚洲成色77777| 哪个播放器可以免费观看大片| 午夜激情久久久久久久| 一本久久精品| 国产在线一区二区三区精| 国产成人精品在线电影| 桃花免费在线播放| 午夜91福利影院| 丝瓜视频免费看黄片| 亚洲,欧美精品.| 国产福利在线免费观看视频| 成年人免费黄色播放视频| 亚洲国产毛片av蜜桃av| 高清视频免费观看一区二区| 欧美最新免费一区二区三区| 天天影视国产精品| 天天影视国产精品| 亚洲成色77777| 欧美日韩成人在线一区二区| 久久久国产一区二区| 日韩欧美精品免费久久| 国产成人欧美| 免费在线观看完整版高清| 人成视频在线观看免费观看| 99热网站在线观看| 卡戴珊不雅视频在线播放| 亚洲av综合色区一区| 18禁国产床啪视频网站| 久久久久久久久久久久大奶| 最近最新中文字幕免费大全7| 亚洲色图综合在线观看| 激情视频va一区二区三区| 久久久欧美国产精品| 麻豆av在线久日| 亚洲内射少妇av| 一区二区三区乱码不卡18| 日韩电影二区| 精品久久久久久电影网| 9热在线视频观看99| 亚洲男人天堂网一区| 精品一区二区免费观看| 日韩一区二区三区影片| 丝袜美腿诱惑在线| 亚洲 欧美一区二区三区| 国产野战对白在线观看| 久久99精品国语久久久| 亚洲精品日本国产第一区| 人妻一区二区av| 亚洲欧美成人综合另类久久久| 国产精品嫩草影院av在线观看| 欧美激情 高清一区二区三区| 国产一区二区在线观看av| 9热在线视频观看99| 久久鲁丝午夜福利片| 一本—道久久a久久精品蜜桃钙片| 最近的中文字幕免费完整| 美女高潮到喷水免费观看| 久久久久国产一级毛片高清牌| 建设人人有责人人尽责人人享有的| 国产男女超爽视频在线观看| 多毛熟女@视频| 亚洲av综合色区一区| av.在线天堂| 青春草亚洲视频在线观看| 纯流量卡能插随身wifi吗| 大话2 男鬼变身卡| 亚洲国产欧美网| 少妇被粗大的猛进出69影院| 桃花免费在线播放| 久久久精品区二区三区| 纵有疾风起免费观看全集完整版| 久久 成人 亚洲| 久久狼人影院| 国产成人91sexporn| 777久久人妻少妇嫩草av网站| 国产成人精品福利久久| 91久久精品国产一区二区三区| 999久久久国产精品视频| 满18在线观看网站| 国产免费一区二区三区四区乱码| 久久精品国产亚洲av涩爱| 国产国语露脸激情在线看| 国产成人免费无遮挡视频| 国产无遮挡羞羞视频在线观看| 久久久久久人妻| 国产在线视频一区二区| 两性夫妻黄色片| 黑人猛操日本美女一级片| 婷婷色综合www| 日本91视频免费播放| 成年美女黄网站色视频大全免费| 一本色道久久久久久精品综合| 日韩欧美精品免费久久| 涩涩av久久男人的天堂| 香蕉丝袜av| 国产一区亚洲一区在线观看| 在线 av 中文字幕| 日韩电影二区| 日韩av在线免费看完整版不卡| 制服人妻中文乱码| 青春草国产在线视频| 亚洲,欧美,日韩| 少妇被粗大猛烈的视频| videossex国产| 天天影视国产精品| 男的添女的下面高潮视频| 国产高清国产精品国产三级| 亚洲精品中文字幕在线视频| 久久久久久久久久久免费av| 久久久久国产一级毛片高清牌| 水蜜桃什么品种好| 亚洲视频免费观看视频| 男女边摸边吃奶| 亚洲精品国产一区二区精华液| 丝袜脚勾引网站| 欧美变态另类bdsm刘玥| 香蕉精品网在线| 美女高潮到喷水免费观看| 亚洲国产av新网站| 精品久久久精品久久久| 久久国产精品男人的天堂亚洲| 午夜福利网站1000一区二区三区| 亚洲欧洲日产国产| 一级黄片播放器| 日本黄色日本黄色录像| 性少妇av在线| av电影中文网址| 午夜免费观看性视频| 亚洲激情五月婷婷啪啪| 制服人妻中文乱码| 免费少妇av软件| 日韩精品有码人妻一区| 欧美 日韩 精品 国产| 母亲3免费完整高清在线观看 | 自线自在国产av| 亚洲人成电影观看| 国产精品久久久久久精品电影小说| 色94色欧美一区二区| 国产精品一区二区在线不卡| 国产av一区二区精品久久| av免费在线看不卡| 日日啪夜夜爽| 少妇被粗大猛烈的视频| 美女视频免费永久观看网站| 久久精品亚洲av国产电影网| 人妻系列 视频| 永久网站在线| 王馨瑶露胸无遮挡在线观看| 黄色配什么色好看| 91在线精品国自产拍蜜月| 亚洲人成网站在线观看播放| 精品午夜福利在线看| 老熟女久久久| 街头女战士在线观看网站| 成年人午夜在线观看视频| 色网站视频免费| 亚洲色图综合在线观看| 老汉色av国产亚洲站长工具| 国产一区二区 视频在线| 在线观看免费视频网站a站| 桃花免费在线播放| 精品99又大又爽又粗少妇毛片| 日产精品乱码卡一卡2卡三| 久久久久网色| 成人毛片60女人毛片免费| 欧美日韩国产mv在线观看视频| 欧美日本中文国产一区发布| 日本爱情动作片www.在线观看| 蜜桃在线观看..| 一级片'在线观看视频| av在线app专区| 精品国产乱码久久久久久小说| 青春草国产在线视频| 欧美精品人与动牲交sv欧美| 美女国产高潮福利片在线看| 久久99一区二区三区| 亚洲综合精品二区| 欧美日韩综合久久久久久| 亚洲,一卡二卡三卡| 黄片小视频在线播放| 一边亲一边摸免费视频| 久久精品熟女亚洲av麻豆精品| 一本—道久久a久久精品蜜桃钙片| 久久久久久人妻| 免费观看a级毛片全部| 建设人人有责人人尽责人人享有的| 老鸭窝网址在线观看| 热99国产精品久久久久久7| 久久国产精品大桥未久av| 啦啦啦视频在线资源免费观看| 亚洲成av片中文字幕在线观看 | 国产精品麻豆人妻色哟哟久久| 国产国语露脸激情在线看| 香蕉精品网在线| 久久99蜜桃精品久久| 少妇人妻精品综合一区二区| 视频在线观看一区二区三区| 在线免费观看不下载黄p国产| 午夜福利网站1000一区二区三区| 曰老女人黄片| 亚洲欧美精品自产自拍| 性色avwww在线观看| 久久久久久久久久久免费av| 大片免费播放器 马上看| www.av在线官网国产| 天美传媒精品一区二区| 久久久久精品人妻al黑| 日日爽夜夜爽网站| 性高湖久久久久久久久免费观看| 精品少妇内射三级| 人人妻人人澡人人爽人人夜夜| 青青草视频在线视频观看| 熟女av电影| 成人影院久久| 欧美黄色片欧美黄色片| 欧美日韩视频高清一区二区三区二| 欧美老熟妇乱子伦牲交| 欧美日韩视频高清一区二区三区二| 黄片小视频在线播放| 久久久久人妻精品一区果冻| 国产毛片在线视频| 人人澡人人妻人| 国产毛片在线视频| xxxhd国产人妻xxx| 女性生殖器流出的白浆| 国产成人a∨麻豆精品| 黄色一级大片看看| 亚洲国产欧美网| 久久女婷五月综合色啪小说| 美国免费a级毛片| 久久国产精品大桥未久av| 精品卡一卡二卡四卡免费| 黄色配什么色好看| 好男人视频免费观看在线| 欧美最新免费一区二区三区| 亚洲欧美清纯卡通| 高清黄色对白视频在线免费看| www日本在线高清视频| 在线看a的网站| 亚洲精品日本国产第一区| 亚洲一码二码三码区别大吗| 一二三四在线观看免费中文在| 91午夜精品亚洲一区二区三区| 97精品久久久久久久久久精品| 亚洲精品视频女| 韩国av在线不卡| 亚洲精品久久久久久婷婷小说| 麻豆乱淫一区二区| 91aial.com中文字幕在线观看| 妹子高潮喷水视频| 国产欧美日韩综合在线一区二区| www.自偷自拍.com| 中国三级夫妇交换| 国产无遮挡羞羞视频在线观看| 高清黄色对白视频在线免费看| 两个人看的免费小视频| 久久久久久久国产电影| 成年美女黄网站色视频大全免费| 男女边吃奶边做爰视频| 在线亚洲精品国产二区图片欧美| 日韩av不卡免费在线播放| 咕卡用的链子| 国产探花极品一区二区| 一边亲一边摸免费视频| 亚洲图色成人| 欧美亚洲 丝袜 人妻 在线| 日韩免费高清中文字幕av| 久久久久精品人妻al黑| 有码 亚洲区| 涩涩av久久男人的天堂| 哪个播放器可以免费观看大片| 国产成人aa在线观看| 丝袜美腿诱惑在线| 亚洲 欧美一区二区三区| 少妇被粗大的猛进出69影院| 亚洲成国产人片在线观看| 国产爽快片一区二区三区| 最新的欧美精品一区二区| 国产精品熟女久久久久浪| 成人亚洲精品一区在线观看| 午夜福利影视在线免费观看| 男的添女的下面高潮视频| 久久久久国产精品人妻一区二区| 国产一区二区三区av在线| 久久久久久久大尺度免费视频| 一区二区av电影网| tube8黄色片| 亚洲精品国产一区二区精华液| av视频免费观看在线观看| 精品第一国产精品| 老汉色av国产亚洲站长工具| av国产久精品久网站免费入址| 中文字幕人妻丝袜制服| 最近2019中文字幕mv第一页| 制服诱惑二区| 久久精品国产鲁丝片午夜精品| 亚洲精品乱久久久久久| 久热这里只有精品99| 91久久精品国产一区二区三区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 男女啪啪激烈高潮av片| 一级a爱视频在线免费观看| 日韩欧美一区视频在线观看| 一级片'在线观看视频| 在线观看免费视频网站a站| 久久久久人妻精品一区果冻| 亚洲精品美女久久av网站| 久久午夜综合久久蜜桃| 久久精品熟女亚洲av麻豆精品| 午夜激情久久久久久久| 国产成人午夜福利电影在线观看| 久久久久国产一级毛片高清牌| 亚洲图色成人| 亚洲一级一片aⅴ在线观看| h视频一区二区三区| 在线观看国产h片| 看免费成人av毛片| 久久精品国产综合久久久| 国产成人精品久久久久久| 少妇 在线观看| 日本色播在线视频| 亚洲精品一区蜜桃| 亚洲综合色网址| 18禁国产床啪视频网站| 啦啦啦在线免费观看视频4| 最近中文字幕2019免费版| 久久久久网色| 国产黄频视频在线观看| 九九爱精品视频在线观看| 精品午夜福利在线看| 18+在线观看网站| 亚洲精品av麻豆狂野| 精品国产露脸久久av麻豆| 亚洲一区中文字幕在线| 国产黄频视频在线观看| 日韩制服丝袜自拍偷拍| 啦啦啦在线观看免费高清www| 1024香蕉在线观看| 国产免费视频播放在线视频| 国产熟女欧美一区二区| 亚洲欧美一区二区三区国产| 国产熟女午夜一区二区三区| 老女人水多毛片| 男女边吃奶边做爰视频| 老汉色av国产亚洲站长工具| av电影中文网址| 亚洲激情五月婷婷啪啪| 亚洲国产成人一精品久久久| 久久久久视频综合| 国产精品不卡视频一区二区| 男男h啪啪无遮挡| 国产成人精品无人区| 久久久精品区二区三区| 性少妇av在线| 一区二区日韩欧美中文字幕| 日本爱情动作片www.在线观看| 精品国产露脸久久av麻豆| 日韩 亚洲 欧美在线| 如日韩欧美国产精品一区二区三区| 日日啪夜夜爽| 男女高潮啪啪啪动态图| av有码第一页| av女优亚洲男人天堂| 亚洲色图 男人天堂 中文字幕| 涩涩av久久男人的天堂| 伦理电影免费视频| 国产免费现黄频在线看| 欧美日韩av久久| 香蕉国产在线看| 新久久久久国产一级毛片| 多毛熟女@视频| 黄色一级大片看看| 国产片特级美女逼逼视频| 如日韩欧美国产精品一区二区三区| 国产日韩欧美亚洲二区| 搡老乐熟女国产| 色94色欧美一区二区| 国产 精品1| 欧美bdsm另类| 男女无遮挡免费网站观看| 久久久久久久精品精品| 欧美激情高清一区二区三区 | a级片在线免费高清观看视频| 美女福利国产在线| 老鸭窝网址在线观看| 在线亚洲精品国产二区图片欧美| 又粗又硬又长又爽又黄的视频| 久久精品国产亚洲av涩爱| 中文字幕精品免费在线观看视频| 一区二区三区精品91| 高清在线视频一区二区三区| 成人影院久久| 色网站视频免费| 韩国av在线不卡| 日本vs欧美在线观看视频| 久久久欧美国产精品| 王馨瑶露胸无遮挡在线观看| 亚洲精品日韩在线中文字幕| 伊人久久国产一区二区| 国产高清国产精品国产三级| 成人亚洲欧美一区二区av| 一二三四在线观看免费中文在| 少妇人妻 视频| 考比视频在线观看| 色播在线永久视频| 视频在线观看一区二区三区| 国产乱人偷精品视频| av天堂久久9| 叶爱在线成人免费视频播放| 亚洲精品在线美女| 大码成人一级视频| 成年人免费黄色播放视频| 国产老妇伦熟女老妇高清| 最黄视频免费看| 日韩一区二区视频免费看| 国产成人a∨麻豆精品| 日韩一本色道免费dvd| 国产深夜福利视频在线观看| 亚洲欧洲精品一区二区精品久久久 | 国产一区有黄有色的免费视频| 久久人人爽av亚洲精品天堂| av.在线天堂| 999精品在线视频| av网站在线播放免费| 美女高潮到喷水免费观看| 最近2019中文字幕mv第一页| 日本vs欧美在线观看视频| 亚洲,欧美精品.| av不卡在线播放| 中文字幕人妻熟女乱码| 久久久精品免费免费高清| 精品卡一卡二卡四卡免费| 国产极品天堂在线| 国产精品麻豆人妻色哟哟久久| xxxhd国产人妻xxx| 电影成人av| 伊人久久国产一区二区| 亚洲情色 制服丝袜| 久久精品国产亚洲av高清一级| 亚洲国产成人一精品久久久| 黑人猛操日本美女一级片| 最近中文字幕2019免费版| 国产无遮挡羞羞视频在线观看| 国产成人精品婷婷| 亚洲经典国产精华液单| 亚洲欧美日韩另类电影网站| 一本色道久久久久久精品综合| 亚洲精品美女久久av网站| 777久久人妻少妇嫩草av网站| 国产一区二区在线观看av| 人妻一区二区av| 成年人午夜在线观看视频| 一区二区av电影网| 人妻一区二区av| 国产片特级美女逼逼视频| 国产精品秋霞免费鲁丝片| 国产伦理片在线播放av一区| 在现免费观看毛片| 国产极品粉嫩免费观看在线| 日韩在线高清观看一区二区三区| 国语对白做爰xxxⅹ性视频网站| 亚洲熟女精品中文字幕| 精品一品国产午夜福利视频| 久久精品国产综合久久久| 制服诱惑二区| 成人午夜精彩视频在线观看| 午夜福利影视在线免费观看| 最黄视频免费看| 国产成人午夜福利电影在线观看| 99久久精品国产国产毛片| 欧美bdsm另类| 久久婷婷青草| 久久青草综合色| 免费观看av网站的网址| 亚洲精品第二区| 最黄视频免费看| 欧美精品高潮呻吟av久久| av又黄又爽大尺度在线免费看| 精品视频人人做人人爽| 成人漫画全彩无遮挡| 看十八女毛片水多多多| 日本欧美国产在线视频|