• <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
    国产伦精品一区二区三区视频9| 日韩成人伦理影院| 亚洲av不卡在线观看| 晚上一个人看的免费电影| 久久精品国产亚洲网站| 狂野欧美白嫩少妇大欣赏| 国产黄色视频一区二区在线观看| 免费黄色在线免费观看| 建设人人有责人人尽责人人享有的 | 下体分泌物呈黄色| 久久久久精品性色| 乱系列少妇在线播放| 精品一区二区三区视频在线| 狂野欧美激情性bbbbbb| 韩国av在线不卡| 亚洲精品第二区| 美女高潮的动态| 中文字幕制服av| 乱码一卡2卡4卡精品| 免费少妇av软件| 国产午夜精品一二区理论片| 人人妻人人看人人澡| 小蜜桃在线观看免费完整版高清| 中国美白少妇内射xxxbb| 日本熟妇午夜| .国产精品久久| 我的老师免费观看完整版| 天天躁夜夜躁狠狠久久av| 亚洲真实伦在线观看| 久久国产乱子免费精品| eeuss影院久久| 草草在线视频免费看| 国产成年人精品一区二区| 国产一区二区亚洲精品在线观看| 深爱激情五月婷婷| 国产成人a∨麻豆精品| 最近手机中文字幕大全| 热99国产精品久久久久久7| 国产精品99久久99久久久不卡 | 天堂俺去俺来也www色官网| 欧美成人精品欧美一级黄| 久久精品国产a三级三级三级| av在线老鸭窝| 国产在线男女| 高清视频免费观看一区二区| 久久久久久久久久久免费av| 三级经典国产精品| 看十八女毛片水多多多| av福利片在线观看| 亚洲国产高清在线一区二区三| 三级国产精品欧美在线观看| 在线观看av片永久免费下载| 日本与韩国留学比较| 亚洲国产色片| 一区二区三区四区激情视频| 国产中年淑女户外野战色| 精品久久久噜噜| 国产精品一区二区三区四区免费观看| 亚洲精品aⅴ在线观看| 日韩免费高清中文字幕av| 91精品伊人久久大香线蕉| 国产高清有码在线观看视频| 国产精品熟女久久久久浪| 国产爽快片一区二区三区| 美女被艹到高潮喷水动态| 在线精品无人区一区二区三 | 日韩免费高清中文字幕av| 国产乱来视频区| 国产又色又爽无遮挡免| 午夜福利高清视频| 日韩av在线免费看完整版不卡| 久久久久久久久久人人人人人人| 亚洲自偷自拍三级| 欧美成人a在线观看| 乱系列少妇在线播放| 国产乱人偷精品视频| 国产成人午夜福利电影在线观看| 亚洲一区二区三区欧美精品 | 啦啦啦啦在线视频资源| 亚洲性久久影院| 久久精品久久久久久噜噜老黄| 欧美3d第一页| 舔av片在线| 日本三级黄在线观看| 国产又色又爽无遮挡免| 精品少妇久久久久久888优播| 伊人久久国产一区二区| 精品久久国产蜜桃| 国精品久久久久久国模美| 男女那种视频在线观看| 国产精品爽爽va在线观看网站| 波野结衣二区三区在线| 欧美少妇被猛烈插入视频| 国产黄色视频一区二区在线观看| 亚洲国产色片| 国产黄频视频在线观看| 国产 一区精品| 亚洲一级一片aⅴ在线观看| 欧美成人午夜免费资源| 国产高潮美女av| av在线播放精品| 麻豆成人午夜福利视频| 欧美三级亚洲精品| a级毛片免费高清观看在线播放| 国产av码专区亚洲av| 丰满乱子伦码专区| 永久网站在线| 精品人妻视频免费看| 乱系列少妇在线播放| 亚洲av福利一区| 黄色配什么色好看| 日韩av免费高清视频| 亚洲av在线观看美女高潮| 欧美区成人在线视频| 日韩国内少妇激情av| 久久久精品94久久精品| 观看免费一级毛片| 老司机影院毛片| 少妇猛男粗大的猛烈进出视频 | 精品久久久久久电影网| 精品国产一区二区三区久久久樱花 | 啦啦啦啦在线视频资源| 精品久久久久久久久亚洲| 一区二区av电影网| 亚洲欧洲日产国产| 午夜亚洲福利在线播放| 大香蕉97超碰在线| 晚上一个人看的免费电影| 成年女人在线观看亚洲视频 | 久久6这里有精品| 国产中年淑女户外野战色| av在线app专区| 97超碰精品成人国产| 精品国产一区二区三区久久久樱花 | 国产精品蜜桃在线观看| 亚洲人成网站高清观看| 国产成人一区二区在线| 成年版毛片免费区| 久久久久久九九精品二区国产| 成人国产麻豆网| 男女国产视频网站| 欧美成人a在线观看| 99热全是精品| 青春草视频在线免费观看| 夫妻性生交免费视频一级片| 免费观看在线日韩| 自拍偷自拍亚洲精品老妇| 色吧在线观看| 观看美女的网站| 18禁裸乳无遮挡免费网站照片| 干丝袜人妻中文字幕| 国产淫语在线视频| 人妻制服诱惑在线中文字幕| 日本黄大片高清| 美女主播在线视频| 午夜激情久久久久久久| 精品一区在线观看国产| 日本av手机在线免费观看| 好男人在线观看高清免费视频| 国产久久久一区二区三区| 中文资源天堂在线| 在线a可以看的网站| 日韩人妻高清精品专区| 一区二区av电影网| 国产高清国产精品国产三级 | 能在线免费看毛片的网站| 51国产日韩欧美| 毛片女人毛片| 精品99又大又爽又粗少妇毛片| 成人午夜精彩视频在线观看| 天天躁夜夜躁狠狠久久av| 亚洲伊人久久精品综合| 美女内射精品一级片tv| 亚洲国产成人一精品久久久| 中文资源天堂在线| 国产成人精品久久久久久| 精品人妻视频免费看| 校园人妻丝袜中文字幕| 热re99久久精品国产66热6| 欧美日本视频| 国产黄色免费在线视频| 久久精品综合一区二区三区| 成年人午夜在线观看视频| 男女那种视频在线观看| 熟女电影av网| 少妇人妻久久综合中文| 丰满少妇做爰视频| videos熟女内射| 在线a可以看的网站| 精品一区二区免费观看| 国产av国产精品国产| 成人二区视频| 亚洲在久久综合| av免费观看日本| 嫩草影院新地址| 国产乱来视频区| 亚洲无线观看免费| 男的添女的下面高潮视频| 亚洲欧洲日产国产| 国产黄色视频一区二区在线观看| 国产成人精品一,二区| 人妻一区二区av| 国产女主播在线喷水免费视频网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 狂野欧美激情性bbbbbb| 欧美潮喷喷水| 日韩av不卡免费在线播放| 又粗又硬又长又爽又黄的视频| 能在线免费看毛片的网站| 国产 一区 欧美 日韩| 夫妻性生交免费视频一级片| 久久久久久久久大av| 在线观看免费高清a一片| 精品99又大又爽又粗少妇毛片| 在线精品无人区一区二区三 | 亚洲av成人精品一二三区| 性插视频无遮挡在线免费观看| 直男gayav资源| 久久鲁丝午夜福利片| av免费观看日本| 成人特级av手机在线观看| 一级黄片播放器| 国产成人a区在线观看| 夫妻性生交免费视频一级片| 久久久久久伊人网av| 国产免费一区二区三区四区乱码| 中文在线观看免费www的网站| 久久久久久久久久人人人人人人| 亚洲欧洲国产日韩| 一级毛片久久久久久久久女| 午夜激情福利司机影院| 亚洲欧美成人精品一区二区| 国产乱人视频| av免费在线看不卡| 久久韩国三级中文字幕| 婷婷色av中文字幕| 边亲边吃奶的免费视频| 婷婷色综合www| 天天一区二区日本电影三级| 91久久精品国产一区二区三区| 国产日韩欧美亚洲二区| 又大又黄又爽视频免费| 亚洲精品456在线播放app| 中文天堂在线官网| 亚洲四区av| 黄色一级大片看看| 欧美性感艳星| 波多野结衣巨乳人妻| 精品久久久久久电影网| 狠狠精品人妻久久久久久综合| 亚洲av中文av极速乱| 欧美日韩国产mv在线观看视频 | 欧美日韩在线观看h| 国产老妇伦熟女老妇高清| 亚洲国产日韩一区二区| 国产综合精华液| 成人一区二区视频在线观看| 我的女老师完整版在线观看| av线在线观看网站| 久久久久精品久久久久真实原创| 免费看日本二区| 在线观看av片永久免费下载| 亚洲精品日本国产第一区| 亚洲成色77777| 久久综合国产亚洲精品| 久久精品熟女亚洲av麻豆精品| 精品久久久久久电影网| 蜜桃久久精品国产亚洲av| 亚洲精品乱码久久久久久按摩| 成人无遮挡网站| 亚洲最大成人中文| 18禁在线播放成人免费| 波多野结衣巨乳人妻| 97超视频在线观看视频| 别揉我奶头 嗯啊视频| 久久99蜜桃精品久久| 国产人妻一区二区三区在| 欧美区成人在线视频| 久久精品久久精品一区二区三区| 免费看不卡的av| 亚洲av福利一区| 尤物成人国产欧美一区二区三区| 亚洲av成人精品一二三区| 国产69精品久久久久777片| 日本黄大片高清| 亚洲高清免费不卡视频| 色5月婷婷丁香| 国产亚洲午夜精品一区二区久久 | 久久久精品欧美日韩精品| 欧美一级a爱片免费观看看| 亚洲国产av新网站| 国产久久久一区二区三区| 国产午夜精品久久久久久一区二区三区| 下体分泌物呈黄色| 国内精品美女久久久久久| 91精品伊人久久大香线蕉| 伦精品一区二区三区| 欧美少妇被猛烈插入视频| 国产成人精品一,二区| 99久久精品热视频| 丝瓜视频免费看黄片| av免费观看日本| 青春草国产在线视频| 我的老师免费观看完整版| 性色avwww在线观看| 白带黄色成豆腐渣| 啦啦啦中文免费视频观看日本| 91精品一卡2卡3卡4卡| 亚洲欧洲日产国产| 国产伦精品一区二区三区四那| 国产一区亚洲一区在线观看| 亚洲人成网站高清观看| 国产欧美另类精品又又久久亚洲欧美| 日本一本二区三区精品| 国产成年人精品一区二区| 亚洲国产精品国产精品| 美女国产视频在线观看| 国产久久久一区二区三区| 五月伊人婷婷丁香| 中文字幕亚洲精品专区| 日韩伦理黄色片| 免费看日本二区| 亚洲精品,欧美精品| 少妇 在线观看| 国内精品美女久久久久久| 亚洲综合精品二区| 欧美三级亚洲精品| 日韩亚洲欧美综合| 成人鲁丝片一二三区免费| 精品国产一区二区三区久久久樱花 | 各种免费的搞黄视频| 国产乱人偷精品视频| 少妇被粗大猛烈的视频| 女的被弄到高潮叫床怎么办| 亚洲自拍偷在线| 三级国产精品片| 欧美日韩一区二区视频在线观看视频在线 | 国产中年淑女户外野战色| 男的添女的下面高潮视频| 国产老妇伦熟女老妇高清| 日本-黄色视频高清免费观看| 六月丁香七月| 国产老妇女一区| 91精品伊人久久大香线蕉| 国产成人精品福利久久| 国产精品久久久久久av不卡| 成人亚洲欧美一区二区av| 免费观看无遮挡的男女| 黄色一级大片看看| 亚洲婷婷狠狠爱综合网| 99视频精品全部免费 在线| 少妇裸体淫交视频免费看高清| 久久亚洲国产成人精品v| 日韩强制内射视频| av在线亚洲专区| 美女xxoo啪啪120秒动态图| 欧美变态另类bdsm刘玥| 国产精品精品国产色婷婷| 亚洲精华国产精华液的使用体验| 噜噜噜噜噜久久久久久91| 欧美97在线视频| 少妇被粗大猛烈的视频| 精品久久久久久久末码| 91精品国产九色| 欧美日韩综合久久久久久| 偷拍熟女少妇极品色| 亚洲无线观看免费| 最近中文字幕高清免费大全6| 搞女人的毛片| 国产黄片视频在线免费观看| 特级一级黄色大片| 最近中文字幕高清免费大全6| 日韩成人伦理影院| 精品少妇久久久久久888优播| 特级一级黄色大片| 亚洲色图综合在线观看| .国产精品久久| videossex国产| 精品国产一区二区三区久久久樱花 | 搡老乐熟女国产| kizo精华| 国产精品女同一区二区软件| 99re6热这里在线精品视频| 青青草视频在线视频观看| 中国美白少妇内射xxxbb| 偷拍熟女少妇极品色| 高清视频免费观看一区二区| 亚洲aⅴ乱码一区二区在线播放| 久久精品国产亚洲av涩爱| 综合色丁香网| 网址你懂的国产日韩在线| 美女内射精品一级片tv| 国产精品久久久久久久电影| 麻豆成人午夜福利视频| 天天躁夜夜躁狠狠久久av| 免费黄色在线免费观看| av在线观看视频网站免费| 青春草国产在线视频| 日韩欧美精品v在线| 干丝袜人妻中文字幕| 日本午夜av视频| 国产91av在线免费观看| 色吧在线观看| 欧美潮喷喷水| 国产精品女同一区二区软件| 女人久久www免费人成看片| 女人十人毛片免费观看3o分钟| 看免费成人av毛片| 亚洲精品影视一区二区三区av| 99久久九九国产精品国产免费| 在线观看免费高清a一片| 久久久久久久久久人人人人人人| 国产精品嫩草影院av在线观看| 国内少妇人妻偷人精品xxx网站| 国产精品熟女久久久久浪| 亚洲av欧美aⅴ国产| 99热6这里只有精品| 亚洲成人精品中文字幕电影| 成人亚洲精品av一区二区| 2021天堂中文幕一二区在线观| 三级国产精品片| 日本午夜av视频| 丝瓜视频免费看黄片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 街头女战士在线观看网站| 亚洲成人一二三区av| 日韩精品有码人妻一区| 乱码一卡2卡4卡精品| 一级片'在线观看视频| 国产欧美日韩精品一区二区| 舔av片在线| 免费看光身美女| 18禁在线无遮挡免费观看视频| 大陆偷拍与自拍| 色婷婷久久久亚洲欧美| 国产亚洲5aaaaa淫片| 晚上一个人看的免费电影| 亚洲怡红院男人天堂| 国产成人精品一,二区| 国国产精品蜜臀av免费| 十八禁网站网址无遮挡 | 欧美性感艳星| 国产黄片美女视频| 一级毛片 在线播放| tube8黄色片| 高清av免费在线| 欧美老熟妇乱子伦牲交| 少妇猛男粗大的猛烈进出视频 | 久久久久性生活片| 国内揄拍国产精品人妻在线| 黄色欧美视频在线观看| 美女xxoo啪啪120秒动态图| 精品酒店卫生间| 婷婷色综合www| 久久热精品热| 欧美日韩亚洲高清精品| 一级毛片 在线播放| av专区在线播放| 在线观看美女被高潮喷水网站| 一个人看的www免费观看视频| 国产亚洲av嫩草精品影院| 久久精品人妻少妇| 久久久久久久国产电影| 在线观看三级黄色| 91狼人影院| 大片免费播放器 马上看| 国内精品美女久久久久久| 日韩一区二区三区影片| 国产精品秋霞免费鲁丝片| 日韩视频在线欧美| 少妇裸体淫交视频免费看高清| 人人妻人人澡人人爽人人夜夜| 99久久精品国产国产毛片| 亚洲av中文av极速乱| 中文字幕制服av| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲综合色惰| 免费人成在线观看视频色| 亚洲av一区综合| 成年版毛片免费区| 亚洲va在线va天堂va国产| 91午夜精品亚洲一区二区三区| 成年免费大片在线观看| 亚洲欧洲国产日韩| 成人欧美大片| 女人十人毛片免费观看3o分钟| 免费观看av网站的网址| 色视频www国产| 久久久成人免费电影| 午夜免费观看性视频| 最近中文字幕高清免费大全6| 在线天堂最新版资源| 99热6这里只有精品| 成人国产麻豆网| 男女无遮挡免费网站观看| 综合色丁香网| 成年女人看的毛片在线观看| 欧美丝袜亚洲另类| 亚洲成人中文字幕在线播放| 免费黄频网站在线观看国产| 欧美性感艳星| 国产精品久久久久久久电影| 天堂俺去俺来也www色官网| 亚洲国产精品专区欧美| av专区在线播放| 久久6这里有精品| 麻豆久久精品国产亚洲av| 青春草国产在线视频| 97人妻精品一区二区三区麻豆| 久久99热6这里只有精品| 交换朋友夫妻互换小说| 日本午夜av视频| 人妻系列 视频| 男女啪啪激烈高潮av片| 色视频在线一区二区三区| 久久鲁丝午夜福利片| 国产精品久久久久久av不卡| 高清av免费在线| 中文字幕久久专区| 欧美国产精品一级二级三级 | 亚洲国产精品成人综合色| 黄色日韩在线| 国产av国产精品国产| 伦理电影大哥的女人| 精品一区在线观看国产| 一级毛片aaaaaa免费看小| 爱豆传媒免费全集在线观看| 亚洲不卡免费看| www.色视频.com| 中文在线观看免费www的网站| 一本色道久久久久久精品综合| 又爽又黄a免费视频| 丰满人妻一区二区三区视频av| 在线观看人妻少妇| 高清av免费在线| 国产午夜精品一二区理论片| 美女国产视频在线观看| 99久久人妻综合| 久久人人爽人人爽人人片va| 中文资源天堂在线| 少妇裸体淫交视频免费看高清| 人妻一区二区av| 视频中文字幕在线观看| 尾随美女入室| 亚洲最大成人av| 99久久九九国产精品国产免费| 久久97久久精品| 99久久九九国产精品国产免费| 日本av手机在线免费观看| 不卡视频在线观看欧美| 91精品伊人久久大香线蕉| 国产免费福利视频在线观看| 婷婷色综合www| 亚洲电影在线观看av| av国产久精品久网站免费入址| 边亲边吃奶的免费视频| 美女cb高潮喷水在线观看| 欧美日韩综合久久久久久| 欧美成人一区二区免费高清观看| 欧美最新免费一区二区三区| 久久97久久精品| 老司机影院毛片| 一本一本综合久久| 免费观看av网站的网址| 黄色视频在线播放观看不卡| 2022亚洲国产成人精品| 美女高潮的动态| 97人妻精品一区二区三区麻豆| 亚洲av在线观看美女高潮| 亚洲自偷自拍三级| 午夜视频国产福利| 91久久精品国产一区二区三区| 人人妻人人看人人澡| 久久精品久久精品一区二区三区| 日韩国内少妇激情av| 男男h啪啪无遮挡| 肉色欧美久久久久久久蜜桃 | 两个人的视频大全免费| 国产精品一区www在线观看| 日韩一本色道免费dvd| 国产精品国产三级国产专区5o| 亚洲欧美日韩无卡精品| 18禁动态无遮挡网站| 国产精品秋霞免费鲁丝片| 中文在线观看免费www的网站| 大码成人一级视频| 精品视频人人做人人爽| 精品一区在线观看国产| 精品酒店卫生间| 99久国产av精品国产电影| 中文天堂在线官网| 97精品久久久久久久久久精品| 成人黄色视频免费在线看| 91精品国产九色| 国产成人91sexporn| 国产色婷婷99| 成人漫画全彩无遮挡| 国产欧美亚洲国产| 日韩欧美精品免费久久| 午夜福利高清视频| 欧美3d第一页| 亚洲伊人久久精品综合| 国产色爽女视频免费观看| 亚洲国产欧美人成| 国产成人freesex在线| 97热精品久久久久久| 亚洲第一区二区三区不卡| 亚洲人与动物交配视频| 欧美日韩在线观看h| av黄色大香蕉| 欧美另类一区| 女人久久www免费人成看片| 久久国内精品自在自线图片| 欧美精品一区二区大全| 建设人人有责人人尽责人人享有的 |