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

    A Novel Disturbance Observer Based Fixed-Time Sliding Mode Control for Robotic Manipulators With Global Fast Convergence

    2024-03-04 07:43:54DanZhangJiabinHuJunChengZhengGuangWuandHuaichengYan
    IEEE/CAA Journal of Automatica Sinica 2024年3期

    Dan Zhang ,,, Jiabin Hu , Jun Cheng , Zheng-Guang Wu ,,, and Huaicheng Yan ,,

    Abstract—This paper proposes a new global fixed-time sliding mode control strategy for the trajectory tracking control of uncertain robotic manipulators.First, a fixed-time disturbance observer (FTDO) is designed to deal with the adverse effects of model uncertainties and external disturbances in the manipulator systems.Then an adaptive scheme is used and the adaptive FTDO (AFTDO) is developed, so that the priori knowledge of the lumped disturbance is not required.Further, a new non-singular fast terminal sliding mode (NFTSM) surface is designed by using an arctan function, which helps to overcome the singularity problem and enhance the robustness of the system.Based on the estimation of the lumped disturbance by the AFTDO, a fixed-time non-singular fast terminal sliding mode controller (FTNFTSMC)is developed to guarantee the trajectory tracking errors converge to zero within a fixed time.The settling time is independent of the initial state of the system.In addition, the stability of the AFTDO and FTNFTSMC is strictly proved by using Lyapunov method.Finally, the fixed-time NFESM (FTNFTSM) algorithm is validated on a 2-link manipulator and comparisons with other existing sliding mode controllers (SMCs) are performed.The comparative results confirm that the FTNFTSMC has superior control performance.

    I.INTRODUCTION

    WITH the further development of modern industry,robotic manipulators can not only assist humans to carry out mechanical automatic production, but also replace labours to complete dangerous industrial tasks, thus they can improve production efficiency and reduce safety problems.In recent years, due to the unique operational flexibility of robotic manipulators, they have received extensive attention in various fields, such as aerospace, medical treatment and logistics transportation [1]-[5].However, complex application scenarios require the robotic manipulators to have superior control performance.As a typical nonlinear system, the manipulator has a complex mechanical structure and it is easily affected by uncertain factors such as parameter perturbations, unmodeled dynamics and external disturbances [6].Therefore, it is necessary to design effective schemes to accurately control manipulators to satisfy practical application demands.

    In the past few decades, various advanced control algorithms were studied to overcome nonlinear and uncertainty problems in robot systems, see, e.g., adaptive control [7], [8],fuzzy control [9], [10], backstepping control [11], [12] and sliding mode control [13], [14].It is known to all that sliding mode control can provide more excellent robustness than other control schemes, which makes it receive extensive attention.However, the traditional linear sliding mode (LSM)design can only ensure the asymptotic stability of the system,and the undesirable chattering is easy to deteriorate the control performance.Those defects have limited the application of sliding mode controllers (SMCs) in practical industry.In order to improve the control performance of dynamic systems and guarantee the finite-time convergence of system states,the terminal sliding mode (TSM) control scheme was proposed in [15].Although TSM control algorithms have gained popularity, it also has some shortcomings, that is, it is prone to singularity problem.By designing an appropriate fractional power in the sliding surface, a new finite-time control method with the non-singular TSM (NTSM) design was proposed in[16] to overcome the undesired singularity.Inspired by the work in [16], the non-singular terminal SMC (NTSMC) was designed in [17] to solve the trajectory tracking problem of manipulators with unknown uncertainties, and the satisfactory tracking performance was achieved.In addition, compared with NTSM control, the non-singular fast TSM (NFTSM)control introduced in [18]-[20] showed faster convergence speed in trajectory tracking.Noteworthily, finite-time sliding mode control for robotic manipulators has been extensively studied in the past few years, but a defect exists in the current finite-time control results, that is, the settling time of the tracking errors is constrained by the initial system states.For example, the finite-time sliding mode control strategy for robotic manipulators presented in [17] and [20] can guarantee a bounded settling time, but the settling time of the error system will increase when the initial state of the system is far from the equilibrium point.This shortcoming affects the global tracking performance of the robotic manipulator and limits the potential of finite-time control strategies in practical engineering.

    Recently, some new fixed-time control methods were proposed and they immediately receive a lot of research attention,which ensures that the convergence of all system states is uniformly bounded and the upper bound of the convergence time is independent of the initial system states, see, e.g., [21]-[28].In [21], a fixed-time feedback control strategy for nonlinear dynamical systems was investigated, which guaranteed the global fixed-time stability.To solve the trajectory tracking problem ofn-link robotic manipulators, an adaptive fuzzy fixed-time control scheme was proposed in [23], in which an error conversion mechanism based on the system performance and a barrier Lyapunov function were designed.This fixed-time control strategy not only helps the manipulator system achieve rapid transient response, but also guarantees the boundedness of its closed-loop signals, thus improving the global convergence performance and confining the trajectory tracking error within the artificially set boundary.In [24], an adaptive fixed-time cooperative control scheme was proposed for multiple manipulator systems with position and velocity constraints, which guaranteed that the consensus tracking error of each intelligent body converged to a small neighborhood near the origin in a fixed time.With help of the above fixed-time control algorithm, the convergence performance of the manipulator is not constrained by the system state, thus eliminating the need for frequent changes in the control structure.Extending the fixed-time scheme to sliding mode control,a robust fixed-time SMC with variable exponential coefficients was recently proposed in [27], which can effectively deal with mismatched disturbances in nonlinear systems.By using nonlinear piecewise functions to overcome the singularity problem, a fixed-time SMC for robotic manipulators was proposed in [28] to ensure that tracking errors converged to a small field near zero within a fixed time.Note that this scheme does not achieve global fixed-time convergence of tracking errors and singularity avoidance process is very complex.In [29], an improved high-order SMC was used to achieve global fixed-time stability and fast state convergence.However, the above fixed-time SMC is not optimal, and unknown parameter perturbations and external disturbances have not been well dealt with.

    As is well known, the system uncertainties and external disturbances are unavoidable in real systems, which will deteriorate the tracking performance of robotic manipulators.At present, many results related to robots mainly focus on how to eliminate the influence of those uncertain factors.The control scheme based on the disturbance observer (DO) is popular among existing research results, which can help strengthen the robustness of the system and improve the tracking accuracy,see, e.g., [30]-[33].More specifically, in [31], an asymptotic disturbance observer for robot systems was designed to estimate the unknown disturbance, which ensured the rapid convergence of tracking errors.In [32], an adaptive fuzzy integral sliding mode control strategy based on the disturbance observer was investigated for robot systems with parameter perturbations, external disturbances and actuator faults.However, the above design schemes cannot guarantee estimation errors converge within a fixed time, which is easy to affect the settling time of system states.Recently, research effort has been devoted onto improving those control performance in literature.In [34], a nonlinear disturbance observer was developed to estimate the lumped disturbance, which ensured that the estimation error converged within a fixed time.In addition, in order to speed up the convergence speed of robot systems, a fixed-time disturbance observer (FTDO) was proposed in [35]-[37] to solve the system uncertainty problem.However, the above design results are also not optimal in terms of estimation accuracy and convergence rate, and there are very few results reported on the FTDO, which motivates the current study.

    Based on the foregoing discussion, a new fixed-time NFTSM control scheme based on an adaptive FTDO(AFTDO) for manipulators is proposed in this paper.Simultaneously, the major contributions of this paper are highlighted as follows:

    1) Compared with the finite-time sliding mode control methods presented in [15]-[20], a new fixed-time NFTSM control strategy for robotic manipulators is proposed in this paper, which ensures that the trajectory tracking error can converge within a fixed time and the upper bound of the convergence time is independent of the initial state.Therefore, a better convergence performance is achieved, which can be found in the simulation study.

    2) Compared with the existing NFTSM control methods introduced in [20], [26] and [28], a new NFTSM surface based on the arctan function is designed to overcome the singularity problem more simply and effectively, and there is no need to worry about the non-differentiability problem that may exist in the piecewise function.In addition, this new sliding surface ensures that the manipulator system can achieve faster trajectory tracking speed.

    3) Compared with the traditional disturbance observers introduced in [31]-[33], which can only achieve asymptotic or finite-time stability, a new AFTDO is designed to ensure that the estimation error converges in a fixed time.Simultaneously, it is required in the existing FTDO designed in[34]-[37] that the prior knowledge of lumped disturbance must be known and it was difficult to provide satisfactory estimation accuracy.Compared with these results, the proposed AFTDO can achieve faster estimation speed and higher estimation accuracy, and does not need to know the perturbation bound in advance with the help of the adaptive law.

    The remaining part of this paper is proceeded as follows.In Section II, the description of the problem and some preparatory knowledge are introduced.The design process of the fixed-time non-singular fast terminal SMC (FTNFTSMC) and the corresponding stability analysis are presented in Section III.In Section IV, the simulation results show that the proposed control strategy is feasible and provides superior control performance.This paper is summarized in Section V.

    II.PROBLEM FORMULATION

    A. Notations

    Inthispaper, thefollowingsymbols are used.Foranyndimensionalvectorz=[z1,z2,...,zn]T∈Rn, ||z||denotesthe Euclidean norm of the vectorz.sig(z) is defined as sig(z)=[sign(z1),sign(z2),...,sign(zn)]T, where sign(*) is symbolic function.λmin{P} is the minimum eigenvalue of matrixP.

    B. Lemmas and Definitions

    Consider the following nonlinear system:

    wherexdenotes the state variable, andf(x(t)) is a continuous nonlinear function.For any initial state, we assume that there is a unique solution for the system (1).

    Definition 1[27]: If the system (1) satisfies: i) It is globally finite-time stable, ii) The settling timeT(z0) is bounded, and there exists a constantTmaxsuch thatT(x0)≤Tmax, then the system (1) is fixed-time stable.

    Lemma 1[21]: For system (1), if there exists a Lyapunov functionV(x) satisfying

    wherepγ <1,kγ >1, α , β,p,kand γ are positive numbers.Then the dynamic system is fixed-time stable.The settling time satisfies

    Lemma 2[38]: Consider the nonlinear system in (1).If there exists a Lyapunov functionV(x) and positive scalars α,β,p,kand γ, such that

    where ? >0.Then the dynamic system is practical fixed-time stable and the residual set of the solution is expressed as

    where μ is a scalar, 0 <μ≤1.The settling time satisfies

    Lemma 3[39]: For the vectorz=[z1,z2,...,zn]T∈Rn, it satisfies

    wherei=1,2,...,n, 0 <ζ1<1 and ζ2>1.

    C. Dynamics of the Robotic Manipulators

    Consider the following dynamic equation for an-link rigid manipulator:

    Property 1[40]:N(q),B(q,q˙) andG(q) are bounded and they satisfy

    Property 2[41]: There are two positive scalars ξ1and ξ2such that the friction term satisfies

    Then dynamic equation (9) is written as

    whereDd=τd-F(q˙)-ΔNq¨-ΔBq˙-ΔG.

    Through simplification, the dynamic model (12) is rearranged as

    whereD=N0-1Dddenote a lumped disturbance.Before proceeding, the following assumptions are necessary.

    Assumption 1[6]: The lumped disturbanceDis bounded,and there is a constant δ1>0 that satisfies

    where ηa1, ηa2and ηa3are unknown positive numbers.

    Assumption 2[42]: The first and second derivatives ofDwith respect to timetexist and are bounded, which is described as

    where δ2and δ3are positive scalars.

    Define the desired reference trajectory asqd, and the main purpose of this paper is: 1) Design a suitable sliding mode control scheme to ensure that the robotic manipulator given in(9) can overcome the system uncertainty and external disturbance, and finally accurately track the reference trajectory within a fixed time; 2) Eliminate the influence of chattering on the control performance.

    Remark 1: From the perspective of control engineering,Assumptions 1 and 2 are reasonable.According to Properties 1 and 2, it is not difficult to infer that ΔN, ΔBand ΔGare bounded.In addition, the unknown disturbance τdis usually bounded in practical applications.Therefore, combined with[6], [19] and [20], the above assumptions can be satisfied.

    III.MAIN RESULTS

    A. Design Approach of a Novel FTDO

    In this subsection, a new FTDO is designed as

    B. Design Approach of AFTDO

    In (17), the sliding componentk3sign(η) is used to compensate for the differential term of the lumped disturbanceD˙.Thereinto, the selection of sliding gaink3is inspired by Assumption 2.However, the disturbance boundary δ2is often unknown.In order to overcome this shortcoming, an equivalent controlueqis used to exactly cancel the derivative of the lumped disturbance, which is the average value that must be taken for the switch signal in (17) to maintain sliding.Interpreting this means that it can obtain |ueqi|=|D˙i| during sliding phase.Note that the equivalent control is conceived as a abstract scheme in the analysis of sliding motion.Therefore,the close approximation ofueqis given in real time by using a nonlinear low-pass filter described below.

    Remark 2: The signal value of |uˉeqi| is estimated by a lowpass filter, which enhances the robustness of the adaptive mechanism.Based on (28) and (29),k3i(t) will gradually increase until sliding motion occurs, and then it starts to decrease and eventually converges to the safety margin.By adjusting parameters βaand υa, the safety margin is changed to meet the requirements of different engineering environments for allowable uncertainties, such as [42].

    Remark 3: The designed adaptive scheme described in (28)and (29) is motivated by [42].Note that the nonlinear lowpass filter mentioned in the above literature is also used as an estimator to give a close approximation of the switching signalueqin real time.The original design approach of the mentioned filter was derived from [43].However, the above design can only achieve finite-time stability.Compared with[42], the proposed adaptive scheme can provide fixed-time stability, so that the dynamic system has better convergence performance.

    Remark 4: In fact, the parameter perturbation and external disturbance in industrial robot systems will adversely affect the control performance.To overcome this problem, the disturbance observer was designed in the [31]-[33].It should be noted that the conventional disturbance observer is difficult to estimate the unknown disturbance within a fixed time and the estimation accuracy is not satisfactory.Therefore, a AFTDO is designed in this paper.Compared with [35]-[37], this design can provide higher estimation accuracy.

    C. Novel Design Approach of FTNFTSMC

    The tracking error of the manipulator is defined as ε1=q-qd, whereqdis the reference position vector.The derivative of ε1is defined as ε˙1=ε2.Subsequently, the error system is given by

    To ensure that the tracking error in (39) quickly converges to zero and overcomes the singularity problem, a new NFTSM surface is designed as

    The specific control block diagram is presented in Fig.1.

    Theorem 3: Consider the manipulator system (9) with system uncertainties and external disturbances.If the AFTDO in(17), (28) and (29), the NFTSM surface in (40) and the FTNFTSMC in (43)-(45) are used, the fixed-time arrival of the sliding surface is achieved and then the tracking error converges to zero within a fixed time.

    Proof:The stability analysis of the reaching phase is described below.

    First, a Lyapunov function is constructed as follows:

    Fig.1.Block diagram of FTNFTSM control strategy.

    DifferentiatingVrwith respect to timet, which has

    The stability analysis of the sliding phase is described below.

    Once the state trajectory converges to the NFTSM surfaces=0, it means that the trajectory tracking enters the sliding phase.Then we have

    Defining Θ =arctan(ε1i), it yields

    By calculating the results in (53), one obtains

    Then the Lyapunov function is constructed as

    The time derivative ofVsis

    Therefore, the error system (39) is fixed-time stable by using FTNFTSMC.The complete convergence time is bounded by

    Remark 6: The NFTSM design introduced in (40) is motivated by [17].However, the above control scheme can only achieve finite-time stability, which results in the convergence time of the system is limited by the initial state.Compared with [17], the proposed control scheme can ensure faster convergence speed and achieve global fixed-time stability.

    Remark 7: The singularity problem is a challenge to be overcome in TSM control, which will not only cause serious overharmonization and increase control loss, but also lead to the collapse of the controlled system, resulting in bad effects.In [19], [20], [26], and [28], the piecewise function method was adopted to eliminate this defect.It should be noted that this method requires careful consideration of the continuity and differentiability of nonlinear functions, which increases the complexity of singularity avoidance.Compared with [19],[20], [26], and [28], the NFTSM surface designed in this paper does not need to combine piecewise functions and can conveniently avoid the singularity problem.

    Remark 8: In brief, the advent of robotic manipulators has accelerated the process of automatic production.However,there are also a variety of uncertain factors in the real factory,which can easily deteriorate the control performance.For example, the nonlinear friction between the individual motion mechanisms of robotic manipulators will affect its dynamics model during operation, which may lead to the uncertainty problem in the robot system.In addition, the wear of mechanical parts can easily cause parameter perturbation in the robotic manipulator system, resulting in modeling uncertainty.Therefore, a new AFTDO is designed to estimate the lumped disturbance.This AFTDO can ensure the fixed-time convergence of estimation errors and has the characteristics of fast estimation speed and high estimation accuracy.Simultaneously, a fixedtime NFTSM (FTNFTSM) control scheme is used to enhance the robustness and improve tracking performance.Furthermore, there is no chattering problem in the actual control input signal.Therefore, the proposed FTNFTSM control strategy is feasible, and it is suitable for practical applications.

    IV.SIMULATION RESULTS

    A. Validation of the Proposed FTNFTSMC

    In this section, the 2-link rigid manipulator shown in Fig.2 is used to visually demonstrate the feasibility of the FTNFTSM control strategy.In addition, numerical simulations are completed on MATLAB/Simulink with the numerical method Ode45.The dynamic model is given by

    Fig.2.Architecture of 2-link robotic manipulators.

    TABLE I CONTROL PARAMETERS OF FTNFTSM METHOD

    Figs.3 and 4 show the position tracking and velocity tracking of the reference trajectory by the two links respectively,and the time responses of the tracking errors are presented in Fig.5.It can be seen that even if manipulators are perturbed by system uncertainties and external disturbances, good tracking performance is achieved by using the FTNFTSMC, that is,the actual position and speed of manipulators are quickly adjusted to the desired signal.The trajectories of control torques are plotted in Fig.6, where it shows that the chattering problem has been overcome.

    Fig.7 depicts the response curves of the sliding surface,which implies that the NFTSM variablesis quickly regulated to zero under the action of the FTNFTSMC.The estimation errors of the AFTDO are presented in Fig.8.Interpreting these estimations results means that the AFTDO proposed in this paper can not only achieve rapid estimation of lumped disturbance, but also provide high estimation accuracy.The satisfactory estimation performance also provides sufficient guarantees to achieve better control performance.As shown in Figs.3-8, it confirms that the FTNFTSM control strategy is feasible.In addition, it is worth mentioning that a priori knowledge about the upper bound of the lumped disturbance is not required to be known in advance by using a disturbance observer.

    Fig.3.The trajectories of position tracking.

    Fig.4.The trajectories of velocity tracking.

    Fig.5.The trajectories of tracking errors.

    Fig.6.The control torques of the two links.

    In order to demonstrate that the proposed control scheme has the property of fixed-time convergence, trajectory tracking simulations are performed for five cases with different initial states of the robotic manipulator.In addition, the control parameters are still selected as shown in Table I.The five different initial states are set as follows.The simulation results are shown in Figs.9-13.

    Fig.7.The trajectories of sliding surface.

    Fig.8.The estimation errors of the disturbance observer.

    Case 1:q(0)=[0.25,2.1]T,q˙(0)=[0,0]T;

    Case 2:q(0)=[0.4,2.25]T,q˙(0)=[0.1,0]T;

    Case 3:q(0)=[-0.1,2.3]T,q˙(0)=[0,-0.1]T;

    Case 4:q(0)=[0.2,2.2]T,q˙(0)=[-0.1,0.15]T;

    Case 5:q(0)=[1,1.5]T,q˙(0)=[0,0]T.

    Figs.9-12 show the tracking response of the manipulator on the reference trajectory in five different initial states.Fig.13 depicts the control torque curves generated by the five different cases.From the results in Figs.10 and 12, it can be observed that the convergence time of the trajectory tracking error is almost the same when different initial values are chosen.Together with the stability analysis in Section III, these simulation results confirm that the proposed FTNFTSM control scheme has the capability of fixed-time convergence and the upper bound of the settling time does not depend on the initial state of the robotic manipulator system.

    B. Comparative Results

    To more intuitively demonstrate the excellent tracking performance of the proposed control strategy, the control schemes introduced in [17], [35], [36], [37], and [44] are selected as comparison examples.These control schemes use adaptive techniques, neural networks and disturbance observers to deal with the uncertainty problem in manipulator systems, respectively.Note that all simulation results are obtained under the same simulation conditions to ensure that the comparison results are reliable.In addition, three performance indices introduced in [19] are used to quantitatively analyze the comparison results, that is, the integral of the absolute value of the tracking error (IAE), the integral of the absolute value of the time multiplied by the error (ITAE) and the integral of the square value of the control input (ISV).The corresponding results are presented in Table II.

    Fig.9.Position tracking trajectory based on different initial states.

    Fig.10.Position tracking error based on different initial states.

    Fig.11.Velocity tracking trajectory based on different initial states.

    Fig.12.Velocity tracking error based on different initial states.

    Fig.13.Control torque trajectories based on different initial states.

    TABLE II QUANTITATIVE ANALYSIS OF DIFFERENT CONTROLLERS

    The comparison results of six different control strategies are plotted in Figs.14-18.Figs.14 and 15 show the response curves for position tracking and speed tracking obtained by using six different controllers.The comparative simulation of position tracking errors is shown in Fig.16.It is easy to observe that Link 1 takes 1.3 s to track the reference position trajectory with the FTNFTSMC, while the other controllers take more than 2 s.For Link 2, it takes 1.2 s to track the reference position trajectory with the FTNFTSMC when the other controllers take about 2.5 s.Interpreting this means that the FTNFTSMC can achieve faster tracking speed.The results of quantitative analysis also confirm this advantage.As shown in Table II, the proposed control algorithm can produce smaller IAE and ITAE values.Fig.17 depicts the comparative results of control inputs given by six different controllers.From Fig.17 and its zoom-in graph, we can see that the FTNFTSMC provides a chattering-free control input signal.However, the chattering problem is not adequately addressed when using FTNTSMC (fixed-time NTSMC) and NFTSMC, which will adversely affect the control performance.In addition,compared with the other five controllers, the proposed controller produces lower ISV values in Table II, which implies that FTNFTSMC can guarantee better control input performance.Fig.18 shows the estimation errors of the lumped disturbance given by the four disturbance observers.From the zoomed-in graphs in Fig.18, it is clearly observed that the proposed AFTDO not only provides faster estimation speed,but also achieves higher estimation accuracy compared to the other three disturbance observers, which confirms its superior estimation performance.Therefore, we can conclude that the FTNFTSM control algorithm designed in this paper is feasible and achieves faster tracking speed, higher tracking accuracy and chattering-free control input.

    Fig.15.Comparative results of velocity tracking.

    Fig.16.Comparative results of position tracking errors.

    Fig.17.Comparative results of control torques.

    Fig.18.Comparative results of disturbance estimation errors.

    V.CONCLUSION

    In this paper, a new global FTNFTSM control strategy has been proposed for robotic manipulators perturbed by modeling uncertainties and external disturbances.To eliminate the adverse effects of lumped disturbance, an AFTDO is designed, which has been proved to have better estimation performance and can achieve fixed-time convergence of estimation errors.By combining fixed-time control, a new NFTSM surface based on an arctan function has been developed, which avoids the singularity problem expediently and enhances the robustness of the system.The fixed-time NFTSM controller based on the AFTDO has been designed to ensure that the tracking errors of manipulators can converge to zero within a fixed time.In addition, the setting time is bounded and its upper bound is not constrained by initial system states.The manipulator system is shown to be globally fixed-time stable and the stability analysis has been rigorously proved by Lyapunov method.Finally, the FTNFTSMC is validated on the 2-link manipulator.The simulation results have proved that the FTNFTSMC has excellent control performance.Future research attention will be focused on the fixedtime cooperative control of multiple robotic manipulators, and the issues of packet losses, delays will be addressed.

    亚洲成人免费电影在线观看| 色综合婷婷激情| 狂野欧美激情性xxxx| 超色免费av| 日日干狠狠操夜夜爽| 久久中文字幕人妻熟女| 免费av毛片视频| 久久国产乱子伦精品免费另类| ponron亚洲| 女同久久另类99精品国产91| 不卡av一区二区三区| 免费一级毛片在线播放高清视频 | 亚洲av电影在线进入| 国产一区在线观看成人免费| 国产av在哪里看| 不卡一级毛片| av电影中文网址| 日日夜夜操网爽| 日韩大尺度精品在线看网址 | av片东京热男人的天堂| 欧美日韩视频精品一区| 亚洲一码二码三码区别大吗| 黄频高清免费视频| 国产精品一区二区三区四区久久 | 视频区欧美日本亚洲| 国产欧美日韩一区二区三| 色在线成人网| 久久国产乱子伦精品免费另类| 99国产极品粉嫩在线观看| 精品一区二区三区视频在线观看免费 | 亚洲欧美精品综合久久99| 中文字幕人妻熟女乱码| a级毛片在线看网站| 欧美在线一区亚洲| 夜夜爽天天搞| 欧美日韩亚洲综合一区二区三区_| 黄片小视频在线播放| 男女高潮啪啪啪动态图| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩大码丰满熟妇| 91国产中文字幕| 18禁裸乳无遮挡免费网站照片 | 久99久视频精品免费| 国产精品偷伦视频观看了| 色精品久久人妻99蜜桃| 天堂中文最新版在线下载| 亚洲性夜色夜夜综合| 久久人妻福利社区极品人妻图片| 大陆偷拍与自拍| 国产精品爽爽va在线观看网站 | 日本vs欧美在线观看视频| 在线观看日韩欧美| 国产国语露脸激情在线看| 欧美精品一区二区免费开放| 午夜福利在线观看吧| 国产成人精品无人区| 99精国产麻豆久久婷婷| 国产av精品麻豆| 亚洲精品美女久久久久99蜜臀| 成人免费观看视频高清| 中文字幕人妻丝袜制服| 国产一区二区在线av高清观看| 老汉色av国产亚洲站长工具| 精品国产美女av久久久久小说| 午夜影院日韩av| 乱人伦中国视频| 一进一出抽搐gif免费好疼 | 天天添夜夜摸| 久久中文字幕一级| 丝袜美腿诱惑在线| 午夜福利免费观看在线| 亚洲精品美女久久久久99蜜臀| 黄片小视频在线播放| 免费高清视频大片| 国产黄a三级三级三级人| 久久精品aⅴ一区二区三区四区| avwww免费| 久久人人精品亚洲av| 国产深夜福利视频在线观看| 亚洲 欧美 日韩 在线 免费| 97人妻天天添夜夜摸| 国产亚洲精品综合一区在线观看 | 欧美日韩国产mv在线观看视频| 国产1区2区3区精品| 咕卡用的链子| av国产精品久久久久影院| 丰满迷人的少妇在线观看| 校园春色视频在线观看| 性欧美人与动物交配| 99久久综合精品五月天人人| 女人精品久久久久毛片| x7x7x7水蜜桃| 日本撒尿小便嘘嘘汇集6| 亚洲一区高清亚洲精品| 女警被强在线播放| 97超级碰碰碰精品色视频在线观看| 久99久视频精品免费| 精品国产一区二区久久| 三级毛片av免费| aaaaa片日本免费| 老司机午夜十八禁免费视频| 无人区码免费观看不卡| 制服诱惑二区| 午夜精品久久久久久毛片777| 99久久国产精品久久久| 精品久久蜜臀av无| 法律面前人人平等表现在哪些方面| √禁漫天堂资源中文www| svipshipincom国产片| 国产高清国产精品国产三级| 麻豆一二三区av精品| 欧美乱码精品一区二区三区| 男男h啪啪无遮挡| 曰老女人黄片| 国产亚洲av高清不卡| 亚洲全国av大片| av欧美777| 亚洲,欧美精品.| 少妇被粗大的猛进出69影院| 国产成年人精品一区二区 | 性欧美人与动物交配| 正在播放国产对白刺激| 亚洲精品av麻豆狂野| 精品高清国产在线一区| 热99国产精品久久久久久7| 女同久久另类99精品国产91| 日韩免费av在线播放| 成人国产一区最新在线观看| 午夜福利欧美成人| 99热只有精品国产| 女人被狂操c到高潮| 老司机午夜福利在线观看视频| 在线播放国产精品三级| 国产精品电影一区二区三区| 窝窝影院91人妻| 好看av亚洲va欧美ⅴa在| 国产精品一区二区精品视频观看| 国产精品亚洲一级av第二区| 女人被躁到高潮嗷嗷叫费观| 亚洲国产精品999在线| 日本精品一区二区三区蜜桃| 五月开心婷婷网| 亚洲国产精品一区二区三区在线| 桃红色精品国产亚洲av| 97超级碰碰碰精品色视频在线观看| www.熟女人妻精品国产| 午夜日韩欧美国产| 热99国产精品久久久久久7| 一区二区三区激情视频| 91成人精品电影| 日韩欧美一区二区三区在线观看| 深夜精品福利| 国产精品爽爽va在线观看网站 | 久久久久久久久免费视频了| 一本大道久久a久久精品| 亚洲专区字幕在线| 热re99久久精品国产66热6| 午夜福利欧美成人| 人人妻人人爽人人添夜夜欢视频| 日本黄色日本黄色录像| 久久99一区二区三区| 嫁个100分男人电影在线观看| www国产在线视频色| 丰满人妻熟妇乱又伦精品不卡| 日日夜夜操网爽| 好看av亚洲va欧美ⅴa在| 亚洲精品国产区一区二| 亚洲片人在线观看| 久久久水蜜桃国产精品网| 在线看a的网站| 在线观看日韩欧美| 精品午夜福利视频在线观看一区| 欧美久久黑人一区二区| 国产精品久久久久成人av| 国产一区二区三区视频了| 亚洲第一av免费看| 99在线人妻在线中文字幕| 一夜夜www| 中文字幕人妻熟女乱码| 国产高清videossex| 日韩国内少妇激情av| 亚洲国产欧美日韩在线播放| 亚洲片人在线观看| 精品国产乱码久久久久久男人| 后天国语完整版免费观看| 午夜久久久在线观看| 精品高清国产在线一区| 日韩中文字幕欧美一区二区| 波多野结衣av一区二区av| 国产欧美日韩一区二区三| 久久99一区二区三区| 热99re8久久精品国产| 欧美另类亚洲清纯唯美| 热re99久久国产66热| 国产无遮挡羞羞视频在线观看| 在线视频色国产色| 亚洲av成人av| 久久久久久久久久久久大奶| 久99久视频精品免费| 国产片内射在线| 亚洲国产欧美网| 亚洲人成网站在线播放欧美日韩| 91成人精品电影| 老汉色av国产亚洲站长工具| 1024视频免费在线观看| 老司机午夜福利在线观看视频| 一区二区三区国产精品乱码| 一级黄色大片毛片| 日本 av在线| 免费观看人在逋| 欧美日韩一级在线毛片| 一级片'在线观看视频| 亚洲一区二区三区不卡视频| 不卡一级毛片| 一进一出好大好爽视频| 精品国产乱码久久久久久男人| 国产熟女午夜一区二区三区| www.熟女人妻精品国产| 热99国产精品久久久久久7| 好男人电影高清在线观看| 欧美激情 高清一区二区三区| 国产有黄有色有爽视频| 中文欧美无线码| 黑人猛操日本美女一级片| 高清欧美精品videossex| 韩国精品一区二区三区| 国产91精品成人一区二区三区| 亚洲av第一区精品v没综合| x7x7x7水蜜桃| 久久青草综合色| 亚洲中文字幕日韩| 国产av一区二区精品久久| 国内毛片毛片毛片毛片毛片| 国产极品粉嫩免费观看在线| xxx96com| 国产亚洲av高清不卡| 日韩欧美一区视频在线观看| 国内毛片毛片毛片毛片毛片| 欧美日韩一级在线毛片| 脱女人内裤的视频| 动漫黄色视频在线观看| 一区福利在线观看| 国产亚洲精品久久久久5区| 国产精品秋霞免费鲁丝片| 国产激情欧美一区二区| 久久久久九九精品影院| www国产在线视频色| 日韩免费高清中文字幕av| 女人高潮潮喷娇喘18禁视频| 黄色视频不卡| 亚洲人成电影免费在线| 黄色a级毛片大全视频| 丰满饥渴人妻一区二区三| 久久天堂一区二区三区四区| av超薄肉色丝袜交足视频| 色综合欧美亚洲国产小说| 一本大道久久a久久精品| 国产男靠女视频免费网站| 黄色视频,在线免费观看| 欧美丝袜亚洲另类 | 精品国产一区二区三区四区第35| 俄罗斯特黄特色一大片| 精品熟女少妇八av免费久了| 国产成人免费无遮挡视频| 精品电影一区二区在线| 国产精品久久久av美女十八| 久久国产精品男人的天堂亚洲| 国产高清视频在线播放一区| 中文字幕色久视频| 国产亚洲欧美在线一区二区| 欧美性长视频在线观看| 午夜福利影视在线免费观看| 一二三四在线观看免费中文在| 波多野结衣高清无吗| 视频区欧美日本亚洲| 18禁黄网站禁片午夜丰满| 成人精品一区二区免费| 午夜视频精品福利| 五月开心婷婷网| 精品国产一区二区久久| 丝袜美足系列| 亚洲伊人色综图| 久久久久久大精品| 久久香蕉激情| 免费在线观看日本一区| 亚洲国产欧美日韩在线播放| 日韩三级视频一区二区三区| 久久久久国内视频| 免费av中文字幕在线| 亚洲欧美精品综合一区二区三区| av在线播放免费不卡| 在线看a的网站| 老司机福利观看| 欧美乱色亚洲激情| 婷婷精品国产亚洲av在线| 亚洲人成77777在线视频| 十分钟在线观看高清视频www| 一本综合久久免费| av中文乱码字幕在线| 亚洲九九香蕉| 成年女人毛片免费观看观看9| 曰老女人黄片| 免费av毛片视频| 波多野结衣一区麻豆| 亚洲色图综合在线观看| 激情视频va一区二区三区| 国产成年人精品一区二区 | www日本在线高清视频| 精品一品国产午夜福利视频| 久久久久国内视频| 国产av一区在线观看免费| 精品少妇一区二区三区视频日本电影| 1024香蕉在线观看| 亚洲一区二区三区欧美精品| 人人澡人人妻人| 久久久久久人人人人人| 一个人免费在线观看的高清视频| 97超级碰碰碰精品色视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 精品高清国产在线一区| 99久久人妻综合| 少妇裸体淫交视频免费看高清 | 十八禁人妻一区二区| 久久热在线av| 妹子高潮喷水视频| 99热只有精品国产| 欧美日韩福利视频一区二区| 亚洲国产欧美网| 久久狼人影院| 电影成人av| 欧美大码av| 欧美日韩亚洲高清精品| 黄色视频,在线免费观看| 色综合欧美亚洲国产小说| 精品少妇一区二区三区视频日本电影| 久久精品影院6| 97超级碰碰碰精品色视频在线观看| 不卡av一区二区三区| 1024香蕉在线观看| 国产av又大| 美国免费a级毛片| 精品人妻在线不人妻| 如日韩欧美国产精品一区二区三区| 免费一级毛片在线播放高清视频 | 亚洲 国产 在线| 亚洲av电影在线进入| 91国产中文字幕| 看片在线看免费视频| 久久婷婷成人综合色麻豆| 丰满的人妻完整版| 欧美日韩福利视频一区二区| 好看av亚洲va欧美ⅴa在| 在线永久观看黄色视频| 欧美日韩视频精品一区| 午夜日韩欧美国产| 久久国产精品男人的天堂亚洲| 99riav亚洲国产免费| 一级毛片精品| 亚洲成人精品中文字幕电影 | 国产精品永久免费网站| 丝袜在线中文字幕| 黄色视频,在线免费观看| 成人18禁在线播放| 很黄的视频免费| 欧美人与性动交α欧美软件| 一a级毛片在线观看| 久9热在线精品视频| 女人爽到高潮嗷嗷叫在线视频| 两个人免费观看高清视频| 欧美日韩视频精品一区| 国产亚洲精品一区二区www| 亚洲国产看品久久| 香蕉国产在线看| 亚洲专区字幕在线| 国产精品免费一区二区三区在线| 欧美日韩福利视频一区二区| 热99国产精品久久久久久7| 国产一卡二卡三卡精品| 人人妻人人澡人人看| 亚洲av五月六月丁香网| 啦啦啦 在线观看视频| 久久精品国产99精品国产亚洲性色 | 99国产精品一区二区蜜桃av| 亚洲欧美日韩高清在线视频| 村上凉子中文字幕在线| 麻豆av在线久日| 亚洲av片天天在线观看| 男人舔女人下体高潮全视频| 欧美日韩视频精品一区| 超碰成人久久| 1024香蕉在线观看| e午夜精品久久久久久久| 777久久人妻少妇嫩草av网站| 五月开心婷婷网| 成人18禁在线播放| 亚洲国产欧美日韩在线播放| 午夜亚洲福利在线播放| av天堂久久9| 色在线成人网| 两性夫妻黄色片| 国产精品久久久久成人av| 中文字幕人妻丝袜一区二区| 一级片'在线观看视频| 日韩免费av在线播放| 性色av乱码一区二区三区2| 久久人人爽av亚洲精品天堂| 国产在线观看jvid| 成人手机av| 欧美精品啪啪一区二区三区| 国产精品免费一区二区三区在线| 性少妇av在线| 亚洲精品一卡2卡三卡4卡5卡| 国产成年人精品一区二区 | 波多野结衣高清无吗| 色精品久久人妻99蜜桃| 久久热在线av| 成人三级黄色视频| 99精品久久久久人妻精品| 国产精品一区二区精品视频观看| 欧美日韩中文字幕国产精品一区二区三区 | 三上悠亚av全集在线观看| 母亲3免费完整高清在线观看| 麻豆一二三区av精品| 男女做爰动态图高潮gif福利片 | 久久久国产精品麻豆| 成熟少妇高潮喷水视频| 欧美一区二区精品小视频在线| 交换朋友夫妻互换小说| 麻豆一二三区av精品| 一区二区日韩欧美中文字幕| 欧美日韩一级在线毛片| 亚洲av五月六月丁香网| 国产一区二区激情短视频| 自拍欧美九色日韩亚洲蝌蚪91| 欧美亚洲日本最大视频资源| а√天堂www在线а√下载| 欧美最黄视频在线播放免费 | 男女下面插进去视频免费观看| 丁香六月欧美| 少妇被粗大的猛进出69影院| 热99re8久久精品国产| 国产亚洲av高清不卡| 久久精品国产亚洲av高清一级| 麻豆一二三区av精品| 久久天堂一区二区三区四区| 极品人妻少妇av视频| 黄色毛片三级朝国网站| 精品人妻在线不人妻| 97超级碰碰碰精品色视频在线观看| 精品一区二区三区视频在线观看免费 | 一个人免费在线观看的高清视频| 亚洲国产中文字幕在线视频| 97碰自拍视频| 成人手机av| 777久久人妻少妇嫩草av网站| 亚洲人成77777在线视频| 9色porny在线观看| 亚洲精品久久成人aⅴ小说| 欧美丝袜亚洲另类 | 亚洲av第一区精品v没综合| 脱女人内裤的视频| 亚洲视频免费观看视频| 少妇的丰满在线观看| 纯流量卡能插随身wifi吗| 在线观看午夜福利视频| 久久精品亚洲熟妇少妇任你| 欧美日韩中文字幕国产精品一区二区三区 | 国产主播在线观看一区二区| 国产麻豆69| 一级a爱视频在线免费观看| a级片在线免费高清观看视频| 久久这里只有精品19| 色婷婷久久久亚洲欧美| av视频免费观看在线观看| 免费在线观看亚洲国产| 999精品在线视频| 久久精品成人免费网站| 99国产精品一区二区三区| 国产免费现黄频在线看| 一级毛片女人18水好多| 亚洲精品一二三| 午夜成年电影在线免费观看| 在线视频色国产色| 日本wwww免费看| 欧美亚洲日本最大视频资源| 丝袜人妻中文字幕| a级毛片黄视频| 俄罗斯特黄特色一大片| 看黄色毛片网站| 啦啦啦 在线观看视频| 精品第一国产精品| 在线观看免费视频网站a站| 少妇的丰满在线观看| 欧美精品啪啪一区二区三区| 十八禁网站免费在线| 亚洲精品国产色婷婷电影| 在线观看一区二区三区| 亚洲avbb在线观看| 日韩精品免费视频一区二区三区| 91麻豆av在线| 国产又色又爽无遮挡免费看| 欧洲精品卡2卡3卡4卡5卡区| 午夜免费观看网址| www.熟女人妻精品国产| 午夜福利一区二区在线看| 欧美大码av| 欧美老熟妇乱子伦牲交| 欧美久久黑人一区二区| 国产精品一区二区三区四区久久 | 免费看a级黄色片| 国产国语露脸激情在线看| 国产精品综合久久久久久久免费 | 亚洲人成77777在线视频| 丁香六月欧美| 精品国产美女av久久久久小说| 熟女少妇亚洲综合色aaa.| 国产野战对白在线观看| 两人在一起打扑克的视频| 亚洲av成人不卡在线观看播放网| 欧美成狂野欧美在线观看| 日本黄色视频三级网站网址| www.www免费av| 一二三四在线观看免费中文在| 亚洲av熟女| 无限看片的www在线观看| 黑丝袜美女国产一区| 亚洲国产欧美网| 90打野战视频偷拍视频| 麻豆一二三区av精品| 欧美激情极品国产一区二区三区| 欧美最黄视频在线播放免费 | 男人的好看免费观看在线视频 | 在线观看一区二区三区激情| aaaaa片日本免费| 亚洲伊人色综图| 在线av久久热| 久久欧美精品欧美久久欧美| av福利片在线| 日韩免费av在线播放| 在线观看日韩欧美| 麻豆久久精品国产亚洲av | 亚洲精品成人av观看孕妇| 久久99一区二区三区| 欧美激情久久久久久爽电影 | 69av精品久久久久久| 欧美国产精品va在线观看不卡| 欧美激情高清一区二区三区| 亚洲一区二区三区不卡视频| 国产有黄有色有爽视频| 美国免费a级毛片| 黑丝袜美女国产一区| 91在线观看av| 午夜福利在线免费观看网站| 黄色毛片三级朝国网站| 琪琪午夜伦伦电影理论片6080| 久久 成人 亚洲| 在线天堂中文资源库| 精品国内亚洲2022精品成人| 日本黄色日本黄色录像| 欧美最黄视频在线播放免费 | 制服诱惑二区| 成人av一区二区三区在线看| 三级毛片av免费| 亚洲精品在线观看二区| 国产亚洲欧美98| 亚洲国产精品999在线| www国产在线视频色| 国产不卡一卡二| 多毛熟女@视频| 在线永久观看黄色视频| 在线播放国产精品三级| 亚洲,欧美精品.| 夜夜躁狠狠躁天天躁| 天堂√8在线中文| 日日摸夜夜添夜夜添小说| 五月开心婷婷网| 精品国产国语对白av| 日韩欧美免费精品| 久久久久国产精品人妻aⅴ院| 一个人观看的视频www高清免费观看 | 黑人巨大精品欧美一区二区mp4| 亚洲国产欧美一区二区综合| 一边摸一边抽搐一进一出视频| 亚洲五月天丁香| 亚洲熟女毛片儿| 欧美色视频一区免费| 大型黄色视频在线免费观看| 免费不卡黄色视频| 美女大奶头视频| √禁漫天堂资源中文www| 精品国产美女av久久久久小说| 9色porny在线观看| 欧美激情极品国产一区二区三区| 在线观看免费高清a一片| 精品久久久久久久毛片微露脸| av电影中文网址| 人成视频在线观看免费观看| 一级a爱视频在线免费观看| 日韩大码丰满熟妇| 老汉色∧v一级毛片| 黄色片一级片一级黄色片| 涩涩av久久男人的天堂| 亚洲第一欧美日韩一区二区三区| 精品人妻在线不人妻| 久久人人爽av亚洲精品天堂| 女生性感内裤真人,穿戴方法视频| 99国产精品99久久久久| 香蕉丝袜av| 欧美日韩av久久| 久久午夜亚洲精品久久| 亚洲全国av大片| 亚洲成国产人片在线观看| 美女午夜性视频免费| 欧美大码av| 日韩三级视频一区二区三区|