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

    基于圓域分析的大型回轉支承初期故障診斷

    2017-05-17 02:35:34黃筱調洪榮晶
    振動與沖擊 2017年9期
    關鍵詞:時域橢圓故障診斷

    封 楊, 黃筱調, 洪榮晶, 陳 捷

    (南京工業(yè)大學 機械與動力工程學院, 南京 210009)

    基于圓域分析的大型回轉支承初期故障診斷

    封 楊, 黃筱調, 洪榮晶, 陳 捷

    (南京工業(yè)大學 機械與動力工程學院, 南京 210009)

    為解決真實工況下大型回轉支承振動信號背景噪聲大、常用故障診斷方法難以適用的問題,提出了一種基于圓域分析的振動信號處理方法。將時域信號進行圓域轉換,并按一定角度將轉換后的圓域信號劃分成多個區(qū)域;判斷各區(qū)域信號鄰域相關離散點擬合橢圓的傾角方向,得到回轉支承整圈對應的多個異常向量;以異常向量的平均向量作為圓域分析的特征向量,分析其均值、方差、歪度和峭度指標的變化情況,實現(xiàn)回轉支承的故障診斷。對某型號回轉支承進行了加速壽命試驗,結果表明,該方法能夠有效診斷出回轉支承滾道的區(qū)域滑移、點蝕等初期故障,相比常見的時域特征、小波分析等方法準確度更高,故障可識別度更強,因此可以用于實際工況下回轉支承的故障診斷。

    回轉支承;圓域分析;圓域重采樣;故障診斷;加速壽命試驗

    大型回轉支承作為重要的回轉連接件被廣泛使用在風力發(fā)電機、港口機械、建筑工程機械等大重型機械中,其滾動體中心圓直徑通常超過600 mm,內圈或外圈會加工出整圈齒輪用于驅動。國家標準對回轉支承的可靠性提出了很高的要求[1],但近年來由于回轉支承的失效導致的機器停機事故屢見不鮮,而回轉支承的制造、運輸、吊裝的周期很長且花費巨大,風機一次的吊裝成本可達數(shù)十萬,港口機械的停機帶來的損失更是難以估量。因此,大型回轉支承的故障診斷尤其是初期故障的判別尤為重要。

    軸承的典型性能退化過程一般如下:初始會有一個磨合過程,然后在較長時間段內運行相對平穩(wěn),到中后期軸承的局部滾道產生裂紋、點蝕等初始故障后,性能很快衰退直至完全失效[2]。大量的專著和文獻對轉速>600 r/min、外圈直徑<500 mm的中高速軸承的故障診斷方法進行了深入地研究和探討,此類軸承故障時的振動信號沖擊成分明顯,信噪比高,因此各種方法在試驗或應用中都得到了有效的驗證。但是,大型回轉支承體積龐大,轉速通常<10 r/min,有時可低至1 r/min,滾道與滾動體接觸時產生的沖擊成分非常微弱,振動信號中有效信息被高頻噪聲完全淹沒,這使得同樣的診斷方法應用在大型低速的回轉支承振動信號中效果并不理想[3],給其故障診斷帶來了很大的挑戰(zhàn)[4]。

    盡管如此,大型回轉支承的故障診斷近年來依然受到了大量學者的關注。Zvokelj等[5-6]先后提出集合經驗模態(tài)分解-多分辨主成分分析(Ensemble Empirical Mode Decomposition-Multi-Scale Principal Component Analysis, EEMD-MSPCA)和集合經驗模態(tài)分解-多分辨核主元分析(EEMD-MSKPCA)的降噪方法用于回轉支承振動信號,并利用PCA/KPCA(Principal Component Analysis/Kernel PCA)模型準確診斷出回轉支承的滾道故障。郭剛濤等[7]將混沌理論用于大口徑天線中的回轉支承故障診斷,利用Lyapunov指數(shù)識別出回轉支承的滾道點蝕故障。趙陽等則采用了經驗模態(tài)分解(Empirical Mode Decomposition,EMD)和獨立分量分析(Independent Component Analysis,ICA)相結合的盲源分離方法,從盲源信號的功率譜中分析出回轉支承故障特征頻率。楊杰等提出了小波-能量模式的回轉支承診斷方法,從加速壽命試驗的振動信號中分離出了故障頻率。然而,文獻[8]中回轉支承的滾道故障均是人工加工出的嚴重缺陷,并不符合回轉支承自然性能退化過程中產生的諸如點蝕、裂紋、剝落等故障,而文獻[9]中的加速壽命試驗超載了90倍,這在工程實際中也較為罕見。此外,上述文獻中的振動信號中均能明顯觀察到高頻沖擊成分,與回轉支承真實工況下振動信號的特征也不甚相符。

    由以上可知,現(xiàn)有的回轉支承故障診斷方法多是針對振動信號的時域、頻域或者時頻域進行分析,此類方法對嚴重故障或是壽命試驗中后期的振動信號的處理非常有效,但是仍然很難用于實際工況下回轉支承的故障尤其是初期故障的診斷。相比之下,振動信號的一些特征在圓域內與在時域、頻域內有著截然不同的呈現(xiàn)方式,圓域采樣及相關的圓域分析近年來越來越多地被用于軸承[10-11]、風機[12-13]、齒輪箱[14]等故障診斷中,而Caesarendra等[15]將其應用于回轉支承狀態(tài)監(jiān)測,取得了較為理想的效果。本文將進一步探討振動信號的圓域分析方法,并將其用于大型回轉支承的初期故障診斷。

    1 背景知識

    1.1 PAA與鄰域相關圖

    軸承振動信號的采樣頻率很高,即使是轉速極低的回轉支承,其采樣率通常>1 kHz,如此龐大的數(shù)據(jù)量將顯著影響故障診斷模型的計算效率,使其很難應用到實時狀態(tài)監(jiān)測系統(tǒng)中。因此,有必要對原始振動信號的數(shù)據(jù)量進行縮減。分段累積近似法(Piecewise Aggregate Approximation, PAA)作為一種數(shù)據(jù)縮減方法,最初由Yi等[16]提出,其特點是在對大量時域數(shù)據(jù)進行縮減的同時盡可能多地保持數(shù)據(jù)的原有特征。PAA實現(xiàn)較為簡單,對于一組有L個數(shù)據(jù)樣本的序列y=(y1,y2,…,yL),首先定義一個常數(shù)w,w應為L的約數(shù),然后將樣本序列y等分成N段,N=L/w,最后將每段序列的代數(shù)平均值連起來便得到縮減后的新序列Y=(Y1,Y2,…,YN),其中每個元素的計算公式為

    (1)

    式中,w又被稱為PAA窗的大小,顯然w越大,縮減后的數(shù)據(jù)量就越小,但是丟失的信息也就越多。

    除了縮減數(shù)據(jù)量外,PAA另一個重要的應用是檢測振動信號中頻率的變化,為此需要引入鄰域相關的概念。鄰域相關即是將數(shù)據(jù)縮減后的時域序列Y以鄰域相關圖的方式呈現(xiàn),圖中每個坐標點的橫坐標為Yn,縱坐標為Yn+1,n∈[1,N-1]。將時域振動信號先進行PAA,然后利用直接橢圓擬合法[17]對鄰域相關離散點進行橢圓擬合的過程稱為PAA過程。為研究其檢測信號中頻率變化的有效性,產生兩組幅值為1、頻率分別為1 Hz和4 Hz的標準正弦波形,其采樣率fs均為48 Hz,采樣時間均為1 s,令w=4,則兩組波形的PAA及鄰域相關圖,如圖1所示。

    (a) 1 Hz正弦信號

    (b) 4 Hz正弦信號

    圖1 兩組仿真信號的PAA及鄰域相關圖

    Fig.1 PAA and neighborhood correlation of two simulated signals

    由圖1可知,正弦信號的鄰域相關圖可擬合出近乎完美的橢圓,但在PAA窗大小w和采樣率fs相同的情況下,不同頻率的正弦波擬合出的橢圓具有不同的傾角,定義圖1(a)中傾角為0~π/2的橢圓為右傾斜橢圓,圖1(b)中傾角為π/2~π的橢圓為左傾斜橢圓。PAA過程正是通過橢圓的傾角判斷出信號中頻率的變化。實際上,回轉支承早期振動信號中的有效成分幾乎完全被背景噪聲淹沒,高頻白噪聲占主導,而隨著故障的產生和加劇,周期性地沖擊信號會越發(fā)明顯,信噪比會提高,此時信號中的中低頻成分占主導。因此,可將PAA過程用于回轉支承的狀態(tài)監(jiān)測,通過擬合橢圓傾角的變化診斷出回轉支承潛在的早期故障。

    1.2 圓域重采樣

    圓域(又稱角度域)分析最初由Luo等[18]提出,用以解決轉動部件變轉速工況下振動信號分析頻譜泄漏的問題。圓域分析的核心思想是將時域振動信號映射到角度域,消除振動信號中的時間屬性,進而消除轉速對振動信號分析的影響。獲取圓域振動信號的方法有兩種:階次追蹤法[19]和圓域重采樣。階次追蹤法通常采用專用儀器直接進行圓域采樣,而圓域重采樣則是通過坐標變換將時域信號轉換成圓域信號。

    假設回轉支承轉動β(0°<β≤360°) 角度需要Ts,時間序列為t=[1,2,…,T],則其對應圓域內的弧度序列為

    (2)

    將轉換后的弧度及其對應的振動幅值在三維空間展開,便可得到時域信號的圓域表達。以圖1(b)中4 Hz的正弦信號為例,若其對應的回轉支承轉動的角度是360°,則圓域重采樣后的信號,如圖2所示。

    (a) 時域信號

    (b) 圓域信號

    由前文可知,PAA過程能夠有效識別振動信號中主要頻率成分的變化,并可用于回轉支承的故障診斷,但實際應用中PAA過程所處理的數(shù)據(jù)量(即序列y的長度)及其對應的物理意義并未給出。在狀態(tài)監(jiān)測系統(tǒng)中,一次PAA過程處理數(shù)據(jù)量越少,總數(shù)據(jù)量一定的情況下需要執(zhí)行PAA過程的次數(shù)就越多,系統(tǒng)計算負荷會大幅提升,而一次PAA過程處理數(shù)據(jù)量越大,信號中的頻率成分就越復雜,擬合橢圓對信號中頻率變化的敏感度會有所降低。因此,確定合適的一次PAA過程處理的數(shù)據(jù)量及其對應的物理意義非常重要。本文將利用圓域重采樣將時域信號轉換成圓域信號,忽略轉速對信號處理的影響,同時將圓域信號按一定角度均分,并對每段信號執(zhí)行PAA過程,最終通過相應的圓域分析實現(xiàn)故障診斷。

    2 圓域分析故障診斷方法

    2.1 故障診斷流程

    本文針對大型回轉支承提出的基于圓域分析的故障診斷流程,如圖3所示,對于一段時域振動信號,首先通過圓域重采樣的方法將其轉換成圓域信號,同時按照一定角度將圓域信號均分成多個區(qū)域,并對每個區(qū)域的信號執(zhí)行PAA過程,然后通過鄰域相關圖中擬合橢圓的傾角方向確定各區(qū)域是否出現(xiàn)異常,最后計算整段圓域信號的相關特征,并監(jiān)測其在時域上的變化趨勢,進而實現(xiàn)回轉支承初期故障診斷。

    圖3 故障診斷流程圖

    具體而言,若振動信號特征頻率為fc、采樣率為fs,回轉支承轉速為qr/min,則回轉支承轉動一周對應著60/qs的振動信號,每秒鐘的振動信號對應圓域內的角度為(6·q)°。因此,可將0~60/qs的時域振動信號轉換成0°~360°的圓域信號,同時每隔(6·q)°將其分割成60/q個角度區(qū)域(每個區(qū)域對應1 s的振動信號),然后對每個區(qū)域的信號執(zhí)行PAA過程。然而,由前文可知,PAA過程中窗大小w與采樣率fs、信號特征頻率fc之間存在某種聯(lián)系,w選取是否合適直接決定了PAA過程識別異常信號的能力,其詳細的優(yōu)化方法將在“2.2”節(jié)中給出。此外,由圖1可知,在同一PAA過程下,高頻信號的擬合橢圓是左傾斜的,當信號頻率降低時,橢圓轉換成右傾斜的,考慮到回轉支承初期振動信號中主要是高頻白噪聲,對應的擬合橢圓應當是左傾斜的,而到中后期振動信號中的中低頻率能量較高,其擬合橢圓應當是右傾斜的。據(jù)此,將PAA過程中擬合橢圓為右傾斜的角度區(qū)域定義為異常區(qū)域并標記,同時作出對應的異常向量,如圖4所示。

    圖4 圓域分析中的異常向量

    圖4將回轉支承一周劃分成8個區(qū)域,對應的轉速q=7.5 r/min,在90°和315°的區(qū)域檢測出異常,作出異常向量V1和V2,然后計算所有異常向量的平均向量得到圓域分析的特征向量Vmean。最后,在狀態(tài)監(jiān)測系統(tǒng)中,回轉支承每轉過一周,可實時計算其振動信號的圓域特征向量Vmean,并求得Vmean在圓域內的均值、方差、歪度和峭度指標作為故障診斷的依據(jù),具體計算過程將在“2.3”節(jié)中給出。

    2.2 PAA窗大小w的優(yōu)化

    由“1.1”節(jié)可知,執(zhí)行PAA過程首先需要確定合適的PAA窗大小w,在同一采樣率fs下,不同的w取值會導致同一段信號擬合出不同傾角方向的橢圓。另一方面,w取定后,不同特征頻率fc的信號也會擬合出不同傾角的橢圓。因此,必然存在某個因變量λ隨著w、fc和fs而改變,控制著擬合橢圓的傾角,當λ達到臨界時橢圓從一個傾角方向轉變?yōu)榱硪粋€方向,其經驗計算公式可定義為

    (3)

    設定fs=2 048 Hz,將fc在0~300 Hz變化,w在2~8取值并對頻率為fc的正弦信號執(zhí)行PAA過程,確定擬合橢圓的傾角方向,最終得到的關系模型如圖5所示。

    由圖5中可知,λ與fc呈線性關系,當λ>0.25時,PAA過程擬合出的橢圓是左傾斜的,當λ<0.25時,擬合橢圓是右傾斜的,當λ=0.25時,PAA過程擬合出的實際上是沒有傾角的圓。因此,λ又可稱為轉換系數(shù),0.25是其閾值。另一方面,當w確定后,每條曲線在λ閾值處都有一個交點,這個點的頻率被稱為轉換頻率,即當信號的頻率高于此頻率時,橢圓是左傾斜的,當信號的頻率低于此頻率時,橢圓是右傾斜的,這與圖1中觀察到的現(xiàn)象是一致的。

    圖5 轉換系數(shù)λ與轉換頻率fc的關系

    如前所述,回轉支承運行初期的振動信號中多是高頻噪聲,此時擬合橢圓是左傾斜的,到中后期振動信號中出現(xiàn)異常時中低頻占主要成分,此時擬合橢圓會出現(xiàn)右傾斜的。因此,工程實際中可先確定信號出現(xiàn)異常時大致的頻率范圍,再由此確定最接近的轉換頻率,最終確定最優(yōu)的w取值。

    2.3 圓域特征計算

    假設信號異常識別步驟中作出的異常向量為V1,V2,…,VM,各向量在圓域內的角度為θ1,θ2,…,θM,則各異常向量可用復數(shù)表達為

    Zm=r·(cosθm+isinθm)=r·eiθm

    (4)

    式中:m∈[1,M];r為各異常向量的模長,此處均為單位1。由此可知,圓域分析的特征向量Vmean的表達式為

    (5)

    其所在圓域的角度θmean應為

    θmean=Angle(Zmean)

    (6)

    據(jù)此,圓域分析的特征向量Vmean的均值、方差、歪度和峭度指標可分別由式(7)~式(10)求得。

    (7)

    (8)

    與輕質薄板相比,ALC板、AS裝配式墻板的保溫、隔熱、防火性能較好。對比ALC板與AS裝配式墻板,ALC板更多地是替代了傳統(tǒng)建筑維護墻體,其缺點是強度低和本身不具備裝飾性。AS裝配式墻板的各項性能均非常優(yōu)越,外墻板為混凝土清水色飾面,質感厚重大氣,可作為墻體和外飾面使用,不需要做額外的外裝飾。由于AS墻板安裝方便,施工周期短和防水防火性能優(yōu)越,也可用于露天防火墻和圍墻。AS墻板全干作業(yè)且施工安裝效率高,耐久性強,維護成本低,更符合裝配式變電站的發(fā)展趨勢。

    (9)

    (10)

    需要注意的是,文獻[15]僅僅分析了圓域0°~180°對應的時域振動信號,這就需要確保每段信號的起始點對應的回轉支承的部位是相同的,否則每次分析結果只能對應回轉支承半圈的區(qū)域,但是每次的區(qū)域可能都是不同的,這將使圓域分析的結果在時域上不具有任何物理意義,而從大量的振動數(shù)據(jù)中提取的每段信號都具有相同的相位顯然是很困難的。仔細觀察式(7)~式(10)可以發(fā)現(xiàn),四個指標的計算結果與異常向量V1,V2,…,VM和特征向量Vmean各自的絕對角度是無關的。以圖4為例,即使將圓域分析得到的異常向量統(tǒng)一轉動一定角度,四個指標的計算結果并不會有所改變。利用這一特性,本文所取的每段時域信號都將對應著圓域0°~360°的范圍,這樣盡管每段信號的起始相位并不相同,但是完全不影響計算結果,且圓域分析的物理意義始終指向回轉支承一整圈的健康狀態(tài),從而有效避免了信號相位同步的問題。

    3 試驗驗證

    3.1 加速壽命試驗

    為驗證本文所提方法,研制出如圖6所示回轉支承試驗臺[21]對型號為QNA-730-22的內齒式回轉支承進行了全壽命疲勞試驗。為模擬實際工況,試驗過程中,陪試回轉支承的內圈(定圈)與底座相連,主試回轉支承的內圈(動圈)通過工裝與陪試回轉支承的外圈(動圈)連接,主試回轉支承的外圈(定圈)與上頂蓋相連。這樣,G1, G2兩個油缸便可通過組合加載的方式將載荷加到上頂蓋,再傳遞到主試回轉支承,然后液壓馬達G3驅動陪試回轉支承動圈,帶動主試回轉支承運行,從而模擬出實際工況下回轉支承的運行。

    根據(jù)某公司此型號回轉支承承載能力曲線,主試回轉支承被施以軸向力Fa=96 kN,傾覆力Mn=240 kNm的極限設計載荷,以4 r/min的轉速進行滿載荷疲勞壽命試驗。試驗過程中采用靈敏度為400 mV/g、響應頻率0~1 kHz的Kistler-8315A010加速度傳感器測量回轉支承徑向振動,并以2 kHz的采樣率記錄了其全壽命的振動信號。

    疲勞壽命試驗共進行了11天,試驗進行到第7天上午時,出現(xiàn)了一根螺栓的疲勞斷裂,此時進行了拆機,發(fā)現(xiàn)外圈(定圈)軟帶附近滾道由于疲勞已經出現(xiàn)了區(qū)域小幅滑移,內圈(動圈)上則是出現(xiàn)了輕微點蝕,分別如圖7(a)和圖7(b)所示。由于圖中的初期故障并不影響回轉支承繼續(xù)服役,因此更換了失效的螺栓并加注了潤滑脂后繼續(xù)裝機試驗。至第11天試驗結束時,回轉支承已經完全失效卡死,外圈滾道出現(xiàn)了嚴重的疲勞剝落,部分滾珠甚至產生了疲勞斷裂,如圖7(c)所示。

    (a) 試驗7天后定圈滾道出現(xiàn)區(qū)域滑移

    (b) 試驗7天后動圈滾道出現(xiàn)點蝕

    (c) 嚴重損壞的定圈滾道和滾珠

    試驗中回轉支承的轉速保持在4 r/min,因此圓域分析中每段信號對應回轉支承轉動一圈的時長應為15 s。據(jù)此,取每天6 h、12 h、18 h和24 h的4段信號,每段時長15 s,得到回轉支承全壽命振動信號,如圖8所示。

    圖8同時給出了試驗第1天、第6天和第11天中各1 s的振動信號??梢钥闯觯駝有盘柗翟谠囼炦^程中基本保持增長趨勢,但自始至終很少出現(xiàn)幅值很高的沖擊成分。此外,在試驗初期,回轉支承振動信號以高頻噪聲為主,幾乎無法看到中低頻成分;第6天信號中的高頻噪聲依然存在,但是已經可見部分中頻成分;而第11天的振動信號中則是以低頻為主,中高頻噪聲已經基本不可見。這些現(xiàn)象與前文的分析是一致的,即回轉支承運行初期振動信號基本被白噪聲覆蓋,當回轉支承產生初始故障后,信號中會出現(xiàn)能量較低的沖擊成分,后期隨著故障進一步加劇,沖擊成分能量會不斷升高,信噪比也會持續(xù)提高。為進一步量化分析,對第6天和第11天的信號進行快速傅里葉變換(Fast Fourier Transform,F(xiàn)FT),結果如圖9所示。

    (a) 全壽命信號

    第1天第6天第11天

    (b) 不同時段1 s信號

    圖8 回轉支承全壽命振動信號

    Fig.8 Life-cycle vibration signal of the tested slewing bearing

    圖9(a)展示的是第6天信號的頻譜分析圖,可以看到在<100 Hz處已經出現(xiàn)了部分的低頻成分,但是>200 Hz的中高頻噪聲仍然占據(jù)主導;相對而言,圖9(b)第11天信號的頻譜中,100 Hz以下的低頻成分能量非常高,而300 Hz以上的頻率成分能量很微弱。此外,信號中頻率為119 Hz的成分在試驗初期能量很低,在第6天時能量達到最高,而在最后階段能量又有所降低,這表明119 Hz可能是回轉支承的某種初期故障產生的。

    (a) 第6天信號的FFT

    (b) 第11天信號的FFT

    據(jù)此,根據(jù)圖5選擇w=4對每段15 s的時域振動信號進行圓域分析,由于每段信號的分析過程是完全相同的,此處僅以第6天和第10天各自的第1段信號為例進行討論。為便于理解,將第6天的分析過程在時域內展開,而將第10天的分析過程在圓域內展開。

    首先,將時域振動信號按秒均分成15個區(qū)域,同時將其映射到0°~360°的圓域中,每秒鐘的時域信號對應24°的圓域信號,區(qū)域劃分結果和轉換后的圓域信號分別如圖10和圖12(a)所示;接著,令w=4,對每個角度區(qū)域的信號進行PAA,并在鄰域相關圖中輸出用于橢圓擬合,分別如圖11和圖13所示;然后,判斷每個擬合橢圓的傾角方向,將右傾斜橢圓所在的角度區(qū)域標記為異常向量或異常區(qū)域分別如圖10和圖12(b)所示,其中的紅色向量即表示了異常向量的平均向量,也就是圓域分析的特征向量Vmean。

    圖10 第6天第1段信號的圓域分析(時域內展開)

    (a)0°~24°(b)48°~72°

    圖11 第6天第1段信號的部分擬合橢圓

    圖12 第10天第1段信號的圓域分析(圓域內展開)

    圖13 第10天第1段信號的部分擬合橢圓

    Fig.13 Partial fitted ellipses of the 10th day

    圖11中的兩個擬合橢圓離心率較小,橢圓更接近圓,且只是分別向左和向右有小幅度的傾斜,而圖13中的兩個擬合橢圓的離心率則大很多,傾斜角度也更大。由此可見,在回轉支承壽命初期,各個角度區(qū)域的左傾斜擬合橢圓離心率和傾斜程度均逐步減小,初期故障產生后,擬合橢圓的傾角方向從左傾斜轉換為右傾斜,隨著故障進一步加劇,右傾斜橢圓的離心率不斷增大。此外,由圖10和圖12(b)中可知,試驗中期出現(xiàn)異常的角度區(qū)域較少,而試驗后期異常區(qū)域較多,其特征向量Vmean的模長也越來越短。

    據(jù)此,得到各段信號圓域分析的特征向量Vmean后,計算其均值、方差、歪度和峭度這4個指標在整個試驗周期中的變化趨勢,如圖14所示。為進行對比分析,圖14同時給出了原始振動信號[22]以及小波分析的相關時域指標。其中,小波分析利用db9小波將原始振動信號進行了3層分解,并選取第3層小波系數(shù)中的高頻細分部分(D3)用于時域指標計算,具體方法此處不再贅述。

    (a)均值(b)方差(c)歪度(d)峭度

    圖14 各項指標在整個壽命周期的變化

    Fig.14 Four features of three different methods on the life cycle

    由圖14可知,部分指標在試驗前兩天有一定的波動,對應回轉支承磨合階段,兩天后趨于平穩(wěn),回轉支承正常運行,直到試驗進行至第6天附近,圓域分析的各項指標均出現(xiàn)了較為明顯的異常幅值,說明回轉支承在試驗進行到第6天時產生了初始故障,結合7天后拆機結果分析,定圈的區(qū)域滑移和動圈的點蝕應當是在此時產生的。相比之下,振動信號相關的時域特征在第6天附近幾乎看不到明顯的異常值,而小波分析結果的均值、方差和峭度指標中能觀察到一定程度的異常產生,但是其可識別度遠不如圓域分析的相關指標。由此可見,相對傳統(tǒng)方法,圓域分析能夠更準確地從高噪聲、低信噪比的振動信號中診斷出回轉支承的初始故障,便于工程人員提前制定有效的維護規(guī)劃,縮短維護周期,從而降低維護成本。

    4 結 論

    回轉支承實際工況下故障診斷方法的相關研究仍處于起步階段,本文提出了一種圓域分析方法,將時域振動信號映射到圓域坐標系中,通過鄰域相關圖中擬合橢圓的傾角方向判斷各角度區(qū)域是否異出現(xiàn)常,并將各異常向量的平均向量作為圓域分析的特征向量,最后以特征向量的均值、方差、歪度和峭度指標在整個壽命周期中的變化趨勢進行故障診斷。加速壽命試驗結果表明:大型回轉支承早期的振動信號以高頻噪聲為主,中后期初始故障產生后中低頻沖擊成分不斷增加;鄰域相關離散點的擬合橢圓的傾角方向能夠準確識別出信號中頻率成分的變化,可以作為判斷信號是否出現(xiàn)異常的依據(jù);相比常用的時域特征和小波分析等方法,本方法能夠更為準確地診斷出回轉支承的初期故障,且故障點處的可識別度更高,分析結果的可解釋性更強;將本方法用于回轉支承狀態(tài)監(jiān)測,能夠為制定及時準確的維護規(guī)劃提供切實的參考依據(jù),具有較強的工程應用價值。

    [1] 回轉支承:JB/T 2300—1999[S]. 北京: 機械科學研究院, 1999.

    [2] WIDODO A, YANG B S. Machine health prognostics using survival probability and support vector machine[J]. Expert Systems with Applications, 2011, 38(7): 8430-8437.

    [3] TAN A C C, KIM Y H, KOSSE V. Condition monitoring of low-speed bearings—a review[J]. Australian Journal of Mechanical Engineering, 2008, 6(1): 61-68.

    [4] 張慧芳, 陳捷. 大型回轉支承故障信號處理方法綜述[J]. 機械設計與制造, 2012(3): 216-218.

    ZHANG Huifang, CHEN Jie. Reserch on signal processing method of large slewing bearing[J]. Machinery Design & Manufacture, 2012(3): 216-218.

    [5] ?VOKELJ M, ZUPAN S, PREBIL I. Multivariate and multiscale monitoring of large-size low-speed bearings using ensemble empirical Mode decomposition method combined with principal component analysis[J]. Mechanical Systems and Signal Processing, 2010, 24(4): 1049-1067.

    [6] ?VOKELJ M, ZUPAN S, PREBIL I. Non-linear multivariate and multiscale monitoring and signal denoising strategy using kernel principal component analysis combined with ensemble empirical Mode decomposition method[J]. Mechanical Systems and Signal Processing, 2011, 25(7): 2631-2653.

    [7] 郭剛濤, 秦文科, 李婧銥, 等. 基于混沌理論的大口徑天線回轉支承故障診斷研究[J]. 機械工程師, 2015(2): 133-135.

    GUO Gangtao, QIN Wenke, LI Jingyi, et al. Fault diagnosis of large-diameter antenna slewing bearing based on chaos theory[J]. Mechanical Engineer, 2015(2): 133-135.

    [8] 趙陽, 陳捷, 洪榮晶, 等. EMD 和獨立分量分析在轉盤軸承故障診斷中的應用[J]. 軸承, 2015 (7): 54-59.

    ZHAO Yang, CHEN Jie, HONG Rongjing, et al. Application of EMD and ICA in fault diagnosis for slewing bearings[J]. Bearing, 2015 (7): 54-59.

    [9] 楊杰, 陳捷, 徐新庭, 等. 基于小波-能量模式的回轉支承故障診斷方法研究與應用[J]. 南京工業(yè)大學學報(自然科學版), 2015, 37(4): 134-140.

    YANG Jie, CHEN Jie, XU Xinting, et al. Slewing bearing analysis on fault diagnosis based on wavelet and energy fault mode and its application[J]. Journal of Nanjing Tech University (Natural Science), 2015, 37(4): 134-140.

    [10] BONNARDOT F, RANDALL R B, ANTONI J. Enhanced unsupervised noise cancellation using angular resampling for planetary bearing fault diagnosis[J]. International Journal of Acoustics and Vibrations, 2004, 9(2): 51-60.

    [11] RENAUDIN L, BONNARDOT F, MUSY O, et al. Natural roller bearing fault detection by angular measurement of true instantaneous angular speed[J]. Mechanical Systems and Signal Processing, 2010, 24(7): 1998-2011.

    [12] VILLA L F, RENONES A, PERAN J R, et al. Angular resampling for vibration analysis in wind turbines under non-linear speed fluctuation[J]. Mechanical Systems and Signal Processing, 2011, 25(6): 2157-2168.

    [13] VILLA L F, RENONES A, PERAN J R, et al. Statistical fault diagnosis based on vibration analysis for gear test-bench under non-stationary conditions of speed and load[J]. Mechanical Systems and Signal Processing, 2012, 29(5): 436-446.

    [14] BONNARDOT F, BADAOUI M E, RANDALL R B, et al. Use of the acceleration signal of a gearbox in order to perform angular resampling (with limited speed fluctuation)[J]. Mechanical Systems and Signal Processing, 2005, 19(4): 766-785.

    [15] CAESARENDRA W, KOSASIH B, TIEU A K, et al. Circular domain features based condition monitoring for low speed slewing bearing[J]. Mechanical Systems and Signal Processing, 2014, 45(1): 114-138.

    [16] YI B K, FALOUTSOS C. Fast time sequence indexing for arbitrary Lp norms[EB/OL].[2015-12-14].http://repository.cmu.edu/compsci/553/.

    [17] FITZGIBBON A, PILU M, FISHER R B. Direct least square fitting of ellipses[J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 1999, 21(5): 476-480.

    [18] LUO M F, MATHEW J. Monitoring machinery with varying speed using angle domain analysis[C]// 1993 Proceedings of 17th Annual Meeting. St Louis: Vibration Institute, 1993: 69-78.

    [19] BOSSLEY K M, MCKENDRICK R J, HARRIS C J, et al. Hybrid computed order tracking[J]. Mechanical Systems and Signal Processing, 1999, 13(4): 627-641.

    [20] PEWSEY A. The large-sample joint distribution of key circular statistics[J]. Metrika, 2004, 60(1): 25-32.

    [21] 黃筱調, 王華, 陳捷, 等. 風電轉盤軸承綜合性能實驗臺:CN100587444C[P]. 2010-02-03.

    [22] 王小玲, 陳進, 從飛云. 基于時頻的頻帶熵方法在滾動軸承故障識別中的應用[J]. 振動與沖擊, 2012, 31(18): 29-33.

    WANG Xiaoling, CHEN Jin, CONG Yunfei. Application of spectral band entropy (SBE) method in rolling bearing fault diagnosis based on time-frequency analysis[J]. Journal of Vibration and Shock, 2012, 31(18): 29-33.

    Incipient fault diagnosis of large-size slewing bearings based on circular domain analysis

    FENG Yang, HUANG Xiaodiao, HONG Rongjing, CHEN Jie

    (School of Mechanical and Power Engineering, Nanjing Tech University, Nanjing 210009, China)

    The background noise of a slewing bearing vibration signal in practical load cases is very high, it makes commonly used fault detection approaches not suitable for slewing bearing fault diagnosis. Therefore, a novel signal processing method was proposed based on circular domain analysis. First of all, the time domain signal was transformed into a circular domain and the transformed signal was divided into several zones according to a certain angle, and then the neighborhood correlation discrete points of each zone were fitted as an ellipse. Afterwards, the ellipses skewing to the right were tagged as abnormalities and the corresponding abnormal vectors were obtained based on the whole cycle of a slewing bearing. Finally, the characteristic vector of circular domain analysis, also the mean vector of all the abnormal vectors was acquired, and its mean, variance, skewness and kurtosis were calculated and taken as the fault indicators. An accelerated life test was conducted on a slewing bearing to validate the proposed method. Results showed that the proposed method has a better performance to detect an incipient fault, such as, slipping and pitting in the raceway than the time domain analysis and the wavelet analysis do, it can be an effective tool for slewing bearing fault diagnosis in engineering practice.

    slewing bearing; circular analysis; circular resampling; fault diagnosis; accelerated life test

    國家自然科學基金項目(51375222;51175242)

    2015-12-17 修改稿收到日期:2016-03-16

    封楊 男,博士生,1988年生

    黃筱調 男,教授,博士生導師,1952年生

    TP206.3

    A

    10.13465/j.cnki.jvs.2017.09.017

    猜你喜歡
    時域橢圓故障診斷
    Heisenberg群上由加權次橢圓p-Laplace不等方程導出的Hardy型不等式及應用
    例談橢圓的定義及其應用
    一道橢圓試題的別樣求法
    基于時域信號的三電平逆變器復合故障診斷
    測控技術(2018年11期)2018-12-07 05:49:02
    橢圓的三類切點弦的包絡
    基于極大似然準則與滾動時域估計的自適應UKF算法
    基于時域逆濾波的寬帶脈沖聲生成技術
    因果圖定性分析法及其在故障診斷中的應用
    基于時域波形特征的輸電線雷擊識別
    電測與儀表(2015年2期)2015-04-09 11:28:50
    基于LCD和排列熵的滾動軸承故障診斷
    av国产精品久久久久影院| 99国产精品免费福利视频| 亚洲国产毛片av蜜桃av| 成年人免费黄色播放视频| 免费看十八禁软件| 人妻久久中文字幕网| 一区二区三区国产精品乱码| www.自偷自拍.com| 99国产综合亚洲精品| 国产欧美日韩一区二区精品| 操出白浆在线播放| 亚洲一区二区三区欧美精品| 国产伦人伦偷精品视频| 精品熟女少妇八av免费久了| 国产av精品麻豆| 麻豆av在线久日| 一进一出抽搐动态| 国产亚洲欧美精品永久| 日本三级黄在线观看| 亚洲欧美日韩高清在线视频| 国产成人av激情在线播放| 色尼玛亚洲综合影院| 国产高清videossex| 国产亚洲欧美精品永久| 日韩欧美在线二视频| 欧美成狂野欧美在线观看| 国产精品成人在线| 最近最新中文字幕大全免费视频| 91av网站免费观看| 午夜福利免费观看在线| 久久久国产精品麻豆| 久久久国产精品麻豆| 欧美成人性av电影在线观看| 亚洲五月天丁香| 1024视频免费在线观看| 黄片大片在线免费观看| 别揉我奶头~嗯~啊~动态视频| 亚洲欧美激情综合另类| 精品午夜福利视频在线观看一区| 亚洲精品国产色婷婷电影| 免费在线观看日本一区| 黑人操中国人逼视频| 久久精品亚洲熟妇少妇任你| 成年版毛片免费区| 激情视频va一区二区三区| 搡老岳熟女国产| 99精国产麻豆久久婷婷| 国产成人欧美在线观看| 亚洲色图av天堂| 亚洲片人在线观看| 久久精品成人免费网站| 级片在线观看| 国产成人欧美| 99国产综合亚洲精品| 桃色一区二区三区在线观看| 欧美 亚洲 国产 日韩一| 亚洲美女黄片视频| 亚洲国产欧美一区二区综合| 最近最新中文字幕大全电影3 | 91国产中文字幕| 久久影院123| 天堂俺去俺来也www色官网| 久久精品国产清高在天天线| 中亚洲国语对白在线视频| 国产成人啪精品午夜网站| 一边摸一边抽搐一进一出视频| 亚洲第一av免费看| 18美女黄网站色大片免费观看| 国产激情久久老熟女| 亚洲一区中文字幕在线| 国产精品香港三级国产av潘金莲| 桃红色精品国产亚洲av| 精品高清国产在线一区| 亚洲成人免费av在线播放| 大陆偷拍与自拍| 日韩大码丰满熟妇| 青草久久国产| 国产熟女午夜一区二区三区| 欧美日韩视频精品一区| 在线观看日韩欧美| 亚洲精品国产一区二区精华液| 91大片在线观看| 搡老熟女国产l中国老女人| 激情在线观看视频在线高清| 亚洲精品成人av观看孕妇| 天堂俺去俺来也www色官网| 亚洲成人免费av在线播放| 欧美午夜高清在线| 亚洲狠狠婷婷综合久久图片| 满18在线观看网站| 免费看十八禁软件| 亚洲人成伊人成综合网2020| 免费av中文字幕在线| 两性午夜刺激爽爽歪歪视频在线观看 | 国产欧美日韩综合在线一区二区| 亚洲成a人片在线一区二区| 国产区一区二久久| 亚洲精品国产一区二区精华液| 成人18禁高潮啪啪吃奶动态图| 级片在线观看| 国产1区2区3区精品| 香蕉久久夜色| 黑人巨大精品欧美一区二区mp4| 日韩免费av在线播放| 精品电影一区二区在线| 亚洲中文字幕日韩| 97人妻天天添夜夜摸| 丰满饥渴人妻一区二区三| 麻豆久久精品国产亚洲av | 精品乱码久久久久久99久播| 国产精品日韩av在线免费观看 | 久久久久亚洲av毛片大全| 亚洲成人免费av在线播放| 国产一区在线观看成人免费| 最好的美女福利视频网| 亚洲一区中文字幕在线| 麻豆一二三区av精品| 日韩欧美在线二视频| 久久青草综合色| 国产精品av久久久久免费| 视频区图区小说| 午夜福利影视在线免费观看| 精品国产美女av久久久久小说| 嫁个100分男人电影在线观看| 精品国内亚洲2022精品成人| 国产xxxxx性猛交| 亚洲aⅴ乱码一区二区在线播放 | 在线天堂中文资源库| 日本一区二区免费在线视频| 国产伦一二天堂av在线观看| 久久天堂一区二区三区四区| 欧美色视频一区免费| 亚洲午夜理论影院| 亚洲精品一区av在线观看| 高清黄色对白视频在线免费看| 伦理电影免费视频| 欧美日本亚洲视频在线播放| 一边摸一边抽搐一进一小说| 窝窝影院91人妻| 久久香蕉精品热| 亚洲精品在线观看二区| 欧美午夜高清在线| 久久人妻福利社区极品人妻图片| a在线观看视频网站| 免费av中文字幕在线| 久久香蕉精品热| www.999成人在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲中文日韩欧美视频| 婷婷精品国产亚洲av在线| 99精国产麻豆久久婷婷| 中出人妻视频一区二区| 国产精品99久久99久久久不卡| 宅男免费午夜| tocl精华| 久久 成人 亚洲| 欧美日本亚洲视频在线播放| 午夜福利在线免费观看网站| 亚洲国产欧美一区二区综合| 欧美日韩视频精品一区| 高清毛片免费观看视频网站 | 又大又爽又粗| 最新在线观看一区二区三区| 窝窝影院91人妻| 一二三四在线观看免费中文在| avwww免费| 中文字幕另类日韩欧美亚洲嫩草| 国产精品 国内视频| 老司机午夜十八禁免费视频| 亚洲av五月六月丁香网| 亚洲人成77777在线视频| 熟女少妇亚洲综合色aaa.| 999久久久精品免费观看国产| 成人免费观看视频高清| 婷婷丁香在线五月| 97碰自拍视频| 大香蕉久久成人网| 色老头精品视频在线观看| a级片在线免费高清观看视频| 亚洲人成电影免费在线| 亚洲色图av天堂| 国产一区二区三区综合在线观看| 精品午夜福利视频在线观看一区| 久久久久久免费高清国产稀缺| 国产精品自产拍在线观看55亚洲| 制服人妻中文乱码| 天天影视国产精品| 操出白浆在线播放| 午夜两性在线视频| 中文字幕另类日韩欧美亚洲嫩草| 三级毛片av免费| 欧美日韩精品网址| 久久久精品欧美日韩精品| 女人被狂操c到高潮| 亚洲三区欧美一区| 欧美成人性av电影在线观看| 一级毛片高清免费大全| 少妇被粗大的猛进出69影院| 美女 人体艺术 gogo| 午夜日韩欧美国产| 欧美激情 高清一区二区三区| 日本vs欧美在线观看视频| 在线观看66精品国产| 淫秽高清视频在线观看| 我的亚洲天堂| 日本欧美视频一区| bbb黄色大片| a级毛片黄视频| 国产欧美日韩精品亚洲av| 亚洲免费av在线视频| 九色亚洲精品在线播放| 免费高清在线观看日韩| 欧美大码av| 村上凉子中文字幕在线| 日本一区二区免费在线视频| 精品国产亚洲在线| 老鸭窝网址在线观看| 亚洲色图综合在线观看| 男女下面进入的视频免费午夜 | 国产成年人精品一区二区 | 黄片小视频在线播放| 亚洲欧美激情综合另类| 午夜视频精品福利| 在线观看免费日韩欧美大片| 777久久人妻少妇嫩草av网站| 国产一区二区在线av高清观看| 自拍欧美九色日韩亚洲蝌蚪91| 悠悠久久av| 亚洲黑人精品在线| 美国免费a级毛片| 在线永久观看黄色视频| 性少妇av在线| 久久久久久人人人人人| 久久精品亚洲av国产电影网| 水蜜桃什么品种好| 女性生殖器流出的白浆| 18禁黄网站禁片午夜丰满| 夜夜躁狠狠躁天天躁| 又黄又粗又硬又大视频| 亚洲人成网站在线播放欧美日韩| 老熟妇仑乱视频hdxx| 美女午夜性视频免费| 黑丝袜美女国产一区| 别揉我奶头~嗯~啊~动态视频| 校园春色视频在线观看| 国产熟女xx| 少妇粗大呻吟视频| 亚洲精品一二三| 一边摸一边抽搐一进一小说| 最近最新中文字幕大全电影3 | 亚洲七黄色美女视频| 一级片'在线观看视频| 亚洲av五月六月丁香网| 久久精品国产清高在天天线| 欧美乱妇无乱码| 久久精品人人爽人人爽视色| 在线观看日韩欧美| 757午夜福利合集在线观看| 老汉色av国产亚洲站长工具| 国产精品久久电影中文字幕| 人人妻,人人澡人人爽秒播| 亚洲五月天丁香| 精品国产国语对白av| 亚洲成av片中文字幕在线观看| 女人被狂操c到高潮| 妹子高潮喷水视频| 老鸭窝网址在线观看| 国产成人精品久久二区二区免费| 美女高潮喷水抽搐中文字幕| 亚洲精品国产色婷婷电影| 天天躁狠狠躁夜夜躁狠狠躁| 免费在线观看黄色视频的| 一边摸一边做爽爽视频免费| 色精品久久人妻99蜜桃| 欧美日韩瑟瑟在线播放| 国产精品av久久久久免费| 黄色a级毛片大全视频| 国产又色又爽无遮挡免费看| 黄色片一级片一级黄色片| 侵犯人妻中文字幕一二三四区| 99国产精品免费福利视频| 国产欧美日韩综合在线一区二区| 色综合婷婷激情| 午夜视频精品福利| 他把我摸到了高潮在线观看| a级毛片在线看网站| 国产高清videossex| 欧美日韩国产mv在线观看视频| 老司机靠b影院| av视频免费观看在线观看| 男女下面插进去视频免费观看| 法律面前人人平等表现在哪些方面| 国产精品99久久99久久久不卡| 黄片小视频在线播放| 制服诱惑二区| 99国产精品99久久久久| 国产男靠女视频免费网站| 大型av网站在线播放| 天堂影院成人在线观看| 少妇粗大呻吟视频| 国内毛片毛片毛片毛片毛片| 美女扒开内裤让男人捅视频| 国产精品亚洲一级av第二区| 日本黄色日本黄色录像| 国产麻豆69| 欧美日韩av久久| 中文字幕人妻丝袜一区二区| bbb黄色大片| 成年版毛片免费区| 亚洲精品美女久久av网站| e午夜精品久久久久久久| 久久国产精品影院| 18美女黄网站色大片免费观看| 免费日韩欧美在线观看| 国产成人精品久久二区二区91| 女人爽到高潮嗷嗷叫在线视频| 99在线人妻在线中文字幕| 国产日韩一区二区三区精品不卡| a级片在线免费高清观看视频| 多毛熟女@视频| 亚洲国产毛片av蜜桃av| 精品欧美一区二区三区在线| 国产xxxxx性猛交| 女人被狂操c到高潮| av天堂在线播放| 99久久国产精品久久久| 18禁观看日本| 久久久精品国产亚洲av高清涩受| 神马国产精品三级电影在线观看 | 三上悠亚av全集在线观看| 久久久久九九精品影院| 亚洲 欧美 日韩 在线 免费| 视频在线观看一区二区三区| 国产精品一区二区精品视频观看| 亚洲黑人精品在线| 一二三四在线观看免费中文在| 在线看a的网站| 国产免费男女视频| 在线看a的网站| 欧美性长视频在线观看| 亚洲av第一区精品v没综合| 午夜福利影视在线免费观看| 51午夜福利影视在线观看| 成人永久免费在线观看视频| 亚洲精品中文字幕一二三四区| 国产伦一二天堂av在线观看| 18禁美女被吸乳视频| 无限看片的www在线观看| 叶爱在线成人免费视频播放| 免费在线观看视频国产中文字幕亚洲| 少妇裸体淫交视频免费看高清 | 色综合站精品国产| 9191精品国产免费久久| 免费观看人在逋| 国产伦一二天堂av在线观看| 亚洲精品中文字幕一二三四区| 国产精品久久久久成人av| 日韩欧美免费精品| 一本综合久久免费| 在线看a的网站| 757午夜福利合集在线观看| 51午夜福利影视在线观看| 欧美 亚洲 国产 日韩一| 性色av乱码一区二区三区2| 精品久久久精品久久久| 日韩欧美免费精品| 美女国产高潮福利片在线看| 亚洲黑人精品在线| 国产视频一区二区在线看| www.www免费av| 久久香蕉精品热| 嫁个100分男人电影在线观看| 一边摸一边做爽爽视频免费| tocl精华| 欧美+亚洲+日韩+国产| 国产精品秋霞免费鲁丝片| 波多野结衣高清无吗| 极品人妻少妇av视频| 亚洲少妇的诱惑av| 欧美成人免费av一区二区三区| 精品国产一区二区三区四区第35| 多毛熟女@视频| 在线永久观看黄色视频| 黄片播放在线免费| 在线天堂中文资源库| 中文字幕av电影在线播放| 视频区欧美日本亚洲| 麻豆一二三区av精品| 一个人免费在线观看的高清视频| 又黄又粗又硬又大视频| 啪啪无遮挡十八禁网站| 黄网站色视频无遮挡免费观看| 18禁裸乳无遮挡免费网站照片 | 我的亚洲天堂| 亚洲av日韩精品久久久久久密| 丝袜美足系列| 国内毛片毛片毛片毛片毛片| 超碰成人久久| 久久精品91无色码中文字幕| 一区二区三区激情视频| 在线观看免费午夜福利视频| 乱人伦中国视频| 国产97色在线日韩免费| 啦啦啦免费观看视频1| 久久久国产一区二区| av在线播放免费不卡| 色尼玛亚洲综合影院| a级片在线免费高清观看视频| 黄片播放在线免费| 久久久水蜜桃国产精品网| 久久国产亚洲av麻豆专区| 久久久精品欧美日韩精品| 一本大道久久a久久精品| 国产成人精品无人区| 一级a爱视频在线免费观看| 高清毛片免费观看视频网站 | 亚洲国产欧美一区二区综合| 亚洲 欧美 日韩 在线 免费| 一进一出抽搐动态| 亚洲国产毛片av蜜桃av| 国产av一区在线观看免费| 亚洲男人天堂网一区| 啦啦啦在线免费观看视频4| 成人免费观看视频高清| 老司机深夜福利视频在线观看| 国产成人精品无人区| 麻豆成人av在线观看| 成年人黄色毛片网站| 亚洲专区中文字幕在线| 国产成人欧美在线观看| 在线永久观看黄色视频| 欧美最黄视频在线播放免费 | 国产高清视频在线播放一区| 国产区一区二久久| 国产高清激情床上av| 一级作爱视频免费观看| 神马国产精品三级电影在线观看 | 美女高潮喷水抽搐中文字幕| 精品高清国产在线一区| 欧美成人性av电影在线观看| av电影中文网址| 90打野战视频偷拍视频| 人妻丰满熟妇av一区二区三区| 丰满的人妻完整版| av片东京热男人的天堂| 日韩精品青青久久久久久| 在线观看午夜福利视频| 国产精品一区二区精品视频观看| 久久99一区二区三区| 久久精品亚洲熟妇少妇任你| 国产亚洲精品久久久久5区| 亚洲午夜理论影院| 黄色视频不卡| 久久国产精品人妻蜜桃| 精品一品国产午夜福利视频| 成人三级做爰电影| 久久人人精品亚洲av| aaaaa片日本免费| 三上悠亚av全集在线观看| 又大又爽又粗| 欧美日韩av久久| 久久99一区二区三区| av电影中文网址| videosex国产| 午夜免费鲁丝| 日韩欧美一区视频在线观看| 色婷婷av一区二区三区视频| 999久久久国产精品视频| 制服诱惑二区| 天堂俺去俺来也www色官网| 交换朋友夫妻互换小说| 老司机深夜福利视频在线观看| 美女高潮到喷水免费观看| 精品久久久久久久久久免费视频 | 国产精品二区激情视频| 国产区一区二久久| 老熟妇仑乱视频hdxx| 香蕉国产在线看| 亚洲全国av大片| 国产精品亚洲av一区麻豆| 亚洲片人在线观看| 亚洲色图综合在线观看| 亚洲七黄色美女视频| 一区二区三区激情视频| 色婷婷av一区二区三区视频| av视频免费观看在线观看| 可以在线观看毛片的网站| 91成年电影在线观看| 亚洲一区中文字幕在线| 国产熟女午夜一区二区三区| 丁香欧美五月| 国产日韩一区二区三区精品不卡| 免费高清在线观看日韩| 在线观看66精品国产| 99热国产这里只有精品6| 天天影视国产精品| 久久久久久久久免费视频了| 欧美日韩亚洲综合一区二区三区_| 黄色成人免费大全| 国产成人av教育| 日本欧美视频一区| 亚洲欧美日韩高清在线视频| 黄频高清免费视频| 欧美日韩黄片免| 首页视频小说图片口味搜索| 亚洲第一青青草原| 女人爽到高潮嗷嗷叫在线视频| 18美女黄网站色大片免费观看| 日本免费一区二区三区高清不卡 | 手机成人av网站| 一区在线观看完整版| 男女午夜视频在线观看| 亚洲欧洲精品一区二区精品久久久| 女性生殖器流出的白浆| 97人妻天天添夜夜摸| 久99久视频精品免费| 国产不卡一卡二| 精品国内亚洲2022精品成人| 高清黄色对白视频在线免费看| 亚洲精品中文字幕在线视频| 大码成人一级视频| 99精品欧美一区二区三区四区| 香蕉久久夜色| 91国产中文字幕| 无遮挡黄片免费观看| 欧美久久黑人一区二区| 国产亚洲精品久久久久5区| 丰满饥渴人妻一区二区三| 十分钟在线观看高清视频www| 久久中文看片网| 91九色精品人成在线观看| a在线观看视频网站| 女人高潮潮喷娇喘18禁视频| 久久国产精品人妻蜜桃| 日本免费a在线| 国产一区二区三区视频了| 成年女人毛片免费观看观看9| 国产精品二区激情视频| 国产精品乱码一区二三区的特点 | 丁香六月欧美| 夜夜爽天天搞| 乱人伦中国视频| 国产高清国产精品国产三级| 国产精品av久久久久免费| 天堂俺去俺来也www色官网| 首页视频小说图片口味搜索| 中文字幕人妻丝袜一区二区| 欧美精品啪啪一区二区三区| av福利片在线| 午夜福利一区二区在线看| 亚洲国产欧美网| 久久天躁狠狠躁夜夜2o2o| 99久久人妻综合| 国产aⅴ精品一区二区三区波| 99国产综合亚洲精品| 亚洲 欧美 日韩 在线 免费| 亚洲第一欧美日韩一区二区三区| 亚洲成av片中文字幕在线观看| 99国产精品一区二区三区| 69精品国产乱码久久久| 久久精品aⅴ一区二区三区四区| 亚洲中文日韩欧美视频| 十八禁人妻一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 久久青草综合色| 母亲3免费完整高清在线观看| 日本免费一区二区三区高清不卡 | 亚洲人成网站在线播放欧美日韩| 在线av久久热| 99在线人妻在线中文字幕| 9热在线视频观看99| 精品无人区乱码1区二区| 欧美老熟妇乱子伦牲交| 久久精品亚洲精品国产色婷小说| 一区在线观看完整版| 日韩欧美一区视频在线观看| 亚洲中文日韩欧美视频| 久久精品成人免费网站| 亚洲欧美日韩无卡精品| 9色porny在线观看| 精品免费久久久久久久清纯| 国产99白浆流出| 亚洲成国产人片在线观看| 黄片小视频在线播放| www日本在线高清视频| 中文字幕人妻熟女乱码| 国产精品国产av在线观看| 亚洲专区字幕在线| 欧美日韩视频精品一区| 国产高清激情床上av| 丁香六月欧美| 高清毛片免费观看视频网站 | 在线观看免费午夜福利视频| 国产精品久久久av美女十八| 久久国产亚洲av麻豆专区| 99riav亚洲国产免费| 日本欧美视频一区| 在线国产一区二区在线| 国产色视频综合| 亚洲片人在线观看| 女人爽到高潮嗷嗷叫在线视频| 国产欧美日韩一区二区精品| 80岁老熟妇乱子伦牲交| 高清欧美精品videossex| 久9热在线精品视频| 丰满的人妻完整版| 国产精品久久久久成人av| bbb黄色大片| 欧美黄色片欧美黄色片| 老熟妇乱子伦视频在线观看| 两人在一起打扑克的视频| 在线永久观看黄色视频| 美女福利国产在线| 黄频高清免费视频| 欧美成人午夜精品|