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

    热99国产精品久久久久久7| 国产欧美日韩一区二区三区在线| 久热这里只有精品99| 天天影视国产精品| 成在线人永久免费视频| 国产区一区二久久| 91九色精品人成在线观看| 亚洲欧美精品综合一区二区三区| 人人妻,人人澡人人爽秒播| 欧美午夜高清在线| 99国产精品一区二区蜜桃av| 成人国语在线视频| 欧美乱码精品一区二区三区| 91成人精品电影| 久久久久久大精品| 中文字幕人妻熟女乱码| 一级毛片高清免费大全| 久久欧美精品欧美久久欧美| 叶爱在线成人免费视频播放| 精品一区二区三区四区五区乱码| 亚洲自拍偷在线| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲国产欧美网| 久久久久久人人人人人| 午夜福利在线免费观看网站| 亚洲成a人片在线一区二区| 真人做人爱边吃奶动态| 免费看a级黄色片| 日韩三级视频一区二区三区| 久久午夜亚洲精品久久| 国产精品野战在线观看 | 国产亚洲欧美98| 丝袜美足系列| 色在线成人网| 国产成人精品久久二区二区91| 首页视频小说图片口味搜索| 亚洲一区中文字幕在线| 日本欧美视频一区| 最近最新免费中文字幕在线| 多毛熟女@视频| 一本综合久久免费| 亚洲精品一区av在线观看| 午夜福利免费观看在线| 国产av又大| 9色porny在线观看| 日韩中文字幕欧美一区二区| 长腿黑丝高跟| 精品一区二区三区四区五区乱码| 日韩国内少妇激情av| 久久国产乱子伦精品免费另类| 成在线人永久免费视频| 成人三级做爰电影| 中文亚洲av片在线观看爽| 一进一出好大好爽视频| 男人操女人黄网站| 丝袜美足系列| 国产真人三级小视频在线观看| 51午夜福利影视在线观看| 国产成人欧美在线观看| 午夜激情av网站| 一区二区日韩欧美中文字幕| 成人手机av| 人人澡人人妻人| 激情视频va一区二区三区| 免费在线观看完整版高清| 日韩高清综合在线| 人妻丰满熟妇av一区二区三区| 999久久久精品免费观看国产| 免费在线观看完整版高清| 成年人黄色毛片网站| 91老司机精品| 老汉色∧v一级毛片| 亚洲第一青青草原| 美女 人体艺术 gogo| 色播在线永久视频| 男人舔女人下体高潮全视频| 午夜福利在线观看吧| 伦理电影免费视频| 老司机靠b影院| 欧美日本中文国产一区发布| e午夜精品久久久久久久| 国产伦一二天堂av在线观看| 日韩精品免费视频一区二区三区| 日韩大码丰满熟妇| 真人做人爱边吃奶动态| 精品久久久久久,| 亚洲七黄色美女视频| 熟女少妇亚洲综合色aaa.| 欧美日韩av久久| 国产免费av片在线观看野外av| 18美女黄网站色大片免费观看| 国产片内射在线| 一级a爱视频在线免费观看| 他把我摸到了高潮在线观看| 99香蕉大伊视频| 中文字幕av电影在线播放| www.自偷自拍.com| 国产亚洲精品久久久久5区| 欧美一级毛片孕妇| 欧美日韩精品网址| 天堂√8在线中文| 免费在线观看黄色视频的| 看黄色毛片网站| 国产亚洲精品第一综合不卡| 中亚洲国语对白在线视频| 级片在线观看| 99久久人妻综合| 精品国产乱码久久久久久男人| 精品久久蜜臀av无| 亚洲av五月六月丁香网| 成年人免费黄色播放视频| 国产国语露脸激情在线看| 在线观看舔阴道视频| 天天躁狠狠躁夜夜躁狠狠躁| 一区二区三区精品91| 免费观看精品视频网站| 国产在线观看jvid| 黄片大片在线免费观看| 91在线观看av| 丰满的人妻完整版| 日本wwww免费看| 国产精品久久久久成人av| 亚洲精品久久成人aⅴ小说| 国产黄色免费在线视频| 午夜91福利影院| 日韩成人在线观看一区二区三区| 69av精品久久久久久| 最新美女视频免费是黄的| 免费在线观看影片大全网站| 日韩欧美免费精品| av天堂久久9| 99久久久亚洲精品蜜臀av| 免费高清在线观看日韩| 久久久久精品国产欧美久久久| 侵犯人妻中文字幕一二三四区| 99热只有精品国产| 亚洲一区高清亚洲精品| 亚洲欧美激情在线| 色综合站精品国产| 国产精品永久免费网站| 99久久人妻综合| 亚洲av熟女| 级片在线观看| 国产精品亚洲一级av第二区| 午夜福利一区二区在线看| 高清黄色对白视频在线免费看| 很黄的视频免费| 黑人欧美特级aaaaaa片| 9热在线视频观看99| 亚洲成人免费av在线播放| 久久久久久久午夜电影 | 亚洲精品中文字幕在线视频| 在线观看一区二区三区| av有码第一页| 制服诱惑二区| 一进一出抽搐动态| 制服人妻中文乱码| 亚洲va日本ⅴa欧美va伊人久久| 精品人妻1区二区| 露出奶头的视频| 国产极品粉嫩免费观看在线| 亚洲在线自拍视频| 激情视频va一区二区三区| 日韩欧美三级三区| 51午夜福利影视在线观看| 操美女的视频在线观看| 侵犯人妻中文字幕一二三四区| 成年人免费黄色播放视频| 国产精品亚洲av一区麻豆| 一级a爱视频在线免费观看| 国产一区二区三区综合在线观看| 国产单亲对白刺激| 国产99白浆流出| 亚洲九九香蕉| 91老司机精品| 99re在线观看精品视频| 一进一出抽搐动态| 又紧又爽又黄一区二区| 淫秽高清视频在线观看| 午夜成年电影在线免费观看| 悠悠久久av| 欧美日韩乱码在线| 人妻丰满熟妇av一区二区三区| 久久人人精品亚洲av| 一个人观看的视频www高清免费观看 | 国产av一区在线观看免费| 不卡一级毛片| 啪啪无遮挡十八禁网站| 久久久久久久久中文| 国产又爽黄色视频| 免费看十八禁软件| 在线播放国产精品三级| 久久香蕉精品热| 美女高潮喷水抽搐中文字幕| 精品一区二区三区四区五区乱码| 久久精品91蜜桃| 国产av精品麻豆| 国产激情久久老熟女| 欧美亚洲日本最大视频资源| av视频免费观看在线观看| www.www免费av| 人人妻人人添人人爽欧美一区卜| 国产精品免费视频内射| 两个人免费观看高清视频| 久久亚洲真实| 精品日产1卡2卡| 九色亚洲精品在线播放| 99re在线观看精品视频| 国产亚洲欧美在线一区二区| 一级作爱视频免费观看| 免费观看精品视频网站| 成人av一区二区三区在线看| 在线天堂中文资源库| 国产av又大| 丁香六月欧美| 欧美+亚洲+日韩+国产| 91成年电影在线观看| 日本vs欧美在线观看视频| 午夜日韩欧美国产| 97碰自拍视频| 99国产综合亚洲精品| 亚洲免费av在线视频| 午夜福利影视在线免费观看| 国产三级在线视频| 中文字幕人妻丝袜制服| 欧美一级毛片孕妇| 9热在线视频观看99| 午夜免费观看网址| 一级作爱视频免费观看| 欧美日韩国产mv在线观看视频| 亚洲欧美日韩无卡精品| 免费搜索国产男女视频| 国产伦一二天堂av在线观看| 国产欧美日韩精品亚洲av| 欧美成人免费av一区二区三区| 国产成人欧美| 十八禁人妻一区二区| 欧美日韩国产mv在线观看视频| 亚洲,欧美精品.| 可以免费在线观看a视频的电影网站| 老汉色∧v一级毛片| 交换朋友夫妻互换小说| 无限看片的www在线观看| 亚洲av成人av| 制服人妻中文乱码| 久热这里只有精品99| 精品一区二区三区四区五区乱码| 成人av一区二区三区在线看| 91国产中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 日韩三级视频一区二区三区| 国产人伦9x9x在线观看| 精品一品国产午夜福利视频| 大码成人一级视频| 最好的美女福利视频网| 国产成人av教育| 一级黄色大片毛片| 啪啪无遮挡十八禁网站| 免费在线观看视频国产中文字幕亚洲| 嫩草影院精品99| 午夜免费鲁丝| 国产精品免费一区二区三区在线| 高清在线国产一区| 午夜a级毛片| 午夜老司机福利片| 精品电影一区二区在线| 国产精品偷伦视频观看了| www日本在线高清视频| 国产精品日韩av在线免费观看 | 亚洲欧美精品综合一区二区三区| 99国产精品一区二区蜜桃av| 不卡av一区二区三区| 琪琪午夜伦伦电影理论片6080| 久久午夜综合久久蜜桃| 午夜福利一区二区在线看| 欧美丝袜亚洲另类 | 免费观看人在逋| 老司机福利观看| 精品少妇一区二区三区视频日本电影| 天堂动漫精品| 国产精品秋霞免费鲁丝片| 免费在线观看视频国产中文字幕亚洲| 久久精品国产亚洲av香蕉五月| 国产精品偷伦视频观看了| 看片在线看免费视频| 欧美日本亚洲视频在线播放| 免费搜索国产男女视频| 欧美久久黑人一区二区| 欧美+亚洲+日韩+国产| 51午夜福利影视在线观看| 欧美大码av| 日韩精品青青久久久久久| 精品免费久久久久久久清纯| 黑人猛操日本美女一级片| 女人被躁到高潮嗷嗷叫费观| 91av网站免费观看| 桃色一区二区三区在线观看| 老司机在亚洲福利影院| www.999成人在线观看| 精品国产一区二区久久| 色综合婷婷激情| 欧美在线黄色| 午夜免费激情av| 黄色毛片三级朝国网站| 久久国产乱子伦精品免费另类| 午夜福利在线观看吧| 可以在线观看毛片的网站| 日韩成人在线观看一区二区三区| 国产精品久久视频播放| 精品国内亚洲2022精品成人| 91成人精品电影| 久久天躁狠狠躁夜夜2o2o| 国产99久久九九免费精品| 亚洲伊人色综图| 夜夜爽天天搞| 看片在线看免费视频| 午夜a级毛片| 别揉我奶头~嗯~啊~动态视频| 午夜福利影视在线免费观看| 欧美激情 高清一区二区三区| 极品人妻少妇av视频| 免费在线观看影片大全网站| 亚洲一卡2卡3卡4卡5卡精品中文| 色哟哟哟哟哟哟| 1024视频免费在线观看| 侵犯人妻中文字幕一二三四区| 美女高潮到喷水免费观看| 亚洲第一欧美日韩一区二区三区| 激情视频va一区二区三区| 美女高潮到喷水免费观看| 91麻豆av在线| 黄色 视频免费看| 国产成人精品无人区| av天堂久久9| 不卡一级毛片| 欧美激情高清一区二区三区| 搡老乐熟女国产| av中文乱码字幕在线| 在线观看午夜福利视频| 精品卡一卡二卡四卡免费| 婷婷六月久久综合丁香| 日韩有码中文字幕| 久久久久国产精品人妻aⅴ院| 亚洲成av片中文字幕在线观看| 最新在线观看一区二区三区| 在线观看66精品国产| 五月开心婷婷网| 午夜成年电影在线免费观看| 亚洲全国av大片| 中亚洲国语对白在线视频| 女性被躁到高潮视频| 免费观看人在逋| av电影中文网址| 国产精品爽爽va在线观看网站 | 久9热在线精品视频| 亚洲精品在线美女| 午夜精品久久久久久毛片777| 欧美亚洲日本最大视频资源| a在线观看视频网站| 97人妻天天添夜夜摸| 欧美在线一区亚洲| 少妇 在线观看| 18美女黄网站色大片免费观看| 黄色成人免费大全| 嫁个100分男人电影在线观看| 女同久久另类99精品国产91| 久久精品亚洲av国产电影网| 国产成人欧美| 日韩大码丰满熟妇| 国产精品乱码一区二三区的特点 | 精品久久久久久久毛片微露脸| 黄色丝袜av网址大全| 91成人精品电影| 一级毛片高清免费大全| 欧美午夜高清在线| 中文字幕人妻丝袜制服| 自拍欧美九色日韩亚洲蝌蚪91| 悠悠久久av| 国产精品久久久久久人妻精品电影| 丁香六月欧美| 国产91精品成人一区二区三区| 国产成人欧美在线观看| 欧美日韩亚洲高清精品| 麻豆av在线久日| 在线观看舔阴道视频| 宅男免费午夜| 亚洲欧美精品综合久久99| 激情视频va一区二区三区| 18禁美女被吸乳视频| av国产精品久久久久影院| 日本 av在线| 性少妇av在线| 999精品在线视频| 亚洲精华国产精华精| 亚洲中文日韩欧美视频| a级片在线免费高清观看视频| 欧美黄色淫秽网站| 老汉色∧v一级毛片| 亚洲国产欧美一区二区综合| 女人爽到高潮嗷嗷叫在线视频| 免费高清视频大片| 多毛熟女@视频| 色在线成人网| 精品熟女少妇八av免费久了| 在线观看66精品国产| 性色av乱码一区二区三区2| 一级毛片精品| a级片在线免费高清观看视频| 婷婷丁香在线五月| 黑人巨大精品欧美一区二区蜜桃| 很黄的视频免费| 无遮挡黄片免费观看| 国产单亲对白刺激| 国产av一区二区精品久久| 91大片在线观看| www.999成人在线观看| 亚洲精品av麻豆狂野| 午夜福利在线免费观看网站| 免费日韩欧美在线观看| 在线观看午夜福利视频| 久久久久久免费高清国产稀缺| 真人一进一出gif抽搐免费| 老司机靠b影院| 久久久精品国产亚洲av高清涩受| 亚洲国产精品sss在线观看 | 在线观看66精品国产| 久久精品亚洲av国产电影网| 久久人人97超碰香蕉20202| 性少妇av在线| 在线观看免费视频网站a站| av天堂久久9| 老汉色av国产亚洲站长工具| 久久久久久久久中文| 日本免费a在线| 欧美不卡视频在线免费观看 | 久久久久久久久中文| 精品国产乱码久久久久久男人| 色老头精品视频在线观看| 免费看十八禁软件| 久久久久久人人人人人| 女人被狂操c到高潮| av在线播放免费不卡| 成人影院久久| av有码第一页| 日韩国内少妇激情av| a级毛片在线看网站| 精品一区二区三区视频在线观看免费 | 国产精品一区二区在线不卡| 91精品三级在线观看| 午夜日韩欧美国产| 69av精品久久久久久| 午夜两性在线视频| 国产成人av教育| 美女大奶头视频| 男女高潮啪啪啪动态图| 啦啦啦 在线观看视频| 精品久久久精品久久久| av超薄肉色丝袜交足视频| 久久久久久免费高清国产稀缺| 国产成人av教育| 国产成人系列免费观看| 精品一区二区三区视频在线观看免费 | 人妻久久中文字幕网| 午夜免费观看网址| 一级毛片精品| 最近最新中文字幕大全免费视频| 天堂俺去俺来也www色官网| 一进一出抽搐动态| 成人三级黄色视频| 天堂影院成人在线观看| 国产亚洲欧美在线一区二区| 757午夜福利合集在线观看| 日韩大码丰满熟妇| 国产深夜福利视频在线观看| 欧美另类亚洲清纯唯美| 一二三四社区在线视频社区8| 午夜91福利影院| 久久久久久久久久久久大奶| 亚洲精品中文字幕在线视频| 首页视频小说图片口味搜索| 亚洲精品一二三| 男男h啪啪无遮挡| 51午夜福利影视在线观看| 麻豆成人av在线观看| 在线观看免费午夜福利视频| 国产真人三级小视频在线观看| 国产日韩一区二区三区精品不卡| 国产一区二区三区在线臀色熟女 | 欧美色视频一区免费| 99久久人妻综合| 午夜久久久在线观看| 国产精品一区二区三区四区久久 | 国产成人欧美在线观看| a级片在线免费高清观看视频| 男男h啪啪无遮挡| 69av精品久久久久久| 不卡一级毛片| 少妇的丰满在线观看| 国产蜜桃级精品一区二区三区| 国产精品电影一区二区三区| 90打野战视频偷拍视频| 一边摸一边抽搐一进一出视频| tocl精华| 亚洲欧美日韩无卡精品| 久久久久久久久免费视频了| 久久中文字幕一级| 亚洲三区欧美一区| e午夜精品久久久久久久| 少妇被粗大的猛进出69影院| 丝袜美足系列| 1024视频免费在线观看| 黑人巨大精品欧美一区二区蜜桃| 久久伊人香网站| 国产午夜精品久久久久久| 超碰成人久久| 后天国语完整版免费观看| 久久欧美精品欧美久久欧美| 性少妇av在线| 波多野结衣av一区二区av| 亚洲精品国产一区二区精华液| 欧美午夜高清在线| 怎么达到女性高潮| 亚洲精品国产精品久久久不卡| 精品国内亚洲2022精品成人| 十八禁网站免费在线| 欧美最黄视频在线播放免费 | 国产国语露脸激情在线看| ponron亚洲| 国产成人av教育| 亚洲精品在线观看二区| 欧美黄色淫秽网站| 啦啦啦免费观看视频1| 国产亚洲欧美精品永久| 国产高清激情床上av| 伦理电影免费视频| 亚洲男人天堂网一区| 两个人看的免费小视频| 久久精品影院6| 精品高清国产在线一区| 亚洲性夜色夜夜综合| 两性午夜刺激爽爽歪歪视频在线观看 | 最近最新中文字幕大全免费视频| 水蜜桃什么品种好| 男人的好看免费观看在线视频 | 19禁男女啪啪无遮挡网站| 久久国产亚洲av麻豆专区| 天天添夜夜摸| 亚洲va日本ⅴa欧美va伊人久久| 69av精品久久久久久| 老熟妇仑乱视频hdxx| 亚洲国产中文字幕在线视频| 琪琪午夜伦伦电影理论片6080| 精品福利观看| 精品国产一区二区久久| 1024香蕉在线观看| 无限看片的www在线观看| 如日韩欧美国产精品一区二区三区| 亚洲第一av免费看| 国产激情欧美一区二区| 免费搜索国产男女视频| 99久久人妻综合| 久久人妻av系列| 90打野战视频偷拍视频| 亚洲三区欧美一区| 久久久久国产一级毛片高清牌| 在线天堂中文资源库| 水蜜桃什么品种好| 免费不卡黄色视频| 中文字幕人妻丝袜一区二区| a级毛片黄视频| 国产精品秋霞免费鲁丝片| 国产91精品成人一区二区三区| 亚洲专区字幕在线| 一级作爱视频免费观看| 80岁老熟妇乱子伦牲交| 韩国av一区二区三区四区| 青草久久国产| 久久狼人影院| 国产亚洲精品综合一区在线观看 | 亚洲少妇的诱惑av| 亚洲情色 制服丝袜| а√天堂www在线а√下载| 精品电影一区二区在线| 欧美大码av| 日本a在线网址| 日韩免费av在线播放| 久久久国产成人精品二区 | xxx96com| 国产一区二区三区在线臀色熟女 | 国产成人精品无人区| 亚洲午夜精品一区,二区,三区| 1024视频免费在线观看| 999精品在线视频| 精品无人区乱码1区二区| 国产成人系列免费观看| 黄色丝袜av网址大全| 日韩三级视频一区二区三区| 美国免费a级毛片| 亚洲一区二区三区欧美精品| 黄色 视频免费看| 欧美丝袜亚洲另类 | 啦啦啦在线免费观看视频4| 日本免费a在线| 亚洲欧美激情在线| 午夜老司机福利片| 国产不卡一卡二| 人人妻人人添人人爽欧美一区卜| 黑丝袜美女国产一区| 国产精品99久久99久久久不卡| 久久九九热精品免费| 新久久久久国产一级毛片| 国产精品永久免费网站| 黄色视频不卡| 国产精品久久久久成人av| 高清毛片免费观看视频网站 |