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

    融合Autogram的共振解調(diào)和1.5維譜的滾動軸承復(fù)合故障診斷方法

    2022-02-15 08:46:08王慧濱剡昌鋒孟佳東陳光億吳黎曉
    振動工程學(xué)報(bào) 2022年6期
    關(guān)鍵詞:峭度特征頻率內(nèi)圈

    王慧濱,剡昌鋒,孟佳東,3,陳光億,吳黎曉

    (1.蘭州理工大學(xué)機(jī)電工程學(xué)院,甘肅 蘭州 730050;2.漳州衛(wèi)生職業(yè)學(xué)院醫(yī)學(xué)技術(shù)系,福建 漳州 363000;3.蘭州交通大學(xué)鐵道技術(shù)學(xué)院,甘肅 蘭州 730030)

    引 言

    在旋轉(zhuǎn)機(jī)械裝備實(shí)際使用過程中,由于外部載荷作用和軸承性能退化,滾動軸承內(nèi)外圈和滾動體的磨損、剝落、裂紋等故障的出現(xiàn)具有隨機(jī)性、并發(fā)性、繼發(fā)性和隱蔽性等,多種或者多點(diǎn)故障發(fā)生并存的情況成為常態(tài),表現(xiàn)為故障特征信號相互耦合。滾動軸承復(fù)合故障的產(chǎn)生給故障確診帶來更大的困難,以往的單故障診斷方法很難適應(yīng)復(fù)合故障模式下的故障診斷問題,因此需要研究滾動軸承復(fù)合故障耦合機(jī)理,分析故障特征變化規(guī)律,提出故障特征分離與識別的方法,對提高故障診斷的準(zhǔn)確性具有重要意義。

    目前復(fù)合故障的診斷方法主要有兩類。第一類方法是將原始信號分解成多個(gè)不同分量信號,通過選取不同分量信號,提取和分析故障特征頻率,實(shí)現(xiàn)復(fù)合故障診斷。Wang 等[1]采用了盲源分離和經(jīng)驗(yàn)?zāi)B(tài)分解相結(jié)合的方法,成功診斷出滾動軸承復(fù)合故障。Tang 等[2]提出了一種基于變分模態(tài)分解(Variational Mode Decomposition,VMD)和盲源分離的滾動軸承復(fù)合故障診斷方法,有效地分離了外圈和滾子的故障。Wan 等[3]采用快速譜峭度和VMD 相結(jié)合的方法對復(fù)合故障進(jìn)行診斷。Miao等[4]提出了參數(shù)優(yōu)化VMD 的復(fù)合故障識別方法,并成功應(yīng)用到機(jī)車軸箱軸承復(fù)合故障診斷中。胡愛軍等[5]采用VMD 結(jié)合1.5 維譜的方法對滾動軸承復(fù)合故障特征進(jìn)行分離。

    采用基于分量信號分析的方法進(jìn)行復(fù)合故障識別,分量信號的選取直接影響診斷的效果。選取有效分量信號分析診斷時(shí),峭度指標(biāo)[6]、互相關(guān)系數(shù)[1]、相關(guān)峭度[5]、頻帶熵[7]、能量比[8]、包絡(luò)譜特征因子[9]等指標(biāo)均得到應(yīng)用。因?yàn)閺?fù)合故障信號的復(fù)雜性,僅根據(jù)單個(gè)指標(biāo)選取分量信號可能會出現(xiàn)分量選取不準(zhǔn)確的問題;同時(shí),復(fù)合故障信號分解后,分量信號有時(shí)不僅僅以單一故障成分存在,所以基于分析不同分量信號的復(fù)合故障診斷方法容易受到分量信號選取的影響,如果選取不準(zhǔn)確將導(dǎo)致診斷錯(cuò)誤或漏診[10]。

    第二類方法是基于解卷積的方法。McDonald等[11]提出了最大相關(guān)峭度解卷積(Maximum Correlated Kurtosis Deconvolution,MCKD)方法,利用已知故障周期從故障信號中解卷積出周期性脈沖成分,具有較好的濾波效果,在齒輪故障診斷中取得了較好的效果。胡愛軍等[12]提出了一種融合譜峭度和MCKD 的滾動軸承復(fù)合故障診斷方法,有效分離故障成分。萬書亭等[13]采用VMD 和MCKD 相結(jié)合的方法,有效分離滾動軸承內(nèi)圈和外圈的故障。

    盡管MCKD 能較好地權(quán)衡周期性脈沖占據(jù)的比重,在復(fù)合故障診斷中得到了應(yīng)用,但是該方法預(yù)設(shè)參數(shù)較多,濾波器長度和周期的選取決定著故障信號分離效果的好壞[14-15]。Tang 等[16]采用布谷鳥搜索算法選取最優(yōu)參數(shù),提出了自適應(yīng)MCKD 的方法,能有效地識別出不同故障類型。Lyu 等[17]采用量子遺傳算法,自適應(yīng)選擇MCKD 中濾波器長度和周期這兩個(gè)關(guān)鍵參數(shù),成功識別出齒輪軸承復(fù)合故障。Mc-Donald 等[18]又提出了多點(diǎn)最優(yōu)最小熵解卷積(Multipoint Optimal Minimum Entropy Deconvolution Adjusted,MOMEDA)方法,適用于周期性故障特征提取。王志堅(jiān)等[19]采用了多點(diǎn)峭度改進(jìn)MOMEDA 方法,能有效提取不同周期故障信號成分,并運(yùn)用到齒輪箱復(fù)合故障診斷中。Wang 等[20]通過優(yōu)化濾波器長度和故障周期范圍改進(jìn)MOMEDA 方法,有效識別了齒輪箱微弱復(fù)合故障。齊詠生等[21]采用自適應(yīng)稀疏共振分解對信號進(jìn)行預(yù)處理,利用MOMEDA進(jìn)行滾動軸承復(fù)合故障診斷。利用解卷積分離復(fù)合故障信號,需要預(yù)先設(shè)置合理的參數(shù)或進(jìn)行預(yù)處理,保證信號分離的準(zhǔn)確性。但是在強(qiáng)噪聲干擾的環(huán)境下,故障的瞬態(tài)沖擊成分容易被軸的轉(zhuǎn)頻和倍頻等諧波分量和噪聲干擾,出現(xiàn)漏診現(xiàn)象。

    共振解調(diào)是提取故障特征的有效方法之一。郭俊等[22]提出基于特征頻率強(qiáng)度選取最優(yōu)共振帶的滾動軸承早期故障診斷,但并沒有應(yīng)用于復(fù)合故障診斷中。Wan 等[23]采用MCKD 對原始信號濾波,結(jié)合快速譜峭度和共振解調(diào)方法,在軸承復(fù)合故障診斷中取得較好效果。該類方法[3,12,23]都是先分離信號,再對不同故障類型的信號分別選取共振帶,其本質(zhì)也是對單故障共振帶的選取。因此,如何在不分離故障信號的前提下有效選取復(fù)合故障的共振頻帶是目前研究的一個(gè)難點(diǎn)。

    針對先分離后診斷的復(fù)合故障識別方法存在的局限性,為提高復(fù)合故障識別的準(zhǔn)確性和直觀性,本文提出了基于Autogram 的共振解調(diào)和1.5 維譜的復(fù)合故障診斷方法,能夠在不分離復(fù)合故障信號的前提下識別故障類型?;谂帕徐?、峭度和相關(guān)系數(shù)三個(gè)指標(biāo),提出了一種新的綜合指標(biāo)Z選取VMD分解后本征模態(tài)函數(shù)(Intrinsic Mode Function,IMF)分量,進(jìn)行信號重構(gòu),提高信噪比,降低背景噪聲對診斷的影響;采用Autogram 算法確定共振頻帶中心頻率和帶寬,對共振信號進(jìn)行包絡(luò)解調(diào),得到包絡(luò)信號的1.5 維譜,根據(jù)1.5 維譜中的故障特征來識別滾動軸承復(fù)合故障,采用滾動軸承3 種不同類型復(fù)合故障的實(shí)測信號驗(yàn)證了所提方法的可行性。

    1 基于Autogram 的共振解調(diào)和1.5 維譜的滾動軸承復(fù)合故障診斷方法

    1.1 選取分量信號的綜合指標(biāo)Z

    排列熵是一種衡量一維時(shí)間序列復(fù)雜度的平均熵參數(shù),根據(jù)Shannon 熵,可定義為:

    式中m代表嵌入維數(shù);J代表符號序列的種類;Pj代表每一種符號序列出現(xiàn)的概率。

    當(dāng)Pj=時(shí),PE(m)達(dá)到最大值ln(m!)。采用ln(m!)將PE(m)進(jìn)行歸一化處理,即:

    式中PE代表時(shí)間序列的規(guī)則程度,PE的取值范圍為(0,1]。

    PE值越小,表示時(shí)間序列越規(guī)則穩(wěn)定,PE值越大,說明時(shí)間序列越復(fù)雜隨機(jī)。

    峭度指標(biāo)描述了振動信號的沖擊特征,常被作為診斷滾動軸承故障的評判標(biāo)準(zhǔn),但是在含有強(qiáng)脈沖噪聲或高重復(fù)率瞬變的情況下,峭度指標(biāo)用于故障診斷的有效性會降低。相關(guān)系數(shù)可以表示VMD分解后各子模態(tài)與原始信號的相關(guān)性,但是易受到噪聲干擾。排列熵是一種檢測時(shí)間序列隨機(jī)性和動力學(xué)突變的方法,可以度量系統(tǒng)內(nèi)部絮亂程度,可用于檢測機(jī)械設(shè)備狀態(tài)變化。但僅采用排列熵單一指標(biāo)很難提取到微弱的故障特征信息。因此,構(gòu)造一種融合峭度、相關(guān)系數(shù)和排列熵的綜合指標(biāo)Z,其具體計(jì)算公式為:

    式中Kurtosis為峭度值,r為相關(guān)系數(shù),i為第i個(gè)分量。

    綜合指標(biāo)Z表示與原始信號相關(guān)的分量信號中含有周期性脈沖成分的豐富程度,Z值越大,表示對應(yīng)的復(fù)合故障信息越豐富。

    式中K為經(jīng)VMD 分解后的模態(tài)分量個(gè)數(shù);k為VMD 的分解層數(shù)。

    各指標(biāo)的計(jì)算公式如表1所示,若Zi≥,則認(rèn)為該模態(tài)分量為有效分量,并將有效分量進(jìn)行重構(gòu)。

    表1 各指標(biāo)計(jì)算公式Tab.1 Calculation formula of each index

    指標(biāo)Z不僅考慮了VMD 各分量信號的沖擊特性,保證重構(gòu)信號與原始信號的相關(guān)性,而且能夠檢測分量信號的隨機(jī)性和動力學(xué)突變行為,使得重構(gòu)信號能有效保留復(fù)合故障信息。

    1.2 滾動軸承復(fù)合故障診斷方法和步驟

    VMD 是一種非遞歸的信號分解方法,對于噪聲表現(xiàn)出較好的魯棒性,該方法涉及分解層數(shù)k和懲罰因子α兩個(gè)重要參數(shù),綜合考慮計(jì)算精度和效率[24],設(shè)置k=4 和α=2000。為準(zhǔn)確選取VMD 分解后的IMF 分量,保證重構(gòu)信號既能保留原信號中的故障信息,又能有效過濾噪聲,提高信噪比,采用提出的綜合指標(biāo)Z選取IMF 分量進(jìn)行信號重構(gòu)。

    Autogram 是一種基于無偏自相關(guān)檢測最佳解調(diào)頻帶的新方法,能夠優(yōu)化解調(diào)頻帶中心頻率和帶寬[25]。

    本文提出了一種基于Autogram 的共振解調(diào)和1.5 維譜的滾動軸承復(fù)合故障診斷方法。為了準(zhǔn)確提取復(fù)合故障信號,采用Autogram 方法確定中心頻率和帶寬,選取有效的復(fù)合故障共振信號,對包絡(luò)信號做1.5 維譜分析,有效突出故障特征,提高復(fù)合故障識別的準(zhǔn)確性和直觀性。流程如圖1所示,具體步驟如下:

    圖1 復(fù)合故障診斷流程圖Fig.1 Flowchart of compound fault diagnosis

    (1)采集滾動軸承振動信號原始數(shù)據(jù)。

    (2)對原始信號進(jìn)行VMD 分解,得到相應(yīng)的模態(tài)分量。

    (3)計(jì)算各IMF 分量的綜合指標(biāo)Z及其均值,選取綜合指標(biāo)在均值之上的分量進(jìn)行信號重構(gòu)。

    (4)采用Autogram 確定濾波器參數(shù),即中心頻率和帶寬,并進(jìn)行帶通濾波。

    (5)計(jì)算共振信號的包絡(luò)信號,對包絡(luò)信號做1.5 維譜分析,提取故障特征頻率。

    (6)將1.5 維譜分析得到的特征頻率與理論上的故障特征頻率進(jìn)行對比,識別出復(fù)合故障類型。

    2 復(fù)合故障診斷方法的相關(guān)基本理論

    2.1 變分模態(tài)分解原理

    VMD 是由Dragomiretskiy 等提出的一種自適應(yīng)的準(zhǔn)正交信號分解的方法[26]。VMD 可將多分量信號非遞歸地分解為多個(gè)帶限本征模態(tài)函數(shù)。IMF分量uk(t)是AM-FM 信號,定義為:

    式中Ak(t)為信號包絡(luò),?k(t)表示相位,為非減函數(shù)。瞬時(shí)頻率ωk(t)=?′k(t)是非負(fù)的,且變化慢于?k(t),VMD 由下式表示:

    式中uk為第k個(gè)模態(tài),ωk為中心頻率,δ(t)為脈沖函數(shù),f表示輸入信號。

    公式(6)由二次懲罰因子和拉格朗日乘子決定。增廣拉格朗日方程如下式所示:

    式中α為二次懲罰因子,λ為拉格朗日乘子。

    公式(7)由乘子交替方向法決定。從頻譜中得到的模態(tài)如下:

    VMD 方法主要包括如下步驟:

    1)更新模態(tài)(ω):

    2)更新中心頻率ωn+1k:

    3)當(dāng)所有的ω≥0,更新為:

    當(dāng)滿足下式時(shí):

    停止更新。

    2.2 Autogram 方法

    Autogram 是由Moshrefzadeh 等[25]提出的一種方法,通過計(jì)算解調(diào)信號平方包絡(luò)的無偏自相關(guān)信號的峭度值來選取包含故障信息的信號頻帶。

    無偏自相關(guān)是基于信號平方包絡(luò)進(jìn)行計(jì)算的,計(jì)算公式如下:

    式中X為信號的平方包絡(luò),N為原始信號總長度,q=0,1,…,N?1。

    為了量化頻帶內(nèi)的信號脈沖,Autogram 方法對傳統(tǒng)的峭度方程進(jìn)行了修改,如下式所示:

    通過尋找最大峭度值對應(yīng)分量所處的頻帶,選出最佳解調(diào)頻帶。該方法克服了主觀判斷的缺點(diǎn),所確定的復(fù)合故障共振信號的中心頻率更準(zhǔn)確,帶寬也更合理,提取的故障信號更加準(zhǔn)確。

    2.3 1.5 維譜

    對于零均值的離散平穩(wěn)隨機(jī)過程x(n),其三階累積量C(m1,m2)定義為:

    式中m1,m2為時(shí)間延遲變量,E[?]表示數(shù)學(xué)期望。

    對于公式(15),令m1=m2=m,則三階累積量的對角切片定義為:

    1.5 維譜是三階累積量對角切片的一維傅里葉變換,即對C(m)進(jìn)行一維傅里葉變換可得到x(n)的1.5 維譜,定義為:

    1.5 維譜具有抑制高斯噪聲,有效檢測二次相位耦合特征,強(qiáng)化信號基頻分量等優(yōu)勢,因此,將1.5維譜運(yùn)用于復(fù)合故障的診斷和分析中。

    3 實(shí)例驗(yàn)證與結(jié)果分析

    在Spectrum Quest Incorporated(SQI)公司生產(chǎn)的MFS 試驗(yàn)臺上測試具有復(fù)合故障的滾動軸承,試驗(yàn)臺如圖2所示。三相交流電機(jī)通過柔性聯(lián)軸器與驅(qū)動軸相聯(lián),轉(zhuǎn)軸兩端各安裝一個(gè)型號為ER-16K 深溝軸承,其中一端是復(fù)合故障軸承,另一端是健康軸承,軸承的具體相關(guān)參數(shù)如表2所示。通過在兩個(gè)軸承中間安裝5.1 kg 轉(zhuǎn)子盤施加50 N 的徑向載荷。采用3 通道數(shù)據(jù)線連接信號采集器和加速度傳感器,并通過USB 接口將信號傳輸?shù)接?jì)算機(jī),以此采集振動信號。

    表2 ER-16K 深溝球軸承參數(shù)Tab.2 Parameters of ER?16K deep groove ball bearing

    圖2 MFS 試驗(yàn)臺Fig.2 MFS test rig

    測試用故障軸承的缺陷是采用激光蝕刻技術(shù)加工得到的,缺陷部位如圖3所示,具體參數(shù)如表3所示。

    表3 復(fù)合故障軸承缺陷參數(shù)Tab.3 Defect parameters of compound fault bearing

    圖3 復(fù)合故障軸承缺陷部位Fig.3 Defect location of compound fault bearing

    實(shí)驗(yàn)中設(shè)置不同轉(zhuǎn)速獲取多組數(shù)據(jù),設(shè)置采樣頻率為15600 Hz,采樣點(diǎn)數(shù)為131072 個(gè),隨機(jī)選擇一組數(shù)據(jù)中連續(xù)的60000 個(gè)數(shù)據(jù)點(diǎn)進(jìn)行分析。

    根據(jù)滾動軸承故障特征頻率計(jì)算公式[27]和相關(guān)參數(shù)得到不同轉(zhuǎn)速下復(fù)合故障軸承的故障特征頻率值如表4所示。以3 種不同類型的滾動軸承復(fù)合故障實(shí)驗(yàn)驗(yàn)證所提出方法的可行性。

    表4 不同轉(zhuǎn)速下內(nèi)圈、外圈、滾動體故障特征頻率值Tab.4 Values of fault feature frequency for inner ring,outer ring and rolling element at different speeds

    3.1 內(nèi)圈和外圈復(fù)合故障

    在1200 r/min 時(shí)得到內(nèi)圈和外圈復(fù)合故障(實(shí)驗(yàn)1)振動信號的時(shí)域圖如圖4所示。從時(shí)域圖中可以看到相應(yīng)的沖擊成分,但不能判定是單故障還是復(fù)合故障。采用VMD 方法把原始信號分解成4 個(gè)IMF 分量,計(jì)算各分量綜合指標(biāo)Z,結(jié)果如圖5所示,選取綜合指標(biāo)Z在均值之上的IMF1 和IMF3 分量進(jìn)行信號重構(gòu)。采用Autogram 方法得到中心頻率為6240 Hz,帶寬為960 Hz,如圖6所示。利用所得中心頻率和帶寬對重構(gòu)信號進(jìn)行帶通濾波,獲取共振信號,并進(jìn)行包絡(luò)解調(diào)和1.5 維譜分析,結(jié)果如圖7所示。從圖7中可以明顯地看到外圈故障特征頻率fBPFO及其2 倍頻2fBPFO成分、內(nèi)圈故障特征頻率fBPFI及其2 倍頻2fBPFI成分,表明滾動軸承同時(shí)存在外圈故障和內(nèi)圈故障。結(jié)果表明所提出的方法能夠?qū)崿F(xiàn)準(zhǔn)確直觀地診斷出復(fù)合故障的目的。

    圖4 內(nèi)圈和外圈復(fù)合故障原始信號時(shí)域圖Fig.4 Compound fault of inner ring and outer ring original signal in time domain

    圖5 內(nèi)圈和外圈復(fù)合故障各分量指標(biāo)圖Fig.5 Component index of compound fault of inner ring and outer ring

    圖6 內(nèi)圈和外圈復(fù)合故障AutogramFig.6 Autogram of compound fault of inner ring and outer ring

    圖7 內(nèi)圈和外圈復(fù)合故障1.5 維包絡(luò)譜Fig.7 1.5-dimension envelope spectrum of compound fault of inner ring and outer ring

    為了說明采用綜合指標(biāo)Z選取分量進(jìn)行信號重構(gòu)的優(yōu)點(diǎn),采用相關(guān)系數(shù)單一指標(biāo)選取分量進(jìn)行分析對比。從圖5中得知,IMF2 和IMF3 兩個(gè)分量的相關(guān)系數(shù)值在均值之上,選取這2 個(gè)分量進(jìn)行信號重構(gòu),通過Autogram 方法得到中心頻率為4320 Hz,帶寬為960 Hz。對重構(gòu)信號濾波并進(jìn)行1.5 維包絡(luò)譜分析,所得結(jié)果如圖8所示。圖中僅能看到外圈故障特征頻率,內(nèi)圈故障特征頻率及其2 倍頻則被淹沒在其他干擾譜線中,不能有效準(zhǔn)確地進(jìn)行故障識別,表明所選分量盡管與原始信號相關(guān)性最大,但是包含的故障信息不夠豐富,不能有效地識別出滾動軸承存在的故障。由此可見,綜合指標(biāo)Z能較為準(zhǔn)確地判別含有復(fù)合故障信息的分量,采用綜合指標(biāo)Z選取分量進(jìn)行信號重構(gòu),能有效地對原始信號進(jìn)行降噪。

    圖8 基于相關(guān)系數(shù)的內(nèi)圈和外圈復(fù)合故障1.5 維包絡(luò)譜Fig.8 1.5-dimension envelope spectrum of compound fault of inner ring and outer ring based on correlation coefficient

    3.2 內(nèi)圈和滾動體復(fù)合故障

    為驗(yàn)證所提方法對含有兩種故障信號強(qiáng)弱差異較大復(fù)合故障的識別效果,選取內(nèi)圈和滾動體復(fù)合故障信號進(jìn)行分析。通過試驗(yàn)發(fā)現(xiàn),在內(nèi)圈和滾動體缺陷寬度和深度一樣的情況下,滾動體的故障信號相對微弱和不穩(wěn)定,這可能是滾動軸承在運(yùn)行過程中滾動體在不同位置時(shí)產(chǎn)生的沖擊振動信號差異較大造成的,這給實(shí)際診斷造成了極大的困難。

    為進(jìn)一步擴(kuò)大兩種故障信號強(qiáng)弱的差異性,設(shè)計(jì)試驗(yàn)軸承滾動體缺陷寬度為0.6 mm,小于內(nèi)圈缺陷寬度1.2 mm(如表3所示),從理論上使?jié)L動體產(chǎn)生的故障信號比外圈產(chǎn)生的故障信號微弱。

    在900 r/min 時(shí)得到內(nèi)圈和滾動體復(fù)合故障(實(shí)驗(yàn)2)的振動信號的時(shí)域圖如圖9所示,存在沖擊成分但是不明顯。對該信號進(jìn)行VMD 分解,計(jì)算各分量指標(biāo),得到IMF1 的綜合指標(biāo)為101.8,在綜合指標(biāo)平均值37.0 之上,采用Autogram 方法對IMF1分量計(jì)算中心頻率和帶寬,結(jié)果如圖10所示。

    圖9 內(nèi)圈和滾動體復(fù)合故障原始信號時(shí)域圖Fig.9 Compound fault of inner ring and rolling element original signal in time domain

    圖10 內(nèi)圈和滾動體復(fù)合故障AutogramFig.10 Autogram of compound fault of inner ring and rolling element

    選擇中心頻率3600 Hz 和帶寬480 Hz 進(jìn)行濾波,對共振信號進(jìn)行包絡(luò)處理,并對包絡(luò)信號進(jìn)行1.5 維譜分析,結(jié)果如圖11所示。從圖11中可以明顯看出轉(zhuǎn)頻、滾動體故障特征頻率、內(nèi)圈的故障特征頻率及其倍頻,滾動體2 倍頻不是特別明顯,但仍然能清楚識別,結(jié)果表明所提出的方法可診斷出滾動軸承中同時(shí)存在內(nèi)圈故障和滾動體故障,沒有遺漏故障信號相對微弱的滾動體故障,能識別出同時(shí)含有兩種故障信號強(qiáng)弱差異較大的復(fù)合故障。

    圖11 內(nèi)圈和滾動體復(fù)合故障1.5 維包絡(luò)譜圖Fig.11 1.5-dimension envelope spectrum of compound fault of inner ring and rolling element

    3.3 內(nèi)圈、外圈和滾動體復(fù)合故障

    兩種故障并存是復(fù)合故障中最為簡單的一種形式,由于滾動軸承正常部位與缺陷部位的接觸和相互撞擊,容易產(chǎn)生新的故障從而演化成三種故障,特別是當(dāng)滾動體存在故障時(shí),與之接觸的內(nèi)圈和外圈極易產(chǎn)生新故障,造成三種故障并存的故障狀態(tài),這也是實(shí)際工程中較為常見的一種復(fù)合故障類型。但是目前對滾動軸承復(fù)合故障診斷的研究較多地集中于兩種故障并存的復(fù)合故障類型,相關(guān)的復(fù)合故障診斷方法應(yīng)用到三種故障并存的復(fù)合故障中還相對較少。

    在1200 r/min 時(shí)得到內(nèi)圈、外圈和滾動體三種故障(實(shí)驗(yàn)3)的振動信號的時(shí)域圖如圖12所示,從時(shí)域圖中可以看出沖擊成分明顯復(fù)雜。利用VMD方法分解原始信號并計(jì)算各分量的綜合指標(biāo)Z,IMF3 和IMF4 的分量指標(biāo)分別為11.5 和16.7,在平均綜合指標(biāo)10.9 之上,選擇IMF3 和IMF4 這兩個(gè)分量進(jìn)行信號重構(gòu),利用Autogram 方法得到中心頻率為4440 Hz,帶寬為240 Hz,結(jié)果如圖13所示。對重構(gòu)信號濾波并進(jìn)行包絡(luò)處理,所得包絡(luò)譜如圖14所示,從圖中雖然可以找出內(nèi)圈、外圈和滾動體故障特征頻率,但是仍然存在部分的干擾譜線,外圈故障特征頻率及2 倍頻不夠明顯,被淹沒在其他干擾譜線中,可見,通過包絡(luò)譜圖還難以直觀地進(jìn)行復(fù)合故障識別。

    圖12 內(nèi)圈、外圈和滾動體復(fù)合故障原始信號時(shí)域圖Fig.12 Compound fault of inner ring,outer ring and rolling element original signal in time domain

    圖13 內(nèi)圈、外圈和滾動體復(fù)合故障AutogramFig.13 Autogram of compound fault of inner ring,outer ring and rolling element original signal in time domain

    圖14 內(nèi)圈、外圈和滾動體復(fù)合故障包絡(luò)譜圖Fig.14 Envelope spectrum of compound fault of inner ring,outer ring and rolling element

    為了強(qiáng)化故障頻率特征,對包絡(luò)信號進(jìn)行1.5維譜分析,所得結(jié)果如圖15所示。從圖中可以明顯看到轉(zhuǎn)頻、內(nèi)圈、外圈和滾動體三種故障的故障特征頻率及其2 倍頻,表明1.5 維譜在抑制高斯噪聲,檢測二次相位耦合特征,強(qiáng)化信號基頻分量方面具有明顯優(yōu)勢,有助于復(fù)合故障診斷。

    圖15 內(nèi)圈、外圈和滾動體復(fù)合故障1.5 維包絡(luò)譜圖Fig.15 1.5-dimension envelope spectrum of compound fault of inner ring,outer ring and rolling element

    通過實(shí)驗(yàn)驗(yàn)證,本文所提方法不僅可以較好地診斷出含有內(nèi)圈故障和外圈故障的復(fù)合故障,也能有效識別兩種故障信號強(qiáng)弱差異較大的復(fù)合故障類型,特別是對同時(shí)存在內(nèi)圈、外圈和滾動體三種故障并存的復(fù)合故障也能成功辨識,具有較好的普遍適用性和工程實(shí)際應(yīng)用價(jià)值。

    4 討 論

    為了進(jìn)一步說明所提方法識別滾動軸承復(fù)合故障的可行性和適用性,選取更貼近實(shí)際磨損缺陷的故障試驗(yàn)數(shù)據(jù)——西安交通大學(xué)雷亞國教授課題組的XJTU-SY Bearing Datasets 實(shí)驗(yàn)數(shù)據(jù)進(jìn)一步驗(yàn)證分析,相關(guān)試驗(yàn)條件詳見文獻(xiàn)[28]。本文隨機(jī)選取Bearing 3_2 這組數(shù)據(jù)中的一個(gè)樣本(第2000 個(gè)樣本傳感器在垂直方向上獲取的數(shù)據(jù))進(jìn)行分析。通過計(jì)算得到內(nèi)圈、外圈和滾動體的故障特征頻率分別為197,123 和83 Hz。

    該試驗(yàn)樣本的時(shí)域信號如圖16所示,能看到有沖擊信號,但是不明顯。對該信號進(jìn)行VMD 分解,計(jì)算各分量綜合指標(biāo),得到IMF1 的綜合指標(biāo)為7.7,在綜合平均值3.1 之上,采用Autogram 方法對IMF1 分量計(jì)算,得到中心頻率為960 Hz,帶寬為1920 Hz,如圖17所示。

    圖16 內(nèi)圈、外圈和滾動體復(fù)合故障原始信號時(shí)域圖Fig.16 Compound fault of inner ring,outer ring and rolling element original signal in time domain

    圖17 內(nèi)圈、外圈和滾動體復(fù)合故障AutogramFig.17 Autogram of compound fault of inner ring,outer ring and rolling element original signal in time domain

    對IMF1 進(jìn)行濾波提取共振信號,對共振信號進(jìn)行1.5 維包絡(luò)譜分析,所得結(jié)果如圖18所示,從圖中可以明顯看到轉(zhuǎn)頻、內(nèi)圈、外圈和滾動體三種故障的故障特征頻率及其2 倍頻,雖仍存在部分干擾譜線,但并不影響復(fù)合故障的診斷。

    圖18 內(nèi)圈、外圈和滾動體復(fù)合故障1.5 維包絡(luò)譜圖Fig.18 1.5-dimension envelope spectrum of compound fault of inner ring,outer ring and rolling element

    采用文獻(xiàn)[13]中的診斷方法對該組試驗(yàn)數(shù)據(jù)進(jìn)行診斷,選取相關(guān)系數(shù)最大的分量IMF1 和峭度最大的分量IMF3 進(jìn)行信號重構(gòu)。通過計(jì)算,內(nèi)圈、外圈和滾動體的解卷積周期分別為TBPFI=130,TBPFO=208,TBSF=308,選取濾波器長度分別為LBPFI=300,LBPFO=180,LBSF=80。通過設(shè)置不同參數(shù),采用MCKD 方法分別提取內(nèi)圈、外圈和滾動體故障信號,并對其包絡(luò)信號進(jìn)行1.5 維分析,所得結(jié)果分別如圖19~21 所示。

    圖19 MCKD 方法提取內(nèi)圈故障信號1.5 維包絡(luò)譜圖Fig.19 1.5-dimension envelope spectrum of inner ring fault signal extracted by MCKD

    從圖19中可以明顯看到內(nèi)圈故障特征頻率及2 倍頻,說明有效地提取了內(nèi)圈故障信號。從圖20中可以看到外圈故障特征頻率及2 倍頻,但2 倍頻不是很明顯。圖21中頻譜線較為雜亂無規(guī)律,滾動體的特征頻率被其他頻譜線所干擾,不能有效識別故障,滾動體故障信號提取失敗,說明采用該方法對此組實(shí)驗(yàn)數(shù)據(jù)進(jìn)行故障診斷出現(xiàn)了漏診現(xiàn)象。這可能是由于設(shè)置的參數(shù)與實(shí)際運(yùn)行工況中的參數(shù)存在不一致所導(dǎo)致的。經(jīng)過對比可知,本文所提方法得出的故障特征頻率更加突出明顯,干擾譜線更少,能夠準(zhǔn)確診斷出不同類型故障,體現(xiàn)了所提方法的優(yōu)點(diǎn)。

    圖20 MCKD 方法提取外圈故障信號1.5 維包絡(luò)譜圖Fig.20 1.5-dimension envelope spectrum of outer ring fault signal extracted by MCKD

    圖21 MCKD 方法提取滾動體故障信號1.5 維包絡(luò)譜圖Fig.21 1.5-dimension envelope spectrum of rolling element fault signal extracted by MCKD

    采用基于MCKD 進(jìn)行診斷時(shí),周期參數(shù)是通過采樣頻率和故障特征頻率計(jì)算得到,通過預(yù)設(shè)周期參數(shù)來驗(yàn)證是否存在該故障,周期和濾波器長度等參數(shù)的準(zhǔn)確性又極大地限制了診斷的效果,這在一定程度上存在著局限性,要提高其在工程實(shí)際中的應(yīng)用效果,還需進(jìn)一步改進(jìn)。

    相比而言,本文所提的方法相對簡單和直觀,改變了傳統(tǒng)的先分離后診斷的模式,無須對故障信號分離就可以有效識別復(fù)合故障類型。既可以避免基于分量信號分析的方法由于分量選取不準(zhǔn)確造成誤診,又可以避免基于解卷積方法由于參數(shù)設(shè)置不合理導(dǎo)致信號分離效果不佳造成漏診,有效克服了先分離后診斷的復(fù)合故障診斷方法的局限性,提高了滾動軸承復(fù)合故障診斷的準(zhǔn)確性和有效性。

    5 結(jié) 論

    本文提出了基于Autogram 共振解調(diào)和1.5 維譜的滾動軸承復(fù)合故障診斷方法,采用滾動軸承3 種不同類型復(fù)合故障的實(shí)測信號驗(yàn)證了所提方法的可行性。

    (1)提出了綜合指標(biāo)Z,可以準(zhǔn)確選擇有效的IMF 分量進(jìn)行信號重構(gòu),降低了背景噪聲的干擾。

    (2)采用基于Autogram 共振解調(diào)和1.5 維譜的方法可以有效識別多種復(fù)合故障,減少了漏診和誤診的可能性。

    該方法克服了先分離信號后診斷方法的局限性,具有一定的工程應(yīng)用價(jià)值。

    猜你喜歡
    峭度特征頻率內(nèi)圈
    基于MCKD和峭度的液壓泵故障特征提取
    特種復(fù)合軸承內(nèi)圈推力滾道磨削用工裝設(shè)計(jì)
    哈爾濱軸承(2021年4期)2021-03-08 01:00:48
    聯(lián)合快速峭度圖與變帶寬包絡(luò)譜峭度圖的輪對軸承復(fù)合故障檢測研究
    瓷磚檢測機(jī)器人的聲音信號處理
    主軸軸承內(nèi)圈鎖緊用臺階套的裝配
    光學(xué)波前參數(shù)的分析評價(jià)方法研究
    基于振動信號特征頻率的數(shù)控車床故障辨識方法
    基于峭度分析的聲發(fā)射故障檢測
    電子世界(2018年12期)2018-07-04 06:34:38
    基于小波去噪和EMD算法在齒輪故障檢測中的應(yīng)用
    基于魯棒性小波包峭度圖的滾動軸承故障診斷*
    国产精华一区二区三区| 欧美av亚洲av综合av国产av| 黄色怎么调成土黄色| 在线观看日韩欧美| 成人免费观看视频高清| 天天躁狠狠躁夜夜躁狠狠躁| 欧美成人午夜精品| 黑人巨大精品欧美一区二区mp4| 伦理电影免费视频| av欧美777| 999久久久国产精品视频| 999久久久国产精品视频| 亚洲美女黄片视频| 免费少妇av软件| 80岁老熟妇乱子伦牲交| 久久久国产欧美日韩av| 久久国产乱子伦精品免费另类| 国产高清视频在线播放一区| 亚洲欧美日韩另类电影网站| 俄罗斯特黄特色一大片| 亚洲国产精品sss在线观看 | 如日韩欧美国产精品一区二区三区| 纯流量卡能插随身wifi吗| 十分钟在线观看高清视频www| 90打野战视频偷拍视频| 久久精品影院6| 久久精品国产清高在天天线| 午夜免费成人在线视频| 97碰自拍视频| 一进一出抽搐gif免费好疼 | 久久婷婷成人综合色麻豆| 麻豆一二三区av精品| 久久久国产成人精品二区 | 亚洲 欧美一区二区三区| 久久久久久久精品吃奶| 88av欧美| 在线永久观看黄色视频| 黄色女人牲交| 精品电影一区二区在线| av免费在线观看网站| 成人永久免费在线观看视频| 久久久久久久久中文| 91精品国产国语对白视频| 丝袜美腿诱惑在线| 国产一区二区三区视频了| 熟女少妇亚洲综合色aaa.| 中国美女看黄片| 在线观看一区二区三区| 悠悠久久av| 美女福利国产在线| 黄色a级毛片大全视频| 欧美黄色淫秽网站| 99国产精品一区二区蜜桃av| 51午夜福利影视在线观看| 国产精品自产拍在线观看55亚洲| 99在线人妻在线中文字幕| 天堂√8在线中文| 88av欧美| 天天影视国产精品| 国产精品一区二区免费欧美| 麻豆av在线久日| 久久久国产欧美日韩av| 最近最新中文字幕大全电影3 | 欧美精品一区二区免费开放| 俄罗斯特黄特色一大片| 欧美黑人欧美精品刺激| 国产一区在线观看成人免费| 欧美日韩黄片免| 高清欧美精品videossex| 亚洲精品国产色婷婷电影| 一区在线观看完整版| 亚洲精品国产精品久久久不卡| 大型黄色视频在线免费观看| 黄色毛片三级朝国网站| 久久久久精品国产欧美久久久| 欧美一级毛片孕妇| 国产精品亚洲av一区麻豆| 欧美黄色淫秽网站| 男女高潮啪啪啪动态图| 制服诱惑二区| 精品一品国产午夜福利视频| 亚洲五月色婷婷综合| 成人国语在线视频| 亚洲免费av在线视频| 高清黄色对白视频在线免费看| 女生性感内裤真人,穿戴方法视频| 国产精品二区激情视频| 欧美黑人欧美精品刺激| 天堂√8在线中文| 99热只有精品国产| 美女高潮喷水抽搐中文字幕| 亚洲久久久国产精品| 国产一区二区三区综合在线观看| 女生性感内裤真人,穿戴方法视频| 久久久久国产精品人妻aⅴ院| 欧美一级毛片孕妇| 免费女性裸体啪啪无遮挡网站| 亚洲第一av免费看| 日韩精品中文字幕看吧| 国产精品 欧美亚洲| 欧美丝袜亚洲另类 | 人成视频在线观看免费观看| 欧美激情高清一区二区三区| av福利片在线| 欧美日本亚洲视频在线播放| 欧美国产精品va在线观看不卡| 欧美丝袜亚洲另类 | 久久久久国内视频| 91av网站免费观看| 亚洲成av片中文字幕在线观看| 日韩成人在线观看一区二区三区| 欧美日韩亚洲高清精品| 亚洲专区中文字幕在线| 男人舔女人的私密视频| 国产av在哪里看| 男女下面进入的视频免费午夜 | 亚洲第一青青草原| 日韩欧美三级三区| 91成年电影在线观看| 日本a在线网址| 99精品欧美一区二区三区四区| 国产主播在线观看一区二区| 日韩欧美三级三区| 嫩草影院精品99| 51午夜福利影视在线观看| 欧美激情极品国产一区二区三区| 国产精品国产av在线观看| 18禁国产床啪视频网站| 美女高潮喷水抽搐中文字幕| 韩国精品一区二区三区| 婷婷精品国产亚洲av在线| 看黄色毛片网站| 国产激情久久老熟女| aaaaa片日本免费| 18禁黄网站禁片午夜丰满| 国产熟女午夜一区二区三区| 日本免费a在线| 亚洲精品国产精品久久久不卡| 国产午夜精品久久久久久| 国产午夜精品久久久久久| 国产一区二区三区综合在线观看| 在线永久观看黄色视频| 成人手机av| 99国产综合亚洲精品| 精品久久久久久久毛片微露脸| 亚洲中文av在线| 国产午夜精品久久久久久| 男女下面进入的视频免费午夜 | 欧美老熟妇乱子伦牲交| 久久亚洲精品不卡| 岛国视频午夜一区免费看| 国产一区二区在线av高清观看| 露出奶头的视频| 1024视频免费在线观看| 精品一区二区三区四区五区乱码| 国产一区二区三区综合在线观看| 一区福利在线观看| 麻豆成人av在线观看| 一级毛片女人18水好多| 久久久精品国产亚洲av高清涩受| 国产成人精品无人区| 久久99一区二区三区| 国产一区二区三区在线臀色熟女 | 正在播放国产对白刺激| 精品国产乱子伦一区二区三区| 热re99久久国产66热| 欧美性长视频在线观看| 黑人巨大精品欧美一区二区mp4| 欧美日韩亚洲高清精品| 三上悠亚av全集在线观看| 成人影院久久| 99精品久久久久人妻精品| 精品福利永久在线观看| 天堂动漫精品| 他把我摸到了高潮在线观看| 日本 av在线| 久久精品国产亚洲av香蕉五月| av视频免费观看在线观看| 91麻豆精品激情在线观看国产 | 69精品国产乱码久久久| 青草久久国产| 日韩人妻精品一区2区三区| 大型av网站在线播放| 啦啦啦 在线观看视频| 欧美激情 高清一区二区三区| 黄片小视频在线播放| 久久久久久大精品| 久久久久久久精品吃奶| 别揉我奶头~嗯~啊~动态视频| 国产97色在线日韩免费| 国产日韩一区二区三区精品不卡| 9色porny在线观看| 熟女少妇亚洲综合色aaa.| 亚洲免费av在线视频| 精品无人区乱码1区二区| 免费在线观看亚洲国产| 老司机福利观看| 国产精品秋霞免费鲁丝片| 国产极品粉嫩免费观看在线| 日韩欧美在线二视频| 一级毛片精品| 侵犯人妻中文字幕一二三四区| 国产亚洲精品久久久久久毛片| 男人舔女人的私密视频| av欧美777| 国产亚洲欧美在线一区二区| 69av精品久久久久久| 在线观看免费午夜福利视频| 国产午夜精品久久久久久| 国产野战对白在线观看| 99国产极品粉嫩在线观看| 亚洲色图av天堂| 黑人猛操日本美女一级片| 在线观看免费日韩欧美大片| 亚洲人成电影免费在线| 1024视频免费在线观看| 久久九九热精品免费| 欧美黄色淫秽网站| 国产人伦9x9x在线观看| 我的亚洲天堂| 久久国产亚洲av麻豆专区| 午夜激情av网站| 午夜两性在线视频| 亚洲av五月六月丁香网| 成年版毛片免费区| 手机成人av网站| 久久久国产一区二区| 丝袜人妻中文字幕| 男人的好看免费观看在线视频 | 最近最新中文字幕大全免费视频| 中国美女看黄片| 国产成人欧美在线观看| 黑人猛操日本美女一级片| 高清在线国产一区| 中文字幕精品免费在线观看视频| 真人做人爱边吃奶动态| 欧美另类亚洲清纯唯美| 伦理电影免费视频| 99精品在免费线老司机午夜| 亚洲国产欧美网| 美国免费a级毛片| 亚洲精品中文字幕在线视频| 日本五十路高清| 人人妻人人爽人人添夜夜欢视频| 亚洲精华国产精华精| 亚洲av熟女| 免费av毛片视频| 国产三级黄色录像| 免费搜索国产男女视频| 国产精品美女特级片免费视频播放器 | 韩国精品一区二区三区| 久久久久久人人人人人| 国产男靠女视频免费网站| 欧美黄色片欧美黄色片| 久久精品人人爽人人爽视色| 日韩人妻精品一区2区三区| 国产精品秋霞免费鲁丝片| 欧洲精品卡2卡3卡4卡5卡区| 亚洲午夜理论影院| 亚洲精品中文字幕在线视频| 丰满的人妻完整版| 交换朋友夫妻互换小说| 中文字幕人妻丝袜制服| 亚洲七黄色美女视频| 黄频高清免费视频| 天堂√8在线中文| 日韩视频一区二区在线观看| 亚洲国产欧美一区二区综合| 国产亚洲欧美98| 成人影院久久| 女人爽到高潮嗷嗷叫在线视频| 亚洲五月色婷婷综合| 一进一出抽搐gif免费好疼 | 日韩有码中文字幕| 亚洲欧洲精品一区二区精品久久久| 搡老熟女国产l中国老女人| 亚洲欧美精品综合一区二区三区| ponron亚洲| 亚洲精品一二三| 免费女性裸体啪啪无遮挡网站| 国产国语露脸激情在线看| 亚洲精品国产色婷婷电影| 欧美人与性动交α欧美软件| 精品久久久久久久久久免费视频 | 久久热在线av| netflix在线观看网站| 久久人人爽av亚洲精品天堂| 国产日韩一区二区三区精品不卡| 国产成人精品久久二区二区免费| 久久人人爽av亚洲精品天堂| 国产亚洲精品一区二区www| 久久久久久久久中文| 国产免费男女视频| 成人特级黄色片久久久久久久| 人人妻,人人澡人人爽秒播| 亚洲久久久国产精品| www.www免费av| 国产97色在线日韩免费| 日日爽夜夜爽网站| 女同久久另类99精品国产91| 美女大奶头视频| 婷婷六月久久综合丁香| 亚洲精品美女久久久久99蜜臀| av有码第一页| 99国产综合亚洲精品| 成年人黄色毛片网站| 法律面前人人平等表现在哪些方面| 香蕉久久夜色| 99国产精品一区二区三区| 一级黄色大片毛片| 午夜两性在线视频| 丰满饥渴人妻一区二区三| 韩国精品一区二区三区| 精品国产美女av久久久久小说| 国产精品九九99| 国产精品久久视频播放| 亚洲一区高清亚洲精品| 9色porny在线观看| 国产亚洲精品一区二区www| 国产精品免费视频内射| 亚洲色图av天堂| 18禁观看日本| 在线观看免费视频网站a站| 日韩欧美免费精品| 99热国产这里只有精品6| 色婷婷av一区二区三区视频| av网站在线播放免费| 啦啦啦在线免费观看视频4| 色老头精品视频在线观看| 精品久久久久久,| 黑丝袜美女国产一区| 在线观看免费午夜福利视频| 亚洲成a人片在线一区二区| www.熟女人妻精品国产| av电影中文网址| 亚洲av成人不卡在线观看播放网| 99国产精品免费福利视频| 又大又爽又粗| 一区二区三区精品91| 免费不卡黄色视频| 久久国产亚洲av麻豆专区| 亚洲成a人片在线一区二区| 久久久国产一区二区| 熟女少妇亚洲综合色aaa.| 亚洲精品av麻豆狂野| 亚洲中文字幕日韩| 国产成人系列免费观看| 免费观看人在逋| 国产精品久久视频播放| 青草久久国产| 亚洲男人的天堂狠狠| 国产人伦9x9x在线观看| 在线av久久热| 91麻豆av在线| 级片在线观看| 亚洲国产欧美网| 一区二区三区激情视频| 国产精品久久久人人做人人爽| 最近最新免费中文字幕在线| 亚洲精品国产一区二区精华液| 精品国产国语对白av| 757午夜福利合集在线观看| 在线观看www视频免费| 99久久久亚洲精品蜜臀av| 午夜福利欧美成人| 亚洲一码二码三码区别大吗| 亚洲成国产人片在线观看| 色综合欧美亚洲国产小说| 亚洲精品av麻豆狂野| 999精品在线视频| 国产精品二区激情视频| 真人一进一出gif抽搐免费| 在线观看免费高清a一片| 国产99白浆流出| 免费av毛片视频| 制服诱惑二区| 欧美大码av| 日本黄色视频三级网站网址| 美女高潮喷水抽搐中文字幕| 成人手机av| 在线观看日韩欧美| 国产精品电影一区二区三区| 精品福利观看| 亚洲精品久久成人aⅴ小说| 99在线人妻在线中文字幕| 真人一进一出gif抽搐免费| 熟女少妇亚洲综合色aaa.| 欧美性长视频在线观看| 免费在线观看影片大全网站| 久久天堂一区二区三区四区| 天堂√8在线中文| 日韩人妻精品一区2区三区| 成人18禁在线播放| 一边摸一边做爽爽视频免费| 久久国产精品男人的天堂亚洲| 精品免费久久久久久久清纯| 欧美丝袜亚洲另类 | 一进一出抽搐动态| 一级片'在线观看视频| 夜夜躁狠狠躁天天躁| 人妻丰满熟妇av一区二区三区| 亚洲欧美一区二区三区久久| 久久香蕉精品热| 一进一出好大好爽视频| 久久天躁狠狠躁夜夜2o2o| 日本三级黄在线观看| 亚洲欧美日韩另类电影网站| 精品国产一区二区三区四区第35| 最近最新中文字幕大全免费视频| 精品欧美一区二区三区在线| 欧美国产精品va在线观看不卡| 国产又爽黄色视频| 国产一区二区三区在线臀色熟女 | 可以在线观看毛片的网站| 黑人巨大精品欧美一区二区mp4| 国产亚洲精品第一综合不卡| 久久久久亚洲av毛片大全| 啦啦啦免费观看视频1| 国产一卡二卡三卡精品| 亚洲国产欧美一区二区综合| 成熟少妇高潮喷水视频| 老司机亚洲免费影院| 一进一出好大好爽视频| 高清黄色对白视频在线免费看| 伦理电影免费视频| 999久久久国产精品视频| 欧美日韩乱码在线| 欧美另类亚洲清纯唯美| 久久精品亚洲精品国产色婷小说| 欧美成人午夜精品| 亚洲自拍偷在线| 午夜福利免费观看在线| 欧美日韩一级在线毛片| 午夜日韩欧美国产| 男人舔女人下体高潮全视频| 香蕉丝袜av| 中文字幕人妻熟女乱码| 男男h啪啪无遮挡| 巨乳人妻的诱惑在线观看| 老鸭窝网址在线观看| 精品国内亚洲2022精品成人| 国产无遮挡羞羞视频在线观看| 欧美老熟妇乱子伦牲交| 亚洲国产精品999在线| 精品乱码久久久久久99久播| av超薄肉色丝袜交足视频| 国产亚洲欧美精品永久| 高清欧美精品videossex| 男女午夜视频在线观看| 亚洲狠狠婷婷综合久久图片| 新久久久久国产一级毛片| 亚洲精华国产精华精| 精品第一国产精品| 日韩欧美国产一区二区入口| 69av精品久久久久久| 亚洲在线自拍视频| 亚洲国产精品999在线| 男人的好看免费观看在线视频 | 久久伊人香网站| 久久久久久人人人人人| 日韩av在线大香蕉| 国产99久久九九免费精品| 欧美av亚洲av综合av国产av| 99riav亚洲国产免费| 中文字幕高清在线视频| avwww免费| 99精品久久久久人妻精品| 手机成人av网站| 国产黄a三级三级三级人| 日本撒尿小便嘘嘘汇集6| 99国产精品一区二区三区| 午夜免费观看网址| 国产国语露脸激情在线看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲成人国产一区在线观看| 三级毛片av免费| 欧美日韩av久久| 19禁男女啪啪无遮挡网站| 精品久久久久久成人av| 老司机靠b影院| 日本免费a在线| 国产人伦9x9x在线观看| 99国产极品粉嫩在线观看| 大香蕉久久成人网| av电影中文网址| 一二三四社区在线视频社区8| 夜夜夜夜夜久久久久| 少妇 在线观看| 国产一区二区激情短视频| 国产亚洲欧美精品永久| 午夜精品国产一区二区电影| av超薄肉色丝袜交足视频| 精品卡一卡二卡四卡免费| 欧美老熟妇乱子伦牲交| 亚洲第一欧美日韩一区二区三区| 亚洲精品中文字幕在线视频| 亚洲一码二码三码区别大吗| 麻豆一二三区av精品| 亚洲精品一二三| 人妻久久中文字幕网| 亚洲国产欧美日韩在线播放| 亚洲aⅴ乱码一区二区在线播放 | 免费av中文字幕在线| 看黄色毛片网站| 岛国视频午夜一区免费看| www.精华液| 日韩成人在线观看一区二区三区| 99精品久久久久人妻精品| 妹子高潮喷水视频| a在线观看视频网站| 嫁个100分男人电影在线观看| 9色porny在线观看| 村上凉子中文字幕在线| 久久香蕉激情| 国产精品99久久99久久久不卡| 国产精品爽爽va在线观看网站 | 伦理电影免费视频| 亚洲国产毛片av蜜桃av| av福利片在线| 侵犯人妻中文字幕一二三四区| 亚洲九九香蕉| 黄色毛片三级朝国网站| 18禁国产床啪视频网站| 久久久国产一区二区| 国产深夜福利视频在线观看| 国产精品成人在线| 精品福利观看| 51午夜福利影视在线观看| 9191精品国产免费久久| 亚洲av成人av| 国产高清国产精品国产三级| 丁香欧美五月| 国产极品粉嫩免费观看在线| 欧美亚洲日本最大视频资源| 久久中文字幕一级| 国产精品一区二区免费欧美| 亚洲精品久久成人aⅴ小说| 午夜精品久久久久久毛片777| 欧美国产精品va在线观看不卡| 中文字幕人妻丝袜一区二区| 亚洲成国产人片在线观看| 三上悠亚av全集在线观看| 欧美在线一区亚洲| 亚洲 欧美一区二区三区| avwww免费| 日本wwww免费看| 身体一侧抽搐| 国产色视频综合| 亚洲一区中文字幕在线| 免费高清在线观看日韩| 精品国内亚洲2022精品成人| 精品国产美女av久久久久小说| 精品日产1卡2卡| 99国产极品粉嫩在线观看| 国产欧美日韩综合在线一区二区| 在线观看www视频免费| 一边摸一边抽搐一进一小说| 亚洲一区二区三区不卡视频| 国产一区二区激情短视频| 免费在线观看日本一区| 亚洲,欧美精品.| 国产av一区在线观看免费| 精品国产一区二区三区四区第35| 十分钟在线观看高清视频www| 欧美+亚洲+日韩+国产| 69精品国产乱码久久久| 夜夜看夜夜爽夜夜摸 | 一级作爱视频免费观看| 亚洲成人免费电影在线观看| 亚洲国产欧美一区二区综合| 亚洲av成人一区二区三| 天堂√8在线中文| 欧美久久黑人一区二区| 免费av毛片视频| √禁漫天堂资源中文www| 免费少妇av软件| 欧美黄色淫秽网站| av天堂在线播放| 久久久精品国产亚洲av高清涩受| 午夜视频精品福利| 欧美人与性动交α欧美精品济南到| 纯流量卡能插随身wifi吗| 久久久久九九精品影院| 黄色 视频免费看| 99精品久久久久人妻精品| 夜夜看夜夜爽夜夜摸 | 可以免费在线观看a视频的电影网站| 99国产极品粉嫩在线观看| 久久久精品国产亚洲av高清涩受| 午夜视频精品福利| 成人手机av| 妹子高潮喷水视频| 成人手机av| 亚洲中文字幕日韩| 十八禁人妻一区二区| 一区二区三区激情视频| 两性夫妻黄色片| 国产一区二区三区视频了| 久久九九热精品免费| 黄网站色视频无遮挡免费观看| 在线观看午夜福利视频| 麻豆成人av在线观看| 中国美女看黄片| 久热爱精品视频在线9| 这个男人来自地球电影免费观看| 侵犯人妻中文字幕一二三四区| 久久亚洲真实| 欧美色视频一区免费| 久热爱精品视频在线9| av超薄肉色丝袜交足视频| a级毛片黄视频| 国产精品久久久久久人妻精品电影| 欧美日韩瑟瑟在线播放| 少妇裸体淫交视频免费看高清 | 成在线人永久免费视频|