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

    基于一維深度卷積生成對(duì)抗網(wǎng)絡(luò)的鋼軌波磨識(shí)別方法

    2022-02-07 03:56:42謝燁趙聞強(qiáng)楊紅運(yùn)包學(xué)海
    鐵道建筑 2022年12期
    關(guān)鍵詞:波磨鋼軌幅值

    謝燁 趙聞強(qiáng) 楊紅運(yùn) 包學(xué)海

    1.浙江省交通投資集團(tuán)有限公司 智慧交通研究分公司,杭州 310016;2.浙江杭海城際鐵路有限公司,浙江 海寧 314499

    軌道交通列車具有載重量大、運(yùn)行密度高等特點(diǎn),頻繁的輪軌相互作用會(huì)導(dǎo)致鋼軌上出現(xiàn)波浪形的磨損,即鋼軌波磨。車輪在經(jīng)過(guò)鋼軌波磨區(qū)域時(shí)引發(fā)的強(qiáng)烈振動(dòng)和噪聲會(huì)大大加快鋼軌和軸箱內(nèi)零件的疲勞損壞速度[1-3]。波磨出現(xiàn)前后駕駛室內(nèi)噪聲波動(dòng)幅值可達(dá)到10 dB 左右,波磨造成的振動(dòng)還會(huì)導(dǎo)致軸向吊耳斷裂、鋼軌扣件彈條斷裂等問題[4-6],大大增加列車運(yùn)行維護(hù)成本。所以,對(duì)鋼軌波磨進(jìn)行識(shí)別具有十分重要的意義。

    由于鋼軌發(fā)生波磨位置不確定,雖然使用波磨測(cè)量?jī)x進(jìn)行檢測(cè)具有很高的準(zhǔn)確性,但需要投入大量的時(shí)間和人工成本。振動(dòng)監(jiān)測(cè)是一種在線狀態(tài)監(jiān)測(cè)方法,具有直接、實(shí)時(shí)的優(yōu)點(diǎn),并且設(shè)置后無(wú)需人工干涉即可對(duì)數(shù)據(jù)進(jìn)行采集,目前已廣泛應(yīng)用于各種設(shè)備的狀態(tài)監(jiān)測(cè)中[7-8]。將振動(dòng)監(jiān)測(cè)應(yīng)用于鋼軌波磨監(jiān)測(cè)中,可根據(jù)輪軌接觸產(chǎn)生的振動(dòng)直觀地觀察到振動(dòng)狀態(tài)的變化;根據(jù)振動(dòng)特征制定辨識(shí)模型,在鋼軌波磨的狀態(tài)監(jiān)測(cè)上更具即時(shí)性。目前許多學(xué)者基于振動(dòng)響應(yīng)對(duì)鋼軌波磨的狀態(tài)監(jiān)測(cè)進(jìn)行了研究。文獻(xiàn)[9]對(duì)鋼軌波磨振動(dòng)信號(hào)進(jìn)行降噪處理,然后使用希爾伯特-黃變換(Hilbert?Huang Transform,HHT)對(duì)信號(hào)進(jìn)行時(shí)頻分析,得到了較為準(zhǔn)確的分析結(jié)果。文獻(xiàn)[10]通過(guò)仿真模型得到鋼軌波磨振動(dòng)仿真信號(hào),計(jì)算小波包能量熵,實(shí)現(xiàn)了較為準(zhǔn)確的鋼軌波磨檢測(cè)。文獻(xiàn)[11]通過(guò)使用時(shí)間-能量信息的波磨指數(shù)和能量因子,設(shè)計(jì)了一種雙指標(biāo)聯(lián)合鋼軌波磨評(píng)價(jià)方法,可以對(duì)鋼軌波磨的波長(zhǎng)和幅值進(jìn)行評(píng)價(jià),并通過(guò)實(shí)測(cè)鋼軌波磨信號(hào)驗(yàn)證了其方法的準(zhǔn)確性。

    隨著大數(shù)據(jù)技術(shù)的發(fā)展與普及,智能狀態(tài)監(jiān)測(cè)與診斷成了新的研究熱點(diǎn),使用設(shè)備的振動(dòng)信號(hào)訓(xùn)練分類模型并對(duì)其狀態(tài)進(jìn)行分類,已經(jīng)廣泛應(yīng)用于滾動(dòng)軸承[12-14]、齒輪箱[15-17]和發(fā)動(dòng)機(jī)[18-20]等機(jī)械設(shè)備的狀態(tài)監(jiān)測(cè)與故障診斷。已有學(xué)者將智能診斷模型引入鋼軌波磨的識(shí)別研究中。文獻(xiàn)[21]對(duì)地鐵鋼軌波磨振動(dòng)信號(hào)進(jìn)行空間域分割后,使用一維卷積神經(jīng)網(wǎng)絡(luò)進(jìn)行了波磨識(shí)別與分類,并使用實(shí)測(cè)鋼軌波磨振動(dòng)信號(hào)進(jìn)行試驗(yàn),得到了很高的鋼軌波磨識(shí)別率。文獻(xiàn)[22-23]通過(guò)經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition,EMD),提取包含故障的本征模函數(shù)(Intrinsic Modal Functions,IMF)的能量和峭度,將其多尺度能量和峭度特征輸入神經(jīng)網(wǎng)絡(luò)進(jìn)行波磨識(shí)別。

    針對(duì)現(xiàn)有研究方法中存在的數(shù)據(jù)模型與真實(shí)鋼軌波磨數(shù)據(jù)有較大偏差,即波磨故障識(shí)別準(zhǔn)確率較低問題,本文提出一種基于一維深度卷積生成對(duì)抗網(wǎng)絡(luò)的鋼軌波磨識(shí)別方法。首先使用一維生成對(duì)抗網(wǎng)絡(luò)生成與實(shí)際波磨振動(dòng)信號(hào)結(jié)構(gòu)相同的偽樣本,對(duì)樣本信號(hào)的數(shù)據(jù)集進(jìn)行擴(kuò)充;然后提取波磨振動(dòng)信號(hào)的時(shí)域統(tǒng)計(jì)指標(biāo)作為波磨數(shù)據(jù)的特征;最后使用分類算法對(duì)不同特征的波磨振動(dòng)數(shù)據(jù)進(jìn)行學(xué)習(xí)與分類,以實(shí)現(xiàn)準(zhǔn)確地提取實(shí)際波磨數(shù)據(jù)故障特征,識(shí)別鋼軌波磨故障的目標(biāo)。

    1 生成對(duì)抗網(wǎng)絡(luò)基本原理

    生成對(duì)抗網(wǎng)絡(luò)(Generative Adversarial Networks,GAN)是一種生成模型,由Goodfellow在2014年首次提出[24],由生成網(wǎng)絡(luò)G和判別網(wǎng)絡(luò)D構(gòu)成。G通過(guò)學(xué)習(xí)真實(shí)樣本的結(jié)構(gòu),生成與真實(shí)樣本結(jié)構(gòu)相似的偽樣本;D用來(lái)分辨輸入數(shù)據(jù)是真實(shí)樣本還是偽樣本。G與D通過(guò)博弈,彼此不斷學(xué)習(xí)提升。直到D無(wú)法判斷輸入數(shù)據(jù)是真實(shí)樣本還是偽樣本時(shí),說(shuō)明G已經(jīng)學(xué)習(xí)到真實(shí)樣本數(shù)據(jù)結(jié)構(gòu)。

    GAN的訓(xùn)練方法如下。

    1)高斯噪聲通過(guò)G生成偽樣本,將偽樣本和真實(shí)樣本輸入最后一層激活函數(shù)為sigmoid 函數(shù)的D,真實(shí)樣本的標(biāo)簽設(shè)定為1,偽樣本的標(biāo)簽設(shè)定為0,轉(zhuǎn)化為一個(gè)二分類問題。使G的權(quán)值不變,通過(guò)反向傳播方法更新D中的權(quán)值,使D的判斷能力提高,用交叉熵作為其損失函數(shù)J(D),表達(dá)式為

    式中:y為樣本標(biāo)簽;x為真實(shí)樣本;G(z)為生成的偽樣本,z為輸入G的高斯噪聲;D(x)和D[G(z)]為真實(shí)樣本和偽樣本通過(guò)判別網(wǎng)絡(luò)的輸出。

    2)將偽樣本輸入D,并將偽樣本的標(biāo)簽記為1,通過(guò)反向傳播方法更新G的權(quán)值,使D將偽樣本判斷為真實(shí)樣本,從而使G可以學(xué)習(xí)真實(shí)樣本數(shù)據(jù)的結(jié)構(gòu)。同樣使用交叉熵作為損失函數(shù)J(G),表達(dá)式為

    3)不斷進(jìn)行迭代交替訓(xùn)練。G的目標(biāo)是生成足夠真的偽樣本,使得D無(wú)法準(zhǔn)確判斷;D的目標(biāo)是盡量準(zhǔn)確判斷真實(shí)樣本和偽樣本。G和D的損失函數(shù)可以等價(jià)為

    式中:Ex~Pdata為x的熵,Pdata為x的概率分布;Ez~Pz為z的熵,Pz為z的概率分布。

    總體損失函數(shù)V(D,G)表達(dá)式為

    網(wǎng)絡(luò)訓(xùn)練的最終理想結(jié)果是使G和D達(dá)到納什均衡,D無(wú)法判別G生成的偽樣本是否為真實(shí)樣本,即D[G(z)]=0.5。深度卷積生成對(duì)抗網(wǎng)絡(luò)(Deep Convolutional Generative Adversarial Networks,DCGAN)是一種基于GAN 的改進(jìn)算法[25]。該算法是將G和D均設(shè)計(jì)為卷積神經(jīng)網(wǎng)絡(luò),以卷積神經(jīng)網(wǎng)絡(luò)強(qiáng)大的特征提取能力,提高GAN 的效果?;诖耍疚氖褂?D?DCGAN 對(duì)一維鋼軌波磨振動(dòng)信號(hào)進(jìn)行偽樣本生成。

    2 基于1D?DCGAN 的鋼軌波磨識(shí)別方法

    DCGAN 最初被應(yīng)用于圖像生成。對(duì)于一維波磨振動(dòng)信號(hào),設(shè)計(jì)合理的1D?DCGAN 有助于網(wǎng)絡(luò)提取波磨振動(dòng)信號(hào)的特征和學(xué)習(xí)波磨振動(dòng)數(shù)據(jù)的結(jié)構(gòu)?;?D?DCGAN 的鋼軌波磨識(shí)別方法分為以下四個(gè)步驟。

    1)波磨振動(dòng)數(shù)據(jù)預(yù)處理;

    2)波磨振動(dòng)數(shù)據(jù)偽樣本生成;

    3)波磨振動(dòng)數(shù)據(jù)時(shí)域特征提取及降維;

    4)分類模型的訓(xùn)練與鋼軌波磨狀態(tài)識(shí)別。

    2.1 波磨振動(dòng)數(shù)據(jù)預(yù)處理

    從實(shí)測(cè)數(shù)據(jù)中篩選出質(zhì)量較好的波磨振動(dòng)數(shù)據(jù),將50%的數(shù)據(jù)作為1D?DCGAN 模型訓(xùn)練數(shù)據(jù),剩余50%的數(shù)據(jù)作為波磨識(shí)別測(cè)試數(shù)據(jù)。將訓(xùn)練數(shù)據(jù)進(jìn)行歸一化處理,然后輸入到設(shè)計(jì)好的1D?DCGAN 模型中。

    2.2 1D?DCGAN的模型結(jié)構(gòu)設(shè)計(jì)

    1D?DCGAN 網(wǎng)絡(luò)參數(shù)見表1。模型中,以高斯噪聲z作為G的輸入,經(jīng)過(guò)反卷積和激活函數(shù)輸出G(z),將G(z)和x隨機(jī)輸入D,經(jīng)過(guò)多層卷積層計(jì)算,輸出標(biāo)簽為0或1,其中0代表偽樣本,1代表真實(shí)樣本。為了將輸出層得到的偽樣本數(shù)值限制在[-1,1],使用tanh作為生成網(wǎng)絡(luò)輸出層的激活函數(shù);為了使判別網(wǎng)絡(luò)的輸出值為0 或1,使用sigmoid 作為判別網(wǎng)絡(luò)的輸出層激活函數(shù)。在GAN 的訓(xùn)練過(guò)程中,由于判別網(wǎng)絡(luò)和生成網(wǎng)絡(luò)訓(xùn)練不平衡,容易出現(xiàn)生成網(wǎng)絡(luò)梯度消失。為減少這一影響,每更新一次判別網(wǎng)絡(luò)參數(shù),訓(xùn)練四次生成網(wǎng)絡(luò)。

    表1 1D?DCGAN網(wǎng)絡(luò)參數(shù)

    2.3 時(shí)域統(tǒng)計(jì)指標(biāo)的特征提取及降維

    對(duì)實(shí)測(cè)鋼軌波磨振動(dòng)數(shù)據(jù)與1D?DCGAN 生成的偽樣本進(jìn)行時(shí)域統(tǒng)計(jì)指標(biāo)特征提取,包括均值、均方根值、方根幅值、絕對(duì)平均值、偏斜度、峭度、方差、最大值、最小值、波形指標(biāo)、峰值指標(biāo)、脈沖指標(biāo)、裕度指標(biāo)、偏斜度指標(biāo)、峭度指標(biāo)共15種時(shí)域特征統(tǒng)計(jì)指標(biāo)。通過(guò)主成分分析(Principal Component Analysis,PCA)降維后,前3 個(gè)主成分分量的累計(jì)貢獻(xiàn)率就可以達(dá)到95%以上。為了便于特征可視化和減少分類模型的計(jì)算量,提取前3 個(gè)主成分作為鋼軌波磨振動(dòng)數(shù)據(jù)的特征。

    2.4 分類模型的訓(xùn)練與鋼軌波磨狀態(tài)識(shí)別

    分別使用K 近鄰算法(K?Nearest Neighbor,KNN)、支持向量機(jī)(Support Vector Machine,SVM)和決策樹(Decision Tree,DT)三種不同的分類模型對(duì)鋼軌波磨特征進(jìn)行分類。為了模擬不同程度的鋼軌波磨狀態(tài),將振動(dòng)信號(hào)的幅值A(chǔ)分別縮小50%和增大50%。

    在鋼軌波磨狀態(tài)識(shí)別過(guò)程中,將正常鋼軌的標(biāo)簽設(shè)定為0,幅值為0.5A、1.0A、1.5A的鋼軌波磨數(shù)據(jù)標(biāo)簽分別設(shè)定為1、2、3。將測(cè)試數(shù)據(jù)輸入到訓(xùn)練好的分類模型中,對(duì)正常鋼軌與不同波磨程度的鋼軌進(jìn)行狀態(tài)識(shí)別。

    3 試驗(yàn)驗(yàn)證

    3.1 數(shù)據(jù)概況

    試驗(yàn)數(shù)據(jù)選用杭海(杭州—海寧)城際鐵路的列車實(shí)測(cè)軸箱振動(dòng)信號(hào)。杭海城際鐵路于2021 年6 月開始進(jìn)行初期運(yùn)營(yíng),全長(zhǎng)46.38 km,全線共12個(gè)站點(diǎn),列車采用DC1500V架空接觸網(wǎng)授電的B2型車,2動(dòng)車 +2 半動(dòng)車,4 輛編組,運(yùn)行速度為80 km/h。使用振動(dòng)采集器對(duì)列車軸箱的振動(dòng)信號(hào)進(jìn)行采集,傳感器的布置見圖1。采樣頻率為5 120 Hz,每條數(shù)據(jù)為1 000 個(gè)采樣點(diǎn)。對(duì)獲取的軸箱振動(dòng)信號(hào)進(jìn)行篩選,共獲得400條品質(zhì)較好的鋼軌波磨樣本數(shù)據(jù)。其中200條數(shù)據(jù)作為訓(xùn)練數(shù)據(jù),輸入到1D?DCGAN 中進(jìn)行模型訓(xùn)練;剩余200 條數(shù)據(jù)作為測(cè)試樣本,用作最后的鋼軌波磨識(shí)別。

    圖1 軸箱振動(dòng)測(cè)試中傳感器的布置

    3.2 數(shù)據(jù)生成

    使用所設(shè)計(jì)的1D?DCGAN 模型進(jìn)行數(shù)據(jù)生成。首先將200條真實(shí)的鋼軌波磨樣本數(shù)據(jù)進(jìn)行歸一化處理,然后輸入到模型中進(jìn)行訓(xùn)練。其中,優(yōu)化器使用Adam,學(xué)習(xí)率設(shè)置為0.000 01,迭代次數(shù)100 萬(wàn)次。為了得到最好的模型,每迭代5 萬(wàn)次對(duì)當(dāng)前模型進(jìn)行保存,然后繼續(xù)迭代。迭代完成后,挑選判別器對(duì)真實(shí)樣本和偽樣本的判別準(zhǔn)確率都在50%左右的模型,這個(gè)模型就是當(dāng)前最優(yōu)模型。使用binary crossentropy作為損失函數(shù),batch_size 設(shè)定為10,padding 選擇same。訓(xùn)練完成后,將1D?DCGAN 生成的數(shù)據(jù)進(jìn)行反歸一化,以得到與真實(shí)樣本相同的幅值。

    通過(guò)訓(xùn)練好的網(wǎng)絡(luò)生成偽樣本,與真實(shí)樣本進(jìn)行對(duì)比,見圖2??芍涸跁r(shí)域上,偽樣本與真實(shí)樣本的幅值基本一致,且結(jié)構(gòu)表現(xiàn)較為一致;在頻域上,偽樣本與真實(shí)樣本均在500 Hz和800 Hz附近幅值較大,二者主要頻段基本一致。

    圖2 真實(shí)樣本與偽樣本的時(shí)頻曲線對(duì)比

    對(duì)真實(shí)樣本、1D?DCGAN 生成的偽樣本及無(wú)損鋼軌數(shù)據(jù),計(jì)算15 個(gè)時(shí)域特征指標(biāo),使用PCA 降維后提取前3 個(gè)主成分分量作為其特征,見圖3??芍簜螛颖九c真實(shí)樣本的特征基本相同,并且與正常鋼軌提取的特征區(qū)別較大。

    圖3 真實(shí)樣本、偽樣本與無(wú)損鋼軌波磨數(shù)據(jù)特征對(duì)比

    3.3 數(shù)據(jù)集的建立

    挑選4 000 個(gè)正常鋼軌的振動(dòng)信號(hào),作為無(wú)波磨故障的振動(dòng)信號(hào);通過(guò)1D?DCGAN 生成1 000 個(gè)偽樣本波磨數(shù)據(jù),作為波磨振動(dòng)信號(hào)。為了對(duì)不同程度的波磨情況進(jìn)行模擬,分別將波磨信號(hào)的幅值A(chǔ)縮小至0.5A和增大至1.5A,即訓(xùn)練數(shù)據(jù)包含4 000 個(gè)無(wú)波磨數(shù)據(jù)和1 000 個(gè)0.5A的波磨數(shù)據(jù)、1 000 個(gè)1.0A的波磨數(shù)據(jù)和1 000個(gè)1.5A的波磨數(shù)據(jù)。將測(cè)試數(shù)據(jù)中的200 條鋼軌波磨數(shù)據(jù)進(jìn)行幅值放大、縮小后,與200 條無(wú)波磨數(shù)據(jù)一起作為測(cè)試數(shù)據(jù)集。

    3.4 試驗(yàn)結(jié)果及分析

    計(jì)算不同幅值波磨數(shù)據(jù)的15種時(shí)域統(tǒng)計(jì)指標(biāo),使用PCA 對(duì)特征進(jìn)行降維后,取前3 個(gè)主成分分量作為其特征,見圖4??芍翰バ盘?hào)的幅值越高,波磨振動(dòng)的信號(hào)越大,使用一維深度卷積生成對(duì)抗網(wǎng)絡(luò)的鋼軌波磨偽樣本與真實(shí)波磨樣本特征越接近。

    圖4 不同幅值波磨數(shù)據(jù)特征對(duì)比

    為了表明研究數(shù)據(jù)量對(duì)模型訓(xùn)練以及分類準(zhǔn)確率的影響,分別用不同數(shù)量的數(shù)據(jù)進(jìn)行模型訓(xùn)練,并對(duì)比其準(zhǔn)確率。使用KNN、SVM 和DT 三種分類模型進(jìn)行模型訓(xùn)練和分類,訓(xùn)練模型的數(shù)據(jù)量分別為樣本總量的20%、40%、60%、80%和100%,其分類準(zhǔn)確率與樣本數(shù)量的關(guān)系見圖5??芍喝N分類模型的平均分類準(zhǔn)確率都隨著訓(xùn)練數(shù)據(jù)樣本數(shù)量的增加而增大,其中平均分類準(zhǔn)確率最高的是SVM,可以達(dá)到96.344%;準(zhǔn)確率隨數(shù)據(jù)規(guī)模變化最明顯的是DT,從89.5%增至92.5%。

    圖5 三種分類模型的平均分類準(zhǔn)確率隨樣本數(shù)量的變化曲線

    4 結(jié)論

    本文提出一種基于1D?DCGAN 的鋼軌波磨識(shí)別方法。首先使用1D?DCGAN 生成的鋼軌波磨偽樣本進(jìn)行數(shù)據(jù)集擴(kuò)充;然后計(jì)算數(shù)據(jù)的15 種時(shí)域統(tǒng)計(jì)指標(biāo),進(jìn)行PCA降維,并取前3個(gè)主成分分量作為鋼軌波磨的特征進(jìn)行分類模型的訓(xùn)練;最后通過(guò)訓(xùn)練好的模型進(jìn)行鋼軌波磨識(shí)別。主要結(jié)論如下:

    1)1D?DCGAN 生成的偽樣本,在時(shí)域、頻域和時(shí)域統(tǒng)計(jì)指標(biāo)方面都與真實(shí)樣本基本一致,可以在樣本不足時(shí)用來(lái)擴(kuò)充數(shù)據(jù)集。

    2)KNN、SVM 和DT 三種分類模型的平均分類準(zhǔn)確率都隨著訓(xùn)練樣本規(guī)模增加而增大,其中SVM 的平均分類準(zhǔn)確率最高。

    猜你喜歡
    波磨鋼軌幅值
    城軌小半徑曲線鋼軌波磨的成因與整治研究綜述
    北京地鐵鋼軌波磨現(xiàn)狀及防治措施
    鋼軌焊后雙中頻感應(yīng)加熱工藝研究
    基于S變換的交流電網(wǎng)幅值檢測(cè)系統(tǒng)計(jì)算機(jī)仿真研究
    電子制作(2017年7期)2017-06-05 09:36:13
    正序電壓幅值檢測(cè)及諧波抑制的改進(jìn)
    低壓電力線信道脈沖噪聲的幅值與寬度特征
    高速鐵路鋼軌疲勞過(guò)程的超聲非線性系數(shù)表征
    國(guó)內(nèi)外高速鐵路鋼軌性能對(duì)比研究
    高速鐵路鋼軌波磨對(duì)車輛—軌道動(dòng)態(tài)響應(yīng)的影響
    基于零序電壓幅值增量的消弧線圈調(diào)諧新方法
    亚洲av免费高清在线观看| 中文字幕制服av| 色吧在线观看| 深夜a级毛片| 99久久久亚洲精品蜜臀av| 免费看光身美女| 18+在线观看网站| 99久久人妻综合| 国产视频内射| 国产精品美女特级片免费视频播放器| 国模一区二区三区四区视频| 岛国在线免费视频观看| 综合色丁香网| 国产一区二区亚洲精品在线观看| 久久久久久久久大av| 搡女人真爽免费视频火全软件| 欧美日韩一区二区视频在线观看视频在线 | 日韩 亚洲 欧美在线| 久久6这里有精品| 国产黄片视频在线免费观看| 乱人视频在线观看| 深夜精品福利| 婷婷色综合大香蕉| 中文欧美无线码| 中文字幕av成人在线电影| 美女脱内裤让男人舔精品视频 | 国产av在哪里看| 亚洲精品日韩在线中文字幕 | 禁无遮挡网站| 六月丁香七月| 国产精品久久视频播放| 一区二区三区四区激情视频 | 国产精品人妻久久久久久| 波多野结衣巨乳人妻| 亚洲丝袜综合中文字幕| 亚洲天堂国产精品一区在线| 国产精品久久电影中文字幕| 男人狂女人下面高潮的视频| 国产精华一区二区三区| 好男人视频免费观看在线| 69人妻影院| 欧美性猛交黑人性爽| 波多野结衣高清作品| 免费不卡的大黄色大毛片视频在线观看 | 九色成人免费人妻av| 久久婷婷人人爽人人干人人爱| 亚洲av中文av极速乱| 国产一级毛片七仙女欲春2| 国产精品乱码一区二三区的特点| 又爽又黄a免费视频| 久久国内精品自在自线图片| 日本一二三区视频观看| 欧美xxxx黑人xx丫x性爽| 亚洲aⅴ乱码一区二区在线播放| 黄色视频,在线免费观看| 色哟哟·www| 国产一区亚洲一区在线观看| 国产一区二区在线观看日韩| 国产一区二区三区av在线 | 性插视频无遮挡在线免费观看| 国产高清有码在线观看视频| 久久精品综合一区二区三区| 成人特级av手机在线观看| 国产成人福利小说| 一边亲一边摸免费视频| 国产精品人妻久久久影院| 免费不卡的大黄色大毛片视频在线观看 | 男插女下体视频免费在线播放| 国产精品久久久久久av不卡| 亚洲av一区综合| 免费看美女性在线毛片视频| 久久久精品94久久精品| 色哟哟·www| 九色成人免费人妻av| 91久久精品国产一区二区三区| 丰满乱子伦码专区| 丰满乱子伦码专区| 深夜精品福利| 亚洲内射少妇av| 波多野结衣高清作品| 亚洲熟妇中文字幕五十中出| 五月玫瑰六月丁香| 长腿黑丝高跟| 国国产精品蜜臀av免费| 舔av片在线| 听说在线观看完整版免费高清| 伦精品一区二区三区| 亚洲国产精品久久男人天堂| 日本在线视频免费播放| 精品久久久噜噜| 可以在线观看毛片的网站| 尾随美女入室| 久久久久国产网址| 夜夜看夜夜爽夜夜摸| 国产精品一区二区性色av| 国产亚洲5aaaaa淫片| 搡女人真爽免费视频火全软件| 久久婷婷人人爽人人干人人爱| 日韩欧美三级三区| 国产男人的电影天堂91| 国产国拍精品亚洲av在线观看| 亚洲国产欧美人成| 欧美一区二区亚洲| 日本av手机在线免费观看| 91久久精品国产一区二区成人| 国产精品久久视频播放| 日本欧美国产在线视频| av女优亚洲男人天堂| 国产69精品久久久久777片| 男人狂女人下面高潮的视频| 成年av动漫网址| 十八禁国产超污无遮挡网站| 午夜视频国产福利| 亚洲欧美日韩东京热| 成人亚洲欧美一区二区av| 12—13女人毛片做爰片一| 日日啪夜夜撸| 哪个播放器可以免费观看大片| 男女下面进入的视频免费午夜| 热99在线观看视频| 亚洲熟妇中文字幕五十中出| 一级二级三级毛片免费看| 在线观看美女被高潮喷水网站| 插阴视频在线观看视频| 久久精品夜色国产| 老司机福利观看| 国产精品麻豆人妻色哟哟久久 | a级毛片免费高清观看在线播放| 天堂√8在线中文| 综合色av麻豆| 丝袜美腿在线中文| 国产精品久久久久久精品电影小说 | 91精品一卡2卡3卡4卡| 伦理电影大哥的女人| av在线蜜桃| 久久久午夜欧美精品| 精品人妻视频免费看| 精品一区二区免费观看| 久久精品影院6| av在线播放精品| 插阴视频在线观看视频| 99热全是精品| 国产精品蜜桃在线观看 | 黄色欧美视频在线观看| 亚洲无线观看免费| 亚洲av二区三区四区| 国产三级中文精品| av免费在线看不卡| 麻豆国产97在线/欧美| а√天堂www在线а√下载| 亚洲婷婷狠狠爱综合网| 欧美激情在线99| 99久国产av精品国产电影| 丰满乱子伦码专区| 午夜精品一区二区三区免费看| 日韩欧美一区二区三区在线观看| 99热全是精品| 国产一级毛片七仙女欲春2| 插阴视频在线观看视频| 亚洲婷婷狠狠爱综合网| 国产淫片久久久久久久久| 欧美丝袜亚洲另类| 桃色一区二区三区在线观看| 麻豆av噜噜一区二区三区| 亚洲第一电影网av| 99精品在免费线老司机午夜| 日韩,欧美,国产一区二区三区 | 亚洲在久久综合| 欧美区成人在线视频| 午夜福利高清视频| 看非洲黑人一级黄片| 婷婷六月久久综合丁香| 久久久久久伊人网av| 精品一区二区免费观看| 国产一区二区激情短视频| 精品免费久久久久久久清纯| 久久国内精品自在自线图片| 乱系列少妇在线播放| 亚洲精品456在线播放app| 国产av一区在线观看免费| 午夜视频国产福利| 国产精品一区二区三区四区免费观看| 寂寞人妻少妇视频99o| 国产乱人偷精品视频| 91精品国产九色| 欧美高清成人免费视频www| 九九在线视频观看精品| 亚洲真实伦在线观看| 美女黄网站色视频| 精品免费久久久久久久清纯| 国产精品一区www在线观看| 久久久久久久久久久丰满| 国产精品野战在线观看| 亚洲最大成人中文| 成人毛片60女人毛片免费| 在线观看一区二区三区| 一级毛片电影观看 | 五月伊人婷婷丁香| 精品午夜福利在线看| 欧美性猛交╳xxx乱大交人| 国产亚洲精品久久久com| 国产v大片淫在线免费观看| 观看免费一级毛片| 日韩国内少妇激情av| 成人毛片a级毛片在线播放| 国产精品无大码| 一边亲一边摸免费视频| 亚洲精品国产av成人精品| 午夜激情欧美在线| 国产久久久一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 国产一级毛片七仙女欲春2| av在线观看视频网站免费| 国产一区亚洲一区在线观看| 一个人免费在线观看电影| 亚洲欧美中文字幕日韩二区| 色视频www国产| 日本在线视频免费播放| 国产精品不卡视频一区二区| 99久久久亚洲精品蜜臀av| 久久精品久久久久久噜噜老黄 | 国产一区二区在线观看日韩| 舔av片在线| 毛片女人毛片| 国内久久婷婷六月综合欲色啪| 国产亚洲精品av在线| 亚洲中文字幕一区二区三区有码在线看| 色吧在线观看| videossex国产| 精品久久久久久成人av| 天天一区二区日本电影三级| 神马国产精品三级电影在线观看| 精品午夜福利在线看| 欧美zozozo另类| 免费看av在线观看网站| 午夜视频国产福利| 12—13女人毛片做爰片一| 亚洲国产欧洲综合997久久,| 最近手机中文字幕大全| 国产av一区在线观看免费| 白带黄色成豆腐渣| 午夜a级毛片| 午夜亚洲福利在线播放| 在线免费观看不下载黄p国产| 亚洲在线观看片| 女同久久另类99精品国产91| 久久这里只有精品中国| 日日摸夜夜添夜夜爱| 精品无人区乱码1区二区| 大又大粗又爽又黄少妇毛片口| 成人特级黄色片久久久久久久| 亚洲欧美成人综合另类久久久 | 久久久久久久久大av| 在线观看美女被高潮喷水网站| 美女高潮的动态| 亚洲无线观看免费| 此物有八面人人有两片| 狂野欧美激情性xxxx在线观看| 男人的好看免费观看在线视频| 亚洲不卡免费看| 国产成人午夜福利电影在线观看| 国产精品久久久久久久电影| 精华霜和精华液先用哪个| 此物有八面人人有两片| av免费观看日本| 草草在线视频免费看| 有码 亚洲区| 免费大片18禁| 嫩草影院新地址| 一区二区三区免费毛片| 欧美色视频一区免费| 午夜激情福利司机影院| 美女被艹到高潮喷水动态| 亚洲国产日韩欧美精品在线观看| 久久久久网色| 18禁裸乳无遮挡免费网站照片| 久久久色成人| 好男人在线观看高清免费视频| 毛片女人毛片| 日本撒尿小便嘘嘘汇集6| 一级黄片播放器| 亚洲国产精品久久男人天堂| 亚洲人成网站高清观看| 国产在视频线在精品| 高清在线视频一区二区三区 | 中文字幕制服av| 色尼玛亚洲综合影院| 小蜜桃在线观看免费完整版高清| 成人特级黄色片久久久久久久| 亚洲成人av在线免费| 亚洲精品乱码久久久久久按摩| 日本色播在线视频| 免费不卡的大黄色大毛片视频在线观看 | 免费人成在线观看视频色| 黄色配什么色好看| 看十八女毛片水多多多| 寂寞人妻少妇视频99o| 99国产极品粉嫩在线观看| 亚洲欧洲国产日韩| 爱豆传媒免费全集在线观看| 久久久精品94久久精品| 日日摸夜夜添夜夜爱| 亚洲精品国产av成人精品| 久久草成人影院| 国内精品一区二区在线观看| 国内少妇人妻偷人精品xxx网站| 国产伦精品一区二区三区视频9| 亚洲精品自拍成人| 成人毛片60女人毛片免费| 热99在线观看视频| 久久久精品94久久精品| 一边摸一边抽搐一进一小说| 久久精品夜色国产| 特大巨黑吊av在线直播| 精品久久久久久久久久免费视频| 丰满乱子伦码专区| 最新中文字幕久久久久| 亚洲精品粉嫩美女一区| 六月丁香七月| 国产美女午夜福利| 干丝袜人妻中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 国内精品美女久久久久久| 久久婷婷人人爽人人干人人爱| 九九热线精品视视频播放| 亚洲av免费高清在线观看| 日韩一区二区视频免费看| 黄色一级大片看看| 久久精品人妻少妇| 久久精品久久久久久噜噜老黄 | 中国美白少妇内射xxxbb| 成人午夜精彩视频在线观看| 精品久久久噜噜| 亚洲欧美精品综合久久99| 成人午夜精彩视频在线观看| 日本黄色视频三级网站网址| 99国产极品粉嫩在线观看| 美女黄网站色视频| 观看美女的网站| 欧美色欧美亚洲另类二区| 久久久久久久亚洲中文字幕| 少妇高潮的动态图| 1000部很黄的大片| www日本黄色视频网| 又粗又爽又猛毛片免费看| 美女国产视频在线观看| 亚洲内射少妇av| 欧美成人a在线观看| 精品一区二区三区视频在线| 国产视频内射| 99热这里只有是精品50| 天天一区二区日本电影三级| 99视频精品全部免费 在线| 最好的美女福利视频网| 色综合亚洲欧美另类图片| 亚洲最大成人av| 国内精品宾馆在线| 亚洲图色成人| 熟妇人妻久久中文字幕3abv| 永久网站在线| 蜜臀久久99精品久久宅男| 男人舔女人下体高潮全视频| 国产视频首页在线观看| 能在线免费看毛片的网站| 性欧美人与动物交配| 日韩一区二区视频免费看| 99视频精品全部免费 在线| 欧美日韩在线观看h| 成年女人永久免费观看视频| 日本一二三区视频观看| av黄色大香蕉| www.av在线官网国产| 欧美一区二区精品小视频在线| 精华霜和精华液先用哪个| 欧美精品国产亚洲| 69人妻影院| 熟女人妻精品中文字幕| 欧美另类亚洲清纯唯美| 免费观看精品视频网站| 高清午夜精品一区二区三区 | 国产成人精品婷婷| 免费搜索国产男女视频| 成人二区视频| 国产精品女同一区二区软件| 亚洲av一区综合| 亚洲人成网站在线播| 国产伦一二天堂av在线观看| 波多野结衣高清无吗| 日韩成人av中文字幕在线观看| 国产美女午夜福利| 少妇高潮的动态图| 美女大奶头视频| 一级毛片久久久久久久久女| av女优亚洲男人天堂| 日韩国内少妇激情av| 久久99蜜桃精品久久| 国产成人91sexporn| 亚洲欧美日韩卡通动漫| 久久人人爽人人片av| 国产精品一区二区三区四区免费观看| 国产 一区精品| 久久久久久大精品| 人妻少妇偷人精品九色| 亚洲最大成人中文| 久久国产乱子免费精品| 99九九线精品视频在线观看视频| 欧美xxxx性猛交bbbb| 伊人久久精品亚洲午夜| 欧美精品一区二区大全| 91麻豆精品激情在线观看国产| 精华霜和精华液先用哪个| 欧美性感艳星| 精品免费久久久久久久清纯| 男人的好看免费观看在线视频| 亚洲在线自拍视频| 欧美最新免费一区二区三区| 欧美日本亚洲视频在线播放| 国产黄a三级三级三级人| 日本欧美国产在线视频| 精品久久久久久成人av| 不卡视频在线观看欧美| 欧美日韩在线观看h| 有码 亚洲区| 男人的好看免费观看在线视频| 久久精品国产亚洲网站| 日韩一区二区视频免费看| 国产精品1区2区在线观看.| 国产一区二区在线观看日韩| 免费观看人在逋| 麻豆精品久久久久久蜜桃| 久久久欧美国产精品| 91久久精品电影网| 国产高清有码在线观看视频| 婷婷色av中文字幕| 伦精品一区二区三区| 亚洲成人中文字幕在线播放| 91在线精品国自产拍蜜月| a级毛色黄片| 久久久久久久久久成人| 日本与韩国留学比较| 亚洲欧洲日产国产| 日韩欧美精品免费久久| 精品久久久久久成人av| 亚洲色图av天堂| 亚洲精品成人久久久久久| 日韩成人av中文字幕在线观看| a级毛色黄片| 小蜜桃在线观看免费完整版高清| 男的添女的下面高潮视频| 又粗又爽又猛毛片免费看| 又粗又硬又长又爽又黄的视频 | 国产一区二区三区在线臀色熟女| 亚洲国产色片| 欧美潮喷喷水| 日本免费一区二区三区高清不卡| 国产亚洲91精品色在线| 欧美色视频一区免费| 美女被艹到高潮喷水动态| 悠悠久久av| 狂野欧美激情性xxxx在线观看| 日本撒尿小便嘘嘘汇集6| 久久久欧美国产精品| 亚洲乱码一区二区免费版| 一区福利在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产成人aa在线观看| 国产成人精品一,二区 | avwww免费| 亚洲av熟女| 啦啦啦韩国在线观看视频| 此物有八面人人有两片| 青青草视频在线视频观看| 国产精品久久电影中文字幕| 国产av在哪里看| 天美传媒精品一区二区| 亚洲在线自拍视频| 日韩人妻高清精品专区| 日本免费a在线| 国产精品美女特级片免费视频播放器| 午夜视频国产福利| 国产精品免费一区二区三区在线| 国产老妇伦熟女老妇高清| 午夜免费男女啪啪视频观看| 色吧在线观看| 亚洲欧美中文字幕日韩二区| 91狼人影院| 久久久成人免费电影| 99在线视频只有这里精品首页| 午夜视频国产福利| 麻豆久久精品国产亚洲av| 亚洲欧美清纯卡通| 精品久久久久久久人妻蜜臀av| 亚洲国产精品sss在线观看| 校园人妻丝袜中文字幕| 中文字幕免费在线视频6| av福利片在线观看| 国产黄色小视频在线观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品乱码久久久久久按摩| 国产老妇伦熟女老妇高清| 中文字幕久久专区| 在线免费观看的www视频| 亚洲第一电影网av| 3wmmmm亚洲av在线观看| a级毛色黄片| 日韩人妻高清精品专区| 禁无遮挡网站| 日韩成人伦理影院| 久久人人爽人人爽人人片va| 午夜激情欧美在线| 国产免费一级a男人的天堂| 99在线人妻在线中文字幕| 美女 人体艺术 gogo| 日本-黄色视频高清免费观看| 色播亚洲综合网| 又爽又黄a免费视频| 国产三级中文精品| 男人和女人高潮做爰伦理| 18+在线观看网站| 中文欧美无线码| 国产黄片视频在线免费观看| 在线国产一区二区在线| 最新中文字幕久久久久| 国产精品久久久久久精品电影小说 | 99riav亚洲国产免费| 欧美丝袜亚洲另类| 尾随美女入室| 欧美bdsm另类| 亚洲精品国产av成人精品| 小蜜桃在线观看免费完整版高清| 日韩欧美精品v在线| 国产精品久久视频播放| 特大巨黑吊av在线直播| h日本视频在线播放| 哪里可以看免费的av片| 欧美丝袜亚洲另类| 亚洲欧美成人精品一区二区| 精华霜和精华液先用哪个| av视频在线观看入口| av在线老鸭窝| 91午夜精品亚洲一区二区三区| 日日干狠狠操夜夜爽| 日韩欧美国产在线观看| 欧美变态另类bdsm刘玥| 岛国毛片在线播放| 长腿黑丝高跟| 最近的中文字幕免费完整| 国产老妇女一区| 97在线视频观看| 亚洲精品456在线播放app| 亚洲成人av在线免费| 国产一区二区亚洲精品在线观看| 如何舔出高潮| 夜夜看夜夜爽夜夜摸| 国产av不卡久久| 非洲黑人性xxxx精品又粗又长| 日本熟妇午夜| 国产精品人妻久久久影院| 51国产日韩欧美| 欧美+日韩+精品| 免费av毛片视频| 啦啦啦啦在线视频资源| 中文亚洲av片在线观看爽| 中文字幕熟女人妻在线| 成年av动漫网址| 我的老师免费观看完整版| 一个人免费在线观看电影| 中文在线观看免费www的网站| 亚洲成人久久性| 日韩成人av中文字幕在线观看| 国产 一区 欧美 日韩| 免费看a级黄色片| 性色avwww在线观看| 国产免费男女视频| 99在线人妻在线中文字幕| 久久午夜亚洲精品久久| 久久亚洲国产成人精品v| 99热网站在线观看| 青春草国产在线视频 | 男插女下体视频免费在线播放| 国产高清三级在线| 91麻豆精品激情在线观看国产| 一夜夜www| 女人被狂操c到高潮| 成年女人看的毛片在线观看| 国产精品女同一区二区软件| 亚洲精品乱码久久久v下载方式| 日本与韩国留学比较| 日本爱情动作片www.在线观看| 国产成人a区在线观看| 久久久久久久久久成人| 国产精品国产三级国产av玫瑰| 亚洲av第一区精品v没综合| 国产久久久一区二区三区| av天堂在线播放| 国产精品嫩草影院av在线观看| 99热这里只有是精品50| av国产免费在线观看| 精品一区二区三区视频在线| 国内精品一区二区在线观看| 亚洲av成人av| 久久精品国产鲁丝片午夜精品| 久久久久久九九精品二区国产| av卡一久久| 欧美激情国产日韩精品一区| 亚洲欧美成人综合另类久久久 | 日韩人妻高清精品专区| 精品免费久久久久久久清纯| 亚洲av成人精品一区久久| 3wmmmm亚洲av在线观看| 简卡轻食公司| 免费无遮挡裸体视频| 日韩高清综合在线| 国产精品99久久久久久久久| 国产高清视频在线观看网站|