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

    喇叭管懸空高度對井筒式泵站運行性能影響研究

    2021-08-25 06:14:42顧雙殷憲柱劉剛許生
    人民長江 2021年6期
    關(guān)鍵詞:懸空漩渦泥沙

    顧雙 殷憲柱 劉剛 許生

    摘要:為了探究喇叭管懸空高度對井筒式泵站運行性能以及泥沙沉積的影響,采用CFX對井筒式泵站內(nèi)流動特性進行了CFD數(shù)值模擬。模擬過程中,將井筒式泵站內(nèi)計算流體視為固液兩相流,設(shè)計了0.3D,0.7D,1.1D和1.5D 4個不同的喇叭管懸空高度方案。計算結(jié)果表明:方案1~4中泵裝置運行最高效率分別為76.4%、79.3%、78.5%和78.8%,方案2中泵裝置的運行效率以及最高高效區(qū)范圍最大;方案1中存在大尺度漩渦,方案3和方案4中存在泥沙沉積現(xiàn)象,因此認為最優(yōu)方案為方案2,該方案的優(yōu)點是沒有泥沙沉積、漩渦尺度小,并且流場結(jié)構(gòu)穩(wěn)定、總水力損失最小。通過分析喇叭管懸空高度對井筒式泵站運行性能和泥沙沉積情況的影響,說明只有合適的喇叭管懸空高度,才能夠保證井筒式泵站的安全穩(wěn)定運行。研究成果可為工程中井筒式泵站設(shè)計提供一定的參考。

    關(guān) 鍵 詞:

    喇叭管懸空高度; 井筒式泵站; CFD; 泥沙淤積; 流場特性; 漩渦尺度

    中圖法分類號: TV732

    文獻標志碼: A

    DOI:10.16232/j.cnki.1001-4179.2021.06.032

    0 引 言

    隨著我國城鎮(zhèn)規(guī)模的不斷擴大,產(chǎn)生了越來越多的城市污水。在城市市政工程中,預(yù)制泵站得到了不斷發(fā)展和應(yīng)用,井筒式泵站則是一種常見的預(yù)制泵站[1]。一體化泵站采用的進水形式是井筒式的進水形式,但這種進水形式泵站中的水流在從進水管進入到井內(nèi)的過程中,往往會比較紊亂,容易形成漩渦[2];城市污水中含有各種泥沙、雜草以及生活垃圾。由于井筒式泵站屬于小型泵站,在實際工程中發(fā)現(xiàn)泵站長時間運行后,在井筒底部會出現(xiàn)泥沙和雜物沉積的問題,既影響泵站的安全運行,又會因泵站內(nèi)沉積生活垃圾產(chǎn)生的腐臭味而影響檢修,因此井筒式泵站的安全運行和自清潔問題一直以來是井筒式泵站面臨的亟需解決的問題。

    查智力等[3-4]針對井筒式泵站內(nèi)水泵數(shù)量、吸水管直徑對泵裝置運行性能的影響進行了CFD數(shù)值模擬計算,得出了井筒式泵站內(nèi)水泵數(shù)量以及吸水管直徑對泵裝置性能有很大影響的結(jié)論。王東進[5]對預(yù)制泵站與傳統(tǒng)混凝土泵站的區(qū)別進行了對比分析和介紹。王卓穎等[6]對預(yù)制泵站進行了優(yōu)化研究。胡凱[7]通過計算得到了預(yù)制泵站優(yōu)化后的有效容積。李文全等[8]對井筒式軸流泵站出水彎管流動特性進行了數(shù)值研究,從而獲得了進水口流量的改變對出水管道水力特性的影響規(guī)律。王默[9]采用CFD數(shù)值模擬方法,研究了不同水泵安裝位置、底部形式和導(dǎo)水錐幾何尺寸對其水力性能的影響。

    目前,針對一體化泵站的研究更多的是集中在結(jié)構(gòu)改進和推廣應(yīng)用方面,而關(guān)于一體化泵站性能優(yōu)化方面的研究卻很少。尤其是應(yīng)當(dāng)重點研究有關(guān)喇叭管的懸空高度對井筒式泵站內(nèi)泥沙沉積和泵裝置運行性能的影響,因為在其他立式泵站中,喇叭管的懸空高度對于水泵的安全運行而言是一個重要的幾何參數(shù),但是關(guān)于井筒式泵站內(nèi)喇叭管懸空高度對于井筒式泵站運行性能及泥沙沉積方面的影響研究目前尚未見到。因此,本文采用CFD數(shù)值模擬的方法來研究喇叭管懸空高度對井筒式泵站性能的影響,研究成果將有助于解決井筒式泵站內(nèi)的自清潔問題以及安裝設(shè)計問題。

    1 計算模型和方法

    1.1 控制方程

    對井筒式泵站進行三維全工況數(shù)值計算。流體力學(xué)的三大基本定律是質(zhì)量守恒定律、能量守恒定律和動量守恒定律[10]。通常情況下,泵裝置內(nèi)的流體可以視其為不可壓縮流體,在整體流動的過程中未考慮熱量傳遞,故可忽略能量守恒定律。流體流動的控制方程組精確地描述了流體的各物理量在空間的分布和隨時間的演化,對于不可壓縮流體,主要是滿足質(zhì)量守恒和動量守恒。

    質(zhì)量守恒方程:

    ρt+(ρui)xi=0(1)

    動量方程:

    (ρui)t+(ρuiuj)xj=-pxi+μ2uixixj(2)

    式中:ρ為流體密度,kg/m3;ui、uj分別為流體在i、j方向上的速度分量,m/s;p為時均壓力,kPa;μ為流體的運動黏性系數(shù),m2/s;xi、xj為空間坐標分量,m。

    1.2 幾何模型的建立

    采用UG 9.0對計算模型進行三維幾何建模,幾何模型部件包括井筒筒體、喇叭管、葉輪、導(dǎo)葉和出水管。基準中心選取葉輪中心為原點,沿進水管進水方向為X方向、與筒體底部平行且垂直進水管方向為Y方向、垂直筒體底部向上為Z方向,如圖1所示。圖1中,葉輪直徑D為120 mm,葉頂間隙為0.1 mm,輪轂直徑為48 mm、葉片數(shù)為4、葉片安放角為0°。喇叭管采用1/4橢圓,長軸和短軸分別為D和0.6D,喇叭管進口直徑為1.6D。井筒筒體直徑為10.0D,出水管直徑為D。為了探究喇叭管懸空高度對井筒式泵站運行性能的影響,分別設(shè)置了4個方案的喇叭管懸空高度,分別為0.3D、0.7D、1.1D和1.5D。

    1.3 網(wǎng)格劃分及無關(guān)性驗證

    采用軟件Workbench中Mesh功能,對井筒泵站筒體、喇叭管和出水管進行幾何模型網(wǎng)格剖分;采用Turbogrid對葉輪和導(dǎo)葉進行網(wǎng)格剖分,如圖2所示。由于黏性的存在,流體邊界層分為黏性底層、過渡層,湍流層,對于高雷諾數(shù)的湍流層y+值取值為30~300[11]。根據(jù)大量類似計算實例,應(yīng)將y+值控制在150以內(nèi),因此,本文計算模型的邊界層網(wǎng)格y+值整體控制在100左右。

    在數(shù)值模擬中,網(wǎng)格數(shù)量和網(wǎng)格質(zhì)量決定了計算結(jié)果精度的可靠性。理論上,隨著網(wǎng)格密度增大,計算結(jié)果的離散誤差減小,計算結(jié)果越準確[12-13]。但是網(wǎng)格密度的增加會加計算工作量,從而造成網(wǎng)格資源的浪費;另外,計算機浮點運算造成的舍入誤差也會隨之增大。為此,對計算模型進行了網(wǎng)格無關(guān)性分析,在此基礎(chǔ)上選取合適的網(wǎng)格數(shù)量,接著采用設(shè)計工況效率來評判網(wǎng)格的無關(guān)性,如圖3所示。當(dāng)網(wǎng)格數(shù)量超過420萬時,設(shè)計工況下的井筒式泵站效率無明顯變化,相對誤差控制在±1%以內(nèi),顯然,此時的網(wǎng)格數(shù)量可滿足計算要求,因此本文總網(wǎng)格數(shù)量采用420萬進行計算。

    1.4 計算參數(shù)及邊界條件設(shè)置

    井筒式泵站主要用于污水工程中,泥沙雜質(zhì)多,本文采用自定義泥沙和水體混合液體模擬工程中的污水,泥沙采用粗砂,密度ρs為1 650 kg/m3,液相為常溫下的清水,密度ρw為997 kg/m3,混合后的密度為

    ρ=1-VsV·ρs+VwV·ρw(3)

    式中:Vs、Vw和V分別為固體體積、液體體積和混合體積,m3。

    邊界條件對流場特性及漩渦流動的數(shù)值計算精度、數(shù)值計算效率和數(shù)值計算穩(wěn)定性都具有明顯的影響[13-14]。選擇適宜的邊界條件處理方法,對精確模擬各類復(fù)雜流態(tài)和漩渦的特性意義重大。

    將進水管進口設(shè)為計算模型液相進口,液相進口采用總壓(Total pressure)進口條件,壓力為1個大氣壓(Total pressure=1atm)。將出水管出口設(shè)置為計算模型出口,出口采用質(zhì)量流量(mass flow rate)出口條件,根據(jù)設(shè)計方案選擇計算流量。在三維建模時,計算模型分為多個部分,不同部件之間需要設(shè)置交界面。除喇叭管出口和葉輪進口、葉輪出口和導(dǎo)葉進口之間的動靜交界面設(shè)置為Stage,將其余靜止部件之間的交界面連接類型設(shè)置為None。計算模型壁面條件采用光滑無滑移壁面條件,包括筒體壁面、喇叭管壁面、葉輪葉片壁面、葉輪室壁面、導(dǎo)葉葉片壁面、導(dǎo)葉室壁面、出水管壁面。葉輪頭導(dǎo)水錐壁面采用旋轉(zhuǎn)壁面(Rotating Wall),壁面轉(zhuǎn)速等于葉輪轉(zhuǎn)速。對計算模型進行初始化設(shè)置,為保證計算流體為固液混合體,將水體初始體積分數(shù)取為0.8,將固體初始體積分數(shù)取為0.2。

    2 計算結(jié)果分析

    2.1 不同方案的能量性能及流場穩(wěn)定性分析

    對不同方案進行全流量工況計算,得到了泵裝置全流量工況性能曲線,如圖4所示。由圖4可知:方案1、方案2、方案3和方案4的泵裝置運行的最高效率點分別為76.4%、79.3%、、78.5%和78.8%,并且方案2的泵裝置運行的高效區(qū)范圍最大。這說明喇叭管的懸空高度為0.7D時,泵裝置的運行效率最高。當(dāng)喇叭管懸空高度為0.3D時,因喇叭管懸空高度過小而使得泵裝置運行效率減小;當(dāng)喇叭管懸空高度為1.1D和1.5D時,喇叭管懸空高度因過大而同樣會使得泵裝置運行效率減小。

    為進一步分析井筒式泵站進口的流動穩(wěn)定性,截取進水管中心高度的X-Y斷面進行流場分析,主要是分析不同方案下X方向的速度Vx的分布及其流線情況,如圖5所示。

    方案1~4中,井筒式泵站喇叭管懸空高度依次增大。從圖5可以看出:當(dāng)懸空高度為0.3D時,水流從進水管進入到井筒泵站,筒體內(nèi)出現(xiàn)了偏流,沿葉輪旋轉(zhuǎn)方向偏向筒體左側(cè),水流高速沖擊到筒體左側(cè)壁,使得筒體右側(cè)流場存在大尺度回流漩渦,從而造成了進水流場結(jié)構(gòu)的不穩(wěn)定,這是由于喇叭管懸空高過低,水泵吸力過大。當(dāng)喇叭管懸空高度為0.7D時,水流從進水管能夠分左右兩股水流平穩(wěn)地進入泵站筒體內(nèi),水流高速進入筒體內(nèi),并且高速水流的范圍幾乎充滿整個斷面,這是由于喇叭管懸空高度較低,水泵喇叭管進口的吸引力非常大,受葉輪旋轉(zhuǎn)作用,筒體左側(cè)的高速流場大于右側(cè)流場使得筒體右側(cè)存在很小的小尺度漩渦。當(dāng)喇叭管懸空高度繼續(xù)增大為1.1D時,由于喇叭管懸空高度增加,水泵喇叭管口吸力減小,從進水管口進入筒體內(nèi)的速度大幅減小,當(dāng)喇叭管懸空高度達到1.5D時,由于喇叭管口與筒體底部之間縱截面過流面積增大,水泵吸力減小使得從進水管口進入筒體內(nèi)的流速進一步減小,在筒體左右兩側(cè)均存在漩渦,左側(cè)漩渦尺度相對于方案3減小,但是右側(cè)漩渦尺度相對于方案3增大非常明顯。說明喇叭管懸空高度過大不利于進水流場的穩(wěn)定分布,易誘發(fā)回流和產(chǎn)生漩渦,喇叭管懸空高度對井筒式泵站進水流場結(jié)構(gòu)的穩(wěn)定性有很大影響。

    2.2 不同方案下的漩渦特性分析

    由于流場分布的不均勻使得在立式泵站喇叭管下方會有附底漩渦存在的可能性,因此,對不同方案的流場采用Q準則三維漩渦進行了可視化處理,從而得到了不同喇叭管懸空高度方案下井筒式泵站內(nèi)可能存在的漩渦,如圖6所示。當(dāng)喇叭管懸空高度為0.3D時,在喇叭管口下方存在著很強的附底漩渦。根據(jù)相關(guān)立式軸流泵喇叭管下方漩渦的有關(guān)文獻可以知道,喇叭管懸空高度越小,漩渦持續(xù)時間越長。當(dāng)喇叭管懸空高度為0.7D時,在喇叭管下方存在尺度很小的附底漩渦,當(dāng)喇叭管懸度空高度增加到1.1D和1.5D時,在喇叭管下方均未發(fā)生附底漩渦。這說明喇叭管的懸空高度增加后,可以有效避免喇叭管下方附底漩渦的發(fā)生。

    通過對井筒式泵站內(nèi)漩渦進行可視化處理,得到了不同喇叭管懸空高度下井筒式泵站喇叭管下方存在漩渦發(fā)生的可能性。為了能夠定量分析喇叭管下方漩渦強度的特性,采用渦量進行評判。渦量是描寫旋渦運動的重要物理量之一[15-16],定義為流體速度矢量的旋度,漩渦通常用渦量來量度其強度和方向,渦量對水頭損失影響較大。對喇叭管進口的流場分布和渦量分布狀況進行了分析,圖7所示為不同方案下的喇叭管進口流線分布和渦量分布的狀況。

    從喇叭管進口的渦量分布和流線分布可以看出:在方案1和方案2中,由于喇叭管懸空高度過低,其周圍水流從喇叭管四周匯聚流入到喇叭管內(nèi),過低的喇叭管懸空高度增加了喇叭管下方流場的混摻,葉輪的旋轉(zhuǎn)作用加強了喇叭管下方流場的剪切運動,誘發(fā)附底漩渦產(chǎn)生,喇叭管下方的紊流強度大。方案3和方案4中的喇叭管進口渦量分布和流線分布在當(dāng)喇叭管懸空高度增大后,喇叭管口與筒體底部之間縱截面面積增大,因此橫向匯集速度減小,周圍水流在喇叭管口匯聚,但是沒有明顯的混摻,水流紊流強度不大。工程中觀察到,喇叭管下方水流紊流強度越大泥沙越容易被攪拌排除,橫向流速越小,水流越穩(wěn)定,泥沙越容易沉積。

    2.3 泥沙沉積含量分析

    對不同方案下的井筒式泵站內(nèi)的泥沙分布情況進行了分析,采用流體體積分數(shù)法顯示泵站內(nèi)的泥沙分布,固體體積分數(shù)的分布代表了井筒式泵站底部的泥沙淤積情況,如圖8所示。在方案1中,井筒底部未見泥沙沉積。在方案2中,在喇叭管附近存在著固體體積分數(shù)為0.7的固液混合體,說明在方案2下尚未形成泥沙沉積的問題。在方案3中,當(dāng)喇叭管懸空高度增加到1.1D時,在井筒底部存在著固體體積分數(shù)為1.0的泥沙分布,說明在井筒底部形成了泥沙沉積;當(dāng)喇叭管懸空高度增加到1.5D時,井筒底部的泥沙沉積明顯增加。方案1~4中的泥沙分布特性說明:喇叭管懸空高度越大,泥沙越容易沉積。

    由圖7和圖8可知:喇叭管懸空高度減小后,行進速度會增大,有利于泥沙的混摻,泥沙被充分攪動,隨主流被水泵抽吸排除;當(dāng)喇叭管懸空高度過大時,喇叭管與筒體底部之間的行進空間過大,使得水流的行進流速過小,泥沙難以被充分攪動與主流混摻,流速過小時,加上水泵的抽吸作用減小使得泥沙產(chǎn)生沉積;同時水中的雜草進入泥沙中,導(dǎo)致淤積情況更加嚴重。

    對喇叭管下方附底漩渦發(fā)生的情況進行分析可以發(fā)現(xiàn):雖然減小喇叭管懸空高度可以減小泥沙沉積,但是喇叭管的懸空高度過低會導(dǎo)致喇叭管下方產(chǎn)生附底漩渦。通過眾多工程案例和相關(guān)資料可知:附底漩渦對機組的安全運行會有很大的危害,將導(dǎo)致機組產(chǎn)生嚴重的振動和葉輪內(nèi)汽蝕,所以在井筒式泵站內(nèi)要做到既能夠有效減少泥沙沉積,又能夠保證泵裝置喇叭管下方不會產(chǎn)生附底漩渦。因此,需要選擇合適的喇叭管懸空高。

    2.4 流動內(nèi)特性分析

    通過對井筒式泵站內(nèi)流動外特性分析可知:井筒式泵站內(nèi)存在漩渦和泥沙沉積的現(xiàn)象,而不同喇叭管懸空高度影響井筒式泵站內(nèi)漩渦和泥沙情況。這些漩渦和泥沙必然對泵裝置的性能產(chǎn)生影響。因此,有必要對不同方案下的井筒泵站內(nèi)的內(nèi)流場特性進行分析。

    為了進一步分析喇叭管懸空高度對井筒式泵站葉輪進口速度分布均勻性的影響,對葉輪進口斷面速度場進行了均勻性分析。分析過程中,流速均勻度采用單元面積為權(quán)重進行加權(quán)計算,以斷面流量為基礎(chǔ)的面積加權(quán)流速均勻度更合理[17],所以本文采用面積加權(quán)流速分布均勻度來評判葉輪進口流速分布情況。面積加權(quán)流速均勻度Vuna計算公式如下:

    Vuna=1-1vani=1νai-va2ΔAiA(4)

    式中:ΔAi為葉輪進口斷面第i單元的面積,m2;

    va為葉輪進口斷面的平均軸向速度,m/s;

    vai為葉輪進口斷面第i單元的軸向速度m/s;

    Q為進口流量,m3/s;

    A為葉輪進口面積,m2;

    va=Q/A。

    圖9為采用公式(4)獲得的設(shè)計工況下不同方案的葉輪進口斷面流速均勻度分布狀況。由圖9可以看出:設(shè)計工況下,方案1~4的葉輪進口流速均勻度分別為88.2%、92.3%、90.50%和89.60%,可以看到在不同喇叭管懸空高度下,葉輪進口流速的均勻度分布差異很大。

    (1) 方案1中,喇叭管懸空高度最低,葉輪進口前水流混摻強烈,紊流強度大和漩渦的存在使得葉輪進口流場分布差,因此,葉輪進口的流速均勻度最差。

    (2) 方案2中,葉輪進口流速均勻度分布最優(yōu),方案中喇叭管懸空高度合適,水流能夠平順地從喇叭管四周匯入喇叭管內(nèi),因此,葉輪進口處流速均勻度高。

    (3) 方案3中,流速均勻度介于方案1和方案2之間,這是由于方案3中的井筒底部存在著泥沙沉積的問題,影響到了進入喇叭管內(nèi)的流態(tài)。

    (4) 在方案4中,葉輪進口流速均勻度低于方案3,這是由于方案四中的喇叭管懸空高度增加了,泥沙沉積更加嚴重,對葉輪進口的流態(tài)影響更強。

    通過對不同方案下的葉輪進口流速的均勻度進行分析,可以得到喇叭管懸空高度對葉輪進口流速均勻度影響很大。

    根據(jù)對不同方案下泵裝置外特性的分析,可以得到井筒式泵站內(nèi)喇叭管懸空高度對機組性能的影響。軸流泵裝置揚程低的特點導(dǎo)致井筒內(nèi)的水頭損失和出水管的損失占裝置揚程很大的比重,對水泵裝置運行效率產(chǎn)生明顯影響。以泵裝置水力損失作為評判指標,更進一步地探究不同喇叭管懸空高對泵裝置的影響。水力損失計算公式如下[18]:

    Δh=Pin-Poutρg(5)

    式中:Δh為水力損失,m;Pin為進口處總壓,kPa;Pout為出口總壓,kPa;ρ為水的密度,kg/m3;g為重力加速度,m/s2。

    根據(jù)公式(5)獲得了不同方案設(shè)計工況下的總水力損失,如圖10所示。由圖10可以看出:方案1~4的總水力損失分別為0.31,0.16,0.23 m和0.26 m;方案1中泵裝置的總水力損失最大,方案2的總水力損失最小。

    通過對不同方案的能量性能、流場特性、漩渦特性以及泥沙含量進行分析,最終確定方案2為最優(yōu)方案;漩渦和泥沙淤積都會對泵裝置的安全運行產(chǎn)生危害,由于泥沙淤積是持續(xù)存在,并且泥沙淤積后喇叭管口與筒體底部之間的距離減小,過流面積減小同樣造成漩渦的產(chǎn)生;而在方案2中泥沙難以淤積,同時漩渦具有不穩(wěn)定性,不是持續(xù)發(fā)生,并且泥沙沉積的危害大于漩渦的危害,因此,最優(yōu)方案為方案2。分析結(jié)果表明,只有合適的喇叭管懸空高度才能夠保證井筒式泵站的安全穩(wěn)定運行。

    3 結(jié) 論

    針對井筒式泵站內(nèi)泥沙沉積導(dǎo)致的自清潔問題,采用軟件CFX對井筒式泵站內(nèi)的流動特性進行了數(shù)值模擬;通過自定義流體介質(zhì)為固液兩相流,分析了井筒式泵站內(nèi)的流量工況、喇叭管懸空高度對泥沙沉積的影響及其對泵裝置性能的影響。得出如下結(jié)論。

    (1) 方案1~4的泵裝置運行的最高效率分別為76.4%、79.3%、、78.5%和78.8%,并且方案2的泵裝置運行的高效區(qū)范圍最大。方案1~4的總水力損失分別為0.31,0.16,0.23 m和0.26 m。葉輪進口流速均勻度分別為88.2%、92.3%、84.6%和88.6%。

    (2) 在方案1和方案2中,由于喇叭管懸空高度過低,周圍水流從喇叭管四周匯聚進入喇叭管內(nèi),過低的喇叭管懸空高度增加了喇叭管下方流場的混摻,葉輪的旋轉(zhuǎn)作用加強了喇叭管下方流場的剪切運動,從而誘發(fā)了附底漩渦產(chǎn)生并使得喇叭管下方的紊流強度大。在方案3和方案4中未見漩渦現(xiàn)象的發(fā)生,表明喇叭管懸空高度增大后可以避免漩渦的發(fā)生。

    (3) 在方案1和方案2中,井筒底部未見明顯泥沙沉積;在方案3和方案4中井筒底部存在泥沙沉積,而且喇叭管懸空高度越大,泥沙越容易沉積。喇叭管越小,喇叭管口下方流場的混摻強度越大,有利于泥沙被及時排除。

    (4) 漩渦和泥沙淤積都會對泵裝置的安全運行產(chǎn)生危害。由于泥沙淤積是持續(xù)存在,漩渦具有不穩(wěn)定性不是持續(xù)發(fā)生,并且泥沙沉積的危害大于漩渦的危害,因此最優(yōu)方案為方案2,這也表明,只有合適的喇叭管懸空高度,才能夠保證井筒式泵站的安全穩(wěn)定運行。

    參考文獻:

    [1] 周佳.一體化污水提升泵站的應(yīng)用探討[J].山西建筑,2014,41(10):130-131.

    [2] LEVEY TOM,JEFFUS,STEPHEN.Prefabricated pumping systems:Sizing up the possibilities[J].Water Engineering and Management,2003,150(2):23-25.

    [3] 查智力,劉超,嚴天序,等.井筒式泵裝置雙泵模型進水特性[J].流體機械,2019,47(5):44-49.

    [4] 查智力,劉超,嚴天序,等.井筒式泵裝置水力特性數(shù)值模擬[J].南水北調(diào)與水利科技,2018,16(2):189-195.

    [5] 王東進.一體化預(yù)制泵站設(shè)計開發(fā)研究及應(yīng)用[J].通用機械,2014(7):77-79

    [6] 王卓穎,吳恩賜,郭德若,等.預(yù)制泵站在中小型泵站應(yīng)用中的優(yōu)化研究[J].中國給水排水,2010,26(22):127.

    [7] 胡凱.一體化預(yù)制泵站的自清潔[J].通用機械,2019(11):15-17.

    [8] 李文全,王凱,陳昆,等.不同運行方式對一體化預(yù)制泵站流動特性和筒體強度的影響[J].中國農(nóng)村水利水電,2019(4):162-167,171.

    [9] 王默.一體化泵站水力性能CFD研究[D].揚州:揚州大學(xué),2016.

    [10] 張自超,王福軍,陳鑫,等.基于改進歐拉算法的雙吸離心泵泥沙磨損特性研究[J].農(nóng)業(yè)機械學(xué)報,2017,48(3):124-133,147.

    [11] 張自超,陳鑫,王福軍,等.考慮顆粒動態(tài)尺度影響的泥沙擴散系數(shù)模型建立及應(yīng)用[J].農(nóng)業(yè)工程學(xué)報,2016,32(22):129-137.

    [12] 許巧哲,崔碩垚.山區(qū)小城鎮(zhèn)污水集中處理設(shè)計[J].給水排水,2016,52(增1):160-162.

    [13] SONG X J,LIU C.Experimental investigation of floor-attached vortex effects on the pressure pulsation at the bottom of the axial flow pump sump[J].Renewable Energy,2019,145(2020):2327-2336.

    [14] 李嘉興,王福軍,陳鑫.基于雙流體模型的尊村泵站雙吸離心泵泥沙磨損[J].科學(xué)通報,2019,64(36):3856-3866.

    [15] 宋希杰,劉超,羅燦,等.軸流泵裝置進水漩渦對壓力脈動的影響[J].農(nóng)業(yè)機械學(xué)報,2018,49(2):113-119.

    [16] 宋希杰,劉超,楊帆,等.進水漩渦誘發(fā)軸流泵壓力脈動的試驗研究[J].振動與沖擊,2018,37(17):25-31.

    [17] 宋希杰,劉超,楊帆,等.水泵進水池底部壓力脈動特性試驗[J].農(nóng)業(yè)機械學(xué)報,2017,48(11):196-203.

    [18] 程煬,孫眾,張慧,等.葉片角度對采用肘形流道的軸流泵裝置特性的影響[J].流體機械,2019,47(6):29-33.

    (編輯:趙秋云)

    Study on influence of suspended height of trumpet tube on operation

    performance of shaft pumping station

    GU Shuang1,YIN Xianzhu2,LIU Gang1,XU Sheng2

    (1.Jiangsu Luoyun Water Conservancy Project Management Office,Suqian 223800,China; 2.Suqian Water Survey,Design and Research Institute Co.,Ltd,Suqian 223800,China)

    Abstract:

    In order to explore the influence of the suspended height of the trumpet tube on the operation performance and sediment deposition of the shaft pumping station,this paper uses CFX to simulate flow characteristics in the shaft pumping station.The solid-liquid two-phase flow is adopted for the calculation fluid in the shaft pumping station.Four schemes of 0.3D,0.7D,1.1D and 1.5D of different suspended height of the trumpet tube are designed.The calculation results show that the maximum efficiency of scheme 1 to scheme 4 is 76.4%,79.3%,78.5% and 78.8% respectively,and scheme 2 has the highest operation efficiency of and the widest efficiency range.There are large-scale vortices in scheme 1,and sediment deposition in schemes 3 and 4.The optimal scheme is scheme 2,and it is characterized as no sediment deposition,small vortex scale,stable flow field structure and minimum total hydraulic loss.By analyzing influence of the suspended height of trumpet tube on the operation performance and sediment deposition of the shaft pumping station,it is found that only the appropriate suspended height of the trumpet tube can ensure the safe and stable operation of the shaft pumping station.The research results can provide some reference for the design of the shaft pumping station in the project.

    Key words:

    suspended height of the trumpet tube;shaft pumping station;CFD;sediment deposition;flow field characteristics;vortex scale

    猜你喜歡
    懸空漩渦泥沙
    泥沙做的父親
    新疆多泥沙河流水庫泥沙處理措施
    土壤團聚體對泥沙沉降速度的影響
    FF陷控制權(quán)爭奪漩渦
    汽車觀察(2018年10期)2018-11-06 07:05:06
    大寶小神探·山澗懸空城
    琴童(2018年2期)2018-03-06 16:36:52
    魚群漩渦
    中外文摘(2017年19期)2017-10-10 08:28:41
    懸空的水
    基于MPU6050六軸傳感器的懸空鼠標設(shè)計與實現(xiàn)
    電子制作(2016年15期)2017-01-15 13:39:11
    中醫(yī)教育陷“量升質(zhì)降”漩渦
    泥沙滅火
    兒童繪本(2015年2期)2015-05-25 18:10:15
    涩涩av久久男人的天堂| 黄色日韩在线| 人妻少妇偷人精品九色| 伊人久久精品亚洲午夜| 超碰97精品在线观看| 少妇的逼好多水| 一级毛片aaaaaa免费看小| 久久精品久久久久久噜噜老黄| 亚洲精华国产精华液的使用体验| 熟女人妻精品中文字幕| 成人漫画全彩无遮挡| 少妇人妻 视频| 黄色日韩在线| 欧美97在线视频| 亚洲av综合色区一区| 一级毛片电影观看| 2022亚洲国产成人精品| 国产精品99久久久久久久久| 如何舔出高潮| 成人漫画全彩无遮挡| 蜜臀久久99精品久久宅男| 欧美老熟妇乱子伦牲交| 久久精品夜色国产| 精品久久久久久久末码| 91在线精品国自产拍蜜月| 2018国产大陆天天弄谢| 三级国产精品片| 天天躁日日操中文字幕| 91精品国产国语对白视频| 国产高清有码在线观看视频| 久久国内精品自在自线图片| 在现免费观看毛片| 最近手机中文字幕大全| 国产精品久久久久久av不卡| 视频中文字幕在线观看| 99久久人妻综合| 精品久久久精品久久久| 亚洲自偷自拍三级| 纵有疾风起免费观看全集完整版| h视频一区二区三区| 51国产日韩欧美| 国产精品一区二区在线不卡| 国产真实伦视频高清在线观看| 亚洲欧美日韩东京热| 国产91av在线免费观看| 欧美成人一区二区免费高清观看| 日本vs欧美在线观看视频 | av不卡在线播放| 久久久久久久精品精品| 国产精品欧美亚洲77777| freevideosex欧美| 精品国产露脸久久av麻豆| 国产在线免费精品| 一个人看视频在线观看www免费| 国产av一区二区精品久久 | 亚洲欧美成人精品一区二区| 欧美变态另类bdsm刘玥| 蜜桃久久精品国产亚洲av| 亚洲人与动物交配视频| 欧美一区二区亚洲| 久久人人爽人人爽人人片va| 亚洲美女视频黄频| 高清午夜精品一区二区三区| 直男gayav资源| 亚洲国产精品成人久久小说| 欧美zozozo另类| 欧美少妇被猛烈插入视频| 久久久久久九九精品二区国产| 波野结衣二区三区在线| 国产精品.久久久| 欧美97在线视频| 久久人人爽人人片av| 久久6这里有精品| 99热国产这里只有精品6| 亚洲国产高清在线一区二区三| 久久久久国产网址| 亚洲国产毛片av蜜桃av| 最后的刺客免费高清国语| 大码成人一级视频| 亚洲美女视频黄频| 亚洲中文av在线| 精品久久久噜噜| 日本黄色片子视频| 国模一区二区三区四区视频| 精品久久国产蜜桃| 少妇被粗大猛烈的视频| 永久免费av网站大全| 久久精品国产亚洲av涩爱| 乱码一卡2卡4卡精品| 99热全是精品| 成人无遮挡网站| 男男h啪啪无遮挡| 99久国产av精品国产电影| 女人十人毛片免费观看3o分钟| 高清av免费在线| 精品久久久久久久久亚洲| 日韩av在线免费看完整版不卡| 欧美三级亚洲精品| 亚洲国产高清在线一区二区三| 欧美三级亚洲精品| av.在线天堂| h日本视频在线播放| 新久久久久国产一级毛片| 亚洲成人一二三区av| 最近2019中文字幕mv第一页| 久久精品夜色国产| a级毛色黄片| 人人妻人人添人人爽欧美一区卜 | 亚洲av在线观看美女高潮| 成人无遮挡网站| 亚洲va在线va天堂va国产| 国产高清国产精品国产三级 | 色5月婷婷丁香| 亚洲av中文字字幕乱码综合| 精品久久久噜噜| 青春草国产在线视频| 插阴视频在线观看视频| 美女xxoo啪啪120秒动态图| 亚洲国产精品专区欧美| 男女啪啪激烈高潮av片| h日本视频在线播放| 久久亚洲国产成人精品v| 亚洲国产毛片av蜜桃av| 国产精品久久久久久久久免| 下体分泌物呈黄色| 亚洲婷婷狠狠爱综合网| 搡老乐熟女国产| 国产成人免费观看mmmm| 亚洲熟女精品中文字幕| 午夜福利网站1000一区二区三区| 亚洲中文av在线| 日日摸夜夜添夜夜添av毛片| av在线老鸭窝| 国产欧美日韩一区二区三区在线 | 亚洲经典国产精华液单| 色5月婷婷丁香| 免费少妇av软件| 精品一区二区三卡| 国产视频内射| 春色校园在线视频观看| 欧美日韩国产mv在线观看视频 | 国产精品99久久久久久久久| 网址你懂的国产日韩在线| .国产精品久久| 亚洲精品自拍成人| 纵有疾风起免费观看全集完整版| 国产男女超爽视频在线观看| 黄色一级大片看看| 国产av精品麻豆| 身体一侧抽搐| 永久免费av网站大全| 亚洲无线观看免费| 欧美日本视频| 久久毛片免费看一区二区三区| 久久精品国产亚洲av涩爱| 一个人看的www免费观看视频| 国产精品.久久久| 国产高清不卡午夜福利| 大话2 男鬼变身卡| 人妻制服诱惑在线中文字幕| 免费大片黄手机在线观看| 熟女电影av网| 色网站视频免费| 内射极品少妇av片p| 蜜桃在线观看..| 国产免费又黄又爽又色| 国产av一区二区精品久久 | 看免费成人av毛片| 美女中出高潮动态图| 99热这里只有是精品50| 精品人妻一区二区三区麻豆| 亚洲美女搞黄在线观看| 国产精品久久久久久久久免| 久久影院123| 啦啦啦中文免费视频观看日本| 高清黄色对白视频在线免费看 | 小蜜桃在线观看免费完整版高清| 亚洲国产毛片av蜜桃av| 97在线人人人人妻| 免费黄网站久久成人精品| 国产精品蜜桃在线观看| 韩国高清视频一区二区三区| 韩国av在线不卡| 欧美精品亚洲一区二区| 男人爽女人下面视频在线观看| 一本久久精品| 国模一区二区三区四区视频| 99久久综合免费| 亚洲va在线va天堂va国产| 丰满少妇做爰视频| 最近中文字幕高清免费大全6| videossex国产| 秋霞在线观看毛片| 精品亚洲成国产av| 日日摸夜夜添夜夜添av毛片| 国产av精品麻豆| 久久99精品国语久久久| 日韩不卡一区二区三区视频在线| 国内揄拍国产精品人妻在线| 在线观看一区二区三区| 成人国产麻豆网| 色哟哟·www| 91久久精品国产一区二区成人| 午夜老司机福利剧场| 久久6这里有精品| 国产成人免费无遮挡视频| 卡戴珊不雅视频在线播放| 国产成人aa在线观看| 国产v大片淫在线免费观看| 丝瓜视频免费看黄片| a级一级毛片免费在线观看| h日本视频在线播放| 男女啪啪激烈高潮av片| 日韩av在线免费看完整版不卡| 国产黄频视频在线观看| 成人国产麻豆网| 久久久久久伊人网av| 国国产精品蜜臀av免费| 噜噜噜噜噜久久久久久91| 亚洲高清免费不卡视频| 99热全是精品| 欧美丝袜亚洲另类| 国产 精品1| 亚洲欧美日韩另类电影网站 | 18禁动态无遮挡网站| 亚洲精品456在线播放app| 中文字幕亚洲精品专区| 成年av动漫网址| 中文在线观看免费www的网站| 国产一区有黄有色的免费视频| 在线观看一区二区三区激情| 亚洲丝袜综合中文字幕| 在线播放无遮挡| 久久久久久久久久成人| 国产精品.久久久| 你懂的网址亚洲精品在线观看| 成人毛片60女人毛片免费| 亚洲av国产av综合av卡| 久久国产乱子免费精品| 在线观看一区二区三区| 亚洲欧美成人综合另类久久久| av又黄又爽大尺度在线免费看| 丰满人妻一区二区三区视频av| 久久综合国产亚洲精品| 一本一本综合久久| 欧美极品一区二区三区四区| 国产av一区二区精品久久 | 尤物成人国产欧美一区二区三区| 欧美高清成人免费视频www| 欧美激情国产日韩精品一区| 午夜免费男女啪啪视频观看| 亚洲国产日韩一区二区| 欧美区成人在线视频| 欧美精品一区二区免费开放| 亚洲一级一片aⅴ在线观看| 亚洲精品乱久久久久久| 亚洲国产欧美人成| 亚洲怡红院男人天堂| 黑人高潮一二区| 啦啦啦视频在线资源免费观看| 国产乱来视频区| 国产精品麻豆人妻色哟哟久久| av在线老鸭窝| 久久久欧美国产精品| 国产免费一级a男人的天堂| 亚洲精品乱码久久久久久按摩| 亚洲怡红院男人天堂| av一本久久久久| 美女cb高潮喷水在线观看| 久久久久久人妻| 国产精品国产av在线观看| 日韩国内少妇激情av| 伊人久久国产一区二区| 国产精品麻豆人妻色哟哟久久| 熟女av电影| 欧美最新免费一区二区三区| 亚洲av欧美aⅴ国产| 日韩中字成人| 国产男人的电影天堂91| 综合色丁香网| 毛片女人毛片| 亚洲成人一二三区av| 国产色爽女视频免费观看| 尾随美女入室| 国产高清不卡午夜福利| 久久久午夜欧美精品| 美女国产视频在线观看| 国产视频首页在线观看| 好男人视频免费观看在线| 97精品久久久久久久久久精品| 免费黄网站久久成人精品| 91精品一卡2卡3卡4卡| 精品国产一区二区三区久久久樱花 | 91精品伊人久久大香线蕉| 乱系列少妇在线播放| 国产高清三级在线| 欧美最新免费一区二区三区| 日韩精品有码人妻一区| 亚洲色图综合在线观看| 亚洲精品456在线播放app| 国产日韩欧美在线精品| 国产成人freesex在线| 少妇被粗大猛烈的视频| 女人十人毛片免费观看3o分钟| 国产乱来视频区| 亚洲精华国产精华液的使用体验| 91午夜精品亚洲一区二区三区| 高清av免费在线| 国产亚洲精品久久久com| 美女高潮的动态| 丝袜喷水一区| 一个人看视频在线观看www免费| 22中文网久久字幕| 高清av免费在线| 男人和女人高潮做爰伦理| 久热久热在线精品观看| 国产毛片在线视频| 只有这里有精品99| 国产精品一区二区在线不卡| 2021少妇久久久久久久久久久| 麻豆乱淫一区二区| 成人午夜精彩视频在线观看| 日韩一本色道免费dvd| 日韩大片免费观看网站| 天天躁日日操中文字幕| 精品99又大又爽又粗少妇毛片| 日本黄色日本黄色录像| 观看av在线不卡| 国精品久久久久久国模美| 男人狂女人下面高潮的视频| 国产成人a区在线观看| 九草在线视频观看| 2022亚洲国产成人精品| 又粗又硬又长又爽又黄的视频| 国产成人a∨麻豆精品| 成年免费大片在线观看| 久久精品熟女亚洲av麻豆精品| 久久久久久久久久久免费av| 国产精品偷伦视频观看了| 国产成人精品婷婷| 丰满少妇做爰视频| 久久久久久久大尺度免费视频| 久久久精品免费免费高清| 在线天堂最新版资源| 国产伦精品一区二区三区视频9| 久久久精品94久久精品| 国产精品熟女久久久久浪| 一级毛片aaaaaa免费看小| 久久热精品热| 舔av片在线| 国产v大片淫在线免费观看| 色综合色国产| 日韩精品有码人妻一区| 久久国产亚洲av麻豆专区| 少妇熟女欧美另类| 久久精品国产a三级三级三级| 伦精品一区二区三区| 日本av免费视频播放| 日本爱情动作片www.在线观看| 中文字幕免费在线视频6| 九草在线视频观看| 欧美xxⅹ黑人| 亚洲国产精品专区欧美| 美女内射精品一级片tv| 内地一区二区视频在线| 18+在线观看网站| 亚洲成人中文字幕在线播放| 91在线精品国自产拍蜜月| 18禁动态无遮挡网站| 黄色日韩在线| 国产精品人妻久久久影院| 精品一区二区三区视频在线| 久久这里有精品视频免费| 成人毛片60女人毛片免费| 免费不卡的大黄色大毛片视频在线观看| 国产精品熟女久久久久浪| tube8黄色片| 天堂中文最新版在线下载| av在线蜜桃| 国产永久视频网站| 搡女人真爽免费视频火全软件| 国产成人91sexporn| 日本爱情动作片www.在线观看| 最近的中文字幕免费完整| 亚洲国产精品国产精品| 大陆偷拍与自拍| 精品一区二区免费观看| 日日摸夜夜添夜夜爱| 欧美日韩亚洲高清精品| 日韩,欧美,国产一区二区三区| 亚洲欧美日韩无卡精品| 偷拍熟女少妇极品色| 少妇熟女欧美另类| 校园人妻丝袜中文字幕| 国产视频内射| 国产av一区二区精品久久 | 男女啪啪激烈高潮av片| 九色成人免费人妻av| a 毛片基地| 另类亚洲欧美激情| 免费人妻精品一区二区三区视频| 日韩中文字幕视频在线看片 | 久久女婷五月综合色啪小说| 亚洲精品乱久久久久久| 最近2019中文字幕mv第一页| 内射极品少妇av片p| 免费不卡的大黄色大毛片视频在线观看| 国产成人aa在线观看| 午夜免费观看性视频| 国产一区二区三区av在线| 日本爱情动作片www.在线观看| av线在线观看网站| 久久人妻熟女aⅴ| 国产亚洲最大av| 国产91av在线免费观看| 波野结衣二区三区在线| 色吧在线观看| 成年人午夜在线观看视频| 亚洲人成网站在线观看播放| 一区二区三区精品91| 啦啦啦中文免费视频观看日本| 国产亚洲av片在线观看秒播厂| 中文天堂在线官网| 激情 狠狠 欧美| 欧美成人午夜免费资源| 亚洲欧美精品自产自拍| 青春草亚洲视频在线观看| 亚洲av不卡在线观看| 麻豆成人av视频| 亚洲怡红院男人天堂| 毛片女人毛片| 在线播放无遮挡| 最黄视频免费看| 欧美成人一区二区免费高清观看| 亚洲图色成人| 两个人的视频大全免费| 黄色怎么调成土黄色| 伊人久久国产一区二区| 国产乱人视频| 久久青草综合色| 久久精品国产a三级三级三级| 日本免费在线观看一区| 亚洲内射少妇av| 三级国产精品片| 久久久久国产网址| 美女国产视频在线观看| 国产精品嫩草影院av在线观看| 夜夜爽夜夜爽视频| 在线观看av片永久免费下载| 成人漫画全彩无遮挡| 男人狂女人下面高潮的视频| 色综合色国产| 一本—道久久a久久精品蜜桃钙片| 亚洲四区av| 国产精品人妻久久久久久| 日本av免费视频播放| 嫩草影院入口| 国产色婷婷99| 在线观看免费日韩欧美大片 | 乱码一卡2卡4卡精品| 亚洲精品aⅴ在线观看| 成人漫画全彩无遮挡| 99热这里只有精品一区| 女性生殖器流出的白浆| 亚州av有码| 麻豆乱淫一区二区| 女性生殖器流出的白浆| 日韩中文字幕视频在线看片 | 全区人妻精品视频| 一级毛片电影观看| 亚洲国产成人一精品久久久| av.在线天堂| 舔av片在线| 亚洲精品一区蜜桃| 精品人妻偷拍中文字幕| 国产精品一区二区在线不卡| 美女cb高潮喷水在线观看| 欧美高清成人免费视频www| 精品久久久久久电影网| 免费在线观看成人毛片| 精品国产露脸久久av麻豆| 久久综合国产亚洲精品| 最近最新中文字幕大全电影3| 国产精品国产三级专区第一集| 99久久综合免费| 午夜视频国产福利| 九九久久精品国产亚洲av麻豆| 在线免费观看不下载黄p国产| 久久久久国产网址| 麻豆国产97在线/欧美| 九九爱精品视频在线观看| 成人毛片60女人毛片免费| 国产av国产精品国产| 小蜜桃在线观看免费完整版高清| 99久久中文字幕三级久久日本| 欧美 日韩 精品 国产| 纵有疾风起免费观看全集完整版| 日韩欧美精品免费久久| 尾随美女入室| 欧美精品国产亚洲| 久久久久久久久久久丰满| 久久久久久九九精品二区国产| 精品久久久久久久末码| 秋霞伦理黄片| 亚洲欧美成人综合另类久久久| 麻豆精品久久久久久蜜桃| 高清午夜精品一区二区三区| 偷拍熟女少妇极品色| 亚洲三级黄色毛片| 菩萨蛮人人尽说江南好唐韦庄| 国产精品国产三级国产专区5o| 黄色日韩在线| 一级毛片久久久久久久久女| 欧美97在线视频| 久久精品久久久久久噜噜老黄| 亚洲四区av| 岛国毛片在线播放| 有码 亚洲区| 高清在线视频一区二区三区| 日日啪夜夜撸| 18禁裸乳无遮挡免费网站照片| 亚洲无线观看免费| 五月玫瑰六月丁香| 亚洲国产成人一精品久久久| 欧美精品人与动牲交sv欧美| 亚洲aⅴ乱码一区二区在线播放| 99久久精品一区二区三区| 99久久精品热视频| 国产亚洲一区二区精品| 制服丝袜香蕉在线| 亚洲欧美日韩卡通动漫| 亚洲av中文字字幕乱码综合| 少妇人妻久久综合中文| 啦啦啦啦在线视频资源| 人妻夜夜爽99麻豆av| 最近手机中文字幕大全| 黄色日韩在线| 午夜免费鲁丝| 日韩电影二区| 精品午夜福利在线看| 亚洲精品中文字幕在线视频 | av免费观看日本| 国产视频内射| 久久精品人妻少妇| 在现免费观看毛片| 中文资源天堂在线| 欧美精品国产亚洲| 高清在线视频一区二区三区| 日本av免费视频播放| 麻豆精品久久久久久蜜桃| 两个人的视频大全免费| 亚洲婷婷狠狠爱综合网| 亚洲中文av在线| 久久久久久久精品精品| 熟女av电影| 精品一区在线观看国产| 国产色婷婷99| 18禁裸乳无遮挡动漫免费视频| 蜜桃久久精品国产亚洲av| 国产午夜精品久久久久久一区二区三区| 麻豆成人午夜福利视频| 亚洲成色77777| 亚洲欧美日韩卡通动漫| 欧美老熟妇乱子伦牲交| 精品久久久久久久久av| 老师上课跳d突然被开到最大视频| 黄片无遮挡物在线观看| 麻豆精品久久久久久蜜桃| 亚洲精品视频女| 久久精品夜色国产| 在线观看av片永久免费下载| 在线 av 中文字幕| 寂寞人妻少妇视频99o| 欧美日韩视频高清一区二区三区二| 免费看不卡的av| 国产精品三级大全| 男女无遮挡免费网站观看| 国产一区亚洲一区在线观看| 乱码一卡2卡4卡精品| 99热国产这里只有精品6| 天堂俺去俺来也www色官网| 久久久亚洲精品成人影院| 简卡轻食公司| 久久久久久九九精品二区国产| 在线 av 中文字幕| 成人毛片60女人毛片免费| 国产片特级美女逼逼视频| 国产免费一级a男人的天堂| 99久久精品国产国产毛片| 精品国产一区二区三区久久久樱花 | 亚洲精品乱码久久久久久按摩| 久久久色成人| 欧美三级亚洲精品| 国语对白做爰xxxⅹ性视频网站| 久久久久久久久久久免费av| 久久6这里有精品| 久久久久久久久久成人| 一级二级三级毛片免费看| 蜜桃在线观看..| 国产亚洲最大av| 午夜免费观看性视频| 在线观看三级黄色| 亚洲精品456在线播放app| 我的老师免费观看完整版| 五月天丁香电影| 永久免费av网站大全| 亚洲欧美一区二区三区国产| 亚洲av国产av综合av卡| 久久人人爽人人片av| 汤姆久久久久久久影院中文字幕| 国产精品一二三区在线看| 中文欧美无线码| 妹子高潮喷水视频| 少妇人妻一区二区三区视频| 久久久久久久大尺度免费视频|