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

    基于RMPC的高超聲速飛行器輸入飽和控制*

    2016-08-09 18:54:09胡超芳任艷麗解倩倩
    航天控制 2016年2期
    關(guān)鍵詞:張量積

    胡超芳 任艷麗 解倩倩

    天津大學(xué)電氣與自動(dòng)化工程學(xué)院,天津 300072

    ?

    基于RMPC的高超聲速飛行器輸入飽和控制*

    胡超芳 任艷麗 解倩倩

    天津大學(xué)電氣與自動(dòng)化工程學(xué)院,天津 300072

    針對(duì)帶有輸入飽和限制的高超聲速飛行器非線性縱向運(yùn)動(dòng)模型,提出了一種魯棒模型預(yù)測(cè)控制器(RMPC)設(shè)計(jì)方法。利用雅克比線性化以及張量積(T-P)模型轉(zhuǎn)換方法,將高超聲速飛行器非線性模型轉(zhuǎn)化為多胞線性參變(LPV)模型。在此基礎(chǔ)上,將輸入飽和表示為實(shí)際反饋控制律與輔助反饋控制律構(gòu)成的凸包,建立飽和RMPC控制器,并通過(guò)引入輔助矩陣來(lái)降低其保守性,利用線性矩陣不等式(LMI)求解,以保證輸入飽和條件下閉環(huán)系統(tǒng)的控制性能和穩(wěn)定性。通過(guò)與其它RMPC控制器的仿真比較,驗(yàn)證了本文方法的有效性。 關(guān)鍵詞 高超聲速飛行器;輸入飽和;魯棒預(yù)測(cè)控制;張量積;線性參變模型

    高超聲速飛行器是指一類飛行速度超過(guò)5倍聲速的飛行器。其機(jī)體/發(fā)動(dòng)機(jī)一體化的設(shè)計(jì),導(dǎo)致高超聲速飛行器的非線性程度較高,對(duì)控制器設(shè)計(jì)提出了極大挑戰(zhàn)。為避免復(fù)雜的非線性設(shè)計(jì)過(guò)程,并實(shí)現(xiàn)高超聲速飛行器大包絡(luò)下的飛行控制,鑒于線性參變(LPV)模型能夠充分描述系統(tǒng)內(nèi)在的非線性和時(shí)變特點(diǎn),可以利用成熟線性控制理論的優(yōu)勢(shì)[1]。近年來(lái)基于LPV系統(tǒng)的高超聲速飛行器控制研究得到了一定發(fā)展,如文獻(xiàn)[2-3]針對(duì)高超聲速飛行器縱向運(yùn)動(dòng)模型,在多胞LPV模型基礎(chǔ)上設(shè)計(jì)了魯棒變?cè)鲆婵刂破?。此外,輸入飽和作為控制系統(tǒng)中較為常見(jiàn)的一種現(xiàn)象,嚴(yán)重影響著飛行控制系統(tǒng)的閉環(huán)性能和穩(wěn)定性[4],對(duì)于高超聲速飛行器的安全穩(wěn)定飛行具有重大現(xiàn)實(shí)意義,引起了廣大學(xué)者的關(guān)注。如Xu等[5]基于動(dòng)態(tài)面控制,利用指令濾波以及跟蹤誤差補(bǔ)償來(lái)解決輸入飽和問(wèn)題;Huang等[6]采用積分二次約束對(duì)高超聲速飛行器的輸入飽和非線性進(jìn)行處理,實(shí)現(xiàn)了平衡點(diǎn)處較小范圍的穩(wěn)定控制。

    魯棒預(yù)測(cè)控制(RMPC)由于其突出的顯式處理約束能力以及良好的魯棒控制性能[7-8],在飛行控制領(lǐng)域也得到了一定的發(fā)展和應(yīng)用。如文獻(xiàn)[9]針對(duì)高超聲速飛行器,設(shè)計(jì)了基于多胞LPV模型的RMPC控制器。而對(duì)于輸入飽和限制要求,目前也已有關(guān)于RMPC的研究成果,如文獻(xiàn)[10]基于多胞LPV模型,設(shè)計(jì)輔助線性反饋律,構(gòu)建RMPC控制器,來(lái)解決輸入飽和問(wèn)題;Huang等[11]對(duì)該飽和方法進(jìn)行了改進(jìn),通過(guò)引入權(quán)重調(diào)節(jié)參數(shù)降低了RMPC控制器的保守性。

    本文針對(duì)帶有輸入飽和限制的高超聲速飛行器縱向運(yùn)動(dòng)模型,設(shè)計(jì)飽和RMPC控制器。首先,基于雅克比線性化和張量積模型轉(zhuǎn)換方法,將非線性系統(tǒng)轉(zhuǎn)化為多胞LPV模型;然后,利用實(shí)際反饋控制律與輔助反饋控制律構(gòu)成的凸集逼近飽和輸入,并引入輔助矩陣來(lái)提高控制器設(shè)計(jì)自由度、降低保守性,以建立飽和RMPC控制器,并通過(guò)線性矩陣不等式(LMI)求解,保證了飛行狀態(tài)的穩(wěn)定收斂。仿真結(jié)果驗(yàn)證了所提方法的有效性。

    1 模型描述

    為了合理簡(jiǎn)化模型的復(fù)雜度,不失一般性,考慮如下假設(shè)[12]:

    1)高超聲速飛行器為理想剛性體,質(zhì)心保持在機(jī)體縱軸上,推力沿縱軸方向;

    2)不考慮地球自轉(zhuǎn),假設(shè)地面局部平坦且重力加速度恒定;

    3)大氣相對(duì)地球靜止,忽略氣流的不對(duì)稱性及壓縮性、陣風(fēng)等因素的影響;

    4)忽略飛行器橫側(cè)向運(yùn)動(dòng)對(duì)縱向運(yùn)動(dòng)的影響,只考慮縱向運(yùn)動(dòng)。

    基于此,高超聲速飛行器的縱向運(yùn)動(dòng)模型[13-14]為:

    (1)

    其中,h,V,α,θ,q分別為高度、速度、攻角、航跡角和俯仰速率,m,g分別為飛行器的質(zhì)量和重力加速度,Iyy為飛行器相對(duì)機(jī)體坐標(biāo)系y軸的轉(zhuǎn)動(dòng)慣量,L,D,T,M分別表示升力、阻力、推力和俯仰力矩。

    (2)

    ρ為空氣密度,氣動(dòng)參數(shù)表達(dá)式為:

    (4)

    (5)

    umax,umin分別為飽和上下限值,ε為所設(shè)計(jì)控制律。

    2 多胞LPV系統(tǒng)建模

    為建立高超聲速飛行器的多胞LPV模型,首先需要構(gòu)建其連續(xù)LPV模型。通過(guò)求解以下方程獲得平衡點(diǎn):

    f1=Vsin(θ-α),

    然后利用雅克比線性化(也就是小擾動(dòng)線性化[15]),在平衡點(diǎn)附近,將縱向運(yùn)動(dòng)模型通過(guò)Taylor展開(kāi)式化簡(jiǎn),略去二階以上無(wú)窮小量,得到近似的線性縱向運(yùn)動(dòng)模型。于是高超聲速飛行器的動(dòng)態(tài)特性就可用平衡點(diǎn)附近的線性系統(tǒng)來(lái)描述。

    為此,對(duì)高超聲速飛行器平衡點(diǎn)進(jìn)行分析,得到雅克比線性化后的LPV模型:

    (6)

    其中,Δx=x-xe,Δu=u-ue,xe,ue為平衡點(diǎn)處的狀態(tài)和輸入,由f(xe,ue)=0求得。p(t)=[hV]T為變參數(shù)向量,h,V為調(diào)度變量,且p(t)∈Γ?R2,Γ=[2286.00m/s,2590.80m/s]×[25908.00m,26517.60m]是由期望飛行速度范圍和高度范圍所確定的有界閉域。其次,使用張量積模型轉(zhuǎn)換方法獲得式(6)在閉域Γ內(nèi)的多胞模型,其主要包括以下步驟:1)在參數(shù)空間Γ內(nèi),分別在高度和速度兩個(gè)維度上選取一定數(shù)量的采樣點(diǎn),將LPV系統(tǒng)離散化,并建立張量以存儲(chǔ)以上離散系統(tǒng);2)通過(guò)高階奇異值分解(HOSVD)[16-17],提取核心線性時(shí)不變(LTI)頂點(diǎn)系統(tǒng);3)構(gòu)建LTI頂點(diǎn)系統(tǒng)的權(quán)系數(shù)函數(shù)。本文中,取樣密度為21×21,通過(guò)以上過(guò)程,時(shí)變矩陣S(p(t))=[A(p(t))B(p(t))]可以近似表達(dá)為變參數(shù)p(t)與有限個(gè)LTI頂點(diǎn)系統(tǒng)的凸組合

    (7)

    3 飽和RMPC控制器設(shè)計(jì)

    Δx(k+1)=A(τ(k))Δx(k)+B(τ(k))Δu(k)

    (8)

    [A(τ(k))B(τ(k))]∈Ω=

    co{[A1B1],[A2B2],...,[A20B20]}

    其中,Ω表示頂點(diǎn)[AnBn]構(gòu)成的凸包,即

    [A(τ(k))B(τ(k))]=

    為考慮控制輸入的飽和約束(5),基于系統(tǒng)(8),采用狀態(tài)反饋來(lái)設(shè)計(jì)飽和RMPC控制器。首先,對(duì)于輸入飽和限制,將Δu表示為凸包的形式[18]

    Δu=sat(KΔx,umin-ur,umax-ur)∈

    (9)

    sat(KΔx,umin-ur,umax-ur)=

    (10)

    然后,在多胞LPV模型(8)的每一頂點(diǎn)處均采用一致的狀態(tài)反饋控制律(10)。Δx(k+i|k)為系統(tǒng)在k時(shí)刻得到的k+i時(shí)刻的預(yù)測(cè)狀態(tài),Δu(k+i|k)為k+i時(shí)刻的預(yù)測(cè)控制輸入。取Lyapunov函數(shù)為V(k+i|k)=Δx(k+i|k)TPΔx(k+i|k),P>0,P為鎮(zhèn)定矩陣。設(shè)k時(shí)刻Lyapunov函數(shù)V(k|k)=Δx(k|k)TPΔx(k|k)的上界值為γ,則目標(biāo)函數(shù)J∞(k)為

    其中,ψ,σ為正定加權(quán)矩陣。

    下面給出定理1來(lái)求解以上飽和RMPC控制器,以保證閉環(huán)系統(tǒng)穩(wěn)定性。

    定理1 針對(duì)多胞LPV模型(8)以及飽和控制律(10),給定初始狀態(tài)Δx(k|k)、正定矩陣ψ,σ以及權(quán)重調(diào)節(jié)參數(shù)β>0,且令γ1=aγ。k時(shí)刻若存在0

    優(yōu)化問(wèn)題1:

    (11)

    (12)

    (13)

    γ-γ1>0

    (14)

    (15)

    (16)

    證明:對(duì)式(11)應(yīng)用Schur補(bǔ)定理[19]可得

    Δx(k|k)TPΔx(k|k)≤γ

    (17)

    本文引入輔助矩陣G來(lái)提高控制器設(shè)計(jì)自由度,降低保守性[20],根據(jù)文獻(xiàn)[21]有G+GT-Q≤GTQ-1G成立,則式(12)和(13)可分別寫(xiě)為式(18)和(19)

    (18)

    (19)

    將式(18)兩邊乘以λ1,式(19)兩邊乘以λj,j=2,3,4,并將兩式左右分別相加可得:

    n= 1,...,20 (20)

    (22)

    式(22)兩邊同時(shí)乘以γ,Δx(k+i|k)T和Δx(k+i|k),令P=γQ-1,K=YG-1,H=ZG-1,則有

    ||AnΔx(k+i|k)+BnΔu(k+i|k)||P-

    (23)

    式(23)保證了Lyapunov函數(shù)遞減,進(jìn)而結(jié)合式(17)可知,Z={Δx(k+i|k)|Δx(k+i|k)TPΔx(k+i|k)≤γ}可作為模型(8)的不變集[8],且該不變集包含原點(diǎn)。

    (24)

    (25)

    (26)

    對(duì)式(25)和(26)應(yīng)用Schur補(bǔ)定理可得:

    (27)

    (28)

    令Hl=ZlG-1,P=γQ-1,并對(duì)式(27)和(28)分別左乘diag{I,GT}和右乘diag{I,G},可得

    (29)

    (30)

    又因?yàn)镚+GT-Q≤GTQ-1G,則有式(15)和(16)成立。此外,若k時(shí)刻存在可行解,則該解對(duì)于k+i,i≥0時(shí)刻仍然可行[8],使得閉環(huán)系統(tǒng)穩(wěn)定。

    4 仿真

    為驗(yàn)證本文控制器的有效性,對(duì)飽和條件下的高超聲速飛行器縱向運(yùn)動(dòng)模型進(jìn)行仿真。飛行參考點(diǎn)取為hr=26243.28m,Vr=2423.16m/s,αr=0.782°,θr=0.782°,qr=0,δer=10.627,Φcr=0.215,初始點(diǎn)取為h0=26060.40m,V0=2286.00m/s,α0=0.844°,θ0=0.844°,q0=0,輸入飽和約束的上下界值分別為δemin=-30°,δemax=30°;Φcmin=0,Φcmax=1.5。

    首先進(jìn)行多胞LPV建模,然后選取合適的矩陣ψ,σ和權(quán)重調(diào)節(jié)參數(shù)β,設(shè)計(jì)飽和RMPC控制器。并將其與文獻(xiàn)[8]中的經(jīng)典約束RMPC以及常規(guī)無(wú)約束RMPC進(jìn)行對(duì)比,仿真結(jié)果如圖1~7所示。

    圖1 舵偏角

    圖2 空氣當(dāng)量比

    圖3 高度

    圖4 速度

    圖5 攻角

    圖6 俯仰角

    圖7 俯仰速率

    從圖1和2可以看出,飽和RMPC與經(jīng)典約束RMPC[8]均能保證輸入滿足給定的飽和限制,而常規(guī)無(wú)約束RMPC控制器使得系統(tǒng)輸入幅值過(guò)大,嚴(yán)重超過(guò)飛行器執(zhí)行機(jī)構(gòu)物理約束范圍。由圖3和4可知,與經(jīng)典約束RMPC[8]相比,本文的飽和RMPC控制器設(shè)計(jì)方法保守性相對(duì)較低,高度超調(diào)量較小,速度上升較快,系統(tǒng)能夠相對(duì)較快地到達(dá)參考點(diǎn)。圖5~7分別表示了攻角α、俯仰角θ和俯仰速率q變化曲線,無(wú)約束RMPC下的狀態(tài)變量變化幅值過(guò)大。而本文飽和RMPC下的α和θ略大于經(jīng)典約束RMPC[8]方法,這是因?yàn)榍罢叩目刂屏肯鄬?duì)于后者較大,但總體變化仍然合理,到達(dá)穩(wěn)態(tài)較快。綜上分析,本文提出的飽和RMPC控制器是有效的,具有良好的穩(wěn)定性和控制性能,實(shí)現(xiàn)了二者的兼顧。

    5 結(jié)論

    針對(duì)帶有輸入飽和的高超聲速飛行器縱向運(yùn)動(dòng)模型,提出了一種飽和RMPC控制器設(shè)計(jì)方法。利用雅克比線性化以及張量積模型轉(zhuǎn)換方法實(shí)現(xiàn)了非線性模型的線性化,建立了多胞LPV系統(tǒng),并通過(guò)構(gòu)造實(shí)際與輔助狀態(tài)反饋律,設(shè)計(jì)了飽和RMPC控制器,其中輔助矩陣的引入一定程度上降低了控制器的保守性,使得系統(tǒng)既滿足實(shí)際輸入飽和限制,又實(shí)現(xiàn)了飛行器快速穩(wěn)定的控制性能,為抗飽和控制方法研究提供了新的思路。

    [1] 王明昊, 劉剛, 趙樹(shù)濤, 等.高超聲速飛行器的LPV變?cè)鲆鏍顟B(tài)反饋H∞控制 [J].宇航學(xué)報(bào), 2013, 34(4): 488-495.(Wang Minghao, Liu Gang, Zhao Shupeng, et al. Variable gain state feedbackH∞control for hypersonic vehicle based on LPV [J]. Journal of Astronautics, 2013, 34(4): 488-495.)

    [2] Sun C Y, Huang Y Q, Qian C S, et al. On modeling and control of a flexible air-breathing hypersonic vehicle based on LPV method [J]. Frontiers of Electrical and Electronic Engineering, 2012, 7(1): 56-68.

    [3] 秦偉偉, 鄭志強(qiáng), 劉剛, 等.高超聲速飛行器的LPV魯棒變?cè)鲆婵刂?[J].系統(tǒng)工程與電子技術(shù), 2011, 33(6): 1327-1331.(Qin Weiwei, Zheng Zhiqiang, Liu Gang, et al. Robust variable gain control for hypersonic vehicles based on LPV [J]. Systems Engineering and Electronics, 2011, 33(6): 1327-1331.)

    [4] Hu T S, Lin Z L. Control systems with actuator saturation: analysis and design [M].Boston: Birkh?user, 2001.

    [5] Xu B, Wang S, Gao D X, et al. Command filter based robust nonlinear control of hypersonic aircraft with magnitude constraints on states and actuators [J]. Journal of Intelligent and Robotic Systems, 2014, 73(1-4): 233-247.

    [6] Huang X L, Ge D M. Robust and anti-windup control of air-breathing hypersonic cruise vehicle [C]. 3rd International Symposium on Systems and Control in Aeronautics and Astronautics, Harbin, 2010, 1515-1520.

    [7] 席裕庚, 李德偉, 林姝.模型預(yù)測(cè)控制-現(xiàn)狀與挑戰(zhàn) [J].自動(dòng)化學(xué)報(bào), 2013, 39(3): 222-236.(Xi Yugeng, Li Dewei, Lin Shu. Model predictive control-status and challenges [J]. Acta Automatica Sinica, 2013, 39(3): 222-236.)

    [8] Kothare M V, Balakrishnan V, Morari M. Robust constrained model predictive control using linear matrix inequalities [J]. Automatica, 1996, 32(10): 1361-1379.

    [9] Qin W W, Zheng Z Q, Zhang L, et al. Robust model predictive control for hypersonic vehicle based on LPV [C]. Proceedings of the IEEE International Conference on Information and Automation, Changsha, 2010, 1012-1017.

    [10] Cao Y Y, Lin Z L. Min-max MPC algorithm for LPV systems subject to input saturation [J]. IEE Proceedings: Control Theory and Applications, 2005, 152(3): 266-272.

    [11] Huang H, Li D W, Lin Z L, et al. An improved robust model predictive control design in the presence of actuator saturation [J]. Automatica, 2011, 47(4): 861-864.

    [12] 程路, 姜長(zhǎng)生, 文杰, 等.近空間飛行器飛/推一體化模糊自適應(yīng)廣義預(yù)測(cè)控制[J].系統(tǒng)工程與電子技術(shù), 2011, 33(1): 127-0133.

    [13] Gao G, Wang J Z. Reference command tracking control for an air-breathing hypersonic vehicle with parametric uncertainties [J]. Journal of the Franklin Institute, 2013, 350(5): 1155-1188.

    [14] Parker J T, Serrani A, Yurkovich S, et al. Control-oriented modeling of an air-breathing hypersonic vehicle[J]. Journal of Guidance, Control, and Dynamics, 2007, 30(3): 856-869.

    [15] 吳森堂, 費(fèi)玉華.飛行控制系統(tǒng)[M].北京航空航天大學(xué)出版社, 2005.

    [16] Petres Z. Polytopic decomposition of linear parameter-varying models by tensor-product model transformation [D]. Budapest, Hungary: University of Technology and Economics, 2006.

    [17] Lathauwer L, Moor B, Vandewalle J. A multi-linear singular value decomposition [J]. SIAM Journal on Matrix Analysis and Applications, 2000, 21(4): 1253-1278.

    [18] Shi T, Su H Y, Chu J. An improved model predictive control for uncertain systems with input saturation [J]. Journal of the Franklin Institute, 2013, 350(9): 2757-2768.

    [19] 俞立.魯棒控制—線性矩陣不等式處理方法 [M]. 北京: 清華大學(xué)出版社, 2002.(Yu Li. Robust control—Linear matrix inequality approach [M].Beijing: Tsinghua University Press, 2002.)

    [20] Wada N, Saito K, Saeki M. Model predictive control for linear parameter varying systems using parameter dependent Lyapunov function [J]. IEEE Transactions on Circuits and Systems II: Express Briefs, 2006, 53(12): 1446-1450.

    [21] Daafouz J, Bernussou J. Parameter dependent Lyapunov functions for discrete time systems with time varying parametric uncertainties [J]. Systems & Control Letters, 2001, 43(5): 355-359.

    Input Saturation Control Based on RMPC for Hypersonic Vehicles

    Hu Chaofang, Ren Yanli, Xie Qianqian

    School of Electrical Engineering and Automation, Tianjin University, Tianjin 300072, China

    Arobustmodelpredictivecontroller(RMPC)designmethodisproposedforthenonlinearlongitudinaldynamicmodelofthehypersonicvehiclewithinputsaturation.Thepoly-topiclinearparametervarying(LPV)modelofthevehiclecanbeobtainedbyusingJacobianlinearizationandtensor-product(TP)modeltransformationapproach.Basedonthismodel,thesaturatedRMPCcontrollerisdesigned,whichisbasedontheconvexhullconstructedbytheactualandauxiliarystatefeedbacklawstoreplacetheinputsaturation,andtheauxiliarymatrixisintroducedtodecreasetheconservatismofthecontroller.Thelinearmatrixinequality(LMI)isusedtosolvethecontroller.Thecontrolperformanceandstabilityoftheclose-loopsystemareguaranteed.TheeffectivenessoftheproposedmethodisdemonstratedbycomparingwithotherRMPCcontrollerssimulation.

    Hypersonicvehicle;Inputsaturation;Robustmodelpredictivecontrol;Tensor-product;Linearparametervaryingmodel

    *天津市自然科學(xué)基金項(xiàng)目(12JCZDJC30300); 國(guó)家留學(xué)基金

    2015-08-27

    胡超芳(1973-),男,河北保定人,副教授,主要從事預(yù)測(cè)控制、飛行器優(yōu)化與控制等研究;任艷麗(1991-),女,山西呂梁人,碩士研究生,主要從事魯棒預(yù)測(cè)控制研究;解倩倩(1990-),女,河北滄州人,碩士研究生,主要從事魯棒預(yù)測(cè)控制研究 。

    TP273

    A

    1006-3242(2016)02-0020-07

    猜你喜歡
    張量積
    B3到B16的一類逆緊全純式
    A3型李代數(shù)的張量積分解
    四種半張量積及其代數(shù)關(guān)系
    Gorenstein投射模的張量積
    Hilbert空間的張量積的連續(xù)性
    非對(duì)稱量子乘積-張量積碼
    余代數(shù)的迭代余扭曲張量積
    次對(duì)合矩陣及其性質(zhì)
    張量積形式的三維延拓Kantorovich法求解彈性動(dòng)力學(xué)問(wèn)題
    有限生成G-投射模的張量積
    精华霜和精华液先用哪个| 天美传媒精品一区二区| 国产午夜精品久久久久久一区二区三区| 精品人妻偷拍中文字幕| 精品酒店卫生间| 国产伦理片在线播放av一区| 国产免费又黄又爽又色| 亚洲精品乱码久久久v下载方式| 国产爽快片一区二区三区| 黄色怎么调成土黄色| 亚洲欧洲日产国产| 久久久国产一区二区| 亚洲综合色惰| 亚洲欧美成人综合另类久久久| 日本与韩国留学比较| 亚洲精品久久久久久婷婷小说| 国内少妇人妻偷人精品xxx网站| 久久亚洲国产成人精品v| 久久这里有精品视频免费| av黄色大香蕉| 中文字幕精品免费在线观看视频 | 日本wwww免费看| av播播在线观看一区| 日本午夜av视频| 国产精品成人在线| 亚洲激情五月婷婷啪啪| av在线蜜桃| 高清午夜精品一区二区三区| 国产熟女欧美一区二区| 最近中文字幕2019免费版| 国产成人freesex在线| 国产毛片在线视频| 亚洲第一av免费看| 亚洲精品一区蜜桃| 国产深夜福利视频在线观看| 亚洲精品国产av蜜桃| 青春草国产在线视频| a 毛片基地| 美女内射精品一级片tv| 久久6这里有精品| 国产黄片视频在线免费观看| 久久久精品94久久精品| 亚洲av男天堂| 色婷婷av一区二区三区视频| 深夜a级毛片| 看非洲黑人一级黄片| 自拍欧美九色日韩亚洲蝌蚪91 | 一级毛片 在线播放| 3wmmmm亚洲av在线观看| 天堂中文最新版在线下载| 久久国产精品大桥未久av | av视频免费观看在线观看| av线在线观看网站| 国产一级毛片在线| 国产精品伦人一区二区| 亚洲欧美一区二区三区国产| 国产亚洲精品久久久com| 一本一本综合久久| 婷婷色综合www| 我要看黄色一级片免费的| 最近中文字幕高清免费大全6| 国产精品无大码| 三级经典国产精品| 久久久a久久爽久久v久久| 久久久久久人妻| 亚洲精品国产av蜜桃| 熟妇人妻不卡中文字幕| 蜜桃久久精品国产亚洲av| 国产精品久久久久久精品古装| 午夜激情福利司机影院| 丝袜脚勾引网站| 国产视频内射| 日韩av不卡免费在线播放| 欧美3d第一页| 色视频www国产| 少妇 在线观看| 国产无遮挡羞羞视频在线观看| 91精品国产九色| 韩国高清视频一区二区三区| 韩国av在线不卡| 青春草国产在线视频| 一级黄片播放器| 超碰97精品在线观看| 女的被弄到高潮叫床怎么办| 香蕉精品网在线| 亚洲精品国产av蜜桃| 久久精品国产a三级三级三级| 中文字幕制服av| 免费黄网站久久成人精品| 亚洲人成网站高清观看| 久久婷婷青草| 午夜免费观看性视频| 赤兔流量卡办理| 国产一区二区三区综合在线观看 | 国产精品国产三级国产专区5o| 亚洲人与动物交配视频| 99热全是精品| 中文在线观看免费www的网站| 精品一区二区三卡| 欧美变态另类bdsm刘玥| 中国美白少妇内射xxxbb| 一个人看视频在线观看www免费| 夫妻性生交免费视频一级片| 黑丝袜美女国产一区| 国产av码专区亚洲av| 国产爱豆传媒在线观看| 亚洲精品日韩av片在线观看| 亚洲av.av天堂| 夫妻午夜视频| 高清不卡的av网站| 春色校园在线视频观看| 色婷婷久久久亚洲欧美| 国产精品久久久久久精品电影小说 | 爱豆传媒免费全集在线观看| 国精品久久久久久国模美| 日本欧美国产在线视频| 一级黄片播放器| 亚洲精品456在线播放app| 伊人久久国产一区二区| 久久精品熟女亚洲av麻豆精品| 日日摸夜夜添夜夜添av毛片| 久久鲁丝午夜福利片| 男女免费视频国产| 国产一区二区三区av在线| 18禁裸乳无遮挡动漫免费视频| 欧美性感艳星| 亚洲国产最新在线播放| 精品一区二区三区视频在线| 国产69精品久久久久777片| 极品少妇高潮喷水抽搐| 秋霞伦理黄片| 欧美人与善性xxx| av福利片在线观看| 亚洲国产高清在线一区二区三| 汤姆久久久久久久影院中文字幕| av国产免费在线观看| www.色视频.com| 国产精品欧美亚洲77777| 黑人猛操日本美女一级片| 一二三四中文在线观看免费高清| 日韩免费高清中文字幕av| 日本黄色片子视频| 亚洲国产精品国产精品| 成人午夜精彩视频在线观看| 成人高潮视频无遮挡免费网站| 久久久久久久久久久免费av| 男女免费视频国产| 国产免费一级a男人的天堂| 国产伦精品一区二区三区视频9| 狂野欧美激情性xxxx在线观看| 男男h啪啪无遮挡| 国产 精品1| 欧美日韩亚洲高清精品| 国产一级毛片在线| 国产成人精品一,二区| 国产成人一区二区在线| 国产成人精品久久久久久| 五月伊人婷婷丁香| 久久人人爽人人爽人人片va| av国产精品久久久久影院| 亚洲高清免费不卡视频| 麻豆精品久久久久久蜜桃| 亚洲图色成人| 在线播放无遮挡| 黑丝袜美女国产一区| 日韩欧美 国产精品| 色吧在线观看| 久久青草综合色| 亚洲精品日本国产第一区| 久久久久久久国产电影| 人妻 亚洲 视频| 欧美日韩综合久久久久久| 午夜福利在线在线| 97超视频在线观看视频| 我要看日韩黄色一级片| 大片免费播放器 马上看| 51国产日韩欧美| 卡戴珊不雅视频在线播放| 国产精品熟女久久久久浪| 超碰97精品在线观看| 成人一区二区视频在线观看| 亚洲最大成人中文| 汤姆久久久久久久影院中文字幕| av免费在线看不卡| 国产av国产精品国产| 亚洲成人一二三区av| 男人添女人高潮全过程视频| 黄色一级大片看看| 久久久久久九九精品二区国产| 日日撸夜夜添| 一区二区三区四区激情视频| 久久精品夜色国产| 日韩视频在线欧美| 网址你懂的国产日韩在线| 只有这里有精品99| 大片电影免费在线观看免费| 久久精品久久精品一区二区三区| 亚洲国产最新在线播放| 内地一区二区视频在线| 91狼人影院| 亚洲精品视频女| 免费不卡的大黄色大毛片视频在线观看| 国产成人免费无遮挡视频| 麻豆乱淫一区二区| 夜夜爽夜夜爽视频| 国产精品.久久久| 免费观看a级毛片全部| 欧美少妇被猛烈插入视频| 一级毛片我不卡| 纯流量卡能插随身wifi吗| 国产探花极品一区二区| 中文字幕av成人在线电影| 色婷婷av一区二区三区视频| 久久99热这里只频精品6学生| 麻豆精品久久久久久蜜桃| 久热这里只有精品99| 插逼视频在线观看| 久久精品久久久久久噜噜老黄| 青春草亚洲视频在线观看| 熟女人妻精品中文字幕| 亚洲欧美日韩另类电影网站 | 亚洲中文av在线| 91aial.com中文字幕在线观看| 爱豆传媒免费全集在线观看| 丰满少妇做爰视频| 久久久久久久大尺度免费视频| 久久精品国产亚洲av涩爱| 男女下面进入的视频免费午夜| 国精品久久久久久国模美| 国产精品久久久久久精品古装| 国产一区二区三区综合在线观看 | 多毛熟女@视频| av不卡在线播放| 亚洲国产毛片av蜜桃av| 国产一区二区三区综合在线观看 | 高清日韩中文字幕在线| 嫩草影院入口| av线在线观看网站| 日韩,欧美,国产一区二区三区| 亚洲熟女精品中文字幕| 黑人猛操日本美女一级片| 狂野欧美激情性bbbbbb| 小蜜桃在线观看免费完整版高清| 国产av码专区亚洲av| 亚洲国产av新网站| 亚洲婷婷狠狠爱综合网| 精品亚洲乱码少妇综合久久| 女人久久www免费人成看片| 深夜a级毛片| 美女xxoo啪啪120秒动态图| 久久久久久久大尺度免费视频| 在线精品无人区一区二区三 | 国产中年淑女户外野战色| 久久久精品免费免费高清| 亚洲精品日韩在线中文字幕| av一本久久久久| 久久久久性生活片| 中文资源天堂在线| av在线观看视频网站免费| 国产成人精品婷婷| 水蜜桃什么品种好| 国产在线一区二区三区精| 国产成人精品福利久久| 日本色播在线视频| 日韩成人av中文字幕在线观看| 一区二区三区乱码不卡18| 亚洲国产精品国产精品| 岛国毛片在线播放| 欧美日韩在线观看h| 精品一区在线观看国产| 亚洲国产色片| 欧美一级a爱片免费观看看| 免费观看在线日韩| 日韩大片免费观看网站| 性色avwww在线观看| 又黄又爽又刺激的免费视频.| 九九在线视频观看精品| 亚洲国产av新网站| 在线观看免费高清a一片| 欧美极品一区二区三区四区| 亚洲在久久综合| 欧美日韩在线观看h| 日本av手机在线免费观看| 久久精品久久久久久久性| 午夜免费男女啪啪视频观看| 亚洲av国产av综合av卡| 激情 狠狠 欧美| 欧美xxxx性猛交bbbb| 免费人成在线观看视频色| 成人毛片a级毛片在线播放| 色哟哟·www| 在线免费观看不下载黄p国产| 欧美日本视频| 蜜臀久久99精品久久宅男| 丰满迷人的少妇在线观看| 韩国高清视频一区二区三区| 丰满少妇做爰视频| 91精品国产九色| 五月开心婷婷网| 成人国产av品久久久| 人人妻人人澡人人爽人人夜夜| 精华霜和精华液先用哪个| 夫妻午夜视频| 久久久久国产精品人妻一区二区| 97超碰精品成人国产| av在线观看视频网站免费| 99九九线精品视频在线观看视频| 国产一级毛片在线| 十八禁网站网址无遮挡 | 亚洲伊人久久精品综合| 午夜老司机福利剧场| 成人综合一区亚洲| 男女啪啪激烈高潮av片| 国产高清国产精品国产三级 | 欧美精品一区二区免费开放| 久久久久久久精品精品| 亚洲精品乱码久久久v下载方式| 肉色欧美久久久久久久蜜桃| 高清av免费在线| 欧美丝袜亚洲另类| 黄色配什么色好看| 少妇人妻 视频| 国产在视频线精品| 国产在线男女| 国产精品偷伦视频观看了| 一级二级三级毛片免费看| 成年av动漫网址| 身体一侧抽搐| av一本久久久久| 欧美成人精品欧美一级黄| 波野结衣二区三区在线| 国产成人91sexporn| 亚洲精品日本国产第一区| 91狼人影院| 亚洲怡红院男人天堂| 性色avwww在线观看| 国产v大片淫在线免费观看| 2022亚洲国产成人精品| 插阴视频在线观看视频| 如何舔出高潮| 我的老师免费观看完整版| 我要看黄色一级片免费的| 久久久久国产精品人妻一区二区| 国产亚洲av片在线观看秒播厂| 永久免费av网站大全| 在线观看国产h片| 中文字幕亚洲精品专区| 国产高清不卡午夜福利| 久久久久视频综合| 干丝袜人妻中文字幕| 看十八女毛片水多多多| 一本色道久久久久久精品综合| 最近的中文字幕免费完整| 亚洲国产精品999| 中文欧美无线码| 亚洲欧美日韩卡通动漫| 老师上课跳d突然被开到最大视频| 夜夜看夜夜爽夜夜摸| 欧美少妇被猛烈插入视频| 免费看av在线观看网站| 国产探花极品一区二区| 99久久中文字幕三级久久日本| 乱系列少妇在线播放| av天堂中文字幕网| 亚洲av成人精品一二三区| 日韩亚洲欧美综合| 日韩成人伦理影院| 久久久成人免费电影| av不卡在线播放| 亚洲,欧美,日韩| 精品久久久久久电影网| 免费观看的影片在线观看| 国产精品国产av在线观看| 最新中文字幕久久久久| 久久久久网色| 麻豆成人av视频| 高清日韩中文字幕在线| 精品一区二区三区视频在线| 国产精品不卡视频一区二区| 菩萨蛮人人尽说江南好唐韦庄| 久久久久国产精品人妻一区二区| 岛国毛片在线播放| av又黄又爽大尺度在线免费看| 啦啦啦啦在线视频资源| av.在线天堂| 夫妻午夜视频| 成年av动漫网址| 在线天堂最新版资源| 午夜激情久久久久久久| 老熟女久久久| 久久久久久九九精品二区国产| 久久99热这里只频精品6学生| 观看免费一级毛片| 亚洲真实伦在线观看| 久久6这里有精品| 人妻 亚洲 视频| 精品久久久久久久久亚洲| videossex国产| 成年女人在线观看亚洲视频| 日本-黄色视频高清免费观看| 国国产精品蜜臀av免费| 久久精品国产亚洲av天美| 中文字幕av成人在线电影| 国产精品精品国产色婷婷| 尤物成人国产欧美一区二区三区| 欧美日韩一区二区视频在线观看视频在线| 色婷婷av一区二区三区视频| 日韩 亚洲 欧美在线| 日韩国内少妇激情av| 国产一区二区在线观看日韩| 男女边吃奶边做爰视频| 韩国高清视频一区二区三区| 久久精品熟女亚洲av麻豆精品| 国产爱豆传媒在线观看| 国产精品精品国产色婷婷| 久久精品国产亚洲av天美| 国内少妇人妻偷人精品xxx网站| 亚洲精品久久午夜乱码| 亚洲国产毛片av蜜桃av| 国产真实伦视频高清在线观看| 性色avwww在线观看| 80岁老熟妇乱子伦牲交| 久久精品久久久久久噜噜老黄| 国产精品熟女久久久久浪| 黄色一级大片看看| 高清av免费在线| 1000部很黄的大片| 久久久久人妻精品一区果冻| 亚洲美女搞黄在线观看| 日产精品乱码卡一卡2卡三| 久久6这里有精品| 国产在线男女| 中国三级夫妇交换| 美女视频免费永久观看网站| 亚洲电影在线观看av| 精品少妇久久久久久888优播| 欧美一级a爱片免费观看看| 视频区图区小说| 亚洲精品中文字幕在线视频 | 国产色爽女视频免费观看| 久久99蜜桃精品久久| 久久久久网色| 99精国产麻豆久久婷婷| 舔av片在线| 99热这里只有是精品50| 精品人妻视频免费看| 日韩 亚洲 欧美在线| 肉色欧美久久久久久久蜜桃| 国产av国产精品国产| av视频免费观看在线观看| 97超视频在线观看视频| 美女视频免费永久观看网站| 在线观看av片永久免费下载| 大片免费播放器 马上看| 99国产精品免费福利视频| 欧美少妇被猛烈插入视频| 中文在线观看免费www的网站| 亚洲熟女精品中文字幕| 最近2019中文字幕mv第一页| av天堂中文字幕网| 免费在线观看成人毛片| 纯流量卡能插随身wifi吗| 亚洲av中文av极速乱| 日韩免费高清中文字幕av| 日韩,欧美,国产一区二区三区| 一级毛片电影观看| 中文字幕人妻熟人妻熟丝袜美| 哪个播放器可以免费观看大片| 丝袜脚勾引网站| 插逼视频在线观看| 91精品国产国语对白视频| kizo精华| 高清av免费在线| 乱码一卡2卡4卡精品| 免费播放大片免费观看视频在线观看| 在线免费观看不下载黄p国产| 久久人妻熟女aⅴ| 最近中文字幕高清免费大全6| kizo精华| 一边亲一边摸免费视频| 极品少妇高潮喷水抽搐| 国产一区二区三区av在线| 在线免费观看不下载黄p国产| 国产成人精品久久久久久| 亚洲高清免费不卡视频| 亚洲精品一二三| 一级a做视频免费观看| 国产乱人偷精品视频| av免费观看日本| 久久久久国产精品人妻一区二区| 1000部很黄的大片| 97热精品久久久久久| 欧美激情极品国产一区二区三区 | 亚洲va在线va天堂va国产| 日韩,欧美,国产一区二区三区| 男女国产视频网站| 青春草视频在线免费观看| 日韩中文字幕视频在线看片 | 亚洲成人中文字幕在线播放| 97精品久久久久久久久久精品| 高清黄色对白视频在线免费看 | av不卡在线播放| 大陆偷拍与自拍| 久久99热这里只频精品6学生| 日日摸夜夜添夜夜添av毛片| 男女国产视频网站| 亚洲av中文字字幕乱码综合| 一级片'在线观看视频| 老师上课跳d突然被开到最大视频| 国内揄拍国产精品人妻在线| 美女视频免费永久观看网站| 亚洲欧美日韩无卡精品| 亚洲欧美日韩另类电影网站 | 黑人高潮一二区| 精品人妻一区二区三区麻豆| 久久精品久久久久久噜噜老黄| 欧美 日韩 精品 国产| 日本av免费视频播放| 国产亚洲一区二区精品| 自拍偷自拍亚洲精品老妇| 亚洲精品成人av观看孕妇| 91午夜精品亚洲一区二区三区| 热99国产精品久久久久久7| 欧美日韩一区二区视频在线观看视频在线| av国产久精品久网站免费入址| 成人高潮视频无遮挡免费网站| 黄色怎么调成土黄色| 久久亚洲国产成人精品v| 亚洲精品国产成人久久av| 少妇的逼水好多| 亚州av有码| 精品人妻一区二区三区麻豆| 精品久久久久久电影网| 久久久久久人妻| 亚洲精品,欧美精品| 亚洲精品日本国产第一区| 你懂的网址亚洲精品在线观看| 美女内射精品一级片tv| 国产在线视频一区二区| 一级毛片久久久久久久久女| 日韩欧美一区视频在线观看 | 99久久中文字幕三级久久日本| 内地一区二区视频在线| 久久久久精品久久久久真实原创| 国产精品久久久久久久久免| 自拍偷自拍亚洲精品老妇| 韩国av在线不卡| 免费人妻精品一区二区三区视频| 99热全是精品| 日韩av不卡免费在线播放| 欧美精品亚洲一区二区| 亚洲精品亚洲一区二区| 久久精品夜色国产| 日韩av在线免费看完整版不卡| 狂野欧美白嫩少妇大欣赏| 国产极品天堂在线| 99热这里只有是精品在线观看| 亚洲真实伦在线观看| 亚洲欧美日韩另类电影网站 | 男女边吃奶边做爰视频| h视频一区二区三区| xxx大片免费视频| 久久婷婷青草| 久久久精品94久久精品| a级一级毛片免费在线观看| 国产亚洲91精品色在线| 成人综合一区亚洲| 亚洲成人av在线免费| 婷婷色综合www| 永久免费av网站大全| 欧美精品亚洲一区二区| 亚洲精品aⅴ在线观看| 日韩 亚洲 欧美在线| 日本与韩国留学比较| 男人狂女人下面高潮的视频| av免费观看日本| 精品99又大又爽又粗少妇毛片| 亚洲熟女精品中文字幕| 国产91av在线免费观看| 男女下面进入的视频免费午夜| 国产在线一区二区三区精| 日本vs欧美在线观看视频 | 国产成人精品一,二区| 中国三级夫妇交换| 国产色爽女视频免费观看| 久久人人爽人人爽人人片va| 另类亚洲欧美激情| 久久ye,这里只有精品| 欧美老熟妇乱子伦牲交| 色视频在线一区二区三区| 国产精品一区二区三区四区免费观看| 男男h啪啪无遮挡| 亚洲成色77777| 大香蕉97超碰在线| 亚洲精品aⅴ在线观看| 少妇 在线观看| 国产午夜精品久久久久久一区二区三区| 极品教师在线视频| 亚洲成人中文字幕在线播放| 亚洲av在线观看美女高潮| 亚州av有码| 亚洲精品一二三| 日本爱情动作片www.在线观看| 欧美最新免费一区二区三区| 日韩欧美精品免费久久| 成人毛片a级毛片在线播放| 久久99热6这里只有精品| 看非洲黑人一级黄片| 乱码一卡2卡4卡精品| 国产精品免费大片| 亚洲欧美成人精品一区二区| 少妇裸体淫交视频免费看高清| 国产成人精品福利久久| 大香蕉97超碰在线| 亚洲av国产av综合av卡|