一种小推力转移轨道快速设计与优化方法技术

技术编号:8980804 阅读:157 留言:0更新日期:2013-07-31 22:52
本发明专利技术提出了一种小推力转移轨道快速设计与优化方法。所述方法引进逆6次多项式法作为初始设计,为后续的精确优化提供高性能的优化参数初值,减少了整个优化过程的计算时间;采用Radau伪光谱法将小推力轨道优化问题参数化为一个非线性规划问题,不仅能减少优化参数的数量还提高了优化结果的精度;无需推导繁琐的一阶最优必要条件,即使是处理含有各种复杂约束的小推力转移轨道,操作起来也方便简单。该方法在太阳帆推进、太阳能电推进等小推力推进航天器的转移轨道设计与优化领域有着极其重要的应用价值。

【技术实现步骤摘要】

本专利技术属于深空探测
,尤其是涉及。
技术介绍
轨道设计与优化是人类开展深空探测活动必须要解决的关键问题之一。小推力轨道设计与优化方法主要分为间接法、直接法和混合法。传统的间接法需要推导繁琐的一阶最优必要条件和猜测没有实际物理意义的协态变量初值,且这种方法的收敛域很窄,对初值的敏感度高。因此,间接法的操作比较困难。传统的混合法同样需要推导繁琐的一阶最优必要条件。尤其,当需要考虑星历、最大推力加速度等约束时,推导过程异常复杂,操作起来同样比较困难。传统的直接法虽然省去了一阶最优必要条件的推导过程,容易操作。但是,由于需要求解一个变量很多的非线性规划问题,计算耗时且计算精度得不到保证。因此,十分有必要发展一种容易操作、计算速度快且计算精度高的小推力轨道设计与优化方法。
技术实现思路
本专利技术所要解决的技术问题在于克服现有技术的不足,针对小推力轨道设计与优化技术求解繁琐、计算量大且精度不高的问题,提出。所述方法引入的改进包括:将目前被认为用于求解最优控制问题最具有潜力的Radau伪光谱法(简称RPM)引入到小推力轨道的优化设计中;增加了存在解析解的小推力初始轨道设计,为后面的精确优化提供高性能的优化参数初值,不仅保留了直接法的易操作性,大幅度地提高了计算效率。为解决上述技术问题,本专利技术所采用的技术方案是:,包括如下步骤:步骤A,小推力转移轨道设计与优化问题描述成如下形式的最优控制问题:目标函数为:J= Φ (X (t0),t0, X (tf),tf)动力学约束为:

【技术保护点】
一种小推力转移轨道快速设计与优化方法,其特征在于,包括如下步骤:步骤A,小推力转移轨道设计与优化问题描述成如下形式的最优控制问题:目标函数为:J=Φ(x(t0),t0,x(tf),tf)动力学约束为:dxdt=f(x(t),u(t),t)边界约束为:φ(x(t0),t0,x(tf),tf)=0路径约束为:C(x(t),u(t),t)≤0其中x(t)表示包含小推力优化问题中的位置矢量、速度矢量、航天器剩余质量的状态变量,t表示时间变量,t0表示发射时刻,tf表示到达时刻,Φ表示目标函数,J表示目标函数值,f表示动力学约束函数,u表示控制变量,φ表示边界约束函数,C表示路径约束函数;步骤B,小推力转移轨道初始设计步骤B?1,采用逆6次多项式法对小推力转移轨道进行初始设计,即利用逆6次多项式来逼近航天器的真实飞行轨迹,从而求得航天器的速度、位移、加速度、剩余质量的近似值;步骤B?2,结合逆6次多项式法和遗传算法完成小推力转移轨道初始设计;步骤C,小推力转移轨道精确优化设计步骤C?1,采用Radau伪光谱法,将连续系统的小推力转移轨道优化问题转换为非线性规划问题,所述非线性规划问题表述如下:minJ=Φ(X(τ0),τ0,X(τN+1),τN+1)满足:DXLGR-tf-t02f(Xk,Uk,τk;t0,tf)=0,(k=1,...,N)φ(X(τ1),τ1,X(τN+1)τN+1)=0tf-t02C(Xk,Uk,τk;t0,tf)≤0,(k=1,...,N)其中,D为微分矩阵;τk指的是第k个离散点值;N表示配点的数量;Xk表示第k个离散点的状态变量值;Uk表示第k个离散点的控制变量值;寻优变量有:状态变量XLGR=(X1,...,XN+1)T、控制变量(U1,...,UN)T以及始末端的时刻分别为t0和tf;(?)T表示矩阵转置;步骤C?2,利用序列二次规划算法SQP求解上述非线性规划问题,得到小推力转移轨道的精确值。...

【技术特征摘要】
1.一种小推力转移轨道快速设计与优化方法,其特征在于,包括如下步骤: 步骤...

【专利技术属性】
技术研发人员:李爽朱永生王宇凯江秀强
申请(专利权)人:南京航空航天大学
类型:发明
国别省市:

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1