LU500286B1

LU500286B1 - An optimal rescue orbital elements online decision-making method based on RBFNN for launch vehicles under thrust drop fault

Info

Publication number: LU500286B1
Application number: LU500286A
Authority: LU
Inventors: Zhigang Wu; Shujun Tan; Xiao He; Liyong Zhang
Original assignee: Univ Dalian Tech
Priority date: 2020-11-12
Filing date: 2021-06-16
Publication date: 2021-12-16
Also published as: CN112455723B; US20220147820A1; US12242962B2; CN112455723A

Abstract

The present invention discloses an optimal rescue orbital elements online decision- making method based on RBFNN for launch vehicles under thrust drop fault? including: establishing the flight dynamic equations of launch vehicles in the second-stage ascending phase in the geocentric inertial coordinate system, to construct a series of optimization problems of maximum semi-major axis of circular orbit under the thrust drop fault; using the adaptive pseudo-spectrum method to solve the optimization problems of maximum semi-major, and using the maximum and minimum method to normalize the sample data to [-1, 1], using the orthogonal least square method to select the data center of the radial basis function neural network (RBFNN), where the Gaussian function is selected as the radial basis function, and the RBFNN is trained offline to establish a nonlinear mapping relationship from the fault states to the optimal rescue orbital elements. The RBFNN is transferred to the actual flight of launch vehicles, and taking the fault states of the flight as the input, the RBFNN determines the optimal rescue orbital elements online.

Description

BL-5254 AN OPTIMAL RESCUE ORBITAL ELEMENTS ONLINE DECISION-MAKING METH&#p00286

BASED ON RBFNN FOR LAUNCH VEHICLES UNDER THRUST DROP FAULT Technical Field The present invention relates to the technical field of launch vehicles under thrust drop fault, in particular to an optimal rescue orbital elements online decision-making method based on RBFNN for launch vehicles under thrust drop fault.

Background As a rocket power plant, the engine is a decisive factor for the reliability and safety of a rocket. Its reliability is related to the success or failure of the entire flight mission. In the actual flight mission, when the thrust drop fault of the launch vehicle engine occurs, it is difficult to complete the mission if the guidance and control scheme under nominal ballistic conditions is continued. In order to avoid the fall of the payload, it is necessary to decide the optimal rescue orbital elements online according to the fault state so that the payload enters into the rescue orbit.

At present, for the problem of online rescue under the thrust drop failure, the main idea is to couple the two problems of rescue trajectory decision and trajectory optimization based on the dynamic model for online optimization [1,2]. SONG [1] proposed an autonomous rescue strategy and algorithm for launch vehicle thrust descent failure, combined with the injection point geocentric angle estimation, convex optimization, adaptive collocation method, etc., to provide a good online solution for the original coupling optimization problem, online providing a good initial value for the original coupling optimization problem to improve the computational efficiency of online trajectory planning. LI [2] proposed an online trajectory optimization algorithm for launch vehicles based on convex optimization to obtain a high-precision and high-efficiency online trajectory optimization method. The above method couples the two problems of rescue trajectory decision and flight trajectory together and optimizes them. Due to the unknown rescue orbit, the search space of the optimal solution is very large, which affects the calculation efficiency of the online part.

Summary

According to the problems exist in the prior art, the present invention provides an optimal rescue orbital elements online decision-making method based on RBFNN for launch vehicles under thrust drop fault, specifically comprising: establishing the dynamic equations of launch vehicles in the second-stage ascending phase in the geocentric inertial coordinate system; setting the boundary condition and the constraint condition with different fault times and thrust drop percentages, to construct a series of optimization problems of maximum semi-major axis of circular orbit under the thrust drop fault; using the adaptive pseudo-spectral method (APM) to solve the optimization problems of maximum semi-major axis offline , to obtain a sample set of the optimal rescue orbital elements under different fault states, wherein the input feature of the sample set is the fault state, and the fault state includes the time of the thrust fault, the percentage of thrust drop, the position, the speed, the mass and the output features of the sample set are the orbital elements of the rescue orbit, the orbital elements of the rescue orbit includes the semi- major axis, the inclination, and the longitude of the ascending node of the rescue orbit; using the maximum-minimum method to normalize the sample data, and all data are normalized to [-1,1]; using the orthogonal least square method to select a data center of a radial basis function neural network (RBFNN), wherein a Gaussian basis function is selected as the radial basis function; the radial basis function neural network is trained offline to establish a nonlinear mapping relationship from the fault state to the optimal rescue orbital elements; transferring the well-trained radial basis function neural network to the actual flight; taking the fault state of as the input of the radial basis function neural network, determining the optimal rescue orbital elements online.

Furthermore, when constructing the optimization problems of maximum semi-major axis: setting the Xi axis to point the primary meridian at the time of launch in an equatorial plane, the Z1 axis vertical to the equatorial plane and to point the North Pole, and the Y1 axis to meet the right-hand rule; establishing the dynamic equations of the launch vehicles in the second-stage ascending phase in the geocentric inertial coordinate system as follows: 17000888 F=V (1) v= A r+ EM om u I 7 (2) i = — Sohn (3) Golsp wherein r and v represent the position and the velocity vector of launch vehicles; u=GM is an earth's gravitational constant; represents a total mass of the launch vehicles, and i; represents a specific impulse of the engine of the launch vehicles; 4 = {it y. ul" is a component of a thrust unit vector of the engine; when the engine fault occurs, a percentage of thrust drop is = ; a thrust magnitude is {1 — #}T um, Wherein Ta is a nominal thrust of the engine; in the case of a thrust drop fault, the specific impulse of the engine remains unchanged, a propellant consumption per second decreased;, and a total flight time exceeds a nominal flight time; assuming the engine thrust drop failure occurs at tq, SO the constraint condition of the starting point is expressed as follow: x{ta) = x, 4) wherein x, is the state of the starting point, a nonlinear relationship from the number of orbital elements to a terminal state is expressed as [aps ep ip 2p, 0017 = Prt) v(t) (5) wherein £; is a terminal moment; €, ef, Er, Ar, wy is the orbital elements of the target orbit including the semi-major axis, eccentricity, inclination, longitude of the ascending node, argument perigee of the terminal point.

a total mass of the launch vehicle and payload after the fuels exhaustion is expressed as my, the radius of the earth is R,, and the minimum safe orbit height is defined as Ê re the terminal mass and the height meets: mt) = mp, À, É v(t) — Ro (6) when the thrust drop fault occurs, because an energy required to send a load into a circular orbit is less than an energy of an elliptical orbit under a same perigee height, so searching for a highest circular orbit within the current orbital plane as the optimal rescue orbital elements; describing, by the solution to the highest circular orbit under the thrust 0200 drop failure, the maximum optimization problem of the semi-major axis: Objective Function: min] = —a(t,) Boundary Conditions: [ey ir, 2, 17 = [0,i1%%, 27°17, Eq.(4), Eq.(6) Dynamics Constraints: Eq. (1)-(3) Control Constraints: ~~ {ul} = 1 The orbital inclination {= and the longitude of the ascending node (2; is in accordance with the fault states of launch vehicles. Due to the above technical scheme, the present invention provides an optimal rescue orbital elements online decision-making method based on RBFNN for launch vehicles under thrust drop fault, which has the following beneficial effects: (1) Benefiting from the invention, the complicated trajectory optimization problem of launch vehicles under the thrust drop fault is decoupled into an optimal rescue orbital elements decision-making and a simplified trajectory optimization problem.

(2) RBFNN can fit the nonlinear relationship from the optimal rescue orbital elements to the fault states accurately. The online decision-making time is almost negligible compared with trajectory optimization time.

(3) The invention can quickly determine the optimal rescue orbital elements online and provide reasonable terminal constraints for subsequent online trajectory optimization. Then the complexity of the online trajectory optimization problem will also be greatly reduced, thereby improving the overall computational efficiency of online rescue.

Brief Description of the Drawings In order to explain more clearly the embodiment of this application or the technical scheme in the prior art, a brief description of the attached drawings required in the specific embodiment or the existing technical description is given below. Obviously, the attached drawings in the following description are only some of the embodiments recorded in this application. For general technical personnel in this field, other drawings can be obtained without creative work.

Fig. 1 shows the scheme of the optimal rescue orbital elements decision-making strategy.

Fig. 2 shows the flow chart of sample set generation. 7500286 Fig. 3 shows the results of the optimal rescue orbital elements decision-making. Detailed Description In order to make the technical scheme and advantages of the invention clearer, a clear and complete description of the technical scheme in the embodiments of the invention will be described in conjunction with drawings attached in the embodiments of the invention hereinafter.

Fig. 1 shows an optimal rescue orbital elements online decision-making method based on RBFNN for launch vehicles under a thrust drop fault, specifically steps are as follows: S1: Taking the two-stage launch vehicle as the research object, it is assumed that the thrust drop failure occurs in the second stage of flight. At this time, the atmosphere is thin and the influence of aerodynamic forces can be ignored. Define the geocentric inertial coordinate system: the origin is at the center of the earth, the X1 axis points to the primary meridian at the time of launch in the equatorial plane, the Z1 axis is vertical to the equatorial plane and points to the North Pole, and the Y1 axis satisfies the right-hand rule, respectively. Establishing the dynamic equations of the launch vehicles in the second-stage ascending phase in the geocentric inertial coordinate system as follows: r=v (1) Ya À. p+ 0 on i bo (2) . T= Tham mp = — Fam Tob = (3) wherein ¥ and” represent the position and the velocity vector of launch vehicles; u=GM is an earth’s gravitational constant; ** represents a total mass of the launch vehicles, and I, represents a specific impulse of the engine of the launch vehicles; # = fre, 1, al‘ isa component of a thrust unit vector of the engine; when the engine fault occurs, a percentage of thrust drop is 1; a thrust magnitude is {1 — 737,,m, Wherein T.……m is à nominal thrust of the engine; in the case of a thrust drop fault, the specific impulse of the engine remains unchanged, a propellant consumption per second decreased n, and a total flight time exceeds a nominal flight time; assuming the engine thrust drop failure occurs at t,, so the constraint condition of the starting point is expressed as follow: x(ty) = x, 4) wherein x, is the state of the starting point, a nonlinear relationship from the number of orbital elements to a terminal state is expressed as fap, ef. ip, 2p, we] = (rte) v(t) (5) wherein tr is a terminal moment; Gr, pip, lp, wp is the orbital elements of the target orbit including the semi-major axis, eccentricity, inclination, longitude of the ascending node, argument perigee of the terminal point.

a total mass of the launch vehicle and payload after the fuels exhaustion is expressed as my, the radius of the earth is R,, and the minimum safe orbit height is defined as À re the terminal mass and the height meets: mt) = mp, À, € TU Ro (6) when the thrust drop fault occurs, because an energy required to send a load into a circular orbit is less than an energy of an elliptical orbit under a same perigee height, so searching for a highest circular orbit within the current orbital plane as the optimal rescue orbital elements; describing, by the solution to the highest circular orbit under the thrust drop failure, the maximum optimization problem of the semi-major axis: Objective Function: min} = —a(t;) Boundary Conditions: fey, ir. {2 1" = [0, ie, een Eq.(4), Eq.(6) Dynamics Constraints: Eq. (1) - (3) Control Constraints: {full = 1 The orbital inclination .** and the longitude of the ascending node £2. is in accordance with the fault states of launch vehicles.

S2: In the actual flight, the position, speed and mass of the rocket are merely near the nominal trajectory due to the noise disturbance. Therefore, the state of the fault is the time when the fault occurs, the size, location, speed, and mass of the thrust drop. Traverse different fault states, construct the optimization problem, and use the adaptive pseudo-spectrum method to solve the optimization problem, as shown in Figure 2. The 786 sample set of optimal rescue orbital elements versus fault states is the training date for RBFNN.

S3: The data of sample set are normalized into the range [-1...1] in order to eliminate the order of the magnitude difference between the data of each dimension and avoid the large prediction error caused by the large magnitude difference between the input and output data. The orthogonal least square method is used to select the RBF network data center, and the nonlinear mapping the relationship from the fault state to the optimal rescue orbital elements is established.

The hidden unit is activated by the basis function, the Gaussian basis function is used in present application. The output of j-th hidden layer is: [ities | Oj = exp he AN | (7) | 20, J wherein XS is the { -th input data. g; is the center of the basis function of the hidden node, and 7; is the expansion speed of the radial basis function. The center of the basis function and the expansion speed are determined by the training of the network. The output layer linearly combines the output of the hidden layer of the radial basis function to generate the expected output. The output of the kth node in the output layer is

A wed wep, (E=L-5) = (8) Wherein w,, is the weight from the jth hidden layer neuron to the à output neuron. In order to achieve appropriate approximation accuracy, the following parameters are determined through training: the number of hidden layer neurons, the center of the basis function of each hidden layer neuron, and the weight of the radial basis function output passed to the summation layer.

S4: The well-trained RBFNN is transferred to online application for decision-making of the optimal rescue orbital elements, taking the fault states of the actual flight as the input,

can quickly determine the optimal rescue orbital elements online. 7500286 Example: In this section, simulation verification is conducted in the entire second stage flight phase of a launch vehicle, and the parameters are from reference [1]. The state distribution of the faults established by the sample set is based on the fault time is incremented by 1s from Os to 375s, and the proportion of thrust-drop is incremented by 1% from 13% to 40%. Excluding the fault conditions that can enter the target orbit and the rescue orbit height is less than 160km, the fault conditions that require rescue are shown in Figure 3. In the sample set, 90% of the data is randomly selected as the training set, and the remaining 10% of the data is used as the test set. A radial basis function neural network is used to establish a nonlinear mapping from the fault state to the optimal rescue orbital elements. The spread factor of RBF neural network training is 1, and the final number of trained hidden layer neurons is 50. In the test set, the optimal rescue orbital elements are determined by the RBF neural network are shown in Figure 3. The error of the semi-major axis decision is [-0.148, 0.774km], and the relative error is within 0.015%. The error of orbital inclination and the longitude of the ascending node is within+10-*deg, and the maximum relative error is within 10.

Table 1 RMSE of different machine learning models No. Mechanism learning model RMSE (km) 1 Linear regression 3.854 2 Decision tree 2.225 3 SVM (Gauss) 1.214 4 SVM (Cubic) 1.111 Gaussian process regression 0.129 6 Radial basis function neural network 0.105 Table 1 is a comparison of the root mean square error (RMSE) of the results of different machine learning methods. The RMSE obtained by linear regression is larger than that of the nonlinear regression method, because the relationship from the fault state to the optimal rescue orbital elements are nonlinear, and it is suitable to use the nonlinear function approximation method. Compared with several other machine models, the RBF neural network model has a smaller RMSE, which means that it has a better function 786 approximation effect in mapping the relationship from the fault state to the optimal rescue orbital elements.

The above is the preferred implement approach of the invention, the protection range is not limited to what mentioned above. Any technician who is familiar with this technical field replaces or revises technical schemes and inventive created on the basis of the present invention within the technical scope disclosed in this invention, schemes and inventive should be involved in the protection range of this invention.

[1] Z. Song, C. Wang, Q. Gong, Joint dynamic optimization of the target orbit and flight trajectory of a launch vehicle based on state-triggered indices, Acta Astronaut., 174 (2020) 82-93.

[2] Y.Li, B.J.Pang, C.Z Wei, N.G.Cui, Y.B.Liu, Online trajectory optimization for power system fault of launch vehicles via convex programming, Aerosp. Sci. Technol., 98(2020) P.1 (1-10).

Claims

Claims LU500286

1. An optimal rescue orbital elements online decision-making method based on RBFNN for launch vehicles under thrust drop fault, comprising: establishing the dynamic equations of launch vehicles in the second-stage ascending phase in the geocentric inertial coordinate system; setting the boundary condition and the constraint condition with different fault times and thrust drop percentages, to construct a series of optimization problems of maximum semi-major axis of circular orbit under the thrust drop fault using the adaptive pseudo-spectral method (APM) to solve the optimization problems of maximum semi-major axis offline , to obtain a sample set of the optimal rescue orbital elements under different fault states, wherein the input features of the sample set are the fault state, and the fault state includes the time of the thrust fault, the percentage of thrust drop, the position, the speed, the mass and the output features of the sample set are the rescue orbital elements, the rescue orbital elements include the semi-major axis, the inclination, and the longitude of the ascending node of the rescue orbit; using a maximum-minimum method to normalize a sample data, and all data are normalized to [-1,1]; using the orthogonal least square method to select a data center of a radial basis function neural network (RBFNN), wherein a Gaussian basis function is selected as the radial basis function; the radial basis function neural network is trained offline to establish a nonlinear mapping relationship from the fault state to the optimal rescue orbital elements; transferring the well-trained radial basis function neural network to the actual flight; taking the fault state of as the input of the radial basis function neural network, to determine the optimal rescue orbital elements online.

2. The method according to claim 1, wherein, when constructing the optimization problems of maximum semi-major axis: setting the Xi axis to point a primary meridian at the time of launch in an equatorial plane, the Z1 axis vertical to the equatorial plane and to point the North Pole, and the Y1 axis to meet the right-hand rule; establishing the dynamic equations of the launch vehicles in the second-stage ascending phase in the geocentric inertial coordinate system as 7500286 follows: PS (1) P= te a ee H I " (2)

i = — Oh 3) wherein 7 and © represent the position and the velocity vector of launch vehicles; u=GM is an earth’s gravitational constant; =: represents a total mass of the launch vehicles, and I, represents a specific impulse of the engine of the launch vehicles; te = fax, [TR is a component of a thrust unit vector of the engine; when the engine fault occurs, a percentage of thrust drop is 1; a thrust magnitude is {1 — #)T.om wherein T,m is à nominal thrust of the engine; in the case of a thrust drop fault, the specific impulse of the engine remains unchanged, a propellant consumption per second decreased and a total flight time exceeds a nominal flight time; assuming the engine thrust drop failure occurs at #; so the constraint condition of the starting point is expressed as follow:

x{te) = Xp (4) wherein Xx; is the state of the starting point, a nonlinear relationship from the number of orbital elements to a terminal state is expressed as lap, ep tp lp wp] = Bre), vit) (5) wherein ¢, is a terminal moment; a;,e., ir, 2; wy is the orbital elements of the target orbit including the semi-major axis, the eccentricity, the inclination, the longitude of the ascending node, the argument perigee of the terminal point. a total mass of the launch vehicle and payload after the fuels exhaustion is expressed as my, the radius of the earth is Æ,, and the minimum safe orbit height is defined as Ë the terminal mass and the height meets: mit) > my, A = Ft) — Ro (6) when the thrust drop fault occurs, because an energy required to send a load into a circular orbit is less than an energy of an elliptical orbit under a same perigee height, so searching for a highest circular orbit within the current orbital plane as the optimal pese 799 orbital elements; describing, by the solution to the highest circular orbit under the thrust drop failure, the maximum optimization problem of the semi-major axis:

Objective Function: min J = —aG,)

Boundary Conditions: Îe,, i, 0:7 = [0 EEE, Eq.(4), Eq.(6)

Dynamics Constraints: Eq. (1) - (3)

Control Constraints: jiujl = 1 determining the orbital inclination à"? and the longitude of the ascending node 12” according to the fault states of launch vehicles.