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

    基于CEEMD和遷移學(xué)習(xí)的滾動(dòng)軸承故障診斷研究

    2024-04-15 20:11:52張潤(rùn)地劉雨暉荊曉遠(yuǎn)韓光信
    河南科技 2024年4期
    關(guān)鍵詞:遷移學(xué)習(xí)滾動(dòng)軸承故障診斷

    張潤(rùn)地 劉雨暉 荊曉 遠(yuǎn)韓光信

    摘 要:【目的】在實(shí)際生產(chǎn)環(huán)境中,由于機(jī)器特征復(fù)雜和工況變化,智能診斷模型在跨機(jī)組遷移時(shí)需要重復(fù)訓(xùn)練,這不僅增加了時(shí)間成本,還加大了算力資源的消耗。為了解決這些問(wèn)題,需要開(kāi)發(fā)出一種能適應(yīng)復(fù)雜機(jī)器特征并在不同工況下保持高準(zhǔn)確度的軸承故障診斷方法,同時(shí),減少模型遷移時(shí)所需的重復(fù)訓(xùn)練,以便實(shí)現(xiàn)更高效的故障識(shí)別和預(yù)測(cè)。【方法】研究提出基于互補(bǔ)集合經(jīng)驗(yàn)?zāi)B(tài)分解(Complementary Ensemble Empirical Mode Decomposition,CEEMD)和遷移學(xué)習(xí)的滾動(dòng)軸承故障診斷方法。首先,采用CEEMD法對(duì)原始信號(hào)進(jìn)行分解,并計(jì)算出對(duì)應(yīng)分量的峭度值。其次,采用多核最大均值差異法對(duì)源域數(shù)據(jù)與目標(biāo)域數(shù)據(jù)進(jìn)行域適應(yīng)處理。最后,在凱斯西儲(chǔ)大學(xué)軸承數(shù)據(jù)集和美國(guó)機(jī)械故障預(yù)防技術(shù)學(xué)會(huì)軸承數(shù)據(jù)集之間進(jìn)行遷移故障診斷試驗(yàn)及對(duì)比分析?!窘Y(jié)果】研究結(jié)果表明,與直接遷移模型算法相比,基于CEEMD改進(jìn)的遷移學(xué)習(xí)網(wǎng)絡(luò)在不同數(shù)據(jù)集上的遷移效果更好,其故障診斷的準(zhǔn)確率最高?!窘Y(jié)論】經(jīng)試驗(yàn)驗(yàn)證,研究所提的方法表現(xiàn)出良好的變工況跨機(jī)組適配能力,具有較高的故障診斷精度,為研究復(fù)雜工況下多機(jī)組相似故障診斷場(chǎng)景提供了非常有價(jià)值的參考。

    關(guān)鍵詞:滾動(dòng)軸承;互補(bǔ)集合經(jīng)驗(yàn)?zāi)B(tài)分解;遷移學(xué)習(xí);故障診斷

    中圖分類(lèi)號(hào):TP391.4? ? ?文獻(xiàn)標(biāo)志碼:A? ? ?文章編號(hào):1003-5168(2024)04-0019-07

    DOI:10.19968/j.cnki.hnkj.1003-5168.2024.04.004

    Research on Rolling Bearing Fault Diagnosis Based on CEEMD and Transfer Learning

    ZHANG Rundi1 LIU Yuhui2 JING Xiaoyuan2 HAN Guangxin1

    (1.College of Information and Control Engineering, Jilin of Chemical Technology, Jilin 132022, China;

    2.College of Computer Science, Guangdong University of Petrochemical Technology, Maoming 525000, China)

    Abstract:[Purposes] In the actual production environment, due to the complexity of machine features and the change of working conditions, the intelligent diagnosis model needs repeated training when migrating across units, which not only increases the time cost, but also increases the consumption of computing resources. In order to solve these problems, it is necessary to develop a bearing fault diagnosis method that can adapt to complex machine features and maintain high accuracy under different working conditions. At the same time, the repeated training required for model migration is reduced to achieve more efficient fault identification and prediction. [Methods] The study proposes a rolling bearing fault diagnosis method based on CEEMD and transfer learning. First, the CEEMD decomposition method is used to decompose the original signal and the kurtosis value of the corresponding component is calculated. Then, the multi-core maximum mean difference method is used for the source domain data and the target domain data. Domain adaptation processing, and finally a migration fault diagnosis test and comparative analysis between the Case Western Reserve University dataset and the American Society for Mechanical Failure Prevention Technology dataset. [Findings]The research results show that compared with the existing direct transfer model algorithm, the improved transfer learning network based on CEEMD has a better transfer effect on different data sets, and its fault diagnosis accuracy is the highest. [Conclusions]It is verified by experiments that the method proposed in the study shows good cross-unit adaptability under variable working conditions, and has high fault diagnosis accuracy, which provides a valuable reference for studying similar fault diagnosis scenarios of multiple units under complex working conditions.

    Keywords: rolling bearing; CEEMD; transfer learning; fault diagnosis

    0 引言

    在實(shí)際工業(yè)生產(chǎn)中,滾動(dòng)軸承是旋轉(zhuǎn)機(jī)械的關(guān)鍵部件,一旦滾動(dòng)軸承發(fā)生故障,極有可能會(huì)對(duì)設(shè)備造成不可修復(fù)的重大損害。若滾動(dòng)軸承長(zhǎng)期處于故障狀態(tài),不僅會(huì)影響設(shè)備的安全性和可靠性,嚴(yán)重時(shí)還會(huì)造成重大的安全事故[1]。因此,研究滾動(dòng)軸承故障診斷技術(shù)對(duì)保障工業(yè)安全生產(chǎn)具有重要意義。

    隨著現(xiàn)代工業(yè)機(jī)械日益精細(xì)化和復(fù)雜化,提取設(shè)備故障特征也變得越來(lái)越困難,如何在減少故障診斷成本的基礎(chǔ)上,對(duì)故障特征實(shí)現(xiàn)準(zhǔn)確高效的提取,成為跨設(shè)備故障診斷工作的重要一環(huán)?,F(xiàn)有的故障特征提取大多是基于恒定工況下的信號(hào)處理方式,而在實(shí)際工業(yè)生產(chǎn)環(huán)境中,不同時(shí)刻各種機(jī)械裝備的轉(zhuǎn)速與工作負(fù)載是不同的,加上各種部件之間的耦合,導(dǎo)致真實(shí)運(yùn)行狀態(tài)的振動(dòng)信號(hào)往往都是非線性、非平穩(wěn)的。目前,常用的故障特征提取方法有小波變換(Wavelet Transform, WT)、經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition, EMD)。石志煒等[2]采用改進(jìn)小波包閾值降噪法,根據(jù)信號(hào)的噪聲比例來(lái)改動(dòng)閾值函數(shù),從而實(shí)現(xiàn)故障特征的提??;Huang等[3]提出的EMD是無(wú)須先驗(yàn)知識(shí)就能自適應(yīng)分解非線性、非平穩(wěn)時(shí)頻信號(hào)的方法,但存在無(wú)法準(zhǔn)確分解信號(hào)的缺陷;針對(duì)EMD分解存在的模態(tài)疊混的問(wèn)題,Wu等[4]提出集合經(jīng)驗(yàn)?zāi)B(tài)分解(Ensemble Empirical Mode Decomposition,EEMD)方法,EEMD方法利用輔助噪聲,有效避免了模態(tài)混疊問(wèn)題,但其存在計(jì)算量大、分解不完備的缺陷。基于此,Yeh等[5]提出了互補(bǔ)集合經(jīng)驗(yàn)?zāi)B(tài)分解(Complementary Ensemble Empirical Mode Decomposition,CEEMD)的方法。與EEMD相比,CEEMD的迭代次數(shù)更少,且信號(hào)分解效果更好。國(guó)內(nèi)研究學(xué)者在上述信號(hào)特征提取理論的基礎(chǔ)上,進(jìn)行滾動(dòng)軸承故障診斷的前期研究工作。張萍等[6]采用互補(bǔ)集合經(jīng)驗(yàn)?zāi)B(tài)分解法對(duì)滾動(dòng)軸承故障進(jìn)行診斷,能有效地彌補(bǔ)EMD分解不準(zhǔn)確的缺陷,但其準(zhǔn)確度還是稍顯不足;王夢(mèng)[7]利用EEMD與峭度譜結(jié)合的方法對(duì)旋轉(zhuǎn)機(jī)械故障進(jìn)行診斷,該方法故障特征提取的準(zhǔn)確度比EEMD法有所提高,但信號(hào)分解的計(jì)算量卻增大;宋治惠[8]提出CEEMD和峭度譜結(jié)合的方法,能準(zhǔn)確有效地提取出滾動(dòng)軸承的早期故障特征,并能有效減少計(jì)算量。

    基于傳統(tǒng)機(jī)器學(xué)習(xí)的故障診斷方法雖然有著較低的計(jì)算成本,但隨著工業(yè)互聯(lián)網(wǎng)的接入與企業(yè)的數(shù)字化轉(zhuǎn)型,已不能適應(yīng)當(dāng)前擁有海量數(shù)據(jù)信息的智慧化工業(yè)生產(chǎn)過(guò)程。由此采用深度學(xué)習(xí)的方法,通過(guò)多層次神經(jīng)元網(wǎng)絡(luò)進(jìn)行特征提取,彌補(bǔ)了傳統(tǒng)機(jī)器學(xué)習(xí)的缺點(diǎn),并取得更好的成效。

    基于深度學(xué)習(xí)的方法具有強(qiáng)大的特征自適應(yīng)提取能力,但其分類(lèi)性能受到網(wǎng)絡(luò)層數(shù)的影響。網(wǎng)絡(luò)層數(shù)越多,所提取到的特征信息越多,計(jì)算復(fù)雜度也就越高。但一味簡(jiǎn)單地疊加神經(jīng)網(wǎng)絡(luò)會(huì)造成梯度彌散或梯度爆炸等問(wèn)題,導(dǎo)致模型的識(shí)別精度不高,訓(xùn)練時(shí)間也會(huì)大大增加。此外,在變工況跨設(shè)備的相似故障診斷場(chǎng)景中,常用的深度學(xué)習(xí)方法仍難以準(zhǔn)確識(shí)別出故障類(lèi)型。引入遷移學(xué)習(xí)能在減少訓(xùn)練成本的同時(shí),大大提高故障診斷的準(zhǔn)確率。因此,一些學(xué)者也對(duì)遷移學(xué)習(xí)應(yīng)用于故障診斷領(lǐng)域進(jìn)行了探索。Yang等[9]提出一種基于特征遷移的故障診斷模型,通過(guò)引入多層最大均值差異來(lái)實(shí)現(xiàn)特征的域適應(yīng),實(shí)現(xiàn)了跨設(shè)備故障遷移診斷;為了解決數(shù)據(jù)稀疏問(wèn)題,Han等[10]提出一種基于生成對(duì)抗網(wǎng)絡(luò)的故障診斷方法,利用多個(gè)并行域鑒別器進(jìn)行數(shù)據(jù)匹配,再單獨(dú)進(jìn)行域適應(yīng)處理,從而完成故障診斷。然而,現(xiàn)有基于遷移學(xué)習(xí)的故障診斷研究仍存在兩個(gè)問(wèn)題:一是用于遷移的模型大都需要多次重復(fù)訓(xùn)練,而隨著網(wǎng)絡(luò)模型的加深,其所需要的標(biāo)記數(shù)據(jù)和訓(xùn)練設(shè)備資源也會(huì)增多,導(dǎo)致模型訓(xùn)練成本有所增加;二是目前的研究主要處理的是一維信號(hào)的輸入數(shù)據(jù),其樣本特征存在單一性問(wèn)題,導(dǎo)致模型無(wú)法提取到足夠多的特征用于故障識(shí)別。

    綜上所述,針對(duì)實(shí)際生產(chǎn)環(huán)境中存在的復(fù)雜機(jī)器特征難以完整提取、變工況跨機(jī)組遷移時(shí)需要模型重復(fù)訓(xùn)練等問(wèn)題,本研究提出一種基于CEEMD和遷移學(xué)習(xí)的滾動(dòng)軸承故障診斷方法。首先,在保證分解效果與EEMD相當(dāng)?shù)那疤嵯?,CEEMD能抑制由白噪聲引起的重構(gòu)誤差。其次,在CEEMD方法的基礎(chǔ)上,結(jié)合計(jì)算每個(gè)固有模態(tài)函數(shù)對(duì)應(yīng)的峭度值,提取出滾動(dòng)軸承故障樣本中更深層次的特征信息。最后,基于特征映射的遷移學(xué)習(xí)網(wǎng)絡(luò)模型,采用多核最大均值差異法來(lái)降低網(wǎng)絡(luò)訓(xùn)練成本,并提高網(wǎng)絡(luò)對(duì)故障識(shí)別的準(zhǔn)確率。因此,本研究基于該方法來(lái)構(gòu)建改進(jìn)的深度神經(jīng)網(wǎng)絡(luò)模型,以獲取更精準(zhǔn)、更全面的故障特征,并通過(guò)試驗(yàn)驗(yàn)證所提出的方法,并證明該方法能更好地適應(yīng)變工況跨機(jī)組的復(fù)雜場(chǎng)景。

    1 軸承故障診斷理論

    1.1 振動(dòng)信號(hào)處理

    在信號(hào)處理方面,有三種典型的處理方法:一是時(shí)域處理方法,如有量綱特征值分析和無(wú)量綱特征值分析;二是頻域處理方法,如傅里葉變換、能量譜和功率譜;三是時(shí)頻域分析,如小波分析。實(shí)際工業(yè)生產(chǎn)中所采集到的信號(hào)通常為時(shí)域信號(hào),為了平衡特征提取工作的成本與準(zhǔn)確性,本研究選擇時(shí)頻域分析方法。

    時(shí)頻處理信號(hào)的方法有小波變換和EMD。小波變換是基于一個(gè)特定小波基的分析方法,一旦確定了小波基,那么在整個(gè)分析過(guò)程中都將無(wú)法對(duì)其進(jìn)行更換,即使該小波基在全局可能是最優(yōu)的,但在一些局部地方卻不是最優(yōu)的。因此,小波分析不能針對(duì)不同工況和不同機(jī)組之間的情況進(jìn)行適應(yīng)調(diào)整。而EMD最突出的特點(diǎn)是解決了小波變換中基函數(shù)不能自適應(yīng)的問(wèn)題,對(duì)一段未知的非線性、非平穩(wěn)信號(hào),EMD不需要去做預(yù)先分析與研究,可以直接開(kāi)始分解。EMD的分解過(guò)程如下:①根據(jù)原始信號(hào)的上下極值點(diǎn),分別畫(huà)出上下包絡(luò)線;②求上下包絡(luò)線的均值,并畫(huà)出均值包絡(luò)線;③原始信號(hào)減去均值包絡(luò)線,得到中間信號(hào);④判斷該中間信號(hào)是否為IMF,如果不是,以該信號(hào)為基礎(chǔ),重新進(jìn)行①到④的操作;⑤在使用上述方法得到第一個(gè)IMF后(即IMF1),用原始信號(hào)減IMF1 作為新的原始信號(hào),再通過(guò)①到④的操作,得到IMF2,以此類(lèi)推,完成EMD分解。

    在實(shí)際的工業(yè)機(jī)器振動(dòng)信號(hào)采集過(guò)程中,一旦出現(xiàn)間歇性信號(hào),采用EMD分解法會(huì)出現(xiàn)頻率混疊的現(xiàn)象,也叫模態(tài)混疊。具體表現(xiàn)為一個(gè)IMF分量中存在多個(gè)尺度成分或一個(gè)尺度成分在多個(gè)IMF分量中存在。

    為了解決這個(gè)問(wèn)題,可以采用噪聲輔助處理的方法,即EEMD。EEMD利用白噪聲均值為0的特性,在信號(hào)分解過(guò)程中,通過(guò)多次引入均勻分布的白噪聲來(lái)掩蓋原信號(hào)本身的噪聲,從而得到更加精準(zhǔn)的上下包絡(luò)線。EEMD方法的分解過(guò)程如下:①設(shè)定原始信號(hào)的處理次數(shù)(m);②給這m個(gè)原始信號(hào)分別添加隨機(jī)白噪聲,并組成一系列新的信號(hào);③對(duì)這一系列新信號(hào)分別進(jìn)行EMD分解,得到一系列的IMF分量;④對(duì)相應(yīng)模態(tài)的IMF分量分別求均值,得到EEMD分解結(jié)果。

    CEEMD是在EMD和EEMD的基礎(chǔ)上,對(duì)信號(hào)分解過(guò)程中存在的模態(tài)混疊問(wèn)題進(jìn)行的改進(jìn)。CEEMD以噪聲輔助數(shù)據(jù)分析方法為基礎(chǔ),在信號(hào)中添加多對(duì)符號(hào)相反的白噪聲后,再應(yīng)用EMD分解,對(duì)分解結(jié)果進(jìn)行集成平均運(yùn)算,從而獲得最終的分解結(jié)果。

    在原始信號(hào)中添加一對(duì)符號(hào)相反的噪聲信號(hào),加入的白噪聲幅值相同,其表現(xiàn)形式見(jiàn)式(1)、式(2)。

    [m+it=xt+n+it] (1)

    [m-it=xt+n-it] (2)

    式中:[xt]為原始信號(hào);[n+it]為正噪聲;[n-it]為負(fù)噪聲。

    首先,對(duì)[m+it]、[m-it]進(jìn)行EMD分解,得到兩組集成的IMF分量IMF1、IMF2。其中,IMF1為加入正噪聲的集成平均結(jié)果,IMF2為加入負(fù)噪聲的集成平均結(jié)果。其次,將每組IMF函數(shù)進(jìn)行組合加權(quán)處理,生成一組總IMF函數(shù)。最后,對(duì)總IMF函數(shù)進(jìn)行EMD分解,得到一組新的IMF函數(shù)。如此反復(fù)迭代,直到得到的IMF函數(shù)的數(shù)量不再增加為止。

    在實(shí)際生產(chǎn)過(guò)程中,機(jī)組不可能維持在一個(gè)恒定的轉(zhuǎn)速、頻率下工作,因此,為了適應(yīng)實(shí)際情況下的機(jī)械故障診斷,本研究還需要對(duì)采集到的時(shí)域信號(hào)進(jìn)行無(wú)量綱指標(biāo)提取。主要選取對(duì)這個(gè)峭度信號(hào)瞬時(shí)沖擊極為敏感的無(wú)量綱指標(biāo),為完善故障信號(hào)的特征提取提供有效補(bǔ)充,進(jìn)而過(guò)濾CEEMD分解后的IMF分量,以提高振動(dòng)信號(hào)故障特征提取效率。峭度值的計(jì)算見(jiàn)式(3)。

    [k4x=Ex-μ4σ4] (3)

    式中:[E]為求期望值;[x]為信號(hào)值;[μ]為信號(hào)均值;[σ]為標(biāo)準(zhǔn)差。

    1.2 遷移學(xué)習(xí)及其網(wǎng)絡(luò)模型

    卷積神經(jīng)網(wǎng)絡(luò)(Convolutional Neural Networks,CNN)是目前深度學(xué)習(xí)理論中應(yīng)用最為成熟的網(wǎng)絡(luò)模型,因其具有強(qiáng)大的特征提取和模型泛化能力而被廣泛應(yīng)用于多個(gè)領(lǐng)域。目前,軸承故障診斷領(lǐng)域也有大量成果是基于CNN網(wǎng)絡(luò)模型演化而來(lái)的。CNN通常由卷積層、池化層和全連接層構(gòu)成。其中,卷積層用于提取特征,是CNN網(wǎng)絡(luò)的核心;池化層對(duì)經(jīng)過(guò)卷積處理得到的特征圖進(jìn)行特征降維,壓縮數(shù)據(jù)和參數(shù)的數(shù)量,同時(shí)減小網(wǎng)絡(luò)參數(shù)規(guī)模,避免過(guò)擬合;全連接層將學(xué)習(xí)到的特征進(jìn)行輸出。

    隨著深度學(xué)習(xí)在軸承故障診斷領(lǐng)域中的廣泛應(yīng)用,在恒定工況下運(yùn)行的單機(jī)器表現(xiàn)良好的CNN網(wǎng)絡(luò)模型,卻在變工況多機(jī)組場(chǎng)景下顯現(xiàn)出了訓(xùn)練耗時(shí)、準(zhǔn)確率低下等弊端。因此,需要將遷移學(xué)習(xí)方法應(yīng)用于軸承故障診斷領(lǐng)域中,來(lái)解決上述問(wèn)題。

    在上述問(wèn)題場(chǎng)景中,已訓(xùn)練好的模型不能很好地適應(yīng)新的目標(biāo)機(jī)器,需要對(duì)不同的目標(biāo)數(shù)據(jù)域重新訓(xùn)練模型,從而導(dǎo)致學(xué)習(xí)成本太高,因此,本研究運(yùn)用已有的知識(shí)或模型將訓(xùn)練好的模型遷移到目標(biāo)域上。根據(jù)遷移方式來(lái)劃分,遷移學(xué)習(xí)方法可分為基于實(shí)例的遷移方法、基于映射的遷移方法、基于模型的遷移方法、基于對(duì)抗的遷移方法等。

    本研究采用基于映射的遷移方法,將深度學(xué)習(xí)與遷移學(xué)習(xí)相結(jié)合,構(gòu)建出深度自適應(yīng)網(wǎng)絡(luò),核心是使用設(shè)計(jì)的自適應(yīng)層來(lái)完成源域與目標(biāo)域數(shù)據(jù)的自適應(yīng),其中,自適應(yīng)層采用多核最大均值差異(Multiple Kernel Maximum Mean Discrepancy,MK-MMD)方法。最大均值差異(Maximum Mean Discrepancy,MMD)方法是遷移學(xué)習(xí)中實(shí)現(xiàn)域適應(yīng)技術(shù)最常用的一種損失函數(shù),其定義見(jiàn)式(4)。

    [MF, p, q?supf?FEpfx-Eqfy]? (4)

    式中:[sup]為求上界,即最大值;[Ep]為求期望值;[f?]為映射函數(shù);[F]為函數(shù)域;[x]為源域數(shù)據(jù);[y]為目標(biāo)域數(shù)據(jù);x和y的分布分別為[p]和[q];[fx]、[fy]分別為源域和目標(biāo)域數(shù)據(jù)經(jīng)過(guò)映射后的值。

    MK-MMD是在MMD基礎(chǔ)上發(fā)展而來(lái)的,用多個(gè)核去構(gòu)造一個(gè)總的核,定義見(jiàn)式(5)。

    [d2kp,q?Ep?xs-Eq?xt2Hk] (5)

    式中:[d2kp,q]為[p]和[q]之間的再生核希爾伯特空間距離,其中,將核(Kernel)定義為[K],表示見(jiàn)式(6)。

    [K?k=u=1mβuki:βu≥0,?u] (6)

    式中:[βu]為權(quán)重。

    本研究使用的損失函數(shù)定義見(jiàn)式(7)。

    [l=lcDs, ys+λMK-MMD2Ds,Dt] (7)

    式中:[l]為總損失;[lcDs,ys]為有標(biāo)注數(shù)據(jù)上的常規(guī)損失;[λMK-MMD2Ds,Dt]為模型的多核最大均值差異損失;[Ds]為源域;[Dt]為目標(biāo)域;[ys]為源域的標(biāo)簽數(shù)據(jù);[λ]為權(quán)重。

    2 網(wǎng)絡(luò)模型的改進(jìn)與故障診斷方法

    2.1 改進(jìn)的AlexNet模型

    AlexNet是在2012年ImageNet圖像分類(lèi)競(jìng)賽中提出的一種經(jīng)典的卷積神經(jīng)網(wǎng)絡(luò),共包含5個(gè)卷積層和3個(gè)全連接層。其中,第一、二、五個(gè)卷積層后面都附加一個(gè)池化層,后三個(gè)層為全連接層。

    本研究所提出的改進(jìn)模型(以下稱(chēng)為FineTune-AlexNet)是在AlexNet基礎(chǔ)上,添加不適用激活函數(shù)的線性輸出層,并將最后幾個(gè)全連接層設(shè)置為域適應(yīng)層。該域適應(yīng)層使用MK-MMD作為域適應(yīng)的正則化項(xiàng),以減少源域和目標(biāo)域之間的分布差異。改進(jìn)后的網(wǎng)絡(luò)模型層級(jí)如圖1所示。

    2.2 基于CEEMD與深度自適應(yīng)網(wǎng)絡(luò)的軸承故障診斷框架

    本研究提出的軸承故障診斷方法流程如圖2所示。包括信號(hào)數(shù)據(jù)預(yù)處理、模型特征映射遷移學(xué)習(xí)、使用訓(xùn)練好的模型進(jìn)行滾動(dòng)軸承故障診斷等。

    2.2.1 信號(hào)預(yù)處理。先使用CEEMD法對(duì)機(jī)器原始時(shí)域振動(dòng)信號(hào)進(jìn)行分解,再求取IMF分量的峭度值,用其組成一組新的數(shù)據(jù)。

    2.2.2 模型特征映射遷移。先對(duì)源域數(shù)據(jù)進(jìn)行模型訓(xùn)練,再利用遷移學(xué)習(xí)方法將訓(xùn)練后的權(quán)重參數(shù)針對(duì)目標(biāo)域的數(shù)據(jù)進(jìn)行特征映射遷移,直到損失收斂完成遷移學(xué)習(xí)為止。

    2.2.3 故障診斷。使用遷移學(xué)習(xí)后的模型進(jìn)行故障診斷和診斷結(jié)果正確率的精度驗(yàn)證。

    3 試驗(yàn)分析

    3.1 數(shù)據(jù)集和相關(guān)參數(shù)

    本研究選取凱斯西儲(chǔ)大學(xué)軸承數(shù)據(jù)集和美國(guó)-機(jī)械故障預(yù)防技術(shù)學(xué)會(huì)的軸承數(shù)據(jù)集(以下分別簡(jiǎn)稱(chēng)CWRU數(shù)據(jù)集和MFPT數(shù)據(jù)集)分別作為試驗(yàn)的源域數(shù)據(jù)和目標(biāo)域數(shù)據(jù)來(lái)源進(jìn)行驗(yàn)證。試驗(yàn)所使用的機(jī)器如圖3所示。主要由風(fēng)扇端軸承、驅(qū)動(dòng)端軸承、電機(jī)、扭矩傳感器、編碼器和負(fù)載組成。

    試驗(yàn)所使用的CWRU源域數(shù)據(jù)集包括四種軸承健康狀態(tài)。即正常狀態(tài)、內(nèi)圈故障、外圈故障和滾動(dòng)體故障,其中,故障直徑有0.017 78 cm、0.035 56 cm、0.053 34 cm這三種尺寸類(lèi)型。試驗(yàn)電機(jī)是在0 kW、0.735 kW、1.470 kW、2.205 kW四種不同負(fù)載和1 730 r/min、1 750 r/min、1 772 r/min、1 797 r/min四種不同轉(zhuǎn)速下工作的。

    試驗(yàn)所使用的MFPT目標(biāo)域數(shù)據(jù)集包括三種軸承健康狀態(tài),即正常狀態(tài)、內(nèi)圈故障、外圈故障。根據(jù)不同的工作狀況,數(shù)據(jù)集包含以采樣率97 656 Hz和負(fù)載122.47 kg的基線條件的3種健康數(shù)據(jù)和3種外圈故障數(shù)據(jù)在11.34 kg、22.68 kg、45.36 kg、68.04 kg、90.72 kg、113.40 kg和136.08 kg負(fù)載下的7種外圈故障數(shù)據(jù)及在0 kg、22.68 kg、45.36 kg、68.04 kg、90.72 kg、113.40 kg和136.08 kg負(fù)載下的7種內(nèi)圈故障數(shù)據(jù)。

    3.2 試驗(yàn)驗(yàn)證及結(jié)果分析

    首先,對(duì)從上述兩個(gè)數(shù)據(jù)集提取的原始振動(dòng)信號(hào)數(shù)據(jù)進(jìn)行CEEMD分解預(yù)處理,分解后的時(shí)域IMF分量圖像示例如圖4所示。計(jì)算出每個(gè)分量所對(duì)應(yīng)的峭度值,再按照表1的數(shù)量以8∶2的比例劃分為訓(xùn)練集和驗(yàn)證集。數(shù)據(jù)信號(hào)樣本分類(lèi)及數(shù)量見(jiàn)表1。

    其次,將經(jīng)過(guò)CEEMD算法處理后的分量信號(hào)及每條分量信號(hào)對(duì)應(yīng)的峭度值按順序排列,并重新構(gòu)建源域和目標(biāo)域的數(shù)據(jù)集。

    最后,分別按照是否進(jìn)行信號(hào)預(yù)處理和是否進(jìn)行改進(jìn)的模型訓(xùn)練進(jìn)行消融試驗(yàn)。試驗(yàn)對(duì)比結(jié)果見(jiàn)表2,模型準(zhǔn)確率與損失迭代的對(duì)比曲線如圖5所示。

    由表2和圖5可知,相較于直接使用AlexNet模型進(jìn)行不同機(jī)器間的遷移訓(xùn)練,將經(jīng)過(guò)CEEMD信號(hào)預(yù)處理的訓(xùn)練模型應(yīng)用到目標(biāo)域的故障數(shù)據(jù)集中,準(zhǔn)確率提升了33%。而對(duì)AlexNet方法、“AlexNet+CEEMD”方法和FineTune-AlexNet方法進(jìn)行對(duì)比分析,不難發(fā)現(xiàn),上述三種方法準(zhǔn)確率均低于“FineTune-AlexNet+CEEMD”方法。這是因?yàn)楸狙芯刻岢龅幕谀繕?biāo)域自適應(yīng)的深度神經(jīng)網(wǎng)絡(luò)模型能有效完成特征映射遷移,從而提高故障特征信息的提取效率和故障診斷準(zhǔn)確率。試驗(yàn)得到的結(jié)果也驗(yàn)證了本研究所提出的方法能更好地適應(yīng)變工況跨機(jī)組的復(fù)雜場(chǎng)景的結(jié)論。

    4 結(jié)論

    針對(duì)復(fù)雜工況下機(jī)器特征提取不完備而導(dǎo)致的故障診斷準(zhǔn)確率不夠高、復(fù)雜工況多機(jī)組需要多次模型訓(xùn)練等問(wèn)題,本研究提出一種基于CEEMD和遷移學(xué)習(xí)的滾動(dòng)軸承故障診斷方法。首先,對(duì)原始振動(dòng)信號(hào)進(jìn)行預(yù)處理,通過(guò)CEEMD分解提取峭度值,從而增強(qiáng)數(shù)據(jù)特征提取精度,并減少數(shù)據(jù)運(yùn)算成本;其次,將處理后的信號(hào)輸入深度遷移網(wǎng)絡(luò)模型進(jìn)行訓(xùn)練,并進(jìn)行故障診斷;最后,經(jīng)過(guò)試驗(yàn)證明,應(yīng)用該方法后,變工況跨數(shù)據(jù)集間遷移精度有很大提升,并降低了數(shù)據(jù)處理的成本。本研究所提出的方案在對(duì)變工況多機(jī)組軸承故障進(jìn)行診斷時(shí),比一些先進(jìn)的算法更為精確。本研究的研究結(jié)果可為實(shí)現(xiàn)智慧工廠的科學(xué)化決策提供理論參考與支持。

    參考文獻(xiàn):

    [1]黃健豪,鄭波,陳國(guó)慶.基于無(wú)監(jiān)督遷移學(xué)習(xí)的核范數(shù)最大化軸承故障診斷方法[J].科學(xué)技術(shù)與工程,2023(11):4638-4646.

    [2]石志煒,張麗萍.基于改進(jìn)小波包閾值降噪的滾動(dòng)軸承故障分析[J].計(jì)算機(jī)測(cè)量與控制,2019(5):58-63.

    [3]HUANG N E,SHEN Z,LONG S R,et al.The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[J].Proceedings of the Royal Society of London Series A,1998,454:903-995.

    [4]WU Z H,HUANG N E.Ensemble empirical mode decomposition:a noise-assisted data analysis method[J].Advances in Adaptive Data Analysis,2009(1):1-41.

    [5]YEH J R,SHIEN J S,HUANG N E.Complementary ensemble empirical mode decomposition:a novel noise enhanced data analysis method[J].Advances in Adaptive Data Analysis,2010(2):135-156.

    [6]張萍,李志農(nóng),陳靜鈴,等.基于互補(bǔ)集合平均經(jīng)驗(yàn)?zāi)B(tài)分解的滾動(dòng)軸承故障診斷方法研究[J].南昌航空大學(xué)學(xué)報(bào)(自然科學(xué)版),2019(2):7-12,49.

    [7]王夢(mèng).基于階次分析的旋轉(zhuǎn)機(jī)械故障診斷方法研究[D].秦皇島:燕山大學(xué),2014.

    [8]宋治惠.基于CEEMDAN的滾動(dòng)軸承早期故障特征提取方法研究[D].天津:天津大學(xué),2017.

    [9]YANG B,LEI Y G,JIA F,et al.An intelligent fault diagnosis approach based on transfer learning from laboratory bearings to locomotive bearings[J].Mechanical Systems and Signal Processing,2019,122:692-706.

    [10]HAN T,LIU C,WU R,et al.Deep transfer learning with limited data for machinery fault diagnosis[J].Applied Soft Computing Journal,2021,103:107150.

    猜你喜歡
    遷移學(xué)習(xí)滾動(dòng)軸承故障診斷
    基于多特征融合的跨域情感分類(lèi)模型研究
    奇異值分解與移移學(xué)習(xí)在電機(jī)故障診斷中的應(yīng)用
    滾動(dòng)軸承裝配方法及注意事項(xiàng)
    滾動(dòng)軸承安裝使用中的幾個(gè)問(wèn)題
    基于小波包絡(luò)分析的滾動(dòng)軸承典型故障診斷技術(shù)研究
    商情(2016年11期)2016-04-15 07:39:54
    因果圖定性分析法及其在故障診斷中的應(yīng)用
    一種基于遷移極速學(xué)習(xí)機(jī)的人體行為識(shí)別模型
    基于包絡(luò)分析法的軸承故障分析
    大數(shù)據(jù)環(huán)境下基于遷移學(xué)習(xí)的人體檢測(cè)性能提升方法
    基于LCD和排列熵的滾動(dòng)軸承故障診斷
    国产一区二区在线观看av| 日日摸夜夜添夜夜爱| 丰满迷人的少妇在线观看| 女性被躁到高潮视频| 18禁观看日本| 巨乳人妻的诱惑在线观看| 国产探花极品一区二区| 成年人免费黄色播放视频| 国产不卡av网站在线观看| 大香蕉久久成人网| 欧美黄色片欧美黄色片| 巨乳人妻的诱惑在线观看| 亚洲一区中文字幕在线| 亚洲一区中文字幕在线| 91国产中文字幕| 欧美变态另类bdsm刘玥| 久久国产精品男人的天堂亚洲| 蜜桃在线观看..| 看十八女毛片水多多多| 亚洲一码二码三码区别大吗| 高清视频免费观看一区二区| 欧美在线黄色| 伊人久久国产一区二区| 亚洲欧美一区二区三区国产| 一级毛片我不卡| 一二三四中文在线观看免费高清| 中文字幕av电影在线播放| 天天躁日日躁夜夜躁夜夜| 久久精品亚洲熟妇少妇任你| 如何舔出高潮| videos熟女内射| 国产一区有黄有色的免费视频| 午夜福利一区二区在线看| 人妻 亚洲 视频| 欧美最新免费一区二区三区| 一本一本久久a久久精品综合妖精| 考比视频在线观看| 久久国产精品男人的天堂亚洲| 伊人久久大香线蕉亚洲五| 男人操女人黄网站| 超色免费av| 只有这里有精品99| 麻豆av在线久日| 男女下面插进去视频免费观看| 91精品国产国语对白视频| 汤姆久久久久久久影院中文字幕| 岛国毛片在线播放| 宅男免费午夜| 久久久国产欧美日韩av| 操出白浆在线播放| 亚洲久久久国产精品| 18禁国产床啪视频网站| 亚洲视频免费观看视频| 国产av一区二区精品久久| 国产av一区二区精品久久| 日韩一区二区视频免费看| 亚洲第一青青草原| 高清不卡的av网站| 日本vs欧美在线观看视频| bbb黄色大片| 成人三级做爰电影| 亚洲av福利一区| 国产亚洲午夜精品一区二区久久| 制服人妻中文乱码| svipshipincom国产片| 日本爱情动作片www.在线观看| 热re99久久国产66热| 老司机深夜福利视频在线观看 | 亚洲精品成人av观看孕妇| 黄色 视频免费看| 天天躁狠狠躁夜夜躁狠狠躁| 久久久国产欧美日韩av| 亚洲第一区二区三区不卡| 最新在线观看一区二区三区 | 99久国产av精品国产电影| 亚洲国产精品国产精品| 在线天堂最新版资源| 久久精品亚洲av国产电影网| 一本色道久久久久久精品综合| 一本大道久久a久久精品| 一个人免费看片子| 满18在线观看网站| 一边亲一边摸免费视频| 日韩 亚洲 欧美在线| 赤兔流量卡办理| 国产黄色免费在线视频| 青春草亚洲视频在线观看| 国产亚洲一区二区精品| 日本黄色日本黄色录像| 日韩伦理黄色片| 国产免费视频播放在线视频| 亚洲精品国产av蜜桃| 久久鲁丝午夜福利片| 欧美日韩一级在线毛片| 男男h啪啪无遮挡| 中文字幕亚洲精品专区| 老司机影院毛片| 少妇人妻 视频| 男女国产视频网站| 极品人妻少妇av视频| 人体艺术视频欧美日本| 波野结衣二区三区在线| 秋霞伦理黄片| 高清不卡的av网站| 91精品国产国语对白视频| 日韩欧美精品免费久久| 欧美久久黑人一区二区| 人妻一区二区av| 尾随美女入室| 精品免费久久久久久久清纯 | 欧美人与性动交α欧美精品济南到| 国产麻豆69| 免费av中文字幕在线| 一区二区三区乱码不卡18| 国产亚洲欧美精品永久| 日韩av在线免费看完整版不卡| 永久免费av网站大全| 国产女主播在线喷水免费视频网站| 视频在线观看一区二区三区| 精品一区二区免费观看| 国产精品女同一区二区软件| 精品久久久精品久久久| 久久久久国产精品人妻一区二区| 九九爱精品视频在线观看| 日本欧美视频一区| 黄色一级大片看看| 午夜日本视频在线| 国产一卡二卡三卡精品 | 最近最新中文字幕大全免费视频 | 乱人伦中国视频| www日本在线高清视频| 亚洲成人一二三区av| 国产免费又黄又爽又色| 叶爱在线成人免费视频播放| 日韩av在线免费看完整版不卡| 国产精品免费大片| 在线 av 中文字幕| 丁香六月天网| 亚洲男人天堂网一区| 亚洲七黄色美女视频| 精品久久蜜臀av无| 高清欧美精品videossex| 丝瓜视频免费看黄片| 久久久欧美国产精品| 在线看a的网站| 一区在线观看完整版| 国产亚洲欧美精品永久| 久久精品国产综合久久久| 婷婷成人精品国产| 国产成人系列免费观看| 国产伦理片在线播放av一区| 老司机影院毛片| 99国产精品免费福利视频| 亚洲欧美精品综合一区二区三区| 亚洲五月色婷婷综合| 涩涩av久久男人的天堂| 精品国产一区二区三区四区第35| 成人亚洲欧美一区二区av| 汤姆久久久久久久影院中文字幕| 桃花免费在线播放| 成人手机av| 性少妇av在线| 青春草亚洲视频在线观看| 在现免费观看毛片| 自线自在国产av| 欧美黄色片欧美黄色片| 国产av精品麻豆| 99热国产这里只有精品6| 久久97久久精品| 日韩一区二区三区影片| 国产一区有黄有色的免费视频| 90打野战视频偷拍视频| 久久精品亚洲av国产电影网| 国产精品欧美亚洲77777| 美国免费a级毛片| 青春草国产在线视频| 秋霞伦理黄片| 亚洲欧美一区二区三区黑人| 99热网站在线观看| 久久韩国三级中文字幕| 99久久人妻综合| 看免费成人av毛片| 日本av免费视频播放| 午夜免费鲁丝| 日本欧美视频一区| 国产黄频视频在线观看| 91精品伊人久久大香线蕉| 天堂8中文在线网| 国产黄频视频在线观看| 亚洲精品国产色婷婷电影| 国语对白做爰xxxⅹ性视频网站| 女性被躁到高潮视频| 午夜福利影视在线免费观看| 国产成人免费无遮挡视频| 欧美成人午夜精品| kizo精华| 亚洲人成电影观看| 日韩熟女老妇一区二区性免费视频| 99精品久久久久人妻精品| 毛片一级片免费看久久久久| 国产精品久久久久久久久免| 久久鲁丝午夜福利片| 在线天堂最新版资源| 七月丁香在线播放| 亚洲国产欧美网| 婷婷成人精品国产| 亚洲综合精品二区| 人妻人人澡人人爽人人| 如何舔出高潮| 国产av一区二区精品久久| 亚洲av欧美aⅴ国产| 欧美日韩一级在线毛片| 欧美激情高清一区二区三区 | av电影中文网址| 午夜影院在线不卡| 人人妻,人人澡人人爽秒播 | 男人舔女人的私密视频| 波野结衣二区三区在线| 日韩大片免费观看网站| 91精品伊人久久大香线蕉| 国产一区有黄有色的免费视频| 在线观看一区二区三区激情| av网站免费在线观看视频| 日本wwww免费看| 99精国产麻豆久久婷婷| 欧美国产精品一级二级三级| 午夜福利影视在线免费观看| 国产深夜福利视频在线观看| 老司机影院毛片| 精品久久蜜臀av无| www日本在线高清视频| 亚洲人成网站在线观看播放| 最新的欧美精品一区二区| 国产成人免费无遮挡视频| 97在线人人人人妻| 久久精品人人爽人人爽视色| 精品一区二区免费观看| av不卡在线播放| 亚洲精品在线美女| 看十八女毛片水多多多| 在线观看三级黄色| 青春草国产在线视频| 日韩精品免费视频一区二区三区| 国产色婷婷99| 国产女主播在线喷水免费视频网站| av福利片在线| 精品少妇内射三级| 在线观看www视频免费| 看十八女毛片水多多多| 欧美人与性动交α欧美精品济南到| 搡老乐熟女国产| 秋霞在线观看毛片| 国产精品久久久人人做人人爽| 日韩人妻精品一区2区三区| 777米奇影视久久| av卡一久久| 国产一级毛片在线| 欧美人与善性xxx| 国产极品天堂在线| 高清av免费在线| 王馨瑶露胸无遮挡在线观看| 婷婷色av中文字幕| 80岁老熟妇乱子伦牲交| 亚洲欧美一区二区三区黑人| av片东京热男人的天堂| 高清欧美精品videossex| 国产成人系列免费观看| 亚洲一码二码三码区别大吗| 欧美少妇被猛烈插入视频| 亚洲欧美激情在线| 最近手机中文字幕大全| 在现免费观看毛片| 亚洲欧美色中文字幕在线| 最近中文字幕2019免费版| 最新在线观看一区二区三区 | 日本午夜av视频| 免费久久久久久久精品成人欧美视频| 九色亚洲精品在线播放| 成年动漫av网址| 亚洲一码二码三码区别大吗| 自拍欧美九色日韩亚洲蝌蚪91| 9热在线视频观看99| 国产av一区二区精品久久| 国产精品 国内视频| 国产成人啪精品午夜网站| 国产精品99久久99久久久不卡 | 在线观看国产h片| 精品人妻一区二区三区麻豆| 日本午夜av视频| 侵犯人妻中文字幕一二三四区| 成年人免费黄色播放视频| 午夜久久久在线观看| 十八禁高潮呻吟视频| 国产精品一国产av| 少妇 在线观看| 免费女性裸体啪啪无遮挡网站| 黄色毛片三级朝国网站| 亚洲精华国产精华液的使用体验| 久久久国产一区二区| 伊人久久国产一区二区| 视频区图区小说| 麻豆精品久久久久久蜜桃| 男女免费视频国产| 一本久久精品| 成人亚洲欧美一区二区av| 免费观看性生交大片5| 另类亚洲欧美激情| 亚洲欧洲日产国产| 精品国产一区二区三区四区第35| 亚洲一级一片aⅴ在线观看| 免费黄网站久久成人精品| 久久性视频一级片| 国产日韩欧美视频二区| 最近最新中文字幕大全免费视频 | 久久精品国产a三级三级三级| 如日韩欧美国产精品一区二区三区| 9热在线视频观看99| 90打野战视频偷拍视频| 一级片'在线观看视频| 亚洲自偷自拍图片 自拍| 国产成人av激情在线播放| 亚洲欧美日韩另类电影网站| 日韩视频在线欧美| 久久免费观看电影| 大片电影免费在线观看免费| 最近的中文字幕免费完整| 亚洲国产精品999| 成人手机av| 高清不卡的av网站| videos熟女内射| 亚洲激情五月婷婷啪啪| 又大又黄又爽视频免费| 一本大道久久a久久精品| 亚洲国产看品久久| 欧美日韩成人在线一区二区| 最近中文字幕高清免费大全6| 日韩制服丝袜自拍偷拍| 99国产综合亚洲精品| 国产视频首页在线观看| 精品亚洲乱码少妇综合久久| 日韩av不卡免费在线播放| 国产成人免费无遮挡视频| 日韩精品免费视频一区二区三区| 中文天堂在线官网| 永久免费av网站大全| 日韩一本色道免费dvd| 国产免费视频播放在线视频| 男女午夜视频在线观看| 亚洲国产欧美一区二区综合| 亚洲成人手机| 精品人妻在线不人妻| 日韩一本色道免费dvd| 久久午夜综合久久蜜桃| 如何舔出高潮| 制服丝袜香蕉在线| 亚洲精品美女久久久久99蜜臀 | 欧美老熟妇乱子伦牲交| bbb黄色大片| 国产老妇伦熟女老妇高清| 另类亚洲欧美激情| 一区二区av电影网| 又黄又粗又硬又大视频| 亚洲精品美女久久av网站| 亚洲精品乱久久久久久| 老汉色∧v一级毛片| 亚洲七黄色美女视频| 国产片内射在线| 亚洲精品国产区一区二| 国产熟女欧美一区二区| 欧美中文综合在线视频| 日韩中文字幕欧美一区二区 | 免费在线观看完整版高清| 国产 一区精品| 久久天躁狠狠躁夜夜2o2o | 亚洲自偷自拍图片 自拍| 精品视频人人做人人爽| 日韩不卡一区二区三区视频在线| 人人澡人人妻人| 黄色 视频免费看| 一级毛片黄色毛片免费观看视频| 一级毛片电影观看| 久久久久久人妻| 最近中文字幕2019免费版| 久久久久精品国产欧美久久久 | 1024视频免费在线观看| 黄色毛片三级朝国网站| 日韩精品免费视频一区二区三区| 亚洲天堂av无毛| 午夜福利影视在线免费观看| 欧美老熟妇乱子伦牲交| 99久久精品国产亚洲精品| 韩国av在线不卡| 国产av码专区亚洲av| 亚洲av电影在线观看一区二区三区| 亚洲成人免费av在线播放| 最近最新中文字幕免费大全7| 精品人妻一区二区三区麻豆| 国产国语露脸激情在线看| 9热在线视频观看99| 国产一卡二卡三卡精品 | 极品少妇高潮喷水抽搐| 亚洲欧美精品自产自拍| 成人漫画全彩无遮挡| 久久精品国产a三级三级三级| 视频在线观看一区二区三区| 看免费av毛片| 国产精品香港三级国产av潘金莲 | 少妇人妻久久综合中文| 国产精品免费视频内射| 成人亚洲精品一区在线观看| 亚洲欧洲日产国产| 男女午夜视频在线观看| 两性夫妻黄色片| 老司机影院成人| 黄色一级大片看看| 国产精品 欧美亚洲| 99久久人妻综合| 久久性视频一级片| 女性被躁到高潮视频| 伦理电影大哥的女人| 中文字幕最新亚洲高清| 亚洲国产毛片av蜜桃av| 午夜av观看不卡| 热99国产精品久久久久久7| 亚洲av日韩精品久久久久久密 | 国产成人免费观看mmmm| 精品久久蜜臀av无| 各种免费的搞黄视频| 亚洲av日韩精品久久久久久密 | 亚洲精品美女久久av网站| 日本一区二区免费在线视频| 亚洲精品日本国产第一区| 欧美日韩一区二区视频在线观看视频在线| 亚洲 欧美一区二区三区| 人人妻人人爽人人添夜夜欢视频| 观看av在线不卡| 国产精品亚洲av一区麻豆 | 青春草国产在线视频| 日韩一区二区视频免费看| 自拍欧美九色日韩亚洲蝌蚪91| 国产亚洲av片在线观看秒播厂| 极品人妻少妇av视频| 国产欧美日韩一区二区三区在线| 亚洲七黄色美女视频| 大片免费播放器 马上看| 亚洲在久久综合| 久久久精品94久久精品| 在线观看一区二区三区激情| 久久精品久久久久久噜噜老黄| 麻豆精品久久久久久蜜桃| 精品福利永久在线观看| 9191精品国产免费久久| 看免费成人av毛片| 久久午夜综合久久蜜桃| 国产淫语在线视频| 国产欧美日韩综合在线一区二区| 亚洲欧美日韩另类电影网站| 亚洲精品aⅴ在线观看| 日本猛色少妇xxxxx猛交久久| 国产激情久久老熟女| 97在线人人人人妻| 免费在线观看黄色视频的| 一边摸一边做爽爽视频免费| 久久精品国产综合久久久| 久久久久精品人妻al黑| av一本久久久久| 精品国产一区二区三区四区第35| 中文精品一卡2卡3卡4更新| 在线观看免费视频网站a站| 亚洲国产精品一区三区| 国产免费现黄频在线看| 在线 av 中文字幕| 久久久亚洲精品成人影院| 99久国产av精品国产电影| 男女床上黄色一级片免费看| www.av在线官网国产| 91老司机精品| 国产一卡二卡三卡精品 | xxxhd国产人妻xxx| 免费在线观看完整版高清| 黑人猛操日本美女一级片| 精品亚洲乱码少妇综合久久| 日韩制服骚丝袜av| 亚洲精品国产区一区二| 交换朋友夫妻互换小说| 久久午夜综合久久蜜桃| 考比视频在线观看| 又粗又硬又长又爽又黄的视频| 涩涩av久久男人的天堂| 操出白浆在线播放| 美女视频免费永久观看网站| 国产精品国产三级专区第一集| 女人久久www免费人成看片| 亚洲情色 制服丝袜| 亚洲成人国产一区在线观看 | 久久人人爽av亚洲精品天堂| 热99久久久久精品小说推荐| 熟女少妇亚洲综合色aaa.| 午夜福利视频在线观看免费| 亚洲av成人精品一二三区| 国产精品一二三区在线看| 十八禁人妻一区二区| 在线免费观看不下载黄p国产| 久久99一区二区三区| 人妻人人澡人人爽人人| 国产av精品麻豆| 91老司机精品| 亚洲国产精品999| 成人漫画全彩无遮挡| 国产精品二区激情视频| 我的亚洲天堂| 日韩,欧美,国产一区二区三区| 亚洲国产欧美在线一区| 日韩一区二区视频免费看| 国产亚洲最大av| 考比视频在线观看| 午夜久久久在线观看| av在线app专区| 视频区图区小说| 亚洲专区中文字幕在线 | 久久久久久免费高清国产稀缺| 丝袜人妻中文字幕| 韩国av在线不卡| 汤姆久久久久久久影院中文字幕| 亚洲欧美日韩另类电影网站| 男女边摸边吃奶| 韩国av在线不卡| 免费人妻精品一区二区三区视频| 国产乱来视频区| 黄片无遮挡物在线观看| 黄网站色视频无遮挡免费观看| 国产精品偷伦视频观看了| 久久精品久久久久久久性| 香蕉国产在线看| 欧美成人精品欧美一级黄| 别揉我奶头~嗯~啊~动态视频 | 永久免费av网站大全| 18禁裸乳无遮挡动漫免费视频| 亚洲国产精品国产精品| 哪个播放器可以免费观看大片| 欧美亚洲日本最大视频资源| 久久精品亚洲av国产电影网| 久久久精品区二区三区| 国产精品国产三级专区第一集| 亚洲欧美日韩另类电影网站| 久久久久久久久久久久大奶| 久久久久网色| 久久久亚洲精品成人影院| 一边摸一边做爽爽视频免费| 国产精品偷伦视频观看了| 超碰97精品在线观看| 99香蕉大伊视频| 晚上一个人看的免费电影| 欧美日韩福利视频一区二区| 狠狠精品人妻久久久久久综合| 又大又黄又爽视频免费| 国产乱人偷精品视频| 18禁裸乳无遮挡动漫免费视频| 男人添女人高潮全过程视频| 亚洲精品国产色婷婷电影| 欧美黑人欧美精品刺激| 日本爱情动作片www.在线观看| 亚洲中文av在线| 国产免费一区二区三区四区乱码| 男女之事视频高清在线观看 | 精品久久蜜臀av无| 欧美黑人精品巨大| 男女国产视频网站| 亚洲精品美女久久久久99蜜臀 | 男女床上黄色一级片免费看| 亚洲精品久久午夜乱码| 90打野战视频偷拍视频| 欧美乱码精品一区二区三区| 夫妻性生交免费视频一级片| 欧美国产精品va在线观看不卡| 国产伦人伦偷精品视频| 欧美日韩av久久| 制服诱惑二区| 嫩草影视91久久| 最新在线观看一区二区三区 | 国产在视频线精品| 欧美日韩一区二区视频在线观看视频在线| 又黄又粗又硬又大视频| 一级毛片黄色毛片免费观看视频| 欧美人与性动交α欧美软件| 纵有疾风起免费观看全集完整版| 97精品久久久久久久久久精品| 青草久久国产| 999精品在线视频| 精品一区二区三区四区五区乱码 | 97精品久久久久久久久久精品| 国产黄频视频在线观看| 午夜精品国产一区二区电影| 日韩制服骚丝袜av| 欧美97在线视频| 中国国产av一级| 大香蕉久久网| 中文字幕精品免费在线观看视频| 天天操日日干夜夜撸| 国产精品一二三区在线看| 日本一区二区免费在线视频| 日本91视频免费播放| 2018国产大陆天天弄谢| 热99久久久久精品小说推荐| 亚洲精品久久午夜乱码| 午夜福利在线免费观看网站| www日本在线高清视频| 最新的欧美精品一区二区| 国产1区2区3区精品| 亚洲av综合色区一区| 久久精品亚洲av国产电影网| 久久久久久人妻| 爱豆传媒免费全集在线观看| 亚洲视频免费观看视频| 日韩成人av中文字幕在线观看|