基于混合法的小推力轨道优化.pdfVIP

基于混合法的小推力轨道优化 尚海滨崔平远栾恩杰 (哈尔滨工业大学深空探测基础研究中心，哈尔滨150001) 摘要针对小推力轨道优化问题，本文给出了一种混合优化算法，混合优化算法结合了庞德里 亚全最大值原理和非线性规划的优点，首先利用最大值原理特问题转化为一个两点边值问题，然后 采用非线性规划方法进行求解。最后利用这种方法优化了星际问的小推力转移轨道，数值计算表 明了该方法的有效性。 关键词 小推力，最太值原理，非线性规划，混合算法。 引言 随着小推力发动机技术的发展，小推力轨道的设计与优化已成为研究的热点。目前对于小推力轨道优 化主要分为两类方法：直接法和间接法。直接法是直接对状态和控制进行离散化后，将寻优问题转化为参数 优化的同题。相对于问接法，直接法的初值猜测较为简单，并且问题的收敛性好。但由于没有引入共轭方 程，直接法得到的控制量仅是依靠某种插值方法得到的，不一定符合实际情况。为保证计算的精度，就要求 取较多的参数，即离散化时采用的时间步长比较小。这样就会增加优化是的计算量，而且如果插值方法选择 不当优化后得到的控制轨线往往不光滑。只能大致反映控制随时问的变化。间接法有很多优点，从理论的 角度上讲庞特里亚金最大值原理堪称完美。通过求解由动力学方程和共轭方程组成的两点边值问题．就能 确定状态和控制随时问的变化。通过分析哈密尔顿函数还能对很多问题给出理论上的解释。但是通过最大 值原理得到的微分方程通常是一个非线性微分方程。这个微分方程的两点边值问题很难得到求解。尤其是 这组方程对共轭状态的初值选取十分敏感，一个不合理的初值就可能导致共轭状态随时间迅速增大。 混合优化算法就是利用直接法和间接法两种方法的优点，不直接求解两点边值问题，而是对这组微分方 程进行离散化，将节点上的状态和共轭状态作为优化参数，将状态和共轭状态在节点上连续作为约束条件。 利用非线性规划的方法使这个两点边值问题得到求解。从本质上讲混合法是计算变分法和参数优化方法的 结合。本文利用混合法优化了日心转移轨道。 1系统模型 对于小推力转移轨道，选择合适的动力学模型可以加快问题的收敛性，目前球坐标系动力学方程由于其 直观性和适用于任何轨道的特点，得到了广泛的应用，球坐标系下的动力学方程如下： r=口． 自：』L rcos‘0 南：坠 证华一号+丢cos例n。 (1) ‘ r 讥=一挚+坐7-tanPT+m!oos卢cos。 246 ”r 沁一学一孚唧+丢sin／／一2一—7。【一7协“中+磊 ． 丁 ”—。瓦 式中，P是中心引力体的引力常数；r为探测器与中心引力体的距离；0和≠是两个方位角；啡和口口分别 是径向、切向以及轨道面法线方向的速度分量；a和口是两个推力方向角；T是小推力发动机推力；m为探测 器的质量，g是地球加速度常数，k是发动机比冲。 对于燃料最省小推力轨道，任务的性能指标可以描述为： J=rain(‘F一‘I) (2) 其中，‘J是探测器从初始轨道变轨时间，tF是与目标轨道交会时间。 利用庞德里亚金最大值原理和式(1)和(2)，我们可以得到小推力轨道优化问题满足的两点边值问题。 共轭状态方程为： ^口=0 护_糍饥(半一