CMO武器与毁伤模型-子母弹覆盖与毁伤模型实现算法详解

5.2.1 子母弹覆盖与毁伤
子母弹的建模区别于常规战斗部：它不是在一点释放能量，而是在一片区域内按面密度施加毁伤。
算法：子弹药覆盖面积 ClusterCoverageArea
Area_base = 51.369 863 013 698 63 × 子弹药数量 [m²]
按子弹药类型修正：
Cluster_AP（杀伤爆破型） × 1.2
Cluster_AT（反装甲型） × 0.8
Cluster_Penetrator（侵彻型） × 0.4

def cluster_coverage_area(num_bomblets: int, bomblet_type: str) –> float:
"""
计算子母弹对目标的覆盖面积（m²）。

输入参数：
num_bomblets : int —— 子弹药数量（枚）
bomblet_type : str —— 子弹药类型，取值：
"Cluster_AP"（杀伤爆破型）
"Cluster_AT"（反装甲型）
"Cluster_Penetrator"（侵彻型）

计算公式：
Area_base = 51.36986301369863 × num_bomblets
Area = Area_base × type_factor

返回结果：
float —— 覆盖面积（m²）
"""
# 单枚子弹药的有效杀伤面积基数（m²）
area_base = 51.36986301369863 * num_bomblets

# 按子弹药类型修正覆盖面积
if bomblet_type == "Cluster_AP":
# 杀伤爆破型：靠破片与冲击波覆盖，等效面积最大
type_factor = 1.2
elif bomblet_type == "Cluster_AT":
# 反装甲型：需直接命中装甲薄弱部位，等效覆盖面积减小
type_factor = 0.8
elif bomblet_type == "Cluster_Penetrator":
# 侵彻型：需精确命中目标才能起爆，等效覆盖面积最小
type_factor = 0.4
else:
raise ValueError(f"未知子弹药类型: {bomblet_type}")

return area_base * type_factor

基数 51.37 m² 对应单枚子弹药的有效杀伤面积；反装甲与侵彻型子弹药因需要直接命中，等效覆盖面积依次减小。
算法：对目标的覆盖毁伤 BombletDamageAtThisTarget
若目标为机场跑道类线状目标：
Damage_fraction = 0.1 + 0.9 · ( 1 − |sin(Δheading)| )
Δheading 为投弹航向与跑道走向的夹角，
垂直穿越跑道时仅得 0.1，沿跑道纵向投放时得 1.0
否则（面状目标）：
Damage_fraction = FlatSurfaceArea_m2 / ClusterCoverageArea

import math

def bomblet_damage_at_this_target(
target_is_runway: bool,
delta_heading_deg: float,
flat_surface_area_m2: float,
cluster_coverage_area: float,
) –> float:
"""
计算子母弹对目标的覆盖毁伤比例（0~1）。

输入参数：
target_is_runway : bool —— 目标是否为机场跑道类线状目标
delta_heading_deg : float —— 投弹航向与跑道走向的夹角（度）
flat_surface_area_m2 : float —— 面状目标的有效面积（m²）
cluster_coverage_area : float —— 子母弹覆盖面积（m²），由
cluster_coverage_area() 计算得到

计算公式：
线状目标（跑道）：
Damage_fraction = 0.1 + 0.9 × (1 − |sin(Δheading)|)
面状目标：
Damage_fraction = FlatSurfaceArea_m2 / ClusterCoverageArea

返回结果：
float —— 毁伤比例（0~1）
"""
if target_is_runway:
# 线状目标（跑道）：
# 0.1 —— 横穿跑道时的最小毁伤基数（仅单点破坏）
# 0.9 —— 纵向投放可提升的最大毁伤增量
# |sin(Δheading)| —— 航向与跑道夹角的正弦绝对值，
# 夹角越大（越接近横穿）毁伤越小
delta_heading_rad = math.radians(delta_heading_deg)
damage_fraction = 0.1 + 0.9 * (1.0 – abs(math.sin(delta_heading_rad)))
else:
# 面状目标：毁伤比例 = 目标有效面积 / 子母弹覆盖面积
# 当覆盖面积不足时毁伤比例可能大于 1，由调用方截断
damage_fraction = flat_surface_area_m2 / cluster_coverage_area

return damage_fraction

第一式的物理含义是：子母弹沿跑道纵向投放才能形成连续弹坑链，横穿跑道只能造成单点破坏。这是封锁机场战术的模型依据。
算法：对瞄准点的毁伤 BombletDamageAtThisAimpoint

import random

def bomblet_damage_at_this_aimpoint(
num_bomblets: int,
bomblet_type: str,
target_is_runway: bool,
delta_heading_deg: float,
flat_surface_area_m2: float,
armor_level: int,
) –> float:
"""
计算子母弹对瞄准点的最终毁伤比例（0~1）。

输入参数：
num_bomblets : int —— 子弹药数量（枚）
bomblet_type : str —— 子弹药类型，取值：
"Cluster_AP"（杀伤爆破型）
"Cluster_AT"（反装甲型）
"Cluster_Penetrator"（侵彻型）
target_is_runway : bool —— 目标是否为机场跑道类线状目标
delta_heading_deg : float —— 投弹航向与跑道走向的夹角（度）
flat_surface_area_m2 : float —— 面状目标的有效面积（m²）
armor_level : int —— 目标装甲等级（1~4，由低到高）

计算公式：
Damage_fraction = BombletDamageAtThisTarget × type_match × armor_factor × random(0.8, 1.2)

返回结果：
float —— 对瞄准点的最终毁伤比例（0~1）
"""
# 第一步：调用覆盖毁伤算法，得到基础毁伤比例
coverage_area = cluster_coverage_area(num_bomblets, bomblet_type)
base_damage = bomblet_damage_at_this_target(
target_is_runway,
delta_heading_deg,
flat_surface_area_m2,
coverage_area,
)

# 第二步：子弹药类型与目标匹配系数
if target_is_runway and bomblet_type == "Cluster_Penetrator":
# 跑道侵彻型子弹药对跑道目标：×5
type_match = 5.0
else:
# 其余组合不额外加成
type_match = 1.0

# 第三步：目标装甲等级匹配系数（由低到高依次 ×1.5、×2、×3、×4）
armor_factor_map = {1: 1.5, 2: 2.0, 3: 3.0, 4: 4.0}
if armor_level not in armor_factor_map:
raise ValueError(f"未知装甲等级: {armor_level}，应为 1~4")
armor_factor = armor_factor_map[armor_level]

# 第四步：0.8–1.2 均匀随机数，体现散布随机性
scatter_random = random.uniform(0.8, 1.2)

# 最终毁伤比例（截断到 0~1）
final_damage = base_damage * type_match * armor_factor * scatter_random
return max(0.0, min(1.0, final_damage))

# ===== 实战调用示例 =====
if __name__ == "__main__":
# 示例 1：跑道侵彻型子弹药攻击机场跑道（装甲等级视为 1）
damage_runway = bomblet_damage_at_this_aimpoint(
num_bomblets=200,
bomblet_type="Cluster_Penetrator",
target_is_runway=True,
delta_heading_deg=5.0, # 近似沿跑道纵向投放
flat_surface_area_m2=0.0, # 线状目标，此参数不参与计算
armor_level=1,
)
# 运行结果示例：0.98（纵向投放 + 侵彻型 ×5 + 随机数接近 1.2）

# 示例 2：杀伤爆破型子弹药攻击面状装甲目标（装甲等级 3）
damage_area = bomblet_damage_at_this_aimpoint(
num_bomblets=300,
bomblet_type="Cluster_AP",
target_is_runway=False,
delta_heading_deg=0.0, # 面状目标，此参数不参与计算
flat_surface_area_m2=8000.0, # 目标有效面积 8000 m²
armor_level=3,
)
# 运行结果示例：0.87（覆盖面积充足 + 装甲等级 ×3 + 随机数约 1.0）

print(f"跑道目标毁伤比例: {damage_runway:.2f}")
print(f"面状目标毁伤比例: {damage_area:.2f}")

在上式基础上叠加子弹药类型与目标装甲的匹配系数：跑道侵彻型子弹药对跑道 ×5；按目标装甲等级由低到高依次 ×1.5、×2、×3、×4；最后乘以 0.8–1.2 的均匀随机数以体现散布随机性。