基于改進毀傷矩陣的人員目標毀傷效能評估

田浩成，盧芳云，李志斌，王 碩

(國防科技大學 文理學院， 長沙 410005)

1 引言

現代戰場已不再遠離城鎮，人員附帶毀傷問題日益突出，準確評估武器對人員的毀傷效能對火力籌劃具有實際意義。武器對目標造成的毀傷概率分布稱為毀傷矩陣；餅切函數和Carlton函數是毀傷概率分布的數學表達，可視為毀傷矩陣的簡化表示。武器毀傷效能一般基于量化指標進行評定，如常見的平均有效毀傷面積(MAE)。精確的毀傷矩陣是武器毀傷效能評估的重要支撐，需基于一系列相對復雜的計算流程來獲取。

對于人員目標而言，常規武器戰斗部主要依靠裝藥爆炸產生的沖擊波和高速破片造成殺傷。由于沖擊波在自由場衰減過快，僅對較近處產生威脅，因此，武器對人員的有效殺傷范圍由破片場決定。破片戰斗部是一種常見的主要依靠破片毀傷目標的武器，其破片場參數通常由靶場試驗獲得，基于理論推導或經驗公式的破片場預估方法也得到了廣泛應用。近年，基于各類商業軟件的數值模擬已經成為破片飛散特性研究的主要方法，各種戰斗部參數化建模平臺的出現也使得破片場模型具有精度可靠、便于獲取的優點。

在人員的毀傷效能評估研究中，之前常用長方體來替代人體或者使用Poisson分布函數直接模擬命中概率，進而結合破片場的分布密度獲得破片命中人體的信息。但是，人員與武器的相對位置不同時，破片命中人體的部位會有很大差異。當今高精度人體模型已經在汽車碰撞、人員跌落和防護等多個研究領域得到了廣泛應用，精細的組織模型支撐了更準確的損傷評估結果。在計算機支持下，采用精細模型可為獲取彈-目交會條件提供更詳細具體的信息，使得進一步提高毀傷矩陣的準確性和效能評估的可靠性成為可能。

彈-目交會涉及到破片場中多枚破片命中目標的情況，需要綜合每一枚破片的殺傷效果。對于破片這類高速投射物，目前廣泛使用的是動能殺傷標準，如對于人員目標，規定質量1 g以上且動能不低于98 J的破片為有效殺傷破片。但是同一破片命中人體位置不同，殺傷效果也會有很大差異，因此該準則只能用于較為粗略的評估，更為準確的評估結果還需結合人員易損性。Allen和Sperrazza基于人體解剖截面進行了人體易損性分析，并擬合出A-S準則，得到了人員目標在進攻、防守、預備和補給4種作戰任務下失能概率與破片打擊條件之間的關系，已經應用于美軍的人員目標毀傷效能評估程序以及國內的一些毀傷評估方法中。

為了更準確獲得單發破片戰斗部對人員目標的毀傷矩陣，本文采用精細的破片場仿真模型和人體模型，改進了毀傷矩陣計算流程中彈-目交會信息的獲取方式，由破片命中人體的部位、數目、質量和速度，基于部件級的A-S準則計算得到了改進的毀傷矩陣。通過引入基于毀傷矩陣的破片平均有效毀傷面積(MAE)來表征毀傷效能，分析了終點姿態對毀傷效能的影響規律，并使用Monte Carlo方法模擬了武器隨機命中情況，進一步評估了單發破片戰斗部對人員的平均毀傷效能。

2 失能概率

2.1 單破片命中

評定單枚破片對人員的損傷是一個復雜的過程，快速計算人員的失能概率需依賴殺傷準則。考慮到人體軀段的易損性不同，一枚破片命中人體不同部位會造成人員能力下降的程度不同，A-S準則在易損性分析結果基礎上以失能概率的形式表征了人員失能情況，可以較好地表示出從無損傷到完全失能的中間階段。因此，本文使用A-S準則來計算單枚破片對人體造成的損傷，其具體形式為：

=1-exp[-(32-)]

(1)

式(1)中：為破片質量；為破片打擊速度；=32也稱為Sperrazza能量；、、為不同殺傷定義下的擬合系數。實際上是造成失能的閾值，質量單位為格令，速度單位為英尺每秒。

由破片的質量和打擊速度可得到該枚破片的Sperrazza能量，根據破片類型和殺傷定義查找相關系數，即可利用式(1)計算出人員失能概率，表1給出了防守30 s和補給12 h殺傷定義下球形鋼質破片命中裸露人體的部件級A-S準則系數。實際上，根據各部位呈現面積對部件級系數加權平均可以得到A-S準則的全身系數，描述了單枚破片命中人體造成的平均能力下降程度，采用該套全身系數得到的結果對于同一破片命中任意部位而言是相同的。因此，通過破片命中位置確定部件級A-S準則系數，有助于提高人員失能概率的準確性。在美軍的作戰任務定義中，人員的能力與四肢狀態緊密關聯，不同作戰任務對四肢的要求也有很大差異，在被相同條件的破片擊中情況下，承擔補給任務的人員會比承擔防守任務能力下降程度更高，因此，補給12 h殺傷定義常被用來評估對平民造成的附帶毀傷。

2.2 多破片命中

單發破片戰斗部的破片場可以包含多達數千枚的破片，人員目標可能面臨不止一枚破片的威脅，還需綜合考慮多枚破片聯合作用的毀傷效果。假設每一枚破片造成的人員失能是相互獨立的，那么根據概率計算方法可以通過生存法則計算多破片命中人體造成的綜合失能概率，具體數學公式可寫為：

(2)

式(2)中：為第枚命中人體破片造成的失能概率，可通過查詢相應系數由式(1)計算得到；為綜合失能概率。

表1 A-S準則系數

3 改進的毀傷矩陣計算

3.1 爆點

武器的最終爆炸位置通常被視為一個點，稱之為爆點，此點距地面的高度即為爆高。爆點在水平面內的坐標由瞄準點和武器精度共同決定，而爆高主要受彈藥引信的影響。

武器精度是武器本身的一種特性，通常使用圓概率偏差CEP來表征水平面內爆點位置與瞄準點之間的偏差，具體定義為：50%彈著點落入以預期平均彈著點為圓心的圓的半徑，即彈著點落入半徑為CEP的圓內的概率為50%。假設這個偏差在和這2個維度上是相互獨立的，并且服從具有同一標準差的正態分布，那么根據概率論知識，與CEP之間的關系可以表示為：

=1177 4

(3)

通過式(3)，可以由武器CEP計算在和上的標準差，進一步隨機生成爆點坐標，從而實現爆點位置的隨機抽樣模擬，圖1展示了對武器CEP為10 m的爆點水平位置進行 1 000次隨機抽樣結果，藍色區域為可能的人員分布區域，綠色圈內抽樣次數為503次，滿足CEP的定義。

圖1 爆點水平位置模擬圖

為了限制人員進入管控區域，需要保證單發武器能夠實現對一定范圍區域的火力壓制，通常會采用定高引信達到較大的威脅范圍。因此，可以假設武器爆高為定值，如5 m、10 m 和15 m，結合上面對武器精度的模擬可以進一步模擬武器爆點。

3.2 彈-目交會

武器在爆點處還會具有一定的姿態，此時以爆點在地面上的投影為原點，建立地面坐標系，武器與地面人員目標的幾何關系如圖2所示，其中為爆高、為落角、為偏航角、為滾轉角。

圖2 武器與人員目標相對位置示意圖

假設不考慮地形起伏，那么所有設定人員目標均處在同一平面，即人員目標的位置由其在坐標系上的坐標確定。實際上，武器爆高決定了武器在圖2地面坐標系中的坐標，而落角、偏航角、滾轉角確定了戰斗部在爆點時的姿態。此時，通過武器和人員目標的一些簡單信息，就可以確定其位置關系。

將特定的破片場仿真模型和精細人體模型放置到指定坐標，按照設定姿態調整模型，并通過附加武器終點落速的方式模擬動爆條件下的破片場，進而追蹤得到每一枚破片與人體的交會情況，最終統計得到命中人體不同部位的破片的初始速度、質量和打擊位置。但是破片速度會在空氣中逐漸衰減，因此還需對打擊速度進行修正，考慮到破片在空中運動的速度衰減規律為：

(4)

式(3)中：為破片質量；為空氣密度；為破片迎風面積；為空氣阻力系數；為破片飛行速度。

對于球形破片，迎風面積取為=π4，取為097，通過式(4)及破片由初始位置到人員目標的距離，即可得到破片擊中目標時的速度。

3.3 計算流程

按照一定的分布規律改變人體模型與武器爆點之間的相對位置，每個位置處均對應了一個失能概率，將所有失能概率存儲為矩陣形式，即得到了單彈對人員目標的毀傷矩陣。本文對毀傷矩陣的改進計算流程主要體現在：① 應用精細模型獲取詳細的彈-目交會信息；② 確定破片命中人體的具體部位，利用部件級A-S準則系數計算人員的失能概率。

因此，可導入特定的戰斗部破片場仿真模型和人體模型得到特定場景下的彈-目交會情況，并基于如表1所示的部件級系數表計算相應的單發破片戰斗部對特定人員目標的毀傷矩陣改進結果。通常情況下會考慮一個較大區域作為目標區域，將區域進行網格劃分并記錄節點處的失能概率，通過編程可以實現改進的毀傷矩陣的便捷計算，具體計算流程如圖3所示。

圖3 單發破片戰斗部的改進毀傷矩陣計算流程框圖

4 實例分析

4.1 模型信息

..破片場仿真模型

本文考慮基于LS-DYNA數值計算得到的203 mm口徑破片戰斗部破片場，破片場由2 072枚預制鋼質球形破片組成，每枚破片質量為1.56 g，破片速度在1 089.6～1 908.5 m/s，平均速度為1 654.3 m/s，具體分布見圖4。

圖4 破片場示意圖

通過對模型平移、旋轉和速度疊加得到多種終點條件下的動態破片場分布。

..精細人體模型

本文選取一個身高為184.5 cm的精細男性人體模型，并按照A-S準則的身體結構將其劃分為頭頸部、胸部、腹部、盆骨、手和腿等6個部位，具體模型如圖5所示。

圖5 人體模型示意圖

4.2 失能概率計算

假設武器爆高10 m，落速500 m/s，落角75°，滾轉角90°，無偏航。將人體模型放置于以武器爆點為中心的地面坐標系上，坐標設為(-8 m，10 m)和(0 m，-15 m)，2個位置處的具體場景分別如圖6(a)和圖6(b)所示，相應的命中人體的破片統計信息及失能概率結果見表2。

由表2可知，圖6(a)中人員目標僅被單枚破片命中胸部，其綜合失能概率即為該枚破片造成的失能概率，需依據表1中胸部對應的A-S準則系數，可得到目標在防守30 s和補給12 h殺傷定義下的失能概率；圖6(b)中共有6枚破片命中人體，同樣地，按照破片命中部位查詢相應的表1系數，得到2種殺傷定義下每枚破片造成的失能概率，進而利用式(2)得到了目標遭到多破片命中的綜合失能概率。

圖6 彈-目交會場景示意圖

表2 命中破片統計信息及失能概率結果

4.3 改進的毀傷矩陣計算

..單發破片戰斗部打擊下的毀傷矩陣計算

假設目標區域為[-150 m，150 m]×[-150 m，150 m]，將區域劃分為尺寸為1.5 m的網格，每個網格節點與武器爆點之間的相對位置是確定的。根據武器終點狀態調整破片場仿真模型得到動態破片場分布，然后按照圖3中循環部分流程，分別計算人體模型在所有網格節點位置處的失能概率，可得到對應圖6武器終點狀態下單發破片戰斗部對人員目標在補給12 h殺傷定義下的毀傷矩陣，如圖7(a)所示。

圖7(a)中網格單元顏色越明亮則表示該位置處失能越嚴重，各網格中的失能概率由破片飛散情況、衰減速度和命中人體部位以及命中破片數目共同決定。從毀傷矩陣給出的概率分布情況可以直觀判斷破片戰斗部對人員目標威脅較大的區域，從而為戰場人員防護及靶場試驗設計提供參考。

..終點姿態對毀傷效能的影響

變化戰斗部的終點狀態可以得到更多打擊場景下的毀傷矩陣。設計單發破片戰斗部的4種終點狀態工況，如表3所示，按照圖3展示的流程，計算得到與4種工況在補給 12 h殺傷定義下對應的改進毀傷矩陣，如圖7所示。

表3 武器終點狀態工況信息

圖7(b)中落角取為85°，此時毀傷區域范圍變得更加集中，由于破片場散布均勻且彈體接近垂直下落，概率分布也相對均勻；圖7(c)中武器發生90°偏航，毀傷矩陣也整體隨之旋轉了90°； 圖7(d)武器滾轉角發生變化，毀傷矩陣并未出現明顯變化，僅在部分區域有細微差別。可見，對于單發破片戰斗部而言，可以通過調整武器終點姿態的方式，實現既滿足一定的火力壓制范圍又避免特定方向上附帶毀傷的多目標優化。

圖7 特定終點姿態下單彈的毀傷矩陣Fig.7 Damage matrix of a single projectile with a specific terminal attitude

..“平均”毀傷矩陣

由圖7不難看出，落角和偏航角對毀傷矩陣影響較大，滾轉角對毀傷矩陣的影響大小與破片場分布有關，武器終點姿態的變化會使得同一型號的單發破片戰斗部對人員目標的毀傷效能也差異很大。如果再加上武器精度給爆點位置帶來的隨機性，評估單彈對人員的毀傷效能僅依賴特定終點狀態下的毀傷矩陣是不夠的，利用Monte Carlo方法可以通過大量迭代得到“平均”的結果。

本文基于Monte Carlo方法對武器精度進行了模擬，引入隨機終點姿態，假設落角范圍為60°～90°，以5°為間隔遞增，只考慮滾轉角為無滾轉、90°、180°和270°等4種情況，偏航角在0°～360°范圍內隨機，設定武器CEP為10 m和50 m，目標區域網格劃分相同，且爆高10 m和落速500 m/s不變。

圖8給出了補給12 h殺傷定義下對2種武器精度進行10 000次隨機抽樣的平均毀傷矩陣，失能概率分布由不規則的散布轉變為以預期彈著點中心為圓心的圓，且圓心附近是對人員威脅相對較大的區域。當CEP較小時，平均毀傷矩陣會受破片飛散產生的中心毀傷盲區影響，在中心區域形成一個概率低谷，精度下降會消除這一影響，但會影響平均毀傷效果。平均毀傷矩陣的失能概率峰值較之單次特定條件的打擊結果大幅降低，且隨著CEP的增大，下降會更為嚴重。

圖8 單彈的“平均”毀傷矩陣

5 毀傷效能量化指標

5.1 平均有效毀傷面積(MAEF)

盡管毀傷矩陣可以直觀表示出武器對不同位置處人員的威脅程度，但是仍需要一種能夠表征一枚彈毀傷效能的量化指標。在武器毀傷效能評估中，面積是一種常見的指標形式，Driels通過破片的平均有效毀傷面積(MAE)表征了破片戰斗部的毀傷效能，可將其變形為：

(5)

式(5)中：MAE為平均有效毀傷面積；(，)為單元處失能概率；、、、為毀傷范圍。

MAE實際上是對毀傷范圍內單元面積的加權求和，權重由單元處的失能概率決定，于是基于毀傷矩陣就可以快速得到平均有效毀傷面積。

5.2 MAEF計算結果

根據式(5)，可計算獲得不同打擊方式下的MAE，基于圖7和圖8所示毀傷矩陣得到的MAE，計算結果見表4。對于不考慮武器精度的單發破片戰斗部特定打擊情況，當落角變大時，更集中的毀傷范圍導致MAE出現明顯下降，而滾轉角和偏航角的改變對MAE影響不大，這與根據毀傷矩陣得到的直觀判斷一致；考慮武器精度和隨機終點狀態，“平均”毀傷矩陣的MAE相比圖7工況的毀傷矩陣下降較多，CEP分別為10 m和50 m的MAE相差也很明顯，兩者MAE之間的差異可能是由于較差的武器精度使得部分零散破片飛出目標區域范圍導致的。

表4的數據表明，特定終點狀態下的單彈毀傷效能與隨機條件下平均的毀傷效能有較大差別，有些計算條件下，MAE下降可達30%，因此實戰條件下的火力規劃單純基于武器的威力參數是不夠的。

表4 MAEF結果

6 結論

通過改進毀傷矩陣的計算方法對人員目標進行了毀傷效能評估，得到以下結論：

1) 本文提出的基于精細模型和A-S準則進行毀傷矩陣計算的改進方法，可為破片戰斗部對人員目標的毀傷效能評估提供更準確的失能概率分布。

2) 改進毀傷矩陣的計算結果表明，單發破片戰斗部打擊下呈現出嚴重威脅人員安全的區域，武器的終點姿態對毀傷矩陣有較大影響，不同終點狀態的毀傷矩陣可以為武器低附帶毀傷規劃運用提供指導。

3) 基于Monte Carlo方法模擬武器的隨機命中以反映實戰條件，所得到的平均毀傷矩陣的毀傷概率峰值會隨著CEP的增大而大幅降低。利用MAE表征單彈的毀傷效能，特定狀態下的毀傷矩陣計算得到的MAE與平均毀傷矩陣得到的MAE差別明顯，說明了實戰打擊效果與武器的理想毀傷性能有差距。所提出的計算流程可以服務于破片戰斗部對人員目標的實戰毀傷效能評估。