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

    基于CEEMD-SA-RNN的柴油機曲軸軸承磨損預測

    2022-08-25 09:36:40李英順田宇左洋張國瑩周通
    車用發(fā)動機 2022年4期
    關(guān)鍵詞:模擬退火曲軸分量

    李英順,田宇,左洋,張國瑩,周通

    (北京石油化工學院信息工程學院,北京 102617)

    步兵戰(zhàn)車是供步兵機動作戰(zhàn)用的裝甲戰(zhàn)斗車輛,具有高機動性。步兵戰(zhàn)車經(jīng)常在惡劣的環(huán)境下工作,為保障車輛的正常運行,定期對機械進行故障檢測是必要的防患措施。傳統(tǒng)的零部件故障檢測方法需要對一些部件進行拆解,在消耗人力物力的同時,對檢測者的要求也很高。而利用計算機通信技術(shù)進行機械故障檢測具有無需拆解、高效等優(yōu)點,相較于傳統(tǒng)的檢測方法更好,成為了國內(nèi)外學者研究的焦點。

    Y. Chen等[1]利用快速傅里葉變換(Fast Fourier transform,F(xiàn)FT)對振動信號進行分析和處理,但其特征具有易聚集特性,不能很好地識別,需要人為介入;張玲玲等[2]采用小波包-AR譜對變速器軸承故障進行特征提取,但該方法有小波基的選取困難缺陷;J. Cheng等[3]利用經(jīng)驗模態(tài)分解(Empirical Mode Decomposition, EMD)對振動信號進行分解,但EMD會出現(xiàn)模態(tài)混淆現(xiàn)象;W. Xia等[4]利用自適應小波和集合經(jīng)驗模態(tài)分解(Ensemble Empirical Mode Decomposition,EEMD)處理信號,用故障樹分析噴油提前角度的大小,但由于實際運算有限,不能很好地消除已添加的白噪聲的影響,會引起重構(gòu)誤差,而補充的集合經(jīng)驗模態(tài)分解方法(Complementary Ensemble Empirical Mode Decomposition,CEEMD),可以很好地消除白噪聲對信號的影響;周小麟等[5]使用補充的集合經(jīng)驗模態(tài)分解、緞藍園丁鳥優(yōu)化算法和最小二乘支持向量回歸(Least Squares Support Vector Regression,LSSVT)實現(xiàn)超短期風電功率組合預測,但此算法的擬合精度和采樣點數(shù)成正比,需要較多數(shù)據(jù),而循環(huán)神經(jīng)網(wǎng)絡(Recurrent Neural Network,RNN)可有效避免此類問題。

    柴油機曲軸軸承的振動信號采集相對困難,很容易混入噪聲。本研究對原始數(shù)據(jù)進行軟閾值處理,使用多種分解方法對已處理的數(shù)據(jù)進行分解,分析各項指標,通過奇異值分解方法對數(shù)據(jù)進行特征提取,并作為神經(jīng)網(wǎng)絡模塊的特征集。循環(huán)神經(jīng)網(wǎng)絡適用于多分類問題,適合對曲軸軸承的磨損程度進行預測。此外,利用模擬退火算法找尋循環(huán)神經(jīng)網(wǎng)絡的最優(yōu)參數(shù),能有效提高系統(tǒng)的預測準確率。

    1 算法設計

    1.1 信號處理算法

    1.1.1 補充的集合經(jīng)驗模態(tài)分解

    補充的集合經(jīng)驗模態(tài)分解是基于集合經(jīng)驗模態(tài)分解(EEMD)的想法,向信號中添加相反白噪聲,實現(xiàn)抵消其影響的新分解方法[6-8]。具體步驟如下。

    (1)

    (2)

    (3)

    (4)

    式中:imfj為最終得到的第j個本征模態(tài)函數(shù);r為分解后的余項分量。利用已得到的本征模態(tài)函數(shù)進行自回歸模型分析,可以觀測系統(tǒng)的特征和工作狀態(tài)[9],在故障診斷和預測中應用廣泛。其定義為

    (5)

    式中:c為常數(shù)項;εt為平均值是0、標準差是σ的隨機誤差,σ值恒不變。

    1.1.2 奇異值分解

    奇異值分解(Singular Value Decom-position,SVD)是根據(jù)矩陣的固有特性對矩陣進行分解的常用方法。其計算步驟如下:

    設矩陣E(e×g)的秩為r′,則存在矩陣U和W及對角陣D,使得[10]

    E=UDWT。

    (6)

    式(6)成立,則式(6)為E(e×g)的奇異值分解,其中:

    U(e×e)=[u1,u2,…ue],

    (7)

    (8)

    Δr′×r′=diag(σ1,σ2,…σr′),

    (9)

    Wg×g=[w1,w2,…wg],

    (10)

    r′=min(e,g)。

    (11)

    1.2 故障預測算法

    1.2.1 循環(huán)神經(jīng)網(wǎng)絡

    曲軸軸承磨損程度需要細致劃分,循環(huán)神經(jīng)網(wǎng)絡(RNN)不僅可以處理任意長度的輸入,還適用于多輸出問題。因此使用RNN作為預測的基本網(wǎng)絡,其結(jié)構(gòu)見圖1。

    圖1 RNN結(jié)構(gòu)

    圖1中x為輸入,h為隱層單元,o為輸出,L為損失函數(shù),y為訓練集的標簽[11-12]。

    ht=f(U·xt+W·ht-1+b),

    (12)

    ot=V·ht+c,

    (13)

    yt=g(ot)。

    (14)

    式中:f()和g()為激活函數(shù)。在RNN的訓練調(diào)參過程中,需要調(diào)優(yōu)參數(shù)只有V,W,U3個參數(shù)。V,W,U的偏導求解和所有歷史時刻的數(shù)據(jù)有關(guān)。各偏導為

    (15)

    (16)

    (17)

    (18)

    (19)

    1.2.2 模擬退火算法

    20世紀80年代,Kirkpatrick等提出了模擬退火算法[13-14](Simulated Annealing,SA),該算法來源于晶體冷卻過程,是一種概率算法。SA原理[15]如下:從某一高溫開始,利用參數(shù)控制溫度下降,結(jié)合概率的突跳性,尋找全局最優(yōu)解。該算法具有時變特性,可以有效避免陷入局部極小問題。

    算法步驟如下。

    1) 初始化:起始溫度T0(足夠高),初始解S(迭代開始),迭代次數(shù)L,終止溫度T1;

    2) 設k=1,2,…L,做第3至第6步計算;

    3) 產(chǎn)生新解S′;

    4) 計算增量ΔC=f(S′)-f(S),f()為目標函數(shù);

    6) 若滿足迭代次數(shù)則繼續(xù),否則轉(zhuǎn)第3步;

    7) 看是否滿足目標溫度,若不滿足則T逐漸減小,且T>T1,然后轉(zhuǎn)第3步。

    2 模型建立

    系統(tǒng)模型由兩部分組成:第一部分是補充的集合經(jīng)驗模態(tài)分解與奇異值分解方法相結(jié)合(CEEMD-SVD)的信號處理模塊;第二部分是模擬退火算法優(yōu)化循環(huán)神經(jīng)網(wǎng)絡(SA-RNN)的故障預測模塊。建立的系統(tǒng)模型主體框架見圖2。

    圖2 模型主體框架

    2.1 信號處理模型

    信號處理流程見圖3。磨損導致軸承游隙增大,運行精度降低,振動噪聲增大。采用連續(xù)加速的方法對曲軸軸承進行診斷,故障時的振動特征相比正常情況更加明顯,因此利用此方法提取振動信號作為依據(jù)對部件進行故障分類。

    圖3 信號處理流程

    對已采集的信號進行軟閾值處理。分別使用無偏似然估計(rigrsure)、最優(yōu)預測變量閾值(heursure)、廣義閾值(sqtwolog)和最小最大閾值(minimaxi)4種方法對信號進行處理。以結(jié)構(gòu)相似性(Structural Similarity,SSIM)作為衡量去噪效果的指標,通過結(jié)果對比選擇rigrsur閾值法對原始信號進行處理。對已處理的信號,使用經(jīng)驗模態(tài)分解、集合經(jīng)驗模態(tài)分解和補充的集合經(jīng)驗模態(tài)分解進行分解,通過對比得出4項參數(shù)指標,選取補充的集合經(jīng)驗模態(tài)分解方法作為最優(yōu)信號分解方法。對信號進行CEEMD分解后,使用自回歸模型對各工況下前6個IMF分量譜圖及累加能量譜圖進行對比分析。觀察到從第7個IMF分量開始,其能量占總能量之比很低,因此選取前6分量作為初始向量矩陣,并對其進行奇異值分解,得到的特征分量作為神經(jīng)網(wǎng)絡的特征集。

    2.2 網(wǎng)絡預測模型

    建立循環(huán)神經(jīng)網(wǎng)絡模對信號進行預測,使用模擬退火優(yōu)化循環(huán)神經(jīng)網(wǎng)絡可以更好地預測曲軸軸承的磨損程度。

    模型步驟如下:

    1) 初始化參數(shù),隨機產(chǎn)生隱含層節(jié)點數(shù)和序列個數(shù)(小于隱含層節(jié)點數(shù)),進行循環(huán)神經(jīng)網(wǎng)絡預測,得出準確率;

    2) 更新隱含層節(jié)點數(shù)和序列個數(shù)為隨機值,且序列個數(shù)小于隱含層節(jié)點數(shù),利用兩參數(shù)創(chuàng)建新的循環(huán)神經(jīng)網(wǎng)絡,并進行預測,得出新網(wǎng)絡的準確率;

    3) 比較新舊準確率的大小,若新值大則接收并存儲新值,否則以概率接受新狀態(tài)(Metropolis準則);

    4) 看是否滿足迭代次數(shù),若不滿足則轉(zhuǎn)第2步;

    5) 看是否滿足目標溫度,若不滿足則T逐漸減小,然后轉(zhuǎn)第2步;滿足則得出最優(yōu)解并結(jié)束。

    模擬退火算法優(yōu)化的循環(huán)神經(jīng)網(wǎng)絡步驟見圖4。

    圖4 SA-RNN步驟

    3 試驗驗證

    以某步兵戰(zhàn)車柴油機第4道曲軸軸承為研究對象。選用PCB601A01型ICP工業(yè)加速度傳感器,將其放置于油底與缸體結(jié)合處3缸左右位置,通過加裝強力磁座吸合于測試部位表面。定義軸承游隙0.05~0.11 mm為正常工況,0.11~0.15 mm為輕度磨損,0.15~0.20 mm為中度磨損,0.20 mm以上為重度磨損。采集轉(zhuǎn)速1 800 r/min正常工況、輕微磨損、中度磨損和嚴重磨損下的曲軸軸承振動信號,傳感器的采樣頻率為12.8 kHz,在運行平穩(wěn)后,每次采樣16 384個點作為一條數(shù)據(jù)。使用該數(shù)據(jù),分別利用CEEMD-SVD算法和SA-RNN算法進行故障預測分析。

    3.1 CEEMD-SVD試驗

    考慮現(xiàn)場采集原因,本研究首先對原始信號進行小波閾值去噪處理,因為硬閾值會引起吉布斯效應,所以選擇軟閾值處理。使用無偏似然估計、最優(yōu)預測變量閾值、廣義閾值、最小最大閾值,并使用SSIM作為衡量去噪效果的指標。正常信號使用不同閾值方法去噪結(jié)果見圖5。

    圖5 正常信號使用不同閾值方法去噪結(jié)果

    由表1中SSIM值可以看出,使用無偏似然估計法去噪效果最好。

    表1 不同閾值方法得到的SSIM值

    對已處理的信號進行傅里葉變換,得到的結(jié)果見圖6。

    圖6 各工況下的FFT結(jié)果

    由圖6可知,經(jīng)過傅里葉變換后的正常工況信號有3個明顯峰值,分別在185 Hz,500 Hz和1 270 Hz左右;輕微磨損下的信號有3個明顯峰值,分別在380 Hz,1 220 Hz和2 050 Hz左右;中度磨損下的信號有2個明顯峰值,分別出現(xiàn)在480 Hz和1 235 Hz附近;嚴重磨損信號有4個明顯峰值,分別出現(xiàn)在185 Hz,490 Hz,1 200 Hz和1 700 Hz附近。

    對已去噪的信號,分別使用EMD、EEMD(0.210 0)、CEEMD(0.110 0)方法進行分解并比較。比較參數(shù)有均方誤差、正交性、分量個數(shù)和時間(見表2)。

    表2 正常去噪信號分解方法比較

    由表2可得出,最優(yōu)分解方法為CEEMD分解,EMD和EEMD算法存在一定的模態(tài)混淆現(xiàn)象,從結(jié)果可以看出,EEMD在多次嘗試改變初始參數(shù)后均方誤差及正交性的參數(shù)值依然較大,運行時間也比EMD長。造成該問題的原因可能是在有限次的運行過程中無法消除白噪聲的影響,致使重構(gòu)誤差較大,完備性差。CEEMD算法經(jīng)過多次改變參數(shù)后得到的均方誤差及正交性都優(yōu)于EMD算法,后續(xù)采用CEEMD算法作為系統(tǒng)模型的信號處理方法。

    進行CEEMD分解,得到10個IMF分量及余項r。信號CEEMD的能量譜見圖7。

    圖7 能量譜(CEEMD)

    由能量譜可得出,從第7個分量開始能量很低。采用前6個分量作為主要研究對象,對前6個IMFs進行傅里葉變換,其頻譜圖見圖8。

    圖8 CEEMD-IMF1~6的傅里葉變換頻譜

    作為類比,本研究也使用了EMD及EEMD方法得到相應的傅里葉變換頻譜圖(見圖9和圖10)。

    圖9 EMD-IMFs1~6的傅里葉變換頻譜

    由圖8至圖10對比發(fā)現(xiàn),前6個分量的頻率范圍有明顯的區(qū)分度, CEEMD方法下的混疊現(xiàn)象更小,再次證明針對曲軸軸承的振動信號CEEMD方法更加適用。

    圖10 EEMD-IMFs1~7的傅里葉變換頻譜

    對正常工況進行CEEMD分解并作出Hilbert邊際譜(見圖11)。

    圖11 CEEMD-Hilbert邊際譜

    由圖11可見,Hilbert邊際譜中3個峰值和之前FFT變換下的峰值相近,在180 Hz,500 Hz及1 300 Hz左右。

    各工況的前6個分量的CEEMD-AR譜對比見圖12,各工況下前6個IMF分量CEEMD-AR譜累加能量對比見圖13。

    圖12 正常與3種磨損情況下的前6個分量CEEMD-AR譜對比

    圖13 各工況下前6個IMF分量CEEMD-AR譜累加能量對比

    對圖12、圖13進行分析可得以下結(jié)論。

    1) CEEMD-AR譜的功率表幅值能量主要集中在6 500 Hz以內(nèi)。

    2) 從分量信號圖中可以看出,能量由正常、輕微磨損、中度磨損和重度磨損依次遞增的特征頻帶有IMF1分量3 850 Hz以內(nèi),IMF2分量1 450 Hz以內(nèi),IMF3分量875 Hz以內(nèi),IMF4分量450 Hz以內(nèi),IMF5分量250 Hz以內(nèi),IMF6分量100 Hz以內(nèi)。曲軸軸承輕微磨損和中度磨損遞增的特征頻帶有IMF1分量4 850~5 550 Hz。

    3) 從累加能量對比圖可以看出,正常磨損能量遞增為75 Hz以內(nèi),輕微磨損能量遞增為50 Hz以內(nèi),中度磨損能量遞增為75 Hz以內(nèi),重度磨損能量遞增范圍為50 Hz以內(nèi)。

    使用SVD對信號進行特征提取,考慮能量占比,使用前6個分量作為柴油機曲軸軸承磨損程度的特征,并作為深度學習的數(shù)據(jù)集。利用CEEMD-SVD方法提取的柴油機曲軸軸承磨損故障特征的數(shù)據(jù)集(部分)見表3。

    表3 訓練數(shù)據(jù)與測試數(shù)據(jù)(部分)

    3.2 SA-RNN試驗

    使用RNN對數(shù)據(jù)進行訓練及預測,參數(shù)設置如下:迭代次數(shù)2 000,批次大小1,學習率0.01 ,輸入6,輸出1,序列長度4,隱含層10。得到120個測試集的預測結(jié)果(見圖14)。

    圖14 測試集預測值與實際值對比

    預測值的準確率為92.307 7%,雖然得到的準確率大于90%,但并不滿足大于95%的條件。本研究嘗試使用模擬退火算法對循環(huán)神經(jīng)網(wǎng)絡進行優(yōu)化,設置隱含層節(jié)點數(shù)和序列長度為可變參數(shù)x1和x2,準確率y作為模擬退火算法的目標函數(shù)。約束條件為隱含層節(jié)點數(shù)小于等于10,并求準確率的最大值、最優(yōu)隱含層節(jié)點數(shù)和序列長度。設置初始溫度為100,截止溫度9,鏈長10,降溫速率0.9。通過迭代,得到每次溫度下降鏈長中的最優(yōu)解(見表4)。

    表4 SA-RNN調(diào)參及準確率

    SA-RNN退火效果見圖15。由圖15可見,使用模擬退火算法優(yōu)化的循環(huán)神經(jīng)網(wǎng)絡在每次降溫時準確率上升速度很快,準確率在達到95%后逐漸趨于穩(wěn)定,這充分說明該算法的效率較高、穩(wěn)定性良好。最后得到最優(yōu)隱含層節(jié)點數(shù)為8,序列長度為2,準確率為97.479%,相比未進行優(yōu)化的準確率為92.307 7%的RNN,優(yōu)化效果明顯,其誤差也在可接受范圍內(nèi)。綜上可見,使用模擬退火算法進行優(yōu)化可以降低局限性,更好獲得全局最優(yōu)解,同時具有良好的魯棒性、穩(wěn)定性,這說明該系統(tǒng)具有可行性。

    圖15 SA-RNN退火效果

    4 結(jié)束語

    利用柴油機曲軸軸承振動信號來預測判斷軸承的磨損程度,提出了一種補充的集合經(jīng)驗模態(tài)分解與奇異值分解相結(jié)合提取信號特征,使用模擬退火算法優(yōu)化循環(huán)神經(jīng)網(wǎng)絡對故障進行預測的方法,建立了CEEMD-SA-RNN模型。利用柴油機曲軸軸承振動信號對提出的算法進行了試驗驗證,結(jié)果表明CEEMD在對曲軸軸承振動信號分解方面更加適用,使用模擬退火算法的循環(huán)神經(jīng)網(wǎng)絡相比常規(guī)循環(huán)神經(jīng)網(wǎng)絡,得到的準確率提高了5%以上。

    猜你喜歡
    模擬退火曲軸分量
    帽子的分量
    淺析曲軸磨床的可靠性改進
    一物千斤
    智族GQ(2019年9期)2019-10-28 08:16:21
    一種曲軸豎置熱處理裝置設計
    模擬退火遺傳算法在機械臂路徑規(guī)劃中的應用
    論《哈姆雷特》中良心的分量
    基于支持向量機的曲軸動平衡控制方法
    分量
    基于模糊自適應模擬退火遺傳算法的配電網(wǎng)故障定位
    SOA結(jié)合模擬退火算法優(yōu)化電容器配置研究
    国内精品久久久久久久电影| 成人永久免费在线观看视频| 夜夜看夜夜爽夜夜摸| ponron亚洲| 久久精品aⅴ一区二区三区四区| av有码第一页| 久久性视频一级片| 国产又黄又爽又无遮挡在线| 国产黄片美女视频| 亚洲国产看品久久| 夜夜爽天天搞| 精品久久久久久成人av| 黑人操中国人逼视频| 妹子高潮喷水视频| 国产欧美日韩一区二区三| 久久久久久久精品吃奶| 欧洲精品卡2卡3卡4卡5卡区| 女性被躁到高潮视频| 色在线成人网| 成人午夜高清在线视频 | 熟女电影av网| 色综合亚洲欧美另类图片| 69av精品久久久久久| 国产成人欧美在线观看| 亚洲国产精品合色在线| 国产精品香港三级国产av潘金莲| 久久人妻福利社区极品人妻图片| 亚洲欧美精品综合久久99| 亚洲一区二区三区不卡视频| 搡老岳熟女国产| 精品福利观看| 亚洲专区国产一区二区| 国产高清激情床上av| 亚洲一码二码三码区别大吗| 国产真实乱freesex| 国产久久久一区二区三区| 亚洲专区字幕在线| 18禁黄网站禁片免费观看直播| 露出奶头的视频| 黄片大片在线免费观看| 级片在线观看| 一进一出抽搐gif免费好疼| 中文字幕人妻丝袜一区二区| 亚洲人成电影免费在线| 搡老熟女国产l中国老女人| 国产精品久久久久久亚洲av鲁大| 免费观看人在逋| 热99re8久久精品国产| 精品一区二区三区av网在线观看| 韩国精品一区二区三区| 久久人妻福利社区极品人妻图片| 法律面前人人平等表现在哪些方面| 岛国在线观看网站| 欧美日韩乱码在线| 美女 人体艺术 gogo| 欧美一级a爱片免费观看看 | av在线播放免费不卡| 精品免费久久久久久久清纯| 日本三级黄在线观看| 欧美精品亚洲一区二区| 88av欧美| 三级毛片av免费| 99国产精品99久久久久| 亚洲国产毛片av蜜桃av| 中国美女看黄片| 最新美女视频免费是黄的| 精品国产国语对白av| www.精华液| 久久伊人香网站| 长腿黑丝高跟| 国产单亲对白刺激| 午夜久久久在线观看| 国产单亲对白刺激| 日韩欧美国产在线观看| 99国产精品99久久久久| 国产乱人伦免费视频| 天天躁夜夜躁狠狠躁躁| 在线观看一区二区三区| 国产亚洲欧美在线一区二区| 日韩av在线大香蕉| 在线观看日韩欧美| 99热这里只有精品一区 | 97人妻精品一区二区三区麻豆 | 一边摸一边抽搐一进一小说| 最近最新中文字幕大全电影3 | 丁香欧美五月| 别揉我奶头~嗯~啊~动态视频| 丝袜人妻中文字幕| 免费在线观看亚洲国产| 国产亚洲欧美精品永久| 色播在线永久视频| 欧美三级亚洲精品| 在线观看一区二区三区| 国产欧美日韩精品亚洲av| 俄罗斯特黄特色一大片| 亚洲精品中文字幕一二三四区| 欧美精品啪啪一区二区三区| 精品国内亚洲2022精品成人| 国产精品久久久久久亚洲av鲁大| 久久亚洲真实| 日本一本二区三区精品| 午夜免费成人在线视频| 久久中文字幕一级| 亚洲天堂国产精品一区在线| 亚洲专区中文字幕在线| 丝袜美腿诱惑在线| 天堂影院成人在线观看| 亚洲精品在线美女| 嫩草影院精品99| 久久久国产成人精品二区| 岛国在线观看网站| 在线天堂中文资源库| 午夜老司机福利片| 欧美精品啪啪一区二区三区| 久久精品影院6| 无人区码免费观看不卡| 老司机在亚洲福利影院| 女性生殖器流出的白浆| 欧美另类亚洲清纯唯美| 动漫黄色视频在线观看| 免费观看人在逋| 欧美中文综合在线视频| 老熟妇仑乱视频hdxx| 国产精品av久久久久免费| 成人三级做爰电影| 悠悠久久av| 午夜福利免费观看在线| 国产久久久一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | 午夜久久久在线观看| 黄色成人免费大全| 久久久国产精品麻豆| 亚洲成国产人片在线观看| 午夜激情av网站| av在线播放免费不卡| 一本久久中文字幕| 十八禁人妻一区二区| 成人国产一区最新在线观看| 久久久久久大精品| 97人妻精品一区二区三区麻豆 | 狠狠狠狠99中文字幕| 国产极品粉嫩免费观看在线| 免费在线观看日本一区| 亚洲国产欧美日韩在线播放| 久久婷婷人人爽人人干人人爱| 一区二区三区高清视频在线| 午夜福利高清视频| 午夜影院日韩av| 日韩三级视频一区二区三区| 久久中文字幕一级| 男人操女人黄网站| 伦理电影免费视频| 天天一区二区日本电影三级| 每晚都被弄得嗷嗷叫到高潮| 草草在线视频免费看| 俄罗斯特黄特色一大片| 色精品久久人妻99蜜桃| 欧美黄色淫秽网站| 制服人妻中文乱码| 香蕉久久夜色| 亚洲精品一卡2卡三卡4卡5卡| 欧美+亚洲+日韩+国产| 成人免费观看视频高清| 亚洲成av人片免费观看| 欧洲精品卡2卡3卡4卡5卡区| 欧美精品啪啪一区二区三区| avwww免费| av视频在线观看入口| 麻豆成人午夜福利视频| 国内久久婷婷六月综合欲色啪| 中文字幕精品免费在线观看视频| 亚洲精品国产区一区二| 国产av一区二区精品久久| 村上凉子中文字幕在线| 给我免费播放毛片高清在线观看| 久久久精品国产亚洲av高清涩受| 日本三级黄在线观看| 免费在线观看日本一区| 动漫黄色视频在线观看| 精品国产美女av久久久久小说| av欧美777| 少妇裸体淫交视频免费看高清 | 巨乳人妻的诱惑在线观看| 亚洲九九香蕉| 成人国语在线视频| 国产精品电影一区二区三区| 老司机在亚洲福利影院| 中文字幕人成人乱码亚洲影| 18禁美女被吸乳视频| 好男人在线观看高清免费视频 | 嫩草影视91久久| 日韩精品免费视频一区二区三区| 亚洲熟妇中文字幕五十中出| 一夜夜www| 99久久无色码亚洲精品果冻| 在线观看免费午夜福利视频| 99久久精品国产亚洲精品| 美女扒开内裤让男人捅视频| 久久国产精品人妻蜜桃| 色播在线永久视频| 亚洲专区字幕在线| 波多野结衣高清无吗| 久久精品影院6| 欧美黑人巨大hd| 久99久视频精品免费| 免费在线观看日本一区| 日本黄色视频三级网站网址| 国产精品久久久久久精品电影 | 99国产精品99久久久久| 久久精品国产亚洲av高清一级| 可以在线观看的亚洲视频| 久久伊人香网站| 国产亚洲精品一区二区www| 国产亚洲精品久久久久久毛片| 色在线成人网| 亚洲中文av在线| 成年人黄色毛片网站| 成熟少妇高潮喷水视频| 久久精品国产综合久久久| 亚洲精品国产区一区二| 欧美黄色片欧美黄色片| 真人一进一出gif抽搐免费| 国产亚洲精品一区二区www| 久久婷婷人人爽人人干人人爱| 大型av网站在线播放| www.www免费av| 制服人妻中文乱码| 黑人巨大精品欧美一区二区mp4| 巨乳人妻的诱惑在线观看| 一二三四社区在线视频社区8| 最近最新免费中文字幕在线| 精品久久久久久成人av| 国产免费男女视频| 在线观看www视频免费| 色综合亚洲欧美另类图片| 少妇裸体淫交视频免费看高清 | 久热这里只有精品99| 99在线人妻在线中文字幕| av视频在线观看入口| 午夜激情福利司机影院| 免费观看精品视频网站| 国产一区二区三区视频了| 欧美黄色片欧美黄色片| 99国产精品一区二区蜜桃av| www国产在线视频色| 国产亚洲精品一区二区www| 国产亚洲精品第一综合不卡| 免费看美女性在线毛片视频| 日韩有码中文字幕| 国产成人精品久久二区二区免费| 国产爱豆传媒在线观看 | 成人三级黄色视频| 亚洲性夜色夜夜综合| 天堂√8在线中文| 中文字幕另类日韩欧美亚洲嫩草| 女人被狂操c到高潮| 激情在线观看视频在线高清| 精华霜和精华液先用哪个| 夜夜爽天天搞| 久99久视频精品免费| 久久人妻av系列| 一二三四在线观看免费中文在| 国产精品,欧美在线| 国产熟女午夜一区二区三区| 亚洲av电影不卡..在线观看| 美女大奶头视频| 午夜福利一区二区在线看| 天天添夜夜摸| 欧美色视频一区免费| 成人欧美大片| 国产不卡一卡二| 免费看美女性在线毛片视频| 国产成人av教育| 啦啦啦观看免费观看视频高清| 精品国产亚洲在线| 亚洲精品美女久久av网站| 日韩大码丰满熟妇| 久久精品国产99精品国产亚洲性色| 国产精品 欧美亚洲| 欧美精品亚洲一区二区| 十八禁网站免费在线| 亚洲一卡2卡3卡4卡5卡精品中文| www.精华液| 免费无遮挡裸体视频| 国产在线观看jvid| 免费观看人在逋| 欧美另类亚洲清纯唯美| 美女国产高潮福利片在线看| 久久精品夜夜夜夜夜久久蜜豆 | 神马国产精品三级电影在线观看 | 88av欧美| 欧美黑人巨大hd| 夜夜夜夜夜久久久久| 亚洲午夜精品一区,二区,三区| 国产蜜桃级精品一区二区三区| 给我免费播放毛片高清在线观看| 看片在线看免费视频| 国产又爽黄色视频| 亚洲欧美精品综合久久99| 女人高潮潮喷娇喘18禁视频| 三级毛片av免费| 悠悠久久av| 在线观看免费日韩欧美大片| 精品国产亚洲在线| 成年人黄色毛片网站| 精品不卡国产一区二区三区| 视频在线观看一区二区三区| 精品不卡国产一区二区三区| 日韩大码丰满熟妇| 亚洲三区欧美一区| 十八禁人妻一区二区| 国产精品一区二区精品视频观看| 色播亚洲综合网| 波多野结衣高清作品| 18禁国产床啪视频网站| 正在播放国产对白刺激| 男女下面进入的视频免费午夜 | 欧美亚洲日本最大视频资源| 久久久久九九精品影院| 国语自产精品视频在线第100页| 亚洲av第一区精品v没综合| 日韩大尺度精品在线看网址| www.999成人在线观看| 最新美女视频免费是黄的| 久久婷婷人人爽人人干人人爱| 成人三级做爰电影| 日韩精品免费视频一区二区三区| 成年女人毛片免费观看观看9| 人人妻,人人澡人人爽秒播| 欧美日韩福利视频一区二区| 亚洲精品中文字幕在线视频| 麻豆成人午夜福利视频| 日韩免费av在线播放| 后天国语完整版免费观看| 午夜两性在线视频| 青草久久国产| 熟女电影av网| 制服诱惑二区| 黄色 视频免费看| 91国产中文字幕| 中文亚洲av片在线观看爽| 99久久无色码亚洲精品果冻| 亚洲专区国产一区二区| 俺也久久电影网| 欧美成人午夜精品| 一区二区三区精品91| 一级a爱视频在线免费观看| 成人午夜高清在线视频 | 一本一本综合久久| 久久婷婷成人综合色麻豆| 一本一本综合久久| ponron亚洲| 久99久视频精品免费| 成人av一区二区三区在线看| 黄色片一级片一级黄色片| 日韩欧美 国产精品| 精品欧美国产一区二区三| 精品一区二区三区视频在线观看免费| 欧美成狂野欧美在线观看| 亚洲中文字幕日韩| 香蕉国产在线看| 国产一级毛片七仙女欲春2 | 国产男靠女视频免费网站| 国产麻豆成人av免费视频| 天天躁夜夜躁狠狠躁躁| xxxwww97欧美| 搡老岳熟女国产| 99久久综合精品五月天人人| 悠悠久久av| 91麻豆精品激情在线观看国产| 熟女少妇亚洲综合色aaa.| 精品国产一区二区三区四区第35| 精品午夜福利视频在线观看一区| 亚洲精品久久成人aⅴ小说| 成熟少妇高潮喷水视频| 国产蜜桃级精品一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 免费电影在线观看免费观看| 精品久久久久久久久久免费视频| 男女床上黄色一级片免费看| 成人18禁在线播放| 宅男免费午夜| 免费观看精品视频网站| 淫秽高清视频在线观看| 中文资源天堂在线| 国产一卡二卡三卡精品| 人人澡人人妻人| 亚洲精品色激情综合| 男女那种视频在线观看| 欧美激情极品国产一区二区三区| 91在线观看av| 亚洲av熟女| 国产精品98久久久久久宅男小说| 一区二区三区精品91| 不卡一级毛片| 中文字幕精品亚洲无线码一区 | 青草久久国产| 天天躁夜夜躁狠狠躁躁| 成人国语在线视频| 国产主播在线观看一区二区| 欧美大码av| 亚洲久久久国产精品| 精品国产乱子伦一区二区三区| 黄色视频不卡| 热re99久久国产66热| 国产精品av久久久久免费| 老司机深夜福利视频在线观看| 国产精品永久免费网站| 亚洲男人天堂网一区| 国内精品久久久久精免费| 欧美午夜高清在线| 精品久久久久久成人av| www.自偷自拍.com| 久久天堂一区二区三区四区| 人人妻人人澡欧美一区二区| 中文字幕人妻熟女乱码| 久久伊人香网站| 一区福利在线观看| 桃红色精品国产亚洲av| 婷婷亚洲欧美| 亚洲欧美一区二区三区黑人| 一本综合久久免费| 国产精品免费一区二区三区在线| 精品久久久久久成人av| 午夜影院日韩av| 老司机午夜十八禁免费视频| 成人国产一区最新在线观看| 可以在线观看的亚洲视频| 国产精品久久久人人做人人爽| 在线观看免费日韩欧美大片| 欧美乱色亚洲激情| 亚洲久久久国产精品| 欧洲精品卡2卡3卡4卡5卡区| 亚洲国产毛片av蜜桃av| 午夜激情福利司机影院| 不卡一级毛片| 免费看十八禁软件| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲va日本ⅴa欧美va伊人久久| 免费搜索国产男女视频| 在线观看www视频免费| 国产极品粉嫩免费观看在线| 美女免费视频网站| 最新美女视频免费是黄的| 19禁男女啪啪无遮挡网站| 91在线观看av| 国产精品 国内视频| 国产一级毛片七仙女欲春2 | 国产主播在线观看一区二区| 午夜免费鲁丝| 久久久水蜜桃国产精品网| 色婷婷久久久亚洲欧美| 欧美日本视频| 亚洲成av片中文字幕在线观看| 少妇的丰满在线观看| 精品少妇一区二区三区视频日本电影| 亚洲欧美一区二区三区黑人| 国产伦一二天堂av在线观看| 久久精品91无色码中文字幕| 国产成人欧美| av片东京热男人的天堂| 日本 av在线| 黄色a级毛片大全视频| 色综合站精品国产| 不卡av一区二区三区| 麻豆av在线久日| 久久久久久人人人人人| 国产v大片淫在线免费观看| 久久精品国产清高在天天线| 久久久久久久久中文| 国产成人精品无人区| 国产亚洲av高清不卡| 久久 成人 亚洲| 日韩欧美国产一区二区入口| 亚洲avbb在线观看| 一进一出好大好爽视频| 麻豆一二三区av精品| 欧美日韩乱码在线| 午夜两性在线视频| 三级毛片av免费| 亚洲国产看品久久| 黄色a级毛片大全视频| av在线天堂中文字幕| 成年版毛片免费区| 最近在线观看免费完整版| 婷婷精品国产亚洲av在线| 国产亚洲欧美精品永久| 国产精品二区激情视频| 国产伦在线观看视频一区| 精品久久久久久久末码| 久久久久久九九精品二区国产 | xxx96com| 久久中文字幕一级| 欧美日韩瑟瑟在线播放| 18禁美女被吸乳视频| 色婷婷久久久亚洲欧美| 啦啦啦 在线观看视频| 桃红色精品国产亚洲av| 日韩精品中文字幕看吧| 亚洲一码二码三码区别大吗| 国产乱人伦免费视频| 免费搜索国产男女视频| 在线十欧美十亚洲十日本专区| 亚洲九九香蕉| 色播在线永久视频| 亚洲人成77777在线视频| 黑人操中国人逼视频| 欧美激情高清一区二区三区| 亚洲国产日韩欧美精品在线观看 | 国产伦人伦偷精品视频| 国产熟女xx| 丝袜美腿诱惑在线| 亚洲专区字幕在线| 精品免费久久久久久久清纯| 少妇熟女aⅴ在线视频| 制服丝袜大香蕉在线| 黄色成人免费大全| 午夜老司机福利片| 色综合亚洲欧美另类图片| 人人妻人人看人人澡| 久久天躁狠狠躁夜夜2o2o| 韩国精品一区二区三区| 久久久久精品国产欧美久久久| 国产三级黄色录像| 成年版毛片免费区| 热99re8久久精品国产| 久久婷婷人人爽人人干人人爱| 精品久久久久久久人妻蜜臀av| 夜夜爽天天搞| 欧美又色又爽又黄视频| 丰满人妻熟妇乱又伦精品不卡| 午夜福利免费观看在线| 国产亚洲av高清不卡| e午夜精品久久久久久久| 国产精品,欧美在线| 黑人操中国人逼视频| 俺也久久电影网| 亚洲美女黄片视频| 黑人欧美特级aaaaaa片| 免费观看人在逋| 欧美亚洲日本最大视频资源| 亚洲av电影不卡..在线观看| 国产精品日韩av在线免费观看| 亚洲av五月六月丁香网| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲自拍偷在线| 熟女少妇亚洲综合色aaa.| 91麻豆精品激情在线观看国产| 国产成人av激情在线播放| 黄色视频不卡| 免费一级毛片在线播放高清视频| 久久久久久国产a免费观看| av电影中文网址| 成人午夜高清在线视频 | 变态另类成人亚洲欧美熟女| 18禁美女被吸乳视频| 亚洲国产精品合色在线| 男女视频在线观看网站免费 | 熟女电影av网| 日韩大码丰满熟妇| 亚洲电影在线观看av| 一边摸一边抽搐一进一小说| 一边摸一边做爽爽视频免费| 精品久久久久久久末码| 国产精品亚洲一级av第二区| 午夜福利视频1000在线观看| 黄色 视频免费看| 免费在线观看成人毛片| 久久久久亚洲av毛片大全| 少妇的丰满在线观看| 国产av一区二区精品久久| 精品久久久久久久久久免费视频| 亚洲第一av免费看| 国产精品永久免费网站| 国产91精品成人一区二区三区| 亚洲天堂国产精品一区在线| 亚洲成av片中文字幕在线观看| 久久精品国产综合久久久| 青草久久国产| 一a级毛片在线观看| 人人妻人人看人人澡| 大香蕉久久成人网| 美女大奶头视频| 免费一级毛片在线播放高清视频| 久久久久国内视频| 成人三级黄色视频| 久久 成人 亚洲| 日韩大码丰满熟妇| 美女免费视频网站| 中文字幕高清在线视频| 男人的好看免费观看在线视频 | 黄频高清免费视频| 亚洲专区中文字幕在线| 91大片在线观看| 亚洲七黄色美女视频| 嫩草影视91久久| 精品一区二区三区av网在线观看| 亚洲黑人精品在线| 国产三级在线视频| 国产精品电影一区二区三区| 这个男人来自地球电影免费观看| 一级a爱片免费观看的视频| 一级毛片精品| 桃红色精品国产亚洲av| 男人舔女人下体高潮全视频| 韩国av一区二区三区四区| 久久午夜亚洲精品久久| 草草在线视频免费看| 91在线观看av| 一级作爱视频免费观看| 久久精品国产清高在天天线| 757午夜福利合集在线观看| 免费在线观看亚洲国产| 长腿黑丝高跟| 国产亚洲av高清不卡|