• <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男天堂| 国产高清有码在线观看视频| 国产日韩欧美视频二区| 亚洲第一av免费看| 亚洲av欧美aⅴ国产| 久久99一区二区三区| 日韩av在线免费看完整版不卡| 国产一区亚洲一区在线观看| 久久国内精品自在自线图片| 久久久国产欧美日韩av| 日本-黄色视频高清免费观看| 久久精品国产自在天天线| 在线观看国产h片| 国产精品一区www在线观看| 在线观看一区二区三区激情| 人妻夜夜爽99麻豆av| 一本—道久久a久久精品蜜桃钙片| 丝袜美足系列| 乱码一卡2卡4卡精品| 十八禁网站网址无遮挡| 午夜福利,免费看| 99视频精品全部免费 在线| 成人黄色视频免费在线看| 日韩av免费高清视频| 国产精品久久久久成人av| 亚洲精品一区蜜桃| 肉色欧美久久久久久久蜜桃| 高清av免费在线| 国产精品久久久久成人av| 亚洲美女搞黄在线观看| av视频免费观看在线观看| 国产女主播在线喷水免费视频网站| 亚洲精品乱久久久久久| 五月天丁香电影| 亚洲av成人精品一区久久| 九九在线视频观看精品| 亚洲av福利一区| 大香蕉97超碰在线| 久久婷婷青草| 国产极品天堂在线| 亚洲精品久久午夜乱码| 欧美成人午夜免费资源| 99国产综合亚洲精品| 精品一区二区免费观看| 高清欧美精品videossex| 91成人精品电影| 国产高清不卡午夜福利| 美女国产视频在线观看| av电影中文网址| 丝瓜视频免费看黄片| 在线观看免费视频网站a站| 考比视频在线观看| 成人免费观看视频高清| 日本wwww免费看| 高清不卡的av网站| 乱人伦中国视频| 免费高清在线观看日韩| 午夜激情久久久久久久| a 毛片基地| 亚洲精品久久成人aⅴ小说 | 久久久久久久精品精品| 午夜精品国产一区二区电影| 大陆偷拍与自拍| 亚洲美女视频黄频| 五月天丁香电影| 亚洲精品久久成人aⅴ小说 | 国产黄色视频一区二区在线观看| 午夜福利,免费看| 亚洲美女黄色视频免费看| 综合色丁香网| av在线app专区| 高清午夜精品一区二区三区| 国产精品三级大全| 999精品在线视频| 伊人久久精品亚洲午夜| 一区二区av电影网| 大香蕉97超碰在线| 午夜免费男女啪啪视频观看| 婷婷色综合大香蕉| 男男h啪啪无遮挡| 亚洲精品乱码久久久v下载方式| 国产一区二区三区综合在线观看 | 久久韩国三级中文字幕| 国产一区有黄有色的免费视频| av.在线天堂| 在线观看三级黄色| 亚洲欧美日韩另类电影网站| 欧美+日韩+精品| 观看av在线不卡| 一二三四中文在线观看免费高清| 久久久久久久亚洲中文字幕| 2022亚洲国产成人精品| 高清不卡的av网站| 欧美97在线视频| 22中文网久久字幕| 亚洲精品国产色婷婷电影| a级毛色黄片| 国产在线视频一区二区| 日日撸夜夜添| 少妇人妻精品综合一区二区| 精品酒店卫生间| 国产一区二区在线观看av| 夜夜骑夜夜射夜夜干| 精品一区二区三卡| 亚洲天堂av无毛| 在线观看免费高清a一片| 高清在线视频一区二区三区| 久久热精品热| 9色porny在线观看| 成年av动漫网址| 欧美日韩综合久久久久久| 街头女战士在线观看网站| 久久久国产一区二区| 欧美日韩av久久| a级片在线免费高清观看视频| 91aial.com中文字幕在线观看| 国产午夜精品久久久久久一区二区三区| 成人影院久久| 久久精品国产a三级三级三级| 在线天堂最新版资源| av福利片在线| 国产免费一区二区三区四区乱码| 欧美变态另类bdsm刘玥| freevideosex欧美| 国产精品欧美亚洲77777| 丰满迷人的少妇在线观看| 91国产中文字幕| 亚洲精品国产色婷婷电影| 最近的中文字幕免费完整| 欧美一级a爱片免费观看看| 国产成人精品久久久久久| 中文精品一卡2卡3卡4更新| 国产免费视频播放在线视频| 国产视频内射| kizo精华| a级毛片在线看网站| av一本久久久久| 最近的中文字幕免费完整| 国产精品人妻久久久久久| 一区二区三区精品91| 欧美性感艳星| 纯流量卡能插随身wifi吗| 九草在线视频观看| 伊人亚洲综合成人网| 亚洲综合色网址| 久久 成人 亚洲| 18禁动态无遮挡网站| 中国美白少妇内射xxxbb| 日产精品乱码卡一卡2卡三| 国产精品一区二区在线不卡| 午夜激情福利司机影院| 欧美变态另类bdsm刘玥| 一级毛片黄色毛片免费观看视频| 精品一区二区免费观看| 精品卡一卡二卡四卡免费| av线在线观看网站| 免费高清在线观看日韩| 亚洲国产av影院在线观看| 亚洲色图综合在线观看| 蜜桃久久精品国产亚洲av| 国产精品国产av在线观看| 国产有黄有色有爽视频| 欧美三级亚洲精品| 婷婷色综合大香蕉| 只有这里有精品99| 又大又黄又爽视频免费| 国产亚洲最大av| 精品久久久精品久久久| 三级国产精品片| 男人添女人高潮全过程视频| 亚洲人成77777在线视频| 观看av在线不卡| 国国产精品蜜臀av免费| 国产亚洲精品久久久com| 亚洲av不卡在线观看| 麻豆成人av视频| 亚洲成色77777| 成人手机av| 精品久久久久久久久av| 高清欧美精品videossex| 丁香六月天网| 欧美最新免费一区二区三区| 人人妻人人澡人人看| 王馨瑶露胸无遮挡在线观看| 极品少妇高潮喷水抽搐| 十八禁网站网址无遮挡| 男女边摸边吃奶| 国产一区亚洲一区在线观看| 亚洲婷婷狠狠爱综合网| 狂野欧美白嫩少妇大欣赏| 日韩欧美一区视频在线观看| 人体艺术视频欧美日本| 观看av在线不卡| 搡老乐熟女国产| 欧美日韩av久久| 国产成人精品一,二区| 亚洲人成77777在线视频| 看非洲黑人一级黄片| 国产综合精华液| 亚洲人与动物交配视频| 亚洲第一av免费看| 国产精品一区二区在线观看99| 亚洲欧美成人综合另类久久久| 观看美女的网站| 日本黄色片子视频| 国产成人午夜福利电影在线观看| 久久免费观看电影| 亚洲欧洲日产国产| a 毛片基地| 午夜免费观看性视频| 亚洲av欧美aⅴ国产| 777米奇影视久久| 国产精品久久久久久av不卡| 国产精品99久久99久久久不卡 | 99热全是精品| 国产色婷婷99| 大陆偷拍与自拍| 3wmmmm亚洲av在线观看| 久久久久人妻精品一区果冻| 欧美精品亚洲一区二区| 亚洲av成人精品一区久久| 免费高清在线观看日韩| 日韩精品有码人妻一区| 亚洲国产精品一区三区| 亚洲欧美日韩另类电影网站| 国产日韩欧美在线精品| av.在线天堂| 亚洲国产av新网站| 国产成人免费无遮挡视频| 国产一级毛片在线| 欧美激情极品国产一区二区三区 | 好男人视频免费观看在线| 国产色婷婷99| av免费在线看不卡| 少妇人妻精品综合一区二区| 国产成人精品无人区| 国产乱来视频区| 国产精品不卡视频一区二区| 国产精品三级大全| 国产免费又黄又爽又色| 亚洲国产av影院在线观看| 久久99热6这里只有精品| 肉色欧美久久久久久久蜜桃| 天美传媒精品一区二区| 亚洲精华国产精华液的使用体验| 久久久久久久久大av| 成人二区视频| 日产精品乱码卡一卡2卡三| 国产在视频线精品| 免费av不卡在线播放| 国产男女内射视频| 亚洲精品456在线播放app| 午夜福利网站1000一区二区三区| 婷婷色综合www| 伊人久久精品亚洲午夜| 国产精品女同一区二区软件| 久久久精品免费免费高清| 亚洲综合色惰| 在线亚洲精品国产二区图片欧美 | 久久久久久久国产电影| 午夜日本视频在线| a 毛片基地| 国产av码专区亚洲av| 综合色丁香网| 中国美白少妇内射xxxbb| 欧美精品国产亚洲| 国产一区有黄有色的免费视频| 3wmmmm亚洲av在线观看| 秋霞在线观看毛片| 男女啪啪激烈高潮av片| 下体分泌物呈黄色| 国产女主播在线喷水免费视频网站| 大片免费播放器 马上看| 欧美+日韩+精品| 国国产精品蜜臀av免费| h视频一区二区三区| 久久久久久伊人网av| 欧美 日韩 精品 国产| 亚洲av中文av极速乱| 欧美 亚洲 国产 日韩一| 亚洲经典国产精华液单| 高清不卡的av网站| 久久鲁丝午夜福利片| 丰满乱子伦码专区| av专区在线播放| 97超视频在线观看视频| 精品亚洲成a人片在线观看| 欧美人与性动交α欧美精品济南到 | 欧美变态另类bdsm刘玥| 亚洲欧洲日产国产| 秋霞伦理黄片| 国产成人精品婷婷| 中文字幕久久专区| 亚洲国产欧美在线一区| 亚洲第一av免费看| av又黄又爽大尺度在线免费看| 亚洲不卡免费看| 国产精品无大码| 少妇人妻精品综合一区二区| 国产精品 国内视频| 大香蕉97超碰在线| 日本-黄色视频高清免费观看| 精品少妇内射三级| 五月开心婷婷网| 天美传媒精品一区二区| 一个人看视频在线观看www免费| 91久久精品国产一区二区三区| 国产成人精品一,二区| 亚洲伊人久久精品综合| 蜜桃久久精品国产亚洲av| 日本色播在线视频| 欧美丝袜亚洲另类| 日本av免费视频播放| 亚洲欧美日韩另类电影网站| 亚洲精品日韩在线中文字幕| 最近的中文字幕免费完整| 久久99一区二区三区| 一本大道久久a久久精品| 成人亚洲精品一区在线观看| 在线观看免费高清a一片| 十八禁网站网址无遮挡| 91在线精品国自产拍蜜月| 国产毛片在线视频| 青青草视频在线视频观看| 国产不卡av网站在线观看| 男女边吃奶边做爰视频| 少妇被粗大的猛进出69影院 | 黑人欧美特级aaaaaa片| 蜜桃在线观看..| 亚洲欧洲日产国产| 另类亚洲欧美激情| 丰满乱子伦码专区| 亚洲精品一二三| √禁漫天堂资源中文www| 青春草视频在线免费观看| 18禁裸乳无遮挡动漫免费视频| 乱人伦中国视频| 日本午夜av视频| 在线 av 中文字幕| 国产一区有黄有色的免费视频| 内地一区二区视频在线| 狠狠婷婷综合久久久久久88av| 国产亚洲欧美精品永久| 美女视频免费永久观看网站| 欧美三级亚洲精品| 人人妻人人添人人爽欧美一区卜| 一个人免费看片子| 欧美日韩精品成人综合77777| 久久精品夜色国产| .国产精品久久| 一级二级三级毛片免费看| 国产色婷婷99| 国产伦精品一区二区三区视频9| 99热这里只有是精品在线观看| 99国产综合亚洲精品| 日韩强制内射视频| 老熟女久久久| 亚洲av日韩在线播放| 色哟哟·www| 国产69精品久久久久777片| 亚洲性久久影院| 国产精品秋霞免费鲁丝片| 女性被躁到高潮视频| 尾随美女入室| 搡女人真爽免费视频火全软件| a级毛片黄视频| 久久久久久久久久成人| 亚洲国产av影院在线观看| 国产精品一区www在线观看| 欧美激情极品国产一区二区三区 | 一区二区av电影网| 在现免费观看毛片| 国产免费又黄又爽又色| 精品久久久久久久久亚洲| 国产亚洲一区二区精品| 欧美人与性动交α欧美精品济南到 | 成人国产麻豆网| 熟妇人妻不卡中文字幕| 美女国产视频在线观看| 天天影视国产精品| 男女边吃奶边做爰视频| 日韩成人伦理影院| 欧美少妇被猛烈插入视频| 热99久久久久精品小说推荐| 视频中文字幕在线观看| 看十八女毛片水多多多| 一区二区三区乱码不卡18| 国产av码专区亚洲av| 婷婷色综合www| 日本-黄色视频高清免费观看| 老司机亚洲免费影院| 日本黄大片高清| 亚洲欧美中文字幕日韩二区| 亚洲精品日韩在线中文字幕| 麻豆精品久久久久久蜜桃| 国产午夜精品久久久久久一区二区三区| 亚洲av综合色区一区| 成人黄色视频免费在线看| 中文字幕制服av| 日韩不卡一区二区三区视频在线| 久久久国产欧美日韩av| 国产精品三级大全| 亚洲,欧美,日韩| av国产精品久久久久影院| 国产在线免费精品| 人妻人人澡人人爽人人| 久久久久久久久久人人人人人人| 一本久久精品| 母亲3免费完整高清在线观看 | 久久ye,这里只有精品| 国产毛片在线视频| 久久韩国三级中文字幕| 久久人妻熟女aⅴ| 久久久久网色| 2022亚洲国产成人精品| 飞空精品影院首页| 国产国语露脸激情在线看| 国产在线一区二区三区精| 最近最新中文字幕免费大全7| 国产日韩欧美在线精品| 成人影院久久| 黄色一级大片看看| 日韩免费高清中文字幕av| 亚洲av日韩在线播放| 亚洲欧美一区二区三区国产| 午夜av观看不卡| 内地一区二区视频在线| 国产欧美另类精品又又久久亚洲欧美| 欧美变态另类bdsm刘玥| av国产精品久久久久影院| 伦理电影大哥的女人| 在线观看美女被高潮喷水网站| 日本91视频免费播放| 美女主播在线视频| 99热这里只有是精品在线观看| 一区在线观看完整版| 丝瓜视频免费看黄片| 免费不卡的大黄色大毛片视频在线观看| 日韩视频在线欧美| 久久影院123| 免费大片黄手机在线观看| 成人无遮挡网站| 寂寞人妻少妇视频99o| 亚洲欧美一区二区三区国产| 免费观看a级毛片全部| 99久久人妻综合| 日韩中文字幕视频在线看片| 热99久久久久精品小说推荐| 日日摸夜夜添夜夜爱| 99视频精品全部免费 在线| 老司机亚洲免费影院| 婷婷色麻豆天堂久久| 国产一级毛片在线| 我的女老师完整版在线观看| 欧美变态另类bdsm刘玥| 免费观看无遮挡的男女| 亚洲欧美精品自产自拍| 国产视频内射| 极品人妻少妇av视频| 亚洲三级黄色毛片| 黄片播放在线免费| 99精国产麻豆久久婷婷| 水蜜桃什么品种好| 欧美日韩精品成人综合77777| 涩涩av久久男人的天堂| 欧美另类一区| 黄色配什么色好看| 在现免费观看毛片| 黄片无遮挡物在线观看| 亚洲精品一区蜜桃| 欧美国产精品一级二级三级| 久久鲁丝午夜福利片| 免费观看性生交大片5| kizo精华| 丰满迷人的少妇在线观看| 欧美变态另类bdsm刘玥| 国产精品国产三级国产av玫瑰| 香蕉精品网在线| 18在线观看网站| 国产 精品1| 满18在线观看网站| 王馨瑶露胸无遮挡在线观看| 久久久久久久大尺度免费视频| 久久99一区二区三区| 我要看黄色一级片免费的| 插逼视频在线观看| av国产久精品久网站免费入址| 亚洲中文av在线| 一二三四中文在线观看免费高清| 亚洲在久久综合| 9色porny在线观看| kizo精华| 亚洲国产精品一区三区| 91精品一卡2卡3卡4卡| 亚洲欧洲精品一区二区精品久久久 | videos熟女内射| 国产精品人妻久久久影院| 在线播放无遮挡| 亚洲婷婷狠狠爱综合网| 91午夜精品亚洲一区二区三区| 亚洲av不卡在线观看| 美女国产视频在线观看| 国产亚洲精品第一综合不卡 | a级毛片免费高清观看在线播放| 亚洲人成网站在线观看播放| 久久影院123| 久久久久久久久久成人| 在线观看美女被高潮喷水网站| 五月开心婷婷网| videossex国产| 肉色欧美久久久久久久蜜桃| 曰老女人黄片| 视频在线观看一区二区三区| 少妇人妻精品综合一区二区| 视频在线观看一区二区三区| 秋霞在线观看毛片| 2022亚洲国产成人精品| 人妻系列 视频| 视频在线观看一区二区三区| 在线免费观看不下载黄p国产| 亚洲精品日韩在线中文字幕| 99热全是精品| 国产精品成人在线| 在线观看www视频免费| 制服诱惑二区| 丝袜在线中文字幕| 婷婷色综合大香蕉| 亚洲图色成人| av免费在线看不卡| 国产日韩欧美在线精品| 一级片'在线观看视频| 亚洲精品久久成人aⅴ小说 | 免费黄色在线免费观看| 只有这里有精品99| 999精品在线视频| 成人国产麻豆网| 免费人成在线观看视频色| 啦啦啦中文免费视频观看日本| 寂寞人妻少妇视频99o| 91精品国产九色| 午夜激情久久久久久久| 国产精品一区二区三区四区免费观看| 国产精品久久久久久精品电影小说| 3wmmmm亚洲av在线观看| 99久久精品一区二区三区| 91精品伊人久久大香线蕉| 久久久亚洲精品成人影院| 亚洲图色成人| 在线观看国产h片| 高清午夜精品一区二区三区| 亚洲成人一二三区av| 亚洲av福利一区| 欧美激情极品国产一区二区三区 | 少妇的逼水好多| 一区二区三区免费毛片| 国产黄频视频在线观看| 黄色欧美视频在线观看| 久久国产亚洲av麻豆专区| 日韩av在线免费看完整版不卡| 亚洲国产精品国产精品| 亚洲av电影在线观看一区二区三区| 国产av码专区亚洲av| 热re99久久精品国产66热6| 久久久久久伊人网av| 国产成人一区二区在线| 国产老妇伦熟女老妇高清| 又大又黄又爽视频免费| 久久毛片免费看一区二区三区| 午夜老司机福利剧场| 91aial.com中文字幕在线观看| 黄色欧美视频在线观看| 91久久精品国产一区二区成人| 亚洲精品视频女| 老女人水多毛片| 免费观看性生交大片5| 国产永久视频网站| 极品少妇高潮喷水抽搐| 黄色怎么调成土黄色| 国产女主播在线喷水免费视频网站| 熟妇人妻不卡中文字幕| 搡老乐熟女国产| 精品熟女少妇av免费看| 狠狠精品人妻久久久久久综合| 亚洲欧洲日产国产| 国产av精品麻豆| 美女视频免费永久观看网站| 夫妻午夜视频| 国产有黄有色有爽视频| 如何舔出高潮| 精品一区在线观看国产| 插逼视频在线观看| 特大巨黑吊av在线直播| 国产亚洲最大av| 免费观看av网站的网址| 中文字幕制服av| 国产高清国产精品国产三级| 精品熟女少妇av免费看| 熟妇人妻不卡中文字幕| 美女xxoo啪啪120秒动态图| 永久网站在线| 99热这里只有精品一区| 久久午夜综合久久蜜桃| 人妻夜夜爽99麻豆av| 国产在线一区二区三区精| 午夜免费鲁丝| 精品人妻偷拍中文字幕| av在线app专区| 国产av国产精品国产| 2022亚洲国产成人精品| 丰满少妇做爰视频| 亚洲精品,欧美精品| 久久精品熟女亚洲av麻豆精品| 国产免费一区二区三区四区乱码| 亚洲色图 男人天堂 中文字幕 |