（海軍航空工程學(xué)院 煙臺(tái) 264001）

1 引言

由于不同目標(biāo)的要害分布以及性質(zhì)差別很大，所以在計(jì)算戰(zhàn)斗部條件殺傷概率時(shí)，要對(duì)目標(biāo)的要害分布類(lèi)型做出假設(shè)。通常的假設(shè)方法有兩種：一是要害艙段殺傷模型，這是針對(duì)給定具體尺寸、位置、性質(zhì)的目標(biāo)時(shí)使用的模型；二是直桿型殺傷模型，它將目標(biāo)假設(shè)為由多條等效直桿組成，便于在不對(duì)特定的目標(biāo)結(jié)構(gòu)進(jìn)行詳細(xì)分析時(shí)使用。本文探究的是利用概率密度積分法，提出了在針對(duì)直桿型殺傷模型時(shí)，防空導(dǎo)彈對(duì)其單發(fā)殺傷概率的計(jì)算方法。

2 單發(fā)殺傷概率一般表達(dá)式

[1]，防空導(dǎo)彈的單發(fā)殺傷概率在由脫靶量ρ，脫靶方位角θ以及xr軸構(gòu)成的圓柱坐標(biāo)系中可以寫(xiě)為

式中Pf（ρ，θ）是脫靶條件為ρ，θ下，引信的啟動(dòng)概率；fg（ρ，θ）是制導(dǎo)誤差為ρ，θ的二維分布密度函數(shù)；ff（xr｜ρ，θ）是引信啟動(dòng)點(diǎn)一維坐標(biāo)xr分布密度函數(shù)。

定義引信戰(zhàn)斗部聯(lián)合條件殺傷概率Pdf（xr｜ρ，θ）為

綜上所述可得：

3 單發(fā)殺傷概率三維數(shù)字積分計(jì)算

在進(jìn)行三位數(shù)字積分之前，可先用簡(jiǎn)化的解析法估算單發(fā)殺傷概率值。進(jìn)行簡(jiǎn)化解析法計(jì)算時(shí)，假設(shè)目標(biāo)為單一直桿型模型，在制導(dǎo)脫靶的可能范圍內(nèi)啟動(dòng)概率Pf＝1，同時(shí)只考慮戰(zhàn)斗部條件殺傷概率，系統(tǒng)誤差為0。通過(guò)參考文獻(xiàn)[1]中的公式，可以估算出：若戰(zhàn)斗部特征半徑與制導(dǎo)誤差的比值不小于4，則可以得到不低于98%的單發(fā)殺傷概率。

而在進(jìn)行三維數(shù)字積分時(shí)，由于彈目交會(huì)情況復(fù)雜，首先要將實(shí)際脫靶量ρ分成兩部分：0～ρ0以及ρ0～ρmax，在0＜ρ＜ρ0時(shí)，戰(zhàn)斗部的威力足以摧毀目標(biāo)，其條件殺傷概率為1，因此，可以得到圓柱坐標(biāo)下的計(jì)算單發(fā)殺傷概率的公式為

其中制導(dǎo)誤差的概率密度函數(shù)fg（ρ，θ）的計(jì)算式為

式中的my、mz為脫靶參數(shù)中的系統(tǒng)誤差，σy、σz為隨機(jī)誤差標(biāo)準(zhǔn)差，兩個(gè)方向上的誤差是相互獨(dú)立的[1]。

3.1 在給定脫靶圓內(nèi)的殺傷概率

在脫靶誤差服從二維正態(tài)分布的條件下，落入給定半徑ρ0內(nèi)的概率可以通過(guò)二重積分或蒙特卡洛法對(duì)上式求得結(jié)果。但是在實(shí)際應(yīng)用中，計(jì)算落入定脫靶圓內(nèi)的概率的情況可以分為圓散布和橢圓散布兩種狀態(tài)，橢圓散布代表著實(shí)際彈道的標(biāo)準(zhǔn)差σy≠σz。簡(jiǎn)便起見(jiàn)，本文假設(shè)彈道為圓散布，即σy＝σz＝σ。在這種情況下，設(shè)t＝ρ／σ，a＝m／σ，最終落入給定脫靶圓的概率可寫(xiě)成：

在實(shí)際彈道的散布中心與質(zhì)心相重合的時(shí)候，系統(tǒng)誤差等于零。在大多數(shù)情況下，系統(tǒng)誤差可以通過(guò)加入校正信號(hào)予以消除，對(duì)于一個(gè)成熟的導(dǎo)彈武器系統(tǒng)而言，都可以認(rèn)為其系統(tǒng)誤差為零，即m＝0，則a＝0，因此零階虛位移貝塞爾函數(shù)I0（aσ／ρ0）＝1，又因?yàn)镮n（aσ／ρ0）＝0，n＞0，則得到最終的落入脫靶圓的概率為[2]

3.2 給定脫靶圓到最大可能脫靶量之間的殺傷概率

3.2.1 直桿型殺傷模型的洞穿概率

在脫靶量ρ大于給定的值ρ0時(shí)，就需要計(jì)算戰(zhàn)斗部的條件殺傷概率。計(jì)算時(shí)要考慮破片對(duì)目標(biāo)要害部位的殺傷作用以及沖擊波對(duì)目標(biāo)的爆破毀傷作用。由于目標(biāo)要害分布的不同則計(jì)算戰(zhàn)斗部條件殺傷概率的方法就不同。本文假設(shè)目標(biāo)為直桿殺傷模型，即目標(biāo)由多條等效直桿組成，每條直桿代表著機(jī)體上的重要部件諸如機(jī)翼、發(fā)動(dòng)機(jī)等。此模型假設(shè)當(dāng)任一根直桿落入戰(zhàn)斗部破片的動(dòng)態(tài)飛散帶時(shí)就認(rèn)為目標(biāo)受到一定概率的殺傷，而不單獨(dú)考慮命中破片數(shù)。在這種假設(shè)下：

其中，m為目標(biāo)可殺傷直桿的總數(shù)，Pdj為戰(zhàn)斗部對(duì)第j個(gè)直桿的殺傷概率。這里不考慮引燃概率，引爆概率與沖擊波毀傷概率，只考慮洞穿殺傷概率。Pdj由Pd1和Pd2兩部分組成，Pd1是駕駛艙等重要單個(gè)設(shè)備艙段的洞穿殺傷概率，Pd2為幾個(gè)艙段組合殺傷概率，即只有在殺傷全部艙段組合才能殺傷目標(biāo)。Pd1表達(dá)式為

其中，j1是這類(lèi)艙段的總數(shù)，Nj是j艙段命中破片的平均數(shù)，P（Ej）是單枚破片對(duì)j艙段的洞穿殺傷概率。

Ej為破片平均比能，即單位破片平均迎風(fēng)面積上的動(dòng)能。

3.2.2 要害艙段命中破片數(shù)的計(jì)算

圖1所示為戰(zhàn)斗部動(dòng)態(tài)飛散過(guò)程中各個(gè)矢量的意義以及之間的關(guān)系。

圖1 戰(zhàn)斗部破片靜態(tài)、動(dòng)態(tài)飛散角及飛散速度

通過(guò)圖1中所標(biāo)注的參數(shù)可以求得j艙段命中破片的平均數(shù)Nj：

式中已知Ωj、ωj以及第j個(gè)要害艙段中心與戰(zhàn)斗部中心的連線矢量Rj，于是可以推導(dǎo)出φj和破片飛散密度κ（φj）。Sj是第j個(gè)要害艙段在破片動(dòng)態(tài)飛散方向上的等效投影面積。在計(jì)算Sj時(shí)，根據(jù)目標(biāo)要害不同位置的不同形狀，做出如下三種不同的等效幾何模型：

1）駕駛艙或油箱等可以等效為球體的要害艙段，其投影面積為常數(shù)，即Sj＝const；

2）發(fā)動(dòng)機(jī)艙、機(jī)身等要害位置，其投影面積為端面與側(cè)面投影面積之和，即Sj＝Sj0｜c(diǎn)osqj｜＋Sj1｜sinqj｜；

3）機(jī)翼、尾翼等位置，其投影面積為這些要害位置在Rj上的投影，即Sj＝Sj0｜c(diǎn)osqj｜。

3.2.3 引信啟動(dòng)區(qū)

在引信啟動(dòng)區(qū)間內(nèi)，引信戰(zhàn)斗部聯(lián)合條件殺傷概率Pdf（xr｜ρ，θ）變?yōu)?/p>

Pf（xr｜ρ，θ）是引信的啟動(dòng)概率，在引信啟動(dòng)區(qū)范圍內(nèi)，計(jì)算式為

在導(dǎo)彈接近目標(biāo)時(shí)，通常近似認(rèn)為引信啟動(dòng)點(diǎn)沿相對(duì)速度坐標(biāo)系oxm軸的分布服從一維正態(tài)分布規(guī)律，其分布密度函數(shù)為

式（14）中的兩個(gè)未知量，mx是啟動(dòng)點(diǎn)的數(shù)學(xué)期望，對(duì)于光學(xué)引信來(lái)講，光學(xué)引信視場(chǎng)較窄，mx與引信延遲時(shí)間τ以及導(dǎo)彈戰(zhàn)斗部中心相對(duì)速度坐標(biāo)系內(nèi)xr軸上的坐標(biāo)有關(guān)，計(jì)算式為

引信啟動(dòng)點(diǎn)的散布σx由四部分的平方和構(gòu)成：延遲時(shí)間散布造成的啟動(dòng)點(diǎn)散布σ1，光學(xué)引信主光傾角散布所造成的啟動(dòng)點(diǎn)散布σ2，引信靈敏度變化引起的啟動(dòng)點(diǎn)散布σ3以及目標(biāo)反射信號(hào)起伏所造成的啟動(dòng)點(diǎn)的散布

由以上的分析，可以依次求得在直桿型殺傷模型的情況下的洞穿概率，引信戰(zhàn)斗部聯(lián)合殺傷概率，進(jìn)而可以求出給定脫靶圓到最大可能脫靶量的殺傷概率，最終得到防空導(dǎo)彈單發(fā)殺傷概率的結(jié)果。

4 仿真計(jì)算

在進(jìn)行仿真計(jì)算時(shí)，以某型防空導(dǎo)彈為例，整個(gè)計(jì)算過(guò)程的流程如圖2所示。

圖2 計(jì)算導(dǎo)彈單發(fā)殺傷概率的流程

由以上流程，防空導(dǎo)彈使用激光引信，假設(shè)該防空導(dǎo)彈破片戰(zhàn)斗部的參數(shù)為：彈目交會(huì)遭遇點(diǎn)高度為22000m，遭遇點(diǎn)斜距為30000m，導(dǎo)彈速度為1000m／s，彈道傾角18°，目標(biāo)運(yùn)動(dòng)速度600m／s，導(dǎo)彈制導(dǎo)系統(tǒng)誤差m＝0，隨機(jī)誤差的標(biāo)準(zhǔn)差σ＝4m，戰(zhàn)斗部破片數(shù)為5000塊，破片的初速度為3000m／s，平均每塊破片的質(zhì)量為6g，破片飛散方位角為30°，飛散方向角為60°。

由之前敘述的計(jì)算方法，假設(shè)模型為直桿型模型，其中1艙段等效為圓柱形，2、3、4、5艙段等效為平板型，計(jì)算出每個(gè)艙段命中的破片數(shù)，其中1艙段命中1000枚，2、3艙段各1500枚，4、5艙段各500枚。目標(biāo)的模型如圖3所示。

圖3 直桿型殺傷模型的目標(biāo)易損構(gòu)件

由以上參數(shù)可以得出當(dāng)系統(tǒng)誤差為0時(shí)單發(fā)防空導(dǎo)彈的落入概率隨脫靶量ρ0變化的關(guān)系，其關(guān)系如圖4所示。

圖4 落入概率隨脫靶量ρ0的變化曲線

通過(guò)前文的論述，當(dāng)脫靶量ρ0為4m時(shí)，由公式：ρmax＝m＋kσ可以得到最大可能脫靶量ρmax為16m。經(jīng)過(guò)計(jì)算，在給定脫靶圓內(nèi)的殺傷概率為39.34%，在給定脫靶圓之外到最大可能脫靶量之間，引信戰(zhàn)斗部聯(lián)合殺傷概率為97.03%，所以殺傷概率為58.80%，兩者求和，最終確定該型防空導(dǎo)彈的單發(fā)殺傷概率為98.14%，驗(yàn)證了前文當(dāng)戰(zhàn)斗部特征半徑與制導(dǎo)誤差的比值大于等于4，則單發(fā)殺傷概率不低于98%這一估算結(jié)果。

5 結(jié)語(yǔ)

防空導(dǎo)彈從探測(cè)目標(biāo)到摧毀目標(biāo)的實(shí)際情況十分復(fù)雜，實(shí)際應(yīng)用中還需要考慮環(huán)境因素，干擾因素等。本文僅選取了落入脫靶圓內(nèi)的分布為圓分布這一特殊情況進(jìn)行了計(jì)算，并且假設(shè)武器系統(tǒng)技術(shù)成熟，系統(tǒng)誤差等于零，在實(shí)際中由于種種因素會(huì)導(dǎo)致單發(fā)殺傷概率降低。另外，針對(duì)不同類(lèi)型的近炸引信，如無(wú)線電引信等，其引信啟動(dòng)區(qū)的各種參數(shù)的不同會(huì)導(dǎo)致單發(fā)殺傷概率的不同。在影響單發(fā)殺傷概率的諸多因素，例如引戰(zhàn)配合效率，戰(zhàn)斗部威力，目標(biāo)易損性等因素。本文針對(duì)這些關(guān)鍵因素進(jìn)行了計(jì)算，得到了在該目標(biāo)類(lèi)型下，針對(duì)激光引信的單發(fā)殺傷概率計(jì)算方法，在一定程度上可以為武器系統(tǒng)的效能評(píng)估提供參考。

參考文獻(xiàn)

[1]張志鴻.防空導(dǎo)彈引信與戰(zhàn)斗部配合效率和戰(zhàn)斗部設(shè)計(jì)[M].北京：宇航出版社，1994：34－88.

[2]李廷杰.導(dǎo)彈武器系統(tǒng)的效能及其分析[M].北京：國(guó)防工業(yè)出版社，2000：166－183.

[3]鄭平泰，楊濤，王寶敏.攔截彈單發(fā)殺傷概率仿真計(jì)算研究[J].系統(tǒng)仿真學(xué)報(bào)，2006，18（7）：1763－1765.

[4]馮融冰，鄭世勇，程延杰，等.基于Simulink的防空導(dǎo)彈單發(fā)殺傷概率仿真研究[J].艦船電子工程，2011，31（11）：82－84.

[5]王征，李言俊，雷蕾.破片型空空彈殺傷概率三維建模與仿真[J].系統(tǒng)仿真學(xué)報(bào)，2009，21（17）：5340－5344.

[6]王淼，宋晗，張麗群.防空導(dǎo)彈單發(fā)殺傷概率仿真計(jì)算研究[J].戰(zhàn)術(shù)導(dǎo)彈技術(shù)，2009（4）：93－96.

[7]孫剛，廖瑛.引戰(zhàn)系統(tǒng)數(shù)學(xué)仿真中戰(zhàn)斗部仿真的初步實(shí)現(xiàn)[J].上海航天，2006（5）：50－54.

[8]王樹(shù)林，段江紅，楊爭(zhēng)，等.用破片覆蓋度進(jìn)行毀傷概率評(píng)估[J].現(xiàn)代防御技術(shù)，2008，36（1）：36－40.

[9]高波，邱立軍，姚躍廷.防空導(dǎo)彈單發(fā)殺傷概率研究[J].彈道學(xué)報(bào)，2011，23（4）：52－55.

[10]王潤(rùn)平，王正向.地空導(dǎo)彈單發(fā)殺傷概率之解析算法[J].系統(tǒng)工程與電子技術(shù)，2001，23（5）：39－40.