基于cube模型的碰撞概率计算方法技术

本公开提供了一种基于cube模型的碰撞概率计算方法，将近地空间划分立方体网格，并对立方体进行编号；计算空间物体的位置及其所处的立方体编号；对空间物体进行配对，将处于相邻和相同立方体内的空间物体作为配对目标；计算配对目标之间的碰撞概率。本公开通过增加对相邻cube空间物体的碰撞可能性的考虑，弥补了原有cube模型可能遗漏碰撞情况的缺陷，使得改进后的cube模型更符合物理真实，分析结果更加合理可信。

【技术实现步骤摘要】
基于cube模型的碰撞概率计算方法
本公开涉及空间物体碰撞概率计算领域，尤其涉及一种基于cube模型的碰撞概率计算方法。
技术介绍
空间碎片的日益增多引发了国际上关于空间环境长期稳定性的研究。而计算机技术的发展使得对N体系统的长期演化仿真成为可能。在空间碎片系统演化中，需要考虑碎片的增加机制和减少机制，与之相应的，需要一种能够利用不断更新的碎片轨道数据进行碰撞概率计算的算法。美国国家宇航局研制了一种称为“Cube”模型的碰撞概率算法。Cube模型的核心思想是，在空间碎片系统长期演化中，每隔一定时间dt对系统采样一次，在采样时刻，将整个近地空间划分为长度h的小立方体，利用更新的轨道数据计算出该时刻每个空间物体所在的位置，然后找出处于同一个立方体的空间物体，计算其两两之间的碰撞概率，若一个立方体中只有一个物体，则不考虑其与其他立方体中物体的碰撞概率。Cube模型的优点在于通过对空间划分立方体可以对N体碎片系统进行快速配对，大大减少了计算量；其缺点则是遗漏了对处于相邻立方体中的空间物体之间碰撞概率的计算，这种碰撞概率甚至会大于同一立方体中两个物体之间的碰撞概率。
技术实现思路
(一)要解决的技术问题本公开旨在弥补Cube模型中没有考虑相邻立方体中空间物体之间碰撞概率的缺陷，提出一种基于cube模型的碰撞概率计算方法。(二)技术方案本公开提供了一种基于cube模型的碰撞概率计算方法，包括：将近地空间划分立方体网格，并对立方体进行编号；计算空间物体的位置及其所处的立方体编号；对空间物体进行配对，将处于相邻和相同立方体内的空间物体作为配对目标；以及计算配对目标之间的碰撞概率。在本公开的一些实施例中，所述将近地空间划分立方体网格，并对立方体进行编号包括：以地心为坐标系原点，建立笛卡尔坐标系，将近地空间划分为立方体网格，并对每个立方体以边长为单位进行编号。在本公开的一些实施例中，所述编号表示为(x，y，z)，x、y、z分别表示立方体在x轴方向、y轴方向和z轴方向的编号。在本公开的一些实施例中，所述计算空间物体的位置及其所处的立方体编号包括：根据每个空间物体更新后的轨道根数计算其所在的位置，查找所述位置所处的立方体，并记录所处立方体的编号。在本公开的一些实施例中，所述对空间物体进行配对，将处于相邻和相同立方体内的空间物体作为配对目标包括：分别在三个坐标轴方向，按序号顺序将空间物体所处立方体的编号进行排序，得到与目标空间物体所处立方体的编号相邻和相同的编号；与所述编号对应的立方体内的空间物体序号构成三个坐标轴方向的空间物体序号集；计算三个坐标轴方向的空间物体序号集的交集，得到与目标空间物体所处立方体相邻和相同的立方体的空间物体序号集；遍历空间物体序号集中的空间物体与目标空间物体的欧氏距离，保留欧式距离在小于阈值的空间物体作为配对目标。在本公开的一些实施例中，利用二分法按序号顺序将空间物体所处立方体的编号进行排序。在本公开的一些实施例中，所述阈值取立方体边长乘以在本公开的一些实施例中，所述计算配对目标之间的碰撞概率包括：对处于相邻和相同立方体内的空间物体，空间物体i和j的平均碰撞数为：c＝SiSjVimpAcdUdt碰撞概率为：pij＝1-exp(-c)其中，dU为体元；dt为时间采样间隔；Si、Sj分别表示空间物体i和j在体元dU内的分布密度；Vimp是空间物体i和j的相对碰撞速度；Ac为空间物体i和j的碰撞截面。在本公开的一些实施例中，所述体元dU是以空间物体i和j之间欧式距离的阈值为半径的球体体积。在本公开的一些实施例中，其中，空间物体在体元dU内均匀分布，Si＝Sj＝1/dU其中，dc为空间物体i和j之间欧式距离的阈值；h为立方体边长。(三)有益效果从上述技术方案可以看出，本公开基于cube模型的碰撞概率计算方法具有以下有益效果：(1)本公开通过增加对相邻cube空间物体的碰撞可能性的考虑，弥补了原有cube模型可能遗漏碰撞情况的缺陷，使得改进后的cube模型更符合物理真实，分析结果更加合理可信。(2)本公开适用于空间碎片系统在长期演化过程中对碰撞概率的评估，有利于提高空间碎片长期演化模型的合理性与可信度。而对空间碎片长期演化模型的研究能够对未来我国进行主动碎片清除提供策略分析和模型支撑。附图说明图1是本公开实施例的基于cube模型的碰撞概率计算方法的流程图。具体实施方式本公开通过对空间物体快速配对，计算空间物体之间的碰撞概率。该方法适用于空间碎片系统在长期演化过程中对碰撞概率的评估。在本公开中，考虑空间目标之间的碰撞，可将其视为空气分子间的碰撞，依据气体动力学来评估空间目标之间的碰撞概率。空气分子间的碰撞概率服从泊松分布。为使本公开的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本公开进一步详细说明。本公开实施例提供了一种基于cube模型的碰撞概率计算方法，参见图1，包括：步骤S1：将近地空间划分立方体网格，并对立方体进行编号。该步骤具体包括：在采样时刻，以地心为坐标系原点，建立笛卡尔坐标系，将地球周围的近地空间划分为边长为h的立方体形成的网格，并对每个立方体按边长h的倍数进行编号，即以边长h为单位进行编号。立方体的编号表示为(x，y，z)，x、y、z分别表示立方体在x轴方向、y轴方向、z轴方向的编号；x轴方向编号表示该立方体处于x轴的第x个单位的位置，同理，y轴方向编号表示该立方体处于y轴的第y个单位的位置，y轴方向编号表示该立方体处于z轴的第z个单位的位置。步骤S2：计算空间物体的位置及其所处立方体的编号。该步骤具体包括：在采样时刻，根据每个空间物体更新后的轨道根数计算其所在的位置，查找该位置所处的立方体，并记录所处立方体的编号。在该步骤中可先对空间物体进行排序，对于序号为i的空间物体，其位置坐标表示为(Xi，Yi，Zi)，该位置所处立方体的编号表示为(xi，yi，zi)。步骤S3：对空间物体进行配对，将处于相邻和相同立方体内的空间物体作为配对目标。该步骤是在采样时刻，通过对所有空间物体所处立方体的编号进行比对，查找出处于相邻或者相同立方体内的空间目标。具体来说，首先，分别在三个坐标轴方向，按序号顺序将空间物体所处立方体的编号进行二分法排序，得到与目标空间物体所处立方体的编号相邻和相同的编号。与上述编号对应的立方体内的空间物体序号构成三个坐标轴方向的空间物体序号集。在一个示例中，可以先按空间物体的序号顺序，将各个空间物体所处立方体的x轴方向编号进行二分法排序，得到空间物体所处立方体的x轴方向编号集，记录x轴方向编号集中与空间物体i所处立方体的x轴方向编号相邻和相同的x轴方向编号，得到与这些x轴方向编号对应的立方体的集合，该集合中的立方体内的空间物体序号构成x轴方向空间物体序号集Mx。此处的二分法排序具体包括：在对第k个空间物体所处立方体的x轴方向编号进行排序时，之前的k-1个空间物体所处立方体的x轴方向编号已经按照x轴排好序并形成数组Nk-1，此时首先利用二分法查找第k个空间物体所处立方体的x轴方向编号在数组Nk-1位置，然后将第k个空间物体所处立方体的x轴方向编号插入Nk-1的对应位置得到数组Nk。同样，按序号顺序将空间物体所处立方体的y轴方向编号进行二分法排序，得到空间物体本文档来自技高网...

【技术保护点】
一种基于cube模型的碰撞概率计算方法，包括：将近地空间划分立方体网格，并对立方体进行编号；计算空间物体的位置及其所处的立方体编号；对空间物体进行配对，将处于相邻和相同立方体内的空间物体作为配对目标；以及计算配对目标之间的碰撞概率。

【技术特征摘要】
1.一种基于cube模型的碰撞概率计算方法，包括：将近地空间划分立方体网格，并对立方体进行编号；计算空间物体的位置及其所处的立方体编号；对空间物体进行配对，将处于相邻和相同立方体内的空间物体作为配对目标；以及计算配对目标之间的碰撞概率。2.如权利要求1所述的碰撞概率计算方法，所述将近地空间划分立方体网格，并对立方体进行编号包括：以地心为坐标系原点，建立笛卡尔坐标系，将近地空间划分为立方体网格，并对每个立方体以边长为单位进行编号。3.如权利要求2所述的碰撞概率计算方法，所述编号表示为(x，y，z)，x、y、z分别表示立方体在x轴方向、y轴方向和z轴方向的编号。4.如权利要求1所述的碰撞概率计算方法，所述计算空间物体的位置及其所处的立方体编号包括：根据每个空间物体更新后的轨道根数计算其所在的位置，查找所述位置所处的立方体，并记录所处立方体的编号。5.如权利要求1所述的碰撞概率计算方法，所述对空间物体进行配对，将处于相邻和相同立方体内的空间物体作为配对目标包括：分别在三个坐标轴方向，按序号顺序将空间物体所处立方体的编号进行排序，得到与目标空间物体所处立方体的编号相邻和相同的编号；与所述编号对应的立方体内的空间物体序号构成三个坐标轴方向的空间物体序号集；计算三个坐标轴方向的空间物体序...

【专利技术属性】
技术研发人员：王晓伟，刘静，崔双星，杨志涛，
申请(专利权)人：中国科学院国家天文台，
类型：发明
国别省市：北京,11