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

    基于K-均值聚類分析的城市道路汽車行駛工況構建方法研究*

    2017-12-15 03:37:43彭育輝楊輝寶李孟良喬學齊
    汽車技術 2017年11期
    關鍵詞:數(shù)據(jù)源特征參數(shù)運動學

    彭育輝 楊輝寶 李孟良 喬學齊

    (1.福州大學,福州 350116;2.中國汽車技術研究中心,天津 300300;3.廈門金龍旅行車有限公司,廈門 361006)

    基于K-均值聚類分析的城市道路汽車行駛工況構建方法研究*

    彭育輝1,2楊輝寶1李孟良2喬學齊3

    (1.福州大學,福州 350116;2.中國汽車技術研究中心,天津 300300;3.廈門金龍旅行車有限公司,廈門 361006)

    以實時采集的乘用車行駛數(shù)據(jù)為數(shù)據(jù)源,進行了城市道路汽車行駛工況構建方法的研究。分別運用運動學片段分析法、主成分分析法和K均值聚類分析法對實測數(shù)據(jù)進行降維和分類,提出以Silhouette函數(shù)實現(xiàn)對聚類結果的篩選,以減少人為選擇的誤差,并根據(jù)聚類中心的大小篩選所需運動學片段構建候選工況。在目標代表工況的遴選方面,提出了綜合6個特征參數(shù)和最大SAFD差異值的評價標準。最后通過試驗驗證了該行駛工況構建方法的有效性和精確性。

    1 前言

    汽車行駛工況(Driving Cycle)又稱車輛測試循環(huán),是描述典型車輛行駛的速度-時間曲線,用于確定車輛污染物排放量、燃油消耗量、新車型技術開發(fā)和評估以及測定交通控制的風險等,是汽車工業(yè)一項共性的核心技術[1]。目前,中國以及大多數(shù)國家是參照歐洲的汽車行駛工況[2]制定并實施本國的排放法規(guī)。中國城市汽車行駛工況采用的是ECE(歐洲城市工況)工況,其加減速過程簡單,與實際車輛行駛狀態(tài)差距較大,已不能滿足車輛開發(fā)測試的需要。針對此問題,我國學者已經(jīng)先后建立了符合本地區(qū)的汽車工況[1,3~5],而在構建行駛工況時,主要采用主成分分析法與聚類分析法相結合的研究方法,其中K-均值聚類(K-Means Cluster)分析法在構建城市汽車行駛工況中應用較多,但在對K-均值聚類分析結果篩選時,普遍通過對比各類對應的汽車行駛狀態(tài)時間分布比例,人為判斷聚類結果是否合理,對聚類結果的合理性存在影響。

    為此,本文以實時采集的乘用車行駛數(shù)據(jù)為研究數(shù)據(jù)源,基于K-均值聚類方法對汽車行駛數(shù)據(jù)進行分析,提出一種以Silhouette函數(shù)篩選聚類結果,并根據(jù)聚類結果構建汽車行駛工況的方法,減少了人為選擇K-均值聚類結果時存在的誤差。

    2 數(shù)據(jù)采集與前期處理

    2.1 數(shù)據(jù)采集

    數(shù)據(jù)采集在廈門市城市道路上進行,為使采集的數(shù)據(jù)真實反映城市道路汽車的行駛特點,利用5輛私家車按照各自目的緊跟行駛車流的方式采集數(shù)據(jù),時間為早上7:00至晚上20∶00,共15天。車輛行駛路線如圖1所示,涵蓋海滄區(qū)的霞飛路、馬青路、陽光路、海滄大橋以及湖里區(qū)湖里大道和仙岳路,體現(xiàn)較為典型的廈門島內(湖里區(qū))與島外(海滄區(qū))的道路交通情況。

    圖1 試驗車輛行駛路線

    數(shù)據(jù)采集車載終端與車輛的OBD接口相連接,如圖2所示。在汽車行駛過程中,以1 Hz的頻率采集包括車速、發(fā)動機轉速、扭矩百分比、瞬時油耗、進氣歧管溫度、發(fā)動機負荷、進氣歧管壓力等參數(shù)。通過信息化數(shù)據(jù)平臺對車輛進行實時監(jiān)測,確保采集數(shù)據(jù)的連續(xù)性和正確性。

    圖2 數(shù)據(jù)采集車載終端及安裝效果

    2.2 數(shù)據(jù)前期處理

    采集數(shù)據(jù)中存在汽車振動導致GPS漂移引起的部分加速度大于2.5 m·s-2的數(shù)據(jù)點,為減少后續(xù)數(shù)據(jù)分析的誤差,前期處理過程中需要將這些數(shù)據(jù)異常點剔出;為保證工況合成結果具有代表性,選用105s左右的數(shù)據(jù)點進行汽車行駛工況的構建分析。經(jīng)前處理后數(shù)據(jù)源的SAFD(Speed Acceleration Frequency Distribution)分布如圖3所示,部分速度-時間曲線和相對應的加速度-時間曲線見圖4。

    圖3 數(shù)據(jù)源的SAFD分布圖

    圖4 部分速度-時間曲線和相對應的加速度-時間曲線

    3 數(shù)據(jù)分析

    3.1 運動學片段分析

    汽車在行駛過程中,由于受到道路交通條件的限制,存在多次怠速、加速、巡航和減速的狀態(tài)。為此將從1個怠速開始到下1個怠速開始定義為1個運動學片段。

    汽車行駛狀態(tài)定義[6]為:怠速為發(fā)動機正在工作但車速為0的連續(xù)過程;加速為汽車加速度大于0.1 m/s2的連續(xù)過程;巡航為汽車加速度的絕對值小于0.1 m/s2但車速不為0的連續(xù)過程;減速為汽車加速度小于-0.1 m/s2的連續(xù)過程。根據(jù)以上定義,可將汽車行駛過程看作如圖5所示的多個運動學片段的組合,通過對每個運動學片段的行駛特征參數(shù)進行分析,得出該車輛的行駛特征。

    在汽車行駛過程中,需要選用一些基本的特征參數(shù)反映每個運動學片段的行駛特征。本文選用如表1所列的15個特征參數(shù)描述運動學片段。將數(shù)據(jù)源分割成297個運動學片段并計算各運動學片段的15個特征參數(shù),部分分析結果見表2。

    圖5 運動學片段

    表1 用于描述運動學片段的15個特征參數(shù)

    表2 運動學片段特征參數(shù)值

    3.2 特征參數(shù)標準化

    基于每個運動學片段中汽車行駛速度的變化可計算相應運動學片段的15個特征參數(shù)。由于各特征參數(shù)的量綱不統(tǒng)一會引起各變量取值的分散程度差異較大,導致后續(xù)進行降維分析和聚類分析過程中優(yōu)先照顧方差較大的變量,使處理結果的穩(wěn)定性差,為此,在對特征參數(shù)進行降維處理前進行特征參數(shù)標準化:

    式中,aij(i=1,2,…;j=1,2,…,n)是第i個運動學片段的第j個參數(shù)。

    對矩陣(1)進行標準化得到標準化矩陣X:

    3.3 特征參數(shù)降維

    利用特征參數(shù)對運動學片段進行描述的過程中,為保證運動學片段的特征信息不丟失,選取的15個特征參數(shù)必然存在所反映的信息重疊,為此采用主成分分析法對特征參數(shù)進行降維處理。

    由矩陣(2)計算協(xié)方差矩陣Σ:

    由矩陣(1)得相關矩陣R:

    用λi表示矩陣R的特征參數(shù)λ1≥λ2≥…≥λn≥0,計算出相應的正交化特征參數(shù)向量為:

    通過對297個運動學片段的特征參數(shù)進行主成分分析,得到如表3所列的主成分貢獻率及累積貢獻率。

    表3 主成分貢獻率及累積貢獻率

    理論上,進行主成分分析時,選取累積貢獻率大于80%的主成分進行后續(xù)的聚類分析,雖然會有一定的信息損失,但對最終的結果影響較小[7]。由表3可知,前8個主成分累積貢獻率達到98.682%,已經(jīng)包含15個特征參數(shù)的所有信息,其中前3個主成分的累積貢獻率之和已經(jīng)達到86.404%,并且這3個主成分的特征值都大于1,因此,用前3個主成分就能反映出15個特征參數(shù)的大部分信息。

    某個參數(shù)在某個主成分上的載荷系數(shù)的絕對值越大,則說明該參數(shù)與這個主成分的相關系數(shù)越大。表4為主成分載荷矩陣,由表4可知,第1主成分M1主要反映減速時間比例、行駛距離、片段持續(xù)時間、加速時間比例、巡航時間比例、平均速度;第2主成分M2主要反映加速度標準偏差、平均加速度、平均減速度;第3主成分M3主要反映怠速時間比例。根據(jù)主成分分析結果,選取第1、第2和第3主成分進一步分析。

    3.4 K-均值聚類分析及結果選擇

    K-均值聚類分析法是根據(jù)對所研究問題的了解程度確定分類數(shù),這樣即可在每一類中選擇1個有代表性的樣品作為聚點(初始聚點),其它樣品根據(jù)與該點的“親疏程度”進行分類,這種“親疏程度”可以用歐氏距離(Euclidian Distance)、歐氏平方距離(Squared Euclidian Distance)、明可夫斯基距離(Minkowski Distance)、切比雪夫距離(Chebyshev Distance)和布洛克距離(Block Distance)5種片段差異程度表示。

    在K-均值聚類分析過程中,參考相關文獻[3,4]采用毆氏平方距離表示各類中的“親疏程度”。其定義如下:設xik表示第i個片段的第k個變量,每個片段定義了p個變量,dik表示不同片段之間的距離,則歐氏平方距離為:

    基于特征參數(shù)的第1、第2和第3主成分,運用K-均值聚類分析方法對所有的運動學片段進行分析,得到分為兩類、3類和4類的聚類分析結果。根據(jù)K-均值聚類分析結果,運用Silhouette函數(shù)繪制輪廓圖,從輪廓圖上判斷每個運動學片段的分類是否合理。其中Silhouette函數(shù)的定義[8]如下:

    式中,a為第i個運動學片段與同類運動學片段之間的平均距離;b為一個向量,其元素是第i個運動學片段與不同類的各運動學片段間的平均距離。

    式(13)中,L(i)的取值范圍為[-1,1],L(i)值越大,說明第i個運動學片段分類越合理;當L(i)<0時,說明第i個運動學片段的分類不合理,應該還有比目前更好的分類。

    K-均值聚類各分類結果的Silhouette函數(shù)值的輪廓見圖6。由圖6a可看出,當分兩類時,各類的Silhouette函數(shù)值均大于0,說明分類時各類型已經(jīng)被很好地區(qū)分;由圖6b可看出,分3類時,第3種類型的Silhouette函數(shù)值出現(xiàn)少量負值;由圖6c可看出,分4類時,第1種類型和第4種類型的Silhouette函數(shù)值出現(xiàn)少量負值的情況,說明聚類分析分為3類和4類時存在未被很好區(qū)分的運動學片段。根據(jù)分析結果,選取分為兩類的結果作為K-均值聚類分析的最終結果。

    圖6 不同分類時Silhouette函數(shù)值的輪廓

    根據(jù)Silhouette函數(shù)值的選擇結果,對分為兩類時各類型中行駛狀態(tài)所占的時間比例進行統(tǒng)計,結果如圖7所示。由圖7可看出,類型1中怠速時間比例所占比重最大,達到37%,其它行駛狀態(tài)的時間比例在20%左右,反映了汽車在城市擁擠道路的行駛過程;類型2中怠速時間比例所占比重不到5%,加減速時間比例和巡航時間比例分布較均勻,反映了汽車在較為順暢道路行駛的過程。

    圖7 各類型行駛狀態(tài)時間比例

    根據(jù)廈門市實際道路交通情況可知,廈門島內湖里區(qū)交通狀況復雜,道路擁擠,因此,第1種類型主要反映試驗車輛在廈門島內的行駛狀態(tài);廈門島外海滄區(qū)道路較為順暢,因此,第2種類型主要反映試驗車輛在廈門海滄區(qū)的行駛過程。

    4 行駛工況的構建

    4.1 工況合成方法

    基于K-均值聚類分析結果,已知各分類中所包含的運動學片段與其聚類中心的大小,由式(12)可知,當dij越小時,說明該運動學片段越能反映該類中的行駛特點。因此,本文的工況合成方法為:在每一種類型中,根據(jù)聚類中心的大小,由小到大分別篩選10個候選運動學片段;當某個類型中所包含的運動學片段不足10個時,將整個類別中的運動學片段作為候選的運動學片段;以各類的候選運動學片段為基礎,從中隨機篩選運動學片段組合成≥1 200 s的候選汽車行駛工況。

    4.2 代表行駛工況的確定

    為了從大量的候選工況中選取具有代表性的汽車行駛工況,國外學者采用最小性能值PV(Performance Value)作為篩選代表性工況的標準[9]。本文提出將加速時間比例、巡航時間比例、怠速時間比例、平均速度、平均加速度、加速度標準偏差以及最大SAFD差異值等7項作為主要篩選參數(shù)。通過仿真模擬計算,確定前6項與數(shù)據(jù)源的相對誤差控制在±5%范圍內,同時最大SAFD差異值的誤差控制在±5%作為代表工況的篩選標準。

    4.3 工況合成結果分析

    根據(jù)前述工況合成方法和篩選標準,運用MATLAB軟件合成如圖8所示的1 514 s的汽車代表行駛工況,其SAFD分布如圖9所示。

    圖8 代表行駛工況的速度-時間曲線

    圖9 代表行駛工況的SAFD分布圖

    將構建的代表行駛工況與實際采集數(shù)據(jù)源的特征參數(shù)進行對比分析,結果如表5所列,表5中主要特征參數(shù)的相對誤差都維持在±5%的范圍內。進一步分析代表行駛工況與數(shù)據(jù)源的SAFD差異值,如圖10所示,由圖10可看出,代表工況和數(shù)據(jù)源對應點的SAFD差異值在±2%的范圍內。通過以上分析表明,所構建的代表工況可以反映試驗車輛的整體行駛特征。

    表5 特征參數(shù)差異率對比

    圖10 數(shù)據(jù)源與代表工況的SAFD差異值

    依據(jù)工況構建過程中的篩選準則,將合成的廈門市乘用車行駛工況與國外主要標準循環(huán)工況進行對比,結果如表6所列。由表6可知,合成的廈門市乘用車行駛工況中,加速時間比例和巡航時間比例與美國城市工況較接近,而平均加速度與日本汽車測試工況(J10-15)和歐盟的NEDC循環(huán)工況最為接近。但廈門工況中怠速時間比例低于美國城市工況,與歐盟提出的輕型車測試循環(huán)工況(WLTC)相比,各參數(shù)都存在較大差距。

    表6 廈門市合成行駛工況與國外主要標準工況的比較

    5 結束語

    基于廈門市乘用車典型交通線路實時采集的行駛數(shù)據(jù),分別運用運動學片段分析法、主成分分析法和K-均值聚類分析法對實測數(shù)據(jù)降維和分類。提出以Silhouette函數(shù)實現(xiàn)對聚類結果的篩選,減少人為選擇的誤差,并根據(jù)聚類中心的大小篩選所需運動學片段構建候選工況。在目標代表工況的遴選標準方面,提出綜合控制6個特征參數(shù)和SAFD最大差異值誤差的方法,并通過實例驗證了所提工況構建方法的有效性和精確性。通過將所構建的廈門工況與國外常用標準循環(huán)工況進行對比表明,所構建的廈門市工況與UDDS工況較為接近,但是怠速時間比例低于UDDS工況。

    1 石琴,鄭與波,姜平.基于運動學片段的城市道路行駛工況的研究.汽車工程,2011,33(3)∶256~261.

    2 張建偉,李孟良,艾國和,等.車輛行駛工況與特征的研究.汽車工程,2005,27(2)∶220~224,245.

    3 蔡鍔,李陽陽,李春明,等.基于K-均值聚類算法的西安市汽車行駛工況合成技術研究.汽車技術,2015(8)∶33~36.

    4 胡志遠,秦艷,譚丕強,等.基于大樣本的上海市乘用車行駛工況構建.同濟大學學報(自然科學版),2015,43(10)∶1523~1527.

    5 李孟良,張建偉,張富,等.中國城市乘用車實際行駛工況的研究.汽車工程,2006,28(6)∶554~557,529.

    6 楊延相,蔡曉林,杜青,等.天津市道路汽車行駛工況的研究.汽車工程,2002,24(3)∶200~204.

    7 范金城,梅長林.數(shù)據(jù)分析.北京:科學出版社,2010.

    8 張德豐.MATLAB概率與數(shù)理統(tǒng)計分析.北京:機械工業(yè)出版社,2010.

    9 Jie Lin,Niemeier D A.Estimating Regional Air Quality VehicleEmission Inventories∶Constructing Robust Driving Cycles.SAE.2002.

    Research on the Construction Method of Driving Cycle for the City Car Based on K-Means Cluster Analysis

    Peng Yuhui1,2,Yang Huibao1,Li Mengliang2,Qiao Xueqi3
    (1.Fuzhou University,Fuzhou,350116;2.China Automotive Technology and Research Center,Tianjin 300300;3.Xiamen Golden Dragon Wagon Bus Co.,Ltd,361006)

    The construction method of city road driving cycle was studied with the real-time collected passenger vehicle driving data as data source.The kinematics fragment analysis,principal component analysis and K-means clustering analysis were applied separately for dimensionality reduction of the measured parameters and classification,then Silhouette Function was proposed to screen the clustering results,to reduce artificial selection error,and the required kinematic fragment candidate cycle was selected according to size of the clustering center.In term of choice of typical driving cycle,an evaluation criterion that integrated six characteristic parameters and the maximum SAFD difference value was proposed.Finally test proved validity and accuracy of the construction method of the driving cycle.

    Driving cycle,Kinematic fragment,Principal component analysis,K-means clustering analysis

    汽車行駛工況 運動學片段 主成分分析 K-均值聚類分析

    U467.1+1

    A

    1000-3703(2017)11-0013-06

    工業(yè)和信息化部“中國新能源汽車產(chǎn)品檢測工況研究開發(fā)”項目資助(FZU201600201603)。

    (責任編輯文 楫)

    修改稿收到日期為2017年3月5日。

    猜你喜歡
    數(shù)據(jù)源特征參數(shù)運動學
    故障診斷中信號特征參數(shù)擇取方法
    基于特征參數(shù)化的木工CAD/CAM系統(tǒng)
    基于MATLAB的6R機器人逆運動學求解分析
    基于D-H法的5-DOF串并聯(lián)機床運動學分析
    Web 大數(shù)據(jù)系統(tǒng)數(shù)據(jù)源選擇*
    基于PSO-VMD的齒輪特征參數(shù)提取方法研究
    基于不同網(wǎng)絡數(shù)據(jù)源的期刊評價研究
    基于運動學原理的LBI解模糊算法
    基于真值發(fā)現(xiàn)的沖突數(shù)據(jù)源質量評價算法
    統(tǒng)計特征參數(shù)及多分類SVM的局部放電類型識別
    電測與儀表(2015年7期)2015-04-09 11:40:04
    亚洲少妇的诱惑av| 精品久久久精品久久久| 狂野欧美激情性xxxx| 后天国语完整版免费观看| 交换朋友夫妻互换小说| 91麻豆精品激情在线观看国产 | 三上悠亚av全集在线观看| 国产免费一区二区三区四区乱码| 老司机午夜十八禁免费视频| 久久久国产精品麻豆| 久久毛片免费看一区二区三区| 久久精品国产a三级三级三级| 亚洲成人手机| 在线观看国产h片| av片东京热男人的天堂| 麻豆av在线久日| 黑丝袜美女国产一区| 巨乳人妻的诱惑在线观看| 亚洲精品第二区| av天堂久久9| 最黄视频免费看| 亚洲一码二码三码区别大吗| 亚洲成av片中文字幕在线观看| 91成人精品电影| 婷婷色av中文字幕| 日本av手机在线免费观看| cao死你这个sao货| 午夜两性在线视频| 久久精品亚洲av国产电影网| 久热这里只有精品99| 一本大道久久a久久精品| 女人精品久久久久毛片| 九草在线视频观看| 欧美日韩黄片免| 老鸭窝网址在线观看| 国产野战对白在线观看| 国产激情久久老熟女| 香蕉丝袜av| 欧美性长视频在线观看| 免费看十八禁软件| 国产福利在线免费观看视频| 高清黄色对白视频在线免费看| 精品久久久久久久毛片微露脸 | 99久久综合免费| 免费一级毛片在线播放高清视频 | 99九九在线精品视频| 亚洲中文日韩欧美视频| 在线观看免费高清a一片| 一级毛片女人18水好多 | 亚洲国产精品国产精品| 纯流量卡能插随身wifi吗| 久久久久久人人人人人| 丝袜在线中文字幕| 三上悠亚av全集在线观看| 999精品在线视频| 999精品在线视频| 激情视频va一区二区三区| 一边摸一边做爽爽视频免费| 9热在线视频观看99| 亚洲av电影在线观看一区二区三区| 少妇人妻 视频| 国产国语露脸激情在线看| 99热国产这里只有精品6| 一区二区av电影网| 国产精品一国产av| 久久性视频一级片| 久9热在线精品视频| 丰满迷人的少妇在线观看| 国产伦人伦偷精品视频| 19禁男女啪啪无遮挡网站| 国产又色又爽无遮挡免| 国产野战对白在线观看| 国产又爽黄色视频| av在线app专区| 巨乳人妻的诱惑在线观看| 国产成人系列免费观看| 一边摸一边抽搐一进一出视频| 一区二区三区四区激情视频| 校园人妻丝袜中文字幕| 亚洲人成电影免费在线| 国产女主播在线喷水免费视频网站| 日本91视频免费播放| 新久久久久国产一级毛片| 午夜免费成人在线视频| 狠狠精品人妻久久久久久综合| 纯流量卡能插随身wifi吗| 又紧又爽又黄一区二区| 日日爽夜夜爽网站| 看免费av毛片| 美女福利国产在线| 女人爽到高潮嗷嗷叫在线视频| 人妻人人澡人人爽人人| 久久久久网色| 狂野欧美激情性bbbbbb| 99久久99久久久精品蜜桃| 国产高清videossex| 国产成人免费观看mmmm| 操出白浆在线播放| 亚洲精品自拍成人| 国产精品av久久久久免费| www.av在线官网国产| 久热这里只有精品99| 午夜免费成人在线视频| 嫩草影视91久久| 午夜两性在线视频| 亚洲国产欧美日韩在线播放| 成年av动漫网址| 精品人妻一区二区三区麻豆| 久久狼人影院| 蜜桃在线观看..| 人妻一区二区av| 久热这里只有精品99| 十八禁高潮呻吟视频| 丝袜脚勾引网站| 在现免费观看毛片| 一级毛片我不卡| 免费在线观看日本一区| 久久久国产一区二区| 老司机午夜十八禁免费视频| 成年人黄色毛片网站| 亚洲精品久久久久久婷婷小说| 国产女主播在线喷水免费视频网站| 老司机影院毛片| 亚洲九九香蕉| 人成视频在线观看免费观看| 在线观看免费午夜福利视频| 亚洲人成电影观看| 青草久久国产| 久久久久久亚洲精品国产蜜桃av| 欧美日韩一级在线毛片| 两性夫妻黄色片| 操出白浆在线播放| 最新在线观看一区二区三区 | 色综合欧美亚洲国产小说| 国产激情久久老熟女| 一级毛片女人18水好多 | 熟女av电影| 国产男女超爽视频在线观看| 亚洲欧美一区二区三区黑人| 国产精品国产三级专区第一集| 69精品国产乱码久久久| 国产一区二区在线观看av| 美女扒开内裤让男人捅视频| 国产1区2区3区精品| 在线观看免费高清a一片| 亚洲精品国产av成人精品| 亚洲午夜精品一区,二区,三区| 超碰97精品在线观看| 嫩草影视91久久| 女人被躁到高潮嗷嗷叫费观| 久久精品国产a三级三级三级| 天堂俺去俺来也www色官网| 亚洲av综合色区一区| 久久久久久久大尺度免费视频| 一二三四在线观看免费中文在| 飞空精品影院首页| 精品亚洲成a人片在线观看| 午夜影院在线不卡| 97在线人人人人妻| 国产免费又黄又爽又色| 美女午夜性视频免费| 日韩视频在线欧美| 亚洲精品一二三| 亚洲成色77777| av网站免费在线观看视频| 汤姆久久久久久久影院中文字幕| 欧美人与性动交α欧美精品济南到| 日韩中文字幕视频在线看片| 国产亚洲午夜精品一区二区久久| 18禁裸乳无遮挡动漫免费视频| 狂野欧美激情性xxxx| 婷婷丁香在线五月| 国产片内射在线| 丰满迷人的少妇在线观看| 啦啦啦中文免费视频观看日本| netflix在线观看网站| 9热在线视频观看99| 久久久久久免费高清国产稀缺| 91字幕亚洲| 在线观看人妻少妇| 欧美黄色片欧美黄色片| 19禁男女啪啪无遮挡网站| 男女免费视频国产| 亚洲熟女精品中文字幕| 少妇人妻 视频| 国产激情久久老熟女| 好男人视频免费观看在线| 97人妻天天添夜夜摸| av在线app专区| 欧美97在线视频| 国产精品麻豆人妻色哟哟久久| 精品人妻1区二区| 天天躁狠狠躁夜夜躁狠狠躁| 涩涩av久久男人的天堂| 国产亚洲av片在线观看秒播厂| 极品少妇高潮喷水抽搐| 国产在线观看jvid| 99国产精品免费福利视频| 老司机影院毛片| 性色av乱码一区二区三区2| 永久免费av网站大全| 日韩大码丰满熟妇| 一区二区av电影网| 亚洲精品国产一区二区精华液| 久久99热这里只频精品6学生| 亚洲国产最新在线播放| 久久人人爽人人片av| 97在线人人人人妻| 人妻一区二区av| 肉色欧美久久久久久久蜜桃| 香蕉国产在线看| av又黄又爽大尺度在线免费看| 成人18禁高潮啪啪吃奶动态图| netflix在线观看网站| 一级黄色大片毛片| 精品亚洲乱码少妇综合久久| 国产真人三级小视频在线观看| 国产不卡av网站在线观看| 久久久久精品人妻al黑| 国产成人精品久久久久久| bbb黄色大片| 亚洲国产av影院在线观看| 男人操女人黄网站| 日本猛色少妇xxxxx猛交久久| 午夜免费鲁丝| 亚洲欧美一区二区三区国产| 国产精品麻豆人妻色哟哟久久| 亚洲av日韩精品久久久久久密 | videosex国产| av欧美777| 人人妻人人澡人人爽人人夜夜| 中文字幕人妻熟女乱码| √禁漫天堂资源中文www| 亚洲精品av麻豆狂野| 另类精品久久| 捣出白浆h1v1| 欧美日韩福利视频一区二区| 最新在线观看一区二区三区 | 一二三四在线观看免费中文在| 久久久国产欧美日韩av| 国产欧美亚洲国产| 男女边摸边吃奶| 最新在线观看一区二区三区 | 国产精品一二三区在线看| 国产xxxxx性猛交| 精品久久久久久久毛片微露脸 | 国产午夜精品一二区理论片| 久久国产精品大桥未久av| 男人操女人黄网站| 亚洲国产中文字幕在线视频| 青青草视频在线视频观看| 一区二区三区激情视频| 国产又色又爽无遮挡免| 男女边摸边吃奶| 七月丁香在线播放| 真人做人爱边吃奶动态| 日韩一区二区三区影片| 日韩电影二区| 国产熟女欧美一区二区| 亚洲av成人精品一二三区| 91精品国产国语对白视频| 丰满饥渴人妻一区二区三| 一区二区日韩欧美中文字幕| 亚洲七黄色美女视频| 十分钟在线观看高清视频www| 亚洲成人手机| 亚洲中文字幕日韩| a级片在线免费高清观看视频| 亚洲欧美一区二区三区黑人| 七月丁香在线播放| netflix在线观看网站| 国产亚洲一区二区精品| 国产欧美日韩一区二区三 | 亚洲精品在线美女| 国产精品久久久av美女十八| 免费观看人在逋| 精品熟女少妇八av免费久了| 波多野结衣一区麻豆| av又黄又爽大尺度在线免费看| 日日夜夜操网爽| 色婷婷久久久亚洲欧美| 下体分泌物呈黄色| 丁香六月天网| 日韩欧美一区视频在线观看| 最近最新中文字幕大全免费视频 | 侵犯人妻中文字幕一二三四区| av在线老鸭窝| 交换朋友夫妻互换小说| 久久久精品94久久精品| 每晚都被弄得嗷嗷叫到高潮| 国产福利在线免费观看视频| 激情五月婷婷亚洲| 三上悠亚av全集在线观看| 国产欧美亚洲国产| 少妇人妻久久综合中文| 黄色毛片三级朝国网站| www.自偷自拍.com| 日韩,欧美,国产一区二区三区| 超色免费av| 成人三级做爰电影| 精品久久久精品久久久| 人人妻,人人澡人人爽秒播 | 国产精品成人在线| 欧美另类一区| 在现免费观看毛片| 日本91视频免费播放| 精品国产乱码久久久久久男人| 黄色片一级片一级黄色片| 中文字幕av电影在线播放| 天堂俺去俺来也www色官网| 日韩视频在线欧美| 男女床上黄色一级片免费看| 如日韩欧美国产精品一区二区三区| 国产激情久久老熟女| 老司机在亚洲福利影院| 午夜免费观看性视频| 久久天躁狠狠躁夜夜2o2o | 亚洲专区中文字幕在线| 大香蕉久久网| 亚洲精品久久久久久婷婷小说| 欧美少妇被猛烈插入视频| 大片免费播放器 马上看| 天天躁夜夜躁狠狠躁躁| 19禁男女啪啪无遮挡网站| 国产一区二区三区av在线| 中国国产av一级| 国产亚洲一区二区精品| 亚洲成人国产一区在线观看 | netflix在线观看网站| 伊人久久大香线蕉亚洲五| 欧美黑人欧美精品刺激| 国产色视频综合| 91老司机精品| 久久久久久免费高清国产稀缺| 久久人人97超碰香蕉20202| 国产一卡二卡三卡精品| 天天操日日干夜夜撸| av片东京热男人的天堂| 啦啦啦中文免费视频观看日本| 晚上一个人看的免费电影| 国产免费现黄频在线看| 18禁黄网站禁片午夜丰满| 久久 成人 亚洲| 精品人妻1区二区| 狂野欧美激情性xxxx| av又黄又爽大尺度在线免费看| 99国产精品一区二区三区| 新久久久久国产一级毛片| 少妇裸体淫交视频免费看高清 | 国产免费现黄频在线看| 免费观看a级毛片全部| 黄色怎么调成土黄色| 侵犯人妻中文字幕一二三四区| 亚洲久久久国产精品| 天天躁夜夜躁狠狠久久av| 9色porny在线观看| 国产野战对白在线观看| 日韩av在线免费看完整版不卡| 多毛熟女@视频| 亚洲欧美激情在线| 国产精品二区激情视频| 美女主播在线视频| 亚洲欧美一区二区三区久久| 午夜精品国产一区二区电影| 色网站视频免费| 爱豆传媒免费全集在线观看| a级毛片黄视频| 80岁老熟妇乱子伦牲交| 国产精品国产三级专区第一集| 久久影院123| 亚洲三区欧美一区| 岛国毛片在线播放| 精品熟女少妇八av免费久了| 国产成人精品久久二区二区免费| 国产成人系列免费观看| 国产伦人伦偷精品视频| 欧美国产精品一级二级三级| 九色亚洲精品在线播放| 国产伦理片在线播放av一区| 别揉我奶头~嗯~啊~动态视频 | 午夜福利影视在线免费观看| 精品一品国产午夜福利视频| 亚洲人成电影免费在线| 777久久人妻少妇嫩草av网站| 欧美激情高清一区二区三区| 精品高清国产在线一区| 日本黄色日本黄色录像| 中文字幕最新亚洲高清| 亚洲国产毛片av蜜桃av| 久久九九热精品免费| 欧美日韩精品网址| 国产亚洲欧美在线一区二区| 男女高潮啪啪啪动态图| 狠狠精品人妻久久久久久综合| 人体艺术视频欧美日本| 国产激情久久老熟女| 精品亚洲成a人片在线观看| 精品熟女少妇八av免费久了| 一级黄片播放器| 欧美黑人欧美精品刺激| 一级毛片电影观看| 十八禁高潮呻吟视频| 亚洲欧美一区二区三区黑人| 久久 成人 亚洲| 尾随美女入室| 爱豆传媒免费全集在线观看| 成人国产一区最新在线观看 | 国产精品久久久久久精品电影小说| 免费在线观看完整版高清| 免费久久久久久久精品成人欧美视频| 亚洲国产成人一精品久久久| 色精品久久人妻99蜜桃| 一级毛片女人18水好多 | 免费久久久久久久精品成人欧美视频| 国产精品久久久久久精品古装| 中文字幕人妻丝袜制服| 亚洲欧美日韩高清在线视频 | 午夜免费鲁丝| 亚洲人成网站在线观看播放| 亚洲av成人精品一二三区| 999精品在线视频| 99精国产麻豆久久婷婷| 日本av手机在线免费观看| 亚洲精品久久成人aⅴ小说| bbb黄色大片| 亚洲精品国产区一区二| 男女边摸边吃奶| 一级毛片我不卡| 国产欧美日韩精品亚洲av| 自拍欧美九色日韩亚洲蝌蚪91| 狂野欧美激情性xxxx| 欧美+亚洲+日韩+国产| 精品国产乱码久久久久久小说| 久热爱精品视频在线9| 亚洲人成77777在线视频| 人人澡人人妻人| 国产日韩欧美亚洲二区| 国产精品一国产av| 欧美精品亚洲一区二区| 色网站视频免费| 久久99热这里只频精品6学生| 女性被躁到高潮视频| 精品亚洲成国产av| 国产女主播在线喷水免费视频网站| 在线观看国产h片| 欧美在线黄色| 各种免费的搞黄视频| 好男人电影高清在线观看| 国产精品偷伦视频观看了| a级毛片黄视频| 欧美日韩亚洲综合一区二区三区_| 亚洲精品久久午夜乱码| 国产成人一区二区在线| 久久影院123| 久热这里只有精品99| 91麻豆av在线| 亚洲伊人色综图| 如日韩欧美国产精品一区二区三区| 亚洲欧美精品综合一区二区三区| 中文字幕人妻丝袜制服| 一本综合久久免费| 国产有黄有色有爽视频| 久久精品国产综合久久久| 国产精品九九99| 欧美日韩视频高清一区二区三区二| 亚洲视频免费观看视频| 久久久精品94久久精品| 国产1区2区3区精品| 久久毛片免费看一区二区三区| 日韩大片免费观看网站| 午夜精品国产一区二区电影| 麻豆国产av国片精品| 一级黄片播放器| 久久天躁狠狠躁夜夜2o2o | 天天影视国产精品| e午夜精品久久久久久久| 老熟女久久久| 人人妻人人澡人人看| 国产熟女欧美一区二区| 午夜福利视频精品| 在现免费观看毛片| 你懂的网址亚洲精品在线观看| 老司机亚洲免费影院| 男女边摸边吃奶| 欧美日韩亚洲高清精品| 91麻豆av在线| 欧美黑人欧美精品刺激| 人妻人人澡人人爽人人| 丝袜人妻中文字幕| 乱人伦中国视频| 午夜激情久久久久久久| 在线观看免费视频网站a站| 黄色怎么调成土黄色| 亚洲九九香蕉| 国产日韩欧美视频二区| 性少妇av在线| 久久毛片免费看一区二区三区| 中文字幕av电影在线播放| 肉色欧美久久久久久久蜜桃| 日韩中文字幕欧美一区二区 | 在线观看免费视频网站a站| 色视频在线一区二区三区| 人体艺术视频欧美日本| 麻豆乱淫一区二区| 色婷婷av一区二区三区视频| 午夜免费男女啪啪视频观看| 精品高清国产在线一区| 日韩人妻精品一区2区三区| 中文字幕亚洲精品专区| av网站在线播放免费| 精品少妇黑人巨大在线播放| 老司机深夜福利视频在线观看 | 黄片播放在线免费| 中文字幕另类日韩欧美亚洲嫩草| 亚洲中文字幕日韩| e午夜精品久久久久久久| 午夜免费男女啪啪视频观看| 精品卡一卡二卡四卡免费| 日韩人妻精品一区2区三区| 午夜精品国产一区二区电影| 国产男女内射视频| 女性生殖器流出的白浆| 国产精品99久久99久久久不卡| 最近手机中文字幕大全| 日韩大码丰满熟妇| 日本五十路高清| 久久久久精品国产欧美久久久 | 国产精品一区二区在线观看99| 精品国产国语对白av| 国产有黄有色有爽视频| 中文乱码字字幕精品一区二区三区| 人成视频在线观看免费观看| 侵犯人妻中文字幕一二三四区| 十八禁人妻一区二区| 热re99久久精品国产66热6| 两性夫妻黄色片| 少妇粗大呻吟视频| 日韩制服骚丝袜av| 波多野结衣av一区二区av| 亚洲国产欧美在线一区| 亚洲第一av免费看| 欧美国产精品一级二级三级| 免费在线观看完整版高清| 精品少妇久久久久久888优播| 国产三级黄色录像| 精品久久久精品久久久| 亚洲一码二码三码区别大吗| 少妇 在线观看| 国产精品偷伦视频观看了| 丝袜在线中文字幕| 国产一区二区在线观看av| 国产亚洲av片在线观看秒播厂| 色综合欧美亚洲国产小说| av视频免费观看在线观看| 欧美久久黑人一区二区| 日韩,欧美,国产一区二区三区| 91精品伊人久久大香线蕉| 色网站视频免费| 不卡av一区二区三区| 交换朋友夫妻互换小说| 午夜两性在线视频| 别揉我奶头~嗯~啊~动态视频 | 日本欧美视频一区| 午夜久久久在线观看| 宅男免费午夜| 亚洲精品久久午夜乱码| 人人妻人人澡人人爽人人夜夜| 亚洲精品国产av成人精品| 欧美精品一区二区大全| 侵犯人妻中文字幕一二三四区| 肉色欧美久久久久久久蜜桃| 最新的欧美精品一区二区| 国产无遮挡羞羞视频在线观看| 在线观看免费午夜福利视频| av天堂在线播放| 亚洲国产欧美在线一区| 精品人妻熟女毛片av久久网站| 2021少妇久久久久久久久久久| 欧美精品一区二区免费开放| 日韩视频在线欧美| 真人做人爱边吃奶动态| 天天躁夜夜躁狠狠久久av| av天堂在线播放| 欧美久久黑人一区二区| 久久精品国产亚洲av高清一级| 在线亚洲精品国产二区图片欧美| svipshipincom国产片| 国产成人精品在线电影| 国产精品一二三区在线看| av又黄又爽大尺度在线免费看| 欧美激情极品国产一区二区三区| 午夜福利在线免费观看网站| 亚洲欧美激情在线| av又黄又爽大尺度在线免费看| 欧美激情极品国产一区二区三区| 亚洲欧美精品综合一区二区三区| 亚洲av男天堂| 十八禁网站网址无遮挡| 国产精品99久久99久久久不卡| 少妇人妻久久综合中文| 中文字幕人妻丝袜一区二区| 成年人午夜在线观看视频| 久久这里只有精品19| 亚洲国产精品999| 中文乱码字字幕精品一区二区三区| 国产成人一区二区三区免费视频网站 | 伊人久久大香线蕉亚洲五| av天堂久久9| 男女之事视频高清在线观看 | 久久精品成人免费网站| 国产野战对白在线观看| 一级毛片女人18水好多 | 亚洲国产欧美日韩在线播放| 亚洲国产中文字幕在线视频| 国产av精品麻豆|