程曉達 曹冬梅 尹緒超 陳海龍

（1.哈爾濱工程大學(xué) 船舶工程學(xué)院 哈爾濱 150001；2.哈爾濱工業(yè)大學(xué) 航天學(xué)院 哈爾濱 150001）

0 引 言

近年來，隨著軍事競爭對海上力量的日益傾斜，水中兵器威力不斷增加，爆炸沖擊波及氣泡脈動、水射流都會對艦船生命力產(chǎn)生巨大的威脅［1-2］。各軍事強國對艦船遭受武器攻擊的生命力問題均開展了深入的探討［3］，我國在該領(lǐng)域的研究也在積極開展。為保持艦船的持續(xù)作戰(zhàn)能力，軍需補給系統(tǒng)的作用不容忽視。然而當(dāng)前國內(nèi)對補給艦艇的抗沖擊問題研究甚少，對于遭受水下武器攻擊時補給艦的強度問題、局部易損區(qū)域、安全性校核等問題均有待探討。

實船水下爆炸實驗費用驚人，目前一些西方發(fā)達國家也往往需要聯(lián)合幾個國家共同實施實船爆炸試驗，且涉及保密，公開的實船試驗數(shù)據(jù)非常稀少。我國在實船試驗方面起步相對較晚，但是近年來也加大了這方面的投入［4］。隨著科技的進步，數(shù)值計算的優(yōu)勢日漸顯著，具有方便、可靠、經(jīng)濟等特點，目前得到了廣泛應(yīng)用［5-7］。本文針對某典型軍用補給油船，通過數(shù)值計算與理論相結(jié)合的方法，對該船的中遠(yuǎn)場抗爆抗沖擊性能進行了分析計算。

ABAQUS由于其優(yōu)異的非線性分析能力，在結(jié)構(gòu)強度領(lǐng)域被廣泛應(yīng)用。本文使用該軟件集成的聲固耦合法可以很好地模擬中遠(yuǎn)場水下爆炸時沖擊波對船體結(jié)構(gòu)的作用，及船體結(jié)構(gòu)與水的相互耦合作用。由于接觸爆炸是個極其復(fù)雜的物理化學(xué)過程，至今仍是各國專家研究的前沿?zé)狳c問題。隨著數(shù)值仿真技術(shù)的不斷發(fā)展，各數(shù)值計算軟件百花爭鳴，互有優(yōu)勢，因而本文不再詳述。

1 有限元模型建立及工況設(shè)置

1.1 船體結(jié)構(gòu)及其流場有限元模型

本文所研究的是某典型補給油船結(jié)構(gòu)，在進行結(jié)構(gòu)抗沖擊分析前，按照設(shè)計圖紙采用通用商業(yè)軟件ANSYS對該船進行建模。主要考察該船抗沖擊性能，為保證模型的計算精度，在建模中全船所有基本結(jié)構(gòu)均詳細(xì)建出；同時，對于便于加工而設(shè)計的倒角、圓孔等結(jié)構(gòu)進行的合理簡化，以提高計算效率。最終建立的油船幾何模型如圖1所示。

對建立的模型進行網(wǎng)格劃分，并對局部結(jié)構(gòu)復(fù)雜以及重點考察區(qū)域的網(wǎng)格進行了細(xì)化處理后，通過實驗室自主開發(fā)的接口程序?qū)⒂邢拊Ｐ蛯?dǎo)入到ABAQUS中進行分析。

為了真實反映艦船在海上的航行狀態(tài)，通過查閱設(shè)計圖紙，將船體內(nèi)部設(shè)備及裝載以MASS質(zhì)量點形式均布在各層甲板，使艦艇排水量達到設(shè)計排水量，并調(diào)整重心位置與設(shè)計值一致。根據(jù)相關(guān)資料［8］，流場的范圍在船體半寬4倍左右可以在精度和計算效率上取得很好的平衡。在導(dǎo)入有限元模型后，在ABAQUS中建立流場模型，并與船體裝配在一起。在船體外表面和流場表面采用TIE約束模擬爆炸后水與結(jié)構(gòu)的耦合作用。船體與流場有限元模型如圖2所示。其中，全船節(jié)點總數(shù)為107746，有限元單元總數(shù)為264040。

1.2 工況設(shè)置

為綜合考慮藥量、爆距等因素的影響，本文采用沖擊因子C描述非接觸爆炸工況［9］，即

式中:W為藥包質(zhì)量（等效TNT質(zhì)量）；

R為爆心位置到艦船結(jié)構(gòu)的最短距離。

在相同的沖擊因子下，艦船結(jié)構(gòu)沖擊響應(yīng)可近似相等。本文擬從三個方面綜合校核該船體結(jié)構(gòu)的抗沖擊性能，即:藥包沿船長方向分布、藥包沿水深分布和不同攻角θ對船體結(jié)構(gòu)的影響。

（1）工況1～工況3為藥包沿船長方向分布，取沖擊因子為0.53，藥包位置分別在船艏處、船舯處和船艉處三種情況，重點考核船舶結(jié)構(gòu)響應(yīng)的變化。

（3）工況4～工況6為不同爆距對船體結(jié)構(gòu)沖擊響應(yīng)的影響，爆點位置均在船舯正下方，分別取沖擊因子0.58、0.63、0.70進行對比分析。工況7～工況9為不同攻角下非接觸爆炸對船體結(jié)構(gòu)的響應(yīng)。爆點位置在船舯正下方，沖擊因子0.58，攻角分別為30°、60°、90°。 藥包布置方式如圖 3 所示:

2 計算結(jié)果分析

2.1 藥包沿船長方向分布

工況1～工況3主要考察爆心在船艏、船舯和船艉三個典型位置時，船體結(jié)構(gòu)的響應(yīng)沿船長的分布規(guī)律，以及在上述工況中船體結(jié)構(gòu)抗沖擊性能，評估其薄弱環(huán)節(jié)。下面給出三種工況下船體結(jié)構(gòu)的應(yīng)力云圖，如圖4所示。

由云圖可以看出，當(dāng)爆炸沖擊波到達結(jié)構(gòu)表面時，由于藥包位置不同，船體結(jié)構(gòu)靠近藥包位置的局部區(qū)域應(yīng)力最大，并沿船長方向向其他結(jié)構(gòu)傳遞，船體結(jié)構(gòu)整體效應(yīng)明顯。為了更直觀的研究船體響應(yīng)沿船長方向分布，將該船沿船長方向平均分成20個區(qū)域，分別提取各分段的無量綱譜速度和無量綱譜加速度，通過樣條曲線繪制不同藥包位置船體沿船長方向無量綱譜速度與譜加速度分布曲線，如圖5、圖6所示。

由圖5可以看出，在爆距為60 m的遠(yuǎn)場爆炸載荷下，隨著藥包沿船長位置的不同，船體結(jié)構(gòu)的無量綱譜速度整體分布均勻，在藥包所在位置響應(yīng)略微偏高，并沒有陡峭的峰值。說明在武器攻擊引起的遠(yuǎn)場爆炸下，艦船主要體現(xiàn)整體響應(yīng)，總體效應(yīng)明顯。圖6中無量綱加速度曲線變化也較平緩，沒有出現(xiàn)劇烈的峰值，也很好地證明了這一點。

2.2 油船結(jié)構(gòu)響應(yīng)沿垂向分布

工況4～工況6主要考察不同爆距對油船結(jié)構(gòu)的沖擊響應(yīng)的影響，由于藥包均布置在船舯正下方，為不失一般性，分別選取船舯部位外底板及主甲板典型部位的垂向加速度時歷曲線，如圖7所示。

由圖7可見，在水下非接觸爆炸發(fā)生后，船體結(jié)構(gòu)的外底板結(jié)構(gòu)加速度快速增大，在0.01～0.02 s內(nèi)達到峰值，這是由于外底板是船體最靠近爆心位置的結(jié)構(gòu)，沖擊波最先到達，隨著爆距的減小，其沖擊波作用產(chǎn)生加速度峰值逐漸遞增。

對比相同工況的主甲板加速度曲線，沖擊波最先到達的外底板結(jié)構(gòu)加速度曲線波形陡峭，包含大量高頻毛刺，這一般是由于沖擊波的前驅(qū)波所致。達到峰值后，外底板結(jié)構(gòu)的加速度迅速減小，并在0.2 s內(nèi)趨于平緩，這是由于船體外板在遭受沖擊波作用時與水產(chǎn)生劇烈的耦合作用，其能量很快損耗，使外板加速度很快趨近于0。而主甲板在空氣中自由振蕩，其加速度減緩的速度明顯低于外底板，在爆炸發(fā)生后0.2 s，其加速度仍在0附近上下振蕩。

沖擊波在船體結(jié)構(gòu)傳遞中能量不斷衰減，通過主甲板加速度曲線可以直觀地觀察到其主甲板加速度峰值要大幅低于外底板，且具有明顯的低頻成分。這是船體結(jié)構(gòu)的低通濾波效應(yīng)引起的，低頻成分來自沖擊波激起的船體內(nèi)部結(jié)構(gòu)的固有頻率振動。

2.3 油船結(jié)構(gòu)響應(yīng)隨攻角分布

油船在遭受非接觸爆炸攻擊時，為考察相同爆距相同藥量下不同攻角對結(jié)構(gòu)的損傷效應(yīng)，定義無量綱攻擊角度為θ/（π/2），其中θ為爆炸攻擊角度，分別取 30°、60°、90°三個典型角度。 工況 7～工況 9對船體結(jié)構(gòu)響應(yīng)受攻角影響進行探究，其無量綱譜速度V及譜加速度A隨無量綱攻擊角的變化如圖8～圖9所示。

由圖可以看出，隨著攻角的增大，船舯部位的無量綱譜速度V逐漸增大，無量綱譜加速度A也略有增加，但整體趨勢平緩，沒有較大的變化，可以認(rèn)為攻角的變化對油船結(jié)構(gòu)沖擊響應(yīng)影響很小。這與中遠(yuǎn)場爆炸艦船局部效應(yīng)不明顯，整體效應(yīng)占主要因素有關(guān)。

3 結(jié) 論

本文采用ANSYS對某油船結(jié)構(gòu)進行有限元建模，通過ABAQUS對油船模型在遭受水下武器攻擊時的非接觸爆炸載荷進行模擬計算，對不同藥包位置不同藥量對船體結(jié)構(gòu)的沖擊響應(yīng)的影響進行分析、總結(jié)規(guī)律，得出以下結(jié)論:

（1）武器攻擊造成的中遠(yuǎn)場水下非接觸爆炸爆心位置對船體結(jié)構(gòu)沖擊響應(yīng)的分布影響很小，船上各部分的響應(yīng)趨于均勻化。油船的沖擊響應(yīng)主要表現(xiàn)為總體響應(yīng)，其局部效應(yīng)并不明顯。

（2）爆距對船體結(jié)構(gòu)的響應(yīng)有直接影響，隨著爆距的減小，艦船結(jié)構(gòu)響應(yīng)逐漸增大，且外底板最靠近爆點，受到非接觸爆炸的威脅最大，需適當(dāng)加強。

（3）沖擊波在船體結(jié)構(gòu)中的傳遞會損耗大量能量，主甲板結(jié)構(gòu)響應(yīng)遠(yuǎn)小于外底板。合理的總體區(qū)劃和艙室布置可以大大提高艦船設(shè)備及人員的安全。

（4）遠(yuǎn)場爆炸時沖擊波的球面效應(yīng)不明顯，其爆炸攻角對結(jié)構(gòu)響應(yīng)的影響不大。

［1］張阿漫，王詩平，白兆宏，等.不同環(huán)境下氣泡脈動特性實驗研究［J］.力學(xué)學(xué)報， 2011， 43（1）:71-83.

［2］張阿漫，姚熊亮.近自由面水下爆炸氣泡的運動規(guī)律研究［J］.物理學(xué)報， 2008， 57（1）:339-353.

［3］Said Michaels.Theory and Practice of Total Ship Survivability Analysis and Design［J］.Naval Engineers Journal，1995，107（4）:191-203.

［4］Calvano Charles N，etc.Establishing the Fundamentals of a Surface Ships Survivability Design Discipline［J］.Naval Engineers Journal，1994，106（1）:71-74.

［5］王官祥，趙玫.用 ANSYS進行抗沖擊仿真［M］//汪玉.艦艇及設(shè)備沖擊響應(yīng)分析技術(shù).北京:海潮出版社，2005:179-181.

［6］王官祥，汪玉.模態(tài)分析在沖擊動力學(xué)分析中的應(yīng)用［M］//汪玉.艦艇及設(shè)備沖擊響應(yīng)分析技術(shù).北京:海潮出版社， 2005:171-174.

［7］吳廣明，沈榮瀛，華宏星.大型艦船模型三維有限元振動及沖擊計算［M］//汪玉.艦艇及設(shè)備沖擊響應(yīng)分析技術(shù).北京:海潮出版社，2005:102-107.

［8］姚熊亮，張阿漫，許維軍.聲固耦合方法在艦船水下爆炸中的應(yīng)用［J］. 哈爾濱工程大學(xué)學(xué)報.2005，26（6）:707-712.

［9］庫爾.水下爆炸［M］.羅耀杰等譯.北京:國防工業(yè)出版社.1960:112-116.