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

    旋轉(zhuǎn)機(jī)械扭振檢測方法的研究現(xiàn)狀及展望

    2021-03-14 12:18:04郭艷玲邱楓李志鵬
    關(guān)鍵詞:檢測方法發(fā)展趨勢

    郭艷玲 邱楓 李志鵬

    摘 要:旋轉(zhuǎn)機(jī)械扭振是影響旋轉(zhuǎn)機(jī)械工作可靠性和使用壽命的主要因素之一,近年來由于扭振引起的大型機(jī)組事故日益增多,因此如何檢測扭振成為機(jī)械行業(yè)一大重點(diǎn)。扭振的檢測過程可分為扭振信號(hào)的測量和扭振信號(hào)的提取分析這兩個(gè)步驟。將扭振信號(hào)的測量方法分為接觸式測量法和非接觸式測量法兩類,將扭振信號(hào)的提取分析方法分為基于傅里葉變化法和基于信號(hào)分解法兩類,并對(duì)這幾類方法中的具體每種方法進(jìn)行了介紹,詳細(xì)地闡述了每種方法的工作原理、研究意義、實(shí)際應(yīng)用以及優(yōu)缺點(diǎn),最后對(duì)這幾類方法的發(fā)展趨勢進(jìn)行了展望,希望對(duì)今后的旋轉(zhuǎn)機(jī)械扭振檢測方法的研究起到一定的參考和指導(dǎo)作用。

    關(guān)鍵詞:旋轉(zhuǎn)機(jī)械;扭振;信號(hào)測量;信號(hào)提取分析;發(fā)展趨勢;檢測方法

    DOI:10.15938/j.jhust.2021.06.007

    中圖分類號(hào): TH825

    文獻(xiàn)標(biāo)志碼: A

    文章編號(hào): 1007-2683(2021)06-0047-08

    Research Status and Prospect of Torsional Vibration

    Detection Methods for Rotating Machinery

    GUO Yan-ling, QIU Feng, LI Zhi-peng

    (School of Electromechanic Engineering, Northeast Forestry University, Harbin 150040, China)

    Abstract:Torsional vibration of rotating machinery is one of the main factors affecting the reliability and service life of rotaing machinery. In recent years, the number of large units accidents caused by torsional vibration is increasing. Therefore, how to detect torsional vibration has become a major focus in the machinery increasing.The process of torsional vibration of detection method can be divided into two steps: the measurement of torsional vibration signals and the extraction and analysis of torsional vibration signal. The measurement methods of torsional vibration signal are divided into contact measurement method and non-contact measurement method. The extraction and analysis methods of torsional vibration signal are divided into Fourier transform method and signal decomposition method. We introduce each method of these methods, espound in detail the woking principle, reaserch significance, practical application, advantages and disadvantages of each methods, and finally forecast the development trend of these methods. It is expected to provide some reference and guidance for the future research on torsional vibration detection methods of rotating machinery.

    Keywords:rotating machinery; torsional vibration; signal measurement; signal extraction analysis; development trend; detection method

    0 引 言

    旋轉(zhuǎn)機(jī)械主要是指電動(dòng)機(jī)、汽輪機(jī)、燃?xì)廨啓C(jī)等機(jī)械設(shè)備。隨著此類機(jī)械的應(yīng)用越來越廣泛,如何確保其安全穩(wěn)定的運(yùn)行,引起越來越多的重視。扭振對(duì)旋轉(zhuǎn)機(jī)械造成很大影響,使旋轉(zhuǎn)機(jī)械疲勞積累不斷加強(qiáng),形成裂紋、切口,并逐漸擴(kuò)散,最終導(dǎo)致此前不容易察覺的軸系的斷裂。同時(shí),扭振會(huì)引起系統(tǒng)噪聲增大、動(dòng)力性能下降。

    近年來國內(nèi)外由于扭振原因引起的大型機(jī)組事故日益增多,2003-2004年,俄羅斯RAO EES Rossii電廠因扭振引起的渦輪轉(zhuǎn)子出現(xiàn)裂紋的事故共12起,修理更換轉(zhuǎn)子的費(fèi)用達(dá)到3.08億盧布[1]。2009年由于扭振造成烏茲別克斯坦塔里曼扎爾州地區(qū)發(fā)電站的800兆瓦渦輪機(jī)組被完全摧毀[1]。2010年廣東茂名某熱電廠200MW機(jī)組鍋爐配備的風(fēng)機(jī)主軸斷裂,2011年該軸再次斷裂,事故原因是扭轉(zhuǎn)疲勞剪切斷裂[2]。此類事故的發(fā)生,引起行業(yè)內(nèi)的廣泛重視。

    扭轉(zhuǎn)振動(dòng)簡稱扭振,其主要原因是喪失了動(dòng)力力矩和負(fù)載扭矩之間的平衡,合成扭矩的方向反復(fù)變化。由于旋轉(zhuǎn)軸不是絕對(duì)的剛體,當(dāng)旋轉(zhuǎn)軸受到動(dòng)力源周期性的扭轉(zhuǎn)激勵(lì)時(shí),旋轉(zhuǎn)軸就會(huì)按動(dòng)力源激振頻率強(qiáng)迫振動(dòng)。對(duì)扭振的檢測可以分為兩個(gè)步驟:扭振信號(hào)的測量和扭振信號(hào)的提取分析。

    1 扭振信號(hào)測量方法

    扭振測量技術(shù)的起源可追溯至19世紀(jì)初,發(fā)展至今,主要可分為兩大類:接觸式和非接觸式扭振測量法。

    1.1 接觸式測量方法

    接觸式測量法主要特點(diǎn)是傳感元件安裝于旋轉(zhuǎn)軸上,通過無線收發(fā)裝置或者滑環(huán)來收集所測得的扭振信號(hào)。其采用的方法主要有電阻應(yīng)變片法、加速度法和光纖光柵應(yīng)變片法。

    1.1.1 電阻應(yīng)變片法

    電阻應(yīng)變片法[3]起源于19世紀(jì)初,最早用于旋轉(zhuǎn)軸扭矩的測量,后也可以用于扭振的測量。該方法將電阻應(yīng)變片互成直角粘貼于被測軸上,如圖1所示。當(dāng)被測旋轉(zhuǎn)軸產(chǎn)生彈性形變時(shí),應(yīng)變片的電阻值會(huì)產(chǎn)生相應(yīng)的變化,阻值的變化量與被測軸形變量成一定的比例,通過這一比例可以計(jì)算出被測軸形變的大小,從而得到扭振信息。Han H S等[4]采用本方法對(duì)某船軸系扭振進(jìn)行了實(shí)測,驗(yàn)證了該方法的有效性。

    電阻應(yīng)變片法可以直接測量被測軸形變量的大小,具有靈敏度高、結(jié)構(gòu)簡單等優(yōu)點(diǎn),適用于中低速工況。但電阻應(yīng)變片的使用壽命有限,不適合在惡劣的條件下工作。

    1.1.2 加速度計(jì)法

    加速度法[5]是一種用壓電加速度計(jì)作為傳感器安裝在旋轉(zhuǎn)軸上來測量扭振的方法。該方法將兩個(gè)加速度計(jì)均勻布置在被測軸的兩端。通過加速度計(jì)來測量被測軸的角加速度信號(hào),對(duì)其進(jìn)行積分得到角速度信號(hào),角速度信號(hào)即為扭振信號(hào)。

    本方法靈敏度高,成本低,適用于中低速工況。相對(duì)于電阻應(yīng)變片法提升了使用壽命且結(jié)構(gòu)更為簡單。缺點(diǎn)是:在工作條件受到限制時(shí),不易于安裝在主軸上。

    1.1.3 光纖光柵應(yīng)變片法

    近年來眾多專家學(xué)者針對(duì)電阻應(yīng)變片不適合惡劣工況,使用壽命短,以及加速度計(jì)不易安裝等缺點(diǎn),致力于研究新材料檢測扭振,其中光纖光柵應(yīng)變片法最為主流。

    本方法是一種基于光纖光柵原理的扭振測量方法,魏莉、劉芹等人設(shè)計(jì)了光纖光柵與扭振之間的應(yīng)變傳遞模型,并根據(jù)此模型研究了光纖光柵傳感器并搭建了實(shí)驗(yàn)臺(tái),利用差分法消除了溫度對(duì)光柵應(yīng)變片的影響,使其只對(duì)應(yīng)變敏感。當(dāng)轉(zhuǎn)軸產(chǎn)生扭振時(shí),通過光柵應(yīng)變量與轉(zhuǎn)軸扭矩變量之間的關(guān)系得到扭振信息[6]。

    由于光纖光柵傳感器,外形和質(zhì)量小,安裝簡單,幾乎不影響軸系的動(dòng)平衡,同時(shí)不受電磁波的干擾,抗腐性強(qiáng) [7],克服了電阻應(yīng)變片壽命短,不宜在惡劣工況下工作和加速度計(jì)不易安裝的缺點(diǎn),但該測量技術(shù)成本高,其測量精度受環(huán)境影響較大。

    綜上3種方法相比較而言,電阻應(yīng)變片法目前使用最為廣泛,光纖光柵應(yīng)變片法目前還處于實(shí)驗(yàn)階段并不成熟,加速度計(jì)法由于加速度計(jì)安裝受工況限制,未能得到廣泛應(yīng)用。

    接觸式測量法的一大優(yōu)點(diǎn)是,可以通過傳感器直接測得被測軸的扭振。但該方法卻有一個(gè)難以改善的缺點(diǎn),即不適合測量高轉(zhuǎn)速工況,因?yàn)楸粶y轉(zhuǎn)軸的轉(zhuǎn)速一旦過高,轉(zhuǎn)軸離心力過大,安裝在被測軸上的傳感元件有脫離被測轉(zhuǎn)軸的危險(xiǎn)。所以接觸式測量技術(shù)漸漸有被非接觸式測量技術(shù)取代的趨勢。

    1.2 非接觸式扭振測量方法

    非接觸測量方法的主要特點(diǎn)是傳感器不與被測軸系直接相連接,適用于高速旋轉(zhuǎn)軸,通常通過測量旋轉(zhuǎn)軸的角速度來體現(xiàn)被測軸的扭振信息。

    1.2.1 脈沖時(shí)序法

    脈沖時(shí)序法起源于19世紀(jì)90年代,屬于最早的非接觸式測量技術(shù),是綜合相位差法和頻數(shù)計(jì)數(shù)法的特點(diǎn)發(fā)展而來的。

    脈沖時(shí)序法有兩大類:一種是磁電式脈沖時(shí)序法[8-10],在旋轉(zhuǎn)軸上安裝等分齒盤,磁電式傳感器充當(dāng)測量工具,當(dāng)轉(zhuǎn)軸工作時(shí),傳感器通過感應(yīng)磁通量的變化,產(chǎn)生一系列關(guān)于電壓的周期性諧波[11-12]。另一種是光電式脈沖時(shí)序法[13-16] ,即利用光電式傳感器,對(duì)旋轉(zhuǎn)軸上粘貼的黑白相間的編碼帶進(jìn)行測量,由于黑白條紋對(duì)于光的反射率不同,在傳感器上產(chǎn)生相應(yīng)高電平信號(hào)和低電平信號(hào),也可測的一系列周期性諧波。

    光電式測量法,靈敏度相對(duì)于磁電式更高,但是不易于應(yīng)用于惡劣的工作條件。磁電式測量適合在惡劣條件下,但在旋轉(zhuǎn)軸上安裝等分齒盤較為麻煩,而且等分齒盤具有一定的質(zhì)量,可能會(huì)引起旋轉(zhuǎn)軸的結(jié)構(gòu)特性改變。

    A. L. Nazolin等[15]基于光電式脈沖時(shí)序法原理建立一個(gè)扭振監(jiān)測系統(tǒng),對(duì)某300MW的汽輪機(jī)組進(jìn)行了扭振監(jiān)測,驗(yàn)證了該方法的有效性,并給出了后續(xù)一些列的改進(jìn)性建議,推動(dòng)了此技術(shù)的發(fā)展。

    張玉皓等[17]提出了一種脈沖信號(hào)輸出的廣義增量編碼器模型,用于測量旋轉(zhuǎn)機(jī)械的瞬時(shí)角速度,優(yōu)化了傳統(tǒng)脈沖時(shí)序法的實(shí)時(shí)性。并在國內(nèi)2臺(tái)1000MW的汽輪機(jī)發(fā)電機(jī)組上實(shí)施了檢測,驗(yàn)證了該方法的有效性。

    脈沖時(shí)序法目前在行業(yè)內(nèi)應(yīng)用最為廣泛,代表產(chǎn)品有英國的TV-1、美國的TVSC型扭振儀等。與此同時(shí)這一測量技術(shù)也存在一些明顯的缺點(diǎn)需要完善:

    1)條紋碼,等分齒輪存在制造誤差[18],或是轉(zhuǎn)速產(chǎn)生波動(dòng)時(shí),會(huì)造成測試結(jié)果失真,進(jìn)而產(chǎn)生測量誤差。

    2)本方法不適合測量旋轉(zhuǎn)軸的低速工況,近幾年為了改善上述缺點(diǎn)眾多專家學(xué)者進(jìn)行了大量的研究和實(shí)驗(yàn)。

    DIAMOND D H等[19]提出了一種計(jì)算方法, 該方法采用了貝葉斯線性回歸來計(jì)算條紋帶的增量距離,并經(jīng)過模擬實(shí)驗(yàn)得到驗(yàn)證,結(jié)果表明,該方法可以準(zhǔn)確地確定轉(zhuǎn)軸上任意不規(guī)則形狀的條紋帶,有效的彌補(bǔ)了光電時(shí)序法由于條紋帶不規(guī)則形狀帶來的誤差。

    田忠旭等[18]針對(duì)磁電式傳感器采集的類似余弦波的電壓信號(hào),提出了一種不斷向零電壓靠近的雙線逼近算法,給出了轉(zhuǎn)速計(jì)算的諧次提取算法,研究了瞬時(shí)轉(zhuǎn)速各階諧波函數(shù)的變化規(guī)律,可以有效去除由等分齒輪誤差所造成的信號(hào)壞點(diǎn)。經(jīng)過在四缸汽油機(jī)的扭振測試結(jié)果表明,該方法具有較高的適應(yīng)性與穩(wěn)定性。LI K等[20]設(shè)計(jì)了一款基于電磁感應(yīng)信號(hào)的角度傳感器。利用了感應(yīng)線圈對(duì)旋轉(zhuǎn)磁場的反饋原理,實(shí)現(xiàn)了同時(shí)滿足同時(shí)適用于高轉(zhuǎn)速和低轉(zhuǎn)速角度測量。并通過實(shí)驗(yàn)驗(yàn)證得出結(jié)論,該角度傳感器可以測量1 r/min到100 r/min的速度范圍,為磁電式脈沖時(shí)序法測量低轉(zhuǎn)速工況邁出了重要的一步。

    FU Y.C等[21]設(shè)計(jì)了一款電容傳感器,其探頭采用差動(dòng)多層環(huán)形結(jié)構(gòu),增加了微板的正面積提高了傳感器的測量精度及范圍。分別進(jìn)行了大量程和小量程的標(biāo)定實(shí)驗(yàn),分析了不同的擬合方法的非線性誤差,驗(yàn)證了設(shè)計(jì)的合理性,為將來磁電式脈沖時(shí)序法傳感器的選用方面提供了有價(jià)值的參考。

    ONA Denis Ijike等[22]通過對(duì)傳統(tǒng)電渦流傳感器探頭中線圈間隙的調(diào)整,提高了電渦流傳感器探頭的靈敏度,這一方法對(duì)今后的磁電式脈沖時(shí)序法精度提高的有著重大的意義。

    1.2.2 激光測扭法

    激光測扭法[23-24]是一種基于激光多普勒技術(shù)的扭振測量方法,將激光投射到被測軸的一個(gè)截面上,反射光形成多普勒頻移,該頻移與旋轉(zhuǎn)軸轉(zhuǎn)速成正比。通過光學(xué)配置如圖2,使相鄰兩個(gè)時(shí)刻所形成的多普勒頻移光信號(hào)在光學(xué)探測器上產(chǎn)生光學(xué)混頻,光電流時(shí)間差正比于轉(zhuǎn)軸在兩個(gè)時(shí)刻的速度差,由于相鄰的光電流時(shí)間差極小,因此可以直接獲得角加速度,這一角加速度反映了軸系角速度的變化,進(jìn)而獲取軸系的扭振的信息。

    1.2.3 CCD扭振測量法

    基于CCD技術(shù)的扭振技術(shù)測量[25]:該方法首先需要在被測旋轉(zhuǎn)軸上標(biāo)記一條螺旋線,作為CCD攝像機(jī)的被測物。旋轉(zhuǎn)軸工作時(shí),螺旋線會(huì)在CCD攝像機(jī)上呈一系列斑點(diǎn),當(dāng)扭振存在時(shí),斑點(diǎn)的分布是非均勻的,即斑點(diǎn)的運(yùn)動(dòng)不處于勻速狀態(tài),通過計(jì)算斑點(diǎn)的瞬時(shí)角加速度,通過數(shù)據(jù)處理,從而獲得旋轉(zhuǎn)軸的扭振信息。

    本方法的優(yōu)點(diǎn)是精度高,實(shí)時(shí)性好,穩(wěn)定性強(qiáng),信息量大,可一次性測量多個(gè)轉(zhuǎn)軸面。但是目前為止,該方法未能應(yīng)用于實(shí)際工程當(dāng)中,沒有相應(yīng)產(chǎn)品問世,可靠性方面沒有得到驗(yàn)證。

    非接觸式測量技術(shù)是當(dāng)今行業(yè)內(nèi)的主流技術(shù),綜上所述三種非接觸式扭振測量技術(shù),脈沖時(shí)序法測量穩(wěn)定,工程上易實(shí)現(xiàn),現(xiàn)有產(chǎn)品最多,應(yīng)用范圍最廣;激光測扭法靈敏度高,響應(yīng)速度快,但是價(jià)格昂貴。CCD測量技術(shù)精度高,實(shí)時(shí)性好,但工程應(yīng)用少,可靠性未得到驗(yàn)證。

    2 扭振測量信號(hào)提取分析方法

    無論采用以上哪一種扭振測量方法,其最終測得的結(jié)果都是一個(gè)輸出信號(hào),得到相應(yīng)的信號(hào)如何消除噪聲,提取有用信息也是研究的一大重點(diǎn)[26]。

    傳統(tǒng)的傅里葉變換法應(yīng)用廣泛,但僅適用于線性穩(wěn)定系統(tǒng),對(duì)于非線性非穩(wěn)定系統(tǒng),傳統(tǒng)傅里葉變換法則受到了很大的局限性[27]。對(duì)于非線性非穩(wěn)定信號(hào),目前多數(shù)使用時(shí)頻分析的方法來進(jìn)行分析。

    常用的時(shí)頻分析法主要包括兩類,一類是以傅里葉變換為基礎(chǔ)的方法,另一種是以經(jīng)驗(yàn)?zāi)B(tài)分解(empirical mode decomposition,EMD)為基礎(chǔ)的方法。

    2.1 基于傅里葉變換法

    此類方法早起最有代表性的是短時(shí)傅里葉變換和小波變換,這些方法能一定程度上解決了傅里葉變換在處理非線性非穩(wěn)定信號(hào)時(shí)的局限性,但是短時(shí)傅里葉變換的窗函數(shù)大小固定,無法同時(shí)滿足較高的時(shí)間分辨率和頻率分辨率,小波變換也面臨著小波基選取的問題,二者都不具有自適應(yīng)性[28-32]。

    近年來諸多專家學(xué)者針對(duì)這類問題展開了深入的研究,S變換,同步擠壓變換等方法相繼問世。

    2.1.1 S變換,廣義S變換

    S變換[33]是Stockwell于1996年提出的時(shí)頻分析算法, S變換采用高斯窗函數(shù),其寬度可以隨信號(hào)頻率變化,相比較其他方法無需選擇窗函數(shù),同時(shí)克服了短時(shí)傅里葉變換窗函數(shù)大小固定的缺陷,所以無論是對(duì)于信號(hào)中的低頻部分還是高頻部分,S變換過濾后的信號(hào)較短時(shí)傅里葉變換及小波變換所得到的信號(hào)都更為準(zhǔn)確,分辨率更高,無損性,降噪性能好。但是S變換中高斯窗口會(huì)隨著頻率的增大而變窄,進(jìn)而導(dǎo)致頻率的分辨率降低[34]。

    為了克服S變換的這一缺點(diǎn),Pinnegar等引入調(diào)節(jié)因子λ對(duì)高斯窗函數(shù)進(jìn)行改進(jìn),提出了廣義S變換,通過調(diào)節(jié)λ使高斯窗函數(shù)隨頻率的變化自適應(yīng)調(diào)整以獲得最佳時(shí)頻分辨率,適合處理非線性,非平穩(wěn)的振動(dòng)信號(hào)[35-36]。

    2.1.2 同步擠壓小波變換

    本方法是2011年Daubechies在小波變換的基礎(chǔ)上提出來的,根據(jù)時(shí)間-尺度平面每個(gè)元素絕對(duì)值的大小,對(duì)平面中的能量進(jìn)行重新分配,最后根據(jù)映射公式將時(shí)間-尺度平面轉(zhuǎn)化為時(shí)間-頻率平面[37-38]。

    本方法的優(yōu)點(diǎn)是信號(hào)時(shí)頻聚集性及分辨性強(qiáng),同時(shí)受不確定準(zhǔn)則的限制較小,因此該方法時(shí)頻聯(lián)合域比其他分析方法的時(shí)頻分析方法的聯(lián)合域要大,受到廣泛使用。與此同時(shí),本方法的結(jié)果也在一定程度上受到小波基的影響,反映高頻弱振幅的信息能力不足。

    2.1.3 同步擠壓S變換

    針對(duì)上述S變換和同步擠壓變換的不足,近年來提出了同步擠壓S變換,該方法將S變換和同步擠壓變換結(jié)合起來,無需選擇小波基,改善了同步擠壓小波變換結(jié)果受小波基影響的缺點(diǎn),同時(shí)也克服了S變換時(shí)頻分辨率不足的缺點(diǎn)。目前,該方法是基于傅里葉變換發(fā)展起來的一種最優(yōu)方法,在處理分析振動(dòng)信號(hào)中極為有效。

    潘高元等[39]運(yùn)用了該方法對(duì)美國Case Western Reserve University電氣工程實(shí)驗(yàn)室中的相關(guān)軸承的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了分析,證明了該方法優(yōu)點(diǎn)的真實(shí)性。

    以上基于傅里葉變換所演化過來的一類方法的共同特點(diǎn)是要選取預(yù)先設(shè)定好的基函數(shù)來表達(dá)原始信號(hào),存在自適應(yīng)不足,時(shí)頻分辨率低的問題,雖然同步擠壓S變換可以在很大程度上克服這些問題 ,但還是存在這方面的弊端。

    2.2 基于EMD信號(hào)分解法

    此類方法是基于經(jīng)驗(yàn)?zāi)B(tài)分解法 [40-41](empirical mode decomposition,EMD)的一類方法,EMD法是依據(jù)數(shù)據(jù)自身的時(shí)間尺度進(jìn)行信號(hào)分解,相對(duì)于基于傅里葉變換法,無需選取基函數(shù),具有自適應(yīng)性強(qiáng)的優(yōu)點(diǎn),非常適合處理非線性非平穩(wěn)信號(hào)。但是存在端點(diǎn)效應(yīng)及模態(tài)混疊現(xiàn)象,嚴(yán)重影響了分解結(jié)果的正確性。近年來,眾多專家學(xué)者在此方法的基礎(chǔ)上進(jìn)行改進(jìn)從而克服以上不足。

    2.2.1 集合經(jīng)驗(yàn)?zāi)B(tài)分解法,完全集合經(jīng)驗(yàn)?zāi)B(tài)分解法

    集合經(jīng)驗(yàn)?zāi)B(tài)分解法(ensemble empirical mode decomposition,EEMD)在EMD法的基礎(chǔ)上在信號(hào)中加入了白噪聲,白噪聲具有頻率均勻分布的統(tǒng)計(jì)特性,使信號(hào)在尺度分布上連續(xù)不斷,改變其在極點(diǎn)上的特性,增強(qiáng)了信號(hào)的內(nèi)捷性,有效地避免了模態(tài)混疊的現(xiàn)象,但是添加的白噪聲殘留會(huì)給信號(hào)帶來噪聲干擾。

    張超等[42]采用本方法對(duì)某鋼鐵公司電機(jī)的滾動(dòng)軸承做了故障診斷,將故障軸承的非穩(wěn)態(tài)原始振動(dòng)信號(hào)進(jìn)行了分析,并與EMD方法做了比較,可以更快的提取軸承損傷性故障特征。驗(yàn)證了該方法的有效性。

    完全集合經(jīng)驗(yàn)?zāi)B(tài)分解法[43](complete ensemble empirical mode decomposition,CEEMD)與EEMD的過程一樣,不同之處在于在添加一種白噪聲之后再加上一組相反符號(hào)的噪聲,再按照EMD分解得到固有模態(tài)函數(shù)(intrinsic mode function,IMF),該方法不僅解決了模態(tài)混疊問題,而且對(duì)原始信號(hào)進(jìn)行了精確的重構(gòu),也可以部分的抵消白噪聲帶來的噪聲干擾,是在EMMD方法的基礎(chǔ)上所做的改進(jìn)。

    2.2.2 希爾伯特-黃變換

    希爾伯特-黃變換(Hilbert-Huang transformation, HHT)方法起源于1998年[44],近年來被廣泛的應(yīng)用于處理振動(dòng)信號(hào)。該方法包含了經(jīng)驗(yàn)?zāi)B(tài)分解法和希爾伯特變換兩部分計(jì)算,其簡要過程是可分為如下步驟:

    1)采用經(jīng)驗(yàn)?zāi)B(tài)分解方法將測得的扭振信號(hào)分解為一些固有模態(tài)函數(shù)(IMF),這些函數(shù)是滿足一定條件的分量;

    2)對(duì)每一個(gè)IMF進(jìn)行Hilbert變換,得到相應(yīng)的 Hilbert譜,換言之就是將每個(gè)IMF表示在聯(lián)合的時(shí)頻域中;

    3)將所有IMF的Hilbert譜歸納在一起就會(huì)得到原信號(hào)的Hilbert譜。

    HHT在處理非線性非平穩(wěn)扭振信號(hào)上優(yōu)勢顯著,本方法由EMD和希爾伯特變換兩部分組成,EMD的模態(tài)混疊缺點(diǎn)一定程度上也會(huì)影響到HHT的精度[45]。

    2.2.3 希爾伯特振動(dòng)分解法

    希爾伯特振動(dòng)分解(Hilbert vibration decomposition,HVD)是Feldman等[46]在小波分析和HHT時(shí)頻分析方法的基礎(chǔ)上總結(jié)出的一種非平穩(wěn)信號(hào)分析的時(shí)頻分析方法。該方法簡要步驟如下[47-48]:

    1)對(duì)扭振原始信號(hào)進(jìn)行Hilbert變換,從而獲得其解析信號(hào)。

    2)需要對(duì)原始信號(hào)的非平穩(wěn)頻率部分進(jìn)行處理,可采用低通濾波的方法對(duì)其進(jìn)行分離,從而獲取扭振信號(hào)中最大幅值分量的瞬時(shí)頻率,利用同步檢測的方法取得此信號(hào)分量的初始相位以及最大振幅值。

    3)利用Hilbert迭代運(yùn)算檢測出原始信號(hào)每個(gè)信號(hào)分量的時(shí)頻信息。

    肖璞等[49]對(duì)HVD法進(jìn)行了優(yōu)化,在HVD法三個(gè)步驟采用鏡像延拓的方法對(duì)信號(hào)進(jìn)行端點(diǎn)預(yù)處理,消除了前兩個(gè)步驟所帶來的邊界效應(yīng),并通過信號(hào)仿真和軸承故障診斷實(shí)驗(yàn)驗(yàn)證了該方法的有效性。

    本方法避免了HHT方法中復(fù)雜的EMD過程,提高了分解速度。與此同時(shí)對(duì)多分量非平穩(wěn)信號(hào)進(jìn)行了細(xì)致分析,分解精度也得到提高。

    表2為各種信號(hào)提取分析方法的特點(diǎn)。

    3 展 望

    由于旋轉(zhuǎn)機(jī)械領(lǐng)域扭振事故的頻發(fā),國內(nèi)外相繼提出多種扭振測量方法。隨著無線遙測技術(shù)的飛速發(fā)展,接觸式測量法在信號(hào)降噪方面會(huì)取得一定程度上的進(jìn)步,但無法應(yīng)用于高速旋轉(zhuǎn)軸的檢測仍是限制其發(fā)展的最大障礙。非接觸式測量為當(dāng)今主流技術(shù),而在非接觸扭振測量技術(shù)當(dāng)中,脈沖時(shí)序計(jì)數(shù)法是當(dāng)前應(yīng)用最廣泛,對(duì)應(yīng)產(chǎn)品最多的一種扭振測量技術(shù),但對(duì)于傳統(tǒng)磁電式傳感器而言,存在漏磁大、磁場強(qiáng)度小、測量誤差大等缺點(diǎn)。目前基于磁聚焦技術(shù)在醫(yī)療、無線電力傳輸、金屬表面探傷等領(lǐng)域開始得到應(yīng)用[50-51],其特點(diǎn)是在很小的區(qū)域范圍內(nèi)形成較大的磁場強(qiáng)度,不存在漏磁問題,所以可以考慮利用聚焦磁場原理設(shè)計(jì)傳感器,提高聚磁效率,提高傳感器精度。同時(shí)用柔性電路板(FPC)粘貼在被測軸上取代傳統(tǒng)齒輪盤,解決了齒輪盤不易安裝這一缺點(diǎn)。對(duì)于扭振信號(hào)提取分析方法,基于傅里葉變換發(fā)展起來的方法,由于基函數(shù)的存在,始終存在自適應(yīng)性不足的缺點(diǎn),今后基于EMD信號(hào)分解法將逐漸成為主流技術(shù),EMMD,CEMMD在一定程度上解決了EMD分解法的模態(tài)混疊問題。對(duì)于HHT法來講,仍需要對(duì)EMD分解法曲線擬合方式、邊界效應(yīng)等問題作進(jìn)一步的改善,以得到更加理想的處理效果。

    參 考 文 獻(xiàn):

    [1] KUMENKOA I. On Ways to Reduce the Risks of Operating Energy Industry Facilities[C]// Proc. 6th Seminar School- Evaluation and Control of Industrial Risks in Industrial Safety. Risk Monitoring of Complex and Unique Objects, Aug. 24-25, 2016 Omsk, Russia, 2016: 107.

    [2] 唐忠順.一次風(fēng)機(jī)主軸斷裂的原因分析及處理[J].風(fēng)機(jī)技術(shù),2011(4):73.

    TANG Zhongshun. Reason Analysis and Disposal for Main Shaft Ruputure of Primary Fan[J]. Chinese Journal of Turbomachinery, 2011(4): 73.

    [3] 周國強(qiáng),甘少煒,雷偉,等. 應(yīng)變測量在軸系扭轉(zhuǎn)振動(dòng)測試中的應(yīng)用[J].船海工程,2014(2):165.

    ZHOU Guoqiang, GAN Shaowei, LEI Wei, et al. Application of Shafting Torsional Vibration Measurement Basedon Strain Measuring[J]. Ship & Ocean Engineering, 2014(2):165.

    [4] HAN H S, LEE K H, PARK S H. Estimate of the Fatigue Life of the Propulsion Shaft from Torsional Vibration Measurement and the Linear Damage Summation Law in Ships[J]. Ocean Engineering, 2015, 107(1): 212.

    [5] 楊卓君,廖明夫.基于加速度傳感器測量扭振方法的研究[J].噪聲與振動(dòng)控制,2008(5):163.

    YANG Zhuojun, LIAO Mingfu. Measurement of Torsional Vibration Using Acceleration Sensors[J]. Noise and Vibration Control, 2008(5): 163.

    [6] 魏莉,劉芹,王兢兢,等.光纖光柵應(yīng)變傳感在扭振測量中的應(yīng)用[J].機(jī)械設(shè)計(jì)與制造,2020(4):204.

    WEI Li, LIU Qin, WANG Jingjing, et al. The Application of Fiber Bragg Grating Strain Sensor in the Torsional Vibration Measurement[J]. Machinery Design & Manufacture, 2020(4): 204.

    [7] 張傳敏,胡國良,張濤.光柵傳感器細(xì)分技術(shù)的研究[J].機(jī)械設(shè)計(jì)與制造,2017(9):159.

    ZHANG Chuanmin, HU Guoliang, ZHANG Tao. Study on Subdivision Technology of Grating Sensor[J]. Machinery Design & Manufacture, 2017(9): 159.

    [8] 蔣云帆,廖明夫,王四季.航空發(fā)動(dòng)機(jī)轉(zhuǎn)子扭振測量新方法[J].振動(dòng)測試與診斷,2013, 33(3):410.

    JIANG Yunfan, LIAO Mingfu, WANG Siji. New Measuring Method for Torsional Vibration of Aeroengine Rotor[J]. Journal of Vibration, Measurement & Diagnosis, 2013, 33(3): 410.

    [9] 郭文新,李志深,李富才,等.基于高頻采樣法的軸系扭振測試新方法[J].振動(dòng)測試與診斷,2017,37(6):1082.

    GUO Wenxin, LI Zhishen, LI Fucai, et al. A New Method for Torsional Vibration Testing Based on the High Frequency Sampling[J]. Machinery Design & Manufacture, 2017,37(6): 1082.

    [10]李建科,何國華,蘭涌森,等.風(fēng)力發(fā)電機(jī)組傳動(dòng)系統(tǒng)軸系扭振測量與抑制[J].船舶工程,2019,41(1):252.

    LI Jianke, HE Guohua, LAN Yongsen, et al. Measurement and Suppression of Shafting Torsional Vibration of Wind Turbine Drive System[J]. Ship Engineering, 2019,41(1): 252.

    [11]WU Z Y, BIAN L X, WANG S X, et al.An Angle Sensor Based on Magnetoelectric Effect[J]. Sensors & Actuators A: Physical, 2017,262:108.

    [12]夏永洪,王善銘,黃劭剛,等.齒諧波繞組感應(yīng)電動(dòng)勢波形快速計(jì)算[J].電機(jī)與控制學(xué)報(bào),2011,15(9):1.

    XIA Yonghong, WANG Shanming, HUANG Shaogang, et al. Fast Calculation of the Induced EMF Waveform of the Tooth Harmonic Windings[J]. Electric Machines and Control, 2011, 15(9): 1.

    [13]杜冬.扭振信號(hào)的測試與半實(shí)物仿真[D].上海:上海交通大學(xué),2015.

    [14]王思文.基于光電式脈沖時(shí)序法的扭振測量技術(shù)研究[D].大連:大連理工大學(xué), 2019.

    [15]NAZOLIN A L, POLYAKOY C V I, GNEZDILOV S G. Diagnostics of Shaft Trains of Large Turbine Sets for Integrity Using Torsional Vibration Measurements[J]. Thermal Engineering, 2020, 67(1): 25.

    [16]LISKA J, JAKL J, KUNKEL S. Measurement and Evaluation of Shaft Torsional Vibrations Using Shaft Instantaneous Angular Velocity [J].Journal of Engineering for Gas Turbines and Power, 2019,141(4):041029.1.

    [17]張玉皓,顧煜炯,馬曉騰,等.廣義增量編碼器瞬時(shí)角速度計(jì)算的扭振在線測量[J].儀器儀表學(xué)報(bào),2020,41(10):187.

    ZHANG Yuhao,GU Yujiong,MA Xiaoteng,et al. On Line Torsional Vibration Measurement Based on Generalized Incrementalencoder Instantaneous Angular Speed Calculation[J]. Chinese Journal of Scientific Instrument,2020,41(10):187.

    [18]田中旭,李廣州,王恒宇,等.軸系扭轉(zhuǎn)振動(dòng)測量的雙線逼近法[J].汽車技術(shù),2020(7):40.

    TIAN Zhongxu, LI Guangzhou, WANG Hengyu, et al. Two-line Approximation for Measurement of Torsional Vibration of Shafting[J]. Automobile Technology, 2020(7): 40.

    [19]DIAMOND D H, HEYNS P S, OBERHOLSTER A J. Online Shaft Encoder Geometry Compensation for Arbitrary Shaft Speed Profiles Using Bayesian Regression[J]. Mechanical Systems and Signal Processing,2016,81: 402.

    [20]LI K, LI Y, HAN Y, et al. An EM Induction Hi-Speed Rotation Angular Rate Sensor[J]. Sensors, 2017,17(3): 610.

    [21]FUY C, FAN W, JIN H X, et al. A New Capacitance Angle Sensor of Concentric Ring Multi-layer Differential[J].Measurement,2020,158:107625.

    [22]ONADenis Ijike, TIAN G Y, SUTTHAWEEKUL Ruslee, et,al. Design and Optimisation of Mutual Inductance Based Pulsed Eddy Current Probe[J]. Measurement, 2019, 144: 402.

    [23]黃震,劉彬,董全林.基于激光多普勒技術(shù)扭振測量的研究[J].光學(xué)學(xué)報(bào),2006,26(3):389.

    HUANG Zhen, LIU Bin, DONG Quanlin. Research on the Torsional Vibration Measurement Based on Laser Doppler Technique[J]. Acta Optica Sinica, 2006, 26(3): 389.

    [24]LING X, YANG S, GAN C. Torsional Vibration Measurements on Rotating Shaft System Using Laser Doppler Vibrometer[J]. Optics and Lasers in Engineering, 2012, 50(11): 1596.

    [25]高欣.基于CCD技術(shù)的汽輪發(fā)電機(jī)組扭振測量方法研究[D].哈爾濱:哈爾濱工業(yè)大學(xué),2006.

    [26]馮曉輝.船舶軸系扭振信號(hào)分析及應(yīng)變測量系統(tǒng)開發(fā)[D].武漢:武漢理工大學(xué),2018.

    [27]GAO W C, LIU W, ZOU J X. Damage Detection Methods Based on Changes of Vibration Parameters: A Summary Review[J]. Journal of Vibration and Shock, 2004, 23(4): 1.

    [28]QIAN S, CHEN D. Joint Time-Frequency Analysis[J]. IEEE Journals and Magazings, 1999, 16(2): 52.

    [29]胡殿剛,馬喜平,趙鳳展,等.一種基于短時(shí)傅里葉變換的電壓閃變信號(hào)的檢測[J].電網(wǎng)與清潔能源,2020,36(3):28.

    HU Diangang, MA Xiping, ZHAO Fengzhan, et al.A New Method of Voltage Flicker Detection Based on Short Time FourierTransform[J]. Power System and Clean Energy, 2020,36(3): 28.

    [30]RUZZENE M, FASANA A, GARIBALDI L, et al. Natural Frequencies and Dampings Identification Using Wavelet Transform Application to Real Data[J]. Mechanical Systemsand Signal Processing, 1997, 11(2): 207.

    [31]WANG Z N, KANG Y, DENG L I, et al. Investigating the Hydrodynamics of Airlift Pumps by Wavelet Packet Transformand the Recurrence Plot[J]. Experimental Thermal & Fluid Science, 2018, 92: 56.

    [32]范志峰,張融,李芳環(huán). 基于小波變換的行星輪系點(diǎn)蝕故障信號(hào)分析[J].煤礦機(jī)械,2020,41(9):193.

    FAN Zhifeng, Zhang Rong, Li Fanghuan. Analysis of Pitting Fault Signal of Planetary Gear Train Based on Wavelet Transform[J]. Coal Mine Machinery, 2020, 41(9):193.

    [33]STOCKWELL R G, MANSINHA L, LOWE R P. Localization of the Complex Spectrum[J]. IEEE Transactions on Signal Processing, 1996, 44(4): 998.

    [34]凌啟輝,閆曉強(qiáng),張義方.基于S變換的熱連軋機(jī)耦合振動(dòng)特征提取[J].振動(dòng),測試與診斷,2016,36(1):115.

    LING Qihui, YAN Xiaoqiang, ZHANG Yifang. Coupling Vibration Feature Extraction of Hot Continuous Rolling Mill Based on Adaptive Frequency Domain Filtering and S-Transform[J]. Journal of Vibration, Measurement & Diagnosis, 1996, 44(4): 998.

    [35]PINNERGAR C R, MANSINHA L. The S-transform with Windows of Arbitrary and Varying Shape [J]. Geophysics, 2003, 68(1): 318.

    [36]陳換過,易永余,陳文華,等.基于廣義S變換的齒輪箱軸承故障診斷方法[J].中國機(jī)械工程,2017,28(1):51.

    CHEN Huanguo, YI Yongyu, CHEN Wenhua. et al. Fault Diagnosis Method of Gearbox Bearings Based on Generalized S-transform[J]. China Mechanical Engineering, 2017, 28(1): 51.

    [37]沈薇,陶新民,高珊,等.基于同步擠壓小波變換的振動(dòng)信號(hào)自適應(yīng)降噪方法[J].振動(dòng)與沖擊,2018,37(14):239.

    SHEN Wei, TAO Xinmin, GAO Shan, et al. Self-adaptive De-noising Algorithm for Vibration Signals Based on Synchrosqueezed Wavelet Transforms[J]. Journal of Vibrations and Shock, 2018, 37(14): 239.

    [38]DAUBECHIES I, LU J F, WU H.T. Synchrosqueezed Wavelet Transforms: An Empirical Mode Decomposition Like Tool[J]. Applied and Computational Harmonic Analysis, 2011, 30(2): 243.

    [39]潘高元,呂舜酩,安增輝,等.滾動(dòng)軸承沖擊特征提取的同步擠壓S變換方法[J].振動(dòng)工程學(xué)報(bào),2020,33(2):433.

    PAN Gaoyuan, LV Shunming, AN Zenghui,et al. Impact Feature Extraction from Rolling Bearing Fault Signal by Synchrosqueezed S-transform[J] Journal of Vibration, 2020,33(2):433.

    [40]LI Y, XU M, LIANG X, et al. Application of Bandwidthemd and Adaptive Multiscale Morphology Analysis for Incipient Fault Diagnosis of Rolling Bearings[J]. IEEE Transactions on Industrial Electronics, 2017, 64(8): 6506.

    [41]陳宗祥,陳明星,焦民勝,等.基于改進(jìn)EMD和雙譜分析的電機(jī)軸承故障診斷和實(shí)現(xiàn)[J].電機(jī)與控制學(xué)報(bào),2018,22(5):78.

    CHEN Zongxiang, CHEN Mingxing, JIAO Minsheng, et.al. Fault Diagnosis of Motor Bearing Using Modified Empirical Mode Decomposition and Bi-spectrum[J]. Electric Machines and Control, 2018, 22(5):78.

    [42]張超,陳建軍,郭迅,等.基于EMMD分解的滾動(dòng)軸承故障診斷[J].機(jī)械強(qiáng)度,2012,34(5):650.

    ZHANG Chao, CHEN Jianjun, Guo Xun, et al. Composite Fault Diagnosis For Bearing of Electric Machine Based on EMMD[J].Journal of Mechanical Strength, 2012, 34(5): 650.

    [43]孫曉娟,王利.基于CEEMD小波包算法的降噪方法研究[J].計(jì)算機(jī)與現(xiàn)代化,2020(9):73.

    SUN Xiaojuan, WANG Li. A Denoising Method Based on CEEMD Wavelet Packet[J]. Computer and Modernization, 2020(9):73.

    [44]HUANG N E.Theempirical Mode Decomposition and the Hilbert Spectrum for Nolinear and Nonstationary Time Series Analysis[J]. Proc.R.Soc Lond A, 1998, 454(1971): 903.

    [45]郝東升.希爾伯特黃變換在風(fēng)機(jī)故障診斷的應(yīng)用[J].內(nèi)燃機(jī)與配件,2016,35(7):44.

    HAO Dongsheng.Application of Hilbert-Huang Transform in Fault Diagnosis of Air Blower[J]. Internal Combustion Engine & Parts,2016,35 (7): 44.

    [46]FELDMAN M. Time-Varying Vibration Decomposition Andanalysis Based on the Hilbert Transform[J]. Journal of Soundand Vibration, 2006, 295(315): 518.

    [47]唐貴基,龐彬.基于改進(jìn)的希爾伯特振動(dòng)分解的機(jī)械故障診斷方法研究[J].振動(dòng)與沖擊,2015,34(3) : 167.

    TANG GuiJi, PANG Bin. A Mechanical Fault Diagnosis Method Based on Improved Hilbert Vibration Decomposition[J]. Journal of Vibration and Shock, 2015, 34(3): 167.

    [48]ZHU X X, YUAN Y M, ZHOU P, et al. An Improved Hilbert Vibration Decomposition Method for Analysis of Rotor Fault Signals[J]. Technical Paper, 2017, 39: 4921.

    [49]肖璞,殷立.改進(jìn)的HVD在滾動(dòng)軸承故障診斷中的應(yīng)用[J].控制工程,2020,27(2):264.

    XIAO Pu, YIN Li.Application of Improved HVD in Fault Diagnosis of Rolling Bearings[J]. Control Engineering of China, 2020, 27(2): 264.

    [50]REN W, MOHAMMED S I, WERELEY S, et al. Magnetic Focus Lateral Flow Sensor for Detection of Cervical Cancer Biomarkers[J]. Anal, Chem, 2019, 91: 2876.

    [51]JUNG Haejoon, LEE Byunghun. Optimization of Magnetic Field Focusing and Null Steering for Selective Wireless Power Transfer [J]. IEEE. Transactions on Power Electronics, 2020, 35(5): 4622.

    (編輯:王 萍)

    收稿日期: 2020-10-13

    基金項(xiàng)目: 國家自然科學(xué)基金(52075090).

    作者簡介:

    郭艷玲(1962—),女,教授,博士研究生導(dǎo)師;

    李志鵬(1963—),男,教授,博士研究生導(dǎo)師.

    通信作者:

    邱 楓(1990—),男,碩士研究生,E-mail:15004528946@163.com.

    3673501908247

    猜你喜歡
    檢測方法發(fā)展趨勢
    關(guān)于食品中氟烯草酸的檢測方法的研究
    淺談瀝青路面施工的非均勻性及檢測方法
    宮頸內(nèi)人乳頭瘤病毒的研究進(jìn)展
    試論當(dāng)前國際金融的形勢與發(fā)展趨勢
    國有企業(yè)經(jīng)濟(jì)市場條件下營銷戰(zhàn)略發(fā)展
    論企業(yè)管理模式與企業(yè)管理現(xiàn)代化
    快速原型技術(shù)在機(jī)械制造中的應(yīng)用及發(fā)展趨勢
    鄉(xiāng)鎮(zhèn)配網(wǎng)規(guī)劃及未來發(fā)展趨勢
    電力計(jì)量裝置異常原因及監(jiān)測方法分析
    小兒氨酚黃那敏顆粒有關(guān)物質(zhì)對(duì)氯苯乙酰胺檢測方法的建立
    中國市場(2016年36期)2016-10-19 05:28:41
    日本精品一区二区三区蜜桃| 亚洲一卡2卡3卡4卡5卡精品中文| 国产亚洲欧美精品永久| 欧美绝顶高潮抽搐喷水| 精品人妻1区二区| 美女高潮喷水抽搐中文字幕| 制服诱惑二区| 高潮久久久久久久久久久不卡| 亚洲精品一卡2卡三卡4卡5卡| 欧美av亚洲av综合av国产av| 18禁美女被吸乳视频| 在线观看免费午夜福利视频| 黄片播放在线免费| 一边摸一边做爽爽视频免费| 亚洲,欧美精品.| 久久精品91无色码中文字幕| 婷婷丁香在线五月| 在线观看舔阴道视频| 亚洲 国产 在线| 长腿黑丝高跟| 后天国语完整版免费观看| 亚洲第一电影网av| 亚洲 欧美 日韩 在线 免费| 中文字幕人妻丝袜一区二区| 精品熟女少妇八av免费久了| 亚洲av成人一区二区三| 涩涩av久久男人的天堂| 51午夜福利影视在线观看| 夜夜夜夜夜久久久久| 久久国产亚洲av麻豆专区| 身体一侧抽搐| 欧美日韩中文字幕国产精品一区二区三区 | 99精品在免费线老司机午夜| 欧美一区二区精品小视频在线| 国产高清激情床上av| 亚洲精品中文字幕一二三四区| 欧美丝袜亚洲另类 | 色在线成人网| 黄色 视频免费看| 中出人妻视频一区二区| 精品一区二区三区四区五区乱码| 禁无遮挡网站| 国产欧美日韩一区二区精品| 国产国语露脸激情在线看| 啦啦啦韩国在线观看视频| 少妇 在线观看| 久久亚洲精品不卡| 亚洲国产精品成人综合色| 亚洲欧美精品综合久久99| 女人爽到高潮嗷嗷叫在线视频| www.自偷自拍.com| av在线天堂中文字幕| 亚洲人成电影免费在线| 脱女人内裤的视频| 日韩免费av在线播放| 欧美日本中文国产一区发布| 热99re8久久精品国产| 麻豆一二三区av精品| 色哟哟哟哟哟哟| 亚洲成av片中文字幕在线观看| 午夜福利视频1000在线观看 | 亚洲国产看品久久| 久久人妻熟女aⅴ| 在线av久久热| 亚洲av五月六月丁香网| 国产精品,欧美在线| 1024视频免费在线观看| 黑人欧美特级aaaaaa片| 亚洲一区二区三区色噜噜| 88av欧美| 真人一进一出gif抽搐免费| 日韩大尺度精品在线看网址 | 免费在线观看亚洲国产| 久久精品亚洲精品国产色婷小说| 亚洲片人在线观看| 成年人黄色毛片网站| 亚洲情色 制服丝袜| 女同久久另类99精品国产91| 午夜福利,免费看| 一区二区三区精品91| 黄片播放在线免费| 乱人伦中国视频| 大陆偷拍与自拍| av片东京热男人的天堂| 女性被躁到高潮视频| 男人舔女人下体高潮全视频| 一级,二级,三级黄色视频| 色尼玛亚洲综合影院| 国产精品98久久久久久宅男小说| 国产精品二区激情视频| 国产极品粉嫩免费观看在线| 无限看片的www在线观看| 久久国产精品男人的天堂亚洲| 国产精品av久久久久免费| av片东京热男人的天堂| 亚洲最大成人中文| 国产一卡二卡三卡精品| 在线观看免费日韩欧美大片| 91在线观看av| 国产精品久久久av美女十八| 大型av网站在线播放| 精品久久久久久成人av| 国产午夜福利久久久久久| 性少妇av在线| 一区二区三区国产精品乱码| 一级a爱片免费观看的视频| 久久人妻熟女aⅴ| 人妻久久中文字幕网| 美女国产高潮福利片在线看| 99热只有精品国产| 国产国语露脸激情在线看| 香蕉久久夜色| 99re在线观看精品视频| 亚洲成av人片免费观看| 国产亚洲精品av在线| 久久人人精品亚洲av| 免费女性裸体啪啪无遮挡网站| 色综合欧美亚洲国产小说| 亚洲一区中文字幕在线| 欧美老熟妇乱子伦牲交| 国产极品粉嫩免费观看在线| 色老头精品视频在线观看| 国产av一区二区精品久久| 亚洲五月婷婷丁香| 韩国av一区二区三区四区| 巨乳人妻的诱惑在线观看| 久久中文字幕一级| 国产人伦9x9x在线观看| 欧美+亚洲+日韩+国产| 大陆偷拍与自拍| 美女扒开内裤让男人捅视频| 丝袜在线中文字幕| 亚洲无线在线观看| 国产精品久久电影中文字幕| 午夜精品国产一区二区电影| 狠狠狠狠99中文字幕| 黄色丝袜av网址大全| 久久久久久免费高清国产稀缺| 黄色视频不卡| 欧美成人免费av一区二区三区| 亚洲一区中文字幕在线| 日本精品一区二区三区蜜桃| 久久精品国产99精品国产亚洲性色 | 曰老女人黄片| 日本 av在线| 嫁个100分男人电影在线观看| 看免费av毛片| 久久亚洲精品不卡| 国产麻豆成人av免费视频| 国产伦一二天堂av在线观看| 亚洲成人国产一区在线观看| bbb黄色大片| 又黄又粗又硬又大视频| 国产成年人精品一区二区| 一进一出抽搐动态| 国产乱人伦免费视频| 深夜精品福利| 久久草成人影院| 变态另类成人亚洲欧美熟女 | 久久精品亚洲精品国产色婷小说| 欧美老熟妇乱子伦牲交| 在线av久久热| 麻豆久久精品国产亚洲av| 黑丝袜美女国产一区| 亚洲avbb在线观看| 国产99白浆流出| 国产一区二区三区在线臀色熟女| 欧美大码av| 999久久久国产精品视频| 黑丝袜美女国产一区| 国产成人啪精品午夜网站| 精品国产超薄肉色丝袜足j| 男人的好看免费观看在线视频 | 亚洲第一电影网av| 精品电影一区二区在线| 日本a在线网址| 国产片内射在线| 一进一出抽搐gif免费好疼| 天天添夜夜摸| 日韩欧美国产在线观看| 人人妻人人澡欧美一区二区 | 久久久久久免费高清国产稀缺| 亚洲av电影在线进入| 一卡2卡三卡四卡精品乱码亚洲| 久久精品91蜜桃| 90打野战视频偷拍视频| 成人18禁在线播放| 成人国产一区最新在线观看| cao死你这个sao货| 精品高清国产在线一区| 免费av毛片视频| 欧美日韩一级在线毛片| 亚洲五月色婷婷综合| 色综合站精品国产| 欧美日韩亚洲国产一区二区在线观看| 日韩免费av在线播放| 操美女的视频在线观看| 人妻丰满熟妇av一区二区三区| 亚洲五月色婷婷综合| 男女做爰动态图高潮gif福利片 | av电影中文网址| 久久精品国产亚洲av香蕉五月| 制服诱惑二区| 曰老女人黄片| 90打野战视频偷拍视频| xxx96com| 国产视频一区二区在线看| 两个人免费观看高清视频| 亚洲无线在线观看| 亚洲精品中文字幕在线视频| 天堂动漫精品| 国产成人一区二区三区免费视频网站| 日韩大尺度精品在线看网址 | 亚洲精品美女久久久久99蜜臀| 十八禁人妻一区二区| 精品久久久久久久久久免费视频| 国产精品98久久久久久宅男小说| 亚洲美女黄片视频| 精品国产超薄肉色丝袜足j| www.自偷自拍.com| 国产精品一区二区三区四区久久 | 日韩免费av在线播放| 久久久久久久午夜电影| 在线观看午夜福利视频| 亚洲天堂国产精品一区在线| 欧美成人午夜精品| 午夜a级毛片| 午夜亚洲福利在线播放| 国产极品粉嫩免费观看在线| 午夜免费鲁丝| 亚洲三区欧美一区| 国产不卡一卡二| 亚洲男人天堂网一区| 免费看美女性在线毛片视频| 69精品国产乱码久久久| 久久热在线av| 精品国产超薄肉色丝袜足j| 欧美黑人精品巨大| 母亲3免费完整高清在线观看| 午夜久久久在线观看| 超碰成人久久| 午夜福利欧美成人| 免费人成视频x8x8入口观看| 日本三级黄在线观看| 一边摸一边抽搐一进一出视频| avwww免费| 久久精品91无色码中文字幕| av天堂久久9| 麻豆成人av在线观看| 国产xxxxx性猛交| 一个人免费在线观看的高清视频| 12—13女人毛片做爰片一| 视频区欧美日本亚洲| 成年人黄色毛片网站| 久久中文字幕人妻熟女| 亚洲中文av在线| 一区二区三区激情视频| 很黄的视频免费| 久久久精品国产亚洲av高清涩受| 高清毛片免费观看视频网站| 青草久久国产| 亚洲国产精品久久男人天堂| 国产精品国产高清国产av| 日韩免费av在线播放| avwww免费| 国内毛片毛片毛片毛片毛片| 精品欧美国产一区二区三| 亚洲免费av在线视频| 国产成人影院久久av| 中亚洲国语对白在线视频| 成人免费观看视频高清| 国产欧美日韩一区二区三| 欧美一区二区精品小视频在线| 亚洲专区中文字幕在线| 国产精品一区二区三区四区久久 | 午夜福利成人在线免费观看| 自线自在国产av| 免费搜索国产男女视频| 男男h啪啪无遮挡| 亚洲国产精品sss在线观看| 在线观看午夜福利视频| 亚洲午夜理论影院| 亚洲第一青青草原| 久久草成人影院| 亚洲欧美日韩另类电影网站| 黄色视频,在线免费观看| 在线播放国产精品三级| 91成人精品电影| 日韩大码丰满熟妇| 无限看片的www在线观看| 国产免费av片在线观看野外av| 国产精品秋霞免费鲁丝片| 在线观看66精品国产| 亚洲人成77777在线视频| √禁漫天堂资源中文www| 身体一侧抽搐| 亚洲黑人精品在线| 亚洲av片天天在线观看| 亚洲男人天堂网一区| 成人国语在线视频| 精品不卡国产一区二区三区| 久久中文字幕一级| 人妻丰满熟妇av一区二区三区| 91成人精品电影| 亚洲精品久久国产高清桃花| 国产精品亚洲美女久久久| 1024香蕉在线观看| 成人三级做爰电影| 精品人妻1区二区| 久久久久九九精品影院| 免费在线观看视频国产中文字幕亚洲| 91精品三级在线观看| 成人国语在线视频| 成人av一区二区三区在线看| 亚洲七黄色美女视频| 久久久久久久久免费视频了| 精品一区二区三区av网在线观看| 国产人伦9x9x在线观看| 亚洲国产欧美网| 久久影院123| 黄频高清免费视频| 狠狠狠狠99中文字幕| 亚洲国产欧美一区二区综合| 在线av久久热| 两性午夜刺激爽爽歪歪视频在线观看 | 大型av网站在线播放| 免费人成视频x8x8入口观看| 母亲3免费完整高清在线观看| 视频区欧美日本亚洲| 亚洲人成77777在线视频| 人成视频在线观看免费观看| 窝窝影院91人妻| 国产主播在线观看一区二区| 最新在线观看一区二区三区| 久久久久久人人人人人| 两性夫妻黄色片| 色婷婷久久久亚洲欧美| 日韩视频一区二区在线观看| 日韩大尺度精品在线看网址 | e午夜精品久久久久久久| 国产在线精品亚洲第一网站| 女人被狂操c到高潮| 久久久久久久午夜电影| 91av网站免费观看| 一级作爱视频免费观看| 中文字幕最新亚洲高清| 中文字幕最新亚洲高清| 国产欧美日韩一区二区三区在线| videosex国产| 操美女的视频在线观看| 在线十欧美十亚洲十日本专区| 亚洲电影在线观看av| 乱人伦中国视频| 国产av一区二区精品久久| 90打野战视频偷拍视频| 久久性视频一级片| 久久久久国产一级毛片高清牌| 精品久久久久久久人妻蜜臀av | 97人妻天天添夜夜摸| 美女免费视频网站| 99国产精品免费福利视频| 91av网站免费观看| 日韩有码中文字幕| 国产成人av激情在线播放| 国产免费男女视频| 好男人电影高清在线观看| 啦啦啦免费观看视频1| 天天添夜夜摸| 午夜福利成人在线免费观看| 少妇 在线观看| av视频免费观看在线观看| 亚洲性夜色夜夜综合| а√天堂www在线а√下载| 欧美色欧美亚洲另类二区 | АⅤ资源中文在线天堂| 亚洲aⅴ乱码一区二区在线播放 | 丁香六月欧美| 少妇 在线观看| 宅男免费午夜| 一区在线观看完整版| 日本vs欧美在线观看视频| 久久人人爽av亚洲精品天堂| 午夜日韩欧美国产| 9热在线视频观看99| 午夜影院日韩av| 99精品欧美一区二区三区四区| 国产一区二区激情短视频| 国产熟女午夜一区二区三区| 免费观看精品视频网站| 国产精品1区2区在线观看.| 黄色女人牲交| 1024视频免费在线观看| 国产高清videossex| 欧美日韩精品网址| 在线视频色国产色| 男女下面插进去视频免费观看| 性少妇av在线| 久久热在线av| 国产视频一区二区在线看| 国产精品二区激情视频| 国产精品 国内视频| 99久久精品国产亚洲精品| 一级毛片精品| 日韩成人在线观看一区二区三区| 色综合站精品国产| 麻豆国产av国片精品| 亚洲第一电影网av| 色婷婷久久久亚洲欧美| 亚洲熟妇中文字幕五十中出| 成人国语在线视频| 曰老女人黄片| 搡老岳熟女国产| 久久午夜亚洲精品久久| 国产成+人综合+亚洲专区| 国产又色又爽无遮挡免费看| 一边摸一边抽搐一进一小说| 亚洲在线自拍视频| 国产精品久久久av美女十八| 首页视频小说图片口味搜索| 日本五十路高清| 在线观看免费午夜福利视频| 日本免费一区二区三区高清不卡 | 波多野结衣高清无吗| 热99re8久久精品国产| 婷婷精品国产亚洲av在线| 极品人妻少妇av视频| 嫩草影视91久久| 男人操女人黄网站| 亚洲国产欧美日韩在线播放| 久久香蕉激情| 黄色 视频免费看| 久久午夜亚洲精品久久| 老汉色∧v一级毛片| 国产熟女午夜一区二区三区| 91字幕亚洲| 久9热在线精品视频| 亚洲一区二区三区色噜噜| 好男人在线观看高清免费视频 | 亚洲男人天堂网一区| 精品久久久久久久久久免费视频| 人人妻人人爽人人添夜夜欢视频| 怎么达到女性高潮| 国产黄a三级三级三级人| 久久青草综合色| 亚洲精品美女久久久久99蜜臀| 一级毛片女人18水好多| 19禁男女啪啪无遮挡网站| 国产99久久九九免费精品| 欧美激情极品国产一区二区三区| 在线观看免费日韩欧美大片| 一个人免费在线观看的高清视频| 两个人看的免费小视频| 婷婷六月久久综合丁香| 免费在线观看黄色视频的| 变态另类成人亚洲欧美熟女 | 无人区码免费观看不卡| 女人高潮潮喷娇喘18禁视频| 美女 人体艺术 gogo| 亚洲美女黄片视频| 国产片内射在线| 日本精品一区二区三区蜜桃| 色综合婷婷激情| 久久人妻熟女aⅴ| 成人永久免费在线观看视频| 啪啪无遮挡十八禁网站| 国产99白浆流出| 精品一区二区三区四区五区乱码| 美女免费视频网站| 亚洲第一电影网av| 欧美日本中文国产一区发布| 黑人操中国人逼视频| 亚洲性夜色夜夜综合| 欧美国产精品va在线观看不卡| 中亚洲国语对白在线视频| 一区二区三区激情视频| 女警被强在线播放| 亚洲自拍偷在线| 精品久久蜜臀av无| av在线播放免费不卡| 人人澡人人妻人| 天天添夜夜摸| 久热爱精品视频在线9| 69精品国产乱码久久久| 欧美日本亚洲视频在线播放| √禁漫天堂资源中文www| 宅男免费午夜| 男男h啪啪无遮挡| 精品久久久久久久毛片微露脸| 天天一区二区日本电影三级 | 丝袜在线中文字幕| 女人精品久久久久毛片| 91麻豆精品激情在线观看国产| 国产成人啪精品午夜网站| 国产蜜桃级精品一区二区三区| 欧美久久黑人一区二区| avwww免费| 国产亚洲欧美98| 男女午夜视频在线观看| 最新在线观看一区二区三区| 国产成人免费无遮挡视频| 90打野战视频偷拍视频| av免费在线观看网站| 亚洲黑人精品在线| 麻豆久久精品国产亚洲av| 又大又爽又粗| 国产高清视频在线播放一区| 色精品久久人妻99蜜桃| 精品电影一区二区在线| 丝袜人妻中文字幕| 欧美+亚洲+日韩+国产| 国产亚洲av嫩草精品影院| 看免费av毛片| 久久精品aⅴ一区二区三区四区| 国产亚洲精品一区二区www| 美女高潮到喷水免费观看| 欧美国产日韩亚洲一区| 无遮挡黄片免费观看| √禁漫天堂资源中文www| 在线观看免费视频网站a站| 中文字幕最新亚洲高清| 亚洲第一青青草原| 亚洲一区二区三区不卡视频| 淫秽高清视频在线观看| 欧美精品啪啪一区二区三区| 国内毛片毛片毛片毛片毛片| 日本黄色视频三级网站网址| 视频在线观看一区二区三区| 午夜两性在线视频| 久久久久国内视频| 波多野结衣高清无吗| 国产极品粉嫩免费观看在线| 一a级毛片在线观看| 久久久久国产一级毛片高清牌| 欧美黑人欧美精品刺激| 久久影院123| 欧美日韩亚洲综合一区二区三区_| 成人精品一区二区免费| 麻豆av在线久日| 亚洲成人免费电影在线观看| 午夜老司机福利片| 国产欧美日韩一区二区三| 亚洲成a人片在线一区二区| 两个人免费观看高清视频| 日日爽夜夜爽网站| 亚洲精品粉嫩美女一区| 91在线观看av| 成人手机av| 亚洲中文字幕日韩| e午夜精品久久久久久久| 黄片小视频在线播放| 欧洲精品卡2卡3卡4卡5卡区| 久久久久久国产a免费观看| 久久精品亚洲精品国产色婷小说| 精品国产乱码久久久久久男人| 国产午夜福利久久久久久| 亚洲免费av在线视频| 国产亚洲精品第一综合不卡| 中文字幕色久视频| 波多野结衣av一区二区av| 一卡2卡三卡四卡精品乱码亚洲| 久久香蕉精品热| 亚洲精品中文字幕一二三四区| 777久久人妻少妇嫩草av网站| 精品国产一区二区久久| 国产一区在线观看成人免费| 久久久久精品国产欧美久久久| av在线天堂中文字幕| 亚洲欧美激情在线| 黄片大片在线免费观看| 精品久久久久久,| 亚洲国产中文字幕在线视频| 国产精品久久久久久精品电影 | 一区二区三区国产精品乱码| 国产精品久久久久久精品电影 | 久久热在线av| 国产高清激情床上av| 欧美成人免费av一区二区三区| av有码第一页| 精品卡一卡二卡四卡免费| 久久久久九九精品影院| 精品国产超薄肉色丝袜足j| 久久精品国产99精品国产亚洲性色 | 成人手机av| 97人妻天天添夜夜摸| 国产精品久久视频播放| 村上凉子中文字幕在线| 久久久久久亚洲精品国产蜜桃av| 亚洲男人天堂网一区| 国内精品久久久久精免费| 久久久久久久久久久久大奶| 日本五十路高清| 女人被狂操c到高潮| 黄色视频不卡| 热re99久久国产66热| 这个男人来自地球电影免费观看| 精品不卡国产一区二区三区| 亚洲免费av在线视频| 午夜影院日韩av| 日日夜夜操网爽| 国产一区在线观看成人免费| 脱女人内裤的视频| www.999成人在线观看| 免费观看精品视频网站| 一级,二级,三级黄色视频| 又紧又爽又黄一区二区| 老司机午夜福利在线观看视频| 久久精品国产亚洲av香蕉五月| 韩国av一区二区三区四区| 午夜亚洲福利在线播放| 免费在线观看影片大全网站| 麻豆av在线久日| 成人三级黄色视频| 一级毛片精品| 国产精品自产拍在线观看55亚洲| 欧美一级a爱片免费观看看 |