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

    一次近月制動的誤差分析與安全關(guān)機策略

    2015-03-13 07:01:16董捷周文艷張熇李驥關(guān)軼峰
    航天器工程 2015年4期
    關(guān)鍵詞:發(fā)動機

    董捷 周文艷 張熇 李驥 關(guān)軼峰

    (1 北京空間飛行器總體設(shè)計部, 北京 100094)(2 北京控制工程研究所, 北京 100190)(3 空間智能控制技術(shù)重點實驗室, 北京 100190)

    ?

    一次近月制動的誤差分析與安全關(guān)機策略

    董捷1周文艷1張熇1李驥2,3關(guān)軼峰2,3

    (1 北京空間飛行器總體設(shè)計部, 北京 100094)(2 北京控制工程研究所, 北京 100190)(3 空間智能控制技術(shù)重點實驗室, 北京 100190)

    實施有效的近月制動進(jìn)入目標(biāo)環(huán)月軌道是實現(xiàn)月球環(huán)繞、安全軟著陸等任務(wù)的前提。文章基于一次近月制動,對點火過程中的主要誤差要素影響進(jìn)行了分析,包括開機與關(guān)機方式、點火姿態(tài)、有限推力制動方式的影響,其中點火姿態(tài)偏差是最大的誤差影響。并分析了一次近月制動條件下,由于推力偏差和制動推質(zhì)比較大,速度增量關(guān)機失效時自主時間關(guān)機和地面支持關(guān)機存在的軌道安全性問題。針對此問題,按近月制動前后的不同階段設(shè)計了3種安全關(guān)機策略,包括相關(guān)的發(fā)動機標(biāo)定、軌控實時監(jiān)視干預(yù)以及應(yīng)急提升近月點。最后介紹了安全關(guān)機策略在嫦娥三號任務(wù)中的應(yīng)用情況。

    近月制動;誤差;點火姿態(tài);發(fā)動機關(guān)機

    1 引言

    月球探測器典型的探測形式包括飛掠、環(huán)繞、著陸、巡視和采樣返回。飛掠與直接著陸通常集中在早期月球探測任務(wù)中,如蘇聯(lián)的月球1號飛掠任務(wù)[1]、蘇聯(lián)的月球2號硬著陸任務(wù)[1]、美國的勘測者號[2]軟著陸任務(wù)等。由于飛掠與直接硬著陸探測成果有限,而直接軟著陸方式對轉(zhuǎn)移軌道發(fā)射窗口要求嚴(yán)格、著陸區(qū)受限[3],所以探測方式逐漸發(fā)展為先在近月點附近執(zhí)行制動完成月球環(huán)繞飛行,再開展后續(xù)遙感、軟著陸等任務(wù)。

    根據(jù)制動推力大小,從降低重力損耗角度考慮,包括一次制動或多次制動。當(dāng)制動初始推質(zhì)比(發(fā)動機推力與探測器質(zhì)量之比)能夠達(dá)到2 m/s2以上時,通常采用一次制動,這類任務(wù)集中在月球軟著陸探測。由于配置了大推力減速發(fā)動機,可以一次制動直接進(jìn)入環(huán)月目標(biāo)軌道,簡化飛行控制程序,如美國的阿波羅系列載人飛船[4]、蘇聯(lián)的“月球”(Luna)系列無人探測器、中國的嫦娥三號、處于規(guī)劃和研制階段的俄羅斯“月球-全球”探測器[5]等。當(dāng)制動初始推質(zhì)比在2 m/s2以內(nèi)時,則采用多次制動方式逐漸降低遠(yuǎn)月點高度,這類任務(wù)通常針對環(huán)繞月球遙感探測,如美國的“月球勘測軌道器”(LRO)、中國的嫦娥一號、二號。

    對于近月制動,如不能及時進(jìn)行減速控制,探測器將飛離地月系統(tǒng)不再返回,若要形成較為穩(wěn)定的環(huán)月軌道,對應(yīng)的減速速度增量至少要達(dá)到200~300 m/s;如減速過大,將有撞月風(fēng)險[6];對于制動前近月點高度為100 km的軌道,如果速度增量超出制動所需速度增量20 m/s,控后軌道近月點將接近月面。這決定了近月制動窗口具有唯一性,其控制時機及控制量直接影響任務(wù)成敗,因此實現(xiàn)近月制動精確控制至關(guān)重要。本文針對一次近月制動中的各項控制誤差影響進(jìn)行了分析,提出了應(yīng)對安全關(guān)機風(fēng)險的策略,最后介紹了嫦娥三號任務(wù)的在軌應(yīng)用情況。

    2 一次近月制動的誤差源分析

    2.1 近月制動控制目標(biāo)

    近月制動后的核心目標(biāo)是保證探測器的軌道高度、軌道傾角滿足任務(wù)要求,在故障情況下能確保軌道的安全性。對于采用一次制動實現(xiàn)的月球軟著陸任務(wù),對近月制動后的軌道高度要求并不嚴(yán)格,一般不超過±50 km。保證軌道傾角的目的,是為了實現(xiàn)月面目標(biāo)區(qū)域的可達(dá)性,只要近月制動后軌道傾角大于著陸區(qū)緯度[7],就可以到達(dá)目標(biāo)著陸區(qū)。只要著陸區(qū)范圍不是嚴(yán)格的月球兩極區(qū)域,綜合軌控精度、環(huán)月軌道攝動對環(huán)月傾角的影響,選擇極軌軌道對兩極區(qū)域具有可達(dá)性,軌道傾角誤差允許范圍通常不大于±5°。

    以2017年5月發(fā)射窗口為例,開展近月制動各種誤差影響分析。選擇標(biāo)稱地月轉(zhuǎn)移軌道轉(zhuǎn)移時間為112 h,到達(dá)月球時的近月點高度為100 km。在地心(第一)赤道坐標(biāo)系[8]下,與運載火箭分離時的標(biāo)稱器箭分離參數(shù)見表1。

    表1 地月轉(zhuǎn)移軌道器箭分離參數(shù)

    本文假設(shè)研究的環(huán)月軌道在月心坐標(biāo)系下主要參數(shù)及偏差如下。

    (1)近月制動后目標(biāo)軌道:半長軸1 837.4 km,相對于月球平均半徑(1 737.4 km)的高度為100 km,制動后允許軌道偏差±30 km(不含最后一次中途修正的誤差影響、近月制動前的軌道預(yù)報誤差);

    (2)制動后目標(biāo)軌道傾角90°,允許偏差±5°;

    (3)近月制動脈沖速度增量843 m/s;

    (4)探測器發(fā)動機標(biāo)稱推力6000 N,推力偏差3%;

    (5)制動前探測器初始質(zhì)量3000 kg,質(zhì)量預(yù)計偏差0.5%;

    (6)發(fā)動機比沖為308 s,比沖偏差取1%。

    對應(yīng)近月制動時探測器的初始推質(zhì)比達(dá)到了2 m/s2,可以采用一次制動策略,這也與阿波羅飛船、嫦娥三號探測器制動初始推質(zhì)比接近。

    2.2 近月制動開關(guān)機影響

    近月制動過程中,發(fā)動機開、關(guān)機時機設(shè)計尤為重要,一方面關(guān)系到制動點火所需速度增量,會影響推進(jìn)劑消耗;另一方面影響控后軌道,如設(shè)計不合理,將造成控后無法環(huán)月或近月點高度過低發(fā)生撞月。

    1)發(fā)動機開機

    在執(zhí)行軌道控制前,地面需提前注入探測器上軌道控制參數(shù)(包括點火開始時刻、速度增量、目標(biāo)姿態(tài)等),探測器上判斷點火時刻到后自主發(fā)出發(fā)動機開機指令。

    為提高有限推力制動的效率,推力弧段中心時刻選擇為地月轉(zhuǎn)移軌道近月點時刻。根據(jù)環(huán)月目標(biāo)軌道、發(fā)動機推力、比沖和探測器質(zhì)量進(jìn)行制動過程有限推力計算,確定點火開始時刻和速度增量。

    發(fā)動機開機的準(zhǔn)確度與星時精度、星上的控制周期相關(guān),通常綜合誤差在百毫秒量級以下。反映了點火中心時刻偏離近月點時刻的偏差。

    2)發(fā)動機關(guān)機

    發(fā)動機的關(guān)機方式包括3種,即由探測器自主執(zhí)行的速度增量關(guān)機、時間關(guān)機,以及在地面判斷關(guān)機無效時的地面支持關(guān)機。為了保證變軌精度,正常條件下是依靠速度增量實現(xiàn)發(fā)動機關(guān)機,并以時間關(guān)機作為加速度計異常等特殊條件下的關(guān)機手段。地面支持關(guān)機指的是在地面判斷出發(fā)動機未正常關(guān)機后,由地面發(fā)送關(guān)機指令。

    (1)速度增量關(guān)機:探測器上根據(jù)加速度計測量數(shù)據(jù)計算速度增量變化,當(dāng)判斷速度增量達(dá)到要求后,自主發(fā)出發(fā)動機關(guān)機指令。這種關(guān)機方式,精度較高,沿推力方向的誤差大小主要包括加速度計誤差、傳輸時延(控制周期量化及串口傳輸時延)和關(guān)機后效三部分。對應(yīng)的誤差模型如下:

    (1)

    (2)時間關(guān)機:一般在加速度計故障等異常條件下,作為速度增量關(guān)機手段的備份。變軌過程中,如果依據(jù)速度增量關(guān)機未正常執(zhí)行,探測器上計算到達(dá)允許的最長點火時間后,將自主發(fā)出發(fā)動機關(guān)機指令。最長點火時間設(shè)計過小,可能在正常速度增量關(guān)機前發(fā)出關(guān)機指令,增大變軌偏差;設(shè)計過大,當(dāng)速度增量關(guān)機失效時,無法保證在“撞月”前發(fā)出關(guān)機指令。時間關(guān)機誤差與加速度計誤差無關(guān),而與發(fā)動機推力偏差直接相關(guān)。則對應(yīng)的誤差模型如下:

    Et=F/m×(ΔTf+td)

    (2)

    式中:Et為自主時間關(guān)機總的速度增量誤差,ΔTf為計算的最長點火時間與實際點火時間之差。

    如果要保證探測器自主時間關(guān)機有效,ΔTf至關(guān)重要。最長點火時間的計算公式如下:

    (3)

    (4)

    式中:Tf為最長點火時間,m0為探測器初始質(zhì)量,mf為點火過程中推進(jìn)劑消耗情況,Isp為比沖,ΔV為速度增量。

    則在ΔV相對固定條件下,最長點火時間與初始質(zhì)量、推力、比沖相關(guān),Tf可以表示為

    Tf=f(m0,F,Isp)

    (5)

    在m0、F、Isp偏差較小時,相對于標(biāo)稱點火時間的偏差可采用如下公式計算:

    (6)

    式中:Δm0、ΔF、ΔIsp分別為m0、F和Isp的偏差量。

    因此,預(yù)計的探測器初始質(zhì)量偏差、發(fā)動機推力偏差和比沖偏差都對預(yù)計的最長點火時間有影響。為了提高自主時間關(guān)機的精度,必須盡可能縮小m0、F和Isp這3個量的偏差值。

    (3)地面支持關(guān)機:由地面介入的備份手段,以進(jìn)一步提高關(guān)機可靠性,保證軌道安全,不發(fā)生撞月。為保證地面支持關(guān)機有效,在經(jīng)歷地面遙測判讀、處置及地面上行執(zhí)行的總時間(設(shè)為ΔTz)后,其帶來的速度增量誤差應(yīng)不大于軌道安全所允許的極限值。對應(yīng)的誤差模型如下:

    Ez=F/m×ΔTz+F/m×td

    (7)

    式中:Ez為地面支持關(guān)機總的速度增量誤差。

    由于ΔTz相對固定,在正常關(guān)機時刻后探測器推質(zhì)比越小,制動速度增量增加越慢,地面支持關(guān)機越能最大程度地發(fā)揮作用。

    2.3 點火姿態(tài)影響

    點火姿態(tài)的設(shè)計包括沿推力方向和垂直于推力方向兩部分。直接影響變軌效果的是沿推力方向,垂直于推力方向用于保證星敏感器的正常使用(姿態(tài)基準(zhǔn)丟失時可以快速獲取慣性姿態(tài))。

    目前,探測器具備在軌控時推力方向沿慣性空間固定和勻速轉(zhuǎn)動兩種能力[9]。固定推力方向較勻速轉(zhuǎn)動制動速度增量大,但對于大推質(zhì)比制動的探測器,速度增量相差僅幾米每秒。而較固定推力控制,勻速轉(zhuǎn)動控制會增大軌道高度偏差,降低軌道安全性,并增加后續(xù)進(jìn)入目標(biāo)軌道修正量;勻速轉(zhuǎn)動控制時探測器慣性姿態(tài)不斷變化,控制策略復(fù)雜,在軌監(jiān)視時也不易進(jìn)行故障診斷。為提高可靠性和降低復(fù)雜度,月球軟著陸探測器近月制動通常設(shè)計采用固定推力方向控制。為提高制動效率,慣性空間的推力方向選擇為地月轉(zhuǎn)移軌道近月點時刻速度的反方向。

    定義標(biāo)稱點火姿態(tài)下,探測器質(zhì)心坐標(biāo)系三軸指向如下:探測器X軸沿標(biāo)稱推力方向,即X軸沿地月轉(zhuǎn)移軌道近月點時刻速度的反方向,Y軸垂直于X軸,與軌道面法線方向一致,Z軸與X軸和Y軸構(gòu)成右手坐標(biāo)系。

    由于繞推力方向的角度誤差不影響制動結(jié)果,因此模型可以定義為沿垂直于推力方向的兩個軸的角度誤差。設(shè)繞Y軸的角度誤差為α,繞Z軸的角度誤差為β。

    根據(jù)前面分析條件,標(biāo)稱點火姿態(tài)四元數(shù)見表2(相對于J2000月心坐標(biāo)系)[8]。

    表2 點火姿態(tài)(慣性空間定向)

    其中,四元數(shù)定義[10]為:設(shè)特征軸單位矢量為e,轉(zhuǎn)角為φ,四元數(shù)q=(q0,q1,q2,q3)定義如下

    (8)

    3 誤差影響仿真分析

    3.1 固定推力方向的控后軌道

    根據(jù)前面計算的輸入條件標(biāo)稱值,計算的控后軌道參數(shù)為:

    (1)點火時間:369.3 s;

    (2)總速度增量:846.9 m/s;

    (3)推進(jìn)劑消耗量:733.6 kg。

    (4)相對于平均月球半徑的近月點高度:92.5 km,遠(yuǎn)月點高度:107.5 km。

    (5)控后探測器質(zhì)量:2 266.4 kg。

    (6)重力損耗;(846.9 m/s-843 m/s)/843 m/s=0.46%。

    可見,在不含其它誤差條件,點火方向沿近月點反方向的點火策略下,由于有限推力控制,近月點高度誤差為7.5 km。

    由于近月制動推質(zhì)比大,實際變軌重力損耗相對于脈沖變軌比例小于1%,一次制動的推進(jìn)劑代價很小。

    3.2 控后軌道的誤差分析

    3.2.1 開機延遲影響

    目前,星載計算機的控制周期通常不大于160 ms,星時誤差不大于20 ms,開機時延不超過20 ms(信號傳輸時間10 ms、閥門動作時間10 ms),則總的開機時間延遲最長為200 ms。在其它參數(shù)沒有誤差的條件下,對應(yīng)的控后軌道高度為近月點92.53 km,遠(yuǎn)月點107.05 km。

    可見在延遲200 ms開機的條件下,對控后軌道的影響相對較小,高度誤差小于1 km。

    3.2.2 速度增量關(guān)機誤差影響

    加速度計測量精度與推質(zhì)比大小相關(guān),對于初始推質(zhì)比為2 m/s2,典型的加速度計測量精度優(yōu)于1×10-3m/s2;傳輸時延包括控制周期160 ms、信號傳輸時間10 ms、閥門動作時間10 ms,則根據(jù)前面介紹的誤差模型,計算的速度增量誤差Ev=0.848 2 m/s。則對應(yīng)的控后軌道參數(shù)為近月點90.09 km,遠(yuǎn)月點105.73 km。

    可見在正常速度增量關(guān)機誤差條件下,對控后軌道的影響也相對較小,近月點高度誤差小于2.5 km,遠(yuǎn)月點高度誤差小于2 km。

    3.2.3 姿態(tài)控制誤差影響

    根據(jù)前面姿態(tài)偏差的定義,分別對兩個軸方向誤差為1°,2°和3°的情況進(jìn)行分析(見表3)。

    表3 姿態(tài)控制誤差對制動后軌道的影響

    注:上述姿態(tài)偏差按相對保守的方法取單方向固定值,實際變軌姿態(tài)偏差會在兩個方向周期性變化。

    可見,隨著軌道平面內(nèi)的姿態(tài)角偏差α不斷增大,軌道近月點高度將快速降低,α角誤差在1°時,軌道高度誤差就達(dá)到25 km;這主要是由于平面內(nèi)誤差影響軌道偏心率,造成軌道高度偏差增大,特別是近月點高度迅速降低,會危及軌道安全。平面外姿態(tài)偏差β主要影響軌道傾角,由于軌控速度增量偏差與探測器相對月球的運動速度是小量,因此合成后對軌道傾角的影響相對較小。

    綜合以上分析,軌控時的姿態(tài)控制誤差應(yīng)該控制在1°以內(nèi)。

    3.3 非正常關(guān)機的近月制動后軌道

    3.3.1 自主時間關(guān)機誤差影響

    根據(jù)前面的分析,自主時間關(guān)機速度增量誤差Et=36.2 m/s。

    不同延遲關(guān)機時間下對應(yīng)的軌道近月點高度變化見表4??梢?,自主時間關(guān)機速度增量誤差明顯大于安全高度允許的速度增量,這樣自主時間關(guān)機無法發(fā)揮作用。只有減小發(fā)動機推力偏差才能提高自主時間關(guān)機的有效性。

    表4 延時關(guān)機時間與軌道近月點高度的變化

    3.3.2 地面支持關(guān)機的誤差影響

    與自主時間關(guān)機相似,按照允許延遲關(guān)機1 min進(jìn)行分析,計算的地面支持關(guān)機速度增量誤差Ez=159.3 m/s。由表4可見,顯然無法滿足近月點高度要求。經(jīng)分析,只有對應(yīng)延遲關(guān)機后的推質(zhì)比小于0.3 m/s2時,才能滿足地面支持關(guān)機的安全性要求。

    針對月球環(huán)繞探測任務(wù),如我國的嫦娥一號[11]、嫦娥二號[12]和美國的LRO[13],這些探測器推質(zhì)比小,只能采用多次近月制動降低重力損耗;但恰恰由于該特點,延遲關(guān)機后至少有1 min以上的時間近月點才處于月球平均半徑以下,地面有相對充足的應(yīng)急處置時間,關(guān)機安全性的風(fēng)險相對較低。

    對于月球軟著陸這類大推質(zhì)比任務(wù),由于延遲關(guān)機后近月點將快速處于月球以下,地面是無法支持再次關(guān)機的。美國阿波羅飛船這類載人任務(wù)從軌道安全性考慮采用了地月自由返回軌道,近月制動時只能處在月球背面,受月球遮擋地面無法支持關(guān)機,但可以由飛船上宇航員實時監(jiān)視支持[14];而對于蘇聯(lián)Luna任務(wù),目前公開的文獻(xiàn)并沒有針對關(guān)機可靠性問題給出明確的策略,只能依賴速度增量關(guān)機。

    3.4 誤差影響分析

    綜合以上誤差分析,在正常速度增量關(guān)機方法下,在開機延遲、速度增量關(guān)機誤差、點火姿態(tài)三項因素中,點火姿態(tài)影響最大,因此需要將姿態(tài)控制誤差控制在1°以內(nèi)。綜合開機延遲、速度增量誤差、有限推力點火制動以及點火姿態(tài)誤差后的軌道誤差如下,滿足要求的控后最大30 km的偏差要求。

    (9)

    式中:ΔEh1為有限推動力制動產(chǎn)生的高度誤差,ΔEh2為開機延遲產(chǎn)生的高度誤差,ΔEh3為速度增量關(guān)機誤差產(chǎn)生的高度誤差,ΔEh4為點火姿態(tài)誤差(1°)產(chǎn)生的高度誤差。

    在速度增量關(guān)機失效的情況下,受較大的推力偏差和制動推質(zhì)比影響,自主時間關(guān)機與地面支持關(guān)機都無法保證軌道安全,因此必須要設(shè)計相應(yīng)的安全關(guān)機策略。

    4 安全關(guān)機策略設(shè)計

    針對自主時間關(guān)機和地面支持關(guān)機存在的問題,可以按階段在近月制動前、中、后期設(shè)計三種策略以提高關(guān)機安全性,避免在速度增量關(guān)機失效的情況下發(fā)生撞月。

    4.1 近月制動前制動發(fā)動機的標(biāo)定

    設(shè)探測器正常關(guān)機時刻至距月面安全高度時刻的時間為T。為使自主時間關(guān)機時,探測器近月點仍在安全高度以上,要求根據(jù)推力偏差計算的最大和最小點火時間差(設(shè)為ΔT)小于T,即以最大推力按最長時間點火后仍能使近月點處于安全高度以上。

    表5對影響自主時間關(guān)機的三個因素進(jìn)行了分析,主要因素是發(fā)動機推力大小,準(zhǔn)確獲取發(fā)動機的常值推力偏差可以有效縮短ΔT的大小。受地面對發(fā)動機推力標(biāo)定手段的限制,地面無法通過熱標(biāo)試驗準(zhǔn)確標(biāo)定在軌發(fā)動機的推力大小,存在較大的常值偏差,通常都在±3%以上。因此近月制動前的策略是利用地月轉(zhuǎn)移過程中途修正的機會,對用于近月制動的大推力發(fā)動機推力大小進(jìn)行標(biāo)定。

    表5 自主時間關(guān)機不同誤差對最大點火時間的影響

    隨著運載火箭制導(dǎo)能力的提升,深空探測器的入軌精度不斷提高,所需要的中途修正量可能很?。欢鴺?biāo)定大推力發(fā)動機,需要使發(fā)動機推力處于穩(wěn)定工作階段,這樣點火時間就需要盡可能長,考慮大推力發(fā)動機流量大,必然會帶來推進(jìn)劑的大量消耗。因此需要在兼顧中途修正需求與優(yōu)化推進(jìn)劑消耗兩方面進(jìn)行平衡開展發(fā)動機推力標(biāo)定策略設(shè)計。

    策略的首要目標(biāo)是先保證發(fā)動機的最短標(biāo)定時間,在此基礎(chǔ)上設(shè)計各次中途修正點火策略,具體如下:

    (1)當(dāng)?shù)?次或第2次中途修正滿足大推力發(fā)動機最短標(biāo)定時間要求時,直接采用大推力發(fā)動機進(jìn)行修正。

    (2)當(dāng)中途修正速度增量小于發(fā)動機最短標(biāo)定時間時,進(jìn)行兩次中途修正組合設(shè)計,即第1次中途修正使用推力器點火進(jìn)行小速度增量軌道拉偏,滿足第2次修正采用大推力發(fā)動機的最短標(biāo)定時間要求。

    通過以上優(yōu)化修正策略,在滿足大推力發(fā)動機點火標(biāo)定最短時間的條件下,使修正速度增量最小。在近月制動前獲得發(fā)動機推力標(biāo)定結(jié)果后,再考慮較小的隨機偏差,計算近月制動最大點火時長并注入探測器上執(zhí)行。

    4.2 近月制動中判斷評估撞月風(fēng)險

    該策略的基本思路是在近月制動軌控中,可以由地面判斷軌控工作狀態(tài)進(jìn)行實時干預(yù),具體策略為:

    (1)首先根據(jù)發(fā)動機最大的推力偏差得到最短點火時間和最長點火時間,結(jié)合點火時刻得到預(yù)期的最早關(guān)機時刻和最晚關(guān)機時刻。

    (2)在點火過程中,探測器上根據(jù)實時計算的剩余點火速度增量等信息,自主計算剩余點火時間。在點火開始一定時間后,根據(jù)已知的點火時刻實時預(yù)測關(guān)機時刻是否處于最早關(guān)機時刻與最晚關(guān)機時刻之間的預(yù)期范圍,如不處于該范圍,說明發(fā)動機推力異?;蛘呤羌铀俣扔嬢敵鼋Y(jié)果異常。這時地面可以發(fā)送指令,實現(xiàn)在指定時刻(即最早關(guān)機時刻)向探測器上注入關(guān)機指令提前關(guān)閉發(fā)動機,以保證探測器處于安全軌道高度。如果關(guān)機時刻始終處于預(yù)期范圍,地面則不加干預(yù),由探測器上自主實現(xiàn)正常的速度增量關(guān)機。

    4.3 近月制動后發(fā)現(xiàn)撞月風(fēng)險時的控制策略

    該策略的基本思路是在正常軌控點火時間之后,快速調(diào)整姿態(tài),調(diào)整點火方向以加速抬高軌道。

    在地面發(fā)現(xiàn)點火時間超過預(yù)期,若僅發(fā)送關(guān)機指令,近月點高度仍將處于月面以下,因此需要盡快調(diào)整姿態(tài)至軌道加速方向,以抬高近月點高度。軟著陸任務(wù)探測器上為了平衡大推力造成的干擾力矩,通常需要推力器和力臂的組合設(shè)計提供足夠的姿態(tài)控制力矩,這種方式可以使姿態(tài)調(diào)整角速度達(dá)到幾度/秒,為及時處置提供了有利的配置保障。

    具體策略如下:

    近月制動時如發(fā)動機未實現(xiàn)正常關(guān)機,制動點火速度增量超出所需速度增量,探測器軌道近月點將低于月面,必須緊急進(jìn)行故障處置。具體故障處置的方法是關(guān)機后根據(jù)探測器上遙測“剩余速度增量(ΔVs)”進(jìn)行是否處置的判斷。

    (1)如額外增加的速度增量大于允許閾值(ΔV)則必須立即進(jìn)行應(yīng)急處置,通過調(diào)整探測器姿態(tài),使探測器由制動減速姿態(tài)調(diào)整為加速姿態(tài)進(jìn)行變軌,提高軌道近月點高度。開機變軌越早撞月風(fēng)險越小,所需的速度增量代價也越小,因此必須盡快完成處置。

    (2)如額外增加的速度增量不大于ΔV,近月點高度滿足安全高度,由于在一天內(nèi)月球攝動對軌道高度的影響不危及安全性,制動后不進(jìn)行應(yīng)急故障處置,進(jìn)行正常軌控后的狀態(tài)恢復(fù)和設(shè)置。并根據(jù)控后定軌數(shù)據(jù)判斷軌道偏差,選擇是否安排第二天進(jìn)行軌道修正。

    5 在軌驗證

    嫦娥三號在軌采用了軟著陸使用的7500 N大推力發(fā)動機進(jìn)行近月制動,對一次近月制動策略進(jìn)行了驗證。

    1)近月制動前

    為了標(biāo)定7500 N發(fā)動機在軌推力,確保近月制動安全,嫦娥三號在飛行程序中安排了用7500 N發(fā)動機進(jìn)行中途修正的飛行事件。探測器發(fā)射入軌精度高,因此于2013年12月2日進(jìn)行了第一次中途修正軌道拉偏,于2013年12月3日第二次中途修正時采用了7500 N推力,修正時長為標(biāo)定所需的最短時間,對穩(wěn)態(tài)推力進(jìn)行了精確標(biāo)定,修正了發(fā)動機地面熱標(biāo)數(shù)據(jù),使近月制動時間關(guān)機參數(shù)的設(shè)置更為準(zhǔn)確。2013年12月6日,嫦娥三號完美完成近月制動。近月制動實際開機時長與標(biāo)稱時長差小于2 s,證明了推力標(biāo)定的有效性。

    2)近月制動過程中

    點火過程中實時預(yù)測的點火時刻始終處于允許誤差范圍內(nèi),如圖1所示,因此未執(zhí)行軌控過程中特定時刻關(guān)機的策略??睾筌壍腊腴L軸1 835.12 km(目標(biāo)半長軸1 835.08 km)、軌道傾角89.92°(目標(biāo)傾角90°),滿足后續(xù)著陸任務(wù)需求。點火過程中姿態(tài)控制精度優(yōu)于0.8°,地面對點火速度增量中的后效誤差進(jìn)行了有效補償,實現(xiàn)關(guān)機精度優(yōu)于0.1 m/s,在考慮近月點軌道預(yù)報誤差后,近月制動綜合的軌道高度偏差不大于10 km。

    圖1 近月制動地面關(guān)機時刻監(jiān)視Fig.1 Monitoring of engine shutdown under braking at perilune

    由于近月制動順利完成,4.3節(jié)的應(yīng)急抬高軌道控制不需要執(zhí)行。但為了保證對該策略進(jìn)行充分驗證,任務(wù)前經(jīng)過了全數(shù)字仿真分析,在此基礎(chǔ)上探測器系統(tǒng)與北京飛行控制中心還開展了多次故障協(xié)同演練。利用地面動力學(xué)計算機仿真在軌動力學(xué)環(huán)境,制導(dǎo)導(dǎo)航與控制(GNC)、數(shù)管等關(guān)鍵分系統(tǒng)采用與星上設(shè)備相同的軟硬件產(chǎn)品開展了1∶1閉環(huán)模飛,充分驗證了對策的有效性。

    6 結(jié)束語

    本文對一次近月制動過程中的各項軌道控制誤差進(jìn)行了分析,歸納了最為關(guān)鍵的誤差源,提出了對主要誤差源點火姿態(tài)的要求。針對速度增量關(guān)機失效條件下,由于推力偏差大及推質(zhì)比大帶來的一次近月制動關(guān)機風(fēng)險,分別在近月制動前、中、后設(shè)計了相應(yīng)策略。近月制動前通過中途修正進(jìn)行推力標(biāo)定,使自主時間關(guān)機參數(shù)設(shè)置更為準(zhǔn)確;近月制動過程中設(shè)計實時預(yù)測并評估判斷點火關(guān)機時刻的方法,使地面有手段在制動過程中及時干預(yù);根據(jù)發(fā)動機地面熱標(biāo)和第一次中途修正的數(shù)據(jù),通過減小點火速度增量補償了關(guān)機后效,從而進(jìn)一步降低了控制誤差;在近月制動后,充分利用了自身大推力器可提供大角速度的優(yōu)勢,設(shè)計調(diào)姿加速對策防止延遲關(guān)機引起的近月點高度快速降低。

    該策略經(jīng)過了嫦娥三號的在軌驗證,在國內(nèi)首次采用一次近月制動進(jìn)入目標(biāo)軌道,優(yōu)化了飛控程序,為圓滿完成任務(wù)奠定了堅實基礎(chǔ)。

    References)

    [1]郗曉寧,曾國強,任萱,等.月球探測器軌道設(shè)計[M].北京:國防工業(yè)出版社,2001:23-34

    Xi Xiaoning, Zeng Guoqiang, Ren Xuan, et al. Orbit design of lunar probe[M]. Beijing: National Defense Industry Press,2001:23-34 (in Chinese)

    [2]褚桂柏,張熇.月球探測器技術(shù)[M].北京:中國科學(xué)技術(shù)出版社,2007:128-129

    Chu Guibai,Zhang He.Lunar probe technology[M].Beijing: China Science and Technology Press,2007:128-129 (in Chinese)

    [3]楊嘉墀.航天器軌道動力學(xué)與控制(下)[M].北京:中國宇航出版社,2005:348

    Yang Jiachi.Spacecraft trajectory dynamics and control[M].Beijing:China Astronautics Press, 2005:348 (in Chinese)

    [4]Interbartolo, Michael A.Apollo mission techniques lunar orbit activities-Part 1a, NASA 20090015331[R].Washington D.C.:NASA,2009

    [5]魯暘筱懿,平勁松.V Shevchenko.俄羅斯“月球-全球”和“月球-資源”探月任務(wù)[J].航天器工程,2013,22(4): 103-108 (in Chinese)

    Lu Yangxiaoyi, Ping Jinsong,V Shevchenko.Russian lunar exploration missions:LUNA-Glob & LUNA-Resource[J]. Spacecraft Engineering, 2013,22(4): 103-108 (in Chinese)

    [6]吳偉仁,周建亮,高薇.繞月探測衛(wèi)星飛行控制[M].北京:中國宇航出版社,2012:134-135

    Wu Weiren,Zhou Jianliang,Gao Wei. Flight control of lunar orbiting satellite[M].Beijing:China Astronautics Press, 2012:134-135 (in Chinese)

    [7]章仁為.衛(wèi)星軌道姿態(tài)動力學(xué)與控制[M].北京:北京航空航天大學(xué)出版社,1998:18

    Zhang Renwei. Dynamics and control of satellite orbit and attitude [M].Beihang University Press,1998:18 (in Chinese)

    [8]張旭輝,余夢倫,茹家欣.QJ 1028B-2008 航天產(chǎn)品常用坐標(biāo)系[S].北京:國防科學(xué)技術(shù)工業(yè)委員會,2008

    Zhang Xuhui,Yu Menglun,Ru Jiaxin. QJ 1028B-2008 General coordinate for astronautics product[S].Beijing:Commission of Science,Technology and Industry for National Defense,2008 (in Chinese)

    [9]李鐵壽.一種減少航天器重力損失的變軌方法:中國,200710301743.2[P].2008-05-28

    Li Tieshou.A maneuver methord for reducing gravity losses:China, 200710301743.2[P].2008-05-28 (in Chinese)

    [10]譚維熾,胡金剛.航天器系統(tǒng)工程[M].北京:中國科學(xué)技術(shù)出版社,2009:44

    Tan Weichi,Hu Jingang.Spacecraft systems engineering[M].Beijing:China Science and Technology Press, 2009:44 (in Chinese)

    [11]楊維廉,周文艷.嫦娥一號月球探測衛(wèi)星軌道設(shè)計[J].航天器工程,2007,16(6):16-24

    Yang Weilian,Zhou Wenyan. Chang’e-1 lunar trajectotry design[J].Spacecraft Engineering,2007,16(6):16-24 (in Chinese)

    [12]周文艷,楊維廉.嫦娥二號衛(wèi)星軌道設(shè)計[J].航天器工程,2010,19(5):24-28

    Zhou Wenyan, Yang Weilian. Chang’e-2 lunar trajectotry design[J].Spacecraft Engineering,2010,19(5):24-28 (in Chinese)

    [13]Ann Nicholson. NASA GSFC Lunar Reconnaissance Orbiter (LRO) orbit estimation and prediction,AIAA 2010-2328[R].Washington D.C.:AIAA,2010

    [14]Berry R L,Wiley R F.Real-time targeting for the Apollo lunar orbit insertion maneuver, NASA-TM-X-69759[R].Washington D.C.:AIAA,1968

    (編輯:李多)

    Error Analysis and Safe Shutdown Strategy Design for One Braking at Perilune

    DONG Jie1ZHOU Wenyan1ZHANG He1LI Ji2,3GUAN Yifeng2,3

    (1 Beijing Institute of Spacecraft System Engineering, Beijing 100094, China (2 Beijing Institute of Control and Engineering, Beijing 100190, China) (3 Science and Technology on Space Intelligent Control Laboratory, Beijing 100190, China)

    For critical missions such as orbiting exploration of and soft landing on the moon, a braking at perilune is the prerequisite. This paper analyzes the main error factors, including firing and shutdown forms of the engine, firing attitude, and influences of finite thrust braking on ignition process under one braking at perilune condition. By comparison and analysis, it concludes that the firing attitude error is the biggest error source. It also indicates the potential safety hazard of the autonomous delay shutdown and ground control shutdown process due to the big thrust bias and thrust-mass ratio of braking, when velocity increment shutdown fails. The safe shutdown strategy in different stages of braking at perilune is designed to resolve the problem mentioned above, including engine calibration, maneuver real-time monitoring and ascending perilune in emergency. Finally the in-orbit utilization of safe shutdown strategy for Chang’e-3 lunar probe is introduced.

    braking at perilune; error; firing attitude; engine shutdown

    2014-12-05;

    2015-05-15

    國家重大科技專項工程

    董捷,男,工程師,從事深空探測器總體設(shè)計工作。Email:donghn13@163.com。

    V476.3

    A

    10.3969/j.issn.1673-8748.2015.04.003

    猜你喜歡
    發(fā)動機
    元征X-431實測:奔馳發(fā)動機編程
    2015款寶馬525Li行駛中發(fā)動機熄火
    2012年奔馳S600發(fā)動機故障燈偶爾點亮
    發(fā)動機空中起動包線擴(kuò)展試飛組織與實施
    RD-180超級火箭發(fā)動機的興衰
    太空探索(2016年8期)2016-07-10 09:21:58
    奔馳E200車發(fā)動機故障燈常亮
    奔馳E260冷車時發(fā)動機抖動
    新一代MTU2000發(fā)動機系列
    2013年車用發(fā)動機排放控制回顧(下)
    2013年車用發(fā)動機排放控制回顧(上)
    国产乱人伦免费视频| 精品99又大又爽又粗少妇毛片 | 亚洲无线观看免费| 国产精品乱码一区二三区的特点| 最新在线观看一区二区三区| 亚洲精品色激情综合| 18禁黄网站禁片免费观看直播| 搡老熟女国产l中国老女人| 免费人成在线观看视频色| xxxwww97欧美| 非洲黑人性xxxx精品又粗又长| 久久久色成人| 久久人人爽人人爽人人片va | 亚洲精品久久国产高清桃花| 香蕉av资源在线| 亚洲国产精品999在线| 国产主播在线观看一区二区| 一进一出抽搐gif免费好疼| 日本五十路高清| 赤兔流量卡办理| 熟妇人妻久久中文字幕3abv| 国产精品1区2区在线观看.| 亚洲 欧美 日韩 在线 免费| 精品久久久久久久人妻蜜臀av| 757午夜福利合集在线观看| 亚洲三级黄色毛片| 一级av片app| 国产精品伦人一区二区| 简卡轻食公司| 精品人妻熟女av久视频| 嫁个100分男人电影在线观看| 国产成人av教育| 国产三级黄色录像| 在线观看av片永久免费下载| 精品人妻1区二区| 国产精品久久久久久亚洲av鲁大| 女同久久另类99精品国产91| 欧美激情久久久久久爽电影| 757午夜福利合集在线观看| 天堂网av新在线| 成人一区二区视频在线观看| 亚洲欧美精品综合久久99| 日韩欧美国产在线观看| 1024手机看黄色片| 国产精品嫩草影院av在线观看 | 久久久精品欧美日韩精品| 久久久精品欧美日韩精品| 国产精品av视频在线免费观看| 亚洲片人在线观看| 成人av在线播放网站| 俺也久久电影网| 18禁黄网站禁片免费观看直播| 每晚都被弄得嗷嗷叫到高潮| 国产精品,欧美在线| 黄色视频,在线免费观看| 欧美激情在线99| 99在线视频只有这里精品首页| 在线播放国产精品三级| 成年版毛片免费区| 搡女人真爽免费视频火全软件 | 免费av毛片视频| 日韩欧美免费精品| 国产精品伦人一区二区| 色5月婷婷丁香| 欧美+亚洲+日韩+国产| 免费高清视频大片| 精品欧美国产一区二区三| 亚洲美女搞黄在线观看 | 国产av麻豆久久久久久久| 国产老妇女一区| 午夜福利18| 国产亚洲欧美在线一区二区| 中文字幕精品亚洲无线码一区| 国产精品一区二区免费欧美| 日本黄大片高清| 国产精品免费一区二区三区在线| 搡老妇女老女人老熟妇| 老司机福利观看| 内射极品少妇av片p| 午夜福利在线在线| 国产精品女同一区二区软件 | 小蜜桃在线观看免费完整版高清| 国产精品久久久久久亚洲av鲁大| 极品教师在线免费播放| 国产精品伦人一区二区| 婷婷精品国产亚洲av| 色噜噜av男人的天堂激情| 噜噜噜噜噜久久久久久91| 一进一出抽搐动态| netflix在线观看网站| 国产乱人视频| 久久久久久大精品| 欧美3d第一页| 激情在线观看视频在线高清| 色综合亚洲欧美另类图片| 夜夜爽天天搞| 观看免费一级毛片| 一进一出抽搐gif免费好疼| 国产高清视频在线观看网站| 熟女人妻精品中文字幕| netflix在线观看网站| av中文乱码字幕在线| 淫秽高清视频在线观看| 又爽又黄无遮挡网站| 国产伦一二天堂av在线观看| 3wmmmm亚洲av在线观看| 18禁裸乳无遮挡免费网站照片| 久久精品国产亚洲av天美| 日本成人三级电影网站| 长腿黑丝高跟| 老司机午夜十八禁免费视频| 老司机深夜福利视频在线观看| 国产三级在线视频| 最近中文字幕高清免费大全6 | 中文在线观看免费www的网站| 变态另类成人亚洲欧美熟女| 精品午夜福利在线看| 亚洲最大成人手机在线| 波多野结衣巨乳人妻| 好男人电影高清在线观看| АⅤ资源中文在线天堂| 亚洲欧美激情综合另类| 欧美精品啪啪一区二区三区| 国产91精品成人一区二区三区| 日韩高清综合在线| 成人一区二区视频在线观看| av中文乱码字幕在线| 天美传媒精品一区二区| 亚洲第一电影网av| 悠悠久久av| 女生性感内裤真人,穿戴方法视频| 精品熟女少妇八av免费久了| 国产伦在线观看视频一区| 亚洲欧美日韩无卡精品| 中国美女看黄片| 国产又黄又爽又无遮挡在线| 成人国产一区最新在线观看| 日本a在线网址| 亚洲av电影不卡..在线观看| 一区二区三区激情视频| 哪里可以看免费的av片| 欧美黄色淫秽网站| 变态另类成人亚洲欧美熟女| 国产亚洲欧美在线一区二区| 欧美另类亚洲清纯唯美| 成年女人看的毛片在线观看| a级毛片a级免费在线| aaaaa片日本免费| 国产伦一二天堂av在线观看| 91麻豆精品激情在线观看国产| 午夜福利视频1000在线观看| 国产亚洲精品综合一区在线观看| 男女下面进入的视频免费午夜| 亚洲精品影视一区二区三区av| 精品久久久久久久久久久久久| 欧美三级亚洲精品| 乱人视频在线观看| 宅男免费午夜| 非洲黑人性xxxx精品又粗又长| 久久久久久久久久成人| 亚洲av.av天堂| 毛片女人毛片| 91在线观看av| 精品一区二区免费观看| 非洲黑人性xxxx精品又粗又长| 国产探花极品一区二区| 免费电影在线观看免费观看| 欧美zozozo另类| 在线播放国产精品三级| 久久久久亚洲av毛片大全| 日韩欧美免费精品| 成人美女网站在线观看视频| 色在线成人网| 一本久久中文字幕| 欧美激情国产日韩精品一区| 深夜精品福利| 成人鲁丝片一二三区免费| 国产精品国产高清国产av| 久久99热这里只有精品18| 毛片女人毛片| 人妻制服诱惑在线中文字幕| 天堂动漫精品| 久久草成人影院| 波野结衣二区三区在线| 琪琪午夜伦伦电影理论片6080| 欧美色欧美亚洲另类二区| 国产精品乱码一区二三区的特点| 男女床上黄色一级片免费看| 又爽又黄无遮挡网站| 变态另类成人亚洲欧美熟女| 听说在线观看完整版免费高清| 亚洲欧美精品综合久久99| 在线观看66精品国产| 久久人人爽人人爽人人片va | 麻豆久久精品国产亚洲av| 中文字幕熟女人妻在线| 国产高清三级在线| eeuss影院久久| 久久久久九九精品影院| 99国产精品一区二区三区| 黄色女人牲交| 在线看三级毛片| 网址你懂的国产日韩在线| 欧美黑人欧美精品刺激| 国产精品伦人一区二区| 成人毛片a级毛片在线播放| 18禁在线播放成人免费| 男女那种视频在线观看| 亚洲欧美日韩无卡精品| 日韩有码中文字幕| 午夜两性在线视频| 91在线精品国自产拍蜜月| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 极品教师在线免费播放| 天堂网av新在线| av在线观看视频网站免费| 午夜福利在线观看吧| 欧美潮喷喷水| 国产私拍福利视频在线观看| 成人鲁丝片一二三区免费| 熟女电影av网| 国产高清有码在线观看视频| 日日摸夜夜添夜夜添小说| 亚洲国产精品sss在线观看| 男人的好看免费观看在线视频| a级一级毛片免费在线观看| 97人妻精品一区二区三区麻豆| 日韩高清综合在线| 国内精品一区二区在线观看| 国产三级黄色录像| 精品人妻一区二区三区麻豆 | 赤兔流量卡办理| 在线观看一区二区三区| 三级毛片av免费| 99久久九九国产精品国产免费| 国产精品一区二区三区四区久久| 成熟少妇高潮喷水视频| 日本五十路高清| 久久天躁狠狠躁夜夜2o2o| 欧美潮喷喷水| 精品乱码久久久久久99久播| 91麻豆av在线| 日韩精品青青久久久久久| 色播亚洲综合网| 真人做人爱边吃奶动态| 免费电影在线观看免费观看| 一卡2卡三卡四卡精品乱码亚洲| 日本黄大片高清| 欧美zozozo另类| 日日摸夜夜添夜夜添小说| av欧美777| 少妇人妻精品综合一区二区 | 男女之事视频高清在线观看| 国产av麻豆久久久久久久| 国产精品av视频在线免费观看| 熟女人妻精品中文字幕| 999久久久精品免费观看国产| 免费大片18禁| 可以在线观看毛片的网站| eeuss影院久久| 日韩 亚洲 欧美在线| 波野结衣二区三区在线| 69人妻影院| 精品国内亚洲2022精品成人| 日韩欧美一区二区三区在线观看| 精品日产1卡2卡| 99久久久亚洲精品蜜臀av| av视频在线观看入口| 成人欧美大片| 丰满人妻一区二区三区视频av| 国产亚洲欧美98| 999久久久精品免费观看国产| 女人十人毛片免费观看3o分钟| 欧美最黄视频在线播放免费| 久久久久久久久大av| 日本撒尿小便嘘嘘汇集6| 一个人看的www免费观看视频| 热99在线观看视频| 精品熟女少妇八av免费久了| 成年女人看的毛片在线观看| 成人无遮挡网站| 国产精品人妻久久久久久| 国内少妇人妻偷人精品xxx网站| 国产一区二区三区在线臀色熟女| 亚洲人成网站高清观看| 亚洲专区中文字幕在线| 老鸭窝网址在线观看| 91九色精品人成在线观看| 97超视频在线观看视频| 免费黄网站久久成人精品 | 欧美另类亚洲清纯唯美| 久久国产乱子免费精品| 在线免费观看不下载黄p国产 | 高清日韩中文字幕在线| 亚洲男人的天堂狠狠| 欧美又色又爽又黄视频| 午夜精品在线福利| 亚洲第一区二区三区不卡| 69人妻影院| netflix在线观看网站| 国产国拍精品亚洲av在线观看| 国产久久久一区二区三区| 最好的美女福利视频网| 国产成人av教育| 一区二区三区激情视频| 亚洲欧美清纯卡通| 在线a可以看的网站| 超碰av人人做人人爽久久| 精品国内亚洲2022精品成人| www.999成人在线观看| 欧美高清性xxxxhd video| 国产男靠女视频免费网站| 3wmmmm亚洲av在线观看| 蜜桃亚洲精品一区二区三区| 美女大奶头视频| 三级毛片av免费| 人妻夜夜爽99麻豆av| 国模一区二区三区四区视频| 精品福利观看| 亚洲五月天丁香| 内射极品少妇av片p| a级毛片a级免费在线| 免费在线观看日本一区| 亚洲熟妇熟女久久| 窝窝影院91人妻| 亚洲自偷自拍三级| 夜夜躁狠狠躁天天躁| 久久99热这里只有精品18| 午夜福利成人在线免费观看| 一级a爱片免费观看的视频| 美女 人体艺术 gogo| 永久网站在线| 久久性视频一级片| 999久久久精品免费观看国产| 午夜老司机福利剧场| 亚洲va日本ⅴa欧美va伊人久久| 哪里可以看免费的av片| 丁香六月欧美| 欧美日韩中文字幕国产精品一区二区三区| 啦啦啦韩国在线观看视频| 色吧在线观看| 国产大屁股一区二区在线视频| 国产主播在线观看一区二区| 别揉我奶头 嗯啊视频| 亚洲成人中文字幕在线播放| 女人被狂操c到高潮| 国产男靠女视频免费网站| 国产中年淑女户外野战色| 欧美一区二区精品小视频在线| 亚洲精品久久国产高清桃花| 成人高潮视频无遮挡免费网站| 长腿黑丝高跟| 国产成人欧美在线观看| 精品国产亚洲在线| 自拍偷自拍亚洲精品老妇| 国产主播在线观看一区二区| 麻豆国产av国片精品| 亚洲国产精品合色在线| 一级黄色大片毛片| 亚洲18禁久久av| 成年免费大片在线观看| 国产毛片a区久久久久| 国产av不卡久久| 亚洲久久久久久中文字幕| 成人国产综合亚洲| 琪琪午夜伦伦电影理论片6080| 亚洲天堂国产精品一区在线| 一区福利在线观看| 国产日本99.免费观看| 我的老师免费观看完整版| 69人妻影院| 窝窝影院91人妻| 丰满的人妻完整版| 桃色一区二区三区在线观看| 亚洲成人免费电影在线观看| 嫩草影院精品99| 国产aⅴ精品一区二区三区波| 免费看a级黄色片| 久久九九热精品免费| 日韩精品中文字幕看吧| 亚洲av日韩精品久久久久久密| 老女人水多毛片| 亚洲国产精品999在线| 欧美成狂野欧美在线观看| 国产白丝娇喘喷水9色精品| 亚洲av电影在线进入| 久久人妻av系列| 午夜a级毛片| 欧美3d第一页| 大型黄色视频在线免费观看| 久久久久国内视频| 亚洲国产精品合色在线| 嫩草影院精品99| 99精品在免费线老司机午夜| 久久这里只有精品中国| 亚洲欧美激情综合另类| 人妻夜夜爽99麻豆av| 中亚洲国语对白在线视频| 在线播放无遮挡| 亚洲熟妇熟女久久| 高清在线国产一区| 动漫黄色视频在线观看| 国产乱人视频| 欧美zozozo另类| 国产高清三级在线| 午夜福利18| 午夜精品久久久久久毛片777| 成人高潮视频无遮挡免费网站| 国产成人福利小说| 一本精品99久久精品77| av中文乱码字幕在线| 久久久国产成人免费| 国产精品一区二区性色av| 18禁裸乳无遮挡免费网站照片| 欧美精品啪啪一区二区三区| 超碰av人人做人人爽久久| 十八禁国产超污无遮挡网站| av在线观看视频网站免费| 亚洲无线在线观看| 久久99热6这里只有精品| 国产毛片a区久久久久| 99国产综合亚洲精品| 成熟少妇高潮喷水视频| 97碰自拍视频| 老鸭窝网址在线观看| 亚洲一区二区三区不卡视频| 91麻豆精品激情在线观看国产| 色综合欧美亚洲国产小说| 欧美一区二区精品小视频在线| 午夜福利在线观看吧| 欧美+日韩+精品| 国产精品三级大全| 国产aⅴ精品一区二区三区波| avwww免费| 久久精品国产99精品国产亚洲性色| 91久久精品电影网| 久久久久九九精品影院| 无人区码免费观看不卡| 午夜精品一区二区三区免费看| 国产男靠女视频免费网站| 最近在线观看免费完整版| 成人国产综合亚洲| 乱人视频在线观看| 国产欧美日韩精品一区二区| 日韩欧美一区二区三区在线观看| 91麻豆av在线| 最好的美女福利视频网| 能在线免费观看的黄片| 精品一区二区三区视频在线观看免费| 老鸭窝网址在线观看| 亚洲va日本ⅴa欧美va伊人久久| 国产一区二区三区在线臀色熟女| 免费观看精品视频网站| 国产蜜桃级精品一区二区三区| 国产成人欧美在线观看| 欧美又色又爽又黄视频| 别揉我奶头~嗯~啊~动态视频| 一级黄色大片毛片| 亚洲aⅴ乱码一区二区在线播放| 91狼人影院| 国产精品日韩av在线免费观看| 国产精品一区二区免费欧美| 免费电影在线观看免费观看| 精品久久久久久久久久久久久| 三级毛片av免费| 国产亚洲精品久久久久久毛片| av福利片在线观看| 舔av片在线| 在线播放无遮挡| 精品乱码久久久久久99久播| 十八禁国产超污无遮挡网站| 国产精品久久久久久精品电影| 亚洲国产精品久久男人天堂| 日本 av在线| 亚洲avbb在线观看| 好男人电影高清在线观看| 在线看三级毛片| 无遮挡黄片免费观看| 欧美3d第一页| 97人妻精品一区二区三区麻豆| 国产精品久久久久久久久免 | 他把我摸到了高潮在线观看| 国产精品久久电影中文字幕| 少妇被粗大猛烈的视频| 女人被狂操c到高潮| 亚洲乱码一区二区免费版| 毛片女人毛片| 亚洲第一区二区三区不卡| 桃红色精品国产亚洲av| 熟女电影av网| 俺也久久电影网| 亚洲中文字幕日韩| 午夜免费激情av| 变态另类成人亚洲欧美熟女| 欧美成人a在线观看| 日韩精品中文字幕看吧| 成人午夜高清在线视频| 又黄又爽又免费观看的视频| 成人午夜高清在线视频| 久久久久国产精品人妻aⅴ院| 午夜亚洲福利在线播放| 日日摸夜夜添夜夜添av毛片 | 精品久久久久久久末码| 99久久九九国产精品国产免费| 99国产精品一区二区蜜桃av| 婷婷亚洲欧美| 欧美成人a在线观看| 国产精品1区2区在线观看.| 色精品久久人妻99蜜桃| av福利片在线观看| 桃红色精品国产亚洲av| 精品一区二区三区视频在线观看免费| 成人av在线播放网站| 成人国产一区最新在线观看| 国产日本99.免费观看| 午夜福利高清视频| 99久久99久久久精品蜜桃| 午夜福利成人在线免费观看| 欧美日韩黄片免| 嫩草影院入口| 欧美日韩国产亚洲二区| 亚洲熟妇中文字幕五十中出| 亚洲人成网站高清观看| 深夜精品福利| 久久人妻av系列| 国产精品一及| 在线播放国产精品三级| 97碰自拍视频| 国产精品av视频在线免费观看| 狂野欧美白嫩少妇大欣赏| 国产乱人视频| 久久午夜亚洲精品久久| 日韩精品中文字幕看吧| 欧美性感艳星| 看黄色毛片网站| 国产私拍福利视频在线观看| 国产成人欧美在线观看| 亚洲专区中文字幕在线| 日韩 亚洲 欧美在线| 欧美日韩黄片免| 国产精品嫩草影院av在线观看 | 99国产精品一区二区蜜桃av| 欧美+日韩+精品| 日本a在线网址| 亚洲人成网站高清观看| 久久久久久久久久成人| 国内毛片毛片毛片毛片毛片| 欧美乱色亚洲激情| 国产精品爽爽va在线观看网站| 欧美丝袜亚洲另类 | 人人妻人人澡欧美一区二区| 又爽又黄a免费视频| 国产亚洲精品久久久久久毛片| 国产黄a三级三级三级人| 免费看日本二区| 成人毛片a级毛片在线播放| 久久人人爽人人爽人人片va | 亚洲国产精品成人综合色| 亚洲天堂国产精品一区在线| 99久久久亚洲精品蜜臀av| 少妇裸体淫交视频免费看高清| 久久人妻av系列| 亚洲精品亚洲一区二区| 亚洲精华国产精华精| 看片在线看免费视频| 国产私拍福利视频在线观看| 美女 人体艺术 gogo| 欧美精品啪啪一区二区三区| 久久人人爽人人爽人人片va | 女同久久另类99精品国产91| 欧美黑人欧美精品刺激| 免费在线观看成人毛片| 亚洲av五月六月丁香网| 国产三级黄色录像| 欧洲精品卡2卡3卡4卡5卡区| 精品一区二区三区人妻视频| 色播亚洲综合网| 免费一级毛片在线播放高清视频| 搡女人真爽免费视频火全软件 | 精品日产1卡2卡| 久久久久亚洲av毛片大全| 精品熟女少妇八av免费久了| 人妻夜夜爽99麻豆av| 成人av在线播放网站| 日韩免费av在线播放| 很黄的视频免费| 久久久久亚洲av毛片大全| 夜夜夜夜夜久久久久| 久久国产乱子免费精品| 亚洲久久久久久中文字幕| 亚洲美女黄片视频| 国产69精品久久久久777片| 2021天堂中文幕一二区在线观| 如何舔出高潮| 女同久久另类99精品国产91| 国内揄拍国产精品人妻在线| 久久精品国产亚洲av天美| 久久婷婷人人爽人人干人人爱| 欧美黄色淫秽网站| 一a级毛片在线观看| 欧美另类亚洲清纯唯美| 国产主播在线观看一区二区| 级片在线观看| 五月玫瑰六月丁香| 国产精品不卡视频一区二区 | 别揉我奶头 嗯啊视频| 悠悠久久av| 国内毛片毛片毛片毛片毛片| 午夜福利视频1000在线观看| 变态另类成人亚洲欧美熟女| 99久久无色码亚洲精品果冻| 三级毛片av免费| 欧美日本亚洲视频在线播放| 亚洲人成网站在线播| 精品一区二区免费观看|