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

    基于特征金字塔卷積循環(huán)神經網絡的故障診斷方法

    2022-03-08 02:37:48劉秀麗徐小力
    上海交通大學學報 2022年2期
    關鍵詞:支路金字塔故障診斷

    劉秀麗, 徐小力

    (北京信息科技大學 現代測控技術教育部重點實驗室,北京 100192)

    行星齒輪和軸承被廣泛應用于風力發(fā)電機組和航空設備中.由于長期工作在變轉速和變載荷的環(huán)境中,部件容易出現磨損等故障[1],對其開展故障診斷方法研究具有重要意義.

    隨著機器學習方法的日漸成熟,傳統(tǒng)智能診斷方法在故障診斷領域應用廣泛.基于機器學習的傳統(tǒng)故障診斷方法分為3個步驟:① 通過對信號的降噪和分解等實現對信號的預處理;② 利用特征提取方法獲得經處理后的信號的時域、頻域和其他特征;③ 將最終提取的特征輸入到基于機器學習的分類器中進行訓練,利用訓練分類器實現對故障類型的識別.在上述步驟中,特征提取方法包括小波變換[1-2]、同步提取變換[3-4]和變分模態(tài)分解[5-6]等,其能夠二次過濾信號的無用特征,使得所需故障特征更明顯.反向傳播神經網絡[7-8]和支持向量機[9-10]等方法已被應用于故障診斷.雖然傳統(tǒng)智能方法具有特征提取操作簡單和分類器參數容易調節(jié)等特點,最終的診斷識別率可滿足大部分場景要求,但其割裂了故障特征提取與診斷識別,在信號特征處理中需要大量的專家經驗.

    深度學習的目的為模仿人腦的學習方式,構建深層網絡模型.其利用大數據訓練模型獲得數據中隱含的特征,是一種集提取特征與模式識別為一體的智能模型.卷積神經網絡(CNN)[11-15]具有特征強大的學習能力和較高的泛化能力,在圖像識別領域應用廣泛.近幾年,一些學者提出將CNN應用于故障診斷領域.李恒等[16]對軸承滾動振動信號進行短時傅里葉變換,得到二維圖的時頻圖樣本,再將其分為訓練集和測試集放入二維CNN中進行故障診斷;韓濤等[17]提出一種基于多小波變換和CNN的滾動軸承智能復合故障診斷方法,將振動信號進行處理后得到多小波系數分支并利用其構造特征圖,輸入二維CNN以實現診斷.但上述研究中的CNN對故障信息提取的方式局限在二維形式,對一維時序原始信息的利用有限,影響故障診斷的準確率,并且對信號特征的處理方法仍依賴人工操作,沒有充分發(fā)揮深度學習模型的優(yōu)勢.Wu等[18]和李東東等[19]將原始一維振動信號直接輸入淺層的一維CNN實現故障診斷,但原始振動信號表現出多尺度復雜性和單一尺度的特征提取,無法充分獲取故障信息,對故障的特征提取效果不佳,診斷準確率和泛化能力有待提升.

    本研究不固定各故障的特征在信號中所占比例和位置,原始振動信號表現出多尺度復雜性,利用特征金字塔網絡(FPN)[20]結構對CNN進行改進,構建基于特征金字塔的卷積循環(huán)神經網絡(CRNN):首先利用CNN框架,結合CNN的卷積層和循環(huán)神經網絡(RNN)中的長短時記憶(LSTM)層,形成新的CRNN以充分利用CNN對空間域信息和RNN對時域信息的學習能力;然后在每一層中權值共享,減少網絡參數;最后利用FPN構建全新特征圖,對傳感器采集的信號進行特征提取,實現故障診斷.該方法可以同時處理一維原始信號和二維信號,不依賴專家經驗和人工提取特征.多尺度的特征提取豐富了診斷信息,更高層級網絡結構使得網絡的特征學習更有效.

    1 基本理論

    1.1 卷積層和池化層

    卷積核是構建卷積層的主要部分,其主要目的為提取信息特征.卷積核主要通過稀疏鏈接和權值共享在輸入信息上提取特征,其與輸入信息以局部鏈接的方式進行鏈接并滑動,滑動時利用權值共享進行計算得到信息特征值:

    (1)

    (2)

    1.2 LSTM層

    RNN是將時序信息納入計算范疇的網絡,對振動信號具有良好的解析能力.LSTM是循環(huán)神經網絡中的特種網絡,具備良好的時域特征提取特性.借助遺忘門,其可以更快收斂,減小過擬合概率.

    LSTM的主要結構如圖1所示,包含輸入門、輸出門和遺忘門,表示為

    圖1 LSTM結構Fig.1 Structure of LSTM

    ft=σ(Wfht-1+Ufxt+bf)

    it=σ(Wiht-1+Uixt+bi)

    at=tanh(Waht-1+Uaxt+ba)

    Ct=Ct-1ft+itat

    ot=σ(Woht-1+Uoxt+bo)

    ht=ottanhCt-1

    (3)

    式中:下標t表示第t輪循環(huán);ft為遺忘門;Wf、Wi、Wa和Wo為各門中的隱藏狀態(tài)權重矩陣的一維分量;Uf、Ui、Ua和Uo為輸入數據的權重矩陣的一維分量;xt為特征值;it和at為輸入門的兩條支路;Ct為細胞狀態(tài);ht為隱藏狀態(tài),由支路ot和Ct構成.

    1.3 FPN

    針對多尺度圖像的檢測共有如下3種方法[20]:

    (1) 利用圖像金字塔構建不同尺度的特征金字塔.其缺點是訓練時間較長,內存較大,并且網絡只對某一個特點的分辨率進行訓練,可能會在訓練與推理時產生矛盾.

    (2) 利用單個高層特征圖進行檢測.

    (3) 利用不同層多尺度的特征圖進行檢測.該方法沒有充分利用低層特征圖中的空間信息.

    FPN可以彌補目標檢測在處理多尺度變化問題時的不足,使得每一層不同尺度的特征圖都具有較強的語義信息.FPN可以融合低分辨率但語義信息較強和高分辨率語義信息較弱但空間信息豐富的特征圖,且僅增加較少計算量.在以特征金字塔為基礎的結構上,其能夠對每一層級的特征圖分別進行預測.具體結構如圖2所示.其中,conv為卷積運算函數,空間分辨率采用2倍采樣(2×up),然后利用逐元素相加融合上采樣的特征圖與相關從下而上的特征圖.迭代此過程,直到生成最精細的分辨率特征圖.通過附加1×1卷積層來產生最粗糙的分辨率特征圖.

    圖2 FPN結構Fig.2 Structure of FPN

    2 基于特征金字塔 CRNN 故障診斷方法

    提出一種特征金字塔的CRNN結構,CRNN結合CNN與RNN中的LSTM層,充分利用兩種網絡對時域特征和空間域特征的強大處理能力,實現對故障特征全方位多尺度特征提取,提高故障診斷的準確率.

    2.1 CRNN模塊的構建

    CRNN模塊的結構如圖3所示.其中,Relu為線性整流函數.其由卷積支路和LSTM支路并聯(lián)構成.卷積支路為卷積層,LSTM支路由LSTM層串聯(lián)池化層構成.為了配合FPN構建特征圖,要求CRNN塊的兩條支路擁有相同的語義尺寸.由于卷積的運算特性,每經過一層CNN,其實際語義強度增加(顏色漸深)同時,語義信息尺寸相應減小.而LSTM算法經過多輪運算后仍保持與原數據相同的語義尺寸,如圖4所示.因此,在LSTM支路中串聯(lián)池化層以達到改變實際語義尺寸輸出的目的.

    圖3 CRNN塊結構Fig.3 Structure of CRNN

    圖4 CNN和LSTM的語義變化Fig.4 Semantic map of CNN and LSTM

    在設定實際參數時,將LSTM支路池化層的池化窗口參數與卷積支路卷積核參數設為相同,可以保證兩支路輸出的語義信息尺寸相等.將兩支路特征進行并聯(lián)運算:

    Y2=UY2+VH+B

    (4)

    式中:Y2為并行運算的輸出;U和V分別為LSTM和卷積兩支路的合并權重;Y2為CNN支路輸出的一維展開;H為多輪LSTM支路輸出的集合;B為偏執(zhí)項.將Y2代入最大池化層后得到CRNN塊的輸出矩陣F.

    2.2 基于特征金字塔的特征圖構建

    利用特征金字塔構建特征圖可以在幾乎不增加預算量的前提下增強特征圖的數量和質量.在利用傳統(tǒng)網絡進行故障診斷時,只有最后一層網絡的輸出被納入分類器的運算輸入范疇.而特征金字塔特征提取方法能夠使用低層級網絡提取的特征,更有效地利用網絡的運算資源.

    每個CRNN塊產生一個F.設共有q個層級的CRNN塊,則其輸出分別為F1,F2,…,Fq.利用FPN構建特征圖的具體實現方法如下:

    (5)

    式中:Fz為重構后第z層的特征圖;Fr為重構前第r層的特征圖.重構特征圖均包含深層特征語義,其中部分特征圖還包含淺層特征語義.

    基于FPN的CRNN模型如圖5所示.其優(yōu)點主要有兩方面:① 在網絡的各個層級同時提取時域信息和空間域信息,使得特征信息更加豐富,有利于提升診斷準確率;② 使用FPN構建特征圖,增強各層級故障特征語義的同時節(jié)約運算資源,有利于加快故障診斷模型的收斂速度.特征金字塔的構建如圖6所示.

    圖5 特征金字塔CRNN結構Fig.5 CRNN structure based on feature pyramid

    圖6 特征金字塔構建示意圖Fig.6 Construction of feature pyramid

    2.3 方法介紹

    2.3.1構建數據集 將采集到的振動信號分為訓練集和測試集.為了取消量綱,加速收斂并增強數據可靠性,對原始振動數據進行歸一化處理:

    (6)

    圖7 順序重疊抽樣法Fig.7 Sequential overlapping sampling

    CNN具有降維特性,導致每經過一次池化運算,數據量均降低至原數據的1/4,因此要求拆分后的單個數據長度可以被22n整除.其中,n為池化層層數,且數據長度大于一個轉動周期內所采集數據的長度,即

    L≥60f/N

    其中:L為抽取數據長度;f為采樣頻率;N為轉速.

    將一維信號堆疊為二維信號,以滿足在減少對原始信號改動的同時,在兩支路分別輸入二維信號和一維信號.該方法可逆,將堆疊后的信號重新展開為一維信號即可還原信號.在卷積支路中使用堆疊后的二維信號,在LSTM支路中使用一維信號.堆疊方法如圖8所示.其中,一維信號長度為mna;二維信號長度為m,寬度為na.

    圖8 信號堆疊法Fig.8 Signal stacking method

    2.3.2網絡訓練 構建如圖9所示的網絡結構并初始化網絡參數.各支路的具體參數包括:卷積核的寬度Wk、高度Hk、數量Kk和步長Sk;LSTM隱藏層數、LSTM細胞隱藏向量維度Dl和循環(huán)次數Kl;池化層窗口寬度Wp、高度Hp和步長Sp;全連接層隱藏層細胞數量以及分類器分類數量.為保證語義信息尺寸一致性,令Wk=Wp、Hk=Hp、Sk=Sp;為保證LSTM中間支路隱藏狀態(tài)作為特征輸出,設

    圖9 基于特征金字塔CRNN的故障診斷流程Fig.9 Fault diagnosis flow based on feature pyramid CRNN

    定Dl等于輸入信號單元素的維度; 為保證兩支路輸出特征數量相等,令Kk=Kl.學習率采用指數衰減法,初始設為0.02,衰減系數設為0.90.共進行80次訓練,每次輸入10個從數據集中隨機抽取的樣本.編譯語言為Python,模型運行框架為Keras+Sklearn,硬件環(huán)境為Intel Core i5 2.5 GHz+Inter Iris Plus Graphics 655.對網絡輸入訓練集樣本進行訓練.

    2.3.3網絡測試 網絡經多次反向傳播更新后,得到最優(yōu)參數,并輸出故障類型和診斷準確率.使用t-分布隨機近鄰嵌入(t-SNE)可視化網絡內部對特征的學習效果.具體流程如圖9所示.

    3 試驗驗證

    3.1 故障模擬試驗

    行星齒輪箱故障數據由試驗臺(見圖10)測得.電機、扭矩儀、行星齒輪箱和磁粉制動器通過聯(lián)軸器相互聯(lián)結實現傳動.行星齒輪箱兩側的扭矩儀對行星齒輪箱的輸入軸和輸出軸施加徑向載荷,磁粉制動器施加軸向載荷.行星齒輪箱參數如表1所示.

    圖10 行星齒輪箱試驗臺Fig.10 Planetary gearbox test bed

    表1 行星齒輪箱基本參數Tab.1 Parameters of planetary gearbox

    3.2 故障信號采集

    在行星齒輪箱中放置三向加速度傳感器,用以收集縱向振動信號構建試驗數據集.利用替換故障件的方法模擬行星齒輪箱故障.采樣頻率設為 20.48 kHz,電機轉速為 1 980 r/min.每組采集到的振動信號點為 3 276 750 個,樣本分割窗尺寸設為 1 024 個點.采用順序重疊抽樣法,共生成 4 000 個樣本.其中,訓練集樣本 3 200 個,測試集樣本800個.數據集基本信息如表2所示.模型分類器的神經元個數與故障標簽數相同.圖11為采集信號時域圖示例.

    表2 行星齒輪箱故障數據集信息Tab.2 Data information of planetary gearbox fault

    圖11 采集時域信號示例Fig.11 Signal acquisition of time domain

    3.3 網絡參數設定

    從采集到的振動信號中抽取每種工況的樣本信號,約80%作為訓練集,約20%作為測試集.抽取信號的步長設為40.選取所抽取樣本信號的每個樣本堆疊為24×128 的二維樣本.使用該信號組進行診斷,網絡參數設定如表3所示.

    表3 網絡參數設定Tab.3 Parameters of network

    3.4 診斷準確率

    將預處理后的數據代入設定好的網絡中,每迭代200輪計算一次故障診斷準確率(ε),所得曲線如圖12所示.學習率設為0.001,經過約 6 000 輪迭代后,學習率設為 0.000 1,網絡的準確率得以收斂.

    圖12 診斷準確率曲線Fig.12 Diagnostic accuracy

    在此基礎上對網絡診斷結果進行K折交叉驗證.K折交叉驗證是評估模型性能的重要驗證方法:將數據集隨機劃分為K組,選擇其中一組作為輸入數據進行測試,其余組依次進行測試,并重復若干次劃分.為實現K折交叉驗證,試驗數據總量需要大于單次試驗所需最小數據量的K倍.設置K=5,通過調整步長獲得更多數據.驗證結果如表4所示.可知準確率的最大值為100%,平均值為99.20%,表明網絡訓練沒有過擬合或欠擬合.

    表4 網絡5折交叉驗證準確率Tab.4 50% cross validation accuracy of network

    3.5 泛用性驗證

    基于特征金字塔的CRNN對試驗臺的故障診斷準確率較高,為驗證其對其他實驗數據的泛用性,利用凱斯西儲大學的滾動軸承數據集進行泛用性驗證,軸承試驗臺如圖13所示.其中,軸承轉速為 1 750 r/min,采樣頻率為12 kHz.軸承故障數據集基本信息如表5所示.對該數據集的診斷結果做5折交叉驗證,診斷準確率結果如圖14所示.

    圖13 凱斯西儲大學軸承試驗臺Fig.13 Bearing test bed of Case Western Reserve University

    圖14 凱斯西儲大學數據診斷準確率Fig.14 Data diagnosis accuracy rate of Case Western Reserve University

    表5 軸承狀態(tài)數據信息Tab.5 Data information of bearing status

    根據驗證結果可知,在該滾動軸承的數據集上,基于特征金字塔的CRNN能很好地實現故障診斷,準確率可達99.5%以上.

    3.6 效果對比

    利用行星齒輪箱試驗臺數據對傳統(tǒng)模型進行故障診斷對比驗證,以驗證基于特征金字塔的CRNN在故障診斷方面的優(yōu)勢.分別對DNN、CNN和LSTM進行5折交叉驗證,重復3次,結果如表6所示.其中,DNN結構為1層輸入層、3層隱藏層和1層輸出層;CNN結構為3層卷積層、3層池化層和1層輸出層;LSTM包括 3 072 個細胞,重復訓練64輪.

    表6 DNN、CNN和LSTM診斷準確率Tab.6 Diagnostic accuracy of DNN, CNN, and LSTM

    計算得到DNN、CNN和LSTM故障診斷準確率的平均值分別為62.66%、86.50%和95.63%.可知,傳統(tǒng)DNN的診斷效果最差,LSTM因為針對時域信息更加有效,所以比CNN診斷準確率更高.基于特征金字塔的CRNN的診斷準確率平均值則比LSTM高3.62%,證明該方法具有優(yōu)越性.

    3.7 t-SNE可視化分析

    t-SNE方法具有可視化效果好、低維數據的分布和原始數據特征空間分布重合度高等特點.對完成訓練的特征金字塔CRNN中卷積池化層和全連接層的輸出特征進行t-SNE降維可視化分析,如圖15所示.在圖15(a)中,由卷積池化層提取的特征,經過降維后具有良好的聚類效果.在圖15(b)中,經過特征金字塔和全連接層篩選后,不同故障類別的特征得到了較好的分類效果,驗證了該方法提取目標域故障特征的有效性.

    圖15 行星齒輪箱試驗臺數據集t-SNE分析Fig.15 t-SNE analysis of data set of planetary gearbox test bench

    4 結論

    (1) 利用CNN和LSTM,特征金字塔CRNN可以同時使用一維和二維信號并行提取信號特征,并使用特征金字塔進行特征構建,實現在較小網絡層級和較低系統(tǒng)內存占用前提下的訓練和故障診斷.

    (2) 利用基于特征金字塔CRNN故障診斷方法對振動信號進行診斷,僅對振動峰值時域信號進行簡單歸一化和堆疊,簡化了人工信號優(yōu)化流程.該方法可以很好地實現對振動信號的識別和分類,經行星齒輪箱故障數據集測試,得到的診斷準確率平均值為99.20%.

    (3) 使用凱斯西儲大學數據集驗證了基于特征金字塔CRNN對振動信號診斷的泛用性,利用t-SNE對所提方法的網絡中全連接層故障特征進行了降維和可視化,驗證了該方法提取目標域故障特征的有效性.

    猜你喜歡
    支路金字塔故障診斷
    “金字塔”
    A Study of the Pit-Aided Construction of Egyptian Pyramids
    基于限流可行方案邊界集的最優(yōu)支路投切
    能源工程(2020年6期)2021-01-26 00:55:22
    海上有座“金字塔”
    神秘金字塔
    童話世界(2017年11期)2017-05-17 05:28:25
    多支路兩跳PF協(xié)作系統(tǒng)的誤碼性能
    電信科學(2016年9期)2016-06-15 20:27:30
    利用支路參數的狀態(tài)估計法辨識拓撲錯誤
    因果圖定性分析法及其在故障診斷中的應用
    多并聯(lián)支路型可控電抗器短路電抗對支路電抗和電流的影響
    基于LCD和排列熵的滾動軸承故障診斷
    色视频在线一区二区三区| 97精品久久久久久久久久精品| 51国产日韩欧美| 韩国av在线不卡| 久久人人爽人人片av| 香蕉精品网在线| av黄色大香蕉| 岛国毛片在线播放| 久久久久国产网址| 一个人看的www免费观看视频| 99久久九九国产精品国产免费| 久久女婷五月综合色啪小说 | 纵有疾风起免费观看全集完整版| 亚洲av中文av极速乱| 深夜a级毛片| 3wmmmm亚洲av在线观看| 免费播放大片免费观看视频在线观看| 亚洲国产精品专区欧美| 国产精品人妻久久久影院| 久久久精品94久久精品| 欧美性猛交╳xxx乱大交人| av国产精品久久久久影院| 日韩欧美精品免费久久| 丝袜脚勾引网站| 又爽又黄无遮挡网站| av线在线观看网站| 国产色婷婷99| 中文资源天堂在线| 99久久九九国产精品国产免费| 日本免费在线观看一区| 午夜福利高清视频| 日本与韩国留学比较| 亚洲天堂av无毛| 亚洲av.av天堂| 精品久久久久久久久av| 免费大片18禁| 亚洲欧美清纯卡通| 99热全是精品| 国产久久久一区二区三区| 亚洲三级黄色毛片| 特级一级黄色大片| 国产熟女欧美一区二区| 日韩av不卡免费在线播放| 精品久久久久久久久亚洲| av又黄又爽大尺度在线免费看| 搡老乐熟女国产| 美女内射精品一级片tv| 午夜爱爱视频在线播放| 一级av片app| 国产高清三级在线| 午夜免费男女啪啪视频观看| 成年版毛片免费区| 午夜福利高清视频| 久久精品久久久久久噜噜老黄| xxx大片免费视频| 午夜免费鲁丝| 禁无遮挡网站| 欧美3d第一页| 免费少妇av软件| 精品少妇黑人巨大在线播放| 黄色一级大片看看| 亚洲精品aⅴ在线观看| 久久国产乱子免费精品| 黑人高潮一二区| 少妇裸体淫交视频免费看高清| 久久久久九九精品影院| 人妻系列 视频| 国产综合懂色| 国产精品国产av在线观看| 成年版毛片免费区| 日韩av不卡免费在线播放| 搡女人真爽免费视频火全软件| 有码 亚洲区| 卡戴珊不雅视频在线播放| 国产淫片久久久久久久久| 如何舔出高潮| 欧美潮喷喷水| h日本视频在线播放| 91午夜精品亚洲一区二区三区| 一区二区av电影网| 亚洲精品国产成人久久av| 王馨瑶露胸无遮挡在线观看| 亚洲精品成人av观看孕妇| 免费黄网站久久成人精品| 国产精品国产av在线观看| 尾随美女入室| 三级男女做爰猛烈吃奶摸视频| 日韩av在线免费看完整版不卡| 久久99热这里只有精品18| 狂野欧美激情性xxxx在线观看| videossex国产| 日韩一区二区三区影片| 涩涩av久久男人的天堂| 国产伦在线观看视频一区| 身体一侧抽搐| 最新中文字幕久久久久| 国产视频首页在线观看| 婷婷色综合www| 亚洲精品中文字幕在线视频 | 国产精品秋霞免费鲁丝片| 丰满乱子伦码专区| 插阴视频在线观看视频| 一级av片app| 午夜免费观看性视频| 亚洲av二区三区四区| av在线观看视频网站免费| 爱豆传媒免费全集在线观看| 欧美zozozo另类| 少妇人妻精品综合一区二区| 中文在线观看免费www的网站| 欧美+日韩+精品| 亚洲av欧美aⅴ国产| 久久久成人免费电影| 国产精品精品国产色婷婷| 日本欧美国产在线视频| 午夜精品国产一区二区电影 | 3wmmmm亚洲av在线观看| 日韩大片免费观看网站| 国产高清有码在线观看视频| 在线免费十八禁| 亚洲不卡免费看| 欧美日本视频| 大话2 男鬼变身卡| 日韩强制内射视频| 女的被弄到高潮叫床怎么办| 欧美区成人在线视频| 亚洲国产欧美在线一区| 久久精品国产亚洲网站| 一个人看视频在线观看www免费| 日韩中字成人| 亚洲欧美清纯卡通| 午夜免费男女啪啪视频观看| 美女主播在线视频| 毛片女人毛片| 日韩在线高清观看一区二区三区| 精品久久久久久久久亚洲| 国产av国产精品国产| 高清在线视频一区二区三区| 欧美精品人与动牲交sv欧美| 少妇裸体淫交视频免费看高清| 18禁动态无遮挡网站| 黄色日韩在线| 国产男女内射视频| 国产精品伦人一区二区| 精品国产露脸久久av麻豆| 国产男人的电影天堂91| 午夜激情福利司机影院| 黄色配什么色好看| xxx大片免费视频| 一级黄片播放器| 国产一区二区三区综合在线观看 | 69av精品久久久久久| 老女人水多毛片| av网站免费在线观看视频| 欧美一区二区亚洲| 国产亚洲av片在线观看秒播厂| 国内精品美女久久久久久| 噜噜噜噜噜久久久久久91| 欧美 日韩 精品 国产| 99视频精品全部免费 在线| 一级片'在线观看视频| 三级国产精品片| 五月开心婷婷网| 97超视频在线观看视频| 午夜激情福利司机影院| 国产黄片美女视频| 麻豆国产97在线/欧美| 天堂俺去俺来也www色官网| 精品酒店卫生间| 亚州av有码| 麻豆乱淫一区二区| 国产老妇女一区| 激情五月婷婷亚洲| 高清毛片免费看| 日本黄大片高清| 欧美性猛交╳xxx乱大交人| 国产爽快片一区二区三区| 日本黄大片高清| 欧美一区二区亚洲| tube8黄色片| 亚洲精品亚洲一区二区| 午夜福利视频精品| 国产精品国产av在线观看| 久久精品国产自在天天线| 大香蕉97超碰在线| 欧美人与善性xxx| 我的女老师完整版在线观看| 免费黄色在线免费观看| 人人妻人人澡人人爽人人夜夜| 国产国拍精品亚洲av在线观看| 欧美3d第一页| 大话2 男鬼变身卡| 国产av国产精品国产| 成人美女网站在线观看视频| 成人黄色视频免费在线看| 一级毛片我不卡| 国产成人午夜福利电影在线观看| 美女国产视频在线观看| 韩国高清视频一区二区三区| 肉色欧美久久久久久久蜜桃 | 欧美最新免费一区二区三区| 亚洲精品国产成人久久av| 日本三级黄在线观看| 国内精品美女久久久久久| 天美传媒精品一区二区| 国产色爽女视频免费观看| 免费高清在线观看视频在线观看| 插阴视频在线观看视频| 亚洲在久久综合| 777米奇影视久久| 晚上一个人看的免费电影| 国产精品久久久久久久电影| 国产精品精品国产色婷婷| 中文在线观看免费www的网站| 日产精品乱码卡一卡2卡三| 美女内射精品一级片tv| 亚洲人成网站在线观看播放| 亚洲国产精品专区欧美| 女人被狂操c到高潮| 国产男女超爽视频在线观看| 一级a做视频免费观看| 成人亚洲精品av一区二区| 成人毛片60女人毛片免费| 91在线精品国自产拍蜜月| 精品99又大又爽又粗少妇毛片| 国产精品国产av在线观看| 狂野欧美白嫩少妇大欣赏| 91狼人影院| 国产伦精品一区二区三区视频9| 欧美xxxx黑人xx丫x性爽| av黄色大香蕉| 在线观看av片永久免费下载| 国产精品女同一区二区软件| 久久久久久久大尺度免费视频| 免费大片黄手机在线观看| eeuss影院久久| 午夜老司机福利剧场| 日韩一区二区视频免费看| 国产精品成人在线| 久久97久久精品| 老女人水多毛片| 亚洲精品日韩在线中文字幕| 久久6这里有精品| 亚洲天堂国产精品一区在线| 国产永久视频网站| 一级毛片久久久久久久久女| 老师上课跳d突然被开到最大视频| 国内少妇人妻偷人精品xxx网站| 青青草视频在线视频观看| 九九爱精品视频在线观看| 欧美区成人在线视频| 午夜激情久久久久久久| 精品人妻熟女av久视频| 亚洲人成网站在线观看播放| 国产精品无大码| 一本久久精品| 久久久精品免费免费高清| 国产乱人视频| 久久久久久久精品精品| 色视频www国产| 国产精品99久久久久久久久| 深夜a级毛片| 精品视频人人做人人爽| 免费av毛片视频| 女人被狂操c到高潮| 午夜福利网站1000一区二区三区| 亚洲精品影视一区二区三区av| 久久亚洲国产成人精品v| 在线观看三级黄色| 欧美极品一区二区三区四区| 黄片无遮挡物在线观看| 国产精品99久久99久久久不卡 | 日本爱情动作片www.在线观看| 可以在线观看毛片的网站| 菩萨蛮人人尽说江南好唐韦庄| 国产69精品久久久久777片| 日韩欧美 国产精品| 久久韩国三级中文字幕| 亚洲在久久综合| 欧美另类一区| 免费人成在线观看视频色| 国产精品av视频在线免费观看| 99热全是精品| 国产永久视频网站| 老司机影院成人| 亚洲av二区三区四区| 一级毛片久久久久久久久女| 国产又色又爽无遮挡免| 日韩av免费高清视频| videos熟女内射| 黄片wwwwww| 国产视频内射| 亚洲精品乱码久久久久久按摩| 免费观看的影片在线观看| 亚洲在久久综合| 91久久精品国产一区二区三区| 亚洲在线观看片| 九九久久精品国产亚洲av麻豆| 成人鲁丝片一二三区免费| 日韩人妻高清精品专区| 免费观看的影片在线观看| 国产黄片美女视频| 国产成人免费无遮挡视频| 亚洲欧美成人精品一区二区| 亚洲国产欧美人成| 在线观看人妻少妇| 精品久久久久久电影网| 亚洲精品乱码久久久v下载方式| 精品人妻视频免费看| 国产精品99久久99久久久不卡 | 一级黄片播放器| 交换朋友夫妻互换小说| 国产亚洲精品久久久com| 成年版毛片免费区| 一二三四中文在线观看免费高清| 天天一区二区日本电影三级| 午夜免费鲁丝| 亚洲在久久综合| 嫩草影院精品99| 99久久精品一区二区三区| 精品午夜福利在线看| 亚洲精品一区蜜桃| 草草在线视频免费看| 男人狂女人下面高潮的视频| 亚洲欧美日韩无卡精品| 免费看av在线观看网站| 亚洲av中文av极速乱| 简卡轻食公司| 国产永久视频网站| 亚洲欧美一区二区三区国产| 在线观看av片永久免费下载| 黄色视频在线播放观看不卡| av网站免费在线观看视频| 日本欧美国产在线视频| 日本黄大片高清| 午夜老司机福利剧场| 三级经典国产精品| 嫩草影院入口| 三级男女做爰猛烈吃奶摸视频| 女人十人毛片免费观看3o分钟| 人妻一区二区av| 狂野欧美激情性xxxx在线观看| 可以在线观看毛片的网站| 欧美精品国产亚洲| 中文字幕久久专区| 高清av免费在线| 国产精品麻豆人妻色哟哟久久| www.色视频.com| 如何舔出高潮| 久久精品国产自在天天线| 亚洲人与动物交配视频| 久久精品人妻少妇| 亚洲精品成人av观看孕妇| 日韩一本色道免费dvd| 久久久色成人| 成人毛片60女人毛片免费| 国产精品嫩草影院av在线观看| 免费观看av网站的网址| 99久久精品国产国产毛片| 最近中文字幕2019免费版| 天天躁日日操中文字幕| 精品酒店卫生间| 69av精品久久久久久| 国产欧美亚洲国产| 在线观看免费高清a一片| 色综合色国产| 亚洲精品影视一区二区三区av| 女人久久www免费人成看片| av国产精品久久久久影院| 久久精品国产亚洲av涩爱| 99九九线精品视频在线观看视频| 亚洲在久久综合| 校园人妻丝袜中文字幕| 色综合色国产| 精品久久久久久久末码| 国产熟女欧美一区二区| 国产亚洲精品久久久com| 日本-黄色视频高清免费观看| 日韩免费高清中文字幕av| 熟女电影av网| 国产精品熟女久久久久浪| 寂寞人妻少妇视频99o| 亚洲av免费高清在线观看| 欧美三级亚洲精品| 亚洲综合精品二区| 日本黄大片高清| 久久国内精品自在自线图片| 欧美精品人与动牲交sv欧美| 三级国产精品片| av播播在线观看一区| 超碰97精品在线观看| 国产真实伦视频高清在线观看| 韩国av在线不卡| 国产黄频视频在线观看| 晚上一个人看的免费电影| 少妇被粗大猛烈的视频| 97在线人人人人妻| 亚洲国产欧美在线一区| 午夜视频国产福利| 青青草视频在线视频观看| 亚洲性久久影院| 久久99热这里只频精品6学生| 亚洲最大成人中文| videossex国产| 欧美日韩一区二区视频在线观看视频在线 | 亚洲最大成人手机在线| av在线蜜桃| 人妻少妇偷人精品九色| 蜜臀久久99精品久久宅男| 久久精品久久久久久久性| 91aial.com中文字幕在线观看| 国产精品伦人一区二区| 日本wwww免费看| 老师上课跳d突然被开到最大视频| 亚洲av免费在线观看| 国产探花在线观看一区二区| 久久ye,这里只有精品| 91狼人影院| av在线播放精品| 亚洲色图av天堂| 成人无遮挡网站| 亚洲国产高清在线一区二区三| 久久久久网色| av在线观看视频网站免费| 免费看日本二区| 亚洲欧洲国产日韩| 最近最新中文字幕大全电影3| 日韩三级伦理在线观看| 小蜜桃在线观看免费完整版高清| 成人漫画全彩无遮挡| 日韩,欧美,国产一区二区三区| 久久97久久精品| 黄色日韩在线| 久久久久网色| 日本与韩国留学比较| 香蕉精品网在线| 久久99热这里只有精品18| tube8黄色片| 有码 亚洲区| 中文资源天堂在线| 少妇猛男粗大的猛烈进出视频 | 国产 精品1| 精品人妻偷拍中文字幕| 婷婷色av中文字幕| 美女国产视频在线观看| 日韩电影二区| 在线观看美女被高潮喷水网站| 不卡视频在线观看欧美| 成人无遮挡网站| 国产免费一级a男人的天堂| 99热6这里只有精品| 黄片无遮挡物在线观看| 国产在线一区二区三区精| 五月玫瑰六月丁香| 免费看不卡的av| 91在线精品国自产拍蜜月| 久久久久久久久久人人人人人人| 亚洲国产最新在线播放| 中文天堂在线官网| 亚洲精品456在线播放app| 九九久久精品国产亚洲av麻豆| 亚洲av成人精品一区久久| 欧美国产精品一级二级三级 | 80岁老熟妇乱子伦牲交| 久久精品久久精品一区二区三区| 熟女电影av网| 在线免费观看不下载黄p国产| 亚洲精品,欧美精品| 亚洲精品第二区| 2022亚洲国产成人精品| 久久人人爽人人片av| 亚洲真实伦在线观看| 国产精品国产三级国产专区5o| 神马国产精品三级电影在线观看| 国产成人aa在线观看| 天天躁夜夜躁狠狠久久av| 国内精品宾馆在线| 国产免费福利视频在线观看| 丝袜脚勾引网站| 精品人妻视频免费看| 又爽又黄a免费视频| 日本猛色少妇xxxxx猛交久久| 欧美高清成人免费视频www| av国产精品久久久久影院| 在线天堂最新版资源| 少妇 在线观看| 精品人妻视频免费看| 乱码一卡2卡4卡精品| 久久久久久久午夜电影| 国产毛片a区久久久久| 日本一二三区视频观看| 亚洲真实伦在线观看| 亚洲图色成人| 一级a做视频免费观看| 日韩国内少妇激情av| 丰满少妇做爰视频| 91午夜精品亚洲一区二区三区| 欧美成人精品欧美一级黄| 中国三级夫妇交换| 中国美白少妇内射xxxbb| 亚洲精品影视一区二区三区av| 纵有疾风起免费观看全集完整版| 国产v大片淫在线免费观看| 我的女老师完整版在线观看| 亚洲精华国产精华液的使用体验| 美女内射精品一级片tv| 久久久亚洲精品成人影院| 中文字幕亚洲精品专区| 一级片'在线观看视频| 99re6热这里在线精品视频| 亚洲综合精品二区| 只有这里有精品99| 精品人妻一区二区三区麻豆| 国产亚洲最大av| 亚洲av二区三区四区| 黄色视频在线播放观看不卡| 久久久久精品久久久久真实原创| 欧美潮喷喷水| 日日摸夜夜添夜夜爱| 国产成人精品久久久久久| 日本黄色片子视频| 18禁裸乳无遮挡动漫免费视频 | 国产亚洲最大av| 三级经典国产精品| 欧美xxxx性猛交bbbb| 久久久久精品久久久久真实原创| 日韩一区二区视频免费看| 免费观看无遮挡的男女| av网站免费在线观看视频| 欧美极品一区二区三区四区| 精品久久国产蜜桃| 久久亚洲国产成人精品v| 午夜亚洲福利在线播放| 国产探花极品一区二区| 免费不卡的大黄色大毛片视频在线观看| 69人妻影院| 十八禁网站网址无遮挡 | 波野结衣二区三区在线| 一二三四中文在线观看免费高清| .国产精品久久| 日韩亚洲欧美综合| 国产永久视频网站| 又爽又黄无遮挡网站| 欧美bdsm另类| 精品亚洲乱码少妇综合久久| 欧美日韩精品成人综合77777| 亚洲国产成人一精品久久久| 日韩不卡一区二区三区视频在线| 欧美成人午夜免费资源| 黑人高潮一二区| 波多野结衣巨乳人妻| 嫩草影院新地址| 精品国产露脸久久av麻豆| 男人和女人高潮做爰伦理| 亚洲精品自拍成人| 香蕉精品网在线| 久久99热这里只有精品18| 国产乱人偷精品视频| 中文精品一卡2卡3卡4更新| 精品久久久久久久人妻蜜臀av| 亚洲精品久久午夜乱码| 在线a可以看的网站| 在线 av 中文字幕| 亚洲真实伦在线观看| 国产精品蜜桃在线观看| 久久久久久伊人网av| 秋霞伦理黄片| 免费大片18禁| 色综合色国产| 日韩欧美精品免费久久| 国产免费福利视频在线观看| 18禁动态无遮挡网站| 久久久久久久午夜电影| 最近中文字幕2019免费版| 成人鲁丝片一二三区免费| 国产男人的电影天堂91| 亚洲色图av天堂| 麻豆乱淫一区二区| 国产欧美另类精品又又久久亚洲欧美| 国产成人freesex在线| 国产人妻一区二区三区在| 少妇人妻 视频| 国产精品人妻久久久久久| 18禁在线无遮挡免费观看视频| 国产精品偷伦视频观看了| 日韩视频在线欧美| 毛片一级片免费看久久久久| 国产黄色免费在线视频| 国产一区有黄有色的免费视频| 国精品久久久久久国模美| 亚洲激情五月婷婷啪啪| 久久久a久久爽久久v久久| 亚洲色图综合在线观看| 午夜日本视频在线| 亚洲精品乱码久久久v下载方式| 青春草亚洲视频在线观看| 国内精品美女久久久久久| 成年av动漫网址| 在线观看国产h片| 三级国产精品欧美在线观看| 亚洲久久久久久中文字幕| 听说在线观看完整版免费高清| 大话2 男鬼变身卡| 欧美高清性xxxxhd video| 一级二级三级毛片免费看| 纵有疾风起免费观看全集完整版| 久久精品国产鲁丝片午夜精品| 青青草视频在线视频观看| 亚洲av欧美aⅴ国产| 亚洲天堂国产精品一区在线| 69人妻影院| 亚洲精品日韩在线中文字幕| 国产一区有黄有色的免费视频| 国产免费视频播放在线视频| 亚洲精品,欧美精品|