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

    利用TLE數(shù)據(jù)判別天宮一號(hào)目標(biāo)飛行器姿態(tài)

    2019-01-10 08:38:48煒,崔文,田鑫,雷園,劉
    載人航天 2018年6期
    關(guān)鍵詞:數(shù)據(jù)服務(wù)天宮機(jī)動(dòng)

    張 煒,崔 文,田 鑫,雷 園,劉 興

    (西安衛(wèi)星測(cè)控中心,西安 710043)

    1 引言

    天宮一號(hào)是我國第一個(gè)目標(biāo)飛行器和空間實(shí)驗(yàn)室,于2011年9月29日發(fā)射入軌,先后與神舟八號(hào)、神舟九號(hào)、神舟十號(hào)成功對(duì)接完成相關(guān)任務(wù)。2016年3月16日天宮一號(hào)目標(biāo)飛行器正式終止數(shù)據(jù)服務(wù),無法繼續(xù)維持正常運(yùn)行軌道,在大氣阻力的影響下軌道不斷衰減,最終于2018年4月2日8時(shí)15分左右再入大氣層。軌道壽命末期天宮一號(hào)重約8 t,是我國非受控再入的第一個(gè)超大型航天器。大質(zhì)量航天器在再入過程中并不能被完全燒毀,仍有10%~40%的殘骸返回地球表面[1],對(duì)地表的生命群體、建筑設(shè)施、生態(tài)環(huán)境等具有很大威脅,因此有必要進(jìn)行精確的再入時(shí)間和落點(diǎn)預(yù)報(bào)[2]。

    天宮一號(hào)的再入過程有以下幾個(gè)特點(diǎn):

    1) 再入過程中飛行器姿態(tài)經(jīng)歷了穩(wěn)定到不穩(wěn)定、緩慢翻滾到快速翻滾的變化;

    2) 再入過程中太陽活動(dòng)基本處于極低水平;

    3) 再入的最后階段地磁環(huán)境處于持續(xù)擾動(dòng)狀態(tài)。

    因此,天宮一號(hào)的再入案例對(duì)于開展軌道動(dòng)力學(xué)和再入預(yù)報(bào)研究具有特殊意義。

    大氣阻力是即將再入航天器所受的最主要非保守?cái)z動(dòng)力,精確的面質(zhì)比及合理的大氣阻力特性建模是準(zhǔn)確計(jì)算大氣阻力加速度、預(yù)報(bào)再入時(shí)間的關(guān)鍵[3]。對(duì)于絕大多數(shù)空間目標(biāo),其準(zhǔn)確外形、質(zhì)量、姿態(tài)和表面材料等都是未知的,分別確定大氣阻力系數(shù)、迎風(fēng)面積和質(zhì)量難度很大,因此通常引入彈道系數(shù)B進(jìn)行統(tǒng)一處理,彈道系數(shù)B的定義為式(1)[4]:

    (1)

    其中,CD為大氣阻力系數(shù),A為迎風(fēng)面積,m為質(zhì)量。

    大部分非受控再入航天器都處于姿態(tài)失控狀態(tài),研究天宮一號(hào)在低軌甚至超低軌道情況下的軌道動(dòng)力學(xué)特性,有必要先判別飛行器姿態(tài),針對(duì)不同研究目的選擇不同姿態(tài)階段的數(shù)據(jù)。傳統(tǒng)的姿態(tài)判別方法是基于目標(biāo)的雷達(dá)特征(包括窄帶雷達(dá)散射截面積、一維距離像、二維/三維雷達(dá)圖像等)[5]或光學(xué)特征[6]進(jìn)行判別,這些判別方法對(duì)觀測(cè)條件具有較高的要求,且一個(gè)地基觀測(cè)弧段往往只有幾分鐘,在航天器翻滾周期大于觀測(cè)弧段時(shí)長(zhǎng)的情況下往往難以準(zhǔn)確判斷飛行器姿態(tài)變化,因此需要其他輔助判別方法。

    本文從姿態(tài)變化引起迎風(fēng)面積變化這一特點(diǎn)出發(fā),提出一種基于彈道系數(shù)的天宮一號(hào)目標(biāo)飛行器姿態(tài)判別方法:使用美國戰(zhàn)略司令部公開發(fā)布的雙行根數(shù)(Two Line Element,TLE)作為數(shù)據(jù)源,首先對(duì)TLE數(shù)據(jù)進(jìn)行篩選,剔除精度較差的異常TLE,然后基于多組TLE數(shù)據(jù)使用最小二乘法計(jì)算彈道系數(shù),對(duì)天宮一號(hào)的飛行姿態(tài)進(jìn)行判別。

    2 TLE數(shù)據(jù)預(yù)處理

    本文選擇美國戰(zhàn)略司令部發(fā)布的TLE[7]作為數(shù)據(jù)源。為了比較數(shù)據(jù)服務(wù)終止前后彈道系數(shù)的變化,選擇2015年1月1日作為數(shù)據(jù)起始時(shí)間。由于數(shù)據(jù)服務(wù)終止前天宮一號(hào)可能發(fā)生軌道機(jī)動(dòng),發(fā)布的TLE可能包括精度差甚至錯(cuò)誤結(jié)果,因此TLE使用前需先進(jìn)行預(yù)處理,檢測(cè)軌道機(jī)動(dòng),并剔除異常TLE。根據(jù)大型空間目標(biāo)TLE根數(shù)精度的特點(diǎn),預(yù)處理按以下幾步進(jìn)行[8]:

    1) 剔除發(fā)布后立即被修正的TLE;

    2) 剔除B*為負(fù)的TLE(B*為TLE中的阻力系數(shù)項(xiàng));

    3) 軌道機(jī)動(dòng)檢測(cè);

    4) 基于平運(yùn)動(dòng)變化剔除異常TLE。

    本節(jié)主要介紹軌道機(jī)動(dòng)檢測(cè)和基于平運(yùn)動(dòng)變化剔除異常TLE的方法。

    低軌航天器的軌道機(jī)動(dòng)一般分為軌道面內(nèi)機(jī)動(dòng)和軌道面外機(jī)動(dòng);軌道面內(nèi)機(jī)動(dòng)會(huì)引起平運(yùn)動(dòng)、半長(zhǎng)軸、能量和(或)偏心率的突然變化,而軌道面外機(jī)動(dòng)改變的則是衛(wèi)星軌道的傾角或升交點(diǎn)赤經(jīng)[9]。比較相鄰兩段數(shù)據(jù)的軌道特征量(半長(zhǎng)軸、能量或傾角等),若特征量出現(xiàn)較大變化,則認(rèn)為兩段數(shù)據(jù)之間發(fā)生軌道機(jī)動(dòng)[10]。圖1展示了基于TLE檢測(cè)軌道機(jī)動(dòng)方法的基本原理,本文將連續(xù)6條TLE數(shù)據(jù)組成數(shù)據(jù)段,使用最小二乘法擬合對(duì)應(yīng)一階多項(xiàng)式,依據(jù)連續(xù)兩段數(shù)據(jù)的多項(xiàng)式變化進(jìn)行軌道機(jī)動(dòng)檢測(cè)。檢測(cè)到2015年1月1日以后天宮一號(hào)共進(jìn)行兩次軌道機(jī)動(dòng),一次發(fā)生在2015年4月10日,另一次發(fā)生在2015年12月16日,均為軌道面內(nèi)軌道機(jī)動(dòng)。

    圖1 軌道機(jī)動(dòng)檢測(cè)的基本原理Fig.1 The basic principle of maneuver detection

    基于平運(yùn)動(dòng)變化剔除異常TLE是軌道機(jī)動(dòng)檢測(cè)方法的簡(jiǎn)化和延伸。使用平運(yùn)動(dòng)n作為特征量,尾部數(shù)據(jù)段長(zhǎng)度仍為6,而頭部數(shù)據(jù)段長(zhǎng)度為1,僅擬合尾部數(shù)據(jù)段的多項(xiàng)式,計(jì)算平運(yùn)動(dòng)的相對(duì)變化量TR如式(2)[8]:

    (2)

    其中,Δnc為頭部數(shù)據(jù)n與外推值之間的差,Δnp為多項(xiàng)式外推至頭部數(shù)據(jù)歷元時(shí)n的變化量。方法如圖2所示。

    圖2 基于平運(yùn)動(dòng)變化的異常TLE剔除方法Fig.2 The outlier filtering method based on mean motion

    本文TR的閾值取為0.6,即認(rèn)為TR大于0.6的TLE是異常數(shù)據(jù)。比較預(yù)處理前后天宮一號(hào)TLE序列的B*,結(jié)果如圖3所示??梢钥闯?,預(yù)處理后B*序列的平滑性有了很大改善,認(rèn)為預(yù)處理后天宮一號(hào)的TLE整體軌道精度有了很大改善。

    圖3 預(yù)處理前后天宮一號(hào)的B*序列Fig.3 The pre and post preprocess of B* of TIANGONG-1

    3 基于TLE的彈道系數(shù)計(jì)算方法

    TLE數(shù)據(jù)并不能直接提供天宮一號(hào)的彈道系數(shù)信息,較為相關(guān)的是B*項(xiàng),B*與彈道系數(shù)B之間的關(guān)系可以近似表示為[11]:B=12.741×B*,但是此結(jié)果與真實(shí)彈道系數(shù)仍存在較大誤差?;赥LE計(jì)算航天器的彈道系數(shù),最簡(jiǎn)單的辦法是根據(jù)兩條TLE的半長(zhǎng)軸的衰減進(jìn)行計(jì)算,即通過調(diào)整彈道系數(shù),使較老TLE預(yù)報(bào)至較新TLE歷元時(shí)預(yù)報(bào)的半長(zhǎng)軸與較新TLE的實(shí)際半長(zhǎng)軸一致。彈道系數(shù)的修正可以表示為式(3)[12]:

    (3)

    其中,ΔaTLE為兩條TLE的半長(zhǎng)軸實(shí)際衰減量,ΔaB為使用初始彈道系數(shù)B0進(jìn)行軌道預(yù)報(bào)的半長(zhǎng)軸衰減量,B′為修正后的彈道系數(shù)。這種方法的問題是TLE精度對(duì)彈道系數(shù)精度影響很大。本文對(duì)此方法進(jìn)行改進(jìn),使用最小二乘方法對(duì)N(N≥3)條TLE進(jìn)行彈道系數(shù)的擬合,使各TLE的實(shí)際半長(zhǎng)軸與計(jì)算半長(zhǎng)軸之間的誤差最小。基本過程如下:

    (4)

    其中,ρ為大氣密度,本文使用MSIS-90大氣模型計(jì)算,V為航天器相對(duì)于大氣的運(yùn)動(dòng)速度,e為偏心率,f為真近點(diǎn)角,ωE為地球自轉(zhuǎn)角速度,i為軌道傾角,μ為引力常數(shù),M為平近點(diǎn)角。

    4) 最后一條TLE計(jì)算結(jié)束后,使用公式(5)計(jì)算彈道系數(shù)修正值;

    (5)

    大氣阻力僅有長(zhǎng)期作用效果,且采用開普勒平根數(shù)作為根數(shù)系統(tǒng),因此進(jìn)行軌道積分時(shí)的積分步長(zhǎng)可取為軌道周期的整數(shù)倍,以提高計(jì)算效率。

    4 姿態(tài)判別結(jié)果

    彈道系數(shù)是空間目標(biāo)的固有屬性,不同空間目標(biāo)的彈道系數(shù)可能會(huì)有較大差異,同一空間目標(biāo)在不同姿態(tài)下彈道系數(shù)也可能不同。

    本文的彈道系數(shù)是基于軌道數(shù)據(jù)的解算結(jié)果,彈道系數(shù)吸收了軌道誤差、大氣環(huán)境建模誤差及航天器姿態(tài)變化等,因此解算的彈道系數(shù)并不是常值。對(duì)于天宮一號(hào)而言,預(yù)處理后的TLE用于解算彈道系數(shù)時(shí)精度已經(jīng)足夠,對(duì)彈道系數(shù)精度影響很小;大氣環(huán)境建模的誤差在大氣環(huán)境平靜時(shí)影響不超過10%,擾動(dòng)期間影響可能超過20%[13],是姿態(tài)穩(wěn)定時(shí)造成彈道系數(shù)波動(dòng)的最主要因素;空間目標(biāo)姿態(tài)穩(wěn)定時(shí)面質(zhì)比不發(fā)生變化,快速翻滾的航天器在翻滾方向不發(fā)生改變的情況下也可以認(rèn)為姿態(tài)對(duì)彈道系數(shù)無影響。天宮一號(hào)外形結(jié)構(gòu)比較復(fù)雜,姿態(tài)變化引起的彈道系數(shù)變化較為明顯。使用前文所述方法對(duì)2015年1月1日起美國戰(zhàn)略司令部公開發(fā)布的TLE進(jìn)行預(yù)處理及彈道系數(shù)計(jì)算,結(jié)果如圖4所示。

    圖4 天宮一號(hào)的彈道系數(shù)結(jié)果Fig.4 The ballistic coefficient estimation results of TIANGONG-1

    其中,綠色虛線為天宮一號(hào)終止數(shù)據(jù)服務(wù)時(shí)間,即2016年3月16日,紅色虛線前后(大約2017年10月8日)天宮一號(hào)的彈道系數(shù)有明顯差異。分段進(jìn)行天宮一號(hào)彈道系數(shù)分析:

    1) 天宮一號(hào)數(shù)據(jù)服務(wù)終止前后

    從圖5可以看到,數(shù)據(jù)服務(wù)終止前彈道系數(shù)的均值約為0.0105 m2/kg,數(shù)據(jù)服務(wù)終止后至紅色虛線所示日期前彈道系數(shù)的均值約為0.0113 m2/kg,前后相差不到8%,差異主要由空間環(huán)境周期性變化引起。數(shù)據(jù)服務(wù)終止前后天宮一號(hào)的彈道系數(shù)變化特點(diǎn)基本一致,即使空間環(huán)境穩(wěn)定也存在周期性變化,變化周期約為5天,初步判斷天宮一號(hào)數(shù)據(jù)服務(wù)終止后姿態(tài)未發(fā)生明顯變化。

    圖5 天宮一號(hào)數(shù)據(jù)服務(wù)終止前后彈道系數(shù)特性Fig.5 The characteristic of ballistic coefficients of TIANGONG-1 before and after data service termination

    2) 2017年10月8日前后

    2017年10月5日至10月13日天宮一號(hào)的彈道系數(shù)如圖6所示,期間天宮一號(hào)的彈道系數(shù)逐漸減小,10月6日前天宮一號(hào)的彈道系數(shù)均值約為0.0113 m2/kg,10月12日后天宮一號(hào)彈道系數(shù)的均值約為0.0069 m2/kg。期間太陽活動(dòng)及地磁環(huán)境均穩(wěn)定,無磁擾、磁暴等異常情況發(fā)生,可以排除彈道系數(shù)變化由大氣環(huán)境建模誤差引起。分析認(rèn)為天宮一號(hào)姿態(tài)變化引起彈道系數(shù)持續(xù)減小,2017年10月6日至10月12日,天宮一號(hào)姿態(tài)逐漸變化。

    圖6 2017年10月5日至10月13日天宮一號(hào)的彈道系數(shù)Fig.6 The ballistic coefficients of TIANGONG-1 between Oct. 5 and Oct. 13, 2017

    圖7 2017年10月12日后天宮一號(hào)的彈道系數(shù)結(jié)果Fig.7 The ballistic coefficients of TIANGONG-1 after Oct. 12, 2017

    3) 2017年10月12日后

    圖7為2017年10月12日后天宮一號(hào)的彈道系數(shù)結(jié)果及地磁指數(shù)情況。可以看出,彈道系數(shù)與地磁環(huán)境具有明顯相關(guān)性,可以認(rèn)為地磁變化是造成彈道系數(shù)變化的主要原因,這是因?yàn)楝F(xiàn)有的大氣密度模型不夠精確,地磁環(huán)境擾動(dòng)情況下的大氣密度計(jì)算誤差大于平靜時(shí)期的大氣密度計(jì)算誤差。期間天宮一號(hào)的翻滾速度可能繼續(xù)加快,但由于翻滾方向不變的情況下平均等效迎風(fēng)面積不變,無法再基于彈道系數(shù)判斷翻滾速度變化情況。

    綜上,可以將天宮一號(hào)的姿態(tài)變化分為以下幾個(gè)階段:

    1) 穩(wěn)定階段:2016年3月16日至2017年10月6日。天宮一號(hào)的彈道系數(shù)變化特點(diǎn)與終止數(shù)據(jù)服務(wù)前基本一致,認(rèn)為天宮一號(hào)處于穩(wěn)定狀態(tài);

    2) 逐漸變化階段:2017年10月6日至2017年10月12日。天宮一號(hào)的姿態(tài)從穩(wěn)定轉(zhuǎn)為緩慢翻滾,彈道系數(shù)逐漸減??;

    3) 翻滾階段:2017年10月12日以后。天宮一號(hào)處于持續(xù)翻滾狀態(tài),翻滾速度變化對(duì)彈道系數(shù)的影響很小,彈道系數(shù)基本穩(wěn)定。

    5 結(jié)論

    本文從姿態(tài)變化引起迎風(fēng)面積變化這一特點(diǎn)出發(fā),使用美國戰(zhàn)略司令部公開發(fā)布的雙行根數(shù)TLE作為數(shù)據(jù)源,提出基于彈道系數(shù)的天宮一號(hào)姿態(tài)判別方法,判別出天宮一號(hào)姿態(tài)穩(wěn)定、逐漸變化及翻滾三個(gè)階段。實(shí)測(cè)的雷達(dá)特征數(shù)據(jù)發(fā)現(xiàn),天宮一號(hào)于2017年10月中旬開始有明顯翻滾現(xiàn)象,結(jié)合本文的分析結(jié)果,可以將天宮一號(hào)姿態(tài)變化的時(shí)間精確到10月6日至10月12日之間。本文提出的方法可以作為傳統(tǒng)姿態(tài)判別方法的有效補(bǔ)充。

    值得指出的是,本文提出的方法也存在一定局限性。在飛行器結(jié)構(gòu)較為簡(jiǎn)單、姿態(tài)失穩(wěn)后彈道系數(shù)變化不明顯的情況下,可能無法進(jìn)行準(zhǔn)確判斷。

    猜你喜歡
    數(shù)據(jù)服務(wù)天宮機(jī)動(dòng)
    地理空間大數(shù)據(jù)服務(wù)自然資源調(diào)查監(jiān)測(cè)的方向分析
    天宮出差樂趣多
    天宮之眼
    裝載機(jī)動(dòng)臂的疲勞壽命計(jì)算
    12萬畝機(jī)動(dòng)地不再“流浪”
    機(jī)動(dòng)三輪車的昨天、今天和明天
    如何運(yùn)用稅收大數(shù)據(jù)服務(wù)供給側(cè)結(jié)構(gòu)性改革
    中國商論(2016年34期)2017-01-15 14:24:18
    基于頻繁子圖挖掘的數(shù)據(jù)服務(wù)Mashup推薦
    天宮二號(hào)蓄勢(shì)待發(fā)
    太空探索(2016年6期)2016-07-10 12:09:06
    天宮二號(hào)發(fā)射成功
    大社會(huì)(2016年8期)2016-05-04 03:42:28
    中文字幕免费在线视频6| 国产精品久久电影中文字幕| 一级av片app| 国产探花在线观看一区二区| 国产亚洲精品久久久久久毛片| 国产不卡一卡二| 极品教师在线免费播放| 久久九九热精品免费| 特大巨黑吊av在线直播| 久久国产精品人妻蜜桃| 国产精品人妻久久久久久| 人人妻人人看人人澡| 老司机深夜福利视频在线观看| 亚洲国产精品合色在线| 国产黄片美女视频| 欧美日韩综合久久久久久 | 中文字幕人妻熟人妻熟丝袜美| 久久国产精品人妻蜜桃| 99久久精品一区二区三区| 99久久99久久久精品蜜桃| 亚洲国产高清在线一区二区三| 久久亚洲精品不卡| 成人无遮挡网站| 成人亚洲精品av一区二区| 国产午夜精品论理片| 久久精品人妻少妇| 嫩草影院精品99| 可以在线观看的亚洲视频| 99热这里只有是精品50| 色哟哟哟哟哟哟| 男女床上黄色一级片免费看| 男人的好看免费观看在线视频| 亚洲不卡免费看| 欧美日韩瑟瑟在线播放| 热99在线观看视频| 十八禁人妻一区二区| 内地一区二区视频在线| 亚洲精品一卡2卡三卡4卡5卡| 91字幕亚洲| 亚洲精品成人久久久久久| eeuss影院久久| 久久精品夜夜夜夜夜久久蜜豆| 欧美午夜高清在线| 国产国拍精品亚洲av在线观看| 久久久久久国产a免费观看| 少妇裸体淫交视频免费看高清| 亚洲七黄色美女视频| 搡老妇女老女人老熟妇| 免费人成视频x8x8入口观看| 在线观看午夜福利视频| 99精品久久久久人妻精品| 男人狂女人下面高潮的视频| 精品久久久久久成人av| 国产av一区在线观看免费| 亚洲人成网站高清观看| 免费一级毛片在线播放高清视频| 欧美日本视频| 成人国产综合亚洲| 精品不卡国产一区二区三区| 午夜福利高清视频| 搡老岳熟女国产| 亚洲激情在线av| 有码 亚洲区| 日日干狠狠操夜夜爽| 97热精品久久久久久| 极品教师在线视频| av天堂中文字幕网| 深夜精品福利| 国产成人av教育| 亚洲中文日韩欧美视频| 亚洲精品成人久久久久久| av福利片在线观看| 精品人妻偷拍中文字幕| 亚洲成人中文字幕在线播放| 美女大奶头视频| 欧美午夜高清在线| 久久久久久大精品| 精品久久久久久久末码| 一进一出抽搐动态| 国产v大片淫在线免费观看| av在线老鸭窝| 久久久久久久久中文| 少妇高潮的动态图| 日本黄大片高清| 村上凉子中文字幕在线| 在线观看美女被高潮喷水网站 | av福利片在线观看| 欧美黄色片欧美黄色片| 久久午夜亚洲精品久久| 91在线观看av| 2021天堂中文幕一二区在线观| 长腿黑丝高跟| 亚洲精华国产精华精| 久久久国产成人免费| 毛片女人毛片| 亚洲av电影不卡..在线观看| 国产色婷婷99| 亚洲va日本ⅴa欧美va伊人久久| 国产亚洲精品久久久久久毛片| 精品人妻偷拍中文字幕| 精品久久久久久久久亚洲 | 我的女老师完整版在线观看| 国产精品三级大全| 在线观看美女被高潮喷水网站 | 午夜精品一区二区三区免费看| 成人永久免费在线观看视频| 黄色日韩在线| 国产乱人视频| 成人午夜高清在线视频| 观看免费一级毛片| 麻豆久久精品国产亚洲av| 亚洲精品粉嫩美女一区| 久久草成人影院| 国产精品98久久久久久宅男小说| 亚洲熟妇熟女久久| 精品午夜福利视频在线观看一区| 国语自产精品视频在线第100页| 国产私拍福利视频在线观看| 亚洲成a人片在线一区二区| 亚洲欧美日韩东京热| 草草在线视频免费看| 性色avwww在线观看| 夜夜躁狠狠躁天天躁| 午夜福利高清视频| 91久久精品国产一区二区成人| 人妻久久中文字幕网| 9191精品国产免费久久| 精品国产亚洲在线| 村上凉子中文字幕在线| 欧美黄色片欧美黄色片| 狠狠狠狠99中文字幕| 性色avwww在线观看| 久久精品国产99精品国产亚洲性色| 亚洲在线观看片| 国产精品国产高清国产av| 此物有八面人人有两片| 国产精品一及| 99国产综合亚洲精品| a在线观看视频网站| 蜜桃亚洲精品一区二区三区| 日本一本二区三区精品| 757午夜福利合集在线观看| 一进一出抽搐gif免费好疼| 国产蜜桃级精品一区二区三区| 淫秽高清视频在线观看| 禁无遮挡网站| 久久精品国产亚洲av天美| 搡老熟女国产l中国老女人| 亚洲成人中文字幕在线播放| 99精品久久久久人妻精品| 91狼人影院| 精品人妻视频免费看| 国产精品久久视频播放| 在线观看舔阴道视频| 久久这里只有精品中国| 亚洲人成电影免费在线| 欧美激情久久久久久爽电影| 国产一区二区激情短视频| 美女高潮的动态| 欧美日韩瑟瑟在线播放| 一a级毛片在线观看| 久久久精品欧美日韩精品| 成人av在线播放网站| 亚洲人成网站在线播| 精品久久久久久久人妻蜜臀av| 一级黄色大片毛片| 亚洲自拍偷在线| 亚洲无线在线观看| 午夜免费激情av| 亚洲男人的天堂狠狠| 久久久久久久久中文| 欧美精品啪啪一区二区三区| 成人av一区二区三区在线看| 国内揄拍国产精品人妻在线| 久久久久久久午夜电影| 高清毛片免费观看视频网站| 欧美黄色淫秽网站| 脱女人内裤的视频| 国产精品98久久久久久宅男小说| 中文字幕人成人乱码亚洲影| 亚洲av第一区精品v没综合| 久99久视频精品免费| 成人美女网站在线观看视频| 国产视频内射| 精品久久久久久久久久免费视频| 国产老妇女一区| 国产精品影院久久| 91麻豆精品激情在线观看国产| 中国美女看黄片| 日本一本二区三区精品| 在线观看av片永久免费下载| 我要看日韩黄色一级片| 校园春色视频在线观看| 亚洲成人免费电影在线观看| 国模一区二区三区四区视频| 男人狂女人下面高潮的视频| 国产免费男女视频| 97人妻精品一区二区三区麻豆| 99国产综合亚洲精品| 免费看日本二区| 黄色一级大片看看| 久久这里只有精品中国| 日韩高清综合在线| 国产成人欧美在线观看| 国产午夜精品论理片| av中文乱码字幕在线| 国产成人福利小说| 成人av在线播放网站| 欧美最黄视频在线播放免费| 欧美三级亚洲精品| 黄色女人牲交| 一个人免费在线观看电影| 久久久久性生活片| netflix在线观看网站| 精品国产三级普通话版| 免费在线观看成人毛片| 在线观看免费视频日本深夜| 精品免费久久久久久久清纯| 国产亚洲精品综合一区在线观看| 99在线人妻在线中文字幕| 亚洲av电影不卡..在线观看| 久久久久亚洲av毛片大全| 欧美成人a在线观看| 18+在线观看网站| 国产蜜桃级精品一区二区三区| 亚洲欧美激情综合另类| 亚洲成人久久性| 两个人视频免费观看高清| 国产色婷婷99| 国产精品一及| 精品久久久久久成人av| 成人三级黄色视频| x7x7x7水蜜桃| 国内少妇人妻偷人精品xxx网站| 男女视频在线观看网站免费| 国产免费男女视频| 亚洲狠狠婷婷综合久久图片| 啦啦啦观看免费观看视频高清| 欧美bdsm另类| 色视频www国产| 身体一侧抽搐| 成人无遮挡网站| 一本综合久久免费| 90打野战视频偷拍视频| 99热这里只有是精品50| 在线观看一区二区三区| 欧美中文日本在线观看视频| 精品人妻1区二区| 午夜免费激情av| 在线十欧美十亚洲十日本专区| 国产一区二区激情短视频| 99久久成人亚洲精品观看| 午夜影院日韩av| 午夜精品在线福利| av天堂中文字幕网| 国产精品久久久久久人妻精品电影| 国产av麻豆久久久久久久| 欧美午夜高清在线| 日本三级黄在线观看| 日本与韩国留学比较| 99在线视频只有这里精品首页| 黄色一级大片看看| 美女xxoo啪啪120秒动态图 | 黄片小视频在线播放| 国产单亲对白刺激| 午夜影院日韩av| 99精品久久久久人妻精品| 日本成人三级电影网站| 制服丝袜大香蕉在线| 天堂av国产一区二区熟女人妻| 高清日韩中文字幕在线| 99热这里只有精品一区| 婷婷色综合大香蕉| 最新在线观看一区二区三区| 全区人妻精品视频| 久久天躁狠狠躁夜夜2o2o| 久久久久免费精品人妻一区二区| 亚洲av日韩精品久久久久久密| 乱码一卡2卡4卡精品| 午夜福利欧美成人| 老司机福利观看| 亚洲真实伦在线观看| 亚洲最大成人手机在线| 久久精品91蜜桃| 毛片一级片免费看久久久久 | .国产精品久久| 天美传媒精品一区二区| 久久久久免费精品人妻一区二区| 色哟哟·www| 一本精品99久久精品77| 黄色女人牲交| 国产精品永久免费网站| 一a级毛片在线观看| 美女高潮的动态| 91午夜精品亚洲一区二区三区 | 欧美潮喷喷水| 日本精品一区二区三区蜜桃| 中文在线观看免费www的网站| 色噜噜av男人的天堂激情| 欧美丝袜亚洲另类 | 国产国拍精品亚洲av在线观看| 51国产日韩欧美| 少妇熟女aⅴ在线视频| 亚洲最大成人av| 国产精品伦人一区二区| 三级国产精品欧美在线观看| 好男人在线观看高清免费视频| 国产精品亚洲一级av第二区| 久久天躁狠狠躁夜夜2o2o| 黄色一级大片看看| 一夜夜www| 国产熟女xx| 精品久久久久久久久av| 国产精品野战在线观看| 久久久久久久午夜电影| 啦啦啦观看免费观看视频高清| 国产精品,欧美在线| 久久精品国产亚洲av天美| av国产免费在线观看| 麻豆av噜噜一区二区三区| 此物有八面人人有两片| 亚洲人成伊人成综合网2020| 人人妻,人人澡人人爽秒播| 熟女人妻精品中文字幕| 又粗又爽又猛毛片免费看| av女优亚洲男人天堂| 999久久久精品免费观看国产| 色精品久久人妻99蜜桃| 在线观看免费视频日本深夜| 一本久久中文字幕| 成年女人看的毛片在线观看| 青草久久国产| 1000部很黄的大片| 男人的好看免费观看在线视频| 波野结衣二区三区在线| 亚洲五月天丁香| 白带黄色成豆腐渣| 一个人观看的视频www高清免费观看| 国产伦在线观看视频一区| 色综合站精品国产| 中文字幕高清在线视频| 国产一区二区三区在线臀色熟女| 午夜日韩欧美国产| 亚洲狠狠婷婷综合久久图片| 一本综合久久免费| 嫩草影院新地址| 久久久久免费精品人妻一区二区| 日本免费a在线| 中文字幕av在线有码专区| 国产亚洲欧美在线一区二区| 波多野结衣巨乳人妻| 国产一区二区亚洲精品在线观看| 天堂av国产一区二区熟女人妻| 亚洲精品粉嫩美女一区| 午夜免费男女啪啪视频观看 | 国产精品免费一区二区三区在线| 日本一本二区三区精品| 九色成人免费人妻av| 亚洲欧美日韩高清在线视频| АⅤ资源中文在线天堂| 赤兔流量卡办理| 美女cb高潮喷水在线观看| 一级黄色大片毛片| 青草久久国产| 国产精品综合久久久久久久免费| 熟女人妻精品中文字幕| 国产男靠女视频免费网站| 三级毛片av免费| 波多野结衣高清无吗| 51国产日韩欧美| 97碰自拍视频| 桃色一区二区三区在线观看| 日韩中字成人| 久久伊人香网站| 国内揄拍国产精品人妻在线| 午夜久久久久精精品| 美女黄网站色视频| 99精品在免费线老司机午夜| 长腿黑丝高跟| 中文字幕av在线有码专区| 噜噜噜噜噜久久久久久91| 国产精品不卡视频一区二区 | 国产精品亚洲一级av第二区| 男人和女人高潮做爰伦理| 国内精品美女久久久久久| 久久久国产成人免费| 91久久精品电影网| 99久久精品国产亚洲精品| 日本五十路高清| 欧美3d第一页| 久久亚洲精品不卡| 九色国产91popny在线| avwww免费| www.熟女人妻精品国产| 国产精品一区二区三区四区久久| 啦啦啦韩国在线观看视频| 国产一区二区亚洲精品在线观看| 久久久久免费精品人妻一区二区| 亚洲av熟女| 黄色配什么色好看| 99riav亚洲国产免费| 午夜福利在线观看免费完整高清在 | 一级av片app| 中文字幕人成人乱码亚洲影| 日韩欧美国产在线观看| 九色成人免费人妻av| 小蜜桃在线观看免费完整版高清| 午夜亚洲福利在线播放| .国产精品久久| 久久性视频一级片| 欧美一区二区精品小视频在线| 简卡轻食公司| 亚洲美女黄片视频| 欧美日韩亚洲国产一区二区在线观看| 亚洲av.av天堂| 最近最新中文字幕大全电影3| 国产高清视频在线观看网站| 99久久99久久久精品蜜桃| 观看美女的网站| 少妇熟女aⅴ在线视频| 嫩草影院入口| 亚洲精品在线美女| bbb黄色大片| 99国产极品粉嫩在线观看| 免费看日本二区| 草草在线视频免费看| 国产色爽女视频免费观看| 欧美黑人欧美精品刺激| 看片在线看免费视频| 亚洲av电影在线进入| 日韩国内少妇激情av| 嫩草影院入口| 悠悠久久av| 少妇人妻精品综合一区二区 | 亚洲人与动物交配视频| 美女大奶头视频| 直男gayav资源| 欧美黑人欧美精品刺激| 在线国产一区二区在线| 久久久久久久精品吃奶| 中文亚洲av片在线观看爽| 真人一进一出gif抽搐免费| 97碰自拍视频| 一卡2卡三卡四卡精品乱码亚洲| 国产精品久久久久久精品电影| 国产成人a区在线观看| 欧美一区二区国产精品久久精品| 欧美激情国产日韩精品一区| 亚洲五月天丁香| 日本黄大片高清| x7x7x7水蜜桃| 久久亚洲精品不卡| 人妻丰满熟妇av一区二区三区| a级毛片a级免费在线| 亚洲最大成人中文| 日本一二三区视频观看| 久久精品人妻少妇| 色av中文字幕| 亚洲,欧美精品.| 欧美潮喷喷水| 嫩草影院新地址| 国产精品亚洲一级av第二区| 51国产日韩欧美| www日本黄色视频网| 国产aⅴ精品一区二区三区波| 亚洲av成人av| 欧美+日韩+精品| 久久天躁狠狠躁夜夜2o2o| 免费在线观看影片大全网站| 在线看三级毛片| 国产不卡一卡二| 久久久久久久亚洲中文字幕 | 伦理电影大哥的女人| 国产亚洲av嫩草精品影院| 亚洲最大成人av| 久久久久久久久久成人| 中亚洲国语对白在线视频| 免费在线观看影片大全网站| 一级av片app| 我的女老师完整版在线观看| 欧美高清成人免费视频www| 国产伦人伦偷精品视频| 亚洲国产高清在线一区二区三| 在线免费观看不下载黄p国产 | 狠狠狠狠99中文字幕| 一区二区三区激情视频| 午夜激情欧美在线| 大型黄色视频在线免费观看| 欧美xxxx性猛交bbbb| 自拍偷自拍亚洲精品老妇| 国产成人啪精品午夜网站| av欧美777| 精品日产1卡2卡| 国产91精品成人一区二区三区| 小蜜桃在线观看免费完整版高清| 99久久无色码亚洲精品果冻| 欧美一区二区精品小视频在线| 蜜桃亚洲精品一区二区三区| 亚洲精品色激情综合| 在线国产一区二区在线| 亚洲精华国产精华精| 国内少妇人妻偷人精品xxx网站| 蜜桃亚洲精品一区二区三区| 亚洲18禁久久av| 韩国av一区二区三区四区| 一区二区三区免费毛片| 久久欧美精品欧美久久欧美| 欧洲精品卡2卡3卡4卡5卡区| 成年女人永久免费观看视频| 欧美日韩瑟瑟在线播放| 久久久国产成人免费| 97人妻精品一区二区三区麻豆| 又黄又爽又刺激的免费视频.| 乱人视频在线观看| 欧美区成人在线视频| 很黄的视频免费| 精品国产三级普通话版| 麻豆一二三区av精品| 国产精品久久久久久人妻精品电影| 国产免费一级a男人的天堂| 麻豆av噜噜一区二区三区| 宅男免费午夜| 国产精品亚洲美女久久久| 男插女下体视频免费在线播放| 此物有八面人人有两片| 99久久成人亚洲精品观看| 综合色av麻豆| 婷婷精品国产亚洲av| 国产人妻一区二区三区在| 搡老岳熟女国产| 夜夜躁狠狠躁天天躁| а√天堂www在线а√下载| 三级国产精品欧美在线观看| 特大巨黑吊av在线直播| 精品久久久久久久久av| 国语自产精品视频在线第100页| 熟女人妻精品中文字幕| 草草在线视频免费看| 少妇的逼水好多| 亚洲五月天丁香| 久久欧美精品欧美久久欧美| 亚洲天堂国产精品一区在线| 我要看日韩黄色一级片| 成人特级黄色片久久久久久久| 亚洲欧美激情综合另类| 宅男免费午夜| 成人永久免费在线观看视频| 男人和女人高潮做爰伦理| 婷婷色综合大香蕉| 国产伦精品一区二区三区四那| 我要搜黄色片| 两个人的视频大全免费| 内地一区二区视频在线| 国产精品爽爽va在线观看网站| 在线观看免费视频日本深夜| 首页视频小说图片口味搜索| 亚洲精品色激情综合| 少妇熟女aⅴ在线视频| 国产真实乱freesex| 热99在线观看视频| 亚洲午夜理论影院| 日韩高清综合在线| 亚洲国产精品sss在线观看| 夜夜躁狠狠躁天天躁| 成人午夜高清在线视频| 久久精品人妻少妇| 国产综合懂色| 69人妻影院| 国产免费一级a男人的天堂| 中文字幕免费在线视频6| 丰满乱子伦码专区| 久久久久久久久大av| 国产色婷婷99| 毛片一级片免费看久久久久 | 日韩欧美一区二区三区在线观看| 日日摸夜夜添夜夜添av毛片 | 亚洲三级黄色毛片| 免费搜索国产男女视频| 亚洲精品在线观看二区| 国产aⅴ精品一区二区三区波| 一夜夜www| 久久热精品热| 国产综合懂色| 国产高清有码在线观看视频| 亚洲精品亚洲一区二区| 日韩欧美国产在线观看| 精品免费久久久久久久清纯| 日韩免费av在线播放| 级片在线观看| 18禁在线播放成人免费| 色哟哟哟哟哟哟| 亚洲国产精品久久男人天堂| 深爱激情五月婷婷| 三级国产精品欧美在线观看| 日本三级黄在线观看| 免费看美女性在线毛片视频| 每晚都被弄得嗷嗷叫到高潮| 丰满人妻熟妇乱又伦精品不卡| 超碰av人人做人人爽久久| 国产一区二区三区视频了| 亚洲美女搞黄在线观看 | 高清日韩中文字幕在线| 国产精品乱码一区二三区的特点| 99精品久久久久人妻精品| 欧美午夜高清在线| 在线十欧美十亚洲十日本专区| 美女 人体艺术 gogo| 国产成年人精品一区二区| 国产乱人伦免费视频| 麻豆av噜噜一区二区三区| 精品久久久久久久久av| 中文亚洲av片在线观看爽| 色av中文字幕| 我要看日韩黄色一级片| 每晚都被弄得嗷嗷叫到高潮|