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

    動(dòng)態(tài)寬卷積殘差網(wǎng)絡(luò)的軸承故障診斷方法

    2023-09-27 12:53:56秦國(guó)浩黃鋒飛丁國(guó)富
    中國(guó)機(jī)械工程 2023年18期
    關(guān)鍵詞:殘差故障診斷軸承

    秦國(guó)浩 張 楷,2 丁 昆 黃鋒飛 鄭 慶,2 丁國(guó)富,2

    1.西南交通大學(xué)機(jī)械工程學(xué)院,成都,6100312.西南交通大學(xué)軌道交通運(yùn)維技術(shù)與裝備四川省重點(diǎn)實(shí)驗(yàn)室,成都,610031

    0 引言

    軸承作為旋轉(zhuǎn)機(jī)械的關(guān)鍵零部件之一,其運(yùn)行狀態(tài)直接影響著機(jī)械裝備的運(yùn)行狀態(tài)[1]。在長(zhǎng)期的運(yùn)行過(guò)程中軸承承受交變載荷和疲勞損傷,易產(chǎn)生磨損、點(diǎn)蝕、剝落等破壞,以致機(jī)械設(shè)備停機(jī),造成巨大經(jīng)濟(jì)損失和人員傷亡等危害[2]。因此,狀態(tài)信號(hào)監(jiān)測(cè)和智能診斷對(duì)軸承長(zhǎng)期健康運(yùn)行具有重要研究意義和應(yīng)用價(jià)值[3]。

    以支持向量機(jī)、淺層神經(jīng)網(wǎng)絡(luò)等為代表的傳統(tǒng)機(jī)器學(xué)習(xí)故障診斷方法是將獲取的軸承振動(dòng)信號(hào)進(jìn)行人工特征提取和篩選,然后通過(guò)分類(lèi)器對(duì)軸承健康狀態(tài)進(jìn)行識(shí)別。如楊宇等[4]使用經(jīng)驗(yàn)?zāi)B(tài)分解方法提取軸承故障特征進(jìn)行診斷,完成了對(duì)軸承工作狀態(tài)和故障類(lèi)別的有效精確識(shí)別,但是此類(lèi)方法存在依賴(lài)專(zhuān)家經(jīng)驗(yàn)、非線性特征提取能力差、故障診斷精度不足等問(wèn)題。隨著計(jì)算能力的提高和人工智能的發(fā)展,神經(jīng)網(wǎng)絡(luò)因其強(qiáng)大的非線性特征自適應(yīng)提取能力受到機(jī)械故障診斷領(lǐng)域的廣泛研究[5]。例如,侯文擎等[6]在自編碼網(wǎng)絡(luò)中引入“稀疏”和“損傷加噪”構(gòu)成了DAE(denosing autoencoder)網(wǎng)絡(luò),提取的數(shù)據(jù)特征魯棒性強(qiáng);趙光權(quán)等[7]使用深度置信網(wǎng)絡(luò)對(duì)原始時(shí)域信號(hào)直接進(jìn)行故障提取和健康狀況識(shí)別,取得了較高的診斷精度。然而,上述神經(jīng)網(wǎng)絡(luò)局部特征提取能力較弱,無(wú)法高效提取數(shù)據(jù)深層表示特征,制約了其在工業(yè)中的應(yīng)用與發(fā)展。卷積神經(jīng)網(wǎng)絡(luò)為解決此問(wèn)題而設(shè)計(jì)的局部感知和權(quán)值共享能有效減少模型參數(shù)量、抑制過(guò)擬合。一些學(xué)者也針對(duì)一維振動(dòng)信號(hào)作為輸入的軸承故障診斷任務(wù)進(jìn)行了適應(yīng)性改進(jìn)[8-9]。

    近年來(lái),卷積神經(jīng)網(wǎng)絡(luò)已成為軸承故障診斷主要智能算法框架之一,許多學(xué)者從不同方面對(duì)卷積神經(jīng)網(wǎng)絡(luò)進(jìn)行了大量改進(jìn),以提高其故障診斷精度。在深度學(xué)習(xí)理論中,通常認(rèn)為網(wǎng)絡(luò)層數(shù)的增加有利于深層非線性抽象特征的提取。趙小強(qiáng)等[10]在軸承故障診斷中引入神經(jīng)網(wǎng)絡(luò)Alexnet,并將一維振動(dòng)信號(hào)按橫向插樣構(gòu)建為二維特征圖以保留軸承故障信息卷積運(yùn)算時(shí)域信號(hào)的關(guān)聯(lián)性;趙宇凱等[11]使用連續(xù)的3×3卷積核替代原本的卷積核,再反復(fù)疊加卷積層和池化層增加網(wǎng)絡(luò)的深度,使得該方法能夠更好地提取滾動(dòng)軸承振動(dòng)信號(hào)的更深層特征;溫江濤等[12]使用堆疊的一維卷積自編碼的殘差結(jié)構(gòu),搭建無(wú)監(jiān)督領(lǐng)域自適應(yīng)的遷移學(xué)習(xí)網(wǎng)絡(luò)模型,較好地實(shí)現(xiàn)了源域到目標(biāo)域特征知識(shí)遷移,提高了軸承故障診斷精度。然而,當(dāng)網(wǎng)絡(luò)層數(shù)超過(guò)一定范圍后,網(wǎng)絡(luò)輸入或梯度信息在傳播過(guò)程中可能會(huì)消失,進(jìn)而導(dǎo)致模型診斷能力下降甚至失效。深度殘差網(wǎng)絡(luò)(deep residual network,ResNet)通過(guò)引入跨層連接擬合殘差項(xiàng)有效解決了深度神經(jīng)網(wǎng)絡(luò)深度增加所帶來(lái)的梯度消失和梯度爆炸等問(wèn)題[13]。趙靖等[14]使用殘差注意力卷積神經(jīng)網(wǎng)絡(luò)方法,有效提取源域故障特征,提高了遷移學(xué)習(xí)方法下的軸承故障診斷精度。除了增加網(wǎng)絡(luò)層數(shù),也可將跨層連接擴(kuò)展至所有層直接相連,以進(jìn)一步確保網(wǎng)絡(luò)中各層之間的最大信息流傳遞[15]。熊鵬等[16]結(jié)合動(dòng)態(tài)加權(quán)密集連接卷積網(wǎng)絡(luò)對(duì)不同頻帶信息自適應(yīng)加權(quán),加強(qiáng)了密集連接卷積網(wǎng)絡(luò)對(duì)變轉(zhuǎn)速工況的特征表征能力,提高了變轉(zhuǎn)速下故障識(shí)別率。此外,為使信息感知多樣化,采用多尺度卷積設(shè)計(jì)是一種可行的方案。如吳靜然等[17]針對(duì)噪聲條件下故障識(shí)別準(zhǔn)確率偏低的問(wèn)題,提出了一種基于注意力機(jī)制的多尺度端到端故障診斷方法;趙小強(qiáng)等[18]也通過(guò)引入Inception網(wǎng)絡(luò)結(jié)構(gòu)和注意力機(jī)制對(duì)卷積神經(jīng)網(wǎng)絡(luò)進(jìn)行改進(jìn),使網(wǎng)絡(luò)能最大化地提取故障數(shù)據(jù)中的信息。

    以上針對(duì)卷積神經(jīng)網(wǎng)絡(luò)及擴(kuò)展網(wǎng)絡(luò)的改進(jìn)方法在軸承故障診斷中取得了較好的效果,但在實(shí)際工程中,往往不能直接對(duì)軸承振動(dòng)進(jìn)行接觸測(cè)量,一些微弱早期故障信號(hào)在傳遞過(guò)程中易受到相關(guān)部件運(yùn)行噪聲干擾,使得其故障特征淹沒(méi)在噪聲中[19]。而現(xiàn)有診斷模型特征提取范圍仍受其卷積核感受野大小的制約,統(tǒng)一尺度的卷積核也難以充分捕捉軸承故障振動(dòng)信號(hào)中的不同沖擊頻率成分,進(jìn)而阻礙故障診斷準(zhǔn)確率的進(jìn)一步提高。本文提出一種基于動(dòng)態(tài)寬卷積殘差網(wǎng)絡(luò)(dynamic wide-kernels residual network,DWResNet)的軸承故障診斷方法,該方法在一維深度殘差網(wǎng)絡(luò)基礎(chǔ)上引入寬殘差核結(jié)構(gòu)形成雙通道并行網(wǎng)絡(luò)結(jié)構(gòu),并通過(guò)注意力機(jī)制對(duì)卷積核進(jìn)行動(dòng)態(tài)加權(quán),自適應(yīng)地充分提取不同尺度特征信息,以實(shí)現(xiàn)強(qiáng)噪聲背景下軸承故障的高精度診斷。

    1 卷積神經(jīng)網(wǎng)絡(luò)基礎(chǔ)理論

    1.1 卷積層

    卷積層是卷積神經(jīng)網(wǎng)絡(luò)中最為核心的組成部分,每一個(gè)卷積層都由多個(gè)卷積核組合而成。相較于全連接網(wǎng)絡(luò),卷積神經(jīng)網(wǎng)絡(luò)具有數(shù)據(jù)局部感知和權(quán)值參數(shù)共享等優(yōu)勢(shì),能夠有效減少網(wǎng)絡(luò)的參數(shù)量。在輸入的振動(dòng)信號(hào)數(shù)據(jù)上,卷積核通過(guò)滑動(dòng)不同的步長(zhǎng)來(lái)提取相應(yīng)位置的局部特征。同時(shí),可以使用一個(gè)卷積核來(lái)提取不同位置的同類(lèi)局部特征進(jìn)行權(quán)值共享,從而實(shí)現(xiàn)卷積網(wǎng)絡(luò)中參數(shù)量的有效利用。卷積的運(yùn)算可表示為

    (1)

    式中,Xl為第l層的輸出;W(c)為第c個(gè)卷積核的權(quán)重矩陣;*為卷積運(yùn)算;b為權(quán)重偏置。

    為了使得卷積層具有不同的感受野,可通過(guò)設(shè)置不同的卷積核大小來(lái)實(shí)現(xiàn)。卷積核感受野的計(jì)算公式為

    (2)

    式中,Ll為第l層對(duì)應(yīng)的感受野大小;kl為第l層的卷積核大小;si為第i層的移動(dòng)步長(zhǎng)。

    1.2 空洞卷積

    空洞卷積就是在標(biāo)準(zhǔn)的卷積中加入空洞,擴(kuò)大卷積層的感受野,且不失對(duì)特征的分辨提取能力。在不增加卷積層層數(shù)的前提下,通過(guò)設(shè)置不同的空洞率,形成不同感受野大小的空洞卷積。

    假設(shè)原來(lái)第l-1層的卷積核大小為kl-1,則卷積核空洞率為d的一層空洞卷積的感受野大小Ll為

    Ll=kl-1+(kl-1-1)(d-1)

    (3)

    假設(shè)空洞卷積的步長(zhǎng)為s,那么經(jīng)過(guò)空洞卷積后的特征圖大小m的計(jì)算公式為

    (4)

    式中,x為當(dāng)前層的輸入;p為填充操作的尺寸。

    通過(guò)設(shè)置不同的空洞率d使得原本的卷積核擁有更大的感受野,同時(shí)配合填充操作,使得輸出的特征圖大小m能夠與原卷積層輸出的特征圖大小一致,從而保證不同卷積核輸出的特征圖能夠進(jìn)行信息融合。

    對(duì)于一維的滾動(dòng)軸承振動(dòng)信號(hào),通過(guò)設(shè)置空洞卷積能夠使得網(wǎng)絡(luò)在提取不同特征區(qū)間有效信息的同時(shí)保證網(wǎng)絡(luò)參數(shù)量的緩慢增長(zhǎng)。

    1.3 殘差網(wǎng)絡(luò)(ResNet)

    ResNet的提出使得在網(wǎng)絡(luò)深度增加的同時(shí)獲得更高的收益成為可能。ResNet不是直接學(xué)習(xí)多層堆疊的函數(shù)擬合的直接映射,而是學(xué)習(xí)它們的殘差映射。圖1展示了一個(gè)標(biāo)準(zhǔn)的殘差結(jié)構(gòu)(x表示輸入,σ表示激活函數(shù)),通過(guò)學(xué)習(xí)F(x)函數(shù)的殘差映射而非直接映射F(x)+x來(lái)有效減少網(wǎng)絡(luò)由于深度增加帶來(lái)的梯度問(wèn)題以及對(duì)直接映射的過(guò)擬合問(wèn)題。

    圖1 標(biāo)準(zhǔn)的殘差塊Fig.1 Residual block

    ResNet模型的核心部分由許多基本的殘差塊堆疊而成,如圖2所示。每一個(gè)基本殘差塊中,均由兩個(gè)卷積層、兩個(gè)ReLU[20]激活函數(shù)、兩個(gè)批量正則化(batch normalization,BN)[21]和一個(gè)恒等映射連接而成。當(dāng)前后通道數(shù)發(fā)生改變時(shí),通過(guò)一個(gè)圖2b所示的卷積層來(lái)適配改變的通道數(shù)。其中,采用BN層的目的是為了加快網(wǎng)絡(luò)收斂,提高網(wǎng)絡(luò)的泛化性,它的計(jì)算公式如下[21]:

    (5)

    式中,E[x]為x的期望;Var[x]為x的方差;γ、β均為可訓(xùn)練的超參數(shù);ε為一個(gè)趨近于0的常量。

    ReLU是常見(jiàn)的激活函數(shù),其主要作用是減小網(wǎng)絡(luò)計(jì)算量,緩解過(guò)擬合,其計(jì)算公式如下:

    (6)

    通過(guò)上述殘差連接的結(jié)構(gòu)可使得網(wǎng)絡(luò)向更深的層次發(fā)展,同時(shí)也能一定程度地緩解網(wǎng)絡(luò)的退化、梯度消失和梯度爆炸等問(wèn)題。

    2 本文所提方法

    為了使得網(wǎng)絡(luò)更易于提取振動(dòng)信號(hào)的特征信息,減少信號(hào)強(qiáng)噪聲干擾所帶來(lái)的影響,本文提出了動(dòng)態(tài)寬卷積殘差網(wǎng)絡(luò)對(duì)軸承振動(dòng)信號(hào)進(jìn)行學(xué)習(xí)。具體來(lái)說(shuō),本文首先在殘差網(wǎng)絡(luò)的基礎(chǔ)上,對(duì)殘差塊的每層卷積進(jìn)行改進(jìn),然后引入可選擇核網(wǎng)絡(luò)(selective kernel networks,SKNet)[22]中可選擇性操作進(jìn)行改進(jìn),使得網(wǎng)絡(luò)能夠在3×1的卷積核和寬核的空洞卷積之間進(jìn)行自適應(yīng)的選擇,最后按照不同的權(quán)重融合兩個(gè)不同卷積核的輸出特征作為每層的輸出特征進(jìn)行學(xué)習(xí)。

    2.1 動(dòng)態(tài)寬卷積殘差網(wǎng)絡(luò)(DWResNet)

    本文所提方法的單層結(jié)構(gòu)如圖3所示,對(duì)于任意一個(gè)給定的特征圖輸入X(X∈RC×H×W,C為輸入數(shù)據(jù)的維度,H、W分別為輸入數(shù)據(jù)的高和寬),將其同時(shí)輸入到兩個(gè)卷積核大小分別為3×1和空洞率d為14的3×1空洞卷積中,使其在擴(kuò)大感受野的同時(shí)減少參數(shù)量的增長(zhǎng)。在經(jīng)過(guò)兩個(gè)卷積變換之后得到兩個(gè)特征圖:

    UA=FA(X)

    (7)

    UB=FB(X)

    (8)

    其中,F為對(duì)應(yīng)的卷積變換函數(shù)。然后,將得到的兩個(gè)卷積核大小不同(其中一個(gè)為寬卷積核)的特征圖UA和UB拼接進(jìn)行信息融合,得到的融合特征表示為

    U=Concat(UA,UB)

    (9)

    類(lèi)似于SE模塊[23]的結(jié)構(gòu),將得到的融合特征通過(guò)全局平均池化(global average pooling,GAP)操作在特征圖上進(jìn)行低維特征融合,以保留特征圖的主要特征,同時(shí)減小網(wǎng)絡(luò)的計(jì)算量,防止網(wǎng)絡(luò)的過(guò)擬合。該過(guò)程的計(jì)算公式為

    (10)

    由于輸入為一維數(shù)據(jù),則高度H為1。

    然后再通過(guò)一個(gè)全連接層(fully connected-layers,FC)將其以r的壓縮率壓縮,使其成為具有全局感受野的實(shí)數(shù)向量。向量的長(zhǎng)度為C/r,S為GAP操作后的輸出,壓縮操作后的輸出為

    (11)

    式中,B為批量正則化層;W為全連接層權(quán)重。

    在壓縮后,通過(guò)兩個(gè)全連接層分別將壓縮的輸出擴(kuò)張為原始維度(ZA和ZB),然后通過(guò)Softmax操作進(jìn)行權(quán)值歸一化,得到不同卷積核和不同通道的權(quán)值A(chǔ)和B:

    (12)

    (13)

    根據(jù)這樣兩個(gè)不同的權(quán)重值,可自適應(yīng)地選擇不同空間尺度上的信息。將歸一化之后的不同的權(quán)重乘以對(duì)應(yīng)的特征圖,可得到最終加權(quán)之后由網(wǎng)絡(luò)自適應(yīng)選擇的輸出:

    Y=AUA+BUB=A′+B′

    (14)

    本文所提方法中的每個(gè)殘差塊的卷積層都使用了可選擇的寬卷積核,使得網(wǎng)絡(luò)能夠在提取特征的每一個(gè)階段都能夠自適應(yīng)地選擇合適的卷積核進(jìn)行卷積,并且通過(guò)寬卷積核的壓縮擴(kuò)張?zhí)幚磉_(dá)到信號(hào)降噪的功能。此外,還修改了網(wǎng)絡(luò)中的通道數(shù)為8、16、32、64,以減小網(wǎng)絡(luò)的參數(shù)量,提高網(wǎng)絡(luò)的訓(xùn)練效率。

    2.2 所提網(wǎng)絡(luò)模型的結(jié)構(gòu)

    圖4展示了所提模型的整體結(jié)構(gòu)。該模型以18層的殘差網(wǎng)絡(luò)為基礎(chǔ),對(duì)其中所有殘差塊中的卷積層進(jìn)行改進(jìn),整個(gè)網(wǎng)絡(luò)由第一層的卷積、16個(gè)改進(jìn)的卷積層和一個(gè)全連接層組成,其具體結(jié)構(gòu)參數(shù)如表1所示。

    圖4 所提方法網(wǎng)絡(luò)模型結(jié)構(gòu)Fig.4 Structure diagram of the proposed method network

    每個(gè)改進(jìn)卷積層的結(jié)構(gòu)如圖3所示。在數(shù)據(jù)輸入網(wǎng)絡(luò)時(shí),首先對(duì)數(shù)據(jù)進(jìn)行卷積操作,再輸入到后續(xù)網(wǎng)絡(luò)層中。后續(xù)網(wǎng)絡(luò)層按照殘差連接的方式分為兩種不同的殘差塊。在輸入輸出特征通道不改變時(shí),使用殘差塊1依次經(jīng)過(guò)卷積層、批量正則化、激活函數(shù)、卷積層和批量正則化后直接與原始輸入特征進(jìn)行連接。在輸入輸出通道改變時(shí),使用殘差塊2經(jīng)過(guò)幾個(gè)與殘差塊1相同的操作之后,將輸入特征進(jìn)行卷積,以匹配改變的通道數(shù)。最后得到提取的特征,輸入全局平均池化層進(jìn)行特征信息的融合,經(jīng)過(guò)全連接層輸出最后分類(lèi)結(jié)果。

    2.3 所提方法故障診斷流程

    本文所提方法的軸承故障診斷具體流程如圖5所示。軸承故障診斷可大致分為4個(gè)步驟:①利用加速度傳感器采集軸承試驗(yàn)臺(tái)數(shù)據(jù),為數(shù)據(jù)集添加不同信噪比的高斯白噪聲以模擬噪聲輸入,并劃分構(gòu)建相應(yīng)的訓(xùn)練集和測(cè)試集;②構(gòu)建DWResNet,并將所有訓(xùn)練樣本進(jìn)行預(yù)處理,按批量輸入到網(wǎng)絡(luò)中進(jìn)行模型訓(xùn)練;③在訓(xùn)練的過(guò)程中,通過(guò)交叉熵?fù)p失函數(shù)進(jìn)行反向傳播,對(duì)網(wǎng)絡(luò)權(quán)重參數(shù)進(jìn)行優(yōu)化迭代;④將測(cè)試數(shù)據(jù)集輸入訓(xùn)練完成的網(wǎng)絡(luò)中,輸出滾動(dòng)軸承故障診斷結(jié)果。

    圖5 所提方法的故障診斷流程Fig.5 Fault diagnosis process of the proposed method

    3 實(shí)驗(yàn)驗(yàn)證

    3.1 實(shí)驗(yàn)數(shù)據(jù)集介紹

    3.1.1實(shí)驗(yàn)一

    本文采用德國(guó)帕德博恩大學(xué)(PBU)軸承數(shù)據(jù)集[24]對(duì)所提方法進(jìn)行驗(yàn)證。該試驗(yàn)臺(tái)的詳細(xì)設(shè)置請(qǐng)參考文獻(xiàn)[24]。試驗(yàn)臺(tái)設(shè)置如圖6所示,所有測(cè)試軸承的型號(hào)均為6203,采樣頻率為64 kHz。

    圖6 帕德博恩大學(xué)軸承試驗(yàn)臺(tái)[24]Fig.6 The test rig of PBU[24]

    本實(shí)驗(yàn)中使用的數(shù)據(jù)集的具體工況設(shè)置如表2所示,本數(shù)據(jù)集包含了5種故障,每個(gè)故障類(lèi)型下的數(shù)據(jù)都處理為表3所示數(shù)據(jù)。當(dāng)軸承的損傷深度為D∈(0,2] mm時(shí),設(shè)置為損傷等級(jí)1;當(dāng)軸承的損傷深度D∈(2.0,4.5] mm時(shí),設(shè)置為損傷等級(jí)2。

    表2 帕德博恩大學(xué)軸承數(shù)據(jù)所選工況參數(shù)

    表3 帕德博恩大學(xué)軸承數(shù)據(jù)詳細(xì)信息

    3.1.2實(shí)驗(yàn)二

    為驗(yàn)證本文方法的有效性,本文還采用了動(dòng)力傳動(dòng)故障診斷綜合實(shí)驗(yàn)臺(tái)(drivetrain diagnostic simulator,DDS)[25]所采集的數(shù)據(jù)進(jìn)行實(shí)驗(yàn)驗(yàn)證。實(shí)驗(yàn)臺(tái)設(shè)置如圖7所示,在平行齒輪箱軸承中,設(shè)定電機(jī)固定頻率為20 Hz,采樣頻率為12.8 kHz。每類(lèi)故障工況重復(fù)實(shí)驗(yàn)4次,每次采集各方向的振動(dòng)加速度信號(hào)共204 800個(gè)點(diǎn)。本數(shù)據(jù)集包含了內(nèi)圈、外圈、滾動(dòng)體故障和復(fù)合故障等4種故障類(lèi)型。

    圖7 動(dòng)力傳動(dòng)故障診斷綜合實(shí)驗(yàn)臺(tái)[25]Fig.7 Power transmission fault diagnosis comprehensive experimental bench[25]

    3.2 實(shí)驗(yàn)環(huán)境及模型參數(shù)設(shè)置

    本實(shí)驗(yàn)中,樣本長(zhǎng)度設(shè)置為4096,通過(guò)滑窗的方式截取樣本。在實(shí)驗(yàn)一中,滑窗步長(zhǎng)與樣本長(zhǎng)度相等,因此樣本之間無(wú)重疊。其中,每類(lèi)樣本1200條,并按1∶1的比例劃分訓(xùn)練集和測(cè)試集,詳細(xì)信息見(jiàn)表3。實(shí)驗(yàn)二中,樣本長(zhǎng)度也設(shè)置為4096,通過(guò)滑窗截取了每類(lèi)樣本600條,滑窗步長(zhǎng)為582,并按1∶1的比例劃分訓(xùn)練集和測(cè)試集。

    分別在信噪比(signal to noise ratio,SNR)為-4,-2,0,2,4 dB的高斯白噪聲以及無(wú)人為添加噪聲情況下,對(duì)卷積神經(jīng)網(wǎng)絡(luò)(CNN)[26]、VGG-16[27]、首層寬卷積深度神經(jīng)網(wǎng)絡(luò)(WDCNN)[28]、ResNet、首層寬卷積的殘差網(wǎng)絡(luò)(WResNet)和本文所提方法進(jìn)行了對(duì)比實(shí)驗(yàn),每組實(shí)驗(yàn)的超參數(shù)設(shè)置相同,但在實(shí)驗(yàn)一、實(shí)驗(yàn)二中設(shè)置初始學(xué)習(xí)率分別為0.01和0.1。并且為保證訓(xùn)練精度和損失函數(shù)收斂,在使用CNN和VGG-16方法時(shí)設(shè)置學(xué)習(xí)率為固定值,詳情見(jiàn)表4,每組實(shí)驗(yàn)重復(fù)10次以保證實(shí)驗(yàn)結(jié)果的準(zhǔn)確性。同時(shí),針對(duì)不同大小的寬卷積核,設(shè)置了空洞率d分別為1、2、14、26和38。信噪比的計(jì)算公式如下:

    (15)

    式中,PSignal、PNoise分別為原始信號(hào)和噪聲的能量。

    以正常狀況為例,經(jīng)過(guò)加噪后的部分信號(hào)如圖8所示。實(shí)驗(yàn)運(yùn)行環(huán)境為Windows系統(tǒng)和Intel(R) Xeon(R) Silver 4210R CPU @2.40GHz處理器,使用了基于Python(3.9)的Pytorch 深度學(xué)習(xí)框架構(gòu)建所提方法和對(duì)比方法模型,并采用NVIDIA GeForce RTX 3090 GPU進(jìn)行訓(xùn)練加速。

    (a)無(wú)人為添加噪聲 (b)RSNR=4 dB

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

    3.3.1實(shí)驗(yàn)一結(jié)果分析

    在實(shí)驗(yàn)一中,根據(jù)實(shí)驗(yàn)結(jié)果得到表5(表中“無(wú)”表示無(wú)人為添加噪聲)。由表5可知,雖然在沒(méi)有人為添加噪聲的情況下,不同方法均能保持較高的診斷準(zhǔn)確率,但是隨著噪聲干擾程度的增加,不同方法的診斷準(zhǔn)確率也隨之降低,說(shuō)明了強(qiáng)噪聲干擾對(duì)模型的抗噪性要求較高,例如CNN的診斷準(zhǔn)確率在RSNR=-4 dB的情況下比無(wú)噪聲情況下降低了11.29個(gè)百分點(diǎn)。區(qū)別于其他的方法,本文所提方法僅降低了1.99個(gè)百分點(diǎn),在各個(gè)不同噪聲情況下,所提方法能夠準(zhǔn)確地提取出滾動(dòng)軸承的振動(dòng)信號(hào)特征,具有更高的穩(wěn)定性。

    表5 實(shí)驗(yàn)一診斷準(zhǔn)確率結(jié)果

    對(duì)比CNN、VGG-16和ResNet在同一噪聲下的診斷精度可以看出,網(wǎng)絡(luò)深度的增加對(duì)滾動(dòng)軸承故障特征的提取能力和穩(wěn)定性均有一定提高。當(dāng)然,層數(shù)并非越深越好,超過(guò)一定數(shù)量后將造成梯度爆炸或梯度消失等問(wèn)題。對(duì)比CNN、WDCNN和WResNet可以發(fā)現(xiàn),卷積核寬度的增大對(duì)網(wǎng)絡(luò)容噪具有較好的效果。同樣地,需要選擇適當(dāng)?shù)木W(wǎng)絡(luò)寬度和網(wǎng)絡(luò)深度,使得方法能夠在軸承故障診斷上有更有效的結(jié)果。

    由對(duì)比實(shí)驗(yàn)結(jié)果可知,本文提出方法具有較優(yōu)的容噪能力,能夠在噪聲環(huán)境下?lián)碛酗@著提高的滾動(dòng)軸承故障分類(lèi)能力和模型的穩(wěn)定性。本文所提方法在不同噪聲水平情況下均有較好的診斷準(zhǔn)確率,特別是在RSNR=-4 dB的強(qiáng)噪聲干擾下,相較于CNN、WDCNN、VGG-16、ResNet和WResNet,診斷準(zhǔn)確率分別提高了17.88、3.65、2.16、2.05和3.35個(gè)百分點(diǎn)。同時(shí),本文所提方法在不同噪聲干擾下的診斷標(biāo)準(zhǔn)差相較于其他方法均較小,表示它在強(qiáng)噪聲下具有較高的模型穩(wěn)定性。

    為了探究強(qiáng)噪聲信號(hào)對(duì)模型分類(lèi)邊界形成的影響,以RSNR=-4 dB為例,繪制了本文所提方法在網(wǎng)絡(luò)不同階段所提取故障特征的t-SNE可視化聚類(lèi)圖,如圖9所示。由圖9不難看出,輸入的強(qiáng)噪聲原始信號(hào)呈現(xiàn)極為嚴(yán)重的混淆、重疊現(xiàn)象,導(dǎo)致難以直接聚類(lèi)辨識(shí)。隨著數(shù)據(jù)經(jīng)過(guò)更深層的階段2、3、4,網(wǎng)絡(luò)對(duì)不同健康狀況的分類(lèi)邊界逐漸清晰。階段4時(shí),同類(lèi)故障特征的類(lèi)內(nèi)聚集程度顯著提高,模型的分類(lèi)邊界也最為明顯。這說(shuō)明經(jīng)過(guò)多個(gè)階段的模型訓(xùn)練,所提方法能夠有效提取強(qiáng)噪聲下不同類(lèi)別的故障特征,較小冗余噪聲的干擾,從而提高對(duì)目標(biāo)軸承的診斷精度。

    (a)輸入 (b)階段2輸出

    同時(shí),為進(jìn)一步展示不同方法對(duì)各個(gè)健康狀況的詳細(xì)辨識(shí)能力,圖10展示了RSNR=-4 dB下不同方法辨識(shí)結(jié)果的混淆矩陣。圖10顯示,不同方法在辨識(shí)分類(lèi)的過(guò)程中對(duì)NC和OR2兩類(lèi)健康狀況是較為容易辨識(shí)的,但是對(duì)OR1、IR1和IR2三類(lèi)故障難以有效辨別,分析其可能的原因是外圈故障和內(nèi)圈故障特征相似,且兩個(gè)不同故障等級(jí)的內(nèi)圈故障特征相似。本文所提出的方法雖然也存在上述情況,但是相較于其他方法,本文方法對(duì)難診斷類(lèi)別的診斷準(zhǔn)確率最高,說(shuō)明了所提方法能夠有效提高不同故障類(lèi)型中判別性特征的辨識(shí)能力,從而提高診斷準(zhǔn)確率。

    (a)CNN (b)WDCNN

    為了進(jìn)一步分析不同空洞率對(duì)模型診斷效果的影響,以強(qiáng)噪聲RSNR=-4 dB為例,分別設(shè)置空洞率d為1、2、14、26和38,并進(jìn)行相應(yīng)的故障診斷實(shí)驗(yàn),結(jié)果如圖11所示。由圖11不難看出,所提方法在d取14時(shí),得到了高達(dá)98%的診斷準(zhǔn)確率,這說(shuō)明本文所選超參數(shù)較為合適,能使模型具有較高的噪聲魯棒性。由于空洞率的大小直接決定了感受野的大小,這也印證了單一地提高空洞率對(duì)噪聲工況下軸承故障診斷準(zhǔn)確率的提高并不是絕對(duì)的,需要選擇合適的空洞率對(duì)數(shù)據(jù)進(jìn)行更為全面的特征提取,以獲得更豐富的分類(lèi)辨識(shí)信息,進(jìn)而提高軸承故障診斷的準(zhǔn)確率。

    圖11 RSNR=-4 dB下不同空洞率下所提方法診斷準(zhǔn)確率Fig.11 Diagnostic accuracy of the proposed method at RSNR=-4dB for different d

    3.3.2實(shí)驗(yàn)二結(jié)果分析

    根據(jù)實(shí)驗(yàn)二實(shí)驗(yàn)結(jié)果繪制相應(yīng)的精度直方圖,見(jiàn)圖12,可以看出,在各個(gè)不同等級(jí)的RSNR下,本文方法相較于其他方法精度均有一定程度的提高,并且具有較好的穩(wěn)定性。RSNR=-4 dB時(shí)的診斷準(zhǔn)確率相較于無(wú)添加噪聲時(shí),在各個(gè)方法下分別降低16.26、8.16、16.91、8.84、9.53和5.39個(gè)百分點(diǎn),由此可見(jiàn)本文所提方法的診斷效果受到噪聲的影響最小,印證了本文所提方法對(duì)強(qiáng)噪聲的魯棒性。并且本文所提方法在RSNR=-4 dB時(shí)的標(biāo)準(zhǔn)差絕對(duì)值相較于其他方法最小,驗(yàn)證了本文所提方法的穩(wěn)定性。

    圖12 實(shí)驗(yàn)二結(jié)果直方圖Fig.12 Histograms of experiment 2 results

    4 結(jié)論

    針對(duì)滾動(dòng)軸承在強(qiáng)噪聲干擾下故障診斷提取特征能力較差、辨識(shí)能力不足的問(wèn)題,提出了一種基于動(dòng)態(tài)寬卷積殘差網(wǎng)絡(luò)的軸承故障診斷方法。該方法通過(guò)可選擇的寬卷積層對(duì)不同噪聲干擾下的滾動(dòng)軸承監(jiān)測(cè)信號(hào)進(jìn)行故障特征提取,通過(guò)生成的自適應(yīng)權(quán)重自動(dòng)地為網(wǎng)絡(luò)選擇使用卷積核,在網(wǎng)絡(luò)深度加深的同時(shí)保證了網(wǎng)絡(luò)的抗噪能力。

    實(shí)驗(yàn)結(jié)果表明,相較于CNN、WDCNN、VGG-16、ResNet和WResNet等網(wǎng)絡(luò),本文所提方法具有更高的診斷準(zhǔn)確率。在無(wú)人為添加噪聲情況下,本文所提方法在實(shí)驗(yàn)一中能夠取得99.99%的診斷準(zhǔn)確率;而在信噪比為-4,-2,0,2,4 dB下,該方法依然可取得99.28%的平均準(zhǔn)確率。并且實(shí)驗(yàn)二的結(jié)果也印證了本文所提方法的有效性。因此,本文所提方法具有較高的辨識(shí)能力、抗噪能力和泛化能力,可為實(shí)際工業(yè)強(qiáng)噪聲干擾條件下的軸承故障診斷提供一定的參考借鑒。

    猜你喜歡
    殘差故障診斷軸承
    軸承知識(shí)
    哈爾濱軸承(2022年2期)2022-07-22 06:39:32
    基于雙向GRU與殘差擬合的車(chē)輛跟馳建模
    軸承知識(shí)
    哈爾濱軸承(2022年1期)2022-05-23 13:13:24
    軸承知識(shí)
    哈爾濱軸承(2021年2期)2021-08-12 06:11:46
    軸承知識(shí)
    哈爾濱軸承(2021年1期)2021-07-21 05:43:16
    基于殘差學(xué)習(xí)的自適應(yīng)無(wú)人機(jī)目標(biāo)跟蹤算法
    基于遞歸殘差網(wǎng)絡(luò)的圖像超分辨率重建
    因果圖定性分析法及其在故障診斷中的應(yīng)用
    平穩(wěn)自相關(guān)過(guò)程的殘差累積和控制圖
    河南科技(2015年8期)2015-03-11 16:23:52
    基于LCD和排列熵的滾動(dòng)軸承故障診斷
    免费黄网站久久成人精品| 成年女人看的毛片在线观看| 狠狠狠狠99中文字幕| 尾随美女入室| 国产成人影院久久av| 最近的中文字幕免费完整| 伦理电影大哥的女人| 国产 一区 欧美 日韩| 我要搜黄色片| 在线国产一区二区在线| 国产日本99.免费观看| 亚州av有码| 日韩精品青青久久久久久| 不卡视频在线观看欧美| 国产成人aa在线观看| 狠狠狠狠99中文字幕| 又粗又硬又长又爽又黄的视频 | 亚洲经典国产精华液单| 老师上课跳d突然被开到最大视频| 麻豆精品久久久久久蜜桃| 日韩av在线大香蕉| 日韩欧美一区二区三区在线观看| 久久久午夜欧美精品| 免费看光身美女| 欧美高清性xxxxhd video| 99久久无色码亚洲精品果冻| 国产精品爽爽va在线观看网站| 亚洲成人av在线免费| 亚洲av中文字字幕乱码综合| 国产亚洲av片在线观看秒播厂 | 日韩在线高清观看一区二区三区| 老师上课跳d突然被开到最大视频| 九九爱精品视频在线观看| 久久精品国产99精品国产亚洲性色| 简卡轻食公司| 婷婷六月久久综合丁香| 神马国产精品三级电影在线观看| 亚洲成人久久性| 插逼视频在线观看| 日本熟妇午夜| 久久久欧美国产精品| 老师上课跳d突然被开到最大视频| a级毛色黄片| 简卡轻食公司| 午夜免费激情av| 精品久久久久久久末码| 美女xxoo啪啪120秒动态图| 国产一区二区在线av高清观看| 丰满的人妻完整版| 男女做爰动态图高潮gif福利片| 国产精品福利在线免费观看| 欧美一区二区精品小视频在线| 欧美日韩综合久久久久久| 亚洲在久久综合| 天堂网av新在线| 高清在线视频一区二区三区 | 九九热线精品视视频播放| 丰满的人妻完整版| eeuss影院久久| 国产极品天堂在线| 1000部很黄的大片| 国产高潮美女av| ponron亚洲| 亚洲精品亚洲一区二区| 久久精品人妻少妇| 网址你懂的国产日韩在线| 国产日韩欧美在线精品| 久久久色成人| 六月丁香七月| 国产成人a∨麻豆精品| 亚洲熟妇中文字幕五十中出| 波多野结衣高清无吗| 美女被艹到高潮喷水动态| av专区在线播放| 亚洲久久久久久中文字幕| 麻豆成人午夜福利视频| 久久午夜亚洲精品久久| 国模一区二区三区四区视频| 青青草视频在线视频观看| 哪里可以看免费的av片| 高清日韩中文字幕在线| 久久人人精品亚洲av| 神马国产精品三级电影在线观看| 在线观看av片永久免费下载| 免费观看在线日韩| 女同久久另类99精品国产91| 黄色配什么色好看| 久久久久国产网址| 亚洲欧美日韩东京热| 亚洲成av人片在线播放无| av专区在线播放| 免费观看的影片在线观看| 一个人免费在线观看电影| 久久精品综合一区二区三区| 一边亲一边摸免费视频| 久久人人爽人人爽人人片va| 久久久精品94久久精品| 久久精品国产清高在天天线| 男女做爰动态图高潮gif福利片| 亚洲成人精品中文字幕电影| 国产三级中文精品| 国产乱人视频| 一进一出抽搐gif免费好疼| 亚洲一区二区三区色噜噜| 亚洲第一区二区三区不卡| 99视频精品全部免费 在线| 国产精品99久久久久久久久| av在线亚洲专区| 小蜜桃在线观看免费完整版高清| 国产精品不卡视频一区二区| 国产成人午夜福利电影在线观看| 亚洲无线观看免费| 久久鲁丝午夜福利片| 成年版毛片免费区| 欧美+日韩+精品| 久久久久国产网址| 色综合色国产| 免费av不卡在线播放| 夫妻性生交免费视频一级片| 国产熟女欧美一区二区| 国产亚洲精品久久久久久毛片| 国产午夜福利久久久久久| 1000部很黄的大片| 一进一出抽搐gif免费好疼| 亚洲成人久久性| 日本黄色视频三级网站网址| 亚洲中文字幕一区二区三区有码在线看| 一级黄色大片毛片| 亚洲一区二区三区色噜噜| 男的添女的下面高潮视频| 国产三级中文精品| 国产av在哪里看| 97人妻精品一区二区三区麻豆| 国产三级中文精品| av国产免费在线观看| 久久99热这里只有精品18| 一边亲一边摸免费视频| 国产精品爽爽va在线观看网站| 乱人视频在线观看| 国产精品日韩av在线免费观看| 中文字幕熟女人妻在线| 一边亲一边摸免费视频| 综合色av麻豆| 美女xxoo啪啪120秒动态图| 日日干狠狠操夜夜爽| 亚洲av熟女| 亚洲aⅴ乱码一区二区在线播放| 欧美在线一区亚洲| 亚洲最大成人手机在线| 国产老妇伦熟女老妇高清| 毛片女人毛片| 久久久精品欧美日韩精品| 天天一区二区日本电影三级| 99国产极品粉嫩在线观看| 国产一区二区在线观看日韩| 热99re8久久精品国产| 卡戴珊不雅视频在线播放| 美女内射精品一级片tv| 欧美成人一区二区免费高清观看| 亚洲欧美日韩卡通动漫| 18禁在线无遮挡免费观看视频| 毛片女人毛片| 久久午夜亚洲精品久久| 亚洲成人久久爱视频| 啦啦啦啦在线视频资源| 51国产日韩欧美| 人妻久久中文字幕网| 成年av动漫网址| 成人一区二区视频在线观看| 久久久久久久亚洲中文字幕| 亚洲人成网站在线观看播放| 中文资源天堂在线| 国产高清不卡午夜福利| 欧美三级亚洲精品| 日日啪夜夜撸| av国产免费在线观看| 婷婷六月久久综合丁香| 99热网站在线观看| 五月玫瑰六月丁香| 一级二级三级毛片免费看| av视频在线观看入口| 一个人观看的视频www高清免费观看| 男人的好看免费观看在线视频| 在线免费观看不下载黄p国产| 亚洲久久久久久中文字幕| av在线老鸭窝| 国产成人一区二区在线| 一个人看的www免费观看视频| av在线观看视频网站免费| 国产蜜桃级精品一区二区三区| 乱人视频在线观看| 久久午夜福利片| 日日摸夜夜添夜夜爱| 在线观看午夜福利视频| 久久九九热精品免费| 免费看日本二区| 国产av麻豆久久久久久久| 校园春色视频在线观看| 最近最新中文字幕大全电影3| 午夜视频国产福利| 黄色一级大片看看| 欧美最新免费一区二区三区| 不卡一级毛片| 哪个播放器可以免费观看大片| 亚洲欧美精品专区久久| 成人国产麻豆网| 亚洲国产欧美人成| 国模一区二区三区四区视频| 99热6这里只有精品| 级片在线观看| 天堂av国产一区二区熟女人妻| av天堂中文字幕网| 国产极品精品免费视频能看的| 又粗又爽又猛毛片免费看| 村上凉子中文字幕在线| 国产精品免费一区二区三区在线| 国产一区二区激情短视频| av在线天堂中文字幕| 听说在线观看完整版免费高清| 精品午夜福利在线看| 亚洲国产欧美在线一区| 国产一区二区三区在线臀色熟女| 色播亚洲综合网| 狂野欧美白嫩少妇大欣赏| 午夜久久久久精精品| 成人午夜精彩视频在线观看| 伦理电影大哥的女人| 欧美区成人在线视频| 中文字幕av成人在线电影| 国产高清激情床上av| 久久久精品欧美日韩精品| 久久午夜亚洲精品久久| 性欧美人与动物交配| 99久久精品一区二区三区| 我的老师免费观看完整版| 午夜a级毛片| 国产视频首页在线观看| 亚洲va在线va天堂va国产| 久久精品久久久久久久性| 免费观看精品视频网站| 在线a可以看的网站| 国产三级中文精品| 国产精品,欧美在线| 久久久精品欧美日韩精品| 能在线免费看毛片的网站| 国产精品一区www在线观看| 亚洲美女搞黄在线观看| 亚洲丝袜综合中文字幕| 国产黄片美女视频| 大型黄色视频在线免费观看| 18禁黄网站禁片免费观看直播| 国产一区二区在线av高清观看| 欧美成人a在线观看| 日本撒尿小便嘘嘘汇集6| 亚洲久久久久久中文字幕| 亚洲国产日韩欧美精品在线观看| 国产熟女欧美一区二区| 国产成人91sexporn| 男女视频在线观看网站免费| 日韩一区二区三区影片| 中文字幕av在线有码专区| 成人毛片60女人毛片免费| 日韩欧美三级三区| 免费观看人在逋| 三级毛片av免费| 一级毛片aaaaaa免费看小| 欧美日韩综合久久久久久| 免费电影在线观看免费观看| 日日摸夜夜添夜夜爱| 国产精品蜜桃在线观看 | 人妻制服诱惑在线中文字幕| 亚洲国产高清在线一区二区三| 五月伊人婷婷丁香| 国产一区二区三区在线臀色熟女| 午夜爱爱视频在线播放| 国产在线精品亚洲第一网站| 国内精品一区二区在线观看| 搞女人的毛片| 成人毛片a级毛片在线播放| 搡女人真爽免费视频火全软件| 长腿黑丝高跟| 国国产精品蜜臀av免费| 亚洲欧美日韩卡通动漫| 国产精华一区二区三区| 一本精品99久久精品77| 免费看光身美女| 亚洲精品久久久久久婷婷小说 | 天堂中文最新版在线下载 | 国内精品美女久久久久久| 日韩制服骚丝袜av| 久久综合国产亚洲精品| 少妇熟女欧美另类| 小蜜桃在线观看免费完整版高清| 欧美一区二区精品小视频在线| 亚洲国产欧洲综合997久久,| 18禁裸乳无遮挡免费网站照片| 一卡2卡三卡四卡精品乱码亚洲| 亚洲最大成人中文| av国产免费在线观看| 嫩草影院新地址| 亚洲成a人片在线一区二区| 成人午夜精彩视频在线观看| 亚洲最大成人av| 久久久色成人| 观看免费一级毛片| 身体一侧抽搐| 亚洲av成人av| 国产美女午夜福利| 国产激情偷乱视频一区二区| 22中文网久久字幕| 丰满乱子伦码专区| 在现免费观看毛片| 天堂网av新在线| 国产69精品久久久久777片| 国产高清三级在线| 欧美日韩一区二区视频在线观看视频在线 | 亚洲欧洲国产日韩| 成人特级av手机在线观看| 久久国内精品自在自线图片| 日本在线视频免费播放| 精品少妇黑人巨大在线播放 | 不卡一级毛片| 亚洲国产欧洲综合997久久,| 日韩欧美精品免费久久| 免费黄网站久久成人精品| 国产精品国产高清国产av| 亚洲欧美日韩无卡精品| 国产一区二区三区在线臀色熟女| 狂野欧美白嫩少妇大欣赏| 亚洲人成网站在线观看播放| 日韩精品青青久久久久久| 91久久精品国产一区二区三区| 少妇熟女欧美另类| 国产伦理片在线播放av一区 | 99在线人妻在线中文字幕| 一本一本综合久久| 久久久a久久爽久久v久久| 在线a可以看的网站| 亚洲av.av天堂| 国产亚洲av嫩草精品影院| 国产久久久一区二区三区| 大香蕉久久网| 欧美日本亚洲视频在线播放| 青青草视频在线视频观看| 国产亚洲欧美98| 国产探花极品一区二区| 亚洲国产精品sss在线观看| 亚洲七黄色美女视频| 桃色一区二区三区在线观看| 欧美激情国产日韩精品一区| 久久久午夜欧美精品| 深爱激情五月婷婷| 一级av片app| 国产毛片a区久久久久| 国产精品久久久久久精品电影| 一本久久中文字幕| 亚洲天堂国产精品一区在线| 久久精品91蜜桃| 麻豆国产97在线/欧美| 国产黄色视频一区二区在线观看 | 日本一本二区三区精品| 两个人的视频大全免费| 哪里可以看免费的av片| 在现免费观看毛片| 可以在线观看毛片的网站| 蜜臀久久99精品久久宅男| 岛国在线免费视频观看| 国产精品麻豆人妻色哟哟久久 | 国产精品一区www在线观看| 国产av麻豆久久久久久久| 男人舔奶头视频| 成人午夜高清在线视频| av在线老鸭窝| 亚洲激情五月婷婷啪啪| 91在线精品国自产拍蜜月| 又爽又黄无遮挡网站| 看十八女毛片水多多多| 日韩中字成人| 天天躁日日操中文字幕| 一本精品99久久精品77| 成年女人看的毛片在线观看| 激情 狠狠 欧美| 午夜精品国产一区二区电影 | 日韩大尺度精品在线看网址| av在线蜜桃| 最近视频中文字幕2019在线8| 九九爱精品视频在线观看| 美女脱内裤让男人舔精品视频 | 日日干狠狠操夜夜爽| 欧洲精品卡2卡3卡4卡5卡区| 国产熟女欧美一区二区| 男女边吃奶边做爰视频| 国产美女午夜福利| 麻豆成人av视频| 久久精品夜夜夜夜夜久久蜜豆| 午夜免费激情av| 国产成人午夜福利电影在线观看| 中文欧美无线码| 校园人妻丝袜中文字幕| 别揉我奶头 嗯啊视频| 日本免费一区二区三区高清不卡| 青青草视频在线视频观看| 国产91av在线免费观看| 身体一侧抽搐| 男女边吃奶边做爰视频| 啦啦啦啦在线视频资源| 欧美日本视频| 国产 一区精品| 久久久久网色| 麻豆成人av视频| 国产精品国产三级国产av玫瑰| 久久久久九九精品影院| 国产精品人妻久久久久久| av国产免费在线观看| 最近手机中文字幕大全| 小蜜桃在线观看免费完整版高清| 97人妻精品一区二区三区麻豆| 久久欧美精品欧美久久欧美| .国产精品久久| 欧美性猛交黑人性爽| 晚上一个人看的免费电影| 日韩一本色道免费dvd| 午夜老司机福利剧场| 亚洲人成网站在线观看播放| 免费电影在线观看免费观看| 成人鲁丝片一二三区免费| 边亲边吃奶的免费视频| 97热精品久久久久久| 不卡视频在线观看欧美| 亚洲精品国产成人久久av| 国产在线精品亚洲第一网站| 亚洲精品日韩av片在线观看| 最近的中文字幕免费完整| 亚洲中文字幕日韩| 91在线精品国自产拍蜜月| 赤兔流量卡办理| 久久久久久久久久黄片| 人妻久久中文字幕网| 久久99精品国语久久久| 亚洲国产色片| 免费观看在线日韩| 日韩成人av中文字幕在线观看| 成人欧美大片| 日韩三级伦理在线观看| 国产久久久一区二区三区| 色噜噜av男人的天堂激情| 精品久久久久久久久久久久久| 亚洲精品成人久久久久久| 九九爱精品视频在线观看| 97超视频在线观看视频| 波多野结衣巨乳人妻| 超碰av人人做人人爽久久| 白带黄色成豆腐渣| 精品欧美国产一区二区三| 亚洲av中文av极速乱| 国产v大片淫在线免费观看| 丰满的人妻完整版| 99久久精品一区二区三区| 久久综合国产亚洲精品| 国产高清不卡午夜福利| 国产美女午夜福利| 欧美日韩精品成人综合77777| 精品人妻一区二区三区麻豆| 国产一区二区亚洲精品在线观看| 国产伦在线观看视频一区| 亚洲五月天丁香| 国产极品精品免费视频能看的| 非洲黑人性xxxx精品又粗又长| 免费观看a级毛片全部| 亚洲国产精品国产精品| 国产一区二区激情短视频| 亚洲在久久综合| 一边摸一边抽搐一进一小说| 日韩人妻高清精品专区| 人妻久久中文字幕网| 免费不卡的大黄色大毛片视频在线观看 | 亚洲av第一区精品v没综合| 在线观看66精品国产| 国产精品久久久久久精品电影小说 | 国模一区二区三区四区视频| 国内少妇人妻偷人精品xxx网站| 女人被狂操c到高潮| 国内揄拍国产精品人妻在线| 国产精品麻豆人妻色哟哟久久 | 亚洲国产欧美在线一区| 欧美日韩国产亚洲二区| 18禁裸乳无遮挡免费网站照片| 精品人妻偷拍中文字幕| 亚洲av电影不卡..在线观看| 国产成人精品婷婷| 国产探花极品一区二区| 给我免费播放毛片高清在线观看| 中文字幕av成人在线电影| 日本黄色片子视频| 国产蜜桃级精品一区二区三区| 少妇丰满av| 99热网站在线观看| 午夜视频国产福利| 毛片一级片免费看久久久久| 一区二区三区高清视频在线| 特大巨黑吊av在线直播| 爱豆传媒免费全集在线观看| 色5月婷婷丁香| 亚洲精品久久国产高清桃花| 男人舔奶头视频| 欧美高清性xxxxhd video| АⅤ资源中文在线天堂| 久久鲁丝午夜福利片| 久久九九热精品免费| 国产精品福利在线免费观看| 午夜久久久久精精品| 淫秽高清视频在线观看| 欧美最新免费一区二区三区| 午夜激情欧美在线| 精品一区二区免费观看| 欧美高清成人免费视频www| 亚洲国产精品成人综合色| 边亲边吃奶的免费视频| 午夜精品国产一区二区电影 | 国产中年淑女户外野战色| 一区福利在线观看| 麻豆一二三区av精品| 卡戴珊不雅视频在线播放| 91久久精品电影网| 国产精品久久电影中文字幕| 国产黄a三级三级三级人| 久久久久国产网址| 成人av在线播放网站| 五月伊人婷婷丁香| 国产精品电影一区二区三区| 亚洲国产精品成人综合色| 国产精品综合久久久久久久免费| 中文字幕免费在线视频6| 麻豆国产97在线/欧美| 99久久九九国产精品国产免费| 成人性生交大片免费视频hd| 看非洲黑人一级黄片| 啦啦啦观看免费观看视频高清| 搡老妇女老女人老熟妇| 两个人的视频大全免费| 精品99又大又爽又粗少妇毛片| 麻豆国产97在线/欧美| 国产精品.久久久| 亚洲人与动物交配视频| 久久久成人免费电影| 在线播放无遮挡| 一本久久精品| 深夜a级毛片| 日韩视频在线欧美| 日日摸夜夜添夜夜爱| 久久热精品热| 高清在线视频一区二区三区 | 美女内射精品一级片tv| 亚洲av中文av极速乱| 天堂影院成人在线观看| av在线天堂中文字幕| 精品熟女少妇av免费看| 99久国产av精品| 少妇的逼好多水| 国语自产精品视频在线第100页| 午夜激情福利司机影院| 精品久久久久久久久亚洲| av专区在线播放| 精品久久久久久久末码| 最好的美女福利视频网| 亚洲在线观看片| 99国产精品一区二区蜜桃av| 在线观看av片永久免费下载| 国产精品永久免费网站| 日韩强制内射视频| 日本黄色视频三级网站网址| 国产精品蜜桃在线观看 | 真实男女啪啪啪动态图| 欧美日韩乱码在线| 一区二区三区高清视频在线| 色尼玛亚洲综合影院| 最近视频中文字幕2019在线8| 免费观看在线日韩| 中文字幕制服av| 欧美性猛交黑人性爽| 高清午夜精品一区二区三区 | 菩萨蛮人人尽说江南好唐韦庄 | www.av在线官网国产| 日韩 亚洲 欧美在线| 亚洲av.av天堂| 免费看光身美女| 日韩 亚洲 欧美在线| 91麻豆精品激情在线观看国产| 亚洲性久久影院| 国产精品久久久久久精品电影小说 | 久久这里只有精品中国| 在线播放无遮挡| 亚洲四区av| 免费看美女性在线毛片视频| 精品少妇黑人巨大在线播放 | 国内精品久久久久精免费| 亚洲aⅴ乱码一区二区在线播放| 又粗又爽又猛毛片免费看| 一本精品99久久精品77| 我要搜黄色片| 成人毛片60女人毛片免费| 久久精品国产亚洲av天美| 一级二级三级毛片免费看| av在线蜜桃| 别揉我奶头 嗯啊视频| 免费看美女性在线毛片视频| 悠悠久久av| 国产乱人偷精品视频| 一个人免费在线观看电影| 国产女主播在线喷水免费视频网站 | 搡老妇女老女人老熟妇| 天天躁夜夜躁狠狠久久av| 国产亚洲av片在线观看秒播厂 | 在线观看一区二区三区| 亚洲欧美成人精品一区二区|