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

    基于影錐特征角的近圓軌道衛(wèi)星進(jìn)出地影時(shí)間的快速判定法

    2020-07-02 09:43:34謝文杰邵圣祥宋易敏
    載人航天 2020年3期
    關(guān)鍵詞:半影陰影軌道

    謝文杰,邵圣祥,王 霞,宋易敏

    (1. 中國人民解放軍63611 部隊(duì),庫爾勒841001; 2. 中國人民解放軍63610 部隊(duì),庫爾勒841001)

    1 引言

    當(dāng)太陽光線受到地球遮擋后會(huì)在地球背后形成陰影區(qū)域,稱為地影。 許多繞地運(yùn)行的衛(wèi)星都要周期性地進(jìn)出地影。 衛(wèi)星位于地影區(qū)時(shí),由于星體受到的太陽輻射大幅減少或消失,不僅影響衛(wèi)星服務(wù)功能的提供[1-3],而且限制航天發(fā)射和測(cè)控任務(wù)的實(shí)施[4-7]。 因此,準(zhǔn)確、快速計(jì)算衛(wèi)星進(jìn)出地影的時(shí)刻及陰影時(shí)間對(duì)制定衛(wèi)星溫控策略、確定太陽輻射壓大小、規(guī)劃航天任務(wù)等至關(guān)重要[8-11]。

    在工程應(yīng)用中主要有2 種地影模型:柱形地影模型和錐形地影模型[12]。 柱形地影模型比較簡(jiǎn)單,只能描述衛(wèi)星被完全照亮或遮擋2 種狀態(tài),適用于低軌衛(wèi)星的陰影預(yù)報(bào)[1],高軌衛(wèi)星建議謹(jǐn)慎采用[2]。 錐形地影模型考慮到太陽被遮掩的程度,構(gòu)造了相對(duì)復(fù)雜和真實(shí)的陰影結(jié)構(gòu),導(dǎo)致陰影判斷條件增多、計(jì)算量加大[13-14],但是卻成為各種高精度應(yīng)用的必然選擇[3]。 目前關(guān)于錐形地影模型陰影判定依據(jù),一種是利用衛(wèi)星的角度得出太陽、地球的視半徑、角距和視面積來判斷[2];另一種是根據(jù)衛(wèi)星到陰影軸的距離與影錐在基本面上的投影半徑進(jìn)行判斷[4]。 由于衛(wèi)星時(shí)刻處于運(yùn)動(dòng)之中,根據(jù)衛(wèi)星位置確定的天體視半徑、投影基本面也隨之不斷改變,所以在迭代判定衛(wèi)星陰影狀態(tài)的過程中需要更新每一時(shí)刻的陰影條件相關(guān)參數(shù)。 現(xiàn)有文獻(xiàn)沒有對(duì)錐形地影模型不同陰影計(jì)算方法的預(yù)報(bào)精度、計(jì)算效率進(jìn)行定量分析。

    本文借鑒賈向華[1]論述的陰影特征角和判定角的概念,把變化率很小的影錐半錐角作為比較的參考基準(zhǔn),提出了一種新的錐形陰影判定條件。 只對(duì)背向太陽側(cè)的衛(wèi)星軌道進(jìn)行地影計(jì)算,并給出了基于時(shí)間窗口數(shù)據(jù)篩選用于近圓軌道衛(wèi)星的地影快速計(jì)算方法;使用STK 和MATLAB 工具對(duì)3 種軌道高度的衛(wèi)星進(jìn)行地影仿真計(jì)算,將不同算法計(jì)算結(jié)果與STK 結(jié)果進(jìn)行比較和分析,以驗(yàn)證本文算法的精度和效率。

    2 錐形地影模型及其判定條件

    將太陽視為有一定大小的發(fā)光球體,忽略地球扁率和大氣的影響,衛(wèi)星為一個(gè)不計(jì)尺寸的點(diǎn)目標(biāo),太陽光以直線形式傳播,這樣因地球遮擋在空間形成的圓錐形陰影如圖1 所示。 在地球的背后形成了太陽光完全照射不到的區(qū)域和受到部分太陽光照射的區(qū)域,前者稱為本影,后者稱為半影。 以日心Os和地心Oe連線為軸,整個(gè)地影由本影圓錐和半影圓錐同軸互相嵌套構(gòu)成。V1為半影圓錐的頂點(diǎn),α為半影錐的半錐角;V2為本影圓錐的頂點(diǎn),β為本影錐的半錐角。 本文結(jié)合圓錐形陰影幾何特征定義陰影特征角和判定角,根據(jù)空間星體幾何關(guān)系推導(dǎo)出它們的求解方法和衛(wèi)星進(jìn)出地影的判定條件。

    圖1 錐形地影模型Fig.1 Conical Earth shadow model

    定義rs、re、rt分別為太陽、地球和衛(wèi)星在慣性坐標(biāo)系(如J2000.0 坐標(biāo)系)中的位置矢量,則太陽、衛(wèi)星相對(duì)于地球的坐標(biāo)s、st分別為式(1)、式(2):

    Rs、Re分別為太陽和地球的半徑,由相似三角形和比例的性質(zhì)可以計(jì)算得到式(3)、式(4):

    半影判定角是影錐頂點(diǎn)V1和衛(wèi)星位置Ot的連線與陰影軸的夾角,角的開口方向與α一致,記作θ;本影判定角是影錐頂點(diǎn)V2和衛(wèi)星位置Ot的連線與陰影軸所夾的非鈍角,記作φ。s1、s2分別為地心Oe指向影錐頂點(diǎn)V1、V2的矢量,φ為輔助變量。 以上夾角和矢量由式(5)~式(9)確定:

    在地球一個(gè)公轉(zhuǎn)周期中,地球與太陽之間的距離在1.47×108km 到1.52×108km 之間變化[8],則α的變化區(qū)間為0.264 8°~0.273 8°,β的變化區(qū)間為0.260 0°~0.268 8°,α的變化率為5.707 8×10-10°/s、β的變化率為5.580 9×10-10°/s,可見α、β隨時(shí)間變化很小,在短期內(nèi)可以視為常值。 由圖1可以看出,地球的本影區(qū)實(shí)際上內(nèi)嵌在半影圓錐空間之中,所以確定衛(wèi)星是否位于地球的陰影內(nèi),首先要判斷衛(wèi)星是否在外層圓錐范圍內(nèi),然后再判斷衛(wèi)星是否在內(nèi)層圓錐區(qū)域內(nèi)。下面給出地影的判定條件,如式(10)所示:

    式中,θ=α或者φ=β時(shí),表示衛(wèi)星正好處于進(jìn)/出半影或者本影位置。

    3 衛(wèi)星地影快速計(jì)算方法

    太陽、地球和衛(wèi)星在空間沿各自的軌道運(yùn)動(dòng),只有在滿足特定的位置關(guān)系時(shí),衛(wèi)星才可能進(jìn)入到地影中。 因此,實(shí)際上僅有一部分衛(wèi)星軌道穿越地球陰影區(qū),顯然對(duì)所有軌道數(shù)據(jù)進(jìn)行計(jì)算不是合理的做法。 本文提出一種數(shù)據(jù)篩選方法,先剔除衛(wèi)星軌道中不可能位于陰影區(qū)的數(shù)據(jù),然后根據(jù)第2 節(jié)提出的錐形陰影判定條件對(duì)剩余的軌道數(shù)據(jù)進(jìn)行計(jì)算,從而大幅度減少數(shù)據(jù)量,以獲得較高的計(jì)算效率。

    太陽、地球和衛(wèi)星的空間幾何關(guān)系可以用太陽地心矢量s和衛(wèi)星地心矢量st之間的夾角σ來描述,如圖2 所示。 衛(wèi)星從地球向陽的一側(cè)朝著地球背陽一側(cè)運(yùn)動(dòng),當(dāng)衛(wèi)星位于位置O′t時(shí)(圖1),σ為90°,衛(wèi)星與半影圓錐母線交于點(diǎn)M,M到地心的距離為d,其計(jì)算公式如式(11)所示:

    圖2 太陽、地球和衛(wèi)星的空間位置關(guān)系Fig.2 Relative position of the Sun,the Earth and the satellite

    圖3 夾角σ 與衛(wèi)星地影的關(guān)系Fig.3 Relationship between angle σ and Earth shadow of the satellite

    把上述星體幾何約束關(guān)系轉(zhuǎn)換為對(duì)應(yīng)的衛(wèi)星軌道數(shù)據(jù)時(shí)間窗口,可以剔除不可能發(fā)生地影的時(shí)段。 如圖2 所示,以地球?yàn)橹行?,太陽在黃道面上運(yùn)動(dòng),衛(wèi)星在其軌道面上運(yùn)動(dòng),太陽地心矢量s在衛(wèi)星軌道面上的投影與衛(wèi)星地心矢量st所成的夾角η和σ的變化趨勢(shì)一致。 根據(jù)幾何中的三垂線定理,當(dāng)σ=90°時(shí),η=90°。 這樣可以用夾角η在同一個(gè)軌道平面中等價(jià)表征太陽、地球和衛(wèi)星的位置關(guān)系,并且方便運(yùn)用天體運(yùn)動(dòng)規(guī)律計(jì)算衛(wèi)星從某一時(shí)刻初始位置到達(dá)特定位置所需經(jīng)歷的時(shí)間。

    為方便計(jì)算,太陽、衛(wèi)星軌道數(shù)據(jù)時(shí)間采用約簡(jiǎn)儒略日表示。 地心衛(wèi)星軌道坐標(biāo)系定義為:坐標(biāo)原點(diǎn)位于地心,x軸在軌道面內(nèi)且指向近地點(diǎn),y軸在軌道面內(nèi)且與x軸垂直,指向由x軸逆時(shí)針旋轉(zhuǎn)90°確定,z軸與x軸、y軸構(gòu)成右手螺旋系。 地心赤道坐標(biāo)系到地心衛(wèi)星軌道坐標(biāo)系的轉(zhuǎn)換如式(12)所示:

    式中Xe為地心赤道坐標(biāo)系矢量,Xo為地心衛(wèi)星軌道坐標(biāo)系矢量,RX、RZ分別為繞x、z軸旋轉(zhuǎn)的算子,T 表示轉(zhuǎn)置矩陣[8],Ω是升交點(diǎn)赤經(jīng),i是傾角,ω是近地點(diǎn)輻角。

    目前大多數(shù)衛(wèi)星都采用近圓軌道[8,15]。 這樣可以近似認(rèn)為衛(wèi)星作勻速圓周運(yùn)動(dòng),同時(shí)把太陽也作勻速運(yùn)動(dòng)處理。 太陽和衛(wèi)星在圓形軌道上的角速度分別如式(13)、(14)所示:

    式中365.25 代表地球公轉(zhuǎn)周期,地球引力系數(shù)GM=398 600.441 5 km3/s2,a為衛(wèi)星軌道半長軸。 在衛(wèi)星軌道坐標(biāo)系中,使用式(15)計(jì)算太陽投影和衛(wèi)星的真近點(diǎn)角ν,以確定某一時(shí)刻太陽、衛(wèi)星在軌道面的相對(duì)位置及其夾角。

    式中x、y分別為太陽或衛(wèi)星在地心衛(wèi)星軌道平面中的位置坐標(biāo)。 假設(shè)衛(wèi)星從某一位置運(yùn)行到另一位置繞地心旋轉(zhuǎn)了Δν弧度,則所需時(shí)間t為式(16):

    當(dāng)衛(wèi)星使用逆行軌道,即與太陽運(yùn)動(dòng)方向相反時(shí),式分母?。?hào)。

    綜上,依據(jù)初始時(shí)刻太陽、衛(wèi)星的相對(duì)位置關(guān)系可以計(jì)算出衛(wèi)星在之后周期運(yùn)動(dòng)中可能會(huì)穿越地影的時(shí)間窗口,并對(duì)這些時(shí)間段內(nèi)的軌道數(shù)據(jù)進(jìn)行陰影判定,此時(shí)陰影判定條件中的第2 個(gè)不等式可以忽略。

    4 仿真計(jì)算與分析

    在適當(dāng)?shù)慕萍僭O(shè)下,基于錐形地影模型,利用衛(wèi)星處的太陽、地球的視半徑和視面積判定衛(wèi)星陰影狀態(tài)的方法為天體視半徑判定法[2,4,8](算法1);利用半影錐、本影錐在基本面上的投影半徑及衛(wèi)星到陰影軸的距離進(jìn)行衛(wèi)星地影解算的方法為影錐基本面投影半徑判定法[4,8](算法2)。 這2 種算法簡(jiǎn)單實(shí)用,輸入高精度太陽和衛(wèi)星星歷,能夠計(jì)算得到滿足工程應(yīng)用要求的解。 定義本文提出的基于影錐特征角的陰影判定條件和快速計(jì)算方法為影錐特征角快速判定法(算法3)。 使用MATLAB R2015b軟件編寫計(jì)算機(jī)程序?qū)πl(wèi)星進(jìn)出地影的時(shí)間進(jìn)行預(yù)報(bào)和分析。

    選擇低軌、中軌和高軌3 顆近圓軌道衛(wèi)星進(jìn)行地影時(shí)間計(jì)算和分析,衛(wèi)星數(shù)據(jù)如表1 所示。首先使用STK11.0.1 計(jì)算2019 年7 月30 日0 時(shí)(UTC)至2019 年8 月2 日0 時(shí)(UTC)期間3 顆衛(wèi)星的地影狀態(tài),然后利用其高精度軌道外推模型獲取該時(shí)間范圍內(nèi)間隔為1 s 的太陽和衛(wèi)星J2000.0 坐標(biāo)系位置矢量數(shù)據(jù)。

    STK(Satellite Tool Kit)是一個(gè)經(jīng)過航天任務(wù)驗(yàn)證的軟件,在仿真分析方面具有極高的準(zhǔn)確性[16]。 表2 為3 種算法計(jì)算得出的衛(wèi)星繞地球飛行一圈,經(jīng)歷一次半影、本影的平均持續(xù)時(shí)間,并以STK 地影計(jì)算結(jié)果為基準(zhǔn),計(jì)算每種算法的相對(duì)誤差。 數(shù)據(jù)表明3 種算法地影時(shí)間計(jì)算精度基本相當(dāng),半影時(shí)間最大誤差5.447%,本影時(shí)間最大誤差0.972%,陰影時(shí)間越短、計(jì)算的相對(duì)誤差也越大。

    表1 衛(wèi)星數(shù)據(jù)Table 1 Satellite data

    表2 衛(wèi)星穿越地影的平均持續(xù)時(shí)間及相對(duì)誤差Table 2 Average duration of satellite through penumbra/umbra and relative errors s

    由圖1 可知,衛(wèi)星在周期運(yùn)動(dòng)中按照半影-本影-半影的順序進(jìn)出地影。 以衛(wèi)星第一次進(jìn)入半影的時(shí)刻和最后一次離開半影的時(shí)刻分別作為衛(wèi)星進(jìn)地影的時(shí)刻和出地影的時(shí)刻,將3 種算法的計(jì)算結(jié)果和STK 預(yù)報(bào)的地影時(shí)刻進(jìn)行比較得到衛(wèi)星進(jìn)出地影時(shí)間差,圖4~圖6 分別是3 顆衛(wèi)星在仿真計(jì)算期間進(jìn)出地影時(shí)間差的變化。 從表2和圖4~圖6 可以看出,算法1 和算法2 對(duì)3 種衛(wèi)星的計(jì)算結(jié)果完全相同;對(duì)于低軌、中軌、高軌衛(wèi)星,算法3與算法1、算法2 預(yù)報(bào)的衛(wèi)星進(jìn)/出地影時(shí)刻的最大時(shí)間差分別為1 s、2 s、4 s,且算法3 預(yù)報(bào)的地影時(shí)刻準(zhǔn)確度稍高于算法1、算法2。

    在處理相同的數(shù)據(jù)量下,比較3 種算法的計(jì)算效率。 計(jì)算機(jī)的軟硬件配置為Intel Core i7-9700 處理器、32 GB 內(nèi)存、Windows 7 Service Pack 1 64 位操作系統(tǒng)。 表3 記錄了輸入3 組時(shí)間間隔為1 s、共計(jì)259 201 點(diǎn)太陽和衛(wèi)星位置矢量數(shù)據(jù)時(shí),MATLAB R2015b 軟件測(cè)量得到的每一種算法的平均完成時(shí)間。 由表3 可知,3 種算法的計(jì)算效率高低依次為:算法3、算法1 和算法2,算法3 比算法1 效率提高25.305%,比算法2 提高52.172%,算法3 具有非常明顯的時(shí)間性能優(yōu)勢(shì)。

    圖4 32060 衛(wèi)星進(jìn)出地影時(shí)間差值變化Fig.4 Changes of time difference of satellite 32060 entry/exit the Earth shadow

    圖5 29061 衛(wèi)星進(jìn)出地影時(shí)間差值變化Fig.5 Changes of time difference of satellite 29061 entry/exit the Earth shadow

    圖6 26880 衛(wèi)星進(jìn)出地影時(shí)間差值變化Fig.6 Changes of time difference of satellite 26880 entry/exit the Earth shadow

    表3 3 種算法計(jì)算時(shí)間的比較Table 3 Comparison of time calculated by three methods s

    由仿真結(jié)果可以看出,算法3 地影判定依據(jù)正確,預(yù)報(bào)時(shí)間窗口設(shè)計(jì)合理,在略優(yōu)于現(xiàn)有其他算法預(yù)報(bào)精度的條件下,具有較高的計(jì)算效率,適用于航天工程中大數(shù)據(jù)量、高強(qiáng)度地影時(shí)間計(jì)算。

    5 結(jié)論

    1)新的錐形地影條件以影錐半錐角為特征角和判定依據(jù),在迭代過程中按時(shí)間段計(jì)算特征角,可大幅減少條件中參數(shù)更新的次數(shù);根據(jù)發(fā)生地影時(shí)太陽、地球和衛(wèi)星存在的幾何邊界條件,確定與這些特定空間位置對(duì)應(yīng)的軌道數(shù)據(jù)的時(shí)間窗口,可大大減少參與地影計(jì)算的數(shù)據(jù)量。

    2)仿真計(jì)算結(jié)果表明,本文建立的地影快速計(jì)算方法預(yù)報(bào)精度略優(yōu)于天體視半徑判定法和影錐基本面投影半徑判定法,并且計(jì)算效率分別提高25.305%和52.172%。

    3)盡管影錐特征角快速判定法以地球作為遮擋物進(jìn)行分析和計(jì)算,但該算法同樣適用于其他天體(如月球)遮擋的情況。 下一步應(yīng)研究橢圓軌道幾何特征和運(yùn)動(dòng)規(guī)律,進(jìn)一步拓寬地影快速計(jì)算方法應(yīng)用范圍。

    猜你喜歡
    半影陰影軌道
    5月天象預(yù)報(bào)
    你來了,草就沒有了陰影
    文苑(2020年11期)2020-11-19 11:45:11
    基于單純形法的TLE軌道確定
    CryoSat提升軌道高度與ICESat-2同步運(yùn)行
    食物相克之說靠譜嗎
    方圓(2020年23期)2020-01-07 00:51:47
    朝美重回“相互羞辱軌道”?
    中國詩歌(2019年6期)2019-11-15 00:26:47
    讓光“驅(qū)走”陰影
    半影月食
    陰影魔怪
    精品久久久久久久人妻蜜臀av| 乱人视频在线观看| 人人妻人人看人人澡| 日韩欧美精品免费久久 | 亚洲无线观看免费| 18禁黄网站禁片午夜丰满| 亚洲国产欧美人成| 久久精品91蜜桃| 婷婷丁香在线五月| 琪琪午夜伦伦电影理论片6080| 九九在线视频观看精品| 国产69精品久久久久777片| 亚洲精品美女久久久久99蜜臀| 一区二区三区高清视频在线| 欧美黄色淫秽网站| 成人无遮挡网站| 97超视频在线观看视频| 欧美日韩瑟瑟在线播放| 午夜激情欧美在线| 久久精品91无色码中文字幕| 成人欧美大片| 内地一区二区视频在线| 91在线精品国自产拍蜜月 | 欧美黑人欧美精品刺激| 欧美日韩一级在线毛片| 丰满乱子伦码专区| 麻豆久久精品国产亚洲av| 午夜日韩欧美国产| 看免费av毛片| 女警被强在线播放| 国产一区二区亚洲精品在线观看| 在线观看免费视频日本深夜| 色av中文字幕| 午夜福利免费观看在线| 精华霜和精华液先用哪个| 久久久久国产精品人妻aⅴ院| 九色国产91popny在线| 欧美3d第一页| 熟女人妻精品中文字幕| 亚洲成人中文字幕在线播放| 天天一区二区日本电影三级| 一卡2卡三卡四卡精品乱码亚洲| 性色avwww在线观看| 中文字幕熟女人妻在线| 日本五十路高清| 欧美三级亚洲精品| 欧美精品啪啪一区二区三区| 天堂影院成人在线观看| 久久草成人影院| 国产成人aa在线观看| 在线观看美女被高潮喷水网站 | 婷婷丁香在线五月| 在线观看美女被高潮喷水网站 | 久久国产精品人妻蜜桃| 亚洲第一欧美日韩一区二区三区| 18禁在线播放成人免费| 19禁男女啪啪无遮挡网站| 岛国在线观看网站| 国产激情偷乱视频一区二区| 欧美黄色淫秽网站| 欧美一区二区亚洲| 成年女人永久免费观看视频| avwww免费| 日本五十路高清| 两性午夜刺激爽爽歪歪视频在线观看| 久久精品91无色码中文字幕| 啦啦啦免费观看视频1| 午夜两性在线视频| 精品久久久久久久人妻蜜臀av| 99热这里只有精品一区| 18禁黄网站禁片午夜丰满| 母亲3免费完整高清在线观看| 久久香蕉国产精品| 美女黄网站色视频| 午夜福利18| 国产精品98久久久久久宅男小说| 国内精品久久久久精免费| 亚洲第一欧美日韩一区二区三区| 免费看十八禁软件| 国产高清视频在线观看网站| 麻豆久久精品国产亚洲av| 手机成人av网站| 亚洲国产欧美网| 成人一区二区视频在线观看| 精品一区二区三区视频在线 | 两人在一起打扑克的视频| 欧美激情在线99| 色综合站精品国产| 久久久成人免费电影| 神马国产精品三级电影在线观看| 午夜精品一区二区三区免费看| 99久久综合精品五月天人人| 琪琪午夜伦伦电影理论片6080| 国产不卡一卡二| 欧美bdsm另类| 亚洲中文字幕日韩| 国产视频一区二区在线看| 国内精品久久久久精免费| 嫁个100分男人电影在线观看| 日韩欧美三级三区| 脱女人内裤的视频| 99国产综合亚洲精品| 亚洲av成人不卡在线观看播放网| av天堂中文字幕网| 好看av亚洲va欧美ⅴa在| 桃色一区二区三区在线观看| 非洲黑人性xxxx精品又粗又长| 久久天躁狠狠躁夜夜2o2o| 欧美性猛交黑人性爽| 一进一出好大好爽视频| 12—13女人毛片做爰片一| 一级毛片女人18水好多| 男女下面进入的视频免费午夜| 国产99白浆流出| 观看免费一级毛片| 亚洲人成伊人成综合网2020| 麻豆国产av国片精品| 性色avwww在线观看| 啪啪无遮挡十八禁网站| 日韩欧美免费精品| 九九热线精品视视频播放| 久久久久久久久大av| 97碰自拍视频| 嫩草影视91久久| 国产亚洲精品综合一区在线观看| 欧美乱妇无乱码| 天堂√8在线中文| 国产精品久久久人人做人人爽| 国产三级中文精品| 国内少妇人妻偷人精品xxx网站| 成人高潮视频无遮挡免费网站| 在线播放无遮挡| 国产精品电影一区二区三区| 两个人看的免费小视频| 美女高潮的动态| 18禁黄网站禁片午夜丰满| 一边摸一边抽搐一进一小说| 欧美性感艳星| АⅤ资源中文在线天堂| 少妇的逼好多水| 少妇人妻一区二区三区视频| 少妇人妻一区二区三区视频| 亚洲精品在线观看二区| 久久久久亚洲av毛片大全| 国产av在哪里看| 精品人妻1区二区| 成人av一区二区三区在线看| 亚洲国产高清在线一区二区三| 国产又黄又爽又无遮挡在线| 国产aⅴ精品一区二区三区波| 51国产日韩欧美| 亚洲五月天丁香| www.www免费av| 乱人视频在线观看| 久久久久久久久久黄片| 少妇丰满av| 五月玫瑰六月丁香| 一级a爱片免费观看的视频| avwww免费| www.色视频.com| 国产精品免费一区二区三区在线| 男女做爰动态图高潮gif福利片| 91麻豆精品激情在线观看国产| 色吧在线观看| АⅤ资源中文在线天堂| 听说在线观看完整版免费高清| 成人特级av手机在线观看| 99精品欧美一区二区三区四区| 97超视频在线观看视频| 国产精品影院久久| 日韩精品青青久久久久久| 波多野结衣高清无吗| 午夜福利在线在线| 男女午夜视频在线观看| 最好的美女福利视频网| 国产真实伦视频高清在线观看 | 丰满人妻一区二区三区视频av | 亚洲国产欧美人成| 老汉色av国产亚洲站长工具| 国产主播在线观看一区二区| 一区福利在线观看| 国产野战对白在线观看| av在线蜜桃| 久久精品人妻少妇| 亚洲欧美日韩无卡精品| 国产av不卡久久| 午夜福利在线观看吧| 真人一进一出gif抽搐免费| 午夜精品久久久久久毛片777| 欧美zozozo另类| 欧美成人性av电影在线观看| 香蕉丝袜av| 香蕉久久夜色| 国产精品国产高清国产av| 三级国产精品欧美在线观看| 亚洲 国产 在线| 91久久精品电影网| 日韩欧美国产一区二区入口| 少妇人妻一区二区三区视频| 欧美在线黄色| 国产精品,欧美在线| 五月玫瑰六月丁香| 淫秽高清视频在线观看| 欧美一区二区亚洲| 亚洲美女视频黄频| 久久婷婷人人爽人人干人人爱| 嫩草影视91久久| 国内少妇人妻偷人精品xxx网站| 亚洲 国产 在线| 999久久久精品免费观看国产| 国产精品99久久99久久久不卡| 激情在线观看视频在线高清| 又粗又爽又猛毛片免费看| 1000部很黄的大片| 午夜影院日韩av| 97超视频在线观看视频| 中文字幕熟女人妻在线| 超碰av人人做人人爽久久 | 一边摸一边抽搐一进一小说| 熟女人妻精品中文字幕| 非洲黑人性xxxx精品又粗又长| 欧美黄色淫秽网站| 尤物成人国产欧美一区二区三区| 男女那种视频在线观看| 人人妻人人看人人澡| 欧美不卡视频在线免费观看| 国产精品1区2区在线观看.| 国产一区二区三区视频了| 国产精品影院久久| 午夜福利欧美成人| 97超视频在线观看视频| or卡值多少钱| 亚洲国产精品成人综合色| 很黄的视频免费| 尤物成人国产欧美一区二区三区| 国产成人a区在线观看| 欧美乱码精品一区二区三区| h日本视频在线播放| 精品不卡国产一区二区三区| 2021天堂中文幕一二区在线观| 亚洲自拍偷在线| 床上黄色一级片| 搡老妇女老女人老熟妇| 1024手机看黄色片| 国产真实乱freesex| 日本三级黄在线观看| 国产精品99久久久久久久久| 老司机深夜福利视频在线观看| 欧美激情久久久久久爽电影| 国产视频一区二区在线看| 成年女人看的毛片在线观看| 18禁黄网站禁片免费观看直播| 97超视频在线观看视频| 亚洲,欧美精品.| 午夜福利免费观看在线| 嫁个100分男人电影在线观看| 一级作爱视频免费观看| 人妻久久中文字幕网| 亚洲在线观看片| 一级黄色大片毛片| 日韩成人在线观看一区二区三区| 老汉色av国产亚洲站长工具| 免费看十八禁软件| 久久久久久久精品吃奶| 黄片大片在线免费观看| 国产精华一区二区三区| 欧美日韩福利视频一区二区| 51国产日韩欧美| aaaaa片日本免费| 一本精品99久久精品77| 国产精品一区二区免费欧美| 真人做人爱边吃奶动态| 日本熟妇午夜| 国产视频一区二区在线看| www日本黄色视频网| 久久久久久国产a免费观看| 日本一二三区视频观看| av片东京热男人的天堂| 国产 一区 欧美 日韩| 欧美乱色亚洲激情| 在线天堂最新版资源| 91麻豆精品激情在线观看国产| 成人午夜高清在线视频| 啪啪无遮挡十八禁网站| av女优亚洲男人天堂| 黄色丝袜av网址大全| 久久久久久久久中文| 三级男女做爰猛烈吃奶摸视频| 精品人妻偷拍中文字幕| 九九久久精品国产亚洲av麻豆| 中文字幕熟女人妻在线| 久久性视频一级片| 波多野结衣巨乳人妻| 午夜免费成人在线视频| 狠狠狠狠99中文字幕| 精品久久久久久成人av| 1000部很黄的大片| 超碰av人人做人人爽久久 | 亚洲五月天丁香| 精品午夜福利视频在线观看一区| 精品国内亚洲2022精品成人| 老熟妇乱子伦视频在线观看| 成人国产一区最新在线观看| 国产av不卡久久| 国产成人影院久久av| 色av中文字幕| 免费观看精品视频网站| 亚洲电影在线观看av| 国产高清激情床上av| 久久6这里有精品| 中文字幕精品亚洲无线码一区| 午夜久久久久精精品| 欧美绝顶高潮抽搐喷水| 1000部很黄的大片| 哪里可以看免费的av片| 欧美成人一区二区免费高清观看| 精品久久久久久成人av| 国产精品亚洲美女久久久| 成人国产一区最新在线观看| 乱人视频在线观看| 国产美女午夜福利| 日韩欧美国产在线观看| 搡老熟女国产l中国老女人| 成人特级黄色片久久久久久久| 婷婷精品国产亚洲av在线| 国产蜜桃级精品一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 国产成年人精品一区二区| 丰满人妻一区二区三区视频av | 国模一区二区三区四区视频| 尤物成人国产欧美一区二区三区| 国产午夜精品论理片| 日韩欧美 国产精品| 精品乱码久久久久久99久播| 久久亚洲真实| 日本在线视频免费播放| 亚洲中文字幕日韩| 午夜福利高清视频| 91麻豆精品激情在线观看国产| 国产精品日韩av在线免费观看| 久9热在线精品视频| 色综合欧美亚洲国产小说| 在线十欧美十亚洲十日本专区| 国产午夜福利久久久久久| 精品乱码久久久久久99久播| 免费看美女性在线毛片视频| 国产精品美女特级片免费视频播放器| 成人欧美大片| 国产午夜福利久久久久久| 日韩欧美三级三区| 给我免费播放毛片高清在线观看| 午夜精品一区二区三区免费看| 欧美三级亚洲精品| 亚洲熟妇熟女久久| 午夜福利视频1000在线观看| 国产综合懂色| 国产午夜福利久久久久久| 99精品欧美一区二区三区四区| 午夜福利在线观看吧| 欧美3d第一页| 国产亚洲精品一区二区www| 欧美又色又爽又黄视频| 香蕉丝袜av| 久9热在线精品视频| 国产探花在线观看一区二区| 欧美日韩综合久久久久久 | 欧美一区二区国产精品久久精品| 99riav亚洲国产免费| 九色国产91popny在线| 日本黄大片高清| 久久99热这里只有精品18| 丝袜美腿在线中文| 中出人妻视频一区二区| www.色视频.com| 成人永久免费在线观看视频| 大型黄色视频在线免费观看| 高清日韩中文字幕在线| 俄罗斯特黄特色一大片| 国产一区二区三区在线臀色熟女| 嫩草影院入口| 法律面前人人平等表现在哪些方面| 国产乱人伦免费视频| 最新中文字幕久久久久| 99在线人妻在线中文字幕| 国产日本99.免费观看| 欧美日韩国产亚洲二区| avwww免费| 日韩高清综合在线| 亚洲欧美日韩高清在线视频| 999久久久精品免费观看国产| 欧美日韩中文字幕国产精品一区二区三区| 中亚洲国语对白在线视频| 国产精品av视频在线免费观看| 亚洲性夜色夜夜综合| 国产视频内射| 精品不卡国产一区二区三区| 国产精品久久久人人做人人爽| 一区二区三区免费毛片| 久久人妻av系列| 国产一级毛片七仙女欲春2| or卡值多少钱| 老司机福利观看| 日本黄色视频三级网站网址| 久久精品影院6| 成年免费大片在线观看| 变态另类成人亚洲欧美熟女| 麻豆国产97在线/欧美| 欧美丝袜亚洲另类 | 精品一区二区三区视频在线 | 一二三四社区在线视频社区8| 免费在线观看亚洲国产| 国产成+人综合+亚洲专区| 亚洲精品成人久久久久久| 一边摸一边抽搐一进一小说| 午夜精品一区二区三区免费看| 日本五十路高清| 观看美女的网站| 国产色婷婷99| 亚洲片人在线观看| 99久久无色码亚洲精品果冻| 国内精品久久久久精免费| 久久国产乱子伦精品免费另类| 亚洲av熟女| 国产色婷婷99| 亚洲成人精品中文字幕电影| 国产精品嫩草影院av在线观看 | 一区福利在线观看| 午夜激情欧美在线| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产欧洲综合997久久,| 国产69精品久久久久777片| 一区二区三区高清视频在线| 亚洲最大成人手机在线| 免费无遮挡裸体视频| 深爱激情五月婷婷| 黑人欧美特级aaaaaa片| 国产蜜桃级精品一区二区三区| 亚洲一区高清亚洲精品| 日韩免费av在线播放| 俺也久久电影网| 国产伦一二天堂av在线观看| 长腿黑丝高跟| 蜜桃久久精品国产亚洲av| 色在线成人网| 日韩有码中文字幕| 熟女电影av网| 欧美激情在线99| 搞女人的毛片| 国产一区二区三区在线臀色熟女| 高清在线国产一区| 午夜福利欧美成人| 最新在线观看一区二区三区| 国产三级中文精品| 亚洲av免费高清在线观看| 午夜两性在线视频| 成人精品一区二区免费| 男人舔奶头视频| 国模一区二区三区四区视频| 欧美日韩精品网址| 19禁男女啪啪无遮挡网站| 色尼玛亚洲综合影院| 少妇丰满av| 欧美一级毛片孕妇| 国产午夜精品久久久久久一区二区三区 | 欧美+日韩+精品| 宅男免费午夜| 中文在线观看免费www的网站| 久久精品国产自在天天线| 给我免费播放毛片高清在线观看| 岛国在线观看网站| or卡值多少钱| 窝窝影院91人妻| 免费高清视频大片| 精品人妻一区二区三区麻豆 | 国内揄拍国产精品人妻在线| 最新在线观看一区二区三区| 神马国产精品三级电影在线观看| 日韩亚洲欧美综合| 国产视频一区二区在线看| 国产精品久久久久久人妻精品电影| 成人永久免费在线观看视频| 舔av片在线| 小蜜桃在线观看免费完整版高清| 一二三四社区在线视频社区8| 日本一本二区三区精品| 宅男免费午夜| 国产 一区 欧美 日韩| 婷婷六月久久综合丁香| 亚洲人成网站在线播| 91av网一区二区| 免费看日本二区| 九九久久精品国产亚洲av麻豆| 免费在线观看日本一区| 99国产综合亚洲精品| 成年版毛片免费区| 亚洲精品久久国产高清桃花| 亚洲专区中文字幕在线| 精品99又大又爽又粗少妇毛片 | 日本 av在线| 精品久久久久久久久久免费视频| 亚洲狠狠婷婷综合久久图片| 免费无遮挡裸体视频| 99在线视频只有这里精品首页| 午夜免费成人在线视频| 国产久久久一区二区三区| 麻豆国产97在线/欧美| 亚洲最大成人中文| 日本一二三区视频观看| 午夜福利在线在线| 国产高清videossex| 久久香蕉国产精品| www日本在线高清视频| 黄片大片在线免费观看| 欧美又色又爽又黄视频| 法律面前人人平等表现在哪些方面| 少妇的逼好多水| 国产精品99久久久久久久久| 久久久久九九精品影院| 最好的美女福利视频网| 99久久成人亚洲精品观看| 国产精品三级大全| 网址你懂的国产日韩在线| 国产亚洲精品久久久com| 神马国产精品三级电影在线观看| 高清在线国产一区| 舔av片在线| 国产精品亚洲一级av第二区| 欧美日韩一级在线毛片| 国产一区在线观看成人免费| 免费人成在线观看视频色| 性色av乱码一区二区三区2| 国产99白浆流出| 日本黄色片子视频| 国产精品免费一区二区三区在线| netflix在线观看网站| 亚洲精品影视一区二区三区av| 国产高清视频在线播放一区| 精品午夜福利视频在线观看一区| 亚洲18禁久久av| 国产亚洲精品一区二区www| 成人无遮挡网站| 又粗又爽又猛毛片免费看| 精品久久久久久成人av| 少妇的逼好多水| 99久久九九国产精品国产免费| 夜夜夜夜夜久久久久| avwww免费| 午夜视频国产福利| 欧美黑人巨大hd| 在线观看av片永久免费下载| 欧美另类亚洲清纯唯美| 19禁男女啪啪无遮挡网站| 亚洲av五月六月丁香网| 在线观看免费午夜福利视频| 国产不卡一卡二| 脱女人内裤的视频| av视频在线观看入口| 麻豆久久精品国产亚洲av| 欧美极品一区二区三区四区| 免费电影在线观看免费观看| 国产熟女xx| 成年人黄色毛片网站| www日本黄色视频网| 怎么达到女性高潮| 日日干狠狠操夜夜爽| 日韩中文字幕欧美一区二区| 日本黄色片子视频| 国产欧美日韩精品亚洲av| 欧美不卡视频在线免费观看| 国产精品美女特级片免费视频播放器| 午夜a级毛片| 看片在线看免费视频| 欧美黑人巨大hd| 99国产极品粉嫩在线观看| 国产高清三级在线| 国产黄色小视频在线观看| 国内精品美女久久久久久| 熟女少妇亚洲综合色aaa.| 成年女人永久免费观看视频| 精品不卡国产一区二区三区| 女人被狂操c到高潮| 99在线人妻在线中文字幕| 国产精品99久久99久久久不卡| 久久久久性生活片| 日本精品一区二区三区蜜桃| 少妇熟女aⅴ在线视频| 99精品久久久久人妻精品| 国产免费av片在线观看野外av| 好看av亚洲va欧美ⅴa在| 国产亚洲精品久久久com| 欧美bdsm另类| 亚洲国产欧美网| 久久久久九九精品影院| 欧美在线黄色| 免费一级毛片在线播放高清视频| 亚洲精品影视一区二区三区av| 国产野战对白在线观看| 亚洲人成网站高清观看| 国产高清有码在线观看视频| 久久久久久九九精品二区国产| 最新在线观看一区二区三区| 国产主播在线观看一区二区| 天堂网av新在线| 极品教师在线免费播放| 欧美绝顶高潮抽搐喷水| 69人妻影院| 国产精品1区2区在线观看.| 亚洲国产色片| 亚洲成人中文字幕在线播放| 久久人妻av系列| 嫩草影院精品99| 欧美+日韩+精品| 欧美成人一区二区免费高清观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产精品久久男人天堂|