趙傳榮，孔德仁

(1.南京理工大學(xué)機(jī)械工程學(xué)院,江蘇南京 210094; 2.安徽工業(yè)大學(xué)電氣與信息工程學(xué)院，安徽馬鞍山 243032)

1 引 言

平面沖擊波在高能密度材料的合成[1-2]、炸藥沖擊起爆[3]及材料的動(dòng)態(tài)特性分析[4]等研究領(lǐng)域有廣泛的應(yīng)用。因此，準(zhǔn)確地獲取研究所需的平面沖擊波至關(guān)重要。

目前實(shí)驗(yàn)室常采用氣體炮[5]驅(qū)動(dòng)飛片撞擊靶板產(chǎn)生所需的平面沖擊波。沖擊波壓力峰值和峰值平臺(tái)壓力保持時(shí)間(脈沖寬度)是表征平面沖擊波強(qiáng)度的兩個(gè)重要特征指標(biāo)。湯文輝等[6]基于流體動(dòng)力學(xué)理論提出了靶板撞擊面沖擊波壓力峰值和脈沖寬度的計(jì)算模型，在一定程度上為沖擊波壓力的獲取提供了理論指導(dǎo)。但是，該模型是在一定的理想假設(shè)條件下根據(jù)飛片撞擊靶板過(guò)程中起支配作用的物理定律而獲得的，需要飛片或靶板材料在高壓下的物態(tài)方程及相關(guān)物態(tài)參數(shù)(如沖擊雨貢紐參數(shù))的準(zhǔn)確資料，才能對(duì)靶板撞擊面的沖擊波壓力做出預(yù)估。然而，當(dāng)飛片或靶板材料為沖擊雨貢紐參數(shù)未知的材料，如99.99%(質(zhì)量分?jǐn)?shù))純鎳、96%(質(zhì)量分?jǐn)?shù))氧化鋁時(shí)，將不能利用該模型預(yù)估靶板撞擊面沖擊波的壓力峰值和脈沖寬度。目前在涉及到相關(guān)研究時(shí)，只能分別利用Ni200[7]和Al300[8]的沖擊雨貢紐參數(shù)進(jìn)行粗略的估算。此外，即使飛片和靶板材料的雨貢紐參數(shù)已知，不同批次的飛片或靶板的雨貢紐參數(shù)也有所差異。因此，有必要對(duì)飛片與靶板撞擊面的壓力峰值和脈沖寬度的理論模型做進(jìn)一步研究。

本工作采用量綱分析法對(duì)影響靶板撞擊面沖擊波壓力峰值和脈沖寬度的因素進(jìn)行分析，以期獲得靶板撞擊面沖擊波壓力峰值和脈沖寬度的經(jīng)驗(yàn)?zāi)Ｐ停徊⑼ㄟ^(guò)實(shí)例分析，檢驗(yàn)該經(jīng)驗(yàn)?zāi)Ｐ偷目尚行浴?/p>

2 平面沖擊波的經(jīng)驗(yàn)?zāi)Ｐ?/h2>

將圓柱形平頭飛片撞擊半無(wú)限厚靶板[9]的過(guò)程看作平面平行正碰撞，影響靶板撞擊面沖擊波壓力峰值pm的因素有：飛片速度v，飛片和靶板的材料密度ρ01、ρ02，飛片和靶板的雨貢紐參數(shù)C01、S1和C02、S2；影響脈沖寬度τ的因素有：飛片速度v，飛片直徑D及長(zhǎng)度L，飛片和靶板的材料密度ρ01、ρ02，飛片和靶板的雨貢紐參數(shù)C01、S1和C02、S2，飛片材料的泊松比ν及Grüneisen系數(shù)γ1。由于γ1、S1、ν及S2為無(wú)量綱量，為便于研究，將它們合并為γ1/S1、S1/ν和S1/S2。pm和τ與上述影響因素之間存在如下函數(shù)關(guān)系

(1)

為簡(jiǎn)化量綱分析過(guò)程，將與兩個(gè)因變量pm和τ相關(guān)聯(lián)的各影響因素合并成一個(gè)函數(shù)關(guān)系式

(2)

在飛片撞擊靶板的物理過(guò)程中，選擇質(zhì)量、長(zhǎng)度和時(shí)間為基本量綱，分別記作M、L、T，由相似理論的Buckingham定理[10-11]，物理方程可變換成無(wú)量綱組合量(π數(shù))的關(guān)系方程式，且無(wú)量綱關(guān)系式中應(yīng)當(dāng)涉及9個(gè)(12-3=9)無(wú)量綱量。由于γ1/S1、S1/ν和S1/S2已是無(wú)量綱量，故只需用量綱分析方法再求取6個(gè)(9-3=6)無(wú)量綱參量即可。無(wú)量綱π項(xiàng)的一般指數(shù)形式為

(3)

式中：x1～x9分別表示各相應(yīng)物理量的指數(shù)。寫出各相應(yīng)物理量的量綱，列于表1。

表1 相應(yīng)物理量的量綱Table 1 Dimensions of the corresponding physical quantities

根據(jù)量綱的齊次性原則可得各指數(shù)的聯(lián)立方程組

(4)

(4)式中有9個(gè)未知數(shù)，需指定6個(gè)才能得到一組解，即一個(gè)無(wú)量綱的組合量。用x4～x9表示x1～x3,有

(5)

若指定一組x4～x9，由(5)式求出相應(yīng)的x1～x3，就可獲得所需的6個(gè)物理量，如表2所示。根據(jù)表2有：π1=L/D，π2=L2ρ01/(pmτ2)，π3=τC01/L，π4=L2ρ02/(pmτ2)，π5=τC02/L，π6=τv/L。參考(1)式，將π1～π66個(gè)量做適當(dāng)變換，得到6個(gè)物理意義較為明確的物理量，即pm/(ρ01C01v)、τv/L、L/D、ρ01/ρ02、C01/C02、v/C01。

表2 相似比求解Table 2 Solution of similar ratio

綜合考慮(1)式及另3個(gè)無(wú)量綱量γ1/S1、S1/ν和S1/S2，可得到靶板撞擊面沖擊波壓力峰值和脈沖寬度的無(wú)量綱經(jīng)驗(yàn)?zāi)Ｐ?/p>

(6)

式中：a1、b1、c1、d1、a2、b2、c2、d2、e2、f2、g2、h2均為常數(shù)。顯然(6)式中各項(xiàng)都有明確的物理意義：ρ01/ρ02、C01/C02、γ1/S1、S1/S2及S1/ν反映了飛片和靶板的材料屬性對(duì)壓力峰值pm和脈沖寬度τ的影響,v/C01反映了飛片速度的影響,L/D反映了飛片結(jié)構(gòu)尺寸對(duì)脈沖寬度τ的影響。

3 實(shí) 例

楊超等[12-13]基于?100 mm一級(jí)輕氣炮分別進(jìn)行了3、4及5 mm厚飛片撞擊靶板的加載實(shí)驗(yàn)，飛片材料是45鋼，飛片直徑為100 mm,靶板材料為純鎳(質(zhì)量分?jǐn)?shù)為99.99%)。在6種沖擊實(shí)驗(yàn)條件下獲得的飛片速度值、靶板撞擊面沖擊波壓力峰值和脈沖寬度數(shù)據(jù)列于表3，將前4組實(shí)驗(yàn)數(shù)據(jù)用于建立經(jīng)驗(yàn)?zāi)Ｐ?，后兩組則用于檢驗(yàn)經(jīng)驗(yàn)?zāi)Ｐ偷目尚行浴?/p>

表3 實(shí)驗(yàn)數(shù)據(jù)Table 3 Experimental data

由于飛片和靶板的材料以及飛片的直徑一定，(6)式中的D、ρ01、C01、ρ02、C02、γ1、S1、S2及ν都是常數(shù)，采用g-cm-μs單位制，則(6)式可簡(jiǎn)化為如下經(jīng)驗(yàn)?zāi)Ｐ?/p>

(7)

式中：s1、s2、t1、t2及λ均為常數(shù)。

采用數(shù)據(jù)回歸的方法對(duì)表3中前4組實(shí)驗(yàn)數(shù)據(jù)進(jìn)行處理，可獲得(7)式中s1、s2、t1、t2及λ的值，建立靶板撞擊面沖擊波壓力峰值和脈沖寬度經(jīng)驗(yàn)?zāi)Ｐ偷亩筷P(guān)系，如下

(8)

用后兩組實(shí)驗(yàn)數(shù)據(jù)對(duì)經(jīng)驗(yàn)?zāi)Ｐ?8)式進(jìn)行檢驗(yàn)，結(jié)果如表4所列。根據(jù)經(jīng)驗(yàn)?zāi)Ｐ?8)式計(jì)算得到的沖擊波壓力峰值的預(yù)測(cè)值與實(shí)驗(yàn)結(jié)果之間的相對(duì)誤差小于9%，壓力脈沖寬度的預(yù)測(cè)值與實(shí)驗(yàn)結(jié)果之間的相對(duì)誤差小于12%。

表4 壓力峰值和脈沖寬度的實(shí)測(cè)值與模型預(yù)測(cè)值的對(duì)比結(jié)果Table 4 Comparison of the measured values of peak pressure and pulse widthwith those predicted by this model

Note：pm,jandτjare the predicted values of peak pressure (pm) and pulse width (τ),respectively.

4 結(jié) 論

采用量綱分析法對(duì)影響壓力峰值和脈沖寬度的因素進(jìn)行分析，分別獲得了靶板撞擊面沖擊波壓力峰值和脈沖寬度的經(jīng)驗(yàn)?zāi)Ｐ?。通過(guò)實(shí)例分析，獲得了45鋼撞擊純鎳產(chǎn)生的平面沖擊波的經(jīng)驗(yàn)?zāi)Ｐ投筷P(guān)系式，并且經(jīng)驗(yàn)?zāi)Ｐ偷念A(yù)測(cè)值與實(shí)驗(yàn)測(cè)量值較為吻合。該經(jīng)驗(yàn)?zāi)Ｐ偷奶岢?，使直接開(kāi)展沖擊加載實(shí)驗(yàn)建立不同沖擊條件下平面沖擊波的預(yù)估模型成為可能。

[1] KIM M,YOO C S.Highly repulsive interaction in novel inclusion D2-N2compound at high pressure:Raman and X-ray evidence [J].J Chem Phys,2011,134(4):44-51.

[2] 羅 寧,胡 強(qiáng),李尚升.人造金剛石的合成技術(shù)分析 [J].應(yīng)用技術(shù),2010(20):257.

LUO N,HU Q,LI S S.Analysis of artificial diamond synthesis technology [J].The Application of Technology,2010(20):257.

[3] PULHAM C R,MILLAR D I A,OSWALD I D H,et al.High-pressure studies of energetic materials [M]//Boldyreva E,Dera P.High-pressure crystallography.Netherlands:Springer Netherlands,2010:447-457.

[4] 曾代朋,陳 軍,譚多望.超壓爆轟產(chǎn)物沖擊絕熱線的實(shí)驗(yàn)研究 [J].高壓物理學(xué)報(bào),2010,24(1):76-80.

ZENG D P,CHEN J,TAN D W.Over pressure experimental study on detonation products shock Hugoniot [J].Chinese Journal of High Pressure Physics,2010,24(1):76-80.

[5] PAVLENKO A V,BALABIN S I,KOZELKOV O E,et al.A one-stage light-gas gun for studying dynamic properties of structural materials in a range up to 40 GPa [J].Instrum Exp Tech+,2013,4(56):482-484.

[6] 湯文輝,張若棋,陳雪芳.LY12-M鋁中沖擊波衰減機(jī)理的實(shí)驗(yàn)研究 [J].高壓物理學(xué)報(bào),1988,2(3):218-226.

TANG W H,ZHANG R Q,CHEN X F.Experimental study on the mechanism of shock wave attenuation in LY12-M aluminum [J].Chinese Journal of High Pressure Physics,1988,2(3):218-226.

[7] CHAN K S,TIAN J W,YANG B,et al.Evolution of slip morphology and fatigue crack initiation in surface grains of Ni200 [J].Metall Mater Trans A,2009,40(11):2545-2556.

[8] 杜曉松,郝建德,楊邦朝,等.薄膜式錳銅傳感器的超高壓力標(biāo)定 [J].儀器儀表學(xué)報(bào),2001,22(增刊2):151-152.

DU X S,HAO J D,YANG B C,et al.Calibration of thin-film ultrahigh pressure manganin gauges[J].Chinese Journal of Scientific Instrument,2001,22(Suppl 2):151-152.

[9] 陳衛(wèi)東,張 忠,劉家良.破片對(duì)屏蔽炸藥沖擊起爆的數(shù)值模擬和分析 [J].兵工學(xué)報(bào),2009,30(9):1187-1190.

CHEN W D,ZHANG Z,LIU J L.Numerical simulation and analysis of fragments impact initiation of shielded explosive [J].Acta Armamentarii,2009,30(9):1187-1190.

[10] 談慶明.量綱分析 [M].合肥：中國(guó)科學(xué)技術(shù)大學(xué)出版社，2007:100-110.

TAN Q M.Dimensional analysis [M].Hefei:University of Science and Technology of China Press,2007:97-98.

[11] KISELEV S P.Numerical simulation of wave formation in an oblique impact of plates by the method of molecular dynamics [J].Journal of Applied Mechanics and Technical Physics,2012,53(6):907-917.

[12] 楊 超,經(jīng)福謙,張萬(wàn)甲.沖擊波壓力和脈沖持續(xù)時(shí)間對(duì)鐵和鎳微結(jié)構(gòu)的影響 [J].兵器材料科學(xué)與工程,1996,19(3):38-44.

YANG C,JING F Q,ZHANG W J.Effects of shock pressure and pulse duration on micro structure of iron and nickel [J].Ordnance Material Science and Engineering,1996,19(3):38-44.

[13] 楊 超,經(jīng)福謙,張萬(wàn)甲,等.沖擊加載作用下鐵和鎳的高應(yīng)變率變形 [J].兵器材料科學(xué)與工程,1996,19(1):49-54.

YANG C,JING F Q,ZHANG W J,et al.High strain rate deformation of Fe and Ni by shock loading [J].Ordnance Material Science and Engineering,1996,19(1):49-54.