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

    Nonsingular Terminal Sliding Mode Control With Ultra-Local Model and Single Input Interval Type-2 Fuzzy Logic Control for Pitch Control of Wind Turbines

    2021-04-16 03:56:26SaberAbrazehAhmadParvareshSaeidRezaMohseniMeisamJahanshahiZeitouniMeysamGheisarnejadandMohammadHassanKhoobanSeniorMemberIEEE
    IEEE/CAA Journal of Automatica Sinica 2021年3期

    Saber Abrazeh, Ahmad Parvaresh, Saeid-Reza Mohseni, Meisam Jahanshahi Zeitouni, Meysam Gheisarnejad, and Mohammad Hassan Khooban, Senior Member, IEEE

    Abstract—As wind energy is becoming one of the fastestgrowing renewable energy resources, controlling large-scale wind turbines remains a challenging task due to its system model nonlinearities and high external uncertainties. The main goal of the current work is to propose an intelligent control of the wind turbine system without the need for model identification. For this purpose, a novel model-independent nonsingular terminal slidingmode control (MINTSMC) using the basic principles of the ultralocal model (ULM) and combined with the single input interval type-2 fuzzy logic control (SIT2-FLC) is developed for non-linear wind turbine pitch angle control. In the suggested control framework, the MINTSMC scheme is designed to regulate the wind turbine speed rotor, and a sliding-mode (SM) observer is adopted to estimate the unknown phenomena of the ULM. The auxiliary SIT2-FLC is added in the model-independent control structure to improve the rotor speed regulation and compensate for the SM observation estimation error. Extensive examinations and comparative analyses were made using a real-time softwarein-the-loop (RT-SiL) based on the dSPACE 1202 board to appraise the efficiency and applicability of the suggested modelindependent scheme in a real-time testbed.

    I. INTRODUCTION AND PRELIMINARIES

    DUE to the shortage of fossil resources and the environmental issues resulting from greenhouse gas emissions, it seems necessary to pay more attention to renewable energy as an alternative energy source [1]. Wind energy has attracted much attention as an invaluable resource and one of the biggest potentials to meet the global energy demand [2]-[4]. Wind turbine technology, which operates under two configurations (i.e., constant and variable speeds)exploits the wind power to generate usable electrical energy for the customers. In the context of wind power production,the implementation of effective production equipment is necessary to optimize the conversion exploitation of such systems. The generation capacity of a wind turbine (WT)system is usually presented by a curve based on the form of power versus the wind velocity characteristics, as illustrated in Fig.1 [5]-[7]. In the first region, the wind speed is lower than the cut-in value and the turbine is not working. This working area is controlled to obtain the maximum power from the available wind. Region II is between the cut-in and the nominal speed power from the available wind, and its control objective is to capture the maximum available power from the wind flow using variable speed operation of the wind turbine.In region III, which is between nominal speeds and cut-out,the wind turbine system power curve remains constant at the rated power. And finally, the wind speed above the cut-out, or region IV, in which the turbine is powered down and stopped to avoid excessive operating loads.

    Fig.1. The operation modes of wind turbine.

    The reliable and optimum operation of a wind power generation system depends on the control systems applied to the wind turbine operating in different regions. In region III,restrictions to the wind power are imposed by applying pitch angle control at the high speeds of this region [8]. Motivated by the real-world problems, developing a proper pitch angle control methodology, which is straightforward to design a wind turbine plant and also robust in terms of nonlinearities,becomes a significant challenge for the control society.Although, being everywhere in the industrial applications of the WT, the proportional integral derivative (PID) does not optimally perform, since necessary modifications such as extra notch filters are required in their structure to handle the nonlinearity of the wind power systems. To eliminate this problem, numerous robust controls such as backstepping controller [9], model predictive control [10], adaptive controller [11], interval type-2 fuzzy logic control (IT2-FLC)[12], and artificial neural network [13] were applied to design a robust scheme for the pitch angle controllers. For instance,aiming to design a controller with the disturbance rejection and reliable operation in the wind turbine systems, Muhando et al. [14] have proposed a pitch angle regulator-based gainscheduling scheme by minimizing the H∞-norm. Since the feedback gains in [14] are designed by the dynamic specifications, the suggested controller needs to be reconfigured against parametric uncertainty and faults. To overcome the deficiency of linear controllers in the wind turbine context, Boukhezzar and Siguerdidjane developed a nonlinear state feedback controller [15] which offers a right balance between optimal power exploitation and train load flexibility. They claimed that the suggested nonlinear scheme offers more efficient performance than the linear quadratic Gaussian, particularly for the extreme wind speeds.

    The non-linear nature of the wind turbine aerodynamics necessitates adopting nonlinear controllers to regulate the power of the uncertain nonlinear plants, and to characterize the stability requirement. The sliding-mode control (SMC) is one of the most effective nonlinear controllers which has the ability to handle the nonlinear plants subjected to external disturbances [16]. To enhance the regulation capability of the standard SMC and improve some control specifications like finite-time stability, higher accuracy, and quick response,terminal SMC (TSMC) has been introduced [17]. To meliorate the convergence speed and overcome the problems associated with the singularity in the traditional TSMC scheme, a nonsingular version of the TSMC scheme, entitled NTSMC, is developed in [18], [19].

    The developments of the model-based control schemes (e.g.,SMC [5], [16], H-infinity [20], fuzzy adaptive control [21],[22], adaptive neural network [23], [24], etc.) have been a topic of intense research in the wind turbine context in the literature. Although the model-based controllers might regulate the wind turbine power effectively, the mathematical model describing the plant dynamics needs to be developed.The main obstacle with these control methodologies is that the presence of un-modeled dynamics and parametric uncertainties in the wind turbine system are inevitable. In this situation these methodologies cannot produce satisfactory outcomes. In addition, a need for the full dynamic states of the system (as the feedback signal) makes the design of the model-based SMC controllers difficult as some turbine models require high-order controllers [25]. These difficulties provoked researchers to develop the model-independent controllers based on the intelligent PI (iPI) [26] and modelindependent SMC (MISMC) [27] for the indeterminate systems to offer a controller with less dependency on the plant dynamics and which compensates the variable parts of such systems.

    Obtaining an accurate mathematical modeling of two-mass WT system is a challenging task due to the complexity of the mechanics, the chaotic and unpredictable conditions in which they work. The main objective of the current work is to design a model-independent controller (MINTSMC), for optimal capturing the power, that considers the nonlinear nature of wind turbines in the presence of unknown disturbances. Since the suggested controller is developed based on an ultra-local model (ULM) scheme, an observer-based sliding mode (SM)is designed for the estimation of the unknown phenomena(e.g., un-modeled dynamics, etc.). For the effective performance of the SM observer, it is essential that the unknown phenomena be bounded, otherwise, an estimation error appears when the observer coefficients are not properly set. To remove this estimation error, an auxiliary single input IT2 fuzzy PI (SIT2-FPI) [28] is added in this work to the observer-based model-independent control. The main contributions can be summarized as:

    i) According to the knowledge of the wind turbine input and output measurements, an SM observer is designed to estimate the unknown dynamics.

    ii) The pitch angle controller is designed based on the model-independent NTSMC scheme to eliminate the prevalent controller's complexity and enter the required control specifications.

    iii) To compensate for the SM observer estimation error in order to ensure the finite-time convergence and ameliorate the system stability, an auxiliary SIT2-FPI is established in the suggested model-independent controller structure.

    The rest of the work is organized as follows. The model of the studied wind turbine is given in Section II. The structure of the suggested model-independent controller is presented in Section III. Afterward, the real-time outcomes and comparison analysis are prepared in Section IV. In Section V,the conclusion and future scope are included.

    II. WIND TURBINE NONLINEAR MODEL

    The general scheme of a wind turbine with variable speed is illustrated in Fig.2, which comprises an aero turbine, mechanical components (gearbox and generator), and electrical devices (power electronic converter and grid connection).

    As is well known, wind energy is first converted into mechanical energy through the wind turbine blades, and then into electrical energy through the generator, which is consequently transmitted into the power grid or consumed locally. For small analysis, a simplified representation of the 2-mass model of the wind turbine plant is presented in Fig.3[8], [15], [25].

    The mechanical power extracted by the wind turbine is expressed by the following equation [15]:

    Fig.2. Global scheme of wind turbines.

    Fig.3. 2-mass dynamic model of aero turbine.

    where ρ is the air density, R is the wind turbine rotor radius, v is the wind speed and CPis the power coefficient of the wind turbines which is a nonlinear function of the pitch angle β and tip-speed ratio λ. The parameter λ is calculated by the blade tip speed and wind speed upstream of the rotor as

    with ωrbeing the rotor angular speed. The curves of CPwith respect to the multiple values of β are sketched in Fig.4.

    The mechanical power Prcan also be rewritten as

    where Tais the aerodynamic torque, given by

    with Cq(λ,β) being the torque coefficient. Furthermore, it can be readily obtained that

    Cq(λ,β) is generally determined experimentally and provided by the manufacturers [29]. In practice, an acceptable numerical approximation can be adopted.

    where

    The nonlinear wind turbine model is shown in a generalized nonlinear form as follows [8]:

    The state vector x, control input uand nonlinear vector G(x)are defined as

    These parameters are listed in Table I.

    TABLE I PARAMETERS OF THE 1.5 MW EXPERIMENTAL WIND TURBINE

    The current work aims to develop a novel modelindependent pitch angle controller for rotor speed regulation at region III of the wind turbine operation by restricting the power derived from the wind turbine.

    In the next sub-section, the overall structure of the suggested model-independent control scheme and the procedure for designing the controller gains are presented.

    III. DESIGN OF MODEL-INDEPENDENT NTSMC BASED IT2-FLC

    Based on the fundamental concept of the ULM control scheme, a combination of SIT2-FPI and model-independent NTSMC scheme is introduced which includes three subcontrollers: i) Firstly, an NTSMC controller is adopted in the model-independent framework to reduce the dependency of the existing controller to the wind turbine model and handle the plant’s nonlinearity; ii) Then, an SM observer is designed to estimate and eliminate the time-varying disturbances and un-modeled terms; iii) And additionally, an auxiliary controller-based SIT2-FPI structure is incorporated into the ULM to compensate the estimation error of SM observer and to enhance the overall transient behavior.

    A. Model-Independent NTSMC Technique

    Definition 1: For a general nonlinear plant, the modelindependent control that is described in [26], [30] and the ULM form which is approximately defined by (11) can be adopted to express the system dynamics.

    where y (t) is the output term, v is the derivative order, ε is the unknown structure of the plant as well as any disturbances which are estimated based on the input signal u(t) and output y(t) , α ∈R is an assumed constant.

    Definition 2: To determine the sliding surface, the new state variables x1and x2are introduced as

    According to (12), the state-space equations are given as

    Based on the SMC framework, a second-order NTSM is presented in (14) with the designed parameter γ , and the p and q are odd defined as1 <p/q <2.

    By differentiating (14), one can obtain

    Theorem 1 [27]: For the state-space equation given in (12),if the NTSMC manifold is considered as (14) and control law is formulated as (16), then one has e →0 in finite time.

    Proof: The Lyapunov function, as usual, is chosen as

    If (17) is differentiated with respect to time, and (15) is substituted into it, one can obtain

    when η1>‖ε‖+μ(μ >0), (19) will be transformed into the following inequality:

    Thus, the stability of the controller will be satisfied by the demand of the Lyapunov theory. ■

    B. Design of SM Observer

    Definition 3: Let us define the observer error as

    In the next step, by subtracting relation (11) from (22), the dynamic error equation can be obtained, given as

    Proof: Let us define the Lyapunov function as

    If (24) is differentiated with respect to time, the following equation is attained:

    Selecting k by k >‖ε‖+η(η >0) yields

    The structure of the established model-independent scheme,which comprises the MINTSM controller with SM observer,is illustrated in Fig.5 [27].

    C. Auxiliary Controller-Based SIT2-FLC Structure

    Fig.5. Block diagram of MINTSMC scheme with SM observer.

    In this sub-section, the performance output of the MINTSMC controller is improved by adding an extra control signal to the system to reduce the observer estimation error.The auxiliary controller which is exploited in this paper, the SIT2-FPI controller, is depicted in Fig.6. In this diagram, kpand kiare the coefficients of the baseline PI controller and kedenotes the input scaling factor (SF). Based on the given scheme, the SIT2-FPI control signal uSIT2is utilized by the output action of the controller (φo) that is shown in (27) [32], [33].

    where χ denotes the learning rate; c1and c2denote acceleration coefficients; R1and R2indicate random numbers within [0 1]. The bounds of the model-independent controller coefficients , , , , are chosen to be 0 and 200.

    Remark 1: In the following, the advantages of the suggested model-independent scheme are summarized.

    i) In comparison with the model-based control strategies(model predictive control, backstepping controllers), the suggested scheme is developed based on the measured input/output (I/O) data of the system which remove the need to model identification,

    ii) In compassion with the conventional SMC and TSMC approaches, the quick time convergence and the elimination of singularity problems are achieved by the NTSMC.

    iii) In compassion with the MINTSMC approach with SM observer, an auxiliary controller-based SIT2-FLC is established in the MINTSMC structure which removes the observer error and ameliorates the pitch angle performance of the WT plant.

    iv) In comparison with the conventional type-1 and interval type-2 FLCs with two inputs [22], [38]-[41], the SIT2-FLC controller does not need to adjust the rule base and membership functions which is valuable in practical applications.

    Ξ={γ η1η2kpki}

    IV. EXPERIMENTAL RESULTS

    The essence of designing ULM control based on the modelindependent scheme is to offer optimal performance for the wind turbine plant under the various level of wind speed variations without the need for model identification. The suggested model-independent scheme was examined by the wind turbine modeled as a 2-mass model using a real-time software-in-the-loop (RT-SiL) simulation. The SiL platform integrates the compiled source code into a functional model simulation to offer a virtual simulation environment for the examination of the control methodologies for complex problems. In comparison with the traditional off-line simulations, the SiL platform provides a real-time simulation to modulate some delays and small errors as well as the effect of the computational burden on the transient performance. The necessary steps for the realization of the SiL setup based on the dSPACE 1202 board are illustrated in Fig.10. In this application, the model of the test system along with the control strategy was run in real-time using the Runge Kutta’s solver with a sampling frequency of 10 kHz [1]. The parameters of the 2-mass WT model in this work are considered from [8].

    Fig.10. Steps for executing RT-SiL setup.

    The transient outcomes of the MINTSMC based SIT2-FPI controller are verified under four typical scenarios of the wind turbine, as i) the step changes of the wind, ii) the random fluctuation of the wind speed iii) parametric uncertainty in the turbine model, and iv) simultaneous change in wind speed and system parameters. Since region III is selected as the operation region of the pitch control, the wind speed for the tested scenarios is considered in the range of [12 m/s 14 m/s].In the comparative studies of the real-time setup, the transient outcomes of the suggested model-independent controller are compared with those of MINTSMC, N-PI [8], and PI [8]controllers.

    Initially, the dynamic performance of the suggested modelindependent controlled wind turbine is investigated by applying the multi-step wind speed variation shown in Fig.11.The optimization issue with the number of particles Np= 50 was solved for the MINTSMC based A-SIT2-FPI (with setting ζIT2= 0.30), the MINTSMC based M-SIT2-FPI (with setting ζIT20.55) and the MINTSMC based S-SIT2-FPID (with setting ζIT2= 0.85). The convergence specifications of the best-explored solutions for the three SIT2-FPI modes over 30 iterations are illustrated in Fig.12. It is seen that the convergence of F(t) under the A-SIT2-FPI action is recognized to be best among the three compared SIT2-FPI controllers.

    Fig.11. Multi-step variation of wind speed.

    Scenario 1 (Dynamic Performance With the Step Wind Speed): The SiL responses for the MINTSMC based SIT2-FPI, MINTSMC, N-PI, and PI controllers under the multi-step speed disturbance of Fig.11, which comprises the outcomes of rotor speed, shaft twist angle, pitch angle and mechanical power are presented in Figs. 13(a), 13(b), 13(c), and 13(d),respectively.

    Fig.12. Convergence specifications of MINTSMC based SIT2-FPI with(a) A-CC IT2, (b) M-CC IT2, (c) S-CC IT2.

    The real-time outcomes reveal that the N-PI and PI controllers experience large overshoot with undamped oscillations. Also, it is seen that, in spite of the MINTSMC approach which exhibits acceptable robust behavior, there remains a non-zero steady-state error in the test system outcomes. In contrast, the dynamic outcomes reveal that with the suggested controller (realized based on the SIT2-FPI), the system outcomes are alleviated with small deviations in a shorter time and better satisfy the steady-state requirements than the other three controllers.

    Scenario 2 (Dynamic Performance With the Wind Speed Fluctuation): To validate the feasibility of the suggested model-independent scheme against the real-world wind, a random fluctuation of wind speed is applied to the concerned 2-mass model wind turbine, and its corresponding profile is sketched in Fig.14(a). The comparative rotor speed outcomes of the four-pitch angle controllers under the real wind speed data are illustrated in Fig.14(b).

    From Fig.14(b), it is seen that the N-PI and PI controllers undergo undesirable overshoot and large settling time.Moreover, the MINTSMC requires a relatively long regulation time to derive the speed rotor deviation to its reference value,while with adding the auxiliary SIT2-FPI controller, superior outcomes are offered by eliminating the SM observer estimation error and improving the steady-state requirements.

    Fig.13. Time-domain behavior comparison of the suggested scheme,MINTSMC, N-PI, and PI according to multi-step wind test. (a) rotor speed,(b) twist angle, (c) pitch angle, (d) mechanical power.

    Scenario 3 (Robustness Against Parametric Uncertainty):In order to examine the influence of the parametric uncertainties, some parameters of the WT are changed as: Rb= +25%, Jr= +50 % and tβ= +75% . The corresponding two standard error measurement criteria including mean square error (MSE) and root MSE (RMSE) after applying the variations to the test system are reported in Table III. Looking at the performance indices presented in Table III, one can observe that with the suggested controller, the calculated values of the error measurement indices in all the test validation cases are best, consequently, superior robustness is reached while comparing with the others. Results of the rotor speed corresponding to the large positive parametric uncertainties for the various designed pitch angle controllers are illustrated in Fig.15(a)-(d).

    Remark 2: From the outcomes of the four scenarios, it is explored that with the application of PI and N-PI controllers,the responses of the WT system experience a large deviation.Moreover, despite the MINTSMC approach that can stabilize the outcomes of studied 2-mass WT plants, the system responses are not optimal under the wind speed changes and system uncertainties. However, with the suggested MINTSMC based SIT-FLC controller, the system outcomes are further alleviated than other designed pitch angle controllers.

    Fig.14. The time-domain response of suggested scheme, MINTSMC, N-PI,and PI according to random fluctuations of wind (a) random wind speed,(b) rotor speed.

    TABLE III COMPASSION OF MSE AND RMSE STANDARDS FOR THE PARAMETRIC UNCERTAINTIES

    Scenario 4 (Robustness Against Wind Speed Changes and Perturbed Condition of Air Density): Additionally, to ascertain the robustness of the suggested model-independent scheme in a more severe condition, the effect of the change in wind speed and uncertainty in the air density are tested simultaneously. Absolute maximum wind speed deviation when the wind speed is changed from 12 m/s to 24 m/s, and the value of ρ is changed from nominal condition to +50%(i.e., it is changed from 1.2250 kg/m3to 1.8375 kg/m3), is shown in Fig.16.

    V. CONCLUSION

    The present work concentrates on developing a combination of MINTSMC and SIT2-FPI controllers for pitch control of a wind turbine plant. The main outcomes of this research project are:

    Fig.15. Time-domain outcomes according to uncertain parameters (a) suggested scheme, (b) MINTSMC, (c) N-PI, (d) PI.

    Fig.16. Absolute maximum wind speed deviation against wind speed change.

    i) An auxiliary SIT2-FL is incorporated into the ULM control scheme to ameliorate the overall performance of an existing MINTSMC controller and eliminate the SM observer estimation error.

    ii) The performance of the designed pitch angle controller in tracking the wind speed rotor under both cases of multi-step and fluctuation wind disturbances is verified.

    iii) The comprehensive tests based on the SiL platform are conducted on a simplified 2-mass model of a wind turbine to validate the applicability of the model-independent controller from the real-time perspective.

    iv) Simultaneous consideration of the wind speed change and uncertainty in the air density of the nonlinear system can be regarded as another outcome of this research.

    Since the chattering adversely affects the system performance which deteriorates the system performance and imposes instabilities; as future work, the chattering attenuation NTSMC can be developed to eliminate this phenomenon. For this purpose, the switching function in the SMC controller should be designed in such a way that mitigates the effect of chattering over time.

    精品一区二区三区视频在线观看免费| 极品教师在线视频| 人妻少妇偷人精品九色| 国产精品爽爽va在线观看网站| av在线观看视频网站免费| 99久久久亚洲精品蜜臀av| 22中文网久久字幕| 国产高清不卡午夜福利| 99热这里只有是精品在线观看| 国产成人a∨麻豆精品| 国内揄拍国产精品人妻在线| 国产国拍精品亚洲av在线观看| 狂野欧美激情性xxxx在线观看| 少妇猛男粗大的猛烈进出视频 | 国产精品乱码一区二三区的特点| 欧美绝顶高潮抽搐喷水| 精品人妻一区二区三区麻豆 | 亚洲国产精品合色在线| 亚洲国产高清在线一区二区三| 一个人看的www免费观看视频| 狂野欧美白嫩少妇大欣赏| 啦啦啦啦在线视频资源| 国产探花在线观看一区二区| 一个人看视频在线观看www免费| 国产白丝娇喘喷水9色精品| 免费无遮挡裸体视频| 亚洲成人av在线免费| 99久久成人亚洲精品观看| 在线观看美女被高潮喷水网站| 少妇裸体淫交视频免费看高清| 91在线精品国自产拍蜜月| 日韩欧美精品免费久久| 成人毛片a级毛片在线播放| 久久人人精品亚洲av| 国产精品av视频在线免费观看| 丰满的人妻完整版| 国产一区二区在线观看日韩| 国产精品av视频在线免费观看| 长腿黑丝高跟| 一区二区三区四区激情视频 | 日韩欧美 国产精品| 神马国产精品三级电影在线观看| 国内揄拍国产精品人妻在线| 99在线人妻在线中文字幕| 国国产精品蜜臀av免费| 久久热精品热| 午夜福利18| 色在线成人网| 久久久久久久久大av| 国产男人的电影天堂91| 国产高清有码在线观看视频| 久久国内精品自在自线图片| 中国美女看黄片| 欧美xxxx性猛交bbbb| 99热这里只有精品一区| 国产激情偷乱视频一区二区| 99久久九九国产精品国产免费| 国产精品一区二区性色av| 最新中文字幕久久久久| 嫩草影院精品99| 久久九九热精品免费| 免费黄网站久久成人精品| 搞女人的毛片| 1000部很黄的大片| 亚洲美女视频黄频| 国产精品日韩av在线免费观看| 成人午夜高清在线视频| 一级黄片播放器| 十八禁国产超污无遮挡网站| 美女被艹到高潮喷水动态| 网址你懂的国产日韩在线| 久久热精品热| 美女免费视频网站| 婷婷精品国产亚洲av在线| 成人综合一区亚洲| 日韩在线高清观看一区二区三区| 色哟哟哟哟哟哟| 亚洲精品亚洲一区二区| 亚洲精品日韩av片在线观看| 美女免费视频网站| 成人漫画全彩无遮挡| 欧美一区二区国产精品久久精品| 久久久久九九精品影院| 人妻久久中文字幕网| 晚上一个人看的免费电影| 可以在线观看毛片的网站| av在线亚洲专区| 亚洲最大成人中文| 久久人人精品亚洲av| 成人特级av手机在线观看| 国产精品女同一区二区软件| 亚洲四区av| 一级av片app| 性欧美人与动物交配| 看片在线看免费视频| 亚洲中文日韩欧美视频| 国产精品国产高清国产av| 深夜a级毛片| 亚洲欧美日韩无卡精品| 久久久精品94久久精品| 亚州av有码| 亚洲av熟女| 美女 人体艺术 gogo| 中出人妻视频一区二区| 久久久久久久久久久丰满| 午夜视频国产福利| 天天一区二区日本电影三级| 12—13女人毛片做爰片一| 99久久精品热视频| 免费观看人在逋| 午夜精品在线福利| 女的被弄到高潮叫床怎么办| 日本一本二区三区精品| 免费av不卡在线播放| 亚洲美女黄片视频| 天堂动漫精品| 我要看日韩黄色一级片| 国模一区二区三区四区视频| 老司机午夜福利在线观看视频| 国产伦精品一区二区三区四那| 老熟妇乱子伦视频在线观看| 国产三级在线视频| 亚洲精品国产成人久久av| 亚洲精品久久国产高清桃花| 欧美中文日本在线观看视频| 国产欧美日韩精品一区二区| av福利片在线观看| 欧美日韩精品成人综合77777| www.色视频.com| 久久久久久久午夜电影| 亚洲婷婷狠狠爱综合网| 成人鲁丝片一二三区免费| 精品久久久噜噜| 亚洲人成网站在线播放欧美日韩| 久久国产乱子免费精品| 亚洲最大成人av| 色视频www国产| 日韩在线高清观看一区二区三区| 永久网站在线| 欧美性猛交╳xxx乱大交人| 成熟少妇高潮喷水视频| 成人av一区二区三区在线看| 国产黄色小视频在线观看| 九九爱精品视频在线观看| 亚洲中文字幕日韩| 我要搜黄色片| 国产精品一区二区免费欧美| 此物有八面人人有两片| 成人无遮挡网站| 免费人成在线观看视频色| 亚洲乱码一区二区免费版| 少妇猛男粗大的猛烈进出视频 | 99热精品在线国产| 国产精品福利在线免费观看| 国产亚洲av嫩草精品影院| 免费一级毛片在线播放高清视频| 一a级毛片在线观看| 性色avwww在线观看| 天堂网av新在线| 大型黄色视频在线免费观看| 亚洲婷婷狠狠爱综合网| 国产精品不卡视频一区二区| 免费电影在线观看免费观看| 亚洲一区高清亚洲精品| 亚洲欧美精品自产自拍| 国产高清有码在线观看视频| 日韩一区二区视频免费看| 无遮挡黄片免费观看| 18禁在线无遮挡免费观看视频 | 九九热线精品视视频播放| 91精品国产九色| 国产精品久久久久久久电影| 免费观看人在逋| 尤物成人国产欧美一区二区三区| 国产三级在线视频| 亚洲欧美日韩高清专用| 大香蕉久久网| 狂野欧美激情性xxxx在线观看| 精品久久久久久成人av| 亚洲av中文av极速乱| 免费人成视频x8x8入口观看| 日韩在线高清观看一区二区三区| 久久中文看片网| 午夜激情福利司机影院| av.在线天堂| 狂野欧美激情性xxxx在线观看| 免费黄网站久久成人精品| 两个人的视频大全免费| 少妇猛男粗大的猛烈进出视频 | 国内精品美女久久久久久| 国产精品精品国产色婷婷| 国产久久久一区二区三区| 人妻夜夜爽99麻豆av| 午夜亚洲福利在线播放| 国产精品电影一区二区三区| 成人国产麻豆网| 精品一区二区三区视频在线| h日本视频在线播放| 变态另类丝袜制服| 如何舔出高潮| 久久热精品热| 在线免费观看的www视频| 青春草视频在线免费观看| 亚洲av五月六月丁香网| 91av网一区二区| 乱人视频在线观看| 日本撒尿小便嘘嘘汇集6| 国产伦精品一区二区三区四那| 老熟妇仑乱视频hdxx| 床上黄色一级片| 全区人妻精品视频| 亚洲综合色惰| 亚洲高清免费不卡视频| 一卡2卡三卡四卡精品乱码亚洲| 人妻制服诱惑在线中文字幕| 日本免费一区二区三区高清不卡| 亚洲精品粉嫩美女一区| 六月丁香七月| 国产在视频线在精品| 欧美三级亚洲精品| 精品乱码久久久久久99久播| 日韩亚洲欧美综合| 97超碰精品成人国产| 亚洲美女黄片视频| a级毛片a级免费在线| 一级黄片播放器| 欧美激情国产日韩精品一区| av国产免费在线观看| 亚洲高清免费不卡视频| 中文字幕久久专区| 久久久午夜欧美精品| 精品久久国产蜜桃| 国产淫片久久久久久久久| 国内精品久久久久精免费| 1000部很黄的大片| 久久精品影院6| 欧美色视频一区免费| 精品久久久噜噜| 91久久精品电影网| 欧美性猛交黑人性爽| 久久人人爽人人片av| 免费av观看视频| 国产一区二区三区av在线 | 免费大片18禁| 日韩大尺度精品在线看网址| 一进一出好大好爽视频| 精品日产1卡2卡| 一区二区三区高清视频在线| 亚洲精品影视一区二区三区av| 国产精品人妻久久久久久| 我的女老师完整版在线观看| 欧美另类亚洲清纯唯美| 人妻夜夜爽99麻豆av| 少妇猛男粗大的猛烈进出视频 | 中文字幕av成人在线电影| 国产精品久久视频播放| 91在线观看av| 有码 亚洲区| 亚洲国产精品国产精品| eeuss影院久久| 51国产日韩欧美| 日本成人三级电影网站| 深爱激情五月婷婷| 久久久a久久爽久久v久久| 天天躁夜夜躁狠狠久久av| 久久综合国产亚洲精品| 老师上课跳d突然被开到最大视频| av在线蜜桃| 男女边吃奶边做爰视频| 国产69精品久久久久777片| 白带黄色成豆腐渣| 久久中文看片网| 亚洲人与动物交配视频| 搡老熟女国产l中国老女人| 亚洲精品日韩av片在线观看| 国产精品久久久久久av不卡| 国产乱人偷精品视频| 亚洲中文字幕日韩| 精品久久国产蜜桃| 国产69精品久久久久777片| 网址你懂的国产日韩在线| 久久中文看片网| av女优亚洲男人天堂| 成人一区二区视频在线观看| 国产精品久久久久久亚洲av鲁大| 嫩草影院精品99| 如何舔出高潮| 国产亚洲精品久久久久久毛片| 亚洲精品色激情综合| 国产中年淑女户外野战色| 熟妇人妻久久中文字幕3abv| 国产私拍福利视频在线观看| 人妻少妇偷人精品九色| 日韩,欧美,国产一区二区三区 | 成人国产麻豆网| 夜夜夜夜夜久久久久| 成人av一区二区三区在线看| 1024手机看黄色片| 三级国产精品欧美在线观看| 久久午夜亚洲精品久久| 国产极品精品免费视频能看的| 嫩草影院入口| 白带黄色成豆腐渣| 少妇猛男粗大的猛烈进出视频 | 日韩一本色道免费dvd| 春色校园在线视频观看| 九九在线视频观看精品| 插阴视频在线观看视频| 黄色日韩在线| 国产精品久久视频播放| 欧美最新免费一区二区三区| av卡一久久| 一个人看视频在线观看www免费| 久久久精品94久久精品| 精品久久久久久久人妻蜜臀av| 深夜精品福利| 99热这里只有是精品50| 97超视频在线观看视频| 真实男女啪啪啪动态图| 丰满人妻一区二区三区视频av| 欧美高清性xxxxhd video| 亚洲美女黄片视频| 女人被狂操c到高潮| 亚洲av中文av极速乱| 美女cb高潮喷水在线观看| 国产精品人妻久久久久久| 男女做爰动态图高潮gif福利片| 97热精品久久久久久| 亚洲在线自拍视频| 寂寞人妻少妇视频99o| 亚洲国产欧美人成| 欧美成人免费av一区二区三区| 国产高清不卡午夜福利| 国产高清视频在线观看网站| 美女 人体艺术 gogo| 伦精品一区二区三区| 菩萨蛮人人尽说江南好唐韦庄 | 黄色日韩在线| 久久久久久国产a免费观看| 成人亚洲精品av一区二区| 三级国产精品欧美在线观看| 激情 狠狠 欧美| 久久草成人影院| 女同久久另类99精品国产91| 丰满乱子伦码专区| 国产亚洲欧美98| 精品熟女少妇av免费看| 精品免费久久久久久久清纯| 别揉我奶头 嗯啊视频| 亚洲中文字幕一区二区三区有码在线看| 亚洲最大成人av| 97在线视频观看| a级毛色黄片| 欧美日韩综合久久久久久| 有码 亚洲区| 九九爱精品视频在线观看| 一进一出抽搐gif免费好疼| 亚洲五月天丁香| 天美传媒精品一区二区| 中文字幕免费在线视频6| 精品一区二区三区视频在线观看免费| 黑人高潮一二区| 欧美色视频一区免费| 国产精品国产高清国产av| 麻豆国产97在线/欧美| 久久久国产成人精品二区| 国产精品一区二区性色av| 美女xxoo啪啪120秒动态图| 又爽又黄a免费视频| 波多野结衣高清无吗| 中国美女看黄片| 国产在线男女| 欧美bdsm另类| 午夜爱爱视频在线播放| 免费无遮挡裸体视频| 丰满乱子伦码专区| 美女高潮的动态| 日本黄色视频三级网站网址| 我要看日韩黄色一级片| 中文字幕人妻熟人妻熟丝袜美| 成人毛片a级毛片在线播放| 听说在线观看完整版免费高清| 国产又黄又爽又无遮挡在线| 又粗又爽又猛毛片免费看| 久久久久国内视频| 国产成人一区二区在线| 在线免费观看的www视频| 免费看av在线观看网站| 午夜福利在线观看免费完整高清在 | 晚上一个人看的免费电影| 精华霜和精华液先用哪个| 人妻久久中文字幕网| 嫩草影视91久久| 免费观看的影片在线观看| 久久鲁丝午夜福利片| 国产极品精品免费视频能看的| 久久欧美精品欧美久久欧美| 亚洲五月天丁香| 国产伦在线观看视频一区| 九九久久精品国产亚洲av麻豆| 麻豆国产av国片精品| 久久久久久久久久成人| 国产亚洲91精品色在线| 少妇的逼好多水| 亚洲欧美清纯卡通| 色5月婷婷丁香| 亚洲国产欧洲综合997久久,| 成年av动漫网址| 人人妻,人人澡人人爽秒播| 99riav亚洲国产免费| 欧美人与善性xxx| 九色成人免费人妻av| 别揉我奶头 嗯啊视频| 亚洲乱码一区二区免费版| 99在线视频只有这里精品首页| 欧美+日韩+精品| 亚洲18禁久久av| 亚洲不卡免费看| 国产精品久久久久久亚洲av鲁大| 毛片女人毛片| 午夜免费男女啪啪视频观看 | 日本撒尿小便嘘嘘汇集6| 永久网站在线| 男女视频在线观看网站免费| 人人妻人人看人人澡| 99久久中文字幕三级久久日本| 91精品国产九色| 男人舔奶头视频| 午夜精品在线福利| 亚洲经典国产精华液单| 国产精品女同一区二区软件| 日日摸夜夜添夜夜添小说| 九九久久精品国产亚洲av麻豆| 亚洲天堂国产精品一区在线| 国产成人福利小说| 欧美人与善性xxx| 少妇猛男粗大的猛烈进出视频 | 男女做爰动态图高潮gif福利片| 国产综合懂色| 欧美一区二区国产精品久久精品| 国产探花极品一区二区| 99久久无色码亚洲精品果冻| 不卡一级毛片| 97超级碰碰碰精品色视频在线观看| 亚洲精品日韩av片在线观看| 久久精品国产鲁丝片午夜精品| 久久久久久久久久黄片| 亚洲av中文字字幕乱码综合| 男女之事视频高清在线观看| 免费av不卡在线播放| 久久国内精品自在自线图片| 国产女主播在线喷水免费视频网站 | 亚洲最大成人中文| 高清毛片免费观看视频网站| 22中文网久久字幕| 99国产精品一区二区蜜桃av| 一级黄色大片毛片| 国产精品一区www在线观看| av中文乱码字幕在线| 成人亚洲精品av一区二区| 91麻豆精品激情在线观看国产| 99热精品在线国产| 日韩成人av中文字幕在线观看 | 亚洲最大成人手机在线| 国产黄片美女视频| 日日摸夜夜添夜夜爱| 亚洲丝袜综合中文字幕| 欧美激情国产日韩精品一区| 成年女人永久免费观看视频| av黄色大香蕉| 婷婷六月久久综合丁香| 日本三级黄在线观看| 色视频www国产| 免费观看精品视频网站| 亚洲自偷自拍三级| 中文字幕免费在线视频6| 久久精品国产自在天天线| 成年免费大片在线观看| 亚洲精品日韩在线中文字幕 | 久久久久国内视频| 一个人看的www免费观看视频| 嫩草影院新地址| 久久精品国产亚洲av涩爱 | 男女视频在线观看网站免费| 国内精品一区二区在线观看| 观看美女的网站| 久久久久久久久久成人| 成人漫画全彩无遮挡| 中国美白少妇内射xxxbb| a级毛片免费高清观看在线播放| 国产欧美日韩精品一区二区| 哪里可以看免费的av片| 亚洲欧美日韩高清在线视频| 亚州av有码| 亚洲欧美日韩高清在线视频| 久久久成人免费电影| 最近最新中文字幕大全电影3| 欧美一区二区国产精品久久精品| 国产精品一区二区三区四区免费观看 | 中文字幕av在线有码专区| 成人午夜高清在线视频| 国产免费一级a男人的天堂| 欧美3d第一页| 免费电影在线观看免费观看| 国产精品不卡视频一区二区| 中国美女看黄片| 91午夜精品亚洲一区二区三区| 欧美成人免费av一区二区三区| 亚洲精品色激情综合| 中文资源天堂在线| 一级毛片我不卡| 日韩精品青青久久久久久| 免费黄网站久久成人精品| 99在线人妻在线中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 99久久精品国产国产毛片| 国产男靠女视频免费网站| 又爽又黄a免费视频| a级毛片免费高清观看在线播放| 99久久九九国产精品国产免费| 国产片特级美女逼逼视频| 国产真实伦视频高清在线观看| 亚洲av电影不卡..在线观看| АⅤ资源中文在线天堂| 黄色配什么色好看| 麻豆国产av国片精品| 18禁黄网站禁片免费观看直播| 又黄又爽又刺激的免费视频.| 久久久精品94久久精品| 亚洲精品一卡2卡三卡4卡5卡| 可以在线观看的亚洲视频| 欧美bdsm另类| 天堂av国产一区二区熟女人妻| 少妇猛男粗大的猛烈进出视频 | 九色成人免费人妻av| 久久久久国产网址| 午夜爱爱视频在线播放| 精品国产三级普通话版| 亚洲中文日韩欧美视频| 97热精品久久久久久| 婷婷精品国产亚洲av在线| av黄色大香蕉| 99九九线精品视频在线观看视频| 熟女电影av网| 成人鲁丝片一二三区免费| 深夜精品福利| 久久久久国产网址| 亚洲成人久久爱视频| 国产精品av视频在线免费观看| 日韩大尺度精品在线看网址| av在线播放精品| 99久久久亚洲精品蜜臀av| 成人亚洲精品av一区二区| 国产一区二区在线av高清观看| 国产精品一区二区三区四区免费观看 | 免费无遮挡裸体视频| 国产黄色视频一区二区在线观看 | 婷婷精品国产亚洲av在线| 午夜久久久久精精品| av在线天堂中文字幕| 夜夜看夜夜爽夜夜摸| 一级av片app| 亚洲熟妇熟女久久| 国产毛片a区久久久久| 国产精品美女特级片免费视频播放器| 久久精品国产亚洲av香蕉五月| 97超碰精品成人国产| 男人的好看免费观看在线视频| 亚洲欧美日韩无卡精品| 亚洲人成网站高清观看| 岛国在线免费视频观看| 日韩人妻高清精品专区| 国产精品永久免费网站| 舔av片在线| 国产午夜福利久久久久久| 老司机午夜福利在线观看视频| 少妇的逼好多水| 亚洲经典国产精华液单| 麻豆精品久久久久久蜜桃| 97超视频在线观看视频| 久久久久国产网址| 欧美xxxx性猛交bbbb| 国产私拍福利视频在线观看| 特级一级黄色大片| 国产人妻一区二区三区在| 久久99热这里只有精品18| 欧美日本视频| 午夜精品在线福利| 久久精品夜夜夜夜夜久久蜜豆| 免费av观看视频| 久久九九热精品免费| 色噜噜av男人的天堂激情| 午夜免费男女啪啪视频观看 | 蜜臀久久99精品久久宅男| 国产高清三级在线| 欧美zozozo另类| 97碰自拍视频| 国产 一区 欧美 日韩| 狠狠狠狠99中文字幕| 日本免费一区二区三区高清不卡| 又爽又黄a免费视频| 校园人妻丝袜中文字幕| 国产高清视频在线播放一区| 九九爱精品视频在线观看| 俄罗斯特黄特色一大片| 悠悠久久av| 精品人妻一区二区三区麻豆 | 少妇的逼水好多| 国产成人福利小说| 在线观看免费视频日本深夜| 蜜桃亚洲精品一区二区三区| 一级毛片电影观看 | 99热这里只有是精品50| 免费人成视频x8x8入口观看| 人人妻人人澡欧美一区二区|