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

    基于H∞理論的交會(huì)對(duì)接模擬器橫向平移控制系統(tǒng)設(shè)計(jì)*

    2009-12-12 05:29:52李貴明
    關(guān)鍵詞:魯棒交會(huì)模擬器

    李貴明,姚 郁

    (1.北京控制工程研究所,北京100190;2.空間智能控制技術(shù)國(guó)家級(jí)重點(diǎn)實(shí)驗(yàn)室,北京100190;3.哈爾濱工業(yè)大學(xué)控制與仿真中心,哈爾濱150001)

    基于H∞理論的交會(huì)對(duì)接模擬器橫向平移控制系統(tǒng)設(shè)計(jì)*

    李貴明1,2,姚 郁3

    (1.北京控制工程研究所,北京100190;2.空間智能控制技術(shù)國(guó)家級(jí)重點(diǎn)實(shí)驗(yàn)室,北京100190;3.哈爾濱工業(yè)大學(xué)控制與仿真中心,哈爾濱150001)

    研究九自由度空間交會(huì)對(duì)接模擬器橫向平移控制系統(tǒng)設(shè)計(jì)問(wèn)題.采用非結(jié)構(gòu)不確定性描述系統(tǒng)中的高頻機(jī)械諧振,給出了系統(tǒng)對(duì)象的不確定性模型.根據(jù)橫向平移系統(tǒng)的性能指標(biāo)和主要不確定性確定了魯棒H∞控制方案,并研究了混合靈敏度問(wèn)題的加權(quán)函數(shù)選取原則和確定方法,給出了基于線性矩陣不等式(LMI)的H∞控制律的存在條件和設(shè)計(jì)方法.實(shí)際系統(tǒng)測(cè)試結(jié)果達(dá)到指標(biāo)要求,說(shuō)明了該方法的有效性.

    空間交會(huì)對(duì)接模擬器;橫向平移控制;H∞控制

    空間交會(huì)對(duì)接是載人航天工程的重大課題[1].導(dǎo)航、制導(dǎo)與控制系統(tǒng)使兩飛行器在對(duì)接前的相對(duì)位置、姿態(tài)和速度都滿足一定要求,對(duì)接機(jī)構(gòu)接觸碰撞并自動(dòng)相互鎖緊,最后完成對(duì)接任務(wù).

    九自由度空間交會(huì)對(duì)接模擬器用于模擬交會(huì)過(guò)程中最后段內(nèi)的相對(duì)運(yùn)動(dòng),對(duì)接任務(wù)決定模擬器應(yīng)具備九個(gè)自由度.模擬器采用6+3的系統(tǒng)方案,其中的追蹤飛行器的六自由度由三個(gè)姿態(tài)運(yùn)動(dòng)和三個(gè)平移運(yùn)動(dòng)組成.要求模擬質(zhì)心運(yùn)動(dòng)的三個(gè)平移運(yùn)動(dòng)具有大行程和嚴(yán)格的垂直度,還要求平移系統(tǒng)滿足低速平穩(wěn)性和穩(wěn)態(tài)精度等指標(biāo).

    自從1981年Zames提出將H∞范數(shù)作為系統(tǒng)設(shè)計(jì)的工具、開(kāi)創(chuàng)H∞控制理論以來(lái),經(jīng)過(guò)二十多年的發(fā)展,H∞控制理論已進(jìn)入完善和推廣時(shí)期,是解決不確定系統(tǒng)控制問(wèn)題的一種綜合方法[2].針對(duì)機(jī)電伺服系統(tǒng),可根據(jù)系統(tǒng)的魯棒穩(wěn)定性和魯棒性能要求直接設(shè)計(jì)出H∞控制器,是一種系統(tǒng)化的設(shè)計(jì)方法.H∞控制問(wèn)題的基于LMI(linearmatrix inequality)的求解方法,需要的假設(shè)條件相對(duì)較少[3],把問(wèn)題轉(zhuǎn)化為求解三個(gè)線性矩陣不等式的凸優(yōu)化問(wèn)題,采用“內(nèi)點(diǎn)法”進(jìn)行求解.

    本文以空間交會(huì)對(duì)接模擬器橫向平移控制系統(tǒng)設(shè)計(jì)為背景,采用了魯棒H∞控制方案,研究了加權(quán)函數(shù)的選取原則,給出了基于LMI的H∞控制律存在條件和求解方法,并給出了實(shí)際系統(tǒng)的測(cè)試結(jié)果.

    1 橫向平移系統(tǒng)結(jié)構(gòu)和不確定性描述

    模擬器橫向平移系統(tǒng)結(jié)構(gòu)如圖1所示.選定哈爾濱工業(yè)大學(xué)自主研制的交流無(wú)刷力矩電機(jī),該電機(jī)設(shè)計(jì)成中空結(jié)構(gòu),絲杠從電機(jī)軸心穿過(guò),電機(jī)轉(zhuǎn)子直接帶動(dòng)螺母旋轉(zhuǎn),采用旋轉(zhuǎn)螺母滾珠絲杠傳動(dòng)為本系統(tǒng)的傳動(dòng)方案.

    圖1 橫向平移系統(tǒng)結(jié)構(gòu)

    平移系統(tǒng)的被控對(duì)象由交流電機(jī)本體、逆變驅(qū)動(dòng)電路和位置傳感器組成,通過(guò)機(jī)理分析推導(dǎo)出其傳遞函數(shù)為

    式中,U(s)為電機(jī)驅(qū)動(dòng)器的輸入電壓,θ(s)為旋轉(zhuǎn)變壓器的輸出位置,τm、τe分別為機(jī)電、電氣時(shí)間常數(shù).

    橫向平移系統(tǒng)的模型可通過(guò)輸入變頻正弦信號(hào)辨識(shí)得到.由于受系統(tǒng)自身?xiàng)l件限制,測(cè)試不能達(dá)到很高的頻段,導(dǎo)致了高頻未建模的動(dòng)態(tài).而作為機(jī)械系統(tǒng)固有特性的機(jī)械諧振在高頻區(qū)表現(xiàn)尤為突出[4].根據(jù)工程實(shí)踐經(jīng)驗(yàn),可將其表述為如下形式:

    式中,ω0為諧振峰頻率,a/b為諧振峰值.諧振特性由機(jī)械結(jié)構(gòu)決定,且難以準(zhǔn)確描述,只可根據(jù)實(shí)際測(cè)試數(shù)據(jù)確定出機(jī)械諧振引入的不確定性的界.考慮到未建模特性是非結(jié)構(gòu)不確定性,故采用乘性攝動(dòng)加以描述,即實(shí)際對(duì)象P(s)為

    式中‖Δ‖∞≤1,界函數(shù)W(s)有待于設(shè)計(jì).乘性不確定性可由圖2來(lái)描述.

    圖2 對(duì)象乘性不確定性描述框圖

    2 基于LM I的H∞控制律

    2.1 混合靈敏度問(wèn)題與標(biāo)準(zhǔn)H∞控制問(wèn)題

    考慮了不確定性、魯棒穩(wěn)定性和魯棒性能的伺服系統(tǒng)的H∞控制問(wèn)題可以表述為如圖3所示的混合靈敏度問(wèn)題[5].

    圖3 混合靈敏度設(shè)計(jì)的系統(tǒng)框圖

    圖3中r、e、u、d和y分別為參考輸入、跟蹤誤差、控制輸入、干擾輸入和系統(tǒng)輸出.從r至e,u和y的閉環(huán)傳遞函數(shù)分別為

    式中,L(s)=P0(s)K(s)為系統(tǒng)開(kāi)環(huán)傳遞函數(shù).通常S(s)和T(s)分別稱為靈敏度函數(shù)和補(bǔ)靈敏度函數(shù).

    控制器設(shè)計(jì)的目標(biāo)為通過(guò)最小化r到z1的傳遞函數(shù)來(lái)滿足跟蹤和抑制干擾的性能;最小化r到z2的傳遞函數(shù)來(lái)防止系統(tǒng)飽和、調(diào)整系統(tǒng)帶寬;最小化r到z3的傳遞函數(shù)來(lái)滿足魯棒穩(wěn)定性,即

    式(7)~式(9)是考慮不確定性上限情況下,系統(tǒng)滿足魯棒穩(wěn)定性的條件.聯(lián)合三式可知,設(shè)計(jì)任務(wù)就是尋找控制器K(s)使

    成立.

    混合靈敏度問(wèn)題便于通過(guò)系統(tǒng)頻域指標(biāo)進(jìn)行設(shè)計(jì),但沒(méi)有統(tǒng)一的解法.將圖3稍加變換可得如圖4所示的標(biāo)準(zhǔn)H∞控制結(jié)構(gòu).因此,可以通過(guò)求解標(biāo)準(zhǔn)H∞控制問(wèn)題來(lái)求解混合靈敏度問(wèn)題[6].

    圖4 標(biāo)準(zhǔn)H∞框圖

    2.2 加權(quán)函數(shù)的選取原則

    下面分別給出各加權(quán)函數(shù)的選取原則[7].

    (1)跟蹤性能加權(quán)函數(shù)We(s)的選取原則

    We(s)代表干擾的頻譜特性.在We(s)的選取中要考慮低頻干擾的頻帶寬度、按系統(tǒng)快速性要求確定的剪切頻率和噪聲抑制比例.因此選取的形式一般為

    式中,k1為比例系數(shù),T和τ為時(shí)間常數(shù).首先通過(guò)綜合性能指標(biāo)確定其階次,然后根據(jù)跟蹤帶寬確定參數(shù)T,根據(jù)噪聲抑制性確定參數(shù)τ,最后根據(jù)We(s)對(duì)高低頻增益的要求以及性能指標(biāo)γ值確定參數(shù)k1.

    We(s)一般采用低通的形式.低通部分是為了保證系統(tǒng)的低頻高增益,以滿足系統(tǒng)的跟蹤性能;高頻段降低增益是為了抑制高頻干擾信號(hào)的影響.為了提高平移系統(tǒng)性能,取We(s)的由低頻高增益到高頻低增益的斜率為-40dB/dec,即

    (2)輸出加權(quán)函數(shù)Wy(s)的選取原則

    因Wy(s)表示乘性攝動(dòng)的范數(shù)界,代表模型中存在的高頻未建模動(dòng)態(tài),因此應(yīng)具有高通特性.其斜率可以取的大些,以保證閉環(huán)系統(tǒng)對(duì)高頻干擾的抑制效果.

    (3)控制量加權(quán)函數(shù)Wu(s)的選擇原則

    因Wu(s)表示加性攝動(dòng)的范數(shù)界,代表模型中存在的中低頻參數(shù)不確定性,因此應(yīng)該具有低通特性,一般取為慣性環(huán)節(jié)或常數(shù).當(dāng)中低頻參數(shù)攝動(dòng)是非主要不確定性時(shí),通過(guò)其得到合適的控制量即可.本系統(tǒng)中Wu(s)的作用是保證控制電壓處于±10V之間,以便實(shí)現(xiàn)數(shù)字控制.

    2.3 基于LM I的H∞控制律

    為便于求解K(s),引入下面的引理.

    引理1[8].設(shè)E,F(xiàn),Q為已知矩陣,且QT=Q,則存在矩陣X,滿足

    的充分必要條件是

    式中,矩陣E⊥,F(xiàn)⊥分別是E和E⊥的垂直補(bǔ).

    設(shè)P0(s)的狀態(tài)空間實(shí)現(xiàn)為(AP0,BP0,CP0,DP0);We的狀態(tài)空間實(shí)現(xiàn)為(AWe,BWe,CWe,DWe);Wu的狀態(tài)空間實(shí)現(xiàn)為(AWu,BWu,CWu,DWu);Wy(s)的模型為Wy(s)=CWy(sI-AWy)-1BWy+Tmsm+…+T1s+T0.

    (3)教師對(duì)作業(yè)的設(shè)置缺乏系統(tǒng)性。作業(yè)的設(shè)計(jì)必須根據(jù)課程標(biāo)準(zhǔn)、教材和學(xué)生的實(shí)際學(xué)習(xí)需要,系統(tǒng)地考慮年級(jí)、單元、課時(shí)的連貫性和一致性,符合語(yǔ)言學(xué)習(xí)的規(guī)律,兼顧聽(tīng)、說(shuō)、讀、看、寫能力的全面發(fā)展。教師對(duì)作業(yè)的設(shè)計(jì)缺乏與生活實(shí)際的聯(lián)系和趣味性,導(dǎo)致學(xué)生產(chǎn)生排斥心理,持敷衍應(yīng)付的態(tài)度完成作業(yè)。

    由圖4可知,廣義輸入r、u,廣義輸出z、e和廣義被控對(duì)象G之間滿足

    其中

    引理2[3].設(shè)廣義對(duì)象G滿足(A,B2)可鎮(zhèn)定、(A,C2)可檢測(cè),則存在使圖4所示的系統(tǒng)內(nèi)部穩(wěn)定、且滿足‖Hzw‖∞<γ的控制器K(s)的充分必要條件為下述LMI:

    存在解矩陣X>0和Y>0.其中n和nk分別表示G和K的次數(shù).

    根據(jù)引理1和引理2及其證明過(guò)程,可得H∞控制律的存在條件和設(shè)計(jì)方法.

    3 橫向平移控制系統(tǒng)設(shè)計(jì)

    3.1 被控對(duì)象測(cè)定

    系統(tǒng)處于開(kāi)環(huán)狀態(tài),給定正弦電壓指令幅值為6 V,從0.1 Hz掃頻至40 Hz得到被控對(duì)象的頻率響應(yīng)曲線族).取的中心線為標(biāo)稱對(duì)象P0(s),通過(guò)調(diào)整式(1)中的三個(gè)參數(shù),采用最小二乘法對(duì)P0(s)進(jìn)行擬合,得P0(s)為

    根據(jù)

    式(22)證實(shí),高頻不確定性為機(jī)械諧振.考慮到魯棒控制器的階次與其有關(guān),并結(jié)合加權(quán)函數(shù)的意義和選取原則,取Wy(s)為

    3.2 性能加權(quán)函數(shù)確定

    橫向平移系統(tǒng)的主要性能指標(biāo)為:

    1)閉環(huán)系統(tǒng)動(dòng)態(tài)響應(yīng)滿足:幅度誤差小于10%、相位滯后小于10°的帶寬高于0.5 Hz;

    2)位置伺服誤差:±1mm;

    3)最低平穩(wěn)速度:1mm/s.

    下面將結(jié)合性能指標(biāo)進(jìn)行基于γ值的We(s)的參數(shù)確定.

    將乘性不確定性加權(quán)函數(shù)Wy(s)重新列寫,如式(24)所示,并選取控制量加權(quán)函數(shù)Wu(s)如式(25)所示

    根據(jù)指標(biāo)要求1),設(shè)定兩個(gè)轉(zhuǎn)折頻率分別為1 Hz和10 Hz,即

    由于此時(shí)We(s)的頻域特性轉(zhuǎn)折段較窄,導(dǎo)致了低頻高增益和高頻低增益只相差25 dB.考慮到高頻低增益是為了抑制高頻干擾,取

    則We(s)的高低增益之間相差60 dB,并且0.5Hz處于慣性環(huán)節(jié)的轉(zhuǎn)折頻率前.由于r=0.125 9,考慮到增大k可更有效抑制10Hz以上的高頻擾動(dòng),故取k=20時(shí),γ=0.706.繼續(xù)增大k,使We(s)低頻增益再高些,且將γ值逼到1附近.最后選定γ值小于1時(shí)的k的最大整數(shù)值33,相應(yīng)的γ=0.983.即最終確定的We(s)為

    綜上,混合靈敏度問(wèn)題的三個(gè)加權(quán)函數(shù)的幅頻特性如圖5所示.

    3.3 H∞控制律與系統(tǒng)魯棒穩(wěn)定性

    利用式(16),可得變換后廣義對(duì)象的系統(tǒng)矩陣G為

    根據(jù)廣義對(duì)象G,利用引理2

    圖5 加權(quán)函數(shù)的幅頻特性

    求得H∞控制律K(s)為

    式中,

    這是一個(gè)6階的線性控制器,相應(yīng)的H∞最優(yōu)性能指標(biāo)為γ=0.983.

    考慮到本系統(tǒng)采用W indows2000與RTX相結(jié)合的定時(shí)機(jī)制以保證計(jì)算機(jī)控制系統(tǒng)采用周期為0.5ms,因此在T=0.5ms的條件下應(yīng)用雙線性變換對(duì)H∞控制器進(jìn)行離散化,得到的離散控制器K(z)為

    式中,

    下面分析閉環(huán)系統(tǒng)的魯棒穩(wěn)定性.

    系統(tǒng)被控對(duì)象模型存在攝動(dòng)時(shí)的開(kāi)環(huán)傳遞函數(shù)的頻域特性如圖6所示.

    圖6 系統(tǒng)對(duì)象攝動(dòng)時(shí)的開(kāi)環(huán)頻率特性

    由圖6可知,被控對(duì)象模型存在攝動(dòng)時(shí)系統(tǒng)的穿越頻率為6 Hz,相角裕度為+45°,滿足穩(wěn)定性條件.

    不確定性攝動(dòng)界函數(shù)1/Wy(s)與補(bǔ)靈敏度函數(shù)T(s)的幅頻特性如圖7所示.

    圖7 系統(tǒng)魯棒性的攝動(dòng)界

    考慮到圖7中不確定性攝動(dòng)界函數(shù)1/Wy(s)表示系統(tǒng)魯棒穩(wěn)定性的界,而補(bǔ)靈敏度函數(shù)T(s)的幅頻特性一直處于1/Wy(s)的幅頻特性之下,即

    由小增益定理知,系統(tǒng)具有魯棒穩(wěn)定性.

    3.4 系統(tǒng)測(cè)試結(jié)果

    將離散化的H∞控制器K(z)應(yīng)用于平移系統(tǒng),分別輸入正弦指令、斜坡指令,得到系統(tǒng)實(shí)際運(yùn)行曲線如圖8~10所示.

    圖8 正弦響應(yīng)

    圖9 正弦指令時(shí)的控制電壓

    圖10 斜坡響應(yīng)

    4 結(jié) 論

    本文以九自由度空間交會(huì)對(duì)接模擬器橫向平移控制系統(tǒng)設(shè)計(jì)為背景,采用了非結(jié)構(gòu)不確定性描述系統(tǒng)中的高頻諧振,并給出了對(duì)象的不確定性模型,確定了魯棒H∞控制方案.研究了混合靈敏度問(wèn)題的加權(quán)函數(shù)選取原則,并結(jié)合系統(tǒng)性能指標(biāo)和主要不確定性給出了確定方法.推導(dǎo)了基于LMI的H∞控制律的存在條件,并給出了控制系統(tǒng)設(shè)計(jì)過(guò)程.系統(tǒng)測(cè)試結(jié)果滿足指標(biāo)要求,并經(jīng)受住了用戶長(zhǎng)時(shí)間測(cè)試和仿真實(shí)驗(yàn)應(yīng)用的考驗(yàn),說(shuō)明了本文方法可以解決實(shí)際系統(tǒng)的魯棒性問(wèn)題.

    [1] 楊嘉墀.航天器軌道動(dòng)力學(xué)與控制(下)[M].北京:宇航出版社,2007

    [2] Zhou K M,Doyle J C,Glover K.魯棒與最優(yōu)控制[M].毛劍琴,鐘宜生,林巖,譯.北京:國(guó)防工業(yè)出版社,2006

    [3] 梅生偉,申鐵龍,劉康志.現(xiàn)代魯棒控制理論與應(yīng)用[M].北京:清華大學(xué)出版社,2003

    [4] 趙希梅.永磁直線交流伺服系統(tǒng)的魯棒跟蹤控制研究[D].沈陽(yáng)工業(yè)大學(xué)碩士學(xué)位論文,2006

    [5] Zames G,F(xiàn)rancis B A.Feedback minimax sensitivity and optimal robustness[J].IEEE Trans Actions Automatic Control,1983,28(5):585-601

    [6] Muhammad EM,Naeem As,Waqars.Weight selection in H∞loop shaping using lead/lag compensators[C].IEEE 2ndInternational Conference on Emerging Technologies,Washington DC,USA,Nov 2006

    [7] 吳旭東,解學(xué)書.H∞魯棒控制中的加權(quán)陣的選擇[J].清華大學(xué)學(xué)報(bào)(自然科學(xué)版),1997,1(37):27-30

    [8] Iwasaki T,Skelton E.All controllers for the general H∞control problem:LM I existence conditions and state space formulas[J].Automatica,1994,30(8):1307-1317

    [9] Englewood C,Gahinet P.Explicit controller formulas for LM I-based H∞synthesis[J].Automatica,1996,32(7):1007-1014

    Design of Space Rendezvous and Docking Transverse Linear Movement Simulator Control System Based on H∞Theory

    LIGuiming1,2,YAO Yu3
    (1.Beijing Institute of Control Engineering,Beijing 100190,China;2.National Key Laboratory of Science and Technology on Space Intelligent Control,Beijing 100190,China;3.Control and Simulation Center,Harbin Institute of Technology,Harbin 150001,China)

    The paper investigates the design of simulator transverse linear movement control system for space rendezvous and docking,and discusses the unstructured uncertainty model.The robust H∞control strategy is designed according to performance indexes and main uncertainties,and the principle and method for choosing weighting functions of the mixed sensitivity optimal design are studied.The robust H∞controller is provided and obtained through solving three LMIs.The physical test resultsmeet the requestments,and show the effectiveness of the proposed method.

    space rendezvous and docking simulator;transverse linearmovement control;H∞control

    2009-02-01

    李貴明(1983—),男,黑龍江人,博士研究生,研究方向?yàn)楹教炱鲗?dǎo)航、制導(dǎo)與控制(e-mail:hitlgm@163.com).

    TP273

    A

    1674-1579(2009)06-0030-06

    猜你喜歡
    魯棒交會(huì)模擬器
    “織交會(huì)”打造大朗樣板
    了不起的安檢模擬器
    盲盒模擬器
    劃船模擬器
    基于學(xué)習(xí)的魯棒自適應(yīng)評(píng)判控制研究進(jìn)展
    2019年,水交會(huì)來(lái)了!
    目標(biāo)魯棒識(shí)別的抗旋轉(zhuǎn)HDO 局部特征描述
    立方體星交會(huì)對(duì)接和空間飛行演示
    太空探索(2016年9期)2016-07-12 09:59:53
    基于Cauchy魯棒函數(shù)的UKF改進(jìn)算法
    目標(biāo)軌跡更新的點(diǎn)到點(diǎn)魯棒迭代學(xué)習(xí)控制
    欧美日本视频| 日本与韩国留学比较| 久久国产乱子免费精品| 一级毛片久久久久久久久女| 婷婷色av中文字幕| a级一级毛片免费在线观看| 麻豆久久精品国产亚洲av| 国产乱人视频| 国产成人精品婷婷| 国产成人aa在线观看| 亚洲av在线观看美女高潮| 免费看光身美女| 永久网站在线| 777米奇影视久久| 日韩欧美精品免费久久| 国产大屁股一区二区在线视频| 亚洲色图av天堂| 午夜日本视频在线| 综合色丁香网| 亚洲精品中文字幕在线视频 | 午夜福利网站1000一区二区三区| 黄片wwwwww| 国产av不卡久久| 啦啦啦中文免费视频观看日本| 国产精品久久久久久精品电影小说 | 日韩国内少妇激情av| 国产成人精品久久久久久| 国产爱豆传媒在线观看| 国产av在哪里看| 日本与韩国留学比较| www.色视频.com| 亚洲欧洲日产国产| 校园人妻丝袜中文字幕| 大片免费播放器 马上看| 久99久视频精品免费| 婷婷色av中文字幕| 精品久久久噜噜| 欧美人与善性xxx| 欧美区成人在线视频| 国产精品一区二区性色av| 国产亚洲5aaaaa淫片| 亚洲av男天堂| 久久综合国产亚洲精品| 亚洲欧美清纯卡通| 麻豆精品久久久久久蜜桃| 免费黄频网站在线观看国产| 国产亚洲精品av在线| 亚洲欧美一区二区三区黑人 | 国产av国产精品国产| 日本免费a在线| 午夜亚洲福利在线播放| 午夜福利视频1000在线观看| 中文乱码字字幕精品一区二区三区 | 高清毛片免费看| av黄色大香蕉| 午夜免费男女啪啪视频观看| 欧美三级亚洲精品| 国产黄a三级三级三级人| 亚洲国产精品成人久久小说| 日日摸夜夜添夜夜添av毛片| 3wmmmm亚洲av在线观看| 国模一区二区三区四区视频| 日韩电影二区| 自拍偷自拍亚洲精品老妇| 激情 狠狠 欧美| 韩国av在线不卡| 97超碰精品成人国产| 国产极品天堂在线| 日韩成人av中文字幕在线观看| 国产69精品久久久久777片| 高清视频免费观看一区二区 | 亚洲天堂国产精品一区在线| 午夜福利网站1000一区二区三区| 偷拍熟女少妇极品色| 日本黄色片子视频| 白带黄色成豆腐渣| 亚洲精品色激情综合| 日本欧美国产在线视频| 免费看不卡的av| 国产av在哪里看| 搡老乐熟女国产| 亚洲欧美中文字幕日韩二区| 乱码一卡2卡4卡精品| 免费看不卡的av| 成人亚洲精品一区在线观看 | 国产精品精品国产色婷婷| 欧美成人一区二区免费高清观看| 最近的中文字幕免费完整| 久久精品国产亚洲av天美| 18+在线观看网站| 欧美区成人在线视频| 国产精品福利在线免费观看| 一级片'在线观看视频| 日日摸夜夜添夜夜爱| 简卡轻食公司| 久久久a久久爽久久v久久| 国产成人精品福利久久| 久久久午夜欧美精品| 国产亚洲午夜精品一区二区久久 | 美女内射精品一级片tv| 免费黄色在线免费观看| 极品少妇高潮喷水抽搐| 国产伦一二天堂av在线观看| 成人毛片60女人毛片免费| 精品久久久久久电影网| 97在线视频观看| 深爱激情五月婷婷| 久久国内精品自在自线图片| 国产v大片淫在线免费观看| 亚洲在线观看片| 亚洲性久久影院| 永久网站在线| eeuss影院久久| 久久久久久久午夜电影| 日韩欧美国产在线观看| 免费看不卡的av| 日本-黄色视频高清免费观看| av免费观看日本| 国产精品日韩av在线免费观看| 好男人视频免费观看在线| 欧美3d第一页| 一级毛片aaaaaa免费看小| 国产老妇女一区| 日韩中字成人| 久久99热这里只频精品6学生| 99热全是精品| 国产亚洲精品av在线| 大话2 男鬼变身卡| 国产伦精品一区二区三区视频9| 国产亚洲午夜精品一区二区久久 | av在线播放精品| 亚洲人与动物交配视频| 青春草视频在线免费观看| 成人亚洲精品一区在线观看 | 最近的中文字幕免费完整| 亚洲av成人精品一区久久| 国产精品一区二区在线观看99 | 丝袜美腿在线中文| 啦啦啦啦在线视频资源| 久久这里有精品视频免费| av.在线天堂| 国产伦精品一区二区三区视频9| 国产黄色小视频在线观看| 一级黄片播放器| 大又大粗又爽又黄少妇毛片口| 欧美区成人在线视频| 欧美bdsm另类| 国产精品精品国产色婷婷| 水蜜桃什么品种好| 精品久久久久久久人妻蜜臀av| 男女啪啪激烈高潮av片| 一边亲一边摸免费视频| 在现免费观看毛片| 国产精品麻豆人妻色哟哟久久 | 欧美成人a在线观看| 高清毛片免费看| 亚洲最大成人中文| 美女大奶头视频| 我要看日韩黄色一级片| 春色校园在线视频观看| h日本视频在线播放| 男女啪啪激烈高潮av片| 又粗又硬又长又爽又黄的视频| 男插女下体视频免费在线播放| 成年免费大片在线观看| 久久久久久久国产电影| 中文字幕av在线有码专区| 久久久亚洲精品成人影院| 黑人高潮一二区| 大片免费播放器 马上看| 亚洲精品国产av蜜桃| 中国美白少妇内射xxxbb| 亚洲激情五月婷婷啪啪| 国产男人的电影天堂91| 禁无遮挡网站| 老女人水多毛片| 99视频精品全部免费 在线| 久久久久精品久久久久真实原创| 日本免费a在线| 成人二区视频| 伊人久久精品亚洲午夜| 黄片无遮挡物在线观看| 偷拍熟女少妇极品色| 一级黄片播放器| 中文字幕av在线有码专区| 91午夜精品亚洲一区二区三区| 777米奇影视久久| av女优亚洲男人天堂| 啦啦啦韩国在线观看视频| 国内精品美女久久久久久| 亚洲熟妇中文字幕五十中出| 日韩亚洲欧美综合| 乱码一卡2卡4卡精品| 成人亚洲精品一区在线观看 | 91精品伊人久久大香线蕉| 国产v大片淫在线免费观看| 国产探花极品一区二区| 91av网一区二区| 两个人视频免费观看高清| 亚洲一区高清亚洲精品| 久久亚洲国产成人精品v| 婷婷色麻豆天堂久久| 亚洲欧洲日产国产| 水蜜桃什么品种好| 小蜜桃在线观看免费完整版高清| 男人舔奶头视频| 欧美一区二区亚洲| 色视频www国产| 国产精品人妻久久久久久| 黄片无遮挡物在线观看| 欧美日韩视频高清一区二区三区二| 国产视频内射| 18+在线观看网站| 天天躁夜夜躁狠狠久久av| 国产国拍精品亚洲av在线观看| 精品欧美国产一区二区三| 99热这里只有是精品在线观看| 少妇的逼水好多| 国产永久视频网站| 久久久久久久久久黄片| 一区二区三区免费毛片| 高清毛片免费看| 大香蕉久久网| 久久精品国产自在天天线| 日本av手机在线免费观看| 韩国高清视频一区二区三区| 久久久久久久久中文| 国产女主播在线喷水免费视频网站 | 777米奇影视久久| 国产淫片久久久久久久久| 日韩一区二区视频免费看| 久久久欧美国产精品| 国产精品一区二区性色av| 最近最新中文字幕免费大全7| 18禁在线无遮挡免费观看视频| 亚洲人成网站在线播| 大香蕉久久网| 日韩中字成人| 国精品久久久久久国模美| 久热久热在线精品观看| 九草在线视频观看| 神马国产精品三级电影在线观看| 国语对白做爰xxxⅹ性视频网站| 国产探花在线观看一区二区| 亚洲精品视频女| 深爱激情五月婷婷| 如何舔出高潮| 网址你懂的国产日韩在线| 一边亲一边摸免费视频| 免费黄频网站在线观看国产| 男女下面进入的视频免费午夜| 成人性生交大片免费视频hd| 国产白丝娇喘喷水9色精品| 一边亲一边摸免费视频| 在线观看人妻少妇| 日韩伦理黄色片| 水蜜桃什么品种好| 成年女人看的毛片在线观看| 久久综合国产亚洲精品| 欧美成人精品欧美一级黄| 人妻一区二区av| 精品人妻偷拍中文字幕| 亚洲18禁久久av| 能在线免费观看的黄片| 精品久久久久久久久久久久久| 久久久久久久久久成人| 亚洲自拍偷在线| 免费观看性生交大片5| 亚洲精品国产av蜜桃| 九九在线视频观看精品| 亚洲av福利一区| 亚洲va在线va天堂va国产| 天堂网av新在线| 日本-黄色视频高清免费观看| av又黄又爽大尺度在线免费看| 精品一区二区三区视频在线| 一级毛片aaaaaa免费看小| 黄片无遮挡物在线观看| 美女cb高潮喷水在线观看| 国产毛片a区久久久久| 丰满乱子伦码专区| 免费人成在线观看视频色| 日韩制服骚丝袜av| 午夜精品国产一区二区电影 | 久久精品国产自在天天线| 免费在线观看成人毛片| freevideosex欧美| 免费黄色在线免费观看| 亚洲激情五月婷婷啪啪| 熟女人妻精品中文字幕| ponron亚洲| 国语对白做爰xxxⅹ性视频网站| 欧美激情久久久久久爽电影| 成人漫画全彩无遮挡| 中文字幕久久专区| 国产亚洲91精品色在线| 午夜精品一区二区三区免费看| 久久99热这里只有精品18| 又爽又黄无遮挡网站| 日韩人妻高清精品专区| 18禁裸乳无遮挡免费网站照片| 免费看光身美女| 国产精品精品国产色婷婷| av专区在线播放| 久久久久久久久久黄片| 日本一本二区三区精品| 亚洲自偷自拍三级| 国产色爽女视频免费观看| 美女大奶头视频| 日本欧美国产在线视频| 国产成人一区二区在线| 少妇丰满av| 高清欧美精品videossex| kizo精华| 亚洲欧洲国产日韩| 国产精品熟女久久久久浪| 国产成人午夜福利电影在线观看| 少妇被粗大猛烈的视频| 午夜精品在线福利| 亚洲国产日韩欧美精品在线观看| 日韩制服骚丝袜av| 午夜爱爱视频在线播放| 久久鲁丝午夜福利片| 久久久久久久久久黄片| a级毛色黄片| 精品不卡国产一区二区三区| 国产视频首页在线观看| 成人亚洲精品av一区二区| 日韩欧美 国产精品| 日韩成人av中文字幕在线观看| 在线观看人妻少妇| 亚洲真实伦在线观看| 能在线免费看毛片的网站| 精品久久久久久久人妻蜜臀av| 日本午夜av视频| 亚洲欧美清纯卡通| 97超碰精品成人国产| 国内揄拍国产精品人妻在线| 五月天丁香电影| 国产成人精品福利久久| 色吧在线观看| 亚洲成人一二三区av| 可以在线观看毛片的网站| 欧美高清性xxxxhd video| 国产成人一区二区在线| 男人舔奶头视频| 六月丁香七月| 日韩在线高清观看一区二区三区| 国产午夜精品久久久久久一区二区三区| 国产精品久久久久久av不卡| 男女那种视频在线观看| 亚洲成人中文字幕在线播放| 亚洲国产精品成人久久小说| 久久久久久久久久久免费av| 淫秽高清视频在线观看| 国产色婷婷99| 国产精品av视频在线免费观看| 亚洲精品中文字幕在线视频 | 免费人成在线观看视频色| www.色视频.com| 亚洲成人一二三区av| 国产男女超爽视频在线观看| 一级毛片久久久久久久久女| 国产日韩欧美在线精品| 免费观看a级毛片全部| 美女cb高潮喷水在线观看| 不卡视频在线观看欧美| 久久久久久久午夜电影| 最近最新中文字幕大全电影3| 两个人的视频大全免费| 淫秽高清视频在线观看| 亚洲三级黄色毛片| 寂寞人妻少妇视频99o| 免费观看性生交大片5| 免费av观看视频| 99热网站在线观看| 色视频www国产| videossex国产| 中文在线观看免费www的网站| 日本黄大片高清| 一个人免费在线观看电影| 乱码一卡2卡4卡精品| 久久99蜜桃精品久久| 免费观看的影片在线观看| 伦精品一区二区三区| 国产精品人妻久久久影院| 久热久热在线精品观看| 国产高潮美女av| 欧美xxxx黑人xx丫x性爽| 麻豆成人av视频| av又黄又爽大尺度在线免费看| 午夜精品一区二区三区免费看| 九色成人免费人妻av| 日韩一区二区视频免费看| 久久久久久久午夜电影| 秋霞伦理黄片| 噜噜噜噜噜久久久久久91| kizo精华| 免费观看性生交大片5| 国产精品.久久久| 91精品伊人久久大香线蕉| 国产高清不卡午夜福利| 老司机影院成人| 中文天堂在线官网| 亚洲精品成人av观看孕妇| 日本与韩国留学比较| 91精品国产九色| 成人国产麻豆网| 直男gayav资源| 美女黄网站色视频| 99久久精品一区二区三区| 久久久久久久午夜电影| 久久久久久久久久成人| 亚洲精品影视一区二区三区av| 国产av码专区亚洲av| 国产亚洲精品久久久com| 麻豆久久精品国产亚洲av| 国产美女午夜福利| 久久99蜜桃精品久久| 18禁动态无遮挡网站| 亚洲国产av新网站| 欧美97在线视频| 亚洲国产欧美人成| 久久久久久久久久成人| 18禁在线无遮挡免费观看视频| 日本免费在线观看一区| 国产麻豆成人av免费视频| 日韩一本色道免费dvd| 色视频www国产| 精品久久久久久久久亚洲| 久久精品国产亚洲av天美| 日韩欧美三级三区| 久久这里有精品视频免费| av在线天堂中文字幕| 免费大片黄手机在线观看| 亚洲高清免费不卡视频| 1000部很黄的大片| 美女高潮的动态| 又大又黄又爽视频免费| 国产色爽女视频免费观看| 精品久久国产蜜桃| 内射极品少妇av片p| 欧美成人a在线观看| 免费观看av网站的网址| 秋霞在线观看毛片| 五月伊人婷婷丁香| av在线老鸭窝| 色综合站精品国产| 黄色配什么色好看| 亚洲四区av| 日韩精品青青久久久久久| 嘟嘟电影网在线观看| 精品一区二区免费观看| 国产欧美日韩精品一区二区| 一个人免费在线观看电影| 免费av不卡在线播放| 日本欧美国产在线视频| 国产精品国产三级专区第一集| 亚洲av不卡在线观看| 美女黄网站色视频| 欧美高清成人免费视频www| 一区二区三区四区激情视频| 免费观看av网站的网址| 久久精品国产亚洲网站| 亚洲美女视频黄频| 亚洲丝袜综合中文字幕| 国产黄片视频在线免费观看| 精品不卡国产一区二区三区| 国产黄色视频一区二区在线观看| 91av网一区二区| 中文字幕人妻熟人妻熟丝袜美| 亚洲内射少妇av| 亚洲性久久影院| 好男人视频免费观看在线| 永久免费av网站大全| 国产淫片久久久久久久久| 亚洲av二区三区四区| 亚洲最大成人手机在线| 国产黄频视频在线观看| 国产乱人偷精品视频| 日本猛色少妇xxxxx猛交久久| 高清视频免费观看一区二区 | 国产毛片a区久久久久| 免费观看的影片在线观看| 久久久午夜欧美精品| 亚洲人成网站在线播| 亚洲国产欧美在线一区| 特级一级黄色大片| 午夜福利视频1000在线观看| 人人妻人人澡欧美一区二区| 色尼玛亚洲综合影院| 国产成人福利小说| 岛国毛片在线播放| 亚洲av免费高清在线观看| 神马国产精品三级电影在线观看| 91精品国产九色| 欧美bdsm另类| 亚洲久久久久久中文字幕| 成人亚洲欧美一区二区av| www.色视频.com| 久久精品国产亚洲av天美| 2022亚洲国产成人精品| 99热全是精品| 国产黄色免费在线视频| 欧美xxxx性猛交bbbb| 2022亚洲国产成人精品| 非洲黑人性xxxx精品又粗又长| 国产av在哪里看| 99热全是精品| 午夜福利高清视频| 最近2019中文字幕mv第一页| 亚洲精品乱码久久久久久按摩| 国产亚洲av嫩草精品影院| 99视频精品全部免费 在线| 嫩草影院入口| 人人妻人人澡人人爽人人夜夜 | 欧美xxxx性猛交bbbb| 亚洲精品亚洲一区二区| 免费在线观看成人毛片| 久久99精品国语久久久| 国产69精品久久久久777片| 日韩欧美精品v在线| 丝袜喷水一区| 一个人免费在线观看电影| 免费观看av网站的网址| 又爽又黄无遮挡网站| 亚洲av免费在线观看| 尾随美女入室| 插阴视频在线观看视频| 国产精品福利在线免费观看| av在线老鸭窝| 日本爱情动作片www.在线观看| 亚洲欧美精品自产自拍| 女人被狂操c到高潮| 色吧在线观看| 熟妇人妻久久中文字幕3abv| 久久久色成人| 国产不卡一卡二| 99久久精品一区二区三区| 91精品伊人久久大香线蕉| 日韩一本色道免费dvd| 别揉我奶头 嗯啊视频| 精品一区在线观看国产| 最近中文字幕2019免费版| 三级国产精品片| 日日撸夜夜添| 亚洲,欧美,日韩| 大香蕉97超碰在线| av国产免费在线观看| 啦啦啦韩国在线观看视频| 亚洲国产精品成人综合色| 久久久久久九九精品二区国产| 熟妇人妻不卡中文字幕| 国产亚洲最大av| 老司机影院毛片| 99热这里只有是精品在线观看| 国产伦精品一区二区三区视频9| 男女那种视频在线观看| 夫妻午夜视频| 午夜视频国产福利| 久久久久久久午夜电影| 午夜免费男女啪啪视频观看| 69人妻影院| 国产高清不卡午夜福利| 亚洲国产色片| 中文欧美无线码| 精品99又大又爽又粗少妇毛片| 亚洲欧美成人综合另类久久久| 别揉我奶头 嗯啊视频| 国产欧美另类精品又又久久亚洲欧美| 国产色婷婷99| 欧美一区二区亚洲| 色综合色国产| 亚洲四区av| 小蜜桃在线观看免费完整版高清| 久久精品综合一区二区三区| 久久这里有精品视频免费| 床上黄色一级片| 97超碰精品成人国产| 一夜夜www| 国产男女超爽视频在线观看| 一本一本综合久久| 亚洲精品亚洲一区二区| 亚洲精品中文字幕在线视频 | 久久99热这里只有精品18| 久久精品熟女亚洲av麻豆精品 | 十八禁网站网址无遮挡 | 男的添女的下面高潮视频| 人妻系列 视频| 国产一区二区在线观看日韩| 又黄又爽又刺激的免费视频.| 国产精品爽爽va在线观看网站| 久久久久免费精品人妻一区二区| 国精品久久久久久国模美| 久久久久久久久久人人人人人人| 久久久色成人| av在线观看视频网站免费| 久久久久久久大尺度免费视频| 美女国产视频在线观看| 国产欧美日韩精品一区二区| 超碰av人人做人人爽久久| av天堂中文字幕网| 最近的中文字幕免费完整| 日韩欧美三级三区| 国产精品.久久久| 国产精品美女特级片免费视频播放器| 少妇丰满av| 国产伦理片在线播放av一区| 在线免费十八禁| 国产成人免费观看mmmm| 精品少妇黑人巨大在线播放| 最近的中文字幕免费完整| 国产成年人精品一区二区| 成人亚洲精品一区在线观看 | 91在线精品国自产拍蜜月|