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

    Emergency Evacuation Plan of the Louvre Based on Cellular Automata Superposition Model

    2020-10-27 08:15:24XIEJiayong謝嘉泳LIQiqi李琪琪LONGJiali

    XIEJiayong(謝嘉泳),LIQiqi(李琪琪),LONGJiali

    School of Physics and Telecommunication Engineering, South China Normal University, Guangzhou 510006, China

    Abstract: The Louvre is a famous French building and a tourist attraction. It is of great significance to make an emergency evacuation plan for tourists. A social force model and cellular automata were established to simulate the crowd evacuation in unidimensional plane, and then the downstairs model was considered to simulate three-dimensional evacuation. Finally the real evacuation situation was simulated, considering the hazard source as well as other factors. The study found that the exit was closer to the corner, and the evacuation distance of descending stairs could be controlled between 0.3 and 0.6 m to maximize evacuation efficiency. The influence of the hazard at the edge of the channel was greater than that at the center of the channel. Therefore, in the process of evacuation, the staff should be guided to the nearest exit, to maintain the optimal distance, the staff to eliminate the danger source in the center of the passage.

    Key words: evacuation; self-propelled particles; cellular automaton; optimal exit

    Introduction

    Increasing terrorist attacks force us to establish a plan for the evacuation of the Louvre in case of emergency[1-2]. In order to let all the tourists leave the building safely as soon as possible, once the evacuation is informed, the tourists need tond the best exit and choose the best escape route according to their own situation at that time,to empty the building as soon as possible.

    The Louvre has four regular exits and many secret exits, which should be fully utilized in the modeling process. It is important to consider emergency evacuees and tourist exits. At the same time, given the low security of secret exports, it is necessary to measure whether and when to use these additional exports. Therefore, we established a social force model[3-7]and considered the influence of staircase factors and the location of the hazard. What’s more, the simulation of cellular automata travel is introduced into the hazard source location model, which is also one of the biggest innovations in our article. In addition, we also need to consider the inuence of language[8], waiting time and other factors on the evacuation plan[9-11].

    1 Model and Methods

    1.1 A social force model

    Figure 1 illustrates the floor plan of the Louvre and the schematic diagram of our model. The system size of the model isL, and the exit width iss, anddis the distance between the corner and the exit.

    Fig. 1 Social force model: (a) floor plan of Louvre and (b) schematic diagram

    Based on the basic principles of the Vicsek model[12]and the leader-follower model, we establish an extended evacuation model of self-propelled particle system with movable exits. The effective self-organizing evacuation strategy is proposed from the closed space.

    Let the size of the space beL×L, the particle moves in this space with a boundary of velocityv, and the position of exits can be moved randomly. The direction of particle movement is determined by the interaction between groups. The position and direction of particle movement are randomly distributed. The details can be expressed as[12]

    θ=cosθi+sinθj,

    (1)

    whereθis the direction of movement;iandjare the unit vectors. In this limited range [-π,π), the update of particle position at the next moment is also random, namely[12]

    xi(t+1)=xi(t)+xi(t)Δt,

    (2)

    wherexi(t+1) is the location att+1 andviis the speed ofi. The direction of motion of particles can be modified as[13]

    (3)

    wheregiis the preferred direction and the weighting factor isω; |〈θi(t)〉r| represents the average direction vector of all particlesjin the interaction radiusrof particlei, which is further revised as[14-15]

    〈θi(t+1)〉=∑θj(t)|∑θj(t)|.

    (4)

    Considering the time everyone displays to the online application when “moving”, in this case, the weight and its preferred direction should be considered as[13]

    (5)

    wherexsign(t) is the location of the exit signs. Furthergi(t) explanation: whenωapproaches 1, the motion direction of the particle is close to its preferred direction. As it approaches 0, the particle has no preferred direction. That is, its direction of motion is random. When the self-propelled particle hits the wall, it will move under the reflection boundary condition.

    1.2 Effect of stair

    On the basis of the model, a simple staircase model[16-17]is added, which is the preliminary overall model of the emergency evacuation of the Louvre. There are no obstacles in the corridor and stairs. During evacuation, tourists on each floor shall take the shortest route to the temporary stairway, and speed up and form a queue in the shortest time. The reaction time of all tourists is the same, and the body thickness of the people in the queue is the same. In the process of evacuation, if there is detention at the door or the stairway because of bottleneck factors, it shall be processed as waiting in line. As is the case on all floors, the phenomenon of congestion, crossing or falling is ignored. The types of stairs in the Louvre are the same, and they are all core tube stairs.

    The speed of movement on a straight path is[16]

    3.617D1+0.95,

    (6)

    D1=NAp/wLs,

    (7)

    whereApis the horizontal projected area of an individual,ωis the traffic density, andLsis the length of stream of people. The rate of movement down the stairs is[16]

    vD=xdownvL,

    (8)

    xdown=0.775+0.44 e-0.39D1down·

    sin(5.16D1down-0.224).

    (9)

    Given that people in emergency situations usually have this great psychological pressure and cause the change of escape speed, so a revision is needed[16]

    v1=u1vL,

    (10)

    v2=u1v0,

    (11)

    u1=1.49-0.36D1,

    (12)

    wherev1,v2, andvLare the normal stair speed, the speed that people go down stairs in the emergency, and the speed when people move on a straight path. The slant length of each stair is known to beLstair. The effective width of floor platform and intermediate platform isB, and the width of ladder section isb. It is easy to get that the evacuation time of stairwell is[16]

    (13)

    whereLstair,B, andbare the dimensions of stairwell.

    When evacuees on top and bottom meet, two situations will occur. In order to ensure evacuation safety, latecomers need to wait for all the first comers to complete evacuation before starting evacuation. The waiting time corresponding to the two situations is[16]

    (14)

    wherev3is the speed that people are in terrain affected emergency,tw1is the waiting time of the personnel of theMfloor when the evacuees of the first floor arrive at the stairway, andtw2is the waiting time when the evacuees arrive at the stairway after thekfloor.

    Considering that the planar evacuationTtime has been calculated in model, the total evacuation time in this case is[16]

    T=T+Tstair, m-1+tw

    .

    (16)

    1.3 Place of the danger source

    The cellular automata model[18-19], namely the dynamic search model, is introduced into the static model of the model, but only the dynamic hazard source is simulated here. The simulation method for other potential factors, such as bottlenecks and contingency, is the same. At the initial moment, pedestrians are randomly distributed in the evacuation space. Judge the evacuation space of the pedestrian, calculate the mobile income matrix of the pedestrian in the row, and select the cell with the largest mobile income matrix. After the calculation of the mobile revenue matrix for all pedestrians is completed, the position of pedestrians will be updated uniformly. When the pedestrians move to the exit, they will move out of the evacuation space. Repeat the above steps until all pedestrians leave the evacuation space, and the simulation ends.

    Near the hazards consider separately the idealized space, namely (M+2)×(N+2), a two-dimensional(2D) discrete cellular network system, and set up evacuation space for walls on both sides, left and right sides for export. What’s more, pedestrians are separated in different directions by a virtual line. Pedestrians cross this line, which centers on dangerous source, and their radius of influence isR, as shown in Fig. 2[18].

    Fig. 2 Schematic diagram of evacuation space

    Here, hazard area is defined as 2D normal distribution[17]

    (17)

    whereRis the radius of influence of hazard source andσ2is the variance of normal distribution.

    Moore type field is shown in Fig. 3(a). Mobile income matrix of mobile income combination in each field, as shown in Fig. 3(b), is defined as[18]

    Aij=[Eij+(1-a)Dij+aFij](1-ηij),

    (18)

    whereAijis the yield parameter andDijis the default direction return. The export earnin with wall isEijandFijis the direction of derivative yield. The area control parameter isa.

    When the pedestrian is in the normal area,a=0; When pedestrians are in the hazard area,a=1.

    Fig. 3 Rule of cellular automata: (a) Moore type eld and (b) mobile income matrix of mobile income combination in each eld

    When the exit is located in the pedestrian field, the attraction of the meta-cell is infinite, and the pedestrian leaves the evacuation space through the exit as far as possible. The exit benefit is defined as[18]

    (19)

    When pedestrians are located in the hazard area, they will move away from the hazard area as fast as possible. To describe the direction of pedestrians’ escape from the hazard area, the following is[18]

    Fij=-gradf(xij,yij),

    (20)

    whereFijis the directional derivative yield.

    When pedestrians are in the normal area, please move in the default direction. At this point, pedestrians are divided into normal pedestrians and pedestrians flee from the source of danger, and move towards the nearest exit. There is[18]

    1.4 Situation of the intersection

    Then, we consider the corridor connecting the intersection, and give the evacuation method when queuing. In the actual evacuation process, due to panic, conformity and other reasons, the evacuees cannot fully follow the guidance and cooperate with the diversion. According to the population distribution rules withdrawn from the organization, the crowd evacuation process is carried out for each source point: the crowd enter the evacuation road from the source node and proceed along the direction of movement of the convergence point to reach the intersection. And they pass the road in sequence, enter through continuous stretching motion and cross nodes, and finally reach the final node.

    The evacuation speed of people is affected by the degree of traffic congestion. When the degree of congestion is small, the evacuation speed of people basically does not interfere with each other and evacuates at the maximum speed. On the contrary, the evacuation speed will decrease with the increase of the degree of traffic congestion[20], which can be described by exponential function. Therefore, at timet, the evacuation speed of people on section (i,j) can be expressed as

    (22)

    Set the decision variablexijto represent the relationship between the road section and evacuation path is[20]

    (23)

    Then the evacuation time of people flow at nodeois[20]

    (24)

    whereUis the number of nodes, and theowithdare the nodes.

    Therefore, the evacuation timeTffor the stream of people is expressed as the maximum value ofTod, and obviously, it is also the final evacuation completion timetfinalof the entire network of the stream of people, that is[20]

    (25)

    whereTfis the crowed evacuation time.

    2 Results and Discussion

    2.1 Exit and visual field radius

    We change the number of groupN, system sizeL, the positiond, width of the exitsand the weighing factor of preferred directionωto investigate the collective evacuation. Our results of MATLAB simulation are shown in Fig. 4.

    Fig. 4 Rules of the system with no preferred direction(no exit sign ω=0)

    Figures 4 (a), (a1), (a2) are the number of individuals left in the spaceNleftas time stepTfor the different distances between the corner and the exit d; The evacuation proportion of the grouppas a function ofdatT=500 was shown in Figs. 4 (b), (b1), (b2) was shown in Figs. 4 (b), (b1), (b2). In Figs. 4 (c), (c1), (c2),pis a function of different width of the exitswhen the exit is located in the corner and middle of the boundary.

    MATLAB is used to simulate the relationship between visual field radius[21-22]and evacuation time. The results are shown in Table 1.

    Table 1 Relationship between visual field radius and evacuation time

    When the exit is located at the edge of the corner and the middle. In Figs. 4(a), (b), and (c),N=400,L=1.0,s=2.0; in Figs. 4(a1), (b1), and (c1),N=800,L=10,s= 2.0; in Figs. 4(a2), (b2), and (c2),N=400,L=1.0,s=0.2. In Figs. 4 (a), (a1), and (a2), the smaller the distancedbetween the corner and the exit, the higher the degree of population evacuation, indicating that the closer the exit is to the corner, the higher the efficiency of crowd evacuation. Figures 4 (b), (b1), and (b2) also verify this conclusion. According to Figs. 4 (c), (c1), and (c2), no matter where the exit is, the evacuation proportion increases with the increase of exit width, so we can know that a larger exit width is conducive to crowd evacuation.

    It is further inferred from the simulation results that increasing the number of exits will reduce the evacuation efficiency of personnel to a certain extent. It can be seen from Fig. 5 that adding a certain mark at the exit can increase the efficiency of evacuation. While some conclusions have been reached, more consideration needs to be given. In the current simulation system, we only consider the two-dimensional rectangular space and ignore some factors, such as visual angle, noise, and oval stadium. More complex models lead to more complex results.

    Fig. 5 Evacuating time needed for all particles Te as a function of weighting factor with one comer exit (with exit sign)

    According to the above model, MATLAB software is used to compile the program. Taking the actual situation of the Louvre as an example, the parameters are set according to the specific situation, and the time for all evacuees to escape from the building can be obtained.

    It can be seen from Fig. 6 that when the total number of evacuations gradually increases, the total time required for complete evacuation of all evacuations increases with a dotted line, which indicates that the total evacuation time and total evacuation are basically positively correlated if the waiting conditions are not considered. According to Figs. 6 (a) -(b) fold points on the polyline occur when the upper and lower evacuation teams meet and when the evacuation teams on the same floor meet. As the waiting time increases, the total evacuation time increases by a larger margin.When the total number of evacuees continues to increase, the detention waiting time will gradually become the most important factor affecting the total evacuation time.

    Fig. 6 Relation chart of the total evacuation time: (a) total evacuation time and distances and (b) total evacuation time and total number of evacuation

    As shown in Fig. 6, evacuation time decreases first and then increases with the increase of evacuation distance. When the evacuation distance is 0.35 m, the evacuation time reaches the minimum value.When the evacuation spacing is between 0.1 m and 0.6 m, the influence of evacuation spacing on evacuation time is small. Considering the best evacuation efficiency and the degree of recognition of the evacuation interval, and to prevent collisions between evacuees, the 0.5 m evacuation interval suggested above is reasonable[23]. Therefore, in the evacuation exercise or actual evacuation process, it is required that the evacuation distance should be controlled within 0. 3-0. 6 m as far as possible to ensure the high efficiency of evacuation.

    2.2 Place of the danger source

    In the simulation study, the density of pedestrian flowing in the evacuation line is defined as the number of pedestrians in the system and the number of cells in the pedestrian movement interval ratio ofMtoN, pedestrian evacuation timeTis defined as the number of time steps required by all pedestrians in the system to leave the room[24]. The evacuation timeTof the system is the average result of 20 samples. Unless otherwise specified, the parameter isM=20 andN=20. Three representative schemes of hazard location distribution are shown in Fig. 7.

    Fig.7 Three representative schemes of hazard location distribution: (a) the hazard at the wall; (b) the hazard at the upper half of evacuation space; (c) the hazard at the central axis of evacuation space

    Figure 8 shows the graph of the change of evacuation timeTwith hazard source radiusρin the case of pedestrian density andρ=0.9. As can be seen in Fig. 8(a), with the increase ofR, the required evacuation time of Case 1 and Case 3 is getting longer and longer, while the evacuation time of Case 2 decreases gradually afterR>6. As can be seen in Fig. 8 (ρ=0.5), when the density is high, with the increase ofR, the evacuation time of Case 3 becomes longer and longer. Case 2 presents a relatively stable state when 48, and the evacuation time decreases with the increase ofR.

    (a) ρ=0.5 (b) ρ=0.9

    Figure 9 shows the typical spatial and temporal evolution patterns of different hazard radiusRand location distribution schemes when pedestrian densityR=0.9. By observing Figs. 9(a), (b) and (c), whenRis small, Case 1 and Case 2 cause the upper pedestrian to cross the dividing line due to the hazard source blocking the upper pedestrian and form a hedge with the lower pedestrian, which affects the evacuation efficiency. In Case 3, there is a gap between the upper and lower sides of the hazard source, through which pedestrians pass, and the congestion is not obvious. However, with the increase ofR, as shown in Figs. 9(c), (d) and (e), the gap between the upper and lower sides becomes smaller and smaller, so that it is completely blocked. Pedestrians fleeing from the hazard source area are completely blocked together with those in the normal area, which greatly affects the evacuation efficiency. Therefore, the evacuation efficiency of Case 3 decreases with the increase of the radius of the hazard source.

    The influence of the hazard source at the edge of the passage is greater than that at the center of the passage. With the increase of the influence radius of the hazard source located in the center of the passage, the required evacuation time increases linearly.

    In most cases, the larger the scope of the danger source, the lower the evacuation efficiency, and the greater the possibility of congestion.

    (a)

    (b)

    (c)

    (d)

    (e)

    2.3 Situation of the intersection and other factors

    According to the actual road network of the Louvre, the source node, intersection node and sink node are numbered respectively to construct the crowd evacuation road network. The flow of people at each source node is evacuated according to the crowd diversion model constructed. The solving process is shown in Fig. 10.

    Fig. 10 Calculation flowchart

    Dijkstra algorithm is an effective algorithm for solving the classical shortest circuit problem[25]. In this paper, Dijkstra algorithm is used to search the shortest path at the intersection node containing each post-sequence section, and to find the length of the post-sequence shortest remaining path.

    According to the real-time congestion degree of evacuation roads and the order of longest remaining shortest roads, the probability of crowd distribution at a road intersection, the traffic input and output relationship of evacuation roads is analyzed. We have formulated rules for evacuation networks for crowds and traffic distribution at road intersections. In addition, a cross-node cross evacuation crowd flow model is set up after the evacuation, which provides reasonable guidance for emergency crowd evacuation. The allocation of the stream of people at each intersection can reduce the evacuation time consumption of the entire road network, guide the balanced distribution of people flow in space and time, effectively prevent congestion, and improve the evacuation efficiency.

    In the process of evacuation, the influence of evacuation speed, visual field radius and evacuation environment of the disabled should be considered[26]. Using the stairs of the barrier-free passage, the disabled staff members for the first time in the field escape from the barrier-free passage after the emergency. On the premise that the rights and needs of all people in the building are balanced, and the elevator is not occupied by ordinary people, the best choice for the disabled to escape is the elevator, which will play an important role in the shelter of the disabled.

    The average speed of visually impaired people is basically the same when there is a guide for each person and when there is no guide for evacuation. However, in the mode of centralized evacuation by a single guide, the average walking speed is much higher than that in the other two cases. Therefore, single guide is adopted for evacuation of visually impaired people.

    3 Conclusions

    In summary, we have proposed a social force model for the group escape considering movable exits and leader-follower effect on the basis of the Vicsek model. We investigate the evacuation behavior in our model with various system parameters via simulations. The optimal location of exit for collective evacuation is the corner of a limited square-shaped cell. Increasing the intensity of exit sign will speed up the evacuation quickly, and then it will take no further effect when the intensity of sign exceeds a certain value. And multi exits will hinder the evacuation.

    In addition, considering the factors, such as the spacing of the stairs and height under the evacuation, the evacuation spacing control within certain range can maximize the evacuation efficiency. Furthermore,the impact of the hazard source at the edge of the channel is greater than the impact of the hazard source at the center of the channel. Other factors such as tourists’ language and the disabled will also affect the flight time of passengers to a certain extent.

    Our study is rather significant to make more safe environment. Besides, our research is likely to be useful in designing a secure passage for a large enclosed place. Although we have drawn some conclusions, more considerations still need to be done. Our model does not take into account all the factors, such as the possible competition and stampede in the process of tourist evacuation. We will explore some complex models to present more realistic effect in future.

    日韩电影二区| 国产成人a∨麻豆精品| 日本一二三区视频观看| 婷婷色麻豆天堂久久| 国产一区二区三区综合在线观看 | 婷婷色综合www| 日韩欧美三级三区| 午夜爱爱视频在线播放| 国产视频首页在线观看| 日韩不卡一区二区三区视频在线| 国产乱人视频| 午夜福利在线观看免费完整高清在| 高清午夜精品一区二区三区| 麻豆成人av视频| 欧美+日韩+精品| 亚洲国产欧美在线一区| 亚洲av日韩在线播放| 午夜免费男女啪啪视频观看| 欧美日韩亚洲高清精品| 性色avwww在线观看| 天天躁日日操中文字幕| 久久久成人免费电影| 色网站视频免费| 日韩欧美一区视频在线观看 | or卡值多少钱| 国产一区亚洲一区在线观看| 男女下面进入的视频免费午夜| 欧美日韩一区二区视频在线观看视频在线 | 亚洲婷婷狠狠爱综合网| 久久久色成人| 国产人妻一区二区三区在| 深爱激情五月婷婷| 国产成人福利小说| 色哟哟·www| 国产成人福利小说| av在线播放精品| 国产亚洲精品久久久com| 午夜福利视频1000在线观看| 国产精品久久久久久av不卡| 嫩草影院新地址| 久久久久久伊人网av| 高清日韩中文字幕在线| 精品久久久久久电影网| 中文欧美无线码| 99久久九九国产精品国产免费| 天天一区二区日本电影三级| 亚洲三级黄色毛片| 九九爱精品视频在线观看| 丰满少妇做爰视频| 激情 狠狠 欧美| 99久久中文字幕三级久久日本| 天堂中文最新版在线下载 | 大片免费播放器 马上看| 男女国产视频网站| 亚洲精品aⅴ在线观看| 日韩,欧美,国产一区二区三区| 人妻一区二区av| 国产午夜精品一二区理论片| 亚洲精品国产av蜜桃| 精品一区二区三区视频在线| 最近的中文字幕免费完整| 国产视频内射| 波多野结衣巨乳人妻| 欧美日韩一区二区视频在线观看视频在线 | 国产淫片久久久久久久久| 嫩草影院精品99| 综合色av麻豆| 日本av手机在线免费观看| 亚洲国产最新在线播放| 精品人妻一区二区三区麻豆| 欧美最新免费一区二区三区| 国产精品福利在线免费观看| 国产精品爽爽va在线观看网站| 美女cb高潮喷水在线观看| 最近的中文字幕免费完整| 又黄又爽又刺激的免费视频.| 内射极品少妇av片p| 天天躁夜夜躁狠狠久久av| 99热网站在线观看| 免费观看无遮挡的男女| 久热久热在线精品观看| a级一级毛片免费在线观看| 成年版毛片免费区| 深夜a级毛片| 色综合站精品国产| 午夜爱爱视频在线播放| 欧美潮喷喷水| 一区二区三区高清视频在线| 日韩大片免费观看网站| 毛片一级片免费看久久久久| 18禁动态无遮挡网站| 赤兔流量卡办理| 夜夜看夜夜爽夜夜摸| av免费观看日本| 国产视频首页在线观看| 女人十人毛片免费观看3o分钟| 亚洲欧美一区二区三区黑人 | 91精品一卡2卡3卡4卡| 搞女人的毛片| 免费看光身美女| 日韩,欧美,国产一区二区三区| 成人午夜精彩视频在线观看| 国产免费又黄又爽又色| 精品一区二区三卡| 亚洲人成网站在线观看播放| 免费黄网站久久成人精品| 久久精品久久久久久久性| 国产美女午夜福利| 一区二区三区免费毛片| 成人午夜高清在线视频| 日本一二三区视频观看| 久久99热6这里只有精品| 亚洲国产精品国产精品| 大又大粗又爽又黄少妇毛片口| 成人欧美大片| 日韩视频在线欧美| 国产探花在线观看一区二区| 国产免费又黄又爽又色| 天堂av国产一区二区熟女人妻| 亚洲丝袜综合中文字幕| 午夜视频国产福利| 天堂俺去俺来也www色官网 | 干丝袜人妻中文字幕| 美女高潮的动态| 国产人妻一区二区三区在| 久久午夜福利片| 一区二区三区免费毛片| 老师上课跳d突然被开到最大视频| 黄色一级大片看看| 国产一区二区在线观看日韩| 六月丁香七月| 成年av动漫网址| 在线播放无遮挡| 午夜日本视频在线| 九色成人免费人妻av| 国产视频首页在线观看| 成人鲁丝片一二三区免费| 人人妻人人看人人澡| 久久久亚洲精品成人影院| 亚洲国产欧美在线一区| 简卡轻食公司| 久久久午夜欧美精品| 3wmmmm亚洲av在线观看| 五月玫瑰六月丁香| 精品人妻熟女av久视频| 成人午夜精彩视频在线观看| 国产精品国产三级专区第一集| 国国产精品蜜臀av免费| 乱系列少妇在线播放| 国产亚洲av片在线观看秒播厂 | 久久精品久久久久久噜噜老黄| 国产精品女同一区二区软件| 麻豆乱淫一区二区| 十八禁网站网址无遮挡 | 亚洲av电影在线观看一区二区三区 | 午夜福利在线在线| 激情五月婷婷亚洲| 三级国产精品片| 青春草国产在线视频| 久久久久久久久中文| 91午夜精品亚洲一区二区三区| 国产国拍精品亚洲av在线观看| 国产亚洲av嫩草精品影院| 老师上课跳d突然被开到最大视频| 国产av码专区亚洲av| 深夜a级毛片| 欧美成人a在线观看| 欧美一级a爱片免费观看看| 日韩精品青青久久久久久| 国内精品美女久久久久久| 亚洲精品日韩在线中文字幕| 久久6这里有精品| 乱码一卡2卡4卡精品| 日韩av在线免费看完整版不卡| 国产黄片视频在线免费观看| 青春草视频在线免费观看| 国产精品无大码| 亚洲aⅴ乱码一区二区在线播放| 中文在线观看免费www的网站| 成人毛片a级毛片在线播放| 国产成年人精品一区二区| 亚洲国产精品成人久久小说| 国产亚洲最大av| 春色校园在线视频观看| 99热6这里只有精品| 精品久久久久久久人妻蜜臀av| 国产爱豆传媒在线观看| 黑人高潮一二区| 成年人午夜在线观看视频 | 国产成人freesex在线| 身体一侧抽搐| 国产精品人妻久久久影院| 卡戴珊不雅视频在线播放| 18禁裸乳无遮挡免费网站照片| 又爽又黄a免费视频| 永久网站在线| 777米奇影视久久| 欧美日韩国产mv在线观看视频 | 午夜激情久久久久久久| 免费观看无遮挡的男女| 最近最新中文字幕大全电影3| 一区二区三区高清视频在线| 国产成人福利小说| 日日摸夜夜添夜夜爱| 国产精品美女特级片免费视频播放器| 尾随美女入室| 国产亚洲5aaaaa淫片| 麻豆国产97在线/欧美| 97热精品久久久久久| 国产精品国产三级国产专区5o| 国内精品一区二区在线观看| 国产精品福利在线免费观看| 夫妻午夜视频| 少妇人妻精品综合一区二区| 亚洲丝袜综合中文字幕| 熟女电影av网| 精品一区二区免费观看| 亚洲综合精品二区| 成年女人在线观看亚洲视频 | 人人妻人人澡欧美一区二区| 少妇裸体淫交视频免费看高清| 99久国产av精品| 白带黄色成豆腐渣| 国内揄拍国产精品人妻在线| 观看美女的网站| 狂野欧美白嫩少妇大欣赏| 久久久精品免费免费高清| 国产永久视频网站| 亚洲内射少妇av| 能在线免费观看的黄片| 美女被艹到高潮喷水动态| 欧美xxⅹ黑人| 国产亚洲91精品色在线| 乱系列少妇在线播放| 国产成人午夜福利电影在线观看| 亚洲精品乱码久久久久久按摩| 欧美精品国产亚洲| 午夜爱爱视频在线播放| 午夜免费观看性视频| 精品国产露脸久久av麻豆 | 国产黄色视频一区二区在线观看| 国产欧美日韩精品一区二区| 午夜福利高清视频| 久久久成人免费电影| 最近2019中文字幕mv第一页| 国产伦精品一区二区三区四那| 国产乱人视频| 国产成人精品婷婷| 干丝袜人妻中文字幕| 国产伦一二天堂av在线观看| 亚洲成人中文字幕在线播放| 青春草亚洲视频在线观看| 99热这里只有精品一区| 亚洲色图av天堂| 看十八女毛片水多多多| 肉色欧美久久久久久久蜜桃 | 黄片无遮挡物在线观看| 成人漫画全彩无遮挡| 国产精品伦人一区二区| 国产精品一区二区三区四区免费观看| 日韩精品有码人妻一区| 少妇猛男粗大的猛烈进出视频 | 久久精品熟女亚洲av麻豆精品 | 欧美性猛交╳xxx乱大交人| 99久国产av精品| 国产精品1区2区在线观看.| 精品人妻视频免费看| 亚洲av电影在线观看一区二区三区 | 2021天堂中文幕一二区在线观| 中国国产av一级| 97精品久久久久久久久久精品| 久久这里有精品视频免费| av网站免费在线观看视频 | 国产淫片久久久久久久久| 精品久久久久久久人妻蜜臀av| 国产三级在线视频| 91在线精品国自产拍蜜月| 亚洲不卡免费看| 看十八女毛片水多多多| 精品久久久久久电影网| 青春草国产在线视频| 纵有疾风起免费观看全集完整版 | 特级一级黄色大片| 国产亚洲最大av| 成人美女网站在线观看视频| 久久久久久久久久黄片| videos熟女内射| 97人妻精品一区二区三区麻豆| 简卡轻食公司| 欧美精品一区二区大全| 成年人午夜在线观看视频 | 成人午夜高清在线视频| 亚洲欧美日韩东京热| 丰满乱子伦码专区| 99热这里只有精品一区| 国产高潮美女av| 欧美激情在线99| 国产精品人妻久久久久久| 午夜激情久久久久久久| 国产在视频线在精品| 亚洲欧美一区二区三区黑人 | 久久精品国产鲁丝片午夜精品| 亚洲精品久久久久久婷婷小说| 99久久精品国产国产毛片| 三级经典国产精品| 久久久精品欧美日韩精品| 成人综合一区亚洲| 亚洲av成人av| 精品国内亚洲2022精品成人| 中文精品一卡2卡3卡4更新| 尤物成人国产欧美一区二区三区| 97人妻精品一区二区三区麻豆| 韩国av在线不卡| 夜夜看夜夜爽夜夜摸| 性插视频无遮挡在线免费观看| 日韩成人伦理影院| 日韩强制内射视频| 小蜜桃在线观看免费完整版高清| 国产男人的电影天堂91| 老女人水多毛片| 婷婷六月久久综合丁香| 国产久久久一区二区三区| av播播在线观看一区| 国产麻豆成人av免费视频| 午夜福利网站1000一区二区三区| 精品一区二区三区视频在线| 国产高潮美女av| 亚洲四区av| 黄色一级大片看看| 少妇猛男粗大的猛烈进出视频 | 国产一区二区在线观看日韩| 亚洲av不卡在线观看| 欧美日韩在线观看h| 少妇熟女欧美另类| 亚洲精品视频女| 能在线免费观看的黄片| 日韩av不卡免费在线播放| 99久国产av精品国产电影| 美女脱内裤让男人舔精品视频| 天天躁日日操中文字幕| h日本视频在线播放| 校园人妻丝袜中文字幕| 99re6热这里在线精品视频| 久久久国产一区二区| 国产日韩欧美在线精品| 中文精品一卡2卡3卡4更新| av在线老鸭窝| 又粗又硬又长又爽又黄的视频| 乱码一卡2卡4卡精品| 菩萨蛮人人尽说江南好唐韦庄| 亚洲人成网站在线播| 黄片wwwwww| 超碰97精品在线观看| 2021天堂中文幕一二区在线观| 在线免费观看不下载黄p国产| 国产免费一级a男人的天堂| 亚洲最大成人手机在线| 永久免费av网站大全| 午夜福利视频1000在线观看| 身体一侧抽搐| 国产成人aa在线观看| 国产高潮美女av| 国产淫片久久久久久久久| 国产精品一区二区性色av| 午夜精品在线福利| 一区二区三区高清视频在线| 简卡轻食公司| 91午夜精品亚洲一区二区三区| 久99久视频精品免费| 寂寞人妻少妇视频99o| 亚洲最大成人手机在线| 91久久精品电影网| 久久精品夜色国产| 亚洲乱码一区二区免费版| 非洲黑人性xxxx精品又粗又长| 午夜日本视频在线| 亚洲一区高清亚洲精品| 免费看a级黄色片| 少妇熟女欧美另类| 91午夜精品亚洲一区二区三区| 久久久久网色| 中文字幕av成人在线电影| 国产黄色视频一区二区在线观看| 身体一侧抽搐| 久久精品夜色国产| 亚洲欧美精品专区久久| 中文字幕制服av| 亚洲熟女精品中文字幕| 男人狂女人下面高潮的视频| 亚洲精品日韩av片在线观看| 蜜桃亚洲精品一区二区三区| 少妇被粗大猛烈的视频| 又爽又黄无遮挡网站| 日韩欧美精品v在线| 免费看日本二区| 欧美日韩精品成人综合77777| 嘟嘟电影网在线观看| 国产黄频视频在线观看| 国产v大片淫在线免费观看| 国产 一区精品| 久久草成人影院| 少妇丰满av| 亚洲久久久久久中文字幕| 少妇丰满av| 99久久人妻综合| 日韩制服骚丝袜av| 丝袜喷水一区| 国产av不卡久久| 午夜视频国产福利| 日本免费a在线| 久久精品国产亚洲av涩爱| 国产av国产精品国产| 国产白丝娇喘喷水9色精品| 国产成人a∨麻豆精品| 在线a可以看的网站| 插阴视频在线观看视频| 国产精品.久久久| 国产黄片视频在线免费观看| 成人高潮视频无遮挡免费网站| 亚洲欧美日韩卡通动漫| 熟女电影av网| 熟女人妻精品中文字幕| 欧美日韩在线观看h| 亚洲av电影不卡..在线观看| av免费在线看不卡| 欧美97在线视频| 国产精品人妻久久久影院| 成人综合一区亚洲| 久久久精品欧美日韩精品| 十八禁网站网址无遮挡 | av在线蜜桃| 国产麻豆成人av免费视频| av在线观看视频网站免费| 蜜桃亚洲精品一区二区三区| 热99在线观看视频| 国产 一区 欧美 日韩| 黑人高潮一二区| 午夜免费观看性视频| 99九九线精品视频在线观看视频| 国产伦在线观看视频一区| 国精品久久久久久国模美| 国产精品.久久久| 一本久久精品| 日韩一本色道免费dvd| 99热全是精品| 久久久久精品久久久久真实原创| 国产高清三级在线| 亚洲av成人av| 国产伦在线观看视频一区| 丝袜美腿在线中文| 亚洲第一区二区三区不卡| 亚洲国产精品成人久久小说| 黄片无遮挡物在线观看| 18禁裸乳无遮挡免费网站照片| 国产女主播在线喷水免费视频网站 | 夫妻性生交免费视频一级片| 丰满少妇做爰视频| 亚洲第一区二区三区不卡| 亚洲人与动物交配视频| 欧美日韩精品成人综合77777| 久热久热在线精品观看| 国产 一区精品| 国产精品国产三级国产专区5o| 精品久久久久久久久亚洲| 久久久久久久久久久丰满| 欧美精品国产亚洲| 麻豆成人av视频| 亚洲精品影视一区二区三区av| 夫妻性生交免费视频一级片| 亚洲精品日本国产第一区| 日韩av在线免费看完整版不卡| 麻豆av噜噜一区二区三区| 搡女人真爽免费视频火全软件| 成人美女网站在线观看视频| 欧美性感艳星| 国产老妇女一区| 亚洲国产色片| 国产午夜精品久久久久久一区二区三区| 蜜桃久久精品国产亚洲av| 亚洲18禁久久av| 久久这里有精品视频免费| 国产黄色小视频在线观看| 国产成人精品福利久久| 一本一本综合久久| 91aial.com中文字幕在线观看| 亚洲不卡免费看| 国产精品国产三级国产av玫瑰| 国产日韩欧美在线精品| 亚洲av一区综合| 亚洲经典国产精华液单| 成人性生交大片免费视频hd| 熟妇人妻不卡中文字幕| 亚洲精品,欧美精品| 狠狠精品人妻久久久久久综合| 美女黄网站色视频| 亚洲国产日韩欧美精品在线观看| 男女视频在线观看网站免费| 天堂中文最新版在线下载 | 亚洲一区高清亚洲精品| 99久久中文字幕三级久久日本| 国产午夜精品一二区理论片| 国产高清国产精品国产三级 | 伦精品一区二区三区| 丰满少妇做爰视频| 神马国产精品三级电影在线观看| av线在线观看网站| 我的女老师完整版在线观看| 91精品一卡2卡3卡4卡| 丝袜美腿在线中文| av网站免费在线观看视频 | 国产黄色免费在线视频| 婷婷色av中文字幕| 国产精品人妻久久久影院| 亚洲精品一二三| 久久精品久久精品一区二区三区| 身体一侧抽搐| 国产精品蜜桃在线观看| 免费电影在线观看免费观看| 色5月婷婷丁香| 精品久久久久久久末码| 久久热精品热| 亚洲不卡免费看| 久久久久久国产a免费观看| 免费看光身美女| 亚洲欧美清纯卡通| 蜜臀久久99精品久久宅男| 亚洲不卡免费看| 久久久久久久午夜电影| 波野结衣二区三区在线| 亚洲精品成人av观看孕妇| 黄片wwwwww| 成人一区二区视频在线观看| 成人午夜高清在线视频| 国产爱豆传媒在线观看| 亚洲av不卡在线观看| 少妇的逼好多水| 亚洲国产日韩欧美精品在线观看| 国精品久久久久久国模美| 永久免费av网站大全| 久久久久久久大尺度免费视频| 久久久国产一区二区| 亚洲无线观看免费| 大香蕉97超碰在线| 亚洲自拍偷在线| 国产综合懂色| 久久综合国产亚洲精品| 你懂的网址亚洲精品在线观看| 99热全是精品| 国产成人91sexporn| 午夜免费男女啪啪视频观看| 69av精品久久久久久| 在线观看av片永久免费下载| 18禁动态无遮挡网站| 麻豆国产97在线/欧美| 亚洲最大成人中文| 亚洲四区av| 婷婷色综合大香蕉| 天堂av国产一区二区熟女人妻| av卡一久久| 天天躁夜夜躁狠狠久久av| 中文字幕免费在线视频6| 国产黄a三级三级三级人| 白带黄色成豆腐渣| 国产成人aa在线观看| 久久久成人免费电影| 国产伦理片在线播放av一区| 99久国产av精品国产电影| 欧美日韩精品成人综合77777| 秋霞伦理黄片| 嫩草影院入口| 日本与韩国留学比较| 十八禁国产超污无遮挡网站| 在线观看美女被高潮喷水网站| 亚洲av日韩在线播放| 精品99又大又爽又粗少妇毛片| 久久久精品94久久精品| 国产av码专区亚洲av| 免费看av在线观看网站| 一边亲一边摸免费视频| 国产一区二区三区综合在线观看 | 午夜福利视频1000在线观看| 精华霜和精华液先用哪个| 久久这里有精品视频免费| av国产久精品久网站免费入址| 国产三级在线视频| 大又大粗又爽又黄少妇毛片口| 毛片一级片免费看久久久久| 男女边摸边吃奶| 中文字幕人妻熟人妻熟丝袜美| 久久综合国产亚洲精品| 亚洲精品aⅴ在线观看| 别揉我奶头 嗯啊视频| 免费看美女性在线毛片视频| 国产免费福利视频在线观看| 成年av动漫网址| 亚洲av电影在线观看一区二区三区 | 午夜日本视频在线| 又黄又爽又刺激的免费视频.| 亚洲精品久久午夜乱码| 久久精品人妻少妇| 真实男女啪啪啪动态图| 欧美激情在线99| 国产亚洲5aaaaa淫片| 亚洲三级黄色毛片| 一区二区三区免费毛片| 亚洲欧美日韩无卡精品| a级一级毛片免费在线观看| 97热精品久久久久久| 亚洲熟女精品中文字幕| 亚洲丝袜综合中文字幕| 亚洲婷婷狠狠爱综合网| 日韩欧美精品v在线| 成人午夜精彩视频在线观看| 欧美成人a在线观看| 婷婷色麻豆天堂久久|