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

    熟妇人妻久久中文字幕3abv| 久久精品国产清高在天天线| 小说图片视频综合网站| 国产伦精品一区二区三区四那| 久久久久久久久久久丰满| 日韩欧美精品免费久久| 一个人免费在线观看电影| 久久精品影院6| 亚洲美女视频黄频| 深夜a级毛片| 国产高潮美女av| 不卡一级毛片| 九九在线视频观看精品| 亚洲av五月六月丁香网| 又黄又爽又免费观看的视频| 久久久欧美国产精品| 99热这里只有是精品50| 亚洲不卡免费看| 国产69精品久久久久777片| 97热精品久久久久久| 在线观看午夜福利视频| av在线播放精品| 成人漫画全彩无遮挡| 变态另类丝袜制服| 日韩欧美三级三区| 免费高清视频大片| 国产一区二区在线av高清观看| 亚洲精品粉嫩美女一区| 亚洲av美国av| 久久精品久久久久久噜噜老黄 | 国产成人91sexporn| eeuss影院久久| 美女免费视频网站| 成人特级黄色片久久久久久久| 不卡视频在线观看欧美| 91狼人影院| 少妇熟女aⅴ在线视频| 亚洲av五月六月丁香网| 三级经典国产精品| 搡老熟女国产l中国老女人| 亚洲丝袜综合中文字幕| 啦啦啦观看免费观看视频高清| 国产精品永久免费网站| 中出人妻视频一区二区| 99热这里只有是精品在线观看| 国产日本99.免费观看| 一进一出抽搐动态| 99riav亚洲国产免费| 中国美女看黄片| 国产精品久久视频播放| 午夜亚洲福利在线播放| 一个人免费在线观看电影| 免费搜索国产男女视频| 校园人妻丝袜中文字幕| 男女边吃奶边做爰视频| 成人漫画全彩无遮挡| 欧美在线一区亚洲| 日韩欧美三级三区| 男女啪啪激烈高潮av片| 日韩av不卡免费在线播放| 欧美+亚洲+日韩+国产| 日日啪夜夜撸| 自拍偷自拍亚洲精品老妇| 亚洲精品在线观看二区| 日日啪夜夜撸| 国产午夜精品久久久久久一区二区三区 | 一进一出好大好爽视频| 日本色播在线视频| 国产精品无大码| 午夜福利在线观看免费完整高清在 | 欧美zozozo另类| 日本精品一区二区三区蜜桃| 搡老岳熟女国产| 亚洲av.av天堂| 亚洲国产精品sss在线观看| 亚洲国产精品国产精品| 色av中文字幕| 国产乱人偷精品视频| 亚洲精品乱码久久久v下载方式| 精品无人区乱码1区二区| 少妇熟女欧美另类| 欧美丝袜亚洲另类| 少妇熟女欧美另类| 女人被狂操c到高潮| 女的被弄到高潮叫床怎么办| 欧美极品一区二区三区四区| 一级av片app| 欧美3d第一页| 亚洲欧美精品综合久久99| 99热这里只有精品一区| 99热全是精品| 中文字幕精品亚洲无线码一区| 97超级碰碰碰精品色视频在线观看| 欧美中文日本在线观看视频| 国产精品伦人一区二区| 国产精品久久电影中文字幕| 国产精品一区二区免费欧美| 十八禁国产超污无遮挡网站| 黄色欧美视频在线观看| 欧美一级a爱片免费观看看| 亚洲中文日韩欧美视频| 亚洲精华国产精华液的使用体验 | 国产亚洲精品久久久com| 久久久久国内视频| 日韩欧美精品免费久久| 亚洲国产精品成人久久小说 | 97超级碰碰碰精品色视频在线观看| 禁无遮挡网站| 精品欧美国产一区二区三| 国产精品嫩草影院av在线观看| 国产成人91sexporn| 草草在线视频免费看| 我要搜黄色片| 成年免费大片在线观看| 赤兔流量卡办理| 久久久久久九九精品二区国产| 天堂动漫精品| 一个人观看的视频www高清免费观看| 狂野欧美白嫩少妇大欣赏| 色尼玛亚洲综合影院| 大香蕉久久网| 亚洲美女黄片视频| 一个人免费在线观看电影| 少妇熟女欧美另类| 亚洲精华国产精华液的使用体验 | av在线蜜桃| 色在线成人网| 欧美精品国产亚洲| 亚洲高清免费不卡视频| 麻豆国产av国片精品| 亚洲成人精品中文字幕电影| 免费观看精品视频网站| 联通29元200g的流量卡| 欧美日本视频| 免费看光身美女| 不卡一级毛片| 又粗又爽又猛毛片免费看| 蜜臀久久99精品久久宅男| 97碰自拍视频| 高清午夜精品一区二区三区 | 极品教师在线视频| 精品久久久久久久久久久久久| 国产精品一区二区免费欧美| 欧美激情国产日韩精品一区| 内射极品少妇av片p| 高清午夜精品一区二区三区 | 你懂的网址亚洲精品在线观看 | 亚洲高清免费不卡视频| 久久午夜亚洲精品久久| 国产精品久久久久久久久免| 国产精品一区二区免费欧美| 久久久久国产网址| 精品熟女少妇av免费看| 国产av在哪里看| 身体一侧抽搐| 91午夜精品亚洲一区二区三区| 一级a爱片免费观看的视频| 亚洲18禁久久av| 超碰av人人做人人爽久久| 99国产精品一区二区蜜桃av| 亚洲国产色片| 亚洲美女搞黄在线观看 | 国产不卡一卡二| 免费不卡的大黄色大毛片视频在线观看 | 免费看光身美女| 亚洲一区二区三区色噜噜| 国产女主播在线喷水免费视频网站 | 欧美日韩国产亚洲二区| 亚洲av五月六月丁香网| 欧美日韩一区二区视频在线观看视频在线 | 久久久久久久久久黄片| 舔av片在线| 毛片女人毛片| 久久精品国产99精品国产亚洲性色| 天堂影院成人在线观看| 国产精品野战在线观看| 国产白丝娇喘喷水9色精品| 久久久久久九九精品二区国产| 国产精品av视频在线免费观看| 成年女人永久免费观看视频| 一个人免费在线观看电影| 成人无遮挡网站| 人妻少妇偷人精品九色| 午夜福利成人在线免费观看| 国产探花在线观看一区二区| 看免费成人av毛片| 1000部很黄的大片| 色噜噜av男人的天堂激情| 精品乱码久久久久久99久播| 嫩草影院新地址| 成人av在线播放网站| 别揉我奶头 嗯啊视频| 看黄色毛片网站| av福利片在线观看| 亚洲欧美日韩无卡精品| 哪里可以看免费的av片| 日日干狠狠操夜夜爽| 亚洲成a人片在线一区二区| 国产成人a区在线观看| 成人高潮视频无遮挡免费网站| 国产高清不卡午夜福利| 97碰自拍视频| 91麻豆精品激情在线观看国产| 国产精品日韩av在线免费观看| 免费人成视频x8x8入口观看| 日韩亚洲欧美综合| 直男gayav资源| 日韩国内少妇激情av| 欧美+亚洲+日韩+国产| 亚洲av.av天堂| 精品久久久久久久末码| 午夜精品国产一区二区电影 | 在线观看一区二区三区| 淫秽高清视频在线观看| 亚洲国产精品成人久久小说 | 亚洲av成人精品一区久久| 欧美最新免费一区二区三区| 精品国产三级普通话版| 国产精品精品国产色婷婷| 久久久午夜欧美精品| 丰满的人妻完整版| 1024手机看黄色片| 久久久久久久久中文| 精品久久久久久久久亚洲| 高清毛片免费观看视频网站| 99热这里只有是精品50| 国内少妇人妻偷人精品xxx网站| 国产精品人妻久久久影院| 一进一出抽搐动态| 欧美绝顶高潮抽搐喷水| 欧美极品一区二区三区四区| 中文字幕免费在线视频6| 在线观看66精品国产| 亚洲乱码一区二区免费版| 一夜夜www| 亚洲最大成人手机在线| 午夜日韩欧美国产| 久久热精品热| 亚洲激情五月婷婷啪啪| 偷拍熟女少妇极品色| av女优亚洲男人天堂| 国产精品亚洲美女久久久| 精品人妻视频免费看| 欧美性感艳星| 午夜福利成人在线免费观看| 日韩精品中文字幕看吧| 日本与韩国留学比较| 国产亚洲91精品色在线| 给我免费播放毛片高清在线观看| 青春草视频在线免费观看| 搡老熟女国产l中国老女人| 亚洲精品色激情综合| av在线天堂中文字幕| 久久久久性生活片| 午夜免费男女啪啪视频观看 | 日本免费a在线| 精品人妻偷拍中文字幕| 男人舔奶头视频| 波多野结衣高清作品| 在现免费观看毛片| 又爽又黄无遮挡网站| 国产精品人妻久久久久久| 日韩av不卡免费在线播放| 午夜a级毛片| 色播亚洲综合网| 99久国产av精品| 欧美一级a爱片免费观看看| 欧美不卡视频在线免费观看| 麻豆久久精品国产亚洲av| 麻豆一二三区av精品| 亚洲国产色片| 又爽又黄a免费视频| 97超碰精品成人国产| 亚洲国产精品合色在线| 国产不卡一卡二| 在线免费观看的www视频| 听说在线观看完整版免费高清| 免费在线观看成人毛片| 国产在视频线在精品| 小蜜桃在线观看免费完整版高清| 国产探花极品一区二区| 国产精品久久久久久亚洲av鲁大| 精品一区二区免费观看| 亚洲性夜色夜夜综合| 麻豆乱淫一区二区| 人妻久久中文字幕网| 丝袜喷水一区| 熟妇人妻久久中文字幕3abv| 亚洲av中文av极速乱| 久久久久久久久大av| 一区福利在线观看| 久久精品久久久久久噜噜老黄 | 国产蜜桃级精品一区二区三区| 国产精品av视频在线免费观看| 欧美一区二区国产精品久久精品| 最后的刺客免费高清国语| 精品一区二区三区av网在线观看| 91在线观看av| 亚洲中文字幕一区二区三区有码在线看| 精品一区二区免费观看| 久久久精品欧美日韩精品| 午夜免费激情av| 能在线免费观看的黄片| 真人做人爱边吃奶动态| h日本视频在线播放| 国内精品美女久久久久久| 亚洲美女黄片视频| 嫩草影院精品99| 欧美xxxx黑人xx丫x性爽| 少妇人妻精品综合一区二区 | 日本精品一区二区三区蜜桃| 日日撸夜夜添| 在线观看午夜福利视频| 夜夜夜夜夜久久久久| 午夜精品国产一区二区电影 | 欧美+日韩+精品| 婷婷色综合大香蕉| 亚洲av免费在线观看| 成人午夜高清在线视频| 中文字幕人妻熟人妻熟丝袜美| 一个人看视频在线观看www免费| 午夜福利高清视频| 91久久精品国产一区二区三区| 美女内射精品一级片tv| 久久精品国产亚洲网站| 搡老熟女国产l中国老女人| 1000部很黄的大片| 国产精品一及| 久久国内精品自在自线图片| 女人被狂操c到高潮| 中文字幕av在线有码专区| 美女大奶头视频| 亚洲欧美清纯卡通| 99热精品在线国产| 蜜臀久久99精品久久宅男| 97热精品久久久久久| 天堂影院成人在线观看| 中文字幕av成人在线电影| 亚洲七黄色美女视频| 尤物成人国产欧美一区二区三区| 三级男女做爰猛烈吃奶摸视频| 国产欧美日韩精品亚洲av| 别揉我奶头~嗯~啊~动态视频| 成人永久免费在线观看视频| 久久精品夜夜夜夜夜久久蜜豆| 春色校园在线视频观看| 看黄色毛片网站| 日韩制服骚丝袜av| 国产精品一区www在线观看| 中文字幕久久专区| 国产亚洲精品久久久久久毛片| 日本精品一区二区三区蜜桃| 日本黄色片子视频| 中文亚洲av片在线观看爽| 少妇的逼好多水| 精品日产1卡2卡| 亚洲国产高清在线一区二区三| av天堂在线播放| 久久综合国产亚洲精品| 国模一区二区三区四区视频| 99热这里只有是精品50| 一级黄色大片毛片| 久久人人精品亚洲av| 亚洲精品影视一区二区三区av| 99久国产av精品| 日韩大尺度精品在线看网址| h日本视频在线播放| 亚洲18禁久久av| 国内精品宾馆在线| 九九爱精品视频在线观看| av在线亚洲专区| 精品人妻视频免费看| 天天躁日日操中文字幕| 人妻少妇偷人精品九色| 精品少妇黑人巨大在线播放 | 久久精品国产99精品国产亚洲性色| 日本一本二区三区精品| 亚洲三级黄色毛片| 淫妇啪啪啪对白视频| 亚洲精品久久国产高清桃花| 日韩强制内射视频| 日产精品乱码卡一卡2卡三| av中文乱码字幕在线| 又粗又爽又猛毛片免费看| 午夜视频国产福利| 日本三级黄在线观看| 欧美xxxx黑人xx丫x性爽| 99riav亚洲国产免费| 日本撒尿小便嘘嘘汇集6| 免费观看精品视频网站| 亚洲成人av在线免费| 亚洲乱码一区二区免费版| 天堂√8在线中文| 99国产精品一区二区蜜桃av| 久久精品国产99精品国产亚洲性色| 看黄色毛片网站| 1000部很黄的大片| 亚州av有码| 男女视频在线观看网站免费| 欧美中文日本在线观看视频| 亚洲人成网站在线观看播放| 免费在线观看影片大全网站| 亚洲av免费在线观看| 美女内射精品一级片tv| 99久国产av精品| 亚洲性久久影院| 国产一区二区三区av在线 | 女的被弄到高潮叫床怎么办| av在线亚洲专区| 此物有八面人人有两片| 床上黄色一级片| 国产淫片久久久久久久久| 日本三级黄在线观看| 三级男女做爰猛烈吃奶摸视频| 久久这里只有精品中国| 国产精品女同一区二区软件| 国产激情偷乱视频一区二区| 亚洲国产精品sss在线观看| 国产精品爽爽va在线观看网站| 久久久久九九精品影院| 国产成人一区二区在线| 久久精品影院6| 一卡2卡三卡四卡精品乱码亚洲| 一区二区三区高清视频在线| 日韩av在线大香蕉| 在线国产一区二区在线| 国产成人91sexporn| 最新在线观看一区二区三区| 亚洲欧美精品自产自拍| 成人午夜高清在线视频| 99久久精品一区二区三区| 婷婷六月久久综合丁香| 少妇人妻精品综合一区二区 | 中文字幕精品亚洲无线码一区| 老熟妇乱子伦视频在线观看| 久久精品综合一区二区三区| 日本三级黄在线观看| 日产精品乱码卡一卡2卡三| 日韩高清综合在线| 久久久久久久久久成人| 色吧在线观看| 国产v大片淫在线免费观看| 国产精品综合久久久久久久免费| 日本爱情动作片www.在线观看 | 精品午夜福利在线看| 欧美极品一区二区三区四区| 人人妻,人人澡人人爽秒播| 男人和女人高潮做爰伦理| 人人妻人人澡人人爽人人夜夜 | 久久中文看片网| 亚洲乱码一区二区免费版| 性色avwww在线观看| 日韩成人伦理影院| 深夜精品福利| 97超级碰碰碰精品色视频在线观看| 国产黄片美女视频| 国产男人的电影天堂91| 国产日本99.免费观看| 国产熟女欧美一区二区| 亚洲无线在线观看| 中文在线观看免费www的网站| 国产三级中文精品| 在线观看午夜福利视频| 国产精品人妻久久久久久| 久久久久久久久久黄片| aaaaa片日本免费| 日本黄大片高清| 亚洲av熟女| 中文字幕av在线有码专区| 国产精品久久久久久精品电影| 97超级碰碰碰精品色视频在线观看| 国产精品亚洲一级av第二区| 舔av片在线| 国产av一区在线观看免费| 日韩欧美精品v在线| 美女内射精品一级片tv| 久久久久免费精品人妻一区二区| 一区福利在线观看| 在线免费观看不下载黄p国产| 国产极品精品免费视频能看的| 偷拍熟女少妇极品色| 国产在线精品亚洲第一网站| 神马国产精品三级电影在线观看| 春色校园在线视频观看| 内射极品少妇av片p| 我要看日韩黄色一级片| 亚洲不卡免费看| av国产免费在线观看| 99热全是精品| 中文亚洲av片在线观看爽| 嫩草影院入口| 久久鲁丝午夜福利片| 综合色av麻豆| 亚洲精品乱码久久久v下载方式| 国产精品久久久久久亚洲av鲁大| 亚洲aⅴ乱码一区二区在线播放| 亚洲丝袜综合中文字幕| 99久久中文字幕三级久久日本| 免费看av在线观看网站| 男人舔奶头视频| 麻豆成人午夜福利视频| 国产真实乱freesex| 欧美+日韩+精品| 国产黄色小视频在线观看| 尾随美女入室| 波多野结衣高清作品| 国产精品久久久久久久久免| 国产在线精品亚洲第一网站| 男女下面进入的视频免费午夜| 久久久国产成人精品二区| 久久热精品热| 国产黄色视频一区二区在线观看 | 亚洲av免费在线观看| 久99久视频精品免费| 99热6这里只有精品| 国产视频一区二区在线看| 免费看美女性在线毛片视频| 亚洲丝袜综合中文字幕| 一a级毛片在线观看| 成人鲁丝片一二三区免费| 精品一区二区三区视频在线观看免费| 国产成人a∨麻豆精品| 日韩在线高清观看一区二区三区| 高清日韩中文字幕在线| 国产精品伦人一区二区| 精品人妻偷拍中文字幕| 能在线免费观看的黄片| 身体一侧抽搐| 国产精品久久久久久av不卡| 黄色配什么色好看| 狂野欧美激情性xxxx在线观看| 免费一级毛片在线播放高清视频| 深夜精品福利| 99热这里只有精品一区| 国内揄拍国产精品人妻在线| 欧美成人免费av一区二区三区| 国产精品久久久久久亚洲av鲁大| 久久九九热精品免费| 波野结衣二区三区在线| 嫩草影院入口| 欧美潮喷喷水| 我的女老师完整版在线观看| 国产精品久久电影中文字幕| 日日啪夜夜撸| 在线看三级毛片| 22中文网久久字幕| 国产麻豆成人av免费视频| 天天躁夜夜躁狠狠久久av| 国产精华一区二区三区| 国产高清不卡午夜福利| 狠狠狠狠99中文字幕| 伦理电影大哥的女人| 日本免费a在线| 成人一区二区视频在线观看| 国产男人的电影天堂91| 精品久久久久久久久久久久久| 精品一区二区免费观看| 久久人人爽人人片av| 自拍偷自拍亚洲精品老妇| a级毛色黄片| 大型黄色视频在线免费观看| 可以在线观看的亚洲视频| 欧美区成人在线视频| 国产亚洲精品久久久com| 日本黄色视频三级网站网址| av卡一久久| 亚洲经典国产精华液单| 少妇裸体淫交视频免费看高清| 午夜免费男女啪啪视频观看 | 俺也久久电影网| 少妇高潮的动态图| 搡老妇女老女人老熟妇| 日日摸夜夜添夜夜添av毛片| 欧美一区二区国产精品久久精品| 国产在视频线在精品| 3wmmmm亚洲av在线观看| 直男gayav资源| 99热网站在线观看| 国产精品久久久久久久久免| avwww免费| 天堂影院成人在线观看| 一a级毛片在线观看| 久久久久久久久久久丰满| 91麻豆精品激情在线观看国产| 国产精品久久电影中文字幕| 久久久久久伊人网av| 最近最新中文字幕大全电影3| 国产不卡一卡二| 中文字幕免费在线视频6| 欧美又色又爽又黄视频| 久久精品国产99精品国产亚洲性色| 久久国产乱子免费精品| 久久人人爽人人片av| 国产av在哪里看| 美女 人体艺术 gogo| 日本与韩国留学比较| 蜜桃亚洲精品一区二区三区| 久久人人精品亚洲av| 午夜福利成人在线免费观看| 在线观看一区二区三区| 亚洲一区高清亚洲精品| 91av网一区二区| 麻豆av噜噜一区二区三区| 一级a爱片免费观看的视频| 极品教师在线视频| 午夜影院日韩av| 亚洲美女搞黄在线观看 | 99九九线精品视频在线观看视频| 伊人久久精品亚洲午夜| 成年女人永久免费观看视频| 国产精品野战在线观看| 床上黄色一级片| 国产伦一二天堂av在线观看| 亚洲七黄色美女视频| 尾随美女入室|