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

    航天器火工沖擊模擬試驗及響應(yīng)預(yù)示方法研究綜述*

    2016-04-20 10:20:31丁繼鋒韓增堯鄒元杰
    爆炸與沖擊 2016年2期
    關(guān)鍵詞:火工火工品航天器

    趙 欣,丁繼鋒,韓增堯,鄒元杰

    (1.中國空間技術(shù)研究院總體部,北京 100094;2.中國空間技術(shù)研究院研究發(fā)展部,北京 100094)

    航天器火工沖擊模擬試驗及響應(yīng)預(yù)示方法研究綜述*

    趙 欣1,丁繼鋒1,韓增堯2,鄒元杰1

    (1.中國空間技術(shù)研究院總體部,北京 100094;2.中國空間技術(shù)研究院研究發(fā)展部,北京 100094)

    航天器火工沖擊力學(xué)環(huán)境是由星箭分離、部組件展開等工作過程中的火工品起爆引起的作用于結(jié)構(gòu)上的高頻、高加速度量級的瞬態(tài)沖擊響應(yīng),能對航天器上含有晶振、脆性材料等的精密電子設(shè)備造成致命損傷,是航天器需要經(jīng)歷的最苛刻的力學(xué)環(huán)境之一。本文中,對國內(nèi)外航天器火工沖擊地面試驗方法和環(huán)境預(yù)示方法做了全面、詳細的介紹,總結(jié)了這兩個方面的研究進展,分析了我國在這兩個方面與航天強國的差距。最后,從我國航天工程實際需求出發(fā),提出了今后航天器火工沖擊領(lǐng)域應(yīng)重點開展的研究方向。

    振動理論;響應(yīng)預(yù)示;火工沖擊;航天器;地面試驗;流體代碼;顯式有限元;統(tǒng)計能量分析;虛模態(tài)綜合法

    航天器在發(fā)射、在軌運行過程中會經(jīng)歷嚴酷復(fù)雜的振動、噪聲和沖擊等力學(xué)環(huán)境?;鸸_擊力學(xué)環(huán)境是沖擊力學(xué)環(huán)境的一種,是由星箭分離、部組件展開等工作過程中的火工品起爆引起的作用于結(jié)構(gòu)上的瞬態(tài)沖擊響應(yīng),具有高頻率、高加速度幅值和高應(yīng)變率的特點。航天器結(jié)構(gòu)及有效載荷可能會在火工沖擊力學(xué)環(huán)境中損傷或失效,這將影響甚至提前終止航天任務(wù)[1]。

    火工沖擊過程中,火工品爆炸產(chǎn)生的載荷以應(yīng)力波的形式在衛(wèi)星結(jié)構(gòu)體中傳播,并在衛(wèi)星結(jié)構(gòu)體產(chǎn)生瞬時、高頻、高加速度幅值的振動響應(yīng)。火工沖擊的加速度響應(yīng)幅值可高達200 000g,作用持續(xù)時間小于20 ms,主要頻率在100 Hz~100 kHz范圍[2-3]?;鸸_擊載荷對航天器主結(jié)構(gòu)通常不會構(gòu)成破壞性影響,但是能夠造成對沖擊敏感的元器件和儀器設(shè)備的損傷,如高頻環(huán)境敏感的高靈敏度設(shè)備、電子產(chǎn)品,他們會發(fā)生故障,輕薄結(jié)構(gòu)、脆性材料會發(fā)生破壞等[4]。美國國家航天局(NASA)早在20世紀70~80年代已經(jīng)認識到該問題,并投入大量人力物力[5]深入研究了火工沖擊產(chǎn)生和傳遞的機理、沖擊響應(yīng)的預(yù)示方法、試驗?zāi)M技術(shù),并制定了相應(yīng)標準[1,6-7],規(guī)范了航天領(lǐng)域火工沖擊相關(guān)的試驗、設(shè)計和預(yù)示方法,目前已經(jīng)能夠?qū)鸸_擊導(dǎo)致的飛行故障進行有效控制。

    隨著我國航天事業(yè)的飛速發(fā)展,航天器上需要采用火工解鎖完成展開任務(wù)的組件不斷增加,精密電子設(shè)備對沖擊載荷也更加敏感[8]。因此,對航天器火工沖擊力學(xué)環(huán)境的預(yù)示、試驗條件的制定、地面的模擬試驗方法提出了更高的要求。

    依據(jù)沖擊響應(yīng)的強度和頻譜范圍,可以將航天器火工沖擊力學(xué)環(huán)境分為近場、中場和遠場[1-2]。在設(shè)計中,要求近場范圍內(nèi)不得布置沖擊敏感設(shè)備,而設(shè)備必須通過沖擊試驗的考核才能布置在中場和遠場上[2]。

    1 地面模擬試驗方法研究

    航天工程型號中,為了獲得火工爆炸時在結(jié)構(gòu)上產(chǎn)生的真實沖擊響應(yīng),一般需利用真實的火工品和結(jié)構(gòu)進行沖擊試驗。采用真實裝備的火工沖擊試驗成本高,進行各種類型試驗往往是現(xiàn)實條件所不允許的。因此,從20世紀70~80年代開始,在美國率先研制了多種火工沖擊模擬試驗方法[9]。這些方法大體可分為火工爆炸式和非火工爆炸式式兩類,應(yīng)用這些方法能夠在極大程度上模擬真實火工品爆炸所產(chǎn)生的沖擊響應(yīng)譜,在航天工程中獲得了廣泛應(yīng)用。此后,歐洲等國家也都大體沿用此方法研制自己的試驗裝備。

    1.1 模擬試驗的考核方法

    為了衡量模擬試驗在多大程度上能夠反應(yīng)真實火工沖擊的載荷環(huán)境,需要找尋一種火工沖擊的描述方式,或者稱為考核方式。用產(chǎn)品的時域響應(yīng)考核模擬試驗的效果是過于苛刻的,在航天工程中,對火工沖擊的沖擊響應(yīng)譜描述獲得了最廣泛的認可[1,6,10]。

    1.2 模擬試驗方法

    1.2.1 火工爆炸式

    火工爆炸式模擬試驗裝置是直接采用火藥以一定的裝藥方式裝載在試驗臺上,火藥通常直接與試驗臺接觸,也可能采用附加的裝藥裝置使火藥與試驗臺不發(fā)生直接接觸。可通過電起爆器或一根簡單的引線將火藥引爆,產(chǎn)生與實際火工品爆炸相似的沖擊激勵。一種典型的火工爆炸式?jīng)_擊模擬試驗裝置如圖1所示,試驗臺由4根尼龍橡膠帶懸掛起來,布置在試驗臺底面的柔性聚能炸藥條位置可以調(diào)節(jié)。被測試的儀器按照實際在航天器上安裝的方式安裝在試驗臺上。通常需要先安裝儀器模擬裝置進行反復(fù)試驗,以確定炸藥條的合適位置,然后安裝實際被測儀器進行正式的試驗[9]。

    采用火工爆炸式的模擬試驗裝置能夠模擬真實的火工品爆炸產(chǎn)生的高頻譜和高加速度量級,因此在航天器系統(tǒng)級試驗中有重要的應(yīng)用價值。歐洲阿利安太空公司精心設(shè)計了火工爆炸式模擬試驗裝置SHOGUN(shock generation unit)裝置[12](見圖2)。該裝置安裝于衛(wèi)星結(jié)構(gòu)和星箭減振轉(zhuǎn)接接頭下方,用于模擬星箭轉(zhuǎn)接接頭與火箭上級分離的點火沖擊。歐洲的SPOT1衛(wèi)星、SPOT5衛(wèi)星及Alphabus等多顆衛(wèi)星都使用該裝置進行與運載火箭的分離沖擊環(huán)境試驗[12-14]。在火工爆炸近場安裝有敏感設(shè)備的情況下,火工爆炸式試驗裝置還被用來模擬火工爆炸近場的沖擊環(huán)境。法國Alcatel ETCA公司研制的火工沖擊模擬試驗裝置可模擬任何類型沖擊響應(yīng)譜[15-16],目前在多臺電子產(chǎn)品的沖擊試驗中得到成功應(yīng)用[17]。北京強度環(huán)境研究所研制了類似的試驗裝置[18],目前該裝置還處于工程應(yīng)用研究階段。

    采用火工爆炸式模擬試驗裝置,在正式試驗前需要進行多次試錯才能得到需要的沖擊響應(yīng)譜,且試驗結(jié)果的重復(fù)性較差。因此,需要特別注意試驗的規(guī)程和細節(jié),比如存放炸藥的溝槽要求經(jīng)過機械加工使表面光滑,以此避免由于表面空穴引起的能量吸收和爆炸方向的改變。這種試驗方法由于使用了炸藥,存在安全隱患,因此需要采取安全措施以確保試驗人員的安全。

    圖1 典型火工爆炸式?jīng)_擊模擬試驗裝置[11]Fig.1 Ordnance-generated pyroshock simulator

    圖2 2010年Alphabus衛(wèi)星SHOGUN試驗現(xiàn)場[12]Fig.2 SHOGUN test spot for Alphabus in 2010

    1.2.2 非火工爆炸式

    非火工爆炸式試驗裝置主要有機械試驗臺、電動振動臺兩種。機械式試驗臺采用擺錘或氣槍等作為沖擊發(fā)生器,作用在一個諧振裝置上激發(fā)其響應(yīng),諧振裝置上安裝有被測儀器,在被測儀器上能夠產(chǎn)生與火工沖擊十分接近的沖擊響應(yīng)譜。實際的諧振結(jié)構(gòu)可以采用諧振桿、諧振梁、諧振板或多層板等結(jié)構(gòu)形式[9,19](見圖3)。

    MIPS(mechanical impulse pyroshock simulator)是在航天領(lǐng)域被廣泛使用的一種試驗裝置,該裝置可以通過調(diào)節(jié)撞擊質(zhì)量塊的質(zhì)量、材料和速度、撞擊位置、被測試儀器的位置及改變夾具和懸吊結(jié)構(gòu)等控制沖擊譜形[9](見圖4~5)。實際試驗過程中為了產(chǎn)生要求的沖擊響應(yīng)譜,需反復(fù)調(diào)試各種參數(shù)。

    通過精心設(shè)計諧振結(jié)構(gòu)可使調(diào)試次數(shù)有效減少??烧{(diào)節(jié)諧振梁裝置可實現(xiàn)沖擊響應(yīng)譜的拐點頻率的連續(xù)變化,它的主體由兩個夾具將一個矩形截面共振梁固定在試驗臺上,通過調(diào)節(jié)兩個夾具的位置來調(diào)節(jié)諧振梁的長度和固有頻率,實現(xiàn)被測量儀器上的沖擊響應(yīng)譜拐點頻率在250~7 000 Hz范圍內(nèi)連續(xù)變化[9](見圖6)。

    圖3 機械式?jīng)_擊模擬試驗裝置示意圖[11]Fig.3 Resonant fixtures for MIPS

    圖4 MIPS試驗臺原理圖[9]Fig.4 MIPS simulator

    圖5 意大利研制的MIPS模擬試驗裝置[20]Fig.5 A mechanical impulse pyroshock simulator in Italy

    圖6 美國研制的可調(diào)節(jié)諧振梁裝置[21]Fig.6 A tunable resonant beam device in America

    機械撞擊式?jīng)_擊試驗一般可以模擬頻率在10 kHz以內(nèi)、幅值不超過5 000g的沖擊響應(yīng)譜,主要用于模擬火工沖擊中、遠場的響應(yīng)。因重復(fù)性好、試驗成本低,航天型號中進行沖擊試驗時,主要采用該方法。但是,采用該方法時也會產(chǎn)生明顯的過試驗問題[2]。

    電動式振動臺或電液式振動臺也可以用來模擬量級不是很高的火工沖擊環(huán)境,經(jīng)過控制系統(tǒng)綜合處理,電動振動臺能夠在被測試設(shè)備上產(chǎn)生接近于火工沖擊響應(yīng)譜的加速度振動響應(yīng)。這種方法具有操作成本低、可控性高等優(yōu)點,但是該設(shè)備僅能實現(xiàn)單向加載,且能夠達到的加速度幅值較低,頻譜范圍較小(3 000 Hz以下),在航天工程型號中多用于元器件篩查和量級較低的沖擊試驗中[9]。

    1.3 數(shù)據(jù)測量、采集和分析

    火工沖擊的典型測量系統(tǒng)通常包含傳感器、導(dǎo)線、信號調(diào)制器、濾波器、模數(shù)轉(zhuǎn)換器和數(shù)據(jù)記錄裝置等。由于火工沖擊產(chǎn)生的加速度量級超過一般加速度傳感器的量程,因而需要采用特制的加速度傳感器,用以防止峰值加速度達到晶體極限而產(chǎn)生的噪聲和零位漂移現(xiàn)象。除此以外,必須特別關(guān)注測量系統(tǒng)的一些細節(jié),如導(dǎo)線的連接方式、數(shù)模轉(zhuǎn)換裝置的頻率響應(yīng)、零位漂移、頻率截斷和動態(tài)量程等[22]。即便試驗的數(shù)據(jù)測量和采集過程特別仔細,測得數(shù)據(jù)仍可能具有很大的離散性,甚至有一些數(shù)據(jù)有本質(zhì)錯誤,因此需要對測量數(shù)據(jù)的有效性進行甄別。一般采用Powers-Piersol程序判斷爆炸沖擊數(shù)據(jù)的有效性,該方法對加速度時間歷程進行一次或兩次積分,并比較SRS的正譜和負譜[23]。如2.1節(jié)中所述,采用SRS對測量結(jié)果進行分析,SRS的計算算法也將影響對測量結(jié)果的評估,采用Smallwood算法可以減少算法引起的誤差[7,24]。另外,在對大量地面測量數(shù)據(jù)進行分析,來確定極限飛行動力學(xué)環(huán)境,或者對結(jié)構(gòu)響應(yīng)進行預(yù)示的情況下,需要對有效數(shù)據(jù)進行必要的統(tǒng)計學(xué)分析[13,25]。

    2 航天器火工沖擊環(huán)境預(yù)示方法

    航天器火工沖擊力學(xué)環(huán)境問題是一個寬頻、瞬態(tài)和強非線性的沖擊動力學(xué)問題。目前,工程上采用的和在研的火工沖擊響應(yīng)預(yù)示方法有兩類:一類是通過工程上大量試驗數(shù)據(jù)分析的結(jié)果估計結(jié)構(gòu)上沖擊載荷的量級,即試驗外推方法;另一類是應(yīng)用有限元方法、統(tǒng)計能量分析及模態(tài)綜合等方法進行結(jié)構(gòu)響應(yīng)預(yù)示分析,即數(shù)值與分析模型方法。

    2.1 試驗外推方法

    由于難以采用傳統(tǒng)的動力學(xué)分析方法準確預(yù)示沖擊載荷的結(jié)構(gòu)響應(yīng),工程中常采用試驗外推方法,即結(jié)合工程經(jīng)驗與相似試驗的數(shù)據(jù)估計結(jié)構(gòu)上沖擊載荷的量級。近年來,沖擊響應(yīng)的試驗外推方法取得了一定的進展,較有效的方法主要有3種:經(jīng)驗?zāi)P头?、?shù)據(jù)外推法和子結(jié)構(gòu)路徑外推法。其中,前2種方法由美國NASA基于經(jīng)驗數(shù)據(jù)提出,第3種方法由法國宇航環(huán)境工程試驗中心(INTESPACE)和法國國家空間中心(CNES)聯(lián)合提出。

    2.1.1 經(jīng)驗?zāi)P头?/p>

    經(jīng)驗?zāi)P头╗1]能夠根據(jù)火工沖擊源的類型以及沖擊傳遞的距離、過連接界面的數(shù)目,快速、粗略地估計航天器結(jié)構(gòu)上某一區(qū)域的沖擊響應(yīng)譜。該方法無需建立結(jié)構(gòu)和火工品的詳細模型,不需要結(jié)構(gòu)設(shè)計的細節(jié),即可快速估計結(jié)構(gòu)上的沖擊響應(yīng),但是預(yù)示精度較低,可以在航天器設(shè)計的方案階段為火工品的選擇、結(jié)構(gòu)構(gòu)型方案的確定提供重要參考。文獻[1]中,給出了不同類型點源附近(10~15 cm)實測沖擊響應(yīng)譜隨距離衰減的關(guān)系。一般沖擊載荷每經(jīng)過一個界面沖擊響應(yīng)譜峰值衰減約40%,但界面總數(shù)不能超過3個,且沖擊譜的斜線部分不衰減。

    2.1.2 數(shù)據(jù)外推法

    應(yīng)用航天工程中已經(jīng)積累的大量火工沖擊數(shù)據(jù),一個新型航天器結(jié)構(gòu)上的沖擊響應(yīng)可以根據(jù)一個參考航天器上的實測響應(yīng)外推得到,此種沖擊響應(yīng)預(yù)示方法即為數(shù)據(jù)外推法。

    文獻[1]中,給出了一種外推法的實施步驟和計算公式。新型航天器結(jié)構(gòu)上的沖擊響應(yīng)譜,可以由一個參考航天器上距爆炸源相同距離的實測沖擊譜計算得到。數(shù)據(jù)外推法的準確性,取決于2種航天器采用的火工品和結(jié)構(gòu)構(gòu)型的相似程度?;鸸て窞橥惍a(chǎn)品、結(jié)構(gòu)構(gòu)型相似時,該方法的準確性較高[2]。歐美各國較多采用該方法進行航天器設(shè)計和沖擊環(huán)境預(yù)示[26-28]。

    2.1.3 子結(jié)構(gòu)路徑外推法

    子結(jié)構(gòu)路徑外推法先對衛(wèi)星整體結(jié)構(gòu)進行特定規(guī)則的子結(jié)構(gòu)分解和歸類,并確定傳遞路徑的定義方法。對衛(wèi)星結(jié)構(gòu)上任意一點,有確定的傳遞路徑。任意一點的響應(yīng),可以用傳遞路徑上的傳感器測點數(shù)據(jù)簡化結(jié)果加上一個連接影響函數(shù)(與傳遞路徑上的連接數(shù)相關(guān))和一個距離影響函數(shù)(與承力筒徑向距離ΔR和軸向距離ΔX相關(guān))得到[13](見圖7)。

    圖7 Walls結(jié)構(gòu)響應(yīng)路徑衰減示意圖[20]Fig.7 Level decrease in the Walls-structure

    子結(jié)構(gòu)路徑外推法由法國等歐洲國家的宇航部門和公司聯(lián)合提出,并實施于SPOT5衛(wèi)星平臺。經(jīng)過SHOGUN試驗數(shù)據(jù)的檢驗,證明該方法可以較準確地估計衛(wèi)星結(jié)構(gòu)與設(shè)備連接界面上任意一點的沖擊響應(yīng)[13]。但是應(yīng)用該方法進行沖擊響應(yīng)預(yù)示時,特定型號衛(wèi)星數(shù)據(jù)的處理結(jié)果(計算公式)只適用于該衛(wèi)星的響應(yīng)預(yù)示;在進行其他衛(wèi)星型號響應(yīng)預(yù)示時,需要對衛(wèi)星結(jié)構(gòu)進行重新分解,連接結(jié)構(gòu)的影響系數(shù)需要重新估算,測量結(jié)果的響應(yīng)譜簡化公式的參數(shù)也有所不同。

    2.1.4 小結(jié)

    采用試驗外推方法能夠給出沖擊環(huán)境的大致估計,但是依據(jù)這種方法進行航天器及其組件沖擊環(huán)境條件設(shè)計往往過于苛刻,或者需要未知型號航天器和已有試驗數(shù)據(jù)的航天器在結(jié)構(gòu)上具有很大的相似性。表1是各種方法特點的對比。

    表1 試驗外推方法比較Table 1 Comparison of different extrapolations

    2.2 數(shù)值與分析模型方法

    火工沖擊載荷具有以高頻為主的寬頻帶、瞬態(tài)和強非線性的特點,用分析方法準確地預(yù)示沖擊載荷引起的結(jié)構(gòu)響應(yīng)具有相當?shù)碾y度。但是相比于試驗方法,應(yīng)用數(shù)值與分析模型方法的硬件成本低、能夠獲取的數(shù)據(jù)量大且詳細、并且能夠先于結(jié)構(gòu)的制造裝配給出響應(yīng)結(jié)果。因而,考慮到航天器及其組件沖擊環(huán)境條件設(shè)計及結(jié)構(gòu)的協(xié)同設(shè)計等,對沖擊響應(yīng)數(shù)值及分析方法的研究仍具有重要的意義。

    2.2.1 Hydrocodes

    Hydrocodes[29]是能夠解決高頻、瞬態(tài)加載條件下流固耦合問題的顯式有限元或有限差分計算程序,應(yīng)用這類程序能夠?qū)Ρ?燃燒過程、非線性結(jié)構(gòu)變形/分離以及應(yīng)力波的產(chǎn)生和傳遞特性,進行精細建模與分析[1]。Hydrocodes的核心是它的計算框架,而不拘泥于某種特定的方法。首先,將表示連續(xù)系統(tǒng)的物理方程在時間和空間上離散,將變量的時間相關(guān)性和空間相關(guān)性分離開來,然后,根據(jù)質(zhì)量、能量和動量守恒建立控制方程求解[30]。Hydrocodes采用傳統(tǒng)的有限元法或有限差分法進行空間離散,為了能夠捕捉應(yīng)力波傳播的細節(jié),在時域上求解時必須采用顯式的差分方程進行。J.A.Zukas[29]系統(tǒng)介紹了Hydrocodes的計算框架以及所包含的理論知識和方法體系。

    典型的Hydrocodes程序有Zeus、EPIC、DYNA和CTH等,基于DYNA程序的商用軟件LS-DYNA已被廣泛應(yīng)用于爆炸沖擊動力學(xué)分析。理論上,該類Hydrocodes程序可求解各類沖擊問題,但是由于精細建模所帶來的時間成本和人力成本較高,該方法應(yīng)用于復(fù)雜航天器火工沖擊響應(yīng)的預(yù)示受到了一定的限制。

    2.2.2 顯式有限元法

    動力有限元法可以對彈性系統(tǒng)的低頻動力問題給出近似的解答,它通過變分原理和空間有限元離散得到一個近似的彈性動力學(xué)微分方程組,然后通過差分法(如顯式的中心差分法和隱式的Newmark法)得到t+Δt、t-Δt、t-2Δt等時刻的運動量之間的遞推關(guān)系[31]。

    隱式有限元法對于穩(wěn)態(tài)、低頻的動力系統(tǒng)能夠給出較準確的解答,但是不適合求解瞬態(tài)、高頻問題。顯式有限元法的臨界步長較小,在每一個時間步中它所需要的計算量也比隱式方法小得多。由于顯式有限元法能夠較好地模擬應(yīng)力波的傳播特點[30],以中心差分法為代表的顯式有限元法被廣泛應(yīng)用于固體中應(yīng)力波傳播等問題的求解。不同于Hydrocodes的是,通常所說的顯式有限元法基于變分原理建立控制方程,適用于求解彈性小變形問題。在航天器火工沖擊力學(xué)環(huán)境問題的求解中,顯式有限元法最適合求解近場響應(yīng),此處的結(jié)構(gòu)響應(yīng)主要由應(yīng)力波傳播決定,且火工品爆炸剩余能量有限,近場結(jié)構(gòu)變形屬于彈性變形。工程上也有應(yīng)用有限元方法進行整星計算的研究,如文獻[32],但是計算結(jié)果與試驗結(jié)果差距較大。

    2.2.3 統(tǒng)計能量分析

    統(tǒng)計能量分析(statistical energy analysis,SEA)方法的提出源于航空航天器工程中的“聲振”問題,它從能量的觀點分析結(jié)構(gòu)振動與聲場間的耦合關(guān)系,并應(yīng)用統(tǒng)計平均的概念,克服了高頻、密模態(tài)造成的響應(yīng)隨機性難以預(yù)示的問題。目前,該方法是解決復(fù)雜系統(tǒng)高頻動力學(xué)問題的一個有力的工具[33]。

    SEA方法將一個復(fù)雜的系統(tǒng)劃分為若干子系統(tǒng),以各子系統(tǒng)輸入、輸出和耗散的“能量”為基本自變量,建立子系統(tǒng)功率流平衡方程。利用功率流平衡方程來描述耦合子系統(tǒng)間的相互作用。它應(yīng)用統(tǒng)計平均的概念,得到每個子系統(tǒng)在頻域上的能量響應(yīng)平均值,該能量平均值可以轉(zhuǎn)化為速度、聲壓級等,用于表征系統(tǒng)的動力學(xué)響應(yīng)。

    瞬態(tài)統(tǒng)計能量分析(transient SEA,TSEA),將傳統(tǒng)的統(tǒng)計能量分析方法在穩(wěn)態(tài)激勵下的響應(yīng)分析擴展至沖擊或瞬態(tài)情況。J.E.Manning等[34]在穩(wěn)態(tài)能量流平衡方程中,加進了隨時間變化的子系統(tǒng)能量和耦合損耗因子,用于表征隨時間變化的子系統(tǒng)間的能量流動。同時,提出了用穩(wěn)態(tài)的耦合損耗因子表征子系統(tǒng)間瞬態(tài)能量流動的不確定性。H.B.Sun等[35]和J.F.Fahy等[36]分別利用修改的統(tǒng)計能量分析預(yù)測隨時間變化的振動,這是統(tǒng)計能量分析在瞬態(tài)領(lǐng)域的應(yīng)用。

    G.Borello[37]提出用子結(jié)構(gòu)模態(tài)信息,將統(tǒng)計能量分析丟失的相位信息補上去,獲得相位重構(gòu)的傳遞函數(shù),再通過傅里葉反變換得到隨時間變化的結(jié)構(gòu)響應(yīng),這種方法稱為局部模態(tài)相位信息重構(gòu)法(LMPR)?;诮y(tǒng)計能量分析的LMPR方法,在歐洲獲得了大量應(yīng)用[38],如SPOT5衛(wèi)星、Smart1衛(wèi)星和Ariane5火箭的星箭分離沖擊力學(xué)環(huán)境分析,該方法被成功集成進商用軟件SEA+中[39]。

    從原理上,統(tǒng)計能量分析是解決高頻振動問題的方法,當結(jié)構(gòu)模態(tài)密度較小(小于2或3)時,該方法將失效[7]。并且,無論是穩(wěn)態(tài)或者瞬態(tài)SEA分析,都是線性分析,理論上僅適合估計距離爆炸裝置較遠位置處的高頻響應(yīng)。另外,由于采用了統(tǒng)計平均的概念,統(tǒng)計能量方法僅能得到在一個區(qū)域上某個頻段內(nèi)的響應(yīng)包絡(luò),無法得到具體位置上的響應(yīng)。

    2.2.4 虛模態(tài)綜合法

    虛模態(tài)綜合法全稱是虛擬模態(tài)綜合與仿真(virtual mode synthesis and simulation,VMSS)。該方法將統(tǒng)計能量分析(SEA)方法的原理引入到經(jīng)典模態(tài)分析格式中,先假定已知結(jié)構(gòu)一個穩(wěn)態(tài)頻率響應(yīng)包絡(luò),如正弦激勵作用下結(jié)構(gòu)上某一點響應(yīng)峰值。初始假定的頻率響應(yīng)既可以是SEA的穩(wěn)態(tài)解,也可以是有限元的結(jié)果,還可以直接通過試驗測量。在高頻段,根據(jù)子結(jié)構(gòu)模態(tài)密度劃分頻段,該頻段上局部模態(tài)響應(yīng)的峰值即可作為每一頻段上的響應(yīng)包絡(luò);虛擬模態(tài)綜合方法依然采用經(jīng)典模態(tài)分析的格式,因此在低頻段可以直接采用實測或者模態(tài)分析結(jié)果,保證較準確的預(yù)示整個頻段上的響應(yīng),更詳細的解釋可參見文獻[1]。

    虛模態(tài)綜合法由E.C.Dalton等[40]提出,闡明了VMSS的原理和基本公式,并將VMSS法集成為通用計算程序MANTA,應(yīng)用于美軍方裝備設(shè)計中高速撞擊等瞬態(tài)問題的求解[41-43]。在航天領(lǐng)域,R.Ullio等[44]在航天器高速沖擊影響分析(HVI)項目中利用Auto SEA的沖擊分析模塊(VMSS方法)進行了GOCE衛(wèi)星沖擊環(huán)境的預(yù)示研究。D.Lee等[45]采用VMSS方法研究了某低軌地球觀測衛(wèi)星與運載火箭分離沖擊響應(yīng)預(yù)示。在我國,王軍評等[46]采用虛擬模態(tài)綜合法完成了運載火箭級間分離沖擊響應(yīng)預(yù)示。

    相比于SEA或TSEA方法,VMSS方法具有無需假設(shè)激勵為穩(wěn)態(tài)或準穩(wěn)態(tài)、可以獲得時域數(shù)據(jù)、當近場結(jié)構(gòu)具有一定的線性特性時可以估計近場響應(yīng)等優(yōu)點,并且 VMSS方法從原理上能夠?qū)崿F(xiàn)包含低頻、中頻、高頻的寬頻帶的計算。因此,應(yīng)用虛模態(tài)綜合法進行航天器火工沖擊的響應(yīng)預(yù)示,從理論具有合理性,在工程上具有一定的研究價值。

    3 總 結(jié)

    綜上所述,在美國和歐洲,對航天器火工沖擊力學(xué)環(huán)境的研究開展較早,發(fā)展較完備,已經(jīng)掌握了航天器火工沖擊的傳遞機理和環(huán)境控制方法,擁有較完備的模擬試驗方法及試驗技術(shù)標準等。切實可行的環(huán)境預(yù)示方法,仍以基于試驗數(shù)據(jù)的外推方法為主,分析和數(shù)值方法則能夠輔助傳遞機理的分析,對環(huán)境控制提供支持。

    在我國,航天部門對于火工沖擊的認識也較早[47],但是火工品爆炸沖擊問題一直沒有引起足夠的關(guān)注,使我國對火工沖擊力學(xué)環(huán)境問題的研究大大滯后于發(fā)達國家。近年來,由于航天器上攜帶火工品和沖擊敏感儀器設(shè)備數(shù)量的大幅增加,火工沖擊問題已經(jīng)引起航天部門的廣泛關(guān)注,相應(yīng)的地面試驗設(shè)備[18,48-49]、測量技術(shù)[50]、環(huán)境防護方法[51]及仿真計算[52-54]的研究均在積極進行中。2009年,頒布了我國現(xiàn)行的爆炸分離沖擊試驗標準[10]。在沖擊試驗和仿真計算方面的研究取得了很大的進步,但仍存在較大的差距,今后應(yīng)重點開展以下方面的研究。

    (1)開展火工爆炸式的地面模擬試驗設(shè)備的研究。目前,系統(tǒng)級或近場地面模擬試驗主要采用火工爆炸式試驗設(shè)備或采用真實火工品進行;設(shè)備級或中、遠場試驗可以采用非火工爆炸式模擬試驗設(shè)備。在我國,目前還沒有能夠進行系統(tǒng)級地面模擬火工沖擊的試驗設(shè)備,非火工爆炸式的試驗設(shè)備的應(yīng)用范圍有限。因此,應(yīng)著力開展火工爆炸式的地面模擬試驗設(shè)備的研究,避免由于地面試驗引入的過試驗或欠試驗的問題。

    (2)建立火工沖擊響應(yīng)數(shù)據(jù)庫。目前,美國NASA和歐洲宇航部門已經(jīng)對大量的試驗數(shù)據(jù)進行了統(tǒng)計學(xué)分析,初步建立了基于試驗數(shù)據(jù)的沖擊力學(xué)環(huán)境預(yù)示方法。在我國,航天工程也已經(jīng)積累了大量的試驗數(shù)據(jù),迫切需要在此基礎(chǔ)上開展相關(guān)的統(tǒng)計分析工作,建立火工沖擊響應(yīng)數(shù)據(jù)庫,為后續(xù)航天器火工沖擊環(huán)境預(yù)示和試驗條件的制定,提供有力的技術(shù)支撐。

    (3)推進航天器火工沖擊環(huán)境預(yù)示方法的研究。對火工沖擊力學(xué)環(huán)境進行準確的預(yù)示,是對試驗數(shù)據(jù)的有力補充,能夠幫助揭示火工沖擊在復(fù)雜結(jié)構(gòu)中的傳遞機理,幫助建立合理的沖擊試驗條件,并且為航天器火工沖擊力學(xué)環(huán)境防護設(shè)計、建立基于火工沖擊力學(xué)環(huán)境防護原則的航天器總體設(shè)計提供理論參考。

    [1] Dynamic environmental criteria: NASA-HDBK-7005[S]. USA, 2001.

    [2] Pyroshock test criteria: NASA-STD-7003A[S]. USA, 2011.

    [3] Environment engineering considerations and laboratory tests: MIL-STD-810F[S]. USA: Department of Defense, 2000.

    [4] Filippi E. Pyroshock simulation using the alcatel etca test facility[C]∥Launch Vehicle Vibrations. First European Conference. Toulouse: CNES, 1999.

    [5] Meoning C J. Pyrotechnic shock flight failures[C]∥Institute of Environmental Sciences Pyrotechnic Shock Tutorial Program, 31st Annual Technical Meeting. 1985.

    [6] Mcgrath M B, Rader W P. Aerospace systems pyrotechnic shock data. Volume Ⅶ: Investigation of mass loading effects[R]. N71-19250, 1970.

    [7] Environment engineering considerations and laboratory tests[M]. USA: Department of Defense, 2000.

    [8] 丁繼鋒,趙欣,韓增堯,等.航天器火工沖擊技術(shù)研究進展[J].宇航學(xué)報,2014,35(12):1339-1349. Ding Jifeng, Zhao Xin, Han Zengyao. Research development of spacecraft pyroshock technique[J]. Journal of Astronautics, 2014,35(12):1339-1349.

    [9] Piersol A G, Paez T L. Harris’ shock and vibration handbook[M]. 6th ed. New York: The McGraw-Hill Companies, Inc, 2010.

    [10] 軍用裝備實驗室環(huán)境試驗方法 第27部分:爆炸分離沖擊實驗:GJB 150.27-2009[S]. 2009.

    [11] Lee J R, Chia C C, Kong C W. Review of pyroshock wave measurement and simulation for space systems[J]. Measurement, 2012,45(4):631-642.

    [12] Http://Telecom.Esa.Int/Telecom/www/Object/Index.Cfm?Fobjectid=30565, Available.

    [13] Girard A, Pont C, Courau E. Very first results form pyroshock database for satellites[C]∥European Conference on Spacecaft Structures, Materials & Mechenical Testing 2005. Noordwijk, The Netherlands, 2005.

    [14] Courau E, Roy P. Ariane 5 shock environment qualification for an earth observation satellite[C]∥European Conference on Spacecraft Structure, Material & Mechanical Testing 2005. Noordwijk, The Netherlands, 2005.

    [15] Filippi E, Dolez F. Development of the alcatel etca pyroshock test facility[C]∥European Conference on Spacecraft Structures, Materials and Mechanical Testing 1999. Braunschweig, Germany, 1999.

    [16] Filippi E, Attouoman H, Dolez B. Pyroshcok simulation using the alcatel etca test facility[C]∥Launch Vehicle Vibrations. First European Conference. Toulouse: CNES, 1999.

    [17] Dilhan D, Piquereau A, Bonnes L. Definition and manufacturing of the pyroshock bench[C]∥7th ESA/CNES International Workshop on Space Pyrotechnics ESTEC. 2008.

    [18] 劉斌,吳江.用火工品產(chǎn)生中高量級沖擊的試驗技術(shù)研究[J].強度與環(huán)境,2007,34(3):8-13. Liu Bin, Wu Jiang. A study of producing high level shock environment by pyrotechnic device[J]. Structure & Environment Engineering, 2007,34(3):8-13.

    [19] Dilhan D, Cipolla V, Grzeskowiak H. Pyroshock generation[C]∥European Conference on Spacecaft Structures, Materials & Mechenical Testing 2005. Noordwijk, The Netherlands, 2005.

    [20] Ali K, Juan F. Tunable beam pyroshock simulation system[C]∥Spacecraft and Launch Vehicle Dynamic Environments Workshop. California, USA, 2007.

    [21] Benedetti M D, Garofalo G, Zumpano M, et al. On the damping effect due to bolted junctions in space structures subjected to pyro-shock[J]. Acta Astronautica, 2007,60(12):947-956.

    [22] Himelblau H, Piersol A G, Wise J H, et al. Handbook for dynamic data acquisition and analysis. Appendix A: Pyroshock data acquisition and analysis[M]. Mt Prospect, IL: Institute of Environmental Sciences, 1994.

    [23] Evans M J, Neubert V H, Bement L J. Measurement, data analysis, and prediction of oyrotechnic shock from pin-pullers and separation joints[J]. Shock and Vibration Bulletin, 1987,57(2).

    [24] Smallwood D O. An improved recursive formula for calculating shock response spectra[J]. Shock and Vibration Bulletin, 1981,51(2):211-217.

    [25] Hughes W O, Mcnelis A M. Statistical analysis of a large sample size pyroshock test data set including post flight data assessment[R]. Cleveland, Ohio: Glenn Research Center, 2009.

    [26] Shi Q, Ando S, Seko H, et al. The summarization of pyroshock testing data and SRS level prediction methodology[C]∥The 5th International Symposium on Environmental Testing for Space Programmes. Noordwijk, The Netherlands, 2004.

    [27] Ullio R, Marucchi-Chierro P C. Utlization of prediction methods in the shock environment evaluation[C]∥The European Conference on Spacecraft Structure, Materials and Mechanical Testing 2001. Noofdwijk, The Netherlands, 2001.

    [28] Girard A, Courau E, Bugeat L P. Pyroshock database for satellites[C]∥European Conferenc on Spacecraft Structures, Materials and Mechanical Testing 2000. Noordwijk, the Netherlands, 2000.

    [29] Zukas J A. Introduction to hydrocodes[M]. Elsevier, 2004.

    [30] 白金澤.LS-DYNA 3D理論基礎(chǔ)與實例分析[M].北京:科學(xué)出版社,2005.

    [31] 張雄,王天舒.計算動力學(xué)[M].北京:清華大學(xué)出版社,2007.

    [32] Mary S, Cipolla V, Courau E, et al. Shock propagation simulation using FEM software[C]∥European Conference on Spacecaft Structures, Materials & Mechenical Testing 2005. Noordwijk, The Netherlands, 2005.

    [33] 姚德源,王其政.統(tǒng)計能量分析原理及其應(yīng)用[M].北京:北京理工大學(xué)出版社,1995.

    [34] Manning J E, Lee K. Predicting mechanical shock transmission[J]. Shock and Vibration Bulletin, 1968,37(4):65-70.

    [35] Sun H B, Sun J C, Richards E J. Prediction of total loss factors of structures. Part Ⅲ: Effective loss factors in quasi-transient conditions[J]. Journal of Sound and Vibration, 1986,106(3):465-479.

    [36] Fahy J F, Yao D Y. Power flow between non-conservatively coupled oscillators[J]. Journal of Sound and Vibration, 1987,114(1):1-11.

    [37] Borello G. SEA航天應(yīng)用[M]. Belgium: Inter AC, 2010.

    [38] Borello G, Courjal A. Modelling payloads using SEA for vibroacoustic and shock prediction[R]. France: Inter AC-L’Union, 2005.

    [39] SEA: Shock module of SEA+ user’s guide[M]. Inter AC, Inc, 2013.

    [40] Dalton E C, Chambers B S. Analysis and validation testing of impulsive load response in complex, multi-compartmented structures[C]∥Proceedings of the 36th AIAA Structures, Structural Dynamics, and Materials Conference. 1995.

    [41] Dalton E C, Frydman A, Li A, et al. High frequency shock predictions in armored vehicles: Ags case study[C]∥Proceedings of the 16th International Symposium on Ballistics. San Francisco, CA, 1996.

    [42] Dalton E C, Loper R B, Frydman A, et al. Simulation of ballistic shock in composite armored vehicles[C]∥Proceedings of the 68th Shock and Vibration Symposium. 1997.

    [43] Dalton E C. High frequency shock prediction, short course notes[M]. Military Technology, Inc, 1999.

    [44] Ullio R, Marucchi-Chierro P C. Auto SEA shock application on shock event simulation: Study case and problematics encountered[R]. Toulouse: Euro PAM, 2006.

    [45] Lee D O, Han J H, Jang H W, et al. Shock respones prediction of a low altitude earth observation satellite during launch vechicle separation[J]. International Journal of Aeronautical & Space Sciences, 2010,11(1):49-57.

    [46] 王軍評,毛勇建,黃含軍,等.統(tǒng)計能量分析法在爆炸分離沖擊響應(yīng)預(yù)示中的應(yīng)用[J].航天器環(huán)境工程,2011,28(5):414-420. Wang Junping, Mao Yongjian, Huang Hanjun, et al. Application of statistical energy analysis method in prediction of pyroshock responses[J]. Spacecraft Environment Engineering, 2011,28(5):414-420.

    [47] 柯受全,金恂叔.衛(wèi)星環(huán)境工程和模擬試驗(下)[M].北京:宇航出版社,1996.

    [48] 焦亮,鄧明.爆炸分離沖擊試驗分析[J].電子產(chǎn)品可靠性與環(huán)境試驗,2013,31(1):31-34. Jiao Liang, Deng Ming. The pyroshock test[J]. Electronic Product Reliability and Environmental Testing, 2013,31(1):31-34.

    [49] 馬斌捷,張建華,吳江.火工品爆炸加載方法在火箭分離沖擊環(huán)境模擬試驗中的應(yīng)用與效果[J].強度與環(huán)境,2007,34(5):1-7. Ma Binjie, Zhang Jianhua, Wu Jiang. Applications and effects of pyrotechnic explode loading in shock environment simulation experiment of rocket separation[J]. Structure & Enviroment Engineering, 2007,34(5):1-7.

    [50] 畢文輝,嚴楠, 何春全,等.航天火工品爆炸沖擊多參量的測試[J].計測技術(shù),2009,29(3):14-16. Bi Wenhui, Yan Nan, He Chunquan, et al. Measurement of pyroshock for explosive initiating device used in spacecraft[J]. Metrology & Measurement Technology, 2009,29(3):14-16.

    [51] 張歡,劉天雄,李長江,等.航天器火工沖擊環(huán)境防護技術(shù)現(xiàn)狀與應(yīng)用[J].航天器工程,2014,23(2):104-113. Zhang Huan, Liu Tianxiong, Li Changjiang, et al. Status and application analysis of spacecraft pyroshock protection techniques[J]. Spacecraft Engineering, 2014,23(2):104-113.

    [52] 毛勇建,李玉龍,陳穎,等.炸藥條加載圓柱殼的數(shù)值模擬(Ⅱ):解耦分析與實驗驗證[J].高壓物理學(xué)報,2013,27(1):76-82. Mao Yongjian, Li Yulong, Chen Ying, et al. Numerical simulation of cylindrical shell loaded by explosive rods (Ⅱ): Decoupling analysis and experimental validation[J]. Chinese Journal of High Pressure Physics, 2013,27(1):76-82.

    [53] 毛勇建,李玉龍,陳穎,等.炸藥條加載圓柱殼的數(shù)值模擬(Ⅰ):流固耦合模擬[J].高壓物理學(xué)報,2012,26(2):155-162. Mao Yongjian, Li Yulong, Chen Ying, et al. Numerical simulation of cylindrical shell loaded by explosive rods (Ⅰ): Fluid-structure interaction simulation[J]. Chinese Journal of High Pressure Physics, 2012,26(2):155-162.

    [54] 王軍評,毛勇建,黃含軍.點式火工分離裝置沖擊載荷作用機制的數(shù)值模擬研究[J].振動與沖擊,2013,32(2):9-13. Wang Junping, Mao Yongjian, Huang Hanjun. Numerical simulation for impulsively loading mechanism of a point pyrotechnic separation device[J]. Journal of Vibration and Shock, 2013,32(2):9-13.

    (責(zé)任編輯 丁 峰)

    Review of pyroshock simulation and response prediction methods in spacecraft

    Zhao Xin1, Ding Jifeng1, Han Zengyao2, Zou Yuanjie1

    (1.InstituteofSpacecraftSystemEngineering,ChinaAcademyofSpaceTechnology,
    Beijing100094,China;2.InstituteofResearchandDevelopment,ChinaAcademyofSpaceTechnology,
    Beijing100094,China)

    The pyroshock environment of satellite-rocket separation is the severest mechanical environment during launching, which is characterized by transient high acceleration and high frequency. While it does not necessarily cause a satellite any structural damage, pyroshock may incur most serious damages on a satellite's precision electronic equipments containing crystals and brittle materials, resulting in either the failure of the entire mission or even catastrophic accidents. Therefore, during the development of a new spacecraft, an accurate prediction of the pyroshock environment and a reasonable specification and determination for components are essential. In this paper, a research review of ground simulation test methods and the pyroshock response prediction is presented, and the technological gap between China and countries highly developed in field is pointed out. In addition to that, according to the requirements of China’s domestic space engineering, the main research directions in the pyroshock are proposed.

    oscillation theory; response prediction; pyroshock; spacecraft; ground simulation test; hydrocodes; explicit finite element analysis; statistical energy analysis; virtual mode synthesis and simulation

    10.11883/1001-1455(2016)02-0259-10

    2014-07-27;

    趙 欣(1982— ),女,博士,工程師,zhaoxinbhu@126.com。

    O32;V414;V416 國標學(xué)科代碼: 1302099

    A

    修回日期: 2014-08-21

    猜你喜歡
    火工火工品航天器
    2022 年第二季度航天器發(fā)射統(tǒng)計
    國際太空(2022年7期)2022-08-16 09:52:50
    關(guān)于火工藥劑發(fā)展的幾點思考
    含能材料(2022年4期)2022-04-16 06:28:04
    關(guān)于“火工煙火技術(shù)‘十四·五’發(fā)展論壇”征文的通知
    火工品(2020年1期)2020-11-28 13:56:51
    電火工品儲存安全評估模型與應(yīng)用研究
    化工管理(2020年16期)2020-10-03 11:26:48
    2019 年第二季度航天器發(fā)射統(tǒng)計
    國際太空(2019年9期)2019-10-23 01:55:34
    2018 年第三季度航天器發(fā)射統(tǒng)計
    國際太空(2018年12期)2019-01-28 12:53:20
    2018年第二季度航天器發(fā)射統(tǒng)計
    國際太空(2018年9期)2018-10-18 08:51:32
    航天器火工引爆鏈路系統(tǒng)設(shè)計驗證方法研究實踐
    淺析火工品安全運輸存儲要求
    化工管理(2017年10期)2017-03-04 09:26:25
    更正
    含能材料(2015年6期)2015-03-27 08:51:45
    少妇人妻久久综合中文| 亚洲激情五月婷婷啪啪| 高清欧美精品videossex| 亚洲国产欧美在线一区| 欧美亚洲 丝袜 人妻 在线| 国产成人freesex在线| 亚洲国产精品国产精品| 2021少妇久久久久久久久久久| 精品久久国产蜜桃| 欧美日韩精品成人综合77777| 黄色日韩在线| av在线播放精品| 亚洲av成人精品一二三区| 国产永久视频网站| 大片电影免费在线观看免费| 熟女电影av网| 欧美老熟妇乱子伦牲交| 99久久人妻综合| 亚洲四区av| 久久久久久久精品精品| 久久久久视频综合| 亚洲精品日韩av片在线观看| 国产av精品麻豆| 91精品一卡2卡3卡4卡| 中文字幕人妻熟人妻熟丝袜美| 亚洲欧美日韩另类电影网站| 91精品伊人久久大香线蕉| 一级毛片我不卡| 简卡轻食公司| 热99国产精品久久久久久7| 在线精品无人区一区二区三| 夫妻性生交免费视频一级片| 人人妻人人看人人澡| 三级经典国产精品| 日韩一本色道免费dvd| 亚洲美女黄色视频免费看| 欧美三级亚洲精品| 久久久亚洲精品成人影院| 十分钟在线观看高清视频www | av天堂久久9| 久久99热这里只频精品6学生| 午夜老司机福利剧场| 国精品久久久久久国模美| 亚洲欧美日韩东京热| 大又大粗又爽又黄少妇毛片口| 欧美丝袜亚洲另类| 成年人午夜在线观看视频| 黄色毛片三级朝国网站 | 特大巨黑吊av在线直播| 下体分泌物呈黄色| 亚洲婷婷狠狠爱综合网| 日韩成人av中文字幕在线观看| 日本猛色少妇xxxxx猛交久久| 精品一区二区免费观看| 亚洲美女黄色视频免费看| 22中文网久久字幕| av天堂中文字幕网| 国产探花极品一区二区| 亚洲丝袜综合中文字幕| 久久久久久久久久久丰满| 免费看av在线观看网站| 日韩伦理黄色片| 亚洲精品视频女| 国产成人精品久久久久久| 性高湖久久久久久久久免费观看| 尾随美女入室| av视频免费观看在线观看| 免费少妇av软件| 久久久久久久久久久丰满| 久久久久久久大尺度免费视频| 在线免费观看不下载黄p国产| 美女cb高潮喷水在线观看| 国产av国产精品国产| 欧美性感艳星| 成人特级av手机在线观看| 国产亚洲91精品色在线| 国产欧美日韩综合在线一区二区 | 欧美国产精品一级二级三级 | 伊人久久国产一区二区| 久久人人爽人人片av| 99久久精品国产国产毛片| 国产免费视频播放在线视频| 18禁在线播放成人免费| 黑人猛操日本美女一级片| 啦啦啦在线观看免费高清www| 桃花免费在线播放| 又大又黄又爽视频免费| 日本午夜av视频| 欧美精品一区二区免费开放| 国产精品国产三级专区第一集| 天天躁夜夜躁狠狠久久av| 三级国产精品欧美在线观看| 国产有黄有色有爽视频| 久久久久国产网址| 精品人妻熟女毛片av久久网站| 亚洲av成人精品一区久久| 又爽又黄a免费视频| 欧美成人午夜免费资源| 如何舔出高潮| 日韩不卡一区二区三区视频在线| 国产极品天堂在线| 日韩人妻高清精品专区| 欧美精品一区二区大全| 国产熟女欧美一区二区| 久久久久久久久久人人人人人人| 亚洲国产日韩一区二区| 色5月婷婷丁香| 午夜免费男女啪啪视频观看| 久久这里有精品视频免费| 亚洲丝袜综合中文字幕| 99热这里只有是精品50| 亚洲人成网站在线播| 特大巨黑吊av在线直播| 国产一区亚洲一区在线观看| 久久99精品国语久久久| 黄片无遮挡物在线观看| 久久99蜜桃精品久久| 日产精品乱码卡一卡2卡三| 亚洲精品自拍成人| 日韩av在线免费看完整版不卡| 男女免费视频国产| 精品熟女少妇av免费看| 夫妻性生交免费视频一级片| 女人久久www免费人成看片| 国产成人一区二区在线| 国产成人午夜福利电影在线观看| 男女无遮挡免费网站观看| 日韩电影二区| 精品一品国产午夜福利视频| 午夜免费观看性视频| 久久久精品免费免费高清| h视频一区二区三区| 一级毛片aaaaaa免费看小| 少妇被粗大的猛进出69影院 | 国产精品福利在线免费观看| 街头女战士在线观看网站| 亚洲精品日韩av片在线观看| 亚洲熟女精品中文字幕| 日本色播在线视频| 在线观看免费日韩欧美大片 | 国内少妇人妻偷人精品xxx网站| 日韩欧美 国产精品| 国产女主播在线喷水免费视频网站| 久久热精品热| 男人舔奶头视频| 中文欧美无线码| 亚洲精品日韩av片在线观看| 人人妻人人澡人人看| 久久久久久久久大av| 国产精品国产三级国产av玫瑰| 自拍欧美九色日韩亚洲蝌蚪91 | 国产免费又黄又爽又色| h视频一区二区三区| 尾随美女入室| 五月开心婷婷网| 亚洲电影在线观看av| 成人毛片60女人毛片免费| 一级,二级,三级黄色视频| 免费人妻精品一区二区三区视频| 一区二区三区精品91| 精品国产一区二区三区久久久樱花| 亚洲av二区三区四区| 日本黄大片高清| 亚洲精品乱码久久久久久按摩| 日本与韩国留学比较| 国内揄拍国产精品人妻在线| 亚洲av.av天堂| 精品久久国产蜜桃| 国产精品女同一区二区软件| 免费观看无遮挡的男女| 国语对白做爰xxxⅹ性视频网站| 国产精品嫩草影院av在线观看| 99久久精品国产国产毛片| 久久久久网色| 51国产日韩欧美| 久久精品久久久久久久性| 丝瓜视频免费看黄片| 夜夜爽夜夜爽视频| 亚洲欧美成人精品一区二区| 不卡视频在线观看欧美| 色94色欧美一区二区| 国产成人精品无人区| 午夜av观看不卡| 少妇人妻 视频| 亚洲精品色激情综合| 美女xxoo啪啪120秒动态图| 黄色毛片三级朝国网站 | av专区在线播放| 内射极品少妇av片p| 色94色欧美一区二区| 伊人亚洲综合成人网| 国产 一区精品| 色吧在线观看| 男人舔奶头视频| 免费播放大片免费观看视频在线观看| av又黄又爽大尺度在线免费看| 我要看日韩黄色一级片| 中文字幕av电影在线播放| 夜夜骑夜夜射夜夜干| 国产日韩欧美视频二区| 免费黄网站久久成人精品| 又黄又爽又刺激的免费视频.| 另类亚洲欧美激情| 少妇高潮的动态图| 日本黄大片高清| 久久国产乱子免费精品| 极品人妻少妇av视频| 国国产精品蜜臀av免费| 国精品久久久久久国模美| 国产老妇伦熟女老妇高清| 香蕉精品网在线| 美女中出高潮动态图| 亚洲欧美日韩东京热| 国产女主播在线喷水免费视频网站| 欧美日韩视频精品一区| 国产成人91sexporn| 成人漫画全彩无遮挡| 性色avwww在线观看| 久久人人爽人人片av| av黄色大香蕉| 国产成人午夜福利电影在线观看| 老司机影院毛片| 午夜福利网站1000一区二区三区| av播播在线观看一区| 国产一区亚洲一区在线观看| 嫩草影院新地址| 亚洲精品国产av蜜桃| 午夜福利,免费看| 日韩av在线免费看完整版不卡| 精品久久久久久久久av| 日韩欧美精品免费久久| 秋霞在线观看毛片| 日韩电影二区| 国产有黄有色有爽视频| 亚洲熟女精品中文字幕| 亚洲欧洲日产国产| 毛片一级片免费看久久久久| 精品亚洲乱码少妇综合久久| 交换朋友夫妻互换小说| 嫩草影院新地址| 国产一区有黄有色的免费视频| 婷婷色麻豆天堂久久| 久久精品久久精品一区二区三区| 精品亚洲成a人片在线观看| 纵有疾风起免费观看全集完整版| 日韩av免费高清视频| 国产日韩欧美视频二区| 内射极品少妇av片p| 亚洲欧美一区二区三区黑人 | 国产极品天堂在线| 亚洲综合色惰| 久久久久人妻精品一区果冻| 久久午夜福利片| 国产日韩欧美在线精品| 亚洲无线观看免费| 欧美xxⅹ黑人| 国产一区二区在线观看日韩| 午夜久久久在线观看| 国产免费一区二区三区四区乱码| 人妻系列 视频| 国产成人午夜福利电影在线观看| 午夜福利在线观看免费完整高清在| 亚洲欧洲国产日韩| 青春草亚洲视频在线观看| 精品久久久精品久久久| 欧美3d第一页| 日韩视频在线欧美| 国产欧美亚洲国产| 成人毛片a级毛片在线播放| 少妇人妻一区二区三区视频| 国产精品人妻久久久影院| 久久人人爽av亚洲精品天堂| 日本欧美国产在线视频| 秋霞伦理黄片| 亚洲欧美精品自产自拍| 免费人妻精品一区二区三区视频| 国产精品免费大片| 人妻少妇偷人精品九色| 亚洲欧美成人精品一区二区| 日韩av不卡免费在线播放| 亚洲国产精品一区三区| 欧美3d第一页| 欧美日本中文国产一区发布| 哪个播放器可以免费观看大片| 中文字幕av电影在线播放| 十八禁高潮呻吟视频 | 美女内射精品一级片tv| av在线app专区| 大话2 男鬼变身卡| 多毛熟女@视频| 国产精品伦人一区二区| 99精国产麻豆久久婷婷| 九色成人免费人妻av| 国产亚洲一区二区精品| 黄色配什么色好看| 极品人妻少妇av视频| 一本大道久久a久久精品| 赤兔流量卡办理| 看非洲黑人一级黄片| 国产精品嫩草影院av在线观看| 久久av网站| 日韩精品有码人妻一区| 一区二区三区免费毛片| 亚洲国产毛片av蜜桃av| 亚洲一区二区三区欧美精品| 久久久久久久精品精品| 国产亚洲av片在线观看秒播厂| 国产熟女欧美一区二区| 亚洲自偷自拍三级| 少妇人妻精品综合一区二区| 国产在线视频一区二区| 久久久久久久久大av| 最近的中文字幕免费完整| 国产探花极品一区二区| 草草在线视频免费看| 插阴视频在线观看视频| 两个人免费观看高清视频 | 亚洲人与动物交配视频| 国产永久视频网站| 偷拍熟女少妇极品色| 国产91av在线免费观看| 欧美日韩综合久久久久久| 大片电影免费在线观看免费| 欧美xxxx性猛交bbbb| 久久精品国产亚洲av涩爱| 亚洲美女视频黄频| 美女cb高潮喷水在线观看| 在线观看免费视频网站a站| 自拍偷自拍亚洲精品老妇| 久久国产亚洲av麻豆专区| 极品少妇高潮喷水抽搐| 久久99一区二区三区| 青春草亚洲视频在线观看| 久久久久久久久久久丰满| 美女大奶头黄色视频| 校园人妻丝袜中文字幕| 中文字幕免费在线视频6| 三级国产精品欧美在线观看| 欧美性感艳星| 高清在线视频一区二区三区| 日韩欧美 国产精品| 国产黄片视频在线免费观看| 国产精品麻豆人妻色哟哟久久| 欧美精品人与动牲交sv欧美| 国产精品女同一区二区软件| 久久精品国产亚洲av天美| 欧美三级亚洲精品| 国产高清国产精品国产三级| 久久久久久久国产电影| 老司机影院毛片| 夫妻午夜视频| 中文字幕人妻熟人妻熟丝袜美| 2022亚洲国产成人精品| 亚洲欧美日韩卡通动漫| 亚洲av在线观看美女高潮| 亚洲av二区三区四区| 国产伦精品一区二区三区视频9| 热re99久久国产66热| 99久久人妻综合| 欧美最新免费一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 六月丁香七月| 自线自在国产av| 精品人妻熟女毛片av久久网站| 边亲边吃奶的免费视频| 一级毛片 在线播放| 久久久久久久精品精品| 久久99热6这里只有精品| 精品一品国产午夜福利视频| 蜜桃在线观看..| 亚洲av日韩在线播放| av天堂久久9| 免费在线观看成人毛片| 少妇的逼水好多| 在线观看av片永久免费下载| 日韩av免费高清视频| 久久久久网色| 国产免费又黄又爽又色| 精品国产露脸久久av麻豆| 国产成人91sexporn| 久久国产亚洲av麻豆专区| 亚洲精品国产av蜜桃| 国产精品.久久久| 这个男人来自地球电影免费观看 | 国产老妇伦熟女老妇高清| 三上悠亚av全集在线观看 | 国产一区亚洲一区在线观看| 午夜91福利影院| 老熟女久久久| 99re6热这里在线精品视频| 国产精品成人在线| 国产高清国产精品国产三级| 亚洲欧美一区二区三区国产| 亚洲三级黄色毛片| 国产成人精品久久久久久| 亚洲av福利一区| 在线 av 中文字幕| 国产免费又黄又爽又色| 成人毛片60女人毛片免费| 久久久国产一区二区| 免费观看性生交大片5| 大香蕉97超碰在线| 中文字幕亚洲精品专区| 狠狠精品人妻久久久久久综合| 国产成人freesex在线| 欧美日韩视频精品一区| 亚洲,一卡二卡三卡| 国产91av在线免费观看| 波野结衣二区三区在线| 亚洲欧美成人综合另类久久久| 一二三四中文在线观看免费高清| 国产精品伦人一区二区| 国产探花极品一区二区| 日韩中字成人| 久久韩国三级中文字幕| 久久精品国产a三级三级三级| 熟女电影av网| 男女边吃奶边做爰视频| 亚洲av不卡在线观看| 久久久久久久久久久免费av| 亚洲精品乱久久久久久| 欧美精品亚洲一区二区| 欧美最新免费一区二区三区| 久久99热这里只频精品6学生| 又爽又黄a免费视频| 亚洲激情五月婷婷啪啪| 啦啦啦视频在线资源免费观看| 99热这里只有是精品50| 蜜臀久久99精品久久宅男| 国产精品久久久久久av不卡| 免费人成在线观看视频色| 国产一区二区三区av在线| 能在线免费看毛片的网站| 欧美变态另类bdsm刘玥| 久久精品久久久久久噜噜老黄| 色婷婷久久久亚洲欧美| 精品少妇黑人巨大在线播放| 各种免费的搞黄视频| 欧美丝袜亚洲另类| 夜夜看夜夜爽夜夜摸| 久久午夜综合久久蜜桃| av黄色大香蕉| 国产精品国产三级专区第一集| 王馨瑶露胸无遮挡在线观看| 久热这里只有精品99| 日韩精品有码人妻一区| 水蜜桃什么品种好| 特大巨黑吊av在线直播| 国产成人精品久久久久久| 少妇裸体淫交视频免费看高清| 国产精品麻豆人妻色哟哟久久| 街头女战士在线观看网站| 高清欧美精品videossex| 亚洲精品色激情综合| 99热6这里只有精品| 精品国产一区二区三区久久久樱花| 搡老乐熟女国产| 久久久国产欧美日韩av| 国产精品无大码| 亚洲av免费高清在线观看| 日本黄色日本黄色录像| 又粗又硬又长又爽又黄的视频| 婷婷色av中文字幕| 午夜福利网站1000一区二区三区| 日日摸夜夜添夜夜爱| 午夜视频国产福利| 日日撸夜夜添| 欧美3d第一页| 我的老师免费观看完整版| 少妇猛男粗大的猛烈进出视频| 赤兔流量卡办理| 日韩大片免费观看网站| 日本黄色片子视频| 三级国产精品欧美在线观看| 91精品国产国语对白视频| 80岁老熟妇乱子伦牲交| 菩萨蛮人人尽说江南好唐韦庄| 国产视频内射| 三级国产精品欧美在线观看| 久久国产亚洲av麻豆专区| kizo精华| 又粗又硬又长又爽又黄的视频| 自线自在国产av| 国产精品成人在线| 亚洲国产精品成人久久小说| 国产白丝娇喘喷水9色精品| 精品一区在线观看国产| 亚洲国产毛片av蜜桃av| 国产欧美日韩综合在线一区二区 | 精品一区二区三卡| 精品一区二区免费观看| 亚洲精品中文字幕在线视频 | 国产男女内射视频| 亚洲av在线观看美女高潮| 国产一区二区三区av在线| 国产 精品1| 欧美日韩一区二区视频在线观看视频在线| 久久99热这里只频精品6学生| 人妻 亚洲 视频| 日韩三级伦理在线观看| 亚洲图色成人| av免费观看日本| 青春草国产在线视频| 成人毛片a级毛片在线播放| 插阴视频在线观看视频| 国产成人精品福利久久| 伦理电影免费视频| 大又大粗又爽又黄少妇毛片口| 国产精品秋霞免费鲁丝片| 亚洲精品中文字幕在线视频 | 一级毛片电影观看| 国产日韩欧美在线精品| 极品人妻少妇av视频| 欧美 日韩 精品 国产| 最新的欧美精品一区二区| 国产淫片久久久久久久久| 黑丝袜美女国产一区| 中文字幕av电影在线播放| 人妻一区二区av| 亚洲精品国产色婷婷电影| 日韩电影二区| 国产精品人妻久久久影院| 又爽又黄a免费视频| 又大又黄又爽视频免费| 亚洲精品中文字幕在线视频 | av天堂久久9| 国产av一区二区精品久久| 国产亚洲欧美精品永久| 久久久a久久爽久久v久久| 日韩制服骚丝袜av| 国产爽快片一区二区三区| 国产熟女午夜一区二区三区 | 免费看av在线观看网站| 汤姆久久久久久久影院中文字幕| 水蜜桃什么品种好| 边亲边吃奶的免费视频| av不卡在线播放| 亚洲电影在线观看av| 午夜老司机福利剧场| 老司机影院毛片| 六月丁香七月| 女性生殖器流出的白浆| 国产av国产精品国产| 十八禁网站网址无遮挡 | 蜜桃在线观看..| 久久久精品94久久精品| 男女无遮挡免费网站观看| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲av男天堂| 高清毛片免费看| 亚洲欧洲国产日韩| 免费黄网站久久成人精品| 欧美bdsm另类| 国产欧美日韩一区二区三区在线 | 最新中文字幕久久久久| 在线亚洲精品国产二区图片欧美 | 日韩精品有码人妻一区| 三上悠亚av全集在线观看 | 亚洲成人一二三区av| 欧美精品一区二区免费开放| 久久精品久久久久久久性| 国产黄频视频在线观看| 韩国高清视频一区二区三区| 自拍偷自拍亚洲精品老妇| 久久99热6这里只有精品| 国产精品福利在线免费观看| 丰满饥渴人妻一区二区三| 亚洲精品久久久久久婷婷小说| 久久人人爽人人片av| 久久99精品国语久久久| 日本黄色片子视频| 国产一区二区三区av在线| 丝袜在线中文字幕| 精品国产乱码久久久久久小说| 亚洲人成网站在线播| 亚洲美女视频黄频| 精品亚洲乱码少妇综合久久| 七月丁香在线播放| 99热这里只有是精品50| 国产成人精品婷婷| 久久久精品免费免费高清| 国产视频首页在线观看| 尾随美女入室| 嫩草影院新地址| 国产白丝娇喘喷水9色精品| videossex国产| 免费观看av网站的网址| 一本一本综合久久| 边亲边吃奶的免费视频| 51国产日韩欧美| 日本与韩国留学比较| 丝袜在线中文字幕| 少妇人妻 视频| 观看免费一级毛片| 国产亚洲91精品色在线| 99久久精品一区二区三区| 久久ye,这里只有精品| 人人妻人人爽人人添夜夜欢视频 | 中文字幕亚洲精品专区| 菩萨蛮人人尽说江南好唐韦庄| 一个人免费看片子| 一级a做视频免费观看| 久久亚洲国产成人精品v| 免费观看性生交大片5| 亚洲av成人精品一区久久| 啦啦啦在线观看免费高清www| 国产伦精品一区二区三区四那| 我要看日韩黄色一级片| 欧美日本中文国产一区发布| 亚洲高清免费不卡视频| 制服丝袜香蕉在线| 在线免费观看不下载黄p国产| 深夜a级毛片| 久久精品国产鲁丝片午夜精品| 大码成人一级视频|