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

    步行通道內(nèi)行人流拉鏈現(xiàn)象的生成機(jī)理與仿真研究*

    2021-03-26 08:43:54武鑫森岳昊劉秋梅張旭邵春福
    物理學(xué)報(bào) 2021年6期
    關(guān)鍵詞:影響模型

    武鑫森 岳昊? 劉秋梅 張旭 邵春福

    1) (北京交通大學(xué), 綜合交通運(yùn)輸大數(shù)據(jù)應(yīng)用技術(shù)交通運(yùn)輸行業(yè)重點(diǎn)實(shí)驗(yàn)室, 北京 100044)

    2) (河南工業(yè)大學(xué)土木建筑學(xué)院, 鄭州 450001)

    1 引 言

    隨著社會(huì)經(jīng)濟(jì)發(fā)展, 公共空間的人群集散問題越來越受到人們關(guān)注, 對(duì)行人動(dòng)力學(xué)的研究也日益增多.揭示行人動(dòng)力學(xué)的運(yùn)動(dòng)特征和相關(guān)規(guī)律, 對(duì)合理指導(dǎo)步行設(shè)施設(shè)計(jì)[1], 科學(xué)制定人群組織和疏散策略[2], 有效誘導(dǎo)設(shè)施內(nèi)行人流移動(dòng), 避免和緩解行人流擁堵等具有重要意義.基于行人微觀行為的仿真研究, 可以定性定量地分析行人流特性[3],再現(xiàn)人群宏觀現(xiàn)象, 避免緊急情況下潛在的安全問題, 以及節(jié)約人力成本、增加效率[4].

    近年來, 基于微觀行為的行人流仿真模型主要包括: 離散仿真模型和連續(xù)仿真模型.離散仿真模型, 即仿真時(shí)間和空間離散, 如元胞自動(dòng)機(jī)模型[5-7]、領(lǐng)域模型[8]、格子氣模型[9]、動(dòng)態(tài)參數(shù)模型[10]等;連續(xù)仿真模型, 即仿真時(shí)間離散, 而仿真空間連續(xù),如社會(huì)力模型[11-23]、離心力模型[24,25]等.從仿真方法的角度, 可以分為基于力和基于啟發(fā)式方法的仿真模型[26-29].研究對(duì)象主要包括步行通道或疏散空間等步行設(shè)施內(nèi), 單向或雙向行人通過流[12]、疏散流[13].疏散場景包括: 行人視線受影響[11], 多安全出口[14]或存在空間障礙物[15]等.仿真研究再現(xiàn)了許多行人流自組織現(xiàn)象, 例如: 自動(dòng)渠化[11,12]、通道形成[11,12]、成拱[11,12]、快即是慢[7]、避讓三角區(qū)[16]、拉鏈效應(yīng)[30,31]等.

    在基于力仿真模型方面, 經(jīng)典模型是社會(huì)力仿真模型, 提出假想社會(huì)力的概念來描述行人之間的相互作用, 行人運(yùn)動(dòng)被定量描述為受自驅(qū)動(dòng)力、排斥力和接觸摩擦力等合力的影響[11].傳統(tǒng)社會(huì)力模型可以克服元胞自動(dòng)機(jī)模型中行人移動(dòng)方向和速度固化等問題, 模擬再現(xiàn)更平滑的行人移動(dòng)軌跡, 更好地定量描述行人之間的動(dòng)態(tài)特性和相互作用; 盡管在模擬行人微觀行為與再現(xiàn)宏觀現(xiàn)象方面有較強(qiáng)的優(yōu)勢(shì), 但仍然有不足, 例如: 參數(shù)復(fù)雜難以準(zhǔn)確合理校準(zhǔn)[18], 會(huì)出現(xiàn)行人抖動(dòng)和重疊現(xiàn)象[19-22]; 仿真少量行人會(huì)出現(xiàn)不合理的行人死鎖現(xiàn)象[16,28]; 缺乏主動(dòng)繞行其他行人的機(jī)制[16,28]等.為克服社會(huì)力模型的缺點(diǎn), 從多種角度進(jìn)行修改完善, 包括: 自驅(qū)動(dòng)力的改進(jìn)[13]和排斥力的改進(jìn)[17];引入其他參數(shù)和參數(shù)標(biāo)定[17]; 構(gòu)建共軛力使行人可以主動(dòng)繞行[16]; 為避免行人重疊引入行人停止機(jī)制[19,20]、引入基于顯式數(shù)值積分的算法[21]、引入減速避讓機(jī)制[22]等.

    在基于啟發(fā)式方法的仿真模型方面[27-29], 很多研究將啟發(fā)式方法整合到行人建模中, 處理復(fù)雜情況并預(yù)測(cè)更平滑的行人軌跡.Moussa?d 等[29]認(rèn)為行人行為模式改變是通過兩種簡單的認(rèn)知過程來完成, 并引入了行人的視覺信息, 提出一種簡單的行為啟發(fā)式模型來捕捉人群行為中的潛在規(guī)律;同時(shí), 局部密度可用于反映行人運(yùn)動(dòng)的舒適程度并作為引發(fā)行人運(yùn)動(dòng)的重要參數(shù)[29].因?yàn)閂oronoi圖特有的幾何特性, 有學(xué)者將Voronoi 圖的概念引入到行人動(dòng)力學(xué)領(lǐng)域, 用于確定行人的個(gè)人空間、領(lǐng)域、安全距離、鄰居, 行人局部密度及基于密度的最佳速度等[27,28].在Voronoi 圖中, 每個(gè)粒子分別位于各自的Voronoi 元胞中; Voronoi 邊位于相鄰粒子連線的中垂線上, 由兩個(gè)相鄰的Voronoi 元胞共享; Voronoi 點(diǎn)在非退化條件下關(guān)聯(lián)3 個(gè)Voronoi 元胞, 是3 個(gè)Voronoi 邊的交叉點(diǎn).在行人仿真過程中, 與行人共享Voronoi 邊的其他行人可視為行人的鄰居, Voronoi 元胞的形狀會(huì)隨著行人位置的變化而實(shí)時(shí)變化.行人的Voronoi 元胞看起來像扭曲的傳統(tǒng)元胞自動(dòng)機(jī), 但又與其不同;Voronoi 圖中的元胞是不規(guī)則的, 且不存在空元胞.

    在行人間相互影響和行人關(guān)注領(lǐng)域研究方面,一般常識(shí)認(rèn)為在行人移動(dòng)過程中行人周圍不同方向和距離的行人會(huì)對(duì)自身產(chǎn)生不同的影響.例如多個(gè)相互不認(rèn)識(shí)的行人一起移動(dòng), 后方行人必須要考慮前方行人的狀態(tài)來決定是否跟隨、避讓或者超越, 而前方行人卻不需要顧及后方行人.同時(shí),Johansson 等[23]基于視頻中的行人運(yùn)動(dòng)軌跡數(shù)據(jù),發(fā)現(xiàn)行人與周圍不同方向行人間的相互作用也表現(xiàn)出各向異性的特征, 行人只受其前方約180°空間范圍內(nèi)其他行人的影響, 而幾乎不受身后行人的影響.通過在社會(huì)力模型的框架下詳細(xì)分析行人間的相互作用, 發(fā)現(xiàn)不同方向行人對(duì)目標(biāo)行人的影響是不同的, 處于正前方的行人對(duì)目標(biāo)行人的作用力最強(qiáng), 處于兩側(cè)的行人對(duì)目標(biāo)行人的作用力最弱.Yu 等[25]在離心力模型中也提到了行人視野范圍的問題, 即行人一般只關(guān)注視野前方180°范圍內(nèi)的情況, 對(duì)兩側(cè)顧及得比較少, 對(duì)身后就幾乎不會(huì)去顧及.Porter 等[32]提出了IM (integrated model)模型框架, 集成了社會(huì)力模型、行為啟發(fā)式模型及材料科學(xué)理論, 充分發(fā)揮了每個(gè)模型的優(yōu)勢(shì), 將材料科學(xué)中“只考慮直接相鄰的分子, 就可以很好地模擬分子間相互作用”的理論應(yīng)用到行人移動(dòng)的相互影響中, 只考慮行人視野范圍中周圍行人的影響, 而不是考慮周圍所有行人, 從而大大簡化了運(yùn)算量.

    在日常觀察中, 雖然作為行人流自組織現(xiàn)象之一的拉鏈現(xiàn)象常被發(fā)現(xiàn), 同時(shí)拉鏈效應(yīng)也會(huì)對(duì)行人的移動(dòng)行為產(chǎn)生影響; 但目前學(xué)者主要是針對(duì)拉鏈現(xiàn)象的行人實(shí)驗(yàn)進(jìn)行了研究, 而缺乏對(duì)拉鏈效應(yīng)形成機(jī)理的研究.在行人實(shí)驗(yàn)中, Moussa?d 等[24],Hoogendoorn 等[30]發(fā)現(xiàn)行人在移動(dòng)過程中有側(cè)向行走的傾向; Hoogendoorn 等[30]發(fā)現(xiàn)通道內(nèi)單向行人流的“拉鏈”效應(yīng)可提高通道的空間利用率和行人的步行速度; Seyfried 等[31]發(fā)現(xiàn)拉鏈效應(yīng)會(huì)對(duì)通道的通行能力產(chǎn)生影響; Ren 等[33]觀察到行人的隊(duì)列層數(shù)會(huì)隨著通道寬度的增加而增加.

    對(duì)拉鏈效應(yīng)的研究, 有利于研究行人拉鏈現(xiàn)象的生成機(jī)理, 以及與步行設(shè)施的關(guān)系, 從而指導(dǎo)步行設(shè)施的規(guī)劃設(shè)計(jì)和人流組織.因此, 本文從拉鏈效應(yīng)的表現(xiàn)形式, 生成拉鏈效應(yīng)的影響因素, 拉鏈效應(yīng)的偏移機(jī)制等方面探究拉鏈效應(yīng)的生成機(jī)理.基于Voronoi 圖, 考慮不同距離和方向周圍行人對(duì)行人移動(dòng)速度的影響, 構(gòu)建速度修正模型, 考慮行人是否有偏移傾向的主觀意愿, 引入偏移規(guī)則, 模擬行人的拉鏈現(xiàn)象.以行人的視野關(guān)注和視野遮擋描述影響行人移動(dòng)過程中偏移行為的因素, 以行人局部密度描述行人步行舒適的情況, 引入拉鏈敏感系數(shù)描述行人客觀偏移的意愿程度, 提出單個(gè)行人側(cè)向偏移的機(jī)制, 獲得行人最佳的偏移位置.仿真研究行人的拉鏈層數(shù)與通道寬度的關(guān)系, 然后, 通過與其他實(shí)證數(shù)據(jù)對(duì)比, 驗(yàn)證模型的可靠性; 研究不同拉鏈效應(yīng)主觀偏移意愿占比與行人平均速度、平均局部密度的關(guān)系等.

    2 拉鏈效應(yīng)機(jī)理

    2.1 拉鏈現(xiàn)象

    行人流自組織現(xiàn)象是指行人流內(nèi)部個(gè)體行人之間在無外界力量指導(dǎo)下, 出于自身因素和作用,自發(fā)組織和形成的行人流宏觀現(xiàn)象.自組織現(xiàn)象之一的拉鏈效應(yīng)[30]是指通道內(nèi)多條單向行人流隊(duì)列之間的相互重疊現(xiàn)象.在通道內(nèi)運(yùn)動(dòng)的單向行人流, 把每一列行人看作一層, 會(huì)自發(fā)產(chǎn)生層或小道,不同層間的行人相互交錯(cuò), 各層行人均占用相鄰層的空間, 形成“拉鏈”現(xiàn)象, 如圖1 所示.該現(xiàn)象的產(chǎn)生主要是由于行人需要盡可能地利用通道中的可行空間, 在確保自身行走空間的同時(shí)提高自己的視野.

    圖1 瓶頸中的行人拉鏈現(xiàn)象[30]Fig.1.Zipper phenomena of pedestrians in a bottleneck[30].

    在拉鏈效應(yīng)行人實(shí)驗(yàn)中[33], 當(dāng)充足行人以正常速度通過不同寬度通道時(shí), 可以觀察到拉鏈層數(shù)會(huì)隨通道寬度增加而增加, 如圖2 所示.當(dāng)通道寬度B 為0.5—0.7 m 時(shí), 因通道寬度較窄, 基本形成1 層行人; 當(dāng)通道寬度B 為0.8—1.1 m 時(shí), 隨通道寬度增加, 可以觀察到形成2—3 層行人; 當(dāng)通道寬度B 為1.2—1.8 m 時(shí), 因通道寬度較寬, 可以觀察到形成3—4 層行人.

    2.2 表現(xiàn)形式

    行人在移動(dòng)過程中, 手腳有節(jié)奏的依次來回?cái)[動(dòng), 其所占動(dòng)態(tài)空間投影區(qū)域類似圓形.在相同方向行走過程中, 兩相鄰行人組成的行人間基本位置形式可分為并排、斜列、直列3 種, 如圖3 所示.

    通過分析行人實(shí)驗(yàn)[33], 研究通道寬度與行人拉鏈現(xiàn)象之間的關(guān)系, 可以發(fā)現(xiàn)當(dāng)通道寬度僅能容納一個(gè)行人通過時(shí), 行人以直列形式排列, 如圖4(a)所示; 當(dāng)通道寬度不足以容納兩個(gè)行人并排通過時(shí), 行人以直列、斜列形式排列, 如圖4(b)所示; 當(dāng)通道寬度足以容納多個(gè)行人并排通過時(shí),行人之間的并排、斜列、直列多種排列形式均會(huì)出現(xiàn), 如圖4(c)所示.圖4(b)和圖4(c)是典型拉鏈現(xiàn)象的表現(xiàn)形式.

    2.3 視野最佳

    圖2 不同通道寬度下的行人拉鏈現(xiàn)象[33]Fig.2.Zipper phenomenon of pedestrians under different channel widths[33].

    圖3 兩相鄰行人同方向行走形態(tài) (a)行人行走形態(tài); (b)兩行人行走并排形式; (c)兩行人行走斜列形式; (d)兩行人行走直列形式Fig.3.Walking pattern of two adjacent pedestrians in the same direction: (a) The real single pedestrian walking pattern; (b) two adjacent pedestrians walking side by side; (c) two adjacent pedestrians walking in a diagonal pattern; (d) two adjacent pedestrians walking in a straight line.

    圖4 行人拉鏈現(xiàn)象表現(xiàn)形式 (a)單列行人截圖[33]與行人直列形式對(duì)比; (b)兩列行人截圖[33]與行人直列、斜列形式對(duì)比;(c)多列行人截圖[33]與行人多種排列形式對(duì)比Fig.4.Manifestation of pedestrian zipper phenomenon: (a) The comparison of single-column pedestrian screenshot[33] with pedestrians in-line pattern; (b) the comparison of two-column pedestrian screenshot[33] with pedestrians in-line and diagonal patterns;(c) the comparison of multi-column pedestrian screenshot[33] with multiple arrangement patterns of pedestrians.

    行人行走時(shí)一般目視前方, 在視野后方的行人基本不會(huì)對(duì)前方行人的行為決策產(chǎn)生影響.在行人有效視野范圍內(nèi), 也不是所有角度和距離都會(huì)對(duì)行人產(chǎn)生同樣的影響.顯然正前方是行人最為關(guān)注的視野區(qū)域, 對(duì)行人的行為決策影響最大; 而從中間向兩側(cè)的視野區(qū)域?qū)π腥诵袨闆Q策的影響程度迅速衰減.以圖5 為例, 移動(dòng)行人 P1的視野被行人 P2,P3, P4阻 擋, 行 人 P1的 可 見 視 野 區(qū) 域 是 Z2, Z4,Z6, 不可見視野區(qū)域是 Z1, Z3, Z5.顯然 Z3區(qū)域阻擋了行人 P1正前方的主要觀察視野和道路, 對(duì)行人 P1行為決策的影響最大, 是促使行人 P1側(cè)向偏離的主要誘因; Z5區(qū)域雖然對(duì)行人 P1觀察周圍視野產(chǎn)生一定遮擋, 但不在行人 P1前進(jìn)道路上, 對(duì)行人 P1行為決策的影響較弱; Z1區(qū)域處于行人P1觀察視野的邊緣區(qū)域, 對(duì)行人 P1行為決策的影響最弱, 幾乎不會(huì)引起行人 P1的特別關(guān)注.

    考慮視野對(duì)行人行為決策的影響, 以行人視野關(guān)注和視野遮擋描述影響行人側(cè)向偏離的因素.視野關(guān)注表示行人對(duì)不同視野區(qū)域的關(guān)注程度.同時(shí), 不同視野區(qū)域?qū)π腥诵袨闆Q策的影響程度也不同.視野遮擋表示在行人有效視野區(qū)域內(nèi), 被遮擋視野區(qū)域面積占行人有效視野區(qū)域面積的比重.正前方區(qū)域是行人視野最為關(guān)注和對(duì)行人視野遮擋影響最大的區(qū)域.

    圖5 有效視野區(qū)域內(nèi)的行人視野關(guān)注和遮擋(橙色遮擋區(qū)域是目標(biāo)行人最為關(guān)注和對(duì)目標(biāo)行人影響最大的區(qū)域,藍(lán)色遮擋區(qū)域?qū)δ繕?biāo)行人的影響較弱)Fig.5.Pedestrian vision attention and occlusion in the effective vision area (the orange shaded area is the area that the target pedestrian pay the most attention to and has the greatest influence on the target pedestrian, while the blue shaded area has a relatively weak influence on the target pedestrian).

    此外, 即使是同一視野方向, 不同距離行人對(duì)目標(biāo)行人視野的遮擋影響程度也不同.距離越近,對(duì)目標(biāo)行人視野遮擋越大; 距離越遠(yuǎn), 對(duì)目標(biāo)行人視野遮擋度越小.以圖6 為例, 行人 P3和 P5相對(duì)目標(biāo)行人 P1在同一視野方向, 與行人 P1的距離分別是 l13和 l15, 對(duì)行人 P1視野的遮擋角度分別是σ13和 σ15, 遮 擋 范 圍 分 別 是 Z5和顯 然 l15>l13,σ15<σ13, 但

    圖6 不同距離行人對(duì)目標(biāo)行人視野的遮擋程度Fig.6.Occlusion degree of the target pedestrian’s vision by surrounding pedestrians at different distances.

    2.4 步行舒適

    當(dāng)通道寬度足以容納多個(gè)行人通過時(shí), 在觀察行人拉鏈現(xiàn)象過程中, 可以引入行人局部密度, 描述行人的步行舒適度.在Voronoi 圖中, 假設(shè)行人Pi所對(duì)應(yīng)Voronoi 元胞的面積是 ai, 則其局部密度為 ρi=1/ai[27].當(dāng)行人局部密度較低時(shí), 行人擁有較大的個(gè)人空間和較好的步行舒適度, 行人之間會(huì)出現(xiàn)并排、斜列、直列等形式, 且行人的步行速度不受影響, 行人以自由流速度移動(dòng), 如圖7(a)所示;當(dāng)行人局部密度較高時(shí), 行人的個(gè)人空間較小, 行人需要充分利用步行設(shè)施空間移動(dòng), 行人的步行舒適度會(huì)受到影響, 且會(huì)影響和限制后面行人的步行速度, 行人之間會(huì)出現(xiàn)各種排列形式, 如圖7(b)所示.

    圖7 不同行人局部密度下的拉鏈現(xiàn)象示意(綠色、黃色、橙色分別表示行人局部密度由低到高) (a)行人局部密度較低時(shí); (b)行人局部密度較高時(shí)Fig.7.The zipper phenomenon of pedestrians under different local density (green, yellow and orange respectively represent the local pedestrian density from low to high):(a) When the local pedestrian density is relatively low;(b) when the local pedestrian density is relatively high.

    同時(shí), 行人是否會(huì)進(jìn)行側(cè)向偏移還與自身的主客觀偏離意愿和習(xí)慣有關(guān).以圖8 為例, 當(dāng)通道寬度正好可以允許兩個(gè)行人并排通過時(shí), 若行人的偏離意愿程度較低, 行人之間有可能形成規(guī)則形狀的直列, 如圖8(a)所示; 當(dāng)通道寬度不足以允許兩個(gè)行人并排通過時(shí), 若行人的偏離意愿程度較高, 行人之間會(huì)自發(fā)形成不規(guī)則形狀的拉鏈, 如圖8(b)所示; 當(dāng)通道寬度可以允許多個(gè)行人并排通過時(shí),若行人的偏離意愿程度較高, 行人在潛意識(shí)驅(qū)動(dòng)下會(huì)自發(fā)進(jìn)行側(cè)向偏離, 追求更寬廣的可見視野和步行舒適度, 形成典型的拉鏈現(xiàn)象, 如圖8(c)所示.

    圖8 不同偏離意愿下的拉鏈現(xiàn)象示意 (a)窄通道行人偏離意愿較低時(shí); (b)窄通道行人偏離意愿較高時(shí); (c) 寬通道行人偏離意愿較高時(shí)Fig.8.Zipper phenomenon of pedestrians under different deviation intentions: (a) When the deviation intention of pedestrians in narrow corridors is relatively low; (b) when the deviation intention of pedestrians in narrow corridors is relatively high; (c) when the deviation intention of pedestrians in wide corridors is relatively high.

    圖9 行人拉鏈效應(yīng)的穩(wěn)定形態(tài)[33]Fig.9.Stable form of pedestrian zipper effect[33].

    在行人追求視野最佳和步行舒適, 處于比較穩(wěn)定的拉鏈效應(yīng)狀態(tài)時(shí), 可以發(fā)現(xiàn)目標(biāo)行人與周圍行人會(huì)形成類似蜂巢的排列, 與視野前方行人形成類似等邊三角形的形態(tài), 如圖9 所示[33].處于穩(wěn)定狀態(tài)的行人, 占用較少的空間, 擁有較好的視野、個(gè)人舒適度和移動(dòng)速度.在現(xiàn)有行人仿真中, 也有學(xué)者采用蜂巢形態(tài)作為行人個(gè)體的穩(wěn)定形態(tài)[34].以圖10 為例, 當(dāng)行人之間彼此相切嚙合, 相鄰的3 個(gè)行人形成夾角為 Φ 的等邊三角形, 處于比較理想的穩(wěn)定結(jié)構(gòu)時(shí), 本文認(rèn)為此時(shí)每個(gè)行人在行人流中處于視野最佳和步行舒適的狀態(tài).

    圖10 行人拉鏈效應(yīng)的最佳形態(tài)Fig.10.The best form of pedestrian zipper effect.

    2.5 偏移機(jī)制

    通過觀察, 分析行人追求視野最佳和步行舒適形成拉鏈現(xiàn)象的行為機(jī)制.假設(shè)不考慮行人個(gè)體和主觀偏移意愿的差異, 認(rèn)為行人均有主觀偏移意愿.當(dāng)通道寬度足夠容納多個(gè)行人并肩通過且行人整體密度較大時(shí), 行人會(huì)呈現(xiàn)各種排列形式.其中,斜列形式在一定程度上是行人側(cè)向偏移的結(jié)果, 是行人拉鏈現(xiàn)象的主要表現(xiàn)形式.在行人行走過程中, 若前方行人阻擋了后方行人的視野和道路, 會(huì)影響和限制后方行人行走的舒適度和步行速度.此時(shí), 為了追求更好的可見視野和行走舒適度, 后方行人會(huì)進(jìn)行側(cè)向偏移, 形成斜列形式.

    以3 個(gè)毗鄰行人作為基本研究單元, 當(dāng)行人處于圖10 所示的視野最佳和步行舒適的拉鏈穩(wěn)定狀態(tài)時(shí), 在有效可見距離L 內(nèi), 在目標(biāo)行人正前方視野遮擋下, 行人可獲得的最佳可見關(guān)注視野為 θs,如圖11 所示.設(shè)定行人半徑為 rh, 頭部半徑為 rs,目標(biāo)行人相對(duì)正前方視線的最佳單側(cè)可見視野范圍標(biāo)準(zhǔn)值為 θh, 具體表達(dá)式為

    引入拉鏈敏感系數(shù) αi, 作為行人個(gè)體的屬性參數(shù), 用以描述行人客觀偏移的意愿程度.設(shè)定目標(biāo)行人 Pi相對(duì)周圍行人 Pj與自身正前方視線的可見范圍為 θij, 若 θij≥θh, 設(shè)定 θij=θh, 行人拉鏈敏感系數(shù)的公式為

    圖11 行人基本單元的拉鏈穩(wěn)定狀態(tài)(陰影區(qū)域?yàn)槟繕?biāo)行人正前方視野遮擋區(qū)域, 非陰影區(qū)域?yàn)槟繕?biāo)行人有效可見視野區(qū)域)Fig.11.Zipper stability state of pedestrian basic unit (the shadow area is the vision occlusion area and the non-shaded area is the effective visible area in front of the target pedestrian).

    當(dāng) αi越接近1 時(shí), 行人越有進(jìn)行偏移的客觀意愿;當(dāng) αi越 接 近0 時(shí), 行 人 越 沒有 進(jìn) 行 客 觀 偏 移 的意愿.

    根據(jù)研究數(shù)據(jù)[35]和日常觀察, 行人的投影形狀接近于橢圓, 如圖12 所示.在傳統(tǒng)的行人仿真中, 行人原橢圓投影區(qū)域通常以圓表示, 且它們的投影面積基本相等.本文將行人橢圓形的投影區(qū)域轉(zhuǎn)換為圓形區(qū)域, 根據(jù)人體尺寸相關(guān)數(shù)據(jù)[35], 可近似得到行人頭部半徑和行人半徑:

    其中, as是橢圓長軸長, bs是橢圓短軸長, rs是圓的半徑.根據(jù)人體尺寸相關(guān)數(shù)據(jù)[35], 近似得到行人頭部半徑 rs為0.08 m, 行人半徑 rh為0.3 m.

    圖12 行人投影區(qū)域Fig.12.Projection area of pedestrians.

    為了方便描述行人的偏移機(jī)制, 設(shè)定目標(biāo)行人Pi的原始位置坐標(biāo)為( Pi(x) , Pi(y) ), 橫向偏移后的位置坐標(biāo)為( Pi′(x) , Pi′(y) ); 周圍行人, 諸如 Pj,Pj2的位置坐標(biāo)分別為( Pj(x) , Pj(y) )和( Pj2(x) ,Pj2(y) ).設(shè)定目標(biāo)行人 Pi的可見視野范圍為 θi, 若目標(biāo)行人 Pi前方有兩個(gè)行人 Pj與 Pj2分別在兩側(cè)阻擋其視野, 則行人 Pi最大可見視野范圍 θi為θij與 θij2之和, 目標(biāo)行人 Pi與周圍行人 Pj的距離為 lij,

    首先, 確定行人拉鏈效應(yīng)的啟動(dòng)機(jī)制.設(shè)定當(dāng)目標(biāo)行人 Pi同時(shí)滿足與視野范圍內(nèi)距離自己最近行人Pj的距離 lij小于臨界值 L , 且行人拉鏈敏感系數(shù), 則目標(biāo)行人 Pi啟動(dòng)拉鏈效應(yīng), 進(jìn)行側(cè)向偏移 βi, 啟動(dòng)條件 γi為

    以圖13 為例, 目標(biāo)行人 P1的視野被左前方行人 P2阻擋(見圖13(a)), 且滿足拉鏈效應(yīng)啟動(dòng)條件γ1, 此時(shí)目標(biāo)行人 P1開始進(jìn)行側(cè)向偏移, 直到新位置處的視野 θ1′2達(dá)到穩(wěn)定視野 θh為止, 如圖13(b)所示.

    圖13 行人拉鏈效應(yīng)啟動(dòng)機(jī)制示意 (a)行人側(cè)向位移前;(b)行人側(cè)向位移后Fig.13.Starting mechanism of pedestrian zipper effect:(a) Before lateral displacement of pedestrians; (b) after lateral displacement of pedestrians.

    在圖13 中, 目標(biāo)行人只考慮了視野范圍內(nèi)距離自己最近的一個(gè)行人, 且符合拉鏈效應(yīng)啟動(dòng)條件, 然后進(jìn)行側(cè)向偏移的情況, 但是目標(biāo)行人也需要處理視野范圍內(nèi)同時(shí)面對(duì)兩側(cè)都有行人的情況, 因此需要進(jìn)一步研究.以圖14 為例, 目標(biāo)行人 P1為避開左前方行人 P2首先側(cè)向偏移到位置P1′(x,y) 處, 此時(shí)左側(cè)視野范圍為 θ1′2, 相對(duì)右前方行 人 P3的 右 側(cè) 視 野 范 圍 為 θ1′3.若 θ1′2=θh, 且θ1′3≥θh, 則 P1′(x,y) 為目標(biāo)行人 P1側(cè)向偏移后的穩(wěn)定位置, 如圖14(a)所示; 若 θ1′2=θh, 且θ1′3<θh, 則 P1′(x,y) 不是目標(biāo)行人 P1側(cè)向偏移后的穩(wěn)定位置, 行人 P1距離右前方行人 P3過近, 如圖14(b)所示; 為了平衡圖14(b)的情況, 目標(biāo)行人 P1在有限的空間內(nèi)稍微左移到更合適的位置, 此時(shí)θ1′2=θ1′3=θs/2 <θh, 如圖14(c)所示.

    定義5 基本單元活動(dòng)制造能力(MCBA)指單個(gè)基本制造活動(dòng)正常生產(chǎn)資料及約束條件運(yùn)作下所能達(dá)到的制造能力水平。

    參照?qǐng)D14(a)的情況, 若目標(biāo)行人 Pi相對(duì)前方一側(cè)行人 Pj最佳偏移后, 與另一側(cè)行人 Pj2的視野范圍 θij2≥θh, 則由(5)—(12)式可得目標(biāo)行人Pi側(cè)向偏移后的穩(wěn)定位置 Pi′(x,y) 及偏離大小 βi, 此時(shí)其與周圍行人 Pj, Pj2的距離分別為 li′j, li′j2; 參照?qǐng)D14(b)的情況, 若 θij2<θh, 則由(5)—(13)式可得目標(biāo)行人 Pi側(cè)向偏移后的均衡位置 Pi′(x,y) 以及偏離大小 βi.

    圖14 考慮兩側(cè)多個(gè)行人的拉鏈效應(yīng)偏移情況示意 (a)兩側(cè)均滿足最佳視野; (b)僅一側(cè)滿足最佳視野; (c)兩側(cè)妥協(xié)后穩(wěn)定視野Fig.14.Deviation of zipper effect when consider multiple pedestrians on both sides: (a) Both sides meet the optimum field of vision;(b) only one side satisfies the optimal field of vision; (c) stabilize visual field after compromise on both sides.

    其中, ni是目標(biāo)行人 Pi的正前方視線方 向, hij是目標(biāo)行人 Pi相對(duì)周圍行人 Pj的豎直距離.

    3 仿真模型

    3.1 速度修正

    行人在正常行走時(shí), 以接近勻速的狀態(tài)移動(dòng).社會(huì)力影響行人行走的本質(zhì)是影響行人的速度變化, 進(jìn)而影響行人的位移變化; 而傳統(tǒng)的社會(huì)力模型更適合仿真擁擠狀態(tài)下的行人相互作用, 在仿真正常狀態(tài)下的行人移動(dòng)時(shí), 會(huì)出現(xiàn)不合理的行人現(xiàn)象[16,27,28,36].為了克服傳統(tǒng)社會(huì)力模型在仿真正常狀態(tài)下行人流的一些不足, 簡化模型參數(shù)校定的復(fù)雜性和不準(zhǔn)確性, 克服行人仿真重疊性和跳動(dòng)性的問題, 提出了一種綜合考慮行人視野范圍內(nèi)周圍不同方向和距離行人或障礙物對(duì)目標(biāo)行人速度影響的速度修正模型, 適合仿真正常行走狀態(tài)下的行人流.

    通常在目標(biāo)行人視野前方一定范圍內(nèi)的周圍行人會(huì)對(duì)目標(biāo)行人的行走決策產(chǎn)生影響, 而目標(biāo)行人視野后方的周圍行人對(duì)目標(biāo)行人的影響較小, 且視野前方不同方向和距離的周圍行人對(duì)目標(biāo)行人行走速度的影響各不相同.結(jié)合Voronoi 圖的幾何特性, 可以有效得到目標(biāo)行人 P0周圍的行人 P1, P2,P3, P4, P5, 如圖15 所示.設(shè)定在目標(biāo)行人 P0正前方視野范圍(180°)內(nèi), 周圍行人可以對(duì)目標(biāo)行人P0速度產(chǎn)生影響的范圍為半徑 r 的半圓.在0 ≤r ≤r1范圍內(nèi), 周圍行人對(duì)目標(biāo)行人 P0的速度影響較大, 如行人 P1; 在 r1<r ≤r2范圍內(nèi), 周圍行人對(duì)目標(biāo)行人 P0的速度影響較小, 如行人 P2; 在r >r2范圍內(nèi), 周圍行人不對(duì)目標(biāo)行人 P0的速度產(chǎn)生影響, 如行人 P3; 同樣, 在目標(biāo)行人 P0視野后方,周圍行人不對(duì)目標(biāo)行人的速度產(chǎn)生影響, 如行人P4和 P5.

    圖15 周圍行人對(duì)目標(biāo)行人的速度影響Fig.15.Velocity effects of surrounding pedestrians on the target pedestrian.

    周圍行人對(duì)目標(biāo)行人速度的影響除了與目標(biāo)行人的相對(duì)距離有關(guān), 也和周圍行人相對(duì)目標(biāo)行人移動(dòng)方向的角度有關(guān).如, 在 0 ≤r ≤r1的范圍內(nèi),正前方行人 P1對(duì)目標(biāo)行人 P0的視野遮擋影響最大,且目標(biāo)行人 P0對(duì)正前方的視野最為關(guān)注, 則行人P1會(huì)對(duì)目標(biāo)行人 P0產(chǎn)生相對(duì)其他行人更大的速度影響, 以使目標(biāo)行人 P0減速甚至靜止.

    基于速度修正描述影響行人移動(dòng)的心理趨勢(shì)和物理行為, 本文提出一種速度修正模型.設(shè)定行人在正常移動(dòng)時(shí)主要受自驅(qū)動(dòng)速度、周圍行人和障礙物排斥速度的影響.自驅(qū)動(dòng)速度提供行人前進(jìn)的速度, 表示行人按照期望速度前進(jìn)的愿望; 排斥速度表示行人受周圍行人和障礙物的排斥心理.該模型將行人的移動(dòng)看成每單位時(shí)間步的勻速運(yùn)動(dòng), 計(jì)算行人在不同速度影響下的位置更新過程為

    v0g((dij-rij)θij)nij表示目標(biāo)行人受周圍行人的排斥心理.在目標(biāo)行人視野范圍內(nèi), 若dijrij≤d1, 當(dāng) θij=0 時(shí), g (x) 函數(shù)取值為 k1, 當(dāng)θij/=0 時(shí), g (x) 函 數(shù) 取 值 為 k2; 若 d1<dij-rij≤d2,g(x) 函 數(shù) 取 值 為 k3; 若 dij-rij>d2, g (x) 函 數(shù) 取值為 k4.其中 dij表示兩行人質(zhì)心間的距離, rij表示兩行人半徑之和, θij表示兩行人位置與目標(biāo)行人速度方向的夾角, nij為周圍行人質(zhì)心指向目標(biāo)行人質(zhì)心的單位向量.

    v0g(diw-ri)niw表示目標(biāo)行人受來自墻壁W的速度影響, 與來自周圍行人的速度影響類似, 在目標(biāo)行人視野范圍(180°)內(nèi), 若 diw-ri≤d3,g(x) 函 數(shù) 取 值 為 k5; 若 diw-ri>d3, g (x) 函 數(shù) 取值為 k6.其中, diw為目標(biāo)行人與墻壁 W 的垂直距離, ri為目標(biāo)行人的半徑, niw為與墻壁垂直且指向目標(biāo)行人的單位向量.通過仿真實(shí)驗(yàn)測(cè)定, 速度修正模型中各參數(shù)的取值如表1 所列.

    表1 速度修正模型參數(shù)取值Table 1.Parameter values of velocity correction model.

    3.2 偏移規(guī)則

    為了真實(shí)再現(xiàn)人群的拉鏈現(xiàn)象, 考慮行人視野最佳和步行舒適的需求, 結(jié)合速度修正模型, 制定行人拉鏈效應(yīng)的側(cè)向偏移規(guī)則, 如圖16 所示.

    首先, 考慮行人是否有進(jìn)行拉鏈效應(yīng)的主觀偏移意愿, 因?yàn)楦鞣N因素影響, 行人有時(shí)會(huì)做出非理性的行為.若行人主觀意愿 κ =1 , 行人有進(jìn)行主觀偏移的意愿; 否則, 行人不會(huì)刻意做出偏移行為.

    其次, 若目標(biāo)行人有進(jìn)行主觀偏移的意愿, 且符合行人進(jìn)行客觀偏移的條件, 即符合(7)式的拉鏈效應(yīng)啟動(dòng)條件 γ1, 則根據(jù)(5)—(13)式可得單個(gè)行人 Pi拉鏈效應(yīng)偏移后的均衡位置 Pi′(x,y).

    圖16 行人拉鏈效應(yīng)偏移規(guī)則流程圖Fig.16.Flow diagram of pedestrian zipper effect deviation rules.

    然后, 根據(jù)目標(biāo)行人 Pi當(dāng)前位置 Pi(x,y) 與其預(yù)測(cè)后的偏移位置 Pi′(x,y) 之間的相對(duì)位置關(guān)系,確定目標(biāo)行人的偏移方向 ni.賦予目標(biāo)行人 Pi一個(gè)主動(dòng)偏移速度 vi, 直到目標(biāo)行人達(dá)到偏移后的均衡位置, (14)式進(jìn)行改進(jìn)為

    在確定行人 Pi的主動(dòng)偏移速度 vi時(shí), 需要考慮模型的仿真精度與效率, 考慮到行人的半徑設(shè)定為0.3 m, 在0.5 s 的反應(yīng)時(shí)間內(nèi), 若取行人 Pi的主動(dòng)偏移速度 vi為0.2 m/s, 則行人在理想狀態(tài)下側(cè)向偏移0.1 m, 比較符合行人的正常行為反應(yīng); 同時(shí),通過仿真實(shí)驗(yàn)測(cè)定, 可以獲得較好的仿真效果, 因此在本文中, 行人 Pi的主動(dòng)偏移速度 vi取0.2 m/s.

    以圖17 為例, 示意有主動(dòng)偏移意愿的目標(biāo)行人 P1產(chǎn)生主動(dòng)偏移速度 v1, 直到其到達(dá)偏移后的均衡位置 P1′(x,y).

    圖17 有偏移速度的行人側(cè)向偏移示意Fig.17.Sideways deviation of pedestrians with deviation velocity.

    4 仿真分析

    為了實(shí)時(shí)區(qū)分每個(gè)行人的局部密度, 根據(jù)行人局部密度的大小, 對(duì)Voronoi 元胞填充不同的顏色, 綠色表示低密度, 黃色表示中密度, 橙色表示中高密度, 紅色表示高密度.為了簡約, 本文均采用一致性的行人進(jìn)行仿真.行人仿真時(shí), 根據(jù)Voronoi 元胞的顏色, 可以觀察到不同位置和時(shí)間行人的擁擠程度, 且Voronoi 元胞的形狀和顏色會(huì)隨著行人位置和密度的變化而實(shí)時(shí)變化.

    4.1 通道循環(huán)流

    在10 m × 3.5 m 的直通道走廊場景中進(jìn)行行人單向循環(huán)流仿真實(shí)驗(yàn), 驗(yàn)證模型的可靠性.走廊下端是行人的產(chǎn)生區(qū)域, 上端是出口, 在走廊起始端2 m × 3.5 m 的范圍內(nèi)隨機(jī)生成20 個(gè)行人, 行人向出口移動(dòng).在實(shí)驗(yàn)中, 行人采用的是周期性循環(huán)流, 行人從走廊下端進(jìn)入并從上端離開后, 從下端重新進(jìn)入走廊; 走廊的左右邊界采用的是封閉性實(shí)體邊界, 上下兩端采用的是開放性截?cái)噙吔?采用開放性截?cái)噙吔鐣r(shí), 在走廊上下端邊緣處行人的Voronoi 元胞, 首先基于開放性邊界生成被截割的Voronoi 元胞, 然后利用上下端的邊界將其截割, 從而生成適應(yīng)邊界的新Voronoi 邊和點(diǎn)[27].

    設(shè)定所有行人都有進(jìn)行偏移的主觀意愿 κ =1 ,行人在走廊移動(dòng)的過程如圖18 所示.在仿真過程中, 行人數(shù)量和整體密度是固定的.仿真開始時(shí),行人隨機(jī)分布在走廊中, 且擁有各自的局部密度;當(dāng)行人從初始生成階段過渡到穩(wěn)定階段, 行人為追求自身的舒適距離和視野, 會(huì)不斷調(diào)整自身的位置, 呈現(xiàn)出一定的拉鏈現(xiàn)象; 最后行人基本均勻分布, 由Voronoi 元胞顏色可見行人的局部密度基本一致, 行人之間分散形成四層隊(duì)列, 且每個(gè)行人的速度都大致相等, 行人勻速向前移動(dòng).同時(shí), 由于在行人Voronoi 圖的計(jì)算中, 仿真場景的上下邊界采用開放性截?cái)噙吔纾?因此, 當(dāng)行人從走廊上端離開并從下端重新進(jìn)入走廊時(shí), 因距離走廊底端邊界較近, 會(huì)增加行人的局部密度.當(dāng)走廊整體行人密度較小時(shí), 重新進(jìn)入走廊的行人與前方行人距離較遠(yuǎn), 雙方行人相互之間的速度影響較小或沒有影響;當(dāng)走廊整體行人密度較大, 重新進(jìn)入走廊的行人與前方行人距離較近, 雙方行人之間會(huì)產(chǎn)生一定的速度影響, 但經(jīng)過短暫時(shí)間行進(jìn)后, 相互影響的行人會(huì)自發(fā)調(diào)節(jié)重新進(jìn)入穩(wěn)定狀態(tài).該仿真真實(shí)再現(xiàn)了常態(tài)行人流, 可以發(fā)現(xiàn)行人維持自身領(lǐng)域, 克服了傳統(tǒng)社會(huì)力模型抖動(dòng)性和重疊性的不足.

    圖18 直通道行人循環(huán)流仿真過程截圖( w =3.5 m ,l =10 m , κ =1 )Fig.18.Screenshot of simulation process of pedestrian flow in the straight corridor ( w =3.5 m , l =10 m , κ =1 ).

    以行人#4 為例, 可以看到在不同仿真時(shí)間步,行人#4 進(jìn)行了明顯的拉鏈效應(yīng)偏移, 以獲得更好的視野和步行舒適度, 如圖19 所示.

    4.2 對(duì)比分析

    在直通道中隨機(jī)產(chǎn)生20 個(gè)行人, 隨著通道寬度 B 的增加, 行人自主選擇的空間逐漸變大, 行人可以充分地進(jìn)行拉鏈效應(yīng)偏移.當(dāng)行人形成穩(wěn)定的行人流時(shí), 可以觀察到不同寬度通道中形成不同的行人隊(duì)列層數(shù), 通道寬度越大, 行人隊(duì)列層數(shù)越多,與行人實(shí)驗(yàn)[33]的結(jié)論相吻合, 如圖20 所示.

    為了獲得行人速度與密度的關(guān)系, 對(duì)不同密度條件下(人/m2)的行人進(jìn)行單向循環(huán)流仿真, 并通過基本圖驗(yàn)證模型的可靠性.在每個(gè)行人密度點(diǎn)的仿真實(shí)驗(yàn)中, 最開始行人隨機(jī)地分布于走廊中, 行人的目標(biāo)點(diǎn)位于走廊上端, 每次仿真實(shí)驗(yàn)重復(fù)10 次,每次持續(xù)90 個(gè)仿真時(shí)間步, 實(shí)驗(yàn)的前30 個(gè)仿真時(shí)間步當(dāng)做熱身時(shí)間, 后60 個(gè)仿真時(shí)間步的數(shù)據(jù)用于統(tǒng)計(jì)分析.圖21 將速度修正模型與其他實(shí)證研究的數(shù)據(jù)進(jìn)行了對(duì)比, 可以發(fā)現(xiàn)速度修正模型的速度-密度關(guān)系與實(shí)證數(shù)據(jù)符合得較好[37-41].

    通過改變行人拉鏈效應(yīng)主觀偏移意愿與無偏移意愿的占比 u , 可以得到不同密度下, 行人主觀偏移意愿占比與行人平均速度的關(guān)系, 如圖22 所示; 也可以獲得不同密度下, 行人主觀偏移意愿占比與行人平均局部密度的關(guān)系, 如圖23 所示.可以發(fā)現(xiàn), 與不考慮拉鏈效應(yīng)相比, 傾向主動(dòng)進(jìn)行拉鏈效應(yīng)偏移的行人占比越大, 越有助于提高通道內(nèi)行人的移動(dòng)速度、舒適度和空間利用率, 與行人實(shí)驗(yàn)[30]的研究結(jié)果相符合.

    圖19 行人#4 在不同仿真時(shí)間步的拉鏈效應(yīng)偏移截圖Fig.19.Screenshot of zipper effect of pedestrian #4 at different simulation time steps.

    圖20 行人拉鏈層數(shù)與通道寬度關(guān)系Fig.20.Relationship between the number of pedestrian zipper layers and the width of the corridor.

    圖21 直通道單向行人流基本圖[37-41](該模型仿真數(shù)據(jù)與實(shí)證數(shù)據(jù)(Older[37], Mori[38], Zhang[39], Weidmann[40], Hankin[41])進(jìn)行對(duì)比)Fig.21.Fundamental diagram of unidirectional pedestrian flow in the corridor[37-41].Actual data (Older[37], Mori[38],Zhang[39], Weidmann[40], Hankin[41]) are gathered to compare with our model.

    圖22 不同密度下行人主觀偏移意愿占比與行人平均速度關(guān)系Fig.22.Relationship between the proportion of pedestrian subjective deviation intention and pedestrian average velocity under different overall densities.

    同時(shí), 為驗(yàn)證Voronoi 圖劃分對(duì)模型仿真效率的影響, 僅分析了針對(duì)Voronoi 圖劃分所需要的時(shí)間.在Win7 64 位操作系統(tǒng), 酷睿雙核i7-6700 CPU 3.4 GHz, 16G 內(nèi)存的臺(tái)式電腦上, 基于Matlab R2018b 的編程環(huán)境, 通過仿真計(jì)算了Voronoi 圖劃分所需的時(shí)間.對(duì)20 人的Voronoi 圖劃分, 需要0.499254 s; 對(duì)200 人的Voronoi 圖劃分, 需要0.992942 s; 對(duì)1000 人的Voronoi 圖劃分, 需要4.952742 s; 對(duì)10000 人的Voronoi 圖劃分, 需要755.423851 s.通過實(shí)際的計(jì)算時(shí)間可以看出: 當(dāng)行人數(shù)量較少時(shí), Voronoi 圖劃分耗時(shí)在可承受的范圍之內(nèi); 當(dāng)行人數(shù)量較大時(shí), Voronoi 圖劃分耗時(shí)急劇增長.因此, 為了減少Voronoi 圖劃分耗時(shí)對(duì)模型仿真效率的影響, 為了提高基于Voronoi 圖模型整體的仿真效率, 使其適用于大規(guī)模行人的仿真, 需要進(jìn)一步研究并行環(huán)境下模型的仿真計(jì)算,包括基于空間分割、基于仿真規(guī)則、基于進(jìn)程的并行計(jì)算等.

    圖23 不同密度下行人主觀偏移意愿占比與行人平均局部密度關(guān)系Fig.23.Relationship between the proportion of pedestrian subjective deviation intention and pedestrian average local density under different overall densities.

    5 結(jié) 論

    本文提出行人拉鏈現(xiàn)象的生成機(jī)理, 并構(gòu)建基于Voronoi 圖的速度修正模型進(jìn)行直通道單向循環(huán)行人流仿真研究.從視野最佳和步行舒適的角度分析拉鏈現(xiàn)象的生成機(jī)理, 引入拉鏈敏感系數(shù)描述行人客觀偏移的意愿程度, 提出行人的偏移機(jī)制.引入Voronoi 圖, 有效確定目標(biāo)行人視野范圍內(nèi)的周圍行人; 基于Voronoi 圖考慮不同距離和方向周圍行人對(duì)目標(biāo)行人移動(dòng)速度的影響, 構(gòu)建速度修正模型, 考慮行人是否有偏移傾向的主觀意愿, 并嵌入偏移規(guī)則, 模擬再現(xiàn)行人的拉鏈現(xiàn)象.仿真結(jié)果真實(shí)再現(xiàn)了直通道常態(tài)行人流, 克服了傳統(tǒng)社會(huì)力模型抖動(dòng)性和重疊性的不足, 可以觀察到均勻分布的自組織行人流和明顯的行人拉鏈效應(yīng), 行人的拉鏈層數(shù)與通道寬度成正比.通過提取行人仿真實(shí)驗(yàn)數(shù)據(jù)與實(shí)證數(shù)據(jù)對(duì)比, 該模型速度密度關(guān)系的基本圖與實(shí)證數(shù)據(jù)符合得較好; 與不考慮拉鏈效應(yīng)相比, 傾向主動(dòng)進(jìn)行側(cè)向偏移的行人占比越大, 越有助于提高通道內(nèi)行人的移動(dòng)速度、舒適度和空間利用率.

    猜你喜歡
    影響模型
    一半模型
    是什么影響了滑動(dòng)摩擦力的大小
    哪些顧慮影響擔(dān)當(dāng)?
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    沒錯(cuò),痛經(jīng)有時(shí)也會(huì)影響懷孕
    媽媽寶寶(2017年3期)2017-02-21 01:22:28
    3D打印中的模型分割與打包
    擴(kuò)鏈劑聯(lián)用對(duì)PETG擴(kuò)鏈反應(yīng)與流變性能的影響
    中國塑料(2016年3期)2016-06-15 20:30:00
    基于Simulink的跟蹤干擾對(duì)跳頻通信的影響
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    日本与韩国留学比较| 91精品国产九色| 国产 一区 欧美 日韩| 一个人看视频在线观看www免费| 国产精品女同一区二区软件| 在现免费观看毛片| 干丝袜人妻中文字幕| 国产精品秋霞免费鲁丝片| 日韩 亚洲 欧美在线| 国语对白做爰xxxⅹ性视频网站| 91久久精品国产一区二区成人| 交换朋友夫妻互换小说| 男女免费视频国产| 亚洲av成人精品一区久久| 91精品国产九色| 日韩免费高清中文字幕av| 一级毛片久久久久久久久女| 91精品伊人久久大香线蕉| 在线观看免费高清a一片| 精品久久久精品久久久| 国产人妻一区二区三区在| 日本黄色片子视频| 亚洲成人一二三区av| 高清毛片免费看| 中文字幕久久专区| 亚洲在久久综合| 熟女电影av网| 日产精品乱码卡一卡2卡三| 国内揄拍国产精品人妻在线| 国产精品三级大全| 大片免费播放器 马上看| 午夜精品国产一区二区电影| av天堂中文字幕网| 亚州av有码| 亚洲精品国产av成人精品| 久久精品国产a三级三级三级| 少妇高潮的动态图| 亚洲欧洲国产日韩| 国产免费视频播放在线视频| 欧美成人午夜免费资源| 国产有黄有色有爽视频| 国产国拍精品亚洲av在线观看| 中文字幕精品免费在线观看视频 | 免费播放大片免费观看视频在线观看| 国产免费一区二区三区四区乱码| 直男gayav资源| 亚洲av二区三区四区| 26uuu在线亚洲综合色| 日本av手机在线免费观看| 干丝袜人妻中文字幕| 久久久久久久久久成人| 亚洲av成人精品一二三区| 少妇猛男粗大的猛烈进出视频| 中文字幕亚洲精品专区| 人体艺术视频欧美日本| 欧美bdsm另类| 一级二级三级毛片免费看| 亚洲国产精品成人久久小说| 18禁动态无遮挡网站| 亚洲国产精品999| 久久99精品国语久久久| 少妇猛男粗大的猛烈进出视频| 国产91av在线免费观看| 熟女电影av网| 2018国产大陆天天弄谢| 人妻系列 视频| 亚洲精品一二三| 这个男人来自地球电影免费观看 | 久久久国产一区二区| 国产老妇伦熟女老妇高清| 欧美97在线视频| 1000部很黄的大片| 中国国产av一级| 少妇熟女欧美另类| 人人妻人人爽人人添夜夜欢视频 | 一本—道久久a久久精品蜜桃钙片| 最后的刺客免费高清国语| 伊人久久国产一区二区| 国产爱豆传媒在线观看| 午夜免费鲁丝| 国产 一区 欧美 日韩| 成年美女黄网站色视频大全免费 | 成人国产av品久久久| 色视频在线一区二区三区| 婷婷色av中文字幕| 老司机影院毛片| 亚洲av.av天堂| 中文字幕精品免费在线观看视频 | 91在线精品国自产拍蜜月| 欧美成人精品欧美一级黄| 黄色一级大片看看| 亚洲自偷自拍三级| 美女高潮的动态| 内射极品少妇av片p| av国产免费在线观看| 91精品伊人久久大香线蕉| 国产日韩欧美亚洲二区| 日韩一本色道免费dvd| 国产高清国产精品国产三级 | 日日撸夜夜添| 狂野欧美白嫩少妇大欣赏| 亚洲欧美成人精品一区二区| 日韩不卡一区二区三区视频在线| 精品久久久精品久久久| 欧美高清成人免费视频www| 久久99蜜桃精品久久| 啦啦啦中文免费视频观看日本| 蜜桃久久精品国产亚洲av| 大陆偷拍与自拍| 中国三级夫妇交换| 国产乱来视频区| 99re6热这里在线精品视频| 26uuu在线亚洲综合色| 97热精品久久久久久| tube8黄色片| 人妻夜夜爽99麻豆av| 久久久精品94久久精品| 亚洲国产成人一精品久久久| 永久免费av网站大全| av一本久久久久| 国产一级毛片在线| 欧美日韩综合久久久久久| 久久久久视频综合| 欧美精品亚洲一区二区| 亚洲av不卡在线观看| 国产成人精品福利久久| 亚洲欧洲日产国产| 亚洲人成网站在线播| 国产av精品麻豆| 国产成人精品婷婷| 中文字幕亚洲精品专区| 国产爱豆传媒在线观看| 午夜福利视频精品| 人妻 亚洲 视频| 在线观看免费日韩欧美大片 | 免费久久久久久久精品成人欧美视频 | 51国产日韩欧美| 国产成人一区二区在线| 亚洲一级一片aⅴ在线观看| 国产人妻一区二区三区在| 久久久欧美国产精品| 五月伊人婷婷丁香| 精品一区二区三区视频在线| 看免费成人av毛片| 日韩大片免费观看网站| 免费看av在线观看网站| 你懂的网址亚洲精品在线观看| 免费久久久久久久精品成人欧美视频 | 人妻一区二区av| 日韩不卡一区二区三区视频在线| 久久精品夜色国产| 精品国产一区二区三区久久久樱花 | 亚洲四区av| 欧美变态另类bdsm刘玥| 七月丁香在线播放| 超碰97精品在线观看| 人人妻人人澡人人爽人人夜夜| 在线观看三级黄色| 狠狠精品人妻久久久久久综合| 色网站视频免费| 国产国拍精品亚洲av在线观看| 国产亚洲一区二区精品| 美女脱内裤让男人舔精品视频| 麻豆国产97在线/欧美| 妹子高潮喷水视频| 国产免费一级a男人的天堂| 一级a做视频免费观看| 婷婷色av中文字幕| 22中文网久久字幕| 六月丁香七月| 亚洲欧美成人精品一区二区| 精品午夜福利在线看| 国产视频内射| 在线看a的网站| 久久人人爽av亚洲精品天堂 | 97在线视频观看| 欧美丝袜亚洲另类| 人妻夜夜爽99麻豆av| 欧美性感艳星| 亚洲国产精品国产精品| 自拍偷自拍亚洲精品老妇| 亚洲婷婷狠狠爱综合网| 国产午夜精品一二区理论片| 亚洲av.av天堂| 午夜激情福利司机影院| 人体艺术视频欧美日本| 国产黄片视频在线免费观看| 成人一区二区视频在线观看| 精品久久久久久久久亚洲| 免费观看的影片在线观看| 日韩在线高清观看一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 亚洲欧美中文字幕日韩二区| 免费av不卡在线播放| 多毛熟女@视频| 国产91av在线免费观看| 国产精品福利在线免费观看| 精品一区二区三区视频在线| 成人亚洲欧美一区二区av| 多毛熟女@视频| 欧美成人午夜免费资源| 日韩一本色道免费dvd| 天天躁日日操中文字幕| 国产伦理片在线播放av一区| 国产精品免费大片| 涩涩av久久男人的天堂| 久久精品国产自在天天线| 亚洲av免费高清在线观看| 久久av网站| 18禁裸乳无遮挡动漫免费视频| 性色av一级| 高清午夜精品一区二区三区| 久久国产亚洲av麻豆专区| 久热这里只有精品99| 国产一区二区在线观看日韩| 建设人人有责人人尽责人人享有的 | www.色视频.com| 一级黄片播放器| 美女脱内裤让男人舔精品视频| 我要看日韩黄色一级片| 久久国产精品男人的天堂亚洲 | 欧美97在线视频| 日韩中字成人| 久久99热6这里只有精品| 肉色欧美久久久久久久蜜桃| 免费大片18禁| 人妻 亚洲 视频| 国产成人91sexporn| 韩国av在线不卡| 精品国产露脸久久av麻豆| 久久影院123| 又黄又爽又刺激的免费视频.| 日本与韩国留学比较| 99久国产av精品国产电影| 日韩,欧美,国产一区二区三区| 少妇的逼好多水| 亚洲国产成人一精品久久久| 国产精品成人在线| 乱系列少妇在线播放| 看十八女毛片水多多多| 超碰av人人做人人爽久久| 国产成人精品一,二区| 免费观看a级毛片全部| 中文字幕精品免费在线观看视频 | 五月天丁香电影| 女人久久www免费人成看片| 免费av中文字幕在线| 亚洲av不卡在线观看| 男女边摸边吃奶| 精品久久久精品久久久| 亚洲av二区三区四区| 亚洲av电影在线观看一区二区三区| 精品亚洲成a人片在线观看 | a级一级毛片免费在线观看| 99热网站在线观看| 狠狠精品人妻久久久久久综合| 亚洲一区二区三区欧美精品| 男人和女人高潮做爰伦理| 看免费成人av毛片| 少妇裸体淫交视频免费看高清| 成人综合一区亚洲| 十分钟在线观看高清视频www | 又黄又爽又刺激的免费视频.| 美女主播在线视频| 欧美丝袜亚洲另类| 午夜福利在线在线| 51国产日韩欧美| 免费观看av网站的网址| av国产久精品久网站免费入址| 国产伦理片在线播放av一区| 男女国产视频网站| 一级爰片在线观看| 99热网站在线观看| 亚洲精品一区蜜桃| 久久精品国产自在天天线| 人妻系列 视频| 成人毛片60女人毛片免费| 成人无遮挡网站| 亚洲在久久综合| a级毛片免费高清观看在线播放| 精品人妻一区二区三区麻豆| 少妇人妻 视频| 性色av一级| 国产爱豆传媒在线观看| 麻豆成人午夜福利视频| av在线蜜桃| 午夜老司机福利剧场| 十八禁网站网址无遮挡 | 一级黄片播放器| 男女免费视频国产| 一本一本综合久久| 成人高潮视频无遮挡免费网站| a级毛片免费高清观看在线播放| 亚洲人成网站高清观看| 色吧在线观看| 乱码一卡2卡4卡精品| 国产一区二区在线观看日韩| 国产v大片淫在线免费观看| 亚洲精品乱码久久久久久按摩| 国产亚洲精品久久久com| 亚洲av二区三区四区| 久久青草综合色| 久久久欧美国产精品| 国产国拍精品亚洲av在线观看| 老师上课跳d突然被开到最大视频| 精品久久久久久久久亚洲| 国产亚洲一区二区精品| 男女啪啪激烈高潮av片| 精品人妻熟女av久视频| 国产女主播在线喷水免费视频网站| 能在线免费看毛片的网站| 国产欧美亚洲国产| 欧美zozozo另类| 人妻一区二区av| 成人影院久久| 97超碰精品成人国产| 欧美国产精品一级二级三级 | 精品人妻偷拍中文字幕| 国产精品久久久久久精品电影小说 | 午夜免费鲁丝| 一个人看视频在线观看www免费| 在线观看国产h片| 免费人成在线观看视频色| 中国三级夫妇交换| 国产精品麻豆人妻色哟哟久久| 亚洲欧洲国产日韩| 亚洲美女搞黄在线观看| 久久久成人免费电影| 中文字幕制服av| 国产真实伦视频高清在线观看| 久久精品夜色国产| 亚洲精品久久久久久婷婷小说| 国产美女午夜福利| 国产乱人视频| 久久人人爽人人片av| 国产乱人视频| 久久久久性生活片| 中文欧美无线码| 亚洲欧美日韩另类电影网站 | 亚洲欧美清纯卡通| av.在线天堂| 国产在视频线精品| 你懂的网址亚洲精品在线观看| 丰满少妇做爰视频| 国产一区二区三区综合在线观看 | 国产精品人妻久久久影院| 一本久久精品| a级一级毛片免费在线观看| 少妇丰满av| 老司机影院成人| 看免费成人av毛片| 一级片'在线观看视频| 日本黄色片子视频| 国产在线免费精品| av免费在线看不卡| a 毛片基地| 黄色配什么色好看| 亚洲av福利一区| 久久精品久久久久久久性| 久久国产精品大桥未久av | 天堂8中文在线网| av在线老鸭窝| 亚洲欧美日韩无卡精品| 亚洲四区av| 日日摸夜夜添夜夜添av毛片| 99热网站在线观看| a级毛色黄片| 国产熟女欧美一区二区| 少妇裸体淫交视频免费看高清| 国产伦理片在线播放av一区| 日本爱情动作片www.在线观看| 一个人看视频在线观看www免费| 黄色怎么调成土黄色| 国产成人freesex在线| 久久久久久久精品精品| 国产白丝娇喘喷水9色精品| 男男h啪啪无遮挡| 亚洲欧美日韩东京热| 女人十人毛片免费观看3o分钟| 精品酒店卫生间| 美女福利国产在线 | 超碰av人人做人人爽久久| 国产精品99久久久久久久久| 啦啦啦中文免费视频观看日本| 国产欧美亚洲国产| 欧美xxxx黑人xx丫x性爽| 亚洲成人av在线免费| 国产在视频线精品| 国产高潮美女av| 麻豆成人午夜福利视频| 久久久午夜欧美精品| 久久人人爽人人片av| 日本与韩国留学比较| 亚洲av综合色区一区| 少妇人妻一区二区三区视频| 在线亚洲精品国产二区图片欧美 | 精品人妻视频免费看| 国产亚洲一区二区精品| 少妇丰满av| 久久久久久久久久久丰满| 国产一区二区在线观看日韩| 久久精品国产亚洲网站| 日本免费在线观看一区| 免费大片18禁| 亚洲第一av免费看| 成人黄色视频免费在线看| 国产精品熟女久久久久浪| 99国产精品免费福利视频| 97超视频在线观看视频| 最近的中文字幕免费完整| www.av在线官网国产| 高清av免费在线| 成年av动漫网址| 91午夜精品亚洲一区二区三区| 97热精品久久久久久| kizo精华| 亚洲精品色激情综合| 久久久久久久大尺度免费视频| 欧美精品人与动牲交sv欧美| 一级毛片我不卡| 国产成人freesex在线| 成人特级av手机在线观看| 丝袜喷水一区| 少妇丰满av| 性色av一级| 亚洲精品自拍成人| 国产精品久久久久久久久免| 亚洲综合色惰| 久久这里有精品视频免费| 久久影院123| 人妻制服诱惑在线中文字幕| av国产久精品久网站免费入址| 永久网站在线| 人妻系列 视频| 婷婷色综合大香蕉| 观看免费一级毛片| 激情五月婷婷亚洲| 人妻制服诱惑在线中文字幕| 一本久久精品| 日韩av不卡免费在线播放| 欧美丝袜亚洲另类| 日韩不卡一区二区三区视频在线| 成人美女网站在线观看视频| a级毛片免费高清观看在线播放| 在线精品无人区一区二区三 | 日本-黄色视频高清免费观看| 国产精品麻豆人妻色哟哟久久| 大片电影免费在线观看免费| 熟女av电影| 久久精品久久精品一区二区三区| .国产精品久久| 国内精品宾馆在线| 午夜福利在线在线| 内地一区二区视频在线| 国产男人的电影天堂91| kizo精华| 午夜福利视频精品| 两个人的视频大全免费| 国产无遮挡羞羞视频在线观看| 日本wwww免费看| 黄色欧美视频在线观看| 日本av免费视频播放| 99热这里只有是精品在线观看| 亚洲精品乱码久久久久久按摩| 国产精品久久久久久精品电影小说 | 哪个播放器可以免费观看大片| 熟女人妻精品中文字幕| 国产熟女欧美一区二区| 欧美xxxx黑人xx丫x性爽| 不卡视频在线观看欧美| 亚洲精品国产色婷婷电影| 亚洲丝袜综合中文字幕| 一本—道久久a久久精品蜜桃钙片| 亚洲av在线观看美女高潮| 精品国产三级普通话版| 蜜桃亚洲精品一区二区三区| 亚洲精品久久久久久婷婷小说| 免费黄色在线免费观看| 两个人的视频大全免费| 亚洲av免费高清在线观看| 在线天堂最新版资源| 最后的刺客免费高清国语| 国产精品一区二区在线观看99| 国产成人a∨麻豆精品| 精品酒店卫生间| 国产高清有码在线观看视频| 久久精品国产亚洲av涩爱| 丰满迷人的少妇在线观看| 国产欧美亚洲国产| 久久久久网色| 观看av在线不卡| 久久久久久久亚洲中文字幕| 欧美性感艳星| 国产av码专区亚洲av| 亚洲欧美一区二区三区国产| 看非洲黑人一级黄片| 亚洲第一av免费看| 色视频在线一区二区三区| 国产熟女欧美一区二区| 麻豆精品久久久久久蜜桃| 日韩制服骚丝袜av| 小蜜桃在线观看免费完整版高清| 美女福利国产在线 | 最后的刺客免费高清国语| 男女啪啪激烈高潮av片| 欧美高清性xxxxhd video| 九色成人免费人妻av| 国产精品99久久久久久久久| 高清欧美精品videossex| 欧美亚洲 丝袜 人妻 在线| 久久精品国产a三级三级三级| 多毛熟女@视频| 色综合色国产| 精品人妻熟女av久视频| 视频区图区小说| h视频一区二区三区| 一级毛片久久久久久久久女| 亚洲国产高清在线一区二区三| 免费在线观看成人毛片| 九九久久精品国产亚洲av麻豆| 久久久精品免费免费高清| 久久国产亚洲av麻豆专区| 交换朋友夫妻互换小说| 久久久久久久国产电影| 亚洲国产欧美人成| 女人久久www免费人成看片| 丝袜脚勾引网站| 最黄视频免费看| 国产淫语在线视频| 91久久精品国产一区二区成人| 交换朋友夫妻互换小说| 精品99又大又爽又粗少妇毛片| 一级片'在线观看视频| 久久鲁丝午夜福利片| 亚洲av免费高清在线观看| 亚洲熟女精品中文字幕| 亚洲三级黄色毛片| 亚洲aⅴ乱码一区二区在线播放| 高清黄色对白视频在线免费看 | 九草在线视频观看| 国产人妻一区二区三区在| 免费播放大片免费观看视频在线观看| a级一级毛片免费在线观看| 国产精品久久久久久av不卡| 在线观看av片永久免费下载| av线在线观看网站| 国产精品一区二区在线不卡| 在线天堂最新版资源| 欧美3d第一页| 国产亚洲91精品色在线| 一个人看视频在线观看www免费| 男人和女人高潮做爰伦理| 婷婷色综合大香蕉| 久热这里只有精品99| 亚洲高清免费不卡视频| 午夜日本视频在线| 涩涩av久久男人的天堂| 国产综合精华液| 亚洲欧美精品自产自拍| 国产欧美日韩一区二区三区在线 | 亚洲自偷自拍三级| 免费大片黄手机在线观看| 下体分泌物呈黄色| 亚洲精品日韩av片在线观看| 尤物成人国产欧美一区二区三区| 五月玫瑰六月丁香| 亚洲精华国产精华液的使用体验| 一级二级三级毛片免费看| 看免费成人av毛片| 18禁在线无遮挡免费观看视频| 一级毛片黄色毛片免费观看视频| 大香蕉97超碰在线| 国产成人精品久久久久久| 精品一区二区三区视频在线| 我的老师免费观看完整版| 欧美精品亚洲一区二区| 少妇人妻久久综合中文| 日韩成人伦理影院| 五月天丁香电影| 亚洲美女视频黄频| 亚洲精品一二三| 99久国产av精品国产电影| 欧美bdsm另类| 蜜桃在线观看..| 毛片女人毛片| 免费播放大片免费观看视频在线观看| 蜜桃在线观看..| 亚洲av免费高清在线观看| 又大又黄又爽视频免费| 熟女av电影| 在线 av 中文字幕| 九九久久精品国产亚洲av麻豆| 国产伦在线观看视频一区| 下体分泌物呈黄色| 久久久久久久久大av| 国产在线视频一区二区| 久久久久国产精品人妻一区二区| 国产精品国产三级国产av玫瑰| 免费大片黄手机在线观看| 一级黄片播放器| xxx大片免费视频| 亚洲色图av天堂| 日本黄色片子视频| 日韩av免费高清视频| 国产av精品麻豆| 久久久久久久大尺度免费视频| 国产精品爽爽va在线观看网站| 亚洲欧美精品自产自拍| 亚洲图色成人| 在线观看美女被高潮喷水网站| 天天躁夜夜躁狠狠久久av| videos熟女内射| 免费观看无遮挡的男女| 国产极品天堂在线| 一级爰片在线观看|