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

    葉片數(shù)對離心泵振動噪聲性能的影響

    2019-02-20 13:47:08談明高陸友東吳澤瑾吳賢芳劉厚林
    農(nóng)業(yè)工程學報 2019年23期
    關鍵詞:蝸殼揚程離心泵

    談明高,陸友東,吳澤瑾,吳賢芳,劉厚林

    葉片數(shù)對離心泵振動噪聲性能的影響

    談明高1,陸友東1,吳澤瑾1,吳賢芳2,劉厚林1

    (1. 江蘇大學流體機械及工程技術研究中心,鎮(zhèn)江 212013; 2. 江蘇大學能源與動力工程學院,鎮(zhèn)江 212013)

    葉片數(shù)是離心泵的主要幾何參數(shù)之一。為研究葉片數(shù)對離心泵振動噪聲性能的影響,以比轉(zhuǎn)速為97的離心泵為例,對比了不同葉片數(shù)下的水力和振動噪聲性能,并采用FEMBEM聲振耦合計算方法對流動激勵下的振動及其聲輻射噪聲進行了數(shù)值模擬,同時與試驗數(shù)據(jù)進行對比分析。結果表明:提出的數(shù)值模擬方法可用于預測泵的流動誘導振動和聲輻射性能,且在模擬中考慮口環(huán)泄漏的影響能夠提高計算精度,有口環(huán)方案預測得到的振幅較無口環(huán)方案的預測精度提高了13.5%。隨著葉片數(shù)的增加,揚程和軸功率均逐漸增大,最大增幅分別為15.9%和14.1%;效率隨葉片數(shù)的增加呈先增大后減小再增大的趨勢。離心泵蝸殼的壓力脈動幅值隨葉片數(shù)的減小而增大。由于葉輪蝸殼動靜干涉的作用,蝸殼隔舌處、第1到第2斷面間和擴壓管壁面等3個區(qū)域的壓力脈動幅值相對較高。隨著葉片數(shù)的減少,蝸殼壁面的振動位移有所增大,最大位移主要發(fā)生蝸殼第8斷面處。振動速度隨著葉片數(shù)的增大后減小,與振動位移的規(guī)律有一定的差異,振動高速區(qū)主要集中在隔舌、蝸殼的第4與第6斷面之間和靠近擴壓管的第8斷面處。設計工況下,泵在葉頻對應的聲壓級和聲強隨著葉片數(shù)的增加先增大后減小,高聲壓級區(qū)域主要出現(xiàn)在泵出口附近的高振動速度引起的垂直方向。綜合考慮水力和振動噪聲性能,確定該模型泵的最佳葉片數(shù)為6。

    離心泵;數(shù)值模擬;壓力;口環(huán)泄漏;葉片數(shù);振動;噪聲

    0 引 言

    離心泵是通過葉輪旋轉(zhuǎn)來驅(qū)動流體的機械,廣泛應用于各種工業(yè)和民用領域[1]。泵在運行過程中產(chǎn)生的振動和噪聲,不僅對設備的使用壽命和系統(tǒng)性能有很大影響[2],而且還會破壞工作環(huán)境,影響人的身心健康。隨著經(jīng)濟和社會的發(fā)展,對泵的振動和噪聲的要求越來越高[3-5]。因此,降低泵振動和噪聲水平了成為研究熱點[6-8]。

    目前,國內(nèi)外對離心泵振動噪聲的研究主要集中在理論[9-10]和試驗研究[11-13]方面。近幾十年來,基于聲學分析的計算流體動力學已取得了一些進展。研究提出了根據(jù)CFD數(shù)值計算[14-15]的方法或者離散渦法[16-17]來預測聲源的方法,然后采用邊界元法(boundary element method,BEM)完成聲輻射計算[15,18-19]。Kato等[20]提出了一種離心泵外表面噪聲預測的流體和結構分析的單向耦合模擬方法?;诖鬁u模擬(large eddy simulation,LES)進行內(nèi)流計算,后運用有限元法計算壓力脈動,最終與測量結果進行對比,研究發(fā)現(xiàn)該方法下葉頻對應的預測振動噪聲和實測吻合較好。針對流體與結構網(wǎng)格的傳輸,Jiang等[21]開發(fā)了一種數(shù)據(jù)接口工具來解決網(wǎng)格匹配的問題,并采用有限元法對結構振動進行了數(shù)值模擬,提取了葉頻對應的振動模態(tài),結果闡明了振動噪聲產(chǎn)生和傳播的機理。

    葉片數(shù)對葉輪流道內(nèi)流體流動的不均勻性有明顯的影響[22]。隨著葉片數(shù)的減少,葉片間距的增大,流體沿圓周流動的不均勻性增加。針對泵性能的不同要求,國內(nèi)外對葉片數(shù)的優(yōu)化[23-24]進行了許多研究,研究表明葉片數(shù)對泵的振動噪聲性能具有較大影響[25-26]。以往的研究大多是在不考慮葉輪和蝸殼中泄漏流場的情況下研究兩者之間的相互作用。事實上,位于定子和轉(zhuǎn)子的流體泄漏對離心泵的能量性能有很大影響,葉輪出口的主流通過間隙流入蝸殼,主流與間隙內(nèi)流體相互作用,引起間隙內(nèi)不均勻的壓力分布。同時,泵內(nèi)泄漏對離心泵的非定常特性也起著重要作用[27-28]。

    本文考慮葉輪和蝸殼間泄漏的影響后,建立了分析離心泵振動和噪聲的預測模型,研究了葉片數(shù)對流體激勵下引起的泵殼振動輻射噪聲的影響,研究結果能夠為離心泵減振降噪設計提供一定的參考。

    1 試驗裝置

    1.1 研究模型

    選取比轉(zhuǎn)速s=97的單級離心泵作為研究對象,其設計參數(shù)d=50 m3/h、轉(zhuǎn)速=2 900 r/min、揚程=30 m。模型泵主要結構參數(shù):葉輪入口直徑1=0.072 m,出口直徑2=0.168m,出口寬度2=0.01 m,葉片出口角2=33°,包角=115°,口環(huán)間隙為0.2 mm。

    1.2 研究方案

    為分析不同葉片數(shù)對離心泵性能的影響,給出了4種不同葉片數(shù)的葉輪,將葉片數(shù)為4、5、6和7的離心泵分別記為泵I、泵II、泵III和泵IV。對圖1所示的葉輪進行快速成型加工,進行性能測試。

    注:z為葉片數(shù)。

    1.3 試驗系統(tǒng)

    離心泵閉式試驗臺如圖2所示,試驗系統(tǒng)主要包括模型泵、渦輪流量計、電機、數(shù)據(jù)采集器、壓阻式壓力傳感器、霍爾傳感器等。泵轉(zhuǎn)矩由傳感器測試,通過計算機采集數(shù)據(jù)。流量計安裝在距出口管1 m處,流量由閥門控制。壓力傳感器測量進出口靜壓,其測量范圍為?100~100 kPa和0~600 kPa。泵的振動加速度由PCB 352A60加速度計測試,安裝在泵進口處,其靈敏度為10 mV/(m/s2)。加速度計的位置如圖3所示。用PXI-6251數(shù)據(jù)采集模塊采集電信號,并由Lab View軟件進行分析。壓力傳感器、渦輪流量計、加速度計和霍爾傳感器的測量不確定度范圍分別為±0.5%、±0.5%、±2%、±1.5%。

    1.真空泵 2.汽蝕筒 3、4、8、10.碟閥 5.渦輪流量計 6.壓力變送器 7.模型泵 9.電機 11.穩(wěn)壓罐 12.球閥

    2 數(shù)值模擬方法

    2.1 三維造型及網(wǎng)格劃分

    研究模型由三維造型軟件Pro/E 5.0生成,并由ICEM進行網(wǎng)格劃分,運用六面體結構化網(wǎng)格劃分確保網(wǎng)格的質(zhì)量,網(wǎng)格域共分為4部分:進口域,葉輪域,蝸殼和出口域,泄露流道域。其中葉輪域與其他3個域之間設置動靜交界面,如圖3所示。

    注:將進口域和泄露流道域間為交界面記為A,泄露流道域和蝸殼及出口域交界面記為F,泄露流道域和葉輪交界面分別交界面B,C,D,E。

    為分析有無口環(huán)泄漏對離心泵噪聲數(shù)值計算的影響,對有無口環(huán)泄漏分別進行了造型,圖4a和4b分別為2種模型的源網(wǎng)格。

    圖4 模型源網(wǎng)格

    2.2 邊界條件設置

    利用CFX軟件進行求解,采用多重參考系模擬葉輪蝸殼的相互作用。葉輪流場處于旋轉(zhuǎn)坐標系中,蝸殼和泄漏流道在靜止坐標系中計算。通過動靜交界面進行數(shù)據(jù)交換,設置一般網(wǎng)格界面,穩(wěn)態(tài)計算采用凍結轉(zhuǎn)子交界面,瞬態(tài)計算采用瞬態(tài)動靜交界面。瞬態(tài)計算初始條件設置為定常計算結果,計算時長為5圈。

    綜合考慮間隙流動的精確性和計算周期,采用SST湍流模型[29],并通過標準壁面函數(shù)計算邊界層變量。相對壓力設為0,計算域中所有表面均采用無滑移壁面條件。收斂殘差設置為10-4。進口條件設為1 atm的恒定總壓,出口條件設為質(zhì)量流量。

    蝸殼材料為鑄鋼,其彈性模量=211 GPa,密度=7 870 kg/m3,泊松比=0.29。蝸殼結構網(wǎng)格如圖5所示。

    注:地腳螺栓孔節(jié)點處的速度為ux= uy= uz=0;進出口法蘭速度為uz=0和ux=0;軸承孔上的節(jié)點速度為ux= uy=uz=0。

    2.3 網(wǎng)格相關性驗證

    通過5套不同數(shù)量的網(wǎng)格檢查網(wǎng)格相關性,分別記為方案A、方案B、方案C、方案D和方案E,如表2所示。根據(jù)不同網(wǎng)格數(shù)下泵揚程的變化確定最佳網(wǎng)格數(shù)量,如圖6所示。隨著網(wǎng)格數(shù)量的增加,泵揚程系數(shù)逐漸接近恒定值,揚程系數(shù)為該網(wǎng)格數(shù)下計算揚程與設計揚程的比值,因此最終選取網(wǎng)格D進行研究分析。

    表1 5種不同網(wǎng)格方案

    圖6 網(wǎng)格數(shù)對揚程的影響

    3 蝸殼結構聲學模擬

    3.1 蝸殼振動模擬方法

    結構動態(tài)響應和周圍空氣中輻射聲壓的控制方程如下

    式中[]為質(zhì)量矩陣;[]為阻尼矩陣;[]為剛度矩陣;{}為節(jié)點結構位移矢量;{P(t)}為施加在節(jié)點結構上的外部激勵力矢量,由CFD計算得到。

    阻尼矩陣[]由Rayleigh理論給出,包含剛度矩陣和質(zhì)量矩陣的線性組合

    式中、分別表示質(zhì)量和剛度比例阻尼常數(shù)。

    式中分別為第、的固有頻率,分別為第、的阻尼比,根據(jù)文獻[30-31],假設2種模式具有相同的阻尼比,==,則阻尼比約為0.04。

    因此式(3)可簡化為

    3.2 蝸殼聲學模擬方法

    運用蝸殼聲學仿真用商用程序SYSNOISE中的邊界元法進行模擬。蝸殼結構振動引起的空氣輻射聲壓的控制方程如下

    提取結構外表面并將其網(wǎng)格化,作為邊界元計算中使用的聲學模型,如圖7所示。聲學網(wǎng)格共包括14 706個元素和13 573個節(jié)點。根據(jù)文獻[32],本研究中所研究模型的最大有效頻率為4 366 Hz,因此,對于葉片通過頻率而言,網(wǎng)格足夠精細。將結構外表面節(jié)點的法向速度轉(zhuǎn)移到聲學模型的表面節(jié)點,將其設置為蝸殼聲學計算的邊界條件,最終采用邊界元法求解聲壓分布。

    圖7 蝸殼聲網(wǎng)格

    4 結果與分析

    4.1 性能試驗結果

    圖8為4臺不同葉片數(shù)的泵在設計工況下的性能曲線。

    圖8 設計工況下不同葉片數(shù)的泵性能曲線

    由圖8可知,數(shù)值計算結果與試驗結果較為吻合。揚程、效率和軸功率誤差值分別在在5.06%、5.34%和5.68%以內(nèi)。隨著葉片數(shù)的增加,揚程總體上呈上升趨勢,軸功率逐漸增大,最大增幅分別為15.9%和14.1%。效率隨葉片數(shù)的增大呈先增大后減小再增大的趨勢,泵Ⅱ(葉片數(shù)為5)的效率達到最大值。這可能是因為當葉片數(shù)減少后,葉片對水流的約束減弱,泵內(nèi)出現(xiàn)流動分離,導致泵效率下降。當葉片數(shù)增加時,泵內(nèi)流體流態(tài)更加均勻,但摩擦損失也隨之增加,效率因此出現(xiàn)變化。就泵的能量性能而言,最佳葉片數(shù)為5。

    4.2 振動試驗

    為了驗證了數(shù)值模擬方法的正確性,對比了設計工況下A1測點振動加速度的試驗測試和數(shù)值預測的頻譜,并分析了口環(huán)泄漏對計算結果的影響,如圖9所示。

    圖9 監(jiān)測點振動加速度頻譜

    從圖9中可以看出,數(shù)值計算能夠很好地預測出特征頻率,振動頻譜中分別出現(xiàn)48和240 Hz 2個峰值。其中48 Hz時的振動由軸旋轉(zhuǎn)產(chǎn)生,是由機械不平衡或水力不平衡引起,240 Hz時的峰值是由轉(zhuǎn)子-定子相互作用激發(fā)的葉頻引起。

    進一步對比圖9中有無口環(huán)泄漏的結果,可以發(fā)現(xiàn)有口環(huán)方案預測得到的葉頻處對應的振幅較無口環(huán)方案更為精確,預測精度提高了13.5%。這表明了離心泵口環(huán)泄漏對水力振動性能的具有顯著影響。因此,考慮口環(huán)泄漏的離心泵振動模擬能夠提高數(shù)值計算的精確性。

    4.3 蝸殼壁面壓力脈動

    圖10為設計工況下泵在葉頻對應的壓力脈動幅值。

    圖10 設計工況蝸殼的壓力脈動幅值

    從圖10可以看出,蝸殼出口處的壓力脈動幅值均較大,壓力脈動幅值隨葉片數(shù)的減小而增大,最大增幅為23.6%。由于葉片與隔舌的相互作用,蝸殼隔舌、第1到第2斷面間和擴壓管壁面等3個區(qū)域的壓力脈動幅值相對較高,這可能是因為葉輪和隔舌間的動靜干涉作用的影響。在葉頻對應的壓力值下,葉片后緣剛好經(jīng)過隔舌前緣,導致了劇烈的壓力脈動。隨著葉片數(shù)減少,葉片間距逐漸增大,葉輪出口流動不均勻度增大,泵的振動和噪聲也隨之增大。

    4.4 蝸殼結構振動

    設計工況下4臺泵的振動位移如圖11所示。圖11表明,隨著葉片數(shù)的減少,振動位移有明顯的上升,泵I(葉片數(shù)為4)的變化最大,振動位移的最大增幅為36.4%,這與壓力脈動的變化相一致。因此,葉片數(shù)對振動位移有很大影響,最大位移主要發(fā)生第8斷面處。

    圖12是設計工況下的振動速度圖。由圖可知,振動速度隨著葉片數(shù)的增加先增大后減小,泵Ⅱ(葉片數(shù)為5)的振動速度幅值最大。根據(jù)圖12中的速度變化,振動的高速區(qū)主要集中在隔舌周圍、蝸殼第4到第6斷面間和擴壓管第8斷面處。但對比圖11可以發(fā)現(xiàn),速度的變化與位移不完全一致,這可能是因為振動速度不僅與位移有關,還與振動頻率有關。總體看葉片數(shù)為5時泵的振動最大。

    圖11 設計流量下的振動位移

    圖12 設計流量下的振動速度

    4.5 蝸殼輻射噪聲

    采用以蝸殼為中心的半徑為0.5 m的球形聲網(wǎng)格,計算了泵輻射聲壓級的方向性分布。圖13給出了4臺泵在設計工況下葉頻對應的聲壓級。從圖13可以看出,隨著葉片數(shù)的增加,泵在設計條件下的聲壓級(sound pressure level,SPL)先增大后減小,泵Ⅱ(葉片數(shù)為5)的SPL最大,與振動速度的變化一致,5葉片較4葉片增幅較大,最大增幅達138.5%。高聲壓級區(qū)域主要出現(xiàn)在泵出口附近的高振動速度引起的垂直方向上。

    圖14給出了設計流量下聲強方向性分布。由圖可以看出,聲強隨葉片數(shù)的增加先增大后減小,與聲壓級變化趨勢一致。5葉片較4葉片增幅較大,最大增幅為237.4%。但泵Ⅱ(葉片數(shù)為5)的聲強明顯高于其他3個方案,這表明了葉片數(shù)的確定應充分考慮振動噪聲的影響。

    圖13 設計流量下的聲壓級

    圖14 設計工況下聲強方向性分布

    綜合以上分析,從水力性能來看,泵I(葉片數(shù)為4)的性能最差,泵III(葉片數(shù)為6)的性能略差于泵Ⅱ(葉片數(shù)為5),但優(yōu)于泵Ⅳ(葉片數(shù)為7)??紤]振動噪聲的影響后,泵III的預測數(shù)據(jù)遠小于泵Ⅱ,接近泵Ⅳ。因此,綜合考慮到離心泵的外特性、壓力脈動和振動噪聲性能后,對于模型泵而言,最佳葉片數(shù)為6。

    5 結 論

    通過數(shù)值仿真和試驗研究的對比分析了葉片數(shù)對離心泵蝸殼的振動噪聲的影響,研究發(fā)現(xiàn):

    1)隨著葉片數(shù)的增加,揚程和軸功率逐漸增大,最大增幅分別為15.9%和14.1%;效率隨葉片數(shù)的增大呈先增大后減小再增大的趨勢。

    2)在數(shù)值模擬中考慮口環(huán)泄漏情況能夠提高離心泵振動仿真的精確度,有口環(huán)方案預測得到的葉頻對應的振幅較無口環(huán)方案更為精確,預測精度提高了13.5%。

    3)離心泵蝸殼的壓力脈動幅值和振動位移隨著葉片數(shù)的減小均有所增大,最大增幅分別為23.6%和36.4%;蝸殼隔舌、蝸殼第1到第2斷面間和擴壓管壁面這3個區(qū)域的壓力脈動幅值相對較大;蝸殼壁面上的最大位移主要發(fā)生第8斷面處。

    4)蝸殼表面振動速度隨著葉片數(shù)的增加先增大后減小,在隔舌周圍、蝸殼的第4與第6斷面之間和靠近擴壓管的第8斷面的振動速度較高。

    5)隨著葉片數(shù)的增加,泵在設計條件下的聲壓級和聲強先增大后減??;高聲壓級區(qū)域主要出現(xiàn)在泵出口附近的高振動速度引起的垂直方向上。

    6)考慮離心泵的外特性、壓力脈動和振動噪聲性能后,模型泵的最佳葉片數(shù)為6。

    [1]張克危. 流體機械原理(上冊)[M]. 北京:機械工業(yè)出版社,2000.

    [2]Birajdar R, Patil R, Khanzode K. Vibration and noise in centrifugal pumps-sources and diagnosis methods[C] Porto: Reliability and Failure, 2009.

    [3]Gulich J F. Centrifugal Pumps[M]. Berlin: Springer, 2010.

    [4]Baun D O, Flack R D. Effects of volute design and number of impeller blades on lateral impeller forces and hydraulic performance[J]. International Journal of Rotating Machinery, 2003, 9(2): 145-152.

    [5]何濤,鐘榮,孫玉東. 離心泵水動力噪聲計算方法研究[J].船舶力學,2012,16(4):449-455.

    He Tao, Zhong Rong, Sun Yudong. Numerical method on hydrodynamic noise of centrifugal pump[J]. CSSRC Reports, 2012, 16(4): 449-455. (in Chinese with English abstract)

    [6]Chakraborty S, Choudhuri K, Dutta P, et al. Performance prediction of Centrifugal Pumps with variations of blade number[J]. Journal of Scientific & Industrial Research, 2013, 72(6): 373-378.

    [7]馮濤,王晶,吳瑞. 蝸舌間隙對離心泵流動噪聲影響的研究[J]. 食品與機械,2012,28(2):75-78.

    Feng Tao, Wang Jing, Wu Rui. Research on the influence to the flow noise of the centrifugal pumps with different gap between impeller and tongue[J]. Food & Machinery, 2012, 28(2): 75-78. (in Chinese with English abstract)

    [8]譚永學,王宏光,楊愛玲,等. 離心泵水動力噪聲預測[J].上海理工大學學報,2011,33(1):89-94.

    Tan Yongxue, Wang Hongguang, Yang Ailing, et al. Numerical prediction of hydrodynamic noise for a centrifugal pumps[J]. Journal of University of Shanghai for Science and Technology, 2011, 33(1): 89-94. (in Chinese with English abstract)

    [9]Simpson H C, Clark T A, Weir G A. A theoretical investigation of hydraulic noise in pumps[J]. Journal of Sound and Vibration, 1967, 5(3): 456-488.

    [10]Chu S, Dong R, Katz J. Relationship between unsteady flow, pressure fluctuations, and noise in a centrifugal pump-part B: Effects of blade-tongue interactions[J]. Journal of Fluids Engineering, 1995, 117(1): 30-35.

    [11]Morgenroth M, Weaver D S. Sound generation by a centrifugal pump at blade passage frequency[J]. Journal of Turbomachinery, 1998, 120(4): 736-743.

    [12]Rzentkowski G, Zbroja S. Experimental characterization of centrifugal pumps as an acoustic source at the blade-passing frequency[J]. Journal of Fluids and Structures, 2000, 14(4): 529-558.

    [13]Choi J S, McLaughlin D K, Thompson D E. Experiments on the unsteady flow field and noise generation in a centrifugal pump impeller[J]. Journal of Sound and Vibration, 2003, 263(3): 493-514.

    [14]董亮,代翠,孔繁余,等. 葉片出口安放角對離心泵作透平噪聲的影響[J]. 農(nóng)業(yè)工程學報,2015(6):77-83.

    Dong Liang, Dai Cui, Kong Fanyu, et al. Effect of blade outlet angle on turbine noise of centrifugal pump[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of CASE), 2015(6): 77-83. (in Chinese with English abstract)

    [15]Kato C, Kaiho M, Manabe A. An overset finite-element large-eddy simulation method with applications to turbomachinery and aeroacoustics[J]. Journal of Applied Mechanics, 2003, 70(1): 32-43.

    [16]Jeon W H, Lee D J. A numerical study on the flow and sound fields of centrifugal impeller located near a wedge[J]. Journal of Sound and Vibration, 2003, 266(4): 785-804.

    [17]Langthjem M A, Olhoff N. A numerical study of flow-induced noise in a two dimensional centrifugal pump. Part I. Hydrodynamics[J]. Journal of Fluids and Structures, 2004, 19: 349-368.

    [18]Langthjem M A, Olhoff N. A numerical study of flow-induced noise in a two-dimensional centrifugal pump. Part II. Hydroacoustics[J]. Journal of Fluids and Structures, 2004, 19: 369-386.

    [19]李躍,施衛(wèi)東,韓笑笑. 不同結構形式對串并聯(lián)離心泵振動特性的影響[J]. 排灌機械工程學報,2015,33(9):744-749.

    Li Yue, Shi Weidong, Han Xiaoxiao. Effects of pump hydraulic structure on vibration characteristics of series-parallel centrifugal pump[J]. Journal of Drainage and Irrigation Machinery Engineering, 2015, 33(9): 744-749. (in Chinese with English abstract)

    [20]Kato C, Yamade Y, Wang H, et al. Numerical prediction of sound generated from flows with a low Mach number[J]. Computers & Fluids, 2007, 36(1): 53-68.

    [21]Jiang Y Y, Yoshimura S, Imai R, et al. Quantitative evaluation of flow-induced structural vibration and noise in turbomachinery by full-scale weakly coupled simulation[J]. Journal of Fluids and Structures, 2007, 23(4): 531-544.

    [22]Adkins D R, Brennen C E. Analysis of hydrodynamic radial forces on centrifugal pump impellers[J]. Journal of Fluids Engineering, 1988, 110(1): 20-28.

    [23]Guinzburg A, Brennen C E, Acosta A J, et al. Experimental results for the rotor dynamic characteristics of leakage flows in centrifugal pumps[J]. Journal of Fluid Engineering, 1994, 116(1): 110-115.

    [24]王勇,劉厚林,袁壽其,等. 葉片數(shù)對離心泵空化誘導振動噪聲的影響[J]. 哈爾濱工程大學學報,2012,33(11):1405-1409.

    Wang Yong, Liu Houlin, Yuan Shouqi, et al. Effects of the blade number on cavitation-induced vibration and noise of centrifugal pumps[J]. Journal of Harbin Engineering University, 2012, 33(11): 1405-1409. (in Chinese with English abstract)

    [25]Tucker P G, Lardeau S. Introduction: Applied large eddy simulation[J]. Phil Trans Roy Soc Lond A-Math Phys Eng Sci, 2009, 367(1899): 2809-2818.

    [26]Bonet J, Peraire J. An alternating digital tree (ADT) algorithm for 3D geometric searching and intersection problems[J]. International Journal for Numerical Methods in Engineering, 1991, 31(1): 1-17.

    [27]Samareh J A. Discrete data transfer technique for fluid-structure interaction[C]//Hampton: NASA Langley Research Center, 2007: 1-12.

    [28]邵春雷,顧伯勤,陳曄. 離心泵內(nèi)部非定常壓力場的數(shù)值研究[J]. 農(nóng)業(yè)工程學報,2009,25(1):87-92.

    Shao Chunlei, Gu Boqin, Chen Ye. Numerical simulation of unsteady pressure field in centrifugal pumps[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of CASE), 2009, 25(1): 87-92

    [29]鄭水華,錢亨,牟介剛,等. 交錯葉片對三通道蝸殼離心泵水動力性能的影響[J]. 農(nóng)業(yè)工程學報,2015,31(23):51-59.

    Zheng Shuihua, Qian Heng, Mou Jiegang, et al. Effect of staggered blades on hydrodynamic performance of three channel volute centrifugal pump[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of CASE), 2015, 31(23): 51-59. (in Chinese with English abstract)

    [30]李德順,王成澤,李銀然,等. 葉片前緣磨損形貌特征對風力機翼型氣動性能的影響[J]. 農(nóng)業(yè)工程學報,2017,33(22):277-283.

    Li Deshun, Wang Chengze, Li Yinran, et al. Effect of blade leading edge wear morphology on aerodynamic performance of wind wing type[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of CASE), 2017, 33(22): 277-283. (in Chinese with English abstract)

    [31]Kaiser T F, Osman R H, Dickau R O. Analysis guide for variable frequency drive operated centrifugal pumps[C]// Texas: Texas A & M University, 2008: 81-106.

    [32]Marburg S. Six boundary elements per wave length: Is that enough[J]. Journal of Computational Acoustics, 2002, 10(1): 25-51.

    Effects of blade number on flow induced vibration and noise in centrifugal pump

    Tan Minggao1, Lu Youdong1, Wu Zejin1, Wu Xianfang2, Liu Houlin1

    (1.,212013,; 2.,,212013,)

    The number of blades is one of the main geometric parameters of centrifugal pump, which is widely used in agricultural machinery, and it has an important influence on the vibration and noise of centrifugal pumps. Both vibration and noise can affect the centrifugal pump performance and its life, and the sources of vibration and noise may lie in hydraulic or mechanical aspects. In fact, most previous works for vibration and noise of centrifugal pumps mostly focused on theoretical and experimental studies. However, these studies seem to simulate the volute and impeller interactions only, without consideration of the leakage flow paths. The leakage flow paths between the rotating impeller and the stationary housing play an important role in centrifugal pumps. Therefore, understanding the influence of the blade number and the leakage flow paths in centrifugal pump is an urgent problem to be solved. In this paper, the vibration and sound radiation of volute under flow excitation was simulated by FEM/BEM acoustic-vibration method. The experiment was carried out to study the effects of blade number on the vibration and noise based on a centrifugal pump with a single entry and a single volute.Comparing the different performances of centrifugal pump with various blade numbers, it was found that with the increase of blade number, the head and shaft power increased gradually, and the efficiency increased first, then decreased and increased with the increase of blade number. In addition, the numerical simulation results of volute with and without leakage flow paths were compared. The vibration and noise induced by inner flow of the pump with different blade number were analyzed under design flow condition. It was found that the results of simulation were validated by the vibration acceleration of the monitoring points on volute compared with the experimental vibration acceleration. The numerical simulation method proposed in this paper could be used to predict flow-induced vibration and acoustic radiation of volutes under design conditions. The errors value of head, efficiency and shaft power between numerical calculation and experiments were within 5.06%, 5.34% and 5.68% respectively. The amplitudes of simulation with the leakage flow paths were coincident with the experimental results than the results without the leakage. The peak error between the simulation amplitude with and without leakage flow paths was 13.5%. To reveal the effects of the blade number on pressure fluctuation and vibration, the contrast with different blade numbers was considered objectively. As the number of blades decreased, the pressure fluctuation and vibration displacement of the volute of centrifugal pump increased. High amplitude regions appeared at the volute tongue, the first and second hydraulic profile of the volute and the eighth hydraulic profile were close to the diffuser. The maximum vibration displacement mainly concentrated at the eighth hydraulic profile. According to the analysis, the significant high levels of vibration velocity can mainly classify in three regions, around the tongue, between the fourth and the sixth hydraulic profile of the volute, and the eighth hydraulic profile was close to the diffuser. However, the variation of the velocity disagreed with the displacement. That meant the vibration velocity was not only relative to the displacement, but also relative to the frequency. In terms of the noise on these impeller with diverse blade number, the vibration speed and radiated sound pressure level of the volute surface first increased and then decreased with the increase of the number of blades, besides, when the impeller was five blades, the vibration speed and noise reached the maximum value. The region of high noise level mainly appears in the vertical direction. The results can provide a reference for the further analysis on vibration and noise reduction design of centrifugal pump.

    centrifugal pump; numerical simulation; pressure; leakage flow paths; number of blades; vibration; noise

    談明高,陸友東,吳澤瑾,吳賢芳,劉厚林. 葉片數(shù)對離心泵振動噪聲性能的影響[J]. 農(nóng)業(yè)工程學報,2019,35(23):73-79.doi:10.11975/j.issn.1002-6819.2019.23.009 http://www.tcsae.org

    Tan Minggao, Lu Youdong, Wu Zejin, Wu Xianfang, Liu Houlin. Effects of blade number on flow induced vibration and noise in centrifugal pump[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(23): 73-79. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2019.23.009 http://www.tcsae.org

    2019-07-20

    2019-10-30

    國家自然科學基金(51679110、51779108);江蘇省自然科學基金B(yǎng)K20161350;江蘇省農(nóng)業(yè)重點研發(fā)計劃(BE2017356)

    談明高,研究員,主要研究高效高可靠性葉片泵水力模型。Email:tmgwxf@ujs.edu.cn

    10.11975/j.issn.1002-6819.2019.23.009

    TH312

    A

    1002-6819(2019)-23-0073-07

    猜你喜歡
    蝸殼揚程離心泵
    焊接蝸殼泵模型開發(fā)
    水泵技術(2022年2期)2022-06-16 07:08:04
    一種改善離心泵運行狀態(tài)的方法
    水泵技術(2021年5期)2021-12-31 05:26:48
    大型立式單級引黃離心泵軸向力平衡的研究
    水泵技術(2021年5期)2021-12-31 05:26:40
    離心泵流場外特性及激勵數(shù)值計算
    防爆電機(2021年5期)2021-11-04 08:16:34
    管路受力誘發(fā)高揚程離心泵振動加劇原因分析
    水泵技術(2021年3期)2021-08-14 02:09:22
    高比速離心泵的高效設計研究
    水泵技術(2021年3期)2021-08-14 02:09:18
    污水泵蝸殼斷裂原因分析
    水泵技術(2021年6期)2021-02-16 01:14:50
    新型中高揚程大流量水錘泵結構技術改進研究
    中國水利(2015年7期)2015-02-28 15:12:58
    蝸殼差壓法在印尼ASAHAN一級水電站的應用
    高揚程長距離輸水管道系統(tǒng)水錘防護的模擬分析
    曰老女人黄片| 国产成人a∨麻豆精品| 精品久久蜜臀av无| 极品人妻少妇av视频| av福利片在线| 国产黄频视频在线观看| 国产精品久久久久久人妻精品电影 | 久久精品亚洲av国产电影网| 国产黄频视频在线观看| 一区在线观看完整版| 精品国产一区二区久久| 亚洲人成电影免费在线| 99精品久久久久人妻精品| 国产熟女午夜一区二区三区| 国产欧美日韩综合在线一区二区| 69精品国产乱码久久久| 女警被强在线播放| 亚洲免费av在线视频| 欧美黄色淫秽网站| 老司机福利观看| 欧美中文综合在线视频| 欧美中文综合在线视频| av一本久久久久| 黄色怎么调成土黄色| 人人妻人人添人人爽欧美一区卜| 我的亚洲天堂| 国产免费一区二区三区四区乱码| 亚洲中文日韩欧美视频| 精品国产国语对白av| 无限看片的www在线观看| 欧美乱码精品一区二区三区| 欧美性长视频在线观看| 少妇 在线观看| 亚洲欧美精品综合一区二区三区| 亚洲精品国产av蜜桃| 亚洲专区中文字幕在线| 中文字幕制服av| 免费黄频网站在线观看国产| a级片在线免费高清观看视频| 日本精品一区二区三区蜜桃| 欧美日韩亚洲综合一区二区三区_| 亚洲欧美精品自产自拍| 精品免费久久久久久久清纯 | 国产精品一区二区精品视频观看| 久热爱精品视频在线9| 久久久精品区二区三区| 久久人人97超碰香蕉20202| 日韩制服丝袜自拍偷拍| 亚洲全国av大片| 亚洲 国产 在线| 免费黄频网站在线观看国产| 考比视频在线观看| 精品人妻在线不人妻| 国产精品久久久人人做人人爽| 99精品欧美一区二区三区四区| 十八禁高潮呻吟视频| 免费黄频网站在线观看国产| av在线老鸭窝| 日韩视频在线欧美| 国产真人三级小视频在线观看| 精品一品国产午夜福利视频| 国产欧美日韩一区二区三 | 精品福利永久在线观看| 久久免费观看电影| 欧美日韩成人在线一区二区| 亚洲色图 男人天堂 中文字幕| 如日韩欧美国产精品一区二区三区| 麻豆国产av国片精品| 精品熟女少妇八av免费久了| 一区二区av电影网| 欧美黄色淫秽网站| 中文字幕色久视频| 久久午夜综合久久蜜桃| 色视频在线一区二区三区| 久久久久久久久免费视频了| 精品少妇内射三级| 亚洲va日本ⅴa欧美va伊人久久 | 波多野结衣av一区二区av| 80岁老熟妇乱子伦牲交| 美女高潮喷水抽搐中文字幕| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲久久久国产精品| 一区二区日韩欧美中文字幕| 一区二区日韩欧美中文字幕| 黄色视频不卡| 欧美日韩亚洲综合一区二区三区_| 国产精品成人在线| 成人国产av品久久久| 美女扒开内裤让男人捅视频| 天堂8中文在线网| 在线观看舔阴道视频| 日韩 亚洲 欧美在线| 嫩草影视91久久| 美女午夜性视频免费| 国产不卡av网站在线观看| avwww免费| 91麻豆精品激情在线观看国产 | 天天躁狠狠躁夜夜躁狠狠躁| 手机成人av网站| 久久久精品区二区三区| 天天躁夜夜躁狠狠躁躁| 男女午夜视频在线观看| 一区福利在线观看| 国产伦理片在线播放av一区| 亚洲av日韩精品久久久久久密| 久久中文字幕一级| 中亚洲国语对白在线视频| 久久九九热精品免费| 免费观看a级毛片全部| 亚洲欧美成人综合另类久久久| 十分钟在线观看高清视频www| 国产成人啪精品午夜网站| 久久久久久久久免费视频了| 91精品伊人久久大香线蕉| 老司机亚洲免费影院| 一级黄色大片毛片| 国产深夜福利视频在线观看| 日本欧美视频一区| 亚洲视频免费观看视频| 亚洲精品国产精品久久久不卡| www.av在线官网国产| 91av网站免费观看| 又大又爽又粗| 中文字幕人妻熟女乱码| 18禁裸乳无遮挡动漫免费视频| 男女国产视频网站| 在线观看免费高清a一片| av免费在线观看网站| 91字幕亚洲| 欧美97在线视频| 国产免费av片在线观看野外av| 大陆偷拍与自拍| 满18在线观看网站| 中文欧美无线码| 少妇的丰满在线观看| 黑人欧美特级aaaaaa片| 91精品国产国语对白视频| 男人爽女人下面视频在线观看| 欧美亚洲日本最大视频资源| 精品亚洲乱码少妇综合久久| 国产不卡av网站在线观看| 人人妻人人澡人人爽人人夜夜| 亚洲一码二码三码区别大吗| 亚洲色图 男人天堂 中文字幕| 操美女的视频在线观看| 国产高清国产精品国产三级| av国产精品久久久久影院| 我要看黄色一级片免费的| 亚洲精品成人av观看孕妇| 国产亚洲av片在线观看秒播厂| 国产欧美日韩一区二区三 | 99香蕉大伊视频| 亚洲国产中文字幕在线视频| 亚洲 欧美一区二区三区| 91精品伊人久久大香线蕉| 美女国产高潮福利片在线看| 性色av一级| 极品少妇高潮喷水抽搐| 国产91精品成人一区二区三区 | 人人妻人人澡人人爽人人夜夜| 亚洲七黄色美女视频| 成人手机av| 久久这里只有精品19| www.熟女人妻精品国产| 国产成人a∨麻豆精品| 韩国高清视频一区二区三区| 十八禁高潮呻吟视频| 热99久久久久精品小说推荐| 国产av精品麻豆| 日本av免费视频播放| 久久久久精品国产欧美久久久 | 999久久久国产精品视频| 一个人免费看片子| 国产91精品成人一区二区三区 | 最近中文字幕2019免费版| 最新的欧美精品一区二区| 亚洲国产av影院在线观看| www.精华液| 精品国产国语对白av| 欧美亚洲日本最大视频资源| av欧美777| 天天添夜夜摸| 操出白浆在线播放| 欧美另类一区| 国产在视频线精品| 久久综合国产亚洲精品| 国产老妇伦熟女老妇高清| 嫩草影视91久久| 国产一区二区激情短视频 | 下体分泌物呈黄色| 丝袜喷水一区| 亚洲精品美女久久av网站| 在线看a的网站| 操美女的视频在线观看| 50天的宝宝边吃奶边哭怎么回事| 久久精品国产综合久久久| 97精品久久久久久久久久精品| 成年av动漫网址| 国产免费av片在线观看野外av| 国产在线一区二区三区精| 国产91精品成人一区二区三区 | 最近中文字幕2019免费版| 青草久久国产| 高清av免费在线| 一级毛片精品| 国产亚洲一区二区精品| 搡老岳熟女国产| 捣出白浆h1v1| 欧美精品人与动牲交sv欧美| 国产精品秋霞免费鲁丝片| 黑人欧美特级aaaaaa片| 精品福利观看| 免费黄频网站在线观看国产| 日韩三级视频一区二区三区| 日韩大码丰满熟妇| 成人国产av品久久久| 亚洲精品乱久久久久久| 在线观看免费高清a一片| 国产日韩欧美视频二区| 国产又爽黄色视频| 美女脱内裤让男人舔精品视频| 亚洲熟女精品中文字幕| 天堂8中文在线网| 国产精品1区2区在线观看. | 中国美女看黄片| 精品久久蜜臀av无| 亚洲自偷自拍图片 自拍| 久久人人97超碰香蕉20202| 日韩熟女老妇一区二区性免费视频| 亚洲人成电影观看| 午夜成年电影在线免费观看| 国产精品久久久久成人av| 欧美 日韩 精品 国产| 午夜福利,免费看| 这个男人来自地球电影免费观看| 首页视频小说图片口味搜索| 91九色精品人成在线观看| 国产在线一区二区三区精| 国产在线免费精品| 国产色视频综合| 国产一区二区 视频在线| 黄网站色视频无遮挡免费观看| 久久av网站| 欧美精品一区二区免费开放| 我的亚洲天堂| 蜜桃国产av成人99| 捣出白浆h1v1| 两性夫妻黄色片| 久久综合国产亚洲精品| 999精品在线视频| 精品人妻熟女毛片av久久网站| 免费在线观看黄色视频的| 久久99热这里只频精品6学生| 精品国产超薄肉色丝袜足j| 午夜福利一区二区在线看| 啪啪无遮挡十八禁网站| 多毛熟女@视频| 黄色视频在线播放观看不卡| 日韩欧美一区二区三区在线观看 | 午夜激情久久久久久久| 欧美日韩国产mv在线观看视频| 男女午夜视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 精品欧美一区二区三区在线| 中文字幕高清在线视频| 日韩 亚洲 欧美在线| 精品人妻1区二区| 丝袜美腿诱惑在线| 少妇的丰满在线观看| www日本在线高清视频| 男人添女人高潮全过程视频| 男女下面插进去视频免费观看| av又黄又爽大尺度在线免费看| 又黄又粗又硬又大视频| 涩涩av久久男人的天堂| 精品国产一区二区三区久久久樱花| 日本撒尿小便嘘嘘汇集6| 天天躁日日躁夜夜躁夜夜| 国产野战对白在线观看| 女性生殖器流出的白浆| 老司机午夜十八禁免费视频| 天天躁日日躁夜夜躁夜夜| 亚洲精品一区蜜桃| 91成年电影在线观看| 国产精品亚洲av一区麻豆| 我要看黄色一级片免费的| 天天躁狠狠躁夜夜躁狠狠躁| 51午夜福利影视在线观看| 国产成人免费观看mmmm| 一区二区三区乱码不卡18| 免费黄频网站在线观看国产| 国产伦人伦偷精品视频| 女人精品久久久久毛片| 夜夜骑夜夜射夜夜干| 日本91视频免费播放| 亚洲国产欧美一区二区综合| 精品一区二区三区av网在线观看 | 精品少妇久久久久久888优播| 亚洲欧美激情在线| 99re6热这里在线精品视频| 久久香蕉激情| 日日摸夜夜添夜夜添小说| 久久久国产欧美日韩av| www.熟女人妻精品国产| 欧美 亚洲 国产 日韩一| 国产在线观看jvid| 亚洲精品中文字幕一二三四区 | 国产av国产精品国产| 女人高潮潮喷娇喘18禁视频| 亚洲一码二码三码区别大吗| 亚洲欧洲精品一区二区精品久久久| 欧美亚洲日本最大视频资源| 亚洲av电影在线进入| 免费不卡黄色视频| 国产亚洲精品一区二区www | 高清在线国产一区| 婷婷色av中文字幕| 天天添夜夜摸| 婷婷丁香在线五月| 啦啦啦中文免费视频观看日本| 国产一级毛片在线| 99国产精品免费福利视频| 少妇人妻久久综合中文| 啦啦啦在线免费观看视频4| 老司机亚洲免费影院| 久久久精品免费免费高清| 国产欧美亚洲国产| 日本av免费视频播放| 亚洲成人手机| 亚洲成人免费av在线播放| 免费观看人在逋| 亚洲国产中文字幕在线视频| 又紧又爽又黄一区二区| 精品乱码久久久久久99久播| 久久久国产成人免费| 精品卡一卡二卡四卡免费| 天天影视国产精品| av网站免费在线观看视频| 黄色a级毛片大全视频| 丝瓜视频免费看黄片| 久久中文字幕一级| 正在播放国产对白刺激| 国产无遮挡羞羞视频在线观看| 视频区欧美日本亚洲| av又黄又爽大尺度在线免费看| 日韩中文字幕欧美一区二区| 日韩,欧美,国产一区二区三区| 一区二区av电影网| 老司机深夜福利视频在线观看 | 欧美精品av麻豆av| 99国产极品粉嫩在线观看| 大型av网站在线播放| 久久99热这里只频精品6学生| 亚洲国产毛片av蜜桃av| 大型av网站在线播放| 我的亚洲天堂| 日韩制服丝袜自拍偷拍| 叶爱在线成人免费视频播放| 大片电影免费在线观看免费| 日韩欧美一区二区三区在线观看 | 日日摸夜夜添夜夜添小说| 亚洲欧美色中文字幕在线| 老鸭窝网址在线观看| 高潮久久久久久久久久久不卡| 天天躁狠狠躁夜夜躁狠狠躁| 免费看十八禁软件| 真人做人爱边吃奶动态| 国产成人系列免费观看| 欧美精品人与动牲交sv欧美| 青春草亚洲视频在线观看| 中文字幕制服av| 国产成人精品无人区| 国产日韩一区二区三区精品不卡| 国产精品影院久久| 大型av网站在线播放| 两个人看的免费小视频| 成人国产一区最新在线观看| 天天操日日干夜夜撸| 久久人人爽av亚洲精品天堂| www.自偷自拍.com| 狠狠狠狠99中文字幕| 人人妻人人添人人爽欧美一区卜| videos熟女内射| 最新的欧美精品一区二区| 亚洲欧美一区二区三区久久| 午夜福利乱码中文字幕| 国产极品粉嫩免费观看在线| a 毛片基地| 国产在线免费精品| 高清av免费在线| 欧美日韩亚洲高清精品| 99国产综合亚洲精品| 高清av免费在线| 热re99久久精品国产66热6| 少妇裸体淫交视频免费看高清 | 免费在线观看影片大全网站| 久久精品国产亚洲av高清一级| 青春草亚洲视频在线观看| 精品亚洲乱码少妇综合久久| 欧美精品一区二区免费开放| 两性午夜刺激爽爽歪歪视频在线观看 | 中文精品一卡2卡3卡4更新| 少妇 在线观看| 99精国产麻豆久久婷婷| 国产日韩一区二区三区精品不卡| 亚洲第一av免费看| 三级毛片av免费| 建设人人有责人人尽责人人享有的| 黄色视频,在线免费观看| 精品少妇久久久久久888优播| 国产精品香港三级国产av潘金莲| 人人妻人人澡人人看| 97人妻天天添夜夜摸| 女性生殖器流出的白浆| 欧美大码av| 人妻 亚洲 视频| 日本vs欧美在线观看视频| 麻豆乱淫一区二区| 黄频高清免费视频| 中国国产av一级| 成年动漫av网址| 桃花免费在线播放| 不卡一级毛片| 亚洲精品自拍成人| 我要看黄色一级片免费的| 黄色视频不卡| www.熟女人妻精品国产| 亚洲三区欧美一区| 国产精品免费视频内射| 久久久水蜜桃国产精品网| 色综合欧美亚洲国产小说| 丰满饥渴人妻一区二区三| 欧美午夜高清在线| 国产麻豆69| 亚洲精品成人av观看孕妇| 俄罗斯特黄特色一大片| 精品第一国产精品| 国产真人三级小视频在线观看| 亚洲色图 男人天堂 中文字幕| 欧美 亚洲 国产 日韩一| 在线精品无人区一区二区三| 人妻久久中文字幕网| 亚洲精品国产av蜜桃| 最近最新免费中文字幕在线| 精品久久久精品久久久| 人妻一区二区av| 法律面前人人平等表现在哪些方面 | 精品人妻在线不人妻| 国产一区二区在线观看av| 久久久精品区二区三区| 久久人妻熟女aⅴ| 久久99热这里只频精品6学生| 欧美黑人精品巨大| 成人国产av品久久久| 精品少妇一区二区三区视频日本电影| 精品国产乱码久久久久久男人| 久久热在线av| 久热这里只有精品99| 亚洲国产欧美在线一区| 亚洲少妇的诱惑av| 亚洲天堂av无毛| 九色亚洲精品在线播放| 91av网站免费观看| 99国产精品99久久久久| 99国产精品免费福利视频| 高清在线国产一区| 欧美av亚洲av综合av国产av| 视频区图区小说| 亚洲 国产 在线| 成人黄色视频免费在线看| 久久人人爽人人片av| 欧美在线一区亚洲| 我要看黄色一级片免费的| 叶爱在线成人免费视频播放| 精品久久久精品久久久| 国产av精品麻豆| 美女福利国产在线| 99国产精品免费福利视频| 黄色a级毛片大全视频| 1024香蕉在线观看| 欧美日韩亚洲国产一区二区在线观看 | 另类亚洲欧美激情| 国产成人啪精品午夜网站| 欧美精品高潮呻吟av久久| e午夜精品久久久久久久| 伊人久久大香线蕉亚洲五| 成人黄色视频免费在线看| 在线观看www视频免费| 人成视频在线观看免费观看| 两性午夜刺激爽爽歪歪视频在线观看 | 老司机在亚洲福利影院| 男女之事视频高清在线观看| 一本综合久久免费| 91成年电影在线观看| 日日摸夜夜添夜夜添小说| 99re6热这里在线精品视频| 国产片内射在线| 精品一品国产午夜福利视频| 18禁裸乳无遮挡动漫免费视频| 亚洲伊人色综图| 精品一区二区三区四区五区乱码| 每晚都被弄得嗷嗷叫到高潮| 国产精品成人在线| 亚洲成人免费电影在线观看| 亚洲激情五月婷婷啪啪| 狂野欧美激情性bbbbbb| 在线看a的网站| 免费日韩欧美在线观看| 免费女性裸体啪啪无遮挡网站| 99国产精品99久久久久| 欧美日韩中文字幕国产精品一区二区三区 | 日韩欧美国产一区二区入口| 久9热在线精品视频| 午夜精品国产一区二区电影| 99香蕉大伊视频| 国产亚洲一区二区精品| 午夜福利在线观看吧| 久热爱精品视频在线9| 18在线观看网站| 久久人人爽人人片av| 国产福利在线免费观看视频| 精品少妇久久久久久888优播| 黄片小视频在线播放| 成人18禁高潮啪啪吃奶动态图| 国产免费av片在线观看野外av| 午夜成年电影在线免费观看| avwww免费| 久久久国产一区二区| 国产又爽黄色视频| 99热网站在线观看| 精品亚洲乱码少妇综合久久| 免费不卡黄色视频| 精品亚洲乱码少妇综合久久| av片东京热男人的天堂| 精品人妻1区二区| 亚洲av男天堂| 久久久国产欧美日韩av| 各种免费的搞黄视频| 国产又爽黄色视频| 成人三级做爰电影| 精品一区二区三卡| av网站在线播放免费| 日韩电影二区| 熟女少妇亚洲综合色aaa.| 精品一品国产午夜福利视频| 国产男女内射视频| 久久久精品免费免费高清| 亚洲欧洲日产国产| 午夜福利一区二区在线看| 欧美另类一区| 久久毛片免费看一区二区三区| www.精华液| 久久久欧美国产精品| 国产一区二区 视频在线| 亚洲男人天堂网一区| 天堂8中文在线网| 亚洲国产日韩一区二区| 搡老乐熟女国产| 精品一品国产午夜福利视频| 国产一区二区在线观看av| 黑人巨大精品欧美一区二区mp4| 黄色视频不卡| 国产亚洲欧美精品永久| 精品少妇一区二区三区视频日本电影| 黄片小视频在线播放| a在线观看视频网站| 久久99热这里只频精品6学生| av在线app专区| 1024视频免费在线观看| 色老头精品视频在线观看| 99精国产麻豆久久婷婷| 精品久久久精品久久久| 亚洲欧美一区二区三区黑人| 老司机亚洲免费影院| 满18在线观看网站| 国产精品1区2区在线观看. | 国产精品99久久99久久久不卡| 亚洲一区二区三区欧美精品| 亚洲欧美日韩另类电影网站| 中国美女看黄片| 黑人巨大精品欧美一区二区mp4| 波多野结衣一区麻豆| 999久久久国产精品视频| 国产麻豆69| 欧美在线黄色| 国产激情久久老熟女| 考比视频在线观看| 久久综合国产亚洲精品| 色94色欧美一区二区| 老司机深夜福利视频在线观看 | 色老头精品视频在线观看| 丝袜人妻中文字幕| 精品国产超薄肉色丝袜足j| 两性夫妻黄色片| 久久人人97超碰香蕉20202| 久久久久久亚洲精品国产蜜桃av| 久久精品熟女亚洲av麻豆精品| 亚洲精品国产色婷婷电影| 日本精品一区二区三区蜜桃| 99国产精品99久久久久| 天堂中文最新版在线下载| 亚洲欧美精品综合一区二区三区| 青草久久国产| 成年av动漫网址| 亚洲精品乱久久久久久| 久久精品国产亚洲av香蕉五月 | 国产男女超爽视频在线观看| 不卡av一区二区三区| 色精品久久人妻99蜜桃| 超碰97精品在线观看| 久久久精品国产亚洲av高清涩受| 午夜影院在线不卡| 亚洲精品国产色婷婷电影| 99国产精品免费福利视频| 美女福利国产在线|