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

    基于前攝效應(yīng)的不耐煩行為建模與雙向行人流動態(tài)*

    2022-04-15 07:32:46張琦渠靜
    物理學(xué)報 2022年7期
    關(guān)鍵詞:水平

    張琦 渠靜

    (北京交通大學(xué)交通運輸學(xué)院,北京 100044)

    不耐煩情緒是一種因人而異的心理因素,且隨環(huán)境實時變化.本研究構(gòu)建了基于元胞自動機的行人微觀模型,并用其研究了雙向行人流在不耐煩行為決策影響下的動態(tài)特征.模型定義了潛在勢能場表達行人對可用空間的占用傾向,合理體現(xiàn)雙向流行人決策的前攝效應(yīng).提出在行人運動過程中,根據(jù)瞬時狀態(tài)測量、記錄和更新行人個體不耐煩水平的方法,并以此觸發(fā)個性化行為.對比研究了3 種不同的不耐煩行為模式的作用效果.仿真試驗結(jié)果顯示,與不考慮不耐煩行為相比,當不耐煩水平用于觸發(fā)橫跨行為時,雙向流在速度、流率和等待時間方面的表現(xiàn)均有改善;當不耐煩用于直接計算行人對有限空間的競爭力時,這3 項指標則在大部分密度水平下均呈惡化.針對不耐煩作用強度與決策閾值的參數(shù)靈敏度分析顯示了模型在不耐煩作用表現(xiàn)方面的靈活性.同時,分場景的仿真試驗還揭示了低比例慢速行人、步道寬度及混合行為模式等對雙向流效率的作用.

    1 引言

    關(guān)于行人雙向流行為的既有研究揭示了行人群集現(xiàn)象源于個體的協(xié)作行為.針對現(xiàn)實世界中公共區(qū)域采集數(shù)據(jù)的調(diào)查和對情景實驗統(tǒng)計數(shù)據(jù)的研究表明:行人流的自組織,尤其是步道生成對雙向流起著至關(guān)重要的作用[1,2].隨著密度的增加,雙向流比單向流效率更高[3],這也可能是自組織效應(yīng)增強的結(jié)果[4,5].行人步道生成有助于有效緩解正面沖突[6],并防止形成僵局[7].因此,無論采用連續(xù)模型[8,9]還是離散模型[10-14],步道生成一直是行人雙向行為模型驗證的主要指標.然而,在不同情況下,行人對協(xié)作的態(tài)度是多變的.Nagatani[15]采用元胞自動機模型的實驗研究表明,行人在低密度時傾向自由移動,而高密度時則試圖生成步道.李明華等[16]構(gòu)建了改進的格子氣模型再現(xiàn)低密度下對向行人流不會出現(xiàn)分層現(xiàn)象的特性.Feliciani和Nishinari[17]對自組織步道生成的研究表明,步道生成后平衡的雙向行人流十分穩(wěn)定,然而步道生成的過程中則要求行人最大程度地橫向移動.這解釋了為什么平衡的雙向流在低密度下效率高,但會快速導(dǎo)致高密度下的死鎖.事實上,行人之間的協(xié)作是基于特定條件的,至少與密度[14,18,19]有關(guān).一些研究探討了高密度的建模解決方案[20,21],然而,問題的核心仍取決于協(xié)作機制是否正確合理地表達.

    在現(xiàn)實世界中,行人個體根據(jù)用戶最佳原則在微觀層面上做出運動決策,在宏觀層面上則體現(xiàn)為可能的系統(tǒng)最佳結(jié)果.例如,行人個體選擇適時跟隨前方的同向行人,避免與對向行人碰撞,這種利己的個體行為對應(yīng)的群體現(xiàn)象即為步道生成,且提高了特定條件下整個群體的效率.因此,行人之間的所謂“協(xié)作”,實際上是個體基于其局部環(huán)境的個性化行為的集成表現(xiàn).目前針對各種個性化行為的大量研究普遍集中在人群疏散場景上.Fu 等[22]引入了模糊理論表現(xiàn)性格各異的行人對危險信息的反應(yīng).Haghani 等[23]通過實驗揭示了適度的競爭壓力有利于中等規(guī)模人群的疏散,并引發(fā)了對傳統(tǒng)的“快即是慢”效應(yīng)的深入討論.Fu 等[24]和Weng等[25]則分別探討了行人差異化行為能力和速度對疏散效率的重要作用.針對行人流的研究也揭示了個性化行為的強大作用,包括隨機減速[26]、右行偏好[27]、偏移傾向[28]以及跟從行為[29-31]等.基于個性化行為的協(xié)作可能提升群體效率,也可能導(dǎo)致群體表現(xiàn)的惡化,尤其是在高密度爭奪有限空間[32,33]的條件下.關(guān)于人群疏散的研究表明,個體過于積極的疏散行為會導(dǎo)致整體效率降低[34,35].可見,群體現(xiàn)象的決定因素是個體的決策原則和行為模式,有必要正視個性化行為對群體表現(xiàn)的重要作用.

    個性化行為視角的引入打破了由統(tǒng)一運動規(guī)則觸發(fā)群體效應(yīng)的常規(guī)思路.在雙向行人流中,協(xié)作表現(xiàn)為基于個性化行為決策的步道生成過程.這里的個性化行為取決于行人個體對同向和對向行人潛在作用的反應(yīng),包括運動狀態(tài)下的碰撞規(guī)避和停滯狀態(tài)下的行為決策.碰撞規(guī)避基于行人對來自其他行人潛在作用的預(yù)判.Suma 等[36]和 Lü等[30]分別改進了勢能場模型和社會力模型用于表現(xiàn)行人的預(yù)見性行為.Jin和Guo[37]針對上下樓梯場景改進了元胞傳輸模型,利用異向行人對元胞勢能的影響驅(qū)動路徑選擇行為.這些研究揭示了前攝效應(yīng)在行人行為研究中的重要性.停滯狀態(tài)則源于排隊或沖突導(dǎo)致的擁堵,并進一步引發(fā)行人的不耐煩情緒.近年來不少行人流動態(tài)研究[38,39]中也考慮到了行人的心理狀況,不耐煩已被證明是影響人群疏散效率的一個關(guān)鍵因素[40,41].本研究假定不耐煩是可能引發(fā)雙向流中個性化行為的一個因素,并進一步導(dǎo)致不同的群體表現(xiàn).根據(jù)日常經(jīng)驗,在不耐煩的情況下,行人會更加渴望移動,更具競爭力或更加積極地尋找可用空間,從而也會影響其預(yù)判和碰撞規(guī)避的相關(guān)決策.值得注意的是,個體之間的不耐煩程度各不相同,而且不耐煩狀態(tài)會隨著行人周圍的情勢發(fā)生持續(xù)變化.因此,前攝效應(yīng)和不耐煩觸發(fā)的行為決策對雙向流的步道生成過程至關(guān)重要,進而影響雙向流的群體表現(xiàn).

    既有研究顯示了基于元胞自動機的模型[13-15,21,24-26,32,41-43]在捕捉與重現(xiàn)個性化行為方面具有強大的表現(xiàn)力和有效性.本文構(gòu)建基于元胞自動機的模型,通過定義潛在勢能場合里表現(xiàn)雙向流中行人決策的前攝效應(yīng).針對不耐煩觸發(fā)的個性化行為及其作用,提出了行人個體不耐煩水平的計量和更新規(guī)則,對比研究了由不耐煩觸發(fā)的3 種行為模式對雙向流的作用效果.

    2 行人雙向流中的前攝效應(yīng)與不耐煩行為

    沖突和死鎖是制約雙向行人行為建模的瓶頸問題之一.在真實世界中,行人并不會盲目前進直至擁堵,而是在行走過程中通過不斷觀察和調(diào)整避免沖突和停滯[14,33],這種前攝效應(yīng)在建模中尤為重要.行人在進行前攝決策時,同向和異向行人的相對數(shù)量發(fā)揮了潛在作用.通常同向行人相對較少的一方更傾向于轉(zhuǎn)向換道,而人數(shù)相對較多的一方則具備更明顯的優(yōu)先權(quán)保持原前進方向.

    如圖1 所示,行人群體Ln和Rn對向行走,行人在行走過程中實時表示自己對前方空間的占用傾向.kLn和kRn分別代表兩個群體的行人釋放占用傾向的作用距離.從時刻t過渡到時刻t+1 時,行人L1(身后沒有同向行人)感知到對向行人的占用傾向后隨即選擇換道至左側(cè)以避免與R1及其身后的人碰撞.

    圖1 雙向行人的前攝效應(yīng)Fig.1.Proactive effect of bi-directional pedestrians.

    本文提出的模型中定義了由行人作用和更新的潛在勢能場來表現(xiàn)上述前攝效應(yīng).每個行人個體對其前方空間的占用傾向以潛在勢能場的方式予以記錄、累加和比較.詳細規(guī)則見3.2 節(jié).

    前攝效應(yīng)的表達機制為緩解面對面沖突提供了途徑,但隨著人群密度的增加,擁堵不可避免.被迫停止和等待的行人可能產(chǎn)生不耐煩心理,這在真實世界中十分普遍.因此,不耐煩的行人可能基于其不耐煩的程度進行行為決策.一方面,不耐煩水平在行走過程中是動態(tài)變化的,當行人停止前進等待較長時間時,不耐煩水平較高,當其繼續(xù)以期望速度行走時不耐煩水平可能降低甚至消失.另一方面,不耐煩對行人的行為決策可能產(chǎn)生不同的影響.一個不耐煩的行人可能橫跨換道,也可能在占用前方有限空間方面更具競爭力.同向和異向不同個體之間差異化的動態(tài)不耐煩水平,使得每個行人面對的環(huán)境更加復(fù)雜和難以預(yù)測.本文針對3 種不同的不耐煩行為模式進行建模和對比研究.不耐煩水平的變化規(guī)則以及行為模式的表達詳見3.3 節(jié).

    3 不耐煩行人運動模型

    3.1 基礎(chǔ)定義

    仿真模型基于一個以方形單元格為基礎(chǔ)結(jié)構(gòu)的二維系統(tǒng).單元格邊長為 0.45 m,每個單元格在任一時間步只z 能由一名行人占用或者空閑.如圖2所示,有兩類行人以相反的方向相向運動,Lefts(L)代表“左向右”行人,Rights(R)代表“右向左”行人.每個行人作用并更新其影響的單元格上的潛在勢能場值,從而生成和更新潛在勢能場.潛在勢能場影響并觸發(fā)行人的移動決策.

    圖2 行人移動的方位及概率Fig.2.Pedestrian movement space and tendency.

    如圖2 所示,每個行人每時間步嘗試移動至其周圍的5 個格子,具體位置與其行走方向相關(guān).Ui,j表示行人移動至格子 (i,j) 的轉(zhuǎn)移概率.對于任意一個格子,行人對其設(shè)定和更新兩個潛在勢能場值:和,分別表示通過右向左行人和左向右行人進行賦值的潛在勢能場值.某個行人向其周圍相鄰格子的 (i,j) 的轉(zhuǎn)移概率Ui,j由該格子的和值確定.每時間步每個行人向其影響距離內(nèi)的所有格子設(shè)定潛在勢能場值,每個格子累計所有對其作用的行人的潛在勢能場值,并更新該時間步內(nèi)格子的和.下個時間步,每個行人依據(jù)自己周圍格子的和及其空閑情況,確定自己移動到哪個格子或原地等待,并進一步再次對其影響距離內(nèi)的所有格子設(shè)定潛在勢能場值.行走方向為d的行人的轉(zhuǎn)移概率如公(1)式所示.

    通過比較轉(zhuǎn)移概率,行人移動到周圍5 個格子中轉(zhuǎn)移概率最大的空閑單元格,轉(zhuǎn)移概率則取決于雙向行人的影響.行人移動可以通過持續(xù)更新潛在勢能場值和轉(zhuǎn)移概率值來實現(xiàn).潛在勢能場值更新的詳細規(guī)則見3.2 節(jié),行人移動規(guī)則見3.3 節(jié).模型的變量和參數(shù)見表1.

    表1 參數(shù)及其物理意義Table 1.Some parameters and the corresponding physical meanings.

    3.2 潛在勢能場值更新

    對于每個格子,兩個潛在勢能場值El和Er分別由可對其施加影響的“左向右”和“右向左”行人進行更新.行人可影響的格子包括其正在占據(jù)的格子以及未來時間步內(nèi)傾向于占據(jù)的格子.

    1)行人影響其正在占據(jù)的格子

    以圖3 中“左向右”行人L為例,一旦L移動到其前方3 個格子之一(圖3(b)和圖3(c)),目標格子的El將被行人L更新為El+1;如果L橫跨至左側(cè)或右側(cè)的格子(圖3(d))或停在目標格子不動(圖3(e)),目標格子的El將不會被行人L更新.因此,El的增加表示曾占用該格子的行人成功前進,該格子值得后方其他與行人L同向的行人考慮選擇.

    圖3 正在占據(jù)的格子的潛在勢能場值更新Fig.3.Potential field value update for occupied cells.

    2)行人影響其未來時間步內(nèi)傾向于占據(jù)的格子

    如果行人前進,前方k個格子為其未來時間步內(nèi)傾向于占據(jù)的格子;如果行人原地等待,將嘗試選擇左前或右前的格子,則其左前和右前格子為其未來時間步內(nèi)傾向于占據(jù)的格子.如圖4 所示,每個時間步內(nèi),一旦行人L移動到其前方3 個格子之一(圖4(b)和圖4(c)),或橫跨至左側(cè)或右側(cè)的格子(圖4(d)),則其前方k個格子的El將被行人L更新為El+1;如果L停在目標格子不動(圖4(e)),則其左前和右前格子的El將被行人L更新為El+1.同時,行人R前方k個格子的Er將被行人R更新為Er+1.因此,El的增加表示該格子對“左向右”行人更具吸引力.

    圖4 傾向于占據(jù)的格子的潛在勢能場值更新Fig.4.Potential field value update for interested cells.

    初始狀態(tài)下無人移動時每個格子的潛在勢能場值均為0.行人更新其正在占據(jù)的和傾向于占據(jù)的格子的潛在勢能場值.每個時間步內(nèi)無人占用的格子的El和Er衰減為El×β和Er×β,β為潛在勢能場值的衰減率.

    每個時間步內(nèi),每個人為自己影響范圍內(nèi)的格子設(shè)定潛在勢能場值,每個格子則依據(jù)所有對其產(chǎn)生影響的行人的作用進行累加并更新其潛在勢能場.對于某個格子,如果更多的“左向右”行人影響該格子,或更多的“左向右”行人占據(jù)并經(jīng)過該格子,則該格子的El應(yīng)當大于Er,對于“左向右”而言其轉(zhuǎn)移概率UEl-Er越大,反之亦然.

    3.3 行人移動規(guī)則

    3.3.1 不考慮不耐煩的行人移動規(guī)則運動過程中,行人依據(jù)轉(zhuǎn)移概率Ui,j向前方、左前和右前3 個格子中的一個移動.如果3 個格子均不空閑(如被他人占用),則原地等待直至這3 個格子有空閑.行人前進方向上前k個格子被視為傾向占用的影響區(qū)域.不考慮不耐煩的基本移動規(guī)則如表2 所列(以“右向左”行人及由其確定的Er為例).

    表2 不考慮不耐煩情況下“右向左”行人R 的基本運動規(guī)則Table 2.Basic movement rules a pedestrian of Rights without considering impatience.

    3.3.2 考慮不耐煩的行人移動規(guī)則

    移動過程中,如果3 個格子均不空閑,行人被迫原地等待,直至3 個格子中有格子空閑.不耐煩情緒在行人原地等待無法移動時產(chǎn)生.不同行人的不耐煩水平存在差異,對于同一個人,不耐煩水平也是變化的.行人的不耐煩水平隨著原地等待時間的增加而提升,一旦其再次開始移動,則不耐煩水平下降甚至消失.

    令it代表行人在時間步t的不耐煩水平.運動過程中,每個行人依據(jù)即時狀態(tài)更新其it.更新規(guī)則如下:

    如果行人不移動,令

    如果行人移動,令

    其中α(α≥ 1)為不耐煩水平增長率,γ(0 ≤γ≤1)為不耐煩水平衰減率.

    行人在不耐煩心理的作用下可能更積極地嘗試搜索空閑格子,可能在停滯不前時選擇橫跨,也可能引發(fā)更強烈的移動意愿.本文針對這3 種可能的不耐煩行為模式展開研究.

    模式1不耐煩水平i作為橫跨行為的觸發(fā)條件(圖5(a));

    模式2不耐煩水平i作為直接參數(shù)用于確定多人競爭同一目標格子時的優(yōu)先權(quán)(圖5(b));

    模式3不耐煩水平i用于計算轉(zhuǎn)移概率,并作為間接參數(shù)用于確定多人競爭同一目標格子時的優(yōu)先權(quán)(圖5(c)).

    令δ≥ 0表示影響決策的不耐煩水平閾值.如圖5(a)所示,“右向左”行人R不耐煩水平為i>δ時,該行人在前方3 個格子不空閑時決定橫跨換道,即不耐煩水平大于閾值作為橫跨行為的觸發(fā)條件(即模式1).圖5(b)中,目標格子的所有競爭者中,“左向右”行人L(不耐煩水平為i>δ)具備最高的不耐煩水平,則L可以優(yōu)先移動至該目標格子,即不耐煩水平i作為直接參數(shù)用于確定多人競爭同一目標格子時的優(yōu)先權(quán)(即模式2).不考慮不耐煩因素時采用隨機選擇規(guī)則處理沖突(表2),模式2 則引入不耐煩水平處理沖突.圖5(c)中,行人依據(jù)i×Ui,j計算轉(zhuǎn)移概率選擇目標格子,目標格子的所有競爭者中,“右向左”行人R(不耐煩水平為i>δ)具備最大的i×Ui,j值,則R可以優(yōu)先移動至該目標格子,即不耐煩水平i作為間接參數(shù)用于確定多人競爭同一目標格子時的優(yōu)先權(quán)(即模式3).

    圖5 3種不耐煩行為模式 (a) 模式1,不耐煩作為回避的觸發(fā)因素(R 不耐煩水平i > δ);(b) 模式2,不耐煩作為估計優(yōu)先權(quán)的直接參數(shù)(目標格子有4 名潛在競爭者,其中,L 的不耐煩水平i 最高且大于δ);(c) 模式3,不耐煩作為估計優(yōu)先級的間接參數(shù)(R 的i×Ui, j 最高且不耐煩水平i 大于δ)Fig.5.Three patterns of impatient behaviors:(a) Pattern 1,impatience as a trigger of sidestep (i of R > δ);(b) Pattern 2,impatience as a direct parameter to estimate priority (4 potential competitors for the target cell,L has the largest i (i of L > δ) among all of the potential competitors for the target cell);(c) Pattern 3,impatience as an indirect parameter to estimate priority (R has the largest i×Ui, j (i of R > δ) among all of the potential competitors for the target cell).

    以一個“右向左”行人為例,其考慮不耐煩行為的行動規(guī)則如表3—表5 所列.

    表3 模式1 的規(guī)則Table 3.Rules of Pattern 1.

    表4 模式2 的規(guī)則Table 4.Rules of Pattern 2.

    表5 模式3 的規(guī)則Table 5.Rules of Pattern 3.

    4 試驗結(jié)果分析

    采用多智能體可編程建模平臺NetLogo 6.0.3[44]搭建仿真試驗環(huán)境.通過設(shè)置不同的參數(shù)組合,在不同條件下實施雙向行人仿真試驗,用于分析不同場景下不耐煩行為模式的作用.N表示行人總數(shù)(雙向數(shù)量對等分別為N/2).定義密度為N/格子總數(shù),以符號ρ表示.T表示仿真時間步長.

    圖6 給出了隨著仿真時間步推進雙向流步道自組織形成過程.當T=0 時,一定數(shù)量的行人(ρ=0.3)隨機分布在22.05 m×22.95 m 的空間內(nèi),雙向行人數(shù)量對等.行人依據(jù)表2 所列的規(guī)則運動.當T=25 時,步道初步成型,T=100 時,步道形態(tài)更加清晰.可見,前攝型潛在勢能場的觸發(fā)和更新對于捕捉雙向行走過程中的自組織效應(yīng)是有效的.

    圖6 基本規(guī)則下的步道形成Fig.6.Lane formation with Basic movement rules (ρ=0.3).

    4.1 前攝型潛在勢能場及不耐煩效應(yīng)的作用分析

    仿真試驗針對3 個不耐煩行為模式及一個不考慮不耐煩的基本模式進行對比分析.圖7 給出了不同密度水平下行人觸發(fā)的潛在勢能場的動態(tài)情形.較深的顏色代表該格子具備較大的潛在勢能場值.圖中Basic 代表基本模式(不考慮不耐煩行為),Pattern 1—3 分別代表模式1—3 (考慮不耐煩行為).總體而言,基本模式下行人相對集中分布,同向行人由于排隊停滯而產(chǎn)生局部的輕微擁堵,并觸發(fā)了較強的潛在勢能場.考慮不耐煩的3 個模式下行人則相對分散,步道之間間隙更小,邊界也相對模糊.考慮不耐煩行為時,不同程度的擁堵主要源自對向行人的面對面沖突.

    圖8 給出了模式1 下不同密度水平行人個體的不耐煩水平.具有較高不耐煩水平的行人呈現(xiàn)較深的顏色.當密度較小時(ρ=0.1),沒有停滯和等待,不耐煩水平較低.密度較大時(ρ=0.5),擁堵較明顯,部分行人開始變得不耐煩,平均不耐煩水平提升.

    4.2 不耐煩行為模式及參數(shù)的作用分析

    4.2.1 不耐煩行為模式的作用分析

    用于結(jié)果分析的統(tǒng)計指標定義如表6 所列.每次仿真試驗前50 個時間步用于啟動試驗,指標的記錄和統(tǒng)計基于后800 時間步的情況.仿真試驗采用的參數(shù)取值為k=3,β=0.8,α=1.5,γ=0.5,δ=0.依據(jù)3.1 節(jié)的定義,試驗參數(shù)取值反映了前攝影響范圍為前方3 個格子,不耐煩水平直接影響決策,并且不耐煩作用力量為中等水平的情況.參數(shù)條件變化及對應(yīng)的靈敏度分析見4.2.2 節(jié),4.2.3節(jié)和4.3 節(jié).

    圖9 給出了速度(Speed)和流率(Flow rate)在不同密度水平下的變化曲線,本研究涉及的基本模式(Basic)和3 個不耐煩行為模式(Pattern 1—3分別代表模式1—3)與經(jīng)典的CA 模型(Blue,2001)進行對比.結(jié)果表明:基本模式與經(jīng)典CA模型的趨勢最為接近.而考慮不耐煩的3 種模式則表現(xiàn)出較大差異.模式1 顯示了與經(jīng)典CA 模型類似的速度曲線趨勢,但具有更高的絕對值.這樣的曲線形態(tài)揭示了合理的橫跨換道行為有助于提高雙向行人流的效率.模式2 在低密度和高密度下的速度和流率均較低,意味著如果不耐煩水平?jīng)Q定行人爭搶優(yōu)先空間的競爭力,則對行走效率產(chǎn)生負面影響.這與圖12 中步道較寬且邊界較為模糊的形態(tài)具有一定關(guān)聯(lián).模式3 在密度較低(ρ=0.2,0.3)時具備最高的速度和流率,但密度大于等于0.4 時卻突然惡化.這與圖12 中呈現(xiàn)的形態(tài)一致,ρ=0.3時行人有條不紊地移動,ρ=0.5 則出現(xiàn)明顯擁塞.

    圖9 不同密度水平下速度和流量的差異Fig.9.Comparison of speed and flow rate under different density levels.

    如圖10 所示,大多數(shù)密度水平下基本模式等待時間最長,表明不耐煩行為對減少行人的平均等待時間有積極的影響.模式1 在大多數(shù)密度下等待時間最短,在高密度時尤為明顯.模式2 在較低和較高密度下均具備3 種不耐煩模式中最長的等待時間.模式3 在較低的密度下獲得最短的等待時間,但隨著密度大于0.4 而快速增加.

    圖10 不同密度水平下的等待時間差異Fig.10.Wait time under different density levels.

    4.2.2 參數(shù)δ 的作用分析

    參數(shù)δ≥ 0表示不耐煩水平影響決策的閾值.δ=0表示只要行人不耐煩(i> 0),就依據(jù)相應(yīng)的模式規(guī)則觸發(fā)其不耐煩行為.δ> 0表示行人只有在不耐煩水平i達到特定閾值δ時,才會按照不耐煩模式規(guī)則觸發(fā)其行為.圖11 給出了在不同δ值作用下速度隨密度的變化關(guān)系.在不同的不耐煩模式下,δ對速度的影響明顯不同.δ=0 時模式1 具備所有密度水平下的最高速度,而δ> 0 時其具體取值對速度的影響不大.對于模式2,當ρ≤ 0.5時,δ> 0 會導(dǎo)致明顯的速度下降,但ρ=0.6 時速度值會隨δ值的增大而不同程度地增大.δ的取值對模式3 的速度影響不大.

    圖11 不同δ 情況下速度和密度的關(guān)系 (a) 模式1;(b) 模式2;(c) 模式3Fig.11.Speed vs.Density under different δ:(a) Pattern 1;(b) Pattern 2;(c) Pattern 3.

    4.2.3 參數(shù)α和γ的作用分析

    參數(shù)α和γ分別表示不耐煩水平的增長率和衰減率.α=1.2和γ=0.2表示不耐煩水平緩慢增加和快速下降,說明不耐煩的作用力量相對較弱.而α=1.8和γ=0.8 則代表不耐煩水平快速增加并緩慢下降,說明不耐煩的作用力量相對較強.

    對于模式1(圖12(a)),更強的不耐煩作用力得到更高的速度,反之,較弱的不耐煩作用力則得到較低的速度.α和γ不同取值時得到的等待時間差異不明顯(圖12(b)).對于模式2(圖12(c)和圖12(d)),最弱的不耐煩作用力(α=1.2和γ=0.2)在所有密度水平下均得到最低的速度和最長的等待時間.而最強的不耐煩作用力(α=1.8和γ=0.8)只在ρ≤ 0.4 時才能獲得最高速度和最短的等待時間.當ρ≥ 0.5 時,中等程度的不耐煩作用力表現(xiàn)更好.模式3 則反映出不同的特征(圖12(e)和圖12(f)).不耐煩作用力的增長和衰減對ρ=0.4—0.5 的速度和等待時間表現(xiàn)出相反的影響,即最弱的不耐煩作用力得到最高的速度和最長的等待時間,而最強的不耐煩作用力則得到最低的速度和最短的等待時間.

    圖12 不同 α和 γ 下的速度和等待時間隨密度的變化趨勢 (a) 模式1 的速度;(b) 模式1 的等待時間;(c) 模式2 的速度;(d) 模式2 的等待時間;(e) 模式3 的速度;(f) 模式3 的等待時間Fig.12.Speed and wait time vs.density under different α and γ:(a) Speed with Pattern 1;(b) wait time with Pattern 1;(c) speed with Pattern 2;(d) wait time with Pattern 2;(e) speed with Pattern 3;(f) wait time with Pattern 3.

    4.3 慢速行人作用分析

    考慮行人群體存在速度差異,部分行人(如老人或行動不便者)速度較低.與正常速度的行人相比,假設(shè)這類慢速行人對前方區(qū)域的影響范圍較小,且其不耐煩的作用力量相對較弱.本模型正常速度行人每時間步移動一格,即0.45 m,3 個時間步相當于1 s,走行速度約為1.35 m/s;慢速行人每2 個時間步移動一格,走行速度約為0.675 m/s.在仿真試驗中,慢速行人的行為參數(shù)取值為k=1,α=1.2,γ=0.2.正常速度行人參數(shù)取值為k=3,α=1.5,γ=0.5.圖13 給出了仿真試驗中連續(xù)時間步內(nèi)捕捉到的正常速度行人對慢速行人的超越過程.

    圖13 仿真試驗中正常速度行人對慢速行人的超越過程(局部圖) 注:其中三角形代表正常速度行人,箭頭形代表慢速行人Fig.13.Overtaking process between pedestrians with normal speed and slower speed (Partial screenshot) Note:the triangles represent pedestrians with normal speed and the arrows represent pedestrians with slower speed.

    設(shè)定慢速行人占總?cè)藬?shù)比例S分別為0%,10%和20%,試驗結(jié)果如圖14 所示.與全體行人均為單一正常速度(S=0%)相比,當10%或20%的行人為慢速行人時,行人的平均速度普遍降低,但慢速行人對不同的不耐煩作用模式的影響存在差異.模式1 的速度衰減在各密度水平下均較為顯著;模式2 的速度衰減在密度大于或等于0.5 時才較為明顯;而模式3 的速度衰減則在各密度水平下均相對平緩.說明在特定密度水平下,適度不耐煩行為對于緩解低比例慢速行人的沖擊具備積極作用.

    圖14 不同慢速行人比例下速度和密度的關(guān)系 (a) 模式1;(b) 模式2;(c) 模式3Fig.14.Speed vs.density under different proportion of slower pedestrians:(a) Pattern 1;(b) Pattern 2;(c) Pattern 3.

    4.4 步道寬度的作用分析

    為了研究不耐煩行為模式在不同步道寬度W條件下的表現(xiàn),選取3 種不同寬度進行仿真試驗.其中W=22.05 m 為理論試驗取值,依據(jù)《地鐵設(shè)計規(guī)范》(GB50157-2013)中規(guī)定島式站臺最小寬度為8 m,以及《鐵路旅客車站設(shè)計規(guī)范》(TB 10100-2018)中規(guī)定特大型站旅客進站、出站通道最小凈寬12 m,結(jié)合本模型的單元格尺寸0.45 m,分別取W=8.1 m,11.7 m;步道長度統(tǒng)一為22.95 m.

    實驗結(jié)果顯示,以W=22.05 m 的試驗結(jié)果為基礎(chǔ),隨著步道寬度的減小,各密度水平下的速度和流率基本呈現(xiàn)減小趨勢,密度高時速度和流率的衰減尤為明顯.當步道寬度減小為11.7 m 時,雙向流率降低了約44%—79%;當步道寬度減小為8.1 m時,雙向流率降低了約63%—86%.對比幾種行為模式,考慮不耐煩行為的模式1 在各密度水平下的平均速度衰減最為明顯,但兩個寬度水平下衰減程度相差不大,衰減率分別為38.4% (W=8.1 m)和36.7% (W=11.7 m);模式3 的速度衰減則對寬度較為敏感,衰減率分別為20.3% (W=8.1 m)和12.9%(W=11.7 m).

    4.5 不耐煩行為模式共存下的作用分析

    行人在產(chǎn)生不耐煩情緒時,可能會采取不同的行為模式.仿真試驗?zāi)M采用不同行為模式的行人共存條件下的雙向流行為,以研究模式共存條件的作用效果.圖15 中“Pattern 1+2”和“Pattern 1+3”分別表示行人中隨機選擇一半的人采用行為模式1,另一半人采用行為模式2 或3 的情形.如圖15 所示,在混合模式下,速度、流率和等待時間的曲線形態(tài)、最大值和臨界值均與單一模式不同.基本上,混合模式的各項統(tǒng)計值是兩個單一模式的中間值.

    圖15 混合模式不同密度下速度、流率、等待時間的變化 (a) 速度與密度關(guān)系;(b)流率與密度關(guān)系;(c) 等待時間與密度關(guān)系Fig.15.Speed,flow rate and wait time under different density levels with mix impatience patterns:(a) Speed;(b) flow rate;(c) wait time.

    5 結(jié)論

    本文構(gòu)建了不耐煩行為模型研究不耐煩行為對雙向行人流的影響.通過定義由個人產(chǎn)生和更新的潛在勢能場表現(xiàn)行人之間的前攝效應(yīng),并觸發(fā)行人的行為決策.在行人運動過程中,根據(jù)瞬時狀態(tài)測量、記錄和更新每個人的不耐煩水平,并通過可控參數(shù)調(diào)節(jié)不耐煩的變化動態(tài).模型支持對兼顧個體多樣性和運動過程動態(tài)性的不耐煩水平的實時表達,并將其作為行人個性化行為的觸發(fā)因素.本文提出并比較了3 種不同的不耐煩行為模式.

    基本圖和相關(guān)統(tǒng)計數(shù)據(jù)揭示了3 種不耐煩行為模式在不同情況下的作用效果.模式1 使雙向流具備最高的速度和流率、最短的等待時間;模式2 具備最低的速度和流率,且在大多數(shù)密度水平下等待時間最長;模式3 在較低的密度水平下具備最高的速度和流率以及最短的等待時間,但在高密度水平下各項指標都急劇惡化.不耐煩水平影響決策的閾值δ,對不同的不耐煩行為模式有不同的影響.由不耐煩水平增長率α和衰減率γ表達的不耐煩水平的作用力量對行人雙向流的效率也具備影響.通過調(diào)整前攝效應(yīng)和不耐煩水平的相關(guān)參數(shù),模型支持對慢速行人差異化不耐煩行為的表現(xiàn).在特定密度水平下,適度不耐煩行為對于緩解低比例慢速行人對雙向流效率的沖擊具備積極作用.同時,不同的不耐煩行為模式在步道寬度減小時的速度衰減程度也體現(xiàn)出差異性.在現(xiàn)實世界中,人群通常由不耐煩水平各異的行人組成,這些行人按照不同的不耐煩行為模式進行行為決策.本文中混合不耐煩模式場景的試驗揭示了雙向流基本圖的不同特征.更加全面的結(jié)論有賴于進一步深入的研究分析.

    本文構(gòu)建的行人模型為探索不耐煩行為對雙向行人流的作用提供了量化方法,試驗結(jié)果揭示了不耐煩行為模式對雙向流基本圖關(guān)鍵因素的作用,為理解個體個性化行為對群體運動效率的作用機理提供了理論支持.由于定量觀測行人的不耐煩行為特征較為困難,本文的研究局限于對特征規(guī)律的理論探討,后續(xù)研究工作中可借助虛擬現(xiàn)實等試驗手段獲取行人不耐煩心理及行為決策的動態(tài)特征,有望進一步標定和完善本文提出的模型.

    猜你喜歡
    水平
    張水平作品
    作家葛水平
    火花(2019年12期)2019-12-26 01:00:28
    深化精神文明創(chuàng)建 提升人大工作水平
    加強上下聯(lián)動 提升人大履職水平
    水平有限
    雜文月刊(2018年21期)2019-01-05 05:55:28
    加強自身建設(shè) 提升人大履職水平
    老虎獻臀
    中俄經(jīng)貿(mào)合作再上新水平的戰(zhàn)略思考
    建機制 抓落實 上水平
    中國火炬(2010年12期)2010-07-25 13:26:22
    做到三到位 提升新水平
    中國火炬(2010年8期)2010-07-25 11:34:30
    午夜福利欧美成人| 色在线成人网| 一个人免费在线观看的高清视频| 日本与韩国留学比较| 国产一级毛片七仙女欲春2| 国产亚洲精品一区二区www| 精品无人区乱码1区二区| 观看美女的网站| 国产免费av片在线观看野外av| 欧美日韩瑟瑟在线播放| 黄色视频,在线免费观看| 国产主播在线观看一区二区| 国产三级中文精品| 日韩成人在线观看一区二区三区| 脱女人内裤的视频| 色尼玛亚洲综合影院| 99国产精品一区二区蜜桃av| 亚洲熟妇中文字幕五十中出| 午夜两性在线视频| 成人一区二区视频在线观看| 久久九九热精品免费| 变态另类丝袜制服| 国内揄拍国产精品人妻在线| 999久久久精品免费观看国产| 熟妇人妻久久中文字幕3abv| 日日夜夜操网爽| 成年女人毛片免费观看观看9| 免费高清视频大片| 国产精品影院久久| 少妇高潮的动态图| 国产精品一区二区三区四区免费观看 | 免费电影在线观看免费观看| 中文字幕久久专区| 免费大片18禁| 我的老师免费观看完整版| 精品久久久久久成人av| 全区人妻精品视频| 99国产精品一区二区蜜桃av| 中文资源天堂在线| 噜噜噜噜噜久久久久久91| 可以在线观看毛片的网站| av欧美777| 欧美中文综合在线视频| 色综合亚洲欧美另类图片| 成人18禁在线播放| 搡女人真爽免费视频火全软件 | 乱人视频在线观看| 人人妻人人澡欧美一区二区| 99久久九九国产精品国产免费| 两个人看的免费小视频| 亚洲成人久久爱视频| 日韩欧美免费精品| 国产精品av视频在线免费观看| 国产av在哪里看| 欧美乱码精品一区二区三区| 久久久久精品国产欧美久久久| 无人区码免费观看不卡| 国产精品女同一区二区软件 | 久久国产精品影院| 岛国在线观看网站| 国产精品电影一区二区三区| 欧美一级毛片孕妇| 久久草成人影院| av片东京热男人的天堂| 一边摸一边抽搐一进一小说| 色综合亚洲欧美另类图片| 成人av一区二区三区在线看| 波野结衣二区三区在线 | 在线天堂最新版资源| 天天一区二区日本电影三级| 色播亚洲综合网| 亚洲av美国av| 午夜福利高清视频| 韩国av一区二区三区四区| 看片在线看免费视频| 精品国内亚洲2022精品成人| 国产单亲对白刺激| 18美女黄网站色大片免费观看| 国产一区二区激情短视频| 精品国产美女av久久久久小说| 婷婷精品国产亚洲av在线| 又紧又爽又黄一区二区| 国产激情欧美一区二区| 高清在线国产一区| 欧美成人性av电影在线观看| 少妇的逼好多水| 搡女人真爽免费视频火全软件 | 色播亚洲综合网| 亚洲av美国av| 国产主播在线观看一区二区| 麻豆一二三区av精品| 日本 欧美在线| 日韩精品青青久久久久久| 天堂影院成人在线观看| 亚洲精品影视一区二区三区av| 久久香蕉精品热| 欧美乱妇无乱码| 女人高潮潮喷娇喘18禁视频| 国产精品 欧美亚洲| 少妇的逼好多水| 国产毛片a区久久久久| 真人一进一出gif抽搐免费| www.999成人在线观看| 成人国产一区最新在线观看| 在线a可以看的网站| 国产精品亚洲美女久久久| 亚洲国产精品久久男人天堂| 可以在线观看的亚洲视频| 免费搜索国产男女视频| 99在线人妻在线中文字幕| 国产精品 国内视频| av欧美777| 日韩欧美精品v在线| 五月伊人婷婷丁香| 少妇人妻一区二区三区视频| 日韩欧美一区二区三区在线观看| 美女大奶头视频| 三级国产精品欧美在线观看| 久久久久国产精品人妻aⅴ院| 天堂av国产一区二区熟女人妻| 午夜日韩欧美国产| 免费av观看视频| 午夜日韩欧美国产| 内地一区二区视频在线| 少妇的丰满在线观看| 99riav亚洲国产免费| 久久中文看片网| 黄片小视频在线播放| 国产av不卡久久| 成年女人永久免费观看视频| 国产三级黄色录像| 99在线人妻在线中文字幕| 男女床上黄色一级片免费看| 十八禁网站免费在线| 国产国拍精品亚洲av在线观看 | 亚洲精品国产精品久久久不卡| 婷婷六月久久综合丁香| 大型黄色视频在线免费观看| 成人av一区二区三区在线看| 亚洲成人久久性| 91久久精品电影网| 熟女电影av网| 国产毛片a区久久久久| 日本熟妇午夜| 一个人看视频在线观看www免费 | 久久99热这里只有精品18| 18禁黄网站禁片免费观看直播| 国产午夜精品久久久久久一区二区三区 | 美女大奶头视频| 12—13女人毛片做爰片一| 国产99白浆流出| 动漫黄色视频在线观看| 国产亚洲精品久久久久久毛片| 变态另类成人亚洲欧美熟女| 精品福利观看| 老司机深夜福利视频在线观看| 人人妻人人看人人澡| 天堂影院成人在线观看| 99热这里只有精品一区| 亚洲一区高清亚洲精品| 18禁黄网站禁片午夜丰满| 1024手机看黄色片| 深夜精品福利| 国产视频一区二区在线看| h日本视频在线播放| 在线十欧美十亚洲十日本专区| 韩国av一区二区三区四区| 国产精品乱码一区二三区的特点| 九色国产91popny在线| 99久久成人亚洲精品观看| 精品午夜福利视频在线观看一区| 看免费av毛片| 欧美日本视频| 99视频精品全部免费 在线| 全区人妻精品视频| 午夜福利在线观看免费完整高清在 | 欧美高清成人免费视频www| 亚洲国产精品久久男人天堂| 国产 一区 欧美 日韩| 国产精品 欧美亚洲| 日本一本二区三区精品| 99久久精品一区二区三区| 欧美绝顶高潮抽搐喷水| 亚洲一区二区三区色噜噜| 白带黄色成豆腐渣| 最新美女视频免费是黄的| 日本黄色片子视频| 美女高潮的动态| 亚洲av免费高清在线观看| 一本久久中文字幕| 精品久久久久久久末码| 午夜老司机福利剧场| 嫁个100分男人电影在线观看| 久久天躁狠狠躁夜夜2o2o| 黄色片一级片一级黄色片| 亚洲成a人片在线一区二区| 国产视频内射| 九色国产91popny在线| 国产亚洲精品av在线| 国产免费男女视频| 夜夜看夜夜爽夜夜摸| 免费av不卡在线播放| 免费在线观看影片大全网站| 午夜福利高清视频| 国产不卡一卡二| 亚洲av成人av| 亚洲av成人av| 乱人视频在线观看| 狠狠狠狠99中文字幕| 日韩高清综合在线| 99在线视频只有这里精品首页| 亚洲人成伊人成综合网2020| 精品午夜福利视频在线观看一区| 久久国产精品人妻蜜桃| 免费一级毛片在线播放高清视频| 手机成人av网站| 国产黄片美女视频| 亚洲最大成人中文| 免费av不卡在线播放| 午夜福利在线观看免费完整高清在 | 51国产日韩欧美| 欧美日韩福利视频一区二区| 成人av在线播放网站| 欧美日韩一级在线毛片| 非洲黑人性xxxx精品又粗又长| 少妇的丰满在线观看| 亚洲人成网站在线播| 国产高清三级在线| 在线播放无遮挡| 成年版毛片免费区| 午夜两性在线视频| 国产精品综合久久久久久久免费| 男女午夜视频在线观看| 最好的美女福利视频网| 久久婷婷人人爽人人干人人爱| 国产精品,欧美在线| 我要搜黄色片| 成人av在线播放网站| www.熟女人妻精品国产| 午夜视频国产福利| 内射极品少妇av片p| 亚洲一区二区三区色噜噜| avwww免费| 我要搜黄色片| 可以在线观看的亚洲视频| 成人鲁丝片一二三区免费| 成年女人毛片免费观看观看9| 精品日产1卡2卡| 国产aⅴ精品一区二区三区波| 国产午夜精品论理片| 网址你懂的国产日韩在线| 精品电影一区二区在线| 国产91精品成人一区二区三区| 极品教师在线免费播放| а√天堂www在线а√下载| 欧美xxxx黑人xx丫x性爽| 无人区码免费观看不卡| 成人18禁在线播放| 国内精品久久久久精免费| 国产精品99久久久久久久久| 日韩欧美 国产精品| 国产精品免费一区二区三区在线| 91久久精品国产一区二区成人 | 国产欧美日韩一区二区精品| 日韩欧美在线二视频| 日日夜夜操网爽| 国产精品野战在线观看| 欧美区成人在线视频| 色在线成人网| 欧美日韩中文字幕国产精品一区二区三区| 久久亚洲真实| 精品一区二区三区视频在线 | 无限看片的www在线观看| 国产视频一区二区在线看| 亚洲专区中文字幕在线| 在线观看av片永久免费下载| 精品久久久久久,| 日本与韩国留学比较| 五月玫瑰六月丁香| 免费观看人在逋| 青草久久国产| 91麻豆av在线| 日日摸夜夜添夜夜添小说| 手机成人av网站| 91字幕亚洲| 精品国产美女av久久久久小说| 国产欧美日韩精品亚洲av| www.www免费av| 亚洲一区高清亚洲精品| 真实男女啪啪啪动态图| 麻豆一二三区av精品| 黄色成人免费大全| 亚洲专区中文字幕在线| 婷婷六月久久综合丁香| 90打野战视频偷拍视频| www.www免费av| 亚洲激情在线av| 18禁在线播放成人免费| 久久九九热精品免费| 久久草成人影院| 国产激情欧美一区二区| 国产精品99久久99久久久不卡| 午夜福利高清视频| 国产精品野战在线观看| 国产在视频线在精品| 亚洲中文日韩欧美视频| av福利片在线观看| 国产三级在线视频| 欧美性猛交黑人性爽| 国产高清有码在线观看视频| e午夜精品久久久久久久| 国产成人啪精品午夜网站| 国产亚洲精品综合一区在线观看| 中文字幕av成人在线电影| 久久精品国产清高在天天线| 18禁在线播放成人免费| 久久久国产成人免费| 欧美3d第一页| 日韩欧美在线乱码| 99国产极品粉嫩在线观看| 亚洲国产欧美人成| 久久性视频一级片| 成人av在线播放网站| 精品国内亚洲2022精品成人| 亚洲成人免费电影在线观看| 亚洲激情在线av| 夜夜爽天天搞| 最新在线观看一区二区三区| 国产av一区在线观看免费| 99久久99久久久精品蜜桃| 亚洲一区二区三区色噜噜| 老汉色av国产亚洲站长工具| 精品一区二区三区视频在线 | 精品福利观看| 又粗又爽又猛毛片免费看| 国产真人三级小视频在线观看| 51午夜福利影视在线观看| 亚洲欧美激情综合另类| 欧美乱色亚洲激情| 久久久久亚洲av毛片大全| 偷拍熟女少妇极品色| 色视频www国产| 日韩精品中文字幕看吧| 欧美在线黄色| 深爱激情五月婷婷| xxxwww97欧美| 欧美+日韩+精品| 久久久久亚洲av毛片大全| 淫秽高清视频在线观看| 一级黄片播放器| 国产成+人综合+亚洲专区| avwww免费| 亚洲aⅴ乱码一区二区在线播放| 色av中文字幕| h日本视频在线播放| 欧美又色又爽又黄视频| 黄色视频,在线免费观看| 欧美丝袜亚洲另类 | 国产视频一区二区在线看| 五月伊人婷婷丁香| 岛国在线免费视频观看| 99热这里只有精品一区| 婷婷精品国产亚洲av| 草草在线视频免费看| 久久天躁狠狠躁夜夜2o2o| 黄色丝袜av网址大全| 国产老妇女一区| 丰满的人妻完整版| 国产精品久久视频播放| 日本一本二区三区精品| 欧美日韩黄片免| 国产欧美日韩一区二区三| 成人无遮挡网站| 亚洲第一电影网av| 在线观看66精品国产| 国产蜜桃级精品一区二区三区| 狠狠狠狠99中文字幕| 在线免费观看不下载黄p国产 | 欧美一级a爱片免费观看看| 婷婷精品国产亚洲av| 黄色片一级片一级黄色片| 久久久色成人| 亚洲久久久久久中文字幕| 国产成人av激情在线播放| 日韩欧美在线乱码| 亚洲精品亚洲一区二区| 亚洲av美国av| 90打野战视频偷拍视频| 制服人妻中文乱码| 亚洲最大成人中文| 无遮挡黄片免费观看| 亚洲五月天丁香| 女人被狂操c到高潮| 我的老师免费观看完整版| 露出奶头的视频| 国产亚洲av嫩草精品影院| 欧美一区二区国产精品久久精品| 十八禁人妻一区二区| 黄色片一级片一级黄色片| 99热这里只有是精品50| 免费在线观看影片大全网站| 国产不卡一卡二| 丝袜美腿在线中文| 老司机午夜十八禁免费视频| 91麻豆av在线| 国产精品1区2区在线观看.| 欧美日韩黄片免| 亚洲美女视频黄频| 日韩高清综合在线| 欧美一级毛片孕妇| 国产精品一区二区三区四区免费观看 | 18+在线观看网站| a在线观看视频网站| 国产一区二区在线av高清观看| 国产熟女xx| 一卡2卡三卡四卡精品乱码亚洲| 久久久久久久午夜电影| 亚洲国产精品999在线| 91字幕亚洲| 级片在线观看| 久久久久精品国产欧美久久久| 99热6这里只有精品| 国产精品嫩草影院av在线观看 | 狂野欧美激情性xxxx| 日韩欧美国产在线观看| 成人一区二区视频在线观看| 久久久久免费精品人妻一区二区| 欧美+日韩+精品| 精品国内亚洲2022精品成人| 日本黄色片子视频| 高清毛片免费观看视频网站| 18禁美女被吸乳视频| 中文字幕人妻丝袜一区二区| 亚洲真实伦在线观看| 日本黄大片高清| 好男人电影高清在线观看| 欧美日韩黄片免| 动漫黄色视频在线观看| av欧美777| 免费观看人在逋| www.色视频.com| 有码 亚洲区| 亚洲电影在线观看av| 精品99又大又爽又粗少妇毛片 | 国产一级毛片在线| 99热6这里只有精品| 国产美女午夜福利| 日韩欧美精品免费久久| 一级二级三级毛片免费看| 久久精品综合一区二区三区| 欧美成人a在线观看| 国产黄片视频在线免费观看| 热99在线观看视频| 久久久午夜欧美精品| 身体一侧抽搐| 久久久久久久久久人人人人人人| 免费av不卡在线播放| 成人午夜精彩视频在线观看| 联通29元200g的流量卡| 亚洲精品乱码久久久v下载方式| 亚洲最大成人手机在线| 啦啦啦中文免费视频观看日本| 深爱激情五月婷婷| 国产av不卡久久| 2021天堂中文幕一二区在线观| 少妇的逼水好多| 欧美高清性xxxxhd video| av在线老鸭窝| 国产午夜精品一二区理论片| 亚洲精品一二三| 中文资源天堂在线| 在线播放无遮挡| 国产精品1区2区在线观看.| 日韩不卡一区二区三区视频在线| 2021少妇久久久久久久久久久| 精品人妻一区二区三区麻豆| 日韩av免费高清视频| av在线播放精品| 在现免费观看毛片| 水蜜桃什么品种好| 禁无遮挡网站| videossex国产| 欧美激情在线99| 日韩欧美精品免费久久| 精品亚洲乱码少妇综合久久| 亚洲色图av天堂| 精品午夜福利在线看| 久久久久精品性色| 久久久欧美国产精品| 国产极品天堂在线| 午夜亚洲福利在线播放| 亚洲精品国产av蜜桃| 三级经典国产精品| 亚洲熟女精品中文字幕| 青青草视频在线视频观看| 国产亚洲一区二区精品| 中文资源天堂在线| 哪个播放器可以免费观看大片| 午夜福利视频精品| 国产免费又黄又爽又色| 亚洲熟女精品中文字幕| 男人爽女人下面视频在线观看| 日韩在线高清观看一区二区三区| 中文天堂在线官网| 婷婷色综合大香蕉| 超碰av人人做人人爽久久| 三级毛片av免费| 久久久久久久久中文| 亚洲成人精品中文字幕电影| 校园人妻丝袜中文字幕| 亚洲av电影不卡..在线观看| 美女脱内裤让男人舔精品视频| 国产亚洲91精品色在线| 99久久精品国产国产毛片| av专区在线播放| 日日啪夜夜撸| 日本av手机在线免费观看| 精品久久久久久成人av| 你懂的网址亚洲精品在线观看| 国产国拍精品亚洲av在线观看| 天天躁日日操中文字幕| 免费人成在线观看视频色| 久久久久久九九精品二区国产| 综合色av麻豆| av线在线观看网站| 我的老师免费观看完整版| 国产综合精华液| 男女那种视频在线观看| 国产乱人视频| 18禁裸乳无遮挡免费网站照片| 免费看光身美女| 有码 亚洲区| 国产综合懂色| 黄片wwwwww| 水蜜桃什么品种好| 免费在线观看成人毛片| 插逼视频在线观看| 美女xxoo啪啪120秒动态图| 色尼玛亚洲综合影院| 国产黄色免费在线视频| 精品人妻熟女av久视频| 欧美激情久久久久久爽电影| 亚洲自偷自拍三级| 精品久久久久久久久久久久久| 日韩欧美 国产精品| 日本爱情动作片www.在线观看| 国产 一区 欧美 日韩| 午夜福利在线在线| 能在线免费看毛片的网站| 国产日韩欧美在线精品| 午夜福利在线在线| 国产成人a∨麻豆精品| a级一级毛片免费在线观看| 五月伊人婷婷丁香| av卡一久久| 久久久a久久爽久久v久久| 中文字幕av成人在线电影| 中国国产av一级| av线在线观看网站| 日本午夜av视频| 成人亚洲精品av一区二区| 亚洲欧美一区二区三区黑人 | 亚洲,欧美,日韩| 少妇熟女欧美另类| 黄色欧美视频在线观看| 蜜臀久久99精品久久宅男| 日本爱情动作片www.在线观看| 人妻少妇偷人精品九色| 亚洲在久久综合| 国产成年人精品一区二区| 九九爱精品视频在线观看| 身体一侧抽搐| 午夜亚洲福利在线播放| 国产大屁股一区二区在线视频| 成年女人看的毛片在线观看| 国产精品爽爽va在线观看网站| 国产在视频线在精品| 久久久久久久久大av| 亚洲图色成人| 亚洲精品日韩av片在线观看| 中文乱码字字幕精品一区二区三区 | 国产在线一区二区三区精| 草草在线视频免费看| 女人被狂操c到高潮| 国产一区二区在线观看日韩| 国产乱来视频区| 国产精品.久久久| 免费观看性生交大片5| 国产高清三级在线| 亚洲真实伦在线观看| 麻豆成人午夜福利视频| 国内精品一区二区在线观看| 日韩三级伦理在线观看| 亚洲人与动物交配视频| 一个人看视频在线观看www免费| 亚洲av福利一区| 久久久久精品性色| 国产伦一二天堂av在线观看| 久久久久久伊人网av| 又爽又黄无遮挡网站| 岛国毛片在线播放| 一边亲一边摸免费视频| 久热久热在线精品观看| 免费无遮挡裸体视频| 亚州av有码| 国产亚洲精品久久久com| 国产在线一区二区三区精| 成人亚洲欧美一区二区av| av在线播放精品| 黄色欧美视频在线观看| 国产黄a三级三级三级人| 久久久久久久久久久丰满| 精品不卡国产一区二区三区| 一级二级三级毛片免费看| 国产亚洲一区二区精品|