• <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)軸承故障診斷
    国产亚洲精品久久久久久毛片 | 亚洲情色 制服丝袜| 国产亚洲精品一区二区www | 免费观看a级毛片全部| 大码成人一级视频| 看黄色毛片网站| 久久久久久免费高清国产稀缺| 亚洲国产毛片av蜜桃av| 天堂中文最新版在线下载| 亚洲人成电影免费在线| 日日夜夜操网爽| 亚洲欧美日韩另类电影网站| av视频免费观看在线观看| 精品久久久久久久毛片微露脸| 免费看十八禁软件| 国产精品欧美亚洲77777| 一本大道久久a久久精品| 久久久久久久国产电影| 视频区欧美日本亚洲| 91国产中文字幕| 夫妻午夜视频| 侵犯人妻中文字幕一二三四区| 国产xxxxx性猛交| 精品国产一区二区三区久久久樱花| 久9热在线精品视频| 成人影院久久| 中文亚洲av片在线观看爽 | www.熟女人妻精品国产| 一级黄色大片毛片| 在线免费观看的www视频| 国产精品久久久久久人妻精品电影| 十分钟在线观看高清视频www| 少妇猛男粗大的猛烈进出视频| 婷婷成人精品国产| 亚洲熟妇熟女久久| 国产成人一区二区三区免费视频网站| 色播在线永久视频| 波多野结衣av一区二区av| 国产精品一区二区在线观看99| 午夜精品在线福利| 怎么达到女性高潮| 伊人久久大香线蕉亚洲五| 天堂动漫精品| 亚洲色图 男人天堂 中文字幕| 国产高清国产精品国产三级| 亚洲av片天天在线观看| 亚洲欧美精品综合一区二区三区| 美女视频免费永久观看网站| 纯流量卡能插随身wifi吗| 好男人电影高清在线观看| 韩国精品一区二区三区| 久久精品aⅴ一区二区三区四区| 久久精品国产亚洲av香蕉五月 | 国产亚洲欧美精品永久| 国产日韩一区二区三区精品不卡| 大型黄色视频在线免费观看| 精品国产一区二区三区四区第35| 黄色视频不卡| 亚洲精华国产精华精| 成人手机av| 国产极品粉嫩免费观看在线| 法律面前人人平等表现在哪些方面| 王馨瑶露胸无遮挡在线观看| 日韩欧美国产一区二区入口| 国产成人精品在线电影| 国产区一区二久久| 免费在线观看视频国产中文字幕亚洲| 久久久精品免费免费高清| 一级片免费观看大全| 法律面前人人平等表现在哪些方面| 老司机在亚洲福利影院| 人妻一区二区av| 91国产中文字幕| 精品国产乱子伦一区二区三区| 丝袜在线中文字幕| 国产精品自产拍在线观看55亚洲 | 精品视频人人做人人爽| 大香蕉久久成人网| 午夜免费鲁丝| 新久久久久国产一级毛片| 日韩视频一区二区在线观看| 成熟少妇高潮喷水视频| 下体分泌物呈黄色| 亚洲成av片中文字幕在线观看| 免费观看a级毛片全部| 黄色毛片三级朝国网站| 十八禁高潮呻吟视频| 大陆偷拍与自拍| 首页视频小说图片口味搜索| 变态另类成人亚洲欧美熟女 | 亚洲成人免费av在线播放| 免费观看人在逋| 可以免费在线观看a视频的电影网站| 丝袜美足系列| 村上凉子中文字幕在线| 亚洲一区二区三区欧美精品| 免费一级毛片在线播放高清视频 | 午夜两性在线视频| 久久精品人人爽人人爽视色| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲午夜精品一区,二区,三区| 亚洲片人在线观看| 欧美日韩视频精品一区| 国产在视频线精品| 日韩精品免费视频一区二区三区| 成人手机av| 国产99久久九九免费精品| 日韩免费av在线播放| 老司机影院毛片| 最新在线观看一区二区三区| 女同久久另类99精品国产91| 99在线人妻在线中文字幕 | av欧美777| 午夜福利在线观看吧| 久久影院123| 悠悠久久av| 视频区欧美日本亚洲| 十分钟在线观看高清视频www| 岛国在线观看网站| 人妻一区二区av| 国产精品美女特级片免费视频播放器 | 免费高清在线观看日韩| 亚洲精品国产一区二区精华液| 一二三四在线观看免费中文在| 他把我摸到了高潮在线观看| 中文字幕人妻熟女乱码| 岛国在线观看网站| 欧美黑人欧美精品刺激| 女人精品久久久久毛片| 一边摸一边做爽爽视频免费| 国产精品自产拍在线观看55亚洲 | 亚洲成人国产一区在线观看| 日本撒尿小便嘘嘘汇集6| 国产精品乱码一区二三区的特点 | tube8黄色片| 欧美精品人与动牲交sv欧美| 国产成人免费无遮挡视频| 极品少妇高潮喷水抽搐| 欧美在线一区亚洲| 日日爽夜夜爽网站| 精品视频人人做人人爽| 女警被强在线播放| 97人妻天天添夜夜摸| 99精品欧美一区二区三区四区| 在线观看免费午夜福利视频| 激情在线观看视频在线高清 | 91九色精品人成在线观看| 亚洲第一av免费看| 男人操女人黄网站| 亚洲国产精品合色在线| 大香蕉久久网| 久久草成人影院| 变态另类成人亚洲欧美熟女 | 丰满的人妻完整版| x7x7x7水蜜桃| 变态另类成人亚洲欧美熟女 | 国产精品免费一区二区三区在线 | 国产亚洲精品久久久久久毛片 | 美女国产高潮福利片在线看| 日本黄色日本黄色录像| 极品人妻少妇av视频| 美女国产高潮福利片在线看| 久久久国产精品麻豆| 777米奇影视久久| 女同久久另类99精品国产91| 日韩欧美免费精品| 法律面前人人平等表现在哪些方面| 国产无遮挡羞羞视频在线观看| 视频区欧美日本亚洲| 91麻豆av在线| 国产高清视频在线播放一区| 亚洲av熟女| 丰满迷人的少妇在线观看| 丝瓜视频免费看黄片| cao死你这个sao货| 一区二区三区国产精品乱码| 欧美老熟妇乱子伦牲交| 国产单亲对白刺激| 每晚都被弄得嗷嗷叫到高潮| 99国产精品一区二区蜜桃av | 黑丝袜美女国产一区| 国产男女内射视频| 久久精品国产a三级三级三级| 天天操日日干夜夜撸| 精品少妇久久久久久888优播| videos熟女内射| 最新在线观看一区二区三区| 亚洲三区欧美一区| 亚洲性夜色夜夜综合| www.精华液| 国产精品综合久久久久久久免费 | 男女午夜视频在线观看| 亚洲五月婷婷丁香| av不卡在线播放| 人妻一区二区av| aaaaa片日本免费| 久久国产精品影院| 久久国产精品大桥未久av| www.精华液| 精品电影一区二区在线| 亚洲人成电影观看| 9191精品国产免费久久| 国产欧美日韩一区二区三| 久99久视频精品免费| 99国产精品一区二区三区| 亚洲男人天堂网一区| 操美女的视频在线观看| 国产成人欧美| 久久久久国内视频| 欧美精品啪啪一区二区三区| 亚洲人成77777在线视频| 日韩欧美在线二视频 | 色综合欧美亚洲国产小说| 999精品在线视频| 中文字幕另类日韩欧美亚洲嫩草| av国产精品久久久久影院| 国产一区有黄有色的免费视频| 久久精品亚洲av国产电影网| 国产精品国产高清国产av | 欧美人与性动交α欧美软件| 一二三四社区在线视频社区8| 中文字幕色久视频| 三上悠亚av全集在线观看| 日本a在线网址| 黄色a级毛片大全视频| 精品亚洲成国产av| 操出白浆在线播放| 亚洲欧美激情在线| 精品福利永久在线观看| 成人永久免费在线观看视频| 亚洲精品粉嫩美女一区| 正在播放国产对白刺激| 十八禁高潮呻吟视频| 大香蕉久久成人网| 国产精品99久久99久久久不卡| 午夜福利免费观看在线| 脱女人内裤的视频| 久久亚洲精品不卡| 亚洲成国产人片在线观看| 精品一区二区三区视频在线观看免费 | 在线观看免费视频网站a站| 成人黄色视频免费在线看| 欧美日韩瑟瑟在线播放| 精品国产乱码久久久久久男人| 50天的宝宝边吃奶边哭怎么回事| 一区福利在线观看| 首页视频小说图片口味搜索| 久久精品国产a三级三级三级| 十八禁人妻一区二区| 精品国产国语对白av| 亚洲第一青青草原| 免费av中文字幕在线| 韩国精品一区二区三区| 午夜精品久久久久久毛片777| 欧美精品亚洲一区二区| 国产伦人伦偷精品视频| 韩国av一区二区三区四区| 亚洲成av片中文字幕在线观看| 国产男女内射视频| 巨乳人妻的诱惑在线观看| 国产欧美日韩综合在线一区二区| 超碰97精品在线观看| 国产成人精品无人区| 一级毛片高清免费大全| 亚洲熟妇熟女久久| 黄色怎么调成土黄色| 亚洲精品久久成人aⅴ小说| 免费不卡黄色视频| 成人国产一区最新在线观看| 欧美色视频一区免费| 午夜激情av网站| 99re6热这里在线精品视频| 中文字幕人妻熟女乱码| 丰满迷人的少妇在线观看| 午夜福利在线观看吧| 99久久精品国产亚洲精品| 精品久久久久久久久久免费视频 | svipshipincom国产片| 精品国产超薄肉色丝袜足j| 国产成人精品在线电影| 美女扒开内裤让男人捅视频| 精品国产亚洲在线| 国产午夜精品久久久久久| 99精品久久久久人妻精品| 欧美性长视频在线观看| 久久精品国产清高在天天线| 青草久久国产| 亚洲精品av麻豆狂野| 精品电影一区二区在线| 俄罗斯特黄特色一大片| 亚洲成国产人片在线观看| 建设人人有责人人尽责人人享有的| 99re在线观看精品视频| 青草久久国产| av天堂久久9| 女人被躁到高潮嗷嗷叫费观| 91精品三级在线观看| 中文字幕人妻熟女乱码| 超碰成人久久| 国产一卡二卡三卡精品| 久久久久国内视频| 91麻豆av在线| 日本黄色视频三级网站网址 | 99热国产这里只有精品6| www.精华液| 每晚都被弄得嗷嗷叫到高潮| 精品久久久久久久毛片微露脸| 国产精品美女特级片免费视频播放器 | 亚洲在线自拍视频| 悠悠久久av| 免费在线观看完整版高清| 夜夜夜夜夜久久久久| 老司机深夜福利视频在线观看| 美女午夜性视频免费| 一夜夜www| 国产1区2区3区精品| 精品视频人人做人人爽| 国产精品成人在线| 亚洲一区二区三区不卡视频| 午夜免费观看网址| 这个男人来自地球电影免费观看| 国产激情欧美一区二区| 搡老岳熟女国产| 国内毛片毛片毛片毛片毛片| 人妻久久中文字幕网| 亚洲熟女毛片儿| 99国产精品免费福利视频| 成人黄色视频免费在线看| 精品国产国语对白av| 国产真人三级小视频在线观看| 飞空精品影院首页| 亚洲第一青青草原| 日本撒尿小便嘘嘘汇集6| 麻豆乱淫一区二区| 免费高清在线观看日韩| 多毛熟女@视频| 12—13女人毛片做爰片一| 亚洲熟妇中文字幕五十中出 | 女警被强在线播放| 岛国在线观看网站| www.999成人在线观看| 少妇的丰满在线观看| 亚洲欧美日韩另类电影网站| 80岁老熟妇乱子伦牲交| 操出白浆在线播放| 欧美人与性动交α欧美软件| 欧美黑人精品巨大| 黄片大片在线免费观看| 欧美最黄视频在线播放免费 | 女人爽到高潮嗷嗷叫在线视频| 日韩 欧美 亚洲 中文字幕| 男人的好看免费观看在线视频 | 国产精华一区二区三区| 丰满迷人的少妇在线观看| 午夜精品国产一区二区电影| 亚洲av电影在线进入| 久久久国产成人精品二区 | a级毛片黄视频| 日本黄色日本黄色录像| 777米奇影视久久| 精品国产亚洲在线| 大香蕉久久网| 大片电影免费在线观看免费| 国产无遮挡羞羞视频在线观看| 国产男女超爽视频在线观看| 午夜成年电影在线免费观看| 不卡一级毛片| 国产伦人伦偷精品视频| 自拍欧美九色日韩亚洲蝌蚪91| 精品无人区乱码1区二区| 高清欧美精品videossex| 日韩欧美一区视频在线观看| 久久九九热精品免费| 无限看片的www在线观看| 色尼玛亚洲综合影院| 三上悠亚av全集在线观看| 男女下面插进去视频免费观看| 精品少妇久久久久久888优播| 妹子高潮喷水视频| 国产亚洲精品久久久久久毛片 | 狂野欧美激情性xxxx| 一级黄色大片毛片| 久久中文字幕人妻熟女| 精品国内亚洲2022精品成人 | 久久精品亚洲av国产电影网| 黄色丝袜av网址大全| 激情视频va一区二区三区| 久久中文字幕一级| 日日爽夜夜爽网站| 日韩精品免费视频一区二区三区| 欧美人与性动交α欧美精品济南到| 久热爱精品视频在线9| 一夜夜www| 99国产精品99久久久久| 日本精品一区二区三区蜜桃| 人妻久久中文字幕网| 日本五十路高清| 中文字幕人妻熟女乱码| 一边摸一边抽搐一进一出视频| 国产成人精品在线电影| 丰满迷人的少妇在线观看| 国产不卡一卡二| 日本撒尿小便嘘嘘汇集6| 久久香蕉精品热| 久久中文字幕一级| 欧美久久黑人一区二区| 日韩成人在线观看一区二区三区| 日韩有码中文字幕| 超碰97精品在线观看| 在线观看午夜福利视频| 99国产精品99久久久久| 国产av一区二区精品久久| 99热国产这里只有精品6| 成年版毛片免费区| 亚洲人成77777在线视频| 丰满人妻熟妇乱又伦精品不卡| 新久久久久国产一级毛片| 亚洲一卡2卡3卡4卡5卡精品中文| 不卡av一区二区三区| 人人妻人人澡人人爽人人夜夜| 亚洲色图综合在线观看| 在线永久观看黄色视频| 最新的欧美精品一区二区| 亚洲色图av天堂| 亚洲欧美一区二区三区黑人| 免费一级毛片在线播放高清视频 | 免费观看精品视频网站| 国产xxxxx性猛交| 久久精品国产99精品国产亚洲性色 | 亚洲情色 制服丝袜| 大型黄色视频在线免费观看| 丝袜美腿诱惑在线| 伦理电影免费视频| tube8黄色片| e午夜精品久久久久久久| 国产一区二区三区视频了| 亚洲aⅴ乱码一区二区在线播放 | 国产97色在线日韩免费| 身体一侧抽搐| 日韩大码丰满熟妇| 一二三四在线观看免费中文在| 亚洲中文av在线| 日韩有码中文字幕| 精品国产超薄肉色丝袜足j| cao死你这个sao货| 免费日韩欧美在线观看| 亚洲人成77777在线视频| 国产精品电影一区二区三区 | www.熟女人妻精品国产| 国产野战对白在线观看| 夜夜夜夜夜久久久久| 黄网站色视频无遮挡免费观看| 50天的宝宝边吃奶边哭怎么回事| 精品久久久久久,| 大香蕉久久成人网| 曰老女人黄片| 丝袜美腿诱惑在线| 日日摸夜夜添夜夜添小说| 80岁老熟妇乱子伦牲交| 免费久久久久久久精品成人欧美视频| 久久精品国产a三级三级三级| 一区在线观看完整版| 亚洲欧美激情在线| 欧美日韩成人在线一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久精品国产欧美久久久| 国产又爽黄色视频| 精品久久久久久电影网| 久久久精品区二区三区| 亚洲精品中文字幕一二三四区| 精品国内亚洲2022精品成人 | 18禁裸乳无遮挡免费网站照片 | 色老头精品视频在线观看| 欧美日韩乱码在线| 99国产精品99久久久久| 成年人免费黄色播放视频| 看免费av毛片| av有码第一页| 精品卡一卡二卡四卡免费| 国产精品一区二区在线观看99| 啦啦啦在线免费观看视频4| 欧美日本中文国产一区发布| 成人18禁高潮啪啪吃奶动态图| 国产高清激情床上av| 可以免费在线观看a视频的电影网站| 成人国产一区最新在线观看| 黄色片一级片一级黄色片| 久久人人97超碰香蕉20202| 成人黄色视频免费在线看| 久久久久久人人人人人| 男人操女人黄网站| ponron亚洲| 中亚洲国语对白在线视频| 久久精品亚洲熟妇少妇任你| 免费观看人在逋| 香蕉久久夜色| 少妇裸体淫交视频免费看高清 | 色94色欧美一区二区| 精品国产乱子伦一区二区三区| 亚洲熟妇熟女久久| 99精品在免费线老司机午夜| 久久久久久久午夜电影 | 欧美日韩亚洲高清精品| 母亲3免费完整高清在线观看| 欧美激情久久久久久爽电影 | 亚洲专区国产一区二区| 亚洲欧美一区二区三区黑人| 国产亚洲精品一区二区www | 男男h啪啪无遮挡| 国产一区二区三区在线臀色熟女 | 亚洲va日本ⅴa欧美va伊人久久| 国产精品自产拍在线观看55亚洲 | 久久香蕉国产精品| 久久人人爽av亚洲精品天堂| 他把我摸到了高潮在线观看| 多毛熟女@视频| 在线观看免费视频网站a站| 久久午夜综合久久蜜桃| 人人妻人人爽人人添夜夜欢视频| 久久ye,这里只有精品| 另类亚洲欧美激情| 伦理电影免费视频| 少妇被粗大的猛进出69影院| 中文字幕色久视频| 亚洲国产精品sss在线观看 | 久久天躁狠狠躁夜夜2o2o| 亚洲专区国产一区二区| 亚洲一区二区三区欧美精品| 日韩中文字幕欧美一区二区| 国产色视频综合| 又紧又爽又黄一区二区| 国产又爽黄色视频| 婷婷丁香在线五月| 亚洲精品av麻豆狂野| 国产精品自产拍在线观看55亚洲 | 91麻豆精品激情在线观看国产 | 一本综合久久免费| 久久精品国产亚洲av高清一级| 久久精品亚洲熟妇少妇任你| 中国美女看黄片| 一进一出抽搐gif免费好疼 | 很黄的视频免费| 天天影视国产精品| 中出人妻视频一区二区| 精品一区二区三区视频在线观看免费 | 欧美日韩亚洲高清精品| 亚洲一卡2卡3卡4卡5卡精品中文| 成人av一区二区三区在线看| aaaaa片日本免费| 精品一区二区三卡| 男女高潮啪啪啪动态图| 亚洲男人天堂网一区| 成人影院久久| 十八禁网站免费在线| 在线国产一区二区在线| 一夜夜www| 国产亚洲精品久久久久久毛片 | 91字幕亚洲| 首页视频小说图片口味搜索| 精品一品国产午夜福利视频| 精品第一国产精品| 久久精品人人爽人人爽视色| 中文字幕精品免费在线观看视频| 99国产精品一区二区蜜桃av | 香蕉国产在线看| 久久精品国产a三级三级三级| 91麻豆av在线| a级片在线免费高清观看视频| 精品国内亚洲2022精品成人 | 免费久久久久久久精品成人欧美视频| 国产欧美日韩一区二区三| 亚洲精品久久午夜乱码| 国产av精品麻豆| 美女国产高潮福利片在线看| 久久九九热精品免费| 久久精品国产a三级三级三级| 99久久综合精品五月天人人| 亚洲成人国产一区在线观看| 99精品在免费线老司机午夜| 国产精品av久久久久免费| 一边摸一边抽搐一进一出视频| 男人的好看免费观看在线视频 | 下体分泌物呈黄色| 欧美乱妇无乱码| 久久久久久人人人人人| 欧美丝袜亚洲另类 | 亚洲七黄色美女视频| 久热爱精品视频在线9| 亚洲一卡2卡3卡4卡5卡精品中文| 日本五十路高清| 国产精品欧美亚洲77777| 深夜精品福利| 亚洲成人手机| 国产成人精品久久二区二区91| 亚洲精品在线美女| 无人区码免费观看不卡| 国产成人系列免费观看| 成人影院久久| 久久久国产精品麻豆| 亚洲午夜精品一区,二区,三区| 在线观看免费日韩欧美大片| 女人被躁到高潮嗷嗷叫费观| 成人国语在线视频| 人妻丰满熟妇av一区二区三区 | 精品国产一区二区久久| 精品一区二区三卡| 9色porny在线观看| 成人亚洲精品一区在线观看| 欧美国产精品va在线观看不卡| 亚洲一区二区三区欧美精品| 国产精品综合久久久久久久免费 | 日本五十路高清| 国产精品久久电影中文字幕 | 免费在线观看黄色视频的| 亚洲 国产 在线|