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

    Optimal allocation model of signalized intersectionswith dynamic use of exit lanes for left turn

    2020-07-20 08:42:40QuZhaoweiQiXingzuChenYonghengBaiQiaowenTaoChuqingLiuFanghong

    Qu Zhaowei Qi Xingzu Chen Yongheng Bai Qiaowen Tao Chuqing Liu Fanghong

    (School of Transportation, Jilin University, Changchun 130022, China)

    Abstract:In order to improve the operational efficiency of heavy left-turn demand intersections, an optimal allocation model of an intersection with dynamic use of exit lanes for left turns (EFL) is proposed. The constraints of setting EFL are analyzed, including the number and length of reverse variable lanes, flow direction constraints, and signal constraints, etc. The constraints and control variables are combined in a unified framework for simultaneous optimization. The objective functions are defined as the average delay and left-turn capacity of an intersection. The model is solved by a non-dominated genetic algorithm (NSGA-Ⅱ). The results show that after the optimal allocation of EFL, the average vehicle delays of the intersection can be reduced by 14.9% and left-turn capacity can be increased by 19.3%. The effectiveness of the optimal allocation model of EFL is demonstrated.

    Key words:traffic engineering; exit lane for left turn; signaled intersection; multi-objective optimization

    Owing to heavy traffic demand, it is difficult to solve the congestion problem of intersections by optimizing signal timing or adding new lanes. Therefore, researchers have been looking for unconventional management methods to increase capacities of signalized intersections. Such designs included superstreet intersections[1-2], displaced left-turn intersections[3-4], roundabouts[5], and others. In these unconventional intersection designs, the exit lanes for left turn (EFL) take full advantages of the intersection capacity, especially under the heavy left-turn traffic. The unique characteristic of the EFL intersection is its mixed-used area. The mixed-used area is used as exit-lane of opposing through vehicles or left-turn lane during different periods of a signal cycle. Also, there is a pre-signal and a median opening at the upstream of the EFL intersection. Therefore, left-turn vehicles can drive into the mixed-used area at the pre-signal and wait for the green light of the main-signal to finish their turning movements.

    In the early development of lane assignment, conventional design methods of lane marking only took all the approach lanes of the signalized intersection into account. Each lane has a specific movement apart from exit lanes. In addition, the design of lane marking should consider the signal timings calculation. The conventional lane assignment considered the lane marking factor first, and then the signal setting is determined based on the lane marking. However, it may not always produce a truly optimal plan of lane markings for the EFL intersection.

    Recently, a lane-based optimization model combining the design of lane marking and signal settings for isolated signal-controlled intersections was developed. The lane marking was always defined as binary variables into an optimization frame. The set of binary variables also include the common flow multiplier, cycle length, traffic flows, maximum reserve capacity as well as starts and durations of greens for traffic movements on approach lanes. In the signal control, Assi et al.[6]proposed a quick method to find the optimum lane group for 3-lane and 4-lane approaches at junctions. Yu et al.[7]proposed an optimal traffic signal control framework to minimize the total travel time. However, the importance of signal timings with lane-based allocation should not be ignored. The integrated optimization model further improved the effectiveness of traffic flows[8-12].

    A lane-based optimization target model was proposed by Wong et al[13-14]. The model optimized lane attributes and signal periods based on the number of approach and exit lanes. The objective function was maximizing the intersection capacity or minimizing cycle length. The multi-objective problems were solved by the binary mixed-integer linear program (BMILP). However, once lane allocations were determined, they can only be used as approach lanes or exit lanes. Zhang et al.[15]proposed a robust optimal fixed-time signal timing and considered the environmental influence factors. However, due to the change of the inlet flow at the intersection, the characteristics of lanes and the timing of signals were not given. Wong and Heydecker[16]further studied the lane optimization problem at a single-signal intersection. They defined new lane attribute variables and introduced them into a new optimization model. However, they ignored the influence of crossing pedestrian signals, especially the phenomenon of high demand for crosswalks in China. Based on the research of Wong and Heydecker[16], Yu et al.[17]improved the single-signal intersection optimization model, in which safety islands were considered. They also proposed a multi-objective lane optimization model, which was not sensitive to traffic flow fluctuations. However, there was a lack of constraints on vehicle lanes change. Li et al.[18]proposed a multi-objective optimization method for improving traffic performance at intersections. The transportation efficiency, energy consumption and road safety were considered as the objective functions in the model. However, the EFL intersection was not considered in the model.

    In order to formulate the intersection design problems into a mathematical program, this study proposes an optimization model of exit lanes for left-turn including lane schedules and signal timings. The design variables can be optimized simultaneously to obtain better intersection capacity. The objective functions are the minimization of the vehicle delay and the maximization of the left-turn capacity. Considering that the non-dominated sorting genetic algorithm-Ⅱ (NSGA-Ⅱ) has a better sorting parameter, incorporates elitism and no sharing parameter needs to be chosenapriori, we select it to solve the model. Also, a set of numerical examples is evaluated to demonstrate the effectiveness of the model.

    1 Model Formulation

    1.1 Optimization objective selection

    In order to improve the efficiency of traffic flows, this study will be conducted from the following two aspects: lanes usage and safety impacts of the EFL design. The intersection traffic capacity and intersection vehicle delay are selected as the optimization objectives.

    1.2 Input data

    In order to obtain the best optimization results, this paper established a lane assignment optimization model including the traffic intersection constraints, lane attributes, signal timing and median open location.

    1) Lane allocation variables

    Fig.1 shows a typical intersection with four arms. If a movement from armito armjis permitted on lanek, thenei,k,j=1; otherwise,ei,k,j=0.

    2) Minimum permitted movement on approach lanes

    Each lane must be an approach lane or an exit lane, and each approach lane must have at least one turning or through movement, which can be specified as

    (1)

    ?j=1,2,…,N-1;i=1,2,…,N

    whereNis the number of traffic arms.

    Fig.1 Numbering convention for the approach lane and exit lane in an intersection

    If the traffic in thek-th lane of approach laneiis prohibited to drive towards exit lanej, its corresponding lane flow should be 0.

    Mei,k,j≥qi,k,j≥0

    (2)

    ?i=1,2,…,N;k=1,2,…,Ni,j=1,2,…,N-1

    whereMis a sufficiently large positive constant andqi,k,jis the lane turning flows assigned on permitted lanes.

    For each traffic flow fromitoj, the number of exit lanes should not be less than that of approach lanes corresponding to this flow direction. It is to avoid impact on each traffic flow on the exit lanes.

    (3)

    ?j=1,2,…,N-1;i=1,2,…,N

    whereNEis the number of exit arms, andNTis the number of traffic lanes in armi.

    3) Cycle length

    The minimum and maximum values of the intersection cycle length are denoted byCminandCmax, respectively. The constraints on the cycle length can be specified as

    Cmin≤C≤Cmax

    (4)

    4) Lane signal setting

    For safety reasons, if an approach lane permits different vehicle movements, it is necessary that they receive the same signal indications, including the green duration and green start time.

    M(1-ei,k,j)≥Θi,k-θi,k≥-M(1-ei,k,j)

    (5)

    M(1-ei,k,j)≥Φi,k-φi,k≥-M(1-ei,k,j)

    (6)

    ?j=1,2,…,NT-1;k=1,2,…,Ni;i=1,2,…,NT

    whereΘi,kis the start of green on approach traffic lanes in armivia lanek;θi,kis the start of green for turning movements on lanekin armi;Φi,kis the durations of green for the approach traffic lanes on lanekin armiandφi,kis the durations of green for turning movements on lanekin armi. If the vehicle movement (i,k,j) is allowed (ei,k,j=1), the green start time and green duration in armivia lanekmust be equal to the signal group (i,j).

    5) Start of the green time setting

    Since the signal period at the intersection is cyclic, the variable of the green start can take any value among a cycle as long as it satisfies other relevant constraints. However, for the convenience of calculation, this paper sets the range of all green start time variables as (0,1), namely,

    1≥θi,j≥0

    (7)

    ?j=1,2,…,N-1;i=1,2,…,N

    6) Duration of green

    The green duration of each approach lane is determined by the minimum green duration of the flow direction. This can be specified by the following constraint set:

    1≥φi,j≥gi,j/C

    (8)

    ?j=1,2,…,NT-1;k=1,2,…,Ni;i=1,2,…,NT

    7) Order of signal displays

    For any two conflicting phase groups (i,j) and (m,n), conflict can be avoided by releasing the sequence of the phase.

    Pij,mn+Pmn,ij=1

    (9)

    ?(i,j),(m,n)∈ψs

    A group of conflicting signal phases is represented asPij,mn. If the green start time signal group (m,n) follows the signal group (i,j), the conflicting signal phasePij,mn=0; otherwisePij,mn=1. If the green start time signal group (i,j) follows the signal group (m,n), the conflicting signal phasePmn,ij=0; otherwisePmn,ij=1.

    8) Flow conservation

    The traffic demand flowsQi,jare known. Lane turning flowsqi,jhave been assigned on the approach lanes. The flow multiplierμ<1 indicates that the traffic volume of the intersection overloads 100(1-μ)%, and the flow multiplierμ>1 indicates that the intersection can also increase the traffic volume of 100(μ-1)%. The following set of constraints can be expressed as

    (10)

    whereμis the common flow multiplier andQi,jis the demand flowing from armito armj.

    9) Mixed-usage-area constraint

    This paper presents a lane-based optimization model combining the exit lanes for the left turn. The vehicles can turn left through the EFL under the pre-signal, indicating that the exit lane becomes a left-turn lane. All the left-turn vehicles that drive into the mixed usage area should be discharged during the green time of the pre-signal. Such design needs a good coordination between the main-signal at the intersection and the pre-signal at the median open.

    For the EFL lanes, the capacity analysis of left-turn vehicles can be more complex. In practice, the vehicles in the EFL lane must be discharged during the left-turn phase of the main signal. Otherwise, the left-turn vehicles will block the traffic coming from the other direction. Thus, the capacity gains of the left-turn movements at an intersection with the EFL design can be considered as the number of left-turn vehicles that enter the EFL. The location of the upstream median opening determines the maximum number of the left-turning vehicles. The basic shockwave diagram of the EFL design is shown in Fig.2.

    (11)

    whereLiis the distance between the main-signal and the pre-signal.

    Fig.2 The basic shockwave diagram of EFL design

    10) Flow factor

    The saturation flow and delay on a pair of adjacent approach lanes with a common signal group must be identical. If there are multiple lanes in the same flow direction, the flow ratio of the same flow direction lanes should be consistent as follows:

    (12)

    M(2-ei,k,j-ei,k+1,j)≥yi,k+1-yi,k≥-M(2-ei,k,j-ei,k+1,j)

    (13)

    ?j=1,2,…,NT-1;k=1,2,…,Ni;i=1,2,…,NT

    11) Average delay of the vehicle

    The calculation method of the average delay at intersections is calculated according to the Highway Capacity Manual HCM (2010)[19].

    Di,k=D1,i,k(PF)+D2,i,k+D3,i,k

    (14)

    wherePFis the adjustment parameter of average delay.

    The calculation method of all vehicle delays in lanekof approach laneican be expressed as

    (15)

    The increased delay calculation method of lanekon approach laneican be expressed as

    (16)

    whereTis the duration of the analysis period;kfis the incremental delay factor;iis the upstream filtering adjustment factor;ci,kis the capacity of lanekat approachiandxi,kis the volume to capacity ratio on lanekin armi.

    1.3 Optimization formulations

    The maximum capacity of the left turn is an important design aspect of intersections. It should consider the requirements of geometric layout and signal timings. The objective functions are defined as the maximization of the left-turn capacity (CAP)Cleftand the minimization of the average vehicle delayD.

    (17)

    (18)

    Subject to the linear constraints (1) to (16), in this nonlinear integer programming model, the decision variables are lane variables and signal variables.

    2 Model Solution Based on Genetic Algorithm

    The above optimization model is a multi-objective mixed integer nonlinear programming problem. In this model,Lis the decision factor. Although it is a multi-objective model, the objective functions can be transformed into a single-objective optimization model by setting weighting parameters. Moreover, the objective functions can be solved by NSGA-Ⅱ.

    The observed average traffic volume is denoted asQ′ and the specific solution steps are as follows:

    Step1Initial population

    Set the evolution counterg=0 and the maximum evolution algebrag=500. Within the range of the established region descriptor,M=100 individuals are randomly generated as the initial populationP(0).

    Step2Fitness function

    The size of fitness indicates that the fitness of individuals and individuals with higher fitness are more likely to inherit to the next generation. According to the characteristics of the objective function established in this paper, the fitness function is selected asf=1/F. Then the fitness value of each individual in populationP(t) is calculated. Two individuals were randomly selected from the population, and the fitness between the two was compared.

    Step3Genetic operations

    The basic genetic operations include selection, crossing and variation.

    1) The selection algorithm selects individuals with a higher fitness value according to the fitness value of each individual in the current population and passes on the selected individuals into the daughter population according to the set rules.

    2) The crossover operation is to randomly match pairs of each individual in the population. Each individual will exchange partial chromosomes with a certain probability to obtain a new generation of individuals. In this paper, the crossover probability is set to be 0.7.

    3) The mutation operation is to generate a new individual by randomly selecting the individual to change the gene value at one or some other loci to other alleles with a certain probability, and the mutation probability is set to be 0.05.

    Step4Algorithm stopping

    If the evolutionary algebragof the genetic algorithm is 500, the algorithm can be terminated by taking the individual with the maximum fitness.

    3 Numerical Studies

    In order to verify the effectiveness of the model in this paper, a four-arm junction as shown in Fig.3 is studied. From the signal timing perspective, the design of the EFL intersection needs a good coordination between the main-signal and the pre-signal as shown in Fig.4. The purpose of the studied case is to demonstrate whether the EFL will improve its efficiency or not. Relevant input data and assumptions are given as follows.

    Fig.3 Layout of a four-arm junction in the example

    For simplicity, the saturation flows for straight-ahead movements are 1 800 pcu/(h·lane). The maximum acceptable degrees of saturation are taken as 90% on all lanes. If the maximum acceptable degrees of saturation

    Fig.4 Signal phases

    is above 90%, the optimization effect of the model is not significant. According to the HCM, the minimum cycle length is 60 s and there is no limitation on the maximum cycle length. In this study, the minimum and maximum cycle lengths are assumed to be 60 and 200 s, respectively. The minimum green durations are 5 s for all traffic movements.

    The traffic demand at the intersection adopts the data during the peak hour. The left-turn phase and the straight-line phase are set to be 5 and 3 vehicles, respectively, in the last cycle. Relevant traffic demand dataQi,jis tabulated in Tab.1. The saturation degrees of all approach traffic lanes are under 0.9. NSGA-Ⅱ can solve the above nonlinear integer programming problems of the mixed usage area.

    Tab.1 Traffic demands of a four-arm junction for the example

    The calculation results shown in Fig.5 show that the EFL intersection has two left-turn lanes in the west approach lanes including the mixed usage area. The flow coefficientμ=1.241 indicates that the lane allocation method can meet the requirements of the traffic demand. In the conventional intersection, it is shown that this lane allocation method can also meet traffic demands under certain saturation conditions. However, the EFL intersection has a greater traffic capacity. Using the VISSIM simulation software, the signal optimization simulation results are shown in Tab.2 and Tab.3. It can be observed from Tabs.4 and 5 that the EFL intersection has a shorter signal cycle. The left-turn effective green time from arm 1 to arm 2 of EFL intersection is shorter than that of conventional intersection. This means that the EFL intersection can not only improve the efficiency of the left turn but also reduce intersection cycle. In the simulation, the default behavior of urban road drivers in VISSIM was adopted. All the models were the cars with the expected speed of 40 to 50 km/h, and the saturation flow rate was about 1 800 veh/h.

    Fig.5 Optimization results of lane arrangements

    Tab.2 Signal optimization simulation results of EFL intersection

    Tab.3 Signal optimization simulation results of conventional intersection

    Fig.6 Signal control process of EFL intersection

    Stage1The left-turn movements from the minor street start to be discharged (movements 5 and 10).

    Stage2The through vehicles (movements 4 and 9) from the minor street start to be discharged. After a few seconds, the pre-signal turns green and the left-turning vehicles (movements 8 and 3) are starting to advance into the exit-lane for left-turn through the median opening.

    Stage3The left-turn phases in the main signals turn green and all the left-turning vehicles (movements 7 and 1) are beginning to be discharged, including those in the contraow and the conventional left-turn lanes.

    Stage4The left-turning vehicles continue to be discharged (movements 1 and 5). The pre-signal turns red, and there are no vehicles left in the exit-lane.

    Stage5The through vehicles (movements 2 and 6) from the major street begin to be discharged, and the cycle is ended. The simulation results are shown in Fig.7 and Tab.4. We developed analytical models to calculate the capacity and average delay. As shown in Fig.7 and Tab.4, when intersection saturation is over 80%, the EFL design can improve the capacity and decrease the average vehicle delay of vehicles at the signal intersection. Moreover, the capacity gains mainly depend on the left-turning vehicles driving through the missed-usage-area. Therefore, the EFL design can significantly improve the capacity of intersections. When the intersection saturation is less than 80%, it needs to be determined whether the intersection needs to set EFL based on the model calculation results or not. The vehicle average delay of the EFL intersection is 14.9% lower than that of the conventional intersection (NEFL). The capacity of left turn is 9.3% higher than that of the conventional intersection. Therefore, in the optimization model, it is necessary to consider the mixed usage area on exit lanes. The EFL can reduce the average vehicle delay caused by obstruction and further improve the capacity of intersections.

    4 Conclusions

    1) The effectiveness of the EFL model and solution algorithms were validated through numerical studies.

    (a)

    Tab.4 Signal optimization simulation results of severe traffic demand

    The results show that both the EFL and conventional intersection can meet the traffic demand under a certain saturation condition. However, the model of EFL has a larger capacity.

    2) In addition, the simulation platform was built via the VISSIM simulation for comparison. The EFL showed significant improvements in traffic operational performance. The average vehicle delays of intersection can be reduced by 14.9% and left-turn capacity can be increased by 19.3%.

    3) Future research work should further pay attention to different traffic demands and ensure that it is adapting to the randomness of traffic demand in each time period. The proposed optimal model has a potential to increase the efficiency of the EFL design. However, more comprehensive research is needed to evaluate the effectiveness and safety impacts of the EFL design. The authors suggest that future studies should focus on these issues.

    成人亚洲欧美一区二区av| 亚洲第一av免费看| 97人妻天天添夜夜摸| 两个人看的免费小视频| 国产成人精品一,二区| 国产免费视频播放在线视频| 宅男免费午夜| 亚洲伊人色综图| 一二三四中文在线观看免费高清| 欧美精品人与动牲交sv欧美| 69精品国产乱码久久久| 亚洲少妇的诱惑av| 黑人猛操日本美女一级片| 亚洲第一青青草原| 啦啦啦视频在线资源免费观看| 日本色播在线视频| 国产成人精品一,二区| 国产在视频线精品| 又大又黄又爽视频免费| 91精品三级在线观看| 欧美 亚洲 国产 日韩一| 免费高清在线观看日韩| 亚洲精品在线美女| 色吧在线观看| 欧美日韩视频精品一区| 这个男人来自地球电影免费观看 | 女人高潮潮喷娇喘18禁视频| 99热全是精品| 亚洲国产色片| 免费观看性生交大片5| 99re6热这里在线精品视频| 九九爱精品视频在线观看| 亚洲四区av| 亚洲av福利一区| 午夜福利在线免费观看网站| 欧美日韩一区二区视频在线观看视频在线| 国产熟女欧美一区二区| 精品国产乱码久久久久久男人| 捣出白浆h1v1| 亚洲,欧美精品.| 久久久欧美国产精品| 午夜福利乱码中文字幕| 女的被弄到高潮叫床怎么办| 美女国产高潮福利片在线看| 国产黄频视频在线观看| 秋霞伦理黄片| 91精品伊人久久大香线蕉| 国产极品粉嫩免费观看在线| 免费在线观看黄色视频的| 婷婷色麻豆天堂久久| 午夜日本视频在线| 天天躁夜夜躁狠狠躁躁| 欧美精品国产亚洲| 久久久久久久久久人人人人人人| 秋霞在线观看毛片| av免费观看日本| 精品视频人人做人人爽| 9热在线视频观看99| 欧美日韩精品网址| 国产精品国产三级国产专区5o| 久热久热在线精品观看| 久久久久久人妻| 日韩制服丝袜自拍偷拍| 99久久精品国产国产毛片| 午夜免费鲁丝| 各种免费的搞黄视频| 人妻 亚洲 视频| www日本在线高清视频| av片东京热男人的天堂| 视频区图区小说| 亚洲综合色惰| 亚洲情色 制服丝袜| 一二三四中文在线观看免费高清| 日韩欧美精品免费久久| 亚洲色图综合在线观看| 亚洲精品视频女| 如何舔出高潮| 久久人人爽av亚洲精品天堂| 纵有疾风起免费观看全集完整版| 国产在线一区二区三区精| 国产黄色视频一区二区在线观看| 永久网站在线| 日韩制服丝袜自拍偷拍| 欧美精品高潮呻吟av久久| 久久亚洲国产成人精品v| 欧美+日韩+精品| 午夜福利在线观看免费完整高清在| 国产熟女欧美一区二区| 色94色欧美一区二区| 亚洲,欧美精品.| 亚洲成人手机| 欧美日韩一级在线毛片| 考比视频在线观看| 精品国产一区二区久久| 久久精品久久久久久噜噜老黄| 日韩制服骚丝袜av| 亚洲精品国产一区二区精华液| 亚洲国产av新网站| 国产av一区二区精品久久| 亚洲国产精品成人久久小说| 亚洲一级一片aⅴ在线观看| 男女无遮挡免费网站观看| 国语对白做爰xxxⅹ性视频网站| 少妇人妻 视频| 中文字幕人妻丝袜一区二区 | 日韩欧美精品免费久久| 午夜91福利影院| 欧美精品亚洲一区二区| 亚洲欧美一区二区三区国产| 2021少妇久久久久久久久久久| 人人妻人人添人人爽欧美一区卜| 多毛熟女@视频| 寂寞人妻少妇视频99o| 丝袜喷水一区| 日韩免费高清中文字幕av| 色播在线永久视频| 我要看黄色一级片免费的| 97精品久久久久久久久久精品| 国产日韩欧美亚洲二区| 精品人妻一区二区三区麻豆| 校园人妻丝袜中文字幕| 亚洲av欧美aⅴ国产| 久久久a久久爽久久v久久| 国产一区二区 视频在线| 免费观看在线日韩| 国产淫语在线视频| www.熟女人妻精品国产| 国产精品三级大全| 亚洲,一卡二卡三卡| av不卡在线播放| 美女视频免费永久观看网站| 欧美激情 高清一区二区三区| 亚洲三级黄色毛片| 午夜影院在线不卡| 亚洲美女搞黄在线观看| 欧美日韩视频高清一区二区三区二| 欧美日韩成人在线一区二区| 秋霞伦理黄片| 国产老妇伦熟女老妇高清| 妹子高潮喷水视频| 丁香六月天网| 精品99又大又爽又粗少妇毛片| 两个人看的免费小视频| 亚洲av成人精品一二三区| 人人妻人人澡人人爽人人夜夜| 黄色视频在线播放观看不卡| 亚洲成国产人片在线观看| 国产精品 国内视频| 免费在线观看完整版高清| 国产乱人偷精品视频| 亚洲第一av免费看| 热99久久久久精品小说推荐| 国产成人精品无人区| 国产精品一国产av| 这个男人来自地球电影免费观看 | 午夜福利乱码中文字幕| 亚洲一区中文字幕在线| 色视频在线一区二区三区| 两性夫妻黄色片| 99re6热这里在线精品视频| 美女脱内裤让男人舔精品视频| 丝袜人妻中文字幕| 亚洲av国产av综合av卡| 亚洲精品一二三| 建设人人有责人人尽责人人享有的| 久久人人爽av亚洲精品天堂| 在线看a的网站| 亚洲精品久久成人aⅴ小说| 黄色视频在线播放观看不卡| 精品少妇久久久久久888优播| 久久这里只有精品19| 777米奇影视久久| 久久99热这里只频精品6学生| 亚洲欧美精品综合一区二区三区 | 国产爽快片一区二区三区| 国产精品久久久久久久久免| 侵犯人妻中文字幕一二三四区| 精品卡一卡二卡四卡免费| 成人国产av品久久久| 久久精品人人爽人人爽视色| 欧美精品一区二区大全| 制服丝袜香蕉在线| 人人妻人人爽人人添夜夜欢视频| 男女午夜视频在线观看| 97精品久久久久久久久久精品| 亚洲精品国产av成人精品| 亚洲国产精品一区三区| 亚洲熟女精品中文字幕| 你懂的网址亚洲精品在线观看| 日本欧美视频一区| 亚洲伊人色综图| 宅男免费午夜| 亚洲精品自拍成人| 一区二区av电影网| 久久久久久久久免费视频了| 麻豆精品久久久久久蜜桃| 晚上一个人看的免费电影| 自线自在国产av| 国产一区有黄有色的免费视频| 少妇的逼水好多| 激情视频va一区二区三区| 超碰成人久久| xxx大片免费视频| 视频在线观看一区二区三区| 久久久久视频综合| 高清不卡的av网站| 人妻系列 视频| 丝袜喷水一区| 黄色一级大片看看| 精品人妻在线不人妻| 国产精品麻豆人妻色哟哟久久| 女的被弄到高潮叫床怎么办| 免费黄色在线免费观看| av电影中文网址| 天美传媒精品一区二区| 国产精品欧美亚洲77777| 少妇人妻久久综合中文| 久久精品亚洲av国产电影网| 日本色播在线视频| 国产男女内射视频| 老司机影院成人| 国产日韩欧美亚洲二区| 精品国产乱码久久久久久男人| 视频区图区小说| 校园人妻丝袜中文字幕| 老司机亚洲免费影院| 天天影视国产精品| 久久久久久久久久久免费av| 男女免费视频国产| 成人免费观看视频高清| 亚洲男人天堂网一区| 老女人水多毛片| 婷婷成人精品国产| 亚洲欧美日韩另类电影网站| 99久国产av精品国产电影| 国产1区2区3区精品| 国产日韩欧美视频二区| 亚洲精品国产av蜜桃| 国产精品香港三级国产av潘金莲 | 中文乱码字字幕精品一区二区三区| 男女无遮挡免费网站观看| 人人妻人人澡人人看| 国产精品久久久久成人av| 熟女av电影| 日本wwww免费看| 日本欧美国产在线视频| av有码第一页| 日韩av免费高清视频| 亚洲国产欧美日韩在线播放| 寂寞人妻少妇视频99o| 中文字幕最新亚洲高清| 国产成人精品一,二区| 婷婷色麻豆天堂久久| 有码 亚洲区| 国产乱人偷精品视频| 永久免费av网站大全| 亚洲国产精品一区三区| 午夜影院在线不卡| 国产一区二区三区av在线| 老司机亚洲免费影院| 少妇被粗大的猛进出69影院| 26uuu在线亚洲综合色| 美女国产高潮福利片在线看| 黄色配什么色好看| 又粗又硬又长又爽又黄的视频| 欧美成人午夜免费资源| 日韩人妻精品一区2区三区| 精品少妇久久久久久888优播| 一区二区日韩欧美中文字幕| 欧美日韩一级在线毛片| 超碰成人久久| 在线精品无人区一区二区三| 一级毛片 在线播放| 如何舔出高潮| 午夜福利网站1000一区二区三区| 9191精品国产免费久久| 美女中出高潮动态图| 在线观看免费日韩欧美大片| 亚洲欧美精品综合一区二区三区 | 熟女少妇亚洲综合色aaa.| 日韩中文字幕欧美一区二区 | 国产一区有黄有色的免费视频| 日本午夜av视频| 国产老妇伦熟女老妇高清| 高清视频免费观看一区二区| √禁漫天堂资源中文www| 女的被弄到高潮叫床怎么办| 高清黄色对白视频在线免费看| 精品卡一卡二卡四卡免费| 亚洲精品久久成人aⅴ小说| 十分钟在线观看高清视频www| 亚洲欧美精品综合一区二区三区 | 欧美国产精品一级二级三级| 999久久久国产精品视频| 午夜福利在线免费观看网站| 国产精品免费大片| 一个人免费看片子| 美女主播在线视频| 成年人午夜在线观看视频| 黄色一级大片看看| 色哟哟·www| 少妇精品久久久久久久| 亚洲天堂av无毛| 久久久久国产精品人妻一区二区| 久久女婷五月综合色啪小说| 黑人猛操日本美女一级片| 天天躁狠狠躁夜夜躁狠狠躁| 日韩伦理黄色片| 免费看av在线观看网站| 国产福利在线免费观看视频| 亚洲色图综合在线观看| 久久鲁丝午夜福利片| 国产精品人妻久久久影院| 乱人伦中国视频| 1024视频免费在线观看| 欧美日韩精品网址| av.在线天堂| 欧美国产精品一级二级三级| 久久精品aⅴ一区二区三区四区 | 少妇人妻精品综合一区二区| 色网站视频免费| 91精品国产国语对白视频| 日韩精品免费视频一区二区三区| 日本wwww免费看| 国产乱来视频区| 熟女电影av网| 高清不卡的av网站| 日本91视频免费播放| 亚洲美女视频黄频| 啦啦啦啦在线视频资源| 日产精品乱码卡一卡2卡三| 性色avwww在线观看| 赤兔流量卡办理| 一本—道久久a久久精品蜜桃钙片| 激情视频va一区二区三区| 国产高清不卡午夜福利| 亚洲成国产人片在线观看| 国产1区2区3区精品| 男的添女的下面高潮视频| 国产亚洲午夜精品一区二区久久| 在现免费观看毛片| 久久av网站| 哪个播放器可以免费观看大片| 婷婷色综合www| 欧美97在线视频| 肉色欧美久久久久久久蜜桃| 极品人妻少妇av视频| 9热在线视频观看99| 久久国产精品男人的天堂亚洲| 国产在线免费精品| 色网站视频免费| 这个男人来自地球电影免费观看 | 肉色欧美久久久久久久蜜桃| 精品一区二区三区四区五区乱码 | 久久青草综合色| 亚洲国产最新在线播放| 美女视频免费永久观看网站| 国产av一区二区精品久久| 中文天堂在线官网| 男女边吃奶边做爰视频| 中文字幕精品免费在线观看视频| 丁香六月天网| 亚洲av国产av综合av卡| 亚洲精品国产色婷婷电影| av女优亚洲男人天堂| 精品国产一区二区三区四区第35| 亚洲精品视频女| 一级毛片电影观看| 欧美av亚洲av综合av国产av | 校园人妻丝袜中文字幕| 久久狼人影院| 各种免费的搞黄视频| 在线观看一区二区三区激情| 亚洲欧美成人综合另类久久久| 九九爱精品视频在线观看| 亚洲美女黄色视频免费看| 啦啦啦啦在线视频资源| 91午夜精品亚洲一区二区三区| 咕卡用的链子| 色94色欧美一区二区| 久久ye,这里只有精品| 亚洲欧美成人精品一区二区| www.精华液| 少妇被粗大的猛进出69影院| av又黄又爽大尺度在线免费看| 欧美av亚洲av综合av国产av | 综合色丁香网| 九草在线视频观看| 久久久久久久精品精品| 欧美av亚洲av综合av国产av | 少妇猛男粗大的猛烈进出视频| 丝袜喷水一区| 日本猛色少妇xxxxx猛交久久| 涩涩av久久男人的天堂| 可以免费在线观看a视频的电影网站 | av在线播放精品| 国产精品国产av在线观看| 青春草亚洲视频在线观看| 国产精品香港三级国产av潘金莲 | 日韩av在线免费看完整版不卡| 水蜜桃什么品种好| 成年美女黄网站色视频大全免费| 曰老女人黄片| 国产成人精品无人区| 日韩 亚洲 欧美在线| 亚洲av日韩在线播放| 老司机影院成人| 免费观看av网站的网址| 一本久久精品| 电影成人av| 欧美国产精品va在线观看不卡| av网站免费在线观看视频| 一本大道久久a久久精品| 最近最新中文字幕大全免费视频 | 又粗又硬又长又爽又黄的视频| 九九爱精品视频在线观看| 99热网站在线观看| 久久久久久久大尺度免费视频| 久久久精品国产亚洲av高清涩受| 狠狠精品人妻久久久久久综合| 菩萨蛮人人尽说江南好唐韦庄| 久热久热在线精品观看| 免费人妻精品一区二区三区视频| 超碰97精品在线观看| 亚洲天堂av无毛| 亚洲三区欧美一区| 久久99一区二区三区| 国产精品熟女久久久久浪| 人妻一区二区av| 久久97久久精品| 久久av网站| 啦啦啦在线免费观看视频4| 人人妻人人爽人人添夜夜欢视频| 男男h啪啪无遮挡| 欧美人与善性xxx| 在线观看免费高清a一片| 日本黄色日本黄色录像| 精品卡一卡二卡四卡免费| 久久久久国产网址| 卡戴珊不雅视频在线播放| 国产成人精品福利久久| 超碰成人久久| 亚洲一区二区三区欧美精品| 精品少妇一区二区三区视频日本电影 | 久久久a久久爽久久v久久| 多毛熟女@视频| 日本欧美视频一区| 久久午夜综合久久蜜桃| 熟女av电影| 国产成人免费无遮挡视频| 捣出白浆h1v1| 国产淫语在线视频| 亚洲av欧美aⅴ国产| 老司机影院成人| 亚洲av电影在线进入| 大片免费播放器 马上看| 精品酒店卫生间| xxxhd国产人妻xxx| 精品一区二区三卡| 一区福利在线观看| 欧美激情高清一区二区三区 | 女的被弄到高潮叫床怎么办| 亚洲经典国产精华液单| 精品一区在线观看国产| 欧美成人精品欧美一级黄| 熟女电影av网| 两个人看的免费小视频| 男女午夜视频在线观看| 免费人妻精品一区二区三区视频| 亚洲综合精品二区| 欧美老熟妇乱子伦牲交| 七月丁香在线播放| 久久精品久久久久久久性| 亚洲人成77777在线视频| 日日爽夜夜爽网站| 午夜福利乱码中文字幕| 久久久久精品久久久久真实原创| 免费黄色在线免费观看| 天天操日日干夜夜撸| 日日摸夜夜添夜夜爱| 另类亚洲欧美激情| 欧美日韩一级在线毛片| 婷婷色av中文字幕| 午夜福利网站1000一区二区三区| 看免费av毛片| 亚洲精品久久成人aⅴ小说| 国产男人的电影天堂91| 中文字幕制服av| 日韩人妻精品一区2区三区| 最近手机中文字幕大全| 日本免费在线观看一区| 看免费av毛片| 国产免费又黄又爽又色| 人人妻人人澡人人看| 亚洲国产成人一精品久久久| 日韩av在线免费看完整版不卡| 一本久久精品| 美女脱内裤让男人舔精品视频| 国产综合精华液| 国产精品香港三级国产av潘金莲 | 免费看av在线观看网站| 成人毛片60女人毛片免费| 男男h啪啪无遮挡| 深夜精品福利| 国产一区二区 视频在线| 国产免费又黄又爽又色| 国产国语露脸激情在线看| 丰满饥渴人妻一区二区三| 丝袜人妻中文字幕| 99久国产av精品国产电影| av有码第一页| 国产高清国产精品国产三级| av视频免费观看在线观看| 日本vs欧美在线观看视频| 另类精品久久| 天天躁夜夜躁狠狠久久av| 一级a爱视频在线免费观看| 亚洲成色77777| 午夜激情久久久久久久| 超碰成人久久| 国产成人一区二区在线| 汤姆久久久久久久影院中文字幕| 国产高清国产精品国产三级| 国产爽快片一区二区三区| 久久久亚洲精品成人影院| 午夜免费鲁丝| 青春草视频在线免费观看| 亚洲av中文av极速乱| 日本爱情动作片www.在线观看| 美女中出高潮动态图| 各种免费的搞黄视频| 丝袜喷水一区| 亚洲欧美成人综合另类久久久| 久久久久国产精品人妻一区二区| 精品国产乱码久久久久久男人| 亚洲国产最新在线播放| 国产av码专区亚洲av| 亚洲欧美中文字幕日韩二区| 69精品国产乱码久久久| 成年动漫av网址| 卡戴珊不雅视频在线播放| 天天躁夜夜躁狠狠久久av| 丝瓜视频免费看黄片| 国产片内射在线| 久久久久久久久久人人人人人人| 麻豆乱淫一区二区| av免费在线看不卡| 在线精品无人区一区二区三| 亚洲av免费高清在线观看| 男人操女人黄网站| 国产淫语在线视频| 毛片一级片免费看久久久久| 亚洲精品日韩在线中文字幕| 日日啪夜夜爽| 久久久久久久久免费视频了| 久久精品国产亚洲av高清一级| 美女视频免费永久观看网站| 一边亲一边摸免费视频| 色94色欧美一区二区| 少妇熟女欧美另类| 亚洲综合色惰| 国产精品二区激情视频| www日本在线高清视频| 欧美bdsm另类| 久久毛片免费看一区二区三区| 欧美日韩国产mv在线观看视频| 一区二区三区精品91| 亚洲国产毛片av蜜桃av| 校园人妻丝袜中文字幕| 亚洲av欧美aⅴ国产| 老熟女久久久| 两个人看的免费小视频| 亚洲av日韩在线播放| 男人舔女人的私密视频| av一本久久久久| 丝袜美腿诱惑在线| 国产精品久久久久久久久免| 国产精品欧美亚洲77777| 大陆偷拍与自拍| 国产无遮挡羞羞视频在线观看| 欧美精品国产亚洲| 久久久国产精品麻豆| 人妻系列 视频| 日韩中文字幕视频在线看片| 大香蕉久久网| 一级黄片播放器| 日韩三级伦理在线观看| 精品酒店卫生间| 韩国av在线不卡| 日韩中文字幕视频在线看片| 精品酒店卫生间| 国产日韩欧美视频二区| 亚洲精品乱久久久久久| 国产在线免费精品| 亚洲精品日本国产第一区| 日韩不卡一区二区三区视频在线| tube8黄色片| www.熟女人妻精品国产| 丝袜美腿诱惑在线| 一级黄片播放器| 亚洲成国产人片在线观看| 国产不卡av网站在线观看| 日本欧美国产在线视频| 大香蕉久久成人网| 久久精品熟女亚洲av麻豆精品| 国产成人av激情在线播放| 老鸭窝网址在线观看| 另类精品久久| 国产精品秋霞免费鲁丝片| 深夜精品福利| 性色avwww在线观看| 最黄视频免费看| 国产成人精品一,二区| 精品少妇内射三级| 国产黄频视频在线观看|