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

    空間軌道交會非線性優(yōu)化方法研究

    2018-03-23 01:59:27周新耀陶勇鵬
    上海航天 2018年1期
    關(guān)鍵詞:交會航天器增量

    周新耀,陶勇鵬,曾 亮,周 亮

    (上海機電工程研究所,上海 201109)

    0 引言

    航天技術(shù)的發(fā)展,以人造衛(wèi)星進入太空為起始點,迄今50年來取得長足的進步。在第一顆人造衛(wèi)星成功發(fā)射后,各國就發(fā)現(xiàn)它在軍事中的巨大潛力,早在20世紀60年代美國總統(tǒng)肯尼迪就說過:誰掌握了太空,誰就能更有效地控制地球。美、俄等超級大國更是不遺余力地進行著空間飛行器的研究,不斷搶占太空資源,力圖掌握絕對的制天權(quán)。

    空間資源爭奪的形勢日趨嚴峻,必然促進各種空間飛行器的發(fā)展,隨著空間格局的變化,空間交會越來越受到各國的重視。交會過程燃料消耗的多少很大程度上決定了交會任務(wù)是否能夠圓滿完成,因此交會軌道的優(yōu)化和控制具有重要的意義[1]。

    目前沖量交會基本都是基于Lambert問題[1-2]和基于微分修正法[3]的單脈沖交會,而Lambert問題是以中心引力體為模型,沒有考慮引力攝動因素,存在一定的誤差。單脈沖交會不存在優(yōu)化問題,但進行遠距離交會時消耗燃料太多。因此采用優(yōu)化后的多脈沖交會,可以有效地節(jié)省燃料和時間[4-11]。

    本文在考慮地球扁率J2攝動影響的情況下,采用非線性規(guī)劃方法對于遠距離情況求解多沖量最優(yōu)空間交會問題,避開了求解兩點邊值問題的復(fù)雜數(shù)學(xué)計算,很方便地得到多沖量交會的最優(yōu)解。在建立普遍適用的空間飛行器數(shù)學(xué)模型的基礎(chǔ)上,通過迭代算法估算交會時間從而確定預(yù)定交會點,進而按照時間、能量以及末端約束條件,建立了用于求解多脈沖最優(yōu)交會的非線性規(guī)劃模型;實現(xiàn)對雙脈沖固定時間交會軌道進行優(yōu)化設(shè)計,進而得到最優(yōu)速度增量的大小、方向和作用時間,但實際過程中,飛行器推力系統(tǒng)并不能提供一個瞬時速度增量,因此采用速度增益導(dǎo)引方法,實現(xiàn)交會軌道優(yōu)化方案[12-14]。最后通過數(shù)學(xué)仿真試驗證明了該方法的有效性。

    1 數(shù)學(xué)模型的建立

    定義地心赤道慣性坐標系(O-xiyizi),坐標系原點O在地球中心,平面xiyi和地球赤道平面重合,xi軸指向春分點,zi軸沿著地球旋轉(zhuǎn)軸(即垂直于赤道平面),指向北極。yi軸與xi軸、zi軸構(gòu)成右手直角坐標系[9]。

    在地心赤道慣性坐標系中,考慮地球扁平J2攝動作用的情況下,航天器運行軌道方程為

    (1)

    式中:r為航天器地心距矢量;r為r的模;v為航天器速度矢量;f為地球J2攝動加速度;μ=3.986×105km3/s2為地球引力常數(shù)。a為發(fā)動機推力的加速度矢量;f在地心慣性坐標系中的分量為

    (2)

    式中:x、y、z為航天器在地心慣性坐標系中的坐標;J2=1.082 64×10-3為地球引力勢的二階帶諧系數(shù);Re=6 378.14 km為地球半徑。

    (3)

    2 預(yù)測交會點確定

    為了達到交會的目的,必須先在目標航天器運行軌道上找到一個合適的預(yù)測交會點,預(yù)測交會點是有條件的,并不是目標將運行到的任何一點都能作為預(yù)測交會點。理論上,在運行時間確定后,空間的兩點間總能找到一條軌道。但預(yù)測交會點的選取必須合理,否則得到的變軌速度增量將超過實際許可,交會無法進行;或者交會軌道不合理,穿過地球表面等[4]。實際上,在交會航天器與目標航天器軌道為近圓的條件下,確定預(yù)測交會點的主要思想是:估算交會航天器和目標航天器運行到預(yù)測交會點的時間,如果兩者能夠大致相等,則此預(yù)測交會點是可取的。具體方法如圖1所示[5]。

    圖1 預(yù)測交會點確定Fig.1 Determination of predictive interception point

    已知一條空間軌道上的兩點r1,r2,則兩點之間的運行時間估算公式為

    (4)

    (5)

    初始條件給定了此刻的交會航天器和目標航天器的位置(分別為r1和r)和速度矢量。先利用數(shù)值方法對目標進行某個步長的軌道遞推,得到目標的一個步長后的位置矢量,設(shè)為r2,再利用剛才的方法分別估算交會航天器和目標航天器從初始位置運行到此位置矢量的大致時間,設(shè)為Δt1,Δt2。如此重復(fù),直到|1-Δt1/Δt2|<ε,便可得到預(yù)測交會點,它可保證比較容易地實現(xiàn)交會航天器與目標航天器同時到達此位置。當(dāng)然,在遞推過程中,必須對遞推步數(shù)作個限制,一般將目標航天器的遞推限制在一個周期內(nèi),如在此之內(nèi)|1-Δt1/Δt2|一直大于ε,則認為當(dāng)前的初始條件不能進行交會。

    3 遠距離交會軌道優(yōu)化

    3.1 最優(yōu)交會的非線性規(guī)劃模型

    首先對航天器多脈沖交會軌道進行非線性規(guī)劃,規(guī)劃過程中的約束條件如下:

    1) 運動狀態(tài)約束

    飛行器在脈沖施加點按照式(1)的規(guī)律運行。以x表示飛行器的狀態(tài)矢量(3個位置和3個速度矢量),令式(1)為Φ(x,t),則在脈沖施加點間(包括最后一次脈沖施加點和交會點間)飛行器的運動狀態(tài)約束為

    (6)

    2) 時間約束

    當(dāng)交會時間確定后,交會過程中的n次脈沖作用的時刻ti(i=1,2,…,n)應(yīng)在初始時刻t0和末端時刻tf之間。即滿足

    t0≤t1

    (7)

    3) 能量約束

    飛行器所攜帶的燃料有限,能夠提供的速度增量有一定的限制,設(shè)能提供的最大速度增量為Δvmax,則交會過程中施加的速度增量要滿足

    Δvtotal≤Δvmax

    (8)

    式中:Δvtotal為n次速度增量總和,表示為

    (9)

    4) 末端狀態(tài)約束

    設(shè)在交會點處飛行器和目標航天器的位置矢量分別為rIp和rtp,速度矢量分別為vIp和vtp。要保證飛行器最后要命中目標,需要滿足

    (10)

    5) 性能指標

    選取能量消耗最少作為性能指標,則性能指標為

    J=min(Δvmax)

    (11)

    綜合以上約束,非線性規(guī)劃模型可表示為

    (12)

    3.2 非線性規(guī)劃算法

    如果目標函數(shù)中或者約束條件中至少有一個是非線性函數(shù),則最優(yōu)化問題就叫做非線性規(guī)劃問題。它的一般形式為

    式中:X=(x1,x2,…,xn)∈Rn;f,gi,hj是定義在Rn上的實值函數(shù)。

    一般來說,解非線性規(guī)劃問題要比求解線性規(guī)劃問題困難得多,而且也不像線性規(guī)劃那樣有統(tǒng)一的數(shù)學(xué)模型及如單純形遺傳法這一通用解法[6]。非線性規(guī)劃的各種算法都有自己特定適用范圍。都有一定的局限性,到目前為止還沒有適合于各種非線性規(guī)劃問題的一般算法[13]。通常求解非線性規(guī)劃問題的方法有罰函數(shù)法和近似規(guī)劃法[8]。本文采用近似規(guī)劃法進行求解:

    式中:Δvmax為發(fā)動機能提供的最大速度增量,為一常值。h(X)=rIp-rtp。以雙脈沖為例X為第一次施加的速度增量和第二次施加的速度增量及其作用時間,為7個量。

    具體迭代算法如下:

    1) 給定初始可行點X1={(x1)1,(x2)1,…,(xn)1},步長限制為(δj)1(j=1,2,…,n),步長縮小系數(shù)β∈(0,1),取0.5。允許誤差ε>0,取1×10-6,令k=1。

    2) 在點Xk處,將f(X),g(X),h(X)按泰勒級數(shù)展開并取一階近似,得到近似線性規(guī)劃問題

    (14)

    3) 在上述近似線性規(guī)劃問題的基礎(chǔ)上增加一組限制步長的線性約束條件,因為線性近似通常只在展開點附近近似程度較高,故需要對變量的取值范圍加以限制,所以增量的約束條件為

    (15)

    求解該線性規(guī)劃問題,得到最優(yōu)解Xk+1。

    4) 檢驗Xk+1對原約束是否可行,若Xk+1對原約束可行,則轉(zhuǎn)化為步驟5);否則,縮小步長限制,令

    (δj)k=β(δj)k(j=1,2,…,n)

    (16)

    返回步驟3),重新解算當(dāng)前的線性規(guī)劃問題。

    5) 判斷精度:若(δj)k<ε(j=1,2,…,n),則點Xk+1為近似最優(yōu)解;否則令(δj)k+1=(δj)k(j=1,2,…,n),k=k+1,返回步驟2)。

    根據(jù)上述最優(yōu)軌道優(yōu)化算法,可以得出雙脈沖最優(yōu)點火時刻t1、t2以及該時刻所需要的速度增量Δv1、Δv2;但實際過程中,飛行器的推力系統(tǒng)并不能夠提供一個瞬時速度增量,而是通過速度增益導(dǎo)引方法來實現(xiàn)軌道轉(zhuǎn)移[10-12],如圖2所示。

    圖2 非線性最優(yōu)交會軌跡Fig.2 Optimal interception trajectory of nonlinearity

    4 仿真計算分析

    交會航天器包括平臺和載荷,其中載荷從平臺發(fā)射,自主變軌完成與目標的交會。設(shè)交會時間為12 000 s,載荷所能提供的最大速度增量為2.0 km/s。初始時刻平臺與目標航天器的參數(shù)如表1所示,目標星為中軌衛(wèi)星,軌道高度16 000 km,三個軌道面參數(shù)如表1所示。

    表1 軌道運行狀態(tài)

    由表2、表3和圖3、圖4對比,可以看出,采用Lambert方法所得的結(jié)果脫靶量較大,同時消耗燃料較多,而本文采取的非線性規(guī)劃方法,不僅速度增量減少了0.08 km/s,而且精度達到85.6 m,遠優(yōu)于Lambert的位置精度18 368 m,如圖5所示。

    表2 雙脈沖最優(yōu)交會結(jié)果

    表3 Lambert交會結(jié)果

    圖3 非線性規(guī)劃彈目相對距離曲線Fig.3 Relative distance curve of missile-target based on nonlinear programming

    圖4 Lambert彈目相對距離變化曲線Fig.4 Relative distance curve of missile-target based on Lambert

    圖5 Lambert交會軌跡Fig.5 Interception trajectory of Lambert

    5 結(jié)束語

    本文對空間軌道的遠距離交會進行了研究,通過迭代算法對預(yù)測交會點進行可行性預(yù)估,進而選取能量消耗和末端精度為控制變量,將軌道交會問題轉(zhuǎn)化為最優(yōu)控制問題,采用非線性規(guī)劃的方法對軌道離散化后所轉(zhuǎn)化的非線性問題進行求解,最終得到最優(yōu)交會軌道。相比于傳統(tǒng)的軌道優(yōu)化設(shè)計,本文將時間、能量以及末端位置、速度狀態(tài)變量同時引入約束中,在節(jié)省能量的前提下,保證了交會的速度需求,極大地提高了末端交會的控制精度。

    非線性規(guī)劃是解決最優(yōu)控制問題的一種有效方法,可以滿足控制中的各項約束,而且在求解過程中避免了復(fù)雜的計算[13]。通過仿真可以看出,本文所采用的非線性規(guī)劃策略大大提高了交會的控制精度,同時當(dāng)飛行時間較長時,單脈沖交會軌跡不再為最優(yōu)軌跡,在交會時間固定的前提下,非線性規(guī)劃得出的雙脈沖策略相比傳統(tǒng)Lambert軌道設(shè)計可有效降低速度增量,減小發(fā)動機燃料消耗。不過非線性規(guī)劃只是一種局部優(yōu)化方法,不能得到全局最優(yōu)解,并且隨著脈沖次數(shù)的增加,計算量較大,因此,如何在本文的研究基礎(chǔ)上,進一步對優(yōu)化算法進行改進,從而更加充分地利用有限燃料,提高對目標的交會精度將是今后需要研究的內(nèi)容。

    [1] 聞新, 楊嘉偉. 軍用衛(wèi)星的發(fā)展趨勢分析[J]. 現(xiàn)代防御技術(shù), 2002, 30(4): 7-11.

    [2] 張莉英, 張啟信, 王輝. 反衛(wèi)星武器技術(shù)及防御措施淺析[J]. 飛航導(dǎo)彈, 2004(3): 28-30.

    [3] 吳勤, 高雁翎. 美國空間對抗裝備的新進展[J]. 航天電子對抗, 2007, 23(3): 1-4.

    [4] 葛之江, 劉品雄, 王乃東. 美俄天基反衛(wèi)星技術(shù)的發(fā)展[J]. 航天電子對抗, 2003, 19(4): 1-5.

    [5] 馮志剛, 方昌華. 世界各國反衛(wèi)星策略綜述[J]. 中國航天, 2006(3): 38-40.

    [6] 周訨元, 劉宇, 譚天樂, 等. 基于遺傳算法的小推力星際轉(zhuǎn)移軌道設(shè)計與控制方法研究[J]. 上海航天, 2016, 33(1): 38-41.

    [7] 周文雅, 劉恒, 李順利, 等. 基于空間交會的天基反衛(wèi)作戰(zhàn)方案[J]. 上海航天, 2009, 26(3): 12-17.

    [8] 方飛, 和興鎖. 空間導(dǎo)彈快速攔截最優(yōu)軌道設(shè)計[J]. 山西大學(xué)學(xué)報, 2006, 29(2): 158-161.

    [9] 王衛(wèi)杰, 沈懷榮. 空間飛行器撞擊交匯模型研究[J]. 電子設(shè)計工程, 2011, 19(8): 66-70.

    [10] 宮穎. 衛(wèi)星空間作戰(zhàn)的基本理論和方法研究[D]. 西安: 西北工業(yè)大學(xué), 2005.

    [11] 羅亞中. 空間最優(yōu)交會路徑規(guī)劃策略研究[D]. 長沙: 國防科技大學(xué), 2007.

    [12] 郭海林, 曲廣吉. 空間交會遠程導(dǎo)引階段最優(yōu)變軌策略研究[J]. 載人航天, 2003(3): 28-33.

    [13] 馬曉兵, 潘騰. 空間交會多沖量變軌非線性規(guī)劃優(yōu)化研究與仿真[J]. 載人航天, 2010(1): 21-25.

    [14] 樂晉, 徐世杰, 張崇峰. 空間交會對接尋的段三次變軌策略分析[J]. 上海航天, 2008, 25(5): 37-41, 48.

    猜你喜歡
    交會航天器增量
    “織交會”打造大朗樣板
    2022 年第二季度航天器發(fā)射統(tǒng)計
    國際太空(2022年7期)2022-08-16 09:52:50
    提質(zhì)和增量之間的“辯證”
    “價增量減”型應(yīng)用題點撥
    2019 年第二季度航天器發(fā)射統(tǒng)計
    國際太空(2019年9期)2019-10-23 01:55:34
    2019年,水交會來了!
    2018 年第三季度航天器發(fā)射統(tǒng)計
    國際太空(2018年12期)2019-01-28 12:53:20
    2018年第二季度航天器發(fā)射統(tǒng)計
    國際太空(2018年9期)2018-10-18 08:51:32
    立方體星交會對接和空間飛行演示
    太空探索(2016年9期)2016-07-12 09:59:53
    基于均衡增量近鄰查詢的位置隱私保護方法
    久久热精品热| 九色成人免费人妻av| 在线观看一区二区三区激情| 成人漫画全彩无遮挡| 国产免费福利视频在线观看| 免费黄色在线免费观看| 日日撸夜夜添| 高清视频免费观看一区二区| 日韩不卡一区二区三区视频在线| 国产精品国产av在线观看| 国产精品久久久久久精品电影小说 | 亚洲内射少妇av| 久久国产精品大桥未久av | 免费av中文字幕在线| 国产成人a区在线观看| 18禁在线无遮挡免费观看视频| 亚洲欧美精品专区久久| 午夜免费男女啪啪视频观看| 国产精品久久久久久精品电影小说 | 久久久久久久久久久免费av| 久久久久久久大尺度免费视频| 97在线视频观看| 欧美精品人与动牲交sv欧美| 精品亚洲成a人片在线观看 | 午夜激情福利司机影院| www.色视频.com| 十八禁网站网址无遮挡 | 99国产精品免费福利视频| 老司机影院成人| 18禁裸乳无遮挡动漫免费视频| 黄色配什么色好看| 久久久久久伊人网av| 九草在线视频观看| 欧美丝袜亚洲另类| 亚洲国产欧美在线一区| 国产白丝娇喘喷水9色精品| a级毛色黄片| 欧美性感艳星| 日本色播在线视频| 中文字幕久久专区| 免费观看的影片在线观看| 婷婷色综合大香蕉| 亚洲中文av在线| 韩国av在线不卡| 狂野欧美激情性xxxx在线观看| 日韩欧美 国产精品| 一级毛片aaaaaa免费看小| 国产精品久久久久久精品电影小说 | 日本黄大片高清| 老师上课跳d突然被开到最大视频| 边亲边吃奶的免费视频| 乱码一卡2卡4卡精品| 日日摸夜夜添夜夜爱| 99久久中文字幕三级久久日本| 久久人妻熟女aⅴ| 国产色婷婷99| 久久久久国产精品人妻一区二区| 精品少妇久久久久久888优播| 自拍偷自拍亚洲精品老妇| 中文字幕亚洲精品专区| 日本色播在线视频| 亚洲性久久影院| 国产无遮挡羞羞视频在线观看| 色哟哟·www| 亚洲激情五月婷婷啪啪| 久久久久久人妻| 国产成人a∨麻豆精品| 日本wwww免费看| 肉色欧美久久久久久久蜜桃| 在线观看一区二区三区激情| 日韩 亚洲 欧美在线| 美女主播在线视频| 色5月婷婷丁香| 91精品国产九色| 欧美成人一区二区免费高清观看| 中国三级夫妇交换| 欧美日韩综合久久久久久| 最新中文字幕久久久久| 如何舔出高潮| 欧美一区二区亚洲| 成人黄色视频免费在线看| 久久久久久久久久成人| 人妻系列 视频| 国产高清国产精品国产三级 | av国产免费在线观看| 精品亚洲成a人片在线观看 | 久久精品国产亚洲av天美| 嫩草影院入口| 少妇人妻久久综合中文| 国内少妇人妻偷人精品xxx网站| 国产毛片在线视频| 99久久综合免费| 精品人妻视频免费看| 亚洲在久久综合| 新久久久久国产一级毛片| 人妻少妇偷人精品九色| 精品少妇久久久久久888优播| 免费久久久久久久精品成人欧美视频 | 2021少妇久久久久久久久久久| 噜噜噜噜噜久久久久久91| 午夜福利在线在线| 亚洲国产日韩一区二区| 亚洲人成网站在线播| 老司机影院成人| 99九九线精品视频在线观看视频| 婷婷色麻豆天堂久久| 女性生殖器流出的白浆| 99热这里只有是精品在线观看| 一区二区三区四区激情视频| 少妇丰满av| 全区人妻精品视频| 校园人妻丝袜中文字幕| 91aial.com中文字幕在线观看| 久久久久久久久久久免费av| 又爽又黄a免费视频| 韩国av在线不卡| 一级毛片电影观看| 激情五月婷婷亚洲| 亚洲av电影在线观看一区二区三区| 亚洲精品国产色婷婷电影| 18禁在线播放成人免费| 校园人妻丝袜中文字幕| 中国国产av一级| 国产亚洲欧美精品永久| av天堂中文字幕网| 亚洲人成网站高清观看| 国产v大片淫在线免费观看| 最近手机中文字幕大全| 亚洲欧洲日产国产| 亚洲欧美日韩无卡精品| 人体艺术视频欧美日本| 最近最新中文字幕大全电影3| 91久久精品电影网| 久久精品国产a三级三级三级| 国精品久久久久久国模美| 菩萨蛮人人尽说江南好唐韦庄| 欧美+日韩+精品| 久久久亚洲精品成人影院| 国产精品久久久久久精品古装| 国产黄片视频在线免费观看| 十分钟在线观看高清视频www | 成人免费观看视频高清| 亚洲国产成人一精品久久久| 国产成人a∨麻豆精品| 99久久中文字幕三级久久日本| 国产永久视频网站| 十八禁网站网址无遮挡 | 久久久国产一区二区| 国产永久视频网站| 欧美人与善性xxx| 亚洲欧美日韩另类电影网站 | 国产精品一二三区在线看| 激情五月婷婷亚洲| 亚洲成人中文字幕在线播放| 日本-黄色视频高清免费观看| 成人美女网站在线观看视频| 少妇熟女欧美另类| 亚洲高清免费不卡视频| 成人特级av手机在线观看| 国产精品秋霞免费鲁丝片| 肉色欧美久久久久久久蜜桃| 联通29元200g的流量卡| 国产69精品久久久久777片| videossex国产| 黑人猛操日本美女一级片| 美女主播在线视频| 欧美三级亚洲精品| 在线精品无人区一区二区三 | 精品久久久久久久末码| 菩萨蛮人人尽说江南好唐韦庄| 熟女电影av网| 国产伦理片在线播放av一区| 国产深夜福利视频在线观看| 亚洲欧美日韩另类电影网站 | 国产一区二区在线观看日韩| 亚洲中文av在线| 99久久精品一区二区三区| 久久久久性生活片| 毛片一级片免费看久久久久| 日韩中文字幕视频在线看片 | 日本黄色日本黄色录像| 国产乱人偷精品视频| av女优亚洲男人天堂| 亚洲婷婷狠狠爱综合网| 有码 亚洲区| 国产精品秋霞免费鲁丝片| 只有这里有精品99| 成人一区二区视频在线观看| 亚洲成人手机| 久久久精品94久久精品| 日韩三级伦理在线观看| 成人黄色视频免费在线看| av.在线天堂| 日韩,欧美,国产一区二区三区| 亚洲人成网站在线观看播放| 国产日韩欧美亚洲二区| 欧美丝袜亚洲另类| 好男人视频免费观看在线| 精品久久国产蜜桃| 久久久精品94久久精品| 亚洲av男天堂| 午夜福利在线在线| 男的添女的下面高潮视频| 国产精品一区二区在线观看99| 97超碰精品成人国产| 91精品伊人久久大香线蕉| 成人漫画全彩无遮挡| 欧美成人一区二区免费高清观看| 国产美女午夜福利| 蜜桃亚洲精品一区二区三区| 亚洲美女搞黄在线观看| 国产v大片淫在线免费观看| 亚洲av在线观看美女高潮| 精品一品国产午夜福利视频| 人妻制服诱惑在线中文字幕| 亚洲av中文av极速乱| 国产 一区精品| 免费观看a级毛片全部| 国产成人91sexporn| 亚洲欧洲日产国产| 美女福利国产在线 | 成人无遮挡网站| 日韩 亚洲 欧美在线| 国产爽快片一区二区三区| 国产免费又黄又爽又色| 亚洲,欧美,日韩| 夜夜看夜夜爽夜夜摸| 亚洲精品久久久久久婷婷小说| 精品国产三级普通话版| 有码 亚洲区| 一二三四中文在线观看免费高清| 色婷婷久久久亚洲欧美| 国产欧美日韩一区二区三区在线 | 亚洲成人中文字幕在线播放| 中国三级夫妇交换| 一级毛片久久久久久久久女| 亚洲欧美清纯卡通| 免费看av在线观看网站| 久久久久久人妻| 在线观看免费高清a一片| 一级片'在线观看视频| 亚洲欧美日韩卡通动漫| 国内精品宾馆在线| 免费看日本二区| 精品视频人人做人人爽| 国产欧美亚洲国产| 亚洲第一区二区三区不卡| 欧美日韩视频精品一区| 2018国产大陆天天弄谢| 丰满乱子伦码专区| 91在线精品国自产拍蜜月| 久久久色成人| 18禁动态无遮挡网站| 国产欧美日韩精品一区二区| 久久久久久久久久久免费av| 亚洲精品视频女| 丰满人妻一区二区三区视频av| 成人影院久久| 大香蕉97超碰在线| 如何舔出高潮| 日韩电影二区| 久久久久久久久久人人人人人人| 午夜福利在线在线| 大码成人一级视频| 久久青草综合色| 亚洲国产精品成人久久小说| 99久久人妻综合| 亚洲性久久影院| 久久久久久九九精品二区国产| 久久久久久久久大av| 97超视频在线观看视频| 亚洲av电影在线观看一区二区三区| 亚洲va在线va天堂va国产| 国产无遮挡羞羞视频在线观看| 免费少妇av软件| 岛国毛片在线播放| freevideosex欧美| 日本vs欧美在线观看视频 | 九九久久精品国产亚洲av麻豆| 国产一区亚洲一区在线观看| 一级毛片久久久久久久久女| 中文字幕精品免费在线观看视频 | 少妇熟女欧美另类| 观看免费一级毛片| 亚洲欧美日韩东京热| 国产乱人视频| 亚洲精品色激情综合| 国产av国产精品国产| 国产淫片久久久久久久久| 伊人久久国产一区二区| 人妻一区二区av| 大话2 男鬼变身卡| 久久久成人免费电影| 18禁在线播放成人免费| 久久久久网色| 涩涩av久久男人的天堂| 国产欧美亚洲国产| 一区在线观看完整版| 26uuu在线亚洲综合色| 国产精品av视频在线免费观看| 成人影院久久| 黄色视频在线播放观看不卡| 乱码一卡2卡4卡精品| 亚洲色图av天堂| 久久99蜜桃精品久久| 欧美精品国产亚洲| 欧美日韩国产mv在线观看视频 | 国产成人精品久久久久久| 日韩伦理黄色片| 色视频在线一区二区三区| 最近2019中文字幕mv第一页| 国产色爽女视频免费观看| 久久女婷五月综合色啪小说| 深夜a级毛片| 国产精品精品国产色婷婷| 国产精品一区www在线观看| av在线app专区| 亚洲av二区三区四区| 少妇丰满av| 女人十人毛片免费观看3o分钟| 美女视频免费永久观看网站| kizo精华| 22中文网久久字幕| 国产精品人妻久久久影院| 美女内射精品一级片tv| 国产在线一区二区三区精| 91精品一卡2卡3卡4卡| 国产在线免费精品| 两个人的视频大全免费| 中国三级夫妇交换| av国产久精品久网站免费入址| 高清毛片免费看| 国产精品国产av在线观看| 欧美区成人在线视频| 欧美亚洲 丝袜 人妻 在线| 国产亚洲欧美精品永久| 毛片女人毛片| tube8黄色片| 一级毛片黄色毛片免费观看视频| 亚洲自偷自拍三级| 久久国内精品自在自线图片| 丝瓜视频免费看黄片| 大香蕉久久网| 日韩精品有码人妻一区| 亚洲欧美日韩卡通动漫| 极品教师在线视频| 亚洲在久久综合| 少妇人妻 视频| 校园人妻丝袜中文字幕| 99热这里只有是精品在线观看| 97超碰精品成人国产| 这个男人来自地球电影免费观看 | 日韩一本色道免费dvd| 日本vs欧美在线观看视频 | 国产精品女同一区二区软件| 人妻系列 视频| 亚洲三级黄色毛片| 欧美xxⅹ黑人| 爱豆传媒免费全集在线观看| 日本欧美国产在线视频| 亚洲第一av免费看| 国产av精品麻豆| 丰满人妻一区二区三区视频av| 激情 狠狠 欧美| a级毛色黄片| 国产一区二区三区综合在线观看 | 国产午夜精品久久久久久一区二区三区| 久久影院123| 久久精品久久久久久久性| 亚洲,欧美,日韩| 亚洲成人中文字幕在线播放| 久久青草综合色| 久久久久久人妻| 插逼视频在线观看| 大码成人一级视频| 成年人午夜在线观看视频| 国产淫语在线视频| 午夜福利网站1000一区二区三区| 久久人人爽人人爽人人片va| 亚洲国产精品成人久久小说| 久久精品熟女亚洲av麻豆精品| 久久久久久九九精品二区国产| 伦理电影免费视频| 国产精品国产三级国产专区5o| 国产成人一区二区在线| 国产 一区精品| 噜噜噜噜噜久久久久久91| 亚洲欧美一区二区三区国产| 国产伦理片在线播放av一区| 五月伊人婷婷丁香| 久久精品久久精品一区二区三区| 欧美日韩亚洲高清精品| 性高湖久久久久久久久免费观看| 99久久精品热视频| 99久久精品国产国产毛片| 小蜜桃在线观看免费完整版高清| 国产高潮美女av| 成人影院久久| 欧美xxxx性猛交bbbb| 久久av网站| 久久久久精品久久久久真实原创| 18禁在线无遮挡免费观看视频| 免费黄频网站在线观看国产| 午夜福利高清视频| 久久影院123| 欧美极品一区二区三区四区| 黄色一级大片看看| 天天躁夜夜躁狠狠久久av| 99九九线精品视频在线观看视频| 能在线免费看毛片的网站| 亚洲人成网站在线播| 精品国产露脸久久av麻豆| 2018国产大陆天天弄谢| 男人狂女人下面高潮的视频| 人人妻人人澡人人爽人人夜夜| 久久久久久伊人网av| 日本欧美国产在线视频| 亚洲av中文字字幕乱码综合| 精华霜和精华液先用哪个| 亚洲高清免费不卡视频| 夜夜爽夜夜爽视频| 一区二区三区精品91| 成人无遮挡网站| 成年女人在线观看亚洲视频| 只有这里有精品99| 亚洲成人手机| 美女中出高潮动态图| 中文欧美无线码| 18禁在线无遮挡免费观看视频| 青春草视频在线免费观看| 久久久午夜欧美精品| av卡一久久| 中文字幕久久专区| 免费人妻精品一区二区三区视频| 女人十人毛片免费观看3o分钟| 久久国产精品男人的天堂亚洲 | 三级经典国产精品| 大码成人一级视频| 青春草国产在线视频| 久久久久久久久久人人人人人人| 国产一区亚洲一区在线观看| 青春草国产在线视频| 亚洲成色77777| 多毛熟女@视频| 亚洲av综合色区一区| 精品国产一区二区三区久久久樱花 | 欧美精品一区二区免费开放| 亚洲av免费高清在线观看| 国产日韩欧美亚洲二区| 午夜福利高清视频| av国产精品久久久久影院| 国产乱人偷精品视频| 久久久久久九九精品二区国产| 91久久精品国产一区二区成人| 日韩视频在线欧美| 久久99蜜桃精品久久| 能在线免费看毛片的网站| 亚洲国产精品999| 国产男女内射视频| 亚洲精品乱码久久久v下载方式| 亚洲欧美一区二区三区黑人 | 精品人妻一区二区三区麻豆| 日韩av不卡免费在线播放| 国产成人免费无遮挡视频| 久久国产精品男人的天堂亚洲 | 久久综合国产亚洲精品| 又爽又黄a免费视频| 国产乱人视频| 男男h啪啪无遮挡| 久久久欧美国产精品| 精品国产三级普通话版| 国产精品久久久久久精品电影小说 | 国产伦理片在线播放av一区| 99热国产这里只有精品6| 女性被躁到高潮视频| 熟女av电影| 草草在线视频免费看| 春色校园在线视频观看| 中文资源天堂在线| 国产成人91sexporn| 青春草国产在线视频| 午夜福利网站1000一区二区三区| 青青草视频在线视频观看| 少妇精品久久久久久久| 亚洲熟女精品中文字幕| 老司机影院成人| 国产毛片在线视频| 高清日韩中文字幕在线| 欧美丝袜亚洲另类| 最近中文字幕2019免费版| 在线观看国产h片| 乱码一卡2卡4卡精品| 自拍偷自拍亚洲精品老妇| 久久国内精品自在自线图片| av福利片在线观看| 亚洲aⅴ乱码一区二区在线播放| 亚洲人成网站高清观看| 少妇人妻一区二区三区视频| a 毛片基地| 成人免费观看视频高清| 高清av免费在线| 欧美日韩亚洲高清精品| 日韩一本色道免费dvd| 美女内射精品一级片tv| 国产成人午夜福利电影在线观看| 亚洲国产av新网站| 国产淫片久久久久久久久| 免费大片黄手机在线观看| 国产男女内射视频| 久久久久久久久久久丰满| 在线观看免费日韩欧美大片 | 久久久久国产网址| 午夜激情久久久久久久| 中文字幕制服av| 国产极品天堂在线| 最后的刺客免费高清国语| 国产伦在线观看视频一区| 在线观看三级黄色| 特大巨黑吊av在线直播| 国内少妇人妻偷人精品xxx网站| 最后的刺客免费高清国语| 成年女人在线观看亚洲视频| 久久久成人免费电影| 亚洲成人一二三区av| 22中文网久久字幕| 成人无遮挡网站| 久久久久国产网址| 亚洲国产欧美人成| 内地一区二区视频在线| 男女啪啪激烈高潮av片| 性色avwww在线观看| 日韩电影二区| 婷婷色麻豆天堂久久| 在现免费观看毛片| 亚洲人成网站在线观看播放| 亚洲国产高清在线一区二区三| 欧美zozozo另类| 97超碰精品成人国产| 国产亚洲午夜精品一区二区久久| 国产一区有黄有色的免费视频| 最黄视频免费看| 成年免费大片在线观看| 国产有黄有色有爽视频| 午夜福利视频精品| 久久精品国产亚洲av涩爱| 成人影院久久| 亚洲丝袜综合中文字幕| 久久人人爽人人爽人人片va| 91精品国产九色| 国产午夜精品久久久久久一区二区三区| 亚洲国产精品成人久久小说| 日本免费在线观看一区| videossex国产| 国产精品.久久久| 天堂8中文在线网| 91狼人影院| 伦理电影免费视频| 99久久精品一区二区三区| 国产一区亚洲一区在线观看| 极品少妇高潮喷水抽搐| 少妇人妻久久综合中文| 免费黄网站久久成人精品| 99热国产这里只有精品6| 成人二区视频| 国产精品无大码| 久久精品久久久久久久性| 国产黄色免费在线视频| 秋霞在线观看毛片| 亚洲激情五月婷婷啪啪| av女优亚洲男人天堂| 亚洲欧美精品自产自拍| 能在线免费看毛片的网站| 亚洲成人av在线免费| 国产精品久久久久久av不卡| 国产在线免费精品| 久久97久久精品| 久久鲁丝午夜福利片| 嫩草影院新地址| 成人高潮视频无遮挡免费网站| 久久久久久久精品精品| 欧美老熟妇乱子伦牲交| 久久ye,这里只有精品| 国产在线男女| 国产视频首页在线观看| 久久久欧美国产精品| 国产精品欧美亚洲77777| 欧美变态另类bdsm刘玥| a级毛片免费高清观看在线播放| 免费看不卡的av| 成人影院久久| 国产精品99久久久久久久久| 少妇裸体淫交视频免费看高清| 大香蕉久久网| 亚洲欧美日韩卡通动漫| 人体艺术视频欧美日本| 一本久久精品| 寂寞人妻少妇视频99o| a级毛色黄片| 亚洲天堂av无毛| 国产v大片淫在线免费观看| av免费观看日本| 精品99又大又爽又粗少妇毛片| 99热国产这里只有精品6| 久久久久久人妻| 亚洲av成人精品一二三区| 蜜桃久久精品国产亚洲av| www.色视频.com| 老女人水多毛片| 欧美精品国产亚洲| 一级毛片aaaaaa免费看小| 久久久久人妻精品一区果冻| 国产精品人妻久久久久久| 视频中文字幕在线观看| 亚洲欧美日韩东京热| 国模一区二区三区四区视频|