基于Held判據(jù)的線(xiàn)型聚能切割器結(jié)構(gòu)優(yōu)化

曹 濤，陳 煒，金建峰，孫 浩，路耀斌

(中國(guó)華陰兵器試驗(yàn)中心， 陜西 華陰 714200)

聚能射流因具有作用方向單一且能量密度高的特點(diǎn)，被越來(lái)越多地應(yīng)用于未爆彈藥銷(xiāo)毀[1-2]。相比較于傳統(tǒng)殉爆的銷(xiāo)毀方法，聚能射流所產(chǎn)生的次生危害更小，對(duì)銷(xiāo)毀場(chǎng)地要求更低，且無(wú)需與未爆彈藥直接接觸，極大地提高了銷(xiāo)毀作業(yè)的安全性。目前英國(guó)、瑞典等歐洲國(guó)家在聚能射流引爆彈藥的研究方面已經(jīng)較為成熟，形成了一批聚能引爆彈產(chǎn)品[3-6]。國(guó)內(nèi)研究人員在聚能銷(xiāo)毀方面也取得了不少成果。宋桂飛等[7]設(shè)計(jì)了一種軸對(duì)稱(chēng)結(jié)構(gòu)的聚能引爆器，并對(duì)某壁厚超過(guò)10 mm的彈藥進(jìn)行了銷(xiāo)毀，成功引爆彈藥內(nèi)部含能物質(zhì)。徐全軍等[8]根據(jù)Held引爆理論，采用正交分析法，設(shè)計(jì)了一種用于銷(xiāo)毀簡(jiǎn)易爆炸裝置的軸對(duì)稱(chēng)聚能彈，并通過(guò)試驗(yàn)驗(yàn)證了聚能彈的微聲作用效果。王敏等[9]對(duì)乳化炸藥線(xiàn)型切割器銷(xiāo)毀薄殼彈藥的過(guò)程進(jìn)行了數(shù)值模擬，并對(duì)藥型罩的材料和炸高作了優(yōu)化。在開(kāi)展了大量彈藥銷(xiāo)毀工作之后，發(fā)現(xiàn)有一類(lèi)彈藥口徑比較小、但數(shù)量往往較大。這類(lèi)彈藥采用殉爆法銷(xiāo)毀時(shí)極易飛散和銷(xiāo)毀不完全，而利用射孔彈進(jìn)行逐發(fā)銷(xiāo)毀時(shí)工作量大，且成本較高。因此為完善聚能銷(xiāo)毀器材系列，本研究提出利用線(xiàn)型聚能切割器對(duì)小口徑彈藥進(jìn)行批量銷(xiāo)毀。

在線(xiàn)型聚能切割器破甲能力試驗(yàn)的基礎(chǔ)上，借助有限元軟件ANSYS/LS-DYNA對(duì)線(xiàn)型聚能切割器射流成型和侵徹靶板的過(guò)程進(jìn)行了數(shù)值模擬，并對(duì)比試驗(yàn)結(jié)果校正了相關(guān)參數(shù)。從試驗(yàn)的結(jié)果來(lái)看，該線(xiàn)型聚能切割器雖然能達(dá)到引爆未爆彈要求，但是在產(chǎn)品長(zhǎng)度、總質(zhì)量等方面還有待進(jìn)一步改進(jìn)。本研究利用正交設(shè)計(jì)原理，基于Held引爆理論對(duì)線(xiàn)型聚能切割器藥型罩開(kāi)口角度、壁厚以及罩頂藥厚進(jìn)行了優(yōu)化，為下一步的改進(jìn)設(shè)計(jì)提供參考。

1 試驗(yàn)與數(shù)值模擬

1.1 線(xiàn)型聚能切割器破甲能力試驗(yàn)

為了確定線(xiàn)型聚能切割器銷(xiāo)毀小口徑彈藥的可行性，對(duì)某型線(xiàn)型聚能切割器進(jìn)行了破甲能力試驗(yàn)。試驗(yàn)設(shè)置如圖1所示，線(xiàn)型聚能切割器平放在兩塊疊加的45#鋼靶板上面，上下靶板厚度均為35 mm。切割器殼體底部直接接觸上靶板表面，內(nèi)置炸高20 mm，采用端面單點(diǎn)起爆。

破甲能力試驗(yàn)共進(jìn)行了3組，除第一組外，該型線(xiàn)型聚能切割器均能成功切開(kāi)35 mm厚的上層靶板，切口中間寬兩頭窄，平均寬度26 mm，上層靶破甲后效如圖2所示。第一組試驗(yàn)上層靶板的侵徹深度分為兩部分，前100 mm未切開(kāi)，后100 mm完全切開(kāi)，分析原因可能是起爆裝置固定不牢靠，與線(xiàn)型切割器之間存在間隙，導(dǎo)致初始起爆能量不充足。試驗(yàn)表明在可靠起爆的情況下該型線(xiàn)型聚能切割器對(duì)45#鋼的破甲能力為35 mm。

圖1 破甲能力試驗(yàn)設(shè)置圖

圖2 破甲能力試驗(yàn)上層靶板后效

1.2 數(shù)值模擬

在A(yíng)NSYS中建立了該型線(xiàn)型聚能切割器切割靶板的等比例模型，模型由炸藥、藥型罩、殼體、靶板和空氣5部分組成，均采用單層實(shí)體網(wǎng)格，并在長(zhǎng)度方向的斷面上施加約束，這樣既可以準(zhǔn)確描述線(xiàn)型射流的特點(diǎn)，又可以大大減少計(jì)算時(shí)間。由于線(xiàn)型聚能切割器是面對(duì)稱(chēng)結(jié)構(gòu)，所以只需建立對(duì)稱(chēng)面一側(cè)的1/2模型。炸藥、藥型罩和空氣采用歐拉網(wǎng)格，殼體和靶板采用拉格朗日網(wǎng)格。聚能射流成型過(guò)程涉及材料的大變形問(wèn)題，為了避免因單元變形量過(guò)大而導(dǎo)致計(jì)算終止，采用多物質(zhì)ALE算法。流體與固體之間采用流固耦合算法，這樣無(wú)需設(shè)置接觸即可實(shí)現(xiàn)各項(xiàng)力學(xué)參數(shù)的相互傳遞。各部分的材料模型和狀態(tài)方程如表1所示，具體材料參數(shù)參見(jiàn)文獻(xiàn)[10]。

表1 數(shù)值模擬的材料模型與狀態(tài)方程

1.3 數(shù)值模擬結(jié)果分析

數(shù)值模擬結(jié)果表明，該線(xiàn)型聚能切割器對(duì)40 mm厚鋼靶板破甲深度約為33 mm，與試驗(yàn)結(jié)果35 mm相差5.7%，驗(yàn)證了數(shù)值模擬的準(zhǔn)確性，為接下來(lái)采用數(shù)值模擬方法進(jìn)行切割器結(jié)構(gòu)優(yōu)化奠定基礎(chǔ)。

圖3 聚能射流成型與侵徹過(guò)程

2 線(xiàn)型聚能切割器結(jié)構(gòu)優(yōu)化

2.1 正交優(yōu)化因子指標(biāo)的選擇

正交試驗(yàn)是工程試驗(yàn)中常用優(yōu)化設(shè)計(jì)方法，對(duì)于因子數(shù)和水平數(shù)不太復(fù)雜的試驗(yàn)，通過(guò)少量的試驗(yàn)即可分析不同因子不同水平的優(yōu)化組合，直觀(guān)分析和方差分析得到各個(gè)因子對(duì)試驗(yàn)指標(biāo)的影響情況，確定因子對(duì)指標(biāo)影響的主次順序[11-12]。

1) 裝藥寬度

線(xiàn)型聚能切割器截面結(jié)構(gòu)如圖4所示，通常情況下，線(xiàn)型聚能切割器裝藥寬度w決定其破甲能力，通過(guò)前期的試驗(yàn)分析，認(rèn)為裝藥寬度為40 mm時(shí)，破甲能力過(guò)剩，且切割器質(zhì)量過(guò)大，影響其適用性，因此將線(xiàn)型聚能切割器寬度設(shè)計(jì)為30 mm。

圖4 線(xiàn)型聚能切割器結(jié)構(gòu)示意圖

2) 藥型罩錐角

線(xiàn)型聚能切割器采用易于加工生產(chǎn)的楔形罩，藥型罩錐角的大小決定侵徹體的形態(tài)和速度。錐角過(guò)小，射流速度顯著提高，但藥型罩材料利用率明顯下降，大部分材料將形成粗大的杵體。錐角過(guò)大，在爆轟波加載過(guò)程中藥型罩頂部發(fā)生翻轉(zhuǎn)，形成線(xiàn)型爆炸成型彈丸。實(shí)踐和理論研究[13]表明，線(xiàn)型聚能射流在作用方向的能量越大，其破甲能力越強(qiáng)，作用方向能量取最大值時(shí)，藥型罩半錐角α與罩材料彎折角θ之間存在關(guān)系：

(1)

通常線(xiàn)型聚能切割器的藥罩質(zhì)量比(裝藥質(zhì)量與藥型罩質(zhì)量之比)R在1.5～4，相應(yīng)的折轉(zhuǎn)角θ=15°～25°，由式(1)可得藥型罩錐角2α=80°～101°，正交優(yōu)化時(shí)取80°、85°、90°、95° 四個(gè)水平。

3) 藥型罩壁厚

The simulation result of a regular patch antenna resonant at 12 GHz was design as below:

常用的楔形罩有等壁厚和變壁厚兩種形式，本研究采用等壁厚楔形罩。最佳藥型罩壁厚隨著藥型罩材料比重的減小而增加，隨著藥型罩口徑的增加而增加，隨著罩頂角的增大而增加，隨著殼體厚度的增加而增加。藥型罩厚度一般約為其裝藥寬度的1～6%，正交優(yōu)化時(shí)取1.4 mm，1.6 mm，1.8 mm，2.0 mm四個(gè)水平。

此外，裝藥作為驅(qū)動(dòng)金屬藥型罩運(yùn)動(dòng)的能源物質(zhì)，在一定范圍內(nèi)增加罩頂藥厚能顯著提高線(xiàn)型切割器的性能，但是當(dāng)罩頂藥厚大于某一臨界值時(shí)切割器性能便不再有明顯改善。藥量過(guò)大時(shí)不僅導(dǎo)致材料浪費(fèi)，而且增大了總體質(zhì)量。因此藥型罩頂部裝藥厚度 也應(yīng)作為影響因子進(jìn)行分析，所設(shè)計(jì)的因子水平如表2。

表2 正交試驗(yàn)因子水平表

4) 優(yōu)化指標(biāo)的選擇

在聚能射流引爆非均質(zhì)炸藥方面，M.Held 發(fā)現(xiàn)了以速度和射流直徑為主要參量的引爆判據(jù)[14]，即：

u2d=Icr

(2)

其中：u是聚能射流撞擊速度；d是聚能射流的直徑；Icr為炸藥的敏感常數(shù)，即臨界起爆常數(shù)。若撞擊截面不是圓時(shí)，上述判據(jù)修正為：

(3)

其中，A為侵徹體對(duì)應(yīng)的截面積。

線(xiàn)型聚能切割器端面單點(diǎn)起爆以后，爆轟波的傳播是一個(gè)復(fù)雜的三維問(wèn)題，不同截面上所形成的聚能射流到達(dá)靶板的時(shí)間不同，因此很難確定射流撞擊靶板的截面面積A。但是Held判據(jù)為確定正交試驗(yàn)指標(biāo)提供了很好的思路，即依舊以u(píng)2d作為考核指標(biāo)，其中u為射流速度，d則為線(xiàn)型聚能射流寬度，u2d值越大，則說(shuō)明線(xiàn)型切割器起爆炸藥的能力越強(qiáng)。

2.2 正交試驗(yàn)方案與結(jié)果

根據(jù)正交試驗(yàn)原則[15]，對(duì)于三因子四水平的試驗(yàn)，選用L16(45)正交表進(jìn)行分析，正交表前三列分別考察罩頂藥厚、藥型罩錐角和藥型罩壁厚，后兩列作為誤差分析。一共建立16組模型，借助LS-DYNA分別計(jì)算射流成型過(guò)程，考察射流頭部速度和射流寬度，計(jì)算結(jié)果如表3所示。

表3 正交試驗(yàn)安排與計(jì)算結(jié)果

2.3 結(jié)果分析

1) 直觀(guān)分析

對(duì)計(jì)算結(jié)果進(jìn)行極差分析，將各列水平數(shù)相同的計(jì)算結(jié)果相加除以4，得到各水平的平均計(jì)算結(jié)果，記為k1、k2、k3、k4，用每個(gè)水平的最大值減去最小值即得到各因子的極差R。根據(jù)每個(gè)極差 的大小可以確定每個(gè)因子對(duì)指標(biāo)影響的重要程度。極差分析結(jié)果如表4，各因子效應(yīng)曲線(xiàn)如圖5。由極差分析的結(jié)果可以看出，線(xiàn)型聚能切割器射流引爆能力隨罩頂藥厚的增加而增加，隨藥型罩錐角的增大而減小，隨壁厚的增大而減小。各因子對(duì)引爆能力影響的重要程度為2α>δ>λ，其中藥型罩角度是影響射流起爆能力的首要因子，其次是藥型罩壁厚，影響程度最小的是罩頂藥厚，極差分析得出的最優(yōu)方案是2α=80°、λ=18 mm、δ=1.4 mm。

(4)

表4 極差分析

圖5 因子效應(yīng)曲線(xiàn)

單個(gè)因子及誤差的偏差平方和：

(5)

方差分析結(jié)果如表5所示，從方差計(jì)算結(jié)果可以看出，罩頂藥厚對(duì)線(xiàn)型聚能切割器引爆能力沒(méi)有影響，實(shí)際上直觀(guān)分析發(fā)現(xiàn)引爆能力隨著罩頂藥厚的增加而增加，可能的原因在于罩頂藥厚水平取值過(guò)大，已經(jīng)接近裝藥飽和值，因此罩頂藥厚λ取13.5 mm。錐角對(duì)線(xiàn)型切割器起爆能力的影響顯著，但不是高度顯著，分析原因在于錐角越小，射流速度越大，但同時(shí)射流質(zhì)量降低，射流寬度變小。藥型罩壁厚對(duì)線(xiàn)型切割器起爆能力的影響顯著，隨著藥型罩壁厚增加，射流速度降低，但同時(shí)射流寬度也增加，所以壁厚的變化對(duì)起爆能力影響有限。方差分析進(jìn)一步排除了誤差的影響，認(rèn)為罩頂藥厚取值接近裝藥飽和值，罩頂藥厚λ取13.5 mm。

表5 方差分析

注：F0.01(3，6)=9.78，F(xiàn)0.05(3，6)=4.76，F(xiàn)0.1(3，6)=3.29。

3 結(jié)論

1) 進(jìn)行了某型線(xiàn)型聚能切割器破甲能力試驗(yàn)，對(duì)45#鋼靶板的破甲深度達(dá)到35 mm。借助ANSYS/LS-DYNA有限元軟件，模擬了該型線(xiàn)型聚能切割器破甲過(guò)程，數(shù)值模擬結(jié)果與試驗(yàn)吻合較好。

2) 采用正交試驗(yàn)的方法對(duì)線(xiàn)型聚能切割器截面參數(shù)進(jìn)行優(yōu)化設(shè)計(jì)，根據(jù)Held引爆理論確定了優(yōu)化指標(biāo)為u2d，數(shù)值模擬得到了16組方案的射流速度和寬度，進(jìn)而確定了u2d值。通過(guò)極差分析和方差分析發(fā)現(xiàn)藥型罩角度是影響射流起爆能力的首要因子，其次是藥型罩壁厚，影響程度最小的是罩頂藥厚，得出最優(yōu)方案：2α=80°、λ=13.5 mm、δ=1.4 mm，為下一步線(xiàn)型聚能切割器改進(jìn)設(shè)計(jì)提供參考。