• <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 μ).

    97超视频在线观看视频| 国内毛片毛片毛片毛片毛片| 国语自产精品视频在线第100页| 级片在线观看| 老司机午夜十八禁免费视频| 成人鲁丝片一二三区免费| 国产伦人伦偷精品视频| 国产真人三级小视频在线观看| 一区福利在线观看| 一进一出好大好爽视频| 无遮挡黄片免费观看| 国产视频一区二区在线看| 每晚都被弄得嗷嗷叫到高潮| 麻豆成人av在线观看| 少妇熟女aⅴ在线视频| 国产精品久久久久久久电影 | 一二三四社区在线视频社区8| 欧美午夜高清在线| 欧美精品啪啪一区二区三区| 日本黄色片子视频| 久久久国产成人免费| 久久精品综合一区二区三区| 国产一区二区亚洲精品在线观看| 90打野战视频偷拍视频| 九九热线精品视视频播放| 天美传媒精品一区二区| 最近最新免费中文字幕在线| 男人和女人高潮做爰伦理| 一卡2卡三卡四卡精品乱码亚洲| 亚洲不卡免费看| 最新美女视频免费是黄的| 亚洲av成人精品一区久久| 啦啦啦观看免费观看视频高清| 无限看片的www在线观看| 国产一级毛片七仙女欲春2| 在线观看免费视频日本深夜| 男人舔奶头视频| 97碰自拍视频| 国产精品香港三级国产av潘金莲| 国产男靠女视频免费网站| 国产不卡一卡二| 波多野结衣高清作品| 亚洲av电影在线进入| 日日干狠狠操夜夜爽| 精品久久久久久久末码| 国产精品亚洲美女久久久| 午夜久久久久精精品| 免费观看精品视频网站| 亚洲无线观看免费| 99视频精品全部免费 在线| 中文字幕av成人在线电影| 亚洲精品日韩av片在线观看 | 母亲3免费完整高清在线观看| 麻豆国产av国片精品| 色视频www国产| 丰满人妻熟妇乱又伦精品不卡| 国产精品电影一区二区三区| 国产精品99久久久久久久久| 女同久久另类99精品国产91| 日本免费一区二区三区高清不卡| 国产男靠女视频免费网站| 18禁国产床啪视频网站| 国产一区二区三区在线臀色熟女| 欧美成人一区二区免费高清观看| 18禁黄网站禁片午夜丰满| 亚洲av不卡在线观看| 国产高清三级在线| 日韩国内少妇激情av| 动漫黄色视频在线观看| 色在线成人网| 97超视频在线观看视频| 日本a在线网址| 真人一进一出gif抽搐免费| 免费看日本二区| 亚洲人成网站高清观看| 亚洲欧美日韩东京热| 12—13女人毛片做爰片一| 一级a爱片免费观看的视频| 亚洲精品色激情综合| 亚洲av一区综合| 国产亚洲精品久久久com| 美女高潮的动态| 男女那种视频在线观看| 亚洲va日本ⅴa欧美va伊人久久| 欧美三级亚洲精品| 欧美bdsm另类| 久久久久久国产a免费观看| 日韩av在线大香蕉| 国产精品美女特级片免费视频播放器| 精品久久久久久成人av| 真人做人爱边吃奶动态| 国产精品99久久99久久久不卡| 午夜福利在线在线| 免费电影在线观看免费观看| 国内毛片毛片毛片毛片毛片| 日韩亚洲欧美综合| 欧美成人免费av一区二区三区| 人妻丰满熟妇av一区二区三区| 欧美日韩瑟瑟在线播放| 精品电影一区二区在线| 九九久久精品国产亚洲av麻豆| 久久久久久久午夜电影| 欧美成人一区二区免费高清观看| 亚洲第一欧美日韩一区二区三区| 美女被艹到高潮喷水动态| 757午夜福利合集在线观看| 91字幕亚洲| 在线免费观看的www视频| 一个人免费在线观看电影| 成人性生交大片免费视频hd| 一卡2卡三卡四卡精品乱码亚洲| 少妇的逼水好多| 欧美三级亚洲精品| 国内揄拍国产精品人妻在线| 九九热线精品视视频播放| 男人的好看免费观看在线视频| 又爽又黄无遮挡网站| 国产精品野战在线观看| 在线a可以看的网站| 日本与韩国留学比较| 亚洲精品美女久久久久99蜜臀| 欧美在线一区亚洲| 老司机午夜十八禁免费视频| 啦啦啦观看免费观看视频高清| 精品福利观看| 男人舔奶头视频| 久久中文看片网| 亚洲欧美日韩卡通动漫| 热99re8久久精品国产| 级片在线观看| 国产精品国产高清国产av| 免费看a级黄色片| 亚洲精品国产精品久久久不卡| 午夜福利在线观看吧| 欧美黄色淫秽网站| 中文字幕av在线有码专区| 亚洲中文字幕一区二区三区有码在线看| 欧美bdsm另类| 国产精品女同一区二区软件 | 高清日韩中文字幕在线| 亚洲欧美一区二区三区黑人| 三级毛片av免费| 人妻丰满熟妇av一区二区三区| 一级黄色大片毛片| 亚洲成人精品中文字幕电影| 国产成人福利小说| 俺也久久电影网| 成熟少妇高潮喷水视频| 日韩欧美在线乱码| ponron亚洲| 国产精品免费一区二区三区在线| 亚洲欧美日韩卡通动漫| 久久久国产成人精品二区| 真人做人爱边吃奶动态| 身体一侧抽搐| www.色视频.com| 91在线精品国自产拍蜜月 | 国产av在哪里看| 久久精品国产清高在天天线| 女人十人毛片免费观看3o分钟| 国产毛片a区久久久久| 婷婷精品国产亚洲av在线| 亚洲av成人av| 亚洲精品久久国产高清桃花| 人妻久久中文字幕网| 又黄又粗又硬又大视频| 欧美+亚洲+日韩+国产| 久久精品影院6| 久久九九热精品免费| 午夜久久久久精精品| 精品福利观看| 在线观看一区二区三区| 免费看美女性在线毛片视频| 无人区码免费观看不卡| 看免费av毛片| 亚洲电影在线观看av| 少妇丰满av| 欧美激情在线99| 欧美中文日本在线观看视频| 波多野结衣巨乳人妻| 我的老师免费观看完整版| 88av欧美| 毛片女人毛片| 亚洲在线自拍视频| 精品国产超薄肉色丝袜足j| 国产伦精品一区二区三区视频9 | 欧美国产日韩亚洲一区| 日韩精品中文字幕看吧| 日韩欧美精品v在线| 精品久久久久久久毛片微露脸| 99精品久久久久人妻精品| 久久国产精品人妻蜜桃| 全区人妻精品视频| 精品国产美女av久久久久小说| 免费人成视频x8x8入口观看| 久久草成人影院| 国产伦在线观看视频一区| 波多野结衣高清无吗| 国产亚洲欧美在线一区二区| 午夜精品久久久久久毛片777| 淫妇啪啪啪对白视频| 中文亚洲av片在线观看爽| 宅男免费午夜| 俺也久久电影网| 真实男女啪啪啪动态图| 女同久久另类99精品国产91| 国产精品一及| 日韩欧美精品免费久久 | 看片在线看免费视频| 亚洲成人久久性| 亚洲18禁久久av| 久9热在线精品视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日本免费a在线| 老鸭窝网址在线观看| 久久午夜亚洲精品久久| 18禁黄网站禁片午夜丰满| 又爽又黄无遮挡网站| 日本在线视频免费播放| 欧美一级a爱片免费观看看| 真人做人爱边吃奶动态| 在线十欧美十亚洲十日本专区| 欧美区成人在线视频| 91字幕亚洲| 一a级毛片在线观看| 美女高潮喷水抽搐中文字幕| 久久久久久大精品| 最近最新中文字幕大全免费视频| 亚洲av免费高清在线观看| 免费在线观看亚洲国产| 亚洲国产欧美网| 国产三级在线视频| 男女之事视频高清在线观看| 蜜桃久久精品国产亚洲av| 免费av不卡在线播放| 香蕉av资源在线| 禁无遮挡网站| 日日摸夜夜添夜夜添小说| 国产亚洲精品久久久久久毛片| 中文字幕av成人在线电影| 国产一区二区亚洲精品在线观看| 久9热在线精品视频| av黄色大香蕉| 婷婷精品国产亚洲av| 91麻豆精品激情在线观看国产| 国模一区二区三区四区视频| xxx96com| 一区二区三区高清视频在线| 成人特级av手机在线观看| 又爽又黄无遮挡网站| 欧美日韩精品网址| 91麻豆精品激情在线观看国产| 真人一进一出gif抽搐免费| 成人精品一区二区免费| 日本五十路高清| 99国产精品一区二区蜜桃av| 内地一区二区视频在线| 久久香蕉精品热| 九色成人免费人妻av| 美女 人体艺术 gogo| 99国产精品一区二区三区| 国产精品久久久人人做人人爽| 久久午夜亚洲精品久久| 人人妻人人看人人澡| 亚洲成人中文字幕在线播放| 母亲3免费完整高清在线观看| 亚洲人成网站高清观看| 又粗又爽又猛毛片免费看| 国产精品精品国产色婷婷| 精品99又大又爽又粗少妇毛片 | 人妻久久中文字幕网| 熟女电影av网| 亚洲av日韩精品久久久久久密| 免费看日本二区| 国产老妇女一区| 亚洲精品国产精品久久久不卡| 欧美黄色淫秽网站| 在线观看66精品国产| 国产精品一区二区三区四区久久| 国产麻豆成人av免费视频| 天天一区二区日本电影三级| 色噜噜av男人的天堂激情| 最后的刺客免费高清国语| 啦啦啦韩国在线观看视频| 成人无遮挡网站| av黄色大香蕉| 好男人电影高清在线观看| 舔av片在线| 欧洲精品卡2卡3卡4卡5卡区| 人妻夜夜爽99麻豆av| 嫩草影院精品99| 精品99又大又爽又粗少妇毛片 | www.熟女人妻精品国产| 一进一出抽搐gif免费好疼| av在线天堂中文字幕| 欧美极品一区二区三区四区| 国产亚洲精品久久久久久毛片| 国产精品日韩av在线免费观看| 国产精品亚洲av一区麻豆| 婷婷精品国产亚洲av| 久久香蕉国产精品| 亚洲av不卡在线观看| 身体一侧抽搐| 免费观看人在逋| 久久人妻av系列| 午夜免费成人在线视频| 免费人成视频x8x8入口观看| 啦啦啦观看免费观看视频高清| 村上凉子中文字幕在线| 18禁裸乳无遮挡免费网站照片| 丁香六月欧美| 国产精品香港三级国产av潘金莲| 麻豆国产av国片精品| 一个人看视频在线观看www免费 | 亚洲av成人精品一区久久| 免费看a级黄色片| 在线观看午夜福利视频| 久久久久久久久中文| 日本 av在线| 久9热在线精品视频| 国产精品乱码一区二三区的特点| 日本精品一区二区三区蜜桃| 欧美色欧美亚洲另类二区| 小说图片视频综合网站| 中文字幕熟女人妻在线| 久久伊人香网站| 日韩免费av在线播放| 色播亚洲综合网| a级一级毛片免费在线观看| 啦啦啦免费观看视频1| 好男人在线观看高清免费视频| 国产精品亚洲美女久久久| 亚洲,欧美精品.| 亚洲一区高清亚洲精品| 午夜两性在线视频| 偷拍熟女少妇极品色| 蜜桃久久精品国产亚洲av| 大型黄色视频在线免费观看| 日韩av在线大香蕉| 天堂av国产一区二区熟女人妻| eeuss影院久久| 久久久国产成人免费| 婷婷亚洲欧美| 欧美日韩中文字幕国产精品一区二区三区| 亚洲av不卡在线观看| 一区二区三区高清视频在线| 亚洲黑人精品在线| 中文字幕高清在线视频| 色综合亚洲欧美另类图片| aaaaa片日本免费| 国产精品久久久久久久电影 | 国产免费男女视频| 一二三四社区在线视频社区8| 精品日产1卡2卡| 在线观看av片永久免费下载| 国产精品 国内视频| 亚洲午夜理论影院| 精品乱码久久久久久99久播| 最新在线观看一区二区三区| 不卡一级毛片| 最近最新中文字幕大全免费视频| 丁香六月欧美| 他把我摸到了高潮在线观看| 一级黄片播放器| 亚洲 欧美 日韩 在线 免费| 精品免费久久久久久久清纯| 亚洲人成网站高清观看| 亚洲乱码一区二区免费版| 舔av片在线| 欧洲精品卡2卡3卡4卡5卡区| 国产av在哪里看| 精品欧美国产一区二区三| 国产精品国产高清国产av| 久久久久久大精品| 黑人欧美特级aaaaaa片| 五月伊人婷婷丁香| 在线免费观看不下载黄p国产 | 欧洲精品卡2卡3卡4卡5卡区| 国产男靠女视频免费网站| 一本一本综合久久| 一个人看的www免费观看视频| av在线天堂中文字幕| 真实男女啪啪啪动态图| 国产乱人视频| 久久亚洲真实| netflix在线观看网站| 国产精品一区二区三区四区久久| 国产精品久久电影中文字幕| 国产精品久久久久久亚洲av鲁大| 美女cb高潮喷水在线观看| 少妇高潮的动态图| 俄罗斯特黄特色一大片| 99热这里只有是精品50| 一a级毛片在线观看| 美女 人体艺术 gogo| 欧美日韩黄片免| 中亚洲国语对白在线视频| 国产亚洲精品一区二区www| 国产精品一区二区三区四区久久| 亚洲不卡免费看| 国产三级在线视频| 熟女人妻精品中文字幕| 日本黄色视频三级网站网址| 色播亚洲综合网| 国产伦一二天堂av在线观看| 国产精品日韩av在线免费观看| 黑人欧美特级aaaaaa片| 久久亚洲精品不卡| 美女高潮的动态| 嫩草影视91久久| 99在线视频只有这里精品首页| 国产精品亚洲一级av第二区| 51午夜福利影视在线观看| 午夜视频国产福利| 美女高潮喷水抽搐中文字幕| 欧美国产日韩亚洲一区| www日本黄色视频网| 国产亚洲欧美98| 97碰自拍视频| 午夜精品久久久久久毛片777| 国产精品久久久久久久电影 | 国产老妇女一区| 日韩av在线大香蕉| 特级一级黄色大片| 色综合站精品国产| 国产成人av激情在线播放| 他把我摸到了高潮在线观看| 欧美极品一区二区三区四区| 久久国产乱子伦精品免费另类| 操出白浆在线播放| 国内精品美女久久久久久| 国产成人影院久久av| 亚洲av二区三区四区| 成人午夜高清在线视频| 免费av观看视频| 99在线人妻在线中文字幕| 男女那种视频在线观看| 国产欧美日韩一区二区三| 成熟少妇高潮喷水视频| 亚洲国产精品999在线| av片东京热男人的天堂| 亚洲专区中文字幕在线| 精品一区二区三区视频在线 | 久9热在线精品视频| 亚洲 欧美 日韩 在线 免费| 国产三级黄色录像| 国内精品久久久久精免费| 亚洲激情在线av| 成人18禁在线播放| 一区福利在线观看| 国产成人啪精品午夜网站| 精品人妻偷拍中文字幕| 午夜精品久久久久久毛片777| 午夜福利在线观看免费完整高清在 | 色av中文字幕| 久久人人精品亚洲av| 国产中年淑女户外野战色| 大型黄色视频在线免费观看| 成人性生交大片免费视频hd| 两个人的视频大全免费| 国产精品嫩草影院av在线观看 | 亚洲成av人片免费观看| 国产精品久久久久久久电影 | 动漫黄色视频在线观看| 十八禁人妻一区二区| 精品福利观看| 老鸭窝网址在线观看| 欧美日韩精品网址| av视频在线观看入口| 欧美精品啪啪一区二区三区| 精品免费久久久久久久清纯| 又黄又粗又硬又大视频| 免费观看人在逋| 无限看片的www在线观看| 超碰av人人做人人爽久久 | 男人舔女人下体高潮全视频| 狂野欧美激情性xxxx| 午夜影院日韩av| 性欧美人与动物交配| 色在线成人网| 免费看日本二区| 亚洲成人中文字幕在线播放| 蜜桃亚洲精品一区二区三区| 美女免费视频网站| 人妻久久中文字幕网| 色噜噜av男人的天堂激情| 精品99又大又爽又粗少妇毛片 | 免费一级毛片在线播放高清视频| 19禁男女啪啪无遮挡网站| 亚洲在线观看片| 亚洲精品色激情综合| 天堂网av新在线| 一二三四社区在线视频社区8| 99热精品在线国产| 久久久久久国产a免费观看| 成人精品一区二区免费| 淫秽高清视频在线观看| 欧美一级a爱片免费观看看| 亚洲av成人av| 丁香六月欧美| 床上黄色一级片| 一区二区三区国产精品乱码| 成年女人永久免费观看视频| 又爽又黄无遮挡网站| 亚洲 国产 在线| 很黄的视频免费| 亚洲一区二区三区色噜噜| 国产伦人伦偷精品视频| 日韩国内少妇激情av| 亚洲国产精品成人综合色| 免费av不卡在线播放| a在线观看视频网站| 国产一级毛片七仙女欲春2| 午夜精品久久久久久毛片777| 亚洲专区中文字幕在线| aaaaa片日本免费| 亚洲国产精品999在线| 午夜福利成人在线免费观看| 国产主播在线观看一区二区| 国产高清videossex| 香蕉av资源在线| 亚洲在线观看片| 琪琪午夜伦伦电影理论片6080| 亚洲国产欧美网| 国产黄色小视频在线观看| 亚洲国产精品合色在线| 最近最新免费中文字幕在线| 国产精品久久久久久久电影 | 大型黄色视频在线免费观看| 欧美色欧美亚洲另类二区| 91av网一区二区| 久久伊人香网站| 岛国视频午夜一区免费看| 久久久精品大字幕| 丰满乱子伦码专区| 特级一级黄色大片| 国产精品 国内视频| 精品国产亚洲在线| 可以在线观看的亚洲视频| 亚洲精品亚洲一区二区| av片东京热男人的天堂| www日本黄色视频网| www.熟女人妻精品国产| 男女床上黄色一级片免费看| 真人做人爱边吃奶动态| 我要搜黄色片| 婷婷精品国产亚洲av| 亚洲最大成人手机在线| 一进一出抽搐gif免费好疼| 97超级碰碰碰精品色视频在线观看| 无遮挡黄片免费观看| 法律面前人人平等表现在哪些方面| 免费电影在线观看免费观看| 啦啦啦观看免费观看视频高清| 亚洲av电影在线进入| 757午夜福利合集在线观看| 日韩欧美精品免费久久 | 欧美在线黄色| 色尼玛亚洲综合影院| 亚洲国产欧洲综合997久久,| 少妇的逼水好多| 欧美成人性av电影在线观看| 在线观看免费视频日本深夜| 欧美成狂野欧美在线观看| 国内少妇人妻偷人精品xxx网站| 色视频www国产| 九九久久精品国产亚洲av麻豆| 夜夜躁狠狠躁天天躁| 天天添夜夜摸| 欧美日韩精品网址| 国产精品嫩草影院av在线观看 | 久久久国产成人精品二区| 久久99热这里只有精品18| 啦啦啦韩国在线观看视频| 免费看美女性在线毛片视频| h日本视频在线播放| 麻豆国产av国片精品| 色在线成人网| 尤物成人国产欧美一区二区三区| 国产精品一区二区三区四区久久| 亚洲国产欧美网| 免费人成视频x8x8入口观看| 动漫黄色视频在线观看| 国产中年淑女户外野战色| 少妇丰满av| 欧美一级a爱片免费观看看| 俄罗斯特黄特色一大片| av福利片在线观看| 欧美成狂野欧美在线观看| av中文乱码字幕在线| 亚洲精品乱码久久久v下载方式 | 真人做人爱边吃奶动态| 国产高清有码在线观看视频| 久久精品91蜜桃| xxxwww97欧美| 国产蜜桃级精品一区二区三区| 窝窝影院91人妻| 91在线精品国自产拍蜜月 | 最近在线观看免费完整版| 亚洲真实伦在线观看| 校园春色视频在线观看| 99精品欧美一区二区三区四区| 色综合婷婷激情| 亚洲精华国产精华精| 日韩人妻高清精品专区| 欧美另类亚洲清纯唯美| 97超视频在线观看视频| 看黄色毛片网站| 嫩草影院入口| 亚洲专区中文字幕在线| 精品久久久久久成人av| 国产一区二区在线av高清观看| 亚洲欧美日韩东京热| 最新在线观看一区二区三区| 手机成人av网站|