張 柱，晉艷娟

（太原科技大學(xué)應(yīng)用科學(xué)學(xué)院，山西 太原030024）

SHPB實(shí)驗(yàn)系統(tǒng)和輕氣炮是目前最常用的動(dòng)態(tài)加載實(shí)驗(yàn)設(shè)備［1］，所能加載的材料試件的尺寸都較小，只能對(duì)符合縮比模型實(shí)驗(yàn)條件的材料進(jìn)行實(shí)驗(yàn)。將兩者應(yīng)用于混凝土和鋼筋混凝土材料的動(dòng)態(tài)力學(xué)性能實(shí)驗(yàn)研究，做縮比模型實(shí)驗(yàn)，材料試件很難滿(mǎn)足等比例縮小和材料相同的條件，所得實(shí)驗(yàn)數(shù)據(jù)的真實(shí)性必然受到影響。要想解決此類(lèi)問(wèn)題，根本的方法是發(fā)展大尺寸試件的實(shí)驗(yàn)加載手段，直接對(duì)真實(shí)混凝土試件進(jìn)行加載，這樣所得的實(shí)驗(yàn)數(shù)據(jù)才可靠。

如果選取炸藥作為驅(qū)動(dòng)源動(dòng)力，一種方法是在保證飛片平面度的前提下，先利用炸藥爆炸驅(qū)動(dòng)大直徑的金屬飛片，再用高速飛片對(duì)真實(shí)混凝土試件進(jìn)行動(dòng)態(tài)加載，完成與輕氣炮相同的加載效果。孫承緯等［2］、趙鋒等［3］利用平面波發(fā)生器實(shí)現(xiàn)了多級(jí)爆轟驅(qū)動(dòng)；欒廣博等［4］和張柱等［5］利用階梯型炸藥驅(qū)動(dòng)錐形金屬飛片，對(duì)大尺寸試件的加載實(shí)驗(yàn)進(jìn)行了初步研究。另一種方法就是利用炸藥直接對(duì)材料試件進(jìn)行加載，但通常只能在材料內(nèi)部產(chǎn)生三角形沖擊波，實(shí)驗(yàn)數(shù)據(jù)很難測(cè)定，且較離散。N.Gebbeken等［6］在點(diǎn)起爆方式下對(duì)真實(shí)混凝土試件實(shí)現(xiàn)了三角形沖擊波加載，得到了不同深度處的壓力－時(shí)間曲線(xiàn)。

在炸藥直接加載情況下，為了得到與輕氣炮加載相同的加載波，本文中提出利用接觸爆炸下的反向起爆模型對(duì)材料試件進(jìn)行加載，在理論上給出沖擊波加載對(duì)應(yīng)的平臺(tái)壓力和平臺(tái)寬度的計(jì)算方法，并結(jié)合數(shù)值模擬對(duì)理論解進(jìn)行驗(yàn)證，以期為爆炸加載沖擊波與輕氣炮加載應(yīng)力波的等效關(guān)系的建立提供理論基礎(chǔ)。

1 爆炸加載沖擊波理論研究

對(duì)于輕氣炮加載實(shí)驗(yàn)而言，理想情況下，飛片與靶板相撞后，會(huì)在材料試件內(nèi)部產(chǎn)生一個(gè)矩形加載應(yīng)力波。炸藥爆炸一般產(chǎn)生的為三角形沖擊波，如果有n個(gè)相同的三角形沖擊波按照到達(dá)的先后順序依次疊加，便可得到帶有平臺(tái)的加載沖擊波，平臺(tái)上帶有明顯的壓力震蕩，如圖1所示。如果相鄰2個(gè)三角形沖擊波的時(shí)間間隔逐漸縮短，甚至連續(xù)，那么疊加所得沖擊波平臺(tái)上的震蕩就會(huì)消失，最終必將獲得與輕氣炮加載相同的矩形沖擊波。

1.1 形成機(jī)理分析

基于圖1所示的疊加原理，采用接觸爆炸反向起爆加載方式，如圖2所示，可以獲得帶平臺(tái)的加載沖擊波。把炸藥看成由無(wú)數(shù)個(gè)矩形微元體組成，分別為1、2、3、4、5、6、7、…。當(dāng)炸藥被起爆后，首先微元體1爆炸，在固體表面施加第1個(gè)三角形沖擊波；然后微元體2爆炸，在固體表面施加第2個(gè)三角形沖擊波；再后微元體3爆炸，在固體表面施加第3個(gè)三角形沖擊波；后續(xù)微元體依次爆炸，產(chǎn)生連續(xù)不斷的三角形沖擊波，依次到達(dá)固體表面。由于這些三角形沖擊波是連續(xù)到達(dá)固體表面的，使固體表面一直維持著平穩(wěn)的壓力供應(yīng)，因此最終必然會(huì)在固體表面獲得帶平臺(tái)的加載沖擊波。

圖1 三角形沖擊波疊加成帶平臺(tái)的沖擊波Fig.1Triangular waves superimposed into shock wave with platform

為驗(yàn)證上述加載沖擊波形成機(jī)理的正確性，利用AUTODYN軟件建立了與圖2相同的反向起爆模型（平面軸對(duì)稱(chēng)問(wèn)題），模型中所用材料均選自軟件自帶的材料庫(kù)，分別為 AIR、COMP B和 AL2024，得到了固體表面的壓力－時(shí)間曲線(xiàn)，如圖3所示。壓力－時(shí)間曲線(xiàn)中存在明顯的加載平臺(tái)，說(shuō)明采用接觸爆炸反向起爆方式得到的加載波形與采用輕氣炮加載得到的加載波形相同。此外，壓力－時(shí)間曲線(xiàn)剛開(kāi)始存在很大的壓力峰值，主要是由起爆瞬間化學(xué)反應(yīng)區(qū)內(nèi)的化學(xué)峰造成的。結(jié)合部分?jǐn)?shù)值模擬發(fā)現(xiàn)，該峰值隨著沖擊波在材料內(nèi)部傳播距離的加大而逐漸衰減。

圖2 平面軸對(duì)稱(chēng)反向起爆模型Fig.2Plane axisymmetric reverse detonation model

圖3 接觸面的壓力－時(shí)間曲線(xiàn)Fig.3Pressure －time curve at contact surface

1.2 沖擊波平臺(tái)壓力的計(jì)算

爆炸加載條件下，固體可以看成流體處理，故稱(chēng)為可壓縮流體。以下理論推導(dǎo)公式涉及的符號(hào)有S、e、p、D、ρ、γ、u、τ和c，分別為熵、內(nèi)能、壓力、沖擊波波速、密度、絕熱指數(shù)、粒子速度、比容和聲速。帶下標(biāo)J的量為爆轟波陣面上的參數(shù)，帶＊的量為固體材料的參數(shù)，帶下標(biāo)0的量為材料的初始參數(shù)。

設(shè)x＜0部分為可壓縮流體，x≥0部分為炸藥柱，炸藥在x＝0處起爆，如圖2所示，此時(shí)有一個(gè)爆轟波在炸藥內(nèi)傳播，流體被壓縮，流體內(nèi)部傳入一個(gè)沖擊波。在爆轟波波后始終跟隨著一個(gè)中心稀疏波，稀疏區(qū)后跟一個(gè)常數(shù)區(qū)，爆炸產(chǎn)物與流體分界面兩邊都是常數(shù)，對(duì)此種情況下的爆炸問(wèn)題，爆炸產(chǎn)物和流體的（u，p）曲線(xiàn)［7］分別為：

在爆轟波陣面上，有：

當(dāng)炸藥和固體材料確定后，對(duì)應(yīng)的各個(gè)材料參數(shù)也就確定，利用逐次逼近法便能求得式（1）確定的2條（u，p）曲線(xiàn)的交點(diǎn)（，），即得到反向起爆模型下可壓縮固體的表面壓力和粒子速度。

顯然，當(dāng)被研究固體材料確定時(shí)，選用不同的炸藥對(duì)應(yīng)不同的平臺(tái)壓力；一旦炸藥和固體材料兩者都確定，平臺(tái)壓力即確定不變，與材料幾何尺寸無(wú)關(guān)。表1中給出了3種常用炸藥（TNT、COMP B和PETN1.77）與3種代表材料（Al、Cu和W）作用后的平臺(tái)壓力和粒子速度。

表1 3種常用炸藥和3種代表材料作用后的平臺(tái)壓力和粒子速度Table 1Shock wave platform pressure and particle velocity for three explosives loading three materials，respectively

1.3 沖擊波平臺(tái)寬度的計(jì)算

從圖2可以看出：在炸藥反向起爆過(guò)程中加載沖擊波的平臺(tái)寬度取決于AD、DE、CF和EF等4條邊上的稀疏波作用。對(duì)于B點(diǎn)的壓力來(lái)說(shuō)，主要是看AD、DE、CF和EF等4條邊上的哪個(gè)反射稀疏波先到達(dá)B點(diǎn)，一旦有稀疏波到達(dá)，壓力平臺(tái)將不能維持，壓力開(kāi)始下降。

設(shè)AD＝BE＝CF＝a1，BC＝EF＝a2，AB＝DE＝a3，取爆轟產(chǎn)物的絕熱指數(shù)γ＝3，結(jié)合波系分區(qū)，爆轟產(chǎn)物包括常數(shù)區(qū)和中心稀疏區(qū)，爆轟產(chǎn)物中［8］：

如果不考慮爆炸產(chǎn)物與流體分界面的運(yùn)動(dòng)，將其看成固壁，此時(shí)，常數(shù)區(qū)中u＝0，c＝DJ/2；中心稀疏區(qū)中，DJ/2≤x/t≤DJ，DJ/2≤c≤3DJ/4。顯然，c＝3DJ/4對(duì)應(yīng)的是爆轟波陣面上的聲速cJ，其值在爆轟產(chǎn)物中最大。

（1）炸藥側(cè)向（EF）稀疏波沿爆轟波陣面上的傳播速度最大，其到達(dá)B點(diǎn)所需時(shí)間：

（2）炸藥的后自由面（CF）稀疏波傳到B點(diǎn)所需的時(shí)間包括ta和tb兩部分，ta為爆炸沖擊波從BE傳到CF所用的時(shí)間，tb為反射稀疏波從CF傳回BE所用的時(shí)間，且有：

當(dāng)爆炸沖擊波傳到CF時(shí)，其后區(qū)域均為爆轟產(chǎn)物，分別為常數(shù)區(qū)和中心稀疏區(qū)，2個(gè)區(qū)對(duì)應(yīng)的長(zhǎng)度均為a2/2。故tb可分2段計(jì)算，tb1為反射稀疏波在常數(shù)區(qū)中的傳播時(shí)間，tb2為反射稀疏波在中心稀疏區(qū)中的傳播時(shí)間，則：

（3）固體中側(cè)向（DE）稀疏波到達(dá)B點(diǎn)所需的時(shí)間：

（4）固體中前自由面（AD）稀疏波傳到B點(diǎn)所需的時(shí)間：

最終，爆炸沖擊波的壓力平臺(tái)寬度：

由式（9）可知，沖擊波平臺(tái)寬度與加載模型幾何尺寸、炸藥的沖擊波波速和固體材料中的聲速有關(guān)。由式（4）和式（6）可知，確定炸藥的幾何尺寸時(shí)，應(yīng)盡量使t1＝t2，即a1/a2＝21/10，此時(shí)兩側(cè)稀疏波同時(shí)達(dá)到，炸藥利用率最高；當(dāng)a1/a2＞21/10時(shí)，炸藥側(cè)向稀疏干擾為主；當(dāng)a1/a2＜21/10時(shí)，炸藥后自由面稀疏干擾為主，此時(shí)固定a1不變，改變a2，可得到不同的沖擊波平臺(tái)寬度，使之與輕氣炮加載波形寬度保持相同。確定固體材料試件的尺寸時(shí)除了考慮t3、t4外，更應(yīng)該考慮輕氣炮加載實(shí)驗(yàn)中樣品的設(shè)計(jì)要求，滿(mǎn)足寬厚比［9］，確保數(shù)據(jù)采集系統(tǒng)能得到有效的實(shí)驗(yàn)結(jié)果。表2以COMP B炸藥和固體材料Al為例，給出了幾種尺寸條件下平臺(tái)寬度的理論計(jì)算值。

表2 COMP B炸藥反向起爆時(shí)，不同尺寸對(duì)應(yīng)的理論平臺(tái)寬度Table 2Theoretical platform width for several sizes in reverse detonation of COMP B

2 反向起爆的數(shù)值模擬及結(jié)果分析

為驗(yàn)證上述理論分析的準(zhǔn)確性，借助AUTODYN軟件對(duì)反向起爆加載模型進(jìn)行了數(shù)值模擬。

2.1 沖擊波平臺(tái)壓力的模擬

討論沖擊波加載平臺(tái)壓力時(shí)，保持?jǐn)?shù)值分析模型的幾何尺寸不變，取a1＝100mm，a2＝100mm，a3＝50mm，不斷改變所用炸藥和固體材料（與表1對(duì)應(yīng)），最終模擬結(jié)果如表3所示。

表3 平臺(tái)壓力和粒子速度的數(shù)值模擬結(jié)果Table 3Numerical simulation results of shock wave platform pressure and particle velocity

觀察表3可知：

（1）無(wú)論是理論值還是數(shù)值模擬值，平臺(tái)壓力和粒子速度的變化規(guī)律均符合客觀事實(shí)，即固體材料相同時(shí)，隨著炸藥密度的增加，炸藥威力變大，平臺(tái)壓力和粒子速度均逐漸變大；炸藥相同時(shí)，隨著固體材料密度的增加，材料可壓縮性變差，粒子速度逐漸變小。

（2）如果以數(shù)值模擬值為準(zhǔn)確值，對(duì)于W，隨著炸藥密度的增加，其相對(duì)誤差逐漸變小。因?yàn)閃本身強(qiáng)度和密度很大，不容易被壓縮，當(dāng)炸藥密度變大時(shí)，對(duì)W的壓力作用也變大，此時(shí)W更近似流體；對(duì)于Cu和Al，理論值與數(shù)值模擬值的相對(duì)誤差都在5%以?xún)?nèi)，兩者基本一致。說(shuō)明平臺(tái)壓力對(duì)應(yīng)的理論推導(dǎo)比較適用于較軟材料的計(jì)算，對(duì)于高強(qiáng)度材料，應(yīng)該盡可能選用高性能炸藥進(jìn)行驅(qū)動(dòng)。

2.2 沖擊波平臺(tái)寬度的模擬

討論沖擊波加載平臺(tái)寬度時(shí)，選定COMP B對(duì)Al作用，研究炸藥的幾何尺寸對(duì)其的影響，具體模擬結(jié)果見(jiàn)表4，表中t′e為沖擊波平臺(tái)寬度的數(shù)值模擬結(jié)果。對(duì)于固體材料幾何尺寸可能引起的影響，不在此處討論，因?yàn)閷?duì)于高爆速炸藥來(lái)說(shuō)，當(dāng)固體材料的幾何尺寸滿(mǎn)足輕氣炮加載實(shí)驗(yàn)樣品的設(shè)計(jì)要求時(shí)（寬厚比約為2），其側(cè)面和自由面的稀疏波到達(dá)分界面中心的時(shí)間都較長(zhǎng)，即t3，t4＞t1，t2。

表4 沖擊波平臺(tái)寬度的數(shù)值模擬結(jié)果Table 4Numerical simulation results of shock wave platform width

觀察表4可知：

（1）理論值與模擬值基本一致，說(shuō)明建立的理論推導(dǎo)公式可以用于實(shí)驗(yàn)前的初步估算。

（2）相對(duì)誤差隨著裝藥量的提高，不斷減小，說(shuō)明建立的理論推導(dǎo)方法更適合大藥量的沖擊波加載。

（3）當(dāng)a1/a2＜21/10時(shí)，若固定a1不變，改變a2便可得到不同的沖擊波平臺(tái)寬度。

（4）對(duì)于炸藥來(lái)說(shuō)，由側(cè)向稀疏波決定（a1）的平臺(tái)寬度（表4中序號(hào)5～10）與由后自由面稀疏波決定（a2）的平臺(tái)寬度（表4中序號(hào)1～4）相比，前者的理論解更接近數(shù)值解，原因是固壁假設(shè)和討論后自

由面稀疏波傳播時(shí)間時(shí)tb2中所用的聲速為一段平均，若選為2段或更多段平均計(jì)算，理論解一定更精確。

3 結(jié) 論

為獲得與輕氣炮相同的加載效果，并突破加載試件的尺寸限制，完成對(duì)真實(shí)混凝土或鋼筋混凝土的動(dòng)態(tài)加載研究，基于疊加原理，提出利用炸藥反向起爆模型實(shí)現(xiàn)對(duì)可壓縮固體材料的沖擊波加載。理論上，通過(guò)聯(lián)立爆炸產(chǎn)物和可壓縮流體的（u，p）曲線(xiàn)，確定了沖擊波平臺(tái)壓力；綜合考慮炸藥與材料各自由邊的反射稀疏波情況，給出沖擊波加載平臺(tái)寬度的計(jì)算方法。并結(jié)合數(shù)值模擬，對(duì)理論分析結(jié)果進(jìn)行了驗(yàn)證，兩者基本一致，說(shuō)明理論分析方法正確有效。如果要建立爆炸加載沖擊波與輕氣炮加載應(yīng)力波的等效關(guān)系，只需對(duì)比相應(yīng)的理論解。

本文在完成過(guò)程中得到北京理工大學(xué)爆炸科學(xué)與技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室寧建國(guó)教授和趙衡陽(yáng)教授的很多指導(dǎo)，在此表示感謝。

［1］ 寧建國(guó)，王成，馬天寶.爆炸與沖擊動(dòng)力學(xué)［M］.北京：國(guó)防工業(yè)出版社，2010：398.

［2］ 孫承緯，文尚剛，趙鋒.多級(jí)炸藥爆轟高速驅(qū)動(dòng)技術(shù)的 Gurney模型優(yōu)化分析［J］.爆炸與沖擊，2004，24（4）：299 －304.Sun Cheng －wei，Wen Shang－gang，Zhao Feng.An optimal analysis of multi －stage explosive accelerated high velocity flyers with the improved Gurney model［J］.Explosion and Shock Waves，2004，24（4）：299 －304.

［3］ 趙鋒，文尚剛，孫承緯，等.多級(jí)串聯(lián)式超高速飛片裝置實(shí)驗(yàn)研究［J］.爆炸與沖擊，2001，21（1）：13 －16.Zhao Feng，Wen Shang－gang，Sun Cheng －wei，et al.Experimental studies on multiple －stage flyers［J］.Explosion and Shock Waves，2001，21（1）：13 －16.

［4］ 欒廣博，郝莉，張柱.平面爆轟波驅(qū)動(dòng)金屬飛片運(yùn)動(dòng)的數(shù)值模擬研究［J］.高壓物理學(xué)報(bào)，2011，25（5）：451 －456.Luan Guang －bo，Hao Li，Zhang Zhu.Numerical study on the acceleration of metallic flyers driven by planar detonation wave of explosives［J］.Chinese Journal of High Pressure Physics，2011，25（5）：451 －456.

［5］ 張柱，欒廣博，寧建國(guó).爆炸載荷作用下錐形飛片的平面度研究［J］.高壓物理學(xué)報(bào)，2011，25（3）：221 －226.Zhang Zhu，Luan Guang －bo，Ning Jian －guo.Flatness study on the metallic flyer driven by explosive loading equipment［J］.Chinese Journal of High Pressure Physics，2011，25（3）：221 －226.

［6］ Gebbeken N，Greulich S，Pietzsch A.Hugoniot properties for concrete determined by full －scale detonation experiments and flyer －plate －impact test［J］.International Journal of Impact Engineering，2006，32（12）：2017 －2031.

［7］ 周毓麟.一維非定常流體力學(xué)［M］.北京：科學(xué)出版社，1998：320.

［8］ 李維新.一維不定常流與沖擊波［M］.北京：國(guó)防工業(yè)出版社，2003：408.

［9］ 經(jīng)福謙.實(shí)驗(yàn)物態(tài)方程導(dǎo)引［M］.2版.北京：科學(xué)出版社，1999：212.