行星借力小推力转移轨道初始设计研究.pdfVIP

果，而对于远行星探测任务，该方法计算得到的 其中：F = mr / r3 转移轨道与最优转移相差较大；另一方面，该方 由控制方程可得最优控制律为： H 法在行星借力处采用了能量匹配策略，这就不可 1 0 u ? (6) 避免的忽略了一些合适的转移机会。 u v 针对以上问题，本节将给出一种基于参考轨 式( 1)~( 6) 给出了小推力轨道两点边值问题的 求解模型。由于该问题不存在解析解，目前主要 道的行星际小推力转移轨道设计方法。该方法利 采用数值解法。然而，由于共轭方程非线性很强， 用参考轨道法计算单推力弧段的两点边值问题， 并且共轭变量没有实际物理意义，很难给出较为 采用有推力借力匹配策略完成轨道的拼接，该方 准确的初值猜测，从而导致数值算法收敛性差， 法具有较好的收敛性，并且相比 Shape-Based 方 计算效率较低。 法，可以搜索到更多可行的转移轨道。 1.2 模型线性化处理与求解 考虑方程( 1) ，如果推进器不工作，即u(t) 0 ， 1 基于参考轨道的两点边值问题 则飞行器的飞行轨迹为开普勒轨道。此时飞行器 满足如下运动方程： T T T T T 3 T 1.1 两点边值问题 z [R ,V ] [V , m R / R ] (7) T T T 在日心惯性坐标系中，采用笛卡尔坐标描述 其中：z = [R ,V ] 为飞行器在开普勒轨道上的状 的飞行器运动方程为 态矢量。 对于开普勒轨道，若已知飞行器从地球发射 r v 3 (1) v m r / r u 的时间t 和到达目标星的时间t ，通过星历计算 0 f 其中：r , v 分别为飞行器的位置和速度矢量，m 为 可得飞行器从地球出发的位置矢量R 和到达目 太阳引力常数，u 为推进器推力加速度。