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

    Modeling and analysis of car-following behavior considering backward-looking effect?

    2021-03-19 03:23:48DongfangMa馬東方YueyiHan韓月一FengzhongQu瞿逢重andShengJin金盛
    Chinese Physics B 2021年3期

    Dongfang Ma(馬東方), Yueyi Han(韓月一), Fengzhong Qu(瞿逢重), and Sheng Jin(金盛)

    1Ocean College,Institute of Marine Information Science and Technology,Zhejiang University,Hangzhou 310058,China

    2College of Civil Engineering and Architecture,Institute of Intelligent Transportation System,Zhejiang University,Hangzhou 310058,China

    Keywords: traffic flow,car-following model,visual angle,backward-looking effect

    1. Introduction

    Car-following models have always been an important research topic in the field of microscopic traffic flow and statistical physics. It has significant application in complex systems,traffic simulation,and intelligent driving.Numerous carfollowing models have been developed in the past decades,and they have important theoretical significance for a profound understanding of complex random characteristics,critical phase transitions,and stability of traffic flow.

    Since the concept of car-following model was proposed in the early 1950s,hundreds of important car-following models have been proposed.[1-8]In recent years, many different car-following models[9-18]are put forward under the modern traffic situation and lots of scholars devote themselves in investigating the behavior of car following.[19]To enrich the traffic flow theory, some scholars proposed extended traffic flow models.[20]However, most models only consider the driver’s perception of the leading vehicle’s motion status information and completely ignore the stimulation of rear vehicle’s motion.In an actual driving environment, the driver will also observe the movement of rear vehicle through the rearview mirror and use it as a basis for adjusting the speed of own vehicle to prevent the rear vehicle from rear-end collision due to its own emergency braking. Especially when traffic is congested or vehicles frequently change lanes,the rear-view effect is more obvious. At the same time, when the distance between two vehicles in the front and rear is small,a severe deceleration of the front vehicle may cause the rear vehicle to brake urgently or even cause a rear-end collision, which may cause congestion and security risks. Therefore, to ensure their own safety and smooth traffic,the drivers will comprehensively judge the state of leading and following vehicles to adjust their driving behavior when driving.The rear-view effect helps the driver to maintain the distance between front and rear vehicles as consistent as possible which effectively stabilizes the traffic flow.

    In recent years,the effect of rear car’s movement on carfollowing behavior has been gradually considered, and relevant models have been established. Sun et al.[21]proposed a car-following model by considering the backward-looking effect and speed difference and showed that the consideration of above factors simultaneously can enhance the stability of traffic flow and effectively suppress the traffic jam. Ge et al.[22]proposed a feedback control car-following model that considers the headway of front and rear vehicles and conducted a stability analysis of the model from a cybernetic perspective, confirming the important role of feedback control signals in suppressing and alleviating traffic congestion. Chen et al.[23]proposed a car-following model that considers driver’s sensory memory and the backward-looking effect, exhibiting good traffic flow stabilization performance. Based on the GIPPS safety distance model, Yang et al.[24]further considered the backward-looking effect and developed a new model.They found that the consideration of backward-looking effect has three types of effects on traffic flow:stabilization,destabilization, and nonphysical phenomena. These models provide important reference value for understanding the effect of rear car on the car-following behavior and have important significance for understanding the factors influencing car-following behavior in a complex environment.

    However,all the above-mentioned models assume that the drivers can accurately perceive the speed, distance, and other variables of the rear car, contrary to the driver’s perception characteristics.Related studies in the field of psychology[25,26]show that the driver cannot accurately perceive the state parameters such as speed and distance; therefore, the following car model based on accurate parameters cannot describe the real driving behavior well. In an actual driving environment,the driver mainly relies on visual information to perceive the vehicle’s movement and then changes his or her own acceleration. In recent years, studies have also been conducted on the visual angle car-following model.[27-30]Jin et al.[28]introduced visual angle as a stimulus for driver to perceive the leading vehicle and built a relevant car-following model. The simulation results show that the model can better describe the characteristics of driver’s asymmetry in acceleration and deceleration and speed fluctuation. Recently, Ma et al.[29]proposed that the solid angle composed of the width and height of the vehicle is a factor influencing the car-following behavior and established the new model using the solid angle as the influencing factor. However,all the car-following models based on visual information focus on the influence of leading vehicle. It is still theoretically significant how to reconstruct the car-following behavior under the combined influence of front and rear vehicles and define more reliable stability conditions of traffic flow parameters.

    Based on the above analysis,this study assumes that during the car-following behavior,the driver judges the movement of leading car and following car depending on visual information and establishes the calculation method of visual angle and its rate of change of front and rear vehicles. On this basis,a fusion mechanism is proposed for the visual information of front and rear vehicles,and then a car-following model is developed under the influence of backward-looking effect. Further analysis of model characteristics and advantages such as stability analysis, numerical simulation, and model comparison provides a basis for building a more realistic and effective microscopic traffic flow model.Using visual angle as influencing factor is a better indicator of human drivers’ awareness.Especially for the rear car, the drivers are harder to perceive the headway and the velocity. The only way to judge the rear car’s behavior is observing the change of the rear car’s size in the driving mirror. Therefore, adopting visual angle in the model can more truly describe the backward-looking effect.Although the theoretical analysis and numerical simulation are similar to those in previous works, there are differences in some aspects such as the asymmetry of the stability and the size of the stability region. So that the results in our paper are more consistent with the actual situation. The backwardlooking effect formed by visual angle can better describe the influence of the oversize vehicles. In reality,the oversize vehicles such as trucks and buses occupy certain proportion which the traditional models cannot reflect the influence of.

    2. Model

    In addition to adjusting the speed and acceleration by observing the leader’s movement,the driver will also observe the change in the movement of rear vehicle through the rearview mirror. Especially when the rear vehicle is a large vehicle,the driver is more willing to increase the distance from the rear car to improve safety.

    Therefore, this study proposes a car-following model framework that considers both front and rear vehicle’s visual information. In this framework,the vehicle size can be introduced into the new model. During the car-following process,the driver’s perception of visual angle information of front and rear vehicles is shown in Fig.1.

    Fig.1. Driver’s sensitivity of forward and backward visual information.

    The formula for calculating the visual angle of leading vehicle (n+1)-th vehicle) to the target vehicle (n-th vehicle)can be given as follows:

    where θnF(t) means the visual angle formed by the n-th car and the car in front of it,wn+1is the width of leading vehicle and ln+1is the length of leading vehicle. Δxn(t)is the distance between n-th car and its leading car at time t. The calculation formula is as follows:

    Similarly, the formula for calculating the visual angle of following vehicle(the(n?1)-th vehicle)relative to the target vehicle can be expressed as follows:

    where θnB(t) means the visual angle formed by the n-th car and the car behind it and wn?1is the width of the following vehicle.

    The driver optimizes the speed according to the visual angle of front and rear vehicles;therefore,the optimized velocity function can be expressed as follows:

    Equation (4) shows the optimal velocity function determined by the forward visual information, while equation (5)shows the optimal velocity function determined by the backward visual information. V1F, V2F, C1F, C2F, V1B, V2B, C1B,and C2Bare the parameters of optimal velocity function to be calibrated. The values are V1F=V1B=6.75 m/s,V2F=V2B=7.91 m/s,C1F=C1B=0.13 m?1,C2F=C2B=1.57.

    Equations(4)and(5)show that the expressions of VBand VFare almost the same. The main difference is the difference in overall value direction because the influence of backwardlooking and forward-looking effects on car-following behavior is opposite.For example,when the front car is close to the current car,the forward visual angle increases,and the current car will respond to deceleration. When the rear car is close to the current car,the rear visual angle will also increase,but the current car will respond to acceleration.The reason why the same optimal velocity function as the function in previous work is adopted is to contrast the new model with the FVD model directly.It is perfectly acceptable to assume the optimal velocity function of other forms rather than the form of headway, but then the parameters must be recalibrated.

    Combining the influence of optimized velocity of front and rear vehicles and the rate of change of visual angle, the car-following model under the influence of backward-looking effect can be expressed as follows:

    where p is the parameter of influencing factors of backwardlooking effect. When p=0, the driver is not influenced by backward-looking effect;α is the sensitivity coefficient,which is the reciprocal of reaction time τ;λ is the response intensity coefficient to the stimulation of visual angle change rate. The value of λ is larger than α, which need to be explained as follows. Firstly,the value of λ is derived according to the calibration of the actual driving data(from NGSIM)and the general range is 40-60, which is not set casually. Secondly, the value of λ is unable to decide whether the changing rate of visual angle(as shown in Eq.(7))can be neglected. Particularly,in the case of the relatively small headway, the changing rate of the forward visual angle can get large, which affected the acceleration deeply. Thirdly, the changing rate of visual angle is used as the parameter influencing acceleration and has a significant influence on the stability analysis result.

    Equation(6)describes the calculation method of acceleration of target vehicle with the front and rear vehicle visual angle information as the stimulus factor. The optimized velocity in the model is obtained from the weighted average of optimized velocities decided by the forward and backward visual angles,i.e.,the driver can obtain the optimal velocity based on the forward and backward visual angle information at the same time and adjust his or her own speed. Essentially, this model indicates that the driver comprehensively considers the visual angle information of front and rear vehicles and tries to maintain the same visual angle forwards and backwards to maintain the stability of headway. When p=0, the backward-looking effect completely disappears, and the model is simplified to only consider the leading car.[28]

    3. Linear stability analysis

    The stability of car-following model is one of its important characteristics. The propagation law of density waves in a traffic flow can be obtained by analyzing the stability of carfollowing model, which has great significance for analyzing the formation mechanism of traffic congestion and proposing related improvement measures.

    Suppose that the initial operating state of vehicle is stable,i.e., the headway is h (h=L/N, L is the total length of road and N is the total number of vehicles), and all the cars maintain the same velocity,the optimized velocity(1?p)VF(θ0F)+pVB(θ0B). Then,the initial position of each vehicle is

    At this time, a small disturbance yn(t) is applied in the traffic flow so that the vehicle position changes from the steady state,written as

    And the headway can be expressed as follows:

    By substituting Eq. (10) into Eqs. (1) and (3), the front and rear visual angles can be calculated as follows:

    Because yn(t)is very small,its quadratic term is ignored.The rate of change of visual angle can be obtained by the derivation of θnF(t)and θnB(t):

    By substituting the above two sets of expressions into Eq.(6)and linearizing the optimized velocity function by Taylor expansion

    Thus,

    The small disturbance yn(t) can be expressed with Fourier mode as follows:

    The first and second derivatives can be calculated and substituted into Eq.(14)

    z can be expanded as follows:

    By substituting z into Eq.(19),we have

    Let z2=0,and the critical condition of linear stability of the model is

    When z2>0, the traffic flow will gradually return to a stable state under the influence of small disturbance, so the condition that makes the model in a linear stable state is

    Figure 2 shows the stability critical curve of the model under different conditions, presenting the relationship between sensitivity coefficient and headway. The four solid lines in Fig.2(a)from the top to the bottom are the critical curves obtained by linear analysis with the conditions λ =20, 40, 60,80,and the values of p are all 0.2. In Fig.2(b),the four solid lines from the top to the bottom are the critical curves obtained by linear analysis with the conditions p=0.2,0.1,0.05,0,and all the λ values are 40. The whole plot is divided into upper and lower regions by the curve. The stable region is above the curve,and its characteristic is that the traffic flow can run stably without obstruction. The unstable area is under the curve,and its characteristic is that a stop-and-go phenomenon will appear in the traffic stream with the gradual expansion of small disturbance.

    Each curve in Fig.2 has a vertex,which is a critical point.Multiple curves indicate that the critical point is different under different parameter conditions. For the critical sensitivity coefficient αc, when α >αc, i.e., the driver’s sensitivity is high,and the sensitivity coefficient is greater than the critical value. Thus, the traffic flow can be smoothly maintained regardless of the headway. If α <αc, the state of traffic flow mainly depends on the headway. Taking the critical headway as a reference,increase or decrease of the headway will reduce the requirement for the sensitivity coefficient, and the traffic flow will more easily enter a stable state. When the headway is near the critical value hc,the sensitivity coefficient must be increased to αcto ensure the stability of traffic flow.

    Figure 2(a)shows that with the increase of λ,the movement direction of critical point is generally downward and slightly to the right, indicating that the critical headway is slightly increased,the critical sensitivity coefficient is reduced to a large extent, and the stable area in the figure becomes larger, i.e., the driver’s attention is more concentrated. When the driver is more sensitive to changes in visual angle, it becomes easier for the traffic flow to reach a stable state. As the value of p decreases in Fig.2(b),the critical point moves to the upper left, the headway decreases, and the critical sensitivity coefficient increases, indicating that the degree of consideration of rear-view effect is gradually reduced. This will make it more difficult for traffic flow to run in a stable state at a high density.

    Fig.2. Linear critical curves of different parameters: different λ (a) and different p(b).

    Figure 3 compares the linear stability of critical curves of different models. Except for the FVD model,[8]all other models consider the influence of backward-looking effect.The stability curves of BLVD model[21]and BL&OVC model[31]almost coincide when the critical headway is consistent, indicating that these models did not substantially improve the stability conditions. When the headway is small, the stability region formed by the stability curve of model proposed in this study is relatively large,indicating that the driver is more concerned about the behavior of rear car in a congested state and uses its behavior as a basis to adjust the speed to avoid collisions.

    Fig.3. Critical curves of different models.

    4. Nonlinear stability analysis

    Nonlinear analysis is widely used in the field of nonlinear systems. In order to analyze the nonlinear characteristic of the new model describing the co-effect of forward and backward visual angle, we research on its gradually changing behavior for long waves of time and space variables around the critical point(αc,hc)using the mKdV equation.

    To do a relatively simple calculation,we use the VF[xn(t)]and VB[xn(t)] to take the place of VF[θn(t)] and VB[θn(t)] respectively and equation(6)can be converted to the following equation:

    then we can gain the expression of the(n+1)-th car further

    Subtract Eq.(24)from Eq.(25)and we can derive

    The condition α =αcis satisfied near the critical point and we define the slow variables X and T as follows using the space variable n and the time variable t in the range 0 <ε ≤1.

    Thereinto, b is a constant coefficient undetermined and we assume that the headway between two cars is

    Expanding Eq. (26) in Taylor’s series and taking to the fifth term of ε,we can get the following equation:

    The expressions of parameters mi(i=1,2,3,4,5,6,7)are

    where

    Here P and Q are constants which are expressed as P=?wn+1/(hc?ln+1)2and Q=?3wn+1/(hc?ln+1)4.

    Thereinto,the expression of gi(i=1,2,3,4,5)can be written as follows by calculation:

    As we can see in the following equation, the variables are supposed to be converted to another form and the purpose of the transformation is to derive the modified KdV equation which is with infinitesimal of higher orders.

    Then we deduce the following regularized equation

    which includes the correction term O(ε). What we do next is ignoring the term O(ε)to get further the mKdV equation with the kink-antikink soliton solution

    Then we assume and reconsider the O(ε)correction term. As a result, the propagation speed c can be composed of the parameters gi

    Therefore, the kink-antikink soliton solution of the headway is expressed in the following form

    the amplitude of the solution is equal to

    The expressions Δx=hc+A and Δx=hc?A represent the headways of two opposite phases, the low-density free flow phase and the high-density congested phase. The two mentioned phases can be combined into one coexisting phase which can be indicated by the kink-antikink soliton solution.Consequently,we achieve our goal of depicting the traffic flow state around the critical point(αc,hc)adopting the solution of the mKdV equation.

    In Fig.4, the black curves are the linear stability critical curves and the red curves with the same relative line style are the nonlinear stability critical curves. Obviously,the addition of nonlinear stability critical curves under the same conditions divides the whole picture into three regions, namely, the stable region(the common region above the two curves under the same conditions), the metastable region (the region between the two curves) and the unstable region (the common region below the two curves). Since the driver’s behavior is asymmetrical in the process of acceleration and deceleration, the nonlinear stability curves show the characteristics of asymmetry which is similar to the linear stability curves.

    Fig.4. Linear critical curves of different parameters: different λ (a) and different p(b).

    5. Numerical simulation

    To further illustrate that the introduction of rear-view effect can more accurately describe the real driving behavior and its spatial and temporal evolution, the model characteristics were analyzed in depth by numerical simulation. In the simulation,all vehicles have the same size: The length is l=5 m,and the width is w=1.8 m. A periodic boundary condition is set as the length of ring road L=1500,and the total number of vehicles N =100. The initial state is stable, i.e., the vehicles are evenly distributed,and the running velocity is the same. In the initial state, the initial position of the first car is 1, indicating that a small disturbance is applied to the first car. The initial state of the vehicles can be expressed as follows:

    The running state parameters of the vehicles can be expressed as follows:

    This section aims to verify the effect of different factors by setting different experimental conditions to prove that the proposed model can further enhance the stability of traffic flow by considering the forward and backward vision effects.

    Figure 5 shows the position distribution of each car under different values of p for a period after the disturbance is applied. The period is selected as t=2000 s-2200 s; p is 0,0.1,0.2, and 0.3; λ is 5, and other parameters remain the same.The aim is to directly analyze the influence of rate of change of visual angle in the model. In the three pictures shown in Figs. 5(a)-5(c), the uneven distribution of vehicles can be clearly observed,showing that the vehicles on the ring road are sparsely or densely distributed on different road sections,i.e.,the stop-and-go phenomena. Beyond that, with the increase of p,the driver’s visual perception of rear car becomes strong,and the degree of sparse or dense lines of the vehicle gradually slows down. The width and number of stripes that can describe the traffic congestion decrease,i.e.,the region of congestion becomes smaller,and the stop-and-go phenomenon is alleviated. When p=0.3, i.e., the stability condition is satisfied, the disturbance of traffic flow disappears after 2000 s,the headway between vehicles remains equal and unchanged,and the traffic flow enters a stable state. This shows that it can effectively enhance the stability of traffic flow for the driver to adjust the observation ratio of forward and backward information to a certain extent.

    Because the model proposed in this study can describe the effect of different car sizes on the driver’s rear-view effect, it was considered to add large cars and set different proportions in the simulation process. The size of a large car was set to 5-m long and 2.5-m wide. In the simulation,the proportion of large vehicles was set to 0.2, 0.1, 0.05, and 0.02; λ is 10; p is 0.2; other parameters were kept the same. Figure 6 shows the change trend of vehicle headway under different proportions of large vehicles after a disturbance is applied. The time points of the two figures are t=300 s and t=2000 s. Figure 6 shows that with the increase in the proportion of large vehicles, the amplitude of up and down fluctuations of headway decreases, and the stop-and-go phenomenon gradually slows down. This indicates that the stability of traffic flow can be improved if there is a certain proportion of large vehicles on the road, whether the large vehicle is observed as the leading car or the rear car because large cars can increase the driver’s attention and sensitivity to changes in their behavior.

    Fig.5. Spatiotemporal vehicle trajectory diagrams for different values of p: (a) p=0,(b) p=0.1,(c) p=0.2,and(d) p=0.3.

    Fig.6. Stop-and-go waves of different truck ratios: (a) t = 300 s, (b)t=2000 s.

    Because there are certain differences in the driving habits of drivers,there will be different decisions on whether to consider the backward-looking effect. Figure 7 shows the spatial and temporal evolution of a headway within a period after the disturbance is applied under different vehicle proportions considering the backward-looking effect. In the simulation, the proportion of vehicles considering the backward-looking effect was set to 0, 0.5, 0.8, and 0.95; the p value is 0.25; the λ value is 5; other parameters were also consistent. All the curves in the three graphs of Figs.7(a)-7(c)fluctuate significantly,indicating that the traffic flow is unstable. In Fig.7(a),all vehicles are only influenced by the front visual angle. The headway greatly fluctuates in the traffic flow. The vehicle is in a stop-and-go state, and the traffic flow is congested. In Fig.7(b),the number of vehicles that only considers the front visual angle and the number of the vehicles that considers both the front and rear visual information make up half of the total number of vehicles. There is still a large fluctuation in the headway. The distribution of headway is discontinuous and varies between two sets of similar data.Such a ratio makes the situation of traffic flow more complex.In Fig.7(c),the ratio of vehicles considering front and rear visual information simultaneously is 80%. The distance between the heads still fluctuates,but the amplitude is significantly reduced compared to the previous two cases. An uneven distribution phenomenon still exists, but it is obviously alleviated. The traffic flow tends to be more stable. In Fig.7(d),95%of the vehicles consider both the front and rear visual information. The head-to-head distance remains unchanged. The small fluctuations are caused by a small number of vehicles that only observe the vehicle in front. The traffic flow is basically stable. Based on the above analysis, the influence of driver’s perception of the rear-view effect on car-following behavior can promote the stability of the traffic flow.This also provides a basis for building an intelligent vehicle control model with the traffic flow stabilization.

    Fig.7. Spatiotemporal evolution of distance headway considering different proportions of visual information forward and backward: (a)rate=0,(b)rate=0.5,(c)rate=0.8,(d)rate=0.95.

    6. Conclusion

    This study fully considers the characteristics of driver’s corresponding car-following behavior by visually acquiring the front and rear vehicle information during driving, introduces a concept of visual angle to describe the visual stimuli of front and rear vehicles,and proposes a car-following model under the influence of rear-view.Through linear stability analysis, a critical stability condition of the model was obtained.The stable area of traffic flow can be expanded by increasing the driver’s sensitivity to the rate of change of visual angle and increasing the ratio of considering the backward-looking effect. The result of nonlinear stability analysis is consistent with that of linear stability analysis. Using numerical simulation to further analyze the characteristics of the model, the results show that different parameter settings, different sizes of large vehicles,and the proportion of backward-looking effects have a significant effect on the spatiotemporal evolution of traffic flow. Therefore, consideration of the visual characteristics of human driver in the model can be closer to the real situation. The driver observes the traffic conditions forwards and backwards with a reasonable allocation ratio and then makes driving decisions,conducive to the stability of traffic flow. This study also provides a basis for studying vehicle control rules in an intelligent networked environment.

    欧美在线黄色| 亚洲av成人精品一区久久| 日本熟妇午夜| 亚洲精品成人久久久久久| 中文字幕熟女人妻在线| 久99久视频精品免费| 久久精品国产亚洲av天美| 日本熟妇午夜| 久久这里只有精品中国| 欧洲精品卡2卡3卡4卡5卡区| 免费电影在线观看免费观看| 久久久久久久久大av| 99久久精品热视频| 精品乱码久久久久久99久播| 亚洲国产欧洲综合997久久,| 久久久久久大精品| 黄色视频,在线免费观看| 日本a在线网址| 97热精品久久久久久| 成人性生交大片免费视频hd| 午夜激情欧美在线| 老司机福利观看| 又黄又爽又刺激的免费视频.| 免费在线观看日本一区| 欧美乱色亚洲激情| 不卡一级毛片| 琪琪午夜伦伦电影理论片6080| 精品欧美国产一区二区三| 日韩欧美国产在线观看| 欧美三级亚洲精品| 哪里可以看免费的av片| 国产欧美日韩一区二区精品| 亚洲中文日韩欧美视频| 国产三级在线视频| 欧美zozozo另类| 国产精品亚洲av一区麻豆| 免费av观看视频| 日本免费一区二区三区高清不卡| 偷拍熟女少妇极品色| 久久99热6这里只有精品| 人人妻人人看人人澡| 婷婷精品国产亚洲av| 少妇熟女aⅴ在线视频| 精品人妻熟女av久视频| x7x7x7水蜜桃| 欧美国产日韩亚洲一区| 淫妇啪啪啪对白视频| 日韩大尺度精品在线看网址| 亚洲人成网站高清观看| 国产精品野战在线观看| 国产探花极品一区二区| 全区人妻精品视频| 亚洲真实伦在线观看| 一进一出好大好爽视频| 亚洲av第一区精品v没综合| 日韩高清综合在线| 国产黄a三级三级三级人| 男人舔奶头视频| 神马国产精品三级电影在线观看| 日韩成人在线观看一区二区三区| 亚洲经典国产精华液单 | 久久99热这里只有精品18| 噜噜噜噜噜久久久久久91| 亚洲精华国产精华精| 欧美xxxx性猛交bbbb| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产精品综合久久久久久久免费| www日本黄色视频网| 色哟哟·www| 国产精品免费一区二区三区在线| 色哟哟·www| 欧美日韩福利视频一区二区| 别揉我奶头 嗯啊视频| 国产aⅴ精品一区二区三区波| 色综合婷婷激情| 大型黄色视频在线免费观看| 婷婷色综合大香蕉| 国产一区二区亚洲精品在线观看| 在线播放无遮挡| 亚洲五月天丁香| 亚洲美女视频黄频| 岛国在线免费视频观看| 国产精品99久久久久久久久| 国产色婷婷99| 午夜亚洲福利在线播放| 欧美日本视频| 国产乱人视频| 波多野结衣高清作品| 好男人在线观看高清免费视频| 在线国产一区二区在线| 亚洲第一欧美日韩一区二区三区| 久久99热这里只有精品18| 成人美女网站在线观看视频| 国产蜜桃级精品一区二区三区| 噜噜噜噜噜久久久久久91| 人人妻人人看人人澡| 给我免费播放毛片高清在线观看| av中文乱码字幕在线| 精品国内亚洲2022精品成人| 欧美丝袜亚洲另类 | 国产伦精品一区二区三区四那| 国产探花在线观看一区二区| 亚洲无线在线观看| 中文字幕熟女人妻在线| 婷婷亚洲欧美| 欧美乱色亚洲激情| 国内精品久久久久久久电影| 日韩有码中文字幕| 国产午夜福利久久久久久| 久久热精品热| 亚洲成人中文字幕在线播放| 亚洲熟妇熟女久久| 欧美性感艳星| 高清毛片免费观看视频网站| 亚洲一区高清亚洲精品| 亚洲天堂国产精品一区在线| 欧美+日韩+精品| 欧美日韩黄片免| 少妇被粗大猛烈的视频| 国产欧美日韩精品亚洲av| 性色av乱码一区二区三区2| 国产欧美日韩一区二区精品| 成人亚洲精品av一区二区| 啦啦啦韩国在线观看视频| 国语自产精品视频在线第100页| 欧美三级亚洲精品| 97人妻精品一区二区三区麻豆| 国产高清视频在线观看网站| 免费看美女性在线毛片视频| 美女cb高潮喷水在线观看| 亚洲一区高清亚洲精品| 国产大屁股一区二区在线视频| 亚洲,欧美,日韩| 每晚都被弄得嗷嗷叫到高潮| 午夜免费成人在线视频| 亚洲无线在线观看| 欧美丝袜亚洲另类 | 老司机深夜福利视频在线观看| 亚洲真实伦在线观看| 久久久精品欧美日韩精品| 一级作爱视频免费观看| 老女人水多毛片| 动漫黄色视频在线观看| 波多野结衣高清作品| 一进一出抽搐动态| 欧美高清成人免费视频www| 99国产极品粉嫩在线观看| 亚洲狠狠婷婷综合久久图片| 国产视频内射| 一个人看视频在线观看www免费| 免费看a级黄色片| 51午夜福利影视在线观看| 成年女人毛片免费观看观看9| 人妻制服诱惑在线中文字幕| 亚洲片人在线观看| 亚洲五月婷婷丁香| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 男女下面进入的视频免费午夜| 亚洲aⅴ乱码一区二区在线播放| 久久久国产成人精品二区| 久久久久免费精品人妻一区二区| 日韩欧美 国产精品| 亚洲综合色惰| 亚洲av熟女| 亚洲第一欧美日韩一区二区三区| 激情在线观看视频在线高清| 丰满的人妻完整版| 高清毛片免费观看视频网站| 欧美潮喷喷水| 波多野结衣高清无吗| 亚洲七黄色美女视频| 欧美高清性xxxxhd video| 精品不卡国产一区二区三区| 一二三四社区在线视频社区8| 美女 人体艺术 gogo| 欧美+日韩+精品| 男女做爰动态图高潮gif福利片| 在线观看午夜福利视频| 精品国内亚洲2022精品成人| 国产av在哪里看| 最近中文字幕高清免费大全6 | 国产v大片淫在线免费观看| 免费av不卡在线播放| 天堂√8在线中文| 两人在一起打扑克的视频| 伦理电影大哥的女人| 香蕉av资源在线| 日韩亚洲欧美综合| 午夜日韩欧美国产| 国产男靠女视频免费网站| 搡老妇女老女人老熟妇| 亚洲国产精品久久男人天堂| 动漫黄色视频在线观看| 亚洲无线在线观看| 成熟少妇高潮喷水视频| 亚洲专区国产一区二区| 在线免费观看的www视频| 久9热在线精品视频| 18+在线观看网站| 在线观看一区二区三区| 毛片女人毛片| 欧美潮喷喷水| 欧美又色又爽又黄视频| a级毛片免费高清观看在线播放| 日韩有码中文字幕| 看黄色毛片网站| 伊人久久精品亚洲午夜| 欧美zozozo另类| 国内精品美女久久久久久| 国产白丝娇喘喷水9色精品| 日本 欧美在线| 男女下面进入的视频免费午夜| 此物有八面人人有两片| 国语自产精品视频在线第100页| 色播亚洲综合网| 日韩成人在线观看一区二区三区| 日本黄色视频三级网站网址| 老司机福利观看| 嫩草影院入口| 级片在线观看| 丝袜美腿在线中文| a在线观看视频网站| 舔av片在线| 一级a爱片免费观看的视频| 国产不卡一卡二| 日韩av在线大香蕉| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 两性午夜刺激爽爽歪歪视频在线观看| 成人国产一区最新在线观看| 国产精品三级大全| 免费在线观看亚洲国产| 看免费av毛片| 十八禁国产超污无遮挡网站| 国产成人影院久久av| a级毛片免费高清观看在线播放| 亚洲中文字幕一区二区三区有码在线看| 99热这里只有是精品50| 男人和女人高潮做爰伦理| 午夜激情福利司机影院| 91字幕亚洲| 国产亚洲欧美在线一区二区| 成人性生交大片免费视频hd| 国产激情偷乱视频一区二区| 亚洲专区中文字幕在线| 精品久久久久久久久av| 少妇人妻精品综合一区二区 | 免费人成在线观看视频色| 久久精品91蜜桃| 亚洲精品色激情综合| 午夜精品在线福利| 小蜜桃在线观看免费完整版高清| 成人国产综合亚洲| 国产午夜精品久久久久久一区二区三区 | 欧美成人a在线观看| 看片在线看免费视频| 久久精品国产自在天天线| a在线观看视频网站| 国产亚洲精品综合一区在线观看| 搡老岳熟女国产| 国产精华一区二区三区| 搡老熟女国产l中国老女人| 色吧在线观看| 99久久精品国产亚洲精品| av天堂中文字幕网| 亚洲av日韩精品久久久久久密| 直男gayav资源| 国产精品久久久久久亚洲av鲁大| 老鸭窝网址在线观看| 精品久久久久久久久久免费视频| 无遮挡黄片免费观看| 桃红色精品国产亚洲av| 内射极品少妇av片p| АⅤ资源中文在线天堂| 亚洲男人的天堂狠狠| 欧美+亚洲+日韩+国产| 一个人观看的视频www高清免费观看| 国产精品电影一区二区三区| 色噜噜av男人的天堂激情| 午夜精品久久久久久毛片777| 欧美不卡视频在线免费观看| 有码 亚洲区| 中文字幕免费在线视频6| 欧美激情久久久久久爽电影| 国产精品三级大全| 午夜福利成人在线免费观看| 精品一区二区三区视频在线观看免费| 深爱激情五月婷婷| 国产综合懂色| 自拍偷自拍亚洲精品老妇| 日韩欧美精品免费久久 | 欧美日韩中文字幕国产精品一区二区三区| 欧美成人免费av一区二区三区| 变态另类成人亚洲欧美熟女| 国产国拍精品亚洲av在线观看| 久久久久性生活片| 国产精品久久久久久久电影| 精品免费久久久久久久清纯| 国产伦精品一区二区三区四那| 久久久成人免费电影| 美女 人体艺术 gogo| 在线天堂最新版资源| 国产伦人伦偷精品视频| 国内精品久久久久精免费| 少妇丰满av| 午夜激情福利司机影院| 久久精品91蜜桃| 欧美在线一区亚洲| 日韩欧美在线乱码| 久久精品国产清高在天天线| 两人在一起打扑克的视频| 国产精品一区二区性色av| 丝袜美腿在线中文| 欧美日本视频| 99在线人妻在线中文字幕| 高清在线国产一区| 午夜福利18| 亚洲成人久久爱视频| 亚洲成人免费电影在线观看| 一个人免费在线观看电影| 怎么达到女性高潮| 国产免费av片在线观看野外av| 女人被狂操c到高潮| 亚洲专区中文字幕在线| 亚洲五月婷婷丁香| 黄色配什么色好看| 亚洲狠狠婷婷综合久久图片| 看免费av毛片| 久久亚洲精品不卡| 亚洲成人久久爱视频| 久久香蕉精品热| 丝袜美腿在线中文| 国产伦在线观看视频一区| 久久精品影院6| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品成人久久久久久| 熟女人妻精品中文字幕| 看十八女毛片水多多多| 女生性感内裤真人,穿戴方法视频| 国产精品综合久久久久久久免费| 午夜免费成人在线视频| 欧美成人性av电影在线观看| 一卡2卡三卡四卡精品乱码亚洲| 欧美三级亚洲精品| 宅男免费午夜| 精品无人区乱码1区二区| av天堂中文字幕网| 国产精品一及| а√天堂www在线а√下载| 大型黄色视频在线免费观看| 精品人妻视频免费看| 高潮久久久久久久久久久不卡| 日本熟妇午夜| 国产毛片a区久久久久| 少妇熟女aⅴ在线视频| 成人国产综合亚洲| 91av网一区二区| www.999成人在线观看| 免费搜索国产男女视频| 国内精品久久久久久久电影| 亚洲五月天丁香| 很黄的视频免费| 欧美日韩乱码在线| 亚洲欧美激情综合另类| 老女人水多毛片| 亚洲国产精品sss在线观看| 午夜精品在线福利| 男人和女人高潮做爰伦理| 中文字幕av在线有码专区| 成年版毛片免费区| www日本黄色视频网| x7x7x7水蜜桃| 成人毛片a级毛片在线播放| 一区二区三区免费毛片| 国产伦在线观看视频一区| 露出奶头的视频| 国产精品伦人一区二区| 熟女电影av网| 免费高清视频大片| 中文字幕免费在线视频6| 国产av在哪里看| 日本一二三区视频观看| 国产精品亚洲美女久久久| 欧美bdsm另类| 最新在线观看一区二区三区| 国产高潮美女av| 亚洲精华国产精华精| 日本熟妇午夜| 99热只有精品国产| 免费在线观看成人毛片| 亚洲经典国产精华液单 | 精品欧美国产一区二区三| 国产色爽女视频免费观看| 中文字幕人成人乱码亚洲影| 欧美色欧美亚洲另类二区| 永久网站在线| 男女之事视频高清在线观看| 蜜桃亚洲精品一区二区三区| 在线国产一区二区在线| 亚洲激情在线av| 免费黄网站久久成人精品 | 校园春色视频在线观看| 国产爱豆传媒在线观看| 久久久久九九精品影院| 亚洲最大成人中文| 天堂√8在线中文| 亚洲不卡免费看| 亚洲经典国产精华液单 | 日韩国内少妇激情av| 久久精品影院6| 亚洲美女搞黄在线观看 | 亚洲欧美日韩东京热| av在线天堂中文字幕| 在线国产一区二区在线| 久久久久久久午夜电影| 国产精品野战在线观看| 麻豆国产97在线/欧美| 日日夜夜操网爽| 国产成人福利小说| 免费观看精品视频网站| 亚洲av成人av| 99热精品在线国产| 99在线人妻在线中文字幕| 欧美最黄视频在线播放免费| 亚洲av美国av| 精品久久久久久久久久免费视频| 久久久久久久久大av| 色尼玛亚洲综合影院| 搡女人真爽免费视频火全软件 | 在线观看一区二区三区| 欧美色欧美亚洲另类二区| 国产精品综合久久久久久久免费| 最近最新免费中文字幕在线| 天天一区二区日本电影三级| 黄色视频,在线免费观看| 天堂av国产一区二区熟女人妻| 久久精品影院6| av专区在线播放| 天堂√8在线中文| 日日夜夜操网爽| 久久精品综合一区二区三区| 免费在线观看日本一区| 哪里可以看免费的av片| 美女被艹到高潮喷水动态| 国产一区二区激情短视频| 精品国内亚洲2022精品成人| 国产av不卡久久| 欧美潮喷喷水| a级毛片免费高清观看在线播放| 国产精品野战在线观看| 久久久久久久午夜电影| 免费av观看视频| 免费大片18禁| 91狼人影院| 99国产综合亚洲精品| 国产伦一二天堂av在线观看| 观看美女的网站| 国产精品乱码一区二三区的特点| 人人妻人人澡欧美一区二区| 美女cb高潮喷水在线观看| 国产白丝娇喘喷水9色精品| 欧美性感艳星| 丁香六月欧美| 黄片小视频在线播放| 欧美日韩亚洲国产一区二区在线观看| 久久欧美精品欧美久久欧美| 午夜精品在线福利| .国产精品久久| 高清日韩中文字幕在线| 色吧在线观看| 久久99热6这里只有精品| 亚洲中文字幕一区二区三区有码在线看| 狠狠狠狠99中文字幕| 亚洲一区高清亚洲精品| 黄色日韩在线| 丝袜美腿在线中文| 国产精品98久久久久久宅男小说| 欧美黑人巨大hd| 亚洲人成网站在线播放欧美日韩| 久久性视频一级片| av在线蜜桃| 中文字幕精品亚洲无线码一区| 午夜福利在线在线| 亚洲精品456在线播放app | 久久国产乱子伦精品免费另类| 亚洲精品456在线播放app | 国产一区二区三区视频了| 色精品久久人妻99蜜桃| 亚洲午夜理论影院| 国产欧美日韩精品亚洲av| 日本一本二区三区精品| 一区福利在线观看| 久久99热6这里只有精品| 人人妻,人人澡人人爽秒播| 精品无人区乱码1区二区| 亚洲精品亚洲一区二区| 最近在线观看免费完整版| 不卡一级毛片| 在现免费观看毛片| 日韩人妻高清精品专区| 国产真实乱freesex| 国产午夜福利久久久久久| 日韩欧美精品免费久久 | 日本免费a在线| 国产亚洲精品综合一区在线观看| 亚洲av电影在线进入| 99精品在免费线老司机午夜| 天堂√8在线中文| 天堂网av新在线| 黄片小视频在线播放| 午夜久久久久精精品| 久久精品久久久久久噜噜老黄 | а√天堂www在线а√下载| 一本综合久久免费| 免费看美女性在线毛片视频| 国产精品亚洲美女久久久| 国产一区二区三区在线臀色熟女| 色综合婷婷激情| 我要看日韩黄色一级片| 亚洲中文日韩欧美视频| 麻豆国产97在线/欧美| 免费一级毛片在线播放高清视频| 国产精品三级大全| 成人精品一区二区免费| 成年女人永久免费观看视频| 啦啦啦韩国在线观看视频| 国产av不卡久久| 国产精品电影一区二区三区| 有码 亚洲区| 赤兔流量卡办理| 露出奶头的视频| 99久久99久久久精品蜜桃| 国产一区二区亚洲精品在线观看| 午夜精品一区二区三区免费看| 十八禁网站免费在线| 国产色爽女视频免费观看| 美女高潮喷水抽搐中文字幕| 人妻丰满熟妇av一区二区三区| 亚洲av免费在线观看| 国产精品一区二区免费欧美| 亚州av有码| 99热这里只有是精品50| 天堂√8在线中文| 99久久精品一区二区三区| 琪琪午夜伦伦电影理论片6080| 成人精品一区二区免费| 国产精品爽爽va在线观看网站| 搡老熟女国产l中国老女人| 国产精品女同一区二区软件 | 看十八女毛片水多多多| 99久久无色码亚洲精品果冻| 亚洲经典国产精华液单 | 亚洲成人免费电影在线观看| 禁无遮挡网站| 国产精品综合久久久久久久免费| 国产成人av教育| 人人妻人人澡欧美一区二区| 桃色一区二区三区在线观看| 日本与韩国留学比较| 国产精品伦人一区二区| 男女下面进入的视频免费午夜| 高清在线国产一区| 欧美黑人欧美精品刺激| 中文字幕久久专区| 国产av不卡久久| 午夜日韩欧美国产| 日本免费a在线| 成年版毛片免费区| 日韩高清综合在线| 两个人的视频大全免费| 18禁黄网站禁片午夜丰满| 色5月婷婷丁香| 亚洲第一区二区三区不卡| 他把我摸到了高潮在线观看| 性插视频无遮挡在线免费观看| 又黄又爽又免费观看的视频| av在线观看视频网站免费| 久久久国产成人免费| 日韩精品中文字幕看吧| 琪琪午夜伦伦电影理论片6080| 久久久久国内视频| 国产伦在线观看视频一区| 日韩欧美精品v在线| 最近在线观看免费完整版| 亚洲三级黄色毛片| 国产免费av片在线观看野外av| 日韩高清综合在线| 老熟妇乱子伦视频在线观看| 国产午夜精品久久久久久一区二区三区 | 精品久久久久久久人妻蜜臀av| 乱人视频在线观看| 亚洲不卡免费看| 亚洲乱码一区二区免费版| 制服丝袜大香蕉在线| 亚洲av中文字字幕乱码综合| av欧美777| 99热这里只有是精品在线观看 | 91av网一区二区| 最新在线观看一区二区三区| 国产伦在线观看视频一区| 国产精品日韩av在线免费观看| 国产成人av教育| 神马国产精品三级电影在线观看| 真人一进一出gif抽搐免费| 久久久久久久亚洲中文字幕 | 免费人成在线观看视频色| 国产精品爽爽va在线观看网站| 国产精品乱码一区二三区的特点| 性色avwww在线观看| 老熟妇仑乱视频hdxx| 最新在线观看一区二区三区| 久久精品国产自在天天线| 国产大屁股一区二区在线视频| 精品久久久久久久久av| 亚洲午夜理论影院| 女人十人毛片免费观看3o分钟|