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

    航天器姿態(tài)穩(wěn)定的自適應(yīng)有限時(shí)間控制

    2016-11-03 02:38:59高計(jì)委蔡遠(yuǎn)利
    固體火箭技術(shù) 2016年2期
    關(guān)鍵詞:控制精度航天器滑模

    高計(jì)委,蔡遠(yuǎn)利

    (西安交通大學(xué) 控制工程所,西安 710049)

    ?

    航天器姿態(tài)穩(wěn)定的自適應(yīng)有限時(shí)間控制

    高計(jì)委,蔡遠(yuǎn)利

    (西安交通大學(xué) 控制工程所,西安710049)

    針對(duì)剛體航天器姿態(tài)控制問(wèn)題,設(shè)計(jì)了有限時(shí)間狀態(tài)反饋控制律??紤]存在慣性不確性和外部擾動(dòng)條件下,采用四元數(shù)將姿態(tài)系統(tǒng)模型轉(zhuǎn)化為L(zhǎng)agrange模型,并且首次將其以狀態(tài)空間形式表示。其次,基于齊次性理論來(lái)設(shè)計(jì)標(biāo)稱(chēng)控制器,在此基礎(chǔ)上利用自適應(yīng)滑模方法抑制總不確定進(jìn)而使姿態(tài)系統(tǒng)能夠在有限時(shí)間內(nèi)收斂到原點(diǎn),并且利用Lyapunov理論對(duì)其進(jìn)行證明。該算法不需要事先獲得總不確定的邊界,對(duì)實(shí)際問(wèn)題具有很好的實(shí)用性和有效性。數(shù)值仿真對(duì)比表明,此方法具有較快的收斂速度、較高的控制精度及較強(qiáng)的魯棒性。

    航天器姿態(tài)控制;有限時(shí)間穩(wěn)定;魯棒性;自適應(yīng)滑模

    0 引言

    航天器的姿態(tài)穩(wěn)定和跟蹤控制是一類(lèi)具有強(qiáng)非線性、強(qiáng)耦合特性的控制問(wèn)題,一直都是航空航天領(lǐng)域的研究熱點(diǎn)與難點(diǎn)。在航天器的控制理論及其應(yīng)用技術(shù)方面,航天器的姿態(tài)穩(wěn)定和跟蹤控制也隨著科學(xué)技術(shù)的發(fā)展而得到密切的關(guān)注和廣泛的研究。近年來(lái),該領(lǐng)域取得許多成果,如PD形式姿態(tài)控制[1]、非線性最優(yōu)化控制[2]、自適應(yīng)控制[3]、Back-stepping控制[4]、滑??刂芠5]等。上述的這些非線性控制方法雖然能保證閉環(huán)控制系統(tǒng)的穩(wěn)定,但是只得到系統(tǒng)在平衡點(diǎn)處的漸近穩(wěn)定。因此,從時(shí)間優(yōu)化的角度來(lái)看,保證控制系統(tǒng)有限時(shí)間收斂的控制方法[6-9]才是時(shí)間最優(yōu)的控制方法。此外,與其他控制方法相比,在有限時(shí)間控制下的閉環(huán)系統(tǒng)還表現(xiàn)出更高的控制精度及更好的抗擾動(dòng)性。由于有限時(shí)間控制算法擁有上述的優(yōu)越性,所以近年來(lái)其在航天器姿態(tài)控制系統(tǒng)中的應(yīng)用也得到了廣泛的重視。

    Jin等[10]較早將有限時(shí)間收斂特性的終端滑模應(yīng)用于姿態(tài)跟蹤控制中,取得較好控制效果。但終端滑模具有奇異性,所以Li等[11]采用非奇異終端滑模對(duì)控制器進(jìn)行設(shè)計(jì),避免在輸入中出現(xiàn)的奇異現(xiàn)象。另外,Hong等[12]基于齊次分離準(zhǔn)則為機(jī)械臂系統(tǒng)設(shè)計(jì)了有限時(shí)間控制器。采用相同方法,Du等[13]利用四元數(shù)將姿態(tài)穩(wěn)定系統(tǒng)轉(zhuǎn)化為L(zhǎng)agrange形式,分別設(shè)計(jì)了局部與全局姿態(tài)有限時(shí)間控制器;而方輝等[14]則針對(duì)修正羅德里格參數(shù)描述的Lagrange模型,提出航天器局部有限時(shí)間控制器。李貴明等[15]則將加冪積分器與齊次理論結(jié)合設(shè)計(jì)姿態(tài)控制器,不但抑制了擾動(dòng)力矩,還使?fàn)顟B(tài)有限時(shí)間收斂到原點(diǎn)的鄰域; Du等[16]利用類(lèi)似方法設(shè)計(jì)全局有限時(shí)間姿態(tài)跟蹤控制器,并且將其拓展到姿態(tài)同步的應(yīng)用中。雖然上述姿態(tài)控制器取得有限時(shí)間特性,但對(duì)于存在慣性不確定與外部擾動(dòng)時(shí),這些控制器的魯棒性能不佳或者需要知道這些不確定的邊界,而這在實(shí)際應(yīng)用中是不現(xiàn)實(shí)的。

    針對(duì)存在慣性不確定及外部干擾的剛體航空航天飛行器,本文提出了一類(lèi)由四元數(shù)作為姿態(tài)描述的飛行器自適應(yīng)有限時(shí)間控制器。首先,將飛行器姿態(tài)動(dòng)力學(xué)與運(yùn)動(dòng)學(xué)模型轉(zhuǎn)化成類(lèi)Lagrange模型,進(jìn)而寫(xiě)成狀態(tài)空間形式。然后,分兩步設(shè)計(jì)控制律:第一步是標(biāo)稱(chēng)控制,它是連續(xù)但不光滑的,而且能在沒(méi)有參數(shù)不確定與外部擾動(dòng)時(shí)有限時(shí)間穩(wěn)定二階姿態(tài)系統(tǒng);第二步為補(bǔ)償控制,采用自適應(yīng)積分滑模方法,補(bǔ)償系統(tǒng)的總不確定,保證控制目標(biāo)的實(shí)現(xiàn)。所設(shè)計(jì)的姿態(tài)穩(wěn)定控制器能有限時(shí)間鎮(zhèn)定系統(tǒng),抑制內(nèi)外擾動(dòng),取得更高的控制精度。另外,不需要事先知道總不確定的邊界,具有很強(qiáng)的實(shí)用意義。

    1 系統(tǒng)與問(wèn)題描述

    1.1有限時(shí)間穩(wěn)定定理

    首先給出有限時(shí)間穩(wěn)定的相關(guān)定義,考慮如下的非線性系統(tǒng):

    (1)

    式中f(x)為連續(xù)函數(shù),且f(0)=0。

    如果系統(tǒng) (1) 是Lyapunov穩(wěn)定的,而且存在一個(gè)時(shí)間函數(shù)T(x),使得對(duì)于所有的t≥T(x),有x(t)=0,則稱(chēng)系統(tǒng) (1) 是有限時(shí)間穩(wěn)定的;如果具有全局性,則稱(chēng)全局有限時(shí)間穩(wěn)定。

    為便于分析,引入齊次函數(shù)的定義,對(duì)于連續(xù)函數(shù)V(x)和(r1,…,rn)∈Rn,ri>0,如果對(duì)于任意的ε>0,存在標(biāo)量k,滿(mǎn)足:

    則稱(chēng)函數(shù)V(x)關(guān)于(r1,…,rn)具有齊次度k。

    令f(x)=(f1(x),…,fn(x))T,若對(duì)于任意的ε>0,存在標(biāo)量k,滿(mǎn)足:

    則稱(chēng)系統(tǒng)(1)關(guān)于(r1,…,rn)具有齊次度k。

    引理1[8]如果系統(tǒng)(1)具有負(fù)的齊次度,即k<0,并且其原點(diǎn)具有漸近穩(wěn)定性,那么系統(tǒng)(1)的原點(diǎn)是有限時(shí)間穩(wěn)定的。

    1.2航天器姿態(tài)系統(tǒng)模型

    飛行器姿態(tài)的運(yùn)動(dòng)學(xué)描述方法有方向余弦、歐拉角、四元數(shù)和修正的羅德里格參數(shù)等。本文采用四元數(shù)來(lái)描述剛體飛行器的姿態(tài)運(yùn)動(dòng)學(xué),其表達(dá)式為

    (2)

    考慮飛行器姿態(tài)動(dòng)力學(xué)方程:

    (3)

    式中J∈R3×3為正定的對(duì)稱(chēng)矩陣,表示飛行器的轉(zhuǎn)動(dòng)慣量矩陣;τ為飛行器的姿態(tài)控制力矩;d為有界的外部干擾向量。

    2 控制系統(tǒng)設(shè)計(jì)

    上述姿態(tài)模型分別采用四元數(shù)與歐拉角速度作為變量來(lái)表示運(yùn)動(dòng)學(xué)與動(dòng)力學(xué)方程,可在此先將其轉(zhuǎn)化為L(zhǎng)agrange模型,再寫(xiě)成狀態(tài)空間的形式,以便于控制器的設(shè)計(jì)。

    2.1基于四元數(shù)的Lagrange模型及狀態(tài)空間

    命題1:根據(jù)假設(shè)q0(t)≠0對(duì)t≥0,姿態(tài)運(yùn)動(dòng)學(xué)與動(dòng)力學(xué)可寫(xiě)成下面Lagrange形式:

    (4)

    其中

    D=MTd+d′

    (5)

    (6)

    其中

    再將J=J0+ΔJ代入式(6)可得式(4),證畢。

    (7)

    其中

    2.2魯棒有限時(shí)間控制器設(shè)計(jì)

    對(duì)式(7)采用如下基本控制器:

    (8)

    其中,w=[w1,w2,w3]T∈R3為實(shí)際控制輸入。w的設(shè)計(jì)分為兩部分,即w=wn+wc(wn為標(biāo)稱(chēng)控制,wc為不確定控制)。標(biāo)稱(chēng)控制設(shè)計(jì)是在無(wú)擾動(dòng)情況下使系統(tǒng)能在有限時(shí)間內(nèi)到達(dá)原點(diǎn),補(bǔ)償控制是在總的不確定擾動(dòng)下也能完成有限時(shí)間控制目標(biāo)。

    定理1考慮如下二階積分鏈系統(tǒng):

    (9)

    對(duì)于任一個(gè)α2∈(0,1),如下?tīng)顟B(tài)反饋控制器能夠使系統(tǒng)(9)在平衡點(diǎn)x=0處有限時(shí)間收斂:

    (10)

    其中,參數(shù)k1,k2>0使得多項(xiàng)式s2+k2s+k1是Hurwitz的,α1,α2滿(mǎn)足α1=α2/(2-α2);Sig(xi)αi=[sig(xi,1)αi,sig(xi,2)αi,sig(xi,3)αi]T(i=1,2),sig(y)v=sign(y)|y|v(v>0)。

    證明:選取一個(gè)Lyapunov函數(shù):

    (11)

    對(duì)上式關(guān)于t求導(dǎo)可得

    由LaSalle不變集定理可知,系統(tǒng)(10)對(duì)于平衡點(diǎn)x=0是漸近穩(wěn)定的。系統(tǒng)(10)對(duì)(2/(1+α1),2/(1+α1),2/(1+α1),1,1,1)的齊次度為α*=α2-1<0,根據(jù)引理1可得,標(biāo)稱(chēng)控制可使系統(tǒng)(10)有限時(shí)間穩(wěn)定。另外,定理1的單輸入單輸出形式在文獻(xiàn)[8]中已有介紹,這里可理解對(duì)其的高維擴(kuò)展。

    上述標(biāo)稱(chēng)控制器(11)在有限時(shí)間內(nèi)使積分鏈系統(tǒng)(10)穩(wěn)定到原點(diǎn),但對(duì)于受內(nèi)外部擾動(dòng)影響的姿態(tài)穩(wěn)定系統(tǒng)(3)或(7)不具有魯棒性。針對(duì)這個(gè)問(wèn)題,本文利用自適應(yīng)滑??刂破髋c上述提出的控制器(11)結(jié)合,在保證姿態(tài)穩(wěn)定系統(tǒng)(7)有限時(shí)間收斂的同時(shí)也使系統(tǒng)具有魯棒性。在設(shè)計(jì)控制器之前,首先定義如下積分滑模流形:

    (12)

    定理2考慮系統(tǒng)(7),若總不確定擾動(dòng)滿(mǎn)足有界假設(shè),則控制律:

    (13)

    可使姿態(tài)系統(tǒng)軌跡(7)在有限時(shí)間內(nèi)穩(wěn)定,其中Sign(s)=[sign(s1),sign(s2),sign(s3)]T。

    注1:由于不連續(xù)控制的原因,閉環(huán)系統(tǒng)會(huì)出現(xiàn)抖動(dòng)現(xiàn)象。這種高頻振動(dòng)會(huì)影響控制性能,嚴(yán)重時(shí)會(huì)導(dǎo)致系統(tǒng)不穩(wěn)定。為了消小抖振,采用飽和函數(shù)代替開(kāi)關(guān)符號(hào)函數(shù):

    Sat(s)=[sat(s1),sat(s2),sat(s3)]

    (14)

    2.3魯棒自適應(yīng)有限時(shí)間控制律設(shè)計(jì)

    控制律(13)可保證姿態(tài)誤差系統(tǒng)(7)的穩(wěn)定性和魯棒性,但增益需要滿(mǎn)足不等式,然而實(shí)際工程中,總的不確定擾動(dòng)邊界是很難獲得,因此下面提出一種自適應(yīng)的補(bǔ)償控制律的設(shè)計(jì)方法,控制增益能夠根據(jù)系統(tǒng)狀態(tài)自動(dòng)調(diào)節(jié),從而不需要事先知道系統(tǒng)不確定的界。

    (15)

    證明:首先選擇一個(gè)Lyapunov函數(shù):

    (16)

    (17)

    (18)

    (19)

    (20)

    其中

    由引理可得系統(tǒng)(6)或(7)能在有限時(shí)間收斂到原點(diǎn)。

    3 仿真與分析

    為了驗(yàn)證提出的有限時(shí)間控制器的有效性,利用MATLAB/SIMULINK軟件,對(duì)姿態(tài)系統(tǒng)進(jìn)行仿真研究。為了對(duì)比,將經(jīng)典PID控制器:

    (21)

    與提出的控制器同時(shí)進(jìn)行仿真。航天器模型和控制器參數(shù)如表1所示,其中不確定慣性量和外部力矩干擾也在考慮當(dāng)中。在仿真中,用d′代替d便于比較不同控制器的控制性能,見(jiàn)表2。

    表1 航天器模型與控制參數(shù)

    表2 不同控制器下仿真對(duì)比結(jié)果

    圖1~圖4中分別描繪四元數(shù)、角速度、滑模變量與輸入控制量隨時(shí)間變化曲線,系統(tǒng)狀態(tài)在約5.2 s以后收斂。通過(guò)圖2~圖4的對(duì)比,可發(fā)現(xiàn)利用自適應(yīng)律和飽和函數(shù)取代開(kāi)關(guān)函數(shù)的方法,都能減少滑模變量的抖振。由圖2、圖3可看出抖振的水平減小一半,并且輸入幅值隨之減小。另外,觀察圖3、圖5可知增益值在滑模變量到達(dá)滑模面后即趨于穩(wěn)定。

    圖1 四元數(shù)響應(yīng)曲線

    (a) 角速度變化曲線

    (b) 滑模曲線

    (c) 轉(zhuǎn)矩輸入

    (a) 角速度變化曲線

    (b) 滑模曲線

    (c) 轉(zhuǎn)矩輸入

    (a) 角速度變化曲線

    (b) 滑模曲線

    (c) 轉(zhuǎn)矩輸入

    圖5 基于定理3和(14)自適應(yīng)增益

    圖6顯示利用經(jīng)典PID控制器對(duì)有相同初值、不同慣性不確定及外界擾動(dòng)的姿態(tài)系統(tǒng)進(jìn)行控制,其穩(wěn)定時(shí)間約為11.3 s。在第一種擾動(dòng)下,設(shè)計(jì)的2種飽和函數(shù)有限時(shí)間控制算法的仿真結(jié)果相似。同樣,由于提出的控制器的參數(shù)沒(méi)有改變,所在第二種擾動(dòng)下的仿真結(jié)果變化也不大。由圖1~圖6中結(jié)果對(duì)比可明顯發(fā)現(xiàn),指數(shù)收斂率緩慢;眾所周知,PID控制器參數(shù)調(diào)節(jié)不易。其次,從圖6可知,當(dāng)初值相同、慣性不確定及外界擾動(dòng)不同時(shí),其控制精度與魯棒性變差,不如提出的控制器 (8)、(10)和(15)那樣具有較強(qiáng)的實(shí)用性。詳細(xì)的對(duì)比結(jié)果在表2中給出,其中收斂時(shí)間指角速度ω進(jìn)入并保持在范圍|ωi|≤0.02(i=1,2,3)內(nèi)所需的時(shí)間,其中魯棒有限時(shí)間控制器在受到2種不確定擾動(dòng)時(shí),基本變化不大,而用飽和函數(shù)代替開(kāi)關(guān)函數(shù)的自適應(yīng)有限時(shí)間控制器的控制精度雖然也發(fā)生變化,但是還能一直保持在較高的范圍內(nèi)。

    (a)第一種情況

    (b) 第二種情況

    4 結(jié)論

    (1)將狀態(tài)空間概念引入以四元數(shù)表示的剛體姿態(tài)穩(wěn)定系統(tǒng)中,有利于各種非線性控制算法的應(yīng)用。

    (2)基于所述轉(zhuǎn)化模型設(shè)計(jì)的魯棒有限時(shí)間控制器使得姿態(tài)穩(wěn)定系統(tǒng)在內(nèi)外部擾動(dòng)下能快速收斂,并且不需要知道擾動(dòng)的邊界,具有很強(qiáng)的實(shí)用性和較高的控制精度。

    (3)用飽和函數(shù)代替符號(hào)函數(shù)可有效減小滑??刂茙?lái)的抖振,提高了控制精度,仿真結(jié)果驗(yàn)證了所提方法的有效性。

    [1]Wen J T Y, Kreutz-Delgado K. The attitude control problem [J]. IEEE Transactions on Automatic Control, 1991, 36(10): 1148-1162.

    [2]Sharma R, Tewari A. Optimal nonlinear tracking of spacecraft attitude maneuvers [J]. IEEE Transactions on Control Systems Technology, 2004, 12(5): 677-682.

    [3]Ahmed J, Coppola V T. Adaptive asymptotic tracking of spacecraft attitude motion with inertia matrix identification [J]. Journal of Guidance, Control, and Dynamics, 1998, 21(5): 684-691.

    [4]孫平, 劉昆. 小型固體運(yùn)載器一級(jí)飛行段姿態(tài)控制方案研究[J]. 固體火箭技術(shù), 2010, 33 (3): 242-246.

    [5]Lu K, Xia Y, Zhu Z, et al. Sliding mode attitude tracking of rigid spacecraft with disturbance [J]. Journal of The Franklin Institutes, 2012, 349(2): 413-440.

    [6]Haimo V T. Finite time controllers [J]. SIAM Journal on Control and Optimization, 1986, 24(4): 760-770.

    [7]Bhat S P, Bernstein D S. Finite-time stability of continu- ous autonomous systems [J]. SIAM Journal on Control Optimization, 2000, 38(3):751-766.

    [8]Bhat S P, Bernstein D S. Finite-time stability of homo- geneous systems [C]//American Control Conference, New Mexico, USA, 1997, 2513-2514.

    [9]Hong Y, Huang J, Xu Y. On an output feedback finite- time stabilization problem [J]. IEEE Transactions on Automatic Control, 2001, 46(2): 305-309.

    [10]Jin E,Sun Z.Robust controllers design with finite time convergence for rigid spacecraft attitude tracking control [J].Aerospace Science and Technology, 2008,12(4):324-330.

    [11]Li S, Wang Z, Fei S. Comments on the paper: Robust controllers design with finite time convergence for rigid spacecraft attitude tracking control [J]. Aerospace Science and Technology, 2011, 15(3): 193-195.

    [12]Hong Y,Xu Y,Huang J.Finite time control for robot manipulators [J].Systems & Control Letters,2002,46(4):243-253.

    [13]Du H, Li S. Finite-time attitude stabilization for a space- craft using homogeneous method [J]. Journal of Guidance, Control and Dynamics, 2012, 35(3): 740-748.

    [14]方輝, 張愛(ài)華, 秦玉華. 剛體航天器姿態(tài)有限時(shí)間穩(wěn)定控制 [J]. 系統(tǒng)工程學(xué)報(bào), 2014, 29 (3): 299-308.

    [15]李貴明, 劉良棟. 剛體衛(wèi)星姿態(tài)的有限時(shí)間控制[J]. 空間控制技術(shù)與應(yīng)用, 2011, 37(3): 1-8.

    [16]Du H, Li S, Qian C. Finite-time attitude tracking control of spacecraft with application to attitude synchronization [J]. IEEE Transactions on Automatic Control, 2011, 56(11): 2711-2717.

    [17]Bhat S P, Bernstein D S. Lyapunov analysis of finite-time differential equations [C]//Proceedings of the American Control Conference. Seattle WA, USA, 1995:1831-1832.

    (編輯:呂耀輝)

    Adaptive finite-time control for spacecraft attitude stabilization

    GAO Ji-wei, CAI Yuan-li

    (Institute of Control Engineering, Xi'an Jiaotong University, Xi'an710049, China)

    In this paper, the finite time feedback control laws were designed for rigid spacecraft attitude control problem. Attitude system model was converted into Lagrange model based on unit quaternion in the presence of inertia uncertainties and external disturbances, and then initially expressed as state-space form. Next, nominal controller was designed based on homogenous theory, and further adaptive sliding mode approach was adopted to steer attitude system to the origin with the lumped uncertainties, while rigorous finite-time stable proofs were achieved in light of the Lyapunov theory. This algorithm requires no information about the lumped uncertainties, and it possesses practicability and effectiveness for realistic systems. Finally, numerical simulation analysis shows that these methodologies can achieve faster convergence, higher control accuracy and stronger robustness.

    spacecraft attitude control;finite-time stability;robustness;adaptive sliding mode

    2014-11-27;

    2015-03-18。

    國(guó)家自然科學(xué)基金(61202128, 61463029);宇航動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室開(kāi)放基金(2011ADL-JD0202)。

    高計(jì)委(1985—),男,博士生,研究方向?yàn)榛?刂萍霸诤教炱髦械膽?yīng)用。E-mail:jwgao2012@163.com

    蔡遠(yuǎn)利(1963—),男,教授/博導(dǎo),研究方向?yàn)轱w行器制導(dǎo)控制與仿真。E-mail:ylicai@mail.xjtu.edu.cn

    V448.2

    A

    1006-2793(2016)02-0275-06

    10.7673/j.issn.1006-2793.2016.02.022

    猜你喜歡
    控制精度航天器滑模
    基于多源異構(gòu)信息融合的采摘機(jī)械臂驅(qū)動(dòng)控制研究
    2022 年第二季度航天器發(fā)射統(tǒng)計(jì)
    2019 年第二季度航天器發(fā)射統(tǒng)計(jì)
    基于組合滑模控制的絕對(duì)重力儀兩級(jí)主動(dòng)減振設(shè)計(jì)
    2018 年第三季度航天器發(fā)射統(tǒng)計(jì)
    測(cè)控技術(shù)(2018年4期)2018-11-25 09:47:26
    并網(wǎng)逆變器逆系統(tǒng)自學(xué)習(xí)滑模抗擾控制
    2018年第二季度航天器發(fā)射統(tǒng)計(jì)
    MW級(jí)太空發(fā)電站微波能量波束指向控制精度分析
    基于安卓的智能車(chē)轉(zhuǎn)速系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)
    51国产日韩欧美| 免费人妻精品一区二区三区视频| 伊人久久国产一区二区| 最近中文字幕高清免费大全6| 久久精品国产亚洲av涩爱| 久久6这里有精品| 亚洲欧美成人精品一区二区| 亚洲欧洲日产国产| 日韩视频在线欧美| 日本猛色少妇xxxxx猛交久久| 黄色一级大片看看| 日韩视频在线欧美| 成年av动漫网址| 99久国产av精品国产电影| 国产精品福利在线免费观看| 性色av一级| 性高湖久久久久久久久免费观看| 爱豆传媒免费全集在线观看| 国产精品蜜桃在线观看| 伦理电影大哥的女人| 亚洲人成网站在线播| 少妇人妻久久综合中文| 免费人成在线观看视频色| 亚洲美女视频黄频| 亚洲av中文字字幕乱码综合| 一区二区三区免费毛片| 亚洲一级一片aⅴ在线观看| 精品一区二区免费观看| 水蜜桃什么品种好| 黄色欧美视频在线观看| 中国国产av一级| av在线app专区| 日韩精品有码人妻一区| 99热这里只有精品一区| 色视频在线一区二区三区| 国产免费一级a男人的天堂| 少妇熟女欧美另类| 日产精品乱码卡一卡2卡三| 国产精品爽爽va在线观看网站| 免费人成在线观看视频色| 日本欧美视频一区| 亚洲熟女精品中文字幕| 欧美区成人在线视频| 久久这里有精品视频免费| 日本黄色片子视频| 18禁裸乳无遮挡免费网站照片| 亚洲欧洲日产国产| 自拍偷自拍亚洲精品老妇| 大又大粗又爽又黄少妇毛片口| 精品久久久久久久末码| 午夜福利高清视频| 国产 一区 欧美 日韩| 久久午夜福利片| 边亲边吃奶的免费视频| 国产人妻一区二区三区在| 欧美精品一区二区大全| 我的老师免费观看完整版| 久久影院123| 又爽又黄a免费视频| 一级a做视频免费观看| 日本av免费视频播放| 最近中文字幕2019免费版| 国产淫语在线视频| 国产成人a区在线观看| 国产爽快片一区二区三区| 国产亚洲一区二区精品| 狂野欧美白嫩少妇大欣赏| 夜夜骑夜夜射夜夜干| 欧美性感艳星| 日韩强制内射视频| 成人黄色视频免费在线看| 亚洲欧洲国产日韩| 亚洲最大成人中文| 乱码一卡2卡4卡精品| 香蕉精品网在线| 黄色视频在线播放观看不卡| 天堂8中文在线网| 舔av片在线| 一二三四中文在线观看免费高清| 一二三四中文在线观看免费高清| 日韩亚洲欧美综合| 成人二区视频| 伊人久久国产一区二区| 成年人午夜在线观看视频| 久久国产亚洲av麻豆专区| 91精品伊人久久大香线蕉| 久久久久久久精品精品| 日韩中文字幕视频在线看片 | 只有这里有精品99| 婷婷色麻豆天堂久久| 欧美三级亚洲精品| 亚洲成人av在线免费| 国产高清三级在线| 亚洲美女黄色视频免费看| 亚洲av日韩在线播放| 亚洲精品乱码久久久久久按摩| 精品久久久久久久久亚洲| 大码成人一级视频| 日韩强制内射视频| 国产精品久久久久久av不卡| 一区二区三区乱码不卡18| 亚洲,一卡二卡三卡| 女人十人毛片免费观看3o分钟| 欧美精品国产亚洲| 国产精品国产三级国产av玫瑰| 国产成人91sexporn| 99热国产这里只有精品6| 国产一区二区三区av在线| 日本免费在线观看一区| 伦理电影大哥的女人| 亚洲精品aⅴ在线观看| 精品人妻视频免费看| 亚洲精品乱久久久久久| 黑人高潮一二区| 激情 狠狠 欧美| 舔av片在线| 18禁裸乳无遮挡动漫免费视频| 黄片无遮挡物在线观看| 久久99精品国语久久久| 久久热精品热| 男女边摸边吃奶| 寂寞人妻少妇视频99o| 国产黄频视频在线观看| 国产欧美亚洲国产| 国产av精品麻豆| av在线蜜桃| 国内揄拍国产精品人妻在线| 国产一区亚洲一区在线观看| 黑人猛操日本美女一级片| xxx大片免费视频| 久久精品国产自在天天线| 夜夜看夜夜爽夜夜摸| 欧美激情极品国产一区二区三区 | 国产成人午夜福利电影在线观看| 亚洲欧美精品专区久久| 日本爱情动作片www.在线观看| 女性被躁到高潮视频| 内地一区二区视频在线| 精品99又大又爽又粗少妇毛片| 直男gayav资源| 久久久久久人妻| 男人添女人高潮全过程视频| 在线观看一区二区三区激情| 蜜桃久久精品国产亚洲av| 国产精品久久久久久精品古装| 老女人水多毛片| 久久精品国产亚洲av涩爱| 高清av免费在线| 欧美老熟妇乱子伦牲交| 亚洲欧美日韩东京热| 一个人看的www免费观看视频| 亚洲av二区三区四区| 波野结衣二区三区在线| 99热国产这里只有精品6| 王馨瑶露胸无遮挡在线观看| 久热这里只有精品99| 亚洲精品视频女| 亚洲国产欧美人成| 国产精品久久久久久av不卡| 97热精品久久久久久| 少妇猛男粗大的猛烈进出视频| 国产成人a∨麻豆精品| 能在线免费看毛片的网站| 在线免费观看不下载黄p国产| 成年免费大片在线观看| 99热这里只有是精品在线观看| 国产精品久久久久久精品古装| 人妻 亚洲 视频| 国产精品一区二区性色av| 免费播放大片免费观看视频在线观看| 久久久久人妻精品一区果冻| 亚洲成人手机| 91精品国产九色| 97超碰精品成人国产| 精品视频人人做人人爽| 午夜免费男女啪啪视频观看| 黄色欧美视频在线观看| 高清视频免费观看一区二区| 女性生殖器流出的白浆| 永久网站在线| 国产黄色免费在线视频| 韩国av在线不卡| 国产中年淑女户外野战色| 乱码一卡2卡4卡精品| 伦精品一区二区三区| 2022亚洲国产成人精品| 欧美国产精品一级二级三级 | 各种免费的搞黄视频| 我要看日韩黄色一级片| 女性被躁到高潮视频| av在线老鸭窝| 国产精品久久久久久久电影| 午夜日本视频在线| 一级毛片 在线播放| 一级二级三级毛片免费看| 老司机影院成人| av在线app专区| 看免费成人av毛片| 欧美亚洲 丝袜 人妻 在线| 亚洲国产欧美人成| 国产片特级美女逼逼视频| 日日摸夜夜添夜夜添av毛片| av又黄又爽大尺度在线免费看| 51国产日韩欧美| 伊人久久精品亚洲午夜| 国产亚洲欧美精品永久| 内地一区二区视频在线| 五月开心婷婷网| av免费观看日本| 黑人高潮一二区| av国产精品久久久久影院| 一本久久精品| 新久久久久国产一级毛片| 最近中文字幕2019免费版| 一区在线观看完整版| 国产精品爽爽va在线观看网站| 国产人妻一区二区三区在| 狠狠精品人妻久久久久久综合| 亚洲国产看品久久| 亚洲av电影在线进入| 国产欧美日韩综合在线一区二区| 欧美日韩视频高清一区二区三区二| 激情视频va一区二区三区| 国产一卡二卡三卡精品| 丝袜人妻中文字幕| 男女高潮啪啪啪动态图| 欧美成人午夜精品| 国产成人精品久久二区二区91| 高潮久久久久久久久久久不卡| 波多野结衣av一区二区av| 日韩大码丰满熟妇| 高清欧美精品videossex| 久热这里只有精品99| 麻豆国产av国片精品| 天堂8中文在线网| 天天躁狠狠躁夜夜躁狠狠躁| 久9热在线精品视频| 色婷婷久久久亚洲欧美| 看免费成人av毛片| 久久久国产一区二区| 久久精品久久久久久久性| 国产高清不卡午夜福利| 久久人人爽av亚洲精品天堂| 欧美黑人精品巨大| 99热网站在线观看| 悠悠久久av| 男女免费视频国产| 国产成人精品久久二区二区免费| 男女床上黄色一级片免费看| www.熟女人妻精品国产| 欧美av亚洲av综合av国产av| 欧美人与性动交α欧美精品济南到| 男人爽女人下面视频在线观看| 久久精品国产亚洲av涩爱| 婷婷丁香在线五月| 欧美+亚洲+日韩+国产| 亚洲激情五月婷婷啪啪| 欧美精品啪啪一区二区三区 | 亚洲第一青青草原| 熟女少妇亚洲综合色aaa.| 国产高清视频在线播放一区 | av片东京热男人的天堂| 黄色 视频免费看| 99国产综合亚洲精品| 欧美久久黑人一区二区| 赤兔流量卡办理| 日韩一卡2卡3卡4卡2021年| 久久人人爽人人片av| 女人久久www免费人成看片| 男女边吃奶边做爰视频| 国产在线观看jvid| 亚洲欧美成人综合另类久久久| 国产精品免费视频内射| 中文字幕人妻丝袜制服| 精品卡一卡二卡四卡免费| 亚洲第一青青草原| 亚洲综合色网址| 丝瓜视频免费看黄片| 丝袜人妻中文字幕| 97人妻天天添夜夜摸| 午夜福利免费观看在线| 校园人妻丝袜中文字幕| 成人午夜精彩视频在线观看| 爱豆传媒免费全集在线观看| 亚洲精品一区蜜桃| 人人妻人人澡人人看| 又紧又爽又黄一区二区| 母亲3免费完整高清在线观看| 9191精品国产免费久久| 性高湖久久久久久久久免费观看| 亚洲精品久久午夜乱码| 国产91精品成人一区二区三区 | 99精品久久久久人妻精品| 欧美日韩精品网址| 51午夜福利影视在线观看| 国产免费视频播放在线视频| 亚洲精品美女久久av网站| www.熟女人妻精品国产| 久久亚洲精品不卡| 精品欧美一区二区三区在线| 丝袜喷水一区| 国产成人免费观看mmmm| 国产在线一区二区三区精| 青草久久国产| 亚洲欧美一区二区三区国产| 亚洲国产精品一区二区三区在线| 中文字幕精品免费在线观看视频| 午夜福利免费观看在线| 国产成人a∨麻豆精品| 亚洲av美国av| www.熟女人妻精品国产| 国产真人三级小视频在线观看| 亚洲国产av影院在线观看| 亚洲精品在线美女| 欧美亚洲 丝袜 人妻 在线| 男人爽女人下面视频在线观看| 黑丝袜美女国产一区| 色视频在线一区二区三区| 另类亚洲欧美激情| 欧美中文综合在线视频| 日韩精品免费视频一区二区三区| 777米奇影视久久| 嫩草影视91久久| 成人黄色视频免费在线看| 秋霞在线观看毛片| 亚洲av国产av综合av卡| 亚洲精品一卡2卡三卡4卡5卡 | 又大又爽又粗| 亚洲精品成人av观看孕妇| 国产免费福利视频在线观看| 美国免费a级毛片| 在线天堂中文资源库| 精品第一国产精品| 免费在线观看完整版高清| 久久久国产一区二区| 午夜福利视频在线观看免费| bbb黄色大片| 黄色 视频免费看| 欧美日韩成人在线一区二区| 建设人人有责人人尽责人人享有的| av网站免费在线观看视频| 最近手机中文字幕大全| 狠狠婷婷综合久久久久久88av| 色播在线永久视频| 成人亚洲精品一区在线观看| 国产1区2区3区精品| 两个人免费观看高清视频| 91麻豆av在线| 国产亚洲av片在线观看秒播厂| 黄色片一级片一级黄色片| 妹子高潮喷水视频| 久久九九热精品免费| 欧美成人午夜精品| 亚洲九九香蕉| 国产视频首页在线观看| 新久久久久国产一级毛片| 日本av免费视频播放| 亚洲一区中文字幕在线| 美女高潮到喷水免费观看| 欧美中文综合在线视频| 亚洲av成人不卡在线观看播放网 | 精品福利永久在线观看| 91老司机精品| 天天影视国产精品| 久久精品亚洲av国产电影网| 精品高清国产在线一区| 宅男免费午夜| 七月丁香在线播放| 亚洲精品一卡2卡三卡4卡5卡 | 成人午夜精彩视频在线观看| 国产成人91sexporn| 精品国产超薄肉色丝袜足j| 免费观看av网站的网址| 啦啦啦中文免费视频观看日本| 两个人免费观看高清视频| 国产97色在线日韩免费| 精品少妇一区二区三区视频日本电影| 成在线人永久免费视频| 一区二区三区四区激情视频| 欧美97在线视频| 国产国语露脸激情在线看| 精品卡一卡二卡四卡免费| 欧美黑人精品巨大| 精品一区二区三区四区五区乱码 | 久久毛片免费看一区二区三区| bbb黄色大片| 亚洲欧美激情在线| 欧美人与性动交α欧美软件| 少妇被粗大的猛进出69影院| 黄色片一级片一级黄色片| 亚洲欧美精品综合一区二区三区| videosex国产| 亚洲情色 制服丝袜| 国产精品久久久av美女十八| 亚洲美女黄色视频免费看| 亚洲一区中文字幕在线| 夫妻性生交免费视频一级片| 国产一区二区 视频在线| 丝袜喷水一区| 婷婷色综合大香蕉| 女性被躁到高潮视频| 国产成人精品在线电影| 制服人妻中文乱码| 国产成人系列免费观看| 国产黄色免费在线视频| 好男人视频免费观看在线| 久热这里只有精品99| 亚洲国产av新网站| 妹子高潮喷水视频| 免费女性裸体啪啪无遮挡网站| 天堂8中文在线网| 男人添女人高潮全过程视频| 最黄视频免费看| 少妇 在线观看| 人体艺术视频欧美日本| 黄色怎么调成土黄色| 国产精品麻豆人妻色哟哟久久| 伊人亚洲综合成人网| 大片免费播放器 马上看| 午夜福利免费观看在线| 人妻人人澡人人爽人人| 精品免费久久久久久久清纯 | 在线观看一区二区三区激情| 日本vs欧美在线观看视频| 国产亚洲一区二区精品| 青草久久国产| 18禁观看日本| 99久久精品国产亚洲精品| 欧美av亚洲av综合av国产av| 国产成人啪精品午夜网站| 成人午夜精彩视频在线观看| 国产精品一区二区免费欧美 | 99香蕉大伊视频| 亚洲 国产 在线| 国产成人91sexporn| 久久久国产欧美日韩av| 自拍欧美九色日韩亚洲蝌蚪91| 美女午夜性视频免费| 2018国产大陆天天弄谢| 热99久久久久精品小说推荐| 人人澡人人妻人| 黄色一级大片看看| 日韩免费高清中文字幕av| 国产成人精品久久久久久| 国产片特级美女逼逼视频| 久久免费观看电影| 婷婷成人精品国产| 老司机午夜十八禁免费视频| 欧美av亚洲av综合av国产av| 久久久精品94久久精品| 深夜精品福利| 亚洲国产欧美网| 成年美女黄网站色视频大全免费| 狂野欧美激情性bbbbbb| 青草久久国产| 久久精品国产a三级三级三级| 操美女的视频在线观看| 久久狼人影院| 中国国产av一级| 免费在线观看日本一区| 国产精品九九99| 咕卡用的链子| 欧美精品人与动牲交sv欧美| 亚洲成人免费电影在线观看 | 少妇精品久久久久久久| 国产精品三级大全| 91精品三级在线观看| 99久久综合免费| 亚洲欧美一区二区三区黑人| 国产精品免费视频内射| 婷婷色综合大香蕉| 一个人免费看片子| 亚洲欧洲精品一区二区精品久久久| 精品少妇黑人巨大在线播放| 亚洲国产精品一区三区| 精品人妻在线不人妻| 中文字幕色久视频| 纯流量卡能插随身wifi吗| 亚洲色图 男人天堂 中文字幕| 欧美日韩av久久| 日本a在线网址| 好男人视频免费观看在线| 久久精品国产亚洲av涩爱| av天堂久久9| 精品人妻在线不人妻| 人人澡人人妻人| 久久久国产精品麻豆| 又黄又粗又硬又大视频| a级毛片在线看网站| 999精品在线视频| 性少妇av在线| 考比视频在线观看| 中文字幕av电影在线播放| 99国产综合亚洲精品| 精品福利观看| 如日韩欧美国产精品一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 高清av免费在线| 久久亚洲国产成人精品v| 日日摸夜夜添夜夜爱| 欧美日韩亚洲综合一区二区三区_| 欧美日韩精品网址| 99香蕉大伊视频| 波多野结衣av一区二区av| 国产片内射在线| 欧美激情极品国产一区二区三区| 不卡av一区二区三区| 欧美日韩成人在线一区二区| 亚洲精品一区蜜桃| 一级毛片电影观看| 女人久久www免费人成看片| 99国产精品一区二区三区| 亚洲精品日韩在线中文字幕| 久久久久久人人人人人| 国产成人欧美在线观看 | 久久国产精品人妻蜜桃| 亚洲av日韩精品久久久久久密 | 人人妻人人澡人人爽人人夜夜| av不卡在线播放| 免费看十八禁软件| av线在线观看网站| svipshipincom国产片| 久久精品国产亚洲av涩爱| 亚洲精品国产av成人精品| 又大又爽又粗| 精品亚洲乱码少妇综合久久| 丰满人妻熟妇乱又伦精品不卡| 91精品三级在线观看| 又黄又粗又硬又大视频| 91国产中文字幕| 91成人精品电影| 国产精品欧美亚洲77777| 欧美激情极品国产一区二区三区| 亚洲精品日本国产第一区| 欧美日韩视频精品一区| 国产精品久久久av美女十八| 黄色a级毛片大全视频| 少妇裸体淫交视频免费看高清 | 日韩中文字幕欧美一区二区 | 中文精品一卡2卡3卡4更新| 人妻一区二区av| 国产在视频线精品| 18在线观看网站| 亚洲精品一区蜜桃| 欧美精品av麻豆av| 成在线人永久免费视频| 麻豆乱淫一区二区| 日韩一卡2卡3卡4卡2021年| 国产精品久久久人人做人人爽| 午夜福利视频精品| 一级,二级,三级黄色视频| 欧美+亚洲+日韩+国产| 多毛熟女@视频| 久久久国产精品麻豆| 亚洲欧美一区二区三区国产| 别揉我奶头~嗯~啊~动态视频 | 乱人伦中国视频| e午夜精品久久久久久久| 香蕉国产在线看| 晚上一个人看的免费电影| 99精品久久久久人妻精品| 亚洲精品国产av蜜桃| 一区二区av电影网| 男人舔女人的私密视频| 午夜老司机福利片| 在线观看人妻少妇| 亚洲精品乱久久久久久| 美女午夜性视频免费| 最近中文字幕2019免费版| 91九色精品人成在线观看| 免费观看人在逋| 女性被躁到高潮视频| 精品福利永久在线观看| 一级黄色大片毛片| 亚洲熟女精品中文字幕| 国产深夜福利视频在线观看| 一本色道久久久久久精品综合| 精品亚洲成国产av| 自拍欧美九色日韩亚洲蝌蚪91| 中文字幕制服av| 欧美在线一区亚洲| 亚洲专区国产一区二区| 亚洲五月婷婷丁香| 少妇粗大呻吟视频| 免费在线观看黄色视频的| 91字幕亚洲| 国产又爽黄色视频| 国产亚洲欧美精品永久| av视频免费观看在线观看| 国产视频首页在线观看| 亚洲久久久国产精品| 免费一级毛片在线播放高清视频 | 男人舔女人的私密视频| 国产麻豆69| 男人添女人高潮全过程视频| 亚洲精品久久久久久婷婷小说| 亚洲精品一卡2卡三卡4卡5卡 | 欧美激情 高清一区二区三区| 午夜免费鲁丝| 色94色欧美一区二区| 国产片特级美女逼逼视频| bbb黄色大片| 宅男免费午夜| 国产精品熟女久久久久浪| 又大又黄又爽视频免费| 久久毛片免费看一区二区三区| 美国免费a级毛片| 最新的欧美精品一区二区| 少妇粗大呻吟视频| 成人国产av品久久久| 亚洲五月婷婷丁香| 欧美日韩亚洲高清精品| 国产成人91sexporn| 日韩电影二区| 伦理电影免费视频|