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

    A Model-Based Unmatched Disturbance Rejection Control Approach for Speed Regulation of a Converter-Driven DC Motor Using Output-Feedback

    2022-01-25 12:51:08LuZhangJunYangandShihuaLi
    IEEE/CAA Journal of Automatica Sinica 2022年2期

    Lu Zhang,Jun Yang,,and Shihua Li,

    Abstract—The speed regulation problem with only speed measurement is investigated in this paper for a permanent magnet direct current (DC) motor driven by a buck converter.By lumping all unknown matched/unmatched disturbances and uncertainties together,the traditional active disturbance rejection control (ADRC) approach provides an intuitive solution for the problem under consideration.However,for such a higher-order disturbed system,the increase of poles for the extended state observer (ESO) therein will lead to drastically growth of observer gains,which causes severe noise amplification.This paper aims to propose a new model-based disturbance rejection controller for the converter-driven DC motor system using output-feedback.Instead of estimating lumped disturbances directly,a new observer is constructed to estimate the desired steady state of control signal as well as errors between the real states and their desired steady-state responses.Thereafter,a controller with only speed measurement is proposed by utilizing the estimates.The performance of the proposed method is tested through experiments on dSPACE.It is further shown via numerical calculations and experimental results that the poles of the observer within the proposed control approach can be largely increased without significantly increasing magnitude of the observer gains.

    I.INTRODUCTION

    DIRECT current (DC) motors are applied extensively in industry,particularly in the systems requiring higherprecision,such as the automotive industry [1],elevator systems [2],robotic manipulators [3],etc.In the early days,the configuration of DC motor drives consists of a resistance network connected in parallel/series with DC motor,but due to the non-smooth regulation,it causes abrupt jerkiness and changes.Recently,DC-DC converters are used to act as smooth starters for DC motors [4],which provides a solution for the aforementioned problem.However,the converterdriven DC motor system usually faces various disturbances/uncertainties,which are briefly summarized as follows:1) unknown constant or time-varying load torques[5];2) unknown parametric uncertainties such as the external source voltage [6],the load resistor,the winding resistor [7],etc.;and 3) unmodeled dynamics such as the equivalent series resistance in converter systems and measurement error [8].Due to the aforementioned challenges,researchers propose various control strategies to solve the angular velocity regulation problem therein,such as proportional-integral (PI)control [9],backstepping control [10],flatness-based tracking control [11],sliding mode control [12],robust control [13],ADRC [14],etc.

    In light of appealing properties including conceptual simplicity and powerful disturbance rejection performance,the active disturbance rejection control (ADRC) method[14]–[16] has gained a lot of attention in the area of speed regulation of buck converter-driven DC motors [14],[17].An observer-based control strategy is proposed in [18] for a permanent magnet motor to reject harmonic torque disturbances,which combines observation with the linear PI controller.In [19],a triple-step controller is proposed to solve the speed regulation problem for a permanent-magnet DC motor subject to uncertainties and disturbances.Furthermore,the authors of [20] propose an observer-based nonlinear triplestep controller to improve the low-speed tracking performance for permanent-magnet DC motor systems.Both the two controller proposed by [19] and [20] utilize a reduced-order mathematical model,which neglects the electrical time constant of the DC motor.In [21],a flatness tracking controller is proposed by using the simple second-order model,but this method would decrease the system bandwidth.In [13] and [22],based on a full fourth-order dynamic model,the ADRC method is applied to address the speed regulation problem for the converter-driven DC motor system.Since the overall system is fourth-order,in the framework of traditional ADRC,it is required to perform fourth-order differentiations to obtain the so-called standard control form,the design procedure of which is clearly shown in [13] and [22].

    For the system under consideration,performing fourth-order differentiations and lumping system dynamics as a part of the unknown input term will lose system dynamic models.It is generally required to assign a conservatively higher observer bandwidth in traditional ADRC to achieve higher estimation and control precision as the lumped uncertainties usually have an unknown complicated structure.However,for such a fourth-order system,the increase of observer bandwidth means large poles,which will lead to a drastic rise of the amplitude of observer gains,causing significant amplification of measurement noises.On the other hand,the derivatives of the desired angular velocity may be unavailable in certain applications such as position control.The traditional ADRC solves this problem by adding tracking differentiators to estimate the unknown information,which considerably increases the complexity in implementation.In [23],a generalized extended state observer (ESO)-based control strategy is proposed for unmatched disturbance compensation in the framework of ADRC,where a static disturbance compensation gain is designed for compensation.It is claimed that the method in [23] can tackle the output-feedback control problem.However,it only provides partial output-feedback control since it usually needs more measurements to satisfy the observability conditions therein.In addition,the number of the extended states is equal to the number of disturbances to be estimated,which leads to a higher dimension of the observer [23].As such,it is imperative to investigate new control strategies for the DC-DC buck converter-driven DC motor system to solve these problems as well as ensure satisfactory control performance and disturbance rejection capability.

    Motivated by the traditional ADRC and the output regulation theory [24],[25],this paper proposes a new disturbance rejection control approach for the buck converterdriven DC motor system to improve the control performance in the presence of various parametric uncertainties and unknown load torque disturbances.The benefits as well as differences from existing works of the proposed method are summarized as follows:

    1) Delicate Analysis of the Benefits of the Proposed Controller:The explicit relationship between the observer gains and the assigned poles has been developed for both the proposed method and the ADRC method.The results show that the observer gains under the proposed observer are much lower than the traditional ADRC method provided that they have the same observer poles.The fundamental difference lies in that the nominal system dynamics are fully exploited in the proposed observer design.This important feature allows for better performance balance between noise attenuation and disturbance rejection of the proposed method.

    2) Comparison With Existing Related Approaches:Since this paper is concerned with unmatched disturbance rejections via output-feedback control,it is worthy to note that most of the existing results focus on state-feedback control design[26]–[28] except for the active disturbance rejection control(ADRC).The potential difficulty for designing an outputfeedback disturbance rejection controller is that it is quite challenging to estimate both the unknown states and unmatched disturbances simultaneously.Consequently,we provide a thorough comparison with the traditional ADRC approach,which is a well-known output-feedback disturbance rejection controller,from different aspects including with a methodological perspective,delicate theoretical analysis and experimental comparisons.

    It should be mentioned that the preliminary version of the paper has been published in [29],where a very basic control design is considered for the buck-converter driven DC motor system subject to only load torque disturbance.The result therein is far from complete due to the following three reasons:1) The parameter uncertainties are not considered therein,which are inevitable for the system under consideration;2) The result is not complete theoretically since no rigorous performance analysis have been presented in that paper to show the relationship between the ultimate bounds and the observer parameters;3) There lacks a thorough theoretical analysis and experimental results to show that the proposed approach exhibits advantages over the traditional control approaches.Compared with [29],this paper enriches the results in both theory and experiment aspects.

    II.PRELIMINARIES AND MOTIVATIONS

    A.Model of Converter-Driven DC Motor System

    Fig.1.Circuit of the buck converter-driven DC motor system.

    The basic circuit of the buck converter-driven DC motor system considered in this paper is shown in Fig.1.In the buck converter circuit,there is a DC input voltageE,a switchVT,a diodeVD,a pulse width modulation (PWM) gate drive,a load resistorR,a capacitorCand an inductorL.The DC motor consists of an induced electromotive forceEa,an rotor winding resistorRaand an armature inductorLa.The nominal moment of the rotor inertia is denoted asJ.Based on Kirchhoff’s law,the system dynamic model is given by ofE,L,C,R,La,RaandJ,respectively.The first two formulas of (1) correspond to the dynamics of the buck converter,in whichiLandvoare the inductor current and capacitor voltage of converter,respectively,μ∈[0,1] is the duty ratio.The last two formulas represent the dynamics of the DC motor,in whichiaandωdenote the armature current and angular velocity of the motor,respectively,kmis motor torque constantkeis counter electromotive force constant,bdenotes the motor viscous friction coefficient,τLis the unknown time-varying load torque.Moreover,d1,d2,d3andd4are lumped disturbances denoted by

    wherefi,i=1,2,3,4 represents the rest unmodeled uncertainties in each channel.

    Remark 1:The unmodeled uncertainties in the converterdriven DC motor systems include the cogging torque,highorder back-electromotive force (EMF) harmonics,dead-time effects,etc.According to the existing results,the cogging torqueTcogis modeled in detail in [30] aswhere θm=θe/npis the mechanical angular position of the rotor,Qis the number of the slots,andTcogiand φiare the corresponding amplitude and phase angle.It is very difficult to model all the unmodeled uncertainties in detail,thus,we treat the unmodeled uncertainties as part of disturbances in this paper.Under the control framework proposed in this paper,the effects of disturbances can be estimated and compensated for in the controller.As a result,the effects from unmodeled uncertainties can be eliminated in the control performance.

    B.Research Motivation and Control Objective

    As mentioned in the introduction of the paper,the traditional ADRC and the generalized ESO-based control approach in [23] provide possible solutions for the problem under consideration.The potential limitations of these two approaches are elaborated in this subsection.

    1) Traditional ADRC:Generally,the traditional ADRC method for system (1) under consideration is designed based on the fourth-order integral chain system [22],where the observer is designed as

    whereκis a constant that depends on the circuit parameters,and ω?is the desired angular velocity.li,i=1,2,3,4,5 is the observer gain to design.ν1=ω?ω?,,i=1,2,3,4,5 denote estimations ofand the lumped disturbance,respectively.The ADRC law is designed as

    where φi,i=1,2,3,4 is the controller gain to be designed.It is assumed that the observer poles are set as ?ωo(ωo>0),and the observer gains in (3) are obtained and given by

    The traditional ADRC method can remove the offset caused by unmatched disturbances asymptotically but through the integral-chain coordinate transformation shown in [22],where system dynamics have not been fully exploited.For the fifthorder observer in (3),the observer gain shows polynomial growth in terms of the increase of the poles for the observer.It is natural that large observer gains will seriously amplify measurement noises and thereby cause abrupt changes in estimates as well as the control signal.

    2) Generalized ESO-Based Control:In [23],a generalized ESO-based control strategy is proposed for unmatched disturbance compensation.The design procedure is revisited as follows.

    Letx=[iL,vo,ia,ω]T.The compact form of the converterdriven DC motor system (1) is given by

    whereA,b,bdandCare matrices with compatible dimensions,andd=[d1,d2,d3,d4]T.Denotexˉ=[iL,vo,ia,ω,d1]T.The generalized ESO is designed as

    Following the guidelines given in [23],for the given converter-driven DC motor system (1),it is verified that the observability condition is not satisfied for the designed observer in (7) and thus the generalized ESO-based control approach is not available for the present control problem.

    Control Objective:This paper aims to design a new outputfeedback speed regulation controller for the converter-driven DC motor system (1) subject to various matched/unmatched uncertainties and disturbances.The controller can fully exploit the nominal model of the system given in (1) and only requires the output measurement.Under the proposed controller,the angular velocityωof the DC motor tracks the desired signal ω?asymptotically even in the presence of unmatched disturbances.

    III.MAIN RESULTS

    A.Controller Design

    By utilizing the basic principle of output regulation theory[24],the steady states of the buck converter-driven DC motor system can be calculated by solving the following regulator equations:

    where πiL,πvo,πia,πωand πμrepresent the desired steadystate responses ofiL,vo,ia,ωandμ,respectively.

    Remark 2:As per the output regulation principle [24],here the actual dynamic system (1) will be driven to its steady-state manifold governed by (9).It is clear that the steady-state manifold is naturally determined by the actual disturbances.This is the reason why the actual disturbances (2) appear in both the actual system (1) and the regulator equations (9).

    Define the following coordinate transformation:

    Taking the time-derivatives of (10) along (1) and (9),one obtains

    Since πμis a function with respect to disturbances,we denote πμas an extended state.Then the augmented system is given by

    where ?=dπμ/dt.The output of system (12) is the measurableeω.

    To begin with the observer design,we denote the state estimations in (12) asAn observer for system (12) is then designed as follows:

    where

    Remark 3:Since both the states and their steady-state responses are unknown,it is generally quite difficult,if not impossible,to estimate the unknown states and their desired steady-state responses simultaneously.Consequently,we propose a new manner to estimate error signals between the unknown system states and their desired steady-state responses directly by using adequate coordinate transformations.By virtue of coordinate transformations,an observer shown in (13) is proposed to identify the estimated parameters in the paper.

    By virtue of the estimations generated in (13),the composite controller is designed as follows:

    Remark 4:The controller can also be designed as

    sinceeωis available.In practice,if the sensor noises are severe and unacceptable,it is advised to choose the controller given in (14) since the observer can be taken as a low-pass filter as well [31].

    Remark 5:The solutions of the regulator equations (9) are not needed in implementation.From the coordinate transformation given in (10),the augmented system (12),the observer,and controller given by (13) and (14),it is known that although we cannot calculate,anddue to the presence of unknown disturbances,the controller can be implemented based on the observation technique.

    B.Performance Analysis

    The stability of the closed-loop buck converter-driven DC motor system is established by the following theorem.

    Theorem 1:For a bounded ?,if the observer gainLin (13)and the feedback control gain Φ in (14) are chosen such thatMandA11are Hurwitz,respectively.Then,both the estimation errors and the angular velocity tracking error will converge to a bounded neighborhood of the origin and the ultimate bounds can be assigned arbitrarily small by adjusting the poles of the observer.

    Proof:The proof of the theorem is divided into the following three steps.

    Step 1:By the definition of input-to-state stability (ISS)[32],there exist a class of KL function β1and a class of Kfunction β2such that for any bounded ?(t) and initial state E(0),the solution E(t) of system (16) satisfies‖E(t)‖≤β1(‖E(0),t‖)+β2(sup0≤τ≤t|?(τ)|).Moreover,ast→∞,there is‖E(∞)‖≤β2(δ)<∞ with |?(t)|≤δ,whereδis a positive constant.Therefore,the speed tracking error and the estimation errors are all bounded.

    Step 2:We assign the eigenvalues ofMandA11as?ωo(ωo>0) and ?ωc(ωc>0).The observer error dynamics are evaluated first,whose formation is depicted by

    To proceed with the analysis,we present the following statement:For a system

    A Lyapunov function is constructed asV(x)=xT Px.The first-order derivative ofV(x) is obtained as

    Through similar analysis,we get the other four estimation errors as follows:

    It is obvious that the ultimate bound of(i=1,2,...,5) can be arbitrarily small by magnifying ωo.

    Step 3:The tracking error subsystem is written as

    Next,the disturbance rejection ability of the closed-loop system under controller (14) is summarized as follows.

    Theorem 2:If ?(t) tends to zero as t tends to infinity,the observer gainLin (13) and the controller gain Φ in (14) are chosen such thatMandA11are Hurwitz,andis invertible withC=[0,0,0,1],then the angular velocity ω in(1) tracks the desired value ω?asymptotically in the presence of disturbances/uncertainties under the control law (14).

    Proof:From (17),under the given condition and Theorem 1,we have

    Fig.2.Control structure of the model-based output-feedback controller.

    IV.ExPERIMENTAL VERIFICATIONS

    To validate the performance of the proposed model-based control strategy,experiments on a buck converter-driven DC motor system are carried out in this section.

    A.Experimental Test Setup

    As shown in Fig.3,the experimental test setup contains a buck converter,a permanent magnet DC motor (DBM100),a tension controller,a magnetic powder brake device for regulating load torque,an incremental photoelectric encoder(A and B two-phase output,resolution of 1 000 lines),a programmable DC power source (EA-PSI9500-20 2U),and dSPACE DS1103 real-time controller.For experimental studies,the sampling frequency and the controller updating frequency is setting as 3.3 kHz.The fixed PWM switching frequency of the converter is 10 kHz.The algorithm is implemented by dSPACE DS1103 control board,which can achieve seamless connection with MATLAB/Simulink.Table I shows nominal values of parameters for the buck converterdriven DC motor circuit.

    Fig.3.Experimental prototype.

    TABLE IPARAMETERS OF BUCK CONVERTER-DRIVEN DC MOTOR SYSTEM

    B.Experimental Results and Analysis

    To demonstrate the efficiency of the proposed controller,the traditional ADRC method and the proportional integral derivative (PID) controller are employed for comparison for experimental studies.The traditional ADRC method is designed along (3) and (4).The PID controller is designed as

    Parameter tuning of different controller is definitely significant in providing fair performance comparisons.It is well known that,there are usually a number of conflicts/constraints in control design and parameter tuning,such as tracking versus disturbance rejection and nominal performance versus robustness.As such,it is generally impossible to design and tune a controller with all performances better than others.For fair comparison,one sensible option is to tune the controllers so that obtain similar performance specifications and then compare the others.Since this paper focuses on disturbance/uncertainties attenuation,for fair performance comparisons among the three controllers,the control parameters of each controller are tuned such that the speed tracking profiles have similar performance specifications,such as settling time and overshoot,in the absence of disturbances and uncertainties.Subsequently,the robustness and disturbance rejection performances are fairly compared.

    For traditional ADRC design,the feedback control and observer parameters are tuned by using pole placement approached [33].The PID controller is tuned by trail and error approach such that one can get satisfactory tracking control performance.To conclude,we aim to assign adequate control parameters of traditional ADRC and PID approaches such that one can obtain desirable static and dynamic speed reference tracking performance,such as short settling time and small overshoot.

    The controller poles of the proposed method and the ADRC are set as the same.In the experiment,the observer poles of the proposed method are chosen as–500,and the observer poles of the ADRC method are chosen as–300.(The reason why we choose different poles for the two observers will be discussed in the next section.) The parameters of the three controllers are shown in Table II.In this section,we present four cases of disturbances and uncertainties,where the detailed experimental scenarios are listed in Table III.

    TABLE IICONTROL PARAMETERS OF THREE CONTROL APPROACHES

    Case I– Time-Varying Reference Speed Tracking:In this case,the reference speed is expressed asω?=80+75sin(2πt)rad/s.For all three controllers,the higher-order timederivations of the reference signal are not utilized.The tracking response curves and duty ratio profiles are shown in Fig.4,which shows that,with similar level of duty ratios,the proposed control approach has much smaller speed tracking error as compared with the traditional ADRC and PID control approaches.

    Case II– Sudden Input Voltage Changes:The input voltage changes from 40 V to 30 V at 1.3 s and recovers to 40 V at 3.3 s.Fig.5 shows the response curves in this case.

    TABLE IIIExPERIMENTAL SETUP SCENARIOS

    Fig.4.Response curves under the proposed method (left),traditional ADRC (middle) and PID (right) for time-varying reference speed tracking (Case I) (top:angular velocity;middle:tracking error;bottom:duty ratio).

    Fig.5.Response curves under the proposed method (left),traditional ADRC (middle) and PID (right) with sudden input voltage change (Case II) (top:angular velocity;bottom:duty ratio).

    Case III– Sudden Load Torque Disturbance:When sudden load torques are imposed,the experimental results are shown in Fig.6.It can be seen from Fig.6 that the proposed control approach has obtained much better transient regulation performance than the other two control approaches.

    Fig.6.Response curves under the proposed method (left),traditional ADRC (middle) and PID (right) with sudden load torque disturbance (Case III) (top:angular velocity;bottom:duty ratio).

    Fig.7.Response curves under the proposed method (left),traditional ADRC (middle) and PID (right) with time-varying load torque disturbance (Case IV)(top:angular velocity;bottom:duty ratio).

    Case IV– Time-Varying Load Torque Disturbance:Different from Cases II and III where constant disturbances are considered,in this case,we aim to evaluate the influences of persistent time-varying external disturbances on the steadystate performance.Response curves of the system output and control signal are shown in Fig.7,from which it can be seen that the proposed control method achieves much smaller steady-state speed fluctuations than the traditional ADRC and PID control schemes.

    In order to compare control performance under different controllers in more detail,we evaluate the performance specification quantitatively.The performance indices including maximum angular velocity drop/raise(MAVD/MAVR) and maximum recovery time (MRT) are utilized for performance validation.The corresponding speed regulation performance indices are summarized in Table IV.

    Remark 6:It is observed from Fig.5 the control performance of ADRC is better than PID controller,but in Fig.6,the control performance of ADRC is worse than that of the PID controller.The reason is that Figs.5 and 6 show experimental results in the presence of sudden input voltage changes and sudden load torque disturbances,respectively,and different disturbances from different channels may have different effects on the control performance.We can actually obtain better control performance by choosing larger poles for ADRC,but larger poles will seriously amplify measurement noises and thereby cause abrupt changes of the estimates as well as the control signal.

    TABLE IVREGULATION PERFORMANCE INDICES IN THREE CASES

    Remark 7:It should be pointed out that there exists a difference between the theoretical analysis of (24) and (25)and the experimental results given in the following section,since there are many factors such as noises and control authorities affecting the control performance in practice.

    V.DISCUSSIONS

    In this Section,we will pay particular attention to observer gains comparisons between the proposed method and the traditional ADRC method.The relationship between the observer poles (set as ?ωo) and observer gains under the proposed model-based observer (13) are given as follows:

    wherea12=?1/L0,a15=?E0/L0,a21=1/C0,a22=1/R0C0,a23=?1/C0,a32=1/La0,a33=?/,a34=?ke/,a43=km/J0,anda44=?b/J0.For the traditional ADRC method,the relationship is shown by (5).By substituting the model parameters of the converter-driven DC motor into (31)and (5),one can obtain the relationships between the observer gains and poles for both methods numerically,which are shown in Fig.8 and Table V.

    As shown by Fig.8 and Table V,the observer gains of the traditional ADRC increase drastically as the poles of the observer increases.However,due to full utilization of the nominal model,when enlarge the poles of the observer,the observer gains of the proposed control approach increase much slower than the traditional ADRC approach.This provides great potential for measurement sensor noise attenuation,which will be elaborated on in the following experimental study.

    We set the speed reference signal as ω?=150 rad/s.Two different observer poles ωo=?500 and ωo=?300 are considered here for the traditional ADRC method,while other controller settings are the same as those in Section V for both ADRC and the proposed method.

    Fig.8.Relationship between observer gains and poles.

    The response curves of the speed,the observer outputs and duty ratio are shown by Fig.9.It can be obtained from the comparison results that with larger observer poles,the traditional ADRC can obtain faster speed tracking performance but larger overshoot and considerable sensor noise amplification.In particular,it can be observed from Fig.9(f) that the duty ratio drastically fluctuates with a traditional ADRC with poles ωo=?500.This supports the previous statement that the proposed strategy has an excellent property where an increase the observer poles does not significantly amplify sensor noise.Therefore,the proposed controller widens the bandwidth without bringing any adverse effects.

    Remark 8:For the experiments in last section,if we set the observer poles as–500 for traditional ADRC,then the duty ratio will have serious chatter especially in the case when the unknown time-varying load torques are imposed on the shaft.Therefore,we set observer poles as–300 in the traditional ADRC method to balance speed tracking performance and duty ratio chatter alleviation.

    VI.CONCLUSIONS

    This paper investigated the speed regulation problem for a buck converter-driven DC motor system with parametric uncertainties and unknown load torque disturbances.Using adequate coordinate transformations,a new observer that cantake full advantage of system nominal dynamics,has been developed to estimate the errors between the system states and their desired steady-state responses.An output-feedback disturbance rejection controller has been proposed thereafter with performance analysis.The benefit with our proposed method is that with the same observer poles,the observer within this method has much smaller observer gains than the traditional ADRC approach.This has been presented with corresponding analysis and discussion.The experiment tests validate that the proposed model-based disturbance rejection control approach has much better control performances including tracking precision and robustness with unknown unmatched disturbances than a PID controller and traditional ADRC method.

    TABLE VOBSERVER GAINS OF TWO CONTROLLERS UNDER DIFFERENT POLES

    Fig.9.Response curves when ω?=150 rad/s (a) Angular velocity;(b) Estimation error;(c) Estimation of ;(d) Estimation of ;(e) Estimation of ;(f) Control input μ).

    久久99热6这里只有精品| 亚洲经典国产精华液单| 成年av动漫网址| 级片在线观看| 中文字幕制服av| 婷婷色av中文字幕| 午夜福利成人在线免费观看| 中文字幕亚洲精品专区| 欧美高清性xxxxhd video| 亚洲不卡免费看| 国产精品永久免费网站| 国产探花在线观看一区二区| 不卡视频在线观看欧美| 久久国产乱子免费精品| 亚洲欧美日韩卡通动漫| 中文字幕制服av| 伊人久久精品亚洲午夜| 精品久久久久久久久亚洲| 欧美性猛交╳xxx乱大交人| 熟女人妻精品中文字幕| 日本三级黄在线观看| 国产精品伦人一区二区| 国产私拍福利视频在线观看| 亚洲18禁久久av| 亚洲美女视频黄频| 国产精品人妻久久久影院| 亚洲美女视频黄频| 五月玫瑰六月丁香| 日韩强制内射视频| 国产黄片视频在线免费观看| 国产成人a∨麻豆精品| 亚洲va在线va天堂va国产| 国产黄色视频一区二区在线观看 | 99视频精品全部免费 在线| 熟女电影av网| 91久久精品电影网| 波多野结衣高清无吗| 国产一区二区三区av在线| 国产私拍福利视频在线观看| 国产av在哪里看| 中文字幕精品亚洲无线码一区| 久久6这里有精品| 毛片一级片免费看久久久久| 性色avwww在线观看| 亚洲欧美日韩高清专用| 国产亚洲一区二区精品| 中文字幕亚洲精品专区| 亚洲在线观看片| 麻豆一二三区av精品| 亚洲无线观看免费| 免费av不卡在线播放| 中文亚洲av片在线观看爽| 亚洲av日韩在线播放| 2022亚洲国产成人精品| 成人国产麻豆网| 国产伦在线观看视频一区| 午夜福利网站1000一区二区三区| 高清av免费在线| 日韩亚洲欧美综合| 国产精品人妻久久久久久| 波多野结衣高清无吗| 国产精品爽爽va在线观看网站| 国产黄色视频一区二区在线观看 | 亚洲av免费在线观看| 亚洲高清免费不卡视频| 久久久a久久爽久久v久久| 精品熟女少妇av免费看| 欧美一区二区国产精品久久精品| 日韩 亚洲 欧美在线| 亚洲人成网站在线观看播放| 人妻夜夜爽99麻豆av| 亚洲av成人av| 中文字幕精品亚洲无线码一区| 国产成人午夜福利电影在线观看| 3wmmmm亚洲av在线观看| 午夜免费激情av| 久久久a久久爽久久v久久| av福利片在线观看| 熟女电影av网| 人妻少妇偷人精品九色| 国产精品一区www在线观看| 午夜精品在线福利| 人妻少妇偷人精品九色| 亚洲久久久久久中文字幕| 国产日韩欧美在线精品| 国产成人午夜福利电影在线观看| 少妇高潮的动态图| 国产女主播在线喷水免费视频网站 | 国产亚洲一区二区精品| 精品人妻一区二区三区麻豆| 免费av观看视频| 免费观看在线日韩| 中文欧美无线码| 超碰av人人做人人爽久久| 免费看日本二区| 久久久精品大字幕| 老女人水多毛片| 亚洲熟妇中文字幕五十中出| 淫秽高清视频在线观看| 一边摸一边抽搐一进一小说| 成人性生交大片免费视频hd| 国产69精品久久久久777片| 国产精品永久免费网站| 日韩欧美在线乱码| 在线天堂最新版资源| 日韩av不卡免费在线播放| 久久99热这里只有精品18| 久久久久久大精品| 精品免费久久久久久久清纯| 国产精华一区二区三区| 久久久国产成人精品二区| 联通29元200g的流量卡| 欧美激情久久久久久爽电影| 国产三级中文精品| 亚洲国产欧洲综合997久久,| 国产精品久久久久久精品电影小说 | 久久草成人影院| 少妇熟女aⅴ在线视频| 久久这里只有精品中国| 在线免费十八禁| 少妇的逼水好多| 村上凉子中文字幕在线| 一级毛片我不卡| 国产成人精品一,二区| 国内少妇人妻偷人精品xxx网站| 国产精品精品国产色婷婷| 一级av片app| 老司机影院毛片| 五月玫瑰六月丁香| 蜜桃亚洲精品一区二区三区| 亚洲第一区二区三区不卡| 精品久久久久久久久av| 精品久久久久久成人av| 91久久精品电影网| 国产成人a∨麻豆精品| 又黄又爽又刺激的免费视频.| 日本免费一区二区三区高清不卡| 日日摸夜夜添夜夜爱| 日韩欧美国产在线观看| 汤姆久久久久久久影院中文字幕 | 丰满少妇做爰视频| 亚洲人成网站在线观看播放| 久久久成人免费电影| 日韩亚洲欧美综合| 久久亚洲国产成人精品v| 久久久久精品久久久久真实原创| 国产精品一区www在线观看| av在线亚洲专区| 免费播放大片免费观看视频在线观看 | 99久久人妻综合| av.在线天堂| 久久久久久国产a免费观看| 国产精品伦人一区二区| 中文亚洲av片在线观看爽| 一级毛片电影观看 | 亚洲av中文av极速乱| 日本免费a在线| 国产在线一区二区三区精 | 嫩草影院入口| 国产色婷婷99| 国产精品乱码一区二三区的特点| 99热这里只有是精品在线观看| 国产久久久一区二区三区| 麻豆成人av视频| 亚洲最大成人中文| 亚洲中文字幕一区二区三区有码在线看| 日日啪夜夜撸| 免费大片18禁| 国产在视频线在精品| 久久久久久九九精品二区国产| 免费观看人在逋| 国产精品永久免费网站| 最近中文字幕2019免费版| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av一区综合| 韩国高清视频一区二区三区| 色网站视频免费| 久久精品综合一区二区三区| 亚洲精品一区蜜桃| 只有这里有精品99| 午夜精品在线福利| 看十八女毛片水多多多| 国产精品精品国产色婷婷| 亚洲国产精品国产精品| 村上凉子中文字幕在线| 国内精品宾馆在线| 日韩国内少妇激情av| 爱豆传媒免费全集在线观看| 如何舔出高潮| 男人舔奶头视频| 69av精品久久久久久| 青春草视频在线免费观看| 中文欧美无线码| 欧美xxxx黑人xx丫x性爽| 亚洲三级黄色毛片| 中文天堂在线官网| 国产精品久久久久久精品电影| 日韩中字成人| 国产大屁股一区二区在线视频| 黄色日韩在线| 高清日韩中文字幕在线| 97热精品久久久久久| 久久精品国产鲁丝片午夜精品| 成年女人看的毛片在线观看| 久久精品影院6| 久久人人爽人人片av| 国产激情偷乱视频一区二区| 国产精品综合久久久久久久免费| 日韩三级伦理在线观看| 春色校园在线视频观看| 精品久久久久久久久久久久久| 熟妇人妻久久中文字幕3abv| 日韩精品青青久久久久久| 国产高清视频在线观看网站| av卡一久久| 一级爰片在线观看| 神马国产精品三级电影在线观看| 亚洲最大成人av| 国产乱来视频区| 欧美潮喷喷水| 国产大屁股一区二区在线视频| 亚洲国产成人一精品久久久| 在线a可以看的网站| 亚洲国产精品专区欧美| 一本久久精品| 久久这里有精品视频免费| av在线观看视频网站免费| 欧美激情在线99| 亚洲国产成人一精品久久久| 小说图片视频综合网站| 亚洲精品国产成人久久av| 亚洲久久久久久中文字幕| 午夜福利在线观看免费完整高清在| 欧美一区二区国产精品久久精品| 中文字幕亚洲精品专区| 国内精品一区二区在线观看| 国产精品美女特级片免费视频播放器| 99久国产av精品| 三级经典国产精品| 国产精品国产三级国产av玫瑰| 好男人在线观看高清免费视频| 亚洲国产日韩欧美精品在线观看| 自拍偷自拍亚洲精品老妇| 69av精品久久久久久| 又黄又爽又刺激的免费视频.| 91久久精品电影网| 青春草亚洲视频在线观看| 亚洲人成网站在线播| 全区人妻精品视频| 亚洲熟妇中文字幕五十中出| 最近的中文字幕免费完整| 麻豆成人av视频| 日日撸夜夜添| 人妻系列 视频| 在线天堂最新版资源| 日本免费a在线| 日韩欧美国产在线观看| 久久欧美精品欧美久久欧美| 免费看美女性在线毛片视频| 亚洲国产成人一精品久久久| 九九在线视频观看精品| 国产极品天堂在线| 欧美最新免费一区二区三区| 精品国产三级普通话版| 亚洲欧美日韩卡通动漫| 干丝袜人妻中文字幕| 少妇高潮的动态图| 精华霜和精华液先用哪个| 尾随美女入室| 亚洲欧洲国产日韩| av女优亚洲男人天堂| 少妇熟女aⅴ在线视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品久久久久久精品电影小说 | 国产色婷婷99| 中文欧美无线码| 噜噜噜噜噜久久久久久91| 内地一区二区视频在线| 欧美最新免费一区二区三区| 菩萨蛮人人尽说江南好唐韦庄 | 级片在线观看| 成人亚洲精品av一区二区| 又黄又爽又刺激的免费视频.| 久久精品91蜜桃| 婷婷六月久久综合丁香| 午夜老司机福利剧场| 又黄又爽又刺激的免费视频.| 麻豆久久精品国产亚洲av| 床上黄色一级片| 成人鲁丝片一二三区免费| 看十八女毛片水多多多| 亚洲精品色激情综合| 亚洲成人久久爱视频| 高清视频免费观看一区二区 | 亚洲精品久久久久久婷婷小说 | 高清毛片免费看| 成人亚洲欧美一区二区av| 99在线人妻在线中文字幕| 男人狂女人下面高潮的视频| .国产精品久久| 床上黄色一级片| 免费看日本二区| 国产免费又黄又爽又色| 精品人妻视频免费看| 国模一区二区三区四区视频| 午夜视频国产福利| 国产不卡一卡二| 日韩亚洲欧美综合| 国产精品乱码一区二三区的特点| 国产亚洲91精品色在线| 欧美另类亚洲清纯唯美| av在线蜜桃| 日本免费a在线| 十八禁国产超污无遮挡网站| 亚洲一级一片aⅴ在线观看| 老司机福利观看| 欧美日韩在线观看h| 少妇人妻精品综合一区二区| 男女啪啪激烈高潮av片| 亚洲成人久久爱视频| 身体一侧抽搐| 99久国产av精品国产电影| 国产老妇女一区| 少妇猛男粗大的猛烈进出视频 | 国产亚洲精品av在线| 国产精品三级大全| 少妇的逼水好多| 国产成人福利小说| 国产精品野战在线观看| 男女国产视频网站| 日本一本二区三区精品| 亚洲欧美日韩无卡精品| 成人午夜高清在线视频| 97超视频在线观看视频| av国产久精品久网站免费入址| 亚洲国产精品国产精品| 男的添女的下面高潮视频| 国产成人a∨麻豆精品| 久久精品熟女亚洲av麻豆精品 | 国产亚洲精品av在线| 五月伊人婷婷丁香| 日本欧美国产在线视频| 国产成人一区二区在线| 欧美一区二区国产精品久久精品| 综合色丁香网| 久久久国产成人精品二区| 永久免费av网站大全| 国模一区二区三区四区视频| 亚洲精品aⅴ在线观看| 久久精品影院6| 久久精品国产99精品国产亚洲性色| 小说图片视频综合网站| 亚洲精品国产成人久久av| 亚洲美女搞黄在线观看| 嫩草影院精品99| 亚洲美女搞黄在线观看| av免费观看日本| 好男人在线观看高清免费视频| 免费看光身美女| АⅤ资源中文在线天堂| 老司机福利观看| 精品人妻偷拍中文字幕| 一个人观看的视频www高清免费观看| 青春草视频在线免费观看| 黄片wwwwww| 国产真实乱freesex| 免费无遮挡裸体视频| 国产一区二区在线av高清观看| 欧美极品一区二区三区四区| 久久热精品热| 国模一区二区三区四区视频| 久久久久免费精品人妻一区二区| 免费av毛片视频| 三级国产精品片| 亚洲电影在线观看av| 一区二区三区四区激情视频| 99热精品在线国产| 日本wwww免费看| 美女内射精品一级片tv| 水蜜桃什么品种好| 国产精品永久免费网站| 色5月婷婷丁香| 免费电影在线观看免费观看| 亚洲一区高清亚洲精品| 久久精品91蜜桃| 两个人视频免费观看高清| 美女被艹到高潮喷水动态| 国产成人91sexporn| 亚洲欧美精品综合久久99| 久久精品人妻少妇| 亚洲国产精品成人综合色| 日日摸夜夜添夜夜添av毛片| 爱豆传媒免费全集在线观看| 女的被弄到高潮叫床怎么办| 22中文网久久字幕| 午夜老司机福利剧场| 91狼人影院| 国产精品人妻久久久影院| 老司机福利观看| 国产伦理片在线播放av一区| www.av在线官网国产| 男女边吃奶边做爰视频| 国产大屁股一区二区在线视频| 国产精品国产三级专区第一集| 亚洲国产精品成人久久小说| 亚洲第一区二区三区不卡| 男女视频在线观看网站免费| 亚洲性久久影院| 国产精品国产三级国产av玫瑰| 女的被弄到高潮叫床怎么办| 日韩精品青青久久久久久| 亚洲四区av| 免费无遮挡裸体视频| 校园人妻丝袜中文字幕| 成人鲁丝片一二三区免费| 国产一区有黄有色的免费视频 | 嫩草影院入口| 中国美白少妇内射xxxbb| 精品国产三级普通话版| 91精品国产九色| 亚洲精品乱码久久久久久按摩| 久99久视频精品免费| 永久免费av网站大全| 亚洲av免费高清在线观看| 日本免费在线观看一区| 能在线免费看毛片的网站| 亚洲国产精品合色在线| 青春草亚洲视频在线观看| 天堂网av新在线| 七月丁香在线播放| 一卡2卡三卡四卡精品乱码亚洲| 日本五十路高清| 中文乱码字字幕精品一区二区三区 | 国产麻豆成人av免费视频| 人体艺术视频欧美日本| 午夜老司机福利剧场| 精品国产露脸久久av麻豆 | 七月丁香在线播放| 亚洲人成网站在线播| 午夜福利在线在线| 一本一本综合久久| 国产精品一区二区在线观看99 | 久久精品夜色国产| 国内精品一区二区在线观看| 乱人视频在线观看| 青春草视频在线免费观看| 日本wwww免费看| 韩国高清视频一区二区三区| 日韩欧美 国产精品| 18禁裸乳无遮挡免费网站照片| 少妇高潮的动态图| 在线观看美女被高潮喷水网站| 秋霞伦理黄片| 国产精品人妻久久久影院| 不卡视频在线观看欧美| 国产视频首页在线观看| 免费人成在线观看视频色| 国产精品一区二区三区四区久久| 高清在线视频一区二区三区 | 欧美日韩综合久久久久久| 国国产精品蜜臀av免费| 十八禁国产超污无遮挡网站| 亚洲伊人久久精品综合 | 亚洲电影在线观看av| 一区二区三区乱码不卡18| 亚洲经典国产精华液单| 女人十人毛片免费观看3o分钟| 赤兔流量卡办理| 亚洲国产欧洲综合997久久,| 天堂影院成人在线观看| 亚洲av成人精品一二三区| 久久草成人影院| 村上凉子中文字幕在线| a级毛色黄片| 欧美又色又爽又黄视频| 亚洲欧美一区二区三区国产| 午夜福利网站1000一区二区三区| 免费观看的影片在线观看| 国产高清不卡午夜福利| 嫩草影院入口| 久久久久久久久久久免费av| 国产高清有码在线观看视频| 亚洲欧美日韩东京热| 视频中文字幕在线观看| 又黄又爽又刺激的免费视频.| 2021天堂中文幕一二区在线观| 蜜桃亚洲精品一区二区三区| 黑人高潮一二区| 欧美日本视频| 精品久久久久久久末码| 久久久久久九九精品二区国产| 亚洲一区高清亚洲精品| 最近最新中文字幕免费大全7| 日韩一区二区视频免费看| 乱人视频在线观看| 欧美日韩综合久久久久久| 国产黄色视频一区二区在线观看 | 男女那种视频在线观看| 亚洲av电影在线观看一区二区三区 | 亚洲精品影视一区二区三区av| 国产熟女欧美一区二区| 日韩一区二区三区影片| 精品久久久久久电影网 | 91久久精品国产一区二区成人| 天堂中文最新版在线下载 | 观看美女的网站| 国产精品国产高清国产av| 日韩亚洲欧美综合| 日本午夜av视频| 久久久精品大字幕| 69人妻影院| 大话2 男鬼变身卡| 99久国产av精品| 国产精品人妻久久久影院| 国产成人一区二区在线| 99久久精品国产国产毛片| 国产精品国产三级专区第一集| 亚洲精品一区蜜桃| 激情 狠狠 欧美| 最近中文字幕2019免费版| 秋霞在线观看毛片| 国产成人免费观看mmmm| 亚洲国产日韩欧美精品在线观看| 国语对白做爰xxxⅹ性视频网站| 午夜视频国产福利| 一个人观看的视频www高清免费观看| 日韩欧美 国产精品| 亚洲av不卡在线观看| 亚洲国产欧洲综合997久久,| 日韩,欧美,国产一区二区三区 | 国产 一区精品| 国产日韩欧美在线精品| 狠狠狠狠99中文字幕| 1000部很黄的大片| 精品99又大又爽又粗少妇毛片| 少妇的逼好多水| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲欧美日韩高清专用| 看非洲黑人一级黄片| av专区在线播放| av卡一久久| 亚洲精品日韩在线中文字幕| 国产在线一区二区三区精 | 免费观看精品视频网站| 欧美zozozo另类| 久久热精品热| 精品久久久久久久久久久久久| 中国国产av一级| 高清午夜精品一区二区三区| 国产单亲对白刺激| 又粗又硬又长又爽又黄的视频| 看免费成人av毛片| 国产成人aa在线观看| 草草在线视频免费看| 中文字幕av在线有码专区| 国产精品久久视频播放| 日本五十路高清| 男女那种视频在线观看| 久久精品久久精品一区二区三区| 1000部很黄的大片| 精品久久久久久成人av| 淫秽高清视频在线观看| 国产精品久久久久久久电影| 又爽又黄无遮挡网站| 亚洲自偷自拍三级| 国产精品国产高清国产av| 欧美xxxx性猛交bbbb| 久久久久久久久大av| 91狼人影院| 老司机影院成人| 免费看a级黄色片| 亚洲久久久久久中文字幕| 在线免费观看的www视频| 秋霞伦理黄片| 看十八女毛片水多多多| 日本av手机在线免费观看| 一级黄色大片毛片| 亚洲欧美日韩高清专用| av播播在线观看一区| av黄色大香蕉| 国产亚洲最大av| 女人久久www免费人成看片 | 伊人久久精品亚洲午夜| 日本一二三区视频观看| 亚洲真实伦在线观看| 欧美高清性xxxxhd video| 国产男人的电影天堂91| 国产伦一二天堂av在线观看| 午夜福利在线观看吧| 日韩av在线免费看完整版不卡| 极品教师在线视频| 日韩大片免费观看网站 | 黄片wwwwww| 亚洲在线观看片| 美女xxoo啪啪120秒动态图| 白带黄色成豆腐渣| 99久国产av精品国产电影| 亚洲精品亚洲一区二区| 日韩一区二区三区影片| 男人舔女人下体高潮全视频| 又黄又爽又刺激的免费视频.| 久久久久久久午夜电影| 亚洲最大成人中文| 午夜免费激情av| 波多野结衣高清无吗| 精品99又大又爽又粗少妇毛片| 亚洲av免费在线观看| 亚洲成av人片在线播放无| 狠狠狠狠99中文字幕| 国产精品久久久久久久电影| 非洲黑人性xxxx精品又粗又长| 少妇丰满av| 99热精品在线国产| 一个人免费在线观看电影| 熟女电影av网|