• <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é)距自動排列
    国产精品久久久av美女十八| 在线看a的网站| 日本午夜av视频| 18禁裸乳无遮挡动漫免费视频| 女人精品久久久久毛片| 精品久久久久久电影网| 一级,二级,三级黄色视频| 亚洲一码二码三码区别大吗| 亚洲av日韩在线播放| 国产精品国产三级专区第一集| 别揉我奶头~嗯~啊~动态视频 | 亚洲精品美女久久av网站| 国产在线视频一区二区| 肉色欧美久久久久久久蜜桃| 大话2 男鬼变身卡| 国产色视频综合| 国产精品秋霞免费鲁丝片| 久久久久久久精品精品| 丝袜在线中文字幕| 国产欧美日韩精品亚洲av| 这个男人来自地球电影免费观看| 两个人看的免费小视频| 久久人人爽人人片av| 精品一品国产午夜福利视频| 亚洲精品日韩在线中文字幕| 国产一区二区 视频在线| 日韩 亚洲 欧美在线| 青春草视频在线免费观看| 天堂中文最新版在线下载| 在线av久久热| 亚洲一区中文字幕在线| 日韩视频在线欧美| 日韩一本色道免费dvd| 熟女少妇亚洲综合色aaa.| 中文字幕精品免费在线观看视频| 新久久久久国产一级毛片| 亚洲av电影在线进入| 国产国语露脸激情在线看| 亚洲,一卡二卡三卡| 婷婷色av中文字幕| 五月开心婷婷网| 国产极品粉嫩免费观看在线| 一级毛片我不卡| 97在线人人人人妻| 国产熟女欧美一区二区| 色综合欧美亚洲国产小说| 我要看黄色一级片免费的| 国产亚洲av片在线观看秒播厂| 亚洲人成77777在线视频| 视频区图区小说| 99国产精品99久久久久| 亚洲专区中文字幕在线| 久久人妻福利社区极品人妻图片 | 日本欧美视频一区| 久久久亚洲精品成人影院| www.精华液| 五月天丁香电影| 国产日韩欧美在线精品| 男女免费视频国产| 国产精品一区二区在线不卡| 亚洲精品国产av蜜桃| 人妻一区二区av| 欧美精品一区二区免费开放| 校园人妻丝袜中文字幕| 嫩草影视91久久| 国产免费现黄频在线看| xxx大片免费视频| www日本在线高清视频| 国产日韩欧美在线精品| 久久久亚洲精品成人影院| 男女高潮啪啪啪动态图| 亚洲av综合色区一区| 99国产精品免费福利视频| 后天国语完整版免费观看| xxxhd国产人妻xxx| 国产野战对白在线观看| 日韩av在线免费看完整版不卡| 免费av中文字幕在线| 久热爱精品视频在线9| 赤兔流量卡办理| 久久久亚洲精品成人影院| 国产人伦9x9x在线观看| 婷婷成人精品国产| 一边摸一边做爽爽视频免费| 日本午夜av视频| 成年人黄色毛片网站| 视频在线观看一区二区三区| 欧美日韩福利视频一区二区| 婷婷丁香在线五月| 日日夜夜操网爽| 亚洲精品一卡2卡三卡4卡5卡 | 久久99热这里只频精品6学生| 中文乱码字字幕精品一区二区三区| 免费看av在线观看网站| 日韩视频在线欧美| 性高湖久久久久久久久免费观看| 老鸭窝网址在线观看| 日本一区二区免费在线视频| 日本欧美国产在线视频| 国产一卡二卡三卡精品| 最黄视频免费看| 国产成人精品久久二区二区免费| 国产精品香港三级国产av潘金莲 | 黄色 视频免费看| 久久天堂一区二区三区四区| 人体艺术视频欧美日本| 大陆偷拍与自拍| 国产国语露脸激情在线看| 亚洲国产欧美在线一区| 成人三级做爰电影| 欧美国产精品一级二级三级| 制服人妻中文乱码| 午夜福利视频精品| 免费日韩欧美在线观看| www日本在线高清视频| 国产精品香港三级国产av潘金莲 | 麻豆乱淫一区二区| 成人亚洲精品一区在线观看| 熟女少妇亚洲综合色aaa.| 国产视频一区二区在线看| 十八禁网站网址无遮挡| 观看av在线不卡| 午夜福利视频精品| 两个人免费观看高清视频| 精品免费久久久久久久清纯 | 性色av一级| 欧美性长视频在线观看| 在线观看一区二区三区激情| 免费高清在线观看视频在线观看| 99九九在线精品视频| 亚洲国产毛片av蜜桃av| 亚洲国产欧美一区二区综合| 高清不卡的av网站| 久久精品熟女亚洲av麻豆精品| 汤姆久久久久久久影院中文字幕| 国产成人a∨麻豆精品| 波多野结衣av一区二区av| 人成视频在线观看免费观看| 制服人妻中文乱码| 亚洲成人免费av在线播放| 永久免费av网站大全| 久久久精品94久久精品| 亚洲伊人色综图| 亚洲一码二码三码区别大吗| 啦啦啦视频在线资源免费观看| 国产日韩欧美亚洲二区| 嫁个100分男人电影在线观看 | 久久精品熟女亚洲av麻豆精品| 制服人妻中文乱码| 熟女少妇亚洲综合色aaa.| 久久青草综合色| 国产成人精品久久二区二区免费| 美国免费a级毛片| 精品一品国产午夜福利视频| 波多野结衣一区麻豆| 国产成人精品久久二区二区免费| 免费少妇av软件| 欧美日韩成人在线一区二区| 免费在线观看视频国产中文字幕亚洲 | 性色av一级| 90打野战视频偷拍视频| 欧美大码av| 日韩制服丝袜自拍偷拍| 亚洲激情五月婷婷啪啪| 老司机午夜十八禁免费视频| 日本vs欧美在线观看视频| 欧美另类一区| 女人爽到高潮嗷嗷叫在线视频| h视频一区二区三区| 国产一区二区三区av在线| 热re99久久国产66热| 午夜影院在线不卡| 欧美日韩福利视频一区二区| 亚洲精品第二区| videosex国产| 亚洲国产精品一区三区| 日本91视频免费播放| 日韩一本色道免费dvd| 日日爽夜夜爽网站| 热99久久久久精品小说推荐| 18禁裸乳无遮挡动漫免费视频| 亚洲精品成人av观看孕妇| 国产欧美亚洲国产| www.熟女人妻精品国产| 免费黄频网站在线观看国产| 只有这里有精品99| 老司机影院成人| 欧美黑人欧美精品刺激| 少妇的丰满在线观看| 丰满饥渴人妻一区二区三| 国产成人精品在线电影| 久久综合国产亚洲精品| 成年动漫av网址| 69精品国产乱码久久久| 国产精品99久久99久久久不卡| 国产在线视频一区二区| 久久久久精品国产欧美久久久 | 男人添女人高潮全过程视频| 一区二区av电影网| 国产精品一二三区在线看| 国产精品一区二区在线观看99| www.精华液| 狠狠精品人妻久久久久久综合| 51午夜福利影视在线观看| 丰满迷人的少妇在线观看| 亚洲精品一区蜜桃| 国产熟女午夜一区二区三区| 一级黄片播放器| 亚洲成色77777| 国产成人啪精品午夜网站| 中文字幕制服av| 免费观看av网站的网址| 亚洲精品av麻豆狂野| 在线观看免费视频网站a站| 精品国产乱码久久久久久小说| 五月天丁香电影| 看免费成人av毛片| 成人18禁高潮啪啪吃奶动态图| videos熟女内射| 日韩av不卡免费在线播放| 男女高潮啪啪啪动态图| 亚洲欧美激情在线| 男女午夜视频在线观看| 亚洲成av片中文字幕在线观看| 18在线观看网站| 久久热在线av| 婷婷成人精品国产| 在线天堂中文资源库| 中文字幕人妻熟女乱码| 一区二区三区四区激情视频| 少妇粗大呻吟视频| 亚洲精品美女久久av网站| 亚洲人成电影免费在线| 真人做人爱边吃奶动态| 欧美+亚洲+日韩+国产| av有码第一页| 99热网站在线观看| 久久精品亚洲熟妇少妇任你| av欧美777| 亚洲黑人精品在线| 久久久久视频综合| 亚洲国产精品国产精品| av福利片在线| 亚洲国产精品成人久久小说| 新久久久久国产一级毛片| 久久久国产精品麻豆| 亚洲欧洲国产日韩| 亚洲国产欧美网| 国产日韩欧美在线精品| 操美女的视频在线观看| 91字幕亚洲| 十八禁人妻一区二区| 热re99久久国产66热| 亚洲第一av免费看| 美女大奶头黄色视频| 久久久久精品国产欧美久久久 | 欧美av亚洲av综合av国产av| 这个男人来自地球电影免费观看| 久久毛片免费看一区二区三区| 免费女性裸体啪啪无遮挡网站| 丰满人妻熟妇乱又伦精品不卡| 午夜免费成人在线视频| 国产国语露脸激情在线看| 侵犯人妻中文字幕一二三四区| 欧美日本中文国产一区发布| 亚洲情色 制服丝袜| 国产精品久久久av美女十八| 亚洲人成电影免费在线| 免费看av在线观看网站| 欧美亚洲 丝袜 人妻 在线| 2021少妇久久久久久久久久久| 九草在线视频观看| 久久99热这里只频精品6学生| 成人午夜精彩视频在线观看| 日韩中文字幕视频在线看片| 极品少妇高潮喷水抽搐| 国产麻豆69| 丰满饥渴人妻一区二区三| 纵有疾风起免费观看全集完整版| 我要看黄色一级片免费的| 国产一区二区激情短视频 | 午夜激情av网站| 精品国产乱码久久久久久男人| 亚洲色图综合在线观看| 老汉色av国产亚洲站长工具| 欧美日韩视频高清一区二区三区二| 日韩 欧美 亚洲 中文字幕| 亚洲国产av影院在线观看| 丰满饥渴人妻一区二区三| 国产精品久久久久久人妻精品电影 | 欧美成狂野欧美在线观看| 9热在线视频观看99| 国产精品av久久久久免费| 国精品久久久久久国模美| 亚洲成人国产一区在线观看 | 欧美日本中文国产一区发布| 黄色怎么调成土黄色| 午夜免费男女啪啪视频观看| 国产精品欧美亚洲77777| 丁香六月天网| 操美女的视频在线观看| 97在线人人人人妻| 性少妇av在线| 亚洲少妇的诱惑av| av片东京热男人的天堂| 精品人妻一区二区三区麻豆| 国产一区二区在线观看av| 晚上一个人看的免费电影| 一级片免费观看大全| 不卡av一区二区三区| 午夜福利免费观看在线| 中文字幕高清在线视频| 久久鲁丝午夜福利片| 亚洲精品国产色婷婷电影| 国产淫语在线视频| 一边亲一边摸免费视频| 亚洲精品久久成人aⅴ小说| 一级毛片 在线播放| 黑人巨大精品欧美一区二区蜜桃| 午夜老司机福利片| 91字幕亚洲| 亚洲一卡2卡3卡4卡5卡精品中文| 晚上一个人看的免费电影| 两性夫妻黄色片| 欧美日韩精品网址| 国产三级黄色录像| 久久久久视频综合| 日韩一区二区三区影片| 国产成人影院久久av| 97人妻天天添夜夜摸| 久久国产精品大桥未久av| 亚洲精品国产一区二区精华液| 国产视频一区二区在线看| 99精国产麻豆久久婷婷| 国产精品一区二区免费欧美 | 丝袜人妻中文字幕| 麻豆乱淫一区二区| 91国产中文字幕| 日韩一区二区三区影片| 91国产中文字幕| 国产精品国产三级国产专区5o| 色精品久久人妻99蜜桃| 国产亚洲av高清不卡| 一级,二级,三级黄色视频| 欧美日韩av久久| 纵有疾风起免费观看全集完整版| 9热在线视频观看99| 男女下面插进去视频免费观看| 高清欧美精品videossex| 国产av精品麻豆| 久久精品国产综合久久久| 成人亚洲精品一区在线观看| 午夜91福利影院| videosex国产| 一区二区三区乱码不卡18| 波野结衣二区三区在线| 咕卡用的链子| 肉色欧美久久久久久久蜜桃| 久久久久网色| 麻豆av在线久日| 777久久人妻少妇嫩草av网站| 欧美精品一区二区大全| 777久久人妻少妇嫩草av网站| 交换朋友夫妻互换小说| 久久国产精品影院| 一区二区三区精品91| 国产精品国产三级国产专区5o| 丝袜人妻中文字幕| 免费看av在线观看网站| 精品国产乱码久久久久久小说| 中文字幕亚洲精品专区| 视频区图区小说| 久久人人97超碰香蕉20202| 国产国语露脸激情在线看| 久久人人97超碰香蕉20202| 人人妻人人澡人人看| 成人国产一区最新在线观看 | 久久亚洲精品不卡| 久久久久精品国产欧美久久久 | 国产成人精品在线电影| 婷婷色综合www| 欧美黑人欧美精品刺激| 又粗又硬又长又爽又黄的视频| 天堂俺去俺来也www色官网| 国产精品 国内视频| 97精品久久久久久久久久精品| 国产黄色免费在线视频| 999精品在线视频| 99精品久久久久人妻精品| 久久久久久久大尺度免费视频| 在线观看免费午夜福利视频| 亚洲国产精品国产精品| 国产精品人妻久久久影院| 美女中出高潮动态图| 国产亚洲午夜精品一区二区久久| 亚洲第一青青草原| 啦啦啦在线免费观看视频4| 女人久久www免费人成看片| 99久久综合免费| 大陆偷拍与自拍| 高清欧美精品videossex| 男女免费视频国产| 97精品久久久久久久久久精品| 十八禁人妻一区二区| 热re99久久精品国产66热6| 国产精品一国产av| 久久久久久久国产电影| 色婷婷av一区二区三区视频| 美女午夜性视频免费| 多毛熟女@视频| 国产91精品成人一区二区三区 | 久久久精品94久久精品| 国产一卡二卡三卡精品| 2018国产大陆天天弄谢| 国产亚洲av片在线观看秒播厂| 91麻豆av在线| 久久国产精品男人的天堂亚洲| 国产国语露脸激情在线看| 国产一级毛片在线| 老司机午夜十八禁免费视频| 欧美97在线视频| av一本久久久久| 久久天堂一区二区三区四区| av国产精品久久久久影院| 国产免费现黄频在线看| 欧美精品啪啪一区二区三区 | 国产欧美日韩综合在线一区二区| 新久久久久国产一级毛片| xxx大片免费视频| 人妻 亚洲 视频| 久久综合国产亚洲精品| 亚洲成色77777| 亚洲自偷自拍图片 自拍| 纵有疾风起免费观看全集完整版| 亚洲成av片中文字幕在线观看| 黑人猛操日本美女一级片| 亚洲av美国av| 王馨瑶露胸无遮挡在线观看| 亚洲欧美一区二区三区黑人| 亚洲精品第二区| 成人黄色视频免费在线看| 91字幕亚洲| 18禁裸乳无遮挡动漫免费视频| 在线观看免费视频网站a站| 青春草亚洲视频在线观看| www.999成人在线观看| 国产精品一区二区免费欧美 | 成人免费观看视频高清| 久久国产精品人妻蜜桃| www.av在线官网国产| 国语对白做爰xxxⅹ性视频网站| 高清黄色对白视频在线免费看| 亚洲精品成人av观看孕妇| 成人黄色视频免费在线看| 国产在线一区二区三区精| 各种免费的搞黄视频| 一个人免费看片子| 国产97色在线日韩免费| 亚洲精品国产av蜜桃| 欧美日韩一级在线毛片| 亚洲欧美色中文字幕在线| 国产片内射在线| 成年人免费黄色播放视频| 国产免费一区二区三区四区乱码| 高清黄色对白视频在线免费看| 日韩一区二区三区影片| 国产成人一区二区三区免费视频网站 | 久久久久久免费高清国产稀缺| 久久天躁狠狠躁夜夜2o2o | 精品免费久久久久久久清纯 | 国产成人系列免费观看| 久久这里只有精品19| 黄色 视频免费看| 18在线观看网站| 国产黄色视频一区二区在线观看| 久久99一区二区三区| 99国产精品一区二区三区| 久久中文字幕一级| 国产免费一区二区三区四区乱码| 日韩熟女老妇一区二区性免费视频| 国产精品一区二区免费欧美 | 国产在线视频一区二区| 久久久久久久久免费视频了| 午夜福利一区二区在线看| 又粗又硬又长又爽又黄的视频| 波多野结衣av一区二区av| 国产成人91sexporn| 久久精品国产a三级三级三级| 国产成人av教育| 在线亚洲精品国产二区图片欧美| 香蕉国产在线看| 黄网站色视频无遮挡免费观看| 亚洲av国产av综合av卡| 亚洲精品第二区| 亚洲一区二区三区欧美精品| 亚洲欧美日韩另类电影网站| 日韩av在线免费看完整版不卡| 久久人人爽人人片av| 亚洲人成电影观看| 亚洲精品成人av观看孕妇| 亚洲精品自拍成人| 777久久人妻少妇嫩草av网站| 久久性视频一级片| 国产亚洲av片在线观看秒播厂| 久久久久精品国产欧美久久久 | 国产成人精品久久久久久| 建设人人有责人人尽责人人享有的| 一区二区av电影网| 少妇被粗大的猛进出69影院| 亚洲色图综合在线观看| 大型av网站在线播放| 国产成人欧美| 又粗又硬又长又爽又黄的视频| 99国产精品免费福利视频| 亚洲色图综合在线观看| 成人亚洲精品一区在线观看| 9191精品国产免费久久| 满18在线观看网站| 女人被躁到高潮嗷嗷叫费观| 大片电影免费在线观看免费| 啦啦啦 在线观看视频| 老司机影院成人| 后天国语完整版免费观看| 国产精品麻豆人妻色哟哟久久| 搡老乐熟女国产| 成人亚洲欧美一区二区av| 日韩 亚洲 欧美在线| 人人妻人人澡人人看| 精品人妻在线不人妻| 亚洲精品av麻豆狂野| 人成视频在线观看免费观看| 午夜精品国产一区二区电影| 国产精品三级大全| 亚洲成色77777| 97精品久久久久久久久久精品| 久久久精品国产亚洲av高清涩受| 久久毛片免费看一区二区三区| 亚洲av片天天在线观看| 国产成人影院久久av| 亚洲少妇的诱惑av| 精品第一国产精品| 天堂中文最新版在线下载| 亚洲精品中文字幕在线视频| 中文精品一卡2卡3卡4更新| 久久鲁丝午夜福利片| 亚洲午夜精品一区,二区,三区| 欧美日韩一级在线毛片| 久久精品aⅴ一区二区三区四区| 真人做人爱边吃奶动态| 1024视频免费在线观看| 精品人妻一区二区三区麻豆| 两个人看的免费小视频| 亚洲欧美精品自产自拍| 纯流量卡能插随身wifi吗| 亚洲欧洲日产国产| videosex国产| 欧美成人午夜精品| 亚洲精品乱久久久久久| 捣出白浆h1v1| 黄色视频在线播放观看不卡| 亚洲国产av影院在线观看| 18禁国产床啪视频网站| 欧美人与性动交α欧美软件| 久9热在线精品视频| 日韩视频在线欧美| 真人做人爱边吃奶动态| 18在线观看网站| tube8黄色片| 啦啦啦中文免费视频观看日本| 欧美精品一区二区大全| av线在线观看网站| 9色porny在线观看| 国产成人精品无人区| www.自偷自拍.com| 久久午夜综合久久蜜桃| 成年美女黄网站色视频大全免费| 每晚都被弄得嗷嗷叫到高潮| √禁漫天堂资源中文www| 国产高清不卡午夜福利| 欧美日韩黄片免| 国产精品免费大片| 黑丝袜美女国产一区| 妹子高潮喷水视频| 狂野欧美激情性xxxx| 黄色一级大片看看| 一二三四在线观看免费中文在| 桃花免费在线播放| 在线精品无人区一区二区三| 亚洲欧美中文字幕日韩二区| 亚洲伊人久久精品综合| 在线精品无人区一区二区三| 黄色视频在线播放观看不卡| 亚洲av男天堂| 18禁国产床啪视频网站| 亚洲成人手机| 国产女主播在线喷水免费视频网站| 黄网站色视频无遮挡免费观看| 这个男人来自地球电影免费观看| 一级毛片电影观看| 18禁国产床啪视频网站| 日日摸夜夜添夜夜爱| 纵有疾风起免费观看全集完整版| 999精品在线视频| 久久久久视频综合| 高清视频免费观看一区二区| av电影中文网址| 只有这里有精品99| av天堂在线播放| 男女高潮啪啪啪动态图| 国产在线免费精品| 高清欧美精品videossex| 国产黄色免费在线视频| 亚洲 国产 在线| 熟女少妇亚洲综合色aaa.| 久久久久精品人妻al黑|