[发明专利]多次抛撒箔条云气动融合轨道再入实时跟踪统计模拟方法

权利要求书

1.一种多次抛撒箔条云气动融合轨道再入实时跟踪统计模拟方法，其特征在于，对于每批次抛撒的箔条个体形成箔条云进行模拟，连续对批次形成的箔条云进行加权统计模拟，以得到对应的批次抛撒整体箔条云覆盖范围，经统计加权得到对应的整体箔条云加权椭球体；

基于每个箔条云加权椭球体均是由有限的组箔条个体团所支撑，故在任意时刻均有组箔条个体团被气动融合轨道再入实时跟踪、随机统计，进而通过一次模拟实现对飞行器上升段多批次抛撒释放箔条云的气动特性、运动性能进行沿轨道再入计算分析；

所述实时跟踪、随机统计的方法被配置为包括以下步骤：

S1、实时计算箔条云中心位置及整体飘落飞行速度；

S2、计算确定整体箔条云雷达电磁散射截面；

S3、循环推进，实时跟踪模拟箔条云；

其中，在S3中，将S1中得出的整体箔条云中心位置、整体飘落飞行速度对应的物理量统一到发射坐标系中，进行弹道计算，进而确定任意时刻任一个箔条个体小组飘落飞行的空间位置与速度，实现实时跟踪；

在S2中，所述整体箔条云雷达电磁散射截面的计算方式被配置为包括：

S21、将第一批抛撒释放的箔条个体所形成的箔条云加权球体作为初始球形分层抛撒释放的所有箔条个体所形成的整体箔条云覆盖范围，进行箔条云切片；

S22、基于整体箔条云被划分的各个切片，对初始抛撒释放的各批次箔条个体层进行循环，并在各层箔条云切片占据两个或两个以上更小的箔条云加权球时，对该箔条云切片中的箔条个体数进行累加，并将该层对应的较小箔条云切片中的箔条数、沿雷达视线方向箔条个体径向速度和箔条个体相对于发射坐标系的俯仰角、偏航角、滚转角，均累加到该切片中，按预定条件进行统计平均处理，以确定最大径向速度、最小径向速度以及平均径向速度分布；

S23、在所有切片层均累加统计完后，得到按要求划分的初始抛撒释放的各批次箔条个体飘落飞行所形成的整体箔条云各切片中箔条个体空间、姿态、速度分布，最终确定整体箔条云雷达电磁散射截面；

在S3中，所述循环推进的方式被配置为包括：

S31、基于给定的头体分离点相关参数，将箔条个体进行轨道计算推进，当计算时间达到首批箔条释放时间，所释放的箔条个体被分成抛撒释放的飞行姿态小组；

S32、使用随机映射方法将与所划分的相应飞行攻角、侧滑角小组建立随机映射关系，确定各箔条小组的飞行攻角与侧滑角(，)；

S33、结合各组箔条个体所具有的速度模拟预测各箔条个体团在下一时刻的空间位置与飞行速度，并按雷达电磁散射特性计算所要求的时间间隔进行循环推进；

S34、在需要统计分析箔条云整体性能的每一时刻，对组箔条个体团中各组箔条数目进行重新调整确定，模拟各时间间隔箔条个体的飞行姿态连续变化；

S35，对于未抛撒释放的箔条个体，仍按飞行器飞行攻角与侧滑角进行计算推进，直到各批次的箔条个体均释放完后；

S36、以时间间隔为步长的箔条云计算推进循环，各自独立地跟踪个组箔条个体团，当计算推进到一定时刻，发现(，)下箔条个体飘落飞行速度已降为零或已达到预先设定的某最低飞行高度，则停止该箔条个体团所有箔条个体沿轨道循环计算，对其余箔条个体团继续进行气动融合轨道计算循环推进，直到所有飞行攻角与侧滑角下箔条个体飘落飞行速度接近或趋于零，循环推进结束，得到该时刻箔条云覆盖范围和不同姿态角箔条个体在空中悬浮分布情况。

2.如权利要求1所述的多次抛撒箔条云气动融合轨道再入实时跟踪统计模拟方法，其特征在于，在S1中，所述箔条云中心位置及整体飘落速度是第L层中箔条个体飘落飞行形成的箔条云形心和速度，再对各批次进行统计加权，可计算得到由飞行器释放出的总数为的所有各层中箔条个体飘落飞行形成的整体箔条云加权形心和加权速度。