李思宇，李曉彬，趙鵬鐸，胡翔

1武漢理工大學(xué)交通學(xué)院，湖北武漢430063

2海軍裝備研究院，北京100073

近爆載荷作用下液艙的吸能研究

李思宇1，李曉彬1，趙鵬鐸2，胡翔1

1武漢理工大學(xué)交通學(xué)院，湖北武漢430063

2海軍裝備研究院，北京100073

為研究近距爆炸載荷作用下液艙各部分的吸能情況，根據(jù)實(shí)驗(yàn)建立數(shù)值仿真模型，研究在有無(wú)液體、不同厚度比和不同水層厚度條件下艙室變形和各部分吸能的占比情況。結(jié)果表明，液體介質(zhì)的存在改變了液艙的吸能模式，液艙總吸能主要受到艙室外壁厚度和水層厚度的影響，厚度比和水層厚度的變化對(duì)艙壁變形模式和爆炸能量在液艙各部分的占比有一定影響。對(duì)液艙各部分吸能機(jī)理的闡述，可作為液艙設(shè)計(jì)的參考依據(jù)。

近距爆炸；防護(hù)液艙；沖擊；吸能

0 引 言

大型水面艦船水下舷側(cè)常設(shè)有防護(hù)水艙，主要作用是針對(duì)魚雷和水雷的局部接觸、近距離非接觸爆炸，以及穿甲式、射流破甲式魚雷在艙室內(nèi)部爆炸的情況，保護(hù)內(nèi)部彈藥艙、動(dòng)力艙等重要艙室免遭破壞。開(kāi)展沖擊波作用下液艙的吸能研究，可為艦船防護(hù)液艙設(shè)計(jì)提供借鑒。

對(duì)液艙抗爆問(wèn)題的研究早在二戰(zhàn)前后就已經(jīng)開(kāi)始，但由于其軍事敏感性，國(guó)外公開(kāi)的文獻(xiàn)不多［1-2］。國(guó)內(nèi)對(duì)此類問(wèn)題的研究開(kāi)展較晚，朱錫等［3-4］對(duì)艦船水下多層防護(hù)結(jié)構(gòu)進(jìn)行了系列的水下接觸爆炸模型的實(shí)驗(yàn)研究。杜志鵬等［5］在適當(dāng)假設(shè)的基礎(chǔ)上，對(duì)液艙壁在近爆載荷作用下的響應(yīng)進(jìn)行了理論推導(dǎo)。張倫平等［6］通過(guò)模型實(shí)驗(yàn)，分析了爆炸載荷能量與結(jié)構(gòu)總吸能的比例關(guān)系，以及藥量、結(jié)構(gòu)參數(shù)對(duì)總吸能分配的影響。張婧等［7］通過(guò)數(shù)值仿真方法對(duì)三層板殼結(jié)構(gòu)在水下爆炸作用下的動(dòng)態(tài)響應(yīng)進(jìn)行了研究，并分析了液面高度的影響。唐廷等［8］采用流固耦合方法對(duì)實(shí)艙模型在水下接觸爆炸載荷作用下的變形和破壞進(jìn)行了研究。

為了研究液艙對(duì)爆炸沖擊波能量的吸收機(jī)理，基于能量守恒，擬對(duì)近爆載荷作用下，液艙吸收的爆炸能量進(jìn)行理論研究。在模型實(shí)驗(yàn)的基礎(chǔ)上，建立液艙在近爆載荷作用下的流固耦合模型，討論有無(wú)液體、不同厚度比和不同水層厚度條件下，液艙艙壁的變形情況和各部分吸能比例的變化，為液艙設(shè)計(jì)提供相關(guān)依據(jù)。

1 液艙吸收的爆炸能量

施加在液艙外壁上的近爆載荷可以表示為［9］

式中：pm為作用在外板的最大壓力；μ為裝藥形狀系數(shù)，球形為3，柱形為2，平面裝藥為1；w為瞬時(shí)爆轟平均壓力；Rw為裝藥半徑；R為炸藥中心與外板的垂直距離；α為爆炸產(chǎn)物噴流角；t為某一時(shí)刻；τ為正壓持續(xù)時(shí)間。

由式（1），根據(jù)沖量定理可得液艙外壁上任一點(diǎn)的比沖量

液艙外壁表面受到爆炸沖擊載荷作用時(shí)，其背面有水，需要考慮水體擾動(dòng)的影響。水體的影響采用等效質(zhì)量來(lái)描述，水的擾動(dòng)隨時(shí)間向縱深方向發(fā)展，因此等效質(zhì)量的增加也隨時(shí)間變化。假設(shè)板后任一質(zhì)點(diǎn)的水體與板同步運(yùn)動(dòng)，考慮厚度為c0τ的水的響應(yīng)，其中c0為水中聲速，則受擾動(dòng)質(zhì)點(diǎn)的質(zhì)量為根據(jù)動(dòng)量守恒原理得到外壁由反射比沖量引起的速度響應(yīng)為

假設(shè)作用在結(jié)構(gòu)表面的壓力被完全吸收且轉(zhuǎn)化為板和水體的初始動(dòng)能，則作用到結(jié)構(gòu)上的總能量可表示為

式中，r為載荷作用半徑，取積分上限l=1/2L，L為板的短邊長(zhǎng)度。式（4）表明，當(dāng)考慮板后水體的擾動(dòng)時(shí)，液艙的總體吸能減少了。將式（4）積分，可得液艙吸收的總能量的表達(dá)式為

2 仿真模型及方法驗(yàn)證

2.1 實(shí)驗(yàn)?zāi)Ｐ?/p>

實(shí)驗(yàn)?zāi)Ｐ褪浅叽鐬?00 mm×500 mm的Q235鋼板。鋼板通過(guò)螺栓和框架夾具面板固定在鋼制立方體爆炸筒的端面，立方體筒的尺寸為400 mm× 400 mm×200 mm。為了模擬艦船液艙，需要在筒內(nèi)注水，因此在鋼板和立方體筒端面之間附加一層橡膠墊，以防止液體泄漏。每個(gè)鋼板通過(guò)32個(gè)螺栓與立方體筒連接，保證四邊固支的邊界條件。夾具面板四邊預(yù)留一定寬度，用以固定鋼板，因此鋼板實(shí)際迎爆面積為400 mm×400 mm。如圖1所示，將液艙模型安置在事先做好的支架上，炸藥通過(guò)繩子固定在液艙外壁的中心位置，爆距為200 mm。炸藥類型選用TNT炸藥，藥量為200 g，板厚為2.5 mm。

圖1 實(shí)驗(yàn)?zāi)Ｐ徒Y(jié)構(gòu)和實(shí)體圖Fig.1 Test model structure and entity graph

2.2 仿真參數(shù)及方法驗(yàn)證

基于實(shí)驗(yàn)?zāi)Ｐ停捎梅蔷€性動(dòng)力分析軟件AUTODYN建立了1/2的對(duì)稱模型（圖2）。理想氣體和水域采用歐拉流體單元建立，液艙內(nèi)外壁用拉格朗日殼單元建立，采用全流固耦合算法計(jì)算。板采用四邊固支邊界條件，空氣域設(shè)置流出邊界條件。在外壁中心到邊界處以20 mm為間距設(shè)置11個(gè)測(cè)點(diǎn)（圖2）。艙壁材料采用Q235鋼，考慮應(yīng)變率強(qiáng)化效應(yīng)，材料的動(dòng)態(tài)流動(dòng)應(yīng)力采用J-C強(qiáng)度模型，具體材料參數(shù)如表1所示。其中，A，B，n，c，m為材料參數(shù)。

圖2 數(shù)值仿真模型Fig.2 Simulation model

表1 模型參數(shù)Table 1 Model parameters

圖3為實(shí)驗(yàn)板實(shí)際變形，圖4為仿真與實(shí)驗(yàn)得到的橫截面撓度曲線的對(duì)比。對(duì)比結(jié)果表明，仿真結(jié)果與實(shí)驗(yàn)曲線的一致性較好，仿真結(jié)果稍大于實(shí)驗(yàn)結(jié)果，最大撓度值誤差為5.8%。靠近邊界處的板撓度差異較大，這是由于實(shí)驗(yàn)的邊界條件不可能完全固支所致。說(shuō)明采用本文的仿真模型研究液艙的吸能問(wèn)題時(shí)，得到的結(jié)論是可靠的。

圖3 液艙外壁變形Fig.3 Deformation of front bulkhead

圖4 仿真與實(shí)驗(yàn)對(duì)比Fig.4 Contrast of simulation and experimental results

3 液艙吸能研究

3.1 有無(wú)液體的影響

圖5為計(jì)算至3 ms時(shí)，有無(wú)液體介質(zhì)液艙艙壁的變形模式對(duì)比?？梢钥闯觯瑹o(wú)液體時(shí)艙室外壁產(chǎn)生了整體塑性大變形，類似空背板的變形模式，而艙室內(nèi)壁變形很小，對(duì)吸能的貢獻(xiàn)很低；有液體時(shí)，艙室外壁僅發(fā)生局部的塑性變形，內(nèi)壁在水體的擾動(dòng)下產(chǎn)生了整體的塑性變形，且變形大于外壁，兩者總體變形量小于無(wú)液體艙室。此時(shí)液艙內(nèi)壁仍有速度，隨著動(dòng)能逐漸向變形能轉(zhuǎn)化，液艙內(nèi)壁最終也將發(fā)展為“中心鼓包”形式［10］。其原因?yàn)闊o(wú)液體介質(zhì)時(shí)，艙室主要靠艙壁的變形吸收爆炸能量，這種情況下，艙室外壁相當(dāng)于對(duì)沖擊載荷起到了“阻斷作用”，容易產(chǎn)生拉伸撕裂破壞［11］；而有液體介質(zhì)時(shí)，水的可壓縮性使得作用在一部分介質(zhì)上的壓力以高速而有限的波動(dòng)形式傳遞給其他質(zhì)點(diǎn)，導(dǎo)致大量能量被水體吸收，并將一部分能量傳遞給液艙內(nèi)壁。結(jié)果表明，液體介質(zhì)通過(guò)自身的吸能降低了艙壁的變形量，并改變了其變形模式，通過(guò)將沖擊載荷的吸收和彌散，將液艙內(nèi)外壁有機(jī)聯(lián)系起來(lái)。

圖5 有無(wú)液體介質(zhì)艙壁變形模式Fig.5 Deformation mode of bulkhead

圖6為有無(wú)液體介質(zhì)時(shí)液艙在近爆載荷作用下的能量轉(zhuǎn)化曲線?？梢钥闯?，無(wú)液體介質(zhì)時(shí)，爆炸能量主要被艙室外壁吸收轉(zhuǎn)化為其動(dòng)能和變形能；有液體介質(zhì)時(shí)，液艙整體的吸能降低了約60%，其中艙壁吸收的爆炸能量?jī)H為無(wú)液體介質(zhì)時(shí)的11%左右，有超過(guò)70%的能量被水體吸收，外壁僅吸收了不到4%的能量，內(nèi)壁約吸收了26%的能量。被吸收的能量大部分轉(zhuǎn)化為水體的動(dòng)能和勢(shì)能，部分轉(zhuǎn)化為艙壁的動(dòng)能和變形能。說(shuō)明艙內(nèi)液體改變了爆炸能量在液艙中的分散效果，由原先主要由外壁承擔(dān)變?yōu)閮?nèi)外壁和液體共同承擔(dān)，并通過(guò)三者之間的相互作用使大量能量通過(guò)水體的振動(dòng)和激蕩，最終以熱能的形式消耗。

圖6 能量轉(zhuǎn)化Fig.6 Energy transformation

3.2 厚度比的影響

文獻(xiàn)［5］和本文的研究都表明，在近爆載荷作用下，液艙內(nèi)壁的變形大于液艙外壁，應(yīng)考慮對(duì)內(nèi)壁進(jìn)行加厚處理。定義液艙內(nèi)壁和外壁厚度的比值為厚度比。壁通過(guò)形變達(dá)到吸能效果，增加內(nèi)壁的厚度勢(shì)必會(huì)影響到爆炸能量在液艙中的彌散效果，因此，需要對(duì)厚度比的影響進(jìn)行探究。圖7為計(jì)算至1 ms時(shí)，不同厚度比液艙艙壁變形的對(duì)比。可以看出，厚度比的變化并未改變艙壁的整體變形趨勢(shì)，雖然該時(shí)刻艙壁尚未達(dá)到最終變形，但隨著厚度比的提升，液艙內(nèi)壁的變形發(fā)展速度逐漸減緩；液艙外壁的總變形量也略有降低，但由于外壁厚度不變，因此對(duì)吸能的影響較小。

圖8為不同厚度比的液艙吸能對(duì)比?？梢钥闯?，兩者的變化規(guī)律基本一致。以厚度比為2∶3的能量曲線為例，在近爆載荷作用下，艙室外壁的能量首先增加，在0.1 ms時(shí)就達(dá)到峰值，約0.18 ms時(shí)外壁的動(dòng)能全部轉(zhuǎn)化，僅剩余變形能；水體中的能量在0.06 ms時(shí)開(kāi)始迅速增加，約在0.18 ms時(shí)到達(dá)峰值，此時(shí)沖擊載荷已到達(dá)液艙內(nèi)壁，內(nèi)壁的能量開(kāi)始上升；0.26 ms時(shí)，水體和內(nèi)壁之間的能量轉(zhuǎn)換結(jié)束，內(nèi)壁中能量基本穩(wěn)定，而水體中能量呈平穩(wěn)下降趨勢(shì)。厚度比為1∶1時(shí)，增加厚度比對(duì)液艙總的吸能和外壁吸能影響很小，但液艙內(nèi)壁吸能增加了約32%，當(dāng)厚度比增加至1∶2時(shí)，液艙內(nèi)壁吸能增加了約54%。結(jié)果表明，增加內(nèi)壁厚度有利于減小液艙內(nèi)壁的變形，同時(shí)增加內(nèi)壁的吸能效果，液艙設(shè)計(jì)時(shí)應(yīng)考慮一定的厚度比。

圖8 不同厚度比液艙吸能對(duì)比Fig.8 Contrast of energy absorption mode of liquid tank

3.3 水層厚度的影響

圖9為計(jì)算至1 ms時(shí)，不同水層厚度時(shí)液艙艙壁的變形對(duì)比?？梢钥闯觯畬雍穸葘?duì)艙壁變形產(chǎn)生了較大的影響。當(dāng)水層厚度由100 mm增加至300 mm時(shí)，液艙外壁變形量有先減小后增大的趨勢(shì)，其變形模式也發(fā)生了改變。當(dāng)水層較薄時(shí)，艙壁呈近似整體塑性大變形趨勢(shì)，隨著水層加厚，艙壁的中心隆起變形越來(lái)越明顯，說(shuō)明增加水層厚度使外壁的變形局部效應(yīng)更加突出。液艙內(nèi)壁在水層較薄時(shí)的變形迅速發(fā)展為中心隆起形式，水層厚度增加后，其變形發(fā)展趨勢(shì)變緩，但隨著水層厚度的繼續(xù)增加，其變形趨勢(shì)沒(méi)有明顯改變。水層較厚時(shí)，作用在液艙內(nèi)壁上的沖擊波近似為平面波，隨著傳播距離減小，載荷作用的局部效應(yīng)逐漸增強(qiáng)，導(dǎo)致液艙內(nèi)壁變形發(fā)展模式改變，進(jìn)而對(duì)其最終變形量產(chǎn)生影響。

圖9 不同水層厚度的液艙艙壁變形對(duì)比Fig.9 Contrast of deformation mode of bulkhead

圖10為不同水層厚度的液艙吸能對(duì)比。當(dāng)水層厚度增加后，水體與液艙內(nèi)壁的能量轉(zhuǎn)換過(guò)程被逐漸推移。當(dāng)水層厚度為100 mm時(shí)，液艙吸收的總吸能增加了約28.4%，其中液艙內(nèi)壁的吸能占比增加了約10%，水體的吸能占比有所下降；當(dāng)水層厚度增加至300 mm時(shí)，液艙的總吸能增加了約35.3%，而對(duì)艙內(nèi)各部分吸能占比影響較小，上下波動(dòng)不超過(guò)3%。薄水層和厚水層雖然都提高了液艙整體吸能，但其吸能增加的原因不同。水體通過(guò)質(zhì)點(diǎn)的運(yùn)動(dòng)吸收爆炸能量，水層厚度較小時(shí)水體擾動(dòng)較為劇烈，且液艙內(nèi)壁變形的迅速發(fā)展導(dǎo)致內(nèi)壁吸能增加。增加水層厚度實(shí)際上是增加了受擾動(dòng)的水體質(zhì)點(diǎn)數(shù)量，進(jìn)而提高了液艙的總吸能。式（5）表明，當(dāng)水層厚度大于c0τ時(shí)，多余的水體將不再對(duì)吸能有所影響。因此，在設(shè)計(jì)液艙時(shí)，應(yīng)綜合考慮依靠艙壁變形和水體擾動(dòng)的吸能效果，選取合適的水層厚度。

4 結(jié) 論

上述研究表明：

1）有液體介質(zhì)時(shí)，艙室吸收的爆炸能量相對(duì)于無(wú)液體介質(zhì)時(shí)有所降低；當(dāng)裝藥和爆距確定時(shí)，液艙吸收的爆炸能量主要和外壁厚度以及水體厚度有關(guān)。

圖10 不同水層厚度的液艙吸能對(duì)比Fig.10 Contrast of energy absorption mode of liquid tank

2）液體介質(zhì)通過(guò)自身的吸能降低了外壁的變形量，并改變了其變形模式；液體介質(zhì)改變了艙室的吸能模式，由單一的依靠外壁變形吸能轉(zhuǎn)變?yōu)橐揽績(jī)?nèi)外壁變形和液體質(zhì)點(diǎn)運(yùn)動(dòng)共同吸能。

3）增加內(nèi)壁厚度不影響液艙的總吸能和外壁的吸能效果，但對(duì)液艙內(nèi)壁的吸能有所提升，同時(shí)減小了其變形發(fā)展的速度，有利于保護(hù)艙室的內(nèi)部結(jié)構(gòu)。

4）水層厚度影響液艙總吸能和艙壁的變形模式，其厚度為c0τ時(shí)可保證液艙的最大吸能效果。水層較薄時(shí)，液艙通過(guò)水體的激蕩和內(nèi)壁的變形提高了整體吸能；水層較厚時(shí)，振動(dòng)水體質(zhì)點(diǎn)的增加是改變液艙吸能的原因。

［1］ 吉田隆.舊海軍艦船の爆彈被害損傷例いつぃて［J］.船の科學(xué)，1990，43（5）：69-73.

［2］ 吉田隆.舊海軍艦船の爆弾被害損傷例について--直撃爆弾および至近爆弾による損傷［J］.船の科學(xué)，1990，43（6）：62-65.

［3］ 朱錫，張振華，劉潤(rùn)泉，等.水面艦艇舷側(cè)防雷艙結(jié)構(gòu)模型抗爆試驗(yàn)研究［J］.爆炸與沖擊，2004，24（2）：133-139. ZHU X，ZHANG Z H，LIU R Q，et al.Experimental study on the explosion resistance of cabin near ship?board of surface warship subjected to underwater con?tact explosion［J］.Explosion and Shock Waves，2004，24（2）：133-139（in Chinese）.

［4］ 張振華，朱錫，黃玉盈，等.水面艦艇舷側(cè)防雷艙結(jié)構(gòu)水下抗爆防護(hù)機(jī)理研究［J］.船舶力學(xué)，2006，10（1）：113-119. ZHANG Z H，ZHU X，HUANG Y Y，et al.Theoreti?cal research on the defendence of cabin near shipboard of surface warship subjected to underwater contact ex?plosion［J］.Journal of Ship Mechanics，2006，10（1）：113-119（in Chinese）.

［5］ 杜志鵬，李曉彬，夏利娟，等.艦船防護(hù)水艙在接近爆炸載荷作用下響應(yīng)的理論研究［J］.船舶力學(xué)，2007，11（1）：119-127. DU Z P，LI X B，XIA L J，et al.Theory research on the response of the warship protective tank under near-by explosion［J］.Journal of Ship Mechanics，2007，11（1）：119-127（in Chinese）.

［6］ 張倫平，張曉陽(yáng)，潘建強(qiáng)，等.多艙防護(hù)結(jié)構(gòu)水下接觸爆炸吸能研究［J］.船舶力學(xué)，2011，15（8）：921-929. ZHANG L P，ZHANG X Y，PAN J Q，et al.Energy re?search about multicamerate defence structure subject?ed to underwater contact explosion［J］.Journal of Ship Mechanics，2011，15（8）：921-929（in Chinese）.

［7］ 張婧，施興華，王善.水下接觸爆炸作用下艦船防護(hù)結(jié)構(gòu)中液艙影響仿真分析［J］.天津大學(xué)學(xué)報(bào)，2008，41（10）：1238-1244. ZHANG J，SHI X H，WANG S.Numerical simulation analysis of liquid cabin of ship defensive structure sub?jected to underwater contact explosions［J］.Journal of Tianjin University，2008，41（10）：1238-1244（in Chinese）.

［8］ 唐廷，朱錫，侯海量，等.大型水面艦艇防雷艙結(jié)構(gòu)防護(hù)機(jī)理數(shù)值仿真［J］.哈爾濱工程大學(xué)學(xué)報(bào)，2012，33（2）：142-149. TANG T，ZHU X，HOU H L，et al.Numerical simula?tion study on the defense mechanism of a cabin near the shipboard for large surface vessels［J］.Journal of Harbin Engineering University，2012，33（2）：142-149（in Chinese）.

［9］ HENRYCH J.The dynamics of explosion and its use［M］.Amsterdam：Elsevier Scientific Pub.Co.，1979：218.

［10］ 徐雙喜.大型水面艦船舷側(cè)復(fù)合多層防護(hù)結(jié)構(gòu)研究［D］.武漢：武漢理工大學(xué)，2010. XU S X.Research on the composite and multilayer de?fensive broadside structure for large surface ship［D］. Wuhan：Wuhan University of Technology，2010（in Chinese）.

［11］ 陳長(zhǎng)海，朱錫，侯海量，等.近距空爆載荷作用下固支方板的變形及破壞模式［J］.爆炸與沖擊，2012，32（4）：368-375. CHEN C H，ZHU X，HOU H L，et al.Deformation and failure modes of clamped square plates under close-range air blast loads［J］.Explosion and Shock Waves，2012，32（4）：368-375（in Chinese）.

Research into energy absorption of liquid cabin subjected to close-range explosion

LI Siyu1，LI Xiaobin1，ZHAO Pengduo2，HU Xiang1

1 School of Transportation，Wuhan University of Technology，Wuhan 430063，China

2 Naval Academy of Armament，Beijing 100073，China

In order to study the energy absorption of different parts of a liquid cabin under a close-range explosion,a fluid-structure coupling model is built on the basis of experiments,and the deformation of the bulkhead and energy absorption ratio of different parts of the liquid cabin are analyzed,in which the influence of the water,bulkhead thickness ratio and water thickness are also discussed.The results show that the existence of a liquid medium can change the energy absorption model of a cabin.The total energy absorption is mainly affected by the front bulkhead thickness and water thickness,and alterations to the bulkhead thickness ratio or water thickness can also affect the deformation model of the bulkhead and energy absorption ratio of different parts of the cabin.A logical explanation of the energy absorption mechanisms of the liquid cabin is proposed,and some useful suggestions for designs are given.

close-range explosion；warship protective tank；shock；energy absorption

U661.4；O38

A

10.3969/j.issn.1673-3185.2017.01.015

2016-07-04

2016-12-28 15:35

國(guó)家自然科學(xué)基金資助項(xiàng)目（11302259）；非線性力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室開(kāi)放基金資助項(xiàng)目（LNM201505）

李思宇，男，1991年生，碩士生。研究方向：結(jié)構(gòu)安全與可靠性。E-mail：596540197@qq.com李曉彬（通信作者），男，1971年生，博士，教授。研究方向：結(jié)構(gòu)安全與可靠性。E-mail：lxbmark@163.com

http://www.cnki.net/kcms/detail/42.1755.TJ.20161228.1535.012.html期刊網(wǎng)址：www.ship-research.com

李思宇，李曉彬，趙鵬鐸，等.近爆載荷作用下液艙的吸能研究［J］.中國(guó)艦船研究，2017，12（1）：101-106，133. LI S Y，LI X B，ZHAO P D，et al.Research into energy absorption of liquid cabin subjected to close-range explosion［J］.Chinese Journal of Ship Research，2017，12（1）：101-106，133.