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

    基于零控脫靶量的衛(wèi)星末端追逃控制方法

    2016-12-19 02:59:22王強葉東范寧軍吳炎烜
    北京理工大學學報 2016年11期
    關鍵詞:航天器矢量控制策略

    王強, 葉東, 范寧軍, 吳炎烜

    (1.北京理工大學 機電學院,北京 100081;2.國防知識產(chǎn)權局,北京 100088;3.哈爾濱工業(yè)大學 衛(wèi)星技術研究所,黑龍江,哈爾濱 150001)

    ?

    基于零控脫靶量的衛(wèi)星末端追逃控制方法

    王強1,2, 葉東3, 范寧軍1, 吳炎烜1

    (1.北京理工大學 機電學院,北京 100081;2.國防知識產(chǎn)權局,北京 100088;3.哈爾濱工業(yè)大學 衛(wèi)星技術研究所,黑龍江,哈爾濱 150001)

    針對攔截衛(wèi)星和目標航天器都可以采用幅值受限連續(xù)推力進行軌道機動時的衛(wèi)星末端攔截控制問題,由于對抗雙方都會采用對自己有利的控制策略來實現(xiàn)追擊和逃逸,因此將該問題可以看作追蹤逃逸問題. 由于兩衛(wèi)星距離較近,因此將非線性攔截逃逸相對動力學簡化為CW方程. 根據(jù)攔截任務終止要求引入零控脫靶矢量將動力學方程降階,采用攔截脫靶量和燃料消耗作為二次最優(yōu)目標函數(shù),推導了衛(wèi)星軌道次優(yōu)控制策略. 仿真表明該控制方法可以在目標衛(wèi)星具有機動逃逸能力時仍能成功將其攔截.

    衛(wèi)星追逃;連續(xù)推力;零控脫靶量;燃料子優(yōu)控制

    衛(wèi)星作為航天器的主體為國家安全提供重要的保障,它可以為指揮者制定作戰(zhàn)策略提供有效戰(zhàn)場信息,為導彈、飛機等傳遞制導導航指令,為地面其他作戰(zhàn)單元提供環(huán)境信息和作戰(zhàn)命令. 因此采用攜帶戰(zhàn)斗部或軟殺傷武器的攔截衛(wèi)星,在戰(zhàn)爭時期對于戰(zhàn)場上方或者本國領土上方空域的偵查衛(wèi)星和預警衛(wèi)星進行攔截可以保證己方在作戰(zhàn)中具有優(yōu)勢.

    當攔截衛(wèi)星和目標之間的距離很近時,認為衛(wèi)星攔截過程進入到終端攔截,此時攔截衛(wèi)星自身攜帶的電子偵測設備可以識別并定位目標衛(wèi)星. 學者們對于目標衛(wèi)星無法進行機動時攔截衛(wèi)星的控制策略進行充分的研究. 高曉光等[1]基于改進高斯法(IGM)和遺傳算法(GA)的混合優(yōu)化算法,為解決空間攔截軌道燃料消耗和轉移時間的綜合最優(yōu)問題,提出一種空間攔截軌道設計方法. Lu S等[2]針對目標航天器運行軌道為橢圓且星間無通信時衛(wèi)星交會問題,采用Lawden方程作為動力學模型,提出了一種幅值有限連續(xù)推力交會的自適應軌道控制方法. 變推力發(fā)動機可進行閉環(huán)控制,精度極高. Singla等[3]針對星載設備測量不確定性,根據(jù)輸出反饋控制提出了結構模型參考自適應控制器,并將其應用在衛(wèi)星自動交會過程中.

    目標軌道靜止時衛(wèi)星交會和攔截的研究技術已經(jīng)在航天工程領域得到應用. 但對于攔截問題,更嚴格的場景中目標航天器也可以實施軌道機動對攔截器進行躲避. 此時該問題可以看作衛(wèi)星非合作交會對接,鄔樹楠等[4]針對快速接近非合作目標航天器任務,提出了一種基于觀測器的相對軌道有限時間控制算法. Gao H等[5]針對衛(wèi)星交會中存在的參數(shù)不確定,未知干擾和控制輸入受限等問題,基于Lyapunov穩(wěn)定理論提出了一種H1狀態(tài)反饋魯棒控制器. 作者將多個控制目標,轉化為線性矩陣不等式,求解得到有效控制參數(shù).

    上述方法在求解控制策略時已經(jīng)考慮了目標的機動能力,但它們并沒有考慮攔截和目標衛(wèi)星機動的最優(yōu)性. 更合理的方法將該問題看作零和微分對策問題,Iaacs首先在該領域展開研究,并將對策雙方都選擇對自己有利的控制策略時的解定義為鞍點解[6]. Menon[7]利用反步法和二次目標函數(shù)得到控制輸出具有飽和特性時航天器追逃問題的反饋控制策略,然而其仿真中所用推力加速度很大不具有工程意義. Mauro等將攔截問題轉化為TPBVP,提出半直接法來求解鞍點解. 該方法被應用在兩航天器遠距離攻防攔截逃[8]問題中. Stupik[9]采用Kriging方法在CW方程的基礎上研究了衛(wèi)星攻防實時反饋解的求取,文章利用打靶法計算不同初值時的雙邊最優(yōu)解,并將這些最優(yōu)解進行插值獲得初態(tài)域和干擾有限時的次優(yōu)反饋控制策略.

    上述方法求解得到的是雙邊最優(yōu)控制開環(huán)解,不能直接為攔截衛(wèi)星提供有效的實時反饋控制. 本文將研究攔截衛(wèi)星和目標航天器都采用幅值受限連續(xù)推力進行軌道機動時的末端攔截問題,由于攔截衛(wèi)星和目標會采用對自己有利的控制策略,因此采用微分對策理論進行分析. 兩衛(wèi)星相對距離較近,因此非線性相對動力學模型將被簡化為CW方程. 通過引入零控脫靶矢量對系統(tǒng)降維,聯(lián)合燃料消耗和攔截脫靶量構建對策目標函數(shù),推導了燃料次優(yōu)反饋控制策略. 通過仿真分析了所提出方法的優(yōu)缺點.

    1 相對動力學模型

    如圖1所示,在攔截衛(wèi)星P附近建立虛擬衛(wèi)星O軌道,以O為原點建立相對軌道坐標系Oxyz,其中Ox軸指向衛(wèi)星O的位置矢量r1方向,Oz軸沿軌道角動量矢量方向,Oy軸在軌道與Ox,Oz成右手系. 假設:① 衛(wèi)星P為質點模型;② 虛擬衛(wèi)星O無機動;③ 忽略地球非球形攝動,大氣阻力攝動,太陽光壓和其他天體產(chǎn)生的引力攝動項的影響. 最后兩項假設可以保證參考衛(wèi)星軌道的形狀不發(fā)生變化.

    令攔截衛(wèi)星P在慣性坐標系下的位置矢量為r2,可得到兩星在赤道慣性坐標系下的動力學方程為

    (1)

    (2)

    虛擬衛(wèi)星軌道坐標系為旋轉動坐標系,則由矢量求導關系可得

    (3)

    式中:δr′表示δr的相對導數(shù);Ω為虛擬軌道角速度.

    (4)

    式中:u=[uxuyuz]為攔截衛(wèi)星推力加速度矢量;

    式(4)稱為CW方程,其線性形式給衛(wèi)星交會和攔截的求解帶來了方便. 文獻[4-5]中給出的軌道動力學模型假定目標航天器靜止或加速度極小,因此其動力學模型與本文有區(qū)別.

    (5)

    其中Φ11,Φ12,Φ21,Φ22可參考文獻[9]給出.

    假設攔截衛(wèi)星P和目標E相對距離遠小于兩星質心到地心的距離,則可在其附近設置虛擬衛(wèi)星,分別對攔截衛(wèi)星和目標,在該虛擬衛(wèi)星軌道坐標系下建立相對軌道運動學方程,可得如下式

    (6)

    式中xP,E分別為攔截衛(wèi)星和目標在虛擬軌道坐標系下的狀態(tài).uP,E為兩星的推力加速度. 它們滿足如式(7)所示幅值限制

    (7)

    為分析簡單,將攔截衛(wèi)星和目標的相對位置和速度作為新狀態(tài)xPE=xP-xE,對xPE求導可得

    (8)

    其中C=-B.

    z方向運動與軌道面內(nèi)運動解耦,后面的控制器設計和分析只針對軌道面內(nèi)運動,此時系統(tǒng)方程與式(4)相似,但參數(shù)矩陣發(fā)生變化.

    2 衛(wèi)星追蹤逃逸的次優(yōu)控制策略

    由于衛(wèi)星自身攜帶燃料的約束,衛(wèi)星攔截時間有限. 而在追蹤逃逸過程中,攔截衛(wèi)星與目標圍繞攔截結束時的距離展開爭奪,因此只考慮兩星的位置矢量,引入零控脫靶矢量(ZEM)來將系統(tǒng)的維度降低[10]

    Z(t)=DΦ(tf,t)xPE(t)=

    (9)

    式中:tf為衛(wèi)星攔截結束時間;Φ(tf,t)為xy方向動力學系統(tǒng)從時間t運動到tf時的狀態(tài)轉移矩陣;Z表示系統(tǒng)從t時刻無控轉移到末端時的脫靶矢量,根據(jù)系統(tǒng)理論可得零控脫靶矢量的狀態(tài)轉移矩陣滿足

    (10)

    對零控脫靶矢量關于時間求導可直接得到動力學方程

    (11)

    式中系數(shù)矩陣BP=-CE=Φ12(tf,t).

    在衛(wèi)星追蹤逃逸過程中,攔截衛(wèi)星和目標盡可能使脫靶量對自己有利,因此定義脫靶量為終端目標函數(shù);為延長攔截衛(wèi)星作用時間,得到反饋控制輸出,因此附加燃料消耗項作為過程控制目標函數(shù). 最終的目標函數(shù)為

    (12)

    式中:Q=qI2>0為二次系數(shù)矩陣;RP和RE都為正定矩陣. 而γ依賴攔截衛(wèi)星和目標的推力,定義了前者機動能力優(yōu)于后者的程度. 系統(tǒng)的哈密頓函數(shù)如下式.

    (13)

    式中λ為系統(tǒng)狀態(tài)的協(xié)態(tài)變量.

    為得到反饋控制器,先假設兩星的控制輸出滿足式(7)給出的幅值約束. 對系統(tǒng)哈密頓函數(shù)使用極大極小值可得攔截衛(wèi)星和目標航天器的控制輸出為

    系統(tǒng)的共軛方程和協(xié)態(tài)變量終端條件如下

    (15)

    則有λ=λ(tf)=QZ(tf). 將兩星控制輸出帶入到系統(tǒng)動力學方程中可得

    對式(16)積分,并帶入?yún)f(xié)態(tài)變量整理后可得終端時刻的零控脫靶矢量為

    Z(tf)=F-1(tf,t)Z(t).

    (17)

    其中

    令RP=RE=RI2,Q=qI4,定義攔截剩余時間為tgo=tf-t,因此F(tf,t)可以如下給出

    (18)

    將式(18)帶入到(14)中,可得反饋控制的顯式形式為

    uP=-q/RF-1Z, uE=-q/(Rγ2)F-1Z.

    (19)

    該控制策略具有最優(yōu)性[11]. 系統(tǒng)閉環(huán)軌跡可將控制輸出帶入方程中積分得到.

    若預期控制指令幅值超過推力上限,則無法得到最優(yōu)控制. 下面給出一個可以實用的次優(yōu)控制策略.

    對于變量x定義它的飽和函數(shù)為[12]

    對于哈密頓函數(shù),利用函數(shù)取極大極小值的必要條件得到兩星最優(yōu)控制輸出,并將其帶入到狀態(tài)方程,可以得到

    此時系統(tǒng)共軛方程和橫截條件與式(15)相同,對式(21)積分,可以得到終端時刻零控脫靶矢量滿足式(22)

    (22)

    式中Z(tf)為飽和函數(shù)的自變量,無法得到解析Z(t),為了給出反饋解,依然假設等式成立. 此時攔截衛(wèi)星和目標航天器的控制輸出可如下給出

    綜上,考慮目標航天器的機動能力時,兩星的反饋控制策略可如式(23)給出,明顯只有當攔截衛(wèi)星和目標的預定控制輸出滿足式(7)的幅值約束時,該策略為最優(yōu)控制,否則非最優(yōu).

    3 仿真和控制性能分析

    假設攔截衛(wèi)星和目標航天器都在地球同步軌道周圍運動,則將虛擬衛(wèi)星選為GEO上軌道與這兩星相近的衛(wèi)星. 軌道角速度為ω=7.272 2×10-5rad/s. 攔截和目標衛(wèi)星的推力加速度幅值上限分別為ρP=0.686 m/s2和ρE=0.343 m/s2,分別對兩組算例進行仿真.

    ① 攔截衛(wèi)星和目標初始時刻的狀態(tài)分別為

    [-2,-1,0.005,0.001]T,[0,0,0,0]T,即目標衛(wèi)星初始時刻與虛擬衛(wèi)星位置相同,其中位置與速度的單位分別為km和km/s. 其他參數(shù)為R=1,Q=1,γ=2. 對攔截結束時間在200~2 000 s變化的攔截逃逸場景進行仿真,統(tǒng)計兩星攔截結束時的距離,可得圖2,可以看到,脫靶量隨攔截結束時間增大而減小,只有當攔截結束時間大于500 s時,脫靶量才能小于50 m. 當攔截時間為200 s時,脫靶量為149.89 m,大于100 m,可認為攔截失敗. 當攔截時間為600 s時,脫靶量為32.49 m. 此時兩星的相對距離如圖3所示. 圖4給出兩星的攔截逃逸軌跡. 可以看到隨著剩余時間減少,兩星的相對位移逐漸減少,在攔截時刻,y方向的相對位移收斂到較小值,而x方向仍有較小位置差異.

    圖5和圖6給出了600 s時兩星在攔截逃逸過程中的推力加速度的變化曲線. 可以看到,兩星在初始時刻都不采用幅值最大的控制策略,直到攔截快要結束時,攔截衛(wèi)星最先開始使用最大幅值輸出控制,而目標衛(wèi)星的控制輸出也逐漸增大到最大幅值.

    ② 攔截衛(wèi)星和目標初始狀態(tài)分別為[-2.000,1.000,0.005,0]T和[2.000,-1.000,-0.005,0]T,即兩星關于狀態(tài)空間零點對稱. 攔截結束時間為1 000 s,可得攔截過程中兩星相對位置的變化曲線和兩星軌跡如圖7和圖8所示.

    此時由于攔截衛(wèi)星初始時刻具有x正方向的速度,因此會先向右運動,而隨后由于目標衛(wèi)星向左運動,又向右運動追擊目標. 圖9和圖10給出了兩星的加速度曲線,明顯此時攔截衛(wèi)星和目標在攔截結束時刻的加速度幅值也迅速變大,直到達到上限.

    上面兩個仿真工況中,攔截衛(wèi)星和目標航天器的控制矢量實時同向,這可以由式(23)得到,具體由于式(20)給出的飽和函數(shù)只改變自變量的幅值大小,而保留了它的方向性質,因此攔截衛(wèi)星和目標的控制向量實時共線,同向可以從式(23)中的負號以及BP=-CE得到.

    4 結 論

    研究了攔截衛(wèi)星和目標航天器都采用幅值有限連續(xù)推力進行變軌時的衛(wèi)星末端攔截軌道控制問題. 由于兩衛(wèi)星相對距離較近,因此可以在其附近建立虛擬衛(wèi)星軌道,進而簡化相對運動方程得到攔截逃逸動力學模型. 通過引入零控脫靶矢量可以將系統(tǒng)方程降維. 利用微分對策理論推導了系統(tǒng)的鞍點平衡解,最后得到了攔截衛(wèi)星在考慮目標最優(yōu)機動下的次優(yōu)反饋控制律. 仿真表明了該控制策略的有效性,并表明采用該策略時,攔截結束時間增加,攔截脫靶量減小;攔截衛(wèi)星和目標在追蹤逃逸過程中任何時刻的控制推力同向,而推力幅值剛開始時變化較小,攔截結束前迅速增加,最后達到衛(wèi)星所能提供的最大推力.

    [1] 高曉光,湯洪,端軍紅.基于混合優(yōu)化算法的空間攔截軌道優(yōu)化設計[J].北京航空航天大學學報,2015,41(9):1574-1581.

    Gao Xiaoguang, Tang Hong, Duan Junhong. Space interception orbit optimization design based on hybrid optimal algorithm[J]. Journal of Beijing University of Aeronautics and Astronautics, 2015,41(9):1574-1581. (in Chinese)

    [2] Lu S, Xu S. Adaptive control for autonomous rende-zvous of spacecraft on elliptical orbit[J]. Acta Mecha-nica Sinica, 2009,25(4):539-545.

    [3] Singla Puneet, Subbarao Kamesh, Junkins John L. Adaptive output feedback control for spacecraft rendezvous and docking under measurementuncer-tainty[J]. Journal of Guidance, Control, and Dynamics, 2006,29(4):892-902.

    [4] 鄔樹楠,吳國強,孫兆偉.接近非合作目標的航天器相對軌道有限時間控制[J].大連理工大學學報,2013,53(6):885-892.

    Wu Shunan, Wu Guoqiang, Sun Zhaowei. Spacecraft relative orbit finite-time control for proximity to non-cooperative strategy[J]. Journal of Dalian University of Technology, 2013,53(6):885-892. (in Chinese)

    [5] Gao H, Yang X, Shi P. Multi-objective robust control of spacecraft rendezvous[J]. IEEE Transactions on Control Systems Technology, 2009,17(4):794-802.

    [6] Isaacs R. Differential games: a mathematical theory with applications to warfare and pursuit, control and optimization[M]. [S.l.]: Courier Dover Publications, 1999.

    [7] Menon P K A, Calise A J, Leung S K M. Guidance laws for spacecraft pursuit-evasion and rendezvous[C]∥AIAA Guidance Navigation and Control. [S.l.]: AIAA, 1988:668-697.

    [8] Conway B A, Pontani M. Numerical solution of the three-dimensional orbital pursuit-evasion game[J]. Journal of Guidance, Control, and Dynamics, 2009,32(2):474-487.

    [9] Stupik J, Pontani M, Conway B. Optimal pursuit-evasion spacecraft trajectories in the hill reference frame[C]∥Proceedings of AIAA/AAS Astrodynamics Specialist Con-ference. Reston, VA: AIAA, 2012:4882.

    [10] Rubinsky S, Gutman S. Vector guidance approach to three-player conflict in exoatmospheric interception[J]. Journal of Guidance, Control, and Dynamics, 2015,38(12):2270-2286.

    [11] Bardhan R, Ghose D. Intercepting maneuvering target with specified impact angle by modified SDRE technique[C]∥Proceedings of 292 American Control Conference (ACC). [S.l.]: IEEE, 2012:4613-4618.

    [12] Rusnak I, Meir L. Optimal guidance for high-order and acceleration constrained missile[J]. Journal of Guidance Control & Dynamics, 1991,14(3):589-596.

    (責任編輯:劉雨)

    Terminal Orbital Control of Satellite Pursuit Evasion Game Based on Zero Effort Miss

    WANG Qiang1,2, YE Dong3, FAN Ning-jun1, WU Yan-xuan1

    (1.School of Electromechanical, Beijing Institute of Technology, Beijing 100081, China;2.National Defense Intellectual Property Office, Beijing 100088, China;3.Institute of Satellite Technology, Harbin Institute of Technology, Harbin, Heilongjiang 150001, China)

    In the satellite interception control strategies, both the interceptor and target could perform orbital maneuver with limited magnitude continuous thrust. Basically, this problem was regarded as a pursuit-evasion game since satellites in both sides would try their best to catch or escape. Assumption that the two player were close such that the nonlinear interception dynamics could be reduced to the CW equations. Zero effort miss was introduced to simplify the system. The controller was constructed by adding the interception miss and the fuel consumption to the optimal objective. Simulation was conducted to analyze the advantage and weakness of proposed control strategy.

    satellite pursuit-evasion game; continuous thrust; zero effort miss; fuel suboptimal control

    2015-11-17

    中央高校基本科研業(yè)務費專項資金資助項目(HIT.NSRIF.2015033);微小型航天器技術國防重點學科實驗室開放基金資助項目(HIT.KLOF.MST.201501)

    王強(1977—),男,博士生,E-mail:yeboyansha@163.com.

    葉東(1985—),男,博士,講師,E-mail:yed@hit.edu.cn.

    V 27

    A

    1001-0645(2016)11-1171-06

    10.15918/j.tbit1001-0645.2016.11.014

    猜你喜歡
    航天器矢量控制策略
    2022 年第二季度航天器發(fā)射統(tǒng)計
    國際太空(2022年7期)2022-08-16 09:52:50
    矢量三角形法的應用
    考慮虛擬慣性的VSC-MTDC改進下垂控制策略
    能源工程(2020年6期)2021-01-26 00:55:22
    2019 年第二季度航天器發(fā)射統(tǒng)計
    國際太空(2019年9期)2019-10-23 01:55:34
    工程造價控制策略
    山東冶金(2019年3期)2019-07-10 00:54:04
    2018 年第三季度航天器發(fā)射統(tǒng)計
    國際太空(2018年12期)2019-01-28 12:53:20
    2018年第二季度航天器發(fā)射統(tǒng)計
    國際太空(2018年9期)2018-10-18 08:51:32
    現(xiàn)代企業(yè)會計的內(nèi)部控制策略探討
    消費導刊(2018年10期)2018-08-20 02:57:02
    基于矢量最優(yōu)估計的穩(wěn)健測向方法
    容錯逆變器直接轉矩控制策略
    波多野结衣高清作品| 中文在线观看免费www的网站 | 宅男免费午夜| 精品久久久久久久末码| 夜夜夜夜夜久久久久| 精品国产乱子伦一区二区三区| 国产精品免费视频内射| 男女做爰动态图高潮gif福利片| 国产一区二区三区在线臀色熟女| 好男人在线观看高清免费视频| 国产久久久一区二区三区| 免费在线观看视频国产中文字幕亚洲| www.自偷自拍.com| 久久久久久国产a免费观看| 两性夫妻黄色片| 亚洲第一电影网av| 99国产精品一区二区三区| ponron亚洲| 国产91精品成人一区二区三区| 日韩av在线大香蕉| 一卡2卡三卡四卡精品乱码亚洲| 在线国产一区二区在线| av天堂在线播放| 亚洲精品美女久久av网站| 久久久久九九精品影院| 高清在线国产一区| 一二三四在线观看免费中文在| 欧美在线黄色| 亚洲欧美一区二区三区黑人| 一卡2卡三卡四卡精品乱码亚洲| 亚洲最大成人中文| 亚洲国产精品999在线| 久久久国产精品麻豆| 午夜福利在线观看吧| www国产在线视频色| 看黄色毛片网站| 国产99白浆流出| 男人舔女人下体高潮全视频| 黄色视频,在线免费观看| 热99re8久久精品国产| 脱女人内裤的视频| 国产精品综合久久久久久久免费| 国产av不卡久久| www.www免费av| 99热这里只有精品一区 | 三级男女做爰猛烈吃奶摸视频| 国产在线观看jvid| 亚洲av美国av| 悠悠久久av| 91老司机精品| 两个人看的免费小视频| 免费在线观看日本一区| 成人永久免费在线观看视频| 亚洲精品av麻豆狂野| 日本撒尿小便嘘嘘汇集6| 国产成人av教育| 国产成人欧美在线观看| www国产在线视频色| 高潮久久久久久久久久久不卡| 91在线观看av| 欧美三级亚洲精品| 国产亚洲精品第一综合不卡| 丰满的人妻完整版| 视频区欧美日本亚洲| 母亲3免费完整高清在线观看| 国产熟女午夜一区二区三区| 狂野欧美激情性xxxx| 九色国产91popny在线| 韩国av一区二区三区四区| 午夜激情av网站| 国产一区二区在线观看日韩 | 午夜视频精品福利| 哪里可以看免费的av片| 18禁黄网站禁片午夜丰满| 日日爽夜夜爽网站| 无遮挡黄片免费观看| 国内毛片毛片毛片毛片毛片| 成人亚洲精品av一区二区| 亚洲一区中文字幕在线| 免费av毛片视频| 欧美精品亚洲一区二区| 国产91精品成人一区二区三区| 国产精华一区二区三区| 国产成人av激情在线播放| 不卡一级毛片| 欧美日韩中文字幕国产精品一区二区三区| 久久久久久亚洲精品国产蜜桃av| 精品福利观看| 中文字幕人妻丝袜一区二区| 老司机午夜十八禁免费视频| 搞女人的毛片| av国产免费在线观看| 亚洲av熟女| 亚洲欧美精品综合久久99| 超碰成人久久| 亚洲美女黄片视频| 欧美性猛交╳xxx乱大交人| 国产熟女午夜一区二区三区| 中文字幕av在线有码专区| 亚洲av成人不卡在线观看播放网| 在线观看美女被高潮喷水网站 | 久久婷婷成人综合色麻豆| 亚洲精品美女久久久久99蜜臀| 久久香蕉激情| 岛国在线观看网站| 亚洲美女视频黄频| 国产av又大| 日日爽夜夜爽网站| 国产精品久久久av美女十八| 国产精品久久视频播放| 一级毛片女人18水好多| 2021天堂中文幕一二区在线观| 日日干狠狠操夜夜爽| 一卡2卡三卡四卡精品乱码亚洲| 一二三四在线观看免费中文在| 99热只有精品国产| 好男人在线观看高清免费视频| 亚洲av中文字字幕乱码综合| 亚洲真实伦在线观看| 国产伦人伦偷精品视频| 亚洲国产精品久久男人天堂| 国产精品亚洲一级av第二区| 亚洲熟妇熟女久久| 亚洲免费av在线视频| 国产精品一区二区三区四区免费观看 | 亚洲中文日韩欧美视频| 亚洲精品一区av在线观看| 日本免费a在线| 欧美日韩亚洲国产一区二区在线观看| 国产亚洲av高清不卡| 五月玫瑰六月丁香| x7x7x7水蜜桃| 欧美不卡视频在线免费观看 | 看片在线看免费视频| 国产av在哪里看| а√天堂www在线а√下载| 国产av一区在线观看免费| 精品国产亚洲在线| 久久人妻福利社区极品人妻图片| 成熟少妇高潮喷水视频| 999久久久精品免费观看国产| 久久久水蜜桃国产精品网| 久久这里只有精品19| 欧美乱色亚洲激情| 国产高清激情床上av| 午夜福利视频1000在线观看| 国产成人精品久久二区二区免费| 99re在线观看精品视频| 99riav亚洲国产免费| 国产av麻豆久久久久久久| 日韩欧美精品v在线| 免费搜索国产男女视频| 一级毛片精品| 中文字幕人成人乱码亚洲影| 1024手机看黄色片| 老司机深夜福利视频在线观看| 草草在线视频免费看| 午夜福利视频1000在线观看| 久久久久久久久中文| 天天躁狠狠躁夜夜躁狠狠躁| 嫁个100分男人电影在线观看| 可以在线观看毛片的网站| 人人妻人人澡欧美一区二区| 日韩高清综合在线| 日日摸夜夜添夜夜添小说| 一进一出抽搐gif免费好疼| 黄色视频,在线免费观看| 法律面前人人平等表现在哪些方面| 久久香蕉精品热| 亚洲中文av在线| 国产单亲对白刺激| 久久香蕉精品热| 久久久久久国产a免费观看| 精品欧美一区二区三区在线| 久久久久久免费高清国产稀缺| 国产成年人精品一区二区| av片东京热男人的天堂| 小说图片视频综合网站| 日韩三级视频一区二区三区| 欧美成人免费av一区二区三区| 在线观看66精品国产| 国产精品日韩av在线免费观看| 欧美色欧美亚洲另类二区| 亚洲欧美日韩东京热| 五月伊人婷婷丁香| 成人av在线播放网站| 欧美乱色亚洲激情| 床上黄色一级片| 成人国语在线视频| 亚洲 欧美 日韩 在线 免费| 免费观看人在逋| 天堂影院成人在线观看| 国产三级在线视频| 精品欧美国产一区二区三| 亚洲午夜理论影院| 免费电影在线观看免费观看| 狂野欧美白嫩少妇大欣赏| 91九色精品人成在线观看| 97碰自拍视频| 国产精品,欧美在线| 成人午夜高清在线视频| 欧美一级毛片孕妇| 伦理电影免费视频| 香蕉久久夜色| 男人舔女人的私密视频| 十八禁网站免费在线| 成在线人永久免费视频| 久久性视频一级片| 少妇熟女aⅴ在线视频| 99精品欧美一区二区三区四区| 热99re8久久精品国产| 国产精品一区二区精品视频观看| 黄色女人牲交| 亚洲真实伦在线观看| 久久草成人影院| 十八禁人妻一区二区| 午夜亚洲福利在线播放| 国产日本99.免费观看| xxx96com| 婷婷精品国产亚洲av| 97人妻精品一区二区三区麻豆| 国产精品av久久久久免费| 国产成人aa在线观看| 啦啦啦观看免费观看视频高清| 每晚都被弄得嗷嗷叫到高潮| 日韩欧美三级三区| 毛片女人毛片| 亚洲人成网站高清观看| 一级黄色大片毛片| 亚洲熟妇中文字幕五十中出| 亚洲国产精品久久男人天堂| 国产精品久久久久久亚洲av鲁大| 久久九九热精品免费| 国产熟女xx| 国产精品一及| 欧美黄色淫秽网站| 亚洲一区二区三区色噜噜| 国内毛片毛片毛片毛片毛片| 人妻丰满熟妇av一区二区三区| 999精品在线视频| 久久中文字幕一级| 欧美一级a爱片免费观看看 | 久久久久性生活片| 日本黄大片高清| 久久精品国产亚洲av高清一级| 一级毛片精品| 国产高清视频在线观看网站| 欧洲精品卡2卡3卡4卡5卡区| 午夜福利在线在线| 国产亚洲av嫩草精品影院| 国产午夜精品论理片| 久久久久久久久久黄片| 特大巨黑吊av在线直播| ponron亚洲| 色播亚洲综合网| 成人欧美大片| 18禁国产床啪视频网站| 老司机福利观看| 日韩欧美国产一区二区入口| 黄片大片在线免费观看| 婷婷亚洲欧美| 亚洲九九香蕉| 免费在线观看日本一区| 99精品欧美一区二区三区四区| 久久草成人影院| 国产精品综合久久久久久久免费| 香蕉国产在线看| 不卡av一区二区三区| 久久久水蜜桃国产精品网| 高清在线国产一区| 无遮挡黄片免费观看| 成年人黄色毛片网站| 日韩精品中文字幕看吧| 久久久久久久午夜电影| 男女那种视频在线观看| 日韩 欧美 亚洲 中文字幕| 亚洲成人久久爱视频| 又紧又爽又黄一区二区| 亚洲无线在线观看| 最新美女视频免费是黄的| 成人三级黄色视频| aaaaa片日本免费| 12—13女人毛片做爰片一| 亚洲国产欧美一区二区综合| 成人高潮视频无遮挡免费网站| a级毛片a级免费在线| 精品久久久久久久久久免费视频| 男女下面进入的视频免费午夜| xxxwww97欧美| 成人国产一区最新在线观看| 久久香蕉激情| 少妇粗大呻吟视频| 国产熟女午夜一区二区三区| 三级国产精品欧美在线观看 | 美女扒开内裤让男人捅视频| 免费看a级黄色片| 国产欧美日韩一区二区精品| 久久精品91蜜桃| 大型av网站在线播放| 中文资源天堂在线| 亚洲人成77777在线视频| 91老司机精品| 五月玫瑰六月丁香| 夜夜躁狠狠躁天天躁| 香蕉av资源在线| 精品欧美一区二区三区在线| 97人妻精品一区二区三区麻豆| 一级作爱视频免费观看| 超碰成人久久| 搡老岳熟女国产| 国产人伦9x9x在线观看| 人人妻人人看人人澡| 亚洲国产精品sss在线观看| 亚洲国产精品成人综合色| 熟女少妇亚洲综合色aaa.| 精品电影一区二区在线| 亚洲性夜色夜夜综合| 很黄的视频免费| 欧美 亚洲 国产 日韩一| 深夜精品福利| 国产成人av教育| 亚洲自偷自拍图片 自拍| 97人妻精品一区二区三区麻豆| 18禁裸乳无遮挡免费网站照片| 国产99久久九九免费精品| 天天躁狠狠躁夜夜躁狠狠躁| 桃红色精品国产亚洲av| 欧美高清成人免费视频www| 国产在线观看jvid| 男女视频在线观看网站免费 | 国产黄片美女视频| 九九热线精品视视频播放| 99精品久久久久人妻精品| 日韩精品中文字幕看吧| 九九热线精品视视频播放| 亚洲精品国产一区二区精华液| 每晚都被弄得嗷嗷叫到高潮| 亚洲一码二码三码区别大吗| 精品一区二区三区av网在线观看| 亚洲精品在线观看二区| 丝袜人妻中文字幕| 国产激情欧美一区二区| 日韩欧美精品v在线| 欧美在线一区亚洲| 淫妇啪啪啪对白视频| 最好的美女福利视频网| 在线观看午夜福利视频| 亚洲av中文字字幕乱码综合| 国产一区在线观看成人免费| 最近在线观看免费完整版| 99国产精品一区二区三区| 国产精品一区二区三区四区久久| 欧美丝袜亚洲另类 | 日韩欧美国产在线观看| 欧美+亚洲+日韩+国产| 黄频高清免费视频| 亚洲精品国产一区二区精华液| 国产伦在线观看视频一区| 久久午夜综合久久蜜桃| 亚洲av片天天在线观看| 2021天堂中文幕一二区在线观| а√天堂www在线а√下载| 首页视频小说图片口味搜索| av福利片在线观看| 午夜免费观看网址| 悠悠久久av| 国产一区二区激情短视频| 成人三级做爰电影| 一区二区三区激情视频| 久久久国产精品麻豆| 国产精品永久免费网站| 天天躁夜夜躁狠狠躁躁| 亚洲国产精品成人综合色| 国内精品久久久久精免费| 国产亚洲精品第一综合不卡| 国产欧美日韩一区二区精品| 变态另类成人亚洲欧美熟女| 中亚洲国语对白在线视频| 久热爱精品视频在线9| 亚洲一码二码三码区别大吗| 一级黄色大片毛片| 午夜福利视频1000在线观看| 老汉色∧v一级毛片| tocl精华| 两个人的视频大全免费| 亚洲欧美日韩高清专用| 国产一区二区三区视频了| 欧美乱妇无乱码| 欧美日韩黄片免| 久久午夜亚洲精品久久| 999精品在线视频| 成人av在线播放网站| 99久久久亚洲精品蜜臀av| 亚洲成人久久爱视频| 天天躁夜夜躁狠狠躁躁| 国产精品亚洲av一区麻豆| 亚洲成人精品中文字幕电影| 精品日产1卡2卡| 国产高清激情床上av| 正在播放国产对白刺激| 亚洲精品色激情综合| 久久国产精品人妻蜜桃| АⅤ资源中文在线天堂| 日日爽夜夜爽网站| 亚洲人成77777在线视频| netflix在线观看网站| 色综合站精品国产| 日韩欧美在线乱码| 露出奶头的视频| 日日爽夜夜爽网站| 成在线人永久免费视频| 不卡av一区二区三区| 国产精品久久久久久人妻精品电影| 欧美绝顶高潮抽搐喷水| 精品国产乱码久久久久久男人| 亚洲国产欧美网| 午夜亚洲福利在线播放| 亚洲成av人片在线播放无| 99国产综合亚洲精品| 亚洲国产日韩欧美精品在线观看 | 18禁黄网站禁片免费观看直播| 国产精品亚洲av一区麻豆| 丁香欧美五月| 亚洲天堂国产精品一区在线| 长腿黑丝高跟| 女人爽到高潮嗷嗷叫在线视频| 俄罗斯特黄特色一大片| 亚洲国产精品合色在线| 又大又爽又粗| 身体一侧抽搐| 国产成人aa在线观看| 欧美日本亚洲视频在线播放| 国产成人啪精品午夜网站| 久久天堂一区二区三区四区| 日日爽夜夜爽网站| 丝袜美腿诱惑在线| 一二三四在线观看免费中文在| 日韩大码丰满熟妇| 大型黄色视频在线免费观看| 午夜激情福利司机影院| 美女大奶头视频| 老司机福利观看| 一进一出好大好爽视频| 免费在线观看日本一区| 精品国产乱码久久久久久男人| 激情在线观看视频在线高清| 男女视频在线观看网站免费 | 香蕉国产在线看| 老熟妇仑乱视频hdxx| 一级作爱视频免费观看| 最好的美女福利视频网| 国产单亲对白刺激| 国产精品一及| 国产激情偷乱视频一区二区| 亚洲精品中文字幕一二三四区| 欧美日韩精品网址| 毛片女人毛片| 神马国产精品三级电影在线观看 | 黄片小视频在线播放| 欧美最黄视频在线播放免费| 欧美性猛交╳xxx乱大交人| 少妇熟女aⅴ在线视频| 日韩高清综合在线| 亚洲成人久久性| 亚洲免费av在线视频| 男人舔女人下体高潮全视频| 国产一区二区激情短视频| 成人欧美大片| 九九热线精品视视频播放| 又黄又粗又硬又大视频| 1024视频免费在线观看| 18禁观看日本| 欧美3d第一页| 国产私拍福利视频在线观看| 国产精品野战在线观看| 俄罗斯特黄特色一大片| 精品不卡国产一区二区三区| 好男人在线观看高清免费视频| av中文乱码字幕在线| 可以免费在线观看a视频的电影网站| 精品国产美女av久久久久小说| 男女那种视频在线观看| 国产精品av久久久久免费| 成人永久免费在线观看视频| 午夜福利在线观看吧| 亚洲成人免费电影在线观看| 99精品久久久久人妻精品| 亚洲欧美日韩无卡精品| 一卡2卡三卡四卡精品乱码亚洲| 国产精品免费一区二区三区在线| 日日夜夜操网爽| 免费看美女性在线毛片视频| 变态另类成人亚洲欧美熟女| 三级毛片av免费| www日本黄色视频网| 91老司机精品| 亚洲九九香蕉| 成人特级黄色片久久久久久久| 女警被强在线播放| 国产成人系列免费观看| 91老司机精品| 精品国产乱码久久久久久男人| 777久久人妻少妇嫩草av网站| 亚洲成av人片免费观看| 最好的美女福利视频网| 伦理电影免费视频| 嫩草影院精品99| 国内精品一区二区在线观看| 两性夫妻黄色片| 色在线成人网| 国产又黄又爽又无遮挡在线| 99国产综合亚洲精品| av免费在线观看网站| 久久久国产成人免费| 亚洲自拍偷在线| 国产高清videossex| 久久久久久国产a免费观看| 曰老女人黄片| 欧美日韩一级在线毛片| 在线观看一区二区三区| 免费一级毛片在线播放高清视频| 亚洲av熟女| 少妇粗大呻吟视频| 欧美成人免费av一区二区三区| 亚洲精品美女久久久久99蜜臀| 欧美 亚洲 国产 日韩一| 男插女下体视频免费在线播放| 人人妻人人澡欧美一区二区| 亚洲欧洲精品一区二区精品久久久| 琪琪午夜伦伦电影理论片6080| 亚洲自偷自拍图片 自拍| 给我免费播放毛片高清在线观看| 国产成人系列免费观看| 亚洲av成人精品一区久久| 欧美日韩乱码在线| 在线观看www视频免费| 欧美中文日本在线观看视频| 88av欧美| 757午夜福利合集在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 在线视频色国产色| 免费在线观看日本一区| 成人午夜高清在线视频| 亚洲专区中文字幕在线| 黄色丝袜av网址大全| 午夜福利免费观看在线| 黄色成人免费大全| 一进一出抽搐动态| 一进一出好大好爽视频| 午夜日韩欧美国产| 欧美成人一区二区免费高清观看 | a级毛片a级免费在线| av免费在线观看网站| 免费在线观看亚洲国产| 国产av一区二区精品久久| 波多野结衣巨乳人妻| 人成视频在线观看免费观看| 国内少妇人妻偷人精品xxx网站 | 麻豆一二三区av精品| 亚洲无线在线观看| 夜夜爽天天搞| 欧美乱色亚洲激情| а√天堂www在线а√下载| 亚洲自偷自拍图片 自拍| 这个男人来自地球电影免费观看| 999久久久精品免费观看国产| 国产熟女午夜一区二区三区| 天堂√8在线中文| 国产av麻豆久久久久久久| 色综合站精品国产| 国产熟女xx| 母亲3免费完整高清在线观看| 日韩欧美国产一区二区入口| 99精品久久久久人妻精品| 欧美+亚洲+日韩+国产| 亚洲黑人精品在线| 又爽又黄无遮挡网站| 美女免费视频网站| 色av中文字幕| 欧美午夜高清在线| 国产成人av激情在线播放| 波多野结衣高清作品| 香蕉久久夜色| 亚洲av第一区精品v没综合| 精品福利观看| 桃红色精品国产亚洲av| 每晚都被弄得嗷嗷叫到高潮| x7x7x7水蜜桃| 国产三级在线视频| 亚洲一区高清亚洲精品| 无遮挡黄片免费观看| e午夜精品久久久久久久| 午夜视频精品福利| 桃色一区二区三区在线观看| 亚洲乱码一区二区免费版| 99久久精品国产亚洲精品| 97人妻精品一区二区三区麻豆| 日韩欧美在线二视频| 男女床上黄色一级片免费看| 免费在线观看日本一区| 欧美一区二区国产精品久久精品 | 天天一区二区日本电影三级| 麻豆av在线久日| 女同久久另类99精品国产91| 毛片女人毛片| 嫩草影视91久久| 成人国语在线视频| 久久亚洲精品不卡| 欧美久久黑人一区二区| 亚洲精品美女久久久久99蜜臀| 看黄色毛片网站| 亚洲七黄色美女视频| 一边摸一边做爽爽视频免费| 国产成人影院久久av| 一边摸一边做爽爽视频免费|