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

    失效航天器姿態(tài)接管的SDRE微分博弈控制

    2020-03-13 08:32:18羅建軍謝劍鋒
    宇航學報 2020年2期
    關鍵詞:微衛(wèi)星微分航天器

    柴 源,羅建軍,韓 楠,謝劍鋒

    (1. 西北工業(yè)大學航天學院,西安 710072;2. 西北工業(yè)大學航天飛行動力學技術重點實驗室,西安 710072;3. 北京航天飛行控制中心,北京 100094)

    0 引 言

    近年來,航天器在軌服務受到世界各國的重視。服務航天器與失效航天器形成組合體后通過自身執(zhí)行機構可接管失效航天器姿軌控系統(tǒng),實現(xiàn)穩(wěn)定控制,失效航天器的接管控制作為一種急需的在軌服務受到廣泛關注[1-2]。由于目前在軌運行的多數(shù)航天器事先未進行可服務和可維修設計,無交會對接輔助設備,因此,采用傳統(tǒng)的交會對接方法進行接管控制存在困難。文獻[3]研究了空間機械臂的失效航天器姿態(tài)接管控制,但空間機械臂難以對失效航天器進行長周期姿態(tài)調整。文獻[4]研究了繩系空間機器人對非合作目標進行抓捕的問題,但是存在繩索被割斷的可能性。此外,針對不同任務定制相應的空間機械臂和繩系機器人,也需要消耗巨大的成本和時間。近年來,微小衛(wèi)星技術得益于微電子技術的進步而快速發(fā)展,具有發(fā)射成本低、研制周期短、可批量生產(chǎn)等優(yōu)勢。鳳凰計劃項目[5]和iboss系列任務[6]利用多顆微小衛(wèi)星來實現(xiàn)空間飛行器的部件再利用和新衛(wèi)星的在軌組裝,是微小衛(wèi)星應用的新方向。文獻[7]通過多顆納星協(xié)同,設計辨識和控制一體化方法,實現(xiàn)失效衛(wèi)星的低成本姿態(tài)接管控制。其中多顆微衛(wèi)星如何通過分布式協(xié)同來實現(xiàn)失效航天器的姿態(tài)接管和長期的姿軌協(xié)同控制,仍有很多理論和方法問題需要解決。

    多顆微衛(wèi)星對失效航天器姿態(tài)接管控制問題的解決方案包括集中式控制和分布式控制等[8]。文獻[9]通過任務分配方法實現(xiàn)了多機器人的集中式控制,但是集中式控制依賴中央控制個體,容錯性差,更適合控制中心可靠的場合。分布式控制通過個體之間信息的交互實現(xiàn)運動行為的協(xié)同,具有靈活性好、容錯性高和容易拓展等優(yōu)勢。傳統(tǒng)的分布式協(xié)同控制[10-12]基于相對運動模型來設計控制器,更適合編隊、集群等空間任務;文獻[13]利用分布式控制分配方法實現(xiàn)空間細胞機器人的姿態(tài)接管,但是頻繁不斷進行的控制分配增加了機器人計算負擔。

    微分博弈方法研究了多個參與者的最優(yōu)決策與控制問題,其中每個個體通過優(yōu)化各自性能指標函數(shù)可得到各自的控制策略[14]。該方法可以視為一種特殊的分布式控制,其通信拓撲結構為可切換的全連通拓撲,允許個體的加入與退出,具有容錯性和靈活性。文獻[15]將微分博弈方法用于多智能體的一致性控制問題中;文獻[16]研究了在軌衛(wèi)星的追逃博弈問題;文獻[17]將二人非零和微分博弈方法用于航天器對接任務中。

    本文將失效航天器姿態(tài)接管控制問題轉化為多顆微衛(wèi)星的微分博弈問題,從而通過微衛(wèi)星的自主決策實現(xiàn)分布式姿態(tài)接管控制。該問題的關鍵在于求解耦合的HJ(Hamilton-Jacobi)方程得到多顆微衛(wèi)星的納什均衡策略[18],但是由于HJ方程的非線性及耦合性,很難得到解析解。傳統(tǒng)的解法大多采用離線計算[19],近年來提出的基于策略迭代的求解算法,利用神經(jīng)網(wǎng)絡近似值函數(shù),但是權值的調整依賴于經(jīng)驗[20]。本文采用狀態(tài)相關黎卡提方程(State-dependent Riccati equation, SDRE)方法,引入狀態(tài)相關系數(shù)(State-dependent coefficient, SDC)矩陣[19-20]將非線性模型轉化為狀態(tài)相關線性模型,從而得到狀態(tài)相關線性二次型微分博弈[21],通過對耦合的狀態(tài)相關黎卡提方程組進行求解,得到微分博弈的納什均衡,從而得到微衛(wèi)星控制策略的閉環(huán)表達式。

    1 問題描述

    1.1 基本假設及定義

    為了實現(xiàn)對失效航天器的姿態(tài)接管控制,N顆微衛(wèi)星需要通過自身的執(zhí)行機構實現(xiàn)組合航天器的姿態(tài)穩(wěn)定,如圖1所示。本文采用基于SDRE的微分博弈方法進行姿態(tài)接管控制?;炯僭O如下:

    1)微衛(wèi)星和失效航天器均為剛體。

    2)微衛(wèi)星已貼附在失效航天器上,且在姿態(tài)穩(wěn)定過程中,各顆微衛(wèi)星的相對位置保持不變。

    3)失效航天器的姿控系統(tǒng)完全失效,所需控制力矩僅由微衛(wèi)星提供。

    4)微衛(wèi)星上安裝有三軸正交的反作用飛輪。

    5)組合航天器的轉動慣量已辨識且保持不變。

    圖1 微衛(wèi)星姿態(tài)接管控制示意圖

    所涉及的坐標系定義如下:

    1)微衛(wèi)星i(i∈N)本體坐標系Oixiyizi:Oi表示微衛(wèi)星i的質心,xi,yi,zi分別表示微衛(wèi)星i的三個慣量主軸,在三個慣量主軸上各安裝一個反作用飛輪。

    2)參考坐標系Oxyz:用來描述每顆微衛(wèi)星的方位。選取微衛(wèi)星1的本體坐標系為參考坐標系,則其他微衛(wèi)星相對于微衛(wèi)星1的方位可由各自的本體坐標系與參考坐標系之間的轉換矩陣來描述。

    1.2 姿態(tài)運動學

    假設微衛(wèi)星和失效航天器均為剛體,因此姿態(tài)接管控制過程,由N顆微衛(wèi)星和失效航天器形成的組合航天器的姿態(tài)運動可由剛體航天器的姿態(tài)運動學方程來描述。本文采用修正羅德里格斯參數(shù)(Modified Rodrigues parameters, MRPs)描述姿態(tài),以避免奇異。組合航天器的姿態(tài)運動學方程為[7]

    (1)

    式中:σ∈R3×1是以MRPs表示的組合航天器姿態(tài)角;ω∈R3×1為組合航天器的姿態(tài)角速度。M(σ)為可逆矩陣,有

    (2)

    式中:σ×=[0,-σ3,σ2;σ3,0,-σ1;-σ2,σ1,0];I3為3×3的單位對角陣。

    1.3 姿態(tài)動力學

    組合航天器的姿態(tài)動力學方程為

    (3)

    式中:J∈R3×3是組合航天器在參考坐標系Oxyz中的轉動慣量矩陣;N是微衛(wèi)星的個數(shù);Cj∈R3×3是從微衛(wèi)星j的本體坐標系Ojxjyjzj到參考坐標系Oxyz的轉換矩陣,用來表示微衛(wèi)星j的方位[7];uj∈R3×1為微衛(wèi)星j在本體坐標系Ojxjyjzj的三軸姿態(tài)控制力矩;ω×=[0,-ω3,ω2;ω3,0,-ω1;-ω2,ω1,0]。

    1.4 微分博弈模型

    根據(jù)組合航天器的姿態(tài)運動學模型和姿態(tài)動力學模型,可得多顆微衛(wèi)星接管失效航天器姿態(tài)運動的微分博弈模型為

    (4)

    式中:x=[σT,ωT]T∈R6×1,且

    (5a)

    (5b)

    式中:03為3×3的全0矩陣。

    微衛(wèi)星的性能指標函數(shù)定義為

    (6)

    定義1.可行控制策略集[17]:若反饋控制策略ui(x)在緊集Ω上連續(xù),ui(0)=0,反饋控制策略集u(x)={u1(x),…,uN(x)}在Ω上能夠穩(wěn)定式(4),且式(6)對于任意的x0∈Ω是有限的,則u(x)在Ω上相對于式(6)定義為可行的,記作ui(x)∈Ψ(Ω)。

    對于給定的可行控制策略,微衛(wèi)星的值函數(shù)可表示為

    (7)

    并滿足Vi(x)>0,Vi(0)=0且連續(xù)可微。

    (8)

    多顆微衛(wèi)星的微分博弈可以由下式來描述

    (9)

    2 基于SDRE的微分博弈控制器設計

    2.1 狀態(tài)相關線性微分博弈模型

    微衛(wèi)星的值函數(shù)的微分等價形式為

    (10)

    定義哈密爾頓函數(shù)為

    (11)

    (12)

    將式(12)代入式(10)中可以得到N個耦合的HJ方程

    (13)

    (14)

    式中:A(x)=[03,M(σ);03,-J-1ω×J]∈R6×6是與狀態(tài)相關的狀態(tài)矩陣,且(A(x),Bj)滿足能控性秩判據(jù)。根據(jù)上述分析,將非線性微分博弈問題轉化為由式(14)和式(6)描述的狀態(tài)相關線性二次型微分博弈問題,能夠在保留了原模型的非線性特性的同時,使得微分博弈問題更加易于求解[19]。

    2.2 姿態(tài)控制器

    針對狀態(tài)相關線性模型(14),獨立優(yōu)化每顆微衛(wèi)星的性能指標函數(shù)(6),可以得到每顆微衛(wèi)星的控制策略,實現(xiàn)失效航天器的姿態(tài)接管控制。

    定理1. 對于狀態(tài)相關線性二次型微分博弈問題,在系統(tǒng)(A(x),Bi)完全能控的條件下,線性狀態(tài)反饋控制策略為

    (15)

    對稱正定矩陣Pi(x)∈R6×6為下面耦合狀態(tài)相關黎卡提方程組的解

    (16)

    證. 為了方便推導,在證明中將Pi(x)記為Pi,將A(x)記為A。

    假設最優(yōu)值函數(shù)(9)在狀態(tài)x(t)下有線性二次型形式的解

    (17)

    則可得

    (18)

    忽略式(18)中的高階項,則微衛(wèi)星i(?j≠i)對應的反饋控制為

    (19)

    將其代入式(10)中可得

    (20)

    整理為二次型形式得

    (21)

    (22)

    將其代入式(21)得

    (23)

    整理得

    0=xT(Qi+PiA+ATPi-

    (24)

    由于存在下式關系

    (25)

    則可以得到

    (26)

    (27)

    取i=1,2,…,N,可以得到耦合的狀態(tài)相關黎卡提方程組(16),對其求解可以得到矩陣Pi,從而得到狀態(tài)反饋控制策略(15)。

    根據(jù)上述分析,系統(tǒng)模型可以表示為如下形式

    (28)

    (29)

    式中:ψ(x)是高階展開項,且滿足

    (30)

    因此,在原點的鄰域內,常系數(shù)矩陣Acl(0)的穩(wěn)定,保證了系統(tǒng)的局部漸近穩(wěn)定性。

    注1.要求系統(tǒng)完全能控是為了保證系統(tǒng)的穩(wěn)定性。由于控制區(qū)間為無限時域,若出現(xiàn)狀態(tài)不能控,無論采取任何控制策略,性能指標都會趨于無窮大。

    注2.上述定理得到的狀態(tài)反饋控制律,在每一時刻更新狀態(tài)x(t),得到對應的狀態(tài)矩陣A(x),并通過求解耦合狀態(tài)相關黎卡提方程組(16)得到對稱正定矩陣Pi(x),在下一時刻重復上述過程,得到時變的狀態(tài)反饋調節(jié)器,能方便地實現(xiàn)閉環(huán)控制,這在工程應用中具有十分重要的意義。

    3 李雅普諾夫迭代

    對于耦合的狀態(tài)相關黎卡提方程組(16),要求在每一步求解耦合的代數(shù)黎卡提方程組。目前對于耦合代數(shù)黎卡提方程組的求解有很多研究成果[14,21],本文采用李雅普諾夫迭代法進行計算。該方法將耦合的代數(shù)黎卡提方程組降階并解耦為李雅普諾夫方程來獨立運算,算法速度快,準確性高。本文將文獻[21]中針對二人微分博弈得到的迭代算法進行歸納,給出一般性的李亞普諾夫迭代公式。

    迭代算法:

    k=0,1,2,…;i= 1,…,N

    (31)

    初值選擇:

    (32)

    通過迭代求解李亞普諾夫方程(31)可以得到狀態(tài)x(t)下的矩陣Pi(x)(i∈N)。

    4 仿真校驗

    為了驗證對失效航天器姿態(tài)接管的SDRE微分博弈控制方法的有效性和容錯性,對組合航天器的姿態(tài)穩(wěn)定控制過程進行仿真研究。

    4.1 有效性仿真

    假設采用三顆微衛(wèi)星對失效航天器進行姿態(tài)接管控制,微衛(wèi)星可提供的控制力矩滿足|u|≤0.25 N·m,其各自的本體坐標系到參考坐標系的轉換矩陣為

    組合航天器的轉動慣量矩陣為

    組合航天器以MRPs表示的初始姿態(tài)角為σ0=[0.015,0.009,0.013]T,初始姿態(tài)角速度為ω0=[0,0,0]T。取兩個微衛(wèi)星的權重矩陣分別為Q1=5I6,Q2=5I6,Q3=5I6,R11=R12=R13=0.01I3,R21=R22=R23=0.01I3,R31=R32=R33=0.01I3。取仿真步長為dt=0.01s。在上述參數(shù)下,仿真驗證SDRE微分博弈控制的有效性。

    圖2和圖3是組合航天器分別以MRPs和歐拉角表示的姿態(tài)角變化曲線,圖4是組合航天器的姿態(tài)角速度變化曲線。從圖中可以看出,組合航天器在40 s到達穩(wěn)定狀態(tài),精度為1×10-5量級,在仿真結束時,可以達到更高的精度。因此,基于微分博弈的姿態(tài)控制器可以實現(xiàn)失效航天器的姿態(tài)穩(wěn)定。

    圖2 姿態(tài)角σ曲線

    圖3 姿態(tài)角歐拉角曲線

    圖4 姿態(tài)角速度ω曲線

    圖5是姿態(tài)接管控制階段三顆微衛(wèi)星的控制力矩隨時間變化曲線。在初始控制階段,由于狀態(tài)量幅值較大,因此所需的控制力矩也較大,隨著狀態(tài)量幅值的不斷減小,所需控制力矩也逐漸減小。

    圖5 微衛(wèi)星控制力矩曲線

    4.2 容錯性仿真

    為了驗證該方法的容錯性,設計如下仿真場景。采用三顆微衛(wèi)星對失效航天器進行姿態(tài)接管控制時,微衛(wèi)星2在第6 s執(zhí)行機構失效,無法提供控制力矩,則前6 s姿態(tài)接管控制所需控制力矩由三顆微衛(wèi)星提供,在6 s之后僅由微衛(wèi)星1和微衛(wèi)星2提供。其他仿真參數(shù)與4.1節(jié)仿真參數(shù)相同。

    圖6為以歐拉角表示的姿態(tài)角和角速度隨時間變化曲線,圖7為三顆微衛(wèi)星的控制力矩變化曲線。從圖6可以看出,在存在一顆微衛(wèi)星失效的情況,其他的微衛(wèi)星將失效微衛(wèi)星視為不參與,更新參與者個數(shù),解算控制力矩使組合航天器的狀態(tài)趨于穩(wěn)定。從圖7可以看出,前6 s的仿真結果和4.1節(jié)中前6 s的仿真結果一致。在第6 s,微衛(wèi)星3無法提供控制力矩,微衛(wèi)星1和微衛(wèi)星2在更新參與者個數(shù)后計算各自的控制力矩,因此在圖7中可以看到6 s時微衛(wèi)星1和微衛(wèi)星2的控制力矩會有突變。6 s之后的微衛(wèi)星3的控制力矩為0,微衛(wèi)星1和微衛(wèi)星2的控制力矩逐漸趨于0。仿真結果表明本文提出的微分博弈方法具有容錯性。

    圖6 部分微衛(wèi)星失效的姿態(tài)角和角速度變化曲線

    圖7 部分微衛(wèi)星失效的微衛(wèi)星控制力矩曲線

    5 結 論

    本文采用SDRE微分博弈方法研究了多顆微衛(wèi)星接管姿態(tài)控制系統(tǒng)失效的航天器的控制問題,通過引入SDC矩陣將多顆微衛(wèi)星的非線性微分博弈轉化為狀態(tài)相關線性二次型微分博弈。多顆微衛(wèi)星通過獨立優(yōu)化各自的性能指標函數(shù)得到各自的控制策略,實現(xiàn)了多顆微衛(wèi)星接管姿態(tài)控制失效的航天器的分布式控制。該方法在保留非線性特性的同時,通過推導得到耦合的狀態(tài)相關黎卡提方程組來逼近微分博弈的納什均衡并得到控制策略的閉環(huán)表達式,采用李亞普諾夫迭代法進行求解,簡化了納什均衡的求解過程。數(shù)值算例的結果表明,SDRE微分博弈方法可以實現(xiàn)多顆微衛(wèi)星對失效航天器的分布式姿態(tài)接管控制,具有較高的控制精度和容錯性。

    猜你喜歡
    微衛(wèi)星微分航天器
    梅花鹿基因組微衛(wèi)星分布特征研究
    2022 年第二季度航天器發(fā)射統(tǒng)計
    國際太空(2022年7期)2022-08-16 09:52:50
    擬微分算子在Hp(ω)上的有界性
    上下解反向的脈沖微分包含解的存在性
    2019 年第二季度航天器發(fā)射統(tǒng)計
    國際太空(2019年9期)2019-10-23 01:55:34
    2018 年第三季度航天器發(fā)射統(tǒng)計
    國際太空(2018年12期)2019-01-28 12:53:20
    2018年第二季度航天器發(fā)射統(tǒng)計
    國際太空(2018年9期)2018-10-18 08:51:32
    林麝全基因組微衛(wèi)星分布規(guī)律研究
    四川動物(2017年4期)2017-07-31 23:54:19
    借助微分探求連續(xù)函數(shù)的極值點
    對不定積分湊微分解法的再認識
    一区二区三区激情视频| 丰满人妻熟妇乱又伦精品不卡| 又大又爽又粗| 色综合站精品国产| 国产蜜桃级精品一区二区三区| 悠悠久久av| 国产亚洲精品一区二区www| 在线播放国产精品三级| 精品欧美国产一区二区三| 黄色视频,在线免费观看| 热99在线观看视频| 久久久久久人人人人人| 欧美激情久久久久久爽电影| 人妻丰满熟妇av一区二区三区| 成熟少妇高潮喷水视频| 亚洲国产色片| 欧美在线一区亚洲| 亚洲精品国产精品久久久不卡| 午夜福利欧美成人| 午夜福利高清视频| 黄色日韩在线| 中文字幕高清在线视频| 99久久成人亚洲精品观看| 99久久国产精品久久久| 久久精品91无色码中文字幕| a级毛片a级免费在线| 在线永久观看黄色视频| 亚洲精品一区av在线观看| 日本黄色视频三级网站网址| 亚洲一区二区三区不卡视频| 亚洲人成网站高清观看| 日本黄色片子视频| 后天国语完整版免费观看| 欧美在线一区亚洲| 最近视频中文字幕2019在线8| 18禁国产床啪视频网站| 日韩欧美精品v在线| 亚洲中文字幕一区二区三区有码在线看 | 搡老妇女老女人老熟妇| 免费高清视频大片| 男女做爰动态图高潮gif福利片| 神马国产精品三级电影在线观看| 午夜两性在线视频| 少妇熟女aⅴ在线视频| 国内精品久久久久精免费| 精品久久久久久成人av| 在线免费观看不下载黄p国产 | 亚洲色图 男人天堂 中文字幕| 久久久久久久精品吃奶| 国产一区二区三区在线臀色熟女| 在线免费观看的www视频| 91av网一区二区| 日本免费一区二区三区高清不卡| 久久国产精品影院| 99国产精品一区二区三区| 两个人看的免费小视频| 99久久综合精品五月天人人| 久99久视频精品免费| 中国美女看黄片| 国产高清有码在线观看视频| 床上黄色一级片| 黄片小视频在线播放| 精华霜和精华液先用哪个| 伊人久久大香线蕉亚洲五| 久久国产精品影院| 欧美激情久久久久久爽电影| 9191精品国产免费久久| 男女那种视频在线观看| 欧美国产日韩亚洲一区| 免费看日本二区| 国产成人一区二区三区免费视频网站| 给我免费播放毛片高清在线观看| 母亲3免费完整高清在线观看| 我要搜黄色片| 淫秽高清视频在线观看| 亚洲国产欧洲综合997久久,| 久久精品国产亚洲av香蕉五月| 免费看a级黄色片| 国产亚洲精品av在线| 成熟少妇高潮喷水视频| 免费大片18禁| 综合色av麻豆| 久久香蕉国产精品| 十八禁网站免费在线| 久久中文看片网| 精品国产美女av久久久久小说| 97人妻精品一区二区三区麻豆| 欧美又色又爽又黄视频| 精品电影一区二区在线| 日韩欧美 国产精品| 动漫黄色视频在线观看| 女同久久另类99精品国产91| 老熟妇乱子伦视频在线观看| 国内精品一区二区在线观看| 无限看片的www在线观看| 亚洲午夜理论影院| 一区二区三区国产精品乱码| 脱女人内裤的视频| 久久久久久久久免费视频了| 国产精品 欧美亚洲| www.精华液| 91麻豆av在线| 最新在线观看一区二区三区| e午夜精品久久久久久久| 99久久国产精品久久久| 可以在线观看的亚洲视频| 性色avwww在线观看| 亚洲五月天丁香| 欧美黑人巨大hd| a级毛片在线看网站| 母亲3免费完整高清在线观看| 91久久精品国产一区二区成人 | 美女免费视频网站| 国产男靠女视频免费网站| 精品电影一区二区在线| 日本免费一区二区三区高清不卡| 老汉色∧v一级毛片| av在线蜜桃| 国产午夜福利久久久久久| 九九在线视频观看精品| 在线a可以看的网站| 一个人看的www免费观看视频| 老司机深夜福利视频在线观看| 天天一区二区日本电影三级| 久久久久国产精品人妻aⅴ院| 国产毛片a区久久久久| 欧美日韩瑟瑟在线播放| 老熟妇乱子伦视频在线观看| 老熟妇乱子伦视频在线观看| 国语自产精品视频在线第100页| 老司机在亚洲福利影院| h日本视频在线播放| 我的老师免费观看完整版| 国产精品,欧美在线| 免费看十八禁软件| 18禁美女被吸乳视频| 亚洲性夜色夜夜综合| 亚洲九九香蕉| 国产精品自产拍在线观看55亚洲| 麻豆国产97在线/欧美| 欧美黑人欧美精品刺激| 99国产综合亚洲精品| 国产高清视频在线观看网站| avwww免费| 少妇人妻一区二区三区视频| 黄色 视频免费看| 成年版毛片免费区| 在线十欧美十亚洲十日本专区| 亚洲中文字幕一区二区三区有码在线看 | 两个人的视频大全免费| 国产探花在线观看一区二区| 午夜激情欧美在线| 伦理电影免费视频| 国产 一区 欧美 日韩| 国产高清三级在线| 99热这里只有是精品50| 特级一级黄色大片| 日韩 欧美 亚洲 中文字幕| 国产不卡一卡二| 一a级毛片在线观看| 搡老岳熟女国产| 亚洲乱码一区二区免费版| 国产精品影院久久| 婷婷精品国产亚洲av在线| 神马国产精品三级电影在线观看| 大型黄色视频在线免费观看| 国产爱豆传媒在线观看| 日本免费一区二区三区高清不卡| 我的老师免费观看完整版| 我的老师免费观看完整版| 亚洲成人久久性| 日韩欧美 国产精品| 午夜福利高清视频| 黄色女人牲交| 美女黄网站色视频| 91麻豆精品激情在线观看国产| 丰满的人妻完整版| 欧美午夜高清在线| 久久精品国产清高在天天线| avwww免费| 久久久成人免费电影| 1024香蕉在线观看| 美女黄网站色视频| 波多野结衣高清作品| 亚洲精品美女久久av网站| 日本 欧美在线| 少妇的丰满在线观看| 成人鲁丝片一二三区免费| 久久久久国产一级毛片高清牌| 国产精品亚洲av一区麻豆| 又爽又黄无遮挡网站| www国产在线视频色| 精品午夜福利视频在线观看一区| 亚洲avbb在线观看| 亚洲天堂国产精品一区在线| 高清在线国产一区| 国产av一区在线观看免费| 国产aⅴ精品一区二区三区波| 久久国产精品人妻蜜桃| 天堂网av新在线| 真人做人爱边吃奶动态| 久久久久免费精品人妻一区二区| 亚洲自偷自拍图片 自拍| 超碰成人久久| 日韩精品中文字幕看吧| 免费无遮挡裸体视频| 九九在线视频观看精品| 国产免费av片在线观看野外av| av天堂中文字幕网| 亚洲片人在线观看| 午夜福利免费观看在线| 久久午夜亚洲精品久久| 国产伦精品一区二区三区四那| av天堂在线播放| 最近视频中文字幕2019在线8| 三级毛片av免费| 国产精品野战在线观看| 亚洲专区中文字幕在线| 男女午夜视频在线观看| av国产免费在线观看| 亚洲美女视频黄频| 亚洲国产欧洲综合997久久,| 亚洲av熟女| 欧美xxxx黑人xx丫x性爽| 一本一本综合久久| 国内精品久久久久久久电影| 成熟少妇高潮喷水视频| 国产 一区 欧美 日韩| 天堂av国产一区二区熟女人妻| 久久久久免费精品人妻一区二区| 又黄又爽又免费观看的视频| av福利片在线观看| 国内精品美女久久久久久| 大型黄色视频在线免费观看| 国产精品久久久av美女十八| 一进一出好大好爽视频| 99在线人妻在线中文字幕| 精品国产乱码久久久久久男人| 欧美黑人巨大hd| 久久久国产成人免费| 成年免费大片在线观看| 99热这里只有精品一区 | 亚洲最大成人中文| 亚洲成人免费电影在线观看| www.自偷自拍.com| 一区福利在线观看| 亚洲精品美女久久av网站| 免费大片18禁| 18禁裸乳无遮挡免费网站照片| 亚洲国产色片| 国产成人欧美在线观看| 伦理电影免费视频| 久久久久久久久中文| 国产毛片a区久久久久| 三级国产精品欧美在线观看 | 观看美女的网站| 性色avwww在线观看| 99热6这里只有精品| 蜜桃久久精品国产亚洲av| 国产探花在线观看一区二区| 欧美绝顶高潮抽搐喷水| 夜夜看夜夜爽夜夜摸| 亚洲真实伦在线观看| 免费av毛片视频| 国产一区二区三区在线臀色熟女| 日本 欧美在线| 非洲黑人性xxxx精品又粗又长| 国产精品久久电影中文字幕| 日日夜夜操网爽| 日韩精品中文字幕看吧| 我的老师免费观看完整版| 美女免费视频网站| 免费高清视频大片| 亚洲七黄色美女视频| 真人一进一出gif抽搐免费| 久久久久性生活片| 18美女黄网站色大片免费观看| 国产麻豆成人av免费视频| 亚洲av电影不卡..在线观看| 啦啦啦韩国在线观看视频| 在线十欧美十亚洲十日本专区| 欧美日韩黄片免| www日本黄色视频网| 久久精品影院6| 亚洲国产精品成人综合色| 操出白浆在线播放| 色尼玛亚洲综合影院| 欧美高清成人免费视频www| 熟女电影av网| 97人妻精品一区二区三区麻豆| 久久精品影院6| 国产精品 国内视频| 欧美色欧美亚洲另类二区| 国产高清三级在线| 亚洲最大成人中文| 大型黄色视频在线免费观看| 日韩欧美一区二区三区在线观看| 国产精品精品国产色婷婷| 欧美成狂野欧美在线观看| 亚洲欧美精品综合久久99| 天堂av国产一区二区熟女人妻| 三级男女做爰猛烈吃奶摸视频| 日韩欧美三级三区| 观看免费一级毛片| 久久国产精品影院| 午夜免费成人在线视频| 宅男免费午夜| 成人午夜高清在线视频| 99热6这里只有精品| 久久国产精品人妻蜜桃| 日韩精品中文字幕看吧| 老汉色∧v一级毛片| avwww免费| 嫩草影院精品99| 后天国语完整版免费观看| 亚洲国产中文字幕在线视频| 一级毛片高清免费大全| 亚洲精品在线美女| 97人妻精品一区二区三区麻豆| 成年女人永久免费观看视频| 俺也久久电影网| 亚洲avbb在线观看| 亚洲一区二区三区不卡视频| 日本黄大片高清| 日韩欧美一区二区三区在线观看| 国产伦人伦偷精品视频| 国产美女午夜福利| 精品久久久久久久久久免费视频| 免费无遮挡裸体视频| 精品久久久久久,| 欧美性猛交黑人性爽| tocl精华| 最新美女视频免费是黄的| 亚洲午夜精品一区,二区,三区| 日韩欧美国产一区二区入口| 亚洲精品美女久久久久99蜜臀| 亚洲精华国产精华精| 欧美一级毛片孕妇| 搡老熟女国产l中国老女人| 天堂动漫精品| 精品久久蜜臀av无| 欧美乱色亚洲激情| 97碰自拍视频| 国产一区二区在线av高清观看| 国产精华一区二区三区| 国产成人啪精品午夜网站| 少妇熟女aⅴ在线视频| 欧美日本亚洲视频在线播放| 亚洲av日韩精品久久久久久密| 青草久久国产| 成人亚洲精品av一区二区| 一本综合久久免费| 国产av在哪里看| 亚洲中文av在线| 精品欧美国产一区二区三| 九色国产91popny在线| 又紧又爽又黄一区二区| 久久婷婷人人爽人人干人人爱| 国产高清有码在线观看视频| 婷婷丁香在线五月| 国产精品99久久久久久久久| 日韩精品青青久久久久久| 成年女人永久免费观看视频| 99在线视频只有这里精品首页| 日本一本二区三区精品| 久久精品91蜜桃| 97人妻精品一区二区三区麻豆| 99国产综合亚洲精品| 久久久国产成人精品二区| 看黄色毛片网站| 亚洲 国产 在线| 久久久久久久午夜电影| 在线免费观看不下载黄p国产 | 欧美三级亚洲精品| 精品久久久久久久毛片微露脸| 亚洲av成人精品一区久久| 一进一出好大好爽视频| 欧美一级毛片孕妇| 国产精品久久久久久亚洲av鲁大| tocl精华| av天堂中文字幕网| 亚洲av第一区精品v没综合| 亚洲中文日韩欧美视频| 长腿黑丝高跟| 久久午夜综合久久蜜桃| 999精品在线视频| 亚洲中文av在线| 亚洲欧美日韩高清在线视频| 久久久成人免费电影| 中出人妻视频一区二区| 人人妻,人人澡人人爽秒播| 在线观看日韩欧美| 搡老熟女国产l中国老女人| 美女免费视频网站| 久久精品91无色码中文字幕| 最近最新中文字幕大全电影3| 我的老师免费观看完整版| 免费观看人在逋| 级片在线观看| 午夜福利欧美成人| 欧美中文日本在线观看视频| 亚洲一区二区三区色噜噜| 黄色视频,在线免费观看| 久久午夜综合久久蜜桃| 久久精品夜夜夜夜夜久久蜜豆| 亚洲第一欧美日韩一区二区三区| 此物有八面人人有两片| 成年人黄色毛片网站| 国语自产精品视频在线第100页| 老汉色av国产亚洲站长工具| 久久久久精品国产欧美久久久| 日本 av在线| 欧美日韩亚洲国产一区二区在线观看| av中文乱码字幕在线| 亚洲人成电影免费在线| 一区二区三区激情视频| 国产精品 国内视频| 欧美精品啪啪一区二区三区| 琪琪午夜伦伦电影理论片6080| 国产在线精品亚洲第一网站| 在线永久观看黄色视频| 欧美丝袜亚洲另类 | 亚洲精品美女久久久久99蜜臀| 欧美中文综合在线视频| 在线视频色国产色| 午夜精品久久久久久毛片777| 色老头精品视频在线观看| 一个人免费在线观看电影 | 精品国产乱码久久久久久男人| 欧美成狂野欧美在线观看| 老司机在亚洲福利影院| 色播亚洲综合网| 91九色精品人成在线观看| 色综合站精品国产| 精品国内亚洲2022精品成人| 在线观看一区二区三区| 精品久久蜜臀av无| 18禁观看日本| 夜夜爽天天搞| 曰老女人黄片| 91在线精品国自产拍蜜月 | 国产伦一二天堂av在线观看| 亚洲精品美女久久av网站| 久久久久精品国产欧美久久久| 国产不卡一卡二| 91麻豆av在线| www日本黄色视频网| 最新在线观看一区二区三区| 身体一侧抽搐| 亚洲欧美一区二区三区黑人| 国内久久婷婷六月综合欲色啪| 亚洲电影在线观看av| 听说在线观看完整版免费高清| 国产成人啪精品午夜网站| 午夜视频精品福利| 性色av乱码一区二区三区2| 中文字幕高清在线视频| 波多野结衣高清作品| 欧美乱色亚洲激情| 一本综合久久免费| 国产欧美日韩一区二区三| 国产精品美女特级片免费视频播放器 | 久久久久免费精品人妻一区二区| 男女下面进入的视频免费午夜| 欧美乱妇无乱码| 欧美日韩国产亚洲二区| 宅男免费午夜| 最新在线观看一区二区三区| 久久香蕉精品热| 亚洲九九香蕉| 一本综合久久免费| 岛国视频午夜一区免费看| 国产精品亚洲av一区麻豆| 久久香蕉精品热| 久久中文字幕人妻熟女| 一本精品99久久精品77| 亚洲精品一区av在线观看| a在线观看视频网站| 好看av亚洲va欧美ⅴa在| 天天添夜夜摸| 真实男女啪啪啪动态图| 村上凉子中文字幕在线| 国产亚洲精品久久久com| 国产成+人综合+亚洲专区| 亚洲av免费在线观看| 欧美3d第一页| 亚洲精品在线美女| 精品免费久久久久久久清纯| 欧美乱色亚洲激情| 久久精品影院6| 热99re8久久精品国产| 色老头精品视频在线观看| 亚洲精品美女久久av网站| 亚洲avbb在线观看| 日韩精品青青久久久久久| www.www免费av| 美女cb高潮喷水在线观看 | av在线天堂中文字幕| 久久人妻av系列| 国产成人系列免费观看| 国产真实乱freesex| 欧美日韩精品网址| 久久久精品欧美日韩精品| 一本精品99久久精品77| 国产真人三级小视频在线观看| 欧美色欧美亚洲另类二区| 动漫黄色视频在线观看| 最近最新中文字幕大全免费视频| 窝窝影院91人妻| 亚洲激情在线av| 婷婷丁香在线五月| 国产精品久久久久久亚洲av鲁大| 黄片大片在线免费观看| 又黄又爽又免费观看的视频| 亚洲黑人精品在线| 国产精品99久久久久久久久| 19禁男女啪啪无遮挡网站| 在线播放国产精品三级| 亚洲乱码一区二区免费版| 人妻夜夜爽99麻豆av| 麻豆av在线久日| 亚洲人成网站在线播放欧美日韩| 制服丝袜大香蕉在线| 一级毛片精品| 99精品欧美一区二区三区四区| 国产成人啪精品午夜网站| 国产av在哪里看| 精品熟女少妇八av免费久了| 黑人巨大精品欧美一区二区mp4| 少妇裸体淫交视频免费看高清| 看片在线看免费视频| 亚洲av电影在线进入| 中文字幕久久专区| 久久精品国产清高在天天线| 人妻久久中文字幕网| 午夜福利高清视频| 国产精品精品国产色婷婷| 伦理电影免费视频| 日本与韩国留学比较| 久久亚洲真实| 中文字幕av在线有码专区| 精品国产乱码久久久久久男人| 精品久久久久久久人妻蜜臀av| 99热只有精品国产| 久久久久精品国产欧美久久久| 丁香欧美五月| 美女 人体艺术 gogo| 悠悠久久av| tocl精华| 久久香蕉精品热| 精品久久久久久久末码| 久久香蕉精品热| 免费看十八禁软件| 欧美激情久久久久久爽电影| 黄色视频,在线免费观看| 人人妻人人看人人澡| 久久久久久国产a免费观看| 欧美大码av| 1024香蕉在线观看| 特大巨黑吊av在线直播| 亚洲国产精品999在线| 亚洲一区高清亚洲精品| 一区福利在线观看| 日本三级黄在线观看| 国产高清激情床上av| 亚洲avbb在线观看| 国产真人三级小视频在线观看| 成人鲁丝片一二三区免费| 2021天堂中文幕一二区在线观| aaaaa片日本免费| 欧美+亚洲+日韩+国产| 一进一出抽搐动态| 精品电影一区二区在线| 亚洲人与动物交配视频| 巨乳人妻的诱惑在线观看| 中出人妻视频一区二区| 啦啦啦免费观看视频1| 欧美一级毛片孕妇| 国产一区在线观看成人免费| 精品久久蜜臀av无| 五月伊人婷婷丁香| 中出人妻视频一区二区| 丁香欧美五月| 一卡2卡三卡四卡精品乱码亚洲| 精华霜和精华液先用哪个| 色尼玛亚洲综合影院| 丰满人妻一区二区三区视频av | 少妇裸体淫交视频免费看高清| 久久久久免费精品人妻一区二区| 90打野战视频偷拍视频| 亚洲av片天天在线观看| 俄罗斯特黄特色一大片| 成人永久免费在线观看视频| 欧美日韩乱码在线| 亚洲欧美精品综合一区二区三区| 久久精品91无色码中文字幕| 九九在线视频观看精品| 黄片大片在线免费观看| 三级国产精品欧美在线观看 | 日本黄色视频三级网站网址| 午夜免费激情av| 午夜福利在线观看免费完整高清在 | 欧美日韩乱码在线| 成人高潮视频无遮挡免费网站| cao死你这个sao货| 变态另类成人亚洲欧美熟女| 一本久久中文字幕| 999精品在线视频| 一进一出抽搐动态| 精品人妻1区二区| 午夜成年电影在线免费观看| 在线免费观看的www视频| 啦啦啦观看免费观看视频高清| 国产一区二区激情短视频| 黑人操中国人逼视频|