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

    航天飛行器在共面多軌道間的機動規(guī)劃方法研究

    2021-09-26 07:08:30劉冰雁于鴻源馬心意房瑩張頌
    動力學與控制學報 2021年4期
    關鍵詞:變軌開普勒沖量

    劉冰雁 于鴻源 馬心意 房瑩 張頌

    (1.軍事科學院,北京100091)(2.解放軍32032部隊,北京100094)

    引言

    航天飛行器在太空軌道間快速開展維護、補給、移除和再利用等在軌服務,都需要在所帶燃料支撐能力下依一定順序訪問各目標軌道.太空軌道是當飛行產生的離心力相當于地球地心引力而在一定高度繞地球做圓周運動的運行軌跡.當這些運行軌跡均位于同一軌道平面內時,為共面軌道.地球周圍太空按高度分高、中、低軌道,范圍可達幾萬公里,飛行器在多軌道間機動,為了節(jié)省燃料通常遵循開普勒定律,多采用雙脈沖、多脈沖的機動形式.但在緊急情況下,飛行器需要以任務為重心,在最短時間內執(zhí)行所有在軌服務,多軌道間的機動時間便成了任務成敗的關鍵.因此,面對相距較遠的多條軌道,需結合飛行器所帶燃料選擇變軌方式,規(guī)劃出用時最短、科學合理的機動方案.

    有關軌道機動規(guī)劃,Zaborsky[1]、Pontan[2]和謝成清[3]研究了兩軌道間的轉移軌道設計問題,對分析不同方式下的軌道間機動效果具有借鑒意義.袁建平[4]、李俊峰[5]和孫沖[6]關于非開普勒軌道機動方面進行了探索研究,對軌道間采取非開普勒徑向機動方式的提出具有啟示作用.歐陽琦[7,8]、余婧[9]和朱嘯宇[10]針對在軌服務進行了任務規(guī)劃研究,以燃料消耗最小為出發(fā)點簡化機動過程,重點研究了任務分配和資源調度問題,但較少考慮多軌道間的機動規(guī)劃問題.

    本文將根據(jù)航天飛行器的變軌機動需求,考慮開普勒的霍曼、高斯以及非開普勒徑向三種典型變軌機動方式,構建共面多軌道間的機動規(guī)劃模型,針對有限時間及燃料通過機動策略尋優(yōu)算法解算,以獲得最短機動時間和最優(yōu)機動方案.

    1 軌道間機動方式分析

    太空軌道是空間任務中的基本要素,是航天飛行器僅受引力作用的穩(wěn)態(tài)飛行路徑.軌道間機動是靠動力裝置改變飛行軌道的行為,是以初始軌道為起點通過轉移軌道抵達目標軌道的過程.共面軌道間機動方式的選取需要考慮初始軌道與目標軌道的共面性、切向速度沖量與轉移軌道形態(tài)的約束性及開普勒機動與非開普勒機動的相異性.

    1.1 共面霍曼變軌機動

    在一定離心力下,航天飛行器只受地球引力和軌道攝動影響,將沿著一定軌道飛行.若想從當前軌道跨越到同一軌道平面內的另一軌道,就需要飛行器自身動力系統(tǒng)提供外力改變原有飛行速度或方向.以最省燃料為出發(fā)點,遵循開普勒運動,只需給予飛行器一個速度增量改變運行軌道使其沿轉移軌道飛行,當此轉移軌道正好與目標軌道相切即可完成軌道間的機動,這就是共面霍曼變軌機動的基本原理.

    共面霍曼變軌機動過程如圖1所示,在半徑為r1的初始軌道C1的任意點P產生第一個速度沖量Δv,改變其原有軌道飛至轉移軌道E,E的近拱點為P、遠拱點為A.經過半個轉移周期達到遠拱點為A,遠拱點正好與半徑為r2的目標軌道C2相切.一般軌道轉移問題在考慮近拱點速度沖量vhm改變原來軌位的同時,還需要考慮到遠拱點順利進入目標軌道問題.但我們研究的軌道間機動任務中,到達目標軌位附近后不一定采取入軌行為,有可能根據(jù)操控任務需要采取伴飛、繞飛和懸停等動作.因此,在軌道間機動過程中論文主要考慮初始變軌過程中的速度沖量、燃料損耗和轉移軌道周期問題.

    圖1 共面軌道間機動方式示意Fig.1 A schematic illustration for maneuver mode between coplanar orbits

    設vc1和vc2分別為初始軌道C1和目標軌道C2的速度,橢圓轉移軌道與初始軌道相切的近拱點P的速度為vEP,μ為天體引力常數(shù).共面霍曼變軌機動方式中:轉移軌道的遠拱點正好與目標軌道相切,據(jù)此可以構建軌道速度與軌道半徑的關系模型[11]:

    由此可知,軌道速度主要跟地球引力常數(shù)、軌道半徑相關,軌道轉移速度可以根據(jù)初始軌道速度、半徑以及目標軌道半徑確定.再根據(jù)開普勒方程,可以證明橢圓軌道的運行周期只與長半軸a有關,而與軌道形狀(偏心率)無關.基于以上分析,可推導出共面霍曼速度沖量模型、長半軸模型以及轉移軌道周期模型:

    1.2 共面高斯變軌機動

    通過對軌道動力學和開普勒運動分析可知,速度沖量Δv在一定范圍內遞增,如圖2所示轉移軌道形態(tài)將由橢圓經長橢圓向拋物線演化,變軌機動時間t也將由半個轉移周期T逐漸縮短.在開普勒定律下,只需知道初始軌道與目標軌道位置以及可提供的速度沖量,便可估算變軌機動時間和各軌道參數(shù)的方法均屬于共面高斯變軌機動方式,霍曼變軌機動可視為高斯方式的一種特例.由速度沖量改變產生的各種共面高斯軌道機動方式,與共面霍曼變軌機動相比,能有效縮減軌道間機動時間,但需以增大速度沖量為代價,這對航天飛行器的動力裝置和燃料儲備提出了更高要求.

    圖2 不同速度下的轉移軌道形態(tài)圖Fig.2 Transition orbit morphologies at different velocities

    由此,共面高斯變軌機動方式中的速度沖量Δv是一個不確定量,需根據(jù)動力裝置及燃料儲備提前確定 .可依據(jù)開普勒運動[11,12]建立以速度沖量為變量的轉移軌道近拱點與初始軌道的速度關系模型:

    接著,根據(jù)活力公式可計算求得轉移軌道的各軌道根數(shù)[11].其中,依據(jù)開普勒定律和運動方程M=E-e?sinE,可構建出飛行器從初始軌道變軌機動至目標軌道的機動時間模型:

    1.3 共面非開普勒徑向機動

    目前太空中的飛行器絕大部分軌道運動都遵循1610年提出的開普勒定律,重點考慮地球中心引力場下的飛行器軌道運動.隨著航天技術的發(fā)展、太空在軌服務需求的不斷拓展,有越來越多緊急、復雜、大范圍和多作業(yè)形式的空間操控任務需要飛行器去完成.這樣,是否能夠依任務需要直接、快速機動至目標軌道便成為當前制約任務完成效率的關鍵因素.基于開普勒定律和開普勒空間軌道理論的機動模型具有機動軌道模式固定、運動時間長以及持續(xù)推力運動難描述等局限,難以適應直接、快速和持續(xù)推力下的空間機動行為.未來太空在軌服務不僅要求飛行器能夠順利抵達目標軌道,更加注重變軌機動的快速性、直接性、可控性和全域可達性,這致使整個運動不再完全滿足開普勒定律.我們定義一切不再符合理性開普勒運動的太空機動為非開普勒運動,即為不再符合開普勒三定律或二體問題解的太空運動.在動力裝置持續(xù)推力作用下,將改變飛行器質心運動的軌跡,即轉移軌道不再是開普勒軌道,而是受不同的瞬時加速度做螺旋、弧線或折線甚至直線軌跡[13,14].本文將拋開開普勒定律限制,考慮軌道間最直接的機動方式,采取在徑向方向給予持續(xù)推力的方式使飛行器在徑向產生加速度,從而更快到達目標軌道.

    設在地球引力場中不考慮攝動因素,飛行器在地球萬有引力下做圓周運動,即只受法向方向的引力FG.現(xiàn)依據(jù)目標軌道方位通過動力系統(tǒng)產生一個持續(xù)發(fā)動機推力F,該F可分解為徑向分力Fh和切向分力 Fs,由此可求得飛行器的合加速度[15,16]:

    式中,m為飛行器質量,rd為法向(徑向)單位矢量,sd為切向單位矢量.設M為地球質量、G為引力常量、r為飛行器初始軌道半徑,把μ=GM、FG=GMm r2[17]代入式(7),得到徑向方向的瞬時加速度模型:

    2 共面多軌道間的機動規(guī)劃方法

    飛行器為完成共面多軌道間機動任務需結合動力系統(tǒng)特性和攜帶燃料量選取合適的變軌機動方式,以最短時間機動到目標軌道.其中,軌道高度不同則速度不同,相互間變軌機動所需的速度沖量也就不同,速度沖量的增加雖可有效縮短機動時間,但卻會消耗更多燃料.因此,飛行器在軌道間采取各變軌機動行為之前,需要將各方面因素統(tǒng)籌考慮,構建規(guī)劃模型,設計適用、高效的尋優(yōu)算法.

    2.1 目標優(yōu)化函數(shù)模型

    面對任務執(zhí)行度、飛行裝置動能轉換率和攜帶燃料量等多優(yōu)化需求,從緊急機動角度選取最重要因素作為優(yōu)化目標.據(jù)上文分析可知,緊急在軌服務任務的完成度主要還是受整個軌道間機動時間影響,時間就是生命,只有在機動上節(jié)省更多的時間才能在每個任務點上留有更充分的操控時間.因此,在共面多軌道間的機動規(guī)劃問題上,選取整個機動時間為優(yōu)化目標并構建目標優(yōu)化函數(shù):

    式中:tij為從i軌道機動至j軌道的時間,單位為s;xij為i軌道與j軌道間轉移軌道的選取情況,為無量綱量;N為需訪問的軌道總數(shù),單位為條.

    2.2 變軌機動約束模型

    變軌機動前需選定變軌方式,依上文分析可大致分為開普勒運動和非開普勒運動,具體選取依動力裝置性能、任務緊急程度和燃料儲量而定.根據(jù)開普勒運動的軌道根數(shù)、速度沖量、軌道轉移周期以及機動時間的動力學關系,建立開普勒高斯變軌機動約束模型:

    式中,aij為從i軌道機動至j軌道的轉移軌道長半軸,單位為km;Eij為ij軌道間轉移軌道的偏近點角,單位為度;eij為ij軌道間轉移軌道的偏心率,為無量綱量;vij為i軌道機動至j軌道的轉移軌道近拱點速度,單位為k∕ms;Δvij為i軌道機動至j軌道的速度沖量,單位為k∕ms;Δυ為基于共面霍曼速度沖量的速度增量,單位為k∕ms.

    非開普勒運動軌跡主要由飛行器在徑向及切向方向上的瞬時加速度共同決定,軌道間機動時間主要受徑向方向的加速度控制.由此,重點考慮徑向方向機動情況,選取均變速徑向機動這一特殊情況建立非開普勒變軌約束模型:

    式中,αij為從i軌道機動至j軌道的徑向合加速度,單位為km∕s2;rij為從i軌道機動至j軌道過程中距地心的徑向距離,單位為km;為i軌道機動至j軌道過程中的徑向發(fā)動機推力,單位為N;vi0為i軌道機動至 j軌道的徑向初始速度,單位為 k∕ms;vij為i軌道機動至 j軌道的機動速度,單位為 k∕ms;sij為i、j軌道間機動的徑向累計路程,單位為km;tij為從i軌道機動至j軌道的時間,單位為s.

    2.3 變軌路徑約束模型

    飛行器面對共面多軌道,需事先明確訪問各軌道的先后順序且每個軌道只訪問一次.這是一類組合優(yōu)化的問題約束,即針對每一目標軌道需明確:僅經由一條轉移軌道進入和一條轉移軌道離開,整個轉移軌道構成一個轉移鏈,中間沒有任何子回路的產生.由此,構建整個變軌機動的路徑約束模型:

    式中,當xij=1時,選取ij軌道之間的轉移軌道,其他則不選取.

    2.4 變軌燃料約束模型

    有的太空軌道間距較遠,航天飛行器自身燃料補給困難,機動過程中必須考慮燃料的消耗情況.可將每次變軌機動前后速度矢量改變量的模,即速度增量的大小作為軌道控制所付出能量代價的一種間接度量.基于動量守恒原理齊奧爾科夫斯基給出了速度與燃料損耗的關系理論[18,19],即航天飛行器通過一個消耗自身燃料質量的反方向推進系統(tǒng),可以在原有運行速度上產生并獲得加速度.由此,構建軌道間機動燃料約束模型:

    式中,vij為i軌道機動至j軌道開普勒運動的近拱點速度或非開普勒運動的機動速度,單位為km∕s2;β為從i軌道機動至j軌道過程中燃料轉換率,是由動力系統(tǒng)決定的常量;λ為i軌道機動至j軌道過程中燃料燃燒后噴出的能量相對火箭的速率,單位為km∕s;mq為飛行器本體的質量,單位為kg;mr為飛行器所攜帶的燃料總量,單位為kg;mij為從i軌道機動至j軌道所消耗的燃料,單位為kg.

    2.5 機動策略尋優(yōu)算法

    根據(jù)所建立的機動規(guī)劃模型,設計基于改進蟻群算法的機動策略尋優(yōu)算法.算法中根據(jù)要訪問的軌道數(shù)量設計螞蟻數(shù),其中每只螞蟻具有特征:根據(jù)以每次變軌機動時間和轉移軌道上外激素的數(shù)量為變量的概率函數(shù),選定下一個訪問軌位;強制螞蟻遵循合法路徑搜索,不重復訪問、不走小環(huán)路,并由禁忌表進行控制;完整一次搜尋后,在曾經訪問的每一條轉移軌道用外激素標記.算法起始時刻,設各軌道間轉移軌道所攜帶的信息素相等,各螞蟻在搜索中根據(jù)轉移軌道上的所需的機動時間決定轉移方向,并按轉移概率模型機動至下一軌道.遍歷完所有軌道后,螞蟻完成一次循環(huán),根據(jù)外激素更新模型[20]修改軌跡強度,并以燃料消耗情況審查進行約束.整個解算流程如圖3所示.

    圖3 基于改進蟻群算法的機動策略尋優(yōu)算法流程圖Fig.3 Flow chart of maneuver strategy optimization algorithm based on improved ant colony algorithm

    3 算例求解與分析

    算例是以航天飛行器在赤道面上的七條軌道間往返執(zhí)行在軌服務任務為背景,將七條軌道簡化為圓形軌道,各軌道的距地球面高度集為:{r1=36300;r2=36700;r3=37200;r4=37800;r5=38500;r6=39300},單位為km.航天飛行器本體質量mq=1× 103kg,攜帶燃料 mr=4×103kg,裝備最大推力為Fmax=39.3kN的肼燃料改進型火箭發(fā)動機,動力系統(tǒng)噴口速率λ=3km∕s,采取非開普勒徑向機動方式時徑向方向最大合加速度αmax=6×10-3km∕s2,地球引力常數(shù) μ =398600.4km3∕s2,地球半徑rε=6378.14km.現(xiàn)需此飛行器在這七條軌道依序開展緊急在軌服務任務,其中預留的軌道間總機動時間只有2.5小時,且要求燃料最少結余50公斤以保障在軌操控任務的順利開展.

    3.1 兩類變軌方式的機動時長及燃料消耗情況

    軌道間采取高斯機動過程中,增大速度沖量將有效提升機動時效性,但燃料消耗也有所增加.圖4仿真展現(xiàn)了r0機動至r1時長隨速度沖量增大呈非線性遞減的情況,及兩軌道間燃料消耗隨速度沖量增大呈非線性遞增的情況.

    圖4 高斯機動中時間、燃料變化情況Fig.4 Time and fuel variations in Gauss maneuver

    當軌道間采取非開普勒徑向機動時,使航天飛行器在徑向方向保持不同合加速度,其機動時間和燃料消耗情況均有所不同.圖5展現(xiàn)了初始軌道r0在徑向方向采取不同合加速度機動至r1軌道的機動時間、燃料消耗情況.其中,機動時長隨徑向合加速度增大呈非線性遞減趨勢,燃料消耗量隨徑向合加速度增大呈非線性遞增趨勢.

    圖5 非開普勒徑向機動中時間、燃料變化情況Fig.5 Time and fuel variations in non-Keplerian radial maneuver

    3.2 四種特殊情況的最優(yōu)機動方案

    如表1所示,為四種特殊變軌機動方式下的最優(yōu)機動方案.其中,采取開普勒霍曼機動時,各軌道間的轉移軌道均為與目標軌道相切的橢圓軌道,例如r3至r4的轉移軌道形態(tài)為長半軸為41.178×103km、短半軸為41.177×103km的橢圓.通過規(guī)劃方法得到時間最優(yōu)的機動策略,整個機動過程所需時間為136.5小時,燃料消耗量為243.77公斤.

    表1 考慮不同變軌方式的最優(yōu)機動方案Table 1 Optimal maneuver scheme for different ways of orbit shifting

    軌道間的速度沖量受初始軌道高度限制,由此可以確定各速度沖量的極限值.當采取極限速度下的高斯機動時,各軌道間的轉移軌道都是在不同程度上與目標軌道相交的橢圓軌道,例如r3至r4的轉移軌道形態(tài)為長半軸為4513.0×103km、短半軸為1920.9×103km的長條形橢圓.由于軌道間機動時間受初始軌道高度、速度沖量影響,當速度沖量改變則軌道間的機動比較優(yōu)勢將可能有所變化.通過規(guī)劃方法可以得到當前環(huán)境下的最優(yōu)機動策略,即由2號軌道經1號、0號、3號、4號和5號機動至6號軌道.整個機動過程所需時間為8.8小時,燃料消耗量為3701.35公斤.

    將軌道間機動分解為徑向運動跟切向運動,只有當徑向方向運動率先與目標軌道交會才能將整個機動時間以徑向時間計算.由此,徑向機動存在一個最低合加速度,只有當大于等于此臨界值時整個機動時間才能以徑向機動時間計算.根據(jù)軌道間切向運動的特性,結合軌道幾何關系及各初始軌道速度,可確定徑向合加速度的最小值.根據(jù)本算例軌道條件,計算得到徑向最低合加速度為αmin=0.25×10-3km∕s2.基于此,軌道間的徑向機動累計路程為兩軌道的軌道高度差,例如r3至r4的機動路程為600km.通過規(guī)劃方法計算得到時間最優(yōu)的機動策略,即由3號軌道經2號、1號、0號、4號和5號機動至6號軌道.整個機動過程所需時間為8.61小時,燃料消耗量為2703.94千克.

    根據(jù)算例所給條件,當在徑向方向保持最大合加速度αmax=6×10-3km∕s2機動時,兩軌道間的徑向機動累計路程為兩軌道的軌道高度差,例如r3至r4的機動路程為600km.通過規(guī)劃方法計算得到時間最優(yōu)的機動策略,即由2號軌道經1號、0號、3號、4號和5號機動至6號軌道.整個機動過程所需時間為1.76小時,燃料消耗量為3983.99千克.

    3.3 時間及燃料約束給定的最優(yōu)機動方案

    根據(jù)算例對機動總時長控制在2.5小時內、燃料結余50公斤的要求,通過本文方法計算,最終得到了如表2所示的最優(yōu)機動方案.軌道間在徑向保持合加速度做徑向均變速機動,采取由2號軌道經1號、0號、3號、4號和5號機動至6號軌道的最優(yōu)訪問策略;整個機動時間為2.22小時,滿足2.5小時內完成軌道間機動的任務要求;燃料消耗3949.12公斤,滿足結余50公斤燃料的在軌服務任務要求.

    表2 時間及燃料限制條件下的最優(yōu)機動方案Table 2 Optimal maneuver scheme under time and fuel constraints

    從整個計算過程知:開普勒機動過程中,速度沖量決定了轉移軌道形態(tài),速度沖量的增加能縮短軌道間機動時間;通過在各初始軌道增加其速度沖量至極限值,可壓縮93.6%的機動時間,但需多消耗15倍的燃料.非開普勒徑向機動中,在徑向保持最低加速度進行機動,機動時間跟極速高斯相近,由于其機動方式直接且所需速度相對較小,從而整體能夠節(jié)省1000kg燃料;提升徑向加速度至算例所給極值,可有效壓縮79.6%機動時長,但需再多消耗3.9倍燃料.針對機動時間2.5小時內及燃料結余50公斤的條件限制,綜合考慮不同機動方式,本文給出了采取非開普勒徑向機動方式的最優(yōu)方案.

    由此,多軌道間機動采取開普勒方式,整體相對節(jié)省燃料,但所需的機動時間較長,適用于攜帶燃料量有限或任務時間較寬松的情況.采取非開普勒機動方式,可進一步縮減機動時間,適用于燃料充足或任務緊急的情況.從應急角度看,是需要飛行器以這類方式按最優(yōu)機動策略,在最短時間內直接完成多軌道間機動的.但非開普勒運動中,合加速度的提升、推力大小的維持及飛行方向的控制均對飛行器的動力系統(tǒng)、燃料儲備和飛行控制提出了更高的要求,這急需航天技術進一步發(fā)展以提升空間全域快速機動能力.

    4 結論

    (1)提出了共面多軌道間機動規(guī)劃方法,能夠根據(jù)具體需求給出最優(yōu)機動策略.

    (2)研究了兩類變軌機動方式四種特殊情況的邊界問題,有助于受時間及燃料限制機動問題的變軌方式選取.

    (3)提出了非開普勒徑向機動方式,能夠有效縮減軌道間機動時間,滿足高時效性要求.

    綜上,本文方法能夠為受時間及燃料條件限制的共面多軌道間機動問題,提出一個綜合考量不同變軌方式的最優(yōu)機動方案.

    猜你喜歡
    變軌開普勒沖量
    連續(xù)3年銷量翻番,2022年欲沖量4000萬,福建這家動保企業(yè)正強勢崛起
    當代水產(2022年3期)2022-04-26 14:27:18
    “朱諾”變軌時間將推至明年2月
    太空探索(2016年12期)2016-07-18 11:13:43
    例析人造衛(wèi)星的圓周運動及變軌問題
    人造衛(wèi)星變軌問題
    固體運載火箭變軌發(fā)動機噴管氣流分離研究
    開普勒的最新發(fā)現(xiàn):第二地球
    “死而復生”的開普勒
    開普勒的大發(fā)現(xiàn)
    孤獨遠行的開普勒
    變力沖量的教學研究
    物理通報(2011年8期)2011-01-24 07:39:38
    一二三四社区在线视频社区8| 久久久国产精品麻豆| 十八禁人妻一区二区| 国产成人欧美在线观看| 免费看光身美女| 两个人的视频大全免费| 波多野结衣高清作品| 亚洲精品久久国产高清桃花| 伦理电影免费视频| 国产精品久久久久久久电影 | 精品久久久久久,| 久久香蕉国产精品| av女优亚洲男人天堂 | 波多野结衣高清作品| 极品教师在线免费播放| 99久久无色码亚洲精品果冻| av片东京热男人的天堂| 一区二区三区国产精品乱码| 波多野结衣高清无吗| 亚洲精品粉嫩美女一区| 精品无人区乱码1区二区| 一级毛片女人18水好多| 国产高清视频在线播放一区| 波多野结衣高清作品| 免费观看的影片在线观看| 精品久久久久久久久久久久久| 激情在线观看视频在线高清| 99re在线观看精品视频| 亚洲国产欧美一区二区综合| 国产激情久久老熟女| 亚洲最大成人中文| 亚洲美女视频黄频| 俺也久久电影网| svipshipincom国产片| 免费无遮挡裸体视频| 巨乳人妻的诱惑在线观看| 俄罗斯特黄特色一大片| 亚洲精品乱码久久久v下载方式 | 大型黄色视频在线免费观看| 舔av片在线| 亚洲人成伊人成综合网2020| 精品久久久久久久人妻蜜臀av| 麻豆久久精品国产亚洲av| 热99在线观看视频| 久久久久九九精品影院| 免费观看精品视频网站| 日韩免费av在线播放| 天堂√8在线中文| 国产亚洲av嫩草精品影院| 法律面前人人平等表现在哪些方面| 窝窝影院91人妻| 白带黄色成豆腐渣| 国产一级毛片七仙女欲春2| 黄色成人免费大全| 全区人妻精品视频| 欧美日韩黄片免| 午夜免费成人在线视频| 亚洲精品一区av在线观看| 国产精品久久久av美女十八| 久久亚洲精品不卡| 精品国产美女av久久久久小说| 久久久久久久久久黄片| 床上黄色一级片| 免费看美女性在线毛片视频| 色精品久久人妻99蜜桃| 精品无人区乱码1区二区| 国产精品 欧美亚洲| 搞女人的毛片| 99国产极品粉嫩在线观看| 岛国视频午夜一区免费看| 少妇裸体淫交视频免费看高清| 国模一区二区三区四区视频 | 美女高潮喷水抽搐中文字幕| 18禁黄网站禁片午夜丰满| 国产1区2区3区精品| 成人av一区二区三区在线看| 久久久久久大精品| 欧美另类亚洲清纯唯美| 精品一区二区三区视频在线观看免费| 长腿黑丝高跟| 97超级碰碰碰精品色视频在线观看| 91在线观看av| 亚洲欧美日韩无卡精品| a级毛片在线看网站| 级片在线观看| 精品国产超薄肉色丝袜足j| 亚洲av美国av| 国内精品久久久久久久电影| 99久久99久久久精品蜜桃| 久久久久亚洲av毛片大全| 亚洲性夜色夜夜综合| 欧美色欧美亚洲另类二区| 久久久久免费精品人妻一区二区| 国产爱豆传媒在线观看| 哪里可以看免费的av片| 女警被强在线播放| 国产精品一区二区免费欧美| 亚洲电影在线观看av| 波多野结衣高清作品| 日韩有码中文字幕| 欧美精品啪啪一区二区三区| 国产欧美日韩一区二区三| 久久性视频一级片| 91九色精品人成在线观看| 国产精品久久久av美女十八| 一边摸一边抽搐一进一小说| 国产精品野战在线观看| 18禁黄网站禁片免费观看直播| 成人午夜高清在线视频| 观看免费一级毛片| 国产男靠女视频免费网站| 国产伦人伦偷精品视频| 午夜福利成人在线免费观看| 免费无遮挡裸体视频| 91麻豆精品激情在线观看国产| 桃红色精品国产亚洲av| 亚洲无线观看免费| 久久天堂一区二区三区四区| 久久午夜亚洲精品久久| 久久精品亚洲精品国产色婷小说| 国产精品av视频在线免费观看| 久久精品91无色码中文字幕| 两个人看的免费小视频| 99久久久亚洲精品蜜臀av| 亚洲狠狠婷婷综合久久图片| 人妻久久中文字幕网| 欧美日韩综合久久久久久 | 久久伊人香网站| 日本黄色视频三级网站网址| 综合色av麻豆| 一卡2卡三卡四卡精品乱码亚洲| 在线观看免费视频日本深夜| 变态另类丝袜制服| 成人亚洲精品av一区二区| 国产视频一区二区在线看| 国产91精品成人一区二区三区| 一本一本综合久久| 国产精品乱码一区二三区的特点| 婷婷六月久久综合丁香| 精品一区二区三区视频在线 | 免费看a级黄色片| 岛国在线免费视频观看| 欧美最黄视频在线播放免费| 人人妻,人人澡人人爽秒播| 午夜亚洲福利在线播放| 国内精品久久久久精免费| 欧美+亚洲+日韩+国产| 男女之事视频高清在线观看| 国产av一区在线观看免费| 国内精品美女久久久久久| 九九热线精品视视频播放| 色综合欧美亚洲国产小说| 国产亚洲精品久久久久久毛片| 哪里可以看免费的av片| 美女扒开内裤让男人捅视频| 久久精品91蜜桃| 露出奶头的视频| 特大巨黑吊av在线直播| 每晚都被弄得嗷嗷叫到高潮| 成在线人永久免费视频| 美女扒开内裤让男人捅视频| 日韩精品青青久久久久久| 制服丝袜大香蕉在线| www.熟女人妻精品国产| 久久人人精品亚洲av| 91久久精品国产一区二区成人 | 国产伦一二天堂av在线观看| 日韩欧美一区二区三区在线观看| 免费在线观看成人毛片| 久久这里只有精品19| 国产精品亚洲美女久久久| 日本 av在线| 久久午夜亚洲精品久久| 日韩国内少妇激情av| 丰满的人妻完整版| 欧美丝袜亚洲另类 | 午夜免费激情av| 久久久久久久久中文| 国产一区在线观看成人免费| 无人区码免费观看不卡| 成人国产综合亚洲| 欧美日韩中文字幕国产精品一区二区三区| 亚洲欧洲精品一区二区精品久久久| 精品电影一区二区在线| 18美女黄网站色大片免费观看| 脱女人内裤的视频| 在线永久观看黄色视频| 在线观看免费午夜福利视频| 三级男女做爰猛烈吃奶摸视频| 国产av不卡久久| 99久久综合精品五月天人人| 91字幕亚洲| cao死你这个sao货| 亚洲午夜理论影院| 美女高潮喷水抽搐中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 身体一侧抽搐| 国产精品久久久久久人妻精品电影| 91老司机精品| 在线观看一区二区三区| 99热6这里只有精品| 亚洲色图av天堂| 久久香蕉国产精品| 国产 一区 欧美 日韩| 国产1区2区3区精品| 小说图片视频综合网站| 久久久国产精品麻豆| 一个人免费在线观看的高清视频| 精品国产乱码久久久久久男人| 国产av在哪里看| 一级a爱片免费观看的视频| 少妇人妻一区二区三区视频| 日本黄大片高清| 1000部很黄的大片| avwww免费| 男女做爰动态图高潮gif福利片| 久久亚洲精品不卡| 欧美性猛交╳xxx乱大交人| 日韩三级视频一区二区三区| 久久亚洲真实| 亚洲欧美日韩高清专用| 狂野欧美激情性xxxx| 国产成人av激情在线播放| 国产真实乱freesex| 99热6这里只有精品| 亚洲av成人精品一区久久| 免费看光身美女| 日日摸夜夜添夜夜添小说| 成人欧美大片| 99视频精品全部免费 在线 | www.精华液| 麻豆成人午夜福利视频| 亚洲黑人精品在线| 国产精华一区二区三区| 熟女人妻精品中文字幕| 给我免费播放毛片高清在线观看| 国产单亲对白刺激| 国产亚洲精品久久久久久毛片| 亚洲熟妇中文字幕五十中出| 成年女人毛片免费观看观看9| 国产视频一区二区在线看| 久久久久久久午夜电影| 亚洲av电影不卡..在线观看| 亚洲真实伦在线观看| 一边摸一边抽搐一进一小说| 亚洲专区字幕在线| 99在线人妻在线中文字幕| 婷婷丁香在线五月| 日韩精品青青久久久久久| 免费搜索国产男女视频| 99久久99久久久精品蜜桃| 久久精品夜夜夜夜夜久久蜜豆| 一级作爱视频免费观看| 亚洲欧美精品综合一区二区三区| 香蕉丝袜av| 久久国产精品人妻蜜桃| 99re在线观看精品视频| 很黄的视频免费| 三级男女做爰猛烈吃奶摸视频| 综合色av麻豆| 国产探花在线观看一区二区| 他把我摸到了高潮在线观看| 亚洲第一欧美日韩一区二区三区| www.自偷自拍.com| 国产亚洲av嫩草精品影院| 国产伦一二天堂av在线观看| 国产一区二区在线av高清观看| 一卡2卡三卡四卡精品乱码亚洲| 一进一出抽搐gif免费好疼| 老鸭窝网址在线观看| 一区二区三区国产精品乱码| 欧美黄色片欧美黄色片| www.精华液| av天堂在线播放| 亚洲五月婷婷丁香| 国语自产精品视频在线第100页| 久久精品影院6| 一卡2卡三卡四卡精品乱码亚洲| 亚洲成人久久性| 久久精品影院6| 首页视频小说图片口味搜索| 身体一侧抽搐| 国产精品自产拍在线观看55亚洲| 99久久成人亚洲精品观看| 人人妻人人澡欧美一区二区| 国产精品美女特级片免费视频播放器 | 欧美日韩综合久久久久久 | 国产精品野战在线观看| 美女高潮喷水抽搐中文字幕| 高清在线国产一区| 国产精品香港三级国产av潘金莲| 成年女人毛片免费观看观看9| 精品电影一区二区在线| 后天国语完整版免费观看| 国产精品一区二区三区四区免费观看 | 日韩欧美免费精品| 国产亚洲精品av在线| 美女高潮喷水抽搐中文字幕| 制服人妻中文乱码| 国产伦精品一区二区三区视频9 | 国产在线精品亚洲第一网站| 国产免费av片在线观看野外av| 欧美中文日本在线观看视频| 日本五十路高清| 最近在线观看免费完整版| 国语自产精品视频在线第100页| 国内毛片毛片毛片毛片毛片| 欧美日本视频| 欧美日韩中文字幕国产精品一区二区三区| 人妻丰满熟妇av一区二区三区| 草草在线视频免费看| e午夜精品久久久久久久| 午夜福利在线观看免费完整高清在 | 一a级毛片在线观看| 国内精品久久久久久久电影| 国产乱人视频| 欧美日韩一级在线毛片| 麻豆久久精品国产亚洲av| 久久精品91无色码中文字幕| 婷婷亚洲欧美| 91在线精品国自产拍蜜月 | 欧美日本亚洲视频在线播放| 99久久精品国产亚洲精品| 久久久久久人人人人人| 最近在线观看免费完整版| 熟女电影av网| 日韩高清综合在线| 小说图片视频综合网站| 欧美xxxx黑人xx丫x性爽| 不卡av一区二区三区| 午夜福利免费观看在线| 日韩欧美在线乱码| 嫩草影院精品99| 窝窝影院91人妻| 97超视频在线观看视频| 久久热在线av| 欧美黑人欧美精品刺激| 91老司机精品| 99精品欧美一区二区三区四区| 色尼玛亚洲综合影院| 搡老岳熟女国产| x7x7x7水蜜桃| 男人舔奶头视频| 国产欧美日韩一区二区精品| 身体一侧抽搐| 99精品久久久久人妻精品| 女同久久另类99精品国产91| 国产成人影院久久av| 欧美日韩精品网址| 夜夜爽天天搞| 国产黄片美女视频| 成人无遮挡网站| 国产精品99久久99久久久不卡| 欧美成人性av电影在线观看| 亚洲va日本ⅴa欧美va伊人久久| 中文字幕熟女人妻在线| 久久精品国产综合久久久| 国产一区二区在线av高清观看| 天堂网av新在线| 久久精品国产清高在天天线| 天堂影院成人在线观看| 免费电影在线观看免费观看| 最近最新中文字幕大全电影3| 久久婷婷人人爽人人干人人爱| 欧美成人免费av一区二区三区| 亚洲精品一区av在线观看| 免费在线观看亚洲国产| 午夜影院日韩av| 国产日本99.免费观看| 神马国产精品三级电影在线观看| 国产日本99.免费观看| 伊人久久大香线蕉亚洲五| 欧美乱妇无乱码| 久久香蕉精品热| 搡老妇女老女人老熟妇| 欧美精品啪啪一区二区三区| 99国产精品一区二区蜜桃av| 手机成人av网站| 亚洲国产精品成人综合色| 757午夜福利合集在线观看| 国产高清有码在线观看视频| 老司机在亚洲福利影院| 免费看十八禁软件| 国产午夜福利久久久久久| 午夜免费激情av| 激情在线观看视频在线高清| 久久亚洲真实| 熟女电影av网| 好男人在线观看高清免费视频| 亚洲一区二区三区色噜噜| 亚洲成av人片免费观看| 色哟哟哟哟哟哟| 最近在线观看免费完整版| 免费在线观看影片大全网站| 国产成人aa在线观看| ponron亚洲| 久99久视频精品免费| 国产爱豆传媒在线观看| 国产日本99.免费观看| 日本与韩国留学比较| 欧美丝袜亚洲另类 | 校园春色视频在线观看| 久久婷婷人人爽人人干人人爱| 高清在线国产一区| 国产成年人精品一区二区| 亚洲av成人精品一区久久| 亚洲精品美女久久久久99蜜臀| 91老司机精品| 亚洲av日韩精品久久久久久密| 可以在线观看毛片的网站| 亚洲国产欧美人成| 9191精品国产免费久久| 香蕉丝袜av| 亚洲最大成人中文| 免费看十八禁软件| 我的老师免费观看完整版| 别揉我奶头~嗯~啊~动态视频| 99久久久亚洲精品蜜臀av| 欧美日韩一级在线毛片| 亚洲精品一区av在线观看| 久久精品国产99精品国产亚洲性色| 啦啦啦免费观看视频1| 国产v大片淫在线免费观看| 国产一区二区在线av高清观看| xxx96com| 国产一区二区在线av高清观看| 亚洲中文日韩欧美视频| 免费大片18禁| h日本视频在线播放| 十八禁人妻一区二区| 欧美日韩黄片免| 精品熟女少妇八av免费久了| 亚洲真实伦在线观看| 老司机午夜福利在线观看视频| 最新中文字幕久久久久 | 97人妻精品一区二区三区麻豆| 欧美成人免费av一区二区三区| 久久这里只有精品中国| 亚洲第一电影网av| 波多野结衣高清作品| 麻豆一二三区av精品| 亚洲九九香蕉| 成年版毛片免费区| 一本一本综合久久| 国产99白浆流出| 成年人黄色毛片网站| 91字幕亚洲| 色视频www国产| 精华霜和精华液先用哪个| www.精华液| 国产成年人精品一区二区| 成人精品一区二区免费| av视频在线观看入口| 久久国产精品影院| 可以在线观看的亚洲视频| 最好的美女福利视频网| 国产激情久久老熟女| 婷婷六月久久综合丁香| 特大巨黑吊av在线直播| 亚洲精品一区av在线观看| 亚洲人成伊人成综合网2020| 日本免费一区二区三区高清不卡| 成人性生交大片免费视频hd| 我的老师免费观看完整版| 亚洲人成电影免费在线| 国产高潮美女av| 美女cb高潮喷水在线观看 | 免费观看的影片在线观看| 国产亚洲精品久久久久久毛片| 成人无遮挡网站| 黄色成人免费大全| 一个人观看的视频www高清免费观看 | 国产伦在线观看视频一区| 天堂√8在线中文| 观看免费一级毛片| 国产不卡一卡二| 老司机午夜福利在线观看视频| av中文乱码字幕在线| 18禁黄网站禁片午夜丰满| 999久久久精品免费观看国产| 国产精品久久视频播放| 免费看光身美女| 日韩欧美免费精品| 搡老熟女国产l中国老女人| 久久久国产精品麻豆| 日日夜夜操网爽| 长腿黑丝高跟| 国内精品久久久久精免费| 免费电影在线观看免费观看| 中文字幕最新亚洲高清| 丰满的人妻完整版| 日本三级黄在线观看| 欧美在线黄色| 少妇的丰满在线观看| 久久久国产欧美日韩av| 色视频www国产| 天天躁狠狠躁夜夜躁狠狠躁| 日本黄色片子视频| 久久欧美精品欧美久久欧美| 国产单亲对白刺激| 亚洲av免费在线观看| 亚洲午夜理论影院| 老鸭窝网址在线观看| 国产高潮美女av| 久久久久国产一级毛片高清牌| 日韩欧美三级三区| 91九色精品人成在线观看| 18美女黄网站色大片免费观看| 两个人的视频大全免费| 免费看日本二区| 国产精品,欧美在线| 亚洲成av人片在线播放无| 中亚洲国语对白在线视频| 99久久久亚洲精品蜜臀av| 中出人妻视频一区二区| 一二三四社区在线视频社区8| 一本综合久久免费| 免费观看的影片在线观看| 最好的美女福利视频网| 国产 一区 欧美 日韩| 国产精品自产拍在线观看55亚洲| 黑人巨大精品欧美一区二区mp4| 中文字幕精品亚洲无线码一区| 最新中文字幕久久久久 | 午夜免费观看网址| 亚洲国产日韩欧美精品在线观看 | 天堂影院成人在线观看| 噜噜噜噜噜久久久久久91| 久久国产精品影院| 日本撒尿小便嘘嘘汇集6| 色综合亚洲欧美另类图片| 欧美高清成人免费视频www| 欧美xxxx黑人xx丫x性爽| 日韩欧美国产一区二区入口| 日本在线视频免费播放| 亚洲av电影在线进入| 一二三四在线观看免费中文在| 熟女人妻精品中文字幕| 婷婷六月久久综合丁香| 亚洲午夜理论影院| 男人的好看免费观看在线视频| 午夜影院日韩av| 国产精品一区二区免费欧美| 少妇裸体淫交视频免费看高清| 成年女人看的毛片在线观看| 免费看日本二区| АⅤ资源中文在线天堂| 亚洲色图av天堂| 可以在线观看的亚洲视频| 亚洲美女视频黄频| netflix在线观看网站| 欧美xxxx黑人xx丫x性爽| 日本五十路高清| 久久热在线av| 一夜夜www| 老汉色∧v一级毛片| 91麻豆av在线| 麻豆国产97在线/欧美| 国产精品一区二区三区四区免费观看 | 亚洲中文日韩欧美视频| 亚洲国产欧美一区二区综合| 嫩草影视91久久| 法律面前人人平等表现在哪些方面| 国产一区二区在线av高清观看| 一个人看的www免费观看视频| 黄色片一级片一级黄色片| 国产精品av视频在线免费观看| 亚洲国产欧美网| 在线观看午夜福利视频| 午夜影院日韩av| 人妻丰满熟妇av一区二区三区| 国产高清有码在线观看视频| 亚洲专区中文字幕在线| 欧美中文日本在线观看视频| 变态另类丝袜制服| 亚洲精品一卡2卡三卡4卡5卡| 欧美色欧美亚洲另类二区| 美女黄网站色视频| 亚洲自偷自拍图片 自拍| av女优亚洲男人天堂 | 看免费av毛片| 亚洲黑人精品在线| av天堂在线播放| 91老司机精品| 99久久无色码亚洲精品果冻| 伊人久久大香线蕉亚洲五| 久久午夜综合久久蜜桃| 亚洲七黄色美女视频| 成人国产一区最新在线观看| www.熟女人妻精品国产| 久久中文字幕一级| 国产高清视频在线播放一区| 国产精品 欧美亚洲| 亚洲九九香蕉| 午夜视频精品福利| 日韩欧美在线二视频| 99国产综合亚洲精品| 午夜福利在线观看免费完整高清在 | 亚洲aⅴ乱码一区二区在线播放| 91字幕亚洲| xxx96com| 99久久无色码亚洲精品果冻| 女人高潮潮喷娇喘18禁视频| 国产成人av激情在线播放| 男女视频在线观看网站免费| 午夜免费成人在线视频| 怎么达到女性高潮| 天天躁日日操中文字幕| 久久久国产欧美日韩av| 免费一级毛片在线播放高清视频| 俺也久久电影网| www日本黄色视频网| 国产激情偷乱视频一区二区| 男插女下体视频免费在线播放| 91字幕亚洲|