晏廷飛，沈志強，張俊剛,2

（1. 北京衛(wèi)星環(huán)境工程研究所；2. 北京衛(wèi)星環(huán)境工程研究所 可靠性與環(huán)境工程技術(shù)重點實驗室：北京100094）

一種氣動式?jīng)_擊響應(yīng)譜試驗系統(tǒng)的設(shè)計與分析

晏廷飛1，沈志強1，張俊剛1,2

（1. 北京衛(wèi)星環(huán)境工程研究所；2. 北京衛(wèi)星環(huán)境工程研究所 可靠性與環(huán)境工程技術(shù)重點實驗室：北京100094）

大多數(shù)航天器單機產(chǎn)品都要求進行沖擊響應(yīng)譜試驗，試驗量級普遍為1000g～2000g，甚至超過3000g。對于較大質(zhì)量（如50 kg以上）產(chǎn)品高量級的沖擊試驗，振動臺、傳統(tǒng)的擺錘式或跌落式?jīng)_擊試驗臺均很難滿足相關(guān)要求。文章研究并設(shè)計了一套能夠進行大質(zhì)量受試產(chǎn)品高量級沖擊響應(yīng)譜試驗的氣動沖擊試驗系統(tǒng)。該系統(tǒng)利用壓縮空氣瞬間釋放膨脹推動質(zhì)量塊加速撞擊具有多階固有頻率的諧振板，通過諧振板被激起的響應(yīng)模擬復(fù)雜的衰減正弦波。測試結(jié)果表明，系統(tǒng)空載時沖擊譜量級達8000g，負載200 kg時可達5000g，時域曲線為振蕩衰減波，持續(xù)時間小于10 ms。文章提出的氣動式?jīng)_擊響應(yīng)譜試驗系統(tǒng)設(shè)計方法可為此類沖擊試驗系統(tǒng)的設(shè)計提供參考和理論依據(jù)。

氣動；沖擊響應(yīng)譜試驗；試驗系統(tǒng)設(shè)計

0 引言

隨著上面級分離等一些新技術(shù)、新結(jié)構(gòu)和新方法在航天器研制中的應(yīng)用，航天器發(fā)射時所經(jīng)歷的沖擊環(huán)境更為嚴酷，部分組件的沖擊量級高達4000g以上，傳統(tǒng)的試驗理論和方法日漸不能滿足需求。加之GJB 1027A新標準的發(fā)布，航天產(chǎn)品在研制階段沖擊試驗要求用沖擊響應(yīng)譜（Shock Response Spectrum, SRS）代替經(jīng)典沖擊譜[1-4]。

當前，航天器試驗所需的 SRS通常采用振動臺激勵方式、擺錘式?jīng)_擊和跌落式?jīng)_擊來實現(xiàn)[5-6]。目前我國比較好的擺錘式?jīng)_擊臺其最大負載不超過100 kg，拐點頻率為300～1000 Hz，滿載時沖擊譜量級不超過2000g。面對嚴酷的SRS試驗要求（如試件質(zhì)量＞100 kg、試驗量級＞2000g），傳統(tǒng)的試驗方式均很難達到。

本文基于壓縮氣體瞬間膨脹產(chǎn)生高速氣動沖擊的理論，通過理論分析、計算并結(jié)合工程實踐經(jīng)驗，給出了一套沖擊試驗的設(shè)計理論和方法，并完成了大質(zhì)量、大尺寸試件以及高量級試驗條件的水平式氣動沖擊響應(yīng)譜試驗系統(tǒng)研制[7-8]。

1 基本構(gòu)成及工作原理

水平式氣動沖擊響應(yīng)譜模擬試驗系統(tǒng)如圖 1所示。其工作原理是利用氣體瞬間釋放膨脹推動質(zhì)量塊加速撞擊具有多階固有頻率的諧振板并激起諧振，調(diào)整諧振板前、后波形墊以及沖擊力的作用效果，使諧振板被激起的響應(yīng)近似于復(fù)雜的衰減正弦波，以模擬爆炸沖擊環(huán)境。該系統(tǒng)的沖擊動力組件由沖擊錘、活塞推桿、儲氣缸、沖擊氣缸等構(gòu)成，其中儲氣缸和沖擊氣缸為沖擊試驗提供動力源。諧振板除了產(chǎn)生滿足譜型要求的響應(yīng)外，還是受試產(chǎn)品的安裝接口[9-11]。質(zhì)量塊是沖擊試驗的運動部件，在壓縮空氣膨脹作用下高速運動撞擊諧振板。

2 理論設(shè)計與分析

2.1 基本理論

將上述復(fù)雜的氣動沖擊試驗系統(tǒng)簡化為若干個不同的單自由度彈簧-質(zhì)量系統(tǒng)，對每個系統(tǒng)進行沖擊響應(yīng)譜分析計算，最后加以合成，即可得到整個試驗系統(tǒng)的沖擊響應(yīng)譜。圖2為水平式氣動沖擊試驗系統(tǒng)中的單自由度質(zhì)量-彈簧-阻尼模型[12]。

該模型的運動方程為

式中：y(t)為系統(tǒng)的響應(yīng)；x(t)為外部輸入激勵；m為諧振板質(zhì)量；k為系統(tǒng)的剛度；fn為系統(tǒng)固有頻率；c為系統(tǒng)阻尼；ζ為系統(tǒng)的阻尼系數(shù)

當諧振板的初始位移、初始速度為0時，式(1)的解為

顯然，y(t)是fn和ζ的函數(shù)，當ζ值一定時，可得到˙( t)的最大值對fn的變化曲線，即加速度沖擊響應(yīng)譜。在fn確定的情況下，通過調(diào)整系統(tǒng)的阻尼和沖擊時的瞬時速度即可得到滿足要求的加速度沖擊響應(yīng)譜曲線。對于氣動沖擊試驗系統(tǒng)而言，關(guān)鍵是諧振板及沖擊動力組件的設(shè)計和分析。

2.2 諧振板設(shè)計與分析

諧振板設(shè)計主要考慮以下因素：試件質(zhì)量、尺寸、頻率范圍、量級、拐點頻率、安裝方式、材料阻尼系數(shù)、幾何形狀等。與擺錘式?jīng)_擊試驗系統(tǒng)的設(shè)計方法基本一致，需要進行模態(tài)分析、頻響分析、強度分析及沖擊響應(yīng)譜分析等。其中關(guān)鍵的是要能實現(xiàn)拐點頻率可調(diào)。根據(jù)沖擊響應(yīng)譜試驗機的設(shè)計原理，只有諧振板在某種確定的邊界條件以及拐點頻率附近具有某階共振頻率時，通過調(diào)整沖擊力的大小及前后波形發(fā)生器才能實現(xiàn)規(guī)定的沖擊響應(yīng)譜試驗規(guī)范。

諧振板設(shè)計及分析方法如下：

1）確定諧振板材料。

2）通過諧振板有限元分析，并根據(jù)其自由狀態(tài)的一階主模態(tài)低于最低拐點頻率的要求確定其外形尺寸。

3）通過施加不同的位移約束邊界，仿真分析確定其與底座的安裝方式。

4）用最大設(shè)計沖擊量級進行強度校核，并要求有足夠的安全余量。

2.3 沖擊動力組件設(shè)計與分析

在沖擊動力組件的設(shè)計中應(yīng)確定的參數(shù)主要有沖擊錘的質(zhì)量，沖擊氣缸、儲氣缸直徑，以及活塞推桿的行程。

設(shè)諧振板的固有頻率為fn，受到?jīng)_擊時產(chǎn)生的峰值加速度為 Anmax，沖擊錘撞擊諧振板產(chǎn)生的脈沖載荷持續(xù)時間近似等于固有頻率周期的一半，即

其中t為激勵的最小脈沖時間。根據(jù)Anmax可確定諧振板的速度為

其中δ為修正因子，一般在2.5～5之間取值。

假設(shè)沖擊氣缸內(nèi)活塞桿的有效運動行程為S，沖擊氣缸內(nèi)徑為d，儲氣缸內(nèi)徑為Φ，開始時刻儲氣缸壓力為 P0，中間任意時刻的壓力為 Pt，對應(yīng)的沖擊活塞運動的距離是St，則氣缸內(nèi)平均氣壓為P，如圖3所示。

在沖擊過程中，假設(shè)儲氣缸內(nèi)氣壓穩(wěn)定，則沖擊氣缸內(nèi)中間任意時刻的壓力為

根據(jù)力的平衡原理則有

式中：mc為沖擊錘的總質(zhì)量；Pa為標準大氣壓；μ為沖擊組件運動過程中的摩擦系數(shù)。

求解式(5)、式(6)，得到氣缸位于 St位置時沖擊錘的速度vt為

則沖擊錘撞擊臺面時（St=S）的速度vc為

假設(shè)諧振板及產(chǎn)品的總質(zhì)量為 mT，沖擊錘撞擊臺面過程遵循彈性碰撞理論，碰撞過程中的能量損耗因子為ξ，則根據(jù)能量守恒定律有

綜合以上各式，假設(shè)沖擊氣缸與儲氣缸內(nèi)徑相同，則可得

因此，根據(jù)最大試驗量級可由式(10)確定沖擊錘質(zhì)量、沖擊氣缸的直徑和沖擊錘運動行程。

3 設(shè)計方案的仿真與試驗驗證

3.1 有限元仿真分析

下面以1.25 m×1.25 m（長×寬）、材料為鋁的諧振板為例進行試驗驗證。設(shè)計指標為：最大負載200 kg，滿負載最大試驗量級 5000g，拐點頻率400～1500 Hz間可調(diào)。

首先進行諧振板有限元優(yōu)化分析（FEA）計算[13]，得到其優(yōu)化厚度為80 mm，采用3根導(dǎo)軌6個滑塊的安裝方式。FEA模型中滑塊的約束條件為運動方向自由，其余方向為固支。諧振板質(zhì)量為216 kg，它的前 15階模態(tài)頻率見表 1，其中一階固有頻率約為389 Hz；一階模態(tài)振型見圖4。

表1 諧振板前15階模態(tài)頻率Table 1 Natural frequencies of the first fifteen orders for the syntonic plate

其次確定計算常量。取δ =2.5，μ =0.2，P0= 0.7 MPa，Pa= 0.1 MPa，ξ=0.20.7 MPa。滿負載時，mT= 216+200 =416 kg。

一般來講，高量級沖擊的拐點頻率比低量級沖擊的要高。假設(shè)沖擊量級為 5000g、拐點頻率為1200 Hz，與沖擊量級為2000g、拐點頻率為400 Hz所需的沖擊能量相當，則式(10)可以簡化為

根據(jù)式(11)進行不同參數(shù)組合分析計算，結(jié)果見表2。

表2 mc、d和S的取值關(guān)系Table 2 Relationship of mc, d and S

如mc=700 kg，沖擊氣缸的直徑d = 0.2 m，則對應(yīng)活塞推桿的運動行程S = 0.781 m，代入式(8)得到?jīng)_擊組件撞擊臺面時的速度約為7.2 m/s，遠高于目前國內(nèi)擺錘沖擊臺的撞擊速度，因此可獲得較高的沖擊能量與沖擊試驗量級。

3.2 試驗驗證

根據(jù)上述設(shè)計參數(shù)研制了一臺諧振板為1.25 m×1.25 m的氣動沖擊響應(yīng)譜試驗系統(tǒng)（見圖5）并進行了試驗測試。測試結(jié)果表明，該氣動沖擊試驗系統(tǒng)可實現(xiàn)400～1500 Hz內(nèi)拐點可調(diào)，空載時試驗量級可達8000g。圖6為該系統(tǒng)在負載200 kg時，試驗條件為 100～1500 Hz、+6 dB/oct、 1500～5000 Hz、5000g的測試曲線。可看出，沖擊響應(yīng)試驗譜符合規(guī)范要求，時域曲線近似爆炸沖擊的振蕩衰減譜。

該沖擊試驗系統(tǒng)目前只能進行水平 2個方向的沖擊試驗，如果條件允許，則可通過設(shè)計翻轉(zhuǎn)工裝或者L型支架實現(xiàn)第3個方向的沖擊試驗。

4 結(jié)束語

本文提出了一種氣動式高量級沖擊響應(yīng)譜試驗系統(tǒng)的設(shè)計分析方法。相較傳統(tǒng)的機械式?jīng)_擊響應(yīng)譜試驗系統(tǒng)，氣動式能獲得較大的瞬態(tài)沖擊速度，可有效提高試驗量級。針對本文的設(shè)計方法與參數(shù)研制了一臺氣動式?jīng)_擊試驗系統(tǒng)并進行了測試試驗。測試結(jié)果表明，時域曲線近似爆炸沖擊響應(yīng)譜的瞬態(tài)振蕩衰減譜，沖擊時間可嚴格控制在10 ms內(nèi)。因此，可根據(jù)此方法及用戶需要，設(shè)計出滿足試驗規(guī)范要求的氣動式水平?jīng)_擊試驗系統(tǒng)。

該系統(tǒng)目前不能用于某些不具備翻轉(zhuǎn)條件的產(chǎn)品3個方向的沖擊試驗，后續(xù)將進一步研究可進行垂直向沖擊譜試驗系統(tǒng)的設(shè)計方法。

（References）

[1] Pyro-shock test criteria: NASA-STD-7003[S], 1999-05-18

[2] Environmental engineering laboratory tests: MIL-STD-810F[S], 2000-01-01

[3] 運載器、上面級和航天器試驗要求: GJB 1027A—2005[S], 2005-12-12

[4] 張建華. 航天產(chǎn)品的爆炸沖擊環(huán)境技術(shù)綜述[J]. 導(dǎo)彈與航天運載技術(shù), 2005(3): 30-35 ZHANG J H. Pyroshock environment of missiles and launch vehicles[J]. Missiles and Space Vehicles,2005(3): 30-35

[5] 薄浩. 擺錘式?jīng)_擊試驗系統(tǒng)開發(fā)[D]. 北京: 北京工業(yè)大學, 2006: 18-25

[6] 徐鵬, 祖靜, 范錦彪. 高 G 值加速度沖擊試驗技術(shù)研究[J]. 振動與沖擊, 2011, 30(4): 241-244 XU P, ZU J, FAN J B. Acceleration shock test technology with higher values of G[J]. Journal of Vibration and Shock, 2011, 30(4): 241-244

[7] 吳斌, 李立名. 雙向式水平?jīng)_擊試驗臺設(shè)計[J]. 振動、測試與診斷, 2002, 22(4): 310-314 WU B, LI L M. Design of dual-horizontal shock test apparatus[J]. Journal of Vibration Measurement &Diagnosis, 2002, 22(4): 310-314

[8] 王飛, 皮本樓, 馮偉干. 氣動式爆炸沖擊環(huán)境模擬裝置試驗技術(shù)研究[J]. 強度與環(huán)境, 2011, 38(4): 16-19 WANG F, PI B L, FENG W G. Experiment research of pneumatic pyroshock environment simulation device[J].Structure & Environment Engineering, 2011, 38(4):16-19

[9] SMITH S, HOLLOWELL B. A proposed method to standardize shock response spectrum (SRS) analysis[J].Journal of the Institute of Environmental Sciences, 1996,39(3): 19-24

[10] 李含忠, 王招霞. 諧振板式大型試件沖擊譜模擬技術(shù)研究[J]. 航天器環(huán)境工程, 2008, 25(2): 159-162 LI H Z, WANG Z X. SRS simulation test using resonant large plates[J]. Spacecraft Environment Engineering,2008, 25(2): 159-162

[11] 馬愛軍, 馮雪梅, 盧來潔, 等. 一種沖擊譜試驗機的研制[J]. 航天醫(yī)學與醫(yī)學工程, 2003, 16(4): 296-299 MA A J, FENG X M, LU L J, et al. Development of a shock response spectrum test facility[J]. Space Medicine& Medical Engineering, 2003, 16(4): 296-299

[12] 柯受全. 衛(wèi)星環(huán)境工程與模擬試驗[M]. 北京: 宇航出版社, 1996: 105-112

[13] 安梅巖, 馬愛軍, 劉洪英. 大載荷諧振板的仿真研究[J].航天器環(huán)境工程, 2007, 24(6): 386-389 AN M Y, MA A J, LIU H Y. Numerical simulation of a resonant plate with heavy load[J]. Spacecraft Environment Engineering, 2007, 24(6): 386-389

（編輯：許京媛）

Design and analysis of a pneumatic test system for shock response spectrum

YAN Tingfei1, SHEN Zhiqiang1, ZHANG Jungang1,2
(1. Beijing Institute of Spacecraft Environment Engineering;2. Science and Technology on Reliability and Environmental Engineering Laboratory,Beijing Institute of Spacecraft Environment Engineering: Beijing 100094, China)

The impact test for spacecraft is generally required with a magnitude between 1000g-2000g,sometime even more than 3000g. In a high-magnitude impact test for large (＞500 kg) products, it is difficult for the pendulum impact test equipment and the drop impact test equipment to meet the test requirements. In order to solve these problems, a pneumatic impact test system is designed for the high-magnitude impact test. The system uses the compressed air as the power source, and the instantaneous release of gas to expand the impact mass block, which is accelerated to impact the resonant plate, so that a spectrum similar to the complex damped sine wave is generated in the resonant plate. The test results show that when the system is unloaded, the impact value is up to 8000g. When the load is 200 kg, the impact value is up to 5000g. The time domain curve is the shock attenuation wave, and the duration is less than 10ms. The method proposed in this paper can serve as a reference and theoretical basis for the design of this kind of impact test systems.

pneumatic; shock response spectrum test; test system design

TH213.3

A

：1673-1379(2017)03-0290-05

10.3969/j.issn.1673-1379.2017.03.011

晏廷飛（1975—），男，碩士學位，高級工程師，主要從事航天器動力學環(huán)境試驗與模擬技術(shù)研究。E-mail:yantingfei@tsinghua.org.cn。

2017-02-21；

2017-05-22

總裝備部共性技術(shù)課題（編號：51334050101）

晏廷飛, 沈志強, 張俊剛. 一種氣動式?jīng)_擊響應(yīng)譜試驗系統(tǒng)的設(shè)計與分析[J]. 航天器環(huán)境工程, 2017, 34(3)：290-294

YAN T F, SHEN Z Q, ZHANG J G, et al. Design and analysis of a pneumatic test system for shock response spectrum[J]. Spacecraft Environment Engineering, 2017, 34(3)： 290-294