金屬約束下定常非理想爆轟的理論研究＊

于 明，孫宇濤，張文宏

（北京應(yīng)用物理與計算數(shù)學(xué)研究所，北京100094）

高能炸藥非理想爆轟問題，是爆轟物理學(xué)領(lǐng)域非常重要的研究內(nèi)容，關(guān)系到對炸藥直徑效應(yīng)、熄爆速度、約束效應(yīng)以及驅(qū)動效率等爆轟規(guī)律的認識，特別是鈍感炸藥（insensitive high explosive，IHE）的非理想爆轟更受關(guān)注。與敏感炸藥爆轟相比，鈍感炸藥爆轟具有2個特點：（1）爆速更低；（2）化學(xué)反應(yīng)區(qū)更寬，因此鈍感炸藥的非理想爆轟現(xiàn)象更顯著。

對炸藥非理想爆轟問題的理論研究開展很早，主要理論有2類［1］：一類是流管方法，另一類是彎曲陣面法。流管方法以F.Wecken等［2］的工作為主，假定前導(dǎo)波陣面是平面，波后產(chǎn)物流線圍成的流管截面是逐漸擴大的，發(fā)散流動被認為是一維的。彎曲陣面法以J.B.Bdzil等［3］的工作為主，假定前導(dǎo)波陣面形狀為已知曲面（如一定半徑的球面等），波后流動是二維的。需指出的是，J.B.Bdzil等［3］的工作是W.W.Wood等［4］的軸線一維理論成果的推廣，初步揭示了二維反應(yīng)區(qū)結(jié)構(gòu)的特征，但在計算聲速面時面臨較大困難。近年來，出現(xiàn)更多的理論方法試圖給出較精確的反應(yīng)區(qū)結(jié)構(gòu)特征，如W.B.Brown等［5］的最大熵增理論和S.D.Watt等［6］的流線理論等。最大熵增理論指出，對定常、自持爆轟，在流線上前導(dǎo)沖擊波與聲速線之間的熵增最大，由此可由變分方法確定前導(dǎo)波陣面后聲速線的位置。流線理論假定波后流線是發(fā)散的直線，給定爆轟速度的徑向變化關(guān)系后沿所有流線求解可以給出波陣面形狀與聲速線位置。以上這些理論方法僅考察了無約束或弱約束條件下的非理想爆轟流動，這些條件下聲速線與前導(dǎo)波陣面在炸藥邊緣處交于一點，而金屬約束（強約束）條件下聲速線與前導(dǎo)波陣面并不相交，并且這些理論方法只給出聲速線的位置，而未給出化學(xué)反應(yīng)結(jié)束線的位置，而且，這些理論方法對未反應(yīng)炸藥和爆轟產(chǎn)物均采用較簡單的多方指數(shù)狀態(tài)方程。

本文中對鋼和鋁約束條件下鈍感炸藥PBX9502的定常非理想爆轟狀態(tài)進行理論研究。理論分析可知，在惰性金屬施加的強約束條件下，爆轟流場中聲速線與前導(dǎo)沖擊波陣面的相對位置同金屬壓縮性相關(guān)，這一點與無約束或弱約束條件下的非理想爆轟流動顯著不同。當(dāng)金屬對炸藥約束作用時，爆轟前導(dǎo)沖擊波折射進入金屬，爆轟沖擊邊緣角由未反應(yīng)炸藥的沖擊波極曲線和金屬的沖擊波極曲線的匹配關(guān)系確定。為了更準(zhǔn)確地分析，對未反應(yīng)炸藥和爆轟產(chǎn)物均采用JWL形式狀態(tài)方程，對金屬采用p（r，T）形式的狀態(tài)方程［7］，炸藥化學(xué)反應(yīng)模型中采用三項式Lee－Tarver反應(yīng)率［8］。采用過爆轟前導(dǎo)沖擊波的流線偏轉(zhuǎn)角在影響域內(nèi)沿高度線性變化的假設(shè)，則由斜沖擊波極曲線關(guān)系式可以求出前導(dǎo)沖擊波陣面的形狀。采用前導(dǎo)沖擊波陣面之后的流場中的流線是直線的假設(shè)，則爆轟流動控制方程由偏微分方程變?yōu)檠亓骶€的常微分方程，沿著所有流線求解便給出聲速線與化學(xué)反應(yīng)結(jié)束線的位置。理論分析同時給出約束金屬折射沖擊波后面的流體的流動狀態(tài)。

1 爆轟前導(dǎo)沖擊波陣面形狀的求解

金屬約束條件下高能炸藥的定常非理想爆轟流動狀態(tài)如圖1所示。如果把參考坐標(biāo)放在爆轟前導(dǎo)沖擊波陣面OA與金屬界面的交點O點上，則該交點附近的流動被定?；?。金屬與炸藥爆轟波相互作用時，爆轟受邊界側(cè)向膨脹的影響呈非理想狀態(tài)，因此炸藥－金屬界面附近的前導(dǎo)沖擊波陣面變彎曲，如果炸藥厚度Y足夠大，變彎曲的前導(dǎo)沖擊波高度被認為是L（圖中OE段），則大于L處的前導(dǎo)沖擊波陣面仍是平面，前導(dǎo)沖擊波后的聲速線為BC，化學(xué)反應(yīng)結(jié)束線為BD。這時前導(dǎo)沖擊波前的流動被認為是速度大小為DCJ的均勻來流。爆轟波與金屬相互作用時，首先是前導(dǎo)沖擊波在金屬中發(fā)生折射，因此將爆轟波極曲線建立于前導(dǎo)沖擊波上。前導(dǎo)沖擊波的沖擊狀態(tài)由未反應(yīng)炸藥物性決定，因此給定未反應(yīng)炸藥狀態(tài)方程和波前來流速度，可以由前導(dǎo)沖擊波的極曲線關(guān)系得到過前導(dǎo)沖擊波的流動偏轉(zhuǎn)角θ。在邊界處，由于前導(dǎo)沖擊波后面的流動是亞聲速的，因此前導(dǎo)沖擊波在炸藥－金屬界面上不會產(chǎn)生任何反射波，這時爆轟沖擊邊緣角θe（前導(dǎo)沖擊波陣面與金屬界面在他們交點處的夾角）由未反應(yīng)炸藥前導(dǎo)沖擊波極曲線與金屬折射沖擊波極曲線決定，并且這2條沖擊波極曲線有共同的初始點（初始壓力和初始流線轉(zhuǎn)折角均為0）。

圖1 金屬約束條件下高能炸藥的定常非理想爆轟流動狀態(tài)Fig.1 Steady non－ideal flow of detonation in high explosives confined by metal

為了求出爆轟前導(dǎo)沖擊波陣面的幾何形狀，假定前導(dǎo)沖擊波陣面由函數(shù)xf＝xf（yf）來描述，并用α表示波陣面與波前速度之間的夾角，則有沖擊波陣面函數(shù)斜率關(guān)系式

考慮到金屬約束時爆轟沖擊邊緣角較小，這里采用過前導(dǎo)沖擊波的流線偏轉(zhuǎn)角在邊界影響域內(nèi)沿高度線性變化的假設(shè)，即在沖擊波陣面OE段上有如下關(guān)系式

對未反應(yīng)炸藥，采用JWL形式的狀態(tài)方程，其斜沖擊波極曲線關(guān)系式為

對金屬，采用p（r，T）形式的狀態(tài)方程，有如下表達式

上述3項分別表示冷壓（冷能）、原子核熱壓（原子核熱能）和電子熱壓（電子熱能），每項具體表達式見文獻［7］，由此有其折射沖擊波極曲線關(guān)系式

式（2）中的爆轟沖擊邊緣角θe由未反應(yīng)炸藥前導(dǎo)沖擊波極曲線（3）和金屬折射沖擊波極曲線（4）的交點決定。圖2中給出了鋼和鋁約束PBX9502炸藥時未反應(yīng)炸藥與金屬的沖擊波極曲線匹配關(guān)系，實線為未反應(yīng)炸藥極曲線，虛線為金屬極曲線，點“Cm”為金屬聲速點狀態(tài)。由圖中極曲線交點可知，鋼約束時θe＝4.20°，鋁約束時θe＝2.28°。

從上述思路可見，確定爆轟沖擊邊緣角以后，前導(dǎo)沖擊波后任意位置處的流動偏轉(zhuǎn)角θ由式（2）給出，并由式（3）迭代求出前導(dǎo)沖擊波后的壓力p和相對體積v－，進而可求出波陣面與波前速度之間的夾角α，最終由式（1）獲得爆轟前導(dǎo)沖擊波陣面的幾何形狀。

圖2 未反應(yīng)炸藥與金屬的沖擊波極曲線Fig.2 Shock polar match of unreacted explosive and metals

2 爆轟化學(xué)反應(yīng)區(qū)聲速線和化學(xué)反應(yīng)結(jié)束線的求解

爆轟化學(xué)反應(yīng)區(qū)流場可由如下計化學(xué)反應(yīng)的二維可壓縮定常無粘流動方程來描述

式中：F為爆轟產(chǎn)物的質(zhì)量分?jǐn)?shù)，Rρ，（p，F(xiàn)）為三項式Lee－Tarver反應(yīng)率［8］。

為方便求解，將方程組（5）寫成如下矢量－矩陣形式

式中：q＝［ρ，u，w，F(xiàn)］T，s（q）＝［0，0，0，R／u］T，且有通量矩陣H（q）如下

從偏微分?jǐn)?shù)學(xué)理論可知，式（6）的通量矩陣有4個特征值：2個為實數(shù)，2個為復(fù)數(shù)，故它為雙曲－橢圓型方程。對于實數(shù)特征值，當(dāng)dy／dx＝w／u時，在特征線上成立

這樣，偏微分方程（5）變成為沿特征線方向的常微分方程式（7）。

這里采用爆轟前導(dǎo)沖擊波陣面之后的流場中的流線是直線的假設(shè)，則對前導(dǎo)沖擊波陣面上的任一點以及此點處的流線偏轉(zhuǎn)角θ，有過該點的流線方程為

由定常流體力學(xué)理論知，定常流動的流線就是特征線，故可得

聯(lián)立式（7）和（9），即沿流線方向求解，可獲得任意x位置處的流動狀態(tài)。

特別地，在u2＋w2＝c2F處，且cF為保持爆轟產(chǎn)物質(zhì)量分?jǐn)?shù)不變時的凍結(jié)聲速，此處即為聲速線的位置。在F＝1處，此處即為化學(xué)反應(yīng)結(jié)束線的位置。

Lee－Tarver反應(yīng)模型要求對化學(xué)反應(yīng)區(qū)采用溫度和壓力平衡以及體積和能量可加的混合原則。在此原則下凍結(jié)聲速的計算式推導(dǎo)如下。由化學(xué)反應(yīng)區(qū)混合原則有

式中：下標(biāo)s表示未反應(yīng)炸藥，g表示氣體產(chǎn)物，T為溫度，Q為單位質(zhì)量炸藥化學(xué)反應(yīng)釋放的能量。

則凍結(jié)聲速可以寫成如下表達式

式（10）和（11）聯(lián)立起來經(jīng)過迭代可以消去中間變量vs從而得到cF。

3 金屬折射沖擊波后的流動狀態(tài)

金屬采用p（r，T）形式的狀態(tài)方程時折射沖擊波后等熵聲速表達式為

由式（4）和式（12）聯(lián)立起來可以確定折射沖擊波后流動是超聲速還是亞聲速。

4 鋼和鋁約束條件下鈍感高能炸藥的結(jié)果分析

PBX9502炸藥的化學(xué)反應(yīng)率參數(shù)見文獻［8］，在金屬約束條件下定常傳播的爆轟波受邊界影響的前導(dǎo)沖擊波高度L可由文獻［1］中相關(guān)公式估計。理論分析同高精度的二維Lagrange流體力學(xué)方法的數(shù)值模擬進行了比較。數(shù)值模擬時，炸藥厚度為2.0cm，金屬厚度為0.5cm，長度均為6.0cm；離散網(wǎng)格劃分為250cm－1，總網(wǎng)格數(shù)約為1 000 000。

圖3中給出了獲得的鋼和鋁約束下物質(zhì)界面附近的爆轟前導(dǎo)沖擊波陣面形狀以及聲速線和化學(xué)反應(yīng)結(jié)束線的位置，其中實線為理論計算結(jié)果，虛線為數(shù)值模擬結(jié)果。從比較結(jié)果可以看出，在2種金屬約束條件下理論結(jié)果同數(shù)值結(jié)果符合得較好。爆轟前導(dǎo)沖擊波陣面形狀符合得好，說明采用的流線偏轉(zhuǎn)角是線性變化的假設(shè)是合理的。聲速線和化學(xué)反應(yīng)結(jié)束線符合得好，說明采用前導(dǎo)沖擊波陣面之后的流場中的流線是直線的假設(shè)是合理的。

從圖中結(jié)果也可以看出，在金屬界面附近，在流動達到聲速以后經(jīng)過一定距離化學(xué)反應(yīng)才結(jié)束，這正是非理想爆轟的特點。在金屬約束條件下，爆轟化學(xué)反應(yīng)區(qū)流場中聲速線與前導(dǎo)沖擊波陣面不相交，這一點與無約束或弱約束條件下的非理想爆轟流動顯著不同。

從圖中結(jié)果還可以看出，由于鋼的聲阻抗（3.658g／（cm2·ms））較大而難于壓縮，因此鋼折射沖擊波后面的流動是超聲速；相反地，鋁的聲阻抗（1.495g／（cm2·ms））較小而易于壓縮，因此鋁折射沖擊波后面的流動是亞聲速。這與圖2中極曲線理論結(jié)果一致，圖2（a）中2曲線交點位于金屬聲速點下方，故流動為超聲速；圖2（b）中2曲線交點位于金屬聲速點上方，故流動為亞聲速?？梢娂s束金屬可壓縮性的不同導(dǎo)致折射沖擊波后流動狀態(tài)不同。

圖3 金屬約束條件下非理想爆轟流動特征狀態(tài)Fig.3 Characteristic flow states of non－ideal detonation under metal confinement

5 結(jié) 論

本文中提出了一種研究高能炸藥定常非理想爆轟的理論方法。該方法采用了過爆轟前導(dǎo)沖擊波的流線偏轉(zhuǎn)角在影響域內(nèi)沿高度線性變化以及前導(dǎo)沖擊波陣面之后的流線是直線這2點假設(shè)，并對未反應(yīng)炸藥和爆轟產(chǎn)物采用JWL形式狀態(tài)方程，對金屬采用p（ρ，T）形式的狀態(tài)方程。本理論方法可以解析地求出炸藥前導(dǎo)沖擊波陣面的形狀，以及給出前導(dǎo)沖擊波后聲速線與化學(xué)反應(yīng)結(jié)束線的位置，同時給出約束金屬內(nèi)的流動狀態(tài)：當(dāng)金屬聲阻抗較大時，金屬折射沖擊波后面的流動通常為超聲速，反之則為亞聲速。從本理論方法獲得的爆轟化學(xué)反應(yīng)區(qū)結(jié)構(gòu)特征和約束金屬內(nèi)的流動特征這些典型結(jié)果看，與高精度Lagrange流體力學(xué)數(shù)值方法的模擬結(jié)果符合良好，說明本文的理論方法具有良好的合理性和適用性。同時，本理論方法的花費遠小于高精度Lagrange流體力學(xué)數(shù)值方法，說明本文的理論方法具有良好的高效性。

［1］孫承緯，衛(wèi)玉章，周之奎.應(yīng)用爆轟物理［M］.北京：國防工業(yè)出版社，2000.

［2］Wecken F.Non－ideal detonation with constant lateral expansion［C］∥Proceedings of the 4th International Symposium on Detonation.Silver Spring，Maryland，1965：107－116.

［3］Bdzil J B，Stewart D S.Modeling two－dimensional detonations with detonation shock dynamics［J］.Physics of Fluids A：Fluid Dynamics，1989，1（7）：1261－1267.

［4］Wood W W，Kirkwood J G.Diameter effect in condensed explosives：The relation between velocity and radius of curvature of the detonation wave［J］.Journal of Chemical Physics，1954，22（11）：1920－1924.

［5］Brown W B.Thermodynamics theory of non－ideal detonation and failure［C］∥Proceedings of the 12th International Symposium on Detonation.San Diego，California，2002：718－723.

［6］Watt S D，Sharpe G J.A streamline approach to steady non－ideal detonation theory［C］∥Proceedings of the 14th International Symposium on Detonation.Coeur d′Alene，Idaho，2010.

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

［8］Tarver C M，McGuire E M.Reactive flow modeling of the interaction of TATB detonation waves with inert materials［C］∥Proceedings of the 12th International Symposium on Detonation.San Diego，California，2002：641－649.