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

    考慮換道壓力增益的高速公路分流區(qū)跟馳模型

    2025-08-15 00:00:00楊龍海車婷婷馬文奎章錫俏陳霖
    重慶大學學報 2025年8期
    關(guān)鍵詞:交通流匝道標準差

    中圖分類號:U491.265 文獻標志碼:A 文章編號:1000-582X(2025)08-040-14

    doi:10.11835/j.issn.1000-582X.2024.005

    Research on car-following model for freeway diversion area considering lane-changing pressure gain

    YANG Longhai,CHE Tingting,MA Wenkui, ZHANG Xiqiao,CHEN Lin (School ofTransportation Science and Engineering,Harbin Institute ofTechnology,Harbin150o90,P.R.China)

    Abstract: To more accurately capture vehicle-following behavior in freeway diversion areas,this study proposes anenhanced car-following model by incorporating a lane-changing pressure gain factor into the full velocity diference model (FVDM).The proposed model accounts for both lane-changing pressure and lane-changing behavior. First,a linear stability analysis reveals that the stability region of the freeway diversion area diminishes as the lane-changing pressure gain factor increases.Second, using trajectory data from 92 lane-changing vehicles extracted from the NGSIM dataset,the improved car-folowing model is calibrated andvalidated.Simulation results demonstrate that the improved model more accurately reproduces vehicle speed and position.Compared with the original FVDM,the simulation error of the proposed pressure-based FVDM (P-FVDM) is reduced by 16% .Compared with the lane pressre FVDM(LP-FVDM),the proposed pressure-lane pressure FVDM(PLPFVDM) reduces error by 12% .Finally, the improved model is used to simulate trafc oscillation in the diversion area.Results show that lane changing behavior can trigger traffic oscilations.Higher driving speeds attenuate oscillation amplitudes,while lower speeds result in a concave growth patrn of speed standard deviation along the vehicle platoon in the upstream direction.Furthermore,lane changes occurring closer to the exit ramp exacerbate oscillations, whereas shorter lane-changing duration help suppress them.

    Keywords: transportation engineering; car-following model; traffic oscilation; lane-changing pressure; freeway diversionarea;numerical simulation

    “交通震蕩\"是一種車輛“走走停停\"的現(xiàn)象[2],對交通流的運行有嚴重的負面影響,如:燃料消耗加劇、有害氣體排放量大、安全風險劇增、增加行車延誤等。換道行為是高速公路中非常普遍的現(xiàn)象,是引發(fā)交通震蕩的一個重要誘因[3-4],具有非常高的研究價值,由于其過程的復(fù)雜性,目前的定量研究相對較少。

    換道行為研究主要是從駕駛員行為和道路交通狀況的角度研究城市和高速公路的變道特征和影響因素。Moridpour等[5描述了車輛換道的過程,主要包括換道動機的產(chǎn)生、換道時機的選擇及換道的執(zhí)行過程。根據(jù)研究動機,可將變道行為分為強制性變道和任意性變道。張穎達等利用Next Generation Simulation(NGSIM)數(shù)據(jù)集中的軌跡數(shù)據(jù),提出半對數(shù)和多項式模型來研究強制型和自由型換道行為的換道時間差異。換道行為中換道位置距離交叉口的距離、行駛速度和換道持續(xù)時間對交通流的穩(wěn)定性有一定的影響8-]。Toledo等[12分析了換道時間的分布范圍,發(fā)現(xiàn)了對換道時間產(chǎn)生影響的最主要因素是具有換道行為車輛與當前車道前后相鄰車輛的速度差值和與目標車道前后相鄰車輛的間距值。

    影響交通震蕩的因素包括縱向相互作用的不穩(wěn)定性、跟馳和換道行為、加減速不對稱特性、駕駛員類型以及道路線形等[13-14]。高速公路分流區(qū)域頻發(fā)的換道行為比交通流縱向不穩(wěn)定性更易誘發(fā)交通震蕩[15]。Mauch等通過研究高速公路交通流軌跡數(shù)據(jù)發(fā)現(xiàn)換道行為對交通震蕩的生成和發(fā)展有著非常大的影響。張娟采用累積移動序列法研究上海市快速路匝道交通震蕩的特性發(fā)現(xiàn),震蕩傳播超越駛?cè)朐训篮?,振幅被削弱,穿過駛出匝道后被放大。黃永現(xiàn)[8為探究交通震蕩的發(fā)展特點,進行了3次實車實驗,大致計算出車隊臨界速度的區(qū)間為 30~40km/h ,并將其用來度量交通流的穩(wěn)定性。Yao等利用現(xiàn)場試驗發(fā)現(xiàn),2輛連續(xù)車輛之間的時間間隔與前車的速度標準差和后車的初始速度呈負相關(guān),且與前車的平均速度和初始間距呈正相關(guān)。這些發(fā)現(xiàn)為更深入的理解交通震蕩現(xiàn)象提供了基礎(chǔ)。

    1 數(shù)據(jù)來源與篩選

    1.1 數(shù)據(jù)來源

    NGSIM數(shù)據(jù)集是對攝像機獲得的視頻數(shù)據(jù)應(yīng)用技術(shù)手段提取而得到的,收集了詳細、高質(zhì)量的交通流數(shù)據(jù),精度較高,多次被研究者們應(yīng)用于微觀駕駛行為分析、模型參數(shù)標定等研究。NGSIM數(shù)據(jù)集包括Peachtree、Lankershim Boulevard、I-80和US-101共4條路段的數(shù)據(jù),其中,前面2條是城市主干道,后面2條是高速公路。

    為研究高速公路分流區(qū)的交通震蕩現(xiàn)象,選取的數(shù)據(jù)采集區(qū)域位于加利福尼亞州的US-101高速公路路段,道路渠化信息如圖1所示。該區(qū)域道路長約2100ft(約 640.08m ),共有8條車道(主路車道5條,減速車道1條,匝道2條),每條車道的寬度為 3.66m ,最高限速 112.65km/h 。數(shù)據(jù)采集于2005年6月15日早晨07:50~08:05、08:05~08:20和08:20~08:35共3個時間段,采集對象為固定時間段和固定路段上的全部車輛,數(shù)據(jù)采集間隔0.1s,即1s產(chǎn)生10條數(shù)據(jù)。

    研究選取分流區(qū)內(nèi)換道車輛的駛出路線為從主路最右側(cè)車道5換道到減速車道6最后再到出口匝道車道8,換道駛出車輛必須在車道5進行強制換道,主要以車道5、車道6和車道8的車輛數(shù)據(jù)為基礎(chǔ)進行研究。

    Fig.1 Schematic diagram ofUS-101collection section

    1.2 數(shù)據(jù)處理

    1.2.1 數(shù)據(jù)篩選

    選擇的數(shù)據(jù)處理工具為python編程語言,對數(shù)據(jù)進行初步預(yù)處理,具體步驟及原則如下。

    1)以從車道5換道駛?cè)胲嚨?時車輛跨越車道線的時刻為依據(jù),提取此時刻前后15s(150條)的數(shù)據(jù)作為換道車輛行為分析的原始數(shù)據(jù)。

    2)提取具有換道行為的車輛軌跡數(shù)據(jù)。車輛的換道行為必須從車道5開始,換道至車道6,最后從出口匝道車道8駛出。剔除既沒有換道行為,也沒有對車輛的換道行為產(chǎn)生影響的車輛。

    3)剔除發(fā)生連續(xù)多次換道和換道失敗等行為的數(shù)據(jù),提取從車道5至車道6時只發(fā)生一次換道行為的車輛數(shù)據(jù)。

    4)換道開始時須保證車道6的上游匝道無車輛,否則將會對換道車輛的換道行為產(chǎn)生影響,進而會對研究分流區(qū)交通震蕩現(xiàn)象產(chǎn)生影響。

    經(jīng)上述篩選后,共提取出符合研究需求的換道車輛92輛。為研究換道駛出車輛的換道位置距出口匝道鼻端的距離和換道時間長度對分流區(qū)交通震蕩的影響,需提取換道車輛的換道時間和換道位置進行分析,如圖2所示。由圖2可知,大部分車輛的換道行為都在3~12 s完成,只有極少數(shù)車輛的換道時間大于12s,同時可以發(fā)現(xiàn)當換道起始點距出口匝道鼻端較近時,駕駛員完成換道行為的時間較短;當距鼻端較遠時,換道時間較長。大部分駕駛員的換道時間大于4s,只有極少數(shù)的駕駛員能夠在較短時間內(nèi)完成換道。因此,可以得出,大多數(shù)駕駛員會盡可能地選擇距出口匝道鼻端較遠的位置進行換道。綜上,出口匝道換道可能會對分流區(qū)內(nèi)駕駛員的心理產(chǎn)生影響,與普通路段的換道行為有明顯的差異。

    圖1 US-101采集路段示意圖

    圖2車輛換道時間-位置圖

    Fig.2Time-locationdiagramofvehicle lanechange

    1.2.2 數(shù)據(jù)修正

    NGSIM數(shù)據(jù)通過視頻拍攝的方式獲取,在拍攝過程中存在被前車遮擋、光線強度變化、周圍環(huán)境顏色突變以及圖像識別不準確(黑色車和陰影、白色車輛和標志標線等現(xiàn)象無法區(qū)分)等現(xiàn)象,造成數(shù)據(jù)中有部分缺失。此外,其速度值是由縱向位置數(shù)據(jù)求一階倒數(shù)獲得,存在突變誤差,研究采取一階差分修正數(shù)據(jù)中的異常值。

    根據(jù)已有研究,滿足車輛性能和人體所能承受極限的加速度值為 -8~5m/s2[20] ,同時,軌跡數(shù)據(jù)集中位移、速度和加速度相鄰幀之間的數(shù)據(jù)差值應(yīng)符合運動學的規(guī)律,相鄰幀的差值大約相等,即

    xt–xt-1≈xt-1-xt-2.

    式中: xt 為當前幀軌跡數(shù)據(jù); xt-1 為前1幀軌跡數(shù)據(jù); xt-2 為前2幀軌跡數(shù)據(jù)。

    利用前1幀和前2幀的數(shù)據(jù)差值來估算當前幀的軌跡數(shù)據(jù)值,計算方式如式(2)所示。

    xt*≈xt-1+(xt-1-xt-2),

    式中, xt* 為當前幀軌跡數(shù)據(jù)的預(yù)測值。

    根據(jù)確定的加速度范圍確定誤差范圍 Δx ,當預(yù)測值和采樣值的差值超出誤差范圍 Δx 時,確定當前時間幀的軌跡數(shù)據(jù)采樣值為異常值,然后將 xt* 重新替換為 x 要規(guī)避異常值為預(yù)測起始點的情況,即要滿足式(3)。

    當連續(xù)出現(xiàn)3個異常值時,應(yīng)重新選取起始點進行預(yù)測。以104號車的速度數(shù)據(jù)為例演示異常值修正效果,如圖3所示。

    圖3104號車速度修正效果圖

    Fig.3Speed correction effectofcar104

    圖3中,藍色線條表示的是車輛的原始軌跡數(shù)據(jù),紅色線條為車輛修正突變誤差之后的軌跡數(shù)據(jù)。從圖中可以看到被黃色圓圈標注的異常值已經(jīng)被剔除并替換,正常值未發(fā)生變化,修復(fù)后的速度數(shù)據(jù)與前后速度變化趨勢一致。數(shù)據(jù)從視頻中提取和誤差修正時,技術(shù)處理會帶來一定的誤差,造成與真實數(shù)據(jù)差距較大,因此,還需將測量誤差和數(shù)據(jù)處理誤差應(yīng)用濾波方法進行去噪處理。

    2考慮換道壓力增益的跟馳模型

    2.1 換道壓力增益因子構(gòu)建

    換道壓力[211是指外部因素對車輛換道行為的影響程度,與離最遲換道位置的距離有關(guān),越靠近最遲換道位置壓力越大,車輛換道的概率越大。

    為研究高速公路分流區(qū)出口匝道對換道行為的影響,提取US-101路段中92輛車換道起始點距出口匝道鼻端的距離。繪制換道起始點距出口匝道鼻端距離 L 的分布頻率圖,如圖4所示。經(jīng)擬合發(fā)現(xiàn)換道起始點距出口匝道鼻端距離 L 的頻率分布圖符合正態(tài)分布曲線,取置信區(qū)間為 95% ,可得出 L 的分布區(qū)間為 120~344m 。定義換道壓力增益因子 K(Ln-1) 衡量換道壓力對駕駛員心理的影響,即

    式中: μ 為換道壓力敏感系數(shù),取值范圍為 [0,1];Ln-1 為換道車輛換道位置距出口匝道鼻端的距離; Lmax 和 Lmin 為換道車輛換道位置距出口匝道鼻端的最遠距離和最近距離。如圖5所示, Lmax 為 344m,Lmin 為 120m ,后文將 K(Ln-1) 簡化為 K(?) 。

    圖4換道起始點距離出口匝道鼻端距離分布頻率直方圖

    Fig.4 Frequency of the distance from the start of the lane change to the nose end of the exit ramp

    圖5換道車輛距出口匝道鼻端距離

    Fig.5Distance from the nose of the exit ramp for lane changevehicles

    2.2 改進跟馳模型

    全速度差模型(fullvelocity difference model,F(xiàn)VDM)[22]是在廣義力模型的基礎(chǔ)上,考慮了前后車正、負速度差對駕駛員的影響所建立的應(yīng)用廣泛的跟馳模型之一。換道車輛到達出口匝道影響范圍之前車輛的行駛采用FVDM,根據(jù)前后2輛車的速度差和間距差調(diào)整自身加速度,改變自身的駕駛狀態(tài)。

    當前車行駛到出口匝道影響范圍時,可能發(fā)生換道駛出行為,會刺激后車集中注意力,有意識的去調(diào)整駕駛狀態(tài),跟馳車輛的駕駛行為受前車和后車的間距差 Δxn,n-1(t), 第 n 輛車在 t 時刻與前車的速度差Δνn,n-1(t) 及 Ln-1 的影響。結(jié)合式FVDM的基本形式,建立考慮換道壓力增益的全速度差跟馳模型(pressure-full velocity difference model,P-FVDM)

    α=1/τ,

    κ=λ/τ

    速度優(yōu)化函數(shù)可表示為

    式中: 為跟馳車輛 n 雞在難 t 時刻的加速度, m/s2;a 為間距差敏感系數(shù);τ為駕駛員反應(yīng)時間, s;νn(t) 為第 n 輛車在 t 時刻的速度, m/s;κ 為速度差敏感系數(shù);λ為敏感系數(shù)常量; lc 為P-FVDM中車輛的最小安全間距, m ;C1 和 C2 為常數(shù)項。

    前車行駛一定距離后,開始發(fā)生換道駛出行為,此時轉(zhuǎn)向燈亮起,開始發(fā)生橫向偏移,如圖5所示。此時跟馳車輛的駕駛行為受 橫向偏移 ΔWn,n+1 及后車距出口匝道鼻端的距離 Ln+1 的影響。應(yīng)用雙車道跟馳模型理論,在FVDM的基礎(chǔ)之上,考慮前車駛離當前車道時產(chǎn)生的橫向偏移對后車的跟馳行為的影響,王玄金[23]建立換道全速度差模型(lane pass-fullvelocity diference model,LP-FVDM)。結(jié)合LP-FVDM的基本形式,建立考慮換道壓力增益和換道駛出間距增益的全速度差模型(pressure lane pass-full velocitydifference model,PLP-FVDM)

    速度優(yōu)化函數(shù)可表示為

    式中: ρ 為橫向間距敏感系數(shù); G(?) 為橫向間距增益因子; ΔWmax 為車道寬度, m 。

    圖6橫向偏移示意圖

    Fig.6Schematicdiagramoflateraloffset

    2.3 線性穩(wěn)定性分析

    假設(shè)在初始狀態(tài)下,車隊的車頭間距都相同為 p,m ,車輛運行速度為優(yōu)化速度 V((1+G(?)+K(?))p) 則得到各車輛的初始位置為

    假設(shè)車輛 n 在 t 時刻由于小的擾動,產(chǎn)生偏離,定義為 yn(t) ,有

    xn({})=xn0(t)+yn(t)

    將式(11)和式(12)帶人到式(10)中,可得

    式中

    Δyn(t)=yn-1(t)-yn(t)

    將式(17)帶人到式(15)中,并將其進行展開,得到關(guān)于 z 的方程

    z2=α[(1+G(?)+K(?))V((1+G(?)+K(?))(e-1)-z]+κ(e-1)z0.

    令:

    z=z1(iκ)+z2(iκ)2+…

    將式(19)帶入式(18)中,計算得 iκ 的一階、二階項的系數(shù)分別為

    當 z2gt;0 時,均勻交通流會逐漸演變?yōu)榉€(wěn)定狀態(tài);當 z2lt;0 時,均勻交通流會逐漸演變?yōu)椴环€(wěn)定狀態(tài)。當z2=0 時擾動衰減,此時可得模型的中性穩(wěn)定條件為

    為更好的分析 K(?) 對模型穩(wěn)定性的影響,分別令 λ=0.5,G(?)=0.3 做出壓力增益因子取不同值時的臨界穩(wěn)定曲線圖,如圖7所示。

    圖7 λ=0.5,G(?)=0.3 的穩(wěn)定性臨界圖

    Fig.7Stabilitycriticalitydiagramfor λ=0.5 and δG(?)=0.3

    圖7中,隨著 K(?) 的增大,臨界曲線向左上方移動,反應(yīng)系數(shù)的最大值升高,穩(wěn)定區(qū)域范圍變小,此時交通流能夠在相對較小的車頭間距下穩(wěn)定運行,但要求駕駛員的反應(yīng)時間足夠小。在臨界曲線附近更容易出現(xiàn)交通流不穩(wěn)定的現(xiàn)象,表明出口匝道對交通流的運行有較大影響,在分流區(qū)內(nèi)交通流的穩(wěn)定性下降,換道駛出的車輛對主路車流的影響較大,極易引發(fā)交通震蕩不斷向上游傳播。

    當 G(?)=0 時,式(19)退化為P-FVDM的穩(wěn)定性條件,同樣可以得出隨著 K(?) 值的增大,模型的穩(wěn)定區(qū)域變小,對駕駛員反應(yīng)能力的要求更高,更易引發(fā)交通震蕩。

    3改進跟馳模型精度驗證

    3.1模型參數(shù)標定

    利用前文得到的換道車輛和其跟馳車輛的軌跡數(shù)據(jù)對改進模型的參數(shù)進行標定,結(jié)合基礎(chǔ)模型的參數(shù)值確定標定參數(shù)的取值范圍,見表1所示。標定方法選取遺傳算法,遺傳算法的參數(shù)設(shè)置見表2。

    表1參數(shù)取值范圍Table1 Parametersrangeofvalues

    表2遺傳算法參數(shù)

    Table2 Parametersof thegeneticalgorithm

    研究選取的性能指標為跟馳車輛的加速度,使用RMSE作為擬合優(yōu)度函數(shù)

    式中: Xobs,i 為實際值; Xmodel,i 為模型計算值; N 為樣本量。

    經(jīng)GA參數(shù)標定后,P-FVDM和PLP-FVDM的模型參數(shù)的均值、中值、最大值、最小值和標準差分別如表3和表4所示。

    Table3 CalibrationresultsoftheP-FVDMmodel

    表4PLP-FVDM標定結(jié)果

    表3P-FVDM標定結(jié)果

    Table4 CalibrationresultsofPLP-FVDMmodel

    3.2 模型仿真

    為了驗證P-FVDM和PLP-FVDM的效果,選取車輛ID為13、1226共2輛車的軌跡數(shù)據(jù),利用P-FVDM和PLP-FVDM對其進行速度和位置的仿真,與實際軌跡數(shù)據(jù)對比分析。用原始軌跡中的初始車輛速度和位置,利用模型計算出加速度,由此可以計算出車輛的仿真速度和位置。速度仿真效果和位置仿真效果分別如圖8和9所示。

    圖8車輛速度仿真效果圖Fig.8Vehicle speed simulation

    由圖8可知,P-FVDM與PLP-FVDM對車輛速度的模擬效果良好,其中1226號車的速度仿真誤差在 1m/s 范圍內(nèi),13號車的速度仿真誤差在 1.5m/s 范圍內(nèi)。綜上,P-FVDM和PLP-FVDM速度仿真誤差都在 1.5m/s 以內(nèi),誤差范圍較小,雖然當車輛速度突變的程度較大時仿真效果相對較差,但總體與車輛的速度變化趨勢是相符的,能夠再現(xiàn)分流區(qū)真實交通流的速度信息。

    圖9車輛位置仿真效果圖

    Fig.9 Vehiclepositionsimulation

    由圖9可知,1226號車位置仿真誤差范圍在 3m 內(nèi);13號車的位置仿真誤差范圍在 4m 內(nèi)。因此,模型對車輛位置的仿真效果總體較好,總體位置仿真誤差范圍在 4m 內(nèi),P-FVDM和PLP-FVDM能夠較好的預(yù)測車輛的位置,能夠再現(xiàn)真實交通流中車輛的位置信息,進一步驗證了模型的準確性。

    3.3 模型有效性驗證

    為檢驗?zāi)P偷挠行裕瑢Ρ萈-FVDM和FVDM、PLP-FVDM和LP-FVDM的速度仿真效果,分別計算13、48、100、359、703和1226號車輛改進前和改進后模型的仿真速度的均方根百分比誤差值,即RMSPE,計算公式為式(24),結(jié)果分別如圖10和11所示。

    由圖10可知,與FVDM相比,在換道車輛進入到出口匝道影響區(qū)域時,P-FVDM所控制的跟馳車輛的仿真誤差比FVDM的誤差減小,其中703號車的仿真誤差減小的最多(約為 32% ),359號車的仿真誤差減小的最少(約為 8% ),模型的總體提升約 16% 。因此,P-FVDM能夠更加精準的反映換道車輛進入出口匝道影響區(qū)發(fā)生換道行為之前主路跟馳車輛的駕駛情況。由圖11可知,與LP-FVDM相比,在換道車輛發(fā)生換道行為過程中,PLP-FVDM所控制的跟馳車輛的速度仿真誤差與真實值的誤差有所減小,模型的精度總體提升約12% ,能夠更加精準的反映換道車輛換道過程中主路跟馳車輛的駕駛情況。

    圖10P-FVDM效果對比

    Fig.10 Comparison of the effect of P-FVDM model

    圖11PLP-FVDM效果對比 Fig.11 Comparison of the effect of PLP-FVDM model

    綜上,提出的模型對分流區(qū)交通流的模擬效果更好。

    4交通震蕩仿真研究

    4.1 仿真條件

    采用城市交通模擬系統(tǒng)(simulation of urban mobility,SUMO)[24]仿真軟件模擬換道車輛駛出時的分流區(qū)內(nèi)交通震蕩,初始狀態(tài)設(shè)定換道車輛后面由20輛車組成的車隊跟馳行駛,所有車輛之間的跟馳模型為FVDM,給定初始速度和初始間距使車隊行駛至穩(wěn)定狀態(tài)。

    t=0 s時,換道駛出車輛到達出口匝道影響區(qū)域的起點,換道駛出車輛和車隊第一輛車之間的控制模型改為P-FVDM,車隊內(nèi)部的車輛控制模型始終為FVDM,保持不變。在 t=21 s時,換道駛出車輛開始換道,發(fā)生橫向偏移,換道時間為 6s 。此時換道駛出車輛和車隊第一輛車之間的控制模型改變?yōu)镻LP-FVDM。換道車輛沿道路縱向的行駛速度不變,橫向偏移的方式由車輛的換道時間和車道寬度( 3.66m 計算得知,在前2.5s內(nèi)車輛勻加速,速度達到最大時車輛橫向位置到達主路車道和減速車道的車道線,在后 2.5s 內(nèi)車輛橫向速度開始均勻減速,當橫向速度減為 0m/s 時,車輛到達減速車道的中心線,此時完成換道行為。在 時,換道駛出車輛完成換道,此時對車隊的第一輛車不產(chǎn)生影響,車隊第一輛車以當前速度繼續(xù)行駛。

    選取速度序列的標準差作為衡量震蕩振幅的指標

    式中: u(t) 為當前車輛的速度采集值; um 為速度序列的平均值; n 為一輛車的仿真數(shù)據(jù)點數(shù); In 為仿真的時間序列。

    4.2交通流速度對交通震蕩的影響

    為研究不同交通流速度對交通震蕩的影響,結(jié)合高速公路出口匝道附近車輛運行特性,設(shè)定換道駛出車輛和車隊的運行速度為 30,40,50km/h ,初始車頭間距由FVDM計算得出,如表5所示。為了使交通震蕩現(xiàn)象更加便于觀察,繪制出了車隊中車輛的速度標準差分布圖,仿真結(jié)果如圖12所示。

    圖12車隊速度標準差圖

    Fig.12Diagramof standarddeviationofvehiclefleetspeed

    由圖12可知,車隊以 30km/h 勻速行駛時,車流內(nèi)部產(chǎn)生擾動,隨著震蕩沿車隊向上游傳播,車隊中車輛的速度標準差呈現(xiàn)凹增長的模式。位于車隊上游的車輛速度標準差較大,即速度變化較大,震蕩振幅增大,震蕩現(xiàn)象更明顯,且隨著傳播過程快速增長。隨著震蕩不斷向上游傳播,車輛的速度標準差的增長幅度逐漸減緩,甚至呈現(xiàn)降低趨勢,即此時震蕩的增長開始趨于飽和。當車隊以 40km/h 勻速行駛時,所有車輛的速度標準差非常小,均小于0.5,且沿車隊向上游傳播過程中迅速減小,即交通震蕩在傳播過程中逐漸消散。當車隊以 50km/h 勻速行駛時,車輛速度標準差更小,均小于0.15,且降低速度更快,即交通震蕩消散的更快。

    車隊以 30km/h 勻速行駛時,交通流內(nèi)部不穩(wěn)定和換道車輛的換道行為共同作用,產(chǎn)生交通震蕩。車隊分別以 40,50km/h 勻速行駛時,換道車輛的換道行為引起了交通震蕩,但此時車輛的車頭間距較大,前方出現(xiàn)擾動時后方駕駛員有更多的時間對車輛進行調(diào)控,因此,交通震蕩在向上游傳播的過程中逐漸消散。

    綜上所述,車輛換道行為能夠引發(fā)交通震蕩,但是當行駛速度較高時震蕩會衰減,當行駛速度較低時速度標準差沿車隊向上游呈凹增長趨勢增長。

    4.3換道位置對交通震蕩的影響

    初始條件下,讓換道駛出車輛和跟馳車隊均以 50km/h 的速度勻速行駛,車頭間距為 28.5m ,所有車輛均以FVDM跟馳行駛,在到達出口匝道影響區(qū)域之前使交通流達到穩(wěn)定狀態(tài)。 t=0s 時,換道駛出車輛到達匝道影響區(qū)域的起點,設(shè)定換道駛出車輛的加速度的值為 -0.2m/s2 ,分別在 s和 t=25 s開始換道,換道時間均為6s,其余仿真設(shè)置與4.2節(jié)設(shè)置相同。

    車隊車輛的速度標準差仿真結(jié)果如圖13所示。換道時刻越晚,表示換道車輛的換道起始點的位置距出口匝道鼻端的距離越近。由圖13可以看出,換道車輛在不同的位置換道駛出時,均發(fā)生了明顯的交通震蕩現(xiàn)象,并向車隊上游不斷傳播。分析發(fā)現(xiàn),換道車輛在不同位置換道駛出時,車輛速度標準差都以凹增長的方式增加,與以往研究相符。隨著換道起始點距出口匝道鼻端距離減小,換道車輛的換道行為對后方跟馳車輛形成的影響越大,同一輛車的速度標準差增大,沿車隊向上游傳播的速度標準差的增加幅度也隨之增大。即換道車輛越早發(fā)生換道駛出行為對交通流的穩(wěn)定越有利。

    圖13車隊速度標準差圖

    Fig.13Diagramofstandarddeviationofvehiclefleetspeed

    圖14車隊速度標準差圖

    Fig.14Diagramof standard deviationofvehicle fleetspeed

    4.4換道時長對交通震蕩的影響分析

    初始條件下,讓換道駛出車輛和跟馳車隊均以 50km/h 的速度勻速行駛,車頭間距為 28.5m ,所有車輛均以FVDM跟馳,達到穩(wěn)定狀態(tài)。

    時,換道駛出車輛到達出口匝道影響區(qū)域的起點,設(shè)定換道駛出車輛的加速度的值為 -0.2m/s2 ,在 時開始換道,產(chǎn)生橫向位移,分別在 s和 時完成換道,即換道時間分別為5s,7s和9s,其余仿真設(shè)置與文中設(shè)置相同。車隊中車輛的速度標準差仿真結(jié)果如圖14所示。

    從圖14中可以發(fā)現(xiàn),換道時間越長,換道行為對后方車輛的影響越大,速度標準差沿車隊向上游傳播增長的越快,同時單個車輛的速度標準差更大,震蕩更加明顯。因此,換道時間越短對交通流的穩(wěn)定性越有利。

    5結(jié)束語

    以NGSIM數(shù)據(jù)集的US-101路段數(shù)據(jù)為數(shù)據(jù)源,對數(shù)據(jù)處理后發(fā)現(xiàn),分流區(qū)內(nèi)換道駛出車輛的換道時間范圍為3~12s,換道位置距出口匝道鼻端的距離為 50~350m ,換道位置和換道時間之間具有一定的關(guān)系。

    通過構(gòu)建高速公路分流區(qū)換道壓力增益因子對模型進行改進,并進行穩(wěn)定性分析,發(fā)現(xiàn)換道壓力敏感系數(shù)越高,模型的穩(wěn)定區(qū)域越小。利用重構(gòu)后的軌跡數(shù)據(jù)對模型進行參數(shù)標定,仿真分析發(fā)現(xiàn),改進后的模型能夠更好地模擬車輛的速度和位置,P-FVDM的仿真精度提升了 16% ,PLP-FVDM的仿真精度提升了 12% 。利用改進后的模型對高速公路分流區(qū)內(nèi)的交通震蕩進行仿真,研究結(jié)果表明,車速為 30km/h 時產(chǎn)生明顯的交通震蕩現(xiàn)象,并呈凹增長趨勢,車速為40和 50km/h 時交通震蕩在向后傳播過程中逐漸消散;隨著換道位置距離出口匝道鼻端距離變小、換道時間變長,速度標準差更大,即交通震蕩振幅更大,交通流的穩(wěn)定性更差。

    研究采用的數(shù)據(jù)集來自美國,不同國家的交通環(huán)境和人文因素相差較大,因此,交通流也存在一定差異,提出的模型仍需要中國的交通流數(shù)據(jù)進行進一步的驗證。再者,研究只考慮了分流區(qū)內(nèi)一輛車換道駛出對主路交通流的影響,后續(xù)可考慮多輛車連續(xù)或間斷駛出時對主路交通流的影響。

    參考文獻

    1」Euie LU.Car-onowmg anu steauy-state meory ior noncongested trac[J].Uperauons Kesearcn,1yol,y(1): o-/0.

    [2]YeoH,SkabardonisA.Understandingstop-and-gotraficinviewofasymmetrictraffic theoryM]/TransportationandTraffic Theory 2009:Golden Jubilee.Boston,MA: Springer US,2009:99-115.

    [3]KernerBS,Rehborn H.Experimentalfeaturesandcharacteristicsof traicjams[J].Physical ReviewE,StatisticalPhysics, Plasmas,F(xiàn)luids,and Related Interdisciplinary Topics,1996,53(2): R1297-R1300.

    [4]Ahn S,CassidyMJ.Freeway traicoscilationsandvehiclelane-change maneuversC]//rocedingsofthe7thIntetioal SymposiumonTraficandTransportationTheory.July23-25,20o7,London,UK.Amsterdam:Elsevier,20o7:691-710.

    [5]Moridpour S,Rose G,SarviM.Modelling the heavy vehicle drivers’lanechanging decisionunder heavy traficconditios[J]. Road amp; Transport Research,2009,18(4):49-57.

    [6]張穎達,邵春福,李慧軒,等.基于NGSIM軌跡數(shù)據(jù)的換道行為微觀特性分析[J].交通信息與安全,2015,33(6):19-24,32. Zhang YD,ShaoCF,LiHX,et al. Microscopic characteristicsof lane-change maneuvers based onNGSIM[J].Joualof Transport Information and Safety, 2015,33(6): 19-24,32.(in Chinese)

    [7]U.S.DepartmentofTransportation,F(xiàn)ederalHighwayAdministration(FHWA).Nextgenerationsimulationfactsheet[B/OL]. (2011-05-27)[2023-12-06]. htp://ps.fhwa.dot.gov/trafficanalysistools/ngsim.htm.

    [8]SusilawatiS,TanKY,SamadKMA,etal.Thestudyoflocationspecificlanechangeimpactsotraficdelayusingetended cell transmissonmodel[C]/Proceedingsof the2018EUREKA International ConferenceonEnginering,Technologyand Construction (EUREKA2018).April24-26,2018,Kuching,Sarawak,Malaysia.Les Ulis:EDP Sciences,2018,203:05008. DOI: 10.1051/matecconf/201820305008.

    [9]LiXP,CuiJX,AnS,etal.Stop-and-gotrafcanalysis:theoreticalproperties,environmentalimpactsandoilation mitigation[J].Transportation Research Part B:Methodological, 2014,70: 319-339.

    [10]鄭淑欣,李振龍,趙曉華.霧天高速公路換道行為特性研究[J].交通信息與安全,2020,38(1):35-42,51. Zheng SX,LiZL,Zhao XH.Behavioral characteristicsof lane-changingon highwayin foggyweather[J].Journalof Transport Information and Safety, 2020,38(1): 35-42, 51.(in Chinese)

    [11]Zheng ST,JiangR,TianJF,etal.Empiricalandexperimentalsudyontegowthpaternof taicoscilationsupstramof fixed bottleneck and model test[J].Transportation Research Part C:Emerging Technologies,2022,140:103729.

    [12]ToledoT,ZoharD.Modelingdurationof lanechanges[J].Transportation Research Record:Journalof theTransportation Research Board,2007,1999(1): 71-78.

    [13]Jiang R,Hu MB,Zhang HM,et al.Onsome experimental featuresof car-folowing behaviorand how to model them[J]. Transportation Research Part B:Methodological, 2015,80:338-354.

    [14]JiangR,JinCJ,ZhangHM,etal.ExperimentalandmpiricalinvestigatiooftraicfloistabilityCProeedingoft 6th International Symposium of the Transport and Trafic ScienceForum (HEART2016).May17-20,2016,Xi'an,China. Amsterdam: Elsevier,2017,23:157-173.DOI: 10.1016/j.trpro.2017.05.011.

    [15]賀玉龍,劉磊,遲佳欣.高速公路車輛換道行為風險研究[J].重慶交通大學學報(自然科學版),2021,40(4):26-33. HeYL,LiuL,ChiJX.Risk of vehiclelanechange behavioronhighway[J].JournalofChongqing Jiaotong University (Natural Science),2021,40(4): 26-33.(in Chinese)

    [16]Mauch ,CassidyM.Freeway traffc:oscilltions observations and predictionsM]/Transportationand TraffcTheorynthe 21stCentury.Amsterdam:Elsevier,2002: 653-673.

    [17]張娟.城市快速路瓶頸交通震蕩特性研究[J].交通信息與安全,2015,33(5):43-48. ZhangJ.Astudyoftraffcoscillationsat thebottlenecksonurbanexpressways[J].JouralofTransportInformationandSafety 2015,33(5): 43-48.(in Chinese)

    [18]黃永現(xiàn).交通震蕩演化特性的實驗研究和建模仿真[D].合肥:中國科學技術(shù)大學,2019. Huang YX.Experimental research and modelingontheevolution of trafic oscilation[D]. Hefei: UniversityofScience and Technology of China,2019.(in Chinese)

    [19]YaoHD,LiQWLiXP.Astudyofrelationshipsintraicoscilltonfeaturesbasdonfieldexperiments[J].Trasprtation ResearchPartA:Policyand Practice,2020,141:339-355.

    [20]馬艷鋒.基于NGSIM數(shù)據(jù)的高速公路出入口強制換道車輛執(zhí)行及后調(diào)整行為研究[D].西安:長安大學,2020. MaYF.Research onvehicles’mandatorylane change behavior in executionand post-adjustment periodonmulti-lane freeway based on NGSIM trajectory data[D]. Xi'an; Chang'an University,2020.(in Chinese)

    [21]劉有軍,曹珊.基于元胞自動機的強制換道模型研究[J].交通信息與安全,2009,27(3):78-80. LiuYJ,Cao S.Compulsorylane-changing traficmodel based on celularautomaton[J].Journalof Tansport Informationand Safety,2009,27(3):78-80.(in Chinese)

    [22]JiangR,WuQS,ZhuZJ.Fullvelocitydierence modelforacar-following heory[J].PhysicalReviewE,2001,64017101.

    [23]王玄金.考慮前車換道過程的跟馳模型研究[D].重慶:重慶大學,2017. Wang X J.Research on car following model considering preceding vehicle’slane-changing processD].Chongqing: Chongqing University,2017. (in Chinese)

    [24]Barcel6J.FundamentalsoftraficsimulationM].NewYork:Springer,20l0.(InternationalSeriesiOperations Researchamp; Management Science,vol 145).

    (編輯 呂建斌)

    猜你喜歡
    交通流匝道標準差
    基于車聯(lián)網(wǎng)的中心導(dǎo)航云性能分析模型
    中國綠色化發(fā)展水平的時空格局及優(yōu)化策略
    智慧收費系統(tǒng)在全椒西收費站的應(yīng)用
    袋栽鹿茸菇菌棒物理力學特性研究
    考慮換道壓力增益的高速公路分流區(qū)跟馳模型
    芻議教材“二次開發(fā)”的若干視角
    亚洲国产最新在线播放| 国产精品福利在线免费观看| www.av在线官网国产| 夜夜爽夜夜爽视频| 九九爱精品视频在线观看| 丝袜在线中文字幕| 91精品伊人久久大香线蕉| 久久久久久久久久久免费av| 水蜜桃什么品种好| 18禁在线无遮挡免费观看视频| 成人特级av手机在线观看| 亚洲国产精品专区欧美| 亚洲精品日韩在线中文字幕| 亚洲成色77777| 看免费成人av毛片| av视频免费观看在线观看| 十八禁高潮呻吟视频 | 国产日韩欧美视频二区| 色婷婷久久久亚洲欧美| 99九九在线精品视频 | 亚洲精品一二三| 国产精品国产av在线观看| 成年女人在线观看亚洲视频| 日韩一区二区三区影片| 亚洲国产色片| 午夜激情福利司机影院| 日韩欧美精品免费久久| 精品少妇内射三级| av国产精品久久久久影院| 亚洲色图综合在线观看| 亚洲色图综合在线观看| 六月丁香七月| 日韩欧美精品免费久久| 一级毛片我不卡| 亚洲综合精品二区| 免费观看的影片在线观看| 嫩草影院入口| 最后的刺客免费高清国语| 王馨瑶露胸无遮挡在线观看| 一二三四中文在线观看免费高清| 热re99久久国产66热| 日本爱情动作片www.在线观看| 精品少妇久久久久久888优播| 一级,二级,三级黄色视频| 亚洲av二区三区四区| 国产精品一区二区三区四区免费观看| 国产伦在线观看视频一区| 日韩电影二区| 亚洲丝袜综合中文字幕| 久久精品熟女亚洲av麻豆精品| 色视频www国产| 三级国产精品片| 国产乱来视频区| 久久精品久久久久久久性| 午夜视频国产福利| 亚洲精品国产av成人精品| 久久久久久久精品精品| 国产日韩一区二区三区精品不卡 | 80岁老熟妇乱子伦牲交| 精品久久久噜噜| 人妻 亚洲 视频| 久久久久久久久久成人| 少妇的逼水好多| 午夜免费男女啪啪视频观看| 国产精品人妻久久久久久| 国产精品人妻久久久久久| 内射极品少妇av片p| 大香蕉久久网| 色94色欧美一区二区| 亚洲成人av在线免费| 欧美日韩综合久久久久久| 亚洲欧美一区二区三区黑人 | 久久免费观看电影| 亚洲精品自拍成人| 国产午夜精品久久久久久一区二区三区| 亚洲av中文av极速乱| 久久久午夜欧美精品| 国产极品天堂在线| 国产午夜精品久久久久久一区二区三区| av福利片在线| 少妇人妻久久综合中文| 一本色道久久久久久精品综合| 精品久久久久久电影网| 一本色道久久久久久精品综合| 亚洲欧洲日产国产| 在线观看美女被高潮喷水网站| 中国国产av一级| av卡一久久| 女性被躁到高潮视频| 亚洲精品第二区| 香蕉精品网在线| 高清在线视频一区二区三区| 在线观看一区二区三区激情| 日韩精品免费视频一区二区三区 | 久久99一区二区三区| 国产免费一级a男人的天堂| 涩涩av久久男人的天堂| 国产伦理片在线播放av一区| 久久久久久久久久成人| 一个人看视频在线观看www免费| 日本wwww免费看| 亚洲一区二区三区欧美精品| 午夜av观看不卡| 美女国产视频在线观看| 啦啦啦在线观看免费高清www| 午夜免费鲁丝| 国产成人免费无遮挡视频| 亚洲婷婷狠狠爱综合网| 日韩免费高清中文字幕av| 亚洲欧美日韩另类电影网站| 国产69精品久久久久777片| 国产无遮挡羞羞视频在线观看| 乱人伦中国视频| 免费少妇av软件| 18+在线观看网站| 十八禁网站网址无遮挡 | 少妇猛男粗大的猛烈进出视频| 亚洲国产最新在线播放| 一本大道久久a久久精品| 国产美女午夜福利| 极品教师在线视频| 极品少妇高潮喷水抽搐| 精品亚洲乱码少妇综合久久| a级一级毛片免费在线观看| 一个人看视频在线观看www免费| 久久久久久久久久人人人人人人| 9色porny在线观看| 国产精品秋霞免费鲁丝片| 国产精品人妻久久久影院| 国产视频首页在线观看| 久久 成人 亚洲| 夜夜骑夜夜射夜夜干| 免费在线观看成人毛片| 亚洲av在线观看美女高潮| 人妻 亚洲 视频| 看十八女毛片水多多多| av网站免费在线观看视频| 一区在线观看完整版| 日日撸夜夜添| 欧美国产精品一级二级三级 | 国产伦在线观看视频一区| 日本黄大片高清| 深夜a级毛片| av免费在线看不卡| 亚洲情色 制服丝袜| 午夜老司机福利剧场| 免费不卡的大黄色大毛片视频在线观看| 最黄视频免费看| 亚洲av.av天堂| 国产亚洲最大av| 毛片一级片免费看久久久久| 麻豆成人午夜福利视频| 超碰97精品在线观看| 免费黄频网站在线观看国产| 精品一品国产午夜福利视频| 五月开心婷婷网| 亚洲精品乱码久久久久久按摩| 亚洲欧美成人综合另类久久久| 97在线视频观看| 国产免费一区二区三区四区乱码| 国产中年淑女户外野战色| 丝瓜视频免费看黄片| 91精品一卡2卡3卡4卡| 肉色欧美久久久久久久蜜桃| 免费人妻精品一区二区三区视频| 爱豆传媒免费全集在线观看| 久久精品久久精品一区二区三区| 亚洲国产色片| 亚洲不卡免费看| 午夜视频国产福利| 国产男女超爽视频在线观看| 亚洲精品亚洲一区二区| 日韩制服骚丝袜av| 日韩一本色道免费dvd| 中文字幕人妻熟人妻熟丝袜美| 80岁老熟妇乱子伦牲交| 国产精品一区二区在线观看99| 男的添女的下面高潮视频| 亚洲美女视频黄频| 成人18禁高潮啪啪吃奶动态图 | 桃花免费在线播放| 少妇裸体淫交视频免费看高清| 久久国内精品自在自线图片| 超碰97精品在线观看| 在线观看免费视频网站a站| 夫妻性生交免费视频一级片| 国产极品天堂在线| 大码成人一级视频| 久久97久久精品| 少妇 在线观看| 国产成人精品一,二区| 99国产精品免费福利视频| 成人毛片60女人毛片免费| 一区二区三区精品91| 日韩三级伦理在线观看| 黄色怎么调成土黄色| 99久国产av精品国产电影| 狠狠精品人妻久久久久久综合| 国产午夜精品一二区理论片| 国产av码专区亚洲av| 亚洲精品日韩av片在线观看| 欧美日韩亚洲高清精品| 乱系列少妇在线播放| 老熟女久久久| 高清视频免费观看一区二区| 伊人亚洲综合成人网| 菩萨蛮人人尽说江南好唐韦庄| 噜噜噜噜噜久久久久久91| 日本欧美国产在线视频| 自线自在国产av| 久久毛片免费看一区二区三区| 插逼视频在线观看| 搡女人真爽免费视频火全软件| 在线观看美女被高潮喷水网站| 国产黄色视频一区二区在线观看| 国产亚洲av片在线观看秒播厂| 亚洲久久久国产精品| 国产 精品1| 天美传媒精品一区二区| 国产片特级美女逼逼视频| 国产91av在线免费观看| 亚洲国产精品999| 国产黄色免费在线视频| av又黄又爽大尺度在线免费看| 卡戴珊不雅视频在线播放| 99热6这里只有精品| 在线观看免费视频网站a站| 精品国产露脸久久av麻豆| 69精品国产乱码久久久| 久久影院123| 最近的中文字幕免费完整| 久久久国产欧美日韩av| 美女内射精品一级片tv| 丰满饥渴人妻一区二区三| 国产色爽女视频免费观看| 国产在线男女| 99视频精品全部免费 在线| 哪个播放器可以免费观看大片| 欧美少妇被猛烈插入视频| 日本欧美国产在线视频| 国产乱人偷精品视频| 亚洲久久久国产精品| 亚洲av国产av综合av卡| 2018国产大陆天天弄谢| 成人无遮挡网站| 欧美日韩亚洲高清精品| 欧美人与善性xxx| 婷婷色麻豆天堂久久| 色哟哟·www| 精品一区二区三区视频在线| 欧美97在线视频| 亚洲电影在线观看av| 亚洲av.av天堂| 精品久久久精品久久久| 丝瓜视频免费看黄片| 国产在线视频一区二区| 国产午夜精品久久久久久一区二区三区| 国产视频首页在线观看| 少妇 在线观看| 美女中出高潮动态图| 国产免费福利视频在线观看| 亚洲精品乱久久久久久| 18禁在线无遮挡免费观看视频| 亚洲av综合色区一区| 久久av网站| 亚洲精品日韩在线中文字幕| 国产亚洲欧美精品永久| 亚洲欧美清纯卡通| 97在线人人人人妻| 国产精品.久久久| 国产一区二区在线观看日韩| 偷拍熟女少妇极品色| 免费在线观看成人毛片| 久久毛片免费看一区二区三区| 日本黄大片高清| 丰满人妻一区二区三区视频av| av国产久精品久网站免费入址| 熟妇人妻不卡中文字幕| 国产白丝娇喘喷水9色精品| 搡女人真爽免费视频火全软件| 久久韩国三级中文字幕| 久久久久久久国产电影| 成人免费观看视频高清| 深夜a级毛片| 啦啦啦中文免费视频观看日本| 婷婷色av中文字幕| 国产精品99久久久久久久久| 免费观看a级毛片全部| 天堂8中文在线网| 国产无遮挡羞羞视频在线观看| 色5月婷婷丁香| 亚洲欧美精品专区久久| 日韩强制内射视频| 亚洲欧美日韩另类电影网站| 97超碰精品成人国产| 欧美日韩一区二区视频在线观看视频在线| 亚洲欧美成人综合另类久久久| 亚洲成人手机| 久久久久久久大尺度免费视频| 色94色欧美一区二区| 日韩中字成人| 99久久综合免费| 精品亚洲成a人片在线观看| 精品国产一区二区三区久久久樱花| 97超视频在线观看视频| 99久久精品国产国产毛片| 美女中出高潮动态图| 永久免费av网站大全| 丰满乱子伦码专区| 麻豆乱淫一区二区| 久久久久国产网址| 欧美bdsm另类| 人妻人人澡人人爽人人| 大又大粗又爽又黄少妇毛片口| 日韩人妻高清精品专区| 一级二级三级毛片免费看| 美女视频免费永久观看网站| 亚洲经典国产精华液单| 成人国产av品久久久| 国产一区二区在线观看日韩| 老司机亚洲免费影院| 熟女av电影| 国产伦精品一区二区三区视频9| 一二三四中文在线观看免费高清| 日本午夜av视频| 久久韩国三级中文字幕| 人人妻人人看人人澡| 亚洲综合色惰| 国产精品久久久久久av不卡| 97在线人人人人妻| 欧美亚洲 丝袜 人妻 在线| 午夜福利影视在线免费观看| 精品一品国产午夜福利视频| 菩萨蛮人人尽说江南好唐韦庄| 国产乱来视频区| 亚洲经典国产精华液单| 午夜免费鲁丝| 久久国产精品男人的天堂亚洲 | 欧美3d第一页| 国内少妇人妻偷人精品xxx网站| 一级毛片久久久久久久久女| 欧美一级a爱片免费观看看| 免费观看a级毛片全部| 简卡轻食公司| 亚洲欧美成人精品一区二区| 国内揄拍国产精品人妻在线| 三级国产精品片| 免费看光身美女| 有码 亚洲区| 99热这里只有是精品在线观看| 精品国产露脸久久av麻豆| 伦理电影免费视频| 国产亚洲av片在线观看秒播厂| 久久99精品国语久久久| 久久久精品94久久精品| 中文字幕免费在线视频6| 精品熟女少妇av免费看| 久久ye,这里只有精品| 涩涩av久久男人的天堂| 亚洲精品日韩av片在线观看| 能在线免费看毛片的网站| 一区二区三区免费毛片| 久久国产精品男人的天堂亚洲 | 一级毛片我不卡| 在线观看一区二区三区激情| a级毛片在线看网站| 日本av手机在线免费观看| 久久精品国产亚洲av天美| 免费观看在线日韩| 国产亚洲最大av| 免费黄网站久久成人精品| 99视频精品全部免费 在线| 久久久久久久精品精品| 中文字幕久久专区| 如日韩欧美国产精品一区二区三区 | 制服丝袜香蕉在线| 亚洲天堂av无毛| 嫩草影院入口| 亚洲三级黄色毛片| 国产精品无大码| 欧美日韩国产mv在线观看视频| 国产又色又爽无遮挡免| 69精品国产乱码久久久| 中文字幕精品免费在线观看视频 | 亚洲精品色激情综合| 少妇熟女欧美另类| 人妻制服诱惑在线中文字幕| 亚洲国产精品专区欧美| 男人和女人高潮做爰伦理| 国产毛片在线视频| 亚洲国产精品专区欧美| 国产成人aa在线观看| 日韩,欧美,国产一区二区三区| 精品午夜福利在线看| 一级av片app| 熟女电影av网| 国产日韩一区二区三区精品不卡 | 婷婷色综合www| 在线播放无遮挡| 99热网站在线观看| 欧美日韩av久久| 免费大片黄手机在线观看| 久久久国产一区二区| 妹子高潮喷水视频| 亚洲av日韩在线播放| 高清不卡的av网站| 成人综合一区亚洲| 26uuu在线亚洲综合色| 国产精品一区www在线观看| 国产伦理片在线播放av一区| 久久久久视频综合| 久久精品国产亚洲网站| 成人亚洲欧美一区二区av| a级毛片免费高清观看在线播放| 免费黄网站久久成人精品| 久久热精品热| 99久久综合免费| 久热这里只有精品99| 搡女人真爽免费视频火全软件| 国产日韩欧美视频二区| 国产av码专区亚洲av| 美女xxoo啪啪120秒动态图| 91成人精品电影| 国产一区二区在线观看日韩| 国产精品熟女久久久久浪| 国产一区二区三区综合在线观看 | 国产淫语在线视频| 国产真实伦视频高清在线观看| 成年人午夜在线观看视频| 伦精品一区二区三区| 日韩中文字幕视频在线看片| 麻豆乱淫一区二区| 成年美女黄网站色视频大全免费 | 中文字幕av电影在线播放| 两个人免费观看高清视频 | 香蕉精品网在线| 2022亚洲国产成人精品| 亚洲四区av| 精品酒店卫生间| 成人特级av手机在线观看| 久久久国产欧美日韩av| 哪个播放器可以免费观看大片| 久久久久网色| 看十八女毛片水多多多| 美女xxoo啪啪120秒动态图| 久久国产精品男人的天堂亚洲 | 久热这里只有精品99| 99久久中文字幕三级久久日本| 如何舔出高潮| 在线观看av片永久免费下载| videossex国产| 精品人妻熟女毛片av久久网站| 中文精品一卡2卡3卡4更新| 亚洲图色成人| 久久99蜜桃精品久久| 人妻少妇偷人精品九色| 不卡视频在线观看欧美| 天堂俺去俺来也www色官网| 高清午夜精品一区二区三区| 久久国产亚洲av麻豆专区| 国产又色又爽无遮挡免| 2018国产大陆天天弄谢| 久久精品国产鲁丝片午夜精品| 国产精品偷伦视频观看了| 欧美日韩av久久| 亚洲av日韩在线播放| 一区二区三区乱码不卡18| 免费大片黄手机在线观看| 狂野欧美激情性xxxx在线观看| 国产免费又黄又爽又色| 极品人妻少妇av视频| 国产精品蜜桃在线观看| 汤姆久久久久久久影院中文字幕| 女性被躁到高潮视频| 亚洲国产日韩一区二区| 国产高清国产精品国产三级| 久久狼人影院| 一本色道久久久久久精品综合| 中文字幕免费在线视频6| 国产黄频视频在线观看| 日韩中文字幕视频在线看片| 免费看日本二区| 日产精品乱码卡一卡2卡三| 国产精品嫩草影院av在线观看| 日韩,欧美,国产一区二区三区| 国产精品一区二区在线不卡| 一级毛片aaaaaa免费看小| 亚洲精品乱久久久久久| 内地一区二区视频在线| 日本色播在线视频| 99热这里只有是精品在线观看| 在线观看av片永久免费下载| 国产老妇伦熟女老妇高清| 高清午夜精品一区二区三区| 亚洲成人手机| 国产成人freesex在线| 天堂8中文在线网| 亚洲精品一二三| 多毛熟女@视频| 丝袜喷水一区| 久久久久久伊人网av| 一级av片app| 亚洲精品久久久久久婷婷小说| 亚洲av在线观看美女高潮| 日韩人妻高清精品专区| 99久久综合免费| 黄色一级大片看看| 99九九线精品视频在线观看视频| av在线app专区| 国产极品粉嫩免费观看在线 | 国产精品一区二区在线不卡| 亚洲国产毛片av蜜桃av| 99久久精品热视频| 中文字幕制服av| 亚洲精品久久久久久婷婷小说| 女性被躁到高潮视频| 欧美三级亚洲精品| 精品卡一卡二卡四卡免费| 久久精品国产亚洲av涩爱| 一级二级三级毛片免费看| 久久毛片免费看一区二区三区| 777米奇影视久久| 国产欧美日韩一区二区三区在线 | 亚洲精品日韩av片在线观看| 性色av一级| 日产精品乱码卡一卡2卡三| 少妇 在线观看| 这个男人来自地球电影免费观看 | 亚洲国产欧美在线一区| av在线老鸭窝| 一级毛片 在线播放| 亚洲精品国产av蜜桃| 日韩中字成人| 久久精品久久久久久久性| 国产精品偷伦视频观看了| 在线观看www视频免费| 自拍偷自拍亚洲精品老妇| 青春草国产在线视频| 男人舔奶头视频| 国产 精品1| 22中文网久久字幕| 亚洲精品国产av成人精品| 美女脱内裤让男人舔精品视频| 国产在视频线精品| 免费观看av网站的网址| 美女cb高潮喷水在线观看| 又黄又爽又刺激的免费视频.| av又黄又爽大尺度在线免费看| 在线观看一区二区三区激情| 亚洲国产欧美在线一区| 男人添女人高潮全过程视频| 2018国产大陆天天弄谢| 久久婷婷青草| av又黄又爽大尺度在线免费看| 成人美女网站在线观看视频| 女人精品久久久久毛片| 久久久国产欧美日韩av| 亚洲精品乱久久久久久| 一级爰片在线观看| 搡女人真爽免费视频火全软件| 国产深夜福利视频在线观看| 亚洲怡红院男人天堂| 亚洲图色成人| av免费在线看不卡| 99热这里只有是精品50| 久久久国产一区二区| 国产精品麻豆人妻色哟哟久久| 欧美精品国产亚洲| 久久久久久久亚洲中文字幕| 中国三级夫妇交换| 欧美精品高潮呻吟av久久| 亚洲内射少妇av| 另类亚洲欧美激情| 一级毛片久久久久久久久女| 久久鲁丝午夜福利片| 日日摸夜夜添夜夜爱| 欧美最新免费一区二区三区| 成人无遮挡网站| 欧美三级亚洲精品| 国产亚洲精品久久久com| 国产在线免费精品| 久久久久久久精品精品| 亚洲中文av在线| 性色avwww在线观看| 精品熟女少妇av免费看| 亚洲av免费高清在线观看| 国产免费一级a男人的天堂| 欧美精品亚洲一区二区| 91午夜精品亚洲一区二区三区| 最近的中文字幕免费完整| 亚洲精品,欧美精品| 男人狂女人下面高潮的视频| 亚洲欧美成人综合另类久久久| 精品亚洲成国产av| 欧美精品国产亚洲| kizo精华| 视频中文字幕在线观看| 九草在线视频观看| 日韩免费高清中文字幕av| 永久网站在线| 最近的中文字幕免费完整| 91aial.com中文字幕在线观看| 久热久热在线精品观看| 国产精品秋霞免费鲁丝片| 自线自在国产av| 精品午夜福利在线看| av专区在线播放| 偷拍熟女少妇极品色| 女性生殖器流出的白浆| 亚洲va在线va天堂va国产| 国内少妇人妻偷人精品xxx网站| 麻豆精品久久久久久蜜桃| 26uuu在线亚洲综合色| 精品亚洲成国产av| 一级a做视频免费观看|