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

    小流量工況下旋渦自吸泵流動降噪優(yōu)化研究

    2017-12-04 09:21:50胡日新湯海濤余昊謙
    農業(yè)機械學報 2017年11期
    關鍵詞:節(jié)距泵體脈動

    王 洋 胡日新 湯海濤 汪 群 余昊謙 賈 靜

    (江蘇大學流體機械工程技術研究中心, 鎮(zhèn)江 212013)

    小流量工況下旋渦自吸泵流動降噪優(yōu)化研究

    王 洋 胡日新 湯海濤 汪 群 余昊謙 賈 靜

    (江蘇大學流體機械工程技術研究中心, 鎮(zhèn)江 212013)

    為了研究旋渦自吸泵內部聲場特性,減少旋渦自吸泵在小流量工況下運行時的噪聲,采用CFD+Lighthill聲類比理論對旋渦泵內部聲場進行求解,并進行降噪優(yōu)化。首先利用CFX軟件提供的RNGk-ε模型,選取0.4Qd、0.8Qd及Qd3個流量點對旋渦泵進行非定常數(shù)值模擬,獲得不同工況下旋渦自吸泵內部壓力脈動情況。提取非定常計算所得的脈動力,導入聲學軟件ACTRAN中進行聲場計算,得到不同工況下旋渦泵出口聲壓級的大小、泵體內部的聲壓分布以及泵內部主要噪聲源分布。采用正弦調制葉片分布方式設計了調制角A為2°、4°、6°、8°的4種不同不等節(jié)距葉輪,通過對模型泵進行聲場計算,最后選取降噪效果較好的調制角A=4°的不等節(jié)距葉輪進行3D打印,并對3D打印樣品進行外特性和噪聲試驗驗證。結果表明:旋渦泵內部流動誘導噪聲與壓力脈動密切相關,主要是由葉輪與泵體的動靜干涉引起,其頻率特性與壓力脈動相似,2倍葉頻時聲壓最大。流量越小,旋渦泵出口聲壓級越大,其噪聲源主要分布在泵體流道及靠近流道出口隔舌處。調制角A=4°的不等節(jié)距葉片在小流量工況下對旋渦泵外特性性能影響不大,且能使小流量工況0.4Qd的噪聲下降2 dB,設計工況噪聲下降4 dB。

    旋渦自吸泵; 壓力脈動; 流動誘導噪聲; 不等節(jié)距葉片; 數(shù)值模擬

    引言

    旋渦自吸泵是一種小流量、高揚程泵,運行時效率較低,噪聲比較大,具有自吸能力,廣泛應用于輸送小流量的酸、堿等腐蝕性液體,汽油、酒精等易揮發(fā)性液體以及船舶供水、鍋爐給水等方面[1]。

    國內外學者對旋渦泵的研究主要關注旋渦泵各參數(shù)對旋渦泵性能的影響,通過研究改變葉片形狀、分布方式以及旋渦泵內部流道等方法來提高旋渦泵揚程和效率、降低能耗等[2-4],對旋渦泵噪聲的研究相對較少。而隨著綠色工業(yè)的推進,旋渦泵運行時噪聲過大的問題開始逐漸制約著其市場的進一步發(fā)展,越來越引起人們重視。因此,從內流機理入手,對旋渦泵聲學特性進行預測,并配合相關的抑制和優(yōu)化手段對旋渦泵進行降噪具有實際意義[5-6]。

    隨著現(xiàn)代聲學仿真技術的發(fā)展,采用數(shù)值模擬技術對泵內部流動誘導噪聲進行預測已成為降噪優(yōu)化的重要手段[7]。目前使用最為廣泛的是混合法[8],它把聲學數(shù)值模擬分為兩步:第一步,利用CFD求解流場。將流體看成不可壓縮流體,采用Reynolds時均法,選擇合適的湍流模型對時均化的N-S方程進行封閉,求解控制方程。目前應用最廣泛的湍流模型為兩方程湍流模型[9-10]。第二步,基于CFD得到的流動信息,根據(jù)Lighthill 聲類比理論定義等價聲源,進行聲學計算[11-12]。

    本文針對浙江省某企業(yè)一款旋渦自吸泵小流量工況運行時噪聲過大的問題,基于Lighthill聲類比理論采用間接混合法著重研究0.4Qd、0.8Qd和Qd3個工況下旋渦泵內流特性及聲場特性,在深入研究其內部流場和聲場基礎上,分析流場對聲場的影響,找到旋渦泵的主要噪聲源,探討旋渦泵內部流動誘導噪聲的傳播規(guī)律,在此基礎上采取一定的降噪優(yōu)化手段,并進行試驗驗證。

    1 計算模型及網(wǎng)格劃分

    1.1 計算模型

    選取浙江省某企業(yè)的一款旋渦自吸泵為研究對象,其主要設計參數(shù):流量Qd=2 m3/h ,揚程H=24 m,轉速n=2 900 r/min,葉片兩側交錯分布,葉片數(shù)Z=40,葉輪外徑D2=76 mm,葉輪出口寬度b2=8 mm,流道截面形狀為矩形,流道寬度B=16 mm。采用Pro/E 5.0軟件對旋渦泵進行三維全流場建模,如圖1所示。其計算域主要分為:泵體、葉輪、進口段、出口段4部分。為了減少進、出口可能存在的回流影響,使流場計算結果更加準確,對進口段水體和泵體出口段水體進行了5倍直徑的延長。

    圖1 泵體和葉輪計算域Fig.1 Computational domains of pump body and impeller

    1.2 計算網(wǎng)格及無關性檢驗

    由于旋渦泵幾何結構復雜,采用ANSYS ICEM對旋渦自吸泵各計算域進行非結構網(wǎng)格劃分。旋渦泵泵體隔舌處間隙較薄,該區(qū)域的網(wǎng)格質量對流動的影響較大,因此對泵體隔舌間隙處進行局部加密處理,如圖2所示,同時為了保證網(wǎng)格質量及計算準確性,對葉片表面、葉輪壁面以及葉輪與泵體交界面進行面加密處理。

    圖2 泵體葉輪網(wǎng)格Fig.2 Grids of pump body and impeller

    圖3 不同網(wǎng)格數(shù)下的揚程Fig.3 Head under different mesh elements

    在保證網(wǎng)格質量的前提下選用了7種不同網(wǎng)格尺寸對網(wǎng)格進行劃分,并進行網(wǎng)格無關性檢驗。如圖3所示,從圖中可以看出隨著總網(wǎng)格數(shù)的增加,模擬結果漸趨于穩(wěn)定,當網(wǎng)格總數(shù)達到150萬附近時,隨著網(wǎng)格數(shù)增加,模擬揚程變化很小,此時可以認為網(wǎng)格對計算結果無影響。綜合考慮計算所需資源以及準確性,最終確定本文采用的網(wǎng)格單元總數(shù)約180萬,其中葉輪網(wǎng)格數(shù)約為53萬,泵體網(wǎng)格數(shù)約為120萬,進出水管路各占約3.5萬。葉輪網(wǎng)格質量在0.36以上,泵體網(wǎng)格質量在0.32以上,葉片表面y+(第一層網(wǎng)格質心到壁面的無量綱距離)分布集中在18~348之間,葉輪壁面y+分布集中在17~320之間,由于近壁區(qū)采用的是Scalable壁面函數(shù)[13],且本文關注的是漩渦泵內部整體流動誘導產(chǎn)生的噪聲,著重研究湍流核心區(qū)內部流動對流動噪聲的影響,壁面剪切力以及粘性子層的數(shù)據(jù)并不是本文研究的重點,因此y+分布滿足計算要求[14]。

    2 數(shù)值模擬方法

    2.1 流場計算

    采用CFX對模型泵進行三維全流場瞬態(tài)場數(shù)值計算,由于RNGk-ε模型對湍動粘度進行了修正,考慮了平均流動中的旋轉及旋流流動情況,同時在ε方程中增加了反映主流的時均變率Eij,較標準k-ε模型能更好處理高應變率及流線彎曲程度較大的復雜流動,對壁面邊界層網(wǎng)格要求較低,考慮到計算能力,湍流模型選擇采用RNGk-ε模型[15],并采用有限體積法對控制方程進行二階精度的空間離散,時間離散采用全隱格式[16]。邊界條件設為壓力進口、質量流量出口,固體壁面采用無滑移邊界條件,近壁區(qū)采用Scalable壁面函數(shù)處理。定常計算時葉輪與泵體交界面設為凍結轉子模型。非定常計算葉輪與泵體交界面設為Transient rotor-state,同時為了加速收斂,在定常計算的基礎上進行非定常數(shù)值計算,共計算6個周期,定義葉輪每旋轉1°計算一次,非定常時間步長為Δt=5.747 126 4×10-5s,總計算時間t=0.124 137 91 s,當每個計算時間步長內的迭代次數(shù)達10次或控制方程變量的絕對殘差均小于10-4時進入下一個時間步迭代。

    為監(jiān)測各個工況下漩渦泵內部壓力脈動變化,在旋渦自吸泵隔舌處每隔5°設置一個監(jiān)測點,流道壁面設置7個監(jiān)測點,共13個監(jiān)測點,監(jiān)測點分布如圖4所示。

    圖4 監(jiān)測點分布圖Fig.4 Monitoring points distribution diagram of pump

    2.2 聲場計算

    提取非定常計算所得的速度、密度、壓力作為聲場計算的時域信號。利用聲學軟件ACTRAN進行聲場計算,圖5為旋渦泵聲學計算模型。

    圖5 聲學計算模型Fig.5 Acoustic computation model1.體源域 2.面源域 3.聲傳播區(qū)域 4.無反射邊界

    計算時,泵體為體聲源;葉輪與進口的聲源信息轉換至葉輪與泵體交界面處,設為面聲源域。在泵體進出口管路內各設置一個場點監(jiān)測聲壓級。為了保證計算準確,需保證聲學網(wǎng)格的單元最大長度不大于最大頻率處聲波波長的1/6[17],根據(jù)本文計算的旋渦泵最大頻率9 000 Hz及聲音在水中的傳播速度1 497 m/s,最終確定本文聲學網(wǎng)格長度為25 mm。由于漢寧窗函數(shù)[18]旁瓣泄漏少,計算時,為了減少頻譜泄漏與畸變,在傅里葉變換中加入漢寧窗函數(shù)。由于人為地對泵體進出口進行了聲源截斷,為了消除截面產(chǎn)生的虛假偶極子聲源影響,計算時加入濾波長度為0.1 mm的余弦濾波函數(shù)[19]。

    3 模擬結果分析

    3.1 外特性對比

    圖6為模擬結果與試驗結果對比圖,從圖中可以看出,模擬揚程和效率與試驗結果有一定的偏差,這主要是由于在數(shù)值計算過程中,忽略了由泄漏引起的容積損失和軸承、摩擦副等引起的機械損失。但總體來說模擬揚程和效率與試驗結果的變化趨勢基本一致,尤其是在小流量工況下吻合度較好,設計工況下計算值與試驗值相比,揚程誤差為3.9%,效率誤差為1.3%,說明流場數(shù)值計算結果準確度較高,具有一定的參考價值,其計算結果可以作為后續(xù)聲場計算的信號源。

    圖6 外特性曲線對比Fig.6 Comparison curves of pump performance

    3.2 壓力脈動分析

    選取0.4Qd、0.8Qd及Qd(Qd表示設計流量)3個流量點對旋渦泵的壓力脈動進行分析。本文研究的旋渦泵葉片數(shù)Z=40,轉速n=2 900 r/min,計算得到軸頻為48.33 Hz,葉頻為1 933.33 Hz。圖7為不同工況下監(jiān)測點的壓力脈動頻域特性。

    圖7 不同工況下監(jiān)測點壓力脈動頻域特性Fig.7 Frequency spectra of pressure fluctuation of monitoring points under different operating conditions

    由圖7可以看出,旋渦泵的壓力脈動由離散頻率和寬頻2部分組成,以離散頻率為主,且隨著流量的增加壓力脈動的幅值逐漸減小。在相同流量下,隔舌及流道壁面處各個監(jiān)測點的壓力脈動幅值均集中在葉頻、2倍葉頻、3倍葉頻、4倍葉頻處,2倍葉頻脈動幅值最大。流道壁面的壓力脈動幅值較小,遠小于隔舌處監(jiān)測的壓力脈動幅值,這主要是由于葉輪周期性地掃過隔舌,引起隔舌處較為劇烈的壓力脈動。從圖7c可以看出,在隔舌處的監(jiān)測點P1~P6,越靠近流道出口隔舌處的監(jiān)測點壓力脈動的幅值越大,P1處的脈動幅值最高,說明在設計工況下,流道出口隔舌處是旋渦泵葉輪與泵體隔舌動靜干涉最強的區(qū)域。通過對比可以發(fā)現(xiàn),流量越小,壓力脈動的程度越劇烈。在小流量工況下,P1~P6監(jiān)測點中,同樣是 P1監(jiān)測點壓力脈動幅值最大,但P1~P6監(jiān)測點不具備像設計工況下的壓力脈動漸變趨勢,這可能是由于小流量工況下旋渦泵內部流動非常紊亂導致規(guī)律性不太明顯。

    3.3 聲場計算結果分析

    3.3.1出口處的聲壓級

    聲壓級的計算公式為

    Spl=20lg(pe/pref)

    (1)

    式中pe——聲壓pref——參考聲壓

    水的參考聲壓一般取10-6Pa。為了研究旋渦泵輻射噪聲頻譜特性,給出了不同流量下旋渦泵出口場點P2的聲壓級頻譜圖,如圖8所示。從圖中可以看出,旋渦泵的流動噪聲主要由寬頻噪聲和離散噪聲2部分組成,不同工況下,最大聲壓級峰值均出現(xiàn)在2倍葉頻處,且在葉頻及其倍頻處聲壓級均有峰值出現(xiàn),呈現(xiàn)較強的離散性,與前文壓力脈動分析的結果一致,說明旋渦泵內部流動誘導噪聲的頻率特性與壓力脈動相似。隨著流量的增加,最大聲壓級的峰值逐漸減小,峰值依次為:186.19、185.12、184.52 dB,這與該泵實際運行時小流量工況下噪聲較大的現(xiàn)實相一致,同時也說明使旋渦泵在設計工況下運行是較好的控制旋渦泵噪聲的方法之一。

    3.3.2聲壓云圖

    由于旋渦泵內部流動誘導噪聲的峰值主要集中

    圖8 不同流量下場點P2聲壓級頻譜Fig.8 Sound pressure of field point P2 at different flow rates

    在葉頻及其倍頻下,而本文研究的主要目的是為了解決小流量工況下旋渦泵噪聲過大問題,因此給出小流量0.4Qd下旋渦泵泵體內葉頻及其倍頻處的聲壓云圖,如圖9所示。從圖中可以看出,在0.4Qd下,聲傳播區(qū)域的聲壓級基本不變,以平面波的形式向前傳播。2倍葉頻處的聲壓級最大,同一頻率下,泵體流道內聲壓級較大,其他位置的聲壓較小,說明泵體流道位置是主要的噪聲輻射源,從其他位置輻射的噪聲較小,說明旋渦泵內部動靜部件的干涉作用可能是旋渦泵內部流動誘導噪聲產(chǎn)生的主要原因之一。

    3.3.3聲源云圖分析

    為了更詳細地了解旋渦泵內部的主要噪聲源分布區(qū)域,對其內部的聲源分布進行分析,聲源場分布即Lighthill方程中聲源項的傅里葉變換模值分布,以聲壓級形式表示(pref=10-6Pa)。由于0.4Qd下2倍葉頻處的聲壓級最大,因此給出2倍葉頻處的體聲源和面聲源分布云圖,如圖10所示。從圖10a可以看出,體聲源主要分布在流道內部、隔舌處、泵體流道進口及出口處,隔舌位置體聲源密度最大,流道內次之,隔舌位置的體聲源主要是由于葉片與蝸殼隔舌的動靜干涉作用引起。流道內體聲源可能是由于液體沖擊隔舌、葉輪與泵體流道內液體能量交換引起的縱向和徑向漩渦等造成的。流道進口處的聲源可能是由于小流量工況下的不穩(wěn)定進流、液體沖擊葉輪導致的。從圖10b可以看出,流道內部面聲源主要集中在隔舌及葉輪流道頂部交接面處。葉輪流道頂部交接面處的聲源主要是由于葉輪與隔舌的動靜干涉引起葉輪流道內流動狀態(tài)發(fā)生較大變化,流體沖擊葉輪以及流出葉輪的液體與流道內的液體相互撞擊引起的。從圖中可以看出,在泵體隔舌處,越靠近流道出口,聲源分布越密集,這與前文壓力脈動分析結果一致,說明葉輪與隔舌的動靜干涉作用產(chǎn)生的不穩(wěn)定脈動力是旋渦泵流動誘導噪聲的主要噪聲源。

    圖9 0.4Qd不同頻率的聲壓級Fig.9 Sound pressure level at different frequencies of 0.4Qd

    圖10 體聲源和面聲源分布云圖Fig.10 Cloud charts of volume and area sound source distribution

    4 降噪優(yōu)化與試驗驗證

    4.1 降噪優(yōu)化

    在風機中經(jīng)常采用不等節(jié)距葉片降低葉輪離散噪聲的方法進行降噪處理。該方法主要是采用不等節(jié)距葉片破壞葉輪葉片均勻分布特性,使相鄰兩葉片撞擊隔舌的相位差不同,從而減小脈動強度,使基頻處的脈動峰值削弱和分散成峰值更小的分布[20-22]?;谇拔姆治?,旋渦泵流動誘導噪聲過大主要也是由于葉輪與隔舌的動靜干涉引起的不穩(wěn)定脈動力導致的。因此本文嘗試參考風機中不等節(jié)距葉片降噪的方法對模型泵進行降噪處理。本文采用經(jīng)典的EWALD頻率調制規(guī)律[23],采用正弦調制葉片分布方式,設計了4種調制角A分別為2°、4°、6°、8°的不等節(jié)距葉輪,A=0°為等節(jié)距葉輪。其葉片分布與周向葉距角分布如圖11、12所示,其分布公式為

    φ′n=φn+Asin(Nφn)

    (2)

    式中φn——等節(jié)距分布時第n個葉片應安裝的位置

    φ′n——不等節(jié)距分布時第n個葉片安裝的位置

    N——調制量的循環(huán)次數(shù),取1

    圖11 非等距葉片分布示意圖Fig.11 Schematic diagram of unequal spacing blade distribution

    圖12 葉片周向葉距角分布圖Fig.12 Blade peripheral intersection angle distribution

    在0.4Qd工況下對調制角為2°、4°、6°和8°的4種不同葉片分布方式的葉輪聲場進行數(shù)值計算,得到其泵體出口場點P2的最高聲壓級分別為:188.23、182.23、184.44、185.14 dB,可以看出,0.4Qd下調制角A為4°的不等節(jié)距葉輪降噪效果最好,泵體出口場點P2的最大聲壓級由186.19 dB降至182.23 dB。

    用3D打印方法做出調制角A=4°的實體葉輪,如圖13所示。在開式試驗臺上進行外特性和噪聲試驗,驗證模擬結果的可靠性。

    圖13 不等節(jié)距葉輪與試驗臺Fig.13 Unequal spacing blade impeller and pump testing platform

    4.2 外特性對比

    外特性對比結果如圖14所示,對比2種不同節(jié)距葉片分布方式的葉輪可以看出,等節(jié)距葉輪和調制角A=4°的不等節(jié)距葉輪,模擬揚程相差不大,隨流量變化趨勢一致,效率在小流量工況下基本一致,大流量工況下有一定誤差,但總體變化趨勢一致,吻合度較好,說明調制角A=4°的不等節(jié)距葉輪對旋渦泵外特性性能影響不大,尤其是在小流量工況下,揚程和效率基本與等節(jié)距葉輪一致。

    圖14 不等節(jié)距葉片對泵外特性的影響Fig.14 Effect of unequal spacing blade on performance of pump

    4.3 降噪效果驗證

    采用型號為日本RION理音NL-22型聲級計,選取0.4Qd和Qd2個工況點對2種不同葉片分布方式的模型泵進行噪聲測試,測試結果如表1所示。

    表1 不同流量下噪聲測試對比Tab.1 Comparison of experiment SPL atdifferent flows dB

    可以看出采用A=4°的不等節(jié)距葉片在設計工況和小流量工況下都有一定的降噪效果,說明將風機上經(jīng)常采用的不等節(jié)距葉片方法應用在旋渦泵上進行降噪是可行的,同樣能起到削弱和分散離散噪聲峰值的作用,達到降噪的目的。

    5 結論

    (1)旋渦泵流動誘導噪聲的產(chǎn)生主要是由葉輪與泵體隔舌之間的動靜干涉作用引起,流量越小,泵體出口聲壓級越大,噪聲源主要集中在泵體流道和隔舌位置,隔舌是主要的噪聲源,且流量越小聲源密度越大,越靠近流道出口隔舌處聲源密度越大。

    (2)旋渦泵內部離散噪聲傳播特性與壓力脈動情況基本一致,噪聲峰值與脈動強度密切相關,噪聲分布主要集中在葉頻及其倍頻處,2倍葉頻處聲壓級最大。

    (3)采用不等節(jié)距葉片對旋渦泵進行降噪的方法是可行的。合理調制角不僅對外特性性能影響不大且能削弱和分散離散噪聲峰值,達到降噪目的。調制角A=4°的不等節(jié)距葉輪對本文研究的旋渦泵外特性性能影響不大,且能使小流量0.4Qd下使噪聲下降2 dB,設計工況下使噪聲下降4 dB。

    1 趙王勇,張凡,王振,等.旋渦泵的研究現(xiàn)狀與發(fā)展方向[J].中國農村水利水電,2010(4):70-73.

    ZHAO Wangyong,ZHANG Fan,WANG Zhen,et al.Research status and research prospect of vortex pumps[J].China Rural Water and Hydropower,2010(4):70-73.(in Chinese)

    2 SONG J W,ENGEDA A,CHUNG M K. A modified theory for the flow mechanism in regenerative flow pump[J].Proc. IMechE, Part A:Journal of Power and Energy,2003,217(3):311-322.

    3 王洋,李亞成,曹璞鈺,等.M型截面流道對旋渦泵性能的影響分析[J/OL].農業(yè)機械學報,2014,45(3):77-82.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20140314amp;flag=1amp;journal_id=jcsam.DOI: 10.6041/j.issn.1000-1298.2014.03.014.

    WANG Yang,LI Yacheng,CAO Puyu,et al.Effects of channel with M-type sectionon the performance of vortex pump[J/OL]. Transactions of the Chinese Society for Agricultural Machinery,2014,45(3):77-82. (in Chinese)

    4 CHOI W C, YOO I S, MU R P, et al. Experimental study on the effect of blade angle on regenerative pump performance[J]. Proc. IMechE, Part A:Journal of Power amp; Energy, 2013, 227(5):585-592.

    5 孔繁余, 周水清, 王志強,等. 漩渦泵內部不穩(wěn)定流場數(shù)值模擬[J]. 農業(yè)機械學報, 2011, 42(1):44-48.

    KONG Fanyu,ZHOU Shuiqing,WANG Zhiqiang,et al. Numerical simulation of internal unsteady flow in vortex pump[J].Transactions of the Chinese Society for Agricultural Machinery, 2011, 42(1):44-48.(in Chinese)

    6 QUAIL F J,SCANLON T,BAUMGARTNER A.Design study of a regenerative pump using one-dimensional and three dimensional numerical techniques[J].European Journal of Mechanics B/Fluids,2012,31:187.

    7 TIMUSHEV S. Development and experimental validation of 3d acoustic-vortex numerical procedure for centrifugal pump noise prediction[C]∥ASME 2009 Fluids Engineering Division Summer Meeting, 2009:389-398.

    8 BAILLY C,BOGEY C,GLOERFELT X. Some useful hybrid approaches for predicting aerodynamic noise[J].Comptes Rendus Mecanique,2005,333(9):666-675.

    9 閻超, 于劍, 徐晶磊,等. CFD模擬方法的發(fā)展成就與展望[J]. 力學進展, 2011, 41(5):562-589.

    YAN Chao,YU Jian, XU Jinglei, et al.On the achievements and prospects for the methods of coputational fluid dynamics[J].Advance in Mechanics, 2011, 41(5):562-589.(in Chinese)

    10 曹璞鈺,印剛,王洋,等.離心泵內雙龍卷風式分離渦數(shù)值分析 [J/OL].農業(yè)機械學報 ,2016, 47(4):22-28.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20160404amp;flag=1amp;journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2016.04.004.

    CAO Puyu,YIN Gang,WANG Yang,et al.Numerical analysis of double tornado-type separation vortices in centrifugal pump[J/OL].Transactions of the Chinese Society for Agricultural Machinery ,2016, 47(4):22-28. (in Chinese)

    11 司喬瑞, 袁建平, 衡亞光,等. 蝸殼形狀對離心泵流動誘導噪聲的影響[J]. 排灌機械工程學報, 2015, 33(3):209-215.

    SI Qiaorui, YUAN Jianping, HENG Yaguang, et al. Effects of cross-section and cut-water shapes of volute on flow induced noise in centrifugal pumps[J]. Journal of Drainage and Irrigation Machinery Engineering, 2015, 33(3):209-215.(in Chinese)

    12 談明高,張景,劉厚林,等.葉片不等間距對離心泵水動力噪聲的影響[J/OL].農業(yè)機械學報,2016,47(2):22-27,34.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20160203amp;flag=1amp;journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2016.02.003.

    TAN Minggao,ZHANG Jing,LIU Houlin, et al.Effect of uneven spacing blade on hydrodynamic noise of centrifugal pump[J/OL].Transactions of the Chinese Society for Agricultural Machinery,2016,47(2):22-27,34. (in Chinese)

    13 ANSYS Inc. ANSYS CFX-solver modeling guide, release 14.5[M]. ANSYS Inc.,2012.

    14 王福軍. 計算流體動力學分析:CFD軟件原理與應用[M]. 北京:清華大學出版社, 2004.

    15 李貴東,王洋,曹璞鈺,等.射流式離心泵非設計工況下內部流動研究 [J/OL]. 農業(yè)機械學報, 2015, 46(8):48-53.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20150808amp;flag=1amp;journal_id=jcsam.DOI:10.6041/j.issn.1000-1298.2015.08.008.

    LI Guidong,WANG Yang,CAO Puyu, et al. Internal flow of flow-ejecting centrifugal pump under off-design conditions [J/OL] .Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(8):48-53.(in Chinese)

    16 胡坤. ANSYS ICEM CFD工程實例詳解[M]. 北京:人民郵電出版社, 2014.

    17 王治國. MSC.ACTRAN工程聲學有限元分析理論與應用[M]. 北京:國防工業(yè)出版社, 2007.

    18 汪偉,謝皓辰,梁光明,等.加窗離散傅里葉變換性能分析和比對[J].現(xiàn)代電子技術,2012,35(3):115-118.

    WANG Wei,XIE Haochen,LIANG Guangming,et al. Performance analysis and comparison of windowed discrete fourier transform[J].Modern Electronics Technique,2012,35(3):115-118. (in Chinese)

    19 趙威. 離心泵的流噪聲數(shù)值計算及聲優(yōu)化[D]. 武漢:華中科技大學, 2013.

    20 馬健峰,袁民建,劉秋洪,等.不等距葉片離心風機氣動噪聲的數(shù)值與實驗研究[J].噪聲與振動控制,2008(3):100-103.

    MA Jianfeng,YUAN Minjian,LIU Qiuhong, et al.Experimental and numerical study of aeroacoustics of a centrifugal fan with non-isometric blade impeller[J].Noise and Vibration Control,2008(3):100-103.(in Chinese)

    21 孫曉峰. 不等節(jié)距葉片風機氣動聲學特性的研究[J]. 北京航空航天大學學報, 1986(4):144-152.

    SUN Xiaofeng. The aeroacoustic nature of unequally spaced fan[J].Journal of Beijing Institute of Aeroacoustic and Astronautics, 1986(4):144-152.(in Chinese)

    22 張瑜. 貫流風機內流及氣動噪聲分析[D]. 武漢:華中科技大學, 2007.

    23 EWALD D, PAVLOVIC A, BOLLINGER J G. Noise reduction by applying modulation principles[J]. Journal of the Acoustical Society of America, 1971, 49(5):1381-1385.

    Flow-inducedNoiseofSelf-primingVortexPumpatLow-flowRate

    WANG Yang HU Rixin TANG Haitao WANG Qun YU Haoqian JIA Jing

    (ResearchCenterofFluidMachineryEngineeringandTechnology,JiangsuUniversity,Zhenjiang212013,China)

    In order to research the internal sound field characteristics and reduce the running noise of self-priming vortex pump at low-flow rate, the CFD and Lighthill sound analogy theory was adopted to solve the internal acoustic field and optimize the noise of self-priming vortex pump. Firstly, the pressure fluctuation of self-priming vortex at 0.4Qd,0.8QdandQdwere obtained from the unsteady simulation with the RNGk-εturbulence model in CFX software. The datum of pressure fluctuation were extracted and input into the ACTRAN acoustics software, the SPL of the pump discharge, internal sound pressure distribution and main noise source in pumps at different flows were got by acoustic field calculation. Then,according to the laws of the Ewald frequency modulation, four different types of unequal spacing impellers were designed with modulation angleAof 2°, 4°,6°and 8°. Finally, the acoustic field calculation of four different impellers was conducted, the impeller with unequal spacing andAvalue of 4°had lower noise, whose sample made by 3D printing was chosen to conduct the performance and noise reduction test. The simulation results showed that the flow-induced noise of self-priming vortex pump was closely related to the pressure fluctuation which mainly caused by rotor-stator interaction,both of them had the same frequency characteristic,the maximum sound pressure appears at two times blade passing frequency. With the decrease of the flow rate, the SPL of the pump discharge was increased, the sound source of the self-priming vortex pump was mainly distributed in the channel and the tongue close to the pump discharge. The noise test results showed that the impeller sample withAvalue of 4°had little effect on the performance of the self-priming vortex pump and can reduce noise at low-flow rate, under the design point, noise can be decreased by 4 dB, and it can be decreased by 2 dB when low flow was 0.4Qd.

    self-priming vortex pump; pressure fluctuation; flow-induced noise; unequal spacing blade; numerical simulation

    10.6041/j.issn.1000-1298.2017.11.023

    TH314

    A

    1000-1298(2017)11-0188-08

    2017-02-27

    2017-04-01

    國家高技術研究發(fā)展計劃(863計劃)項目(2011AA100506)、國家自然科學基金項目(51379091)和江蘇高校優(yōu)勢學科建設工程項目(PAPD)

    王洋(1955—),男,研究員,博士生導師,主要從事流體機械及工程研究,E-mail: pgwy@ujs.edu.cn

    猜你喜歡
    節(jié)距泵體脈動
    新學期,如何“脈動回來”?
    家教世界(2023年25期)2023-10-09 02:11:56
    基于ANSYS Workbench三螺桿軸頭泵泵體設計與試驗
    RBI在超期服役脈動真空滅菌器定檢中的應用
    圓柱螺旋扭轉彈簧節(jié)距對剛度與應力的影響研究*
    機械制造(2021年9期)2021-10-12 12:28:08
    輪胎花紋節(jié)距噪聲研究及其應用
    單級蝸殼泵強度有限元分析
    泵體水力模具逆向工程技術研究
    地球脈動(第一季)
    一種可移動的充氣橡膠制品裝置
    基于Autolisp的輪胎花紋節(jié)距自動排列
    99精品久久久久人妻精品| 黑人巨大精品欧美一区二区蜜桃| 国产99久久九九免费精品| 最黄视频免费看| 亚洲成人国产一区在线观看 | 黄色视频在线播放观看不卡| 国产精品一国产av| 日韩一区二区视频免费看| 午夜激情久久久久久久| 十八禁人妻一区二区| 蜜桃在线观看..| 在线免费观看不下载黄p国产| 91精品国产国语对白视频| 日本午夜av视频| 成人亚洲欧美一区二区av| 欧美日韩亚洲综合一区二区三区_| 午夜福利一区二区在线看| netflix在线观看网站| 卡戴珊不雅视频在线播放| 亚洲国产欧美一区二区综合| 欧美人与善性xxx| 中文字幕制服av| 最近2019中文字幕mv第一页| 老司机深夜福利视频在线观看 | 中文字幕av电影在线播放| 最新在线观看一区二区三区 | 777久久人妻少妇嫩草av网站| 韩国高清视频一区二区三区| 青春草亚洲视频在线观看| 街头女战士在线观看网站| 欧美日韩一级在线毛片| 黄片无遮挡物在线观看| 少妇精品久久久久久久| 免费少妇av软件| 秋霞在线观看毛片| 亚洲欧美中文字幕日韩二区| 色婷婷久久久亚洲欧美| 大片电影免费在线观看免费| 欧美激情高清一区二区三区 | 在线观看国产h片| 一级毛片电影观看| 亚洲国产av影院在线观看| 欧美日韩亚洲综合一区二区三区_| 老鸭窝网址在线观看| 国产精品久久久人人做人人爽| 国产伦理片在线播放av一区| 久热这里只有精品99| 国产伦理片在线播放av一区| 亚洲精品中文字幕在线视频| 韩国精品一区二区三区| 亚洲欧美中文字幕日韩二区| 熟妇人妻不卡中文字幕| 欧美 亚洲 国产 日韩一| 欧美成人午夜精品| 亚洲av成人精品一二三区| 日日撸夜夜添| 中国国产av一级| 精品一区二区三区四区五区乱码 | videosex国产| 亚洲成人国产一区在线观看 | 亚洲国产精品国产精品| av在线老鸭窝| 青春草视频在线免费观看| 少妇人妻精品综合一区二区| 下体分泌物呈黄色| av电影中文网址| 一边摸一边抽搐一进一出视频| 精品一区二区三区四区五区乱码 | 99国产综合亚洲精品| av有码第一页| 国产不卡av网站在线观看| 亚洲精品国产区一区二| 一级片'在线观看视频| 免费av中文字幕在线| 午夜激情久久久久久久| 久久久久久久久久久免费av| 赤兔流量卡办理| 久久国产精品大桥未久av| 亚洲,欧美,日韩| 欧美日本中文国产一区发布| 国产一区二区激情短视频 | av不卡在线播放| 日韩大片免费观看网站| 男女高潮啪啪啪动态图| 交换朋友夫妻互换小说| 别揉我奶头~嗯~啊~动态视频 | 亚洲四区av| 国产精品一区二区精品视频观看| 男女国产视频网站| 中文字幕av电影在线播放| 欧美日韩视频高清一区二区三区二| av女优亚洲男人天堂| 大片免费播放器 马上看| 日韩 亚洲 欧美在线| 悠悠久久av| 免费黄色在线免费观看| 悠悠久久av| 一边摸一边抽搐一进一出视频| 婷婷色综合www| 婷婷成人精品国产| 国产精品一二三区在线看| 国产爽快片一区二区三区| 国产亚洲av片在线观看秒播厂| 日韩一区二区视频免费看| 蜜桃在线观看..| 久久久久久久久久久免费av| 赤兔流量卡办理| 精品国产露脸久久av麻豆| 久久久久精品人妻al黑| 久久久亚洲精品成人影院| 大陆偷拍与自拍| 日韩人妻精品一区2区三区| 免费在线观看完整版高清| 黄色一级大片看看| 性高湖久久久久久久久免费观看| 老汉色av国产亚洲站长工具| 在线看a的网站| 日本wwww免费看| av又黄又爽大尺度在线免费看| 久久久精品国产亚洲av高清涩受| 成人漫画全彩无遮挡| 可以免费在线观看a视频的电影网站 | 国产黄色免费在线视频| 久久女婷五月综合色啪小说| 在线天堂最新版资源| 亚洲在久久综合| 国产av精品麻豆| 老司机影院毛片| 99国产综合亚洲精品| 欧美亚洲 丝袜 人妻 在线| 亚洲人成电影观看| xxx大片免费视频| 久久精品久久精品一区二区三区| 我要看黄色一级片免费的| 高清欧美精品videossex| 一边摸一边抽搐一进一出视频| 精品国产超薄肉色丝袜足j| 日韩制服丝袜自拍偷拍| 一二三四在线观看免费中文在| 国产精品国产三级国产专区5o| 天堂俺去俺来也www色官网| 午夜激情av网站| 老熟女久久久| a级毛片黄视频| 十八禁网站网址无遮挡| 亚洲成人av在线免费| 男女下面插进去视频免费观看| 国产黄色免费在线视频| 一区二区三区精品91| 亚洲一码二码三码区别大吗| 免费高清在线观看视频在线观看| 欧美老熟妇乱子伦牲交| 香蕉丝袜av| 久久国产精品大桥未久av| 精品少妇黑人巨大在线播放| 免费观看av网站的网址| 久久天躁狠狠躁夜夜2o2o | 欧美日韩福利视频一区二区| 最新在线观看一区二区三区 | 国产精品香港三级国产av潘金莲 | 巨乳人妻的诱惑在线观看| 午夜福利乱码中文字幕| 久久精品aⅴ一区二区三区四区| 人妻一区二区av| 久久人人爽人人片av| 伊人久久大香线蕉亚洲五| 日本91视频免费播放| 精品人妻熟女毛片av久久网站| 亚洲欧美成人精品一区二区| 精品国产一区二区三区四区第35| 19禁男女啪啪无遮挡网站| 国产成人精品在线电影| 日本黄色日本黄色录像| 纵有疾风起免费观看全集完整版| 午夜日本视频在线| 国产一区二区三区av在线| 777米奇影视久久| 青青草视频在线视频观看| 一级毛片我不卡| 97精品久久久久久久久久精品| 亚洲精品在线美女| 亚洲成av片中文字幕在线观看| 日韩熟女老妇一区二区性免费视频| 丝袜美腿诱惑在线| 国产精品.久久久| 一区二区三区精品91| 伊人亚洲综合成人网| 精品国产一区二区久久| 又大又黄又爽视频免费| 精品一区二区三区av网在线观看 | 亚洲精品乱久久久久久| av网站在线播放免费| 少妇人妻久久综合中文| 国产成人91sexporn| 亚洲精品久久午夜乱码| 狂野欧美激情性bbbbbb| 午夜福利影视在线免费观看| 欧美成人午夜精品| av女优亚洲男人天堂| 中文字幕制服av| 久久久久精品人妻al黑| 天天操日日干夜夜撸| 天堂中文最新版在线下载| 中国三级夫妇交换| 青青草视频在线视频观看| 亚洲av电影在线观看一区二区三区| 一区二区三区乱码不卡18| 日日摸夜夜添夜夜爱| 99久久人妻综合| 欧美最新免费一区二区三区| 亚洲av日韩在线播放| 99久久精品国产亚洲精品| 一级毛片电影观看| 狂野欧美激情性bbbbbb| 久久久国产欧美日韩av| 精品国产露脸久久av麻豆| 老鸭窝网址在线观看| 丰满乱子伦码专区| 夜夜骑夜夜射夜夜干| 飞空精品影院首页| 日韩熟女老妇一区二区性免费视频| 18在线观看网站| 国产亚洲午夜精品一区二区久久| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品久久精品一区二区三区| 97精品久久久久久久久久精品| 国产精品久久久久久精品电影小说| 极品少妇高潮喷水抽搐| 亚洲欧美色中文字幕在线| 免费观看人在逋| 久久久久久久久久久免费av| 在线观看免费午夜福利视频| 国产人伦9x9x在线观看| 久久久久网色| 一个人免费看片子| 人人妻人人澡人人爽人人夜夜| 波多野结衣一区麻豆| 中文字幕色久视频| 哪个播放器可以免费观看大片| av国产精品久久久久影院| 黑人欧美特级aaaaaa片| av在线播放精品| 国产精品久久久久久精品电影小说| 美女国产高潮福利片在线看| 亚洲视频免费观看视频| 久久亚洲国产成人精品v| 大陆偷拍与自拍| 另类亚洲欧美激情| 国产av一区二区精品久久| 狠狠婷婷综合久久久久久88av| 国产麻豆69| 亚洲国产欧美在线一区| 国产精品熟女久久久久浪| 18禁观看日本| 久久亚洲国产成人精品v| 国产黄频视频在线观看| 我的亚洲天堂| 国产老妇伦熟女老妇高清| tube8黄色片| 免费看不卡的av| 中文欧美无线码| 青青草视频在线视频观看| 嫩草影院入口| 国产亚洲av片在线观看秒播厂| 一级爰片在线观看| 成人毛片60女人毛片免费| 日本色播在线视频| 观看美女的网站| 丝袜在线中文字幕| 男女无遮挡免费网站观看| 亚洲一区二区三区欧美精品| 久久这里只有精品19| 中文字幕av电影在线播放| 久久精品久久久久久久性| 熟女av电影| 欧美精品一区二区大全| 精品久久久精品久久久| 精品久久久精品久久久| 免费在线观看视频国产中文字幕亚洲 | 黄色 视频免费看| 91精品国产国语对白视频| 午夜福利视频在线观看免费| av在线app专区| 成人亚洲欧美一区二区av| 久久国产精品大桥未久av| 婷婷色综合大香蕉| 久久久久精品国产欧美久久久 | av片东京热男人的天堂| 卡戴珊不雅视频在线播放| 卡戴珊不雅视频在线播放| 精品亚洲成a人片在线观看| 999精品在线视频| 国产一区二区激情短视频 | 最近手机中文字幕大全| 国产野战对白在线观看| 午夜福利视频在线观看免费| 国产精品一区二区在线观看99| 亚洲一码二码三码区别大吗| 中文字幕人妻熟女乱码| 黄片播放在线免费| 久久精品久久久久久久性| 深夜精品福利| 日韩熟女老妇一区二区性免费视频| 国产精品亚洲av一区麻豆 | 午夜福利在线免费观看网站| 狠狠婷婷综合久久久久久88av| 欧美另类一区| 亚洲国产欧美网| 18禁观看日本| 两个人免费观看高清视频| 亚洲成色77777| 欧美精品一区二区大全| 大香蕉久久成人网| 大码成人一级视频| 久久久久网色| 亚洲国产欧美日韩在线播放| av有码第一页| 美女高潮到喷水免费观看| 各种免费的搞黄视频| 日韩欧美精品免费久久| 亚洲一区中文字幕在线| 日韩欧美一区视频在线观看| 啦啦啦 在线观看视频| 欧美日韩视频精品一区| 美女福利国产在线| 免费看av在线观看网站| 国产免费现黄频在线看| 亚洲色图 男人天堂 中文字幕| 9色porny在线观看| 久久 成人 亚洲| 日本91视频免费播放| 看免费成人av毛片| 精品国产乱码久久久久久小说| 国产精品亚洲av一区麻豆 | 日本91视频免费播放| 亚洲av中文av极速乱| 久久99热这里只频精品6学生| 高清欧美精品videossex| 免费少妇av软件| 一区福利在线观看| 亚洲av国产av综合av卡| 桃花免费在线播放| 国产精品免费大片| 日本av手机在线免费观看| 国产欧美日韩一区二区三区在线| 女人久久www免费人成看片| 少妇人妻久久综合中文| 美女国产高潮福利片在线看| 欧美黑人欧美精品刺激| 少妇被粗大猛烈的视频| av不卡在线播放| 黄片小视频在线播放| 欧美人与性动交α欧美精品济南到| 亚洲成人手机| 亚洲国产看品久久| 男女边摸边吃奶| 中文字幕高清在线视频| 搡老岳熟女国产| 国产精品一国产av| 亚洲国产欧美网| 欧美变态另类bdsm刘玥| 波多野结衣一区麻豆| 91精品国产国语对白视频| 久久精品aⅴ一区二区三区四区| 这个男人来自地球电影免费观看 | 尾随美女入室| 日本wwww免费看| 美女国产高潮福利片在线看| 尾随美女入室| 精品久久久久久电影网| 亚洲欧美一区二区三区国产| 五月开心婷婷网| 欧美在线一区亚洲| 国产激情久久老熟女| 美女视频免费永久观看网站| av视频免费观看在线观看| 波多野结衣一区麻豆| 欧美国产精品va在线观看不卡| 美女中出高潮动态图| 亚洲av国产av综合av卡| 国产成人av激情在线播放| 午夜福利视频在线观看免费| 国产精品成人在线| 亚洲综合精品二区| 成年人免费黄色播放视频| 日韩av免费高清视频| 亚洲伊人色综图| 汤姆久久久久久久影院中文字幕| 少妇猛男粗大的猛烈进出视频| 日韩一区二区视频免费看| av片东京热男人的天堂| 亚洲欧美激情在线| 中文字幕人妻丝袜制服| 精品第一国产精品| 在线观看人妻少妇| 一区二区三区精品91| 亚洲精品日本国产第一区| www.av在线官网国产| 激情五月婷婷亚洲| 男女午夜视频在线观看| 久久热在线av| 亚洲专区中文字幕在线 | 免费日韩欧美在线观看| 国产欧美亚洲国产| 国产精品欧美亚洲77777| 亚洲国产精品国产精品| 大香蕉久久成人网| 在现免费观看毛片| 午夜福利在线免费观看网站| 国产精品女同一区二区软件| 国产精品香港三级国产av潘金莲 | 午夜激情久久久久久久| 如何舔出高潮| 久久天躁狠狠躁夜夜2o2o | 少妇 在线观看| 男女之事视频高清在线观看 | 精品一区二区三区四区五区乱码 | 亚洲一区中文字幕在线| 久久久久久久久久久免费av| 人人澡人人妻人| 日韩精品免费视频一区二区三区| 99精品久久久久人妻精品| 久久人人爽人人片av| 叶爱在线成人免费视频播放| 久久影院123| 亚洲欧美色中文字幕在线| 国产麻豆69| 亚洲精品中文字幕在线视频| 亚洲欧美一区二区三区黑人| 精品国产国语对白av| 不卡视频在线观看欧美| 国产人伦9x9x在线观看| 国产免费又黄又爽又色| e午夜精品久久久久久久| 国产成人欧美| 日韩人妻精品一区2区三区| 在线 av 中文字幕| 色网站视频免费| 丁香六月欧美| 亚洲精品一二三| 日韩大片免费观看网站| 亚洲国产精品国产精品| 婷婷色综合大香蕉| 免费看av在线观看网站| 麻豆av在线久日| 男人爽女人下面视频在线观看| 99国产综合亚洲精品| 欧美亚洲 丝袜 人妻 在线| 一本久久精品| 性高湖久久久久久久久免费观看| 综合色丁香网| 免费看不卡的av| 欧美另类一区| 国产成人系列免费观看| 国产成人免费观看mmmm| 如何舔出高潮| 日本av免费视频播放| 亚洲欧美日韩另类电影网站| 老鸭窝网址在线观看| 人妻人人澡人人爽人人| 国产成人精品在线电影| 久久99热这里只频精品6学生| 午夜免费观看性视频| 国产一区亚洲一区在线观看| 国产精品久久久久久精品古装| 久久99一区二区三区| 亚洲精品av麻豆狂野| 不卡视频在线观看欧美| 少妇猛男粗大的猛烈进出视频| 国产在线免费精品| 日日摸夜夜添夜夜爱| 老司机影院毛片| 久久韩国三级中文字幕| 啦啦啦 在线观看视频| 日本爱情动作片www.在线观看| 国产精品国产av在线观看| 国产精品三级大全| 国产免费又黄又爽又色| 一本色道久久久久久精品综合| 中文字幕制服av| 免费观看a级毛片全部| 黄网站色视频无遮挡免费观看| 亚洲一区二区三区欧美精品| 中文乱码字字幕精品一区二区三区| 色94色欧美一区二区| av国产精品久久久久影院| 国产伦理片在线播放av一区| 男女高潮啪啪啪动态图| 久久毛片免费看一区二区三区| 久久久久精品国产欧美久久久 | √禁漫天堂资源中文www| 无遮挡黄片免费观看| 亚洲欧美成人综合另类久久久| 纯流量卡能插随身wifi吗| 亚洲一区二区三区欧美精品| 一级a爱视频在线免费观看| 欧美精品一区二区大全| 欧美国产精品va在线观看不卡| 国产有黄有色有爽视频| 精品少妇内射三级| 捣出白浆h1v1| 久久人人爽av亚洲精品天堂| 一区二区日韩欧美中文字幕| 久久午夜综合久久蜜桃| 少妇被粗大的猛进出69影院| 成人影院久久| 国产成人午夜福利电影在线观看| 青青草视频在线视频观看| 国产探花极品一区二区| 黄色怎么调成土黄色| 精品少妇黑人巨大在线播放| 午夜福利视频精品| 国产成人精品在线电影| 欧美久久黑人一区二区| 黄色视频在线播放观看不卡| 伦理电影大哥的女人| 一级毛片我不卡| 美女脱内裤让男人舔精品视频| 亚洲自偷自拍图片 自拍| 天堂中文最新版在线下载| 亚洲精品日本国产第一区| 国产免费福利视频在线观看| 亚洲精品日本国产第一区| 日韩av不卡免费在线播放| 日韩一区二区三区影片| 免费高清在线观看视频在线观看| 日韩一区二区三区影片| 1024视频免费在线观看| 黄网站色视频无遮挡免费观看| 国产黄色免费在线视频| 777久久人妻少妇嫩草av网站| 日韩一本色道免费dvd| a级毛片黄视频| 精品福利永久在线观看| 女人爽到高潮嗷嗷叫在线视频| 欧美日韩亚洲综合一区二区三区_| 亚洲欧美一区二区三区久久| 国产有黄有色有爽视频| 精品卡一卡二卡四卡免费| 久久久国产精品麻豆| 亚洲 欧美一区二区三区| 美女高潮到喷水免费观看| 最近中文字幕高清免费大全6| 成年人免费黄色播放视频| 91精品国产国语对白视频| 国产成人精品在线电影| 久久97久久精品| 亚洲精品久久成人aⅴ小说| 欧美国产精品一级二级三级| 999久久久国产精品视频| 美女主播在线视频| 欧美 日韩 精品 国产| 欧美日韩一级在线毛片| 新久久久久国产一级毛片| 成人国语在线视频| h视频一区二区三区| 日韩av免费高清视频| 国产精品久久久人人做人人爽| 免费久久久久久久精品成人欧美视频| 欧美变态另类bdsm刘玥| 久久久久精品国产欧美久久久 | 嫩草影视91久久| 精品一区二区三区四区五区乱码 | 国产精品一区二区精品视频观看| 1024视频免费在线观看| 国产亚洲av高清不卡| 亚洲色图 男人天堂 中文字幕| 国产成人欧美| 国产麻豆69| 免费黄频网站在线观看国产| 成人国语在线视频| 久热这里只有精品99| 欧美另类一区| 男女床上黄色一级片免费看| 黄片播放在线免费| 欧美激情极品国产一区二区三区| 亚洲国产欧美网| 大话2 男鬼变身卡| 成年人午夜在线观看视频| 亚洲人成电影观看| 婷婷成人精品国产| 成人国产麻豆网| tube8黄色片| 欧美中文综合在线视频| 欧美黑人欧美精品刺激| 久久青草综合色| 黑丝袜美女国产一区| 免费在线观看视频国产中文字幕亚洲 | 国产亚洲欧美精品永久| 考比视频在线观看| 9色porny在线观看| 午夜福利免费观看在线| 亚洲精品国产av成人精品| 久久人妻熟女aⅴ| 欧美日韩av久久| 波多野结衣一区麻豆| 十八禁高潮呻吟视频| 国产成人免费观看mmmm| 免费少妇av软件| 美女扒开内裤让男人捅视频| 1024香蕉在线观看| 观看美女的网站| 99精国产麻豆久久婷婷| 不卡av一区二区三区| 亚洲精品视频女| 婷婷成人精品国产| 日韩一卡2卡3卡4卡2021年| 欧美精品高潮呻吟av久久| 多毛熟女@视频| 亚洲成人国产一区在线观看 | 国产av国产精品国产| 国产男女超爽视频在线观看| 老司机亚洲免费影院| 亚洲精品自拍成人| 久久久久久人人人人人|