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

    智能車自主換道超車控制方法研究

    2023-06-07 00:22:18許小偉崔金融嚴運兵馬霽旻王亞瑋
    武漢科技大學學報 2023年3期
    關(guān)鍵詞:主車勢場車道

    許小偉,崔金融,嚴運兵,馬霽旻,王亞瑋

    (1.武漢科技大學汽車與交通工程學院,湖北 武漢,430065; 2.智新科技股份有限公司,湖北 武漢,430056;3.華中科技大學電氣與電子工程學院,湖北 武漢,430074)

    智能車換道超車是一種常見的駕駛行為,其中換道時刻的選擇和高速換道行駛過程中的安全性需要重點關(guān)注,必須解決換道決策、路徑規(guī)劃以及軌跡跟蹤等多個關(guān)鍵技術(shù)問題[1]。

    換道行為涉及到多個車道,由于人工勢場理論簡單且直觀,因此被廣泛用于量化道路的動態(tài)環(huán)境。張朋[2]利用人工勢場法對道路的動態(tài)環(huán)境復雜度進行建模,并對提出的量化模型進行實驗標定。李林恒等[3]為了直觀地表達智能車在換道過程中面臨的風險,選用人工勢場法建立換道模型,并利用國內(nèi)道路環(huán)境進行模型標定。朱乃宣等[4]借助人工勢場理論量化前車產(chǎn)生的風險場場強,依據(jù)駕駛風格來設計不同的規(guī)劃閾值而達到個性化換道的目的,但是換道觸發(fā)時刻僅考慮前車的影響,未考慮換道對相鄰車道正常行駛車輛造成的影響。

    軌跡跟蹤的目的是使車輛按照規(guī)劃的路徑行駛。目前智能車軌跡跟蹤方法主要有PID控制、模型預測控制和LQR控制等,其中LQR控制能考慮到多項性能指標且實時性好,因此得到廣泛應用。陳亮等[5]構(gòu)建了加入前饋控制的LQR控制器,有效減少了跟蹤誤差。胡杰等[6]根據(jù)車速變化對LQR權(quán)重系數(shù)進行模糊調(diào)節(jié),并驗證了控制器在實車環(huán)境中的控制性能。高琳琳等[7]對LQR控制器的權(quán)重矩陣進行了自適應調(diào)整。然而,模糊調(diào)節(jié)和自適應調(diào)整權(quán)重都是依靠經(jīng)驗設計,并不能保證輸出最優(yōu)的控制量。

    針對上述問題,本文基于人工勢場理論量化道路風險場,建立換道規(guī)則,并借助最小縱向安全距離來判斷換道的安全性,同時通過改進的LQR控制器對換道軌跡進行跟蹤。本文最后利用Carsim和Simulink搭建聯(lián)合仿真平臺,來驗證所提出的換道超車決策的有效性以及軌跡跟蹤控制器的效果。

    1 道路風險場設計

    本文基于人工勢場理論計算車輛正常行駛過程中遇到的風險水平,并建立道路風險場,包括障礙車輛風險場和引導風險場。障礙車輛風險場考慮了障礙車輛的自身屬性,引導風險場考慮了主車與周圍車輛的相對位置和相對速度。

    1.1 障礙車輛風險場

    障礙車輛風險場涉及換道車輛與周圍車輛的相對運動狀態(tài)。由實際交通流可知,車輛換道過程中,在縱向上,換道車輛接近周圍車輛時,風險場會不斷增強,且車輛相對距離越近,場強的增長幅度越大;在橫向上,允許換道車輛以較小的換道距離超越前車,且前車的風險場應盡可能地不影響相鄰車道。綜上,選用二次高斯函數(shù)模型[8]來表示周圍車輛的場強,建立障礙車輛勢場Uobst。

    Uobst=

    (1)

    ?x,obst=kxLobst

    (2)

    ?y,obst=kyWobst

    (3)

    式中:x、y分別為換道車輛質(zhì)心的縱、橫向坐標;xobst、yobst分別為障礙車輛質(zhì)心的縱、橫向坐標;Aobst為障礙車輛的場強幅值,本文取值為1;?x,obst、?y,obst分別為障礙車輛產(chǎn)生的勢場在縱向和橫向上的影響范圍;Lobst、Wobst分別為障礙車輛的長度和寬度;kx、ky分別為障礙車輛在x、y方向上的縮放系數(shù),考慮到實際超車過程中障礙車輛只對同車道后車產(chǎn)生影響,而對同向相鄰的車道影響很小甚至沒有影響,故取kx=1,ky=0.5。

    構(gòu)建圖1所示的換道超車場景,道路環(huán)境為同向行駛的三車道,規(guī)定向右行駛為正方向,文中均以A車為主車構(gòu)建行駛風險場,其余車輛視為周圍的障礙車輛,各車輛行駛狀態(tài)參數(shù)見表1。

    表1 車輛狀態(tài)參數(shù)

    圖1 智能車換道超車場景示意圖

    依據(jù)圖1所示場景形成的勢場如圖2所示,可以看出,障礙車輛風險場在縱向上影響范圍較大,而在橫向上影響范圍較小;其次,隨著主車與障礙車相對距離的減小,主車所受到的風險也越大。由此可知,使用二次高斯函數(shù)建立障礙車輛勢場滿足換道場景的要求。

    (a)障礙車輛勢場3D圖

    1.2 引導風險場

    為了避免兩車發(fā)生碰撞事故,增加一個額外的引導勢場,更好地引導在慢速行駛的障礙車后方的車輛提前換道行駛。選用的引導勢場函數(shù)表達式為

    Uext=Aobst+Fx(x-xobst)-Fy(y-yobst)2,

    ?x∈[-S,0]

    (4)

    式中:S為引導勢場在縱向上的影響范圍;Fx、Fy分別為引導勢場在縱向和橫向上的影響系數(shù),計算方式如下:

    (5)

    (6)

    引導勢場只有在車輛靠近障礙車輛的時候才會產(chǎn)生實際意義,故引導勢場的強度必然不小于零,小于零的勢場值用零代替。影響范圍公式為:

    S=0.5Lobst+Smin+tavrel

    (7)

    式中:Smin為引導勢場產(chǎn)生影響的最短距離;vrel為兩車之間的相對速度;ta為相對速度減小時的接近時距。

    圖1所示換道場景的引導勢場見圖3,可以看出,引導勢場影響范圍隨著兩車相對速度絕對值的增加而擴大,且距離障礙車輛質(zhì)心越近,引導勢場值越大。

    (a)引導勢場3D圖

    1.3 總的道路風險場

    綜上所述,車輛總的道路風險場為障礙車輛勢場與引導勢場的整合:

    Utotal=Uobst+Uext

    (8)

    2 基于道路風險場的換道決策模型

    2.1 基于風險場的換道意圖

    智能車的換道動機是為了追求更高的通行效率以及更優(yōu)的行駛狀態(tài),即智能車換至目標車道可獲得足夠的行駛空間以加速到期望車速。首先,比較主車當前所處位置與相鄰車道中心線位置的風險場,若是相鄰車道中心線位置的風險場強度低于主車所處位置的場強,則認為主車有向相鄰車道的換道動機,相反則認為本車保持當前車道行駛更佳。但是,在實際情況下,若兩處的風險場場強的差值很小,一旦車輛狀態(tài)發(fā)生微小變化,就會引起場強變化,導致車輛頻繁產(chǎn)生換道意圖,反而影響通行效率。因此,本文定義一個規(guī)劃閾值P,只有當主車所處位置與相鄰車道中心線位置的場強差值超過預設的規(guī)劃閾值,才認為換道后可獲取更優(yōu)行駛狀態(tài),即滿足式(9)時,車輛產(chǎn)生向相鄰車道的換道超車意圖。

    Utotal(x,yi)>Utotal(x,yj)+P

    (9)

    式中:x為主車所處位置縱向坐標,yi為主車所處位置橫向坐標,yj為相鄰車道中心線橫向坐標。

    2.2 換道安全性分析

    換道意圖產(chǎn)生后,為了保證換道安全性,本文以最小縱向安全距離作為判據(jù),其中包括本車道內(nèi)最小縱向安全距離以及目標車道前、后車最小縱向安全距離。由文獻[9]統(tǒng)計分析可知,車輛切入場景的縱向速度變化很小,因此本文假設換道過程中主車與周圍車輛都保持恒定的縱向速度行駛。

    (1)本車道最小縱向安全距離

    在換道過程中,智能車可通過傳感器設備獲得與主車相同車道的前車的速度以及二者之間的縱向車距,保持合理的車距進行換道。與前車的最小縱向安全距離MSS1計算如下:

    (10)

    式中:v為主車當前車速,vf為前車速度,tC為換道臨界時間,dmin為最小安全距離的臨界值。

    (2)目標車道前、后車最小縱向安全距離

    在更為復雜的交通環(huán)境中,車輛換道至有正常行駛車輛的目標車道,因此要保證換道車輛與目標車道上車輛之間的安全性。換道起始時刻,智能車獲取目標車道車輛的位置與速度信息,進行分析計算,保證換道結(jié)束時滿足安全車距的要求。換道車輛與目標車道前、后車之間的最小縱向安全距離MSS2、MSS3分別為:

    (11)

    (12)

    式中:vi,f、vi,r分別代表目標車道前、后車的縱向速度。

    2.3 自主換道觸發(fā)流程

    在車輛產(chǎn)生換道意圖并判斷換道安全后,決策系統(tǒng)發(fā)出觸發(fā)換道指令,換道行為觸發(fā)流程如圖4所示。

    圖4 換道行為觸發(fā)流程

    3 軌跡跟蹤控制器

    3.1 軌跡跟蹤誤差力學模型

    本文利用LQR控制器對智能車規(guī)劃好的軌跡進行跟蹤。首先建立一個適用于控制器設計的二自由度動力學跟蹤誤差模型[10],如圖5所示。

    圖5 二自由度動力學跟蹤誤差模型

    圖5中:δ為前輪轉(zhuǎn)角;lf、lr分別為車輛質(zhì)心到前、后軸的距離;ed為車輛質(zhì)心到參考軌跡的最短距離,即位置偏差;φ為車輛航向角,φdes為參考航向角;ω為橫擺角速度;vx、vy分別為車輛縱、橫向速度。

    根據(jù)動力學模型得到二自由度微分方程為:

    (13)

    式中:m為車輛質(zhì)量,Iz為轉(zhuǎn)動慣量,Cf、Cr分別為車輛前、后輪側(cè)偏剛度。

    定義航向誤差為

    eφ=φ-φdes

    (14)

    設參考軌跡的半徑為Rdes,可得到理想的橫向加速度為

    (15)

    智能車產(chǎn)生的橫向加速度為

    (16)

    (17)

    (18)

    綜上所述,得出如下公式:

    (19)

    (20)

    將式(19)和式(20)轉(zhuǎn)換為狀態(tài)空間表達形式,則為

    (21)

    3.2 預瞄模型

    軌跡跟蹤控制器是基于誤差模型計算出前輪轉(zhuǎn)角的輸入,存在一定的滯后性,故在本文方法中對車輛狀態(tài)進行一定采樣時間ts的預瞄。設車輛在預瞄時間內(nèi)勻速直線行駛,預瞄后車輛的航向角為φpre、坐標為(Xpre,Ypre),有

    (22)

    式中:(X0,Y0)為軌跡規(guī)劃點的坐標。

    3.3 LQR控制器

    本文利用LQR[11-12]對智能車換道軌跡進行跟蹤控制。首先明確線性系統(tǒng)的性能指標函數(shù)J為

    (23)

    式中:Q和R為控制器的權(quán)重矩陣。

    為了使性能指標J達到最優(yōu),采用變分法進行求解。計算得到反饋矩陣為

    K=(R+BTPB)-1BTPA

    (24)

    其中P為如下Riccati方程的正定解:

    P=Q+ATPA-ATPB(R+BTPB)-1BTPA

    (25)

    最終得到LQR控制器的最優(yōu)控制律為

    U(t)=-KX(t)

    (26)

    式中:K=[k1,k2,k3,k4]為LQR控制器增益。

    3.4 前饋控制器

    將式(26)帶入式(21)可得

    (27)

    由此可見,系統(tǒng)存在穩(wěn)態(tài)誤差,即存在無法消除的位置偏差和航向偏差,為此設計前饋控制器消除該誤差。設前饋控制輸出的前輪轉(zhuǎn)角為δf,即控制量輸入U(t)=-KX(t)+δf。

    根據(jù)終值定理可得系統(tǒng)的穩(wěn)態(tài)誤差為

    e=-(A-BK)-1(Bδf+B1C)

    (28)

    為保證橫向偏差穩(wěn)態(tài)值趨于零,設計前饋控制律為

    (29)

    式中:kref為道路參考曲率。

    將前饋控制律與LQR反饋控制律相結(jié)合,得到跟蹤算法的前輪轉(zhuǎn)角輸出量。

    3.5 遺傳算法優(yōu)化權(quán)重矩陣參數(shù)

    式(23)中的權(quán)重矩陣Q和R在車輛跟蹤控制中起著非常重要的作用。尋找最優(yōu)的Q和R是一項復雜的工作,往往需要大量試錯才能得出可行解。本文利用遺傳算法來優(yōu)化LQR控制器參數(shù)[12-13]。

    鑒于二自由度動力學誤差模型有且僅有前輪轉(zhuǎn)角一個控制量,因此確定控制權(quán)重矩陣R=1。Q中的對角線參數(shù)代表了各項誤差的權(quán)值,在車輛軌跡跟蹤過程中更加關(guān)注位置誤差和航向誤差,因此將Q矩陣設置為

    (30)

    下面利用遺傳算法對矩陣Q中參數(shù)a、b進行優(yōu)化。為了提高跟蹤的精準性,同時考慮位置誤差和航向誤差,因此設目標適應度函數(shù)為

    (31)

    根據(jù)目標適應度函數(shù)最小化原則,LQR控制器參數(shù)優(yōu)化流程如圖6所示。具體步驟如下:

    圖6 采用遺傳算法優(yōu)化LQR 控制器的流程

    步驟1在0~10的約束范圍內(nèi)產(chǎn)生N個初始個體。

    步驟2將初始種群中的個體賦值給權(quán)重矩陣Q中的a、b,得到反饋增益矩陣K,求出最優(yōu)控制律,并輸入給車輛仿真模型,運行模型,計算適應度函數(shù)值。

    步驟3選取適應度數(shù)值小的個體進行迭代。

    步驟4通過保優(yōu)、交叉、變異等遺傳操作,保留父代和產(chǎn)生的子代,組成新的種群。

    步驟5重復步驟2~4,直至滿足終止條件。

    4 仿真分析

    為了驗證本文方法的有效性,利用Simulink和Carsim搭建聯(lián)合仿真平臺,針對車輛在快速路上自主換道行為進行仿真測試,車輛的部分結(jié)構(gòu)參數(shù)如表2所示。

    表2 車輛結(jié)構(gòu)參數(shù)

    4.1 軌跡跟蹤控制器性能分析

    在帶有前饋的LQR控制器基礎之上,分別采用經(jīng)驗法和遺傳算法確定控制器參數(shù),在相同的車速和路況條件下進行雙移線仿真測試,驗證軌跡跟蹤控制效果。

    在雙移線測試中,路面附著系數(shù)設為0.85,車速恒定為60 km/h,雙移線的軌跡方程[7]為

    (32)

    方案一 按照經(jīng)驗法選取控制效果最好的一組控制參數(shù):a=5,b=5。

    方案二 采用遺傳算法對控制器參數(shù)優(yōu)化,遺傳算法的終止代數(shù)設為20,初始種群個體數(shù)量為100,交叉和變異的概率分別取為0.75和0.15。多次運行程序,雖然每次結(jié)果不同,但是非常接近,最后得出a=9.9608,b=0.1233。

    兩個方案的測試結(jié)果如圖7所示。采用經(jīng)驗法確定參數(shù)的LQR控制器可以將軌跡跟蹤的位置誤差控制在0.9 m以內(nèi),航向誤差控制在0.12 rad以內(nèi);利用遺傳算法優(yōu)化后的LQR控制器可以將位置誤差控制在0.6 m以內(nèi),航向誤差控制在0.1 rad以內(nèi)??梢钥闯?通過遺傳算法優(yōu)化LQR控制器矩陣參數(shù)能減少軌跡跟蹤系統(tǒng)的超調(diào)量,提升跟蹤精度。

    (a)側(cè)向位置和位置誤差

    4.2 換道超車仿真與結(jié)果分析

    本文主要利用兩類場景、四種工況來驗證換道決策方法的有效性和跟蹤的準確性。第一類為簡單的換道場景,僅本車道前面存在低速障礙車輛;第二類為復雜的換道場景,存在多車道車輛影響。

    (1)第一類場景

    此類場景設計了兩種工況,通過設定不同的主車速度來驗證換道超車觸發(fā)時刻的合理性。

    工況一仿真場景如圖8所示,主車A行駛在快速路中間車道上,其前方有車輛B,而相鄰車道無任何車輛,該工況下各車輛的狀態(tài)參數(shù)如表3所示。

    表3 工況一的仿真參數(shù)設定

    圖8 工況一場景

    車輛換道超車過程如圖9所示。結(jié)合表3以及圖8~圖9可知,本車道上車輛B的行駛速度低于主車A,且主車A的期望速度高于其當前速度,為了追求更好的通行效率,主車A不斷靠近前車B,當主車所處位置的勢場值與相鄰道路中心勢場值之差大于閾值時,在t1時刻主車A產(chǎn)生換道意圖,由于左側(cè)車道沒有障礙車,選擇該車道為目標車道。下一步進行安全性分析,左側(cè)沒有障礙車且主車A與前車B的相對距離為41 m,大于本車道最小縱向安全距離MSS1=17 m的要求,故滿足換道的安全性要求。然后規(guī)劃出換道軌跡,t1時刻開始觸發(fā)換道指令,進行超車,t2時刻主車A駛?cè)肽繕塑嚨?換道結(jié)束。

    圖9 工況一的換道過程

    圖10所示為工況一的仿真結(jié)果,包括換道過程中車輛實際軌跡與規(guī)劃出的目標軌跡的對比、軌跡跟蹤誤差以及主車A的側(cè)向加速度和橫擺角。從圖10可以看出,整個換道過程相對平穩(wěn),軌跡跟蹤誤差小于0.3 m,側(cè)向加速度和橫擺角的峰值均在約束范圍之內(nèi)。

    (a)側(cè)向位置和位置誤差

    工況二仿真場景與圖8相同。相較于工況一,主車A的速度有所提升,車輛的具體位置與速度如表4所示。

    表4 工況二的仿真參數(shù)設定

    工況二的車輛換道超車過程如圖11所示。與工況一相比,主車A的速度有所提升,即主車與障礙車的速度差增加了,導致引導勢場的增強,兩車產(chǎn)生追尾的可能性也提高了,因此換道觸發(fā)時刻有所提前。在t1時刻車輛A產(chǎn)生換道意圖,且此時前、后車距離為53 m,滿足本車道最小縱向安全距離MSS1=26 m的要求。車輛A在t1時刻觸發(fā)換道指令進行換道,于t2時刻完成換道。

    圖11 工況二的換道過程

    工況二的仿真結(jié)果如圖12所示。與工況一相比,在車輛A行駛速度增加后,只有側(cè)向加速度與跟蹤誤差的峰值有所增加,且增幅很小,同時,換道過程中橫擺角的峰值幾乎沒有變化,這體現(xiàn)出本文提出的軌跡跟蹤控制方法具有良好的適應性。

    (a)側(cè)向位置和位置誤差

    (2)第二類場景

    下面通過兩種工況來驗證本文換道決策方法在復雜道路場景中的適用性。

    工況三仿真場景參照圖1,主車A行駛在快速路中間車道上,其前方有車輛B,且左右相鄰車道上均有車輛。該工況下各車輛的狀態(tài)參數(shù)如表5所示。

    表5 工況三的仿真參數(shù)設定

    工況三的換道過程如圖13所示。車輛A和B的位置及速度參數(shù)與工況二一致,依據(jù)工況二測試結(jié)果可知,A車在與前車B的距離為53 m處產(chǎn)生換道意圖,且C車所在車道速度更快,故被定為目標車道。此刻A車與C車相距36 m,滿足與目標車道后車之間的最小縱向安全距離MSS3=17 mm的要求,t1時刻觸發(fā)換道指令,t2時刻A車駛?cè)肽繕塑嚨劳瓿蓳Q道。

    圖13 工況三的換道過程

    工況四仿真場景與工況三相同,只是車輛C的速度由100 km/h提升到120 km/h,各車輛的狀態(tài)參數(shù)見表6。

    表6 工況四的仿真參數(shù)設定

    圖14所示為工況四的換道超車過程。相較于工況三,車輛C速度提升后,主車A受環(huán)境中其他車輛的位置和行駛速度的影響,換道觸發(fā)時刻有明顯的延后,這是車輛綜合考慮換道意圖和換道安全性的結(jié)果。首先主車A在中間車道行駛,不斷靠近障礙車輛B,由工況三可知,A車在距B車53 m時產(chǎn)生換道意圖,此時A車與C車的縱向距離為21 m,低于最小安全換道距離MSS3=33 mm的要求,因此等待C車超越A車后,并且兩車的相對距離大于前車最小縱向安全距離MSS2=8 mm的條件時,A車開始換道。由工況四的仿真結(jié)果可知,本文換道超車控制方法充分考慮了復雜道路環(huán)境下的換道安全性。

    圖14 工況四的換道過程

    5 結(jié)語

    為了保障智能車自主換道超車的安全性,本文首先根據(jù)人工勢場理論量化車輛換道風險,建立了基于道路風險場的換道觸發(fā)決策模型,然后采用改進的LQR控制器對軌跡進行跟蹤,其中LQR控制參數(shù)通過遺傳算法進行優(yōu)化。本文最后通過Simulink和Carsim聯(lián)合仿真驗證了改進后控制器的精度提升,并且通過不同道路環(huán)境下的換道超車仿真試驗,證明了所提出的換道觸發(fā)決策模型的有效性和安全性。然而,本研究未考慮換道過程中車輛縱向速度的改變,且目前還都只是基于聯(lián)合仿真測試,下一步工作將針對本文方法進行實車驗證。

    猜你喜歡
    主車勢場車道
    北斗+手機實現(xiàn)車道級導航應用
    基于三次B樣條曲線擬合的主車軌跡預測算法的研究
    基于三次B樣條曲線擬合的主車軌跡預測算法的研究
    基于Frenet和改進人工勢場的在軌規(guī)避路徑自主規(guī)劃
    基于改進人工勢場方法的多無人機編隊避障算法
    避免跟車闖紅燈的地面車道線
    科學大眾(2021年9期)2021-07-16 07:02:50
    淺談MTC車道改造
    歐曼牽引車制動系統(tǒng)異常的故障處理
    時代汽車(2018年4期)2018-05-31 02:53:08
    主車與掛車分別設立第三者責任保險的賠償額如何確定
    山東青年(2017年7期)2018-01-11 16:09:15
    庫車坳陷南斜坡古流體勢場對陸相油氣運聚的控制
    美女主播在线视频| 极品少妇高潮喷水抽搐| 性色av一级| 在线观看一区二区三区激情| 在线播放无遮挡| 国产高清不卡午夜福利| 2021天堂中文幕一二区在线观| 高清av免费在线| 在线观看免费高清a一片| 日本黄色片子视频| kizo精华| 精品人妻偷拍中文字幕| 2018国产大陆天天弄谢| 亚洲最大成人手机在线| 国产v大片淫在线免费观看| kizo精华| 蜜臀久久99精品久久宅男| 高清日韩中文字幕在线| 嫩草影院新地址| 简卡轻食公司| 久久精品熟女亚洲av麻豆精品| 最近中文字幕2019免费版| 久久精品久久久久久噜噜老黄| 18禁在线播放成人免费| 亚洲欧美一区二区三区国产| 国产大屁股一区二区在线视频| 国产精品.久久久| 国产精品无大码| 啦啦啦在线观看免费高清www| 嫩草影院入口| 国产一级毛片在线| 日韩国内少妇激情av| 国产黄色免费在线视频| 少妇的逼好多水| 国产片特级美女逼逼视频| av黄色大香蕉| 三级经典国产精品| 国产一区二区三区av在线| 中国国产av一级| 51国产日韩欧美| 欧美潮喷喷水| 国产乱人视频| 少妇的逼好多水| 精品一区二区三区视频在线| 大又大粗又爽又黄少妇毛片口| 欧美日韩一区二区视频在线观看视频在线 | 亚洲精品日本国产第一区| 久久久久久久大尺度免费视频| 在线播放无遮挡| 中文字幕av成人在线电影| 国产成人精品婷婷| 人妻少妇偷人精品九色| 人妻 亚洲 视频| 永久免费av网站大全| 国产精品人妻久久久影院| 亚洲精品国产av成人精品| av黄色大香蕉| 日本-黄色视频高清免费观看| 99热网站在线观看| 国产在视频线精品| 国产一区有黄有色的免费视频| 插逼视频在线观看| 22中文网久久字幕| 国产精品三级大全| 亚洲四区av| 久久精品国产鲁丝片午夜精品| 人人妻人人澡人人爽人人夜夜| 视频区图区小说| 一个人观看的视频www高清免费观看| 免费黄频网站在线观看国产| 人妻一区二区av| 免费看a级黄色片| 99九九线精品视频在线观看视频| 亚洲精品中文字幕在线视频 | 免费黄色在线免费观看| 欧美一级a爱片免费观看看| 日本黄大片高清| 午夜激情久久久久久久| 亚洲精品日本国产第一区| 成年av动漫网址| 熟妇人妻不卡中文字幕| 亚洲精品亚洲一区二区| 久热这里只有精品99| 欧美一级a爱片免费观看看| 亚洲人成网站在线播| 国产成人a区在线观看| 亚洲丝袜综合中文字幕| 国产一区有黄有色的免费视频| 亚洲精华国产精华液的使用体验| 亚洲丝袜综合中文字幕| 成人免费观看视频高清| 白带黄色成豆腐渣| 日本免费在线观看一区| 国产精品.久久久| 99精国产麻豆久久婷婷| 美女国产视频在线观看| 亚洲四区av| 成人美女网站在线观看视频| 国产女主播在线喷水免费视频网站| 欧美精品国产亚洲| av女优亚洲男人天堂| 久久久精品免费免费高清| 日韩av在线免费看完整版不卡| 激情五月婷婷亚洲| 97在线视频观看| 美女内射精品一级片tv| 99热这里只有是精品50| 在线观看av片永久免费下载| 听说在线观看完整版免费高清| 99热全是精品| 欧美+日韩+精品| av免费观看日本| 69av精品久久久久久| av线在线观看网站| 国产女主播在线喷水免费视频网站| 欧美另类一区| 国产一级毛片在线| av天堂中文字幕网| 2018国产大陆天天弄谢| 日本黄色片子视频| 欧美激情久久久久久爽电影| 亚洲精品自拍成人| 免费观看a级毛片全部| 国产精品人妻久久久影院| 五月天丁香电影| 国产欧美日韩精品一区二区| 亚洲欧美精品专区久久| 女人十人毛片免费观看3o分钟| 99热这里只有精品一区| 免费观看av网站的网址| 综合色av麻豆| 成人亚洲精品av一区二区| 成人亚洲欧美一区二区av| 最近最新中文字幕免费大全7| 99久久精品国产国产毛片| 久久精品熟女亚洲av麻豆精品| 亚洲av国产av综合av卡| 欧美激情在线99| 中文字幕人妻熟人妻熟丝袜美| 国产午夜福利久久久久久| 日本猛色少妇xxxxx猛交久久| 秋霞在线观看毛片| 亚洲色图av天堂| 亚洲一级一片aⅴ在线观看| 国产片特级美女逼逼视频| 欧美老熟妇乱子伦牲交| 中文字幕av成人在线电影| 免费不卡的大黄色大毛片视频在线观看| av国产精品久久久久影院| 久久久久精品性色| 九九爱精品视频在线观看| 成年女人在线观看亚洲视频 | 亚洲成人中文字幕在线播放| 国产亚洲最大av| 久久97久久精品| a级毛色黄片| 99热这里只有是精品在线观看| 中文在线观看免费www的网站| 直男gayav资源| 成人欧美大片| 免费av不卡在线播放| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 熟女av电影| 91午夜精品亚洲一区二区三区| 在线观看国产h片| 女的被弄到高潮叫床怎么办| 一级毛片黄色毛片免费观看视频| 麻豆国产97在线/欧美| 少妇人妻精品综合一区二区| 亚洲国产日韩一区二区| 日本午夜av视频| 青春草亚洲视频在线观看| 亚洲国产欧美人成| 好男人视频免费观看在线| 国产 精品1| 内地一区二区视频在线| 国产伦精品一区二区三区四那| 成人亚洲精品一区在线观看 | 国产探花在线观看一区二区| 伦精品一区二区三区| 在线观看av片永久免费下载| 亚洲综合色惰| 国产成人福利小说| 综合色丁香网| 99热6这里只有精品| 国产精品99久久99久久久不卡 | 三级国产精品欧美在线观看| 69人妻影院| 国产成人a区在线观看| 国产亚洲5aaaaa淫片| 色播亚洲综合网| 精品久久久噜噜| 久久99蜜桃精品久久| 亚洲欧洲日产国产| 在现免费观看毛片| av国产久精品久网站免费入址| 视频中文字幕在线观看| 色播亚洲综合网| 久久精品久久精品一区二区三区| 狂野欧美激情性bbbbbb| 午夜精品国产一区二区电影 | 亚洲自拍偷在线| 大码成人一级视频| 91在线精品国自产拍蜜月| 亚洲精品中文字幕在线视频 | 亚洲精品日韩av片在线观看| 我要看日韩黄色一级片| 69人妻影院| 极品教师在线视频| 亚洲欧美成人精品一区二区| 亚洲av.av天堂| 在线 av 中文字幕| 欧美人与善性xxx| 亚洲国产精品专区欧美| 久久久久国产精品人妻一区二区| 网址你懂的国产日韩在线| 九色成人免费人妻av| 嘟嘟电影网在线观看| 亚洲成人av在线免费| 国产精品三级大全| 极品少妇高潮喷水抽搐| 午夜日本视频在线| 国产成人a∨麻豆精品| 秋霞伦理黄片| 久久久久久伊人网av| 国产精品一区www在线观看| 国产欧美亚洲国产| 亚洲内射少妇av| 亚洲精品国产av蜜桃| 久久久久久九九精品二区国产| 少妇猛男粗大的猛烈进出视频 | 青青草视频在线视频观看| av黄色大香蕉| 国产午夜精品久久久久久一区二区三区| 丰满少妇做爰视频| 99re6热这里在线精品视频| av又黄又爽大尺度在线免费看| 建设人人有责人人尽责人人享有的 | 国产精品蜜桃在线观看| 1000部很黄的大片| 丰满少妇做爰视频| 国产成人freesex在线| 国产成人精品婷婷| 欧美成人午夜免费资源| 中文字幕av成人在线电影| 好男人视频免费观看在线| 欧美高清成人免费视频www| 99久久精品热视频| 成人鲁丝片一二三区免费| 国产高清不卡午夜福利| 欧美日韩视频高清一区二区三区二| 欧美97在线视频| 黄片wwwwww| 国产精品一区www在线观看| 亚洲成色77777| 我要看日韩黄色一级片| 一级毛片aaaaaa免费看小| 99久久精品热视频| 极品教师在线视频| 只有这里有精品99| 国产精品一区www在线观看| 国产男女超爽视频在线观看| 嫩草影院新地址| 人妻 亚洲 视频| 一级毛片 在线播放| 干丝袜人妻中文字幕| 国内精品美女久久久久久| 亚洲aⅴ乱码一区二区在线播放| 联通29元200g的流量卡| 亚洲美女视频黄频| 亚洲精品456在线播放app| 边亲边吃奶的免费视频| 日本一二三区视频观看| 成人亚洲精品一区在线观看 | 国产精品久久久久久av不卡| 精品少妇黑人巨大在线播放| av一本久久久久| 国产精品麻豆人妻色哟哟久久| 久久国产乱子免费精品| 99热6这里只有精品| 国产午夜福利久久久久久| 国产免费视频播放在线视频| 亚洲精品一二三| 看非洲黑人一级黄片| 爱豆传媒免费全集在线观看| 国产精品成人在线| 又大又黄又爽视频免费| 日韩制服骚丝袜av| 亚洲色图综合在线观看| 一级av片app| 国产成人freesex在线| eeuss影院久久| 欧美zozozo另类| 建设人人有责人人尽责人人享有的 | 欧美国产精品一级二级三级 | 免费电影在线观看免费观看| 欧美zozozo另类| 在线亚洲精品国产二区图片欧美 | 男的添女的下面高潮视频| 久久影院123| 成人亚洲精品av一区二区| 只有这里有精品99| 秋霞在线观看毛片| 男人舔奶头视频| 国产亚洲精品久久久com| 亚洲成色77777| 最近中文字幕高清免费大全6| 亚洲av中文字字幕乱码综合| 国产亚洲午夜精品一区二区久久 | 欧美少妇被猛烈插入视频| 日本一本二区三区精品| 久久久久久伊人网av| 直男gayav资源| 亚洲精品日韩在线中文字幕| 少妇人妻精品综合一区二区| 91午夜精品亚洲一区二区三区| 亚洲国产最新在线播放| 老司机影院毛片| 免费观看无遮挡的男女| 久久久久久久久久成人| 久久综合国产亚洲精品| 麻豆成人av视频| 亚洲最大成人av| 在线观看人妻少妇| 亚洲一级一片aⅴ在线观看| av线在线观看网站| 亚洲伊人久久精品综合| 亚洲图色成人| 菩萨蛮人人尽说江南好唐韦庄| 晚上一个人看的免费电影| 2021天堂中文幕一二区在线观| 国产精品一区www在线观看| 全区人妻精品视频| 久久久久九九精品影院| 亚洲第一区二区三区不卡| 丰满少妇做爰视频| 国产av码专区亚洲av| 王馨瑶露胸无遮挡在线观看| 在线a可以看的网站| 啦啦啦中文免费视频观看日本| 91狼人影院| av一本久久久久| 久久这里有精品视频免费| 国产亚洲91精品色在线| 国产色婷婷99| 国产成人精品婷婷| 午夜免费观看性视频| 亚洲av福利一区| 老师上课跳d突然被开到最大视频| 久久99热这里只有精品18| 中文字幕制服av| 日本色播在线视频| 如何舔出高潮| 又大又黄又爽视频免费| 欧美成人午夜免费资源| 美女国产视频在线观看| av黄色大香蕉| 黑人高潮一二区| 国产精品麻豆人妻色哟哟久久| www.色视频.com| 丝瓜视频免费看黄片| 久久久久精品久久久久真实原创| 亚洲图色成人| 中国三级夫妇交换| 男人添女人高潮全过程视频| 日韩av免费高清视频| 男人添女人高潮全过程视频| 欧美国产精品一级二级三级 | 亚洲精品一区蜜桃| 亚洲一区二区三区欧美精品 | 街头女战士在线观看网站| 一级黄片播放器| 久久热精品热| 最近2019中文字幕mv第一页| 国产乱人视频| 国产亚洲精品久久久com| 国产熟女欧美一区二区| 国产亚洲精品久久久com| 欧美xxⅹ黑人| 在线 av 中文字幕| 亚洲激情五月婷婷啪啪| 国产精品国产三级国产专区5o| 午夜日本视频在线| 又爽又黄无遮挡网站| 日韩欧美精品v在线| 国产亚洲一区二区精品| 精品久久久久久久久亚洲| 欧美极品一区二区三区四区| 3wmmmm亚洲av在线观看| 欧美极品一区二区三区四区| 99热网站在线观看| 国产成人精品一,二区| 国产成年人精品一区二区| 亚洲av成人精品一区久久| 丰满人妻一区二区三区视频av| 日本一本二区三区精品| 国国产精品蜜臀av免费| 赤兔流量卡办理| 亚洲国产av新网站| 三级国产精品片| 国产精品久久久久久久电影| 在线观看三级黄色| 国产综合懂色| 国产精品人妻久久久影院| 亚洲三级黄色毛片| 亚洲av福利一区| 精品一区二区三卡| 男女国产视频网站| 国产精品99久久久久久久久| 亚洲,一卡二卡三卡| 黄色一级大片看看| videossex国产| 精品熟女少妇av免费看| 欧美xxxx性猛交bbbb| av在线观看视频网站免费| 成人二区视频| 久久99热这里只有精品18| 国产综合懂色| 欧美国产精品一级二级三级 | 国产成人精品一,二区| 久久人人爽人人爽人人片va| 亚洲精品一二三| 日韩成人av中文字幕在线观看| 91狼人影院| 女人十人毛片免费观看3o分钟| 黄色日韩在线| 好男人在线观看高清免费视频| 交换朋友夫妻互换小说| 天堂俺去俺来也www色官网| 国产真实伦视频高清在线观看| 极品教师在线视频| 七月丁香在线播放| 久热这里只有精品99| 青春草亚洲视频在线观看| 十八禁网站网址无遮挡 | 国产综合精华液| 大香蕉久久网| 毛片女人毛片| 日日啪夜夜撸| 亚洲av不卡在线观看| 亚洲美女视频黄频| 中文字幕久久专区| 久久综合国产亚洲精品| 国产精品一区二区在线观看99| 亚洲综合精品二区| 交换朋友夫妻互换小说| 欧美3d第一页| 国产精品久久久久久精品电影| 天堂网av新在线| 国产视频首页在线观看| 大香蕉久久网| 在现免费观看毛片| 亚洲真实伦在线观看| 极品教师在线视频| 中国三级夫妇交换| 国产又色又爽无遮挡免| 国产成人免费无遮挡视频| 高清午夜精品一区二区三区| 久久精品夜色国产| eeuss影院久久| 国产亚洲av片在线观看秒播厂| 国产成人福利小说| 国产成人freesex在线| 欧美成人一区二区免费高清观看| 男的添女的下面高潮视频| 天堂俺去俺来也www色官网| 国内揄拍国产精品人妻在线| 亚洲第一区二区三区不卡| 国内精品宾馆在线| 午夜日本视频在线| 久久久久久久大尺度免费视频| av免费在线看不卡| 人人妻人人看人人澡| 美女xxoo啪啪120秒动态图| 国产欧美亚洲国产| 六月丁香七月| 制服丝袜香蕉在线| 男人爽女人下面视频在线观看| 国产永久视频网站| 嘟嘟电影网在线观看| 亚洲av.av天堂| 亚洲av中文字字幕乱码综合| 久久97久久精品| 波野结衣二区三区在线| 男人爽女人下面视频在线观看| 国产精品成人在线| videos熟女内射| 精品午夜福利在线看| 亚洲精品一二三| 丝袜美腿在线中文| 国产女主播在线喷水免费视频网站| 日本黄色片子视频| 一级毛片 在线播放| 80岁老熟妇乱子伦牲交| 国产免费视频播放在线视频| 国产老妇女一区| 国产精品久久久久久精品电影小说 | 国产 一区精品| 亚洲国产精品国产精品| 亚洲国产欧美人成| www.av在线官网国产| 亚洲国产欧美人成| 黄色配什么色好看| 男人和女人高潮做爰伦理| 2022亚洲国产成人精品| 亚洲av在线观看美女高潮| 国产成人福利小说| 欧美激情久久久久久爽电影| 国产片特级美女逼逼视频| 国产真实伦视频高清在线观看| 亚洲国产精品国产精品| 国产亚洲午夜精品一区二区久久 | 午夜福利视频1000在线观看| 蜜臀久久99精品久久宅男| 亚洲精品一区蜜桃| 国产精品爽爽va在线观看网站| 亚洲色图av天堂| 午夜爱爱视频在线播放| 最新中文字幕久久久久| 直男gayav资源| 国国产精品蜜臀av免费| www.av在线官网国产| 成年人午夜在线观看视频| 街头女战士在线观看网站| 中文在线观看免费www的网站| 久久ye,这里只有精品| av天堂中文字幕网| 少妇 在线观看| 视频区图区小说| 美女被艹到高潮喷水动态| 欧美激情在线99| 午夜福利网站1000一区二区三区| 天堂俺去俺来也www色官网| 高清日韩中文字幕在线| 亚洲国产精品999| 丰满少妇做爰视频| 欧美亚洲 丝袜 人妻 在线| av在线app专区| 听说在线观看完整版免费高清| 大陆偷拍与自拍| 人妻系列 视频| 少妇人妻久久综合中文| 久热这里只有精品99| 国产成人a∨麻豆精品| 久久久亚洲精品成人影院| 中文字幕免费在线视频6| 男女边吃奶边做爰视频| 成人亚洲精品一区在线观看 | 蜜桃亚洲精品一区二区三区| 国产成年人精品一区二区| 十八禁网站网址无遮挡 | 啦啦啦在线观看免费高清www| 日韩国内少妇激情av| 天美传媒精品一区二区| 国产精品一区二区性色av| 亚洲国产成人一精品久久久| 日韩av免费高清视频| 国产精品.久久久| 欧美xxxx性猛交bbbb| 免费av毛片视频| 色5月婷婷丁香| 高清欧美精品videossex| 亚洲熟女精品中文字幕| 免费少妇av软件| 欧美激情在线99| 日韩一本色道免费dvd| 内地一区二区视频在线| 亚洲内射少妇av| 一个人看的www免费观看视频| 国产精品熟女久久久久浪| 日本午夜av视频| 国产欧美日韩精品一区二区| 亚洲精品影视一区二区三区av| 汤姆久久久久久久影院中文字幕| 极品少妇高潮喷水抽搐| 少妇裸体淫交视频免费看高清| 三级国产精品欧美在线观看| 国产v大片淫在线免费观看| 精品久久久久久久久av| 夫妻性生交免费视频一级片| 少妇 在线观看| 中国三级夫妇交换| 在线观看一区二区三区激情| 亚洲,一卡二卡三卡| 成年女人看的毛片在线观看| 18禁动态无遮挡网站| 丰满乱子伦码专区| 成人亚洲精品一区在线观看 | 女的被弄到高潮叫床怎么办| 日本猛色少妇xxxxx猛交久久| 黄色日韩在线| 啦啦啦啦在线视频资源| 在线观看三级黄色| 久久久久国产网址| 国产美女午夜福利| 色吧在线观看| 国产一区二区三区综合在线观看 | 在线观看一区二区三区激情| 久久6这里有精品| 久久久午夜欧美精品| 又黄又爽又刺激的免费视频.| 69人妻影院| 亚洲一区二区三区欧美精品 | 欧美一区二区亚洲| 成人毛片60女人毛片免费| 高清毛片免费看| 久久综合国产亚洲精品| 欧美精品一区二区大全| 中文字幕人妻熟人妻熟丝袜美| 久久久久久久精品精品| 国产爱豆传媒在线观看| 国产大屁股一区二区在线视频| 国产av国产精品国产| 汤姆久久久久久久影院中文字幕| 亚洲在线观看片|