堅硬目標毀傷效能評估的模擬模型

摘 要：在堅硬目標毀傷效能評估的理論模型的基礎上，建立計算機模擬模型，依據防護結構毀傷概率的數學模型，采用蒙特卡洛模擬方法模擬其毀傷概率。

關鍵詞：堅硬目標 毀傷效能 毀傷模型

中圖分類號：TP761 文獻標識碼：A 文章編號：1672-3791（2013）04（a）-0036-02

在建立堅硬目標特性數據庫[1]、堅硬目標防護層等效素混凝土靶[2]的基礎上，討論了毀傷評估的理論方法[3]，依據防護工事的特點，以及戰斗部的毀傷方式，初步定義了堅硬目標防護結構的毀傷等級，并在此基礎上確定了不同等級的毀傷判據，下面據此建立計算機模擬模型，采用蒙特卡洛模擬方法模擬其毀傷概率。

1 計算機模擬模型

根據系統誤差的正態分布規律，在區域內進行隨機抽樣，模擬初系統誤差，采用數值積分法計算出，于是的數學期望為：

當隨機抽樣樣本容量足夠大時，導彈落入有效幅員概率的估計值為：

同理可知，導引誤差為正態分布隨機變量z分布密度函數，因此的數學期望：

當隨機抽樣樣本容量足夠大時，導彈在落入有效幅員條件下的毀傷概率的估計值為：

求出單發毀傷概率，進而求出在不考慮毀傷積累情況下的多發毀傷概率。即：

精度檢驗及樣本容量的確定：的估計誤差與樣本容量之間的關系為：，其中：為置信度相對應的分位數，為隨機變量的方差，計算中以次預試驗得到的樣本方差作為的估計值代入計算。同理可得，的估計誤差與樣本容量的關系。當估計誤差，具有置信度的樣本數為：。

蒙特卡洛法：在計算機模型中，產生隨機誤差的方法主要是運用蒙特卡洛法。本文中采用數學方法計算隨機數列。由于計算需要，需要分別了解均勻分布、正態分布、和互相獨立的二維正態分布的隨機變量抽樣值。（1）均勻分布的隨機抽樣：設是在[0，1]區間上的均勻分部隨機量的一次抽樣值，則×即為在[a，b]區間上均勻分布隨機量的一次抽樣值；（2）正態分布的隨機抽樣：采用極限近似法，設是[0，1]區間上均勻分布隨機量的12次抽樣值，則由產生一個標準正態分布隨機數，如果需要正態隨機數滿足，則×；（3）互相獨立的二維正態分布隨機抽樣：設落點坐標x與y分別服從一維正態分布：，而且相互獨立，則由下式對x，y抽樣：，，其中，為按瑞利分布的隨機變量的抽樣值，為均勻分布的隨機相位的抽樣值：，其中：為（0，1）區間上均勻分布隨機數的兩次抽樣值。

由試驗數據衡量樣本的均方差：本文利用殘余誤差得出樣本均方差的估計值。設對某量進行n次測量，得到n個測量值，被測量值的真值為，隨機誤差的抽樣值為：

其中：為測量值的樣本均值，且為殘余誤差，。于是，求得樣本方差：，因為：

根據隨機誤差的低償性，當n充分大時：，這樣可以將（4.46）式改寫為：，于是得：，所以：，此式稱為貝賽爾（Bessel）公式，根據此式可以由殘余誤差得到單次測量的樣本均方差。

2 防護結構震塌厚度計算

關于侵徹深度的計算，詳見文章[2]。下面討論如何計算防護結構的震塌厚度問題。在爆炸或沖擊載荷作用下，防護結構將產生彈坑、貫穿和震塌三種形態的局部破壞。由于應力波在復雜的鋼筋混凝土結構中難以確定的傳播過程，使得震塌破壞問題成為相當復雜的理論問題，這里采用漏斗坑半徑表示震塌厚度，公式為：，式中：w為射彈裝藥重量，K為標準拋擲爆破系數，e為炸藥換算系數，為孔口堵塞系數，W為炸藥中心至地表距離。

3 實例計算

戰斗部參數如下表1，靶板尺寸為20×18×6 m，彈體垂直侵徹，著速為400 m/s。導引誤差均方差為0.1 m，起爆點誤差帶0.2 m，圓概率偏差為6。

由空腔膨脹理論[2]計算最大侵深，并用最小漏斗坑半徑計算最小震塌厚度。根據兩者和與靶板厚度的對比，確定靶板是否能貫穿。可以計算得最大侵深=3.823668，最小震塌厚度=3.20478，因為且，所以毀傷等級為B，毀傷概率：。本例中，由于等效靶的尺寸相對來講不是很大，所以從計算結果可以看出，散布圓概率偏差對目標毀傷概率的影響較大，其變化規律為散布圓概率偏差越大，單發毀傷概率越小。當等效靶的范圍變大為100×100×6m，由于靶的易損性范圍變大，散布圓概率偏差在0～30之間基本上沒有變化。

4 結論

為量化毀傷效能評估，分別建立了由射擊誤差決定的命中概率和命中條件下的毀傷概率計算方法，進而求出了受射擊誤差影響的毀傷概率[3]。根據上述理論基礎建立計算機模擬模型，為下一步的毀傷軟件結構的搭建奠定基礎。

參考文獻

[1]徐輝，李加財.堅硬目標特性數據庫的建立[J].今日科苑，2007（4）.

[2]徐輝，李加財.堅硬目標防護層等效素混凝土靶的方法[J].科技風，2008.

[3]徐輝，潘陽.堅硬目標毀傷效能評估的理論模型[J].科技風，2012.