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

    考慮狀態(tài)約束的彈性高超聲速飛行器自適應(yīng)飽和容錯控制

    2021-09-07 06:48:44譚立國魏毅寅段廣仁
    宇航學(xué)報 2021年7期
    關(guān)鍵詞:李雅普諾夫超聲速

    陳 峣,譚立國,魏毅寅,段廣仁

    (1. 哈爾濱工業(yè)大學(xué)控制理論與制導(dǎo)技術(shù)研究中心,哈爾濱 150001;2. 哈爾濱工業(yè)大學(xué)空間基礎(chǔ)科學(xué)研究中心,哈爾濱 150001)

    0 引 言

    高超聲速飛行器具有時變快、非線性強、耦合強及不確定性大等特點[1-2]。通常采用細(xì)長體外形和輕結(jié)構(gòu)設(shè)計,在發(fā)生彈性形變時可能會對飛行器的控制量產(chǎn)生影響,進而加大了執(zhí)行器發(fā)生故障的概率[3]。另外,高超聲速飛行器主要以超燃沖壓發(fā)動機作為動力推進系統(tǒng),飛行狀態(tài)變量的變化對其工作性能影響很大,只有當(dāng)飛行狀態(tài)滿足相關(guān)約束條件時,超燃沖壓發(fā)動機才可以正常工作。所以,研究執(zhí)行器故障和狀態(tài)約束下的彈性高超聲速飛行器的魯棒控制方法,對提升高超聲速飛行器的可靠性和安全性具有重要意義。

    近年來,隨著現(xiàn)代控制理論的快速發(fā)展,非線性控制方法在高超聲速飛行器跟蹤控制領(lǐng)域中得到了廣泛應(yīng)用,并取得豐富的研究成果。文獻[4-5]利用滑??刂坪妥赃m應(yīng)控制,針對高超聲速飛行器剛體模型,設(shè)計了自適應(yīng)滑??刂破?。文獻[6]采用跟蹤微分器,設(shè)計了高超聲速飛行器動態(tài)面控制器。文獻[7]針對彈性高超聲速飛行器,將滑??刂婆c反步法相結(jié)合,設(shè)計了反步滑??刂破?。為增強系統(tǒng)的魯棒性,文獻[8-9]將外界干擾、模型參數(shù)不確定性視為未知有界的系統(tǒng)干擾,利用非齊次觀測器對干擾進行估計,設(shè)計了高超聲速飛行器自適應(yīng)快速終端跟蹤控制器。進一步考慮彈性模態(tài)對高超聲速飛行器的影響,文獻[10-11]利用超螺旋滑模控制理論,針對彈性高超聲速飛行器設(shè)計了自適應(yīng)超螺旋滑模跟蹤控制器。文獻[12-15]利用智能控制算法,對彈性高超聲速飛行器跟蹤控制問題進行了研究分析。另外,高超聲速飛行器的動力推進系統(tǒng)為超燃沖壓發(fā)動機,只有當(dāng)飛行狀態(tài)滿足一定約束時,該發(fā)動機才能正常工作[12-14]。為解決上述問題,文獻[16-17]通過引入屏障李雅普諾夫函數(shù),針對非線性系統(tǒng)設(shè)計了滿足狀態(tài)約束條件的自適應(yīng)控制器。文獻[18-19]利用屏障李雅普諾夫函數(shù)、反步法和自適應(yīng)技術(shù),設(shè)計了能同時滿足跟蹤性能和飛行過程中狀態(tài)約束的控制器。

    在實際控制系統(tǒng)中,由于物理機構(gòu)的限制使得執(zhí)行機構(gòu)提供的控制力(力矩)是有限大小的,忽略執(zhí)行機構(gòu)的飽和會使得所設(shè)計的控制器魯棒性降低。因此,在設(shè)計控制器時必須考慮輸入飽和問題。文獻[20-21]通過引入輔助系統(tǒng),利用滑??刂评碚摵妥赃m應(yīng)控制算法,設(shè)計了剛體高超聲速飛行器飽和跟蹤控制器。文獻[22-24]針對彈性高超聲速飛行器設(shè)計了抗飽和的自適應(yīng)反步跟蹤控制器。文獻[25]在切換控制理論的基礎(chǔ)上,利用線性矩陣不等式方案,設(shè)計了抗飽和切換控制器。此外,由于一體化結(jié)構(gòu)設(shè)計以及高溫,高速等復(fù)雜飛行條件的影響,使得飛行器執(zhí)行機構(gòu)容易發(fā)生飽和,這將進一步提高執(zhí)行器發(fā)生故障的頻率,從而導(dǎo)致系統(tǒng)性能下降。文獻[26]在二階滑模趨近律和有限時間觀測器的基礎(chǔ)上,設(shè)計了彈性高超聲速飛行器故障容錯控制器。文獻[27] 針對高超聲速飛行器,根據(jù)被動容錯思想,設(shè)計了具有容錯功能的控制器。文獻[28-29] 設(shè)計了魯棒自適應(yīng)容錯控制器,該控制器可以保證系統(tǒng)模型參數(shù)不確定性和執(zhí)行機構(gòu)故障情況下高超聲速飛行器的穩(wěn)定性。文獻[30] 設(shè)計了抗飽和容錯控制器,該控制器以滑模控制和自適應(yīng)控制器為基礎(chǔ),可以保證控制輸出滿足執(zhí)行器的物理約束條件。

    為進一步解決帶有多種約束條件下的彈性高超聲速飛行器跟蹤控制問題,本文采用被動容錯控制的思想,結(jié)合自適應(yīng)控制、反步控制、設(shè)計了跟蹤控制器。與上述相關(guān)文獻相比,本文的主要創(chuàng)新點如下:

    1) 在控制器設(shè)計過程中,通過引入新型正切型屏障李雅普諾夫函數(shù)來確保所設(shè)計的控制器能夠滿足高超聲速飛行器的狀態(tài)約束和跟蹤性能。

    2) 與文獻[8]相比,本文通過引入雙曲正切函數(shù)處理輸入飽和問題,保證高超聲速飛行器在實際控制系統(tǒng)執(zhí)行過程中滿足執(zhí)行器機構(gòu)物理約束條件。

    3)與文獻[20]相比,本文同時考慮了輸入飽和、執(zhí)行器故障、狀態(tài)約束等條件,使得所設(shè)計的控制策略,具有更好的工程意義。

    本文的主要內(nèi)容如下:首先,給出了彈性高超聲速飛行器控制模型;其次,利用反步法、雙曲正切函數(shù)和自適應(yīng)控制技術(shù),分別針對速度子系統(tǒng)和高度子系統(tǒng)設(shè)計了自適應(yīng)抗飽和故障容錯跟蹤控制器,且借助Lyapunov理論證明了所設(shè)計的跟蹤控制器的速度和高度有限時間收斂性;再次,結(jié)合仿真實驗分析了所設(shè)計控制器的性能;最后,給出本文的結(jié)論。

    1 模型描述及相關(guān)引理

    1.1 彈性高超聲速飛行器模型

    考慮到高超聲速飛行器剛體動力學(xué)和彈性動力學(xué)氣動耦合,系統(tǒng)模型[2]:

    (1)

    推力T、阻力D、升力L、俯仰角力矩M和廣義力Ni表達式為:

    (2)

    式中:氣動力和氣動力矩系數(shù)CT,φ(α,Δτ1,M∞)、CT(α,Δτ1,M∞,Ad)、CD(α,δe,δc,Δτ1,Δτ2)、CL(α,δe,δc,Δτ1,Δτ2)、CM(α,δe,δc,Δτ1,Δτ2)、CNi(α,δe,δc,Δτ1,Δτ2)為關(guān)于M∞、Ad、α、Δτ1和Δτ2的非線性函數(shù),具體的表達式如下所示:

    (3)

    1.2 控制面模型

    考慮氣動參數(shù)的不確定性,可得到如下的氣動力、氣動力矩以及廣義力的不確定部分的表達式:

    (4)

    其中:ΔCT,φ、 ΔCT、 ΔCD、 ΔCL、 ΔCM及ΔCNi的表達式如下:

    (5)

    在高超聲速飛行器剛體模型中將彈性模態(tài)部分視為系統(tǒng)擾動進行處理,考慮到空氣參數(shù)的不確定性,則進一步控制模型為:

    (6)

    其中:dV、dγ、dα和dq為氣動參數(shù)和彈性模態(tài)引起的系統(tǒng)不確定項,則dV、dγ、dα和dq的具體表述式如下所示:

    其中,S為參數(shù)特征面積。

    同時考慮輸入飽和及執(zhí)行器故障模型如下:

    φ=ρVsat(φ),δe=ρhsat(δe)

    (7)

    其中:ρV和ρh為失效故障因子,且滿足0<ρV<1, 0<ρh<1, sat(φ)為飽和函數(shù)。

    注1:在執(zhí)行器故障模型(8)中,φ和δe分別高超聲速飛行器控制系統(tǒng)中的發(fā)動機節(jié)流閥和舵偏角的參考輸入;ρV和ρh分別表示為速度及高度執(zhí)行機構(gòu)中發(fā)生的故障,且滿足于0≤ρi<1(i=V,h)。

    根據(jù)式(7),式(6)可重寫為:

    (8)

    控制目標(biāo):針對彈性高超聲速飛行控制系統(tǒng)式(8),在控制限幅、執(zhí)行機構(gòu)失效的情況下,在所設(shè)計的抗飽和故障容錯控制策略作用下,能夠控制飛行器對速度參考信號Vd和高度參考信號hd進行快速穩(wěn)定高精度跟蹤,同時保證攻角、俯仰角、俯仰角速率滿足給定狀態(tài)約束條件。

    2 控制器設(shè)計

    針對速度和高度兩個子系統(tǒng),根據(jù)被動容錯思想、結(jié)合反步法和自適應(yīng)控制方法,分別設(shè)計了自適應(yīng)抗飽和的故障容錯控制器,通過引入新型屏障李雅普諾夫函數(shù),能夠保證系統(tǒng)狀態(tài)在飛行過程中滿足狀態(tài)約束條件。

    為方便控制器的設(shè)計,給出以下引理。

    引理1[27].對于任意實數(shù)x和非零實數(shù)y,下面不等式成立:

    0≤|x|(1-tanh(|x/y|))≤α|y|

    (9)

    其中:α>0,其最小值α*滿足α*=x*(1-tanhx*),則x*滿足方程e-2x*+1-2x*=0。

    2.1 速度子系統(tǒng)控制器設(shè)計

    定義速度跟蹤誤差zV:

    zV=V-Vd

    (10)

    其中:Vd為速度參考信號。

    對式(10)求導(dǎo):

    (11)

    為處理執(zhí)行器輸入飽和問題,引入光滑雙曲正切函數(shù)如下:

    sat(φ)=h(φ)=φ+Δφ

    (12)

    h(φ)=φmaxtanh(φ/φmax)

    (13)

    其中:φmax為正常數(shù)。

    根據(jù)式(12),則(11)可重寫為:

    (14)

    假設(shè)1:在系統(tǒng)(14)中,干擾gVρVΔφ+dV(t)存在未知上界,則滿足下列不等式:

    |gVρVΔφ+dV(t)|≤εV

    (15)

    其中:εV為未知正常數(shù)。

    為了保證速度期望的跟蹤性能|zV|≤AzV,其中AzV>0,定義屏障李雅普諾夫函數(shù)如下:

    (16)

    對式(16)求導(dǎo)可得:

    (17)

    根據(jù)式(17),利用反步法和自適應(yīng)算法,設(shè)計了帶有狀態(tài)約束的自適應(yīng)飽和故障容錯控制器:

    (18)

    (19)

    (20)

    其中:kV1、ξV、λV1和λV 2為正常數(shù)。

    定理1.考慮系統(tǒng)(8),且滿足假設(shè)條件1,在控制器(18)和自適應(yīng)律(19)~(20)作用下,則跟蹤誤差變量zV收斂到如下區(qū)域:

    (21)

    其中:cV和cV 0是與控制器參數(shù)相關(guān)的常數(shù)。

    證明:選取李雅普諾夫函數(shù)如下:

    (22)

    利用式(18)對V1求導(dǎo)整理可得:

    (23)

    由引理1可知,下列不等式成立:

    (24)

    將式(24)代入式(23)整理可得:

    (25)

    cV=min{kV1,λV1,λV 2}

    (26)

    (27)

    根據(jù)式(26)~(27),式(25)可以整理為:

    (28)

    對式(28)求解可得:

    (29)

    由式(29)可得V1(t)≤V1(0),即V1(t)是有界的,聯(lián)合屏障李雅普諾夫函數(shù)式(16),可得不等式如下:

    (30)

    進一步可得:

    (31)

    定理1證畢。

    2.2 高度子系統(tǒng)控制器設(shè)計

    高度誤差變量zh定義為:

    zh=h-hd

    (32)

    其中:hd為速度參考信號。

    對式(32)求導(dǎo):

    (33)

    進一步可得期望的航跡角信號:

    (34)

    其中:kh1>0是常數(shù)。

    為了系統(tǒng)的反饋形式,下式將式(8)進行轉(zhuǎn)化:

    (35)

    假設(shè)2.在高度系統(tǒng)模型(35)中,總的干擾項dγ和dα有界,且滿足式:

    |dγ|≤εγ, |dα|≤εα

    (36)

    其中:εγ>0和εα>0且未知。

    Step1:航跡角跟蹤誤差變量zγ定義:

    zγ=γ-γc

    (37)

    對式(37)求導(dǎo)可得:

    (38)

    為了保證航跡角期望的跟蹤性能|zγ|≤Azγ,其中:Azγ>0,定義屏障李雅普諾夫函數(shù)如下:

    (39)

    對式(39)求導(dǎo)可得:

    (40)

    根據(jù)式(40),設(shè)計虛擬控制器αc如下:

    (41)

    (42)

    其中:kγ1、ξγ和λγ為正常數(shù)。

    為了克服對虛擬控制器αc多次微分所導(dǎo)致計算復(fù)雜問題,引入一階命令濾波器如下:

    (43)

    其中:αc和αd分別為一階命令濾波器的輸入及輸出信號,且τ1為正常數(shù)。

    定義濾波器微分跟蹤誤差yα為:

    yα=αd-αc

    (44)

    將虛擬控制器式(41)代入式(40)整理可得:

    (45)

    Step2:定義攻角誤差zα=α-αc,并求導(dǎo)得:

    (46)

    為保證攻角期望的跟蹤性能|zα|≤Azα,其中Azα>0,定義屏障李雅普諾夫函數(shù)如下:

    (47)

    對式(47)求導(dǎo)可得:

    (48)

    根據(jù)式(46),設(shè)計虛擬控制器qc如下:

    (49)

    (50)

    其中:kα1,ξα和λα為正常數(shù)。

    為克服對虛擬控制器qc多次微分所導(dǎo)致計算復(fù)雜,引入一階命令濾波器如下:

    (51)

    其中:qc和qd分別為命令濾波器的輸入和輸出信號,且τ2為正常數(shù)。

    定義濾波器微分跟蹤誤差yq為:

    yq=qd-qc

    (52)

    將式(49)代入(48)整理可得:

    (53)

    Step3:定義俯仰角速率跟蹤誤差zq=q-qd,并求導(dǎo)可得:

    (54)

    為了處理執(zhí)行輸入飽和問題,引入光滑雙曲正切函數(shù)如下:

    sat(δe)=h(δe)=δe+Δδe

    (55)

    h(δe)=δemaxtanh(δe/δemax)

    (56)

    其中:δemax為正常數(shù)。則式(54)可重寫為:

    (57)

    假設(shè)3:在系統(tǒng)(57)中,干擾gqρhΔδe+dq存在未知上界,則滿足下列不等式:

    |gqρhΔδe+dq|≤εh

    (58)

    其中,εh為未知正常數(shù)。

    為保證俯仰角速率期望的跟蹤性能|zq|≤Azq,其中Azq>0,定義屏障李雅普諾夫函數(shù)如下:

    (59)

    對式(59)求導(dǎo)可得:

    (60)

    根據(jù)式(57),結(jié)合自適應(yīng)控制算法,設(shè)計了帶有狀態(tài)約束的飽和故障容錯控制器如下:

    (61)

    (62)

    (63)

    其中:kq1、ξq、λq1和λq2為正常數(shù)。

    定理2.考慮系統(tǒng)(35)且滿足假設(shè)條件2,在控制器式(61)和自適應(yīng)律式(62)~(63)作用下為制導(dǎo)律,誤差變量zh、zγ、zα和zq漸近收斂到如下區(qū)域。

    (64)

    其中:ch0和ch是與控制器參數(shù)相關(guān)的常數(shù)。

    證明.考慮李雅普諾夫函數(shù)

    (65)

    利用式(40)、(53)、(60)、(61)對式(65)求導(dǎo)可得:

    (66)

    由引理1可知,下列不等式成立:

    (67)

    將式(67)代入式(66)整理可得:

    0.2785ξαδα+0.2785ξqδq

    (68)

    根據(jù):

    (69)

    (70)

    將式(69)和式(70)代入式(68)整理可得:

    0.2785ξαδα+0.2785ξqδq

    (71)

    (72)

    0.2785ξγδγ+0.2785ξαδα+0.2785ξqδq

    (73)

    利用式(72)和式(73),式(71)可整理為:

    (74)

    對式(72)求解可得:

    (75)

    由式(75)可得V2(t)≤V2(0),即V2(t)是有界的,聯(lián)合屏障李雅普諾夫函數(shù)式(39)、(47)和(59),可得不等式如下:

    (76)

    (77)

    (78)

    進一步可得:

    (79)

    (80)

    (81)

    定理2證畢。

    3 仿真分析

    為了定量分析本文所設(shè)計控制策略的性能,對非線性運動方程(1)和氣動模型(2)~(3)進行仿真,參考文獻[2]中氣動參數(shù),見表1。

    表1 高超聲速飛行器參數(shù)

    3.1 執(zhí)行器無故障的仿真分析

    高超聲速飛行器期望速度指令為Vd=4650.3 m/s,期望高度指令為hd=34328 m,控制參數(shù)選取如下:AzV=3、kV1=0.54、ξV=0.01、λV1=0.02、λV 2=0.02、kh1=0.1、kγ1=0.5、ξγ=ξq=ξα=0.01、λγ=λα=λq1=λq2=0.02、τ1=τ2=0.02、kα1=0.3、kq1=0.46和Azγ=Azα=Azq=0.15。則仿真結(jié)果如圖1~圖5所示。

    圖1 速度跟蹤曲線

    圖2 高度跟蹤曲線

    圖3 控制輸入φc和δe曲線

    圖4 狀態(tài)變量γ,α和q曲線

    圖5 彈性模態(tài)η1和η2曲線

    從圖1~圖2分別給出了速度和高度的跟蹤曲線仿真結(jié)果可知,即使存在外界干擾、模型參數(shù)不確定性和輸入飽和情形下,速度誤差和高度誤差在有限時間內(nèi)收斂到零的附近區(qū)域,能夠滿足跟蹤性能。從圖3給出控制輸入曲線可知,控制輸入在整個控制過程中是有界的。由圖4給出飛行器的其他狀態(tài)曲線可知,航跡角、攻角和俯仰率在短時間內(nèi)趨于穩(wěn)態(tài)值且滿足預(yù)期的狀態(tài)約束條件。圖5描述了彈性模型的曲線,表明彈性模態(tài)變量η1和η2在經(jīng)歷較短暫態(tài)后趨于穩(wěn)定值。

    3.2 帶有執(zhí)行器故障的仿真分析

    為驗證所設(shè)計控制器對執(zhí)行器發(fā)生故障時的有效性,控制參數(shù)和參考信號與3.1節(jié)相同,在仿真中,假設(shè)故障形式如式(82)所示,其仿真結(jié)果如圖6~圖11所示。

    圖6 速度跟蹤曲線

    圖7 高度跟蹤曲線

    圖8 控制輸入φc和δe曲線

    圖9 狀態(tài)變量γ和α曲線

    圖10 俯仰角速率q曲線

    圖11 彈性模態(tài)η1和η2曲線

    (82)

    圖6~圖7分別給出了兩個狀態(tài)的跟蹤結(jié)果,從圖中可以看出,當(dāng)飛行器發(fā)生執(zhí)行機構(gòu)失效故障時,所設(shè)計的容錯控制器能夠快速自動調(diào)整控制增益,以實現(xiàn)對故障影響的有效處理,使得速度跟蹤誤差與高度跟蹤誤差均在較短時間內(nèi)收斂到平衡點,且其精度達到了控制系統(tǒng)的性能要求。圖8所產(chǎn)生的控制變化曲線,從圖中可以看出,輸入曲線變化平緩,并能夠滿足執(zhí)行器物理受限的要求。從圖9~圖10所示的高超聲速飛行器其它狀態(tài)量的變化曲線可知,當(dāng)控制器發(fā)生故障時,由于控制器對故障進行有效處理,使得較短的時間內(nèi)航跡角γ,攻角α和俯仰角θ迅速趨于穩(wěn)定,且各個誤差變量zγ、zθ和zq分別快速的收斂至|zγ|≤0.15、|zθ|≤0.15和|zq|≤0.15。圖11為執(zhí)行機構(gòu)故障下的彈性狀態(tài)變化曲線,彈性模態(tài)變量η1和η2在較短時間內(nèi)趨于穩(wěn)態(tài)值,并能保持在一定范圍內(nèi)變化。

    4 結(jié)論

    本文針對帶有狀態(tài)約束的高超聲速飛行器跟蹤控制問題進行了深入的研究,同時考慮到彈性耦合、執(zhí)行器輸入受限和執(zhí)行器故障的工程實際需求,在構(gòu)造一種新型屏障李雅普諾夫函數(shù)基礎(chǔ)上,利用反步法、雙曲正切函數(shù)、一階濾波器和自適應(yīng)控制技術(shù),提出了一種自適應(yīng)抗飽和故障容錯控制策略,通過屏障李雅普諾夫函數(shù)的有界性,使得所設(shè)計的控制策略在滿足飛行狀態(tài)約束的同時具有良好的控制性能。借助李雅普諾夫函數(shù)對所設(shè)計的控制器進行了穩(wěn)定性證明,并通過仿真分析進一步驗證了所設(shè)計控制策略的性能。

    猜你喜歡
    李雅普諾夫超聲速
    李雅普諾夫:彼得堡數(shù)學(xué)學(xué)派的健將
    基于增廣Lyapunov 泛函的時變時滯T-S模糊系統(tǒng)穩(wěn)定性分析
    高超聲速出版工程
    高超聲速飛行器
    系統(tǒng)H∞范數(shù)計算:Lyapunov函數(shù)的直接優(yōu)化方法
    超聲速旅行
    Novoselova瑪莎·諾夫斯洛娃【英國】
    博客天下(2014年24期)2014-09-25 02:47:29
    高超聲速大博弈
    太空探索(2014年5期)2014-07-12 09:53:28
    采用李雅普諾夫函數(shù)的電液伺服系統(tǒng)反饋線性化控制
    幽默與漫畫
    海外文摘(2009年11期)2009-12-31 08:33:26
    51国产日韩欧美| 777米奇影视久久| 国产日韩欧美视频二区| 亚洲四区av| 亚洲欧美清纯卡通| 建设人人有责人人尽责人人享有的| 免费黄色在线免费观看| 亚洲婷婷狠狠爱综合网| 最近2019中文字幕mv第一页| 久久精品久久久久久久性| 美女福利国产在线| 王馨瑶露胸无遮挡在线观看| 如日韩欧美国产精品一区二区三区| 亚洲中文av在线| 日本黄色日本黄色录像| 伦理电影大哥的女人| 肉色欧美久久久久久久蜜桃| 日本与韩国留学比较| 国产亚洲欧美精品永久| 亚洲精品国产色婷婷电影| 国产男女内射视频| 如日韩欧美国产精品一区二区三区| 高清黄色对白视频在线免费看| 日本黄色日本黄色录像| 黑丝袜美女国产一区| 少妇高潮的动态图| 99精国产麻豆久久婷婷| 插逼视频在线观看| 亚洲激情五月婷婷啪啪| 国产精品成人在线| 国产女主播在线喷水免费视频网站| 久久国产精品大桥未久av| 国产成人精品久久久久久| 久久人妻熟女aⅴ| 精品国产露脸久久av麻豆| 久久亚洲国产成人精品v| 91精品三级在线观看| 国产白丝娇喘喷水9色精品| 大香蕉久久网| 99国产精品免费福利视频| 国产不卡av网站在线观看| 午夜福利在线观看免费完整高清在| av线在线观看网站| 日韩精品免费视频一区二区三区 | 麻豆乱淫一区二区| 久久精品国产a三级三级三级| 男人舔女人的私密视频| 九色亚洲精品在线播放| 免费黄网站久久成人精品| 国产色婷婷99| 色5月婷婷丁香| 国产精品一区二区在线观看99| 国产av码专区亚洲av| 亚洲av日韩在线播放| 国产精品一国产av| 国产精品成人在线| 国产欧美日韩综合在线一区二区| 最近最新中文字幕免费大全7| 亚洲欧洲日产国产| 免费观看性生交大片5| 我要看黄色一级片免费的| 99热6这里只有精品| 韩国精品一区二区三区 | 亚洲精品色激情综合| 精品第一国产精品| 欧美日韩国产mv在线观看视频| 人妻少妇偷人精品九色| 久久久精品区二区三区| 久久狼人影院| 美女大奶头黄色视频| 一级a做视频免费观看| 丝袜脚勾引网站| 亚洲精品久久成人aⅴ小说| 纵有疾风起免费观看全集完整版| 搡老乐熟女国产| 久久精品夜色国产| 91精品三级在线观看| 日本午夜av视频| av国产精品久久久久影院| 三上悠亚av全集在线观看| 女性生殖器流出的白浆| 人人妻人人澡人人看| 卡戴珊不雅视频在线播放| 中文字幕免费在线视频6| 欧美97在线视频| 久久人妻熟女aⅴ| 又粗又硬又长又爽又黄的视频| 成年人午夜在线观看视频| 久久国产精品男人的天堂亚洲 | 国产日韩欧美视频二区| 国产白丝娇喘喷水9色精品| 国产精品秋霞免费鲁丝片| 丝袜美足系列| 精品国产一区二区三区四区第35| 亚洲精品视频女| 中国国产av一级| a级毛色黄片| 人妻一区二区av| 大香蕉久久成人网| 久久 成人 亚洲| 久久久久久久国产电影| 亚洲中文av在线| 草草在线视频免费看| 精品国产国语对白av| 国产亚洲精品第一综合不卡 | 亚洲精品第二区| 亚洲国产精品999| 街头女战士在线观看网站| 2018国产大陆天天弄谢| 99热6这里只有精品| 久久av网站| 一区二区三区精品91| 亚洲国产精品一区三区| 欧美日韩亚洲高清精品| 免费大片黄手机在线观看| 日韩中字成人| 丰满少妇做爰视频| 久热这里只有精品99| 国产永久视频网站| 久久午夜综合久久蜜桃| 天堂俺去俺来也www色官网| 久久亚洲国产成人精品v| 制服诱惑二区| 美女内射精品一级片tv| 国产xxxxx性猛交| 高清黄色对白视频在线免费看| 日韩av免费高清视频| 国产精品久久久久久精品电影小说| h视频一区二区三区| av网站免费在线观看视频| 色哟哟·www| 亚洲人与动物交配视频| 看免费成人av毛片| 我的女老师完整版在线观看| 一二三四在线观看免费中文在 | 最近最新中文字幕免费大全7| 我要看黄色一级片免费的| 精品一品国产午夜福利视频| 日本与韩国留学比较| 免费观看无遮挡的男女| av在线app专区| 免费观看在线日韩| 制服人妻中文乱码| 新久久久久国产一级毛片| 亚洲一码二码三码区别大吗| 精品一区在线观看国产| 一级a做视频免费观看| 日本爱情动作片www.在线观看| 男人舔女人的私密视频| 美女大奶头黄色视频| 国产精品一区二区在线不卡| 午夜精品国产一区二区电影| 久久99一区二区三区| 激情五月婷婷亚洲| 自线自在国产av| 一级毛片黄色毛片免费观看视频| 中国国产av一级| 亚洲精品国产av蜜桃| 欧美日韩综合久久久久久| 国产男人的电影天堂91| 国产精品99久久99久久久不卡 | 日韩电影二区| 只有这里有精品99| 夜夜骑夜夜射夜夜干| 日韩视频在线欧美| 哪个播放器可以免费观看大片| 中文字幕av电影在线播放| 黑人巨大精品欧美一区二区蜜桃 | 人妻系列 视频| 十八禁网站网址无遮挡| 亚洲精品乱久久久久久| 99热全是精品| 九九在线视频观看精品| 国国产精品蜜臀av免费| 中文字幕免费在线视频6| 熟女电影av网| 成人漫画全彩无遮挡| 久久久久久久大尺度免费视频| 人人妻人人澡人人爽人人夜夜| 亚洲成国产人片在线观看| 最近2019中文字幕mv第一页| 亚洲国产日韩一区二区| 国产激情久久老熟女| 性色av一级| 看十八女毛片水多多多| 日本猛色少妇xxxxx猛交久久| 国产乱来视频区| 在线天堂中文资源库| 亚洲国产精品专区欧美| 国产成人aa在线观看| 亚洲性久久影院| 亚洲伊人久久精品综合| 久久久久国产精品人妻一区二区| 五月开心婷婷网| 国产精品人妻久久久影院| 亚洲精品第二区| 国产欧美亚洲国产| 2018国产大陆天天弄谢| 亚洲av在线观看美女高潮| 边亲边吃奶的免费视频| 毛片一级片免费看久久久久| 巨乳人妻的诱惑在线观看| 亚洲精品乱码久久久久久按摩| 午夜免费鲁丝| 少妇精品久久久久久久| 香蕉国产在线看| av国产久精品久网站免费入址| 精品国产一区二区三区四区第35| 国产精品成人在线| 丝袜脚勾引网站| 免费看光身美女| 色吧在线观看| a级毛色黄片| 亚洲av电影在线进入| 2018国产大陆天天弄谢| 欧美亚洲日本最大视频资源| 在线看a的网站| 午夜免费鲁丝| 久久这里有精品视频免费| 国产精品无大码| 精品久久久精品久久久| 亚洲精品美女久久久久99蜜臀 | 亚洲精品一二三| 亚洲美女视频黄频| 秋霞在线观看毛片| 肉色欧美久久久久久久蜜桃| 人人妻人人添人人爽欧美一区卜| 精品亚洲成国产av| 制服诱惑二区| 欧美 亚洲 国产 日韩一| 免费在线观看完整版高清| 亚洲国产精品专区欧美| 婷婷色av中文字幕| 国产免费一级a男人的天堂| 两个人看的免费小视频| 天天影视国产精品| 中文字幕人妻熟女乱码| 亚洲欧美一区二区三区国产| 在线天堂中文资源库| xxx大片免费视频| 亚洲精品456在线播放app| 国产精品久久久久久久电影| 中文字幕av电影在线播放| 青春草视频在线免费观看| 91aial.com中文字幕在线观看| 一本久久精品| 免费大片黄手机在线观看| 久久久国产一区二区| 女的被弄到高潮叫床怎么办| 亚洲av中文av极速乱| 美女内射精品一级片tv| 久久亚洲国产成人精品v| 日日撸夜夜添| 在线亚洲精品国产二区图片欧美| 亚洲av中文av极速乱| 久久精品国产综合久久久 | 日日摸夜夜添夜夜爱| 日韩精品免费视频一区二区三区 | xxxhd国产人妻xxx| 天天操日日干夜夜撸| 只有这里有精品99| 成人18禁高潮啪啪吃奶动态图| 国产亚洲一区二区精品| 五月开心婷婷网| 久久精品久久久久久噜噜老黄| 日韩欧美精品免费久久| 日本欧美视频一区| 久久人妻熟女aⅴ| 国产片内射在线| 人妻系列 视频| 日韩成人伦理影院| 亚洲精品久久久久久婷婷小说| 国产精品人妻久久久影院| 久久99热6这里只有精品| 人体艺术视频欧美日本| 久久97久久精品| 久久狼人影院| 我的女老师完整版在线观看| 18禁观看日本| 视频在线观看一区二区三区| 国产精品嫩草影院av在线观看| 日韩熟女老妇一区二区性免费视频| 免费久久久久久久精品成人欧美视频 | 啦啦啦视频在线资源免费观看| 久久精品熟女亚洲av麻豆精品| 少妇被粗大猛烈的视频| 欧美成人午夜精品| 亚洲欧美精品自产自拍| 一二三四中文在线观看免费高清| 赤兔流量卡办理| 日韩三级伦理在线观看| 精品99又大又爽又粗少妇毛片| 汤姆久久久久久久影院中文字幕| 在线 av 中文字幕| 久久久国产精品麻豆| 91久久精品国产一区二区三区| 精品人妻一区二区三区麻豆| 视频区图区小说| 国产乱来视频区| 寂寞人妻少妇视频99o| xxxhd国产人妻xxx| 在线观看一区二区三区激情| 永久网站在线| 一本色道久久久久久精品综合| 中文乱码字字幕精品一区二区三区| 午夜免费男女啪啪视频观看| av播播在线观看一区| 日韩免费高清中文字幕av| 又黄又爽又刺激的免费视频.| 国产男女超爽视频在线观看| 男女高潮啪啪啪动态图| 久久精品aⅴ一区二区三区四区 | 男女国产视频网站| 亚洲欧美色中文字幕在线| 日韩不卡一区二区三区视频在线| 国产无遮挡羞羞视频在线观看| 在线看a的网站| 大片免费播放器 马上看| 国产毛片在线视频| 男女国产视频网站| 日韩欧美一区视频在线观看| 久久久国产一区二区| 老女人水多毛片| 欧美人与性动交α欧美精品济南到 | 91精品伊人久久大香线蕉| 亚洲欧美成人精品一区二区| 一边亲一边摸免费视频| videos熟女内射| 亚洲美女黄色视频免费看| 国产精品久久久久久精品电影小说| 国产男女超爽视频在线观看| 亚洲欧美日韩卡通动漫| 日本黄色日本黄色录像| 中国国产av一级| 成年动漫av网址| 精品久久蜜臀av无| 久久国产精品男人的天堂亚洲 | 亚洲精品一二三| 亚洲欧洲日产国产| 亚洲丝袜综合中文字幕| 国产亚洲午夜精品一区二区久久| 久久婷婷青草| 日韩大片免费观看网站| 91成人精品电影| 日韩不卡一区二区三区视频在线| 欧美成人午夜免费资源| 久久这里只有精品19| 香蕉精品网在线| 亚洲精品,欧美精品| a级毛片黄视频| 香蕉丝袜av| 丰满少妇做爰视频| 中文字幕另类日韩欧美亚洲嫩草| videossex国产| 少妇高潮的动态图| 黄色毛片三级朝国网站| 岛国毛片在线播放| 丝袜美足系列| 国产精品一区二区在线观看99| 91在线精品国自产拍蜜月| kizo精华| 天堂8中文在线网| 18禁国产床啪视频网站| 中国国产av一级| 午夜视频国产福利| 日本猛色少妇xxxxx猛交久久| 国产免费又黄又爽又色| 亚洲国产欧美在线一区| 亚洲丝袜综合中文字幕| 国产激情久久老熟女| 少妇精品久久久久久久| 国产成人精品婷婷| 免费人成在线观看视频色| 久久久久视频综合| 天天影视国产精品| 欧美精品国产亚洲| 男人添女人高潮全过程视频| 日韩av免费高清视频| 久久狼人影院| 亚洲欧洲日产国产| 国产av一区二区精品久久| 啦啦啦在线观看免费高清www| 久久久国产欧美日韩av| 桃花免费在线播放| 欧美最新免费一区二区三区| 校园人妻丝袜中文字幕| 人妻一区二区av| 中国国产av一级| 免费看不卡的av| 丰满乱子伦码专区| 少妇被粗大的猛进出69影院 | 97在线人人人人妻| 亚洲精品色激情综合| 国产色爽女视频免费观看| 爱豆传媒免费全集在线观看| 亚洲第一av免费看| 免费久久久久久久精品成人欧美视频 | 欧美日本中文国产一区发布| 国产av国产精品国产| 成人免费观看视频高清| 国产免费一区二区三区四区乱码| 你懂的网址亚洲精品在线观看| 又黄又粗又硬又大视频| 日本黄大片高清| 精品一区二区三区视频在线| 午夜福利网站1000一区二区三区| 久久久久国产精品人妻一区二区| 大香蕉久久成人网| 欧美日韩亚洲高清精品| 菩萨蛮人人尽说江南好唐韦庄| 在线观看人妻少妇| 亚洲av福利一区| 成年女人在线观看亚洲视频| 亚洲精品国产色婷婷电影| 熟女av电影| 色视频在线一区二区三区| 看免费成人av毛片| 午夜老司机福利剧场| 99视频精品全部免费 在线| 美国免费a级毛片| 午夜福利在线观看免费完整高清在| 色婷婷久久久亚洲欧美| 美女大奶头黄色视频| 一边摸一边做爽爽视频免费| 日本午夜av视频| 亚洲美女搞黄在线观看| 十分钟在线观看高清视频www| 国产高清三级在线| av一本久久久久| 亚洲性久久影院| 国产av国产精品国产| 久久精品人人爽人人爽视色| 国产精品久久久久成人av| 亚洲欧美精品自产自拍| 精品一区二区三区视频在线| 母亲3免费完整高清在线观看 | 欧美成人精品欧美一级黄| 人人妻人人澡人人看| 伊人久久国产一区二区| av又黄又爽大尺度在线免费看| 国产精品久久久av美女十八| 国产精品三级大全| 亚洲国产精品999| 校园人妻丝袜中文字幕| 成人亚洲欧美一区二区av| 99久国产av精品国产电影| 国产精品 国内视频| 黄网站色视频无遮挡免费观看| av免费观看日本| 国产在线一区二区三区精| 中国三级夫妇交换| 成年动漫av网址| 亚洲精品国产色婷婷电影| 亚洲国产精品999| xxxhd国产人妻xxx| 欧美最新免费一区二区三区| 日本欧美国产在线视频| 男女啪啪激烈高潮av片| 天堂俺去俺来也www色官网| 看非洲黑人一级黄片| 国产亚洲精品第一综合不卡 | 热99久久久久精品小说推荐| 亚洲 欧美一区二区三区| 黑丝袜美女国产一区| av电影中文网址| 国产成人精品无人区| 精品酒店卫生间| 日日摸夜夜添夜夜爱| 亚洲人成77777在线视频| 一级毛片我不卡| 成年美女黄网站色视频大全免费| 9热在线视频观看99| 国产极品天堂在线| 自线自在国产av| 精品午夜福利在线看| 午夜免费观看性视频| 丰满迷人的少妇在线观看| 十八禁网站网址无遮挡| 一本久久精品| 伊人亚洲综合成人网| 久久99热6这里只有精品| 日韩中字成人| 老女人水多毛片| 国产激情久久老熟女| 熟女av电影| 两个人看的免费小视频| 少妇高潮的动态图| 精品久久久精品久久久| 1024视频免费在线观看| 女的被弄到高潮叫床怎么办| 成人国语在线视频| 在线观看免费视频网站a站| 亚洲国产精品一区三区| 久久精品人人爽人人爽视色| 亚洲经典国产精华液单| 久久久亚洲精品成人影院| 精品少妇久久久久久888优播| 久久婷婷青草| 日本vs欧美在线观看视频| 久久精品国产亚洲av天美| 亚洲欧美成人综合另类久久久| 国产一级毛片在线| 交换朋友夫妻互换小说| 精品一区二区三区视频在线| 日本黄大片高清| 一本大道久久a久久精品| 插逼视频在线观看| 王馨瑶露胸无遮挡在线观看| 日本猛色少妇xxxxx猛交久久| 亚洲久久久国产精品| 婷婷色麻豆天堂久久| 亚洲综合精品二区| 色94色欧美一区二区| 一二三四中文在线观看免费高清| 久久久久精品人妻al黑| 男女国产视频网站| 精品视频人人做人人爽| 国产成人精品婷婷| 亚洲少妇的诱惑av| 亚洲欧洲国产日韩| 欧美成人午夜免费资源| 亚洲美女黄色视频免费看| 韩国高清视频一区二区三区| 制服诱惑二区| 欧美亚洲日本最大视频资源| 日韩中文字幕视频在线看片| 91精品三级在线观看| 2021少妇久久久久久久久久久| 街头女战士在线观看网站| 午夜福利视频在线观看免费| 丰满迷人的少妇在线观看| 精品国产露脸久久av麻豆| 黄网站色视频无遮挡免费观看| 精品一区二区免费观看| 宅男免费午夜| 最后的刺客免费高清国语| 日韩av在线免费看完整版不卡| 国产深夜福利视频在线观看| 男女高潮啪啪啪动态图| 成人亚洲欧美一区二区av| 人妻系列 视频| www.色视频.com| kizo精华| 婷婷色综合www| 国产黄频视频在线观看| 高清在线视频一区二区三区| 欧美少妇被猛烈插入视频| 看十八女毛片水多多多| 在线 av 中文字幕| av网站免费在线观看视频| 久久久久久久国产电影| 久久精品国产亚洲av涩爱| 高清不卡的av网站| 成人免费观看视频高清| 欧美国产精品va在线观看不卡| 欧美人与性动交α欧美精品济南到 | 国产女主播在线喷水免费视频网站| 亚洲国产色片| 亚洲欧美成人综合另类久久久| 久久精品国产鲁丝片午夜精品| 中文字幕最新亚洲高清| 欧美日本中文国产一区发布| 亚洲精品日本国产第一区| 伊人亚洲综合成人网| 老熟女久久久| 亚洲一级一片aⅴ在线观看| 午夜福利,免费看| 亚洲精品久久午夜乱码| 在线看a的网站| 在线免费观看不下载黄p国产| 亚洲欧洲国产日韩| 97超碰精品成人国产| 久久97久久精品| 久久久亚洲精品成人影院| 久久精品夜色国产| av播播在线观看一区| 老熟女久久久| 97在线视频观看| 国产精品国产三级国产专区5o| 中文字幕另类日韩欧美亚洲嫩草| 久久久久精品性色| 汤姆久久久久久久影院中文字幕| 亚洲国产毛片av蜜桃av| 亚洲 欧美一区二区三区| 美女脱内裤让男人舔精品视频| 欧美丝袜亚洲另类| 国产极品粉嫩免费观看在线| kizo精华| 少妇猛男粗大的猛烈进出视频| 国产永久视频网站| 欧美xxxx性猛交bbbb| 国产麻豆69| av黄色大香蕉| 久久久久视频综合| 亚洲中文av在线| 99久久精品国产国产毛片| 国产欧美日韩综合在线一区二区| 天天躁夜夜躁狠狠躁躁| av天堂久久9| 大话2 男鬼变身卡| 国产精品国产三级国产专区5o| 精品午夜福利在线看| 熟女av电影| av不卡在线播放| 成人毛片a级毛片在线播放| 丝瓜视频免费看黄片| 亚洲国产精品一区三区| 五月开心婷婷网| 天堂俺去俺来也www色官网| 久久这里只有精品19| 欧美激情国产日韩精品一区| 国产成人aa在线观看| 婷婷色综合大香蕉| 日韩欧美一区视频在线观看| 久久青草综合色|