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

    航天器再入隕落解體模型及分析預(yù)報(bào)策略研究

    2020-10-31 06:47:20唐小偉李四新石衛(wèi)波黨雷寧李志輝
    載人航天 2020年5期
    關(guān)鍵詞:隕落解體氣動(dòng)力

    唐小偉,李四新,石衛(wèi)波,黨雷寧,李志輝

    (中國空氣動(dòng)力研究與發(fā)展中心超高速空氣動(dòng)力研究所, 綿陽621000)

    1 引言

    隨著人類在大氣層外航天活動(dòng)的逐漸增多,航天器離軌隕落再入問題逐漸受到關(guān)注和重視。這類問題屬于典型的非常規(guī)再入問題[1-2]。 所謂非常規(guī)再入問題即指各種人造飛行物如失效衛(wèi)星、軌道艙、空間站、重復(fù)使用航天器(故障情況)、多級(jí)運(yùn)載火箭等在無控飛行或受控變軌后軌道衰降并再入地球大氣層的過程。 這些再入隕落的航天器往往不是用于在大氣層中飛行的,或者處于非設(shè)計(jì)飛行狀態(tài),沒有專門的氣動(dòng)布局外形及熱防護(hù)措施。 航天器超高速飛行進(jìn)入大氣層后,在強(qiáng)烈氣動(dòng)力/熱作用下,其原始構(gòu)型整器會(huì)出現(xiàn)解體,金屬材料軟化熔融,復(fù)合材料熱解/燒蝕以及期間或伴隨燃燒或爆炸等劇烈反應(yīng)現(xiàn)象。對(duì)航天器再入隕落解體過程的把握和定性理解是進(jìn)行建模及分析預(yù)報(bào)的前提。

    國外關(guān)于航天器再入隕落解體分析預(yù)測(cè)及地面風(fēng)險(xiǎn)評(píng)估的研究已有20 余年,建立了相對(duì)成熟的軟件系統(tǒng)[3-5],比如美國NASA 的DAS 和ORSAT[6-8]等,歐空局(ESA)的DRAMMA、SESAM 和SCARAB 等[9-11]。 國內(nèi)開發(fā)的有隕落預(yù)測(cè)軟件如DRAPS,受眾小、應(yīng)用力度不夠[12-14]。

    面向物體法和面向航天器法是由德國學(xué)者Lips 和Fritsche 提出的關(guān)于航天器再入隕落分析的分類方法[3]。 面向物體法的基本思路是將復(fù)雜航天器結(jié)構(gòu)簡(jiǎn)化為簡(jiǎn)單形狀物體,比如圓球、圓柱、箱體等,建模簡(jiǎn)單并且計(jì)算速度快,適宜于大量碎片再入的快速預(yù)測(cè),DAS 和ORSAT 均屬于這一類方法。 面向航天器法則盡可能模擬真實(shí)的航天器外形,采用基于表面網(wǎng)格的飛行器和碎片模型,理論上講具有更高的預(yù)測(cè)精度,但建模較復(fù)雜并需要較長的計(jì)算時(shí)間,SCARAB 是目前唯一屬于面向航天器法的軟件。

    唐小偉等[15-16]開展了再入隕落相關(guān)工程實(shí)例的應(yīng)用分析,包括對(duì)多級(jí)運(yùn)載火箭發(fā)動(dòng)機(jī)殘骸高速墜入大氣層、失效小衛(wèi)星隕落再入損毀情況及某衛(wèi)星使用的放射性同位素?zé)嵩唇M件再入進(jìn)行了計(jì)算分析。 基于彈道-氣動(dòng)力-氣動(dòng)熱的綜合計(jì)算分析的技術(shù)途徑,針對(duì)LTG(類似天舟、天宮等)軌道艙再入隕落解體過程進(jìn)行了分析預(yù)報(bào)。經(jīng)驗(yàn)表明,航天器再入隕落解體過程中各種現(xiàn)象具有很大隨機(jī)性,從解決工程問題的角度,復(fù)雜的航天器精細(xì)建模不一定能夠反映系統(tǒng)的不確定性。 航天器再入隕落解體過程的分析預(yù)報(bào),從技術(shù)途徑概述就是要得出氣動(dòng)力/熱和航天器物理(解體/運(yùn)動(dòng))化學(xué)(熱解/燒蝕)狀態(tài)交互作用的過程。 從技術(shù)邏輯分析,航天器解體及其導(dǎo)致的物形變化是氣動(dòng)力/熱分析的基礎(chǔ);對(duì)航天器及其解體后對(duì)象氣動(dòng)力/熱的準(zhǔn)確分析,是研究對(duì)象飛行運(yùn)動(dòng)、殘骸存活性、殘骸落區(qū)及地面風(fēng)險(xiǎn)評(píng)估的關(guān)鍵。

    為此,本文首先針對(duì)航天器再入隕落解體分析預(yù)報(bào)的研究對(duì)象,綜合以往工作經(jīng)驗(yàn)歸納提煉出一種新型簡(jiǎn)潔的航天器再入隕落解體模型——三層級(jí)模型;然后針對(duì)航天器再入隕落解體分析預(yù)報(bào)的技術(shù)途徑,提出基于條件邊界的參數(shù)統(tǒng)計(jì)方法;最后通過相關(guān)軟件研制和一例大型航天器再入隕落解體過程分析預(yù)報(bào)表明文章提出的模型和方法的適用性。

    2 再入隕落解體過程技術(shù)簡(jiǎn)析

    如前所述,航天器再入隕落解體過程分析是典型的非常規(guī)再入問題[1-2],這類航天器再入過程中的主要技術(shù)問題包括飛行運(yùn)動(dòng)、氣動(dòng)力、氣動(dòng)熱和結(jié)構(gòu)解體4 個(gè)方面。

    對(duì)航天器再入隕落過程進(jìn)行分析時(shí),對(duì)飛行運(yùn)動(dòng)的把握是貫穿全過程的主要線索,也是分析預(yù)報(bào)建模關(guān)注的主要內(nèi)容之一。 在航天器整器解體前,對(duì)其進(jìn)行精確的氣動(dòng)力/熱及運(yùn)動(dòng)分析具有重大意義,有助于對(duì)整器解體時(shí)刻的飛行高度、速度等參數(shù)進(jìn)行準(zhǔn)確推斷;由于隕落解體后的航天器部件或碎片外形一般都比較復(fù)雜,質(zhì)量特性及外形仍然處于持續(xù)變化之中,導(dǎo)致飛行運(yùn)動(dòng)姿態(tài)及相關(guān)的氣動(dòng)力/熱作用存在極大不確定性。 對(duì)于此類問題,要獲得隕落體部件或碎片的姿態(tài)變化歷程是困難的,且不具備充分的工程意義。

    氣動(dòng)力問題是隕落體再入/進(jìn)入大氣層時(shí)面臨的最重要和最復(fù)雜的技術(shù)問題之一。 力是運(yùn)動(dòng)變化的直接原因;隕落體受到的重力相對(duì)較為確定,因而氣動(dòng)力的準(zhǔn)確分析是隕落體運(yùn)動(dòng)動(dòng)力學(xué)評(píng)估的基礎(chǔ)。 隕落體再入飛行過程中,穿越了自由分子流區(qū)域(高度100~120 km 以上高空)及過渡流區(qū)域(約高度70~120 km 區(qū)間),其部分殘骸碎片可能會(huì)穿過連續(xù)流區(qū)域(約高度70 km 以下)直至地球表面。 隕落體飛行過程中涉及到氣動(dòng)力分析,以及跨流域、復(fù)雜外形、氣動(dòng)物理等方面的多種手段。 對(duì)這些氣動(dòng)力相關(guān)科學(xué)問題的探索,可牽引發(fā)展以求解Boltzmann 方程可計(jì)算建模為核心的跨流域氣體動(dòng)理論統(tǒng)一算法基礎(chǔ)研究、稀薄氣體動(dòng)力學(xué)DSMC 數(shù)值模擬研究、過渡流區(qū)域N-S/DSMC 耦合算法研究及考慮熱化學(xué)非平衡效應(yīng)的近連續(xù)流區(qū)、連續(xù)流區(qū)N-S 方程數(shù)值模擬研究。

    隕落體飛行穿越大氣層過程中,高超聲速運(yùn)動(dòng)必然會(huì)帶來強(qiáng)烈的氣動(dòng)熱作用。 氣動(dòng)熱分析主要反饋其對(duì)隕落體的外形影響及物性變化,作為解體分析判斷的依據(jù),并作為氣動(dòng)力分析的重要基礎(chǔ)。 氣動(dòng)熱問題可細(xì)分為熱環(huán)境、金屬材料軟化/熔融和復(fù)合材料熱解/燒蝕3 個(gè)方面,劇烈的爆炸或燃燒反應(yīng)也和氣動(dòng)熱作用密切相關(guān)。

    熱環(huán)境分析的主要目標(biāo)是獲得隕落體表面的對(duì)流換熱情況,這是后續(xù)材料結(jié)構(gòu)破壞及解體分析的基礎(chǔ)。 隕落體主要是一些壽命末期的大型軌道飛行器、人造衛(wèi)星或上面級(jí)運(yùn)載火箭分離拋棄的殘骸等,它們絕大部分由金屬或合金材料(統(tǒng)稱為金屬材料)構(gòu)成。 因此,對(duì)金屬材料軟化/熔融導(dǎo)致的破壞情況進(jìn)行分析評(píng)估是非常重要的方面。 在高速氣流氣動(dòng)加熱作用下,金屬材料構(gòu)件會(huì)逐漸升溫;當(dāng)溫度上升到一定程度,金屬材料會(huì)出現(xiàn)軟化現(xiàn)象;當(dāng)溫度達(dá)到熔點(diǎn),金屬材料將發(fā)生熔融。 航天器再入隕落飛行過程中,對(duì)金屬材料軟化/熔融的分析目的在于獲得其由于軟化或熔融而導(dǎo)致的結(jié)構(gòu)破壞和解體情況。 雖然再入大氣層內(nèi)的人工天體類隕落體一般絕大部分由金屬材料構(gòu)成,但是仍然有部分部件是由復(fù)合材料構(gòu)成的,主要是一些氣瓶、貯箱之類。 復(fù)合材料在氣動(dòng)熱作用下一般要經(jīng)歷升溫、熱解和燒蝕等復(fù)雜的物理化學(xué)過程,這是航天器再入隕落分析預(yù)報(bào)中需要重點(diǎn)關(guān)注的技術(shù)問題。

    進(jìn)入大氣層內(nèi)的隕落體可能呈現(xiàn)出的最明顯表象即是結(jié)構(gòu)失效崩潰及激烈的解體。 結(jié)構(gòu)解體形成的部件或碎片形狀是氣動(dòng)力、氣動(dòng)熱評(píng)估分析的幾何基礎(chǔ),解體后部件或碎片的質(zhì)量特性及運(yùn)動(dòng)參數(shù)同時(shí)也是后續(xù)彈道預(yù)測(cè)的主要依據(jù)。 如果要剖析結(jié)構(gòu)解體的細(xì)節(jié),須基于固體力學(xué)為基礎(chǔ)的嚴(yán)格的內(nèi)應(yīng)力分析,融合氣動(dòng)力/熱對(duì)結(jié)構(gòu)產(chǎn)生的物理化學(xué)作用,并考慮其它影響因素(如重力、自旋等)的作用。 原則上,結(jié)構(gòu)解體屬于典型的固體力學(xué)問題,屬于固體力學(xué)中材料的應(yīng)變應(yīng)力分析技術(shù)問題。 通過數(shù)值模擬方法可獲得每一處結(jié)構(gòu)微元體的應(yīng)變、應(yīng)力,結(jié)合材料當(dāng)?shù)貤l件下的物性參數(shù)和破壞判據(jù)給出定量仿真結(jié)果。 航天器再入隕落飛行過程中,隕落體的力學(xué)環(huán)境分析是非常困難的技術(shù)難點(diǎn),而且力學(xué)環(huán)境和氣動(dòng)加熱作用及材料物性密切相關(guān)。

    3 三層級(jí)模型

    為了對(duì)航天器再入隕落解體過程分析預(yù)報(bào)提供必要物形幾何前提,提出了一種新型簡(jiǎn)潔的航天器再入隕落解體模型——三層級(jí)模型(Three Phases Model,TPM)。 該模型的基本思想是面向物體法,同時(shí)融入了面向航天器法的思維。

    航天器再入隕落解體三層級(jí)模型的3 個(gè)層級(jí)是:①系統(tǒng)/子系統(tǒng)(Whole);②部件(Part);③碎片/微粒(Debris)。 因此,三層級(jí)模型又可稱為WPDM(Whole-Part-Debris Model),如圖1 所示。

    圖1 航天器再入隕落解體三層級(jí)模型Fig.1 Three phases model for spacecraft disintegration during reentry fall

    航天器再入隕落解體三層級(jí)模型的系統(tǒng)/子系統(tǒng),即指航天器整器形式。 對(duì)于中小型在軌運(yùn)行航天器,可以是包括太陽能電池帆板或不包括太陽能電池帆板的形式;對(duì)于特大型在軌運(yùn)行航天器,如國際空間站,顯然在某些隕落階段可以分為若干子系統(tǒng)。 通過實(shí)例分析可知,對(duì)航天器整器解體時(shí)刻的狀態(tài)參數(shù),尤其是高度、速度、彈道傾角、飛行方位角的精確程度,對(duì)后續(xù)存活殘骸的落區(qū)范圍估計(jì)起到?jīng)Q定性作用,對(duì)殘骸存活性分析和地面風(fēng)險(xiǎn)評(píng)估也有影響。 這是三層級(jí)模型中特別指定系統(tǒng)/子系統(tǒng)層級(jí)的最大意義。

    航天器再入隕落解體三層級(jí)模型的部件指的是航天器內(nèi)部具有一定功能性及結(jié)構(gòu)外形確定性的內(nèi)容物,包括各種功能部件、氣瓶及燃料貯箱等。 當(dāng)隕落至一定高度在氣動(dòng)力/熱作用下,航天器整器外殼及附屬連接固定結(jié)構(gòu)失效,其內(nèi)部各部件在氣流沖擊下將飛散開來。 整器解體時(shí)刻參數(shù)即是各部件承受嚴(yán)苛氣動(dòng)力/熱及運(yùn)動(dòng)分析的初始參數(shù)。 原則上,可以對(duì)具有確定幾何特性和質(zhì)量特性的部件進(jìn)行精細(xì)的氣動(dòng)力/熱和運(yùn)動(dòng)分析,但是在熔融或燒蝕過程中,任何部件的幾何特性和質(zhì)量特性的確定都是暫時(shí)性的。

    4 塊條片模型

    航天器再入隕落解體三層級(jí)模型的碎片/微粒指的是不可辨識(shí)原本幾何形態(tài)的殘留物。 碎片/微粒的來源有3 個(gè)方面:①航天器整器或部件熔融燒蝕后殘留物;②航天器整器或部件在氣流剪切撕扯下產(chǎn)生的碎片;③熔融物再凝結(jié)并得以存活的殘留物。 碎片/微粒的特征尺度下限根據(jù)關(guān)注的問題進(jìn)行界定。 碎片/微粒的技術(shù)分析不再考慮氣動(dòng)熱作用,即它們是被作為落地實(shí)物來看待的。 碎片/微粒的總量需封閉與航天器的質(zhì)量守恒及材料屬性,同時(shí)考慮扣除燃燒或霧化湮沒的部分(一般情況下是絕大部分)。

    碎片/微粒層級(jí)的外形建模過程中,對(duì)其幾何特征進(jìn)行高度抽象。 幾何建模抽象基于完全歸納法原則,即可以覆蓋不同的幾何特征,從而覆蓋體現(xiàn)相應(yīng)的氣動(dòng)特性及飛行運(yùn)動(dòng)特性。 碎片/微粒按幾何特征分為3 類:①塊(Block);②條(Column);③片(Slice)。 因此,碎片/微粒層級(jí)的這種劃分模式可稱為塊條片模型(Block-Column-Slice Model, BCSM)。

    3 類碎片/微粒的幾何形狀中(圖2),塊以球或橢球指代,由短長軸比定形;片以半圓邊緣圓片指代,由厚度直徑比定形;條以半球端面圓柱指代,由長細(xì)比定形。 3 類碎片/微粒的尺度范圍(上下限)根據(jù)需要設(shè)定,塊由特征尺度長軸確定;條由特征尺度柱長確定;片由特征尺度直徑確定。 碎片/微粒的幾何形狀采用表面圓潤的設(shè)定,基于隕落殘骸歷經(jīng)燒蝕及氣流沖刷后一般不再具有棱角分明的特征。

    圖2 碎片/微粒層級(jí)的幾何建模Fig.2 Geometric modeling at debris/particle phase

    碎片/微粒需滿足材料(密度)及對(duì)應(yīng)總質(zhì)量約束條件。 根據(jù)需要對(duì)3 類碎片/微粒的相關(guān)參數(shù)具體進(jìn)行設(shè)定或覆蓋。 這些參數(shù)的設(shè)定和控制依賴于基礎(chǔ)研究成果對(duì)隕落解體過程的認(rèn)識(shí)程度。

    某些情況下,若系統(tǒng)/子系統(tǒng)及整器發(fā)生劇烈解體過程(如爆破等),則部件級(jí)層級(jí)可能不存在,直接由系統(tǒng)/子系統(tǒng)層級(jí)演化至碎片/微粒層級(jí)。 在部件級(jí)存在情況下,某些碎片/微粒實(shí)體也可能由整器(主要是殼體)解體直接產(chǎn)生,這也是圖1 上部虛線的含義。

    表1 定義了碎片/微粒的塊條片模型相關(guān)的幾何及材料屬性參數(shù)。

    表1 碎片/微粒層級(jí)塊條片模型的幾何及材料屬性參數(shù)Table 1 Geometrical and material parameters of BCSM

    對(duì)于每一個(gè)具體的碎片/微粒,其幾何及材料屬性確定的個(gè)體數(shù)量可表示為式(1)。

    上述參量中,T 和M 為非連續(xù)性變量,R 和L為連續(xù)性變量,特別地,塊、條、片類型對(duì)應(yīng)的個(gè)體數(shù)量可分別表示為式(2)。

    連續(xù)性變量R 和L 一般可表示為在某種范圍內(nèi)的分布函數(shù)的形式實(shí)際應(yīng)用中根據(jù)需要進(jìn)行離散化處理。 基于隕落過程中熔融及燒蝕作用具有使碎片外形圓潤化的趨勢(shì),建議定形參數(shù)R 的取值下限為1,上限根據(jù)基礎(chǔ)研究和觀測(cè)結(jié)果的經(jīng)驗(yàn)積累進(jìn)行評(píng)估確定。 此外,3 種碎片類型的定形參數(shù)為1 時(shí),實(shí)際都?xì)w化為基準(zhǔn)的圓球形。 特征尺度L 必須具有明確的上下限,這是根據(jù)所關(guān)注的問題及實(shí)際研究對(duì)象的情況評(píng)估確定的。 具有相同材料屬性的任何類型碎片,應(yīng)該滿足同一材料的質(zhì)量約束條件。

    碎片數(shù)量的關(guān)鍵依賴參數(shù)是特征尺度(或碎片質(zhì)量),根據(jù)類似研究[17-19]表明,碎片數(shù)量和特征尺度近似于指數(shù)變化規(guī)律,即尺度越小的碎片,其數(shù)量越多。 為此,構(gòu)建特征尺度大于L 的碎片數(shù)量的表達(dá)式如式(3)所示。

    式中,N0為納入統(tǒng)計(jì)的碎片數(shù)量總數(shù), b 為解體強(qiáng)度系數(shù),L1和L2分別為碎片特征尺度的下限和上限,L 取值介于L1和L2之間。

    以上碎片/微粒層級(jí)的塊條片模型,根據(jù)實(shí)際問題的工程要求情況,針對(duì)類別、定形、特征尺度和材料等方面,可以進(jìn)行必要的簡(jiǎn)并。

    5 三層級(jí)模型的技術(shù)特點(diǎn)

    航天器再入隕落解體三層級(jí)模型能夠涵蓋航天器解體研究對(duì)象的結(jié)構(gòu)組成及碎片各種幾何特征,這些幾何外形特征能夠反映研究對(duì)象的氣動(dòng)力/熱及飛行運(yùn)動(dòng)特性,通過調(diào)整參數(shù)的應(yīng)用,針對(duì)航天器不同解體情況具有普遍適應(yīng)性。

    航天器再入隕落解體三層級(jí)模型的系統(tǒng)/子系統(tǒng)層級(jí)和部件層級(jí)具有面向航天器法的性質(zhì),碎片/微粒層級(jí)具有面向物體法的性質(zhì)。 參照國內(nèi)外現(xiàn)有方法模型,三層級(jí)模型具有如下特點(diǎn):

    1)融合了面向物體法和面向航天器法的思想,是該兩類方法的綜合體現(xiàn);

    2)三層級(jí)模型對(duì)研究對(duì)象結(jié)構(gòu)定義清晰明確,更精煉適用,便于工程操作;

    3)對(duì)系統(tǒng)/子系統(tǒng)的明確關(guān)注有利于精確分析航天器整器解體時(shí)刻的狀態(tài)參數(shù);

    4)對(duì)部件層級(jí)的定義有利于對(duì)隕落解體過程的細(xì)節(jié)情況進(jìn)行探索;

    5)碎片/微粒層級(jí)的塊條片模型能夠代表不同幾何特征及相應(yīng)的氣動(dòng)力/熱和飛行運(yùn)動(dòng)特性;

    6)三層級(jí)模型能夠滿足工程分析需要,也有助于挖掘科學(xué)問題推進(jìn)基礎(chǔ)研究。

    6 基于條件邊界的參數(shù)統(tǒng)計(jì)方法

    基于航天器再入隕落解體三層級(jí)模型,原則上對(duì)航天器再入隕落解體過程的分析預(yù)報(bào)提供了具體的研究對(duì)象,但是航天器再入隕落解體過程中的氣動(dòng)力/熱現(xiàn)象實(shí)際上是外形及環(huán)境不斷變化的非定常問題,由此導(dǎo)致隕落過程必然是一個(gè)隨機(jī)過程。 換句話說,假設(shè)初始狀態(tài)完全相同的2 個(gè)隕落體,其隕落解體過程、殘骸存活情況及實(shí)際落區(qū)等都會(huì)有所不同。 對(duì)隨機(jī)過程的分析不能局限于確定性問題的解決思路,為此提出一種基于條件邊界的參數(shù)統(tǒng)計(jì)方法,作為航天器再入隕落分析預(yù)報(bào)的基本策略。

    基于條件邊界的參數(shù)統(tǒng)計(jì)方法即是針對(duì)航天器再入隕落過程仿真建模分析中存在不確定性的氣動(dòng)力/熱及飛行狀態(tài)參數(shù),根據(jù)現(xiàn)有基礎(chǔ)研究成果評(píng)估確定其上下限即條件邊界;在隕落過程分析時(shí),在這些條件邊界值域內(nèi),運(yùn)用合適的參數(shù)分布及統(tǒng)計(jì)模型,對(duì)關(guān)注的目標(biāo)參數(shù)進(jìn)行定量分析。由此獲得統(tǒng)計(jì)意義上的航天器再入隕落過程及結(jié)果,這些結(jié)果應(yīng)該是具有一定置信度的分布帶形式。

    從飛行運(yùn)動(dòng)方面分析,由于質(zhì)量特性對(duì)于隕落體往往是變化不定,很難把握的,且在高度100 km以上大氣極其稀薄、動(dòng)壓非常小,航天器飛行過程的姿態(tài)轉(zhuǎn)動(dòng)運(yùn)動(dòng)情況對(duì)平動(dòng)運(yùn)動(dòng)幾乎不產(chǎn)生影響,一般情況下采用三自由度彈道方程即可。 但當(dāng)需要把握解體前運(yùn)動(dòng)細(xì)節(jié)時(shí)可采用六自由度彈道方程。 基于三自由度橫向氣動(dòng)力的需要,可把變形或局部破壞后的隕落體以及殘骸碎片的運(yùn)動(dòng)姿態(tài)參數(shù)納入合適的分布統(tǒng)計(jì)模型考慮。

    從氣動(dòng)力方面分析,在高超聲速條件下隕落體及其殘骸外形復(fù)雜且多變。 如果存在部分殘骸體接近地球表面且速度減至高超聲速以下,初步評(píng)估獲得的經(jīng)驗(yàn)表明,這類較低速度飛行的殘骸碎片不會(huì)在氣動(dòng)力/熱作用下發(fā)生破壞,在這種情況下其當(dāng)?shù)貜椀纼A角一般已接近-90°,即其落點(diǎn)經(jīng)過該類模式飛行后不會(huì)發(fā)生顯著變化。 因此,采用氣動(dòng)力數(shù)據(jù)表或快速工程方法是沿彈道計(jì)算的主要選擇,必要的配平姿態(tài)評(píng)估可通過數(shù)值模擬分析,任何先進(jìn)的大氣模型對(duì)大氣參數(shù)的推算也存在不確定度,故需要納入條件邊界參數(shù)統(tǒng)計(jì)的參量包括氣動(dòng)阻力系數(shù)及大氣密度。

    從氣動(dòng)熱方面分析,由于對(duì)變形隕落體及殘骸碎片形狀的不確定性,對(duì)氣動(dòng)分析的外形可以通過分類抽象模式化處理。 據(jù)此,對(duì)氣動(dòng)熱直接相關(guān)參數(shù)進(jìn)行參數(shù)分布統(tǒng)計(jì)已不存在邏輯上的必要性;實(shí)際上氣動(dòng)熱的影響可轉(zhuǎn)化歸于對(duì)氣動(dòng)力不確定度的影響。 對(duì)金屬材料部件的傳熱熔融計(jì)算及復(fù)合材料部件的熱解/燒蝕計(jì)算可采用零維或一維計(jì)算模型。

    結(jié)構(gòu)解體的分析是航天器再入隕落分析預(yù)報(bào)技術(shù)中最困難和最缺乏針對(duì)性基礎(chǔ)研究成果支撐的部分。 技術(shù)途徑是結(jié)合航天器再入隕落解體三層級(jí)模型,通過建立解體條件準(zhǔn)則進(jìn)行分析,即根據(jù)若干基礎(chǔ)研究的成果提煉相關(guān)的解體破壞判據(jù)。 可以考慮的解體條件準(zhǔn)則包括:①高度準(zhǔn)則,即設(shè)定在某一高度值發(fā)生某種形式的解體過程;②熔融準(zhǔn)則,即設(shè)定金屬材料在達(dá)到熔融溫度發(fā)生某種形式的解體過程;③熱解準(zhǔn)則,即設(shè)定復(fù)合材料在達(dá)到熱解溫度時(shí)發(fā)生某種形式的解體過程;④溫度準(zhǔn)則,即認(rèn)為物面溫度達(dá)到某一設(shè)定值將發(fā)生某種形式的解體過程;⑤熱流準(zhǔn)則,即認(rèn)為累計(jì)熱流或瞬時(shí)熱流達(dá)到某一設(shè)定值將發(fā)生某種形式的解體過程;⑥動(dòng)壓準(zhǔn)則,即認(rèn)為動(dòng)壓達(dá)到某一設(shè)定值將發(fā)生某種形式的解體過程;⑦時(shí)間準(zhǔn)則,即設(shè)定在某一時(shí)刻發(fā)生某種形式的解體過程;⑧綜合準(zhǔn)則,上述兩種或多種準(zhǔn)則的綜合,可包括權(quán)重和、任一準(zhǔn)則達(dá)成、全體準(zhǔn)則達(dá)成等。

    上述熔融準(zhǔn)則及熱解準(zhǔn)則可以認(rèn)為是溫度準(zhǔn)則和熱流準(zhǔn)則的某種復(fù)合形式,但由于金屬材料在目前的航天器結(jié)構(gòu)材料中占據(jù)最大百分比,故把熔融準(zhǔn)則單列便于應(yīng)用分析;復(fù)合材料則是目前航天器上有特殊要求的容器類的主要材料,因此也把熱解準(zhǔn)則單列。

    在DAS 和ORSAT 中,采用最簡(jiǎn)單的高度準(zhǔn)則預(yù)測(cè)解體;DRAPS 除了高度準(zhǔn)則外還提供溫度準(zhǔn)則、燒蝕準(zhǔn)則(燒蝕到一定程度發(fā)生解體)和綜合準(zhǔn)則(達(dá)到上述任何一個(gè)解體條件就認(rèn)為解體)。 本文擬建立結(jié)構(gòu)解體的較為系統(tǒng)完整的解體條件準(zhǔn)則及解體形式的框架體系。 這些準(zhǔn)則的選用根據(jù)實(shí)際的隕落過程特點(diǎn);準(zhǔn)則的量化建立依賴于實(shí)際研究的隕落體對(duì)象相關(guān)的基礎(chǔ)支撐成果的積累。

    原則上,解體形式可以包括:①部件級(jí)解體,整器被肢解為部件;②碎片級(jí)解體,整器或部件被解體為碎片或微粒;③變形,外形發(fā)生變化。

    由于外形變化在解體過程中是一種連續(xù)不斷的持續(xù)性過程,因此要對(duì)此進(jìn)行精細(xì)分析或仿真是非常困難的;從關(guān)注解體結(jié)果的角度,可以把變形解體形式的影響納入其它解體條件的值域方面進(jìn)行考慮。 結(jié)構(gòu)解體分析中需要納入條件邊界參數(shù)統(tǒng)計(jì)的最重要的參量包括:解體碎片的尺度分布、解體碎片的形狀模式分布、解體高度范圍等。

    隕落飛行過程中物體形狀及氣動(dòng)力/熱局部參數(shù)的隨機(jī)性使得采用面向隕落體的完全有限元力/熱及結(jié)構(gòu)應(yīng)力耦合數(shù)值模擬代價(jià)巨大,需要大規(guī)模高性能并行計(jì)算,且其個(gè)別條件下的確定性數(shù)值分析結(jié)果并不能代表解體過程參數(shù)的隨機(jī)性。 但是,必須通過基礎(chǔ)研究成果來提煉工程上需要的支撐數(shù)據(jù),并不斷深入研究逐步完善工程適用的解體判據(jù)。 如ORSAT 中默認(rèn)的解體高度為78 km,然而本文研究表明,解體高度和研究對(duì)象材質(zhì)及飛行模式關(guān)系密切,有必要進(jìn)行精細(xì)化分析并設(shè)定合適的量化解體準(zhǔn)則。

    初步歸納,航天器再入隕落解體分析預(yù)報(bào)的基于條件邊界的參數(shù)統(tǒng)計(jì)方法涉及的統(tǒng)計(jì)分析參量如下:①隕落分析飛行起始參數(shù),含初始位置三分量(坐標(biāo)值x/y/z 或經(jīng)度/緯度/高度)和初始速度三分量(速率/初始彈道傾角/初始方位角),共6 個(gè)參量;②阻力系數(shù);③大氣密度;④解體高度;⑤解體碎片尺度;⑥解體碎片形狀;⑦解體碎片姿態(tài)角。 上述12 個(gè)參數(shù)可根據(jù)情況進(jìn)一步簡(jiǎn)并,必要時(shí)也可擴(kuò)展。 為了使用基于條件邊界的參數(shù)統(tǒng)計(jì)方法,需把握上述參數(shù)的基準(zhǔn)值、不確定度及其分布形式,其中對(duì)結(jié)構(gòu)解體相關(guān)參數(shù)的設(shè)定是最關(guān)鍵和最困難的。

    7 實(shí)例應(yīng)用分析

    基于上述航天器再入隕落解體三層級(jí)模型及基于條件邊界的參數(shù)統(tǒng)計(jì)方法的分析策略及思路,開發(fā)研制了一套針對(duì)性的軟件系統(tǒng)。 軟件專門用于大型航天器再入隕落解體過程分析預(yù)報(bào)及地面風(fēng)險(xiǎn)評(píng)估,采用數(shù)據(jù)庫技術(shù)對(duì)軟件系統(tǒng)研究對(duì)象相關(guān)的素材資源、方法流程調(diào)度及結(jié)果數(shù)據(jù)進(jìn)行綜合管理,并集成了各類相關(guān)基礎(chǔ)研究算法模塊。 圖3 為軟件系統(tǒng)基本邏輯流程圖,圖4 為軟件系統(tǒng)運(yùn)行時(shí)的典型界面。

    圖3 航天器再入隕落解體分析預(yù)報(bào)軟件系統(tǒng)基本流程圖Fig.3 Basic flow chart of analysis and prediction software system for spacecraft disintegration during reentry fall

    圖4 航天器再入隕落解體分析預(yù)報(bào)軟件系統(tǒng)典型界面Fig.4 Typical interface of analysis and prediction software system for spacecraft disintegration during reentry fall

    基于上述模型方法及開發(fā)的軟件系統(tǒng),對(duì)一個(gè)大型航天器實(shí)例進(jìn)行了初步分析。 由于相關(guān)基礎(chǔ)研究支撐工作是一個(gè)長期持續(xù)的過程,對(duì)解體模型中一些參數(shù)的邊界及分布設(shè)定還需要不斷細(xì)化和充實(shí)完善。

    圖5(a)為用于隕落分析的大型航天器(簡(jiǎn)稱LTG),該航天器整器由一兩艙結(jié)構(gòu)及兩側(cè)太陽能電池帆板組成。 艙體材質(zhì)為鋁合金,部分組件外露艙體,包括中繼天線、軌控和姿控發(fā)動(dòng)機(jī)噴管、各種相機(jī)光學(xué)窗口等。 圖5(b)為無太陽帆板的兩艙組合體,其與帶太陽帆板的整器共同構(gòu)成算例應(yīng)用分析中LTG 隕落解體三層級(jí)模型的系統(tǒng)/子系統(tǒng)層級(jí)。

    圖5 LTG 隕落解體三層級(jí)模型之系統(tǒng)/子系統(tǒng)層級(jí)Fig.5 System/sub-system phase example of three phases model for a LTG spacecraft

    LTG 再入隕落解體三層級(jí)模型的部件級(jí)大致包括了圖6 中的十多種類型,它們是組成該航天器的核心功能組件。 大部分內(nèi)部部件由鋁合金材料構(gòu)成,少量不銹鋼部件;軌控和姿控發(fā)動(dòng)機(jī)噴管由鈮合金構(gòu)成;還有部分由復(fù)合材料或碳纖維構(gòu)成的瓶罐類容器;此外還有蓄電池模塊。

    圖6 LTG 隕落解體三層級(jí)模型之部件層級(jí)Fig.6 Part phase example of three phases model for a LTG spacecraft

    本例中,LTG 再入隕落解體三層級(jí)模型的碎片/微粒層級(jí),簡(jiǎn)化設(shè)定為球體、長方體、片類、桿類(長細(xì)比=10)和桿類(長細(xì)比=5),并在此基礎(chǔ)上采用分組方法進(jìn)行了統(tǒng)計(jì)計(jì)算分析。 圖7 為某設(shè)定條件下碎片落區(qū)的計(jì)算結(jié)果。

    圖7 LTG 隕落解體三層級(jí)模型之碎片/微粒層級(jí)計(jì)算結(jié)果Fig.7 Debris/particle phase example of three phases model for a LTG spacecraft

    針對(duì)LTG 再入隕落解體過程的系統(tǒng)建模和模擬,需要把握如下情況:隕落再入初始條件的確定;太陽電池帆板脫落的時(shí)機(jī)及對(duì)稱性;兩艙組合體逐漸熔融伴隨解體的復(fù)雜性;部件存活期間形變及氣動(dòng)力熱局部參數(shù)的隨機(jī)性;燃料貯箱及電池模塊物理或化學(xué)爆破的可能性及影響。 上述情況的隨機(jī)性使得對(duì)整器解體的隕落分析的條件參數(shù)設(shè)置原則是盡量采取偏于保守的設(shè)定,包括對(duì)條件狀態(tài)參數(shù)基準(zhǔn)值及其值域范圍,都需要在工程對(duì)象可能的范圍內(nèi)進(jìn)行綜合評(píng)估分析。

    計(jì)算分析涉及的結(jié)果數(shù)據(jù)眾多,本文給出部分結(jié)果。 圖8 為LTG 兩艙組合體于高度100 km、迎角170°時(shí)表面溫度分布,此時(shí)其艙體迎風(fēng)端面溫度已超過鋁合金材質(zhì)熔點(diǎn),熔融速率取決于當(dāng)?shù)責(zé)崃鞔笮 ?/p>

    圖8 LTG 兩艙組合體表面溫度分布(高度100 km、迎角170°)Fig.8 Surface temperature distribution of LTG two cabin combination (hight is 100 km,angle of incidence is 170°)

    圖9 給出了艙體解體后一些典型部件在飛行過程中的溫度、壓力、熱流等表面參數(shù)分布。 通過對(duì)這些部件組成材質(zhì)(金屬類)的熔融情況分析,可以得出其最終存活性結(jié)論及落區(qū)范圍。

    圖9 LTG 典型部件表面參數(shù)分布Fig.9 Surface parameter distribution of LTG typical components

    根據(jù)計(jì)算結(jié)果分析表明,該大型航天器飛行至高度120 km 以下時(shí),氣動(dòng)力/熱即產(chǎn)生明顯作用。 在高度約105 km 處,太陽電池帆板在氣動(dòng)力熱作用下脫落解體。 在高度約100 ~95 km 區(qū)間,由鋁合金構(gòu)成的兩艙組合體殼體熔融解體,內(nèi)部部件開始承受強(qiáng)烈氣動(dòng)力/熱作用。 在高度約95~70 km 區(qū)間,鋁合金及不銹鋼部件分別熔融解體;碳纖維及復(fù)合材料部件熱解充分而解體。耐高溫合金組件如軌控發(fā)動(dòng)機(jī)組和姿控發(fā)動(dòng)機(jī)組,其合金材料在高度約60 km 處基本完全熔融,留存可辨識(shí)原本初始形態(tài)的殘骸落至地面的可能性極小,通常是經(jīng)多次解體或熔融凝固殘骸、碎片墜落地面。 此外,兩艙解體后電池模塊在氣動(dòng)熱作用下將很快因溫升發(fā)生爆破解體。

    8 結(jié)論

    1) 提出并闡述了航天器再入隕落解體三層級(jí)模型(TPM)(含碎片/微粒模型BCSM)以及基于條件邊界的參數(shù)統(tǒng)計(jì)方法。 該三層級(jí)模型融合了面向物體法涉及碎片/微粒層級(jí)和面向航天器法設(shè)計(jì)系統(tǒng)/子系統(tǒng)和部件層級(jí)的思想。

    2) 針對(duì)航天器再入隕落解體過程仿真建模存在不確定性氣動(dòng)力/熱及飛行狀態(tài)參數(shù),發(fā)展了基于條件邊界的參數(shù)統(tǒng)計(jì)方法,通過相關(guān)軟件研制及對(duì)一個(gè)大型航天器實(shí)例的再入隕落解體過程分析,說明提出的模型和分析預(yù)報(bào)基本策略是合理可行的。

    3) 后續(xù)工作將繼續(xù)細(xì)化量化三層級(jí)模型、條件邊界及參數(shù)分布統(tǒng)計(jì)等設(shè)定。 航天器再入隕落解體過程分析預(yù)報(bào)的精準(zhǔn)度依賴于基礎(chǔ)研究成果的支撐程度,這是針對(duì)航天器再入隕落解體過程空氣動(dòng)力學(xué)相關(guān)科學(xué)問題深入研究的驅(qū)動(dòng)力。

    猜你喜歡
    隕落解體氣動(dòng)力
    行星UC2的隕落
    飛行載荷外部氣動(dòng)力的二次規(guī)劃等效映射方法
    《芳華》:事關(guān)理想主義的隕落
    蘇聯(lián)1991年解體前的最后時(shí)光
    中外文摘(2017年14期)2017-07-31 16:16:48
    側(cè)風(fēng)對(duì)拍動(dòng)翅氣動(dòng)力的影響
    “娃娃親”因兩家發(fā)展不同而解體
    美空軍又一退役氣象衛(wèi)星在軌解體
    太空探索(2016年12期)2016-07-18 11:13:43
    “東方之星”隕落長江全媒體報(bào)道體現(xiàn)大愛
    新聞前哨(2015年8期)2015-12-23 03:22:39
    隕落的夏代城市:大師姑城址發(fā)掘記
    大眾考古(2014年10期)2014-06-21 07:12:08
    高速鐵路接觸線覆冰后氣動(dòng)力特性的風(fēng)洞試驗(yàn)研究
    久久精品熟女亚洲av麻豆精品 | 三级国产精品片| 国产精品女同一区二区软件| 亚洲成人中文字幕在线播放| 一个人看的www免费观看视频| 色综合色国产| 国产成人午夜福利电影在线观看| 欧美激情国产日韩精品一区| 99热这里只有是精品50| 成人一区二区视频在线观看| 欧美xxxx黑人xx丫x性爽| 观看美女的网站| 国产免费视频播放在线视频 | 欧美bdsm另类| 亚洲电影在线观看av| av国产久精品久网站免费入址| 黄片wwwwww| 国产精品国产高清国产av| 97超碰精品成人国产| 中文字幕av在线有码专区| 午夜视频国产福利| 干丝袜人妻中文字幕| 男人狂女人下面高潮的视频| 麻豆乱淫一区二区| 国产乱来视频区| 99久国产av精品| 韩国高清视频一区二区三区| 秋霞在线观看毛片| 久久久久久久国产电影| 欧美xxxx性猛交bbbb| 春色校园在线视频观看| 日本一二三区视频观看| 亚洲电影在线观看av| 亚洲精品日韩在线中文字幕| 人人妻人人澡欧美一区二区| 国产极品天堂在线| 精品人妻视频免费看| 26uuu在线亚洲综合色| 97超视频在线观看视频| 国产成人福利小说| 中文字幕熟女人妻在线| 日本爱情动作片www.在线观看| 亚洲五月天丁香| 亚洲欧美中文字幕日韩二区| 欧美丝袜亚洲另类| 婷婷色av中文字幕| 99久国产av精品| 国产大屁股一区二区在线视频| 一级黄片播放器| 免费看日本二区| 乱码一卡2卡4卡精品| 成人午夜精彩视频在线观看| 国产成人午夜福利电影在线观看| 亚洲在久久综合| 国语自产精品视频在线第100页| 精品国产露脸久久av麻豆 | 国产高清三级在线| 成人综合一区亚洲| 亚洲精华国产精华液的使用体验| 熟妇人妻久久中文字幕3abv| 欧美激情久久久久久爽电影| 波野结衣二区三区在线| 夜夜看夜夜爽夜夜摸| 亚洲精品乱码久久久久久按摩| 免费播放大片免费观看视频在线观看 | 国产色爽女视频免费观看| 成人特级av手机在线观看| 亚洲最大成人av| 欧美成人一区二区免费高清观看| 成人av在线播放网站| 亚洲不卡免费看| 国产精品野战在线观看| 午夜免费激情av| 久久久久久国产a免费观看| 亚洲av一区综合| 免费观看性生交大片5| 好男人视频免费观看在线| 简卡轻食公司| 色视频www国产| av线在线观看网站| 成人综合一区亚洲| 欧美人与善性xxx| 国产高清不卡午夜福利| 人妻夜夜爽99麻豆av| 国产午夜精品久久久久久一区二区三区| 自拍偷自拍亚洲精品老妇| 男女下面进入的视频免费午夜| 国产av不卡久久| 国产精品女同一区二区软件| 丰满少妇做爰视频| 国产精品久久久久久久久免| 中文字幕亚洲精品专区| 夜夜爽夜夜爽视频| 日韩一区二区视频免费看| 国产伦理片在线播放av一区| 欧美最新免费一区二区三区| av天堂中文字幕网| 国产精品一区二区三区四区免费观看| 99视频精品全部免费 在线| 永久免费av网站大全| 精华霜和精华液先用哪个| 免费av毛片视频| 嫩草影院精品99| 能在线免费看毛片的网站| 亚洲成人久久爱视频| 成年免费大片在线观看| 99久久成人亚洲精品观看| 免费人成在线观看视频色| 日本熟妇午夜| 免费观看在线日韩| 亚洲精品aⅴ在线观看| 亚洲乱码一区二区免费版| 国产精品av视频在线免费观看| 美女国产视频在线观看| 有码 亚洲区| 性插视频无遮挡在线免费观看| 国产精华一区二区三区| 国产亚洲一区二区精品| 国产综合懂色| 欧美日本视频| 三级男女做爰猛烈吃奶摸视频| 22中文网久久字幕| 精品人妻视频免费看| 国产精品,欧美在线| 欧美性感艳星| 免费搜索国产男女视频| 麻豆成人午夜福利视频| 欧美另类亚洲清纯唯美| 一级毛片aaaaaa免费看小| 99九九线精品视频在线观看视频| 免费看a级黄色片| 日韩精品有码人妻一区| 床上黄色一级片| 最近中文字幕高清免费大全6| 插阴视频在线观看视频| 秋霞在线观看毛片| 午夜久久久久精精品| 日韩一区二区视频免费看| 三级男女做爰猛烈吃奶摸视频| 一区二区三区高清视频在线| 日本色播在线视频| 国产国拍精品亚洲av在线观看| 久久99热这里只频精品6学生 | 一本久久精品| 秋霞在线观看毛片| 观看免费一级毛片| 亚洲无线观看免费| 国产成人精品久久久久久| 91av网一区二区| 国产亚洲91精品色在线| 听说在线观看完整版免费高清| 一边摸一边抽搐一进一小说| 国产精品国产三级国产专区5o | 深夜a级毛片| 七月丁香在线播放| 成人午夜高清在线视频| 欧美色视频一区免费| 蜜桃久久精品国产亚洲av| 联通29元200g的流量卡| 伦精品一区二区三区| 精品一区二区三区人妻视频| 国产午夜精品久久久久久一区二区三区| 精品午夜福利在线看| 日韩av在线大香蕉| 国产黄色视频一区二区在线观看 | 精品久久久久久久久亚洲| 又粗又硬又长又爽又黄的视频| 日韩av在线大香蕉| 非洲黑人性xxxx精品又粗又长| 乱码一卡2卡4卡精品| 天天一区二区日本电影三级| 91久久精品电影网| av免费观看日本| 国产一区二区在线av高清观看| 最后的刺客免费高清国语| 免费看日本二区| 成人欧美大片| 三级男女做爰猛烈吃奶摸视频| 国产三级在线视频| 有码 亚洲区| 男女视频在线观看网站免费| 午夜激情欧美在线| 99热这里只有是精品50| 大香蕉久久网| 日本黄大片高清| 婷婷六月久久综合丁香| 中国国产av一级| 国产精品久久久久久精品电影小说 | 国产不卡一卡二| 深爱激情五月婷婷| 午夜爱爱视频在线播放| 亚洲国产精品国产精品| 看免费成人av毛片| 亚洲成人精品中文字幕电影| 国产亚洲91精品色在线| av国产久精品久网站免费入址| 美女高潮的动态| 国产在线男女| 我的老师免费观看完整版| 国产私拍福利视频在线观看| 岛国在线免费视频观看| 日韩一区二区三区影片| 夜夜爽夜夜爽视频| 国产白丝娇喘喷水9色精品| 中文欧美无线码| 成人特级av手机在线观看| 日韩亚洲欧美综合| 亚洲人成网站在线播| 久久精品国产亚洲av涩爱| 在线观看66精品国产| 国产精品无大码| 中文字幕熟女人妻在线| 成人性生交大片免费视频hd| 国产 一区 欧美 日韩| h日本视频在线播放| 少妇丰满av| 两个人视频免费观看高清| 男人狂女人下面高潮的视频| 午夜免费男女啪啪视频观看| 我的老师免费观看完整版| 在线观看一区二区三区| 亚洲av电影在线观看一区二区三区 | 亚洲在线观看片| 夜夜看夜夜爽夜夜摸| 国产精品av视频在线免费观看| 少妇裸体淫交视频免费看高清| 黄片无遮挡物在线观看| www.色视频.com| 嫩草影院入口| 一卡2卡三卡四卡精品乱码亚洲| 少妇的逼水好多| 日本免费一区二区三区高清不卡| 欧美精品国产亚洲| 男插女下体视频免费在线播放| 级片在线观看| 白带黄色成豆腐渣| 久久精品国产自在天天线| 国产黄片视频在线免费观看| 在线免费观看不下载黄p国产| 免费观看精品视频网站| 寂寞人妻少妇视频99o| 久久久久久久久久黄片| 一个人观看的视频www高清免费观看| 久久久国产成人精品二区| 国产伦在线观看视频一区| 天天躁夜夜躁狠狠久久av| 国产国拍精品亚洲av在线观看| 国产精华一区二区三区| 国产又黄又爽又无遮挡在线| 色网站视频免费| 草草在线视频免费看| 日韩av在线大香蕉| 成人亚洲精品av一区二区| 国产伦在线观看视频一区| 日韩av在线免费看完整版不卡| 99久久中文字幕三级久久日本| 国产成人a∨麻豆精品| 国产一级毛片在线| 欧美变态另类bdsm刘玥| 国产高潮美女av| 九色成人免费人妻av| 好男人在线观看高清免费视频| 中文在线观看免费www的网站| 亚洲国产高清在线一区二区三| 日韩欧美三级三区| 国内精品宾馆在线| 欧美97在线视频| 欧美精品一区二区大全| 国产av不卡久久| 人妻少妇偷人精品九色| 可以在线观看毛片的网站| 精品人妻偷拍中文字幕| 亚洲精品国产成人久久av| 建设人人有责人人尽责人人享有的 | 国国产精品蜜臀av免费| 好男人在线观看高清免费视频| 国产一区有黄有色的免费视频 | 看黄色毛片网站| 久久99精品国语久久久| 69人妻影院| 午夜福利成人在线免费观看| 欧美不卡视频在线免费观看| 少妇被粗大猛烈的视频| 亚洲婷婷狠狠爱综合网| 色尼玛亚洲综合影院| 欧美区成人在线视频| 又粗又硬又长又爽又黄的视频| 成人性生交大片免费视频hd| 欧美97在线视频| 两个人的视频大全免费| 少妇的逼好多水| 狂野欧美白嫩少妇大欣赏| 91狼人影院| 中文欧美无线码| 国产黄片视频在线免费观看| 亚洲国产成人一精品久久久| 伦理电影大哥的女人| 亚洲av.av天堂| 18禁在线无遮挡免费观看视频| 亚洲精品一区蜜桃| 精品熟女少妇av免费看| 高清日韩中文字幕在线| 大香蕉97超碰在线| 日韩一区二区视频免费看| 亚洲精品亚洲一区二区| 一级二级三级毛片免费看| 如何舔出高潮| 日韩,欧美,国产一区二区三区 | 美女高潮的动态| 欧美激情国产日韩精品一区| 哪个播放器可以免费观看大片| 午夜亚洲福利在线播放| 午夜精品国产一区二区电影 | 国产伦精品一区二区三区视频9| 亚洲国产最新在线播放| 亚洲在线观看片| 汤姆久久久久久久影院中文字幕 | 久久热精品热| a级毛片免费高清观看在线播放| 秋霞在线观看毛片| 国产精品.久久久| 99久久中文字幕三级久久日本| 久久久国产成人免费| 亚洲成av人片在线播放无| 国产一区二区三区av在线| 99久久无色码亚洲精品果冻| 国产成人精品一,二区| 又粗又爽又猛毛片免费看| 我要搜黄色片| 内地一区二区视频在线| 亚洲第一区二区三区不卡| 欧美性猛交黑人性爽| 国产精品久久久久久久久免| 波野结衣二区三区在线| 一边亲一边摸免费视频| 熟女人妻精品中文字幕| 亚洲精华国产精华液的使用体验| 在线观看一区二区三区| 99久久精品国产国产毛片| 欧美日本视频| 国产不卡一卡二| 美女cb高潮喷水在线观看| 亚洲av中文字字幕乱码综合| 黄片无遮挡物在线观看| videossex国产| 可以在线观看毛片的网站| 91aial.com中文字幕在线观看| 亚洲欧美日韩卡通动漫| 国产高清三级在线| 神马国产精品三级电影在线观看| 一夜夜www| 日韩欧美国产在线观看| 国产精品综合久久久久久久免费| 一本久久精品| 男女那种视频在线观看| 欧美一区二区精品小视频在线| 日韩成人av中文字幕在线观看| 亚洲精品一区蜜桃| 纵有疾风起免费观看全集完整版 | 久久99蜜桃精品久久| 一个人观看的视频www高清免费观看| 国产成人a区在线观看| 婷婷色麻豆天堂久久 | 亚洲在久久综合| 国产精品一及| 日韩 亚洲 欧美在线| 天天一区二区日本电影三级| 爱豆传媒免费全集在线观看| 99热全是精品| 三级男女做爰猛烈吃奶摸视频| 九色成人免费人妻av| 女人被狂操c到高潮| 中文亚洲av片在线观看爽| 九九热线精品视视频播放| 亚洲人成网站在线播| 永久网站在线| 欧美一区二区亚洲| 日本三级黄在线观看| 亚洲欧美精品自产自拍| 国产黄色小视频在线观看| 男女那种视频在线观看| 自拍偷自拍亚洲精品老妇| 精品久久久久久久末码| 国产精品伦人一区二区| 永久网站在线| 美女国产视频在线观看| 长腿黑丝高跟| 国产精品女同一区二区软件| 色哟哟·www| 日本av手机在线免费观看| 国产成人福利小说| 国产伦精品一区二区三区视频9| 亚洲精品国产成人久久av| 在线免费观看不下载黄p国产| 91狼人影院| 久久婷婷人人爽人人干人人爱| 亚洲国产精品成人久久小说| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲av男天堂| 久久精品国产亚洲网站| 国产精品一区二区在线观看99 | 全区人妻精品视频| 国国产精品蜜臀av免费| 亚洲国产成人一精品久久久| 日本-黄色视频高清免费观看| 欧美性猛交╳xxx乱大交人| 亚洲精品一区蜜桃| 又粗又爽又猛毛片免费看| 91久久精品国产一区二区成人| 国产激情偷乱视频一区二区| 五月伊人婷婷丁香| 国产精品熟女久久久久浪| 国产伦在线观看视频一区| 国产私拍福利视频在线观看| 青春草国产在线视频| 成人亚洲欧美一区二区av| av在线天堂中文字幕| 亚洲人成网站高清观看| 久久欧美精品欧美久久欧美| 天天躁日日操中文字幕| 成人欧美大片| 国产一区二区三区av在线| 中文乱码字字幕精品一区二区三区 | 女人十人毛片免费观看3o分钟| 精品久久久噜噜| 五月伊人婷婷丁香| 好男人在线观看高清免费视频| 国产伦理片在线播放av一区| 久久这里只有精品中国| 在线免费观看不下载黄p国产| 亚洲成人av在线免费| 国产亚洲最大av| 中文在线观看免费www的网站| 日韩人妻高清精品专区| 国产精品久久久久久av不卡| av专区在线播放| 少妇被粗大猛烈的视频| 一本一本综合久久| 日韩,欧美,国产一区二区三区 | 亚洲av不卡在线观看| 亚洲精品aⅴ在线观看| av在线播放精品| 赤兔流量卡办理| 亚洲国产精品专区欧美| 51国产日韩欧美| 看十八女毛片水多多多| 亚洲最大成人av| videos熟女内射| 日韩一区二区三区影片| 美女被艹到高潮喷水动态| 26uuu在线亚洲综合色| 99热6这里只有精品| 午夜福利在线观看免费完整高清在| 精品免费久久久久久久清纯| 成人一区二区视频在线观看| 色综合亚洲欧美另类图片| 成年免费大片在线观看| 日本黄色视频三级网站网址| 身体一侧抽搐| 日韩制服骚丝袜av| 桃色一区二区三区在线观看| 边亲边吃奶的免费视频| 亚洲精品影视一区二区三区av| 亚洲成av人片在线播放无| 亚洲av成人av| av黄色大香蕉| 精华霜和精华液先用哪个| a级毛色黄片| 成人毛片a级毛片在线播放| 一级毛片我不卡| 伦理电影大哥的女人| 欧美性感艳星| 精品久久久久久久人妻蜜臀av| 丰满人妻一区二区三区视频av| 晚上一个人看的免费电影| 色综合站精品国产| 午夜久久久久精精品| 少妇丰满av| 国产亚洲一区二区精品| 三级国产精品欧美在线观看| 日本av手机在线免费观看| 小蜜桃在线观看免费完整版高清| 26uuu在线亚洲综合色| 色视频www国产| 午夜福利高清视频| 国模一区二区三区四区视频| 精品酒店卫生间| 床上黄色一级片| 精品久久久久久久久久久久久| 97人妻精品一区二区三区麻豆| 午夜爱爱视频在线播放| 亚洲国产高清在线一区二区三| 亚洲国产精品成人综合色| 日韩在线高清观看一区二区三区| 亚洲av免费在线观看| 内地一区二区视频在线| 欧美日本亚洲视频在线播放| 国产视频内射| 神马国产精品三级电影在线观看| 中文乱码字字幕精品一区二区三区 | 国产在线一区二区三区精 | 六月丁香七月| 国产成人freesex在线| 乱人视频在线观看| 欧美97在线视频| av又黄又爽大尺度在线免费看 | 黄色日韩在线| 日日摸夜夜添夜夜爱| 欧美人与善性xxx| av在线天堂中文字幕| 国产成人午夜福利电影在线观看| 3wmmmm亚洲av在线观看| 蜜桃亚洲精品一区二区三区| 欧美又色又爽又黄视频| 久久久国产成人免费| .国产精品久久| 国产免费视频播放在线视频 | 日日摸夜夜添夜夜添av毛片| 一边亲一边摸免费视频| 日韩视频在线欧美| 亚洲av免费在线观看| 尤物成人国产欧美一区二区三区| 婷婷六月久久综合丁香| 蜜桃久久精品国产亚洲av| 婷婷色综合大香蕉| 又爽又黄无遮挡网站| 国产亚洲5aaaaa淫片| 精品无人区乱码1区二区| 亚洲精品自拍成人| 51国产日韩欧美| 波多野结衣高清无吗| 国产精品国产三级国产专区5o | 一本一本综合久久| 中文亚洲av片在线观看爽| 亚洲综合精品二区| 国产av码专区亚洲av| av国产久精品久网站免费入址| av天堂中文字幕网| 欧美成人一区二区免费高清观看| 国产精品一区二区在线观看99 | 免费电影在线观看免费观看| 高清在线视频一区二区三区 | 国模一区二区三区四区视频| 在线a可以看的网站| 精品久久久久久久末码| 国语自产精品视频在线第100页| 91精品伊人久久大香线蕉| 国产毛片a区久久久久| 毛片一级片免费看久久久久| 波野结衣二区三区在线| 国语对白做爰xxxⅹ性视频网站| 欧美三级亚洲精品| 国产亚洲av嫩草精品影院| 你懂的网址亚洲精品在线观看 | 白带黄色成豆腐渣| 国产精品久久视频播放| 色网站视频免费| 最近手机中文字幕大全| 在线天堂最新版资源| 午夜免费激情av| 国产极品精品免费视频能看的| 一级毛片我不卡| 国产av码专区亚洲av| 日本欧美国产在线视频| 午夜福利网站1000一区二区三区| 99在线视频只有这里精品首页| 精品一区二区三区人妻视频| 国产片特级美女逼逼视频| 日韩欧美 国产精品| 亚洲欧美一区二区三区国产| 欧美又色又爽又黄视频| 久久精品久久久久久久性| av免费观看日本| 成年av动漫网址| 精品一区二区免费观看| 久久久久久久亚洲中文字幕| 国产欧美另类精品又又久久亚洲欧美| 国产高清三级在线| 国产av不卡久久| 成人午夜精彩视频在线观看| 国产高清国产精品国产三级 | 亚洲精华国产精华液的使用体验| 我要看日韩黄色一级片| 内射极品少妇av片p| 99国产精品一区二区蜜桃av| 小说图片视频综合网站| 日本一二三区视频观看| 三级国产精品欧美在线观看| 国产精品美女特级片免费视频播放器| 午夜日本视频在线| 又爽又黄无遮挡网站| 热99在线观看视频| av在线老鸭窝| 97在线视频观看| 高清视频免费观看一区二区 | 91久久精品国产一区二区成人| 免费av不卡在线播放| 99久久精品国产国产毛片| 日本免费一区二区三区高清不卡| 久久久久性生活片| 26uuu在线亚洲综合色| 久久精品国产亚洲av天美| 我的女老师完整版在线观看| 国产91av在线免费观看| 日本五十路高清| 又爽又黄无遮挡网站| 中文欧美无线码| 国产片特级美女逼逼视频| 亚洲第一区二区三区不卡| 最近2019中文字幕mv第一页| 国产精品国产三级专区第一集| 亚洲精品456在线播放app| 精品一区二区免费观看| 久久精品熟女亚洲av麻豆精品 | 91午夜精品亚洲一区二区三区|