• <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é)改進
    巧選參考系簡化解題過程
    国产亚洲精品av在线| av又黄又爽大尺度在线免费看 | 亚洲欧美精品自产自拍| 亚洲人成网站在线播| 日韩大尺度精品在线看网址| 女同久久另类99精品国产91| 久久亚洲精品不卡| 欧美在线一区亚洲| 成熟少妇高潮喷水视频| 亚洲国产精品合色在线| 久久人人精品亚洲av| 亚洲av成人av| 久久久色成人| 亚洲第一电影网av| 小说图片视频综合网站| 丰满的人妻完整版| 亚洲精品久久国产高清桃花| 人妻久久中文字幕网| 国产午夜精品论理片| ponron亚洲| 国产精品一区二区三区四区久久| 午夜精品在线福利| 大又大粗又爽又黄少妇毛片口| 国产精品1区2区在线观看.| 亚洲成人中文字幕在线播放| 欧美性感艳星| 极品教师在线视频| 成人高潮视频无遮挡免费网站| 2021天堂中文幕一二区在线观| 九色成人免费人妻av| 在现免费观看毛片| 99精品在免费线老司机午夜| 国产女主播在线喷水免费视频网站 | 国产成人午夜福利电影在线观看| 亚洲精品日韩在线中文字幕 | 日本熟妇午夜| 久久精品久久久久久噜噜老黄 | 午夜福利成人在线免费观看| 免费无遮挡裸体视频| 欧美日韩精品成人综合77777| 一级毛片久久久久久久久女| 啦啦啦啦在线视频资源| 国产精品99久久久久久久久| 亚洲精品日韩av片在线观看| 少妇猛男粗大的猛烈进出视频 | 少妇的逼好多水| 日本欧美国产在线视频| 国产欧美日韩精品一区二区| 日日啪夜夜撸| 校园人妻丝袜中文字幕| 亚洲av电影不卡..在线观看| 免费av不卡在线播放| 亚洲欧美成人综合另类久久久 | 亚洲第一电影网av| 永久网站在线| 国产一区二区在线av高清观看| 亚洲欧美日韩高清专用| 可以在线观看毛片的网站| av在线蜜桃| 亚洲av成人av| 麻豆精品久久久久久蜜桃| 国产精品不卡视频一区二区| 亚洲精品久久久久久婷婷小说 | 午夜爱爱视频在线播放| 久久久欧美国产精品| 日韩精品青青久久久久久| a级毛片免费高清观看在线播放| 能在线免费观看的黄片| 97超碰精品成人国产| 波多野结衣高清无吗| 亚洲精品国产成人久久av| 亚洲人与动物交配视频| 91精品一卡2卡3卡4卡| 2022亚洲国产成人精品| av视频在线观看入口| 免费黄网站久久成人精品| 国产毛片a区久久久久| eeuss影院久久| 亚洲av熟女| 国产精品久久电影中文字幕| 全区人妻精品视频| 免费不卡的大黄色大毛片视频在线观看 | 日日撸夜夜添| 久久久久九九精品影院| 啦啦啦韩国在线观看视频| ponron亚洲| or卡值多少钱| 成熟少妇高潮喷水视频| 亚洲国产高清在线一区二区三| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 色综合色国产| 精品一区二区免费观看| 男人和女人高潮做爰伦理| 亚洲欧美日韩高清在线视频| 悠悠久久av| 久久精品综合一区二区三区| 国产中年淑女户外野战色| 亚洲成av人片在线播放无| 最好的美女福利视频网| 亚洲欧美日韩卡通动漫| 久久人妻av系列| 寂寞人妻少妇视频99o| 亚洲内射少妇av| 在线免费观看不下载黄p国产| 97超碰精品成人国产| 深夜a级毛片| 免费搜索国产男女视频| 亚洲在线自拍视频| 蜜臀久久99精品久久宅男| 欧美日韩乱码在线| 校园春色视频在线观看| 国产精品久久久久久久电影| 国产一区二区三区在线臀色熟女| 久久久国产成人精品二区| 精品日产1卡2卡| 免费无遮挡裸体视频| 99久久中文字幕三级久久日本| 亚洲欧美日韩高清专用| 成人欧美大片| 国产精品乱码一区二三区的特点| 我的女老师完整版在线观看| 人妻系列 视频| 内射极品少妇av片p| 国产不卡一卡二| 精品久久久久久久末码| 欧美一级a爱片免费观看看| 国产精品综合久久久久久久免费| 69人妻影院| 亚洲国产精品国产精品| 波多野结衣高清无吗| 狠狠狠狠99中文字幕| 在线播放国产精品三级| 国产精品久久久久久亚洲av鲁大| 国产精品久久久久久av不卡| 丝袜喷水一区| 亚洲欧美日韩无卡精品| 国产亚洲av嫩草精品影院| 大又大粗又爽又黄少妇毛片口| 国产极品天堂在线| 丝袜美腿在线中文| 国产探花在线观看一区二区| 国产精品无大码| 18禁在线无遮挡免费观看视频| 久久久a久久爽久久v久久| ponron亚洲| 天美传媒精品一区二区| 国产综合懂色| 99国产极品粉嫩在线观看| 日日摸夜夜添夜夜爱| 国产亚洲av片在线观看秒播厂 | 欧美高清性xxxxhd video| 在线天堂最新版资源| 精品久久久久久久久av| 99热这里只有精品一区| 国产乱人视频| 亚洲精品影视一区二区三区av| 啦啦啦观看免费观看视频高清| 久久欧美精品欧美久久欧美| 免费搜索国产男女视频| 免费av观看视频| 亚洲精品成人久久久久久| 国产成人影院久久av| 99精品在免费线老司机午夜| 精华霜和精华液先用哪个| 一区二区三区免费毛片| 亚洲欧美成人精品一区二区| 久久久久久久久大av| 久99久视频精品免费| 日韩国内少妇激情av| 婷婷六月久久综合丁香| 国产中年淑女户外野战色| 成年版毛片免费区| 国产国拍精品亚洲av在线观看| 日本撒尿小便嘘嘘汇集6| 久久久久久久久久久免费av| 淫秽高清视频在线观看| 九九爱精品视频在线观看| 亚洲一区二区三区色噜噜| 欧美性感艳星| 国产探花在线观看一区二区| 国产不卡一卡二| 午夜福利视频1000在线观看| 国产午夜精品论理片| 亚洲av.av天堂| kizo精华| 国产精品人妻久久久影院| 亚洲天堂国产精品一区在线| 日日摸夜夜添夜夜添av毛片| 亚洲aⅴ乱码一区二区在线播放| 免费一级毛片在线播放高清视频| 免费搜索国产男女视频| 在线观看av片永久免费下载| 精品国内亚洲2022精品成人| 简卡轻食公司| .国产精品久久| 久久婷婷人人爽人人干人人爱| 99热这里只有是精品50| 久久中文看片网| 亚洲第一电影网av| 人妻夜夜爽99麻豆av| 国产精品嫩草影院av在线观看| 国产伦理片在线播放av一区 | 2022亚洲国产成人精品| 久久久久久久久大av| 在线播放国产精品三级| 久久人人爽人人爽人人片va| 国产精品麻豆人妻色哟哟久久 | 2021天堂中文幕一二区在线观| 久久久久九九精品影院| 亚洲av免费在线观看| 黄色一级大片看看| av视频在线观看入口| 午夜视频国产福利| 欧美日韩乱码在线| 国产成人91sexporn| 天堂中文最新版在线下载 | 欧美精品一区二区大全| 中文字幕av在线有码专区| 国产精品电影一区二区三区| 国产精品一及| 啦啦啦观看免费观看视频高清| 国产精品免费一区二区三区在线| 久久久久久久久久成人| 我要看日韩黄色一级片| 黄色日韩在线| 亚洲四区av| 亚洲美女视频黄频| 日韩欧美 国产精品| 在线播放无遮挡| 级片在线观看| 中文字幕精品亚洲无线码一区| а√天堂www在线а√下载| 男女做爰动态图高潮gif福利片| 夫妻性生交免费视频一级片| 熟妇人妻久久中文字幕3abv| 一个人看的www免费观看视频| 欧美人与善性xxx| 蜜桃亚洲精品一区二区三区| 久久久久久久久久黄片| 亚洲丝袜综合中文字幕| 久久精品影院6| 久久草成人影院| 国产成人影院久久av| 女人十人毛片免费观看3o分钟| 国产麻豆成人av免费视频| 国产色婷婷99| 日本黄色视频三级网站网址| 少妇人妻精品综合一区二区 | 国产69精品久久久久777片| 精品久久久久久成人av| 国产成人一区二区在线| 国产爱豆传媒在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 在线免费观看的www视频| 日日摸夜夜添夜夜爱| 成年免费大片在线观看| 美女被艹到高潮喷水动态| 全区人妻精品视频| 五月伊人婷婷丁香| 深爱激情五月婷婷| 非洲黑人性xxxx精品又粗又长| 免费人成在线观看视频色| 国产亚洲精品久久久久久毛片| 亚洲欧洲国产日韩| 综合色丁香网| 亚洲人与动物交配视频| 别揉我奶头 嗯啊视频| 成熟少妇高潮喷水视频| 大又大粗又爽又黄少妇毛片口| 91狼人影院| 在线观看av片永久免费下载| 久久中文看片网| 日韩av在线大香蕉| 亚洲精品456在线播放app| 国产午夜精品一二区理论片| 亚洲一区二区三区色噜噜| 深夜精品福利| 男人的好看免费观看在线视频| 亚洲精品乱码久久久v下载方式| 两性午夜刺激爽爽歪歪视频在线观看| 国产熟女欧美一区二区| av国产免费在线观看| 国产精品人妻久久久影院| 最近视频中文字幕2019在线8| 嫩草影院新地址| 久久韩国三级中文字幕| videossex国产| 一级黄片播放器| 亚洲丝袜综合中文字幕| 国产精品久久久久久久电影| 哪个播放器可以免费观看大片| 成年女人永久免费观看视频| 亚洲四区av| 国产 一区 欧美 日韩| 女同久久另类99精品国产91| 少妇被粗大猛烈的视频| 国产精品日韩av在线免费观看| 99在线视频只有这里精品首页| 亚洲欧美日韩高清在线视频| 国产午夜福利久久久久久| 久久久久久久久久黄片| 日韩欧美国产在线观看| 亚洲美女视频黄频| 内地一区二区视频在线| 97人妻精品一区二区三区麻豆| 人妻少妇偷人精品九色| 国产精品久久久久久久电影| 18禁裸乳无遮挡免费网站照片| 国产探花极品一区二区| 18禁在线播放成人免费| kizo精华| 我要看日韩黄色一级片| 看黄色毛片网站| 亚洲自偷自拍三级| av在线老鸭窝| 91麻豆精品激情在线观看国产| 国产精品伦人一区二区| 久久人人精品亚洲av| 国产人妻一区二区三区在| 国产老妇女一区| 欧美在线一区亚洲| 熟妇人妻久久中文字幕3abv| 18禁在线播放成人免费| 日韩成人伦理影院| 中国国产av一级| 欧美潮喷喷水| 国产成人福利小说| 亚洲三级黄色毛片| 有码 亚洲区| 又黄又爽又刺激的免费视频.| 搡女人真爽免费视频火全软件| 国产单亲对白刺激| 91av网一区二区| 日本撒尿小便嘘嘘汇集6| 在线观看av片永久免费下载| 最近的中文字幕免费完整| 三级国产精品欧美在线观看| 亚洲精品国产av成人精品| 成年版毛片免费区| 天天躁日日操中文字幕| av天堂中文字幕网| videossex国产| 婷婷精品国产亚洲av| 狠狠狠狠99中文字幕| 熟妇人妻久久中文字幕3abv| 日日干狠狠操夜夜爽| 亚洲性久久影院| 狂野欧美激情性xxxx在线观看| 精品久久久久久久人妻蜜臀av| 亚洲精品粉嫩美女一区| 精品国内亚洲2022精品成人| 国产精品日韩av在线免费观看| 可以在线观看毛片的网站| 黑人高潮一二区| 午夜久久久久精精品| 免费av毛片视频| 九九爱精品视频在线观看| 婷婷色综合大香蕉| 欧美极品一区二区三区四区| avwww免费| 内地一区二区视频在线| 久久久久久伊人网av| 亚洲国产精品国产精品| 91麻豆精品激情在线观看国产| 精品人妻一区二区三区麻豆| 我的女老师完整版在线观看| 欧美一级a爱片免费观看看| 美女黄网站色视频| 免费观看精品视频网站| 波多野结衣高清作品| 又黄又爽又刺激的免费视频.| 亚洲一级一片aⅴ在线观看| 久久久a久久爽久久v久久| 美女国产视频在线观看| 国产亚洲av片在线观看秒播厂 | 白带黄色成豆腐渣| 亚洲精品久久国产高清桃花| 国产极品精品免费视频能看的| 亚洲欧美清纯卡通| 悠悠久久av| 麻豆国产av国片精品| 日韩欧美 国产精品| 能在线免费观看的黄片| 日韩精品青青久久久久久| 亚洲成人久久爱视频| 亚洲精品自拍成人| 夜夜夜夜夜久久久久| 亚洲成人精品中文字幕电影| 中文字幕免费在线视频6| 欧美bdsm另类| 国产69精品久久久久777片| 身体一侧抽搐| 最近中文字幕高清免费大全6| 成年av动漫网址| 亚洲国产高清在线一区二区三| 激情 狠狠 欧美| 亚洲av电影不卡..在线观看| 在线a可以看的网站| 欧美性感艳星| 亚洲人成网站高清观看| 精华霜和精华液先用哪个| 在线免费十八禁| 欧美最新免费一区二区三区| 深爱激情五月婷婷| 成人欧美大片| 菩萨蛮人人尽说江南好唐韦庄 | 精品久久久久久久末码| 中文资源天堂在线| 国产高清视频在线观看网站| 中文字幕免费在线视频6| 蜜臀久久99精品久久宅男| 麻豆国产av国片精品| 色综合色国产| 国产成人影院久久av| 国内精品美女久久久久久| 久久久久网色| 国产爱豆传媒在线观看| 啦啦啦啦在线视频资源| 久久这里只有精品中国| 久久久精品大字幕| 久久久精品94久久精品| 国产成人一区二区在线| 嫩草影院入口| 我的老师免费观看完整版| 最近2019中文字幕mv第一页| 欧美日本亚洲视频在线播放| 老司机福利观看| 国产精品福利在线免费观看| 校园春色视频在线观看| 亚洲av.av天堂| 国产色爽女视频免费观看| 精品久久久久久久久av| 中国美白少妇内射xxxbb| 日日干狠狠操夜夜爽| av福利片在线观看| 男人和女人高潮做爰伦理| 国产精品久久久久久久电影| 欧美xxxx性猛交bbbb| 国产亚洲5aaaaa淫片| 久久精品国产99精品国产亚洲性色| av免费观看日本| 国产91av在线免费观看| 久久久久久伊人网av| 啦啦啦观看免费观看视频高清| АⅤ资源中文在线天堂| 欧美变态另类bdsm刘玥| 嘟嘟电影网在线观看| 91av网一区二区| 国产精品永久免费网站| 亚州av有码| 美女被艹到高潮喷水动态| 综合色av麻豆| 国产免费男女视频| 免费搜索国产男女视频| 永久网站在线| 韩国av在线不卡| eeuss影院久久| 丝袜喷水一区| 91久久精品国产一区二区三区| 日韩在线高清观看一区二区三区| 欧美zozozo另类| 三级男女做爰猛烈吃奶摸视频| 久久国内精品自在自线图片| 91久久精品电影网| 国产麻豆成人av免费视频| 只有这里有精品99| 欧洲精品卡2卡3卡4卡5卡区| 爱豆传媒免费全集在线观看| 午夜视频国产福利| 一夜夜www| 99视频精品全部免费 在线| 老司机影院成人| 久久欧美精品欧美久久欧美| 精品99又大又爽又粗少妇毛片| 亚洲av免费在线观看| 在线国产一区二区在线| 精品免费久久久久久久清纯| 中文亚洲av片在线观看爽| 国产真实乱freesex| 最好的美女福利视频网| 亚洲久久久久久中文字幕| av又黄又爽大尺度在线免费看 | 婷婷精品国产亚洲av| 好男人在线观看高清免费视频| 波多野结衣高清无吗| 一卡2卡三卡四卡精品乱码亚洲| 国产精品一及| 亚洲成av人片在线播放无| 两个人视频免费观看高清| 成人亚洲欧美一区二区av| 麻豆国产av国片精品| 亚洲欧美精品专区久久| 日本在线视频免费播放| 有码 亚洲区| 少妇熟女aⅴ在线视频| 99久久精品国产国产毛片| 男的添女的下面高潮视频| 91精品国产九色| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成人国产麻豆网| 国内精品久久久久精免费| 久久久精品欧美日韩精品| 黄色配什么色好看| 麻豆国产av国片精品| 久久99蜜桃精品久久| 级片在线观看| 久久久精品欧美日韩精品| 婷婷精品国产亚洲av| 久99久视频精品免费| 少妇被粗大猛烈的视频| 99热精品在线国产| 久久综合国产亚洲精品| 夜夜看夜夜爽夜夜摸| 午夜激情欧美在线| 男女做爰动态图高潮gif福利片| 久久这里有精品视频免费| 一夜夜www| 国产极品精品免费视频能看的| 色视频www国产| 99久久精品热视频| 在线国产一区二区在线| 一区二区三区免费毛片| 国产乱人视频| 免费黄网站久久成人精品| 又粗又硬又长又爽又黄的视频 | 麻豆国产97在线/欧美| 国产探花极品一区二区| 老司机影院成人| 五月伊人婷婷丁香| 色尼玛亚洲综合影院| a级毛片免费高清观看在线播放| 久久久久国产网址| 国产老妇伦熟女老妇高清| 国产真实伦视频高清在线观看| 亚洲第一电影网av| 嫩草影院精品99| 人体艺术视频欧美日本| 久久精品久久久久久久性| 免费不卡的大黄色大毛片视频在线观看 | 精品久久国产蜜桃| 三级经典国产精品| 精品无人区乱码1区二区| 国产精品久久久久久精品电影| 波多野结衣巨乳人妻| av视频在线观看入口| 国产女主播在线喷水免费视频网站 | 国产成人精品婷婷| 亚洲国产日韩欧美精品在线观看| 少妇猛男粗大的猛烈进出视频 | 日本三级黄在线观看| 精品久久国产蜜桃| 免费av不卡在线播放| 亚洲人与动物交配视频| 99久久精品热视频| 久久久久免费精品人妻一区二区| 久久国内精品自在自线图片| 成人三级黄色视频| 成人一区二区视频在线观看| 欧美激情在线99| 男女啪啪激烈高潮av片| 高清毛片免费看| 久久久久久久午夜电影| 精品不卡国产一区二区三区| 成人鲁丝片一二三区免费| 亚洲熟妇中文字幕五十中出| 亚洲欧美日韩东京热| 国产老妇伦熟女老妇高清| 最新中文字幕久久久久| 成人午夜精彩视频在线观看| 亚洲五月天丁香| 此物有八面人人有两片| av黄色大香蕉| 日韩成人av中文字幕在线观看| 我的老师免费观看完整版| 亚洲人成网站在线播| 国国产精品蜜臀av免费| 国产精品一区二区在线观看99 | 亚洲av一区综合| 亚洲人成网站在线播放欧美日韩| 成人二区视频| 久久精品影院6| 看十八女毛片水多多多| 一级黄片播放器| 午夜a级毛片| 国产 一区 欧美 日韩| 久久久久久国产a免费观看| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品久久久久久精品电影小说 | 亚洲成人久久性| 精品99又大又爽又粗少妇毛片| 国产精品爽爽va在线观看网站| 亚洲国产高清在线一区二区三| 女的被弄到高潮叫床怎么办| 色哟哟·www| 夫妻性生交免费视频一级片| 干丝袜人妻中文字幕| 国产一区二区在线av高清观看| 一级黄色大片毛片| 免费看日本二区| 欧美一区二区国产精品久久精品| 高清午夜精品一区二区三区 | 日韩欧美三级三区| www.av在线官网国产| 少妇熟女欧美另类| 免费不卡的大黄色大毛片视频在线观看 | 老司机福利观看| 噜噜噜噜噜久久久久久91| 国产熟女欧美一区二区| 人妻久久中文字幕网| 我要看日韩黄色一级片| 国产免费一级a男人的天堂| 欧美高清性xxxxhd video| 嫩草影院精品99| 亚洲成人久久性| 中文在线观看免费www的网站|