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

    Extended state observer?based pressure control for pneumatic actuator servo systems

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

    John Mohorcic · Lili Dong

    Abstract

    Keywords Pneumatic actuator · Pressure control · ESO-based SMC · ADRC · Impedance control module · Stability proof ·Lyapunov approach

    1 Introduction

    A pneumatic actuator converts energy (typically in the form of compressed air) into linear or rotary mechanical motions.In this paper, we focus on a linear actuator which mainly consists of a cylinder, a piston (or diaphragm attached to a rod), and valves. The piston separates the cylinder into two chambers. When compressed air at high pressure moves into the back chamber and the front chamber is exhausted, the net force on the piston pushes the assembly forward, and vice versa. As the same force (induced by air pressure) is applied to both chambers, there is no movement [1]. Valves connected to the actuator allows air to flow into or out of each chamber. An external controller can turn on or off the valves to control the air flow and thus the pressure. Being a highly reliable, efficient, and safe source of motion control, a pneumatic actuator has been one of the most popular engineering devices in industry. They are broadly applied to automation and manufacturing process such as pick and place motion,robot manipulation, packing, clamping, broaching, and drilling etc. [1]. These applications typically require fast and accurate positioning. Improving positioning control has thus remained an enduring research issue. As mentioned above, the positioning control of pneumatic actuators can be implemented through regulating and directing air flow in chambers via valves. As air mass flow changes, air pressure varies, which in turn changes the force on the piston [2], and hence alters the positions of piston incrementally. It has been confirmed the effectiveness of pressure control on enhancing the performance of pneumatic position servo systems [3].

    However, the nonlinearity, system uncertainties, and disturbances pose challenges on position servo control of pneumatic actuators. Due to the compressibility of air and delay of valves, the governing dynamic equation of the pneumatic actuator is nonlinear [4]. In addition, intermittent loss of air and pressure (caused by imperfect sealing)and variations of load and friction [5] introduce uncertainty and disturbance to the pneumatic actuator system. To overcome the issues, different pressure control strategies have been applied to pneumatic actuators. Typically, pressure control is combined with position control to form a position servo system as given in [6-9]. Pressure control can also act as an independent control method such as the one in [10-13]. In [6, 7], a PI controller is designed to control the air pressure, while a position controller is constructed based on the controlled pressure. Extra feedforward loops and disturbance observers are utilized in [6] and [7], respectively, to “counterbalance the friction force and parameter changes”. In [8], a PID controller with feedback linearization is used in the pressure control inner loop, and another PID controller is employed in the position control outer loop. In addition to [8], feedback linearization is reported in[9-11] on position tracking control of pneumatic actuators.The linearization method provides “theoretic guidance at the initial stage of controller design and leads to a simpler tracking control strategy” [11]. However, its effectiveness is dependent on accurate mathematical model of the actuator and is not robust enough against modeling uncertainties.As an independent control approach, both open-loop and closed-loop pressure control have been generated. The openloop pressure control, as introduced in [12], is based on the working curve between input voltage and output pressure of a proportional valve. Open-loop pressure control has simple structure and is easy to implement. However, it cannot compensate for disturbances and system parameter variations.Among closed-loop robust control loops which involve pressure sensors, a linear quadratic Gaussian (LQG) self-tuning pressure regulator is utilized to control the pressure based on online estimation of system parameters [13]. Compared with open-loop pressure control, the LQG self-tuning controller is more adaptive to variable parameters such as the volume of the chamber and the setting pressure. In addition,intelligent controllers such as hybrid fuzzy PID [14] and fuzzy force feedback controllers [15] are applied to adjust the pressure of pneumatic systems where fuzzy controller handles structured uncertainty. Although adaptive and intelligent controllers are robust against parameter variations,they do not cope well with external disturbances. To reduce the effect of disturbances, sliding mode controller (SMC)has been employed to regulate the pressure of pneumatic actuators [3, 16, 17]. Sliding mode control is also an effective method for positioning control [18-21]. Constructed on the assumption that the unknown terms in the actuator model are bounded, SMC provides an efficient way to maintain system stability in the presence of modeling imprecision and disturbance [18]. As indicated in [3, 16-21], SMC is a proven method to increase the performance of pneumatic position servo system.

    In [22-24], an extended state observer (ESO) is created for active disturbance rejection controller (ADRC), which integrates ESO with a feedback controller for motion control. The ESO estimates not only system states, but also a generalized disturbance which includes both external disturbance and unknown internal dynamics. Thus, it is a robust estimator for the uncertain systems [24]. Recently, ESO has been embedded in adaptive controller and SMC [23, 24]due to its simple configuration with only one tuning parameter (i.e., observer bandwidth) and good observation performance. Especially in [24], an ESO-based SMC successfully controlled the piston position of a pneumatic actuator.However, the dynamics for air pressure is disregarded in [24]because “it will complicate the control loop” for positioning.“For simplicity” of controller design, only the dynamics of piston position is considered in [24].

    In this paper, an ESO-based SMC is developed based on both piston-load dynamics and chamber pressure model. A single closed-loop position servo system is constructed for pressure and position controls. Compared to double control loops in [6-9], the single control loop has simpler structure.Our control objective is to drive the piston to reach a goal position through regulating the air pressure in chambers.To this end, an impedance control module, which regulates the relationship between force and position of the piston,is originally designed and combined with the ESO-based SMC. While the former is employed to produce a reference pressure level based on desired position, the latter regulates the pressure output of chambers to follow the setpoint. In addition, the extended state of the ESO can approximate the external disturbances and model uncertainties (or generalized disturbance) caused by an air leak, friction, and environmental changes. A SMC then compensates for the estimated generalized disturbance. Therefore, the robustness and stability of the control system are achieved. Unlike the adaptive and intelligent control algorithms in [13-15] where there are at least five tuning parameters, the proposed control system has three tuning parameters (one observer bandwidth for ESO and two controller gains for SMC), making it computationally tractable and easier to implement. In this paper, the ESO-based SMC is compared with classic ADRC and SMC in terms of position tracking and robustness via simulations. The stability of the proposed control system is proved through Lyapunov approach.

    The rest of this paper is organized as follows. The dynamics of load-piston and chamber pressure model are presented in Sect. 2. The ESO-based SMC, ADRC and classic SMC are developed on a pneumatic actuator in Sect. 3. The simulation results are demonstrated in Sect. 4. Concluding remarks and future research are given in Sect. 5.

    2 Modeling of pneumatic actuators

    A linear pneumatic actuator typically “includes a force element (cylinder), a command device (valve), and position and pressure sensors” [25]. It is driven by compressed air through valves and outputs displacement.

    2.1 Piston?load dynamics

    The cylinder consists of a piston, rod, and two chambers.The difference between the air pressure in two chambers determines the position, force, and compliance of the piston.The force diagram of piston and load is shown in Fig. 1.

    Fig. 1 Force diagram of piston and load [25]

    Fig. 2 Pneumatic valves connected to chamber 1

    In Fig. 1,P1(t) andP2(t) are the absolute air pressures for chambers 1 and 2 respectively,FL(t) is payload force,Ff(t) is coulomb friction,βis viscous friction coefficient,MPis mass of the piston,MLis mass of the load, andPa(t)is the atmospheric pressure. The parameterAprefers to the area of the smaller rod portion of the piston, whileA1represents the area of larger disc portion of the piston that divides chambers 1 and 2, andA2is the difference betweenA1andAp. From [25], the piston-load dynamics can be modeled as

    wherex(t) is piston position. The right side of (1) “represents the actuator active force, produced by the different pressures acting on the opposite sides of the piston” [25]. To control the piston position, one must “finely tune the pressure levels in the cylinder chambers using the pneumatic valve” [25].

    2.2 Chamber pressure model

    We suppose the compressed air (or gas) enters chamber 1 and exists chamber 2. According to [25-27], the heat exchange within the respective chambers is represented by

    where˙Q(t) is the heat rate for the air volume,˙mithe mass flow rate for each chamber (where subscriptiindicates chamber number andi=1 or 2 ),˙mlthe mass flow rate of air leakage due to imperfect sealing of cylinder,˙E(t) the rate of change of total energy,˙W(t) the rate of work delivered to

    2.3 Valve model

    As shown in Fig. 1, there are two chambers in one cylinder.“One chamber is connected to the pressure tank through the supply path, and the compressed air will flow inward.The other chamber is connected to the atmosphere through the exhaust path, and the air will flow outward” [25]. Each cylinder chamber is attached to a pair of valves. Figure 2 shows the two valves tied to chamber 1. The pair of valves for each chamber is mutually exclusive so that air cannot be injected and expelled simultaneously. In Fig. 2,Kirepresents the on/off state of each valve andPs(t) is supply pressure. As valve 1 is on (K1= 1), the compressed air will be supplied to chamber 1. As valve 2 is on (K2= 1), the air in chamber 1 will flow outward.

    Therefore, the mass flow rates for chambers 1 and 2 can be represented by (13) and (14) whereqm(defined in [25])is a function ofPs(t), and/orPi(t), and/orPa(t).

    As one turns on/off different valves, the air mass flow rate will change accordingly. We can control the air pressure for each chamber through varying mass flow rate.

    3 Controller design

    The control objective is to regulate the piston to a desired position through adjusting the air pressure in the presence of external disturbance and parameter variations. A block diagram of the pneumatic position servo control system is given in Fig. 3, wherexdis desired position,PDis desired pressure,u(t) is control signal,P(t) is the pressure output of each chamber, and˙x(t) is velocity. The purpose of impedance control module is to producePDaccording toxd. The control system makesP(t) trackPD. In this section, ESO and SMC are introduced first. Then, the ESO-based SMC is derived on the grounds of ESO and SMC.

    3.1 Impedance control module

    Mechanical impedance is defined as force divided by velocity. Thus, if a motion is taken as an input and force as an output, an impedance equation can be established to derive the force from motion. The impedance control module merely represents a mass-spring-damper behavior that regulates force to displacement, velocity, and acceleration. In this section, such a module will be created to yield a desired pressure from a desired positionxd. Equation (1) can be rewritten as

    Fig. 3 Block diagram of pneumatic servo control system

    In (18) and (19),P1D(t) andP2D(t) are the desired pressures for chambers 1 and 2, respectively,Psum(t) is a chosen parameter for the total supplied pressure, andPa(t) is known.Thus, solving (18) and (19) will produceP1D(t) andP2D(t) ,which can be utilized as reference signals for pressure control. For simplicity, the subscripti(i=1 or 2) is disregarded in the rest of Sect. 3 as the same controller can be designed on each chamber.

    3.2 Extended state observer

    3.3 Sliding mode controller

    whereΦ(t) is one of the tuning parameters for SMC.

    3.4 ESO?based SMC

    In Sect. 3.3, the SMC is developed based on the assumption that the feedback signals of position, velocity, and pressure,and the bound of external disturbance are available. In practice, velocity is typically “obtained by differentiation of the position signal with a digital Butterworth low-pass filter”[19]. Thus, the measurement of velocity signal is prone to errors and bandwidth fluctuations [2]. In addition, the bound of disturbance is unknown. Thus, the accurate information of velocity and disturbance are absent in the real world. In this section, to make the SMC more practical, ESO is added to the control system to estimate generalized disturbance which contains both velocity and external disturbance. Based on the outputs of ESO, SMC can effectively reduce the effect of disturbance. In addition, ESO only has one tuning parameterωo. Therefore, it does not add much complexity to the controller design.

    The ESO in (22) is rewritten as

    To reduce chattering effect, the sign function in (53) can be replaced by a saturation function as given by (43), wheresis replaced byz1-PD.

    4 Simulation results

    In this section, ESO-based SMC is compared with classic ADRC and SMC. According to [22-24], the control law for ADRC is designed as

    whereωcis the controller bandwidth and a positive real number. The ADRC is a state feedback controller based on estimated statesz1(t) andz2(t) from ESO. In (54), the difference betweenPDandz1is tracking error. Both (23) and(54) constitute the ADRC control system. The SMC control law is given by (42). As shown in Fig. 3, ESO-based SMC is represented by (44) and (53). Three controllers are simulated on the pneumatic actuator (described by (22)), respectively.Our control objectives are to drive the piston to desired position through adjusting the air pressure to required level. To change air pressure, control commands need to turn on or off valves to direct the air flow into or out of each chamber.Frequent on and off commands will wear out the valves.Thus, we also aim at producing the least excitation signals to valves. The performances of these three controllers are compared in terms of position tracking, pressure regulation,amount of control effort (required to excite valves), as well as robustness against system uncertainties and disturbances.From [25-28], the system parameter values are listed in Table 1, which shows the reference model for impedance control module is critically damped. The controller parameters are selected as in Table 2. For fair comparison, three controllers are tuned such that the magnitudes of their control efforts are comparable. Although ADRC has a simpler structure than ESO-based SMC, the controller parameters of the former have much larger values than that of the latter (asshown in Table 2). This is because ADRC is derived from tracking error while ESO-based SMC is built on the bound of tracking error. The reduced controller gains (for ESObased SMC) are advantageous for practical exercises where small gains are easy to implement.

    Table 1 System parameters

    Table 2 Controller parameters

    4.1 Pressure control and position tracking

    The goal of the pneumatic position servo control system is to command the piston to any position within the cylinder.While position is not directly regulated by the controllers, it is a direct outcome of the controlled pressure in each chamber. Figure 4 shows the position outputs of three controllers and reference trajectory. As seen in Fig. 4, the position changes att=0.2 s and 0.4 s , respectively. Figure 5 shows the steady-state errors (SSE) in the position outputs. The average SSE for ADRC is 4.2 × 10-9m, for ESO-based SMC is 2.5 × 10-9m, and for SMC is - 1.7 × 10-7m. Although all three controllers can drive the piston to reference position,the ESO-based SMC provides the smallest SSE. The SSE of SMC is consistently below zero. Therefore, the tracking performance of SMC is worse than the other two. The chamber pressures under three controllers are shown in Fig. 6.As reference position varies, the chamber pressures change accordingly to push the piston to the goal position. Among three controllers, only SMC controlled pressure has an overshoot att=0.4 s . The performances of ADRC and ESObased SMC are equivalent in pressure regulation.

    Fig. 4 Position outputs

    Figure 7 shows the SSEs for chamber pressures. Specifically, it displays the error between desired and measured pressures. In Fig. 7, “Ch. 1” represents chamber 1, and “Ch.2” represents chamber 2. The average SSEs for ESO-based SMC and ADRC are almost 0.05 kPa and 1 kPa, respectively. However, the SSE for SMC is about 5 kPa. SMC shows the largest SSE. As discussed in Sect. 2, positive mass flow rate turns on valve 1 to provide supply pressure.Negative mass flow rate turns on valve 2 to exhaust air. If the mass flow rate is zero, both valves are off. The excitations of these valves must be scrutinized since they will degrade the mechanical parts gradually. Figure 8 shows the valve states in chamber 1 under the controls of ESO-based SMC, ADRC, and SMC. From simulation, we can see that SMC excites the valves most frequently, and hence has the worst performance among the three controllers. ESO-based SMC performs slightly better than ADRC, particularly at the beginning of the simulation. During the simulation period,ESO-based SMC switches the valves every 1.2 ms on average, ADRC switches the valves every 0.6 ms on average, and SMC switches the valves every 27 μs on average. From [27],the switching valve can change its state on a microsecond timescale to adapt the rapid change of pressure. Thus, both ESO-based SMC and ADRC are feasible. The control efforts of three controllers are shown in Fig. 9. While the average magnitudes of the controllers are nearly the same, SMC exhibits spikes and chattering in the effort. Thus, ESO-based SMC and ADRC outperform SMC.

    Fig. 6 Desired and actual chamber pressure for a ESO-based SMC, b ADRC, and c SMC

    Fig. 7 SSEs for chamber pressures under the control of a ESO-based SMC, b ADRC, and c SMC

    Fig. 8 Valve states in chamber 1 with a ESO-based SMC, b ADRC,and c SMC

    Fig. 9 Control effort of a ESO-based SMC, b ADRC, and c SMC

    Fig. 10 Position outputs with parameter variations

    Figure 4 through Fig. 9 show that ADRC, ESO-based SMC and SMC successfully drive the piston to desired positions through adjusting chamber pressures to required levels using comparable amount of control efforts. Thus, the control goal is achieved by the three controllers. With that being said, ESO-based SMC outperforms the other two controllers in that it produces the least frequently changing signals to excite valves. In addition, it has smaller controller gains than ADRC and SMC. Thus, it is pragmatic and feasible. On the other hand, SMC shows chattering in its control effort, and excites the valves most frequently. Therefore, it is not an ideal choice compared to the other two. In the next section,we aim to test the robustness of the three controllers.

    4.2 Robust analysis

    Fig. 11 Pressure outputs with parameter variations: a ESO-based SMC, b ADRC, and c SMC

    First, we vary the system parameter values to introduce structured uncertainty to the pneumatic actuator system.The impedance parameters are varied withincreased by 1%,reduced by 80%, andreduced by 95%. Therefore, the impedance system remains critically damped. The position outputs of three control systems are illustrated in Fig. 10. Despite the presence of parameter variations, all three controllers drive the piston to desired positions but with a bit longer settling time than the one in Fig. 5. This is expected because friction force is increased with the change of impedance parameters. The resulting chamber pressure and control efforts are shown in Figs. 11 and 12, respectively, where ADRC and ESO-based SMC show comparable performances, while SMC displays overshoots in pressure outputs and chattering effects in its control effort. The average magnitudes of three control signals are almost same.

    Subsequently, we need to test the robustness of the three controllers against external disturbance. An air leak is a common problem in a pneumatic actuator system. It could be caused by a loose connection, a hole in the tubing, or wear to a seal and so on. Here, a step disturbance representing a sudden air leak is added to the chamber pressure model att=0.5 s . Figures 13 through 16 show how three control systems respond to the disturbance. The position outputs,the chamber pressure, and control efforts are depicted in Figs. 13, 14, 15, and 16, respectively. Despite the presence of the disturbance, the three controllers command the piston to the desired positions. Both ESO-based SMC and ADRC effectively compensate for the disturbance. However, the SMC yields SSE in position output aftert=0.5 s (as shown in Fig. 14).

    In addition, ADRC and ESO-based SMC successfully control the chamber pressures to desired levels with equivalent amount of control efforts (see Figs. 15 and 16). There are no noticeable fluctuations in the position and pressure responses of ADRC and ESO-based SMC, even with the existing disturbance. However, SMC shows evident chattering in the control input.

    4.3 Conclusions

    Fig. 12 Control efforts of a ESO-based SMC, b ADRC, and c SMC in the presence of parameter variations

    In this paper, an ESO-based SMC is developed to improve the performance of a pneumatic position servo system through regulating the air pressure of each chamber in the pneumatic cylinder. ESO, an accurate observer that was originally used in ADRC, is applied to estimate the velocity and generalized disturbance. SMC, a proven method for positioning pneumatic actuators, is employed to control the air pressure based on the estimated states of ESO. ESO and SMC form a robust controller that can successfully drive the piston to desired position through regulating the air pressure,despite the presence of system uncertainty and external disturbance. The Lyapunov approach proves the stability of the control system. In addition, for comparison purposes, both ADRC and SMC are designed and applied to control the air pressure, respectively. The performances of three controllers(i.e., ADRC, ESO-based SMC and SMC) are evaluated in simulation. For a fair assessment, the three controllers are tuned in a way such that the average magnitudes of control efforts are nearly the same. It is shown that both ADRC and ESO-based SMC outperform SMC in terms of robustness and control effort. When it comes to pressure regulation and positioning, ADRC and ESO-based SMC demonstrate equivalent results. Nevertheless, ESO-based SMC produces less valve excitation than ADRC to accomplish the control goal, particularly at the beginning of the simulation. In addition, ESO-based SMC requires smaller controller gains than both ADRC and SMC.

    Fig. 13 Position outputs with input disturbance

    Fig. 14 Zoom-in view of position outputs as t = 0.5-0.55 s

    One contribution of the paper is the original and successful development of ESO-based SMC for the pressure control of a pneumatic position servo system. Another contribution is the establishment of impedance control module, on which the ESO-based SMC is constructed. The impedance control module translates the desired position to required air pressure level for each chamber. With the module, piston-load dynamics and chamber pressure model are connected in such a way that it enables the controller to position the actuator through regulating air pressure.

    In the future, we plan to implement the ADRC and ESObased SMC on a real pneumatic actuator. The experimental results are expected to match the simulation results presented in the paper.

    Fig. 15 Chamber pressure in the presence of disturbance under the control of a ESO-based SMC, b ADRC, and c SMC

    Fig. 16 Control efforts of a ESO-based SMC, b ADRC, and c SMC in the presence of disturbance

    亚洲精品乱码久久久v下载方式| 99久久99久久久精品蜜桃| 国产精品1区2区在线观看.| 精品人妻1区二区| 亚洲18禁久久av| 久久6这里有精品| 桃色一区二区三区在线观看| 看片在线看免费视频| 成年女人看的毛片在线观看| 国产精品av视频在线免费观看| 偷拍熟女少妇极品色| 可以在线观看的亚洲视频| 中文在线观看免费www的网站| 国产精品国产高清国产av| 精品久久久久久,| 麻豆国产97在线/欧美| 午夜精品在线福利| 99久久成人亚洲精品观看| 亚洲av电影不卡..在线观看| 国产成人a区在线观看| 久久久久国产精品人妻aⅴ院| 国产三级在线视频| 精品午夜福利在线看| 欧美成人一区二区免费高清观看| 色综合欧美亚洲国产小说| 91狼人影院| 男人舔奶头视频| 国产精品三级大全| 国产淫片久久久久久久久 | 欧美最新免费一区二区三区 | 午夜激情福利司机影院| av专区在线播放| 亚洲色图av天堂| 免费观看精品视频网站| 脱女人内裤的视频| 高清在线国产一区| 国产精品久久电影中文字幕| 国产在视频线在精品| 99热这里只有是精品50| 天堂av国产一区二区熟女人妻| 亚洲国产日韩欧美精品在线观看| 又黄又爽又免费观看的视频| 亚洲精品一卡2卡三卡4卡5卡| 亚洲中文字幕日韩| 亚洲专区中文字幕在线| 无遮挡黄片免费观看| 久久精品国产清高在天天线| 国产精品三级大全| 日韩欧美免费精品| 麻豆成人午夜福利视频| 免费观看精品视频网站| 1000部很黄的大片| 欧美精品国产亚洲| 在现免费观看毛片| 国产成人啪精品午夜网站| 欧美丝袜亚洲另类 | 淫妇啪啪啪对白视频| 日本熟妇午夜| 国产老妇女一区| 国产乱人视频| 国产高潮美女av| 亚洲国产高清在线一区二区三| 日本免费一区二区三区高清不卡| 深夜a级毛片| 高清日韩中文字幕在线| 国产黄a三级三级三级人| 久久久久久大精品| 搡老妇女老女人老熟妇| 最后的刺客免费高清国语| 精品一区二区三区视频在线| 99热这里只有是精品在线观看 | 精品人妻视频免费看| 国产私拍福利视频在线观看| 国产黄色小视频在线观看| 国产成人欧美在线观看| 久久久久久久亚洲中文字幕 | 亚洲欧美清纯卡通| 色综合站精品国产| 亚洲精品影视一区二区三区av| 99热6这里只有精品| 中文字幕久久专区| 麻豆成人av在线观看| 成年女人看的毛片在线观看| 亚洲av成人不卡在线观看播放网| 在线观看一区二区三区| 国产黄a三级三级三级人| 91在线精品国自产拍蜜月| 男人和女人高潮做爰伦理| 非洲黑人性xxxx精品又粗又长| 亚洲欧美日韩高清专用| aaaaa片日本免费| 国产精品永久免费网站| 国产激情偷乱视频一区二区| 欧美激情久久久久久爽电影| 欧美日韩黄片免| 久久人人精品亚洲av| 色吧在线观看| 九色成人免费人妻av| 搡老妇女老女人老熟妇| 成年人黄色毛片网站| 欧美3d第一页| 级片在线观看| 国产伦一二天堂av在线观看| 国产精品一区二区性色av| 中文资源天堂在线| 91久久精品电影网| 久久久精品欧美日韩精品| 国产单亲对白刺激| 少妇高潮的动态图| 九九久久精品国产亚洲av麻豆| 国模一区二区三区四区视频| 观看免费一级毛片| 身体一侧抽搐| 变态另类丝袜制服| 午夜激情欧美在线| 国产亚洲精品av在线| 国产午夜精品久久久久久一区二区三区 | 老司机午夜十八禁免费视频| 国产伦在线观看视频一区| 美女大奶头视频| 欧美在线一区亚洲| 国产精品永久免费网站| 999久久久精品免费观看国产| 国产一区二区在线av高清观看| 很黄的视频免费| 国产成+人综合+亚洲专区| 日韩中字成人| 99riav亚洲国产免费| 精品久久久久久久久久免费视频| 日韩欧美 国产精品| 嫩草影院精品99| 老司机午夜福利在线观看视频| 婷婷精品国产亚洲av| 我的女老师完整版在线观看| 1024手机看黄色片| 可以在线观看毛片的网站| 熟女电影av网| 欧美色欧美亚洲另类二区| 国内精品久久久久久久电影| 老司机午夜福利在线观看视频| 国产激情偷乱视频一区二区| 亚洲色图av天堂| 精品一区二区三区人妻视频| 亚洲中文字幕日韩| 亚洲 国产 在线| 久久精品久久久久久噜噜老黄 | 亚洲人成网站在线播放欧美日韩| 亚洲乱码一区二区免费版| 久久人人爽人人爽人人片va | 直男gayav资源| 欧美激情国产日韩精品一区| 91久久精品电影网| 我的女老师完整版在线观看| 一区二区三区四区激情视频 | 少妇的逼水好多| 国产免费一级a男人的天堂| 国产精品精品国产色婷婷| 精品日产1卡2卡| 夜夜看夜夜爽夜夜摸| 国产三级黄色录像| 亚洲精品456在线播放app | 国产成年人精品一区二区| 精品熟女少妇八av免费久了| 国产探花极品一区二区| 我的女老师完整版在线观看| 美女 人体艺术 gogo| 久久久久国产精品人妻aⅴ院| 97超视频在线观看视频| 免费在线观看亚洲国产| 欧美激情久久久久久爽电影| 日本精品一区二区三区蜜桃| 精品无人区乱码1区二区| 国内精品一区二区在线观看| 深夜精品福利| 99久久精品国产亚洲精品| 伦理电影大哥的女人| 美女高潮喷水抽搐中文字幕| 国产中年淑女户外野战色| 黄色一级大片看看| 久久精品影院6| 成年免费大片在线观看| 亚洲专区国产一区二区| 欧美一区二区国产精品久久精品| 老司机福利观看| 国产极品精品免费视频能看的| 国产一区二区三区在线臀色熟女| 免费黄网站久久成人精品 | 中出人妻视频一区二区| 国产午夜福利久久久久久| 国产精品自产拍在线观看55亚洲| 色尼玛亚洲综合影院| 久久久久久九九精品二区国产| 性色av乱码一区二区三区2| 级片在线观看| 亚洲精品在线观看二区| 久久天躁狠狠躁夜夜2o2o| 亚洲av中文字字幕乱码综合| 国产主播在线观看一区二区| 国产淫片久久久久久久久 | 乱人视频在线观看| 又爽又黄a免费视频| 大型黄色视频在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲av成人精品一区久久| 久久国产乱子伦精品免费另类| 午夜老司机福利剧场| netflix在线观看网站| 久久精品国产亚洲av香蕉五月| 国产成人aa在线观看| 日韩精品中文字幕看吧| 有码 亚洲区| 蜜桃久久精品国产亚洲av| 久久国产精品影院| av欧美777| 天美传媒精品一区二区| 最近在线观看免费完整版| 非洲黑人性xxxx精品又粗又长| 搞女人的毛片| 他把我摸到了高潮在线观看| 亚洲精品久久国产高清桃花| 欧洲精品卡2卡3卡4卡5卡区| 国产伦一二天堂av在线观看| 欧美又色又爽又黄视频| a级毛片免费高清观看在线播放| 免费观看的影片在线观看| 久久伊人香网站| 免费无遮挡裸体视频| 婷婷精品国产亚洲av在线| 亚洲av电影在线进入| 国产中年淑女户外野战色| 国产一区二区三区在线臀色熟女| 日韩 亚洲 欧美在线| 丁香六月欧美| 国产高清视频在线播放一区| 九色成人免费人妻av| 日韩亚洲欧美综合| 日日摸夜夜添夜夜添av毛片 | 亚洲五月婷婷丁香| 国产精品亚洲美女久久久| 亚洲18禁久久av| 成人亚洲精品av一区二区| 91午夜精品亚洲一区二区三区 | 国产在线精品亚洲第一网站| 搡老岳熟女国产| 别揉我奶头~嗯~啊~动态视频| 黄色一级大片看看| 亚洲无线在线观看| 中文在线观看免费www的网站| 精品久久久久久,| 一级作爱视频免费观看| h日本视频在线播放| 亚洲性夜色夜夜综合| 国内久久婷婷六月综合欲色啪| 久久香蕉精品热| 国内揄拍国产精品人妻在线| 两个人的视频大全免费| 国产白丝娇喘喷水9色精品| www.色视频.com| 国产中年淑女户外野战色| 伊人久久精品亚洲午夜| 亚洲自拍偷在线| 男人的好看免费观看在线视频| 国产精品电影一区二区三区| 两个人的视频大全免费| 亚洲无线在线观看| 中出人妻视频一区二区| 欧美三级亚洲精品| 成年女人永久免费观看视频| 精品久久久久久成人av| 亚洲中文字幕一区二区三区有码在线看| 淫秽高清视频在线观看| 尤物成人国产欧美一区二区三区| 亚洲精品在线美女| 久久久精品欧美日韩精品| 成人亚洲精品av一区二区| 国产精品久久久久久久电影| 欧美黄色片欧美黄色片| 午夜精品在线福利| 成人亚洲精品av一区二区| 国语自产精品视频在线第100页| 99久久精品热视频| 18美女黄网站色大片免费观看| 欧美日韩亚洲国产一区二区在线观看| 精品不卡国产一区二区三区| 国产精品一区二区三区四区久久| 国产精品久久视频播放| 国产精品爽爽va在线观看网站| 一本一本综合久久| 可以在线观看毛片的网站| 美女免费视频网站| 老鸭窝网址在线观看| 亚洲av成人精品一区久久| 搡老岳熟女国产| 国产精品嫩草影院av在线观看 | 国产成人啪精品午夜网站| 特大巨黑吊av在线直播| 女同久久另类99精品国产91| a级毛片免费高清观看在线播放| 丰满人妻熟妇乱又伦精品不卡| 一进一出抽搐gif免费好疼| 88av欧美| 色视频www国产| 看十八女毛片水多多多| 国产成人福利小说| 亚洲一区高清亚洲精品| 在线播放无遮挡| 少妇丰满av| 男女视频在线观看网站免费| 精品人妻一区二区三区麻豆 | 亚洲欧美精品综合久久99| 一级黄片播放器| 九九热线精品视视频播放| 两个人视频免费观看高清| 最新在线观看一区二区三区| 韩国av一区二区三区四区| 亚洲精品日韩av片在线观看| 欧美激情在线99| 国产精品久久电影中文字幕| 在线天堂最新版资源| 波多野结衣高清作品| 99热6这里只有精品| 欧美日韩国产亚洲二区| 精品国产三级普通话版| 禁无遮挡网站| netflix在线观看网站| 热99在线观看视频| 麻豆成人午夜福利视频| 欧美乱色亚洲激情| 香蕉av资源在线| 麻豆国产av国片精品| 丁香六月欧美| 日韩成人在线观看一区二区三区| 夜夜爽天天搞| 窝窝影院91人妻| 五月伊人婷婷丁香| 国模一区二区三区四区视频| 看片在线看免费视频| 国产免费一级a男人的天堂| 99国产综合亚洲精品| 国产av一区在线观看免费| 此物有八面人人有两片| 久久6这里有精品| 国内精品久久久久精免费| 国产三级在线视频| 日韩高清综合在线| 国产在线男女| 久久精品久久久久久噜噜老黄 | 一个人免费在线观看电影| 男人和女人高潮做爰伦理| 国产午夜精品论理片| 久久中文看片网| 欧美黄色淫秽网站| 我的老师免费观看完整版| 搡老妇女老女人老熟妇| 人妻丰满熟妇av一区二区三区| 亚洲精品成人久久久久久| 99精品久久久久人妻精品| 网址你懂的国产日韩在线| 久久精品综合一区二区三区| 校园春色视频在线观看| 悠悠久久av| 久久6这里有精品| 小蜜桃在线观看免费完整版高清| 国产蜜桃级精品一区二区三区| 非洲黑人性xxxx精品又粗又长| 小说图片视频综合网站| 亚洲性夜色夜夜综合| 国产av一区在线观看免费| 欧美日韩黄片免| 日韩欧美免费精品| 午夜精品一区二区三区免费看| 国产男靠女视频免费网站| 老鸭窝网址在线观看| 天堂av国产一区二区熟女人妻| 亚洲精品日韩av片在线观看| 好男人在线观看高清免费视频| 天天一区二区日本电影三级| 99在线人妻在线中文字幕| 欧美国产日韩亚洲一区| 日本五十路高清| 国产亚洲av嫩草精品影院| 国产精品乱码一区二三区的特点| 久久人人爽人人爽人人片va | 在线十欧美十亚洲十日本专区| 久久人人爽人人爽人人片va | 国产一区二区三区视频了| 最近中文字幕高清免费大全6 | 精品日产1卡2卡| ponron亚洲| 亚洲五月天丁香| 搡老熟女国产l中国老女人| 美女高潮喷水抽搐中文字幕| 蜜桃亚洲精品一区二区三区| 国产高清视频在线播放一区| 久久午夜福利片| 久久久久久久久久成人| 哪里可以看免费的av片| 日韩欧美在线二视频| 免费在线观看影片大全网站| 亚洲av熟女| 午夜免费男女啪啪视频观看 | 国产乱人伦免费视频| 亚洲最大成人av| 国产色爽女视频免费观看| 五月伊人婷婷丁香| 国产日本99.免费观看| 色哟哟哟哟哟哟| 亚洲国产精品成人综合色| 成人欧美大片| 日韩欧美一区二区三区在线观看| 女人被狂操c到高潮| 国产高潮美女av| 十八禁人妻一区二区| 欧美乱色亚洲激情| 九九热线精品视视频播放| 欧美一级a爱片免费观看看| 久久亚洲精品不卡| 欧美黄色淫秽网站| 亚洲精品亚洲一区二区| 国产精华一区二区三区| 国产精品国产高清国产av| av在线天堂中文字幕| 脱女人内裤的视频| 欧美午夜高清在线| 老司机福利观看| 自拍偷自拍亚洲精品老妇| 99在线视频只有这里精品首页| 在线a可以看的网站| 成人午夜高清在线视频| netflix在线观看网站| 国产精品三级大全| 亚洲在线自拍视频| 久久精品久久久久久噜噜老黄 | 午夜老司机福利剧场| 黄色配什么色好看| 中文字幕熟女人妻在线| 国产探花在线观看一区二区| 国产欧美日韩一区二区三| 欧美成人免费av一区二区三区| 国产三级在线视频| 99热6这里只有精品| 91久久精品电影网| 久久九九热精品免费| 欧美最新免费一区二区三区 | 亚洲av成人精品一区久久| 他把我摸到了高潮在线观看| 1024手机看黄色片| 在线十欧美十亚洲十日本专区| 国产精品女同一区二区软件 | 亚洲经典国产精华液单 | 欧美国产日韩亚洲一区| 亚洲乱码一区二区免费版| 无人区码免费观看不卡| 国产高清视频在线播放一区| 一本一本综合久久| 日韩中文字幕欧美一区二区| 夜夜爽天天搞| 婷婷亚洲欧美| av视频在线观看入口| 一本综合久久免费| 听说在线观看完整版免费高清| a级一级毛片免费在线观看| 亚洲av成人av| 欧美国产日韩亚洲一区| 欧美在线黄色| 99久国产av精品| 亚洲国产精品sss在线观看| 日本撒尿小便嘘嘘汇集6| 亚洲专区中文字幕在线| 亚洲人成电影免费在线| 国产极品精品免费视频能看的| 性色avwww在线观看| 亚洲av二区三区四区| 99热这里只有是精品在线观看 | 亚洲欧美激情综合另类| 日韩精品青青久久久久久| 成人性生交大片免费视频hd| 三级国产精品欧美在线观看| 美女cb高潮喷水在线观看| 黄色日韩在线| 成人午夜高清在线视频| 在线十欧美十亚洲十日本专区| 亚洲av第一区精品v没综合| 亚洲一区二区三区不卡视频| or卡值多少钱| 赤兔流量卡办理| 亚洲国产欧美人成| 美女被艹到高潮喷水动态| 国产一区二区在线av高清观看| 免费在线观看日本一区| 国产高清视频在线观看网站| 国产色爽女视频免费观看| 久久久久久久久久黄片| 午夜久久久久精精品| 他把我摸到了高潮在线观看| 中文字幕人妻熟人妻熟丝袜美| 精品一区二区免费观看| 丰满人妻熟妇乱又伦精品不卡| 丁香欧美五月| 好男人在线观看高清免费视频| av在线蜜桃| 欧美色视频一区免费| 国产一区二区三区视频了| 最近在线观看免费完整版| 12—13女人毛片做爰片一| 男人和女人高潮做爰伦理| 两个人的视频大全免费| 麻豆国产av国片精品| 午夜福利高清视频| 亚洲人与动物交配视频| 久久热精品热| 99热6这里只有精品| 亚洲精品一区av在线观看| 99久久精品国产亚洲精品| 人人妻,人人澡人人爽秒播| 国内精品久久久久久久电影| 国产精品久久久久久亚洲av鲁大| 三级国产精品欧美在线观看| 亚洲人成电影免费在线| 嫩草影视91久久| 老司机午夜十八禁免费视频| 人人妻人人澡欧美一区二区| 日本黄色片子视频| 五月伊人婷婷丁香| 亚洲国产精品久久男人天堂| 高清日韩中文字幕在线| 又黄又爽又刺激的免费视频.| 欧美黑人欧美精品刺激| 欧美一级a爱片免费观看看| 久久久久久国产a免费观看| 亚洲av中文字字幕乱码综合| 淫秽高清视频在线观看| 91在线观看av| 日本免费a在线| 亚洲国产精品合色在线| 日韩av在线大香蕉| 欧美日韩综合久久久久久 | 久久久久久久精品吃奶| 国产精品98久久久久久宅男小说| 18禁裸乳无遮挡免费网站照片| 夜夜看夜夜爽夜夜摸| 午夜精品久久久久久毛片777| 1000部很黄的大片| 久久久色成人| 国内少妇人妻偷人精品xxx网站| 无人区码免费观看不卡| 午夜福利18| 校园春色视频在线观看| 国产精品一及| 日韩大尺度精品在线看网址| 亚洲内射少妇av| 久久久久久久久中文| 日韩欧美免费精品| 国产黄片美女视频| 十八禁人妻一区二区| 两个人的视频大全免费| 免费在线观看成人毛片| 国产真实伦视频高清在线观看 | a级一级毛片免费在线观看| 18禁在线播放成人免费| 又黄又爽又免费观看的视频| 成人无遮挡网站| 精品一区二区三区视频在线观看免费| 精品人妻熟女av久视频| 国产野战对白在线观看| 99久国产av精品| 91麻豆精品激情在线观看国产| 午夜精品一区二区三区免费看| 露出奶头的视频| 内射极品少妇av片p| av视频在线观看入口| 小说图片视频综合网站| 久久午夜亚洲精品久久| 此物有八面人人有两片| 精品午夜福利视频在线观看一区| 久久久成人免费电影| 99热精品在线国产| 亚洲狠狠婷婷综合久久图片| 天天一区二区日本电影三级| www.www免费av| 免费在线观看影片大全网站| 国产日本99.免费观看| 一区福利在线观看| av在线老鸭窝| 国产一区二区激情短视频| 好看av亚洲va欧美ⅴa在| 女人十人毛片免费观看3o分钟| 成人特级黄色片久久久久久久| 亚洲精华国产精华精| 白带黄色成豆腐渣| 欧美成人一区二区免费高清观看| 12—13女人毛片做爰片一| 日韩亚洲欧美综合| 亚洲电影在线观看av| 99精品久久久久人妻精品| 日韩成人在线观看一区二区三区| 高清在线国产一区| 国产国拍精品亚洲av在线观看| 亚洲无线观看免费| 亚洲aⅴ乱码一区二区在线播放| 久久精品91蜜桃| 在线看三级毛片| 久久亚洲精品不卡| 免费在线观看日本一区| 舔av片在线| 国内精品一区二区在线观看| 久久天躁狠狠躁夜夜2o2o| av在线蜜桃| 亚洲精品456在线播放app | or卡值多少钱| 网址你懂的国产日韩在线| 国产真实乱freesex| 亚洲成人精品中文字幕电影| 国产久久久一区二区三区| 91午夜精品亚洲一区二区三区 |