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

    基于Hilbert-Huang變換的混流泵流動誘導(dǎo)振動試驗

    2018-02-28 05:55:46季磊磊施衛(wèi)東楊勇飛平元峰張文全
    農(nóng)業(yè)工程學(xué)報 2018年2期
    關(guān)鍵詞:方向振動信號

    李 偉,季磊磊,施衛(wèi)東,楊勇飛,平元峰,張文全

    ?

    基于Hilbert-Huang變換的混流泵流動誘導(dǎo)振動試驗

    李 偉,季磊磊,施衛(wèi)東,楊勇飛,平元峰,張文全

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

    混流泵水力誘導(dǎo)的機(jī)組振動是混流泵運(yùn)行失穩(wěn)的重要因素之一,為了研究混流泵水力激振誘導(dǎo)的機(jī)組振動情況,基于本特利408數(shù)據(jù)采集系統(tǒng),測量獲得了空載和負(fù)載工況下混流泵泵體和泵體基座不同位置處的振動信號,通過希爾伯特-黃變換對原始振動信號進(jìn)行經(jīng)驗篩分分解,獲得了不同模函數(shù)分量的頻譜分布。研究結(jié)果表明,相比空載運(yùn)行,混流泵負(fù)載工況運(yùn)行時水力誘導(dǎo)的機(jī)組振動明顯加劇,但在不同方向上,水力激振引起的振動各不相同。方向上2個工況下的振動頻譜分布基本相似,而在方向、方向和混流泵底座上,負(fù)載工況下波形的頻帶分布變窄,能量分布較為集中,且方向的原始振幅要明顯大于方向,約為方向原始振幅的2倍。混流泵負(fù)載工況運(yùn)行時,低頻振動占據(jù)主要振動能量分布,使得不同模函數(shù)分量的主頻向低頻方向移動,水力誘導(dǎo)混流泵機(jī)組的振動以中低頻振動為主。該研究可為有效降低或防止混流泵水力誘導(dǎo)的機(jī)組振動惡化提供參考。

    泵;振動;信號分析;混流泵;空載;希爾伯特-黃變換

    0 引 言

    混流泵廣泛應(yīng)用于農(nóng)業(yè)排灌、城市供排水、礦山、大型水利工程等領(lǐng)域。在混流泵機(jī)組的運(yùn)行過程中,總是存在著不同程度的振動,一方面是由于軸向平衡共振、轉(zhuǎn)子不平衡、安裝原因?qū)е碌臋C(jī)械振動,另一方面很多已經(jīng)確定的振動問題都與流動的大尺度振蕩有關(guān)。當(dāng)這些振動超過一定限度時,就會對設(shè)備造成不同程度的危害,不僅降低水泵效率,而且縮短了機(jī)組零部件的使用壽命,嚴(yán)重時甚至導(dǎo)致機(jī)組被迫停機(jī)[1-5]。

    隨著水力機(jī)械單機(jī)容量和尺寸的增加,人們對流動誘導(dǎo)機(jī)組振動的研究越來越重視,許多文獻(xiàn)探討了非定常流動誘導(dǎo)的水力機(jī)械系統(tǒng)穩(wěn)定性問題[6-9]。在水輪機(jī)領(lǐng)域,張松松[10]研究了不同工況下額定功率為500 kW的小型混流式水輪發(fā)電機(jī)組的振動問題,研究發(fā)現(xiàn)在負(fù)荷變動的過程中,發(fā)電機(jī)內(nèi)部磁拉力不平衡,機(jī)組的軸線不直和對中不良,以及導(dǎo)葉開度變化造成的水力不平衡是導(dǎo)致軸心軌跡出現(xiàn)大范圍偏移的主要原因。Shi等[11]采用數(shù)值模擬的方法研究了非定常流動下導(dǎo)葉式水輪機(jī)內(nèi)部的壓力脈動特性,發(fā)現(xiàn)渦室內(nèi)壓力脈動的時域特性呈明顯的周期性變化;壓力脈動的主頻集中在低頻區(qū)域,并且是葉片通過頻率的2倍。孟龍等[12]發(fā)現(xiàn)機(jī)組同時存在轉(zhuǎn)子質(zhì)量不平衡問題及間隙過大問題,根據(jù)軸心軌跡特性逐步調(diào)整上導(dǎo)及水導(dǎo)間隙。

    然而,混流泵作為水力機(jī)械的一種,由于水力不穩(wěn)定誘導(dǎo)機(jī)組振動的研究文獻(xiàn)較少,大部分學(xué)者僅單純的關(guān)注了泵內(nèi)非定常流場結(jié)構(gòu)或僅單一的研究了泵不穩(wěn)定運(yùn)行時的振動特征[13-15]。施衛(wèi)東等[16]研究了高比轉(zhuǎn)速混流泵不同工況下的葉輪進(jìn)出口等位置的壓力脈動時域和頻域特性,研究結(jié)果表明混流泵最大壓力脈動發(fā)生在葉輪進(jìn)口前,泵運(yùn)行偏離最優(yōu)工況越遠(yuǎn),葉輪進(jìn)口處壓力系數(shù)幅值越大。劉建瑞等[17]研究發(fā)現(xiàn)在相同流量下,徑向間隙越小,監(jiān)測點頻域振幅越大。李偉等[18]發(fā)現(xiàn)流量工況的改變影響了轉(zhuǎn)子系統(tǒng)的不平衡量和不對中程度,并且隨著遠(yuǎn)離設(shè)計流量點,不平衡量引起的軸系工頻振動和不對中引起的水平方向振動不斷加劇。Wang等[19]研究發(fā)現(xiàn),隨著流量的增大,壓力脈動的峰值逐漸減小,不同監(jiān)測點處壓力脈動的主頻為葉片通過頻率和其倍頻,而不同監(jiān)測點處的振動頻率為軸頻和其倍頻。

    在振動信號的處理方面,希爾伯特-黃變換(Hilbert-Huang transform,HHT)是在傅立葉變換、小波變換等方法的基礎(chǔ)上構(gòu)建的一種信號時頻分析理論,它引入了固有模態(tài)和經(jīng)驗篩分(empirical mode decomposition,EMD)概念,能夠?qū)⒄駝有盘柗纸鉃橛邢薜木哂袑嶋H物理意義的本征模函數(shù),并對每個本征模函數(shù)進(jìn)行Hilbert變換,得到每個本征模函數(shù)的瞬時頻譜,常用來對機(jī)械故障進(jìn)行診斷[20-22]。本文以導(dǎo)葉式混流泵為研究對像,在前期研究其內(nèi)部流動和軸系振動的基礎(chǔ)上[23-24],基于本特利408數(shù)據(jù)采集系統(tǒng)和希爾伯特-黃變換,對比分析了混流泵空載和負(fù)載運(yùn)行時的振動特性,探索混流泵水力非定常特性誘導(dǎo)的機(jī)組振動,為降低或防止混流泵運(yùn)行中的機(jī)組振動惡化提供參考。

    1 試驗對象及試驗裝置

    1.1 試驗對象

    本文研究的導(dǎo)葉式混流泵模型的參數(shù)如下:流量opt=380 m3/h,揚(yáng)程=6m,轉(zhuǎn)速=1 450 r/min,比轉(zhuǎn)速s=480。葉片數(shù)=4,導(dǎo)葉葉片數(shù)d=7。葉輪進(jìn)口直徑in=92 mm,出口直徑out=121 mm試驗用混流泵模型泵如圖1所示。

    圖1 混流泵模型

    1.2 泵外特性試驗裝置

    試驗測量在江蘇大學(xué)國家水泵及系統(tǒng)工程技術(shù)研究中心的250 mm不銹鋼軸(混)流泵閉式試驗臺上進(jìn)行,該試驗臺專門用于混流泵、軸流泵模型泵段及模型裝置試驗,試驗裝置系統(tǒng)如圖2所示。

    1.增壓泵 2.渦輪流量計 3.出口閘閥 4.水箱 5.進(jìn)口閘閥 6.伸縮管 7.排氣孔 8.進(jìn)口測壓段 9.試驗泵段 10.扭矩儀 11.電機(jī) 12.排氣孔 13.出口測壓段 14.橡膠軟接頭

    在試驗泵段進(jìn)、出口位置各安置一個麥克公司生產(chǎn)壓力變送器以便進(jìn)行揚(yáng)程的測量,壓力變送器為WT-1151型電容式壓力變送器;進(jìn)口測量用壓力變送器測量范圍為±100 kPa,精度為0.2;出口測量用壓力變送器測量范圍為0~600 kPa,精度為0.2。流量測量采用上海自儀九生產(chǎn)的公稱壓力為1.6 MPa,精度為0.5的LWGY-250型渦輪流量計。采用上海良標(biāo)智能終端股份有限公司生產(chǎn)的精度為0.2級的ZJ型轉(zhuǎn)矩轉(zhuǎn)速測量儀測量模型泵的軸轉(zhuǎn)速、轉(zhuǎn)矩和軸功率,轉(zhuǎn)矩轉(zhuǎn)速測量儀測量誤差為±0.2%。測試系統(tǒng)達(dá)到1級精度要求。

    2 振動數(shù)據(jù)采集系統(tǒng)與信號分析方法

    2.1 振動數(shù)據(jù)采集系統(tǒng)

    按照試驗要求搭建好實驗臺并連接好數(shù)據(jù)采集和檢測系統(tǒng)。將本特利加速度傳感器接入本特利408數(shù)據(jù)采集系統(tǒng)。試驗采用本特利公司生產(chǎn)的200350加速度傳感器進(jìn)行振動信號的采集,該傳感器靈敏度系數(shù)為100,頻率范圍0.5~10 kHz,傳感器頭部有電磁鐵,使用時直接將傳感器貼于待測位置即可。為了獲得混流泵泵體在不同方向上的振動特性,在泵體環(huán)形渦室的對稱中心處、靠近轉(zhuǎn)軸處和泵體基座上設(shè)置了傳感器進(jìn)行監(jiān)測,傳感器的現(xiàn)場布置如圖3所示。傳感器輸出的信號由本特利408數(shù)據(jù)采集系統(tǒng)進(jìn)行采集和后處理,該系統(tǒng)由408動態(tài)信號處理儀器(DSPi)以及ADRE Sxp軟件構(gòu)成。其中,ADRE Sxp軟件可進(jìn)行各種信號的幅值分析、頻譜分析、變化趨勢分析等,還可通過自選擇濾波繪制伯德圖和頻譜圖以及軸心軌跡等各種試驗結(jié)果分析曲線。

    1.模型泵 2.電機(jī) 3. Y方向振動傳感器(CH2) 4.底座上振動傳感器(CH4) 5. X方向振動傳感器(CH1) 6. Z方向振動傳感器(CH3) 7.本特利408數(shù)據(jù)采集系統(tǒng) 8.壓力變送器

    首先,進(jìn)行額定負(fù)載工況下的試驗。在試驗開始階段,將出口閥門全開,先啟動泵類產(chǎn)品測試系統(tǒng)并調(diào)試好軟件,開始記錄數(shù)據(jù),再啟動電機(jī),調(diào)節(jié)出口閥門,使流量計讀數(shù)達(dá)到設(shè)計工況點并且穩(wěn)定在380 m3/h,轉(zhuǎn)速穩(wěn)定在1 450 r/min時,分別獲取能量性能參數(shù),隨后啟動并觸發(fā)本特利408數(shù)據(jù)采集系統(tǒng),開始記錄泵體的振動信號。關(guān)閉電源停機(jī),待管路內(nèi)流體重新穩(wěn)定,重新將出口閥門全開,進(jìn)行上述3次重復(fù)性試驗,并記錄試驗數(shù)據(jù)。在相同運(yùn)行條件下進(jìn)行空載試驗,關(guān)閉混流泵上下游管道閥門,并將泵體內(nèi)的水排盡,振動信號采集方法與負(fù)載試驗時相同。

    2.2 振動信號的分解和變換

    HHT變換是基于信號局部特征,自適應(yīng)地篩選出模函數(shù)分量,克服了小波變換中選擇小波基的困難[25]。同時,它可以在時間和頻率同時達(dá)到很高的精度,克服了傳統(tǒng)的傅里葉變換只能得到信號某一段時間內(nèi)頻率的均值,無法準(zhǔn)確描述頻率-時間變化的缺陷。HHT變換通過生成復(fù)解析信號,得到復(fù)平面上具有明確解析意義的瞬時頻率[26-28]。故本文采用Hilbert-Huang變換對獲得的振動信號進(jìn)行處理。

    EMD是用波動上、下包絡(luò)平均值確定“瞬時平衡位置”的方法,分解過程相當(dāng)于對原始信號進(jìn)行分頻濾波,最終得到不同頻段的特征信號。對于任意時間序列()經(jīng)過EMD分解,可以得到一系列模函數(shù)分量,其表達(dá)式如式(1)所示。

    式中本征模式函數(shù)1()、2()、···、c()分別表示原始信號中所有頻率下由高頻到低頻的信號成分;余量r()表示信號的趨勢項或漂移。

    在任意的某一時間序列()時,Hilbert-Huang變換()可表示為

    構(gòu)造解析函數(shù)

    式中()和()分別稱為信號()的瞬時振幅和瞬時相位。其中

    由瞬時相位可得信號的瞬時頻率()

    在進(jìn)行Hilbert-Huang變換之后,各模函數(shù)分量可以轉(zhuǎn)化為信號對應(yīng)的以時間為變量的瞬時振幅和瞬時頻率。由于Hilbert-Huang變換對局部特性的要求較高,所以應(yīng)去除所分析信號在頻域中的虛假成分。以時間和頻率為自變量,振幅為因變量,就能得到Hilbert-Huang幅值譜(,),即

    式中表示取實部;a(),ω()是以時間為變量函數(shù),可以構(gòu)成時間、頻率、幅值的三維時頻譜圖。

    3 試驗結(jié)果與分析

    3.1 能量性能試驗結(jié)果

    獲得混流泵模型3組試驗外特性數(shù)據(jù),如圖4所示。通過試驗結(jié)果可知,所進(jìn)行的3次試驗所得的揚(yáng)程和效率曲線基本一致,揚(yáng)程的最大誤差在3%以內(nèi)而效率誤差在2%以內(nèi),證明了本次試驗可重復(fù)性好,試驗所得結(jié)果可靠性高。

    3.2 機(jī)組振動信號的EMD分解

    測試得到了混流泵在設(shè)計流量工況和空載運(yùn)行下泵體的(CH1)、(CH2)、(CH3)3個方向上的振動信號以及泵體底座上的振動信號(CH4),通過對4個加速度傳感器所得信號進(jìn)行希爾伯特-黃變換,得到了不同傳感器所對應(yīng)的EMD分解圖(未全部標(biāo)出),如圖5所示。

    注:Q為試驗流量工況,m3·h-1;Qopt為設(shè)計流量工況,m3·h-1。

    由圖5可知,從原始信號可以看出,在負(fù)載工況下,各個傳感器所采集的信號的幅值均有所增加,尤其是泵體上方向、方向和底座上的振動明顯增大,而在方向上,振動增加不明顯,從2個流量工況下方向傳感器對應(yīng)的EMD分解圖可以看出,和空載運(yùn)行時相似,混流泵在設(shè)計流量工況下運(yùn)行時,方向上各個模函數(shù)分量的曲線較為相似,說明混流泵在運(yùn)行時,水力激振對方向上的振動影響較小,這和混流泵的安裝形式和出口方向有關(guān)。從2個流量工況下方向和方向傳感器對應(yīng)的原始振動信號和其EMD分解圖可以看出,混流泵在空載運(yùn)行時,其泵體在2個方向上的振動幅值較小,而當(dāng)混流泵運(yùn)行在負(fù)載工況時,2個方向上的振幅均明顯增大,呈條帶狀分布,并且方向的原始振幅要明顯大于方向,約為方向原始振幅的2倍。同時,方向EMD分解圖上對應(yīng)的各模函數(shù)分量和剩余項的振幅也均要大于相應(yīng)的方向上的各分量幅值。

    在泵體底座上,雖然傳感器所測方向和方向相同,但其原始振動的幅值卻略小于方向的振幅,這是由于泵體底座和試驗臺的接觸面積較大,泵體上振動傳遞到底座上有所減弱,因而振幅相對較小。由于在負(fù)載工況下方向上的振幅明顯增大,因此,對比方向上空載工況和負(fù)載工況下的各模函數(shù)分量和剩余項的振幅可知,隨時間的變化,負(fù)載工況下各模函數(shù)分量的振幅在一段時間內(nèi)存在周期性,各個波峰值和波谷值出現(xiàn)的頻次較高且相鄰波峰值和波谷值的差值相差不大,而空載工況下各模函數(shù)分量的振幅在不同時間段內(nèi)會出現(xiàn)極值,并且波峰值和波谷值相差較大。因此,在負(fù)載工況下,混流泵泵體的整體振動幅值明顯增加,而在空載運(yùn)行下,雖然整泵的振幅較小,但會出現(xiàn)振動極值。綜上所述,混流泵在負(fù)載運(yùn)行時,其水力誘導(dǎo)的機(jī)組振動明顯加劇,EMD分解圖上各模函數(shù)分量上和剩余項的振幅也明顯增加,但在不同方向上,水力激振引起的振動幅值各不相同。

    圖5 振動信號EMD分解圖

    3.3 設(shè)計工況和空載工況下頻譜分析

    獲得了不同模函數(shù)分量的頻譜分布(未全部標(biāo)出),如圖6所示。混流泵空載和額定工況載荷下振動原始信號頻譜較為復(fù)雜,除主頻信號及前幾階分頻信號外,其余信號大多被寬頻信號所覆蓋,難以看出其頻譜信號特征。經(jīng)EMD分解后,各階模函數(shù)的頻譜分布由模函數(shù)1到殘余項11,主頻呈逐漸降低趨勢。結(jié)合圖5可知EMD成功的將振動信號按頻率高低按降序分解出來,其中殘余項11與模函數(shù)11頻率分布主要集中于零附近,其所對應(yīng)模函數(shù)信號周期性也不明顯,可判定此兩項為分解中產(chǎn)生的虛假分量[29],對于混流泵振動信號分析不具有參考意義。模函數(shù)9附近模量的頻譜主要成分為軸頻(24.16 Hz)以下的低頻信號,這部分信號主要對應(yīng)于葉輪與導(dǎo)葉內(nèi)部液體回流以及漩渦結(jié)構(gòu)等低頻流動現(xiàn)象。模函數(shù)5~模函數(shù)7模量的頻譜能量較高且分布較為集中,對應(yīng)信號主要來源于轉(zhuǎn)子的轉(zhuǎn)動以及葉輪與導(dǎo)葉葉片之間的動靜干涉。模函數(shù)1~模函數(shù)3模量的頻譜主要反應(yīng)機(jī)械的高頻信號,根據(jù)文獻(xiàn)[30]中結(jié)論,這部分信號對機(jī)組的不穩(wěn)定以及故障較為敏感。

    混流泵在空載時,在不同位置處的頻譜分布基本相似,原始信號頻譜分布的幅值約為0。而混流泵在負(fù)載工況時,各個傳感器的頻譜分布各不相同。在方向上,2個工況下混流泵的振動頻譜分布基本相似,不同模函數(shù)分量的頻譜分布也基本一致,說明水力激振對方向上的振動影響較小。相比混流泵空載時的頻譜分布,在負(fù)載工況下,方向上原始信號對應(yīng)的頻譜幅值波動較大,頻譜分布較為雜亂,分頻成分較多,隨著模函數(shù)分量的增大,各分量下波形的主頻從高頻向低頻移動。同時,對比2個工況相對應(yīng)的模函數(shù)分量可知,在不同模函數(shù)分量下,負(fù)載工況波形的頻域分布范圍略有縮減。對比方向上空載工況和負(fù)載工況下的各模函數(shù)分量和剩余項的頻譜分布可知,在模函數(shù)1~模函數(shù)3分量下,各個模函數(shù)分量的振動頻譜出現(xiàn)明顯區(qū)別,在空載工況下,3個分量下的振動幅值在600 Hz附近有最高值,而在負(fù)載工況下,這3個模函數(shù)分量下的振動極值從900向300 Hz方向偏移,在其他模函數(shù)分量下,振動頻譜的極值均出現(xiàn)在200 Hz內(nèi),說明混流泵在負(fù)載工況下,水力誘導(dǎo)振動以中、低頻振動為主。在方向上的振動頻譜分布和方向趨勢類似,均是隨著模函數(shù)分量的增大,頻域波形的主頻從高頻向低頻移動。同時,相比空載工況,在相同分量下,波形的頻帶分布變窄,能量分布較為集中,說明水力誘導(dǎo)的機(jī)組振動加劇了泵體的振動,并且水力振動使得振動主頻向低頻移動,這是水力誘導(dǎo)泵體振動最直觀的體現(xiàn)。在原始振動波形的頻譜分布中,在200~600 Hz內(nèi)振動頻譜的幅值較高,說明該頻段內(nèi),振動的能量分布較多,水力誘導(dǎo)泵體振動較強(qiáng)。在泵體底座上,從原始信號的頻譜分布可以看出,在0~400 Hz內(nèi),負(fù)載工況下原始波形頻譜分布的幅值較高,約為其他頻率下幅值的2倍,說明水力因素導(dǎo)致的混流泵底座的振動基本是低頻振動。綜上所述,混流泵在負(fù)載運(yùn)行時,水力誘導(dǎo)振動對泵體的影響較為明顯,低頻振動占據(jù)主要振動能量分布,使得不同模函數(shù)分量的主頻向低頻方向移動,水力誘導(dǎo)混流泵機(jī)組的振動以中、低頻振動為主。

    圖6 不同模函數(shù)分量的頻譜分布

    4 結(jié) 論

    1)通過對原始振動信號進(jìn)行希爾伯特-黃變換可知,混流泵在負(fù)載運(yùn)行時,水力誘導(dǎo)的機(jī)組振動明顯加劇,但在不同方向上,水力激振引起的振動各不相同,方向、方向和底座上的振動明顯增大,而在方向上,振動增加不明顯。

    2)混流泵空載運(yùn)行時,在不同位置處的頻譜分布基本相似,原始信號頻譜分布的幅值約為0。而在負(fù)載工況運(yùn)行時,方向上2個工況下的振動頻譜分布基本相似,在方向、方向和混流泵底座上,波形的頻帶分布變窄,能量分布較為集中,水力誘導(dǎo)的機(jī)組振動加劇了泵體的振動。

    3)混流泵在負(fù)載運(yùn)行時,低頻振動占據(jù)主要振動能量分布,使得不同模函數(shù)分量的主頻向低頻方向移動。水力誘導(dǎo)混流泵機(jī)組的振動以中低頻振動為主。因此,減弱混流泵負(fù)載運(yùn)行時的中低頻振動是提高混流泵運(yùn)行穩(wěn)定性的關(guān)鍵。

    [1] 呂海平. 臥式離心泵受力與軸的振動對泵體破壞分析[J]. 才智,2012,2012(17):49-50.

    Lü Haiping. Failure analysis of horizontal centrifugal pump force and shaft vibration on pump body[J].Intelligence, 2012, 2012(17): 49-50. (in Chinese with English abstract)

    [2] Christopher S, Kumaraswamy S. Study of noise and vibration signal for a radial flow pump during performance test[J]. Fluid Mechanics and Fluid Power-Contemporary Research, 2017,12: 853-861.

    [3] Li W, Shi W, Xu Y, et al. Effects of guide vane thickness on pressure pulsation of mixed-flow pump in pumped-storage power station[J]. Journal of Vibroengineering, 2013, 15(3): 1177-1185.

    [4] 李偉, 季磊磊, 施衛(wèi)東,等.不同流量工況下混流泵壓力脈動試驗[J]. 農(nóng)業(yè)機(jī)械學(xué)報, 2016, 47(12):70-76.

    Li Wei, Ji Leilei, Shi Weidong, et al. Experiment on pressure fluctuation in mixed-flow pump under different flow rate conditions[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(12): 70-76. (in Chinese with English abstract)

    [5] Li W, Ji L, Shi W, et al. Vibration characteristics of the impeller at multi-conditions in mixed-flow pump under the action of fluid-structure interaction[J]. Journal of Vibroengineering, 2016, 18(5): 3213-3224 .

    [6] 潘羅平,高明. 水輪機(jī)水力穩(wěn)定性的分析[J]. 長春工程學(xué)院學(xué)報:自然科學(xué)版,2002,3(4):41-43.

    Pan Luoping, Gao Ming. Analysis of hydraulic stability of hydraulic turbine[J]. Journal of Changchun Institute of Technology: Natural Science Edition, 2002, 3(4): 41-43. (in Chinese with English abstract)

    [7] 陶星明,劉光寧. 高比速混流式水輪機(jī)的水力穩(wěn)定性問題[J]. 大電機(jī)技術(shù),2003(4):46-48.

    Tao Xingming, Liu Guangning. Hydraulic stability problem of high specific speed mixed flow turbine[J]. Large motor technology, 2003(4): 46-48. (in Chinese with English abstract)

    [8] 劉麗娜,潘偉峰,麻志成. 混流式水輪機(jī)低負(fù)荷運(yùn)行水力穩(wěn)定性的研究[J]. 水電自動化與大壩監(jiān)測,2015(2):20-23. Liu Lina, Pan Weifeng, Ma Zhicheng. Study on hydraulic stability of mixed flow turbine at low load[J]. Hydropower Automation and Dam Monitoring, 2015(2): 20-23. (in Chinese with English abstract)

    [9] 成立,吳璐璐,劉超. 大型軸流泵水力不穩(wěn)定區(qū)研究[J]. 灌溉排水學(xué)報,2010,29(2):102-104.

    Cheng Li, Wu Lulu, Liu Chao. Study on hydraulic instability zone of large axial flow pump[J]. Journal of Irrigation and Drainage, 2010, 29(2): 102-104. (in Chinese with English abstract)

    [10] 張松松. 小型混流式水輪發(fā)電機(jī)組振動試驗與分析[D]. 邯鄲:河北工程大學(xué),2014.

    Zhong Songsong. Vibration Test and Analysis of Small Mixed Flow Turbine[D]. Handan:Hebei University of Engineering, 2014. (in Chinese with English abstract)

    [11] Shi F X, Yang J H, Wang X H. Analysis on characteristic of pressure fluctuation in hydraulic turbine with guide vane[J]. International Journal of Fluid Machinery and Systems, 2016, 9(3): 237-244.

    [12] 孟龍,劉孟,支發(fā)林,等. 機(jī)械不平衡及軸瓦間隙對水輪機(jī)運(yùn)行穩(wěn)定性的影響分析[J]. 機(jī)械工程學(xué)報,2016,52(3):49-55.

    Meng Long, Liu Meng, Zhi Faling, et al. Analysis of the effect of mechanical imbalances and bush clearance on the stability of hydraulic turbines[J]. Journal of Mechanical Engineering, 2016, 52(3): 49-55. (in Chinese with English abstract)

    [13] 李偉,季磊磊,施衛(wèi)東,等. 混流泵非均勻輪緣間隙流場數(shù)值計算[J]. 農(nóng)業(yè)機(jī)械學(xué)報,2016,47(10):66-72.

    Li Wei, Ji Leilei, Shi Weidong, et al. Numerical calculation of the flow field in the nonuniform rim clearance of mixed flow pump[J]. Transactions of the Chinese Society for Agricultural Machinery 2016, 47(10): 66-72. (in Chinese with English abstract)

    [14] 李偉. 導(dǎo)葉式混流泵多工況內(nèi)部流場的PIV測量[J]. 農(nóng)業(yè)工程學(xué)報,2016,32(24):82-88.

    Li Wei. PIV measurement of the internal flow field in multiple conditions of guide vane mixed flow pump[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(24): 82-88. (in Chinese with English abstract)

    [15] 季磊磊,李偉,施衛(wèi)東,等. 導(dǎo)葉式混流泵內(nèi)部非定常流動特性數(shù)值模擬[J]. 農(nóng)業(yè)機(jī)械學(xué)報,2016,47(增刊1):155-162,188.

    Ji Leilei, Li Wei, Shi Weidong, et al. Numerical simulation of unsteady flow in guide vane mixed flow pump[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016,47(Supp.1):155-162,188. (in Chinese with English abstract)

    [16] 施衛(wèi)東,鄒萍萍,張德勝,等. 高比轉(zhuǎn)速斜流泵內(nèi)部非定常壓力脈動特性[J]. 農(nóng)業(yè)工程學(xué)報,2011,27(4):147-152. Shi Weidong, Zou Pingping, Zhang Deisheng, et al. Unsteady pressure pulsation characteristics in high specific speed slanted flow pump[J].Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011, 27(4): 147-152. (in Chinese with English abstract)

    [17] 劉建瑞,鄭俊峰,付登鵬,等. 混流泵徑向間隙對內(nèi)部非定常流場影響的分析[J]. 流體機(jī)械,2014(3):19-23.

    Liu Jianrui, Zhen Junfemg, Fu Dengpeng, et al. Analysis of the influence of radial clearance on unsteady flow field in mixed flow pump[J]. Fluid machinery, 2014(3): 19-23. (in Chinese with English abstract)

    [18] 李偉,季磊磊,施衛(wèi)東,等. 變流量工況對混流泵轉(zhuǎn)子軸心軌跡的影響[J]. 農(nóng)業(yè)工程學(xué)報,2016,32(4):91-97.

    Li Wei, Ji Leilei, Shi Weidong, et al. Effects of variable flow conditions on rotor axis orbit of mixed flow pump[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(4): 91-97. (in Chinese with English abstract)

    [19] Wang K, Liu H, Zhou X, et al. Experimental research on pressure fluctuation and vibration in a mixed flow pump[J]. Journal of Mechanical Science and Technology, 2016, 30(1): 179-184.

    [20] 薛延剛,羅興錡,王瀚,等. 基于改進(jìn)HHT方法提取水輪機(jī)動態(tài)特征信息研究[J]. 水力發(fā)電學(xué)報,2011,30(4):214-221.

    Xue Yangang, Luo Xingqi, Wang Han, et al. Study on dynamic feature extraction of hydraulic turbine based on improved HHT method[J]. Journal of Hydroelectric Power Generation, 2011, 30(4): 214-221. (in Chinese with English abstract)

    [21] 梁興,劉梅清,劉志勇,等. 立式混流泵異常振動測試分析[J]. 排灌機(jī)械工程學(xué)報,2013,31(5):373-378.

    Liao Xing, Liu Meiqing, Liu Zhiyong, et al. Test and analysis of abnormal vibration of vertical mixed flow pump[J]. Drainage and Irrigation Machinery Engineering, 2013, 31(5): 373-378. (in Chinese with English abstract)

    [22] 楊國安,王澤棟. 基于改進(jìn)希爾伯特-黃的泵閥故障診斷新方法[J]. 北京化工大學(xué)學(xué)報:自然科學(xué)版,2008,35(4):81-85.

    Yang Guoan, Wang Zeidong. A new method for fault diagnosis of pump valve based on improved Hilbert-Huang[J]. Journal of Beijing University of Chemical Technology: Natural Science Edition, 2008, 35(4): 81-85. (in Chinese with English abstract)

    [23] Li W, Zhou L, Shi W, et al. PIV experiment of the unsteady flow field in mixed-flow pump under part loading condition[J]. Experimental Thermal and Fluid Science, 2017, 83(2017): 191-199.

    [24] 李偉,季磊磊,施衛(wèi)東,等. 混流泵起動過程轉(zhuǎn)子軸心軌跡的試驗研究[J]. 機(jī)械工程學(xué)報,2016, 52(22):168-177.

    Li Wei, Ji Leilei, Shi Weidong, et al. Experimental study on the rotor axis orbit in the starting process of mixed flow pump[J]. Journal of Mechanical Engineering, 2016, 52(22): 168-177. (in Chinese with English abstract)

    [25] Zhu K P,Wong Y S,Hong G S. Wavelet analysis of sensor signals for tool condition monitoring: A review and some new results[J]. International Journal of Machine Tools & Manufacture, 2009, 49(2009): 537-553.

    [26] Zhang Y, Tang B, Xiao X. Time-frequency interpretation of multi-frequency signal from rotating machinery using an improved Hilbert–Huang transform[J]. Measurement, 2016, 82(2016): 221-239.

    [27] Interpretation of mechanical signals using an improved Hilbert–Huang transform[J]. Mechanical Systems & Signal Processing, 2008, 22(5): 1061-1071.

    [28] Babu T R, Srikanth S, Sekhar A S. Hilbert–Huang transform for detection and monitoring of crack in a transient rotor[J]. Mechanical Systems & Signal Processing, 2008, 22(4): 905-914.

    [29] 劉美汝. 基于HHT的主泵飛輪振動監(jiān)測系統(tǒng)研究[D]. 哈爾濱:哈爾濱工程大學(xué),2015.

    Liu Meiru. Research on the Vibration Monitoring System of the Main Pump Flywheel based on HHT[D]. Harbin:Harbin Engineering University, 2015. (in Chinese with English abstract)

    [30] 王若. 火箭發(fā)動機(jī)渦輪泵實時故障檢測方法研究[D]. 成都:電子科技大學(xué),2013.

    Wang Ruo. Research on Real Time Fault Detection Method for Turbopump of Rocket Engine[D]. Chengdu: University of Electronics Technology, 2013. (in Chinese with English abstract)

    Experiment of flow induced vibration of mixed-flow pump based on Hilbert-Huang transform

    Li Wei, Ji Leilei, Shi Weidong, Yang Yongfei, Ping Yuanfeng, Zhang Wenquan

    (,,212013,)

    The vibration induced by flow is one of the important factors to the instability of the mixed-flow pump. With the increase of the capacity and the size of mixed-flow pumps, the vibration of the hydraulic components attracts more attention from researchers and engineers. In order to study the vibration of mixed-flow pump induced by hydrodynamic force, based on the Bentley 408 data acquisition system, vibration signals in 3 directions (,and) on the base under unloaded and designed flow conditions of the mixed-flow pump are tested and then analyzed using Hilbert-Huang transform (HHT). The original vibration signal is decomposed by EMD (empirical mode decomposition) using the Hilbert-Huang transform, and the spectral distribution of the different mode function components is obtained. The decomposed signal contains intrinsic modulus with 11 different orders and one residual. The acquired intrinsic modulus represents vibration signal with different frequencies, except the Intrinsic mode function 11 and the residual, which show no periodic characteristics. The results show that the vibration under loaded condition is significantly increased compared with the unloaded operating condition, but the vibrations caused by the hydraulic excitation differ in different directions. When the mixed-flow pump is working under designed operating condition, the low frequency vibration occupies the main vibration energy distribution, making the main mode of the different mode function components move to the low frequency direction, and the hydraulic induced vibration is dominated by the middle and low frequency vibration. In thedirection, the vibration spectrum distribution is similar in the 2 cases, while in theanddirection and on the base of mixed-flow pump the frequency distribution of the waveform is narrowed and the energy distribution is concentrated under loaded condition. Under both loaded and unloaded conditions, the main frequency of the spectrum moves toward the low frequency region with the increase of the modulus order. The difference of the frequency spectrum under loaded and unloaded conditions is more obvious for the intrinsic modulus with lower order, namely Intrinsic mode function 1-3. Under unloaded condition, frequency spectrum of the Intrinsic mode function 1-3 reaches the highest value near 600 Hz; while under loaded operating condition, the position for the peak values moves from 900 to 300 Hz for the spectrum of the intrinsic modulus of the first 3 orders. For the other intrinsic moduli, the peak of the frequency spectrum is distributed below 200 Hz, which indicates that the vibration induced by flow is mainly composed of low and middle frequency vibrations. Compared with the frequency spectrum of vibration signal under unloaded condition, the distribution of vibration spectrum under loaded operating condition is narrower, and the power distribution of the vibration is more concentrated, indicating that the vibration induced by flow increases the vibration of the pump and the hydraulic vibration makes the main frequency of the pump vibration move toward the low frequency region in the frequency spectrum. In the original spectrum under loaded operating condition, the amplitude of the vibration spectrum is higher in the region of 200-600 Hz, which indicates that the vibration power is high in this frequency section, and the vibration caused by hydraulic force is greater there. According to the original vibration signal on the base of the pump, it was found that in the frequency region from 0 to 400 Hz, the amplitude of the vibration spectrum under loaded operating condition is 2 times higher than that of other frequencies, which indicates that low frequency vibration occupies the vibration on the base of the mixed-flow pump. The research results have important engineering application value and theoretical guidance for effectively reducing or preventing the vibration of mixed-flow pumps.

    pumps; vibrations; signal analysis; mixed-flow pump; unloaded condition; Hilbert-Huang transform

    10.11975/j.issn.1002-6819.2018.02.007

    TH313

    A

    1002-6819(2018)-02-0047-08

    2017-07-02

    2017-12-08

    國家自然科學(xué)基金項目(51679111、51579118);江蘇省自然科學(xué)基金項目(BK20161472);江蘇高校優(yōu)勢學(xué)科建設(shè)工程資助項目(PAPD)

    李 偉,博士,研究員,博士生導(dǎo)師,主要研究方向為流體機(jī)械(泵)的優(yōu)化設(shè)計、流場計算和動力學(xué)特性研究。Email:lwjiangda@ujs.edu.cn

    李 偉,季磊磊,施衛(wèi)東,楊勇飛,平元峰,張文全. 基于Hilbert-Huang變換的混流泵流動誘導(dǎo)振動試驗[J]. 農(nóng)業(yè)工程學(xué)報,2018,34(2):47-54. doi:10.11975/j.issn.1002-6819.2018.02.007 http://www.tcsae.org

    Li Wei, Ji Leilei, Shi Weidong, Yang Yongfei, Ping Yuanfeng, Zhang Wenquan. Experiment of flow induced vibration of mixed-flow pump based on Hilbert-Huang transform[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(2): 47-54. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2018.02.007 http://www.tcsae.org

    猜你喜歡
    方向振動信號
    振動的思考
    2022年組稿方向
    2021年組稿方向
    信號
    鴨綠江(2021年35期)2021-04-19 12:24:18
    2021年組稿方向
    完形填空二則
    振動與頻率
    基于FPGA的多功能信號發(fā)生器的設(shè)計
    電子制作(2018年11期)2018-08-04 03:25:42
    中立型Emden-Fowler微分方程的振動性
    基于LabVIEW的力加載信號采集與PID控制
    啦啦啦中文免费视频观看日本| 日本欧美国产在线视频| 一级毛片aaaaaa免费看小| av不卡在线播放| 久久99蜜桃精品久久| 国产成人freesex在线| 麻豆精品久久久久久蜜桃| a级毛片免费高清观看在线播放| 制服丝袜香蕉在线| 久久久久久久亚洲中文字幕| 麻豆成人午夜福利视频| 九色成人免费人妻av| 亚洲精品日韩在线中文字幕| 精品一品国产午夜福利视频| 色婷婷久久久亚洲欧美| 国产片特级美女逼逼视频| 国产综合精华液| 亚洲欧美成人综合另类久久久| 国产熟女欧美一区二区| 黄色配什么色好看| 伦精品一区二区三区| 久久久久久人妻| 精品卡一卡二卡四卡免费| 亚洲内射少妇av| 国产黄色免费在线视频| 伦理电影大哥的女人| 亚洲在久久综合| 亚洲av成人精品一区久久| 国产精品蜜桃在线观看| 天美传媒精品一区二区| 一边亲一边摸免费视频| 美女主播在线视频| a级毛片在线看网站| 黄色一级大片看看| 免费少妇av软件| 特大巨黑吊av在线直播| 少妇人妻精品综合一区二区| 夜夜骑夜夜射夜夜干| 91精品伊人久久大香线蕉| 综合色丁香网| 亚洲欧洲国产日韩| 久久人妻熟女aⅴ| 亚洲第一av免费看| 午夜91福利影院| 亚洲欧洲日产国产| 久久99精品国语久久久| 成人亚洲精品一区在线观看| 99精国产麻豆久久婷婷| 日韩成人伦理影院| 特大巨黑吊av在线直播| 午夜av观看不卡| 麻豆成人av视频| av黄色大香蕉| 国产精品蜜桃在线观看| 国产亚洲精品久久久com| 噜噜噜噜噜久久久久久91| 2018国产大陆天天弄谢| 婷婷色综合大香蕉| 国产成人精品久久久久久| 91久久精品电影网| 伊人亚洲综合成人网| 高清视频免费观看一区二区| 成年av动漫网址| 亚洲精品456在线播放app| 青春草国产在线视频| 啦啦啦啦在线视频资源| 丝袜在线中文字幕| 三上悠亚av全集在线观看 | 新久久久久国产一级毛片| 久久精品国产a三级三级三级| 又爽又黄a免费视频| 中文字幕亚洲精品专区| 欧美丝袜亚洲另类| 精品久久久精品久久久| 成人18禁高潮啪啪吃奶动态图 | 亚洲欧美精品自产自拍| 99精国产麻豆久久婷婷| 欧美亚洲 丝袜 人妻 在线| 亚洲美女视频黄频| 亚洲精品国产色婷婷电影| 成人国产av品久久久| 一区二区三区精品91| 久久久亚洲精品成人影院| av福利片在线观看| 国产伦理片在线播放av一区| 亚洲无线观看免费| 精品人妻熟女毛片av久久网站| 日日摸夜夜添夜夜添av毛片| 成人午夜精彩视频在线观看| 中文字幕免费在线视频6| 久久人人爽人人片av| 热99国产精品久久久久久7| av一本久久久久| 亚洲一区二区三区欧美精品| kizo精华| 夫妻午夜视频| 女人精品久久久久毛片| 22中文网久久字幕| 成人免费观看视频高清| 国产成人精品一,二区| 亚洲成人一二三区av| 午夜老司机福利剧场| 成人国产麻豆网| 亚洲美女视频黄频| 久久99精品国语久久久| 午夜福利在线观看免费完整高清在| 在线观看一区二区三区激情| 秋霞在线观看毛片| 亚洲欧美清纯卡通| 婷婷色综合www| 91精品伊人久久大香线蕉| 在线观看美女被高潮喷水网站| 亚洲av电影在线观看一区二区三区| 国产精品久久久久久精品电影小说| 另类精品久久| 欧美精品亚洲一区二区| 亚洲精品国产av成人精品| 最近中文字幕高清免费大全6| 黑人猛操日本美女一级片| 久久人人爽人人爽人人片va| 一本色道久久久久久精品综合| 国产 一区精品| 日本av免费视频播放| 人妻夜夜爽99麻豆av| 日本-黄色视频高清免费观看| 国产av精品麻豆| 99国产精品免费福利视频| 午夜免费男女啪啪视频观看| 亚洲av.av天堂| 一本色道久久久久久精品综合| 国产高清有码在线观看视频| 国产av码专区亚洲av| 久久午夜综合久久蜜桃| 亚洲欧美精品自产自拍| 91精品伊人久久大香线蕉| 亚洲成人av在线免费| 亚洲国产色片| 最近最新中文字幕免费大全7| 欧美激情国产日韩精品一区| 欧美精品亚洲一区二区| 2018国产大陆天天弄谢| 全区人妻精品视频| 精品亚洲成a人片在线观看| 男男h啪啪无遮挡| 如日韩欧美国产精品一区二区三区 | 国产成人freesex在线| 国产精品三级大全| 女人精品久久久久毛片| 色吧在线观看| 久久久精品94久久精品| 女的被弄到高潮叫床怎么办| 人体艺术视频欧美日本| 亚洲精品国产色婷婷电影| 成人二区视频| av福利片在线观看| 91久久精品国产一区二区三区| 韩国高清视频一区二区三区| 国产av精品麻豆| 国产精品人妻久久久影院| 女性被躁到高潮视频| 久久婷婷青草| av福利片在线| 曰老女人黄片| 黄色配什么色好看| 精品久久久久久久久av| 免费黄网站久久成人精品| 国产av码专区亚洲av| 国产精品欧美亚洲77777| 另类精品久久| 欧美xxⅹ黑人| 日韩成人av中文字幕在线观看| 亚洲精品国产av蜜桃| 欧美区成人在线视频| 久久久国产一区二区| 99热6这里只有精品| 国产探花极品一区二区| 在线观看三级黄色| 亚洲国产欧美日韩在线播放 | 在现免费观看毛片| 不卡视频在线观看欧美| 国产成人精品婷婷| 纯流量卡能插随身wifi吗| 熟女电影av网| 有码 亚洲区| 亚洲精品456在线播放app| 久久久久久久久久人人人人人人| 成人影院久久| 黑人巨大精品欧美一区二区蜜桃 | 国产黄片美女视频| 麻豆成人av视频| 极品人妻少妇av视频| 国产精品嫩草影院av在线观看| 人人妻人人看人人澡| 最黄视频免费看| 最近手机中文字幕大全| 国产精品99久久久久久久久| 国产精品人妻久久久影院| 韩国高清视频一区二区三区| 男女免费视频国产| 久久6这里有精品| 日本91视频免费播放| 国产永久视频网站| 久久午夜福利片| 成年女人在线观看亚洲视频| 少妇人妻 视频| 老司机亚洲免费影院| 日韩成人伦理影院| 亚洲精品色激情综合| 麻豆精品久久久久久蜜桃| 黄色日韩在线| 一级二级三级毛片免费看| 亚洲不卡免费看| 两个人的视频大全免费| videossex国产| 久久这里有精品视频免费| 亚洲熟女精品中文字幕| 久久久久精品性色| .国产精品久久| 有码 亚洲区| 中文字幕免费在线视频6| 大话2 男鬼变身卡| 国产精品久久久久久av不卡| 亚洲欧美成人精品一区二区| 亚洲自偷自拍三级| 亚洲美女搞黄在线观看| 欧美亚洲 丝袜 人妻 在线| 亚洲婷婷狠狠爱综合网| 色5月婷婷丁香| 亚洲av男天堂| 国产精品久久久久久久久免| 国内少妇人妻偷人精品xxx网站| 嫩草影院新地址| 又黄又爽又刺激的免费视频.| 亚洲欧美精品专区久久| 亚洲国产色片| 一级毛片黄色毛片免费观看视频| 高清不卡的av网站| 97超碰精品成人国产| 亚洲美女视频黄频| 在线观看www视频免费| 国产精品久久久久久久电影| av国产精品久久久久影院| 久久狼人影院| av免费在线看不卡| 欧美日韩亚洲高清精品| 午夜日本视频在线| 在现免费观看毛片| 99精国产麻豆久久婷婷| 国产高清有码在线观看视频| 欧美老熟妇乱子伦牲交| 又粗又硬又长又爽又黄的视频| 亚洲av综合色区一区| 大话2 男鬼变身卡| 亚洲综合精品二区| av有码第一页| 久久综合国产亚洲精品| 日韩一区二区视频免费看| 久久久久精品久久久久真实原创| 亚洲成人一二三区av| 国产成人精品无人区| 高清不卡的av网站| 亚洲,欧美,日韩| 大码成人一级视频| 乱码一卡2卡4卡精品| 国产精品无大码| av女优亚洲男人天堂| 国产在线免费精品| 亚洲情色 制服丝袜| 久久久久久久亚洲中文字幕| 精品久久久久久电影网| 香蕉精品网在线| 国产女主播在线喷水免费视频网站| av在线播放精品| 精品一品国产午夜福利视频| 日韩大片免费观看网站| 免费黄网站久久成人精品| 欧美亚洲 丝袜 人妻 在线| 日本欧美国产在线视频| 99视频精品全部免费 在线| 黄色怎么调成土黄色| 精品一品国产午夜福利视频| xxx大片免费视频| 欧美 日韩 精品 国产| 日韩av在线免费看完整版不卡| 丰满少妇做爰视频| 免费人成在线观看视频色| 亚洲真实伦在线观看| 亚洲四区av| 精品一区二区三卡| 久久精品久久久久久噜噜老黄| 少妇熟女欧美另类| 伊人亚洲综合成人网| 人妻人人澡人人爽人人| 五月开心婷婷网| 久久ye,这里只有精品| 精品亚洲成国产av| 亚洲国产av新网站| 中文字幕亚洲精品专区| 午夜av观看不卡| 99热6这里只有精品| 亚洲国产精品一区二区三区在线| 高清欧美精品videossex| www.av在线官网国产| 成人漫画全彩无遮挡| 夜夜爽夜夜爽视频| 爱豆传媒免费全集在线观看| av在线app专区| 观看av在线不卡| 亚洲人与动物交配视频| 成人毛片a级毛片在线播放| videossex国产| av有码第一页| 国产免费视频播放在线视频| 免费大片黄手机在线观看| av天堂久久9| 免费少妇av软件| 日韩制服骚丝袜av| 自线自在国产av| 午夜av观看不卡| 国产视频内射| 国产男女超爽视频在线观看| av免费在线看不卡| 99热全是精品| 伦精品一区二区三区| av福利片在线| 免费大片18禁| 亚洲不卡免费看| 少妇高潮的动态图| 日韩精品有码人妻一区| 性色avwww在线观看| 日韩精品免费视频一区二区三区 | 美女cb高潮喷水在线观看| 国产淫语在线视频| 91午夜精品亚洲一区二区三区| 成人漫画全彩无遮挡| 99热这里只有精品一区| 欧美三级亚洲精品| 亚洲欧洲日产国产| a级毛色黄片| 欧美日韩综合久久久久久| 各种免费的搞黄视频| 亚洲精品第二区| 看免费成人av毛片| 一本—道久久a久久精品蜜桃钙片| 一本一本综合久久| 亚洲国产精品999| 免费人妻精品一区二区三区视频| 久久久久久久久久成人| 街头女战士在线观看网站| 插阴视频在线观看视频| 高清视频免费观看一区二区| 中文欧美无线码| 一本大道久久a久久精品| 色5月婷婷丁香| 国产精品秋霞免费鲁丝片| 少妇 在线观看| 99精国产麻豆久久婷婷| 777米奇影视久久| 欧美 日韩 精品 国产| 亚洲精品乱码久久久久久按摩| 91精品伊人久久大香线蕉| 中文字幕制服av| 亚洲欧美精品自产自拍| 一区二区三区四区激情视频| 精品国产国语对白av| 亚洲,欧美,日韩| av一本久久久久| 丝袜喷水一区| www.色视频.com| 丰满饥渴人妻一区二区三| 九九在线视频观看精品| 观看美女的网站| 午夜免费男女啪啪视频观看| 日本黄大片高清| 一二三四中文在线观看免费高清| 欧美人与善性xxx| 精品久久久久久电影网| 男人狂女人下面高潮的视频| 日韩人妻高清精品专区| 一区二区三区免费毛片| 在线观看国产h片| 在线观看免费视频网站a站| 精品国产一区二区久久| 丝袜脚勾引网站| 夫妻性生交免费视频一级片| 国产精品一区二区在线不卡| 国产精品麻豆人妻色哟哟久久| 观看av在线不卡| 99热6这里只有精品| 午夜福利在线观看免费完整高清在| 亚洲欧美日韩卡通动漫| 80岁老熟妇乱子伦牲交| 国产男人的电影天堂91| 乱人伦中国视频| 中文欧美无线码| 狂野欧美激情性bbbbbb| 最近中文字幕2019免费版| 欧美日韩视频高清一区二区三区二| 国产淫语在线视频| av免费在线看不卡| 极品教师在线视频| 2022亚洲国产成人精品| 亚洲国产色片| 国产免费福利视频在线观看| 亚洲婷婷狠狠爱综合网| 啦啦啦在线观看免费高清www| 国产精品久久久久久久久免| 这个男人来自地球电影免费观看 | 女的被弄到高潮叫床怎么办| 国产高清不卡午夜福利| 亚洲av成人精品一二三区| 午夜激情福利司机影院| 天堂中文最新版在线下载| av国产精品久久久久影院| 久久精品国产亚洲网站| 国产亚洲av片在线观看秒播厂| 我的老师免费观看完整版| 91精品国产国语对白视频| 国产成人一区二区在线| 777米奇影视久久| 老熟女久久久| 一级,二级,三级黄色视频| 婷婷色综合www| 国产精品一区二区三区四区免费观看| 青春草亚洲视频在线观看| 国产免费又黄又爽又色| 日日撸夜夜添| 国产91av在线免费观看| 国产 一区精品| 狂野欧美白嫩少妇大欣赏| 99视频精品全部免费 在线| 偷拍熟女少妇极品色| 精品一区二区三卡| 亚洲精品日韩在线中文字幕| 汤姆久久久久久久影院中文字幕| 嫩草影院入口| 人妻 亚洲 视频| 一本—道久久a久久精品蜜桃钙片| av黄色大香蕉| 欧美bdsm另类| 18禁动态无遮挡网站| 汤姆久久久久久久影院中文字幕| 精品卡一卡二卡四卡免费| kizo精华| 久久久久精品久久久久真实原创| 免费看av在线观看网站| 街头女战士在线观看网站| 国产乱来视频区| av在线观看视频网站免费| 最近的中文字幕免费完整| 天天操日日干夜夜撸| 91久久精品国产一区二区成人| 国产探花极品一区二区| 国产黄片视频在线免费观看| 中国美白少妇内射xxxbb| 又粗又硬又长又爽又黄的视频| 久久久国产欧美日韩av| 欧美日韩视频精品一区| 日韩一区二区三区影片| 国产成人精品久久久久久| 午夜福利影视在线免费观看| 亚洲成色77777| 欧美日韩视频高清一区二区三区二| 日韩欧美一区视频在线观看 | 亚洲av成人精品一区久久| 欧美精品一区二区大全| 亚洲国产最新在线播放| 亚洲,欧美,日韩| 国国产精品蜜臀av免费| 在线观看人妻少妇| 天天躁夜夜躁狠狠久久av| 我要看黄色一级片免费的| 亚洲无线观看免费| 五月玫瑰六月丁香| 亚洲av在线观看美女高潮| 伦理电影大哥的女人| 天堂8中文在线网| 另类亚洲欧美激情| 中文精品一卡2卡3卡4更新| 精品99又大又爽又粗少妇毛片| 制服丝袜香蕉在线| 卡戴珊不雅视频在线播放| 女性生殖器流出的白浆| 久久精品久久精品一区二区三区| 国产伦在线观看视频一区| 国产精品欧美亚洲77777| 美女中出高潮动态图| 777米奇影视久久| 亚洲国产精品国产精品| 哪个播放器可以免费观看大片| 少妇熟女欧美另类| 美女大奶头黄色视频| 亚洲国产最新在线播放| 久久鲁丝午夜福利片| 美女xxoo啪啪120秒动态图| 青春草视频在线免费观看| 亚洲欧美成人综合另类久久久| 韩国av在线不卡| 成人亚洲精品一区在线观看| 久久青草综合色| 欧美+日韩+精品| 18禁在线播放成人免费| 欧美高清成人免费视频www| 国产精品免费大片| 成年女人在线观看亚洲视频| 日韩大片免费观看网站| 免费在线观看成人毛片| 欧美精品人与动牲交sv欧美| 热99国产精品久久久久久7| 全区人妻精品视频| 能在线免费看毛片的网站| 九草在线视频观看| av国产精品久久久久影院| 国产精品人妻久久久久久| av一本久久久久| 免费人成在线观看视频色| 久久青草综合色| 日本黄大片高清| 美女视频免费永久观看网站| 成人国产av品久久久| 91精品国产国语对白视频| 一个人看视频在线观看www免费| 亚洲怡红院男人天堂| 亚洲av综合色区一区| 亚洲天堂av无毛| 亚洲欧美日韩另类电影网站| 国产黄色免费在线视频| 久久久久久人妻| 精品久久久久久久久av| 大香蕉久久网| 亚洲欧美一区二区三区国产| 九九久久精品国产亚洲av麻豆| 国产欧美日韩精品一区二区| 毛片一级片免费看久久久久| 欧美精品人与动牲交sv欧美| 免费黄网站久久成人精品| √禁漫天堂资源中文www| 国产精品一区二区三区四区免费观看| 熟女人妻精品中文字幕| 国产极品粉嫩免费观看在线 | 97超碰精品成人国产| 欧美日韩国产mv在线观看视频| 国产精品国产三级专区第一集| videos熟女内射| 狂野欧美激情性xxxx在线观看| 亚洲成人av在线免费| 一级黄片播放器| 亚洲美女视频黄频| 亚洲精品一二三| 亚洲欧美清纯卡通| 国产精品偷伦视频观看了| tube8黄色片| 日本爱情动作片www.在线观看| 天堂中文最新版在线下载| 香蕉精品网在线| 日本猛色少妇xxxxx猛交久久| www.av在线官网国产| 欧美精品亚洲一区二区| 少妇人妻久久综合中文| 国产伦在线观看视频一区| 在现免费观看毛片| 久久精品久久久久久噜噜老黄| 午夜视频国产福利| 精品少妇内射三级| av网站免费在线观看视频| 久久精品国产a三级三级三级| 国产高清国产精品国产三级| 日本免费在线观看一区| 久热久热在线精品观看| 色哟哟·www| 老司机影院毛片| 国产乱来视频区| 99久久中文字幕三级久久日本| 久久久欧美国产精品| 国模一区二区三区四区视频| 国产午夜精品久久久久久一区二区三区| 观看免费一级毛片| 成人18禁高潮啪啪吃奶动态图 | 亚洲精品视频女| 久热久热在线精品观看| 国产精品欧美亚洲77777| 熟女av电影| 国产免费一区二区三区四区乱码| 国产亚洲5aaaaa淫片| 久久久国产欧美日韩av| av国产精品久久久久影院| av网站免费在线观看视频| 97在线人人人人妻| 黄色怎么调成土黄色| 在线观看av片永久免费下载| 老熟女久久久| 观看美女的网站| 交换朋友夫妻互换小说| 国产深夜福利视频在线观看| 在线观看国产h片| 欧美三级亚洲精品| 如何舔出高潮| 国国产精品蜜臀av免费| 国产精品久久久久久久电影| 在线观看www视频免费| 各种免费的搞黄视频| 伦理电影大哥的女人| 国产精品一二三区在线看| 亚洲精品乱久久久久久| 国产探花极品一区二区| 麻豆成人午夜福利视频| 亚洲av男天堂| 国产欧美亚洲国产| 麻豆乱淫一区二区| 中文字幕免费在线视频6| 成人影院久久| 国产成人精品无人区| 十八禁高潮呻吟视频 | 老司机影院毛片| 国产无遮挡羞羞视频在线观看|