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

    改進LSTM 的脫軌系數(shù)預(yù)測方法

    2022-03-07 06:57:50劉怡伶張文靜劉學(xué)文
    軟件導(dǎo)刊 2022年2期
    關(guān)鍵詞:短時記憶集上輪軌

    張 卜,劉怡伶,張文靜,劉學(xué)文

    (上海工程技術(shù)大學(xué)機械與汽車工程學(xué)院,上海 201620)

    0 引言

    隨著社會經(jīng)濟的不斷發(fā)展和鐵路技術(shù)的不斷提升,輪軌列車的運行速度至今已達到每小時500 多公里,運載重量也在逐漸增加,但隨之而來的列車運行安全成為最棘手的問題之一,其中最為危險的工況就是脫軌。

    在軌道交通建設(shè)領(lǐng)域,輪軌之間的作用研究自始至終是車輛—基礎(chǔ)耦合動力學(xué)的核心所在。理清輪軌在不同工況、不同運行狀態(tài)等條件下的接觸狀態(tài),是研究輪軌作用的焦點問題。隨著深度學(xué)習(xí)理論的迅速發(fā)展,除直接以理論方法研究脫軌機理外,也可以通過深度學(xué)習(xí)方法挖掘大數(shù)據(jù)中的潛在規(guī)律,Sitton 等建立多個BP 網(wǎng)絡(luò)分別加以訓(xùn)練,網(wǎng)絡(luò)收斂后將預(yù)測結(jié)果與驗證集進行比較,根據(jù)網(wǎng)絡(luò)性能表現(xiàn)對它們賦予權(quán)重,將多個帶權(quán)重的BP 網(wǎng)絡(luò)組成投票系統(tǒng)以預(yù)測貨沖擊對軌、橋帶來的破壞;Samira 等和Serdar 等在神經(jīng)網(wǎng)絡(luò)基礎(chǔ)上,建立一種新的預(yù)測模型,用于列車脫軌事故風(fēng)險評估;Dindar 等設(shè)計一種由貝葉斯分類器(NB)組成的網(wǎng)絡(luò),該網(wǎng)絡(luò)上層NB 的分類結(jié)果會成為下一層NB 的輸入,可對因極端天氣造成的脫軌事故進行準(zhǔn)確分類;Zhang 等利用深層信念網(wǎng)絡(luò)(Deep Belief Network,DBN)對列車運行安全性進行預(yù)測,并且考慮乘客舒適性;張俊甲等利用附加權(quán)重法對BP 神經(jīng)網(wǎng)絡(luò)進行改進,克服了傳統(tǒng)BP 網(wǎng)絡(luò)學(xué)習(xí)率調(diào)節(jié)困難的缺點,得到一種能夠有效預(yù)測脫軌系數(shù)的網(wǎng)絡(luò)模型;陳皓等建立基于BR網(wǎng)絡(luò)的NARX 網(wǎng)絡(luò)預(yù)測模型以預(yù)測列車脫軌系數(shù),解決了傳統(tǒng)BP 網(wǎng)絡(luò)預(yù)測脫軌系數(shù)成本高、精度低的問題;楊桐等采用T-snake 模型分析輪軌接觸圖像,推算出輪軌相對橫向位移量,結(jié)合遺傳算法優(yōu)化小波神經(jīng)網(wǎng)絡(luò),對相對位移、速度、加速度以及輪重減載率和脫軌系數(shù)間的映射進行建模。

    列車運行時所產(chǎn)生的各項數(shù)據(jù)按照被采樣的先后順序記錄下來,是典型的時間序列。時間序列的特殊性在于其中的各樣例不是孤立的,之前的樣例會影響后面樣例的狀態(tài),長短時記憶神經(jīng)網(wǎng)絡(luò)(LSTM)所具有的循環(huán)結(jié)構(gòu)使它能夠掌握樣例間的關(guān)聯(lián)性。陳渝等采用LSTM 對醫(yī)院門診量數(shù)據(jù)進行建模,取得了較好的預(yù)測精度;李彬等借助隨機梯度下降法對LSTM 進行優(yōu)化,對風(fēng)力發(fā)電功率概率進行預(yù)測,得出一種能夠預(yù)測未來200h 風(fēng)電功率概率函數(shù)的神經(jīng)網(wǎng)絡(luò);Ugurlu 等以均值絕對誤差(MAE)和均方根誤差(RMSE)作為性能度量,使用LSTM 精確預(yù)測了電價的價格滯后現(xiàn)象。

    要想保證深度學(xué)習(xí)效果,需要準(zhǔn)備大量的訓(xùn)練數(shù)據(jù),然而列車實車試驗成本高昂,鐵路日常運營中收集到的數(shù)據(jù),一般研究機構(gòu)也難以獲得,通過仿真方法獲取數(shù)據(jù)集是解決該難題的方法之一。近年來,仿真軟件的性能逐步提升,其中SIMPACK 在多體動力學(xué),尤其是輪軌接觸方面的仿真效果非常好。陳楊在多體系統(tǒng)動力學(xué)軟件SIMPACK 中建立了整車動力學(xué)仿真模型,并對輪軌接觸面和懸掛系統(tǒng)做了處理,加入了軌道不平順激勵模型,得到了較為真實的仿真場景,利用仿真研究了各向不平順對蛇形失穩(wěn)臨界速度的影響;王海濤在SIMPACK 中建立了車—線—隧剛?cè)狁詈蟿恿W(xué)模型,利用仿真與實測數(shù)據(jù)研究了車輪抬升力,車輪踏面磨耗數(shù)與隧道壁振幅等參數(shù)對列車運行的影響;何銀川在SIMPACK 中建立了馱背運輸車輛系統(tǒng)振動模型,根據(jù)仿真分析了臨界速度、曲線通過能力和懸架對垂向穩(wěn)定性的影響,并且做了實際試驗,與仿真結(jié)果是相吻合的;姚常偉通過FLUENT、ANSYS 和SIMPACK 等仿真軟件對受側(cè)風(fēng)影響的列車模型進行了仿真,仿真結(jié)果與通過理論計算得到的風(fēng)載荷對車輛的影響是等效的。

    上述科研成果都能證明SIMPACK 多體動力學(xué)仿真軟件在研究輪軌接觸方面的可靠性,因此借助仿真工具得到數(shù)據(jù)集,再利用深度學(xué)習(xí)方法分析數(shù)據(jù),從而得到關(guān)鍵參數(shù)與列車輪軌狀態(tài)之間的關(guān)系是可行的。

    本研究采用SIMPACK 仿真軟件收集列車通過彎道時的各項參數(shù),包括最大脫軌系數(shù)、最大橫向振動加速度、橫向重心偏移量、車速、彎道半徑和軌道超高,組成訓(xùn)練數(shù)據(jù)集。

    循環(huán)神經(jīng)網(wǎng)絡(luò),特別是其中的長短時記憶神經(jīng)網(wǎng)絡(luò)已被諸多研究證實適用于處理時間序列數(shù)據(jù),因此本研究也將以長短時記憶神經(jīng)網(wǎng)絡(luò)作為基礎(chǔ)模型。但是一方面,長短時記憶網(wǎng)絡(luò)訓(xùn)練過程對學(xué)習(xí)率非常敏感,學(xué)習(xí)率的調(diào)整往往要根據(jù)學(xué)者的經(jīng)驗或反復(fù)嘗試才得以確定,這樣調(diào)參不僅要耗費大量時間,也未必能達到最優(yōu)效果;另一方面,得益于長短時記憶網(wǎng)絡(luò)強大的學(xué)習(xí)能力,網(wǎng)絡(luò)可以高效地從訓(xùn)練集中找到屬性與標(biāo)記之間的映射,但是訓(xùn)練后的網(wǎng)絡(luò)在測試集上的表現(xiàn)往往不如訓(xùn)練集,也即發(fā)生了過擬合現(xiàn)象。針對這兩種情況,本研究提出使用簡便、有效的動態(tài)學(xué)習(xí)率和Dropout 方法針對學(xué)習(xí)率和網(wǎng)絡(luò)結(jié)構(gòu)進行優(yōu)化,將經(jīng)過改進的長短時記憶神經(jīng)網(wǎng)絡(luò)用于列車通過彎道時的脫軌系數(shù)預(yù)測。實驗結(jié)果表明,改進的長短時記憶神經(jīng)網(wǎng)絡(luò)能夠較好地預(yù)測脫軌系數(shù)隨工況變化趨勢,為進一步的脫軌系數(shù)預(yù)測研究提供了支持。

    1 長短時記憶神經(jīng)網(wǎng)絡(luò)(LSTM)

    1.1 LSTM 神經(jīng)網(wǎng)絡(luò)

    LSTM 神經(jīng)網(wǎng)絡(luò)改良于循環(huán)神經(jīng)網(wǎng)絡(luò)(RNN),除了與RNN 相同的外部循環(huán)外,還有本身節(jié)點間的自循環(huán)。自循環(huán)由遺忘門、輸入門和輸出門控制,這三者可根據(jù)當(dāng)前時刻的輸入從隱層中刪除、添加、提取信息,改善了RNN 長期依賴導(dǎo)致梯度消失和梯度爆炸的問題,圖1 展示了LSTM 的單元結(jié)構(gòu)。

    Fig.1 Unit structure of LSTM圖1 LSTM 的單元結(jié)構(gòu)

    輸入門根據(jù)

    x

    、輸入門的輸入偏置

    b

    h

    決定計算當(dāng)前時刻輸出的狀態(tài)

    c

    (,可用式(2)計算。其中,

    U

    為單元體輸入門的輸入權(quán)重,

    W

    為單元體輸入門的循環(huán)權(quán)重,Sigmoid 為激活函數(shù),則:

    1.2 動態(tài)學(xué)習(xí)率方法

    神經(jīng)網(wǎng)絡(luò)的權(quán)重更新離不開梯度下降法,傳統(tǒng)的梯度下降方法如式(5)所示。其中,

    W

    為第

    i

    次更新時的權(quán)重,

    η

    為學(xué)習(xí)率,

    L

    為損失函數(shù)的值。

    學(xué)習(xí)率

    η

    的取值對網(wǎng)絡(luò)訓(xùn)練效果影響很大,取值過大會導(dǎo)致更新跳過最優(yōu)權(quán)重,或是損失函數(shù)震蕩,網(wǎng)絡(luò)難以收斂,取值過小會導(dǎo)致網(wǎng)絡(luò)收斂慢、訓(xùn)練效率差。為解決該問題可以引入動態(tài)學(xué)習(xí)率方法,每輪訓(xùn)練都根據(jù)當(dāng)前損失函數(shù)的梯度Δ

    L

    重新分配下一輪訓(xùn)練的學(xué)習(xí)率,其計算方法如式(6)所示。

    動態(tài)學(xué)習(xí)率可按式(7)推導(dǎo),可以看出,隨著訓(xùn)練輪數(shù)的增加,動態(tài)學(xué)習(xí)率會越來越小,防止越過最優(yōu)解,但當(dāng)Δ

    L

    ≤0 時,即網(wǎng)絡(luò)尚未出現(xiàn)收斂趨勢時,動態(tài)學(xué)習(xí)率的降低速度會變慢,以保證網(wǎng)絡(luò)訓(xùn)練效率。若學(xué)習(xí)率過低,會引發(fā)梯度消失問題,因此

    η

    不會小于0.01 。

    1.3 Dropout 方法

    Fig.2 Principle of dropout圖2 Dropout 原理

    2 數(shù)據(jù)集準(zhǔn)備

    2.1 SIMPACK 仿真場景

    該仿真場景中包含模型與線路兩個部分,其中模型中有17個剛體、1個車身、2個構(gòu)架、4 組輪對、8個軸箱和2個搖枕。其中,車身、構(gòu)架和輪對考慮到了前后伸縮、左右晃動、上下浮動以及點頭、搖頭和翻滾6個自由度,軸箱考慮點頭,只有1個自由度,搖枕與車身剛性連接,不考慮其自由度,故該模型共有50個自由度。該模型包含54個力元,因只考慮單節(jié)車輛,故沒有車鉤力,其中軸套與構(gòu)架間使用43 號力元,一系懸掛中的彈簧使用85 號力元,一二系懸掛中的阻尼以及蛇形減振器選用6 號力元,二系懸掛的兩個剪切彈簧使用79 號力元,車架與搖枕之間用13 號力元作為緩沖器以防止二者碰撞,輪軌接觸力元采用78 號力元并考慮到了它的非線性特性,該模型如圖1 所示。

    Fig.3 SIMPACK simulation model圖3 SIMPACK 仿真模型

    車輛行駛線路包括:直線—過渡曲線—彎道—過渡曲線—直線,列車在直線的行駛狀態(tài)不是研究重點,故設(shè)置為10m,過渡曲線用線性過渡,長度為10m,彎道的長度為500m,半徑可調(diào)。引入的軌道水平不平順激勵參考美國6 級譜,式(8)是其譜密度分析式,各參數(shù)如表1 所示。

    Table 1 Parameter values of the six-order spectrum of the US orbit表1 美國軌道六級譜的參數(shù)值

    2.2 數(shù)據(jù)集及其預(yù)處理

    列車在彎道上的脫軌系數(shù),會受到線路、車速等眾多因素的影響,本文將重點放在車速、重心橫向偏移量、半徑和超高上,研究這三者對脫軌系數(shù)的影響。通過不同條件下的仿真實驗,提取車速、半徑、超高、重心橫向偏移量作為屬性,對應(yīng)的脫軌系數(shù)作為標(biāo)記,共得到如表2 所示的1 269 條樣例。其中,超高按照式(9)計算,一般應(yīng)該按5mm的整倍數(shù)設(shè)置,且按照國家標(biāo)準(zhǔn)一般不能超過120~150mm,但為了在仿真中能得到盡可能大的脫軌系數(shù),暫不考慮這些限制。

    Table 2 Derailment coefficient data set表2 脫軌系數(shù)數(shù)據(jù)集

    3 脫軌系數(shù)預(yù)測

    3.1 評價標(biāo)準(zhǔn)

    為了評估模型最終在測試集上的表現(xiàn),采用均方根誤差(RMSE)作為評價標(biāo)準(zhǔn),如式(10)所示,其為各條樣例預(yù)測值和真實值偏差的平方之和的平方根,反映了預(yù)測值偏離真實值的程度,當(dāng)預(yù)測值出現(xiàn)較大誤差時,會導(dǎo)致RMSE 激增,故RMSE 能夠較為全面地反映預(yù)測的精確程度。

    3.2 實驗結(jié)果與分析

    將數(shù)據(jù)集的前890條作為訓(xùn)練集,其余作為測試集。初始學(xué)習(xí)率為0.01,循環(huán)次數(shù)500次,batch_size 為20,keep_prob 參數(shù)決定了dropout 層保留輸出權(quán)重的概率,設(shè)置為0.9、0.8、0.7,即單元有10%、20%、30%的概率失活,分別進行訓(xùn)練,選取測試結(jié)果較好的一種作為最終參數(shù)。3 種使用不同keep_prob 參數(shù)的LSTM 在訓(xùn)練集上都有不錯的表現(xiàn),圖4 展示了500 輪訓(xùn)練中,損失函數(shù)值變化曲線。

    Fig.4 Loss function in training圖4 訓(xùn)練中的損失函數(shù)

    幾種不同Keep_Prob 參數(shù)下的LSTM 在測試集上的RMSE 如表3 所示。傳統(tǒng)的LSTM 雖然在訓(xùn)練集上表現(xiàn)優(yōu)異,但在測試集上的誤差反而較大,出現(xiàn)了過擬合現(xiàn)象。引入Dropout 方法,設(shè)置keep_prob 為0.9 和0.8 時,過擬合現(xiàn)象有所改善,但進一步減少保留輸出權(quán)重的概率時,測試集上的誤差再次增加,出現(xiàn)了欠擬合趨勢。因此,最后選取Keep_Prob 參數(shù)為0.8。

    Table 3 Prediction errors of different keep_prob parameter networks表3 不同keep_prob 參數(shù)網(wǎng)絡(luò)的預(yù)測誤差

    最終測試結(jié)果如圖5 所示,模型經(jīng)過優(yōu)化后有著較好的擬合精度,預(yù)測曲線雖然未和真實值完全重合,但是模型準(zhǔn)確預(yù)測出了脫軌系數(shù)隨工況的變化趨勢。其中,在113 號樣例附近,模型準(zhǔn)確地預(yù)測出了因為車速突然下降而導(dǎo)致的脫軌系數(shù)斷崖式下跌,在其他位置模型也把握住了車速越高、彎道半徑越小,則脫軌系數(shù)越大的特點,這一特點在測試集和訓(xùn)練集上都相同,說明模型具有了較好的泛化性能??傮w而言,雖然預(yù)測值和真實值之間有一定誤差,但均方根誤差控制在了0.01 以內(nèi),與未經(jīng)優(yōu)化的LSTM預(yù)測結(jié)果相比,改善較大,RMSE 減小近24%,對于預(yù)測脫軌系數(shù)的變化趨勢有一定參考價值。

    Fig.5 Prediction results of 2 kinds of LSTM圖5 兩種LSTM 預(yù)測結(jié)果

    4 結(jié)語

    列車脫軌系數(shù)是評價列車運行安全的重要參考,影響著軌道運輸?shù)脑O(shè)計、實驗、運維等方方面面。本研究針對在彎道上受重心偏移量影響的列車脫軌系數(shù),使用深度學(xué)習(xí)方法,建立長短時記憶神經(jīng)網(wǎng)絡(luò)進行脫軌系數(shù)預(yù)測。為了解決模型學(xué)習(xí)率調(diào)整困難的問題,引入了動態(tài)學(xué)習(xí)率方法,為了解決過擬合問題,使用了Dropout 方法。驗證集表明,雖然在預(yù)測精度上仍有可以改善的地方,但通過優(yōu)化,模型在驗證集上的RMSE 減小了23.9%,不僅較好地預(yù)測出脫軌系數(shù)變化趨勢,也把握住了脫軌系數(shù)驟然變化情況??傮w而言,這種優(yōu)化的長短時記憶網(wǎng)絡(luò)時間序列預(yù)測方法能應(yīng)用在脫軌系數(shù)預(yù)測上,能夠為鐵路安全運營提供數(shù)據(jù)支持。下一步研究將考慮更多因素對脫軌系數(shù)的影響,比如車輪抬升量、車鉤力及橫風(fēng)等,并嘗試采用其他優(yōu)化方法以建立更加準(zhǔn)確的預(yù)測方法,進一步挖掘列車運行數(shù)據(jù)的實用價值。

    猜你喜歡
    短時記憶集上輪軌
    基于長短時記憶神經(jīng)網(wǎng)絡(luò)的動力電池剩余容量預(yù)測方法
    Cookie-Cutter集上的Gibbs測度
    中低速磁浮道岔與輪軌道岔的差異
    鏈完備偏序集上廣義向量均衡問題解映射的保序性
    復(fù)扇形指標(biāo)集上的分布混沌
    短時記憶、長時記憶對英語聽力的影響
    中低速磁浮與輪軌交通信號系統(tǒng)的差異
    短時記憶理論的影響
    非線性穩(wěn)態(tài)曲線通過時輪軌滾動接觸的數(shù)值求解方法
    不同輪軌沖角下高速輪軌穩(wěn)態(tài)滾動接觸的蠕滑特性
    最后的刺客免费高清国语| 亚洲欧美日韩卡通动漫| 日本a在线网址| 国产精品嫩草影院av在线观看| 老司机午夜福利在线观看视频| 欧美潮喷喷水| 波野结衣二区三区在线| 熟女电影av网| 国内揄拍国产精品人妻在线| 97超视频在线观看视频| 欧美xxxx黑人xx丫x性爽| 熟女人妻精品中文字幕| 精品久久久久久久久久久久久| 91在线观看av| 国产精品久久久久久亚洲av鲁大| 亚洲av成人av| 久久久精品欧美日韩精品| 嫩草影院精品99| 国产亚洲精品久久久com| 麻豆精品久久久久久蜜桃| 国产高清三级在线| 欧美最新免费一区二区三区| 热99在线观看视频| 亚洲五月天丁香| 久久草成人影院| 国产激情偷乱视频一区二区| 日韩人妻高清精品专区| 国产国拍精品亚洲av在线观看| 欧美成人一区二区免费高清观看| 成年av动漫网址| 亚洲,欧美,日韩| 91午夜精品亚洲一区二区三区| 亚洲av.av天堂| 亚洲国产精品成人综合色| 成年免费大片在线观看| 伦精品一区二区三区| 日韩人妻高清精品专区| 亚洲无线在线观看| 在线看三级毛片| 成人精品一区二区免费| 日韩一本色道免费dvd| 波多野结衣高清作品| 成人无遮挡网站| 成人亚洲精品av一区二区| 国产午夜精品论理片| 国产三级中文精品| 久久精品综合一区二区三区| 97超级碰碰碰精品色视频在线观看| 成人欧美大片| 大型黄色视频在线免费观看| 久久人人爽人人爽人人片va| 在现免费观看毛片| 亚洲在线观看片| 国产亚洲精品久久久久久毛片| 精品国产三级普通话版| 在线观看午夜福利视频| 又爽又黄a免费视频| 看十八女毛片水多多多| 亚洲美女视频黄频| 日日摸夜夜添夜夜爱| 亚洲精品456在线播放app| 久久久久免费精品人妻一区二区| 一级a爱片免费观看的视频| 超碰av人人做人人爽久久| 99热这里只有是精品在线观看| 麻豆一二三区av精品| 麻豆久久精品国产亚洲av| 久久久久久久久久黄片| av.在线天堂| 菩萨蛮人人尽说江南好唐韦庄 | 麻豆乱淫一区二区| 少妇的逼水好多| 免费在线观看影片大全网站| 久久人人爽人人爽人人片va| 国产片特级美女逼逼视频| 精品午夜福利视频在线观看一区| 少妇熟女aⅴ在线视频| 国产成人精品久久久久久| 日韩欧美一区二区三区在线观看| 一边摸一边抽搐一进一小说| 亚洲在线自拍视频| 身体一侧抽搐| 91久久精品国产一区二区三区| 久久精品夜色国产| 久久久久免费精品人妻一区二区| 一a级毛片在线观看| 国产高清有码在线观看视频| 乱系列少妇在线播放| 亚洲三级黄色毛片| 亚洲精品一区av在线观看| 99九九线精品视频在线观看视频| 精品99又大又爽又粗少妇毛片| 欧美日韩国产亚洲二区| 婷婷色综合大香蕉| 欧美成人a在线观看| 久久精品国产鲁丝片午夜精品| 一级黄色大片毛片| 亚洲综合色惰| 国产高清视频在线观看网站| 免费电影在线观看免费观看| 小说图片视频综合网站| 亚洲一区二区三区色噜噜| 久久久午夜欧美精品| 长腿黑丝高跟| 色综合站精品国产| 亚洲va在线va天堂va国产| 国产在视频线在精品| 精品少妇黑人巨大在线播放 | 日本免费a在线| 国产麻豆成人av免费视频| 国产aⅴ精品一区二区三区波| 久久精品国产亚洲av涩爱 | 91久久精品电影网| 亚洲av.av天堂| 99久久无色码亚洲精品果冻| 中出人妻视频一区二区| 欧美高清成人免费视频www| 国产爱豆传媒在线观看| 99久久中文字幕三级久久日本| 性色avwww在线观看| 内地一区二区视频在线| 免费无遮挡裸体视频| 一级a爱片免费观看的视频| a级毛片a级免费在线| 春色校园在线视频观看| 成年av动漫网址| 搡老熟女国产l中国老女人| 久久欧美精品欧美久久欧美| 国产一区二区亚洲精品在线观看| 日日摸夜夜添夜夜爱| 久久久久性生活片| 在线看三级毛片| 国产综合懂色| 日韩,欧美,国产一区二区三区 | 欧美成人精品欧美一级黄| 国产色婷婷99| 成年版毛片免费区| 男人舔奶头视频| 色吧在线观看| 欧美激情在线99| 99热全是精品| 亚洲av成人精品一区久久| 亚洲av五月六月丁香网| 人人妻,人人澡人人爽秒播| av免费在线看不卡| aaaaa片日本免费| 两性午夜刺激爽爽歪歪视频在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 少妇的逼好多水| 三级毛片av免费| 久久亚洲国产成人精品v| 男人舔女人下体高潮全视频| 亚洲欧美日韩高清专用| 久久精品国产亚洲av涩爱 | 午夜免费激情av| 午夜久久久久精精品| 人妻夜夜爽99麻豆av| 最近视频中文字幕2019在线8| 99久国产av精品| 波野结衣二区三区在线| 噜噜噜噜噜久久久久久91| 91精品国产九色| 少妇高潮的动态图| 亚洲精品粉嫩美女一区| 中文字幕免费在线视频6| 熟女电影av网| 日本免费a在线| 亚洲不卡免费看| 亚洲中文字幕一区二区三区有码在线看| 欧美一区二区国产精品久久精品| 精品久久久久久久久久久久久| 91久久精品电影网| 美女内射精品一级片tv| 97超碰精品成人国产| 久久精品国产自在天天线| 美女cb高潮喷水在线观看| 真人做人爱边吃奶动态| 国产精品久久久久久精品电影| 别揉我奶头~嗯~啊~动态视频| 91狼人影院| 日韩欧美一区二区三区在线观看| 欧美丝袜亚洲另类| 亚洲欧美精品自产自拍| 岛国在线免费视频观看| 午夜视频国产福利| 高清日韩中文字幕在线| 极品教师在线视频| 日本黄大片高清| 免费看美女性在线毛片视频| 啦啦啦观看免费观看视频高清| 久久人人精品亚洲av| 一级毛片aaaaaa免费看小| 国产免费一级a男人的天堂| 国产精品精品国产色婷婷| 色尼玛亚洲综合影院| 精品人妻偷拍中文字幕| 99久久无色码亚洲精品果冻| av在线亚洲专区| 精品一区二区免费观看| 蜜桃亚洲精品一区二区三区| 亚洲第一区二区三区不卡| 国产精品精品国产色婷婷| 尤物成人国产欧美一区二区三区| 精品人妻偷拍中文字幕| 中国美女看黄片| 少妇的逼水好多| 久久99热这里只有精品18| 日韩 亚洲 欧美在线| 亚洲av五月六月丁香网| 亚洲国产欧美人成| 国产高清三级在线| 国产精品av视频在线免费观看| 国产成人aa在线观看| 国产精品一区二区三区四区免费观看 | 国产蜜桃级精品一区二区三区| 亚洲人与动物交配视频| 永久网站在线| 亚洲五月天丁香| 人人妻人人看人人澡| 草草在线视频免费看| 国产美女午夜福利| 欧美三级亚洲精品| av卡一久久| 国产亚洲精品综合一区在线观看| 五月玫瑰六月丁香| 久久精品国产亚洲av天美| 日韩精品青青久久久久久| 在线观看午夜福利视频| 欧美激情久久久久久爽电影| 舔av片在线| 中文字幕人妻熟人妻熟丝袜美| 波多野结衣高清无吗| 亚洲va在线va天堂va国产| 欧美+日韩+精品| 亚洲自偷自拍三级| 亚洲人成网站在线播| 亚洲激情五月婷婷啪啪| 给我免费播放毛片高清在线观看| 日韩国内少妇激情av| 深夜精品福利| 日本在线视频免费播放| 3wmmmm亚洲av在线观看| 欧美xxxx性猛交bbbb| av福利片在线观看| 一个人免费在线观看电影| 亚洲三级黄色毛片| 亚州av有码| 成熟少妇高潮喷水视频| 亚洲美女黄片视频| 在线免费观看不下载黄p国产| 国产亚洲精品久久久久久毛片| 一个人免费在线观看电影| 成人一区二区视频在线观看| www.色视频.com| 女人被狂操c到高潮| 麻豆久久精品国产亚洲av| 欧美成人a在线观看| 中国美白少妇内射xxxbb| 天堂动漫精品| 精品人妻一区二区三区麻豆 | 久久午夜亚洲精品久久| 精品不卡国产一区二区三区| 卡戴珊不雅视频在线播放| 少妇的逼水好多| 97超视频在线观看视频| 日韩欧美国产在线观看| 成人国产麻豆网| a级毛色黄片| 三级国产精品欧美在线观看| 男女视频在线观看网站免费| 国产成人一区二区在线| 国产伦精品一区二区三区视频9| 一本久久中文字幕| 国产aⅴ精品一区二区三区波| 中文亚洲av片在线观看爽| 99久久中文字幕三级久久日本| 春色校园在线视频观看| 亚洲性久久影院| 国产精品,欧美在线| 最近2019中文字幕mv第一页| 又黄又爽又免费观看的视频| 天堂√8在线中文| 成人亚洲精品av一区二区| 插逼视频在线观看| 久久久精品欧美日韩精品| 国产精品久久电影中文字幕| 国产一级毛片七仙女欲春2| 免费大片18禁| 亚洲在线自拍视频| 99国产极品粉嫩在线观看| 亚洲va在线va天堂va国产| 亚洲中文日韩欧美视频| 男人舔女人下体高潮全视频| 亚洲国产色片| 欧美激情国产日韩精品一区| 亚洲av美国av| 男人狂女人下面高潮的视频| 亚洲成人久久性| 欧美zozozo另类| 亚洲欧美日韩高清在线视频| 久久久午夜欧美精品| 亚洲aⅴ乱码一区二区在线播放| 亚洲在线观看片| www日本黄色视频网| 国产亚洲精品久久久com| 日韩在线高清观看一区二区三区| 日本 av在线| 亚洲国产精品成人久久小说 | 亚洲无线观看免费| 日韩精品有码人妻一区| 久久久色成人| 欧美绝顶高潮抽搐喷水| 丰满人妻一区二区三区视频av| 51国产日韩欧美| 人妻制服诱惑在线中文字幕| 精品午夜福利在线看| 99热这里只有是精品在线观看| 小说图片视频综合网站| 日韩强制内射视频| 精品久久久久久久久久久久久| 蜜臀久久99精品久久宅男| 嫩草影院入口| 久久精品国产亚洲av天美| 久久久久久久亚洲中文字幕| 在线观看免费视频日本深夜| 俄罗斯特黄特色一大片| 亚洲久久久久久中文字幕| 自拍偷自拍亚洲精品老妇| 毛片女人毛片| 久久久久国产网址| 欧美性猛交黑人性爽| 国产又黄又爽又无遮挡在线| 亚洲真实伦在线观看| 国产淫片久久久久久久久| 白带黄色成豆腐渣| 有码 亚洲区| 国产精品永久免费网站| 麻豆久久精品国产亚洲av| 国产高清有码在线观看视频| 精品久久久久久久人妻蜜臀av| 色综合亚洲欧美另类图片| 美女高潮的动态| 欧美一区二区精品小视频在线| 一区二区三区免费毛片| 尾随美女入室| 精品午夜福利在线看| 别揉我奶头~嗯~啊~动态视频| 一级黄色大片毛片| 尤物成人国产欧美一区二区三区| 成人亚洲欧美一区二区av| 淫秽高清视频在线观看| 国内久久婷婷六月综合欲色啪| 成人亚洲欧美一区二区av| 男女之事视频高清在线观看| 天堂av国产一区二区熟女人妻| 国产私拍福利视频在线观看| 如何舔出高潮| 国产三级中文精品| 插阴视频在线观看视频| 国产免费男女视频| 亚洲精品色激情综合| 成年av动漫网址| 91午夜精品亚洲一区二区三区| 亚洲精品影视一区二区三区av| 青春草视频在线免费观看| 美女大奶头视频| 高清毛片免费观看视频网站| 天天一区二区日本电影三级| 一级毛片电影观看 | 国产三级中文精品| 亚洲四区av| 在线国产一区二区在线| 久久精品国产清高在天天线| 亚洲一区高清亚洲精品| 亚洲中文字幕一区二区三区有码在线看| 欧美zozozo另类| 国产精品不卡视频一区二区| 精品午夜福利在线看| 成人亚洲精品av一区二区| 成人毛片a级毛片在线播放| 久久天躁狠狠躁夜夜2o2o| 看黄色毛片网站| 国产私拍福利视频在线观看| 中出人妻视频一区二区| 成人高潮视频无遮挡免费网站| www.色视频.com| 色视频www国产| 麻豆久久精品国产亚洲av| 噜噜噜噜噜久久久久久91| 午夜精品一区二区三区免费看| 一区福利在线观看| 非洲黑人性xxxx精品又粗又长| 亚洲专区国产一区二区| 我要搜黄色片| 国产探花在线观看一区二区| 亚洲欧美成人综合另类久久久 | 蜜臀久久99精品久久宅男| 午夜福利视频1000在线观看| 亚洲欧美日韩东京热| 99在线人妻在线中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 日韩精品青青久久久久久| 综合色丁香网| 欧美潮喷喷水| 久久久久久伊人网av| 又爽又黄无遮挡网站| 成人av一区二区三区在线看| 卡戴珊不雅视频在线播放| 精品久久久久久久人妻蜜臀av| 国产私拍福利视频在线观看| 久久精品国产亚洲网站| 国产老妇女一区| 观看美女的网站| 麻豆一二三区av精品| 国产精品国产三级国产av玫瑰| 精品午夜福利视频在线观看一区| 男女啪啪激烈高潮av片| 最后的刺客免费高清国语| 亚洲欧美日韩高清在线视频| avwww免费| 色视频www国产| 午夜福利高清视频| 色播亚洲综合网| 尾随美女入室| 精品久久国产蜜桃| 久久久a久久爽久久v久久| 无遮挡黄片免费观看| 国产片特级美女逼逼视频| 大香蕉久久网| 成人漫画全彩无遮挡| 啦啦啦韩国在线观看视频| 久久99热这里只有精品18| 97超级碰碰碰精品色视频在线观看| 不卡一级毛片| 高清毛片免费观看视频网站| 亚洲欧美日韩东京热| 狂野欧美激情性xxxx在线观看| 亚洲在线观看片| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲精品一区av在线观看| 国产片特级美女逼逼视频| 亚洲av二区三区四区| 男女之事视频高清在线观看| 亚洲中文日韩欧美视频| 亚洲av免费在线观看| 欧美国产日韩亚洲一区| 欧美日韩在线观看h| 可以在线观看的亚洲视频| 午夜福利在线观看免费完整高清在 | 久久精品91蜜桃| 婷婷六月久久综合丁香| 97热精品久久久久久| 亚洲欧美日韩东京热| 嫩草影院精品99| 日日摸夜夜添夜夜添小说| 亚洲自偷自拍三级| 国产av一区在线观看免费| 看非洲黑人一级黄片| 美女高潮的动态| av在线蜜桃| 国产精华一区二区三区| 亚洲色图av天堂| 搡老妇女老女人老熟妇| 亚洲av中文字字幕乱码综合| 亚洲18禁久久av| 精品久久久久久久久久久久久| 国产精品无大码| aaaaa片日本免费| 亚洲国产高清在线一区二区三| 亚洲高清免费不卡视频| 大香蕉久久网| 精品免费久久久久久久清纯| 久久久久免费精品人妻一区二区| 欧美一区二区国产精品久久精品| 午夜精品一区二区三区免费看| 我的女老师完整版在线观看| 成人二区视频| 国产av在哪里看| 久久精品夜色国产| 久久精品人妻少妇| 成人av在线播放网站| 不卡一级毛片| 日本一二三区视频观看| 国产精品1区2区在线观看.| 亚洲av熟女| 真人做人爱边吃奶动态| 亚洲内射少妇av| 日韩大尺度精品在线看网址| 国产成年人精品一区二区| 国产色爽女视频免费观看| 精品一区二区免费观看| 高清毛片免费看| 男人舔女人下体高潮全视频| 精品日产1卡2卡| 深夜a级毛片| 亚洲一区高清亚洲精品| 亚洲乱码一区二区免费版| 国产视频内射| 欧美一区二区亚洲| 欧美高清成人免费视频www| 丰满乱子伦码专区| 国产伦精品一区二区三区四那| 亚洲自偷自拍三级| 亚洲人成网站高清观看| 九九久久精品国产亚洲av麻豆| 一卡2卡三卡四卡精品乱码亚洲| 国产午夜精品久久久久久一区二区三区 | 性色avwww在线观看| 日本成人三级电影网站| 男女视频在线观看网站免费| a级一级毛片免费在线观看| 国产国拍精品亚洲av在线观看| 成人美女网站在线观看视频| 国产高清视频在线观看网站| 成人美女网站在线观看视频| 在线免费观看的www视频| 老师上课跳d突然被开到最大视频| 精品一区二区免费观看| av黄色大香蕉| 香蕉av资源在线| 精品久久久久久久久久久久久| 日韩欧美一区二区三区在线观看| 国内久久婷婷六月综合欲色啪| 久久久久久久久大av| 天堂影院成人在线观看| 日韩欧美一区二区三区在线观看| 精品99又大又爽又粗少妇毛片| 一区福利在线观看| 亚洲图色成人| 菩萨蛮人人尽说江南好唐韦庄 | 在线天堂最新版资源| 午夜免费激情av| 亚洲第一区二区三区不卡| 国产精品嫩草影院av在线观看| 精品一区二区三区人妻视频| 最近在线观看免费完整版| 12—13女人毛片做爰片一| 国产美女午夜福利| 久久精品人妻少妇| 久久精品夜色国产| av中文乱码字幕在线| 欧美中文日本在线观看视频| 国产真实伦视频高清在线观看| 干丝袜人妻中文字幕| 精品人妻视频免费看| 极品教师在线视频| 成人国产麻豆网| 最新在线观看一区二区三区| 午夜影院日韩av| 成人二区视频| 五月伊人婷婷丁香| 别揉我奶头 嗯啊视频| 91麻豆精品激情在线观看国产| 国产又黄又爽又无遮挡在线| 亚洲av美国av| 成人一区二区视频在线观看| 九九热线精品视视频播放| 国产精品不卡视频一区二区| 我要搜黄色片| 国产毛片a区久久久久| 久久午夜福利片| 你懂的网址亚洲精品在线观看 | 精品人妻视频免费看| 99热这里只有精品一区| 国产三级中文精品| 最新在线观看一区二区三区| 精品久久久久久久久亚洲| 久久精品国产鲁丝片午夜精品| 亚洲av免费在线观看| 欧美潮喷喷水| 国产精品爽爽va在线观看网站| 欧美激情国产日韩精品一区| 99国产精品一区二区蜜桃av| 网址你懂的国产日韩在线| 亚洲欧美清纯卡通| 成人av一区二区三区在线看| 久久精品国产亚洲网站| 搡老岳熟女国产| 亚洲美女搞黄在线观看 | 中文在线观看免费www的网站| 在线免费观看的www视频| 99精品在免费线老司机午夜| 舔av片在线| 国产白丝娇喘喷水9色精品| 超碰av人人做人人爽久久| 国产亚洲精品av在线| 最近视频中文字幕2019在线8| 国产亚洲精品综合一区在线观看| 日本爱情动作片www.在线观看 | 久久天躁狠狠躁夜夜2o2o| 国产一级毛片七仙女欲春2| 欧美xxxx性猛交bbbb| 国产三级在线视频| 人妻久久中文字幕网| 亚洲精品色激情综合| 日日摸夜夜添夜夜添av毛片| 亚洲高清免费不卡视频| 日本免费一区二区三区高清不卡| 中出人妻视频一区二区| 日日撸夜夜添| 国产在线男女| 亚洲av熟女| 身体一侧抽搐| 性欧美人与动物交配| 高清日韩中文字幕在线| 午夜免费男女啪啪视频观看 | 亚洲一区高清亚洲精品| 精品久久久久久久久亚洲| 成人二区视频| 亚洲av美国av| 欧美+日韩+精品| 99久久精品国产国产毛片| 国产三级在线视频| 欧美国产日韩亚洲一区| 久久人人爽人人爽人人片va|