殼單元和實體單元模擬爆炸荷載作用鋼板動力響應(yīng)的差異性比較

朱俊杰，盧永剛，馮星奎，陶俊林

(1. 西南科技大學(xué) 土木工程與建筑學(xué)院，四川 綿陽 621010；2. 中國工程物理研究院 總體工程研究所，四川 綿陽 621000；3. 火箭軍駐中國工程物理研究院 型號辦事處，四川 綿陽 621900)

殼單元和實體單元模擬爆炸荷載作用鋼板動力響應(yīng)的差異性比較

朱俊杰1，盧永剛2，馮星奎3，陶俊林1

(1. 西南科技大學(xué) 土木工程與建筑學(xué)院，四川 綿陽 621010；2. 中國工程物理研究院 總體工程研究所，四川 綿陽 621000；3. 火箭軍駐中國工程物理研究院 型號辦事處，四川 綿陽 621900)

以艦船艙室鋼板為研究對象，利用動力分析有限元程序 Ls-dyan，分別采用殼單元和實體單元 2 種建模方式對爆炸荷載作用下鋼板的動力響應(yīng)進行數(shù)值模擬研究。對比鋼板的應(yīng)力分布云圖、中線上各點最大位移曲線圖，發(fā)現(xiàn) 2 種方法計算結(jié)果與試驗值大致吻合，采用殼單元建模在保證計算結(jié)果準確性的前提下可以大大提高計算效率。如何正確地選擇單元進行分析是數(shù)值模擬很重要的一個環(huán)節(jié)，通過對殼單元和實體單元計算結(jié)果的初步對比，為艦船艙室數(shù)值模擬研究如何選取單元提供參考。

船用鋼板；爆炸荷載；數(shù)值模擬；殼單元

0 引 言

艦艇是由縱橫構(gòu)架和殼板結(jié)構(gòu)組成，其中殼板主要承受水平荷載。艦船艙室結(jié)構(gòu)承受的主要載荷之一是戰(zhàn)斗部爆炸產(chǎn)生的沖擊波載荷[1]，鋼板在受到具有傳播速度快、峰值荷載大、作用時間短等特征的爆炸沖擊波載荷作用時，其受力特性和動態(tài)響應(yīng)顯著異于靜載情況下，且十分復(fù)雜，故而研究鋼板在爆炸荷載作用下的動力響應(yīng)很有必要，是研究艙室結(jié)構(gòu)的基礎(chǔ)。目前國內(nèi)對艦船艙室方面的研究大多集中在艙室結(jié)構(gòu)的仿真分析上，朱建方等[2]對艦艇艙室在爆炸荷載下的動態(tài)響應(yīng)過程進行了數(shù)值模擬，得到了艦艇艙室在爆炸沖擊波作用下的響應(yīng)特性和破壞機理。鄢順偉等[3]仿真計算了不同壁厚艙室結(jié)構(gòu)在反艦導(dǎo)彈戰(zhàn)斗部內(nèi)爆作用下的毀傷效應(yīng)。但是，對船用鋼板在爆炸荷載作用下的響應(yīng)研究較少。本文采用數(shù)值模擬的方法對爆炸載荷下四邊固支的鋼板進行研究。在利用數(shù)值分析研究時，單元類型的選取對計算效率有很大的影響。鋼板屬于板結(jié)構(gòu)，實體單元與殼單元均可以對其進行模擬，但是 2 種建模方法得到結(jié)果的差異值得探討。通過分析兩者的差別，可以更好地為鋼板使用何種單元在保證其結(jié)果準確性的基礎(chǔ)上又能提高計算效率提供參考。為此，本文以文獻中試驗為基礎(chǔ)[4–5]，在驗證數(shù)值計算方法可靠的基礎(chǔ)上，分別采用殼單元與實體單元對四邊固支的方形板在爆炸荷載作用下的動力響應(yīng)進行數(shù)值模擬研究，考察使用不同的建模單元類型對計算結(jié)果的影響[6]，并對結(jié)果的差異性進行了比較與分析。

1 Shell163 單元介紹

Shell163 是一個具有 4 節(jié)點單元，有彎曲和膜特征，可施加平面和法向載荷。單元在每個節(jié)點上有 12個自由度，分別為在節(jié)點 x，y 和 z 方向的平動、加速度、速度和繞 x，y 和 z 軸的轉(zhuǎn)動。Shell163 單元一共有 12 種算法，這里選第 10 種算法，對于翹曲情況可以得到正確結(jié)果。和實體單元一樣，積分點的個數(shù)直接影響著 CPU 的時間。典型地，對于彈性材料沿厚度方向需要 2 個積分點，而對于塑性材料則需要 3 個或更多的積分點[7]。

2 數(shù)值計算模型

2.1 算法選取

分析主要涉及到流體流動、固體大變形以及流體與固體相互作用等問題，Ls-dyna 程序中提供了 Lagrange，Euler 和 ALE（任意拉格郎日歐拉算法）3 種不同的算法，其中 ALE 算法兼具 Lagrange 算法和 Euler算法的優(yōu)點，可以描述流體與固體的相互耦合作用；同時在求解過程中網(wǎng)格還可以根據(jù)定義的參數(shù)適當調(diào)整位置，克服網(wǎng)格嚴重畸變的問題[8]。因此，本文中選取 Ls-dyna 程序提供的 ALE 算法建立流體與固體間的聯(lián)系，實現(xiàn)流固耦合動態(tài)分析，模擬方形鋼板在爆炸荷載作用下的響應(yīng)和破壞。

2.2 材料本構(gòu)與狀態(tài)方程

結(jié)構(gòu)材料在爆炸沖擊載荷作用下會產(chǎn)生巨大的變形，與受靜載荷作用的情況相比，動態(tài)屈服應(yīng)力比靜態(tài)屈服應(yīng)力有明顯提高，材料力學(xué)性能表現(xiàn)相對也更為復(fù)雜。固支方板采用低碳鋼，鋼板選用與應(yīng)變率相關(guān)的各向同性塑性隨動硬化模型，材料參數(shù)見表 1。

表 1 鋼板材料參數(shù)Tab. 1 Steel plate material parameters

炸藥在空氣中爆炸時，會在極短的時間內(nèi)把化學(xué)能全部轉(zhuǎn)變?yōu)闊崮?，?dǎo)致爆轟產(chǎn)物向周圍迅速膨脹而做功，形成空氣沖擊波。文中采用 *HIGH_EXPLOSIVE_BURN 模型作為炸藥的材料模型，采用 Jones-Wilkins-Lee（JWL）狀態(tài)方程來描述，即

式中：A ，B，R1，R2和 ω 為 JWL 狀態(tài)方程參數(shù)；E為炸藥的內(nèi)能；V 為當前相對體積。

對空氣采用 *MAT_NULL 材料模型以及線性多項式 *EOS_LINEAR_POLYNOMI 狀態(tài)方程加以描述，即

式中：C0～C6為常量；e 為初始體積內(nèi)能；μ = 1/V0–1；V0為相對體積；C1～C4均為 0。

2.3 有限元模型

1） 空氣網(wǎng)格的驗證

在數(shù)值模擬中，網(wǎng)格的劃分至關(guān)重要，網(wǎng)格的大小不僅關(guān)系著數(shù)值模擬計算的精度，也關(guān)系著數(shù)值模擬計算所用的時間。一般而言，網(wǎng)格越小則計算的結(jié)果越接近于真實情況，但隨著網(wǎng)格的減少計算時間也急劇增大，所以選擇合理的網(wǎng)格尺寸很有必要。網(wǎng)格驗證幾何模型尺寸為空氣 120 mm × 120 mm × 60 mm，炸藥為 27 mm × 27 mm × 27 mm。利用對稱性取 1/4 模型進行構(gòu)建，模型如圖 2 所示。

表 2 炸藥與狀態(tài)方程參數(shù)Tab. 2 Explosive and equation of state parameters

表 3 多項式狀態(tài)方程參數(shù)Tab. 3 Polynomial equation of state parameters

空氣與炸藥均采用 8 節(jié)點實體單元，TNT 當量為 20 g，炸藥中心位置起爆。在對稱面上施加對稱約束，四周（除對稱面）施加固定邊界條件。為了模擬出無限空氣域，空氣模型邊界施加透射邊界。工況如表 4 所示。

表 4 空氣網(wǎng)格驗證工況表Tab. 4 Air grid authentication working table

上述工況在距離炸藥中心 6 cm 和 8 cm 位置處的超壓如圖 3～圖 4 所示。

由圖 5～圖6 可以看出，隨著網(wǎng)格的減少，壓力的峰值越來越大，當網(wǎng)格取 0.3 cm 和 0.2 cm 時壓力峰值的差值約為 5%，考慮到計算精度和時間的因素空氣網(wǎng)格取 0.3 cm。

2） 炸藥網(wǎng)格的驗證

根據(jù)上面網(wǎng)格確定炸藥網(wǎng)格的驗證工況如表 5 所示。

表 5 炸藥網(wǎng)格驗證工況表Tab. 5 Explosive grid authentication working table

由圖 5～圖 6 可知，隨著炸藥網(wǎng)格的減少，壓力峰值有所增加，當網(wǎng)格取為 0.2 cm 和 0.1 cm 時的壓力峰值差值約為 2%，炸藥網(wǎng)格選取 0.2 cm 較為合適。

3） 2 種單元計算模型

采用 Ansys/Ls-dyna 建立有限元模型，對于固支鋼板在爆炸載荷的作用過程進行了數(shù)值模擬研究。固支鋼板尺寸為 63 mm × 63 mm，厚度為 1.6 mm，建立 1/4模型。方板尺寸如圖 7 所示。

數(shù)值模型由炸藥、空氣、鋼板 3 部分組成，采用殼單元和實體單元分別建立板模型，2 種模擬方式除了方板的單元類型不同，其他參數(shù)均一致?？諝鈫卧叽鐬?3 mm，炸藥單元尺寸為 2 mm，板單元為 0.4 mm，有限元模型如圖 8～圖9 所示。

2.4 失效準則

采用 *MAT_ADD_EROSION 對鋼板材料失效準則進行定義，塑性失效應(yīng)變 εmax= 0.30，考慮到在沖擊荷載作用下薄壁材料容易出現(xiàn)斷裂、絕熱剪切等破壞情況，在計算模型中增加最大剪應(yīng)變失效準則，剪切失效應(yīng)變 τmax= 0.20[9]。利用單元刪除技術(shù)來模擬鋼板的破壞情況，即當單元變形引起的塑性應(yīng)變或剪切應(yīng)變超過臨界值時，就認為該單元破壞。

3 殼單元與實體單元模擬結(jié)果比較

3.1 方形鋼板的動態(tài)響應(yīng)

1）t = 0.15 ms 時刻，方形鋼板 Von Mises 應(yīng)力云圖對比。

從圖 10 可知，無論是采用殼單元還是實體單元，方形板中的應(yīng)力分布大致相同，最大值均出現(xiàn)在炸藥最上方中心位置處。同一時刻，實體單元中最大應(yīng)力值為 503.6 MPa，殼單元中最大應(yīng)力值為 504.7 MPa。

從圖 11 可知，方形鋼板的變形大致呈現(xiàn)為中心處的外鼓，呈“金字塔”狀，開始先發(fā)生中心處的翹曲，隨著時間的增加，翹曲越來越嚴重，而后慢慢下降，2 種方法模擬的鋼板變形發(fā)展趨勢與試驗結(jié)果變形趨勢相同。

2）鋼板中線上各點處最大位移值對比

圖 12給出了在 20 g 炸藥量產(chǎn)生的爆炸荷載作用下鋼板中線上各點的最大位移曲線。由圖中可以看出，2種方法得到的各點處最大位移曲線基本一致，與試驗值較為吻合，各點處最大位移差值均在 10% 以內(nèi)。最大的位移發(fā)生在鋼板中心位置處，2 種模擬方法得到的結(jié)果與試驗值差值均在 5% 以內(nèi)，其他位置處的差值較小。

3.2 兩種方法計算時間對比

2 種模型的炸藥和空氣域單元數(shù)目相同，而使用殼單元建模使得鋼板單元數(shù)目減少 50%。由圖 13 可以看出使用殼單元建模計算時間降低約 50%，顯著提高計算效率。尤其對于單元數(shù)目比較多的模型，殼單元的優(yōu)勢更為突出。

無論是從鋼板的動力響應(yīng)歷程還是鋼板的變形特征層面上的對比驗證，其結(jié)構(gòu)均表明：本次數(shù)值模擬的建模技術(shù)和參數(shù)選取合理、可靠，其數(shù)值模擬的結(jié)果可以艦船艙室內(nèi)爆選取單元提供依據(jù)。

4 結(jié) 語

對方形鋼板在爆炸荷載作用下的響應(yīng)，分別采用殼單元和實體單元進行模擬。對比 2 種方法中板的應(yīng)力分布和中線處各點的最大位移值，得到的結(jié)果大致相同，差值在 10% 以內(nèi)。結(jié)果表明使用實體單元建模和使用殼單元建模對于等厚度方板的計算結(jié)果影響不大。而采用殼單元建模相比采用實體單元建模可以顯著提高計算效率。

數(shù)值分析方法在爆炸荷載作用下艦船艙室的動力響應(yīng)研究中有著廣泛的應(yīng)用，但如何正確地選擇單元進行分析是數(shù)值模擬中很重要的一個環(huán)節(jié)。本文通過對殼單元和實體單元模擬方形鋼板計算結(jié)果的初步對比，為艦船艙室模型中如何選擇單元提供參考。

[1]陳攀, 劉志忠. 艙室內(nèi)爆沖擊波載荷特性及影響因素[J]. 艦船科學(xué)技術(shù), 2016, 38(2): 43–48. CHEN Pan, LIU Zhi-zhong. Research on loading of explosive and influencing factors inside closed cabin [J]. Ship Science and Technology, 2016, 38(2): 43–48.

[2]朱建方, 王偉力, 曾亮. 艦艇艙室內(nèi)爆毀傷的建模與仿真分析[J]. 系統(tǒng)仿真學(xué)報, 2009, 21(22): 7066–7068.

[3]杜志鵬, 李曉彬, 夏利娟, 等. 反艦導(dǎo)彈攻擊艦船舷側(cè)防護結(jié)構(gòu)過程數(shù)值仿真[J]. 哈爾濱工程大學(xué)學(xué)報, 2006. 27(4): 484–487.

[4]LONGDON G S, CHUNG Ki Y S, NURICK G N. Experimental and numercial studiesion the response of quadrangular stiffened plates. Part I: subjected to uniform blast load [J]. In-ternational Journal of Impact Engineering. 2005, (31): 55–83.

[5]LONGDON G S, CHUNG Ki Y S, NURICK G N. Experimental and numercial studiesion the response of quadrangular stiffened plates. Part II: localist blast loading [J]. International Journal of Impact Engineering. 2005, (31): 85–11.

[6]徐新光. 梁單元和實體單元模擬框架結(jié)構(gòu)的差異性比較[J].四川建材, 2014, 40(5): 46–48.

[7]ANSYS工程結(jié)構(gòu)數(shù)值分析[M]. 北京. 人民交通出版社, 2012: 11–13.

[8]ANSYS/LSDYNA在爆炸與沖擊領(lǐng)域內(nèi)的工程應(yīng)用[M]. 北京:中國建筑工業(yè)出版社, 2011.

[9]李帆. 船艦室內(nèi)爆炸破壞的數(shù)值模擬和試驗研究[D]. 北京: 國防科技大學(xué), 2012.

Comparision of differentiation with the steel plates' dynamic response subjected to blast loads using shell element and solid element modeling methods

ZHU Jun-jie1, LU Yong-gang2, FENG Xing-kui3, TAO Jun-lin1
(1. School of Civil Engineering and Architecture, Southwest University of Science and Technology, Mianyang 621010, China; 2. Institute of Structural Mechanics,China Academy of Engineering Physics, Mianyang 621900, China; 3. The Office of Rocket Force, China Academy of Eugineering Physicies, Mianyang 621900, China)

Take the steel plates of warship cabins as research objects, numerical simulation study on the response of the quadrangular steel plates subjected to blast loads is presented by dynamic analysis finite element procedure Ls-dyan3d, using modeling methods of shell element and solid element respectively. Comparing pressure distribution cloud picture and each point of midcourt line miximum displacement curve graph, it can be found that transformation trend of steel plates is similar for different thickness while each point of midcourt line miximum displacement is discriminating. However results from two methods are alike for equal thickness. Meanwhile, it can improve calculative effciency. How to choose which kind of element types correctly to analysis is an important link. Through the preliminary contrast study of calculative results of shell element and solid element, it can provide reference for choosing element type in warship cabins numerical simulation.

marine steel plate；blast loads；numerical simulation；shell element

TJ03

A

1672–7619(2017)03–0018–05

10.3404/j.issn.1672–7619.2017.03.004

2016–09–18；

2016–10–27

國家自然科學(xué)基金資助項目(11502258)

朱俊杰(1991–)，女，碩士研究生，研究方向為結(jié)構(gòu)沖擊與爆炸。

陶俊林(1972–)，男，博士，教授，E-mail:junlintao@126.com。