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

    Four Wheel Independent Drive Electric Vehicle Lateral Stability Control Strategy

    2020-11-05 09:37:22YantaoTianXuanhaoCaoXiaoyuWangandYanboZhao
    IEEE/CAA Journal of Automatica Sinica 2020年6期

    Yantao Tian, Xuanhao Cao, Xiaoyu Wang, and Yanbo Zhao

    Abstract—In this paper, a kind of lateral stability control strategy is put forward about the four wheel independent drive electric vehicle. The design of control system adopts hierarchical structure. Unlike the previous control strategy, this paper introduces a method which is the combination of sliding mode control and optimal allocation algorithm. According to the driver’s operation commands (steering angle and speed), the steady state responses of the sideslip angle and yaw rate are obtained. Based on this, the reference model is built. Upper controller adopts the sliding mode control principle to obtain the desired yawing moment demand. Lower controller is designed to satisfy the desired yawing moment demand by optimal allocation of the tire longitudinal forces. Firstly, the optimization goal is built to minimize the actuator cost. Secondly, the weighted least-square method is used to design the tire longitudinal forces optimization distribution strategy under the constraint conditions of actuator and the friction oval. Beyond that, when the optimal allocation algorithm is not applied, a method of axial load ratio distribution is adopted. Finally, CarSim associated with Simulink simulation experiments are designed under the conditions of different velocities and different pavements. The simulation results show that the control strategy designed in this paper has a good following effect comparing with the reference model and the sideslip angle β is controlled within a small rang at the same time.Beyond that, based on the optimal distribution mode, the electromagnetic torque phase of each wheel can follow the trend of the vertical force of the tire, which shows the effectiveness of the optimal distribution algorithm.

    I. Introduction

    THE road traffic accidents occur frequently in recent years,much attention has been paid on the vehicle active safety studies [1]–[3]. Electric vehicle lateral stability, as an important part of vehicle active safety, has obtained the extensive and in-depth research by domestic and foreign scholars [4]–[8].

    Tabtiet al. put forward a new control method for the independent drive electric vehicle stability system, which described the vehicle dynamics as a linear hybrid system, and designed mixed model predictive controller. The experiment verifies that the controller makes the yaw angular velocity well follow the reference, ensures the vehicle driving safety and improves the vehicle stability [9]. Oharaet al. put forward a strategy of stability control based on the estimated information, the model reference control method is used for the composition of the controller, and this controller is able to make the yaw rate follow the reference yaw rate robustly.Numerical simulation is implemented with the consideration of a gust of wind and vehicle modeling error. Also, a vehicle experiment on the road is conducted [10]. Canaleet al.designed a reference generator to improve vehicle handling and provide the desired yaw rate value. After that, a closed loop controller was designed using the second-order sliding mode (SOSM) methodology to guarantee robust stability in front of disturbances and model uncertainties [11]. Peeie proposed a direct yaw moment control (DYC) using the front tire steer angle to increase the stability of the small electric vehicle with HMBS, and a kind of state feedback controller was presented. But several conditions must be met in the simulation to increase the stability of the vehicle: 1) The velocity of the vehicle must be constant; 2) The changes of the lateral velocity must be small [12]. Nasriet al.demonstrated the possibility of an improved four wheel vehicle stability which uses four independent drive wheels for motion by using the sliding mode control [13]. Namet al. put forward a method to estimate vehicle side-slip angle based on the lateral tire force sensor and improve the wheel motor drive electric vehicle stability. The extended Kalman filter (EKF) state observer is designed to estimate vehicle side-slip angle. In addition, a vehicle stability controller was designed, and a simulation experiment was implemented to verify the effectiveness [14]. Nahidiet al. proposed a modular control structure for optimal integrated vehicle longitudinal and lateral stability control. The optimal torque distribution adopted the HCC strategy, which minimizes these CG discrepancies by considering the errors between the target and actual vehicle CG forces and yaw moment [15]. Zhanget al.proposed a kind of hierarchical framework to reduce the complexity for design and analysis. A high-level controller was designed to generate the desired longitudinal dynamics,and at the lower level, an adaptive sliding-mode controller was designed. Numerical simulations demonstrated the robustness of the controller [16]. Xiaet al.presented a tracking control method for the lateral motion of an autonomous land vehicle (ALV), which is based on active disturbance rejection control (ADRC) scheme and differential flatness theory [17]. Khosravaniet al. proposed a new torque distribution method which can contain multi-objective functions such as force tracking at CG, minimum actuator adjustments, etc. And several physical constraints such as constraint for avoiding wheel overspin/lock were incorporated into the allocation task. Numerical simulations demonstrated that it could improve the vehicle stability [18]. Rafailaet al.proposed a nonlinear model based predictive control strategy for the autonomous control of the steering system for ground vehicles. According to the simulation results, this approach could provide good performances in practical use [19].

    To sum up, the study of stability control system of the vehicle mainly concentrates on the upper controller development and design at present. The allocation algorithm of the lower controller is widely used by axial load distribution ratio method, but the utilization rate of each wheel is not high, which affects the vehicle’s stability and safety.Therefore, in this paper, firstly, the control allocation algorithm is studied as the research key, which is used to distribute longitudinal force of each tire reasonably. The weighted least-square method is adopted to design the tire longitudinal forces optimization distribution strategy. Beyond that, when the optimal allocation algorithm does not apply, the method of axial load ratio distribution is adopted. Secondly, a new hierarchical structure of lateral stability controller is designed which adopts the method of sliding mode control combined with optimal allocation algorithm.

    The purpose of this paper is to design a kind of lateral stability control strategy and improve the vehicle stability in more critical conditions. And the control strategy is double layered system. Upper controller adopts the sliding mode control principle to obtain the desired yawing moment demand. Lower controller is designed to follow the desired yawing moment demand by optimal allocation of the tire longitudinal forces. Finally,CarSimassociated withSimulinksimulation experiments are designed under the conditions of different velocities and different pavements.

    The main contributions of this paper can be described as two points: 1) A kind of vehicle tire longitudinal forces optimization distribution strategy is designed, and the optimization goal is built to minimize the actuators cost.According to this approach, the utilization rate of each wheel can be increased, and vehicle stability and safety can be improved. 2) The optimal allocation algorithm will not apply when the vehicle is travelling on the low adhesion pavement,because when the coefficient of friction is very low, the longitudinal force constraint may appear empty set. So, in this case, the distribution mode contains the switch of the optimal allocation and the axle load distribution. According to this approach, the effectiveness of actuator can be guaranteed.

    This paper is organized as follows. In Section II, a lateral dynamics model is designed based on the vehicle yawing plane model which is used to the design of upper controller. In Section III, firstly, the lateral stability system control strategy and basic frame are described. Secondly, a kind of vehicle reference model is calculated. After that, the sliding mode controller is designed to control the vehicle tracking vehicle reference model. In Section IV, two different control distribution methods are designed, the first one uses the quadratic programming method, and the other one is the axle load proportional allocation algorithm. In Section V,Carsim/Simulinksimulation experiments are conducted under the condition of different velocities and different pavements.According to the lateral stability control system simulation experiments, the validity of the control strategy and the distribution algorithm designed are demonstrated.

    II. Vehicle Dynamics Modeling

    The vehicle system is a kind of high degree of freedom and strong nonlinear system which are difficult for control design.Therefore, in order to simplify the vehicle model, it is necessary to make the following assumptions [20]: 1) the vehicle runs on a horizontal road; 2) the impact of vehicle vertical load transfer is not considered; 3) the two front wheel angles are equal and the front and rear axle wheelbases are approximately equal; 4) the vehicle tire and suspension remain relatively vertical state to the ground.

    A. Lateral Dynamics Model

    Lateral dynamics model of this paper is designed based on the vehicle yawing plane model. The main difference from the traditional vehicle dynamic model is that this paper takes the direct yaw moment as an additional input variable, which is based on the independent output torque of each wheel [21].The vehicle’s yawing plane is shown in Fig. 1, and the parameters of vehicle dynamics are given in Table I. The vehicle lateral dynamics model and the yawing moment balance equation at the barycenter position are established as follows:

    Fig. 1. Vehicle yaw plane four wheel model.

    TABLE I Parameters of Vehicle Dynam ics

    Fig.2.Vehicle linear bicycle dynamicmodel.

    According to the above vehicle linear bicycle dynam ic model,the simplified vehicle lateral and the yaw dynam ics equation are described as follows:

    B.Wheel Dynamics Model

    Fig.3.Thewheel torque balance diagram.

    2)The Wheel Vertical Load

    Ignoring the vehicle suspension system,the vehicle can be taken as a rigid body.According to the model vehicle dynam ic model and the rigid mechanics principles,the calculation formulas of wheel vertical load are shown as follows[20]:

    Fig.4.Theelectric vehicle lateral stability control system structure diagram.

    C.Tire Model

    III.Direct ly Yaw ing Moment Con t rol ler

    A.The Lateral Stability System Control Strategy

    Lateral stability control goal is to enhance the vehicle steering stability and ensure the vehicle to track the expected trajectory at the same time.Thewhole system isdivided into a two-layer structure, the upper controller is designed for the yaw rate and side-slip angle tracking control.Based on the yaw moment requirement and considering constraints at the same time,a lower controller is designed to optimize the motor torque distribution and ensure the electric vehicle stability.Fig.4 is the lateral stability control system structure adopted in thisarticle.

    B. Reference Model

    According to the sideslip angleand the yaw rate steady-state response,which are produced by the driver demands(the steering angle and speed), the vehicle reference model can be built. Normally,the vehicle steady state responseβ˙=γ˙=0 is used as the referencemodel.According to(4),(5)and the tire cornering stiffness values,the desired vehicle sideslip angle and yaw rateare defined as follows:

    Generally speaking,the maximum value of the stability vehicle sideslip angle changes w ith the road adhesion conditions.On the good road condition,the sideslip angle stability value is greater than thaton theworse road adhesion conditions[25].According to this experience,this article selects road adhesion coefficient to lim it the valueof

    To sum up,the desired vehicle sideslip angle is revised as follows:

    C.Sliding Mode Controller

    For the vehicle stability control system, the controleffectof general linear controllers usually can not satisfy the requirementof improving the vehicle steering control stability and safety.According to this,this paper chooses the sliding mode controller to control the vehicle to track the desired yaw rate and side-slip angle,and the sliding mode surface[26]is chosen as follows:

    IV.Con t rol Torque A l loca tion A lgor ithm

    The controlof electric vehicles has been conducted inmany research reports,such as using linear quadratic regulator or state feedback controller to maintain the vehicle driving stability.However,in order to make full use of the vehicle’s performance potential,a well-designed control distribution method is required.

    A.Optimization Goal

    The vehicle total tire longitudinal force and total yaw moment equationsareas follows:

    The equation can use the quadratic programmingmethod to solve.Based on the quadratic programm ing theory,firstly,the optimal objective function and the constraint conditions should be determined.In this paper,we choose the sum of squaresof four wheel longitudinal forcesand lateral forcesdivided by the square of the product of the vertical loadand the friction coefficient as the optimization objective function,which can be used to express the vehicle stability margin.Generally, the higher of the tire utilization,the lower the vehicle stabilitymargin When the objective function value is close to 1, the adhesion between tire and road w illbe close to its lim it,and the vehicle is in a critical state of losing stability at this time.So in this paper,the optimization goal is to minimize the objective function and ensure the maximization of the vehicle stability margin[27].The objective function isexpressed as[28]

    As(34)shows,in this paper, the optim ization objective function is the variantof theminimum sum of squares of the tire utilization.The lateral force is ignored,and the tire utilization minimum is produced by the selection of the longitudinal forcewhich is used in the objective function.To be clear,in the actual situation, the longitudinal and lateral forces cannot be completely independent.We can think about it in a different way,the optimization objective function adopted in this paper can minim ize the actuator control cost,and ensure the vehicle driving stability.The optim ization goal is to distribute the longitudinal force and torque,which makes the motor to retain a bigger margin[29]in sudden unstable situation and makes the vehicle havegreater ability to own the initiative to maintain or change the vehicle state and enhance the vehicle safety and stability.Italso has certain significance for thestability of control systems.

    B.Constraint Conditions

    C.The Optimization Allocation Algorithm

    According to(31),(34),and(36), the follow ing quadratic programm ing standard isobtained:

    D.The Axle Load Proportional Allocation Algorithm

    The optimal allocation algorithm w ill be no longer applicablewhen the vehicle is running on the pavement w ith poor adhesion,because the coefficient of friction is very low,the longitudinal force constraint(39) may result in an empty set.In this case,for the vehicle stability control system,the vehicle tire longitudinal force distribution adopts themethod of axial load ratio distribution.The total driving force and the total yaw ing moment have been distributed according to the proportion of front and rear axle load,axle load valuesas

    So each wheel longitudinal force by axial load proportion distribution can be obtained as follows:

    V.La tera l Stabil ity Con t rol System Simu la tion Resu l ts

    This section is the study and analysis of the lateral stability control system simulation experiments.Compared w ith the open-loop experiment,the closed-loop driver model can better reflect the driver’s real driving process.In this simulation experiment, weadopt double line conditions in different speed and road surface conditions, respectively.The parameters of theCarsimvehiclemodel are shown in Table II.

    TABLE II Vehic le Pa rameters

    A.The Simulation Results With Initial Speed 60 km/h Under High Adhesion and Low Adhesion Road Surface Condition

    Figs.5 and 6 show the simulation results w ith the vehicle initialspeed 60 km/h when the road friction coefficient is 0.9.Fig.6 shows the vehicle yaw rate and sideslip angle tracking errors.And according to Fig.6,the tracking errors based on the optimal distribution mode is smaller than the tracking errors based on the average distribution mode.Figs.5(a),5(b),and Fig.6 show that the valuesof the vehicle yaw rate and the vehicle sideslip angle have good follow ing performances for the reference values,which demonstrates that the designed controller has good control results.Fig.5(c)shows the simulation result of the vehicle four wheel vertical force;Fig.5(d)shows the simulation results of the yaw ing moment under the optim ization distribution mode,the average distribution mode and the PI controller,respectively.Figs.5(e)and 5(f)show the electromagnetic torques under the optimization distribution mode and the average distribution mode,respectively.According to Figs.5(e)and 5(f),it can be seen that the optim ization distribution mode and the average distribution mode both have good torque output results.However,according to contrast Figs.5(e)and 5(f)w ith 5(c),it can be seen that in the optimal distribution mode,the electromagnetic torque phase of each wheel can follow the trend of the vertical force of the tiremore than in the average distribution mode.Generally,under the same road surface condition,the larger the vertical force of the wheel is,the greater longitudinal driving force can be produced,so the larger longitudinal force output isalso assigned in the optimal distributionmode.This results also show the effectiveness of theoptimal distribution algorithm and itsgood performance.

    Figs.7 and 8 show the simulation results w ith the vehicle initial speed 60 km/hwhen the road friction coefficient is0.4.According to Figs.7(a),7(b),and Fig.8,the designed control system has good follow ing effect for the desired yaw rate and sideslip angle.Fig.7(c)is the vehicle four wheel vertical force;From Fig.7(d),we can see that the yaw moment has sudden changes in neighborhood of 3 s to 4 s and 7 s to 8 s,and the optim ization distribution mode needs smaller yaw moment,so the vehicle has more ability to deal w ith the sudden unstable situation,which enhances the vehicle stability to a certain extent.It shows that the yaw moment control resultof the optimization distribution mode is superior to that of the average distribution mode.According to contrast Figs.7(e)and 7(f)w ith Fig.7(c), the torque output of the optimal distribution mode can also follow the trend of the vertical force of the tire.In this respect,the optim ization distribution mode is still superior to the average distribution mode.Figs.7(e)and 7(f)show that, when the pavement condition is poor,the vehicle right front wheel electromagnetic torque suddenly drops to zero in the average distribution mode in the 2–4 s and the 6–8 s, which indicates that the right front wheel is locked,and it is harm ful to the stability of the vehicle and the com fort of driving.

    From Figs.5–7,under the same speed,vehicle stability control system based on the optimization distribution mode can have good control effect no matter the tire-road friction condition is good or poor.So according to the simulation results,the control result based on the optim ization distribution mode is better than that of the average distribution mode,which improves the vehicle stability on the low adhesion road.

    B.The Simulation Results With Initial Speed 120 km/h Under High Adhesion and Low Adhesion Road Surface Condition

    Fig.5.Resultsw ith initial velocity 60 and tire-road friction coefficient 0.9.

    Fig.6.Yaw rate and sideslip angle tracking errors.

    Figs.9 and 10 show the simulation results w ith the vehicle initial speed 120 km/h when the road friction coefficient is 0.9.Fig.10 shows the vehicle yaw rate and sideslip angle tracking errors.And according to Fig.10,the tracking errors based on the optimal distribution mode are smaller than the tracking errors based on the average distribution mode.Fig.9(c)denotes the vehicle four wheel vertical force;Fig.9(d)is the yaw moment figure of the optim ization distribution mode,average distribution mode and the PI control.From Fig.9(d), we can see that the required yaw ing moment based on the optimaldistribution algorithm isslightly less than thatof the average distributionmode no matterwhat controller isadopted (PIcontroller or sliding mode controller),which shows that the optimization distribution mode has better controleffectand helps to improve the vehicle stability.Figs.9(e)and 9(f)denote the electromagnetic torque inputs under the optim ization distribution mode and average distribution mode,respectively.The optimization distribution mode contains the sw itch of the optimal allocation and the axle load distribution,so there are some larger amplitude changes but it can still follow the trend of the vertical force of the tire and maintain the stability of the vehicle.

    Fig.7.Resultsw ith initial velocity 60 and tire-road friction coefficient 0.4.

    Fig.8.Yaw rate and sideslip angle tracking errors.

    Figs.11 and 12 show the simulation resultsw ith the vehicle initial speed 120 km/h when the road friction coefficient is 0.4.From Figs.11(a)and 11(b), we can see that the vehicle control effect is inferior to the effect of the high tire-road friction coefficient but yaw rate and vehicle sideslip angle can still follow the changes of the references.Fig.11(c)is the vehicle four wheel vertical forces.Fig.11(d)denotes the vehicle yaw ing moment control result figure.From Fig.11(d),it can be seen that yaw ing moment provided by the wheel is smaller and the vehicle has better stability and security no matter what controller is adopted(PI controller or sliding mode controller).So the optim ization distribution has obvious advantages over the average distribution.According to contrast Figs.11(e)and 11(f)w ith 11(c)the torque output of the optimal distributionmode can also follow the trend of the vertical forceof the tire when the vehicle speed ishigh and the road friction coefficient is low.Finally,the optim ization distributionmode by sw itching the optimal allocation and the axle load distribution mode has better ability to maintain vehicle stability.

    Fig.9.Resultsw ith initial velocity 120 and tire-road friction coefficient 0.9.

    Fig.10. Yaw rateand sideslip angle tracking errors.

    According to contrast Figs.5(b)w ith 11(b),the optimization control algorithm proposed in this paper can still have a good control effect when the sideslip angle is larger and vehicle speed is higher.Under the constraints of road conditions and motor torque, the vehicle stability and security are improved based on the reasonable distribution of the tire longitudinal force.Finally,the simulation experiments verify the effectiveness of the proposed control strategy and the distribution algorithm.

    VI.Conc lusion

    In this paper,a kind of electric vehicle lateral stability control strategy is put forward.It can be concluded as follows.

    1)The hierarchical control structure is adopted,the upper controller and the lower controller are designed,respectively.

    2)The upper controller is themotion controller.In order to guarantee the robustness and stability of the system,the sliding mode controller is applied.The vehicle sideslip angleand yaw rateare selected as control variables to calculate the yaw ing moment needed to stabilize thevehicle.

    Fig.11.Resultsw ith initial velocity 120 and tire-road friction coefficient 0.4.

    Fig.12. Yaw rateand sideslip angle tracking errors.

    3)The lower controller adopts the sum of the squaresof tire utilization as the stability control target.Considering the road adhesion condition andmotor torque constraints,theweighted least squares method has been adopted to optim ize the tire longitudinal force distribution,which is used to guarantee the vehicle stability.

    Finally,CarSimassociated w ithSimulinksimulation experiments are designed under the condition of different velocities and different pavements.The experimental results show the control strategy designed in this paper has a good follow effect in comparison w ith the referencemodel and the sideslip angleis controlled in a small range at the same time.In addition, based on the optimal distributionmode,the electromagnetic torque phase of each wheel can follow the trend of the vertical force of the tire.Generally,under the same road surface condition,the larger the vertical force of the wheel,the greater the longitudinal driving force can be produced,so the larger longitudinal force output is also assigned in the optimal distribution mode.The optimal distribution algorithm designed in this paper can also prevent wheels from locking under poor pavements and needs smaller yaw ing moments,which shows the effectiveness of the optimal distribution algorithm.Beyond that,when the optimal allocation algorithm is no longer applied,a method of axial load ratio distribution is adopted.Sometimes the optim ization distribution mode contains the sw itch of theoptimalallocation and the axle load distribution,although there are some larger amplitude changes,it can still maintain the stability of the vehicle.

    In the next work, the further research and effort can be devoted to the follow ing two points:1)Online estimation of the road friction coefficient.In this paper,the road friction coefficient is a constant,and the effectivenessof the system is proved under different pavements whose road friction coefficients are fixed.But in order to carry out real vehicle experiments,it is necessary to estimate the road friction coefficient on line.2)Based on the online estimation of the road friction coefficient,some real vehicle experiments can be carried at the next work.And different experimental conditions can be considered,for example snake-like condition and double lane change test.

    久热久热在线精品观看| 成人漫画全彩无遮挡| 十分钟在线观看高清视频www | 国产精品一区二区在线不卡| 一本色道久久久久久精品综合| 日本wwww免费看| 亚洲怡红院男人天堂| www.色视频.com| 国产欧美日韩综合在线一区二区 | 久久久久久久精品精品| 天美传媒精品一区二区| 欧美日韩国产mv在线观看视频| 美女内射精品一级片tv| 久久97久久精品| 在线观看免费高清a一片| 一本久久精品| 大香蕉97超碰在线| 日韩 亚洲 欧美在线| 午夜激情福利司机影院| 妹子高潮喷水视频| 女的被弄到高潮叫床怎么办| 亚洲不卡免费看| 老熟女久久久| 老女人水多毛片| 国产色婷婷99| 日韩电影二区| 在线播放无遮挡| 麻豆成人av视频| 亚洲中文av在线| 国产高清不卡午夜福利| 亚洲欧美精品自产自拍| 最新的欧美精品一区二区| 欧美日韩视频精品一区| 日韩精品有码人妻一区| 亚洲一级一片aⅴ在线观看| 九色成人免费人妻av| 欧美国产精品一级二级三级 | 亚洲人成网站在线观看播放| 精品少妇黑人巨大在线播放| 久久女婷五月综合色啪小说| 亚洲欧美一区二区三区黑人 | 中国三级夫妇交换| 中国国产av一级| 亚洲欧美日韩卡通动漫| 成人无遮挡网站| 国产淫片久久久久久久久| 在线 av 中文字幕| 少妇人妻精品综合一区二区| 国产精品久久久久久精品古装| 亚洲不卡免费看| 国产免费又黄又爽又色| 97超碰精品成人国产| 国产在线视频一区二区| 免费观看在线日韩| 2022亚洲国产成人精品| 九九久久精品国产亚洲av麻豆| 丰满人妻一区二区三区视频av| 丰满乱子伦码专区| 欧美日韩视频高清一区二区三区二| 中文精品一卡2卡3卡4更新| 黑人猛操日本美女一级片| 一级二级三级毛片免费看| 久久韩国三级中文字幕| 2021少妇久久久久久久久久久| 国产中年淑女户外野战色| 国产精品久久久久久久久免| 日本91视频免费播放| 日韩欧美一区视频在线观看 | 久久午夜综合久久蜜桃| 国产免费视频播放在线视频| 老司机亚洲免费影院| 色视频www国产| 亚洲不卡免费看| 午夜免费观看性视频| 亚洲在久久综合| 看十八女毛片水多多多| 男女免费视频国产| 欧美精品一区二区免费开放| 色婷婷av一区二区三区视频| 日本猛色少妇xxxxx猛交久久| 久久午夜综合久久蜜桃| av卡一久久| 熟女人妻精品中文字幕| 视频区图区小说| 美女国产视频在线观看| 日韩av免费高清视频| 亚洲情色 制服丝袜| 国语对白做爰xxxⅹ性视频网站| 97精品久久久久久久久久精品| 纯流量卡能插随身wifi吗| 夜夜爽夜夜爽视频| 亚洲性久久影院| 99视频精品全部免费 在线| 欧美亚洲 丝袜 人妻 在线| 日韩大片免费观看网站| 97在线视频观看| 最近中文字幕高清免费大全6| 大香蕉97超碰在线| 另类精品久久| 日本91视频免费播放| 伊人久久精品亚洲午夜| 亚洲电影在线观看av| 午夜日本视频在线| 亚洲欧美成人综合另类久久久| 国产一区亚洲一区在线观看| 9色porny在线观看| 午夜91福利影院| 国产亚洲一区二区精品| 99久久人妻综合| 99久久精品一区二区三区| 亚洲在久久综合| 美女脱内裤让男人舔精品视频| 两个人的视频大全免费| 久久国产乱子免费精品| 97超视频在线观看视频| 日韩三级伦理在线观看| 国产黄片视频在线免费观看| 一本一本综合久久| 欧美3d第一页| 国产一区有黄有色的免费视频| 国产成人午夜福利电影在线观看| 欧美人与善性xxx| 制服丝袜香蕉在线| 深夜a级毛片| 国产精品国产三级国产av玫瑰| 激情五月婷婷亚洲| 久久韩国三级中文字幕| 简卡轻食公司| 自线自在国产av| 久久久精品免费免费高清| 国产白丝娇喘喷水9色精品| av网站免费在线观看视频| 久久久久精品久久久久真实原创| 国产精品免费大片| 综合色丁香网| 丰满少妇做爰视频| 韩国av在线不卡| 极品少妇高潮喷水抽搐| 国产成人午夜福利电影在线观看| 国产成人精品福利久久| 国产精品久久久久久久电影| 男女啪啪激烈高潮av片| 国产欧美日韩综合在线一区二区 | 美女福利国产在线| 曰老女人黄片| 国产爽快片一区二区三区| 日本欧美视频一区| 亚洲在久久综合| 肉色欧美久久久久久久蜜桃| 免费观看a级毛片全部| xxx大片免费视频| 国产黄片视频在线免费观看| 国产老妇伦熟女老妇高清| 狠狠精品人妻久久久久久综合| 青春草视频在线免费观看| 国产在线免费精品| 在线观看av片永久免费下载| 69精品国产乱码久久久| 久久ye,这里只有精品| 最近手机中文字幕大全| 麻豆精品久久久久久蜜桃| 高清毛片免费看| 久久av网站| 亚洲欧美精品自产自拍| 欧美日韩精品成人综合77777| 日本猛色少妇xxxxx猛交久久| 91精品伊人久久大香线蕉| a级毛色黄片| 黄色欧美视频在线观看| 大香蕉久久网| 日韩一本色道免费dvd| 色婷婷av一区二区三区视频| 日韩成人av中文字幕在线观看| av女优亚洲男人天堂| 18禁动态无遮挡网站| 国产av一区二区精品久久| 人人妻人人澡人人爽人人夜夜| av在线观看视频网站免费| 日韩三级伦理在线观看| 高清不卡的av网站| 国产美女午夜福利| 狠狠精品人妻久久久久久综合| 另类亚洲欧美激情| 日本免费在线观看一区| 五月伊人婷婷丁香| 国产视频首页在线观看| 亚洲va在线va天堂va国产| 青春草国产在线视频| 亚洲第一av免费看| 成人亚洲精品一区在线观看| 欧美亚洲 丝袜 人妻 在线| 国产成人免费观看mmmm| 欧美日韩视频高清一区二区三区二| 成年人午夜在线观看视频| 午夜精品国产一区二区电影| 人人妻人人看人人澡| 午夜福利网站1000一区二区三区| 欧美bdsm另类| 国产日韩一区二区三区精品不卡 | 少妇人妻久久综合中文| 久久人人爽av亚洲精品天堂| 亚洲综合精品二区| 亚洲精品456在线播放app| 两个人免费观看高清视频 | 精品国产国语对白av| 国产69精品久久久久777片| 男女无遮挡免费网站观看| 伦理电影大哥的女人| 婷婷色麻豆天堂久久| 日韩三级伦理在线观看| 国产免费视频播放在线视频| a级毛色黄片| 久久久精品免费免费高清| 午夜福利在线观看免费完整高清在| 亚洲av成人精品一二三区| 综合色丁香网| 久久久久国产精品人妻一区二区| 久久久a久久爽久久v久久| 看免费成人av毛片| 午夜日本视频在线| 午夜激情久久久久久久| 精品卡一卡二卡四卡免费| 欧美日韩亚洲高清精品| 久久久欧美国产精品| 日本-黄色视频高清免费观看| 精品久久久久久久久亚洲| 色哟哟·www| 大香蕉97超碰在线| 高清毛片免费看| 人妻一区二区av| 天天躁夜夜躁狠狠久久av| 最近最新中文字幕免费大全7| 青青草视频在线视频观看| 精品一区二区三区视频在线| 国产亚洲av片在线观看秒播厂| 国产欧美亚洲国产| 蜜臀久久99精品久久宅男| 啦啦啦在线观看免费高清www| 亚洲成人手机| 美女视频免费永久观看网站| 国产av一区二区精品久久| 日本黄大片高清| 成人国产av品久久久| 久久久久久伊人网av| 中国三级夫妇交换| 美女xxoo啪啪120秒动态图| 国产片特级美女逼逼视频| 国产高清三级在线| 欧美区成人在线视频| 国产成人91sexporn| 91在线精品国自产拍蜜月| 99热这里只有是精品在线观看| 中文字幕久久专区| 欧美日韩一区二区视频在线观看视频在线| 亚洲国产精品专区欧美| 国产免费视频播放在线视频| 在线亚洲精品国产二区图片欧美 | 看十八女毛片水多多多| 七月丁香在线播放| 亚洲国产最新在线播放| 精品99又大又爽又粗少妇毛片| 国产精品一区二区三区四区免费观看| 久久国产乱子免费精品| 久久青草综合色| 久久久久久人妻| 亚洲av二区三区四区| 亚洲欧美成人综合另类久久久| 啦啦啦视频在线资源免费观看| 中文字幕免费在线视频6| 亚洲国产日韩一区二区| 国产欧美日韩综合在线一区二区 | 嘟嘟电影网在线观看| 国产精品.久久久| 亚洲怡红院男人天堂| 久久 成人 亚洲| 蜜桃久久精品国产亚洲av| 欧美日韩综合久久久久久| 少妇人妻久久综合中文| 亚洲高清免费不卡视频| 亚洲精品一区蜜桃| 最近中文字幕2019免费版| 99热网站在线观看| 亚洲四区av| 中文字幕久久专区| 美女大奶头黄色视频| 亚洲精品国产成人久久av| 99热全是精品| 亚洲在久久综合| 精华霜和精华液先用哪个| av国产精品久久久久影院| 欧美 日韩 精品 国产| 伊人久久精品亚洲午夜| 成人黄色视频免费在线看| 亚洲成人手机| 亚洲国产精品国产精品| 精品熟女少妇av免费看| 又粗又硬又长又爽又黄的视频| 国产淫片久久久久久久久| 精品一品国产午夜福利视频| 成人毛片a级毛片在线播放| 欧美精品国产亚洲| 成年人午夜在线观看视频| 久久婷婷青草| 午夜日本视频在线| 我要看日韩黄色一级片| 韩国高清视频一区二区三区| 国产伦理片在线播放av一区| 国产精品成人在线| 国产综合精华液| 久热久热在线精品观看| h日本视频在线播放| 最新中文字幕久久久久| 六月丁香七月| 国产色爽女视频免费观看| 精品熟女少妇av免费看| 国产亚洲91精品色在线| 免费在线观看成人毛片| 亚洲国产成人一精品久久久| 人人妻人人澡人人爽人人夜夜| 最新的欧美精品一区二区| 免费大片18禁| 亚洲人与动物交配视频| 久久这里有精品视频免费| 国产免费一级a男人的天堂| 色哟哟·www| 熟女人妻精品中文字幕| 国产成人午夜福利电影在线观看| 免费观看a级毛片全部| 成人无遮挡网站| 老女人水多毛片| 国内少妇人妻偷人精品xxx网站| 亚洲国产精品一区三区| 如何舔出高潮| 国产高清不卡午夜福利| 爱豆传媒免费全集在线观看| 久久久久国产精品人妻一区二区| 亚洲久久久国产精品| 桃花免费在线播放| 乱系列少妇在线播放| 中文资源天堂在线| 搡女人真爽免费视频火全软件| 人妻人人澡人人爽人人| 午夜免费男女啪啪视频观看| 国产 一区精品| 国产av国产精品国产| 中文在线观看免费www的网站| 91久久精品电影网| 国产精品女同一区二区软件| 亚洲美女搞黄在线观看| 亚洲欧美精品自产自拍| av天堂久久9| 亚洲国产成人一精品久久久| 精品国产一区二区三区久久久樱花| 少妇的逼好多水| 国产精品一二三区在线看| 性高湖久久久久久久久免费观看| 丰满人妻一区二区三区视频av| 亚洲人与动物交配视频| 亚洲在久久综合| 亚洲国产精品专区欧美| 亚洲欧美一区二区三区国产| 亚洲欧美日韩东京热| 亚洲av男天堂| 亚洲欧美一区二区三区国产| 少妇高潮的动态图| 国产一区亚洲一区在线观看| 亚洲欧美日韩另类电影网站| 少妇高潮的动态图| 亚洲精品色激情综合| 热re99久久精品国产66热6| 欧美最新免费一区二区三区| 国产精品免费大片| 国产无遮挡羞羞视频在线观看| 亚洲成人av在线免费| 大片电影免费在线观看免费| a级毛色黄片| 色网站视频免费| 免费av中文字幕在线| 国产精品伦人一区二区| 亚洲成人av在线免费| 国产国拍精品亚洲av在线观看| 国产69精品久久久久777片| 国产av码专区亚洲av| 日本猛色少妇xxxxx猛交久久| 青春草国产在线视频| 黄色一级大片看看| 一本色道久久久久久精品综合| 97精品久久久久久久久久精品| 日韩,欧美,国产一区二区三区| 成人国产麻豆网| 亚洲国产精品成人久久小说| 久久av网站| 久久久久久伊人网av| 简卡轻食公司| 99re6热这里在线精品视频| 七月丁香在线播放| 最近的中文字幕免费完整| 亚洲欧美日韩东京热| 丰满人妻一区二区三区视频av| 亚洲国产色片| 亚洲,欧美,日韩| 久久精品久久久久久噜噜老黄| 成人无遮挡网站| 日韩成人伦理影院| 一级片'在线观看视频| 青春草国产在线视频| 亚洲精品aⅴ在线观看| 亚洲va在线va天堂va国产| 欧美高清成人免费视频www| 日本午夜av视频| 一级黄片播放器| 插逼视频在线观看| 国产黄片视频在线免费观看| 男女无遮挡免费网站观看| 一区在线观看完整版| 在现免费观看毛片| www.av在线官网国产| 丰满乱子伦码专区| 国产深夜福利视频在线观看| 亚洲第一区二区三区不卡| 噜噜噜噜噜久久久久久91| 午夜精品国产一区二区电影| 日韩人妻高清精品专区| 国产av国产精品国产| 精品久久久噜噜| 日韩精品有码人妻一区| 一级黄片播放器| 欧美xxⅹ黑人| 久久国产亚洲av麻豆专区| 80岁老熟妇乱子伦牲交| 桃花免费在线播放| 亚洲av在线观看美女高潮| 中国三级夫妇交换| 精品一品国产午夜福利视频| 国产日韩一区二区三区精品不卡 | 九草在线视频观看| 精品久久久精品久久久| 纵有疾风起免费观看全集完整版| 国产乱来视频区| 纵有疾风起免费观看全集完整版| 麻豆乱淫一区二区| 人人妻人人澡人人看| av在线app专区| 亚洲精品色激情综合| 精品一品国产午夜福利视频| 人人妻人人添人人爽欧美一区卜| 成人亚洲精品一区在线观看| 日本欧美国产在线视频| 国产极品天堂在线| 欧美成人精品欧美一级黄| 国产精品秋霞免费鲁丝片| 久久久a久久爽久久v久久| 国国产精品蜜臀av免费| 亚洲性久久影院| 中文字幕av电影在线播放| 免费av不卡在线播放| 女性被躁到高潮视频| 又粗又硬又长又爽又黄的视频| 在线观看人妻少妇| 免费av中文字幕在线| 成年女人在线观看亚洲视频| 中文精品一卡2卡3卡4更新| 乱人伦中国视频| 国产伦理片在线播放av一区| 亚洲国产欧美在线一区| 亚洲国产精品成人久久小说| 精品酒店卫生间| 国产综合精华液| 一级毛片aaaaaa免费看小| 日本av免费视频播放| 下体分泌物呈黄色| 夜夜爽夜夜爽视频| 亚洲国产av新网站| 秋霞在线观看毛片| 国产免费福利视频在线观看| 在线观看免费高清a一片| 内射极品少妇av片p| 2018国产大陆天天弄谢| 欧美精品高潮呻吟av久久| 国产午夜精品久久久久久一区二区三区| 久久免费观看电影| 97超视频在线观看视频| 丝瓜视频免费看黄片| 久久久久精品久久久久真实原创| 亚洲av综合色区一区| 在线观看免费日韩欧美大片 | 久久影院123| 观看免费一级毛片| 国产成人免费观看mmmm| 久久精品国产亚洲av天美| 在线观看免费日韩欧美大片 | 亚洲图色成人| 极品人妻少妇av视频| 熟妇人妻不卡中文字幕| 在线观看免费日韩欧美大片 | 国产一级毛片在线| 国产极品天堂在线| a 毛片基地| 亚洲第一区二区三区不卡| 色哟哟·www| 亚洲三级黄色毛片| 七月丁香在线播放| 欧美 日韩 精品 国产| 成人黄色视频免费在线看| 婷婷色av中文字幕| 国产精品福利在线免费观看| 狂野欧美白嫩少妇大欣赏| 两个人的视频大全免费| 国产免费一区二区三区四区乱码| 国产在线视频一区二区| 成人18禁高潮啪啪吃奶动态图 | 麻豆成人午夜福利视频| 成人亚洲欧美一区二区av| 视频中文字幕在线观看| 80岁老熟妇乱子伦牲交| 成人国产av品久久久| 精品一区二区三区视频在线| 免费人成在线观看视频色| 3wmmmm亚洲av在线观看| 51国产日韩欧美| 交换朋友夫妻互换小说| 99九九在线精品视频 | 午夜日本视频在线| 日日爽夜夜爽网站| 久久久久久久久久成人| 欧美性感艳星| 国产精品一区www在线观看| 精品国产露脸久久av麻豆| 国产成人a∨麻豆精品| 在线播放无遮挡| 9色porny在线观看| 亚洲欧美成人精品一区二区| 久久影院123| 国产淫语在线视频| 久久久国产一区二区| 99热这里只有精品一区| 欧美日韩精品成人综合77777| 女的被弄到高潮叫床怎么办| 亚洲国产精品国产精品| 亚洲综合色惰| 欧美成人午夜免费资源| 人人妻人人澡人人看| 极品教师在线视频| 欧美精品一区二区免费开放| 亚洲精品456在线播放app| 日韩伦理黄色片| 欧美成人午夜免费资源| 女人精品久久久久毛片| av在线老鸭窝| 啦啦啦在线观看免费高清www| av在线老鸭窝| 亚洲精品色激情综合| 亚洲国产毛片av蜜桃av| 五月玫瑰六月丁香| 欧美一级a爱片免费观看看| 国产男人的电影天堂91| 国产色婷婷99| 亚洲精品国产成人久久av| 午夜福利,免费看| 啦啦啦中文免费视频观看日本| 精品人妻熟女毛片av久久网站| 欧美国产精品一级二级三级 | av又黄又爽大尺度在线免费看| 在线观看一区二区三区激情| 亚洲一级一片aⅴ在线观看| 涩涩av久久男人的天堂| 久久久久视频综合| 乱码一卡2卡4卡精品| 国产精品人妻久久久影院| 美女中出高潮动态图| √禁漫天堂资源中文www| 丝袜喷水一区| 国产男人的电影天堂91| 午夜影院在线不卡| 国产一区亚洲一区在线观看| 久久精品国产鲁丝片午夜精品| 亚洲欧洲国产日韩| 高清毛片免费看| 国产一区二区在线观看av| 老女人水多毛片| 一边亲一边摸免费视频| kizo精华| 久久青草综合色| 国产黄片美女视频| 久久午夜福利片| 欧美日韩一区二区视频在线观看视频在线| 夜夜骑夜夜射夜夜干| 亚洲av欧美aⅴ国产| 久久99热6这里只有精品| 最黄视频免费看| 免费黄频网站在线观看国产| 国产精品三级大全| 亚洲久久久国产精品| 亚洲国产最新在线播放| 91久久精品电影网| 亚洲成人一二三区av| 国产成人91sexporn| 99久久精品一区二区三区| 青春草亚洲视频在线观看| 人妻制服诱惑在线中文字幕| 热re99久久国产66热| 在线看a的网站| 青春草国产在线视频| 少妇丰满av| 在线观看人妻少妇| 国产精品三级大全| 精品亚洲成a人片在线观看| 观看免费一级毛片| 观看av在线不卡| 国产深夜福利视频在线观看| av又黄又爽大尺度在线免费看| 99久久精品一区二区三区| 成人18禁高潮啪啪吃奶动态图 | 在线观看www视频免费| 日本wwww免费看|