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

    Adaptive Backstepping Control Design for Semi-Active Suspension of Half-Vehicle With Magnetorheological Damper

    2021-04-16 03:56:42KhalidElMajdoubFouadGiriandFatimaZahraChaoui
    IEEE/CAA Journal of Automatica Sinica 2021年3期

    Khalid El Majdoub, Fouad Giri, and Fatima-Zahra Chaoui

    Abstract—This paper investigates the problem of controlling half-vehicle semi-active suspension system involving a magnetorheological (MR) damper. This features a hysteretic behavior that is presently captured through the nonlinear Bouc-Wen model.The control objective is to regulate well the heave and the pitch motions of the chassis despite the road irregularities. The difficulty of the control problem lies in the nonlinearity of the system model, the uncertainty of some of its parameters, and the inaccessibility to measurements of the hysteresis internal state variables. Using Lyapunov control design tools, we design two observers to get online estimates of the hysteresis internal states and a stabilizing adaptive state-feedback regulator. The whole adaptive controller is formally shown to meet the desired control objectives. This theoretical result is confirmed by several simulations demonstrating the supremacy of the latter compared to the skyhook control and passive suspension.

    I. INTRODUCTION

    AUTOMOTIVE control has become a driving factor in automotive innovation over the last twenty-five years. In order to meet the enhanced requirement for lower fuel consumption, lower exhaust emission, improved safety as well as comfort and convenience function, automotive control had to be applied. Faced with the need to reduce gas emissions,the hybrid vehicle or the electric vehicle is now a credible option [1]-[5], and already commercialized despite the intrinsic constraints associated with dual powering and extra equipment costs.

    Vehicle suspension system control aims at ensuring passenger’s safety and comfort by reducing vibrations caused by the onboard engine and the road irregularities [6]-[9].Nowadays, magnetorheological (MR) dampers have become among the most effective means toward the achievement of this objective, due to versatility, low-cost, small-size, high achievable forces. They are constituted of MR fluids that reversibly change their rheological properties in the presence of a varying magnetic field. Then, the fluid force can be acted on by changing the magnetic field amplitude, an action that necessitates low energy. These appealing features explain why MR dampers have found, in recent years, several applications in so various domains such as vibration isolation and damping,earthquake and civil engineering, automotive and vehicle industry [10]-[12].

    The complexity of MR dampers dynamics lies in the nonlinearity and the hysteretic nature of the force-velocity relationship [13]. Different, more or less complex, dynamic models have been proposed to capture this relationship[14]-[16]. Most of them are in fact variants of the Bouc-Wen hysteresis model and the LuGre friction model.

    For the MR dampers to be used as actuators in half-vehicle semi-active system suspension, they need to be appropriately controlled. To this end, several control approaches have been proposed, ranging from simple on-off control techniques to much more advanced linear and nonlinear control techniques including the optimal control and the linear quadratic control(LQR) technique [17]-[19], linear parameter varying (LPV)technique [20], adaptive backstepping active or semi-active suspension for a half-vehicle or quarter vehicle [21]-[24]. In[25], [26], a skyhook damper control algorithm has been proposed for full-vehicle suspension. A hybrid control technique, combining skyhook control and the Groundhook control, has been proposed in [27]. The disturbance (road profile) was assumed to be not a measurable and predictable signal. The control techniques were based on the on-off switched damping strategy, switching between high and low damping coefficients in order to achieve body comfort or road holding specifications.

    Another research direction is that using neural network control [28] and fuzzy logic control [29], which are two techniques that do not require the perfect knowledge of the model, but the expertise of the system is mandatory. Other existing control techniques include H2optimal control [30].Optimal and model predictive controllers were combined in[31] leading to semi-active hybrid controllers switching between different control laws.

    A fraction of the existing research activity has been devoted to adaptive control design for vehicle suspension. In [32] an adaptive fault-tolerant control design method was proposed for nonlinear pure-feedback systems (with dead-zone actuators and stochastic failures) and validated by simulation on a quarter-car suspension system. The method is not applicable to the problem studied in this work because the considered half-car suspension system does not fit the required pure-feedback structure, due to the hysteresis dynamics of the rheological actuators. The problem of adaptive (finite time) fault-tolerant control was also investigated in [33] for a class of switched nonlinear systems in the lower-triangular form (under arbitrary switching signals) where the unknown functions were approximated by the radial basis function neural networks (RBFNNs). The proposed method applies to single control input systems but it is not applicable to the present two-input two-output system.Furthermore, it is not judicious to use RBFNNs when the system model is available (even in the presence of parameter uncertainty), which is the case in the present paper.

    A common shortcoming of all mentioned control systems is that all system states have been supposed to be accessible to measurements and some only achieved a part of the desired performances. Some works involved switching between different control laws, making complex the formal analysis of the closed-loop performances [28]. Others used fully active suspension actuators [21] which are not yet available on most of today’s vehicles because of their high costs. As a matter of fact, semi-active actuators are generally preferred in vehicle industry applications.

    In the light of the above description, it becomes clear that most proposed (model-based) control methods for semi-active suspension have based on linear or LPV models. The point is that these models are approximate representations limiting the achievable performances.

    The present paper constitutes a significant progress with respect to the existing works because the control design is presently based on a model that accounts for the system nonlinearity, the model parameter uncertainty, and the hysteretic nature of MR damper dynamics. The dynamics of nonlinearity and hysteresis have been captured using the Bouc-Wen model. The Bouc-Wen model is more accurate, but more complex, than the Dahl model. Of course, more complex hysteresis models could be considered but the point is that not all models allow the design of controllers using the available nonlinear control design techniques. Therefore, the control design model should be a compromise between model accuracy and complexity. From this viewpoint, Bouc-Wen model is quite convenient and so is presently based upon.

    This paper proposes a new (model-based) adaptive backstepping semi-active suspension control for vehicles involving a semi-active component, i.e., MR damper. The aim is to manage the tradeoff between conflicting requirements, i.e.,ride comfort and driving safety. On the other hand, stabilizing the attitude of the vehicle in the presence of parameter uncertainties is a principal requirement, where suspension spaces and dynamic tire loads are considered to be guaranteed within their allowable ranges. The complexity of the considered control problem is threefold: i) the involved MR dampers dynamics are nonlinear and hysteretic; ii) the state variables associated to MR dampers hysteresis are not accessible to measurements; iii) some parameters of the overall half-vehicle system are uncertain. The control problem is dealt with using the adaptive backstepping design technique and Lyapunov analysis tools [34], [35]. The control design is based on the complete system model that accounts not only for the linear halfvehicle vertical motion but also for angular movement and the hysteretic behavior of the semi-active MR dampers. The adaptive controller includes three main parts: i) state observers providing online estimates of the damper internal state variables; ii) a parameter adaptive law providing online estimates of the uncertain parameters; iii) a nonlinear control law ensuring stability and regulation of the bouncing and pitching motions. The adaptive controller thus constructed is formally shown to meet its objectives.

    The paper is organized as follows: Section II is devoted to the modeling of the half-vehicle suspension system with MR damper; the controller design is presented in Section III; the controller performances are illustrated by numerical simulation in Section IV. To alleviate the paper body, some technical parts are placed in Appendices.

    II. NONLINEAR HALF-VEHICLE MODEL AND SEMI-ACTIVE COMPONENT MODELING

    A. Half-Vehicle Model

    Vehicle suspension systems are mainly composed of springs and dampers that filter effort transmission between the vehicle body and the road. Driving safety is the consequence of a harmonious suspension design in terms of wheels suspension,springing, steering, and braking. The driving comfort is closely related to the vibration level: the lower vehicle vibration level causes the better driving comfort quality.Considering the vehicle’s symmetry, the study can be dealt with based on the so-called half-vehicle model shown in Fig.1.This model is characterized by the sprung mass msand the unsprung masses for the front and rear axles mus1and mus2,respectively. As the damping coefficient is negligible, the front and rear tires are simply modeled by springs kus1and kus2linked to the road contact points zr1and zr2.

    Fig.1. Model of half vehicle with MR dampers.

    The passive suspension between chassis and the front, and rear axles is modeled by two dampers Cs1, Cs2(linearized viscous damping coefficients for the front and rear dampers)and two springs ks1, ks2(spring stiffness for the front and rear dampers). The non-linear half-vehicle model used in this study has four Degree of Freedom due to the heave and pitch motions of the sprung mass and the vertical motions of the unsprung masses. In this model, J stands for the mass moment of inertia for the pitch motion. For the vehicle body, zsand θ represent the heave and pitch motions, respectively. zs1and zs2are the displacements for the front and rear corners of the chassis vehicle, respectively. zus1and zus2are the unsprung masses displacements for the front and rear axles,respectively. In Fig.1, a and b show the distances of the suspensions to the center of the vehicle body and ? is the velocity of the vehicle.

    The numerical values of the vehicle parameters are given in Table I. Analytically, the nonlinear passive model based upon later, is defined by the following equation:

    Bouncing movement:

    Pitching movement:

    Vertical front and rear axles movements:

    TABLE I PARAMETER NUMERICAL VALUES OF THE HALF-VEHICLE MODEL

    B. MR Damper Modeling

    Presently, the Bouc-Wen model is retained because it represents a good accuracy/complexity ratio.

    Fig.2 illustrates how the hysteresis effect in MR dampers is captured using this model. Analytically, the Bouc-Wen model consists of the following couple of equations for the front and rear dampers, indexed by ( i=1,2), respectively:

    where the different notations are defined as follows: ωiis the internal state variable, viis the input voltage, ko1,ko2are the linear spring stiffness, αa1, αa2are the stiffness of ωi, αbiis the stiffness of ωiinfluenced by vi, Coaiis the viscous damping coefficient, Cobiis the viscous damping coefficient influenced by viand ρi, βi, -λiare the positive parameters characterizing the shape and size of hysteresis loop.

    Fig.2. Bouc-Wen model for MR damper.

    The Bouc-Wen model internal state variable ωiis not accessible to measurements but is formally proved to be bounded. In effect, it is shown in [36], if

    Then, the signal ωiremains bounded as follows:

    C. Half-Vehicle With MR Damper Model

    For convenience, the following vector state is introduced:

    where x1and x2are the vertical position and bouncing velocity of the centre of the gravity of the vehicle,respectively. x3and x4denote the pitch position and pitch velocity of the vehicle, respectively. On the other hand, x5, x6and x7, x8denote the vertical displacements and vertical speeds for the front and rear unsprung masses, respectively. In coherence with the notations (6), the state variables of the Bouc-Wen models are denoted: x9=ω1, x10=ω2. Let us introduce the novel notation for the MR damper forces,ui=Fmri. Then, the whole controlled system, including the motion equations (1a-d) and the damper equations (3) and (4),turns out to be described by the following state-space model:

    The parameters ki(i=1,10) and hi(i=1,8) are defined in Table II.

    TABLE II THE PARAMETERS k i(i=1,10) AND hi(i=1,8)

    In (7k-l), the control inputs are v1and v2. Observe that they are multiplied by states that are not accessible to measurements, namely x9for v1and x10for v2. It is assumed that all state variables x1to x8are accessible to measurements,except for the damper states x9and x10.

    In the system suspension with parametric uncertainty there are two types of nonlinearities. The first is smooth nonlinearity introduced by the trigonometric functions and the second is the Bouc-Wen hysteresis nonlinearity introduced by the MR damper actuator. Interestingly, the parameters kiare not required to be exactly known. This is coherent with the fact that the physical quantities ms,J,ks1,ks2,Cs1, and Cs2are likely to change with the operation conditions. It is only required that

    where the lower and upper bounds (k5min,k5max,k10min,k10max)are supposed to be a priori known. This assumption is a realistic one since k5and k10are related to the chassis mass and inertia respectively, whose variation ranges are a priori known for each given car. From Table I, one gets that k5min=1/msmaxand k10min=1/Jmax. In (7f-h), the parameters hi(i=1,8) are exactly known. All the parameters in hysteresis Bouc-Wen model are exactly known in (7i-l).

    III. ADAPTIVE CONTROL DESIGN

    A. Control Objectives

    Consider a half-vehicle equipped by two MR dampers (Fig.1)and described by (7) and (8). We seek a controller that meets the following objective.

    1) Bouncing and Pitching Regulation: The bouncing and pitching positions x1and x3must track their reference values,denoted x1refand x3ref, presently supposed to be constant. The road contact point zr1and zr2are considered to be known bounded external disturbance.

    2) Ride Comfort: The vertical acceleration of the car body must be limited and improved using the considered controller within an admissible range.

    3) Good Road Holding: The dynamic tires loads should not exceed the static ones for both of the front and rear wheels[21], i.e.,

    where the static tire loads Fi(i=1,2) are given by

    4) Suspension Space Limits: Because of the mechanical structure (front and rear axles), the suspension spaces should not exceed the allowable maxima, i.e.,

    B. Controller Design

    The controller design is performed in two stages. The first stage is relatively short and will consist in designing an observer structure providing estimates of the hysteresis internal variables x9and x10. Then, the second stage will be devoted to designing an adaptive regulator that meets the control objectives of Section III-A.

    1) Observer Structure of the Damper Internal States

    The internal states x9and x10are not accessible to measurements. To get online estimates of these variables, the following observer structure, suggested by (7i-l), is considered

    In the above observer structure, the observation of the state variables x9and x10relies on the measurements of signals u1and u2. Since these signals enter the sub-model (7i-l) through multiplication by the quantities ( αa1+αb1v1) and ( αa2+αb2v2),it follows that observability of the sub-model (7i-l) is only ensured unless the mentioned quantities remain away from zero. Therefore, we make the following observability assumptions:

    An observer analysis will be performed later as a part of the whole adaptive controller analysis. In particular, the way the observer gains should be online tuned will be determined bearing in mind the achievement of the output-reference tracking objective. To prepare that analysis, introduce the following observation errors:

    Then, it is readily seen, subtracting side-to-side (13a-d)from (7i-j), that the observation errors undergo the following equations:

    2) Adaptive Backstepping Controller Design

    We seek a bouncing-pitching position regulator that enforces the following position errors to asymptotically vanish, whatever the initial conditions:

    x1refand x3refdenote the bouncing and pitching position reference values supposed to be constant.

    To compensate for parameter uncertainty, the position regulator will be provided with a parameter adaptation capability. To this end, introduce the following notations:

    The control law, the parameter adaptation law and the observer gain will now be designed on the basis of (7a-d),using the backstepping technique [35]. This is performed in three steps. For convenience, the following new compact form of the system (7a-d) is introduced:

    Design Step 1 (Stabilization Function Design for Subsystem(18a) and (18c)):

    Consider the following Lyapunov function candidate for the subsystem (18a) and (18c):

    The corresponding time-derivative is

    In (18a) and (18c), the variables x2and x4stands for the virtual control signals. Let α1and α3denote the corresponding stabilization functions. Equation (20) suggests that

    with this notation, (20) and (22a-b) can be rewritten as

    Design Step 2 (Control Law Design for the System (18b)and (18d)):

    From (22a-b) it follows, using (23a-b), (18b), and (18d),that the errors e2and e4undergoe the following equations:

    Consider the Lyapunov function candidate

    Fig.3. Vehicle suspension adaptive backstepping control.

    TABLE IV THE NUMERICAL VALUES OF THE MR DAMPERS PARAMETERS

    As a semi-active control approach, the skyhook control strategy supposes that the sprung body is connected to a virtual sky by a damper [25]. The front and rear dampers are indexed by (i=1,2), respectively. The skyhook control law is mathematically described by

    The passive suspension is a simpler variant of the damping force described by (3).

    where the index “p” refers to “passive”. The parameters Coaiand koi( i=1,2) are given in Table IV.

    The root mean square (RMS) value of the vehicle body acceleration is often used as a measure of passenger drive comfort. That is, the vertical acceleration RMS is considered in this study to evaluate the quality of ride comfort. Recall that the RMS value of an n-dimensional vector x is defined as follows:

    where ‖·‖ denotes the usual Euclidian norm. The acceleration peak value is also of some interest. It refers to the maximal magnitude of vehicle body or passenger acceleration. The peak value is

    where ‖ ·‖∞denotes the infinity norm.

    The static tire loads for both front and rear wheels can be calculated by (11a-b) as F1=5718N and F2=4768N.

    The control design parameters are given the values c1=1,c2=0.01, c3=1, c4=0.01, δ1=0.01, δ2=0.01, γ=1,λo1=10, and λo2=10, which proved to be suitable.

    A. Time Domain Analysis With Deterministic Road

    The vehicle speed produces an effect on the performances of the half vehicle suspension model. Indeed, the speed determines the frequency of the disturbances caused by road irregularities which in turn affect the suspension frequency behavior. Two considered road profiles are referred to bump road input and road input with limited ramp. The input zr2for the rear wheels is the same as that for the front ones but delayed by a certain time τ =(a+b)/?.

    Bump Road Input: This profile is the most encountered in practice; it is analytically described by

    where hb=4 cm represents the height of the bump road input,T1=1s and T2=2s. This profile is shown by Fig.4 and the system responses, obtained with the proposed semi-active control and with the passive suspension, are plotted in Figs. 5(a)(chassis heave displacement) and 5(b) (heavy acceleration),Figs. 6(a) (pitch motion) and 6(b) (pitch acceleration).

    Fig.4. The bump road inputs.

    The variation of the front damper force is shown in Fig.7(a),and the rear damper force is plotted in Fig.7(b). It is seen from Figs. 5(a), 5(b) and 6(a), 6(b) that the amplitude for the vertical position, the pitch motion, the bouncing acceleration and the pitch acceleration vanish much faster with the semiactive suspension control than with skyhook control and passive suspension. To quantify the ride comfort improvement, the RMS values for the vertical acceleration of the vehicle body are computed using (41). The RMS values for the bouncing and pitch accelerations of the vehicle body are presented in Table V. It is seen that these values are greatly reduced with the adaptive controller resulting in improved ride comfort.

    Fig.5. (a) The heave of the sprung mass for the bump road; (b) The heave acceleration of the sprung mass for the bump road.

    Fig.6. (a) The pitch of the sprung mass for the bump road; (b) The pitch acceleration of the sprung mass for the bump road.

    Fig.7. (a) The front damper force for the bump road; (b) The rear damper force for the bump road.

    TABLE V COMPARATIVE RESULT BETWEEN DIFFERENT CONTROL LAW

    The decrease in the RMS values is 77% for the bouncing acceleration and 95% for the pitch acceleration.

    The comparison between suspensions techniques can also be checked by analyzing the peak values (42). As seen, the peak values are substantially decreased by the adaptive controller, that is, the maximum accelerations sensed by the passengers are decreased by 54% for the heave acceleration and 90% for the pitch acceleration. This results in increased ride comfort. In the adaptive suspension control, the limitations of the suspension spaces should be taken into account, it can be observed from Figs. 8(a) (front working space) and 8(b) (rear working space) that the controlled suspension spaces are all below the limitations (12). The road holding is computing by (10). Figs. 9(a) (dynamic front tire load) and 9(b) (dynamic rear tire load) show that the dynamic tires loads are all within the bounds.

    Fig.8. (a) The front working space for the bump road; (b) The rear working space for the bump road.

    Limited Ramp: A road with limited ramp, depicted by Fig.10,is considered. Analytically, the road shape is defined as follows, where hr=0.02 m is the final road surface elevation,T1=1s and T2=1.1s:

    Fig.9. (a) The dynamic front tire load for the bump road; (b) The dynamic rear tire load for the bump road.

    Fig.10. The limited ramp road inputs.

    The obtained control performances are illustrated by Figs. 11(a), 11(b), 12(a), and 12(b). These show the halfvehicle body vertical displacement (Fig.11(a)), the vertical acceleration (Fig.11(b)), the pitch motion (Fig.12(a)), and the pitch acceleration (Fig.12(b)). It is seen that the vehicle settles smoothly to the final height value of the road input and the vertical displacement magnitude and acceleration are not large. Let us evaluate the RMS value of the vertical motion acceleration and the pitch acceleration of the sprung mass(Table V). The ride comfort is significantly improved with the adaptive controller; the RMS decreases by 85% for the bouncing acceleration and decreases by 90% in the case of the pitch acceleration. The peak value of the vertical motion acceleration is decreases by 60%. On the other hands, the peak value is decreases by 88% in the case of the pitch acceleration.That is, the decrease is quite significant. The front and rear MR dampers forces produced by adaptive controller, with respect to time are depicted in Figs. 13(a) and 13(b). It is seen from Figs. 14(a) (front working space) and 14(b) (rear working space) that the controlled suspension spaces are respected. The road holding is checked by Figs. 15(a)(dynamic front tire load) and 15(b) (dynamic rear tire load).That shows that the dynamic tires loads are all within the bounds.

    Fig.12. (a) The pitch of the sprung mass for the limited ramp road; (b) The pitch acceleration of the sprung mass for the limited ramp road.

    Fig.13. (a) The front damper force for the limited ramp road; (b) The rear damper force for the limited ramp road.

    Fig.14. (a) The front working space for the limited ramp road; (b) The rear working space for the limited ramp road.

    Fig.15. (a) The dynamic front tire load for the limited ramp road; (b) The dynamic rear tire load for the limited ramp road.

    Fig.16. (a) The heave of the sprung mass for the random road; (b) The heave acceleration of the sprung mass for the random road.

    B. Time Domain Analysis With Random Road With C Level, ISO 8608

    Now, a random road profile is considered to illustrate the validity and effectiveness of the presently proposed controller.In this test, a zero reference is considered, i.e., zsref=0.Random road profile is generated by linear filtering methods and harmonic superposition method [38]. The main idea is to generate random number road profile with road class roughness-C according to the inter-national standard organization (ISO) 8608. The equation of the random road profile is as follows:

    where f0=0.011 Hz stands for the offline cutoff frequency, zriis the road roughness amplitude, ? is a white noise whose intensity is 1, ? is the vehicle speed, n0is the frequency of reference space whose value is 0.1 Hz, G=256×106m3is the road roughness coefficient for road of C level.

    The responses of a half vehicle model using adaptive controller, on one hand, and skyhook and passive control to random road excitation, on the other hand, are presented in the subsequent lines. The sprung mass vertical displacement, the bouncing acceleration, the pitch motion and the pitch acceleration are plotted in Figs. 16(a), 16(b), 17(a), and 17(b),respectively. It is seen that the magnitudes of the heave, pitch and corresponding accelerations are significantly decreased with the adaptive controller, confirming the efficiency of the latter in the presence of known road. The decreases in the RMS values of the heave and pitch accelerations (of the vehicle body) are 80% and 95%, respectively. Such large decrease of RMS values entails an improved ride comfort. In turn, the Peak values for the heave and pitch accelerations are substantially decreased with the adaptive controller.Specifically, the maximum accelerations sensed by the passengers are decreased by 90% and 97%, respectively,which results in substantially increased ride comfort.

    V. CONCLUSION

    This paper presents a new approach of controlling the vertical movement of a half car vehicle facing deterministic with a constant velocity. The main component of the proposed suspension system is an adaptive controller designed by the adaptive backstepping technique on the basis of a model that accounts for the hysteresis effect in the MR dampers. The controller consists of two observers estimating the damper hysteresis internal states, a parameter adaptive algorithm estimating the system uncertain parameters and an adaptive state-feedback control law stabilizing the suspension system.The adaptive controller is formally shown to meet its control objectives. This theoretical result is confirmed by several simulations that also illustrate the high supremacy of the semiactive suspension over the passive suspension and the semiactive skyhook suspension. Compared to previous works on semi-active suspension, our solution presents several features:i) this solution involves a new component such as MR dampers modeled by the Bouc-Wen model; ii) it does not assume all system states and all system parameters to be known; iii) it does not involve switching between different control laws; iv) it enjoys a formal closed-loop stability analysis. Finally, it is checked with the simulation study of Section IV that the backstepping adaptive controller involving MR dampers improves a good ride comfort and road holding.

    Fig.17. (a) The pitch of the sprung mass for the random road; (b) The pitch acceleration of the sprung mass for the random road.

    This study can be continued in many directions, e.g.,designing a higher supervision loop computing online optimal values of the control design parameters; optimality is to be defined in terms of control performances.

    APPENDIX A MODEL AND CONTROLLER PARAMETERS AND FUNCTIONS

    国产av一区二区精品久久| 热99re8久久精品国产| www.精华液| 国产精品亚洲美女久久久| 欧美亚洲日本最大视频资源| 国产欧美日韩一区二区三| 欧美激情高清一区二区三区| tocl精华| 国产在线精品亚洲第一网站| 精品国产超薄肉色丝袜足j| av在线播放免费不卡| 99国产综合亚洲精品| 国产国语露脸激情在线看| 久久久久久久久免费视频了| 精品国产超薄肉色丝袜足j| 国产欧美日韩一区二区三| 电影成人av| 91麻豆av在线| 亚洲av五月六月丁香网| 国产真人三级小视频在线观看| 黄色视频,在线免费观看| 波多野结衣一区麻豆| 最近最新中文字幕大全电影3 | 窝窝影院91人妻| 日本黄色视频三级网站网址| 日本精品一区二区三区蜜桃| 男人操女人黄网站| 日日爽夜夜爽网站| 久久久久亚洲av毛片大全| 久久久久久国产a免费观看| 亚洲人成电影观看| 久久青草综合色| 国产一级毛片七仙女欲春2 | 午夜精品久久久久久毛片777| 国产精品一区二区在线不卡| 精品不卡国产一区二区三区| 亚洲av美国av| 男男h啪啪无遮挡| 亚洲成人久久性| 天堂动漫精品| 女警被强在线播放| 美国免费a级毛片| 人人妻,人人澡人人爽秒播| 琪琪午夜伦伦电影理论片6080| 欧美日韩福利视频一区二区| 丝袜美足系列| 大型黄色视频在线免费观看| 露出奶头的视频| 黄色片一级片一级黄色片| 可以在线观看毛片的网站| 不卡一级毛片| 久久久久久大精品| 成年女人毛片免费观看观看9| 色播在线永久视频| 成人三级黄色视频| 国产精品久久视频播放| 高潮久久久久久久久久久不卡| 亚洲国产高清在线一区二区三 | 色综合站精品国产| 女人被狂操c到高潮| 日韩国内少妇激情av| 在线视频色国产色| 日本 欧美在线| √禁漫天堂资源中文www| 一边摸一边抽搐一进一出视频| 97人妻精品一区二区三区麻豆 | 欧美乱妇无乱码| 精品久久久久久久久久免费视频| 国产高清视频在线播放一区| 国产成人精品在线电影| 精品午夜福利视频在线观看一区| 国产精品99久久99久久久不卡| 精品人妻在线不人妻| 亚洲国产欧美一区二区综合| 亚洲情色 制服丝袜| 在线观看www视频免费| 亚洲av美国av| 亚洲欧美激情综合另类| 女性被躁到高潮视频| 国产一区在线观看成人免费| 在线观看午夜福利视频| 一级毛片精品| 亚洲 国产 在线| 黑丝袜美女国产一区| 免费不卡黄色视频| 亚洲va日本ⅴa欧美va伊人久久| 亚洲第一av免费看| 十八禁人妻一区二区| 亚洲情色 制服丝袜| 亚洲国产欧美日韩在线播放| 黑丝袜美女国产一区| 国产精品亚洲av一区麻豆| 久久精品亚洲熟妇少妇任你| 午夜免费观看网址| 国产麻豆成人av免费视频| 色播在线永久视频| 一区二区日韩欧美中文字幕| 香蕉国产在线看| 成人欧美大片| 丁香欧美五月| 国产精华一区二区三区| 日日夜夜操网爽| 日韩精品青青久久久久久| 日韩视频一区二区在线观看| 一卡2卡三卡四卡精品乱码亚洲| 国产乱人伦免费视频| 久久天堂一区二区三区四区| 女警被强在线播放| 一本大道久久a久久精品| 巨乳人妻的诱惑在线观看| 免费搜索国产男女视频| 啦啦啦韩国在线观看视频| 嫩草影视91久久| 欧美日韩亚洲综合一区二区三区_| 久久久国产成人精品二区| 国产亚洲精品一区二区www| 亚洲第一av免费看| 欧美色视频一区免费| 久久国产亚洲av麻豆专区| av天堂在线播放| 欧美成人性av电影在线观看| 中文字幕最新亚洲高清| av中文乱码字幕在线| 狠狠狠狠99中文字幕| 成人18禁在线播放| 国产精品久久久久久亚洲av鲁大| 国产精品久久久久久人妻精品电影| 日本在线视频免费播放| 日韩欧美国产一区二区入口| 一级片免费观看大全| 91大片在线观看| 纯流量卡能插随身wifi吗| 黄色视频,在线免费观看| 亚洲av电影不卡..在线观看| 妹子高潮喷水视频| 搡老妇女老女人老熟妇| 欧美一级a爱片免费观看看 | 国产男靠女视频免费网站| 国产伦人伦偷精品视频| 久久人人爽av亚洲精品天堂| 国产男靠女视频免费网站| 人成视频在线观看免费观看| 日韩大尺度精品在线看网址 | 亚洲国产精品久久男人天堂| 男人舔女人的私密视频| 亚洲av熟女| 午夜久久久久精精品| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品中文字幕一二三四区| 级片在线观看| 国产av在哪里看| 亚洲国产精品成人综合色| 亚洲情色 制服丝袜| 神马国产精品三级电影在线观看 | 国产精品爽爽va在线观看网站 | 啦啦啦韩国在线观看视频| 国产三级在线视频| 国产色视频综合| 99国产精品99久久久久| 丰满的人妻完整版| av福利片在线| 一级a爱视频在线免费观看| 亚洲久久久国产精品| x7x7x7水蜜桃| 日本黄色视频三级网站网址| 日本免费a在线| 高清在线国产一区| 午夜影院日韩av| 在线观看午夜福利视频| 国产成+人综合+亚洲专区| 男女做爰动态图高潮gif福利片 | 国产一级毛片七仙女欲春2 | 成人18禁在线播放| 国产三级在线视频| 欧美人与性动交α欧美精品济南到| 咕卡用的链子| 亚洲熟妇中文字幕五十中出| 欧美一级毛片孕妇| 欧美亚洲日本最大视频资源| 亚洲黑人精品在线| 搡老熟女国产l中国老女人| 在线观看免费视频网站a站| 欧美色视频一区免费| 免费在线观看视频国产中文字幕亚洲| 亚洲视频免费观看视频| 欧美黄色淫秽网站| 一区二区三区高清视频在线| 久久天躁狠狠躁夜夜2o2o| 精品久久久久久久人妻蜜臀av | 国产精品亚洲av一区麻豆| 一边摸一边做爽爽视频免费| 香蕉丝袜av| 天堂影院成人在线观看| 欧美av亚洲av综合av国产av| 国产精华一区二区三区| 欧美不卡视频在线免费观看 | 宅男免费午夜| 欧美一级毛片孕妇| 神马国产精品三级电影在线观看 | 国产成人精品无人区| 啦啦啦观看免费观看视频高清 | 九色国产91popny在线| 久久久久亚洲av毛片大全| 久久精品91无色码中文字幕| 午夜福利影视在线免费观看| 日韩三级视频一区二区三区| 国产一级毛片七仙女欲春2 | 久久久久九九精品影院| 日韩国内少妇激情av| 欧美黄色片欧美黄色片| 免费久久久久久久精品成人欧美视频| 国产三级在线视频| 女警被强在线播放| 不卡一级毛片| 日本 欧美在线| 久久久久久久久中文| 色综合亚洲欧美另类图片| 人人妻人人澡人人看| 97人妻天天添夜夜摸| 91精品三级在线观看| 可以免费在线观看a视频的电影网站| 国产午夜福利久久久久久| 男女之事视频高清在线观看| 色播在线永久视频| 国产黄a三级三级三级人| 日韩欧美在线二视频| 久久人人97超碰香蕉20202| 久久天躁狠狠躁夜夜2o2o| 咕卡用的链子| 女性生殖器流出的白浆| 日韩大尺度精品在线看网址 | 啦啦啦韩国在线观看视频| 精品一品国产午夜福利视频| 亚洲成国产人片在线观看| 老司机靠b影院| av网站免费在线观看视频| 91精品国产国语对白视频| av在线播放免费不卡| 国产成人一区二区三区免费视频网站| 丝袜在线中文字幕| 国产精华一区二区三区| 麻豆成人av在线观看| 欧美激情高清一区二区三区| 啦啦啦韩国在线观看视频| 久久精品aⅴ一区二区三区四区| 国产成人av教育| 又紧又爽又黄一区二区| 国产成人精品久久二区二区91| 亚洲 欧美一区二区三区| 一区二区日韩欧美中文字幕| 欧美成人免费av一区二区三区| 天天躁夜夜躁狠狠躁躁| 麻豆国产av国片精品| 老司机在亚洲福利影院| 中文字幕另类日韩欧美亚洲嫩草| 在线国产一区二区在线| 国产一区二区激情短视频| 制服人妻中文乱码| 国产欧美日韩综合在线一区二区| 在线观看66精品国产| 91av网站免费观看| 在线观看舔阴道视频| 99国产精品免费福利视频| 多毛熟女@视频| 女人高潮潮喷娇喘18禁视频| 美女 人体艺术 gogo| 激情视频va一区二区三区| 两个人视频免费观看高清| 中文字幕高清在线视频| 亚洲精品美女久久av网站| 亚洲国产精品久久男人天堂| 国产不卡一卡二| 免费不卡黄色视频| 人人妻,人人澡人人爽秒播| 真人一进一出gif抽搐免费| 国产精品免费视频内射| 搡老妇女老女人老熟妇| 欧美乱色亚洲激情| 午夜福利免费观看在线| 91精品国产国语对白视频| 日韩欧美三级三区| 午夜久久久久精精品| 国产xxxxx性猛交| 国产精品秋霞免费鲁丝片| 亚洲精品国产区一区二| 国产精品久久视频播放| 国产精品1区2区在线观看.| 日本免费a在线| 中文字幕高清在线视频| 琪琪午夜伦伦电影理论片6080| 91老司机精品| 人人妻人人爽人人添夜夜欢视频| 国产精品永久免费网站| 日本 av在线| 可以免费在线观看a视频的电影网站| 韩国精品一区二区三区| 国产精品美女特级片免费视频播放器 | 久久人人97超碰香蕉20202| 一区二区三区精品91| 国产精品99久久99久久久不卡| 少妇裸体淫交视频免费看高清 | 两人在一起打扑克的视频| 少妇粗大呻吟视频| 欧美日韩一级在线毛片| 每晚都被弄得嗷嗷叫到高潮| 亚洲一卡2卡3卡4卡5卡精品中文| 999久久久精品免费观看国产| 一级片免费观看大全| 少妇熟女aⅴ在线视频| 亚洲人成77777在线视频| 亚洲精品av麻豆狂野| 国产成年人精品一区二区| 亚洲精品国产色婷婷电影| 日韩 欧美 亚洲 中文字幕| e午夜精品久久久久久久| 国产精品综合久久久久久久免费 | 国产精品1区2区在线观看.| 亚洲欧美一区二区三区黑人| 欧美性长视频在线观看| 亚洲av第一区精品v没综合| 精品久久久精品久久久| 91av网站免费观看| 不卡av一区二区三区| 国产精品精品国产色婷婷| 黄色毛片三级朝国网站| www.999成人在线观看| 久久热在线av| 欧美黑人欧美精品刺激| 国产精品一区二区免费欧美| www.999成人在线观看| 久久人妻av系列| 亚洲片人在线观看| 国内久久婷婷六月综合欲色啪| 日韩 欧美 亚洲 中文字幕| 久久草成人影院| 无遮挡黄片免费观看| 免费在线观看亚洲国产| 女人高潮潮喷娇喘18禁视频| 国产精品亚洲av一区麻豆| 国产男靠女视频免费网站| 麻豆国产av国片精品| 国产一级毛片七仙女欲春2 | 欧美一级a爱片免费观看看 | 国产精品一区二区三区四区久久 | 波多野结衣巨乳人妻| 国产一区二区三区综合在线观看| 啦啦啦韩国在线观看视频| 天堂动漫精品| 中文字幕av电影在线播放| www.999成人在线观看| 久久人人爽av亚洲精品天堂| АⅤ资源中文在线天堂| 波多野结衣巨乳人妻| 久久久久久久午夜电影| 中文字幕av电影在线播放| 国产黄a三级三级三级人| 久久青草综合色| 亚洲熟妇熟女久久| 真人做人爱边吃奶动态| 91成年电影在线观看| 亚洲中文日韩欧美视频| 天天添夜夜摸| 欧美性长视频在线观看| 久久久久久久久久久久大奶| 一边摸一边做爽爽视频免费| 如日韩欧美国产精品一区二区三区| 亚洲中文字幕日韩| 久久午夜亚洲精品久久| 久久久久久人人人人人| 老司机午夜十八禁免费视频| 人人妻人人澡欧美一区二区 | 91麻豆精品激情在线观看国产| 国内久久婷婷六月综合欲色啪| 国产免费av片在线观看野外av| 多毛熟女@视频| 婷婷六月久久综合丁香| 成人18禁高潮啪啪吃奶动态图| 免费av毛片视频| 黄色丝袜av网址大全| 亚洲激情在线av| 91大片在线观看| 午夜福利欧美成人| 非洲黑人性xxxx精品又粗又长| 欧美日韩黄片免| 国产色视频综合| 色老头精品视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品中文字幕在线视频| 可以在线观看的亚洲视频| 亚洲专区中文字幕在线| aaaaa片日本免费| av中文乱码字幕在线| 亚洲成国产人片在线观看| 97人妻精品一区二区三区麻豆 | 国产激情欧美一区二区| 久久久国产精品麻豆| 亚洲国产日韩欧美精品在线观看 | 欧美日韩亚洲国产一区二区在线观看| 在线永久观看黄色视频| 人人澡人人妻人| 男女床上黄色一级片免费看| 窝窝影院91人妻| 91成年电影在线观看| 一本综合久久免费| 色播亚洲综合网| 侵犯人妻中文字幕一二三四区| 黄片大片在线免费观看| 少妇粗大呻吟视频| 夜夜夜夜夜久久久久| 啦啦啦免费观看视频1| 欧美中文综合在线视频| 日韩欧美三级三区| 看片在线看免费视频| 90打野战视频偷拍视频| 免费在线观看视频国产中文字幕亚洲| av超薄肉色丝袜交足视频| 亚洲欧美日韩另类电影网站| 久久午夜亚洲精品久久| 亚洲伊人色综图| 女人被狂操c到高潮| av在线播放免费不卡| 级片在线观看| 亚洲精品国产精品久久久不卡| 亚洲国产欧美网| 中文字幕另类日韩欧美亚洲嫩草| 国产精品99久久99久久久不卡| 黄片大片在线免费观看| 18禁观看日本| 女人被狂操c到高潮| av欧美777| 最近最新中文字幕大全电影3 | 大型av网站在线播放| 两性午夜刺激爽爽歪歪视频在线观看 | 国产三级在线视频| 法律面前人人平等表现在哪些方面| 亚洲视频免费观看视频| 亚洲精品av麻豆狂野| 欧美日韩一级在线毛片| 亚洲全国av大片| 97超级碰碰碰精品色视频在线观看| 美女大奶头视频| 欧洲精品卡2卡3卡4卡5卡区| 十八禁网站免费在线| 在线国产一区二区在线| 亚洲av第一区精品v没综合| 一级a爱视频在线免费观看| 香蕉国产在线看| 日韩精品青青久久久久久| 中文字幕色久视频| 国产91精品成人一区二区三区| 欧美一级a爱片免费观看看 | av有码第一页| 亚洲最大成人中文| 国产亚洲av嫩草精品影院| 男人舔女人的私密视频| 婷婷丁香在线五月| 18禁国产床啪视频网站| 国产精品野战在线观看| 免费女性裸体啪啪无遮挡网站| 亚洲av五月六月丁香网| 波多野结衣一区麻豆| 国产男靠女视频免费网站| 亚洲精品国产一区二区精华液| 精品乱码久久久久久99久播| 精品卡一卡二卡四卡免费| 国产成人精品久久二区二区91| 久久香蕉精品热| 亚洲一卡2卡3卡4卡5卡精品中文| 精品久久久久久,| 色尼玛亚洲综合影院| 天堂√8在线中文| 香蕉久久夜色| 一边摸一边抽搐一进一小说| 久久人人97超碰香蕉20202| 亚洲第一欧美日韩一区二区三区| 久久天堂一区二区三区四区| 精品高清国产在线一区| 日日干狠狠操夜夜爽| 亚洲国产高清在线一区二区三 | 亚洲国产精品久久男人天堂| 99在线人妻在线中文字幕| 国产精品亚洲美女久久久| 午夜福利在线观看吧| 一区二区三区高清视频在线| 精品欧美一区二区三区在线| 国产精品,欧美在线| 免费一级毛片在线播放高清视频 | 夜夜看夜夜爽夜夜摸| 精品国产一区二区久久| av在线播放免费不卡| 国内精品久久久久精免费| 黑丝袜美女国产一区| 九色亚洲精品在线播放| av福利片在线| 一a级毛片在线观看| 50天的宝宝边吃奶边哭怎么回事| 欧美激情高清一区二区三区| 亚洲中文字幕日韩| 夜夜夜夜夜久久久久| 少妇被粗大的猛进出69影院| 亚洲中文日韩欧美视频| 日本 欧美在线| 午夜精品国产一区二区电影| 大码成人一级视频| 亚洲熟妇中文字幕五十中出| 不卡av一区二区三区| 给我免费播放毛片高清在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲欧洲精品一区二区精品久久久| 色综合站精品国产| 欧美精品亚洲一区二区| 国产国语露脸激情在线看| 一a级毛片在线观看| 免费观看精品视频网站| 亚洲av熟女| 亚洲成人久久性| 婷婷丁香在线五月| 亚洲熟妇熟女久久| 欧美不卡视频在线免费观看 | 国产亚洲av嫩草精品影院| 最好的美女福利视频网| 正在播放国产对白刺激| 亚洲性夜色夜夜综合| 国产精品亚洲一级av第二区| 日韩一卡2卡3卡4卡2021年| 国产精品久久久人人做人人爽| 99国产精品99久久久久| 色精品久久人妻99蜜桃| 国产精品影院久久| 欧美大码av| 国内毛片毛片毛片毛片毛片| 精品国产亚洲在线| 亚洲第一欧美日韩一区二区三区| 国产欧美日韩精品亚洲av| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜福利,免费看| 欧美乱码精品一区二区三区| 一本久久中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久久久久免费视频了| 叶爱在线成人免费视频播放| 亚洲欧美激情在线| 国产单亲对白刺激| 亚洲精品美女久久av网站| 又黄又爽又免费观看的视频| 欧美日韩一级在线毛片| 婷婷六月久久综合丁香| 精品卡一卡二卡四卡免费| 18禁裸乳无遮挡免费网站照片 | 18禁国产床啪视频网站| 国产亚洲精品一区二区www| 午夜福利视频1000在线观看 | 国产在线精品亚洲第一网站| 一夜夜www| 成人av一区二区三区在线看| 欧美日韩亚洲国产一区二区在线观看| 女生性感内裤真人,穿戴方法视频| 老司机在亚洲福利影院| 亚洲人成伊人成综合网2020| 亚洲激情在线av| 国产精品亚洲美女久久久| 精品欧美国产一区二区三| 亚洲第一av免费看| 国产精品野战在线观看| 一边摸一边抽搐一进一小说| 亚洲自拍偷在线| 午夜免费鲁丝| 啪啪无遮挡十八禁网站| 午夜免费成人在线视频| 精品久久久久久久久久免费视频| 亚洲av熟女| 欧美+亚洲+日韩+国产| 亚洲天堂国产精品一区在线| 午夜亚洲福利在线播放| 老熟妇仑乱视频hdxx| √禁漫天堂资源中文www| 国产精品自产拍在线观看55亚洲| 亚洲欧美激情在线| 99久久综合精品五月天人人| 琪琪午夜伦伦电影理论片6080| 女同久久另类99精品国产91| 女人被狂操c到高潮| 99在线人妻在线中文字幕| 九色亚洲精品在线播放| 最新在线观看一区二区三区| 国产精品二区激情视频| 岛国在线观看网站| 黄色片一级片一级黄色片| 国产精品,欧美在线| 女性被躁到高潮视频| 美女 人体艺术 gogo| 午夜福利,免费看| 亚洲一区中文字幕在线| 美女 人体艺术 gogo| 88av欧美| 777久久人妻少妇嫩草av网站| 亚洲欧美激情综合另类| 十八禁网站免费在线| 黄色成人免费大全| 女人精品久久久久毛片| 日韩三级视频一区二区三区| 精品少妇一区二区三区视频日本电影| 可以免费在线观看a视频的电影网站| 国产精品九九99| 黑人巨大精品欧美一区二区mp4| 极品人妻少妇av视频| 亚洲熟妇中文字幕五十中出| 国产精品98久久久久久宅男小说| 十分钟在线观看高清视频www| 国产一区二区在线av高清观看| 又紧又爽又黄一区二区| 99久久综合精品五月天人人| 国产成人精品在线电影| 亚洲一区二区三区色噜噜| 国产精品永久免费网站|