張友君，谷 巖，裴曉陽，杜金梅，楊 佳，喻 寅

（中國工程物理研究院流體物理研究所沖擊波物理與爆轟物理重點實驗室，四川 綿陽 621900）

高能炸藥具有較高的化學釋放能和較快的反應(yīng)速度，是一種常用的高壓沖擊波能源。在化爆加載中，利用高能炸藥爆轟釋放的化學能轉(zhuǎn)變?yōu)轱w片的動能，已有了一系列的化爆加載裝置，如接觸爆炸、空腔增壓、多級飛片增壓、二維會聚增壓等方法［1］。目前的化爆加載裝置大都采用平面波透鏡傳爆主炸藥的方式驅(qū)動飛片，這種裝置獲得的飛片擊靶速度為數(shù)千米每秒，相應(yīng)的沖擊加載壓力范圍為幾十至幾百吉帕。而在研究材料的損傷與斷裂、低壓沖擊相變以及彈塑性形變特性等動態(tài)響應(yīng)時，需要使待測樣品受到的沖擊應(yīng)力位于較低的區(qū)間。目前，獲得數(shù)百米每秒擊靶速度或產(chǎn)生幾吉帕壓力大多使用一級氣體炮，而爆轟驅(qū)動技術(shù)結(jié)構(gòu)相對簡單、體積小、容易與測試儀器同步工作（如原位實時測量、高速攝影）等［2］，因此研究化爆驅(qū)動飛片在低速和穩(wěn)定狀態(tài)下的加載技術(shù)十分必要。

在平面波透鏡結(jié)構(gòu)的加載裝置中，為了在材料中獲得較低的沖擊壓力，可采用低爆速炸藥或加惰性材料衰減層等方法［3］。В.И.Цыпкцнн［4］對TNT炸藥爆轟產(chǎn)物經(jīng)過不同厚度空氣隙驅(qū)動10mm銅飛片的運動過程進行了實驗研究，表明隨著空氣隙的增大，飛片速度下降，但空氣隙增加到一定厚度時飛片速度難以再繼續(xù)下降。趙劍衡等［5］結(jié)合實驗和數(shù)值模擬研究了在炸藥透鏡和飛片之間引入空腔對爆轟驅(qū)動飛片的影響，結(jié)果表明引入空腔的爆轟加載裝置可以推動6mm銅飛片達到約1km／s的速度。張萬甲［6］、G.I.Kanel等［7］設(shè)計了組合飛片作為低壓爆轟加載裝置，實驗結(jié)果初步表明組合飛片作為低壓爆轟加載的原理是可行的，但設(shè)計較復(fù)雜。金柯等［8］對爆轟產(chǎn)物與飛片之間預(yù)留空氣隙和使用組合飛片的方式分別進行了數(shù)值模擬研究，結(jié)果表明增加空氣隙或使用組合飛片可以使中等阻抗飛片降至1～2km／s。綜上所述，采用平面波透鏡結(jié)構(gòu)的爆轟加載裝置，很難獲得低于1km／s飛片速度的加載能力。

本文中，通過網(wǎng)絡(luò)爆轟技術(shù)設(shè)計一種新型低速低沖擊應(yīng)力的加載裝置。結(jié)合數(shù)值模擬、實驗測量與標定，對網(wǎng)絡(luò)爆轟驅(qū)動飛片的速度、平面性及裝置結(jié)構(gòu)參數(shù)進行系統(tǒng)的研究，擬為低壓動態(tài)加載實驗研究提供技術(shù)手段。

1 實驗原理與設(shè)計

1.1 基本原理

炸藥爆轟驅(qū)動飛片的過程中，可用不可壓縮飛片模型估算飛片受到壓縮波驅(qū)動的運動規(guī)律，得到飛片速度u的關(guān)系式［3］

式中：m為裝藥質(zhì)量，M為飛片質(zhì)量，l為裝藥長度，D為爆轟波速度，t為飛片飛行時間。

當t趨近于無窮時，可以得到飛片的極限速度umax

由式（2）可知，η越小，飛片速度越低。在炸藥成分與飛片質(zhì)量一定的情況下，要降低飛片速度，必須降低裝藥質(zhì)量或裝藥長度。在傳統(tǒng)的平面波透鏡化爆加載裝置中，需要一定的裝藥長度才能保證點源發(fā)出的散心爆轟波改造成為平面爆轟波，這個總長度一般都在幾十毫米。同時，在平面波透鏡的加載裝置中，通常還要再加一段主炸藥，厚度在數(shù)十毫米范圍，主裝藥長度過短就會影響飛片擊靶的平面性。因此使用平面波透鏡結(jié)構(gòu)驅(qū)動飛片的裝置很難降低飛片速度。

1.2 裝置設(shè)計

為了得到較低的飛片擊靶速度，本文中提出并使用平面網(wǎng)絡(luò)爆轟技術(shù)。網(wǎng)絡(luò)爆轟是在大口徑裝藥的端面上布置多個起爆點，這些起爆點由炸藥網(wǎng)絡(luò)線連接，當雷管引爆時，通過炸藥網(wǎng)絡(luò)線的傳爆使多個起爆點同時起爆，由此獲得平面范圍內(nèi)的爆轟波輸出。利用網(wǎng)絡(luò)起爆可以降低裝藥厚度，通常主炸藥厚度在1mm以下。

在一維應(yīng)變動態(tài)加載實驗中，對加載裝置有兩個基本要求。首先，壓力可調(diào)，有一定的壓力調(diào)節(jié)范圍；其次，所產(chǎn)生的沖擊波有一定的平面尺寸和平面度，以便物理力學參數(shù)的測量［3］。本文的網(wǎng)絡(luò)爆轟驅(qū)動飛片技術(shù)滿足這兩個基本要求。

根據(jù)上述基本原理，網(wǎng)絡(luò)爆轟驅(qū)動飛片加載的裝置如圖1所示。它由雷管、炸藥網(wǎng)絡(luò)板、衰減層、有機玻璃板、飛片、空腔和靶系統(tǒng)等組成。此處的炸藥網(wǎng)絡(luò)板包括了多點起爆網(wǎng)絡(luò)和主裝藥。

炸藥網(wǎng)絡(luò)板中裝藥厚度0.8mm，邊長95mm。衰減層采用沖擊阻抗較高的45鋼板，厚度可視所需飛片速度調(diào)節(jié)。為使多點爆轟波被改造成為準一維平面波，經(jīng)過計算取有機玻璃板厚10mm。飛片為?100mm的圓板，厚度3mm，材料為無氧銅。

圖1 網(wǎng)絡(luò)爆轟驅(qū)動飛片的實驗裝置示意圖Fig.1 Schematic diagram of network detonation－loading device

1.3 數(shù)值模擬

使用數(shù)值計算檢驗了網(wǎng)絡(luò)爆轟驅(qū)動飛片技術(shù)的可行性。裝置材料采用彈塑性流體模型，炸藥爆轟產(chǎn)物的描述選用JWL狀態(tài)方程，表達式為［9］

式中：A、B、ω、R1、R2為常數(shù)，e為比內(nèi)能。炸藥 GI－920的CJ參數(shù)為：ρ＝1.540g／cm3，D＝7.35km／s，e＝4.29MJ／kg，狀態(tài)方程參數(shù)為：A＝2 714.0GPa，B＝17.930GPa，R1＝7.00，R2＝1.60，ω＝0.35。

圖2（a）給出了飛片后自由面粒子速度歷史隨45鋼板厚度變化的計算結(jié)果，炸藥在時刻零起爆。圖2（b）給出了飛片穩(wěn)定飛行后的擊靶速度隨鋼板厚度的變化趨勢。結(jié)果表明，當鋼板厚度從1mm增加到10mm時，飛片的擊靶速度相應(yīng)從330m／s減小到210m／s，并且飛片飛行距離在2mm內(nèi)即可達到穩(wěn)定的速度平臺。應(yīng)力波在飛片中來回反射可造成飛片自由面速度振蕩；鋼板厚度與飛片速度之間并不完全成線性關(guān)系，隨著鋼板厚度持續(xù)增加，擊靶速度繼續(xù)減小的幅度十分有限。如果使用該飛片撞擊相同材料的無氧銅樣品可以達到3.8～6.2GPa沖擊應(yīng)力，因此初步計算認為該結(jié)構(gòu)的裝置具有可行性，并且通過改變衰減層鋼板的厚度可以使飛片的速度發(fā)生連續(xù)性的變化，得到不同的沖擊應(yīng)力。

圖2 飛片速度隨衰減層45鋼板厚度的變化規(guī)律Fig.2 Simulation results of flyer velocity as a function of steel thickness

2 實驗測量與標定

為了對飛片飛行狀態(tài)進行實時測量，將靶板設(shè)計為可安裝測量探針的基板。探針基板為?80mm的LY12鋁，可布局多路多圈的測量探針。飛片與探針基板間的空腔距離為約2mm。為了準確測量和標定，使用DPS測速系統(tǒng)［10］。通過多路多圈的DPS探針布置（見圖3），可以得出飛片直徑40mm內(nèi)的平面性分布。其中，探針1～9從內(nèi)至外分別分布于圓心、?20mm和?40mm，主要是測量飛片速度和評估飛片擊靶的平面性。探針10、11和12分別對應(yīng)炸藥網(wǎng)絡(luò)爆轟的一波點、二波點和四波點，主要測量網(wǎng)絡(luò)起爆點對爆轟波傳播均勻性的影響。在位于?40mm的圓上設(shè)計了4個?6mm的圓孔用于排出飛片飛行時的壓縮氣流，避免對DPS測量的影響。

進行了兩發(fā)典型加載狀態(tài)的實驗測量和標定，實驗參數(shù)與部分結(jié)果見表1，DPS測量結(jié)果見圖4。表中，h1、h2、h3和h4分別為飛片厚度、有機玻璃厚度、45鋼板厚度和空腔距離?？梢园l(fā)現(xiàn)，在2mm空腔距離內(nèi)，飛片已經(jīng)達到一個穩(wěn)定的最大速度飛行狀態(tài)。實驗1、2的飛片速度分別為v1＝（299±4）m／s，v2＝（221±3）m／s。多點測量表明，?40mm內(nèi)飛片的擊靶速度相對誤差在1.5%以內(nèi)，飛片飛行速度均勻性很好。

圖3 DPS探針布局Fig.3 Position of DPS probes

表1 標定實驗的加載狀態(tài)和部分實驗結(jié)果Table1 Experimental parameters and results of the calibrated equipment

圖4 不同測點處飛片后自由面速度歷史Fig.4 Free surface velocity histories for flyers at different points

為了分析飛片的擊靶平面性，對飛片后自由面速度歷史進行積分處理。在空腔的某個面上飛片最早和最晚到達的時間差即為這個面上飛片的擊靶波形差。實驗1、2中當飛片分別到達1.9和1.7mm的空腔距離時，各測點的時間點分布見圖5所示。實驗1、2中飛片擊靶波形差在?20mm范圍內(nèi)分別為146和224ns，到?40mm時，飛片擊靶波形差分別為224和468ns，各個測點的位移差在0.04～0.08mm之間，表明飛片在?40mm的范圍內(nèi)具有較好的平面性。利用最小二乘法對實驗點進行曲線擬合，時間t隨半徑r呈拋物線分布，實驗1、2的t（r）關(guān)系分別為t1（r）＝3.1×10－4r2＋48.795 4、t2（r）＝8.7×10－4r2＋51.519 5，t和r的單位分別為μs和mm。由擊靶波形差可計算得?40mm范圍內(nèi)實驗1、2飛片的擊靶傾角分別為0.189°和0.298°，與一級輕氣炮約0.5°和二級輕氣炮約1°相比［11］，本文的飛片驅(qū)動裝置的擊靶波形完全可以滿足一維平面沖擊加載的要求。

綜上所述，實驗表明：網(wǎng)絡(luò)爆轟驅(qū)動飛片的速度可以在一定范圍內(nèi)可調(diào)，飛片擊靶均勻性和平面性較好，飛片在2mm的空腔距離內(nèi)能夠達到穩(wěn)定飛行狀態(tài)。

圖5 飛片擊靶波形Fig.5 The flyer front surface

將實驗結(jié)果與數(shù)值模擬進行對比。首先，對飛片穩(wěn)定時的平均速度幅值進行比較，實驗1、2的飛片擊靶速度實驗值分別為299和221m／s，相應(yīng)的計算值為302和220m／s，實驗測量與計算的相對誤差小于2%。其次，對飛片飛行過程的速度歷史進行比較，圖6給出了飛片中心測點1的后自由面速度－時間曲線的實驗結(jié)果與計算結(jié)果，兩者符合程度很好。

圖6 后自由面速度剖面的實驗結(jié)果與數(shù)值模擬結(jié)果對比Fig.6 Comparison between experimental and calculated free－surface velocity profile

3 結(jié) 論

基于網(wǎng)絡(luò)爆轟的原理，提出并設(shè)計了一種網(wǎng)絡(luò)爆轟驅(qū)動飛片的加載技術(shù)，并對典型加載裝置進行了實驗測量和標定，結(jié)合實驗和數(shù)值模擬對該技術(shù)進行了系統(tǒng)的分析和研究，得到以下主要結(jié)論：

使用網(wǎng)絡(luò)爆轟驅(qū)動飛片技術(shù)，可以有效降低傳統(tǒng)平面波透鏡加載裝置中的主炸藥裝藥厚度，極大地降低飛片運動速度，使低壓下的材料動態(tài)響應(yīng)特性研究得到新的加載手段；網(wǎng)絡(luò)爆轟驅(qū)動飛片在一定范圍內(nèi)的擊靶均勻性和平面性較好，完全滿足一維應(yīng)變平面沖擊波加載的需求；通過改變衰減層的厚度，飛片速度可以實現(xiàn)在200～350m／s間連續(xù)可調(diào)。

［1］經(jīng)福謙，陳俊祥.動高壓原理與技術(shù)［M］.北京：國防工業(yè)出版社，2006：128－155.

［2］趙峰.炸藥強爆轟驅(qū)動高速金屬飛片的實驗和理論研究［D］.綿陽：中國工程物理研究院，2005.

［3］經(jīng)福謙.實驗物態(tài)方程導引［M］.2版.北京：科學出版社，1999：125－163.

［4］ЦыпкцннВ И.Исспедованиепроцессаразгонамедньгхпродуктамивзръгвачереззазор［J］.ЖурналТехничекий Физики，1975，45（3）：624－629.

［5］趙劍衡，趙峰，文尚剛，等.帶空腔爆轟加載裝置對驅(qū)動飛片的影響［J］.爆炸與沖擊，2001，21（4）：307－310.ZHAO Jian－h(huán)eng，ZHAO Feng，WEN Shang－gang，et al.The effects of explosive device with a cavity on flyer driving［J］.Explosion and Shock Waves，2001，21（4）：307－310.

［6］張萬甲，劉倉理.可調(diào)控靶中應(yīng)力波形的吉帕壓力范圍平面爆炸加載技術(shù)［J］.爆炸與沖擊，2006，26（1）：91－96.ZHANG Wan－jia，LIU Cang－li.A planar explosive loading technique of GPa pressures and capable of adjusting the stress waves in targets［J］.Explosion and Shock Waves，2006，26（1）：91－96.

［7］Kanel G I，Razorenov S V，F(xiàn)ortov V E.Shock－wave phenomena and the properties of condensed matter［M］.New York：Springer，2003：15－17.

［8］金柯，李平，吳強，等.爆轟產(chǎn)物驅(qū)動飛片運動數(shù)值模擬研究［J］.爆炸與沖擊，2004，24（5）：419－424.JIN Ke，LI Ping，WU Qiang，et al.Numerical simulation on acceleration process of flyer driven by detonation product through an air cushion［J］.Explosion and Shock Waves，2004，24（5）：419－424.

［9］張寶坪，張慶明，黃風雷.爆轟物理學［M］.北京：兵器工業(yè)出版社，1997：157－164.

［10］王德田，李澤仁，吳建榮，等.光纖位移干涉儀在爆轟加載飛片速度測量中的應(yīng)用［J］.爆炸與沖擊，2009，29（1）：105－108.WANG De－tian，LI Ze－ren，WU Jian－rong，et al.An optical－fiber displacement interferometer for measuring velocities of explosively－driven metal plates［J］.Explosion and Shock Waves，2009，29（1）：105－108.

［11］王金貴.氣體炮原理與技術(shù)［M］.北京：國防工業(yè)出版社，2001：165－171.