仲 倩,王伯良,黃 菊,惠君明

(南京理工大學(xué)化工學(xué)院,江蘇 南京 210094)

引 言

爆炸威力研究中,炸藥爆炸產(chǎn)生的空氣沖擊波超壓常常作為評價爆炸威力的一個重要參數(shù)。通常,研究人員根據(jù)實測的爆炸沖擊波超壓值與 TNT爆炸沖擊波超壓值進(jìn)行比較,分別求出各對應(yīng)點的等效 TNT藥量,從而求出該戰(zhàn)斗部爆炸威力 TNT當(dāng)量比[1]。因此,TNT炸藥爆炸沖擊波超壓值的準(zhǔn)確性直接影響到戰(zhàn)斗部爆炸威力的評估。Brode[2]及Henrych[3]等總結(jié)的超壓經(jīng)驗公式計算結(jié)果與試驗數(shù)據(jù)偏差較大,無法直接應(yīng)用到戰(zhàn)斗部威力評估中。於津[4]及楊鑫等[5]總結(jié)的經(jīng)驗公式精度有所提高,但由于缺乏大量的試驗數(shù)據(jù),其經(jīng)驗公式的適用范圍較小。所以現(xiàn)有的TN T沖擊波超壓經(jīng)驗公式需要進(jìn)一步完善。

沖擊波超壓測試系統(tǒng)的配置是影響擬合結(jié)果的一個重要因素,根據(jù)美國軍用標(biāo)準(zhǔn)的要求建立了一套沖擊波超壓測試系,精度為 14bit,采樣頻率為65MHz,與以往的測壓系統(tǒng)相比,本測試系統(tǒng)更加符合爆炸特性的要求[6-7]。

本研究采用該測試系統(tǒng)所獲得數(shù)據(jù)與理論研究相結(jié)合的方法對 TNT的沖擊波超壓公式進(jìn)行分析研究,得到 TNT爆炸沖擊波超壓經(jīng)驗公式,并將經(jīng)驗公式與試驗數(shù)據(jù)及文獻(xiàn)試驗數(shù)據(jù)進(jìn)行對比,新擬合的經(jīng)驗公式系統(tǒng)偏差較小,具有一定的適用性,為戰(zhàn)斗部威力評估提供可靠的依據(jù)。

1 實 驗

1.1 沖擊波超壓測試系統(tǒng)簡介

空中爆炸超壓有自由場超壓和壁面超壓兩種。壁面超壓是指初始沖擊波到達(dá)壁面(主要是地面)產(chǎn)生反射后壁面所承受的壓力,該壓力比自由場壓力高。由于爆炸的破壞作用等原因,要測定自由場超壓很困難,尤其對于大藥量試驗基本上是不可能的,而且實際戰(zhàn)斗部也是經(jīng)反射后的沖擊波作用毀傷目標(biāo)。因此,本研究測試的是壁面壓力。

根據(jù)美國軍用標(biāo)準(zhǔn)和爆炸相關(guān)特性,建立了一套測試系統(tǒng),如圖1所示。由聯(lián)動控制儀同步引爆爆源,爆源產(chǎn)生的沖擊波經(jīng)過各壓力傳感器的工作面時,其晶體在壓力作用下發(fā)生形變,內(nèi)部產(chǎn)生“極化現(xiàn)象”,在垂直于電軸的表面上產(chǎn)生電荷并由其內(nèi)置的電荷放大器對信號進(jìn)行放大,經(jīng)過信號調(diào)理儀傳送到多通道的VXI數(shù)據(jù)采集儀,然后由計算機(jī)系統(tǒng)對所獲信號進(jìn)行數(shù)據(jù)處理。

傳感器選用美國 PCB公司 102A系列,信號調(diào)理儀選用美國 PCB公司 481A型號;數(shù)據(jù)采集儀是成都微測公司VXI-1115型,14位的精度,65M的采樣率,符合美國軍用標(biāo)準(zhǔn)。

圖1 測試系統(tǒng)示意圖Fig.1 Measurement system chart

1.2 靶場布置和試驗裝置

將 TNT炸藥放在距地面1m的支架上,以藥柱在地面上的垂直投影點為爆心,在爆心相對垂直的兩條射線上布置地面?zhèn)鞲衅?傳感器分別布置在距爆心 3、 4、 5、 6、 8、 12、 14、 17、 20、 22和 24m 處。

1.3 樣 品

試驗樣品為熔鑄 TNT藥柱,藥量分別為 1、8、10、30和 35kg,裝藥的長徑比為 1∶1,以鈍化 RDX為傳爆藥,用 8號雷管上端起爆。

2 結(jié)果與分析

2.1 沖擊波超壓的經(jīng)驗公式

Brode[2]用模型相似理論給出 TNT爆炸沖擊波峰值超壓的表達(dá)式為:

式中:Δp為沖擊波超壓,MPa;R為比例距離(R=r/3r為測試點到爆心的距離,m;W為藥量,kg)。

Henrych[3]用試驗的方法得到?jīng)_擊波峰值超壓公式為:

Mills用相似理論及數(shù)值模擬相結(jié)合的方法得到 TNT爆炸沖擊波超壓的表達(dá)式:

M.A.Sadovskyi[9]根據(jù)模型相似律理論建立比例距離關(guān)系比較式,由試驗確定系數(shù),得到高爆炸藥沖擊波峰值超壓的表達(dá)式為:

對于大藥量 TNT炸藥爆炸,文獻(xiàn)[10]中給出TNT在近地面爆炸時沖擊波超壓的計算公式:

文獻(xiàn)[5]給出 TNT爆炸的沖擊波超壓計算公式:

式(6)只是對 M.A.Sadovskyi經(jīng)驗公式在比例距離小于 1m/kg1/3情況下的經(jīng)驗公式進(jìn)行了修正,不具有代表性。

2.2 經(jīng)驗公式計算值與測試結(jié)果的對比分析

利用沖擊波超壓測壓系統(tǒng)測試了兩組 TNT的超壓,每組試驗用不同藥量分別做3次平行試驗。為直觀地分析經(jīng)驗公式計算值與試驗數(shù)據(jù)之間的差別,繪出不同方法計算沖擊波超壓比例距離的關(guān)系曲線,如圖 2所示。

由圖2可知,Brode和 Henrych方法的計算值偏低,Mills經(jīng)驗公式的計算值和文獻(xiàn) [5]的計算值較其他公式相對偏高,尤其是在-R 小 于 2.5m/kg1/3范圍內(nèi),其計算結(jié)果與試驗結(jié)果偏差較大。隨著比例距離-R的 增大,經(jīng)驗公式的計算結(jié)果與試驗數(shù)據(jù)的偏差逐漸減小,當(dāng)比例距離大于3m/kg 時,兩者基本比較接近。

分析認(rèn)為,產(chǎn)生上述偏差的主要原因:一是超壓峰值具有一定的離散性;二是不少研究者的實驗研究進(jìn)行得比較早,測試儀器相對簡單,而超壓峰值的獲取對實驗裝置的靈敏度、精度及采樣頻率的要求都很高。

圖 2 測試值與經(jīng)驗公式峰值超壓對比Fig.2 Comparison of the measured veluesΔp with eiteratureonesΔp

因此,在本測試系統(tǒng)所測得數(shù)據(jù)基礎(chǔ)上對 TNT的經(jīng)驗公式重新擬合,可為戰(zhàn)斗部威力評估提供客觀依據(jù)。

2.3 公式擬合

2.3.1 計算原理

根據(jù)爆炸相似規(guī)律,超壓與比例距離的關(guān)系式為:

根據(jù)試驗測試數(shù)據(jù),以公式(7)為基礎(chǔ),運用最小二乘法原理,對 TNT進(jìn)行經(jīng)驗公式擬合。

2.3.2 擬合公式

借助爆炸相似率和最小二乘法原理,對TNT超壓測試數(shù)據(jù)進(jìn)行經(jīng)驗公式擬合,擬合公式如下:

2.4 分析與討論

2.4.1 相對偏差分析

經(jīng)驗公式的擬合過程是一個理想化的爆炸相似律的模擬,從每次測試數(shù)據(jù)來看,數(shù)據(jù)曲線并不像擬合的結(jié)果那樣光滑。為驗證上述擬合公式相對理論計算偏差,表 1列出了用擬合公式(8)計算的 TNT超壓預(yù)測值與本次試驗測試值。

從表 1可以看出,試驗測試數(shù)據(jù)與公式擬合值之間的相對偏差均小于 5.30%,其絕對值的平均值為2.65%,正向偏差的平均值為6.17%,負(fù)向偏差平均值為5.62%。該公式可用于TNT沖擊波超壓的計算,能夠滿足后續(xù)的威力評價與理論分析的要求。

表 1 TNT超壓值的試驗數(shù)據(jù)和計算值偏差分析Table 1 Error analysis of the experimental data and caluclated ones of the TNT overpressure

2.4.2 與參考文獻(xiàn)數(shù)據(jù)的比較

為驗證上述擬合公式適用的普遍性,選取文獻(xiàn)[11]中的 TNT實測數(shù)據(jù)帶入擬合公式進(jìn)行比較分析,文獻(xiàn)[11]中的測試方法及系統(tǒng)配置與本研究相似,分析結(jié)果如表 2所示。

由表 2可見,文獻(xiàn)中測試值與擬合值之間的相對偏差絕對值的平均值為5.61%,說明該擬合公式可普遍適用于 TNT沖擊波超壓的計算。

表2 與文獻(xiàn)試驗數(shù)據(jù)對比分析Table2 ComparisonoftheexperimentalvaluesΔp withtheexperimentalonesΔpinliterature

3 結(jié) 論

(1)根據(jù)美國軍用標(biāo)準(zhǔn)和爆炸相關(guān)特性,建立了一套高精度和高分辨率的沖擊波超壓測試系統(tǒng)。

(2)總結(jié)出能較好描述沖擊波超壓峰值與對比距離相似律關(guān)系的表達(dá)式。通過與試驗數(shù)據(jù)分析比較,可知該公式的預(yù)測值與試驗值相對偏差均小于5.30%,其絕對值的平均值為2.65%。

(3)通過對文獻(xiàn)試驗數(shù)據(jù)與本研究表達(dá)式預(yù)測值的結(jié)果比較,其相對偏差絕對值的平均值為5.61%,表明該擬合公式可適用于TNT沖擊波超壓的計算。

[1] 北京工業(yè)學(xué)院八系.爆炸及其作用 (下冊)[M].北京:國防工業(yè)出版社,1979.

[2] Brode H L.Blast wave from a spherical charge [J].Phys Fluids,1959(2):217.

[3] Henrych J.The Dynamics of Explosion and Its Use [M].Amsterdam:Elsevier,1979.

[4] 於津,彭金華,張?zhí)?等.基于 MATLAB的 FAE戰(zhàn)斗部超壓計算的公式擬合 [J].彈箭與制導(dǎo)學(xué)報,2004,24(3):306-308.YU Jin, PENG Jin-hua,ZHANG Tao,etal.Formula fitting to overpressure calculation of FAE w arhead in MATLAB[J]. Journal of Projectiles,Rockets,Missiles and Guidance,2004,24(3):306-308.

[5] 楊鑫,石少卿.空氣中 TNT爆炸沖擊波超壓峰值的預(yù)測及數(shù)值模擬[J].爆破,2008(2):15-19.YANG Xin,SHI Shao-qing.Forecast and simulation of peak overpressure of TNT explosion shockwave in air [J].Blasting,2008(2):15-19.

[6] 王建玲,郭煒 ,馮曉軍.TNT、PBX和 Hexel空中爆炸沖擊波參數(shù)的實驗研究[J].火炸藥學(xué)報,2008,31(6):42-44.WANG Jian-ling, GUO Wei, FENG Xiao-jun.Experimental research on the air explosion shock wave parameters of TNT,PBX and hexel[J].Chinese Journal of Explosives and Propellants,2008,31(6):42-44.

[7] Muzychuk E,May seless M,Belsky I.Influence of thepost-detonation burning process on blast wave parameters in air [C]∥22nd International Symposium on Ballistics.Canada:Destech Publications,2005:536-543.

[8] Ah med Fah my Farag Tolba.Response of FRP-Retrofitted Reinforced Concrete Panels to Blast Loading [D].Canada:Carleton University,2001.

[9] Sadovskyi M A.Mechanical action of air shock waves of explosion,based on experimental data [M].Moscow:Izd Akad Nauk SSSR,1952.

[10]惠君明,陳天云.炸藥爆炸理論 [M].南京:江蘇科學(xué)技術(shù)出版社,1995.

[11]張?zhí)?惠君明.FAE爆炸場超壓與威力的實驗研究[J].爆炸與沖擊,2004(3):176-181.ZHANG Tao,HHI Jun-ming.Experimental research on the overpressure and power in the FAEblast field [J].Explosive and Blast,2004(3):176-181.

致謝:感謝范能躍、房瑩瑩、王鳳丹等對本實驗的幫助。