• <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
    91九色精品人成在线观看| 一进一出抽搐gif免费好疼| 久久人妻av系列| 亚洲最大成人手机在线| 免费人成视频x8x8入口观看| 国产成人福利小说| 国产午夜福利久久久久久| 小说图片视频综合网站| 一卡2卡三卡四卡精品乱码亚洲| 十八禁人妻一区二区| 国产三级黄色录像| 无遮挡黄片免费观看| 成人av一区二区三区在线看| 亚洲av中文字字幕乱码综合| 可以在线观看毛片的网站| 怎么达到女性高潮| 久久久久免费精品人妻一区二区| 欧美乱色亚洲激情| 99精品在免费线老司机午夜| 亚洲欧美精品综合久久99| 一本久久中文字幕| 最近在线观看免费完整版| 欧美性猛交╳xxx乱大交人| 欧美黑人欧美精品刺激| 波多野结衣高清作品| 午夜福利免费观看在线| 偷拍熟女少妇极品色| 丝袜美腿在线中文| 一区二区三区免费毛片| 免费观看的影片在线观看| 日韩 亚洲 欧美在线| 国产成人欧美在线观看| 久久精品国产亚洲av天美| 国产又黄又爽又无遮挡在线| 窝窝影院91人妻| 欧美绝顶高潮抽搐喷水| av天堂中文字幕网| 18禁黄网站禁片免费观看直播| 久久久精品欧美日韩精品| 精品熟女少妇八av免费久了| 看黄色毛片网站| 国产免费一级a男人的天堂| 国产精品一及| 精品欧美国产一区二区三| av黄色大香蕉| 亚洲色图av天堂| 亚洲精品成人久久久久久| 黄色视频,在线免费观看| 18禁黄网站禁片午夜丰满| 亚洲三级黄色毛片| 日本免费一区二区三区高清不卡| 亚洲性夜色夜夜综合| 国产成人av教育| 窝窝影院91人妻| 免费在线观看影片大全网站| 久久久久久久午夜电影| 国产精品美女特级片免费视频播放器| 色综合婷婷激情| 久久久久久九九精品二区国产| 亚洲国产欧洲综合997久久,| 狂野欧美白嫩少妇大欣赏| 免费人成在线观看视频色| 1000部很黄的大片| 亚洲精品日韩av片在线观看| 婷婷色综合大香蕉| 久久久精品大字幕| 看片在线看免费视频| 国产白丝娇喘喷水9色精品| 麻豆成人av在线观看| 亚洲中文字幕一区二区三区有码在线看| 永久网站在线| 三级国产精品欧美在线观看| 日本黄色视频三级网站网址| 午夜免费成人在线视频| 国产蜜桃级精品一区二区三区| 国产黄a三级三级三级人| 18禁在线播放成人免费| 一级毛片久久久久久久久女| 免费高清视频大片| 一二三四社区在线视频社区8| 在线免费观看的www视频| 一夜夜www| 欧美日韩福利视频一区二区| 99久久99久久久精品蜜桃| 亚洲av成人不卡在线观看播放网| 一级毛片久久久久久久久女| 日韩欧美在线二视频| 婷婷亚洲欧美| 一区福利在线观看| 欧美xxxx黑人xx丫x性爽| 直男gayav资源| 精品人妻1区二区| 亚洲熟妇熟女久久| 国产中年淑女户外野战色| 在线国产一区二区在线| 精品久久久久久,| 亚洲av一区综合| 免费在线观看成人毛片| 成人高潮视频无遮挡免费网站| 中文字幕人成人乱码亚洲影| 18美女黄网站色大片免费观看| 亚洲久久久久久中文字幕| 国产精品99久久久久久久久| 久久欧美精品欧美久久欧美| 欧美午夜高清在线| 搞女人的毛片| 国产精品不卡视频一区二区 | a级毛片免费高清观看在线播放| 亚洲成av人片免费观看| 国产单亲对白刺激| 国产aⅴ精品一区二区三区波| 色综合婷婷激情| 午夜福利在线观看免费完整高清在 | 欧美高清成人免费视频www| 国产午夜福利久久久久久| 色综合亚洲欧美另类图片| 日韩av在线大香蕉| 欧美黄色淫秽网站| 欧美潮喷喷水| 日韩欧美 国产精品| 91麻豆精品激情在线观看国产| 国产精品久久久久久人妻精品电影| 99热这里只有是精品在线观看 | 十八禁网站免费在线| 国产高清有码在线观看视频| 国产人妻一区二区三区在| 高清在线国产一区| 日本熟妇午夜| 又粗又爽又猛毛片免费看| 久久精品国产清高在天天线| 床上黄色一级片| 极品教师在线免费播放| 禁无遮挡网站| 成人亚洲精品av一区二区| 18+在线观看网站| 欧美日韩黄片免| 欧美在线黄色| 亚洲av日韩精品久久久久久密| 高清毛片免费观看视频网站| 首页视频小说图片口味搜索| 可以在线观看毛片的网站| 99视频精品全部免费 在线| 91av网一区二区| 亚洲av日韩精品久久久久久密| 亚洲av免费高清在线观看| 午夜久久久久精精品| 国产熟女xx| 大型黄色视频在线免费观看| 国产伦人伦偷精品视频| 精品久久久久久成人av| 成人特级黄色片久久久久久久| 精品午夜福利在线看| 亚洲五月婷婷丁香| 天堂av国产一区二区熟女人妻| 亚洲真实伦在线观看| 亚洲七黄色美女视频| 亚洲第一电影网av| 久久热精品热| 欧美成人性av电影在线观看| 亚洲在线观看片| 又粗又爽又猛毛片免费看| 少妇裸体淫交视频免费看高清| 别揉我奶头 嗯啊视频| 欧美日韩黄片免| 亚洲性夜色夜夜综合| 日韩欧美在线二视频| 国产精品一区二区性色av| 黄色视频,在线免费观看| 亚洲国产精品久久男人天堂| 老熟妇仑乱视频hdxx| 日本免费a在线| 亚洲av电影不卡..在线观看| 中文字幕av成人在线电影| 91在线精品国自产拍蜜月| 国产综合懂色| 日本五十路高清| av视频在线观看入口| 国产精品一区二区免费欧美| 18禁黄网站禁片免费观看直播| 国产亚洲av嫩草精品影院| 18+在线观看网站| 热99在线观看视频| 能在线免费观看的黄片| 亚洲最大成人中文| 国产一区二区激情短视频| 人妻久久中文字幕网| 欧美日韩黄片免| 少妇裸体淫交视频免费看高清| 老女人水多毛片| 日本五十路高清| 在线观看午夜福利视频| 特级一级黄色大片| 亚洲熟妇熟女久久| 亚洲av不卡在线观看| 麻豆久久精品国产亚洲av| 免费看美女性在线毛片视频| 亚洲三级黄色毛片| 丰满的人妻完整版| 国产成人欧美在线观看| 午夜福利高清视频| 久久久久久久精品吃奶| 国产精品不卡视频一区二区 | 日韩大尺度精品在线看网址| 五月伊人婷婷丁香| 一本一本综合久久| 欧美日韩国产亚洲二区| 三级男女做爰猛烈吃奶摸视频| 免费观看的影片在线观看| 成人无遮挡网站| 亚洲av电影在线进入| 在线观看午夜福利视频| 99精品久久久久人妻精品| 99久国产av精品| 99久久99久久久精品蜜桃| 中文字幕免费在线视频6| 亚洲中文日韩欧美视频| 欧美绝顶高潮抽搐喷水| 久久香蕉精品热| 欧美在线黄色| 内射极品少妇av片p| 激情在线观看视频在线高清| 在线观看一区二区三区| www.www免费av| 亚洲精品乱码久久久v下载方式| 级片在线观看| 国产伦一二天堂av在线观看| 国产精品一区二区性色av| 欧美一区二区精品小视频在线| 久久国产乱子伦精品免费另类| 亚洲第一区二区三区不卡| 欧美在线黄色| 亚洲18禁久久av| 搡老妇女老女人老熟妇| 午夜精品在线福利| 精品无人区乱码1区二区| 亚洲不卡免费看| 搡老熟女国产l中国老女人| 夜夜看夜夜爽夜夜摸| 午夜免费成人在线视频| 国产精品免费一区二区三区在线| 日本a在线网址| 在线看三级毛片| 欧美日韩瑟瑟在线播放| 国产一区二区三区在线臀色熟女| 最新在线观看一区二区三区| 国产免费一级a男人的天堂| 国产精品爽爽va在线观看网站| 国内久久婷婷六月综合欲色啪| 精品福利观看| 免费av毛片视频| 久久久色成人| 中文字幕久久专区| 免费av毛片视频| 午夜免费激情av| 成人高潮视频无遮挡免费网站| 丰满的人妻完整版| 俄罗斯特黄特色一大片| 亚洲天堂国产精品一区在线| 国产男靠女视频免费网站| 天美传媒精品一区二区| 又爽又黄a免费视频| 国产午夜精品久久久久久一区二区三区 | 成人一区二区视频在线观看| 日韩中字成人| 一本久久中文字幕| 欧美精品啪啪一区二区三区| 国产国拍精品亚洲av在线观看| 99热这里只有精品一区| 我的老师免费观看完整版| 国产亚洲欧美在线一区二区| 日本三级黄在线观看| 99国产极品粉嫩在线观看| 丝袜美腿在线中文| 日日夜夜操网爽| 在线a可以看的网站| 国产黄a三级三级三级人| 国产单亲对白刺激| 身体一侧抽搐| 在线观看美女被高潮喷水网站 | 亚洲中文字幕日韩| 欧美一区二区国产精品久久精品| 99热这里只有精品一区| 嫩草影视91久久| 亚洲最大成人手机在线| 亚洲成av人片在线播放无| 精品国产三级普通话版| 国产精品嫩草影院av在线观看 | 国产高清激情床上av| 在线观看午夜福利视频| 国产色婷婷99| 不卡一级毛片| 久久久久久久精品吃奶| 国产欧美日韩精品一区二区| 美女大奶头视频| 91久久精品国产一区二区成人| 18+在线观看网站| 男女下面进入的视频免费午夜| 国产亚洲精品久久久com| 欧美中文日本在线观看视频| 国产爱豆传媒在线观看| 最近中文字幕高清免费大全6 | 一个人看的www免费观看视频| 久久婷婷人人爽人人干人人爱| 久久久国产成人免费| 好看av亚洲va欧美ⅴa在| 国产av不卡久久| 国产视频内射| 熟女电影av网| 少妇的逼好多水| 久久久久国内视频| 国产免费男女视频| 色吧在线观看| 夜夜夜夜夜久久久久| 在线观看av片永久免费下载| 国产精品久久久久久人妻精品电影| 夜夜躁狠狠躁天天躁| 3wmmmm亚洲av在线观看| 18禁裸乳无遮挡免费网站照片| 日日摸夜夜添夜夜添小说| 亚洲精品成人久久久久久| 国产精品嫩草影院av在线观看 | 中文在线观看免费www的网站| 国产高潮美女av| 久久午夜亚洲精品久久| 伊人久久精品亚洲午夜| 国产精品免费一区二区三区在线| 婷婷丁香在线五月| 久久精品国产亚洲av涩爱 | 国产蜜桃级精品一区二区三区| 亚洲国产欧美人成| 免费搜索国产男女视频| 少妇人妻精品综合一区二区 | 久久久久久久久久成人| 在线观看免费视频日本深夜| 2021天堂中文幕一二区在线观| 国产91精品成人一区二区三区| 亚洲最大成人av| 真人一进一出gif抽搐免费| 久久久久精品国产欧美久久久| 午夜精品久久久久久毛片777| 97超级碰碰碰精品色视频在线观看| 啦啦啦韩国在线观看视频| 欧美乱色亚洲激情| 在线观看66精品国产| 宅男免费午夜| 国产淫片久久久久久久久 | 美女被艹到高潮喷水动态| 婷婷亚洲欧美| 黄色日韩在线| 国产乱人视频| 国产精品,欧美在线| 成人午夜高清在线视频| 校园春色视频在线观看| 欧美日韩国产亚洲二区| 亚洲精品在线观看二区| 嫩草影院精品99| 国产视频内射| 亚洲av中文字字幕乱码综合| 动漫黄色视频在线观看| 色哟哟·www| 最近在线观看免费完整版| 少妇的逼水好多| 欧美区成人在线视频| 内射极品少妇av片p| 91av网一区二区| 观看美女的网站| 99热6这里只有精品| 中文字幕高清在线视频| 欧美bdsm另类| 国产中年淑女户外野战色| 日韩人妻高清精品专区| 一区二区三区四区激情视频 | 日韩人妻高清精品专区| 午夜老司机福利剧场| 国产伦精品一区二区三区四那| 久久精品国产自在天天线| 国产美女午夜福利| 在线观看免费视频日本深夜| 欧美成人性av电影在线观看| 精品国内亚洲2022精品成人| 亚洲熟妇熟女久久| 国产av一区在线观看免费| 国产三级中文精品| 一本久久中文字幕| 一级黄色大片毛片| 在线观看美女被高潮喷水网站 | 大型黄色视频在线免费观看| 亚洲熟妇熟女久久| 精品久久久久久久人妻蜜臀av| 美女大奶头视频| 69av精品久久久久久| 免费在线观看亚洲国产| 精品久久久久久久末码| 可以在线观看的亚洲视频| 亚洲国产精品久久男人天堂| 欧美三级亚洲精品| 草草在线视频免费看| 欧美在线一区亚洲| 蜜桃亚洲精品一区二区三区| a级一级毛片免费在线观看| 两个人的视频大全免费| 床上黄色一级片| 国产中年淑女户外野战色| 国产欧美日韩精品一区二区| eeuss影院久久| 精品熟女少妇八av免费久了| 9191精品国产免费久久| 大型黄色视频在线免费观看| 床上黄色一级片| 色精品久久人妻99蜜桃| 首页视频小说图片口味搜索| 国产aⅴ精品一区二区三区波| 日韩精品中文字幕看吧| 99国产极品粉嫩在线观看| 欧美一区二区国产精品久久精品| 嫁个100分男人电影在线观看| 亚洲第一欧美日韩一区二区三区| 91狼人影院| 一级a爱片免费观看的视频| 琪琪午夜伦伦电影理论片6080| 可以在线观看毛片的网站| 美女免费视频网站| 成年人黄色毛片网站| 国产伦在线观看视频一区| 久久中文看片网| 国产精品久久电影中文字幕| 日韩人妻高清精品专区| 欧美在线一区亚洲| 在线国产一区二区在线| 我要看日韩黄色一级片| 夜夜躁狠狠躁天天躁| www日本黄色视频网| 欧美日韩乱码在线| 欧美日韩瑟瑟在线播放| 俄罗斯特黄特色一大片| 国产精品不卡视频一区二区 | 精品午夜福利视频在线观看一区| 动漫黄色视频在线观看| а√天堂www在线а√下载| 国产伦一二天堂av在线观看| 少妇的逼好多水| 人妻久久中文字幕网| 最好的美女福利视频网| 亚洲av电影不卡..在线观看| 午夜福利18| 美女 人体艺术 gogo| 精品人妻视频免费看| 午夜精品久久久久久毛片777| 在线天堂最新版资源| 亚洲最大成人av| 人妻制服诱惑在线中文字幕| 国产一区二区三区在线臀色熟女| 欧美bdsm另类| 成人性生交大片免费视频hd| 免费av不卡在线播放| 最近最新免费中文字幕在线| 少妇熟女aⅴ在线视频| 久久久久久大精品| 夜夜躁狠狠躁天天躁| 久久久久久国产a免费观看| 日日摸夜夜添夜夜添小说| 99热只有精品国产| 午夜视频国产福利| 十八禁网站免费在线| 国产在线精品亚洲第一网站| 乱人视频在线观看| 欧美高清成人免费视频www| 丰满乱子伦码专区| 亚洲欧美日韩东京热| 在线观看一区二区三区| www.www免费av| .国产精品久久| 久久久久久久精品吃奶| 婷婷亚洲欧美| 久久九九热精品免费| 亚洲精品在线美女| 日本 欧美在线| 久久久久久久久久成人| 国产91精品成人一区二区三区| 精品日产1卡2卡| 精品人妻视频免费看| 欧美色视频一区免费| 身体一侧抽搐| 国产成人影院久久av| 国产亚洲精品久久久久久毛片| 国产麻豆成人av免费视频| 国产免费一级a男人的天堂| 真人一进一出gif抽搐免费| 女同久久另类99精品国产91| 哪里可以看免费的av片| 69av精品久久久久久| 最后的刺客免费高清国语| 欧美极品一区二区三区四区| 欧美区成人在线视频| 久久久精品欧美日韩精品| 又爽又黄a免费视频| 一个人看的www免费观看视频| 在线国产一区二区在线| 欧美一区二区精品小视频在线| 日日干狠狠操夜夜爽| 丁香六月欧美| av国产免费在线观看| 麻豆成人午夜福利视频| 国内精品久久久久精免费| 人人妻,人人澡人人爽秒播| 亚洲七黄色美女视频| 精品久久久久久久久久免费视频| 国产精品久久久久久精品电影| 欧美xxxx性猛交bbbb| 日韩国内少妇激情av| 亚洲自偷自拍三级| 日韩 亚洲 欧美在线| 看十八女毛片水多多多| 久久人人爽人人爽人人片va | 欧美日韩中文字幕国产精品一区二区三区| 1024手机看黄色片| 国产欧美日韩精品一区二区| 国产乱人视频| av视频在线观看入口| 国产精品永久免费网站| 欧美日韩综合久久久久久 | 高潮久久久久久久久久久不卡| 午夜福利在线在线| 色精品久久人妻99蜜桃| 波多野结衣巨乳人妻| 97人妻精品一区二区三区麻豆| 男插女下体视频免费在线播放| 国产高清三级在线| 国产视频一区二区在线看| 亚洲国产欧美人成| 在线播放国产精品三级| 老熟妇仑乱视频hdxx| 啪啪无遮挡十八禁网站| 十八禁国产超污无遮挡网站| 午夜a级毛片| 午夜老司机福利剧场| 变态另类成人亚洲欧美熟女| 性插视频无遮挡在线免费观看| 国产三级黄色录像| 日韩国内少妇激情av| 亚洲欧美日韩高清专用| 国产不卡一卡二| 免费观看精品视频网站| 91在线精品国自产拍蜜月| 亚洲成人久久性| 如何舔出高潮| 女生性感内裤真人,穿戴方法视频| 伦理电影大哥的女人| 国产美女午夜福利| 国产伦精品一区二区三区视频9| 中文字幕av在线有码专区| 久久亚洲真实| 人人妻,人人澡人人爽秒播| 午夜久久久久精精品| 欧美一区二区亚洲| 欧美日韩黄片免| 老司机福利观看| 国产一区二区亚洲精品在线观看| 白带黄色成豆腐渣| 少妇人妻精品综合一区二区 | 噜噜噜噜噜久久久久久91| 国产精品精品国产色婷婷| 亚洲国产精品久久男人天堂| av欧美777| 天堂影院成人在线观看| 在线观看美女被高潮喷水网站 | 国产精品一区二区三区四区久久| 欧美绝顶高潮抽搐喷水| 嫩草影院入口| 此物有八面人人有两片| 美女 人体艺术 gogo| 动漫黄色视频在线观看| 特大巨黑吊av在线直播| 久久婷婷人人爽人人干人人爱| 亚洲三级黄色毛片| 久久久精品欧美日韩精品| 757午夜福利合集在线观看| 一级作爱视频免费观看| 国产高清视频在线播放一区| 极品教师在线视频| 在线观看一区二区三区| 色尼玛亚洲综合影院| or卡值多少钱| 国产亚洲精品综合一区在线观看| 国产一区二区激情短视频| 两性午夜刺激爽爽歪歪视频在线观看| 91久久精品国产一区二区成人| 看十八女毛片水多多多| 99久久精品国产亚洲精品| 2021天堂中文幕一二区在线观| 99精品久久久久人妻精品| 色在线成人网| 亚洲人成电影免费在线| 国产男靠女视频免费网站| 国内精品久久久久精免费| 淫秽高清视频在线观看| 国语自产精品视频在线第100页| 久久午夜亚洲精品久久| 中文字幕av成人在线电影| 国产综合懂色| 成人特级av手机在线观看| 成人特级黄色片久久久久久久| 热99re8久久精品国产| 久久久久久久久久黄片| 真人做人爱边吃奶动态| 亚洲五月婷婷丁香| 日韩精品青青久久久久久| 久久午夜福利片| 一个人看视频在线观看www免费| 欧美日韩中文字幕国产精品一区二区三区| 观看免费一级毛片| 动漫黄色视频在线观看| 国产伦人伦偷精品视频| 亚洲无线观看免费| 偷拍熟女少妇极品色|