CN109190155B - 一种采用电推进/太阳帆推进的混合连续小推力轨道设计方法 - Google Patents
一种采用电推进/太阳帆推进的混合连续小推力轨道设计方法 Download PDFInfo
- Publication number
- CN109190155B CN109190155B CN201810828084.6A CN201810828084A CN109190155B CN 109190155 B CN109190155 B CN 109190155B CN 201810828084 A CN201810828084 A CN 201810828084A CN 109190155 B CN109190155 B CN 109190155B
- Authority
- CN
- China
- Prior art keywords
- thrust
- gravitational field
- virtual
- orbit
- propulsion
- 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
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
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)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种采用电推进/太阳帆推进的混合连续小推力轨道设计方法,首先建立虚拟引力场动力学模型;其次建立在虚拟引力场中的切向推力轨道动力学模型,求解其解析解;然后对所求得的解析解进行精度分析;最后建立轨道优化模型,采用电推进/太阳帆推进的混合连续小推力方法进行轨道优化。本发明在保证精度的条件下,能够大大的减小计算量,能够减小实现轨道转移所需的推力,扩大了小推力轨道应用的范围。
Description
技术领域
本发明属于航天技术开发领域,具体涉及一种采用电推进/太阳帆推进的混合连续小推力轨道设计方法。
背景技术
连续小推力轨道设计在空间机动轨道设计方面有广泛的应用。与传统脉冲推力或大推力相比,连续小推力轨道机动容易控制,机动能力强;另外其比冲小,实现轨道机动所需要的燃料很少,因而成为研究热点。Mcinnes首次提出了使用太阳帆来改变引力场的大小,设计了悬浮轨道,并推导了实现悬浮轨道的所需推力的解析解。进一步,Malcolm等采用太阳帆产生切向固定推力,设计了从地球轨道到太阳极轨的转移轨道,并进一步解决了轨道优化问题。然而,由于单独采用电推进或者太阳帆产生的推力较小,轨道转移常常需要很长时间。针对这个问题,Mengali和Quarta提出了太阳帆和太阳帆电推进结合的混合推进方法,并采用间接优化方法进行了轨道优化。其中,太阳帆产生径向推力,用于改变引力场的引力常数,而于此同时电推进产生固定切向推力。在假定太阳帆能够产生50%的中心引力(μ2=0.5μ1)的条件下,混合推力能够减少燃料的消耗并且缩短空间机动任务时间。然而该方法是作为一种间接优化方法,并没有考虑推力约束问题,很难应用于实际空间机动轨道设计中。
发明内容
本发明的目的在于提供一种采用电推进/太阳帆推进的混合连续小推力轨道设计方法,以克服上述现有技术存在的缺陷,本发明在保证精度的条件下,能够大大的减小计算量,能够减小实现轨道转移所需的推力,扩大了小推力轨道应用的范围。
为达到上述目的,本发明采用如下技术方案:
一种采用电推进/太阳帆推进的混合连续小推力轨道设计方法,包括如下步骤:
1)建立虚拟引力场动力学模型;
2)建立在虚拟引力场中的切向推力轨道动力学模型,求解其解析解;
3)对所求得的解析解进行精度分析;
4)采用电推进/太阳帆推进的混合连续小推力方法进行轨道优化。
进一步地,步骤1)具体为:
其中,v1为航天器在实际引力场中的速度,v2为航天器在虚拟引力场中的速度,由式(2),式(3),在虚拟中心引力场坐标系下,航天器运动方程表示为如下形式:
其中μ2为虚拟引力场常数,实现虚拟中心引力场所需推力由下式得到,
即通过调整推力,能够形成一种虚拟中心引力场,在该虚拟中心引力场中,航天器运动轨迹为虚拟圆锥曲线,若该虚拟圆锥曲线轨道满足轨道转移约束,则能够实现轨道机动。
进一步地,步骤2)具体为:
在切向推力作用下,平面轨道高斯摄动方程为,
其中,h为轨道角动量,r为轨道高度,fr为径向推力;a为轨道半长轴,e为偏心率,ω为近地点俯角,θ为真近点角,E为偏近点角,p=a(1-e2)为半通径,为瞬时时刻的平均轨道角速度,μ为引力场数,假定径向加速度为零,则有:
径向推力fr以及周向推力fθ和切向推力f之间的关系为:
而真近点角和偏近点角的关系为:
因此,将式(7)中自变量转换为偏近点角,则轨道半长轴和偏心率的变化为:
其中E为偏近点角,E0、Ef分别为初始时刻和终点时刻的偏近点角,ft为连续切向推力;
在虚拟引力场中的切向推力轨道近似解析解为:
avg、Δavg、μvg、Δevg分别为虚拟引力场轨道半长轴,虚拟引力场轨道半长轴增量,虚拟引力场数,虚拟引力场轨道偏心率增量;当初始点和目标点已知时,上式中的参数E已知,边界约束则转化为:
其中(ai)vg为虚拟引力场中的初始轨道半长轴,(at)vg为虚拟引力场中的目标轨道半长轴,(Δa)vg为虚拟引力场中轨道机动产生的半长轴增量,其中(ei)vg为虚拟引力场中的初始轨道偏心率,(et)vg为虚拟引力场中的目标轨道偏心率,(Δa)vg为虚拟引力场中轨道机动产生的偏心率增量;由此,采用虚拟引力场下的连续推力轨道设计方法,能够将所有满足边界约束的轨道用一组参数x=(μvg,rvg,ft)表示,即转移轨道参数化,rvg为虚拟引力场的引力中心偏移量,这样转移轨道设计与优化问题就转化成为参数设计与优化问题,采用虚拟引力场下的混合连续小推力轨道设计,需要给航天器施加推力,一部分推力用于形成虚拟引力场,另一部分推力用于形成虚拟引力场中的切向力,具体推力分配:利用太阳帆形成虚拟引力场,利用电推进来产生虚拟引力场中的切向推力;当太阳帆和电推进同时作用时,实际的轨道相当于航天器在虚拟引力场中受到固定切向推力作用的机动轨道。
进一步地,步骤3)具体为:
为描述解析解与数值解的误差,在此定义误差因子参数,如式(14)所示:
其中Δθ为轨道转移的转移角,通过对比式(10)及式(12),解析解的误差因子参数会随着切向推力、轨道半长轴、偏近点角的增大而增大,当e<0.7时,偏心率越大,误差因子参数越大;
在采用虚拟引力场常值切向力解析方法设计转移轨道时,只需要保证误差因子k<k*,此处k*为容许误差因子。
进一步地,步骤4)具体为:
假定混合径向推力加速度为Tr,太阳帆和电推进分别产生的为assr,aSEPr;切向推力为Tc,太阳帆和电推进分别产生的为assc,aSEPc,太阳帆推进和电推进的数学模型如下所示:
在此假定太阳帆推进的姿态角α是真近点角θ的二阶多项式,假定多项式系数为ai,i=0,1,2,如下所示:
α=a2·θ2+a1·θ+a0
通过对多项式的参数进行优化,能够最大限度地增大太阳帆的使用效率,减小电推进的燃料消耗,进而减小混合推力的燃料消耗,电推进所需要提供的推力加速度函数Fr(·),Fr(·)表示为:
假定虚拟引力场参数上下限为μvgL,μvgU,切向力参数上下限为fsepL,fsepU,引力场中心位置偏移量的上下限为rvgL,rvgU,则有:
μvg∈[μvgL,μvgU],rvg∈[rvgL,rvgU],fsep∈[fsepL,fsepU]
由此可见,轨道优化问题即转化为参数优化问题,优化的目标函数为转移轨道消耗燃料质量:
J=min(mfuel)=max(massf) (15)
其中mfuel为燃料消耗量,massf为任务结束时航天器的质量。
与现有技术相比,本发明具有以下有益的技术效果:
本发明方法能够以解析的形式求出大量的满足轨道边界约束的转移轨道,并求解出转移轨道的解析解。同时,在给定目标函数的条件下,可对满足推力约束的转移轨道进行优化。与其它的小推力轨道设计方法相比,该方法有以下的特点:1、该轨道设计方法是一种解析方法,在保证精度的条件下,能够大大的减小计算量,因而该方法可以为数值优化方法(直接法或间接法)提供轨道优化初值;2、与常值切向推力相比,该方法能够减小实现轨道转移所需的推力,扩大了小推力轨道应用的范围;3、与常值切向推力相比,该方法不仅能够减小轨道转移所消耗的能量,而且能够减小轨道转移的飞行时间。
附图说明
图1为虚拟中心引力场示意图;
图2为混合小推力轨道交会示意图;
图3为混合小推力转移轨道优化示意图;
图4为速度误差、位置误差随误差因子的变化。
具体实施方式
下面对本发明作进一步详细描述:
本发明一种采用电推进/太阳帆推进的混合连续小推力轨道设计方法,包含以下步骤:a、该方法结合轨道平均方法和虚拟引力势场方法,将小推力作用下的机动轨道参数化,并能够求解出实现轨道机动的解析解;b、在分析了轨道误差特性的基础上,确定轨道解析解的误差因子,在保证满足误差因子约束的条件下,得到高精度的解析解;基于上述设计方法的优化算法,包含以下步骤:a1、采用形状方法,将太阳帆的姿态角表示为转移角的多项式形式;b1、给定以能量消耗为目标函数,对该多项式的系数进行优化,以减小混合推进推进器的燃料质量。
具体步骤如下:
1、建立虚拟引力场动力学模型,求解动力学解析解;
其中,v1为航天器在实际引力场中的速度,v2为航天器在虚拟引力场中的速度。由式(2),式(3)可知,在虚拟中心引力场坐标系下,航天器运动方程可表示为如下形式:
其中μ2为虚拟引力场常数,实现虚拟中心引力场所需推力可由下式得到,
由上述分析可知,通过调整推力,可形成一种虚拟中心引力场。在该虚拟中心引力场中,航天器运动轨迹为虚拟圆锥曲线。若该虚拟圆锥曲线轨道满足轨道转移约束,则可实现轨道机动。
2、建立在虚拟引力场中的切向推力轨道动力学模型,求解其解析解;
在切向推力作用下,平面轨道高斯摄动方程为,
其中,h为轨道角动量,r为轨道高度,fr为径向推力;a为轨道半长轴,e为偏心率,ω为近地点俯角,θ为真近点角,E为偏近点角,p=a(1-e2)为半通径,为瞬时时刻的平均轨道角速度,μ为引力场数。假定径向加速度为零,则有
径向推力fr以及周向推力fθ和切向推力f之间的关系为:
而真近点角和偏近点角的关系为,
因此,将公式(7)中自变量转换为偏近点角,则轨道半长轴和偏心率的变化为,
其中E为偏近点角,E0 Ef分别为初始时刻和终点时刻的偏近点角,ft为连续切向推力,c为常数;
在给定轨道交会的初始点A和终点B的条件下,采用混合推力方法设计机动轨道就是求解引力场参数(μ2,rvg)以及切向推力fvg。如图2所示,在虚拟引力场下,航天器从A点可以转移到B’点,而若选择合适参数x=(μ2,rvg,fvg)时,在虚拟引力场中和切向推力的同时作用下,航天器在虚拟引力场可从A点转移到目标点B点,从而满足轨道交会的边界约束。这就是采用基于虚拟引力场方法的混合连续小推力轨道设计的基本思路。
采用虚拟引力场方法的混合连续小推力轨道设计方法来设计轨道,实际上是设计在虚拟引力场中的常值切向推力转移轨道。转移轨道需要满足两点边界值约束(位置约束和速度约束),以及推力约束。当初始轨道的偏心率很小,推力较小的条件下,在虚拟引力场中的切向推力轨道近似解析解为:
avgΔavgμvgΔevg分别为虚拟引力场轨道半长轴,虚拟引力场轨道半长轴增量,虚拟引力场数,虚拟引力场轨道偏心率增量。当初始点和目标点已知时,上式中的参数E已知。边界约束可以转化为,
其中(ai)vg为虚拟引力场中的初始轨道半长轴,(at)vg为虚拟引力场中的目标轨道半长轴,(Δa)vg为虚拟引力场中轨道机动产生的半长轴增量。其中(ei)vg为虚拟引力场中的初始轨道偏心率,(et)vg为虚拟引力场中的目标轨道偏心率,(Δa)vg为虚拟引力场中轨道机动产生的偏心率增量。由此可见,采用虚拟引力场下的连续推力轨道设计方法,可将所有满足边界约束的轨道用一组参数x=(μvg,rvg,ft)表示,即转移轨道参数化。(此处rvg为虚拟引力场的引力中心偏移量)这样转移轨道设计与优化问题就可以转化成为参数设计与优化问题。
由上分析可知,采用虚拟引力场下的混合连续小推力轨道设计方法,需要给航天器施加推力,一部分推力用于形成虚拟引力场,另一部分推力用于形成虚拟引力场中的切向力。具体推力分配:利用太阳帆形成虚拟引力场,利用电推进来产生虚拟引力场中的切向推力。当太阳帆和电推进同时作用时,实际的轨道相当于航天器在虚拟引力场中受到固定切向推力作用的机动轨道。
3、解析解精度分析;
为描述解析解与数值解的误差,在此定义一个误差因子参数,如公式(14)所示,
其中Δθ为轨道转移的转移角,对比公式(10)及(12)可知,解析解的误差因子会随着切向推力、轨道半长轴、轨道转移时间(偏近点角)的增大而增大,当(e<0.7)时,偏心率越大,误差参数越大。
如图4所示,采用解析方法可以得到大量的转移轨道,同时采用四阶隆格库塔数值方法,就可以求解并分析这些解析解的速度误差、位置误差以及误差因子之间的关系。
以地球-火星转移轨道为例,采用解析方法可以得到大量的转移轨道,在其中任意选择n条转移轨道(n=157)。采用数值积分方法,可以得出这些解析解的速度误差和位置误差与误差因子的关系,如图4所示。由图4可知,速度误差、位置误差和误差因子是同步的,当误差因子很小时,其对应的速度误差和位置误差也很小;相反,当误差因子比较大时,其速度误差和位置误差都很大,可以看出速度误差和位置误差都随着误差因子的增大而增大。当保证误差因子小于0.0025时,速度误差小于0.001,而位置误差小于0.002。由以上分析可知,在采用虚拟引力场常值切向力解析方法设计转移轨道时,只需要保证误差因子k<k*,此处k*为容许误差因子,就可以得到满足相应的精度要求的转移轨道。
4、建立轨道优化模型,采用混合小推力方法进行轨道优化
采用太阳帆/电推进混合推力,其推力由两部分组成,第一部分为太阳帆产生的推进,第二部分为电推进产生,假定混合径向推力加速度为Tr,太阳帆和电推进分别产生的为assr,aSEPr;切向推力为Tc,太阳帆和电推进分别产生的为assc,aSEPc;如下所示:
为太阳帆推进和电推进的数学模型。在此假定太阳帆推进的姿态角α是真近点角θ的二阶多项式,假定多项式系数为ai(i=0,1,2),如下所示,
α=a2·θ2+a1·θ+a0
通过对多项式的参数进行优化,可以最大限度地增大太阳帆的使用效率,这样就减小了电推进的燃料消耗,进而减小了混合推力的燃料消耗。因此在整个任务中,电推进所需要提供的推力加速度函数Fr(·),Fr(·)可以表示为:
采用基于混合连续小推力轨道设计方法,可求出大量的满足边界约束的转移轨道,如图3所示。假定虚拟引力场参数上下限μvgL,μvgU,切向力参数上下限为fsepL,fsepU,引力场中心位置偏移量的上下限位rvgL,rvgU,则有:
μvg∈[μvgL,μvgU],rvg∈[rvgL,rvgU],fsep∈[fsepL,fsepU]
由此可见,轨道优化问题即可转化为参数优化问题,优化的目标函数可以为转移轨道消耗燃料质量:
J=min(mfuel)=max(massf) (15)
其中mfuel为燃料消耗量,massf为任务结束时航天器的质量。
综上所述,本发明提出一种基于混合连续小推力的轨道设计与优化方法,该方法能够给出连续小推力的解析解,并且求出大量的满足约束条件的轨道。通过对太阳帆姿态角的优化,可以得到最节省燃料的转移轨道。与传统的单独电推进或者单独太阳帆推进相比,混合推进方法能够降低实现轨道机动的推力大小,在节省燃料消耗的同时缩短飞行时间。该方法能够很大程度上扩展小推力方法在轨道设计中的应用范围,为小推力轨道设计提出了另外一种设计思路。最后以地球-火星轨道转移为算例,对比了单纯使用太阳帆或者单独使用电推进的方法,验证了混合推力轨道的可行性,以及在燃料消耗、任务时间上的优越性,对比结果如表1和表2所示。
表1地球火星轨道转移虚拟引力场方法的优化结果
表2混合推力与太阳帆推进、电推进的结果比较
由表2所示的推力方法对比可知,无论从飞行时间,还是推力大小,本发明提出的方法均优于常值推力方法。
Claims (3)
1.一种采用电推进/太阳帆推进的混合连续小推力轨道设计方法,其特征在于,包括如下步骤:
1)建立虚拟引力场动力学模型;
具体为:
其中μ2为虚拟引力场常数,实现虚拟中心引力场所需推力由下式得到,
即通过调整推力,能够形成一种虚拟中心引力场,在该虚拟中心引力场中,航天器运动轨迹为虚拟圆锥曲线,若该虚拟圆锥曲线轨道满足轨道转移约束,则能够实现轨道机动;
2)建立在虚拟引力场中的切向推力轨道动力学模型,求解其解析解;
具体为:
在切向推力作用下,平面轨道高斯摄动方程为,
其中,h为轨道角动量,r为轨道高度,fr为径向推力;a为轨道半长轴,e为偏心率,ω为近地点俯角,θ为真近点角,E为偏近点角,p=a(1-e2)为半通径,为瞬时时刻的平均轨道角速度,μ为引力场数,假定径向加速度为零,则有:
径向推力fr以及周向推力fθ和切向推力f之间的关系为:
而真近点角和偏近点角的关系为:
因此,将式(7)中自变量转换为偏近点角,则轨道半长轴和偏心率的变化为:
其中E为偏近点角,E0、Ef分别为初始时刻和终点时刻的偏近点角,ft为连续切向推力;
在虚拟引力场中的切向推力轨道近似解析解为:
avg、Δavg、μvg、Δevg分别为虚拟引力场轨道半长轴,虚拟引力场轨道半长轴增量,虚拟引力场数,虚拟引力场轨道偏心率增量;当初始点和目标点已知时,上式中的参数E已知,边界约束则转化为:
其中(ai)vg为虚拟引力场中的初始轨道半长轴,(at)vg为虚拟引力场中的目标轨道半长轴,(Δa)vg为虚拟引力场中轨道机动产生的半长轴增量,其中(ei)vg为虚拟引力场中的初始轨道偏心率,(et)vg为虚拟引力场中的目标轨道偏心率,(Δe)vg为虚拟引力场中轨道机动产生的偏心率增量;由此,采用虚拟引力场下的连续推力轨道设计方法,能够将所有满足边界约束的轨道用一组参数x=(μvg,rvg,ft)表示,即转移轨道参数化,rvg为虚拟引力场的引力中心偏移量,这样转移轨道设计与优化问题就转化成为参数设计与优化问题,采用虚拟引力场下的混合连续小推力轨道设计,需要给航天器施加推力,一部分推力用于形成虚拟引力场,另一部分推力用于形成虚拟引力场中的切向力,具体推力分配:利用太阳帆形成虚拟引力场,利用电推进来产生虚拟引力场中的切向推力;当太阳帆和电推进同时作用时,实际的轨道相当于航天器在虚拟引力场中受到固定切向推力作用的机动轨道;
3)对所求得的解析解进行精度分析;
4)采用电推进/太阳帆推进的混合连续小推力方法进行轨道优化。
3.根据权利要求2所述的一种采用电推进/太阳帆推进的混合连续小推力轨道设计方法,其特征在于,步骤4)具体为:
假定混合径向推力加速度为Tr,太阳帆和电推进分别产生的为assr,aSEPr;切向推力加速度为Tc,太阳帆和电推进分别产生的为assc,aSEPc,太阳帆推进和电推进的数学模型如下所示:
在此假定太阳帆推进的姿态角α是真近点角θ的二阶多项式,假定多项式系数为ai,i=0,1,2,如下所示:
α=a2·θ2+a1·θ+a0
通过对多项式的参数进行优化,能够最大限度地增大太阳帆的使用效率,减小电推进的燃料消耗,进而减小混合推力的燃料消耗,电推进所需要提供的径向推力加速度函数Fr(·),切向加速度函数Ft(·)分别表示为:
假定虚拟引力场参数上下限为μvgL,μvgU,切向力参数上下限为fsepL,fsepU,引力场中心位置偏移量的上下限为rvgL,rvgU,则有:
μvg∈[μvgL,μvgU],rvg∈[rvgL,rvgU],fsep∈[fsepL,fsepU]
由此可见,轨道优化问题即转化为参数优化问题,优化的目标函数为转移轨道消耗燃料质量:
J=min(mfuel)=max(massf) (15)
其中mfuel为燃料消耗量,massf为任务结束时航天器的质量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810828084.6A CN109190155B (zh) | 2018-07-25 | 2018-07-25 | 一种采用电推进/太阳帆推进的混合连续小推力轨道设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810828084.6A CN109190155B (zh) | 2018-07-25 | 2018-07-25 | 一种采用电推进/太阳帆推进的混合连续小推力轨道设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109190155A CN109190155A (zh) | 2019-01-11 |
CN109190155B true CN109190155B (zh) | 2022-09-16 |
Family
ID=64937357
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810828084.6A Active CN109190155B (zh) | 2018-07-25 | 2018-07-25 | 一种采用电推进/太阳帆推进的混合连续小推力轨道设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109190155B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111324137B (zh) * | 2020-03-05 | 2022-08-05 | 中国工程物理研究院总体工程研究所 | 一种基于常值连续推力的区域悬停轨道控制方法 |
CN112393835B (zh) * | 2020-11-03 | 2021-12-21 | 西北工业大学深圳研究院 | 一种基于扩展卡尔曼滤波的小卫星在轨推力标定方法 |
CN113602534B (zh) * | 2021-06-26 | 2023-02-28 | 山东航天电子技术研究所 | 一种微型电推进推力大小的在轨标定方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7913953B2 (en) * | 2005-12-28 | 2011-03-29 | Frank Werner Ellinghaus | Solar sail launch system and solar sail attitude control system |
CN106202640B (zh) * | 2016-06-28 | 2018-03-27 | 西北工业大学 | 日‑地三体引力场中的晕轨道航天器偏置轨道设计方法 |
CN107310752B (zh) * | 2017-05-23 | 2020-09-08 | 西北工业大学 | 一种太阳帆航天器行星圆悬浮轨道之间的转移方法 |
-
2018
- 2018-07-25 CN CN201810828084.6A patent/CN109190155B/zh active Active
Non-Patent Citations (2)
Title |
---|
Displaced Geostationary Orbit Design Using Hybrid Sail Propulsion;Jeannette Heiligers et al;《Guidance, Control, and Dynamics》;20111130;全文 * |
混合推力作用下基于傅里叶级数展开的轨道设计方法;方群 等;《西北工业大学学报》;20171031;第35卷(第5期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109190155A (zh) | 2019-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110794863B (zh) | 一种控制性能指标可定制的重型运载火箭姿态控制方法 | |
CN109190155B (zh) | 一种采用电推进/太阳帆推进的混合连续小推力轨道设计方法 | |
CN109911249B (zh) | 低推重比飞行器的星际转移有限推力入轨迭代制导方法 | |
CN110989644B (zh) | 一种考虑目标点多终端约束的飞行器轨迹规划方法 | |
CN110733670B (zh) | 一种短航程低过载的再入轨迹设计方法 | |
CN112572835B (zh) | 一种具有姿态切换的卫星在轨角动量管理及控制方法 | |
CN109613928B (zh) | 一种用于多矢量螺旋桨组合浮空器的复合控制系统及方法 | |
CN112061424B (zh) | 一种基于融合目标姿态的机动过程能源角动态跟踪方法 | |
CN108536009B (zh) | 一种混合小推力航天器的轨道保持和扰动抑制方法 | |
CN114637312A (zh) | 一种基于智能变形决策的无人机节能飞行控制方法及系统 | |
Wang et al. | Six-DOF trajectory optimization for reusable launch vehicles via Gauss pseudospectral method | |
CN113581494A (zh) | 一种geo卫星化电混合推进变轨方法 | |
Yu et al. | Analytical solutions to aeroassisted orbital transfer problem | |
CN114265420B (zh) | 适于高动态、控制慢响应的制导控制一体化设计方法 | |
CN113830332B (zh) | 一种电推轨道转移的点火姿态建立与动态跟踪方法 | |
CN113343561B (zh) | 一种求解最优月球飞掠转移轨道的方法及系统 | |
CN115993777A (zh) | 基于轨道摄动模型反演的径切联控解耦迭代标定方法 | |
CN108984934B (zh) | 一种基于周向或切向混合连续小推力的轨道设计方法 | |
CN113741551A (zh) | 一种基于代理模型的全过程轨迹优化方法及装置 | |
Lee | Mission and trajectory optimization of the air-launching rocket system using MDO techniques | |
Cheng et al. | Cross-cycle iterative unmanned aerial vehicle reentry guidance based on reinforcement learning | |
Aprovitola et al. | Optimal design of a next generation high-lift reusable re-entry vehicle | |
CN113968360B (zh) | 一种静止轨道卫星星上自主电推进变轨方法 | |
Dai et al. | Trajectory Optimization and Tracking Control of a Morphing Aircraft | |
Shi et al. | A pseudospectral approach to ascent trajectory optimization for hypersonic air-breathing vehicles |
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 |