CN105511493A - 一种基于火星大气辅助的低轨星座部署方法 - Google Patents
一种基于火星大气辅助的低轨星座部署方法 Download PDFInfo
- Publication number
- CN105511493A CN105511493A CN201511000976.XA CN201511000976A CN105511493A CN 105511493 A CN105511493 A CN 105511493A CN 201511000976 A CN201511000976 A CN 201511000976A CN 105511493 A CN105511493 A CN 105511493A
- Authority
- CN
- China
- Prior art keywords
- aircraft
- orbit
- detector
- constellation
- air
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 76
- 230000008569 process Effects 0.000 claims abstract description 20
- 230000005484 gravity Effects 0.000 claims abstract description 4
- 238000012546 transfer Methods 0.000 claims description 29
- 238000013459 approach Methods 0.000 claims description 23
- 230000001174 ascending effect Effects 0.000 claims description 18
- 238000005457 optimization Methods 0.000 claims description 17
- 230000035807 sensation Effects 0.000 claims description 9
- 206010048669 Terminal state Diseases 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 3
- 101150036841 minJ gene Proteins 0.000 claims description 3
- 230000000541 pulsatile effect Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 5
- 230000008859 change Effects 0.000 description 7
- 239000000446 fuel Substances 0.000 description 6
- 230000009466 transformation Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/10—Simultaneous control of position or course in three dimensions
- G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
- G05D1/104—Simultaneous control of position or course in three dimensions specially adapted for aircraft involving a plurality of aircrafts, e.g. formation flying
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开的一种基于火星大气辅助的低轨星座部署方法,涉及一种火星大气与其引力系统下的火星星座部署方法,属于航空航天技术领域。本发明通过优化得到满足气动力要求的控制率来求解出所需初始轨道进入大气施加的速度脉冲和飞行器进入目标轨道施加的速度脉冲。探测器通过施加所需的飞行器从初始轨道进入大气速度脉冲将载有的飞行器从远火点位置释放并进入大气,在大气内通过优化给出的控制率进行气动力辅助轨道转移,并通过施加所需的飞行器进入目标轨道施加的速度脉冲将飞行器定轨到目标轨道上,将多颗星座飞行器分别部署到各自的目标轨道上,实现对整个星座的部署。本发明部署过程所消耗的能量低,对目标轨道没有严苛的要求,部署过程灵活。
Description
技术领域
本发明涉及一种火星大气与其引力系统下的火星星座部署方法,尤其涉及适用于火星探测器给火星部署小飞行器星座过程的低轨火星星座部署方法,属于航空航天技术领域。
背景技术
星座是部署入轨能正常工作的卫星集合,火星低轨星座是分布在火星低轨道上为实现某些特定任务而存在的卫星群。火星是太阳系中的八大行星之一,由于其距离地球较近,并且与地球有相近的自转周期,因此被认为是最有可能未来人来移居的且最有探测价值的行星。在八大行星中,除了地球,火星是最早被探测的行星。通过给火星部署星座,不仅可以长时间观测火星,而且对于未来的火星探测任务,如果有火星星座的存在,可以给着陆器提供导航以及通信服务,因此火星星座的存在对于未来火星探测以及相关的任务都有重大的意义。火星和地球一样,都具有大气层,由于借助大气进行轨道机动相比于传统的霍曼转移能节省大量的燃料,所以借助火星大气层进行星座的部署能有效的节省部署过程的能量消耗,从而为后续星座的维持节省更多的燃料。
在已发展的关于将飞行器从初始轨道部署到目标轨道的方法中在先技术[1](参见McCUEG.A.OptimumTwo-impulseOrbitalTransferandRendezvousBetweenInclinedEllipticalOrbits[J].AIAAJournal,Vol.1,No.8,1963,pp.1865-1872.)给出采用能量等高线方法,对飞行器轨道转移过程进行优化,通过能量等高线图选取能量最优的转移轨道,从而将飞行器从初始轨道部署到目标轨道上去。对于这种传统的轨道转移过程,轨道机动所消耗的燃料总量非常大,对于存在轨道面改变的轨道转移过程,在现有飞行器燃料配置的前提下几乎很难实现。
在先技术[2](参见BeererJ.,etal.AerobrakingatMars:theMGSMission[J].AIAAJournal,1996,10.2514/6.)给出采用大气阻力实现飞行器轨道转移的方法,也即气动辅助转移方法。对于火星探测器捕获任务,通过选择飞行器进入大气时的初始状态以及在大气内飞行过程中控制角的大小,使卫星在穿越大气之后,到达目标轨道,并通过一次较小脉冲机动实现飞行器转移至目标轨道的目的。该方法在轨道转移过程中消耗的燃料较少,但并没有给出如何通过优化的方法将多颗飞行器部署到火星低轨从而形成星座。
在星座部署过程中,星座飞行器全部由探测器释放,将星座飞行器精确的的部署到对应的目标轨道上是整个星座构建过程中最为关键的过程之一。
发明内容
本发明公开的一种基于火星大气辅助的低轨星座部署方法,要解决的技术问题是实现将星座飞行器精确的的部署到对应的目标轨道上,且可实现部署过程所消耗的能量低,对目标轨道的高度和倾角并没有严苛的要求,部署过程灵活。
本发明的目的是通过下述技术方案实现的:
本发明公开的一种基于火星大气辅助的低轨星座部署方法,通过优化得到满足气动力要求的控制率来求解出所需初始轨道进入大气施加的速度脉冲和飞行器进入目标轨道施加的速度脉冲。探测器通过施加所需的飞行器从初始轨道进入大气速度脉冲将载有的飞行器从远火点位置释放并进入大气,在大气内通过优化给出的控制率进行气动力辅助轨道转移,并通过施加所需的飞行器进入目标轨道施加的速度脉冲将飞行器定轨到目标轨道上,将多颗星座飞行器分别部署到各自的目标轨道上,从而实现对整个星座的部署。本发明公开的一种基于火星大气辅助的低轨星座部署方法,结合行星-星座卫星二体系统和气动力辅助变轨需求,利用火星大气阻力辅助作用代替直接的脉冲机动,部署过程所消耗的能量低,且对目标轨道的高度和倾角并没有严苛的要求,部署过程灵活。
本发明公开的一种基于火星大气辅助的低轨星座部署方法,包括如下步骤:
步骤一:在火星赤道惯性坐标系下建立探测器运动方程,确定探测器的轨道根数,从而确定探测器远火点的位置。所述的探测器载有需要低轨星座部署的多个星座飞行器。载有星座飞行器的数量根据星座任务需求而定。
步骤1.1:在火星赤道惯性坐标系下建立探测器运动方程;
选择火星的质心作为原点建立赤道惯性坐标系,选择X轴为行星赤道面内并指向春分点,Z轴指向火星自转轴的方向,并与火星旋转角速度的方向一致,Y轴在火星赤道面内垂直于X轴和Z轴,并与X轴和Z轴构成右手直角坐标系。
探测器在赤道惯性坐标系下的运动方程可以表示为,
其中μ表示火星的引力常数,xyz分别表示探测器在惯性坐标系OXYZ下的位置坐标,r表示探测器到引力中心的距离:
步骤1.2:定义探测器轨道六根数。
探测器轨道六根数分别定义为半长轴a、偏心率e、轨道倾角i、升交点赤经Ω、近地点幅角ω和真近点角θ。
I半长轴a
确定圆锥曲线轨道大小的参数,其几何意义是椭圆轨道半长轴。
II偏心率e
确定圆锥曲线轨道形状的参数,其几何意义是椭圆轨道偏心率。
III轨道倾角i
赤道平面与轨道平面的夹角。方向由右手法则确定,即由赤道平面向轨道平面绕交线矢量逆时针测量,i∈[0°,180°]。
IV升交点赤经Ω
由春分点沿着赤道至升交点的角度,Ω∈[0°,360°]。所述的升交点为探测器由南半球至北半球穿过赤道平面的点。
V近地点幅角ω
自轨道升交点,在轨道平面内沿探测器运动方向度量至近火点的角度,即近火点矢径与交线矢径之间的夹角,ω∈[0°,360°]。
VI真近点角θ
自近火点沿探测器运动方向度量至探测器某时刻所在位置的角度,是随着时间变化的,θ∈[0°,360°]。
步骤1.3:根据步骤1.2定义的探测器轨道六根数,给出赤道惯性系下的坐标到轨道根数的求解方法,并确定探测器远火点的位置。
首先,分别给出探测器在赤道惯性坐标系下坐标r=[X,Y,Z],v=[vX,vY,vZ]。
径向速度vr:
如果vr>0,则探测器正飞离近火点;若vr<0,则探测器正飞向近火点。
比角动量h:
轨道平面和赤道面交线矢量N:
I根据公式(5)轨道半长轴a
II根据公式(6)轨道偏心率e:
III轨道倾角i:
此处如果i∈[0,90°]则为顺行轨道,否则为逆行轨道。
IV升交点赤经Ω:
V近火点幅角ω:
VI真近点角θ:
根据式(5)和(6)可以得到探测器在远火点的位置为:
ra=a(1+e)(11)
步骤二:确定火星大气边缘高度以及从初始椭圆轨道进入大气所需脉冲Δv1,并确定飞行器转移到大气边缘时的速度v0与航迹角γ0。
给定火星大气边缘位置矢径为r0。
由于在远火点加机动进入大气所需脉冲量Δv1最小,所以在远火点给探测器载有的飞行器施加一次脉冲机动使得飞行器进入大气,根据椭圆轨道能量方程(12)得到在远火点探测器的速度v。
远火点探测器的速度v为:
当给出进入大气的转移轨道的近地点为rpe,则可以求得飞行器从初始轨道进入大气所需的速度脉冲Δv1为:
根据初始轨道远火点高度ra、进入大气的转移轨道的近地点rpe和火星大气边缘位置矢径r0可以求得飞行器转移到大气边缘时的速度v0与航迹角γ0:
步骤三:确定气动力辅助转移过程的运动方程、控制量、终端约束和定轨速度脉冲Δv2。
飞行器在火星大气内的运动如方程(17)所示:
其中,V为飞行器速度,r为飞行器矢径,γ为飞行航迹角,ψ为飞行航向角,θ为飞行器相对火星经度,为飞行器相对火星纬度。m为飞行器质量,μe为火星引力常数,Isp,g0分别为发动机比冲和重力加速度。α为攻角,σ为滚转角,T为发动机推力,所述的三个变量攻角α、滚转角σ、发动机推力T均属于控制量。若考虑无推力的气动力辅助轨道转移过程,则取T=0。
给出每个飞行器出大气之后的倾角if和目标轨道高度hf,其中目标轨道高度hf对应的矢径为rf,从而确定出大气时刻飞行器终端状态需满足的约束,即为:
此时,当飞行器飞出大气之后,飞行器沿转移轨道转移至目标轨道所在高度,之后通过施加第二次速度脉冲Δv2使飞行器进入目标轨道,所需施加第二次速度脉冲Δv2大小为:
步骤四:给出飞行器部署过程的优化性能指标,并根据优化性能指标和优化方法给确定优化后的所需初始轨道进入大气施加的速度脉冲Δv1和飞行器进入目标轨道施加的速度脉冲Δv2,并给出飞行器在大气内飞行的控制率。
给出飞行器部署过程优化性能指标:
minJ=Δv1+Δv2(21)
根据优化性能指标和优化方法给确定优化后的所需初始轨道进入大气施加的速度脉冲Δv1和飞行器进入目标轨道施加的速度脉冲Δv2。
对初始轨道进入大气施加的速度脉冲Δv1和飞行器进入目标轨道施加的速度脉冲Δv2进行优化,可通过求解两点边值问题进行求解,所述的优化求解方法可采用直接法或间接法。
步骤五:根据步骤四的优化结果实现低轨星座部署。
通过步骤四的优化可以得到步骤二需确定的飞行器从初始轨道进入大气施加的速度脉冲Δv1、步骤三需确定的飞行器进入目标轨道施加的速度脉冲Δv2和飞行器在大气内飞行的控制率。探测器通过施加所需的飞行器从初始轨道进入大气速度脉冲Δv1将载有的飞行器从远火点位置释放并进入大气,在大气内通过优化给出的控制率进行气动力辅助轨道转移,并通过施加所需的飞行器进入目标轨道施加的速度脉冲Δv2将飞行器定轨到目标轨道上,从而实现低轨星座部署。
有益效果:
1、本发明公开的一种基于火星大气辅助的低轨星座部署方法,通过优化得到满足气动力要求的控制率来求解出Δv1和Δv2。相比于传统的直接利用发动机推力进行飞行器部署所需的Δv1和Δv2,借助气动力求解得到的Δv1和Δv2可大幅减小。即由于本发明的低轨星座部署方法采用气动力代替或部分代替发动机的推力,相比直接施加机动的方法,星座部署过程所需速度增量Δv1和Δv2很小,进而节省燃料。
2、本发明公开的一种基于火星大气辅助的低轨星座部署方法,由于经过优化给出控制率,气动辅助过程灵活性高,对目标轨道的适用范围广。
3、本发明公开的一种基于火星大气辅助的低轨星座部署方法,鲁棒性强,可重复性高,受星座构型限制影响小,适用范围广。
附图说明
图1是本发明步骤1火星赤道惯性坐标系的示意图。
图2是基于火星大气辅助的低轨星座部署流程图。
图3是本发明实施例中飞行器纬度变化曲线。
图4是本发明实施例中飞行器航迹角变化曲线。
图5是本发明实施例中飞行器高度变化曲线。
图6是本发明实施例中飞行器航向角变化曲线。
具体实施方式
为了更好的说明本发明的目的和优点,下面结合附图和实例对发明内容做进一步说明。
实施例1:
如图2所示,本实施例公开的一种基于火星大气辅助的低轨星座部署方法,包括如下步骤:
步骤一:在火星赤道惯性坐标系下建立探测器运动方程,确定探测器的轨道根数,从而确定探测器远火点的位置。所述的探测器载有需要低轨星座部署的多个星座飞行器。载有星座飞行器的数量根据星座任务需求而定。
步骤1.1:在火星赤道惯性坐标系下建立探测器运动方程;
如图2所示,选择火星的质心作为原点建立赤道惯性坐标系,选择X轴为行星赤道面内并指向春分点,Z轴指向火星自转轴的方向,并与火星旋转角速度的方向一致,Y轴在火星赤道面内垂直于X轴和Z轴,并与X轴和Z轴构成右手直角坐标系。
探测器在赤道惯性坐标系下的运动方程可以表示为,
其中μ=42828km3/s2表示火星的引力常数,xyz分别表示探测器在惯性坐标系OXYZ下的位置坐标,r表示探测器到引力中心的距离:
步骤1.2:定义探测器轨道六根数。
探测器轨道六根数分别定义为半长轴a、偏心率e、轨道倾角i、升交点赤经Ω、近地点幅角ω和真近点角θ。
I半长轴a
确定圆锥曲线轨道大小的参数,其几何意义是椭圆轨道半长轴。
II偏心率e
确定圆锥曲线轨道形状的参数,其几何意义是椭圆轨道偏心率。
III轨道倾角i
赤道平面与轨道平面的夹角。方向由右手法则确定,即由赤道平面向轨道平面绕交线矢量逆时针测量,i∈[0°,180°]。
IV升交点赤经Ω
由春分点沿着赤道至升交点的角度,Ω∈[0°,360°]。所述的升交点为探测器由南半球至北半球穿过赤道平面的点。
V近地点幅角ω
自轨道升交点,在轨道平面内沿探测器运动方向度量至近火点的角度,即近火点矢径与交线矢径之间的夹角,
VI真近点角θ
自近火点沿探测器运动方向度量至探测器某时刻所在位置的角度,是随着时间变化的,θ∈[0°,360°]。
步骤1.3:根据步骤1.2定义的探测器轨道六根数,给出赤道惯性系下的坐标到轨道根数的求解方法,并确定探测器远火点的位置。
首先,分别给出探测器在赤道惯性坐标系下坐标r=[X,Y,Z],v=[vX,vY,vZ]。
径向速度vr:
如果vr>0,则探测器正飞离近火点;若vr<0,则探测器正飞向近火点。
比角动量h:
轨道平面和赤道面交线矢量N:
I根据公式(5)轨道半长轴a
II根据公式(6)轨道偏心率e:
III轨道倾角i:
此处如果则为顺行轨道,否则为逆行轨道。
IV升交点赤经Ω:
V近火点幅角ω:
VI真近点角θ:
根据式(5)和(6)可以得到探测器在远火点的位置为:
ra=a(1+e)(11)
步骤二:确定火星大气边缘高度以及从初始椭圆轨道进入大气所需脉冲Δv1,并确定飞行器转移到大气边缘时的速度v0与航迹角γ0。
给定火星大气边缘位置矢径为r0。
由于在远火点加机动进入大气所需脉冲量Δv1最小,所以在远火点给探测器载有的飞行器施加一次脉冲机动使得飞行器进入大气,根据椭圆轨道能量方程(12)得到在远火点探测器的速度v。
远火点探测器的速度v为:
当给出进入大气的转移轨道的近地点为rpe,则可以求得飞行器从初始轨道进入大气所需的速度脉冲Δv1为:
根据初始轨道远火点高度ra、进入大气的转移轨道的近地点rpe和火星大气边缘位置矢径r0可以求得飞行器转移到大气边缘时的速度v0与航迹角γ0:
步骤三:确定气动力辅助转移过程的运动方程、控制量、终端约束和定轨速度脉冲Δv2。
飞行器在火星大气内的运动如方程(17)所示:
其中,V为飞行器速度,r为飞行器矢径,γ为飞行航迹角,ψ为飞行航向角,θ为飞行器相对火星经度,为飞行器相对火星纬度。飞行器质量m为50kg,火星引力常数μ为42828km3/s2,Isp,g0分别为发动机比冲和重力加速度。α为攻角,σ为滚转角,T为发动机推力,所述的三个变量攻角α、滚转角σ、发动机推力T均属于控制量。此处考虑无推力的气动力辅助轨道转移过程,所以取T=0。
给出每个飞行器出大气之后的倾角if和目标轨道高度hf,其中目标轨道高度hf对应的矢径为rf,从而确定出大气时刻飞行器终端状态需满足的约束,即为:
此时,当飞行器飞出大气之后,飞行器沿转移轨道转移至目标轨道所在高度,之后通过施加第二次速度脉冲Δv2使飞行器进入目标轨道,所需施加第二次速度脉冲Δv2大小为:
步骤四:给出飞行器部署过程的优化性能指标,并根据优化性能指标和优化方法给确定优化后的所需初始轨道进入大气施加的速度脉冲Δv1和飞行器进入目标轨道施加的速度脉冲Δv2,并给出飞行器在大气内飞行的控制率。
给出飞行器部署过程优化性能指标:
minJ=Δv1+Δv2(21)
根据优化性能指标和优化方法给确定优化后的所需初始轨道进入大气施加的速度脉冲Δv1和飞行器进入目标轨道施加的速度脉冲Δv2。
对初始轨道进入大气施加的速度脉冲Δv1和飞行器进入目标轨道施加的速度脉冲Δv2进行优化,可通过求解两点边值问题进行求解,所述的优化求解方法可采用直接法或间接法。
步骤五:根据步骤四的优化结果实现低轨星座部署。
通过步骤四的优化可以得到步骤二需确定的飞行器从初始轨道进入大气施加的速度脉冲Δv1、步骤三需确定的飞行器进入目标轨道施加的速度脉冲Δv2和飞行器在大气内飞行的控制率。探测器通过施加所需的飞行器从初始轨道进入大气速度脉冲Δv1将载有的飞行器从远火点位置释放并进入大气,在大气内通过优化给出的控制率进行气动力辅助轨道转移,并通过施加所需的飞行器进入目标轨道施加的速度脉冲Δv2将飞行器定轨到目标轨道上,从而实现低轨星座部署。
为了验证方法的可行性,选择探测器的轨道为椭圆轨道,考虑探测器在远火点释放3颗飞行器,并经过火星大气辅助变轨之后分别部署到各自的目标轨道上,形成三星星座。假设火星的半径为3389km。
选择探测器所在轨道半长轴为32910km,偏心率为0.8887,轨道倾角为0°。
选择三颗星座飞行器的目标轨道为:
飞行器1:高度为400km,倾角为10°的圆轨道;
飞行器2:高度为400km,倾角为5°的圆轨道;
飞行器3:高度为400km,倾角为1°的圆轨道;
根据探测器的轨道根数可以求得探测器的远火点高度为58762.736km。
根据直接法求解方程(17)可以得到每个飞行器在大气内飞行对应状态变量随时间的变化关系,并得到相应的第一次轨道转移的近火点高度hpe和离轨和最后定轨所施加的速度脉冲Δv1,Δv2,以及总消耗的速度脉冲Δv。
表1部署过程飞行器对应的参数与燃耗
尽管已经示出和描述了本发明的实施例,本领域的普通技术人员可以理解:在不脱离本发明的原理和宗旨的情况下,可对这些实施例进行多种变化、修改、替换和变形,本发明的范围由权利要求及其等同物限定。
Claims (8)
1.一种基于火星大气辅助的低轨星座部署方法,其特征在于:包括如下步骤,
步骤一:在火星赤道惯性坐标系下建立探测器运动方程,确定探测器的轨道根数,从而确定探测器远火点的位置;
步骤二:确定火星大气边缘高度以及从初始椭圆轨道进入大气所需脉冲△v1,并确定飞行器转移到大气边缘时的速度v0与航迹角γ0;
给定火星大气边缘位置矢径为r0;
由于在远火点加机动进入大气所需脉冲量△v1最小,所以在远火点给探测器载有的飞行器施加一次脉冲机动使得飞行器进入大气,根据椭圆轨道能量方程(12)得到在远火点探测器的速度v;
远火点探测器的速度v为:
当给出进入大气的转移轨道的近地点为rpe,则可以求得飞行器从初始轨道进入大气所需的速度脉冲△v1为:
根据初始轨道远火点高度ra、进入大气的转移轨道的近地点rpe和火星大气边缘位置矢径r0可以求得飞行器转移到大气边缘时的速度v0与航迹角γ0:
步骤三:确定气动力辅助转移过程的运动方程、控制量、终端约束和定轨速度脉冲△v2;
飞行器在火星大气内的运动如方程(17)所示:
其中,V为飞行器速度,r为飞行器矢径,γ为飞行航迹角,ψ为飞行航向角,θ为飞行器相对火星经度,为飞行器相对火星纬度;m为飞行器质量,μe为火星引力常数,Isp,g0分别为发动机比冲和重力加速度;α为攻角,σ为滚转角,T为发动机推力,所述的三个变量攻角α、滚转角σ、发动机推力T均属于控制量;若考虑无推力的气动力辅助轨道转移过程,则取T=0;
给出每个飞行器出大气之后的倾角if和目标轨道高度hf,其中目标轨道高度hf对应的矢径为rf,从而确定出大气时刻飞行器终端状态需满足的约束,即为:
此时,当飞行器飞出大气之后,飞行器沿转移轨道转移至目标轨道所在高度,之后通过施加第二次速度脉冲△v2使飞行器进入目标轨道,所需施加第二次速度脉冲△v2大小为:
步骤四:给出飞行器部署过程的优化性能指标,并根据优化性能指标和优化方法给确定优化后的所需初始轨道进入大气施加的速度脉冲△v1和飞行器进入目标轨道施加的速度脉冲△v2,并给出飞行器在大气内飞行的控制率;
给出飞行器部署过程优化性能指标:
minJ=△v1+△v2(21)
根据优化性能指标和优化方法给确定优化后的所需初始轨道进入大气施加的速度脉冲△v1和飞行器进入目标轨道施加的速度脉冲△v2;
步骤五:根据步骤四的优化结果实现低轨星座部署;
通过步骤四的优化可以得到步骤二需确定的飞行器从初始轨道进入大气施加的速度脉冲△v1、步骤三需确定的飞行器进入目标轨道施加的速度脉冲△v2和飞行器在大气内飞行的控制率;探测器通过施加所需的飞行器从初始轨道进入大气速度脉冲△v1将载有的飞行器从远火点位置释放并进入大气,在大气内通过优化给出的控制率进行气动力辅助轨道转移,并通过施加所需的飞行器进入目标轨道施加的速度脉冲△v2将飞行器定轨到目标轨道上,从而实现低轨星座部署。
2.如权利要求1所述的一种基于火星大气辅助的低轨星座部署方法,其特征在于:
步骤一的具体实现步骤为,
步骤1.1:在火星赤道惯性坐标系下建立探测器运动方程;
选择火星的质心作为原点建立赤道惯性坐标系,选择X轴为行星赤道面内并指向春分点,Z轴指向火星自转轴的方向,并与火星旋转角速度的方向一致,Y轴在火星赤道面内垂直于X轴和Z轴,并与X轴和Z轴构成右手直角坐标系;
探测器在赤道惯性坐标系下的运动方程可以表示为,
其中μ表示火星的引力常数,xyz分别表示探测器在惯性坐标系OXYZ下的位置坐标,r表示探测器到引力中心的距离:
步骤1.2:定义探测器轨道六根数;
探测器轨道六根数分别定义为半长轴a、偏心率e、轨道倾角i、升交点赤经Ω、近地点幅角ω和真近点角θ;
I半长轴a
确定圆锥曲线轨道大小的参数,其几何意义是椭圆轨道半长轴;
II偏心率e
确定圆锥曲线轨道形状的参数,其几何意义是椭圆轨道偏心率;
III轨道倾角i
赤道平面与轨道平面的夹角;方向由右手法则确定,即由赤道平面向轨道平面绕交线矢量逆时针测量,i∈[0°,180°];
IV升交点赤经Ω
由春分点沿着赤道至升交点的角度,Ω∈[0°,360°];所述的升交点为探测器由南半球至北半球穿过赤道平面的点;
V近地点幅角ω
自轨道升交点,在轨道平面内沿探测器运动方向度量至近火点的角度,即近火点矢径与交线矢径之间的夹角,ω∈[0°,360°];
VI真近点角θ
自近火点沿探测器运动方向度量至探测器某时刻所在位置的角度,是随着时间变化的,θ∈[0°,360°];
步骤1.3:根据步骤1.2定义的探测器轨道六根数,给出赤道惯性系下的坐标到轨道根数的求解方法,并确定探测器远火点的位置;
首先,分别给出探测器在赤道惯性坐标系下坐标r=[X,Y,Z],v=[vX,vY,vZ];
径向速度vr:
如果vr>0,则探测器正飞离近火点;若vr<0,则探测器正飞向近火点;
比角动量h:
轨道平面和赤道面交线矢量N:
I根据公式(5)轨道半长轴a
II根据公式(6)轨道偏心率e:
III轨道倾角i:
此处如果i∈[0,90°]则为顺行轨道,否则为逆行轨道;
IV升交点赤经Ω:
V近火点幅角ω:
VI真近点角θ:
根据式(5)和(6)可以得到探测器在远火点的位置为。
ra=a(1+e)(11)。
3.如权利要求1或2所述的一种基于火星大气辅助的低轨星座部署方法,其特征在于:所述的探测器载有需要低轨星座部署的多个星座飞行器;载有星座飞行器的数量根据星座任务需求而定。
4.如权利要求1或2所述的一种基于火星大气辅助的低轨星座部署方法,其特征在于:步骤四中对初始轨道进入大气施加的速度脉冲△v1和飞行器进入目标轨道施加的速度脉冲△v2进行优化,可通过求解两点边值问题进行求解,所述的优化求解方法采用直接法或间接法。
5.如权利要求3所述的一种基于火星大气辅助的低轨星座部署方法,其特征在于:步骤四中对初始轨道进入大气施加的速度脉冲△v1和飞行器进入目标轨道施加的速度脉冲△v2进行优化,可通过求解两点边值问题进行求解,所述的优化求解方法采用直接法或间接法。
6.一种基于火星大气辅助的低轨星座部署方法,其特征在于:通过优化得到满足气动力要求的控制率来求解出所需初始轨道进入大气施加的速度脉冲和飞行器进入目标轨道施加的速度脉冲;探测器通过施加所需的飞行器从初始轨道进入大气速度脉冲将载有的飞行器从远火点位置释放并进入大气,在大气内通过优化给出的控制率进行气动力辅助轨道转移,并通过施加所需的飞行器进入目标轨道施加的速度脉冲将飞行器定轨到目标轨道上,将多颗星座飞行器分别部署到各自的目标轨道上,从而实现对整个星座的部署。
7.如权利要求6所述的一种基于火星大气辅助的低轨星座部署方法,其特征在于:所述的探测器载有需要低轨星座部署的多个星座飞行器;载有星座飞行器的数量根据星座任务需求而定。
8.如权利要求6或7所述的一种基于火星大气辅助的低轨星座部署方法,其特征在于:对目标轨道的适用范围广,故受星座构型限制影响小,适用范围广。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511000976.XA CN105511493B (zh) | 2015-12-28 | 2015-12-28 | 一种基于火星大气辅助的低轨星座部署方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511000976.XA CN105511493B (zh) | 2015-12-28 | 2015-12-28 | 一种基于火星大气辅助的低轨星座部署方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105511493A true CN105511493A (zh) | 2016-04-20 |
CN105511493B CN105511493B (zh) | 2018-04-03 |
Family
ID=55719553
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511000976.XA Active CN105511493B (zh) | 2015-12-28 | 2015-12-28 | 一种基于火星大气辅助的低轨星座部署方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105511493B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106335655A (zh) * | 2016-10-31 | 2017-01-18 | 上海航天控制技术研究所 | 一种火星器器分离轨道设计方法 |
CN106383994A (zh) * | 2016-09-05 | 2017-02-08 | 北京理工大学 | 一种基于脉冲和气动辅助结合的低轨轨道面转移方法 |
CN107609267A (zh) * | 2017-09-12 | 2018-01-19 | 北京理工大学 | 一种月球有限推力多次捕获轨道实现方法 |
CN108491650A (zh) * | 2018-03-29 | 2018-09-04 | 北京理工大学 | 一种基于智能学习的火星进入终端状态高效评估方法 |
CN111301715A (zh) * | 2020-01-17 | 2020-06-19 | 哈尔滨工业大学 | 基于霍曼变轨的同轨道特定相位分布的星座布局与轨道调整方法、装置及计算机存储介质 |
CN113911398A (zh) * | 2021-11-09 | 2022-01-11 | 中国人民解放军火箭军工程大学 | 飞行器单脉冲躲避策略确定方法及系统 |
CN114139277A (zh) * | 2021-11-05 | 2022-03-04 | 中山大学 | 一种火星大气电离层穿刺探测轨迹优化方法、装置及设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2520629C1 (ru) * | 2012-11-19 | 2014-06-27 | Федеральное государственное унитарное предприятие "Центральный научно-исследовательский институт машиностроения" (ФГУП ЦНИИмаш) | Способ управления выведением космического аппарата на орбиту искусственного спутника планеты |
CN104019818A (zh) * | 2014-06-19 | 2014-09-03 | 北京理工大学 | 一种基于预测轨迹的行星导航轨道器布局优化方法 |
CN104794311A (zh) * | 2015-05-15 | 2015-07-22 | 哈尔滨工业大学 | 一种时间自由条件下的全局最优双脉冲轨道转移方法 |
CN105115512A (zh) * | 2015-09-23 | 2015-12-02 | 北京理工大学 | 一种火星大气进入段侧向预测校正制导方法 |
-
2015
- 2015-12-28 CN CN201511000976.XA patent/CN105511493B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2520629C1 (ru) * | 2012-11-19 | 2014-06-27 | Федеральное государственное унитарное предприятие "Центральный научно-исследовательский институт машиностроения" (ФГУП ЦНИИмаш) | Способ управления выведением космического аппарата на орбиту искусственного спутника планеты |
CN104019818A (zh) * | 2014-06-19 | 2014-09-03 | 北京理工大学 | 一种基于预测轨迹的行星导航轨道器布局优化方法 |
CN104794311A (zh) * | 2015-05-15 | 2015-07-22 | 哈尔滨工业大学 | 一种时间自由条件下的全局最优双脉冲轨道转移方法 |
CN105115512A (zh) * | 2015-09-23 | 2015-12-02 | 北京理工大学 | 一种火星大气进入段侧向预测校正制导方法 |
Non-Patent Citations (2)
Title |
---|
符俊 等: "气动力辅助双脉冲最优异面变轨问题研究", 《弹道学报》 * |
闫旭: "航天器气动力辅助变轨方法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106383994A (zh) * | 2016-09-05 | 2017-02-08 | 北京理工大学 | 一种基于脉冲和气动辅助结合的低轨轨道面转移方法 |
CN106335655A (zh) * | 2016-10-31 | 2017-01-18 | 上海航天控制技术研究所 | 一种火星器器分离轨道设计方法 |
CN106335655B (zh) * | 2016-10-31 | 2019-06-07 | 上海航天控制技术研究所 | 一种火星器器分离轨道设计方法 |
CN107609267A (zh) * | 2017-09-12 | 2018-01-19 | 北京理工大学 | 一种月球有限推力多次捕获轨道实现方法 |
CN107609267B (zh) * | 2017-09-12 | 2019-07-16 | 北京理工大学 | 一种月球有限推力多次捕获轨道实现方法 |
CN108491650A (zh) * | 2018-03-29 | 2018-09-04 | 北京理工大学 | 一种基于智能学习的火星进入终端状态高效评估方法 |
CN108491650B (zh) * | 2018-03-29 | 2021-05-18 | 北京理工大学 | 一种基于智能学习的火星进入终端状态高效评估方法 |
CN111301715A (zh) * | 2020-01-17 | 2020-06-19 | 哈尔滨工业大学 | 基于霍曼变轨的同轨道特定相位分布的星座布局与轨道调整方法、装置及计算机存储介质 |
CN111301715B (zh) * | 2020-01-17 | 2022-02-25 | 哈尔滨工业大学 | 基于霍曼变轨的同轨道特定相位分布的星座布局与轨道调整方法、装置及计算机存储介质 |
CN114139277A (zh) * | 2021-11-05 | 2022-03-04 | 中山大学 | 一种火星大气电离层穿刺探测轨迹优化方法、装置及设备 |
CN113911398A (zh) * | 2021-11-09 | 2022-01-11 | 中国人民解放军火箭军工程大学 | 飞行器单脉冲躲避策略确定方法及系统 |
CN113911398B (zh) * | 2021-11-09 | 2023-11-10 | 中国人民解放军火箭军工程大学 | 飞行器单脉冲躲避策略确定方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN105511493B (zh) | 2018-04-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Xie et al. | Spacecraft dynamics and control | |
CN105511493A (zh) | 一种基于火星大气辅助的低轨星座部署方法 | |
CN109774974B (zh) | 一种用于空间碎片接近的轨道设计方法 | |
CN106383994A (zh) | 一种基于脉冲和气动辅助结合的低轨轨道面转移方法 | |
Dutta et al. | Aerocapture as an enhancing option for ice giants missions | |
CN104309822A (zh) | 一种基于参数优化的航天器单脉冲水滴形绕飞轨迹悬停控制方法 | |
Lee et al. | Preliminary design of the guidance, navigation, and control system of the Altair Lunar lander | |
Omar et al. | Satellite formation control using differential drag | |
CN108082538B (zh) | 一种考虑始末约束的多体系统低能量捕获轨道方法 | |
Kang et al. | Nanosat formation flying design for SNIPE mission | |
Varma | Control of satellites using environmental forces: aerodynamic drag/solar radiation pressure | |
Xie et al. | Guidance, navigation, and control for spacecraft rendezvous and docking: theory and methods | |
CN104354877A (zh) | 一种基于地球-火星循环轨道的载人火星探测系统及方法 | |
Gao et al. | Analytical design methods for transfer trajectories between the Earth and the lunar orbital station | |
Nan et al. | Global 4-D trajectory optimization for spacecraft | |
CN113741551A (zh) | 一种基于代理模型的全过程轨迹优化方法及装置 | |
Newman | Drift recovery and station keeping results for the historic canx-4/canx-5 formation flying mission | |
Abilleira | 2011 Mars Science Laboratory Mission Design Overview | |
Zhang et al. | Theory and Design Methods of Special Space Orbits | |
Xie et al. | Autonomous guidance, navigation, and control of spacecraft | |
Liu et al. | Orbit transfer from Earth to Mars using modified Lambert algorithm | |
BLIDERAN et al. | RE-ENTRY STRATEGY ANALYSIS OF REUSABLE LAUNCH VEHICLES | |
Bettinger | Spacecraft Demand Tasking and Skip Entry Responsive Maneuvers | |
Ehricke | Orbital operations | |
Di Campli Bayard de Volo et al. | A Gliding Vehicle for ISS Crew Rescue-Mission Operational Concept |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |