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

    高超聲速飛行器線性變參數(shù)一體化式控制律設(shè)計

    2022-12-01 08:04:36錢云霄
    上海交通大學(xué)學(xué)報 2022年11期
    關(guān)鍵詞:參考系超聲速航跡

    楊 庶, 錢云霄, 楊 婷,2

    (1. 西北工業(yè)大學(xué) 自動化學(xué)院,西安 710072; 2. 北京空天技術(shù)研究所,北京 100074)

    高超聲速飛行器可進行跨空域、跨速域飛行,其復(fù)雜的飛行環(huán)境和時變的飛行動態(tài)給控制系統(tǒng)的設(shè)計帶來嚴(yán)峻挑戰(zhàn).針對上述問題,目前常采用動態(tài)逆控制、反演控制、滑??刂啤⑶袚Q控制、自適應(yīng)控制、線性變參數(shù)(LPV)控制等方法進行控制器設(shè)計.文獻[1]針對制導(dǎo)控制一體化問題,提出一種基于魯棒動態(tài)逆的動態(tài)面一體化算法,通過引入動態(tài)面控制,避免反演控制的“計算膨脹”問題.文獻[2]設(shè)計了一種無需在線求導(dǎo)的高超聲速飛行器反步法控制律方法,用于實現(xiàn)高超聲速飛行器縱向通道的航跡傾斜角跟蹤控制.模型的強非線性和強耦合性特點使得動態(tài)逆、反步法等控制器設(shè)計過程復(fù)雜,增加了實際應(yīng)用難度[1,3].對于高超聲速飛行器縱向控制通道的速度和高度控制問題,文獻[4]設(shè)計干擾補償器和狀態(tài)反饋控制律,在實現(xiàn)指令跟蹤的同時,提高飛行器閉環(huán)系統(tǒng)相對于外界擾動的魯棒性.文獻[5]給出一種基于模糊控制理論的高超聲速飛行器縱向最優(yōu)控制律,能夠在飛行器未建模動力學(xué)特性的影響下,實現(xiàn)對飛行器高度和速度指令的跟蹤.文獻[6]采用有限時間滑??刂评碚?,設(shè)計高超聲速飛行器在俯沖段的協(xié)同制導(dǎo)律,實現(xiàn)多枚高超聲速飛行器攻擊同一固定目標(biāo)的控制任務(wù).文獻[7-8]分別采用滑模和自適應(yīng)控制理論實現(xiàn)高超聲速飛行器俯沖段的控制.文獻[9]針對吸氣式高超聲速飛行器的縱向飛行控制,建立多胞剛性LPV模型,進而基于區(qū)域極點配置實現(xiàn)系統(tǒng)的H∞/H2多目標(biāo)魯棒控制.文獻[10]將切換控制思想與LPV控制結(jié)合,通過減小子系統(tǒng)的參數(shù)變化范圍,增強系統(tǒng)的魯棒性.

    上述控制方法中,LPV控制利用參數(shù)依賴的思想,既可以有效適應(yīng)高超聲速飛行器飛行包線大、參數(shù)快時變的特點,又可以基于增益調(diào)度技術(shù)保證閉環(huán)系統(tǒng)的穩(wěn)定性和控制器的實用性,在高超聲速飛行器控制中得到廣泛應(yīng)用.需要指出的是,現(xiàn)有的飛行器LPV控制研究一般僅考慮飛行器縱向或橫航向動力學(xué)特性[11-12],忽略縱向和橫航向通道之間的耦合.基于飛行器完整模型的LPV控制研究通常引入模型假設(shè)(如速度傾側(cè)角為0),降低控制律的設(shè)計難度.

    本文基于高度-水平航跡控制概念,將高超聲速飛行器三維航跡分解為豎直面航跡和水平面航跡,通過高度、速度、側(cè)滑角、偏航角控制,實現(xiàn)飛行器的三維航跡控制.基于高度-水平航跡控制概念,采用具有L2誘導(dǎo)范數(shù)控制性能的LPV控制方法,設(shè)計高超聲速飛行器的一體化式LPV控制律.在地心地固(ECEF)參考系下建立高超聲速飛行器的數(shù)學(xué)模型,通過數(shù)值仿真檢驗一體化式LPV控制律在測量噪聲影響下的控制性能.仿真結(jié)果表明,該控制律能夠跟蹤給定控制指令,并具有較好的噪聲抑制性能.

    1 高超聲速飛行器數(shù)學(xué)模型

    基于WGS84地球模型參數(shù),在ECEF參考系(見圖1)下建立飛行器的非線性數(shù)學(xué)模型,在模型中考慮地球自轉(zhuǎn)、地球扁率、地球引力二階簡諧效應(yīng)的影響.圖中:l和μ為大地經(jīng)度和大地緯度;ωE為地球自轉(zhuǎn)速度;lE為地球自轉(zhuǎn)角;ECI代表地心慣性參考系;B代表機體參考系.為便于控制律設(shè)計,基于平板地球假設(shè),建立飛行器的簡化非線性數(shù)學(xué)模型,進而得到高超聲速飛行器的LPV模型.在兩組數(shù)學(xué)模型中均假設(shè)飛行器為剛體,并忽略飛行器的質(zhì)量變化.

    圖1 ECEF參考系與機體參考系示意圖Fig.1 ECEF reference frame and body reference frame

    1.1 ECEF參考系下的運動方程

    由于地球自轉(zhuǎn)的影響,圖1所示的ECEF參考系是相對于ECI參考系的非慣性參考系,其角速率向量為

    (1)

    (2)

    在ECEF參考系中,高超聲速飛行器的平動運動受到地球自轉(zhuǎn)角速度、Coriolis效應(yīng)、地球扁率等因素影響,根據(jù)牛頓第二定律,在機體參考下可將高超聲速飛行器的平動動力學(xué)方程表達為

    (3)

    ge=

    (4)

    式中:GM為地心引力常數(shù);J2為地球引力二階簡諧效應(yīng)常數(shù);a為地球半長軸;μc為地心緯度;rk(k=1,2,3)為飛行器位移向量r的第k個分量;‖r‖為飛行器相對于地心的距離.GM,J2,a均根據(jù)WGS84地球模型參數(shù)確定,式(3)和(4)的詳細(xì)推導(dǎo)過程參見文獻[13].

    根據(jù)文獻[14]中的氣動參數(shù)進行計算,有

    Fb=

    (5)

    式中:α為迎角;T為發(fā)動機推力,可通過燃油流量比δt進行控制;D,L和Y為飛行器的阻力、升力和側(cè)力,可分別通過副翼舵偏角δa、升降舵舵偏角δe和方向舵舵偏角δr進行控制.

    高超聲速飛行器的轉(zhuǎn)動動力學(xué)方程為

    (6)

    ωb=ωbn+Rbnωne+Rbeωe

    (7)

    (8)

    高超聲速飛行器的位移和四元數(shù)姿態(tài)運動學(xué)方程為

    (9)

    1.2 平板地球假設(shè)下的運動方程和LPV模型

    為便于控制律設(shè)計,可忽略地球自轉(zhuǎn)和引力二階簡諧效應(yīng),根據(jù)式(3)、式(6)和式(9)進行簡化,在平板地球假設(shè)下建立高超聲速飛行器的簡化運動方程[13],即

    (10)

    根據(jù)基于線性化的LPV模型建模方法[16],可將式(10)轉(zhuǎn)化為如下高超聲速飛行器LPV模型,即

    (11)

    2 一體化式LPV控制律

    針對高超聲速飛行器的三維航跡控制需求,基于高度-水平航跡控制概念設(shè)計高超聲速飛行器一體化式LPV控制律,進而采用線性矩陣不等式(LMI)方法求解該LPV輸出反饋控制律的參數(shù).

    2.1 高度-水平航跡控制概念

    根據(jù)笛卡爾參考系中基向量相互正交的性質(zhì),可根據(jù)北東地參考系,將高超聲速飛行器的三維航跡控制問題分解為水平面和豎直面航跡控制問題(見圖2).豎直面航跡控制可通過高度控制實現(xiàn),水平面航跡控制問題可通過速度和偏航角控制實現(xiàn).由于高超聲速飛行器通常要求實現(xiàn)協(xié)調(diào)轉(zhuǎn)彎性能,即水平面航跡控制時側(cè)滑角為0°,所以側(cè)滑角控制應(yīng)作為水平面航跡控制目標(biāo)之一.

    根據(jù)上述討論,控制指令可選:速度V、側(cè)滑角β、高度h、偏航角ψ.由于這些控制指令既包含常規(guī)飛行控制的內(nèi)回路姿態(tài)控制指令,也包含外回路控制指令,所以通常需要進行內(nèi)外回路控制律的分離設(shè)計與綜合,進而實現(xiàn)飛行器三維航跡控制.由于內(nèi)外回路控制律分離設(shè)計方法有可能導(dǎo)致控制律之間的耦合問題,所以本文采用LPV輸出反饋控制理論,設(shè)計高超聲速飛行器一體化式LPV飛行控制律,使用LMI方法求解具有L2誘導(dǎo)范數(shù)性能的LPV魯棒控制律參數(shù).

    圖2 高度-水平航跡控制概念示意圖Fig.2 Altitude-horizontal trajectory control concept

    2.2 LPV控制律設(shè)計

    高超聲速飛行器一體化式LPV控制律如圖3所示,通過參考模型和L2誘導(dǎo)范數(shù)意義下的模型匹配方法[17],實現(xiàn)速度、側(cè)滑角、高度、偏航角通道的模型跟蹤控制性能.圖中:下標(biāo)c代表控制指令;向量n為反饋噪聲.各控制通道的參考模型均為阻尼比0.7的二階系統(tǒng),則有傳遞函數(shù):

    (12)

    式中:s為傳遞函數(shù)的復(fù)頻域自變量.

    誤差權(quán)值函數(shù)設(shè)計為

    (13)

    式中:WV(s)和Wh(s)的低頻增益為0.01,設(shè)計目標(biāo)為使速度和高度的穩(wěn)態(tài)控制誤差不超過1%;Wβ(s)和Wψ(s)的低頻增益為0.001,設(shè)計目標(biāo)為使側(cè)滑角和偏航角的穩(wěn)態(tài)控制誤差不超過0.1%.

    控制權(quán)值函數(shù)設(shè)計為

    (14)

    其設(shè)計目標(biāo)為在單位指令信號的作用下,高超聲速飛行器燃油流量的變化量不超過20%,副翼、升降舵、方向舵的變化量不超過10°.式(13)和式(14)中各權(quán)值函數(shù)的輸出分別對應(yīng)于圖3中的eV、eβ、eh、eψ和eu.

    反饋信號設(shè)計為

    xf=

    (15)

    各反饋信號均受測量噪聲影響.

    圖3 LPV控制律與系統(tǒng)連接結(jié)構(gòu)Fig.3 System interconnections for LPV control law

    噪聲權(quán)值函數(shù)設(shè)計為

    Wn=

    (16)

    式中:I3和I6分別為三維單位矩陣和六維單位矩陣.式(16)代表速度、位移、角度和角速率反饋信號的噪聲幅值分別不超過0.01 m/s、0.01 m、0.01° 和 1.75×10-4rad/s.

    根據(jù)式(11)~(16),用于控制律設(shè)計的LPV模型可表達為

    (17)

    式中:xL為狀態(tài)變量,包括式(11)~(13)的狀態(tài)變量;d為廣義輸入,包括控制指令和測量噪聲;e為廣義誤差,由式(13)的輸出信號給定;y為輸出信號,由控制指令和反饋信號給定.對于式(17)中的線性系統(tǒng),A為系統(tǒng)狀態(tài)矩陣,B1為廣義輸入的輸入矩陣,B2為控制輸入的輸入矩陣,C1為廣義誤差的輸出矩陣,C2為輸出信號的輸出矩陣,D11為從廣義輸入到廣義誤差的前饋矩陣,D12為從控制輸入到廣義誤差的前饋矩陣,D21為從廣義輸入到輸出信號的前饋矩陣.

    控制律在馬赫數(shù)包線內(nèi)進行設(shè)計,所采用的馬赫數(shù)包線為4≤Ma≤6.為保證LMI的約束條件為有限維,采用如下的包線離散網(wǎng)格進行控制律參數(shù)求解:

    P=

    {4, 4.25, 4.5, 4.75, 5, 5.25, 5.5, 5.75, 6}

    (18)

    LPV控制器參數(shù)通過求解如下的LMI進行確定[18],即

    (19)

    (20)

    式(19)中的LMI使閉環(huán)系統(tǒng)具有D-穩(wěn)定性[18-19],并且保證閉環(huán)系統(tǒng)滿足如下的L2誘導(dǎo)范數(shù)性能,即輸入-輸出L2穩(wěn)定性:

    (21)

    LPV控制律可表達為

    (22)

    控制律的參數(shù)可通過下式計算:

    (23)

    式中:矩陣N為中間變量.式(23)的具體推導(dǎo)過程參見文獻[18].式(19)中的待求矩陣采用如下的參數(shù)化格式,即

    (24)

    式中:0, 1, 2代表階次;ρ定義為ρ=(Ma-4)/2.矩陣R(ρ)表達為馬赫數(shù)的二次方程形式,使控制律能夠?qū)Ω叱曀亠w行器在不同飛行速度下的運動特性非線性變化進行有效控制.其余待求矩陣均設(shè)為常數(shù)矩陣,降低控制律求解和實現(xiàn)所需計算量.

    3 仿真結(jié)果

    基于高超聲速飛行器平板地球假設(shè)下的運動方程,分析飛行器開環(huán)系統(tǒng)和閉環(huán)系統(tǒng)的穩(wěn)定性.采用ECEF參考系下的高超聲速飛行器運動方程,通過豎直面機動和水平面機動的數(shù)值仿真,檢驗一體化式LPV控制律的控制性能和魯棒性.

    3.1 穩(wěn)定性分析

    根據(jù)平板地球假設(shè)下的運動方程,高超聲速飛行器在馬赫數(shù)包線4≤Ma≤6中的極點分布如圖4所示.圖中:zb為飛行器的極點;ωn和ξ分別為飛行器極點的自然頻率和阻尼比.由圖可見,飛行器是開環(huán)不穩(wěn)定系統(tǒng),不穩(wěn)定極點的最高自然頻率接近于10 rad/s,說明該飛行器的開環(huán)不穩(wěn)定性較高,而本文設(shè)計的一體化式LPV控制律能夠保證閉環(huán)系統(tǒng)為穩(wěn)定系統(tǒng),并且所有極點在復(fù)平面上均分布在阻尼比大于0.5、自然頻率小于30 rad/s的扇形區(qū)域內(nèi),實現(xiàn)式(19)給定的閉環(huán)系統(tǒng)D-穩(wěn)定性要求.同時,閉環(huán)系統(tǒng)的大多數(shù)極點的阻尼比大于0.7,說明一體化式LPV控制律能夠為高超聲速飛行器閉環(huán)系統(tǒng)提供良好的阻尼特性.

    圖4 高超聲速飛行器開環(huán)和閉環(huán)系統(tǒng)極點分布Fig.4 Open-loop and closed-loop poles of a hypersonic vehicle

    3.2 躍升加速機動

    采用ECEF參考系下的高超聲速飛行器運動方程進行數(shù)值仿真,使飛行器在豎直面內(nèi)進行高度躍升和水平加速機動,檢驗基于飛行器簡化模型設(shè)計的LPV控制律在地球自轉(zhuǎn)、地球扁率、地球引力二階簡諧效應(yīng)等擾動下的控制性能.仿真中,在反饋信號中增加測量噪聲,檢驗該控制律相對于噪聲的魯棒性.速度、位移、角度和角速率反饋信號的噪聲幅值分別設(shè)置為0.005 m/s、0.005 m、0.005°和8.726×10-5rad/s.

    如圖5所示,一體化式LPV控制律能夠使高超聲速飛行器跟蹤給定的高度和速度指令,機動過程中側(cè)滑角和偏航角沒有出現(xiàn)明顯的控制誤差.同時,測量噪聲對飛行器閉環(huán)系統(tǒng)控制性能的影響很小,沒有出現(xiàn)控制性能惡化的現(xiàn)象.圖中:t為時間.從圖6和圖7所示的高超聲速飛行器在ECEF參考系下的狀態(tài)變量響應(yīng)可以看出,高超聲速飛行器具有閉環(huán)穩(wěn)定性,說明該控制律在地球自轉(zhuǎn)、地球扁率、地球引力二階簡諧效應(yīng)等擾動和測量噪聲的影響下仍具有良好的控制性能和魯棒性.

    圖5 控制指令與系統(tǒng)響應(yīng)(躍升加速機動)Fig.5 Control commands and system responses (ascent-acceleration maneuver)

    圖6 ECEF參考系下高超聲速飛行器運動方程的狀態(tài)響應(yīng)(躍升加速機動)Fig.6 State responses of hypersonic vehicle equations of motion in ECEF reference frame (ascent-acceleration maneuver)

    圖7 ECEF參考系下高超聲速飛行器的姿態(tài)四元數(shù)響應(yīng)(躍升加速機動)Fig.7 Quaternion responses of hypersonic vehicle equations of motion in ECEF reference frame (ascent-acceleration maneuver)

    從圖8所示的控制輸入可以看出,在標(biāo)稱條件下高超聲速飛行器完成躍升加速機動所需的控制輸入量較小.但在測量噪聲的影響下,LPV控制律輸出的控制輸入中包含較為明顯的噪聲,不利于實際應(yīng)用.為消除該噪聲影響,在LPV控制律的反饋信號中增加二階低通濾波器:

    (25)

    式中:濾波器自然頻率ωf=40 rad/s;阻尼比ξf=1.

    圖8 控制輸入(躍升加速機動)Fig.8 Control inputs (ascent-acceleration maneuver)

    從圖5~7的仿真曲線中可以看出,對反饋信號進行濾波后,高超聲速飛行器閉環(huán)系統(tǒng)在測量噪聲的影響下仍然具有良好的響應(yīng)特性,并且LPV控制律生成的控制輸入得到明顯的降噪和平滑化.

    3.3 水平轉(zhuǎn)彎機動

    采用ECEF參考系下的高超聲速飛行器運動方程進行數(shù)值仿真,使飛行器在水平面內(nèi)進行水平轉(zhuǎn)彎機動.仿真中,在反饋信號中增加測量噪聲,噪聲幅值設(shè)置與3.2節(jié)相同.

    從圖9可以看出,一體化式LPV控制律能夠使高超聲速飛行器跟蹤給定的偏航角指令,使飛行器完成45°的水平轉(zhuǎn)彎機動,機動過程中側(cè)滑角的控制誤差小于0.2°,說明飛行器實現(xiàn)了協(xié)調(diào)轉(zhuǎn)彎控制性能.同時,測量噪聲對飛行器閉環(huán)系統(tǒng)的控制性能影響很小,沒有出現(xiàn)控制性能惡化的現(xiàn)象.

    圖9 控制指令與系統(tǒng)響應(yīng)(水平轉(zhuǎn)彎機動)Fig.9 Control commands and system responses (level turn maneuver)

    從圖10和圖11所示的高超聲速飛行器在ECEF參考系下狀態(tài)變量響應(yīng)可以看出,高超聲速飛行器具有閉環(huán)穩(wěn)定性,說明該控制律在地球自轉(zhuǎn)、地球扁率、地球引力二階簡諧效應(yīng)等擾動和測量噪聲的影響下仍具有良好的控制性能和魯棒性.從圖12所示的控制輸入可以看出, 高超聲速飛行器完成水平轉(zhuǎn)彎機動所需的控制輸入量較小.

    圖10 ECEF參考系下高超聲速飛行器運動方程的狀態(tài)響應(yīng)(水平轉(zhuǎn)彎機動)Fig.10 State responses of hypersonic vehicle equations of motion in ECEF reference frame (level turn maneuver)

    圖11 ECEF參考系下高超聲速飛行器的姿態(tài)四元數(shù)響應(yīng)(水平轉(zhuǎn)彎機動)Fig.11 Quaternion responses of hypersonic vehicle equations of motion in ECEF reference frame (level turn maneuver)

    圖12 控制輸入(水平轉(zhuǎn)彎機動)Fig.12 Control inputs (level turn maneuver)

    從圖9~12中的仿真曲線可知,在應(yīng)用式(25)中的二階低通濾波器對反饋信號進行濾波后,測量噪聲對控制輸入的影響得到明顯的抑制,同時LPV控制律仍能夠保持良好的控制性能.

    4 結(jié)論

    針對高超聲速飛行器三維航跡控制問題,基于高度-水平航跡控制概念,設(shè)計了一體化式LPV控制律.該控制律不區(qū)分常規(guī)飛行控制的內(nèi)外控制回路,在L2誘導(dǎo)范數(shù)意義下具有良好的控制精度和魯棒性,具體結(jié)論如下:

    (1) 高度-水平航跡控制概念基于笛卡爾參考系基向量相互正交的性質(zhì),將三維航跡控制問題轉(zhuǎn)化為水平面和豎直面航跡控制問題,等價的控制指令為速度、側(cè)滑角、高度和偏航角.

    (2) 一體化式LPV控制律在馬赫數(shù)包線內(nèi)進行設(shè)計,在4≤Ma≤6包線內(nèi)使高超聲速飛行器閉環(huán)系統(tǒng)具有L2誘導(dǎo)范數(shù)控制性能,在保證控制指令跟蹤精度的同時,實現(xiàn)相對于噪聲的魯棒性.

    (3) 基于極點配置和LMI方法設(shè)計的一體化式LPV控制律使高超聲速飛行器閉環(huán)系統(tǒng)的極點具有阻尼比不低于0.5、自然頻率不高于30 rad/s的性能.

    (4) 采用ECEF參考系下建立的高超聲速飛行器運動方程進行數(shù)值仿真,一體化式LPV控制律能夠在地球自轉(zhuǎn)、地球扁率、地球引力二階簡諧效應(yīng)等擾動下保持良好的控制性能,并且相對于反饋信號中的測量噪聲具有良好的魯棒性.

    猜你喜歡
    參考系超聲速航跡
    為何點不能作為參考系
    ——兼談參考系與坐標(biāo)系的關(guān)聯(lián)關(guān)系
    高超聲速出版工程
    高超聲速飛行器
    夢的航跡
    青年歌聲(2019年12期)2019-12-17 06:32:32
    超聲速旅行
    自適應(yīng)引導(dǎo)長度的無人機航跡跟蹤方法
    視覺導(dǎo)航下基于H2/H∞的航跡跟蹤
    靈活選擇參考系簡化多物體相對運動問題お
    平面轉(zhuǎn)動參考系的教學(xué)改進
    巧選參考系簡化解題過程
    欧美3d第一页| 热re99久久精品国产66热6| av片东京热男人的天堂| 亚洲第一av免费看| 国产日韩欧美视频二区| 国产乱来视频区| 在线免费观看不下载黄p国产| 色婷婷久久久亚洲欧美| 爱豆传媒免费全集在线观看| 国产免费一级a男人的天堂| 亚洲丝袜综合中文字幕| 飞空精品影院首页| 日韩欧美精品免费久久| 男女边吃奶边做爰视频| 日本与韩国留学比较| 久久久国产精品麻豆| 又黄又爽又刺激的免费视频.| 欧美日韩亚洲高清精品| 亚洲,一卡二卡三卡| 日韩大片免费观看网站| 这个男人来自地球电影免费观看 | 成人毛片a级毛片在线播放| 国产精品不卡视频一区二区| 国产日韩欧美亚洲二区| 国产日韩欧美亚洲二区| 国产老妇伦熟女老妇高清| 久久热在线av| 日韩制服骚丝袜av| 美女大奶头黄色视频| 免费看光身美女| 国产日韩欧美在线精品| 激情视频va一区二区三区| av国产久精品久网站免费入址| 国产精品免费大片| 免费观看a级毛片全部| 欧美xxⅹ黑人| 91精品三级在线观看| 乱码一卡2卡4卡精品| 建设人人有责人人尽责人人享有的| 少妇猛男粗大的猛烈进出视频| 国产精品欧美亚洲77777| 夫妻午夜视频| 成人漫画全彩无遮挡| 久久97久久精品| 好男人视频免费观看在线| 精品一区在线观看国产| 最新的欧美精品一区二区| 国产av码专区亚洲av| 狠狠婷婷综合久久久久久88av| 日本与韩国留学比较| 久久精品国产自在天天线| 亚洲伊人久久精品综合| xxx大片免费视频| 日本wwww免费看| 美女主播在线视频| 国产成人91sexporn| 中国三级夫妇交换| 777米奇影视久久| 丰满乱子伦码专区| 在线观看三级黄色| 满18在线观看网站| 一区二区日韩欧美中文字幕 | 男人添女人高潮全过程视频| 久久婷婷青草| 成年动漫av网址| 丰满饥渴人妻一区二区三| 最后的刺客免费高清国语| 视频中文字幕在线观看| 夫妻午夜视频| 伦理电影免费视频| 亚洲高清免费不卡视频| 人妻少妇偷人精品九色| 香蕉精品网在线| 在线观看一区二区三区激情| videos熟女内射| 女人被躁到高潮嗷嗷叫费观| 高清在线视频一区二区三区| 丝瓜视频免费看黄片| 性色av一级| 国产爽快片一区二区三区| 日韩欧美一区视频在线观看| 久久精品国产亚洲av天美| 日本欧美国产在线视频| av电影中文网址| 一边亲一边摸免费视频| 搡女人真爽免费视频火全软件| 久久精品久久久久久噜噜老黄| 国产男女内射视频| 丰满乱子伦码专区| 国产av国产精品国产| 日本vs欧美在线观看视频| 18在线观看网站| 久久免费观看电影| 日韩精品免费视频一区二区三区 | 最近2019中文字幕mv第一页| 国产精品久久久久久av不卡| 在线观看美女被高潮喷水网站| 国产亚洲欧美精品永久| 免费久久久久久久精品成人欧美视频 | 久久久精品区二区三区| 国产亚洲精品久久久com| 国产精品久久久久久av不卡| 人人澡人人妻人| 激情五月婷婷亚洲| xxx大片免费视频| 久久久久久伊人网av| 看非洲黑人一级黄片| 日本与韩国留学比较| 母亲3免费完整高清在线观看 | 三上悠亚av全集在线观看| 黄网站色视频无遮挡免费观看| 丝袜脚勾引网站| 亚洲婷婷狠狠爱综合网| 五月天丁香电影| 韩国av在线不卡| 国产精品久久久久久精品电影小说| 国精品久久久久久国模美| 久久人人爽av亚洲精品天堂| 亚洲欧洲日产国产| 9191精品国产免费久久| 在线观看www视频免费| 久久久国产欧美日韩av| 免费看不卡的av| 国产日韩欧美亚洲二区| 久久精品国产鲁丝片午夜精品| 久久人人97超碰香蕉20202| 国产免费又黄又爽又色| 九草在线视频观看| 91午夜精品亚洲一区二区三区| 最新中文字幕久久久久| 69精品国产乱码久久久| 欧美精品一区二区大全| 91国产中文字幕| 免费不卡的大黄色大毛片视频在线观看| 9色porny在线观看| 中文乱码字字幕精品一区二区三区| 青春草国产在线视频| 最近2019中文字幕mv第一页| 亚洲,一卡二卡三卡| 亚洲精品久久午夜乱码| 啦啦啦中文免费视频观看日本| 日韩免费高清中文字幕av| 国产成人aa在线观看| 亚洲国产欧美日韩在线播放| 久久精品国产综合久久久 | 亚洲欧美色中文字幕在线| 男女无遮挡免费网站观看| 成年动漫av网址| 久久亚洲国产成人精品v| 国产av精品麻豆| 丝袜人妻中文字幕| 男人操女人黄网站| 日本色播在线视频| 亚洲国产欧美在线一区| 一级毛片电影观看| 亚洲欧洲日产国产| 亚洲 欧美一区二区三区| 亚洲精品日本国产第一区| 激情视频va一区二区三区| 自线自在国产av| 三级国产精品片| 亚洲欧美日韩另类电影网站| 秋霞在线观看毛片| av视频免费观看在线观看| 黑丝袜美女国产一区| 日本午夜av视频| 欧美精品av麻豆av| 大片电影免费在线观看免费| 久久久久视频综合| 午夜激情久久久久久久| 免费久久久久久久精品成人欧美视频 | 精品人妻偷拍中文字幕| freevideosex欧美| 制服诱惑二区| 搡女人真爽免费视频火全软件| 九色成人免费人妻av| 美国免费a级毛片| av播播在线观看一区| 人妻一区二区av| 久久av网站| 成年人免费黄色播放视频| 国产一区二区激情短视频 | 校园人妻丝袜中文字幕| 国精品久久久久久国模美| videossex国产| 欧美人与性动交α欧美软件 | av线在线观看网站| 飞空精品影院首页| 9热在线视频观看99| 下体分泌物呈黄色| 一级黄片播放器| 黑人巨大精品欧美一区二区蜜桃 | 精品一区在线观看国产| 午夜激情av网站| 日日啪夜夜爽| av黄色大香蕉| 久久狼人影院| 日韩欧美一区视频在线观看| av有码第一页| 成人黄色视频免费在线看| 亚洲内射少妇av| 久久国内精品自在自线图片| 久久久精品区二区三区| 欧美日韩亚洲高清精品| 国产淫语在线视频| 久久久久久久久久久免费av| 两性夫妻黄色片 | 99九九在线精品视频| 好男人视频免费观看在线| 日本爱情动作片www.在线观看| 亚洲成人手机| 男人爽女人下面视频在线观看| 伊人久久国产一区二区| 80岁老熟妇乱子伦牲交| 国产精品一区二区在线观看99| videosex国产| 18在线观看网站| 亚洲经典国产精华液单| 婷婷色综合大香蕉| www.熟女人妻精品国产 | 亚洲国产精品国产精品| 亚洲第一av免费看| 人人妻人人爽人人添夜夜欢视频| 在线天堂最新版资源| 午夜日本视频在线| 国产一区二区在线观看av| 亚洲综合色网址| 日本av手机在线免费观看| 一级毛片我不卡| 一级黄片播放器| 国产色爽女视频免费观看| 国产精品人妻久久久久久| 亚洲中文av在线| 男女啪啪激烈高潮av片| 久久人人爽av亚洲精品天堂| 国产成人欧美| 成年人免费黄色播放视频| 观看av在线不卡| 亚洲婷婷狠狠爱综合网| 乱人伦中国视频| 国产有黄有色有爽视频| 日韩三级伦理在线观看| 精品久久久精品久久久| av黄色大香蕉| 国产白丝娇喘喷水9色精品| 精品人妻偷拍中文字幕| √禁漫天堂资源中文www| 亚洲国产色片| 久久久a久久爽久久v久久| 亚洲国产精品一区二区三区在线| 高清av免费在线| 免费黄色在线免费观看| 日韩精品免费视频一区二区三区 | 亚洲av.av天堂| 亚洲色图综合在线观看| 美女视频免费永久观看网站| 女人久久www免费人成看片| 汤姆久久久久久久影院中文字幕| 爱豆传媒免费全集在线观看| 性色avwww在线观看| 久热久热在线精品观看| 日韩精品免费视频一区二区三区 | 人人妻人人澡人人看| 色视频在线一区二区三区| 99热全是精品| 两个人免费观看高清视频| 午夜福利网站1000一区二区三区| 又粗又硬又长又爽又黄的视频| 一本大道久久a久久精品| 夜夜爽夜夜爽视频| 日本黄大片高清| 大码成人一级视频| 亚洲精品国产av成人精品| 久久久a久久爽久久v久久| 少妇高潮的动态图| 一区二区三区精品91| 国产亚洲欧美精品永久| av女优亚洲男人天堂| 久热这里只有精品99| 在线看a的网站| 777米奇影视久久| 美女国产高潮福利片在线看| 久久精品熟女亚洲av麻豆精品| 精品福利永久在线观看| www.熟女人妻精品国产 | 国产免费福利视频在线观看| 五月玫瑰六月丁香| 国内精品宾馆在线| 精品少妇内射三级| 欧美丝袜亚洲另类| 欧美人与善性xxx| 国产精品国产三级专区第一集| 亚洲国产精品专区欧美| 亚洲色图综合在线观看| 中文字幕亚洲精品专区| 亚洲人成77777在线视频| 日韩大片免费观看网站| 久久久久网色| 免费看不卡的av| 亚洲天堂av无毛| 国产精品久久久久久av不卡| 内地一区二区视频在线| 亚洲精华国产精华液的使用体验| 天堂中文最新版在线下载| 国产精品 国内视频| 亚洲国产av影院在线观看| 亚洲精品久久午夜乱码| 9热在线视频观看99| 国产伦理片在线播放av一区| 欧美人与性动交α欧美精品济南到 | 一边摸一边做爽爽视频免费| 免费黄频网站在线观看国产| 美女脱内裤让男人舔精品视频| 最近最新中文字幕免费大全7| 我要看黄色一级片免费的| 成年女人在线观看亚洲视频| 亚洲av国产av综合av卡| 成人漫画全彩无遮挡| 亚洲高清免费不卡视频| 久久99热6这里只有精品| 免费人成在线观看视频色| 99国产综合亚洲精品| av免费观看日本| 国产高清国产精品国产三级| 国产精品熟女久久久久浪| 亚洲精品一区蜜桃| 免费看光身美女| 免费高清在线观看视频在线观看| 啦啦啦在线观看免费高清www| 王馨瑶露胸无遮挡在线观看| 如何舔出高潮| 色5月婷婷丁香| 久久免费观看电影| 国产视频首页在线观看| 99热6这里只有精品| 亚洲人成网站在线观看播放| 久久久久久久久久久久大奶| 欧美bdsm另类| 国产成人精品婷婷| 日韩免费高清中文字幕av| 伊人久久国产一区二区| 欧美xxⅹ黑人| 亚洲综合色网址| 一区二区三区四区激情视频| 国产老妇伦熟女老妇高清| 天天影视国产精品| 久久人妻熟女aⅴ| 国产亚洲午夜精品一区二区久久| 婷婷成人精品国产| 久久女婷五月综合色啪小说| 最近最新中文字幕免费大全7| 狂野欧美激情性bbbbbb| 国产成人午夜福利电影在线观看| 久久女婷五月综合色啪小说| 亚洲精品久久久久久婷婷小说| 亚洲av电影在线观看一区二区三区| 一级黄片播放器| 高清视频免费观看一区二区| 日本色播在线视频| 亚洲欧美日韩卡通动漫| 精品亚洲成国产av| 满18在线观看网站| 搡老乐熟女国产| 亚洲图色成人| 成年美女黄网站色视频大全免费| 亚洲综合色惰| 午夜91福利影院| 亚洲国产看品久久| 欧美日韩精品成人综合77777| 久久午夜福利片| 97在线人人人人妻| 波野结衣二区三区在线| 亚洲精品自拍成人| 欧美日韩视频高清一区二区三区二| 天天躁夜夜躁狠狠躁躁| 日本-黄色视频高清免费观看| 制服人妻中文乱码| 在线亚洲精品国产二区图片欧美| 女性生殖器流出的白浆| 街头女战士在线观看网站| 制服人妻中文乱码| 麻豆乱淫一区二区| 人人妻人人爽人人添夜夜欢视频| 街头女战士在线观看网站| 午夜免费鲁丝| 亚洲国产成人一精品久久久| 女性生殖器流出的白浆| 精品一区二区三区四区五区乱码 | 夜夜骑夜夜射夜夜干| 亚洲第一av免费看| 国产免费一级a男人的天堂| 国产亚洲欧美精品永久| 亚洲国产色片| 午夜福利视频在线观看免费| 色94色欧美一区二区| 亚洲精品乱久久久久久| 狂野欧美激情性bbbbbb| 国产成人精品在线电影| 欧美精品国产亚洲| 中文字幕亚洲精品专区| 午夜福利网站1000一区二区三区| 国产激情久久老熟女| 下体分泌物呈黄色| tube8黄色片| 高清视频免费观看一区二区| 波多野结衣一区麻豆| 色婷婷久久久亚洲欧美| 亚洲精品美女久久久久99蜜臀 | xxx大片免费视频| 韩国av在线不卡| 不卡视频在线观看欧美| 亚洲成国产人片在线观看| 校园人妻丝袜中文字幕| av在线老鸭窝| 日本午夜av视频| 亚洲人与动物交配视频| 在线天堂最新版资源| 黄片播放在线免费| 中文字幕免费在线视频6| 在线精品无人区一区二区三| av在线观看视频网站免费| 亚洲 欧美一区二区三区| 高清在线视频一区二区三区| 永久免费av网站大全| 日日摸夜夜添夜夜爱| 国产一区有黄有色的免费视频| 纵有疾风起免费观看全集完整版| 午夜久久久在线观看| 少妇熟女欧美另类| 国产一区二区激情短视频 | 韩国精品一区二区三区 | 一级毛片我不卡| 人妻 亚洲 视频| 搡女人真爽免费视频火全软件| av电影中文网址| 国产成人午夜福利电影在线观看| 中国三级夫妇交换| 久久久久久人妻| 中文欧美无线码| 国产免费一区二区三区四区乱码| 亚洲少妇的诱惑av| 国产永久视频网站| xxxhd国产人妻xxx| 国产亚洲精品第一综合不卡 | 王馨瑶露胸无遮挡在线观看| 成人综合一区亚洲| 免费高清在线观看视频在线观看| 亚洲精品av麻豆狂野| av在线老鸭窝| 丝袜人妻中文字幕| 一二三四中文在线观看免费高清| 少妇精品久久久久久久| 亚洲av中文av极速乱| 日日爽夜夜爽网站| 巨乳人妻的诱惑在线观看| 超色免费av| 国产免费一区二区三区四区乱码| 青春草视频在线免费观看| 黄片无遮挡物在线观看| 嫩草影院入口| 大香蕉久久网| 中国国产av一级| 男人添女人高潮全过程视频| 国产精品成人在线| videossex国产| 午夜视频国产福利| 又大又黄又爽视频免费| 国产片特级美女逼逼视频| 99国产精品免费福利视频| 十分钟在线观看高清视频www| 在线 av 中文字幕| 99re6热这里在线精品视频| 亚洲三级黄色毛片| 五月开心婷婷网| 国产精品免费大片| 内地一区二区视频在线| 只有这里有精品99| 曰老女人黄片| 日本免费在线观看一区| 黄色毛片三级朝国网站| 97超碰精品成人国产| 精品午夜福利在线看| 在线亚洲精品国产二区图片欧美| 日韩一本色道免费dvd| 王馨瑶露胸无遮挡在线观看| 久久亚洲国产成人精品v| 少妇的丰满在线观看| 久久久久久人人人人人| 交换朋友夫妻互换小说| 又黄又粗又硬又大视频| 五月开心婷婷网| 国产精品女同一区二区软件| 秋霞伦理黄片| 欧美精品国产亚洲| 精品亚洲成a人片在线观看| 中文字幕人妻熟女乱码| 国产精品国产三级国产专区5o| 亚洲欧美一区二区三区黑人 | 国产亚洲午夜精品一区二区久久| 91成人精品电影| 最近最新中文字幕免费大全7| 国产成人欧美| av网站免费在线观看视频| 久久精品国产综合久久久 | 久久综合国产亚洲精品| 国产不卡av网站在线观看| 免费久久久久久久精品成人欧美视频 | 熟妇人妻不卡中文字幕| 精品国产乱码久久久久久小说| 婷婷成人精品国产| 美女大奶头黄色视频| 春色校园在线视频观看| 人人妻人人爽人人添夜夜欢视频| 一本色道久久久久久精品综合| 黄色 视频免费看| 这个男人来自地球电影免费观看 | 一区二区三区精品91| 欧美成人午夜精品| 亚洲精品久久久久久婷婷小说| 亚洲精品成人av观看孕妇| 亚洲一区二区三区欧美精品| 97人妻天天添夜夜摸| 一级,二级,三级黄色视频| 狂野欧美激情性bbbbbb| 日本黄大片高清| 最后的刺客免费高清国语| 欧美激情极品国产一区二区三区 | 国产精品人妻久久久久久| 激情五月婷婷亚洲| 国产不卡av网站在线观看| 在线观看国产h片| 18禁国产床啪视频网站| 国产亚洲最大av| 亚洲欧美成人精品一区二区| 国产视频首页在线观看| 成年动漫av网址| 久久久久视频综合| 乱人伦中国视频| 哪个播放器可以免费观看大片| 欧美xxxx性猛交bbbb| 午夜福利在线观看免费完整高清在| 卡戴珊不雅视频在线播放| 色5月婷婷丁香| 黑人猛操日本美女一级片| 丝袜喷水一区| 国产麻豆69| 亚洲国产成人一精品久久久| 大香蕉久久成人网| 最近2019中文字幕mv第一页| 亚洲第一区二区三区不卡| 美女内射精品一级片tv| 晚上一个人看的免费电影| 纵有疾风起免费观看全集完整版| 美国免费a级毛片| 90打野战视频偷拍视频| 久久精品久久精品一区二区三区| 日韩成人伦理影院| 97在线视频观看| 男女无遮挡免费网站观看| 亚洲精品乱久久久久久| 精品亚洲乱码少妇综合久久| 婷婷色综合大香蕉| 亚洲av综合色区一区| 免费高清在线观看视频在线观看| 国产 精品1| 午夜福利网站1000一区二区三区| av电影中文网址| 黑丝袜美女国产一区| 热re99久久国产66热| 黄色一级大片看看| 久久99热这里只频精品6学生| 人人妻人人爽人人添夜夜欢视频| 久久热在线av| 亚洲欧美成人精品一区二区| 在线观看国产h片| 欧美激情国产日韩精品一区| 亚洲精品久久午夜乱码| www.av在线官网国产| 精品少妇内射三级| 亚洲成人av在线免费| 欧美人与性动交α欧美软件 | 男男h啪啪无遮挡| 欧美激情极品国产一区二区三区 | 日韩精品免费视频一区二区三区 | 日韩欧美精品免费久久| 人体艺术视频欧美日本| 国产亚洲精品久久久com| 亚洲一码二码三码区别大吗| 丝袜脚勾引网站| 亚洲av中文av极速乱| 婷婷色av中文字幕| 久久久久精品久久久久真实原创| videos熟女内射| 精品午夜福利在线看| 久久狼人影院| 夫妻午夜视频| 国产成人午夜福利电影在线观看| 在线观看国产h片| 伦精品一区二区三区| 久久97久久精品| 久久精品人人爽人人爽视色| 王馨瑶露胸无遮挡在线观看| 美女脱内裤让男人舔精品视频| av不卡在线播放| 22中文网久久字幕| 丰满饥渴人妻一区二区三| 久久精品国产自在天天线| 乱人伦中国视频| 国产成人精品婷婷| kizo精华| 少妇猛男粗大的猛烈进出视频| 肉色欧美久久久久久久蜜桃| 国产黄色视频一区二区在线观看| 男人操女人黄网站| 国产亚洲av片在线观看秒播厂|