崔金雷，侯俊亮，肖輝朗，王修益，曲利峰，何小旭

(1. 四川航天系統(tǒng)工程研究所，四川 成都 610100; 2. 四川航天電子設(shè)備研究所，四川 成都 610100;3. 陸軍裝備部駐成都第三軍代室，四川 成都 610000)

0 引 言

破片是殺爆戰(zhàn)斗部的主要毀傷元，破片的特征參數(shù)有破片質(zhì)量、速度和飛散角等。破片的空間飛散特性決定破片是否作用于目標(biāo)，而速度則決定破片能否對目標(biāo)產(chǎn)生一定程度的殺傷。作為戰(zhàn)斗部殺傷威力的重要參數(shù)，破片速度的準(zhǔn)確計算對戰(zhàn)斗部威力的評估具有相當(dāng)重要的意義。

許多學(xué)者以經(jīng)典的能量守恒方程為基礎(chǔ)提出了破片驅(qū)動理論，如[1]：

式中：C——裝藥質(zhì)量；

Ec——裝藥能量密度；

D——裝藥爆速；

r0——裝藥初始半徑；

rf——殼體破裂半徑；

Ms——殼體質(zhì)量；

Ap——其斷裂比能；

Pf——殼體破裂時爆轟產(chǎn)物壓力；

Vf——殼體破裂時爆轟產(chǎn)物的體積；

β=C/M，M——被驅(qū)動物質(zhì)量。

彈藥工程專業(yè)學(xué)者通過對以上公式實用改進獲得了諸多解析公式，如最著名、應(yīng)用最廣的Gurney公式[2]：

其中Eg為Gurney能。式中只考慮了裝藥質(zhì)量、性能和被驅(qū)動質(zhì)量的影響，獲得的往往是更符合自然破片驅(qū)動的極限速度。應(yīng)用計算不同結(jié)構(gòu)的戰(zhàn)斗部破片初速時，其精度越來越不能滿足威力評估的需要。尤其是對預(yù)制破片戰(zhàn)斗部破片初速計算時，由于內(nèi)襯破裂和預(yù)制破片排布間隙的存在，爆轟產(chǎn)物提前泄露，計算與試驗結(jié)果有較大誤差，見表1。

表1 經(jīng)典公式計算與試驗結(jié)果對比

陳興旺等基于 Gurney 假定和相鄰層顆粒之間力和力的波動量等概率傳遞假定，忽略排列方式引起的孔隙率變化，建立了多層破片初速場計算，其研究過程中忽略了孔隙率的影響[7]。隨著計算機技術(shù)的發(fā)展，數(shù)值計算逐漸成為研究破片驅(qū)動和飛散的重要手段，苗春壯等利用三維有限元軟件ANSYS/LS-DYNA對不同形狀預(yù)控破片驅(qū)動過程進行了數(shù)值計算，獲得了破片初速與破片形狀的影響規(guī)律[8]。凌琦采用數(shù)值計算方法對定向破片的驅(qū)動過程進行了仿真，獲得了破片加載過程、初速和空間分布情況[9]。

本文通過理論和量綱分析，對影響破片初速的內(nèi)襯材料特性、破片形狀、破片尺寸等與爆轟產(chǎn)物泄漏相關(guān)的參數(shù)進行了分析，采用數(shù)值計算和試驗相結(jié)合的方法對不同工況預(yù)制破片戰(zhàn)斗部的爆炸驅(qū)動過程進行了研究，以Gurney公式為基礎(chǔ)構(gòu)建出考慮爆轟氣體泄漏因子的預(yù)制破片初速修正模型。該公式可為預(yù)制破片初速準(zhǔn)確計算和戰(zhàn)斗部毀傷評估提供參考依據(jù)。

1 理論模型

在Gurney公式基礎(chǔ)上，本文針對柱形戰(zhàn)斗部提出一種預(yù)制破片初速修正計算模型，其他外形戰(zhàn)斗部可采取軸向分段近似參考計算：

式中：fd——起爆方式修正系數(shù)；

fe——端面效應(yīng)修正系數(shù)；

fl——氣體泄露系數(shù)；

fh——爆轟產(chǎn)物密度及速度分布修正系數(shù)。

有些學(xué)者對其中的系數(shù)已經(jīng)進行了修正，如蔣建偉、馮順山等對起爆方式對破片的影響進行了試驗修正擬合出了單端起爆、兩端起爆和軸線起爆方式下的fd表達式[10]，Charron對端面效應(yīng)進行了大量試驗研究，給出了fe的表達式[11]。fh系數(shù)的研究需要高精尖的儀器設(shè)備和大量的試驗經(jīng)費支撐，本文將研究重點著眼于fl，即戰(zhàn)斗部破裂造成的氣體泄露對破片驅(qū)動的影響規(guī)律，研究內(nèi)容包括不同的預(yù)制破片類型和尺寸，不同的內(nèi)襯材料等。

氣體泄露系數(shù)為純數(shù)，經(jīng)分析形成fl的函數(shù)關(guān)系如下：

式中：re——裝藥半徑；

ρt——內(nèi)襯密度；

σs——內(nèi)襯強度；

Ap——內(nèi)襯斷裂比能；

δ——內(nèi)襯厚度；

λ——破片尺寸。

函數(shù)中有6個獨立的有量綱物理量，選擇為re、ρt、σs為基本量，根據(jù)Π定理，可組成3個無量綱組合，可得到氣體泄露系數(shù)的相似準(zhǔn)數(shù)。

假定內(nèi)襯破裂半徑僅僅與內(nèi)襯材料相關(guān)，忽略厚度的影響。以此為基礎(chǔ)，式（5）可簡化為：

令x1=Ap/σs，x2=λ/re，從式（6）可見，x1表征材料延展性，數(shù)值越小延展性越好，爆轟產(chǎn)物對破片的加速就越充分；而x2則表征破片周向排布情況，數(shù)值越小排布越緊密爆轟產(chǎn)物泄漏情況越輕微，破片獲得的能量也越多。

2 數(shù)值計算

2.1 計算模型

選取柱形戰(zhàn)斗部為研究目標(biāo)，利用AUTODYN-3D對三種內(nèi)襯材料的預(yù)制破片戰(zhàn)斗部作用過程進行數(shù)值模擬，以獲得不同形狀、不同尺寸的預(yù)制破片在爆轟產(chǎn)物驅(qū)動下的初速規(guī)律。為氣體泄露系數(shù)fl的數(shù)值擬合提供數(shù)據(jù)支撐。

柱形戰(zhàn)斗部一般由主裝藥、前后端板、蒙皮、內(nèi)襯、預(yù)制破片及灌封膠等組成，數(shù)值模擬模型示意見圖1。主裝藥為Ф36 mm的TNT藥柱，前后端板均為5 mm厚的45#鋼板，蒙皮為鋁或非金屬，內(nèi)襯為1100鋁、45#鋼或聚氨酯泡沫材料，灌封膠選用nylon材料，而預(yù)制破片選用30CrMnSi材料，破片特征尺寸選取2～6 mm，戰(zhàn)斗部裝藥由前端面中心點起爆。

圖1 數(shù)值計算模型

根據(jù)幾何縮比，數(shù)值模擬模型中蒙皮厚度小于0.5 mm?？紤]到計算蒙皮厚度方向尺寸小，網(wǎng)格劃分時將導(dǎo)致計算規(guī)模的大幅度增加，計算中將蒙皮轉(zhuǎn)化為內(nèi)襯一并建模以提高效率?？紤]到結(jié)構(gòu)的對稱性，計算模型采用1/4模型。

建立網(wǎng)格模型過程中，綜合考慮計算精度和效率，尺寸規(guī)定如下：模擬主裝藥Euler區(qū)軸向徑向網(wǎng)格尺寸1 mm，前后端板網(wǎng)格尺寸1.66 mm，破片特征尺寸(直徑、邊長)方向網(wǎng)格數(shù)6，內(nèi)襯厚度方向網(wǎng)格數(shù)7，軸向網(wǎng)格數(shù)240。

空氣域采用Euler算法，除對稱面外，各邊界設(shè)置flow_out邊界條件，實現(xiàn)Euler場邊界物質(zhì)和能量的流出。灌封膠及炸藥則采用填充方式fill到空氣域中，炸藥瞬態(tài)爆轟后在空氣域中流動。多物質(zhì)Euler場及邊界條件如圖2所示。

圖2 多物質(zhì)填充及Euler邊界條件

2.2 材料模型

空氣采用 IDEALGAS狀態(tài)方程 (equation of state，EOS)描述，方程為：

式中：P0——空氣壓力；

e——單位體積空氣熱力學(xué)能；

γ——多方指數(shù)；

Pshift——壓力偏移量。

對于空氣模型： γ=1.4，ρ取 0.001225 g/cm3，e在mm/mg/ms單位制下取2.068×105。

TNT炸藥采用Jones-Wilkins-Lee狀態(tài)方程，方程及相關(guān)參數(shù)如下：

其中，P為爆轟壓力；V為相對體積；E為單位體積爆轟產(chǎn)物熱力學(xué)能，其余為材料參數(shù)。

對于破片、前后端板、灌封膠和內(nèi)襯，選用沖擊波速度和材料質(zhì)點速度相關(guān)的沖擊波狀態(tài)方程(shock)，即 Mie-Gruneisen 方程；

Mie-Gruneisen方程為：

其中，μ = ρ/ρ0?1，ρ0為材料初始密度；ρ為材料密度；c0為材料聲速。

質(zhì)點速度方程為：

其中，up為沖擊波速度；C0、S1和S2為試驗常數(shù)。

破片和內(nèi)襯強度模型選取Johnson-cook彈塑性方程，該模型一般用于描述材料在大應(yīng)變、高應(yīng)變率下的強度狀態(tài)，非常適合材料在高速沖擊、爆炸載荷環(huán)境中的應(yīng)力計算，屈服強度方程為：

TH——對比溫度；

A——材料低應(yīng)變時的屈服強度；

B和n——表征材料應(yīng)變硬化的常數(shù)；

C和m——試驗常數(shù)。

其中TH=(T?Tr)/(Tm?Tr)，T為材料溫度，Tm為材料熔點，Tr為環(huán)境室溫。

聚氨酯泡沫由無強度材料描述，在沖擊波作用時瞬間破裂。材料的Gruneisen狀態(tài)方程參數(shù)見表2，Johnson-cook強度模型參數(shù)見表3。

表2 Gruneisen狀態(tài)方程參數(shù)

表3 Johnson-cook強度模型材料參數(shù)

2.3 結(jié)果及分析

對建立的數(shù)值模擬模型進行計算，獲得了不同算例的破片初速，從破片初速軸向分布來看，戰(zhàn)斗部兩端速度較低，尤其是前端面，這與Charron、馮順山和蔣建偉等對破片軸向速度分布的研究結(jié)果一致[2,8]。圖3為鋁內(nèi)襯典型工況典型時刻的破片速度分布云圖。

圖3 尺寸5 mm的球形和立方體破片典型時刻的速度云圖

從不同預(yù)制破片在典型時刻的飛散情況可以看出，隨著爆轟波軸向傳播，戰(zhàn)斗部呈喇叭狀膨脹，最后形成前端直徑偏小的鼓狀，破片逐漸分離徑向飛散。邊長5 mm的立方體破片算例中，不同破片的速度-時間曲線如圖4所示。

圖4 邊長5 mm的立方體破片算例中軸向破片速度曲線

可見，戰(zhàn)斗部起爆40 μs后，破片速度已達到最大值，繼續(xù)計算只是破片在空氣中無阻尼飛行，對初速計算結(jié)果已無影響。本文研究的是內(nèi)襯材料，破片類型、尺寸對破片初速的影響規(guī)律，因此需排除端面效應(yīng)、起爆方式等影響因素。根據(jù)破片初速軸向分布特點，取戰(zhàn)斗部軸向60～75 mm范圍內(nèi)破片作重點研究，取其中破片的最高速度為初速數(shù)據(jù)。

爆轟產(chǎn)物驅(qū)動下，不同內(nèi)襯材料的破裂半徑不同。根據(jù)能量守恒理論，破裂半徑越大，爆轟產(chǎn)物泄露的能量越小，破片也就獲得更多動能，從破片初速曲線看內(nèi)襯約在15～25 μs時破裂。以直徑為2 mm的球形破片不同內(nèi)襯材料的工況為例，對20 μs時刻破片上的壓力值進行對比，計算結(jié)果見圖5。

圖5 直徑2 mm球形破片不同內(nèi)襯工況下破片壓力云圖

從壓力峰值分布看，高壓力區(qū)均位于與裝藥接觸一面的破片底部。同質(zhì)量裝藥不同內(nèi)襯材料條件下，同尺寸破片受到的壓力差別較大，峰值差距達94%。可知20 μs時刻三種內(nèi)襯破裂程度各異：聚氨酯泡沫爆炸瞬間破裂，爆轟產(chǎn)物很早泄露，破片受到的壓力最??；鋁與鋼材料相比破裂半徑小，受到的壓力值相對鋼內(nèi)襯低。圖6給出了三種工況破片速度-時間曲線，由動量定理可知，破片受到的壓力值與破片獲得的動量正相關(guān)，壓力值大的工況破片獲得的初速也高。

圖6 兩種典型破片不同內(nèi)襯工況下破片速度曲線

對于2 mm球形，內(nèi)襯45#鋼對應(yīng)的破片速度大于1100鋁對應(yīng)的破片速度，而對于3 mm立方體破片，兩種內(nèi)襯材料下速度接近，分析認(rèn)為由于材料延展性好、膨脹破裂半徑大，內(nèi)襯45#鋼時破片體應(yīng)能獲得更大動能，破片初速更高，但立方體破片與內(nèi)襯貼合性相對較差，較球形破片氣體泄漏更嚴(yán)重，導(dǎo)致多加載到破片體的能量更多被內(nèi)襯材料或泄漏的氣體獲得，因此45#鋼內(nèi)襯時破片初速未見明顯差別。

從圖5中可知破片初速從高到低的順序為：(內(nèi)襯材料) 45#鋼>1100鋁>聚氨酯泡沫；而裝填比大小順序卻是(內(nèi)襯材料)45#鋼<1100鋁<聚氨酯泡沫。若采用經(jīng)典的Gurney公式計算，破片速度基本取決于裝填比，計算結(jié)果與圖6完全相反，與實際情況不符。可見選取不同內(nèi)襯材料進行預(yù)制破片初速研究具有實際應(yīng)用價值。聚氨酯泡沫和45#鋼內(nèi)襯工況破片初速理論及數(shù)值模擬結(jié)果見圖7。

圖7 不同內(nèi)襯材料下破片初速理論及數(shù)值模擬結(jié)果

可見排布和驅(qū)動過程中存在氣體泄漏等問題，導(dǎo)致現(xiàn)有的理論模型計算結(jié)果明顯偏高。

2.4 模型建立

根據(jù)氣體泄露系數(shù)的量綱分析結(jié)果，對數(shù)值模擬得到的破片初速結(jié)果進行處理，便得到數(shù)據(jù)擬合所需的數(shù)據(jù)，見表4和表5。

表4 球形破片氣體泄漏系數(shù)擬合數(shù)據(jù)

表5 立方體破片氣體泄漏系數(shù)擬合數(shù)據(jù)

對氣體泄露系數(shù)fl進行多項式擬合，公式形式如下式所示，a1～a3、b1～b3、d為待定系數(shù)。

對兩種不同形狀的破片進行擬合計算，獲得曲線見圖8，擬合得到的參數(shù)見表6。

圖8 不同形狀破片fl系數(shù)擬合曲線

表6 不同形狀破片初速公式參數(shù)擬合結(jié)果

考慮到數(shù)值模擬中內(nèi)襯材料范圍及破片與裝藥半徑的相對尺寸，給出預(yù)制破片初速模型的定義域為 0.56≤x1≤10，0.05≤x2≤0.32。

3 試驗研究

設(shè)計的預(yù)制破片裝置主要由前后端蓋、裝藥、破片和內(nèi)襯等件組成，預(yù)制破片放置于單端加工有凹槽的聚乙烯圓筒中，結(jié)構(gòu)示意和實物如圖9所示。選取的裝藥為JH-2藥柱，裝藥半徑18 mm，長度50 mm，工況 1、3、5選取Ф5 mm 破片內(nèi)襯分別選取45#鋼、聚乙烯和1100鋁三種材料；工況2、4、6選取Ф6 mm破片內(nèi)襯分別選取45#鋼、聚乙烯和1100鋁三種材料，內(nèi)襯尺寸與數(shù)值計算工況相同。

圖9 試驗裝置

試驗在爆炸洞中進行，預(yù)制破片裝置軸線向上放置于中間支撐桿之上，裝置中心距離地面0.5 m，與測速靶中心平齊。為了獲得準(zhǔn)確的破片初速結(jié)果，分別在距離裝置2 m和4 m處布置了兩組梳狀測速靶。試驗現(xiàn)場布置如圖10所示。

圖10 靜爆試驗現(xiàn)場布置

試驗中取測得數(shù)據(jù)的最高值作為破片速度值。對于不同工況兩公式計算與以上試驗結(jié)果及文獻試驗結(jié)果的對比情況匯總?cè)绫?所示，其中立方體破片數(shù)據(jù)引用文獻[11]的試驗結(jié)果。

表7 試驗結(jié)果與計算結(jié)果對比

從以上結(jié)果及分析可知，經(jīng)典的預(yù)制破片初速公式計算鋼內(nèi)襯工況時誤差相對小，其他內(nèi)襯工況誤差均超過20%；與經(jīng)典的預(yù)制破片初速公式相比，本文提出的預(yù)制破片初速修正模型由于考慮了內(nèi)襯材料、破片形狀和尺寸的影響，與試驗結(jié)果符合度更好，模型大大降低了經(jīng)典公式的計算誤差，精度平均提高了14%以上，可用于預(yù)制破片初速的工程計算。

4 結(jié)束語

1)爆轟產(chǎn)物泄漏對炸藥驅(qū)動預(yù)制破片速度有一定影響，隨著內(nèi)襯材料特性、破片形狀和尺寸等參數(shù)的調(diào)整，預(yù)制破片初速呈現(xiàn)相關(guān)性變化。

2) 通過在經(jīng)典公式基礎(chǔ)上增加爆轟產(chǎn)物泄漏系數(shù)，本文構(gòu)建了一種預(yù)制破片初速的修正計算方法并通過了試驗驗證，結(jié)果顯示相比原公式精度平均提高14%以上。

3)本文結(jié)論對于預(yù)制破片戰(zhàn)斗部設(shè)計和毀傷威力準(zhǔn)確評估具有實際的指導(dǎo)意義。