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

    日地L1點(diǎn)探測任務(wù)設(shè)計與結(jié)果分析

    2021-12-21 09:19:38鄒樂洋高珊趙晨喬德治李曉光孟占峰
    中國空間科學(xué)技術(shù) 2021年6期
    關(guān)鍵詞:測控姿態(tài)坐標(biāo)系

    鄒樂洋,高珊,趙晨,喬德治,李曉光,孟占峰

    1. 北京空間飛行器總體設(shè)計部,北京 100094 2. 上海宇航系統(tǒng)工程研究所,上海 201109 3. 北京控制工程研究所,北京 100080

    1 引言

    日地L1點(diǎn)具有器-地-日相對位置關(guān)系較為固定、能源條件穩(wěn)定、測控通信連續(xù)性好、軌道維持速度增量需求小等特點(diǎn),是開展太陽觀測任務(wù)的理想觀測點(diǎn)。目前國際上已實(shí)施的日地L1點(diǎn)探測任務(wù)包括歐美合作的ISEE-3任務(wù)和SOHO任務(wù)、美國的Wind任務(wù)和Ace任務(wù)等[1],中國正在論證的“夸父”任務(wù)也選擇日地L1點(diǎn)開展太陽觀測任務(wù)。

    日地L1點(diǎn)對于我國后續(xù)太陽觀測任務(wù)具有重大意義。我國目前尚未實(shí)施過日地L1點(diǎn)探測任務(wù),對日地L1點(diǎn)的軌道設(shè)計、測控鏈路環(huán)境、太陽輻照環(huán)境、三體動力學(xué)環(huán)境、空間輻射環(huán)境等方面的研究缺少在軌實(shí)際飛行數(shù)據(jù)的支撐,因此,開展一次日地L1點(diǎn)飛行驗(yàn)證任務(wù),獲得在軌實(shí)測結(jié)果,對于支撐我國后續(xù)太陽觀測任務(wù)設(shè)計與實(shí)施具有重要意義。

    在嫦娥五號任務(wù)的既定目標(biāo)全部實(shí)現(xiàn)后,利用軌道器能力和剩余推進(jìn)劑,在拓展任務(wù)階段設(shè)計并實(shí)施了我國首次日地L1點(diǎn)探測,對日地L1點(diǎn)李薩茹軌道(Lissajous orbit)特性以及日地L1點(diǎn)的測控鏈路環(huán)境、太陽輻照環(huán)境等飛行環(huán)境進(jìn)行了探測,對太陽翼撓性動力學(xué)參數(shù)設(shè)計、環(huán)境力擾動影響和工業(yè)級元器件在空間輻照環(huán)境下的工作性能進(jìn)行了實(shí)際飛行驗(yàn)證,為我國未來太陽觀測等深空探測任務(wù)的設(shè)計提供在軌實(shí)際飛行數(shù)據(jù)。

    本文對日地L1點(diǎn)探測任務(wù)的設(shè)計與結(jié)果分析情況進(jìn)行介紹,總結(jié)我國首次日地L1點(diǎn)探測任務(wù)的技術(shù)成果,為未來深空探測任務(wù)設(shè)計提供參考。

    2 日地L1點(diǎn)探測實(shí)施過程

    嫦娥五號軌道器在2020年12月17日與返回器分離后,通過規(guī)避機(jī)動和軌道入射機(jī)動進(jìn)入了日地L1點(diǎn)轉(zhuǎn)移軌道,開始實(shí)施拓展任務(wù)。軌道器經(jīng)過約3個月的轉(zhuǎn)移軌道飛行,期間執(zhí)行了2次中途修正,于2021年3月15日到達(dá)日地L1點(diǎn),進(jìn)入環(huán)繞日地L1點(diǎn)運(yùn)行的李薩茹軌道。軌道器在日地L1點(diǎn)環(huán)繞期間完成了日地L1點(diǎn)李薩茹軌道飛行試驗(yàn)、日地L1點(diǎn)飛行環(huán)境試驗(yàn)、太陽翼撓性在軌辨識試驗(yàn)、光學(xué)相機(jī)成像試驗(yàn)等試驗(yàn)項(xiàng)目,試驗(yàn)結(jié)果正常,獲取了預(yù)期的試驗(yàn)數(shù)據(jù)。日地L1點(diǎn)探測實(shí)施過程的關(guān)鍵事件如表1所示。

    表1 日地L1點(diǎn)探測實(shí)施過程關(guān)鍵事件

    3 日地L1點(diǎn)探測設(shè)計

    嫦娥五號拓展任務(wù)的主要目的是獲取在軌飛行試驗(yàn)數(shù)據(jù),為未來深空探測任務(wù)和產(chǎn)品的設(shè)計提供參考。在拓展任務(wù)設(shè)計中,避免與嫦娥二號、嫦娥五號飛行試驗(yàn)器等拓展任務(wù)目標(biāo)發(fā)生重復(fù)[2-5],獲得更多全新的探測結(jié)果。因此,選擇了我國之前從未開展過探測活動的日地L1點(diǎn)作為嫦娥五號拓展任務(wù)的探測目標(biāo)。嫦娥五號拓展任務(wù)的具體內(nèi)容包括:

    1)驗(yàn)證日地L1點(diǎn)(距離地球約1.5×106km)轉(zhuǎn)移軌道的設(shè)計與控制技術(shù);

    2)在日地L1點(diǎn)附近開展環(huán)繞飛行任務(wù),驗(yàn)證L1點(diǎn)環(huán)繞軌道設(shè)計與控制技術(shù);

    3)對日地L1點(diǎn)附近測控鏈路、太陽輻照、三體動力學(xué)、空間輻射等飛行環(huán)境進(jìn)行監(jiān)測,獲取實(shí)際在軌飛行數(shù)據(jù)。

    3.1 日地L1點(diǎn)李薩茹軌道飛行試驗(yàn)

    日地L1點(diǎn)李薩茹軌道飛行試驗(yàn)的目的是通過設(shè)計并實(shí)施日地L1點(diǎn)轉(zhuǎn)移和環(huán)繞飛行,驗(yàn)證軌道設(shè)計和軌道控制策略的正確性、可行性。具體飛行過程可劃分為以下兩個階段:

    1)地球-日地L1點(diǎn)轉(zhuǎn)移段。此階段從軌道器規(guī)避機(jī)動后開始,至達(dá)到日地L1點(diǎn)附近過日地L1點(diǎn)旋轉(zhuǎn)坐標(biāo)系XOZ面為止。規(guī)避機(jī)動后2 h實(shí)施一次轉(zhuǎn)移軌道入射機(jī)動,目的是通過一次變軌使軌道器進(jìn)入飛往日地L1點(diǎn)的轉(zhuǎn)移軌道,進(jìn)入飛往日地L1點(diǎn)的穩(wěn)定流形,軌道器飛行約90 d到達(dá)日地L1點(diǎn)附近(第一次過日地L1點(diǎn)旋轉(zhuǎn)坐標(biāo)系XOZ面),直接進(jìn)入環(huán)繞L1點(diǎn)的李薩茹軌道。

    2)日地L1點(diǎn)環(huán)繞飛行段。軌道器第一次過日地L1點(diǎn)旋轉(zhuǎn)坐標(biāo)系XOZ面后即進(jìn)入環(huán)繞日地L1點(diǎn)的李薩茹軌道運(yùn)行,一圈周期約為180 d,每半圈(約3個月)安排一次軌道維持,軌道器在該軌道上運(yùn)行并開展相關(guān)的驗(yàn)證試驗(yàn)任務(wù)。

    (1)日地L1點(diǎn)轉(zhuǎn)移軌道設(shè)計

    從地球出發(fā)的日地平動點(diǎn)轉(zhuǎn)移軌道設(shè)計主要考慮地球、月球、太陽的引力以及其它大天體的三體引力,此外還應(yīng)考慮太陽光壓的影響。轉(zhuǎn)移軌道的設(shè)計采用數(shù)值方法[3],根據(jù)嫦娥五號軌道器規(guī)避機(jī)動完成后的軌道參數(shù),通過一次變軌控制調(diào)整出發(fā)速度來滿足軌道器到達(dá)平動點(diǎn)時的速度需求。

    定義初始狀態(tài)為軌道器的出發(fā)速度矢量v0:

    p=v0

    (1)

    將終端狀態(tài)矩陣定義為:

    q=v1x

    (2)

    式中:v1x為軌道器第一次到達(dá)日地L1點(diǎn)旋轉(zhuǎn)坐標(biāo)系XOZ面時刻相對旋轉(zhuǎn)坐標(biāo)系X方向的速度大小,目標(biāo)值為0。日地L1旋轉(zhuǎn)坐標(biāo)系的定義為:原點(diǎn)在日地L1點(diǎn),X軸始終沿太陽指向地球方向,Z軸沿地球公轉(zhuǎn)平面法線方向,Y軸和Z軸、X軸構(gòu)成右手螺旋坐標(biāo)系。

    根據(jù)初始狀態(tài)與終端狀態(tài)間的關(guān)系q=q(p),兩者的誤差關(guān)系可線性化表達(dá)為:

    (3)

    利用上述關(guān)系,通過若干次迭代即可獲得滿足第一次過XOZ面時刻速度需求的軌道器初始出發(fā)速度。再以該初始出發(fā)狀態(tài)作為下一步迭代的初值,采用相同方法求解滿足軌道器第二次到達(dá)日地L1點(diǎn)旋轉(zhuǎn)坐標(biāo)系XOZ面時刻,相對旋轉(zhuǎn)坐標(biāo)系X方向的速度大小v1x=0的初始狀態(tài),以此類推,計算至第五次過旋轉(zhuǎn)坐標(biāo)系XOZ面時刻停止[3]。

    利用上述方法計算獲得軌道器入射進(jìn)入轉(zhuǎn)移軌道所需的速度增量。軌道器在完成規(guī)避機(jī)動后2 h,通過補(bǔ)充實(shí)施一次轉(zhuǎn)移軌道入射機(jī)動,即可進(jìn)入飛往日地L1點(diǎn)的穩(wěn)定流形。軌道器飛行約90 d到達(dá)日地L1點(diǎn)附近(第一次過日地L1點(diǎn)旋轉(zhuǎn)坐標(biāo)系XOZ面),進(jìn)入環(huán)繞L1點(diǎn)的李薩茹軌道。在轉(zhuǎn)移軌道飛行過程中安排了三次中途修正,第一次安排在規(guī)避后約10 h,用于修正規(guī)避控制殘差;第二次安排在規(guī)避機(jī)動后約40 d,用于消除前一次修正的殘差;第三次安排在第一次過日地L1點(diǎn)旋轉(zhuǎn)坐標(biāo)系XOZ面時,用于消除前一次修正的殘差,保證進(jìn)入李薩茹軌道后的穩(wěn)定性[6-8]。

    中途修正策略均為通過調(diào)整軌控速度增量的三方向分量瞄準(zhǔn)控后第五次過日地L1點(diǎn)旋轉(zhuǎn)坐標(biāo)系XOZ面時的X向速度為零。

    (2)日地L1點(diǎn)李薩茹軌道設(shè)計

    在圓限制性三體問題下,共線平動點(diǎn)附近的軌道運(yùn)行一般表達(dá)為:

    (4)

    式中:x、y和z分別為探測器在旋轉(zhuǎn)坐標(biāo)系下的三方向位置分量;Ax和Az分別表示平動點(diǎn)軌道在旋轉(zhuǎn)坐標(biāo)系下x向和z向的軌道振幅;λ為x向振動頻率;ν為z向振動頻率;Φ和φ為振動的初始相位。

    考慮到進(jìn)入Halo軌道相比李薩茹軌道需要消耗更多的推進(jìn)劑,同時拓展任務(wù)軌道本身沒有特定的構(gòu)形要求,因此在軌道設(shè)計中選擇了經(jīng)過穩(wěn)定流形到達(dá)后自然形成的李薩茹軌道作為使命軌道。

    軌道器第一次過日地L1點(diǎn)旋轉(zhuǎn)坐標(biāo)系XOZ面后即進(jìn)入環(huán)繞日地L1點(diǎn)的李薩茹軌道運(yùn)行,一圈周期約為180 d,每半圈(約3個月)安排一次軌道維持,維持策略為通過調(diào)整軌控速度增量的三方向分量瞄準(zhǔn)控后第五次過日地L1點(diǎn)旋轉(zhuǎn)坐標(biāo)系XOZ面時的X向速度為零。軌道器可在該軌道上運(yùn)行一圈或數(shù)圈,開展相關(guān)的探測任務(wù)[9]。

    日地L1點(diǎn)探測的飛行軌跡如圖1所示,李薩茹軌道各向振幅如圖2所示。李薩茹軌道X方向振幅約4.5×105km、Y方向振幅約1.4×106km、Z方向振幅約1.5×105km。

    圖1 日地L1點(diǎn)探測飛行軌跡Fig.1 Flight trajectory of Sun-Earth L1 point exploration

    圖2 日地L1點(diǎn)李薩茹軌道各向振幅Fig.2 The orbit amplitude on each direction of the Lissajous orbit around the Sun-Earth L1 point

    3.2 日地L1點(diǎn)飛行環(huán)境探測

    日地L1點(diǎn)飛行環(huán)境探測的主要目的是通過日地L1點(diǎn)環(huán)繞飛行,對日地L1點(diǎn)附近的測控鏈路環(huán)境、太陽輻照環(huán)境等飛行環(huán)境進(jìn)行探測與驗(yàn)證。

    在測控鏈路環(huán)境探測與驗(yàn)證方面,主要通過日地L1點(diǎn)飛行期間的長期測控跟蹤,驗(yàn)證日地L1點(diǎn)飛行環(huán)境條件下的測控鏈路性能,通過低SEP角(即太陽-地球-探測器夾角,一般認(rèn)為SEP角低于5°時會發(fā)生日凌現(xiàn)象,太陽噪聲影響地面接收下行信號)條件下的測控跟蹤以及擇機(jī)開展高碼速率模式下行試驗(yàn),驗(yàn)證測控鏈路對日地L1點(diǎn)飛行環(huán)境的適應(yīng)性。

    在太陽輻照環(huán)境探測與驗(yàn)證方面,主要使用太陽翼輸出電流作為太陽輻照強(qiáng)度的表征,將日地L1點(diǎn)環(huán)繞飛行期間實(shí)際太陽翼輸出電流的變化情況與基于日地距離變化的日地因子模型的預(yù)示結(jié)果進(jìn)行比對,驗(yàn)證日地因子模型的正確性。

    3.3 三體動力學(xué)環(huán)境探測

    日地L1點(diǎn)是太陽-地球的引力平動點(diǎn),繞日地L1點(diǎn)的飛行環(huán)境是典型的三體動力學(xué)環(huán)境。在日地L1點(diǎn)探測期間,設(shè)計了太陽翼撓性在軌辨識試驗(yàn),并通過動量輪角動量累積情況對太陽光壓等環(huán)境力產(chǎn)生的姿態(tài)擾動情況進(jìn)行了探測。

    軌道器以三軸穩(wěn)定控制模式在軌運(yùn)行,巡航期間主要使用動量輪實(shí)現(xiàn)高精度姿態(tài)維持,而周期性的動量輪控制會激發(fā)太陽翼撓性振動,成為影響姿態(tài)維持精度的主要因素之一。由于地面試驗(yàn)時存在重力影響,難以真實(shí)模擬在軌力學(xué)環(huán)境,因此很難通過地面試驗(yàn)獲得真實(shí)的太陽翼撓性振動數(shù)據(jù),一般采用仿真分析的方式獲得相關(guān)參數(shù),用于控制系統(tǒng)參數(shù)設(shè)計[10-11]。利用日地L1點(diǎn)三體動力學(xué)環(huán)境的失重、天體引力場作用穩(wěn)定、引力梯度影響小等特點(diǎn),設(shè)計了太陽翼撓性在軌辨識試驗(yàn),通過修改動量輪控制參數(shù)主動激發(fā)軌道器太陽翼撓性振動,根據(jù)軌道器整器姿態(tài)角速度變化情況反演0°、315°、270°等不同太陽翼轉(zhuǎn)角下的太陽翼撓性振動情況,實(shí)現(xiàn)太陽翼的動力學(xué)阻尼特性在軌辨識。

    另一方面,在三軸穩(wěn)定控制條件下,太陽光壓等環(huán)境力產(chǎn)生的姿態(tài)擾動力矩是通過動量輪加速轉(zhuǎn)動產(chǎn)生的控制力矩進(jìn)行抵消的,因此姿態(tài)擾動力矩的作用可以由動量輪角動量累積情況進(jìn)行表征。通過對動量輪角動量累積情況進(jìn)行分析,可以對環(huán)境力產(chǎn)生的姿態(tài)擾動情況進(jìn)行探測。

    3.4 日地L1點(diǎn)輻射環(huán)境對工業(yè)元器件影響試驗(yàn)

    軌道器配備的2臺光學(xué)相機(jī)(分別為寬視場監(jiān)視攝像機(jī)和監(jiān)視傳感器B)為實(shí)現(xiàn)輕小型化設(shè)計使用了部分工業(yè)級元器件。為驗(yàn)證工業(yè)級元器件在日地L1點(diǎn)輻射環(huán)境下的工作性能[12-13],設(shè)計了日地L1點(diǎn)輻射環(huán)境對工業(yè)元器件影響試驗(yàn),使用2臺光學(xué)相機(jī)實(shí)施地月合影拍攝、太陽拍攝和太陽翼狀態(tài)圖像監(jiān)視[14]。

    4 日地L1點(diǎn)探測結(jié)果分析

    4.1 日地L1點(diǎn)李薩茹軌道飛行試驗(yàn)

    日地L1點(diǎn)李薩茹軌道飛行試驗(yàn)期間共實(shí)施了1次軌道入射機(jī)動、2次轉(zhuǎn)移軌道中途修正(設(shè)計中的第3次中途修正因速度增量極小而取消)和1次李薩茹軌道維持,各次軌控的控制參數(shù)如表2所示。

    表2 日地L1點(diǎn)李薩茹軌道飛行試驗(yàn)軌道控制參數(shù)

    根據(jù)實(shí)際飛行獲得的定軌數(shù)據(jù),對日地L1點(diǎn)的飛行情況進(jìn)行了外推分析,李薩茹軌道的各向位置振幅偏差結(jié)果見圖3~圖5,圖中的歷元時刻為軌道器到達(dá)日地L1點(diǎn)時刻。

    圖3 旋轉(zhuǎn)坐標(biāo)系X向位置偏差Fig.3 Position deviation in the X direction of rotating coordinate system

    由圖3~圖5可以看出,在日地L1點(diǎn)環(huán)繞飛行期間,與李薩茹軌道X方向設(shè)計振幅約4.5×105km,Y方向設(shè)計振幅約1.4×106km,Z方向設(shè)計振幅約1.5×105km相比,實(shí)際飛行的李薩茹軌道在日地旋轉(zhuǎn)坐標(biāo)系X向和Y向的位置偏差均不超過2 000 km,Z向偏差不超過200 km,各向位置的相對偏差不到0.5%,與設(shè)計的一致性良好,天體引力、太陽光壓等環(huán)境力的作用符合預(yù)期。

    圖4 旋轉(zhuǎn)坐標(biāo)系Y向位置偏差Fig.4 Position deviation in the Y direction of rotating coordinate system

    圖5 旋轉(zhuǎn)坐標(biāo)系Z向位置偏差Fig.5 Position deviation in the Z direction of rotating coordinate system

    4.2 日地L1點(diǎn)飛行環(huán)境探測

    (1)測控鏈路環(huán)境探測與驗(yàn)證

    嫦娥五號軌道器上測控系統(tǒng)總接收功率由式(5)給出:

    EIRP|g+LTP+LSP|u+La+LP+

    (5)

    式(5)中,除EIRP|g和LSP|u兩項(xiàng)外,其它項(xiàng)均不受軌道器運(yùn)行位置的影響。軌道器運(yùn)行于月球軌道和日地L1點(diǎn)兩處的接收功率不同,是由于EIRP|g和LSP|u兩項(xiàng)造成的。其中:

    LSP|u=20lgR+20lgfup+32.44

    (6)

    式中:R為軌道器和地面天線之間的距離(km);fup為上行載波頻率(MHz)[15]。經(jīng)計算,器地距離為4.0×105km和1.5×106km時,LSP|u分別為221.7 dB和233.2 dB,即軌道器運(yùn)行至日地L1點(diǎn)時,空間損失較月球軌道增加了11.5 dB。作為補(bǔ)償,地面站提高了上行功率,日地L1點(diǎn)探測期間的有效發(fā)射功率EIRP|g比環(huán)月期間增加了18 dB。綜合考慮EIRP|g和LSP|u兩項(xiàng),由上述公式計算得到的軌道器運(yùn)行于日地L1點(diǎn)處的接收功率應(yīng)比月球軌道高6.5 dB。實(shí)際飛行數(shù)據(jù)分析結(jié)果表明,軌道器環(huán)月飛行階段的應(yīng)答機(jī)接收通道信號強(qiáng)度約為-104~-108 dBm,日地L1點(diǎn)探測期間應(yīng)答機(jī)接收通道信號強(qiáng)度約為-100 dBm,比環(huán)月飛行階段高出4~8 dB,實(shí)際飛行結(jié)果與理論計算結(jié)果基本相符,存在誤差是由于器上接收天線增益隨方位角不同有一定波動。

    日地L1點(diǎn)飛行期間SEP角變化情況如圖6所示,在一圈環(huán)繞軌道飛行過程中,SEP角變化范圍為0°~50°。日地L1點(diǎn)飛行期間,地面站對軌道器進(jìn)行了測控跟蹤。軌道器接收地面站上行信號正常,應(yīng)答機(jī)接收通道信號強(qiáng)度高于-108 dBm,高于解調(diào)門限-128 dBm,地面站對軌道器上行跟蹤鎖定正常。地面站接收軌道器下行信號正常,在低SEP角(低于5°)條件下下行信號接收正常。日地L1點(diǎn)飛行期間開展了高碼速率模式下行試驗(yàn),地面站能夠正常接收軌道器下行數(shù)據(jù),表明測控鏈路性能能夠支持日地L1點(diǎn)飛行環(huán)境下的測控通信,為我國未來日地L1點(diǎn)探測任務(wù)測控鏈路設(shè)計提供了在軌實(shí)際飛行數(shù)據(jù)。

    圖6 日地L1點(diǎn)飛行期間SEP角變化情況Fig.6 Variation of SEP angle during the flight around the Sun-Earth L1 point

    (2)太陽輻照環(huán)境探測與驗(yàn)證

    嫦娥五號軌道器太陽翼輸出電流采用式(7)計算:

    I=I0×Xs×Ns×[βp×(T-25℃)+1]×S

    (7)

    其中,I0為標(biāo)準(zhǔn)條件下測試電流值,Xs為太陽季節(jié)性變化因子,Ns為電池片并聯(lián)數(shù),βp為溫度系數(shù),T為溫度,S為單片電池片有效面積。式(7)中,除Xs外,其余參數(shù)為能源系統(tǒng)自身的設(shè)計值,I0為按照地球附近平均太陽常數(shù)對地面測試結(jié)果進(jìn)行修正后獲得的基準(zhǔn)電流值。在地月系統(tǒng)實(shí)施探測任務(wù)的探測器的能源系統(tǒng)設(shè)計中,太陽常數(shù)一般取平均值,即1 353 W/m2,由于地月距離相比日地距離為小量,因此一般不考慮Xs受探測器與太陽之間的距離變化的影響。對于日地L1點(diǎn)探測任務(wù),器地距離約為1.5×106km,距離變化對Xs造成的影響不可忽略,同時,由于日地距離的季節(jié)性變化,Xs也相應(yīng)季節(jié)性變化,在能源系統(tǒng)設(shè)計中需要加以考慮[16-17]。

    根據(jù)Xs與日地距離之間的反比關(guān)系建立了日地因子模型,如式(8)所示,其中XsL1為日地L1點(diǎn)處的太陽季節(jié)性變化因子,L0為太陽與地球的平均距離,LL1為日地L1點(diǎn)探測期間太陽與探測器之間的距離。根據(jù)模型修正后的太陽季節(jié)性變化因子對太陽翼輸出電流進(jìn)行了理論計算,并將太陽翼輸出電流的預(yù)示值與實(shí)際在軌飛行的實(shí)測值進(jìn)行了比較,如圖7所示。模型預(yù)示值與實(shí)測值之間的最大偏差約為2.5%,表明日地因子模型對于日地L1點(diǎn)的太陽輻照環(huán)境具有良好的預(yù)示作用,能夠應(yīng)用于未來日地L1點(diǎn)探測任務(wù)的能源系統(tǒng)設(shè)計。

    XsL1=Xs×L0/LL1

    圖7 太陽翼輸出電流預(yù)示值與實(shí)測值比較Fig.7 Comparison between the predictive value and the measured value of solar array output current

    4.3 三體動力學(xué)環(huán)境探測

    (1)太陽翼撓性在軌辨識試驗(yàn)

    在太陽翼撓性在軌辨識試驗(yàn)中,注入高帶寬動量輪控制參數(shù)后,軌道器偏航方向姿態(tài)角速度震蕩幅值從小于0.000 3(°)/s增加至最大0.012(°)/s,說明太陽翼內(nèi)彎模態(tài)被動量輪控制激發(fā)。觀察20 min后將太陽翼轉(zhuǎn)至315°轉(zhuǎn)角位置,觀察到滾動和偏航姿態(tài)角速度均出現(xiàn)小于0.01(°)/s幅值的震蕩,說明太陽翼內(nèi)彎模態(tài)同時對滾動方向和偏航方向的姿態(tài)造成耦合影響。再觀察20 min后,將太陽翼轉(zhuǎn)至270°轉(zhuǎn)角位置,觀察到滾動姿態(tài)角速度出現(xiàn)震蕩,說明太陽翼內(nèi)彎模態(tài)對滾動方向姿態(tài)造成影響。再過20 min后,將太陽翼轉(zhuǎn)回0°轉(zhuǎn)角位置,并恢復(fù)低帶寬動量輪控制參數(shù),觀察到偏航姿態(tài)角速度震蕩幅值逐漸減小趨于收斂,整器姿態(tài)角速度恢復(fù)至小于0.000 3(°)/s,說明太陽翼內(nèi)彎模態(tài)振動趨于收斂。試驗(yàn)過程中的本體角速度震蕩情況如圖8所示。

    圖8 試驗(yàn)過程中軌道器三軸姿態(tài)角速度震蕩情況Fig.8 The vibration of angular velocity in three axis on the orbiter during the test

    通過對注入低帶寬參數(shù)后的系統(tǒng)響應(yīng)進(jìn)行分析,對太陽翼振動的真實(shí)阻尼比進(jìn)行了評估。結(jié)果表明,太陽翼阻尼比介于0.002~0.005之間,且與太陽翼振幅正相關(guān)。

    (2)環(huán)境力姿態(tài)擾動情況探測

    嫦娥五號軌道器在日地L1點(diǎn)飛行期間,受到的環(huán)境力主要包括大天體引力和太陽光壓,環(huán)境力作用對軌道器產(chǎn)生姿態(tài)擾動力矩,由動量輪加速轉(zhuǎn)動產(chǎn)生的控制力矩抵消。日地L1點(diǎn)是太陽和地球之間的引力平動點(diǎn),引力場穩(wěn)定,引力梯度造成的姿態(tài)擾動較小,可以忽略,主要的擾動力是太陽光壓。

    太陽光壓主要作用在太陽翼上,壓心可近似為太陽翼的幾何中心。由于軌道器自身存在質(zhì)心偏移,因此會產(chǎn)生姿態(tài)干擾力矩,可以通過軌道器動量輪角動量累積情況進(jìn)行表征。軌道器配備了4臺動量輪,其中包括3臺與軌道器坐標(biāo)系方向相同的正交動量輪和1臺斜裝動量輪。如圖9所示,實(shí)際在軌飛行數(shù)據(jù)表明,角動量累積主要發(fā)生在軌道器-Y方向,角動量累計周期約為2周,與預(yù)期結(jié)果一致。

    圖9 動量輪角動量累積情況Fig.9 The angular momentum accumulation of the momentum wheels

    4.4 日地L1點(diǎn)輻射環(huán)境對工業(yè)元器件影響試驗(yàn)

    2021年4月15日,使用寬視場監(jiān)視攝像機(jī)和監(jiān)視傳感器B完成了地月合影拍攝、太陽拍攝和太陽翼撓性在軌辨識試驗(yàn)過程中的太陽翼狀態(tài)監(jiān)視。

    試驗(yàn)獲得的圖像數(shù)據(jù)質(zhì)量良好,圖像顯示清晰,表明使用了工業(yè)級元器件的光學(xué)相機(jī)在日地L1點(diǎn)輻射環(huán)境下工作狀態(tài)正常。試驗(yàn)過程圖像如圖10~圖12所示。

    圖10 寬視場監(jiān)視攝像機(jī)拍攝的地月合影圖像Fig.10 Group photo of the Earth and the Moon taken by wide-field camera

    圖11 寬視場監(jiān)視攝像機(jī)拍攝的太陽圖像Fig.11 Photo of the Sun taken by wide-field camera

    圖12 監(jiān)視傳感器B拍攝的太陽翼圖像Fig.12 Photo of the solar array taken by monitoring camera B

    5 結(jié)論

    充分利用嫦娥五號軌道器平臺能力和剩余推進(jìn)劑,設(shè)計并實(shí)施了中國首次日地L1點(diǎn)探測任務(wù),完成了日地L1點(diǎn)李薩茹軌道飛行試驗(yàn)、日地L1點(diǎn)飛行環(huán)境探測、三體動力學(xué)環(huán)境探測以及日地L1點(diǎn)輻射環(huán)境對工業(yè)元器件影響試驗(yàn)等項(xiàng)目,獲取了大量技術(shù)成果,進(jìn)一步豐富了嫦娥五號任務(wù)的成果產(chǎn)出,對于中國未來深空探測任務(wù)的設(shè)計具有重要參考價值。通過日地L1點(diǎn)探測的設(shè)計與實(shí)施,主要獲得了以下技術(shù)成果:

    1)實(shí)現(xiàn)了中國首次日地L1點(diǎn)環(huán)繞探測,驗(yàn)證了日地L1點(diǎn)轉(zhuǎn)移軌道、環(huán)繞軌道設(shè)計與控制技術(shù),獲得了日地L1點(diǎn)軌道演化特性數(shù)據(jù),環(huán)繞軌道實(shí)際飛行振幅與設(shè)計值相比偏差不超過0.5%,與設(shè)計的一致性良好,天體引力、太陽光壓等環(huán)境力的作用符合預(yù)期,可以作為未來日地L1點(diǎn)探測任務(wù)的參考。

    2)實(shí)現(xiàn)了對日地L1點(diǎn)測控鏈路環(huán)境、太陽輻照環(huán)境等飛行環(huán)境的就位探測,獲取了在軌實(shí)測數(shù)據(jù),測控系統(tǒng)實(shí)際接收增益與理論計算值基本一致,太陽翼輸出電流實(shí)測值與模型預(yù)測值之間偏差不超過2.5%,為中國未來日地L1點(diǎn)探測任務(wù)的測控系統(tǒng)、能源系統(tǒng)設(shè)計提供了重要參考。

    3)實(shí)現(xiàn)了三體動力學(xué)環(huán)境探測,完成了太陽翼撓性動力學(xué)參數(shù)在軌辨識和環(huán)境力姿態(tài)擾動情況探測,試驗(yàn)結(jié)果與仿真分析結(jié)果基本匹配一致,在未來深空探測器動力學(xué)分析與建模、控制系統(tǒng)參數(shù)設(shè)計中具有較高的參考價值。

    4)實(shí)現(xiàn)了工業(yè)級元器件在日地L1點(diǎn)輻射環(huán)境下的成功應(yīng)用,驗(yàn)證了工業(yè)級圖像傳感器等元器件在日地L1點(diǎn)空間輻射環(huán)境下的工作性能,對未來深空探測器的設(shè)備輕小型、長壽命設(shè)計以及低成本空間產(chǎn)品開發(fā)有重要借鑒意義。

    猜你喜歡
    測控姿態(tài)坐標(biāo)系
    攀爬的姿態(tài)
    《測控電路》實(shí)踐教學(xué)改革探討
    電子測試(2018年22期)2018-12-19 05:12:58
    全新一代宋的新姿態(tài)
    汽車觀察(2018年9期)2018-10-23 05:46:40
    跑與走的姿態(tài)
    中國自行車(2018年8期)2018-09-26 06:53:44
    解密坐標(biāo)系中的平移變換
    坐標(biāo)系背后的故事
    基于重心坐標(biāo)系的平面幾何證明的探討
    基于現(xiàn)代測控技術(shù)及其應(yīng)用分析
    向著新航程進(jìn)發(fā)——遠(yuǎn)望7號測控船首航記錄
    太空探索(2016年12期)2016-07-18 11:13:43
    基于USB2.0協(xié)議的通用測控通信接口設(shè)計
    av女优亚洲男人天堂| 在线观看免费视频网站a站| 我要看黄色一级片免费的| 91精品国产国语对白视频| 中国国产av一级| 国产精品人妻久久久影院| 美女国产视频在线观看| 久久这里只有精品19| 日韩中字成人| 久久久国产精品麻豆| 热99久久久久精品小说推荐| 国产成人午夜福利电影在线观看| 成年女人在线观看亚洲视频| 欧美日韩av久久| 秋霞在线观看毛片| 国产欧美日韩一区二区三区在线| 在线观看www视频免费| 建设人人有责人人尽责人人享有的| 91在线精品国自产拍蜜月| 午夜福利视频精品| 高清欧美精品videossex| 国产免费一区二区三区四区乱码| 亚洲国产最新在线播放| 色播在线永久视频| 精品卡一卡二卡四卡免费| 色视频在线一区二区三区| 日日摸夜夜添夜夜爱| 人妻人人澡人人爽人人| 亚洲精品日本国产第一区| 九色亚洲精品在线播放| 免费在线观看视频国产中文字幕亚洲 | 国产麻豆69| 少妇人妻 视频| 波多野结衣一区麻豆| 国产成人免费观看mmmm| 性少妇av在线| 亚洲欧洲国产日韩| 亚洲一区二区三区欧美精品| 午夜日本视频在线| 亚洲色图 男人天堂 中文字幕| 久久99蜜桃精品久久| 欧美bdsm另类| av免费在线看不卡| 亚洲综合色网址| 精品第一国产精品| 亚洲熟女精品中文字幕| 久久精品国产亚洲av高清一级| 久久久亚洲精品成人影院| 日韩av不卡免费在线播放| 亚洲精品久久久久久婷婷小说| 视频区图区小说| 亚洲成国产人片在线观看| 国产97色在线日韩免费| 一边亲一边摸免费视频| 国产亚洲最大av| 午夜免费鲁丝| videossex国产| 亚洲成人手机| 日本av手机在线免费观看| 99九九在线精品视频| 国产在视频线精品| 国产片特级美女逼逼视频| 国产精品免费视频内射| 超碰97精品在线观看| 三上悠亚av全集在线观看| 亚洲av在线观看美女高潮| 久久精品国产鲁丝片午夜精品| 欧美精品高潮呻吟av久久| 91成人精品电影| 成年动漫av网址| 亚洲,欧美,日韩| 亚洲经典国产精华液单| 免费高清在线观看日韩| 中文字幕另类日韩欧美亚洲嫩草| 亚洲av欧美aⅴ国产| 只有这里有精品99| 亚洲美女视频黄频| 一区在线观看完整版| 国产一区亚洲一区在线观看| 免费av中文字幕在线| 国产精品一区二区在线不卡| 丰满迷人的少妇在线观看| 美女福利国产在线| 精品少妇内射三级| 国产精品久久久久成人av| 欧美国产精品一级二级三级| 久久这里只有精品19| 啦啦啦在线观看免费高清www| 久久人人爽av亚洲精品天堂| 国产老妇伦熟女老妇高清| 亚洲精品第二区| 久久精品国产亚洲av天美| 一区二区日韩欧美中文字幕| 不卡视频在线观看欧美| 亚洲人成77777在线视频| 高清黄色对白视频在线免费看| 日本av免费视频播放| av福利片在线| 大片电影免费在线观看免费| 人人妻人人添人人爽欧美一区卜| 超碰成人久久| 久久久久久免费高清国产稀缺| 亚洲精品第二区| 高清在线视频一区二区三区| 精品一品国产午夜福利视频| 午夜福利乱码中文字幕| 精品酒店卫生间| 90打野战视频偷拍视频| 中文天堂在线官网| 亚洲一级一片aⅴ在线观看| 亚洲伊人久久精品综合| 精品午夜福利在线看| 9色porny在线观看| 亚洲色图综合在线观看| 欧美人与性动交α欧美软件| 国产精品久久久久久精品古装| 精品亚洲成国产av| 18禁国产床啪视频网站| 深夜精品福利| 精品国产乱码久久久久久男人| a级毛片在线看网站| 精品久久久久久电影网| 欧美精品av麻豆av| 亚洲欧美一区二区三区久久| 国产成人aa在线观看| 亚洲国产精品国产精品| 国产毛片在线视频| 香蕉丝袜av| 亚洲av欧美aⅴ国产| 日韩一卡2卡3卡4卡2021年| 男女国产视频网站| 久久午夜福利片| 五月天丁香电影| 女人被躁到高潮嗷嗷叫费观| 新久久久久国产一级毛片| 欧美变态另类bdsm刘玥| 久久女婷五月综合色啪小说| 91午夜精品亚洲一区二区三区| 日本色播在线视频| 欧美中文综合在线视频| 男女边吃奶边做爰视频| 天堂8中文在线网| 老司机亚洲免费影院| 激情视频va一区二区三区| 美女午夜性视频免费| 老汉色av国产亚洲站长工具| 永久免费av网站大全| 欧美人与性动交α欧美精品济南到 | 一级毛片黄色毛片免费观看视频| 国产成人精品久久二区二区91 | 在线天堂中文资源库| 精品福利永久在线观看| 精品国产一区二区久久| www.精华液| 一二三四中文在线观看免费高清| 人人妻人人添人人爽欧美一区卜| 午夜福利视频精品| 啦啦啦啦在线视频资源| 赤兔流量卡办理| 香蕉丝袜av| 看非洲黑人一级黄片| 青青草视频在线视频观看| 人妻少妇偷人精品九色| 国产亚洲精品第一综合不卡| 岛国毛片在线播放| 国产精品 国内视频| 桃花免费在线播放| 国产免费视频播放在线视频| a级毛片黄视频| 国产日韩欧美亚洲二区| 国产在线视频一区二区| 国产精品 国内视频| 日本av免费视频播放| 成人国产av品久久久| 考比视频在线观看| 在线精品无人区一区二区三| 永久网站在线| 久久久精品国产亚洲av高清涩受| 国产视频首页在线观看| 欧美人与性动交α欧美软件| 欧美成人午夜免费资源| 考比视频在线观看| 久久精品国产亚洲av高清一级| 少妇人妻久久综合中文| 综合色丁香网| 一边摸一边做爽爽视频免费| 97在线人人人人妻| 成人免费观看视频高清| videosex国产| 亚洲精品自拍成人| 日本猛色少妇xxxxx猛交久久| 中文字幕精品免费在线观看视频| 免费黄网站久久成人精品| 亚洲,欧美,日韩| 亚洲人成电影观看| 欧美精品国产亚洲| 午夜日韩欧美国产| 国产成人欧美| 熟女av电影| 只有这里有精品99| 国产精品一区二区在线不卡| 国产老妇伦熟女老妇高清| 日韩精品免费视频一区二区三区| 久久久久久久久久人人人人人人| 国产精品一国产av| 国产日韩一区二区三区精品不卡| 国产老妇伦熟女老妇高清| 在线观看一区二区三区激情| 最近中文字幕2019免费版| 中文欧美无线码| 国产成人av激情在线播放| 97精品久久久久久久久久精品| 巨乳人妻的诱惑在线观看| 黑丝袜美女国产一区| 久久人人爽av亚洲精品天堂| 欧美精品亚洲一区二区| 国产免费又黄又爽又色| 精品酒店卫生间| 久久精品久久久久久噜噜老黄| 久久综合国产亚洲精品| 久久精品夜色国产| 国产精品久久久av美女十八| 欧美激情高清一区二区三区 | 一本色道久久久久久精品综合| 91精品伊人久久大香线蕉| 亚洲国产欧美日韩在线播放| 日韩在线高清观看一区二区三区| 久久精品国产亚洲av高清一级| 亚洲成国产人片在线观看| 久久国内精品自在自线图片| 极品少妇高潮喷水抽搐| 国产成人精品在线电影| 亚洲美女视频黄频| av在线播放精品| 建设人人有责人人尽责人人享有的| 国产av一区二区精品久久| 亚洲精品aⅴ在线观看| 一本—道久久a久久精品蜜桃钙片| 91精品伊人久久大香线蕉| 欧美成人午夜免费资源| 国产白丝娇喘喷水9色精品| 久久久国产精品麻豆| 女人高潮潮喷娇喘18禁视频| 搡女人真爽免费视频火全软件| av在线老鸭窝| 日本猛色少妇xxxxx猛交久久| 日韩视频在线欧美| 中文字幕人妻丝袜一区二区 | 国产女主播在线喷水免费视频网站| 日韩三级伦理在线观看| av不卡在线播放| 免费黄频网站在线观看国产| 99国产综合亚洲精品| 美女国产视频在线观看| 又黄又粗又硬又大视频| 欧美精品高潮呻吟av久久| 久久精品国产亚洲av高清一级| 热re99久久国产66热| 丰满饥渴人妻一区二区三| 亚洲,欧美精品.| 亚洲精品日韩在线中文字幕| www.精华液| 国产精品99久久99久久久不卡 | 国产毛片在线视频| 欧美日韩精品成人综合77777| 啦啦啦中文免费视频观看日本| 如何舔出高潮| 免费在线观看完整版高清| 美女大奶头黄色视频| 日本色播在线视频| 丰满迷人的少妇在线观看| 一级黄片播放器| 久久久欧美国产精品| 黄色配什么色好看| 国产 一区精品| 大码成人一级视频| 一区二区三区四区激情视频| 免费看av在线观看网站| 中文字幕人妻丝袜一区二区 | av片东京热男人的天堂| 一本久久精品| 亚洲国产毛片av蜜桃av| 香蕉丝袜av| 国产成人精品婷婷| 人妻少妇偷人精品九色| 国产精品国产av在线观看| 成人手机av| 久久 成人 亚洲| 国产乱人偷精品视频| 午夜福利在线观看免费完整高清在| av网站在线播放免费| 97精品久久久久久久久久精品| 尾随美女入室| 欧美日韩一级在线毛片| 欧美人与性动交α欧美软件| 久久婷婷青草| 欧美亚洲 丝袜 人妻 在线| 天堂8中文在线网| 熟妇人妻不卡中文字幕| 精品国产一区二区三区久久久樱花| 在线观看三级黄色| 久久精品夜色国产| 高清视频免费观看一区二区| 美女主播在线视频| 色播在线永久视频| 免费少妇av软件| 一区二区三区精品91| 最近中文字幕高清免费大全6| 麻豆乱淫一区二区| 99国产精品免费福利视频| 亚洲av成人精品一二三区| 又黄又粗又硬又大视频| 赤兔流量卡办理| 男人操女人黄网站| 嫩草影院入口| 丝袜喷水一区| 这个男人来自地球电影免费观看 | av在线播放精品| 亚洲少妇的诱惑av| 日韩人妻精品一区2区三区| 国产一区二区 视频在线| 欧美日韩亚洲高清精品| 国产欧美亚洲国产| 国产97色在线日韩免费| 777久久人妻少妇嫩草av网站| 成人午夜精彩视频在线观看| 日本av免费视频播放| 色播在线永久视频| 男女无遮挡免费网站观看| 男女边摸边吃奶| 国产精品一区二区在线不卡| 香蕉精品网在线| 人体艺术视频欧美日本| 777米奇影视久久| 久久精品夜色国产| 91精品三级在线观看| 搡女人真爽免费视频火全软件| 婷婷色综合www| 欧美在线黄色| 亚洲精品日韩在线中文字幕| 啦啦啦中文免费视频观看日本| av线在线观看网站| 天堂中文最新版在线下载| 国产精品三级大全| 超碰97精品在线观看| 国产白丝娇喘喷水9色精品| 成人国产av品久久久| 永久免费av网站大全| 成人漫画全彩无遮挡| 丝袜在线中文字幕| 国产成人精品在线电影| 午夜激情久久久久久久| 亚洲欧美中文字幕日韩二区| 久久久欧美国产精品| 亚洲美女视频黄频| 五月天丁香电影| 国产成人一区二区在线| 成人国语在线视频| 久久精品国产综合久久久| 天天影视国产精品| 天天操日日干夜夜撸| 免费观看性生交大片5| av线在线观看网站| 中文天堂在线官网| 亚洲色图综合在线观看| 国产免费又黄又爽又色| 亚洲人成电影观看| 久久 成人 亚洲| 天天影视国产精品| 哪个播放器可以免费观看大片| av国产久精品久网站免费入址| 青春草视频在线免费观看| 99香蕉大伊视频| 中文字幕精品免费在线观看视频| 亚洲精品一二三| 日韩一卡2卡3卡4卡2021年| 最近的中文字幕免费完整| 看免费成人av毛片| 国产一区二区在线观看av| 午夜激情av网站| 性高湖久久久久久久久免费观看| 99久久中文字幕三级久久日本| 97精品久久久久久久久久精品| 国产免费现黄频在线看| 有码 亚洲区| 久久久久久人人人人人| 免费女性裸体啪啪无遮挡网站| 国产精品女同一区二区软件| 国产乱人偷精品视频| 18在线观看网站| av卡一久久| 亚洲精品一二三| 99热网站在线观看| 亚洲精品成人av观看孕妇| 精品国产露脸久久av麻豆| 欧美精品av麻豆av| 天天躁夜夜躁狠狠躁躁| 成人毛片a级毛片在线播放| 日日爽夜夜爽网站| 国产色婷婷99| 少妇人妻 视频| 精品少妇一区二区三区视频日本电影 | 久久99热这里只频精品6学生| 国产日韩欧美亚洲二区| 亚洲精品久久久久久婷婷小说| 丰满少妇做爰视频| 国产在线视频一区二区| 中国三级夫妇交换| 国产欧美日韩一区二区三区在线| 欧美精品人与动牲交sv欧美| 免费观看无遮挡的男女| 赤兔流量卡办理| 女人精品久久久久毛片| 一区二区三区激情视频| 欧美xxⅹ黑人| 久久97久久精品| 日日爽夜夜爽网站| 99久久人妻综合| 亚洲激情五月婷婷啪啪| 天天躁日日躁夜夜躁夜夜| 99国产精品免费福利视频| 日韩中文字幕视频在线看片| 亚洲av综合色区一区| 大话2 男鬼变身卡| 久久97久久精品| 国产精品不卡视频一区二区| 午夜激情久久久久久久| 99久久人妻综合| 免费人妻精品一区二区三区视频| 国产精品嫩草影院av在线观看| 老汉色av国产亚洲站长工具| 亚洲精品国产一区二区精华液| 国产成人欧美| 免费女性裸体啪啪无遮挡网站| 国产麻豆69| 九色亚洲精品在线播放| 亚洲成人av在线免费| 欧美日韩视频高清一区二区三区二| 九九爱精品视频在线观看| 国产黄频视频在线观看| 亚洲四区av| 久久久久久久久久人人人人人人| 青青草视频在线视频观看| 亚洲精品久久成人aⅴ小说| 搡女人真爽免费视频火全软件| 国产一级毛片在线| 日韩人妻精品一区2区三区| 性色avwww在线观看| 日本黄色日本黄色录像| 久久久精品94久久精品| 国产成人精品一,二区| 日韩精品有码人妻一区| 久久精品久久久久久久性| 大话2 男鬼变身卡| 毛片一级片免费看久久久久| 久久毛片免费看一区二区三区| 日韩三级伦理在线观看| 三上悠亚av全集在线观看| 两性夫妻黄色片| 九色亚洲精品在线播放| 高清黄色对白视频在线免费看| 成人黄色视频免费在线看| 日韩三级伦理在线观看| 国产爽快片一区二区三区| 中文欧美无线码| 97精品久久久久久久久久精品| 久久人妻熟女aⅴ| 亚洲成人手机| 看免费成人av毛片| 成人漫画全彩无遮挡| 在线观看免费视频网站a站| av免费在线看不卡| 久久人人97超碰香蕉20202| a级毛片在线看网站| 亚洲av免费高清在线观看| 大香蕉久久网| 国产成人av激情在线播放| av视频免费观看在线观看| av又黄又爽大尺度在线免费看| 纯流量卡能插随身wifi吗| 边亲边吃奶的免费视频| av卡一久久| 国产在线视频一区二区| 国产日韩一区二区三区精品不卡| 久久精品熟女亚洲av麻豆精品| 中文精品一卡2卡3卡4更新| 一个人免费看片子| 欧美日韩精品成人综合77777| 激情视频va一区二区三区| 欧美亚洲 丝袜 人妻 在线| 麻豆精品久久久久久蜜桃| 久久免费观看电影| 女性被躁到高潮视频| 午夜福利在线观看免费完整高清在| 国产成人a∨麻豆精品| 十八禁网站网址无遮挡| 欧美日韩视频高清一区二区三区二| 蜜桃在线观看..| 久久久久人妻精品一区果冻| 精品人妻在线不人妻| 中文字幕色久视频| 狠狠精品人妻久久久久久综合| 久久久久精品久久久久真实原创| 亚洲国产成人一精品久久久| 久久99蜜桃精品久久| 国产黄色视频一区二区在线观看| 天天操日日干夜夜撸| 老司机影院成人| 日本-黄色视频高清免费观看| 免费在线观看视频国产中文字幕亚洲 | 国产日韩一区二区三区精品不卡| 国产av一区二区精品久久| 国产激情久久老熟女| 精品一区二区三卡| av在线观看视频网站免费| 最近2019中文字幕mv第一页| 男的添女的下面高潮视频| 亚洲欧美中文字幕日韩二区| 大话2 男鬼变身卡| 在现免费观看毛片| 欧美xxⅹ黑人| 国产精品久久久久成人av| 丝袜美腿诱惑在线| videossex国产| 免费高清在线观看日韩| 久久久国产一区二区| 伊人久久大香线蕉亚洲五| 国产欧美日韩综合在线一区二区| 韩国高清视频一区二区三区| 最新的欧美精品一区二区| 777久久人妻少妇嫩草av网站| 国产成人一区二区在线| 亚洲av福利一区| 天天躁狠狠躁夜夜躁狠狠躁| 91aial.com中文字幕在线观看| 9191精品国产免费久久| 老汉色av国产亚洲站长工具| 久久99热这里只频精品6学生| 国产在线视频一区二区| 综合色丁香网| 亚洲国产日韩一区二区| 一二三四中文在线观看免费高清| 五月开心婷婷网| 高清不卡的av网站| 国产探花极品一区二区| 人人妻人人添人人爽欧美一区卜| 久久精品国产亚洲av涩爱| 精品少妇久久久久久888优播| 欧美人与善性xxx| 看免费av毛片| a级片在线免费高清观看视频| 夫妻午夜视频| 另类亚洲欧美激情| 国产高清国产精品国产三级| 国产伦理片在线播放av一区| 亚洲av电影在线观看一区二区三区| 午夜福利,免费看| 免费观看无遮挡的男女| 九草在线视频观看| 天美传媒精品一区二区| 另类亚洲欧美激情| 中国国产av一级| 人人妻人人爽人人添夜夜欢视频| 在线观看美女被高潮喷水网站| 人妻一区二区av| 亚洲国产精品一区三区| 午夜福利,免费看| 免费观看性生交大片5| 伦理电影大哥的女人| 男男h啪啪无遮挡| 亚洲人成网站在线观看播放| 9色porny在线观看| 少妇的逼水好多| 永久免费av网站大全| 国产精品一二三区在线看| 中文字幕人妻熟女乱码| 9色porny在线观看| 九草在线视频观看| 毛片一级片免费看久久久久| 在线精品无人区一区二区三| 观看美女的网站| 美女高潮到喷水免费观看| 国产成人a∨麻豆精品| 国产熟女欧美一区二区| 午夜福利影视在线免费观看| 成年美女黄网站色视频大全免费| 91在线精品国自产拍蜜月| 深夜精品福利| 在线观看免费高清a一片| av.在线天堂| 水蜜桃什么品种好| 国产黄频视频在线观看| 亚洲国产欧美日韩在线播放| 亚洲欧洲国产日韩| 国产亚洲精品第一综合不卡| 亚洲av福利一区| 国产又色又爽无遮挡免| 精品第一国产精品| 精品人妻熟女毛片av久久网站| 毛片一级片免费看久久久久| 中文字幕人妻丝袜制服| 亚洲欧美一区二区三区国产| 国产精品无大码| 精品亚洲乱码少妇综合久久| 中文精品一卡2卡3卡4更新| 日韩在线高清观看一区二区三区| 最近中文字幕高清免费大全6| 大陆偷拍与自拍| 欧美人与性动交α欧美软件| 纵有疾风起免费观看全集完整版| 一级黄片播放器| 日韩一本色道免费dvd| 老汉色∧v一级毛片| 国产精品秋霞免费鲁丝片| 青春草视频在线免费观看|