CN107066641B - 大规模空间碎片分布演化的数值计算方法及系统 - Google Patents
大规模空间碎片分布演化的数值计算方法及系统 Download PDFInfo
- Publication number
- CN107066641B CN107066641B CN201611036194.6A CN201611036194A CN107066641B CN 107066641 B CN107066641 B CN 107066641B CN 201611036194 A CN201611036194 A CN 201611036194A CN 107066641 B CN107066641 B CN 107066641B
- Authority
- CN
- China
- Prior art keywords
- space debris
- space
- debris
- fragments
- distribution
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 94
- 238000000034 method Methods 0.000 title claims abstract description 44
- 239000012634 fragment Substances 0.000 claims abstract description 196
- 238000004880 explosion Methods 0.000 claims abstract description 34
- 238000012545 processing Methods 0.000 claims abstract description 11
- 238000004364 calculation method Methods 0.000 claims description 56
- 230000010354 integration Effects 0.000 claims description 49
- 238000005315 distribution function Methods 0.000 claims description 19
- 230000001133 acceleration Effects 0.000 claims description 18
- 238000012937 correction Methods 0.000 claims description 13
- 239000013598 vector Substances 0.000 claims description 11
- 239000002360 explosive Substances 0.000 claims description 10
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 claims description 8
- 238000013467 fragmentation Methods 0.000 claims description 7
- 238000006062 fragmentation reaction Methods 0.000 claims description 7
- 230000005484 gravity Effects 0.000 claims description 6
- 230000014509 gene expression Effects 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 5
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 claims description 4
- 238000011423 initialization method Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 238000004134 energy conservation Methods 0.000 claims description 3
- 238000009827 uniform distribution Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 230000001351 cycling effect Effects 0.000 claims description 2
- 238000009795 derivation Methods 0.000 claims description 2
- 238000009877 rendering Methods 0.000 claims description 2
- 230000007774 longterm Effects 0.000 abstract description 5
- 230000000116 mitigating effect Effects 0.000 abstract description 3
- 238000009472 formulation Methods 0.000 abstract description 2
- 239000000203 mixture Substances 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 8
- 230000008859 change Effects 0.000 description 7
- 230000009471 action Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 241000282414 Homo sapiens Species 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000001932 seasonal effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
本发明提供一种大规模空间碎片分布演化的数值计算方法及系统,系统包括:空间碎片状态推演子系统、空间碎片碰撞/爆炸解体产生新碎片子系统、空间碎片运动状态三维显示子系统、空间碎片运动状态输出子系统和空间碎片演化结果处理子系统。本发明提供的大规模空间碎片分布演化的数值计算方法及系统具有以下优点:本发明提供的大规模空间碎片分布演化的数值计算方法及系统,能够对大规模空间碎片分布状态进行长期演化计算,能够动态三维演示碎片分布状态,提供了碎片分布状态数据的输出和处理分析功能,可以为航天器在轨运行空间碎片碰撞规避、空间碎片减缓措施制定以及确保空间资源可持续利用的重要的保障。
Description
技术领域
本发明属于航天动力学技术和计算机技术领域,具体涉及一种大规模空间碎片分布演化的数值计算方法及系统。
背景技术
随着人类空间活动的日益增加,人为发射进入空间的物体不断增多,使得地球外部空间的物体数量迅速增加,在轨航天器受到的碰撞威胁也愈加严重。
据NASA空间监视网(SSN)的统计,地球同步轨道高度以下的大于10cm尺寸的空间物体数量达到15000多个,小于10cm尺寸的空间物体数量更多,粗略估计至少为400000个,致使地球外部空间达到了前所未有的拥挤状态。这些空间物体中,只有很少的一部分是在轨正常工作的航天器,大部分是由于发射任务、在轨航天器爆炸分解、在轨航天器碰撞分解产生的空间碎片。因此,这些在轨道动力学的约束下沿轨道绕地球自由运动的空间碎片,会挤占有限的空间资源,使得诸如近地轨道和静止轨道等区域可用的资源减少。同时,在大气阻力、地球非球形摄动以及太阳、月球三体引力等摄动力的作用下,自由无控的空间碎片轨道会随时间不断变化,当空间碎片轨道与正常工作航天器的轨道存在交点时,就会存在与航天器发生碰撞的可能,对正常工作航天器产生威胁。另外,航天器发生碰撞后,会进一步使空间碎片的数量增加,从而使得产生新的碰撞的可能性增加,进而产生更多碎片,对正常运行航天器产生更大威胁,这种链式的雪崩式效应,最终会导致航天器在轨运行的成本迅速增加,甚至会出现近地或同步轨道范围内航天器均无法生存的情况。
由此可见,大量空间碎片对航天器造成的碰撞威胁,将是人类空间活动面临的主要安全问题,如不加以管理,将导致航天器运行成本大大增加,甚至会使空间资源无法利用。有效的对当前空间碎片的运动分布状态进行预测分析,是航天器在轨运行碰撞规避、空间碎片减缓措施制定以及确保空间资源可持续利用的重要的技术保障。现有技术中,尚未出现有效的对当前空间碎片的运动分布状态进行预测分析的相关内容。
发明内容
针对现有技术存在的缺陷,本发明提供一种大规模空间碎片分布演化的数值计算方法及系统,可有效解决上述问题。
本发明采用的技术方案如下:
本发明提供一种大规模空间碎片分布演化的数值计算方法,包括以下步骤:
S1,空间碎片状态推演子系统获取空间碎片的初始运动状态数据以及空间碎片受到的空间摄动力,并将所述初始运动状态数据以及所述空间摄动力输入到数值积分模型中;所述数值积分模型按照给定的时间步长进行积分推演,得到空间碎片未来运动状态数据;然后,采用所述空间碎片未来运动状态数据更新空间碎片当前运动状态数据;
S2,判断是否达到演化时间,如果未达到,将状态更新后的空间碎片运动状态和属性数据传输给空间碎片碰撞/爆炸解体产生新碎片子系统;
S3,所述空间碎片碰撞/爆炸解体产生新碎片子系统将所述状态更新后的空间碎片运动状态和属性数据输入到碎片解体模型中;所述碎片解体模型进行解体计算,得到空间碎片解体后产生新空间解体碎片的属性数据;并将所述新空间解体碎片的属性数据反馈传输到所述空间碎片状态推演子系统;
S4,所述空间碎片状态推演子系统将所述新空间解体碎片的属性数据以及已有的空间碎片的运动状态数据输入到数值积分模型中,再次按照给定的时间步长进行积分推演,得到空间碎片未来运动状态数据,如此不断循环S2-S4;直到达到演化时间,结束流程。
优选的,S1中,所述空间摄动力包括大气阻力摄动力、地球非球形引力摄动力、太阳光压摄动力以及第三体引力摄动力;所述数值积分模型为4阶阿达姆斯预估矫正积分模型。
优选的,S3中,所述空间碎片碰撞/爆炸解体产生新碎片子系统具体用于:
首先,根据能量定律确定空间碎片解体后产生新空间解体碎片的数量N和新空间解体碎片特征尺寸l的分布;
然后,以新空间解体碎片特征尺寸l作为独立变量,利用概率分布模型确定新空间解体碎片的面质比A/md、新空间解体碎片的迎风截面积A、新空间解体碎片相对于解体前空间碎片的速度增量Δv;
最后,根据A/md、A确定碎片质量md。
优选的,S3中,所述空间碎片碰撞/爆炸解体产生新碎片子系统具体用于:
S3.1,确定解体碎片数量:
其中,解体碎片区分为爆炸解体碎片和碰撞解体碎片;
1)对于爆炸解体碎片,根据能量守恒定律,解体碎片特征尺寸l大于被研究解体碎片最小特征尺寸lc的碎片数量Nf满足函数关系:
其中:为解体碎片的最小特征尺寸;
cs是修正系数,与爆炸类型有关,对于历史爆炸事件,修正系数cs的表达式:
其中:
其中:H为爆炸解体碎片的高度;
2)对于碰撞解体碎片,碰撞解体碎片数量满足函数关系:
式中:为两碰撞目标的等效质量:
其中:
其中:Ep为等效碰撞能量,为灾难性碰撞临界能量;mp为相互碰撞中较小质量;mt为相互碰撞中较大质量;vr是碰撞时两碎片的相对速度,即碰撞速度;
S3.2,确定解体碎片尺寸分布律:
爆炸解体碎片尺寸分布律F:
碰撞解体碎片尺寸分布律F:
S3.3,确定解体碎片面质比参数:
本步骤中,爆炸解体碎片和碰撞解体碎片的碎片面质比参数满足同样的分布函数;
(1)对于尺寸大于11cm解体碎片,定义解体碎片面质比的对数值γ=log(A/md),解体碎片特征长度的对数值θ=log(l),解体碎片的面质比分布由如下二元正态分布决定:
p(γ,θ)=ε(θ)p1(γ)+(1-ε(θ))p2(γ) (19)
其中:p1(γ)和p2(γ)均是正态分布概率密度函数:
其中:权重参数ε(θ)∈[0,1],是θ的饱和函数;正态分布概率密度函数的均值μk和正态分布概率密度函数的方差σk均是θ的函数;
对于航天器解体碎片,其面质比分布函数中的参数由下式决定:
对于火箭上面级解体碎片,其面质比分布函数的参数由下式决定
(2)对于尺寸小于8cm的航天器解体碎片和尺寸小于1.7cm的火箭上面级解体碎片,碎片的面质比满足正态分布律p1(γ),即权重函数ε=1,
p(γ,θ)=p1(γ) (23)
此时,分布函数参数由下式统一决定:
(3)对于尺寸介于8cm~11cm之间的航天器解体碎片和尺寸介于1.7cm~11cm之间的火箭上面级碎片,碎片的面质比通过随机采样的方法确定,即:首先生成服从均匀分布且取值区间在[0.0,1.0]内的随机变量ζ,然后将ζ与利用(25)计算得到的进行对比,若则利用公式(19)计算面质比参数,反之则利用公式(23)计算面质比参数;
步骤3.4,确定解体碎片的速度增量:
解体碎片相对于解体前目标的速度增量满足正态分布规律,定义v=lg(Δv),γ=lg(A/md),则v的分布函数p(ν)满足
其中:v为解体碎片的速度;Δv为解体碎片的速度增量;
解体碎片的均值μν和解体碎片的方差σν通过下式计算:
优选的,还包括:
S5,空间碎片运动状态三维显示子系统将所述空间碎片状态推演子系统推演到的空间碎片未来运动状态数据进行三维显示,即:动态的显示所有空间碎片在空间中的位置分布状态;具体包括:
步骤S5.1,接收空间碎片状态推演子系统推演到的空间碎片未来运动状态数据;
步骤S5.2,调用三维显示接口函数,将空间碎片的位置和形状在三维模型中进行渲染;
步骤S5.3,驱动三维模型,对空间碎片运动状态进行动态更新。
优选的,还包括:
S6,在固定时刻或者设定的时刻,空间碎片运动状态输出子系统将所述空间碎片状态推演子系统推演到的空间碎片未来运动状态输出;其中,所述空间碎片未来运动状态包括:空间碎片的位置、速度矢量和面值比、质量、来源属性参数;此外,空间碎片运动状态输出子系统还输出不同高度层内空间碎片状态数据。
优选的,还包括:
S7,空间碎片演化结果处理子系统根据空间碎片状态推演子系统推演到的空间碎片未来运动状态,计算空间碎片状态量,进一步计算到空间碎片密度;根据空间碎片密度得到空间碎片沿轨道高度方向和纬度方向上的分布状态。
优选的,所述空间碎片密度通过以下步骤计算得到:
步骤S7.1,确定空间碎片在高度区间[r,r′]内的平均空间密度分布其中,r和r′分别为高度区间的最小地心距和最大地心距;包括以下步骤:
步骤S7.1.1,当目标轨道为椭圆轨道,近地点地心距和远地点地心距分别为ra和rp,当r′≤ra且r≥rp时,通过下式计算:
其中:a为目标轨道半长轴;
步骤S7.1.2,当目标轨道为椭圆轨道,且r′>ra或r<rp时,则令r′=ra或r=rp,带入(28)进行计算,得到
步骤S7.1.3,当目标轨道为圆轨或近圆轨道,即rp→ra∈[r,r′]时,则令r′=ra或r=rp,公式(28)变为公式(29),通过公式(29)计算
步骤S7.1.4,其他情况,
步骤S7.2,确定空间碎片在纬度区间[β,β′]内的平均空间密度分布函数其中,β和β′分别是纬度区间的最小值和最大值;
具体的,若β′>0,则直接执行步骤S7.2.1;
若β′≤0,则令β′=-β,β=-β′,均转换为关于原点对称的正值区间后,再执行步骤S7.2.1;
步骤S7.2.1,碎片轨道倾角记为i,当β′≤i且β≥-i时
步骤S7.2.2,当β≤-i且β′≤i时
步骤S7.2.3,当β′≥i且-i≤β≤i时
步骤S7.2.4,当β′≥i且β≤-i时
步骤S7.2.5,当β>i时,
本发明还提供一种大规模空间碎片分布演化的数值计算系统,包括:
空间碎片状态推演子系统,用于获取空间碎片的初始运动状态数据以及空间碎片受到的空间摄动力,并将所述初始运动状态数据以及所述空间摄动力输入到数值积分模型中;所述数值积分模型按照给定的时间步长进行积分推演,得到空间碎片未来运动状态数据;然后,采用所述空间碎片未来运动状态数据更新空间碎片当前运动状态数据,得到空间碎片新运动状态数据;
空间碎片碰撞/爆炸解体产生新碎片子系统,用于将所述空间碎片新运动状态数据输入到碎片解体模型中;所述碎片解体模型进行解体计算,得到空间碎片解体后产生新空间解体碎片的属性数据;并将所述新空间解体碎片的属性数据反馈传输到所述空间碎片状态推演子系统;
空间碎片运动状态三维显示子系统,用于将所述空间碎片状态推演子系统推演到的空间碎片未来运动状态数据进行三维显示,即:动态的显示所有空间碎片在空间中的位置分布状态;
空间碎片运动状态输出子系统,用于将所述空间碎片状态推演子系统推演到的空间碎片未来运动状态输出;其中,所述空间碎片未来运动状态包括:空间碎片的位置、速度矢量和面值比;此外,空间碎片运动状态输出子系统还输出不同高度层内空间碎片状态数据;
空间碎片演化结果处理子系统,用于根据空间碎片状态推演子系统推演到的空间碎片未来运动状态,计算空间碎片状态量,进一步计算到空间碎片密度;根据空间碎片密度得到空间碎片沿轨道高度方向和纬度方向上的分布状态。
本发明提供的大规模空间碎片分布演化的数值计算方法及系统具有以下优点:
本发明提供的大规模空间碎片分布演化的数值计算方法及系统,能够对大规模空间碎片分布状态进行长期演化计算,能够动态三维演示碎片分布状态,提供了碎片分布状态数据的输出和处理分析功能,可以为航天器在轨运行空间碎片碰撞规避、空间碎片减缓措施制定以及确保空间资源可持续利用的重要的保障。
附图说明
图1为本发明提供的大规模空间碎片分布演化的数值计算方法的流程示意图;
图2为本发明提供的空间碎片状态推演子系统的结构示意图;
图3是大规模空间碎片分布演化数值计算系统运行时,空间碎片初始分布状态的三维显示图;
其中:在图3中,积分步数为0步;当前空间碎片数量为17000个;已发生碰撞次数为0次;
图4是大规模空间碎片分布演化数值计算系统运行过程中,产生新的解体碎片的分布状态示意图。
具体实施方式
为了使本发明所解决的技术问题、技术方案及有益效果更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
本发明提供一种大规模空间碎片分布演化的数值计算方法,应用于由5个子系统构成的大规模空间碎片分布演化的数值计算系统中,大规模空间碎片分布演化的数值计算系统包括:空间碎片状态推演子系统、空间碎片碰撞/爆炸解体产生新碎片子系统、空间碎片运动状态三维显示子系统、空间碎片运动状态输出子系统和空间碎片演化结果处理子系统。
大规模空间碎片分布演化的数值计算方法包括以下步骤:
S1,空间碎片状态推演子系统获取空间碎片的初始运动状态数据以及空间碎片受到的空间摄动力,并将所述初始运动状态数据以及所述空间摄动力输入到数值积分模型中;所述数值积分模型按照给定的时间步长进行积分推演,得到空间碎片未来运动状态数据;然后,采用所述空间碎片未来运动状态数据更新空间碎片当前运动状态数据。
本步骤中,所述空间摄动力包括大气阻力摄动力、地球非球形引力摄动力、太阳光压摄动力以及第三体引力摄动力;所述数值积分模型为4阶阿达姆斯(Adams-Bashforth)预估矫正积分模型。
空间碎片状态推演子系统的结构如图2所示,具体实现过程如下:
步骤S1.1,计算空间碎片受到的大气阻力摄动力。
大气阻力作为唯一清除在轨空间碎片的因素,其作用随轨道高度变化较为明显。计算表明,250km的典型空间碎片,在大气阻力的作用下,2个月内会再入到大气层内;而600km高度的空间碎片,再入到大气层的时间达到15年,高于850km的空间碎片,再入到大气层的时间将需要数个世纪之久,而地球同步轨道处的空间碎片在通常情况下永远不会再入到大气层中,除非空间环境发生了大规模变化。
在轨空间碎片受到的大气阻力随着大气状态的不同而变化,阻力加速度计算式为
其中:是大气阻力摄动力;CD是空间碎片的阻力系数;A是空间碎片的迎风截面积;md是空间碎片的质量;ρ是空间碎片所在位置的大气密度;是空间碎片与当地大气的相对速度;
假设大气相对地球静止,即大气相对于地球的速度其中为地球自转角速度,为地心矢径;若空间碎片的速度为则相对速度阻力加速度的计算式可进一步改写为
步骤S1.2,计算空间碎片受到的地球非球形引力摄动力。
在地固坐标系中,地球引力位函数作为拉普拉斯方程的解,其非球形部分U为:
其中:和为勒让德多项式相关项,定义为:
其中:
G为引力常数,Me是地球质量,r为地心距;
U表示地球引力场位函数;
λ和分别表示单位质点在地固坐标系中的地心经和纬度;
ae表示地球平均半径;
为规格化的缔合勒让德多项式;
和为规格化的地球引力位系数;
n和m为多项式的阶和次,N为取的最高阶数。
步骤S2.3,计算空间碎片受到的太阳光压摄动力。
太阳光压摄动力加速度可以表示为
其中:为太阳光压摄动力加速度;pSP=4.5605×10-6N/m2为地球附近太阳光压常数;CR为太阳光压系数;ASR是与太阳垂直的目标横截面积;AU=1.49597870×1011m是IAU1976天文常数系统给出的天文单位;为太阳位置矢量。
步骤S2.4,计算空间碎片受到的第三体引力摄动力。
第三体引力摄动力加速度为
其中:为第三体引力摄动力加速度;S和L分别代表太阳和月球;Mj为太阳或月球质量;代表太阳或月球的地心矢径。太阳和月球的位置可以通过平根公式计算,也可以采用JPL星历确定。
步骤S2.5,利用4阶阿达姆斯预估矫正方法对空间碎片的运动状态进行积分更新。该积分方法利用4阶龙格-库塔方法进行参数初始化,充分利用历史计算结果,每步积分仅计算一次摄动力,减少积分计算量;采用4阶积分方法,比高阶积分具有跟高的计算稳定性。
4阶龙格-库塔参数初始化方法为:
针对初值问题:
其中:t是时刻;y是状态量,如空间碎片的位置和速度矢量的坐标;是状态变量y的导数;f(t,y)是约束函数,如空间碎片受到大气阻力加速度、地球非球形引力摄动力加速度以及太阳光压摄动力加速度等;tn表示离散化后的不同时刻;yn是第n个时刻状态变量的取值。
龙格-库塔积分公式为:
其中:h=tn+1-tn,为积分步长;k为积分式所取的阶数,当取值为4的时候即为4阶龙格-库塔初始化方法;Ci,ai,bij均为已知的常数项。
设已知4个时刻的函数值分别为fi-3,fi-2,fi-1,fi,则i+1时刻y的近似估计值为
S2,判断是否达到演化时间,如果未达到,将状态更新后的空间碎片运动状态和属性数据传输给空间碎片碰撞/爆炸解体产生新碎片子系统;
S3,S3,所述空间碎片碰撞/爆炸解体产生新碎片子系统将所述状态更新后的空间碎片运动状态和属性数据输入到碎片解体模型中;所述碎片解体模型进行解体计算,得到空间碎片解体后产生新空间解体碎片的属性数据;并将所述新空间解体碎片的属性数据反馈传输到所述空间碎片状态推演子系统;
其中,所述空间碎片碰撞/爆炸解体产生新碎片子系统,具体用于:
首先,根据能量定律确定空间碎片解体后产生新空间解体碎片的数量N和新空间解体碎片特征尺寸l的分布;其中,碎片的特征尺寸l=(lx+ly+lz)/3,即为碎片3个轴向尺寸的平均。
然后,以新空间解体碎片特征尺寸l作为独立变量,利用概率分布模型确定新空间解体碎片的面质比A/md、新空间解体碎片的迎风截面积A、新空间解体碎片相对于解体前空间碎片的速度增量Δv;
最后,根据A/md、A确定碎片质量md。
更具体的,所述空间碎片碰撞/爆炸解体产生新碎片子系统具体用于:
S3.1,确定解体碎片数量:
其中,解体碎片区分为爆炸解体碎片和碰撞解体碎片;
1)对于爆炸解体碎片,根据能量守恒定律,解体碎片特征尺寸l大于被研究解体碎片最小特征尺寸lc的碎片数量Nf满足函数关系:
其中:为解体碎片的最小特征尺寸;
cs是修正系数,与爆炸类型有关,对于历史爆炸事件,修正系数cs的表达式:
其中:修正系数cs与目标探测直径门限dSSN直接相关,dSSN与目标高度H的关系式如下:
其中:H为爆炸解体碎片的高度;
2)对于碰撞解体碎片,解体碎片的数量和分布与相互碰撞的能量有关。碰撞分为灾难性碰撞和非灾难性碰撞,两种碰撞条件下解体碎片的数量和分布规律满足不同的表达式。
一次碰撞是否为灾难性碰撞可通过下式判断:
其中:Ep为等效碰撞能量,为灾难性碰撞临界能量。mp为相互碰撞中较小质量;mt为相互碰撞中较大质量;vr是碰撞时两碎片的相对速度,即碰撞速度;当上式成立时,碰撞为灾难性碰撞,反之则为非灾难性碰撞。
因此,给出下面的碰撞解体碎片数量满足函数关系:
式中:为两碰撞目标的等效质量:
其中:
S3.2,确定解体碎片尺寸分布律:
爆炸解体碎片尺寸分布律F:
碰撞解体碎片尺寸分布律F:
计算时,首先确定要研究解体碎片的最小尺寸lc,然后根据关系式(11)和(15)计算爆炸或解体产生大于lc尺寸碎片的总数目,再根据分布律(17)或(18)确定这些碎片具有的尺寸。利用分布律确定解体碎片尺寸时,需要对随机变量抽样,可采用反函数的方法进行随机变量的抽样。
S3.3,确定解体碎片面质比参数:
得到解体碎片的特征尺寸以后,碎片的面质比参数满足确定的分布函数,利用随机变量采样法可以得到每一个碎片具有的面质比。
本步骤中,爆炸解体碎片和碰撞解体碎片的碎片面质比参数满足同样的分布函数;
(1)对于尺寸大于11cm解体碎片,定义解体碎片面质比的对数值γ=log(A/md),解体碎片特征长度的对数值θ=log(l),解体碎片的面质比分布由如下二元正态分布决定:
p(γ,θ)=ε(θ)p1(γ)+(1-ε(θ))p2(γ) (19)
其中:p1(γ)和p2(γ)均是正态分布概率密度函数:
其中:权重参数ε(θ)∈[0,1],是θ的饱和函数;正态分布概率密度函数的均值μi和正态分布概率密度函数的方差σi均是θ的函数;
对于航天器解体碎片,其面质比分布函数中的参数由下式决定:
对于火箭上面级解体碎片,其面质比分布函数的参数由下式决定
(2)对于尺寸小于8cm的航天器解体碎片和尺寸小于1.7cm的火箭上面级解体碎片,碎片的面质比满足正态分布律p1(γ),即权重函数ε=1,
p(γ,θ)=p1(γ) (23)
此时,分布函数参数由下式统一决定:
(3)对于尺寸介于8cm~11cm之间的航天器解体碎片和尺寸介于1.7cm~11cm之间的火箭上面级碎片,碎片的面质比通过随机采样的方法确定,即:首先生成服从均匀分布且取值区间在[0.0,1.0]内的随机变量ζ,然后将ζ与利用(25)计算得到的进行对比,若则利用公式(19)计算面质比参数,反之则利用公式(23)计算面质比参数;
步骤3.4,确定解体碎片的速度增量:
解体碎片相对于解体前目标的速度增量满足正态分布规律,定义v=lg(Δv),γ=lg(A/md),则v的分布函数p(ν)满足
其中:v为解体碎片的速度;Δv为解体碎片的速度增量;
解体碎片的均值μν和解体碎片的方差σν通过下式计算:
S4,所述空间碎片状态推演子系统将所述新空间解体碎片的属性数据以及已有的空间碎片的运动状态数据输入到数值积分模型中,再次按照给定的时间步长进行积分推演,得到空间碎片未来运动状态数据,如此不断循环S2-S4;直到达到演化时间,结束流程。
还包括:
S5,空间碎片运动状态三维显示子系统将所述空间碎片状态推演子系统推演到的空间碎片未来运动状态数据进行三维显示,即:利用DirectX等三维显示技术,动态的显示所有空间碎片在空间中的位置分布状态;
具体包括:
步骤S5.1,接收空间碎片状态推演子系统推演到的空间碎片未来运动状态数据;
步骤S5.2,调用三维显示接口函数(API),将空间碎片的位置和形状在三维模型中进行渲染;
步骤S5.3,驱动三维模型,对空间碎片运动状态进行动态更新。
空间碎片的分布状态如图3和图4所示。
还包括:
S6,在固定时刻或者设定的时刻,空间碎片运动状态输出子系统将所述空间碎片状态推演子系统推演到的空间碎片未来运动状态输出;其中,所述空间碎片未来运动状态包括:空间碎片的位置、速度矢量和面值比、质量、来源属性参数;此外,空间碎片运动状态输出子系统还输出不同高度层内空间碎片状态数据。
还包括:
S7,空间碎片演化结果处理子系统根据空间碎片状态推演子系统推演到的空间碎片未来运动状态,计算空间碎片状态量,进一步计算到空间碎片密度;根据空间碎片密度得到空间碎片沿轨道高度方向和纬度方向上的分布状态。
所述空间碎片密度通过以下步骤计算得到:
步骤S7.1,确定空间碎片在高度区间[r,r′]内的平均空间密度分布其中,r和r′分别为高度区间的最小地心距和最大地心距;包括以下步骤:
步骤S7.1.1,当目标轨道为椭圆轨道,近地点地心距和远地点地心距分别为ra和rp,当r′≤ra且r≥rp时,通过下式计算:
其中:a为目标轨道半长轴;
步骤S7.1.2,当目标轨道为椭圆轨道,且r′>ra或r<rp时,则令r′=ra或r=rp,带入(28)进行计算,得到
步骤S7.1.3,当目标轨道为圆轨或近圆轨道,即rp→ra∈[r,r′]时,则令r′=ra或r=rp,公式(28)变为公式(29),通过公式(29)计算
步骤S7.1.4,其他情况,
步骤S7.2,确定空间碎片在纬度区间[β,β′]内的平均空间密度分布函数其中,β和β′分别是纬度区间的最小值和最大值;
具体的,若β′>0,则直接执行步骤S7.2.1;
若β′≤0,则令β′=-β,β=-β′,均转换为关于原点对称的正值区间后,再执行步骤S7.2.1;
步骤S7.2.1,碎片轨道倾角记为i,当β′≤i且β≥-i时
步骤S7.2.2,当β≤-i且β′≤i时
步骤S7.2.3,当β′≥i且-i≤β≤i时
步骤S7.2.4,当β′≥i且β≤-i时
步骤S7.2.5,当β>i时,
本发明还提供一种大规模空间碎片分布演化的数值计算系统,包括空间碎片状态推演子系统、空间碎片碰撞/爆炸解体产生新碎片子系统、空间碎片运动状态三维显示子系统、空间碎片运动状态输出子系统和空间碎片演化结果处理子系统。下在对这5个子系统详细介绍:
(1)空间碎片状态推演子系统
空间碎片状态推演子系统,用于获取空间碎片的初始运动状态数据以及空间碎片受到的空间摄动力,并将所述初始运动状态数据以及所述空间摄动力输入到数值积分模型中;所述数值积分模型按照给定的时间步长进行积分推演,得到空间碎片未来运动状态数据;然后,采用所述空间碎片未来运动状态数据更新空间碎片当前运动状态数据,得到空间碎片新运动状态数据。
也就是说,空间碎片状态推演子系统以空间碎片的初始运动状态为输入,综合考虑空间碎片受到的大气阻力、非球形引力、太阳光压、三体引力等摄动力的作用,采用高效稳定的积分方法,积分推演碎片的未来的运动状态。
其中,在计算大气阻力的过程中,要首先确定高层大气密度。太阳辐射作用使大气密度呈现周日变化、季节变化等周期性变化规律;太阳活动会导致大气密度发生剧烈变化,严重影响大气阻力对目标的衰减作用。计算系统选取了Harris-Priester大气密度模型,该模型能够反映太阳活动对大气密度的影响,模型复杂性低,对计算能力要求不高,适宜进行大规模的长期演化计算。
计算空间碎片受到的非球形摄动时,计算系统对其量级进行了分析,并进行了合理的取舍。地球扁率J2项是主要的摄动项,导致空间碎片半长轴和升交点赤经周期性变化,影响着物体沿经度方向上的分布规律;在扁率J2项和J2,2项影响下,GEO轨道带物体在地球赤道短轴上空,即东经75°和西经105°附近,呈现聚集分布规律;J3项和J4项的量级与J2,2项量级相当,均为10-6,计算空间碎片受到的非球形摄动时也加以考虑。
计算系统中的积分方法是经过优选确定的。对空间碎片的运动进行动力学积分时,要求积分算法精度高且具有有一定稳定性,避免长期计算时结果发散;由于空间碎片数量规模大,需要开发高效的并行计算方法。系统采用4阶阿达姆斯预估矫正方法(Adams-Bashforth)进行积分,该积分方法利用4阶龙格-库塔方法进行参数初始化,充分利用历史计算结果,每步积分仅计算一次摄动力,减少积分计算量;采用4阶积分方法,比高阶积分具有跟高的计算稳定性;模型基于MPI的并行计算框架,利用多进程和多核处理技术,设计了高效的并行演化系统。
(2)空间碎片碰撞/爆炸解体产生新碎片子系统
空间碎片碰撞/爆炸解体产生新碎片子系统,用于将所述空间碎片新运动状态数据输入到碎片解体模型中;所述碎片解体模型进行解体计算,得到空间碎片解体后产生新空间解体碎片的属性数据;并将所述新空间解体碎片的属性数据反馈传输到所述空间碎片状态推演子系统。
空间碎片碰撞/爆炸解体产生新碎片子系统的功能包括两个方面:首先确定空间碎片之间发生碰撞的概率,然后根据碎片解体模型生成大量解体碎片。然后,根据能量定律确定解体后产生碎片的数量N和碎片特征尺寸l的分布;然后以特征尺寸l作为独立变量,利用概率分布模型确定碎片的面质比A/md、有效横截面积A和速度增量Δv;最后根据A/md、A确定碎片质量md。
(3)空间碎片运动状态三维显示子系统
空间碎片运动状态三维显示子系统,用于利用DirectX等三维显示技术,将所述空间碎片状态推演子系统推演到的空间碎片未来运动状态数据进行三维显示,子系统能够根据物体的推演状态,动态的显示所有物体在空间中的位置分布状态。
(4)空间碎片运动状态输出子系统
空间碎片运动状态输出子系统,用于在固定时刻或者设定的时刻,将所述空间碎片状态推演子系统推演到的空间碎片未来运动状态输出,具体的,输出包块解体碎片在内的当前所有碎片的位置、速度矢量和物体的面值比等属性参数;此外,空间碎片运动状态输出子系统还输出不同高度层内空间碎片状态数据;
(5)空间碎片演化结果处理子系统
空间碎片演化结果处理子系统,用于根据空间碎片状态推演子系统推演到的空间碎片未来运动状态,分析碎片沿轨道高度方向、纬度方向等空间维度方向上的变化情况,分析碎片规模随时间的变化规律。计算空间碎片状态量,进一步计算到空间碎片密度;根据空间碎片密度得到空间碎片沿轨道高度方向和纬度方向上的分布状态。
本发明提供的大规模空间碎片分布演化数值计算方法及系统,能够根据碎片的初始运动状态,综合考虑大气阻力、地球非球形摄动力以及太阳/月球三体引力等复杂摄动力的作用,碎片相互碰撞产生新碎片的作用,对碎片的分布状态进行长期演化计算。系统能够三维动态的显示空间碎片在轨分布状态随时间的变化,能够对碎片分布状态的演化计算结果进行统计分析,得到空间碎片在不同轨道高度、不同纬度等方向上的分布规律,从而为航天器在轨运行碰撞规避、空间碎片减缓措施制定以及确保空间资源可持续利用的重要的技术保障。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视本发明的保护范围。
Claims (8)
1.一种大规模空间碎片分布演化的数值计算方法,其特征在于,包括以下步骤:
S1,空间碎片状态推演子系统获取空间碎片的初始运动状态数据以及空间碎片受到的空间摄动力,并将所述初始运动状态数据以及所述空间摄动力输入到数值积分模型中;所述数值积分模型按照给定的时间步长进行积分推演,得到空间碎片未来运动状态数据;然后,采用所述空间碎片未来运动状态数据更新空间碎片当前运动状态数据;
空间碎片状态推演子系统具体实现过程如下:
步骤S1.1,计算空间碎片受到的大气阻力摄动力;
在轨空间碎片受到的大气阻力随着大气状态的不同而变化,阻力加速度计算式为
其中:是大气阻力摄动力;CD是空间碎片的阻力系数;A是空间碎片的迎风截面积;md是空间碎片的质量;ρ是空间碎片所在位置的大气密度;是空间碎片与当地大气的相对速度;
假设大气相对地球静止,即大气相对于地球的速度其中为地球自转角速度,为地心矢径;若空间碎片的速度为则相对速度阻力加速度的计算式可进一步改写为
步骤S1.2,计算空间碎片受到的地球非球形引力摄动力;
在地固坐标系中,地球引力位函数作为拉普拉斯方程的解,其非球形部分U为:
<mrow>
<mi>U</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>2</mn>
</mrow>
<mi>N</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>n</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>C</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>n</mi>
<mi>m</mi>
</mrow>
</msub>
<msubsup>
<mover>
<mi>U</mi>
<mo>&OverBar;</mo>
</mover>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mo>+</mo>
<msub>
<mover>
<mi>S</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>n</mi>
<mi>m</mi>
</mrow>
</msub>
<msubsup>
<mover>
<mi>V</mi>
<mo>&OverBar;</mo>
</mover>
<mi>n</mi>
<mi>m</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:和为勒让德多项式相关项,定义为:
其中:
G为引力常数,Me是地球质量,r为地心距;
U表示地球引力场位函数;
λ和分别表示单位质点在地固坐标系中的地心经和纬度;
ae表示地球平均半径;
为规格化的缔合勒让德多项式;
和为规格化的地球引力位系数;
n和m为多项式的阶和次,N为取的最高阶数;
步骤S1.3,计算空间碎片受到的太阳光压摄动力;
太阳光压摄动力加速度可以表示为
其中:为太阳光压摄动力加速度;pSP=4.5605×10-6N/m2为地球附近太阳光压常数;CR为太阳光压系数;ASR是与太阳垂直的目标横截面积;AU=1.49597870×1011m是IAU1976天文常数系统给出的天文单位;为太阳位置矢量;
步骤S1.4,计算空间碎片受到的第三体引力摄动力;
第三体引力摄动力加速度为
其中:为第三体引力摄动力加速度;S和L分别代表太阳和月球;Mj为太阳或月球质量;代表太阳或月球的地心矢径;太阳和月球的位置可以通过平根公式计算,也可以采用JPL星历确定;
步骤S1.5,利用4阶阿达姆斯预估矫正方法对空间碎片的运动状态进行积分更新;该积分方法利用4阶龙格-库塔方法进行参数初始化,充分利用历史计算结果,每步积分仅计算一次摄动力,减少积分计算量;采用4阶积分方法,比高阶积分具有跟高的计算稳定性;
4阶龙格-库塔参数初始化方法为:
针对初值问题:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>y</mi>
<mi>n</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:t是时刻;y是状态量,如空间碎片的位置和速度矢量的坐标;是状态变量y的导数;f(t,y)是约束函数,包括空间碎片受到大气阻力加速度、地球非球形引力摄动力加速度以及太阳光压摄动力加速度;tn表示离散化后的不同时刻;yn是第n个时刻状态变量的取值;
龙格-库塔积分公式为:
<mrow>
<msub>
<mi>y</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>y</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<mi>h</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>k</mi>
</munderover>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mi>h</mi>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msub>
<mi>b</mi>
<mn>10</mn>
</msub>
<msub>
<mi>hf</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mi>k</mi>
</msub>
<mi>h</mi>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mi>k</mi>
</msub>
<mi>h</mi>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>k</mi>
</munderover>
<msub>
<mi>b</mi>
<mrow>
<mi>k</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:h=tn+1-tn,为积分步长;k为积分式所取的阶数,当取值为4的时候即为4阶龙格-库塔初始化方法;Ci,ai,bij均为已知的常数项;
设已知4个时刻的函数值分别为fi-3,fi-2,fi-1,fi,则i+1时刻y的近似估计值为
<mrow>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<mfrac>
<mi>h</mi>
<mn>24</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mn>9</mn>
<msub>
<mi>f</mi>
<mrow>
<mo>-</mo>
<mn>3</mn>
</mrow>
</msub>
<mo>+</mo>
<mn>37</mn>
<msub>
<mi>f</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<mn>59</mn>
<msub>
<mi>f</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<mn>55</mn>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
S2,判断是否达到演化时间,如果未达到,将状态更新后的空间碎片运动状态和属性数据传输给空间碎片碰撞/爆炸解体产生新碎片子系统;
S3,所述空间碎片碰撞/爆炸解体产生新碎片子系统将所述状态更新后的空间碎片运动状态和属性数据输入到碎片解体模型中;所述碎片解体模型进行解体计算,得到空间碎片解体后产生新空间解体碎片的属性数据;并将所述新空间解体碎片的属性数据反馈传输到所述空间碎片状态推演子系统;
S4,所述空间碎片状态推演子系统将所述新空间解体碎片的属性数据以及已有的空间碎片的运动状态数据输入到数值积分模型中,再次按照给定的时间步长进行积分推演,得到空间碎片未来运动状态数据,如此不断循环S2-S4;直到达到演化时间,结束流程。
2.根据权利要求1所述的大规模空间碎片分布演化的数值计算方法,其特征在于,S1中,所述空间摄动力包括大气阻力摄动力、地球非球形引力摄动力、太阳光压摄动力以及第三体引力摄动力;所述数值积分模型为4阶阿达姆斯预估矫正积分模型。
3.根据权利要求1所述的大规模空间碎片分布演化的数值计算方法,其特征在于,S3中,所述空间碎片碰撞/爆炸解体产生新碎片子系统具体用于:
首先,根据能量定律确定空间碎片解体后产生新空间解体碎片的数量N和新空间解体碎片特征尺寸l的分布;
然后,以新空间解体碎片特征尺寸l作为独立变量,利用概率分布模型确定新空间解体碎片的面质比A/md、新空间解体碎片的迎风截面积A、新空间解体碎片相对于解体前空间碎片的速度增量△v;
最后,根据A/md、A确定碎片质量md。
4.根据权利要求3所述的大规模空间碎片分布演化的数值计算方法,其特征在于,S3中,所述空间碎片碰撞/爆炸解体产生新碎片子系统具体用于:
S3.1,确定解体碎片数量:
其中,解体碎片区分为爆炸解体碎片和碰撞解体碎片;
1)对于爆炸解体碎片,根据能量守恒定律,解体碎片特征尺寸l大于被研究解体碎片最小特征尺寸lc的碎片数量Nf满足函数关系:
<mrow>
<msub>
<mi>N</mi>
<mi>f</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>&GreaterEqual;</mo>
<msub>
<mi>l</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>6</mn>
<msub>
<mi>c</mi>
<mi>s</mi>
</msub>
<msubsup>
<mover>
<mi>l</mi>
<mo>^</mo>
</mover>
<mi>c</mi>
<mrow>
<mo>-</mo>
<mn>1.6</mn>
</mrow>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:为解体碎片的最小特征尺寸;
cs是修正系数,与爆炸类型有关,对于历史爆炸事件,修正系数cs的表达式:
<mrow>
<msub>
<mi>c</mi>
<mi>s</mi>
</msub>
<mo>=</mo>
<msup>
<mn>10</mn>
<mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mi>exp</mi>
<mrow>
<mo>{</mo>
<mrow>
<mo>-</mo>
<mn>2.464</mn>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mi>lg</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>d</mi>
<mrow>
<mi>S</mi>
<mi>S</mi>
<mi>N</mi>
</mrow>
</msub>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>m</mi>
<mi>d</mi>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>1.22</mn>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mo>}</mo>
</mrow>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:
其中:H为爆炸解体碎片的高度;
2)对于碰撞解体碎片,碰撞解体碎片数量满足函数关系:
<mrow>
<msub>
<mi>N</mi>
<mi>f</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>&GreaterEqual;</mo>
<msub>
<mi>l</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0.1</mn>
<msup>
<mover>
<mi>m</mi>
<mo>^</mo>
</mover>
<mn>0.75</mn>
</msup>
<msubsup>
<mover>
<mi>l</mi>
<mo>^</mo>
</mover>
<mi>c</mi>
<mrow>
<mo>-</mo>
<mn>1.71</mn>
</mrow>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:为两碰撞目标的等效质量:
其中:
<mrow>
<mover>
<mi>m</mi>
<mo>^</mo>
</mover>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>m</mi>
<mi>t</mi>
</msub>
<mo>+</mo>
<msub>
<mi>m</mi>
<mi>p</mi>
</msub>
</mrow>
<mrow>
<mo>&lsqb;</mo>
<mi>k</mi>
<mi>g</mi>
<mo>&rsqb;</mo>
</mrow>
</mfrac>
</mtd>
<mtd>
<mrow>
<msub>
<mi>E</mi>
<mi>p</mi>
</msub>
<mo>&GreaterEqual;</mo>
<msubsup>
<mi>E</mi>
<mi>p</mi>
<mo>*</mo>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mrow>
<msub>
<mi>m</mi>
<mi>p</mi>
</msub>
<msub>
<mi>v</mi>
<mi>r</mi>
</msub>
</mrow>
<mrow>
<mn>1000</mn>
<mo>&lsqb;</mo>
<mi>k</mi>
<mi>g</mi>
<mo>&CenterDot;</mo>
<mi>m</mi>
<mo>/</mo>
<mi>s</mi>
<mo>&rsqb;</mo>
</mrow>
</mfrac>
</mtd>
<mtd>
<mrow>
<msub>
<mi>E</mi>
<mi>p</mi>
</msub>
<mo><</mo>
<msubsup>
<mi>E</mi>
<mi>p</mi>
<mo>*</mo>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:Ep为等效碰撞能量,为灾难性碰撞临界能量;mp为相互碰撞中较小质量;mt为相互碰撞中较大质量;vr是碰撞时两碎片的相对速度,即碰撞速度;
S3.2,确定解体碎片尺寸分布律:
爆炸解体碎片尺寸分布律F:
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>{</mo>
<mtable>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msubsup>
<mi>l</mi>
<mi>c</mi>
<mn>1.6</mn>
</msubsup>
<msup>
<mi>l</mi>
<mrow>
<mo>-</mo>
<mn>1.6</mn>
</mrow>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>,</mo>
<mi>l</mi>
<mo>&GreaterEqual;</mo>
<msub>
<mi>l</mi>
<mi>c</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>,</mo>
<mi>l</mi>
<mo><</mo>
<msub>
<mi>l</mi>
<mi>c</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
碰撞解体碎片尺寸分布律F:
<mrow>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>{</mo>
<mtable>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msubsup>
<mi>l</mi>
<mi>c</mi>
<mn>1.71</mn>
</msubsup>
<msup>
<mi>l</mi>
<mrow>
<mo>-</mo>
<mn>1.71</mn>
</mrow>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>,</mo>
<mi>l</mi>
<mo>&GreaterEqual;</mo>
<msub>
<mi>l</mi>
<mi>c</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>,</mo>
<mi>l</mi>
<mo><</mo>
<msub>
<mi>l</mi>
<mi>c</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>18</mn>
<mo>)</mo>
</mrow>
</mrow>
S3.3,确定解体碎片面质比参数:
本步骤中,爆炸解体碎片和碰撞解体碎片的碎片面质比参数满足同样的分布函数;
(1)对于尺寸大于11cm解体碎片,定义解体碎片面质比的对数值γ=log(A/md),解体碎片特征长度的对数值θ=log(l),解体碎片的面质比分布由如下二元正态分布决定:
p(γ,θ)=ε(θ)p1(γ)+(1-ε(θ))p2(γ) (19)
其中:p1(γ)和p2(γ)均是正态分布概率密度函数:
<mrow>
<msub>
<mi>p</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&gamma;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>&sigma;</mi>
<mi>k</mi>
</msub>
<msqrt>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>&gamma;</mi>
<mo>-</mo>
<msub>
<mi>&mu;</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mn>2</mn>
<msubsup>
<mi>&sigma;</mi>
<mi>k</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>20</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:权重参数ε(θ)∈[0,1],是θ的饱和函数;正态分布概率密度函数的均值μk和正态分布概率密度函数的方差σk均是θ的函数;
对于航天器解体碎片,其面质比分布函数中的参数由下式决定:
<mrow>
<mi>&epsiv;</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>&theta;</mi>
<mo>&le;</mo>
<mo>-</mo>
<mn>1.95</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.3</mn>
<mo>+</mo>
<mn>0.4</mn>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>+</mo>
<mn>1.2</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>1.95</mn>
<mo><</mo>
<mi>&theta;</mi>
<mo><</mo>
<mn>0.55</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow>
<mi>&theta;</mi>
<mo>&GreaterEqual;</mo>
<mn>0.55</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<msub>
<mi>&mu;</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.6</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>&theta;</mi>
<mo>&le;</mo>
<mo>-</mo>
<mn>1.1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.6</mn>
<mo>-</mo>
<mn>0.318</mn>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>+</mo>
<mn>1.1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>1.1</mn>
<mo><</mo>
<mi>&theta;</mi>
<mo><</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.95</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>&theta;</mi>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<msub>
<mi>&sigma;</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>0.1</mn>
</mtd>
<mtd>
<mrow>
<mi>&theta;</mi>
<mo>&le;</mo>
<mo>-</mo>
<mn>1.3</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.1</mn>
<mo>+</mo>
<mn>0.2</mn>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>+</mo>
<mn>1.3</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>1.3</mn>
<mo><</mo>
<mi>&theta;</mi>
<mo><</mo>
<mo>-</mo>
<mn>0.3</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0.3</mn>
</mtd>
<mtd>
<mrow>
<mi>&theta;</mi>
<mo>&GreaterEqual;</mo>
<mo>-</mo>
<mn>0.3</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<msub>
<mi>&mu;</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>1.2</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>&theta;</mi>
<mo>&le;</mo>
<mo>-</mo>
<mn>0.7</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>1.2</mn>
<mo>-</mo>
<mn>1.333</mn>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>+</mo>
<mn>0.7</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.7</mn>
<mo><</mo>
<mi>&theta;</mi>
<mo><</mo>
<mo>-</mo>
<mn>0.1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>2.0</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>&theta;</mi>
<mo>&GreaterEqual;</mo>
<mo>-</mo>
<mn>0.1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<msub>
<mi>&sigma;</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>0.5</mn>
</mtd>
<mtd>
<mrow>
<mi>&theta;</mi>
<mo>&le;</mo>
<mo>-</mo>
<mn>0.5</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.5</mn>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>+</mo>
<mn>0.5</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.5</mn>
<mo><</mo>
<mi>&theta;</mi>
<mo><</mo>
<mo>-</mo>
<mn>0.3</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0.3</mn>
</mtd>
<mtd>
<mrow>
<mi>&theta;</mi>
<mo>&GreaterEqual;</mo>
<mo>-</mo>
<mn>0.3</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>21</mn>
<mo>)</mo>
</mrow>
</mrow>
对于火箭上面级解体碎片,其面质比分布函数的参数由下式决定
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&epsiv;</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow>
<mi>&theta;</mi>
<mo>&le;</mo>
<mo>-</mo>
<mn>1.4</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mn>0.3571</mn>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>+</mo>
<mn>1.4</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>1.4</mn>
<mo><</mo>
<mi>&theta;</mi>
<mo><</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0.5</mn>
</mtd>
<mtd>
<mrow>
<mi>&theta;</mi>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&mu;</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>0.45</mn>
</mtd>
<mtd>
<mrow>
<mi>&theta;</mi>
<mo>&le;</mo>
<mo>-</mo>
<mn>0.5</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.45</mn>
<mo>-</mo>
<mn>0.9</mn>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>+</mo>
<mn>0.5</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.5</mn>
<mo><</mo>
<mi>&theta;</mi>
<mo><</mo>
<mo>-</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.9</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>&theta;</mi>
<mo>&GreaterEqual;</mo>
<mo>-</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&sigma;</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0.55</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&mu;</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<mn>0.9</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&sigma;</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>0.28</mn>
</mtd>
<mtd>
<mrow>
<mi>&theta;</mi>
<mo>&le;</mo>
<mo>-</mo>
<mn>1.0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.28</mn>
<mo>-</mo>
<mn>0.1636</mn>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>1.0</mn>
<mo><</mo>
<mi>&theta;</mi>
<mo><</mo>
<mn>0.1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0.1</mn>
</mtd>
<mtd>
<mrow>
<mi>&theta;</mi>
<mo>&GreaterEqual;</mo>
<mn>0.1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>22</mn>
<mo>)</mo>
</mrow>
</mrow>
(2)对于尺寸小于8cm的航天器解体碎片和尺寸小于1.7cm的火箭上面级解体碎片,碎片的面质比满足正态分布律p1(γ),即权重函数ε=1,
p(γ,θ)=p1(γ) (23)
此时,分布函数参数由下式统一决定:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&mu;</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.3</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>&theta;</mi>
<mo>&le;</mo>
<mo>-</mo>
<mn>1.75</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>0.3</mn>
<mo>-</mo>
<mn>1.4</mn>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>+</mo>
<mn>1.75</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>1.75</mn>
<mo><</mo>
<mi>&theta;</mi>
<mo><</mo>
<mo>-</mo>
<mn>1.25</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>1.0</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>&theta;</mi>
<mo>&GreaterEqual;</mo>
<mo>-</mo>
<mn>1.25</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&sigma;</mi>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>0.2</mn>
</mtd>
<mtd>
<mrow>
<mi>&theta;</mi>
<mo>&le;</mo>
<mo>-</mo>
<mn>3.5</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.2</mn>
<mo>+</mo>
<mn>0.1333</mn>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mo>+</mo>
<mn>3.5</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>&theta;</mi>
<mo>></mo>
<mo>-</mo>
<mn>3.5</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>24</mn>
<mo>)</mo>
</mrow>
</mrow>
(3)对于尺寸介于8cm~11cm之间的航天器解体碎片和尺寸介于1.7cm~11cm之间的火箭上面级碎片,碎片的面质比通过随机采样的方法确定,即:首先生成服从均匀分布且取值区间在[0.0,1.0]内的随机变量ζ,然后将ζ与利用(25)计算得到的进行对比,若则利用公式(19)计算面质比参数,反之则利用公式(23)计算面质比参数;
步骤3.4,确定解体碎片的速度增量:
解体碎片相对于解体前目标的速度增量满足正态分布规律,定义v=lg(△v),γ=lg(A/md),则v的分布函数p(ν)满足
<mrow>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>&sigma;</mi>
<mi>v</mi>
</msub>
<msqrt>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>-</mo>
<msub>
<mi>&mu;</mi>
<mi>v</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mn>2</mn>
<msubsup>
<mi>&sigma;</mi>
<mi>v</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>26</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:v为解体碎片的速度;为解体碎片的速度增量;
解体碎片的均值μν和解体碎片的方差σν通过下式计算:
σv=0.4。
5.根据权利要求1所述的大规模空间碎片分布演化的数值计算方法,其特征在于,还包括:
S5,空间碎片运动状态三维显示子系统将所述空间碎片状态推演子系统推演到的空间碎片未来运动状态数据进行三维显示,即:动态的显示所有空间碎片在空间中的位置分布状态;具体包括:
步骤S5.1,接收空间碎片状态推演子系统推演到的空间碎片未来运动状态数据;
步骤S5.2,调用三维显示接口函数,将空间碎片的位置和形状在三维模型中进行渲染;
步骤S5.3,驱动三维模型,对空间碎片运动状态进行动态更新。
6.根据权利要求5所述的大规模空间碎片分布演化的数值计算方法,其特征在于,还包括:
S6,在固定时刻或者设定的时刻,空间碎片运动状态输出子系统将所述空间碎片状态推演子系统推演到的空间碎片未来运动状态输出;其中,所述空间碎片未来运动状态包括:空间碎片的位置、速度矢量和面值比、质量、来源属性参数;此外,空间碎片运动状态输出子系统还输出不同高度层内空间碎片状态数据。
7.根据权利要求6所述的大规模空间碎片分布演化的数值计算方法,其特征在于,还包括:
S7,空间碎片演化结果处理子系统根据空间碎片状态推演子系统推演到的空间碎片未来运动状态,计算空间碎片状态量,进一步计算到空间碎片密度;根据空间碎片密度得到空间碎片沿轨道高度方向和纬度方向上的分布状态。
8.根据权利要求7所述的大规模空间碎片分布演化的数值计算方法,其特征在于,所述空间碎片密度通过以下步骤计算得到:
步骤S7.1,确定空间碎片在高度区间[r,r′]内的平均空间密度分布其中,r和r′分别为高度区间的最小地心距和最大地心距;包括以下步骤:
步骤S7.1.1,当目标轨道为椭圆轨道,近地点地心距和远地点地心距分别为ra和rp,当r′≤ra且r≥rp时,通过下式计算:
<mrow>
<mover>
<mi>s</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<msup>
<mi>r</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mn>3</mn>
<mi>a</mi>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>a</mi>
<mi>r</mi>
<mi>c</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mo>(</mo>
<mn>2</mn>
<mfrac>
<mrow>
<msup>
<mi>r</mi>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<mi>a</mi>
</mrow>
<mrow>
<msub>
<mi>r</mi>
<mi>a</mi>
</msub>
<mo>-</mo>
<msub>
<mi>r</mi>
<mi>p</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>-</mo>
<mi>arcsin</mi>
<mo>(</mo>
<mn>2</mn>
<mfrac>
<mrow>
<mi>r</mi>
<mo>-</mo>
<mi>a</mi>
</mrow>
<mrow>
<msub>
<mi>r</mi>
<mi>a</mi>
</msub>
<mo>-</mo>
<msub>
<mi>r</mi>
<mi>p</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mn>3</mn>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msqrt>
<mrow>
<mo>(</mo>
<msup>
<mi>r</mi>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<msub>
<mi>r</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>r</mi>
<mi>a</mi>
</msub>
<mo>-</mo>
<msup>
<mi>r</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
</msqrt>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msqrt>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>-</mo>
<msub>
<mi>r</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
<mo>(</mo>
<msub>
<mi>r</mi>
<mi>a</mi>
</msub>
<mo>-</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
</msqrt>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<mn>4</mn>
<msup>
<mi>&pi;</mi>
<mn>2</mn>
</msup>
<mi>a</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>r</mi>
<mrow>
<mo>&prime;</mo>
<mn>3</mn>
</mrow>
</msup>
<mo>-</mo>
<msup>
<mi>r</mi>
<mn>3</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>28</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:a为目标轨道半长轴;
步骤S7.1.2,当目标轨道为椭圆轨道,且r′>ra或r<rp时,则令r′=ra或r=rp,带入(28)进行计算,得到
步骤S7.1.3,当目标轨道为圆轨或近圆轨道,即rp→ra∈[r,r′]时,则令r′=ra或r=rp,公式(28)变为公式(29),通过公式(29)计算
<mrow>
<mover>
<mi>s</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<msup>
<mi>r</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mn>3</mn>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mi>a</mi>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>8</mn>
<msup>
<mi>&pi;</mi>
<mn>2</mn>
</msup>
<mi>a</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>r</mi>
<mrow>
<mo>&prime;</mo>
<mn>3</mn>
</mrow>
</msup>
<mo>-</mo>
<msup>
<mi>r</mi>
<mn>3</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>29</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤S7.1.4,其他情况,
步骤S7.2,确定空间碎片在纬度区间[β,β′]内的平均空间密度分布函数其中,β和β′分别是纬度区间的最小值和最大值;
具体的,若β′>0,则直接执行步骤S7.2.1;
若β′≤0,则令β′=-β,β=-β′,均转换为关于原点对称的正值区间后,再执行步骤S7.2.1;
步骤S7.2.1,碎片轨道倾角记为i,当β′≤i且β≥-i时
<mrow>
<mover>
<mi>f</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>&beta;</mi>
<mo>,</mo>
<msup>
<mi>&beta;</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mrow>
<mi>arcsin</mi>
<mfrac>
<mrow>
<msup>
<mi>sin&beta;</mi>
<mo>&prime;</mo>
</msup>
</mrow>
<mrow>
<mi>sin</mi>
<mi>i</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mi>arcsin</mi>
<mfrac>
<mrow>
<mi>sin</mi>
<mi>&beta;</mi>
</mrow>
<mrow>
<mi>sin</mi>
<mi>i</mi>
</mrow>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<mi>sin&beta;</mi>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<mi>sin</mi>
<mi>&beta;</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>30</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤S7.2.2,当β≤-i且β′≤i时
<mrow>
<mover>
<mi>f</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>&beta;</mi>
<mo>,</mo>
<msup>
<mi>&beta;</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mrow>
<mi>arcsin</mi>
<mfrac>
<mrow>
<msup>
<mi>sin&beta;</mi>
<mo>&prime;</mo>
</msup>
</mrow>
<mrow>
<mi>sin</mi>
<mi>i</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mi>&pi;</mi>
<mn>2</mn>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<mi>sin&beta;</mi>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<mi>sin</mi>
<mi>&beta;</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>31</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤S7.2.3,当β′≥i且-i≤β≤i时
<mrow>
<mover>
<mi>f</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>&beta;</mi>
<mo>,</mo>
<msup>
<mi>&beta;</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mi>&pi;</mi>
<mn>2</mn>
</mfrac>
<mo>-</mo>
<mi>arcsin</mi>
<mfrac>
<mrow>
<mi>sin</mi>
<mi>&beta;</mi>
</mrow>
<mrow>
<mi>sin</mi>
<mi>i</mi>
</mrow>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<mi>sin&beta;</mi>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<mi>sin</mi>
<mi>&beta;</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>32</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤S7.2.4,当β′≥i且β≤-i时
<mrow>
<mover>
<mi>f</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>&beta;</mi>
<mo>,</mo>
<msup>
<mi>&beta;</mi>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
<mrow>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>sin&beta;</mi>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&beta;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>33</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤S7.2.5,当β>i时,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611036194.6A CN107066641B (zh) | 2016-11-23 | 2016-11-23 | 大规模空间碎片分布演化的数值计算方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611036194.6A CN107066641B (zh) | 2016-11-23 | 2016-11-23 | 大规模空间碎片分布演化的数值计算方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107066641A CN107066641A (zh) | 2017-08-18 |
CN107066641B true CN107066641B (zh) | 2018-05-11 |
Family
ID=59618417
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611036194.6A Active CN107066641B (zh) | 2016-11-23 | 2016-11-23 | 大规模空间碎片分布演化的数值计算方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107066641B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107871047A (zh) * | 2017-11-21 | 2018-04-03 | 中国人民解放军战略支援部队航天工程大学 | 一种复杂空间系统安全管理平行计算方法 |
CN107992669B (zh) * | 2017-11-28 | 2021-02-02 | 中国人民解放军战略支援部队航天工程大学 | 一种航天器解体事件的类型判定方法和系统 |
CN108021765B (zh) * | 2017-12-18 | 2021-03-02 | 北京卫星环境工程研究所 | 激光烧蚀驱动不规则三维目标力学行为的计算方法 |
CN110175431B (zh) * | 2019-06-05 | 2022-05-20 | 哈尔滨工业大学 | 一种地固系下空间碎片空间密度确定方法 |
CN111707274B (zh) * | 2020-05-29 | 2022-01-18 | 南京航空航天大学 | 能量最优的航天器连续动态避障轨迹规划方法 |
CN112307615A (zh) * | 2020-10-27 | 2021-02-02 | 西北工业大学 | 一种空间碎片轨道快速演化方法 |
CN112363410B (zh) * | 2020-11-13 | 2022-09-30 | 浙江大学 | 航天器智能自主控制研究与验证系统 |
CN113642785B (zh) * | 2021-07-28 | 2023-10-20 | 中国测绘科学研究院 | 基于先验信息的空间碎片轨道长期预报方法、系统及设备 |
CN115285381B (zh) * | 2022-10-09 | 2023-01-20 | 北京开运联合信息技术集团股份有限公司 | 一种太空碎片的碰撞预警方法及装置 |
CN115662193A (zh) * | 2022-10-19 | 2023-01-31 | 中国民航大学 | 一种面向空中交通管制的亚轨道碎片危险区生成方法 |
CN117744502A (zh) * | 2024-02-07 | 2024-03-22 | 中国人民解放军战略支援部队航天工程大学 | 一种基于兵棋的轨道碎片演化方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20120089517A (ko) * | 2010-12-15 | 2012-08-13 | 한국항공우주연구원 | 인공위성과 우주파편 사이의 충돌위험관리 시스템 및 방법 |
CN105223170A (zh) * | 2014-05-30 | 2016-01-06 | 中国科学院空间科学与应用研究中心 | 一种模拟微小空间碎片撞击诱发放电的装置及方法 |
-
2016
- 2016-11-23 CN CN201611036194.6A patent/CN107066641B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20120089517A (ko) * | 2010-12-15 | 2012-08-13 | 한국항공우주연구원 | 인공위성과 우주파편 사이의 충돌위험관리 시스템 및 방법 |
CN105223170A (zh) * | 2014-05-30 | 2016-01-06 | 中国科学院空间科学与应用研究中心 | 一种模拟微小空间碎片撞击诱发放电的装置及方法 |
Non-Patent Citations (5)
Title |
---|
An analytic method of space debris cloud evolution and its collision evaluation for constellation satellites;Zhang Binbin 等;《ADVANCES IN SPACE RESEARCH》;20160316;第903-913页 * |
LEGEND – a three-dimensional LEO-to-GEO debris evolutionary model;J.-C. Liou等;《ADVANCES IN SPACE RESEARCH》;20041231;第981-986页 * |
NASA’S NEW BREAKUP MODEL OF EVOLVE 4.0;N. L. Johnson等;《ADVANCES IN SPACE RESEARCH》;20011231;第28卷(第9期);第1377-1384页 * |
空间物体解体碎片云的长期演化建模与分析;张斌斌等;《中国空间科学技术》;20160825;第36卷(第4期);第1-8页 * |
空间碎片环境模型研究;王若璞;《中国博士学位论文全文数据库工程科技Ⅱ辑》;20110715;C031-15 * |
Also Published As
Publication number | Publication date |
---|---|
CN107066641A (zh) | 2017-08-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107066641B (zh) | 大规模空间碎片分布演化的数值计算方法及系统 | |
Wang | Analysis of Debris from the Collision of the Cosmos 2251 and the Iridium 33 Satellites | |
Reiland et al. | Assessing and minimizing collisions in satellite mega-constellations | |
CN106709145B (zh) | 大规模空间碎片分布状态数值演化的并行计算方法 | |
Armellin et al. | End-of-life disposal of high elliptical orbit missions: The case of INTEGRAL | |
Letizia et al. | Collision probability due to space debris clouds through a continuum approach | |
Dell’Elce et al. | Probabilistic assessment of the lifetime of low-earth-orbit spacecraft: uncertainty characterization | |
Pu et al. | Optimal small satellite orbit design based on robust multi-objective optimization method | |
Khoury | Orbital rendezvous and spacecraft loitering in the earth-moon system | |
Ray et al. | A drag coefficient modeling approach using spatial and temporal Fourier expansions for orbit determination | |
Wen et al. | Natural landing dynamics near the secondary in single-tidal-locked binary asteroids | |
Woo et al. | On the planar motion in the full two-body problem with inertial symmetry | |
Navabi et al. | Close approach analysis of space objects and estimation of satellite-debris collision probability | |
Rossi et al. | Collision risk against space debris in Earth orbits | |
Popova et al. | Planetary dynamics in the system α Centauri: The stability diagrams | |
Liu et al. | Analytical propagation of space debris density for collisions near sun-synchronous orbits | |
Zhang et al. | Discrete evolution model based on mean spatial density for space debris environment | |
Anderson et al. | Methodology for characterizing high-risk orbital debris in the geosynchronous orbit regime | |
Guardabasso et al. | Cislunar debris mitigation: Development of a methodology to assess the sustainability of lunar missions | |
Zhang et al. | Collision risk investigation for an operational spacecraft caused by space debris | |
Awad et al. | The method of multiple scales for orbit propagation with atmospheric drag | |
Sanchez et al. | Accessibility of the resources of near Earth space using multi-impulse transfers | |
Wang et al. | The assessment of the semi-analytical method in the long-term orbit prediction of Earth satellites | |
Machuca et al. | Robust optimization of descent trajectories on irregular-shaped bodies in the presence of uncertainty | |
Rich | Investigating analytical and numerical methods to predict satellite orbits using two-line element sets |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |