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

    亚洲真实伦在线观看| 搡老妇女老女人老熟妇| 在线观看免费午夜福利视频| 色尼玛亚洲综合影院| 日本熟妇午夜| 久久国产乱子伦精品免费另类| 曰老女人黄片| 成人精品一区二区免费| 自线自在国产av| 国产精品香港三级国产av潘金莲| av在线天堂中文字幕| 国产色视频综合| e午夜精品久久久久久久| 精品国产乱子伦一区二区三区| 亚洲五月色婷婷综合| 久久人人精品亚洲av| 老汉色av国产亚洲站长工具| 亚洲国产精品成人综合色| 国产极品粉嫩免费观看在线| 日日夜夜操网爽| 亚洲中文字幕日韩| 日韩欧美一区二区三区在线观看| 亚洲狠狠婷婷综合久久图片| 免费在线观看成人毛片| 久久人妻福利社区极品人妻图片| 麻豆国产av国片精品| 久久热在线av| 精品国产超薄肉色丝袜足j| 久久亚洲精品不卡| 色播亚洲综合网| 成年免费大片在线观看| 中文字幕精品免费在线观看视频| 最好的美女福利视频网| 欧美日韩亚洲综合一区二区三区_| 在线观看日韩欧美| 欧美最黄视频在线播放免费| 国产成人欧美| 亚洲精品美女久久av网站| 一个人观看的视频www高清免费观看 | 最近最新免费中文字幕在线| www.精华液| 老司机深夜福利视频在线观看| 久久欧美精品欧美久久欧美| 亚洲性夜色夜夜综合| 日韩av在线大香蕉| 中出人妻视频一区二区| 成人国产综合亚洲| 香蕉av资源在线| 国产精品久久视频播放| 国产av一区在线观看免费| 白带黄色成豆腐渣| 两个人免费观看高清视频| 日韩av在线大香蕉| 国产在线精品亚洲第一网站| 老熟妇仑乱视频hdxx| 国产欧美日韩精品亚洲av| 禁无遮挡网站| 欧美久久黑人一区二区| 婷婷精品国产亚洲av| 一a级毛片在线观看| 国产成人精品久久二区二区91| 欧美成人免费av一区二区三区| 久久精品亚洲精品国产色婷小说| 一区二区日韩欧美中文字幕| 一本一本综合久久| 真人做人爱边吃奶动态| 亚洲 国产 在线| 在线天堂中文资源库| a在线观看视频网站| 伦理电影免费视频| 亚洲成人免费电影在线观看| 黄色视频不卡| 成人国语在线视频| 变态另类成人亚洲欧美熟女| 国产主播在线观看一区二区| 亚洲色图 男人天堂 中文字幕| 精品一区二区三区av网在线观看| 嫩草影视91久久| 国产精品九九99| 18禁黄网站禁片午夜丰满| 热re99久久国产66热| 成人国产一区最新在线观看| 18禁美女被吸乳视频| 老司机福利观看| 久久精品国产综合久久久| 神马国产精品三级电影在线观看 | 欧美日韩一级在线毛片| 中文字幕av电影在线播放| 自线自在国产av| 极品教师在线免费播放| 国产精品久久电影中文字幕| 日本黄色视频三级网站网址| 色综合站精品国产| 99久久国产精品久久久| 欧美日本视频| 美女免费视频网站| 男人的好看免费观看在线视频 | 色哟哟哟哟哟哟| 在线天堂中文资源库| 非洲黑人性xxxx精品又粗又长| 热re99久久国产66热| 人人妻,人人澡人人爽秒播| 国产成年人精品一区二区| 极品教师在线免费播放| 99热这里只有精品一区 | 亚洲精品久久国产高清桃花| 国产精品野战在线观看| 免费看a级黄色片| а√天堂www在线а√下载| 国产亚洲欧美精品永久| 国产午夜福利久久久久久| АⅤ资源中文在线天堂| 首页视频小说图片口味搜索| 亚洲国产精品合色在线| 亚洲久久久国产精品| 操出白浆在线播放| 国产精品亚洲美女久久久| 99re在线观看精品视频| 99riav亚洲国产免费| 日韩欧美在线二视频| 亚洲va日本ⅴa欧美va伊人久久| 日韩 欧美 亚洲 中文字幕| 欧美成人一区二区免费高清观看 | www.熟女人妻精品国产| 欧美+亚洲+日韩+国产| 99在线视频只有这里精品首页| 满18在线观看网站| 国产成人欧美在线观看| 欧美日韩精品网址| 日韩三级视频一区二区三区| 少妇的丰满在线观看| 亚洲一区中文字幕在线| 亚洲精品美女久久久久99蜜臀| 欧美日韩福利视频一区二区| 老司机在亚洲福利影院| 黄色视频不卡| 午夜久久久在线观看| 亚洲三区欧美一区| 听说在线观看完整版免费高清| 国产精品亚洲av一区麻豆| 成人国产综合亚洲| 久久精品影院6| xxxwww97欧美| 日韩中文字幕欧美一区二区| 嫩草影视91久久| 国产伦人伦偷精品视频| 国产精品一区二区免费欧美| 亚洲真实伦在线观看| 国产精品久久久久久人妻精品电影| 国产人伦9x9x在线观看| 久久人人精品亚洲av| av超薄肉色丝袜交足视频| 伦理电影免费视频| 欧美日本亚洲视频在线播放| 欧美成人一区二区免费高清观看 | aaaaa片日本免费| 日本熟妇午夜| 美女高潮到喷水免费观看| 十分钟在线观看高清视频www| 婷婷精品国产亚洲av| 岛国视频午夜一区免费看| x7x7x7水蜜桃| 一边摸一边做爽爽视频免费| 日韩欧美免费精品| xxx96com| 精品国产一区二区三区四区第35| 可以在线观看的亚洲视频| 老司机深夜福利视频在线观看| 他把我摸到了高潮在线观看| 久99久视频精品免费| 一本大道久久a久久精品| 日韩欧美一区二区三区在线观看| 人人妻人人澡人人看| 亚洲最大成人中文| 国产黄a三级三级三级人| 成年女人毛片免费观看观看9| 欧美成人午夜精品| 欧美激情高清一区二区三区| 亚洲欧美日韩高清在线视频| 美女午夜性视频免费| 91国产中文字幕| 国产亚洲精品av在线| 午夜a级毛片| 麻豆国产av国片精品| 欧美乱色亚洲激情| 欧美黑人欧美精品刺激| 日韩欧美在线二视频| 黄色成人免费大全| 男女下面进入的视频免费午夜 | 国产熟女午夜一区二区三区| 国产99久久九九免费精品| 一二三四在线观看免费中文在| 免费av毛片视频| 桃红色精品国产亚洲av| 久久伊人香网站| 国产亚洲精品综合一区在线观看 | 老熟妇乱子伦视频在线观看| 熟女电影av网| 精品久久久久久,| 亚洲狠狠婷婷综合久久图片| 欧美色欧美亚洲另类二区| 精品日产1卡2卡| 看片在线看免费视频| 午夜视频精品福利| 亚洲成人国产一区在线观看| 欧美成人性av电影在线观看| 一进一出抽搐动态| 啦啦啦观看免费观看视频高清| 久久香蕉精品热| 两个人免费观看高清视频| 51午夜福利影视在线观看| 免费在线观看亚洲国产| 丁香六月欧美| 国产精品综合久久久久久久免费| 在线天堂中文资源库| 欧美日韩福利视频一区二区| 女同久久另类99精品国产91| 两个人视频免费观看高清| 怎么达到女性高潮| 在线永久观看黄色视频| 午夜老司机福利片| 欧美日韩黄片免| 黄频高清免费视频| 国产精品一区二区免费欧美| 久久精品aⅴ一区二区三区四区| 黄色片一级片一级黄色片| 777久久人妻少妇嫩草av网站| 亚洲精品粉嫩美女一区| 最新美女视频免费是黄的| 亚洲自拍偷在线| 国产精华一区二区三区| 欧美精品亚洲一区二区| 国产黄色小视频在线观看| 国产av一区在线观看免费| 男人舔奶头视频| 十八禁人妻一区二区| 国产亚洲精品一区二区www| 91麻豆精品激情在线观看国产| 亚洲欧美激情综合另类| 18禁观看日本| 中亚洲国语对白在线视频| 色综合婷婷激情| 免费av毛片视频| e午夜精品久久久久久久| 亚洲精品在线美女| 中国美女看黄片| 欧美乱码精品一区二区三区| 国产一区二区激情短视频| 在线播放国产精品三级| 日韩成人在线观看一区二区三区| 很黄的视频免费| 成人国产综合亚洲| 男女之事视频高清在线观看| 国产亚洲精品av在线| 中文字幕人妻丝袜一区二区| 人妻久久中文字幕网| 神马国产精品三级电影在线观看 | 一级毛片女人18水好多| 久久久久久久久久黄片| 国产片内射在线| 久久人人精品亚洲av| 99久久无色码亚洲精品果冻| 亚洲精品色激情综合| 亚洲第一欧美日韩一区二区三区| 亚洲精品在线观看二区| 亚洲成人久久爱视频| 午夜两性在线视频| 99久久无色码亚洲精品果冻| 757午夜福利合集在线观看| 久久久久久大精品| 国产精品电影一区二区三区| 欧美激情 高清一区二区三区| 成熟少妇高潮喷水视频| 美女免费视频网站| 亚洲欧美激情综合另类| 日韩欧美国产一区二区入口| 免费在线观看黄色视频的| 欧美精品亚洲一区二区| 国产精品久久久久久人妻精品电影| 中文资源天堂在线| a级毛片a级免费在线| 久久久久久久久免费视频了| 久久精品aⅴ一区二区三区四区| 最新在线观看一区二区三区| 男人舔奶头视频| 麻豆国产av国片精品| 女同久久另类99精品国产91| 亚洲人成77777在线视频| 免费在线观看黄色视频的| 亚洲国产高清在线一区二区三 | 日本在线视频免费播放| 久久亚洲精品不卡| 国产精品久久电影中文字幕| 一级a爱片免费观看的视频| x7x7x7水蜜桃| 午夜福利在线在线| 18禁黄网站禁片午夜丰满| 亚洲 欧美 日韩 在线 免费| 久久久久国内视频| 午夜激情av网站| 精品国产乱码久久久久久男人| 亚洲 欧美一区二区三区| 国产一卡二卡三卡精品| 国产精品 欧美亚洲| 色在线成人网| x7x7x7水蜜桃| 精品人妻1区二区| 欧美色视频一区免费| 麻豆成人av在线观看| 久久精品国产亚洲av高清一级| 可以在线观看毛片的网站| 久久99热这里只有精品18| 精品欧美一区二区三区在线| 搡老岳熟女国产| 桃色一区二区三区在线观看| 国产爱豆传媒在线观看 | 成人亚洲精品av一区二区| 日本五十路高清| 国产极品粉嫩免费观看在线| 12—13女人毛片做爰片一| 一区二区三区精品91| 一级a爱视频在线免费观看| av视频在线观看入口| 国产伦在线观看视频一区| videosex国产| 亚洲午夜理论影院| 女生性感内裤真人,穿戴方法视频| 精品一区二区三区av网在线观看| 天堂动漫精品| 91大片在线观看| 色综合亚洲欧美另类图片| 亚洲国产中文字幕在线视频| 亚洲片人在线观看| 国产久久久一区二区三区| 在线观看一区二区三区| 久久伊人香网站| 久久狼人影院| 狂野欧美激情性xxxx| 91成人精品电影| 久久这里只有精品19| 亚洲全国av大片| 国产欧美日韩一区二区精品| 老司机靠b影院| 欧美乱妇无乱码| 国产亚洲精品久久久久久毛片| 一个人观看的视频www高清免费观看 | 久久天躁狠狠躁夜夜2o2o| 手机成人av网站| 久久精品国产亚洲av高清一级| 国产亚洲精品一区二区www| 久久久精品欧美日韩精品| 男男h啪啪无遮挡| 波多野结衣高清作品| xxx96com| 久久精品影院6| 看免费av毛片| 视频区欧美日本亚洲| 99精品在免费线老司机午夜| av有码第一页| 看黄色毛片网站| 久久香蕉国产精品| 麻豆一二三区av精品| 啦啦啦 在线观看视频| 午夜免费鲁丝| 黄色视频不卡| 精品国产国语对白av| 久久久久久久午夜电影| 麻豆一二三区av精品| 色综合站精品国产| 好看av亚洲va欧美ⅴa在| 欧美日韩亚洲综合一区二区三区_| 欧美乱码精品一区二区三区| 亚洲欧美一区二区三区黑人| 国产精品亚洲美女久久久| 中文字幕另类日韩欧美亚洲嫩草| 超碰成人久久| 国产精品 国内视频| 免费无遮挡裸体视频| 精品免费久久久久久久清纯| 不卡av一区二区三区| 51午夜福利影视在线观看| www日本在线高清视频| 亚洲aⅴ乱码一区二区在线播放 | 免费电影在线观看免费观看| 中文字幕久久专区| a级毛片a级免费在线| netflix在线观看网站| 亚洲人成77777在线视频| 91九色精品人成在线观看| 亚洲九九香蕉| 在线观看舔阴道视频| 久久中文字幕一级| 欧美不卡视频在线免费观看 | 亚洲精品国产一区二区精华液| 久久精品aⅴ一区二区三区四区| www.熟女人妻精品国产| 男女床上黄色一级片免费看| 波多野结衣巨乳人妻| 听说在线观看完整版免费高清| 国产伦在线观看视频一区| 一夜夜www| 成人国产综合亚洲| 高清毛片免费观看视频网站| 丝袜在线中文字幕| 老司机午夜十八禁免费视频| 啦啦啦免费观看视频1| 国产成人啪精品午夜网站| 在线永久观看黄色视频| 啦啦啦观看免费观看视频高清| av有码第一页| 国产免费男女视频| 亚洲国产欧美一区二区综合| 精品一区二区三区视频在线观看免费| 亚洲人成电影免费在线| 欧美日本亚洲视频在线播放| 国产精品久久电影中文字幕| 99精品欧美一区二区三区四区| 国产一级毛片七仙女欲春2 | 变态另类成人亚洲欧美熟女| 啦啦啦 在线观看视频| 免费在线观看影片大全网站| 18美女黄网站色大片免费观看| 50天的宝宝边吃奶边哭怎么回事| 国产精品久久久人人做人人爽| 亚洲欧美精品综合久久99| www.自偷自拍.com| 中国美女看黄片| 亚洲午夜精品一区,二区,三区| 听说在线观看完整版免费高清| 成人精品一区二区免费| 国产av不卡久久| 欧美一级a爱片免费观看看 | 亚洲激情在线av| 欧美中文综合在线视频| 97碰自拍视频| 99精品久久久久人妻精品| 久热爱精品视频在线9| 人成视频在线观看免费观看| 欧美乱色亚洲激情| 美女高潮到喷水免费观看| svipshipincom国产片| 国产成人系列免费观看| 在线看三级毛片| 亚洲国产欧美日韩在线播放| 精品久久久久久久末码| 午夜激情av网站| cao死你这个sao货| 一本大道久久a久久精品| 免费高清在线观看日韩| 亚洲av中文字字幕乱码综合 | 日本一本二区三区精品| 日韩国内少妇激情av| 亚洲人成伊人成综合网2020| 欧美日韩黄片免| 国产精品 国内视频| 香蕉久久夜色| 色综合站精品国产| 国产aⅴ精品一区二区三区波| 久久狼人影院| 久久午夜亚洲精品久久| 日韩精品中文字幕看吧| 国产三级黄色录像| 精品国内亚洲2022精品成人| 免费高清视频大片| 丁香六月欧美| 91麻豆精品激情在线观看国产| 久久午夜综合久久蜜桃| 亚洲成人国产一区在线观看| 久久久国产成人精品二区| 亚洲专区中文字幕在线| 嫁个100分男人电影在线观看| 国产亚洲av高清不卡| 18美女黄网站色大片免费观看| 一级黄色大片毛片| 亚洲精品中文字幕一二三四区| 亚洲最大成人中文| 91av网站免费观看| 亚洲国产日韩欧美精品在线观看 | 麻豆成人午夜福利视频| 每晚都被弄得嗷嗷叫到高潮| 1024手机看黄色片| 一进一出抽搐gif免费好疼| 一区二区三区精品91| 亚洲精品国产区一区二| 99在线视频只有这里精品首页| 欧美日韩黄片免| 欧美日韩福利视频一区二区| 亚洲全国av大片| 免费在线观看黄色视频的| 精品国产亚洲在线| 欧美黑人巨大hd| 欧美国产精品va在线观看不卡| 精品久久久久久久毛片微露脸| 欧美最黄视频在线播放免费| 哪里可以看免费的av片| 可以在线观看毛片的网站| 久久精品影院6| 日本一本二区三区精品| 国产精品免费一区二区三区在线| www.精华液| 久久精品国产亚洲av香蕉五月| 操出白浆在线播放| 久久久久久九九精品二区国产 | 俺也久久电影网| 亚洲 欧美 日韩 在线 免费| 嫁个100分男人电影在线观看| 成年版毛片免费区| 99久久无色码亚洲精品果冻| 国产精品99久久99久久久不卡| 久久精品91无色码中文字幕| 亚洲色图av天堂| 欧美黑人巨大hd| 一级a爱视频在线免费观看| cao死你这个sao货| 亚洲色图av天堂| 夜夜爽天天搞| 欧美日韩精品网址| 久久久久国内视频| 精品少妇一区二区三区视频日本电影| 在线观看午夜福利视频| 亚洲熟女毛片儿| 一区二区三区精品91| 老司机靠b影院| 三级毛片av免费| 夜夜看夜夜爽夜夜摸| 久久精品国产亚洲av高清一级| 精品久久久久久久久久免费视频| 精品国产一区二区三区四区第35| 99久久综合精品五月天人人| 精品人妻1区二区| aaaaa片日本免费| www.自偷自拍.com| 亚洲成人久久爱视频| 麻豆久久精品国产亚洲av| 日韩大码丰满熟妇| 欧美成人午夜精品| 给我免费播放毛片高清在线观看| 成人特级黄色片久久久久久久| 母亲3免费完整高清在线观看| 久久精品亚洲精品国产色婷小说| 精品久久久久久久毛片微露脸| 脱女人内裤的视频| 亚洲精品色激情综合| 成人18禁高潮啪啪吃奶动态图| 日本 av在线| 国产又黄又爽又无遮挡在线| 免费高清视频大片| 成人精品一区二区免费| 亚洲av成人av| 久久 成人 亚洲| 黑丝袜美女国产一区| www日本在线高清视频| 视频在线观看一区二区三区| 国产亚洲精品久久久久久毛片| av福利片在线| av中文乱码字幕在线| 久久精品国产99精品国产亚洲性色| 久久人人精品亚洲av| 老司机福利观看| 1024香蕉在线观看| 日韩大尺度精品在线看网址| 国产亚洲欧美在线一区二区| 精品国产国语对白av| 免费女性裸体啪啪无遮挡网站| 女人爽到高潮嗷嗷叫在线视频| 午夜亚洲福利在线播放| 99国产精品一区二区三区| 香蕉丝袜av| 搞女人的毛片| 精品国产乱码久久久久久男人| 老司机深夜福利视频在线观看| 亚洲av第一区精品v没综合| 国产精品国产高清国产av| 老鸭窝网址在线观看| 成人国产综合亚洲| 一边摸一边抽搐一进一小说| 动漫黄色视频在线观看| 好男人在线观看高清免费视频 | 国产黄色小视频在线观看| 欧美精品啪啪一区二区三区| 一个人观看的视频www高清免费观看 | 精品少妇一区二区三区视频日本电影| 桃红色精品国产亚洲av| 国产午夜福利久久久久久| 欧美日韩福利视频一区二区| 亚洲全国av大片| 亚洲免费av在线视频| 美女大奶头视频| 日本在线视频免费播放| 十八禁网站免费在线| 看黄色毛片网站| 亚洲精品中文字幕一二三四区| 精品不卡国产一区二区三区| 中国美女看黄片| 欧美亚洲日本最大视频资源| 欧美日韩福利视频一区二区| 婷婷精品国产亚洲av在线| 国产视频内射| 最好的美女福利视频网| 亚洲精品国产一区二区精华液| 少妇被粗大的猛进出69影院| 黄色a级毛片大全视频| 免费高清视频大片| 看片在线看免费视频| 香蕉丝袜av| 国产乱人伦免费视频| 男女视频在线观看网站免费 | 亚洲一区中文字幕在线| 母亲3免费完整高清在线观看| 免费无遮挡裸体视频| 亚洲五月天丁香| 中文亚洲av片在线观看爽| 桃红色精品国产亚洲av| 欧美亚洲日本最大视频资源| 国产精品美女特级片免费视频播放器 |