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

    基于GADF-TL-ResNeXt的滾動(dòng)軸承故障診斷方法

    2023-11-03 11:34:38侯東曉周子安程榮財(cái)
    計(jì)量學(xué)報(bào) 2023年10期
    關(guān)鍵詞:殘差分組故障診斷

    侯東曉, 周子安, 程榮財(cái), 閻 爽

    (1.東北大學(xué)秦皇島分校 控制工程學(xué)院 河北 秦皇島,066004;2.京能秦皇島熱電有限公司 河北 秦皇島,066004)

    1 引 言

    在現(xiàn)代化設(shè)備制造業(yè)中,機(jī)械制造正朝著自動(dòng)化、智能化的方向不斷發(fā)展。其中,旋轉(zhuǎn)機(jī)械結(jié)構(gòu)更加精密、功能更為復(fù)雜,在各個(gè)生產(chǎn)領(lǐng)域發(fā)揮著重要作用[1]。滾動(dòng)軸承作為旋轉(zhuǎn)機(jī)械的代表零件之一,被用作支撐回轉(zhuǎn)件、減小摩擦以及保證回轉(zhuǎn)精度等,對機(jī)械設(shè)備的壽命與穩(wěn)定起著決定性的影響。在實(shí)際工況下,軸承往往需要承受變速、變載等嚴(yán)苛的工作條件,導(dǎo)致故障也時(shí)常發(fā)生。因此,需要及時(shí)、高效地查明軸承的故障原因。

    當(dāng)前傳統(tǒng)的軸承故障診斷方法有2種,一是以信號(hào)分析為主的診斷方法。Ye X等[2]將改進(jìn)的經(jīng)驗(yàn)?zāi)B(tài)分解(EMD)方法應(yīng)用到軸承故障診斷中,對振動(dòng)信號(hào)的微弱噪聲和突發(fā)脈沖具有良好的識(shí)別效果。Hou J等[3]提出將基于聚類算法的集合經(jīng)驗(yàn)?zāi)B(tài)分解(EEMD)應(yīng)用到軸承故障診斷中,結(jié)果表明所提的方法能夠獲得較高的聚類指標(biāo)。此方法的局限性是過分依賴人工參與,并且過程非常的耗費(fèi)時(shí)力。另一種則是以數(shù)據(jù)驅(qū)動(dòng)為主的機(jī)器學(xué)習(xí)方法。如Chen F F等[4]將預(yù)處理的振動(dòng)信號(hào)直接輸入到改進(jìn)的支持向量機(jī)(SVM)中,用于軸承早期故障的識(shí)別,并取得了良好的效果。在文獻(xiàn)[5]中,研究人員先將軸承一維振動(dòng)信號(hào)進(jìn)行預(yù)處理,再輸入到一維卷積神經(jīng)網(wǎng)絡(luò)(1D-CNN)中完成故障的識(shí)別與分類。此方法能夠自動(dòng)提取特征,且準(zhǔn)確率和魯棒性也大大提高,因此近些年被廣泛應(yīng)用于故障診斷領(lǐng)域。CNN因其有強(qiáng)大的圖像特征提取能力也受到了許多學(xué)者的青睞。Li G等[6]通過ST變換將振動(dòng)信號(hào)直接變?yōu)槎S矩陣并輸入到CNN中進(jìn)行分類。Yang J等[7]通過生成對抗網(wǎng)絡(luò)來擴(kuò)充小樣本數(shù)據(jù)集,并轉(zhuǎn)換為灰度圖像輸入到SE-CNN自動(dòng)分類識(shí)別。Han T等[8]向二維圖像中添加紅色,將時(shí)域彩色特征圖輸入到CNN中以增強(qiáng)CNN的識(shí)別能力。但此類方法仍存在以下幾個(gè)問題:1) 直接以原始信號(hào)輸入模型中可能會(huì)導(dǎo)致部分特征提取不完整;2) 是隨著網(wǎng)絡(luò)模型的深度的加深,模型存在退化問題;3) 是模型需要成百上千次的訓(xùn)練才能達(dá)到一個(gè)好的效果。

    基于以上問題,本文提出了格拉姆角場(Gramian angular field,GAF)、遷移學(xué)習(xí)(transfer learning, TL)、分組殘差網(wǎng)絡(luò)(ResNeXt)相結(jié)合的數(shù)據(jù)驅(qū)動(dòng)故障診斷模型(GADF-TL-ResNeXt)。首先通過GAF將一維時(shí)序信號(hào)轉(zhuǎn)換成二維特征圖,再將這些特征圖輸入到與遷移學(xué)習(xí)結(jié)合的ResNeXt網(wǎng)絡(luò)中進(jìn)行自動(dòng)的識(shí)別、分類。為了驗(yàn)證該方法的有效性,使用凱斯西儲(chǔ)大學(xué)數(shù)據(jù)與其他方法進(jìn)行對比,最后使用軋機(jī)模擬實(shí)驗(yàn)平臺(tái)上收集的軸承故障數(shù)據(jù)進(jìn)行變工況下的軸承故障分類。

    2 數(shù)據(jù)到圖像的轉(zhuǎn)換方法

    格拉姆角場,是一種將時(shí)序信號(hào)轉(zhuǎn)換為二維圖像的編碼方式[9],它的原理是將時(shí)間序列信號(hào)在極坐標(biāo)系統(tǒng)內(nèi)進(jìn)行編碼,極坐標(biāo)系統(tǒng)保留了絕對的時(shí)間關(guān)系,GAF圖像編碼將一維時(shí)間序列1個(gè)點(diǎn)對應(yīng)的時(shí)間、幅度轉(zhuǎn)換到極坐標(biāo)系下的半徑、角度,這樣數(shù)據(jù)便從一維數(shù)據(jù)轉(zhuǎn)換到二維空間上,保留了對時(shí)間的依賴性。GAF在編碼時(shí)通過三角函數(shù)的和或差運(yùn)算生成GASF/GADF編碼圖,相同的數(shù)據(jù)點(diǎn)下對應(yīng)的圖像如圖1所示。

    由圖1可知,格拉姆角差場(GADF)相較于GASF在圖像色彩、交叉邊界以及細(xì)節(jié)刻畫方面效果均更好[10]。因此采用GADF進(jìn)行編碼計(jì)算(見圖2),假設(shè)有n個(gè)點(diǎn)的時(shí)間序列X={x1,x2,…,xn},GADF對此序列X轉(zhuǎn)化成二維空間,具體步驟為:

    圖2 GADF圖像編碼流程Fig.2 GADF images encoding process

    1) 對n個(gè)時(shí)間序列信號(hào)進(jìn)行歸一化處理,目的是將一維數(shù)據(jù)壓縮至區(qū)間[-1,1],處理式

    (1)

    2) 將縮放后的序列數(shù)據(jù)轉(zhuǎn)換到極坐標(biāo)系,即將數(shù)值看作夾角余弦值,時(shí)間戳看作半徑,從而在極坐標(biāo)中表示重新調(diào)整的時(shí)間序列X,轉(zhuǎn)換過程中,隨著時(shí)間的增加,相應(yīng)的數(shù)值在跨度圓上的不同角度點(diǎn)之間發(fā)生扭曲,從而得到極坐標(biāo)系下的曲線,具體可定義為

    (2)

    式中:φ表示角余弦的極坐標(biāo);ti表示時(shí)間戳;N表示正則化極坐標(biāo)系跨度的常數(shù)因子。

    3) 縮放后的時(shí)間序列轉(zhuǎn)換到極坐標(biāo)系統(tǒng)后,對于存儲(chǔ)時(shí)間信息的極坐標(biāo)系統(tǒng)φ,通過計(jì)算系統(tǒng)中每個(gè)極坐標(biāo)的三角函數(shù)差,可以很容易地進(jìn)行角度視角來識(shí)別出不同時(shí)間間隔內(nèi)的時(shí)間相關(guān)性??蓪⑵渚幋a到矩陣的幾何結(jié)構(gòu)中,定義為

    AGADF=[sin(φi-φj)]

    (3)

    (4)

    式中:I=[1,1,…,1],它是單位行向量;X′和X代表著不同行向量。

    4) 最后將AGADF矩陣表示為二維GADF圖像,并通過圖像建立極坐標(biāo)保持絕對時(shí)間關(guān)系。

    3 網(wǎng)絡(luò)結(jié)構(gòu)及原理

    3.1 殘差神經(jīng)網(wǎng)絡(luò)(ResNet)

    ResNet[11]屬于CNN的一種,其引入的殘差結(jié)構(gòu)可將網(wǎng)絡(luò)底層的特征通過捷徑直接傳遞到高層,通過引入快捷連接來擬合殘差映射函數(shù),只學(xué)習(xí)輸入和輸出的差別。ResNet由許多殘差塊堆疊組成,不僅加拓寬了網(wǎng)絡(luò)深度,還成功的解決了網(wǎng)絡(luò)退化問題,提高了模型訓(xùn)練的穩(wěn)定性,結(jié)構(gòu)如圖3所示。

    圖3 殘差結(jié)構(gòu)Fig.3 Residual structure

    圖3中x為殘差塊的輸入;F(x)為殘差映射函數(shù);H(x)=F(x,{Wi})+x為恒等映射函數(shù),{Wi}為輸入經(jīng)過殘差塊中第i個(gè)卷積層時(shí)獲得的權(quán)重。

    3.2 分組殘差結(jié)構(gòu)(GRS)

    分組殘差結(jié)構(gòu)是在原有的殘差模塊的基礎(chǔ)上引入了Inception[12]模塊思想,將殘差塊由原來的瓶頸結(jié)構(gòu)拆分成多分支結(jié)構(gòu),并且每條支路上都遵循相同的拓補(bǔ)結(jié)構(gòu),分別進(jìn)行相同的卷積操作,最后再進(jìn)行拼接求和匯總,具體可以定義為:

    (5)

    式中:X表示輸入數(shù)據(jù);Y表示輸出數(shù)據(jù);表示相同的拓補(bǔ)結(jié)構(gòu);G表示分支數(shù),經(jīng)實(shí)驗(yàn)驗(yàn)證,增加分組數(shù)量可以提高模型的精度,本文中G取32。

    當(dāng)殘差塊處于不同位置的時(shí)候,輸出數(shù)據(jù)的維度也不同,需要作出調(diào)整,其主要結(jié)構(gòu)分為2種:恒等映射圖4(a)和投影映射圖4(b)。當(dāng)輸入x與Conv3+BN層輸出的數(shù)據(jù)維度相同時(shí),二者可直接經(jīng)過捷徑(shortcut)相加。反之,則需要在捷徑分支上加一層卷積操作調(diào)整輸入數(shù)據(jù)x的維度以方便于Conv3+BN層輸出數(shù)據(jù)相加。

    圖4 2種分組殘差結(jié)構(gòu)Fig.4 Two group residual structures

    分組殘差結(jié)構(gòu)的優(yōu)點(diǎn)是不增加參數(shù)復(fù)雜度的條件下提高模型的精確度,對于未進(jìn)行分組操作的殘差塊,所提方法可以明顯降低殘差塊的參數(shù)量。假設(shè)輸入特征矩陣的通道數(shù)為C,卷積核大小為M×M,輸出特征矩陣的通道數(shù)為N,分為G組。

    未分組的殘差塊部分參數(shù)量Tw:

    Tw=C×M×M×N

    (6)

    分組的殘差塊部分參數(shù)量Tf:

    (7)

    由式(6)和(7)可知,相比較未分組的殘差塊而言,分組以后的殘差塊部分參數(shù)量要減少G倍,對參數(shù)處理方面的誤差會(huì)有一定的減小,圖5為本文ResNeXt網(wǎng)絡(luò)模型結(jié)構(gòu)。

    圖5 ResNeXt網(wǎng)絡(luò)框架結(jié)構(gòu)Fig.5 ResNeXt network framework structure

    3.3 遷移學(xué)習(xí)

    遷移學(xué)習(xí)給定源域Ds、相應(yīng)的原任務(wù)Ts及目標(biāo)域Dt和目標(biāo)任務(wù)Tt,其目的是使用Ds、Ts中的知識(shí)來改進(jìn)Dt中的目標(biāo)預(yù)測函數(shù)ft(·)。源域Ds、目標(biāo)域Dt可定義為:Ds={X,P(X)}、Dt={Y,ft(·)}。X代表特征空間,P(X)代表邊際概率分布。由此獲得Ds={(xs1,ys1),(xs2,ys2),…,(xsn,ysn)},其中數(shù)據(jù)樣本滿足xsi∈Xs,相應(yīng)的類標(biāo)簽滿足ysi∈Ys,Dt可以由數(shù)據(jù)樣本Xt和標(biāo)簽空間Yt表示。遷移學(xué)習(xí)還可以在相似或不同領(lǐng)域?qū)嵤┛焖賹W(xué)習(xí)[13]。由于故障診斷中相對的特征樣本量較少,ResNeXt模型很難在短時(shí)間內(nèi)達(dá)到一個(gè)完美訓(xùn)練效果,所以可加入遷移學(xué)習(xí)方法。此外,訓(xùn)練模型參數(shù)具備較強(qiáng)的特征遷移能力,對其他數(shù)據(jù)集進(jìn)行特征提取時(shí)可直接引入,有助強(qiáng)化模型性能,以學(xué)習(xí)特定的功能[14]。通過對初始模型在ImageNet數(shù)據(jù)集上先進(jìn)行預(yù)訓(xùn)練,可得到具有提煉抽象特征的模型權(quán)重參數(shù),即使加入樣本量較少的數(shù)據(jù)集,也可以有較高的識(shí)別準(zhǔn)確率。

    4 基于GADF與TL-ResNeXt故障診斷流程

    本文所用方法結(jié)合了GAF編碼可以保留時(shí)間序列的相關(guān)性以及加入遷移學(xué)習(xí)模塊的ResNeXt可以快速的識(shí)別時(shí)間序列特征能力的優(yōu)點(diǎn),提出了基于GADF-TL-ResNeXt的軸承故障診斷方法,具體流程如圖6所示。模型訓(xùn)練部分引入了ImageNet數(shù)據(jù)集上預(yù)訓(xùn)練好的模型權(quán)重,訓(xùn)練過程中改變了softmax層的分類節(jié)點(diǎn)個(gè)數(shù),以適用于軸承故障分類的幾種情況,再將訓(xùn)練過程中模型準(zhǔn)確率最大的權(quán)重文件引入到測試集模型中,最后進(jìn)行準(zhǔn)確率最高的故障分類,以達(dá)到最好的分類效果。

    圖6 基于GADF-TL-ResNeXt的故障診斷方法流程圖Fig.6 Fault diagnosis process based on GADF-TL-ResNeXt

    5 實(shí)驗(yàn)研究與分析

    5.1 CWRU軸承數(shù)據(jù)庫

    5.1.1 數(shù)據(jù)劃分

    為了評估所建模型的性能,使用來自美國凱斯西儲(chǔ)大學(xué)(Case Western Reserve University,CWRU)軸承數(shù)據(jù)庫的數(shù)據(jù)與其他方法進(jìn)行對比。以驅(qū)動(dòng)端深溝球軸承SKF6205為例,軸承的損傷是用電火花加工的單點(diǎn)損傷。分別對軸承的外圈(3、6、12點(diǎn)鐘方向)、內(nèi)圈、滾動(dòng)體設(shè)置0.17、0.35、0.53 mm 共3種故障類型損傷,在驅(qū)動(dòng)端軸承座上放置加速度信號(hào)傳感器,采樣頻率為12 kHz。實(shí)驗(yàn)中電機(jī)負(fù)載分為4類,選功率為0.75 kW,轉(zhuǎn)速為1772 r/min的情況。分別選用內(nèi)圈、外圈、滾動(dòng)體(6點(diǎn)鐘方向)3種類型的3種不同損傷直徑作為故障數(shù)據(jù)樣本。另設(shè)一種健康狀態(tài)的數(shù)據(jù)樣本,共10個(gè)類別,每種類別構(gòu)造400份樣本,按照9:1劃分為訓(xùn)練集和驗(yàn)證集。使用表1數(shù)據(jù),驗(yàn)證所提模型分類的有效性。

    表1 故障類型樣本Tab.1 Sample fault types

    5.1.2 樣本生成方法參數(shù)設(shè)置

    一維時(shí)序信號(hào)通過GADF編碼方法生成二維圖像,采用滑動(dòng)窗口在時(shí)序上進(jìn)行平移滑動(dòng)采集數(shù)據(jù)點(diǎn),同時(shí)生成對應(yīng)數(shù)據(jù)點(diǎn)下的二維圖像,為了確定滑動(dòng)窗口的大小,首先應(yīng)該計(jì)算軸承在旋轉(zhuǎn)一周內(nèi)采樣點(diǎn)的個(gè)數(shù),由公式(8)確定:

    (8)

    式中:L代表軸承一個(gè)周期內(nèi)采樣點(diǎn)的個(gè)數(shù);S代表軸承的轉(zhuǎn)速;Fq代表采樣頻率。

    為確保每個(gè)圖片編碼包含1個(gè)周期內(nèi)的采樣點(diǎn),經(jīng)過計(jì)算可知,軸承每轉(zhuǎn)采樣大約400個(gè)數(shù)據(jù)點(diǎn),所以滑動(dòng)窗口大小至少要為400,本文數(shù)據(jù)以最小的滑動(dòng)窗口采樣?;瑒?dòng)步長為滑動(dòng)窗口每次平移采樣的距離,當(dāng)步長小于窗口時(shí),樣本之間存在重疊,相同數(shù)據(jù)點(diǎn)下可以獲得更多的樣本,為保證數(shù)據(jù)利用率及信號(hào)的完整性,實(shí)驗(yàn)數(shù)據(jù)樣本處理方法為半額重疊采樣。為保證輸入圖片的分辨率適中,每一個(gè)滑動(dòng)窗口生成1張256×256的GADF二維編碼圖像,采樣所生成的故障圖如圖7所示。

    圖7 不同故障類型GADF編碼圖Fig.7 GADF coding diagram of different fault types

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

    為了驗(yàn)證模型的可行性與穩(wěn)定性,對比了其它幾類不同網(wǎng)絡(luò),分別是未加入遷移學(xué)習(xí)模塊的ResNeXt、遷移學(xué)習(xí)但未分組處理的殘差網(wǎng)絡(luò)TL-ResNet以及淺層網(wǎng)絡(luò)GoogLeNet、AlexNet。分別迭代140輪,由圖8可知,本模型相比于其他模型而言,訓(xùn)練第一輪就有87%的準(zhǔn)確率,在迭代到30輪時(shí)趨于穩(wěn)定,最高驗(yàn)證準(zhǔn)確率為99%,訓(xùn)練損失最低達(dá)到0.01。

    圖8 不同模型的準(zhǔn)確率和損失曲線Fig.8 Accuracy and loss curves for different models

    以上數(shù)據(jù)說明本模型相較于其他模型,魯棒性比較強(qiáng),并且趨于完全收斂,而其他未進(jìn)行遷移學(xué)習(xí)處理以及淺層網(wǎng)絡(luò)模型,很難在很少的輪數(shù)達(dá)到一個(gè)理想的效果,訓(xùn)練損失及穩(wěn)定性要遠(yuǎn)遠(yuǎn)低于本模型。

    5.1.4 模型評估指標(biāo)

    評估指標(biāo)是判斷診斷算法性能的標(biāo)準(zhǔn),在數(shù)據(jù)分析中具有重要意義。在深度學(xué)習(xí)故障診斷領(lǐng)域,準(zhǔn)確率Ac、精確率Pr、召回率Re、F1分?jǐn)?shù)都是評判模型性能的標(biāo)準(zhǔn),表達(dá)式如式(9)~式(12),通過實(shí)驗(yàn)驗(yàn)證,本模型對比了其他幾種流行網(wǎng)絡(luò)的這幾個(gè)參數(shù)值,如表2所示。

    表2 不同網(wǎng)絡(luò)模型各個(gè)參數(shù)Tab.2 Parameters of different network models (%)

    (9)

    (10)

    (11)

    (12)

    式中(假設(shè)分類器分類預(yù)測x、y共2種樣本):TP為x樣本被正確識(shí)別的數(shù)量;TN為y樣本被正確識(shí)別的數(shù)量;FP為分類器預(yù)測結(jié)果是x樣本,實(shí)際是y樣本,即誤報(bào)的y樣本數(shù)量;FN為分類器預(yù)測結(jié)果是y樣本,實(shí)際是x樣本,即漏報(bào)的x樣本數(shù)量。

    為了更加直觀地感受TL-ResNeXt的學(xué)習(xí)能力與分類效果,繪制了圖9所示的分類混淆矩陣,橫縱坐標(biāo)每個(gè)標(biāo)簽代表一類故障類型,對角線上的元素代表正確標(biāo)簽和預(yù)測標(biāo)簽的重合度,以驗(yàn)證集為例,每類故障樣本數(shù)為40,共計(jì)400份故障樣本。由圖可知,有4個(gè)誤分類樣本,分別為3個(gè)內(nèi)圈0.35 mm數(shù)據(jù)誤分成1個(gè)外圈0.35 mm數(shù)據(jù)和2個(gè)滾動(dòng)體0.53 mm數(shù)據(jù),1個(gè)滾動(dòng)體0.53 mm數(shù)據(jù)被誤分成外圈0.35 mm數(shù)據(jù),其余種類均達(dá)到滿分類效果,由此證明本模型是有效、可行的。

    圖9 軸承分類混淆矩陣Fig.9 Classification confusion matrix

    5.2 與其他診斷算法對比

    為進(jìn)一步說明模型的優(yōu)勢,將GADF-TL-ResNeXt與一些常用的方法進(jìn)行對比,結(jié)果如表3所示。與參考文獻(xiàn)[15~17]方法相比,所提的GADF-TL-ResNeXt方法可識(shí)別更多的故障類型,且提高了分類準(zhǔn)確度。

    表3 不同算法準(zhǔn)確率對比Tab.3 Comparison of accuracy of different algorithms

    5.3 變工況條件下軸承故障診斷

    5.3.1 實(shí)驗(yàn)平臺(tái)介紹

    實(shí)驗(yàn)數(shù)據(jù)來自軋鋼機(jī)模擬實(shí)驗(yàn)平臺(tái),圖10為軋機(jī)監(jiān)控系統(tǒng),包括軋機(jī)試驗(yàn)臺(tái)、傳感器采集系統(tǒng)和實(shí)驗(yàn)控制臺(tái)。它們分別如圖10中(Ⅰ)、(Ⅱ)、(Ⅲ)所示。2個(gè)加速計(jì)傳感器位于上工作輥的軸向和徑向位置,以收集軸承振動(dòng)數(shù)據(jù)。加載驅(qū)動(dòng)電機(jī)模擬軋機(jī)工作時(shí)軋輥所受到的載力(N0、N1、N2),由電機(jī)經(jīng)過傳動(dòng)系統(tǒng)給到軋輥不同的驅(qū)動(dòng)轉(zhuǎn)速(W0—600 r/min、W1—1 200 r/min、W2—1 800 r/min)來模擬不同工況下軋輥的轉(zhuǎn)速,采樣時(shí)間設(shè)置為20 s,采樣率為10 240 /s。

    圖10 軋鋼機(jī)模擬實(shí)驗(yàn)平臺(tái)Fig.10 Rolling mill simulation experimental platform

    根據(jù)實(shí)驗(yàn)設(shè)計(jì)的故障場景,這些數(shù)據(jù)用于軋機(jī)軸承的不同健康狀態(tài);收集的數(shù)據(jù)文件標(biāo)記為“正常”、“內(nèi)圈故障”、“外圈故障”和“滾動(dòng)體故障”;軸承故障如圖11所示。

    圖11 軸承各類健康狀態(tài)圖Fig.11 Bearing health charts

    5.3.2 數(shù)據(jù)處理及劃分

    采用控制變量法,即用固定轉(zhuǎn)速下的不同載力以及相同載力下不同轉(zhuǎn)速, 來模擬實(shí)際工況下的故障00000000000000000000000000000000000000000,考慮到采樣頻率為104/s,為了覆蓋3種轉(zhuǎn)速下軸承轉(zhuǎn)一圈所采集到的點(diǎn),根據(jù)式(8)可以將滑動(dòng)窗口大小設(shè)置為1 024,再使用半額重疊采樣法,制作相對應(yīng)的GADF編碼樣本。

    在轉(zhuǎn)速相同(W0、W1、W2)的情況下,分別采用2種電機(jī)加載(N0、N1)的軸承徑向振動(dòng)數(shù)據(jù)制作相應(yīng)的數(shù)據(jù)集,每類制作400份故障樣本,同樣按9:1劃分為訓(xùn)練集和驗(yàn)證集,每種轉(zhuǎn)速下包含6類故障狀態(tài)和1類健康狀態(tài),共7種類別,在加載相同(N0、N1、N2)的情況下,分別采用2種轉(zhuǎn)速(W1、W2)的軸承軸向振動(dòng)數(shù)據(jù)制作相應(yīng)的數(shù)據(jù)集,數(shù)據(jù)劃分結(jié)果與轉(zhuǎn)速相同情況下的結(jié)果保持一致,用來驗(yàn)證本模型的有效性。

    5.3.3 模型實(shí)驗(yàn)結(jié)果

    分別為相同轉(zhuǎn)速下施加不同載力、相同加載下給出不同轉(zhuǎn)速的2種情況,將振動(dòng)信號(hào)經(jīng)GADF編碼成二維圖像后,劃分訓(xùn)練集和驗(yàn)證集,輸入到TL-ResNeXt網(wǎng)絡(luò)中,得到如圖12所示的準(zhǔn)確率、訓(xùn)練損失曲線,圖13所示的分類混淆矩陣。

    圖12 驗(yàn)證準(zhǔn)確率及損失曲線Fig.12 Verification accuracy and loss curve

    圖13 分類混淆矩陣Fig.13 Classification confusion matrix

    由圖12(a)~圖12(c)可知,軸承在相同轉(zhuǎn)速下施加不同載力,模型識(shí)別最高準(zhǔn)確率分別為98.2%,98.6%,98.2%,平均準(zhǔn)確率為98.3%;由圖12(d)~圖12(f)可知,軸承在施加相同載力,再給出不同轉(zhuǎn)速的情況下,模型識(shí)別最高準(zhǔn)確率分別為97.5%,97.5%,97.9%,平均準(zhǔn)確率為97.6%;模型的總體平均準(zhǔn)確率為97.95%。最低訓(xùn)練損失平均值為0.03,模型達(dá)到收斂,由此可得出本模型在變工況的條件下依然可以達(dá)到一個(gè)很好的分類效果。

    圖13清晰地展示了軸承各種情況下的分類效果,由圖13(a)~圖13(c)可知,僅有個(gè)別故障被誤分為其他類別,其中軸承正常類型與軸承故障類型均被完全分開,由圖13(d)~圖13(f)可知,僅圖13(e)的2類標(biāo)簽分類效果不佳,原因是2號(hào)標(biāo)簽與1號(hào)標(biāo)簽的某段振動(dòng)數(shù)據(jù)較為相似;其余的情況僅個(gè)別數(shù)據(jù)發(fā)生誤判,由此說明本模型在變工況下分類是有效的。

    6 結(jié) 論

    本研究提出了一種以數(shù)據(jù)驅(qū)動(dòng)為主的CNN(TL-ResNeXt)軸承故障診斷方法,主要得出了以下幾點(diǎn)結(jié)論:

    1) 首先利用了一種振動(dòng)信號(hào)到圖像的轉(zhuǎn)換方式-GAF編碼,拋棄了以原始振動(dòng)信號(hào)作為輸入,既保留了信號(hào)時(shí)間序列的相關(guān)性,又可以充分提取軸承不同故障的特征,無需人工手動(dòng)提取。

    2) 利用了層級更深的殘差網(wǎng)絡(luò)作為模型主體,殘差塊部分采用分組卷積情況,通過更深的網(wǎng)絡(luò)層和更好的特征提取層,ResNeXt將提高故障診斷的最終預(yù)測精度,并結(jié)合了遷移學(xué)習(xí)模塊,在較少的訓(xùn)練輪數(shù)下模型可以達(dá)到更好的穩(wěn)定性和準(zhǔn)確率。

    3) 在西儲(chǔ)大學(xué)軸承數(shù)據(jù)集上對比了其他幾種流行的網(wǎng)絡(luò)模型,本方法表現(xiàn)的更好;在軋機(jī)模擬試驗(yàn)臺(tái)上模擬軸承在變工況條件下情況,依然能達(dá)到一個(gè)良好的分類效果。在軸承給出相同轉(zhuǎn)速下施加不同載力的情況下,模型識(shí)別平均準(zhǔn)確率為98.3%,在施加相同載力下給出不同轉(zhuǎn)速的情況下,模型識(shí)別平均準(zhǔn)確率為97.6%,僅有個(gè)別類型發(fā)生誤判,因此可見本模型在故障診斷領(lǐng)域有著良好的性能。

    猜你喜歡
    殘差分組故障診斷
    基于雙向GRU與殘差擬合的車輛跟馳建模
    基于殘差學(xué)習(xí)的自適應(yīng)無人機(jī)目標(biāo)跟蹤算法
    基于遞歸殘差網(wǎng)絡(luò)的圖像超分辨率重建
    分組搭配
    怎么分組
    分組
    因果圖定性分析法及其在故障診斷中的應(yīng)用
    平穩(wěn)自相關(guān)過程的殘差累積和控制圖
    河南科技(2015年8期)2015-03-11 16:23:52
    基于LCD和排列熵的滾動(dòng)軸承故障診斷
    基于WPD-HHT的滾動(dòng)軸承故障診斷
    免费黄网站久久成人精品| 久久午夜福利片| 国产精品一及| 成人毛片a级毛片在线播放| 亚洲av中文字字幕乱码综合| 久久99热6这里只有精品| 69av精品久久久久久| 大码成人一级视频| 欧美日韩在线观看h| 国内少妇人妻偷人精品xxx网站| 免费看不卡的av| 国产一区有黄有色的免费视频| 欧美激情在线99| 免费av毛片视频| 亚洲av日韩在线播放| 欧美 日韩 精品 国产| 人妻少妇偷人精品九色| eeuss影院久久| 国产成人freesex在线| 性色avwww在线观看| 91久久精品国产一区二区三区| 日韩欧美精品v在线| 人体艺术视频欧美日本| 国产在线男女| 又爽又黄无遮挡网站| 亚洲va在线va天堂va国产| 男女边吃奶边做爰视频| 亚洲成人精品中文字幕电影| 精品国产一区二区三区久久久樱花 | 国产探花在线观看一区二区| 亚洲国产精品国产精品| 亚洲美女搞黄在线观看| 在线播放无遮挡| av女优亚洲男人天堂| 综合色av麻豆| 九九久久精品国产亚洲av麻豆| 色综合色国产| 黄色一级大片看看| 国产高清三级在线| 日本欧美国产在线视频| 18禁在线播放成人免费| 欧美bdsm另类| 草草在线视频免费看| 国产一区二区三区av在线| 五月开心婷婷网| 国产精品不卡视频一区二区| 成年女人在线观看亚洲视频 | 亚洲经典国产精华液单| 在线精品无人区一区二区三 | 国产精品麻豆人妻色哟哟久久| 中文乱码字字幕精品一区二区三区| 亚洲婷婷狠狠爱综合网| 色哟哟·www| 亚洲国产欧美在线一区| 大片电影免费在线观看免费| 亚洲av成人精品一二三区| 亚洲成人av在线免费| 18+在线观看网站| 伊人久久国产一区二区| 国产成人aa在线观看| 国产成人a∨麻豆精品| 亚洲精品国产色婷婷电影| 激情 狠狠 欧美| 亚洲真实伦在线观看| av.在线天堂| 波野结衣二区三区在线| 国产午夜精品一二区理论片| 看黄色毛片网站| 亚洲av成人精品一区久久| 中文在线观看免费www的网站| 婷婷色综合www| 久久久精品欧美日韩精品| 久久鲁丝午夜福利片| 欧美极品一区二区三区四区| 亚洲人与动物交配视频| 国产精品嫩草影院av在线观看| 国产免费视频播放在线视频| 久久久久久久精品精品| 麻豆国产97在线/欧美| 免费播放大片免费观看视频在线观看| 日韩av在线免费看完整版不卡| 在线天堂最新版资源| 欧美极品一区二区三区四区| 国产精品99久久99久久久不卡 | 午夜福利视频1000在线观看| 视频中文字幕在线观看| xxx大片免费视频| 亚洲av免费在线观看| 国产一区二区三区av在线| 亚洲精品自拍成人| 国内精品美女久久久久久| www.av在线官网国产| 中文在线观看免费www的网站| 人妻系列 视频| 亚洲国产精品成人久久小说| 久久国内精品自在自线图片| 大片免费播放器 马上看| 卡戴珊不雅视频在线播放| 一个人观看的视频www高清免费观看| 亚洲精品久久午夜乱码| 亚洲色图综合在线观看| 亚洲精品影视一区二区三区av| 欧美性猛交╳xxx乱大交人| 又粗又硬又长又爽又黄的视频| 国产男人的电影天堂91| 日本与韩国留学比较| 国产一区二区三区av在线| 免费观看a级毛片全部| 婷婷色综合大香蕉| 国产精品成人在线| 国产爱豆传媒在线观看| 九九在线视频观看精品| 亚洲天堂av无毛| 成人免费观看视频高清| 人人妻人人爽人人添夜夜欢视频 | 国产又色又爽无遮挡免| 精品久久久噜噜| 亚洲精品一区蜜桃| 久热久热在线精品观看| 人妻系列 视频| 少妇人妻精品综合一区二区| 插阴视频在线观看视频| av在线播放精品| 亚洲四区av| 国产精品嫩草影院av在线观看| 国模一区二区三区四区视频| 91精品伊人久久大香线蕉| av在线播放精品| 91精品伊人久久大香线蕉| 丰满少妇做爰视频| 色播在线永久视频| 美女福利国产在线| 国语对白做爰xxxⅹ性视频网站| 我要看黄色一级片免费的| av国产精品久久久久影院| 久久久久久久精品精品| 精品一区二区三区av网在线观看 | 老鸭窝网址在线观看| 人成视频在线观看免费观看| 男女边摸边吃奶| 亚洲精品日韩在线中文字幕| 看非洲黑人一级黄片| 少妇被粗大猛烈的视频| 91国产中文字幕| 91国产中文字幕| 国产精品成人在线| 街头女战士在线观看网站| 超色免费av| 国产精品国产三级专区第一集| 黄片无遮挡物在线观看| 国产 一区精品| 人人妻人人澡人人看| 国产精品秋霞免费鲁丝片| 国产成人免费观看mmmm| 亚洲国产av新网站| 亚洲国产欧美在线一区| 我的亚洲天堂| 日日摸夜夜添夜夜爱| 啦啦啦在线免费观看视频4| 最近中文字幕高清免费大全6| 校园人妻丝袜中文字幕| 亚洲精品久久午夜乱码| 一区二区三区四区激情视频| 中文字幕av电影在线播放| av在线播放精品| 中文字幕高清在线视频| 久久精品国产综合久久久| 人妻人人澡人人爽人人| 我要看黄色一级片免费的| 不卡av一区二区三区| 在线 av 中文字幕| 久久久久精品性色| 一二三四在线观看免费中文在| 最近的中文字幕免费完整| 亚洲五月色婷婷综合| 韩国高清视频一区二区三区| www.自偷自拍.com| 亚洲精品一区蜜桃| 在线观看免费视频网站a站| 国产一区二区三区综合在线观看| 搡老岳熟女国产| 91成人精品电影| av在线app专区| 欧美国产精品va在线观看不卡| 9191精品国产免费久久| 狂野欧美激情性xxxx| 丰满乱子伦码专区| 黄片播放在线免费| 啦啦啦在线免费观看视频4| 一级毛片我不卡| 99久国产av精品国产电影| 国产精品国产三级国产专区5o| 天天躁夜夜躁狠狠躁躁| 欧美97在线视频| 久热这里只有精品99| 国产欧美亚洲国产| 中文天堂在线官网| 亚洲精品日韩在线中文字幕| 看非洲黑人一级黄片| 成年女人毛片免费观看观看9 | 中文字幕最新亚洲高清| 国产日韩一区二区三区精品不卡| 黄色怎么调成土黄色| 在线观看免费高清a一片| 人妻人人澡人人爽人人| 成年人午夜在线观看视频| 丝袜在线中文字幕| 成年av动漫网址| 国产野战对白在线观看| 一区福利在线观看| 欧美人与善性xxx| 性高湖久久久久久久久免费观看| 色网站视频免费| 亚洲国产精品一区二区三区在线| 亚洲一区二区三区欧美精品| 日韩大码丰满熟妇| 在线观看免费午夜福利视频| 国产精品亚洲av一区麻豆 | 亚洲国产精品国产精品| 欧美黑人欧美精品刺激| 最新的欧美精品一区二区| 亚洲成国产人片在线观看| av网站在线播放免费| 国产高清国产精品国产三级| 免费观看性生交大片5| 好男人视频免费观看在线| 在线免费观看不下载黄p国产| 女人被躁到高潮嗷嗷叫费观| 成年动漫av网址| 大陆偷拍与自拍| 成人黄色视频免费在线看| 欧美中文综合在线视频| 欧美日韩视频高清一区二区三区二| 9色porny在线观看| 高清不卡的av网站| 女性被躁到高潮视频| 欧美亚洲日本最大视频资源| 国产精品久久久av美女十八| 中文字幕人妻丝袜一区二区 | 综合色丁香网| 不卡av一区二区三区| 中文字幕亚洲精品专区| 亚洲第一av免费看| 黄色一级大片看看| 亚洲七黄色美女视频| 91国产中文字幕| 精品国产一区二区久久| 激情五月婷婷亚洲| 日韩欧美精品免费久久| 在线免费观看不下载黄p国产| 波多野结衣av一区二区av| 亚洲国产日韩一区二区| 老鸭窝网址在线观看| 亚洲伊人久久精品综合| 亚洲情色 制服丝袜| 亚洲四区av| 久久久久久久大尺度免费视频| 国产免费又黄又爽又色| 国产成人一区二区在线| 大陆偷拍与自拍| 日韩大码丰满熟妇| 99久久精品国产亚洲精品| 久久ye,这里只有精品| 好男人视频免费观看在线| 国产黄色视频一区二区在线观看| 欧美日韩亚洲高清精品| 国产精品麻豆人妻色哟哟久久| 人人妻人人爽人人添夜夜欢视频| 国产有黄有色有爽视频| 免费人妻精品一区二区三区视频| 99热网站在线观看| 男女高潮啪啪啪动态图| 亚洲婷婷狠狠爱综合网| 一级黄片播放器| 亚洲成人免费av在线播放| 中文字幕人妻丝袜一区二区 | 老司机亚洲免费影院| 一级a爱视频在线免费观看| 熟妇人妻不卡中文字幕| 国产精品久久久人人做人人爽| 女性生殖器流出的白浆| 日韩av在线免费看完整版不卡| 日韩不卡一区二区三区视频在线| 午夜激情久久久久久久| 毛片一级片免费看久久久久| 看免费av毛片| 9热在线视频观看99| 亚洲国产精品国产精品| 国产欧美亚洲国产| 波多野结衣一区麻豆| 丰满乱子伦码专区| 9热在线视频观看99| 亚洲国产毛片av蜜桃av| 国产有黄有色有爽视频| 在线观看免费午夜福利视频| 国产精品女同一区二区软件| 国产精品国产三级国产专区5o| 亚洲国产毛片av蜜桃av| 美女脱内裤让男人舔精品视频| 操出白浆在线播放| 国产99久久九九免费精品| 性色av一级| 国产成人午夜福利电影在线观看| 叶爱在线成人免费视频播放| 国产 一区精品| 国产精品欧美亚洲77777| 精品少妇内射三级| 另类精品久久| 亚洲精品久久成人aⅴ小说| 99国产综合亚洲精品| 别揉我奶头~嗯~啊~动态视频 | 久久精品久久精品一区二区三区| 一级毛片我不卡| 亚洲国产成人一精品久久久| 97在线人人人人妻| 少妇被粗大的猛进出69影院| 欧美日本中文国产一区发布| 亚洲精品aⅴ在线观看| 日本猛色少妇xxxxx猛交久久| av不卡在线播放| 97在线人人人人妻| 精品国产乱码久久久久久男人| 如日韩欧美国产精品一区二区三区| 高清av免费在线| 欧美少妇被猛烈插入视频| 王馨瑶露胸无遮挡在线观看| 成年av动漫网址| 最近最新中文字幕大全免费视频 | 哪个播放器可以免费观看大片| 中国国产av一级| 日韩中文字幕欧美一区二区 | videos熟女内射| 最近中文字幕2019免费版| 日韩伦理黄色片| 综合色丁香网| 国产乱来视频区| 午夜福利一区二区在线看| av.在线天堂| 色婷婷av一区二区三区视频| av国产精品久久久久影院| 黄网站色视频无遮挡免费观看| 丝袜人妻中文字幕| 在线观看三级黄色| 日本欧美国产在线视频| 国产精品蜜桃在线观看| 欧美黄色片欧美黄色片| 亚洲av国产av综合av卡| 多毛熟女@视频| 亚洲欧美中文字幕日韩二区| 丝袜喷水一区| 美女脱内裤让男人舔精品视频| 美女主播在线视频| 国产激情久久老熟女| 大片免费播放器 马上看| 伊人久久国产一区二区| 精品国产露脸久久av麻豆| 日韩制服骚丝袜av| 中文字幕精品免费在线观看视频| 九九爱精品视频在线观看| 伦理电影大哥的女人| 亚洲一级一片aⅴ在线观看| 色视频在线一区二区三区| 久久久久久久国产电影| 亚洲国产日韩一区二区| 亚洲欧洲国产日韩| 91精品国产国语对白视频| 欧美亚洲日本最大视频资源| 人人妻人人爽人人添夜夜欢视频| 超色免费av| 成人漫画全彩无遮挡| 一边摸一边做爽爽视频免费| 亚洲欧美成人精品一区二区| 视频在线观看一区二区三区| 狠狠精品人妻久久久久久综合| 亚洲国产日韩一区二区| 69精品国产乱码久久久| 国精品久久久久久国模美| 三上悠亚av全集在线观看| av有码第一页| 成人国产麻豆网| 精品卡一卡二卡四卡免费| 人人妻,人人澡人人爽秒播 | 侵犯人妻中文字幕一二三四区| 国产又色又爽无遮挡免| 免费人妻精品一区二区三区视频| 久久99热这里只频精品6学生| 亚洲av在线观看美女高潮| 久久久久精品人妻al黑| 欧美日韩福利视频一区二区| 久久精品国产综合久久久| 日韩视频在线欧美| 国产黄频视频在线观看| 久久精品久久精品一区二区三区| 999精品在线视频| 无遮挡黄片免费观看| 天天添夜夜摸| 亚洲免费av在线视频| av卡一久久| 一级a爱视频在线免费观看| 亚洲精品久久成人aⅴ小说| 国产精品一区二区在线观看99| 日本欧美国产在线视频| 成年av动漫网址| av.在线天堂| 日韩av免费高清视频| 日韩 亚洲 欧美在线| 人妻一区二区av| 日韩制服骚丝袜av| 中文字幕精品免费在线观看视频| 秋霞伦理黄片| 亚洲人成电影观看| 国产成人a∨麻豆精品| 美女午夜性视频免费| av在线观看视频网站免费| 国产免费视频播放在线视频| 久久人人97超碰香蕉20202| 国产成人精品在线电影| 极品人妻少妇av视频| xxxhd国产人妻xxx| 亚洲图色成人| 精品免费久久久久久久清纯 | 国产亚洲一区二区精品| 亚洲精品美女久久av网站| 菩萨蛮人人尽说江南好唐韦庄| 无限看片的www在线观看| 久久婷婷青草| 精品一区二区三区四区五区乱码 | 亚洲国产看品久久| 亚洲av中文av极速乱| 观看美女的网站| 三上悠亚av全集在线观看| 丝袜美腿诱惑在线| 国产男女内射视频| 日韩中文字幕欧美一区二区 | 人体艺术视频欧美日本| 香蕉国产在线看| 欧美乱码精品一区二区三区| 波多野结衣av一区二区av| 一级,二级,三级黄色视频| 久久 成人 亚洲| 日韩一区二区视频免费看| 国产熟女午夜一区二区三区| 日韩免费高清中文字幕av| 人人妻人人爽人人添夜夜欢视频| 久久这里只有精品19| 亚洲自偷自拍图片 自拍| 国产成人午夜福利电影在线观看| 精品一品国产午夜福利视频| 中文字幕人妻丝袜一区二区 | 国产精品女同一区二区软件| 丝瓜视频免费看黄片| 日本欧美视频一区| 亚洲激情五月婷婷啪啪| 99久久综合免费| 午夜福利一区二区在线看| 视频区图区小说| 欧美日韩精品网址| 午夜福利,免费看| 中文字幕色久视频| 免费看av在线观看网站| 黄频高清免费视频| www.熟女人妻精品国产| 黑人巨大精品欧美一区二区蜜桃| 中文天堂在线官网| 看免费av毛片| 国产精品麻豆人妻色哟哟久久| 国产极品天堂在线| 国产免费一区二区三区四区乱码| 一区福利在线观看| av福利片在线| 亚洲国产精品成人久久小说| 精品人妻熟女毛片av久久网站| 日本vs欧美在线观看视频| 一二三四在线观看免费中文在| 黄色一级大片看看| 亚洲精品国产一区二区精华液| 欧美 日韩 精品 国产| 满18在线观看网站| 啦啦啦中文免费视频观看日本| 在线观看一区二区三区激情| 18在线观看网站| 人妻人人澡人人爽人人| 亚洲精品一二三| 99久久精品国产亚洲精品| 韩国av在线不卡| 美女脱内裤让男人舔精品视频| 国产精品国产av在线观看| 十八禁人妻一区二区| 天天躁日日躁夜夜躁夜夜| 人妻一区二区av| 午夜老司机福利片| 少妇人妻精品综合一区二区| 香蕉丝袜av| 国精品久久久久久国模美| 男女午夜视频在线观看| 亚洲成人av在线免费| 久久青草综合色| 亚洲伊人色综图| 国产成人精品在线电影| 久久午夜综合久久蜜桃| tube8黄色片| 亚洲一卡2卡3卡4卡5卡精品中文| 男人爽女人下面视频在线观看| 美女高潮到喷水免费观看| 国产片内射在线| 国产一区有黄有色的免费视频| 99精品久久久久人妻精品| 97精品久久久久久久久久精品| 午夜福利免费观看在线| 国产精品 欧美亚洲| videos熟女内射| 欧美人与性动交α欧美精品济南到| 国产精品久久久久成人av| www.av在线官网国产| 久久久久精品国产欧美久久久 | 婷婷成人精品国产| 国产免费福利视频在线观看| 最近中文字幕高清免费大全6| 日韩精品免费视频一区二区三区| 国产一区二区在线观看av| 侵犯人妻中文字幕一二三四区| 最近中文字幕2019免费版| 日韩视频在线欧美| 国产精品无大码| 欧美日韩精品网址| 成人亚洲欧美一区二区av| 人人妻人人添人人爽欧美一区卜| 国产精品香港三级国产av潘金莲 | 视频在线观看一区二区三区| 国产在视频线精品| 黄色视频不卡| 丰满饥渴人妻一区二区三| 在线天堂中文资源库| 国产人伦9x9x在线观看| 免费少妇av软件| 亚洲国产av影院在线观看| 亚洲精品国产区一区二| 国产精品 国内视频| 又黄又粗又硬又大视频| 国产成人精品福利久久| 久久热在线av| 80岁老熟妇乱子伦牲交| 日韩中文字幕视频在线看片| 国产深夜福利视频在线观看| 一区在线观看完整版| 亚洲精品第二区| 满18在线观看网站| 欧美乱码精品一区二区三区| 国产欧美亚洲国产| 丰满迷人的少妇在线观看| 久热这里只有精品99| 久久久久视频综合| 国产乱人偷精品视频| 精品国产乱码久久久久久小说| 久久久久久人人人人人| 中文欧美无线码| 国产精品嫩草影院av在线观看| 波野结衣二区三区在线| 丝袜美足系列| 大香蕉久久网| 久久综合国产亚洲精品| 少妇被粗大的猛进出69影院| av网站在线播放免费| 日韩精品免费视频一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 久久久久久人人人人人| 中文欧美无线码| 日韩中文字幕欧美一区二区 | 午夜福利免费观看在线| 在线观看免费高清a一片| 岛国毛片在线播放| 色视频在线一区二区三区| 欧美激情 高清一区二区三区| 一级毛片电影观看| 欧美黑人精品巨大| 国产精品免费大片| 秋霞伦理黄片| 国产精品欧美亚洲77777| 久热这里只有精品99| 多毛熟女@视频| 久热这里只有精品99| 一级毛片电影观看| 另类亚洲欧美激情| 搡老乐熟女国产| a级片在线免费高清观看视频| 夫妻午夜视频| 久久久精品免费免费高清| 男女国产视频网站| 午夜免费男女啪啪视频观看| 丰满乱子伦码专区| 国产精品久久久人人做人人爽| 久久久精品94久久精品| 久久精品aⅴ一区二区三区四区| 国产激情久久老熟女| 成年动漫av网址| 男女边吃奶边做爰视频| 秋霞伦理黄片| 久久久精品94久久精品| 亚洲欧美一区二区三区黑人| 超碰成人久久| 高清黄色对白视频在线免费看| 老司机影院成人| 午夜福利网站1000一区二区三区| 一边摸一边抽搐一进一出视频| 交换朋友夫妻互换小说| 中文乱码字字幕精品一区二区三区| 久久99一区二区三区| 久久久久网色| 一级片免费观看大全| 秋霞在线观看毛片| 国产精品成人在线| 国产精品久久久久久久久免| 看非洲黑人一级黄片| 韩国精品一区二区三区| 99久久精品国产亚洲精品|