999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于自然駕駛數據的危險評估算法研究?

2017-09-15 05:04:32朱西產沈劍平
汽車工程 2017年8期
關鍵詞:駕駛員

吳 斌,朱西產,沈劍平,李 霖

基于自然駕駛數據的危險評估算法研究?

吳 斌1,2,朱西產1,沈劍平2,李 霖3

(1.同濟大學汽車學院,上海 201804;2.國家機動車產品質量監督檢驗中心(上海),上海 201805;3.上海國際汽車城(集團)有限公司,上海 201805)

本文中提出了一種基于自然駕駛數據的汽車碰撞危險評估算法。首先分析了在自然駕駛工況下的駕駛員制動和轉向的輸入特性,研究了表征駕駛員的制動避撞和轉向避撞緊急程度的行為特征參數,建立了緊急制動和緊急轉向的駕駛員數學模型;在此基礎上,分析了不同車速下制動避撞和轉向避撞的特征,提出了基于自然駕駛數據的危險評估算法。仿真結果表明,本文中提出的危險評估算法可同時兼顧制動避撞和轉向避撞,并與駕駛員在真實交通環境下的緊急避撞行為相對應,可用于汽車避撞系統控制策略的開發。

自然駕駛數據;緊急避撞;駕駛員行為;危險評估

前言

汽車碰撞危險估計是指通過分析當前駕駛環境,估計本車及周圍駕駛環境的未來變化,判斷車輛的潛在碰撞危險程度。危險估計算法是智能汽車安全系統制定控制策略的基礎,決定了系統發出預警或主動干預的時刻及可靠性。

目前主要有兩類危險估計算法:確定性算法和概率性算法。其中,概率性算法通過計算本車與不同目標發生碰撞的概率分布來估算危險程度,該方法用于復雜的交通環境,能夠考慮多個目標之間的相互作用,但計算量巨大,目前還處在研究中[1]。確定性危險估計方法具有計算效率高、便于開發控制策略等優點,從而得到廣泛應用。文獻[2]中對于確定性危險估計算法的表征參數進行了總結,并將這些參數分成基于時間、基于距離和基于減速度3類。文獻[3]中對于時間等各種估計測量參數進行了分析和比較。文獻[4]中通過實際交通工況數據對基于距離的危險估計算法與TTC(time-to-collision,兩車相對距離與相對速度計算得到的碰撞時間)進行了評價和比較。應用比較廣泛的危險估算方法,主要集中在縱向碰撞預警策略,如Mazda模型、Honda模型、伯克利模型和Jaguar模型等。但是,目前這些模型未能充分考慮我國駕駛員的行為特征,且無法同時兼容考慮轉向避撞的可行性。

本文中提出了一種基于我國自然駕駛工況的汽車碰撞危險估計算法。分析了在自然駕駛工況下駕駛員的制動和轉向避撞行為特征,建立了駕駛員緊急制動和轉向的數學模型;在此基礎上,提出了基于我國自然駕駛工況的危險估計算法,對汽車碰撞的危險區域進行分級劃分。

1 自然駕駛工況的數據采集

為了研究駕駛員在真實交通環境中的駕駛行為,本文中所使用的數據全部來自中國大型實車路試(China field operational test,China FOT)的先行實驗China-Pilot FOT。所有實驗車都統一安裝了數據采集設備,包括4個攝像頭,全面記錄駕駛過程中駕駛員操作信息和道路環境信息,如圖1所示。在為期2個月的數據采集中,收集了來自不同駕駛員的160個駕駛片段,得到5 600min的自然駕駛數據。

圖1 自然駕駛工況的采集

2 駕駛行為輸入特性的研究

對于自動駕駛輔助系統,核心問題是確定合適的介入時刻以達到安全效用和用戶接受度的最優化。傳統的危險估計算法中大多以車輛物理極限來定義碰撞無法避免時刻,但這種方法并不是最為合理。故本文中引入駕駛員行為數據,通過制動和轉向參數來定義工況的避撞操作的緊急程度,作為危險估計算法的輸入。

2.1 駕駛員的制動輸入特性

在研究駕駛員的制動輸入特性時,本文中利用制動減速度和制動減速度梯度來表征制動工況操作緊急程度。

圖2為自然駕駛工況下車輛的縱向加速度隨速度變化的分布圖,負值為制動減速度,表征了制動的強烈程度,代表了駕駛員在危險工況下對車輛制動力的利用水平。通過統計分析,自然駕駛工況中95%的駕駛員利用的最大制動減速度絕對值小于5m/s2;99%的駕駛員利用的最大制動減速度絕對值小于6m/s2;駕駛員最大的制動減速度為8m/s2,與地面的最大附著力相關。

圖2 縱向加速度隨速度變化的分布圖

縱向制動減速度梯度是指制動時減速度值的上升快慢,反映駕駛員踩制動踏板的速度,也是駕駛員對碰撞危險感知的直觀反應。圖3為自然駕駛工況下車輛縱向加速度梯度隨速度變化的分布圖,負值為制動減速度梯度。通過統計,在自然駕駛工況中95%可操作的最大制動減速度梯度小于10m/s3;99%的可操作的最大制動減速度梯度小于15m/s3;駕駛員可操作的最大制動減速度梯度為20m/s3。

利用制動減速度和制動減速度梯度來表達制動操作的緊急程度,根據上述分析將制動操作緊急程度分成3級,如表1所示。

2.2 駕駛員的轉向輸入特性

在研究駕駛員的轉向輸入特性時,側向加速度、側向加速度梯度、轉向盤轉角輸入、轉向盤轉速等參數都與駕駛員的轉向行為密切相關。其中轉向盤轉角輸入和轉向盤轉速與駕駛員轉向輸入線性相關,故選取轉向盤轉速作為駕駛員轉向操作劇烈程度的參數;同時考慮轉向過程中,車輛對地面附著力利用程度,故選取側向加速度作為轉向強烈程度的參數。

圖4為駕駛員的側向加速度隨速度變化的分布圖,圖中緊急程度1、緊急程度2和緊急程度3分別表示95%,99%范圍線和最大側向加速度線。

轉向側向加速度隨速度變化比較明顯:車速在40~80 km/h之間時,駕駛員使用的側向加速度較大;車速小于40km/h和大于80km/h時,駕駛員使用的側向加速度明顯降低。

圖5為駕駛員轉向盤角速度隨車速的變化分布圖,圖中緊急程度1、緊急程度2和緊急程度3分別表示95%,99%范圍線和最大轉向盤角速度線。隨著車速的增加,轉向盤角速度逐漸減小。在較低車速行駛時,駕駛員傾向于大幅度轉向盤轉角和角速度進行轉向避撞;在較高車速行駛時,駕駛員傾向于小幅度轉向盤轉角和角速度進行轉向避撞。

利用側向加速度和轉向盤角速度來表達轉向工況的操作緊急程度,將轉向操作緊急程度分成3級,如表2所示。

3 駕駛員緊急制動和轉向模型建立

在駕駛過程中車輛受駕駛員控制,在某些工況下車輛運動雖然從車輛動力學上是可行的,但駕駛員通過操作實現這種車輛運動的可能性很低。本節利用自然駕駛工況數據,分析駕駛員緊急工況下的制動轉向行為特征,并建立相應的駕駛員數學模型,作為危險估算方法的輸入,可提高算法的安全效用。

表2 轉向操作緊急程度的劃分

3.1 駕駛員緊急制動模型的建立

以表1中的緊急程度1最大制動減速度絕對值>5m/s2、最大制動減速度梯度絕對值>10m/s3為篩選閾值,并結合視頻錄像篩選出制動危險工況400多例。緊急制動工況下駕駛員遵循的原則是保證駕駛的平順性,即緊急制動是不斷調整優化控制的過程,符合人體行為學規律,優化目標為

式中j(t)為制動減速度梯度。利用五次多項式來表示制動過程中的車輛縱向位移變化S(t),并推導出其他運動參數:

假設車輛初始狀態:S(0)=0,v(0)=v0,a(0)=a0;車輛的最終狀態:S(T)=Sf,v(T)=vf,a(T)=0;其中T為緊急制動的時間,將車輛的初始狀態和最終狀態代入多項式表征的運動參量中,可求解出五項式系數的表達式。

約束條件:(1)避免碰撞的發生,以目標車靜止的追尾工況為例,如圖6所示約束條件為Sf<D;(2)最大制動減速度絕對值小于8m/s2。

圖6 制動工況約束條件的示意圖

通過優化算法可最終計算出車輛的運動狀態參量。為了驗證上述駕駛員緊急制動數學模型的有效性,隨機選取2個駕駛員緊急制動工況進行擬合,圖7和圖8分別為制動工況1#的擬合情況(兩車的相對車速為28km/h,碰撞時間TTC=0.8s),圖9和圖10為制動工況2#的擬合情況(兩車的相對速度為54km/h,TTC=1.8s)。

圖7 制動工況1#的制動減速度曲線

圖8 制動工況1#的制動減速度梯度曲線

圖9 制動工況2#的制動減速度曲線

圖10 制動工況2#的制動減速度梯度曲線

通過對比,建立的駕駛員緊急制動模型可以很好地復現駕駛員在緊急工況下的制動操作。

3.2 駕駛員緊急轉向模型的建立

對緊急轉向避撞變道過程進行分析,如圖11所示,整個緊急轉向變道過程可分成3個階段:避撞階段、側移階段以及穩定階段。本文中研究集中在避撞階段(從駕駛員開始轉向避撞行為的時刻T0開始,以轉向盤轉角達到最大值的時刻T1為終點)的轉向行為。

圖11 緊急轉向變道避撞階段劃分

為了研究避撞階段駕駛員緊急轉向特征,圖12為33例緊急變道工況避撞階段的轉向盤最高角速度和轉向盤轉角變化的線性關系圖,線性相關系數為0.98,線性比例約為2.21。

圖12 轉向盤最高角速度和轉角的線性關系

上述線性關系表明,避撞階段的駕駛員緊急轉向避撞可用開環模型來表示,本文中利用高斯函數建立轉向盤角速度的表達式為

駕駛員在緊急轉向變道的避撞階段,同樣是以保證轉向穩定性,優化目標為

上已述及,駕駛員的轉向盤角速度可通過高斯函數來表達;而利用車輛自行車模型,可以推導車身的其他運動狀態參量。

不失一般性,考慮圖13所示將兩車都視作相同大小矩形,有如下兩個約束條件。

圖13 緊急轉向避撞的示意圖

(1)避免與前方障礙物發生碰撞,如圖13所示,以向左轉向避撞為例,在本車通過目標車時,目標車左后角點(x0,y0)到本車右側邊界距離大于一定的安全距離Dsm,即

式中:c為本車右側邊界所在直線斜率,c=tan(?);(xrf,yrf)為本車矩形邊界的右上頂點,可對車身的狀態參數進行推導。

(2)緊急轉向的最大側向加速度絕對值小于6m/s2。

通過優化算法可最終計算出車輛的運動狀態參量。為了驗證上述駕駛員緊急轉向數學模型的有效性,隨機選取2個駕駛員緊急轉向工況進行擬合,圖14為轉向工況1#的擬合情況(兩車的相對車速為57km/h,TTC=1.35s),圖15為轉向工況2#的擬合情況(兩車的相對速度為70km/h,TTC=1.3s)。通過對比可知,建立的駕駛員緊急轉向模型可以很好地復現駕駛員在緊急工況下的轉向操作。

圖14 轉向工況1#的轉向盤角速度和側向加速度曲線

圖15 轉向工況2#的轉向盤角速度和側向加速度曲線

4 危險估計算法的研究

本文中以目標車靜止、后車追尾的工況為例,進行危險估計算法的研究。通過上節駕駛員緊急制動和轉向模型,可模擬出從車輛初始危險狀態到安全狀態的駕駛員緊急操作。模型預測的駕駛員參考操作,可判斷該操作是否超出了駕駛員正常駕駛的可行域范圍,并通過表1和表2中的表征制動和轉向操作緊急程度的閾值參數來估計目前車輛行駛狀態的危險程度。

根據文獻[6]和文獻[7]中的研究成果,駕駛員的緊急操作行為與TTC倒數TTC-1相關,本文中采用TTC-1進行危險估計的分析。以本車30和70km/h速度接近靜止的目標車為例進行計算,表3為達到表1中制動操作緊急程度1~3條件下的緊急制動工況。

表3 達到表1的緊急程度_1~3條件的緊急制動工況

圖16和圖17為30與70km/h情況下,不同工況下模擬車輛緊急制動減速度和制動減速度梯度的模擬曲線,分別為達到表1中制動操作緊急程度的3個等級的工況操作曲線。

圖16 70km/h車速下各工況的制動減速度和梯度曲線

圖17 30km/h車速下各工況的制動減速度和梯度曲線

隨著TTC-1增加,兩車之間的距離D減少,車輛最大制動減速度及其梯度的絕對值都逐漸變大,危險等級不斷增加。進一步分析可知,在低速情況下危險程度主要受最大制動減速度梯度的影響,在高速情況下危險程度主要是受最大制動減速度的影響。

用同樣的方法分析兩種車速下緊急轉向的駕駛行為,表4示出達到表2中的操作緊急程度1~3條件下的緊急轉向工況。

表4 達到表2的緊急程度_1~3條件的緊急轉向工況

圖18和圖19為在表4中不同工況下自車轉向盤角速度和側向加速度的模擬曲線。分別為達到表2中轉向操作緊急程度的3個等級的工況操作曲線。

圖18 70km/h各工況的轉向盤角速度和側向加速度曲線

通過分析隨著TTC-1的不斷增加,相對距離不斷減少,所需的轉向盤最高角速度和車輛最大側向加速度不斷增加,緊急轉向危險程度也不斷增加。

圖19 30km/h各工況的轉向盤角速度和側向加速度曲線

綜合比較緊急制動和緊急轉向操作,在速度為70km/h時,緊急轉向避撞所需的縱向距離更小;在低速30km/h時,制動避撞所需的縱向距離更小。

為了更好地說明制動和轉向操作在不同車速下的區別,以最高操作緊急等級緊急程度3作為約束條件,計算并對比緊急制動和轉向的臨界避撞TTC-1,如圖20所示。

圖20 緊急制動和轉向的臨界TTC-1曲線

在相對車速較高(大于40km/h)時,轉向避撞的臨界TTC-1更高,說明此時雖然制動無法避免碰撞,但通過轉向還能夠避免碰撞,因此緊急轉向避撞在相對車速較高時相比于緊急制動可以提供更好的避撞效果。而在相對車速較低(小于40km/h)時,緊急制動相對于緊急轉向避撞可以提供更好的避撞效果。

通過駕駛員緊急工況的制動轉向模型,分別估算緊急制動和轉向操作,并通過制動和轉向緊急程度的參數來估計目前車輛行駛狀態的危險程度。利用表1和表2建立的制動和轉向操作不同緊急程度的劃分,分別計算不同車速在不同危險等級下的臨界TTC-1值,圖21為不同車速下緊急制動和緊急轉向的臨界TTC-1值的曲線圖。

圖21 制動和轉向操作不同緊急程度的臨界TTC-1曲線

根據上述研究,結合圖21,在速度較低(小于40km/h)時,緊急制動有更好的避撞效果,故本文中危險估計算法在低速情況下以緊急制動的臨界TTC-1作為危險等級區域劃分的依據;在速度較高(大于40km/h)時,緊急轉向有更好的避撞效果,且隨著車速增加駕駛員通過轉向避撞有增多的趨勢,故在高速階段以緊急轉向的臨界TTC-1作為危險等級區域劃分的依據,同時考慮到高速階段發生碰撞的危險程度較高,故增加兩檔以緊急制動的臨界TTC-1來劃分的危險等級。最終基于自然駕駛工況的危險估計算法的區域劃分如圖22所示。

由圖22可見,按照TTC-1與車速的關系劃分危險等級區域。I0區域,未能檢測到碰撞發生的危險;I1區域,碰撞危險程度低,提示駕駛員進行制動或轉向來避撞;I2區域,碰撞危險程度低,提示駕駛員進行轉向來避撞;II區域,危險程度較高,一級預警;III區域,危險等級極高,二級預警;IV區域,碰撞危險極高,碰撞即將發生甚至無法避免。

圖22 基于自然駕駛工況的危險估計的區域劃分

5 結論

本文中基于我國自然駕駛工況數據,提出了一種汽車碰撞危險估計算法。首先分析了緊急制動和轉向的駕駛員輸入特性,對緊急操作的等級進行了劃分;篩選了緊急工況,建立了駕駛員緊急制動和轉向的駕駛員模型,并進行了驗證。通過駕駛員模型,并結合制動和轉向操作的緊急程度,分析了不同速度下駕駛員的緊急制動和轉向避撞操作的特征。結果表明,在低速階段緊急制動操作所需避撞的縱向距離小,而在高速階段緊急轉向操作的避撞效果更好。結合不同速度情況下的避撞特征,提出了汽車碰撞危險估計算法,以TTC-1作為危險等級的參數,對汽車碰撞的危險程度進行分級區域的劃分。

本文中建立的危險估計算法中的閾值是基于自然駕駛工況中的數據,與我國駕駛員駕駛行為特征相符,并同時考慮制動和轉向兩種避撞措施,可用于設計針對我國用戶的汽車避撞系統。

[1] EIDEHALL A.Tracking and threat assessment for automotive collision avoidance[D].Linkoping:Department of Electrical Engineering of Linkoping University,2007.

[2] ZHANG Y,ANTONSSON E K,GROTE K.A new threat assessment measure for collision avoidance systems[C].Proceeding of the IEEE ITSC 2006.Toronto:IEEE Press,2006:968-975.

[3] KITAJIMA S,MARUMO Y,HIRAOKA T,et al.Comparison of evaluation indices concerning estimation of driver's risk perceptionrisk perception of rear-end collision to a preceding vehicle[J].Review of Automotive Engineering,2009,30(2):191.

[4] LEE K,PENG H.Evaluation of automotive forward collision warning and collision avoidance algorithms[J].Vehicle System Dynamics,2005,43(10):735.

[5] NAJM W G,SMITH J D,YANAGISAWA M.Pre-crash scenario typology for crash avoidance research[R].Washington,DC:National Highway Traffic Safety Administration,2007.

[6] MOON S,YI K.Human driving data-based design of a vehicle adaptive cruise control algorithm[J].Vehicle System Dynamics,2008,46(8):661.

[7] KIEFER R J,LEBLANC D J,FLANNAGAN C A.Developing an inverse time-to-collision crash alert timing approach based on drivers last-second braking and steering judgments[J].Accident Analysis&Prevention,2005,37(2):295.

[8] BIRAL F,LIO M D,LOT R,et al.Vehicle and driver modeling and threat assessment for driving support functions[C].Atti XX Congresso dell’Associazione Italiana di Meccanica Teorica e Applicata,2011.

[9] MATTIAS M,COELINGH E,et al.Model-based threat assessment for avoiding arbitrary vehicle collisions[J].IEEE Transactions on Intelligent Transportation Systems,2010,11(3):658-669.

[10] BROWN T L,LEE J D,MCGEHEE D V.Human performance models and rear-end collision avoidance algorithms[J].Human Factors: The Journal of the Human Factors and Ergonomics Society,2001,43(3):462.

[11] BERTOLAZZI E,BIRAL F,LIO M D,et al.Supporting drivers in keeping safe speed and safe distance:The SASPENCE subproject within the european framework programme 6 integrating project PReVENT[J].Intelligent Transportation Systems IEEE Transactions on,2010,11(3):525-538.

[12] KIEFER R J,CASSAR M T,FLANNAGAN C A,et al.Forward collision warning requirements project final report-task 1[R]. Washington D C:National Highway Traffic Safety Administration,2003.

[13] HUGEMANN D W,NICKEL D M,MORAWSKI I,et al.Longitudinal and lateral accelerations in normal day driving[C].The 6thInternational Conference of The Institute of Traffic Accident Investigators.London:The Institute of Traffic Accident Investigator,2003:1-8.

[14] BRACKSTONE M,MCDONALD M.Driver headway:how close is too close on a motorway?[J].Ergonomics,2007,50(8): 1183.

[15] HIRAOKA T,TANAKA M,KUMAMOTO H,et al.Collision risk evaluation index based on deceleration for collision avoidance (First Report)[J].Review of Automotive Engineering,2009,30 (4):429.

A Study on Risk Assessment Algorithm Based on Natural Driving Data

Wu Bin1,2,Zhu Xichan1,Shen Jianping2&Li Lin3
1.College of Automotive Studies,Tongji University,Shanghai 201804;2.National Center of Supervision and Inspection on Motor Vehicle Products Quality(Shanghai),Shanghai 201805;3.Shanghai International Automobile City(Group)Co.Ltd.,Shanghai 201805

A vehicle collision risk assessment algorithm based on natural driving data.is proposed in this paper.Firstly the driver input characteristics for braking and steering under natural driving condition are analyzed,the feature parameters of driver behavior characterizing the emergency extent of collision avoidance by braking and by steering are studied,and the driver models for emergency braking and emergency steering are established.Then on this basis,the features of collision avoidance by braking and by steering are analyzed.Finally a risk assessment algorithm based on natural driving data is proposed.The simulation results show that the risk assessment algorithm proposed can concurrently consider collision avoidance by both braking and steering and correspond to emergency collision avoidance behavior of driver under real traffic environment,so can be used for the development of control strategy for vehicle collision avoidance system.

natural driving data;collision avoidance;driver behavior;risk assessment

10.19562/j.chinasae.qcgc.2017.08.009

?上海市科委研發平臺項目(16DZ2291000)資助。

原稿收到日期為2016年11月23日,修改稿收到日期為2017年2月16日。

朱西產,教授,博士生導師,E-mail:xczhu@163.com。

猜你喜歡
駕駛員
基于高速公路的駕駛員換道意圖識別
基于眼動的駕駛員危險認知
駕駛員安全帶識別方法綜述
如何提高汽車駕駛員預防事故的能力
活力(2019年22期)2019-03-16 12:46:58
基于matlab的純電動客車駕駛員行為分析
汽車電器(2018年1期)2018-06-05 01:22:54
駕駛員擅自偷開公車肇事公司仍需擔責
紅土地(2016年6期)2017-01-15 13:45:55
美國駕駛證申請與新駕駛員培訓簡介
汽車與安全(2016年5期)2016-12-01 05:22:10
起步前環顧四周是車輛駕駛員的義務
公民與法治(2016年4期)2016-05-17 04:09:26
基于車輛軌跡的駕駛員自由換道行為研究
總裝某汽車營:想當駕駛員直接“上車”考試
解放軍生活(2013年5期)2013-04-29 21:51:50
主站蜘蛛池模板: 一级毛片免费观看久| 免费人成黄页在线观看国产| 国产精品伦视频观看免费| 国产福利在线观看精品| 亚洲成人高清无码| 麻豆精品视频在线原创| 中日无码在线观看| 久久精品无码一区二区日韩免费| 无码内射在线| 国产成人h在线观看网站站| 亚洲九九视频| 91成人在线观看| 国产日韩精品一区在线不卡| 成人午夜天| 中文字幕丝袜一区二区| 夜夜爽免费视频| 国产清纯在线一区二区WWW| 国产a网站| 伊大人香蕉久久网欧美| 91青草视频| 婷婷伊人五月| 亚洲国产日韩在线观看| 无码日韩视频| 少妇精品网站| 中文字幕欧美成人免费| 亚洲伊人久久精品影院| 99久久成人国产精品免费| 狂欢视频在线观看不卡| 欧美日韩亚洲综合在线观看| 五月天香蕉视频国产亚| 国产在线视频自拍| 99re精彩视频| 亚洲婷婷六月| 日本午夜精品一本在线观看| 国产第四页| 久久免费观看视频| 欧美中文字幕一区| 天天综合网色中文字幕| 色视频国产| 97视频在线精品国自产拍| 伊人大杳蕉中文无码| 四虎国产永久在线观看| 国产美女无遮挡免费视频网站 | 国产精品第一区在线观看| 日韩欧美高清视频| 国产综合无码一区二区色蜜蜜| 成人亚洲国产| 欧美色丁香| 国产免费a级片| 在线播放91| 欧美一区二区自偷自拍视频| 欧美成人看片一区二区三区 | 色综合久久无码网| 国产日本视频91| 久久狠狠色噜噜狠狠狠狠97视色| 国产精品男人的天堂| 狠狠色丁婷婷综合久久| 国产a在视频线精品视频下载| 欧美一区二区三区欧美日韩亚洲| 米奇精品一区二区三区| 精品一区二区三区视频免费观看| 亚洲中文精品人人永久免费| 91九色最新地址| 日韩在线观看网站| 国产成人综合网| 亚亚洲乱码一二三四区| 欧美日韩中文国产va另类| 男女男免费视频网站国产| 片在线无码观看| 午夜福利在线观看入口| 毛片网站在线看| 呦女亚洲一区精品| 99精品国产自在现线观看| 欧美高清国产| 精品久久高清| 国产成人精品高清不卡在线| 中国一级特黄视频| 亚洲综合18p| 久久这里只有精品国产99| 精品一区二区三区波多野结衣| 国产办公室秘书无码精品| 免费女人18毛片a级毛片视频|