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

    泵葉輪出口寬度對蝸殼內(nèi)壓力脈動強(qiáng)度的影響

    2017-08-31 01:50:49顧延?xùn)|袁壽其裴吉張金鳳黃茜王文杰
    關(guān)鍵詞:蝸殼離心泵脈動

    顧延?xùn)|,袁壽其,裴吉,張金鳳,黃茜,王文杰

    (江蘇大學(xué) 國家水泵及系統(tǒng)工程技術(shù)研究中心,江蘇 鎮(zhèn)江 212013)

    泵葉輪出口寬度對蝸殼內(nèi)壓力脈動強(qiáng)度的影響

    顧延?xùn)|,袁壽其,裴吉,張金鳳,黃茜,王文杰

    (江蘇大學(xué) 國家水泵及系統(tǒng)工程技術(shù)研究中心,江蘇 鎮(zhèn)江 212013)

    為了研究低比轉(zhuǎn)數(shù)離心泵葉輪出口寬度對蝸殼內(nèi)壓力脈動強(qiáng)度的影響,將葉輪出口寬度分別設(shè)計(jì)為5、6和7 mm,并保持蝸殼和葉輪其他幾何參數(shù)不變。應(yīng)用ANSYS CFX 14.5進(jìn)行數(shù)值計(jì)算,對比了定常數(shù)值計(jì)算與外特性、粒子圖像測速(particle image velocimetry, PIV)試驗(yàn)結(jié)果,驗(yàn)證了數(shù)值計(jì)算的可信度,3個(gè)方案的揚(yáng)程、效率均滿足設(shè)計(jì)要求。在非定常數(shù)值計(jì)算中,引入壓力脈動強(qiáng)度表征3個(gè)旋轉(zhuǎn)周期內(nèi)壓力脈動的累積效應(yīng)。結(jié)果表明:對于同一葉輪出口寬度,小流量下蝸殼內(nèi)壓力脈動強(qiáng)度最大,大流量次之,而在設(shè)計(jì)工況下,壓力脈動強(qiáng)度及其梯度均是最小的;由于葉輪和蝸殼的時(shí)序干涉作用,隔舌處出現(xiàn)壓力脈動強(qiáng)度的大梯度變化;蝸殼進(jìn)口寬度與葉輪出口寬度之比為3時(shí),該低比轉(zhuǎn)數(shù)離心泵的性能滿足設(shè)計(jì)要求,且壓力脈動強(qiáng)度較低,綜合性能較優(yōu)。

    低比轉(zhuǎn)數(shù);離心泵;數(shù)值模擬;壓力脈動;多工況;葉輪;揚(yáng)程;蝸殼

    低比轉(zhuǎn)數(shù)離心泵具有流量小、揚(yáng)程高等特點(diǎn),廣泛應(yīng)用于農(nóng)業(yè)灌溉、化工生產(chǎn)等領(lǐng)域。由于低比轉(zhuǎn)數(shù)離心泵的空間非對稱性結(jié)構(gòu)、葉輪和蝸殼的動靜干涉作用、偏離設(shè)計(jì)工況運(yùn)行等原因,造成其內(nèi)部流動呈現(xiàn)出復(fù)雜的非定常特性,導(dǎo)致泵內(nèi)部流體的壓力隨時(shí)間快速波動。壓力脈動是離心泵產(chǎn)生振動和噪聲的重要原因,因此研究低比轉(zhuǎn)數(shù)離心泵運(yùn)行時(shí)的壓力脈動對提高泵的穩(wěn)定性和可靠性有著重要作用[1-4]。目前,研究離心泵壓力脈動的方法主要是數(shù)值模擬與試驗(yàn),國內(nèi)外學(xué)者的研究主要集中在壓力脈動的起因、壓力脈動特性分析和壓力脈動與振動噪聲之間的關(guān)系等。Khalifa等通過試驗(yàn)方法測量了雙吸泵內(nèi)部壓力脈動,研究了壓力脈動與流動誘導(dǎo)振動之間的關(guān)系,并指出優(yōu)化葉輪與隔舌間的間隙能有效降低流動誘導(dǎo)振動[5-6]。Barrio等采用數(shù)值模擬的方法計(jì)算了單級單吸離心泵的非定常流動,葉輪出口流動的脈動頻率以葉頻為主[7]。Benra等研究了多級泵第一級葉輪、導(dǎo)葉內(nèi)的壓力脈動,采用結(jié)構(gòu)網(wǎng)格進(jìn)行全流場非定常計(jì)算,并進(jìn)行時(shí)、頻域分析[8]。Spence等數(shù)值模擬了不同葉輪雙吸泵內(nèi)的非定常不穩(wěn)定流動,并對15個(gè)關(guān)鍵位置進(jìn)行壓力脈動監(jiān)測,分析了各個(gè)位置的壓力脈動特性[9]。施衛(wèi)東等分析了蝸殼基圓大小對超厚葉片離心泵隔舌附近壓力脈動的影響[10]。劉厚林等運(yùn)用FEMBEM聲振耦合計(jì)算和試驗(yàn)測量方法,分析了葉輪出口寬度對離心泵噪聲輻射的影響,發(fā)現(xiàn)葉輪出口寬度存在一個(gè)較好的范圍使得能量性能與外場噪聲等綜合性能較優(yōu)[11]。談明高等分析不同葉輪外徑下,固體顆粒對蝸殼內(nèi)壓力脈動的影響[12]。

    從壓力脈動時(shí)頻域特性中,只能看出壓力值隨時(shí)間的變化過程、或者某一點(diǎn)的壓力脈動頻率,很難得到葉輪和蝸殼內(nèi)每個(gè)網(wǎng)格節(jié)點(diǎn)在整個(gè)旋轉(zhuǎn)周期內(nèi)累積的壓力脈動結(jié)果[13-14]。本文引入壓力脈動標(biāo)準(zhǔn)差對離心泵蝸殼內(nèi)各個(gè)網(wǎng)格節(jié)點(diǎn)在3個(gè)旋轉(zhuǎn)周期內(nèi)非定常壓力脈動強(qiáng)度進(jìn)行表征并用葉輪出口處圓周速度U2進(jìn)行無量綱處理[15-19],對比研究了葉輪出口寬度對各個(gè)工況下蝸殼內(nèi)壓力脈動強(qiáng)度的影響,以期為低比轉(zhuǎn)數(shù)離心泵的設(shè)計(jì)和提高運(yùn)行可靠性提供重要的理論支撐。

    1 模型泵與數(shù)值模擬

    1.1 模型泵

    本文以單級單吸低比轉(zhuǎn)數(shù)離心泵為研究對象,其結(jié)構(gòu)形式為懸臂式,設(shè)計(jì)參數(shù)為:流量QDES=6.3 m3/h,揚(yáng)程HDES=8.0 m,轉(zhuǎn)速n=1 450 r/min,比轉(zhuǎn)數(shù)ns=46.5。主要幾何參數(shù)為:葉輪進(jìn)口直徑Dj=50 mm,葉輪出口直徑D2=160 mm,葉片進(jìn)口安放角β1=24°,葉片出口安放角β2=29°,葉片包角φ=156°,葉片數(shù)Z=6,蝸殼基圓直徑D3=170 mm,蝸殼進(jìn)口寬度b3=18 mm,隔舌安放角α3=27°,蝸殼出口直徑寬度D4=32 mm,蝸殼截面形狀為矩形。為了研究葉輪出口寬度(如圖1(a)所示)對離心泵內(nèi)部非定常壓力脈動強(qiáng)度的影響,將葉輪出口寬度b2設(shè)計(jì)為5、6和7 mm(蝸殼進(jìn)口寬度與葉輪出口寬度之比b3/b2分別為3.6、3和2.6),分別與原蝸殼匹配,并對比分析3個(gè)不同葉輪出口寬度在0.7QDES、1.0QDES、1.3QDES工況下的蝸殼內(nèi)壓力脈動強(qiáng)度分布。

    1.2 數(shù)值模擬

    運(yùn)用Creo Parameter三維造型軟件對葉輪和蝸殼中的流域進(jìn)行造型,在葉輪進(jìn)口處增加了長度為5倍葉輪進(jìn)口直徑的進(jìn)口段,并在蝸殼出口增加了出口段,長度為5倍蝸殼出口直徑,其目的在于減小設(shè)置的邊界條件對泵內(nèi)部流動的影響,保證計(jì)算精度。采用ICEM CFD軟件對整個(gè)計(jì)算域進(jìn)行多塊結(jié)構(gòu)化網(wǎng)格劃分,對整個(gè)流域的壁面進(jìn)行邊界層網(wǎng)格加密,并進(jìn)行了網(wǎng)格無關(guān)性分析,最終確定3個(gè)葉輪的網(wǎng)格數(shù)為1 500 000左右,蝸殼網(wǎng)格數(shù)為1 165 264,進(jìn)口段網(wǎng)格數(shù)為561 235,圖1(b)是結(jié)構(gòu)網(wǎng)格示意圖。應(yīng)用ANSYS CFX 14.5商業(yè)軟件對離心泵進(jìn)行定常和非定常計(jì)算,湍流模型選擇SSTk-ω模型。

    圖1 低比轉(zhuǎn)數(shù)離心泵模型及結(jié)構(gòu)網(wǎng)格Fig.1 Model and structured grids of low-specific-speed centrifugal pump

    在定常計(jì)算設(shè)置中,進(jìn)、出口邊界條件設(shè)置為總壓進(jìn)口和質(zhì)量流量出口,動、靜計(jì)算域的交界面選擇Frozen rotor,收斂殘差RMS設(shè)置為1×10-5,最大迭代步數(shù)為1 000,收斂結(jié)果較好。非定常計(jì)算以定常計(jì)算結(jié)果為初始值,邊界條件設(shè)置與定常計(jì)算設(shè)置相同,只在動、靜計(jì)算域之間的交界面改成Transient frozen rotor,計(jì)算16圈,前6圈以3°為1個(gè)時(shí)間步長,在后10圈中以1°為1個(gè)時(shí)間步長,在分析蝸殼中壓力脈動強(qiáng)度時(shí)采用非定常計(jì)算的最后3個(gè)周期結(jié)果,壓力脈動強(qiáng)度系數(shù)的定義式為

    (1)

    為了準(zhǔn)確地分析數(shù)值模擬結(jié)果,定義了描述位置的變量,如圖1(b)所示,首先定義了x-y直角坐標(biāo)系,點(diǎn)m為蝸殼周向上任意位置,θ表示ox至om的角度,逆時(shí)針為正。

    2 試驗(yàn)驗(yàn)證

    為了保證數(shù)值模擬的準(zhǔn)確性,將葉輪出口寬度6 mm模型泵的外特性、粒子圖像測速(particle image velocimetry, PIV)實(shí)驗(yàn)結(jié)果與定常計(jì)算結(jié)果進(jìn)行對比,在具有二級精度的開式試驗(yàn)臺上進(jìn)行試驗(yàn),試驗(yàn)臺如圖2所示。由于數(shù)值模擬計(jì)算的是水力效率,而試驗(yàn)臺測試的是泵系統(tǒng)的總效率,2個(gè)效率沒有可比性。因此從揚(yáng)程的角度驗(yàn)證數(shù)值模擬的準(zhǔn)確性,如圖3所示,在0.7QDES~1.4QDES工作范圍內(nèi),葉輪出口寬度6 mm模型泵的模擬揚(yáng)程與試驗(yàn)具有很高的一致性,設(shè)計(jì)工況下的揚(yáng)程偏差為3.5%,小流量0.7QDES工況下偏差3.8%,大流量工況1.3QDES相差2.5%。從圖3可以看出,泵的揚(yáng)程隨著葉輪出口寬度的增大而升高,主要是因?yàn)槌隹诘脑龃髮p小工作介質(zhì)的軸面速度,從而提高理論揚(yáng)程。3個(gè)方案最佳效率點(diǎn)的流量均大于設(shè)計(jì)流量,設(shè)計(jì)工況下的水力效率隨著出口葉輪出口寬度的增加而升高,增大葉輪出口寬度也體現(xiàn)了低比轉(zhuǎn)數(shù)離心泵加大流量設(shè)計(jì)法的基本思路[20-23]。該低比轉(zhuǎn)數(shù)離心泵的設(shè)計(jì)揚(yáng)程為8 m,葉輪出口寬度為6 mm和7 mm模型泵的模擬揚(yáng)程達(dá)到了9.2 m,如果再增大葉輪出口寬度會導(dǎo)致?lián)P程進(jìn)一步增大,導(dǎo)致水頭過剩、浪費(fèi)能源,而葉輪出口寬度5 mm模型泵的揚(yáng)程為8.5 mm,如果再減小葉輪出口寬度,會導(dǎo)致水頭低于設(shè)計(jì)要求,因此滿足揚(yáng)程設(shè)計(jì)要求的葉輪出口寬度區(qū)間是5~7 mm,3個(gè)方案的揚(yáng)程均滿足設(shè)計(jì)要求;葉輪出口寬度為6和7 mm的水力效率比較接近,而葉輪出口寬度為5 mm的水力效率相對較低。采用PIV技術(shù)拍攝葉輪出口寬度6 mm模型泵蝸殼中的內(nèi)部流動[24-25],選取蝸殼內(nèi)r=91 mm、區(qū)間為[-20°,50°]這段70°圓弧上的速度作為試驗(yàn)、模擬的對比對象,計(jì)算這段圓弧上的速度平均值并用葉輪出口的圓周速度U2對平均線速度進(jìn)行無量綱化處理[26],無量綱線平均速度如圖4所示。在3個(gè)工況下,PIV試驗(yàn)與數(shù)值模擬的平均線速度誤差均在8.5%以內(nèi),一定程度上驗(yàn)證了數(shù)值模擬的正確性。

    圖2 外特性及粒子圖像測速試驗(yàn)臺Fig.2 External characteristics and PIV test rig

    圖3 不同葉輪出口寬度方案的性能曲線Fig.3 Performance curves of three different impeller outlet width schemes

    圖4 葉輪出口寬度6 mm蝸殼中的粒子圖像測速試驗(yàn)與數(shù)值模擬結(jié)果對比Fig.4 Comparisons between PIV test and numerical simulation results in volute of impeller outlet width 6 mm

    3 結(jié)果與分析

    圖5為不同葉輪出口寬度情況下蝸殼中開面上的壓力脈動強(qiáng)度分布,由圖5可知:在小流量0.7QDES工況下,從隔舌(72°)~150°內(nèi),不同葉輪出口寬度的壓力脈動強(qiáng)度均較大并呈現(xiàn)出大梯度變化特征,而在150°~180°減弱,在180°~200°又再次變大,在200°~隔舌(72°)壓力脈動強(qiáng)度出現(xiàn)反復(fù)、不規(guī)則變化,主要原因是不同初相位、不同頻率的脈動源在各個(gè)網(wǎng)格節(jié)點(diǎn)的壓力波疊加,使得某些區(qū)域的脈動幅值加強(qiáng),某些區(qū)域的脈動幅值減弱。在設(shè)計(jì)工況1.0QDES下,葉輪出口寬度為5和6 mm時(shí),蝸殼中開面上的壓力脈動強(qiáng)度分布較相似,并在相近的位置出現(xiàn)變化;但是在隔舌(72°)~85°內(nèi),葉輪出口寬度5 mm的壓力脈動強(qiáng)度較大,而葉輪出口寬度為6 mm時(shí)則較弱,葉輪出口寬度為7 mm時(shí)在整個(gè)蝸殼中開面上的壓力脈動強(qiáng)度的波動較小,主要是在150°處出現(xiàn)梯度變化,其他方位變化不明顯。在大流量1.3QDES工況下,葉輪出口寬度為5 mm時(shí),在隔舌(72°)~105°區(qū)間內(nèi)出現(xiàn)較大的壓力脈動強(qiáng)度,隨著出口寬度的增大,這個(gè)區(qū)間大小逐漸減小,說明葉輪和隔舌的相互作用所引起的壓力脈動在減弱;葉輪出口寬度為5和7 mm的壓力脈動強(qiáng)度分布特征較相似,其幅值和梯度較接近,而葉輪出口寬度6 mm則以較弱的壓力脈動強(qiáng)度分布為主。在各個(gè)工況下,蝸殼擴(kuò)散段上的壓力脈動強(qiáng)度均較大,主要是因?yàn)閿U(kuò)散段的過流斷面從矩形變化到圓形,過流斷面面積變化梯度大,造成動能到壓能的劇烈變化,而擴(kuò)散段上較大的壓力脈動強(qiáng)度反應(yīng)了能量劇烈轉(zhuǎn)化的過程。對于同一運(yùn)行工況,葉輪出口寬度6 mm模型泵蝸殼中開面上的平均壓力脈動強(qiáng)度最小,在小流量工況下,葉輪出口寬度7 mm的平均壓力脈動強(qiáng)度明顯大于葉輪出口寬度5 mm,而在大流量工況下,兩者比較接近。對于同一葉輪出口寬度,小流量下的平均壓力脈動強(qiáng)度最大,大流量次之,而設(shè)計(jì)工況下的平均壓力脈動強(qiáng)度最小。

    圖5 不同工況下3個(gè)方案蝸殼中開面上的壓力脈動強(qiáng)度對比(QDES=6.3 m3·h-1)Fig.5 Comparisons of pressure fluctuation intensity distribution on center section of three schemes under multiple operating conditions(QDES=6.3 m3·h-1)

    圖6為不同葉輪出口寬度情況下蝸殼72°(隔舌)處流道截面的壓力脈動強(qiáng)度分布。從圖中可以看到,葉輪出口寬度為7 mm方案在小流量0.7QDES工況下的壓力脈動強(qiáng)度的幅值最大,且壓力脈動強(qiáng)度的梯度也是最大的,而葉輪出口寬度6 mm在設(shè)計(jì)流量1.0QDES工況下的壓力脈動強(qiáng)度是最小的,整個(gè)流道截面分布均勻,僅在隔舌附近存在較小的壓力脈動強(qiáng)度梯度。在不同葉輪出口寬度情況下,隔舌處出現(xiàn)的大強(qiáng)度壓力波動沿著徑向擴(kuò)散并不斷減小直至平穩(wěn)。隨著流量的增大,葉輪出口寬度7 mm的壓力脈動強(qiáng)度分布趨于對稱,而葉輪出口寬度5 mm和6 mm在各個(gè)工況下均比較對稱。在大流量工況下,葉輪出口寬度對壓力脈動強(qiáng)度的影響較小,分布較相似。對于同一葉輪出口寬度,設(shè)計(jì)流量下的平均壓力脈動強(qiáng)度最小,偏工況下的平均壓力脈動強(qiáng)度均較大。在同一流量下,葉輪出口寬度6 mm的隔舌處流道截面上的平均壓力脈動強(qiáng)度最小,而葉輪出口寬度5 mm和7 mm在小流量下差異較大,大流量下較接近。

    圖6 不同工況下3個(gè)方案蝸殼72°截面上的壓力脈動強(qiáng)度對比Fig.6 Comparisons of pressure fluctuation intensity distribution on tongue section of three schemes under multiple operating conditions

    低比轉(zhuǎn)數(shù)離心泵葉輪出口的射流-尾跡與非對稱對數(shù)螺旋線型蝸殼中的流動相互影響,造成了葉輪和蝸殼交界面非定常、不均勻的壓力分布,圖7為不同葉輪出口寬度情況下的蝸殼進(jìn)口圓周面上的壓力脈動強(qiáng)度分布對比圖。在小流量0.7QDES工況下,不同的葉輪出口寬度在72°(隔舌)~120°內(nèi)均存在較大的壓力脈動強(qiáng)度,但其壓力脈動強(qiáng)度的梯度較小,葉輪出口寬度7 mm在180°區(qū)域再次出現(xiàn)劇烈的壓力脈動,而葉輪出口寬度5和6 mm在180°區(qū)域內(nèi)的壓力脈動較為平穩(wěn)。葉輪出口寬度7 mm在小流量0.7QDES工況下壓力脈動強(qiáng)度分布呈現(xiàn)出不均勻性,分布紊亂,在其他工況下,在中開線兩側(cè)壓力脈動強(qiáng)度分布比較對稱,而葉輪出口寬度5和6 mm在各個(gè)工況下均比較對稱。在設(shè)計(jì)工況1.0QDES和大流量1.3QDES工況下,3個(gè)不同葉輪出口寬度的壓力脈動強(qiáng)度分布比較接近,在相近的位置出現(xiàn)壓力脈動強(qiáng)度變化,變化幅值也較小。而且在各個(gè)工況下,均在72°(隔舌)出現(xiàn)壓力脈動強(qiáng)度梯度的極大值,劇烈的壓力波動也造成了下游出現(xiàn)較大壓力脈動強(qiáng)度分布??傮w而言,在各個(gè)工況下,葉輪出口寬度6 mm的壓力脈動強(qiáng)度更小,分布更均勻。

    圖7 不同工況下3個(gè)方案蝸殼進(jìn)口圓周面壓力脈動強(qiáng)度對比Fig.7 Comparisons of pressure fluctuation intensity distribution on volute inlet section of three schemes under multiple operating conditions

    4 結(jié)論

    1)沿著蝸殼內(nèi)的主流方向,壓力脈動強(qiáng)度出現(xiàn)強(qiáng)弱交替、不規(guī)則變化。由于葉輪和蝸殼的動靜干涉作用,導(dǎo)致了隔舌部位的壓力脈動強(qiáng)度出現(xiàn)大梯度變化,并在隔舌下游區(qū)域出現(xiàn)壓力脈動強(qiáng)度極大值。在蝸殼擴(kuò)散段,壓力脈動強(qiáng)度較大,但分布比較均勻。

    2)在相同的葉輪出口寬度下,小流量工況下的蝸殼內(nèi)壓力脈動強(qiáng)度及其梯度均是最大的,大流量次之,而設(shè)計(jì)工況下壓力脈動強(qiáng)度及其梯度較小。

    3)低比轉(zhuǎn)數(shù)離心泵葉輪出口寬度存在一個(gè)合適的取值范圍,使得各個(gè)工況下蝸殼內(nèi)的壓力脈動強(qiáng)度較小。綜合考慮該低比轉(zhuǎn)數(shù)離心泵的性能與壓力脈動特性,蝸殼進(jìn)口寬度與葉輪出口寬度之比約為3時(shí),其綜合性能較優(yōu)。

    [1]袁壽其,施衛(wèi)東,劉厚林,等.泵理論與技術(shù)[M]. 北京:機(jī)械工業(yè)出版社,2014.

    [2]楊軍虎.低比轉(zhuǎn)數(shù)離心泵葉輪內(nèi)的流動機(jī)理和葉輪設(shè)計(jì)[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2002,33(2): 47-49. YANG Junhu. Flow mechanism inside the impeller of a centrifugal pump with low specific speed and design of the impeller [J]. Transactions of the Chinese society of agricultural machinery, 2002, 33(2): 47-49.

    [3]王維軍,王洋,劉洋,等.低比轉(zhuǎn)數(shù)變曲率葉型離心泵非定常流動分析[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2013,44(11): 61-66. WANG Weijun, WANG Yang, LIU Yang, et al.Analysis on unsteady flow in low specific speed centrifugal pump with curvature blade [J]. Transactions of the Chinese society of agricultural machinery, 2013, 44(11): 61-66.

    [4]KELDER J D H, DIJKERS R J H, VAN ESCH B P M, et al. Experimental and theoretical study of the flow in the volute of a low specific-speed pump [J]. Fluid dynamics research, 2001, 28(4): 267-280.

    [5]KHALIFA A E, ALQUTUB A M, BEN-MANSOUR R. Study of pressure fluctuations and induced vibration at blade-passing frequencies of a double volute pump [J]. Arabian journal for science and engineering, 2011, 36(7): 1333-1345.

    [6]KHALIFA A E. Performance and vibration of a double volute centrifugal pump: effect of impeller trimming[C]//ASME 2014 International Mechanical Engineering Congress and Exposition. American Society of Mechanical Engineers. Montreal, Quebec, Canada, 2014.

    [7]BARRIO R, PARRONDO J, BLANCO E. Numerical analysis of the unsteady flow in the near-tongue region in a volute-type centrifugal pump for different operating points [J]. Computers & fluids, 2010, 39(5): 859-870.

    [8]BENRA F K, DOHMEN H J. Numerical investigation of the transient flow in a centrifugal pump stage[C]//ASME 2005 Fluids Engineering Division Summer Meeting. American Society of Mechanical Engineers. Houston, Texas, USA, 2005: 71-76.

    [9]SPENCE R, AMARAL-TEIXEIRA J. Investigation into pressure pulsations in a centrifugal pump using numerical methods supported by industrial tests [J]. Computers & fluids, 2008, 37(6): 690-704.

    [10]施衛(wèi)東,張磊,陳斌,等.離心泵間隙對壓力脈動及徑向力的影響[J].排灌機(jī)械工程學(xué)報(bào),2012,30(3): 260-264. SHI Weidong, ZHANG Lei, CHEN Bin, et al.Influence of gap on pressure pulsation and radial force of centrifugal pumps [J]. Journal of drainage and irrigation machinery engineering, 2012, (30)3: 260-264.

    [11]劉厚林,丁劍,談明高,等.葉輪出口寬度對離心泵噪聲輻射影響的分析與試驗(yàn)[J].農(nóng)業(yè)工程學(xué)報(bào),2013,29(16): 66-72. LIU Houlin, DING Jian, TAN Minggao, et al. Analysis and experimental of centrifugal pump noise based on outlet width of impeller[J]. Transactions of the Chinese society of agricultural engineering, 2013, 29(16): 66-72.

    [12]談明高,丁榮,劉厚林,等.不同葉輪外徑下雙流道污水泵壓力脈動特性分析[J].農(nóng)業(yè)工程學(xué)報(bào),2015,31(13): 53-59. TAN Minggao, DING Rong, LIU Houlin, et al. Analysis on pressure pulsation under different impeller diameters in double channel sewage pump[J]. Transactions of the Chinese society of agricultural engineering, 2015, 31(13): 53-59.

    [13]WANG H, TSUKAMOTO H. Fundamental analysis on rotor-stator interaction in a diffuser pump by vertex method [J]. Journal of fluids engineering, 2001, 123(4): 737-747.

    [14]KELDER J D H, DIJKERS R J H, ESCH B P M VAN, et al. Experimental and theoretical study of the flow in the volute of a low specific pump [J]. Fluid dynamics research, 2001, 28(4): 267-280.

    [15]SPENCE R, AMARAL-TEIXEIRA J. A CFD parametric study of geometrical variations on the pressure pulsations and performance characteristics of a centrifugal pump [J]. Computers & fluids, 2009, 38(6): 1243-1257.

    [16]GAETANI P, BOCCAZZI A, SALA R. Low field in the vaned diffuser of a centrifugal pump at different vane setting angles [J]. Journal of fluids engineering, 2012, 134(3):031101.

    [17]袁建平,朱鈺雯,周幫倫,等.余熱排出泵內(nèi)部壓力脈動特性分析[J].排灌機(jī)械工程學(xué)報(bào),2015,33(6): 475-480. YUAN Jianping, ZHU Yuwen, ZHOU Banglun, et al.Pressure fluctuation in residual heat removal pumps [J]. Journal of drainage and irrigation machinery engineering, 2015, 33(6): 475-480.

    [18]王洋,彭帥,劉瑞華,等.旋渦自吸泵內(nèi)部流場壓力脈動數(shù)值模擬[J].排灌機(jī)械工程學(xué)報(bào),2015,33(7): 583-588. WANG Yang, PENG Shuai, LIU Ruihua, et al.Numerical simulation of pressure fluctuation in self-priming vortex pump [J]. Journal of drainage and irrigation machinery engineering, 2015, 33(7): 583-588.

    [19]LUCIUS A, BRENNER G. Numerical simulation and evaluation of velocity fluctuations during rotating stall of a centrifugal pump [J]. Journal of fluids engineering, 2011, 133(8):081102.

    [20]袁壽其,低比速離心泵理論與設(shè)計(jì)[M].北京:機(jī)械工業(yè)出版社,1997.

    [21]王洋,何文俊.基于Fluent的無過載離心泵改型設(shè)計(jì)[J].農(nóng)業(yè)機(jī)械學(xué)報(bào),2009,40(9): 85-88. WANG Yang, HE Wenjun. Improved attempt of non-overload centrifugal pumps based on fluent[J]. Transactions of the Chinese society of agricultural machinery, 2009, 40(9): 85-88.

    [22]倪永燕,袁壽其,袁建平,等.低比轉(zhuǎn)速離心泵加大流量設(shè)計(jì)模型[J].排灌機(jī)械,2008,26(1): 21-24. NI Yongyan, YUAN Shouqi, YUAN Jianping, et al.Model of enlarged flow design for low specific speed centrifugal pump [J]. Drainage and irrigation mchinery, 2008, 26(1): 21-24.

    [23]GüLICH J F. Centrifugal pumps [M]. Berlin: Springer, 2008.

    [24]WESTRA R W, BROERSMA L, VAN ANDEL K, et al. PIV measurements and CFD computations of secondary flow in a centrifugal pump impeller [J]. Journal of fluids engineering, 2010, 132(6):061104.

    [25]WUIBAUT G, BOIS G, DUPONT P, et al. PIV measurements in the impeller and the vaneless diffuser of a radial flow pump in design and off-design operating conditions [J]. Journal of fluids engineering, 2002, 124(3): 791-797.

    [26]冒杰云.低比速離心泵小流量工況內(nèi)部非穩(wěn)態(tài)流動特性研究[D]. 鎮(zhèn)江:江蘇大學(xué), 2015. MAO Jieyun. Research on Unstable flow in a low-specific speed centrifugal pump at partload conditions[D]. Zhenjiang: Jiangsu University, 2015.

    本文引用格式:

    顧延?xùn)|,袁壽其,裴吉,等. 泵葉輪出口寬度對蝸殼內(nèi)壓力脈動強(qiáng)度的影響[J]. 哈爾濱工程大學(xué)學(xué)報(bào), 2017, 38(7): 1023-1030.

    GU Yandong, YUAN Shouqi, PEI Ji, et al. Effects of the outlet width of pump impeller on pressure fluctuation intensity in volute[J]. Journal of Harbin Engineering University, 2017, 38(7): 1023 -1030.

    Effects of the outlet width of pump impeller on pressure fluctuation intensity in volute

    GU Yandong, YUAN Shouqi, PEI Ji, ZHANG Jinfeng, HUANG Xi, WANG Wenjie

    (National Research Center of Water Pump And System Engineering Technology, Jiangsu University, Zhenjiang 212013, China)

    To investigate the effects of impeller outlet width on pressure fluctuation intensity in the volute of a pump with low-specific-speed, the impeller outlet widths were designed as 5, 6, and 7 mm. Other geometric parameters, including impeller and volute, were unchanged. The commercial software ANSYS CFX 14.5 was applied to solve the steady and unsteady Reynolds-averaged Navier-Stokes equations. The head curves and dimensionless mean velocity in the volute calculated by the steady numerical simulations displayed good agreement with those obtained by external characteristics and PIV experiments; therefore, the numerical simulation results were reliable to a certain extent. Furthermore, the external characteristics of three schemes met the design requirements. In the unsteady simulations, a dimensionless standard deviation coefficient of pressure fluctuation for three entire revolution periods was defined as the pressure fluctuation intensity. For the same impeller outlet width, the magnitude and gradient of pressure fluctuation intensity in the volute were the largest under the partial load operating condition, moderate under the overload operating condition, and minimum under the design operating condition. The gradient of pressure fluctuation intensity in the tongue region was larger than that in any other region because of the rotor-stator interaction. When the ratio of the volute inlet width to impeller outlet width was around 3, the performance and operating reliability of this low-specific-speed pump were in good balance, and the fluctuating intensity was low and comprehensive performance was high.

    low specific speed; centrifugal pump; numerical simulation; pressure fluctuation; multiple operating conditions; impeller; head; volute

    2016-08-16.

    日期:2017-04-28.

    中國博士后科學(xué)基金項(xiàng)目(2015T80507);江蘇高校優(yōu)勢學(xué)科建設(shè)工程項(xiàng)目(PAPD).

    顧延?xùn)|(1989-), 男, 博士研究生; 袁壽其(1963-), 男, 研究員, 博士生導(dǎo)師.

    袁壽其,E-mail:shouqiy@ujs.edu.cn.

    10.11990/jheu.201608030

    TH311

    A

    1006-7043(2017)07-1023-08

    網(wǎng)絡(luò)出版地址:http://www.cnki.net/kcms/detail/23.1390.u.20170428.0803.006.html

    猜你喜歡
    蝸殼離心泵脈動
    新學(xué)期,如何“脈動回來”?
    家教世界(2023年25期)2023-10-09 02:11:56
    RBI在超期服役脈動真空滅菌器定檢中的應(yīng)用
    焊接蝸殼泵模型開發(fā)
    一種改善離心泵運(yùn)行狀態(tài)的方法
    大型立式單級引黃離心泵軸向力平衡的研究
    離心泵流場外特性及激勵(lì)數(shù)值計(jì)算
    高比速離心泵的高效設(shè)計(jì)研究
    污水泵蝸殼斷裂原因分析
    地球脈動(第一季)
    蝸殼差壓法在印尼ASAHAN一級水電站的應(yīng)用
    波多野结衣巨乳人妻| 99久久99久久久精品蜜桃| 免费av观看视频| 精品一区二区三区视频在线观看免费| 久久这里只有精品中国| 露出奶头的视频| 久久性视频一级片| 国产aⅴ精品一区二区三区波| 国产精品国产高清国产av| 亚洲真实伦在线观看| 99在线人妻在线中文字幕| 日本 欧美在线| 99久国产av精品| 韩国av一区二区三区四区| 91久久精品电影网| 国产单亲对白刺激| а√天堂www在线а√下载| 别揉我奶头~嗯~啊~动态视频| 精品人妻一区二区三区麻豆 | 婷婷丁香在线五月| 少妇人妻精品综合一区二区 | 一区二区三区四区激情视频 | 搡老岳熟女国产| 国产高清激情床上av| 99热只有精品国产| 丰满乱子伦码专区| 亚洲无线在线观看| 欧美成人性av电影在线观看| 久久久久九九精品影院| 男插女下体视频免费在线播放| 亚洲成人久久爱视频| 免费看日本二区| 丰满人妻熟妇乱又伦精品不卡| 看免费av毛片| 黄色配什么色好看| 国产视频一区二区在线看| 在线观看一区二区三区| 搡老妇女老女人老熟妇| 午夜福利成人在线免费观看| 精品福利观看| 欧美午夜高清在线| 亚洲av成人不卡在线观看播放网| 欧美高清成人免费视频www| 12—13女人毛片做爰片一| 99热只有精品国产| 日韩有码中文字幕| 真人做人爱边吃奶动态| 欧美乱色亚洲激情| 日韩欧美在线乱码| 波多野结衣巨乳人妻| 欧美日韩国产亚洲二区| 日本成人三级电影网站| 麻豆久久精品国产亚洲av| 国产成+人综合+亚洲专区| 怎么达到女性高潮| 欧美激情国产日韩精品一区| 在线观看免费视频日本深夜| 91午夜精品亚洲一区二区三区 | 一a级毛片在线观看| 香蕉av资源在线| 久久久久久久午夜电影| 欧美激情在线99| 真人一进一出gif抽搐免费| 国产精品伦人一区二区| 国产综合懂色| 中文资源天堂在线| 女同久久另类99精品国产91| 午夜福利在线在线| 国产精品一及| 极品教师在线免费播放| 超碰av人人做人人爽久久| 99在线人妻在线中文字幕| 日韩欧美国产一区二区入口| 成人特级黄色片久久久久久久| 国产精品一区二区性色av| 亚洲 国产 在线| 国产精品三级大全| 2021天堂中文幕一二区在线观| 麻豆国产97在线/欧美| 国产欧美日韩一区二区精品| 九色国产91popny在线| 99久久成人亚洲精品观看| 亚洲av美国av| 久久久成人免费电影| 久久久精品欧美日韩精品| 日本a在线网址| 国内毛片毛片毛片毛片毛片| 欧美日韩乱码在线| 乱人视频在线观看| 午夜福利视频1000在线观看| 免费观看的影片在线观看| 亚洲久久久久久中文字幕| 亚洲男人的天堂狠狠| 国产在线男女| 一区二区三区激情视频| 少妇被粗大猛烈的视频| 精品久久久久久,| 天堂动漫精品| 亚洲精华国产精华精| 成年免费大片在线观看| 午夜久久久久精精品| 国产亚洲精品综合一区在线观看| 亚洲精品一卡2卡三卡4卡5卡| 又紧又爽又黄一区二区| 成人特级黄色片久久久久久久| 麻豆久久精品国产亚洲av| 国产伦人伦偷精品视频| 一区二区三区免费毛片| 欧美日韩福利视频一区二区| 国产精品亚洲av一区麻豆| 国产aⅴ精品一区二区三区波| 好看av亚洲va欧美ⅴa在| 永久网站在线| 亚洲精品成人久久久久久| 国产一区二区在线观看日韩| 国产黄a三级三级三级人| 一夜夜www| 国产淫片久久久久久久久 | 色尼玛亚洲综合影院| 亚洲国产日韩欧美精品在线观看| 九九在线视频观看精品| 午夜免费成人在线视频| 一进一出好大好爽视频| 一个人免费在线观看的高清视频| 一边摸一边抽搐一进一小说| 亚洲一区高清亚洲精品| 一本精品99久久精品77| 亚洲午夜理论影院| 在线天堂最新版资源| 午夜福利18| 亚洲av美国av| 免费观看人在逋| 18禁黄网站禁片午夜丰满| 亚洲第一电影网av| 91九色精品人成在线观看| 国产精品亚洲一级av第二区| 日本与韩国留学比较| 亚洲18禁久久av| 日本a在线网址| 三级毛片av免费| 国产av一区在线观看免费| 好看av亚洲va欧美ⅴa在| 日韩欧美 国产精品| 一进一出好大好爽视频| 69av精品久久久久久| 亚洲成人久久爱视频| av在线老鸭窝| 欧美激情久久久久久爽电影| 国内精品美女久久久久久| 日本 av在线| 国产一区二区三区在线臀色熟女| 亚洲五月婷婷丁香| 亚洲avbb在线观看| 亚洲午夜理论影院| 久久久久久九九精品二区国产| 国产三级在线视频| 久久久久久久午夜电影| 久久久国产成人免费| 久久国产精品人妻蜜桃| 757午夜福利合集在线观看| 国产三级在线视频| 免费看美女性在线毛片视频| 亚洲av.av天堂| 超碰av人人做人人爽久久| 国产精品女同一区二区软件 | 成人无遮挡网站| xxxwww97欧美| 日日夜夜操网爽| 一级黄色大片毛片| 亚洲美女黄片视频| 亚洲熟妇熟女久久| 性插视频无遮挡在线免费观看| 亚洲人成电影免费在线| 国产精品永久免费网站| 欧美激情国产日韩精品一区| 欧美色欧美亚洲另类二区| 人人妻人人看人人澡| 精品一区二区三区视频在线观看免费| 久久婷婷人人爽人人干人人爱| 听说在线观看完整版免费高清| 午夜福利视频1000在线观看| 国产久久久一区二区三区| 好看av亚洲va欧美ⅴa在| 又粗又爽又猛毛片免费看| 亚洲18禁久久av| 色精品久久人妻99蜜桃| 久久草成人影院| 91九色精品人成在线观看| 成人国产一区最新在线观看| 日本精品一区二区三区蜜桃| 午夜福利在线观看吧| 亚洲不卡免费看| 成人亚洲精品av一区二区| 99国产精品一区二区三区| 中文字幕av成人在线电影| 精品久久久久久久久av| 禁无遮挡网站| 亚洲天堂国产精品一区在线| 国产精品1区2区在线观看.| 婷婷六月久久综合丁香| 日韩欧美免费精品| 99热这里只有是精品在线观看 | 长腿黑丝高跟| 欧美成人性av电影在线观看| 欧美中文日本在线观看视频| 757午夜福利合集在线观看| 色噜噜av男人的天堂激情| 高清毛片免费观看视频网站| 一卡2卡三卡四卡精品乱码亚洲| 赤兔流量卡办理| 欧美日韩黄片免| 91麻豆精品激情在线观看国产| 男女视频在线观看网站免费| 亚洲一区高清亚洲精品| 亚洲第一欧美日韩一区二区三区| 人妻丰满熟妇av一区二区三区| 亚洲在线自拍视频| 欧美精品啪啪一区二区三区| 亚洲精品成人久久久久久| 亚洲av一区综合| 国产中年淑女户外野战色| 赤兔流量卡办理| 97碰自拍视频| 90打野战视频偷拍视频| 一进一出抽搐动态| 精品久久久久久久人妻蜜臀av| 欧美日韩黄片免| 精品乱码久久久久久99久播| 十八禁人妻一区二区| 国产在视频线在精品| 波野结衣二区三区在线| 亚洲成人中文字幕在线播放| 一级毛片久久久久久久久女| 男人狂女人下面高潮的视频| 91狼人影院| av天堂在线播放| 午夜亚洲福利在线播放| 99久国产av精品| 成人高潮视频无遮挡免费网站| 少妇丰满av| 一个人免费在线观看电影| 老司机午夜福利在线观看视频| 欧美激情国产日韩精品一区| 男人舔女人下体高潮全视频| 精品久久久久久久末码| 成人亚洲精品av一区二区| 久久国产精品影院| 97超视频在线观看视频| 亚洲人与动物交配视频| 国产精品亚洲美女久久久| 日韩中字成人| 久久亚洲真实| 国产大屁股一区二区在线视频| 全区人妻精品视频| 亚洲熟妇中文字幕五十中出| av女优亚洲男人天堂| 可以在线观看毛片的网站| 日本五十路高清| 岛国在线免费视频观看| 日本三级黄在线观看| 成年版毛片免费区| 免费电影在线观看免费观看| 午夜福利视频1000在线观看| 亚洲欧美日韩高清在线视频| 国产极品精品免费视频能看的| 99热只有精品国产| 好男人在线观看高清免费视频| 欧美最黄视频在线播放免费| 亚洲精品在线观看二区| 成人午夜高清在线视频| 天堂影院成人在线观看| 日本黄大片高清| 国产一区二区在线av高清观看| 一进一出抽搐动态| 色5月婷婷丁香| www日本黄色视频网| 高清毛片免费观看视频网站| 男女床上黄色一级片免费看| 国产精华一区二区三区| 国产探花在线观看一区二区| 中文字幕av成人在线电影| 亚洲熟妇熟女久久| 久久久久九九精品影院| 日韩欧美国产在线观看| 成人毛片a级毛片在线播放| 哪里可以看免费的av片| 简卡轻食公司| 小蜜桃在线观看免费完整版高清| 亚洲精品一区av在线观看| а√天堂www在线а√下载| av黄色大香蕉| 亚洲 欧美 日韩 在线 免费| 淫秽高清视频在线观看| 老司机午夜福利在线观看视频| 在线观看66精品国产| 天堂√8在线中文| 美女高潮的动态| 日韩精品中文字幕看吧| 日韩欧美一区二区三区在线观看| 日韩精品中文字幕看吧| 国产不卡一卡二| 在现免费观看毛片| 日日摸夜夜添夜夜添小说| 18禁黄网站禁片免费观看直播| 国模一区二区三区四区视频| 天堂网av新在线| 99国产极品粉嫩在线观看| 一个人免费在线观看电影| 久久伊人香网站| 两个人的视频大全免费| 欧美日韩中文字幕国产精品一区二区三区| 最新在线观看一区二区三区| 一个人看视频在线观看www免费| 久久久久久久久大av| 一级作爱视频免费观看| 久久久久久久久中文| 搡老岳熟女国产| 日本一二三区视频观看| 波多野结衣巨乳人妻| 亚洲成人中文字幕在线播放| 国产亚洲av嫩草精品影院| 国产精品电影一区二区三区| 国产三级在线视频| 精品久久久久久久久av| 日本a在线网址| 一夜夜www| 午夜精品在线福利| 久久精品国产99精品国产亚洲性色| 午夜精品一区二区三区免费看| 在线十欧美十亚洲十日本专区| 成人av一区二区三区在线看| 国产午夜精品久久久久久一区二区三区 | 欧美xxxx性猛交bbbb| 免费人成在线观看视频色| av在线观看视频网站免费| 亚洲黑人精品在线| 看免费av毛片| 村上凉子中文字幕在线| av视频在线观看入口| 精品久久久久久久久久久久久| 国产高潮美女av| 一夜夜www| 99热这里只有是精品50| 日本与韩国留学比较| 欧美成人一区二区免费高清观看| 精品国内亚洲2022精品成人| 日韩欧美精品v在线| 日本 欧美在线| 观看免费一级毛片| 99久久久亚洲精品蜜臀av| 免费看美女性在线毛片视频| 一卡2卡三卡四卡精品乱码亚洲| 身体一侧抽搐| 一个人免费在线观看电影| 国产色爽女视频免费观看| 久久国产乱子伦精品免费另类| 蜜桃久久精品国产亚洲av| a级毛片a级免费在线| 首页视频小说图片口味搜索| 免费av不卡在线播放| 成人特级av手机在线观看| 校园春色视频在线观看| 在线a可以看的网站| 亚洲内射少妇av| 特大巨黑吊av在线直播| av在线蜜桃| 午夜亚洲福利在线播放| 国产又黄又爽又无遮挡在线| 亚洲国产精品999在线| 最后的刺客免费高清国语| 看黄色毛片网站| 日韩欧美一区二区三区在线观看| 日韩有码中文字幕| 国产毛片a区久久久久| 亚洲综合色惰| 国产美女午夜福利| 精品福利观看| 丰满乱子伦码专区| 黄色丝袜av网址大全| 国产av在哪里看| 在现免费观看毛片| 伊人久久精品亚洲午夜| 高清日韩中文字幕在线| 亚洲一区高清亚洲精品| 国产精品久久久久久久电影| 国产亚洲欧美在线一区二区| 久久精品国产99精品国产亚洲性色| 夜夜看夜夜爽夜夜摸| 俄罗斯特黄特色一大片| 国产 一区 欧美 日韩| 国产野战对白在线观看| 18禁裸乳无遮挡免费网站照片| 亚洲午夜理论影院| 免费电影在线观看免费观看| 国产精品久久久久久久久免 | 一进一出好大好爽视频| 国产日本99.免费观看| 动漫黄色视频在线观看| 观看美女的网站| 国产激情偷乱视频一区二区| 日韩欧美精品v在线| 午夜精品在线福利| 少妇熟女aⅴ在线视频| 国产精品一区二区三区四区久久| 亚洲精品456在线播放app | 国产精品综合久久久久久久免费| 国产高清有码在线观看视频| 欧美一级a爱片免费观看看| 在线观看一区二区三区| 欧美xxxx性猛交bbbb| 国产乱人伦免费视频| 久久久久久大精品| 天堂影院成人在线观看| 色在线成人网| 在现免费观看毛片| 伦理电影大哥的女人| a级毛片免费高清观看在线播放| www.熟女人妻精品国产| 欧美日本亚洲视频在线播放| 久久午夜亚洲精品久久| 日本熟妇午夜| 久久久久久久久久黄片| 无人区码免费观看不卡| 变态另类成人亚洲欧美熟女| 哪里可以看免费的av片| 别揉我奶头 嗯啊视频| 亚洲成人精品中文字幕电影| 好男人在线观看高清免费视频| 熟妇人妻久久中文字幕3abv| 亚洲精品粉嫩美女一区| 色av中文字幕| 亚洲精品亚洲一区二区| 午夜老司机福利剧场| 欧美绝顶高潮抽搐喷水| 国产一区二区激情短视频| 亚洲精品乱码久久久v下载方式| 又粗又爽又猛毛片免费看| 精品人妻视频免费看| 亚洲欧美日韩卡通动漫| 亚洲天堂国产精品一区在线| 精品人妻1区二区| 成人av在线播放网站| 亚洲电影在线观看av| 淫秽高清视频在线观看| 看免费av毛片| 高清日韩中文字幕在线| 国产久久久一区二区三区| 男女下面进入的视频免费午夜| 亚洲成a人片在线一区二区| 中文字幕av成人在线电影| 成年人黄色毛片网站| 日韩精品中文字幕看吧| 偷拍熟女少妇极品色| 免费人成在线观看视频色| 亚洲av.av天堂| 直男gayav资源| 国产不卡一卡二| 国产欧美日韩一区二区精品| 日本成人三级电影网站| 我要看日韩黄色一级片| 日日摸夜夜添夜夜添小说| 精品午夜福利在线看| 老女人水多毛片| 99久久精品热视频| 人妻久久中文字幕网| 欧美在线黄色| 美女 人体艺术 gogo| 久久精品夜夜夜夜夜久久蜜豆| а√天堂www在线а√下载| 免费av观看视频| 久久久久久国产a免费观看| 国内精品久久久久精免费| 国产精品一及| 国产v大片淫在线免费观看| 美女大奶头视频| 在线看三级毛片| 午夜福利在线在线| 国产美女午夜福利| 99久久久亚洲精品蜜臀av| 日本一二三区视频观看| .国产精品久久| 两个人视频免费观看高清| 久久久久免费精品人妻一区二区| 国产一区二区亚洲精品在线观看| 国内精品久久久久精免费| 亚洲18禁久久av| 亚洲人成网站高清观看| 久久99热这里只有精品18| 国产成人啪精品午夜网站| 亚洲精品成人久久久久久| 最新在线观看一区二区三区| 性色avwww在线观看| 18禁在线播放成人免费| 国产伦精品一区二区三区四那| 亚州av有码| 国产aⅴ精品一区二区三区波| 男女做爰动态图高潮gif福利片| 亚洲av电影不卡..在线观看| 欧美日本亚洲视频在线播放| 欧美色欧美亚洲另类二区| 人人妻,人人澡人人爽秒播| 三级毛片av免费| 日韩欧美国产一区二区入口| 51国产日韩欧美| 久久久久久九九精品二区国产| 成年人黄色毛片网站| 国产午夜福利久久久久久| 免费观看的影片在线观看| 亚洲电影在线观看av| 一本久久中文字幕| 午夜两性在线视频| 日韩高清综合在线| 精品不卡国产一区二区三区| 91在线精品国自产拍蜜月| 色5月婷婷丁香| 国产一区二区三区在线臀色熟女| 国产激情偷乱视频一区二区| 女生性感内裤真人,穿戴方法视频| 又黄又爽又刺激的免费视频.| 日韩欧美三级三区| 十八禁国产超污无遮挡网站| 精品国产三级普通话版| 亚洲男人的天堂狠狠| 搡老岳熟女国产| 人妻久久中文字幕网| avwww免费| 欧美高清性xxxxhd video| 男女做爰动态图高潮gif福利片| 欧美极品一区二区三区四区| 国产三级中文精品| 成人av在线播放网站| 又粗又爽又猛毛片免费看| 一级黄片播放器| 亚洲 国产 在线| 在线看三级毛片| 精品人妻视频免费看| 91麻豆av在线| 国产伦精品一区二区三区视频9| 99久久成人亚洲精品观看| 香蕉av资源在线| 天堂网av新在线| 国产黄a三级三级三级人| 国产伦精品一区二区三区四那| 国产一区二区在线av高清观看| 国产欧美日韩精品一区二区| 久久精品综合一区二区三区| 99精品在免费线老司机午夜| 九色成人免费人妻av| 免费一级毛片在线播放高清视频| 九九热线精品视视频播放| 国产伦在线观看视频一区| 欧美激情在线99| 女同久久另类99精品国产91| 美女cb高潮喷水在线观看| 国产爱豆传媒在线观看| 亚洲欧美激情综合另类| 亚洲自拍偷在线| 校园春色视频在线观看| 国产蜜桃级精品一区二区三区| 久久久久精品国产欧美久久久| 国产精品免费一区二区三区在线| 国产 一区 欧美 日韩| 免费高清视频大片| 亚洲在线观看片| 哪里可以看免费的av片| 中文字幕av在线有码专区| 舔av片在线| 欧美精品啪啪一区二区三区| 夜夜爽天天搞| 亚洲成人免费电影在线观看| 一级作爱视频免费观看| 午夜福利在线在线| 亚洲av不卡在线观看| 久久久成人免费电影| 国产色爽女视频免费观看| 国产国拍精品亚洲av在线观看| 免费一级毛片在线播放高清视频| 色综合站精品国产| 老司机深夜福利视频在线观看| 免费搜索国产男女视频| 国产视频一区二区在线看| 最近在线观看免费完整版| 99久久久亚洲精品蜜臀av| 日韩 亚洲 欧美在线| 国产av在哪里看| 成人亚洲精品av一区二区| 亚洲精品456在线播放app | 久久久国产成人免费| 婷婷精品国产亚洲av在线| 欧美3d第一页| 国产精品av视频在线免费观看| 亚洲人成网站高清观看| 琪琪午夜伦伦电影理论片6080| 熟妇人妻久久中文字幕3abv| 欧美成人a在线观看| 夜夜看夜夜爽夜夜摸| 亚洲人与动物交配视频| 精品久久久久久,| 精品久久久久久成人av| 日韩欧美一区二区三区在线观看| 国产在线精品亚洲第一网站| 成人特级av手机在线观看| 亚洲va日本ⅴa欧美va伊人久久| 久久午夜福利片| 久久久久久九九精品二区国产| 国产成人av教育| 美女大奶头视频| 国产在视频线在精品| 长腿黑丝高跟| 亚洲,欧美,日韩| 亚洲黑人精品在线| 不卡一级毛片| 久久99热6这里只有精品| 男女下面进入的视频免费午夜| 国产高清三级在线|