[发明专利]一种火箭推力下降故障下椭圆救援轨道的轨迹重规划方法

权利要求书

1.一种火箭推力下降故障下椭圆救援轨道的轨迹重规划方法，其特征在于，包括：

在地心惯性坐标系中建立火箭的上升段二级飞行动力学方程，以不同的故障时刻、推力下降百分比设置边界条件和约束条件，构建各种推力下降故障下椭圆救援轨道的轨迹优化问题；

采用自适应伪谱法离线求解椭圆救援轨道的轨迹优化问题，得到故障状态-入轨参数的样本集，其中样本集的输入特征为故障状态，所述故障状态包括推力故障的时刻、推力下降的大小、位置、速度和质量，其中样本集的输出特征为入轨参数，所述入轨参数包括轨道半长轴、偏心率、轨道倾角、升交点赤经、终端的控制变量；

采用最大最小法对样本数据进行归一化处理，将所有数据规范化到[-1,1]之间，通过正交最小二乘法选择径向基神经网络数据中心，其中径向基函数为高斯基函数，离线训练径向基神经网络，从而建立故障状态到入轨参数非线性映射关系；

将径向基神经网络迁移到火箭实际飞行中，以飞行的故障状态作为输入，该径向基神经网络在线决策出入轨参数；利用决策出的救援轨道根数和终端控制变量，为轨迹规划设置合理的初始猜测值；利用决策出的救援轨道根数将椭圆救援轨道的轨迹优化问题转化为最优推进剂优化问题，采用自适应伪谱法在线求解最优推进剂优化问题，即可得到飞行轨迹。

2.根据权利要求1所述一种火箭推力下降故障下椭圆救援轨道的轨迹重规划方法，其特征在于，构建各种推力下降故障下椭圆救援轨道的轨迹优化问题时：

设X₁轴在赤道平面内指向发射时刻本初子午线方向，Z₁轴垂直赤道平面指向北极，Y₁轴满足右手定则，在地心惯性坐标系中建立火箭的上升段二级飞行动力学方程如下：

式中r，v分别为运载火箭的位置、速度向量，为r的一阶导数，为v的一阶导数；μ＝GM为地球引力常数，G为引力常量，M为地球质量；m为火箭的总质量，为m的一阶导数，I_sp为火箭的发动机比冲；u＝[u_x,u_y,u_z]^T为发动机的推力单位矢量分量；当发动机发生故障，推力下降的比例是η，推力大小是(1-η)T_nom，T_nom是发动机标称推力；g₀为重力加速度；

在推力下降故障情况下，发动机比冲不变，推进剂的秒耗量下降η，总飞行时间超过标称的飞行时间，假定发动机推力下降故障出现时刻为t₀，运载火箭需要以t₀时刻的状态为起始点优化救援轨迹，因此起始点等式约束条件表示为：

x(t₀)＝x₀ (4)

式中，x₀是起始点的状态。

用m_f表示剩余燃料耗尽后运载火箭与有效载荷的总质量，地球半径是R₀；定义最低安全轨道高度为h_safe，终端质量和高度应当满足：

m(t_f)≥m_f,h_safe≤r(t_f)-R₀ (5)

发生推力下降故障，为椭圆救援轨道类型，在调整轨道平面的同时，尽量提高轨道高度；椭圆救援轨道的轨迹优化问题描述为：

其中J是优化目标，λ_a,λ_e,λ_i,λ_Ω为权重因子，hp_f对应入轨轨道的近地点高度，是目标轨道的近地点高度,t_f为终端时刻，a_f,i_f,Ω_f分别为入轨轨道的半长轴、倾角、升交点赤经，分别为目标轨道的半长轴、倾角、升交点赤经。