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

    Sliding mode control design for oblique wing aircraft in wing skewing process

    2019-02-27 08:59:06TingYUEZijinXULixinWANGTongWANG
    CHINESE JOURNAL OF AERONAUTICS 2019年2期
    關(guān)鍵詞:信息

    Ting YUE,Zijin XU,Lixin WANG,*,Tong WANG

    aSchool of Aeronautic Science and Engineering,Beihang University,Beijing 100083,China

    bChina Ship Development and Design Center,Wuhan 430064,China

    Abstract When the wing of Oblique Wing Aircraft(OWA)is skewed,the center of gravity,inertia and aerodynamic characteristics of the aircraft all significantly change,causing an undesirable flight dynamic response,affecting the flying qualities,and even endangering the flight safety.In this study,the dynamic response of an OWA in the wing skewing process is simulated,showing that the three axis movements of the OWA are highly coupled and present nonlinear characteristics during the wing skewing.As the roll control efficiency of the aileron decreases due to the shortened control arm in an oblique configuration,the all-moving horizontal tail is used for additional roll and the control allocation is performed based on minimum control energy.Given the properties of pitchroll-yaw coupling and control input and state coupling,and the difficulty of establishing an accurate aerodynamic model in the wing skewing process due to unsteady aerodynamic force,a multi-loop sliding mode controller is formulated by the time-scale separation method.The closed-loop simulation results show that the asymmetric aerodynamics can be balanced and that the velocity and altitude of the aircraft maintain stable,which means that a smooth transition is obtained during the OWA's wing skewing.

    KEYWORDS Decoupling;Dynamic response;Oblique wing aircraft;Sliding mode control;Wing skewing process

    1.Introduction

    Oblique Wing Aircraft(OWA)can rotate the wing at different flight velocities to form various wing sweep configurations.Representing a simplified minimum drag solution in supersonic flow,the oblique elliptic wing can achieve efficient flight under subsonic cruise/loiter conditions while providing excellent supersonic dash/cruise capability.1,2Therefore,OWA can adapt to multi-mission flight with extended flight envelope,and demonstrate higher combat effectiveness than conventional fixed-wing aircraft.

    A significant amount of early theoretical work was conducted by Jones starting from the 1950s.3In the past decades,a number of aircrafts were flown to prove the feasibility of flight control of OWA,including small gliders,small remotely piloted aircraft and the manned AD-1 aircraft flown by NASA.4-8Numerous conceptual design studies and research papers have addressed oblique wing-body-tail or oblique flying wing designs,and a number of wind tunnel tests and flight tests have been performed on both oblique wing and oblique flying wing designs.9-17Recently,there has been renewed interest in developing OWA due to its aerodynamic and structural advantages.

    However,a variety of technical difficulties associated with the oblique wing configuration have prevented its application to operational aircraft.A critical technical challenge of OWA arises from its nonlinear and strongly coupled aerodynamic characteristics,and the inertial couplings further complicate the flight dynamics of such aircraft. These complicated characteristics are shown during the wing skewing of OWA,when a relative movement occurs between the wing and the fuselage.The aircraft in various configurations can be described as a time-varying dynamic system.Wing skewing largely changes the aircraft's configuration parameters,including the moment of inertia,the location of the center of gravity,etc.It not only causes changes of lift,drag and pitching moment,but also generates side force,and rolling and yawing moments.This results in the aircraft's highly coupled threeaxis movement and presents strong nonlinearity.In addition,wing skewing causes longitudinal and lateral couplings of aileron control,and lateral control inefficiency of the aileron at a large angle of the oblique wing.As a result,OWA will have an undesirable dynamic response in the presence of wing skews,and the flying qualities and flight safety will deteriorate.It is thus necessary to investigate the flight dynamic characteristics of OWA during wing skewing and design the flight control law to ensure a smooth transition between different configurations.

    The wing skewing process of OWA is similar to that of variable-sweep wing aircraft and morphing aircraft.Several studies have been carried out on the flight dynamic modelling and control of variable-sweep wing aircraft and morphing aircraft in configuration variations.18-23However,unlike the aircraft designed in a symmetric configuration,OWA is characterized by the asymmetry configuration with strong three-axis aerodynamic and inertial couplings,resulting in a significant variation in three-axis motion parameters during wing skewing.Currently,there have been a number of studies mainly focusing on the analysis of flight dynamic characteristics24-26and flight control design of OWA in static configurations(wing- fixed configurations at different oblique angles).27-32However,the existing literature does not consider the influence of the wing skewing process on aircraft's dynamic characteristics;additionally,there is lack of published literature about dynamic characteristics of and flight control design for the process of wing skewing.

    Due to the inertial and pitch-yaw-roll aerodynamic couplings in the wing skewing process,OWA is a non-linear time-varying system.This causes a difficulty in modeling precisely the aerodynamics of wing skewing which has uncertain dynamic behaviors.In addition,due to the inertial and aerodynamic cross-couplings,it is difficult for conventional control methods such as proportional-Integral-Derivative(PID)control to obtain acceptable control performance.A feasible approach for flight control design is sliding mode control theory,which can guarantee the dynamic qualities of the sliding mode motion and approaching motions by designing the sliding mode function and approaching rate,respectively.33,34

    In this study,OWA is regarded as a rigid body system composed of the fuselage and rotating wing,and a time-varying and nonlinear dynamic model is built in order to examine the flight dynamic characteristics of OWA in the wing skewing process.The dynamic responses are numerically simulated and analyzed.According to the highly pitch-roll-yaw coupling,control input and state coupling,and nonlinear properties of OWA,a multi-loop sliding mode controller is designed based on the time-scale separation method.The all-moving horizontal tail is used for roll control due to the insufficient control effect of the aileron,and the control allocation of the control surfaces is performed based on minimum control energy.The simulation results show that the controller can ensure the good dynamic characteristics of OWA during wing skewing.The work introduced in this paper can provide a theoretical reference for flight control design for OWA.

    2.Dynamic modeling and analysis

    2.1.Layout and aerodynamic characteristics of OWA

    The OWA investigated in this study is presented in Fig.1.The oblique wing is designed to rotate from 0°to 60°with the right wing forward.

    When the wing is skewed,the configuration parameters of the OWA change significantly,which mainly include the center of gravity position and the moment of inertia.Since the rotation axis of the wing is in front of the wing gravity center and parallel to thez-axis in the body axes,the aircraft's gravity center shifts forward to the right when the wing is skewed,while the gravity center position along thez-axis remains unchanged.Meanwhile,the inertia moment and product of inertia also change significantly.As the symmetry of the conventional layout is destroyed,the cross inertial productsIxyandIyzappear.The inertia coupling is strong between the longitudinal and lateral axes,which obviously affects the flight dynamic characteristics of the OWA.32

    The asymmetric layout formed by the skewed wing also has a great influence on aerodynamic characteristics,which are stated as follows.32The drag gradually decreases with the increase of the skew angle Λ.Since the effective wing span reduces after the wing is oblique,the lift decreases.In addition,the asymmetric wing layout generates side force,which is rather large and cannot be ignored.The oblique wing generates not only additional nose-down pitching moment,but also rolling moment and yawing moment.The rolling moment sinks the front wing and the yawing moment makes the aircraft rotate to the right.

    Fig.1 Layout of the oblique wing aircraft.

    After the wing is oblique,the ailerons are arranged back and forth.When the ailerons are deflected,the difference in lift will generate additional pitching moment.32This should be considered in designing flight control law for the wing skewing process.In addition,the control arm of the aileron is shortened after the wing is oblique,which gradually decreases the roll control efficiency of the aileron with the increase of the wing skew angle.This may result in insufficient roll control ability of the ailerons.Therefore,it is necessary to use the all moving horizontal tail for additional roll control of the OWA in the large skew angled configuration.27

    2.2.Dynamic modeling

    During the oblique wing skewing,the OWA's gravity center keeps changing.In this paper,the origin of the body axes is located at the center of gravity of the OWA at 0°skew angle.According to the momentum theorem and the theorem of moment of momentum for any moving point,18,35we get the nonlinear dynamic equations of the OWA in the wing skewing process in a scalar form as below:

    where[X,Y,Z]is the total force on the aircraft,[L,M,N]the total moment relative to the origin of the body axes,mthe mass of the aircraft,[u,v,w]the flight velocity,[p,q,r]the angular velocity in the body axes relative to the inertial axes,[Sx,Sy,Sz]the static moment of aircraft,[S1x,S1y,S1z]the static moments of the wing,mthe mass of the wing,Ithe rotational inertia,and[ω1x, ω1y, ω1z]the projections of the angular velocity of the wing in the body axes.It can be seen from Eqs.(1)and(2)that,unlike the dynamic equations of conventional fixed-wing aircraft,there are additional terms caused by the center of gravity position shift and wing skew rate.

    In the variable-sweep process of swept-wing aircraft and the wing-folding process of morphing aircraft,when the changing rate of configurations is slow,the influence of unsteady aerodynamic characteristics of the aircraft is small.18,21,22Analogously,when the angular velocity of the OWA's wing skew is comparatively small,the quasi-steady assumption can be used to simplify the calculation of aerodynamic force during the wing skewing process.It is considered that the aerodynamic force of the whole aircraft depends on the static configuration and the flight status of the aircraft at that time.In this study,we just investigate the case in which the angular velocity of the OWA's wing skew is comparatively slow,and in which the aerodynamic force is simplified using the quasi-steady assumption.The aerodynamic uncertainty caused by simplification is solved by robust sliding mode control design.Therefore,here we simplify that the aerodynamic model during the wing skewing is the same with the static configuration.32

    When the OWA is in an asymmetric configuration,the pitching motion will generate roll and yaw motions which will then generate pitching motion.Moreover,the deflection of the aileron will generate pitching moment.The aerodynamic crosscouplings are dependent on the wing skew angle and the angle of attack.As a result,the nonlinear equation of motion can be expressed as

    where f is a nonlinear time-varying function about x and u.The changes of the state variables affect the efficiency of the control input,and the control input can lead to changes of the state variables,i.e.,coupling exists between x and u and cannot be overlooked.This is the main difference of the six degrees of freedom nonlinear model with that of conventional aircraft.

    2.3.Analysis of dynamic characteristics

    The flight dynamics of the OWA in the wing skewing process without control is mathematically simulated.The aircraft's CG position of the straight wing configuration is selected as the origin of the body axes.The wing skewing process and returning to straight wing process are simulated separately.The initial condition is a straight and level flight atH=4000 m andMa=0.4.The dynamic responses are presented in Figs.2 and 3,where γ,φ,andVare climb angle,roll angle and velocity respectively.

    It can be seen from Fig.2 that the initial responses of the aircraft are mainly the changes of the three-axis angular velocities and attitudes when the wing is skewed.This is mainly caused by the asymmetric aerodynamic layout,resulting in a nose down pitching moment and positive rolling and yawing moments.The rolling moment sinks the front wing and the yawing moment makes the aircraft rotate to the right.As the skew angle continues to increase to 30°,the asymmetric rolling moment still exists.Then the roll angle increases rapidly and the lift direction keeps tilting.Meanwhile,the lift coefficient decreases with the increase of the skewing angle,so that the lift cannot balance the gravity of the aircraft.Therefore,the aircraft dives with acceleration and cannot reach a new steady state.

    Fig.2 Dynamic responses in the wing skewing process.

    During the changes of configurations,the attitude,altitude and speed of the OWA are all obviously changed.Therefore,a flight control law should be designed to ensure the flight safety and flying qualities of the aircraft in the dynamic process.

    3.Flight control law design based on sliding mode control method

    In the wing skewing process,the three-axis movements of the OWA are highly coupled and strongly nonlinear,and it is thus difficult to establish an accurate aerodynamic model.This leads to a certain error between the simulation results and the real flight response,meaning that there exists uncertainty in the motion model of the aircraft.In addition,the range of the motion parameters is large,which requires the flight control system to be more robust.Since the sliding mode control method can ensure the dynamic quality through the design of the sliding mode function and the approach rate,the designed controller can adapt to the parameter changes of the OWA during the wing skewing.The sliding mode controller has strong robustness and can obtain precise control effect,more suitable for the OWA with a large range of parameter variations,and with time-varying and nonlinear characteristics.Therefore,the flight control law using sliding mode control method is designed in this study.The control objective is to keep the altitude and velocity of the OWA constant during the wing skewing,so that a smooth transition of the configuration changes can be achieved.

    Fig.3 Dynamic responses during the returning to the straight wing.

    3.1.Flight control system design based on sliding mode control

    If a nonlinear system˙x=f( x,u,t)can be transformed to an affine system˙x=f( x,t)+g( x,t)u,the sliding mode control can be applied conveniently.However,the OWA's dynamic equations of motion during the wing skewing exhibit inertia and aerodynamic cross couplings,and cannot be directly transformed to the affine form.According to the characteristics of the aircraft motions,we propose that the OWA's variables of movement are divided into fast and slow variables by using the time-scale separation method.Then the complicated nonlinear equations can be simplified to an inner-loop and an outer-loop affine systems,and the sliding mode control design can be easily applied to the design of the inner-loop and outer loop sliding mode controllers.

    The sliding mode flight control system introduced in this paper mainly includes five modules:the modules of control command generation,outer-loop control law,inner-loop control law,control surface allocation,and automatic throttle.The inner and the outer loops are separated by time-scale separation method for the aircraft motions.The control structure is shown in Fig.4.

    Fig.4 Flight control structure.

    The function of each module is introduced as follows:

    (1)The control command generation module generates control commands to the outer-loop control law module and the automatic throttle module based on the states of the aircraft.In order to ensure that the aircraft's altitude and speed remain constant without sideslip,the roll angle is used to generate force to balance the side force generated by the asymmetric configuration.The output control commands are formulated as

    where ΔV, αc, βcand μcare the output command of speed,angle of attack,side slip angle and roll angle in the speed axis,respectively.HcandVcare the target altitude and speed,respectively.KΔH,K˙HandK˙yare the gains corresponding to the altitude error ΔH,the altitude varying rateand the yaw rate,tuned to be 2 rad/m,-0.5 rad s/m and-6 rad s/m according to the control performance,respectively.

    本文中我們采用哈希算法,主要的過程為:發(fā)送方采用單向哈希函數(shù)對消息進行計算,得到摘要并發(fā)送消息和摘要。接收方將接收到的消息,按同樣方式進行哈希函數(shù)計算,并將新得出的結(jié)果與發(fā)送方的原摘要結(jié)果進行比對。如結(jié)果一致,說明消息完整。在本系統(tǒng)中,摘要信息的不可變,保證了需要存證信息的完整性和真實性。將需要存證的電子數(shù)據(jù)放在區(qū)塊鏈中,避免數(shù)據(jù)被惡意篡改。

    (2)The outer-loop control model generates the input command [pc,qc,rc]Tfor the inner-loop model based on the control reference[αc,βc,μc]T.According to the time-scale separation method for the state variables,the dynamic responses of the fast variables of the inner-loop are considered to have reached steady states and are neglected in the outer-loop design.Thus,the inner and outer-loop control laws can be designed separately.The differential equation can be expressed as

    where xs=is the 3 × 3 matrix in Eq.(5),us(t)=[p,q,r]T,ys(t)is the output.

    For a system such as Eq.(6),according to the sliding mode control theory,36the sliding mode surface can be designed as

    where the values of the sliding surface control parameterscis(i=1,2,3)should makeHurwitz stable.The system can reach the sliding mode surface within a limited time and the error tends to be zero when the control law takes the following form37,38:

    In order to suppress the vibration of the system,the signfunction sgn of the control law is replaced by the following saturation function

    where ε is a small positive number.According to the control effect,the parameters are tuned as below:Rs=diag{2,2,2},Γs=diag{0.1,0.1,0.1},and ε=0.1.

    (3)The inner loop calculates the required three-axis control moments to track the command [pc,qc,rc]T.The rotational dynamics Eq.(2)can be rewritten as

    Similar to the nonlinear MIMO system of the outer loop,Mccan also be obtained by using the control law of Eq.(8).Here and∑f=also replaced by a saturation function wherea1=a2=a3=2/3,and ε=0.05.

    (4)The control allocation module calculates the required deflections of the control surfaces based on the demand of three-axis control moment Mc.Since the roll control efficiency of the OWA's aileron is low in the large skew angled configuration,the all-moving horizontal tails are used to assist roll control.27The control surface allocation can be expressed as

    where G is the matrix of the control coefficient,and u= [δa,δr,δeL,δeR],δa,δr,δeL,δeRare deflections of the aileron,rudder,left elevator and right elevator,gij(i=p,q,randj= δa,δr,δeL,δeR) are the three-axis moment control efficiencies.

    In order to avoid the full or extremely small deflection of a certain control surface,we use Δ =]to weigh the allocation.Then Eq.(11)can be rewritten as

    where GΔ=GΔ,and^u=[δa/δamax,δr/δrmax,δeL/δeLmax,δeR/δeRmax]T.The control inputucan be obtained by a pseudoinverse of Eq.(12),an expenditure of minimum control energy,39and can be expressed as

    (5)Since the drag and flight velocity of the OWA both change significantly in the wing skewing process,the throttle is used to control the speed of the aircraft.Designed by the traditional PID control method,40the automatic throttle control law is expressed as

    Fig.5 Closed-loop system responses(H=4 km and Ma=0.4).

    wherek1,k2,andk3are the gains of PID,and the parameters are tuned ask1=-6 s/m,k2=-5 m-1andk3=-1.5 s2/m.

    3.2.Closed-loop simulation and analysis

    The sliding mode control law is verified by closed-loop simulations for the wing skewing process and the process of wing returning to the straight position.The initial flight conditions areH=4 km andMa=0.4.The dynamic responses of the parameters are shown in Fig.5.

    The closed-loop simulation results for the wing skewing process are shown in Fig.5(a).It can be seen that since the lift coefficient of the aircraft gradually decreases with the increase of the skewing angle,the angle of attack gradually increases to achieve the balance of the longitudinal force.The sideslip angle is maintained as 0°,which is consistent with the target command.The aircraft rolls to the left and the component of gravity is used to balance the asymmetric side force;the aircraft eventually maintains a certain roll angle.In the entire process,the speed and altitude of the aircraft basically remain unchanged.We can see that the three-axis control surfaces deflect collaboratively to balance the asymmetric aerodynamic moments and that the throttle is decreased to balance the reduced drag.The all-moving horizontal tail participates in roll control and involves in the control allocation based on minimal control energy,and the deflections of the control surfaces and the throttle are smooth.

    The closed-loop simulation results for the process of wing returning to the straight position are shown in Fig.5(b).It can be seen that the attitude of the aircraft changes exactly in the opposite direction of the wing skewing process.The lift coefficient increases,the lateral force decreases to zero,and the angle of attack decreases.Moreover,the roll angle for balancing the side force returns to zero,while the sideslip angle remains zero.After returning to the symmetrical configuration of the straight wing,the asymmetric moment disappears and the drag increases.As a result,the lateral control surfaces(rudder,ailerons,and differential horizontal tail)decrease to zero and the throttle increases accordingly.During the whole process,the flight velocity and altitude of the aircraft basically remain unchanged.It can be seen that the sliding mode flight control law can fully guarantee the OWA flight at the required altitude and speed during the processes of wing skewing and returning to the straight position.

    Fig.6 Closed-loop system responses(H=5 km and Ma=0.5).

    To validate its robustness in the presence of the oblique wing,the sliding mode flight controller is applied to control the OWA in another flight state.Fig.6 presents the closedloop dynamic responses of the OWA during both the wing skewing and the returning to the straight wing atH=5 km andMa=0.5;the responses are controlled by the sliding mode flight controller designed for the state ofH=4 km andMa=0.4.We can see that the flight controller can still ensure that the altitude and speed in the wing skewing process are almost unchanged.This means that the sliding mode controller designed in this study has good robustness.

    4.Conclusions

    (1)The three-axis motions of the OWA are significantly coupled and highly nonlinear in the wing skewing process.The initial dynamic responses of the aircraft in the processes of wing skewing and wing returning to the straight position are mainly the changes of threeaxis angular velocities and attitudes.With the rapid increase of the roll angle,the lift direction of the aircraft continuously changes,thus being unable to balance the gravity of the aircraft.The aircraft continues to dive with acceleration in the two processes mentioned above,and cannot reach a new balance state.

    (2)The moment arm of the aileron is shortened when the wing is skewed,which leads to an insufficient rolling control effectiveness.Therefore,the all-moving horizontal tail is used for assistant roll control.According to the three-axis coupling characteristics of aileron control,the control surfaces are allocated based on minimum control energy,achieving an accurate control effect.

    (3)According to the characteristics of pitch-roll-yaw coupling,control input and state coupling,and difficulty to accurately establish the aerodynamic model,a multi-loop sliding mode flight controller is designed by the time-scale separation method for the OWA to ensure the smooth transition during the wing skewing.The simulation results show that the robust sliding mode flight controller can smoothly adjust the attitude of the OWA to balance the asymmetric aerodynamic forces generated in the processes of wing skewing and returning.In addition,the speed and altitude of the aircraft can be kept constant.

    Acknowledgement

    This work was supported by the National Natural Science Foundation of China(No.11402010).

    猜你喜歡
    信息
    訂閱信息
    中華手工(2017年2期)2017-06-06 23:00:31
    展會信息
    中外會展(2014年4期)2014-11-27 07:46:46
    信息超市
    展會信息
    展會信息
    展會信息
    展會信息
    展會信息
    信息
    健康信息
    祝您健康(1987年3期)1987-12-30 09:52:32
    av天堂久久9| 亚洲天堂av无毛| 国产午夜精品一二区理论片| 超碰97精品在线观看| 国产精品女同一区二区软件| 26uuu在线亚洲综合色| 在线天堂最新版资源| 午夜福利,免费看| 极品少妇高潮喷水抽搐| 久久99热6这里只有精品| 精品久久久精品久久久| 精品少妇久久久久久888优播| 男人狂女人下面高潮的视频| 精品久久久精品久久久| 一级毛片久久久久久久久女| 久久久国产一区二区| 深夜a级毛片| 免费观看在线日韩| 亚洲国产精品999| 亚洲av在线观看美女高潮| 韩国高清视频一区二区三区| 亚洲欧美中文字幕日韩二区| 久久午夜福利片| 人人澡人人妻人| 22中文网久久字幕| 精品人妻熟女毛片av久久网站| 久久精品国产鲁丝片午夜精品| 欧美日韩亚洲高清精品| 天堂俺去俺来也www色官网| 国产成人免费无遮挡视频| 日韩三级伦理在线观看| 亚洲欧美日韩东京热| 极品少妇高潮喷水抽搐| 在线免费观看不下载黄p国产| 亚洲欧美日韩卡通动漫| 一级片'在线观看视频| av免费观看日本| 人体艺术视频欧美日本| 久久婷婷青草| 美女视频免费永久观看网站| 日本午夜av视频| 亚洲成人av在线免费| 尾随美女入室| 婷婷色av中文字幕| 美女cb高潮喷水在线观看| 视频中文字幕在线观看| 婷婷色综合大香蕉| 高清黄色对白视频在线免费看 | 国产av一区二区精品久久| 久久人人爽人人片av| 视频中文字幕在线观看| 亚洲欧美日韩东京热| 少妇人妻精品综合一区二区| 美女内射精品一级片tv| av又黄又爽大尺度在线免费看| 在线观看av片永久免费下载| 黄色欧美视频在线观看| 春色校园在线视频观看| 国产亚洲欧美精品永久| 国产极品粉嫩免费观看在线 | av一本久久久久| av免费观看日本| 国产免费一级a男人的天堂| 欧美日本中文国产一区发布| a级一级毛片免费在线观看| 亚洲精品国产av蜜桃| 免费黄色在线免费观看| 午夜免费鲁丝| 爱豆传媒免费全集在线观看| 日韩av在线免费看完整版不卡| 不卡视频在线观看欧美| 久久久精品免费免费高清| 久久精品国产a三级三级三级| 大香蕉97超碰在线| 免费大片18禁| 日本午夜av视频| 纵有疾风起免费观看全集完整版| 国产黄片美女视频| 亚洲无线观看免费| 一级av片app| 99久久综合免费| 国产精品女同一区二区软件| 色婷婷av一区二区三区视频| 欧美日韩视频精品一区| 亚洲电影在线观看av| 在线播放无遮挡| 日本欧美视频一区| 嘟嘟电影网在线观看| 亚洲欧美一区二区三区国产| 国模一区二区三区四区视频| 成人影院久久| 日韩中文字幕视频在线看片| 26uuu在线亚洲综合色| 一级,二级,三级黄色视频| 18+在线观看网站| 久久久午夜欧美精品| 黄色一级大片看看| 国产美女午夜福利| 777米奇影视久久| 亚洲美女黄色视频免费看| 亚洲经典国产精华液单| 久久狼人影院| 在线免费观看不下载黄p国产| 人妻少妇偷人精品九色| 人人澡人人妻人| 69精品国产乱码久久久| 一区二区三区精品91| 久久久午夜欧美精品| 久久热精品热| 精品国产露脸久久av麻豆| 亚洲精华国产精华液的使用体验| 国产亚洲一区二区精品| 国产精品一区www在线观看| 一二三四中文在线观看免费高清| a级毛片免费高清观看在线播放| 国产一区二区在线观看av| 亚洲精品,欧美精品| 人妻夜夜爽99麻豆av| 人体艺术视频欧美日本| 国产精品99久久久久久久久| 日韩欧美一区视频在线观看 | 国产一区二区三区综合在线观看 | 免费av不卡在线播放| 中文资源天堂在线| 亚洲av成人精品一区久久| 日本黄色片子视频| 99热全是精品| 曰老女人黄片| 婷婷色综合www| 六月丁香七月| av卡一久久| 国精品久久久久久国模美| 视频中文字幕在线观看| 亚洲国产精品专区欧美| 91久久精品电影网| 国产熟女欧美一区二区| 三级经典国产精品| 免费看光身美女| 日日摸夜夜添夜夜添av毛片| 成年女人在线观看亚洲视频| 一级爰片在线观看| 妹子高潮喷水视频| 亚洲成人手机| 国精品久久久久久国模美| 国产黄片视频在线免费观看| 美女xxoo啪啪120秒动态图| 亚洲国产精品一区三区| 男女国产视频网站| 日韩三级伦理在线观看| 精品国产一区二区久久| 亚洲精品日韩在线中文字幕| 久久热精品热| 在线观看av片永久免费下载| 日本黄大片高清| 成人免费观看视频高清| 久热这里只有精品99| 亚洲人成网站在线观看播放| 日韩电影二区| 久久精品国产鲁丝片午夜精品| 国产精品蜜桃在线观看| 美女视频免费永久观看网站| 晚上一个人看的免费电影| 熟妇人妻不卡中文字幕| 新久久久久国产一级毛片| 精品亚洲乱码少妇综合久久| 午夜精品国产一区二区电影| 国产精品国产三级国产专区5o| 在线观看www视频免费| 夜夜骑夜夜射夜夜干| 国产亚洲最大av| 九色成人免费人妻av| 色哟哟·www| 国产精品欧美亚洲77777| 2018国产大陆天天弄谢| 国产精品三级大全| 国产精品一区二区在线不卡| 久久精品久久精品一区二区三区| 国产免费视频播放在线视频| 成人无遮挡网站| 国产精品久久久久久av不卡| 久久 成人 亚洲| 亚洲婷婷狠狠爱综合网| av卡一久久| 18禁在线无遮挡免费观看视频| 久久97久久精品| 久久久久久久精品精品| 国产色爽女视频免费观看| 超碰97精品在线观看| 狂野欧美激情性bbbbbb| 啦啦啦视频在线资源免费观看| 精品少妇久久久久久888优播| 久久青草综合色| 日本爱情动作片www.在线观看| 国产视频首页在线观看| 黄色怎么调成土黄色| 亚洲精品国产成人久久av| 国产综合精华液| 卡戴珊不雅视频在线播放| 久久综合国产亚洲精品| 日韩免费高清中文字幕av| 国语对白做爰xxxⅹ性视频网站| 国产黄色视频一区二区在线观看| 蜜桃久久精品国产亚洲av| 一级av片app| 久久久久久人妻| 亚洲国产成人一精品久久久| 欧美成人午夜免费资源| 午夜影院在线不卡| 在线观看一区二区三区激情| 在线观看一区二区三区激情| 亚洲三级黄色毛片| 免费少妇av软件| h视频一区二区三区| 国产又色又爽无遮挡免| 人人妻人人爽人人添夜夜欢视频 | 精品一区二区三卡| 女的被弄到高潮叫床怎么办| 91精品国产九色| 十八禁高潮呻吟视频 | 国产av码专区亚洲av| www.av在线官网国产| 国产精品一区二区在线观看99| 国产av精品麻豆| 久久久久久久久久久久大奶| 美女大奶头黄色视频| 日本色播在线视频| 日本av免费视频播放| 99久久精品热视频| 美女主播在线视频| 99热网站在线观看| 大片电影免费在线观看免费| 欧美国产精品一级二级三级 | 久久韩国三级中文字幕| 日韩一区二区三区影片| 晚上一个人看的免费电影| 亚洲国产欧美日韩在线播放 | 色视频www国产| 在线亚洲精品国产二区图片欧美 | 纵有疾风起免费观看全集完整版| 22中文网久久字幕| 少妇被粗大猛烈的视频| 99久久中文字幕三级久久日本| 精品久久久久久久久av| 亚洲av电影在线观看一区二区三区| 有码 亚洲区| 久久97久久精品| 99九九在线精品视频 | 国产中年淑女户外野战色| 久久99热6这里只有精品| 亚洲人成网站在线播| 蜜桃久久精品国产亚洲av| 中文在线观看免费www的网站| 国产亚洲av片在线观看秒播厂| 成人综合一区亚洲| 日韩精品免费视频一区二区三区 | 久久影院123| 黄色毛片三级朝国网站 | 亚洲精品456在线播放app| 老女人水多毛片| 日韩一区二区视频免费看| 亚洲无线观看免费| 色网站视频免费| 亚洲电影在线观看av| 亚洲国产精品一区二区三区在线| 在线看a的网站| 国产精品成人在线| 亚洲美女黄色视频免费看| 日韩大片免费观看网站| 国产精品欧美亚洲77777| 啦啦啦视频在线资源免费观看| 最近中文字幕高清免费大全6| 99九九在线精品视频 | 亚洲精品一二三| 成人亚洲精品一区在线观看| 肉色欧美久久久久久久蜜桃| 黄片无遮挡物在线观看| 少妇人妻 视频| 日韩中文字幕视频在线看片| 国产成人aa在线观看| 麻豆精品久久久久久蜜桃| 亚洲精品一区蜜桃| 秋霞伦理黄片| 国产成人精品福利久久| 亚洲欧美精品自产自拍| 国产亚洲5aaaaa淫片| 亚洲国产av新网站| 大片电影免费在线观看免费| 青青草视频在线视频观看| 大香蕉97超碰在线| 热re99久久精品国产66热6| 精品亚洲成国产av| 少妇的逼水好多| 成人亚洲精品一区在线观看| 久久热精品热| 免费观看a级毛片全部| 日韩av免费高清视频| 国产精品久久久久久精品电影小说| 久久精品夜色国产| 黑人巨大精品欧美一区二区蜜桃 | 男人爽女人下面视频在线观看| 亚洲av综合色区一区| 精品视频人人做人人爽| 人妻制服诱惑在线中文字幕| 欧美日韩一区二区视频在线观看视频在线| 久久99精品国语久久久| 亚洲精品乱久久久久久| 视频中文字幕在线观看| 亚洲四区av| 久久久久久久大尺度免费视频| 在线观看一区二区三区激情| 久久久精品免费免费高清| 看十八女毛片水多多多| 国内少妇人妻偷人精品xxx网站| 少妇人妻 视频| 日日撸夜夜添| 国产日韩一区二区三区精品不卡 | 国产精品熟女久久久久浪| 性色av一级| 亚洲第一av免费看| 亚洲经典国产精华液单| 欧美bdsm另类| 国产精品一区二区在线不卡| 久久6这里有精品| 成人二区视频| 午夜精品国产一区二区电影| 国产毛片在线视频| 街头女战士在线观看网站| 99九九线精品视频在线观看视频| 丰满饥渴人妻一区二区三| 国产探花极品一区二区| 亚洲天堂av无毛| 亚洲精品久久午夜乱码| 国产亚洲一区二区精品| 久久久精品免费免费高清| 九九在线视频观看精品| 欧美日韩视频精品一区| 国产亚洲午夜精品一区二区久久| 免费观看的影片在线观看| 亚洲欧美一区二区三区国产| 99热这里只有精品一区| 国产精品不卡视频一区二区| 三级国产精品欧美在线观看| 久久精品久久久久久噜噜老黄| 搡老乐熟女国产| 久久久久人妻精品一区果冻| 曰老女人黄片| 极品少妇高潮喷水抽搐| 国产探花极品一区二区| 国产精品国产三级专区第一集| 久久ye,这里只有精品| 一本久久精品| 性色av一级| av视频免费观看在线观看| 最近中文字幕2019免费版| 欧美精品一区二区大全| 18禁裸乳无遮挡动漫免费视频| 男男h啪啪无遮挡| 热re99久久精品国产66热6| 国产精品人妻久久久久久| 国产免费一级a男人的天堂| 97超视频在线观看视频| 一级毛片电影观看| 人妻人人澡人人爽人人| 欧美成人午夜免费资源| 欧美激情国产日韩精品一区| 亚洲精品乱码久久久久久按摩| 18禁裸乳无遮挡动漫免费视频| 人人澡人人妻人| a 毛片基地| 黄色欧美视频在线观看| 插逼视频在线观看| 一区二区av电影网| 激情五月婷婷亚洲| 人妻少妇偷人精品九色| 丝袜喷水一区| 亚洲av二区三区四区| 最新中文字幕久久久久| 久久久久久久精品精品| 亚洲精品456在线播放app| 亚洲性久久影院| 久久久亚洲精品成人影院| 又粗又硬又长又爽又黄的视频| 在现免费观看毛片| 国产乱人偷精品视频| 国产av一区二区精品久久| 男女免费视频国产| av又黄又爽大尺度在线免费看| 天堂8中文在线网| 国产69精品久久久久777片| 久久毛片免费看一区二区三区| 免费观看在线日韩| 天堂中文最新版在线下载| 少妇人妻久久综合中文| av卡一久久| 日韩中字成人| 亚洲va在线va天堂va国产| 午夜免费男女啪啪视频观看| 国产精品国产三级专区第一集| 亚洲欧美日韩卡通动漫| 又爽又黄a免费视频| 我的老师免费观看完整版| 一本—道久久a久久精品蜜桃钙片| 最近中文字幕高清免费大全6| 中文乱码字字幕精品一区二区三区| 久久久国产欧美日韩av| 老司机影院毛片| 亚洲av成人精品一二三区| 久久97久久精品| 乱人伦中国视频| 三级经典国产精品| 久久久久精品久久久久真实原创| 青青草视频在线视频观看| 国产日韩欧美视频二区| 欧美变态另类bdsm刘玥| 亚洲国产精品一区二区三区在线| 精品午夜福利在线看| 9色porny在线观看| 少妇人妻精品综合一区二区| 大又大粗又爽又黄少妇毛片口| 午夜福利网站1000一区二区三区| 亚洲av.av天堂| 人妻人人澡人人爽人人| 日日啪夜夜爽| 高清午夜精品一区二区三区| 桃花免费在线播放| 人妻人人澡人人爽人人| 十八禁高潮呻吟视频 | 久久热精品热| 激情五月婷婷亚洲| 亚洲精品成人av观看孕妇| 欧美高清成人免费视频www| 嫩草影院新地址| 高清欧美精品videossex| 免费黄频网站在线观看国产| 婷婷色麻豆天堂久久| 亚洲高清免费不卡视频| 国产高清三级在线| 亚洲精品日韩在线中文字幕| 美女福利国产在线| av在线老鸭窝| 最新的欧美精品一区二区| 91精品伊人久久大香线蕉| 另类精品久久| 伦理电影大哥的女人| 日韩不卡一区二区三区视频在线| 欧美性感艳星| 新久久久久国产一级毛片| 汤姆久久久久久久影院中文字幕| 免费观看无遮挡的男女| 精品视频人人做人人爽| 精品久久久久久久久亚洲| 2021少妇久久久久久久久久久| 在线精品无人区一区二区三| av免费在线看不卡| 三级国产精品片| 人体艺术视频欧美日本| 亚洲色图综合在线观看| 青春草视频在线免费观看| 成人国产麻豆网| 丝袜脚勾引网站| 国产白丝娇喘喷水9色精品| 有码 亚洲区| 乱人伦中国视频| 天天躁夜夜躁狠狠久久av| 校园人妻丝袜中文字幕| 成年女人在线观看亚洲视频| 美女主播在线视频| 人妻少妇偷人精品九色| 黄色日韩在线| 一级毛片电影观看| av在线老鸭窝| 黄片无遮挡物在线观看| 久久99一区二区三区| 国产亚洲最大av| 观看免费一级毛片| 久久99蜜桃精品久久| 国产精品熟女久久久久浪| 国产成人freesex在线| 一级片'在线观看视频| 亚洲精品乱码久久久v下载方式| 日韩人妻高清精品专区| 毛片一级片免费看久久久久| 国产午夜精品一二区理论片| 日本vs欧美在线观看视频 | 国产精品.久久久| 久久久久人妻精品一区果冻| 亚洲av成人精品一区久久| 久久精品久久久久久噜噜老黄| 国语对白做爰xxxⅹ性视频网站| 日韩成人伦理影院| 自拍欧美九色日韩亚洲蝌蚪91 | 天天躁夜夜躁狠狠久久av| 免费观看性生交大片5| 中文字幕久久专区| 少妇被粗大猛烈的视频| 韩国av在线不卡| a级毛片在线看网站| 插阴视频在线观看视频| 最近中文字幕2019免费版| 久久久国产一区二区| 亚洲精品一区蜜桃| 欧美亚洲 丝袜 人妻 在线| 亚洲四区av| 成人免费观看视频高清| 99九九线精品视频在线观看视频| 中文字幕久久专区| 亚洲av成人精品一二三区| 男人狂女人下面高潮的视频| 天天操日日干夜夜撸| 欧美日韩国产mv在线观看视频| 国产成人一区二区在线| 国产精品三级大全| 欧美成人午夜免费资源| 最新的欧美精品一区二区| 国产老妇伦熟女老妇高清| 女性生殖器流出的白浆| 一级毛片久久久久久久久女| 成人午夜精彩视频在线观看| 女人精品久久久久毛片| 国产精品一区www在线观看| 一级毛片黄色毛片免费观看视频| 亚洲色图综合在线观看| a级片在线免费高清观看视频| av又黄又爽大尺度在线免费看| 亚洲美女搞黄在线观看| 亚洲怡红院男人天堂| 欧美成人精品欧美一级黄| 亚洲av国产av综合av卡| 下体分泌物呈黄色| 国产精品久久久久久久电影| 18禁在线无遮挡免费观看视频| 欧美+日韩+精品| 国产高清国产精品国产三级| 日韩视频在线欧美| 九九久久精品国产亚洲av麻豆| 妹子高潮喷水视频| 亚洲精品视频女| 国产精品久久久久久精品电影小说| 久久国产精品大桥未久av | 精品一区二区三卡| 97超视频在线观看视频| 熟妇人妻不卡中文字幕| 哪个播放器可以免费观看大片| 免费在线观看成人毛片| 男女啪啪激烈高潮av片| 大香蕉97超碰在线| 国产精品久久久久久久电影| 亚洲欧美日韩卡通动漫| 内射极品少妇av片p| 国产亚洲欧美精品永久| tube8黄色片| 中国三级夫妇交换| 十八禁高潮呻吟视频 | 欧美xxxx性猛交bbbb| 欧美xxⅹ黑人| 亚洲综合精品二区| 亚洲av日韩在线播放| 女性被躁到高潮视频| 夜夜看夜夜爽夜夜摸| 精品少妇内射三级| 国产精品.久久久| 久久青草综合色| 伦精品一区二区三区| 国产精品一区二区在线不卡| 国产av一区二区精品久久| 又粗又硬又长又爽又黄的视频| 丝袜喷水一区| 极品人妻少妇av视频| 亚洲图色成人| 精品久久久久久久久亚洲| 嘟嘟电影网在线观看| 一区在线观看完整版| 18禁动态无遮挡网站| 国产精品成人在线| 韩国高清视频一区二区三区| 男人舔奶头视频| 97精品久久久久久久久久精品| 我的老师免费观看完整版| 热99国产精品久久久久久7| 在线天堂最新版资源| 亚洲婷婷狠狠爱综合网| 亚洲中文av在线| 99九九在线精品视频 | 一级毛片久久久久久久久女| 少妇的逼水好多| av天堂久久9| 97超碰精品成人国产| 国产日韩一区二区三区精品不卡 | 丰满少妇做爰视频| 欧美日韩一区二区视频在线观看视频在线| 能在线免费看毛片的网站| 国产精品久久久久久久电影| 精品久久久精品久久久| 色网站视频免费| 插阴视频在线观看视频| 色网站视频免费| 欧美 亚洲 国产 日韩一| 日本免费在线观看一区| 最近手机中文字幕大全| 永久网站在线| 老司机影院成人| 交换朋友夫妻互换小说| 国产精品一区www在线观看| 亚洲精品乱码久久久久久按摩| 97精品久久久久久久久久精品| 91午夜精品亚洲一区二区三区| 欧美xxxx性猛交bbbb| 亚洲精品,欧美精品| 美女视频免费永久观看网站| av免费观看日本| 自线自在国产av| 亚洲伊人久久精品综合| 伊人久久国产一区二区| 桃花免费在线播放| 黄色欧美视频在线观看|