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

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

高能炸藥具有較高的化學(xué)釋放能和較快的反應(yīng)速度，是一種常用的高壓沖擊波能源。在化爆加載中，利用高能炸藥爆轟釋放的化學(xué)能轉(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銅飛片的運(yùn)動過程進(jìn)行了實驗研究，表明隨著空氣隙的增大，飛片速度下降，但空氣隙增加到一定厚度時飛片速度難以再繼續(xù)下降。趙劍衡等［5］結(jié)合實驗和數(shù)值模擬研究了在炸藥透鏡和飛片之間引入空腔對爆轟驅(qū)動飛片的影響，結(jié)果表明引入空腔的爆轟加載裝置可以推動6mm銅飛片達(dá)到約1km／s的速度。張萬甲［6］、G.I.Kanel等［7］設(shè)計了組合飛片作為低壓爆轟加載裝置，實驗結(jié)果初步表明組合飛片作為低壓爆轟加載的原理是可行的，但設(shè)計較復(fù)雜。金柯等［8］對爆轟產(chǎn)物與飛片之間預(yù)留空氣隙和使用組合飛片的方式分別進(jìn)行了數(shù)值模擬研究，結(jié)果表明增加空氣隙或使用組合飛片可以使中等阻抗飛片降至1～2km／s。綜上所述，采用平面波透鏡結(jié)構(gòu)的爆轟加載裝置，很難獲得低于1km／s飛片速度的加載能力。

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

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

1.1 基本原理

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

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

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

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

1.2 裝置設(shè)計

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

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

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

炸藥網(wǎng)絡(luò)板中裝藥厚度0.8mm，邊長95mm。衰減層采用沖擊阻抗較高的45鋼板，厚度可視所需飛片速度調(diào)節(jié)。為使多點(diǎn)爆轟波被改造成為準(zhǔn)一維平面波，經(jīng)過計算取有機(jī)玻璃板厚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)方程，表達(dá)式為［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é)果表明，當(dāng)鋼板厚度從1mm增加到10mm時，飛片的擊靶速度相應(yīng)從330m／s減小到210m／s，并且飛片飛行距離在2mm內(nèi)即可達(dá)到穩(wěn)定的速度平臺。應(yīng)力波在飛片中來回反射可造成飛片自由面速度振蕩；鋼板厚度與飛片速度之間并不完全成線性關(guān)系，隨著鋼板厚度持續(xù)增加，擊靶速度繼續(xù)減小的幅度十分有限。如果使用該飛片撞擊相同材料的無氧銅樣品可以達(dá)到3.8～6.2GPa沖擊應(yīng)力，因此初步計算認(rèn)為該結(jié)構(gòu)的裝置具有可行性，并且通過改變衰減層鋼板的厚度可以使飛片的速度發(fā)生連續(xù)性的變化，得到不同的沖擊應(yīng)力。

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

2 實驗測量與標(biāo)定

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

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

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

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

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

為了分析飛片的擊靶平面性，對飛片后自由面速度歷史進(jìn)行積分處理。在空腔的某個面上飛片最早和最晚到達(dá)的時間差即為這個面上飛片的擊靶波形差。實驗1、2中當(dāng)飛片分別到達(dá)1.9和1.7mm的空腔距離時，各測點(diǎn)的時間點(diǎn)分布見圖5所示。實驗1、2中飛片擊靶波形差在?20mm范圍內(nèi)分別為146和224ns，到?40mm時，飛片擊靶波形差分別為224和468ns，各個測點(diǎn)的位移差在0.04～0.08mm之間，表明飛片在?40mm的范圍內(nèi)具有較好的平面性。利用最小二乘法對實驗點(diǎn)進(jìn)行曲線擬合，時間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)能夠達(dá)到穩(wěn)定飛行狀態(tài)。

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

將實驗結(jié)果與數(shù)值模擬進(jìn)行對比。首先，對飛片穩(wěn)定時的平均速度幅值進(jìn)行比較，實驗1、2的飛片擊靶速度實驗值分別為299和221m／s，相應(yīng)的計算值為302和220m／s，實驗測量與計算的相對誤差小于2%。其次，對飛片飛行過程的速度歷史進(jìn)行比較，圖6給出了飛片中心測點(diǎn)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ù)，并對典型加載裝置進(jìn)行了實驗測量和標(biāo)定，結(jié)合實驗和數(shù)值模擬對該技術(shù)進(jìn)行了系統(tǒng)的分析和研究，得到以下主要結(jié)論：

使用網(wǎng)絡(luò)爆轟驅(qū)動飛片技術(shù)，可以有效降低傳統(tǒng)平面波透鏡加載裝置中的主炸藥裝藥厚度，極大地降低飛片運(yùn)動速度，使低壓下的材料動態(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］趙峰.炸藥強(qiáng)爆轟驅(qū)動高速金屬飛片的實驗和理論研究［D］.綿陽：中國工程物理研究院，2005.

［3］經(jīng)福謙.實驗物態(tài)方程導(dǎo)引［M］.2版.北京：科學(xué)出版社，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］金柯，李平，吳強(qiáng)，等.爆轟產(chǎn)物驅(qū)動飛片運(yùn)動數(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］張寶坪，張慶明，黃風(fēng)雷.爆轟物理學(xué)［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.