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

    基于振動信號的柴油機NOx排放虛擬傳感研究

    2022-07-22 14:31:34胡桂誠李國興和超亮楊甜甜
    關鍵詞:缸蓋缸內瞬態(tài)

    胡桂誠,李國興,沈 亮,和超亮,楊甜甜

    (1.太原理工大學 機械與運載工程學院, 太原 030024; 2.中國重汽集團濟南特種車有限公司,濟南 250000)

    0 引言

    基于NEDC循環(huán)等穩(wěn)態(tài)工況開展的汽車尾氣排放性能檢測無法對內燃機瞬態(tài)排放特性進行有效評估,因此實際行駛污染物排放(real drive emission,RDE)檢測方法得到了重視與推廣。輕型柴油車實際行駛中的氮氧化物(NOx)排放約為車輛認證時的3~4倍[1]。RDE檢測對車輛瞬態(tài)排放監(jiān)測技術提出了更高要求,然而目前尚未有兼顧成本和瞬態(tài)測量性能的測量技術得到推廣。PEMS車載瞬態(tài)測試設備昂貴且笨重,無法滿足推廣需要。商用NOx傳感器則無法滿足較高的瞬態(tài)測量要求,不利于閉環(huán)控制和排放改善的實現(xiàn)[2]。而虛擬傳感技術為NOx瞬態(tài)監(jiān)測提供了更多可能,獲得了越來越多的關注與研究。

    已有的虛擬NOx傳感器多是基于實測缸內壓力信號作為輸入變量開發(fā)的。這是因為NOx的產(chǎn)生主要與缸內燃燒行為有關,通常以缸內壓力作為缸內燃燒行為的評價指標[3]。

    2011年,Guardiola等[4]開發(fā)了一個柴油發(fā)動機快速NOx預測模型,選擇缸內壓力和電子控制單元中平均變量作為預測模型的輸入信號,實現(xiàn)了NOx的實時預測。2014年,F(xiàn)ormentin等[5]提出了一種基于發(fā)動機轉速和缸壓的NOx預測方法。通過主成分分析構建預測變量,并且使用L2正則化技術建立了NOx預測模型。Christian等[6]基于靜態(tài)多項式黑盒建模為重型越野柴油機設計了虛擬NOx傳感器,使用較少的回歸量實現(xiàn)了高精度的NOx預測。2019年,Paul等[7]提出一種基于現(xiàn)象學模型的NOx虛擬傳感器,與物理NOx傳感器相比,其頻率響應更好。Dipankar等[8]提出的基于神經(jīng)網(wǎng)絡的虛擬NOx傳感器進一步證實了將基于缸壓的虛擬NOx傳感器應用于發(fā)動機排放實時控制的可行性。然而,缸內壓力傳感器比商用NOx傳感器昂貴數(shù)10倍,使得基于壓力的虛擬NOx傳感器依然不便于推廣應用。有必要探索更具成本效益的虛擬NOx傳感器。

    缸內燃燒行為與NOx的產(chǎn)生密切相關,同時又是發(fā)動機產(chǎn)生振動和噪聲的最重要激勵源之一。通過分析結構振動可以評估缸內燃燒情況,進而有望對燃燒產(chǎn)生的污染物排放水平加以預測。唐娟等[9]通過有限元模擬得到在僅施加燃燒壓力時的缸蓋振動響應,研究表明,缸蓋表面位移可以有效地反映缸內壓力變化。程勇等[10]研究表明,缸蓋振動速度和缸壓升高率具有相似的趨勢。Zhao等[11]采用已測加速度信號來識別時域中的缸內燃燒事件。Polonowski等[12]分析了柴油機的加速度信號和燃燒放熱引起的壓力信號之間的相干性。昝晶等[13]使用結構相似性(SSIM)分析方法從不同的數(shù)字圖像中識別和提取信息。基于SSIM算法可以獲得2個三維矩陣之間更多的相關信息,這意味著可以直接從振動時頻圖中提取與燃燒事件相關的局部振動事件。楊甜甜[14]分析了結構振動與缸內壓力的相似性,研究了基于振動信號的缸內壓力重構技術,為基于振動的排放監(jiān)測奠定了基礎。

    本文中提出一種基于振動信號重構缸壓,通過構建預測模型來預測NOx排放的方法。首先,對缸蓋振動時頻圖譜和缸壓二階導時頻圖譜進行相似性分析,獲得相似性評價結果。將振動信號作為輸入?yún)?shù),以相似性評價結果作為控制參數(shù),通過形態(tài)學重構計算出缸壓二階導信號。對重構的缸壓信號進行分段提取預測變量,基于穩(wěn)態(tài)工況下提取的預測變量與排放數(shù)據(jù)構造主成分回歸(PCR)預測模型,對比分析排放預測效果。最后,通過瞬態(tài)工況預測的排放與實測對比驗證了所提預測模型的可行性。結果表明:基于振動信號重構的缸壓信號所建立的回歸模型能夠很好地預測柴油機各種運行工況下的NOx排放。預測值和測量值之間的確定系數(shù)R2為0.971~0.995,誤差小于10%。瞬態(tài)試驗表明,預測NOx值能夠較好地跟隨實測NOx排放變化,虛擬NOx傳感器響應速度更快。

    1 試驗設計

    為了便于提取振動特征信號,選用一臺單缸柴油機開展相關試驗,具體參數(shù)見表1。

    表1 發(fā)動機技術參數(shù)

    試驗設備如圖1所示,由立式單缸柴油機、測功機、測控系統(tǒng)、數(shù)據(jù)采集儀、缸內壓力傳感器、振動傳感器組成。缸套與缸蓋上分別安裝有振動加速度傳感器,用以測量柴油機結構振動。NOx排放用AVL FTIR多成分測試儀進行測量。此外,同步采集缸內壓力、發(fā)動機轉速和曲軸位置信息用于數(shù)據(jù)結果分析。

    圖1 試驗臺架示意圖

    排放試驗分為2個部分:建模試驗與驗證試驗。建模試驗是由15個工況組成的穩(wěn)態(tài)試驗,包括1 000~1 800 r/min范圍內以200 r/min為間隔的5個轉速工況,每種轉速施加3種負荷:10、30和50 N·m。驗證試驗為定負荷變轉速瞬態(tài)試驗。在恒定的發(fā)動機負荷下(10、30、50 N·m),發(fā)動機轉速階梯上升(以200 r/min為間隔從1 000 r/min增加到1 800 r/min),隨后階梯下降(以200 r/min為間隔從1 800 r/min降低到1 000 r/min)。排放數(shù)據(jù)采樣頻率為1 Hz,其他數(shù)據(jù)采樣頻率為96 kHz,每個運行工況的采樣時間為30 s。為避免單次測量引起誤差,每個工況下采集2次數(shù)據(jù)。

    2 基于振動信號的缸壓重構

    基于振動信號構建虛擬NOx傳感器的關鍵是從測量的振動信號中準確提取缸內壓力相關信息。已有關于缸內壓力信息提取的研究通常是根據(jù)振動信號中預先選定的峰值或拐點來識別缸壓曲線的變化趨勢,能夠有效評估部分燃燒事件的相位信息,但無法精確評估缸內壓力的振幅或能量信息。本文中提出的基于振動信號的燃燒壓力時頻重構的方法實現(xiàn)了缸壓信號的相位與振幅的同步重構,其步驟主要分為以下4步。

    2.1 數(shù)據(jù)選擇

    內燃機的燃燒室由缸蓋的底部,缸套的內壁和活塞頂面組成。燃氣沖擊和燃氣驅動下的活塞側擊是內燃機振動響應主要的2種激勵源。相較于缸蓋,同時受到2種激勵的缸套往往呈現(xiàn)出更為復雜的振動響應特征。

    實測振動信號包含隨機噪聲并呈現(xiàn)顯著非平穩(wěn)性,因此對原始振動數(shù)據(jù)進行了均值濾波和連續(xù)小波變換(continuous wavelet transform,CWT)。圖2是在1 600 r/min-30 N·m運行工況下缸套振動、缸蓋振動和缸壓二階導的CWT時頻圖。

    相較于缸套振動,缸蓋振動時頻圖的能量分布規(guī)律更接近于缸壓二階導時頻圖。因此,后續(xù)選取缸蓋振動響應作為缸壓重構的數(shù)據(jù)源。

    2.2 歸一化

    對缸蓋振動和缸壓二階導時頻圖進行相似性分析之前,首先需要對信號的時頻矩陣進行歸一化處理,將圖2中不同物理量的時頻圖轉化為亮度峰值一致的RGB圖像。結構相似性(structural similarity, SSIM)分析的適用輸入為灰度圖像,為此需要對RGB圖像進行灰度化處理,去除色相和飽和度信息。由圖3可看出,灰度轉化后圖像顯示的特征信息少于時頻圖譜,實質上轉換過程并不會造成信息缺失,而是由顯示模式造成的。為了優(yōu)化微弱特征的顯示和相似性對比效果,需要使用伽瑪校正方法進行灰度圖像的修正處理。所謂伽瑪校正就是對圖像的伽瑪曲線進行編輯,從而實現(xiàn)對圖像進行非線性色調編輯的方法,檢測出圖像信號中的深色部分和淺色部分,并使兩者比例增大,從而提高圖像對比度效果。如圖4所示,通過伽瑪校正后,更多局域特征得以在灰度圖中凸顯。

    圖2 1 600 r/min-30 N·m工況下缸套和缸蓋振動

    圖3 時頻圖譜的灰度顯示

    圖4 伽瑪校正后的灰度圖

    信號與壓力二階導的CWT時頻

    2.3 結構相似性分析

    SSIM(structural similarity)是一種用來衡量圖像相似度指標的算法,基本原理由亮度對比、對比度對比、結構對比3部分組成,SSIM指數(shù)如式(1)所示。

    (1)

    對伽馬校正前后的缸蓋振動和缸壓二階導灰度圖分別進行結構相似性分析,獲得SSIM指數(shù)均值。從圖5中可以看出,經(jīng)過伽馬校正后的缸壓二階導和缸蓋振動的SSIM均值(0.359 8)低于伽馬校正前的SSIM均值(0.657 1)。這主要是由于經(jīng)過伽馬校正后,非平穩(wěn)響應特征得以凸顯,與燃燒事件無關的局域相似性虛高得到了有效抑制,導致整個工作循環(huán)的SSIM均值下降。事實上,在燃燒上止點(TDC)附近,經(jīng)過伽馬校正后的SSIM圖譜顯示出更加明顯的與燃燒有關的相似結構,如三角區(qū)域所示。

    圖5 伽瑪校正前、后SSIM圖譜

    2.4 形態(tài)學重構

    基于振動信號重構缸壓的實質是依據(jù)相似性矩陣對振動信號中與燃燒事件相關的局域響應特征進行濾波提取。將缸蓋振動灰度圖作為輸入?yún)?shù),以SSIM分析計算得出的相似矩陣作為控制參數(shù),可通過形態(tài)學重構計算出壓力二階導的灰度圖。如圖6所示,重構的壓力二階導灰度圖與實測灰度圖之間有著較好的相似性。

    圖6 壓力二階導的灰度圖

    從圖7(a)與(b)可看出重構的壓力二階導準確地保留了與燃燒事件相關的波動信息。由于重構過程導致的計算誤差,重構的壓力二階導時頻曲面比原始壓力二階導曲面更加光滑。然后通過逆小波變換將重構的壓力二階導時頻圖轉化為曲線圖。如圖8所示,重構的壓力曲線與原始曲線的相位吻合度較好,并且幅值也基本接近。尤其是燃燒上止點附近的燃燒波動特征得到了較好地重構復現(xiàn),如圖8中局域放大圖所示。

    圖7 壓力二階導的時頻顯示

    圖8 重構的燃燒壓力二階導曲線

    3 基于缸壓重構的NOx排放虛擬傳感

    本文中提出一種NOx排放的主成分回歸(PCR)預測方法,通過主成分分析(PCA)獲取預測特征值與目標NOx值之間對應關系并構建了NOx回歸預測模型。PCR是在預測變量較多且預測變量高度相關甚至共線的情況下,對響應變量進行預測建模的一種方法。通過主成分分析,由原始預測變量線性組合出新的預測變量稱為主成分,以各主成分作為新自變量代替原來預測變量與目標值進行回歸分析,構建回歸預測模型。

    首先,需要從重構的缸壓曲線中選取合適的預測變量用于模型構建。不同于依賴壓力曲線輪廓信息拾取的預測變量構造方法,本文中以缸壓曲線分段的均方根值、峰值和峭度作為預測變量。NOx排放的產(chǎn)生與缸內燃燒行為關系緊密,為縮減輸入維度,只選取燃燒過程相關的缸壓片段來構造預測變量。選擇進氣門關閉角(IVC)和排氣門打開角(EVO)之間變化劇烈的缸壓曲線作為預測變量的選取區(qū)間,根據(jù)不同燃燒階段以噴油開始角(ISA)、上止點(TDC)和燃燒結束角(EOC)將其分割為4段,對應的曲軸轉角分別為:218~338° CA,338~360° CA,360~445° CA和445~485° CA。分別求取4個區(qū)間中壓力曲線的均方根值(RMS)、峰值和峭度,完成12個預測變量的選取。

    3.1 基于穩(wěn)態(tài)試驗的模型構建

    首先,對12個預測變量和排放結果進行主成分分析,然后分別采用2個主成分和5個主成分建立用于NOx排放的PCR預測模型,以預測15種不同工況下的排放水平。如圖9所示,由2個主成分和5個主成分分別構建的PCR模型預測值與測量值有良好的一致性。

    圖9 NOx排放測量值和預測值

    由5個主成分構成的PCR模型的回歸分析相關性系數(shù)R2=0.955,略高于由2個主成分構成的PCR模型的回歸分析相關性系數(shù)R2=0.971??梢钥闯?,從重構的缸壓曲線中選擇的預測變量與排放值之間存在密切的相關性。

    如圖10所示,主成分分析方法表明前2個主成分解釋了99.53%的變量,前5個主成分解釋了99.98%的變量,可以看出,僅由2個主成分構成的回歸模型可以很好地預測NOx排放。

    圖10 主成分的累計方差貢獻率

    如圖11所示,預測和實測的NOx之間的相對誤差在多數(shù)工況下誤差均低于10%,其中前2個主成分平均相對誤差4.73%,最大相對誤差11.56%,前5個主成分平均相對誤差3.39%,最大相對誤差10.25%。在實際應用中,基于2個主成分建立回歸模型能夠準確高效地預測NOx排放水平。

    圖11 NOx排放實測值和預測值的相對誤差

    3.2 基于瞬態(tài)試驗的模型驗證

    物理傳感器的響應時間受到排放物采樣傳輸速率和傳感器分析耗時的綜合影響,響應時間往往超過1 s,難以滿足發(fā)動機瞬態(tài)排放控制需要。臺架試驗中NOx測量設備的最小采樣間隔為1 s,而基于虛擬傳感器預測的排放刷新頻率與發(fā)動機轉速頻率一致。對預測排放結果進行了降頻重采樣,以增強預測數(shù)據(jù)與實測數(shù)據(jù)之間的可比性。

    由圖12可以看出,在定負荷變轉速瞬態(tài)試驗中,隨著發(fā)動機轉速的階躍變化(黑色點劃線),預測NOx值能夠較好地跟隨實測NOx排放的瞬態(tài)變化。在高負荷瞬態(tài)工況下(50 N·m瞬態(tài)工況),預測的NOx與實測數(shù)據(jù)相比有更多的峰值偏差,這主要是由于參與模型構建的試驗工況偏少(15個)造成的,后續(xù)可以通過增大建模數(shù)據(jù)庫來改善。

    圖12 隨速度變化的NOx測量值和預測值

    圖13和圖14分別表示轉速在1 200 r/min和1 000 r/min之間降低和升高過程中實測與預測的NOx排放結果。由圖可知,預測NOx的幅度變化更為尖銳,變化趨勢與實測NOx較為吻合。相較于預測NOx,實測NOx與轉速瞬變輸入之間存在更長的響應延遲,這主要是由排放儀的采樣延遲導致的。排氣經(jīng)導流管傳輸至排放儀分析單元和樣本分析所需的時間決定了實測NOx的延遲程度。基于振動信號開發(fā)的虛擬NOx傳感器響應更快,可以實現(xiàn)NOx排放預測,為瞬態(tài)排放監(jiān)測與控制提供實時排放數(shù)據(jù)。

    圖13 1 200~1 000 r/min-30 N·m工況下NOx測量值和預測值

    圖14 1 000~1 200 r/min-30 N·m工況下NOx測量值和預測值

    圖15顯示了瞬態(tài)試驗過程中全部實測與預測的NOx數(shù)據(jù)。如圖所示,多數(shù)預測 NOx與實測 NOx的相對誤差小于±10%(10%誤差限值標為綠線),其中平均相對誤差為6.5%,最大相對誤差為14.3%,表明所開發(fā)的虛擬傳感器可以較為準確地實現(xiàn)NOx排放水平的瞬態(tài)預測。

    圖15 瞬態(tài)NOx測量值和預測值的關系圖

    4 結論

    1) 缸內壓力曲線重構的相位特性與實測曲線相吻合,幅值基本接近,證實了實測振動的缸內壓力重構方法的合理性。

    2) 通過對比實測與預測NOx排放量,可以得出基于振動信號重構的缸壓信號所建立的回歸模型能夠很好地預測柴油機各種運行工況下的NOx排放。預測值和測量值之間的相關性系數(shù)(R2)為0.971~0.995,誤差小于±10%。瞬態(tài)試驗表明,在與實測NOx變化趨勢較好吻合的前提下,虛擬傳感器的響應速度快于物理傳感器,且瞬態(tài)平均相對誤差低于7%。

    本文中提出的虛擬NOx傳感器驗證了基于實測振動信號預測NOx排放的可行性?;谡駝拥腘Ox排放虛擬傳感器有望為實時駕駛車輛提供一種經(jīng)濟高效的非嵌入式在線監(jiān)測排放技術。

    猜你喜歡
    缸蓋缸內瞬態(tài)
    高壓感應電動機斷電重啟時的瞬態(tài)仿真
    防爆電機(2020年3期)2020-11-06 09:07:36
    缸內直噴發(fā)動機的燃燒和爆震仿真分析
    北京汽車(2018年6期)2019-01-03 11:27:38
    模擬缸蓋工藝在當代內燃機制造業(yè)中的應用
    高性能連體缸蓋關鍵鑄造工藝研究及應用
    缸蓋氣門導管座圈的加工及改進
    十億像素瞬態(tài)成像系統(tǒng)實時圖像拼接
    中國光學(2015年5期)2015-12-09 09:00:39
    鑄造缺陷對A356 鋁合金缸蓋疲勞性能的影響
    上海金屬(2015年5期)2015-11-29 01:13:43
    基于瞬態(tài)流場計算的滑動軸承靜平衡位置求解
    DC/DC變換器中的瞬態(tài)特性分析
    支持直接噴射汽油機歐6開發(fā)目標的缸內和循環(huán)可辨顆粒生成的評估
    少妇丰满av| 久久草成人影院| 青春草亚洲视频在线观看| 免费av毛片视频| 免费av毛片视频| 韩国高清视频一区二区三区| 汤姆久久久久久久影院中文字幕 | 亚洲伊人久久精品综合 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品嫩草影院av在线观看| 老司机影院成人| 精品无人区乱码1区二区| 非洲黑人性xxxx精品又粗又长| 美女高潮的动态| 国产一级毛片在线| 天堂影院成人在线观看| 97超碰精品成人国产| 国产精品99久久久久久久久| 午夜日本视频在线| 国产精品人妻久久久影院| 能在线免费看毛片的网站| 免费在线观看成人毛片| 伦精品一区二区三区| 2022亚洲国产成人精品| 久久久久精品久久久久真实原创| 亚洲国产日韩欧美精品在线观看| 欧美丝袜亚洲另类| 99久国产av精品国产电影| 免费观看的影片在线观看| 免费黄色在线免费观看| 天堂影院成人在线观看| 欧美高清成人免费视频www| 亚洲丝袜综合中文字幕| 免费黄网站久久成人精品| 综合色av麻豆| 亚洲av一区综合| 嫩草影院新地址| 国产免费视频播放在线视频 | 亚洲国产欧美在线一区| 国产成人精品婷婷| 好男人在线观看高清免费视频| 最近最新中文字幕大全电影3| 国产日韩欧美在线精品| 国产精品乱码一区二三区的特点| 91午夜精品亚洲一区二区三区| 欧美日韩精品成人综合77777| 少妇裸体淫交视频免费看高清| 青春草国产在线视频| 亚洲人成网站在线观看播放| av在线观看视频网站免费| 亚洲最大成人av| 久久精品夜色国产| 一个人看视频在线观看www免费| 自拍偷自拍亚洲精品老妇| 听说在线观看完整版免费高清| 中文字幕精品亚洲无线码一区| 成人午夜精彩视频在线观看| 亚洲精品,欧美精品| 久久久精品欧美日韩精品| 午夜激情福利司机影院| 日韩成人伦理影院| 成人三级黄色视频| 欧美zozozo另类| 特大巨黑吊av在线直播| 大香蕉久久网| 97超视频在线观看视频| 久久精品夜色国产| 人妻制服诱惑在线中文字幕| 亚洲欧洲国产日韩| 久久国内精品自在自线图片| 联通29元200g的流量卡| 亚洲乱码一区二区免费版| 一级二级三级毛片免费看| 日本黄色片子视频| 久久久久久久久久久免费av| 在线免费观看的www视频| 国产亚洲精品久久久com| 99热这里只有是精品50| 中文字幕av成人在线电影| 久久久久久久午夜电影| 免费看日本二区| 黄色配什么色好看| 最近最新中文字幕免费大全7| av卡一久久| 日本三级黄在线观看| av福利片在线观看| 亚洲国产欧美人成| 三级毛片av免费| 97超碰精品成人国产| 国产午夜精品久久久久久一区二区三区| 深爱激情五月婷婷| 别揉我奶头 嗯啊视频| 在线a可以看的网站| 久久久久久久久久黄片| 在线免费观看不下载黄p国产| 97超碰精品成人国产| 亚洲自拍偷在线| av在线播放精品| 欧美高清性xxxxhd video| 国产精品1区2区在线观看.| 日日撸夜夜添| 久久人人爽人人片av| 91狼人影院| 亚洲成人中文字幕在线播放| 亚洲精品影视一区二区三区av| 亚洲国产精品国产精品| 亚洲精品自拍成人| av视频在线观看入口| 国产久久久一区二区三区| 亚洲色图av天堂| 色综合亚洲欧美另类图片| 午夜免费激情av| 亚洲三级黄色毛片| 亚洲真实伦在线观看| 日本爱情动作片www.在线观看| 久久精品国产亚洲av天美| 欧美区成人在线视频| 国产欧美另类精品又又久久亚洲欧美| 亚洲一区高清亚洲精品| 精品国产一区二区三区久久久樱花 | 亚洲av成人av| 午夜免费激情av| 国产精品电影一区二区三区| 少妇丰满av| 午夜免费激情av| 99久久精品一区二区三区| 国产亚洲5aaaaa淫片| 69人妻影院| 久久久久久久久久久丰满| 亚洲电影在线观看av| 日本猛色少妇xxxxx猛交久久| 桃色一区二区三区在线观看| 插逼视频在线观看| 国产国拍精品亚洲av在线观看| 国产精品一区www在线观看| 欧美日本视频| 日本wwww免费看| 免费观看的影片在线观看| 国产伦精品一区二区三区四那| 亚洲精品,欧美精品| 色播亚洲综合网| 国产综合懂色| 精品人妻偷拍中文字幕| 亚洲aⅴ乱码一区二区在线播放| 最近最新中文字幕大全电影3| 欧美区成人在线视频| 在线观看av片永久免费下载| 日韩一区二区视频免费看| 天天一区二区日本电影三级| 精品一区二区三区人妻视频| 国产亚洲5aaaaa淫片| 亚洲自偷自拍三级| 国模一区二区三区四区视频| 韩国高清视频一区二区三区| 欧美潮喷喷水| 精品熟女少妇av免费看| 菩萨蛮人人尽说江南好唐韦庄 | 最近视频中文字幕2019在线8| 丰满少妇做爰视频| 久久久国产成人精品二区| 亚洲三级黄色毛片| 欧美日韩综合久久久久久| 91av网一区二区| 我要搜黄色片| 秋霞在线观看毛片| 老师上课跳d突然被开到最大视频| 亚洲国产精品合色在线| 美女脱内裤让男人舔精品视频| 国产爱豆传媒在线观看| 成人国产麻豆网| 亚洲精品国产成人久久av| 丝袜美腿在线中文| 在线免费观看不下载黄p国产| 国产探花在线观看一区二区| 国产精品av视频在线免费观看| 亚洲最大成人手机在线| 国产伦一二天堂av在线观看| 亚洲18禁久久av| 久久99热6这里只有精品| 国产白丝娇喘喷水9色精品| 国产男人的电影天堂91| 国产v大片淫在线免费观看| 久久精品91蜜桃| 97热精品久久久久久| 久久精品国产亚洲网站| 久久久国产成人精品二区| 中文乱码字字幕精品一区二区三区 | 国产在线一区二区三区精 | 午夜福利在线观看吧| 美女被艹到高潮喷水动态| 久久99蜜桃精品久久| 97热精品久久久久久| 国产午夜精品一二区理论片| 国产精品国产三级国产av玫瑰| 中文字幕精品亚洲无线码一区| 99久久九九国产精品国产免费| av又黄又爽大尺度在线免费看 | 永久免费av网站大全| 亚洲av一区综合| 丝袜喷水一区| 国产不卡一卡二| 麻豆成人午夜福利视频| .国产精品久久| 黑人高潮一二区| 亚洲成人久久爱视频| 国产精品女同一区二区软件| 久久精品久久精品一区二区三区| 男人狂女人下面高潮的视频| 少妇的逼水好多| 国产精品综合久久久久久久免费| 精品午夜福利在线看| 直男gayav资源| 天天一区二区日本电影三级| 2021天堂中文幕一二区在线观| 久久久精品欧美日韩精品| 一级黄片播放器| 小蜜桃在线观看免费完整版高清| 两个人的视频大全免费| 床上黄色一级片| 亚洲人与动物交配视频| 免费人成在线观看视频色| 精品一区二区三区人妻视频| 岛国在线免费视频观看| 精品久久久久久久人妻蜜臀av| 亚洲欧洲国产日韩| 免费看av在线观看网站| 久久久久久九九精品二区国产| 丝袜美腿在线中文| 精品国内亚洲2022精品成人| 永久网站在线| av又黄又爽大尺度在线免费看 | 午夜激情欧美在线| 黄色欧美视频在线观看| 国产精品1区2区在线观看.| 91aial.com中文字幕在线观看| 男女国产视频网站| 麻豆成人午夜福利视频| 啦啦啦韩国在线观看视频| 日韩欧美 国产精品| 床上黄色一级片| 国产麻豆成人av免费视频| 国产精品人妻久久久久久| 午夜福利在线在线| 成人漫画全彩无遮挡| 久久久国产成人精品二区| 看免费成人av毛片| 久久婷婷人人爽人人干人人爱| 成人无遮挡网站| 观看免费一级毛片| 大话2 男鬼变身卡| 久久精品综合一区二区三区| 麻豆乱淫一区二区| 国产精品一区二区三区四区免费观看| 国产一区亚洲一区在线观看| 国产成人aa在线观看| www.av在线官网国产| 国产精品三级大全| 久久久精品大字幕| 我的女老师完整版在线观看| 99热这里只有精品一区| 村上凉子中文字幕在线| 三级毛片av免费| 日本免费在线观看一区| 成人无遮挡网站| 特级一级黄色大片| 嫩草影院精品99| 欧美一级a爱片免费观看看| 国产一区二区在线观看日韩| 欧美丝袜亚洲另类| 国产爱豆传媒在线观看| 精品欧美国产一区二区三| 看免费成人av毛片| 淫秽高清视频在线观看| 天堂av国产一区二区熟女人妻| 高清日韩中文字幕在线| 中文字幕制服av| 亚洲av中文av极速乱| 亚洲最大成人手机在线| 一个人看的www免费观看视频| 深爱激情五月婷婷| 国产精品不卡视频一区二区| 亚洲欧洲日产国产| 免费人成在线观看视频色| 亚洲综合色惰| 亚洲精品乱码久久久v下载方式| 日本五十路高清| 国产老妇伦熟女老妇高清| .国产精品久久| 精品免费久久久久久久清纯| 看黄色毛片网站| 黑人高潮一二区| 色尼玛亚洲综合影院| 亚洲乱码一区二区免费版| 精品人妻熟女av久视频| 国产一级毛片七仙女欲春2| 亚洲精品影视一区二区三区av| 欧美激情国产日韩精品一区| 丰满乱子伦码专区| 麻豆一二三区av精品| 高清日韩中文字幕在线| 久久人人爽人人片av| 熟妇人妻久久中文字幕3abv| 日韩欧美精品v在线| 91精品伊人久久大香线蕉| 尤物成人国产欧美一区二区三区| 身体一侧抽搐| 我要看日韩黄色一级片| 亚洲国产成人一精品久久久| 99视频精品全部免费 在线| 国产精品人妻久久久影院| 韩国av在线不卡| 插逼视频在线观看| 噜噜噜噜噜久久久久久91| 蜜桃久久精品国产亚洲av| 大话2 男鬼变身卡| 久久综合国产亚洲精品| 亚洲国产成人一精品久久久| 久久精品91蜜桃| 色综合站精品国产| 国产精品电影一区二区三区| 有码 亚洲区| av黄色大香蕉| 婷婷六月久久综合丁香| av线在线观看网站| 色综合亚洲欧美另类图片| 简卡轻食公司| 看十八女毛片水多多多| 一个人免费在线观看电影| 99久久成人亚洲精品观看| 边亲边吃奶的免费视频| av女优亚洲男人天堂| 午夜福利在线观看免费完整高清在| 六月丁香七月| 在线a可以看的网站| 99久久无色码亚洲精品果冻| av福利片在线观看| 日韩高清综合在线| 你懂的网址亚洲精品在线观看 | 国产av码专区亚洲av| 91精品伊人久久大香线蕉| 黄色配什么色好看| 色播亚洲综合网| 欧美激情久久久久久爽电影| 69av精品久久久久久| 欧美+日韩+精品| 欧美性感艳星| 国产精品99久久久久久久久| 日日干狠狠操夜夜爽| 亚洲在久久综合| 我要看日韩黄色一级片| 亚洲,欧美,日韩| 乱人视频在线观看| 久久久久免费精品人妻一区二区| 亚洲av一区综合| 亚洲性久久影院| 舔av片在线| 97超碰精品成人国产| 91久久精品国产一区二区三区| 亚洲丝袜综合中文字幕| 亚洲精品国产成人久久av| 亚洲人成网站在线观看播放| 18禁在线播放成人免费| 日本熟妇午夜| 少妇猛男粗大的猛烈进出视频 | 一级毛片久久久久久久久女| 亚洲欧美成人精品一区二区| 九九热线精品视视频播放| 久久久精品大字幕| 亚洲欧洲日产国产| 欧美日本视频| 久久99热这里只频精品6学生 | 国产亚洲最大av| 少妇的逼水好多| 国产视频首页在线观看| 麻豆成人av视频| 午夜亚洲福利在线播放| 黄色一级大片看看| 亚洲高清免费不卡视频| 九九在线视频观看精品| 日韩欧美精品免费久久| 免费观看在线日韩| 国产老妇伦熟女老妇高清| 边亲边吃奶的免费视频| 少妇熟女欧美另类| 91av网一区二区| 18禁动态无遮挡网站| www日本黄色视频网| 亚洲欧美精品综合久久99| 久久精品国产自在天天线| 欧美另类亚洲清纯唯美| 国产精品一及| 中文字幕久久专区| or卡值多少钱| 亚洲最大成人中文| 99久久无色码亚洲精品果冻| 在线a可以看的网站| 一本久久精品| 中文天堂在线官网| 秋霞在线观看毛片| 日本午夜av视频| 亚洲久久久久久中文字幕| 男女下面进入的视频免费午夜| 高清毛片免费看| 日韩制服骚丝袜av| 欧美一级a爱片免费观看看| 九色成人免费人妻av| 欧美变态另类bdsm刘玥| 亚洲国产欧洲综合997久久,| 亚洲不卡免费看| 精品人妻视频免费看| 男人狂女人下面高潮的视频| 国产精品av视频在线免费观看| 亚洲av不卡在线观看| 日韩av在线免费看完整版不卡| 亚洲av免费高清在线观看| 国产高潮美女av| 熟女人妻精品中文字幕| 久久99蜜桃精品久久| 国产精品久久久久久久久免| 婷婷六月久久综合丁香| 亚洲在线自拍视频| 久久久精品94久久精品| 国产免费男女视频| 午夜福利在线观看吧| 日本黄色视频三级网站网址| 国产麻豆成人av免费视频| 久久99精品国语久久久| 欧美成人免费av一区二区三区| 在线观看一区二区三区| 国产精品麻豆人妻色哟哟久久 | 亚洲人与动物交配视频| 少妇人妻一区二区三区视频| 亚洲精品乱码久久久v下载方式| av免费在线看不卡| 久久这里有精品视频免费| 中国国产av一级| 搞女人的毛片| 国模一区二区三区四区视频| 纵有疾风起免费观看全集完整版 | 亚洲久久久久久中文字幕| 国产精品一及| 成人午夜精彩视频在线观看| 精品免费久久久久久久清纯| 欧美区成人在线视频| 国产片特级美女逼逼视频| 日日啪夜夜撸| 久久久久久久国产电影| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国语对白做爰xxxⅹ性视频网站| 国内精品美女久久久久久| 热99re8久久精品国产| 国产午夜精品久久久久久一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 日韩制服骚丝袜av| 美女被艹到高潮喷水动态| 欧美一区二区亚洲| 国产高清不卡午夜福利| 精品人妻熟女av久视频| 亚洲成人中文字幕在线播放| 老女人水多毛片| 九色成人免费人妻av| 91狼人影院| 毛片一级片免费看久久久久| 1000部很黄的大片| 能在线免费观看的黄片| 小蜜桃在线观看免费完整版高清| 爱豆传媒免费全集在线观看| 人人妻人人看人人澡| 久久精品久久精品一区二区三区| 人体艺术视频欧美日本| 最近的中文字幕免费完整| 亚洲国产精品sss在线观看| 在线播放国产精品三级| 有码 亚洲区| 欧美激情国产日韩精品一区| 18禁裸乳无遮挡免费网站照片| 边亲边吃奶的免费视频| 亚洲精品456在线播放app| 国产成人91sexporn| 一级毛片aaaaaa免费看小| 一级毛片久久久久久久久女| 人体艺术视频欧美日本| 午夜激情福利司机影院| 久久亚洲精品不卡| 国产精品电影一区二区三区| 最近中文字幕高清免费大全6| 九九在线视频观看精品| 99热全是精品| 精品久久久久久电影网 | 久久精品影院6| 中文资源天堂在线| 能在线免费看毛片的网站| 精品人妻偷拍中文字幕| 老司机福利观看| 国产欧美另类精品又又久久亚洲欧美| 久久久久免费精品人妻一区二区| 欧美成人一区二区免费高清观看| ponron亚洲| 极品教师在线视频| 亚洲国产欧美在线一区| 亚洲在久久综合| 久久久久九九精品影院| av在线观看视频网站免费| 亚洲图色成人| a级毛片免费高清观看在线播放| 国内精品一区二区在线观看| 1024手机看黄色片| 日韩欧美精品v在线| 日韩欧美国产在线观看| 中国国产av一级| 三级男女做爰猛烈吃奶摸视频| 国产精品美女特级片免费视频播放器| 色噜噜av男人的天堂激情| 成人综合一区亚洲| 国产黄色小视频在线观看| 国产高清不卡午夜福利| 久久人人爽人人片av| 久久婷婷人人爽人人干人人爱| 日韩一区二区视频免费看| 久久6这里有精品| 我的老师免费观看完整版| 婷婷色av中文字幕| 国产成人免费观看mmmm| 国产免费又黄又爽又色| 亚洲熟妇中文字幕五十中出| 久久久久性生活片| av在线老鸭窝| 国产一区二区在线观看日韩| 亚洲国产成人一精品久久久| av国产免费在线观看| 免费看av在线观看网站| 午夜福利在线在线| 国产精品国产三级国产专区5o | 久久韩国三级中文字幕| 波野结衣二区三区在线| 六月丁香七月| 亚洲国产高清在线一区二区三| 2021少妇久久久久久久久久久| 国产午夜福利久久久久久| 国产片特级美女逼逼视频| 久久热精品热| 日韩一区二区三区影片| 国产 一区精品| 一卡2卡三卡四卡精品乱码亚洲| 日韩精品有码人妻一区| 99久久人妻综合| 免费av毛片视频| 亚洲在线观看片| 国产一区二区在线av高清观看| 少妇的逼水好多| 欧美性感艳星| 国产毛片a区久久久久| 一区二区三区免费毛片| 欧美日韩一区二区视频在线观看视频在线 | 欧美另类亚洲清纯唯美| 国产一区二区亚洲精品在线观看| 白带黄色成豆腐渣| 午夜激情福利司机影院| 免费电影在线观看免费观看| 国产成人一区二区在线| 亚洲精品乱码久久久久久按摩| 成人鲁丝片一二三区免费| av国产免费在线观看| 国产v大片淫在线免费观看| 国产伦精品一区二区三区视频9| 精品人妻一区二区三区麻豆| 日韩强制内射视频| 国产熟女欧美一区二区| 美女高潮的动态| 亚洲av成人av| 一个人看的www免费观看视频| 99热全是精品| 亚洲精品色激情综合| 人妻制服诱惑在线中文字幕| 国产黄色视频一区二区在线观看 | 小蜜桃在线观看免费完整版高清| 丰满少妇做爰视频| 国产 一区精品| 国产欧美日韩精品一区二区| 国产精品久久视频播放| 国产三级在线视频| 中文字幕久久专区| 亚洲五月天丁香| av在线观看视频网站免费| 精品久久久噜噜| 99国产精品一区二区蜜桃av| 日韩精品青青久久久久久| 久久久a久久爽久久v久久| av在线播放精品| 午夜日本视频在线| 久久久久久久亚洲中文字幕| 69人妻影院| 日韩一本色道免费dvd| 在线播放国产精品三级| 观看美女的网站| 在线观看av片永久免费下载| 老师上课跳d突然被开到最大视频| 天美传媒精品一区二区| 日韩人妻高清精品专区| 最近最新中文字幕大全电影3| 免费观看性生交大片5| 日韩av在线免费看完整版不卡| 观看免费一级毛片| 欧美zozozo另类| 亚洲真实伦在线观看| 亚洲三级黄色毛片| av国产久精品久网站免费入址| 亚洲欧美精品专区久久| 亚洲三级黄色毛片| videossex国产| 精华霜和精华液先用哪个| 久久久欧美国产精品| 男人和女人高潮做爰伦理| 免费播放大片免费观看视频在线观看 | 亚洲av中文av极速乱|