周學(xué)濱 張永坤

（91439部隊(duì) 大連 116041）

1 引言

沖擊是指一個系統(tǒng)在相當(dāng)短的時間內(nèi)，受到瞬態(tài)激勵，其位移、速度或加速度突然發(fā)生變化的物理現(xiàn)象。在沖擊激振力作用下，系統(tǒng)所產(chǎn)生的運(yùn)動為瞬態(tài)運(yùn)動，而當(dāng)激振力作用停止后，系統(tǒng)將按其自身固有頻率作自由振動。艦船設(shè)備的沖擊環(huán)境主要表現(xiàn)為設(shè)備基礎(chǔ)所經(jīng)受的沖擊加速度、沖擊速度和位移，是設(shè)備抗沖擊研究的基礎(chǔ)。研究大中型艦船設(shè)備的沖擊環(huán)境方法有兩種：實(shí)船抗沖擊試驗(yàn)和沖擊平臺試驗(yàn)。由于實(shí)船試驗(yàn)的組織實(shí)施難度大，費(fèi)用高，因此多采用浮動沖擊平臺水下爆炸試驗(yàn)方法，通常需要選擇中型浮動沖擊平臺作為載體［1～2］。借助浮動沖擊平臺進(jìn)行艦船設(shè)備抗沖擊性能考核，國內(nèi)作了很多研究［3～8］?？己诉^程中，浮動沖擊平臺主要作用是為安裝于其內(nèi)部的艦船設(shè)備提供符合要求的沖擊環(huán)境，研究中型浮動沖擊平臺的沖擊環(huán)境對于建設(shè)一型標(biāo)準(zhǔn)平臺，為設(shè)備考核提供更加真實(shí)的模擬實(shí)船環(huán)境，具有重要的意義。

文中選取一型標(biāo)準(zhǔn)浮動沖擊平臺為計算對象，采用Abaqus軟件，對其標(biāo)準(zhǔn)工況下的沖擊環(huán)境及毀傷參數(shù)進(jìn)行計算，分析了不同爆炸距離、不同爆炸攻角、不同沖擊因子作用下的沖擊環(huán)境參數(shù)的變化規(guī)律，將毀傷量值與理論分析算法進(jìn)行比對，沖擊環(huán)境結(jié)果與BV043/85標(biāo)準(zhǔn)規(guī)定的算法進(jìn)行比對，對爆炸作用下標(biāo)準(zhǔn)浮動沖擊平臺的毀傷程度及沖擊環(huán)境進(jìn)行了分析研究。

2 仿真計算過程

將中型浮動沖擊平臺置于MIL-S-901D中的標(biāo)準(zhǔn)工況下進(jìn)行數(shù)值計算，檢查塑性形變的大小以考核平臺的強(qiáng)度。加載時使用Geer and hunter模型計算壓力載荷［9～10］。

2.1 仿真方法驗(yàn)證

結(jié)合某實(shí)船爆炸數(shù)據(jù)，對仿真方法進(jìn)行了驗(yàn)證［11］。計算工況為1000kg TNT爆源距左舷70m，水深50m沉底爆炸，龍骨沖擊因子為0.337。根據(jù)試驗(yàn)時測點(diǎn)的實(shí)際情況選取進(jìn)行比對的測點(diǎn)位置，分別選擇了3個應(yīng)變測點(diǎn)和3個加速度測點(diǎn)。測點(diǎn)具體位置分別為01甲板55#肋位迎爆面中縱桁縱向應(yīng)變測點(diǎn)、主甲板25#肋位迎爆面縱桁縱向應(yīng)變測點(diǎn)、內(nèi)底45#肋位迎爆面第一扶強(qiáng)材邊中間根部板格應(yīng)變測點(diǎn)、01甲板27#肋位艦長室垂向加速度測點(diǎn)、主甲板76#肋位垂向加速度測點(diǎn)、內(nèi)底47#肋位垂向加速度測點(diǎn)。對比結(jié)果見表1～2。

表1 應(yīng)變峰值對比

表2 加速度峰值對比

通過對峰值的相對誤差計算，得出應(yīng)變峰值平均精度77.48%，加速度峰值平均精度77.52%，平均相對誤差均在30%以內(nèi)，計算精度符合工程要求。

2.2 仿真計算對象

仿真計算對象模型如圖1所示。

2.3 仿真工況設(shè)計

為能夠較全面地分析平臺的工況，設(shè)置了80個工況進(jìn)行計算。設(shè)置方法如下：

1）選取4個平面。這4個平臺為x=6.5，x=3.25，x=0，及y=3.25，4個平面位置在下文中簡稱1號位置、2號位置、3號位置、4號位置，如圖2（a）。

圖1 中型浮動沖擊平臺數(shù)值模型

2）在各平面內(nèi)選取4個扇形區(qū)域。以美國MIL-S-901D為基礎(chǔ)，參考4.1，3.1中4個工況選取4個沖擊因子，分別為C1＝0.51，C2＝0.63，C3＝0.73，C4＝0.86，對應(yīng) 爆 距為 R1＝14.39m，R2＝11.61m，R3＝10.01m，R4＝8.51m。以坐標(biāo)原點(diǎn)為圓心，以上述4個爆距為半徑，作4個扇形區(qū)域如圖2（b）。

圖2 工況設(shè)置示意圖

3）作5條射線。在各扇形內(nèi)分別取α1＝0°，α2＝30°，α3＝45°，α4＝60°，α5＝90°的5種攻角（圖中只標(biāo)出α2＝30°），作5條射線，如圖2（b）。

4）設(shè)置爆源位置。在各平面內(nèi)的90°扇形內(nèi)各有20個交點(diǎn)，共4個平面位置，則共有80個交點(diǎn)。將爆源分別放在這80個交點(diǎn)上，就獲得80個不同的爆炸工況。

2.4 測點(diǎn)選取

浮動沖擊平臺內(nèi)底是待測設(shè)備的主要安裝基礎(chǔ)，其沖擊環(huán)境即為設(shè)備（或基座）的沖擊輸入。本文在平臺內(nèi)底中央選擇36個點(diǎn)作為測點(diǎn)，這些點(diǎn)所占的面積接近1m2，并在平臺上建立平面直角坐標(biāo)系，如圖3，原點(diǎn)位于圖中平臺外底左下角。

圖3 平臺內(nèi)底測點(diǎn)的選取

2.5 仿真計算結(jié)果

2.5.1 平臺抗沖擊強(qiáng)度計算結(jié)果

按2.3所設(shè)的工況，對中型浮動沖擊平臺結(jié)構(gòu)模型在Abaqus中進(jìn)行仿真計算。根據(jù)結(jié)構(gòu)毀傷特征，沖擊因子最大，所受到的沖擊也最為強(qiáng)烈，因此若平臺在該工況下強(qiáng)度滿足要求，可正常工作，即可認(rèn)為為該平臺的結(jié)構(gòu)強(qiáng)度滿足設(shè)計要求。圖4為工況4下（沖擊因子為0.86）的平臺受沖擊后的位移云圖，從該圖上可以看出：輕質(zhì)舷墻的頂部為平臺發(fā)生最大形變的區(qū)域，形變最大的位置為迎爆面的兩個頂角。各工況下平臺最大變形情況如表3所示。

圖4 中型浮動沖擊平臺在工況4時的橫向位移云圖

表3 平臺最大變形隨沖擊因子變化關(guān)系

由表3中可知，四種沖擊因子（0.51，0.63，0.73，0.86）狀況下，底部中間的最大變形量值分別為0.026cm，0.04cm，0.056cm，0.075cm。

2.5.2 沖擊環(huán)境計算結(jié)果

對沖擊平臺在上述共計80個工況下進(jìn)行仿真計算后，得到各測點(diǎn)處的沖擊環(huán)境參數(shù)，給出將每個工況各測點(diǎn)的沖擊環(huán)境相同參數(shù)進(jìn)行平均處理。為了消除位置參數(shù)的影響，選擇最危險的截面對沖擊響應(yīng)特征進(jìn)行分析，給出1號位置x=6.5m處的沖擊參數(shù)，如表4所示。

3 計算結(jié)果比對分析

3.1 沖擊環(huán)境特征分析

毀傷的量值在2.5.1的計算結(jié)果中已經(jīng)給出，此處根據(jù)表4沖擊環(huán)境的計算結(jié)果對沖擊環(huán)境的特征進(jìn)行分析，重點(diǎn)說明不同攻角情況下的變化特征。比較表4的數(shù)據(jù)可知，0°攻角，沖擊因子增加，譜位移增大，譜速度增大，譜加速度增大；30°攻角，沖擊因子增加，譜位移增大，譜速度增大，譜加速度增大；45°攻角，沖擊因子增加，譜位移增大，譜速度增大，譜加速度增大；60°攻角，沖擊因子增加，譜加速度增大；譜位移和譜速度在0.73出現(xiàn)最大值；90°攻角，沖擊因子增加，譜位移減小，譜速度和譜加速度變化不規(guī)律。

表4 x=6.5m位置處各工況的沖擊環(huán)境參數(shù)

可總結(jié)為沖擊因子相同情況下，沖擊加速度隨攻角增加而增大，在90°時達(dá)到最大值；沖擊因子相同情況下，沖擊速度及譜位移隨攻角變化不規(guī)律，需要綜合沖擊因子量值及攻角情況綜合考慮。

3.2 BV043/85規(guī)范計算

當(dāng)設(shè)備實(shí)際重量大于5t時，需考慮其受沖擊時對所處環(huán)境產(chǎn)生的反作用，必需對沖擊譜載荷進(jìn)行折減計算，折減僅針對譜加速度和譜速度，譜位移不用折減。對于100t設(shè)備的部分沖擊環(huán)境參數(shù)［11］，如表5所示。

表5 折減后BV043/85規(guī)范中的部分沖擊環(huán)境參數(shù)

說明：I區(qū)為船殼基座和主甲板以下隔艙壁；II區(qū)為下甲板和主甲板隔壁；III區(qū)為上層建筑。

3.3 毀傷等級理論計算

采用文獻(xiàn)中的理論計算方法［12］對平臺的毀傷等級進(jìn)行計算。水中爆炸造成的毀傷一般分為四級，見表6。

表6 毀傷的四個級別

毀傷因子和毀傷等級的關(guān)系曲線，如表7所示，表7中的毀傷等級與表6中毀傷等級相對應(yīng)，五級毀傷可當(dāng)作無毀傷或彈性變形狀態(tài)。

表7 水面艦艇的毀傷等級劃分

表8 浮動沖擊平臺毀傷因子

由表8計算結(jié)果可知，浮動沖擊平臺的毀傷因子在設(shè)計工況下的量值為0.029～0.049，對應(yīng)于水面艦艇毀傷等級為五級；對應(yīng)表6可知，中型浮動沖擊平臺為無毀傷或彈性變形狀態(tài)。

4 計算結(jié)果比對

4.1 毀傷因子比對分析

將表3仿真計算結(jié)果以及表8理論計算結(jié)果進(jìn)行比較，結(jié)果如表9所示。

由表可知，仿真計算與理論預(yù)報計算的仿真因子的誤差最大為3.48%，在工程許可范圍內(nèi)。證明了仿真計算方法可行。

4.2 沖擊環(huán)境計算結(jié)果比較

比對表4中的仿真計算結(jié)果以及表5中的BV043/85標(biāo)準(zhǔn)預(yù)估結(jié)果，可以得到：

1）取沖擊加速度參數(shù)進(jìn)行分析，對于重量小于5t的設(shè)備，中型浮動沖擊平臺垂向譜加速度為320g，沖擊加速度環(huán)境基本滿足使用要求，只是在垂直90°，沖擊因子為0.86工況下，沖擊加速度值偏大，試驗(yàn)實(shí)施過程中應(yīng)該采取一定的防護(hù)措施。

2）取沖擊速度參數(shù)進(jìn)行分析，對于重量小于5t的設(shè)備，中型浮動沖擊平臺譜垂向速度為7m/s，選取沖擊速度的情況下，各種工況均滿足要求。

由于大多數(shù)情況都是以沖擊速度值作為判據(jù)，在此對載重量為100t情況下的沖擊速度進(jìn)行分析，由表可中折減后的譜速度為2.11m/s，只有在角度很小的情況下滿足使用要求。仿真計算結(jié)果中最大的速度為4.44m/s，假定該值為折減后譜速度，根據(jù)折減公式可反推出此種沖擊環(huán)境下平臺的最大載重約為15.6t，此時最大載重量下降主要是由于爆源位置變化到平臺下方所致。

5 結(jié)語

通過本文研究，得到一些有意義的結(jié)論：

1）中型浮動沖擊平臺在標(biāo)準(zhǔn)工況下的毀傷等級為五級，可以看作為無毀傷或彈性變形狀態(tài)。

2）采用仿真計算及理論分析方法對毀傷因子的量值進(jìn)行計算并比對，兩種計算方法最大誤差在工程誤差范圍之內(nèi)，結(jié)果比較吻合。

3）對平臺底部沖擊環(huán)境進(jìn)行計算分析，沖擊加速度隨攻角增加而增大，在90°時達(dá)到最大值；沖擊因子相同情況下，沖擊速度及譜位移隨攻角變化不規(guī)律。

4）將仿真計算的沖擊環(huán)境數(shù)值與BV043/85標(biāo)準(zhǔn)量值進(jìn)行比對，探索了中型浮動沖擊平臺可承載設(shè)備的量值。平臺在進(jìn)行正式試驗(yàn)時其爆源布放位置對試驗(yàn)結(jié)果將有較大影響，應(yīng)嚴(yán)格按標(biāo)準(zhǔn)工況實(shí)施。