• <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)用
    基于魯棒性小波包峭度圖的滾動軸承故障診斷*
    午夜两性在线视频| 亚洲avbb在线观看| 亚洲一区中文字幕在线| 日本欧美视频一区| 久久久久久久精品吃奶| 久久久国产成人精品二区| 久久热在线av| 在线观看66精品国产| 免费少妇av软件| 极品教师在线免费播放| 久久精品国产99精品国产亚洲性色 | av天堂久久9| 成人国产综合亚洲| 久热这里只有精品99| 黑丝袜美女国产一区| 午夜精品久久久久久毛片777| 丰满人妻熟妇乱又伦精品不卡| 757午夜福利合集在线观看| 成人国产一区最新在线观看| 变态另类丝袜制服| 可以在线观看的亚洲视频| 午夜免费鲁丝| 黄色 视频免费看| 一进一出好大好爽视频| 黄片大片在线免费观看| 国内精品久久久久精免费| 在线观看免费视频日本深夜| 中文字幕人妻丝袜一区二区| 免费在线观看黄色视频的| 一级毛片女人18水好多| 国产精品综合久久久久久久免费 | 首页视频小说图片口味搜索| 国产国语露脸激情在线看| av片东京热男人的天堂| 精品久久久精品久久久| 亚洲国产精品sss在线观看| 久久久精品欧美日韩精品| 国产欧美日韩综合在线一区二区| 亚洲午夜理论影院| 美女大奶头视频| 中文字幕另类日韩欧美亚洲嫩草| 老司机福利观看| 久久久久久久久中文| 精品久久久久久久久久免费视频| 欧美乱妇无乱码| 中文字幕色久视频| 精品福利观看| 免费在线观看完整版高清| 又大又爽又粗| 50天的宝宝边吃奶边哭怎么回事| 此物有八面人人有两片| 欧美日韩乱码在线| 国产97色在线日韩免费| 91在线观看av| 亚洲第一av免费看| 18美女黄网站色大片免费观看| www.精华液| 激情在线观看视频在线高清| 亚洲第一av免费看| 十八禁人妻一区二区| 级片在线观看| 搡老熟女国产l中国老女人| 黄色 视频免费看| 国产黄a三级三级三级人| 午夜精品久久久久久毛片777| 欧美最黄视频在线播放免费| 欧美+亚洲+日韩+国产| 一个人免费在线观看的高清视频| 久久精品国产99精品国产亚洲性色 | 黄色 视频免费看| 天堂√8在线中文| 久9热在线精品视频| 免费在线观看日本一区| 一进一出抽搐gif免费好疼| АⅤ资源中文在线天堂| 一级作爱视频免费观看| 国产亚洲精品av在线| 18禁国产床啪视频网站| av网站免费在线观看视频| www.精华液| 黑人巨大精品欧美一区二区mp4| 1024香蕉在线观看| 操出白浆在线播放| 少妇熟女aⅴ在线视频| 啦啦啦免费观看视频1| 欧美在线一区亚洲| 黄色a级毛片大全视频| 欧美亚洲日本最大视频资源| 丝袜美腿诱惑在线| 久久伊人香网站| 国产亚洲av高清不卡| 国产成年人精品一区二区| 好男人在线观看高清免费视频 | 国产欧美日韩一区二区三区在线| 精品国产乱码久久久久久男人| 亚洲第一欧美日韩一区二区三区| 国产精品亚洲一级av第二区| 亚洲av日韩精品久久久久久密| 美女 人体艺术 gogo| 国产亚洲精品综合一区在线观看 | 欧美一区二区精品小视频在线| 国产精品日韩av在线免费观看 | 亚洲视频免费观看视频| 禁无遮挡网站| 国产亚洲精品av在线| 免费少妇av软件| 母亲3免费完整高清在线观看| √禁漫天堂资源中文www| 变态另类成人亚洲欧美熟女 | 在线观看舔阴道视频| 亚洲国产精品sss在线观看| 男女下面插进去视频免费观看| 久久伊人香网站| 精品一区二区三区视频在线观看免费| 免费高清在线观看日韩| 欧洲精品卡2卡3卡4卡5卡区| 国产99白浆流出| 90打野战视频偷拍视频| 亚洲欧美精品综合久久99| 777久久人妻少妇嫩草av网站| 国产一卡二卡三卡精品| 久久久精品国产亚洲av高清涩受| 成人国产综合亚洲| 变态另类丝袜制服| 亚洲av电影不卡..在线观看| 久久久久久人人人人人| 久久人妻熟女aⅴ| 午夜成年电影在线免费观看| 亚洲人成网站在线播放欧美日韩| 亚洲无线在线观看| 日韩欧美在线二视频| 国产视频一区二区在线看| 长腿黑丝高跟| 桃色一区二区三区在线观看| 亚洲人成网站在线播放欧美日韩| 黑丝袜美女国产一区| 人妻久久中文字幕网| 国产亚洲av高清不卡| 国产黄a三级三级三级人| 两人在一起打扑克的视频| 久久久久国产精品人妻aⅴ院| 久久精品国产亚洲av香蕉五月| 黄片大片在线免费观看| 又黄又粗又硬又大视频| 国产精品一区二区三区四区久久 | 狂野欧美激情性xxxx| 亚洲国产精品合色在线| 天堂√8在线中文| 两人在一起打扑克的视频| 曰老女人黄片| 国产一区二区在线av高清观看| 很黄的视频免费| 亚洲av五月六月丁香网| 久久久国产精品麻豆| 老司机午夜福利在线观看视频| 久久国产精品影院| 黑丝袜美女国产一区| 欧美中文综合在线视频| 好男人电影高清在线观看| 日日干狠狠操夜夜爽| 国产精品一区二区在线不卡| 久久国产乱子伦精品免费另类| 国产精品久久视频播放| 99久久99久久久精品蜜桃| 国产xxxxx性猛交| 成人手机av| 一级毛片高清免费大全| 少妇被粗大的猛进出69影院| 亚洲人成伊人成综合网2020| 亚洲欧美日韩另类电影网站| 亚洲专区中文字幕在线| 国产成人精品在线电影| 国产aⅴ精品一区二区三区波| www.自偷自拍.com| 啪啪无遮挡十八禁网站| 在线国产一区二区在线| 久久中文看片网| 精品无人区乱码1区二区| 成人国产综合亚洲| 亚洲人成77777在线视频| 精品久久久久久成人av| 波多野结衣高清无吗| 久久久久久久久久久久大奶| 色播在线永久视频| 美女高潮喷水抽搐中文字幕| 国产亚洲欧美98| 在线观看66精品国产| 色综合站精品国产| 久久国产精品人妻蜜桃| 国产真人三级小视频在线观看| 午夜免费激情av| 精品国产国语对白av| 在线观看免费视频网站a站| 一区二区三区精品91| 亚洲国产欧美网| 国产精品久久久人人做人人爽| 亚洲人成电影免费在线| 国产亚洲av高清不卡| 亚洲狠狠婷婷综合久久图片| 欧美一级a爱片免费观看看 | 免费看a级黄色片| 中国美女看黄片| 久久久久久久久中文| 久久婷婷人人爽人人干人人爱 | 一区福利在线观看| 两人在一起打扑克的视频| 国产视频一区二区在线看| 成在线人永久免费视频| 欧美一级a爱片免费观看看 | 美女午夜性视频免费| tocl精华| 久久精品亚洲熟妇少妇任你| 国产精品爽爽va在线观看网站 | 看黄色毛片网站| 一区福利在线观看| 久久久国产成人精品二区| 久久久久久免费高清国产稀缺| 少妇裸体淫交视频免费看高清 | 男人操女人黄网站| 欧美激情久久久久久爽电影 | 欧美 亚洲 国产 日韩一| 亚洲第一青青草原| 久久草成人影院| 免费在线观看影片大全网站| 这个男人来自地球电影免费观看| 欧美 亚洲 国产 日韩一| www.999成人在线观看| 可以免费在线观看a视频的电影网站| 精品一区二区三区视频在线观看免费| 日本 欧美在线| 色综合婷婷激情| 久久天堂一区二区三区四区| 亚洲av电影不卡..在线观看| 午夜福利18| 91国产中文字幕| 成人永久免费在线观看视频| 成人亚洲精品一区在线观看| 看免费av毛片| 亚洲第一青青草原| 国产精品永久免费网站| 亚洲少妇的诱惑av| 国产精品美女特级片免费视频播放器 | 久久久久久大精品| 午夜免费鲁丝| 国产精品久久久久久精品电影 | 国产精品乱码一区二三区的特点 | 亚洲国产精品合色在线| 日韩欧美免费精品| 久久久久久久久中文| 黄片小视频在线播放| 精品欧美一区二区三区在线| 人人妻人人澡欧美一区二区 | 波多野结衣av一区二区av| 国产熟女午夜一区二区三区| 亚洲电影在线观看av| 黄频高清免费视频| 午夜日韩欧美国产| 19禁男女啪啪无遮挡网站| 日本三级黄在线观看| 人人妻,人人澡人人爽秒播| 一区二区三区精品91| 亚洲国产精品合色在线| 久久人妻熟女aⅴ| 看片在线看免费视频| 久久中文看片网| 欧美日韩亚洲国产一区二区在线观看| av中文乱码字幕在线| 欧美乱妇无乱码| 成年版毛片免费区| 久久精品国产亚洲av高清一级| 黄色毛片三级朝国网站| 欧美午夜高清在线| 日韩欧美在线二视频| 亚洲美女黄片视频| 亚洲成a人片在线一区二区| 悠悠久久av| 免费看十八禁软件| 黄色成人免费大全| 老鸭窝网址在线观看| 国产亚洲精品av在线| 天堂√8在线中文| 亚洲欧洲精品一区二区精品久久久| 国产野战对白在线观看| 亚洲人成伊人成综合网2020| 自线自在国产av| 成年人黄色毛片网站| 国语自产精品视频在线第100页| 久久中文字幕一级| 精品卡一卡二卡四卡免费| 日本精品一区二区三区蜜桃| 亚洲精品国产色婷婷电影| 免费少妇av软件| 日本撒尿小便嘘嘘汇集6| 999精品在线视频| 又大又爽又粗| 在线观看免费视频日本深夜| 亚洲黑人精品在线| 9热在线视频观看99| 丝袜美腿诱惑在线| 香蕉久久夜色| 黄色视频,在线免费观看| 免费在线观看日本一区| 欧美老熟妇乱子伦牲交| 国产片内射在线| 国产一区在线观看成人免费| 精品欧美一区二区三区在线| 久久国产精品男人的天堂亚洲| 美女 人体艺术 gogo| 亚洲国产欧美日韩在线播放| 国产欧美日韩一区二区精品| 每晚都被弄得嗷嗷叫到高潮| 欧美色视频一区免费| 国产亚洲av高清不卡| 一个人观看的视频www高清免费观看 | 日本欧美视频一区| 欧美成人免费av一区二区三区| 此物有八面人人有两片| 国内精品久久久久久久电影| 久久久久亚洲av毛片大全| 国产精品美女特级片免费视频播放器 | 香蕉丝袜av| 久久国产亚洲av麻豆专区| 亚洲久久久国产精品| 欧美另类亚洲清纯唯美| 一a级毛片在线观看| 久久精品国产99精品国产亚洲性色 | 一个人免费在线观看的高清视频| 午夜日韩欧美国产| 国产av精品麻豆| 午夜日韩欧美国产| av在线天堂中文字幕| 男男h啪啪无遮挡| 女生性感内裤真人,穿戴方法视频| 精品久久久久久久久久免费视频| 十八禁网站免费在线| 9色porny在线观看| 美女高潮喷水抽搐中文字幕| 啦啦啦观看免费观看视频高清 | 亚洲一区二区三区不卡视频| 精品国产乱码久久久久久男人| 99久久99久久久精品蜜桃| 国产精品日韩av在线免费观看 | 亚洲欧美一区二区三区黑人| 少妇熟女aⅴ在线视频| 国产亚洲精品久久久久5区| 国产熟女xx| 免费在线观看完整版高清| 很黄的视频免费| 女人被躁到高潮嗷嗷叫费观| 高清毛片免费观看视频网站| 欧美绝顶高潮抽搐喷水| av有码第一页| 国产亚洲精品综合一区在线观看 | 国产欧美日韩精品亚洲av| 精品久久久久久久久久免费视频| 欧美午夜高清在线| 精品第一国产精品| 亚洲精品美女久久久久99蜜臀| 国产麻豆成人av免费视频| 9191精品国产免费久久| 日韩欧美免费精品| 女性被躁到高潮视频| 亚洲精品中文字幕一二三四区| 中文字幕精品免费在线观看视频| 满18在线观看网站| 亚洲午夜精品一区,二区,三区| 搡老妇女老女人老熟妇| 欧美不卡视频在线免费观看 | 制服丝袜大香蕉在线| 国产伦一二天堂av在线观看| 欧美成人性av电影在线观看| 99国产精品一区二区三区| 男女之事视频高清在线观看| 欧美不卡视频在线免费观看 | 精品免费久久久久久久清纯| 天天一区二区日本电影三级 | 亚洲精品国产区一区二| 欧美乱妇无乱码| 亚洲一区二区三区色噜噜| 亚洲成av人片免费观看| 久久婷婷成人综合色麻豆| 波多野结衣一区麻豆| 悠悠久久av| 亚洲人成电影免费在线| av天堂在线播放| 97超级碰碰碰精品色视频在线观看| 午夜免费成人在线视频| 精品欧美一区二区三区在线| 国产av一区二区精品久久| 国内久久婷婷六月综合欲色啪| 日韩欧美一区二区三区在线观看| 久久久久久久久免费视频了| 国产黄a三级三级三级人| 看片在线看免费视频| 久久久国产精品麻豆| 久久精品国产综合久久久| 久久人妻福利社区极品人妻图片| 欧美日韩瑟瑟在线播放| 18美女黄网站色大片免费观看| 国产精品亚洲一级av第二区| 好男人在线观看高清免费视频 | 露出奶头的视频| 久久精品国产综合久久久| 精品国产乱码久久久久久男人| 免费在线观看视频国产中文字幕亚洲| 久久久久久久久久久久大奶| 亚洲国产高清在线一区二区三 | 大码成人一级视频| 黄色女人牲交| 日韩大尺度精品在线看网址 | 久久久久国产精品人妻aⅴ院| 一个人观看的视频www高清免费观看 | 欧美黄色片欧美黄色片| 成人欧美大片| 久久人妻福利社区极品人妻图片| 激情在线观看视频在线高清| 一区二区三区精品91| 每晚都被弄得嗷嗷叫到高潮| 99久久精品国产亚洲精品| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品久久电影中文字幕| 欧美激情 高清一区二区三区| 国内精品久久久久精免费| 久久久久久国产a免费观看| 国产亚洲精品av在线| 18禁美女被吸乳视频| 手机成人av网站| 国产精品久久视频播放| 亚洲国产欧美日韩在线播放| 两人在一起打扑克的视频| 国产xxxxx性猛交| 在线播放国产精品三级| 日本精品一区二区三区蜜桃| 午夜日韩欧美国产| 99国产精品一区二区蜜桃av| 黄网站色视频无遮挡免费观看| 男人操女人黄网站| 国产三级在线视频| 国产成人啪精品午夜网站| 亚洲,欧美精品.| 怎么达到女性高潮| 久久中文字幕人妻熟女| 欧美不卡视频在线免费观看 | 日本免费a在线| 色婷婷久久久亚洲欧美| 大香蕉久久成人网| 国产欧美日韩一区二区三区在线| 性欧美人与动物交配| 一区福利在线观看| 亚洲欧美日韩高清在线视频| 亚洲性夜色夜夜综合| 久久久久国产精品人妻aⅴ院| 久久久久国内视频| 国产欧美日韩一区二区三区在线| 中文字幕色久视频| 18禁国产床啪视频网站| 男男h啪啪无遮挡| а√天堂www在线а√下载| 亚洲欧美精品综合一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲自拍偷在线| 日韩欧美在线二视频| 亚洲成av片中文字幕在线观看| 色综合婷婷激情| av视频在线观看入口| 亚洲人成网站在线播放欧美日韩| 黄色片一级片一级黄色片| 成人三级黄色视频| 欧美国产日韩亚洲一区| 99在线人妻在线中文字幕| 精品国产国语对白av| 色播亚洲综合网| 成人特级黄色片久久久久久久| 国产麻豆成人av免费视频| 亚洲第一av免费看| 黄片播放在线免费| 久久中文字幕一级| 国产视频一区二区在线看| 国产精品美女特级片免费视频播放器 | 亚洲人成电影免费在线| 国产成人系列免费观看| 两个人看的免费小视频| 90打野战视频偷拍视频| 又黄又爽又免费观看的视频| 精品欧美国产一区二区三| 涩涩av久久男人的天堂| 国产高清videossex| 国产三级黄色录像| 亚洲最大成人中文| 在线观看免费日韩欧美大片| 免费看a级黄色片| 亚洲第一电影网av| 国产精品日韩av在线免费观看 | 欧美色视频一区免费| 国产野战对白在线观看| 国产精品一区二区在线不卡| 免费高清在线观看日韩| 在线观看66精品国产| 国产麻豆成人av免费视频| 久久狼人影院| 久久青草综合色| 亚洲九九香蕉| 日韩中文字幕欧美一区二区| 美女 人体艺术 gogo| 亚洲精品久久成人aⅴ小说| 中文字幕久久专区| 欧美成人性av电影在线观看| 999精品在线视频| 久久久国产成人精品二区| 欧美激情 高清一区二区三区| 亚洲av片天天在线观看| 18美女黄网站色大片免费观看| 男人舔女人的私密视频| 多毛熟女@视频| 久久婷婷人人爽人人干人人爱 | 黄色女人牲交| 国产亚洲欧美精品永久| 亚洲国产精品999在线| 操美女的视频在线观看| 日韩精品青青久久久久久| 一进一出抽搐动态| 韩国精品一区二区三区| 桃色一区二区三区在线观看| 久久天堂一区二区三区四区| 在线观看舔阴道视频| 99riav亚洲国产免费| 国产av在哪里看| 正在播放国产对白刺激| av在线天堂中文字幕| 国产精品亚洲一级av第二区| 操美女的视频在线观看| 免费在线观看完整版高清| 一本综合久久免费| 欧美乱妇无乱码| 国产在线观看jvid| 欧美最黄视频在线播放免费| 色播亚洲综合网| 亚洲 国产 在线| 国内久久婷婷六月综合欲色啪| 色综合站精品国产| 侵犯人妻中文字幕一二三四区| 女警被强在线播放| 久久精品国产亚洲av香蕉五月| 看黄色毛片网站| 一边摸一边抽搐一进一小说| 人人妻人人澡欧美一区二区 | 欧美日韩中文字幕国产精品一区二区三区 | 在线观看一区二区三区| 国产精品野战在线观看| 老鸭窝网址在线观看| 久久欧美精品欧美久久欧美| av天堂久久9| 91av网站免费观看| 在线观看日韩欧美| 一区二区日韩欧美中文字幕| 亚洲国产欧美网| 国产成人影院久久av| 国产三级在线视频| 黄色a级毛片大全视频| 久久精品国产亚洲av高清一级| 99久久久亚洲精品蜜臀av| 亚洲男人的天堂狠狠| 色综合站精品国产| 一边摸一边做爽爽视频免费| 久久久国产成人精品二区| 一二三四社区在线视频社区8| 99国产精品一区二区三区| 别揉我奶头~嗯~啊~动态视频| 国产主播在线观看一区二区| 色播在线永久视频| 亚洲五月婷婷丁香| 国产在线观看jvid| 国产成人免费无遮挡视频| 中国美女看黄片| 亚洲成人精品中文字幕电影| 久久青草综合色| 最近最新免费中文字幕在线| 国产亚洲精品第一综合不卡| 怎么达到女性高潮| 99国产精品一区二区蜜桃av| 欧美一区二区精品小视频在线| 久久伊人香网站| 极品教师在线免费播放| 色哟哟哟哟哟哟| 欧洲精品卡2卡3卡4卡5卡区| 黄色视频,在线免费观看| 国产一级毛片七仙女欲春2 | 国产精品免费视频内射| 桃色一区二区三区在线观看| 一二三四社区在线视频社区8| 国产乱人伦免费视频| 免费在线观看亚洲国产| 日本精品一区二区三区蜜桃| 手机成人av网站| 97超级碰碰碰精品色视频在线观看| 中文字幕人妻丝袜一区二区| 午夜激情av网站| 国产区一区二久久| 男女床上黄色一级片免费看| 欧美日韩精品网址| 国产区一区二久久| 19禁男女啪啪无遮挡网站| 国产视频一区二区在线看| 精品欧美一区二区三区在线| 成人亚洲精品一区在线观看| 亚洲精品在线美女| 熟女少妇亚洲综合色aaa.| 露出奶头的视频| 国产一区二区三区视频了| 国产精品亚洲av一区麻豆| 久久精品aⅴ一区二区三区四区| 国内精品久久久久精免费| 极品教师在线免费播放| 黄色成人免费大全| 首页视频小说图片口味搜索|