某彈簧分離裝置沖擊環(huán)境統(tǒng)計分析及天地差異性研究

高 慶， 商 霖， 蘇大亮， 朱 辰， 蔡毅鵬， 趙永輝

(1.南京航空航天大學(xué) 機械結(jié)構(gòu)力學(xué)及控制國家重點實驗室,南京 210016;2.中國運載火箭技術(shù)研究院，北京 100076; 3.中國航天科技集團公司，北京 100048)

某彈簧分離裝置沖擊環(huán)境統(tǒng)計分析及天地差異性研究

高 慶1,2， 商 霖2， 蘇大亮3， 朱 辰2， 蔡毅鵬2， 趙永輝1

(1.南京航空航天大學(xué) 機械結(jié)構(gòu)力學(xué)及控制國家重點實驗室,南京 210016;2.中國運載火箭技術(shù)研究院，北京 100076; 3.中國航天科技集團公司，北京 100048)

在分析某彈簧分離裝置沖擊環(huán)境的過程中，發(fā)現(xiàn)設(shè)備支架相同位置處的沖擊數(shù)據(jù)在多次地面和飛行試驗之間存在較大差異。統(tǒng)計分析表明兩組沖擊環(huán)境近似服從對數(shù)正態(tài)分布；一致性檢驗結(jié)果證明兩沖擊環(huán)境存在顯著差異，但散差量級基本一致；初步分析認(rèn)為邊界條件和測量系統(tǒng)是導(dǎo)致天地差異性的主要原因。這些結(jié)論能為沖擊環(huán)境的統(tǒng)計分析、條件制定以及試驗設(shè)計等提供參考。

高頻沖擊；彈簧分離裝置；對數(shù)正態(tài)分布；天地差異性

飛行器在發(fā)射和飛行期間，必須承受由于起飛、分離、關(guān)機和入軌等各種時序動作引起的噪聲、振動、沖擊等動力學(xué)環(huán)境[1]。美國航空航天局(NASA)的統(tǒng)計數(shù)據(jù)表明：高頻沖擊環(huán)境幾乎是所有航天器最嚴(yán)酷的力學(xué)環(huán)境之一,是導(dǎo)致飛行故障的主要環(huán)境因素，對飛行成功有著重要影響[2]；國內(nèi)也多次發(fā)生由于沖擊環(huán)境引起的電子設(shè)備故障和局部結(jié)構(gòu)損傷的情況，甚至導(dǎo)致飛行失敗。因此國內(nèi)外航天界對高頻瞬態(tài)沖擊開展了大量的研究，取得了豐富的研究成果[3]，歐美等國逐步建立了MIL-STD-7003、NASA-HDBK-7005、MIL-STD-810、火工品沖擊設(shè)計指南、ECSS、ACTEP等系列標(biāo)準(zhǔn)規(guī)范以及設(shè)計手冊，國內(nèi)也建立了GJB150、GJB/Z126等沖擊數(shù)據(jù)分析和試驗方法等方面的標(biāo)準(zhǔn)規(guī)范。

20世紀(jì)70年代后，隨著航空航天事業(yè)的發(fā)展，歐美逐步積累了一定數(shù)量的力學(xué)環(huán)境數(shù)據(jù)，開展了大量的數(shù)據(jù)分析研究工作。Piersol[4]在大量數(shù)據(jù)分析和工程經(jīng)驗的基礎(chǔ)上，提出在無其他有力證據(jù)的情況下，可認(rèn)為飛行器振動、噪聲和沖擊等動力學(xué)環(huán)境服從對數(shù)正態(tài)分布假設(shè)，并已被NASA-HDBK-7005動力學(xué)環(huán)境準(zhǔn)則采用[5]。William等[6]利用EOS Terra航天器地面星箭分離試驗數(shù)據(jù)，開展了大樣本爆炸沖擊環(huán)境的統(tǒng)計分析，確定了火工品爆炸沖擊環(huán)境的統(tǒng)計分布特性和分布規(guī)律，驗證了爆炸螺栓產(chǎn)生的高頻沖擊基本服從對數(shù)正態(tài)分布，并對比分析了該沖擊環(huán)境的天地差異性。Thomas[7]對比分析了Polaris導(dǎo)彈飛行試驗和地面試驗沖擊環(huán)境的相關(guān)性，驗證了地面試驗沖擊環(huán)境大于或等于飛行沖擊環(huán)境，并建議采用地面試驗獲得沖擊環(huán)境。國內(nèi)的沖擊環(huán)境研究也主要集中在切割索、爆炸螺栓、包帶等含火工品的分離裝置及其產(chǎn)生的沖擊環(huán)境[8～11]，研究彈簧等非火工品分離裝置產(chǎn)生的沖擊環(huán)境較少，同時在沖擊環(huán)境的天地差異性分析方面也少見報道。

本文針對某典型彈簧分離裝置及其產(chǎn)生的沖擊環(huán)境，搜集整理了多次地面試驗和飛行試驗中相同位置的沖擊環(huán)境數(shù)據(jù)，發(fā)現(xiàn)兩者之間存在一定的差異。因此，首先采用統(tǒng)計分析方法，分析該分離彈簧沖擊環(huán)境的分布規(guī)律，研究兩組樣本(地面試驗和飛行試驗)的分布形式是否相同，統(tǒng)計特性是否相似，散差量級是否一致；然后開展了地面試驗和飛行試驗沖擊環(huán)境的相關(guān)性研究，初步分析了天地差異性的主要原因。上述工作及其結(jié)果能為飛行器沖擊環(huán)境的統(tǒng)計分析、沖擊環(huán)境條件制定以及地面試驗設(shè)計等提供參考。

1 某彈簧分離裝置原理及其試驗

某分離方案中，采用彈簧作為分離儲能元件，當(dāng)分離前后體運動至一定相對位置時，即滿足分離條件，分離彈簧約束自動解除，預(yù)壓縮的彈簧開始動作，釋放儲存的能量，推動前后體產(chǎn)生Y向相對運動，實現(xiàn)兩體的側(cè)向分離。分離前體的結(jié)構(gòu)形式為薄壁圓柱殼體，中間部分艙段內(nèi)部安裝儀器設(shè)備，其殼體直徑為700 mm，壁厚5 mm，沖擊測點位于距分離彈簧安裝點200 mm的設(shè)備支架上(圖1)；分離后體為柱段結(jié)構(gòu)，內(nèi)部無儀器設(shè)備；彈簧預(yù)緊力為12 000 N±10%，標(biāo)稱預(yù)壓縮量為30 mm，安裝時根據(jù)預(yù)壓縮量對預(yù)緊力微調(diào)。

圖1 某彈簧分離裝置示意圖

飛行試驗中分離后體后端面為簡支狀態(tài)，分離前體為自由狀態(tài),僅通過彈簧分離裝置與分離后體連接，且分離前體為多艙段連接結(jié)構(gòu)。飛行中的分離過程為，分離前體在外界推力作用下，沿X軸向前運動，當(dāng)運動至某一相對位置時，彈簧約束解除，實現(xiàn)兩體側(cè)向分離。飛行試驗中受各種因素限制，使用8位數(shù)字采集的無線遙測系統(tǒng)，測量精度稍差。

而地面試驗中，受試驗設(shè)備、技術(shù)方法以及試驗規(guī)模等的限制，僅取分離前體的中間部分艙段(設(shè)備艙)，并固定該艙段的后端面，保證無剛體運動；利用地面試驗設(shè)備，沿X軸負(fù)向移動分離后體，模擬兩體的相對運動，保證分離時的相對運動速度與實際飛行狀態(tài)一致，最終模擬實際飛行狀態(tài)的兩體側(cè)向分離。試驗中，在該艙段內(nèi)儀器設(shè)備安裝支架上布置沖擊測點，利用測量精度較高的24位模數(shù)轉(zhuǎn)換編碼的數(shù)字采集有線測量系統(tǒng)，得到分離過程的沖擊數(shù)據(jù)。

2 各組沖擊環(huán)境數(shù)據(jù)的統(tǒng)計分析

2.1沖擊數(shù)據(jù)的統(tǒng)計分析方法

常用檢驗樣本是否服從正態(tài)分布的方法的包括圖方法(正態(tài)概率圖)、χ2檢驗、雅克比檢驗、Lilliefors檢驗、峰度/偏度檢驗、SHAPIRO-WILK檢驗等，各種檢驗方法均較為成熟，但檢驗結(jié)果存在一定的差別，詳細(xì)計算公式和過程見統(tǒng)計分析方面的相關(guān)文獻(xiàn)[12-13]。

沖擊環(huán)境的假設(shè)檢驗，是對描述沖擊環(huán)境的沖擊響應(yīng)譜進(jìn)行的，即在各頻點上均進(jìn)行假設(shè)檢驗，由于采用對數(shù)正態(tài)分布檢驗，因此需要對沖擊響應(yīng)譜值進(jìn)行對數(shù)計算。在對數(shù)正態(tài)分布的基礎(chǔ)上，可利用統(tǒng)計方法計算得到?jīng)_擊環(huán)境的P95/50(95%概率，50%置信度)預(yù)計的上下限，制定沖擊環(huán)境條件時，通常只使用P95/50上限。

2.2沖擊環(huán)境數(shù)據(jù)的統(tǒng)計分析

在彈簧分離裝置研制過程中，獲得了設(shè)備支架上2個沖擊測點的各6次沖擊數(shù)據(jù)，進(jìn)而可計算得到各沖擊數(shù)據(jù)的沖擊響應(yīng)譜(SRS)(圖2)，分析得出各次試驗數(shù)據(jù)之間的重復(fù)性較好，離散性較小。對于飛行試驗中獲得的設(shè)備支架6次沖擊數(shù)據(jù)，也可分析得到其的沖擊響應(yīng)譜(圖3)，可見各次飛行試驗數(shù)據(jù)的離散性稍大。將所有數(shù)據(jù)按照其最大值進(jìn)行歸一化后，分別如圖2和圖3所示，對比兩組沖擊環(huán)境數(shù)據(jù)(地面試驗與飛行試驗)，可見在量級和譜形方面還是存在一定的差別。注：歸一化方法為SRSnorm=20lg(SRSi/SRSmax)。

圖2 地面試驗測量數(shù)據(jù)的沖擊響應(yīng)譜

圖3 飛行試驗測量數(shù)據(jù)的沖擊響應(yīng)譜

分別利用χ2檢驗、雅克比檢驗、Lilliefors檢驗、峰度/偏度檢驗、SHAPIRO-WILK檢驗等方法，對上述兩組沖擊響應(yīng)譜數(shù)據(jù)進(jìn)行對數(shù)正態(tài)分布假設(shè)檢驗，檢驗結(jié)果如表1所示，可見大部分頻點的沖擊響應(yīng)譜均服從對數(shù)正態(tài)分布。根據(jù)兩組沖擊響應(yīng)譜數(shù)據(jù)的假設(shè)檢驗結(jié)果，大部分頻點的沖擊響應(yīng)譜均服從對數(shù)正態(tài)分布假設(shè)，少部分頻點不服從對數(shù)正態(tài)分布，因此可認(rèn)為沖擊環(huán)境近似服從對數(shù)正態(tài)分布假設(shè)。

3 沖擊環(huán)境的天地差異性分析

在確定兩組沖擊環(huán)境數(shù)據(jù)均近似服從對數(shù)分布之后，可開展地面與飛行沖擊環(huán)境的對比分析，研究兩組數(shù)據(jù)的統(tǒng)計特性是否一致，分析地面試驗沖擊環(huán)境與飛行沖擊環(huán)境之間的相關(guān)性，初步分析兩者之間的差異性原因。

3.1兩組沖擊環(huán)境數(shù)據(jù)一致性的分析方法

研究兩組沖擊環(huán)境數(shù)據(jù)相關(guān)性的方法有很多，在統(tǒng)計角度上常常對比均值、方差以及P95/50最大預(yù)計環(huán)境等參數(shù)，Thomas采用的兩組沖擊環(huán)境數(shù)據(jù)一致性的判斷標(biāo)準(zhǔn)為平均值相差20%以內(nèi)，方差量級相同等。William等采用統(tǒng)計參數(shù)進(jìn)行一致性判斷，包括分布形式相同、統(tǒng)計特性(均值和方差)相似(顯著水平取0.05)，即在服從對數(shù)正態(tài)分布的情況下，能通過均值T檢驗和方差F檢驗。

3.2兩組沖擊環(huán)境數(shù)據(jù)一致性的統(tǒng)計分析結(jié)果

3.2.1 均值、方差以及P95/50對比結(jié)果

基于對數(shù)正態(tài)分布，分別計算兩組數(shù)據(jù)的均值、P95/50上下限、方差等統(tǒng)計參數(shù)，并進(jìn)行對比(圖4)，可見2 000 Hz之前兩者的均值相差較大，方差均為2～3 dB；2 000 Hz之后兩者的均值基本吻合，飛行試驗稍大，且飛行試驗的方差顯著增大(4～6 dB)；此外兩者的P95/50統(tǒng)計上限相差稍大，主要原因是地面試驗樣本數(shù)稍多，重復(fù)性較好，而飛行試驗的樣本數(shù)較少，離散性較大。根據(jù)均值和方差判定準(zhǔn)則，不能認(rèn)為這兩組沖擊環(huán)境是一致的。

圖4 兩組數(shù)據(jù)的均值、方差及P95/50對比

3.2.2 均值T檢驗和方差F檢驗結(jié)果

由于兩組數(shù)據(jù)均可認(rèn)為近似服從對數(shù)正態(tài)分布，因此可通過檢驗兩組數(shù)據(jù)的均值、方差是否一致，以判斷兩組數(shù)據(jù)是否一致。檢驗結(jié)果(圖5)表明，在0.05顯著水平條件下，大部分頻點的均值T檢驗是拒絕兩組數(shù)據(jù)均值一致的假設(shè)，大部分頻點的方差F檢驗是接受兩組數(shù)據(jù)方差一致的假設(shè)，這表明兩組樣本存在顯著差異，即飛行試驗的沖擊環(huán)境與地面試驗的沖擊環(huán)境存在顯著差異，不能認(rèn)為兩者是一致的。但需要注意的是，兩組樣本的方差基本一致。

3.3地面試驗沖擊環(huán)境的覆蓋性分析

由于地面試驗與飛行試驗的沖擊環(huán)境存在顯著差異，因此需要確定地面試驗的沖擊環(huán)境是否能夠覆蓋飛行沖擊環(huán)境的散布范圍，對比了飛行試驗沖擊環(huán)境與地面試驗P95/50統(tǒng)計結(jié)果(圖6)，低頻部分低于地面試驗P95/50統(tǒng)計上限，高頻部分存在大于P95/50統(tǒng)計上限的情況，因此地面試驗的沖擊環(huán)境并不能完全覆蓋飛行沖擊環(huán)境，特別是高頻部分，這表明采用該地面試驗方案存在一定風(fēng)險，未能夠完全考核對儀器設(shè)備對沖擊環(huán)境的適應(yīng)性。

圖5 兩組數(shù)據(jù)的T檢驗和F檢驗結(jié)果

圖6 飛行試驗沖擊環(huán)境與地面試驗P95/50上下限對比

3.4沖擊環(huán)境天地差異的原因初步分析

在確定沖擊環(huán)境天地差異性較大后，需要研究分析導(dǎo)致差異的原因，為后續(xù)環(huán)境條件制定、試驗設(shè)計等提供參考。地面試驗和飛行試驗之間主要差異如表2所示，列出了產(chǎn)品技術(shù)狀態(tài)、彈簧蓄能、分離時間、邊界條件、測量系統(tǒng)等，可見主要影響因素是分離前體在地面試驗和飛行試驗中邊界條件差異較大，影響了其內(nèi)部儀器設(shè)備及支架的沖擊環(huán)境；測量系統(tǒng)精度不同，可能影響測量結(jié)果的準(zhǔn)確性。

表2 地面試驗和飛行試驗狀態(tài)的異同點

4 結(jié)論和建議

本文以某彈簧分離裝置產(chǎn)生的沖擊環(huán)境為例，開展了地面和飛行試驗沖擊環(huán)境的分布規(guī)律分析，同時利用統(tǒng)計特性、假設(shè)檢驗等方法，研究了地面試驗和飛行試驗沖擊環(huán)境的一致性，最后對影響沖擊環(huán)境天地差異性的主要因素進(jìn)行了初步分析，結(jié)果表明在沖擊環(huán)境方面，彈簧分離裝置與火工分離裝置之間即存在相同之處，也存在明顯差異：

(1) 分布特性檢驗結(jié)果表明，與火工分離裝置爆炸沖擊環(huán)境一致，該彈簧分離裝置的沖擊環(huán)境也近似服從對數(shù)正態(tài)分布，且散布較小，方差基本為3 dB左右，僅少數(shù)頻段達(dá)到6 dB；

(2) 天地一致性分析結(jié)果表明，與火工分離裝置爆炸沖擊環(huán)境不同，該彈簧分離裝置沖擊環(huán)境存在顯著的天地差異性，地面試驗沖擊環(huán)境未能完全覆蓋飛行沖擊環(huán)境，這說明不能經(jīng)驗性的人為地面沖擊環(huán)境不小于飛行試驗，因此在產(chǎn)品的環(huán)境適應(yīng)性、可靠性分析時，需要重視地面試驗的設(shè)計方案及對沖擊環(huán)境的影響；

(3) 影響該沖擊環(huán)境天地一致性的主要因素為分離前體的邊界條件以及測量系統(tǒng)的性能，建議開展更詳細(xì)的研究分析和對比試驗，同時也建議搜集更多其他分離裝置的沖擊環(huán)境，研究其分布特性和天地一致性。

[1] 馬興瑞，于登云，韓增堯，等. 星箭力學(xué)環(huán)境分析與試驗技術(shù)研究進(jìn)展[J]. 宇航學(xué)報，2006，27(3): 323-331.

MA Xingrui, YU Dengyun, HAN Zengyao, et al. Research evolution on the satellite-rocket mechanical environment analysis & test technology[J]. Journal of Astronautics, 2006, 27(3): 323-331.

[2] MOENING C J. Pyrotechnic shock flight failures[C]//Institute of Environmental Sciences Pyrotechnic Shock Tutorial Program, 31st Annual Technical Meeting, Inst. Envir. Sc., 985-04～05.

[3] 張建華. 航天產(chǎn)品的爆炸沖擊環(huán)境技術(shù)綜述[J].導(dǎo)彈與航天運載技術(shù), 2005(3):30-36.

ZHANG Jianhua. Pyroshock environment of missiles and launch vehicles[J]. Missiles and Space Vehicles, 2005(3): 30-36.

[4] PIERSOL A G. Review of procedures to compute maximum structural response from predictions or measurements at selected points[J]. Shock and Vibration, 1996, 3(3):211-221.

[5] NASA-HDBK-7005, NASA technical standard: Dynamic Environmental Criteria[S]. 2001.

[6] WILLIAM O H, ANNE M M. Statistical analysis of a large sample size pyroshock test data set including post flight data assessment[R]. NASA-TM-1998-206621.

[7] THOMAS J K. Correlation between Polaris ground and flight test shock and vibration data[R]. AIAA Paper NO. 65-206.

[8] 李會娜，黃友仙，肖健，等.某系列導(dǎo)彈爆炸沖擊環(huán)境統(tǒng)計[J]. 強度與環(huán)境, 2010, 37(6):12-17.

LI Huina, HUANG Youxian, XIAO Jian, et al. Statistical estimate of pyroshock environment for certain missiles[J]. Structure & Environment Engineering, 2010, 37(6):12-17.

[9] 丁繼鋒，趙欣，韓增堯. 航天器火工沖擊技術(shù)研究進(jìn)展[J]. 宇航學(xué)報, 2014, 35(12): 1339-1349.

DING Jifeng, ZHAO Xin, HAN Zengyao. Research development of spacecraft pyroshcok technique[J]. Journal of Astronautics, 2014,35(12):1339-1349.

[10] 段家希，李志來，曹乃亮，等. 包帶式解鎖支座沖擊環(huán)境研究[J].振動與沖擊, 2013,32(9):16-20.

DUAN Jiaxi, LI Zhilai, CAO Nailiang, et al. Shock environment of clamp-band unlocking support[J]. Journal of Vibration and Shock, 2013,32(9):16-20.

[11] 何春全,嚴(yán)楠,葉耀坤. 導(dǎo)彈級間火工分離裝置綜述[J]. 航天返回與遙感,2009, 30(3):47-49.

HE Chunquan, YAN Nan, YE Yaokun. Study of stage separation device for missile[J]. Spacecraft Recovery & Remote Sensing, 2009, 30(3):47-49.

[12] 吳贛昌. 概率論和數(shù)理統(tǒng)計[M]. 北京：中國人民大學(xué)出版社，2011.

[13] 數(shù)據(jù)的統(tǒng)計處理和解釋-正態(tài)性檢驗:GB/T 4882—2001[S]. 北京:中國標(biāo)準(zhǔn)出版社,2001.

Shockenvironmentstatisticalanalysisforaspringseparatedeviceandshockdatadifferencebetweenflighttestandgroundone

GAO Qing1,2， SHANG Lin2， SU Daliang3， ZHU Chen2， CAI Yipeng2， ZHAO Yonghui1

(1. State Key Laboratory of Mechanics and Control of Mechanical Structures， Nanjing University of Aeronautics and Astronautics， Nanjing 210016， China; 2. China Academy of Launch Vehicle Technology, Beijing 100076, China;3. China Aerospace Science and Technology Corporation, Beijing 100076, China)

In analyzing shock environment of a spring separate device, the shock data difference between flight test and ground one was found at the same position of one equipment bracket. It was shown through statistical analysis that the two shock environments obey the log-normal distribution; there is a significant difference between flight shock environment and ground shock one using a consistency check, but their variance levels are almost the same; boundary conditions and the measurement system are considered to be the major factors causing the difference. These conclusions provided a reference for shock environment statistical analysis, defining test specification and designing shock environment tests.

high frequency transient shock; spring separate device; log-normal distribution; difference between flight test and ground one

國家重點實驗室2015年開放課題(MCMS-0115G01)；國防技術(shù)基礎(chǔ)科研項目(JSZL2015203B002)

2016-03-30 修改稿收到日期：2016-07-13

高慶 男，碩士，高級工程師，1982年4月生

趙永輝 男，博士,教授,博士生導(dǎo)師，1969年11月生

V216.5

： A

10.13465/j.cnki.jvs.2017.17.031