• <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
    精品少妇久久久久久888优播| 99热网站在线观看| 成人国产av品久久久| 制服诱惑二区| 黄色片一级片一级黄色片| 欧美精品一区二区免费开放| 精品熟女少妇八av免费久了| 亚洲国产中文字幕在线视频| 久久精品国产99精品国产亚洲性色 | 欧美精品一区二区免费开放| 免费观看a级毛片全部| 免费观看a级毛片全部| 搡老乐熟女国产| 9191精品国产免费久久| 人人妻,人人澡人人爽秒播| 免费在线观看完整版高清| 黄频高清免费视频| 精品熟女少妇八av免费久了| 中文字幕高清在线视频| 国产成人欧美| 午夜福利,免费看| 亚洲av欧美aⅴ国产| 菩萨蛮人人尽说江南好唐韦庄| 国产aⅴ精品一区二区三区波| 王馨瑶露胸无遮挡在线观看| 不卡一级毛片| 啦啦啦 在线观看视频| 高清毛片免费观看视频网站 | 女性生殖器流出的白浆| 免费观看a级毛片全部| 国产在线精品亚洲第一网站| 亚洲精品中文字幕在线视频| 午夜福利欧美成人| 中文字幕人妻丝袜一区二区| 最近最新免费中文字幕在线| 免费在线观看黄色视频的| 无人区码免费观看不卡 | 亚洲人成伊人成综合网2020| 国产成人免费观看mmmm| 久久国产精品人妻蜜桃| 91精品国产国语对白视频| 我要看黄色一级片免费的| 9色porny在线观看| 国产精品亚洲av一区麻豆| 大码成人一级视频| av国产精品久久久久影院| 国产成人av激情在线播放| 2018国产大陆天天弄谢| 正在播放国产对白刺激| 亚洲中文日韩欧美视频| 啦啦啦在线免费观看视频4| 热re99久久国产66热| 色尼玛亚洲综合影院| av一本久久久久| 欧美成人免费av一区二区三区 | 两个人免费观看高清视频| 亚洲成人免费av在线播放| 国产精品免费一区二区三区在线 | 亚洲国产欧美网| 黑人操中国人逼视频| 国产老妇伦熟女老妇高清| 亚洲色图av天堂| 婷婷成人精品国产| 欧美精品高潮呻吟av久久| 国产一区二区三区综合在线观看| 欧美亚洲日本最大视频资源| 亚洲国产欧美日韩在线播放| 亚洲av成人不卡在线观看播放网| 亚洲自偷自拍图片 自拍| 亚洲欧美一区二区三区久久| 老司机在亚洲福利影院| 一边摸一边抽搐一进一出视频| 久久婷婷成人综合色麻豆| 男人舔女人的私密视频| 亚洲,欧美精品.| 国产成人精品无人区| 欧美激情极品国产一区二区三区| 十八禁高潮呻吟视频| 免费黄频网站在线观看国产| 久久狼人影院| 国产在线视频一区二区| 一个人免费在线观看的高清视频| 国内毛片毛片毛片毛片毛片| 视频在线观看一区二区三区| 久久久久视频综合| 亚洲国产欧美在线一区| 窝窝影院91人妻| 香蕉久久夜色| 窝窝影院91人妻| 一二三四在线观看免费中文在| av免费在线观看网站| 国产精品熟女久久久久浪| 999久久久国产精品视频| 国产日韩欧美亚洲二区| 男女无遮挡免费网站观看| 99在线人妻在线中文字幕 | 老司机在亚洲福利影院| 每晚都被弄得嗷嗷叫到高潮| 在线观看66精品国产| 久久精品国产a三级三级三级| 日韩欧美一区二区三区在线观看 | 亚洲成a人片在线一区二区| 精品免费久久久久久久清纯 | 免费在线观看日本一区| 国产精品1区2区在线观看. | 韩国精品一区二区三区| 国产av一区二区精品久久| 欧美日韩亚洲综合一区二区三区_| 久久久久精品人妻al黑| 高清黄色对白视频在线免费看| 日本精品一区二区三区蜜桃| 1024香蕉在线观看| 中国美女看黄片| 最黄视频免费看| 色尼玛亚洲综合影院| 人人妻人人添人人爽欧美一区卜| 女警被强在线播放| 日韩一卡2卡3卡4卡2021年| 欧美精品av麻豆av| 激情在线观看视频在线高清 | 两个人看的免费小视频| 老司机福利观看| 少妇猛男粗大的猛烈进出视频| 一二三四在线观看免费中文在| 亚洲少妇的诱惑av| 国产精品一区二区在线观看99| 99国产综合亚洲精品| 青草久久国产| 亚洲人成电影观看| 国产欧美日韩一区二区精品| 国产精品久久久久久人妻精品电影 | 亚洲色图av天堂| 老熟妇仑乱视频hdxx| 精品人妻在线不人妻| 男女床上黄色一级片免费看| 一边摸一边抽搐一进一小说 | 亚洲成人国产一区在线观看| 人妻 亚洲 视频| 日本精品一区二区三区蜜桃| 91国产中文字幕| 丝袜人妻中文字幕| 亚洲欧洲日产国产| 午夜91福利影院| 午夜视频精品福利| 国产高清激情床上av| 日韩大码丰满熟妇| 中文字幕高清在线视频| 一个人免费看片子| 少妇裸体淫交视频免费看高清 | 涩涩av久久男人的天堂| 久久99一区二区三区| 亚洲精品乱久久久久久| 一本综合久久免费| 国产高清激情床上av| 制服诱惑二区| 色播在线永久视频| 国产xxxxx性猛交| 男女免费视频国产| 9色porny在线观看| 国产深夜福利视频在线观看| 亚洲性夜色夜夜综合| 成人黄色视频免费在线看| 91大片在线观看| 国产日韩欧美视频二区| 99re6热这里在线精品视频| 亚洲精品自拍成人| 青青草视频在线视频观看| 超碰97精品在线观看| 欧美黑人欧美精品刺激| 精品国产一区二区久久| 夜夜骑夜夜射夜夜干| 久久久久久久大尺度免费视频| 国产亚洲欧美在线一区二区| 国产精品亚洲av一区麻豆| 欧美性长视频在线观看| 午夜福利欧美成人| 午夜老司机福利片| 国产精品国产高清国产av | 亚洲成av片中文字幕在线观看| 国产成人一区二区三区免费视频网站| 一区二区三区乱码不卡18| 国产亚洲精品第一综合不卡| 亚洲av第一区精品v没综合| 婷婷丁香在线五月| 亚洲三区欧美一区| 成人亚洲精品一区在线观看| 制服人妻中文乱码| 久久久久久久国产电影| 最近最新中文字幕大全免费视频| 丁香六月欧美| 欧美在线黄色| 一区在线观看完整版| 日日摸夜夜添夜夜添小说| 9热在线视频观看99| 国产精品二区激情视频| 国产欧美日韩综合在线一区二区| 免费一级毛片在线播放高清视频 | 在线观看66精品国产| 宅男免费午夜| 国产精品美女特级片免费视频播放器 | 午夜福利影视在线免费观看| 久久亚洲精品不卡| 欧美乱妇无乱码| 咕卡用的链子| av在线播放免费不卡| 久久久久国内视频| 美女福利国产在线| 久久精品熟女亚洲av麻豆精品| 国产av又大| 热re99久久国产66热| 亚洲国产毛片av蜜桃av| 最近最新免费中文字幕在线| 一本—道久久a久久精品蜜桃钙片| 久久久久网色| 国产av精品麻豆| av又黄又爽大尺度在线免费看| 老司机在亚洲福利影院| 深夜精品福利| 欧美黄色淫秽网站| 欧美成人午夜精品| 老汉色av国产亚洲站长工具| 日本wwww免费看| 国产老妇伦熟女老妇高清| 热99久久久久精品小说推荐| 99精品久久久久人妻精品| 国产精品熟女久久久久浪| 一区二区日韩欧美中文字幕| 大香蕉久久成人网| netflix在线观看网站| 一边摸一边抽搐一进一出视频| 国产精品免费一区二区三区在线 | 国产高清视频在线播放一区| 国产男女超爽视频在线观看| 高清av免费在线| 曰老女人黄片| 国产在线免费精品| 欧美一级毛片孕妇| 亚洲 国产 在线| 免费观看av网站的网址| 午夜激情久久久久久久| 亚洲少妇的诱惑av| 免费人妻精品一区二区三区视频| 香蕉久久夜色| 日本精品一区二区三区蜜桃| 免费av中文字幕在线| 亚洲精品在线美女| 男女下面插进去视频免费观看| 最新美女视频免费是黄的| 一个人免费看片子| 欧美日韩视频精品一区| 我要看黄色一级片免费的| 亚洲精品一二三| 视频区图区小说| 色综合婷婷激情| 在线永久观看黄色视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲午夜精品一区,二区,三区| 精品少妇一区二区三区视频日本电影| 免费高清在线观看日韩| 国产精品久久电影中文字幕 | 国产亚洲欧美精品永久| 亚洲国产欧美日韩在线播放| 欧美日本中文国产一区发布| 免费久久久久久久精品成人欧美视频| 国产在线一区二区三区精| 少妇裸体淫交视频免费看高清 | 人妻 亚洲 视频| av片东京热男人的天堂| 欧美精品啪啪一区二区三区| 人妻 亚洲 视频| 欧美日韩福利视频一区二区| av免费在线观看网站| 久久中文字幕一级| 久久久久久久久免费视频了| 久热爱精品视频在线9| 波多野结衣av一区二区av| 1024视频免费在线观看| 男男h啪啪无遮挡| 美女午夜性视频免费| 久久九九热精品免费| 国产高清videossex| 国产精品1区2区在线观看. | 岛国在线观看网站| 国产亚洲一区二区精品| 国产野战对白在线观看| 老汉色av国产亚洲站长工具| 精品福利永久在线观看| 18禁裸乳无遮挡动漫免费视频| 午夜日韩欧美国产| a级毛片黄视频| 制服人妻中文乱码| 少妇 在线观看| 天天躁夜夜躁狠狠躁躁| 精品乱码久久久久久99久播| 波多野结衣一区麻豆| 两个人看的免费小视频| 建设人人有责人人尽责人人享有的| 亚洲国产看品久久| 亚洲av第一区精品v没综合| 女人精品久久久久毛片| 美女扒开内裤让男人捅视频| 欧美性长视频在线观看| 久久天堂一区二区三区四区| 中亚洲国语对白在线视频| 欧美另类亚洲清纯唯美| 热99国产精品久久久久久7| 在线观看人妻少妇| 最近最新中文字幕大全免费视频| 国产单亲对白刺激| 亚洲va日本ⅴa欧美va伊人久久| 免费黄频网站在线观看国产| 国产精品自产拍在线观看55亚洲 | 大型黄色视频在线免费观看| 视频区欧美日本亚洲| 久久久久久亚洲精品国产蜜桃av| 国产高清视频在线播放一区| 免费人妻精品一区二区三区视频| 黑人巨大精品欧美一区二区mp4| 极品教师在线免费播放| 精品国产乱子伦一区二区三区| 999精品在线视频| 性色av乱码一区二区三区2| 久久久久国产一级毛片高清牌| 757午夜福利合集在线观看| 深夜精品福利| 国产精品 欧美亚洲| 欧美在线黄色| 69精品国产乱码久久久| 中文字幕高清在线视频| 久久中文看片网| 欧美一级毛片孕妇| 97在线人人人人妻| 亚洲免费av在线视频| 日本黄色视频三级网站网址 | 99久久精品国产亚洲精品| 亚洲精品一卡2卡三卡4卡5卡| 国产福利在线免费观看视频| 亚洲午夜理论影院| 极品教师在线免费播放| 免费在线观看日本一区| 捣出白浆h1v1| 正在播放国产对白刺激| 亚洲第一青青草原| 大片电影免费在线观看免费| 日本欧美视频一区| 国产一区二区在线观看av| 中文字幕人妻丝袜制服| 18禁裸乳无遮挡动漫免费视频| a级毛片黄视频| 亚洲国产精品一区二区三区在线| 在线观看舔阴道视频| 亚洲成人国产一区在线观看| 91大片在线观看| 波多野结衣av一区二区av| 色婷婷av一区二区三区视频| 下体分泌物呈黄色| 国产精品自产拍在线观看55亚洲 | 99热网站在线观看| 亚洲全国av大片| 久久人妻福利社区极品人妻图片| 一区二区三区国产精品乱码| 大陆偷拍与自拍| av又黄又爽大尺度在线免费看| 国产一区二区三区视频了| cao死你这个sao货| 99久久99久久久精品蜜桃| 久久中文看片网| 国产高清国产精品国产三级| 国产真人三级小视频在线观看| 国产97色在线日韩免费| 欧美精品亚洲一区二区| 国产精品九九99| 人人妻人人添人人爽欧美一区卜| 18禁国产床啪视频网站| 国产精品一区二区免费欧美| 欧美午夜高清在线| 两个人看的免费小视频| 日韩欧美一区二区三区在线观看 | 国产av又大| 999久久久精品免费观看国产| 电影成人av| 亚洲精品一卡2卡三卡4卡5卡| 免费观看av网站的网址| 成年人午夜在线观看视频| 精品国产一区二区三区久久久樱花| 大码成人一级视频| 国产男靠女视频免费网站| 制服诱惑二区| 久久久久网色| 一边摸一边抽搐一进一小说 | 9热在线视频观看99| 日韩欧美免费精品| 一本久久精品| 日韩视频一区二区在线观看| 又黄又粗又硬又大视频| 欧美日韩亚洲综合一区二区三区_| 天天躁日日躁夜夜躁夜夜| 国产三级黄色录像| 丝袜美足系列| 五月开心婷婷网| 黑人猛操日本美女一级片| 国产在线观看jvid| 国产亚洲欧美精品永久| 国产日韩欧美在线精品| 欧美黄色淫秽网站| 新久久久久国产一级毛片| 桃花免费在线播放| 亚洲av成人不卡在线观看播放网| 久久这里只有精品19| 久久久久久久国产电影| e午夜精品久久久久久久| 午夜精品国产一区二区电影| 黄色片一级片一级黄色片| 黄色 视频免费看| 别揉我奶头~嗯~啊~动态视频| 男人舔女人的私密视频| 亚洲精品乱久久久久久| 午夜福利一区二区在线看| 久久人妻av系列| 国产在线视频一区二区| 欧美日韩国产mv在线观看视频| 免费观看av网站的网址| 一级a爱视频在线免费观看| 汤姆久久久久久久影院中文字幕| 女人被躁到高潮嗷嗷叫费观| 桃花免费在线播放| 久久精品aⅴ一区二区三区四区| 久久久精品94久久精品| 一区二区日韩欧美中文字幕| 性少妇av在线| 国产日韩一区二区三区精品不卡| 12—13女人毛片做爰片一| 国产精品av久久久久免费| 怎么达到女性高潮| 涩涩av久久男人的天堂| 91国产中文字幕| 午夜激情av网站| 亚洲精品粉嫩美女一区| 色综合欧美亚洲国产小说| 美女国产高潮福利片在线看| 国产午夜精品久久久久久| 人妻久久中文字幕网| 黄色 视频免费看| 性色av乱码一区二区三区2| 亚洲午夜理论影院| 精品视频人人做人人爽| 亚洲色图 男人天堂 中文字幕| 欧美日韩av久久| 国产欧美亚洲国产| 18禁裸乳无遮挡动漫免费视频| 麻豆成人av在线观看| 成年动漫av网址| 1024视频免费在线观看| 一边摸一边抽搐一进一出视频| 国产男女超爽视频在线观看| 久久久久精品国产欧美久久久| 天天躁狠狠躁夜夜躁狠狠躁| 久久性视频一级片| 精品卡一卡二卡四卡免费| 宅男免费午夜| 十分钟在线观看高清视频www| 黄网站色视频无遮挡免费观看| 国产成人精品久久二区二区91| 国产一区二区 视频在线| 波多野结衣av一区二区av| 精品一区二区三区四区五区乱码| 亚洲欧美日韩另类电影网站| 成在线人永久免费视频| 一级毛片女人18水好多| 99re6热这里在线精品视频| 高清在线国产一区| 啦啦啦 在线观看视频| 亚洲精品久久午夜乱码| 一个人免费在线观看的高清视频| 黑人操中国人逼视频| 欧美精品一区二区大全| 又紧又爽又黄一区二区| 国产国语露脸激情在线看| 91字幕亚洲| av又黄又爽大尺度在线免费看| 国产又爽黄色视频| 久久国产精品人妻蜜桃| 91精品三级在线观看| netflix在线观看网站| 久久中文字幕人妻熟女| 欧美人与性动交α欧美软件| 99国产精品一区二区蜜桃av | 亚洲avbb在线观看| 91麻豆精品激情在线观看国产 | 欧美亚洲 丝袜 人妻 在线| 国产视频一区二区在线看| 欧美精品一区二区免费开放| 狠狠精品人妻久久久久久综合| 中文字幕人妻熟女乱码| 亚洲伊人色综图| 日本撒尿小便嘘嘘汇集6| 精品一区二区三卡| 日韩中文字幕视频在线看片| av在线播放免费不卡| 50天的宝宝边吃奶边哭怎么回事| 亚洲熟妇熟女久久| www日本在线高清视频| a在线观看视频网站| 亚洲三区欧美一区| 亚洲成a人片在线一区二区| 无遮挡黄片免费观看| 91大片在线观看| 久久中文看片网| 精品国产国语对白av| 99精国产麻豆久久婷婷| 国产视频一区二区在线看| 国产国语露脸激情在线看| 丁香六月天网| 久久热在线av| 日本黄色日本黄色录像| 美女福利国产在线| 最近最新免费中文字幕在线| 动漫黄色视频在线观看| 99精品久久久久人妻精品| 久久久久久久国产电影| 国产男女内射视频| 宅男免费午夜| 啦啦啦在线免费观看视频4| 啦啦啦 在线观看视频| 亚洲av成人一区二区三| 国产伦人伦偷精品视频| 女性被躁到高潮视频| 一区在线观看完整版| 黑人巨大精品欧美一区二区蜜桃| 久久久久久亚洲精品国产蜜桃av| 国产伦理片在线播放av一区| 国产1区2区3区精品| 俄罗斯特黄特色一大片| 激情视频va一区二区三区| 中文字幕人妻熟女乱码| 老熟女久久久| 日本vs欧美在线观看视频| 他把我摸到了高潮在线观看 | 一本一本久久a久久精品综合妖精| 最近最新中文字幕大全电影3 | 国产一区二区三区视频了| 亚洲人成电影免费在线| 午夜福利乱码中文字幕| 国产日韩一区二区三区精品不卡| 欧美激情极品国产一区二区三区| 亚洲精品久久午夜乱码| 日本a在线网址| 少妇 在线观看| 青草久久国产| 国产精品久久久av美女十八| √禁漫天堂资源中文www| 日本欧美视频一区| 亚洲精品国产色婷婷电影| 亚洲av日韩精品久久久久久密| 午夜福利视频精品| 亚洲,欧美精品.| 国产免费福利视频在线观看| 国产97色在线日韩免费| 欧美日韩亚洲高清精品| 国产精品偷伦视频观看了| 色老头精品视频在线观看| 在线天堂中文资源库| 少妇裸体淫交视频免费看高清 | 无限看片的www在线观看| 欧美精品高潮呻吟av久久| 亚洲国产av新网站| 国产高清视频在线播放一区| 中文字幕人妻熟女乱码| 免费在线观看影片大全网站| 精品一区二区三区四区五区乱码| 好男人电影高清在线观看| 亚洲色图 男人天堂 中文字幕| 天堂俺去俺来也www色官网| av超薄肉色丝袜交足视频| 在线十欧美十亚洲十日本专区| 欧美激情高清一区二区三区| a级毛片黄视频| 18禁美女被吸乳视频| 国产成人精品久久二区二区免费| 国产亚洲精品久久久久5区| 免费不卡黄色视频| 脱女人内裤的视频| 757午夜福利合集在线观看| 电影成人av| 热re99久久国产66热| 日本精品一区二区三区蜜桃| 国产99久久九九免费精品| 91精品三级在线观看| 国产97色在线日韩免费| 免费在线观看影片大全网站| 人人妻人人澡人人看| www.自偷自拍.com| 国产成人一区二区三区免费视频网站| 国产成人欧美| 国产极品粉嫩免费观看在线| 国产精品99久久99久久久不卡| a级毛片在线看网站| 777久久人妻少妇嫩草av网站| 18禁美女被吸乳视频| 咕卡用的链子| 成年人免费黄色播放视频| 精品一区二区三区av网在线观看 | 99久久人妻综合| 亚洲国产av影院在线观看| 狠狠精品人妻久久久久久综合| 久久精品亚洲av国产电影网| 国产日韩欧美视频二区| 欧美日韩成人在线一区二区| 国产视频一区二区在线看| 黄片大片在线免费观看| 激情视频va一区二区三区| 久久精品亚洲av国产电影网| av线在线观看网站| 国产国语露脸激情在线看| www.自偷自拍.com|