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

    基于機器學習的柴油機納米級微粒預測模型

    2020-04-08 13:33:48李加強王艷艷譚建偉
    科學技術與工程 2020年3期
    關鍵詞:顆粒物模型

    鄒 浪, 何 超*, 李加強, 王艷艷, 譚建偉

    (1.西南林業(yè)大學機械與交通學院,昆明 650224; 2.西南林業(yè)大學云南省高校高原山區(qū)機動車環(huán)保與安全重點實驗室,昆明 650224; 3.北京理工大學機械與車輛學院,北京 100081)

    近年來,機動車尾氣排放已成為大氣污染主要來源,尤其是大氣顆粒物的增加對人體健康影響顯著。研究表明,顆粒物粒徑與呼吸道和心血管疾病的增加以及死亡率之間存在關聯(lián)[1-2]。據報道,柴油車尾氣已成為碳質氣溶膠來源之一[3-5]。因此,減少大氣顆粒物的排放成為當前研究的熱點。

    隨著柴油機排放法規(guī)發(fā)展到歐VI階段,在原來顆粒物質量排放測試的基礎上增加了顆粒物數(shù)量的測量[6]。然而,納米技術的進步使得試驗儀器能夠精確測量1~1 000 nm顆粒物粒徑及數(shù)量,從而降低偶然誤差和系統(tǒng)誤差帶來的影響。研究發(fā)現(xiàn),機動車輛燃燒產生的顆粒對人體的危害主要來自細顆粒物和超細顆粒物[7-10],尤其是直徑小于100 nm的超細顆粒物對人和動物的危害比細顆粒物更加突出[11]。相比汽油車,柴油車具有使用壽命長久、燃油經濟性好、壓縮比高、低污染等優(yōu)點,被作為動力裝置安裝在許多歐美轎車及重型車輛上。由此可見,改善柴油車的污染物排放對控制大氣污染起著重要的作用。

    針對發(fā)動機轉速、空燃比、氣缸壓力、廢棄再循環(huán)(exhaust gas re-circulation,EGR)率等燃燒特性與納米級PM的研究已有大量報道[12-16]。因此,利用發(fā)動機燃燒產生的缸壓循環(huán)及空燃比等物理參數(shù)預測污染物的排放必不可少。Henningsson等[17]利用柴油機在不同速度、噴油定時及EGR率等條件下氣缸壓力對NOx和尾氣中的氧氣濃度λ進行了有效預測;Pu等[14]利用氣缸壓力作為神經網絡的輸入,預測了粒徑為15~1 000 nm的顆粒物分布;Weymann等[18]建立燃燒壓力軌跡的神經網絡模型,結果表明該模型的預測值與實際值的相關系數(shù)為0.992;余林嘯[19]研究不同海拔高度對柴油機的燃燒特性和排放特性的影響,結果表明隨著海拔的升高,柴油機缸內最高燃燒壓力逐漸減小,核模態(tài)和積聚模態(tài)顆粒物數(shù)量增加。

    然而,在不同海拔、轉速及轉距情況下,利用柴油機燃燒產生的氣缸壓力預測納米級顆粒物(particulate matter,PM)在相關領域的文獻非常有限。由于實驗條件非常復雜,需要測試大量原始數(shù)據,導致模型的構建十分困難。基于此,運用機器學習領域的方法,采用主成分分析(principal component analysis, PCA)研究不同海拔下氣缸壓力的相關特征。隨后訓練神經網絡,以主成分訓練數(shù)據作為輸入,顆粒物排放的粒徑數(shù)量濃度作為輸出,從而比較不同主成分對于顆粒物濃度預測情況,選擇合適的主成分建立模型是研究的重點。

    1 試驗設備

    試驗用柴油機為錫柴CA6DF3,試驗發(fā)動機規(guī)格如表1所示, 燃料為0號柴油。

    表1 試驗發(fā)動機參數(shù)

    采用CW440D型電渦流測功機實時測量發(fā)動機的輸出轉矩和轉速。為了得到精確的燃燒特性參數(shù),使用DEWETRON 5000型燃燒分析儀對角標儀與缸壓傳感器所測試的數(shù)據進行實時現(xiàn)場測量,按每0.2曲柄角間隔采集1個氣缸壓力數(shù)據。同時,使用低壓靜電沖擊式采樣器ELPI對不同粒徑的顆粒物數(shù)量濃度進行測量,第一級沖擊器中加入濾紙,使測量范圍擴展到7 nm粒徑,粒徑分布如表2所示。

    表2 ELPI粒徑分級

    2 試驗方法

    2.1 試驗設計

    試驗對280種工況進行測試,其中包括4種海拔(0、1 608、2 408、3 300 m),每種海拔下8種轉速(1 030、1 200、1 400、1 600、1 800、2 100、2 300 r/min及怠速700 r/min),以及每種轉速下10個不同的轉矩(怠速除外),相同海拔下轉速與轉矩的關系如圖1所示,其中轉矩為47~790 nm。每個工況點的轉矩和轉速保持一致,待發(fā)動機工作穩(wěn)定后,再進行數(shù)據采集。經過多次測試,對每個工況隨機選取50個周期取平均值,收集大量氣缸壓力及不同粒徑的顆粒物濃度數(shù)據。

    圖1 相同海拔下轉速與轉矩關系

    2.2 主成分分析原理

    主成分分析是一種大數(shù)據降維的方法[20],通過查找變量之間的聯(lián)系,使原始變量重新線性組合得到一組不相關的綜合指標,這為降低數(shù)據的維數(shù)提供了方法。因此,主城分分析在圖像識別、模型預測及計算機等領域運用廣泛。主成分分析原理公式如式(1)所示:

    (1)

    式(1)中:A∈RN是原始數(shù)據集,N為原始變量的維數(shù),即不同海拔、轉速和轉矩條件下,N變量的氣缸壓力數(shù)據組合成原始數(shù)據集;X=[x1,x2,…,xL]為主成分組成的基向量矩陣;C=[c1,c2,…,cL]為是由L個系數(shù)組成的向量矩陣;e為誤差。

    當L?N,使誤差e最小化時,可由式得到式(2),其中T為轉置。這表明C的每個分量即主成分系數(shù)能表示相對應原始數(shù)據的特征。

    C=AXT或X=CTA

    (2)

    由于奇異值分解不僅是譜分析的推廣運用,還能消除干擾信號帶來的影響。因此,原始數(shù)據使用MATLAB svd函數(shù)進行奇異值分解,其原理如式(3)所示:

    A=USVT

    (3)

    式(3)中:A為任意實維矩陣;U和V為正交矩陣;S為奇異值矩陣,奇異值按列從左到右逐漸減小;U和V的列向量分別是左奇異向量和右奇異向量,其中U的列向量是奇異值從大到小依次對應的特征向量。若奇異值矩陣的某部分對角線元素遠遠大于其他奇異值之和,則這部分元素對應的特征向量代表原始數(shù)據的大部分信息。

    2.3 神經網絡結構

    對影響指標A=[a1,a2,…,an]進行主成分分析,其結果X=[X1,X2,…,Xn]作為神經網絡的輸入向量,得到對應的輸出結果y,如圖2所示。wij和wih分別為輸入層與隱含層、隱含層與輸出層的連接權值。通過設定期望誤差來評估實模型的精確程度,如果實際輸出與預測輸出之間的誤差大于期望誤差,則進行多次逆向傳播修正,直至小于期望誤差終止,最終建立神經網絡模型。

    圖2 神經網絡結構原理

    3 試驗結果與分析

    3.1 主成分分析氣缸壓力

    發(fā)動機燃燒過程集中在活塞上止點附近,這時氣缸壓力變化趨勢最能代表整個周期的燃燒特性。由此可見,選擇上止點-30~45 ℃曲柄角進行數(shù)據研究是降低實驗難度的可行方案[14]。每循環(huán)采樣1次,采樣間隔為0.2曲柄角,可獲得375個樣本點,最終對10 500個數(shù)據進行分析。

    通過MATLAB奇異值分解并計算氣缸壓力的前10個主成分貢獻率如表3所示。圖3所示為X1、X4、X7、X10主成分的氣缸壓力變化曲線,其中橫坐標表示曲柄角,縱坐標表示數(shù)據預處理后的缸壓。由表3、圖3可知,第1主成分提取原始數(shù)據62.02%的特征信息,而主成分10主要捕捉活塞上止點附近壓力變化的頻率。圖4所示為不同海拔下X1、X4、X7、X10主成分的重構數(shù)據與原始數(shù)據的分布,表明主成分維數(shù)越高,奇異值分解重構數(shù)據與原始數(shù)據信息重合度越好,說明奇異值分解降低了原始數(shù)據特征信息的損失。這是由于PCA的本質是空間坐標的變換且數(shù)據結構保持不變,獲得的主成分集是原始變量的線性組合造成的結果。

    表3 主成分貢獻率

    圖3 部分主成分的缸壓變化曲線

    對280組原始數(shù)據和重構數(shù)據進行均方根誤差(RMSE)檢驗,結果如圖5所示。在圖5中,第1主成分因貢獻率小于80%造成均方根誤差過大,導致信息提取過少,因此不能用于分析模型的訓練數(shù)據。圖3~圖5表明,前4個主成分代表缸內壓力總體變化趨勢,且隨著主成分貢獻率的增大,均方根誤差急劇減小,當達到一定程度后誤差緩慢降低。當達到10個主成分時,均方根誤差為0.014~0.082 MPa,信息提取效果較好。綜上所述,不同工況下氣缸壓力變化趨勢能由主成分系數(shù)表示。

    圖4 主成分重構數(shù)據與原始數(shù)據的缸壓分布

    3.2 神經網絡預測模型分析

    3.2.1 神經網絡模型的設計

    由于柴油機的顆粒物排放受負荷、溫度、轉速、噴油時間、火花正時、燃料及EGR率等物理參數(shù)的影響。因此建立氣缸壓力與顆粒物的排放模型屬于非線性映射。顯然,一個多層反向傳播的前饋神經網絡是非線性映射模型仿真的最佳選擇之一。鑒于原始數(shù)據涉及不同的量綱,為保證預測模型精確度,必須對數(shù)據進行標準化處理。

    選擇單隱含層的神經網絡模型,輸入是經過前4、7、10個主成分分析的280組標準化數(shù)據,輸出為具有6個神經元的粒徑為7~990 nm顆粒物濃度(即第一組粒徑為7~29 nm;第二組粒徑為29~57 nm;第三組粒徑為57~101 nm;第四組粒徑為101~255 nm;第五組粒徑為255~990 nm,第六組粒徑為7~990 nm)。其中隨機產生238組數(shù)據用于神經網絡訓練,剩余42組數(shù)據用于神經網絡測試。主成分分析與神經網絡模型流程如圖6所示。

    圖5 氣缸壓力均方根誤差曲線

    圖6 主成分分析與神經網絡模型流程

    由于神經網絡極易陷入局部最優(yōu)化,而貝葉斯正則化反向傳播算法能使神經網絡模型具有最佳整體性能。因此,選擇MATLAB定義的貝葉斯函數(shù)(trainbr)對該模型進行訓練優(yōu)化。最大訓練次數(shù)為1 000次,訓練目標為0.001,學習率設置為0.01。

    3.2.2 模型預測結果分析

    為研究不同工況下不同主成分的氣缸壓力與PM粒徑濃度模型的最佳結合點,需要對神經網絡模型進行反復測試和訓練。其中平均絕對誤差(MAE)和均方根誤差(RSME)越小,回歸系數(shù)R2越接近1,則表明該模型預測效果越好。

    經過多次訓練后,X1~X4主成分預測PM粒徑數(shù)量濃度的神經網絡采取4個輸入層神經元,隱含層節(jié)點為5時,該模型能準確預測7~990 nm的顆粒物粒徑濃度,測試的R2達到0.85,整體性能R2達到0.81,如圖7所示,其中橫坐標從左到右依次表示為訓練目標、測試目標和全局輸出數(shù)據,縱坐標表示神經網絡對應的預測輸出。此外,X1~X7主成分選用7個輸入神經元,8個隱含層節(jié)點時模型預測穩(wěn)定高效,測試的R2達到0.86,整體性能R2達到0.91,如圖8所示;X1~X10主成分采用10個輸入變量,9個隱含層神經元,測試的R2達到0.92,整體性能R2達到0.95,如圖9所示。

    圖7 X1~X4回歸分析

    圖8 X1~X7回歸分析

    對X1~X4、X1~X7和X1~X10三個模型進行仿真測試,引入平均絕對誤差(MAE)和均方根誤差(RSME)對樣本預測值與實測值進行檢驗,得到結果如表4所示。由表4可知,模型的平均絕對誤差為90.74~2.284 29×103cm3,均方根誤差為1.612×104~9.769×104cm3,這表明柴油機燃燒產生的氣缸壓力能準確預測7~990 nm每組顆粒物的數(shù)量濃度分布。由表4可知,平均絕對誤差和均方根誤差隨著主成分貢獻率的增大而逐漸降低,表明主成分占比與模型的預測精度成正相關,這是由于新加入的主成分相比舊主成分對于原始數(shù)據的特征提取更全面,引起主成分與不同粒徑的顆粒物濃度的相關性增加的結果。其中7~990 nm PM數(shù)濃度在X1~X10模型的樣本預測值與試驗值如圖10所示。

    圖9 X1~X10回歸分析圖

    研究表明在不同行駛工況下,當主成分達到一定比例時,利用柴油機燃燒產生的氣缸壓力提取主成分來預測PM的排放是一種有效方法。此外,前10主成分的缸壓預測7~990 nm顆粒物的平均絕對誤差為90.74 cm3,均方根誤差為1.612×104cm3,回歸系數(shù)R2為0.95,表明在某些復雜實驗條件下測量顆粒物的粒徑濃度排放時,可以獨立測量氣缸壓力進行預測。

    表4 每組粒徑濃度的試驗值和預測值的誤差

    圖10 PM數(shù)濃度在X1~X10模型的試驗值與樣本預測值

    4 結論

    通過主成分分析與神經網絡組合模型,得到以下結論。

    (1)奇異值分解能有效減少原始數(shù)據信息的損失。

    (2)在運用主成分方法分析氣缸壓力時,不同工況下的氣缸壓力變化能夠由主成分系數(shù)表示。

    (3)利用主成分分析和神經網絡建模,當主成分達到一定比例時,實測數(shù)據與模型的預測值具有良好的一致性,其回歸系數(shù)R2達到0.95。證明主成分分析與神經網絡結合是一種降低數(shù)據維度和模型構造難度的可行方法之一。

    (4)柴油機氣缸壓力與在7~990 nm顆粒物的數(shù)量濃度預測分布擬合效果較好,表明發(fā)動機燃燒產生缸壓力可以用來標定不同工況下車輛顆粒物的排放。

    猜你喜歡
    顆粒物模型
    一半模型
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權M-估計的漸近分布
    南平市細顆粒物潛在來源分析
    固定源細顆粒物監(jiān)測技術現(xiàn)狀分析與思考
    3D打印中的模型分割與打包
    錯流旋轉填料床脫除細顆粒物研究
    化工進展(2015年3期)2015-11-11 09:18:15
    FLUKA幾何模型到CAD幾何模型轉換方法初步研究
    多層介質阻擋放電處理柴油機尾氣顆粒物
    借鑒歐洲經驗加快我國顆粒物污染防治
    亚洲欧美成人精品一区二区| 午夜免费鲁丝| 欧美精品国产亚洲| av免费在线看不卡| 大片电影免费在线观看免费| 多毛熟女@视频| 国产视频首页在线观看| 在线观看国产h片| av在线播放精品| 少妇猛男粗大的猛烈进出视频| 综合色丁香网| av免费观看日本| 中文字幕人妻丝袜制服| 一个人免费看片子| 丰满饥渴人妻一区二区三| 精品一区二区免费观看| 婷婷色麻豆天堂久久| 成年美女黄网站色视频大全免费| 日韩,欧美,国产一区二区三区| 久久毛片免费看一区二区三区| 久久久久久久亚洲中文字幕| 999精品在线视频| 中国国产av一级| 男女高潮啪啪啪动态图| 亚洲av福利一区| 美女国产视频在线观看| 亚洲伊人久久精品综合| 99热网站在线观看| www.色视频.com| 久久这里只有精品19| 久久这里只有精品19| 亚洲成国产人片在线观看| 亚洲国产精品专区欧美| 亚洲美女视频黄频| 91aial.com中文字幕在线观看| 免费人妻精品一区二区三区视频| 成人毛片a级毛片在线播放| 久久鲁丝午夜福利片| 99国产精品免费福利视频| 大香蕉久久网| 日韩制服丝袜自拍偷拍| 婷婷成人精品国产| 少妇 在线观看| 大香蕉97超碰在线| 欧美另类一区| 我的女老师完整版在线观看| 日韩制服骚丝袜av| 国产精品成人在线| 国产精品国产三级专区第一集| 亚洲美女视频黄频| 91在线精品国自产拍蜜月| 欧美人与性动交α欧美软件 | 亚洲伊人色综图| 午夜免费观看性视频| a级片在线免费高清观看视频| 亚洲欧美精品自产自拍| 成人黄色视频免费在线看| 亚洲内射少妇av| 天美传媒精品一区二区| 国产精品久久久久久精品电影小说| 欧美xxxx性猛交bbbb| 内地一区二区视频在线| 色5月婷婷丁香| 一级毛片黄色毛片免费观看视频| 母亲3免费完整高清在线观看 | 亚洲国产精品专区欧美| 国产欧美另类精品又又久久亚洲欧美| 两个人看的免费小视频| av一本久久久久| 久久精品国产综合久久久 | 在线免费观看不下载黄p国产| 久久精品国产鲁丝片午夜精品| 精品一区二区三卡| 91国产中文字幕| 9热在线视频观看99| 国产精品三级大全| 少妇被粗大猛烈的视频| 久久女婷五月综合色啪小说| 国产亚洲av片在线观看秒播厂| 国产精品一区www在线观看| 亚洲国产欧美日韩在线播放| 亚洲精品一区蜜桃| 一区二区日韩欧美中文字幕 | 中文天堂在线官网| 中文字幕人妻熟女乱码| 国产精品成人在线| 一级毛片黄色毛片免费观看视频| 宅男免费午夜| 女人被躁到高潮嗷嗷叫费观| 欧美人与性动交α欧美软件 | 午夜免费观看性视频| 欧美3d第一页| 国产精品久久久久久久久免| 亚洲精品国产av蜜桃| 亚洲精品一区蜜桃| 国产伦理片在线播放av一区| 90打野战视频偷拍视频| 尾随美女入室| 在线精品无人区一区二区三| 亚洲丝袜综合中文字幕| 纯流量卡能插随身wifi吗| 热re99久久国产66热| 人人澡人人妻人| 99久久综合免费| av在线播放精品| 国产成人av激情在线播放| 97在线视频观看| 午夜久久久在线观看| 一级片'在线观看视频| 自拍欧美九色日韩亚洲蝌蚪91| www.色视频.com| 国产熟女午夜一区二区三区| 国产永久视频网站| 亚洲精品aⅴ在线观看| 男女免费视频国产| 成年美女黄网站色视频大全免费| 久久国产精品大桥未久av| 女人久久www免费人成看片| 大香蕉久久网| 欧美激情极品国产一区二区三区 | 在线观看国产h片| 美女xxoo啪啪120秒动态图| 激情五月婷婷亚洲| 免费av中文字幕在线| 久久久国产一区二区| 精品久久蜜臀av无| 亚洲av电影在线观看一区二区三区| 九色成人免费人妻av| 亚洲欧美清纯卡通| 久久久精品区二区三区| 王馨瑶露胸无遮挡在线观看| 国产1区2区3区精品| 免费大片18禁| av国产精品久久久久影院| av国产精品久久久久影院| 卡戴珊不雅视频在线播放| 美女福利国产在线| 高清欧美精品videossex| 国产片特级美女逼逼视频| 熟女人妻精品中文字幕| 亚洲欧美一区二区三区国产| 欧美xxxx性猛交bbbb| 秋霞在线观看毛片| 91精品三级在线观看| 国产片特级美女逼逼视频| 国产黄色视频一区二区在线观看| 高清欧美精品videossex| 成年人午夜在线观看视频| a 毛片基地| 欧美日韩av久久| a 毛片基地| 久久人人爽人人片av| 亚洲精品成人av观看孕妇| 一本大道久久a久久精品| 亚洲,一卡二卡三卡| 日日爽夜夜爽网站| 亚洲,一卡二卡三卡| 久久精品国产综合久久久 | 日本黄色日本黄色录像| 人人妻人人澡人人爽人人夜夜| 国产av精品麻豆| 肉色欧美久久久久久久蜜桃| 国产成人精品久久久久久| a级毛片在线看网站| av在线观看视频网站免费| 美女主播在线视频| 日韩制服骚丝袜av| 国产男人的电影天堂91| 一级毛片电影观看| 综合色丁香网| 大片免费播放器 马上看| av电影中文网址| 久久久久网色| 亚洲欧美日韩另类电影网站| 亚洲成人一二三区av| 亚洲伊人色综图| 日本91视频免费播放| 午夜免费观看性视频| 高清视频免费观看一区二区| 成人亚洲精品一区在线观看| 亚洲欧美精品自产自拍| 国产精品 国内视频| 欧美精品一区二区免费开放| 99热国产这里只有精品6| 亚洲一级一片aⅴ在线观看| 国产男女内射视频| 青青草视频在线视频观看| 好男人视频免费观看在线| 国产片特级美女逼逼视频| 激情视频va一区二区三区| 国产爽快片一区二区三区| 欧美97在线视频| 亚洲精品一二三| 成人免费观看视频高清| 狠狠婷婷综合久久久久久88av| 美女大奶头黄色视频| av播播在线观看一区| 亚洲国产精品专区欧美| 一边亲一边摸免费视频| 9色porny在线观看| 一级毛片我不卡| 侵犯人妻中文字幕一二三四区| 狠狠精品人妻久久久久久综合| 天天影视国产精品| 亚洲,欧美,日韩| 你懂的网址亚洲精品在线观看| 亚洲精品日韩在线中文字幕| 成人国产麻豆网| 国产精品免费大片| 精品少妇久久久久久888优播| 夜夜骑夜夜射夜夜干| 熟女人妻精品中文字幕| 纯流量卡能插随身wifi吗| 精品视频人人做人人爽| 免费在线观看完整版高清| 亚洲av电影在线进入| 亚洲一区二区三区欧美精品| 亚洲美女视频黄频| 色吧在线观看| 青春草亚洲视频在线观看| 伦理电影免费视频| 免费观看无遮挡的男女| 国产精品久久久久成人av| 人妻系列 视频| 国产乱人偷精品视频| 少妇被粗大的猛进出69影院 | 国产熟女午夜一区二区三区| 777米奇影视久久| av国产精品久久久久影院| 曰老女人黄片| 在线亚洲精品国产二区图片欧美| 国产淫语在线视频| 国产白丝娇喘喷水9色精品| 日韩av不卡免费在线播放| 人人妻人人澡人人看| av一本久久久久| 国产精品久久久久成人av| 亚洲精品国产av成人精品| 最近2019中文字幕mv第一页| 黄网站色视频无遮挡免费观看| 女人被躁到高潮嗷嗷叫费观| 日本91视频免费播放| 精品人妻熟女毛片av久久网站| 亚洲天堂av无毛| 热re99久久国产66热| 亚洲少妇的诱惑av| 女性生殖器流出的白浆| 成人国产麻豆网| 亚洲国产最新在线播放| 亚洲国产精品国产精品| 亚洲经典国产精华液单| 久久精品久久久久久噜噜老黄| 男女高潮啪啪啪动态图| 日韩一区二区三区影片| av视频免费观看在线观看| 欧美xxⅹ黑人| www.熟女人妻精品国产 | 国产精品成人在线| 人人澡人人妻人| 在线 av 中文字幕| 18禁观看日本| 日韩一区二区视频免费看| 亚洲精品自拍成人| 成人国语在线视频| 亚洲情色 制服丝袜| 亚洲欧美清纯卡通| 涩涩av久久男人的天堂| 日韩人妻精品一区2区三区| 十分钟在线观看高清视频www| 有码 亚洲区| 久久97久久精品| 亚洲一码二码三码区别大吗| 街头女战士在线观看网站| 精品一区二区免费观看| 亚洲av成人精品一二三区| 久久人人爽人人爽人人片va| 91在线精品国自产拍蜜月| 亚洲国产最新在线播放| 日韩一区二区三区影片| 极品少妇高潮喷水抽搐| 欧美变态另类bdsm刘玥| 老司机影院毛片| 日韩三级伦理在线观看| av在线播放精品| 国产老妇伦熟女老妇高清| 国产成人欧美| 日韩视频在线欧美| 久久久久精品人妻al黑| 寂寞人妻少妇视频99o| 99香蕉大伊视频| 一级片'在线观看视频| 免费观看无遮挡的男女| 男女午夜视频在线观看 | 男男h啪啪无遮挡| 午夜激情av网站| 极品少妇高潮喷水抽搐| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 精品少妇久久久久久888优播| 一区二区三区乱码不卡18| 一边摸一边做爽爽视频免费| 十八禁高潮呻吟视频| 色5月婷婷丁香| 永久网站在线| 最黄视频免费看| 国产视频首页在线观看| 婷婷色综合www| 欧美 日韩 精品 国产| 精品人妻在线不人妻| 国产高清不卡午夜福利| 2022亚洲国产成人精品| av又黄又爽大尺度在线免费看| 中文字幕av电影在线播放| 成年av动漫网址| 国产有黄有色有爽视频| 美女国产视频在线观看| 最后的刺客免费高清国语| 亚洲国产成人一精品久久久| 在线观看免费视频网站a站| 欧美激情极品国产一区二区三区 | 免费黄色在线免费观看| 亚洲第一av免费看| 成人免费观看视频高清| 黄片无遮挡物在线观看| 免费久久久久久久精品成人欧美视频 | 高清视频免费观看一区二区| 日本免费在线观看一区| 日韩精品有码人妻一区| 欧美精品av麻豆av| 性色avwww在线观看| 欧美 日韩 精品 国产| 久久这里有精品视频免费| 国产av精品麻豆| 国产极品粉嫩免费观看在线| 亚洲精品国产色婷婷电影| av不卡在线播放| av黄色大香蕉| 性高湖久久久久久久久免费观看| 日产精品乱码卡一卡2卡三| 亚洲内射少妇av| 亚洲国产最新在线播放| 精品99又大又爽又粗少妇毛片| 边亲边吃奶的免费视频| 亚洲精品色激情综合| 少妇猛男粗大的猛烈进出视频| 最近2019中文字幕mv第一页| 日本91视频免费播放| 亚洲欧美一区二区三区黑人 | 男女啪啪激烈高潮av片| 欧美精品一区二区免费开放| 夫妻午夜视频| 亚洲国产色片| 男女午夜视频在线观看 | 七月丁香在线播放| 毛片一级片免费看久久久久| 如何舔出高潮| 巨乳人妻的诱惑在线观看| 熟女电影av网| 在线观看三级黄色| 又粗又硬又长又爽又黄的视频| 观看美女的网站| 一二三四在线观看免费中文在 | 欧美激情极品国产一区二区三区 | 国产精品久久久久久精品电影小说| 99久久人妻综合| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 高清欧美精品videossex| 宅男免费午夜| 这个男人来自地球电影免费观看 | 欧美97在线视频| 视频在线观看一区二区三区| 大香蕉97超碰在线| 大香蕉久久网| 99热国产这里只有精品6| 最近中文字幕高清免费大全6| 另类精品久久| 久久精品国产亚洲av涩爱| 如日韩欧美国产精品一区二区三区| 激情五月婷婷亚洲| 侵犯人妻中文字幕一二三四区| 精品人妻在线不人妻| 日韩av不卡免费在线播放| 精品久久国产蜜桃| www日本在线高清视频| 巨乳人妻的诱惑在线观看| 在线看a的网站| xxx大片免费视频| 中文乱码字字幕精品一区二区三区| 日本欧美视频一区| 满18在线观看网站| 爱豆传媒免费全集在线观看| 欧美 亚洲 国产 日韩一| 在线观看人妻少妇| 午夜福利在线观看免费完整高清在| 美女国产视频在线观看| 亚洲精品国产av成人精品| 国产高清三级在线| 男女免费视频国产| 国产片内射在线| 91成人精品电影| 亚洲天堂av无毛| 男男h啪啪无遮挡| 日日啪夜夜爽| 欧美日韩一区二区视频在线观看视频在线| 国产成人免费观看mmmm| 亚洲欧美清纯卡通| 成人综合一区亚洲| 久久免费观看电影| 夫妻午夜视频| 人妻系列 视频| 91精品三级在线观看| 精品国产乱码久久久久久小说| av卡一久久| 午夜福利视频在线观看免费| 寂寞人妻少妇视频99o| 亚洲内射少妇av| 久久久久精品性色| 蜜桃在线观看..| 波多野结衣一区麻豆| 91精品伊人久久大香线蕉| 精品亚洲成a人片在线观看| 在线看a的网站| av一本久久久久| 成人二区视频| 黄色配什么色好看| 美女国产视频在线观看| 欧美国产精品va在线观看不卡| 久久精品aⅴ一区二区三区四区 | 巨乳人妻的诱惑在线观看| 哪个播放器可以免费观看大片| 久久久精品区二区三区| 国产在线免费精品| 亚洲av男天堂| 曰老女人黄片| 老司机影院毛片| 日本爱情动作片www.在线观看| 亚洲图色成人| 国产成人精品无人区| 国产日韩一区二区三区精品不卡| 18禁动态无遮挡网站| 欧美日韩综合久久久久久| 精品久久蜜臀av无| 亚洲精品第二区| 国产在线免费精品| 免费不卡的大黄色大毛片视频在线观看| 曰老女人黄片| 欧美+日韩+精品| www.色视频.com| 9热在线视频观看99| 深夜精品福利| 美女国产视频在线观看| 亚洲欧美中文字幕日韩二区| 深夜精品福利| 观看美女的网站| 国产成人午夜福利电影在线观看| 一区在线观看完整版| 国产在视频线精品| 日韩一区二区视频免费看| 亚洲色图 男人天堂 中文字幕 | 亚洲av综合色区一区| 久久久久久久久久人人人人人人| 九色成人免费人妻av| 亚洲精品456在线播放app| 日本欧美国产在线视频| 午夜福利视频在线观看免费| 亚洲av男天堂| 乱码一卡2卡4卡精品| 亚洲综合色惰| 久久久精品区二区三区| 国产精品一国产av| 黄色 视频免费看| 精品国产乱码久久久久久小说| 久久精品国产亚洲av天美| 人人妻人人澡人人看| 成人二区视频| 日本-黄色视频高清免费观看| 亚洲婷婷狠狠爱综合网| 国产欧美亚洲国产| 中文字幕av电影在线播放| 成人亚洲精品一区在线观看| 久久精品国产综合久久久 | 免费日韩欧美在线观看| 亚洲国产最新在线播放| 最近中文字幕2019免费版| 啦啦啦视频在线资源免费观看| 精品一品国产午夜福利视频| 草草在线视频免费看| 青春草国产在线视频| 国产免费福利视频在线观看| 咕卡用的链子| 久久婷婷青草| √禁漫天堂资源中文www| 午夜免费观看性视频| 中文字幕最新亚洲高清| 精品人妻熟女毛片av久久网站| 久久久久人妻精品一区果冻| 久久久久久久亚洲中文字幕| 国产一区二区在线观看日韩| 亚洲精品日本国产第一区| 99精国产麻豆久久婷婷| 丰满少妇做爰视频| av.在线天堂| 亚洲精品国产av成人精品| 日本-黄色视频高清免费观看| 亚洲国产精品成人久久小说| 国产精品 国内视频| 欧美成人午夜精品| 永久免费av网站大全| 97人妻天天添夜夜摸| 久久久亚洲精品成人影院| 亚洲人与动物交配视频| 亚洲精品自拍成人| 黄网站色视频无遮挡免费观看| 大香蕉久久成人网| 国产成人欧美| 日本vs欧美在线观看视频| 丝袜喷水一区| 美女国产视频在线观看| 午夜福利视频在线观看免费| 少妇的逼好多水| 免费看不卡的av| 在线观看国产h片| 捣出白浆h1v1| 99久久综合免费| 成年美女黄网站色视频大全免费| 曰老女人黄片| 亚洲国产精品一区二区三区在线| 久久久久网色| freevideosex欧美| 国产乱来视频区| 制服诱惑二区| 国产在线一区二区三区精| 久久久久精品久久久久真实原创| 亚洲伊人色综图| 亚洲综合色网址| 久久亚洲国产成人精品v| 一级毛片电影观看| 18禁国产床啪视频网站| 国产综合精华液| 五月开心婷婷网| 咕卡用的链子| 性高湖久久久久久久久免费观看| 国产免费一区二区三区四区乱码| 蜜桃在线观看..| 成人18禁高潮啪啪吃奶动态图| 久久精品国产亚洲av涩爱| 国产免费视频播放在线视频| 肉色欧美久久久久久久蜜桃| 欧美精品亚洲一区二区| 久久99精品国语久久久| 香蕉精品网在线| 少妇人妻久久综合中文| 男女午夜视频在线观看 | 自线自在国产av| 人人妻人人澡人人看| 欧美亚洲日本最大视频资源| 日本午夜av视频| 免费黄网站久久成人精品| 亚洲综合色惰| 日本欧美视频一区| 成人影院久久| 国产一区二区三区综合在线观看 | 精品人妻熟女毛片av久久网站| 美女福利国产在线| 你懂的网址亚洲精品在线观看| 成年人午夜在线观看视频| 街头女战士在线观看网站| 亚洲精品国产av蜜桃| 国产精品秋霞免费鲁丝片| 国产欧美日韩一区二区三区在线| 欧美精品人与动牲交sv欧美| 国产片特级美女逼逼视频| 大话2 男鬼变身卡| 亚洲四区av| 亚洲欧美一区二区三区黑人 | 婷婷色综合www| 日本免费在线观看一区| 在线免费观看不下载黄p国产| 一级片免费观看大全| 熟女人妻精品中文字幕| 成人毛片a级毛片在线播放| av女优亚洲男人天堂| 国产男女超爽视频在线观看| av在线播放精品| 国产在线免费精品| av在线播放精品| www.av在线官网国产| 中文字幕av电影在线播放| 欧美xxⅹ黑人| 狂野欧美激情性bbbbbb| 亚洲精品视频女| 韩国av在线不卡| 欧美精品一区二区免费开放| 热re99久久精品国产66热6| 欧美日韩视频精品一区| 久热久热在线精品观看| 久久精品久久久久久噜噜老黄| 欧美 日韩 精品 国产| 免费播放大片免费观看视频在线观看| 午夜精品国产一区二区电影| 日本色播在线视频| 在线观看一区二区三区激情| 亚洲成国产人片在线观看| 精品国产一区二区三区四区第35| 日韩av免费高清视频| 国产成人午夜福利电影在线观看| 欧美精品av麻豆av| 亚洲,欧美,日韩| 高清不卡的av网站| 99精国产麻豆久久婷婷| 99久久精品国产国产毛片| 交换朋友夫妻互换小说| 韩国av在线不卡| 亚洲av.av天堂| 亚洲丝袜综合中文字幕| 搡老乐熟女国产|