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

    新型高超聲速飛行器耦合姿態(tài)控制系統(tǒng)設(shè)計

    2017-03-31 02:21:08郭宗易郭建國
    宇航學(xué)報 2017年3期
    關(guān)鍵詞:超聲速觀測器飛行器

    郭宗易,周 軍,郭建國

    (西北工業(yè)大學(xué)精確制導(dǎo)與控制研究所,西安710072)

    新型高超聲速飛行器耦合姿態(tài)控制系統(tǒng)設(shè)計

    郭宗易,周 軍,郭建國

    (西北工業(yè)大學(xué)精確制導(dǎo)與控制研究所,西安710072)

    針對具有強耦合和強不確定性特點的高超聲速飛行器控制問題,提出一種考慮耦合特性的新型魯棒姿態(tài)控制律。首先,建立高超聲速飛行器的姿態(tài)動力學(xué)方程,并基于反步法思想轉(zhuǎn)化為誤差動態(tài)模型;然后結(jié)合干擾觀測器提出基于耦合特性的飛行器魯棒控制方法,并設(shè)計滑模補償項保證姿態(tài)角漸近收斂到跟蹤指令;最后基于絕對誤差積分(IAE)指標(biāo)從理論上證明本文方法相對于未考慮耦合特性的方法具有更好的動態(tài)性能。本文提出的控制方法不僅保證系統(tǒng)漸近穩(wěn)定性和強魯棒性,而且將耦合特性考慮進控制器設(shè)計中,進而提升了系統(tǒng)的動態(tài)性能。數(shù)值仿真考慮氣動參數(shù)拉偏以及實際的執(zhí)行機構(gòu)幅值和速率限制,仿真結(jié)果校驗了本文控制系統(tǒng)設(shè)計的有效性。

    高超聲速飛行器;耦合特性;干擾觀測器;滑模補償項;絕對誤差積分(IAE)

    0 引 言

    高超聲速飛行器以其快速響應(yīng)能力、高突防能力和軍事戰(zhàn)略意義成為近些年各國的研究熱點[1]。高超聲速飛行器的高速、高空域飛行特點,導(dǎo)致其存在強不確定性、強非線性和強耦合等控制難題,尤其耦合問題目前仍是亟待研究和解決的關(guān)鍵。

    解決耦合問題的思路之一是進行解耦控制。由于飛行器模型的復(fù)雜性,常規(guī)的線性解耦方法已經(jīng)難以應(yīng)用到飛行器控制器設(shè)計中,所以近些年常用的方法主要是微分幾何動態(tài)逆控制策略,其基本思想是利用微分幾何方法實現(xiàn)系統(tǒng)的精確反饋線性化,實現(xiàn)解耦目的,再結(jié)合各種方法設(shè)計動態(tài)逆控制器。文獻[2]針對再入飛行器進行反饋線性化,得到其解耦模型,然后設(shè)計滑模干擾觀測器對不確定性進行估計和補償;文獻[3]利用動態(tài)逆與μ綜合結(jié)合的方法設(shè)計飛行器縱向通道魯棒控制器;文獻[4]采用反饋線性化得到解耦模型,然后基于H∞方法設(shè)計控制器。然而,動態(tài)逆對模型精確度十分敏感,飛行器的非線性動態(tài)特性非常復(fù)雜,難以建立精確的數(shù)學(xué)模型,設(shè)計的控制器魯棒性也較差。

    另一種處理耦合的思路是在設(shè)計控制器時考慮耦合因素。文獻[5]模型中考慮了飛行器姿態(tài)模型的通道耦合,在進行解耦后定義殘余耦合量,將之考慮進系統(tǒng)的不確定性并采用變結(jié)構(gòu)控制保證系統(tǒng)穩(wěn)定;文獻[6]將高超聲速飛行器耦合系統(tǒng)表示成具有非匹配不確定性的關(guān)聯(lián)大系統(tǒng)形式,基于廣義關(guān)聯(lián)大系統(tǒng)Lyapunov穩(wěn)定理論及Riccati方程設(shè)計變結(jié)構(gòu)控制方法;文獻[7]建立考慮控制通道間耦合因素的制導(dǎo)與控制一體化設(shè)計模型,設(shè)計非光滑擴張狀態(tài)觀測器進行觀測補償,然后結(jié)合塊動態(tài)逆方法設(shè)計保證系統(tǒng)的全局有限時間穩(wěn)定;文獻[8]針對飛行器耦合系統(tǒng),采用反步法保證系統(tǒng)狀態(tài)的一致最終有界。此類方法本質(zhì)上是將耦合視之為不確定性,通過魯棒控制方法(比如變結(jié)構(gòu)控制)或者擾動估計補償方法(比如擴張狀態(tài)觀測器)進行處理,缺乏對系統(tǒng)耦合特性的詳細分析。目前對于系統(tǒng)耦合對系統(tǒng)性能影響的研究較少,也沒有在高超聲速飛行器控制中的應(yīng)用。文獻[9]首次在理論上提出基于能量指標(biāo)的耦合利弊評價方法,并分析耦合對系統(tǒng)性能的影響,但其僅基于傳遞函數(shù)形式,缺乏對常用狀態(tài)空間形式的分析。

    針對以上問題,本文針對高超聲速飛行器姿態(tài)控制系統(tǒng),提出一種新型的基于耦合特性的魯棒控制方法(Coupling-based control, CBC),并結(jié)合干擾觀測器和滑模補償項保證其漸近收斂。本文與其他文獻研究的區(qū)別在于:1)提出基于耦合特性的飛行器姿態(tài)閉環(huán)控制器,控制器中反映了耦合對系統(tǒng)性能的影響;2)針對干擾觀測器估計不準導(dǎo)致跟蹤存在誤差的問題,結(jié)合魯棒補償項實現(xiàn)系統(tǒng)漸近收斂;3)基于絕對誤差積分指標(biāo)(Integral absolute error,IAE),從理論上證明本文提出的控制方法相對于未考慮耦合特性的控制具有更好的動態(tài)性能。

    1 飛行器數(shù)學(xué)模型

    1.1 高超聲速飛行器姿態(tài)系統(tǒng)模型

    飛行器姿態(tài)系統(tǒng)模型如下所示:

    (1)

    式中:α、β和γV分別是攻角、側(cè)滑角和傾側(cè)角,ωx、ωy和ωz分別是滾轉(zhuǎn)、偏航和俯仰角速度,θ為速度傾角,L和N代表升力和側(cè)向力,V是飛行器速度,g是重力加速度。Δα、Δβ和ΔγV表達式為:

    (2)

    式中:ψV為航跡偏航角,φ是飛行器緯度,R為飛行器相對地心的距離。由于實際中存在R>>V,所以一般在模型簡化處理中可以將Δα、Δβ和ΔγV考慮到不確定性中。另外,關(guān)于坐標(biāo)系定義、飛行器模型中各量的詳細定義和方程推導(dǎo)過程,本文不再贅述,具體請參考文獻[8]。

    下面給出飛行器的繞質(zhì)心動力學(xué)方程??紤]飛行器為面對稱,繞質(zhì)心動力學(xué)方程可以寫為

    (3)

    式中:Mx、My和Mz是飛行器受到的氣動力矩,Jx、Jy和Jz分別是飛行器相對機體坐標(biāo)系三軸的轉(zhuǎn)動慣量。

    (4)

    令ξ=[α,β,γV]T,ω= [ωz,ωy,ωx]T,u=[δz,δy,δx]T, 那么將式(4)代入式(1)和式(3)經(jīng)過簡化可得如下高超聲速飛行器姿態(tài)系統(tǒng)模型

    (5)

    式中:DF和DM代表飛行器受到的集總不確定性,包括未建模動態(tài)、氣動系數(shù)不確定性和外界擾動等。各系數(shù)矩陣的表達式見式(7),而gξ表達式為

    (6)

    (7)

    控制目標(biāo)是使得α、β和γV跟蹤參考指令αc、βc和γVc。

    對于飛行器模型,可以認為ω+gξ項是ξ系統(tǒng)的耦合項,其中ω代表姿態(tài)角速度對各通道角度的狀態(tài)耦合,而gξ代表通道耦合項,反映了各個通道之間的交聯(lián)影響。在角速度動力學(xué)方程中,Aωξξ項是ω系統(tǒng)的耦合項,Aωξξ代表姿態(tài)角對各自通道的角速度的狀態(tài)耦合。模型(5)中Aωω(ω)項包含的交叉耦合項是復(fù)雜的非線性耦合,為了簡便將之視為已知項,所以本文主要考慮狀態(tài)耦合和各通道耦合影響。

    接下來,需要把式(5)轉(zhuǎn)化為誤差系統(tǒng)。指令ξc=[αc,βc,γVc]T,設(shè)誤差狀態(tài)eξ=ξ-ξc,對其求導(dǎo)有

    (8)

    (9)

    進一步對eω求導(dǎo)可得

    (10)

    (11)

    如果控制器使得eξ→0,那么有[α,β,γV]T→ [αc,βc,γVc]T。下面,針對誤差模型(11),本文提出一種基于耦合特性的新型魯棒控制器。

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

    高超聲速飛行器的快速響應(yīng)特點和受到各種不確定性影響的特性,要求設(shè)計的控制系統(tǒng)不僅達到良好的動態(tài)性能,而且能夠具有抵抗不確定性的強魯棒性。為了滿足這種要求,本部分首先引入一種干擾觀測器設(shè)計方法;然后結(jié)合干擾觀測器提出考慮耦合特性的魯棒控制器,并證明系統(tǒng)的有界穩(wěn)定;最后結(jié)合滑模補償器實現(xiàn)了整個系統(tǒng)的漸近穩(wěn)定。

    2.1 干擾觀測器設(shè)計

    令x=[ξT,ωT]T,那么系統(tǒng)(5)變?yōu)槿缦滦问?/p>

    (12)

    其中:

    (13)

    受文獻[11]啟發(fā),本文提出新型干擾觀測器:

    (14)

    式中:z∈R6×1是觀測器狀態(tài),σ0和σ1是觀測器增益,其中矩陣Λ∈R6×6是待設(shè)計矩陣。

    引理 1. 針對系統(tǒng)(12)設(shè)計如式(14)的干擾觀測器,如果滿足假設(shè)1,而且σ0=2Λ,σ1=Λ2,矩陣Λ取為Λ=diag(λ,λ,λ,λ,λ,λ),其中λ> 0,那么干擾的估計誤差滿足

    (15)

    式中:r=μ/λ2。

    (16)

    (17)

    下面參考文獻[11]中引理1的求解方法可得結(jié)論(15)。

    2.2 基于耦合特性的魯棒控制

    定理 1. 考慮姿態(tài)誤差系統(tǒng)(11),設(shè)計如下的控制器

    (18)

    (19)

    式中:(°)表示Hadamard乘積,函數(shù)sgn(x) = [sgn(x1), sgn(x2), …, sgn(xn)]T(x=[x1,x2, …,xn]T)。Jξ1= [Jξ11,Jξ12,Jξ13]T,Jξ2= [Jξ21,Jξ22,Jξ23]T和Jω= [Jω1,Jω2,Jω3]T均是3×1矩陣。而Θ函數(shù)定義為

    Θ(J)=diag(sgn(J1),sgn(J2),…,sgn(Jn))

    (20)

    式中:矩陣J=[J1J2…Jn]T。那么姿態(tài)系統(tǒng)(11)的跟蹤誤差漸近穩(wěn)定到原點附近鄰域內(nèi)。

    證. 將控制器(18)代入系統(tǒng)(11)得到閉環(huán)系統(tǒng)

    (21)

    (22)

    容易證明不論eξieωi正負,必然有

    (23)

    1.形象思維和抽象思維的統(tǒng)一符合兒童思維和認知發(fā)展特點。兒童主要以形象思維來觀察世界。情境教學(xué)抓住兒童的這一思維特點,將“形真、情切、意遠、理寓其中”融為一個整體,一以貫之,將知識從具體到抽象,再從抽象到具體,讓兒童真正理解,這符合兒童思維和認知的發(fā)展特點。[1]

    (24)

    (25)

    注1. 控制器中的Θ函數(shù)即考慮耦合特性的控制器附加項。如果符號函數(shù)判斷為-1,控制律中將耦合進行抵消;符號函數(shù)判斷為1,則對耦合加強,最終保證耦合始終對系統(tǒng)性能有利。

    注2. 耦合特性量Jξ1、Jξ2和Jω是定義在誤差狀態(tài)基礎(chǔ)上的,所以耦合特性是在時刻變化的。這和經(jīng)典的耦合評價方法相比在本質(zhì)上是不同的。耦合評價大多針對開環(huán)系統(tǒng)耦合評估,但本文提出的耦合特性量可以設(shè)計閉環(huán)控制器。

    注3. 值得注意的是,本文提出的控制器中存在切換項,所以穩(wěn)定性要在Filippov意義[13]下理解。

    2.3 補償器設(shè)計

    定理 2. 考慮如式(11)的誤差姿態(tài)系統(tǒng),設(shè)計如下的控制器

    (26)

    證. 仍取Lyapunov函數(shù)V,將控制器(26)代入系統(tǒng)(11),對V求導(dǎo)并變換可得

    (27)

    3 基于絕對誤差積分的性能分析

    絕對誤差積分指標(biāo)是控制系統(tǒng)中一種常用的性能分析指標(biāo)[12]:

    (28)

    式中:e代表系統(tǒng)的跟蹤誤差。IAE指標(biāo)能夠綜合反映出系統(tǒng)跟蹤控制的動態(tài)特性。下面將基于此指標(biāo)分析本文方法的優(yōu)點。針對系統(tǒng)(11),不考慮耦合特性的常規(guī)控制器(Non-coupling-based control, NCBC)為如下形式

    u′=B-1(-Kωeω-eξ-Aωω(ω)-Aωξeξ-h+

    (29)

    定理 3. 不考慮耦合特性的控制器(29)的系統(tǒng)絕對誤差積分指標(biāo)為ΦIAE1,而采用基于耦合特性控制器(26)的指標(biāo)為ΦIAE2,那么有ΦIAE2≤ΦIAE1。

    證. 將控制器(29)代入系統(tǒng)(11)得到閉環(huán)系統(tǒng)

    (30)

    根據(jù)式(30)并結(jié)合定義(28)計算ΦIAE1,將式(26)代入式(11),并結(jié)合定義(28)可計算ΦIAE2,分別如下所示

    (31)

    ΦIAE2-ΦIAE1=

    (32)

    注意到

    (33)

    將式(33)代入式(32)可得ΦIAE2-ΦIAE1≤ 0。證畢。

    4 仿真校驗

    仿真校驗中飛行器的基本參數(shù)如表1所示[8]。

    表1 飛行器基本參數(shù)

    氣動系數(shù)取值參考文獻[14]。文中飛行器速度馬赫數(shù)為6,飛行高度為20km。跟蹤指令如表2所示。

    表2 姿態(tài)角跟蹤指令

    初始值分別為α0=2°,β0=0.1°,γV0=1°。執(zhí)行機構(gòu)限幅為[-20°, 20°],限速率為[-300°/s, 300°/s]。干擾觀測器參數(shù)λ=20,濾波器參數(shù)T=diag(30, 30, 30)??刂破鲄?shù)取值為

    (34)

    4.1 有效性校驗及對比仿真

    為了檢驗本文控制方法的有效性和魯棒性,t=6s時在仿真中加入外界擾動為

    (35)

    分別采用本文的控制器(26)和不考慮耦合特性的控制器(29)進行仿真對比,結(jié)果如圖1~4所示。

    圖1~3分別是姿態(tài)角、干擾估計和舵偏的對比變化曲線。

    從圖1可以看出,本文的控制方法保證了系統(tǒng)的漸近穩(wěn)定跟蹤,而且在相同的控制器參數(shù)下,相對于未考慮耦合特性的控制器NCBC,飛行器姿態(tài)系統(tǒng)實現(xiàn)了更快的上升時間和更小的超調(diào)量,獲得了更好的動態(tài)性能。

    由圖1可知,在加入外部擾動的情況下,本文提出的CBC控制方法依然能夠保證姿態(tài)角對跟蹤指令的精確跟蹤;而圖2反映了干擾估計值對真值的跟蹤情況,可以看出干擾估計值能夠較快跟蹤上真值,從而校驗了本文所設(shè)計干擾觀測器的有效性。仿真結(jié)果說明本文的高超聲速飛行器姿態(tài)控制方法具有較強的魯棒性。

    表3列出了本文設(shè)計方法CBC和未考慮耦合特性的NCBC之間IAE性能指標(biāo)的對比。從表3可以看出,三個通道中CBC方法相對于NCBC都具有更小的IAE指標(biāo),說明本文設(shè)計的控制方法使得姿態(tài)系統(tǒng)獲得了更好的動態(tài)性能。

    表3 IAE性能指標(biāo)對比

    從圖3可以看出,系統(tǒng)的舵偏出現(xiàn)了一定程度的抖動。這同樣是由于控制舵偏的速率限制和控制器中采用了切換項而導(dǎo)致的。實際仿真中可以考慮采用邊界層方法將切換函數(shù)sgn(·)替換為飽和函數(shù)sat(·),從而減小舵偏的抖振。

    4.2 氣動系數(shù)拉偏情況

    采用本文的控制律,考慮氣動系數(shù)拉偏,拉偏范圍為±20%。為簡便起見,仿真時間取為6s,仿真結(jié)果如圖4~6所示。

    圖4~6分別是俯仰、偏航和滾動通道的姿態(tài)角以及舵偏的變化曲線。從圖4可以看出,在正負拉偏情況下攻角均能夠保證對跟蹤指令的精確漸近跟蹤,且具有良好的動態(tài)性能。同樣地,圖5和圖6反映出偏航和滾動通道的良好動態(tài)特性。這說明本文設(shè)計的姿態(tài)控制器能夠適應(yīng)氣動系數(shù)的不確定性變化,具有較強魯棒性。

    5 結(jié) 論

    本文提出了一種新型的基于耦合特性的高超聲速飛行器姿態(tài)控制方法,然后基于干擾觀測器技術(shù)和滑模補償增強了控制器的魯棒性能。仿真校驗了本文方法具有較好魯棒性和動態(tài)性能。本文是作者對耦合性質(zhì)及控制器設(shè)計的初步研究,在以后的工作中將進行進一步的理論拓展和深入研究。

    [1] Duan H B, Li P. Progress in control approaches for hypersonic vehicle [J]. Science China - Technological Sciences, 2012, 55(10): 2965-2970.

    [2] Charles E H, Yuri B S. Sliding mode disturbance observer-based control for a reusable launch vehicle[J]. Journal of Guidance,Control, and Dynamics, 2006, 29(6): 1315-1328.

    [3] Georgie J, Valasek J. Evaluation of longitudinal desired dynamics for dynamic inversion controlled generic reentry vehicles [J]. Journal of Guidance, Control, and Dynamics, 2003, 26(5): 811-819.

    [4] Su X F, Jia Y M. Self-scheduled robust decoupling control withH∞performance of hypersonic vehicles [J]. Systems Control & Letters, 2014, 70: 38-48.

    [5] 周軍, 劉玉璽, 周鳳岐. 動態(tài)解耦和變結(jié)構(gòu)魯棒補償?shù)膶?dǎo)彈再入控制[J]. 系統(tǒng)仿真學(xué)報, 2009, 21(10): 3071-3073, 3083. [Zhou Jun, Liu Yu-xi, Zhou Feng-qi. Study of missile reentry dynamic decoupling control based on variable structure robustness compensation [J]. Journal of System Simulation, 2009, 21(10): 3071-3073, 3083.]

    [6] 周鳳岐, 王延, 周軍, 等. 高超聲速飛行器耦合系統(tǒng)變結(jié)構(gòu)控制設(shè)計[J]. 宇航學(xué)報, 2011, 32(1): 66-71. [Zhou Feng-qi, Wang Yan, Zhou Jun, et al. Design of variable structure controller for hypersonic coupled vehicle system [J]. Journal of Astronautics, 2011, 32(1): 66-71.]

    [7] 孫向宇, 晁濤, 王松艷, 等. 考慮通道耦合因素的制導(dǎo)控制一體化設(shè)計方法 [J]. 宇航學(xué)報, 2016, 37(8): 936-945. [Sun Xiang-yu, Chao Tao, Wang Song-yan, et al. Integrated guidance and control design method considering channel coupling [J]. Journal of Astronautics, 2016, 37(8): 936-945.]

    [8] 王建華, 劉魯華, 湯國建. 高超聲速飛行器俯沖段制導(dǎo)與姿控系統(tǒng)設(shè)計[J]. 宇航學(xué)報, 2016, 37(8): 964-973. [Wang Jian-hua, Liu Lu-hua, Tang Guo-jian. Guidance and attitude control system design for hypersonic vehicle in dive phase [J]. Journal of Astronautics, 2016, 37(8): 964-973.]

    [9] Gigi S, Tangirala A K. Quantification of interaction in multiloop control system using directed spectral decomposition [J]. Automatica, 2013, 49(5): 1174-1183.

    [10] Krstic M, Kanellakopoulos I, Kokotovic P V. Nonlinear and adaptive control design [M]. New York: John Wiley & Sons, 1995.

    [11] 張?zhí)煲? 周軍, 郭建國. 基于干擾觀測器的高速飛行器預(yù)測控制律設(shè)計[J]. 航空學(xué)報, 2013, 35(1): 215-222. [Zhang Tian-yi, Zhou Jun, Guo Jian-guo. Design of predictive controller for hypersonic vehicle based on disturbance observer[J]. Acta Aeronautica et Astronautica Sinica, 2013, 35(1): 215-222.]

    [12] Astrom K J, Hagglund T. PID controllers: theory, design and tuning [M]. Research Triangle Park, North Carolina: Instrument Society of America, 1995.

    [13] Filippov A F. Differential equations with discontinuous righthand sides [M]. Dordrecht, the Netherlands: Kluwer Academic Publishers, 1988.

    [14] Colgren R, Keshmiri S, Mirmirani M. Nonlinear ten- degree-of-freedom dynamics model of a generic hypersonic vehicle [J]. Journal of Aircraft,2009, 46(3) : 800-813.

    通信地址:陜西省西安市碑林區(qū)友誼西路西北工業(yè)大學(xué)234郵箱(710072)

    電話:(029)88493113

    E-mail:Eason.work@qq.com

    (編輯:牛苗苗)

    Novel Coupling Based Attitude Control System Design for Hypersonic Vehicles

    GUO Zong-yi, ZHOU Jun, GUO Jian-guo

    (Institute of Precision Guidance and Control, Northwestern Polytechnical University, Xi’an 710072,China)

    This paper presents a novel coupling characterization based robust control scheme for hypersonic vehicles. First of all, the attitude equations of a hypersonic vehicle are denoted, and the error dynamic model of the hypersonic vehicle is obtained based on the backstepping technique. Then, a robust coupling-based control scheme is proposed combining with the disturbance observer, and the asymptotic convergence is guaranteed through the sliding mode compensation. Finally, the integral absolute error (IAE) performance index is introduced, and the better performance of the proposed method is proven theoretically. The proposed method not only guarantees the asymptotic stability and the strong robustness, but also improves the dynamic performance via considering the coupling characterization into the controller design. The aerodynamic parameters perturbation and actual actuator constraints are considered in the numerical simulations, and the results demonstrate the validity of the control method in this paper.

    Hypersonic vehicle; Coupling property; Disturbance observer; Sliding mode compensation; Integral absolute error (IAE)

    2016-12-01;

    2017-01-17

    國家自然科學(xué)基金(61473226)

    V448

    A

    1000-1328(2017)03-0270-09

    10.3873/j.issn.1000-1328.2017.03.007

    郭宗易(1988-),男,博士生,主要從事飛行器制導(dǎo)控制與仿真研究。

    猜你喜歡
    超聲速觀測器飛行器
    高超聲速出版工程
    高超聲速飛行器
    超聲速旅行
    復(fù)雜飛行器的容錯控制
    電子制作(2018年2期)2018-04-18 07:13:25
    基于觀測器的列車網(wǎng)絡(luò)控制
    基于非線性未知輸入觀測器的航天器故障診斷
    神秘的飛行器
    基于干擾觀測器的PI控制單相逆變器
    高超聲速大博弈
    太空探索(2014年5期)2014-07-12 09:53:28
    采用干擾觀測器PI控制的單相SPWM逆變電源
    亚洲七黄色美女视频| 国产精品久久久久久精品电影| 我的老师免费观看完整版| 91av网一区二区| 免费av毛片视频| 一个人免费在线观看电影| 国产亚洲av嫩草精品影院| 热99re8久久精品国产| 亚洲美女搞黄在线观看 | 国产精品一区二区三区四区久久| 国产精品久久久久久av不卡| 欧美成人精品欧美一级黄| 国产av麻豆久久久久久久| 欧美日韩综合久久久久久| 日本与韩国留学比较| 天堂√8在线中文| 免费人成视频x8x8入口观看| 国产精品一区www在线观看| 亚洲天堂国产精品一区在线| 亚洲国产精品合色在线| 国产一区二区在线观看日韩| 日韩欧美免费精品| aaaaa片日本免费| 一卡2卡三卡四卡精品乱码亚洲| 桃色一区二区三区在线观看| 日韩欧美国产在线观看| 人妻制服诱惑在线中文字幕| 久久国产乱子免费精品| 男女之事视频高清在线观看| 麻豆乱淫一区二区| 91在线观看av| 国产爱豆传媒在线观看| 波多野结衣高清作品| 久久午夜福利片| 老司机福利观看| 尾随美女入室| 白带黄色成豆腐渣| 波野结衣二区三区在线| 精品不卡国产一区二区三区| 婷婷亚洲欧美| 一级黄色大片毛片| 国产精品一区二区性色av| 成人av一区二区三区在线看| 亚洲av中文av极速乱| 亚洲精品成人久久久久久| 久久久久久伊人网av| 天堂√8在线中文| 日韩欧美免费精品| 国模一区二区三区四区视频| 最近在线观看免费完整版| 欧美丝袜亚洲另类| 免费黄网站久久成人精品| 观看免费一级毛片| 亚洲第一区二区三区不卡| 久久婷婷人人爽人人干人人爱| av在线播放精品| 久久精品国产亚洲网站| 青春草视频在线免费观看| 可以在线观看毛片的网站| 午夜日韩欧美国产| 亚洲欧美日韩东京热| 人妻制服诱惑在线中文字幕| 麻豆久久精品国产亚洲av| 免费在线观看影片大全网站| 亚洲自拍偷在线| 欧美日韩精品成人综合77777| 97人妻精品一区二区三区麻豆| 国产精品久久久久久亚洲av鲁大| 久久欧美精品欧美久久欧美| 草草在线视频免费看| av国产免费在线观看| 五月伊人婷婷丁香| 国产高清视频在线观看网站| 十八禁网站免费在线| 日本精品一区二区三区蜜桃| 成人午夜高清在线视频| 亚洲久久久久久中文字幕| 欧美人与善性xxx| 欧美日本亚洲视频在线播放| 一级毛片久久久久久久久女| 亚洲精品456在线播放app| 91狼人影院| 91午夜精品亚洲一区二区三区| 天堂√8在线中文| 亚洲四区av| 色尼玛亚洲综合影院| 少妇高潮的动态图| 国产精品av视频在线免费观看| 综合色丁香网| 成人亚洲精品av一区二区| 亚洲,欧美,日韩| 国产av在哪里看| 久久精品综合一区二区三区| 久99久视频精品免费| 日韩大尺度精品在线看网址| 中文字幕av成人在线电影| 深夜精品福利| 国模一区二区三区四区视频| 国产又黄又爽又无遮挡在线| 日本黄色视频三级网站网址| 国产精品,欧美在线| 性插视频无遮挡在线免费观看| 午夜激情欧美在线| 国产中年淑女户外野战色| 亚洲国产高清在线一区二区三| 欧美日本视频| 麻豆精品久久久久久蜜桃| 欧美+日韩+精品| 波野结衣二区三区在线| 成人亚洲欧美一区二区av| 亚洲第一区二区三区不卡| 日产精品乱码卡一卡2卡三| 国产不卡一卡二| 黑人高潮一二区| 一进一出好大好爽视频| 男女边吃奶边做爰视频| 美女cb高潮喷水在线观看| 亚洲一区高清亚洲精品| 99热这里只有精品一区| 中国美白少妇内射xxxbb| 中文亚洲av片在线观看爽| 日韩精品青青久久久久久| 久久中文看片网| 黄色视频,在线免费观看| 国产黄色小视频在线观看| 亚洲va在线va天堂va国产| 国产视频内射| 精品一区二区三区av网在线观看| 国产精品乱码一区二三区的特点| 国产极品精品免费视频能看的| 日本色播在线视频| 国产精品三级大全| 亚洲国产精品成人综合色| 热99re8久久精品国产| 男女那种视频在线观看| 国产探花极品一区二区| 免费在线观看影片大全网站| eeuss影院久久| 99久久中文字幕三级久久日本| 看黄色毛片网站| av天堂中文字幕网| 国产高清视频在线播放一区| 国产精品国产三级国产av玫瑰| 超碰av人人做人人爽久久| 久久久国产成人免费| 校园春色视频在线观看| 十八禁国产超污无遮挡网站| 内射极品少妇av片p| 国产极品精品免费视频能看的| 少妇丰满av| 天天一区二区日本电影三级| 嫩草影视91久久| 黑人高潮一二区| 性插视频无遮挡在线免费观看| 啦啦啦啦在线视频资源| 日韩欧美在线乱码| 免费看a级黄色片| 欧美一区二区国产精品久久精品| 内地一区二区视频在线| 亚洲精品亚洲一区二区| 美女内射精品一级片tv| 少妇的逼水好多| 亚洲成人精品中文字幕电影| 少妇猛男粗大的猛烈进出视频 | 久久人妻av系列| 国产高清不卡午夜福利| 国产毛片a区久久久久| 国产av一区在线观看免费| 免费一级毛片在线播放高清视频| 免费看光身美女| av.在线天堂| 日韩精品青青久久久久久| 3wmmmm亚洲av在线观看| 毛片一级片免费看久久久久| 亚洲av不卡在线观看| 亚洲欧美成人精品一区二区| 精品一区二区三区人妻视频| 国产伦精品一区二区三区视频9| 男女之事视频高清在线观看| 日本撒尿小便嘘嘘汇集6| 亚洲在线自拍视频| 成年女人看的毛片在线观看| 一进一出抽搐gif免费好疼| 亚洲精品成人久久久久久| 欧美一区二区亚洲| 韩国av在线不卡| 国产在线男女| 99久久久亚洲精品蜜臀av| 国产精品国产三级国产av玫瑰| 亚洲高清免费不卡视频| 一本精品99久久精品77| 欧美激情久久久久久爽电影| 一级毛片电影观看 | 欧美区成人在线视频| 精品熟女少妇av免费看| 在线天堂最新版资源| 美女cb高潮喷水在线观看| 久久精品国产99精品国产亚洲性色| 级片在线观看| 免费观看在线日韩| 久久久久久伊人网av| 91久久精品国产一区二区三区| 亚洲四区av| 真实男女啪啪啪动态图| 国产白丝娇喘喷水9色精品| 美女免费视频网站| 综合色丁香网| 色视频www国产| 97超碰精品成人国产| 97碰自拍视频| 亚洲无线在线观看| 久久九九热精品免费| 亚洲人成网站高清观看| 插逼视频在线观看| 美女黄网站色视频| 亚洲婷婷狠狠爱综合网| 日本爱情动作片www.在线观看 | 国产精品久久视频播放| 最近手机中文字幕大全| 晚上一个人看的免费电影| 国产伦在线观看视频一区| 亚洲激情五月婷婷啪啪| 国产成人福利小说| 日韩 亚洲 欧美在线| 99热精品在线国产| 22中文网久久字幕| 97在线视频观看| 精品人妻熟女av久视频| av在线老鸭窝| 神马国产精品三级电影在线观看| 高清毛片免费观看视频网站| 久久鲁丝午夜福利片| 亚洲av成人av| 别揉我奶头 嗯啊视频| 此物有八面人人有两片| 久久九九热精品免费| 天天一区二区日本电影三级| 女生性感内裤真人,穿戴方法视频| 午夜精品国产一区二区电影 | 亚洲经典国产精华液单| 国产伦一二天堂av在线观看| ponron亚洲| 午夜a级毛片| 天美传媒精品一区二区| 91精品国产九色| 淫妇啪啪啪对白视频| 亚洲中文字幕一区二区三区有码在线看| 99热只有精品国产| 亚洲成人久久爱视频| 国产伦精品一区二区三区四那| 亚洲激情五月婷婷啪啪| 日韩欧美精品v在线| 日韩欧美免费精品| 成人av在线播放网站| 国产日本99.免费观看| 99热精品在线国产| 日韩一本色道免费dvd| 黄色视频,在线免费观看| 性欧美人与动物交配| 日韩,欧美,国产一区二区三区 | 国产精品永久免费网站| 精品不卡国产一区二区三区| 我的女老师完整版在线观看| 亚洲最大成人手机在线| 久久久精品大字幕| 12—13女人毛片做爰片一| 长腿黑丝高跟| 欧美成人精品欧美一级黄| 午夜福利高清视频| videossex国产| 熟妇人妻久久中文字幕3abv| 九九热线精品视视频播放| 国产又黄又爽又无遮挡在线| 少妇丰满av| 日本撒尿小便嘘嘘汇集6| 在线天堂最新版资源| 91在线精品国自产拍蜜月| 赤兔流量卡办理| 亚洲av免费在线观看| 午夜精品国产一区二区电影 | 在线观看一区二区三区| 蜜桃亚洲精品一区二区三区| 人人妻,人人澡人人爽秒播| 久久午夜亚洲精品久久| 你懂的网址亚洲精品在线观看 | 噜噜噜噜噜久久久久久91| 国产 一区 欧美 日韩| 看非洲黑人一级黄片| 最近的中文字幕免费完整| 中文亚洲av片在线观看爽| 亚洲av熟女| 中国美白少妇内射xxxbb| 国产国拍精品亚洲av在线观看| 久久久久久久久久久丰满| 午夜福利18| 最近2019中文字幕mv第一页| 我的老师免费观看完整版| 欧美成人一区二区免费高清观看| h日本视频在线播放| 久久综合国产亚洲精品| 尾随美女入室| 成年女人永久免费观看视频| 成人一区二区视频在线观看| 国产高清视频在线观看网站| 亚洲在线自拍视频| 久久久久久久久大av| 女的被弄到高潮叫床怎么办| 国产午夜精品久久久久久一区二区三区 | 欧美日韩乱码在线| 欧美+日韩+精品| 日韩强制内射视频| 18+在线观看网站| 国产成人福利小说| 国产av麻豆久久久久久久| 五月玫瑰六月丁香| 黄色欧美视频在线观看| 亚洲熟妇中文字幕五十中出| 国产在线精品亚洲第一网站| 深夜精品福利| 国内精品久久久久精免费| 亚洲av.av天堂| 精品久久久久久成人av| 久久精品国产自在天天线| 最新在线观看一区二区三区| 亚洲精品在线观看二区| 少妇熟女欧美另类| 人人妻人人澡欧美一区二区| 在线看三级毛片| 91久久精品电影网| 亚洲成人久久爱视频| 国产免费一级a男人的天堂| 亚洲精品久久国产高清桃花| 欧美区成人在线视频| 男女视频在线观看网站免费| 成年版毛片免费区| 美女内射精品一级片tv| 一级黄色大片毛片| 插阴视频在线观看视频| 插逼视频在线观看| 久久国内精品自在自线图片| .国产精品久久| 一区二区三区四区激情视频 | 亚洲真实伦在线观看| 最近的中文字幕免费完整| 老熟妇仑乱视频hdxx| 99久久精品国产国产毛片| 国产精品亚洲一级av第二区| 性色avwww在线观看| 在线a可以看的网站| 欧美激情在线99| 天天一区二区日本电影三级| 日韩av在线大香蕉| 欧美最新免费一区二区三区| 国产精品久久视频播放| 在线天堂最新版资源| 日韩高清综合在线| av福利片在线观看| 春色校园在线视频观看| 久久精品91蜜桃| 亚洲精品国产av成人精品 | 日韩 亚洲 欧美在线| 国产精品综合久久久久久久免费| 简卡轻食公司| 亚洲精品日韩在线中文字幕 | 日韩在线高清观看一区二区三区| 中国国产av一级| 午夜视频国产福利| 91精品国产九色| 91av网一区二区| 蜜桃亚洲精品一区二区三区| 男女下面进入的视频免费午夜| 欧美极品一区二区三区四区| 99久久精品热视频| 久久久a久久爽久久v久久| 丰满的人妻完整版| 亚洲人成网站高清观看| 床上黄色一级片| 联通29元200g的流量卡| 麻豆国产97在线/欧美| 美女被艹到高潮喷水动态| 久久精品影院6| 国产精品电影一区二区三区| or卡值多少钱| 色5月婷婷丁香| 超碰av人人做人人爽久久| 欧美+日韩+精品| 校园人妻丝袜中文字幕| 大香蕉久久网| 99国产精品一区二区蜜桃av| 亚洲经典国产精华液单| 春色校园在线视频观看| 国产高潮美女av| 亚洲成人久久爱视频| 亚洲精品亚洲一区二区| a级毛色黄片| 午夜福利高清视频| 亚洲va在线va天堂va国产| 国产精品一区www在线观看| 高清毛片免费观看视频网站| 伊人久久精品亚洲午夜| or卡值多少钱| 国产激情偷乱视频一区二区| 国产精品一及| 国产v大片淫在线免费观看| 精品人妻偷拍中文字幕| 午夜精品在线福利| av天堂中文字幕网| 91在线精品国自产拍蜜月| 欧美绝顶高潮抽搐喷水| 我要搜黄色片| 欧美日本视频| 成人漫画全彩无遮挡| 中文字幕av在线有码专区| 变态另类丝袜制服| 女的被弄到高潮叫床怎么办| 亚洲国产精品sss在线观看| 国产淫片久久久久久久久| 在线天堂最新版资源| 夜夜爽天天搞| 少妇的逼水好多| 波多野结衣高清作品| 欧美极品一区二区三区四区| 成人三级黄色视频| 麻豆乱淫一区二区| 久久久久久九九精品二区国产| 99riav亚洲国产免费| 久久久久国产精品人妻aⅴ院| 国产精品日韩av在线免费观看| 亚洲在线自拍视频| 婷婷六月久久综合丁香| 精品午夜福利视频在线观看一区| 国产一区二区三区av在线 | 国产精品爽爽va在线观看网站| 免费人成在线观看视频色| 亚洲18禁久久av| 亚洲在线观看片| 久久婷婷人人爽人人干人人爱| 国产精品久久久久久久久免| 亚洲国产精品sss在线观看| 99热6这里只有精品| 麻豆久久精品国产亚洲av| 亚洲高清免费不卡视频| 国产白丝娇喘喷水9色精品| 一个人免费在线观看电影| 成人精品一区二区免费| 国产高清视频在线观看网站| 亚洲第一电影网av| eeuss影院久久| 久久久国产成人免费| 亚洲精品一卡2卡三卡4卡5卡| 日本a在线网址| 精品久久久久久久人妻蜜臀av| 久久久久久伊人网av| 亚洲成人久久性| 一区二区三区四区激情视频 | 一个人看视频在线观看www免费| 国产精品一及| 精品久久久久久成人av| 精品国产三级普通话版| 国产在线精品亚洲第一网站| 成人特级av手机在线观看| 精品少妇黑人巨大在线播放 | 午夜亚洲福利在线播放| 成人亚洲欧美一区二区av| 好男人在线观看高清免费视频| 日韩一本色道免费dvd| 一级av片app| www.色视频.com| 国内精品久久久久精免费| 久久久精品欧美日韩精品| 内射极品少妇av片p| 卡戴珊不雅视频在线播放| 亚洲第一区二区三区不卡| 午夜激情福利司机影院| 99久久久亚洲精品蜜臀av| 69人妻影院| 真人做人爱边吃奶动态| 99在线人妻在线中文字幕| 一本精品99久久精品77| 国内揄拍国产精品人妻在线| 中出人妻视频一区二区| 色在线成人网| 中出人妻视频一区二区| 亚洲成人中文字幕在线播放| 色av中文字幕| 男人舔女人下体高潮全视频| 免费看av在线观看网站| 日日啪夜夜撸| 在线免费观看不下载黄p国产| 人妻夜夜爽99麻豆av| 国产精品嫩草影院av在线观看| 亚洲欧美精品自产自拍| 亚洲av免费高清在线观看| 成年女人看的毛片在线观看| 最后的刺客免费高清国语| 三级经典国产精品| 搡老妇女老女人老熟妇| 亚洲国产日韩欧美精品在线观看| 日韩强制内射视频| 午夜视频国产福利| 高清日韩中文字幕在线| 国产v大片淫在线免费观看| 午夜福利18| a级毛色黄片| 日韩精品青青久久久久久| 日韩 亚洲 欧美在线| 大型黄色视频在线免费观看| 亚洲国产精品成人久久小说 | 不卡一级毛片| 99久国产av精品国产电影| 亚洲人成网站在线播放欧美日韩| 一边摸一边抽搐一进一小说| 小说图片视频综合网站| 日韩高清综合在线| 九九在线视频观看精品| 97热精品久久久久久| 欧美日韩综合久久久久久| 亚洲精品国产av成人精品 | 久久九九热精品免费| 悠悠久久av| 欧美成人精品欧美一级黄| 国产高清视频在线观看网站| 日韩欧美在线乱码| 俄罗斯特黄特色一大片| 国产乱人偷精品视频| 国产成人一区二区在线| 久久久精品欧美日韩精品| 麻豆一二三区av精品| 成年女人永久免费观看视频| 亚洲丝袜综合中文字幕| av在线天堂中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 一进一出抽搐动态| avwww免费| 免费一级毛片在线播放高清视频| 日韩三级伦理在线观看| 中文字幕久久专区| 两个人视频免费观看高清| 两个人的视频大全免费| 中文字幕免费在线视频6| 亚洲熟妇熟女久久| 黄色配什么色好看| 国产精品美女特级片免费视频播放器| 日本 av在线| 欧美一级a爱片免费观看看| 最近中文字幕高清免费大全6| 亚洲一区二区三区色噜噜| 国产单亲对白刺激| 国产高清视频在线观看网站| 黄色配什么色好看| 亚洲成人久久爱视频| 亚洲美女视频黄频| 亚洲中文日韩欧美视频| 国产成人影院久久av| 亚洲真实伦在线观看| 国产高清三级在线| 一个人免费在线观看电影| 麻豆久久精品国产亚洲av| 老司机福利观看| 色哟哟·www| 一个人看的www免费观看视频| 美女 人体艺术 gogo| 国产精品一及| 精品少妇黑人巨大在线播放 | 狂野欧美激情性xxxx在线观看| 亚洲精品一卡2卡三卡4卡5卡| 免费在线观看成人毛片| 久久精品国产清高在天天线| 成人永久免费在线观看视频| 毛片女人毛片| 色噜噜av男人的天堂激情| 中文字幕av在线有码专区| 插阴视频在线观看视频| 日本色播在线视频| 国产精品野战在线观看| 精品久久久久久久末码| 最好的美女福利视频网| 啦啦啦观看免费观看视频高清| 免费搜索国产男女视频| 日产精品乱码卡一卡2卡三| 男女那种视频在线观看| 女人十人毛片免费观看3o分钟| 老熟妇仑乱视频hdxx| 亚洲第一电影网av| 少妇裸体淫交视频免费看高清| 我的老师免费观看完整版| 亚洲精品乱码久久久v下载方式| 狂野欧美激情性xxxx在线观看| 久久精品91蜜桃| 午夜福利成人在线免费观看| 国产精品爽爽va在线观看网站| 久久精品91蜜桃| 午夜精品国产一区二区电影 | 一级毛片我不卡| 国产一区二区三区av在线 | 波多野结衣高清无吗| 观看免费一级毛片| 欧美激情在线99| 99国产极品粉嫩在线观看| 久久久久久久亚洲中文字幕| 成人午夜高清在线视频| 日日啪夜夜撸| 亚洲在线观看片| 嫩草影院入口| 午夜影院日韩av| 成人欧美大片| 日韩欧美精品免费久久| av国产免费在线观看| 看黄色毛片网站| 一本一本综合久久| 九九久久精品国产亚洲av麻豆| 欧美激情国产日韩精品一区| 可以在线观看的亚洲视频| 午夜a级毛片| 国产精品亚洲美女久久久|