• <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)化電容器配置研究
    免费av中文字幕在线| 欧美亚洲 丝袜 人妻 在线| 久久久久久免费高清国产稀缺| 中文字幕最新亚洲高清| 国产一卡二卡三卡精品| 三级毛片av免费| 丁香欧美五月| 中国美女看黄片| 久久久国产欧美日韩av| 身体一侧抽搐| 久久久精品免费免费高清| 亚洲精品国产一区二区精华液| 天堂动漫精品| 久久久精品区二区三区| 亚洲avbb在线观看| 午夜老司机福利片| 69精品国产乱码久久久| 深夜精品福利| 老熟妇乱子伦视频在线观看| 久久久国产成人精品二区 | 国产在线一区二区三区精| 午夜91福利影院| 成年女人毛片免费观看观看9 | 国产亚洲欧美98| 亚洲九九香蕉| 亚洲少妇的诱惑av| 久久国产亚洲av麻豆专区| 久久久久视频综合| 久久亚洲真实| 天天躁狠狠躁夜夜躁狠狠躁| 免费在线观看视频国产中文字幕亚洲| 正在播放国产对白刺激| 久久精品国产亚洲av高清一级| 在线视频色国产色| 少妇 在线观看| 久久ye,这里只有精品| 69精品国产乱码久久久| 亚洲一码二码三码区别大吗| 国产欧美日韩一区二区三| 国产亚洲一区二区精品| 国产成人免费无遮挡视频| 久久久久国产一级毛片高清牌| 高清黄色对白视频在线免费看| 一区二区日韩欧美中文字幕| 纯流量卡能插随身wifi吗| 欧美精品啪啪一区二区三区| 成人18禁在线播放| 别揉我奶头~嗯~啊~动态视频| 亚洲中文av在线| 国产免费av片在线观看野外av| 亚洲va日本ⅴa欧美va伊人久久| 国产精品一区二区免费欧美| av福利片在线| 51午夜福利影视在线观看| 日韩欧美一区视频在线观看| 欧美 日韩 精品 国产| 精品国产国语对白av| √禁漫天堂资源中文www| 18禁黄网站禁片午夜丰满| 老鸭窝网址在线观看| 久久久国产精品麻豆| 18禁观看日本| 欧美 亚洲 国产 日韩一| av一本久久久久| 久久久国产成人免费| 麻豆av在线久日| 精品人妻熟女毛片av久久网站| 日韩欧美国产一区二区入口| 欧美国产精品va在线观看不卡| 亚洲欧洲精品一区二区精品久久久| 国产成人av教育| 18禁裸乳无遮挡免费网站照片 | 丝袜美足系列| 久久久久精品国产欧美久久久| 18禁裸乳无遮挡动漫免费视频| cao死你这个sao货| 精品一品国产午夜福利视频| 亚洲精品一卡2卡三卡4卡5卡| 欧美中文综合在线视频| 99香蕉大伊视频| 狠狠狠狠99中文字幕| 黄色女人牲交| 在线观看一区二区三区激情| 12—13女人毛片做爰片一| 精品少妇一区二区三区视频日本电影| netflix在线观看网站| 成人18禁高潮啪啪吃奶动态图| 欧美日韩视频精品一区| 极品教师在线免费播放| 精品久久久久久,| 老汉色∧v一级毛片| 国产欧美亚洲国产| 黑人巨大精品欧美一区二区蜜桃| 国产视频一区二区在线看| 亚洲av成人一区二区三| 精品一区二区三区av网在线观看| 欧洲精品卡2卡3卡4卡5卡区| 99精品久久久久人妻精品| 一夜夜www| 欧美黑人欧美精品刺激| 午夜福利,免费看| 国产野战对白在线观看| 亚洲色图av天堂| 久久 成人 亚洲| 高清av免费在线| 欧美日韩一级在线毛片| 大型av网站在线播放| 亚洲精品粉嫩美女一区| 91成人精品电影| 女同久久另类99精品国产91| 美女福利国产在线| 欧美另类亚洲清纯唯美| 久久99一区二区三区| 日本精品一区二区三区蜜桃| 三上悠亚av全集在线观看| 天天躁日日躁夜夜躁夜夜| 日韩欧美在线二视频 | 成人18禁高潮啪啪吃奶动态图| 中文字幕人妻熟女乱码| 欧美激情极品国产一区二区三区| 亚洲国产欧美日韩在线播放| 色播在线永久视频| 大片电影免费在线观看免费| 一级,二级,三级黄色视频| 一本一本久久a久久精品综合妖精| 韩国精品一区二区三区| 一级黄色大片毛片| 亚洲专区字幕在线| 色尼玛亚洲综合影院| 日本wwww免费看| 高清毛片免费观看视频网站 | 99国产精品一区二区三区| 久久精品国产99精品国产亚洲性色 | 亚洲av成人一区二区三| 国产97色在线日韩免费| 欧美精品人与动牲交sv欧美| 男女床上黄色一级片免费看| 韩国精品一区二区三区| 成年人午夜在线观看视频| 看黄色毛片网站| 久久久久精品国产欧美久久久| 国产单亲对白刺激| 免费不卡黄色视频| 国产熟女午夜一区二区三区| 免费女性裸体啪啪无遮挡网站| 捣出白浆h1v1| 久久精品aⅴ一区二区三区四区| 国产高清视频在线播放一区| 国产精品久久久久久精品古装| a级毛片黄视频| 美女福利国产在线| 男女高潮啪啪啪动态图| 国产亚洲欧美98| 99热网站在线观看| 最近最新免费中文字幕在线| 欧美日韩乱码在线| 国产精品一区二区在线观看99| 欧美激情 高清一区二区三区| 国产精品影院久久| av免费在线观看网站| 中文字幕精品免费在线观看视频| 亚洲欧美日韩高清在线视频| 午夜老司机福利片| 91国产中文字幕| 视频在线观看一区二区三区| 亚洲精华国产精华精| 亚洲精品中文字幕一二三四区| 狠狠婷婷综合久久久久久88av| 99在线人妻在线中文字幕 | 国产成人免费观看mmmm| 亚洲色图 男人天堂 中文字幕| 两性夫妻黄色片| 嫁个100分男人电影在线观看| 国产色视频综合| 久99久视频精品免费| 99精品久久久久人妻精品| av天堂久久9| 国产av又大| 亚洲精品中文字幕一二三四区| 成在线人永久免费视频| 露出奶头的视频| 一进一出好大好爽视频| 亚洲熟妇熟女久久| 黄色女人牲交| 国产精品香港三级国产av潘金莲| 国产野战对白在线观看| 女人高潮潮喷娇喘18禁视频| 婷婷成人精品国产| 国产亚洲精品第一综合不卡| 国产野战对白在线观看| 久久精品国产a三级三级三级| 每晚都被弄得嗷嗷叫到高潮| 亚洲全国av大片| 18禁观看日本| 国产精品久久久久成人av| 亚洲第一av免费看| 他把我摸到了高潮在线观看| 亚洲中文av在线| 王馨瑶露胸无遮挡在线观看| 久久草成人影院| 国产精品美女特级片免费视频播放器 | 国产精品国产高清国产av | 久久精品国产a三级三级三级| ponron亚洲| 美女视频免费永久观看网站| 热re99久久精品国产66热6| 下体分泌物呈黄色| 一边摸一边做爽爽视频免费| 精品国产超薄肉色丝袜足j| 成人av一区二区三区在线看| 亚洲精品在线美女| 亚洲精品在线美女| 欧美 日韩 精品 国产| 亚洲精品在线观看二区| 男人舔女人的私密视频| 热re99久久精品国产66热6| 黄网站色视频无遮挡免费观看| 国产精品乱码一区二三区的特点 | 亚洲片人在线观看| 午夜视频精品福利| 亚洲专区中文字幕在线| 成人国产一区最新在线观看| 日本欧美视频一区| 99久久人妻综合| 亚洲中文字幕日韩| 精品一区二区三区视频在线观看免费 | 女人高潮潮喷娇喘18禁视频| 国产深夜福利视频在线观看| 午夜视频精品福利| 丰满饥渴人妻一区二区三| 亚洲中文字幕日韩| 不卡一级毛片| 不卡av一区二区三区| 精品国产国语对白av| 在线国产一区二区在线| 精品少妇久久久久久888优播| 久久青草综合色| 亚洲一区二区三区不卡视频| 777久久人妻少妇嫩草av网站| 岛国在线观看网站| 午夜激情av网站| 欧美日韩精品网址| 又黄又爽又免费观看的视频| 国产无遮挡羞羞视频在线观看| 日本五十路高清| av在线播放免费不卡| 久久久久精品人妻al黑| 亚洲精品美女久久久久99蜜臀| videosex国产| 亚洲av电影在线进入| cao死你这个sao货| 亚洲精品中文字幕一二三四区| 啪啪无遮挡十八禁网站| 亚洲精华国产精华精| 久久性视频一级片| bbb黄色大片| 国产97色在线日韩免费| 久久这里只有精品19| 日本vs欧美在线观看视频| 老司机午夜福利在线观看视频| 99久久精品国产亚洲精品| 亚洲精品国产一区二区精华液| 国产精品欧美亚洲77777| 黑人猛操日本美女一级片| 亚洲熟妇中文字幕五十中出 | 午夜精品国产一区二区电影| 自线自在国产av| 国产又爽黄色视频| 亚洲av美国av| 精品国产乱码久久久久久男人| x7x7x7水蜜桃| 看免费av毛片| 9191精品国产免费久久| 亚洲熟女毛片儿| 村上凉子中文字幕在线| 国产淫语在线视频| 在线永久观看黄色视频| 91精品三级在线观看| 一个人免费在线观看的高清视频| a级毛片在线看网站| 亚洲色图 男人天堂 中文字幕| 在线av久久热| 久久人妻av系列| 一本综合久久免费| 国产亚洲av高清不卡| 国产在线一区二区三区精| 啦啦啦免费观看视频1| 黄色视频不卡| 一级片免费观看大全| 久久久精品国产亚洲av高清涩受| 欧美成人免费av一区二区三区 | 王馨瑶露胸无遮挡在线观看| a级毛片黄视频| 黄色 视频免费看| av在线播放免费不卡| 久久久久精品人妻al黑| 亚洲av熟女| 国产国语露脸激情在线看| 狂野欧美激情性xxxx| 91大片在线观看| 黄色毛片三级朝国网站| 国产精品免费视频内射| 久久人人97超碰香蕉20202| 国产欧美日韩一区二区三| 欧美日本中文国产一区发布| 国产成人精品无人区| 中文欧美无线码| 亚洲成国产人片在线观看| 国产人伦9x9x在线观看| 脱女人内裤的视频| 免费日韩欧美在线观看| 老熟女久久久| 国产精品国产av在线观看| 欧美激情高清一区二区三区| 欧美精品人与动牲交sv欧美| 中文字幕色久视频| 91大片在线观看| 在线观看舔阴道视频| 少妇 在线观看| 久久久久久久国产电影| 国产精品香港三级国产av潘金莲| 日韩三级视频一区二区三区| 国产av精品麻豆| 男人的好看免费观看在线视频 | 国产精品1区2区在线观看. | 午夜福利乱码中文字幕| 欧美老熟妇乱子伦牲交| 黄色视频,在线免费观看| 十八禁人妻一区二区| 亚洲免费av在线视频| 国产片内射在线| 欧美日韩精品网址| 亚洲av成人不卡在线观看播放网| 精品一品国产午夜福利视频| 午夜福利免费观看在线| 欧美亚洲 丝袜 人妻 在线| 99久久精品国产亚洲精品| 久久青草综合色| 亚洲中文字幕日韩| 男人舔女人的私密视频| 高清视频免费观看一区二区| 国产成人精品无人区| 精品无人区乱码1区二区| 亚洲精品国产色婷婷电影| av电影中文网址| 国产精品久久久久久精品古装| 国产欧美日韩一区二区三| 一个人免费在线观看的高清视频| 这个男人来自地球电影免费观看| 免费在线观看黄色视频的| 久久人人爽av亚洲精品天堂| 欧美精品啪啪一区二区三区| 久久婷婷成人综合色麻豆| 日韩 欧美 亚洲 中文字幕| 日韩三级视频一区二区三区| 国产成人影院久久av| 日韩一卡2卡3卡4卡2021年| 制服人妻中文乱码| 色老头精品视频在线观看| 精品电影一区二区在线| 久久中文看片网| 69av精品久久久久久| 在线av久久热| 看免费av毛片| 老司机靠b影院| 国产精品免费一区二区三区在线 | 91av网站免费观看| 久久精品91无色码中文字幕| 极品人妻少妇av视频| 午夜成年电影在线免费观看| 久久 成人 亚洲| 一二三四社区在线视频社区8| 五月开心婷婷网| av网站免费在线观看视频| 成人三级做爰电影| 久久国产亚洲av麻豆专区| 国产乱人伦免费视频| 99久久人妻综合| 亚洲久久久国产精品| 亚洲黑人精品在线| 怎么达到女性高潮| 精品电影一区二区在线| 欧美黄色片欧美黄色片| 黄片大片在线免费观看| 69精品国产乱码久久久| 正在播放国产对白刺激| 一级毛片精品| 亚洲国产欧美日韩在线播放| 99热只有精品国产| 久久精品成人免费网站| 成人影院久久| 国产视频一区二区在线看| 两个人免费观看高清视频| 黄色视频,在线免费观看| 极品人妻少妇av视频| 黄色怎么调成土黄色| 国产无遮挡羞羞视频在线观看| 久久久水蜜桃国产精品网| 久久香蕉激情| 久久国产亚洲av麻豆专区| 国产高清视频在线播放一区| 成年人免费黄色播放视频| 精品国产一区二区久久| 久久精品国产综合久久久| 丝瓜视频免费看黄片| 午夜福利在线免费观看网站| 一级片'在线观看视频| 97人妻天天添夜夜摸| 人人妻,人人澡人人爽秒播| 久久人妻熟女aⅴ| 午夜久久久在线观看| 丰满饥渴人妻一区二区三| 19禁男女啪啪无遮挡网站| 很黄的视频免费| av不卡在线播放| 国产一卡二卡三卡精品| 如日韩欧美国产精品一区二区三区| 国产高清videossex| 老司机靠b影院| 久久人妻福利社区极品人妻图片| 熟女少妇亚洲综合色aaa.| 麻豆成人av在线观看| 黄色视频不卡| 亚洲免费av在线视频| 欧美精品啪啪一区二区三区| 美女高潮喷水抽搐中文字幕| bbb黄色大片| 在线十欧美十亚洲十日本专区| 在线永久观看黄色视频| 免费在线观看黄色视频的| 色婷婷久久久亚洲欧美| 久久中文看片网| x7x7x7水蜜桃| 黄色成人免费大全| 五月开心婷婷网| 在线av久久热| 天堂俺去俺来也www色官网| 视频区欧美日本亚洲| 国产亚洲欧美精品永久| 一边摸一边做爽爽视频免费| 另类亚洲欧美激情| 丝袜美腿诱惑在线| 19禁男女啪啪无遮挡网站| 国产91精品成人一区二区三区| 好男人电影高清在线观看| av有码第一页| www.精华液| 亚洲全国av大片| 麻豆av在线久日| 搡老熟女国产l中国老女人| 久久亚洲精品不卡| 亚洲人成电影免费在线| 欧美日韩av久久| 国产欧美日韩精品亚洲av| av网站在线播放免费| 12—13女人毛片做爰片一| 亚洲七黄色美女视频| 一级毛片精品| 一区二区日韩欧美中文字幕| 99久久国产精品久久久| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av日韩精品久久久久久密| √禁漫天堂资源中文www| www.自偷自拍.com| 在线观看午夜福利视频| 香蕉久久夜色| 亚洲熟妇熟女久久| 成人三级做爰电影| 日韩欧美国产一区二区入口| 欧美亚洲 丝袜 人妻 在线| 亚洲一区二区三区不卡视频| 日本精品一区二区三区蜜桃| 国产成人欧美在线观看 | 欧美在线一区亚洲| 亚洲av片天天在线观看| 国产主播在线观看一区二区| 黄色片一级片一级黄色片| 又紧又爽又黄一区二区| 大陆偷拍与自拍| 精品国产乱子伦一区二区三区| 欧美av亚洲av综合av国产av| 欧美日韩一级在线毛片| 黄色成人免费大全| 99国产综合亚洲精品| 人人澡人人妻人| 欧美黑人精品巨大| 亚洲一区高清亚洲精品| 美女视频免费永久观看网站| 一级片免费观看大全| 成人永久免费在线观看视频| 捣出白浆h1v1| av一本久久久久| 日韩成人在线观看一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 变态另类成人亚洲欧美熟女 | 久久久久国产精品人妻aⅴ院 | 啦啦啦在线免费观看视频4| 国产精品国产av在线观看| 婷婷丁香在线五月| 91字幕亚洲| 男人的好看免费观看在线视频 | 一级毛片女人18水好多| 亚洲国产看品久久| 老司机午夜十八禁免费视频| ponron亚洲| 欧美亚洲 丝袜 人妻 在线| 老熟妇乱子伦视频在线观看| 久久中文字幕一级| 少妇猛男粗大的猛烈进出视频| av视频免费观看在线观看| 久久精品国产99精品国产亚洲性色 | 亚洲成人免费电影在线观看| 成人免费观看视频高清| 两个人看的免费小视频| 国产淫语在线视频| 女人被狂操c到高潮| 男人舔女人的私密视频| 欧美乱妇无乱码| 精品少妇一区二区三区视频日本电影| 精品国产乱码久久久久久男人| 国产1区2区3区精品| 免费在线观看完整版高清| 色尼玛亚洲综合影院| 一区二区日韩欧美中文字幕| 久久精品国产99精品国产亚洲性色 | 国产成人欧美| 亚洲五月婷婷丁香| 美女国产高潮福利片在线看| 91av网站免费观看| 真人做人爱边吃奶动态| 十八禁网站免费在线| 久久影院123| 国产精品久久久人人做人人爽| 国产区一区二久久| 好男人电影高清在线观看| 国产色视频综合| 成人特级黄色片久久久久久久| 国产精品久久视频播放| 手机成人av网站| 狠狠婷婷综合久久久久久88av| 99久久99久久久精品蜜桃| 国产一区二区激情短视频| 亚洲精品久久午夜乱码| 精品乱码久久久久久99久播| 亚洲国产中文字幕在线视频| 在线观看一区二区三区激情| 日本撒尿小便嘘嘘汇集6| 99久久人妻综合| 人人妻人人澡人人看| 18禁观看日本| 操出白浆在线播放| 精品电影一区二区在线| 99香蕉大伊视频| 老司机靠b影院| 亚洲熟女精品中文字幕| 亚洲精品国产区一区二| 久久九九热精品免费| 精品国产超薄肉色丝袜足j| 嫁个100分男人电影在线观看| 国产亚洲av高清不卡| 十八禁人妻一区二区| 亚洲精品一卡2卡三卡4卡5卡| 人成视频在线观看免费观看| 老司机亚洲免费影院| 18禁美女被吸乳视频| 色婷婷久久久亚洲欧美| 最近最新中文字幕大全免费视频| 成人永久免费在线观看视频| 香蕉丝袜av| 中文字幕精品免费在线观看视频| 免费观看人在逋| 老司机深夜福利视频在线观看| 少妇裸体淫交视频免费看高清 | 国产不卡av网站在线观看| 高清在线国产一区| 一进一出好大好爽视频| 成人影院久久| 人妻一区二区av| 国产99久久九九免费精品| 国产一区二区三区综合在线观看| 精品一区二区三区av网在线观看| 亚洲av成人av| 中文字幕高清在线视频| 亚洲av电影在线进入| 一本综合久久免费| 久久99一区二区三区| 国产精品98久久久久久宅男小说| 久久亚洲真实| 久久香蕉精品热| 欧美日韩视频精品一区| 三上悠亚av全集在线观看| 日本一区二区免费在线视频| 国产在线一区二区三区精| 99久久精品国产亚洲精品| 免费黄频网站在线观看国产| 成人三级做爰电影| 国产精品久久久人人做人人爽| 久久国产精品人妻蜜桃| 久久中文看片网| 国产99久久九九免费精品| 精品久久久久久久毛片微露脸| 好看av亚洲va欧美ⅴa在| 亚洲国产精品一区二区三区在线| 王馨瑶露胸无遮挡在线观看| 免费观看人在逋| 手机成人av网站| 大香蕉久久成人网| 婷婷成人精品国产| 一级作爱视频免费观看| 国产野战对白在线观看| 人妻久久中文字幕网| 80岁老熟妇乱子伦牲交| 人人澡人人妻人| 视频区欧美日本亚洲| 高清欧美精品videossex|