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

    Analysis of dynamic features in intersecting pedestrian flows?

    2017-08-30 08:26:56HaiRongDong董海榮QiMeng孟琦XiuMingYao姚秀明XiaoXiaYang楊曉霞andQianLingWang王千齡
    Chinese Physics B 2017年9期

    Hai-Rong Dong(董海榮),Qi Meng(孟琦),Xiu-Ming Yao(姚秀明), Xiao-Xia Yang(楊曉霞),and Qian-Ling Wang(王千齡)

    1 State Key Laboratory of Rail Traffic Control and Safety,Beijing Jiaotong University,Beijing 100044,China

    2 School of Electronic and Information Engineering,Beijing Jiaotong University,Beijing 100044,China

    Analysis of dynamic features in intersecting pedestrian flows?

    Hai-Rong Dong(董海榮)1,Qi Meng(孟琦)1,Xiu-Ming Yao(姚秀明)2,?, Xiao-Xia Yang(楊曉霞)1,and Qian-Ling Wang(王千齡)1

    1 State Key Laboratory of Rail Traffic Control and Safety,Beijing Jiaotong University,Beijing 100044,China

    2 School of Electronic and Information Engineering,Beijing Jiaotong University,Beijing 100044,China

    This paper focuses on the simulation analysis of stripe formation and dynamic features of intersecting pedestrian flows. The intersecting flows consist of two streams of pedestrians and each pedestrian stream has a desired walking direction. The model adopted in the simulations is the social force model,which can reproduce the self-organization phenomena successfully.Three scenarios of different cross angles are established.The simulations confirm the empirical observations that there is a stripe formation when two streams of pedestrians intersect and the direction of the stripes is perpendicular to the sum of the directional vectors of the two streams.It can be concluded from the numerical simulation results that smaller cross angle results in higher mean speed and lower level of speed fluctuation.Moreover,the detailed pictures of pedestrians’moving behavior at intersections are given as well.

    intersecting pedestrian flows,social force model,stripes formation,cross angle

    1.Introduction

    Thousands of people died and more were injured in crowd stampedes over the past 100 years around the world.[1]Cities are densely populated nowadays and the collective behaviors seem to occur more frequently,such as traveling by subway,shopping in malls,and watching sports games in stadia. Hence,how to guarantee the safety and comfort of pedestrian motion and improve the efficiency of public space has become an urgent problem to be solved.

    Pedestrian dynamics refers to the study of how and where crowds form and move,which usually studies the dynamic features of pedestrian flows.The study on pedestrian dynamics has attracted considerable attention of researchers infields of physics,transportation,and more.Some good pedestrian models have been put forward with the advance of research.In general,pedestrian models are divided into three categories:(i)macroscopic models,where pedestrians are studied as a group without considering the interactions among pedestrians and pedestrian flows are analogous to gases and fluids in some sense;[2](ii)microscopic models,where every pedestrian is treated as an individual and the interactions among pedestrians are considered;[3](iii)mesoscopic models,which inherit from both macroscopic models and microscopic models.[4]The overview of some representative pedestrian models can be found in Refs.[5]–[7].Pedestrian models can describe the individual and group behaviors and reproduce the empirical scenarios,so they can help researchers to analyze the walking behaviors and dynamic features of pedestrian flows.In the previous investigations,arch-like blockings in bottlenecks,lane formation in bidirectional flows,and other self-organization phenomena have been observed and studied by various approaches.[1,8]Furthermore,the escape panic[9,10]and the emergency evacuation process under normal conditions and from a smoke-filled or afire room[11–15]have been studied by using multiple pedestrian models.

    Recently,intersecting flows become one of the topics of most concern because they are practically unavoidable and likely to cause clogging and accidents.The stripe formation is observed when two pedestrian flows intersect,but no stable pattern exists when three or more pedestrian flows intersect.[1,8]Zhang et al.proposed a two-dimensional model based on the cellular automata model in order to simulate pedestrians’crossing behaviors and the simulation results showed that this model can be applied successfully.[16]Saegusa et al.studied the intersecting flows under a periodic boundary condition by using the lattice gas model,and found that the overshooting in the fundamental diagram is attributed to the unstable congestion.[17]Guo et al.put forward a new microscopic model which is similar to the discrete choice model,and conducted several real pedestrian experiments with different intersecting angles to calibrate it.[18]Besides,pedestrians’avoidance behaviors in crossing flow were investigated in Ref.[19].

    On the basis of the previous studies,little research has been reported about the dynamic features of intersecting flows with different cross angles,so we aim to find the direction of the stripes,the influence of the cross angles on dynamic fea-tures,and the moving details of pedestrians at intersections in this paper.The model we choose is the social force model, which is a continuous microscopic model.It is universally accepted that self-organization phenomena are the consequence of interactions,which are considered in the social force model. By using the social force model,we can reproduce the stripes formation successfully.To reflect the effect of the cross angles, we set up several simulation scenarios in which only the cross angles are different,and then investigate the mean speeds and the fluctuations of speed.By comparing the mean speeds and the standard deviations,we can obtain whether the cross angles affect the motion of pedestrians and what the effects are. As for the moving details,we adopt the approach that picks some individuals and then snapshots them for each second in a time period when they are at intersections,then we find the common regularities.

    The organization of this paper is as follows.In Section 2, we will introduce the social force model.Section 3 presents the simulation scenarios.Section 4 shows the simulation results and their analysis.Section 5 is the conclusion summarizing the whole paper.

    2.The social force model

    Based on socio-psychological studies that both sociopsychology and physical forces affect the motion of pedestrians,Helbing et al.introduced the social force model to investigate the pedestrian movement,which is particularly suited for describing the occurrence of pressure in crowds.[9]

    The social force model considers the influence of the neighboring pedestrians and the environment on individuals’movement in crowds.Pedestrians are driven by:the desired force,the interaction forcefijbetween pedestriansiand j, the interaction forcefiWbetween pedestrianiand wall W.The mathematical formula of the social force model is expressed as

    where miis the mass of pedestrian i,and viis his or her actual velocity at time t.

    Generally,pedestrians have certain destinations,so each pedestrian has a desired velocity and a desired direction.In practice,the actual speed and the desired speed are different due to the interference of obstacles and other neighboring pedestrians.The desired force expresses pedestrians’willingness to achieve the desired velocity.Pedestrianilikes to move with a certain desired speedin a certain directionand therefore tends to adapt his or her actual velocity viwith a certain characteristic time τi.So the desired force is given by

    Pedestrian i prefers to keep a distance from other pedestrian j,so there exists an interaction forcefij.This force describes not only the psychological tendency to stay away from others,but also the physical force that occurs when the pedestrians touch each other.The interaction force between pedes-triansiandj is given by

    where ridenotes the position of the center of mass of pedestrian i.In two-dimensional space,ri=(xi,yi),xiand yiare respectively the horizontal and vertical coordinates of the center of mass of pedestrian i.denotes the distance between the centers of the two pedestrians.rijdenotes the sum of the radii of the two pedestrians,namely,rij=ri+rj,where riand rjare the radii of pedestriansi and j,respectively.nijrepresents the unit vector pointing from pedestrian j to pedestrian i,tijrepresents the tangential direction,

    g(x)is a piecewise function

    The term Aiexp[(rij?dij)/Bi]nijreflects the repulsive force describing the tendency that each pedestrian would like to keep a distance from others to avoid collisions,where Aiand Biare constants.

    When dijis smaller than rij,the pedestrians touch each other and compressions emerge.In this situation,there exist two extra forces.The term kg(rij?dij)nijreflects the body compression force,which can neutralize the compression of the bodies.The term κg(rij?dij)Δvtjitijreflects the sliding friction force,which blocks the tangential motion.k and κ are large positive constants.

    The interaction force between pedestrianiand wall W is treated analogously as

    where diWdenotes the distance between pedestriani and wall W,niWdenotes the direction perpendicular to the wall pointing from the wall to pedestrian i,and tiWreflects the tangential direction.Since the speed of the wall is zero,the tangential velocity can be expressed as

    3.Scenario setup

    The subway stations,malls,stadia,or other public places are likely to show intersecting pedestrian flows.Intersections of pedestrian flows are one of the major problems that affect the safety and efficiency of pedestrian movement and impossible to avoid in practice,so it is meaningful to know how people behave at intersections and the effect of the cross angles.

    Real-life experiment is unexpected and dangerous,so there is an urgent need to setup a reliable model to simulate the pedestrian motion.The parameters in the social force model are specified in Table 1,which is referred to in Ref.[9].

    Table 1.Parameters of the social force model.

    We study the intersecting flows which consist of two components of pedestrians moving into an area whose size is 50 m×50 m.The schematic diagram of the simulation scenarios is shown in Fig.1.Three simulation scenarios are built with different cross angles,which are 60°,90°,and 120°, named scenario 1,scenario 2,and scenario 3,respectively.For each scenario,the flow consists of two components of pedestrians;one component moves from west to east,and the other component moves from south to northeast or north or northwest.We use red hollow circles to represent the pedestrians going eastward and blue hollow circles to represent the pedestrians who walk in the other direction.

    In each scenario,the initial speed of pedestrians on both sides is 0.On each side,the distance between pedestrians’initial positions and the entrance is 1 m,which is designed in order to allow the pedestrians to enter the area with their normal walking state.The width of pedestrians’initial position line is 4 m,that is,the initial positions are(?1 m,13–17 m), (13–17 m,?1 m)in scenario 1,(?1 m,23–27 m),(23–27 m,?1 m)in scenario 2,and(?1 m,13–17 m),(23–27 m,?1 m) in scenario 3.Also,the initial number of pedestrians and the inflow rate of the two streams are exactly the same in the three scenarios.The initial number is 50,and the inflow rate is set to 5 p/s.The positions and speeds of pedestrians are updated after each time step,and time step ΔT is set to 0.01 s.

    Fig.1.(color online)Schematic diagram of the simulation scenarios. Hollow circles of different colors represent pedestrians walking in different directions,which are e1 and e2.θ represents the cross angle. The dashed square in the center of the area is the statistical area.

    The parameters of the social force model in these scenarios,namely,the desired speed,reaction time,and body size, are identical.According to Ref.[20],the parameters are set as follows.Reaction time τ satisfies a uniform distribution between 0.2 s and 0.5 s with interval Δτ=0.05 s.The desired speed is subject to a normal distribution,the mean of which is 1.35 m/s and the standard deviation is 0.2.The radius uniformly distributes between 0.2 m and 0.3 m with interval Δr=0.025 m.The parameters are presented in Table 2.

    Table 2.Parameter settings in different scenarios.

    The simulation time of each scenario is 150 s,and each scenario runs 10 times.

    To ensure the accurate analysis of dynamic characteristics at intersections,only the data collected when the two streams of pedestrians intersect completely are used,so we define a statistical area of size 10 m×10 m,which is the central dashed square shown in Fig.1.The selection of the location of the statistical area is essential.The coordinates of the statistical area’s vertexes are selected as follows:In scenario 1,S1(23 m, 22 m),S2(33 m,22 m),S3(33 m,12 m),S4(23 m,12 m);in scenario 2,S1(20 m,30 m),S2(30 m,30 m),S3(30 m,20 m), S4(20 m,20 m);in scenario 3,S1(11 m,20 m),S2(21 m, 20 m),S3(21 m,10 m),S4(11 m,10 m).Only the data in the statistical area are considered.The dynamic characteristics of intersecting pedestrian flows can be analyzed by the mean speed and the fluctuation of speed.According to the observations,stripes are formed after 20 s,so we calculate the mean speed in the statistical area from 40 s to 150 s when the stripes are basically stable.

    4.Analysis of the intersecting pedestrian flows

    4.1.Direction of stripes

    When two wide pedestrian flows intersect,stripe formation can be observed.This can be understood as a selforganized behavior,which means that the patterns are not planned or organized artificially,instead,the self-organization phenomena are considered as the consequence of interactions among pedestrians.What is more,the empirical investigations showed that the stripes are perpendicular to the sum of the directional vectors of the two pedestrian flows,[1]for example, the stripes are expected to have an angle of 45°in the case of a perpendicular intersection.

    We try to find the collective pattern of intersecting flows in our simulations.As shown in Fig.2,in order to find the regularity noticeably,the vector diagram is given in the northwest.The red vector is the directional vector of pedestrians moving from west to east,and the blue vector is the directional vector of pedestrians moving from south to northeast. It should be noted that the magnitudes of the two vectors are the same.The black vector is the sum direction.We draw a magenta line which is perpendicular to the black vector.It can be seen clearly that the direction of the stripes is the same as the magenta line.

    Figures 3 and 4 are analyzed analogously.It is apparent that the direction of the stripes is the same as the magenta line, which is perpendicular to the sum direction.

    The conclusion drawn from our simulations is identical to the empirical investigations,that is,the stripes are perpendicular to the sum of the directional vectors of the two pedestrian flows.

    Fig.2.(color online)Stripes formation in scenario 1.The vector diagram in the northwest is a contrast from the perspective of geometry. The coordinates of the statistical area’svertexes are(23 m,22 m),(33 m, 22 m),(33 m,12 m),and(23 m,12 m).

    Fig.3.(color online)Stripes formation in scenario 2.The vector diagram in the northwest is a contrast from the perspective of geometry. The coordinates of the statistical area’svertexes are(20 m,30 m),(30 m, 30 m),(30 m,20 m),and(20 m,20 m).

    Fig.4.(color online)Stripes formation in scenario 3.The vector diagram in the northeast is a contrast from the perspective of geometry. The coordinates of the statistical area’s vertexes are(11 m,20 m),(21 m, 20 m),(21 m,10 m),and(11 m,10 m).

    4.2.Analysis of speed with different angles

    The performance of pedestrian flows is characterized by the mean speed and the fluctuation of speed.When the simulation is conducted once,we obtain a mean speed which is calculated from 40 s to 150 s.After running simulations of each scenario for 10 times,we obtain 10 mean speeds for each scenario,the results are shown in Table 3.In order to show the difference intuitively,we draw a line chart(see Fig.5)with the data from Table 3.Then we calculate the average value of 10 mean speeds and obtain the corresponding standard deviation for each scenario(see Table 4).The standard deviations are all less than 0.01,which means that the results are close to each other,thus the stability and accuracy of the results can be confirmed.

    Table 3.Mean speeds(in units of m/s)of 10 simulations.

    Table 4.Average value and standard deviation of 10 mean speeds.

    Fig.5.(color online)Line chart of mean speeds of 10 simulations.The data are from Table 3.

    As for the fluctuation of speed,we choose one from 10 simulations of each scenario as the example.The fluctuation diagrams are shown in Fig.6,which reflect the change of the speeds with the elapse of time.It is worth noting that each speed point in the diagram is the average value of all pedestrians’speeds in the statistical area at each time step.The mean speeds and standard deviations of instant speeds are shown in Table 5.

    Table 5.One representation from 10 simulations.

    By analyzing Tables 4 and 5 and Figs.5 and 6,we canfind that greater cross angle results in lower mean speed and higher level of fluctuation.Since three simulation areas are of the same size,higher mean speed re flects shorter time spent to pass through this area.

    These results indicate that greater cross angle leads to stronger interactions among pedestrians in two directions, which means more con flicts at intersections.In order to get a precise and direct comprehension of the results,we draw a sample force diagram in which the cross angle is 120°(see Fig.7).Pedestrians who go eastward are regarded as group m and pedestrians moving in the other desired direction are regarded as group n.denotes the desired direction of pedestrian group m,is the perpendicular vector to it.fmnis the interaction force exerted on pedestrian group m by pedestrian group n.The interaction force can be decomposed in the directionsandThe component force in the directionopposes pedestrians’motion,while the component force in the directiondeviates the pedestrians’direction of motion.In scenario 1,on one hand,the interaction force exerted on pedestrian group m deflects the moving direction,on the other hand,it increases the speed of motion.In scenario 2, the interaction force deviates pedestrians’direction of motion. So the conclusion that greater cross angle leads to low mean speed and higher level of fluctuation is related to the interactions between pedestrians in two directions.When all pedestrians move in the same desired direction,that is,unidirectional pedestrian flows,pedestrians walk forward with rare interactions.In bidirectional pedestrian flows,pedestrians walking in opposite desired directions face each other with a great deal of conflicts,so the encountering pedestrians move a little aside to pass each other and the lane formation arises.Pedestrians moving in uniform lanes will have weak interactions,which is regarded as optimal self-organization.[21]Greater cross angle means bigger deviation between two directions,therefore pedestrians need to overcome more resistance to reach the desired location.

    Fig.7.(color online)Sample force diagram of scenario 3.Pedestrians who go eastward are regarded as group m and pedestrians moving in the other desired direction are regarded as group n.denotes the desired direction of pedestrian group m,is the perpendicular vector to it.fmn is the interaction force exerted on pedestrian group m by pedestrian group n.The interaction force can be decomposed in the two directions.

    Moreover,because of the dramatic speed changes in scenario 3,we infer preliminarily from experience that the fluctuation of speed is related to the process of congestion to some extent and a higher level of fluctuation is dangerous when the pedestrian density is large,while this inference should be verified by future work.

    4.3.Analysis of individual motion details

    In order to know how pedestrians move at intersections, we use different colors to mark several individual pedestrians, then we choose a time period and take snapshots for each second,the pictures are shown in Figs.8–10.

    Figure 8 exhibits the snapshots of scenario 1 from 48 s to 51 s.We use magenta,red,and orange solid circles to mark three pedestrians who move from west to east and green,blue, and black solid circles to mark three pedestrians who moves from south to northeast.The blue solid circle is at a stripe precisely.In order to show clearly,we draw a magenta ellipse to highlight the stripe.It can be observed that this pedestrian moves with his or her neighboring pedestrians who have the same desired direction.Figure 8(d)shows that the black solid circle approaches to the stripe continually,because they have the same direction.

    Fig.8.(color online)Snapshots of scenario 1 from 48 s to 51 s.Magenta,red,and orange solid circles represent three marked pedestrians who go eastward;green,blue,and black solid circles represent three marked pedestrians moving in the other direction. Magenta ellipse is used to show the stripe clearly.

    Fig.9.(color online)Snapshots of scenario 2 from 47 s to 50 s.The marks are the same as those in Fig.8.Magenta and black ellipses are used to show the stripes clearly.

    Figure 9 shows the snapshots of scenario 2 from 47 s to 50 s.Analogously,we use magenta,red,and orange solid circles to mark three pedestrians who move eastward and green, blue,and black solid circles to mark three pedestrians who move northward.At the beginning,the magenta solid circle is at a stripe,and we draw a magenta ellipse to highlight this stripe.With the elapse of time,the green solid circle approaches the intersection,and we use a black ellipse to highlight the stripe.It is obvious that they move with a group of pedestrians who have the same desired direction.What is more,the length of the stripe marked by the magenta ellipse has changed,which means that the pedestrians can extend along the stripe.

    Figure 10 shows the snapshots of scenario 3 from 46 s to 49 s.The marks are made in the same way as scenario 1.In Fig.10(a),the green solid one is at a stripe,we use a magenta rectangle to highlight it.The magenta and red solid circles are at the same stripe and we draw a black rectangle to highlight the stripe.It should be noted that the width of the stripes is narrow,so we use rectangles to highlight the stripes in pictures. We can see clearly that the pedestrians move forward with the stripes.

    The above observations and analysis lead us to the following conclusions.(i)Pedestrians always move forward with neighboring pedestrians who have the same desired direction, and these pedestrians form a group in the form of a stripe.(ii) As time goes by,it is possible that some other pedestrians who are in the same direction with the stripe and not far from the stripe will approach the stripe gradually(see Fig.8,the black one approaches the marked stripe gradually),and then join the pedestrian groups.(iii)The length of the stripe is not constant (see Fig.9),pedestrians can extend along the stripe.

    As we all know,pedestrians moving in two different directions occupy the same space when these pedestrians meet at intersections,so they are not able to choose the direct way to go forward.By forming stripes,they are allowed to penetrate each other with less conflicts and can continue walking to their destinations without stopping.In this way,the stripe formation is an optimal way to pass the intersection.

    Fig.10.(color online)Snapshots of scenario 3 from 46 s to 49 s.The marks are the same as those in Fig.8.Magenta and black rectangles are used to show the stripes clearly.

    5.Conclusion

    By using the social force model,we have investigated the direction of self-organized stripes,presented and analyzed the dynamic features of the intersecting flows with three typical cross angles,and found the way that pedestrians move at intersections.We have confirmed the empirical observations in simulations that pedestrians form stripes in a self-organized way when two pedestrian streams intersect and the direction of the stripes is perpendicular to the sum of the directional vectors of the two pedestrian flows.The simulation results show that greater cross angle results in more conflicts which are characterized by low mean speed and higher level of fluctuation.Greater cross angle is more likely to cause blockings and have adverse effects on the passing of pedestrians.What is more,we found that pedestrians move with stripes and can extend along the stripes.This phenomenon is regarded as an optimal way to pass the intersections quickly.

    [1]Helbing D,Buzna L,Johansson A and Werner T 2005 Transportation Science 39 1

    [2]Helbing D and Molnar P 1995 Phys.Rev.E 51 4282

    [3]Teknomo K,Takeyama Y and Inamura H 2000 Proceedings Japan Society of Civil Engineering Conference,March,2000,Morioka,Japan, p.1

    [4]Mehran R,Oyama Aand Shah M 2009 IEEE Computer Society Conference on Computer Vision and Pattern Recognition,June 20-26,2009, Miami,USA,p.935

    [5]Zheng X P,Zhong T K and Liu M T 2009 Building and Environment 44 437

    [6]Schadschneider A,Klüpfel H,Kretz T,Rogsch C and Seyfried A 2009 Multi-AgentSystems for Traffic and Transportation Engineering(Pennsylvania:IGI Global)pp.125–154

    [7]Duives D C,Daamen W and Hoogendoorn S 2013 Transportation Research Part C 37 193

    [8]Hoogendoorn S and Bovy P H L 2003 Optimal Control Applications and Methods 24 153

    [9]Helbing D,Farkas I and Vicsek T 2000 Nature 407 487

    [10]Zou Y,Xie J R and Wang B H 2016 Plos One 11 e0153388

    [11]Yuan W F and Tan K H 2011 Physica A 390 4210

    [12]Shi J Y,Ren A Z and Chen C 2009 Automation in Construction 18 338

    [13]Zheng Y,Jia B,Li X G and Zhu N 2011 Physica A 390 3147

    [14]Wang L,Zheng J H,Zhang X S,Zhang J L Wang Q Z and Zhang Q 2016 Chin.Phys.B 11 118901

    [15]Yang X X,Dong H R,Yao X M and Sun X B 2016 Chin.Phys.B 25 048902

    [16]Zhang J,Wang H and Li P 2004 Journal of Zhejiang University Science 5 835

    [17]Saegusa T,Mashiko T and Nagatani T 2008 Physica A 387 4119

    [18]Guo R Y,Wong S C and Huang H J 2010 Physica A 389 515

    [19]Das S,Manski C F and Manuszak M D 2005 Journal of Applied Econometrics 20 529

    [20]Yang X X,Winnie D,Serge P H,Dong H R and Yao X M 2016 Chin. Phys.B 25 028901

    [21]Helbing D and Vicsek T 1999 New J.Phys.1 13

    25 February 2017;revised manuscript

    22 May 2017;published online 18 July 2017)

    10.1088/1674-1056/26/9/098902

    ?Project supported by the National Natural Science Foundation of China(Grant No.61233001)and the Fundamental Research Funds for the Central Universities, China(Grant No.2017JBM014).

    ?Corresponding author.E-mail:xmyao@bjtu.edu.cn

    ?2017 Chinese Physical Society and IOP Publishing Ltd http://iopscience.iop.org/cpb http://cpb.iphy.ac.cn

    好看av亚洲va欧美ⅴa在| 国产精品美女特级片免费视频播放器| 国产单亲对白刺激| 国产探花极品一区二区| 人人妻,人人澡人人爽秒播| 哪里可以看免费的av片| 国产亚洲精品久久久com| 日本三级黄在线观看| 欧美高清成人免费视频www| 全区人妻精品视频| 久久精品国产清高在天天线| 欧美激情在线99| 美女 人体艺术 gogo| 免费在线观看日本一区| 精品人妻偷拍中文字幕| 九九热线精品视视频播放| 国产欧美日韩一区二区精品| 丁香六月欧美| 老熟妇仑乱视频hdxx| 搡女人真爽免费视频火全软件 | 欧美一级a爱片免费观看看| 久久国产精品影院| 久久中文看片网| 女人十人毛片免费观看3o分钟| 欧美中文日本在线观看视频| 亚洲欧美日韩东京热| 日本五十路高清| a级毛片免费高清观看在线播放| 麻豆国产av国片精品| 国产又黄又爽又无遮挡在线| 国产极品精品免费视频能看的| 免费av毛片视频| 毛片一级片免费看久久久久 | 午夜激情欧美在线| 日日摸夜夜添夜夜添av毛片 | 最近视频中文字幕2019在线8| 亚洲国产色片| 午夜免费男女啪啪视频观看 | 永久网站在线| 成年女人永久免费观看视频| 美女 人体艺术 gogo| 麻豆一二三区av精品| 国内精品美女久久久久久| 精品乱码久久久久久99久播| 日韩亚洲欧美综合| 真人做人爱边吃奶动态| 亚洲真实伦在线观看| 午夜福利免费观看在线| 男人狂女人下面高潮的视频| 老熟妇乱子伦视频在线观看| 亚洲激情在线av| 午夜福利在线观看吧| 亚洲不卡免费看| 白带黄色成豆腐渣| 宅男免费午夜| 高清毛片免费观看视频网站| 久久久久久久亚洲中文字幕 | 两人在一起打扑克的视频| 中文字幕高清在线视频| 国产av在哪里看| 国产精品自产拍在线观看55亚洲| 成年女人毛片免费观看观看9| 欧美激情国产日韩精品一区| 在线观看一区二区三区| 69av精品久久久久久| 成人亚洲精品av一区二区| 国产精品乱码一区二三区的特点| 久久婷婷人人爽人人干人人爱| 69av精品久久久久久| 国产成人啪精品午夜网站| 白带黄色成豆腐渣| 国产精品,欧美在线| 中国美女看黄片| 久久久久精品国产欧美久久久| .国产精品久久| 丰满的人妻完整版| 俺也久久电影网| 日本一本二区三区精品| 亚洲 国产 在线| a级毛片a级免费在线| 国产中年淑女户外野战色| 少妇的逼水好多| 精品人妻1区二区| 美女xxoo啪啪120秒动态图 | 国产探花极品一区二区| 精品国产三级普通话版| 一本久久中文字幕| 小蜜桃在线观看免费完整版高清| 亚洲精品一区av在线观看| 免费看美女性在线毛片视频| 国产精品永久免费网站| 亚洲成人免费电影在线观看| 波多野结衣高清作品| 久99久视频精品免费| 小蜜桃在线观看免费完整版高清| 少妇人妻一区二区三区视频| 69人妻影院| 美女cb高潮喷水在线观看| 精品人妻一区二区三区麻豆 | 国产精品嫩草影院av在线观看 | 欧美日韩瑟瑟在线播放| 在线免费观看的www视频| 久久午夜亚洲精品久久| 91在线观看av| 亚洲成人中文字幕在线播放| 精品久久久久久久末码| 国内精品久久久久精免费| 久久久久久九九精品二区国产| 丰满乱子伦码专区| 中文字幕免费在线视频6| 亚洲精品一卡2卡三卡4卡5卡| 九九久久精品国产亚洲av麻豆| 午夜老司机福利剧场| 高清在线国产一区| 91字幕亚洲| 俄罗斯特黄特色一大片| 直男gayav资源| 久久精品91蜜桃| 亚洲熟妇熟女久久| 成年免费大片在线观看| 亚洲专区国产一区二区| 成人特级av手机在线观看| 桃色一区二区三区在线观看| av福利片在线观看| 欧美性猛交黑人性爽| 欧美黄色淫秽网站| 亚洲 欧美 日韩 在线 免费| 亚洲性夜色夜夜综合| 国产v大片淫在线免费观看| 欧美不卡视频在线免费观看| 97超视频在线观看视频| 午夜激情福利司机影院| 国语自产精品视频在线第100页| 国产精品av视频在线免费观看| 天美传媒精品一区二区| 夜夜爽天天搞| www.999成人在线观看| 最好的美女福利视频网| www.999成人在线观看| 国产真实乱freesex| 午夜a级毛片| 日本 欧美在线| 十八禁国产超污无遮挡网站| 精品人妻1区二区| 欧美在线一区亚洲| 窝窝影院91人妻| 欧美日本亚洲视频在线播放| 欧美黑人巨大hd| 日本一本二区三区精品| 亚洲 国产 在线| 69人妻影院| 日韩欧美国产在线观看| av天堂在线播放| 久久6这里有精品| 亚洲天堂国产精品一区在线| 好男人电影高清在线观看| 亚洲av一区综合| 欧美国产日韩亚洲一区| 久久欧美精品欧美久久欧美| 91在线精品国自产拍蜜月| 在线播放国产精品三级| 身体一侧抽搐| 偷拍熟女少妇极品色| 日本三级黄在线观看| 亚洲最大成人中文| 麻豆国产97在线/欧美| 性欧美人与动物交配| 亚洲专区国产一区二区| 俄罗斯特黄特色一大片| av黄色大香蕉| 精品久久久久久,| 一区二区三区免费毛片| 在线观看舔阴道视频| 热99在线观看视频| 日本撒尿小便嘘嘘汇集6| 亚洲经典国产精华液单 | 99久国产av精品| 亚洲人成网站在线播放欧美日韩| 国产精品女同一区二区软件 | 亚洲精品一卡2卡三卡4卡5卡| АⅤ资源中文在线天堂| 99热这里只有精品一区| 午夜福利在线在线| 久久九九热精品免费| 9191精品国产免费久久| 久久精品国产自在天天线| 亚洲成av人片免费观看| 亚洲国产精品合色在线| 少妇的逼好多水| 如何舔出高潮| 91久久精品电影网| 日本 欧美在线| 麻豆一二三区av精品| 在线观看一区二区三区| 久9热在线精品视频| 日日摸夜夜添夜夜添av毛片 | 男插女下体视频免费在线播放| 免费av观看视频| 久久这里只有精品中国| 国产69精品久久久久777片| 欧美日韩综合久久久久久 | 午夜福利欧美成人| 欧美激情在线99| 日韩欧美国产在线观看| 永久网站在线| 亚州av有码| 一本久久中文字幕| 亚洲第一电影网av| 在线观看66精品国产| 午夜a级毛片| 日韩av在线大香蕉| 欧美另类亚洲清纯唯美| 国内精品一区二区在线观看| 午夜福利18| 免费观看人在逋| 一级av片app| 国产免费男女视频| 欧美日本亚洲视频在线播放| 美女xxoo啪啪120秒动态图 | 午夜福利18| 69人妻影院| 日本三级黄在线观看| 欧美激情久久久久久爽电影| 国内揄拍国产精品人妻在线| 欧美xxxx性猛交bbbb| 欧美精品啪啪一区二区三区| 亚洲国产欧美人成| 久久伊人香网站| 波多野结衣高清无吗| 国模一区二区三区四区视频| 久久精品国产亚洲av涩爱 | 搞女人的毛片| 久久久久久大精品| 丰满人妻熟妇乱又伦精品不卡| 99国产极品粉嫩在线观看| 色精品久久人妻99蜜桃| 亚洲最大成人中文| 欧美一区二区亚洲| 天堂av国产一区二区熟女人妻| 午夜a级毛片| 一本一本综合久久| 韩国av一区二区三区四区| 3wmmmm亚洲av在线观看| 69av精品久久久久久| 亚洲精品乱码久久久v下载方式| 久久香蕉精品热| 精品午夜福利在线看| 亚洲成av人片免费观看| 久久人妻av系列| 动漫黄色视频在线观看| 长腿黑丝高跟| 国产成人aa在线观看| 中文资源天堂在线| 日本五十路高清| 欧美日韩乱码在线| 欧美极品一区二区三区四区| 免费看光身美女| 免费黄网站久久成人精品 | 在线观看美女被高潮喷水网站 | 夜夜看夜夜爽夜夜摸| 久久精品综合一区二区三区| 成人鲁丝片一二三区免费| 麻豆久久精品国产亚洲av| 久久亚洲真实| 国产毛片a区久久久久| 国产亚洲欧美在线一区二区| 99热这里只有是精品50| 国产主播在线观看一区二区| 91麻豆av在线| 精品久久久久久,| 亚洲人成电影免费在线| 精品一区二区三区视频在线观看免费| 亚洲aⅴ乱码一区二区在线播放| 99riav亚洲国产免费| 尤物成人国产欧美一区二区三区| 色噜噜av男人的天堂激情| 中亚洲国语对白在线视频| 久久人妻av系列| 97人妻精品一区二区三区麻豆| 国产真实乱freesex| 女同久久另类99精品国产91| 日日干狠狠操夜夜爽| 日韩欧美在线乱码| 国产亚洲欧美在线一区二区| 国产色婷婷99| 亚洲av免费在线观看| 亚洲av不卡在线观看| 国产乱人视频| 成人av在线播放网站| 网址你懂的国产日韩在线| 欧美性感艳星| 日日摸夜夜添夜夜添av毛片 | 精品一区二区三区视频在线| 不卡一级毛片| 久久九九热精品免费| 国产精品1区2区在线观看.| 国产精品爽爽va在线观看网站| 久久久精品大字幕| 国产精品免费一区二区三区在线| 丰满人妻熟妇乱又伦精品不卡| 精品久久久久久,| 精品人妻熟女av久视频| 1000部很黄的大片| 在线观看午夜福利视频| 又紧又爽又黄一区二区| 嫩草影院入口| 一个人看的www免费观看视频| 深爱激情五月婷婷| 国产av麻豆久久久久久久| 欧美黄色片欧美黄色片| 国产色爽女视频免费观看| 久久精品国产亚洲av香蕉五月| 国产精品久久久久久久久免 | 久久午夜亚洲精品久久| 国产一区二区三区视频了| 丝袜美腿在线中文| 日韩国内少妇激情av| 久久久久久久午夜电影| 99精品在免费线老司机午夜| 亚洲熟妇中文字幕五十中出| 国产 一区 欧美 日韩| av国产免费在线观看| 美女 人体艺术 gogo| 国产一级毛片七仙女欲春2| 男女视频在线观看网站免费| 国产乱人伦免费视频| 免费看美女性在线毛片视频| 国产成+人综合+亚洲专区| 九九热线精品视视频播放| 国产欧美日韩精品一区二区| 久久久久久国产a免费观看| 夜夜躁狠狠躁天天躁| 国产91精品成人一区二区三区| 亚洲不卡免费看| eeuss影院久久| 女人被狂操c到高潮| 午夜免费成人在线视频| 麻豆av噜噜一区二区三区| 男人舔奶头视频| 日韩 亚洲 欧美在线| 每晚都被弄得嗷嗷叫到高潮| 可以在线观看的亚洲视频| 日本黄大片高清| 一本久久中文字幕| 国产一级毛片七仙女欲春2| 亚洲五月婷婷丁香| 一区二区三区激情视频| 久久性视频一级片| 不卡一级毛片| 欧美日韩福利视频一区二区| 亚洲 国产 在线| 欧美日韩福利视频一区二区| 男人舔奶头视频| 三级毛片av免费| 久久热精品热| aaaaa片日本免费| 亚洲精品在线观看二区| 成人鲁丝片一二三区免费| 人人妻人人澡欧美一区二区| 性欧美人与动物交配| 亚洲美女黄片视频| 如何舔出高潮| 国产伦精品一区二区三区视频9| 精品午夜福利视频在线观看一区| 婷婷色综合大香蕉| 亚洲国产精品久久男人天堂| 香蕉av资源在线| 精品人妻一区二区三区麻豆 | 女同久久另类99精品国产91| 美女cb高潮喷水在线观看| 日本黄色视频三级网站网址| 精品一区二区三区视频在线| 成人国产一区最新在线观看| 免费在线观看成人毛片| 欧美日本亚洲视频在线播放| 亚洲欧美精品综合久久99| 搡老岳熟女国产| 欧美成人免费av一区二区三区| 亚洲国产精品sss在线观看| 亚洲一区高清亚洲精品| 美女免费视频网站| 亚洲人成网站在线播| 在线观看一区二区三区| 级片在线观看| 亚洲av成人av| 精品久久久久久久人妻蜜臀av| 九九热线精品视视频播放| 免费av观看视频| 听说在线观看完整版免费高清| 在线观看舔阴道视频| 国产精品女同一区二区软件 | 精品久久久久久久久亚洲 | 变态另类成人亚洲欧美熟女| 禁无遮挡网站| 中文字幕av成人在线电影| 国产中年淑女户外野战色| 亚洲真实伦在线观看| 桃色一区二区三区在线观看| 少妇被粗大猛烈的视频| 级片在线观看| 国产国拍精品亚洲av在线观看| 中文字幕久久专区| 淫秽高清视频在线观看| 欧美绝顶高潮抽搐喷水| 性插视频无遮挡在线免费观看| 国产精品1区2区在线观看.| www.www免费av| 久久精品91蜜桃| 99在线视频只有这里精品首页| 两个人的视频大全免费| 狂野欧美白嫩少妇大欣赏| 午夜免费成人在线视频| 精品日产1卡2卡| 男女做爰动态图高潮gif福利片| 久久伊人香网站| 免费电影在线观看免费观看| av专区在线播放| 一个人免费在线观看的高清视频| 两性午夜刺激爽爽歪歪视频在线观看| 色精品久久人妻99蜜桃| 亚洲一区二区三区色噜噜| 精品熟女少妇八av免费久了| 波多野结衣高清无吗| 性欧美人与动物交配| 男女下面进入的视频免费午夜| 又爽又黄a免费视频| 波多野结衣高清作品| 麻豆一二三区av精品| 欧美性感艳星| 免费无遮挡裸体视频| 日韩高清综合在线| 国产精品一区二区性色av| 久久午夜亚洲精品久久| 老熟妇仑乱视频hdxx| 欧美日韩福利视频一区二区| 亚洲av成人不卡在线观看播放网| 97超级碰碰碰精品色视频在线观看| ponron亚洲| 一级作爱视频免费观看| 1000部很黄的大片| 在线看三级毛片| 国产野战对白在线观看| 午夜影院日韩av| 可以在线观看的亚洲视频| 日本黄大片高清| 一级作爱视频免费观看| 一区二区三区高清视频在线| 99热6这里只有精品| 精品人妻1区二区| 亚洲成a人片在线一区二区| 亚洲乱码一区二区免费版| 久久精品国产99精品国产亚洲性色| 精品人妻视频免费看| 欧美中文日本在线观看视频| 久久精品久久久久久噜噜老黄 | 国产亚洲精品久久久com| 国产视频一区二区在线看| 一本一本综合久久| 精品国产三级普通话版| 亚洲欧美精品综合久久99| 日韩中字成人| 国内久久婷婷六月综合欲色啪| 黄色女人牲交| 久久国产精品人妻蜜桃| 十八禁网站免费在线| 日日摸夜夜添夜夜添小说| 欧美日韩中文字幕国产精品一区二区三区| 在线播放无遮挡| 欧美成人性av电影在线观看| av天堂中文字幕网| 欧美xxxx黑人xx丫x性爽| 免费av不卡在线播放| www.www免费av| 国产午夜福利久久久久久| 99在线视频只有这里精品首页| 久久久久久九九精品二区国产| 中文字幕精品亚洲无线码一区| 国产高清视频在线播放一区| 欧美色欧美亚洲另类二区| av中文乱码字幕在线| 国产单亲对白刺激| 少妇被粗大猛烈的视频| 又黄又爽又刺激的免费视频.| 久久久久久国产a免费观看| 国产亚洲精品av在线| 精品免费久久久久久久清纯| 小说图片视频综合网站| 午夜激情福利司机影院| 精品国产三级普通话版| 国产精品电影一区二区三区| 性色av乱码一区二区三区2| 亚洲人成网站在线播放欧美日韩| 非洲黑人性xxxx精品又粗又长| 少妇丰满av| 亚洲在线自拍视频| 在线观看美女被高潮喷水网站 | 两个人视频免费观看高清| 亚洲中文字幕日韩| 国产视频内射| 国产精品久久久久久精品电影| 99热精品在线国产| 成人午夜高清在线视频| 我的老师免费观看完整版| 亚洲一区高清亚洲精品| 精品乱码久久久久久99久播| 久久久久免费精品人妻一区二区| 午夜老司机福利剧场| 神马国产精品三级电影在线观看| 亚洲国产精品久久男人天堂| 国产精品日韩av在线免费观看| 欧美激情久久久久久爽电影| 久久久国产成人免费| 亚洲自拍偷在线| 观看免费一级毛片| 国产精品女同一区二区软件 | 深夜精品福利| 一本综合久久免费| 午夜福利视频1000在线观看| 自拍偷自拍亚洲精品老妇| 国产av一区在线观看免费| a级毛片免费高清观看在线播放| 白带黄色成豆腐渣| 欧美xxxx性猛交bbbb| 亚洲片人在线观看| 亚洲av中文字字幕乱码综合| 嫩草影院入口| 久久欧美精品欧美久久欧美| 尤物成人国产欧美一区二区三区| 我要看日韩黄色一级片| 在线观看66精品国产| 精品无人区乱码1区二区| 黄色女人牲交| 在线播放无遮挡| 最后的刺客免费高清国语| 精品人妻1区二区| 中文字幕人妻熟人妻熟丝袜美| 亚洲最大成人av| 又紧又爽又黄一区二区| 日本撒尿小便嘘嘘汇集6| 亚洲激情在线av| 国产熟女xx| 国产精华一区二区三区| 亚洲成人精品中文字幕电影| 亚洲在线观看片| 91在线精品国自产拍蜜月| 免费在线观看影片大全网站| 成人亚洲精品av一区二区| 国产一区二区在线av高清观看| 亚洲av五月六月丁香网| 日韩 亚洲 欧美在线| 久久99热这里只有精品18| 小蜜桃在线观看免费完整版高清| 欧美成人一区二区免费高清观看| 亚洲精品在线观看二区| 老女人水多毛片| 国产伦人伦偷精品视频| 色哟哟哟哟哟哟| 欧美日韩亚洲国产一区二区在线观看| 国产精品三级大全| 一个人看的www免费观看视频| 99热6这里只有精品| 色吧在线观看| 99国产精品一区二区三区| 国产淫片久久久久久久久 | 99久久精品国产亚洲精品| 我的女老师完整版在线观看| 国产免费男女视频| 成人性生交大片免费视频hd| 亚洲美女视频黄频| netflix在线观看网站| 国产精华一区二区三区| 日韩中字成人| 小蜜桃在线观看免费完整版高清| 久久久久国内视频| 亚洲国产精品sss在线观看| 国产精品电影一区二区三区| 色精品久久人妻99蜜桃| 午夜影院日韩av| 国产色婷婷99| 嫩草影院入口| 亚洲美女黄片视频| 怎么达到女性高潮| 欧美成人a在线观看| 自拍偷自拍亚洲精品老妇| 夜夜躁狠狠躁天天躁| 91久久精品国产一区二区成人| 中文字幕av成人在线电影| 天堂影院成人在线观看| 一夜夜www| 日韩免费av在线播放| 国内久久婷婷六月综合欲色啪| 18禁在线播放成人免费| 悠悠久久av| 黄色丝袜av网址大全| 一个人免费在线观看电影| 午夜福利高清视频| 天堂影院成人在线观看| 最近中文字幕高清免费大全6 | 亚洲美女黄片视频| 欧美一级a爱片免费观看看| 夜夜爽天天搞| 欧美又色又爽又黄视频| 乱码一卡2卡4卡精品| 搡老妇女老女人老熟妇| 亚洲 国产 在线| 又黄又爽又刺激的免费视频.| 狠狠狠狠99中文字幕| 一进一出抽搐gif免费好疼| 久久伊人香网站| 偷拍熟女少妇极品色| 成人毛片a级毛片在线播放| 国产精品永久免费网站| 免费看美女性在线毛片视频| 国内少妇人妻偷人精品xxx网站| 国产又黄又爽又无遮挡在线|