• <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| 五月天丁香电影| 最近的中文字幕免费完整| 韩国av在线不卡| 建设人人有责人人尽责人人享有的| 亚洲精品在线美女| 男男h啪啪无遮挡| 丝袜美足系列| 97精品久久久久久久久久精品| 国产日韩欧美视频二区| 男女午夜视频在线观看| 久久亚洲国产成人精品v| 国产免费现黄频在线看| 欧美精品人与动牲交sv欧美| 人人妻人人澡人人看| 免费不卡黄色视频| 久久ye,这里只有精品| 日日撸夜夜添| 亚洲精品久久午夜乱码| 亚洲一区中文字幕在线| 国产免费一区二区三区四区乱码| 亚洲国产精品成人久久小说| 老司机影院毛片| 久久精品aⅴ一区二区三区四区| 亚洲图色成人| 性色av一级| www日本在线高清视频| 日本av免费视频播放| 亚洲男人天堂网一区| 久久久久久免费高清国产稀缺| 一边摸一边做爽爽视频免费| 少妇精品久久久久久久| 99re6热这里在线精品视频| 久久鲁丝午夜福利片| 国产男女超爽视频在线观看| 最新的欧美精品一区二区| 日本午夜av视频| 操美女的视频在线观看| 久久综合国产亚洲精品| 又黄又粗又硬又大视频| 国产乱来视频区| 国产成人午夜福利电影在线观看| 大香蕉久久网| 亚洲国产精品999| 自线自在国产av| 狂野欧美激情性bbbbbb| 国产成人av激情在线播放| 亚洲国产中文字幕在线视频| 久久韩国三级中文字幕| 国产精品一二三区在线看| 国产精品麻豆人妻色哟哟久久| 深夜精品福利| 电影成人av| 亚洲欧美精品自产自拍| av在线观看视频网站免费| 男人操女人黄网站| 欧美日韩福利视频一区二区| 性色av一级| 午夜福利网站1000一区二区三区| 亚洲第一av免费看| 亚洲av综合色区一区| 在线观看人妻少妇| 丝瓜视频免费看黄片| 美女视频免费永久观看网站| 韩国精品一区二区三区| 亚洲男人天堂网一区| 美女脱内裤让男人舔精品视频| 久久久久久久久久久久大奶| 亚洲国产精品一区三区| 国产男女内射视频| 一区二区av电影网| 国产亚洲最大av| 99久久综合免费| 97人妻天天添夜夜摸| 在线亚洲精品国产二区图片欧美| 日韩视频在线欧美| 国产男女超爽视频在线观看| 欧美日韩成人在线一区二区| 多毛熟女@视频| 观看美女的网站| 无限看片的www在线观看| 国产乱来视频区| 国产一区有黄有色的免费视频| 最近最新中文字幕大全免费视频 | 热99久久久久精品小说推荐| 欧美日韩一区二区视频在线观看视频在线| 操美女的视频在线观看| 亚洲欧美一区二区三区国产| 亚洲欧美色中文字幕在线| 久久鲁丝午夜福利片| 99精品久久久久人妻精品| √禁漫天堂资源中文www| 精品久久久精品久久久| 黑人巨大精品欧美一区二区蜜桃| 国产毛片在线视频| 又大又爽又粗| 亚洲国产中文字幕在线视频| 亚洲精品日本国产第一区| 狠狠精品人妻久久久久久综合| 久久精品熟女亚洲av麻豆精品| videosex国产| bbb黄色大片| 久久ye,这里只有精品| 国产在视频线精品| 亚洲精品国产av蜜桃| 久久精品aⅴ一区二区三区四区| 大陆偷拍与自拍| 国产精品无大码| 久久久国产一区二区| 黄色毛片三级朝国网站| 国产一级毛片在线| 日本黄色日本黄色录像| 黄片无遮挡物在线观看| 久久热在线av| 最近的中文字幕免费完整| 国产一区有黄有色的免费视频| 午夜福利免费观看在线| 日本av免费视频播放| 一级片'在线观看视频| 亚洲av福利一区| 中文字幕高清在线视频| 三上悠亚av全集在线观看| 熟女av电影| 黄片无遮挡物在线观看| 成人影院久久| 国产亚洲av高清不卡| 国产精品麻豆人妻色哟哟久久| 久久精品人人爽人人爽视色| 一区二区av电影网| 亚洲色图综合在线观看| 99热全是精品| 波野结衣二区三区在线| 哪个播放器可以免费观看大片| 大片电影免费在线观看免费| 欧美乱码精品一区二区三区| 水蜜桃什么品种好| avwww免费| 在线天堂最新版资源| 精品国产一区二区久久| 少妇猛男粗大的猛烈进出视频| 亚洲四区av| 下体分泌物呈黄色| h视频一区二区三区| 国产免费福利视频在线观看| 欧美日韩一级在线毛片| 亚洲国产中文字幕在线视频| 日本av手机在线免费观看| 少妇猛男粗大的猛烈进出视频| 午夜福利影视在线免费观看| 蜜桃国产av成人99| 1024香蕉在线观看| 亚洲精品,欧美精品| 国产精品三级大全| 2021少妇久久久久久久久久久| 久久久久久人妻| 亚洲成av片中文字幕在线观看| 国产精品久久久久久人妻精品电影 | www日本在线高清视频| 免费黄频网站在线观看国产| 人妻一区二区av| 日韩一本色道免费dvd| 在线观看三级黄色| 亚洲欧美中文字幕日韩二区| 精品一区二区三卡| 精品福利永久在线观看| 又大又黄又爽视频免费| 一级片'在线观看视频| 中文字幕最新亚洲高清| 日日摸夜夜添夜夜爱| av卡一久久| 国产精品久久久久久精品古装| 欧美黑人精品巨大| 国产精品久久久久久人妻精品电影 | 国产无遮挡羞羞视频在线观看| 亚洲在久久综合| 麻豆av在线久日| 欧美久久黑人一区二区| 晚上一个人看的免费电影| 国产伦理片在线播放av一区| 制服诱惑二区| 毛片一级片免费看久久久久| 亚洲伊人色综图| 伦理电影大哥的女人| 亚洲熟女精品中文字幕| 男男h啪啪无遮挡| 久久久久久免费高清国产稀缺| 80岁老熟妇乱子伦牲交| 亚洲精品久久久久久婷婷小说| 少妇的丰满在线观看| 欧美另类一区| 亚洲av综合色区一区| 日本黄色日本黄色录像| 日韩一卡2卡3卡4卡2021年| 丝袜美足系列| 精品视频人人做人人爽| 国产成人精品无人区| 丰满迷人的少妇在线观看| 国产精品国产三级国产专区5o| 秋霞伦理黄片| 成年人免费黄色播放视频| 别揉我奶头~嗯~啊~动态视频 | 波野结衣二区三区在线| 亚洲国产av新网站| 一本大道久久a久久精品| 男人舔女人的私密视频| 欧美日韩亚洲综合一区二区三区_| 成人国语在线视频| 热99国产精品久久久久久7| 国产又爽黄色视频| 精品福利永久在线观看| 久久 成人 亚洲| 欧美亚洲日本最大视频资源| 国产精品香港三级国产av潘金莲 | 亚洲欧美成人综合另类久久久| 亚洲精品久久午夜乱码| 亚洲国产精品国产精品| 国产伦人伦偷精品视频| av天堂久久9| 精品国产超薄肉色丝袜足j| 亚洲第一青青草原| 黄频高清免费视频| 两性夫妻黄色片| 天堂8中文在线网| 超碰97精品在线观看| 欧美人与善性xxx| 黄片无遮挡物在线观看| 啦啦啦视频在线资源免费观看| 丝袜在线中文字幕| 国产精品av久久久久免费| 一级毛片 在线播放| 美国免费a级毛片| 久久久国产欧美日韩av| 99国产综合亚洲精品| 日韩大片免费观看网站| 国产精品嫩草影院av在线观看| 巨乳人妻的诱惑在线观看| 亚洲少妇的诱惑av| 久久精品国产亚洲av高清一级| 亚洲国产欧美在线一区| 777久久人妻少妇嫩草av网站| 国产精品.久久久| 一区二区av电影网| 9热在线视频观看99| 国产精品.久久久| 免费不卡黄色视频| av天堂久久9| av电影中文网址| 久久国产精品男人的天堂亚洲| 男的添女的下面高潮视频| 欧美精品av麻豆av| 在线亚洲精品国产二区图片欧美| 亚洲第一区二区三区不卡| 一本色道久久久久久精品综合| 99久国产av精品国产电影| 亚洲国产看品久久| 又大又爽又粗| 国产av国产精品国产| 久久久久久久大尺度免费视频| 国产欧美亚洲国产| 亚洲欧美成人综合另类久久久| 美国免费a级毛片| 成人午夜精彩视频在线观看| 日韩一本色道免费dvd| av天堂久久9| 亚洲av在线观看美女高潮| 久久天堂一区二区三区四区| a 毛片基地| 亚洲成人手机| av卡一久久| 黄色怎么调成土黄色| 国产成人a∨麻豆精品| 国产极品粉嫩免费观看在线| 久久精品久久久久久噜噜老黄| 少妇的丰满在线观看| 国产激情久久老熟女| 国产1区2区3区精品| www.熟女人妻精品国产| 精品一区二区三区av网在线观看 | 尾随美女入室| 精品国产乱码久久久久久男人| 亚洲美女黄色视频免费看| 亚洲av福利一区| 狂野欧美激情性bbbbbb| 欧美日本中文国产一区发布| 中文字幕色久视频| 2021少妇久久久久久久久久久| 在线观看免费视频网站a站| 最新的欧美精品一区二区| 亚洲一区中文字幕在线| 黄色一级大片看看| 国产男人的电影天堂91| 欧美精品亚洲一区二区| 尾随美女入室| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲情色 制服丝袜| 青春草国产在线视频| 熟女av电影| 街头女战士在线观看网站| 亚洲精品乱久久久久久| 一级毛片电影观看| 久久97久久精品| 我要看黄色一级片免费的| 亚洲自偷自拍图片 自拍| 99热全是精品| 国产乱人偷精品视频| 亚洲精品美女久久久久99蜜臀 | 热re99久久精品国产66热6| 大香蕉久久网| 免费在线观看黄色视频的| 国产精品人妻久久久影院| 97人妻天天添夜夜摸| 深夜精品福利| 成人国产av品久久久| 久久99精品国语久久久| 综合色丁香网| 欧美精品高潮呻吟av久久| 黄色视频不卡| 国产97色在线日韩免费| 免费黄色在线免费观看| 日本黄色日本黄色录像| 国产97色在线日韩免费| 一区二区三区激情视频| 男女之事视频高清在线观看 | 免费观看人在逋| 亚洲精品日韩在线中文字幕| 丰满少妇做爰视频| 亚洲,欧美精品.| 亚洲欧美成人综合另类久久久| 999精品在线视频| 亚洲欧美一区二区三区黑人| 黄片播放在线免费| 波多野结衣av一区二区av| 久久99精品国语久久久| 不卡av一区二区三区| 午夜福利一区二区在线看| 国产熟女午夜一区二区三区| 色婷婷久久久亚洲欧美| av有码第一页| 成人三级做爰电影| 亚洲av电影在线观看一区二区三区| 国产精品99久久99久久久不卡 | 在线亚洲精品国产二区图片欧美| 午夜福利网站1000一区二区三区| 国产1区2区3区精品| 亚洲国产欧美日韩在线播放| 制服丝袜香蕉在线| 日韩成人av中文字幕在线观看| 欧美av亚洲av综合av国产av | 日韩欧美精品免费久久| 飞空精品影院首页| 国产成人精品无人区| 人人妻人人添人人爽欧美一区卜| 国产xxxxx性猛交| a级毛片在线看网站| 日日爽夜夜爽网站| 国产av一区二区精品久久| 观看av在线不卡| 久久精品亚洲熟妇少妇任你| 欧美黄色片欧美黄色片| 久久久久视频综合| 亚洲精品中文字幕在线视频| 中文字幕另类日韩欧美亚洲嫩草| 丰满迷人的少妇在线观看| 伊人久久国产一区二区| 欧美日韩成人在线一区二区| 欧美日韩国产mv在线观看视频| 飞空精品影院首页| 一二三四在线观看免费中文在| 日韩av免费高清视频| 午夜91福利影院| 亚洲av成人不卡在线观看播放网 | 欧美日韩亚洲综合一区二区三区_| 自拍欧美九色日韩亚洲蝌蚪91| 一边摸一边抽搐一进一出视频| 国产日韩欧美在线精品| 大陆偷拍与自拍| 女人精品久久久久毛片| 久久午夜综合久久蜜桃| 久久久久精品人妻al黑| 亚洲美女黄色视频免费看| 中文字幕最新亚洲高清| 九色亚洲精品在线播放| av在线app专区| 国产又色又爽无遮挡免| 狂野欧美激情性xxxx| 久久久久久免费高清国产稀缺| 制服丝袜香蕉在线| 丝袜美足系列| 午夜福利,免费看| 成人国语在线视频| 国产精品久久久久久久久免| 欧美在线一区亚洲| 9191精品国产免费久久| 深夜精品福利| 美女扒开内裤让男人捅视频| 水蜜桃什么品种好| 看十八女毛片水多多多| 午夜久久久在线观看| 丝袜脚勾引网站| 大香蕉久久网| 精品一区二区三区av网在线观看 | 精品人妻熟女毛片av久久网站| 亚洲国产欧美在线一区| 久久久久久免费高清国产稀缺| 国产成人精品无人区| 日日爽夜夜爽网站| 蜜桃在线观看..| 最新在线观看一区二区三区 | 黄片无遮挡物在线观看| 欧美精品亚洲一区二区| 夫妻性生交免费视频一级片| 我的亚洲天堂| 男女无遮挡免费网站观看| bbb黄色大片| 国产精品一二三区在线看| 51午夜福利影视在线观看| 亚洲 欧美一区二区三区| 老汉色∧v一级毛片| 午夜福利网站1000一区二区三区| 欧美日韩国产mv在线观看视频| 亚洲在久久综合| 一级毛片电影观看| 美女中出高潮动态图| 精品国产乱码久久久久久小说| 9热在线视频观看99| 日韩熟女老妇一区二区性免费视频| 日韩精品免费视频一区二区三区| 男女午夜视频在线观看| 国产不卡av网站在线观看| 在线观看免费午夜福利视频| 国产 精品1| 国产av码专区亚洲av| 亚洲成人手机| 岛国毛片在线播放| 丰满少妇做爰视频| 成人毛片60女人毛片免费| 老司机深夜福利视频在线观看 | 精品一区二区三区四区五区乱码 | 美女中出高潮动态图| 免费久久久久久久精品成人欧美视频| 搡老岳熟女国产| 亚洲av福利一区| 精品一品国产午夜福利视频| av一本久久久久| 欧美黄色片欧美黄色片| 日日摸夜夜添夜夜爱| 午夜福利在线免费观看网站| 2018国产大陆天天弄谢| 一本久久精品| 青春草视频在线免费观看| 国产探花极品一区二区| 亚洲人成电影观看| 中国国产av一级| 秋霞在线观看毛片| 无遮挡黄片免费观看| 午夜福利网站1000一区二区三区| 99国产综合亚洲精品| 国产av一区二区精品久久| 亚洲av中文av极速乱| 精品少妇内射三级| 免费少妇av软件| 国产精品久久久久久人妻精品电影 | 亚洲美女黄色视频免费看| 一级毛片我不卡| 国产xxxxx性猛交| 91国产中文字幕| 欧美久久黑人一区二区| 新久久久久国产一级毛片| 久久久久久久久久久久大奶| 我的亚洲天堂| 考比视频在线观看| 日韩大片免费观看网站| 悠悠久久av| 亚洲一码二码三码区别大吗| 色播在线永久视频| 精品一品国产午夜福利视频| 欧美日韩一区二区视频在线观看视频在线| 伊人亚洲综合成人网| 亚洲精品第二区| 亚洲色图 男人天堂 中文字幕| 久久久久久久国产电影| 少妇的丰满在线观看| 午夜免费男女啪啪视频观看| 久久97久久精品| 人体艺术视频欧美日本| 午夜激情av网站| 亚洲人成电影观看| 老熟女久久久| 午夜激情久久久久久久| 免费在线观看完整版高清| 中文字幕av电影在线播放| 黄频高清免费视频| 一级黄片播放器| 别揉我奶头~嗯~啊~动态视频 | 国产精品嫩草影院av在线观看| 一级片免费观看大全| 亚洲成人一二三区av| 国产一区二区在线观看av| 精品一区二区三卡| 人人澡人人妻人| 黑丝袜美女国产一区| 中文乱码字字幕精品一区二区三区| 又大又爽又粗| 免费黄网站久久成人精品| 精品久久蜜臀av无| 欧美激情 高清一区二区三区| 国产成人免费无遮挡视频| 免费观看av网站的网址| 欧美成人午夜精品| 国产成人91sexporn| 一区福利在线观看| 看免费av毛片| 国产伦人伦偷精品视频| 国产日韩欧美视频二区| 久久这里只有精品19| 男女无遮挡免费网站观看| 欧美久久黑人一区二区| 国产精品av久久久久免费| 欧美精品亚洲一区二区| 99精品久久久久人妻精品| 无遮挡黄片免费观看| 成人18禁高潮啪啪吃奶动态图| 伦理电影大哥的女人| 欧美国产精品va在线观看不卡| 亚洲人成电影观看| 精品亚洲成国产av| 久久ye,这里只有精品| 欧美最新免费一区二区三区| 纵有疾风起免费观看全集完整版| 精品久久蜜臀av无| 亚洲成人一二三区av| 午夜福利网站1000一区二区三区| 99久久99久久久精品蜜桃| 人妻 亚洲 视频| 亚洲欧洲国产日韩| 水蜜桃什么品种好| 亚洲七黄色美女视频| 国产av国产精品国产| 国产色婷婷99| 亚洲国产av影院在线观看| 天堂俺去俺来也www色官网| 看十八女毛片水多多多| 国产人伦9x9x在线观看| 国产又爽黄色视频| 日本wwww免费看| 亚洲精品在线美女| av网站免费在线观看视频| 久久这里只有精品19| svipshipincom国产片| 亚洲国产成人一精品久久久| 久久天躁狠狠躁夜夜2o2o | 中文欧美无线码| 国产成人a∨麻豆精品| 成人午夜精彩视频在线观看| 国产麻豆69| 在现免费观看毛片| 精品国产露脸久久av麻豆| 欧美久久黑人一区二区| 熟女av电影| 日本爱情动作片www.在线观看| 国产一区二区三区综合在线观看| 亚洲国产欧美在线一区| 97在线人人人人妻| 日韩一卡2卡3卡4卡2021年| 黑人巨大精品欧美一区二区蜜桃| 亚洲一区二区三区欧美精品| 最近2019中文字幕mv第一页| 国产极品天堂在线| 国产精品免费视频内射| 国产黄色免费在线视频| 韩国高清视频一区二区三区| av在线观看视频网站免费| 97精品久久久久久久久久精品| 男的添女的下面高潮视频| 韩国av在线不卡| 各种免费的搞黄视频| 亚洲成av片中文字幕在线观看| 2021少妇久久久久久久久久久| 制服诱惑二区| 男女无遮挡免费网站观看| 丝袜美足系列| 最新的欧美精品一区二区| 涩涩av久久男人的天堂| 看免费av毛片| 男女边摸边吃奶| 午夜福利视频在线观看免费| 天美传媒精品一区二区| 电影成人av| 91成人精品电影| 色视频在线一区二区三区| 欧美精品一区二区大全| 免费黄网站久久成人精品| 国产日韩一区二区三区精品不卡| 久久人妻熟女aⅴ| 熟妇人妻不卡中文字幕| svipshipincom国产片| 欧美激情高清一区二区三区 | 日韩一卡2卡3卡4卡2021年| 人人妻人人爽人人添夜夜欢视频| 国产老妇伦熟女老妇高清| 亚洲欧美成人综合另类久久久| 热99久久久久精品小说推荐| 欧美 亚洲 国产 日韩一| 观看av在线不卡| 精品人妻在线不人妻| 咕卡用的链子| 国产精品一国产av| av视频免费观看在线观看| 最近手机中文字幕大全| 看非洲黑人一级黄片| 久久毛片免费看一区二区三区| 亚洲四区av|