爆炸載荷下艙板與加強筋抗爆性能分析

梅宗書，石成英，蔡星會，王　爭

(第二炮兵工程大學(xué)，西安　710025)

?

爆炸載荷下艙板與加強筋抗爆性能分析

梅宗書，石成英，蔡星會，王爭

(第二炮兵工程大學(xué)，西安710025)

摘要：以救生艙殼體的抗爆性能為研究對象，通過控制變量法，分別研究了殼體中艙板與加強筋兩部分的抗爆性能。通過有限元分析軟件LS-DYNA，對艙板厚度一定，加強筋厚度不同和加強筋厚度一定，艙板厚度不同兩種組合進行數(shù)值仿真計算，得出了整體救生艙殼體的塑性變形與艙板和加強筋厚度之間的關(guān)系曲線。研究結(jié)論表明：在一定載荷范圍內(nèi)，加強筋在厚度與整體殼體塑性變形上成雙段線性曲線規(guī)律；艙板在厚度與整體殼體塑性變形上成拋物線曲線規(guī)律。研究成果為救生艙艙板與加強筋厚度設(shè)計提供有益借鑒。

關(guān)鍵詞：沖擊；LS-DYNA；艙板；加強筋；塑性變形

本文引用格式：梅宗書，石成英，蔡星會，等.爆炸載荷下艙板與加強筋抗爆性能分析[J].兵器裝備工程學(xué)報，2016(1):55-58.

Citation format:MEI Zong-shu, SHI Cheng-ying, CAI Xing-hui, et al.Antiknock Performance Analysis of Tank Plate and Stiffener Under Explosion Loading[J].Journal of Ordnance Equipment Engineering,2016(1):55-58.

作為煤礦、軍事等領(lǐng)域的人員防護裝備，救生艙的抗爆性能是研究設(shè)計其需要考慮的最重要的性能之一。救生艙根據(jù)功能有諸多分類，但在殼體結(jié)構(gòu)上大體相同，均為在救生艙艙板上焊接加強筋。以往在救生艙艙板和加強筋研究領(lǐng)域中，大多數(shù)為對復(fù)合材料加強筋的抗爆性能和不同材質(zhì)、夾層結(jié)構(gòu)的艙板抗爆性能的研究[1.2]。我國學(xué)者丁曉紅、李國杰等通過對自然界分支系統(tǒng)形態(tài)的最優(yōu)性研究[3]，提出一種能使加強筋沿結(jié)構(gòu)整體剛度最大方向生長的自適應(yīng)成長法。但是將艙壁與加強筋放在一起，對其在厚度維度上的抗爆性能進行研究分析，公開文獻還比較少。

目前，在救生艙的抗爆性能研究方面，對于艙板和加強筋厚度的設(shè)計，通常做法是艙板厚度由材質(zhì)和抗爆性能指標(biāo)確定，加強筋厚度參照相關(guān)行業(yè)標(biāo)準(zhǔn)，一般不考慮艙板厚度與加強筋尺寸方面的相互聯(lián)系。本研究則分別探究兩者厚度與救生艙殼體抗爆性能之間的關(guān)系，從而為艙板與加強筋厚度的確定提供依據(jù)。

本研究通過有限元分析軟件LS-DYNA，對艙板厚度一定，加強筋厚度不同和加強筋厚度一定，艙板厚度不同兩種組合下，救生艙殼體的最大塑性變形進行研究，得到相關(guān)線性曲線，為艙板與加強筋的厚度設(shè)計提供有益參考。

1模型與參數(shù)設(shè)置

本研究對救生艙殼體進行建模，殼體由高強度鋼質(zhì)大板和內(nèi)附加強筋組成，艙板與加強筋材質(zhì)相同，材料模型采用塑性隨動模型*MAT_PLASTIC_KINEMATIC，單位制為cm-g-us，參數(shù)設(shè)置如表1所示。

表1　參數(shù)設(shè)置

模型采用Shell163薄殼單元，加強筋截面形狀設(shè)計為半圓形[4]，圓面半徑為6 cm，通過在*SECTION_SHELL中輸入厚度參數(shù)設(shè)置加強筋與艙板厚度。由關(guān)鍵字*LOSD_BLAST施加爆炸載荷，爆炸點設(shè)置在相對艙板對稱的位置，沖擊載荷大小為0.8 MPa，模型尺寸如圖1所示。

圖1　艙板示意圖

艙板在某一時刻的應(yīng)力云圖如圖2所示。

圖2　艙板應(yīng)力原圖

由圖2可知，在救生艙殼體的中間位置，加強筋承受的應(yīng)力值最大，塑性變形也最為嚴(yán)重，為研究殼體最大形變量，本研究選取艙板中心位置的節(jié)點進行殼體法線方向的位移取樣[5]，并進行數(shù)據(jù)分析。

本研究采用控制變量法，將實驗工況分為兩組，具體參數(shù)設(shè)置如表2所示。通過對表2中兩種組合進行計算，得到艙板與加強筋各自的抗爆性能。

表2　實驗工況參數(shù)

2結(jié)果分析

通過計算，并在后處理器LS-PREPOST中選取艙板中心位置節(jié)點為研究對象，輸出兩種組合下該節(jié)點隨時間的位移變化曲線，如圖3所示。

由圖3可知，沖擊載荷作用下，在艙板的動力響應(yīng)初始階段，不同厚度的艙板與加強筋，救生艙殼體在法線方向的變形同步性較好。通過比較曲線最低點對應(yīng)形變值可知，組合2比組合1在變形幅度上要?。挥钟汕€之間的疏密程度可知，在增加相同厚度時，組合2比組合1應(yīng)變值增加輻度大。隨著曲線波動區(qū)域平緩，最后圍繞一條中心線上下波動，如圖4所示。

圖3　艙板中心節(jié)點初始位移時程曲線

圖4　艙板中心點完整位移時程曲線

由圖4中曲線的變化規(guī)律可以預(yù)測，假定時間足夠長，則所有曲線都將趨于一條縱坐標(biāo)不等的水平線，稱該水平線為塑性變形線[6]，如圖5所示。本研究以塑性變形線所對應(yīng)的縱坐標(biāo)大小作為救生艙殼體抗爆性能的指標(biāo)，將計算出的中心節(jié)點位移的平均值作為殼體塑性變形量的大小。

圖5　加強筋為1.4 cm時變形曲線

根據(jù)理論推導(dǎo)可知，變形曲線中的直線即為艙體最大塑性變形曲線，縱坐標(biāo)值即為塑性變形值——即撓度的大小，不同厚度下，艙板最大撓度值如圖6所示。

圖6　艙體塑性變形隨結(jié)構(gòu)厚度的變化曲線

由圖6可知，在厚度越大，艙板塑性變形量越小的同時，厚度與塑性變形并不滿足一階線性關(guān)系。兩條曲線中組合1中的曲線可看成分段雙直線，組合2表現(xiàn)為拋物線的左支。進一步計算每一段的斜率可得表3所示的參數(shù)。

由表3可知，加強筋在厚度遞增時，救生艙殼體撓度呈現(xiàn)分段式遞減規(guī)律，即對于加強筋而言其抗爆性能呈現(xiàn)分段式遞減規(guī)律，而對于艙板而言，當(dāng)厚度增加到一定厚度時，繼續(xù)增加其厚度，對抗爆性能的提高效果并不明顯。同時，從1.7 cm厚度開始組合2中的斜率下降幅度較大，基本和組合1中的斜率在一個區(qū)間，因此，此時通過增加加強筋厚度來提高救生艙殼體的抗爆性能是比較合適的。依據(jù)本文中的模型與材料設(shè)置，本文中艙板與加強筋的最佳的厚度組合為1.6+1.5 cm。

表3　分段斜率參數(shù)

3結(jié)論

本研究通過分別對救生艙殼體中的艙體大板和加強筋的抗爆性能進行模擬實驗和數(shù)據(jù)分析，以曲線圖的方式給出了加強筋和艙板在抗爆性能上的不同。

當(dāng)艙板厚度達到一定值時，繼續(xù)增加其厚度并不能較有效地提高救生艙的抗爆性能，同時，增加艙板厚度會大幅增加艙體的總重和制造成本，此時因結(jié)合救生艙抗爆性能指標(biāo)，選擇艙板與加強筋的塑性應(yīng)變與厚度之比較大時各自對應(yīng)的厚度，這樣才能最大程度縮減材料整體厚度和自重，有效節(jié)約生產(chǎn)成本。

參考文獻：

[1]項晨,陳浩然,郭兆璞.復(fù)合材料加筋結(jié)構(gòu)的流固耦合振動及動力響應(yīng)分析[J].復(fù)合材料學(xué)報,1996,13(3):100-110.

[2]肖鋒,諶勇,章振華,等.夾層結(jié)構(gòu)沖擊動力學(xué)研究綜述[J].振動與沖擊,2013,32(18):1-7.

[3]丁曉紅,李國杰,蔡戈堅,等.薄板結(jié)構(gòu)的加強筋自適應(yīng)成長設(shè)計法[J].中國機械工程,2005,16(12):1057-1060.

[4]平玉城,金龍哲,栗婧,等.不同加強筋對救生艙艙體抗爆性能影響的研究[J].煤,2014,23(10):4-7.

[5]MAIN J A,GAZONAS G A.Uniaxial crushing of sandwich plates under air blast:Influence of mass distribution[J].International Journal of Solids and Structures,2008,45(7/8):2297-2321.

[6]張旭紅,王志華,趙隆茂.爆炸載荷作用下鋁蜂窩夾芯板的動力響應(yīng)[J].爆炸與沖擊,2009,29(4):356-360.

(責(zé)任編輯周江川)

石成英(1964—)，男，教授，博士生導(dǎo)師，主要從事戰(zhàn)斗部工程研究。

【裝備理論與裝備技術(shù)】

Antiknock Performance Analysis of Tank Plate and Stiffener

Under Explosion Loading

MEI Zong-shu, SHI Cheng-ying, CAI Xing-hui, WANG Zheng

(The Second Artillery Engineering University, Xi’an 710025, China)

Abstract:With the capsule shell antiknock performance as the research object, two parts of the capsule shell, plate and stiffener, whose antiknock performance were researched through the control variable method. Through finite element analysis software LS-DYAN, two kinds of operating condition, that the thickness of the stiffener is different with a certain thickness of plate and that the thickness of the plate is different with a certain thickness of stiffener, were under numerical simulation calculation, which discovered the relationship between the plastic deformation of the whole capsule shell and the thickness of the plate and the stiffener. The result show that, within a certain range of load, the plastic deformation of the overall shell and the thickness of stiffener forms into double linear curve law; the plastic deformation of the overall shell and the thickness of plate forms into a parabolic curve. This research provides beneficial reference for the thickness design of capsule plate and stiffener.

Key words:impact; LS-DYNA; deck stiffener; plastic deformation

文章編號：1006-0707(2016)01-0055-04

中圖分類號：TD774

文獻標(biāo)識碼：A

doi:10.11809/scbgxb2016.01.013

作者簡介：梅宗書(1991—)，男，碩士研究生，主要從事兵器科學(xué)與技術(shù)研究；

收稿日期：2015-06-03；修回日期：2015-06-25