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

    Adaptive Trajectory Tracking Control for Nonholonomic Wheeled Mobile Robots:A Barrier Function Sliding Mode Approach

    2024-04-15 09:37:30YunjunZhengJinchuanZhengKeShaoHanZhaoHaoXieandHaiWang
    IEEE/CAA Journal of Automatica Sinica 2024年4期

    Yunjun Zheng , Jinchuan Zheng ,,, Ke Shao ,,, Han Zhao ,Hao Xie , and Hai Wang ,,

    Abstract—The trajectory tracking control performance of nonholonomic wheeled mobile robots (NWMRs) is subject to nonholonomic constraints, system uncertainties, and external disturbances.This paper proposes a barrier function-based adaptive sliding mode control (BFASMC) method to provide high-precision, fast-response performance and robustness for NWMRs.Compared with the conventional adaptive sliding mode control,the proposed control strategy can guarantee that the sliding mode variables converge to a predefined neighborhood of origin with a predefined reaching time independent of the prior knowledge of the uncertainties and disturbances bounds.Another advantage of the proposed algorithm is that the control gains can be adaptively adjusted to follow the disturbances amplitudes thanks to the barrier function.The benefit is that the overestimation of control gain can be eliminated, resulting in chattering reduction.Moreover, a modified barrier function-like control gain is employed to prevent the input saturation problem due to the physical limit of the actuator.The stability analysis and comparative experiments demonstrate that the proposed BFASMC can ensure the prespecified convergence performance of the NWMR system output variables and strong robustness against uncertainties/disturbances.

    I.INTRODUCTION

    IN recent decades, the control design of high-order nonlinear multiple input multiple output (MIMO) system has attracted much attention due to their wide applications in many fields, such as autonomous vehicles, unmanned aerial vehicles and intelligent robots [1].As typical high-order nonlinear MIMO systems, the nonholonomic wheeled mobile robots (NWMRs) play a vital role in environment exploration,industrial transportation and intelligent inspection and other fields owing to the characteristics of fast speed, high energy utilization, high interactivity and integration [2].However, the control performance of NWMRs may be considerably affected by time-varying uncertainties, such as parameter variation,external disturbance, and unmodeled dynamics.Besides,another main challenge faced by the NWMRs control design is the inevitable nonholonomic constraint which is not integrable and can not be eliminated in the mathematical model.Fortunately, in recent years, researchers have designed a variety of robust or adaptive control methods for NWMRs to deal with the time-varying uncertainties and nonholonomic constraint, such as sliding mode control (SMC) [3], [4], backstepping control [5], [6], model predictive control [7], [8], fuzzy control [9], learning control [10], [11] and adaptive control[12], [13], etc.Although these control methods achieve some outstanding performance for NWMRs, for instance, fast convergence rate, high accuracy, and robustness against uncertainties, they also have some drawbacks.For example, only asymptotic stability is guaranteed; estimating the finite time of error convergence is complicated; prior knowledge of uncertainties is employed, and a tremendous amount of computing is required.

    Due to the excellent robustness, fast response and good transient performance, the SMC is considered as an effective tool for nonlinear controller design, which has been widely used in robot systems, such as robotic manipulators [14], actuated exoskeletons [15], bioinspired robots [16], and mobile robots [17].However, the SMC is not flawless.Once a more considerable switching gain is selected to suppress the unknown uncertainties and disturbances, the severe control chattering will emerge, even leading to the actuator saturation and energy loss.Furthermore, if uncertainties and disturbances increase to some extent that exceeds the ability of the selected switching gain, then the sliding mode will be lost, and errors will increase accordingly.The control chattering in mechatronic systems may cause severe deterioration resulting in undesired overshoots, oscillations, and instabilities.To alleviate the control chattering, several techniques have been developed, such as boundary layer technique (BLT) [18], filter technique [19], high-order sliding mode control (HSMC)[20], and adaptive SMC (ASMC) [21].The BLT is a classic approach to eliminate chattering by adopting a boundary layer around the switching surface, yet it compromises performance and chattering.Besides, a prior knowledge of the uncertainties bounds is still needed to ensure the stability of SMC.The filter technique alleviates the control chattering by softening the discontinuous control input and generating a continuous equivalent control input [22].However, similarly to BLT, it also requires the information of uncertainties bounds, and the control performance is sensitive to filter time constants.The HSMC solves the chattering problem by redesigning the sliding mode surface.In [23], a fast nonsingular terminal sliding mode (FNTSM) is proposed to achieve the tracking control of a linear motor with a finite-time zero-convergence property.In FNTSM, the switching element is replaced by a continuous function, which makes the FNTSM control (FNTSMC) chattering-free in nature [24].However,the converging neighborhood size relies on the prior knowledge of the uncertainties bounds.In ASMC, an adaptive control algorithm (ACA) is introduced to adapt the switching gains with respect to uncertainties and disturbances.Hence,the relatively small control gains are generated to suppress uncertainties and disturbances, which ensures the SMC with little control chattering.Nevertheless, the overestimation issue may occur when ACA estimates the unknown uncertainties and disturbances.

    Since the ASMC method combines the attractive properties of both SMC and ACA, it has attracted a great deal of researchers’ attention recently, and some interesting research results have been applied to the trajectory-tracking control of NWMRs.For example, an extended state observer-based ASMC is proposed in [25] to estimate the aggregate uncertainty of NWMRs online with a PID-type sliding surface.Nevertheless, it can only guarantee asymptotic trajectory convergence for NWMRs as time tends to infinity.In [26], a fuzzy integral ASMC is employed to compensate for the lumped uncertainty of NWMRs with fuzzy logic.However,the complicated structure using membership functions and rules must be evaluated, resulting in a higher computational load for execution.A neural dynamic model-based ASMC for NWMRs is proposed in [27] to ensure the global asymptotic stability of the tracking errors without the prior knowledge of uncertainties bounds.Although the neural network scheme significantly alleviates the control chattering in SMC, the complexity issue of the controller and learning algorithm severely limits its practical application.In [28], an adaptive FNTSMC is developed to achieve the finite-time stabilization of the trajectory tracking errors for NWMRs.However, the tracking errors under this method can only be guaranteed with uniformly ultimate boundedness.The ASMCs mentioned above share the idea of using ACA to estimate the lumped uncertainty and applying SMC to counteract time-varying uncertainties and disturbances.Furthermore, some things could be improved for these controllers in reality.First, the ACAs in these controllers are designed in a gain-increasing form.More specifically, the control gains increase monotonically to be sufficient for solving the reaching phase and are locked at fixed values.When the amplitudes of uncertainties and disturbances become small, the control gains can not follow the variations.Thus, the control gains may be overestimated, which results in larger chattering and control magnitude.Second, the settling time and convergence domain are unpredictable, though some methods theoretically guarantee finite-time zero-convergence of the tracking errors.Moreover,the convergence properties rely on the initial conditions and the knowledge of the lumped uncertainty, which is unknown in reality.

    Several solutions have already been proposed to solve the overestimation problem in ASMC with the first/second order sliding mode.Talebet al.[29], [30] and Utkin and APoznyak[31] have designed the first/second-order ASMC for the uncertain MIMO nonlinear systems with time-varying control gains defined in piecewise function form.The piecewise function-based ASMC (PFASMC) algorithms allow control gains to increase (or decrease) versus the variations of the uncertainties and disturbances without overestimating the switching gain and with guaranteeing the finite-time convergence of sliding mode variables to a neighborhood of zero.However,these approaches could be further enhanced.The adaptive schemes of these controllers have only been applied to the first/second order sliding mode, namely, no results of the adaptive scheme for high-order sliding mode.Besides, the size of the convergence domain and time rely on the bounds of the uncertainties and perturbations.Recently, the barrier function(BF) has been applied for the first-order SMC as a novel modality of adapting control gains [32], [33].In [33], a BFbased ACA makes the control gains follow the absolute values of uncertainties and drives the sliding mode variables to reach a prescribed neighborhood of zero in a finite time.Nevertheless, it is also imperfect that when the sliding mode variables approach the zone boundary (i.e., the vertical asymptote of BF), the control inputs step to infinity, leading to an actuator saturation problem.

    Motivated by [29], [33], we propose a BF-based ASMC(BFASMC) strategy for the trajectory tracking control of NWMRs subject to uncertainties and perturbations.By adopting the BFASMC, the chattering and overestimation problem can be effectively mitigated without compromising the control performance, and the actuator saturation problem can also be avoided.The comparative experiments are conducted on a real NWMR platform between the BFASMC and compared controllers, i.e., FNTSMC, PFASMC and adaptive sliding mode dynamic control (ASMDC) [34].The experimental results demonstrate the superior performance of the proposed BFASMC in terms of faster convergence rate, higher tracking accuracy and better robustness.The main contributions of this paper are as follows:

    1) The dynamic model of the NWMR formulated with uncertainties and disturbances is presented with novel error dynamics to serve the BFASMC design.Owing to the advantages of the novel error dynamics, the under-actuated problem and singularity problem are solved for the NWMR control system.

    2) A BFASMC is presented to ensure that the NWMR tracking errors fast converge to a predefined neighborhood of zero with a predefined reaching time despite uncertainties and disturbances.Thanks to the advantages of BF, the requirements for the uncertainties and disturbances bounds and the overestimation of control gain are eliminated during the reaching phase and adaptive phase.

    3) A modified BF-like control gain is designed to prevent the actuator saturation problem and retain the specified transient and steady-state performance, i.e., predefined reaching time and error bound, for the NWMR.

    Notations: Fors,f∈Rm,|s| d enotes [|s1|,...,|sm|]T,s>f(ors≤f) denotes all elements ofsare greater than (less than or equal to) those off.For the sake of brevity, the arguments of system functions may be omitted in the following sections.

    II.PLANT MODELING

    A typical NWMR platform applied in our experiment is shown in Fig.1(a).The simplified model of the NWMR platform is shown in Fig.1(b).In this model, suppose the mass center coordinate of the NWMR is located at the mid-point of the connection between the driving wheels, i.e.,Oxy.The position of the NWMR in the reference coordinateOXYis defined as χ=[x,y,θ]T, whereθis the angle between the forward direction andX-axis.vandwdenote the forward velocity and angular velocity of NWMR, respectively.τland τrdenote the torques of the left and right driving wheels generated by DC motors, respectively.Besides, 2bdenotes the distance between two driving wheels, andris the radius of the driving wheel.In this work, the anti-clockwise direction is defined as the positive direction.The kinematic model and dynamic model of NWMR are given by [35]

    Due to the pure rolling and non-slipping assumptions of NWMR, a nonholonomic constraint is imposed on it, i.e.,Aχ˙=0.Taking the time derivative of (1), it yields

    Substituting (3) into the dynamic model (2), we obtain

    Multiplying both sides of (4) byHTyields

    Obviously, theHsatisfies the following equations:

    By (5) and (6), the dynamic model (2) can be rewritten as

    III.TRAJECTORY ERROR DYNAMICS

    To formulate the trajectory tracking control problem, a reference trajectory for the NWMR is defined as follows:

    Remark 1: The novelties of new error dynamics (17) are that it translates the under-actuated original dynamical model(14) into a fully-actuated one to facilitate the controller design, and avoids the singularity problem of the controller law which will be discussed in Remark 5.

    Assumption 1: Suppose the lumped uncertaintydis timevarying and unknown, but the bound ofdexists under the following property:

    IV.BFASMC DESIGN

    This section proposes a BFASMC strategy for the uncertain NWMR system to achieve accurate position tracking with predefined reaching time and prescribed precision.

    A. Preliminaries of BF

    Fig.2.Schematic diagram of BF.

    B. Controller Design

    First, design a nonsingular terminal sliding manifold(NTSM) [23] as

    C. Time Scale Transformation

    D. Stability Analysis

    Remark 8: The results of Theorem 1 introduce two features of the proposed controller.First, the designed control gain matrix ?Aguarantees the finite-time stability of the sliding dynamics and error dynamics.It ensures the system trajectory converges to a predefined region with a prescribed reaching time, which suits many mobile systems that demand specified tracking accuracy and settling time.This feature makes the BFASMC different from the control method of [39], which only guarantees the time-varying bound of system trajectory with indistinct reaching time.

    V.MODIFICATION AND IMPROVEMENT OF BFASMC

    Based on Assumption 3, it can be inferred that σ1,2satisfy the finite-time stabilization condition [38].Hence, theσ1,2will converge to the region (-hi,hi) in a finite time.Furthermore, it can be obtained that

    Remark 9: The foregoing discussion indicates that the BFASMC (27) with the ?A(22) and the modified ?B(64) possesses the capability of enforcing the σ1,2and ζ1,2to converge to region {|σ1,2|≤ρ1,2<η} and {|ζ1,2|≤2ρ1,2<2η},respectively.Even if σ1,2leave the region { |σ1,2|≤ρ1,2} due to severe disturbances, the modified BF-like control gain ?Bican still guarantee the finite time stability and prevent the actuator saturation problem.

    Remark 10: From the above analysis, the designed control scheme (27) with the adaptive strategy (22)-(26) and the modified BF-like control gain (64) has the following advantages: i) removal of the requirement for uncertainties and disturbances bounds during the control design, ii) elimination of overestimation of adaptive control gain leading to less control chattering and vibrations, iii) explicitly predefined bounds for the tracking error and reaching time, and iv) prevention of the input actuator saturation.

    Furthermore, the designed control scheme is worthy of going into further development.For example, the convergence time and tracking accuracy can be further optimized using the heuristic optimization approach [41]; the undesirable chattering in control input could be alleviated by integrating fuzzy control scheme.

    VI.PARAMETERS DESIGN OF BFASMC

    The anticipated control performance in the practical application involves accuracy, response speed, robustness, signal chattering, and control effort.In order to acquire a more balanced system performance, this paper provides a design guideline of the BFASMC parameters, which also facilitates the usage of this control method.

    1)Selections ofλ1,λ2andμ1,μ2: From the sliding mode function (19), we know that λ1,2and μ1,2critically dominate the sliding mode dynamics performance, such as convergence rate, tracking accuracy and control chattering.The smallerλ1,2and larger μ1,2lead to a faster convergence rate and higher tracking accuracy ofebut at the cost of amplifying the measurement noise and control chattering.In experiment,λ1=4.1, λ2=4.2, μ1=1.9, μ2=1.9.

    VII.EXPERIMENTAL RESEARCH

    To demonstrate the capability of the proposed BFASMC,we implement the verification experiments on a real NWMR test platform displayed in Figs.1(a) and 4.The NWMR system contains a main body supported by two driving wheels,two auxiliary wheels and some essential electrical devices,including a real-time microcontroller, a Li-polymer battery,two motor drivers (Maxon ESCON-36/2), and two DC motors(Maxon 118391) (equipped with encoder to measure the rotary angle) driving wheels to achieve autonomous motion.The controller (27) is first programmed in MATLAB 2021a environment to verify the validity with the numerical simulation of the NWMR system.Afterward, the proposed controller is programmed in the LabVIEW software and implemented on a real-time NI microcontroller with a sampling time of 10 ms to drive the real NWMR.The system inputs,i.e., the posture signals of the NWMR, are calculated using a posture estimator with the motor encoder signals.The nominal parameters of the NWMR arem0= 45 kg,I0=0.56 kg·m2,b=0.188m,r=0.042m.Since the payload is the main parameter uncertainty for the NWMR system, a3 kg payload is fixed at the NWMR platform to test the robustness against parameter uncertainty.The recently presented control methods, i.e., FNTSMC [23], PFASMC [30], and a widely adopted NWMR control strategy, i.e., ASMDC [34], are also carried out as the compared controllers.

    A. Comparative Controllers

    Fig.4.Test platform.

    For a fair comparison, the sliding manifolds for PFASMC and FNTSMC are the same as that in BFASMC (27), which results in the equivalent control inputs of PFASMC and FNTSMC are consistent with the τ0(20).This subsection introduces only the reaching control inputs of PFASMC and FNTSMC.All parameter values of the compared controllers are carefully tuned to achieve the ideal tracking performance for the NWMR system.

    1)PFASMC: The reaching control of PFASMC is

    B. Experimental Results and Analyses

    To illustrate the performance metrics of the proposed BFASMC for NWMRs, we conduct two experiments with respective scenarios, i.e., with initial position offset and external shock disturbance.Besides, the NWMR platform is mounted with a 3 kg payload for all experiments.A circular reference trajectory for the NWMR is designed as follows:

    namely, the radius is 0.6 m and angular speed is 0.7 r ad/s.The initial posture is set asx(0)=-0.1,y(0)=-0.1, θ(0)=0,which makes the initial conditions of σ1,2satisfy|σ1,2(0)|≥η/2 , and experiment duration is set as 20 s.In order to analyze the experimental results more accurately, it is important to define some nomenclatures, such as the root-mean-square of error (RMSe) , the maximum value of error (MAXe) and accumulative error (Ae).The mathematical forms of these nomenclatures are as follows:

    wherendenotes the sample number,eidenotes theith sampled error and ?tis the sampling period.The statistical results ofRMSe,MAXe, andAeare summarized by Table I in which some superior values are highlighted in boldface.

    1)Trajectory Tracking Control With Initial Position Offset:This experiment verifies the tracking performance and robustness of the controllers with an initial position offset and payload.In Fig.5, it can be found that the outputs of the sliding mode variables σ1,2under BFASMC meet the desired tracking performance illustrated in Theorem 1.Specifically, the σ1,2converge to the set {|σ1,2(t)|<0.025} before the time upper boundTp=1.5 s, and the trajectories of σ1,2are contained in the set {|σ1,2(t)|<0.05} for all future timet≥1.5 s.The tracking results of position profiles and tracking errors under the three controllers are shown in Figs.6 and 7.It can be found that the BFASMC has the fastest convergence rate,and the settling times of tracking errors under BFASMC are 0.95 s, 1.19 s and 1.35 s which are 185.2%,583.3%,551.3%;213.7%,712.4%,541.2% and 375.7%,273.9%,167.4% faster than those of FNTSMC, PFASMC and ASMDC, respectively.The rapid convergence rate of BFASMC can be attributed to the sliding mode function (19), the time-varying gain φ (t) (25)and the proportional termk? with the adaptive laws (24).Besides, the steady-state errors (SSEs) under BFASMC are bounded in 0.0045 m, 0.0015 m, and 0.006 rad, which are smaller than those of FNTSMC, PFASMC and ASMDC.This is because the adaptive algorithm in BFASMC can estimate the uncertainties and disturbances bounds and avoid overestimating based on the BF structure.A similar conclusion can also be observed in Table I.TheRMSe,MAXeandAevalues of BFASMC are smaller than those of FNTSMC, PFASMC and ASMDC.Fig.8 illustrates the voltage inputs of driving wheels, which reveal the control input levels of the NWMR under the four controllers.It can be observed that the PFASMC and ASMDC have severe chattering instead of the less chattering in BFASMC and FNTSMC.The reason lies in the switching elements in the adaptive laws and reaching controls of PFASMC and ASMDC.In addition, the chattering of FNTSMC is smaller than that of BFASMC.This is because the usage of continuous function makes the FNTSMC chattering-free theoretically; however, the chattering due to noise is inevitable in experiment.

    TABLE I TRACKING PERFORMANCE EVALUATION OF THE TWO EXPERIMENTS

    Fig.5.The trajectories of the sliding mode variables σ1,2 under BFASMC with an initial position offset.

    Fig.6.Comparative results of the trajectory tracking with an initial position offset.

    Fig.7.Comparative results of the tracking error e with an initial position offset.

    Fig.8.Comparative results of the control input with an initial position offset.

    2)Trajectory Tracking Control With External Shock Disturbance: To verify the robustness performance of the proposed controller against the external disturbance, two half-sine waveform electronic disturbances δ1=δ2=5|sin(2πt)| are imposed in the control inputτ.The disturbances duration is 0.5 s, occurring at the 12th s when the NWMR system is stable.The trajectories of the sliding mode variables σ1,2under BFASMC are shown in Fig.9.It can be seen that the tracking responses of σ1,2escape from the steady state when the external shock disturbances occur.Since the control gain matrixΦ has been switched to the BF-based term ?B, the control inputs increase rapidly to pull back the σ1,2to the origin due to the characteristics of BF.After the |σ1,2| decrease, the control inputs decrease accordingly to avoid overestimating control gains.This phenomenon verifies the performance features of the design adaptive mechanism and BF.From Figs.10 and 11,it can be found that the BFASMC takes a shorter time to make the NWMR return to steady state operation and retains the tracking errors in a desired level, i.e.,±0.0045 m,±0.0015 m,±0.006rad.Nevertheless, FNTSMC, PFASMC and ASMDC all have a significant degradation of the trajectory tracking and error tracking.The statistical results in Table I also confirm the effectiveness of the BFASMC to counteract the external shock disturbance.For instance, the increment ofRMSeunder BFASMC in this experimental case is smaller than those of the other controllers.The control signals illustrated in Fig.12 also verify the feature of BF, i.e., rapidly changing with respect to the external shock disturbance.In addition, it can be seen that the peak values of the control inputs after the shock disturbance under BFASMC are 3.17977 V and 3.17996V, which are 7 7.4%,75.9%; 7 4.6%,76.8% and85.6%,81.6%smaller than those of FNTSMC, PFASMC and ASMDC, respectively.This phenomenon demonstrates that the modified BF-like control gain can effectively prevent the actuator saturation; meanwhile, the specified tracking precision is guaranteed.

    Fig.9.The trajectories of the sliding mode variables σ1,2 under BFASMC with external shock disturbance.

    Fig.10.Comparative results of the trajectory tracking with external shock disturbance.

    Fig.11.Comparative results of the tracking error e with external shock disturbance.

    Fig.12.Comparative results of the control input with external shock disturbance.

    VIII.CONCLUSION

    This paper developed a BFASMC strategy for the trajectory tracking control of NWMRs in situations of external disturbances and parameter uncertainties.The proposed approach differs from the existing ASMC methods, which apply the monotonically increasing adaptive mechanism to suppress uncertainties and disturbances.The usage of BF in the BFASMC guarantees that the control gains can be adaptively adjusted to follow the change of disturbances, i.e., control gains increase (decrease) as disturbances increase (decrease),resulting in the elimination of gain overestimation.On the other hand, the suggested approach ensures that the sliding mode variables converge to a prescribed neighborhood of zero in a predefined time, independent of the initial condition and uncertainties bound.As a result, the NWMR system output variables, i.e., tracking errors, can converge to a prescribed region near zero with a shorter settling time.Besides, the developed modified BF-like control gain resolves the actuator saturation problem for the existing BF-based control methods.Finally, the comparative experimental results on the NWMR platform confirm the superior properties of the developed BFASMC.

    The future work of our research will include further improving the control performance of the proposed BFASMC by integrating with the heuristic optimization approach; integrating the proposed BFASMC with fuzzy control scheme to alleviate the chattering in control input.

    APPENDIX A

    APPENDIX B

    国产伦一二天堂av在线观看| 国产伦精品一区二区三区视频9| 搡女人真爽免费视频火全软件| 久久综合国产亚洲精品| 一级av片app| 中国美白少妇内射xxxbb| 色哟哟哟哟哟哟| 亚洲真实伦在线观看| 精品一区二区免费观看| 亚洲欧洲日产国产| 国产黄色小视频在线观看| 亚洲欧美清纯卡通| 最近最新中文字幕大全电影3| 国语自产精品视频在线第100页| 婷婷色综合大香蕉| av卡一久久| 成熟少妇高潮喷水视频| 久久久久国产网址| 婷婷六月久久综合丁香| 最近手机中文字幕大全| 黄色欧美视频在线观看| 日韩一本色道免费dvd| 九色成人免费人妻av| 99久国产av精品国产电影| 午夜久久久久精精品| .国产精品久久| 国内揄拍国产精品人妻在线| 夜夜爽天天搞| 国产91av在线免费观看| 三级男女做爰猛烈吃奶摸视频| avwww免费| 日日干狠狠操夜夜爽| 国产精品伦人一区二区| 天天一区二区日本电影三级| 在线观看一区二区三区| 午夜福利在线观看吧| 亚洲美女搞黄在线观看| 免费观看在线日韩| 国产亚洲91精品色在线| 男插女下体视频免费在线播放| 欧美另类亚洲清纯唯美| 一级黄色大片毛片| 久久人人爽人人爽人人片va| 99热6这里只有精品| 国产av麻豆久久久久久久| 国产免费一级a男人的天堂| 亚洲熟妇中文字幕五十中出| 三级毛片av免费| 日韩视频在线欧美| av在线蜜桃| АⅤ资源中文在线天堂| 国产日本99.免费观看| 久久亚洲国产成人精品v| 亚洲精品乱码久久久久久按摩| 国产一区二区激情短视频| 精品一区二区三区视频在线| 日韩成人伦理影院| 黄片wwwwww| 国产黄色小视频在线观看| 最近2019中文字幕mv第一页| 欧美人与善性xxx| 亚洲自拍偷在线| 一个人免费在线观看电影| 美女高潮的动态| 国产伦在线观看视频一区| 亚洲成人精品中文字幕电影| 日韩成人av中文字幕在线观看| 干丝袜人妻中文字幕| 欧美人与善性xxx| 久久国产乱子免费精品| 成人无遮挡网站| 高清在线视频一区二区三区 | 黄色日韩在线| 亚洲国产欧洲综合997久久,| 免费在线观看成人毛片| 国产伦在线观看视频一区| 深夜a级毛片| а√天堂www在线а√下载| 女人被狂操c到高潮| 国产精品乱码一区二三区的特点| 美女国产视频在线观看| 最后的刺客免费高清国语| 99在线视频只有这里精品首页| 日韩av不卡免费在线播放| av在线蜜桃| 亚洲欧美精品专区久久| 一边亲一边摸免费视频| 中文字幕熟女人妻在线| 国产精品爽爽va在线观看网站| or卡值多少钱| 亚洲av熟女| 国产老妇伦熟女老妇高清| 国产单亲对白刺激| 天堂网av新在线| 日日摸夜夜添夜夜添av毛片| 一个人看的www免费观看视频| 一进一出抽搐gif免费好疼| 国内揄拍国产精品人妻在线| 国产精品三级大全| 色播亚洲综合网| 国产精品久久久久久精品电影小说 | 一本久久精品| 99久久精品一区二区三区| 国产成人精品久久久久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美色视频一区免费| 国产在线精品亚洲第一网站| 日本五十路高清| 精品免费久久久久久久清纯| 久久久a久久爽久久v久久| 国产亚洲av片在线观看秒播厂 | 又爽又黄无遮挡网站| 天天一区二区日本电影三级| 美女黄网站色视频| 简卡轻食公司| 九草在线视频观看| 99久久中文字幕三级久久日本| 欧美色视频一区免费| 99久久精品一区二区三区| 亚洲精品久久久久久婷婷小说 | 日本一二三区视频观看| 中国美白少妇内射xxxbb| 精品久久久久久久久亚洲| 性插视频无遮挡在线免费观看| 观看免费一级毛片| 在线观看av片永久免费下载| 舔av片在线| 99久国产av精品| 亚洲经典国产精华液单| 欧美性猛交黑人性爽| 性欧美人与动物交配| 色综合色国产| 国产精品三级大全| 亚洲美女搞黄在线观看| 在线观看av片永久免费下载| 在线观看美女被高潮喷水网站| 亚洲久久久久久中文字幕| 高清午夜精品一区二区三区 | 国产精品人妻久久久久久| 欧美日韩在线观看h| 毛片一级片免费看久久久久| 小蜜桃在线观看免费完整版高清| 深夜a级毛片| 美女内射精品一级片tv| 国产色婷婷99| 3wmmmm亚洲av在线观看| 简卡轻食公司| 亚洲经典国产精华液单| 综合色av麻豆| 久久人妻av系列| 国产乱人视频| 日韩成人av中文字幕在线观看| 亚洲不卡免费看| 午夜福利在线观看免费完整高清在 | 国产美女午夜福利| h日本视频在线播放| 国产男人的电影天堂91| 亚洲av电影不卡..在线观看| 亚洲欧美精品自产自拍| 狠狠狠狠99中文字幕| 老师上课跳d突然被开到最大视频| 国产女主播在线喷水免费视频网站 | av.在线天堂| 国产精品爽爽va在线观看网站| 蜜桃亚洲精品一区二区三区| 12—13女人毛片做爰片一| av在线天堂中文字幕| 国产久久久一区二区三区| 午夜福利在线观看免费完整高清在 | 人妻少妇偷人精品九色| 精品少妇黑人巨大在线播放 | 日韩高清综合在线| 91狼人影院| 丰满的人妻完整版| 一进一出抽搐动态| 欧美不卡视频在线免费观看| 国产亚洲91精品色在线| 性欧美人与动物交配| 国产亚洲精品久久久久久毛片| 成人特级av手机在线观看| 成人亚洲精品av一区二区| 免费观看在线日韩| a级毛片a级免费在线| 亚洲人成网站在线播| 两个人的视频大全免费| 狠狠狠狠99中文字幕| 国产精品麻豆人妻色哟哟久久 | 在线天堂最新版资源| 九九热线精品视视频播放| 99热这里只有精品一区| 日本五十路高清| 亚洲精品影视一区二区三区av| 成人欧美大片| 亚洲av.av天堂| 能在线免费观看的黄片| 国产精品国产高清国产av| 亚洲熟妇中文字幕五十中出| 蜜桃久久精品国产亚洲av| 免费观看在线日韩| 久久综合国产亚洲精品| 99热精品在线国产| 久久久久久久久大av| 日韩 亚洲 欧美在线| 欧美激情久久久久久爽电影| 淫秽高清视频在线观看| 日本一本二区三区精品| 成人国产麻豆网| 可以在线观看的亚洲视频| 黄色一级大片看看| 噜噜噜噜噜久久久久久91| 国内久久婷婷六月综合欲色啪| 欧美+亚洲+日韩+国产| 毛片一级片免费看久久久久| 麻豆成人av视频| 黑人高潮一二区| 婷婷色av中文字幕| 午夜福利在线观看免费完整高清在 | 国产单亲对白刺激| 久久久久久久久久成人| 美女大奶头视频| 我要搜黄色片| 成人二区视频| 91狼人影院| 99久久精品热视频| 给我免费播放毛片高清在线观看| 男人狂女人下面高潮的视频| 国产亚洲av片在线观看秒播厂 | 久久国内精品自在自线图片| 国产成人一区二区在线| 精品午夜福利在线看| 一级av片app| 国产精品一区二区性色av| av视频在线观看入口| 欧美潮喷喷水| av国产免费在线观看| 99久国产av精品国产电影| 免费看光身美女| 长腿黑丝高跟| 在线免费观看不下载黄p国产| 一级毛片我不卡| 在线观看美女被高潮喷水网站| 亚洲成a人片在线一区二区| 少妇熟女aⅴ在线视频| 久久精品久久久久久噜噜老黄 | 麻豆久久精品国产亚洲av| 免费电影在线观看免费观看| 亚洲aⅴ乱码一区二区在线播放| 国产69精品久久久久777片| 国产伦在线观看视频一区| 日本五十路高清| 青春草亚洲视频在线观看| 亚洲成a人片在线一区二区| 免费av观看视频| 日本五十路高清| 日本免费a在线| 日韩一区二区三区影片| 91狼人影院| 哪个播放器可以免费观看大片| 欧美性猛交╳xxx乱大交人| 性欧美人与动物交配| 国产精品99久久久久久久久| 日日撸夜夜添| 国产v大片淫在线免费观看| 寂寞人妻少妇视频99o| 免费搜索国产男女视频| 欧美激情国产日韩精品一区| 哪里可以看免费的av片| 级片在线观看| 国产 一区精品| 在现免费观看毛片| 成人特级黄色片久久久久久久| 亚洲av一区综合| 美女被艹到高潮喷水动态| 丰满的人妻完整版| 直男gayav资源| 久久久午夜欧美精品| 久久九九热精品免费| 三级经典国产精品| 一区二区三区四区激情视频 | 免费观看人在逋| 欧美潮喷喷水| 日韩三级伦理在线观看| 午夜久久久久精精品| 精品久久久久久久久亚洲| 99热6这里只有精品| 永久网站在线| 久久国内精品自在自线图片| 亚洲人成网站在线播放欧美日韩| 免费搜索国产男女视频| 黄色一级大片看看| 最近中文字幕高清免费大全6| 午夜视频国产福利| 国产精品三级大全| 可以在线观看毛片的网站| 免费电影在线观看免费观看| 精品久久久久久久久久久久久| 91av网一区二区| 国产中年淑女户外野战色| 成年女人看的毛片在线观看| 性色avwww在线观看| 国产精品国产高清国产av| 国产一区二区三区在线臀色熟女| 日韩视频在线欧美| 亚洲av.av天堂| 中文字幕免费在线视频6| 久久99蜜桃精品久久| 亚洲欧美中文字幕日韩二区| 久久精品91蜜桃| 六月丁香七月| 亚洲av二区三区四区| 嫩草影院新地址| 日韩在线高清观看一区二区三区| 少妇的逼好多水| 国产又黄又爽又无遮挡在线| 色综合亚洲欧美另类图片| 成熟少妇高潮喷水视频| 日韩中字成人| 成人二区视频| 午夜福利在线在线| 好男人在线观看高清免费视频| 夫妻性生交免费视频一级片| 亚洲精品影视一区二区三区av| 免费看光身美女| 美女xxoo啪啪120秒动态图| 成人欧美大片| 国产一区二区亚洲精品在线观看| 免费一级毛片在线播放高清视频| 亚洲一级一片aⅴ在线观看| 色尼玛亚洲综合影院| 欧美bdsm另类| 欧美高清性xxxxhd video| 免费观看a级毛片全部| 久久99热6这里只有精品| 男人舔奶头视频| 午夜爱爱视频在线播放| 寂寞人妻少妇视频99o| 高清毛片免费观看视频网站| 欧美日韩一区二区视频在线观看视频在线 | 18禁裸乳无遮挡免费网站照片| 国产探花极品一区二区| 嫩草影院新地址| 熟妇人妻久久中文字幕3abv| 日本爱情动作片www.在线观看| 亚洲性久久影院| 嫩草影院精品99| 久久久欧美国产精品| 亚洲国产欧美人成| 免费一级毛片在线播放高清视频| 国国产精品蜜臀av免费| 精品久久久久久久人妻蜜臀av| 九色成人免费人妻av| 亚洲欧美日韩高清在线视频| 97超视频在线观看视频| 国产av一区在线观看免费| 国产精品av视频在线免费观看| 国产91av在线免费观看| 久久精品久久久久久噜噜老黄 | 天天一区二区日本电影三级| 欧美在线一区亚洲| 婷婷精品国产亚洲av| 天天躁日日操中文字幕| 国产精品1区2区在线观看.| 69av精品久久久久久| 国产成人freesex在线| 搡老妇女老女人老熟妇| 国产精品国产高清国产av| 欧美激情在线99| www.av在线官网国产| 一区二区三区四区激情视频 | 亚洲最大成人手机在线| 我的老师免费观看完整版| 日本欧美国产在线视频| 久久久久久九九精品二区国产| 麻豆成人午夜福利视频| 边亲边吃奶的免费视频| 51国产日韩欧美| 亚洲av中文av极速乱| 国产中年淑女户外野战色| 日韩精品青青久久久久久| 在线免费观看不下载黄p国产| 男人舔奶头视频| 欧美色欧美亚洲另类二区| 亚洲精品久久久久久婷婷小说 | 欧美高清性xxxxhd video| 久久精品夜夜夜夜夜久久蜜豆| 99热这里只有精品一区| 国产私拍福利视频在线观看| 乱码一卡2卡4卡精品| 精品国产三级普通话版| 久久精品久久久久久久性| 我的老师免费观看完整版| 日韩欧美在线乱码| 黄片无遮挡物在线观看| 日韩高清综合在线| 可以在线观看毛片的网站| 在线天堂最新版资源| 91久久精品国产一区二区三区| 极品教师在线视频| 悠悠久久av| 亚洲最大成人手机在线| 国产爱豆传媒在线观看| 久久99热6这里只有精品| 精品人妻视频免费看| 特大巨黑吊av在线直播| 一区福利在线观看| 国产亚洲av片在线观看秒播厂 | 尾随美女入室| 亚洲欧美日韩卡通动漫| .国产精品久久| 久久久久九九精品影院| 久久久久久九九精品二区国产| 最近2019中文字幕mv第一页| 只有这里有精品99| 午夜精品国产一区二区电影 | av免费观看日本| 人妻久久中文字幕网| 18禁黄网站禁片免费观看直播| 中出人妻视频一区二区| 菩萨蛮人人尽说江南好唐韦庄 | 老司机福利观看| 黄色日韩在线| 免费av不卡在线播放| a级毛片免费高清观看在线播放| 亚洲欧美精品综合久久99| 成人毛片a级毛片在线播放| 国产亚洲精品久久久com| 亚洲中文字幕日韩| 免费看美女性在线毛片视频| 欧美极品一区二区三区四区| 黄色日韩在线| 激情 狠狠 欧美| 老司机福利观看| 女同久久另类99精品国产91| 久久久久久国产a免费观看| 99热这里只有精品一区| 国产高潮美女av| av在线天堂中文字幕| 国产高清三级在线| 欧美日韩在线观看h| 一本久久中文字幕| 日韩欧美在线乱码| 看非洲黑人一级黄片| 日韩 亚洲 欧美在线| 一个人看视频在线观看www免费| 国产亚洲av嫩草精品影院| 日韩国内少妇激情av| 丰满乱子伦码专区| 日韩视频在线欧美| 你懂的网址亚洲精品在线观看 | 亚洲第一区二区三区不卡| 99热只有精品国产| 麻豆av噜噜一区二区三区| 免费观看在线日韩| 午夜免费激情av| www.av在线官网国产| 久久久精品欧美日韩精品| 成熟少妇高潮喷水视频| 久久99精品国语久久久| 人妻少妇偷人精品九色| 99视频精品全部免费 在线| 女人十人毛片免费观看3o分钟| 中文字幕人妻熟人妻熟丝袜美| 三级国产精品欧美在线观看| 18禁在线播放成人免费| 国产精品野战在线观看| 3wmmmm亚洲av在线观看| 免费观看的影片在线观看| 国产黄片美女视频| 天天躁日日操中文字幕| 不卡视频在线观看欧美| 深爱激情五月婷婷| 久久亚洲国产成人精品v| 精品午夜福利在线看| 国产成人a∨麻豆精品| 男人狂女人下面高潮的视频| 久久欧美精品欧美久久欧美| 国产成人影院久久av| 国产探花极品一区二区| 99九九线精品视频在线观看视频| 午夜精品一区二区三区免费看| 欧美bdsm另类| 麻豆av噜噜一区二区三区| 午夜福利在线在线| 日本黄大片高清| 亚洲高清免费不卡视频| 赤兔流量卡办理| 麻豆国产97在线/欧美| 久久久a久久爽久久v久久| 黄片无遮挡物在线观看| 国产麻豆成人av免费视频| 欧美+日韩+精品| 国产免费一级a男人的天堂| 青青草视频在线视频观看| 日韩一区二区视频免费看| 五月玫瑰六月丁香| av黄色大香蕉| 18禁黄网站禁片免费观看直播| 秋霞在线观看毛片| 99久久九九国产精品国产免费| 黄色视频,在线免费观看| 国产 一区 欧美 日韩| 色噜噜av男人的天堂激情| 热99re8久久精品国产| 天堂中文最新版在线下载 | 3wmmmm亚洲av在线观看| 狂野欧美白嫩少妇大欣赏| 高清午夜精品一区二区三区 | 三级男女做爰猛烈吃奶摸视频| 最后的刺客免费高清国语| 免费av观看视频| 亚洲自偷自拍三级| 久久草成人影院| 精品人妻偷拍中文字幕| 嫩草影院入口| 欧美高清成人免费视频www| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 色噜噜av男人的天堂激情| 午夜激情欧美在线| 永久网站在线| 18禁在线无遮挡免费观看视频| 一个人看的www免费观看视频| 日韩av在线大香蕉| 国产精品国产高清国产av| 床上黄色一级片| 国产激情偷乱视频一区二区| 淫秽高清视频在线观看| 天堂网av新在线| 久久这里只有精品中国| 久久99热6这里只有精品| 国产精品一区二区性色av| 国产精品久久久久久av不卡| 久久久久免费精品人妻一区二区| 成人漫画全彩无遮挡| 国产又黄又爽又无遮挡在线| 哪个播放器可以免费观看大片| 国产69精品久久久久777片| 欧美精品一区二区大全| 最近的中文字幕免费完整| 国产一区二区三区av在线 | 久久人人精品亚洲av| 国产精品久久电影中文字幕| 高清日韩中文字幕在线| 日韩视频在线欧美| 99九九线精品视频在线观看视频| 婷婷亚洲欧美| 春色校园在线视频观看| 91午夜精品亚洲一区二区三区| 麻豆av噜噜一区二区三区| 久久精品久久久久久久性| 亚洲欧美日韩无卡精品| 99久国产av精品| 99久久中文字幕三级久久日本| 成人一区二区视频在线观看| 最近视频中文字幕2019在线8| 久久久久国产网址| 最后的刺客免费高清国语| av在线蜜桃| 久久99热6这里只有精品| 午夜a级毛片| 观看免费一级毛片| 国产一区二区激情短视频| 特级一级黄色大片| 欧美成人一区二区免费高清观看| 三级国产精品欧美在线观看| 国产精品嫩草影院av在线观看| ponron亚洲| 成人av在线播放网站| 久久欧美精品欧美久久欧美| 如何舔出高潮| 日韩欧美三级三区| 亚洲三级黄色毛片| 国产精品精品国产色婷婷| 白带黄色成豆腐渣| 亚洲最大成人手机在线| 1000部很黄的大片| 久久精品夜色国产| 午夜福利高清视频| 麻豆乱淫一区二区| 男人舔女人下体高潮全视频| 最近2019中文字幕mv第一页| 久久99蜜桃精品久久| 国产精品一二三区在线看| a级一级毛片免费在线观看| 哪里可以看免费的av片| 日韩成人伦理影院| 国产成人a∨麻豆精品| 一个人看的www免费观看视频| 春色校园在线视频观看| 在线观看av片永久免费下载| 99热全是精品| 日韩人妻高清精品专区| 国产视频内射| 日本色播在线视频| 内地一区二区视频在线| 国产男人的电影天堂91| 日韩欧美在线乱码| 日本免费a在线| 大香蕉久久网| 午夜福利视频1000在线观看| 国产高潮美女av| 久久久国产成人免费| 国产亚洲5aaaaa淫片| 成熟少妇高潮喷水视频| 免费av毛片视频| 乱码一卡2卡4卡精品| 亚洲精品色激情综合| 国产探花在线观看一区二区| 亚洲欧洲日产国产| 亚洲欧美成人精品一区二区| 免费观看人在逋| 免费无遮挡裸体视频| 最后的刺客免费高清国语| 国产极品精品免费视频能看的| 熟女人妻精品中文字幕| 亚洲乱码一区二区免费版|