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

    Design of a Robust Optimal Decentralized PI Controller Based on Nonlinear Constraint Optimization For Level Regulation: An Experimental Study

    2020-02-29 14:20:14SoumyaRanjanMahapatroBidyadharSubudhiandSandipGhosh
    IEEE/CAA Journal of Automatica Sinica 2020年1期

    Soumya Ranjan Mahapatro,, Bidyadhar Subudhi,, and Sandip Ghosh,

    Abstract—This paper presents the development of a new robust optimal decentralized PI controller based on nonlinear optimization for liquid level control in a coupled tank system. The proposed controller maximizes the closed-loop bandwidth for specified gain and phase margins, with constraints on the overshoot ratio to achieve both closed-loop performance and robustness. In the proposed work, a frequency response fitting model reduction technique is initially employed to obtain a first order plus dead time (FOPDT) model of each higher order subsystem. Furthermore, based on the reduced order model, a proposed controller is designed. The stability and performance of the proposed controller are verified by considering multiplicative input and output uncertainties. The performance of the proposed optimal robust decentralized control scheme has been compared with that of a decentralized PI controller. The proposed controller is implemented in real-time on a coupled tank system. From the obtained results,it is shown that the proposed optimal decentralized PI controller exhibits superior control performance to maintain the desired level, for both the nominal as well as the perturbed case as compared to a decentralized PI controller.

    NOMENCLATURE

    ACross sectional area of the tanks.

    a12Cross-sectional area of the pipe from tank 1 to tank 2.

    h1water level in Tank 1

    h2Water level in Tank 2

    a1Outlet area of tank 1

    a2Outlet area of tank 2

    gGravitational constant

    η Constant relating the control input with water inflow from the pumps

    ωciPhase crossover frequency

    τiiOverall time constant

    AmGain margin

    φmPhase margin

    PTMaximum closed loop amplitude ratio

    om*Upper bound of maximum closed loop amplitude ratio

    I. INTRODUCTION

    AUTOMATIC regulation of precise liquid level is an important issue in many process industries such as pharmaceutical, chemical, and food processing industries. In general,most industrial processes are multi input multi output (MIMO)systems. The coupled tank system (CTS) is a special class of an unstable non-minimum phase system. Control of a CTS is one of the many challenging problems in this area owing to the complexities in controlling a nonlinear MIMO system. In practice, implementing a controller for a multivariable system is more complex than implementing one for a single input single output (SISO) system because of variations in process dynamics which arise due to changes in operating points and coupling effects. In multivariable control design, often, a decoupler is designed to disintegrate a multiloop control system into a set of single individual loops [1]-[3]. Due to the presence of coupling effects between inputs and outputs in a MIMO system, many well developed closed-loop control methods for a SISO system with time delays can hardly be extended to a MIMO system with delay. To handle coupling effect in a multivariable systems, usually two types of control strategies are adopted i.e. a decentralized controller or a centralized controller. Out of these two, decentralized control structure is widely preferred due to its superior performance,simplicity, and robustness.

    Over the past few decades, many methods have been reported for developing multiloop controller design including the autotuning method [4], [5], independent loop tuning method [3], [4]and detuning method [6]. Although, the detuning process is simple, the loop performance and stability of the closed loop system cannot be adequately defined through this approach. In the independent loop tuning method, the response may be inadequate, as the exact information about controller dynamics of the other loops are not used. Besides these multiple loop control design approaches, many other decoupling approaches have been reported in literature i.e. decentralized control[7]-[9], adaptive neural decentralized control [10], and adaptive decoupling control [11]-[13]. Despite the availability of several other advanced controllers, SISO PID controllers with decouplers are the preferred controllers to use in industry owing to their simple structure and ease in implementation.Despite there are extensive applications of PID controller in industry, it can still be demanding to find an adequate tuning method to determine the appropriate PID controller parameters.A considerable amount of research work has been carried out to find out an efficient tuning formula for the PID controller in[4], [14]-[17]. Robustness is an important aspect of controller design, because the process model used is often an approximation of actual industrial process. Robustness plays a vital role in the control design and is often measured by the gain margin (GM) and phase margin (PM) [18]-[20]. Developing a PID controller to meet GM and PM specifications is an important design approach. However, control design to satisfy GM and PM criteria approach is not new. The solution is usually achieved by numerical methods or graphically by using a Bode plot [21], [22]. In [1], [15], [23], a decentralized PID controller has been proposed based on gain and phase margin (GPM)specifications PID for MIMO system. In general, PI control design based on GM and PM specifications normally leads to four equations with four unknowns ( ωg, ωp,Kcland τi). To solve this problem, numerical methods can be directly employed, or an approximation of arctan function can be adopted to obtain these parameters [24], [25]. In [26], a GPM specification PID controller incorporates an internal model control (IMC) algorithm which significantly simplifies the problem as there is only one tuning parameter.

    Although PID controller tuning always requires a good closed loop performance controller along with robustness,there is a tradeoff between performance and robustness. The maximum closed-loop amplitude ratio (PT) and bandwidth( ωc) are an essential closed-loop performance criteria. In general, a largerPTleads to a faster response. However, it provides a larger overshoot with less robustness. Thus, it should be properly bounded. The bandwidth ( ωc) can be made as large as possible, since it results in a faster closed loop response and lowers the settling time. Even though several well-developed GPM specifications based PID controllers have been proposed [1], [11], [20] and [27], there is no record found on directly imposing constraints on the performance criteriaPTand ωc, along with robustness criteria on GPM PID tuning approach or any experimental validation for the TITO system. Therefore, the motivation of the proposed work is to design an optimal robust PID controller while imposing constraints on both performance and robustness criteria. In this proposed work, constraints onPTbounds peak overshoot and constraints on GPM assure robustness. The settling time is thus minimized as the response is made as fast as possible through maximization of the bandwidth. The contribution of this work lies in the design of an optimal robust PID controller based on nonlinear optimization to precisely regulate the level in a coupled tank system in the face of disturbance as well as uncertainties. Furthermore, the performance of the proposed controller is validated experimentally on a coupled-tank liquid level MIMO system to verify the efficacy of the proposed controller.

    The rest of the paper is organized as follows. In Section II,modeling of the coupled tank system is presented. Decoupler design is presented in Section II. In Section III, a frequency response fitting based model reduction technique is presented to obtain a first order plus dead time (FOPDT) model of each higher order decoupled subsystem. The proposed controller design is described in Section IV. Section V presents robustness analysis of the proposed controller. Simulation as well as experimental results with analysis are presented in Section VI.Finally, conclusions are drawn in Section VII.

    II. MODELING OF A COUPLED TANK SYSTEM

    A schematic diagram of a coupled tank system (CTS) is shown in Fig 1. The objective of controlling the CTS is to maintain the desired liquid level in the individual tanks. By using the mass balance equation and Bernoulli’s principle, the model of the CTS can be written as

    Fig. 1. Schematic diagram of a coupled tank system.

    Some parameters of the CTS such as η andAare not known exactly. Hence, there is a need to obtain an accurate dynamic model of the CTS. So, in this work, FOPDT model is obtained for the CTS from the set of experimental input-output data by employing the system identification technique [28]. A set of input output data of 2500 samples were collected from open loop experiments on the coupled tank system with a sampling time 0.1 s.

    Fig. 2. Model validation based on residual analysis.

    We consider pump control voltage as the input and water level as the output for conducting system identification. By employing the system identification tool box, we obtain a model of the CTS as given in (28). The identified model is validated by using residual analysis. The residual analysis can be performed by evaluating the autocorrelation of output residuals and calculating the cross correlation between the input of the process and to the residual output. The parameters of the identified model in (28) are estimated with prediction error method (PEM). The corresponding results are presented in Figs. 2(a)-2(d). From Figs. 2(a)-2(d), it is observed that the correlation functions remain inside the confidence interval.Hence it is confirmed that, the identified model is suitable for the decoupler design. The obtained identified model is given by

    The above model given in (1) is used for designing the decoupling controllers. In the next section, we present the design of the decoupling controller.

    III. DECOUPLING CONTROLLER DESIGN

    Fig. 3 depicts the block diagram of the control structure of a TITO system. In this structure, two PID controllers are used and the outputs of the PID controllers are fed as inputs to the decouplers. Also, the outputs from the decoupler are used as the two inputs to the CTS. Therefore, in this structure, the control problem is separated into two parts, namely the first part, which deals with decoupling and the second part, which deals with the control of the decoupled loops. Fig. 4 represents the schematic of real-time controller implementation of proposed algorithm. The objective of designing a decoupler is to reduce control loop interactions. In general, two cases are considered for design of decouplers [14] as follows

    Case 1: Off diagonal elements of (1) have no right hand pole (RHP) poles and diagonal elements do not have RHP zeros. Then a decoupling control can be designed as follows.

    Fig. 3. Schematic decentralized control structure for TITO system.

    wherep1(s),p2(s),d12(s) andd21(s) are defined as

    Case 2: Diagonal elements of (1) have no RHP pole and off diagonal elements do not have RHP zeros.

    The decoupled processH(s) with processGp(S) and decouplerD(s) is given by

    The decoupled elementshii(s) are to be controlled by the proposed decentralized controllerCii(s) where i = 1, 2.

    A. Model Order Reduction

    As elementshii(s) of (2) consist of large time delays, it is quite complicated to implement and design a decoupling controller. Therefore, a model reduction approach is considered which is further simplified as (4). In this section, a FOPDT modelPii(s) of each elementshii(s) ofH(s) is obtained as follows [29].

    This reduced order FOPDT model in (4) is derived by employing the frequency response fitting at two points i.e.ω=0 , ω=ωcito obtain the FOPDT model ofhii(s). where ωcidenotes the phase crossover frequency [30].

    By considering the above conditions, the parameters of the FOPDT model given in (4) can be calculated as follows:

    Further, this reduced order decoupled model is used in controller design. The proposed control design approach is presented in the next section.

    IV. PROPOSED OPTIMAL DECENTRALIZED PID CONTROLLER DESIGN

    The proposed PI controller tuning method is designed based on nonlinear optimization. In this work, the overshoot is bounded by the constraint onPTand simultaneously, the robustness is guaranteed by the constraint on GPM. The proposed approach given below is employed to obtain the control laws for the reduced order decoupled subsystem as given in (4), i.e.,P11(s) andP22(s). Often, a FOPDT model is used in industry because it reasonably represents the process gain, time constant, and deadtime of a higher order process.One can write the transfer function of a FOPDT process as follows.

    The transfer function of the PID controller can be written as

    Further, from (6) and (7), the open loop transfer function can be obtained as

    With frequency analysis of each term of (8), the amplitude ratioAROand phase change φOare written as

    where

    From the open loop systemGO, one can obtained the closed loop transfer function as follows:

    The amplitude ratio of closed loop system can be calculated as

    The amplitude ratio of the closed loop system can be obtained from open loop amplitude ratioAROand φO. Further,ωccan be obtained by solving (12)

    The maximum closed loop amplitude ratioPTcan be obtained as

    using definitions of GPM one can write the following set of equations

    Fig. 4. Real-time controller implementation of proposed control.

    where

    By substituting (9) and (10) into (14)-(15), one obtains

    In (16)-(19), there are five unknowns ( ωg, ωp,Kcl, τiand τd).Hence, these equations can not be solved directly. However,the following lemma is defined to formulate an optimization problem by imposing constraints on gain margin, phase margin, and maximum closed loop bandwidth.

    Lemma 1:The optimization problem is defined as follows:

    whereom*is the upper bound of maximum amplitude ratio.

    Remark 1:

    1) The maximum closed loop amplitude ratio is often related to gain and phase margin bounds [30]

    Equation (20) gives the relation betweenPTand GPM bounds and it is identified that ifom*is low, then the actual value ofAmandom*may be higher thanAm*andom*. Thus, a more robust and less aggressive response may be achieved.

    3) In the proposed controller, to solve the unknown variables in (16)-(19) the “fmincon” function is used from the MATLAB optimization tool box.

    V. ROBUSTNESS AND STABILITY ANALYSIS

    In general, there usually exists an unmodeled process dynamics in real-time systems. Therefore, robustness analysis of the designed controller is an important aspect of the control system. The dynamics of the real-time system often has many sources of uncertainties which may cause poor performance and even degrades the stability of the system. Therefore, in this paper to analyze the stability of the decoupled control system multiplicative input and multiplicative output,uncertainties are incorporated into the nominal model system.The corresponding schematic structure of the multiplicative input and multiplicative output uncertainty and itsT-Δ form is shown in Fig. 5 and Fig. 6. By representing the perturbed system in the form ofT-Δ structure for stability analysis, the transfer function from the output to input i.e. ‘u’ to ‘v’ can be written as

    According to the small gain theorem in [30], a perturbed system with uncertainty holds robust stability well if and only if the mentioned robust constraints conditions below are satisfied, i.e.

    Fig. 5. Schematic diagram of multiplicative input uncertainty and itsT-Δ structure.

    Fig. 6. Schematic diagram of multiplicative output uncertainty and itsT-Δ structure.

    But the stability constraints shown in (23) and (24) results in more time needed to evaluate the norm. Hence, to overcome this issue, an equivalent relationship is developed using the small gain theorem and the spectral radius stability[30] which is as follows

    With the aforementioned stability constraints, (25) can be further reformulated as

    VI. RESULTS AND DISCUSSION

    A. Simulation Results

    In this section, simulation results are incorporated to verify the efficacy and performances of the proposed controller. The results of the proposed controller is compared with that of [1].The proposed controller is simulated on the identified model of CTS which is given in Section II. The decoupler for the identified model is obtained by using Case 1 of the decoupling controller design that is described in Section III. The obtained identified model for the CTS that is described in Section II as follows.

    The decouplerD(s) is obtained as

    The resulting diagonal decoupled subsystems are as follows.

    For controller realization,h11(s) andh22(s) should be expressed as FOPDT process. Hence, using (4), the FOPDT models for ofh11(s) andh22(s) are determined as follows.

    Solving (20) forAm= 3 and φm38°, the proposed decentralized PI controller is obtained as

    The suggested ranges of GM and PM are 2-5 and 3 0°-60°,(Am) and PM (φm) to be 3 and 38°respectively. The proposed respectively [31]. To pursue the controller design, we set GM controller parameters are obtained by solving (16)-(20) by using “fmincon” function from MATLAB optimization toolbox. The obtained controller parameters are given in Table I.

    TABLE I PROPOSED CONTROLLER PARAMETER FOR SIMULATION

    To verify the efficacy of the proposed optimal GPM PI controller, its performance is compared to a GPM PI controller reported in [1]. In simulation, first a step change inh1is given att= 0 s andh2demand is set tot= 500 s respectively. Similarly, to illustrate the disturbance rejection performance, a step disturbanced= 0.8 att= 1700 s is introduced to both loops. The corresponding results of both set-point tracking and disturbance rejection are given in Figs. 7-8. From the above simulation results, it is observed that both tanks attain the desired level smoothly without affecting other tank’s levels. However, GPM PI [1] needs more control action to archive steady state and also yields steady state error when there is a load disturbance as compared to the proposed optimal GPM PI controller.Figs. 9-10 show the set-point tracking and disturbance attenuation responses of the proposed optimal GPM PI with a set of differentom*values. From Figs. 9-10, it is observed that, smallerom*values lead to better set-point tracking as well as disturbance attenuation performances.

    Fig. 7. Set-point tracking for nominal system.

    Fig. 8. Disturbance rejection for nominal system.

    Fig. 9. Set point tracking with different om* value.

    Fig. 10. Disturbance rejection with different om* value.

    To study the robustness of the proposed controller in face of parameter uncertainty, the delay, gain, and time constants are changed by ±10%. The corresponding results are presented in Figs. 11-13. To show the practicality and feasibility of the proposed controller, both the performance of the proposed optimal GPM PI and GPM PI [1] are also investigated based on the integral absolute error (IAE), total variance (TV), and RMS tracking error (eti). These are presented in Table II.Apart from this, the robust stability of the proposed controller is analyzed while considering the multiplicative input and multiplicative output uncertainties in Section V. The corresponding magnitude plot of the spectral radius of(26)-(27) are shown in Fig. 14. From Fig. 14, it is observed that the proposed controller also provides robust stability.

    Fig. 11. Set-point tracking for perturbed system with 10% parameter uncertainty.

    Fig. 12. Set-point tracking for perturbed system with -10% parameter uncertainty.

    Fig. 13. Simulation results of disturbance rejection for perturbed system.

    Fig. 14. Magnitude plots of spectral radius.

    TABLE II PERFORMANCE ANALYSIS

    B. Experimental Results

    In order to demonstrate practicality of the proposed decoupled optimal GPM PI controller described in Section IV,it is validated in real-time. Fig. 15 depicts the sketch of the coupled tank system that is used for real-time implementation.In this set-up, each tank is fitted with an outlet pipe to transmit overflow water into the storage tank. Also, there are two submersible pumps to drive water from bottom to top of the tanks. An air pressure sensor is attached to the base of each container to measure the water level of the corresponding tank. A PC with an Advantech card and MATLAB/SIMULINK is used to implement the algorithms. The PSUPA (Power Supply Power Amplify) unit amplifies the water pressurelevel signals and passes them as an analogue signal to the PCI1711 DAQ card. The pumps control signal can be sent from the PC through the DAQ (PCI1711) card and PSUPA unit. The control signals, which are between 0 V-5 V, are transferred to the PSUPA unit where these are transformed into 24 V PWM signals to drive the pumps. The water level information is also sent to the PC via the same PSUPA unit.The model of the CTS is obtained based on the set of experimental input-output data of 2000 samples. Each sample is then collected from the open loop experiment on the CTS with a sampling time of 0.1 s by using the System Identification Tool Box in MATLAB/SIMULINK. The response of model output and actual output using pump 1(pump 2 off) is shown in Fig. 16(a) for tank 1 and Fig. 16(b)for tank 2. Similarly, using pump 2 (pump 1 off) is shown in Fig. 17(a) for tank 2 and Fig. 17(b) for tank 1. Figs. 18-19 exhibit the set-point tracking and disturbance rejection results of the proposed controller. From Fig. 19, it is observed that,the proposed controller exhibits good disturbance rejection capabilities when a load disturbance is applied at 1200 s by opening of valve 4 and valve 5 for 100 s as shown in Fig. 15.Fig. 20 depicts the efficacy of the proposed controller in the face of uncertainty, which is introduced by 50 % opening of the coupling valve whenever both tanks are at their steady state. The corresponding control signals are shown in Figs.21-22. In practice, the opening of the coupling valve not only affects the level interaction, but also alters the operating point of the system. In general, the changes introduces parametric uncertainty. However, from Fig. 20 it is observed that, the proposed controller provides appropriate control action to reduce the interaction for achieving the desired level in the face of uncertainty.

    In the proposed approach, to determine the unknown variables (ωg, ωp,Kcl, τi, τd), a nonlinear optimization margin, phase margin, and maximum amplitude ratio (om*).problem is formulated, imposing constraints on the gain This nonlinear optimization problem is solved by using Lemma 1 with the “fmincon” function in MATLAB. In the proposed controller, the value ofom*influences the dynamic response (both set-point tracking and disturbance rejection) of the closed loop system as seen from Figs. 9-10, in which, the effect on the closed loop performance of the system with variation ofom*are studied. Thus, the solution of the above optimization problem with the “fmincon” function, shows that the value ofom*plays an essential role in obtaining appropriate parameters ( ωg, ωp,Kcl, τi, τd) of the PI controller. From the obtained results shown in Figs. 9-10, it is clear that selection of an appropriateom*is a trade-off between the closed loop performance and robustness. From a computational view point, it is observed that the proposed controller involves less computational complexity for the choice of an appropriate value forom*. However, as compared to the existing GPM PI controller [1], the proposed optimal GPM PI controller shows a trade-off between good closed loop dynamic performance and robustness.

    VII. CONCLUSIONS

    In this paper, a robust optimal decentralized PID controller is proposed based on nonlinear constraint optimization to simultaneously meet both robustness and loop performance requirements. The maximum closed loop amplitude ratio (PT)and bandwidth ( ωc) are considered as the closed-loop performance criteria. On the other hand, GPM serves as a robustness criterion. Furthermore, the unique advantage of the proposed controller, is the flexibility that is brought by imposing constraints on the maximum closed-loop amplitude.With proper tuning of the bound onPT, one can easily obtain the parameters of the PID to achieve good set-point tracking as well as a good load disturbance attenuation response. To investigate the performance of the proposed controller, the simulation results and different performance indices are compared in Table II with the work reported in [1]. To illustrate the robustness of the proposed controller,±10%parametric uncertainty is added to gain and time delays. Also,robust stability is analysed while considering multiplicative input and multiplicative output uncertainties. The obtained results reveal that the proposed controller exhibits improved performance with minimal interaction and also provides robust performance in face of uncertainty as compared to [1].To represent the practical applicability of the proposed controller, an experiment is conducted with a real-time coupled-tank liquid level system in the laboratory. From the obtained simulation and experimental results, it is determined that the proposed controller provides robust stability as well as excellent loop performance.

    Fig. 15. Sketch of the experimental set-up.

    Fig. 16. The step response of model from input 1 to output 1 and output 2.

    Fig. 17. The step response of model from input 2 to output 2 and output 1.

    Fig. 18. Experimental response of set-point tracking.

    Fig. 19. Experimental response of disturbance rejection.

    Fig. 20. Experimental response of both tanks for 50% opening of coupling valve.

    Fig. 21. Experimental response of control input of both tanks for 50% opening of coupling valve.

    Fig. 22. Experimental response of control input of both tanks for set-point tracking.

    APPENDIX

    The open loop transfer function of a FOPDT system that is given in (8) is rewritten here

    In general the solution and closed-loop analysis of a third order system is complex. Therefore, the Routh stability criterion model reduction method [32] is adopted for simplification of the obtained model (34). The reduced order transfer function is determined directly from the elements of the higher order denominator and numerator of the Routh Array. The obtained reduced model is yielded as follows

    where κ1and κ2are defined as

    Now considering the characteristics equation by substituting(35) in (36) gives

    where

    Equation (37) is in the form ofs2+2ξωns+ωn2. After simplification of (37), one can write the following set of equations

    After model order reduction, the model given in (34) is approximated as a second order model. Therefore, one can write the closed loop resonant frequency, resonant peak,maximum overshoot, and bandwidth as follows:

    The damping ratio and phase margin relation can be written as

    From (40), it is clear that, bothMrandMpare the functions of ξ. Also from Fig. 23, it is seen that as ξ increase, bothMrandMpdecrease. Fig. 24 depicts the frequency domain performance specification of a second order system. Also from Fig. 23 it is observed that when ξ decreases, bothMrandMpincrease. Thus, whenMrincreases, the range of frequency with gain becomes more than 0.707, i.e., the bandwidth also increases.

    Fig. 23. Response of ξ vs bandwith, Mp and Mr.

    Fig. 24. Frequency domain performance specification.

    免费高清在线观看视频在线观看| 极品教师在线视频| 欧美变态另类bdsm刘玥| 欧美精品一区二区大全| 午夜激情久久久久久久| 亚洲欧美中文字幕日韩二区| 国产一区有黄有色的免费视频| 欧美激情国产日韩精品一区| 亚洲天堂av无毛| 一边亲一边摸免费视频| 人妻少妇偷人精品九色| 亚洲精品国产av蜜桃| 欧美人与善性xxx| 国内少妇人妻偷人精品xxx网站| 久久久久视频综合| 色婷婷久久久亚洲欧美| 免费人成在线观看视频色| 国产有黄有色有爽视频| 国产男人的电影天堂91| 在线免费十八禁| 精品亚洲成国产av| 久久久午夜欧美精品| 又粗又硬又长又爽又黄的视频| 一级毛片我不卡| 亚洲精品乱码久久久v下载方式| 成人高潮视频无遮挡免费网站| 国产黄频视频在线观看| av国产精品久久久久影院| 热re99久久精品国产66热6| 国产在线男女| 免费高清在线观看视频在线观看| 国产黄片视频在线免费观看| 国产爽快片一区二区三区| 看免费成人av毛片| 亚洲怡红院男人天堂| 久久6这里有精品| 美女视频免费永久观看网站| 女性生殖器流出的白浆| 一区二区三区精品91| 搡老乐熟女国产| 午夜激情福利司机影院| 亚洲国产精品999| 久久久午夜欧美精品| 成人特级av手机在线观看| a 毛片基地| 联通29元200g的流量卡| freevideosex欧美| 在线观看美女被高潮喷水网站| 国产精品99久久99久久久不卡 | 大香蕉97超碰在线| 国产淫语在线视频| 日韩欧美 国产精品| 在线观看三级黄色| 国产淫片久久久久久久久| 亚洲人成网站高清观看| 两个人的视频大全免费| 我要看日韩黄色一级片| 91精品国产九色| 伊人久久精品亚洲午夜| 久久久成人免费电影| 欧美成人a在线观看| 国产91av在线免费观看| 大陆偷拍与自拍| 色婷婷av一区二区三区视频| 久久久久网色| 另类亚洲欧美激情| 欧美日韩视频精品一区| 能在线免费看毛片的网站| 久久久久久久久久久丰满| 黄色日韩在线| 日韩制服骚丝袜av| 五月天丁香电影| 联通29元200g的流量卡| 亚洲精品乱码久久久久久按摩| 777米奇影视久久| 高清欧美精品videossex| 久久国产精品男人的天堂亚洲 | 美女高潮的动态| 麻豆成人av视频| 日本av免费视频播放| 人人妻人人添人人爽欧美一区卜 | 日韩av在线免费看完整版不卡| 国产精品三级大全| 国产一级毛片在线| 精品国产露脸久久av麻豆| 菩萨蛮人人尽说江南好唐韦庄| 丰满少妇做爰视频| 免费黄频网站在线观看国产| 亚洲精品乱码久久久久久按摩| 噜噜噜噜噜久久久久久91| 韩国av在线不卡| 人人妻人人爽人人添夜夜欢视频 | 丰满人妻一区二区三区视频av| 最近手机中文字幕大全| 亚洲人成网站在线观看播放| 少妇人妻一区二区三区视频| 王馨瑶露胸无遮挡在线观看| 精品亚洲乱码少妇综合久久| 1000部很黄的大片| 激情 狠狠 欧美| 成人二区视频| 国产熟女欧美一区二区| 黑人高潮一二区| 老司机影院成人| 赤兔流量卡办理| 免费人妻精品一区二区三区视频| 亚洲成人手机| 大又大粗又爽又黄少妇毛片口| 男女啪啪激烈高潮av片| 91在线精品国自产拍蜜月| 蜜桃久久精品国产亚洲av| 免费av不卡在线播放| 一本色道久久久久久精品综合| 五月天丁香电影| 赤兔流量卡办理| 街头女战士在线观看网站| 黄色怎么调成土黄色| 成年av动漫网址| 全区人妻精品视频| 亚洲国产精品一区三区| 中文字幕亚洲精品专区| 高清黄色对白视频在线免费看 | 人妻少妇偷人精品九色| www.av在线官网国产| 久久99热这里只有精品18| 中文字幕精品免费在线观看视频 | 一区二区三区四区激情视频| 久久精品国产a三级三级三级| 国产有黄有色有爽视频| 18禁在线播放成人免费| 又爽又黄a免费视频| 1000部很黄的大片| 久久6这里有精品| 国产 一区 欧美 日韩| 我的老师免费观看完整版| 激情 狠狠 欧美| 中文字幕亚洲精品专区| 欧美精品人与动牲交sv欧美| 国产高清有码在线观看视频| 免费少妇av软件| 2022亚洲国产成人精品| 大片免费播放器 马上看| 成人美女网站在线观看视频| 亚洲美女黄色视频免费看| 大片免费播放器 马上看| 国产免费又黄又爽又色| 亚洲成人中文字幕在线播放| 国产精品熟女久久久久浪| 国产午夜精品一二区理论片| 免费播放大片免费观看视频在线观看| 中文乱码字字幕精品一区二区三区| av视频免费观看在线观看| 亚洲久久久国产精品| 26uuu在线亚洲综合色| 日韩免费高清中文字幕av| 国产在线男女| 欧美bdsm另类| 久久久久国产精品人妻一区二区| 黑丝袜美女国产一区| 精品久久国产蜜桃| 亚洲人成网站在线播| 精品熟女少妇av免费看| 久久影院123| 国产深夜福利视频在线观看| 这个男人来自地球电影免费观看 | 搡老乐熟女国产| 精品99又大又爽又粗少妇毛片| 精品熟女少妇av免费看| 免费观看性生交大片5| videos熟女内射| 日韩av在线免费看完整版不卡| av在线播放精品| 亚洲国产精品国产精品| 男女边摸边吃奶| 啦啦啦视频在线资源免费观看| 国产精品麻豆人妻色哟哟久久| 亚洲精品日韩在线中文字幕| 日韩欧美精品免费久久| 亚洲无线观看免费| 国产亚洲91精品色在线| 亚洲精品国产av蜜桃| 在线看a的网站| 美女xxoo啪啪120秒动态图| 蜜桃久久精品国产亚洲av| 国产精品一区二区性色av| 99九九线精品视频在线观看视频| 国产色婷婷99| 国产精品av视频在线免费观看| 女性被躁到高潮视频| 国产伦在线观看视频一区| 99热网站在线观看| 91精品国产九色| 午夜福利在线观看免费完整高清在| 亚洲精品国产色婷婷电影| 国产欧美日韩一区二区三区在线 | 亚洲av二区三区四区| 国产成人精品婷婷| 亚洲不卡免费看| 亚洲精品aⅴ在线观看| 亚洲精品久久午夜乱码| 国产免费又黄又爽又色| 精品久久久久久电影网| 美女视频免费永久观看网站| 中文天堂在线官网| 欧美xxxx黑人xx丫x性爽| 男女无遮挡免费网站观看| 婷婷色综合大香蕉| 亚洲欧美日韩另类电影网站 | 久久影院123| 尤物成人国产欧美一区二区三区| 哪个播放器可以免费观看大片| 亚洲欧美清纯卡通| 精品酒店卫生间| 国产精品三级大全| 少妇精品久久久久久久| 26uuu在线亚洲综合色| 在线免费十八禁| 五月开心婷婷网| 色哟哟·www| 在线亚洲精品国产二区图片欧美 | 精品国产一区二区三区久久久樱花 | 丰满乱子伦码专区| 国产成人精品一,二区| 亚洲熟女精品中文字幕| 久久人人爽人人片av| 人体艺术视频欧美日本| 少妇猛男粗大的猛烈进出视频| 99久久精品一区二区三区| 日韩不卡一区二区三区视频在线| 国产69精品久久久久777片| 亚洲国产高清在线一区二区三| 国产精品一区二区在线不卡| 男女无遮挡免费网站观看| 99热国产这里只有精品6| 最近2019中文字幕mv第一页| 午夜福利在线观看免费完整高清在| 久久99热这里只频精品6学生| 夜夜爽夜夜爽视频| 亚洲av男天堂| 午夜精品国产一区二区电影| 亚洲精品日韩在线中文字幕| av女优亚洲男人天堂| 性高湖久久久久久久久免费观看| 青春草亚洲视频在线观看| 中文在线观看免费www的网站| 91aial.com中文字幕在线观看| 亚洲av免费高清在线观看| 免费大片黄手机在线观看| 97在线视频观看| 国产在线一区二区三区精| av在线观看视频网站免费| 毛片女人毛片| av福利片在线观看| 天美传媒精品一区二区| 亚洲精品一二三| 黑人高潮一二区| 久久久久久久国产电影| 国产午夜精品久久久久久一区二区三区| 高清欧美精品videossex| 久久鲁丝午夜福利片| 日本猛色少妇xxxxx猛交久久| 丝瓜视频免费看黄片| 久久99精品国语久久久| 狂野欧美白嫩少妇大欣赏| 黑丝袜美女国产一区| 色视频在线一区二区三区| 久久精品熟女亚洲av麻豆精品| 国产伦精品一区二区三区视频9| 秋霞在线观看毛片| 欧美三级亚洲精品| 国语对白做爰xxxⅹ性视频网站| 欧美+日韩+精品| 一边亲一边摸免费视频| 成人美女网站在线观看视频| 国产成人91sexporn| 国产无遮挡羞羞视频在线观看| 午夜福利高清视频| 一本一本综合久久| 精品国产一区二区三区久久久樱花 | 免费在线观看成人毛片| av网站免费在线观看视频| av卡一久久| 91狼人影院| 日本猛色少妇xxxxx猛交久久| av在线app专区| 国产一区二区在线观看日韩| 亚洲av综合色区一区| 身体一侧抽搐| 欧美日韩国产mv在线观看视频 | 春色校园在线视频观看| 久久久久久久久久人人人人人人| 91精品国产国语对白视频| 欧美老熟妇乱子伦牲交| 日韩强制内射视频| 中文字幕久久专区| 日韩中文字幕视频在线看片 | 精品久久久久久久久av| 老师上课跳d突然被开到最大视频| 欧美日韩精品成人综合77777| 韩国av在线不卡| 国产女主播在线喷水免费视频网站| 亚洲国产最新在线播放| 不卡视频在线观看欧美| 日韩强制内射视频| a级一级毛片免费在线观看| 色吧在线观看| 久久久a久久爽久久v久久| 亚洲av成人精品一二三区| 啦啦啦啦在线视频资源| 免费高清在线观看视频在线观看| 2021少妇久久久久久久久久久| 亚洲国产日韩一区二区| 欧美性感艳星| 高清视频免费观看一区二区| 久久人人爽av亚洲精品天堂 | 一区二区三区四区激情视频| 不卡视频在线观看欧美| 久久精品久久精品一区二区三区| 内射极品少妇av片p| 丰满乱子伦码专区| 人人妻人人看人人澡| 国产亚洲91精品色在线| 欧美97在线视频| 日日撸夜夜添| 精品亚洲乱码少妇综合久久| 亚洲国产av新网站| 亚洲图色成人| 日韩制服骚丝袜av| 亚洲图色成人| 亚洲伊人久久精品综合| 国产乱人偷精品视频| 成人黄色视频免费在线看| 国产精品一二三区在线看| 亚洲av日韩在线播放| 国产在视频线精品| 国产精品久久久久久久电影| 日本一二三区视频观看| 女人久久www免费人成看片| 免费观看性生交大片5| 麻豆成人av视频| 丰满迷人的少妇在线观看| 日韩在线高清观看一区二区三区| 免费av中文字幕在线| 日本欧美视频一区| 一本色道久久久久久精品综合| 精品人妻一区二区三区麻豆| 我的女老师完整版在线观看| 一个人免费看片子| 久久人人爽人人爽人人片va| 精品一品国产午夜福利视频| 久久久国产一区二区| 欧美97在线视频| 国产男人的电影天堂91| 日韩一区二区视频免费看| 一级片'在线观看视频| 国产色婷婷99| 国国产精品蜜臀av免费| h视频一区二区三区| 午夜激情久久久久久久| av专区在线播放| 一级片'在线观看视频| 黄色欧美视频在线观看| 夫妻午夜视频| 欧美日本视频| 午夜免费鲁丝| 国产精品.久久久| 在线观看免费高清a一片| 国产成人精品福利久久| 成人无遮挡网站| av福利片在线观看| 99久久精品国产国产毛片| 欧美三级亚洲精品| 亚洲国产精品成人久久小说| 男人舔奶头视频| 免费黄色在线免费观看| 亚洲美女黄色视频免费看| 免费看av在线观看网站| 3wmmmm亚洲av在线观看| 伊人久久国产一区二区| 18禁裸乳无遮挡动漫免费视频| 精品亚洲成a人片在线观看 | 久久久欧美国产精品| 男女国产视频网站| 国产伦在线观看视频一区| 日韩成人av中文字幕在线观看| 一个人看视频在线观看www免费| 综合色丁香网| 网址你懂的国产日韩在线| 男男h啪啪无遮挡| 国产黄片美女视频| 能在线免费看毛片的网站| 国产有黄有色有爽视频| 在线免费十八禁| 日韩一本色道免费dvd| 少妇 在线观看| 婷婷色麻豆天堂久久| 韩国av在线不卡| 日本与韩国留学比较| 青春草国产在线视频| 国产午夜精品久久久久久一区二区三区| 国产精品久久久久久精品古装| 国产 一区精品| 美女视频免费永久观看网站| 国产伦精品一区二区三区视频9| 尤物成人国产欧美一区二区三区| 777米奇影视久久| av专区在线播放| 日韩强制内射视频| 免费看光身美女| 在线观看免费高清a一片| 亚洲成人av在线免费| 我要看黄色一级片免费的| 我的老师免费观看完整版| 国产成人aa在线观看| 欧美亚洲 丝袜 人妻 在线| 亚洲精品中文字幕在线视频 | 一级毛片 在线播放| 亚洲av成人精品一区久久| 大片电影免费在线观看免费| 在线播放无遮挡| 校园人妻丝袜中文字幕| 夫妻性生交免费视频一级片| 男人和女人高潮做爰伦理| 国产女主播在线喷水免费视频网站| 最近2019中文字幕mv第一页| 国产 一区精品| 九九爱精品视频在线观看| 日韩制服骚丝袜av| 日韩电影二区| 99九九线精品视频在线观看视频| 国产精品.久久久| 国语对白做爰xxxⅹ性视频网站| 日本黄大片高清| 亚洲欧美一区二区三区国产| 黄片无遮挡物在线观看| 久久精品国产a三级三级三级| 亚洲精品乱久久久久久| 午夜免费男女啪啪视频观看| av在线播放精品| 欧美日韩综合久久久久久| 大香蕉97超碰在线| 汤姆久久久久久久影院中文字幕| 日韩成人伦理影院| 国产亚洲午夜精品一区二区久久| 亚洲无线观看免费| 美女内射精品一级片tv| 国产免费一级a男人的天堂| 国产中年淑女户外野战色| 哪个播放器可以免费观看大片| 国产深夜福利视频在线观看| 国产亚洲91精品色在线| 国产毛片在线视频| 亚洲,一卡二卡三卡| 色网站视频免费| 在线免费十八禁| 自拍偷自拍亚洲精品老妇| 久久久久久人妻| 久久久久视频综合| 国产精品不卡视频一区二区| 国产精品久久久久久av不卡| 成人毛片60女人毛片免费| 国产精品一区www在线观看| 免费人妻精品一区二区三区视频| 久久精品国产亚洲av天美| 亚洲自偷自拍三级| 熟女av电影| 久久6这里有精品| 国产精品久久久久久久电影| 99热全是精品| 国产伦精品一区二区三区视频9| 国产日韩欧美在线精品| 国产视频首页在线观看| 少妇高潮的动态图| 久久久久网色| 亚洲精品日韩av片在线观看| av在线老鸭窝| 亚洲国产精品一区三区| 国产男女超爽视频在线观看| 熟女人妻精品中文字幕| 五月伊人婷婷丁香| 日本欧美视频一区| 欧美日韩视频高清一区二区三区二| 久久久久国产网址| 亚洲四区av| 自拍偷自拍亚洲精品老妇| 国产深夜福利视频在线观看| 色哟哟·www| 伊人久久精品亚洲午夜| 国产亚洲精品久久久com| 亚洲精品乱码久久久久久按摩| 18禁动态无遮挡网站| 黑人高潮一二区| 久久亚洲国产成人精品v| 国产精品久久久久久av不卡| 久久久久久九九精品二区国产| 国产精品久久久久久久久免| 爱豆传媒免费全集在线观看| 91狼人影院| 成人无遮挡网站| 日韩av在线免费看完整版不卡| 中文字幕制服av| 啦啦啦在线观看免费高清www| 18+在线观看网站| 91精品国产国语对白视频| 1000部很黄的大片| 极品教师在线视频| 啦啦啦在线观看免费高清www| 伦精品一区二区三区| 一区二区三区乱码不卡18| 久久久亚洲精品成人影院| 亚洲国产成人一精品久久久| 国产亚洲欧美精品永久| 在线精品无人区一区二区三 | 日韩精品有码人妻一区| 大又大粗又爽又黄少妇毛片口| 黄色怎么调成土黄色| 亚洲国产高清在线一区二区三| 青春草国产在线视频| 精品一区二区免费观看| 国产在线男女| 18禁在线无遮挡免费观看视频| 久久久亚洲精品成人影院| 精品99又大又爽又粗少妇毛片| 男女无遮挡免费网站观看| 超碰av人人做人人爽久久| 全区人妻精品视频| 妹子高潮喷水视频| 最近最新中文字幕免费大全7| 熟女av电影| 一区二区三区乱码不卡18| 婷婷色麻豆天堂久久| 久久精品国产鲁丝片午夜精品| 97在线视频观看| 国产黄片美女视频| 午夜激情福利司机影院| 日韩制服骚丝袜av| kizo精华| 国产精品欧美亚洲77777| 免费大片黄手机在线观看| 亚洲国产色片| 26uuu在线亚洲综合色| av专区在线播放| 午夜福利视频精品| 午夜日本视频在线| 尤物成人国产欧美一区二区三区| 中国国产av一级| 色网站视频免费| 国产精品一二三区在线看| 日韩大片免费观看网站| 99国产精品免费福利视频| 大香蕉97超碰在线| 观看av在线不卡| 在线观看免费视频网站a站| 午夜福利视频精品| 亚洲精品视频女| 精品人妻一区二区三区麻豆| 久久综合国产亚洲精品| 你懂的网址亚洲精品在线观看| 亚洲精品成人av观看孕妇| 色婷婷av一区二区三区视频| 人妻一区二区av| 我要看黄色一级片免费的| 哪个播放器可以免费观看大片| 亚洲三级黄色毛片| 亚洲欧美日韩卡通动漫| 久久精品人妻少妇| 日日啪夜夜爽| 18禁在线播放成人免费| 大话2 男鬼变身卡| 久久久久久人妻| 成人影院久久| 欧美日韩一区二区视频在线观看视频在线| 超碰97精品在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 日韩伦理黄色片| 少妇的逼水好多| 成人无遮挡网站| 国产午夜精品久久久久久一区二区三区| 国产成人a区在线观看| 久久青草综合色| 久久久久久九九精品二区国产| 直男gayav资源| 激情 狠狠 欧美| 亚洲欧美精品自产自拍| 亚洲人成网站在线播| 狂野欧美激情性xxxx在线观看| 夜夜看夜夜爽夜夜摸| 国产精品蜜桃在线观看| 国产精品免费大片| 色网站视频免费| 免费看日本二区| 夫妻午夜视频| 精品国产露脸久久av麻豆| 蜜桃亚洲精品一区二区三区| 夫妻午夜视频| 一级毛片黄色毛片免费观看视频| 国产精品蜜桃在线观看| 亚洲精品国产成人久久av| 肉色欧美久久久久久久蜜桃| 国产v大片淫在线免费观看| 免费av中文字幕在线| 全区人妻精品视频| 国产欧美日韩一区二区三区在线 | 一级片'在线观看视频| 黄色日韩在线| 欧美3d第一页| 亚洲欧洲日产国产| 精品酒店卫生间| 久久久色成人| 只有这里有精品99| 99久久人妻综合| 成人亚洲欧美一区二区av| 精品一区在线观看国产| 小蜜桃在线观看免费完整版高清| 久久97久久精品|