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

    滾動(dòng)軸承振動(dòng)信號(hào)處理方法綜述

    2016-02-07 12:20:28胡智勇胡杰鑫謝里陽(yáng)佟安時(shí)張嘯塵
    關(guān)鍵詞:故障診斷振動(dòng)故障

    胡智勇,胡杰鑫,謝里陽(yáng),佟安時(shí),張嘯塵

    (1.東北大學(xué) 航空動(dòng)力裝備振動(dòng)及控制教育部重點(diǎn)實(shí)驗(yàn)室,遼寧 沈陽(yáng) 110819;2.東北大學(xué) 現(xiàn)代設(shè)計(jì)與分析研究所,遼寧 沈陽(yáng) 110819)

    滾動(dòng)軸承振動(dòng)信號(hào)處理方法綜述

    胡智勇,胡杰鑫,謝里陽(yáng),佟安時(shí),張嘯塵

    (1.東北大學(xué) 航空動(dòng)力裝備振動(dòng)及控制教育部重點(diǎn)實(shí)驗(yàn)室,遼寧 沈陽(yáng) 110819;2.東北大學(xué) 現(xiàn)代設(shè)計(jì)與分析研究所,遼寧 沈陽(yáng) 110819)

    針對(duì)軸承振動(dòng)信號(hào)處理的各種方法進(jìn)行了全面綜述,研究了軸承振動(dòng)信號(hào)處理的時(shí)域分析法和頻域分析法,指出時(shí)域參數(shù)的具體用法和適用范圍,以及各種頻域分析方法的優(yōu)缺點(diǎn);闡述了時(shí)頻分析方法,該方法克服了時(shí)域與頻域分析不能反映信號(hào)局部特征的缺陷,且對(duì)非平穩(wěn)的故障信號(hào)有較強(qiáng)的處理能力;介紹了一些特殊的分析方法,對(duì)一些具體的軸承振動(dòng)信號(hào)問(wèn)題有較好的處理效果.分析現(xiàn)有的軸承振動(dòng)信號(hào)處理方法,表明各種分析方法取長(zhǎng)補(bǔ)短綜合運(yùn)用以及新技術(shù)、新理論的應(yīng)用,將會(huì)是軸承振動(dòng)信號(hào)處理方法未來(lái)的發(fā)展方向.

    滾動(dòng)軸承; 信號(hào)處理; 振動(dòng); 時(shí)域分析; 頻域分析

    滾動(dòng)軸承是一種將轉(zhuǎn)軸與軸座之間的滑動(dòng)摩擦轉(zhuǎn)變?yōu)闈L動(dòng)摩擦的精密機(jī)械元件.由于加工或安裝誤差,以及服役過(guò)程中的疲勞、磨損、腐蝕等,都將會(huì)使軸承產(chǎn)生振動(dòng),直接影響整個(gè)轉(zhuǎn)動(dòng)系統(tǒng)的轉(zhuǎn)動(dòng)精度與平穩(wěn)性,甚至?xí)l(fā)嚴(yán)重的機(jī)械事故.因此,有必要對(duì)軸承的運(yùn)行狀態(tài)進(jìn)行監(jiān)測(cè),以便及時(shí)發(fā)現(xiàn)問(wèn)題,避免產(chǎn)生嚴(yán)重的后果,從而可以降低故障發(fā)生率,減少維修費(fèi)用,提高經(jīng)濟(jì)效益.由此可見(jiàn),研究軸承故障監(jiān)測(cè)與診斷的振動(dòng)信號(hào)處理方法十分必要.

    軸承發(fā)生局部缺陷時(shí),會(huì)以一定頻率產(chǎn)生一系列帶寬沖擊和沖擊衰減響應(yīng),表現(xiàn)為復(fù)雜的調(diào)幅振動(dòng);由于缺陷位置不同,振動(dòng)特性也不完全相同[1,2].對(duì)軸承振動(dòng)信號(hào)的處理來(lái)研究軸承運(yùn)行狀態(tài),是一種非常有效的檢測(cè)、診斷方法.

    國(guó)際上,美國(guó)學(xué)者是最早進(jìn)行滾動(dòng)軸承故障診斷研究的,提出了用加速度傳感器來(lái)采集振動(dòng)信號(hào)進(jìn)行檢測(cè).英國(guó)人Dyer采用峭度系數(shù)法、Taylor提出尺寸故障計(jì)算方法、Randall采用倒頻譜分析、Mathew采用共振解調(diào)法和AR模型等對(duì)軸承故障監(jiān)測(cè)進(jìn)行的研究也都取得了成果.小波變換在信號(hào)處理方面的優(yōu)越性,使其在軸承故障診斷中得到應(yīng)用廣泛;最小熵反卷積與譜峭度結(jié)合,可有效地對(duì)軸承微弱故障進(jìn)行檢測(cè)[3-5].一些新的方法,如Hilbert-Huang變換(經(jīng)驗(yàn)?zāi)B(tài)分解方法)、分型理論、局部值分解法等,也都應(yīng)用在滾動(dòng)軸承故障中.

    國(guó)內(nèi)軸承診斷技術(shù)起步較晚,起初主要是學(xué)習(xí)國(guó)外的先進(jìn)理論與方法,近年來(lái)滾動(dòng)軸承信號(hào)處理方法也有了較大的發(fā)展.李昊[6]在研究小波變換消噪性能的基礎(chǔ)上,提出了一種小波循環(huán)自相關(guān)函數(shù)分析算法,應(yīng)用在了滾動(dòng)軸承振動(dòng)信號(hào)的處理.陳曉平等[7]闡述了利用小波閾值消噪技術(shù)和小波包特征提取方法在滾動(dòng)軸承振動(dòng)信號(hào)處理中的應(yīng)用,大幅度提高了信噪比,同時(shí)也提高了特征信息的靈敏度和分辨率.李毓輝等[8]結(jié)合滾動(dòng)軸承故障信號(hào)譜結(jié)構(gòu)的先驗(yàn)知識(shí),提出了一種半盲均衡方法,實(shí)現(xiàn)了盲均衡器參數(shù)的優(yōu)化,實(shí)驗(yàn)數(shù)據(jù)驗(yàn)證了該方法的有效性.

    1 滾動(dòng)軸承振動(dòng)機(jī)理

    滾動(dòng)軸承振動(dòng)是由外部振源或軸承損傷引起的,潤(rùn)滑劑在軸承運(yùn)轉(zhuǎn)時(shí)產(chǎn)生的流體動(dòng)力也可以是振源.在研究滾動(dòng)軸承振動(dòng)信號(hào)過(guò)程中,主要研究軸承因損壞而產(chǎn)生的振動(dòng),損傷原因主要有磨損、腐蝕、疲勞等.

    軸承的磨損主要發(fā)生在滾道與滾動(dòng)體之間,原因可能有異物進(jìn)入、載荷沖擊、潤(rùn)滑不良等.磨損會(huì)使軸承的游隙、表面粗糙度增大,旋轉(zhuǎn)精度降低,同時(shí)也會(huì)加劇軸承磨損,在振動(dòng)方面主要表現(xiàn)為振動(dòng)幅度增大、產(chǎn)生噪聲.工作表面損傷后,損傷點(diǎn)受力時(shí)會(huì)產(chǎn)生突變的沖擊力,其較寬的信號(hào)頻帶會(huì)覆蓋軸承系統(tǒng)高頻處的固有頻率,引起系統(tǒng)簡(jiǎn)諧振動(dòng)[9].

    損傷引起的振動(dòng)有兩種類型:一種是由軸承損傷點(diǎn)處受到反復(fù)撞擊而產(chǎn)生的,其頻率與軸承的性質(zhì)有關(guān);另一種是因損傷沖擊而誘發(fā)的高頻固有振動(dòng)[10].

    2 信號(hào)處理方法

    2.1 時(shí)域分析法

    時(shí)域分析是根據(jù)信號(hào)的時(shí)間歷程記錄波形,分析信號(hào)的組成和特征量,從而在時(shí)間坐標(biāo)軸上顯示出振動(dòng)信息的方法,包括時(shí)域波形分析和時(shí)域參數(shù)值分析.滾動(dòng)軸承的振動(dòng)監(jiān)測(cè),一般采用加速度計(jì)來(lái)采集振動(dòng)信號(hào)的時(shí)域波形.

    時(shí)域波形作為原始波形,保留了最全面、最完整的信息,為防止有用信息的丟失,在進(jìn)行分析時(shí)不加處理,直接進(jìn)行觀察.

    在時(shí)域分析方法中,信號(hào)的基本統(tǒng)計(jì)特征和信號(hào)的概率分布特征在軸承的故障監(jiān)測(cè)與診斷中起到很大的作用.其中,應(yīng)用較為廣泛的時(shí)域指標(biāo)包括峰值、有效值、峰值指標(biāo)、脈沖指標(biāo)、裕度指標(biāo)、峭度系數(shù)等.

    峰值:反映振幅的最大值,可能是由于點(diǎn)蝕而引起振動(dòng)沖擊,是用于監(jiān)測(cè)軸承點(diǎn)蝕的時(shí)域指標(biāo);同樣,比值形式的脈沖指標(biāo)、裕度指標(biāo)和峰值指標(biāo)分子均為峰值,同樣可以根據(jù)其值的變化來(lái)監(jiān)測(cè)點(diǎn)蝕.

    有效值:反映振幅隨時(shí)間的緩慢變化,而磨損也是一種由緩及重的變化過(guò)程,因此可以用來(lái)檢測(cè)軸承的磨損程度.

    峰值指標(biāo):峰值與有效值的比,當(dāng)出現(xiàn)點(diǎn)蝕產(chǎn)生損傷時(shí),峰值加大,而有效值變化緩慢,其指標(biāo)值增加,后又因有效值的增大,比值減小,由此可以判斷軸承是否發(fā)生故障.

    峰值、脈沖指標(biāo)、裕度指標(biāo)等是對(duì)軸承早期故障敏感的參數(shù),在早期軸承監(jiān)測(cè)中,其值變化較大;有效值因變化不太明顯,不適合早期的軸承故障診斷,但其較好的穩(wěn)定性,與前三個(gè)參數(shù)相結(jié)合,也可以很好地進(jìn)行軸承的故障監(jiān)測(cè)與診斷.

    峭度系數(shù)K,表示樣本密度函數(shù)的頂峰凸平度和觀測(cè)值的分散程度,可以分為正峭度(K>3)、零峭度(K=3)和負(fù)峭度(K<3),此參數(shù)不適合軸承晚期的故障診斷.

    時(shí)域分析方法中,峰值與峭度系數(shù)等均為無(wú)量綱特征參數(shù),克服了軸承及其工作特性的影響.整體來(lái)說(shuō),時(shí)域分析法的處理過(guò)程不會(huì)使信號(hào)發(fā)生畸變或損失,使用最原始資料進(jìn)行軸承的故障診斷;時(shí)域分析法不受信號(hào)的絕對(duì)水平的影響,不論測(cè)點(diǎn)是否有變動(dòng),對(duì)信號(hào)處理的結(jié)果影響不大.時(shí)域法適合于軸承故障監(jiān)測(cè)的早期診斷,并且時(shí)域分析法的精密程度不夠好.

    2.2 頻域分析法

    滾動(dòng)軸承的頻域分析法較時(shí)域分析法精密度要高,其原理是將時(shí)域信號(hào)變換到頻域中,根據(jù)頻譜分析來(lái)對(duì)軸承進(jìn)行監(jiān)測(cè).在軸承振動(dòng)信號(hào)提取過(guò)程中會(huì)混入大量的背景噪聲,這會(huì)對(duì)頻譜分析產(chǎn)生不利的影響,因此需要先進(jìn)行濾波處理,提高信噪比.頻域分析法的理論基礎(chǔ)是傅里葉變換法,常用的具體方法有倒頻譜分析、共振解調(diào)技術(shù)和譜峭度法.

    倒頻譜分析法也稱為二次頻譜分析,是在倒頻譜上識(shí)別軸承振動(dòng)信號(hào)頻率組成分量的方法.具體方法是將功率譜的對(duì)數(shù)值進(jìn)行傅里葉變換,把復(fù)雜的卷積關(guān)系轉(zhuǎn)變?yōu)楹?jiǎn)單的線性疊加.倒頻譜分析抗傳遞途徑的干擾能力強(qiáng),即所測(cè)得的信號(hào)受信號(hào)的傳播路徑影響小,可以較為真實(shí)地反映振源的特性.但倒頻譜分析對(duì)于軸承出現(xiàn)多點(diǎn)損傷時(shí)的監(jiān)測(cè)有一定的盲目性,可能會(huì)降低診斷精度.

    共振解調(diào)技術(shù)的具體原理為:當(dāng)軸承的工作表面出現(xiàn)損傷時(shí),損傷點(diǎn)的接觸動(dòng)作會(huì)產(chǎn)生很大的沖擊脈沖力,脈沖的帶寬覆蓋于軸承系統(tǒng)的整個(gè)固有頻率,會(huì)引起系統(tǒng)的固有振動(dòng),也會(huì)使軸承在各個(gè)固有頻率上發(fā)生諧振,并伴有幅值的調(diào)制現(xiàn)象;根據(jù)實(shí)際需要,選擇某一頻率作為研究對(duì)象,利用高頻帶通濾波器將故障信息解調(diào)出來(lái),再經(jīng)過(guò)低通濾波器處理,得到相對(duì)低頻沖擊并放大擴(kuò)展的共振諧波,處理后的信號(hào)包含故障特征信息的低頻包絡(luò)信號(hào),通過(guò)對(duì)這一包絡(luò)信號(hào)進(jìn)行頻譜分析就可以對(duì)軸承進(jìn)行故障監(jiān)測(cè)與分析.共振解調(diào)技術(shù)在很大程度上提高了信噪比,對(duì)軸承微弱的沖擊信號(hào)有很強(qiáng)的敏感性,是軸承故障診斷較為有效的方法;但該方法也存在局限性,即帶通濾波器參數(shù)的選擇.

    譜峭度法可以為共振解調(diào)技術(shù)提供準(zhǔn)確的帶通濾波器參數(shù)[11].Dwyer最早提出譜峭度這一概念,目的是為了解決功率譜無(wú)法檢測(cè)和提取信號(hào)中瞬態(tài)現(xiàn)象等問(wèn)題.起初該方法并未得到普遍應(yīng)用,主要因?yàn)槠涠x不完善、估計(jì)程序不完整.J.Antoni對(duì)譜峭度法進(jìn)行了深入研究,給出了譜峭度的正式定義,才使其在故障診斷領(lǐng)域的應(yīng)用變得成熟.

    其他的頻域分析法還有傅里葉變換法、功率譜分析法等.在信號(hào)處理中,傅里葉變換的典型用途是將信號(hào)分解成幅值譜,用來(lái)顯示與頻率對(duì)應(yīng)的幅值大小.傅里葉變換法不會(huì)造成信息的丟失,通過(guò)觀察各頻率的譜線來(lái)進(jìn)行故障的監(jiān)測(cè)與診斷.功率譜分析包含兩種方法,自功率譜密度分析與互功率譜密度分析,都是研究信號(hào)在頻域中的各種特征,從頻域內(nèi)提取淹沒(méi)在噪聲中有用的信號(hào).

    2.3 時(shí)頻分析法

    軸承振動(dòng)信號(hào)時(shí)域與頻域分析僅僅能從整體對(duì)信號(hào)進(jìn)行處理,無(wú)法反映信號(hào)局部的特征.時(shí)頻分析方法則可以反映信號(hào)頻率與時(shí)間之間的相互關(guān)系,以及信號(hào)的局部特征,并且對(duì)非平穩(wěn)的故障信號(hào)有較強(qiáng)的處理能力.時(shí)頻分析法主要包括小波變換、短時(shí)傅里葉變換、威格納時(shí)頻分布和EMD方法等.

    在信號(hào)提取過(guò)程中不可避免會(huì)帶來(lái)噪聲的干擾,如果不加必要的處理,會(huì)使特征的提取出現(xiàn)問(wèn)題,正確性降低,必須采取相應(yīng)的措施進(jìn)行消噪的處理.采集到的非平穩(wěn)振動(dòng)信號(hào),利用傳統(tǒng)消噪方法處理獲得的結(jié)果不理想,因此可以利用小波分析來(lái)進(jìn)行消噪.小波分析是一種多分辨率的方法,在時(shí)頻平面不同位置具有不同的分辨率,其克服了傳統(tǒng)信號(hào)處理方法只有單一分辨率的局限.

    小波分解法還可以進(jìn)一步改進(jìn)為小波包分解.小波包分解是將小波分解得到的高頻成分進(jìn)一步分解,克服了小波分解只對(duì)低頻成分進(jìn)行分解的缺點(diǎn),是一種更加細(xì)化的分解.當(dāng)滾動(dòng)軸承出現(xiàn)故障時(shí),某些頻帶內(nèi)的振動(dòng)信號(hào)能量將減少,而另一些頻帶內(nèi)振動(dòng)信號(hào)能量將增加,各頻率成分都含有有用的特征信息,利用小波包分解提取后,直接根據(jù)故障頻率的映射關(guān)系便可判別出軸承的故障原因.因此,小波包分解是一種十分有效的軸承振動(dòng)信號(hào)處理方法.

    短時(shí)傅里葉變換不能同時(shí)兼顧頻率與時(shí)間分辨率:當(dāng)頻率增高時(shí),時(shí)間分辨率就會(huì)降低;反之,時(shí)間分辨率就會(huì)增高.而當(dāng)信號(hào)變化劇烈時(shí),要求窗函數(shù)具有較高的時(shí)間分辨率;波形變化平緩時(shí),則要求窗函數(shù)具有較高的頻率分辨率.因此,短時(shí)傅里葉變換僅適用于分段平穩(wěn)信號(hào)或近似平穩(wěn)信號(hào),并不適用于非平穩(wěn)信號(hào).

    威格納(Wigner-Ville)時(shí)頻分布有利于處理非平穩(wěn)和時(shí)變信號(hào).這種方法由Wigner提出,由Ville將其引入信號(hào)分析領(lǐng)域.威格納時(shí)頻分布比短時(shí)傅里葉變換有更好的分辨率,但其卻存在交叉干擾項(xiàng),迄今這一問(wèn)題也未能很好地解決.

    經(jīng)驗(yàn)?zāi)B(tài)分解(EMD)方法是由美國(guó)NASA的黃鍔博士提出的,與傅里葉分解與小波分解方法相比,該方法不需要設(shè)定任何基函數(shù),是一種直觀的、自適應(yīng)的、后驗(yàn)的方法.EMD方法可以處理任何類型的信號(hào),尤其適用于非平穩(wěn)信號(hào).

    局部均值分解(LMD)[12]方法是由Smith提出的,是一種新的非線性和非平穩(wěn)信號(hào)分析方法.其核心思想是將復(fù)雜的多分量信號(hào)分解為若干個(gè)PF(Product function,簡(jiǎn)稱PF)分量之和,每個(gè)分量都是包絡(luò)信號(hào)和純調(diào)頻信號(hào)的乘積,較好地保持了原始信號(hào)的幅值和頻率的變化特性,將所有的分量進(jìn)行一定的組合便可以得到原始信號(hào)的時(shí)頻分布.這種方法與經(jīng)驗(yàn)?zāi)B(tài)分解法有相似之處,并且具有較高的信噪比,適用于非平穩(wěn)、非線性的信號(hào)處理.

    2.4 其他信號(hào)處理方法

    G.L.McDonald等人提出了一種新的信號(hào)分析處理方法——最大相關(guān)峭度反卷積(MCKD)方法,此方法通過(guò)選取一個(gè)有限沖擊響應(yīng)濾波器,使已知周期信號(hào)濾波后相關(guān)峭度最大,迭代求取最大相關(guān)峭度反卷積,此方法在診斷周期性的沖擊和抑制噪聲方面有很好的效果.在MCKD方法的基礎(chǔ)上,徐繼剛提出了高階位移相關(guān)峭度反卷積方法,在對(duì)采樣信號(hào)進(jìn)行降噪的同時(shí)強(qiáng)化周期性的沖擊成分[13].該方法可以提高故障監(jiān)測(cè)能力,但高階也增加了計(jì)算的復(fù)雜性.在信號(hào)處理時(shí),正確選取T(周期參數(shù))與M(移位數(shù))是十分關(guān)鍵的,如果M取值較大時(shí),迭代法會(huì)降低數(shù)值的精度.

    李毓輝等綜合運(yùn)用盲均衡技術(shù)和計(jì)算階次跟蹤技術(shù)提出一種滾動(dòng)軸承振動(dòng)信號(hào)半盲處理方法,該方法結(jié)合了滾動(dòng)軸承故障信號(hào)譜結(jié)構(gòu)的先驗(yàn)知識(shí),利用計(jì)算階次跟蹤分析得到盲均衡后振動(dòng)信號(hào)的階次譜圖,然后根據(jù)故障信號(hào)的譜結(jié)構(gòu)實(shí)現(xiàn)盲均衡濾波器參數(shù)N(每個(gè)通道的樣本數(shù))和τ(時(shí)間延遲數(shù))的優(yōu)化.該方法基于任意旋轉(zhuǎn)機(jī)械的信號(hào)都是具有循環(huán)平穩(wěn)的特性,通過(guò)數(shù)據(jù)關(guān)聯(lián)和數(shù)據(jù)級(jí)信息融合,得到的信號(hào)要優(yōu)于原始傳感器獲得的信號(hào),這在工程中有利于對(duì)軸承故障的診斷.

    唐國(guó)明等提出了基于虛擬儀器,運(yùn)用LabVIEW軟件,應(yīng)用自回歸模型功率譜分析AR(p)方法對(duì)滾動(dòng)軸承信號(hào)進(jìn)行處理[14].虛擬儀器由軟、硬件平臺(tái)組成,用LabVIEW語(yǔ)言編寫(xiě)程序,而且AR(p)方法易于在計(jì)算機(jī)上實(shí)現(xiàn).在軸承信號(hào)的處理方法中,頻域分析應(yīng)用較為廣泛,頻域參數(shù)的獲得大多數(shù)是通過(guò)傅里葉變換,其缺點(diǎn)是計(jì)算量非常大.而AR(p)功率譜處理信息具有過(guò)程簡(jiǎn)單、計(jì)算量小、模型參數(shù)估計(jì)較為準(zhǔn)確,得到的譜線光滑、譜峰突出等優(yōu)點(diǎn),準(zhǔn)確的AR(p)能夠較好的反映出動(dòng)態(tài)系統(tǒng)的客觀規(guī)律,為研究軸承振動(dòng)的本質(zhì)提供可靠的數(shù)據(jù),但AR(p)對(duì)噪聲干擾也較為敏感.虛擬儀器的使用,能實(shí)現(xiàn)傳統(tǒng)儀器的功能,隨著科技的不斷發(fā)展,虛擬儀器可以更廣泛的運(yùn)用到故障監(jiān)測(cè)與診斷當(dāng)中.

    數(shù)學(xué)形態(tài)學(xué)信號(hào)處理方法[15]:形態(tài)學(xué)最初是應(yīng)用在圖像處理領(lǐng)域,由法國(guó)地質(zhì)學(xué)家Matheron和Serra創(chuàng)立,后來(lái)在信號(hào)處理方面也得到了應(yīng)用.相對(duì)于小波變換、短時(shí)傅里葉變換、威格納時(shí)頻分布等非線性非平穩(wěn)信號(hào)的處理方法,數(shù)學(xué)形態(tài)法的濾波效果更好.利用其處理滾動(dòng)軸承振動(dòng)信號(hào)時(shí),根據(jù)軸承振動(dòng)特點(diǎn)引入形態(tài)學(xué)方法,構(gòu)造出有降噪功能的形態(tài)學(xué)濾波器對(duì)信號(hào)進(jìn)行處理,還原振動(dòng)的真實(shí)信號(hào),從而對(duì)軸承的故障進(jìn)行監(jiān)測(cè)與診斷.數(shù)學(xué)形態(tài)學(xué)在振動(dòng)故障診斷中的應(yīng)用才剛剛起步,還有待于進(jìn)一步地完善.由于數(shù)學(xué)形態(tài)學(xué)的濾波特性,小波變換的多分辨率特性,二者的結(jié)合使用是未來(lái)重點(diǎn)深入研究的內(nèi)容.

    模糊識(shí)別方法:模糊診斷的前期工作是建立隸屬矩陣R,即故障與征兆之間的模糊關(guān)系矩陣,這種矩陣可以通過(guò)經(jīng)驗(yàn)數(shù)據(jù)、統(tǒng)計(jì)數(shù)據(jù)、實(shí)驗(yàn)數(shù)據(jù)得到.當(dāng)提取軸承的振動(dòng)信號(hào)X以后,經(jīng)過(guò)運(yùn)算Y=XR,得到Y(jié)值,通過(guò)Y值的特征來(lái)進(jìn)行軸承故障的監(jiān)測(cè)與診斷.這種方法適用于軸承的早期故障診斷,并且有較好的效果.

    神經(jīng)網(wǎng)絡(luò)方法:以神經(jīng)元為信息處理的基本單元,信息通道為神經(jīng)元之間的連接弧,由多個(gè)神經(jīng)元共同組成網(wǎng)絡(luò)結(jié)構(gòu),特點(diǎn)是具有聯(lián)想、記憶和學(xué)習(xí)的功能.在機(jī)械故障診斷中,BP神經(jīng)網(wǎng)絡(luò)及其改進(jìn)算法應(yīng)用較為廣泛.BP算法的基本思想是把相似模式的特征值提取出來(lái),并映射到連接權(quán)值上,這是網(wǎng)絡(luò)學(xué)習(xí)的一個(gè)過(guò)程.當(dāng)網(wǎng)絡(luò)遇到一個(gè)新模式時(shí),將其特征值與所學(xué)習(xí)到的各類特征值相比較,從而判斷其類別.由于神經(jīng)網(wǎng)絡(luò)的函數(shù)具有逼近和記憶的處理能力,其在處理軸承振動(dòng)信號(hào)時(shí),將振動(dòng)的特征參數(shù)轉(zhuǎn)化為特征向量,進(jìn)而實(shí)現(xiàn)對(duì)滾動(dòng)軸承的故障監(jiān)測(cè)與診斷.BP的改進(jìn)算法優(yōu)點(diǎn)是網(wǎng)絡(luò)學(xué)習(xí)的速度較快,并且很大程度避免發(fā)生陷入局部極小的情況.

    遺傳算法(GA)是一種模擬自然界中自然選擇和種群遺傳的算法,具有隨機(jī)、迭代及進(jìn)化等特點(diǎn).這種算法可使問(wèn)題的解不斷優(yōu)化,并逼近最優(yōu)的值,適用于處理復(fù)雜非線性問(wèn)題,在信息處理中得到廣泛的應(yīng)用.此外,遺傳算法的改進(jìn)方法使得收斂更快,并避免陷入局部極小.

    灰色理論是由華中理工大學(xué)教授鄧聚龍?jiān)?982年提出的,陳家驥等基于此提出了灰色系統(tǒng)關(guān)聯(lián)診斷法,其特點(diǎn)為不需要大的樣本量,不受振動(dòng)信號(hào)參數(shù)分布影響,計(jì)算量較小,結(jié)果與定性結(jié)論一致等.通過(guò)建立故障與正常的模型,將實(shí)測(cè)模型與其相關(guān)聯(lián),找出關(guān)聯(lián)程度最大的模型,即可判斷是否發(fā)生故障.

    奇異值分解方法[16-18]:在軸承損傷處于早期或比較輕微時(shí),軸承工作時(shí)受到的沖擊并不明顯,由其引起的調(diào)幅現(xiàn)象也是十分微弱的,這樣就被軸承其它振動(dòng)頻率以及參入的噪聲淹沒(méi)了.這時(shí)如果使用包絡(luò)法,其處理的振動(dòng)信號(hào)的效果并不明顯,而奇異值分解方法則解決了這一問(wèn)題.同時(shí),奇異值分解方法還可以根據(jù)其處理得到的完整調(diào)幅信號(hào),得出載波和被調(diào)制的低頻信號(hào)頻率;SUV消噪方法具有零相移、波形失真小、信噪比高等優(yōu)點(diǎn);奇異值分解方法使處理結(jié)果調(diào)幅特征十分明顯,要好于帶通濾波器的效果.

    Teager能量算子是一種非線性信號(hào)處理方法,實(shí)質(zhì)是運(yùn)用了差分運(yùn)算.這種算法計(jì)算較為簡(jiǎn)單,且時(shí)間分辨率高.基于Teager能量算子,王天金等提出了一種的頻譜分析方法[19],具體應(yīng)用步驟:第一步,運(yùn)用Teager能量算子進(jìn)行運(yùn)算;第二步,對(duì)瞬時(shí)Teager能量序列進(jìn)行傅里葉變換,獲得Teager能量譜;第三步,將Teager能量譜中的主導(dǎo)頻率與已知軸承各種故障的頻率進(jìn)行對(duì)比,從而得出故障原因.

    隱馬爾可夫模型(HMM)作為一種統(tǒng)計(jì)分析模型,創(chuàng)立于20世紀(jì)70年代.80年代得以傳播、發(fā)展,成為信號(hào)處理的一個(gè)重要研究方向.有些學(xué)者將其運(yùn)用在軸承故障狀態(tài)辨識(shí)中,也取得了較好的效果.

    沖擊脈沖法(SPM)是根據(jù)故障軸承產(chǎn)生的低頻沖擊脈沖信號(hào)激起SPM傳感器的共振,經(jīng)過(guò)窄帶濾波器和脈沖形成電路得到一定頻段的脈沖序列.SPM方法就是依據(jù)這些脈沖序列來(lái)對(duì)軸承進(jìn)行故障的監(jiān)測(cè)與診斷,這種方法被認(rèn)為是一種有效的軸承故障診斷方法,但當(dāng)背景噪聲很強(qiáng)或有其它的沖擊源時(shí),SPM方法的判斷能力就會(huì)下降,失去了診斷能力.

    集合經(jīng)驗(yàn)?zāi)7纸?EEMD)方法是針對(duì)EMD方法的不足,提出的一種噪聲輔助數(shù)據(jù)分析方法,有效解決了EMD方法的混頻現(xiàn)象.蘇文勝提出了基于EMD降噪和譜峭度法的滾動(dòng)軸承振動(dòng)信號(hào)處理方法[20],突出了振動(dòng)信號(hào)的特征成分,有利于軸承故障的診斷.同時(shí)提出了一種基于小波包樣本熵的滾動(dòng)軸承故障診斷方法,通過(guò)計(jì)算,得出該算法可以較好地預(yù)測(cè)軸承的運(yùn)行狀態(tài).還提出了一種基于最優(yōu)Morlet小波濾波器和自相關(guān)增強(qiáng)算法,并通過(guò)實(shí)驗(yàn)證明了這種信號(hào)處理方法比基于EMD降噪和譜峭度法的滾動(dòng)軸承振動(dòng)信號(hào)處理方法能診斷更早期的軸承故障.

    Kalman濾波的信號(hào)處理方法,計(jì)算量較小,可以用于非平穩(wěn)狀態(tài).通過(guò)仿真模擬和實(shí)驗(yàn),可以看出Kalman濾波的降噪功能非常好,得出的譜線十分光滑,有用的特征信息十分明顯.

    S Abdusslam等[21]人采用了一種小波變換、包絡(luò)分析與時(shí)間編碼信號(hào)處理識(shí)別(TESPAR)相結(jié)合的方法進(jìn)行軸承的輕微故障診斷,這種方法有很好的故障識(shí)別能力.

    LIU Qing-jie等人比較了消除多項(xiàng)式趨勢(shì)法和五點(diǎn)三次平滑方法[22].實(shí)驗(yàn)結(jié)果表明,前者效果不明顯,而后者消除了高頻噪聲,同時(shí)信號(hào)曲線變得光滑.

    Li Changning等[23]提出一種基于α穩(wěn)態(tài)分布的新的非高斯統(tǒng)計(jì)模型,可以精確描述軸承故障信號(hào)的統(tǒng)計(jì)特征.

    Wei zhou等[24]描述了維納濾波信號(hào)處理方法在軸承振動(dòng)信號(hào)處理中的應(yīng)用,指出其具有信號(hào)的降噪功能,不受軸承參數(shù)的影響.Yang Ming等[25]提出了一種基于循環(huán)維納濾波和包絡(luò)譜的信號(hào)處理方法,并通過(guò)仿真與實(shí)驗(yàn)真實(shí)數(shù)據(jù)證明了方法的有效性.

    Li Meng等[26]研究了支持向量機(jī)(SVM)在軸承故障診斷中的應(yīng)用.SVM對(duì)小訓(xùn)練樣本有很好的分類性能,診斷精度對(duì)內(nèi)核函數(shù)和參量依賴小,適合于工程應(yīng)用.SVM與主成分分析(PCA)法的結(jié)合可以有效地提取故障特征,準(zhǔn)確識(shí)別故障模式.Wu Bin等[27]提出了SVM與HMM相結(jié)合的方法,并驗(yàn)證了該方法的有效性.Du Yi等[28]提出了一種概率盒與登普斯特謝佛結(jié)構(gòu)(DSS)相結(jié)合方法,并利用SVM方法減少計(jì)算量,這種方法具有快速、對(duì)噪聲不敏感、識(shí)別率高和準(zhǔn)確度高等優(yōu)點(diǎn).

    此外,分形盒維數(shù)方法、維譜分析法及其相關(guān)方法等也被不少學(xué)者應(yīng)用于軸承振動(dòng)信號(hào)的處理中,并取得了一些成效.

    3 結(jié)論

    3.1 現(xiàn)有信號(hào)處理技術(shù)的綜合運(yùn)用

    到目前為止,軸承振動(dòng)信號(hào)故障診斷方法已發(fā)展到可觀的數(shù)量,每一種信號(hào)處理方法都有其優(yōu)點(diǎn),但也存在其局限性.例如,峰值、峭度等參數(shù)對(duì)沖擊信號(hào)敏感,但不能分辨出沖擊信號(hào)的來(lái)源,可能有盲目性,產(chǎn)生誤判.有些信號(hào)處理方法較為準(zhǔn)確,但算法復(fù)雜,計(jì)算量很大,對(duì)設(shè)備的要求較高.有些信號(hào)處理方法對(duì)工作人員的經(jīng)驗(yàn)水平和操作能力要求較高,有一定的先驗(yàn)性,如譜分析法.此外,新出現(xiàn)的智能信號(hào)處理方法理論不夠完善,有待于進(jìn)行改進(jìn)和優(yōu)化.

    既然各種信號(hào)處理方法各有優(yōu)缺點(diǎn),如果在它們之中進(jìn)行選取,結(jié)合起來(lái)綜合運(yùn)用,會(huì)帶來(lái)好的效果.可以先采用時(shí)域法進(jìn)行信號(hào)的預(yù)處理,提取出有用信息,如特征參數(shù),故障特點(diǎn)等;再結(jié)合預(yù)處理結(jié)果通過(guò)頻域法或時(shí)頻法進(jìn)行精確診斷,會(huì)使結(jié)果較為準(zhǔn)確.有的信號(hào)處理計(jì)算量小,有的方法分析準(zhǔn)確,但計(jì)算量大,可以用計(jì)算量小的方法進(jìn)行計(jì)算處理,再用分析準(zhǔn)確的方法進(jìn)行分析,使得方法各自發(fā)揮其優(yōu)點(diǎn),互相彌補(bǔ)不足,從而實(shí)現(xiàn)軸承故障監(jiān)測(cè)診斷又好又快.

    3.2 信號(hào)處理技術(shù)的發(fā)展趨勢(shì)

    軸承振動(dòng)信號(hào)的處理不是一個(gè)簡(jiǎn)單的問(wèn)題,振動(dòng)過(guò)程中的噪聲參入,振動(dòng)狀態(tài)的非線性,非高斯性,非平穩(wěn)性等,都會(huì)給軸承振動(dòng)信號(hào)的處理帶來(lái)困難.如何監(jiān)測(cè)軸承的運(yùn)行狀態(tài),準(zhǔn)確找出軸承的故障的根源,是一個(gè)必須解決的問(wèn)題.隨著計(jì)算機(jī)技術(shù)及虛擬網(wǎng)絡(luò)技術(shù)的發(fā)展,診斷儀器硬件、軟件診斷設(shè)備的日趨精良,技術(shù)人員水平的不斷提高,軸承故障監(jiān)測(cè)與診斷技術(shù)將不斷發(fā)展、進(jìn)步.軸承故障診斷信號(hào)處理手段的可能發(fā)展趨勢(shì)歸納如下:

    (1) 多種信號(hào)處理技術(shù)的完善與相互結(jié)合應(yīng)用.單一的信號(hào)處理方法并不十分完善,可靠性有待提高,信號(hào)的處理過(guò)程也是有待于改善的.多種信號(hào)處理方法的結(jié)合,克服單一方法的不足,以提高信號(hào)的診斷效率與準(zhǔn)確率.

    (2) 智能的信號(hào)處理方法將會(huì)極大地發(fā)展.智能化是當(dāng)今時(shí)代的主流,這一現(xiàn)象也會(huì)反映在軸承的振動(dòng)信號(hào)處理方法方面.計(jì)算機(jī)及網(wǎng)絡(luò)的發(fā)展,虛擬神經(jīng)系統(tǒng)的開(kāi)發(fā)與完善,人工智能技術(shù)的進(jìn)步,開(kāi)拓了信號(hào)處理的另一個(gè)領(lǐng)域.自動(dòng)識(shí)別信號(hào)的特征,對(duì)其進(jìn)行分析、歸納與整理,加以判斷等功能將會(huì)不斷發(fā)展,降低技術(shù)人員的勞動(dòng)強(qiáng)度.這是與計(jì)算機(jī)的發(fā)展密不可分的.

    (3) 信號(hào)處理技術(shù)在工程中的應(yīng)用.所有的技術(shù)理論發(fā)展都會(huì)應(yīng)用于工程中來(lái).理論結(jié)合工程實(shí)際,將會(huì)引領(lǐng)信號(hào)處理技術(shù)的發(fā)展方向.工程問(wèn)題一般較為復(fù)雜,影響因素多,也有其固有的特點(diǎn)與規(guī)律,熟悉與掌握了這些特點(diǎn)與規(guī)律,將會(huì)使信號(hào)處理方法更加有效,更加有實(shí)用性.

    (4) 新的信號(hào)處理方法不斷出現(xiàn),其它理論的結(jié)合應(yīng)用.信號(hào)處理的方法趨于多樣化,復(fù)合化.其它的理論正逐漸應(yīng)用于信號(hào)的處理之中,同時(shí),又會(huì)有新的理論提出并應(yīng)用于信號(hào)的處理.這將會(huì)給信號(hào)處理技術(shù)注入新鮮的血液,有利于這一領(lǐng)域快速地發(fā)展.

    [1] GAO Q,DUAN C,FAN H.Rotating machine fault diagnosis using empirical mode decomposition[J].Mechanical Systems and Signal Processing,2008,22:1072-1081.

    [2] WANG Jiangping,HU Hongtao.Vibration-based fault diagnosis of pump using fuzzy technique[J].Measurement,2006,39:176-185.

    [3] Sawalhi N,Randall R B,Endo H.The enhancement of fault detection and diagnosis in rolling element bearings using minimum entropy deconvolution combined with spectral kurtosis [J].Mechanical Systems and Signal Processing 2007,(21):2616-2633.

    [4] Antoni Jeromef,The spectral kurtosis.A useful tool for characterizing non-stationary signals[J].Mechanical System and Signal Processing,2006,20:282-307.

    [5] Antoni Jeromef,Randall R B.The spectral kurtosis:application to the vibratory surveillance and diagnosis of rotationmachines[J].Mechanical System and Signal Processing,2006,20:308-331.

    [6] 李昊.基于小波循環(huán)自相關(guān)分析方法的滾動(dòng)軸承振動(dòng)信號(hào)處理技術(shù)的研究[D].長(zhǎng)沙:中南大學(xué),2005.

    LI Hao.Method of rolling bearing based on wavelet cyclic autocorrelation analysis of vibration signal processing research [D].Changsha:Central South University,2007.

    [7] 陳曉平,和衛(wèi)星,陳季云.基于小波分析的滾動(dòng)軸承振動(dòng)信號(hào)處理[J].儀器儀表學(xué)報(bào),2008,29(4):197-200.

    CHEN Xiaoping,HE Weixing,CHEN Jiyun.Wavelet based analysis for rolling bearing vibration signal processing [J].Chinese Journal of Scientific Instrument,2008,29(4):197-200.

    [8] 李毓輝,盧如平,方丹,等.故障滾動(dòng)軸承振動(dòng)信號(hào)的半盲處理方法研究[J].振動(dòng)與沖擊,2009,28(12):80-83.

    LI Yuhui,LU Ruping,FANG Dan,et al.Research on semi-blindsignal processing method for faulty rolling bearing vibration signal [J].Journal of vibration and shock,2009,28(12):80-83.

    [9] Nan Chyuan Tsai,Yueh Hsun King,Rong Mao Lee.Fault diagnosis for magnetic bearing systems[J].Mechanical Systems and Signal Processing,2009,23:1339-1351.

    [10] Alonsoa F J,Salgado D R.Analysis of the structure of vibration signals for tool wear detection [J].Mechanical Systems and Signal Processing,2008,22:735-748.

    [11] Daubechies I.The wavelet transform time-frequency localization and signal analysis[J].IEEE Trans on Information Theory,1990,36(5):961-1006.

    [12] Jonathan S Smith.The local mean decomposition and its application to EEG perception data [J].Journal of The Royal Society Interface,2005,2:433-454.

    [13] 徐繼剛.滾動(dòng)軸承振動(dòng)信號(hào)處理與故障識(shí)別方法研究[D].蘭州:蘭州理工大學(xué),2012.

    XU Jigang.The research of rolling bearing vibration signal processing and fault recognition method [D].Lanzhou:Lanzhou University of Technology,2012.

    [14] 唐國(guó)明,楊曙年.基于虛擬儀器的軸承振動(dòng)信號(hào)處理技術(shù)[J].軸承,2005(12):24-27.

    TANG Guoming,YANG Shunian.Bearing vibration signal processing technique based on virtual instrument [J].Bearing,2005(12):24-27.

    [15] 郭曉鑫,劉立,馮志鵬.數(shù)學(xué)形態(tài)學(xué)在信號(hào)處理中的應(yīng)用綜述[J].儀器儀表學(xué)報(bào),2009,30(10):187-191.

    GUO Xiaoxin,LIU Li,FENG Zhipeng.Mathematical morphology review and its applicaitons in signal porcessing [J].Chinese Journal of Scientific Instrument,2009,30(10):187-191.

    [16] Lehtola L,Karsikas M,Koskinen M,et al.Effects of noise and filtering on SVD-based morphological parameters of the T wave in the ECG[J].Journal of Medical Engineering and Technology,2008,32(5):400- 407.

    [17] Mark B,Jiang D.Noise reduction in multiple-echo data sets using singular value decomposition[J].Magnetic Resonance Imaging,2006,24(2):849- 856.

    [18] Cempel C.Generalized singular value decomposition in multidimensional condition monitoring of machines- A proposal of comparative diagnostics[J].Mechanical Systems and Signal Processing,2009,23(3):701- 711.

    [19] 王天金,馮志鵬,郝如江,等.基于Teager能量算子的滾動(dòng)軸承故障診斷研究[J].振動(dòng)與沖擊,2012,31(2):1-6.

    WANG Tianjin,FENG Zhipeng,HAO Rujiang,et al.Fault diagnosis of rolling element bearings based on Teager energy operator [J].Journal of vibration and shock,2012,31(2):1-6.

    [20] 蘇文勝.滾動(dòng)軸承振動(dòng)信號(hào)處理及特征提取方法研究[D].大連:大連理工大學(xué),2010.

    SU Wensheng.Research on rolling element bearing vibration signal processing and feature extraction method [D].Dalian:Dalian University of Technology,2010.

    [21] ABDUSSLAM S,RAHARJO P,GU P,et al.Bearing defect detection and diagnosis using a time encoded signal processing and pattern recognition method[C]∥Comadem2012,25th International Congress on Condition Monitoring and Diagnostic Engineering.Bristol,UK:Institute of Physics Publishing,2012:1-12.

    [22] LIU Qingjie,LIU Xiaofang,ZHAN Jun.Study on pre-processing methods of precision bearing vibration signal[C]∥Esiat,2009 International Conference on Environmental Science and Information Application Technology.Piscataway,NJ:IEEE Press,2009:572-574.

    [23] LI Changning,YU Gang.A new statistical modelfor rolling element bearing fault signal based on alpha-stable distribution[C]∥Iccms2010,Proceedings of the 2010 Second International Conference on Computer Modeling and Simulation.Piscataway,NJ:IEEE Press,2010:386-390.

    [24] ZHOU W,HABETLER T G,HARLEY R G,et al.Incipient bearing fault detection via stator current noise cancellation using Wiener filter[C]∥Sdemped,IEEE Symposium on Diagnostics for Electric Machines Power Electronics and Drives.Piscataway,NJ:IEEE Press,2007:11-16.

    [25] YANG Ming,JIN Chen,DONG Guang-ming.Weak fault feature extraction of rolling bearing based on cyclic Wiener filter and envelopespectrum[J],Mechanical Systems and Signal Processing,2011,25(5):1773-1785.

    [26] LI Meng,WANG Miao,WANG Changchun.Research on SVM classification performance in rolling bearing diagnosis[C]∥Icicta2010,Proceedings of the 2010 International Conference on Intelligent Computation Technology and Automation.Los Alamitos,CA:IEEE Computer Society,2010:132-135.

    [27] WU Bin,YU Shanping,LUO Yuegang,et al.Rolling bearing faults diagnosis method based on SVM-HMM[C]∥Icmtma2010,Proceedings of the 2010 International Conference on Measuring Technology and Mechatronics Automation.Piscataway,NJ:IEEE Press,2010:295-298.

    [28] DU Y,CHI Y L,WU X.Bearing fault diagnosis based on probability boxes theory and SVM method[J].Applied Mechanics and Materials,2011,79:93-98.

    Review on signal processing for rolling bearing vibrations

    HU Zhi-yong,HU Jie-xin,XIE Li-yang,TONG An-shi,ZHANG Xiao-chen

    (1.Northeastern University,Key Laboratory of Vibration and Control of Aero-Propulsion Systems,Shenyang 110819,China; 2.Northeastern University,Institute of Modern Design and Analysis,Shenyang 110819,China)

    Firstly, various methods for bearing vibration signal processing were summarized including the time-domain and frequency-domain analyses. Then, applications of time-domain parameters, along with advantages and disadvantages of frequency domain analysis methods, are then pinpointed. Finally, it is found that the time-frequency analysis method reflects local characteristics of non-stationary signals. Specifically for bearing vibration signal processing, the method combination and new technology application ride on the future trends.

    rolling bearings; signal processing; vibration; time domain; frequency domain

    中央高校基本科研業(yè)務(wù)費(fèi)項(xiàng)目“復(fù)雜裝備關(guān)鍵子系統(tǒng)冗余配置與維修策略的綜合決策研究”(課題編號(hào):N150304013),遼寧省博士啟動(dòng)基金(項(xiàng)目編號(hào):201501154)

    胡智勇(1976-),男,博士,講師.E-mail:neujasonhu@163.com

    TH 113.1; TH 165+.3

    A

    1672-5581(2016)06-0525-07

    猜你喜歡
    故障診斷振動(dòng)故障
    振動(dòng)的思考
    振動(dòng)與頻率
    故障一點(diǎn)通
    中立型Emden-Fowler微分方程的振動(dòng)性
    奔馳R320車ABS、ESP故障燈異常點(diǎn)亮
    因果圖定性分析法及其在故障診斷中的應(yīng)用
    故障一點(diǎn)通
    江淮車故障3例
    UF6振動(dòng)激發(fā)態(tài)分子的振動(dòng)-振動(dòng)馳豫
    基于LCD和排列熵的滾動(dòng)軸承故障診斷
    日本wwww免费看| 少妇精品久久久久久久| 亚洲全国av大片| 亚洲国产毛片av蜜桃av| 亚洲av美国av| 99国产精品免费福利视频| 久久精品国产亚洲av香蕉五月 | 国产人伦9x9x在线观看| 午夜日韩欧美国产| 国产在视频线精品| 狂野欧美激情性bbbbbb| 少妇粗大呻吟视频| 国产精品影院久久| 国产av一区二区精品久久| 久久久久精品人妻al黑| 精品福利永久在线观看| 国产黄频视频在线观看| 中文字幕色久视频| 国产1区2区3区精品| 国产精品久久久久久精品电影小说| 侵犯人妻中文字幕一二三四区| 18禁观看日本| 巨乳人妻的诱惑在线观看| 中文字幕人妻丝袜一区二区| 黄频高清免费视频| 久久九九热精品免费| 亚洲欧美成人综合另类久久久| 欧美黄色淫秽网站| 正在播放国产对白刺激| 精品少妇内射三级| 在线 av 中文字幕| 12—13女人毛片做爰片一| 国产日韩欧美亚洲二区| 久久狼人影院| 人人妻人人爽人人添夜夜欢视频| 亚洲男人天堂网一区| 成人亚洲精品一区在线观看| 久久女婷五月综合色啪小说| 久久国产精品大桥未久av| 黄色视频,在线免费观看| 国产在线视频一区二区| 日韩欧美一区二区三区在线观看 | 一本综合久久免费| 十八禁高潮呻吟视频| 国产精品一区二区在线观看99| 男人操女人黄网站| 精品高清国产在线一区| 精品国产超薄肉色丝袜足j| 久久久久久人人人人人| 美女高潮到喷水免费观看| 一个人免费在线观看的高清视频 | 国产亚洲av高清不卡| 久久精品aⅴ一区二区三区四区| 在线看a的网站| 中文精品一卡2卡3卡4更新| 又黄又粗又硬又大视频| av欧美777| 男女免费视频国产| 欧美精品av麻豆av| 动漫黄色视频在线观看| 久久久久久免费高清国产稀缺| 91字幕亚洲| 岛国毛片在线播放| 18禁黄网站禁片午夜丰满| av在线播放精品| 悠悠久久av| 亚洲成av片中文字幕在线观看| a在线观看视频网站| 久久国产精品人妻蜜桃| 18禁观看日本| 亚洲欧美精品自产自拍| 高清av免费在线| 男人爽女人下面视频在线观看| tube8黄色片| 精品国产乱码久久久久久男人| 亚洲免费av在线视频| 一本久久精品| 国产av一区二区精品久久| 视频在线观看一区二区三区| 91老司机精品| 国产在线一区二区三区精| 国产免费现黄频在线看| 大码成人一级视频| 下体分泌物呈黄色| 一本一本久久a久久精品综合妖精| 国产色视频综合| 黄片播放在线免费| 亚洲av成人不卡在线观看播放网 | 丝袜脚勾引网站| 日韩精品免费视频一区二区三区| 又大又爽又粗| 一级毛片精品| 丝瓜视频免费看黄片| 国产成人精品在线电影| 欧美另类亚洲清纯唯美| 一二三四在线观看免费中文在| 久热这里只有精品99| 久久久久久久久久久久大奶| 午夜免费鲁丝| 男女午夜视频在线观看| 国产精品 国内视频| 欧美大码av| 激情视频va一区二区三区| 啦啦啦 在线观看视频| 国产极品粉嫩免费观看在线| 91九色精品人成在线观看| 极品少妇高潮喷水抽搐| 亚洲国产欧美一区二区综合| 国产免费视频播放在线视频| 天天影视国产精品| 欧美精品一区二区大全| 日本av手机在线免费观看| 午夜视频精品福利| 亚洲一码二码三码区别大吗| 9色porny在线观看| av在线app专区| 51午夜福利影视在线观看| 亚洲avbb在线观看| 美女视频免费永久观看网站| 免费观看人在逋| 欧美另类亚洲清纯唯美| 国产高清视频在线播放一区 | 亚洲午夜精品一区,二区,三区| 久久毛片免费看一区二区三区| 成年人黄色毛片网站| 日日夜夜操网爽| 人人妻,人人澡人人爽秒播| 丝袜美腿诱惑在线| 国产男女内射视频| 成人国产一区最新在线观看| 老鸭窝网址在线观看| 国产一区二区 视频在线| 久久久久久人人人人人| 欧美精品一区二区大全| 91精品国产国语对白视频| 午夜免费成人在线视频| 啪啪无遮挡十八禁网站| 91老司机精品| 久久国产精品男人的天堂亚洲| 啪啪无遮挡十八禁网站| 免费在线观看日本一区| 天天影视国产精品| 别揉我奶头~嗯~啊~动态视频 | 视频在线观看一区二区三区| 久久青草综合色| av天堂在线播放| av网站免费在线观看视频| 成人18禁高潮啪啪吃奶动态图| 国产无遮挡羞羞视频在线观看| 久久久水蜜桃国产精品网| 亚洲伊人久久精品综合| 天堂俺去俺来也www色官网| 岛国毛片在线播放| 男女国产视频网站| 亚洲伊人色综图| 51午夜福利影视在线观看| 国产高清国产精品国产三级| 嫩草影视91久久| 日本a在线网址| 亚洲欧洲精品一区二区精品久久久| 超碰97精品在线观看| 一个人免费在线观看的高清视频 | 丁香六月天网| 欧美日韩成人在线一区二区| 午夜福利影视在线免费观看| 母亲3免费完整高清在线观看| 最新在线观看一区二区三区| 99热全是精品| 女人久久www免费人成看片| 十分钟在线观看高清视频www| 亚洲欧美精品综合一区二区三区| 最近中文字幕2019免费版| 国产成人啪精品午夜网站| 性色av一级| 91成年电影在线观看| 亚洲精品国产一区二区精华液| 国产亚洲一区二区精品| 99精国产麻豆久久婷婷| 两个人看的免费小视频| 在线观看一区二区三区激情| 国产av又大| 女人爽到高潮嗷嗷叫在线视频| 成年动漫av网址| 制服诱惑二区| 曰老女人黄片| 在线精品无人区一区二区三| 一级片'在线观看视频| 韩国高清视频一区二区三区| 亚洲性夜色夜夜综合| 欧美中文综合在线视频| 国产在线一区二区三区精| 国产精品 国内视频| 日日爽夜夜爽网站| 另类亚洲欧美激情| 成人国产av品久久久| 久久久久久久国产电影| 最近最新免费中文字幕在线| 无限看片的www在线观看| 久久久欧美国产精品| 久久这里只有精品19| 老熟妇乱子伦视频在线观看 | 久久综合国产亚洲精品| 乱人伦中国视频| 电影成人av| 中文字幕高清在线视频| 人人妻人人爽人人添夜夜欢视频| 国产男女内射视频| 国产亚洲av片在线观看秒播厂| 最近最新免费中文字幕在线| 热re99久久国产66热| 在线永久观看黄色视频| 国产精品亚洲av一区麻豆| 一边摸一边做爽爽视频免费| 日韩三级视频一区二区三区| 老司机福利观看| 亚洲精品中文字幕在线视频| 男女边摸边吃奶| 99久久99久久久精品蜜桃| 考比视频在线观看| 一本色道久久久久久精品综合| 一区二区三区激情视频| 精品福利永久在线观看| 9191精品国产免费久久| 亚洲久久久国产精品| 亚洲成国产人片在线观看| 一区二区三区激情视频| 亚洲欧美一区二区三区黑人| 亚洲av片天天在线观看| 亚洲国产日韩一区二区| 亚洲国产精品一区三区| 在线观看免费视频网站a站| 国产成人一区二区三区免费视频网站| 国产亚洲一区二区精品| 久久精品国产亚洲av高清一级| 91老司机精品| 少妇被粗大的猛进出69影院| 午夜激情久久久久久久| 亚洲精品国产av蜜桃| 国产精品 国内视频| 9191精品国产免费久久| 亚洲精品国产精品久久久不卡| 一级片'在线观看视频| 成年av动漫网址| 亚洲精品自拍成人| 老司机亚洲免费影院| 99国产精品免费福利视频| 自拍欧美九色日韩亚洲蝌蚪91| 老司机靠b影院| 国产成人精品无人区| 中文字幕高清在线视频| 在线天堂中文资源库| 丰满人妻熟妇乱又伦精品不卡| 乱人伦中国视频| 国产精品偷伦视频观看了| 高清欧美精品videossex| 日韩,欧美,国产一区二区三区| 2018国产大陆天天弄谢| 国产精品久久久av美女十八| 麻豆国产av国片精品| 老司机在亚洲福利影院| 精品免费久久久久久久清纯 | 国产精品一区二区在线不卡| 日韩 欧美 亚洲 中文字幕| 日本一区二区免费在线视频| 亚洲精品中文字幕在线视频| 桃红色精品国产亚洲av| 不卡av一区二区三区| 五月开心婷婷网| 精品一区在线观看国产| 一区二区日韩欧美中文字幕| 一级片'在线观看视频| 乱人伦中国视频| 国产一区二区三区av在线| 美女国产高潮福利片在线看| 欧美激情高清一区二区三区| 欧美日韩国产mv在线观看视频| 久久亚洲国产成人精品v| 操美女的视频在线观看| 国产高清视频在线播放一区 | 麻豆乱淫一区二区| 高清黄色对白视频在线免费看| 国产av国产精品国产| 80岁老熟妇乱子伦牲交| 精品国产一区二区久久| 悠悠久久av| 一级,二级,三级黄色视频| 国产免费一区二区三区四区乱码| 99国产精品免费福利视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲国产成人一精品久久久| 丝袜人妻中文字幕| 在线 av 中文字幕| 精品国产乱码久久久久久男人| 人人妻人人添人人爽欧美一区卜| 欧美变态另类bdsm刘玥| 国产精品免费大片| avwww免费| 日韩大片免费观看网站| 91成人精品电影| 欧美日韩一级在线毛片| 国产高清视频在线播放一区 | av天堂在线播放| 99精品久久久久人妻精品| 日韩中文字幕欧美一区二区| 桃红色精品国产亚洲av| 国产精品国产三级国产专区5o| 纵有疾风起免费观看全集完整版| 亚洲精品国产av成人精品| 亚洲精品国产精品久久久不卡| 久久国产精品人妻蜜桃| 国产伦理片在线播放av一区| 亚洲国产精品一区三区| 午夜福利,免费看| 夫妻午夜视频| avwww免费| 精品久久蜜臀av无| tube8黄色片| 99久久精品国产亚洲精品| 狠狠精品人妻久久久久久综合| 成年人黄色毛片网站| 亚洲欧美日韩另类电影网站| 久久精品成人免费网站| 亚洲精品一卡2卡三卡4卡5卡 | 国产亚洲精品一区二区www | 国产免费现黄频在线看| 热re99久久国产66热| 19禁男女啪啪无遮挡网站| 人人澡人人妻人| 宅男免费午夜| 久久ye,这里只有精品| 久久狼人影院| 国产亚洲午夜精品一区二区久久| 日韩电影二区| 熟女少妇亚洲综合色aaa.| 亚洲成人免费av在线播放| 少妇人妻久久综合中文| av天堂在线播放| 最黄视频免费看| 国产亚洲一区二区精品| 国产精品成人在线| 欧美亚洲日本最大视频资源| 中文字幕人妻熟女乱码| 丝袜在线中文字幕| 99国产综合亚洲精品| 人人妻人人澡人人看| 老鸭窝网址在线观看| 久久久久久亚洲精品国产蜜桃av| av电影中文网址| 美女高潮喷水抽搐中文字幕| 久久国产精品人妻蜜桃| 99国产精品免费福利视频| 黄片播放在线免费| 中文字幕人妻丝袜一区二区| 亚洲 欧美一区二区三区| 国产精品二区激情视频| 久久热在线av| 日本精品一区二区三区蜜桃| av在线播放精品| 日韩,欧美,国产一区二区三区| 黄片播放在线免费| 19禁男女啪啪无遮挡网站| 久久青草综合色| 捣出白浆h1v1| 精品少妇内射三级| 精品国产乱码久久久久久男人| 一区二区三区乱码不卡18| 免费高清在线观看日韩| 一区二区三区乱码不卡18| 香蕉丝袜av| 高清欧美精品videossex| 狠狠婷婷综合久久久久久88av| 国产亚洲精品一区二区www | 少妇被粗大的猛进出69影院| 日本91视频免费播放| 欧美日韩福利视频一区二区| 91字幕亚洲| 亚洲欧美成人综合另类久久久| www.自偷自拍.com| 国产成人一区二区三区免费视频网站| av在线app专区| 欧美黄色淫秽网站| 国产淫语在线视频| 中文字幕最新亚洲高清| 中文字幕人妻丝袜一区二区| av在线播放精品| 性高湖久久久久久久久免费观看| 免费一级毛片在线播放高清视频 | 中文字幕人妻熟女乱码| 国产精品影院久久| av又黄又爽大尺度在线免费看| 精品一区二区三卡| 欧美国产精品一级二级三级| 美女高潮喷水抽搐中文字幕| 中文字幕av电影在线播放| 19禁男女啪啪无遮挡网站| 亚洲 欧美一区二区三区| 国产在线免费精品| 在线观看免费高清a一片| 欧美av亚洲av综合av国产av| 香蕉国产在线看| 国产无遮挡羞羞视频在线观看| 久9热在线精品视频| 亚洲精品国产av蜜桃| 精品久久久久久电影网| 波多野结衣av一区二区av| 精品久久蜜臀av无| 国产精品国产av在线观看| 国产一区二区 视频在线| 十八禁网站免费在线| 国产成人免费无遮挡视频| 亚洲av电影在线进入| 国产在视频线精品| 亚洲国产欧美一区二区综合| 啦啦啦啦在线视频资源| 热99re8久久精品国产| 国产野战对白在线观看| 色综合欧美亚洲国产小说| 日韩电影二区| 成人影院久久| 精品国产一区二区三区四区第35| av天堂久久9| 十八禁网站网址无遮挡| 午夜91福利影院| 精品国产乱码久久久久久男人| 女性被躁到高潮视频| 精品福利永久在线观看| 性色av一级| 久久精品久久久久久噜噜老黄| 黄色 视频免费看| 电影成人av| 丝袜美足系列| 老汉色av国产亚洲站长工具| 中国国产av一级| 久久av网站| 男女高潮啪啪啪动态图| 国产精品香港三级国产av潘金莲| 少妇的丰满在线观看| 黑人巨大精品欧美一区二区蜜桃| 十八禁网站网址无遮挡| 69精品国产乱码久久久| 又紧又爽又黄一区二区| 91av网站免费观看| 啦啦啦在线免费观看视频4| 日本精品一区二区三区蜜桃| 丰满迷人的少妇在线观看| videos熟女内射| 国产欧美日韩一区二区精品| 在线观看www视频免费| 美女大奶头黄色视频| 老司机深夜福利视频在线观看 | 国产av一区二区精品久久| 人成视频在线观看免费观看| 又黄又粗又硬又大视频| 精品久久久久久久毛片微露脸 | 91字幕亚洲| 麻豆国产av国片精品| 丰满饥渴人妻一区二区三| 久久天躁狠狠躁夜夜2o2o| 国产老妇伦熟女老妇高清| 国产亚洲精品久久久久5区| 亚洲国产精品成人久久小说| 一区二区三区精品91| 好男人电影高清在线观看| 国产欧美日韩综合在线一区二区| 蜜桃国产av成人99| 亚洲精品美女久久av网站| 最新的欧美精品一区二区| 国产精品免费视频内射| 欧美xxⅹ黑人| 久久精品国产综合久久久| 美女高潮到喷水免费观看| 黄色怎么调成土黄色| 国产亚洲欧美在线一区二区| 久久久久网色| 国产av又大| 国产野战对白在线观看| 中文字幕最新亚洲高清| 国产成人精品在线电影| 777久久人妻少妇嫩草av网站| 男人操女人黄网站| 少妇精品久久久久久久| 欧美久久黑人一区二区| 99精国产麻豆久久婷婷| 又紧又爽又黄一区二区| 国产成人一区二区三区免费视频网站| 久久狼人影院| 国产99久久九九免费精品| 欧美激情高清一区二区三区| 日韩一区二区三区影片| 日本vs欧美在线观看视频| 久热这里只有精品99| av在线播放精品| 精品免费久久久久久久清纯 | e午夜精品久久久久久久| 天天躁日日躁夜夜躁夜夜| 91精品伊人久久大香线蕉| 国产一区二区激情短视频 | 丝袜喷水一区| 色婷婷av一区二区三区视频| 久久人人97超碰香蕉20202| 女警被强在线播放| 飞空精品影院首页| 在线观看免费日韩欧美大片| av免费在线观看网站| 久久免费观看电影| 国产精品免费视频内射| 欧美亚洲 丝袜 人妻 在线| 美女视频免费永久观看网站| 最新的欧美精品一区二区| 九色亚洲精品在线播放| 亚洲精品久久久久久婷婷小说| 麻豆乱淫一区二区| 亚洲天堂av无毛| 在线亚洲精品国产二区图片欧美| 国产成人一区二区三区免费视频网站| 国产日韩欧美视频二区| 最近最新免费中文字幕在线| 性色av一级| 两人在一起打扑克的视频| 午夜免费观看性视频| 久久99热这里只频精品6学生| 天堂中文最新版在线下载| 热re99久久精品国产66热6| 性少妇av在线| 国产又色又爽无遮挡免| 在线观看舔阴道视频| 电影成人av| 自线自在国产av| 欧美另类一区| 日韩视频一区二区在线观看| 久热爱精品视频在线9| 亚洲欧美成人综合另类久久久| tube8黄色片| 国产精品免费大片| 久久亚洲精品不卡| 丰满迷人的少妇在线观看| 国产一区二区激情短视频 | 国产精品二区激情视频| 国产男女内射视频| 成人av一区二区三区在线看 | 亚洲av日韩精品久久久久久密| 欧美激情久久久久久爽电影 | 黄色毛片三级朝国网站| 亚洲国产欧美一区二区综合| 丝瓜视频免费看黄片| www.熟女人妻精品国产| 久久女婷五月综合色啪小说| 麻豆乱淫一区二区| 18禁裸乳无遮挡动漫免费视频| 国产一卡二卡三卡精品| 91精品国产国语对白视频| 伊人亚洲综合成人网| 丝瓜视频免费看黄片| 成人影院久久| 国产亚洲精品久久久久5区| 蜜桃国产av成人99| 亚洲avbb在线观看| 一个人免费在线观看的高清视频 | 99国产精品一区二区蜜桃av | 建设人人有责人人尽责人人享有的| 久久久久久亚洲精品国产蜜桃av| 99re6热这里在线精品视频| 高清欧美精品videossex| 汤姆久久久久久久影院中文字幕| 日本五十路高清| 亚洲国产精品一区三区| 国产免费福利视频在线观看| 午夜精品国产一区二区电影| 国产片内射在线| 99国产综合亚洲精品| 久久毛片免费看一区二区三区| 91国产中文字幕| 欧美成人午夜精品| 亚洲色图 男人天堂 中文字幕| 日韩欧美国产一区二区入口| 欧美一级毛片孕妇| 国产又色又爽无遮挡免| 黄色怎么调成土黄色| 女人久久www免费人成看片| 无限看片的www在线观看| 亚洲精华国产精华精| 国产精品久久久久久精品古装| 日本一区二区免费在线视频| tube8黄色片| 黑人猛操日本美女一级片| 久久av网站| 黄色毛片三级朝国网站| 91av网站免费观看| 国产精品99久久99久久久不卡| 香蕉国产在线看| 国产亚洲精品第一综合不卡| 大码成人一级视频| 欧美黑人精品巨大| 欧美精品av麻豆av| 美女高潮到喷水免费观看| 99久久综合免费| 亚洲性夜色夜夜综合| 欧美乱码精品一区二区三区| 国产精品久久久久久人妻精品电影 | 男人添女人高潮全过程视频| 亚洲欧洲日产国产| 在线天堂中文资源库| 免费在线观看视频国产中文字幕亚洲 | av国产精品久久久久影院| 99精国产麻豆久久婷婷| 水蜜桃什么品种好| 中文字幕人妻熟女乱码| 乱人伦中国视频| 亚洲欧洲日产国产| 在线天堂中文资源库| 亚洲精品国产色婷婷电影| 19禁男女啪啪无遮挡网站| 老司机深夜福利视频在线观看 | 黑人巨大精品欧美一区二区mp4| 精品一区二区三区四区五区乱码|