趙傳榮，孔德仁

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

1 引 言

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

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

本工作采用量綱分析法對影響靶板撞擊面沖擊波壓力峰值和脈沖寬度的因素進行分析，以期獲得靶板撞擊面沖擊波壓力峰值和脈沖寬度的經(jīng)驗?zāi)Ｐ?；并通過實例分析，檢驗該經(jīng)驗?zāi)Ｐ偷目尚行浴?/p>

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

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

(1)

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

(2)

在飛片撞擊靶板的物理過程中，選擇質(zhì)量、長度和時間為基本量綱，分別記作M、L、T，由相似理論的Buckingham定理[10-11]，物理方程可變換成無量綱組合量(π數(shù))的關(guān)系方程式，且無量綱關(guān)系式中應(yīng)當(dāng)涉及9個(12-3=9)無量綱量。由于γ1/S1、S1/ν和S1/S2已是無量綱量，故只需用量綱分析方法再求取6個(9-3=6)無量綱參量即可。無量綱π項的一般指數(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個未知數(shù)，需指定6個才能得到一組解，即一個無量綱的組合量。用x4～x9表示x1～x3,有

(5)

若指定一組x4～x9，由(5)式求出相應(yīng)的x1～x3，就可獲得所需的6個物理量，如表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個量做適當(dāng)變換，得到6個物理意義較為明確的物理量，即pm/(ρ01C01v)、τv/L、L/D、ρ01/ρ02、C01/C02、v/C01。

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

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

(6)

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

3 實 例

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

表3 實驗數(shù)據(jù)Table 3 Experimental data

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

(7)

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

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

(8)

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

表4 壓力峰值和脈沖寬度的實測值與模型預(yù)測值的對比結(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é) 論

采用量綱分析法對影響壓力峰值和脈沖寬度的因素進行分析，分別獲得了靶板撞擊面沖擊波壓力峰值和脈沖寬度的經(jīng)驗?zāi)Ｐ?。通過實例分析，獲得了45鋼撞擊純鎳產(chǎn)生的平面沖擊波的經(jīng)驗?zāi)Ｐ投筷P(guān)系式，并且經(jīng)驗?zāi)Ｐ偷念A(yù)測值與實驗測量值較為吻合。該經(jīng)驗?zāi)Ｐ偷奶岢?，使直接開展沖擊加載實驗建立不同沖擊條件下平面沖擊波的預(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] 羅 寧,胡 強,李尚升.人造金剛石的合成技術(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)物沖擊絕熱線的實驗研究 [J].高壓物理學(xué)報,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].高壓物理學(xué)報,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é)報,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)東,張 忠,劉家良.破片對屏蔽炸藥沖擊起爆的數(shù)值模擬和分析 [J].兵工學(xué)報,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].合肥：中國科學(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)福謙,張萬甲.沖擊波壓力和脈沖持續(xù)時間對鐵和鎳微結(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)福謙,張萬甲,等.沖擊加載作用下鐵和鎳的高應(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.