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

    Hardware-in-loop adaptive neural control for a tiltable V-tail morphing aircraft

    2023-05-06 01:59:58FuxingQioJingpingShiXioboQuYongxiLyu
    Defence Technology 2023年4期

    Fu-xing Qio ,Jing-ping Shi ,b,* ,Xio-bo Qu ,b ,Yong-xi Lyu ,b

    a School of Automation,Northwestern Polytechnical University,Xi'an,710129,China

    b Shaanxi Province Key Laboratory of Flight Control and Simulation Technology,Xi'an,710129,China

    Keywords:Morphing aircraft Back-stepping control Adaptive control Neural networks Radial basis function

    ABSTRACT This paper proposes an adaptive neural control(ANC)method for the coupled nonlinear model of a novel type of embedded surface morphing aircraft which has a tiltable V-tail.A nonlinear model with sixdegrees-of-freedom is established.The first-order sliding mode differentiator (FSMD) is applied to the control scheme to avoid the problem of “differential explosion”.Radial basis function neural networks are introduced to estimate the uncertainty and external disturbance of the model,and an ANC controller is proposed based on this design idea.The stability of the proposed ANC controller is proved using Lyapunov theory,and the tracking error of the closed-loop system is semi-globally uniformly bounded.The effectiveness and robustness of the proposed method are verified by numerical simulations and hardware-in-the-loop (HIL) simulations.

    1.Introduction

    For the past two decades,morphing aircraft has become an interesting topic due to its superior flight performance.Compared to the traditional fixed-wing aircraft,this aircraft has multi-mission adaptability by altering the aerodynamic structure of the aircraft[1—3].With the development of several research fields such as actuators,materials,structure,and control theory,morphing aircraft can easily and effectively implement multiple missions in different flight environments [4—8].As a complicated system,the control of morphing aircraft is faced with the huge variations of aerodynamic parameters,the existence of external disturbances,and internal uncertainties,which is one of the major hot spots of research on morphing aircraft [9].

    In previous research,most of the extant control schemes mainly focus on the decoupled longitudinal model,and most of the control objects in this literature are the sweep-back wings aircraft which is similar to the F-14 Tomcat.The aircraft with morphing wing technology mentioned in many studies are the Firebee [10]and the MFX-2 [11].Based on the unmanned aerial vehicle (UAV) Firebee,Yan et al.[12]designed a morphing aircraft that has a pair of variable-sweep wings.And designed a controller by proportionalintegral-derivative (PID) with different gains for the longitudinal model of this aircraft.Then,an adaptive super-twisting algorithm sliding mode control method was presented based on their existing research work [13].According to the longitudinal dynamics of the MFX-2,Wu et al.[14]proposed an adaptive approach based on neural network and high-order integral chained differentiator.Then they proposed a control scheme by composite switching neural prescribed-performance [15].And for the longitudinal dynamics with input-output constraints and unknown parameters,a modified control method by adaptive neural dynamic surface has been developed based on their previous work [16].Based on the same longitudinal model as Wu et al.,Wang et al.[17]modeled the longitudinal altitude motion as a switched nonlinear system in the form of a lower triangle.And the control scheme was proposed for it based on adaptive dynamic programming and switched nonlinear systems.

    In addition to the morphing aircraft with swept-back wings,some scholars have designed different types of morphing aircraft and proposed corresponding control methods.Inspired by a seagull's wing structure,Grant et al.[18]designed a multiple-joint morphing-sweep aircraft,which is independent on the wings along with on the inner and outer segments.And they analyzed the flight dynamics under the different swept-back wing configurations in symmetric and asymmetric cases.In Ref.[19],they advanced their work by studying the longitudinal dynamics of this aircraft and presented a time-varying characteristic equation for the influence of various morphing schemes.Based on X—24B,Jiao et al.[20]researched a hypersonic morphing vehicle with retracted winglets.According to the longitudinal model of this aircraft,they proposed a new learning approach for the control law of mode switching based on a type-2 Takagi-Sugeno-Kang fuzzy neural network.Yue et al.[21]investigated a folding-wing morphing aircraft and used the linear parameter varying (LPV) approach to build a longitudinal model of the aircraft and simulated the longitudinal dynamic response during wing folding.Based on the LPV model,they designed a pair of linear controllers by the selfscheduled gain H∞r(nóng)obust control approach to ensure the control effect of aircraft altitude and speed during the process of wing folding [22].

    Despite the designs and control schemes of morphing aircraft have made significant progress in recent years,some phenomena and issues need to be considered:

    (1) Many researchers have focused on the sweep-back-wing morphing aircraft as the control object,and the transformation process of this type of morphing aircraft only has notable changes in the longitudinal dynamics.

    (2) And most of the control schemes and performance analyses are based on the decoupled longitudinal model,however,there is not much analysis and research on the lateral dynamics of the morphing aircraft.Essentially,the controlled object is a single input single output model,it is difficult to reflect the complexity of morphing aircraft and the actual performance of the controller.And most previous researches have only considered the control of the longitudinal system,rather than considering longitudinal and lateral dynamics as a whole system.

    (3) In the field of morphing aircraft control,few scholars carry out HIL simulations after numerical simulations,they did not verify the control method by HIL simulation in a complex environment.However,it is very common to use HIL simulation to verify the feasibility and reliability of the algorithm in some fields such as navigation and vehicle control[23—25].

    Inspired by the mentioned discussion,this paper presents a new type of embedded surface morphing aircraft with an adjustable Vtail and proposes an ANC method for the coupled nonlinear model of this novel morphing aircraft.The main contributions of this work are listed as follows:

    (1) A novel morphing aircraft which is a tailless blended wing/body aircraft with an adjustable V-tail is designed to satisfy the requirements of different flight environments and missions.Based on a large number of aerodynamic data from different V-tail configurations,a nonlinear model with six degrees of freedom is established.

    (2) The FSMD is introduced into the traditional back-stepping control (TBC) scheme to avoid the “differential explosion”.To solve the problem of the external disturbance and uncertainty of the model,radial basis function (RBF) neural networks (NNs) are introduced,and an ANC controller is proposed based on this design idea.

    (3) By contrast to studies that focus only on a linearized model with decoupled longitudinal dynamics,the controller presented in this paper is designed for the longitudinal and lateral dynamics of the whole system.And the stability of the closed-loop system is proved by using Lyapunov synthesis based on stability analysis.

    (4) After numerical simulations,the robustness and effectiveness of the proposed method are verified by HIL simulations.

    The layout of this paper is organized as follows.Section 2 introduces the preliminaries and problem formulation.And section 3 presents an adaptive neural controller and the corresponding stability analysis.The numerical simulations and HIL simulations are given in section 4.Finally,section 5 concludes the work of this study.

    2.Preliminaries and problem formulation

    2.1.Morphing aircraft specification

    A novel type of embedded surface morphing aircraft,named XQ-8A which is shown in Fig.1(a),has been devised by a research group at the Northwestern Polytechnical University,Xi'an.The takeoff weight of the XQ-8A is 28 kg and the cruise speed is designed to be 150 km/s at an altitude of 1000 m.The key parameters of this aircraft are given in Table 1.And the aerodynamic data of the aircraft are shown in Fig.2.

    Table 1 Main parameters of XQ-8A.

    Table 2 Parameters of the ADRC.

    All control surfaces on the XQ-8A are shown in Fig.1(b).It has four empennage surfaces and a pair of variable flight control surfaces (VFCSs),which can form a V-tail or be embedded in the fuselage.And there are twin all-moving fins on the VFCSs.

    The control surfaces are embedded in the fuselage and the aircraft has a flying-wing configuration when the angle of the V-tail is zero (Λ=0°).Because of the mechanical dead zone,the allmoving fins cannot deflect when the angle of the V-tail is less than 15°(Λ<15°).Thus,only when the angle of the V-tail is greater than or equal to 15°(Λ ≥15°) do the all-moving fins work,operating as a rudder.The other two pairs of empennage surfaces can be divided into an inner control surface (CS1) and an outer control surface (CS2) according to their positions.The range of the V-tail angle (Λ) is 0°-90°.And the range of flight control surfaces CS1(δcs1R,δcs1L),and CS2 (δcs2R,δcs2L) is limited to -20°-20°.The position limit of the all-moving fins (δvcsR,δvcsL) is -15°~15°when the angle of the V-tail exceeds 15°(Λ ≥15°).

    Fig.1.Morphing aircraft XQ-8A: (a) Different morphing configurations;(b) Definition of control surface.

    Fig.2.Some aerodynamic data of XQ-8A: (a) Lift coefficient;(b) Drag coefficient;(c) Pitching moment coefficient;(d) Yawing moment coefficient.

    2.2.Dynamic model

    The control-oriented model of the dynamics of this aircraft are as follows:

    For simplicity,the dynamic model in Eq.(2) and Eq.(3) can be converted as

    Assumption 1.The deflection of the control surface has negligible effects on the aerodynamic force component [26,27],i.e.,g1u(xm)u≈0.

    Considering Assumption 1,the modeling uncertainty and external interference,the dynamic model can be written as

    where parts(a)and(b)in Eq.(8)denote two separate subsystems.The functions fi(?)and gi(?)are given in Appendix A.Δ1?R3and Δ2?R3denote the approximation errors in each subsystem,respectively.To deal with the system uncertainty in each subsystem effectively,RBFNNs are adopted.

    Assumption 2.The command signal y1d=x1d=[αdβdμd]Tis bounded and continuously differentiable,and its first derivative(1d) is also bounded.

    Assumption 3.The matrixes g1(xm)?R3×3and g2(xn)?R3×6are invertible.

    Remark 1.Note that g1u(xm)denotes the aerodynamic force component caused by the control surface deflection and numerical analysis finds that the control surfaces have a negligible influence on the aerodynamic force,so the size of g1u(xm)is negligible.Assumption 2 ensures that the command signal for aircraft control satisfies the conditions of continuous differentiability and a bounded first derivative.Because the control surfaces of the aircraft are designed to control each axis angular rate of aircraft independently,g2(xn)is invertible for all cases.The morphing aircraft in this paper has six control surfaces,so g2(xn)?R3×6is a matrix of three rows and six columns,and g2(xn)-1?R6×3is a pseudo-inverse matrix.

    2.3.Neural networks

    RBFNNs are used to approximate the continuous function [28].An arbitrary continuous function ⊿?R can be denoted as follows:

    where ω?RNis the weight matrix of the NN (N>1 denotes the number of neural network nodes,andN=50 in this paper),ρ?RNis the RBF matrix,and ε?Rexpresses the approximation error of the NN.

    Definition 1.The function trv(A)is defined to obtain the vectors of diagonal elements of the matrix A?RM×M.trv(A)and A are defined as follows:

    In Eq.(8),Δ1and Δ2denote the approximation errors of the two subsystems,and are defined as follows:

    The RBF ρi,j(xini)usually takes the form of a Gaussian function

    where ζiniis the center vector of the RBF and σi,jis thejth RBF bandwidth.The specific definitions of xiniand ζiniare as follows:

    Taking (ω1)Tρ1(xin1)as an example,the structure of the NN is shown in Fig.3.

    2.4.Active disturbance rejection control

    Active disturbance rejection control (ADRC) has achieved satisfactory performance for the uncertain nonlinear system [29,30].The ADRC approach is applied to design the velocity controller in this paper.ADRC comprises 3 parts: the extended state observer(ESO),the tracking differentiator(TD),and the nonlinear state error feedback(NLSEF).ADRC approach is applied to design the velocity controller,the structural diagram of velocity control is shown in Fig.4.

    The algorithm of TD is expressed as follow:

    ESO is used to observe the internal state and estimate the unpredictable state of the system.It is the core of ADRC.Form as follows:

    3.Controller design and stability analysis

    3.1.Controller design

    The proposed controller is designed for the attitude subsystem,which contains a two-step recursive design procedure.The FSMD is introduced into the TBC scheme to avoid the problem of “differential explosion”and RBFNNs are used to solve the problem of the external disturbance and uncertainty of the model.The diagram of the control scheme is shown in Fig.5.

    Step 1:Define e1=[e1,e2,e3]T?R3as the tracking error and y1dis the reference signal.The derivative of e1can be formulated as

    Define x2das the desired signal of the inner subsystem.The virtual control law of the outer system can then be designed as

    where K1?R3×3is a diagonal positive-definite matrix anddenote the estimations of W1.

    Fig.3.The structural diagram of NN.

    Lemma 1.The first-order sliding mode differentiator (FSMD) is designed as

    Step 2:To prevent the problem of the differential explosion,according to Lemma 1,we have

    Define e2=[e4,e5,e6]T?R3as the tracking error of the inner subsystem.So2can be expressed as

    where ηiis a positive design constant and Γi?RN×Nis a diagonal positive-definite matrix.

    3.2.Stability analysis

    Define the errors

    Then,note that the following inequalities hold

    Lemma 2.For any vectoraandb,the following inequality holds:

    Lemma 3.(Young’s inequality[32,33])The following inequality holds for any matrixAwhich‖A‖≤amax,amax>0.

    For convenience,in the following,abbreviations such as f1are used for f1(xm).According to the above equations,the time derivative of e1and e2are given by

    According to the matrix operation principle,the following expressions can be obtained

    Thus,Eq.(32) and Eq.(33) can be rewritten as

    Note thatg1min≤‖g1‖,ηi≤ηmaxand λmin(K1),λmin(K2),and λmin(Γi)denote the minimum eigenvalues of the corresponding matrices.Additionally,λmin(Γi)≥λminΓ(i=1,…,6).

    Therefore,the derivative ofLwis obtained as

    To ensure closed-loop stability,the design parametersC1should be selected such thatC1>0.According to Eq.(36),it can be shown that the tracking error of the closed-loop system is semi-globally uniformly bounded.

    4.Simulations and results

    4.1.Numerical simulations

    To illustrate the performance of the proposed method,two setsof numerical simulations are devised,and corresponding comparison and analysis are performed.

    The control parameters of the controllers are shown in Tables 2 and 3.For the RBFNNs component of the proposed controller,each control channel includes 50 NN nodes (N=50).The RBF bandwidths σi,j(i=1,2,3,j=1,…,50)in the first subsystem are evenly spaced in the scope [3,4].And the corresponding NN weight elements ωi,j(i=1,2,3,j=1,…,50)are initially set to random values in the scope[0,1].In the second subsystem,the initial values of the NN weight elements ωi,j(i=4,5,6,j=1,…,50)are also random in the range of[0,1].And the corresponding RBF bandwidths σi,j(i=4,5,6,j=1,…,50) are evenly spaced in the scope [20,21].

    Table 3 Parameters of the ANC.

    4.1.1.Verification under the static V-tail

    In this scenario,the angle of the V-tail is fixed at 45°,and the initial trim conditions are set as [V0,h0,α0,β0,μ0]=[41.67 m/s,1000 m,4.849°,0°,0°].The change of the reference signal is as follows:

    The tracking performance in this simulation is revealed in Fig.6(a).Fig.6(b)shows the root mean square errors (RMSEs)and the settling time.The results suggest that the proposed method decreases the RMSE and reduces the settling time under the condition of the static V-tail.

    4.1.2.Verification under the dynamic V-tail with wind and aerodynamics uncertainties

    In this scenario,a new numerical simulation is carried out to verify the efficiency and robustness of the proposed method.The Vtail angle is variable and the wind disturbance and aerodynamic coefficient variation are taken into consideration.The wind model is the Dryden Wind Turbulence Model in Simulink,and the wind velocities denoted in the three-body axes areuw,vw,andww,as revealed in Fig.7(b).The aerodynamic coefficient variation is expressed as:

    The V-tail command signal is shown in Fig.7(a) and the initial trim conditions are set as [V0,h0,α0,β0,μ0]=[41.67 m/s,1000 m,4.774°,0°,0°].As can be seen from Fig.8(a),the proposed controller has better robustness and performance and the RMSAEs revealed in Fig.8(b) support this conclusion.

    4.2.HIL simulations

    4.2.1.HIL simulation setup

    To further validate the performance and robustness of the proposed method,the corresponding methods are transferred to the HIL simulation.And three sets of HIL simulation scenarios are devised.

    The HIL simulation system contains the following components:(1) an autopilot with sensors for controller implementation and measuring angular velocity,(2)a real-time simulation computer for implement the morphing aircraft model,(3) a 3-axis turntable for driving the attitude motion of the autopilot,(4) a visual display software for displaying real-time flight data of morphing aircraft.Fig.9 Shows the structural diagram of the HIL simulation platform.

    Fig.6.Simulation results comparison in scenario 1: (a) Tracking performance;(b) RMSEs and settling time.

    Fig.7.V-tail command signal and wind velocity: (a) V-tail command signal;(b) Wind velocity.

    Fig.8.Simulation results comparison in scenario 2: (a) Tracking performance;(b) RMSEs.

    The autopilot was developed by our laboratory,which integrated a high-precision IMU integrated module for measuring the angular velocities of morphing aircraft.The real-time simulation computer runs the morphing aircraft model by MATLAB/xPC Target,a PC-based real-time control system from MathWorks [34,35].It receives the control signals from the autopilot and sends the states of morphing aircraft to autopilot.And it also sends the states to the visual display software via the UDP for observing the kinematic characteristics of the aircraft intuitively.The autopilot is fixed on a 3-axis turntable which was driving by the signals which are the angular velocities of morphing aircraft from the real-time simulation computer.Considering the real physical data source,the angular velocities used by autopilot comes from the data measured by IMU (sensor data,p,q,r) instead of the data calculated by the real-time computer (xPC data,p,q,r),and the comparison of the two data is shown in Fig.10.And Fig.11 shows the snapshot of the HIL simulation platform.

    4.2.2.Verification under the static V-tail

    In this scenario,the angle of the V-tail is fixed at 45°,and the initial conditions are set as[V0,h0,α0,β0,μ0]=[41.67 m/s,1000 m,4.849°,0°,0°].The reference signals are the same as scenario 1.The tracking performance is revealed in Fig.12(a).And the RMSEs are revealed in Fig.12(b).

    Compared with the previous numerical simulation results,the effect of HIL experiment becomes worse under the same control command.The main reason is the change of aircraft angular velocity data source of the controller.As shown in Fig.9,the angular velocities (xPC data,p,q,r) calculated by the real-time simulation computer drives the turntable to move,and the IMU of the autopilot is measuring the angular velocities (sensor data,p,q,r).Fig.10(a)shows the angular velocities under different information sources.It can be seen that there are errors and delays in the angular velocities used by the controller.Under such harsh signal interference,the two controllers can also track the reference signal well,which shows that the robustness of the controllers is good.As expected,the proposed method has a better tracking effect,especially on the α and μ channels.The comparison of RMSEs in Fig.12(b) and Table 4 supports this conclusion.

    4.2.3.Verification under the dynamic V-tail

    Fig.9.The architecture of the HIL simulation.

    Fig.10.Angular velocities under different information sources: (a) Data in scenario 3;(b) Data in scenario 4.

    In scenario 4,the V-tail command signal is the same as in scenario 2,and the initial trim conditions are set as [V0,h0,α0,β0,μ0]=[41.67 m/s,1000 m,4.774°,0°,0°].The angular velocities under different information sources are shown in Fig.10(b).The tracking performance and the RMSEs are shown in Fig.13.Table 4 shows the detailed comparison of RMSEs in this scenario.Scenario 4 adds the V-tail control command,which increases the interference to the system,so the control effect becomes worse compared with scenario 3.However,the results still support the previous conclusions,and it can be concluded from the detailed data in Table 4 that the greater the interference to the system,the better the effect of ANC controller compared to the TBS controller.

    Table 4 Comparison of the RMSEs of different HIL simulations.

    Although the obtained outcomes still support the previous conclusions,it can be observed through Fig.13(a)that around 1.7 s,the effects of both controllers become very poor.By comparing the V-tail control command in Fig.7(a),it can be found that the V-tail angle is less than 15°at 1.7 s before.Through analysis,it is found that this phenomenon is also following the objective laws of physics.Due to the mechanical dead zone,the all-moving fins cannot work when the V-tail is less than 15°.In this situation,the lateral stability and controllability of the aircraft become worse.And from the simulation in the 5 s of Fig.13(a),it can be noticed that the performance of the proposed controller is preferred over the TBS controller when Λ<15°.

    Fig.11.The snapshots of the HIL simulation platform.

    Fig.12.Simulation results comparison in scenario 3: (a) Tracking performance;(b) RMSEs.

    Leaving aside the effect of V-tail on the aircraft,a new experiment (scenario 5) is designed.In this scenario,the simulation conditions are the same as scenario 4 except for the V-tail command signal as shown in Fig.14(a).The angular velocities under different information sources are shown in Fig.14(b).The tracking performance and the RMSEs which are shown in Fig.15 and Table 4 also support the conclusion that the proposed method is superior.The next step to improve the simulation effect of the HIL experiment should start with the data source of the sensor and filter the state required by the controller.This will increase the effectiveness of the controller,which is also a way of our research in the future.

    5.Conclusions

    This paper presented a novel morphing aircraft that has an adjustable V-tail.An adaptive neural control method using RBFNNs was proposed and the longitudinal and lateral channel flight controller was designed by the proposed method.Using Lyapunov theory,the parameter ranges of the proposed method were derived and the Lyapunov stability for the proposed methods was proved.To verified the performance of the method in terms of rapidity,accuracy,and robustness,numerical and HIL simulations across five scenarios were designed.The results show that the proposed method works well with good robustness and tracking performance.

    Fig.13.Simulation results comparison in scenario 4: (a) Tracking performance;(b) RMSEs.

    Fig.14.New V-tail command signal and angular velocities: (a) New V-tail command signal;(b) Angular velocity.

    Fig.15.Simulation results comparison in scenario 5: (a)Tracking performance;(b)RMSEs.

    Declaration of competing interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgments

    This work was funded by the National Natural Science Foundation of China (No.61573286),the Natural Science Foundation of Shaanxi Province (2019JM-163,2020JQ-218),the Fundamental Research Funds for the Central Universities (3102019ZDHKY07),and supported by Shaanxi Province Key Laboratory of Flight Control and Simulation Technology.

    Appendix A

    where

    Note that abbreviations are used for aerodynamic coefficients,such asClsforCls(Λ).

    国产精品一区二区在线不卡| 一边亲一边摸免费视频| 亚洲美女视频黄频| 一本—道久久a久久精品蜜桃钙片| 99热国产这里只有精品6| 人妻 亚洲 视频| 日日摸夜夜添夜夜添av毛片| 色94色欧美一区二区| 91久久精品电影网| 草草在线视频免费看| 国产成人免费观看mmmm| 97超视频在线观看视频| 男女高潮啪啪啪动态图| 久久韩国三级中文字幕| 中文字幕最新亚洲高清| 久久青草综合色| 老司机亚洲免费影院| 秋霞在线观看毛片| 熟妇人妻不卡中文字幕| 18禁观看日本| 久久久久久人妻| 国产乱人偷精品视频| 国产高清三级在线| 男人添女人高潮全过程视频| 两个人免费观看高清视频| 午夜免费男女啪啪视频观看| 尾随美女入室| 波野结衣二区三区在线| 全区人妻精品视频| 中文天堂在线官网| 久久精品久久久久久噜噜老黄| 视频中文字幕在线观看| 亚洲av不卡在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久久久国产精品人妻一区二区| 亚洲精华国产精华液的使用体验| 交换朋友夫妻互换小说| 中文字幕最新亚洲高清| 麻豆乱淫一区二区| 亚洲在久久综合| 91国产中文字幕| 狂野欧美激情性xxxx在线观看| 国产 精品1| 久久精品熟女亚洲av麻豆精品| 日本黄大片高清| 一本色道久久久久久精品综合| 性色avwww在线观看| 亚洲精品一区蜜桃| 永久网站在线| 亚洲欧洲日产国产| 国产视频内射| 国产成人精品福利久久| 亚洲综合色惰| 夜夜骑夜夜射夜夜干| videos熟女内射| 欧美日韩综合久久久久久| 日韩不卡一区二区三区视频在线| 插阴视频在线观看视频| 夜夜爽夜夜爽视频| 久久免费观看电影| 色视频在线一区二区三区| 国产69精品久久久久777片| 人妻一区二区av| 成人手机av| 免费大片18禁| 国产深夜福利视频在线观看| 欧美三级亚洲精品| 国产精品人妻久久久影院| 国产亚洲精品久久久com| 在线观看人妻少妇| 纵有疾风起免费观看全集完整版| 国产成人精品婷婷| 欧美丝袜亚洲另类| 色94色欧美一区二区| 色婷婷av一区二区三区视频| 亚洲国产精品成人久久小说| 久久av网站| 欧美精品一区二区免费开放| 王馨瑶露胸无遮挡在线观看| 日韩视频在线欧美| 国产精品国产三级专区第一集| 国产日韩一区二区三区精品不卡 | 婷婷色综合www| 国产亚洲av片在线观看秒播厂| 99国产综合亚洲精品| 各种免费的搞黄视频| 日韩亚洲欧美综合| 91精品一卡2卡3卡4卡| 成人手机av| 久久青草综合色| 中文字幕免费在线视频6| 女性生殖器流出的白浆| 成人国语在线视频| 狂野欧美白嫩少妇大欣赏| kizo精华| 亚洲国产av新网站| a 毛片基地| 考比视频在线观看| 99热国产这里只有精品6| 日韩成人av中文字幕在线观看| 特大巨黑吊av在线直播| 天堂中文最新版在线下载| 又粗又硬又长又爽又黄的视频| 亚洲国产色片| 国产成人午夜福利电影在线观看| 欧美bdsm另类| 日韩亚洲欧美综合| 久久av网站| 美女主播在线视频| 一级毛片我不卡| 国产精品麻豆人妻色哟哟久久| 寂寞人妻少妇视频99o| 午夜久久久在线观看| 最近的中文字幕免费完整| 在线播放无遮挡| 免费观看a级毛片全部| 黄色毛片三级朝国网站| 成年女人在线观看亚洲视频| 亚洲天堂av无毛| 亚洲高清免费不卡视频| 天天影视国产精品| 成人手机av| 国产免费视频播放在线视频| 中文字幕久久专区| 18在线观看网站| xxxhd国产人妻xxx| 亚洲三级黄色毛片| 2022亚洲国产成人精品| 丰满少妇做爰视频| 日本91视频免费播放| 国产探花极品一区二区| 最新的欧美精品一区二区| 成人影院久久| 少妇 在线观看| 久久久精品免费免费高清| 成人毛片a级毛片在线播放| 亚洲欧美成人综合另类久久久| 黄色毛片三级朝国网站| 免费日韩欧美在线观看| 另类精品久久| 精品熟女少妇av免费看| 少妇被粗大猛烈的视频| 精品少妇久久久久久888优播| 一级黄片播放器| 亚洲人成网站在线观看播放| av网站免费在线观看视频| 久久久国产一区二区| 欧美 亚洲 国产 日韩一| 97在线人人人人妻| 成人亚洲精品一区在线观看| 国产 一区精品| 国产国拍精品亚洲av在线观看| 美女内射精品一级片tv| 最近的中文字幕免费完整| 免费黄频网站在线观看国产| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | av一本久久久久| 中文字幕精品免费在线观看视频 | 精品亚洲乱码少妇综合久久| 国产成人午夜福利电影在线观看| 美女大奶头黄色视频| 国产亚洲欧美精品永久| 一本一本综合久久| 国产乱人偷精品视频| 女性生殖器流出的白浆| 啦啦啦在线观看免费高清www| 曰老女人黄片| 看十八女毛片水多多多| av天堂久久9| 99热这里只有精品一区| 亚洲精品日韩在线中文字幕| 91精品伊人久久大香线蕉| 国产成人精品一,二区| 亚洲精品自拍成人| 亚洲精品国产av成人精品| 久久这里有精品视频免费| 两个人的视频大全免费| 最黄视频免费看| 国产又色又爽无遮挡免| 老司机影院毛片| 日本vs欧美在线观看视频| 国产免费一区二区三区四区乱码| 亚洲人成网站在线播| 夫妻性生交免费视频一级片| 日本av免费视频播放| 久久精品久久精品一区二区三区| 99久久人妻综合| 欧美最新免费一区二区三区| 亚洲国产精品专区欧美| 久久97久久精品| 亚洲精品成人av观看孕妇| 久久精品熟女亚洲av麻豆精品| 91久久精品国产一区二区三区| 男男h啪啪无遮挡| 亚洲成人一二三区av| 欧美变态另类bdsm刘玥| 国产成人精品久久久久久| 日日撸夜夜添| 高清不卡的av网站| 啦啦啦啦在线视频资源| 亚洲精品一二三| 久久久午夜欧美精品| 成人18禁高潮啪啪吃奶动态图 | 如何舔出高潮| www.av在线官网国产| 欧美精品人与动牲交sv欧美| 欧美国产精品一级二级三级| 日本黄色片子视频| 久久久久久久亚洲中文字幕| 日本91视频免费播放| 久久亚洲国产成人精品v| 免费观看的影片在线观看| 久久精品国产亚洲av天美| 另类亚洲欧美激情| 国产免费福利视频在线观看| 国产国拍精品亚洲av在线观看| 免费黄色在线免费观看| 午夜视频国产福利| 汤姆久久久久久久影院中文字幕| 制服人妻中文乱码| 三上悠亚av全集在线观看| 欧美最新免费一区二区三区| 亚洲欧洲国产日韩| 18禁在线播放成人免费| 欧美xxxx性猛交bbbb| 在线观看一区二区三区激情| 夫妻性生交免费视频一级片| 午夜av观看不卡| 黑人欧美特级aaaaaa片| 亚洲欧美日韩另类电影网站| 在线观看一区二区三区激情| 大陆偷拍与自拍| 草草在线视频免费看| 一区二区三区四区激情视频| 超碰97精品在线观看| kizo精华| 亚洲熟女精品中文字幕| 日韩免费高清中文字幕av| 伦精品一区二区三区| 国产精品人妻久久久影院| 如何舔出高潮| 亚洲av国产av综合av卡| 曰老女人黄片| 国产av码专区亚洲av| 国产成人av激情在线播放 | 久久久久久久久久久久大奶| 丝袜美足系列| 免费大片黄手机在线观看| 亚洲欧洲日产国产| 国产色婷婷99| 日韩大片免费观看网站| 亚洲国产精品国产精品| 中文欧美无线码| 国产片特级美女逼逼视频| 大香蕉久久网| 久久青草综合色| 80岁老熟妇乱子伦牲交| 久久久国产一区二区| 黑人巨大精品欧美一区二区蜜桃 | 99九九线精品视频在线观看视频| av视频免费观看在线观看| 日本wwww免费看| 久久久精品区二区三区| 亚洲av国产av综合av卡| 九色亚洲精品在线播放| 国产成人精品无人区| 亚洲高清免费不卡视频| 欧美激情国产日韩精品一区| 日本猛色少妇xxxxx猛交久久| av线在线观看网站| 婷婷色综合大香蕉| 欧美丝袜亚洲另类| 蜜臀久久99精品久久宅男| 国产亚洲精品久久久com| 国产在线免费精品| 观看美女的网站| 在线观看免费视频网站a站| 久久久久久伊人网av| 晚上一个人看的免费电影| 亚洲av日韩在线播放| 99九九线精品视频在线观看视频| 久久精品久久久久久久性| 国产免费现黄频在线看| 中文精品一卡2卡3卡4更新| 少妇人妻精品综合一区二区| 日本黄色日本黄色录像| 国产欧美另类精品又又久久亚洲欧美| 熟女人妻精品中文字幕| 成人亚洲精品一区在线观看| 欧美另类一区| 97在线人人人人妻| 国产精品无大码| 春色校园在线视频观看| a级片在线免费高清观看视频| 国产极品粉嫩免费观看在线 | 成年女人在线观看亚洲视频| 成人手机av| 亚洲婷婷狠狠爱综合网| 亚洲精华国产精华液的使用体验| 亚洲国产精品专区欧美| 少妇 在线观看| 伦理电影大哥的女人| 国产男人的电影天堂91| 国产69精品久久久久777片| .国产精品久久| av一本久久久久| 亚洲精品国产av蜜桃| 国产免费一级a男人的天堂| 天堂8中文在线网| 久久精品国产鲁丝片午夜精品| 熟女电影av网| 国产国拍精品亚洲av在线观看| 青春草国产在线视频| 久久av网站| 国产伦精品一区二区三区视频9| 国产精品久久久久成人av| 久久这里有精品视频免费| 中文欧美无线码| 大片免费播放器 马上看| 欧美一级a爱片免费观看看| 2022亚洲国产成人精品| 黄色毛片三级朝国网站| 国产免费又黄又爽又色| 在线观看免费日韩欧美大片 | 免费大片黄手机在线观看| 亚洲av电影在线观看一区二区三区| 男女无遮挡免费网站观看| 免费不卡的大黄色大毛片视频在线观看| 99久久人妻综合| 又大又黄又爽视频免费| 免费黄网站久久成人精品| 街头女战士在线观看网站| 在线免费观看不下载黄p国产| 中文天堂在线官网| 欧美xxxx性猛交bbbb| 交换朋友夫妻互换小说| 九色亚洲精品在线播放| 亚洲精品日韩在线中文字幕| 一边摸一边做爽爽视频免费| 国产黄频视频在线观看| 性色av一级| 欧美亚洲 丝袜 人妻 在线| 97精品久久久久久久久久精品| 天天操日日干夜夜撸| 国产亚洲精品久久久com| 插逼视频在线观看| av不卡在线播放| 美女脱内裤让男人舔精品视频| 丝袜美足系列| 熟女电影av网| 日韩亚洲欧美综合| av网站免费在线观看视频| av女优亚洲男人天堂| 亚洲欧美精品自产自拍| a级毛片黄视频| 午夜福利在线观看免费完整高清在| 丰满少妇做爰视频| 日本-黄色视频高清免费观看| 午夜老司机福利剧场| 99精国产麻豆久久婷婷| 欧美 亚洲 国产 日韩一| 肉色欧美久久久久久久蜜桃| 国产亚洲精品久久久com| 免费观看av网站的网址| 男人爽女人下面视频在线观看| 午夜福利在线观看免费完整高清在| 一级a做视频免费观看| 国产乱人偷精品视频| 极品人妻少妇av视频| 久久精品国产鲁丝片午夜精品| www.色视频.com| 丝袜在线中文字幕| 天天操日日干夜夜撸| 精品国产一区二区三区久久久樱花| 亚洲精品色激情综合| 婷婷色av中文字幕| 精品国产露脸久久av麻豆| 亚洲内射少妇av| 成年av动漫网址| 岛国毛片在线播放| 欧美国产精品一级二级三级| 国产精品蜜桃在线观看| 寂寞人妻少妇视频99o| 美女xxoo啪啪120秒动态图| 国产亚洲精品第一综合不卡 | 我的老师免费观看完整版| 久久 成人 亚洲| 国产不卡av网站在线观看| 日韩三级伦理在线观看| 亚洲图色成人| 国产精品偷伦视频观看了| 精品久久久精品久久久| 亚洲无线观看免费| 免费不卡的大黄色大毛片视频在线观看| 国产精品一二三区在线看| 日本免费在线观看一区| 国产视频首页在线观看| 精品亚洲成a人片在线观看| 亚洲美女黄色视频免费看| 日韩亚洲欧美综合| 丝袜喷水一区| 在线播放无遮挡| 国产精品国产三级国产av玫瑰| 中文字幕制服av| 亚洲精品自拍成人| 日本欧美国产在线视频| 精品人妻熟女毛片av久久网站| 最新中文字幕久久久久| 九色亚洲精品在线播放| 看免费成人av毛片| 国产亚洲最大av| 国产日韩欧美亚洲二区| 欧美人与性动交α欧美精品济南到 | 插阴视频在线观看视频| 色5月婷婷丁香| 亚州av有码| 欧美日韩亚洲高清精品| 精品人妻偷拍中文字幕| 国模一区二区三区四区视频| 一级毛片aaaaaa免费看小| 精品一区二区免费观看| 国产日韩欧美视频二区| 欧美日韩一区二区视频在线观看视频在线| 日韩成人av中文字幕在线观看| 丁香六月天网| 制服诱惑二区| a级毛片黄视频| 夫妻性生交免费视频一级片| 国产亚洲最大av| 国产色婷婷99| 精品久久久久久电影网| 国产无遮挡羞羞视频在线观看| 纵有疾风起免费观看全集完整版| av电影中文网址| 亚洲第一区二区三区不卡| 黄色视频在线播放观看不卡| 日本av免费视频播放| 成年av动漫网址| 色吧在线观看| 两个人免费观看高清视频| 久久综合国产亚洲精品| 亚洲人成77777在线视频| 免费观看的影片在线观看| 久久久久久久久久久久大奶| 亚洲精品中文字幕在线视频| 久久久久久伊人网av| 男女无遮挡免费网站观看| 久久午夜综合久久蜜桃| 久久久久久久久久久免费av| 国产精品嫩草影院av在线观看| 男女边摸边吃奶| 热re99久久精品国产66热6| 久久久久精品久久久久真实原创| 亚洲国产毛片av蜜桃av| av天堂久久9| 黄色一级大片看看| 一本大道久久a久久精品| 在线观看免费高清a一片| 国产精品一区二区在线观看99| 免费观看的影片在线观看| 人妻夜夜爽99麻豆av| 亚洲精品成人av观看孕妇| 3wmmmm亚洲av在线观看| 又黄又爽又刺激的免费视频.| 婷婷色综合大香蕉| 亚洲精品乱码久久久v下载方式| 色吧在线观看| 蜜桃在线观看..| 国产黄片视频在线免费观看| 99视频精品全部免费 在线| 丝袜脚勾引网站| 日本与韩国留学比较| 午夜精品国产一区二区电影| xxx大片免费视频| 久久毛片免费看一区二区三区| av国产精品久久久久影院| 国产精品久久久久久久电影| 少妇人妻精品综合一区二区| 国产国语露脸激情在线看| 免费少妇av软件| 99九九线精品视频在线观看视频| 波野结衣二区三区在线| 国产精品99久久99久久久不卡 | 国产乱人偷精品视频| 国产男女内射视频| 精品人妻熟女毛片av久久网站| 人妻人人澡人人爽人人| 日韩av在线免费看完整版不卡| 18禁动态无遮挡网站| 亚洲精品色激情综合| 日韩成人伦理影院| 久久韩国三级中文字幕| 亚洲av欧美aⅴ国产| 久热久热在线精品观看| 国产一区亚洲一区在线观看| 亚洲,欧美,日韩| 精品久久久精品久久久| 汤姆久久久久久久影院中文字幕| 久久精品国产a三级三级三级| 精品久久久久久久久av| 黄色配什么色好看| av女优亚洲男人天堂| 亚洲天堂av无毛| 亚洲国产毛片av蜜桃av| 欧美日韩视频高清一区二区三区二| 国产av精品麻豆| 亚洲精品第二区| 草草在线视频免费看| 午夜免费观看性视频| 国产av国产精品国产| 制服丝袜香蕉在线| 亚洲精品456在线播放app| 久久久久视频综合| 国产成人一区二区在线| 美女cb高潮喷水在线观看| 欧美日韩成人在线一区二区| videos熟女内射| 中国国产av一级| 精品国产一区二区三区久久久樱花| 少妇人妻精品综合一区二区| 亚洲国产av新网站| 热99久久久久精品小说推荐| 午夜91福利影院| 18在线观看网站| 亚洲少妇的诱惑av| 一级二级三级毛片免费看| 99精国产麻豆久久婷婷| 建设人人有责人人尽责人人享有的| 亚洲美女视频黄频| 国产又色又爽无遮挡免| 最后的刺客免费高清国语| 汤姆久久久久久久影院中文字幕| av黄色大香蕉| 久久综合国产亚洲精品| 天堂中文最新版在线下载| 婷婷成人精品国产| 久久久国产一区二区| 免费黄色在线免费观看| 99热国产这里只有精品6| 日韩精品免费视频一区二区三区 | 中文乱码字字幕精品一区二区三区| 高清黄色对白视频在线免费看| 亚洲成人一二三区av| 五月玫瑰六月丁香| 看免费成人av毛片| 国产精品99久久99久久久不卡 | av在线播放精品| 午夜91福利影院| .国产精品久久| 丰满少妇做爰视频| 男女免费视频国产| 国产黄色视频一区二区在线观看| 久久热精品热| 国产成人aa在线观看| 交换朋友夫妻互换小说| 久久久久精品性色| 亚洲一区二区三区欧美精品| 99热全是精品| 日本午夜av视频| 黑人巨大精品欧美一区二区蜜桃 | 69精品国产乱码久久久| 三上悠亚av全集在线观看| 久久久久久久久大av| 精品一区二区三卡| 亚洲精品国产色婷婷电影| 欧美激情国产日韩精品一区| 久久国产亚洲av麻豆专区| 插逼视频在线观看| 曰老女人黄片| 只有这里有精品99| 纯流量卡能插随身wifi吗| 亚洲av.av天堂| 男人添女人高潮全过程视频| 午夜久久久在线观看| 亚洲无线观看免费| 国产免费视频播放在线视频| 高清午夜精品一区二区三区| 一区二区三区四区激情视频| 亚洲av日韩在线播放| 精品99又大又爽又粗少妇毛片| 国产片内射在线| 国产精品久久久久久久电影| 欧美日韩亚洲高清精品| a级毛片免费高清观看在线播放| 人体艺术视频欧美日本| 嘟嘟电影网在线观看| 夫妻午夜视频| 九九在线视频观看精品| 国产精品久久久久久久久免| 欧美另类一区| 免费观看av网站的网址| 亚洲精品国产av蜜桃| 成人影院久久| 国产免费福利视频在线观看| 九九爱精品视频在线观看| 日韩强制内射视频| 国产免费福利视频在线观看| 亚洲不卡免费看| 久久人人爽av亚洲精品天堂| 中文天堂在线官网| 久久精品国产鲁丝片午夜精品| 天堂8中文在线网| 成人毛片a级毛片在线播放| 高清在线视频一区二区三区| 校园人妻丝袜中文字幕| 91午夜精品亚洲一区二区三区| 久久久a久久爽久久v久久| 免费观看无遮挡的男女| 欧美日韩成人在线一区二区| 大又大粗又爽又黄少妇毛片口| 亚洲婷婷狠狠爱综合网| 亚洲色图 男人天堂 中文字幕 | 国产精品.久久久| 99热这里只有是精品在线观看| 2018国产大陆天天弄谢|