• <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
    日韩大片免费观看网站| 涩涩av久久男人的天堂| 99国产精品免费福利视频| 老鸭窝网址在线观看| 国产成人一区二区三区免费视频网站| 男男h啪啪无遮挡| 国产无遮挡羞羞视频在线观看| 精品亚洲成国产av| 欧美激情久久久久久爽电影 | 日本av手机在线免费观看| 男女高潮啪啪啪动态图| 亚洲综合色网址| 日韩中文字幕视频在线看片| 三级毛片av免费| 麻豆国产av国片精品| 日本一区二区免费在线视频| 久9热在线精品视频| 色婷婷av一区二区三区视频| 99香蕉大伊视频| 操美女的视频在线观看| 亚洲精品国产色婷婷电影| 亚洲一区中文字幕在线| 午夜福利一区二区在线看| 精品国产乱码久久久久久男人| 视频区图区小说| 看免费av毛片| 999久久久国产精品视频| 国产不卡av网站在线观看| 91成年电影在线观看| 欧美av亚洲av综合av国产av| 在线av久久热| 老司机福利观看| 淫妇啪啪啪对白视频 | 一级毛片精品| 免费在线观看日本一区| 正在播放国产对白刺激| 精品第一国产精品| 高清在线国产一区| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品久久午夜乱码| 国产高清国产精品国产三级| 国产日韩一区二区三区精品不卡| 啦啦啦中文免费视频观看日本| 亚洲中文日韩欧美视频| 国产高清视频在线播放一区 | 菩萨蛮人人尽说江南好唐韦庄| 国产av国产精品国产| 欧美 日韩 精品 国产| 亚洲久久久国产精品| 亚洲人成电影观看| av在线播放精品| 亚洲精品中文字幕在线视频| 成人免费观看视频高清| 日本一区二区免费在线视频| 国产成+人综合+亚洲专区| 极品少妇高潮喷水抽搐| 满18在线观看网站| 午夜激情久久久久久久| 国产精品一区二区精品视频观看| 免费在线观看完整版高清| 国产成人影院久久av| √禁漫天堂资源中文www| 精品福利观看| 人妻一区二区av| 国产黄频视频在线观看| 日日摸夜夜添夜夜添小说| 黄色 视频免费看| 国产欧美日韩精品亚洲av| 成人手机av| 精品第一国产精品| 91字幕亚洲| 女人久久www免费人成看片| 亚洲人成电影免费在线| 精品高清国产在线一区| 欧美精品高潮呻吟av久久| 狠狠婷婷综合久久久久久88av| 亚洲精品乱久久久久久| 亚洲国产欧美网| 黄色a级毛片大全视频| 久久亚洲精品不卡| 男人添女人高潮全过程视频| 日韩视频在线欧美| 国产一卡二卡三卡精品| 99国产精品免费福利视频| 国内毛片毛片毛片毛片毛片| 纵有疾风起免费观看全集完整版| 亚洲av欧美aⅴ国产| tocl精华| 最新在线观看一区二区三区| 久久久国产一区二区| 久久久精品94久久精品| 国产99久久九九免费精品| 不卡av一区二区三区| 欧美日韩黄片免| 久久精品熟女亚洲av麻豆精品| 一级,二级,三级黄色视频| 国产日韩一区二区三区精品不卡| 中文字幕制服av| 日本黄色日本黄色录像| 桃花免费在线播放| 欧美激情久久久久久爽电影 | 精品少妇一区二区三区视频日本电影| 亚洲第一青青草原| 亚洲美女黄色视频免费看| 色视频在线一区二区三区| 韩国高清视频一区二区三区| 女人久久www免费人成看片| 久久精品久久久久久噜噜老黄| 久久99热这里只频精品6学生| 成人av一区二区三区在线看 | 天天躁日日躁夜夜躁夜夜| 日本五十路高清| 久久国产精品大桥未久av| 精品国产一区二区久久| 无限看片的www在线观看| 啦啦啦免费观看视频1| 国产免费福利视频在线观看| 久久久精品94久久精品| 国产欧美日韩一区二区三区在线| 黑丝袜美女国产一区| av片东京热男人的天堂| 不卡av一区二区三区| tube8黄色片| 日本a在线网址| 亚洲欧洲日产国产| 亚洲精品国产区一区二| 久久人妻福利社区极品人妻图片| 午夜成年电影在线免费观看| 欧美日韩黄片免| 欧美日本中文国产一区发布| 少妇的丰满在线观看| 在线 av 中文字幕| 欧美老熟妇乱子伦牲交| 日韩制服骚丝袜av| 国产成+人综合+亚洲专区| 天天影视国产精品| 桃红色精品国产亚洲av| 色视频在线一区二区三区| 亚洲专区中文字幕在线| 51午夜福利影视在线观看| 精品久久蜜臀av无| 亚洲中文字幕日韩| av有码第一页| 日本猛色少妇xxxxx猛交久久| 国产精品免费视频内射| 岛国毛片在线播放| 日韩中文字幕视频在线看片| 国产无遮挡羞羞视频在线观看| 岛国毛片在线播放| 99久久精品国产亚洲精品| 日韩免费高清中文字幕av| 在线观看一区二区三区激情| 一区二区三区精品91| 国产成人av激情在线播放| 国产欧美日韩一区二区三区在线| 99国产精品一区二区蜜桃av | 十八禁人妻一区二区| 黄色a级毛片大全视频| 狠狠精品人妻久久久久久综合| 五月开心婷婷网| 性少妇av在线| 日韩 亚洲 欧美在线| 久久天堂一区二区三区四区| 色婷婷av一区二区三区视频| 水蜜桃什么品种好| 国产在线免费精品| 多毛熟女@视频| 在线观看www视频免费| 国产精品九九99| 国产精品亚洲av一区麻豆| 黄色视频在线播放观看不卡| www.999成人在线观看| videos熟女内射| 国产一区二区三区综合在线观看| 亚洲精品乱久久久久久| 欧美成人午夜精品| 精品久久久久久久毛片微露脸 | 丝袜喷水一区| 男人操女人黄网站| 老熟妇乱子伦视频在线观看 | 男女国产视频网站| 丝袜喷水一区| 免费人妻精品一区二区三区视频| 我要看黄色一级片免费的| 欧美激情 高清一区二区三区| 99久久国产精品久久久| 欧美日韩精品网址| 亚洲中文av在线| 大片电影免费在线观看免费| 国产深夜福利视频在线观看| 亚洲av欧美aⅴ国产| 窝窝影院91人妻| 午夜影院在线不卡| 欧美成狂野欧美在线观看| 精品久久久精品久久久| 国产欧美日韩一区二区三 | 婷婷色av中文字幕| 最近最新中文字幕大全免费视频| 真人做人爱边吃奶动态| 精品视频人人做人人爽| 亚洲国产av影院在线观看| 日韩 亚洲 欧美在线| 成人亚洲精品一区在线观看| 色视频在线一区二区三区| 国产精品麻豆人妻色哟哟久久| 日韩大码丰满熟妇| 一边摸一边抽搐一进一出视频| 欧美日本中文国产一区发布| 国产亚洲午夜精品一区二区久久| 精品亚洲成国产av| bbb黄色大片| 国产精品久久久久成人av| 丝袜喷水一区| 黄色视频,在线免费观看| cao死你这个sao货| 性少妇av在线| 激情视频va一区二区三区| 国产精品自产拍在线观看55亚洲 | 亚洲精品乱久久久久久| 国产无遮挡羞羞视频在线观看| 成年美女黄网站色视频大全免费| 中文字幕av电影在线播放| 久热这里只有精品99| 精品国产一区二区久久| 久久久久久免费高清国产稀缺| 狠狠婷婷综合久久久久久88av| 天堂8中文在线网| 天天影视国产精品| 操美女的视频在线观看| 久久久久久人人人人人| 99国产精品一区二区蜜桃av | 777米奇影视久久| 少妇粗大呻吟视频| 超色免费av| 日韩免费高清中文字幕av| 中文字幕人妻丝袜一区二区| 女人精品久久久久毛片| 999久久久精品免费观看国产| 狠狠狠狠99中文字幕| 日日爽夜夜爽网站| 亚洲精品久久成人aⅴ小说| 99久久国产精品久久久| 少妇 在线观看| 免费高清在线观看日韩| 亚洲精品国产色婷婷电影| 91成人精品电影| 国产黄频视频在线观看| 免费日韩欧美在线观看| 久久久久精品人妻al黑| 国产成人啪精品午夜网站| 黑人巨大精品欧美一区二区蜜桃| 久久久久国内视频| 丰满少妇做爰视频| 99精品欧美一区二区三区四区| 少妇的丰满在线观看| 中文字幕制服av| 国产黄色免费在线视频| 老司机深夜福利视频在线观看 | 一本久久精品| 欧美黑人欧美精品刺激| 亚洲成国产人片在线观看| 12—13女人毛片做爰片一| 纵有疾风起免费观看全集完整版| 一级毛片精品| 91麻豆精品激情在线观看国产 | 欧美激情 高清一区二区三区| 可以免费在线观看a视频的电影网站| 精品国产超薄肉色丝袜足j| 久久99一区二区三区| 国产99久久九九免费精品| 国产免费av片在线观看野外av| 久久这里只有精品19| 午夜两性在线视频| 大型av网站在线播放| 9191精品国产免费久久| 久久av网站| 国产97色在线日韩免费| 国产无遮挡羞羞视频在线观看| 啪啪无遮挡十八禁网站| 女人爽到高潮嗷嗷叫在线视频| 久久久水蜜桃国产精品网| 法律面前人人平等表现在哪些方面 | 久久人妻熟女aⅴ| 啦啦啦视频在线资源免费观看| 久久精品熟女亚洲av麻豆精品| 日本一区二区免费在线视频| 秋霞在线观看毛片| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品av麻豆狂野| 热99久久久久精品小说推荐| 777米奇影视久久| 国产精品影院久久| 一区二区三区精品91| a 毛片基地| 国产精品av久久久久免费| 久久久久久久久免费视频了| 久久久精品免费免费高清| 国产97色在线日韩免费| 老熟妇仑乱视频hdxx| 91国产中文字幕| 国产免费一区二区三区四区乱码| 亚洲伊人色综图| 亚洲一区中文字幕在线| 亚洲精品一区蜜桃| 男女下面插进去视频免费观看| 搡老乐熟女国产| 悠悠久久av| 久久精品国产亚洲av高清一级| 国产99久久九九免费精品| 涩涩av久久男人的天堂| 丝袜人妻中文字幕| 黑人操中国人逼视频| 久久精品aⅴ一区二区三区四区| 国产亚洲精品久久久久5区| 精品一区在线观看国产| 欧美 日韩 精品 国产| 永久免费av网站大全| 91麻豆av在线| 国产一级毛片在线| 亚洲精品中文字幕一二三四区 | 69av精品久久久久久 | 午夜激情av网站| 亚洲精品美女久久av网站| 老司机午夜十八禁免费视频| 久久久久国产精品人妻一区二区| 老司机影院毛片| 国产高清国产精品国产三级| 午夜激情av网站| 19禁男女啪啪无遮挡网站| 18禁黄网站禁片午夜丰满| 亚洲情色 制服丝袜| 老熟妇乱子伦视频在线观看 | tocl精华| 高清欧美精品videossex| 久久精品aⅴ一区二区三区四区| 啪啪无遮挡十八禁网站| 正在播放国产对白刺激| 欧美日韩黄片免| 人人妻人人爽人人添夜夜欢视频| 亚洲色图 男人天堂 中文字幕| 日日爽夜夜爽网站| 女人高潮潮喷娇喘18禁视频| 69av精品久久久久久 | 91老司机精品| av线在线观看网站| 老鸭窝网址在线观看| 日日夜夜操网爽| 少妇精品久久久久久久| 电影成人av| 黄色毛片三级朝国网站| 欧美精品高潮呻吟av久久| 久久九九热精品免费| 亚洲av片天天在线观看| 欧美av亚洲av综合av国产av| 老司机影院毛片| 青春草视频在线免费观看| 国产成人精品久久二区二区91| 丰满人妻熟妇乱又伦精品不卡| 婷婷色av中文字幕| 亚洲九九香蕉| 午夜久久久在线观看| 动漫黄色视频在线观看| 黄色视频,在线免费观看| 大片电影免费在线观看免费| 悠悠久久av| 精品国产一区二区三区久久久樱花| 中文字幕色久视频| 建设人人有责人人尽责人人享有的| 色94色欧美一区二区| 黑人猛操日本美女一级片| 国产欧美日韩一区二区三区在线| 99国产精品一区二区蜜桃av | 在线天堂中文资源库| 十分钟在线观看高清视频www| 又紧又爽又黄一区二区| 啪啪无遮挡十八禁网站| 一边摸一边抽搐一进一出视频| 男女下面插进去视频免费观看| 欧美国产精品va在线观看不卡| 亚洲国产精品一区二区三区在线| 人妻人人澡人人爽人人| 国产伦人伦偷精品视频| 亚洲va日本ⅴa欧美va伊人久久 | 黑人欧美特级aaaaaa片| 久久精品aⅴ一区二区三区四区| 高清av免费在线| 日韩 欧美 亚洲 中文字幕| 国产男人的电影天堂91| 一区福利在线观看| 午夜久久久在线观看| 80岁老熟妇乱子伦牲交| 黑人猛操日本美女一级片| 国产欧美日韩一区二区三 | 精品一品国产午夜福利视频| av又黄又爽大尺度在线免费看| 亚洲成人免费av在线播放| 桃花免费在线播放| 国产亚洲精品第一综合不卡| 狂野欧美激情性xxxx| 国产成人系列免费观看| 首页视频小说图片口味搜索| 亚洲国产毛片av蜜桃av| av又黄又爽大尺度在线免费看| 日韩电影二区| 午夜成年电影在线免费观看| 蜜桃在线观看..| 国产一区二区 视频在线| 亚洲精品一二三| 人人澡人人妻人| 亚洲男人天堂网一区| 国产成人av激情在线播放| 十八禁网站免费在线| 看免费av毛片| av有码第一页| 搡老熟女国产l中国老女人| 91av网站免费观看| 欧美精品啪啪一区二区三区 | 久久精品亚洲熟妇少妇任你| 男女下面插进去视频免费观看| 国产精品一区二区精品视频观看| 午夜福利在线观看吧| 在线观看www视频免费| 一区在线观看完整版| www.自偷自拍.com| 日本vs欧美在线观看视频| 欧美av亚洲av综合av国产av| 欧美一级毛片孕妇| 少妇被粗大的猛进出69影院| 亚洲国产精品一区三区| 免费在线观看完整版高清| 久久青草综合色| 国产欧美日韩精品亚洲av| 啦啦啦在线免费观看视频4| 亚洲精品一二三| av国产精品久久久久影院| 久久久水蜜桃国产精品网| 91精品伊人久久大香线蕉| 黄色视频,在线免费观看| 国产成人精品久久二区二区免费| xxxhd国产人妻xxx| 中文字幕人妻丝袜一区二区| 久久国产亚洲av麻豆专区| 美女高潮到喷水免费观看| 一本综合久久免费| 国产福利在线免费观看视频| 97精品久久久久久久久久精品| 欧美成狂野欧美在线观看| 老汉色∧v一级毛片| 亚洲国产av影院在线观看| 免费少妇av软件| 亚洲国产中文字幕在线视频| 亚洲成av片中文字幕在线观看| 精品少妇久久久久久888优播| 91成年电影在线观看| 欧美国产精品va在线观看不卡| 一级片免费观看大全| 亚洲专区国产一区二区| 久久九九热精品免费| 日韩欧美国产一区二区入口| 亚洲精品乱久久久久久| 国产一卡二卡三卡精品| 91精品国产国语对白视频| 亚洲人成77777在线视频| 国产精品熟女久久久久浪| 日韩大码丰满熟妇| 狠狠精品人妻久久久久久综合| 亚洲精品乱久久久久久| 国产精品国产av在线观看| 久久久久久久精品精品| 91成年电影在线观看| 老司机深夜福利视频在线观看 | 不卡一级毛片| 免费不卡黄色视频| 久久亚洲精品不卡| 精品少妇久久久久久888优播| 国产av精品麻豆| 手机成人av网站| 亚洲美女黄色视频免费看| 国产成人欧美| 日韩大片免费观看网站| a级片在线免费高清观看视频| 国产欧美日韩综合在线一区二区| 女人高潮潮喷娇喘18禁视频| 高清在线国产一区| 欧美激情极品国产一区二区三区| 热99国产精品久久久久久7| 超色免费av| 亚洲全国av大片| 久久久久视频综合| 宅男免费午夜| netflix在线观看网站| av网站免费在线观看视频| 999久久久精品免费观看国产| 少妇人妻久久综合中文| 国产主播在线观看一区二区| 淫妇啪啪啪对白视频 | 国产免费视频播放在线视频| 精品免费久久久久久久清纯 | 9热在线视频观看99| 国产精品成人在线| 久久久国产欧美日韩av| 欧美另类一区| 五月开心婷婷网| 亚洲va日本ⅴa欧美va伊人久久 | 久久久久国产精品人妻一区二区| 成人黄色视频免费在线看| 男女午夜视频在线观看| 国产亚洲午夜精品一区二区久久| 美女高潮到喷水免费观看| 不卡av一区二区三区| 国产成人影院久久av| 久久人人爽人人片av| 91国产中文字幕| 妹子高潮喷水视频| 制服诱惑二区| 在线亚洲精品国产二区图片欧美| bbb黄色大片| 中文字幕制服av| 亚洲欧美成人综合另类久久久| 成人av一区二区三区在线看 | 妹子高潮喷水视频| 桃红色精品国产亚洲av| 亚洲精品在线美女| 国产黄频视频在线观看| 婷婷色av中文字幕| 亚洲天堂av无毛| 免费在线观看影片大全网站| 欧美日韩视频精品一区| 国产一级毛片在线| 久久香蕉激情| 丝瓜视频免费看黄片| 三级毛片av免费| videos熟女内射| 可以免费在线观看a视频的电影网站| 久久国产精品男人的天堂亚洲| av片东京热男人的天堂| 欧美av亚洲av综合av国产av| 亚洲成国产人片在线观看| 国产成人系列免费观看| 亚洲国产欧美在线一区| 在线观看免费视频网站a站| 久久午夜综合久久蜜桃| 在线观看免费视频网站a站| 欧美国产精品一级二级三级| 巨乳人妻的诱惑在线观看| 成人三级做爰电影| 亚洲少妇的诱惑av| 亚洲国产av新网站| 成年动漫av网址| 久久99热这里只频精品6学生| 满18在线观看网站| 成人av一区二区三区在线看 | a级片在线免费高清观看视频| 黑人操中国人逼视频| 亚洲人成电影观看| 日本91视频免费播放| 正在播放国产对白刺激| 欧美激情高清一区二区三区| 韩国高清视频一区二区三区| av又黄又爽大尺度在线免费看| 亚洲精品美女久久av网站| 欧美精品啪啪一区二区三区 | 一级黄色大片毛片| 老司机午夜十八禁免费视频| 国产精品国产av在线观看| 亚洲国产欧美网| 亚洲精品日韩在线中文字幕| 十八禁高潮呻吟视频| 欧美大码av| 香蕉丝袜av| 欧美成人午夜精品| 精品久久蜜臀av无| 日韩三级视频一区二区三区| 欧美日韩av久久| 少妇人妻久久综合中文| 国产欧美日韩一区二区三 | av欧美777| 亚洲 欧美一区二区三区| 日韩人妻精品一区2区三区| 丝袜在线中文字幕| 亚洲精品国产一区二区精华液| 90打野战视频偷拍视频| 一边摸一边做爽爽视频免费| 日韩中文字幕视频在线看片| 国产一区二区三区综合在线观看| 国产精品熟女久久久久浪| 国产国语露脸激情在线看| 久久精品aⅴ一区二区三区四区| 精品高清国产在线一区| 在线观看免费日韩欧美大片| 美女午夜性视频免费| av天堂久久9| 肉色欧美久久久久久久蜜桃| 纯流量卡能插随身wifi吗| svipshipincom国产片| 久久国产精品大桥未久av| 欧美日韩黄片免| 大片电影免费在线观看免费| 色94色欧美一区二区| 性色av一级| 成人三级做爰电影| 日韩熟女老妇一区二区性免费视频| 麻豆乱淫一区二区| 中文字幕色久视频| 久久久国产精品麻豆| 久久久精品区二区三区| 可以免费在线观看a视频的电影网站| 五月开心婷婷网| 天堂俺去俺来也www色官网| 中文字幕色久视频| 国产精品九九99| 少妇 在线观看| 国产熟女午夜一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 精品一区二区三区av网在线观看 |