龐桂鳳，張 琦，李 巍，李 磊

0 引 言

隨著彈道導(dǎo)彈精度的日益提高，整體類戰(zhàn)斗部攻擊硬目標(biāo)的需求日益強烈。對整體侵徹類戰(zhàn)斗部而言，終點效應(yīng)試驗[1,2]是研究戰(zhàn)斗部侵徹強度、侵徹能力、裝藥安定性和引戰(zhàn)配合等性能必不可少的項目。

本文在分析采用地面火箭橇和滑膛火炮完成整體類戰(zhàn)斗部終點效應(yīng)試驗的優(yōu)、缺點的基礎(chǔ)上，結(jié)合本文研究的發(fā)射筒的特點，提出了一種利用導(dǎo)彈發(fā)射筒進(jìn)行大口徑、大質(zhì)量、高速度的整體類戰(zhàn)斗部終點效應(yīng)試驗的方案。仿真結(jié)果表明，方案具有可行性，并可使廢棄的導(dǎo)彈發(fā)射筒得到充分利用。

本文研究的整體戰(zhàn)斗部終點速度馬赫數(shù)大于2，戰(zhàn)斗部試驗件約數(shù)百公斤，具有復(fù)雜的結(jié)構(gòu)約束邊界。采用地面火箭橇進(jìn)行這類戰(zhàn)斗部的終點效應(yīng)試驗，存在因振動環(huán)境惡劣過考核、研制周期長、試驗成本高等缺點；若采用滑膛火炮進(jìn)行戰(zhàn)斗部終點效應(yīng)試驗，則存在炮管加工難度大、生產(chǎn)周期長、制造成本高，因膛內(nèi)壓力過大過考核，結(jié)構(gòu)約束邊界模擬不真實等缺點。

彈道導(dǎo)彈發(fā)射筒口徑大，戰(zhàn)斗部作為導(dǎo)彈的一部分與發(fā)射筒本身具有較好的相容性。初步分析，通過調(diào)整彈射壓力和彈射行程，在可接受的彈射過載條件下，可以獲得理想的彈射出筒速度，達(dá)到戰(zhàn)斗部終點效應(yīng)試驗的要求。而且一般導(dǎo)彈發(fā)射筒規(guī)定使用若干次以后，就廢棄不再使用，因此利用發(fā)射筒進(jìn)行地面試驗是完全可行的。利用廢棄不用的導(dǎo)彈發(fā)射筒進(jìn)行戰(zhàn)斗部終點效應(yīng)試驗，還具有試驗成本低、研制周期短的優(yōu)點。

1 系統(tǒng)及方案

試驗系統(tǒng)由導(dǎo)彈發(fā)射筒、試驗件、高速攝像裝置和終點效應(yīng)試驗靶 4大部分組成，其中導(dǎo)彈發(fā)射筒[3]由發(fā)射筒本體、彈射動力裝置和適配器組成，如圖 1所示。試驗時，導(dǎo)彈發(fā)射筒水平放置，試驗件與適配器固定好后安裝于導(dǎo)彈發(fā)射筒內(nèi)，發(fā)射筒出口處架設(shè)高速攝像裝置，用于錄制試驗件出筒及撞擊試驗靶時刻的姿態(tài)、速度及相關(guān)景像[4]。試驗靶根據(jù)整體戰(zhàn)斗部打擊目標(biāo)的特點進(jìn)行設(shè)計、施工，滿足試驗要求。

圖1 試驗系統(tǒng)組成示意Fig.1 Composition Diagram of Test System

因試驗件在發(fā)射筒內(nèi)的滑行時間很短[5]，為毫秒量級，為便于計算，假設(shè)試驗件在發(fā)射筒內(nèi)所受的摩擦力和空氣阻力恒定不變。根據(jù)能量守恒定律，建立基本方程（1）和輔助方程（2）。

基本方程：

輔助方程：

式中 t為試驗件在發(fā)射筒內(nèi)的運動時間；tv為試驗件t時刻在發(fā)射筒內(nèi)的速度；tF為試驗件t時刻在發(fā)射筒內(nèi)受到的推力；xF為試驗件在發(fā)射筒內(nèi)受到的空氣阻力；ρ為空氣密度，為試驗件t時刻在發(fā)射筒內(nèi)的阻力系數(shù)， )(txvC =0.2；S為試驗件的受力面積；0F為試驗件在筒內(nèi)所受到的摩擦阻力；m為試驗件的質(zhì)量；tL為試驗件t時刻在筒內(nèi)的行程；tΡ為試驗件t時刻在筒內(nèi)受到的壓強；φ為試驗件的底部直徑；μ為滑動摩擦力系數(shù)，μ=0.1。

2 計算方法

本文采用龍格-庫塔法[6]，數(shù)值求解上述方程組，仿真流程如圖2所示。

圖2 仿真流程示意Fig.2 The Simulation Flow Chart

3 仿真結(jié)果分析

3.1 某典型工況仿真結(jié)果分析

試驗件終點速度取700～750 m/s，通過仿真計算可知，所需的彈射壓強分別為8.7 MPa和10 MPa，即若發(fā)射筒內(nèi)的彈射動力裝置產(chǎn)生的壓強處于8.7～10 MPa之間，則試驗件可獲得滿足要求的終點速度。某次典型工況仿真結(jié)果如圖3和圖4所示。

圖3 壓強-時間變化曲線Fig.3 Intensity of Pressure-time Changes Curve

圖4 速度-時間變化曲線Fig.4 Velocity-time Changes Curve

由圖3可知，發(fā)射筒內(nèi)所需的壓力不超過10 MPa，而發(fā)射筒結(jié)構(gòu)可以承受12 MPa的壓強。為分析在試驗過程中試驗件的受載情況，仿真計算了試驗件運動過程中的軸向過載nx，計算結(jié)果如圖5所示。

圖5 軸向過載-時間變化曲線Fig.5 Axial Overload-time Changes Curve

由圖 5可知，在試驗過程中試驗件要承受的最大軸向過載為3500，強度計算結(jié)果表明，試驗件強度滿足試驗要求。因此，該方案可行。

3.2 其它工況仿真結(jié)果分析

為進(jìn)一步分析這種終點試驗方案的可行性及適用范圍，仿真計算了另外2種工況：

a）在發(fā)射筒長度（9 m）和試驗件質(zhì)量（500 kg）一定的條件下，試驗件終點速度與發(fā)射筒內(nèi)壓強的關(guān)系，仿真結(jié)果如圖6所示。

圖6 終點速度-時間變化曲線Fig.6 Terminal Speed-time Changes Curve

由圖 6可知，在發(fā)射筒長度和試驗件質(zhì)量一定的條件下，發(fā)射筒內(nèi)壓強越大，試驗件達(dá)到的終點速度越大，當(dāng)發(fā)射筒壓強在7.4～8.3 MPa之間時，試驗件的終點速度范圍是660～700 m/s。

b）在發(fā)射筒長度（9 m）和發(fā)射筒內(nèi)壓強（8.3 MPa）一定的條件下，試驗件終點速度與試驗件質(zhì)量的關(guān)系，仿真結(jié)果如圖7所示。

圖7 終點速度-時間變化曲線Fig.7 The Terminal Speed-time Changes Curve

由圖7可知，在發(fā)射筒長度和發(fā)射筒內(nèi)壓強一定的條件下，試驗件質(zhì)量越小，試驗件可以達(dá)到的終點速度越大，當(dāng)試驗件質(zhì)量在350～500 kg之間時，試驗件的終點速度范圍是700～839 m/s。

4 結(jié)束語

分析認(rèn)為，若根據(jù)內(nèi)彈道計算結(jié)果重新設(shè)計發(fā)射筒[7,8]的彈射動力裝置，進(jìn)行500 kg以內(nèi)、終點速度小于750 m/s的戰(zhàn)斗部終點效應(yīng)試驗，方案具有可行性。

仿真結(jié)果表明，試驗件質(zhì)量增加或發(fā)射筒壓強減少，都將引起試驗件終點速度降低，反之，則試驗件的終點速度會升高。故本試驗方案是否適用其它質(zhì)量、結(jié)構(gòu)的戰(zhàn)斗部終點效應(yīng)試驗，需從終點速度及結(jié)構(gòu)強度等方面進(jìn)行分析，以確定其適用性。

[1] 隋樹元, 王樹山. 終點效應(yīng)學(xué)[M]. 長沙： 國防工業(yè)出版社, 2000.

Sui Shuyuan, Wang Shushan. The effect of terminal point[M]. Changsha：Defence Industry Book Concern, 2000.

[2] 王儒策, 趙國志. 彈丸終點效應(yīng)學(xué)[M]. 北京： 北京理工大學(xué)出版社,1993.

Wang Ruce, Zhao Guozhi. Projectile end point effect[M]. Beijing： Beijing Institute of Technology Press, 1993.

[3] 高明坤, 宋廷倫. 火箭導(dǎo)彈發(fā)射裝置構(gòu)造[M]. 北京： 北京理工大學(xué)出版社, 1996.

Gao Mingkun, Song Tinglun. The construction of rocket’s or missile’s launcher[M]. Beijing： Beijing Institute of Technology Press, 1996.

[4] 張國偉. 終點效應(yīng)及其應(yīng)用技術(shù)[M]. 北京： 國防工業(yè)出版社, 2006.

Zhang Guowei. End point effect and application technique[M]. Beijing：Defence Industry Book Concern, 2006.

[5] 姚昌仁. 唐國良, 宋廷倫. 火箭導(dǎo)彈發(fā)射動力學(xué)[M]. 北京： 北京理工大學(xué)出版社, 1996.

Yao Changren, Tang Guoliang, Song Tinglun. The launch dynamics of rocket or missile[M]. Beijing： Beijing Institute of Technology Press, 1996.

[6] 封建湖, 車剛明, 聶玉峰. 數(shù)值分析原理[M]. 北京： 科學(xué)出版社, 2002.

Feng Jianhu, Che Gangming, Nie Yufeng. The principle of numerical analysis[M]. Beijing： Science Book Concern, 2002.

[7] 姚昌仁. 導(dǎo)彈發(fā)射裝置設(shè)計[M]. 北京： 國防工業(yè)出版社, 1981.

Yao Changren. The design of missile’s launcher[M]. Beijing： Defence Industry Book Concern, 1981.

[8] 姚昌仁, 張波. 火箭導(dǎo)彈發(fā)射裝置設(shè)計[M]. 北京： 北京理工大學(xué)出版社, 1998.

Yao Changren, Zhang Bo. The design of rocket’s or missile’s launcher[M].Beijing： Beijing Institute of Technology Press,1998.