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

    衛(wèi)星編隊連續(xù)推力控制的燃料平衡方法

    2016-04-25 09:12:10段晨陽趙亞飛張世杰孔憲仁
    航天器工程 2016年1期
    關鍵詞:主星相位角編隊

    段晨陽 趙亞飛 張世杰 孔憲仁

    (哈爾濱工業(yè)大學衛(wèi)星技術研究所,哈爾濱 150080)

    衛(wèi)星編隊連續(xù)推力控制的燃料平衡方法

    段晨陽 趙亞飛 張世杰 孔憲仁

    (哈爾濱工業(yè)大學衛(wèi)星技術研究所,哈爾濱 150080)

    在衛(wèi)星編隊飛行中,編隊重構等機動過程會導致整個編隊衛(wèi)星之間燃料消耗不均勻,甚至出現(xiàn)某一成員衛(wèi)星燃料消耗完,而導致整個編隊構型提前結束乃至任務失敗。針對該問題,文章提出了在衛(wèi)星編隊軌道重構過程中可采用的一種燃料平衡方法,即基于連續(xù)推力控制,以燃料最優(yōu)為控制目標,通過建立燃料消耗函數(shù),推導了不同相位角及重構半徑時的最優(yōu)控制加速度,通過減小各從星之間的燃料消耗函數(shù)的差異,使得不同成員衛(wèi)星燃料消耗差別最小。編隊衛(wèi)星燃料平衡程度取決于初始相位角,文章給出了最佳初始相位角的表達式。最后,對以一主二從的三星編隊在從星軌道重構中的從星燃料平衡問題進行了仿真,分別驗證了衛(wèi)星編隊連續(xù)推力控制方法和編隊衛(wèi)星燃料平衡方法的正確性和有效性。

    衛(wèi)星編隊;軌道重構;燃料平衡;連續(xù)推力

    1 引言

    衛(wèi)星編隊飛行是人類對天空探索及實現(xiàn)衛(wèi)星協(xié)同工作的重要技術。通過編隊飛行,可實現(xiàn)干涉測量[12]、立體成像[3-4]、空間監(jiān)測等多種任務。通過衛(wèi)星間的協(xié)調(diào)工作,可實現(xiàn)單星無法完成的任務,對未來宇宙探索和空間服務具有重要意義。隨著衛(wèi)星推進技術的發(fā)展,衛(wèi)星推進方式實現(xiàn)了從化學推進—化學/電推進—電推進的過渡,在未來編隊飛行任務中,電推進等新型推進技術有助于減少衛(wèi)星對燃料的依賴,具有廣闊的發(fā)展前景。在ESA發(fā)布的下一代衛(wèi)星平臺(Next-generation Satellite Platform)研究中,電推進方式以其高效節(jié)能的特點被列為重點發(fā)展目標之一。

    編隊飛行中各成員衛(wèi)星之間協(xié)同工作,在軌以某種構型存在,實現(xiàn)任務的需求和功能。編隊飛行與單星飛行相比,可實現(xiàn)的功能更為復雜。但是,編隊飛行也有其自身的缺點,如軌道機動導致衛(wèi)星之間燃料消耗不一致,最終可能使得某一成員衛(wèi)星燃料耗盡,進而導致編隊性能降低,構型提前結束,甚至無法完成預期的在軌任務。美國、俄羅斯等國家已經(jīng)開始了對編隊衛(wèi)星重構中的燃料平衡問題的研究[5]。一些學者提出在工程上采用控制分配方法,約束軌道重構中成員衛(wèi)星之間的燃料消耗[6-8]。Susan A.Frost將燃料平衡指標轉(zhuǎn)化為范數(shù)最優(yōu)問題,達到了均衡燃料消耗的目的[9-10]。在重力場和深空環(huán)境兩種約束下,Amirreza R.等對雙星和三星編隊機動中的燃料平衡問題進行了仿真驗證[11]。Beard R.W.等開展了相關的工作,為燃料平衡問題的研究提供了一定的理論依據(jù)[12]。

    本文對采用連續(xù)推力軌道控制模式下的編隊重構的燃料平衡問題進行研究,給出了響應的理論和方法,可供未來編隊衛(wèi)星的燃料平衡控制參考。

    2 編隊燃料平衡優(yōu)化設計

    在確定燃料平衡問題模型之前,首先須要定義衛(wèi)星編隊重構中采用的坐標系,并介紹衛(wèi)星編隊重構的基本原理。在此基礎上,逐步分析編隊燃料平衡優(yōu)化設計中涉及到的關鍵技術。

    2.1 坐標系定義

    本節(jié)主要介紹本文涉及的兩種坐標系:慣性坐標系和軌道坐標系,如圖1所示。

    (1)慣性坐標系:原點O位于地心,OX軸指向春分點方向,OZ軸垂直于赤道面,并指向北極方向,OY軸與OX和OZ軸成右手系。

    (2)軌道坐標系:原點o位于主星質(zhì)心,ox軸由地心指向衛(wèi)星方向,oz軸垂直于軌道面,并指向軌道角動量方向,oy軸和ox、oz軸構成右手系。

    圖1 坐標系Fig.1 Coordinate frame

    2.2 衛(wèi)星編隊重構可行性分析

    假設編隊衛(wèi)星主星在圓軌道上運行,在主星軌道坐標系下,各顆從星其相對于主星的位置變化規(guī)律可用CW方程表述。令x,y,z為主星軌道坐標系的坐標分量,為主星軌道坐標系內(nèi)表達的速度分量??紤]進行連續(xù)推力軌道控制,令三個方向的控制加速度分別用Tx,Ty,Tz表示,它們可以是關于時間的函數(shù),則相應的運動方程可以寫成

    式中:ω為軌道坐標系相對慣性坐標系的轉(zhuǎn)動角速度。

    在衛(wèi)星編隊飛行重構的過程中,總的燃料消耗可以表示為各個方向控制加速度的函數(shù),假設編號為i的編隊成員衛(wèi)星初始相對運動狀態(tài)和期望相對運動狀態(tài)分別表述為ri0,rif,那么在編隊重構過程中,燃料消耗表示為兩個狀態(tài)量的函數(shù)。在滿足編隊飛行總體構型性能的前提下,根據(jù)每一個編隊成員初始狀態(tài),合理確定其期望狀態(tài),以達到每顆衛(wèi)星的剩余燃料差別最小,實現(xiàn)燃料平衡的目的。

    2.3 核心技術分析

    1)連續(xù)推力模型

    在對編隊重構中的推力器工作情況未知的情況下,可以將Tx,Ty,Tz寫成關于時間的任意連續(xù)函數(shù),本文將其寫成傅里葉函數(shù)形式[13]:

    式中:ax0,ay0,az0,bxn,byn,bzn為常值量;ωf=2π/tf。為了簡化運算,只取前5項的表達式作為加速度的表達形式。

    2)編隊重構定義

    編隊重構涉及到兩個方面:x-y軌道平面內(nèi)機動和軌道面外機動。給定軌道機動時間tf,通過連續(xù)推力控制,衛(wèi)星從初始位置和速度,到達目標狀態(tài)。在主星軌道坐標系下,o點與從星在y-z平面的投影位置的連線,按逆時針方向與z軸的夾角α稱為相位角。在主星軌道周期內(nèi),從星在y-z平面投影位置與主星位置連線掃過的圓的半徑定義為投影圓半徑。圖2給出了投影圓編隊衛(wèi)星的相位角關系。

    圖2 相位角Fig.2 Phase angle

    衛(wèi)星重構包含以下兩個方面的內(nèi)容:①從初始編隊投影圓半徑R,經(jīng)過軌道重構后,編隊重構后新的投影圓半徑為R+ΔR;②從初始相位角α,經(jīng)過軌道重構后,相位角變?yōu)棣粒う痢?/p>

    3)最優(yōu)控制技術

    最優(yōu)控制方法從每個從星的軌道重構出發(fā),分別建立各自的燃料目標函數(shù),通過對每個從星燃料消耗的合理調(diào)度,達到機動后編隊衛(wèi)星之間彼此燃料平衡的要求,并限制總?cè)剂舷?,從而在完成衛(wèi)星編隊燃料平衡問題的同時,使得編隊重構總?cè)剂舷淖顑?yōu)。

    3 編隊衛(wèi)星重構方法

    3.1 編隊衛(wèi)星重構問題研究

    在連續(xù)推力假設下,相對軌道面外和軌道面內(nèi)的燃料目標函數(shù)可以寫成

    對每一個編隊成員衛(wèi)星,在編隊重構過程中,最優(yōu)控制可以寫成(相關推導見參考文獻[13])

    式中:T0,T1,λ0,λ1是從星初始相位角、投影圓半徑及其變化函數(shù);Λ,Φ,ζ為推導過程中產(chǎn)生的中間變量(見參考文獻[13])。

    根據(jù)z方向和x-y方向的燃料消耗最優(yōu)目標函數(shù),經(jīng)過優(yōu)化目標函數(shù)J=Jxy+Jz,可以得到最優(yōu)相位角α*的兩個解:

    根據(jù)取得極小值條件,驗證哪一個是取得極小值的最優(yōu)相位角。這也證實了在這種機動策略下,總的燃料消耗僅僅取決于從星在主星軌道坐標系的相位角初值,因此通過選擇合理的從星初始相位角,可以達到節(jié)省燃料的目的。

    3.2 實例驗證

    假設主星運行在高度為1000km圓軌道,初始編隊投影圓半徑為30km,軌道重構后的軌道半徑為50km,相位角變化30°,轉(zhuǎn)移時間1600s。根據(jù)式(5)的結論,可以得到α*=72.352 6°時燃料最省。一般來說,使得燃料平衡的解并不是唯一的,總是趨向于選取燃料消耗最小的那一組解作為最優(yōu)解。

    圖3給出了編隊軌道重構過程中的相對軌道變化情況,細圓圈“o”表示無控情況下,從星經(jīng)過時間tf之后的位置,黑線表示在主星軌道坐標系下的從星相對轉(zhuǎn)移軌道。從圖3中可以看出,從星從初始相對軌道半徑為30km的投影圓軌道出發(fā),經(jīng)過轉(zhuǎn)移時間tf,重構后的相對軌道半徑為50km,能夠完成軌道半徑變化的要求。在軌道重構過程中,完成軌道重構需要的速度增量為0.056 8km/s。

    圖3 主星軌道坐標系下的軌道重構Fig.3 Orbital reconfiguration described in orbital frame of chief satellite

    假設從星自初始相對位置經(jīng)過,自由飛行時間tf后的相對位置為(yf0,zf0),采用連續(xù)推力飛行時間tf后的位置為(yf,zf),則有

    相應地,y-z平面投影圓相位角的變化為

    這與期望的相位變化角一致,在完成對相對軌道半徑大小控制的同時完成了相位角變化機動。

    圖4給出了軌道轉(zhuǎn)移過程中的燃料消耗隨初始相位角變化的情況,從圖4中可以看出,燃料消耗變化的周期是180°,使得燃料消耗最少的初始相位角為αmin=72.352 6°,這和計算得到的結果是一致的,從而證實了上述原理的正確性。由于解的周期性,可見只要相位角的最優(yōu)解并不是唯一的,這里只給出了0~180°的結果。

    圖4 燃料消耗-初始相位角變化曲線Fig.4 Variation of fuel consumption as a function of initial phase angle

    4 編隊燃料平衡方法

    兩個從星的軌道重構所對應的三星燃料平衡問題。假設兩顆從星的初始質(zhì)量和攜帶燃料的狀況相同,那么一個簡單的想法是在同等的轉(zhuǎn)移條件下(轉(zhuǎn)移時間,起始時間等),使得兩顆星的燃料消耗相等,并且能夠完成相應的軌道重構半徑和相位角的改變。

    4.1 三星燃料平衡問題

    本節(jié)討論的三星燃料平衡問題,是指在主從式衛(wèi)星編隊中,從星#1和從星#2同一時間從不同的初始位置進行軌道機動,達到一個新的編隊構型。經(jīng)過這樣的軌道重構后,從星#1和從星#2相對主星不僅相對距離發(fā)生變化,而且在主星軌道坐標系下,其相位角也發(fā)生改變,相位角的改變量分別用Δα1,Δα2表示,如圖5所示。對于具有兩顆從星的編隊衛(wèi)星重構問題,如果其初始燃料相同,并且在軌道重構過程各顆從星之間彼此燃料消耗差別最小,則滿足編隊燃料平衡的需求,其流程如圖6所示。

    引入如下中間變量

    式中:ω為軌道坐標系相對慣性坐標系的角速度;θ為從星間的初始相位差;Δα為從星相位變化量。

    假設兩顆星初始相位角差θ已知,因此,只要能夠求得其中一顆星的初始相位角,就可以確定另一顆衛(wèi)星的初始相位角,經(jīng)過理論推導和構建燃料平衡的目標函數(shù),滿足燃料平衡時,相位角須滿足:

    圖5 軌道重構示意圖Fig.5 Depiction of orbital reconfiguration

    圖6 編隊重構燃料平衡方法流程Fig.6 Flow chart of balancing formation reconfiguration

    4.2 實例仿真

    假設主星在軌道高度為1000km的圓軌道上運行,在主星軌道坐標系下,兩顆從星在y-z平面半徑為30km的投影圓軌道上運行,初始相位角相差60°。在t0時刻,兩顆從星開始機動,經(jīng)過時間tf,兩顆從星到達新的投影圓相對運行軌道,新的軌道y-z平面軌道半徑為50km,假設兩顆從星相位角變化相同,即Δα1=Δα2=30°,對4.1節(jié)中提出的燃料平衡算法進行仿真驗證。

    基于上述假設,從星#2的初始相位角與從星#1的初始相位角之間恒定差一個角度θ,因此,兩顆從星的燃料消耗都可以看作從星#1初始相位角的函數(shù),如圖7所示,實線表示從星#1的燃料消耗-初始相位角變化曲線,虛線表示從星#2的燃料消耗-從星#1初始相位角變化曲線??梢宰C明,燃料-相位角變化曲線的周期為180°,因此,我們這里只關注初始相位角在-90°≤α0≤90°范圍內(nèi)的結果,從圖7中可以看出,兩條曲線共有兩個交點,分別為,并且在從星#1的初始相位角時,達到燃料平衡,并且燃料消耗最優(yōu)。在軌道重構過程中,完成軌道重構需要的速度增量為0.064 8km/s。

    圖7 兩顆從星的燃料消耗隨從星#1初始相位角變化曲線Fig.7 Variations of fuel consumptions of two deputy satellites with different initial phase angles of deputy#1

    通過施加連續(xù)推力,從星#1和從星#2從初始位置出發(fā),進入轉(zhuǎn)移軌道,脫離原來半徑為30km的投影圓軌道,經(jīng)過設定的軌道重構時間tf,同時到達預定的半徑為50km的投影圓軌道,整個過程在y-z平面的投影如圖8所示,其中粗線圈表示從星#1和#2在t0時刻的位置以及在t0+tf時刻的位置,從圖8中可以看出采用4.1節(jié)給出的方法能夠完成軌道重構,滿足相對軌道半徑變化的要求。與此同時,軌道重構另一個目標要求是完成相位角的改變,初始時刻兩顆星的相位角差為θ=60°。對于從星#1,機動前后的相位變化角為Δα,同樣地,對于從星#2有一樣的結論,經(jīng)計算,Δα=30°,這與任務要求的兩顆星相位角變化30°的要求是一致的。因此,此方法既可以完成相對軌道位置的變化,又可以完成相位角變化的要求,從而證明了該方法的有效性。

    圖9給出了3個主星軌道坐標系下三個坐標軸方向的變化曲線,從推力函數(shù)的構造不難看出,推力函數(shù)是有三角函數(shù)構成的初等函數(shù),因此表現(xiàn)為不同峰值、不同頻率的周期變化曲線的疊加形式,與PD控制率等不同,其連續(xù)變化曲線在逐漸接近目標位置時,并不表現(xiàn)為逐漸趨近于0,因此,在本文的理論框架下,在設定好轉(zhuǎn)移時間后,須要主動關閉發(fā)動機。

    圖8 y-z平面投影圓相位角變化圖Fig.8 An increment in phase angle illustrated in projected circle in y-z plane

    圖9 從星#2推力變化情況Fig.9 Working status of thrusts of deputy#2

    5 結論

    本文對衛(wèi)星編隊重構過程的燃料平衡問題進行了研究,在給定初始軌道情況及目標軌道參數(shù)的情況下,給出了衛(wèi)星編隊重構方法,在此基礎上,對編隊重構中的燃料平衡問題進行了研究。經(jīng)研究,得出了以下結論:

    (1)在給定初始編隊狀態(tài)和目標編隊狀態(tài)的情況下,可以通過連續(xù)推力控制實現(xiàn)軌道重構。本文給出了編隊重構方法,并對其有效性進行了仿真驗證。

    (2)在初始燃料相同的情況下,研究了雙從星編隊重構的燃料平衡問題,通過構造燃料目標函數(shù),給出了求解最優(yōu)初始相位角的求解方法,在設定相對運動軌道半徑和相位角變化的情況下,實例仿真證明該方法能夠有效地完成機動過程,并達到燃料平衡。

    (References)

    [1]郝繼剛,張育林.SAR干涉測高分布式小衛(wèi)星編隊構型優(yōu)化設計[J].宇航學報,2006,27(4):654-669 Hao Jigang,Zhang Yulin.Formation optimized design for the height measurement of InSAR using distributed micro-satellites[J].Journal of Astronautics,2006,27(4):654-669(in Chinese)

    [2]Alberto Moreira,Gerhard Krieger.Spaceborne synthetic aperture radar(SAR)systems:state of the art and future developments[C]//11thGAAS Symposium.Munich:Global Automotive Aftermarket Symposium,2003:385-388

    [3]Jianying Wang,Zhaowei Sun.6-DOF robust adaptive terminal sliding mode control for spacecraft formation flying[J].Acta Astronautica,2012,73:76-87

    [4]Wu Y H,Cao X B,Xing Y J,et al.Relative motion coupled control for formation flying spacecraft via convex optimization[J].Aerospace Science and Technology,2010,14(6):415-428

    [5]Sung Moon Yoo,Sang Young Park,Kyu Hong Choi.A fuel balancing method for reconfiguration of satellite formation flying[C]//International Conference on Control,Automation and Systems.Seoul:International Mechanical and Mechatronics Engineering Committee,2007:1426-1429

    [6]W C Durham.Constrained control allocation[J].Journal of Guidance,Control,and Dynamics,1993,16(4):717-725

    [7]Johansen T A,F(xiàn)ossen T I.Control allocation—a survey[J].Automatica,2013,49(5):1087-1103

    [8]Shengyong T,Zhang S,Zhang Y.A modified direct allocation algorithm with application to redundant actuators[J].Chinese Journal of Aeronautics,2011,24(3):299-308

    [9]Susan A Frost,Marc Bodson.Resource balancing control allocation[C]//American Control Conference,Washington D.C.:AIAA,2010:1326-1331

    [10]Bodson M,F(xiàn)rost S A.Load balancing in control allocation[J].Journal of Guidance,Control,and Dynamics,2011,34(2):380-387

    [11]Rahmani A,Mesbahi M,Hadaegh F Y.Optimal balancedenergy formation flying maneuvers[J].Journal of guidance,control,and dynamics,2006,29(6):1395-1403

    [12]Beard R W,McLain T W,Hadaegh F Y.Fuel optimization for constrained rotation of spacecraft formations[J].Journal of Guidance,Control,and Dynamics,2000,23(2):339-346

    [13]P Palmer.Optimal relocation of satellites flying in nearcircular-orbit formations[J].Journal of Guidance,Control,and Dynamics,2006,29(3):519-526

    (編輯:張小琳)

    Fuel-balancing Using Continuous Thrust in Satellite Formation Flying

    DUAN Chenyang ZHAO Yafei ZHANG Shijie KONG Xianren
    (Research Center of Satellite Technology,Harbin Institute of Technology,Harbin 150080,China)

    In satellite formation flying,maneuvering such as formation reconfiguration would cause unbalanced fuel distribution among the formation satellites.It may be worse if one of the satellites fails,which leads to a short lifetime formation,or complete loss of its basic function.Aiming at the problem the paper proposes a fuel-balance approach using continuous thrust,which can be adopted during orbital reconfiguration of satellite formation.The approach uses continuous thrust,brings up a fuel-consumption objective function,and derives different optimal control acceleration with different phase angles and radius changes.We reduce the fuel difference in the formation by minimizing the difference of fuel-consumption functions between each two satellites.Fuel-balancing depends on the different initial phase angle.This paper gives a expression of initial phase angle.Eventually,we use a formation of three satellites in simulation,the simulation proves the proposed methods,continuous thrust control and fuel-balancing method are valid and effective.

    satellite formation;orbit reconfiguration;fuel-balance;continuous thrust

    V474.3

    :ADOI:10.3969/j.issn.1673-8748.2016.01.003

    2015-11-19;

    :2016-01-09

    國家自然科學基金(91438202),民用航天“十二五”預研項目

    段晨陽,男,博士研究生,研究方向為航天器動力學與控制。Email:hitsat@gmail.com。

    猜你喜歡
    主星相位角編隊
    2023年1月25日,美軍一次演習期間,空軍正在進行編隊飛行
    軍事文摘(2023年5期)2023-03-27 08:56:26
    開封地區(qū)健康成人相位角及其影響因素
    相位角對容性耦合電非對稱放電特性的影響?
    物理學報(2018年22期)2018-12-18 05:58:24
    歐洲嘗試用魚叉和網(wǎng)來清理太空垃圾
    首顆繞褐矮星旋轉(zhuǎn)的巨行星現(xiàn)身
    科學導報(2018年22期)2018-05-14 12:05:58
    基于事件驅(qū)動的多飛行器編隊協(xié)同控制
    共軸剛性旋翼直升機旋翼控制相位角問題分析
    輻射開環(huán)繩系衛(wèi)星編隊自旋展開動力學與控制策略
    航空學報(2016年9期)2016-12-06 07:07:51
    精確檢測電網(wǎng)電壓矢量相位角方法的研究
    基于蟻群算法的分布式衛(wèi)星光網(wǎng)絡波長路由分配技術研究
    欧美+日韩+精品| 一区福利在线观看| 国产精品电影一区二区三区| 好男人在线观看高清免费视频| 一区福利在线观看| 日韩中字成人| 亚洲精品日韩在线中文字幕 | 欧美最黄视频在线播放免费| 天堂影院成人在线观看| 91精品国产九色| 桃色一区二区三区在线观看| 非洲黑人性xxxx精品又粗又长| 精品一区二区三区视频在线| 高清毛片免费观看视频网站| 男女之事视频高清在线观看| 国产精品精品国产色婷婷| 日本免费一区二区三区高清不卡| 少妇高潮的动态图| 中文资源天堂在线| 波多野结衣巨乳人妻| av在线老鸭窝| 又黄又爽又刺激的免费视频.| 国产高清三级在线| 日本一本二区三区精品| 国产又黄又爽又无遮挡在线| 欧美成人精品欧美一级黄| 久久精品国产亚洲网站| 在线免费十八禁| 国产成人影院久久av| 国产亚洲av嫩草精品影院| 欧美xxxx黑人xx丫x性爽| 国产精品久久久久久av不卡| 亚洲婷婷狠狠爱综合网| 国产高清有码在线观看视频| 国产黄色小视频在线观看| 久久精品综合一区二区三区| 综合色av麻豆| 精品免费久久久久久久清纯| 欧美性猛交╳xxx乱大交人| 麻豆国产97在线/欧美| a级毛片免费高清观看在线播放| 深夜精品福利| 国产精品人妻久久久久久| 亚洲国产精品久久男人天堂| 国产久久久一区二区三区| 亚洲图色成人| 男女那种视频在线观看| 亚洲不卡免费看| 久久久久精品国产欧美久久久| 国产精品一区二区免费欧美| 别揉我奶头 嗯啊视频| 亚洲久久久久久中文字幕| 亚洲欧美精品综合久久99| 亚洲成a人片在线一区二区| 18禁在线无遮挡免费观看视频 | 两性午夜刺激爽爽歪歪视频在线观看| 国产精品电影一区二区三区| av专区在线播放| 国产精品野战在线观看| 少妇的逼水好多| 国产午夜精品论理片| 尾随美女入室| 97超碰精品成人国产| 在线国产一区二区在线| 国产亚洲精品av在线| 内射极品少妇av片p| 少妇熟女aⅴ在线视频| 色综合色国产| 婷婷精品国产亚洲av| 欧美日韩一区二区视频在线观看视频在线 | 97热精品久久久久久| 国产成人a∨麻豆精品| 美女内射精品一级片tv| 国模一区二区三区四区视频| 性色avwww在线观看| 成人无遮挡网站| 毛片女人毛片| 国产精品嫩草影院av在线观看| 久久午夜亚洲精品久久| 给我免费播放毛片高清在线观看| videossex国产| 美女cb高潮喷水在线观看| 51国产日韩欧美| 在线看三级毛片| 人妻久久中文字幕网| 久久久久精品国产欧美久久久| 午夜精品一区二区三区免费看| 美女被艹到高潮喷水动态| 亚洲成av人片在线播放无| 国产精品国产高清国产av| 欧美日韩综合久久久久久| 成年女人看的毛片在线观看| 精品人妻偷拍中文字幕| 无遮挡黄片免费观看| 男女下面进入的视频免费午夜| 成人国产麻豆网| 免费av毛片视频| 国产精品一区二区三区四区久久| 久久久久久大精品| 亚洲七黄色美女视频| 久久韩国三级中文字幕| 亚洲婷婷狠狠爱综合网| 99久久精品国产国产毛片| 日韩 亚洲 欧美在线| 成人二区视频| 欧美成人免费av一区二区三区| 久久久久精品国产欧美久久久| 97超碰精品成人国产| 成人永久免费在线观看视频| 在线天堂最新版资源| 中文字幕av成人在线电影| 久久欧美精品欧美久久欧美| 亚洲美女视频黄频| 精品国产三级普通话版| 国产一区二区在线观看日韩| 精品乱码久久久久久99久播| 长腿黑丝高跟| 观看美女的网站| 亚洲婷婷狠狠爱综合网| av女优亚洲男人天堂| 51国产日韩欧美| 亚洲一区高清亚洲精品| 最近最新中文字幕大全电影3| 亚洲精华国产精华液的使用体验 | 真人做人爱边吃奶动态| 日韩精品中文字幕看吧| 欧美日韩一区二区视频在线观看视频在线 | 亚洲精品一卡2卡三卡4卡5卡| ponron亚洲| 亚洲精品日韩av片在线观看| 亚洲美女视频黄频| 99在线人妻在线中文字幕| 精品免费久久久久久久清纯| 99riav亚洲国产免费| 99热6这里只有精品| 日本黄色片子视频| 十八禁国产超污无遮挡网站| 欧美日韩国产亚洲二区| 国产伦在线观看视频一区| 不卡视频在线观看欧美| 亚洲中文字幕日韩| 嫩草影视91久久| 日本与韩国留学比较| 午夜a级毛片| 日日撸夜夜添| 97在线视频观看| 国产成人a∨麻豆精品| 一本一本综合久久| 欧美色欧美亚洲另类二区| 中文字幕免费在线视频6| 国产精品一区二区免费欧美| 日韩大尺度精品在线看网址| 18+在线观看网站| 亚洲,欧美,日韩| 全区人妻精品视频| 村上凉子中文字幕在线| 五月玫瑰六月丁香| 国国产精品蜜臀av免费| 97在线视频观看| 欧美xxxx黑人xx丫x性爽| 国产单亲对白刺激| 国产在视频线在精品| 又爽又黄a免费视频| 欧美日韩国产亚洲二区| 亚洲中文字幕日韩| 人人妻人人澡欧美一区二区| 国产精品人妻久久久影院| 亚洲av一区综合| 国产高清有码在线观看视频| 亚洲欧美日韩高清专用| 91精品国产九色| 成人一区二区视频在线观看| 日韩国内少妇激情av| 国产乱人视频| 精品国产三级普通话版| 深爱激情五月婷婷| 日韩 亚洲 欧美在线| 丝袜喷水一区| 国产 一区 欧美 日韩| 国产亚洲精品av在线| 小说图片视频综合网站| 蜜臀久久99精品久久宅男| 日韩精品有码人妻一区| 国产精品一及| 欧美一区二区国产精品久久精品| а√天堂www在线а√下载| 老司机影院成人| 国产单亲对白刺激| av在线老鸭窝| 亚洲国产日韩欧美精品在线观看| 国产毛片a区久久久久| 神马国产精品三级电影在线观看| 国产精品av视频在线免费观看| 91麻豆精品激情在线观看国产| 天美传媒精品一区二区| 午夜精品在线福利| 特大巨黑吊av在线直播| 日日摸夜夜添夜夜添小说| 亚洲经典国产精华液单| 午夜福利在线在线| 最新中文字幕久久久久| 日韩强制内射视频| 性色avwww在线观看| 成年女人毛片免费观看观看9| 99热精品在线国产| 国产精品一区二区三区四区免费观看 | 小说图片视频综合网站| 免费观看在线日韩| 亚洲成人中文字幕在线播放| 有码 亚洲区| 97热精品久久久久久| 又爽又黄无遮挡网站| 国内精品久久久久精免费| 国产伦在线观看视频一区| 成人毛片a级毛片在线播放| 少妇被粗大猛烈的视频| 九九热线精品视视频播放| 男人的好看免费观看在线视频| 国产aⅴ精品一区二区三区波| 国产精品电影一区二区三区| 三级经典国产精品| 亚洲真实伦在线观看| or卡值多少钱| 亚洲欧美中文字幕日韩二区| 真实男女啪啪啪动态图| 欧美激情久久久久久爽电影| 2021天堂中文幕一二区在线观| 床上黄色一级片| 免费人成视频x8x8入口观看| 熟女人妻精品中文字幕| 国产在线精品亚洲第一网站| 国产黄片美女视频| 亚洲精品在线观看二区| 日韩精品有码人妻一区| 国产一区二区亚洲精品在线观看| 国产亚洲av嫩草精品影院| 欧美另类亚洲清纯唯美| 日韩在线高清观看一区二区三区| av在线蜜桃| 国产私拍福利视频在线观看| 我要搜黄色片| 美女 人体艺术 gogo| 人妻丰满熟妇av一区二区三区| 人人妻人人澡人人爽人人夜夜 | 国产成人a区在线观看| 国产精品国产高清国产av| 深夜a级毛片| 久久久久久久午夜电影| 成人欧美大片| 国产成人91sexporn| 国模一区二区三区四区视频| 久久久精品94久久精品| 成人无遮挡网站| 亚洲高清免费不卡视频| 韩国av在线不卡| 人人妻人人看人人澡| 12—13女人毛片做爰片一| 搡女人真爽免费视频火全软件 | 精品一区二区三区视频在线| 亚洲天堂国产精品一区在线| 在线观看66精品国产| www日本黄色视频网| av在线天堂中文字幕| 99热6这里只有精品| 又爽又黄a免费视频| 欧美区成人在线视频| 久久亚洲精品不卡| 两个人视频免费观看高清| 夜夜夜夜夜久久久久| 免费av观看视频| 日韩欧美 国产精品| 在线免费观看的www视频| 国产精品久久视频播放| 亚洲欧美精品综合久久99| 日韩国内少妇激情av| 永久网站在线| 国内揄拍国产精品人妻在线| 成人漫画全彩无遮挡| 成人性生交大片免费视频hd| 国产aⅴ精品一区二区三区波| 国语自产精品视频在线第100页| 国产伦精品一区二区三区视频9| 岛国在线免费视频观看| 又爽又黄无遮挡网站| 亚洲最大成人中文| 一级a爱片免费观看的视频| 日本免费a在线| 欧美中文日本在线观看视频| 小蜜桃在线观看免费完整版高清| 国产爱豆传媒在线观看| 亚洲无线在线观看| 热99re8久久精品国产| 成人一区二区视频在线观看| 国内少妇人妻偷人精品xxx网站| 国产国拍精品亚洲av在线观看| 亚洲欧美日韩高清专用| 观看免费一级毛片| 在线天堂最新版资源| 久久国产乱子免费精品| 色哟哟·www| 国产视频内射| 亚洲最大成人中文| 看非洲黑人一级黄片| or卡值多少钱| 欧美xxxx性猛交bbbb| 日韩欧美在线乱码| 国产 一区精品| 91在线精品国自产拍蜜月| 欧美另类亚洲清纯唯美| 国产白丝娇喘喷水9色精品| 99热网站在线观看| 欧美激情在线99| 男人狂女人下面高潮的视频| 赤兔流量卡办理| 免费在线观看影片大全网站| 99在线人妻在线中文字幕| 国产精品一区www在线观看| 国产又黄又爽又无遮挡在线| 亚洲成人久久性| 成年女人永久免费观看视频| 精品欧美国产一区二区三| 国内精品久久久久精免费| 日本五十路高清| 性插视频无遮挡在线免费观看| 欧美国产日韩亚洲一区| 精品不卡国产一区二区三区| 免费av毛片视频| 日韩高清综合在线| 久久精品国产鲁丝片午夜精品| 亚洲中文字幕日韩| 国产国拍精品亚洲av在线观看| 国产黄色视频一区二区在线观看 | 91在线精品国自产拍蜜月| 亚洲不卡免费看| 精品熟女少妇av免费看| or卡值多少钱| 99热全是精品| 小蜜桃在线观看免费完整版高清| 午夜福利在线在线| 久久综合国产亚洲精品| 成年女人毛片免费观看观看9| 黄色一级大片看看| АⅤ资源中文在线天堂| 精品无人区乱码1区二区| 国产精品亚洲美女久久久| 欧美bdsm另类| 一本久久中文字幕| 波多野结衣巨乳人妻| 久久精品久久久久久噜噜老黄 | 变态另类成人亚洲欧美熟女| 亚洲国产欧洲综合997久久,| 秋霞在线观看毛片| 欧美绝顶高潮抽搐喷水| 男人舔女人下体高潮全视频| 久久久久国内视频| 国产精品电影一区二区三区| 久久精品国产自在天天线| 毛片女人毛片| 亚洲中文日韩欧美视频| 男人舔女人下体高潮全视频| 一本久久中文字幕| 一级毛片aaaaaa免费看小| 日产精品乱码卡一卡2卡三| 免费大片18禁| av黄色大香蕉| 91午夜精品亚洲一区二区三区| 91久久精品国产一区二区成人| 亚洲av成人精品一区久久| 男人舔奶头视频| 91久久精品国产一区二区成人| 国内久久婷婷六月综合欲色啪| 国产一级毛片七仙女欲春2| 亚洲国产欧美人成| 国产不卡一卡二| 国产淫片久久久久久久久| 亚洲,欧美,日韩| 长腿黑丝高跟| 欧美丝袜亚洲另类| 日韩国内少妇激情av| 亚洲综合色惰| 久久久久国内视频| 日韩一区二区视频免费看| av国产免费在线观看| 男女视频在线观看网站免费| 丝袜喷水一区| 欧美又色又爽又黄视频| 夜夜爽天天搞| 精品日产1卡2卡| 亚洲电影在线观看av| 亚洲国产欧洲综合997久久,| 99热6这里只有精品| 伊人久久精品亚洲午夜| 久久久久久久久大av| 精品一区二区三区视频在线| 欧美人与善性xxx| 久久精品国产清高在天天线| 精品无人区乱码1区二区| 白带黄色成豆腐渣| 亚洲熟妇熟女久久| 男人狂女人下面高潮的视频| 少妇熟女aⅴ在线视频| 丰满的人妻完整版| 日本撒尿小便嘘嘘汇集6| 亚洲第一区二区三区不卡| 舔av片在线| 丝袜喷水一区| 婷婷精品国产亚洲av在线| 午夜福利在线观看吧| 亚洲国产日韩欧美精品在线观看| 久久久久久伊人网av| 国产黄色视频一区二区在线观看 | 欧美日韩乱码在线| 天堂影院成人在线观看| 两个人的视频大全免费| 乱人视频在线观看| 久久久久久久久久久丰满| 国产av在哪里看| 最后的刺客免费高清国语| 国产亚洲精品av在线| 精品乱码久久久久久99久播| 99热全是精品| 九九在线视频观看精品| 久久久久久大精品| 99久久成人亚洲精品观看| 亚洲精品亚洲一区二区| 最后的刺客免费高清国语| 搞女人的毛片| 日韩三级伦理在线观看| 国产91av在线免费观看| 悠悠久久av| 嫩草影院入口| 欧美高清成人免费视频www| 精品午夜福利在线看| 成年av动漫网址| 夜夜看夜夜爽夜夜摸| 我要搜黄色片| 国产伦精品一区二区三区四那| 搡女人真爽免费视频火全软件 | 国产精品电影一区二区三区| 男女做爰动态图高潮gif福利片| 欧美色欧美亚洲另类二区| 欧美一区二区精品小视频在线| 欧美激情在线99| 久久久欧美国产精品| 久久中文看片网| 精品午夜福利视频在线观看一区| 男女啪啪激烈高潮av片| 亚洲精品国产av成人精品 | 我要搜黄色片| 国产精品野战在线观看| 香蕉av资源在线| 久久草成人影院| 日韩欧美 国产精品| 最近手机中文字幕大全| 国产伦一二天堂av在线观看| 国产成年人精品一区二区| 欧美极品一区二区三区四区| 亚洲一区二区三区色噜噜| 99久久无色码亚洲精品果冻| 国产在线男女| 深夜精品福利| 欧美另类亚洲清纯唯美| 亚洲激情五月婷婷啪啪| 老司机影院成人| 少妇高潮的动态图| 听说在线观看完整版免费高清| 男女之事视频高清在线观看| 变态另类丝袜制服| 人妻久久中文字幕网| 国产三级在线视频| 国产人妻一区二区三区在| 99视频精品全部免费 在线| 一级毛片久久久久久久久女| 99热精品在线国产| 丰满的人妻完整版| 成年免费大片在线观看| 国产精品综合久久久久久久免费| 一本久久中文字幕| 一个人看视频在线观看www免费| 激情 狠狠 欧美| 啦啦啦观看免费观看视频高清| 日韩欧美一区二区三区在线观看| a级毛色黄片| 日本免费一区二区三区高清不卡| 最近手机中文字幕大全| 亚洲经典国产精华液单| 亚洲七黄色美女视频| a级毛色黄片| 老熟妇仑乱视频hdxx| 春色校园在线视频观看| 国产精品国产高清国产av| 欧美一区二区亚洲| 真人做人爱边吃奶动态| 亚洲精品粉嫩美女一区| 真实男女啪啪啪动态图| 少妇被粗大猛烈的视频| 女人被狂操c到高潮| 人人妻人人澡欧美一区二区| 成人性生交大片免费视频hd| 亚洲中文字幕日韩| 嫩草影院入口| 久久久久国产精品人妻aⅴ院| 搡老妇女老女人老熟妇| 真人做人爱边吃奶动态| 长腿黑丝高跟| 淫妇啪啪啪对白视频| 亚洲av免费在线观看| 日本在线视频免费播放| 欧美日韩精品成人综合77777| a级毛色黄片| 在线观看av片永久免费下载| 午夜精品在线福利| 国产熟女欧美一区二区| 亚洲av第一区精品v没综合| 亚洲中文日韩欧美视频| 久久久久久大精品| 99久久精品一区二区三区| 99在线人妻在线中文字幕| 久久久a久久爽久久v久久| 色吧在线观看| 乱码一卡2卡4卡精品| 搡女人真爽免费视频火全软件 | 久久久久国产网址| 日本免费a在线| av在线亚洲专区| 黄色视频,在线免费观看| 午夜免费男女啪啪视频观看 | 国产探花在线观看一区二区| 亚洲av免费在线观看| 最近最新中文字幕大全电影3| 日韩欧美一区二区三区在线观看| 国产精品永久免费网站| 看非洲黑人一级黄片| 蜜桃久久精品国产亚洲av| 国产一区二区在线观看日韩| 国产精品久久久久久久久免| 国产精品久久久久久亚洲av鲁大| 国产真实乱freesex| 成人毛片a级毛片在线播放| 国产片特级美女逼逼视频| 国内精品美女久久久久久| 午夜精品在线福利| 不卡一级毛片| 久久久国产成人免费| 久久热精品热| 久久精品国产亚洲网站| 久久久午夜欧美精品| 国产亚洲欧美98| 一级a爱片免费观看的视频| 男人狂女人下面高潮的视频| 亚洲国产精品sss在线观看| 极品教师在线视频| 久久精品国产自在天天线| 人妻久久中文字幕网| 成人国产麻豆网| 久久精品国产鲁丝片午夜精品| 久久久久久久久大av| 老师上课跳d突然被开到最大视频| 国产成人91sexporn| 午夜久久久久精精品| 五月玫瑰六月丁香| 可以在线观看毛片的网站| 男女边吃奶边做爰视频| 人人妻人人看人人澡| 听说在线观看完整版免费高清| 精品少妇黑人巨大在线播放 | ponron亚洲| 成人性生交大片免费视频hd| av天堂在线播放| 亚洲国产色片| 亚洲真实伦在线观看| 女人被狂操c到高潮| 91av网一区二区| 91精品国产九色| 99久久精品热视频| 老女人水多毛片| 五月玫瑰六月丁香| 真人做人爱边吃奶动态| 最近视频中文字幕2019在线8| 国产成人一区二区在线| 综合色丁香网| 又黄又爽又免费观看的视频| 亚洲精品久久国产高清桃花| 午夜福利在线观看免费完整高清在 | 搡老岳熟女国产| 干丝袜人妻中文字幕| 国产一区二区在线观看日韩| 日本三级黄在线观看| 久久婷婷人人爽人人干人人爱| 在线播放国产精品三级| 国产精品,欧美在线| 亚洲丝袜综合中文字幕| 成人特级av手机在线观看| 日韩成人伦理影院| 久久6这里有精品| 高清午夜精品一区二区三区 | 日本撒尿小便嘘嘘汇集6| 日韩欧美精品v在线| 欧美国产日韩亚洲一区| 免费不卡的大黄色大毛片视频在线观看 | 尤物成人国产欧美一区二区三区| 亚洲av美国av| 在线国产一区二区在线| 日韩成人av中文字幕在线观看 | 91在线精品国自产拍蜜月| 欧美最新免费一区二区三区| 中国美女看黄片| 99久久久亚洲精品蜜臀av| 国产中年淑女户外野战色| 国产 一区精品| 国产av一区在线观看免费| 十八禁网站免费在线| 在线国产一区二区在线| 观看美女的网站| 变态另类成人亚洲欧美熟女|