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

    Data?driven disturbance observer?based control: an active disturbance rejection approach

    2021-05-19 04:05:42HarveyRojasCubidesJohnCortRomeroJaimeArcosLegarda
    Control Theory and Technology 2021年1期

    Harvey Rojas?Cubides · John Cortés?Romero · Jaime Arcos?Legarda

    Abstract

    Keywords Active disturbance rejection control · Data-driven disturbance observer · Disturbance observer-based control ·Signal prediction

    1 Introduction

    External disturbances, measurement noise, parameter uncertainty and unmodeled dynamics are the most common issues in the design of robust control strategies. For decades, each of these aspects has been handled independently. However,starting in the 90s, based on the work of J. Han, a unified approach to deal with all these adverse effects, which lumps internal and external disturbances, including the unmodeled dynamics, into a concept called total disturbance was developed [1, 2]. This concept joined to the mechanisms for its estimation and compensation is the core of the active disturbance rejection control (ADRC) approach [3, 4].Besides, total disturbance idea is closely related to lumped disturbance or equivalent input disturbance in the context of disturbance/uncertainty estimation and attenuation techniques [5].

    Under this approach, external disturbances, uncertainties and unmodeled dynamics are considered as a lumped disturbance, which is attenuated with a disturbance rejectionbased control. This technique is widely used to improve the robustness of feedback controllers against external disturbances and model uncertainties. The methods to perform disturbance rejection control can be grouped, according to their design approach, into three main categories: (1) timedomain, (2) frequency domain, and (3) data-driven.

    The first method, the time domain-based design, uses the state-space representation of the system model to assign additional states that represent the internal dynamics of the disturbance signal. In this way, a state observer is designed to estimate all the states of the model, including the states that represent the disturbance signal and its internal model.There are several observer strategies that are used this approach: the extended state observer (ESO) [6, 7], active disturbance rejection control (ADRC) [8], generalized extended state observer (GESO), also known as generalized proportional integral (GPI) observer [9, 10], unknown input observers (UIO) [11], and high-gain observer (HGO) [12].

    The second method uses the frequency domain analysis. The main structure of this method is the disturbance observer (DOB), which usually uses a plant inverse and a Q-filter design that ensures a causal realization of the observer [13]. The main drawbacks of this method is the use of the plant inverse, which might be an unstable inverse if the system is non-minimum phase. This problem has been addressed for several researchers that have proposed the use of stable inverse [14]. To mitigate the design efforts of the plant inverse and the Q-filter necessary into the conventional DOB’s structure, Zheng et al. [15] proposed a generalized design method that minimizes the norm of the dynamics from disturbance to its estimation error. The digital implementation of the DOB entails the Q-filter implementation, which is generally designed as a digital versions of basic analog filters [16, 17], zero-phase low-pass FIR filter[18, 19], discrete binomial filters [20] and infinite impulse response filters [21-23], among others.

    The third method uses data-driven techniques. These methods do not necessarily consider the explicit system models or the disturbance internal model. Instead, it uses the input-output data to predict the trend of the disturbance signal [24]. This paper explores the third method with a novel data-driven disturbance observer. Despite of the huge number of scientific contributions on disturbance observers, there are very few reports exploring the opportunities that the data analysis and signal processing could produce into the area of disturbance estimation. This issue has been approached using the disturbance history data to solve a least square regression problem that fits the disturbance signal history to a polynomial basis function, which, in the sense of Taylor’s series, is the most effective and simplest methods for curve fitting [25-27]. Moreover, the time series expansion and the Moore-Penrose pseudoinverse has been used to estimate higher order disturbances in nonlinear systems [28].Besides the curve fitting, data information has been used in learning strategies to enhance the DOB’s performance under quasi-periodic disturbances [15]. In this way, techniques, such as iterative learning control [29, 30] have been developed to attenuate the effect of periodic disturbances;however, the learning process hardly guarantees stability in the transient responses.

    This paper proposes a control based on a data-driven disturbance observer (3DOB), which guarantees stability in the transient response through an inner control loop. Additionally, the approach used by the 3DOB to estimates the disturbance does not require explicit design of the plant inverse or a Q-filter. Thus, this approach eliminates the structural constraints of the conventional DOB. The remaining of this paper is organized as follows: Sect. 2 describes the general outline of the proposed approach. Section 3 presents the data-driven disturbance observer structure. Section 4 shows the proposed predictor scheme that is the essential part of the data-driven disturbance estimation and rejection. Section 5 presents design criteria and stability-performance analysis.Section 6 shows the experimental validation of the control strategy on a switched DC-DC synchronous buck converter.Finally, Sect. 7 summarizes the accomplishments of this work and draws recommendations for future developments.

    2 Proposed approach

    A general outline of the proposed approach is presented in Fig. 1. Consider a minimum phase plant described by the following discrete time model that relates the measurable outputy(k) and the control inputu(k):

    wherenis the system order, andNG(z) andDG(z) are polynomials in thezdomain. In addition, systemG(z) is perturbed by an equivalent input (total) disturbanceξ(k) which includes all endogenous and exogenous adverse effects.Therefore, the goal is to design a control system that allows robust tracking of the referencey*(k).

    The proposal is defined under the general disturbance estimation and attenuation approach, and particularly in the context of the active disturbance rejection control framework. Likewise, the control strategy is completely presented in discrete time, avoiding the use of discretization methods for its digital implementation, e.g., a bilinear transformations.

    In general terms, the proposed strategy is divided into two stages: inner or stabilization loop and outer or filter-prediction loop. The purpose of the inner loop is to stabilize the system dynamics while adjusting its internal model. This allows, on the one hand, to design a complementary disturbance rejection scheme and, on the other, to adjust a robust tracking according to the internal model considered. A key point in the inner loop is the definition of a disturbance signal equivalent to the output(ξ1) that can be synthesized thorough knowledge of the available output and an auxiliary control input (u1). This equivalent disturbance signal will be processed as part of the synthesis of the outer loop controller.

    Fig. 1 General control scheme for the proposed approach

    The outer loop is inspired by the guidelines of a control structure based on a disturbance observer, where for this case,the inverse of the equivalent system resulting from the inner loop is required. As usual, this inverse corresponds to an improper transfer function causing implementation problems due to dependence on future values. Precisely, those related predictions associated with the equivalent disturbance can be estimated online based on a polynomial approximation of the trend associated with the data-signal history in the Taylor sense.

    2.1 Stabilization using the inner loop controller

    Consider the plantG(z) described in (1) and a unified equivalent disturbance to the control inputξ(t) such that The degree and structure of the controllerC1(z) is such that it allows, on the one hand, the arbitrary location of the roots of the polynomialD1(z)=DC1(z)DG(z)+NY(z)NG(z) and, on the other hand, the rejection of structured disturbance signals whose internal models have been considered in theC1(z)design. As initial requirement, the numerator resulting from the controller,NC1(z) , must have its roots within the unit circle of the complex plane. In the event thatNC1(z) exhibits unstable roots, an additional procedure must be performed.This restriction can be easily manipulated from the algebraic point of view through a two-degree of freedom control structure [31]. For this purpose, the following subsection proposes an adaptation of the two-parameter compensator.

    At this point, it is important to highlight that the system obtained from the internal loop is stable and it also has the ability to track or reject the signals that are included in the internal model considered inC1(z).

    2.2 Two?parameter compensator issue

    As previously indicated, a requirement in the inner loop is that theNC1(z) polynomial has stable roots for a stable model-inverse when formulating the control law. This problem can be solved by including the classic two-parameter compensator presented in Fig. 2.

    The proposed approach follows, in the first stage of the design, the model matching methodology to impose a stable numerator in the resulting control system.

    Under the current methodology, a two input-one output dynamics is considered as in Eq. (6):

    Hence, the intent of using the two-parameter compensator is to match the transfer function of the control system to a desirable numerator and denominator model, both with stable roots by construction. The first is to implement a stable inversion and the second is to achieve arbitrary pole placement.

    A supplementary polynomialNU1(z) is imposed through the following control law:

    Fig. 2 Two-parameter controller structure

    Remark 1In the following,K1(z) can be used instead ofC1(z) when a two-parameter controller procedure was not necessary, and the design is based on the inner loop diagram presented in Fig. 1.

    Once a suitable stabilization has been carried out by means of the internal loop, an outer loop based on a datadriven disturbance observer is proposed.

    3 Data?driven disturbance observer

    To find a simplified equivalent disturbance signal in the outer loop, let us defineΞ1(z)=Ξ(z)∕K1(z) as a filtered version of the equivalent disturbanceΞ(z) . This is,

    Then a useful characteristic definition of the new filtered disturbance, justifying the application of the stabilization loop, can be obtained online by a simpley(k),u1(k) processing, i.e.,

    Hence, an ideal and non-implementable control law forU1(z)is given by

    whereT-1(z) is an improper transfer function that demands future values ofY(z) causing implementation problems. This general math definition, bears similarity with the idea used in the classical disturbance observer design (in an inner loop); however, in this work, a different way to deal with this problem is presented as described below.

    Definingras the relative degree ofT(z) and to obtain an implementable version of (18), a delayed version of the disturbanceΞ1(z) is proposed as

    Nevertheless, up to this point only ar-delayed version (not the current one) of the disturbanceξ1(k-r) is available.For this reason, a data-driven predictorP(z) that estimates the current version of the disturbance(k) is included as follows:

    The proposed control scheme is shown in Fig. 3. This scheme is composed of an inner loop together with a data-driven disturbance observer in the outer loop. The inner loop includes a Two-parameter compensatorK(z)=[K1(z) -K2(z)]designed to obtain a desiredT(z), based on the nominal/known model of the plantG(z) and assuming a perfect disturbance rejection. On the other hand, the disturbance observer in the outer loop uses the reference transfer functionT(z) and itsr-delayed inversez-rT-1(z) (winch is a proper transfer function) to obtain ar-delayed disturbanceξ1(k-r) . Finally, a predictor schemeP(z) is included to achieve an estimation of the current disturbance(k).

    Fig. 3 Proposed control scheme

    Fig. 4 Prediction process

    4 Predictor scheme

    An excess in the number of samplespover the order of the polynomial regressionmis required to obtain a low-pass filter effect. The criteria for this adjustment involve the knowledge of noise level of the signals, sampling time, and bandwidth.

    According to the previous analysis, it is possible to define the family of transfer functions for the predictorP(z) given by

    5 Design criteria and stability?performance analysis

    To analyze the effect of the data-driven disturbance observer over the system performance and stability, the control law(22) is replaced in the output dynamics (15) as follows:

    whereQ(z)=z-rP(z) . From this expression, it can be inferred thatTQn=Q(z) , andSQn=1-Q(z) represent the complementary sensitive and sensitive functions of the estimation process, respectively. Now, using the predictor transfer function (26), we have

    5.1 Closed?loop stability

    With the aim to explore the closed-loop stability analysis,equation (29) is rewritten as

    which allows to draw the following observations about closed-loop stability.

    On the one hand, the inner loop process design guarantees the stability of the transfer functionT(z). On the other hand,given the structure of the obtained predictor, the equivalent function 1-Q(z) places its poles at zero, which also assure its inherent stability. Therefore, (35) presents a linear dynamic system with two inputs one output, predominantly governed by two stable transfer functions.

    The inner loop, called the stabilizer loop, is in charge of the stable priming structure of the overall control system.As conceived in the proposed two-stage ADRC approach,the outer loop performs an estimation and consequent rejection of the EID while preserving the stable structure imposed by the inner loop. This property remains invariant only for adequate performance of the equivalent DOB. In fact, the proposed strategy can be viewed as robust feedback linearization scheme. The validity of the resulting linearization is directly proportional to the effectiveness of the rejection of the unified disturbance. This can be expressed simply (frequency domain) since the DOB must comply with a minimum bandwidth. Once this requirement is fulfilled, the linearization of the system is achieved and one could even speak of separability.

    It is worth to note that according the ADRC framework,every mismatch of the model based of the design with respect to the real behaviour of the system can be included into the unified total disturbance. Therefore, missing nonlinear dynamics or related uncertainties are also considered as part of the lumped disturbance and can be handled by the proposed method.

    On the other hand, a deeper non linear robust stability analysis can be difficult issue and it is out of the scope of the present work, see, e.g., [32]. However, a simple stability analysis, based on the linearization process, like the one presented in this work could be favorable in an practical engineering context [33]. In this sense, even though in this paper has not given an exhaustive result related with the specific effect of nonlinear uncertainties, owing to mathematical complexity, the main features of the proposed approach and its efficacy has been illustrated using experimental results.

    5.2 Design parameters and performance

    According to the equation (32) and the block diagram presented in Fig. 3, it is possible to establish that transfer functionQ(z)=z-rP(z) fulfills a similar role as the Q-filter in the classical disturbance observer-based control (DOBC) [20],or an extended state observer in the ADRC framework [34].However, the proposed approach focuses on the design of a data-driven observer based on a predictor. Nonetheless,this work takes advantage of the similarities with DOBC,using some well-known frequency domain analysis tools.This facilitates the robustness-performance analysis and it makes the choice of the design parameters clearer.

    The tuning parameters of the data-driven disturbance observer are: the polynomial orderm, the number of samplespin the predictor, and the sampling period (in some cases). Thus, for each combination in the design parameters,the predictor is represented by a different transfer function(P(z)) and the frequency response of the observer changes.This feature shows the flexibility of the proposed approach.To present the properties of this strategy, Figs. 5 and 6 show the frequency response ofQ(z) forr=1 and some combination ofmandp.

    Fig. 5 Frequency response of Q(z) varying the polynomial order for 5, 10 and 20 samples p, and h=20 μs

    An additional feature shown in Fig. 6 is the effect of number of samples over robustness and performance.Here, a low number of samples achieve better performance in low frequency, but the bandwidth is higher causing issues with the noise. While, a larger number of samples produces a low-pass filter effect; and however, this also produces a lower performance.

    5.3 Equivalence of Q(z) in the S?domain

    It is possible to establish an equivalent in the S-domain for the final results obtained in this work. For instance,using a bilinear transformation ofQ(z), where for a given sampling periodhwe have

    As an example, Fig. 7 presents the frequency response ofQ(z) and its equivalent in the s-domainQ(s) (based on the bilinear transformation) for three combinations of design parametersmandp, forr=1 andh=20 μs.

    The results show that it is conceivable that under a certain design and with an appropriate tuning, a Q-filter could be obtained with a similar performance as presented in this work, at least in the zone of low and medium frequencies.However, it is clear that this continuous-time design would not correspond to a data-driven approach, which is the core of this work.

    6 Experimental validation

    With the aim to validate the proposed control strategy, the case of a switched DC-DC synchronous buck converter (SBC) is considered. This kind of converter produces a regulated output voltage that is lower than its input voltage, and can providehigh currents while improving the efficiency. As shown in Fig. 8, a generic SBC is comprised mainly of an input voltage sourcevi, two power MOSFETs controlled with PWM (with their gate drivers), an LC output filter, and the load resistance.The simplified average dynamic model of a buck converter for voltage-mode control is described as

    Fig. 6 Frequency response of Q(z) varying the number of samples for polynomial order m=[0,1,2] , and h=20 μs

    whereRis the load resistance of the converter,Lis a filter inductance,Cis an output capacitor,vi(t) is the input voltage,v0(t) is the average output capacitor voltage, and the duty cycleu(t)∈[0,1] represents the control signal.

    Some elements that motivate the use of this platform to test the control proposed approach are (1) synchronous buck DC-DC converter is a hybrid system, given the switching mode of the circuit. Moreover, there are several commonly unmodeled nonlinear effects, such as magnetic characteristics of the inductor, electro-magnetic interference caused by switching and dead time, among others, which have influence in the control performance [35]. (2) The system components values can undergo changes during its lifetime,which cause a relevant uncertainty factor for output voltage regulation [36]. (3) During its normal operation, the system is affected both by matched disturbances at the control channel (input voltage changes), as well as, by mismatched disturbances caused by current load changes. This fact, is often a challenge for the control strategies under the ADRC and DOB approaches [37]. (4) The voltage-mode for buck converters is characterized both by its simplicity due it no need to sense the inductor current, and also its low openloop output impedance. However, it has some disadvantages due its poor line rejection and because control-to-output transfer function changes between operating modes; therefore, a disturbance estimation and rejection scheme could be useful [38]. (5) Taking into account its real applications,such as motor drivers, chargers for electric vehicles and DC transmission, among others, it is desirable to achieve a direct digital controller with an easier and cheaper implementation in standard electronic devices, e.g., a microcontroller [39].

    Considering the above, the main challenge is to find a control law such that the output voltagev0(t) is forced to track the given referencev*o(t) , even in presence of disturbances, such as changes in the load, input voltage variations,uncertainties in electronic components and non modeled dynamics, among others.

    Fig. 7 Frequency response of Q(z) and its equivalent in the s-domain varying the design parameters, for r=1 and h=20 μs

    Fig. 8 Circuit diagram of DC-DC synchronous buck converter

    Table 1 Parameter of the DC-DC synchronous buck converter

    6.1 Experimental setup

    Fig. 9 Experimental setup including (1) DC-DC synchronous buck converter with active load, (2) 32-bit floating-point microcontroller,(3) Programmable DC power supply, (4)-(6) system to generate the input voltage disturbances, (7) and (8) measuring instruments

    The experimental test setup configuration and the prototype are depicted in Fig. 9. In general terms, the platform includes: (1) the power stage composed by a DC-DC synchronous buck converter with active load,(2) a 32-bit floating-point microcontroller unit (MCUTMS320F28377S) to implement the control algorithms using Code Composer Studio v10.1 and Matlab v2020a,(3) a programmable desktop laboratory DC power supply(HM305) with the measurement of voltage, current and power in the input side of the converter; a system to generate the input voltage disturbances that includes: (4) an arbitrary waveform generator (UTG932E), (5) a generic power supply and (6) electronic circuits for conditioning;and some measuring instruments, such as (7) digital oscilloscope (DS1045Z) and (8) digital multimeter (UT70B).On the hole, the MCU controls the synchronous buck operation by driving the high-side and low-side switches of the half bridge using on-chip PWM signals with a frequency of 200 kHz and a dead time of 10 ns. For his part, a resistive voltage divider together with a passive first order low pass filter are used to measure the output signal. Moreover, an analog-to-digital converter with a resolution of 12 bits and a sample time for the controller updating is set up to 20 μ s. Finally, the nominal parameters of the power converter are listed in Table 1.

    6.2 Controller design

    Based into the general design procedure presented in Sects.2 and 3, the nominal transfer functionsG(z),T(z), andz-rT(z)-1(withr=1 ) are

    Table 2 Predictor design parameters

    Regarding to the predictor design parameters, three combinations of polynomial order (m) and number of samples (p),were selected. The design criteria for the controllers used in the experimental tests are presented in Table 2.

    6.3 Performance?robustness test

    In this section, the robustness and performance of the proposed approach were tested on the DC-DC synchronous buck converter affected by disturbances in the input voltage and the load.

    6.3.1 Case 1: Robustness against load changes

    In this first scenario, the load resistance is reduced during the operating process. To achieve this, we used a load permanently connected in the converter output, (Rn) , in addition to an active load to increase the output current up to 4.75 times the nominal value. It provides an active load feature at run time and allows transient performance tests.

    The experimental results of the output voltage with the proposed strategy using the three different design criteria and the classical-DOB are shown in Fig. 10. To compare the performances obtained by each controller, some metrics index,such as IAE, ISE, MAE and MSE are used, as proposed in [40]. Finally, the quantitative comparison using the performance indices are shown in Table 4.

    Table 3 Parameters used for DOB and HODOB, where k is the number of extended states

    Fig. 10 Experimental results against load changes

    Table 4 Closed-loop performance criteria against load changes

    6.3.2 Case 2: Robustness against input voltage fluctuations

    The second scenario test the robustness of the proposed strategy against input voltage variations. This test has special relevance in buck converters due its poor line rejection and because the control-to-output transfer function changes between operating modes. Likewise, the input voltage disturbance acts as a matched perturbation input in the input-output system model. During the test, the input voltage has smooth changes around its nominal value (Vin), with an amplitude of 22% and a frequency of 10 Hz, according to the following profile:

    The experimental curves of the output voltage under the proposed 3DOB-based controller and the classical DOB are shown in Fig. 11. Moreover, the performances indices IAE,ISE, MAE and MSE are presented in Table 5.

    According to the above experimental tests, although both 3DOB and classical DOB approaches are able to reject the mismatched disturbances caused by hard sudden load changes, the closed-loop performance is improved by the proposed 3DOB-based control as is shown in Table 4. Furthermore, the performance results in 3DOB further improves with a large polynomial order since the local approximation of the disturbance is better, even for smaller bandwidth.On the other hand, results presented in Fig. 11 and Table 4 show that the results of the two approaches are comparable; however, the proposed control strategy minimizes all criteria compared with the classical strategy in the case of input voltage fluctuations. Besides, significant improvements occur in 3DOB with polynomials with order greater than zero.

    Finally, it is worth noting that according with practical implementation issues, such as noise and sampling period,it is possible to find different tuning rules for both the proposed data-driven disturbance observer and classical DOB and HODOB, which allow to achieve satisfactory performances. However, from the results obtained in this work, it is viable to conclude that the data-driven proposed approach is capable of achieving similar and even superior performances (under certain conditions) compared to another disturbance/uncertainty estimation and attenuation techniques,such as disturbance observer-based control (DOBC). Likewise, it is important to note that, the inherent phase lag isreduced in the disturbance observer mechanism, leading to potential for faster response, which is needed in applications like DC-DC converter.

    Fig. 11 Experimental results against input voltage fluctuations

    Table 5 Closed-loop performance criteria against input voltage fluctuations

    7 Conclusions

    A data-driven disturbance observer-based control strategy with an active disturbance rejection approach is developed.The general control scheme is composed by two stages: first a inner stabilisation loop based on an unity feedback controller or a two parameters compensator is developed to set the desired closed-loop dynamics. Then to estimate and reject the equivalent disturbance in the outer loop, an observer that includes a prediction scheme, based on a least-square regression, is proposed.

    A key feature of the proposed approach is its flexibility to select the design parameters, such as the polynomial orderm, and the number of samplespin the predictor. This feature allows to adjust the best performance and at the same time to reduce the effect of noise. Furthermore, the proposed predictor scheme fulfills the same effects of the classical-DOB, but without the explicit use of an inverse of the plant dynamics or the design of a digital Q-filter. Additionally, it is possible to use the frequency response analysis to help the parameter tuning and to achieve the desired performance.

    In other matters, the proposed approach not only has a performance and a robustness comparable at least with those of another related methods, but also it presents a new point of view that takes advantage of the use of data. This aspect that could be considered as a superiority and an opportunity to open a line of research regarding the use of data-driven techniques in disturbance estimation and rejection schemes.

    Finally, the robustness and performance of the proposed strategy is validated in a case study for the voltage control of a DC-DC synchronous buck converter affected by disturbances in the input voltage and the load. The experimental results show the data-driven disturbance observer-based control is able to reject matched and mismatched disturbances and the performance is improved in comparison to the classical disturbance observer-based control.

    av电影中文网址| 国产精品乱码一区二三区的特点| 久久久水蜜桃国产精品网| 久久婷婷成人综合色麻豆| 精品国产乱子伦一区二区三区| 三级毛片av免费| 亚洲九九香蕉| 亚洲片人在线观看| 国产一级毛片七仙女欲春2 | 欧美不卡视频在线免费观看 | 欧美另类亚洲清纯唯美| 久久性视频一级片| 亚洲狠狠婷婷综合久久图片| 婷婷六月久久综合丁香| 午夜激情av网站| 波多野结衣巨乳人妻| 亚洲精品国产一区二区精华液| 在线观看66精品国产| 亚洲国产欧美网| 成人亚洲精品av一区二区| 欧美另类亚洲清纯唯美| 在线av久久热| 别揉我奶头~嗯~啊~动态视频| 久久天躁狠狠躁夜夜2o2o| 亚洲精品美女久久av网站| 变态另类丝袜制服| 久久精品aⅴ一区二区三区四区| 日韩成人在线观看一区二区三区| 村上凉子中文字幕在线| av中文乱码字幕在线| 久久精品人妻少妇| 精品一区二区三区视频在线观看免费| 一二三四在线观看免费中文在| 白带黄色成豆腐渣| 久久精品国产综合久久久| 国产精品野战在线观看| 国产人伦9x9x在线观看| 精品欧美一区二区三区在线| 亚洲色图 男人天堂 中文字幕| 亚洲国产精品成人综合色| 精品国产乱子伦一区二区三区| 首页视频小说图片口味搜索| 最近在线观看免费完整版| 俄罗斯特黄特色一大片| 黄片小视频在线播放| 1024香蕉在线观看| 免费电影在线观看免费观看| 宅男免费午夜| 国产精品,欧美在线| 草草在线视频免费看| 男女做爰动态图高潮gif福利片| 欧美黑人欧美精品刺激| 免费看a级黄色片| 日本黄色视频三级网站网址| 色在线成人网| 一a级毛片在线观看| 国产精品久久久久久亚洲av鲁大| 亚洲专区中文字幕在线| av中文乱码字幕在线| 成人特级黄色片久久久久久久| 国产乱人伦免费视频| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美激情久久久久久爽电影| 久久国产精品人妻蜜桃| 成人18禁在线播放| xxxwww97欧美| 久久久久久亚洲精品国产蜜桃av| 丝袜人妻中文字幕| 色播亚洲综合网| 免费女性裸体啪啪无遮挡网站| 波多野结衣高清作品| 欧美午夜高清在线| 国产精品一区二区三区四区久久 | 成人亚洲精品av一区二区| 一级毛片女人18水好多| 免费在线观看视频国产中文字幕亚洲| 天天一区二区日本电影三级| 免费一级毛片在线播放高清视频| 黄色丝袜av网址大全| 成年女人毛片免费观看观看9| 精品久久蜜臀av无| www.www免费av| 99精品欧美一区二区三区四区| 精品国产一区二区三区四区第35| 波多野结衣高清无吗| 美女免费视频网站| 久久中文看片网| 97碰自拍视频| 国产一级毛片七仙女欲春2 | 亚洲无线在线观看| 国产成年人精品一区二区| 日本 av在线| 免费看a级黄色片| 免费在线观看影片大全网站| 国产亚洲av高清不卡| 国产免费男女视频| 波多野结衣巨乳人妻| 91大片在线观看| x7x7x7水蜜桃| www.www免费av| 国内毛片毛片毛片毛片毛片| 国产精品,欧美在线| 国产又色又爽无遮挡免费看| av中文乱码字幕在线| 亚洲全国av大片| 亚洲中文字幕一区二区三区有码在线看 | 国产精品1区2区在线观看.| 99久久99久久久精品蜜桃| 国产成人av激情在线播放| 国产成人av激情在线播放| 中文字幕最新亚洲高清| 法律面前人人平等表现在哪些方面| 狂野欧美激情性xxxx| 亚洲免费av在线视频| 精品一区二区三区四区五区乱码| 久久久国产成人免费| 久久午夜亚洲精品久久| 高清毛片免费观看视频网站| 久久天躁狠狠躁夜夜2o2o| 99久久无色码亚洲精品果冻| 欧美zozozo另类| 国产成人精品无人区| 亚洲精品中文字幕一二三四区| 国产精品av久久久久免费| 露出奶头的视频| 一个人免费在线观看的高清视频| 国产一区二区激情短视频| 亚洲第一青青草原| 日韩欧美三级三区| 亚洲成人精品中文字幕电影| 欧美av亚洲av综合av国产av| 波多野结衣高清无吗| 日本五十路高清| av中文乱码字幕在线| 亚洲国产毛片av蜜桃av| 亚洲欧美日韩无卡精品| 巨乳人妻的诱惑在线观看| 在线观看一区二区三区| 国产真人三级小视频在线观看| 日韩欧美一区二区三区在线观看| 久久中文看片网| 岛国视频午夜一区免费看| 欧美黑人欧美精品刺激| 在线天堂中文资源库| 淫秽高清视频在线观看| 亚洲精品中文字幕一二三四区| 成人午夜高清在线视频 | 91在线观看av| 一本大道久久a久久精品| 曰老女人黄片| 久热爱精品视频在线9| 又紧又爽又黄一区二区| 精品不卡国产一区二区三区| 丝袜在线中文字幕| 国产午夜福利久久久久久| 亚洲精品在线美女| 老汉色∧v一级毛片| 每晚都被弄得嗷嗷叫到高潮| 无遮挡黄片免费观看| 丰满的人妻完整版| 国产欧美日韩精品亚洲av| 国内毛片毛片毛片毛片毛片| 51午夜福利影视在线观看| a在线观看视频网站| 男人操女人黄网站| 99精品久久久久人妻精品| 欧美丝袜亚洲另类 | 国内精品久久久久久久电影| 成人亚洲精品av一区二区| 亚洲 欧美一区二区三区| 亚洲专区字幕在线| 亚洲一区二区三区色噜噜| 波多野结衣av一区二区av| 亚洲久久久国产精品| 亚洲av熟女| 黄片大片在线免费观看| 精品无人区乱码1区二区| √禁漫天堂资源中文www| 欧美乱码精品一区二区三区| 国产高清视频在线播放一区| 国产成人系列免费观看| АⅤ资源中文在线天堂| a在线观看视频网站| 欧美黄色淫秽网站| 夜夜躁狠狠躁天天躁| 亚洲色图av天堂| 2021天堂中文幕一二区在线观 | 欧美成人免费av一区二区三区| 可以在线观看毛片的网站| av免费在线观看网站| 午夜福利高清视频| 中文字幕人成人乱码亚洲影| 丰满的人妻完整版| 亚洲狠狠婷婷综合久久图片| 黄片大片在线免费观看| 性欧美人与动物交配| 国产精品永久免费网站| 两个人看的免费小视频| 久久精品91蜜桃| 99国产极品粉嫩在线观看| 亚洲片人在线观看| 国产精品九九99| 一级片免费观看大全| 午夜激情福利司机影院| 成人18禁高潮啪啪吃奶动态图| 免费在线观看视频国产中文字幕亚洲| 99久久无色码亚洲精品果冻| 午夜激情av网站| 少妇裸体淫交视频免费看高清 | 日韩有码中文字幕| 国产免费av片在线观看野外av| 亚洲性夜色夜夜综合| 亚洲婷婷狠狠爱综合网| 一区二区三区免费毛片| 美女xxoo啪啪120秒动态图| 精品欧美国产一区二区三| 男人狂女人下面高潮的视频| 在线播放国产精品三级| 日韩国内少妇激情av| 女同久久另类99精品国产91| 我的老师免费观看完整版| 亚洲欧美成人综合另类久久久 | 欧美最黄视频在线播放免费| 悠悠久久av| 国产私拍福利视频在线观看| 狂野欧美激情性xxxx在线观看| 国产成人影院久久av| h日本视频在线播放| 悠悠久久av| 久久精品影院6| 久久久久九九精品影院| 欧美丝袜亚洲另类| 99久国产av精品国产电影| 一卡2卡三卡四卡精品乱码亚洲| 久久鲁丝午夜福利片| 精品久久久久久久久亚洲| 黑人高潮一二区| 一本一本综合久久| 欧美国产日韩亚洲一区| 人妻制服诱惑在线中文字幕| 欧美+日韩+精品| 亚洲电影在线观看av| 不卡一级毛片| 少妇丰满av| 国产一区二区三区在线臀色熟女| 国产成人a∨麻豆精品| 亚洲一区高清亚洲精品| 亚洲av中文字字幕乱码综合| 国产成人精品久久久久久| 男人舔女人下体高潮全视频| 欧美日韩乱码在线| 麻豆精品久久久久久蜜桃| 又黄又爽又刺激的免费视频.| 天天一区二区日本电影三级| 三级国产精品欧美在线观看| 精品国产三级普通话版| 亚洲熟妇熟女久久| 一本久久中文字幕| 国产视频内射| 国内精品久久久久精免费| 夜夜看夜夜爽夜夜摸| 亚洲欧美成人综合另类久久久 | 伦精品一区二区三区| 九九在线视频观看精品| 无遮挡黄片免费观看| 欧美性感艳星| 国产麻豆成人av免费视频| 亚洲av成人精品一区久久| 一个人看的www免费观看视频| 国产一区二区在线观看日韩| 在线观看美女被高潮喷水网站| 1024手机看黄色片| 久久久精品大字幕| 国产亚洲精品久久久久久毛片| 婷婷色综合大香蕉| 美女xxoo啪啪120秒动态图| 少妇熟女欧美另类| 在线国产一区二区在线| 欧美+亚洲+日韩+国产| 国产一区二区激情短视频| 久久天躁狠狠躁夜夜2o2o| 国产av一区在线观看免费| 老熟妇乱子伦视频在线观看| 亚洲美女黄片视频| 国产成人福利小说| 乱系列少妇在线播放| 十八禁国产超污无遮挡网站| 国产精品国产高清国产av| 国产精品av视频在线免费观看| 丝袜美腿在线中文| 亚洲久久久久久中文字幕| 日韩精品中文字幕看吧| 少妇丰满av| 一区二区三区免费毛片| 一本精品99久久精品77| 国产一级毛片七仙女欲春2| 少妇的逼好多水| 国产精品美女特级片免费视频播放器| 午夜视频国产福利| 国产乱人视频| 中文字幕熟女人妻在线| 久久久国产成人精品二区| 少妇熟女aⅴ在线视频| 18禁黄网站禁片免费观看直播| 日韩欧美免费精品| 如何舔出高潮| 伊人久久精品亚洲午夜| 18禁黄网站禁片免费观看直播| 精品久久久久久久末码| 天堂影院成人在线观看| 夜夜看夜夜爽夜夜摸| 97在线视频观看| 国产精品一二三区在线看| 亚洲,欧美,日韩| 亚洲国产精品sss在线观看| 欧美zozozo另类| 国产乱人偷精品视频| 亚洲专区国产一区二区| 精品午夜福利在线看| 国产美女午夜福利| 国产午夜精品论理片| 日本色播在线视频| 一个人观看的视频www高清免费观看| 久久久久久久久久久丰满| 熟女电影av网| 免费高清视频大片| 日本一二三区视频观看| 亚洲在线观看片| 国产乱人视频| 日韩成人av中文字幕在线观看 | 欧美日韩一区二区视频在线观看视频在线 | 草草在线视频免费看| 欧美激情在线99| 真实男女啪啪啪动态图| 午夜福利在线观看免费完整高清在 | 又粗又爽又猛毛片免费看| 中国国产av一级| 蜜桃亚洲精品一区二区三区| 欧美区成人在线视频| 欧美性猛交黑人性爽| 免费看日本二区| 欧美区成人在线视频| 国产一区二区激情短视频| 国产91av在线免费观看| 中文亚洲av片在线观看爽| 免费高清视频大片| 精品午夜福利视频在线观看一区| 欧美+日韩+精品| 在线播放无遮挡| 少妇人妻一区二区三区视频| 嫩草影院入口| 成人一区二区视频在线观看| 欧美日本亚洲视频在线播放| 一边摸一边抽搐一进一小说| 看片在线看免费视频| 国产真实伦视频高清在线观看| 亚洲精品久久国产高清桃花| 美女被艹到高潮喷水动态| 亚洲成人精品中文字幕电影| 少妇猛男粗大的猛烈进出视频 | 高清毛片免费看| 国产精品久久久久久亚洲av鲁大| 久久人人爽人人片av| 国产三级中文精品| 国产aⅴ精品一区二区三区波| 成年免费大片在线观看| 亚洲性久久影院| 一级a爱片免费观看的视频| 插阴视频在线观看视频| 天天一区二区日本电影三级| 日本精品一区二区三区蜜桃| 99久国产av精品| 亚洲国产高清在线一区二区三| 91久久精品国产一区二区三区| 久久精品国产亚洲av天美| 国产女主播在线喷水免费视频网站 | 免费人成视频x8x8入口观看| 亚洲图色成人| 91久久精品电影网| 色播亚洲综合网| 联通29元200g的流量卡| 精华霜和精华液先用哪个| 美女cb高潮喷水在线观看| 亚洲人成网站在线观看播放| 久久精品夜夜夜夜夜久久蜜豆| 久久国内精品自在自线图片| 特级一级黄色大片| 亚洲av.av天堂| 欧美3d第一页| 久久久久久九九精品二区国产| 日本a在线网址| 人妻丰满熟妇av一区二区三区| 久久精品国产自在天天线| 久久久国产成人免费| 网址你懂的国产日韩在线| 狂野欧美激情性xxxx在线观看| 亚洲天堂国产精品一区在线| av天堂在线播放| 国产亚洲av嫩草精品影院| 欧美日韩国产亚洲二区| 一进一出好大好爽视频| 美女高潮的动态| 1000部很黄的大片| 亚洲av熟女| 联通29元200g的流量卡| 国产男靠女视频免费网站| 嫩草影院入口| 亚洲成av人片在线播放无| 日日撸夜夜添| 干丝袜人妻中文字幕| 国内精品美女久久久久久| 日韩av在线大香蕉| 联通29元200g的流量卡| 国产不卡一卡二| 可以在线观看的亚洲视频| 三级经典国产精品| 久久精品综合一区二区三区| 国产av一区在线观看免费| 色在线成人网| 亚洲精品久久国产高清桃花| 特级一级黄色大片| 哪里可以看免费的av片| 国产久久久一区二区三区| 国产三级中文精品| 欧美成人a在线观看| 天堂av国产一区二区熟女人妻| 搡女人真爽免费视频火全软件 | 啦啦啦观看免费观看视频高清| 3wmmmm亚洲av在线观看| 国产美女午夜福利| 成人亚洲欧美一区二区av| 少妇裸体淫交视频免费看高清| 啦啦啦韩国在线观看视频| www.色视频.com| 麻豆av噜噜一区二区三区| 国内精品美女久久久久久| 69av精品久久久久久| 免费观看精品视频网站| 国产蜜桃级精品一区二区三区| 久久99热6这里只有精品| 亚洲一区二区三区色噜噜| 天天躁日日操中文字幕| 中文资源天堂在线| 十八禁网站免费在线| 蜜桃久久精品国产亚洲av| 2021天堂中文幕一二区在线观| 国产淫片久久久久久久久| 三级男女做爰猛烈吃奶摸视频| 日韩成人av中文字幕在线观看 | 少妇的逼好多水| 国产在线精品亚洲第一网站| 欧美bdsm另类| a级毛色黄片| 国产爱豆传媒在线观看| 插阴视频在线观看视频| 色综合色国产| 舔av片在线| 国产单亲对白刺激| av黄色大香蕉| 亚洲精品一区av在线观看| 国产淫片久久久久久久久| 老熟妇乱子伦视频在线观看| 久久综合国产亚洲精品| 99久久精品一区二区三区| 香蕉av资源在线| 蜜桃久久精品国产亚洲av| 无遮挡黄片免费观看| 韩国av在线不卡| 久久久久久九九精品二区国产| 午夜免费激情av| 国内精品一区二区在线观看| 18禁在线播放成人免费| 午夜视频国产福利| 91麻豆精品激情在线观看国产| 亚洲不卡免费看| 精华霜和精华液先用哪个| 久久精品国产亚洲网站| 国产爱豆传媒在线观看| 1024手机看黄色片| 在线天堂最新版资源| 天天躁日日操中文字幕| 天天一区二区日本电影三级| 欧美性感艳星| 看黄色毛片网站| 国产亚洲91精品色在线| 老师上课跳d突然被开到最大视频| 日韩精品有码人妻一区| 又爽又黄无遮挡网站| 最好的美女福利视频网| 日韩 亚洲 欧美在线| 身体一侧抽搐| 久久久色成人| 国内精品美女久久久久久| 99久久精品国产国产毛片| 一级黄色大片毛片| 一夜夜www| 久久久久久久久中文| 97在线视频观看| 国产精品野战在线观看| 十八禁网站免费在线| 亚洲在线自拍视频| 国产精品一及| 国产精品永久免费网站| 国产精品av视频在线免费观看| 欧美另类亚洲清纯唯美| 成人三级黄色视频| 啦啦啦啦在线视频资源| 国产三级中文精品| 一进一出抽搐动态| 国产高清视频在线播放一区| 久久午夜亚洲精品久久| 国产精品久久久久久久电影| 久久久久国内视频| 波野结衣二区三区在线| 一个人观看的视频www高清免费观看| 免费看光身美女| 嫩草影院新地址| 色5月婷婷丁香| 99九九线精品视频在线观看视频| 18+在线观看网站| or卡值多少钱| 91在线观看av| 乱人视频在线观看| 少妇的逼水好多| 久久久久久大精品| 看黄色毛片网站| 国产精品爽爽va在线观看网站| 18+在线观看网站| 淫秽高清视频在线观看| 免费人成在线观看视频色| 99热只有精品国产| 观看免费一级毛片| 国产蜜桃级精品一区二区三区| 亚洲第一电影网av| 99热网站在线观看| 日本一本二区三区精品| 99热只有精品国产| 久久人妻av系列| 国产精品久久久久久精品电影| 尤物成人国产欧美一区二区三区| 欧美日本亚洲视频在线播放| 男人狂女人下面高潮的视频| av国产免费在线观看| 亚洲久久久久久中文字幕| 久久精品国产自在天天线| 精品久久久久久成人av| 日韩精品有码人妻一区| 欧美国产日韩亚洲一区| 色综合亚洲欧美另类图片| 乱码一卡2卡4卡精品| 日本精品一区二区三区蜜桃| 国产探花在线观看一区二区| 国国产精品蜜臀av免费| 97超碰精品成人国产| 欧美性猛交黑人性爽| 精品久久久久久久人妻蜜臀av| 少妇的逼好多水| 日韩一区二区视频免费看| 久久久久久国产a免费观看| 真人做人爱边吃奶动态| 亚洲国产欧洲综合997久久,| 中文字幕久久专区| 国产精品日韩av在线免费观看| 真人做人爱边吃奶动态| 欧美性猛交╳xxx乱大交人| av在线亚洲专区| 亚洲av第一区精品v没综合| 国产精品一区二区三区四区免费观看 | 亚洲四区av| 免费在线观看影片大全网站| 人妻少妇偷人精品九色| 12—13女人毛片做爰片一| 三级毛片av免费| 久久久久久九九精品二区国产| 日本撒尿小便嘘嘘汇集6| 精品人妻视频免费看| 国产91av在线免费观看| 香蕉av资源在线| 日韩大尺度精品在线看网址| 九九热线精品视视频播放| av在线亚洲专区| 久久韩国三级中文字幕| 插逼视频在线观看| 男插女下体视频免费在线播放| 国产成人精品久久久久久| 最近2019中文字幕mv第一页| 欧美日本亚洲视频在线播放| 欧美成人免费av一区二区三区| 久久久久久伊人网av| 大又大粗又爽又黄少妇毛片口| 国产一区二区亚洲精品在线观看| 中文字幕av在线有码专区| 免费看美女性在线毛片视频| 老熟妇乱子伦视频在线观看| 日本一二三区视频观看| 日本免费一区二区三区高清不卡| 国产视频内射| avwww免费| 国内精品久久久久精免费| 精品久久国产蜜桃| 国产伦一二天堂av在线观看| 丝袜美腿在线中文| 成年女人看的毛片在线观看| 午夜激情福利司机影院| 精华霜和精华液先用哪个| 国产乱人偷精品视频| 亚洲av成人av| 久久久久久久久久成人| 一区二区三区免费毛片| 两个人的视频大全免费| 啦啦啦啦在线视频资源| 日本 av在线| 亚洲av成人精品一区久久| 成人特级av手机在线观看| 搞女人的毛片|