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

    A phase?locked loop using ESO?based loop filter for grid?connected converter: performance analysis

    2021-05-19 05:48:16BaolingGuoSeddikBachaMazenAlamirJulienPouget
    Control Theory and Technology 2021年1期

    Baoling Guo · Seddik Bacha · Mazen Alamir · Julien Pouget

    Abstract

    Keywords Grid-connected converters · Grid disturbances · Phase-locked loop · Loop filter · Extended state observer ·Tuning approach · Stability analysis · Phase margin · Noises attenuation · Power hardware-in-the-loop (PHIL) test

    1 Introduction

    Grid-connected converters (GcCs) have various applications such as grid-integrated renewable energy generation systems[1], active power filters [2], uninterruptible power supplies[3], etc. The phase-locked loop (PLL) is widely used for the GcCs to obtain the information of the voltage phase and the grid frequency [4, 5]. A PLL is a closed-loop feedback control system, which consists of three basic elements: phase detector, loop filter (LF), and voltage-controlled oscillator[5]. The LF is used to suppress disturbances inside the PLL control loop. The tracking characteristics, response dynamics, and stability properties of the PLL are also mainly determined by the LF [5]. The PLLs commonly use a PI controller as the LF [4-7]. However, the PLL control systems involve various disturbances such as phase jumps, frequency deviations, voltage magnitude variations, and nonlinear dynamics[6-10]. The PLL systems may loose its synchronization after a sudden large disturbance [11, 12].

    Large efforts are devoted to enhancing the steady-state performance and the quality of the transient of conventional PLLs. The author of [4] reviews various types of LFs including Type-1, quasi-Type-1, quasi-Type-2, and Type-3, which are applied for particular purposes. Moreover, the grid voltage measurements involve switching noises and ripples,and sensor noises [13, 14]. Under such circumstances, the expected results of PLL would be not reliable. Some works attempt to use the prefilter before its input or add the lowpass filler inside the control loop to address such issues [13,14]. Furthermore, some adaptive LF-based PLLs are proposed such as the adaptive frequency estimation loop design[15] and the adaptive loop gain mechanism [16]. However,the strategy proposed in [15] has a high requirement on the PLL control input. The adaptive loop gain may cause oscillations problems [16].

    The extended state observer (ESO) is a partial-modelbased observer, which is first defined in the frame of active disturbance rejection control (ADRC) [17]. Although it is originally proposed in a nonlinear structure [17], a linear ESO has been increasingly used as it is easier to design and simpler to implement [18]. The ESO-based control treats all disturbed factors affecting the plant such as system nonlinearities, uncertainties, and external disturbances as a generalized disturbance [17]. The generalized disturbance is then observed by ESO and compensated into the closedloop dynamics. Such structure enables excellent capabilities in dealing with disturbances or uncertainties [19]. Moreover, the previous work indicates that the ESO-based control achieves better ability to attenuate high-frequency noises[20]. This work attempts to exploit the ESO’s potential advantages to achieve an efficient and robust PLL system,and therefore, improve control performances of the disturbed GcCs.

    The previous works commonly compare the ESO-based controller with PI controllers to verify its effectiveness[21-23]. The critical point is to determine under which conditions, the comparison might be viewed as fair. The conventional PI is a single-degree-of-freedom controller[24]. The proportional term and integral term are selected by considering the stability conditions and the quality of the transient [7]. An ESO-based controller, however, is a two-degree-of-freedom (2DOF) controller [25]. The controller bandwidth and the observer bandwidth can be adjusted separately. The observer bandwidth is mainly limited by sensor sensitivity [26]. In [8], a linear ESO is used to enhance control dynamics of the PLL. The natural frequency of PI closed-loop control system and the control bandwidth of ESO-based control are tuned to be the same for the comparison study [8]. However, such comparison is not really fair[20]. The authors of [20] formulate the linear ADRC by a feedback compensator with a prefilter. The compensator can be generated with a PI controller filtered by a low-pass filter[20]. The tuning approach discussed in [20], which is based on the well-tuned parameters of PI controller, is considered in this work. The feedback compensator is tuned to have almost the same gain at cross frequency with the PI control system, which results in a relatively fair comparison between these two different types of controllers.

    This work can be viewed as an extension of our previous work [8], which particularly deals with the following three aspects:

    1. First, this work introduces a different tuning approach.The parameters of ESO-based controller are generated based on the well-designed PI controller, which results in a fair comparison with conventional PI-type PLLs.

    2. Then, the frequency domain properties are quantitatively analysed with respect to the control stabilities and the noises rejection.

    3. Finally, more complete Matlab/Simulink simulation and power hardware-in-the-loop (PHIL)-based experimental results are provided to analyse the performances of the ESO-based PLLs and check its applicability to the GcCs.

    The paper is organized as follows: the ESO-based PLLs are briefly reviewed in Sect. 2. The tuning approach and the parameter selection guidelines are presented in Sect. 3. The frequency domain analysis is followed in Sect. 4. Simulation studies are conducted before the experimental validation in Sect. 5. The PHIL tests and experimental results are presented in Sect. 6. Conclusions are provided in Sect. 7.

    2 ESO?based SRF?PLL design

    The proposed ESO-based LF can be incorporated into SRF-PLLs equipped with additional filters for abnormal grid conditions, whose filters are placed before its input [27] or inside the control loop [28]. Two representative designs under unbalanced grid, decoupled PLL (DPLL) [27] and moving average filter (MAF)-PLL [28], as illustrated in Figs. 2 and 3, are particularly studied for the performance analysis. The design details of DPLL/MAF-PLL can be found in previous works [27, 28].

    Fig. 2 Schematic diagram of ESO-based DPLL, with low-pass filter (LPF), vα , vβ the voltages in a αβ rotating reference frame, Tdq+ and Tdq- are the positive and negative sequences transforming matrices, respectively, superscript +∕- represent the positive and negative sequences components, respectively, and the upper index ‘ - ’indicates variables after filters

    Fig. 3 Schematic diagram of ESO-based MAF-PLL, where the upper index ‘ - ’ indicates variables after filters

    2.1 SRF?PLL modelling

    The modelling of SRF-PLL is briefly reviewed (see Fig. 1).The grid voltages are generally represented by

    2.2 Extended state observer

    A key procedure is to reformulate a practical control issue to the canonical model of cascade integrators with total disturbances [29]. Based on expression (4) and Fig. 4, the manipulation is conducted as follows:

    A proper disturbance definition is critical to enable a highefficiency ESO [29]. Since, in the PLL, thevdeis the realtime estimation ofVm, the control gainb=Vm∕vdeapproximates 1. Therefore, the estimation of control gain is set tob0=1 . Besides, the nominal grid frequencyωgis feed-forwarded in the control input. Such operation helps improving the ESO estimation’s efficiency [29]. Finally, the total disturbancedis formulated by

    2.3 Control law formulation

    Based on Fig. 4, the control law is represented by

    Fig. 4 Linearised control diagram of ESO-based SRF-PLL,β1 and β2 are the positive gains of ESO, δw is the unknown input disturbances, and b0 is the estimation of control gain

    3 Tuning approach

    The previous research [20] has proved that, under some conditions, a first-order ESO-based controller can stabilize the plant if it can be stabilized by a PI controller. Therefore,we can configure an ESO-based controller based on a welltuned PI controller. More importantly, such design enables a relatively fair comparison between the PI and the ESO-based controller.

    3.1 Formulation of 2DOF block diagram

    First, the transfer functions can be obtained from (5)-(10)as follows:

    the ESO-based controller is formulated to a 2DOF closedloop system as shown in Fig. 5a. Now, we have a similar structure as a classical PI control system Fig. 5b except the prefilter. The authors of [32] indicate that the prefilter does not affect the feedback properties.

    Fig. 5 a ESO-based control 2DOF diagram. b Classical PI closedloop control diagram. δw : the unknown input disturbances. δn : measurement noises

    3.2 Updating ESO?based controller with PI controller

    where the gain correction factorNis introduced to ensure that the ESO-based control has the similar behaviour with PI control in low frequency. By solving Eq. (18), we have

    Note that the denominator (KPωo-2KI) in Eq. (19) must be positive to satisfy the stability condition, leading to the following expression:

    Once the parametersξandωoare determined, the gain correction factorNcan be obtained with Eq. (18). As shown in Fig. 6, the transfer functions (14) and (16) present the same bode behaviours by introducing the gain correction factor. The gain correction factorNis introduced to ensure that the ESO-based control has the similar behaviour with PI control in low frequency. Otherwise, there would exist a constant gain difference between them. The parameters tuning guidelines will be detailed as follows.

    Fig. 6 Bode diagram of C(s), CPI(s) , and

    3.3 Tuning guidelines

    Now, we can generate the ESO-based controller with the PI controller. Based on the above discussions, the tuning guidelines are summarized as follow:

    · First, we choose the well-tuned PI parameters that have been validated in previous works.

    · Regarding the parameterb0, since, in the PLL, thevdeis the real-time estimation ofVm, the control gainb=Vm∕vdeapproximates 1. It is reasonable to set the estimation of control gain tob0=1.

    · The coefficientξ=2 is commonly chosen for the Hurwitz stable consideration [18]. Starting fromωo >ξKI∕KP, we increaseωofor Eqs. (16)-(19) until reaching noise limitations [26].

    · Then, the parametersNandωcare the solutions of (18)-(19), respectively.

    · Note that, for a PLL with in-loop filter, the coefficientξhas to be corrected to ensure the control stability.

    · Finally, it is mandatory to validate the applied parameters of ESO-based controller in the physical tests.

    The parameters of the chosen PI controllers for PLLs are given in Table 1. The parametersKPandKIfor SRF-PLL are selected based on the design method described in [7]. Besides,the well-designed parameters for the DPLL and the MAF-PLL provided in [27, 28] are considered, as well.

    Table 1 Well-tuned PI parameters for different types of PLLs

    4 Frequency domain analysis

    The stability index is a crucial factor for evaluating the quality of a controller [10, 20, 33]. An ESO-based PLL is described in [8]; however, no detailed stabilities analyses are provided. The stability properties are, therefore, quantitatively analysed in this work. Besides, the high-frequency noise rejection ability and the closed-loop performance are analysed, as well. The noise suppression effect of ESO itself is also presented. Note that the tuning guidelines discussed in Sect. 3 will be followed in the forthcoming analysis.

    4.1 Stabilities analysis

    The stability of the prefilterH(s) (15) is guaranteed, as all coefficients of its denominator are positive [33]. Indeed, the prefilter does not affect the feedback properties [32]. Therefore, the closed-loop stability can be determined by the forward gain transfer function, which is given by

    where the plant model isGp(s)=b∕s, andb=Vm∕vde. The stability will be analysed with respect to the observer bandwidthωo, the robustness for plant model uncertainness, and the effects of extra filters.

    Fig. 7 Bode plots of L(s) for ωo ∈{3ωn,5ωn,7ωn}

    The bode diagram ofL(s)=C(s)Gp(s) is shown in Fig. 7,withωo∈{3ωn,5ωn,7ωn} rad/s. The remaining parameters are selected by following the tuning guidelines summarized in Sect. 3. Figure 7 shows that, in low-frequency range, the ESO-based control behaves closely to the bode feature of the PI controller. Compared to the PI controller, however, its stability margin decreases due to the effect of low-pass filter.The phase margin (PM) within a range of 30?~60?is commonly recommended regarding the control stability [28].Therefore, the phase margin remains acceptable ( PM>50?).It can also be observed that the PM increases with largerωo,which can therefore enhance the control stabilities. However,theωocannot be too large, because it is limited by sensor noises and practical considerations. The parameterωomust be validated in practical conditions.

    The stability of the control system can be affected by the plant model discrepancies. The main uncertain parameter for the SRF-PLL is the plant gainbas indicated by Eq. (22).The control robustness are therefore assessed under different values ofb:

    The control gainb=Vm∕vdeapproximates 1, and its estimation is set tob0=1 . The bode plots of the loop transfer functionLb(s) forb∈{0.5,1.0,2.0,N,3.0} are provided in Fig. 8. The phase margin decreases when theNb0is far from the real value ofb. Though the phase margins are always acceptable.

    Now, the stabilities of ESO-based PLLs with extra fillers are discussed. Regarding the DPLL, the prefilter has no direct effect on the control stability [27, 32]. However, the stability of the ESO-based MAF-PLL could be affected,because a Moving Average Filter (MAF) is placed inside the control loop (see Fig. 3). The parameters have to be corrected; otherwise, the oscillation problems would occur in the practical tests. The loop gain transfer function is modified as follows:

    Fig. 8 Bode plots of Lb(s) for b ∈{0.5,1.0,2.0,N,3.0}

    with the parameterTw=0.01 s [28].

    The bode plots ofLM(s) for differentωoandξare, respectively, shown in Figs. 9 and 10. We use the following settings:ωo∈{3ωn,5ωn,7ωn} rad/s,ξ=2 for Fig. 9 andωo=5ωnrad/s,ξ∈{2,3,4} for Fig. 10. It can be seen that the PM and the gain margin (GM) decrease much due to the effect of the MAF. The PM and GM increase with a largerωo. However, theωois limited by sensor noises. The PM and GM can be enhanced using largerξ. While too largeξenables slow dynamics response [18]. Finally, the coefficientξ=4 ,ωo=5ωnis chosen for its application into the MAF-PLL.

    Fig. 9 Bode plots of LM(s) for ωo ∈{3ωn,5ωn,7ωn}

    Fig. 10 Bode plots of LM(s) for coefficient ξ ∈{2,3,4}

    4.2 Noises rejection analysis

    As illustrated in Fig. 5, the measurement noises rejection ability is determined by the transfer function as follows:

    wherey(s) andδn(s) are the Laplace transforms of the outputyand the measurement noisesδn, respectively.

    Figure 11 shows the bode diagram ofGn(s) . The comparative results indicate that ESO-based controller behaves better to attenuate the measurements noise compared to the PI controller. As indicted in (16), that is the low-pass filter which enhances its ability to reduce high-frequency noises. Besides,the noise rejection ability decreases when the parameterωoincreases.

    Now, we discuss the noise attenuation effect of ESO itself.Based on Eqs. (7)-(10), the transfer function of the estimated state variable1is obtained as follows:

    Then, theresulting from the measurement noiseδnand the input disturbanceδwcan be given by

    Fig. 11 Bode plots of Gn(s) for ωo ∈{3ωn,5ωn,7ωn}

    The bode plots offδn(s) andfδw(s) are shown in Figs. 12 and 13, respectively, withωo∈{3ωn,5ωn,7ωn} rad/s. It can be observed from Fig. 12 that the ESO performs as a low-pass filter, which can reject the high-frequency signals. Besides,the largerωocan increase the high-frequency gain. Namely,the ESO would be more sensitive to the high-frequency measurement noise when we increase theωo. It can be seen from Fig. 13 that the phase lag ofis decreased with largerωo; however, there is no influence on the high-frequency gain; the ESO can more efficiently attenuate the input disturbance in the low-frequency region. Consequently, there exists a trade-off between the input disturbance rejection and the measurement noise sensitivity [34].

    Fig. 12 Bode diagrams of the transfer function fδn(s) for ωo ∈{3ωn,5ωn,7ωn}

    Fig. 13 Bode diagrams of the transfer function fδw(s) for ωo ∈{3ωn,5ωn,7ωn}

    Based on (9), the transfer function of disturbances estimation error becomes

    The bode plots offe2(s) are shown in Fig. 14, withωo∈{3ωn,5ωn,7ωn} rad/s. The ESO can estimate the unknown disturbance very well only when the disturbances are composed of low-frequency components [10]. The disturbance estimation errors can be reduced by increasing theωo; however, theωois limited by the noises sensitivity.

    4.3 Closed?loop performance analysis

    Based on the formulated 2DOF diagram (see Fig. 5), the closed-loop transfer function of the ESO-based PLL is represented as follows:

    The bode plots ofGcl(s) are shown in Fig. 15. The closedloop performances of PI and ESO-based control system behave closely at the low-frequency range. It can also be noticed that there is a peak greater than 0 dB for the PI control system, which indeed corresponds a time-domain overshoot of the step response [20]. However, for the ESO-based control, the overshoot is attenuated due to the prefilter H(s).

    Fig. 14 Bode diagram of the transfer function fe2(s) for ωo ∈{3ωn,5ωn,7ωn}

    Fig. 15 Bode plot of closed-loop transfer function Gcl(s)

    5 Simulation studies and results

    Simulation studies are conducted with the Matlab/Simulink.We use the solver Ode3 (Bogacki-Shampine) with a fixed step size 0.1 ms [8]. The balanced three-phase voltages are set toVm=179 V,f=50 Hz. A white noise is manually added to simulate the voltage measurement noises. The ESO-based SRF-PLL/DPLL/MAF-PLL are, respectively,simulated by considering different disturbed regimes. The parameters of ESO-based PLLs are selected by following the tuning guidelines in Sect. 3; meanwhile, the control stability is taken into account, as well. The control parameters of ESO-based SRF-PLL and DPLL are set toωo=5ωnrad/s,ωn=157 rad/s,ξ=2 , andb0=1 . Regarding the MAF-PLL,the coefficientξ=4 is applied to ensure the control stability,and the remaining parameters are set to the same as those of the SRF-PLL.

    Fig. 16 Estimated frequencies of ESO-based SRF-PLL

    Fig. 17 Phase tracking errors of ESO-based SRF-PLL

    The unbalanced case is tested to assess the effects of the extra filters on the ESO-based LF. Figures 18 and 19 show the performances of ESO-based DPLL and MAFPLL. The voltage unbalances happen att=0.2 s, and the unbalanced voltages are set toVam=Vm,Vbm=0.8Vm,Vcm=0.8Vm. A phase jump of +20?is induced att=0.4 s.The grid frequency varies from 50 Hz to 52 Hz att=0.8 s.The simulation results for PI controllers are also provided for the comparisons. Some conclusions are suggested by the results, as shown in Figs. 18 and 19. First, the ESO-based DPLL and MAF-PLL present almost the same dynamic responses as the corresponding PI control for tracking the information of the frequency and the phase. Then, as the MAF itself could help attenuate the high-frequency noises,consequently, the advantage of ESO-based LF when applied to the MAF-PLL is not very significant. Moreover, the DPLL performs faster response dynamics compared to the MAF-PLL. Therefore, in the following part, only the ESObased DPLL is implemented into the closed-loop control of the GcC.

    Fig. 18 Estimated frequencies of ESO-based DPLL and MAF-PLL

    Fig. 19 Phase tracking errors of ESO-based DPLL and MAFPLL

    6 Experimental tests and results

    The closed-loop control performances of GcC equipped with ESO-based PLL are experimentally validated in this section.Experiments are conducted in a real-time hybrid PHIL test benchmark, as schematically illustrated in Fig. 20, which consists of the three-phase grid-connected converter, the disturbed grid emulator [35], and the dSPACE hardware (DS1005) [36].The three-phase GcC is a back-to-back Pulse Width Modulation (PWM) voltage source converter. The control signal is converted to the PWM command for the converters in real time[36]. The disturbed grid is emulated by the RT-LAB real-time digital simulator [8]. Key parameters of the testing benchmark are listed in Table 2.

    The control of GcCs mainly has three functions: the grid synchronization [8], the DC-link voltage control [37],and the grid current control [38]. First, the outer voltage loop maintains the constant DC-link voltage, where a PIcontroller is used in this work. Then, the proposed ESObased PLL is used to obtain the grid voltage phase to synchronize the injected current with the grid voltage. Furthermore, the injected current control is achieved through three Proportional Resonant (PR) controllers upon theabcthree-phase frame. More design details with respect to the control designs of GcCs can be found in the previous literatures [37-39]. Regarding the unbalanced conditions, various control strategies have been proposed before [40, 41].In our tests, we use the positive sequence control strategy proposed in [40]. The PLL is used to detect the phase angle that follows the positive sequence of the grid voltage. For the ESO-based PLLs, we keep the same control parameters as configured in the simulation studies. The well-tuned parameters of the PI controller for the DC-link voltage and the PR controllers for grid-injected currents are provided in Table 3, and the detailed definitions of these parameters can be found in [9]. Since there is no difference between the synchronization phase angle during the unbalances and that during normal grid conditions, the same control parameters values are thus used for both balanced and unbalanced cases.

    Fig. 20 PHIL experimental benchmark

    Table 2 Parameters of the experimental benchmark

    Figure 22 presents the closed-loop control performances of GcC equipped with ESO-based PLL/DPLL. Figure 22a shows the control performances under balanced grid. Testing results for unbalanced grid are provided in Fig. 22b. To highlight, we suppose that the unbalance is not affecting the phases but the magnitudes only. Then, the positive sequence is synchronized with the phase “a” of the original threephase voltages. Therefore, the estimated phase of the positive sequence follows the original phase angle.

    The test conditions are set as follows: a phase jump of+20?is manually added in Fig. 22a1 and b1, a frequency step of 2 Hz happens in Fig. 22a2 and b2, and an magnitude drop of 0.2Vmis considered in Fig. 22a3 and b3. First, theESO-based SRF-PLL/DPLL can correctly follow the phase,the frequency, and the voltage amplitude information. Then,the results indicate that the injected grid currentiasynchronizes with the grid voltageva. Furthermore, the GcC’s controls equipped with the ESO-based SRF-PLL/DPLL enable efficient performance regarding the control dynamics and the robustness in the presence of various disturbances.

    Table 3 Control parameters of grid-connected converter

    Fig. 22 Grid-connected control performances with ESO-based PLL/DPLL under balanced and unbalanced grid conditions, the grid current of phase a ia (2A/div), the grid voltage of phase a va (100 V/div), for a1 and b1, a phase jump of Δθ =+20π∕180 rad, the actual voltage phase θ(2rad/div), the estimated voltage phase (2rad/div); for a2 and b2, a frequency step of Δf =+2 Hz, the actual grid frequency f (1Hz/div), the estimated grid frequency (1Hz/div); for a3 and b3, a voltage magnitude variation of ΔVm =0.2Vm , the d-axis voltage vde (20 V/div), and the q-axis voltage vqe (20 V/div). (a1)-(a3) Balanced grid voltages. (b1)-(b3) Unbalanced grid voltages

    7 Conclusions

    This work is an extension for our previous work [8]. A linear ESO-based controller is designed as the LF for the SRFPLL. A tuning approach based on the well-tuned parameters of PI controller is provided, which results in a fair comparison with PI-type PLLs. Some conclusions can be drawn from the frequency domain analysis and the simulation results:the performances of the generated ESO-based controller are close to those of PI control at low frequency, while achieving a better ability of attenuating measurement noises. The phase margin decreases slightly, but remains acceptable. For the MAF-PLL with the in-loop filter, the equipped MAF itself could help attenuating the effects of high-frequency noises. Consequently, the advantage of ESO-based LF for its application to the MAF-PLL is not very significant. The PHIL experimental results suggest that the ESO-based PLLs can correctly obtain the phase, the frequency, and the voltage amplitude information. The obtained results show the effectiveness of ESO-based PLLs when applied to the GcC in the presence of various disturbances.

    It is worth mentioning that this work has discussed a tuning approach based on the well-designed PI control. However, for the implementations of ESO-based control, it is not mandatory to tune them based on a PI controller. The linear ESO-based controller can be tuned simply by following the parameters selection guidelines provided in the previous works [8, 18, 21].

    AcknowledgementsThis paper was supported by G2elab, Grenoble INP, University Grenoble Alpes, France and School of Engineering,HES-SO, Valais, Switzerland.

    FundingOpen Access funding provided by Haute Ecole Specialisée de Suisse occidentale (HES-SO).

    Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source,provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.

    免费人成在线观看视频色| 久久久久国产精品人妻一区二区| 亚洲人成网站在线观看播放| 青春草亚洲视频在线观看| 欧美精品一区二区免费开放| 一级,二级,三级黄色视频| 大话2 男鬼变身卡| 成年人免费黄色播放视频| 丁香六月天网| 中文字幕久久专区| 黄色欧美视频在线观看| 日韩欧美一区视频在线观看| 久久精品熟女亚洲av麻豆精品| 精品少妇黑人巨大在线播放| 国产高清不卡午夜福利| 永久免费av网站大全| 汤姆久久久久久久影院中文字幕| 久久精品久久久久久噜噜老黄| 亚洲精品乱久久久久久| 狂野欧美激情性xxxx在线观看| 一区二区三区四区激情视频| 大香蕉97超碰在线| 99热全是精品| av黄色大香蕉| 蜜桃在线观看..| 久久久久久久久久人人人人人人| 亚洲欧美一区二区三区黑人 | 999精品在线视频| 秋霞伦理黄片| 久久久午夜欧美精品| 一级a做视频免费观看| 亚洲av免费高清在线观看| 日日啪夜夜爽| 成年美女黄网站色视频大全免费 | 日韩一区二区视频免费看| 寂寞人妻少妇视频99o| 一本—道久久a久久精品蜜桃钙片| 99久久综合免费| 欧美日韩亚洲高清精品| 不卡视频在线观看欧美| 搡女人真爽免费视频火全软件| 啦啦啦中文免费视频观看日本| 啦啦啦啦在线视频资源| 欧美激情极品国产一区二区三区 | 欧美性感艳星| 欧美成人精品欧美一级黄| 日日摸夜夜添夜夜添av毛片| 亚洲精品日韩av片在线观看| 校园人妻丝袜中文字幕| 久久久久久人妻| 成年人免费黄色播放视频| 男人操女人黄网站| av国产精品久久久久影院| 99热全是精品| 欧美激情国产日韩精品一区| av免费观看日本| 欧美日韩在线观看h| 亚洲av不卡在线观看| av有码第一页| 亚洲少妇的诱惑av| 久久久久久久久久久丰满| 18禁观看日本| 亚洲精华国产精华液的使用体验| 国产极品粉嫩免费观看在线 | 制服丝袜香蕉在线| 亚洲欧美一区二区三区黑人 | 热99国产精品久久久久久7| 最近的中文字幕免费完整| 熟妇人妻不卡中文字幕| 人妻系列 视频| 在线观看www视频免费| 国产片特级美女逼逼视频| 2021少妇久久久久久久久久久| 亚洲情色 制服丝袜| 精品国产乱码久久久久久小说| av在线观看视频网站免费| 亚洲精品自拍成人| freevideosex欧美| 少妇被粗大猛烈的视频| 18禁在线播放成人免费| 久久ye,这里只有精品| 亚洲精品久久午夜乱码| 国产午夜精品久久久久久一区二区三区| 在线天堂最新版资源| 亚洲无线观看免费| kizo精华| 热99久久久久精品小说推荐| 97超碰精品成人国产| 日日摸夜夜添夜夜添av毛片| 91精品一卡2卡3卡4卡| 国产精品.久久久| www.av在线官网国产| 国产高清不卡午夜福利| av有码第一页| 亚洲图色成人| 91午夜精品亚洲一区二区三区| 国产毛片在线视频| 91精品一卡2卡3卡4卡| 97超碰精品成人国产| 新久久久久国产一级毛片| 简卡轻食公司| 91精品国产国语对白视频| 搡老乐熟女国产| 精品一区二区免费观看| 国产男女超爽视频在线观看| 久久久欧美国产精品| 欧美日韩亚洲高清精品| 精品国产国语对白av| 久久久久久久久大av| 一级毛片aaaaaa免费看小| 中文字幕免费在线视频6| 久久午夜福利片| av在线播放精品| 好男人视频免费观看在线| 成人18禁高潮啪啪吃奶动态图 | 大陆偷拍与自拍| 欧美激情极品国产一区二区三区 | 夜夜爽夜夜爽视频| 亚洲国产精品国产精品| 丝瓜视频免费看黄片| 欧美精品亚洲一区二区| 老司机影院成人| 91精品伊人久久大香线蕉| xxxhd国产人妻xxx| 国产综合精华液| 五月天丁香电影| 日韩 亚洲 欧美在线| 另类精品久久| 欧美成人午夜免费资源| 国产无遮挡羞羞视频在线观看| 麻豆乱淫一区二区| 少妇精品久久久久久久| 另类亚洲欧美激情| av黄色大香蕉| 精品亚洲成a人片在线观看| 久久久久国产网址| 日日撸夜夜添| 美女内射精品一级片tv| 纯流量卡能插随身wifi吗| 夫妻性生交免费视频一级片| 26uuu在线亚洲综合色| 成人影院久久| 国产精品蜜桃在线观看| 亚洲,欧美,日韩| 久久99一区二区三区| 国产高清国产精品国产三级| 午夜激情福利司机影院| 国产在线一区二区三区精| 成年人免费黄色播放视频| 大香蕉久久成人网| 91在线精品国自产拍蜜月| 亚洲精品第二区| 十八禁高潮呻吟视频| 男女免费视频国产| 精品久久久久久久久av| 国产精品蜜桃在线观看| 人人妻人人爽人人添夜夜欢视频| 极品人妻少妇av视频| 日本爱情动作片www.在线观看| 美女国产高潮福利片在线看| 五月开心婷婷网| 99久久精品国产国产毛片| 亚洲精品国产av蜜桃| 欧美xxⅹ黑人| 亚洲精品国产色婷婷电影| 熟女人妻精品中文字幕| 91精品三级在线观看| 国产精品国产三级专区第一集| 激情五月婷婷亚洲| 欧美日韩av久久| 色5月婷婷丁香| 只有这里有精品99| 免费观看a级毛片全部| 精品国产国语对白av| 婷婷成人精品国产| 亚洲精华国产精华液的使用体验| 少妇的逼好多水| 九色亚洲精品在线播放| 午夜免费观看性视频| 国产成人精品在线电影| 亚洲国产欧美在线一区| 五月天丁香电影| 亚洲av综合色区一区| 夫妻性生交免费视频一级片| 亚洲av男天堂| 久久久久国产精品人妻一区二区| 人妻一区二区av| 少妇人妻 视频| 最近中文字幕高清免费大全6| 99九九在线精品视频| av视频免费观看在线观看| 三级国产精品欧美在线观看| 亚洲精品日韩av片在线观看| 男人爽女人下面视频在线观看| 国产高清三级在线| 国产 一区精品| 午夜免费鲁丝| 亚洲国产成人一精品久久久| 天天躁夜夜躁狠狠久久av| 一级片'在线观看视频| 精品卡一卡二卡四卡免费| 日韩电影二区| 亚洲精品视频女| 精品国产一区二区久久| 亚洲成人手机| 欧美激情 高清一区二区三区| 中文天堂在线官网| 性色avwww在线观看| 国产乱来视频区| 亚洲精华国产精华液的使用体验| 亚洲av男天堂| 在线亚洲精品国产二区图片欧美 | 2018国产大陆天天弄谢| 在线亚洲精品国产二区图片欧美 | 婷婷成人精品国产| 日韩,欧美,国产一区二区三区| 亚洲精品乱久久久久久| 91精品国产九色| 欧美激情 高清一区二区三区| 亚洲av男天堂| 九色亚洲精品在线播放| 夫妻性生交免费视频一级片| 久久久国产一区二区| 纯流量卡能插随身wifi吗| 女人久久www免费人成看片| 日韩中字成人| 亚洲国产精品成人久久小说| 99久久精品一区二区三区| 亚洲美女黄色视频免费看| 精品少妇久久久久久888优播| 18禁动态无遮挡网站| 人成视频在线观看免费观看| 制服人妻中文乱码| 日韩亚洲欧美综合| 成人免费观看视频高清| 爱豆传媒免费全集在线观看| 久久99热6这里只有精品| 尾随美女入室| 日韩伦理黄色片| 亚洲国产色片| 秋霞伦理黄片| 欧美日韩成人在线一区二区| 各种免费的搞黄视频| 欧美日韩视频精品一区| 国产极品粉嫩免费观看在线 | 国产极品天堂在线| 国产一区二区三区av在线| 91久久精品电影网| 国产69精品久久久久777片| 成人二区视频| 亚洲欧美色中文字幕在线| 亚洲无线观看免费| 欧美三级亚洲精品| 91aial.com中文字幕在线观看| 美女主播在线视频| av国产精品久久久久影院| 久久99热这里只频精品6学生| 制服诱惑二区| 爱豆传媒免费全集在线观看| av国产精品久久久久影院| 美女cb高潮喷水在线观看| 三上悠亚av全集在线观看| 亚洲图色成人| 美女中出高潮动态图| 日韩,欧美,国产一区二区三区| 伊人久久精品亚洲午夜| 亚洲精品乱码久久久久久按摩| 国产欧美另类精品又又久久亚洲欧美| 午夜福利视频精品| 丝袜脚勾引网站| 天堂俺去俺来也www色官网| 国产精品一二三区在线看| 一级黄片播放器| 久久久久久久大尺度免费视频| 精品人妻在线不人妻| 国产精品一区二区在线观看99| 亚洲av中文av极速乱| 久久99一区二区三区| 国产熟女午夜一区二区三区 | 麻豆乱淫一区二区| 精品亚洲乱码少妇综合久久| 亚洲欧美色中文字幕在线| 我的老师免费观看完整版| 午夜久久久在线观看| 人妻人人澡人人爽人人| 黄色一级大片看看| 高清av免费在线| 熟妇人妻不卡中文字幕| 高清黄色对白视频在线免费看| 欧美国产精品一级二级三级| 少妇熟女欧美另类| 久久国产精品大桥未久av| 黄片播放在线免费| 国产日韩欧美亚洲二区| 日韩不卡一区二区三区视频在线| 欧美丝袜亚洲另类| 久久99蜜桃精品久久| 国产精品99久久99久久久不卡 | 午夜福利视频精品| 国产在线免费精品| 亚洲精品久久午夜乱码| 十八禁网站网址无遮挡| 校园人妻丝袜中文字幕| 伦精品一区二区三区| 日本猛色少妇xxxxx猛交久久| 亚洲av在线观看美女高潮| 哪个播放器可以免费观看大片| 亚洲精品成人av观看孕妇| 亚洲成人av在线免费| 午夜免费鲁丝| 国产成人一区二区在线| 国产男女内射视频| h视频一区二区三区| 极品人妻少妇av视频| 热re99久久国产66热| 少妇被粗大的猛进出69影院 | 如日韩欧美国产精品一区二区三区 | 欧美人与性动交α欧美精品济南到 | 午夜激情久久久久久久| 国产深夜福利视频在线观看| 久久久亚洲精品成人影院| 国产精品三级大全| 久久婷婷青草| 久久久久久伊人网av| 美女国产高潮福利片在线看| 精品久久久噜噜| 欧美日韩在线观看h| 成人无遮挡网站| www.av在线官网国产| 如日韩欧美国产精品一区二区三区 | 亚洲精品av麻豆狂野| 美女视频免费永久观看网站| 五月伊人婷婷丁香| 一本久久精品| 黑丝袜美女国产一区| av又黄又爽大尺度在线免费看| 九九爱精品视频在线观看| 日本免费在线观看一区| 精品少妇内射三级| 欧美97在线视频| 久久99一区二区三区| 人人澡人人妻人| 久久国产精品男人的天堂亚洲 | 亚洲av中文av极速乱| 中文乱码字字幕精品一区二区三区| 日韩成人伦理影院| 日本av免费视频播放| 韩国高清视频一区二区三区| 亚洲中文av在线| 国产高清有码在线观看视频| 日本欧美视频一区| 我的女老师完整版在线观看| 少妇 在线观看| 国产一区二区三区综合在线观看 | 精品一区在线观看国产| 亚洲av不卡在线观看| 欧美人与性动交α欧美精品济南到 | 内地一区二区视频在线| 国产国语露脸激情在线看| a 毛片基地| 欧美日本中文国产一区发布| 91精品国产国语对白视频| 久久精品熟女亚洲av麻豆精品| 七月丁香在线播放| 欧美亚洲日本最大视频资源| 精品国产露脸久久av麻豆| 日韩中文字幕视频在线看片| 99久久中文字幕三级久久日本| 制服丝袜香蕉在线| 日韩中字成人| 丰满迷人的少妇在线观看| 高清不卡的av网站| 欧美亚洲日本最大视频资源| 两个人免费观看高清视频| 久久免费观看电影| 91久久精品国产一区二区三区| 水蜜桃什么品种好| 中文字幕人妻熟人妻熟丝袜美| 波野结衣二区三区在线| 精品久久久久久久久亚洲| 亚洲欧美日韩卡通动漫| 免费黄色在线免费观看| 国产成人一区二区在线| 国产精品一国产av| 人妻人人澡人人爽人人| 国产精品蜜桃在线观看| 寂寞人妻少妇视频99o| 老司机影院成人| 国产精品一区二区三区四区免费观看| 午夜福利视频精品| av视频免费观看在线观看| 亚洲精品国产色婷婷电影| 亚洲av综合色区一区| 成人亚洲欧美一区二区av| 嫩草影院入口| 久久精品国产亚洲av涩爱| 观看美女的网站| 久久狼人影院| 肉色欧美久久久久久久蜜桃| 一本一本综合久久| 久久亚洲国产成人精品v| 欧美性感艳星| 夜夜爽夜夜爽视频| 一级a做视频免费观看| 亚洲精品一区蜜桃| 国产一区二区在线观看av| 久久久精品区二区三区| 国产精品成人在线| 久久国内精品自在自线图片| 男女无遮挡免费网站观看| 多毛熟女@视频| 亚洲精品乱码久久久v下载方式| 男人操女人黄网站| 女性被躁到高潮视频| 精品熟女少妇av免费看| 91久久精品国产一区二区三区| 国产精品麻豆人妻色哟哟久久| 中文字幕久久专区| 99久久中文字幕三级久久日本| 亚洲国产欧美日韩在线播放| 日韩不卡一区二区三区视频在线| 一边摸一边做爽爽视频免费| 超色免费av| 春色校园在线视频观看| 亚洲欧美一区二区三区黑人 | 国产精品一区二区三区四区免费观看| 丰满迷人的少妇在线观看| 不卡视频在线观看欧美| 久久午夜综合久久蜜桃| 欧美日本中文国产一区发布| 夫妻午夜视频| 女的被弄到高潮叫床怎么办| 菩萨蛮人人尽说江南好唐韦庄| 免费看av在线观看网站| 国产一区二区在线观看日韩| 在线观看www视频免费| av福利片在线| 亚洲成人一二三区av| 欧美日韩精品成人综合77777| 在线 av 中文字幕| 边亲边吃奶的免费视频| 久久午夜福利片| 久久人妻熟女aⅴ| 国产亚洲最大av| 一个人看视频在线观看www免费| 97超视频在线观看视频| 久久这里有精品视频免费| 中国三级夫妇交换| 欧美变态另类bdsm刘玥| 日韩中文字幕视频在线看片| 欧美变态另类bdsm刘玥| 草草在线视频免费看| 欧美精品高潮呻吟av久久| 哪个播放器可以免费观看大片| 亚洲精品国产av蜜桃| 街头女战士在线观看网站| av视频免费观看在线观看| 成人国语在线视频| av线在线观看网站| 亚洲av成人精品一二三区| 亚洲美女搞黄在线观看| av免费观看日本| 亚洲不卡免费看| 免费人妻精品一区二区三区视频| 欧美成人精品欧美一级黄| 欧美亚洲 丝袜 人妻 在线| 91精品国产九色| 成人毛片a级毛片在线播放| 久久精品国产亚洲网站| 亚洲精品亚洲一区二区| 在线亚洲精品国产二区图片欧美 | 三上悠亚av全集在线观看| 国产精品国产三级专区第一集| 亚洲国产成人一精品久久久| 日本与韩国留学比较| 蜜桃国产av成人99| av播播在线观看一区| 亚洲图色成人| 久久精品久久久久久噜噜老黄| 18禁动态无遮挡网站| 国产极品天堂在线| 一本久久精品| 国产极品天堂在线| 亚洲综合精品二区| 亚洲精品久久久久久婷婷小说| 成年人午夜在线观看视频| 国产精品偷伦视频观看了| 亚洲欧美日韩另类电影网站| 又大又黄又爽视频免费| 91精品国产九色| 亚洲一级一片aⅴ在线观看| 久久久久国产精品人妻一区二区| av网站免费在线观看视频| 精品久久久精品久久久| 看免费成人av毛片| 天堂8中文在线网| 精品视频人人做人人爽| 中文欧美无线码| 老司机亚洲免费影院| 看免费成人av毛片| 五月天丁香电影| 欧美精品一区二区免费开放| 国产无遮挡羞羞视频在线观看| 日韩一本色道免费dvd| 亚洲精品av麻豆狂野| 亚洲,欧美,日韩| 狠狠婷婷综合久久久久久88av| 欧美日韩视频高清一区二区三区二| 人妻一区二区av| 观看美女的网站| 色5月婷婷丁香| 黄片播放在线免费| 十八禁高潮呻吟视频| a级毛片免费高清观看在线播放| 欧美日本中文国产一区发布| 久久久久视频综合| 18+在线观看网站| 男女国产视频网站| 成人免费观看视频高清| 黄片播放在线免费| 欧美日韩一区二区视频在线观看视频在线| 夜夜骑夜夜射夜夜干| 亚洲av中文av极速乱| 午夜激情av网站| 欧美bdsm另类| 在线观看免费高清a一片| 美女国产高潮福利片在线看| 永久网站在线| 男女边吃奶边做爰视频| 在线观看www视频免费| 久久久午夜欧美精品| 国产精品久久久久成人av| 成人免费观看视频高清| 国产男人的电影天堂91| 成年人午夜在线观看视频| 午夜福利,免费看| 免费大片18禁| 国产日韩欧美亚洲二区| 日日撸夜夜添| 亚洲精品国产av成人精品| kizo精华| 亚洲美女黄色视频免费看| av国产精品久久久久影院| 日韩,欧美,国产一区二区三区| 最近最新中文字幕免费大全7| 国产精品久久久久久精品电影小说| 少妇被粗大的猛进出69影院 | 久久国内精品自在自线图片| 午夜影院在线不卡| 久久久久久久精品精品| 亚洲在久久综合| 女性被躁到高潮视频| 高清黄色对白视频在线免费看| 亚洲国产成人一精品久久久| 永久免费av网站大全| 久久久久久久国产电影| 成人手机av| 狂野欧美激情性bbbbbb| 亚洲无线观看免费| 国产男人的电影天堂91| 九草在线视频观看| 欧美精品国产亚洲| 五月开心婷婷网| 99九九线精品视频在线观看视频| 免费观看av网站的网址| 热99国产精品久久久久久7| 美女主播在线视频| 纯流量卡能插随身wifi吗| 成人国语在线视频| 十八禁网站网址无遮挡| 91在线精品国自产拍蜜月| 男的添女的下面高潮视频| 国产成人aa在线观看| 亚洲伊人久久精品综合| 中文字幕人妻熟人妻熟丝袜美| 久久综合国产亚洲精品| 黄色欧美视频在线观看| 最近最新中文字幕免费大全7| 91精品三级在线观看| 欧美国产精品一级二级三级| 欧美日韩综合久久久久久| 国产av国产精品国产| 一级,二级,三级黄色视频| 丰满少妇做爰视频| 如日韩欧美国产精品一区二区三区 | 18禁观看日本| 91午夜精品亚洲一区二区三区| 黑人巨大精品欧美一区二区蜜桃 | 国产一区二区在线观看av| 超碰97精品在线观看| 在线观看三级黄色| 国产一级毛片在线| 九九久久精品国产亚洲av麻豆| 亚洲无线观看免费| 欧美国产精品一级二级三级| 熟女电影av网| 中文字幕av电影在线播放| 国产av精品麻豆| 91精品伊人久久大香线蕉| 欧美激情 高清一区二区三区| av卡一久久| 亚洲中文av在线| 国产精品嫩草影院av在线观看| 欧美日本中文国产一区发布| 中文精品一卡2卡3卡4更新| 久久久欧美国产精品| 午夜久久久在线观看| 精品亚洲成国产av| 亚洲精品乱久久久久久| 久久精品国产亚洲av天美| 男女边摸边吃奶| 婷婷色av中文字幕| 哪个播放器可以免费观看大片| 精品亚洲成a人片在线观看| 一边亲一边摸免费视频|