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

    基于精細微觀交通流模型的信號交叉口人-車相互干擾研究*

    2019-12-24 08:21:22陶亦舟韋艷芳高慶飛董力耘
    物理學報 2019年24期
    關鍵詞:人行橫道交通燈元胞

    陶亦舟 韋艷芳 高慶飛 董力耘?

    1) (上海應用技術大學理學院, 上海 201418)

    2) (玉林師范學院物理科學與工程技術學院, 玉林 537000)

    3) (上海大學力學與工程科學學院, 上海市應用數(shù)學和力學研究所, 上海 200072)

    4) (上海市力學在能源工程中的應用重點實驗室, 上海 200072)

    采用微觀離散模型研究信號交叉口的人-車相互干擾機理, 其中車輛運動描述基于細化NaSch模型并考慮了司機對交通燈信號切換的預期效應; 行人運動描述則基于多步格子氣模型并考慮了過街綠燈周期內(nèi)行人過街速度逐漸增加的特征.機動車與行人之間相互干擾表現(xiàn)為交通燈信號切換時由于車輛(行人)滯留沖突區(qū)而造成行人(車輛)運動的延誤.通過數(shù)值模擬得到了車輛運動的基本圖、行人等待時間以及相應的相圖, 并給出行人(車輛)滯留沖突區(qū)內(nèi)造成車輛(行人)延遲的定量特征.模擬結果發(fā)現(xiàn)存在一個臨界綠信比.當綠信比小于臨界值時, 隨著密度的增大, 車流出現(xiàn)自由流相、飽和流相和擁堵流相; 而當綠信比大于臨界值, 則出現(xiàn)自由流相、共存相、飽和流相和擁堵流相.由人-車相互干擾所導致的延誤與車流和行人流的運行狀態(tài)密切相關.當行人到達率和綠信比都比較大時, 等待區(qū)內(nèi)的行人無法在行人綠燈周期內(nèi)一次清空.文中詳細地討論了人-車相互干擾的定性和定量特征.發(fā)現(xiàn)通過合理設置綠信比, 可以在保證車流運行效率的同時, 減少行人過街的等待時間.

    1 引 言

    交叉路口是城市道路網(wǎng)絡中的重要節(jié)點, 也是瓶頸所在.在流量較大的交叉路口通常都安裝了交通燈, 屬于信號交叉口.交叉口機動車與行人之間的沖突主要發(fā)生在人行橫道處.由于交通燈的調(diào)控作用, 行人過街的安全性顯著提高, 但仍然存在著機動車與行人的相互干擾.即使不考慮行人綠燈期間過街行人與右轉車輛之間的沖突, 當直行車輛在綠閃期間越過停車線, 若交通燈由綠轉紅時車輛仍繼續(xù)通過路口則會導致行人運動延誤.另一方面,由于大城市中不少主干道的車道數(shù)偏多, 路面過寬致使人行橫道過長, 當行人綠燈相位配時不足時,往往當車輛綠燈亮起時, 行人還沒有走出人行橫道, 從而導致車輛運動的延誤.因此很有必要研究信號交叉口處車輛和行人之間的相互干擾.

    交叉口處機動車和行人沖突不但會降低交叉口通行效率, 還會導致交通安全事故, 國內(nèi)外許多學者對此開展了多方面的研究:對人行橫道上人群過街行為進行觀測, 所得到的行人運動定量特征可用于行人流建模和模型參數(shù)的標定[1-4]; 從安全的角度對行人過街的安全性進行評估[5,6]; 對人行橫道上的行人過街行為進行動力學建模和模擬[7-11];對信號交叉口的機動車流進行動力學建模和模擬[12-15].這些模型多從微觀角度描述機動車或行人個體的運動, 如行人流的社會力模型[16]、格子氣模型[17], 機動車的優(yōu)化速度模型[18]以及行人和機動車的元胞自動機模型[19-21], 它們可以細致地反映車輛之間、行人之間以及人-車之間的復雜相互作用.以往的研究中多以機動車流或行人流中的一種為主要研究對象.如Zeng等[10]基于行人心理,建立了考慮群體規(guī)避和個體規(guī)避、轉彎車輛避碰等典型交叉口行人行為的改進社會力模型; Belbasi和Foulaadvand[12]采用NaSch模型研究單個信號交叉口的機動車流特征, 并考慮了固定時長的交通燈和自適應交通燈.

    最近對機動車與行人之間干擾的建模和模擬越來越受到重視[22-30], 其中元胞自動機模型是近年來被廣泛用于交通流模擬的一類模型.由于在信號交叉口行人和機動車之間存在著相互干擾, 因此對機動車流和行人流同步模擬才能更好地反映人-車相互作用的機理.Xie等[24]采用元胞自動機模型研究了信號人行橫道的混合交通特性, 考慮了人-車之間的相互作用機制, 并引入謹慎、正常和冒險三種不同類型的行人, 發(fā)現(xiàn)三種類型的行人在頻率曲線上均存在兩個局部最大值, 這一現(xiàn)象與真實的交通狀況是一致的.Li等[25]提出了一種元胞自動機模型, 用于模擬信號交叉口行人和車輛的過馬路行為, 研究了不同行人信號周期和人行橫道寬度對人行橫道通行能力、行人與車輛之間的沖突數(shù)量以及沖突導致的行人延誤.孫澤和賈斌[26]通過引入行人和機動車的沖突干擾規(guī)則, 構建了能夠描述人行橫道處機動車和行人相互干擾行為的元胞自動機模型, 并通過數(shù)值模擬研究了行人到達率、機動車到達率、行人冒險穿越機動車道的閾值以及行人對機動車敏感系數(shù)等因素對機動車和行人流量的影響.余艷等[27]基于 VDR (velocity-dependentrandomization)模型[21], 建立行人與機動車相互干擾的元胞自動機模型, 在開放邊界條件下研究有、無交通燈控制時車輛產(chǎn)生概率、消失概率和交通燈綠信比對車輛交通流和行人流的影響.此外, 學者們對無信號交叉口的人-車相互干擾也開展了大量的研究.值得一提的是, 近年來自動駕駛車輛的出現(xiàn)也進一步推動了人-車相互干擾的研究[31].

    總體而言, 目前關于信號交叉口的人-車相互作用已經(jīng)有不少研究成果, 但對于間斷交通流的模擬, 雖然在模型中考慮了交通燈的影響, 然而未能充分考慮司機根據(jù)交通燈切換進行的速度調(diào)整, 往往出現(xiàn)加速度過大的情況, 因此如何準確刻畫司機對交通燈的反應仍然是一個值得研究的問題.此外, 行人過街的一個引人注意的特征是隨著行人綠燈時間的減少, 尚未過街的行人會加快速度, 因而在行人綠燈周期內(nèi)行人平均過街速度并非一個常數(shù), 而是越來越快[3].這一特征尚未在已有的工作中得到充分的展現(xiàn).基于上述分析, 本文根據(jù)信號交叉口的機動車和行人流運動特征, 分別建立了相應的微觀離散模型, 對機動車和行人相互作用進行了細致的研究, 并深入分析人-車相互干擾的機理.

    2 行人-車輛運動的微觀離散模型

    在已有的機動車-行人微觀離散模型中, 從模擬效率的角度考慮, 常采用較為簡單的模型, 如Nagel和 Schreckenberg[19]提出的 NaSch模型以及Blue和Adler[20]提出的行人流元胞自動機模型.NaSch模型[19]原用于描述高速公路交通流(屬于連續(xù)交通流), 雖然目前已經(jīng)提出了很多機動車元胞自動機模型, 但由于NaSch模型規(guī)則簡單, 計算效率高, 因此得到廣泛應用.盡管NaSch模型不能再現(xiàn)同步流, 但是由于城市地面道路通常較短, 受到交通燈的制約, 難以出現(xiàn)同步流等自組織現(xiàn)象.基于上述考慮, 本文所提出機動車運動模型仍是基于NaSch模型, 并充分考慮交通燈對司機駕駛行為的影響.本文假設人-車不能共存于沖突區(qū)(見圖1), 即在人行橫道內(nèi)車輛不會作為障礙物而出現(xiàn), 這種行人運動場景相對簡單, 因此行人過街可用較為簡單的格子氣模型來描述, 相比于元胞自動機模型, 其優(yōu)點是采用隨機串行更新規(guī)則, 避免了并行更新時發(fā)生的行人沖突[32,33].

    2.1 研究場景和基本假設

    圖1 信號交叉口人-車相互作用示意圖Fig.1.Schematic diagram for vehicle-pedestrian interaction at a signalized crossing.

    將信號交叉口處行人過街和車輛通行的場景簡化如下:長為L的車道中間設有一個與其垂直的人行橫道.車道與人行橫道的重疊部分稱為沖突區(qū).在沖突區(qū)的下游邊緣設置一個交通燈, 如圖1所示.圖中位于停車線后的紫色矩形區(qū)域代表一輛車, 箭頭所示為其運動方向, 即車輛由左向右行駛;三角形代表向上運動的行人.停車線附近的車道可分為三部分:沖突區(qū)、沖突區(qū)上游和下游.在沖突區(qū)上游(即停車線左側)設置了一個減速區(qū), 位于減速線和停車線之間, 其長度為LD.車輛由沖突區(qū)的下游離開.人行橫道包括兩個部分:沖突區(qū)和等待區(qū).沖突區(qū)長為WC, 寬為W(即車道的寬度).沖突區(qū)的上游為行人等待區(qū), 其長度為LW.

    本文模型的基本假設如下:

    1) 行人從固定位置(即人行橫道)過街, 行人的運動限制在人行橫道(包括沖突區(qū))內(nèi), 且行人是從下到上單向運動, 即不考慮雙向的人群運動;

    2) 僅考慮單車道的情況, 車輛從左到右運動;

    3) 行人和車輛在沖突區(qū)內(nèi)是分離的, 即行人和車輛不允許同時出現(xiàn)在沖突區(qū)內(nèi);

    4) 交通燈周期分為兩個相位, 車輛綠燈相位時, 車輛可以通行, 行人需在路邊等待, 其周期為TG; 車輛紅燈相位即行人綠燈相位, 其周期為TR,此時行人可以過街, 車輛則需在停車線后等待.

    2.2 行人運動模型

    本文僅考慮同向行人過街的情況.由于人行橫道內(nèi)沒有障礙物, 且行人運動區(qū)域的幾何形狀簡單, 因此行人的期望運動方向明確, 這里我們采用最基本的格子氣模型來描述行人的運動.格子氣模型可以視為一種特殊的元胞自動機模型, 通常采取隨機串行更新規(guī)則.將人行橫道離散成由正方形元胞組成的二維網(wǎng)格.一個元胞的大小就是典型行人所占的面積, 即 0.4 m × 0.4 m.

    李文勇等[2]的觀測表明行人步行速度在過街綠燈周期的不同階段會有所變化:當行人綠燈信號啟亮, 平均步行速度為 1.00—1.35 m/s; 當行人綠燈時間過半, 平均步行速度為 1.28—1.56 m/s; 當行人綠燈閃爍, 平均步行速度為1.56—2.12 m/s.可見行人在過街綠燈周期內(nèi)的平均步行速度是隨著時間的增加而逐步增大, 并非恒定不變.此結果與日常經(jīng)驗相符, 應當在行人流模型中有所體現(xiàn).由上述實測數(shù)據(jù)可知, 行人最大速度Umax可 達2 m/s左右, 即每秒運動 5 個元胞.本文采用多步模型來實現(xiàn), 即將 1 s分為 5個子時間步, 每個子時間步行人最多移動1個元胞.

    行人運動規(guī)則如下.

    1)根據(jù)交通燈周期確定行人的運動概率

    a) 在車輛綠燈周期, 設行人直行概率p=p0,起動概率q=q0, 且 0<p0<1 , 0<q0<1 .如果行人尚在沖突區(qū)內(nèi), 則直接設p=1 ,q=1 , 即行人將會用確定性的方式以最快速度離開沖突區(qū).如果行人尚未進入沖突區(qū), 則到達路邊(即等待區(qū)上邊界)后, 將會停下等待.其中, 行人直行概率p反映了行人選擇出行目的方向(即下游出口方向)概率的權重, 起動概率q則反映行人運動過程中可能出現(xiàn)的延遲, 與VDR模型[21]中車輛的慢啟動類似.

    b)在行人綠燈周期, 行人可以穿越人行橫道.為反映行人在過街周期內(nèi)的運動特征, 假設行人直行概率p和起動概率q隨時間t線性增加:

    其中TR為行人綠燈周期長度; 時間t∈[0,TR] , 當行人綠燈啟亮瞬間, 有t=0 .

    表1 各種情況下行人向鄰近元胞的行走概率Table 1.Probabilities for a pedestrian to move towards his/her neighboring cells.

    2)行人根據(jù)相鄰元胞的狀態(tài)確定目標元胞

    對于行人而言, 當r<q時才可以運動.其中r是在 [0 ,1] 之間均勻分布的隨機數(shù).如不特別說明,后文中出現(xiàn)的r總是指這種類型的隨機數(shù).

    如果行人i可以運動, 則有四種可能選擇, 即向前、向左或向右運動, 還可以停在原地不動.行人i左側、前方和右側元胞的狀態(tài)Si用一個數(shù)組表示 (sL,sF,sR) , 其中每個分量可取 0 或 1, 表示對應的元胞為空或被占據(jù), 如 (0 ,1,0) 表示行人左右兩側元胞無人, 而其前方的元胞被占據(jù).行人根據(jù)概率 (pL,pF,pR) 隨機選擇目標元胞, 然后進入該元胞或停留在原位, 如表1 所列.其中pL,pF和pR分別表示行人向左、向前和向右運動的概率.

    采用隨機串行更新方式逐一更新行人的位置,即在每一子時間步, 都會對所有行人進行隨機排序.然后依次逐一確定行人的運動并更新.注意在1 s內(nèi)行人運動要更新5次.人行橫道采用開放邊界條件.假設行人的到達率λ滿足泊松分布.行人等待區(qū)最多可以容納LW×WC個人, 每隔 1 s會在等待區(qū)的下邊界隨機產(chǎn)生一些新的行人, 并向上游運動.本文取λmax=2 , 相當于 1 h 內(nèi)有 7200 人通過路口(對于一般路口而言, 這已是足夠大的行人流量).行人在邊界以一定概率γ消失.為簡單起見,本文中取γ=1 , 即行人達到?jīng)_突區(qū)的上邊界處可自由離開.

    2.3 車輛運動模型

    為了準確描述車輛的起動和剎車過程, 本文在網(wǎng)格細化的NaSch模型上考慮交通燈的影響.為與行人流模型匹配, 車道也離散成長度為0.4 m的元胞.根據(jù)假設3), 暫不考慮車輛的橫向運動, 因此沖突區(qū)以外的車道不必離散成二維的網(wǎng)格.當交通燈切換時, 需要考慮司機的預期反應, 為此在停車線上游設置長為LD的減速區(qū).當交通燈由綠轉紅時, 減速區(qū)內(nèi)的車輛會采取相應的減速措施, 避免出現(xiàn)交通事故.

    對于車輛i而言, 在t時刻它的位置為xi(t) , 速度為vi(t) , 車長均為lc.其運動規(guī)則如下.

    1)加速

    其中a+為車輛加速時的加速度值.

    2)出于安全考慮對于車速的限制

    其中di(t) 為車輛i與前車的間距,

    這意味著在本模型中相鄰兩輛車的最小間距為1.

    3)減速

    a) 對于車輛綠燈周期的車輛或車輛紅燈周期內(nèi)但不在減速區(qū)內(nèi)的車輛i, 考慮司機駕駛行為中的不確定因素而導致的隨機減速:

    其中pd為隨機減速概率,a-為隨機減速時的加速度值.

    b)對于車輛紅燈周期內(nèi)處于減速區(qū)內(nèi)的車輛i:

    其中?i=xS-xi(t) 是車輛i到停車線的距離, 且?i≤LD,xS是停車線的位置,是正常減速時允許的最大加速度值.

    c)當交通燈由綠轉紅時, 對于停車線后的第一輛車j, 如果它即將到達停車線, 則根據(jù)當前速度預測該車是否可以順利在停車線前減速停下.如果

    則車輛j將會在停車線后停下, 否則它將繼續(xù)駛過路口.

    4)運動

    車輛運動采用并行更新方式, 而車道采用周期性邊界條件.

    2.4 人-車相互作用

    由于交通燈的調(diào)控, 直行車輛和過街行人的運動基本分離, 而當交通燈切換時, 車輛與行人之間將發(fā)生相互干擾.

    1)當交通燈由綠轉紅時, 即將駛過停車線的車輛如果來不及減速, 會在紅燈期開始后繼續(xù)駛過沖突區(qū), 此時行人會預判車輛的行為(在模擬中雖然已是紅燈, 沖突區(qū)仍將設置為行人不可通行), 因此并不會立即進入沖突區(qū), 只有當車輛駛出沖突區(qū)后, 行人才可以過街, 因此會造成行人運動的延誤.此外, 如果車輛滯留在沖突區(qū)也會導致行人運動的延誤.每小時內(nèi)由車輛導致的行人平均時間延誤用表示.

    2)當交通燈由紅轉綠時, 尚未進入沖突區(qū)的行人將停在路邊等待, 而在沖突區(qū)內(nèi)尚未過街的行人將繼續(xù)前行, 直到?jīng)_突區(qū)內(nèi)無人時, 車輛才可以通行(在模擬中雖然已是綠燈, 沖突區(qū)仍將設置為車輛不可通行).這種情況下, 由于行人滯留在沖突區(qū)而造成車輛運動的延誤.每小時內(nèi)由行人導致的車輛平均時間延誤用表示.

    3 數(shù)值模擬與分析

    3.1 參數(shù)設置

    模擬時車道長度L=3000 , 車道寬度W=10 ,減速區(qū)長度LD=120 , 人行橫道寬度WC=W, 等待區(qū)長度LW=20 .假設所有車輛都是同樣的轎車,車長lc=10 , 最大速度Vmax=40 , 加速時的加速度值a+=4 , 減速時的最大加速度值這樣車輛可以在10 s加速到最大速度, 并可以在4 s以內(nèi)減速到零, 可以較好地反映車輛的實際運動情況.車輛的隨機減速概率pd=0.3 .

    假設每個行人都只占據(jù)人行橫道上的一個元胞.行人到達符合泊松分布, 到達率為λ, 且0< λ≤2.設p0=0.8 ,q0=0.6 , 據(jù) 此可以估 算,在行人綠燈周期開始時行人平均速度約為1.2 m/s,即將結束時約為2.0 m/s.

    設交通信號燈周期時長T= 100 s, 即車輛綠燈周期TG與行人綠燈周期TR之和:T=TG+TR.綠信比μ定義為車輛綠燈周期TG與交通燈周期時長T的比值, 即

    可知TR=(1-μ)T.本文中設μ≥0.5 .

    車輛密度k定義為

    其中Nc為車輛總數(shù).

    所有車輛的平均速度V定義為

    其中T0是系統(tǒng)達到穩(wěn)態(tài)的時間,T1是計算速度平均值所用的總時間步長.

    平均車流量Jc定義為

    3.2 典型參數(shù)下模型的基本性質

    首先研究不同綠信比μ時車輛運動的基本圖.模擬結果顯示, 存在一個臨界綠信比μc≈ 0.8 .當μ≤μc時, 可以觀察到兩個臨界密度k1和k2, 見圖2(a).當k<k1時, 車流處于自由流相, 其特征為車流量Jc隨密度k的增大而線性增加, 平均車速為常數(shù),見圖2(b); 當k1≤k≤k2時, 車流處于飽和流相, 在流密曲線中出現(xiàn)一個流量平臺, 此時車流量Jc約為一常數(shù)J0, 表明綠燈期間通過的車輛數(shù)大致相等.這是交通燈控制下間斷交通流的普遍特征; 當k>k2時, 是擁堵流相, 此時流量隨密度增大而線性減少, 這是由于上游路段的剩余空間限制了綠燈期間可以通過的數(shù)量.當k>k1時, 平均速度隨密度的增加而單調(diào)減少.由圖2可見, 不同行人到達率λ對車輛運動基本圖的影響并不明顯, 表明車輛的運動主要取決于綠信比, 與實際情況相符.隨著綠信比μ的增大, 車輛通行時間相應增加, 飽和流量值J0隨之增大, 同時k1逐漸增大, 而k2逐漸減小, 即流量平臺出現(xiàn)的密度范圍逐漸減小, 見圖2(c).這表明, 當μ較大時, 雖然車流量更大, 但也更容易發(fā)生交通擁堵.值得注意的是, 在圖2(d) 中, 當μ=0.9> μc時, 此時未飽和流相可以分為兩段, 當k<k0時, 車流量隨密度線性增加, 且V≈Vmax=40 , 為自由流相;當k>k0時, 平均速度隨密度的增加而單調(diào)減小,但當k0≤k<k1時, 車流量仍隨密度線性增加, 見圖2(c).

    圖2 不同綠信比下車輛運動基本圖 (a)流量-密度曲線 (μ = 0.5); (b)速度-密度曲線 (μ = 0.5); (c)流量-密度曲線 (μ = 0.9);(d)速度-密度曲線 (μ = 0.9)Fig.2.Fundamental diagrams for vehicles under different splits:(a) Flux-density relation (μ = 0.5); (b) speed-density relation (μ =0.5); (c) flux-density relation (μ = 0.9); (d) speed-density relation (μ = 0.9).

    圖3給出了典型密度下的車輛運行時空斑圖,是模擬 1 h 后取兩個交通燈周期、處于 [xS— 600,xS+ 600]范圍內(nèi)的車輛運動軌跡, 其中每個點代表一輛車, 綠色和紅色豎線分別表示處于車輛綠燈和紅燈周期內(nèi).圖3(a)為自由流相, 車輛在紅燈周期結束前開始在停車線后排隊, 在綠燈周期結束前全部通過路口.這表明車輛不會直接影響行人過街.從圖中可以觀察到明顯的剎車波和起動波, 小尺度的交通擁堵向上游傳播并消散.圖3(b)為飽和流相, 紅燈期間出現(xiàn)車輛排隊, 由于此時排隊的車輛足夠多, 在綠燈期間通過路口的車輛數(shù)主要取決于綠燈周期的長度.在飽和流相, 綠燈期間通過路口的車輛數(shù)基本不變.圖3(c)為擁堵流相, 可以發(fā)現(xiàn)飽和流與擁堵流之間的差別.擁堵流相時下游車輛的擁堵不但會向上游傳播, 還會影響到上游路口車輛的通行.圖3(d)給出了k= 0.15 (k0<k<k1)時的車輛運動斑圖.可以發(fā)現(xiàn)既有車輛自由通過路口, 也有紅燈期間形成小尺度擁堵的消散過程, 因此本文將圖2(c)和圖2(d)中k0<k<k1范圍內(nèi)的交通流相稱之為共存相, 其平均車速隨著車輛密度增加而減小, 見圖2(d).當k<k0時, 車輛均可自由通過路口, 因此V≈Vmax.

    3.3 行人等待時間的影響

    行人等待時間TW定義為從行人進入等待區(qū)到行人離開等待區(qū)(即進入沖突區(qū))的時間之差的平均值.圖4反映了四種典型行人到達率下車輛運動對于行人等待時間的影響, 可以觀察到車輛運動基本圖(參考圖2)與行人等待時間兩者之間的相關性.當μ≤μc時, 發(fā)現(xiàn)車流處于自由流相(k<k1)時,TW可以視為不隨車輛密度而變化的常數(shù), 從圖3(a)可知, 這是由于車輛綠燈周期結束之前, 所有車輛已經(jīng)通過路口, 因而不會影響行人過街.當車流是飽和流狀態(tài) (k1≤k≤k2)時, 車隊在綠燈期間的消散過程類似, 對行人過街的影響大致相同, 可以預期會出現(xiàn)一個等待時間平臺,TW也近似為一個常數(shù).由于行人綠燈周期足夠長, 在車輛綠燈期間集聚的行人基本可以一次通過, 但是若交通燈由綠轉紅時仍有車輛通過路口, 行人只能等其離開沖突區(qū)再過街, 因而導致TW略有增大.當k>k2時, 車輛處于擁堵流狀態(tài)時, 車輛可能滯留在沖突區(qū), 因而導致TW明顯增加, 特別是行人到達率λ很大時, 見圖4(b).當μ=0.9> μc時, 由于行人綠燈周期過短, 即使車輛密度很小, 當行人到達率λ=2時因為等待過街行人很多而無法一次清空, 因此行人等待時間很長, 見圖4(c).隨著μ的增加,TW也相應增大.當行人到達率λ還不夠大時(如λ=1),此時可將人行橫道視為十車道情形, 行人以最大速度Umax=5 運動, 因此在較短時間內(nèi), 可以通過相當多的行人.在車輛綠燈期間積累的人群在行人綠燈周期都可以清空(不考慮在車輛綠燈周期即將結束時剛進入等待區(qū)的行人).當λ=2 時, 行人綠燈啟亮時, 等待區(qū)已經(jīng)充滿了行人.當紅轉綠時, 仍有不少行人滯留在等待區(qū).這也使得在車輛綠燈周期內(nèi), 等待區(qū)很容易被人群充滿.

    圖3 車輛運動時空演化斑圖 (λ =1.0) (a) k = 0.08, μ =0.5; (b) k = 0.40, μ = 0.5; (c) k = 0.70, μ = 0.5; (d) k =0.15, μ = 0.9Fig.3.Temporal-spatial patterns of vehicles in motion(λ =1.0):(a) k = 0.08, μ = 0.5; (b) k = 0.40, μ = 0.5;(c) k = 0.70, μ = 0.5; (d) k = 0.15, μ = 0.9.

    圖4 不同綠信比下行人等待時間隨車輛密度的變化(a) μ = 0.5; (b) μ = 0.8; (c) μ = 0.9Fig.4.Waiting time of pedestrian varies with vehicle density under different splits:(a) μ = 0.5; (b) μ = 0.8; (c) μ =0.9.

    圖5 不同綠信比下行人等待時間隨行人到達率的變化(a) μ = 0.5; (b) μ = 0.8; (c) μ = 0.9Fig.5.Waiting time of pedestrians varies with the arrival rate of pedestrian under different splits:(a) μ = 0.5; (b) μ =0.8; (c) μ = 0.9.

    圖5給出了四種典型車輛密度下行人等待時間TW隨行人到達率λ的變化趨勢.由圖5 可見, 對于給定的車輛密度k, 隨著行人到達率λ的增加,行人等待時間TW相應增大.對于同一綠信比和行人到達率, 車輛密度越大, 行人等待時間越長.當μ=0.5,0.6 時, 車輛處于自由流與飽和流(k=0.1,0.3,0.5<k2, 參考圖2(a)), 由于車輛運動對行人過街的干擾很少, 且行人綠燈周期夠長, 因此TW隨λ的增加而緩慢增大.對于擁擠流狀態(tài)k=0.7>k2, 由于車輛會暫時滯留在沖突區(qū), 對行人過街影響很大, 但是因為紅燈周期足夠長, 行人仍可順利過街.當μ=0.7,0.8 時, 對于k= 0.1, 0.3 <k2, 車流尚未超出飽和流相范圍,TW隨λ緩慢增長.但由于車輛綠燈周期長, 行人的等待時間明顯增加.k=0.5,0.7>k2時車流處于擁擠流相.當λ < λc時,TW隨λ緩慢增長, 由于此時人數(shù)較少,雖然行人綠燈周期較短, 仍可基本一次清空.λc依賴于μ和k.當λ(λ > λc≈ 1.0) 足夠大時,TW會較快增長, 如圖5(b) 所示.當λ足夠大時,TW>T(明顯大于交通燈周期), 這表明有些行人在一個交通燈周期內(nèi)沒能過街.值得注意的是,μ=0.9 時,由于行人綠燈周期很短, 即使車輛密度很小(如k=0.1) 時, 當λ足夠大時,TW會快速增長, 見圖5(c),這是由于車輛綠燈周期過長, 在等待區(qū)積累的大量行人無法在很短的行人綠燈周期內(nèi)通過人行橫道,而滯留的行人又容易導致等待區(qū)滿員, 但與高車輛密度情況相比, 仍有明顯的差別, 這表明滯留在沖突區(qū)的車輛對行人過街影響很大.此外, 當λ較小時, 不同車輛密度下TW隨λ變化曲線基本重合.這表明當行人來流較小時, 絕大部分行人還可以順利過街, 即使行人綠燈周期很短.

    3.4 人-車相互作用

    下面研究人-車相互干擾所導致的延誤, 即由于車輛紅燈(綠燈)期間車輛(行人)占據(jù)沖突區(qū)所導致的行人(車輛)的延誤.首先研究車輛對行人的影響, 包括車輛紅燈周期內(nèi)滯留在沖突區(qū)的車輛或者交通燈剛由綠轉紅時才進入沖突區(qū)的車輛.出于安全考慮, 即使行人綠燈已經(jīng)啟亮, 此時行人也不會進入沖突區(qū), 因此這兩種情況下車輛都會占用行人過街時間.

    圖6 不同綠信比時車輛所導致的行人延誤 (a) μ = 0.5;(b) μ = 0.8; (c) μ = 0.9Fig.6.Pedestrian delay time caused by vehicles under different splits:(a) μ = 0.5; (b) μ = 0.8; (c) μ = 0.9.

    下面研究由于行人未能過街而滯留對車輛通行造成的影響.從圖7可以看出, 當行人到達率較小時(或者行人綠燈周期長), 行人對于車輛的影響主要是陸續(xù)到達的行人直接通過路口(車輛綠燈周期在路口累積的等待過街行人已經(jīng)在行人綠燈亮起后率先過街了), 當車輛紅燈已轉綠, 此時尚未完全通過路口的行人會造成車輛運動的延誤, 因而隨行人到達率的增加而增大.圖7(b)中較高車輛密度時, 由于上游車輛擁堵延伸到路口, 阻礙了行人過街, 加上紅燈周期時間不夠長, 在路口聚集的大量行人在較少的通行時間內(nèi)集中過街而導致車輛的延誤明顯增大.

    由圖7(a)可知, 當μ較小時, 行人有足夠的時間過馬路, 影響車輛通行的主要因素是行人綠燈周期內(nèi)滯留在沖突區(qū)內(nèi)的行人.隨著行人到達率λ增加,也相應增大, 但增速逐漸減緩.當車流處于擁堵流狀態(tài)(k=0.7)時, 行人對車輛的影響也略有增大.由圖7(b) 可知, 當?shù)竭_率較大時, 雖然行人較多, 過街時間還足夠, 當k<k2時, 行人對車輛的影響與μ=0.5 時基本一致.由于此時車輛不會占據(jù)沖突區(qū), 行人可以順利過街(雖然行人綠燈周期已經(jīng)明顯變短).當k>k2時(車輛處于擁堵狀態(tài)), 由于車輛擠占行人過街時間, 會造成密集人群集中過街, 因此交通燈紅轉綠時, 行人滯留沖突區(qū)的概率增大, 對車輛產(chǎn)生更大的影響.當行人達到率λ較大時, 圖中出現(xiàn)一個平臺.由于等待區(qū)容量有限, 此λ區(qū)間內(nèi)行人都將充滿等待區(qū), 人群的消散過程類似, 對車輛的影響基本相同, 表現(xiàn)為圖中的平臺.由圖7(c)可知, 當交通燈由紅轉綠時, 仍有大批行人集中進入沖突區(qū).當達到率相當大時,上游來流可以認為近似不變.即使車輛密度很小,由于行人來流足夠大時, 聚集的人群也無法在很短的行人綠燈周期內(nèi)清空.在等待區(qū)存在需要二次清空的人群, 當行人入流較大(λ>1), 等待區(qū)幾乎總是處于滿員狀態(tài), 因此出現(xiàn)了一個的平臺.如果繼續(xù)擴大等待區(qū)的范圍, 由行人滯留沖突區(qū)所造成的車輛延誤還會有所增大.

    圖7 不同綠信比時行人造成車輛的延誤 (a) μ = 0.5;(b) μ = 0.8; (c) μ = 0.9Fig.7.Vehicle delay time caused by pedestrians under different splits:(a) μ = 0.5; (b) μ = 0.8; (c) μ = 0.9.

    3.5 相 圖

    下面研究人-車相互作用系統(tǒng)的全局性質.圖8和圖9分別給出了車流量Jc和行人等待時間TW隨車輛密度k和行人到達率λ的三維分布圖.從圖8可見, 在不同λ下車流量隨k的變化曲線具有類似的特征.當μ< μc時, 車流都可以分為三種狀態(tài),即自由流相、飽和流相和擁堵流相, 見圖8(a).當μ>μc時, 車流都可以分為四種狀態(tài), 即自由流相、共存相、飽和流相和擁堵流相, 見圖8(b).值得說明的是, 當μ> μc時, 從流量上看, 自由流相和共存相的區(qū)分不太明顯, 可以參考圖2(d).隨著綠信比μ的增加, 臨界密度k1增大,k2減小, 飽和流量值相應增大.由圖9可見, 行人等待時間呈現(xiàn)出復雜的空間分布情況.在μ< μc的情況下, 當k<k2時, 行人的等待時間較短, 通常等待時間不會超過TG/2 , 行人綠燈開啟后已經(jīng)在等待區(qū)的行人可以一次過街; 而當k>k2且行人到達率較大時, 等待區(qū)內(nèi)行人很多, 且因過街時間被車輛占用較多, 因此行人綠燈開啟后已經(jīng)在等待區(qū)的人群無法清空, 當μ> μc時尤為明顯, 由于行人綠燈周期很短, 即使車輛密度不大, 當行人來流很大時, 也不可能一次清空等待過街的人群, 見圖9(b).

    圖8 車流量隨車輛密度和行人到達率的變化 (a) μ =0.5; (b) μ = 0.9Fig.8.Flux as a function of vehicle density and pedestrian arrival rate:(a) μ = 0.5; (b) μ = 0.9.

    綜合考慮圖9和圖10, 可以繪制出在 (λ,k) 空間的相圖.當μ< μc時相圖可以分為 4 個區(qū):I區(qū)內(nèi)車輛運動屬于自由流, 行人可順利過街(一次清空); II區(qū)內(nèi)車輛運動屬于飽和流, 行人可順利過街(一次清空); III區(qū)內(nèi)車輛運動屬于擁堵流, 行人仍可順利過街(一次清空); IV區(qū)內(nèi)車輛運動屬于擁堵流, 部分行人需要二次清空.隨著綠信比的增大, II區(qū)減小, 而 I, III和 IV 區(qū)則相應增大.當μ>μc時相圖可以分為8個區(qū):I區(qū)—IV區(qū)內(nèi)車輛運動分別屬于自由流相、共存相、飽和流相和擁堵流相, 行人可一次清空; V區(qū)—VIII區(qū)車輛運動分別屬于自由流相、共存相、飽和流相和擁堵流相,而行人需二次清空.隨著綠信比的增大, 相圖中行人需二次清空所占的區(qū)域也會相應增大.值得說明的是, 上述相圖是在信號周期T= 100 s時給出的.根據(jù)本文所采用的道路長度, 當綠信比較大時,自由流時車輛可以在其綠燈周期內(nèi)進入該路段并通過, 而當綠信比較小時, 則在綠燈周期開始進入該路段的車輛無法在綠燈結束時通過該路段.如果采取滿足上述條件的信號周期, 可以預期將會得到與圖10類似的相圖.如果信號周期過短或者過長時, 則相圖會有較大的變化.本文所采用的信號周期長度屬于真實交通最常見的一種情形.

    圖9 行人等待時間隨車輛密度和行人到達率的變化(a) μ = 0.5; (b) μ = 0.9Fig.9.Waiting time of pedestrian as a function of vehicle density and pedestrian arrival rate:(a) μ = 0.5; (b) μ =0.9.

    圖10 (λ, k)空間中的相圖 (a) μ = 0.5; (b) μ = 0.9Fig.10.Phase diagram in (λ, k) space:(a) μ = 0.5, (b) μ =0.9.

    4 結 論

    基于信號交叉口司機對交通燈切換的預期效應以及行人過街時平均速度隨著行人綠燈時間減少而增大的特征, 建立了描述機動車和行人運動的精細模型.當車輛在遇到紅燈時的剎車過程更加平滑, 其加速度控制在合理的范圍內(nèi), 而對行人過街的刻畫也更接近于實際.人-車相互干擾主要體現(xiàn)在當交通燈切換時, 由于車輛或行人滯留在沖突區(qū)內(nèi)對行人過街或車輛通行所造成的延誤, 我們定義了相應的延誤時間.通過數(shù)值模擬, 給出了人-車相互干擾的定量特征, 并對其機理進行了詳細的討論, 可為提高交叉口的通行效率提供參考.作為人-車相互作用的一種簡化描述, 本文尚未考慮機動車和行人混雜出現(xiàn)在沖突區(qū)內(nèi)的情況.在現(xiàn)實交通中, 即使信號交叉口發(fā)生車輛擁堵甚至溢流的情況, 當行人綠燈啟亮時, 行人仍可穿過車輛的間隙而過街.這種類型的人-車相互干擾將在下一步的工作中繼續(xù)研究.

    猜你喜歡
    人行橫道交通燈元胞
    “機動車禮讓斑馬線”,到底怎么“讓”
    新傳奇(2021年32期)2021-09-22 09:59:20
    基于單片機的交通燈控制系統(tǒng)設計
    電子測試(2018年23期)2018-12-29 11:11:34
    基于元胞自動機下的交通事故路段仿真
    智富時代(2018年5期)2018-07-18 17:52:04
    放學了
    小小藝術家(2018年5期)2018-07-01 17:05:42
    一直飄紅的交通燈
    基于單片機的LED模擬交通燈設計
    電子制作(2016年21期)2016-05-17 03:52:44
    基于元胞數(shù)據(jù)的多維數(shù)據(jù)傳遞機制
    北京測繪(2016年2期)2016-01-24 02:28:28
    世界各地創(chuàng)意人行橫道
    基于AIS的航道移動瓶頸元胞自動機模型
    中國航海(2014年1期)2014-05-09 07:54:25
    基于protues仿真的交通燈即時控制的一種實現(xiàn)
    河南科技(2014年11期)2014-02-27 14:09:57
    成年动漫av网址| 精品人妻在线不人妻| 性高湖久久久久久久久免费观看| 69精品国产乱码久久久| 在线观看国产h片| 精品国产国语对白av| 黄色视频在线播放观看不卡| 成人国产麻豆网| 少妇精品久久久久久久| 欧美丝袜亚洲另类| 国产亚洲av片在线观看秒播厂| 国产成人精品一,二区| 欧美激情 高清一区二区三区| 免费av不卡在线播放| 视频在线观看一区二区三区| av片东京热男人的天堂| 久久午夜综合久久蜜桃| 日韩成人av中文字幕在线观看| 91精品国产国语对白视频| 久久国内精品自在自线图片| 日本爱情动作片www.在线观看| 国产成人av激情在线播放| 最近中文字幕高清免费大全6| 欧美xxxx性猛交bbbb| 免费人成在线观看视频色| 国产精品久久久久久久电影| 国产国语露脸激情在线看| 黄色怎么调成土黄色| 亚洲精品第二区| 秋霞伦理黄片| av网站免费在线观看视频| 七月丁香在线播放| 午夜91福利影院| 国产精品久久久久久久久免| 国产一区二区激情短视频 | 亚洲国产精品一区二区三区在线| 777米奇影视久久| av电影中文网址| 十八禁高潮呻吟视频| 看十八女毛片水多多多| 色吧在线观看| 中文乱码字字幕精品一区二区三区| 又黄又爽又刺激的免费视频.| 最近中文字幕2019免费版| 精品一区二区三区四区五区乱码 | 久热这里只有精品99| 男的添女的下面高潮视频| 久久人人爽人人爽人人片va| 久久久久国产精品人妻一区二区| 免费在线观看完整版高清| 超色免费av| 国产视频首页在线观看| 免费看av在线观看网站| 性高湖久久久久久久久免费观看| 亚洲av男天堂| 99久久精品国产国产毛片| 精品亚洲成国产av| 日本黄大片高清| 亚洲av福利一区| 欧美 亚洲 国产 日韩一| 久久久a久久爽久久v久久| 久久久久久久久久人人人人人人| 黑人高潮一二区| 人妻人人澡人人爽人人| 国产精品麻豆人妻色哟哟久久| 国产有黄有色有爽视频| 一级a做视频免费观看| 国产精品一区www在线观看| 久久婷婷青草| av在线观看视频网站免费| 国产高清国产精品国产三级| 久久久久国产精品人妻一区二区| 日本黄色日本黄色录像| 亚洲美女视频黄频| 国内精品宾馆在线| 18+在线观看网站| 婷婷色麻豆天堂久久| av一本久久久久| 国产日韩欧美视频二区| 国产欧美亚洲国产| 菩萨蛮人人尽说江南好唐韦庄| videossex国产| 国产在线一区二区三区精| 亚洲经典国产精华液单| 国产男人的电影天堂91| 青春草国产在线视频| 国产成人一区二区在线| 最新中文字幕久久久久| 99久国产av精品国产电影| 夜夜爽夜夜爽视频| 蜜臀久久99精品久久宅男| 国产日韩一区二区三区精品不卡| 久久人人97超碰香蕉20202| 色5月婷婷丁香| 寂寞人妻少妇视频99o| 99香蕉大伊视频| av免费在线看不卡| 午夜激情av网站| 免费看不卡的av| 成人手机av| 精品少妇久久久久久888优播| 丝袜在线中文字幕| av有码第一页| 日韩不卡一区二区三区视频在线| 中文字幕精品免费在线观看视频 | 欧美精品av麻豆av| 免费看av在线观看网站| 欧美最新免费一区二区三区| 国产男女内射视频| 麻豆乱淫一区二区| 亚洲成色77777| 成人国产av品久久久| 777米奇影视久久| 熟女人妻精品中文字幕| 黑人巨大精品欧美一区二区蜜桃 | 成人二区视频| 欧美3d第一页| 日韩制服丝袜自拍偷拍| 亚洲欧美成人精品一区二区| 最近中文字幕2019免费版| 国产精品成人在线| 精品午夜福利在线看| 亚洲色图综合在线观看| 五月天丁香电影| 男人舔女人的私密视频| 欧美变态另类bdsm刘玥| 免费av不卡在线播放| 久久99蜜桃精品久久| 少妇的丰满在线观看| 两个人看的免费小视频| 亚洲成人一二三区av| 成人免费观看视频高清| 日韩三级伦理在线观看| 欧美3d第一页| 久久久久久久久久人人人人人人| 久久久久人妻精品一区果冻| 在线亚洲精品国产二区图片欧美| 亚洲一码二码三码区别大吗| 免费观看性生交大片5| 亚洲国产精品国产精品| 免费高清在线观看视频在线观看| 午夜福利视频精品| 国产av国产精品国产| 18禁观看日本| 22中文网久久字幕| 高清不卡的av网站| h视频一区二区三区| 久久99蜜桃精品久久| 精品少妇内射三级| 久久这里有精品视频免费| 一级,二级,三级黄色视频| xxxhd国产人妻xxx| 性色av一级| 国产精品.久久久| 女人精品久久久久毛片| 激情视频va一区二区三区| 性高湖久久久久久久久免费观看| 精品午夜福利在线看| 麻豆乱淫一区二区| 熟女电影av网| 激情五月婷婷亚洲| 乱码一卡2卡4卡精品| 99视频精品全部免费 在线| 777米奇影视久久| 天堂俺去俺来也www色官网| 岛国毛片在线播放| 97在线视频观看| 国产一级毛片在线| 国产永久视频网站| 亚洲高清免费不卡视频| 久久狼人影院| 男人操女人黄网站| 国产 精品1| 国产探花极品一区二区| 亚洲欧洲日产国产| 中文字幕亚洲精品专区| 22中文网久久字幕| 久久97久久精品| 国产色爽女视频免费观看| 亚洲精品视频女| 亚洲色图综合在线观看| 考比视频在线观看| 一本色道久久久久久精品综合| 99热全是精品| 天堂中文最新版在线下载| 狂野欧美激情性xxxx在线观看| 欧美日韩亚洲高清精品| 天天躁夜夜躁狠狠久久av| 久久人人爽av亚洲精品天堂| 男人添女人高潮全过程视频| av免费观看日本| 春色校园在线视频观看| 欧美日本中文国产一区发布| 日韩三级伦理在线观看| 欧美3d第一页| 永久网站在线| 99久久综合免费| h视频一区二区三区| 日本午夜av视频| 日本欧美国产在线视频| av有码第一页| 久久热在线av| 青春草视频在线免费观看| 亚洲天堂av无毛| 亚洲四区av| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美日韩视频精品一区| 欧美日韩av久久| 国产在线免费精品| 狂野欧美激情性xxxx在线观看| 乱人伦中国视频| 91午夜精品亚洲一区二区三区| 午夜福利视频精品| 赤兔流量卡办理| 亚洲伊人久久精品综合| 亚洲精品中文字幕在线视频| 成年美女黄网站色视频大全免费| 美女福利国产在线| 亚洲,欧美,日韩| 精品午夜福利在线看| 亚洲精品456在线播放app| 久久人妻熟女aⅴ| 精品少妇久久久久久888优播| 亚洲第一av免费看| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品第二区| 久久国产精品大桥未久av| 久久ye,这里只有精品| 亚洲伊人色综图| 超色免费av| 99久久精品国产国产毛片| 精品少妇黑人巨大在线播放| 97在线视频观看| 十八禁高潮呻吟视频| 国产黄频视频在线观看| 亚洲欧美中文字幕日韩二区| 国产av码专区亚洲av| 青春草亚洲视频在线观看| 一个人免费看片子| 日韩制服骚丝袜av| 18禁动态无遮挡网站| 男女边摸边吃奶| av国产久精品久网站免费入址| 国产成人午夜福利电影在线观看| 国产亚洲欧美精品永久| 乱码一卡2卡4卡精品| 国产精品一国产av| 久久久久国产精品人妻一区二区| 永久免费av网站大全| 亚洲精品av麻豆狂野| 少妇的逼水好多| 亚洲av男天堂| 国产精品久久久久久久久免| 在线观看免费日韩欧美大片| 美女国产高潮福利片在线看| 精品熟女少妇av免费看| 咕卡用的链子| av有码第一页| 最近中文字幕高清免费大全6| 秋霞伦理黄片| 大片电影免费在线观看免费| 久久久a久久爽久久v久久| 欧美人与性动交α欧美精品济南到 | 在线看a的网站| 欧美另类一区| 97精品久久久久久久久久精品| 久久久久久人妻| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产精品国产三级国产av玫瑰| 飞空精品影院首页| 欧美精品一区二区免费开放| 一边摸一边做爽爽视频免费| 久久久久久久久久久久大奶| 国产精品久久久久成人av| 五月开心婷婷网| 欧美性感艳星| 18禁国产床啪视频网站| 91久久精品国产一区二区三区| av在线老鸭窝| 一级毛片 在线播放| 亚洲av综合色区一区| av网站免费在线观看视频| 美国免费a级毛片| 久久鲁丝午夜福利片| 免费av中文字幕在线| 91成人精品电影| 国产极品粉嫩免费观看在线| 欧美 日韩 精品 国产| 最新的欧美精品一区二区| 在线 av 中文字幕| 精品国产露脸久久av麻豆| 免费观看无遮挡的男女| 欧美激情国产日韩精品一区| 久久人人爽av亚洲精品天堂| 中国美白少妇内射xxxbb| 午夜激情久久久久久久| 国产xxxxx性猛交| 美女国产视频在线观看| 亚洲精品一区蜜桃| 一区二区av电影网| 国产69精品久久久久777片| 日韩一本色道免费dvd| 日本爱情动作片www.在线观看| 久久婷婷青草| 国产精品免费大片| 美女中出高潮动态图| 亚洲av.av天堂| 午夜精品国产一区二区电影| 69精品国产乱码久久久| 免费看光身美女| 啦啦啦在线观看免费高清www| 老熟女久久久| 国产精品秋霞免费鲁丝片| 超碰97精品在线观看| 日韩av在线免费看完整版不卡| 中文字幕精品免费在线观看视频 | 国产 精品1| 制服丝袜香蕉在线| 热99国产精品久久久久久7| 最近中文字幕2019免费版| 一二三四中文在线观看免费高清| 亚洲精品美女久久av网站| av天堂久久9| 日韩精品免费视频一区二区三区 | 欧美成人午夜免费资源| 亚洲美女视频黄频| 国产在线一区二区三区精| 中文字幕人妻熟女乱码| 亚洲精品国产av蜜桃| 成人黄色视频免费在线看| 韩国av在线不卡| 一级,二级,三级黄色视频| 美女国产高潮福利片在线看| 中文天堂在线官网| a 毛片基地| 天堂中文最新版在线下载| 日韩欧美精品免费久久| 久久久久久久久久人人人人人人| 午夜老司机福利剧场| 伦理电影大哥的女人| 校园人妻丝袜中文字幕| 黄色配什么色好看| 卡戴珊不雅视频在线播放| 久久久久精品人妻al黑| 高清不卡的av网站| 亚洲经典国产精华液单| 最后的刺客免费高清国语| 久久久久久久久久久免费av| 日日摸夜夜添夜夜爱| 久久久久国产网址| 精品少妇内射三级| 九九在线视频观看精品| 婷婷色综合www| 久久久久久久久久人人人人人人| 国语对白做爰xxxⅹ性视频网站| 看十八女毛片水多多多| 欧美日韩国产mv在线观看视频| 性高湖久久久久久久久免费观看| 国产免费一区二区三区四区乱码| 人体艺术视频欧美日本| 大香蕉97超碰在线| 欧美精品一区二区大全| 国产成人aa在线观看| 国产精品 国内视频| 国产精品99久久99久久久不卡 | 国产av码专区亚洲av| 日韩欧美精品免费久久| 午夜福利在线观看免费完整高清在| 中文字幕精品免费在线观看视频 | 日韩一区二区三区影片| 香蕉国产在线看| 亚洲精品第二区| 18在线观看网站| 少妇熟女欧美另类| 亚洲色图 男人天堂 中文字幕 | 精品人妻偷拍中文字幕| 久久人人爽人人片av| 亚洲激情五月婷婷啪啪| 国产亚洲一区二区精品| 国产熟女午夜一区二区三区| 狠狠精品人妻久久久久久综合| 日韩中字成人| 久久久久久人妻| 一区二区av电影网| av在线老鸭窝| 亚洲欧洲日产国产| 欧美激情 高清一区二区三区| 美女xxoo啪啪120秒动态图| 国产不卡av网站在线观看| 狂野欧美激情性xxxx在线观看| a级片在线免费高清观看视频| 中文字幕最新亚洲高清| 人妻一区二区av| 观看美女的网站| 成年人免费黄色播放视频| 啦啦啦视频在线资源免费观看| 国产在线免费精品| 毛片一级片免费看久久久久| 你懂的网址亚洲精品在线观看| 婷婷色综合www| 亚洲国产精品一区二区三区在线| 51国产日韩欧美| 亚洲av中文av极速乱| 黄网站色视频无遮挡免费观看| 91国产中文字幕| 美女视频免费永久观看网站| 国产国语露脸激情在线看| 亚洲精品一区蜜桃| 777米奇影视久久| 精品一品国产午夜福利视频| 精品久久国产蜜桃| 国产成人一区二区在线| 91精品伊人久久大香线蕉| 亚洲,欧美精品.| 自拍欧美九色日韩亚洲蝌蚪91| 老女人水多毛片| 爱豆传媒免费全集在线观看| 亚洲伊人色综图| 中文欧美无线码| 久久久久久人人人人人| 精品一区二区免费观看| 欧美97在线视频| 国产精品三级大全| 国产成人一区二区在线| 精品少妇内射三级| 大话2 男鬼变身卡| 久久久久国产精品人妻一区二区| 波多野结衣一区麻豆| 九九爱精品视频在线观看| 十分钟在线观看高清视频www| 欧美97在线视频| 国产精品欧美亚洲77777| 亚洲,欧美精品.| 国精品久久久久久国模美| 国产男女超爽视频在线观看| 国产精品熟女久久久久浪| a级毛片在线看网站| 性色av一级| 国产欧美日韩综合在线一区二区| 国产色爽女视频免费观看| 免费久久久久久久精品成人欧美视频 | 精品福利永久在线观看| 欧美丝袜亚洲另类| 亚洲国产av新网站| 日韩视频在线欧美| 黑人欧美特级aaaaaa片| 亚洲av电影在线进入| 国产精品熟女久久久久浪| 亚洲激情五月婷婷啪啪| 久久久久久人人人人人| 一本色道久久久久久精品综合| 国产色爽女视频免费观看| 色婷婷av一区二区三区视频| 欧美日韩亚洲高清精品| 高清黄色对白视频在线免费看| 亚洲色图综合在线观看| kizo精华| 久久精品久久精品一区二区三区| 国产av精品麻豆| 国产片特级美女逼逼视频| 国产极品天堂在线| 久久人人97超碰香蕉20202| 爱豆传媒免费全集在线观看| 蜜桃在线观看..| 免费观看无遮挡的男女| 制服丝袜香蕉在线| 韩国高清视频一区二区三区| 精品久久久久久电影网| 久久精品aⅴ一区二区三区四区 | 丰满少妇做爰视频| 男女啪啪激烈高潮av片| 性色avwww在线观看| 国产1区2区3区精品| 欧美变态另类bdsm刘玥| 人人妻人人澡人人看| 国产乱人偷精品视频| 国产成人免费观看mmmm| 国产精品久久久av美女十八| 精品国产一区二区久久| 日韩制服丝袜自拍偷拍| 日日啪夜夜爽| 五月伊人婷婷丁香| 国产亚洲一区二区精品| 亚洲精品,欧美精品| 国产乱来视频区| 国产成人精品久久久久久| 国产免费一区二区三区四区乱码| 亚洲成色77777| 精品国产国语对白av| 国产白丝娇喘喷水9色精品| 国产1区2区3区精品| 伦理电影大哥的女人| 久久99热这里只频精品6学生| 一区二区日韩欧美中文字幕 | 国产高清国产精品国产三级| 久久久精品免费免费高清| 一级a做视频免费观看| 久久精品熟女亚洲av麻豆精品| 日韩三级伦理在线观看| 免费黄频网站在线观看国产| av免费在线看不卡| 免费黄频网站在线观看国产| 99久久中文字幕三级久久日本| 亚洲美女搞黄在线观看| 亚洲欧洲国产日韩| 自线自在国产av| 草草在线视频免费看| 亚洲av在线观看美女高潮| 日本wwww免费看| 97在线人人人人妻| 色吧在线观看| 国产一区有黄有色的免费视频| 亚洲av.av天堂| 91成人精品电影| 久久久久视频综合| 日本欧美视频一区| 日本av手机在线免费观看| 美女脱内裤让男人舔精品视频| 97精品久久久久久久久久精品| 日韩精品有码人妻一区| 久久精品aⅴ一区二区三区四区 | 色婷婷久久久亚洲欧美| av免费观看日本| 国产日韩欧美在线精品| 最近中文字幕高清免费大全6| 亚洲欧美精品自产自拍| 黄色 视频免费看| 精品午夜福利在线看| 在线精品无人区一区二区三| 中文字幕精品免费在线观看视频 | 精品人妻一区二区三区麻豆| 男女国产视频网站| 亚洲成人av在线免费| 国产黄色视频一区二区在线观看| 草草在线视频免费看| 少妇的逼水好多| 嫩草影院入口| 亚洲精品日韩在线中文字幕| 18禁观看日本| 男女高潮啪啪啪动态图| 久久国产精品大桥未久av| 99热全是精品| 成人国产av品久久久| 亚洲国产欧美日韩在线播放| 大码成人一级视频| 精品人妻熟女毛片av久久网站| 高清欧美精品videossex| 91aial.com中文字幕在线观看| 少妇精品久久久久久久| 亚洲精品久久午夜乱码| 纵有疾风起免费观看全集完整版| 欧美 日韩 精品 国产| 2018国产大陆天天弄谢| 岛国毛片在线播放| 51国产日韩欧美| 日本vs欧美在线观看视频| 中国美白少妇内射xxxbb| 高清欧美精品videossex| 国产 一区精品| 国产又爽黄色视频| 久久青草综合色| 男的添女的下面高潮视频| 高清毛片免费看| 国产精品一区二区在线不卡| 99国产综合亚洲精品| 国产亚洲一区二区精品| av又黄又爽大尺度在线免费看| 欧美激情国产日韩精品一区| 国产老妇伦熟女老妇高清| 看非洲黑人一级黄片| 欧美变态另类bdsm刘玥| 国语对白做爰xxxⅹ性视频网站| 亚洲美女黄色视频免费看| 日韩视频在线欧美| 午夜福利视频精品| 2022亚洲国产成人精品| 又黄又爽又刺激的免费视频.| 久久精品久久精品一区二区三区| 母亲3免费完整高清在线观看 | 精品人妻在线不人妻| 9191精品国产免费久久| 一区二区三区精品91| 丝袜喷水一区| av国产精品久久久久影院| 久久免费观看电影| 91成人精品电影| 久久青草综合色| 久久免费观看电影| 在线观看人妻少妇| 午夜免费鲁丝| 国产精品 国内视频| 国产亚洲精品第一综合不卡 | 免费女性裸体啪啪无遮挡网站| 欧美成人精品欧美一级黄| 久久久精品94久久精品| 插逼视频在线观看| 晚上一个人看的免费电影| www.熟女人妻精品国产 | 侵犯人妻中文字幕一二三四区| 中国三级夫妇交换| 九草在线视频观看| 99热6这里只有精品| 国产一区亚洲一区在线观看| 九草在线视频观看| 国产日韩欧美视频二区| 国产男人的电影天堂91| 人人妻人人澡人人爽人人夜夜| 国产69精品久久久久777片| 日韩av免费高清视频| 精品国产一区二区三区久久久樱花| 丝袜在线中文字幕| 亚洲精品国产色婷婷电影| 大话2 男鬼变身卡| 久久精品国产a三级三级三级| 人人妻人人添人人爽欧美一区卜|