CN107310752A - 一种太阳帆航天器行星圆悬浮轨道之间的转移方法 - Google Patents

一种太阳帆航天器行星圆悬浮轨道之间的转移方法 Download PDF

Info

Publication number
CN107310752A
CN107310752A CN201710369368.9A CN201710369368A CN107310752A CN 107310752 A CN107310752 A CN 107310752A CN 201710369368 A CN201710369368 A CN 201710369368A CN 107310752 A CN107310752 A CN 107310752A
Authority
CN
China
Prior art keywords
planet
suspension railway
solar sail
oval
railway
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
Application number
CN201710369368.9A
Other languages
English (en)
Other versions
CN107310752B (zh
Inventor
袁建平
陈建林
袁静
孙冲
高琛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201710369368.9A priority Critical patent/CN107310752B/zh
Publication of CN107310752A publication Critical patent/CN107310752A/zh
Application granted granted Critical
Publication of CN107310752B publication Critical patent/CN107310752B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Chemical & Material Sciences (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Combustion & Propulsion (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

本发明公开了一种太阳帆航天器行星圆悬浮轨道之间的转移方法,该方法首先分析了行星圆悬浮轨道和椭圆悬浮轨道的动力学模型,并给出了实现这两类轨道的太阳帆姿态参数和反射率参数求解方法;同时分析了初始、目标行星圆悬浮轨道与行星椭圆悬浮轨道之间进行拼接和转移的边界条件;然后根据任务要求,确定航天器进行转移的初始圆悬浮轨道和最终圆悬浮轨道轨道参数,并根据边界条件求解用于轨道转移的行星椭圆悬浮轨道,进而结合行星椭圆悬浮轨道的动力学平衡点和太阳光压力模型解算太阳帆姿态参数和反射率参数。该发明可针对行星圆悬浮轨道之间的转移问题,给出合适的行星椭圆悬浮轨道。

Description

一种太阳帆航天器行星圆悬浮轨道之间的转移方法
技术领域
本发明属于航空航天技术领域,涉及一种太阳帆航天器行星圆悬浮轨道之间的转移方法。
背景技术
近年来,随着先进材料技术、航天技术的发展和突破,太阳帆推进技术、电推进技术等新型的推进技术飞速发展。通常来说,太阳帆推进技术可借助太阳光产生推力而不消耗航天器本身的化学燃料;电推进装置则具有较大的比冲而在推进航天器的过程中可以减少化学燃料的消耗。因此,相比于传统的推进装置,这两类推进系统在长时间、高能量的航天任务中具有显著的优势,特别适用于需要采用强非开普勒轨道作为任务轨道的任务,以及实现强非开普勒轨道之间的转移任务。
相比于传统开普勒轨道,强非开普勒轨道在行星的极点观测任务、土星环和地磁场的观测任务,以及空间环境监测等任务中有着广泛的潜在应用价值。因此,大量的文献对基于太阳帆推进技术的日心圆悬浮轨道、行星圆悬浮轨道、人工拉格朗日点轨道建模理论进行了广泛的研究。对这类技术的在轨验证实验也相继展开,特别是2010年,日本宇航局首次将名为IKAROS的试验太阳帆航天器成功送入了太空并进行了大量的技术验证,这次具有划时代意义的轨道实验正式拉开了人类利用太阳帆装置推进航天器的序幕。
目前而言,基于太阳帆推进技术的任务轨道研究主要集中于日心、行星圆悬浮轨道,拉格朗日平动点轨道,以及日心椭圆悬浮轨道等,对行星椭圆悬浮轨道的研究仍较为少见。作为太阳帆能够实现的轨道的一部分,行星椭圆悬浮轨道在行星附近环境监测、实现行星圆悬浮轨道之间的转移等任务中具有潜在的应用价值。本发明针对行星圆悬浮轨道之间的转移问题,提出了一种利用行星椭圆悬浮轨道进行拼接和转移的方法。
发明内容
本发明的目的在于提供一种太阳帆航天器行星圆悬浮轨道之间的转移方法;通过利用太阳帆推进技术设计行星椭圆悬浮转移轨道,从而实现了行星圆悬浮轨道的转移任务。
本发明的目的是通过以下技术方案来实现的:
这种太阳帆航天器行星圆悬浮轨道之间的转移方法,包括以下步骤:
步骤1,确定太阳帆航天器行星椭圆悬浮轨道模型和行星圆悬浮轨道模型,并且确定两类轨道下太阳帆的唯一描述参数;
步骤2,通过行星圆悬浮轨道上的太阳帆航天器的位置矢量和速度矢量,确定行星椭圆悬浮轨道和行星圆悬浮轨道进行拼接的边界条件;
步骤3,将圆悬浮轨道之间的转移作为目标,得到实现行星圆悬浮轨道之间转移的椭圆悬浮轨道;
步骤4,通过边界条件确定椭圆悬浮轨道的轨道参数,根据椭圆悬浮轨道的平衡条件和太阳光压加速度模型,得到实现椭圆悬浮转移轨道的太阳帆姿态角参数和反射率参数;
步骤5,判定太阳帆姿态参数和反射率参数是否可实现,若能实现,则椭圆悬浮轨道可作为转移轨道;否则重复步骤1-4;
步骤6,完成太阳帆航天器在行星圆悬浮轨道之间的转移。
更进一步的,本发明的特点还在于:
其中步骤1中通过建立惯性参考坐标系非一致旋转相合坐标系和旋转坐标系确定两类轨道下太阳帆的唯一描述参数。
其中步骤2中行星椭圆悬浮轨道和行星圆悬浮轨道进行拼接的位置为椭圆悬浮轨道的近焦点或远焦点。
其中步骤3中两条行星圆悬浮轨道位于同一平面内时,得到实现行星圆悬浮轨道之间转移的椭圆悬浮轨道。
其中步骤3两条行星圆悬浮轨道不在同一平面内时,得到实现行星圆悬浮轨道之间转移的椭圆悬浮轨道。
其中步骤4中在非一致旋转相合坐标系下建立太阳帆航天器的动力学模型其中Ri是太阳帆相对于行星质心的位置矢量;R是太阳帆相对于坐标系原点的位置矢量;ω表示坐标系本身的旋转角速度;μp表示行星的引力常数;太阳光压加速度a表示为:a=ac(ns·n)·[2(1-u)ρs(ns·n)·n+(1-ρs+uρs)ns],单位矢量ns表示太阳光方向;n表示太阳帆的法向;u为镜面反射部分面积占太阳帆总面积的比值;ρs表示照射在镜面反射上的太阳光完全反射的比率;ac表示行星附近的最大太阳光压加速度。
其中步骤5中太阳帆姿态约束c(α,θ)和反射率参数u满足以下条件为太阳帆姿态参数和反射率参数能够实现:c=n·ns∈(0,1]:太阳帆的法向方向n必须指向太阳的反方向;u∈[0,umax]:太阳帆的反射率不能超过其最大反射率。
与现有技术相比,本发明的有益效果是:本发明的技术方案提出了太阳帆航天器行星圆悬浮轨道之间的转移轨道设计方法。首先给出行星椭圆悬浮轨道和圆悬浮轨道的动力学模型,明确描述这两类轨道的轨道参数和实现这两类轨道的太阳帆姿态参数和反射率参数解算方法;在此基础上,分析初始、目标行星圆悬浮轨道与行星椭圆悬浮轨道之间进行拼接和转移的边界条件;然后根据任务要求,确定航天器进行转移的初始圆悬浮轨道和目标圆悬浮轨道,并按照初始和目标轨道之间的位置关系进行分类,针对不同分类,根据边界条件求解用于轨道转移的行星椭圆悬浮轨道;最后,结合行星椭圆悬浮轨道的动力学平衡点和太阳光压力模型解算太阳帆姿态参数和反射率参数,并验证两者是否满足约束。
附图说明
图1为本发明的流程示意图;
图2为本发明中惯性参考坐标系示意图;
图3为本发明中非一致旋转相合坐标系示意图;
图4为本发明中共平面转移示意图;
图5为本发明中一般轨道转移情况一的坐标示意图;
图6为本发明中一般轨道转移情况二的坐标示意图;
图7为本发明中太阳帆姿态参数θ的变化示意图;
图8为本发明中太阳帆姿态参数α的变化示意图;
图9为本发明中太阳帆镜面反射率u的变化示意图。
具体实施方式
下面结合附图对本发明做进一步详细描述:
本发明提供了一种太阳帆航天器行星圆悬浮轨道之间的转移方法,包括以下步骤:
步骤1确定太阳帆航天器行星椭圆悬浮轨道模型和行星圆悬浮轨道模型,通过建立惯性参考坐标系非一致旋转相合坐标系和旋转坐标系从而确定两类轨道下太阳帆的唯一描述参数。
步骤2,通过行星圆悬浮轨道上的太阳帆航天器的位置矢量和速度矢量,确定行星椭圆悬浮轨道和行星圆悬浮轨道进行拼接的位置为椭圆悬浮轨道的近焦点或远焦点。
步骤3,将圆悬浮轨道作为转移目标,在两条行星圆悬浮轨道位于同一平面或不同平面内时,得到实现行星圆悬浮轨道之间转移的椭圆悬浮轨道。
步骤4,在非一致旋转相合坐标系下建立太阳帆航天器的动力学模型其中Ri是太阳帆相对于行星质心的位置矢量;R是太阳帆相对于坐标系原点的位置矢量;ω表示坐标系本身的旋转角速度;μp表示行星的引力常数;太阳光压加速度a表示为:a=ac(ns·n)·[2(1-u)ρs(ns·n)·n+(1-ρs+uρs)ns],单位矢量ns表示太阳光方向;n表示太阳帆的法向;u为镜面反射部分面积占太阳帆总面积的比值;ρs表示照射在镜面反射上的太阳光完全反射的比率;ac表示行星附近的最大太阳光压加速度。
步骤5,根据椭圆悬浮轨道的轨道参数和太阳光压加速度模型,得到作为转移轨道的椭圆悬浮轨道的太阳帆姿态角参数θ,α和反射率参数u;太阳帆姿态约束c(α,θ)和反射率参数u满足以下条件:c=n·ns∈(0,1]:太阳帆的法向方向n必须指向太阳的反方向;u∈[0,umax]:太阳帆的反射率不能超过其最大反射率;当太阳帆姿态约束c(α,θ)和反射率参数u满足上述条件时,太阳帆姿态参数和反射率参数能够实现,将椭圆悬浮轨道可作为转移轨道;不满足时,重复步骤1-4,重新确定椭圆悬浮轨道。
步骤6,通过以上步骤,实现太阳帆航天器在行星圆悬浮轨道之间的转移。
本发明的具体实施例为:
太阳帆航天器利用行星椭圆悬浮轨道实现圆悬浮轨道之间的转移的设计流程如图1所示。在行星椭圆悬浮轨道设计之前,首先分析初始、目标行星圆悬浮轨道与行星椭圆悬浮轨道之间进行拼接的边界条件;在此基础上,根据任务要求和航天器测量设备性能,确定航天器进行转移的初始圆悬浮轨道和目标圆悬浮轨道,并根据边界条件求解用于轨道转移的行星椭圆悬浮轨道轨道参数。在此基础上,建立一个新的旋转非一致相合坐标系,并利用该坐标系下太阳帆航天器的动力学平衡点条件,结合太阳光压力模型,解算太阳帆姿态参数和反射率参数。最后,结合约束条件,判断太阳帆参数是否具有物理意义,如果有意义,则该行星椭圆悬浮轨道可以作为初始和目标圆悬浮轨道之间的转移轨道,如果没有意义,则重新设计目标行星椭圆悬浮轨道。其具体过程为:
建立三个坐标系。惯性参考坐标系如图2所示,该坐标系原点o1位于行星质心,kz轴沿着太阳行星连线并指向行星;ix轴垂直于kz轴,且位于由kz轴和行星椭圆悬浮轨道主焦点组成的平面内,jy满足右手螺旋坐标系。非一致旋转相合坐标系如图3所示,该坐标系原点F位于行星椭圆悬浮轨道的主焦点,初始的iξ轴指向参考方向,该参考方向位于平面o1ixkz内;kζ轴垂直于行星椭圆悬浮轨道面;jη满足右手螺旋坐标系。该坐标系以kζ轴为旋转轴,具有与参考轨道上航天器瞬时角速度相同的旋转角速度iξ轴和参考方向之间的夹角为f。3)旋转坐标系该坐标系原点os位于太阳帆的质心,k轴垂直于行星椭圆悬浮轨道平面;i轴沿着矢径r的方向;j满足右手螺旋坐标系。太阳帆的姿态角在该坐标系内可由两个姿态角和α确定。
在行星圆悬浮轨道和椭圆悬浮轨道拼接过程中,假设太阳帆的姿态角和反射率可以进行瞬时切换,另外,在拼接点处不允许出现额外的推力。在上述条件下,为了实现行星圆悬浮轨道和行星椭圆悬浮轨道之间的拼接,要求拼接点处位置连续、速度连续,即rCDO=rEDO,vCDO=vEDO。另外,太阳帆航天器运行在圆悬浮轨道上时,其瞬时速度和瞬时位置垂直。因此,考虑到上述位置和速度连续条件,拼接完成后,运行在椭圆悬浮轨道上的航天器在拼接点的瞬时位置矢量和瞬时速度矢量也必须垂直。综上可知,行星圆悬浮轨道和椭圆悬浮轨道的拼接只能发生在椭圆悬浮轨道的近焦点处或远焦点处,即fEDO=0,π。
在非一致旋转相合坐标系下建立太阳帆航天器的动力学模型为:
其中Ri是太阳帆相对于行星质心的位置矢量;R是太阳帆相对于坐标系原点的位置矢量;ω表示坐标系本身的旋转角速度;μp表示行星的引力常数;太阳光压加速度a可表示为:
a=ac(ns·n)·[2(1-u)ρs(ns·n)·n+(1-ρs+uρs)ns] (2)
其中,单位矢量ns表示太阳光方向;n表示太阳帆的法向;u为镜面反射部分面积占太阳帆总面积的比值;ρs表示照射在镜面反射上的太阳光完全反射的比率;ac表示行星附近的最大太阳光压加速度,即太阳帆帆面垂直太阳光方向。
为了方便描述行星椭圆悬浮轨道,假设一个半长轴a、偏心率e的参考轨道,该参考轨道位于行星引力场。假设行星椭圆悬浮轨道与该参考轨道具有相同的瞬时角速度和轨道周期T。因此,参考轨道的轨道方程s、轨道周期T,瞬时角速度和瞬时角加速度分别为:
以变化的距离S(f)为单位长度,因此,R=s·r,Ri=s·ri,将此连同坐标系角速度和方程(3)代入方程(1),得到太阳帆航天器在非一致旋转相合坐标系中的无量纲动力学方程为:
其中,r为航天器相对于主焦点F无量纲位置矢量,ri为航天器相对于行星质心的无量纲矢量,f表示轨道真近点角,e表示轨道偏心率,无量纲的太阳光压加速度aN可表示为:
通过强制r″=0和r′=0,方程(4)给出了太阳帆航天器在非一致旋转相合坐标系下保持静止的条件:
通过化简和坐标变化,可得其代数形式的平衡条件为:
该条件只包含了三个自由变量,因此,结合方程(2),在给定真近点处可计算太阳帆的两个姿态角和其反射率。并且太阳帆姿态约束和反射率参数u必须满足以下两个条件:c=n·ns∈(0,1],太阳帆的法向方向n必须执行太阳的反方向;u∈[0,umax],太阳帆的反射率不能超过其最大反射率。
为了得到旋转坐标系下的行星圆悬浮轨道,假设悬浮半径是ρ,悬浮高度为z,轨道角速度为ωc,太阳帆相应的倾斜角α和轻系数β分别为:
其中
将航天器正在运行的圆悬浮轨道作为转移的初始圆悬浮轨道,并根据任务要求,猜测转移的目标圆悬浮轨道。根据初始轨道和目标轨道的相对位置,可将其分为以下两种情况讨论:
共平面转移:如图4所示,两条行星圆悬浮轨道位于同一平面内,椭圆悬浮轨道的轨道参数分别记为ξ·a,e和L,其中a表示参考轨道的半长轴,ξ表示椭圆悬浮轨道半长轴和参考轨道半长轴之间的比值,e表示偏心率,L表示悬浮高度。因此
L=z1=z2 (10)
使用位置连续和速度连续条件,可以得到:
ρ1=ξa(1-e),ρ2=ξa(1+e) (11)
ρ1=ξa(1-e),ρ2=ξa(1+e) (12)
联合求解方程(11)和(12),可以得到:
由此可得到,用于实现两个圆悬浮轨道之间的椭圆悬浮轨道的轨道参数。
一般轨道转移:如图5所示,在情况一下,点SP1和点SP2的坐标可以表示为:
rSP1=[ρ1,0,z1],rSP2=[-ρ2,0,z2] (14)。
根据几何关系,可以求得椭圆悬浮轨道的悬浮高度为:
与上述情况相仿,点SP1和点SP2处速度连续条件和位置连续条件将被应用于求解椭圆悬浮轨道的轨道参数a,e和ξ:
更加一般的情况,如图6所示,在情况二下,椭圆悬浮轨道的参数求解除了计算点SP1和点SP2的坐标上有所区别外,其他计算与图5情况类似。
根据椭圆悬浮轨道的平衡条件(方程7),椭圆悬浮轨道的轨道参数和太阳光压加速度模型(方程2),计算能够实现作为转移轨道的椭圆悬浮轨道的太阳帆姿态参数和反射率参数。
在本例中,以共平面转移为例,假设任务在地心引力场中,且ac=6mm/s-2。初始地球圆悬浮轨道的轨道参数为:
ρ1=18RE,z1=30RE,Ω1=6×10-6rad/s-1
目标地球圆悬浮轨道的轨道参数为:
ρ2=22RE,z2=30RE,Ω2=4×10-6rad/s-1
通过设计,可得到作为椭圆悬浮轨道的轨道参数为:
a=40RE,e=0.1,ξ=0.5,L=30RE
为了实现初始圆悬浮轨道和目标圆悬浮轨道之间的转移,经过设计,得到了作为转移轨道的椭圆悬浮轨道。图7-9描述了在地球引力场中实现该椭圆悬浮轨道需要的太阳帆姿态变化和反射率参数变化。

Claims (7)

1.一种太阳帆航天器行星圆悬浮轨道之间的转移方法,其特征在于,包括以下步骤:
步骤1,确定太阳帆航天器行星椭圆悬浮轨道模型和行星圆悬浮轨道模型,并且确定两类轨道下太阳帆的唯一描述参数;
步骤2,通过行星圆悬浮轨道上的太阳帆航天器的位置矢量和速度矢量,确定行星椭圆悬浮轨道和行星圆悬浮轨道进行拼接的边界条件;
步骤3,将圆悬浮轨道之间的转移作为目标,得到实现行星圆悬浮轨道之间转移的椭圆悬浮轨道;
步骤4,通过边界条件确定椭圆悬浮轨道的轨道参数,根据椭圆悬浮轨道的平衡条件和太阳光压加速度模型,得到实现椭圆悬浮转移轨道的太阳帆姿态角参数和反射率参数;
步骤5,判定太阳帆姿态参数和反射率参数是否可实现,若能实现,则椭圆悬浮轨道可作为转移轨道;否则重复步骤1-4;
步骤6,完成太阳帆航天器在行星圆悬浮轨道之间的转移。
2.根据权利要求1所述的太阳帆航天器行星圆悬浮轨道之间的转移方法,其特征在于,所述步骤1中通过建立惯性参考坐标系非一致旋转相合坐标系和旋转坐标系确定两类轨道下太阳帆的唯一描述参数。
3.根据权利要求1所述的太阳帆航天器行星圆悬浮轨道之间的转移方法,其特征在于,所述步骤2中行星椭圆悬浮轨道和行星圆悬浮轨道进行拼接的位置为椭圆悬浮轨道的近焦点或远焦点。
4.根据权利要求1所述的太阳帆航天器行星圆悬浮轨道之间的转移方法,其特征在于,所述步骤3中两条行星圆悬浮轨道位于同一平面内时,得到实现行星圆悬浮轨道之间转移的椭圆悬浮轨道。
5.根据权利要求1所述的太阳帆航天器行星圆悬浮轨道之间的转移方法,其特征在于,所述步骤3两条行星圆悬浮轨道不在同一平面内时,得到实现行星圆悬浮轨道之间转移的椭圆悬浮轨道。
6.根据权利要求1所述的太阳帆航天器行星圆悬浮轨道之间的转移方法,其特征在于,所述步骤4中在非一致旋转相合坐标系下建立太阳帆航天器的动力学模型其中Ri是太阳帆相对于行星质心的位置矢量;R是太阳帆相对于坐标系原点的位置矢量;ω表示坐标系本身的旋转角速度;μp表示行星的引力常数;太阳光压加速度a表示为:a=ac(ns·n)·[2(1-u)ρs(ns·n)·n+(1-ρs+uρs)ns],单位矢量ns表示太阳光方向;n表示太阳帆的法向;u为镜面反射部分面积占太阳帆总面积的比值;ρs表示照射在镜面反射上的太阳光完全反射的比率;ac表示行星附近的最大太阳光压加速度。
7.根据权利要求1所述的太阳帆航天器行星圆悬浮轨道之间的转移方法,其特征在于,所述步骤5中太阳帆姿态约束c(α,θ)和反射率参数u满足以下条件为太阳帆姿态参数和反射率参数能够实现:c=n·ns∈(0,1]太阳帆的法向方向n必须指向太阳的反方向;u∈[0,umax]太阳帆的反射率不能超过其最大反射率。
CN201710369368.9A 2017-05-23 2017-05-23 一种太阳帆航天器行星圆悬浮轨道之间的转移方法 Expired - Fee Related CN107310752B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710369368.9A CN107310752B (zh) 2017-05-23 2017-05-23 一种太阳帆航天器行星圆悬浮轨道之间的转移方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710369368.9A CN107310752B (zh) 2017-05-23 2017-05-23 一种太阳帆航天器行星圆悬浮轨道之间的转移方法

Publications (2)

Publication Number Publication Date
CN107310752A true CN107310752A (zh) 2017-11-03
CN107310752B CN107310752B (zh) 2020-09-08

Family

ID=60181885

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710369368.9A Expired - Fee Related CN107310752B (zh) 2017-05-23 2017-05-23 一种太阳帆航天器行星圆悬浮轨道之间的转移方法

Country Status (1)

Country Link
CN (1) CN107310752B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108303874A (zh) * 2017-12-28 2018-07-20 北京航空航天大学 一种针对绳系空间拖船系统摆振的小推力切换控制方法
CN108860662A (zh) * 2018-06-24 2018-11-23 中国人民解放军战略支援部队航天工程大学 一种基于太阳光压的中轨道卫星离轨方法
CN109032176A (zh) * 2018-07-25 2018-12-18 西北工业大学 一种基于微分代数的地球同步轨道确定和参数确定方法
CN109190155A (zh) * 2018-07-25 2019-01-11 西北工业大学 一种采用电推进/太阳帆推进的混合连续小推力轨道设计方法
CN109213190A (zh) * 2018-07-23 2019-01-15 西北工业大学 一种日心悬浮轨道上电动帆航天器编队飞行协同控制方法
CN109299493A (zh) * 2018-07-23 2019-02-01 西北工业大学 一种日心悬浮轨道上太阳帆航天器相对运动可达域确定方法
CN114013676A (zh) * 2022-01-10 2022-02-08 南充职业技术学院 一种城市规划测绘用无人机起飞云台

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106096204A (zh) * 2016-06-28 2016-11-09 西北工业大学 一种基于太阳帆推进技术的航天器日心椭圆悬浮轨道设计方法
CN106200376A (zh) * 2016-06-28 2016-12-07 西北工业大学 一种航天器日心非开普勒悬浮轨道拼接方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106096204A (zh) * 2016-06-28 2016-11-09 西北工业大学 一种基于太阳帆推进技术的航天器日心椭圆悬浮轨道设计方法
CN106200376A (zh) * 2016-06-28 2016-12-07 西北工业大学 一种航天器日心非开普勒悬浮轨道拼接方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
钱航: "太阳帆航天器轨道和姿态耦合设计与优化", 《中国博士学位论文全文数据库-工程科技II辑》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108303874A (zh) * 2017-12-28 2018-07-20 北京航空航天大学 一种针对绳系空间拖船系统摆振的小推力切换控制方法
CN108860662A (zh) * 2018-06-24 2018-11-23 中国人民解放军战略支援部队航天工程大学 一种基于太阳光压的中轨道卫星离轨方法
CN109213190A (zh) * 2018-07-23 2019-01-15 西北工业大学 一种日心悬浮轨道上电动帆航天器编队飞行协同控制方法
CN109299493A (zh) * 2018-07-23 2019-02-01 西北工业大学 一种日心悬浮轨道上太阳帆航天器相对运动可达域确定方法
CN109032176A (zh) * 2018-07-25 2018-12-18 西北工业大学 一种基于微分代数的地球同步轨道确定和参数确定方法
CN109190155A (zh) * 2018-07-25 2019-01-11 西北工业大学 一种采用电推进/太阳帆推进的混合连续小推力轨道设计方法
CN114013676A (zh) * 2022-01-10 2022-02-08 南充职业技术学院 一种城市规划测绘用无人机起飞云台
CN114013676B (zh) * 2022-01-10 2022-03-15 南充职业技术学院 一种城市规划测绘用无人机起飞云台

Also Published As

Publication number Publication date
CN107310752B (zh) 2020-09-08

Similar Documents

Publication Publication Date Title
CN107310752A (zh) 一种太阳帆航天器行星圆悬浮轨道之间的转移方法
CN102880193B (zh) 航天器开发试验台系统
CN106096204B (zh) 一种基于太阳帆推进技术的航天器日心椭圆悬浮轨道设计方法
Bassetto et al. Sliding mode control for attitude maneuvers of Helianthus solar sail
CN104090612A (zh) 一种基于偏航导引的倾斜轨道航天器获取能源的方法
CN106200376B (zh) 一种航天器日心非开普勒悬浮轨道拼接方法
CN107247825A (zh) 一种基于太阳帆航天器的行星椭圆悬浮轨道设计方法
CN105974822A (zh) 一种航天器自主绕飞交会控制系统验证装置及其验证方法
CN105253329B (zh) 一种基于弱稳定边界的两脉冲行星捕获轨道方法
CN107450309B (zh) 一种绳系组合体姿态机动协同控制方法
CN106125745A (zh) 一种对空间合作目标跟踪成像的卫星姿态控制方法
Adeli et al. A scalable bus-based attitude control system for Solar Sails
CN105631099A (zh) 一种小天体探测器着陆动力学模拟系统
CN107589756A (zh) 一种奔月卫星编队初始化方法
CN101186236A (zh) 一种减少航天器重力损失的变轨方法
CN106679674A (zh) 基于星历模型的地月L2点Halo轨道阴影分析方法
CN113641182B (zh) 星间激光通信系统的高精度瞄准指向方法及系统
CN104123461B (zh) 一种用于空间物体光度分析的光照可视关系计算方法
CN112061424B (zh) 一种基于融合目标姿态的机动过程能源角动态跟踪方法
CN113310496A (zh) 一种确定月地转移轨道的方法及装置
CN105035371A (zh) 一种基于osg三维引擎的经典轨道三维空间关系构建方法
CN108082538A (zh) 一种考虑始末约束的多体系统低能量捕获轨道方法
CN105301958B (zh) 一种基于气动力辅助的平衡点周期轨道捕获方法
CN112945242B (zh) 一种卫星在轨自主规划任务最佳时间与姿态的方法
Calaon et al. Momentum management of a spacecraft equipped with a dual-gimballed electric thruster

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200908