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

    Adaptive Sensor-Fault Tolerant Control of Unmanned Underwater Vehicles With Input Saturation

    2024-04-15 09:36:54XueraoWangQinglingWangYanxuSuYunchengOuyangandChangyinSun
    IEEE/CAA Journal of Automatica Sinica 2024年4期

    Xuerao Wang , Qingling Wang , Yanxu Su , Yuncheng Ouyang , and Changyin Sun ,,

    Abstract—This paper investigates the tracking control problem for unmanned underwater vehicles (UUVs) systems with sensor faults, input saturation, and external disturbance caused by waves and ocean currents.An active sensor fault-tolerant control scheme is proposed.First, the developed method only requires the inertia matrix of the UUV, without other dynamic information,and can handle both additive and multiplicative sensor faults.Subsequently, an adaptive fault-tolerant controller is designed to achieve asymptotic tracking control of the UUV by employing robust integral of the sign of error feedback method.It is shown that the effect of sensor faults is online estimated and compensated by an adaptive estimator.With the proposed controller, the tracking error and estimation error can asymptotically converge to zero.Finally, simulation results are performed to demonstrate the effectiveness of the proposed method.

    I.INTRODUCTION

    OVER the past few years, considerable progress has been made in the area of underwater robotics, which has received extensive attention.Among them, the control of unmanned underwater vehicles (UUVs) is a popular research field with a significant number of positive results [1]–[6].Benefiting from these results, UUVs are widely applied in the ocean environment, such as ocean forecasting, underwater inspection, and pipeline tracking.However, although the development of the oceanic industries has grown rapidly, the research of UUVs is limited by the complexity and variability of the marine environment [7].One of the core limitations is that UUVs are prone to failure in a complex environment.Therefore, the ability to handle various faults is essential to the achievement of vehicle’s mission.

    The most common failures in UUV systems are of two kinds: actuator faults and sensor faults.To solve these problems, the fault-tolerant control (FTC) method has been well investigated to handle the potential faults and maintain system stability [8]–[13].In [14], a unified actuator fault model with the consideration of four different types of faults was first proposed for an unmanned marine vehicle.Based on the established model, the authors gave a quantized sliding mode FTC design procedure, in which the quantization parameters can be adjusted by the lower bound of the fault information.A dynamic output feedback FTC scheme was proposed in [15].The scheme includes a dynamic compensator and an adaptive learning observer to deal with the compound uncertainties so that the tracking errors can be bounded by a small value.However, while there has been significant study on the actuator FTC, less exploration has been carried out on the sensor FTC because it is more complicated and changeable [16].In reality, the closed-loop control of the UUV depends heavily on the measurement of position and velocity information by the shipboard sensors.Hence, even a small measured deviation may seriously deteriorate the control performance.Therefore, the research of sensor FTC schemes for UUV is indispensable and meaningful.

    The most serious problem caused by sensor faults is the inability to obtain accurate state information, which may result in fatal errors in feedback control.To better compensate for the biased measurement information, fault estimation is essential, where the estimate information is utilized to compensate for the effect of unknown possible faults on the control performance.To tackle this problem, the design and analysis of fault estimation and compensation have attracted lots of research attention [17]–[19].A unified FTC architecture including sensor fault detection and isolation was proposed in[20], which consists of a bank of nonlinear adaptive estimators.Specifically, the estimator in [20] can detect sensor faults and isolate the biased dynamics after fault detection.In [21],the system states, sensor faults, and disturbances were grouped into a new state which can be estimated by a novel robust observer simultaneously.Then the developed fault system can be stabilized with the estimation results.However,there are some issues with the aforementioned studies that require further investigation: 1) These developed FTC schemes do not provide sufficiently satisfactory behavior, e.g.,the asymptotic convergence of tracking errors; and 2) Most of the sensor FTC control schemes for the UUV system consider only one or two kinds of sensor faults, which reduces their applicability in the complicated marine environment.

    A significant problem for all unmanned systems is how to ensure the tracking performance under autonomous operation[22]–[24].Motivated to maintain high tracking performance in the presence of failures, varieties of control techniques are investigated.In [25], an observer-based FTC for an underwater vehicle was proposed with a finite-time adaptive sliding mode observer, which can ensure the estimation errors converge to zero in finite time.A composite FTC scheme was designed in [26] with three subsystems for fault detection, isolation, and control, where the tracking errors of the UUV were guaranteed to meet the preset performance.However, the existence of possible faults makes it difficult to achieve asymptotic tracking results, a more effective high-performance control scheme is required.Recently, an important result, that is,the robust integral of the sign of error (RISE) feedback control technique was first developed in [27].The basic principle of this method is to utilize a continuous high-gain controller with an integral robust term, to guarantee the controlled system asymptotically stable, which provides a new solution to the aforementioned problem.Whereas, RISE control also brings the problem of excessively high control gain.For UUVs, since the energy of the actuators is limited, the highgain controller cannot be realized, and the input saturation is taken into account in this paper.To implement the trajectory tracking control with input saturation for UUVs, several results have been obtained in [28]–[30].But few studies have considered the simultaneous existence of both sensor faults and input saturation.Therefore, it is still challenging to develop a sensor FTC scheme to ensure asymptotic stability for UUVs subject to input saturation and external disturbance.

    In this paper, an active sensor FTC scheme is proposed for UUVs subject to sensor faults, input saturation, and external disturbances to achieve asymptotic tracking performance.The main contributions of this paper are summarized as follows:

    1) The proposed FTC in this paper can deal with both additive and multiplicative types of sensor faults, which is different from existing results [26], [31], [32] that only consider one or two types of sensor faults.Then, an online estimator is given to estimate the sensor faults without traditional fault detection and isolation (FDI) module.

    2) In contrast to [7], [28], [33], the nonlinear dynamic knowledge of the UUV or the sensor fault information is not required in this paper.The increased uncertainty makes it more challenging to deal with sensor faults, but enhances the availability of the developed FTC scheme in practical implementation.

    3) The asymptotic tracking performance is guaranteed with the proposed RISE-based active FTC scheme, which is different from the existing UUV control methods [7], [25], [34].To the best of our knowledge, this is the first asymptotic tracking result for UUVs that are capable of tolerating various sensor faults.

    The rest of this paper is organized as follows: Section II introduces the problem statement and preliminaries.The main results, including the design of error systems and the FTC controller, are given in Section III.The stability analysis of the closed-loop system is detailed in Section IV.Section V presents the simulation examples on a UUV platform.Finally,Section VI includes conclusions.

    II.PROBLEM FORMULATION AND PRELIMINARIES

    A. System Description

    Referring to [35], the nonlinear UUV model moving in the horizontal planes subject to ocean current disturbances and input saturation can be expressed in three DOFs as

    where η =[x,y,ψ]T∈R3describes the pose of the vehicle represented in the inertial reference frame, and ν=[u,v,r]T∈R3denotes the surge, sway, and yaw velocities represented in the body reference frame (see Fig.1).The rotation matrix between two frames is

    Fig.1.UUV model with two reference frames.

    In view of the physical limitation of the engines of UUVs,the command control force may exceed the action range of the actuators.To decrease the influence induced by this constraint, the input saturation of the actuators is considered in the controller design.Therefore, the control force and torque τsproduced by actuators suffers from asymmetric saturation nonlinearities as

    wherei=1, 2, 3, τ =[τ1,τ2,τ3]∈R3is the designed input signal to the saturated actuator.τbiis the known constant denoting constraints on thei-th element of τs.

    Using the conversion relation betweenηandν, we have

    Then, the motion equation of the UUV in inertial reference frame is rewritten as [26]

    On account of the sensor faults, some state variables cannot be utilized directly in the control design.Hence, we define the measured output by the sensor as(Rν)T]T.The subsequent section is based on the fact that the pose and velocity information of the UUV is obtainable using sensors (e.g., long baseline, magnetic compass, gyros), that is,ηandνare measurable [36].In practice, it can be known that the accurate value of the pose informationηis easily obtained via precision measuring instruments.However, due to the limitations of underwater measurement technology, the measured velocity valueνoften contains faulty information.Therefore,in this paper, we mainly investigate the case when the faults occurred in the velocity sensors, to guarantee the reliability and safety of UUVs.

    The control objective of this paper is to develop an asymptotic tracking control scheme only based on measurable output, so that the UUV model in (1) can track the desired trajectory in the presence of input saturation, external disturbance,and different types of sensor faults.The active sensor FTC structure for the UUV is shown in Fig.2.As illustrated in the diagram, only the imprecise measurements of the system output are used for feedback control.

    Assumption 1([37],[38]): The external disturbance τdis bounded by τd≤d, wheredis a positive constant.Moreover,the first two time derivatives of τdare bounded.

    B. Sensor Faults

    To distinguish the fault status and fault-free status, the output of velocity sensors with faults is defined as νf, and the sensors are supposed to be faulty at timet≥Tf.The additive faults and multiplicative faults are both considered in this paper, and the definitions of which are given as follows.

    Fig.2.Adaptive sensor fault tolerant control scheme.

    Definition 1(Sensor Additive Fault): The sensor for measuring variable δ(t)∈R is identified to have an additive fault at timeTfif the measured output

    wherebf(t) denotes the accuracy coefficient parameter.It is noted thatbf(t) and its time derivative are bounded.

    Definition 2(Sensor Multiplicative Fault): The sensor for measuring variable ι(t)∈R is identified to have a multiplicative fault at timeTfif the measured output

    where σf(t)∈(0,1] represents the effective factor.

    According to Definitions 1 and 2, the possible sensor faults can be modeled as

    where σ(t)=diag(σ1(t),σ2(t),σ3(t)) denotes the effectiveness matrix with σi(i=1,2,3)∈(0,1], andb=[b1,b2,b3]Tis the accuracy coefficient vector.On the basis of Definition 1, it can be assumed that the accuracy coefficient vectorb(t) and its time derivative are bounded by 0 ≤//b(t)//≤b1, and 0 ≤//b˙(t)//≤b2, respectively, whereb1andb2are unknown constants [21], [39].

    Remark 1: Definition 2 implies the following categories, in which three typical types of additive failures are included.

    1)Bias:

    2)Drift:

    3)Loss of Accuracy:

    These cases correspond to the classical type of sensor faults given in [40], [41].

    According to (8), the measured output subject to faults is redefined as

    wherefs=R(σ(t)-I3)ν+Rb(t) is an auxiliary term containing the overall information of sensor faults.According to the practical operation of the UUV and the limitation of the actuators, the acceleration of the UUV will be limited to a certain range.Adding that into the consideration with Definitions 1 and 2, it can be obtained thatf˙sis bounded by //f˙s//≤f, wherefis a known constant.

    III.CONTROLLER DESIGN

    In this section, an adaptive sensor-fault tolerant controller based on the RISE feedback method is proposed to guarantee the UUV can track the desired trajectory asymptotically.

    A. Error Systems

    To facilitate the subsequent design, we multiply both sides of (5) by an auxiliary multiplierR-T.Then, system (5) is equivalent to

    According to (12) and (13), we can get the dynamics of output variables with sensor faults as

    Letηd(t)∈R3be the reference trajectory such that

    Definey1d=ηdas the desired trajectory ofyf1.Then, the tracking error is defined as

    Now, we further design a virtual control signaly2d∈R3as

    Substituting (14) into (18), we have

    Next, the second filtering error variable is defined as

    whereβ ∈R3×3is the filter parameter matrix to be designed.Substituting (18) and its time derivative into (20) yields

    It is worth noting thatsis unavailable since the definition in(20) is based on η¨.Thus, the variablesis used only to analyze and prove the stability results.

    Now, after multiplying both sides of (21) byM(η), the open-loop tracking error system is developed as

    B. Controller Design

    Before designing the controller directly, an auxiliary dynamic system is proposed to eliminate the effect of input saturation.Defining an auxiliary input signal=R-Tτ, the difference between the designed input signal and the actual control input is denoted as

    An input saturation auxiliary dynamic system is designed as

    Invoking the developed open-loop system (22) and considering the input saturation effect, the FTC controller is developed as

    whereks∈R3and γ ∈R are positive control gain parameters,μ′=μ0sgn(s(ξ)) is a small adjustment term.Furthermore, differentiating (25) with respect to time, one has

    However, the input control signalincludes an auxiliary errorz, which consists of the sensor fault estimation value.To handle this problem, an adaptive estimator is presented to obtain the estimation value as

    where ?fis a positive constant estimator gain,?f=yf2-y˙f1-f?s-zis an auxiliary term.

    Motivated for the stability proof, we take the derivative of(22)

    Substituting the controller (26) into (28), we can obtain the closed-loop system as follows:

    ?N?D

    where and are given by

    and

    Assumption 2:The first two time derivatives of the saturated difference signal ?τare bounded.

    SinceN?(t) is continuously differentiable, taking a similar approach as in [27] and [38], the mean value theorem is employed to suppose thatN?(t) is bounded by

    φ(t)∈R9

    where the auxiliary error variable is designed as

    and the bounded function ρ(·):R →R is a positive globally invertible nondecreasing function.Due to the property ofηdgiven in (15), and according to Assumptions 1 and 2, it can be obtained that

    Remark 3: Many practical systems such as underwater vehicles always operate in some specified regions where they are controllable under the input saturation [47].Furthermore, the energy of UUV system is limited, which is inadequate to allow the overlarge change of control torque and system state.Consequently, it can be obtained from (25) and (26) that the control input signal and its derivative are bounded, and Assumption 2 is reasonable.

    Remark 4: In the control design, an adaptive estimator is developed to estimate the unknown term about sensor faults.It is highlighted that the derivative of measured statecan not be obtained directly.To solve this problem, the substitution algorithm given in [41] is utilized to calculate the current state.Then, the proposed adaptive estimator (27) can be integrated as

    IV.STABILITY ANALYSIS

    In this section, the stability of the resulting closed-loop system is rigorously proved.Before giving the main results, a useful lemma is shown and proved as follows, which has an important role in the proof of main results.

    Lemma 1:Define an auxiliary functionL(t)∈R as

    With (35), the integral ofL(t) satisfies the following inequality:

    whereLis a positive constant given by

    with the following condition holding:

    Proof: Substituting (21) into (35) and integrating in time,we have

    Considering (38), the upper bound of the right-hand of (39)is obtained as

    Thus, (36) holds.■

    Remark 5: The design ofL(t) in Lemma 1 is different from the traditional design of the RISE feedback method in [11],[48].The error systems are redesigned on account of the consideration of sensor faults, so thatL(t) and its property are adjusted accordingly.

    The main results are summarized as follows.

    Theorem 1: Consider the nonlinear and uncertain dynamics of the UUV illustrated by (1) with input saturation and sensor faults (8).Under Assumptions 1 and 2, the controller given in(24) and (25) with estimator (27) ensures the asymptotic tracking performance of UUV with all system states bounded,and the tracking error satisfies

    provided the control gainksis chosen sufficiently large, and the remaining control parameters are chosen with the following conditions:

    Proof:Define an auxiliary functionP(t) as

    It can be obtained from (36) in Lemma 1 thatP(t) is positive.Then, a Lyapunov candidate function is designed as

    where ζ (t) is defined as

    wherem1is a positive constant andm2(·):R →R denotes a positive nondecreasing function.Further, we can obtain the bounds ofVas

    Substituting (24), (27) and the time derivative of (43) into(48) yields

    Substituting (35) into (49) yields

    Based on (32), an auxiliary inequality can be derived as

    Then,V˙ can be further bounded by

    V.SIMULATION

    In this section, the proposed FTC scheme is applied to the trajectory tracking problem of UUV subject to ocean currents and external disturbances.The model considered in this simulation is omni directional intelligent navigator (ODIN) weighing 125 kg [49], [50].The actuating system consists of four horizontal thrusters and four vertical thrusters, and each thruster is a brushless DC motor that can provide a maximum thrust of approximately 27 N.Further details can be found in[51].It is noted that the horizontal movement is based on four horizontal thrusters [25].

    For the ODIN vehicle, the system matricesM,C(ν),D(ν) in(1) are given as

    Through the description of the ODIN structure in [51], there are four identical thrusters in the horizontal plane.The control configuration matrix is given by

    with θ=1/4π andRz=0.508.The reference trajectories are chosen as follows:

    To illustrate the effectiveness of the designed control algorithm, three tests with different fault types are conducted.The fault situations are specified as

    Case 1:The first test is supposed to perform a base result for comparison.Therefore, the UUV system is under fault-free condition as

    Case 2: In this test, we suppose that the velocity sensor bias fault occurs att=50 s.The bias fault is represented by

    Case 3: In the third test, we suppose that the velocity sensor bias fault and the multiplicative fault (loss of effectiveness of 25%) occur simultaneously att=80 s.The composite fault is represented by

    To provide a comparison, and verify the effectiveness of the proposed RISE-based sensor FTC method, we compare it with the other comparable integral terminal sliding mode (ITSM)control scheme.It is worth mentioning that the ITSM control scheme can guarantee the finite-time convergence of the error dynamics, which is commonly used in the control of underwater vehicles [52]–[55].The design of the ITSM-based controller refers to the results in [55], and the input saturation constraints are removed.In this case, we have

    wheresηis the kinematic ITSM surface defined as

    and the parameters are given as Γs=diag(1,1,1), αs=1,βs=diag(30,30,55),m= 3,n= 7.

    The simulation results are shown respectively in Figs.3-7.From Fig.3, it can be observed that the UUV can track the reference trajectory under the proposed controller, and the xyplane trace is a complete “8” type.The tracking performance in the fault-free case is shown in Fig.4, from which we can get that the tracking errors asymptotically converge to zero.The control signals given by the controller are depicted in Fig.4(c), and the control forces assigned to each thruster are shown in Fig.4(d).It can be obtained that the control forces and torque are bounded by the given restrained conditions.In the second case, the bias faults are considered, and the occurrence time is set to 50 s to verify the capability of the proposed scheme to cope with sudden faults.The tracking performance of the UUV under the bias sensor faults, input saturation, and external disturbance is illustrated in Fig.5.It can be observed from Figs.5(a) and 5(b) that even if the fault occurs,the system can still be restored to a stable state within 5 s.Fig.5(e) gives the measured and actual velocities to show the effect of sensor bias faults on the data measurement.Fig.5(f)depicts that the lumped faults can be estimated by the designed estimator.The third test verifies the ability of the proposed scheme to deal with composite faults, especially if there is already a bias fault.The tracking performance of the UUV under the composite sensor faults, input saturation, and external disturbance is illustrated in Fig.6.It can be observed from Figs.6(a)-6(d) that although both bias faults and multiplicative faults exist at the same time after 80 s, the UUV can still guarantee asymptotic tracking performance with the given input constraints satisfied.Fig.6(e) gives the measured and actual velocities to show the effect of sensor composite faults,Fig.6(f) presents the estimation results for the composite faults.

    The results in Fig.7 give the comparisons between the proposed controller and ITSM-based controller.First, it can be seen from Fig.7(a) both controllers can well track the desired trajectory in a short period, but the proposed controller provides smoother tracking curves.Moreover, Fig.7(b) shows that the tracking performance of the proposed RISE-based FTC scheme outperforms the ITSM-based controller in terms of overshoot, settling time, and steady-state error.Furthermore, the proposed RISE-based FTC scheme shows superior fault tolerance compared to the ITSM-based controller when a system failure occurs.From the above analysis, we can conclude that the designed controller has a good fault compensation effect and asymptotic tracking performance, which illustrates the effectiveness of the proposed FTC scheme.

    Fig.5.Case 2: Control performance under sensor bias faults.

    Remark 6: The tuning rules for the control parameters are built based on (42).First, the conditions given in (42) present a fundamental range of the parameter selection such as α and κ.Some parameter tuning skills are further used to select the exact value by analyzing the effect of each parameter.For example, α,ks, β are responsible for regulating the scale of filtering error.The parameters α andksjointly determine the proportional and differential error coefficients, which are recommended to be chosen as large as possible within the appropriate range.The parameter β affects the integration coefficient, which should be selected depending on the residuals.Furthermore, for coupling parameters such as β, ?f, the parameter that is easier to adjust can be selected preferentially,and the range of the remaining one will be determined easily.

    VI.CONCLUSION

    An adaptive sensor FTC scheme has been investigated for a UUV subject to sensor bias, drift, loss of accuracy, and loss of effectiveness with input saturation.To deal with the problem of input saturation, an auxiliary dynamic system is designed to eliminate the effect of constraints.The proposed sensor faulttolerant controller can compensate for the effects of sensor faults via an adaptive fault estimator.Moreover, by employing the RISE feedback method, the proposed scheme can achieve asymptotic tracking performance of the faulty UUV system.Simulation results are reported to validate the proposed control scheme has fault-tolerance capability and excellent tracking performance.

    Fig.6.Case 3: Control performance under sensor composite faults.

    Fig.7.Comparison results.Control performance under the action of RISE-based and ITSM-based controllers.

    a级一级毛片免费在线观看| 亚洲怡红院男人天堂| 一级毛片久久久久久久久女| 少妇人妻精品综合一区二区| 老女人水多毛片| 日韩国内少妇激情av| 日韩av在线免费看完整版不卡| 美女黄网站色视频| 大片免费播放器 马上看| 欧美丝袜亚洲另类| 80岁老熟妇乱子伦牲交| 我的老师免费观看完整版| 26uuu在线亚洲综合色| 日韩强制内射视频| 青春草国产在线视频| 免费人成在线观看视频色| 国产精品.久久久| 91aial.com中文字幕在线观看| 色综合站精品国产| 午夜激情福利司机影院| 在线播放无遮挡| 1000部很黄的大片| 亚洲内射少妇av| 国内精品宾馆在线| av专区在线播放| av在线蜜桃| 亚洲精品色激情综合| 久久午夜福利片| 狠狠精品人妻久久久久久综合| 亚洲无线观看免费| 日韩欧美精品v在线| 欧美性猛交╳xxx乱大交人| 亚洲成人久久爱视频| 午夜精品一区二区三区免费看| 九九在线视频观看精品| 一级片'在线观看视频| 春色校园在线视频观看| 亚洲在线自拍视频| ponron亚洲| 九色成人免费人妻av| 国内精品宾馆在线| 秋霞伦理黄片| 天天躁夜夜躁狠狠久久av| 一级片'在线观看视频| 久久久色成人| 国产在线男女| 日韩欧美精品v在线| 国产综合精华液| 建设人人有责人人尽责人人享有的 | 成人午夜精彩视频在线观看| 人妻制服诱惑在线中文字幕| 人妻一区二区av| 男女视频在线观看网站免费| 纵有疾风起免费观看全集完整版 | 久久久精品94久久精品| 男女边摸边吃奶| 日韩精品有码人妻一区| 中国国产av一级| 亚洲国产精品成人综合色| 91在线精品国自产拍蜜月| 嫩草影院新地址| 久久韩国三级中文字幕| 国产人妻一区二区三区在| 欧美丝袜亚洲另类| videos熟女内射| 中文资源天堂在线| 免费观看在线日韩| 日本午夜av视频| 爱豆传媒免费全集在线观看| 免费看美女性在线毛片视频| 亚洲综合精品二区| 国产中年淑女户外野战色| 国产乱人偷精品视频| 免费av观看视频| 国产熟女欧美一区二区| 欧美日韩在线观看h| 国产男人的电影天堂91| 久久99热这里只频精品6学生| 久久久色成人| 韩国av在线不卡| 国产乱人偷精品视频| 波多野结衣巨乳人妻| 精品不卡国产一区二区三区| 我的老师免费观看完整版| 亚洲精品日韩av片在线观看| 日韩av免费高清视频| 成人午夜高清在线视频| 亚洲欧美清纯卡通| 亚洲精品456在线播放app| 亚洲av成人精品一区久久| 亚洲性久久影院| 少妇高潮的动态图| 亚洲激情五月婷婷啪啪| 亚洲综合色惰| 日本免费a在线| 熟女人妻精品中文字幕| 久久精品国产亚洲网站| 国产综合懂色| 亚洲人成网站在线播| 国产亚洲av片在线观看秒播厂 | 亚洲最大成人中文| 久久久午夜欧美精品| 亚州av有码| 国产色爽女视频免费观看| .国产精品久久| 国产精品一及| 97精品久久久久久久久久精品| 在线播放无遮挡| 精品人妻视频免费看| 亚洲国产精品成人综合色| 天堂影院成人在线观看| 午夜日本视频在线| 嫩草影院新地址| 色尼玛亚洲综合影院| 精品国产露脸久久av麻豆 | 国产黄片视频在线免费观看| 日日摸夜夜添夜夜爱| 一个人观看的视频www高清免费观看| 久久精品国产亚洲网站| 高清在线视频一区二区三区| 国产色婷婷99| 永久免费av网站大全| 国产一区二区亚洲精品在线观看| 亚洲一区高清亚洲精品| 老司机影院毛片| 亚洲av福利一区| 51国产日韩欧美| 搞女人的毛片| 五月玫瑰六月丁香| 久久99热6这里只有精品| 亚洲欧美日韩东京热| 久久这里只有精品中国| 日本一二三区视频观看| 亚洲,欧美,日韩| 草草在线视频免费看| 一本一本综合久久| 国内精品一区二区在线观看| 中文精品一卡2卡3卡4更新| 三级毛片av免费| 天天一区二区日本电影三级| 美女国产视频在线观看| 99热全是精品| 老司机影院毛片| 国产av不卡久久| 青春草国产在线视频| 伦理电影大哥的女人| 日韩欧美三级三区| 超碰97精品在线观看| 免费av不卡在线播放| 午夜激情久久久久久久| 国产一级毛片在线| 午夜福利网站1000一区二区三区| 成人漫画全彩无遮挡| 亚洲久久久久久中文字幕| 国产精品久久久久久久电影| 搡老妇女老女人老熟妇| 亚洲三级黄色毛片| 亚洲欧洲国产日韩| 久久久久性生活片| 亚洲电影在线观看av| 国产视频首页在线观看| 老女人水多毛片| 亚洲精品日韩在线中文字幕| 国产精品久久视频播放| kizo精华| 国产精品人妻久久久影院| 久久精品国产亚洲网站| 国产午夜福利久久久久久| 91久久精品国产一区二区三区| 国产一级毛片在线| 在线观看av片永久免费下载| 精品久久久久久久末码| av专区在线播放| 欧美日韩精品成人综合77777| 性色avwww在线观看| 看黄色毛片网站| 男女下面进入的视频免费午夜| 午夜免费激情av| 亚洲精品一二三| 草草在线视频免费看| 亚洲欧美中文字幕日韩二区| 免费av不卡在线播放| 亚洲国产欧美人成| 亚洲国产成人一精品久久久| 久久精品人妻少妇| 亚洲欧美清纯卡通| 少妇高潮的动态图| 听说在线观看完整版免费高清| 亚洲精品国产av蜜桃| 精品一区二区三区人妻视频| av在线天堂中文字幕| 日韩亚洲欧美综合| 日本黄色片子视频| 日韩欧美国产在线观看| 十八禁网站网址无遮挡 | 床上黄色一级片| 精品欧美国产一区二区三| 日本三级黄在线观看| 久久鲁丝午夜福利片| 亚洲欧美日韩无卡精品| 在线观看美女被高潮喷水网站| 久久99热6这里只有精品| 国产精品不卡视频一区二区| 亚洲真实伦在线观看| 亚洲va在线va天堂va国产| 久久亚洲国产成人精品v| www.色视频.com| 免费黄色在线免费观看| 免费观看av网站的网址| 日韩成人av中文字幕在线观看| 国产一区有黄有色的免费视频 | 人妻少妇偷人精品九色| 精品不卡国产一区二区三区| 一区二区三区四区激情视频| 小蜜桃在线观看免费完整版高清| 黑人高潮一二区| 韩国高清视频一区二区三区| 午夜福利高清视频| 国产亚洲精品久久久com| 国产伦一二天堂av在线观看| 97在线视频观看| 蜜桃久久精品国产亚洲av| 天天躁夜夜躁狠狠久久av| 亚洲精品影视一区二区三区av| 国产伦精品一区二区三区四那| 成人国产麻豆网| 乱系列少妇在线播放| 国产精品.久久久| 男女国产视频网站| 国产成人a区在线观看| 欧美日韩在线观看h| 精品欧美国产一区二区三| 国产毛片a区久久久久| 亚洲av福利一区| 亚洲国产色片| 久久精品熟女亚洲av麻豆精品 | 婷婷六月久久综合丁香| av卡一久久| 老师上课跳d突然被开到最大视频| 熟女电影av网| 国产精品一区二区三区四区久久| 国产片特级美女逼逼视频| 青春草视频在线免费观看| 午夜福利高清视频| 青春草亚洲视频在线观看| 亚洲精品456在线播放app| 白带黄色成豆腐渣| 日本爱情动作片www.在线观看| 日本黄色片子视频| 国产一区二区三区综合在线观看 | 亚洲欧美清纯卡通| 午夜久久久久精精品| 成人亚洲欧美一区二区av| 丰满人妻一区二区三区视频av| 久久国内精品自在自线图片| 久久久欧美国产精品| av国产免费在线观看| 嫩草影院入口| 亚洲va在线va天堂va国产| 一区二区三区乱码不卡18| 97精品久久久久久久久久精品| 午夜久久久久精精品| 身体一侧抽搐| 欧美bdsm另类| 偷拍熟女少妇极品色| 国产亚洲av片在线观看秒播厂 | 精品一区二区免费观看| 日本一二三区视频观看| 免费观看精品视频网站| 亚洲精品久久久久久婷婷小说| 天美传媒精品一区二区| 久久久久久久大尺度免费视频| 成人亚洲精品一区在线观看 | 国产精品一二三区在线看| 精品少妇黑人巨大在线播放| 久久精品夜夜夜夜夜久久蜜豆| av国产久精品久网站免费入址| 亚洲av成人精品一区久久| 国产精品综合久久久久久久免费| 成人漫画全彩无遮挡| 女人被狂操c到高潮| 麻豆成人av视频| 久久久成人免费电影| 久久久久网色| 中文字幕久久专区| 男女边摸边吃奶| 少妇裸体淫交视频免费看高清| 精品99又大又爽又粗少妇毛片| 国产免费福利视频在线观看| 国产精品三级大全| 嘟嘟电影网在线观看| 人人妻人人澡人人爽人人夜夜 | 国产精品久久久久久精品电影小说 | 欧美日韩精品成人综合77777| 欧美xxxx黑人xx丫x性爽| 99热这里只有是精品50| 日韩av在线免费看完整版不卡| 日本av手机在线免费观看| 国产在视频线在精品| 寂寞人妻少妇视频99o| 久久精品综合一区二区三区| 99久国产av精品| 一个人看的www免费观看视频| 亚洲精品色激情综合| 麻豆精品久久久久久蜜桃| 夫妻性生交免费视频一级片| 黑人高潮一二区| 久久人人爽人人片av| 99re6热这里在线精品视频| 国产成人免费观看mmmm| 亚洲国产色片| 国产精品99久久久久久久久| 美女主播在线视频| 国产人妻一区二区三区在| 欧美激情国产日韩精品一区| 亚洲av男天堂| 熟女人妻精品中文字幕| 久久久久久久久久成人| 高清日韩中文字幕在线| 亚洲成人久久爱视频| 成人高潮视频无遮挡免费网站| 丰满人妻一区二区三区视频av| 亚洲精品第二区| 老司机影院毛片| 91久久精品国产一区二区三区| 内地一区二区视频在线| 国产免费一级a男人的天堂| 色综合站精品国产| 91精品国产九色| 全区人妻精品视频| 亚洲av电影在线观看一区二区三区 | 婷婷色综合www| 国产又色又爽无遮挡免| 国产精品av视频在线免费观看| 午夜爱爱视频在线播放| 国产一区二区三区综合在线观看 | 国产精品麻豆人妻色哟哟久久 | 中文字幕制服av| 精品一区二区三区人妻视频| 午夜亚洲福利在线播放| 国产成人精品一,二区| 精品人妻视频免费看| av女优亚洲男人天堂| 成人高潮视频无遮挡免费网站| 日韩av免费高清视频| 欧美成人午夜免费资源| 日韩欧美 国产精品| 一区二区三区四区激情视频| 丝袜喷水一区| 久久国产乱子免费精品| 久久国内精品自在自线图片| 80岁老熟妇乱子伦牲交| 国内揄拍国产精品人妻在线| 午夜激情欧美在线| 人妻制服诱惑在线中文字幕| 亚洲av福利一区| 晚上一个人看的免费电影| 永久免费av网站大全| 日本免费a在线| 成人综合一区亚洲| 少妇人妻一区二区三区视频| 春色校园在线视频观看| 亚洲国产日韩欧美精品在线观看| 国产精品久久久久久精品电影| 青春草国产在线视频| 你懂的网址亚洲精品在线观看| 听说在线观看完整版免费高清| 极品少妇高潮喷水抽搐| 久久99精品国语久久久| 97热精品久久久久久| 欧美高清成人免费视频www| 久久99热这里只有精品18| 中文乱码字字幕精品一区二区三区 | 男女边摸边吃奶| 亚洲aⅴ乱码一区二区在线播放| 联通29元200g的流量卡| 欧美bdsm另类| 国产麻豆成人av免费视频| 婷婷六月久久综合丁香| 搡老妇女老女人老熟妇| 大又大粗又爽又黄少妇毛片口| eeuss影院久久| 亚洲人与动物交配视频| 网址你懂的国产日韩在线| 国产精品日韩av在线免费观看| 国产91av在线免费观看| 91午夜精品亚洲一区二区三区| 18禁在线无遮挡免费观看视频| 免费av观看视频| 国产精品国产三级国产av玫瑰| a级毛色黄片| 天堂中文最新版在线下载 | 欧美精品国产亚洲| 久久久久久久大尺度免费视频| 久久久久久九九精品二区国产| 少妇猛男粗大的猛烈进出视频 | 免费大片黄手机在线观看| 亚洲自偷自拍三级| 亚洲久久久久久中文字幕| 国产三级在线视频| 亚洲va在线va天堂va国产| 伊人久久精品亚洲午夜| 97在线视频观看| 亚洲精品一区蜜桃| 久久99蜜桃精品久久| 男人爽女人下面视频在线观看| 亚洲欧美日韩东京热| 日本猛色少妇xxxxx猛交久久| 久久精品国产鲁丝片午夜精品| 国产av在哪里看| 男插女下体视频免费在线播放| 特级一级黄色大片| 亚洲欧洲日产国产| 乱码一卡2卡4卡精品| 2021天堂中文幕一二区在线观| 成人午夜精彩视频在线观看| 麻豆精品久久久久久蜜桃| 日本欧美国产在线视频| 夜夜爽夜夜爽视频| 午夜福利网站1000一区二区三区| 国产精品女同一区二区软件| 中文字幕制服av| 国产成年人精品一区二区| 国模一区二区三区四区视频| 免费黄频网站在线观看国产| 男人舔奶头视频| 2021少妇久久久久久久久久久| 超碰97精品在线观看| 99久久九九国产精品国产免费| 国产男人的电影天堂91| 老师上课跳d突然被开到最大视频| 亚洲精品一区蜜桃| 欧美bdsm另类| 国产综合懂色| 一级毛片 在线播放| 国产日韩欧美在线精品| 天堂√8在线中文| 久久久精品免费免费高清| 中文字幕人妻熟人妻熟丝袜美| 一个人看的www免费观看视频| 久久久久免费精品人妻一区二区| 美女xxoo啪啪120秒动态图| 久久精品国产亚洲网站| 日韩欧美 国产精品| 久久午夜福利片| 日韩,欧美,国产一区二区三区| 爱豆传媒免费全集在线观看| a级毛色黄片| 男女边摸边吃奶| 国产精品一区www在线观看| 高清午夜精品一区二区三区| 搡老妇女老女人老熟妇| 夜夜爽夜夜爽视频| av国产免费在线观看| 国产黄色视频一区二区在线观看| 色综合色国产| 亚洲最大成人av| 精品人妻一区二区三区麻豆| 国产av在哪里看| 在线 av 中文字幕| 免费无遮挡裸体视频| 国产精品精品国产色婷婷| 一级毛片aaaaaa免费看小| 亚洲怡红院男人天堂| 如何舔出高潮| 91aial.com中文字幕在线观看| 亚洲av免费在线观看| 中文乱码字字幕精品一区二区三区 | eeuss影院久久| 色尼玛亚洲综合影院| 一级毛片 在线播放| 欧美潮喷喷水| 免费观看性生交大片5| 老司机影院毛片| 国产精品熟女久久久久浪| 亚洲av中文av极速乱| 极品教师在线视频| 一级毛片 在线播放| 午夜视频国产福利| 两个人视频免费观看高清| 三级男女做爰猛烈吃奶摸视频| 五月天丁香电影| 日日啪夜夜撸| 高清毛片免费看| 亚洲熟妇中文字幕五十中出| 久久久欧美国产精品| 嫩草影院精品99| 国产成人精品福利久久| 久久久色成人| 国产一区二区三区av在线| 精品人妻一区二区三区麻豆| 成人特级av手机在线观看| 九草在线视频观看| 国内精品宾馆在线| 黄片无遮挡物在线观看| 最近最新中文字幕大全电影3| 一级毛片电影观看| 99久久九九国产精品国产免费| 一本一本综合久久| 亚洲自拍偷在线| av黄色大香蕉| 狂野欧美激情性xxxx在线观看| 国产精品一区二区三区四区免费观看| 在现免费观看毛片| 欧美日韩在线观看h| 中文字幕人妻熟人妻熟丝袜美| 激情五月婷婷亚洲| 成年版毛片免费区| 成人av在线播放网站| 美女大奶头视频| 午夜久久久久精精品| 久久久精品94久久精品| 久久久久久久久久久免费av| 国产人妻一区二区三区在| 美女大奶头视频| 99热这里只有是精品50| 午夜亚洲福利在线播放| 亚洲第一区二区三区不卡| 韩国高清视频一区二区三区| 亚洲av.av天堂| 91精品一卡2卡3卡4卡| 美女内射精品一级片tv| 麻豆久久精品国产亚洲av| 成人鲁丝片一二三区免费| 一级毛片 在线播放| 91精品一卡2卡3卡4卡| 能在线免费看毛片的网站| 中文资源天堂在线| 成人国产麻豆网| 成人综合一区亚洲| 久久6这里有精品| 1000部很黄的大片| 18禁动态无遮挡网站| 欧美精品国产亚洲| 久久鲁丝午夜福利片| 亚洲av日韩在线播放| 国产亚洲一区二区精品| 色吧在线观看| 国产午夜精品论理片| 高清欧美精品videossex| 日本与韩国留学比较| 国产综合懂色| 日本免费a在线| 日本黄大片高清| 91精品国产九色| 亚洲av成人av| 欧美激情久久久久久爽电影| 简卡轻食公司| 亚洲怡红院男人天堂| kizo精华| 最近手机中文字幕大全| 亚洲自拍偷在线| 午夜精品一区二区三区免费看| 校园人妻丝袜中文字幕| 纵有疾风起免费观看全集完整版 | 国产高潮美女av| 中文字幕制服av| 夫妻午夜视频| .国产精品久久| 青春草国产在线视频| 亚洲av.av天堂| 久久久久久久久久成人| 国产精品嫩草影院av在线观看| 日本午夜av视频| 国产淫语在线视频| 国产免费视频播放在线视频 | 精品99又大又爽又粗少妇毛片| 一边亲一边摸免费视频| 白带黄色成豆腐渣| 少妇人妻精品综合一区二区| 精品一区二区免费观看| 日日摸夜夜添夜夜添av毛片| 美女cb高潮喷水在线观看| 少妇人妻一区二区三区视频| 联通29元200g的流量卡| 毛片一级片免费看久久久久| freevideosex欧美| 一级爰片在线观看| 最近中文字幕高清免费大全6| 永久免费av网站大全| 成年av动漫网址| 久久久精品94久久精品| 你懂的网址亚洲精品在线观看| 国产成年人精品一区二区| 国产成人freesex在线| 国产人妻一区二区三区在| 久久精品夜夜夜夜夜久久蜜豆| 日韩在线高清观看一区二区三区| 婷婷色综合www| 免费观看a级毛片全部| 亚洲精品一区蜜桃| 青青草视频在线视频观看| 国产女主播在线喷水免费视频网站 | 午夜福利高清视频| 久久精品国产鲁丝片午夜精品| av在线天堂中文字幕| 日本黄大片高清| 丝瓜视频免费看黄片| 免费看美女性在线毛片视频| 国产免费福利视频在线观看| 欧美成人午夜免费资源| av在线天堂中文字幕| 边亲边吃奶的免费视频| 亚洲精品一区蜜桃| 一夜夜www| 久久精品人妻少妇| 在线观看人妻少妇| 免费观看无遮挡的男女| 国产极品天堂在线| 国产高清三级在线| 91精品伊人久久大香线蕉| 亚洲真实伦在线观看| 啦啦啦韩国在线观看视频| 国产成人a∨麻豆精品| 国产三级在线视频|