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

    潛水貫流泵裝置過流部件水力性能分析與優(yōu)化

    2018-04-11 01:55:46夏臣智蔣紅櫻
    農(nóng)業(yè)工程學報 2018年7期
    關(guān)鍵詞:進線過流揚程

    夏臣智,成 立,蔣紅櫻,辛 健

    ?

    潛水貫流泵裝置過流部件水力性能分析與優(yōu)化

    夏臣智1,成 立1,蔣紅櫻2,辛 健3

    (1. 揚州大學水利與能源動力工程學院,揚州 225127; 2. 江蘇省水利工程科技咨詢有限公司,南京 210029; 3. 江蘇亞太泵閥有限公司,泰州 225400)

    為研究潛水貫流泵裝置過流部件的水力性能,該文采用CFD方法對潛水貫流泵裝置進行數(shù)值計算,分析了不同過流部件形式對泵裝置內(nèi)水力性能的影響,并對計算結(jié)果進行試驗驗證。結(jié)果表明:潛水貫流泵裝置燈泡體支撐片的數(shù)量會影響導葉與支撐片之間的水流流態(tài),支撐片的數(shù)量應(yīng)與導葉片數(shù)一致。潛水貫流泵裝置宜采用橢球體的燈泡體尾部形式,能避免回流、脫流等不良流態(tài)的產(chǎn)生。采用流線形進線孔,并且將進線孔與支撐片結(jié)合在一起,能改善出水燈泡體的流態(tài),提高裝置效率。在閘門槽間的進水流道過渡形式宜采用漸縮方管式。優(yōu)化后泵裝置在最優(yōu)工況點的效率提高2.5%,達到78.0%。在最高效率點,數(shù)值計算預(yù)測揚程和流量的不確定度均小于1%,試驗與數(shù)值計算結(jié)果吻合較好。該研究可為潛水貫流泵裝置在實際工程中推廣應(yīng)用提供參考。

    泵;計算機仿真;優(yōu)化;潛水貫流泵;過流部件;模型試驗

    0 引 言

    貫流泵是一種用于低揚程、大流量泵站的經(jīng)濟型泵,具有水力損失小、裝置效率高等特點。南水北調(diào)東線一期工程新建泵站中有7座選用貫流泵[1-4]。潛水貫流泵是一種采用潛水電機作為燈泡體的貫流泵型,與其他貫流泵型相比,具有優(yōu)良的水力性能、緊湊的結(jié)構(gòu)形式、良好的抗噪性和低廉的運行費用等特點,在城市雨水集放、防洪排澇和灌溉方面有良好的應(yīng)用前景[5-6]。

    潛水貫流泵裝置常被用于低揚程、大流量的泵站中,其水泵形式為軸流泵。目前關(guān)于軸流泵水力性能和內(nèi)部流動特性的文獻較為豐富。施衛(wèi)東[7]對高比轉(zhuǎn)速低揚程軸流泵的水力模型進行研究,提出一系列提高效率和抗汽蝕的設(shè)計方法。王福軍等[8]率先采用瞬態(tài)流理論和大渦模擬方法,對軸流泵內(nèi)部非定常流動進行了深入地研究,得到了多工況下軸流泵內(nèi)壓力脈動特性。Qian等[9-10]研究了可調(diào)式導葉對軸流泵的水力性能的影響。鄭源等[11]采用數(shù)值模擬方法和模型試驗對低揚程大流量泵裝置馬鞍區(qū)的流動特性進行研究。張德勝等[12-15]采用大渦模擬方法對模型軸流泵內(nèi)部非定常流動進行了數(shù)值模擬研究,獲得了軸流泵葉頂區(qū)湍流特性,并得到試驗驗證。Li等[16-17]分別采用數(shù)值模擬方法,對軸流泵內(nèi)部空化流動特性進行了研究。Liu等[18]研究了停止工況下軸流泵裝置內(nèi)的瞬態(tài)流動。Wu等[19]基于RBF神經(jīng)網(wǎng)格方法對軸流泵系統(tǒng)的動力模型進行了研究。唐學林等[20]采用流固耦合方法,對燈泡貫流泵葉輪強度進行研究,并與試驗結(jié)果比較。Wang等[21]采用數(shù)值模擬方法對豎井貫流泵裝置內(nèi)的流固耦合特性進行了研究。Kan等[22]采用雙向流固耦合方法,分析了軸伸式貫流泵裝置葉片的壓力動力特性。

    上述文獻對軸流泵水力模型的能量性能、空化特性和葉片壓力特性做了大量的研究工作,為貫流泵裝置的水力模型設(shè)計提供豐富的資料。但是,影響貫流泵裝置水力性能的因素,除了水泵的水力模型之外,還有貫流泵裝置的過流部件。研究表明[23],貫流泵裝置的過流部件對整個裝置的水力性能有較大的影響。目前對于潛水貫流泵裝置過流部件水力性能分析優(yōu)化的研究,未見報道。

    本文采用CFD軟件ANSYS CFX,對設(shè)計工況下潛水貫流泵裝置的過流部件水力性能進行分析研究,并通過模型試驗對計算結(jié)果進行驗證,以期潛水貫流泵裝置在實際工程中推廣應(yīng)用提供參考。

    1 研究對象及模擬方法

    1.1 研究對象

    潛水貫流泵裝置的過流部件包括進水流道、葉輪、導葉、燈泡體和出水流道,其中燈泡體由支撐片、進線孔和燈泡體尾部組成,如圖1所示。葉輪水力模型型號為ZM25,葉輪直徑為300 mm,葉頂間隙為0.15 mm,輪轂比為0.36,葉片數(shù)為3,配套導葉葉片數(shù)為5,轉(zhuǎn)速為1 067 r/min。設(shè)計流量為 270 L/s。

    圖1 潛水貫流泵裝置全流道示意圖

    1.2 數(shù)值模擬

    采用分塊策略對潛水貫流泵裝置進行六面體網(wǎng)格剖分。進、出水流道選用O型拓撲結(jié)構(gòu),便于對邊界層進行加密。借鑒相關(guān)文獻的方法[24],葉輪采用H/J/L-Grid拓撲結(jié)構(gòu),導葉選用H-Grid拓撲結(jié)構(gòu)。為獲得更為精確的結(jié)果,對葉片周邊區(qū)域的網(wǎng)格進行加密處理,整個計算域內(nèi)+值在30~500之間[24-27]。進水流道和葉輪的網(wǎng)格如圖2所示。

    圖2 進水流道和葉輪網(wǎng)格

    為提高計算精度,排除網(wǎng)格數(shù)量對計算結(jié)果的影響,同時適當減少計算時長,對整個計算域進行了網(wǎng)格敏感性分析。網(wǎng)格節(jié)點數(shù)達到1 258 795個后,泵裝置揚程和葉片扭矩波動值均小于1%。

    相關(guān)文獻研究表明[28-29],RNG湍流模型考慮分離流動和渦旋的效應(yīng),適用于泵裝置的數(shù)值計算。采用Scalable wall function處理近壁區(qū)流動。進口條件采用質(zhì)量進口,出口邊界設(shè)置為壓力出口。

    2 數(shù)值模擬方案

    針對潛水貫流泵裝置的結(jié)構(gòu)特點,本文選取3個關(guān)鍵過流部件的5種參數(shù)進行研究,具體為燈泡體中支撐片的數(shù)量和尾部形狀、電纜進線孔斷面形狀和位置,以及進水流道的過渡形狀。由于研究參數(shù)較多,本文采用單因數(shù)控制方法,依次對各過流部件的結(jié)構(gòu)參數(shù)進行研究,首先研究過流部件結(jié)構(gòu)參數(shù)是支撐片的數(shù)量,后續(xù)研究均以前面研究較優(yōu)結(jié)果為基礎(chǔ)。共研究了9種結(jié)構(gòu)方案,其中方案1、2和3研究燈泡體支撐片數(shù)量,分別對應(yīng)為3片、5片和7片;方案2、4和5研究燈泡體尾部形狀,分別對應(yīng)為圓臺體、半圓球體和半橢球體;方案5、6、7和8研究電纜進線孔的型式,分為斷面形狀和位置,其中方案5不設(shè)置進線孔,方案6和方案7的斷面形狀分別為圓形和流線形,進線孔的位置在支撐片之間,方案8的進線孔斷面形狀為流線形,位置在支撐片上;方案8和方案9研究進水流道的過渡形狀,方案3為最初方案。

    3 結(jié)果與分析

    3.1 燈泡體支撐片數(shù)量對泵裝置水力性能影響

    圖3和圖4分別為設(shè)計工況下,支撐片上的靜壓分布和燈泡體的流線圖。由圖3與圖4可知,方案1、2和3在各支撐片間的靜壓分布均勻,順水流方向靜壓逐步降低,未見局部靜壓突變區(qū)域,支撐片之間的流線變化平順,沒有回流、漩渦等不良流態(tài),表明支撐片數(shù)量的增加,對燈泡體內(nèi)的水流流動擾動較小。

    圖3 不同方案支撐片靜壓分布

    圖4 不同方案燈泡體流線圖

    表1為不同支撐片數(shù)方案的水力損失。由表1可知,方案1、2和3的燈泡體水力損失占整個泵裝置水力損失的比重約在30%左右,表明燈泡體的水力性能優(yōu)化對潛水貫流泵裝置能量特性提高有重要影響。隨著支撐片數(shù)量的增加,燈泡體水力損失隨之增大,但支撐片數(shù)量為5時(方案2),泵裝置的水力損失最小。

    表1 不同支撐片數(shù)方案的水力損失

    為了解導葉和支撐片之間水流擴散狀態(tài),對導葉和支撐片之間的軸向速度均勻度[30]變化進行了分析。圖5a為軸向速度分布均勻度的取樣斷面示意圖,圖5b為方案1、2和3的導葉和支撐片之間軸向速度均勻度的變化。如圖5b所示,方案1、2和3的軸向速度均勻度在導葉出口(斷面1)上相近,在之后的斷面上均逐漸增大,其中方案2的軸向速度均勻度在各個斷面上均大于方案1和方案3。這表明,方案2(5片支撐片)中水流在導葉和支撐片之間的擴散要優(yōu)于方案1和方案3。

    圖5 導葉與支撐片間的流動特性

    3.2 燈泡體尾部形式對泵裝置水力性能影響

    圖6為燈泡體尾部縱斷面流態(tài)。由圖6可知,方案2燈泡體尾部有漩渦,由于方案2的尾部形式結(jié)構(gòu)上的特點:底部與燈泡不是平順的連接以及頭部的凸臺狀形式,這兩處水流的穩(wěn)定狀態(tài)被破壞,水流擾動劇烈,產(chǎn)生漩渦。方案4作了相應(yīng)的改進,采用半圓球體的尾部形式,使得水流能夠平順的擴散,得到較好的流態(tài)。由于橢球的斷面形式較圓形更具有流線形的特征,能改變流體在尾部脫離點位置,獲得更好的流態(tài), 方案5采用半橢球體。表3為不同尾部形式方案的水力損失。由表2可知,方案5(半橢球體)的燈泡體水力損失為31 mm,較方案2的減少14 mm,較方案4減少了3 mm,且方案5裝置的水力損失較方案2減少了17 mm,較方案4減少了6 mm。這表明燈泡體尾形式的改變,不僅會影響燈泡體的水力性能,對泵裝置其他部件的性能也有所改善,也與上述流態(tài)的結(jié)果相吻合。

    圖6 燈泡體尾部縱斷面流態(tài)

    表2 不同尾部形式方案的水力損失

    綜上,潛水貫流泵裝置宜采用半橢球體的燈泡體尾部形式。

    3.3 電纜進線孔斷面形狀和位置對泵裝置水力性能影響

    圖7和圖8分別為方案5、6、7和8的燈泡體表面和燈泡體支撐片末端斷面靜壓分布。表4為上述方案燈泡段和裝置的水力損失。由圖7和圖8可知,無電纜進線孔(方案5),燈泡體表面靜壓沿水流流動方向逐漸減小,沒有局部靜壓突變的區(qū)域,各個支撐片以及支撐片之間的燈泡體表面靜壓變化一致。支撐片末端斷面上,5個扇形面上靜壓分布也是近乎一致的。在2個支撐片之間設(shè)置電纜進線孔時(方案6),燈泡體表面的靜壓分布與方案5不同,在進線孔兩側(cè)形成局部低壓區(qū),所在位置的支撐片前端靜壓較其他支撐片高。在支撐片末端斷面上,進線孔后的扇形面上存在局部高壓區(qū),但其他扇形面的靜壓分布與方案5近乎一致。由表3可知,方案5的水力損失為31 mm,方案6的水力損失為36 mm,相對與方案5增大了16%。表明進線孔會影響燈泡體的水流流態(tài),增大燈泡體和泵裝置的水力損失。

    圖7 燈泡體表面靜壓分布

    圖8 支撐片末端斷面靜壓分布

    表3 不同進線孔型式方案的水力損失

    進線孔橫斷面形狀由圓形(方案6)變?yōu)榱骶€形(方案7)時,支撐片末端斷面上的局部高壓區(qū)變小,其他區(qū)域的靜壓分布變化較小,方案7燈泡體的水力損失為35 mm,與方案6相比,僅僅減少了1 mm,表明與圓形進線孔相比,流線形進線孔略微改善進線孔后方的水流流態(tài),但對燈泡體和泵裝置的水力損失影響較小。

    進線孔的位置由支撐片之間(方案7)改成支撐片之上(方案8),對燈泡體表面的靜壓分布有明顯的改善。方案7的進線孔前后均有局部高壓區(qū),兩側(cè)有局部低壓區(qū),水流流動變化急劇。與此相比,方案8進線孔周圍的靜壓變化較為平順,支撐片末端斷面上未見局部高壓區(qū),減少了進線孔由于阻水分流和尾跡帶來流場擾動,因而水力損失較?。ㄒ姳?)。

    潛水貫流泵裝置的進線孔宜采用流線形的斷面,并與支撐片結(jié)合在一起。

    3.4 閘門槽間進水流道過渡形狀對泵裝置水力性能影響

    圖9方案8和方案9進水流道近壁面流線圖。表5為方案8和方案9的進水流道水力損失和進水流道出口流速分布均勻度。由圖9可知,方案8進水流道近壁面流線較為平順,方案9進水流道近壁面流線在第二道門槽上流線紊亂。由表4可知,方案8的水力損失較小,流速分布均勻度較高,能提供較好的水泵入流條件。表明在兩道閘門槽之間的過渡形式宜采用漸縮方管,既能避免因閘門槽部位過流斷面急劇變化而造成的不良流態(tài),又能方便工程施工。

    圖9 進水流道流態(tài)

    表4 不同進水流道過渡形狀方案的水力損失及進水流道出口流速分布均勻度

    4 最優(yōu)方案與優(yōu)化試驗驗證

    4.1 過流部件優(yōu)化前后水力性能比較

    方案3為最初方案,方案8為優(yōu)化方案。與方案3相比,方案8對3種過流部件的5種結(jié)構(gòu)參數(shù)進行了優(yōu)化,優(yōu)化后的過流部件參數(shù)包括支撐片數(shù)量(5片)、燈泡體尾部形狀(半橢球體)、進線孔斷面形狀(流線形)和位置(支撐片上),以及進水流道的過渡形式(漸縮方管)。

    圖10為方案8和方案3的泵裝置能量性能對比。由圖可知,當流量系數(shù)/BEP為1時(為泵裝置流量,L/s;BEP為試驗的泵裝置最高效率點對應(yīng)的流量,L/s),方案3的揚程系數(shù)/BEP為0.89(為泵裝置揚程,m;BEP為試驗的泵裝置最高效率點對應(yīng)的流量,m),效率為75.5%;方案8的揚程系數(shù)/BEP為1.01,效率為78.0%。方案8較方案3在最優(yōu)點效率提高了2.5%。

    注:Q為泵裝置流量,L·s-1;QBEP為試驗最高效率點對應(yīng)的流量,L·s-1;H為泵裝置揚程,m;HBEP為試驗最高效率點對應(yīng)的揚程,m。

    4.2 優(yōu)化試驗驗證

    為驗證潛水貫流泵裝置數(shù)值計算結(jié)果的可靠性,采用物理模型進行能量性能測試。物理模型的葉輪直徑為300 mm,葉頂間隙為0.15 mm,輪轂比為0.36,葉片數(shù)為3,配套導葉葉片數(shù)為5,試驗轉(zhuǎn)速為1 067 r/min。物模試驗在經(jīng)過中國計量認證的水力試驗臺上進行,試驗參數(shù)測量參照SL140-2006《水泵模型及裝置模型驗收試驗規(guī)程》。潛水貫流泵裝置物理模型如圖11所示。試驗結(jié)果與計算結(jié)果比較如圖12所示。試驗和計算得到的流量-揚程曲線在計算最高效率點(/BEP=1)處交叉,在大流量區(qū)(/BEP>1.05),計算預(yù)測的揚程大于試驗值,在小流量區(qū)(/BEP<1)計算預(yù)測的揚程小于試驗值。在最高效率點,揚程和流量預(yù)測的不確定度在1%以內(nèi)。流量系數(shù)/BEP在0.98~1.16范圍內(nèi),計算預(yù)測的效率在70%以上,與試驗結(jié)果吻合較好。數(shù)值計算結(jié)果是可信的。

    圖11 潛水貫流泵裝置物理模型

    圖12 泵裝置外特性曲線對比

    5 結(jié) 論

    本文比較了不同支撐片數(shù)量對泵裝置水力性能的影響,結(jié)果表明,隨著支撐片數(shù)量的增加,燈泡體的水力損失隨之增大,但當支撐片數(shù)目與導葉一致時,整個泵裝置的水力損失最小。

    比較了圓臺體、圓球體和橢圓體3種尾部形式對泵裝置水力性能影響。采用半橢球體的燈泡體尾部形式,能避免回流、脫流等不良流態(tài)的產(chǎn)生,減小泵裝置的水力損失。潛水貫流泵裝置電纜進線孔的形狀和位置對該處水流態(tài)有明顯的影響。采用流線形的進線孔,并且將進線孔與支撐片結(jié)合在一起,有助于改善出水燈泡體的流態(tài),提高裝置效率。

    潛水貫流泵裝置防洪事故門槽和檢修門槽之間的進水流道過渡形式,對進水流道內(nèi)的流態(tài)有明顯影響。采用漸縮方管的形式,能獲得較好的進水流態(tài),并且也方便工程施工。優(yōu)化后潛水貫流泵裝置在最優(yōu)工況點的效率提高2.5%,達到78.0%。

    采用模型試驗對優(yōu)化后的潛水貫流泵裝置的水力性能進行驗證,在最高效率點,揚程和流量預(yù)測的不確定度小于1%,試驗與數(shù)值計算結(jié)果吻合較好。

    [1] 張仁田,單海春,卜舸,等. 南水北調(diào)東線一期工程燈泡貫流泵結(jié)構(gòu)特點[J]. 排灌機械工程學報,2016,34(9):774-782,789. Zhang Rentian, Shan Haichun, Bu Ge, et al. Structural features of bulb tubular pumps in first phase of South-to- NorthWater Diversion Eastern Route Project in China[J]. Journal of Drainage and Irrigation Machinery Engineering, 2016, 34(9): 774-782, 789 (in Chinese with English abstract)

    [2] 徐磊,陸林廣,陳偉,等. 南水北調(diào)工程邳州站豎井貫流泵裝置進出水流態(tài)分析[J]. 農(nóng)業(yè)工程學報,2012,28(6):50-56. Xu Lei, Lu Linguang, Chen Wei, et al. Flow pattern analysis on inlet and outlet conduit of shaft tubular pump system of Pizhou pumping station in South-to-North Water Diversion Project[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(6): 50-56. (in Chinese with English abstract)

    [3] 金燕. 貫流泵內(nèi)部流動的數(shù)值模擬與三維LDV測量研究[D]. 揚州:揚州大學,2010. Jin Yan. Study on the Numerical Simulation and 3D-LDV Measurement of Tubular Pump[D]. Yangzhou: Yangzhou University, 2010. (in Chinese with English abstract)

    [4] 馮旭松,關(guān)醒凡,井書光,等. 南水北調(diào)東線燈泡貫流泵水力模型及裝置研究開發(fā)與應(yīng)用[J]. 南水北調(diào)與水利科技,2009,7(6):32-35.Feng Xusong, Guan Xingfan, Jing Shuguang, et al. Development and application on hydraulic model and equipment of bulb tubular pumps in the eastern route of the south-to-north water transfer project[J]. South-to-North Water Transfers and Water Science & Technology, 2009, 7(6): 32-35. (in Chinese with English abstract)

    [5] 楊帆,金燕,劉超,等. 雙向潛水貫流泵裝置性能試驗與數(shù)值分析[J]. 農(nóng)業(yè)工程學報,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 Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(16): 60-67. (in Chinese with English abstract)

    [6] Dijksman J F. Hydrodynamics of small tubular pumps [J]. Journal of Fluid Mechanics, 1984, 139: 173-191.

    [7] 施衛(wèi)東. ZM931高比轉(zhuǎn)數(shù)軸流泵水力模型的設(shè)計[J]. 農(nóng)業(yè)機械學報,1998,29(2):49-53. Shi Weidong. Design of Axial Flow Pump Hydraulic model Zm931 on high specific speed [J]. Transactions of the Chinese Society for Agricultural Machinery, 1998, 29(2): 49-53. (in Chinese with English abstract)

    [8] 王福軍,張玲,張志民. 軸流泵不穩(wěn)定流場的壓力脈動特性研究[J]. 水利學報,2007,38(8):1003-1009. Wang Fujun, Zhang Ling, Zhang Zhimin. Analysis on pressure fluctuation of unsteady flow in axial-flow pump [J]. Journal of Hydraulic Engineering, 2007, 38(8): 1003-1009. (in Chinese with English abstract)

    [9] Qian Z, Wang F, Guo Z, et al. Performance evaluation of an axial-flow pump with adjustable guide vanes in turbine mode[J]. Renewable Energy, 2016, 99: 1146-1152.

    [10] 錢忠東,王凡,王志遠,等. 可調(diào)導葉式軸流泵馬鞍區(qū)水力特性試驗研究[J]. 排灌機械工程學報,2013,31(6):461-465. Qian Zhongdong, Wang Fan, Wang Zhiyuan, et al. Experimental study on hydraulic performance of saddle zone in axial flow pump with adjustable guide vane[J]. Journal of Drainage and Irrigation Machinery Engineering, 2013, 31(6): 461-465. (in Chinese with English abstract)

    [11] 鄭源,茅媛婷,周大慶,等. 低揚程大流量泵裝置馬鞍區(qū)的流動特性[J]. 排灌機械工程學報,2011,29(5):369-373. Zheng Yuan, Mao Yuanting, Zhou Daqing, et al. Flow characteristics of low-lift and large flow rate pump installationin saddle zone[J]. Journal of Drainage and Irrigation Machinery Engineering, 2011, 29(5): 369-373. (in Chinese with English abstract)

    [12] 張德勝,石磊,陳健,等. 基于大渦模擬的軸流泵葉頂泄漏渦瞬態(tài)特性分析[J]. 農(nóng)業(yè)工程學報,2015,31(11):74-80. Zhang Desheng, Shi Lei, Chen Jian, et al. Analysis on transient characteristics of tip leakage vortex in axial flow pump using large eddy simulation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(11): 74-80. (in Chinese with English abstract)

    [13] Zhang D, Shi L, Zhao R, et al. Study on unsteady tip leakage vortex cavitation in an axial-flow pump using an improved filter-based model[J]. Journal of Mechanical Science & Technology, 2017, 31(2): 659-667.

    [14] Zhang D, Shi W, Pan D, et al. Numerical and experimental investigation of tip leakage vortex cavitation patterns and mechanisms in an axial flow pump[J]. Journal of Fluids Engineering, 2015, 137(12): 815-816.

    [15] Zhang D, Shi W, Esch B P M V, et al. Numerical and experimental investigation of tip leakage vortex trajectory and dynamics in an axial flow pump[J]. Computers & Fluids, 2013, 112(1): 61-71.

    [16] Li Z, Zhang N, Hong B, et al. Simulation research on cavitation flow in tip clearance of axial-flow pump[J]. Journal of Engineering Thermophysics, 2011, 588(8): 1255-1258.

    [17] Feng W, Cheng Q, Guo Z, et al. Simulation of cavitation performance of an axial flow pump with inlet guide vanes[J]. Advances in Mechanical Engineering, 2016, 8(6): 1-8.

    [18] Liu Y, Zhou J, Zhou D. Transient flow analysis in axial-flow pump system during stoppage[J]. Advances in Mechanical Engineering, 2017, 9(9): 1-8.

    [19] Wu Q, Wang X, Shen Q. Research on dynamic modeling and simulation of axial-flow pumping system based on RBF neural network[J]. Neurocomputing, 2016, 186: 200-206.

    [20] 唐學林,王秀葉,賈玉霞. 基于流固耦合的燈泡貫流泵葉輪強度分析[J]. 排灌機械工程學報,2014,32(11):921-926. Tang Xuelin, Wang Xiuye, Jia Yuxia. Strength analysis of bulb tubular pump impeller based on fluid-structure interaction[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 32(11): 921-926. (in Chinese with English abstract)

    [21] Wang S, Zhang L, Yin G. Numerical investigation of the fsi characteristics in a tubular pump[J]. Mathematical Problems in Engineering, 2017(17): 1-9.

    [22] Kan K, Zheng Y, Fu S, et al. Dynamic stress of impeller blade of shaft extension tubular pump device based on bidirectional fluid-structure interaction[J]. Journal of Mechanical Science & Technology, 2017, 31(4): 1561-1568.

    [23] 仇寶云,裴蓓,申劍,等. 燈泡貫流泵機組支撐形式比較[J]. 排灌機械工程學報,2011,29(1):61-66. Qiu Baoyun, Pei Bei, Shen Jian, et al. Camparison of bulb tubular pump unit’s support forms[J]. Journal of Drainage and Irrigation Machinery Engineering, 2011, 29(1): 61-66. (in Chinese with English abstract)

    [24] 楊帆. 低揚程泵裝置水動力特性及多目標優(yōu)化關(guān)鍵技術(shù)研究[D]. 揚州:揚州大學,2013. Yang Fan. Research on Hydraulic Performance and Multi- objective optimization Design of Low-Lift Pump system [D]. Yangzhou : Yangzhou University, 2013.

    [25] 資丹,王福軍,陶然,等. 邊界層網(wǎng)格尺度對泵站流場計算結(jié)果影響研究[J]. 水利學報,2016,47(2):139-149. Zi Dan, Wang Fujun, Tao Ran, et al. Research for impacts of boundary layer grid scale on flow field simulation results in pumping station[J]. Journal of Hydraulic Engineering, 2016, 47(2): 139-149. (in Chinese with English abstract)

    [26] Salim S M, Ariff M, Cheah S C. Wall+ approach for dealing with turbulent flows over a wall mounted cube[J]. Progress in Computational Fluid Dynamics An International Journal, 2010, 10(8): 341-351.

    [27] Tu J, Yeoh G H, Liu C. Computational fluid dynamics: A practical approach[J]. Artificial Organs, 2013, 33(9): 727-32.

    [28] Yang G. Pump Turbine Working Condition of Numerical Simulation and Optimization of the RNGModel[M]. Intelligent Computing Methodologies. Springer International Publishing, 2016: 684-693.

    [29] Mostafa N H, Boraey M A. Numerical and experimental investigation of cavitation in axial pumps[J]. International Journal of Manufacturing, 2013, 4(1): 1225-1225.

    [30] 劉超. 水泵及水泵站[M].北京:中國水利水電出版社,2009.

    Hydraulic performance analysis and optimization on flow passage components of diving tubular pumping system

    Xia Chenzhi1, Cheng Li1, Jiang Hongying2, Xin Jian3

    (1.225127; 2.210029; 3.225400,)

    Diving tubular pump is a kind of low-head large-flow horizontal pump. This type is postpositive bulb pump and the motor is coaxial with pump, which is suitable for irrigation or storm drainage pumping station. For the advantages of excellent hydraulic performance, compact structure, good noise resistance and low operating cost, diving tubular pump is widely applied in low head pumping station. The internal structure of diving tubular pumping system is more complex than other pumping systems for the combination of motor and rotor. Bulb, support slices, cable hole, inlet passage and outlet passage constitute the flow passage components of diving tubular pumping system, which have influence on the internal flow pattern and hydraulic loss of diving tubular pumping system. For low head condition, small increase of hydraulic loss will result in large efficiency loss of pumping system. So the analysis on hydraulic performance of flow passage components of diving tubular pumping system needs to be carried out. In this paper, to acquire hydraulic performance of flow passage components, the inner flow field of diving tubular pumping system is calculated by CFD (computational fluid dynamics). The influence on hydraulic performance of pumping system with different flow passage components is analyzed and the experiment is performed to verify the result of numerical calculation. The head and efficiency difference between model test and CFD simulation are less than 1% in the best efficiency point condition, and the high efficiency range predicted by CFD coincides with the results of the experiment. In order to acquire the influence on hydraulic performance for different numbers of support slices, 3 kinds of slices’ numbers are simulated to show that the number of support slices affects the flow pattern between guide vane and support slice. The hydraulic loss of bulb part increases with the number of support slices. The hydraulic loss of pumping system is relatively lower while the support slices’ number is equal to guide vane’s, which is 5 in this pumping system. Also the shape of bulb tail of pumping system has influence on the flow pattern in outlet passage. Three kinds of bulb tail shapes are calculated, including cone, semi-ellipsoid and semi-sphere. Cone is bad shape for bulb tail, due to the bad flow pattern caused by blunt tail. Semi-sphere is better than cone for less hydraulic loss and better flow expansion in outlet passage. Compared to cone and semi-sphere, semi-ellipsoid is the best shape for bulb tail, which can prevent bad flow pattern like backflow and flow separation. Cable hole is a small pipe which locates at bulb, and its influence on flow pattern and hydraulic loss is ignored by manufactures for its small volume. This paper compares 2 different section shapes of cable hole, specific circle and streamline shape. The result shows that cable hole causes 16% hydraulic loss of bulb, and streamline shape acquires better pressure distribution on cable hole and slightly enhances the flow pattern behind cable hole. Besides, the combination of cable hole and support slice can improve the flow pattern in bulb part and improve efficiency of pumping system. The transition shape of inlet passage between 2 gate slots is also a key factor on hydraulic performance of diving tubular pumping system. Two transition shapes are calculated, including converging square pipe and converging square-circle pipe. Compared with converging square-circle pipe, converging square pipe decreases hydraulic loss of pumping system by 5% and acquires better flow pattern around gate slots. In addition, the other advantage of converging square pipe is the convenience for construction. In conclusion, the efficiency of optimized diving tubular pumping system increases by 2.5%, which reaches 78%. In best efficiency point, the uncertainty of hydraulic performance between simulation and test result less than 1%.

    pumps; computer simulation; optimization; diving tubular pump; flow passage components; model test

    夏臣智,成 立,蔣紅櫻,辛 健.潛水貫流泵裝置過流部件水力性能分析與優(yōu)化[J]. 農(nóng)業(yè)工程學報,2018,34(7):45-51. doi:10.11975/j.issn.1002-6819.2018.07.006 http://www.tcsae.org

    Xia Chenzhi, Cheng Li, Jiang Hongying, Xin Jian. Hydraulic performance analysis and optimization on flow passage components of diving tubular pumping system[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(7): 45-51. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2018.07.006 http://www.tcsae.org

    2017-10-22

    2018-01-03

    國家自然科學基金面上(51779214);江蘇省高校優(yōu)勢學科建設(shè)工程資助項目(PAPD);江蘇省研究生培養(yǎng)創(chuàng)新工程(KYLX16_1395);江蘇水利科技項目(2016035)

    夏臣智,博士生,主要從事泵站工程研究。Email:xiachenzhi@foxmail.com

    成 立,教授,博士生導師,主要從事泵站工程研究。Email:chengli@yzu.edu.cn

    10.11975/j.issn.1002-6819.2018.07.006

    TH31;TV131

    A

    1002-6819(2018)-07-0045-07

    猜你喜歡
    進線過流揚程
    管路受力誘發(fā)高揚程離心泵振動加劇原因分析
    變壓器相間過流保護不正確動作原因的探討
    電子制作(2017年19期)2017-02-02 07:08:58
    小型PLC在變配電站電源自動切換控制中的應(yīng)用
    國產(chǎn)500kW短波發(fā)射機過流保護電路的應(yīng)用缺陷與改進
    電子制作(2016年11期)2016-11-07 08:43:36
    水輪機過流部件改造與節(jié)能增效
    配網(wǎng)雙端閉鎖式過流保護的實現(xiàn)方法
    新型中高揚程大流量水錘泵結(jié)構(gòu)技術(shù)改進研究
    中國水利(2015年7期)2015-02-28 15:12:58
    高揚程長距離輸水管道系統(tǒng)水錘防護的模擬分析
    高速泵與旋殼泵在小流量高揚程情況下的適用性分析
    河南科技(2014年4期)2014-02-27 14:07:13
    基于PLC實現(xiàn)的低壓三進線 三母聯(lián)備自投系統(tǒng)
    国产精品一区二区三区四区免费观看 | 欧美日韩福利视频一区二区| 国产又黄又爽又无遮挡在线| 亚洲av五月六月丁香网| 99re在线观看精品视频| 亚洲九九香蕉| 一本一本综合久久| 看免费av毛片| 999久久久国产精品视频| 黄片小视频在线播放| 亚洲aⅴ乱码一区二区在线播放| 无人区码免费观看不卡| 亚洲av成人一区二区三| 999久久久国产精品视频| 亚洲国产欧洲综合997久久,| 香蕉国产在线看| www.精华液| 99国产精品一区二区三区| 国产乱人视频| 亚洲av美国av| 免费人成视频x8x8入口观看| 日本在线视频免费播放| 久久香蕉精品热| 99re在线观看精品视频| 亚洲精华国产精华精| 精品日产1卡2卡| 国产亚洲精品综合一区在线观看| 国产av在哪里看| 视频区欧美日本亚洲| 好男人在线观看高清免费视频| 成人18禁在线播放| 国产精品久久久久久久电影 | 亚洲国产精品成人综合色| 欧美3d第一页| 不卡av一区二区三区| av中文乱码字幕在线| 琪琪午夜伦伦电影理论片6080| 99在线视频只有这里精品首页| 淫妇啪啪啪对白视频| 久久精品91蜜桃| 精品国产美女av久久久久小说| 日日干狠狠操夜夜爽| 这个男人来自地球电影免费观看| 欧美在线一区亚洲| 黄色视频,在线免费观看| a级毛片a级免费在线| 精品久久久久久久毛片微露脸| 亚洲aⅴ乱码一区二区在线播放| 亚洲成人久久爱视频| 色综合欧美亚洲国产小说| 中文在线观看免费www的网站| 啦啦啦韩国在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看| 日本成人三级电影网站| 日韩欧美国产一区二区入口| 亚洲自拍偷在线| 国产精品久久久av美女十八| 国产精品,欧美在线| 亚洲aⅴ乱码一区二区在线播放| 欧美成人免费av一区二区三区| 超碰成人久久| av国产免费在线观看| 国产成人精品无人区| 老司机福利观看| 99国产精品一区二区三区| 18禁黄网站禁片免费观看直播| 亚洲九九香蕉| 日本在线视频免费播放| 国产一级毛片七仙女欲春2| 99精品欧美一区二区三区四区| 国产亚洲精品久久久com| xxx96com| 国内少妇人妻偷人精品xxx网站 | 国产成+人综合+亚洲专区| 久久久国产欧美日韩av| 久久亚洲真实| 亚洲av成人精品一区久久| 欧美性猛交黑人性爽| 美女cb高潮喷水在线观看 | 国产精品九九99| 别揉我奶头~嗯~啊~动态视频| 免费在线观看日本一区| 婷婷精品国产亚洲av在线| 国产精品自产拍在线观看55亚洲| 中文字幕精品亚洲无线码一区| 麻豆一二三区av精品| 久久久久久九九精品二区国产| 99在线人妻在线中文字幕| 精品一区二区三区视频在线 | 亚洲成人久久性| 亚洲精品美女久久av网站| 久99久视频精品免费| 小蜜桃在线观看免费完整版高清| 亚洲乱码一区二区免费版| 日韩欧美免费精品| 色av中文字幕| 熟女电影av网| 久久香蕉国产精品| 女人被狂操c到高潮| www日本在线高清视频| 91av网一区二区| 日本免费a在线| 狂野欧美激情性xxxx| 亚洲国产精品合色在线| 黄色女人牲交| 高清毛片免费观看视频网站| 中文字幕久久专区| 天天躁狠狠躁夜夜躁狠狠躁| 搡老妇女老女人老熟妇| 国产综合懂色| 我的老师免费观看完整版| 女生性感内裤真人,穿戴方法视频| bbb黄色大片| 中文资源天堂在线| 黄色女人牲交| 国产成人av教育| 国产熟女xx| 美女免费视频网站| 男插女下体视频免费在线播放| 亚洲在线观看片| 每晚都被弄得嗷嗷叫到高潮| 日本免费a在线| 国产成人啪精品午夜网站| 熟女电影av网| 国产精品精品国产色婷婷| 黄片大片在线免费观看| 亚洲国产日韩欧美精品在线观看 | 日韩欧美免费精品| 久久国产精品人妻蜜桃| 免费在线观看影片大全网站| 18禁美女被吸乳视频| 成人18禁在线播放| 老汉色∧v一级毛片| 国产av不卡久久| 国产精品美女特级片免费视频播放器 | 日本 av在线| 桃红色精品国产亚洲av| 国产熟女xx| 男女下面进入的视频免费午夜| 嫩草影院入口| 小蜜桃在线观看免费完整版高清| 少妇人妻一区二区三区视频| 日本一二三区视频观看| 欧美乱码精品一区二区三区| 午夜免费观看网址| 桃色一区二区三区在线观看| 久久久久国产精品人妻aⅴ院| 丰满的人妻完整版| 亚洲色图 男人天堂 中文字幕| 在线国产一区二区在线| 日韩欧美一区二区三区在线观看| 美女 人体艺术 gogo| 亚洲午夜理论影院| 在线永久观看黄色视频| 可以在线观看的亚洲视频| 国产精品亚洲美女久久久| 久久精品国产亚洲av香蕉五月| 日本一二三区视频观看| 免费在线观看视频国产中文字幕亚洲| 亚洲熟妇中文字幕五十中出| 狂野欧美激情性xxxx| 在线a可以看的网站| 亚洲国产欧洲综合997久久,| 色综合欧美亚洲国产小说| 全区人妻精品视频| 综合色av麻豆| 日韩欧美国产在线观看| 亚洲 国产 在线| 搡老妇女老女人老熟妇| 制服人妻中文乱码| www.自偷自拍.com| 夜夜爽天天搞| 国产成人精品久久二区二区免费| 91麻豆精品激情在线观看国产| 午夜精品久久久久久毛片777| 欧美乱码精品一区二区三区| 精华霜和精华液先用哪个| 日本免费a在线| av女优亚洲男人天堂 | 国产亚洲av高清不卡| 国产乱人视频| 久久精品国产清高在天天线| 天天躁日日操中文字幕| 狂野欧美激情性xxxx| 黑人操中国人逼视频| 久久精品人妻少妇| www.999成人在线观看| 亚洲精品久久国产高清桃花| 欧美中文综合在线视频| 亚洲最大成人中文| 国产av在哪里看| 听说在线观看完整版免费高清| 在线观看日韩欧美| 国产探花在线观看一区二区| 免费在线观看影片大全网站| 长腿黑丝高跟| 国产久久久一区二区三区| 精品国产超薄肉色丝袜足j| a级毛片a级免费在线| 午夜成年电影在线免费观看| av天堂在线播放| 在线观看午夜福利视频| 搞女人的毛片| 亚洲精品美女久久av网站| 久久草成人影院| 国产99白浆流出| 香蕉丝袜av| 国产 一区 欧美 日韩| 人人妻人人澡欧美一区二区| 精品一区二区三区视频在线 | 日韩成人在线观看一区二区三区| 熟女人妻精品中文字幕| 亚洲熟女毛片儿| 极品教师在线免费播放| 99精品在免费线老司机午夜| av欧美777| 可以在线观看的亚洲视频| 亚洲午夜精品一区,二区,三区| 亚洲欧美日韩无卡精品| 成人av在线播放网站| 国产伦精品一区二区三区四那| 亚洲成人免费电影在线观看| 欧美激情在线99| 精品日产1卡2卡| 一个人看的www免费观看视频| 国产三级黄色录像| 最新美女视频免费是黄的| 午夜福利18| 男人和女人高潮做爰伦理| 在线十欧美十亚洲十日本专区| 国产欧美日韩精品亚洲av| 色综合欧美亚洲国产小说| 国产高潮美女av| 99久久无色码亚洲精品果冻| 国产精品99久久久久久久久| 午夜福利免费观看在线| 中文字幕久久专区| 美女黄网站色视频| 日韩欧美在线二视频| 久久久久亚洲av毛片大全| 日韩高清综合在线| 日韩精品中文字幕看吧| 免费看美女性在线毛片视频| 亚洲av中文字字幕乱码综合| 亚洲,欧美精品.| 欧美成人免费av一区二区三区| 久久久久久人人人人人| 又大又爽又粗| 国产黄片美女视频| 亚洲av第一区精品v没综合| 啦啦啦韩国在线观看视频| 久久这里只有精品中国| 中出人妻视频一区二区| 男女午夜视频在线观看| 无人区码免费观看不卡| 99久久久亚洲精品蜜臀av| 嫩草影院入口| 日韩精品中文字幕看吧| 午夜亚洲福利在线播放| 9191精品国产免费久久| АⅤ资源中文在线天堂| 午夜久久久久精精品| 午夜福利18| 99精品久久久久人妻精品| 成年版毛片免费区| 一夜夜www| 国产伦精品一区二区三区四那| 日韩欧美一区二区三区在线观看| 熟女少妇亚洲综合色aaa.| 久久久久亚洲av毛片大全| 欧美xxxx黑人xx丫x性爽| 香蕉久久夜色| АⅤ资源中文在线天堂| 淫妇啪啪啪对白视频| 日本 av在线| 最新在线观看一区二区三区| 欧美+亚洲+日韩+国产| 国产精品久久视频播放| 日韩国内少妇激情av| 国产午夜精品久久久久久| 国产成人啪精品午夜网站| 午夜福利高清视频| 男插女下体视频免费在线播放| 亚洲中文日韩欧美视频| 国产一区二区激情短视频| 最新美女视频免费是黄的| 亚洲国产看品久久| 欧美成人免费av一区二区三区| 久久欧美精品欧美久久欧美| 女警被强在线播放| 一区二区三区高清视频在线| 国产欧美日韩一区二区精品| 国产极品精品免费视频能看的| 国产成人精品无人区| 亚洲国产日韩欧美精品在线观看 | 99在线视频只有这里精品首页| 国产精品98久久久久久宅男小说| 亚洲黑人精品在线| 不卡av一区二区三区| 国产精品爽爽va在线观看网站| 亚洲av成人精品一区久久| 久久精品aⅴ一区二区三区四区| 日本撒尿小便嘘嘘汇集6| 三级国产精品欧美在线观看 | 国产1区2区3区精品| 婷婷六月久久综合丁香| 特大巨黑吊av在线直播| 18禁黄网站禁片免费观看直播| 精品一区二区三区四区五区乱码| 成人永久免费在线观看视频| 97超视频在线观看视频| 国产人伦9x9x在线观看| 99国产极品粉嫩在线观看| 性色av乱码一区二区三区2| 国产伦精品一区二区三区四那| 欧美黄色片欧美黄色片| 久久国产精品人妻蜜桃| 黄片小视频在线播放| 国产精品 欧美亚洲| 村上凉子中文字幕在线| 欧美丝袜亚洲另类 | 一边摸一边抽搐一进一小说| 成人高潮视频无遮挡免费网站| 国产亚洲欧美在线一区二区| 国产激情偷乱视频一区二区| 欧美黄色片欧美黄色片| 国产在线精品亚洲第一网站| 丰满人妻一区二区三区视频av | 两个人看的免费小视频| 国产精品久久久久久精品电影| 成人亚洲精品av一区二区| 色在线成人网| 国产精品亚洲av一区麻豆| 国产精品久久久久久亚洲av鲁大| 狂野欧美白嫩少妇大欣赏| 一级黄色大片毛片| 色综合站精品国产| 首页视频小说图片口味搜索| 少妇的逼水好多| 亚洲成人久久性| 久久中文字幕人妻熟女| 91麻豆av在线| 亚洲 欧美 日韩 在线 免费| 欧美性猛交╳xxx乱大交人| 丰满人妻熟妇乱又伦精品不卡| 999久久久国产精品视频| 搡老妇女老女人老熟妇| 真实男女啪啪啪动态图| 久久久国产欧美日韩av| 俺也久久电影网| 欧美高清成人免费视频www| 国产精品久久久人人做人人爽| 狂野欧美白嫩少妇大欣赏| 国产精品一区二区精品视频观看| 国产亚洲欧美在线一区二区| 午夜成年电影在线免费观看| 男人舔女人的私密视频| 国产精品 国内视频| 亚洲av中文字字幕乱码综合| 国内揄拍国产精品人妻在线| 可以在线观看的亚洲视频| 亚洲国产精品合色在线| 搡老妇女老女人老熟妇| 日韩精品中文字幕看吧| 男女视频在线观看网站免费| 日本成人三级电影网站| 禁无遮挡网站| 999精品在线视频| 真人一进一出gif抽搐免费| 精品国产亚洲在线| 欧美成狂野欧美在线观看| 国产亚洲av高清不卡| 欧美成人一区二区免费高清观看 | 久久性视频一级片| 亚洲一区二区三区色噜噜| 男人的好看免费观看在线视频| 女同久久另类99精品国产91| 禁无遮挡网站| 国产高潮美女av| 精品熟女少妇八av免费久了| 国产精品久久久久久亚洲av鲁大| 国产97色在线日韩免费| 欧美在线黄色| 国产成人欧美在线观看| 国产不卡一卡二| 亚洲avbb在线观看| 精品国产美女av久久久久小说| 精品久久久久久久久久久久久| 国产成人欧美在线观看| 少妇人妻一区二区三区视频| 美女扒开内裤让男人捅视频| 日韩欧美一区二区三区在线观看| 真实男女啪啪啪动态图| 亚洲熟妇熟女久久| 精品久久久久久久人妻蜜臀av| 女生性感内裤真人,穿戴方法视频| 丁香六月欧美| av在线天堂中文字幕| netflix在线观看网站| 国产精品爽爽va在线观看网站| 亚洲av成人一区二区三| 男女之事视频高清在线观看| 国语自产精品视频在线第100页| 嫩草影院入口| 午夜福利在线观看免费完整高清在 | 搡老岳熟女国产| 免费看日本二区| 国产不卡一卡二| www.精华液| 草草在线视频免费看| 黄色丝袜av网址大全| 在线永久观看黄色视频| 欧美一级a爱片免费观看看| 波多野结衣高清作品| 精品国内亚洲2022精品成人| 九九热线精品视视频播放| 韩国av一区二区三区四区| 性色av乱码一区二区三区2| 国内精品一区二区在线观看| 18美女黄网站色大片免费观看| 亚洲熟女毛片儿| 久久久色成人| 欧美一区二区精品小视频在线| 午夜激情欧美在线| 村上凉子中文字幕在线| 亚洲国产中文字幕在线视频| 1024香蕉在线观看| 日日夜夜操网爽| 一个人免费在线观看的高清视频| 免费在线观看视频国产中文字幕亚洲| 母亲3免费完整高清在线观看| 国产视频一区二区在线看| 午夜精品在线福利| 亚洲成人久久性| 非洲黑人性xxxx精品又粗又长| 日本与韩国留学比较| 青草久久国产| 国产一级毛片七仙女欲春2| 成人性生交大片免费视频hd| 人人妻人人澡欧美一区二区| 校园春色视频在线观看| 午夜免费观看网址| 国产欧美日韩精品一区二区| 在线观看美女被高潮喷水网站 | 成人鲁丝片一二三区免费| 午夜福利欧美成人| 国产又黄又爽又无遮挡在线| 亚洲中文av在线| 国产亚洲精品av在线| 欧美日韩乱码在线| 国产黄片美女视频| 亚洲欧美日韩高清专用| 午夜免费成人在线视频| 亚洲国产中文字幕在线视频| 国产高清视频在线观看网站| 亚洲av片天天在线观看| 国产成人精品久久二区二区91| 国产精品99久久99久久久不卡| 午夜免费激情av| 国产爱豆传媒在线观看| 狂野欧美激情性xxxx| 亚洲精品一区av在线观看| 国产美女午夜福利| 日本与韩国留学比较| 欧美另类亚洲清纯唯美| 波多野结衣高清作品| 亚洲av成人一区二区三| 夜夜爽天天搞| 12—13女人毛片做爰片一| 桃色一区二区三区在线观看| 国产伦精品一区二区三区视频9 | 国产伦在线观看视频一区| 男女视频在线观看网站免费| 1000部很黄的大片| 久久久国产成人精品二区| 九九久久精品国产亚洲av麻豆 | 久久这里只有精品中国| 脱女人内裤的视频| 少妇裸体淫交视频免费看高清| 亚洲av五月六月丁香网| 麻豆国产av国片精品| 99国产极品粉嫩在线观看| 一卡2卡三卡四卡精品乱码亚洲| 婷婷六月久久综合丁香| 成人欧美大片| а√天堂www在线а√下载| 亚洲人成网站在线播放欧美日韩| 无人区码免费观看不卡| 午夜a级毛片| 啦啦啦观看免费观看视频高清| 天天添夜夜摸| 一二三四在线观看免费中文在| 精品一区二区三区视频在线 | 免费观看人在逋| 日本成人三级电影网站| 午夜激情福利司机影院| 国产主播在线观看一区二区| 成人欧美大片| 欧美zozozo另类| 动漫黄色视频在线观看| 国产日本99.免费观看| 草草在线视频免费看| 日日夜夜操网爽| 日本a在线网址| 亚洲电影在线观看av| x7x7x7水蜜桃| 日韩欧美一区二区三区在线观看| 一级毛片高清免费大全| 老鸭窝网址在线观看| 成人av一区二区三区在线看| 欧美成人免费av一区二区三区| 亚洲av免费在线观看| 亚洲成av人片在线播放无| 最近最新免费中文字幕在线| 精品电影一区二区在线| 99热这里只有精品一区 | 伦理电影免费视频| 亚洲av电影在线进入| 一个人看的www免费观看视频| 麻豆成人午夜福利视频| 中文字幕人成人乱码亚洲影| 欧美最黄视频在线播放免费| 岛国在线免费视频观看| 女同久久另类99精品国产91| 亚洲精品久久国产高清桃花| 麻豆成人av在线观看| 夜夜爽天天搞| 亚洲午夜精品一区,二区,三区| 狠狠狠狠99中文字幕| 国产精品99久久久久久久久| 听说在线观看完整版免费高清| 91在线观看av| 久久中文字幕一级| 久久精品亚洲精品国产色婷小说| 免费观看人在逋| 不卡av一区二区三区| 亚洲av熟女| 免费无遮挡裸体视频| 18美女黄网站色大片免费观看| 在线观看免费视频日本深夜| 性色avwww在线观看| 亚洲成人免费电影在线观看| 国产精华一区二区三区| 精品一区二区三区视频在线观看免费| 亚洲成av人片免费观看| 一进一出抽搐gif免费好疼| 国产美女午夜福利| 国产av在哪里看| 中出人妻视频一区二区| 国产精品香港三级国产av潘金莲| 老司机福利观看| 亚洲,欧美精品.| 特级一级黄色大片| 国产欧美日韩一区二区精品| 午夜福利欧美成人| 看片在线看免费视频| 最近最新免费中文字幕在线| 人人妻,人人澡人人爽秒播| 国产av一区在线观看免费| 中文字幕熟女人妻在线| 国产精品亚洲一级av第二区| 午夜激情福利司机影院| 色吧在线观看| 国产午夜精品论理片| 成年人黄色毛片网站| 五月伊人婷婷丁香| 99久久精品一区二区三区| 亚洲五月天丁香| 日本五十路高清| 国产高清三级在线| 最近最新中文字幕大全免费视频| 午夜福利高清视频| 国产av一区在线观看免费| 黄色视频,在线免费观看| 动漫黄色视频在线观看| 成人三级做爰电影| 色综合站精品国产| 国产精品亚洲一级av第二区| 国产精品亚洲av一区麻豆| 一二三四社区在线视频社区8| 激情在线观看视频在线高清| 精品无人区乱码1区二区| 国产欧美日韩一区二区三| 国产午夜精品论理片| 一本一本综合久久| 久久精品夜夜夜夜夜久久蜜豆| 日韩欧美在线乱码| 久久性视频一级片| 午夜福利欧美成人| 日韩欧美在线乱码| 毛片女人毛片| 男女下面进入的视频免费午夜| 色综合婷婷激情| 国产精华一区二区三区| 午夜福利欧美成人| 亚洲av成人一区二区三| 亚洲成av人片免费观看| 国内揄拍国产精品人妻在线| 欧美日本视频| 给我免费播放毛片高清在线观看| 麻豆av在线久日| 亚洲av成人不卡在线观看播放网| 久久久久精品国产欧美久久久| 免费观看人在逋| 国产爱豆传媒在线观看| 午夜久久久久精精品| 亚洲精品中文字幕一二三四区| 久久精品人妻少妇| 男女那种视频在线观看| 亚洲国产精品成人综合色| 国产av在哪里看| 国产探花在线观看一区二区| 午夜福利在线观看吧| 国产精品久久久av美女十八| 真人做人爱边吃奶动态|