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

    Energy Control of Plug-In Hybrid Electric Vehicles Using Model Predictive Control With Route Preview

    2021-11-07 02:23:12YangZhaoYanguangCaiandQiwenSong
    IEEE/CAA Journal of Automatica Sinica 2021年12期

    Yang Zhao,Yanguang Cai,and Qiwen Song

    Abstract—The paper proposes an adoption of slope,elevation,speed and route distance preview to achieve optimal energy management of plug-in hybrid electric vehicles (PHEVs).The approach is to identify route features from historical and realtime traffic data,in which information fusion model and traffic prediction model are used to improve the information accuracy.Then,dynamic programming combined with equivalent consumption minimization strategy is used to compute an optimal solution for real-time energy management.The solution is the reference for PHEV energy management control along the route.To improve the system’s ability of handling changing situation,the study further explores predictive control model in the realtime control of the energy.A simulation is performed to model PHEV under above energy control strategy with route preview.The results show that the average fuel consumption of PHEV along the previewed route with model predictive control (MPC)strategy can be reduced compared with optimal strategy and base control strategy.

    I.INTRODUCTION

    T HE increasing fuel price and pollution problem have aroused the concern for electric vehicle in transportation area.One of the most promising technologies for fuel economy and emission reduction is the plug-in hybrid electric vehicle,in which the vehicle power-train is built with an engine enhanced by a battery and electric machine.The challenge faced by vehicle power-train is the control strategy design of energy management system(EMS),which determines the power split between the engine and battery.In the past decades,there are various studies focusing on the EMS to maximize the energy efficiency for plug-in hybrid electric vehicles (PHEVs).The researches generally believe that a near-optimal fuel economy can be achieved only if the vehicle knows the future traffic condition [1],[2].A common way to get future information is to adopt a GPS-based navigation system with the driver inputting the destination,which provides speed and topography data to the EMS along the route.However,there are limitations for the method.The predicted traffic condition is not always the actual driving condition.In addition,the existing energy control strategies do not consider the information error,the computational demand should allow for the reasonable limitation,adding constraint on the vehicle model.

    However,the recent development trend of automotives is to connect the vehicle with the traffic information center through 3G network.The traffic information center will provide realtime and predicted traffic condition information for the vehicle.Then,it is possible for the vehicle to achieve dynamic energy management.When the user inputs the destination and applies for the dynamic navigation,these data will be transmitted into traffic information center server,which identifies further traffic condition in the route and send back to the vehicle.An optimal strategy will be precomputed based on future information for the vehicle,as showed in Fig.1.

    There are two main constructing blocks of such a system,including route prediction and optimization of the EMS.Several studies have researched different short-term prediction techniques that are used to predict route characteristic from purely historical or real-time traffic data and relate the route to the vehicle operation[3],[4],however this prediction based approach has a shortcoming that it cannot reflect the actual traffic condition accurately in a required prediction horizon.And the information error will affect the performance of energy optimization.A PHEV’s excellent energy economy is not only determined by the absorbed information,but also by the energy control strategy.A number of optimization strategies have been proposed for the EMS combining with traffic condition information[5].Given the preview of vehicle speed,slope,and drive power demand in a drive cycle,dynamic programming can be used to identify the optimal solution on the powertrain operation globally or in a certain time horizon.The previous study in[6]has showed that know ledge of future trip information can contribute to fuel economy improvement in PHEVs.However,the requirement for future trip information makes it difficult to operate at real time for dynamic control strategy.Existing work lacks integrating route preview information with real-time implementable algorithm[7].The dynamic programming can calculate the estimated power split based on only the future power demand regardless of past or current state of vehicle.The presumption for dynamic programming computation that may be useful for a route is that the immediate future power demand does not change greatly.

    Fig.1.The workflow of energy management with route preview.

    Based on the above survey,this paper makes contributions to the literature in two aspects:1)it provides a route preview with information fusion of real-time and historic data to improve the accuracy of prediction;2)proposes a model predictive control model to predict the energy split dynamics with optimal control while retaining the linearity of state evolution equation.

    The methodology used in the paper is to identify routes from GPS driving data,using historical and real-time speed data.The data fusion model and short-term prediction model are used for the route speed prediction.As the slope and elevation are static,only speed is necessary to be predicted.A fter the route character is predicted,the route is identified,which is used to model the driving condition along the route.Then,dynamic programming combining with an equivalent consumption minimization strategy (ECMS) as the real-time control strategy is proposed to compute an optimal solution based on a simplified vehicle model.A fter the optimal trace is determined,a model predictive control model is introduced to optimize the energy system at real time within a certain prediction horizon while tracing the optimal trajectory.The objective of model predictive control (MPC) control is to improve the system’s ability of handling changing situation by measuring and analyzing current state while adhering to constraints of vehicle character and future traffic condition.Recent scholars have put emphasis on the MPC model for the hybrid electric vehicles energy control.The MPC model is used for working with different operating conditions.Through calculating the future control methodology with quadratic programming,the manipulated variable is at real time as well as the changing MPC.The proposed MPC can improve the frequency stabilization of the system [2].Some scholars also adopt the nonlinear model predictive control for optimizing the energy consumption.The fast and efficient control orientation model integrated with nonlinear MPC model is developed to achieve the global optimum energy control[3].However,there is a lack of MPC control integrating with the previewed route information,which also affects the energy control for PHEV.

    Outline: The paper is organized as follows.Following the introduction,the route identification procedure and the corresponding route parameter prediction are explained.The next section describes the simplified vehicle model used during the precomputations.Then,the section introduces how dynamic programming and ECMS can be used to precompute an optimal control strategy for a route.The detailed MPC model used for real-time energy control is explained in the following section.In the final part of the paper,the simulation study and its results are presented before the paper ends with a discussion and conclusion.

    II.ROUTE PREV IEW

    In this paper,the concept of a route preview is defined as preview of slope,elevation,vehicle speed and trip distance for the future route.The future power demand is determined by the upcoming road slope and future speed.The road terrain information can be obtained from the onboard 3-D maps and the GPS navigation system [8].The future speed profile of the route is estimated through using real-time data and historical data.

    Trip features: To facilitate the estimation,each trip is associated with a finite number of features,e.g.,driving length,starting time,slope,elevation,and GPS coordinates at different positions.

    A.Route Prediction

    The process of route preview can be viewed as the task of estimates the route features in n links in the specific route.For each link,they have the same feature with the route,and all links get together to make the route [9].

    In the study,we adopt the traffic data collected by Guangzhou Transportation Department.The route prediction based on historic and real-time data is achieved with the following procedure.

    Speed prediction.The system archives historic and real-time data,the prediction of the speed based on these data is separate.The ground traffic in the same site over several logged data (a week or month) has similar intra-day trend,which has been used as approximate predicted trend.Researchers have proposed several different methods to identify the trend.We adopt a relatively practical method,that is moving average method to predict the historic speed trend [10].In the realtime prediction,we emphasize on improving the timeliness of the prediction,a real-time recurrent learning(RTRL)algorithm is proposed to predict future speed in a certain horizon with different steps [11].

    Information fusion.In order to take into account the timeliness of the speeds as well the periodic trend,the realtime prediction model is considered to integrate the historic trend.In the study,we relate historic trend with current traffic state through building the information fusion model,and then predicting future steps of the link speed.Mathematically,the information fusion methodology is represented by the following equations:

    wheref(t) is fused speed,λ(t) is the weight of history data in the prediction at timet,h(t) is the historic trend,r(t) is the predicted speed with RTRL.

    Equation (1) is a convex combination which is a standard information fusion technique used for speed estimation.

    B.Route Modeling

    In order to optimize the energy management of PHEV in a specified route,the traffic condition along the route should be modeled in a reasonable way.The route modeling will be processed with following steps:

    Step 1:For a defined route,the origin,destination and departure time are determined by the user.

    Step 2:Route operation features assignment.It is assumed that there areklinks in the route.A ll links are assigned with q features,and a link feature data matrix is formed.The row vector corresponds to the feature vectors of links,including the terrain,speed and length,etc.And the column vectors reflect the observations of the features.As the speed feature varies with time,the speed observation will be a matrix,which represents different speed profiles in k link at different time.

    Step 3:The time consumption for each link at any time point can be obtained based on the speed prediction and link length.Estimate the predicted leaving timel1at first link and the speed of the second links2with timel1,which is used to calculate the travel time in the second link.Repeat the process until arriving the destination and get the timelkand the speedskin the final link.

    Step 4:Upgrade the feature vectors.

    III.VEHICLE MODELING

    A post-transmission parallel PHEV is selected for this study.In the vehicle,the electric motor connects with the front axle,the engine is coupled by a clutch and an automatic transmission [12].In the vehicle model,an inverse method is used to determine the force acting on powertrain,which means that the required wheel torqueTrfollows with the input speed and road slope,it is calculated by the following equation

    wherevis the vehicle speed,a is the acceleration,θis the road slope,ρis the density of air,Cd Ais the air drag resistance,fris the rolling resistance,mis the vehicle mass,andmeis the equivalent vehicle mass.The traction torque of the powertrainTpis determined by

    whereTmis the motor torque,Teis the engine torque,andTbis the torque of the friction brakes.rfdenotes the ratio of the final gear,and the corresponding efficiencyηfdepends on the sign of the torque demand at the wheels.The gears,i= 1,...,5 are represented by a drive ratiorgb,iand a mechanical efficiencyηgb,i.It is assumed that torque responses are instantaneous and that gear shifts and the enginestate transitions are both instantaneous and lossless [8].The instantaneous fuel mass rate is mapped from the engine torque,and the instantaneous fuel cost g is thus given by

    wherecfdenotes the price of fuel,andeonis the engine state.In addition,the output electrical power of the battery is related with motor torque.The resulting battery powerPbis given by

    wherePais the auxiliary electrical load.The coefficientsc0:1andd0:2are determined by linear least squares from the engine and motor maps.An equivalent circuit is considered to model a Li-ion battery consisting ofncbattery cells connected in series,the parameters are listed in Table I.The battery with an open-circuit voltage and the series connection with a constant resistanceRidepends on the battery state of charge (SOC)[13].The SOC dynamics are described by

    whereQis the battery capacity.As for the control signaluin the vehicle model,it contains the engine-state decision information,the torques of the engine and the electric motor,i.e.,

    whereUis the torque constraints of the engine and motor.In addition,the control signal should be selected,for example,Tpshould be equal toTr.However,the choice of gear is excluded in the control signal,which is determined by a separate gearbox controller.

    TABLE IVEHICLE DATA

    IV.ENERGY MANAGEMENT STRATEGY

    A PHEV can operate in charge-depleting (CD) mode and charge-sustaining(CS)mode.When the battery SOC is nearly depleted,the PHEV will change into CS mode through blending operation of the engine and the electric motor [8].

    At the beginning of the trip,the battery SOC is high,the PHEV operates in CD mode: the battery SOC is nearly depleted with either all-electric mode or mode of blending the electric motor and combustion engine.The electric grid energy is generally economical compared with gasoline fuel energy.Thus,it is better to operate the PHEV with all-electric mode in a short trip.For the trips out of the all-electric range,hybrid operation of the electric motor and the combustion engine is presented to achieve fuel saving compared with all-electric mode.With preview of future driving condition,the vehicle can optimize the operation suggestion in all-electric CD mode or blended CD mode.

    The energy management control problem is to maximize fuel saving for a PHEV along a route,while adhering to the constraint and dynamics of the vehicle model.The optimal control problem can be formulated as a minimization of the overall expected energy cost

    whereJx t(SOC(t))is the terminal cost at time t from positionxtto the end position at timetf; current SOC isSOC(t);uis the power-split ratio,which is also the control vector.The termφ(·) represents the recharge cost of the finalSOCffromSOC(t).As the paper focuses on the charge depletion optimization within a charge cycle,φ(·) is considered to remove in this study.In addition,the solution of the optimization problem should adhere to constraints including powertrain model,physical limitation and SOC operation.

    Due to the constraints and nonlinearities,there is no direct analytical solution for the optimization problem.To solve the optimal energy control problem,a two-scale DP solution is provided.In a higher layer,the dynamic programming is used to plan the battery usage with preview traffic information.In the lower layer,an EMCS real-time control strategy tracks the trajectory in the higher layer while subject to the real-time control of energy.

    With preview of future velocity and future power demand,the control problem illustrated in (8) that subjects to model equations and constraints can be solved according to Bellman’s optimality principle and dynamic programming.In the vehicle model,the power-split controller adopts an ECMS strategy,the optimization problem is transferred to minimize instantaneous equivalent fuel ratedefined as

    wherePe(t,u) is the net power charged to the battery;s(t) is the ECMS equivalence factor;andHfis the fuel lower heating value.

    The Hamiltonian for the cost function (8) is as follows:

    whereλ(t) is the co-state,from which the optimal value is determined by future power demands.Applying Pontryagin’s minimum principle,the co-stateλ(t) has dynamics

    As the optimal co-state depends on the future power demands,it is obtained based on future route conditions and current battery’s SOC.However,the future power demand can be obtained because of the uncertainness.We provide an estimation of optimal co-state ?λ(t) through using preview information.

    From the Hamilton-Jacobi-Bellman equation,the optimal co-stateλ?(t) can be written as

    whereλ?(SOC(t),xt) is the optimal co-state determined bySOC(t) and positionxt,and(SOC(t)) is the optimized cost-to-go in (8).

    Through estimating the upcoming power demand based on the traffic preview,a dynamic programming can be used to obtain the estimation of the optimal cost-to-go(SOC(t)).Then,the estimation of optimal co-state is written as

    V.MODEL PREDICTIVE CONTROL FOR ENERGY CONTROL

    A.Model Predictive Control

    Model predictive control (MPC) is a methodology to repeatedly determine optimal values online in a rolling horizon for control measures,which uses optimization and prediction model for the optimization problem control [14].The MPC can handle different constraints on system input and state.MPC also builds a feedback mechanism as it adopts a receding horizon approach.

    Fig.2.Representation of MPC.

    Fig.3.Conceptual MPC control structure.

    The work mechanism of MPC is as follows (see Fig.2):the control time intervalTctrlis the time interval,defined ast=kTctrl.At each control stepk,the current statex(k) of the system is measured and the predictive state overkstep horizon is determined.Next,an optimization model combining with system model is used to derive the control inputsu(k),...,u(k+Np ?1) that optimize a given performance criterion over a given horizon [kTctrl,(k+Np)Tctrl]subject to constraints,whereNpis the prediction horizon.To improve the computation speed,we introduce the control horizonNc(< Np),and there is a form ofu(k+j) =u(k+j ?1)forj=Nc,...,Np ?1 [14].

    The optimal control values are imputed in the system with a moving horizon concept: at each control time step,only the first control vectoru?(k) of the optimal control sequenceu?(k),...,u?(k+N c ?1) is applied to the system.Then,the prediction horizon shifts one step forward.In the next state,the optimal control problem is solved again based on the new state and the new optimal control inputs.The process is repeated until the final prediction horizon is got.As the moving horizon concept includes a feedback mechanism,it allows the system to reduce the effect of disturbances and improve the accuracy.

    B.MPC Strategy for the Route

    MPC can repeatedly solve optimization problems online in a rolling horizon to derive a sequence of control decisions[15].We focus on how MPC minimizes the cost and improves vehicle energy control.

    The process of MPC,seen in Fig.3,can be described in the following steps:

    Step 1:Prediction model.The prediction model of energy demand can be generally described,for

    whereSOC(k+1) is the predicted state,SOC(k) is current state,u(k) is the control input,fis the shifted function.

    Step 2:Control Inputs.At each control step,the MPC controller obtains or predicts the current state of the vehicle energy state.Since the future slope,speed and elevation are the basic components affecting the energy consumption,control signalufor the MPC problem of control stepkincludes slope,speed and elevation in the route over the control horizon.To take the future demand of energy into account in the optimization process,we can use autonomous control for vehicle energy control betweent= (k+N c)Tctrlandt= (k+Np)Tctrl.With MPC energy control,the vehicle would be assigned with the energy model of its immediate predecessor energy splitting [14].

    Step 3:Optimization problem.Assuming that the control intervalThasNsteps,the optimization problem of MPC can be given as

    whereg(kctrl) is the future control input vector for control steps.The performance criteriaJperf(k) for MPC is the total fuel comsumption in the route,which is evaluated during[kTctrl,(k+Np)Tctrl].Moreover,to prevent great shifting from the DP based optimal value,the performance criteria adds a penalty on variations between MPC prediction value and DP based optimal value,which results in the total performance function

    wherefoptimal(k+j) is the look-up value at timek+jfor optimal control strategy.

    Step 4:Rolling horizon.After the control input is obtained from the optimization,the first optimal vector is transferred to the process and implemented.When the prediction model comes to the next control step,it is updated with realtime traffic data,vehicle state,and the prediction horizon is shifted one step forward,the optimization process begins again.The rolling horizon process makes the control loop close,and enables the control system obtain feedback from the real driving system.Thus,the MPC controller is robust to disturbances [15].

    C.Optimization Methods

    One of the challenges in MPC approach is the extensive demand on computation.The prediction model is nonlinear,and the optimization of MPC is a nonconvex optimization problem solved online.Thus,it is essential to select a suitable optimization technology for deriving feasible optimal control vectors.In the study,global or multistart local optimization methods are necessary,including multistart sequential quadratic programming [16],multistart pattern search [17],genetic algorithms [18],or simulated annealing.

    D.Performance Measures

    Prediction error indices such as mean absolute relative error(MARE) and root relative square error (RRSE) are computed with equations,respectively.

    To compare and analyze the prediction error indices in different time intervals,the real velocity and the predicted velocity will be analyzed with the prediction error analysis,which indicate the performance of MPC strategy.

    VI.SIMULATION AND VEHICLE TEST

    To test the proposed MPC strategy in more realistic driving scenarios,a simulation study is implemented using real-world driving condition provided by Guangzhou traffic analysis system,whose database contains geographic information system(GIS) data and the floating car data from about 16 000 taxies,gathered since 2013 in Guangzhou area.The study simulates a specified route in the system.

    The data for driving condition are classified into the training period and the validation period.The training period is four weeks while the validation duration is 1 week.In the training period,the route data is used to represent the route feature as described in Section II.In the route modeling,the trip is assigned with link length,speed profile,terrain,slope,i.e.,the validation period is reserved for simulation.

    Only know ledge on trip distance is assumed and the simulation results provide key insights to the fuel energy (FE) improvement potentials of MPC toward real-world applications.The traffic condition information is described in Fig.4.

    The simulation is studied with federal test cycle and realworld traffic profile.A default rule-based strategy without incorporating route information is used for comparison purpose.The default strategy applies a fixedλcalibration that results a 35 mile PER in HWYFET drive cycle at a maximal rate of charge depletion.In the MPC control,the prediction horizon is defined as 120 s.About 2%?7%FE improvement has been observed in the simulations for MPC strategy.For example,the simulation result from a simulation case with a 120 km trip distance and 75% initial battery SOC is plotted in Fig.5.By incorporating the trip information,the fuel consumption for MPC is 5.4704 kg while fuel economy for basic control strategy is 5.931 kg.An 8.65% FE improvement is achieved in comparison with the basic control process.

    Fig.4.Speed,elevation,and grade in the route.

    Fig.5.Battery SOC and control setpoint for MPC and basic control strategy.

    The future study will explore the improvement of MPC strategies compared with optimal control strategy.In Fig.6,the optimal battery SOC EMCS process is plotted together with that from the MPC control process for comparison.The MPC controlledλsetpoint converges toλ?adaptively.Given the know ledge of driving conditions in a trip,simulation validation is carried out using the optimized SOC profile as reference and the recommendedλ?as feed forwardλsetpoint.The fuel economy for optimal strategy is 5.731 kg.Fig.6 shows that the fuel consumption for MPC is 5.4704 kg,there is a 5.35%FE improvement achieved for MPC strategy comparing with optimal strategy.

    Fig.6.Battery SOC and control setpoint for MPC and optimal control strategy.

    One of the advantages of MPC is its robustness to the uncertainties.To test the robustness of the operation strategy in terms of inaccuracy of the estimated velocity profile,the value of the traffic limitation as well as the configuration parameters of the driver model were disturbed with variables subject to a normal distribution and a standard deviation of 20% of the original value.Then,simulations over the same driving trip are repeatedly performed with these uncertainties and the results are statistically analyzed.The deviations between predicted and actual drive trip result in a far more frequent re-calculation of the variable.However,the fuel consumptions of the disturbed system deviate slightly from the consumptions of the undisturbed system.Fig.7 shows an example of the simulation.The fuel consumption in disturbed system is 5.5801 kg,which is 1.9% larger than the normal system and 2.7% less than the optimal control system.The performance of the disturbed system (20% variance) with MPC strategy is worse than the normal system with MPC strategy while it is still better than the optimal control system.The result shows that the MPC control system has a good robustness in the disturbed environment.In addition,the result indicates that MPC control system has a good robustness in the energy consumption control.

    With the increasing prediction horizon of the route,the vehicle can obtain more future operation state information,the fuel saving effect will be better.However,when the forecast horizon increases to a certain distance,predictive control computation efficiency is very low.Therefore,it is necessary to select a proper horizon,making the computed horizon forecast with a better fuel economy.

    Fig.7.Battery SOC for MPC,optimal control and MPC (20% variance)strategy.

    We can see from Table II,with the increase of the prediction horizon,the equivalent fuel consumption is reduced gradually.With prediction horizon increasing,the vehicle can obtain more future traffic information,and modify to special work conditions of the engine earlier.When the forecast horizon is more than 160 s,equivalent fuel consumption slightly increased,this is led by the accumulation calculation error of MPC.Considering a reasonable fuel economy,we select a length of 140 s as the forecast horizon.Comparing with the previous control strategies,it is obvious that the fuel economy is better no matter how long the prediction horizon is.Thus,in the actual energy control,it is necessary to combine the MPC model and the route preview information,thus optimizing the system and obtaining the maximum operation effectiveness.

    TABLE IITHE FUEL ECONOMY WITH DIFFERENT PREDICTION HORIZON

    V II.CONCLUSIONS

    This paper presents real-time control strategies for PHEV energy management that takes advantage of traffic information preview.The proposed strategies aim at fuel minimization for PHEV.Initially,the route features are identified with information fusion model and prediction model for historical and realtime traffic data.The objective of fuel saving is achieved at three levels:1)the route preview information is used for global optimization estimation of fuel-electricity usage solution; 2)the ECMS is applied to achieve real-time optimal energy control with information of instantaneous power demand,current value of batter’s SOC and the global optimization parameters; 3) the MPC is designed to control the electricity consumption ratio along the route based on future power demand,historical and current value of battery’s SOC and instantaneous optimization parameters.The simulation study shows that there is fuel saving by applying the energy control strategies with route preview.MPC control indicates 8.65%improvement compared with base control strategy and 5.35%improvement compared with optimal control combined with ECMS.Another simulation indicates that the MPC strategy has a good robustness.In addition,the prediction horizon of MPC affects the energy control effect,a 140 s prediction horizon leads to a better control result.

    亚洲五月色婷婷综合| 天堂俺去俺来也www色官网| 99久久99久久久精品蜜桃| 一区二区三区乱码不卡18| 窝窝影院91人妻| 青青草视频在线视频观看| 狂野欧美激情性xxxx| 性少妇av在线| 两性夫妻黄色片| 亚洲欧美成人综合另类久久久| 别揉我奶头~嗯~啊~动态视频 | 真人做人爱边吃奶动态| 男男h啪啪无遮挡| 国产成人免费无遮挡视频| 亚洲欧美日韩高清在线视频 | xxxhd国产人妻xxx| 亚洲综合色网址| 最新在线观看一区二区三区| 亚洲伊人久久精品综合| 久久精品久久久久久噜噜老黄| 国产成人一区二区三区免费视频网站| av超薄肉色丝袜交足视频| 亚洲精品美女久久av网站| 亚洲av美国av| 巨乳人妻的诱惑在线观看| 国产野战对白在线观看| 男女午夜视频在线观看| 亚洲成人免费av在线播放| 老熟妇仑乱视频hdxx| 婷婷成人精品国产| 精品亚洲乱码少妇综合久久| avwww免费| 亚洲专区中文字幕在线| 丝袜脚勾引网站| 免费日韩欧美在线观看| 亚洲av片天天在线观看| 老司机亚洲免费影院| 亚洲欧美精品综合一区二区三区| 九色亚洲精品在线播放| www日本在线高清视频| 亚洲精品国产av成人精品| 老司机深夜福利视频在线观看 | 精品人妻熟女毛片av久久网站| av视频免费观看在线观看| 黄片小视频在线播放| 久久女婷五月综合色啪小说| 国产欧美日韩精品亚洲av| 亚洲男人天堂网一区| 80岁老熟妇乱子伦牲交| 免费在线观看完整版高清| 亚洲精品一二三| 美女午夜性视频免费| 免费少妇av软件| 男女午夜视频在线观看| 男人舔女人的私密视频| 丝袜喷水一区| 丝袜喷水一区| av视频免费观看在线观看| 久久狼人影院| 伊人久久大香线蕉亚洲五| 婷婷色av中文字幕| 亚洲熟女毛片儿| 热99国产精品久久久久久7| 97人妻天天添夜夜摸| 高清欧美精品videossex| 亚洲欧美一区二区三区久久| 国产免费一区二区三区四区乱码| 亚洲专区字幕在线| 国产老妇伦熟女老妇高清| 国产av国产精品国产| 国产日韩一区二区三区精品不卡| 亚洲av电影在线观看一区二区三区| 色播在线永久视频| 欧美国产精品va在线观看不卡| av免费在线观看网站| 亚洲少妇的诱惑av| 男女午夜视频在线观看| 搡老岳熟女国产| 亚洲第一欧美日韩一区二区三区 | 性高湖久久久久久久久免费观看| 国产av精品麻豆| 国产男女超爽视频在线观看| 成年美女黄网站色视频大全免费| 亚洲国产欧美一区二区综合| 国产成人影院久久av| 美女主播在线视频| 在线观看免费日韩欧美大片| 十八禁人妻一区二区| 最近最新免费中文字幕在线| 亚洲欧美一区二区三区黑人| 免费看十八禁软件| 亚洲精品中文字幕一二三四区 | 国产精品二区激情视频| 亚洲第一青青草原| 欧美日韩亚洲综合一区二区三区_| 高清黄色对白视频在线免费看| 成人手机av| 亚洲国产成人一精品久久久| 黑人巨大精品欧美一区二区mp4| 美女国产高潮福利片在线看| 中文字幕高清在线视频| 91精品伊人久久大香线蕉| av不卡在线播放| 欧美日韩黄片免| 好男人电影高清在线观看| 亚洲中文字幕日韩| 韩国精品一区二区三区| 国产精品久久久人人做人人爽| 亚洲狠狠婷婷综合久久图片| 两个人视频免费观看高清| av视频在线观看入口| 午夜精品久久久久久毛片777| 中出人妻视频一区二区| 亚洲欧美精品综合久久99| 校园春色视频在线观看| 精品久久久久久成人av| 亚洲男人的天堂狠狠| 国产精品野战在线观看| 日韩精品中文字幕看吧| 亚洲 国产 在线| 男女下面进入的视频免费午夜| 美女高潮喷水抽搐中文字幕| 18美女黄网站色大片免费观看| 午夜福利成人在线免费观看| 波多野结衣巨乳人妻| 激情在线观看视频在线高清| 亚洲一区二区三区色噜噜| 91字幕亚洲| 午夜免费成人在线视频| 精华霜和精华液先用哪个| 亚洲av美国av| a级毛片在线看网站| 久久久国产成人免费| 日本在线视频免费播放| 色综合婷婷激情| 久久精品91蜜桃| 久久中文看片网| 给我免费播放毛片高清在线观看| 久久久国产成人精品二区| 最新美女视频免费是黄的| 国产午夜精品论理片| 国产爱豆传媒在线观看 | 好男人在线观看高清免费视频| 国产av又大| 久久草成人影院| 国产探花在线观看一区二区| 女人爽到高潮嗷嗷叫在线视频| 少妇粗大呻吟视频| av天堂在线播放| 亚洲美女视频黄频| 十八禁网站免费在线| 欧美日韩乱码在线| 美女 人体艺术 gogo| 淫妇啪啪啪对白视频| 亚洲自偷自拍图片 自拍| 岛国在线观看网站| 国产午夜福利久久久久久| 91字幕亚洲| 亚洲18禁久久av| 女生性感内裤真人,穿戴方法视频| 亚洲一码二码三码区别大吗| 丝袜人妻中文字幕| 久久精品国产99精品国产亚洲性色| 亚洲在线自拍视频| 日韩精品中文字幕看吧| 欧美又色又爽又黄视频| 麻豆国产97在线/欧美 | 熟女少妇亚洲综合色aaa.| 精品一区二区三区av网在线观看| 在线十欧美十亚洲十日本专区| 麻豆成人av在线观看| 一区福利在线观看| 九九热线精品视视频播放| 日韩高清综合在线| 国产成人欧美在线观看| 日韩欧美在线乱码| 97碰自拍视频| 中文字幕人妻丝袜一区二区| 给我免费播放毛片高清在线观看| 亚洲国产精品久久男人天堂| 久久久精品大字幕| 欧美精品亚洲一区二区| 亚洲美女黄片视频| 中文亚洲av片在线观看爽| 又粗又爽又猛毛片免费看| 欧美av亚洲av综合av国产av| 床上黄色一级片| 久久久久国内视频| 成熟少妇高潮喷水视频| 香蕉久久夜色| 欧美+亚洲+日韩+国产| 国产亚洲欧美在线一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 精品一区二区三区四区五区乱码| 蜜桃久久精品国产亚洲av| 一级毛片高清免费大全| 免费在线观看视频国产中文字幕亚洲| 精品国产超薄肉色丝袜足j| 欧美+亚洲+日韩+国产| 老汉色av国产亚洲站长工具| www.www免费av| bbb黄色大片| 夜夜躁狠狠躁天天躁| 中文亚洲av片在线观看爽| 中文字幕人成人乱码亚洲影| videosex国产| 亚洲五月婷婷丁香| 欧洲精品卡2卡3卡4卡5卡区| 无限看片的www在线观看| 日韩精品青青久久久久久| 国产亚洲欧美在线一区二区| 国产真实乱freesex| 亚洲 欧美一区二区三区| 老熟妇乱子伦视频在线观看| 精华霜和精华液先用哪个| 九九热线精品视视频播放| 国产熟女午夜一区二区三区| 午夜成年电影在线免费观看| 岛国在线观看网站| 日韩成人在线观看一区二区三区| 制服丝袜大香蕉在线| 国产69精品久久久久777片 | 一本精品99久久精品77| 国产av不卡久久| 中亚洲国语对白在线视频| 看黄色毛片网站| 色综合亚洲欧美另类图片| 午夜老司机福利片| 国产午夜福利久久久久久| 亚洲精品中文字幕一二三四区| 中文字幕久久专区| 18禁国产床啪视频网站| 99久久综合精品五月天人人| 97碰自拍视频| 欧美人与性动交α欧美精品济南到| 国产精品一区二区三区四区免费观看 | 国产精品免费视频内射| 亚洲无线在线观看| 黄色a级毛片大全视频| 国产一区在线观看成人免费| 两人在一起打扑克的视频| 欧美三级亚洲精品| 亚洲黑人精品在线| 99国产精品一区二区三区| АⅤ资源中文在线天堂| 成人18禁高潮啪啪吃奶动态图| 人人妻人人澡欧美一区二区| 99在线视频只有这里精品首页| 手机成人av网站| a级毛片a级免费在线| 成人18禁高潮啪啪吃奶动态图| 成人三级黄色视频| 精品一区二区三区四区五区乱码| 久久久久久人人人人人| 色精品久久人妻99蜜桃| 又黄又爽又免费观看的视频| 亚洲九九香蕉| 一区二区三区国产精品乱码| 亚洲熟妇中文字幕五十中出| 国产激情久久老熟女| 午夜免费成人在线视频| 黄色片一级片一级黄色片| 国产精品亚洲一级av第二区| 亚洲国产精品成人综合色| 精品高清国产在线一区| 床上黄色一级片| 精品久久久久久久久久久久久| 日本五十路高清| 亚洲欧美日韩高清专用| 一边摸一边做爽爽视频免费| 夜夜爽天天搞| www日本在线高清视频| 亚洲av日韩精品久久久久久密| 国内久久婷婷六月综合欲色啪| 婷婷精品国产亚洲av在线| 欧美成狂野欧美在线观看| 免费无遮挡裸体视频| 成在线人永久免费视频| 中文字幕av在线有码专区| 黄片小视频在线播放| 熟女少妇亚洲综合色aaa.| 婷婷精品国产亚洲av在线| 1024手机看黄色片| 国产午夜福利久久久久久| 夜夜爽天天搞| 亚洲av五月六月丁香网| 国产成人精品无人区| 91九色精品人成在线观看| 免费高清视频大片| 高潮久久久久久久久久久不卡| 亚洲片人在线观看| 老司机午夜十八禁免费视频| 久久欧美精品欧美久久欧美| 美女 人体艺术 gogo| 亚洲一区二区三区色噜噜| 1024视频免费在线观看| 久热爱精品视频在线9| 日本免费a在线| 看片在线看免费视频| 国产精品,欧美在线| 人妻久久中文字幕网| 国产亚洲精品第一综合不卡| 精品一区二区三区四区五区乱码| 国产三级中文精品| 午夜老司机福利片| 欧美+亚洲+日韩+国产| 亚洲欧美日韩高清专用| 老汉色∧v一级毛片| www.www免费av| 激情在线观看视频在线高清| 99国产精品一区二区蜜桃av| 亚洲一区中文字幕在线| av视频在线观看入口| 国产亚洲精品一区二区www| 国产一区二区激情短视频| 精品一区二区三区av网在线观看| 亚洲av美国av| 黄色a级毛片大全视频| 久久天堂一区二区三区四区| 日韩欧美国产在线观看| 两个人免费观看高清视频| 中出人妻视频一区二区| 国产真人三级小视频在线观看| 国产97色在线日韩免费| 日韩成人在线观看一区二区三区| 国产97色在线日韩免费| 久久国产乱子伦精品免费另类| 午夜激情福利司机影院| 国产三级中文精品| 老司机午夜十八禁免费视频| 最好的美女福利视频网| 久久久国产成人精品二区| 女人高潮潮喷娇喘18禁视频| 国产视频内射| 精品电影一区二区在线| 99久久无色码亚洲精品果冻| 久久精品91无色码中文字幕| 久久久久免费精品人妻一区二区| 叶爱在线成人免费视频播放| 亚洲自偷自拍图片 自拍| av在线天堂中文字幕| 亚洲av片天天在线观看| 极品教师在线免费播放| 熟妇人妻久久中文字幕3abv| 亚洲欧美一区二区三区黑人| 亚洲午夜精品一区,二区,三区| а√天堂www在线а√下载| 韩国av一区二区三区四区| a级毛片a级免费在线| 波多野结衣巨乳人妻| 19禁男女啪啪无遮挡网站| 中文字幕久久专区| 亚洲国产欧美人成| 9191精品国产免费久久| 免费看美女性在线毛片视频| 国产精品乱码一区二三区的特点| 色播亚洲综合网| 欧美精品啪啪一区二区三区| 一个人观看的视频www高清免费观看 | 在线观看舔阴道视频| 90打野战视频偷拍视频| 搡老熟女国产l中国老女人| 黄色视频不卡| 成人一区二区视频在线观看| 久久久久久久午夜电影| 在线观看舔阴道视频| 国产精品精品国产色婷婷| 中文资源天堂在线| 一夜夜www| 757午夜福利合集在线观看| 午夜亚洲福利在线播放| 超碰成人久久| 99久久无色码亚洲精品果冻| 国产视频一区二区在线看| 国模一区二区三区四区视频 | 黄片小视频在线播放| 久久国产精品人妻蜜桃| 国产av不卡久久| 欧美一级a爱片免费观看看 | 成人国产综合亚洲| 午夜福利在线观看吧| 小说图片视频综合网站| 亚洲成av人片免费观看| 日韩免费av在线播放| 国产麻豆成人av免费视频| 久久久久久九九精品二区国产 | 啦啦啦免费观看视频1| 午夜福利成人在线免费观看| 亚洲美女视频黄频| 国产在线精品亚洲第一网站| 亚洲熟妇熟女久久| 国产亚洲精品综合一区在线观看 | 日本成人三级电影网站| 黄片大片在线免费观看| 久久婷婷成人综合色麻豆| 别揉我奶头~嗯~啊~动态视频| 国内少妇人妻偷人精品xxx网站 | 一本大道久久a久久精品| 一夜夜www| 亚洲午夜精品一区,二区,三区| 男女视频在线观看网站免费 | 午夜激情福利司机影院| 少妇人妻一区二区三区视频| 亚洲免费av在线视频| 老司机靠b影院| 日本成人三级电影网站| АⅤ资源中文在线天堂| 欧美极品一区二区三区四区| 18禁观看日本| 在线观看美女被高潮喷水网站 | 女人被狂操c到高潮| 亚洲精品一卡2卡三卡4卡5卡| 一个人免费在线观看电影 | 亚洲精品一卡2卡三卡4卡5卡| 99国产精品一区二区蜜桃av| 成人欧美大片| 亚洲,欧美精品.| 久久精品国产99精品国产亚洲性色| 国产伦在线观看视频一区| 免费电影在线观看免费观看| 国产成人av教育| 午夜两性在线视频| 日韩欧美免费精品| 精品高清国产在线一区| 国产精品一及| 欧美一级毛片孕妇| av天堂在线播放| 亚洲人与动物交配视频| 怎么达到女性高潮| 黄色 视频免费看| 身体一侧抽搐| 国产精品爽爽va在线观看网站| 亚洲avbb在线观看| 亚洲欧美日韩无卡精品| 国产成人aa在线观看| 国产一区二区三区视频了| 日本一二三区视频观看| 久久国产精品影院| 免费搜索国产男女视频| netflix在线观看网站| 九九热线精品视视频播放| 怎么达到女性高潮| 91麻豆精品激情在线观看国产| 国产成人欧美在线观看| 欧美+亚洲+日韩+国产| 久久性视频一级片| 999久久久精品免费观看国产| 久久中文字幕一级| 免费在线观看日本一区| 黄色女人牲交| 精品欧美一区二区三区在线| 色在线成人网| 精品欧美一区二区三区在线| 日韩欧美国产一区二区入口| 欧美日韩一级在线毛片| 不卡一级毛片| 中国美女看黄片| 欧美不卡视频在线免费观看 | 熟女少妇亚洲综合色aaa.| 日韩 欧美 亚洲 中文字幕| 99久久精品国产亚洲精品| 好男人电影高清在线观看| 久久久国产成人免费| 黄色 视频免费看| 日本一二三区视频观看| 国产伦一二天堂av在线观看| 天天添夜夜摸| netflix在线观看网站| 久久精品影院6| 中文在线观看免费www的网站 | 亚洲国产精品久久男人天堂| 色在线成人网| 成年女人毛片免费观看观看9| 波多野结衣高清无吗| 国产成人av教育| 真人做人爱边吃奶动态| 国产91精品成人一区二区三区| 成人18禁在线播放| avwww免费| 97超级碰碰碰精品色视频在线观看| 老司机靠b影院| 桃红色精品国产亚洲av| 国产黄片美女视频| 岛国在线免费视频观看| 国产亚洲精品一区二区www| 国产一区二区三区在线臀色熟女| 麻豆一二三区av精品| 国产一区二区激情短视频| 脱女人内裤的视频| 天堂√8在线中文| 欧美乱码精品一区二区三区| a级毛片在线看网站| 午夜视频精品福利| 成熟少妇高潮喷水视频| 国产欧美日韩精品亚洲av| 国产成+人综合+亚洲专区| 亚洲精华国产精华精| 国产精品综合久久久久久久免费| 91麻豆av在线| 亚洲男人天堂网一区| 韩国av一区二区三区四区| 麻豆av在线久日| 精品电影一区二区在线| 国产激情久久老熟女| 国产av不卡久久| 不卡av一区二区三区| 禁无遮挡网站| 无遮挡黄片免费观看| 国产aⅴ精品一区二区三区波| 淫秽高清视频在线观看| 久热爱精品视频在线9| 亚洲人与动物交配视频| 久久久国产欧美日韩av| 国产在线精品亚洲第一网站| 欧美日韩乱码在线| 999久久久精品免费观看国产| 欧美丝袜亚洲另类 | 亚洲欧美一区二区三区黑人| 久久久久久免费高清国产稀缺| 亚洲免费av在线视频| 宅男免费午夜| 国产伦在线观看视频一区| 俺也久久电影网| 人人妻,人人澡人人爽秒播| 国产成人啪精品午夜网站| 久久香蕉国产精品| 亚洲最大成人中文| 1024手机看黄色片| 午夜福利欧美成人| 脱女人内裤的视频| 色播亚洲综合网| 亚洲全国av大片| 久久伊人香网站| 久久香蕉国产精品| 欧美精品亚洲一区二区| 午夜老司机福利片| 欧美午夜高清在线| 搡老熟女国产l中国老女人| 国产成人影院久久av| 国产三级在线视频| 最近最新中文字幕大全电影3| 叶爱在线成人免费视频播放| 动漫黄色视频在线观看| 精品久久蜜臀av无| 亚洲熟妇中文字幕五十中出| 久久久久国产精品人妻aⅴ院| 国产真实乱freesex| 日韩精品免费视频一区二区三区| 天堂av国产一区二区熟女人妻 | 免费在线观看视频国产中文字幕亚洲| 国产亚洲av嫩草精品影院| 男女午夜视频在线观看| 淫妇啪啪啪对白视频| 欧美成人性av电影在线观看| 一区二区三区高清视频在线| 亚洲国产精品久久男人天堂| 国产成年人精品一区二区| 天堂av国产一区二区熟女人妻 | 大型黄色视频在线免费观看| 欧美日韩福利视频一区二区| 日韩精品青青久久久久久| 日本撒尿小便嘘嘘汇集6| 欧美av亚洲av综合av国产av| 欧美另类亚洲清纯唯美| 国产精品久久久久久精品电影| 最新美女视频免费是黄的| АⅤ资源中文在线天堂| 午夜日韩欧美国产| 亚洲精品中文字幕一二三四区| 又紧又爽又黄一区二区| 国产亚洲av嫩草精品影院| 久久久久久人人人人人| 99热这里只有精品一区 | 国产高清有码在线观看视频 | 香蕉国产在线看| 美女午夜性视频免费| 99国产精品99久久久久| 国产成人影院久久av| 国产成人系列免费观看| 黄色毛片三级朝国网站| 一边摸一边抽搐一进一小说| 亚洲美女黄片视频| 国产精品电影一区二区三区| 麻豆av在线久日| 亚洲成人精品中文字幕电影| 99精品在免费线老司机午夜| 中国美女看黄片| 精品久久久久久,| 日日爽夜夜爽网站| www.www免费av| 美女免费视频网站| 成人精品一区二区免费| 一本久久中文字幕| 国产三级在线视频| 成人18禁在线播放| 亚洲国产欧洲综合997久久,| 黄色丝袜av网址大全| 久久香蕉精品热| 久久午夜亚洲精品久久| 久久人人精品亚洲av| 高清在线国产一区| 麻豆国产97在线/欧美 | 啦啦啦免费观看视频1| 精品久久久久久久久久免费视频| 黄色视频,在线免费观看| 在线观看66精品国产| 国产成人一区二区三区免费视频网站| 国产又色又爽无遮挡免费看| 高潮久久久久久久久久久不卡| 亚洲九九香蕉| 亚洲熟女毛片儿| 老司机午夜十八禁免费视频| 国产精品香港三级国产av潘金莲| 国产精品亚洲av一区麻豆| 18禁国产床啪视频网站|