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

    步行設(shè)施內(nèi)疏散行人擁擠踩踏仿真研究

    2020-07-02 06:55:42劉秋梅武鑫森
    關(guān)鍵詞:倒地元胞摩擦力

    岳 昊,劉秋梅,武鑫森

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

    0 引 言

    隨著越來越多大型步行設(shè)施的建設(shè),行人擁擠誘發(fā)的踩踏事件時有發(fā)生,不僅帶來巨額的經(jīng)濟損失,還嚴重威脅人們的生命安全.擁擠踩踏事件具有突發(fā)性、危險性、不可再現(xiàn)性等特征,很難獲取踩踏事件發(fā)生時行人微觀動力學方面數(shù)據(jù),且組織行人擁擠踩踏實驗同樣具有一定的危險性.因此,基于微觀仿真的行人擁擠踩踏研究,成為行人流領(lǐng)域研究的熱點與難點.基于微觀仿真,探究擁擠踩踏的生成機理與時空分布等特征,有助于加強行人組織管理,優(yōu)化空間設(shè)施設(shè)計,提高行人疏散的安全性與效率性,強化行人疏散安全感.

    行人流疏散仿真分為:以社會力為代表的空間連續(xù)模型[1]和以元胞自動機為代表的空間離散模型[2].原始社會力模型存在“重疊”“群體特征偏少”等缺陷,后人從隨機行為[3]、從眾效應(yīng)[4]、信息傳遞等角度對其進行改進,但計算復雜、高密度人群仿真精度低的問題依然存在,不適合直接模擬擁擠人群踩踏現(xiàn)象.元胞自動機模型因空間分割相對任意性,用于疏散空間特征、擁擠致傷[5]、人群擠壓[6]和行人移動概率[7]等研究.元胞離散空間可以模擬擁擠時人群擠壓形變現(xiàn)象,卻很難體現(xiàn)行人之間力的微觀作用機制.本文采用相互力作用下的空間細化元胞自動機模型,利用元胞自動機空間離散特點,將空元胞作為行人倒地空間,即行人擁擠倒地誘發(fā)踩踏的初始條件,構(gòu)建行人間作用力的傳遞機制,再現(xiàn)離散空間下作用力的傳遞過程,達到與連續(xù)空間相同的作用力傳遞效果.在行人擁擠踩踏領(lǐng)域,與以往研究事故觸發(fā)機理、預(yù)防措施等不同,本文從微觀作用力角度探究行人踩踏內(nèi)在機理.

    1 踩踏機理

    以出口前疏散空間為例.疏散過程中,行人沿著指向安全出口的期望方向移動,在安全出口前產(chǎn)生擁擠排隊,每個行人迫切離開危險空間的意愿轉(zhuǎn)化為期望方向上的擁擠力.擁擠力造成行人之間的空間縮小和彈性擠壓,還在行人之間傳遞、累積,形成強大的破壞力;由于行人移動空間急劇縮小,在行人間接觸的條件下,擁擠力還會衍生出行人之間的摩擦力,從而降低行人的移動效率,增加行人的擁擠程度.在行人擁擠且競爭移動空間的過程中,擁擠力催促行人移動,摩擦力阻礙行人前進,在擁擠力與摩擦力相互作用下,行人無法按其意愿移動,移動速度大小和方向發(fā)生改變.當兩者作用力明顯加劇時,易使行人移動失敗,甚至失控倒地;此時,擁擠行人間相互作用力的平衡態(tài)被打破,使倒地者周邊行人的受力瞬間失衡,發(fā)生相繼失控倒地現(xiàn)象;由于倒地信息傳遞的滯后性,后方行人繼續(xù)向前擁擠,繼而發(fā)生連續(xù)踩踏現(xiàn)象.

    2 仿真模型

    2.1 空間與行人

    為模擬行人之間的擠壓形變現(xiàn)象,仿真模型基于疏散空間細分的元胞自動機構(gòu)建.根據(jù)文獻[8],行人在密集人群中占據(jù)0.4 m×0.4 m 的方形空間,近似為3×3 個元胞.由于行人可承受的擠壓形變程度有限,人體被視為由剛體部分和彈性部分構(gòu)成,如圖1(a)所示.行人剛體部分不可被他人侵占,保證行人在疏散過程中不會發(fā)生重疊;行人受力超過一定閾值時,彈性空間可與他人共享,如圖1(b)所示.基于空間細分,單位時間步長內(nèi),行人的移動領(lǐng)域為以行人為圓心半徑為3 個元胞的范圍,如圖1(c)所示.

    圖1 行人邊界與移動領(lǐng)域Fig.1 Pedestrian boundary and moving range

    當行人移動領(lǐng)域內(nèi)無空元胞時,移動期望方向指向安全出口;當行人移動領(lǐng)域內(nèi)有空元胞,且距離出口更近時,為尋求空間舒適性,行人移動期望方向指向空元胞,如圖2所示.

    2.2 擁擠力產(chǎn)生與更新

    由于安全出口的瓶頸效應(yīng),行人產(chǎn)生期望方向的移動意愿在擁擠阻滯情況下轉(zhuǎn)化為擁擠力Fn.行人移動選擇由移動收益決定,總是朝收益最大的元胞移動,由方向參數(shù)和空格參數(shù)兩部分構(gòu)成[2],即

    式中:Pij是移動收益;Dij是方向參數(shù);Sxy是點(x,y)到門口距離;Eij是空格參數(shù);(x0,y0)是門口坐標;(x,y)是行人當前坐標;(a,b)是行人目標點坐標;i,j分別表示行人橫縱坐標變化量,i=a-x,j=b-y.

    圖2 行人期望方向Fig.2 Expected direction

    由于行人的異質(zhì)性,假設(shè)擁擠力在0~Fmax間隨機分布,F(xiàn)max為最大擁擠力.擁擠力的作用效果具有連續(xù)性,需要在元胞離散化的基礎(chǔ)上再現(xiàn)擁擠力的連續(xù)作用效果.在仿真步T=t時刻,行人n的擁擠力為;在t→t+1 過程中,行人n的擁擠力變化量為ΔF;在T=t+1 時刻,理論上行人擁擠力既保留t時刻對受力對象的作用效果,也包含擁擠力變化量的實際作用效果的更新規(guī)則為

    2.3 擁擠力傳遞

    (1)受力對象.

    行人期望方向總是指向出口或距出口更近的元胞,故由行人移動意愿轉(zhuǎn)化而來的擁擠力的受力對象更接近安全出口.模型將行人編號及距安全出口的距離一一對應(yīng),按照與安全出口距離的大小從大到小進行遍歷.初始時刻第n個行人隨機賦值擁擠力,隨時間的推進,根據(jù)行人擁擠力的變化,計算下一時刻.根據(jù)擁擠力方向及行人位置找尋受力對象.受力對象滿足兩個基本條件:與施力主體接觸,距安全出口更近.

    (2)擁擠力動態(tài)傳遞分解.

    首先,尋找每個行人的受力對象,計算行人當前所受的擁擠力;然后,進行受力分解,直到遍歷空間內(nèi)所有行人.如圖3所示.

    圖3 擁擠力傳遞分解過程Fig.3 Push-force transmission process

    圖3中,吸收系數(shù)α為行人對擁擠力的吸收能力;抗擠傷系數(shù)β為行人對擁擠力的承受能力.α和β由行人自身體質(zhì)、年齡、穿著等決定.為行人輸出的擁擠力;為行人n所受擁擠力,i為對n施加擁擠力的行人;為行人n分解傳遞到行人m上的分力;θm為行人m與擁擠力的夾角.

    擁擠力有多個受力對象時,采用距離擁擠力方向最近的原則,將其分解傳遞給所有受力對象,如圖4 所示.擁擠力在傳遞過程中,由于行人自身彈性吸收作用,會被部分消化.其吸收、傳遞和分解規(guī)則[5]為

    圖4 擁擠力分解規(guī)則Fig.4 Decomposition rule

    2.4 擁擠踩踏

    (1)移動阻滯.

    行人在擁擠力的作用下產(chǎn)生擠壓形變,在移動過程中與他人產(chǎn)生相對速度,進而產(chǎn)生摩擦力.在擁擠狀態(tài)下,摩擦力不僅阻滯行人向前移動,還可能造成行人移動失敗甚至踩踏.行人摩擦力的大小由其所受擁擠力決定,其方向與行人運動方向呈180°,即

    式中:μ為摩擦系數(shù);為行人n所受擁擠力.

    (2)踩踏發(fā)生機理.

    為探究踩踏現(xiàn)象機理,引入抗摩擦系數(shù)γ,描述行人對摩擦力的克服作用.根據(jù)擁擠力與摩擦力的矛盾程度,判斷踩踏事故是否發(fā)生.當摩擦力與擁擠力矛盾達到一定值時,人群受力平衡點被打破,弱勢群體失控倒地,誘發(fā)踩踏.T=t時刻:當時,摩擦力不足以影響行人移動意愿,行人正常移動;當時,行人被迫移動失敗.T=t+1 時刻,若行人n再一次被選擇進行移動,且所受擁擠力,則該行人失控倒地.行人n的施力主體瞬間受力失衡,相繼倒地,直到施力主體受力在自身可承受范圍內(nèi)時,倒地結(jié)束,踩踏發(fā)生.踩踏事故發(fā)生過程如圖5所示.

    圖5 踩踏事故發(fā)生過程Fig.5 Process of stampede

    行人移動狀態(tài)為

    式中:S為行人移動狀態(tài),0和1分別表示正常移動和踩踏.

    2.5 仿真規(guī)則

    (1)移動沖突.

    若空間內(nèi)空元胞周圍行人滿足:移動后不與其他行人重疊,該空元胞在行人可移動范圍內(nèi),行人不會移動到距離相等或更遠的空元胞處,同一時間步內(nèi)同一行人最多移動一次,則其構(gòu)成該元胞的待移動行人集,如圖6所示.

    根據(jù)概率選擇模型確定最終移動行人.待移動行人集中所有行人競爭元胞空間,并阻礙彼此移動.行人移動概率由自驅(qū)動力、擁擠力和摩擦力共同決定.行人所受合力和移動概率為

    圖6 待移動行人集Fig.6 Pedestrians waiting to move

    式中:Fn為所受合力;為自驅(qū)動力;pn為選擇概率.

    (2)移動仿真流程.

    在行人疏散仿真過程中,首先計算疏散行人受力情況,然后根據(jù)行人受力情況進行移動選擇,Ns為踩踏行人,Nb為連續(xù)倒地行人,Nh為擠傷行人,仿真流程如圖7所示.

    圖7 疏散仿真過程Fig.7 Simulation process of evacuation

    3 仿真結(jié)果分析

    3.1 仿真場景

    擁擠踩踏往往發(fā)生在疏散瓶頸處,疏散人數(shù)與瓶頸寬度是主要的預(yù)防調(diào)控因素.行人擁擠踩踏是有效評價行人組織方案或瓶頸寬度設(shè)計合理性的主要指標.在預(yù)防擁擠踩踏的實時監(jiān)控中,主要監(jiān)控瓶頸前的疏散人數(shù).擁擠踩踏具有局部性、瓶頸處等特點,可在整體設(shè)施的基礎(chǔ)上,抽象出局部的、單瓶頸、非全局動態(tài)的仿真模型.本文以局部的、單出口疏散瓶頸為例,探究行人擁擠踩踏的生成機理,達到預(yù)防與保障安全的作用.

    仿真場景設(shè)置為尺寸8 m×8 m,出口寬度L=0.8 m,占據(jù)6個元胞的空間,如圖2所示.初始時刻行人呈拱形排列,待疏散行人數(shù)量N=202 人,人群密度為0.16 m2/人. 從行人屬性出發(fā),假設(shè)|Fmax|=100 單位,吸收系數(shù)α取[0,160]單位,抗擠傷系數(shù)β=170 單位,抗踩踏系數(shù)γ取[0,1 800]單位.擁擠疏散時行人步行速度為0.5 m/s,仿真步長為0.8 s.在此基礎(chǔ)上,在疏散空間設(shè)置4個柱狀障礙物,研究其是否影響疏散空間的安全性,如圖8所示.

    圖8 疏散空間布局Fig.8 Layout of obstacles

    3.2 行人踩踏情況

    圖9為N=202,L=6,不同α時,踩踏行人數(shù)量Ns隨γ的變化情況.當γ增加時,Ns先快速減小,然后緩慢達到趨于0的穩(wěn)定狀態(tài).隨著α增加,曲線簇由密集逐漸稀疏,故存在臨界吸收系數(shù)α1.當α≤α1時,此時α較小,吸收擁擠力的能力較弱,對踩踏人數(shù)變化的影響不明顯,α的增加不會帶來明顯的Ns的變化;當α>α1時,此時Ns的大小受α影響十分明顯,隨著α的增大,Ns下降的速度明顯增加.

    與無障礙物疏散空間相比,設(shè)置障礙物后,踩踏行人的數(shù)量明顯減少,如圖10所示.

    圖9 Ns 隨α 變化情況Fig.9 Change of Ns with α

    當γ在[0,1 800]變化時,總踩踏行人數(shù)量sumNs及踩踏與失衡倒地總?cè)藬?shù)sum(Ns+Nb)隨著α的變化情況如圖11(a)所示.sumNs隨α在[0,160]之間變化時,呈現(xiàn)先上升后下降趨勢;當α介于20~40時,達到最大值.sum(Ns+Nb)隨著α在[0,160]變化時,呈現(xiàn)單調(diào)遞減趨勢.

    圖10 空間布局對Ns 的影響Fig.10 Influence of spatial layout on Ns

    Ns隨著α的改變呈現(xiàn)先增加后減小趨勢,α=20 時出現(xiàn)峰值,如圖11(b)所示.此時,擁擠力有一定積聚時間,累積效果達到最大,爆發(fā)性更強.當γ在[0,1 800]變化時,Ns對γ的敏感程度逐漸減小,即存在臨界值γ1,決定Ns的變化幅度.當γ≤γ1時,隨著γ增加Ns迅速減小;當γ>γ1時,隨著γ增加Ns保持平穩(wěn)狀態(tài).

    圖11 sum Ns 、sum(Ns+Nb) 和Ns 隨α 變化情況Fig.11 Changes of sum Ns,sum(Ns+Nb) and Ns with α

    γ=700 時,行人擠傷次數(shù)隨著α增大不斷減少,整體分布在[0,50],如圖12(a)所示.高擠傷次數(shù)出現(xiàn)頻率受α變化影響較大,而低擠傷次數(shù)出現(xiàn)頻率對α變化不敏感.α=40 時,多次仿真得出行人平均擠傷情況,近80%的行人擠傷少于21次,隨擠傷次數(shù)增加,頻率呈現(xiàn)緩慢遞減狀態(tài),如圖12(b)所示.

    3.3 踩踏行人時空分布

    吸收系數(shù)α=120,抗擠傷系數(shù)β=170,抗踩踏系數(shù)γ=200 時,緊急疏散過程中踩踏行人分布情況如圖13 所示.踩踏行人均分布在疏散出口附近區(qū)域,時間t=8~20 s 是事故高發(fā)時段,其后僅偶爾出現(xiàn)踩踏現(xiàn)象.空元胞需要一定時間才能從門口傳遞到人群深處,故在疏散初期,后方行人分布還保留初始均勻狀態(tài),擁擠力均勻向前傳遞.隨著時間的推進,空元胞逐漸向后移動,行人排列呈現(xiàn)越來越明顯的無序化,擁擠力傳遞呈現(xiàn)不均勻性.局部行人聚集極易造成擁擠力的高度累積,處于擁擠累積鏈高處的行人有很大概率發(fā)生踩踏事故.在t=20 s 以后,部分行人疏離,以及行人踩踏倒地,使疏散空間內(nèi)行人數(shù)量大幅減少,閑暇空間增多,保證行人順利離開,踩踏事故較少發(fā)生.

    圖12 行人擠傷次數(shù)統(tǒng)計Fig.12 Statistics of pedestrian injury times

    圖13 踩踏行人分布圖Fig.13 Pedestrian distribution map

    4 結(jié) 論

    通過分析疏散人群間擁擠力和摩擦力的產(chǎn)生與作用機制探究踩踏發(fā)生機理.將連續(xù)與離散模型相結(jié)合,搭建力作用下的元胞自動機模型,再現(xiàn)行人擁擠踩踏現(xiàn)象.通過建立擁擠力的吸收、傳遞、摩擦力的產(chǎn)生與作用、行人概率選擇、事故發(fā)生等規(guī)則進行行人踩踏現(xiàn)象仿真.研究發(fā)現(xiàn):當多個行人競爭同一疏散空間時,局部人群矛盾沖突達到最大;當所受擁擠力和摩擦力超過安全閾值時,行人容易失控使局部受力失衡進而發(fā)生擁擠踩踏.在疏散瓶頸前設(shè)置障礙物,能有效中斷擁擠力的傳遞,減少擁擠踩踏發(fā)生.行人安全閾值越高外界擁擠力作用效果越低,擁擠人群穩(wěn)定性越高.因此,步行設(shè)施可以通過提高疏散空間安全系數(shù)和行人安全閾值避免踩踏,為事故防范提供建議.仿真結(jié)果發(fā)現(xiàn),疏散出口前是擁擠力高度聚集區(qū),擁擠踩踏空間上呈現(xiàn)環(huán)出口集中分布現(xiàn)象,時間上具有一定滯后性.本文模擬再現(xiàn)了行人擁擠踩踏現(xiàn)象,但未能有效再現(xiàn)實際擁堵中行人晃動現(xiàn)象,后續(xù)研究進一步挖掘行人擁擠力間的作用機理,再現(xiàn)擁堵晃動行人作用力的聚集效應(yīng).

    猜你喜歡
    倒地元胞摩擦力
    離奇的交通事故
    『摩擦力』知識鞏固
    理順摩擦力
    透析摩擦力
    騎車撞狗倒地受傷 飼養(yǎng)者有責賠四萬
    基于元胞自動機下的交通事故路段仿真
    智富時代(2018年5期)2018-07-18 17:52:04
    小偷被追身亡,失主該擔責嗎?
    公民導刊(2017年2期)2017-05-30 10:48:04
    神奇的摩擦力
    基于元胞數(shù)據(jù)的多維數(shù)據(jù)傳遞機制
    北京測繪(2016年2期)2016-01-24 02:28:28
    “土地爺”退休后“倒地”
    清風(2014年10期)2014-09-08 13:11:04
    另类亚洲欧美激情| 少妇精品久久久久久久| 九色亚洲精品在线播放| 老司机深夜福利视频在线观看 | 人妻 亚洲 视频| 国产亚洲精品第一综合不卡| 一级黄色大片毛片| 亚洲欧美精品综合一区二区三区| 丰满少妇做爰视频| 黄色a级毛片大全视频| av又黄又爽大尺度在线免费看| 成人国产av品久久久| 999久久久国产精品视频| 又紧又爽又黄一区二区| 99久久99久久久精品蜜桃| 男男h啪啪无遮挡| 首页视频小说图片口味搜索| 亚洲国产中文字幕在线视频| 欧美精品一区二区免费开放| 午夜福利在线观看吧| 男女边摸边吃奶| 亚洲av欧美aⅴ国产| 国产av精品麻豆| 美女高潮到喷水免费观看| 欧美激情 高清一区二区三区| 久久毛片免费看一区二区三区| 日本五十路高清| 国产精品免费视频内射| 国产免费一区二区三区四区乱码| 两个人看的免费小视频| 欧美人与性动交α欧美精品济南到| 丝袜人妻中文字幕| 亚洲精品一卡2卡三卡4卡5卡 | √禁漫天堂资源中文www| 国产一区二区三区av在线| 日韩视频一区二区在线观看| 免费看十八禁软件| 免费人妻精品一区二区三区视频| 国产男女超爽视频在线观看| 大陆偷拍与自拍| 亚洲av日韩在线播放| 操出白浆在线播放| 国产成人啪精品午夜网站| 久久影院123| 黄色片一级片一级黄色片| 国产免费福利视频在线观看| 国产免费av片在线观看野外av| 水蜜桃什么品种好| 久久午夜综合久久蜜桃| 桃花免费在线播放| 日韩 欧美 亚洲 中文字幕| 一个人免费看片子| 日韩 亚洲 欧美在线| 黑人猛操日本美女一级片| av福利片在线| 国产成人精品无人区| 精品乱码久久久久久99久播| 亚洲专区中文字幕在线| 午夜日韩欧美国产| 人妻久久中文字幕网| 一区在线观看完整版| a在线观看视频网站| 久久天躁狠狠躁夜夜2o2o| 亚洲av电影在线进入| 一二三四社区在线视频社区8| 亚洲国产看品久久| 国产精品影院久久| 飞空精品影院首页| 97精品久久久久久久久久精品| 我要看黄色一级片免费的| netflix在线观看网站| 丁香六月欧美| 欧美少妇被猛烈插入视频| 亚洲国产中文字幕在线视频| 操美女的视频在线观看| 欧美人与性动交α欧美精品济南到| 亚洲av国产av综合av卡| 精品免费久久久久久久清纯 | 伦理电影免费视频| 夜夜夜夜夜久久久久| 久久性视频一级片| 51午夜福利影视在线观看| 精品亚洲成a人片在线观看| 另类精品久久| 国产日韩欧美视频二区| 亚洲国产精品成人久久小说| 国产亚洲精品一区二区www | 高清视频免费观看一区二区| 国产精品香港三级国产av潘金莲| videos熟女内射| 在线亚洲精品国产二区图片欧美| 欧美成人午夜精品| 欧美日韩av久久| 在线 av 中文字幕| 丝袜美足系列| 亚洲精品久久午夜乱码| 亚洲国产精品一区三区| 黑人猛操日本美女一级片| 午夜久久久在线观看| 成年女人毛片免费观看观看9 | 精品久久久精品久久久| 在线观看免费高清a一片| 精品一区二区三区av网在线观看 | 1024视频免费在线观看| 亚洲国产精品一区二区三区在线| 精品国产乱码久久久久久男人| 91成年电影在线观看| 日韩熟女老妇一区二区性免费视频| 精品人妻熟女毛片av久久网站| 国产免费av片在线观看野外av| 久久狼人影院| 91成人精品电影| 精品人妻1区二区| 91字幕亚洲| 在线亚洲精品国产二区图片欧美| 国产精品久久久久久精品古装| 国产精品一二三区在线看| 狂野欧美激情性bbbbbb| 欧美人与性动交α欧美精品济南到| 视频在线观看一区二区三区| 欧美在线一区亚洲| 在线看a的网站| 亚洲国产精品一区二区三区在线| 成年美女黄网站色视频大全免费| 叶爱在线成人免费视频播放| 丝袜美腿诱惑在线| 欧美日韩亚洲综合一区二区三区_| 免费高清在线观看日韩| 老司机影院毛片| 国产男女超爽视频在线观看| 俄罗斯特黄特色一大片| 大陆偷拍与自拍| 国产有黄有色有爽视频| 国产av精品麻豆| 国产亚洲精品久久久久5区| 男女之事视频高清在线观看| 性色av乱码一区二区三区2| 欧美日韩亚洲综合一区二区三区_| 亚洲国产精品一区二区三区在线| 久久精品国产亚洲av高清一级| 国产免费福利视频在线观看| 国产福利在线免费观看视频| 亚洲欧美精品自产自拍| 天天操日日干夜夜撸| cao死你这个sao货| 最近最新免费中文字幕在线| 亚洲精品自拍成人| 久久精品成人免费网站| 国产黄色免费在线视频| 中文欧美无线码| 久久久国产成人免费| 99精品久久久久人妻精品| 国产欧美日韩综合在线一区二区| 热99re8久久精品国产| 国产区一区二久久| 老司机福利观看| 大片电影免费在线观看免费| 美女扒开内裤让男人捅视频| 国产精品欧美亚洲77777| 国产片内射在线| 久久久久久久久免费视频了| 啦啦啦中文免费视频观看日本| 男男h啪啪无遮挡| 一本综合久久免费| 欧美激情久久久久久爽电影 | 国产欧美亚洲国产| 国产在线视频一区二区| 免费高清在线观看日韩| 国产伦人伦偷精品视频| 天天操日日干夜夜撸| 桃红色精品国产亚洲av| 欧美日韩视频精品一区| 午夜老司机福利片| 99久久精品国产亚洲精品| 老司机午夜福利在线观看视频 | 视频在线观看一区二区三区| 男人操女人黄网站| 老司机午夜十八禁免费视频| 中文精品一卡2卡3卡4更新| 久久久精品94久久精品| 淫妇啪啪啪对白视频 | 蜜桃国产av成人99| 97在线人人人人妻| 亚洲国产看品久久| 高清视频免费观看一区二区| 久久久国产精品麻豆| 免费观看av网站的网址| 1024视频免费在线观看| 美女主播在线视频| 夜夜夜夜夜久久久久| 黑人猛操日本美女一级片| 亚洲欧美日韩高清在线视频 | 99久久综合免费| 曰老女人黄片| 妹子高潮喷水视频| 久久精品国产a三级三级三级| 男人添女人高潮全过程视频| 深夜精品福利| 99国产综合亚洲精品| 母亲3免费完整高清在线观看| 欧美成狂野欧美在线观看| 免费不卡黄色视频| 男女高潮啪啪啪动态图| 窝窝影院91人妻| 又黄又粗又硬又大视频| 国产精品一区二区精品视频观看| 亚洲精品国产区一区二| 精品亚洲成a人片在线观看| 欧美变态另类bdsm刘玥| 捣出白浆h1v1| 天堂中文最新版在线下载| 免费高清在线观看视频在线观看| 精品国产一区二区久久| 熟女少妇亚洲综合色aaa.| 久久久久久久精品精品| 欧美另类一区| 男人舔女人的私密视频| 无限看片的www在线观看| 午夜福利影视在线免费观看| 啦啦啦中文免费视频观看日本| 夜夜骑夜夜射夜夜干| 亚洲五月色婷婷综合| 美女视频免费永久观看网站| av在线老鸭窝| 深夜精品福利| 狠狠婷婷综合久久久久久88av| 久久精品成人免费网站| 久久这里只有精品19| www.熟女人妻精品国产| 午夜福利视频精品| 日韩 亚洲 欧美在线| 成年av动漫网址| 国产精品香港三级国产av潘金莲| 在线天堂中文资源库| 99国产精品一区二区蜜桃av | 亚洲欧美激情在线| 极品少妇高潮喷水抽搐| 日韩欧美一区二区三区在线观看 | 日韩制服骚丝袜av| 午夜久久久在线观看| 脱女人内裤的视频| 两性午夜刺激爽爽歪歪视频在线观看 | 国产日韩欧美视频二区| 日韩制服丝袜自拍偷拍| 亚洲自偷自拍图片 自拍| 欧美午夜高清在线| 80岁老熟妇乱子伦牲交| 在线看a的网站| 亚洲国产欧美在线一区| 免费在线观看日本一区| 亚洲专区国产一区二区| 国产亚洲av片在线观看秒播厂| 亚洲男人天堂网一区| 欧美精品一区二区大全| 亚洲综合色网址| 黑人猛操日本美女一级片| 一级毛片女人18水好多| av电影中文网址| 久久国产精品影院| 亚洲精品中文字幕在线视频| 男女午夜视频在线观看| 午夜免费观看性视频| 亚洲中文av在线| 亚洲一区二区三区欧美精品| 欧美人与性动交α欧美精品济南到| 爱豆传媒免费全集在线观看| 国产成人av激情在线播放| 黄片播放在线免费| 精品亚洲成国产av| 一二三四社区在线视频社区8| 亚洲av电影在线进入| 丰满人妻熟妇乱又伦精品不卡| 色播在线永久视频| 亚洲精品国产色婷婷电影| av电影中文网址| 每晚都被弄得嗷嗷叫到高潮| 80岁老熟妇乱子伦牲交| 日韩欧美免费精品| 美女主播在线视频| 欧美大码av| 精品国产国语对白av| 大片免费播放器 马上看| 亚洲精品日韩在线中文字幕| 一级毛片精品| 亚洲国产中文字幕在线视频| 91九色精品人成在线观看| 纯流量卡能插随身wifi吗| avwww免费| 免费观看av网站的网址| 日本wwww免费看| 一区二区三区四区激情视频| 黄网站色视频无遮挡免费观看| 精品一品国产午夜福利视频| 黄色视频,在线免费观看| 午夜精品久久久久久毛片777| 中国国产av一级| 国产精品免费视频内射| 啦啦啦免费观看视频1| 欧美日韩亚洲综合一区二区三区_| 99热国产这里只有精品6| 色视频在线一区二区三区| a 毛片基地| av不卡在线播放| 亚洲国产精品成人久久小说| 亚洲少妇的诱惑av| 日韩熟女老妇一区二区性免费视频| tocl精华| 国产高清国产精品国产三级| 欧美精品亚洲一区二区| 日本欧美视频一区| 少妇裸体淫交视频免费看高清 | 亚洲精品国产区一区二| 1024香蕉在线观看| 亚洲精华国产精华精| 熟女少妇亚洲综合色aaa.| 国产精品熟女久久久久浪| 97在线人人人人妻| 90打野战视频偷拍视频| 天堂中文最新版在线下载| 制服诱惑二区| 免费在线观看黄色视频的| 青春草亚洲视频在线观看| 91老司机精品| 精品福利观看| 国产淫语在线视频| 狠狠狠狠99中文字幕| av又黄又爽大尺度在线免费看| 美国免费a级毛片| 男女午夜视频在线观看| 国产成人精品在线电影| 午夜老司机福利片| 丝袜美足系列| 亚洲 国产 在线| 19禁男女啪啪无遮挡网站| 亚洲熟女精品中文字幕| 我的亚洲天堂| 啦啦啦啦在线视频资源| 亚洲av美国av| 成年人黄色毛片网站| 国产精品一区二区免费欧美 | 电影成人av| 亚洲欧洲日产国产| 水蜜桃什么品种好| 国产免费av片在线观看野外av| 最近最新免费中文字幕在线| 欧美另类一区| 亚洲精品第二区| 亚洲成人国产一区在线观看| 亚洲人成电影观看| 精品第一国产精品| 国产精品久久久久久精品古装| 日本黄色日本黄色录像| 99精品久久久久人妻精品| 国产成人免费无遮挡视频| 午夜激情av网站| 久久精品aⅴ一区二区三区四区| 母亲3免费完整高清在线观看| 高清在线国产一区| 亚洲五月色婷婷综合| 精品久久久久久久毛片微露脸 | 精品国产乱码久久久久久男人| 黑人巨大精品欧美一区二区蜜桃| 青青草视频在线视频观看| 99久久99久久久精品蜜桃| 亚洲av片天天在线观看| 国产精品一区二区在线不卡| 久久久精品免费免费高清| 午夜免费观看性视频| 在线精品无人区一区二区三| 国产在线免费精品| 这个男人来自地球电影免费观看| 男女高潮啪啪啪动态图| 亚洲欧美日韩另类电影网站| 国产xxxxx性猛交| 国产精品av久久久久免费| 69精品国产乱码久久久| 看免费av毛片| 久久青草综合色| 狠狠精品人妻久久久久久综合| 亚洲精品美女久久av网站| 国产国语露脸激情在线看| 91精品国产国语对白视频| www.自偷自拍.com| 91字幕亚洲| 亚洲视频免费观看视频| 人成视频在线观看免费观看| 国产精品免费大片| 国产片内射在线| 脱女人内裤的视频| 国产精品一区二区在线不卡| 欧美日韩亚洲国产一区二区在线观看 | 少妇 在线观看| 日本a在线网址| 午夜福利影视在线免费观看| 欧美激情高清一区二区三区| 亚洲精品国产av成人精品| 亚洲成人免费av在线播放| 国产亚洲av高清不卡| 日本av免费视频播放| www.999成人在线观看| 国产淫语在线视频| 欧美黑人精品巨大| 丝袜脚勾引网站| 在线观看舔阴道视频| 十八禁高潮呻吟视频| 国产成人一区二区三区免费视频网站| 91字幕亚洲| 日本wwww免费看| 精品人妻1区二区| 亚洲av男天堂| 美女午夜性视频免费| av又黄又爽大尺度在线免费看| 欧美av亚洲av综合av国产av| 精品欧美一区二区三区在线| 成年动漫av网址| 午夜久久久在线观看| 一区福利在线观看| 伦理电影免费视频| 十八禁网站免费在线| 亚洲欧美日韩高清在线视频 | 国产av精品麻豆| 久久精品国产综合久久久| 一个人免费在线观看的高清视频 | 高清视频免费观看一区二区| 亚洲精品av麻豆狂野| www.av在线官网国产| 丁香六月欧美| 日韩视频一区二区在线观看| 欧美人与性动交α欧美软件| 国产成人精品无人区| 亚洲国产精品成人久久小说| 久久天躁狠狠躁夜夜2o2o| 一二三四在线观看免费中文在| 久久天躁狠狠躁夜夜2o2o| 亚洲成人免费电影在线观看| tocl精华| 日韩制服丝袜自拍偷拍| 亚洲免费av在线视频| 欧美av亚洲av综合av国产av| 国产精品 欧美亚洲| 国产又色又爽无遮挡免| 看免费av毛片| 俄罗斯特黄特色一大片| 色婷婷久久久亚洲欧美| 一级,二级,三级黄色视频| 国产亚洲av高清不卡| 后天国语完整版免费观看| 欧美性长视频在线观看| 久久人人爽av亚洲精品天堂| 三上悠亚av全集在线观看| 久久av网站| 国产有黄有色有爽视频| tocl精华| 丰满人妻熟妇乱又伦精品不卡| 欧美少妇被猛烈插入视频| av网站在线播放免费| 亚洲欧美成人综合另类久久久| 久久久久视频综合| 欧美另类一区| 热re99久久国产66热| 夜夜骑夜夜射夜夜干| 热99国产精品久久久久久7| 国产97色在线日韩免费| 久久av网站| 国产免费av片在线观看野外av| 成人国产一区最新在线观看| 午夜精品国产一区二区电影| 亚洲国产毛片av蜜桃av| 国产av一区二区精品久久| 国产av一区二区精品久久| 国产av一区二区精品久久| 中文字幕最新亚洲高清| 女人精品久久久久毛片| 亚洲国产av新网站| 国产精品偷伦视频观看了| 日韩人妻精品一区2区三区| 日韩制服骚丝袜av| 国产精品久久久久成人av| 99精品欧美一区二区三区四区| 日韩大片免费观看网站| 亚洲精品一卡2卡三卡4卡5卡 | 在线av久久热| 青春草视频在线免费观看| 多毛熟女@视频| 99国产精品免费福利视频| 丁香六月欧美| 男女之事视频高清在线观看| 一二三四在线观看免费中文在| 国产精品成人在线| 国产日韩欧美亚洲二区| 一级片'在线观看视频| 蜜桃在线观看..| 黄色a级毛片大全视频| 国产在线免费精品| 久久久久网色| 久久九九热精品免费| 久久人人97超碰香蕉20202| 一本—道久久a久久精品蜜桃钙片| 在线观看免费视频网站a站| 搡老乐熟女国产| 欧美激情极品国产一区二区三区| 不卡av一区二区三区| 日韩人妻精品一区2区三区| 12—13女人毛片做爰片一| 精品国产乱码久久久久久小说| 国产野战对白在线观看| 午夜久久久在线观看| 捣出白浆h1v1| 咕卡用的链子| 亚洲激情五月婷婷啪啪| 中文字幕制服av| 777米奇影视久久| 日韩欧美一区二区三区在线观看 | 国产精品成人在线| 成人国产一区最新在线观看| 亚洲国产中文字幕在线视频| 少妇精品久久久久久久| 美女主播在线视频| 秋霞在线观看毛片| 国产在线免费精品| 69av精品久久久久久 | 中国美女看黄片| 亚洲国产欧美网| 国产xxxxx性猛交| 国产97色在线日韩免费| 日本五十路高清| 精品国产乱码久久久久久小说| 欧美日韩亚洲高清精品| 如日韩欧美国产精品一区二区三区| 国产一级毛片在线| 十八禁人妻一区二区| 午夜激情av网站| 色94色欧美一区二区| netflix在线观看网站| 成人免费观看视频高清| av国产精品久久久久影院| 日韩,欧美,国产一区二区三区| 国产片内射在线| 在线观看人妻少妇| 中文字幕人妻熟女乱码| 夜夜骑夜夜射夜夜干| 日韩三级视频一区二区三区| 久久午夜综合久久蜜桃| 久热这里只有精品99| 婷婷丁香在线五月| 亚洲精品国产一区二区精华液| 极品少妇高潮喷水抽搐| 亚洲色图 男人天堂 中文字幕| 曰老女人黄片| 日韩视频一区二区在线观看| 99久久国产精品久久久| 午夜两性在线视频| 男人操女人黄网站| 免费在线观看影片大全网站| 一进一出抽搐动态| 亚洲精品国产av蜜桃| 亚洲人成电影免费在线| 国产成人欧美| 国产熟女午夜一区二区三区| 性色av一级| 男人操女人黄网站| 欧美日韩亚洲高清精品| 一区二区av电影网| 男女午夜视频在线观看| 视频区图区小说| av天堂在线播放| 免费在线观看黄色视频的| 午夜视频精品福利| 如日韩欧美国产精品一区二区三区| 超碰97精品在线观看| 99香蕉大伊视频| 999久久久精品免费观看国产| 久久精品亚洲av国产电影网| 亚洲精品第二区| 少妇裸体淫交视频免费看高清 | 国产黄频视频在线观看| 欧美成人午夜精品| 丝袜人妻中文字幕| 欧美日韩黄片免| 这个男人来自地球电影免费观看| 一二三四在线观看免费中文在| 亚洲av国产av综合av卡| 午夜激情av网站| 亚洲人成电影观看| 国产精品国产三级国产专区5o| 十八禁人妻一区二区| 老汉色∧v一级毛片| 亚洲美女黄色视频免费看| 日韩中文字幕欧美一区二区| 黄片大片在线免费观看| 国产无遮挡羞羞视频在线观看| 亚洲天堂av无毛| 妹子高潮喷水视频| 香蕉国产在线看| 纯流量卡能插随身wifi吗| 国产成人一区二区三区免费视频网站| 午夜免费成人在线视频| 久久这里只有精品19| 涩涩av久久男人的天堂| 一本—道久久a久久精品蜜桃钙片| 99久久人妻综合| 日本撒尿小便嘘嘘汇集6| 午夜免费成人在线视频| 久久久国产一区二区| 精品一区二区三卡| 欧美日韩亚洲国产一区二区在线观看 | 国产免费av片在线观看野外av| 精品一区在线观看国产| 久久毛片免费看一区二区三区| 熟女少妇亚洲综合色aaa.| 精品人妻在线不人妻| 欧美另类亚洲清纯唯美| 99国产精品免费福利视频| 女人爽到高潮嗷嗷叫在线视频| 欧美精品人与动牲交sv欧美|