• <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)
    精品福利观看| 亚洲国产高清在线一区二区三| 人人妻,人人澡人人爽秒播| 国产91精品成人一区二区三区| 18禁黄网站禁片免费观看直播| 在线天堂最新版资源| 国产精品电影一区二区三区| 夜夜夜夜夜久久久久| 桃红色精品国产亚洲av| 成人特级av手机在线观看| 18禁在线播放成人免费| 可以在线观看的亚洲视频| 亚洲人与动物交配视频| 欧美成人性av电影在线观看| 亚洲成人精品中文字幕电影| 欧美日韩国产亚洲二区| 国产精品永久免费网站| 女警被强在线播放| 十八禁人妻一区二区| 久久精品影院6| 天天躁日日操中文字幕| 婷婷丁香在线五月| 成人鲁丝片一二三区免费| 免费搜索国产男女视频| 精品乱码久久久久久99久播| 美女大奶头视频| 手机成人av网站| 99久久成人亚洲精品观看| 精品日产1卡2卡| 欧美大码av| 国产在视频线在精品| 老鸭窝网址在线观看| 亚洲av美国av| 国产三级在线视频| 免费在线观看影片大全网站| 亚洲一区二区三区不卡视频| 国产av在哪里看| 少妇人妻一区二区三区视频| 一区二区三区激情视频| 欧洲精品卡2卡3卡4卡5卡区| 亚洲人与动物交配视频| 嫁个100分男人电影在线观看| 国产国拍精品亚洲av在线观看 | 中文字幕人妻丝袜一区二区| 18禁国产床啪视频网站| 国产蜜桃级精品一区二区三区| 国产真实伦视频高清在线观看 | 亚洲av成人av| 在线观看av片永久免费下载| 国内精品美女久久久久久| 免费观看精品视频网站| 欧美日韩福利视频一区二区| 男女之事视频高清在线观看| 男女下面进入的视频免费午夜| 99国产极品粉嫩在线观看| 日韩欧美国产一区二区入口| 黄色成人免费大全| 99热精品在线国产| 99热这里只有是精品50| 免费在线观看成人毛片| 精品熟女少妇八av免费久了| 精品欧美国产一区二区三| 一进一出抽搐动态| 欧美xxxx黑人xx丫x性爽| 天堂动漫精品| 少妇人妻一区二区三区视频| bbb黄色大片| 国产亚洲av嫩草精品影院| 国产精品三级大全| 亚洲国产欧洲综合997久久,| 黄色视频,在线免费观看| 老鸭窝网址在线观看| 在线观看日韩欧美| 十八禁人妻一区二区| 国产一区二区三区在线臀色熟女| 国产欧美日韩精品亚洲av| 久久亚洲精品不卡| 欧美+亚洲+日韩+国产| ponron亚洲| 成人三级黄色视频| 精品一区二区三区人妻视频| 亚洲精品一卡2卡三卡4卡5卡| 女同久久另类99精品国产91| 国产成+人综合+亚洲专区| 搡老熟女国产l中国老女人| 人人妻人人看人人澡| 夜夜看夜夜爽夜夜摸| 日韩欧美三级三区| 黄色片一级片一级黄色片| 一区二区三区国产精品乱码| 日本熟妇午夜| 99国产精品一区二区三区| 国产伦精品一区二区三区四那| 久久九九热精品免费| 免费人成视频x8x8入口观看| 制服丝袜大香蕉在线| 偷拍熟女少妇极品色| 三级毛片av免费| 色综合婷婷激情| 99视频精品全部免费 在线| 12—13女人毛片做爰片一| 伊人久久精品亚洲午夜| 国内精品一区二区在线观看| 亚洲午夜理论影院| 亚洲欧美日韩高清在线视频| 人妻丰满熟妇av一区二区三区| 中文在线观看免费www的网站| 99久久综合精品五月天人人| 国产精品久久久久久久久免 | 男女之事视频高清在线观看| а√天堂www在线а√下载| 免费看美女性在线毛片视频| 嫁个100分男人电影在线观看| 午夜精品在线福利| 禁无遮挡网站| 在线观看免费午夜福利视频| 女人被狂操c到高潮| 91在线观看av| 午夜福利在线观看免费完整高清在 | 啦啦啦观看免费观看视频高清| xxx96com| 丁香六月欧美| 国产午夜精品论理片| 亚洲美女黄片视频| 狠狠狠狠99中文字幕| 国产精品 欧美亚洲| 国产亚洲精品久久久久久毛片| 国产高清视频在线观看网站| 国产精品 国内视频| 亚洲人成伊人成综合网2020| 亚洲欧美日韩高清专用| 久99久视频精品免费| 精品欧美国产一区二区三| 男插女下体视频免费在线播放| 国产欧美日韩一区二区精品| 亚洲成av人片免费观看| 国产亚洲精品综合一区在线观看| 国产野战对白在线观看| 男人的好看免费观看在线视频| 精品国内亚洲2022精品成人| 国产真人三级小视频在线观看| 日日干狠狠操夜夜爽| 成年人黄色毛片网站| 操出白浆在线播放| 国产精品久久久久久久久免 | 国产精品久久久久久精品电影| 三级国产精品欧美在线观看| 狂野欧美激情性xxxx| 黄色丝袜av网址大全| 真实男女啪啪啪动态图| 色老头精品视频在线观看| 一本精品99久久精品77| 香蕉丝袜av| 在线观看免费午夜福利视频| 欧美日韩亚洲国产一区二区在线观看| 在线播放国产精品三级| 国产熟女xx| 成人av一区二区三区在线看| 午夜a级毛片| 麻豆国产97在线/欧美| av片东京热男人的天堂| 久久精品91蜜桃| 国产亚洲欧美在线一区二区| 中出人妻视频一区二区| 在线十欧美十亚洲十日本专区| 日韩人妻高清精品专区| 国产不卡一卡二| 有码 亚洲区| 老鸭窝网址在线观看| 99久久精品国产亚洲精品| 亚洲av成人精品一区久久| 成人高潮视频无遮挡免费网站| 夜夜看夜夜爽夜夜摸| 神马国产精品三级电影在线观看| www日本黄色视频网| 一二三四社区在线视频社区8| 国产亚洲欧美在线一区二区| 夜夜爽天天搞| 在线看三级毛片| 国产激情偷乱视频一区二区| 日韩 欧美 亚洲 中文字幕| 欧美乱妇无乱码| 在线观看舔阴道视频| 久久亚洲精品不卡| 亚洲欧美日韩卡通动漫| 性欧美人与动物交配| 天堂影院成人在线观看| 午夜免费男女啪啪视频观看 | 动漫黄色视频在线观看| 偷拍熟女少妇极品色| 国产久久久一区二区三区| 人妻夜夜爽99麻豆av| 51午夜福利影视在线观看| 国产97色在线日韩免费| 麻豆一二三区av精品| 九九在线视频观看精品| 国产真实伦视频高清在线观看 | 夜夜夜夜夜久久久久| 色av中文字幕| 国产伦一二天堂av在线观看| 亚洲av免费在线观看| 99久久精品国产亚洲精品| 亚洲专区国产一区二区| 亚洲精品粉嫩美女一区| bbb黄色大片| 啦啦啦观看免费观看视频高清| 又黄又爽又免费观看的视频| 亚洲精品在线美女| 国产高潮美女av| 他把我摸到了高潮在线观看| 国产精品久久视频播放| 国产精品亚洲美女久久久| 91麻豆精品激情在线观看国产| aaaaa片日本免费| 久久国产精品影院| 我要搜黄色片| 少妇的逼水好多| 久久久成人免费电影| 久9热在线精品视频| 亚洲中文日韩欧美视频| 每晚都被弄得嗷嗷叫到高潮| 香蕉av资源在线| 久久久久久大精品| 久久性视频一级片| 99久久综合精品五月天人人| 一本综合久久免费| 亚洲熟妇熟女久久| 国产在视频线在精品| 2021天堂中文幕一二区在线观| 亚洲精品久久国产高清桃花| 岛国在线观看网站| 久久精品综合一区二区三区| 日本一本二区三区精品| 久久99热这里只有精品18| 1024手机看黄色片| 9191精品国产免费久久| 人人妻,人人澡人人爽秒播| 美女高潮喷水抽搐中文字幕| 国产黄色小视频在线观看| 一本一本综合久久| 99热这里只有是精品50| 国产午夜福利久久久久久| 亚洲片人在线观看| 亚洲一区高清亚洲精品| 在线观看av片永久免费下载| 亚洲av免费高清在线观看| 在线十欧美十亚洲十日本专区| 国产成人福利小说| 亚洲人成网站高清观看| 亚洲国产精品sss在线观看| 免费观看人在逋| 一级黄片播放器| 成人欧美大片| 网址你懂的国产日韩在线| 12—13女人毛片做爰片一| 国产精品一区二区三区四区久久| 一本一本综合久久| 婷婷亚洲欧美| 亚洲av一区综合| 免费人成在线观看视频色| 男女做爰动态图高潮gif福利片| 可以在线观看的亚洲视频| 久久国产精品影院| 国内毛片毛片毛片毛片毛片| 精品午夜福利视频在线观看一区| 亚洲无线观看免费| 欧美一区二区精品小视频在线| 国产不卡一卡二| 露出奶头的视频| 国产成人av教育| 国产熟女xx| 亚洲国产欧洲综合997久久,| 真实男女啪啪啪动态图| 亚洲激情在线av| 亚洲第一欧美日韩一区二区三区| av天堂在线播放| 亚洲欧美日韩高清专用| 国产熟女xx| 亚洲精品成人久久久久久| 欧美黑人巨大hd| 在线天堂最新版资源| 99视频精品全部免费 在线| 国产高清激情床上av| 欧美黑人巨大hd| 天堂影院成人在线观看| 国产主播在线观看一区二区| 99久国产av精品| 三级国产精品欧美在线观看| 久久九九热精品免费| 国产乱人伦免费视频| 午夜亚洲福利在线播放| 国产午夜精品论理片| 18禁国产床啪视频网站| 美女 人体艺术 gogo| 国产精品亚洲美女久久久| 久久人人精品亚洲av| 亚洲狠狠婷婷综合久久图片| 一本久久中文字幕| 国产高清videossex| 99久久九九国产精品国产免费| 美女大奶头视频| 国产中年淑女户外野战色| 老汉色∧v一级毛片| 国模一区二区三区四区视频| 免费无遮挡裸体视频| 中文字幕熟女人妻在线| 人妻丰满熟妇av一区二区三区| 无遮挡黄片免费观看| 欧美成人a在线观看| 在线播放国产精品三级| 日本黄色视频三级网站网址| 国产精品,欧美在线| 成人无遮挡网站| 两个人看的免费小视频| 国产高潮美女av| 51午夜福利影视在线观看| 日韩欧美在线乱码| or卡值多少钱| 亚洲国产日韩欧美精品在线观看 | 免费在线观看日本一区| 长腿黑丝高跟| 日本三级黄在线观看| 日本免费一区二区三区高清不卡| 免费看美女性在线毛片视频| 成人鲁丝片一二三区免费| 国内精品一区二区在线观看| 欧美又色又爽又黄视频| 99在线人妻在线中文字幕| 热99re8久久精品国产| 美女被艹到高潮喷水动态| 91字幕亚洲| 国产精品女同一区二区软件 | 午夜激情福利司机影院| 午夜福利高清视频| 老鸭窝网址在线观看| 丁香六月欧美| 国产高清视频在线播放一区| 欧美性猛交黑人性爽| 午夜激情欧美在线| 免费av观看视频| 人人妻,人人澡人人爽秒播| 成人特级av手机在线观看| bbb黄色大片| 免费在线观看日本一区| 91久久精品电影网| 国产av不卡久久| 在线看三级毛片| 白带黄色成豆腐渣| 最近最新免费中文字幕在线| 成年版毛片免费区| 天天躁日日操中文字幕| 女同久久另类99精品国产91| 九九热线精品视视频播放| 人妻久久中文字幕网| 国产高清视频在线播放一区| 国产一区二区三区在线臀色熟女| 性欧美人与动物交配| 波野结衣二区三区在线 | 51午夜福利影视在线观看| 老司机午夜十八禁免费视频| 成人永久免费在线观看视频| 黄色女人牲交| 99精品久久久久人妻精品| 好看av亚洲va欧美ⅴa在| 国产黄a三级三级三级人| 校园春色视频在线观看| 一区二区三区高清视频在线| АⅤ资源中文在线天堂| 亚洲最大成人手机在线| 亚洲不卡免费看| 欧美日韩黄片免| 成人特级av手机在线观看| 99久久综合精品五月天人人| 激情在线观看视频在线高清| 免费av毛片视频| 男女视频在线观看网站免费| 国产在线精品亚洲第一网站| 1000部很黄的大片| 国产高清视频在线观看网站| 中文字幕精品亚洲无线码一区| 精品久久久久久久人妻蜜臀av| 无遮挡黄片免费观看| 国产精品 国内视频| 老司机深夜福利视频在线观看| 女生性感内裤真人,穿戴方法视频| 色尼玛亚洲综合影院| 国产精品影院久久| 久久亚洲真实| 真实男女啪啪啪动态图| 老熟妇仑乱视频hdxx| av在线天堂中文字幕| 欧美成人免费av一区二区三区| 亚洲专区国产一区二区| 99精品欧美一区二区三区四区| 夜夜夜夜夜久久久久| 麻豆成人av在线观看| 麻豆国产97在线/欧美| 婷婷丁香在线五月| 久久精品国产亚洲av香蕉五月| 久久精品国产99精品国产亚洲性色| 老熟妇仑乱视频hdxx| 色播亚洲综合网| 亚洲av成人精品一区久久| 91字幕亚洲| av在线天堂中文字幕| xxx96com| 亚洲成人精品中文字幕电影| 久久久国产成人免费| 久久国产精品影院| 1024手机看黄色片| 精品人妻1区二区| 亚洲精品456在线播放app | 国产99白浆流出| 国产真实乱freesex| 97碰自拍视频| 成人av在线播放网站| 白带黄色成豆腐渣| 国产精品一区二区三区四区久久| 日本撒尿小便嘘嘘汇集6| 黄色日韩在线| 99在线人妻在线中文字幕| 国产v大片淫在线免费观看| 欧美一区二区精品小视频在线| 天美传媒精品一区二区| 淫妇啪啪啪对白视频| 亚洲成人精品中文字幕电影| tocl精华| 黄色日韩在线| 国产精品嫩草影院av在线观看 | 高清在线国产一区| АⅤ资源中文在线天堂| 韩国av一区二区三区四区| 欧美性猛交╳xxx乱大交人| 亚洲在线自拍视频| 国产精品久久久人人做人人爽| 别揉我奶头~嗯~啊~动态视频| 无限看片的www在线观看| 国产伦精品一区二区三区视频9 | 成熟少妇高潮喷水视频| 男人舔奶头视频| 午夜精品在线福利| 日韩免费av在线播放| 国产日本99.免费观看| 欧美在线黄色| 亚洲久久久久久中文字幕| 精品乱码久久久久久99久播| 一区二区三区国产精品乱码| 一二三四社区在线视频社区8| 一本精品99久久精品77| 国产真实伦视频高清在线观看 | 国产老妇女一区| 波多野结衣巨乳人妻| 丰满人妻一区二区三区视频av | 午夜a级毛片| 久久人妻av系列| 亚洲国产欧美人成| 免费看光身美女| 久久久久国内视频| 日本免费一区二区三区高清不卡| 最近在线观看免费完整版| 亚洲精品在线美女| 日本 欧美在线| 亚洲午夜理论影院| 波多野结衣高清作品| 美女 人体艺术 gogo| 国产精品国产高清国产av| 亚洲内射少妇av| 国产 一区 欧美 日韩| 精品人妻1区二区| 亚洲成av人片免费观看| 亚洲国产精品合色在线| 麻豆国产97在线/欧美| 天堂网av新在线| 亚洲av免费在线观看| 中文字幕av在线有码专区| 国产高清有码在线观看视频| 国产精品亚洲一级av第二区| 国产成人av教育| 国产真人三级小视频在线观看| 国产 一区 欧美 日韩| 欧美日韩中文字幕国产精品一区二区三区| 成人高潮视频无遮挡免费网站| 高清毛片免费观看视频网站| 久久精品人妻少妇| 亚洲avbb在线观看| 九色成人免费人妻av| 久久精品亚洲精品国产色婷小说| 久久天躁狠狠躁夜夜2o2o| 精品福利观看| 嫩草影院入口| 老司机在亚洲福利影院| 69人妻影院| 搡老熟女国产l中国老女人| 淫秽高清视频在线观看| 丁香六月欧美| АⅤ资源中文在线天堂| 国产aⅴ精品一区二区三区波| 日本黄大片高清| 亚洲精品久久国产高清桃花| 一个人看视频在线观看www免费 | 久久久国产成人精品二区| 久久欧美精品欧美久久欧美| 成年人黄色毛片网站| 国产精品99久久99久久久不卡| 一夜夜www| 久久久久久久久中文| 91麻豆精品激情在线观看国产| 亚洲欧美精品综合久久99| 亚洲国产精品sss在线观看| 亚洲人与动物交配视频| 亚洲精品在线美女| 女人十人毛片免费观看3o分钟| 一区二区三区国产精品乱码| 一二三四社区在线视频社区8| 亚洲久久久久久中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 亚洲欧美精品综合久久99| 亚洲av美国av| 熟妇人妻久久中文字幕3abv| 亚洲五月婷婷丁香| АⅤ资源中文在线天堂| 性色av乱码一区二区三区2| 成人欧美大片| 天堂√8在线中文| 看免费av毛片| 亚洲狠狠婷婷综合久久图片| 在线观看66精品国产| 非洲黑人性xxxx精品又粗又长| 国产精品av视频在线免费观看| 国产激情偷乱视频一区二区| 免费人成视频x8x8入口观看| ponron亚洲| 人妻久久中文字幕网| 99在线人妻在线中文字幕| 好看av亚洲va欧美ⅴa在| 午夜亚洲福利在线播放| 国内精品一区二区在线观看| 午夜亚洲福利在线播放| a级一级毛片免费在线观看| АⅤ资源中文在线天堂| 亚洲天堂国产精品一区在线| 国产真实伦视频高清在线观看 | 国产精品1区2区在线观看.| 男女之事视频高清在线观看| 国产亚洲精品一区二区www| 国产av不卡久久| 怎么达到女性高潮| 性色av乱码一区二区三区2| 不卡一级毛片| 精品日产1卡2卡| 国产成人福利小说| 一个人看的www免费观看视频| 无人区码免费观看不卡| 久久99热这里只有精品18| 国产亚洲精品综合一区在线观看| www.www免费av| 午夜福利18| 高潮久久久久久久久久久不卡| 中文字幕熟女人妻在线| 亚洲欧美日韩东京热| 99久久无色码亚洲精品果冻| 757午夜福利合集在线观看| 最近最新中文字幕大全免费视频| 国产亚洲精品一区二区www| 国产精品女同一区二区软件 | a级毛片a级免费在线| 韩国av一区二区三区四区| 一边摸一边抽搐一进一小说| 在线观看av片永久免费下载| 亚洲男人的天堂狠狠| 90打野战视频偷拍视频| 日韩中文字幕欧美一区二区| 蜜桃亚洲精品一区二区三区| 亚洲精品粉嫩美女一区| 一级黄色大片毛片| 女人被狂操c到高潮| 国产国拍精品亚洲av在线观看 | 国产黄色小视频在线观看| av在线蜜桃| 一卡2卡三卡四卡精品乱码亚洲| 有码 亚洲区| 麻豆成人av在线观看| 首页视频小说图片口味搜索| 变态另类丝袜制服| 三级男女做爰猛烈吃奶摸视频| 成年免费大片在线观看| www.999成人在线观看| 老司机在亚洲福利影院| 欧美色欧美亚洲另类二区| 99精品欧美一区二区三区四区| 成人国产综合亚洲| 好男人在线观看高清免费视频| 嫩草影视91久久| 免费看日本二区| 五月伊人婷婷丁香| 少妇的丰满在线观看| 国产伦精品一区二区三区视频9 | 手机成人av网站| 美女cb高潮喷水在线观看| 一区二区三区高清视频在线| 亚洲人与动物交配视频| 91字幕亚洲| 精品午夜福利视频在线观看一区| 欧美av亚洲av综合av国产av| 欧美成人免费av一区二区三区| 国产高清激情床上av| 九色国产91popny在线| 97超视频在线观看视频| 美女被艹到高潮喷水动态| 偷拍熟女少妇极品色| 亚洲专区国产一区二区| 一本一本综合久久| 男女下面进入的视频免费午夜| 韩国av一区二区三区四区| 久久久国产成人免费|