裝備典型艦載平臺(tái)振動(dòng)環(huán)境嚴(yán)酷度分析

蔡健平，張萌，趙婉

（中國(guó)航天標(biāo)準(zhǔn)化研究所，北京 100071）

在滿足我軍實(shí)戰(zhàn)化要求和建設(shè)海洋強(qiáng)國(guó)的背景下，當(dāng)前各種艦載武器裝備正在從過去經(jīng)常近海待命和庫(kù)房貯存的狀態(tài)轉(zhuǎn)變?yōu)殚L(zhǎng)期在海上待機(jī)、值班和備戰(zhàn)的狀態(tài)，由于海風(fēng)、海浪以及發(fā)動(dòng)機(jī)等導(dǎo)致的艦船平臺(tái)誘發(fā)力學(xué)（如振動(dòng)、沖擊、搖擺）環(huán)境對(duì)武器裝備產(chǎn)生作用，從而對(duì)武器裝備振動(dòng)環(huán)境敏感的薄弱環(huán)節(jié)的有很大的影響[1—3]。例如仿真結(jié)果表明[2—3]：由于艦載誘發(fā)的力學(xué)環(huán)境長(zhǎng)期作用下導(dǎo)致的累積損傷，使固體發(fā)動(dòng)機(jī)粘接面在海上戰(zhàn)備1年的壽命比倉(cāng)庫(kù)貯存至少降低8.62%，而導(dǎo)彈的艦載貯存壽命只相當(dāng)于庫(kù)房壽命的15%。另外試驗(yàn)和分析也表明[4—5]：振動(dòng)或者熱導(dǎo)致的疲勞對(duì)發(fā)動(dòng)機(jī)藥柱和推進(jìn)劑會(huì)造成累積損傷，影響其壽命；另外振動(dòng)環(huán)境也會(huì)影響電子產(chǎn)品其可靠性，從而影響裝備的使用[7]。為了評(píng)價(jià)艦載裝備例如發(fā)動(dòng)機(jī)的貯存壽命，需要通過實(shí)驗(yàn)?zāi)M艦載振動(dòng)環(huán)境的作用[7]，或者以實(shí)際的振動(dòng)環(huán)境為依據(jù)設(shè)計(jì)振動(dòng)耐久性試驗(yàn)[8]。這就需要正確識(shí)別艦載環(huán)境條件，文中實(shí)測(cè)了裝備艦載典型測(cè)點(diǎn)的振動(dòng)環(huán)境，并對(duì)其環(huán)境嚴(yán)酷度進(jìn)行了評(píng)估，從而為相關(guān)武器裝備貯存試驗(yàn)和壽命評(píng)估提供依據(jù)。

1 實(shí)驗(yàn)和研究方法

1.1 典型艦載平臺(tái)選取

考慮到實(shí)測(cè)和搜集數(shù)據(jù)的可行性以及便利性，要求選擇常年在南海航行的補(bǔ)給船代表裝備實(shí)測(cè)艦載平臺(tái)。該艦共有4層甲板，包括主甲板、二甲板、艇甲板及駕駛和羅徑甲板，發(fā)動(dòng)機(jī)主機(jī)頻率為750～720 r/min。該艦正常工作航速較為穩(wěn)定，并主要在1～4級(jí)的海況下航行，因此該艦提供了一種長(zhǎng)期的、典型的艦載平臺(tái)環(huán)境，而非極端的力學(xué)環(huán)境。

1.2 艦載平臺(tái)振動(dòng)環(huán)境數(shù)據(jù)采集設(shè)備

采用定制的INV3062T2數(shù)據(jù)采集系統(tǒng)（如圖1所示），配合ICP加速度傳感器進(jìn)行振動(dòng)環(huán)境測(cè)量。數(shù)據(jù)采集分析系統(tǒng)具有下列技術(shù)指標(biāo)：頻率范圍為0.5～8 kHz；量程為50g；安裝諧振頻率＞25 kHz；分辨率為0.002 m/s2。

圖1 INV3062T2數(shù)據(jù)采集系統(tǒng)Fig.1 INV3062T2 data collecting system

1.3 艦載平臺(tái)振動(dòng)環(huán)境采集和分析方法

選取裝備服役的典型部位如艙室進(jìn)行振動(dòng)環(huán)境采集，采用ICP加速度傳感器同時(shí)采集船運(yùn)動(dòng)方向、船側(cè)向和豎直方向的加速度數(shù)據(jù)，采集頻率為12.8 kHz?？紤]到時(shí)間和成本，經(jīng)過初測(cè)后，選取振動(dòng)環(huán)境較為惡劣的點(diǎn)進(jìn)行測(cè)量，以便于進(jìn)行嚴(yán)酷度評(píng)價(jià)。這些點(diǎn)每次測(cè)量持續(xù)10 min或10 min以上，共進(jìn)行5次測(cè)量。獲得振動(dòng)數(shù)據(jù)后，先剔除異常值，再參照相關(guān)振動(dòng)分析的方法和標(biāo)準(zhǔn)進(jìn)行下列分析[9—11]。

1）振動(dòng)數(shù)據(jù)的時(shí)域分析。對(duì)振動(dòng)數(shù)據(jù)進(jìn)行時(shí)域分析時(shí)，將測(cè)得的加速度數(shù)據(jù)按照時(shí)間進(jìn)行繪制，就得到時(shí)域曲線，由于艦載振動(dòng)環(huán)境基本上接近穩(wěn)態(tài)，因此不需要進(jìn)行特殊處理，直接觀察時(shí)域曲線，就可以對(duì)數(shù)據(jù)的特性進(jìn)行判斷。

2）振動(dòng)數(shù)據(jù)的頻域分析。振動(dòng)數(shù)據(jù)的頻域分析采用快速傅立葉變換（FFT）進(jìn)行分析。FFT的基本思想就是，將1個(gè)長(zhǎng)序列依次分解為2個(gè)短序列來進(jìn)行離散傅里葉變換（DFT）。這里采樣頻率為12.8 KHz，分析頻率為6400 Hz，F(xiàn)FT點(diǎn)數(shù)取1024點(diǎn)，譜線條數(shù)512，采用Hanning窗，采用線性平均分析，分析主要通過DASP專業(yè)分析軟件完成。

3）概率分析[12—13]。分析不同加速度值出現(xiàn)的概率密度和概率。對(duì)在某范圍內(nèi)的加速度，其概率密度為：

式中：P|d1-d2|為d1到d2范圍內(nèi)加速度出現(xiàn)的概率密度；k|d1-d2|為d1到d2范圍內(nèi)加速度出現(xiàn)的個(gè)數(shù)；K為總的加速度個(gè)數(shù)。

同理，小于某加速度的概率為：

式中：P|-∞-d|為小于d的加速度出現(xiàn)的概率密度；k|-∞-d|為小于d的加速度出現(xiàn)的個(gè)數(shù)；K為總的加速度個(gè)數(shù)。

進(jìn)行概率分析時(shí)，讀取所有的時(shí)域數(shù)據(jù)，在Matlab環(huán)境按照（1），（2）進(jìn)行計(jì)算。

4）隨機(jī)振動(dòng)的加速度頻譜密度（ASD）。其定義為：

式中：αr.m.s·Δf是Δf頻率范圍內(nèi)加速度均方根值。

1.4 艦載環(huán)境嚴(yán)酷度評(píng)價(jià)方法

參照針對(duì)電子產(chǎn)品和機(jī)械產(chǎn)品的嚴(yán)酷度分級(jí)方法[14—15]，主要以GB/T 14091《機(jī)械產(chǎn)品環(huán)境參數(shù)分類及其嚴(yán)酷度分級(jí)》進(jìn)行振動(dòng)環(huán)境評(píng)級(jí)分類。對(duì)于隨機(jī)振動(dòng)可分成兩類：I類——有顯著的低頻量，頻率范圍為10～2000 Hz；II類——有較平坦分布的振動(dòng)能量，頻率范圍為20～2000 Hz。

隨機(jī)振動(dòng)的典型加速度頻譜分類見表1。

表1 隨機(jī)振動(dòng)的典型頻譜分類Table 1 Classification of typical frequency spectra of stochastic vibrations

2 結(jié)果和討論

2.1 船向振動(dòng)采集和分析結(jié)果

典型的船向上的振動(dòng)時(shí)域曲線如圖2所示，從45 ms的振動(dòng)時(shí)域曲線可見，船向上的振動(dòng)曲線是光滑過渡的，沒有尖刺或者較大的突變，因此測(cè)得的時(shí)域曲線是可信的隨機(jī)振動(dòng)曲線。

經(jīng)過FFT分析得到的船向振動(dòng)頻域曲線如圖3所示，對(duì)該圖中的頻率峰進(jìn)行進(jìn)一步分析，得到表2。從圖3和表2可見，船向振動(dòng)的加速度主峰大約在50Hz左右，加速度峰值約為0.412 72 m/s2，而在562.5 Hz和1112.5 Hz下各出現(xiàn)一個(gè)頻率峰，其加速度分別為0.152 34 m/s2和0.100 8 m/s2。

50 Hz的加速度主峰與交流電頻率一致，但在非艦載條件下，測(cè)量系統(tǒng)本身未見這個(gè)量級(jí)的加速度信號(hào)，可排除該處信號(hào)完全由測(cè)量系統(tǒng)的電源干擾引起，可能與測(cè)點(diǎn)附近使用的交流電的工作裝置如電機(jī)有關(guān)。562.5 Hz和1112.5 Hz的頻率與船上的旋轉(zhuǎn)裝置如油水分離器、轉(zhuǎn)軸等相關(guān)。

船向振動(dòng)加速度概率統(tǒng)計(jì)結(jié)果如圖4所示，從圖4中可見，船向振動(dòng)加速度接近正態(tài)分布，其數(shù)學(xué)期望或者均值為：

圖2 船向振動(dòng)時(shí)域曲線Fig.2 Time domain vibration curves in the ship direction

圖3 船向振動(dòng)頻域曲線Fig.3 Frequency domain vibration curve in the ship direction

表2 船向振動(dòng)頻域分布Table 2 Frequency domain distribution in the ship direction

式中：xi為第i個(gè)加速度值；p（xi）為xi出現(xiàn)的概率。正態(tài)分布的方差為：

從公式（4），（5）計(jì)算得到船向振動(dòng)加速度數(shù)學(xué)期望為 0.045 56 m/s2，方差為 0.267 57（m/s2）2。因此船向振動(dòng)加速度是一個(gè)正態(tài)隨機(jī)變量N（μ，σ2），其概率密度函數(shù)為：

為確定某x值的概率?？赏ㄟ^下列公式標(biāo)準(zhǔn)化為N（0，1）的標(biāo)準(zhǔn)正態(tài)分布后，查正態(tài)分布表即可：

圖4 船向振動(dòng)加速度概率統(tǒng)計(jì)Fig.4 Statistics of vibration acceleration in the ship direction

船向振動(dòng)加速度譜密度如圖5所示，船向振動(dòng)加速度譜密度峰值情況見表3。可以看出，船向振動(dòng)加速度譜密度峰值主要集中在50 Hz和562.5 Hz。其中50 Hz的振動(dòng)加速度譜密度為4.53×10-3（m/s2）2/Hz，562.5 Hz的振動(dòng)加速度譜密度為1.64×10-3（m/s2）2/Hz。

圖5 船向振動(dòng)加速度譜密度Fig.5 Spectral density of vibration acceleration in the ship direction

2.2 橫向振動(dòng)

采用2.1同樣的分析方法對(duì)橫向上的振動(dòng)進(jìn)行頻域和加速度譜密度分析，得到的結(jié)果見表4、表5。從表4和表5可見，橫向振動(dòng)的加速度主峰在562.5 Hz，加速度峰值約為0.812 m/s2，而在橫向主峰50 Hz處的加速度峰值較小，小于0.1 m/s2，1112.5 Hz頻率處的加速度很小。這和船向振動(dòng)的峰值也有不同，因此可以看出，在振動(dòng)在船向和橫向上有一定的擇優(yōu)性。另外在487.5 Hz到587.5 Hz區(qū)間還有其他小的頻率峰出現(xiàn)，和船向的頻率峰相比，這些頻率峰的量值與船向的頻率峰大致相當(dāng)。

表3 船向振動(dòng)加速度譜密度峰值情況Table 3 Acceleration spectral density peaks of vibration in the ship direction

表4 橫向振動(dòng)頻域分布Table 4 Frequency domain distribution in the transverse direction

表5 橫向振動(dòng)加速度譜密度峰值情況Table 5 Acceleration spectral density peaks of vibration in the transverse direction

根據(jù)方程（4），（5）得到橫向振動(dòng)加速度，其概率密度函數(shù)為：

表5中，562.5Hz的振動(dòng)加速度譜密度為17.59×10-3（m/s2）2/Hz，大約是船向數(shù)值的10倍。

2.3 垂直方向振動(dòng)

采用2.1同樣的分析方法對(duì)橫向上的振動(dòng)進(jìn)行頻域和加速度譜密度分析，得到的結(jié)果見表6、表7。從表6和表7可見，垂直振動(dòng)的加速度主峰在562.5 Hz，其加速度峰值約為1.560 24 m/s2，而在625，550，275，125，50 Hz下各有小的頻率峰，同橫向和船向的相比，562.5 Hz主峰的加速度量值最大。垂直振動(dòng)加速度譜密度峰值主要集中在562.5 Hz，其振動(dòng)加速度譜密度為97.37×10-3（m/s2）2/Hz，這個(gè)數(shù)值是橫向、船向數(shù)值中最大的。

表6 垂直方向振動(dòng)頻域分布Table 6 Frequency domain distribution in the vertical direction

表7 垂直振動(dòng)加速度譜密度峰值情況Table 7 Acceleration spectral density peaks of vibration in the vertical direction

根據(jù)方程（4），（5）得到垂直振動(dòng)加速度概率密度函數(shù)為：

2.4 振動(dòng)環(huán)境嚴(yán)酷度評(píng)價(jià)

艙室振動(dòng)環(huán)境的船向、橫向和垂直向的加速度譜密度、頻率分布均有較大差別，在垂直方向上主峰的加速度量值和加速度譜密度均大于船向、橫向上的相關(guān)數(shù)值。在艙室振動(dòng)環(huán)境的船向、橫向和垂直向上小于200 Hz下，其加速度譜密度均小于1（m/s2）2/Hz，而 200 Hz以上的加速度譜密度均小于 0.3（m/s2）2·Hz。因此艙室在三個(gè)方向上典型振動(dòng)環(huán)境屬于I型。由于垂直方向上加速度量值和加速度譜密度較大，所以考慮艦載平臺(tái)振動(dòng)環(huán)境時(shí)，應(yīng)把垂直方向作為考慮的重點(diǎn)。

3 結(jié)論

1）文中采用INV3062T2數(shù)據(jù)采集分析系統(tǒng)，配合ICP加速度傳感器對(duì)艦載平臺(tái)艙室的振動(dòng)環(huán)境進(jìn)行了測(cè)量。

2）艙室內(nèi)測(cè)得的振動(dòng)曲線是光滑過渡的，沒有尖刺或者較大的突變，是可信的隨機(jī)振動(dòng)曲線。各方向上的振動(dòng)加速度是一個(gè)正態(tài)隨機(jī)變量，并得到了各個(gè)方向上振動(dòng)加速度的概率密度函數(shù)。

3）船向振動(dòng)的加速度主峰大約在50 Hz左右，加速度峰值約為0.412 72 m/s2，而在562.5 Hz和1112.5 Hz下各出現(xiàn)1個(gè)頻率峰；橫向振動(dòng)的加速度主峰在562.5 Hz，加速度峰值約為0.812 m/s2；垂直振動(dòng)的加速度主峰在562.5 Hz，其加速度峰值約為1.560 24 m/s2，而在625，550，275，125，50 Hz下各有小的頻率峰。在垂直方向上的主峰的加速度量值和加速度譜密度均大于船向、橫向上的相關(guān)數(shù)值。

4）艙室在三個(gè)方向上典型振動(dòng)環(huán)境屬于I型。由于垂直方向上加速度量值和加速度譜密度較大，所以考慮艦載平臺(tái)振動(dòng)環(huán)境時(shí)，應(yīng)把垂直方向作為考慮的重點(diǎn)。

[1] 方書甲.海洋環(huán)境對(duì)海軍裝備性能的影響分析[J].艦船科學(xué)技術(shù)，2004，26（2）：5—10.FANG Shu-jia.Analysis of the Ocean Environment Influence on the Performance of the Navy Equipment[J].Ship Science and Technology，2004，26（2）：5—10.

[2] 曲凱，張杰，張旭東.艦船運(yùn)動(dòng)隊(duì)固體發(fā)動(dòng)機(jī)粘接面疲勞損傷研究[J].兵工學(xué)報(bào)，2012，33（8）：986—990.QU Kai，ZHANG Jie，ZHANG Xu-dong.Research on Effects of Ship Motion on Interface Fatigue Damage of Solid Rocket Motor[J].ACTA Armamentarii，2012，33（8）：986—990.

[3] 刑耀國(guó)，曲凱，許俊松，等.艦船搖擺條件下固體發(fā)動(dòng)機(jī)艦載壽命預(yù)估[J].推進(jìn)技術(shù)，2011，32（1）：32—35.XING Yao-guo，QU Kai，XU Jun-song，et al.Life Prediction of Shipboard Solid Rocket Motor under the Ship Swing Motion[J].Journal of Propulsion Technology，2011，32（1）：32—35.

[4] 劉著卿，李高春，丁彪，等.定應(yīng)力幅值往復(fù)拉伸過程中推進(jìn)劑的損傷擴(kuò)展[J].海軍航空工程學(xué)院學(xué)報(bào)，2011，26（3）：309—313.LIU Zhu-qing，LI Gao-chun，DING Biao，et al.Damage Characteristics of Propellant in the Process of Cyclic Tension under Certain Stress Amplitude[J].Journal of Navy Aeronautical University，2011，26（3）：309—313.

[5] 曲凱，刑耀國(guó)，張旭東.搖擺載荷作用下艦載固體發(fā)動(dòng)機(jī)藥柱疲勞損傷[J].航空動(dòng)力學(xué)報(bào)，2011，26（1）：2636—2640.QU Kai，XING Yao-guo，ZHANG Xu-dong.Fatigue Damage of Shipborne Solid Rocket Motor Propollent under Loading[J].Journal of Aerospace Power，2011，26（1）：2636—2640.

[6]靳哲峰.環(huán)境綜合應(yīng)力作用下航天電連接器可靠性分析和建模[D].杭州：浙江大學(xué)，2003.JIN Zhe-feng.Failure Analysis and Reliability Modeling for Electrical Connector under Action of Combined Environmental Stresses[D].Hangzhou：Zhejiang University，2003.

[7] 徐金洲，賈耀明.艦上導(dǎo)彈固體發(fā)動(dòng)機(jī)貯存壽命分析方法研究[J].固體火箭技術(shù)，2009，32（3）：271—273 XU Jin-zhou，JIA Yao-ming.Research of Analysis Method for Storage Life of SRM[J].Journal of Solid Rocket Technology，2009，32（3）：271—273

[8]GJB 150A，軍用裝備實(shí)驗(yàn)室環(huán)境試驗(yàn)方法，第16部分：振動(dòng)試驗(yàn)[S].GJB 150A，Laboratory Environmental Test Method for Military Materiel-Part 16：Vibration Test[S].

[9] 王孚懋，任勇生，韓寶坤，等.機(jī)械振動(dòng)與噪聲分析基礎(chǔ)[M].北京：國(guó)防工業(yè)出版社，2009.WANG Fu-mao，REN Yong-sheng，HAN Bao-kun，et al.Mechanical Vibration and Noise Analysis[M].Beijing：Defense Industry Press，2009.

[10]徐天成，谷亞林，錢玲.信號(hào)與系統(tǒng)[M].北京：電子工業(yè)出版社，2012.XU Tian-cheng，GU Ya-lin，QIAN Ling.Signal and System[M].Beijing：Electronic Industry Press，2012.

[11]GJB 10593.3，電工電子產(chǎn)品環(huán)境參數(shù)測(cè)量方法——振動(dòng)數(shù)據(jù)處理與歸納[S].GJB 10593.3，Methods of Measuring Environmental Parameters for Electrical and Electronic Products，Treatment and Induction for Vibration Data[S].

[12]吳贛昌.概率論和數(shù)理統(tǒng)計(jì)[M].北京：中國(guó)人民大學(xué)出版社，2011.WU Gan-chang.Probability and Mathematical Statistics[M].Beijing：Press of Renmin University，2011.

[13]峁詩松，湯銀才，王玲玲.可靠性統(tǒng)計(jì)[M].北京：高等教育出版社，2008.MAO Shi-song，TANG Yin-cai，WANG Ling-ling.Reliability Statistics[M].Beijing：Advanced Education Press，2008.

[14]GB/T 14091，機(jī)械產(chǎn)品環(huán)境參數(shù)分類及其嚴(yán)酷度分級(jí)[S].GB/T 14091，Classification and Severity of Mechanical Product Environmental Factors[S].

[15]GB/T 4769，電子電工產(chǎn)品環(huán)境條件分類[S].GB/T 4769，Classification of Environmental Conditions of Electric and Electronic Products[S].