CN113944570A - Track control method based on electric pump engine and computing equipment - Google Patents

Track control method based on electric pump engine and computing equipment Download PDF

Info

Publication number
CN113944570A
CN113944570A CN202111015246.2A CN202111015246A CN113944570A CN 113944570 A CN113944570 A CN 113944570A CN 202111015246 A CN202111015246 A CN 202111015246A CN 113944570 A CN113944570 A CN 113944570A
Authority
CN
China
Prior art keywords
time length
preset
track
charging time
working time
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111015246.2A
Other languages
Chinese (zh)
Other versions
CN113944570B (en
Inventor
徐�明
左小玉
白雪
丁纪昕
陈琳
刘佳琦
和星吉
李林澄
郑亚茹
张丹阳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN202111015246.2A priority Critical patent/CN113944570B/en
Publication of CN113944570A publication Critical patent/CN113944570A/en
Application granted granted Critical
Publication of CN113944570B publication Critical patent/CN113944570B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02KJET-PROPULSION PLANTS
    • F02K9/00Rocket-engine plants, i.e. plants carrying both fuel and oxidant therefor; Control thereof
    • F02K9/42Rocket-engine plants, i.e. plants carrying both fuel and oxidant therefor; Control thereof using liquid or gaseous propellants
    • F02K9/44Feeding propellants
    • F02K9/46Feeding propellants using pumps
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02KJET-PROPULSION PLANTS
    • F02K9/00Rocket-engine plants, i.e. plants carrying both fuel and oxidant therefor; Control thereof
    • F02K9/42Rocket-engine plants, i.e. plants carrying both fuel and oxidant therefor; Control thereof using liquid or gaseous propellants
    • F02K9/44Feeding propellants
    • F02K9/56Control
    • F02K9/563Control of propellant feed pumps
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02KJET-PROPULSION PLANTS
    • F02K9/00Rocket-engine plants, i.e. plants carrying both fuel and oxidant therefor; Control thereof
    • F02K9/42Rocket-engine plants, i.e. plants carrying both fuel and oxidant therefor; Control thereof using liquid or gaseous propellants
    • F02K9/60Constructional parts; Details not otherwise provided for

Abstract

The application discloses a track control method and computer equipment based on an electric pump transmitter, wherein the method comprises the following steps: constructing a track dynamics model, obtaining a first final track state according to a preset first working time length set, a first charging time length set, a preset initial track state and the track dynamics model recursion, and judging whether a preset expected track state and the first final track state meet preset constraint conditions or not; if the working time length and the charging time length do not meet the preset requirements, the first working time length set and the first charging time length set are adjusted to obtain a second working time length set and a second charging time length set, the working time length and the charging time length corresponding to the derailment control are respectively determined from the second working time length set and the second charging time length set, and the track control is carried out according to the working time length and the charging time length. The technical problem how to carry out track control through the electric pump engine is solved in this application.

Description

Track control method based on electric pump engine and computing equipment
Technical Field
The application relates to the technical field of aerospace, in particular to an orbit control method based on an electric pump transmitter and computer equipment.
Background
In order to meet the requirements of low cost, quick launching and high reliability of the launching market of the micro spacecraft, the spacecraft design based on the electric pump engine is produced. The electric pump propulsion system drives the motor through a battery carried by the spacecraft, and then drives the fuel pump and the oxidizer pump to pump the propellant in the storage tank, so that the pressure of the propellant reaching a combustion chamber is increased, and the propellant is conveyed. Compared with the traditional turbo pump propellant supply mode, the electric pump propellant supply mode has the advantages of simple structure, convenience in thrust ratio adjustment, high energy utilization efficiency, low cost and the like, and is highly concerned by researchers at home and abroad.
For a spacecraft, after the pressure of a combustion chamber and the initial pressure of a high-pressure gas storage tank are determined, the working time of an engine needs to ensure the normal work of the spacecraft, for example, the spacecraft needs to meet the requirement of multiple orbit control operations or position maintenance operations. The working time of the transmitter is related to the battery capacity, which is related to the battery quality, which in turn affects the overall quality of the supply system; in addition, because the battery cannot be charged during the launching process of the spacecraft, the weight of the system cannot be reduced from the perspective of the battery in order to keep the spacecraft in normal operation. However, if the electric pump engine is used as the orbit control engine of the spacecraft, the solar wing of the spacecraft can be used for charging the battery for multiple times, the capacity of the battery only needs to meet the requirement of the orbit control machine for operation at every time, the battery can be charged for multiple times by using the solar wing of the spacecraft, the multiple orbit control work of the battery is guaranteed, the quality of the battery can be reduced by using the capacity of the battery, and therefore the improvement of the quality of an effective load or the mobility of the battery is reduced. Therefore, how to realize the orbit control of the spacecraft by the electric pump engine becomes a problem to be solved urgently.
Disclosure of Invention
The technical problem that this application was solved is: aiming at how to realize the orbit control of the spacecraft through the electric pump engine, the application provides an orbit control method and computer equipment based on the electric pump engine.
In a first aspect, an embodiment of the present application provides a method for controlling a track based on an electric pump transmitter, where the method includes:
constructing a track dynamics model, and recursion according to a preset first working time length set, a preset first charging time length set, a preset initial track state and the track dynamics model to obtain a first final track state, wherein the first working time length set is a set of initial working time lengths corresponding to at least one orbit control, and the first charging time length set is a set of charging working time lengths corresponding to at least one orbit control;
judging whether the expected track state and the first end track state meet preset constraint conditions or not according to the preset expected track state and the first end track state;
if not, respectively adjusting the first working time length set and the first charging time length set in a preset value range to obtain a second working time length set and a second charging time length set, and recurrently obtaining a second final orbit state again according to the second working time length set, the second charging time length set, a preset initial orbit state and an orbit dynamics model;
and determining a second working time length set and a second charging time length set corresponding to the preset expected track state and the second last track state meeting the preset constraint condition, determining the working time length from the second working time length set, determining the working time length and the charging time length corresponding to the derailment control from the second charging time length set, and performing track control according to the working time length and the charging time length.
Optionally, before constructing the orbit dynamics model, the method further includes:
the method comprises the steps of respectively determining initial working time and initial charging time corresponding to each time of rail control of the electric pump propeller according to preset rail control times, obtaining a first working time according to the initial working time corresponding to all the secondary rail controls of the electric pump propeller, and obtaining a first charging time set according to the initial charging time corresponding to all the secondary rail controls of the electric pump propeller.
Optionally, the initial working duration and the initial charging duration corresponding to each orbit control of the electric pump propeller are respectively determined according to the preset orbit control times, and the method comprises the following steps:
determining the initial working time length corresponding to each orbit control within the range not greater than the preset working time length;
and determining the initial charging time length corresponding to each orbit control within the range of not less than the preset charging time length.
Optionally, constructing an orbital dynamics model, comprising:
the orbital dynamics model is represented by the following formula:
Figure BDA0003240108180000031
Figure BDA0003240108180000032
wherein a represents a semi-major axis of a spacecraft orbit; e represents the eccentricity of the spacecraft orbit; v represents the flight velocity vector of the spacecraft relative to the atmosphere; μ represents a gravitational constant; f. oftRepresenting a perturbed acceleration tangential component in a current track velocity direction; θ represents a true proximal angle; r represents spacecraft and earthThe distance of (d); f. ofmRepresenting the normal component of the perturbed acceleration in the plane of the orbit away from the center of curvature.
Optionally, the preset constraint condition includes:
presetting that the absolute value of the difference between the long half shaft of the expected track and the long half shaft of the tail track is not greater than a preset first threshold;
the absolute value of the difference between the eccentricity of the preset expected track and the eccentricity of the last track is not greater than a preset second threshold.
Optionally, judging whether the first initial orbit state and the first final orbit state satisfy a preset constraint condition according to the first initial orbit state and the first final orbit state includes: respectively determining a first long half shaft and a first eccentricity corresponding to the last track according to the state of the first last track; calculating a first absolute value of a difference between the first major axis and a value of a major axis of the preset desired track, and a second absolute value of a difference between a value of the first eccentricity and an eccentricity of the preset desired track; and judging whether the first absolute value is not larger than a preset first threshold value or not, and whether the second absolute value is not larger than a preset second threshold value or not.
Optionally, determining the operating time length corresponding to the derailment control from the second operating time length set and determining the charging time length corresponding to the derailment control from the second charging time length set includes:
determining the working time length and the charging time length according to the following formula:
Figure BDA0003240108180000041
where x represents a second set of charging time periods, x ═ x1,x2,…,xi,…xn](ii) a q denotes a second set of operating durations, q ═ q1,q2,…,qi,…qn](ii) a f (x, q) represents a function with the second set of charging time periods and the second set of operating time periods as variables; k represents a preset weight coefficient; n represents the preset track control times; q. q.siRepresenting the working time length corresponding to the ith track control in the second working time length; x is the number ofiAnd the charging time length corresponding to the ith rail control in the second charging time length is represented.
In a second aspect, the present application provides a computer device comprising:
a memory for storing instructions for execution by at least one processor;
a processor for executing instructions stored in a memory to perform the method of the first aspect.
Compared with the prior art, the application has the following beneficial effects:
in the scheme provided by the embodiment of the application, the electric pump propeller is adopted and used as an orbit control engine, and the transfer problem based on continuous thrust is converted into the optimal design problem aiming at transfer time and fuel consumption by utilizing the charge-discharge characteristics of the battery, so that reference is provided for the actual spacecraft to apply the electric pump propeller to carry out orbital maneuver.
Drawings
Fig. 1 is a schematic flowchart of a method for controlling a track based on an electric pump transmitter according to an embodiment of the present disclosure;
FIG. 2 is a schematic diagram of a continuous thrust transfer track based on an electric pump engine provided in an embodiment of the present application;
FIG. 3 is a schematic diagram of a continuous thrust transfer trajectory based on an electric pump engine provided by an embodiment of the present application;
fig. 4 is a schematic structural diagram of a computer device according to an embodiment of the present application.
Detailed Description
In the solutions provided in the embodiments of the present application, the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
The track control method based on the electric pump transmitter provided by the embodiment of the present application is further described in detail with reference to the drawings in the specification, and a specific implementation manner of the track control method may include the following steps (a flow of the method is shown in fig. 1):
step 101, constructing a track dynamics model, and obtaining a first last track state according to a preset first working time length set, a preset first charging time length set, a preset initial track state and a track dynamics model recursion, wherein the first working time length set is a set of initial working time lengths corresponding to at least one orbit control, and the first charging time length set is a set of charging working time lengths corresponding to at least one orbit control.
For the traditional continuous thrust rail transfer problem, the working time of the engine is almost not limited, and the rail design only needs to meet a rail dynamics equation and realize the working mode of the engine with the optimal fuel. However, for an electric pump engine, the operational time after each ignition is limited by the battery capacity, and a period of unpowered coasting is required after ignition, during which the electric pump propeller battery is charged until the battery state SOC reaches 100%, ready for the next engine operation. Referring to fig. 2, a schematic diagram of a continuous thrust transfer orbit based on an electric pump engine is provided for an embodiment of the present application.
Further, the solution provided in the embodiments of the present application requires setting basic parameters before constructing the orbit dynamics model. As an example, the basic parameter is that the satellite runs on a circular orbit with an altitude of 200km, and the altitude is expected to be raised to 500km after the orbit transfer. Taking the data of the aerospace vehicle X-37B as an example for analysis, the engine thrust is 490N, the specific impulse is 315.5s, the empty weight is 3500kg, and the load is 227 kg.
In order to realize attitude and orbit control of the spacecraft based on the electric pump engine, the scheme provided by the embodiment of the application converts the problem of an engine mode into the problem of working time and charging time of the electric pump engine. The track dynamics problem is converted into a problem of optimizing the working time and the charging time of the motor of the electric pump. In order to optimize the operation time and the charging time of the electric pump engine, in a possible implementation manner, before constructing the track dynamics model, the method further comprises the following steps: and respectively determining the initial working time and the initial charging time corresponding to each orbit control of the electric pump propeller according to the preset orbit control times, and respectively obtaining a first working time set and a first charging time set according to the initial working time and the initial charging time.
In a possible implementation manner, the method includes the steps of determining the initial working time and the initial charging time corresponding to each tracking control of the electric pump propeller according to the preset tracking control times, respectively, and including: determining the initial working time length corresponding to each orbit control within the range not greater than the preset working time length; and determining the initial charging time length corresponding to each orbit control within the range of not less than the preset charging time length.
In the solution provided in the embodiments of the present application, there is a maximum thrust duration constraint t due to the electric pump propeller battery capacity limitationoptI.e. the preset operating time range is xi≤toptI is 1,2, … n. In addition, the battery needs to be charged after each ignition is finished, the spacecraft is in the unpowered sliding stage in the period, and the constraint t of the minimum charging time duration is consideredmaxChargeThen the preset charging time length range q must be satisfiedi≥tmaxCharge,i=1,2,…n。
And (4) considering perturbation influence of the earth oblateness and establishing an orbit recursion model. In practical situations, the specific impulse and the thrust of the rail-controlled engine are constant and limited, and the thrust direction of the engine is frequently adjusted to be unfavorable for improving the reliability of tasks, so that the tangential speed increment is adopted to realize the rail maneuvering. The second set of orbital perturbation equations is as follows:
Figure BDA0003240108180000061
Figure BDA0003240108180000062
Figure BDA0003240108180000063
Figure BDA0003240108180000064
Figure BDA0003240108180000065
Figure BDA0003240108180000066
wherein a, e, omega, i, omega and u respectively represent a track semi-major axis, eccentricity, ascension at a rising intersection point, a track inclination angle, a perigee argument and a latitude argument; μ represents a gravitational constant of the earth; r represents the distance of the satellite from the earth, v represents the flight velocity vector of the spacecraft relative to the atmosphere; f. oftRepresenting a perturbed acceleration tangential component in a current track velocity direction; f. ofmRepresenting the normal component of the perturbed acceleration in the plane of the orbit away from the center of curvature.
Aiming at the problem of the height lifting of the circular track, only the state changes of the semi-major axis and the eccentricity of the track need to be considered, and the track dynamics equation obtained by simplification is as follows:
Figure BDA0003240108180000071
Figure BDA0003240108180000072
wherein the power fmOnly at J2The items are related. In the continuous thrust section, shown in FIG. 2, the perturbation force fttMainly composed of the oblateness of the earth J2The perturbation and the rail control engine thrust are formed; in the unpowered gliding section, ftOnly include and J2Perturbation of the relevant terms. When the track recursion time reaches tcharge+tengineThen the iteration terminates and outputs the track state at that time.
Further, in the solution provided in this embodiment of the present application, the first working duration set is a set of initial working durations corresponding to at least one tracking control, and the first charging duration set is a set of initial charging durations corresponding to at least one tracking controlThe set is a set of charging working time lengths corresponding to at least one rail control. As an example, if the number of times of the spacecraft orbit control is set to three, the first charging time period x is ═ x1,x2,x2]Wherein x is1,x2,x3Charging time periods corresponding to the first time, the second time and the third time of rail control respectively, wherein a second working time period set q is [ q ═ q [ [ q ]1,q2,q3]Wherein q is1,q2,q3The working time lengths corresponding to the first, second and third tracking control.
Further, after a preset first working time length set and a preset first charging time length set are obtained, a first last track state is obtained according to the preset first working time length set, the preset first charging time length set, a preset initial track state and a track dynamics model, wherein the first working time length set is a set of initial working time lengths corresponding to at least one track control, and the first charging time length set is a set of charging working time lengths corresponding to at least one track control.
And 102, judging whether the preset expected track state and the first end track state meet preset constraint conditions or not according to the preset expected track state and the first end track state.
In one possible implementation, the preset constraint includes: presetting that the absolute value of the difference between the long half shaft of the expected track and the long half shaft of the tail track is not greater than a preset first threshold; the absolute value of the difference between the eccentricity of the preset expected track and the eccentricity of the last track is not greater than a preset second threshold.
Further, in a possible implementation manner, determining whether the first initial track state and the first final track state satisfy a preset constraint condition according to the first initial track state and the first final track state includes: respectively determining a first long half shaft and a first eccentricity corresponding to the last track according to the state of the first last track; calculating a first absolute value of a difference between the first major axis and a value of a major axis of the preset desired track, and a second absolute value of a difference between a value of the first eccentricity and an eccentricity of the preset desired track; and judging whether the first absolute value is not larger than a preset first threshold value or not, and whether the second absolute value is not larger than a preset second threshold value or not.
Specifically, in the solution provided in the embodiment of the present application, the preset constraint condition is represented by the following formula:
|af-a0|≤1m,|ef-e0|≤10-6
xi≤topt,qi≥tMaxCharge
wherein, a0And afSemi-major axes, e, representing respectively a preset desired trajectory and a last state trajectory0And efRespectively representing the eccentricity of a preset expected track and the eccentricity of a final state track; optimum operating time toptDetermined by the battery properties, expressed as the ratio of the battery energy density to the power density; t is tMaxchargeMay be predetermined. By way of example, the deviation of the semi-major axis and the eccentricity sets the error tolerance: tolerance of 1 meter (m) for semi-major axis and 10 for eccentricity-6. In addition, in the scheme provided in the embodiment of the present application, the number of times of tracking n is set to 3.
Step 103, if the preset expected track state and the first last track state are not satisfied, adjusting the first operating duration set and the first charging duration set within a preset value range to obtain a second operating duration set and a second charging duration set, and recurrently obtaining a second last track state according to the second operating duration set, the second charging duration set, the preset initial track state and the track dynamics model.
In one possible implementation manner, determining the operating time length corresponding to the derailment control from the second operating time length set and determining the charging time length corresponding to the derailment control from the second charging time length set includes:
determining the working time length and the charging time length according to the following formula:
Figure BDA0003240108180000081
where x represents a second set of charging time periods, x ═ x1,x2,…,xi,…xn](ii) a q denotes a second set of operating durations, q ═ q1,q2,…,qi,…qn](ii) a f (x, q) represents a function with the second set of charging time periods and the second set of operating time periods as variables; k represents a preset weight coefficient; n represents the preset track control times; q. q.siRepresenting the working time length corresponding to the ith track control in the second working time length; x is the number ofiAnd the charging time length corresponding to the ith rail control in the second charging time length is represented.
Specifically, in the solution provided in the embodiment of the present application, the circular orbit height lift problem is converted into the following optimization model:
Figure BDA0003240108180000091
and 104, determining a second working time length set and a second charging time length set corresponding to the preset expected track state and the second last track state meeting the preset constraint condition, determining the working time length from the second working time length set and determining the working time length and the charging time length corresponding to the derailment control from the second charging time length set, and performing track control according to the working time length and the charging time length.
In the solution provided in the embodiment of the present application, for solving the problem of the optimization of the orbit control, as an example, a combined optimization method of a genetic algorithm and a function set is used for solving. An optimal continuous thrust transfer trajectory considering the charging conditions is shown in fig. 3. The working time of the propeller of the two-time electric pump is x1=551.07s,x2=809.52s,x3362.81 s; the charging time of the propeller of the electric pump twice is q1=1862.58s,q2=1816.96s,q31816.41 s. The total maneuver duration is 1723s, and the fuel consumed by the optimal transfer trajectory is 169.23kg, compared to the Hulman transfer savings of 100.45 kg.
Further, in the solution provided in this embodiment of the application, if the preset expected track state and the first last track state satisfy the preset constraint condition, the operating time set corresponding to the derailment control is determined from the first operating time set, the charging time set corresponding to the derailment control is determined from the first charging time set, and the track control is performed according to the operating time set and the charging time set.
In the scheme provided by the embodiment of the application, the electric pump propeller is adopted and used as an orbit control engine, and the transfer problem based on continuous thrust is converted into the optimal design problem aiming at transfer time and fuel consumption by utilizing the charge-discharge characteristics of the battery, so that reference is provided for the actual spacecraft to apply the electric pump propeller to carry out orbital maneuver.
Referring to fig. 4, the present application provides a computer device comprising:
a memory 401 for storing instructions for execution by at least one processor;
a processor 402 for executing instructions stored in memory to perform the method described in fig. 1.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
It will be apparent to those skilled in the art that various changes and modifications may be made in the present application without departing from the spirit and scope of the application. Thus, if such modifications and variations of the present application fall within the scope of the claims of the present application and their equivalents, the present application is intended to include such modifications and variations as well.

Claims (8)

1. A method of rail control based on an electric pump dispenser, comprising:
constructing a track dynamics model, and recursion according to a preset first working time length set, a preset first charging time length set, a preset initial track state and the track dynamics model to obtain a first final track state, wherein the first working time length set is a set of initial working time lengths corresponding to at least one orbit control, and the first charging time length set is a set of charging working time lengths corresponding to at least one orbit control;
judging whether the expected track state and the first end track state meet preset constraint conditions or not according to the preset expected track state and the first end track state;
if not, respectively adjusting the first working time length set and the first charging time length set in a preset value range to obtain a second working time length set and a second charging time length set, and recurrently obtaining a second final orbit state again according to the second working time length set, the second charging time length set, a preset initial orbit state and an orbit dynamics model;
and determining a second working time length set and a second charging time length set corresponding to the preset expected track state and the second last track state meeting the preset constraint condition, determining the working time length and the charging time length corresponding to the derailment control from the second working time length set and the second charging time length set, and performing track control according to the working time length and the charging time length.
2. The method of claim 1, wherein prior to constructing the orbital dynamics model, further comprising:
the method comprises the steps of respectively determining initial working time and initial charging time corresponding to each time of rail control of the electric pump propeller according to preset rail control times, obtaining a first working time according to the initial working time corresponding to all the secondary rail controls of the electric pump propeller, and obtaining a first charging time set according to the initial charging time corresponding to all the secondary rail controls of the electric pump propeller.
3. The method of claim 2, wherein determining the initial operating time and the initial charging time for each tracking of the electric pump thruster according to the preset number of tracking respectively comprises:
determining the initial working time length corresponding to each orbit control within the range not greater than the preset working time length;
and determining the initial charging time length corresponding to each orbit control within the range of not less than the preset charging time length.
4. The method of claim 3, wherein constructing an orbital dynamics model comprises:
the orbital dynamics model is represented by the following formula:
Figure FDA0003240108170000021
Figure FDA0003240108170000022
wherein a represents a semi-major axis of a spacecraft orbit; e represents the eccentricity of the spacecraft orbit; v represents the flight velocity vector of the spacecraft relative to the atmosphere; μ represents a gravitational constant; f. oftRepresenting a perturbed acceleration tangential component in a current track velocity direction; θ represents a true proximal angle; r represents the distance of the spacecraft from the earth; f. ofmRepresenting the normal component of the perturbed acceleration in the plane of the orbit away from the center of curvature.
5. The method of claim 4, wherein the preset constraints comprise:
presetting that the absolute value of the difference between the long half shaft of the expected track and the long half shaft of the tail track is not greater than a preset first threshold;
the absolute value of the difference between the eccentricity of the preset expected track and the eccentricity of the last track is not greater than a preset second threshold.
6. The method of claim 5, wherein determining whether the first initial track state and the first end track state satisfy a predetermined constraint condition comprises:
respectively determining a first long half shaft and a first eccentricity corresponding to the last track according to the state of the first last track;
calculating a first absolute value of a difference between the first major axis and a value of a major axis of the preset desired track, and a second absolute value of a difference between a value of the first eccentricity and an eccentricity of the preset desired track;
and judging whether the first absolute value is not larger than a preset first threshold value or not, and whether the second absolute value is not larger than a preset second threshold value or not.
7. The method of claim 6, wherein determining the operating time period for the derailment control from the second set of operating time periods and determining the charging time period from the second set of charging time periods comprises:
determining the working time length and the charging time length according to the following formula:
Figure FDA0003240108170000023
where x represents a second set of charging time periods, x ═ x1,x2,…,xi,…xn](ii) a q denotes a second set of operating durations, q ═ q1,q2,…,qi,…qn](ii) a f (x, q) represents a function with the second set of charging time periods and the second set of operating time periods as variables; k represents a preset weight coefficient; n represents the preset track control times; q. q.siRepresenting the working time length corresponding to the ith track control in the second working time length; x is the number ofiAnd the charging time length corresponding to the ith rail control in the second charging time length is represented.
8. A computer device, comprising:
a memory for storing instructions for execution by at least one processor;
a processor for executing instructions stored in a memory to perform the method of any of claims 1 to 7.
CN202111015246.2A 2021-08-31 2021-08-31 Track control method based on electric pump engine and computing equipment Active CN113944570B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111015246.2A CN113944570B (en) 2021-08-31 2021-08-31 Track control method based on electric pump engine and computing equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111015246.2A CN113944570B (en) 2021-08-31 2021-08-31 Track control method based on electric pump engine and computing equipment

Publications (2)

Publication Number Publication Date
CN113944570A true CN113944570A (en) 2022-01-18
CN113944570B CN113944570B (en) 2023-02-03

Family

ID=79327670

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111015246.2A Active CN113944570B (en) 2021-08-31 2021-08-31 Track control method based on electric pump engine and computing equipment

Country Status (1)

Country Link
CN (1) CN113944570B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN203849915U (en) * 2014-05-05 2014-09-24 昆明市官渡区第六中学 Double-layer-track electric pump water rocket launcher
US20160208670A1 (en) * 2015-01-21 2016-07-21 Deere & Company Def injector cooling system and method
CN106697333A (en) * 2017-01-12 2017-05-24 北京理工大学 Robustness analysis method for spacecraft orbit control strategy
CN107031868A (en) * 2017-03-23 2017-08-11 北京空间飞行器总体设计部 A kind of low rail remote sensing satellite autonomous Orbit control method
CN109630317A (en) * 2018-12-06 2019-04-16 上海空间推进研究所 Rail attitude control integration space propulsion system based on electrodynamic pump
CN112412660A (en) * 2020-12-03 2021-02-26 西安航天动力研究所 Space power system combining extrusion and electric pump auxiliary pressurization

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN203849915U (en) * 2014-05-05 2014-09-24 昆明市官渡区第六中学 Double-layer-track electric pump water rocket launcher
US20160208670A1 (en) * 2015-01-21 2016-07-21 Deere & Company Def injector cooling system and method
CN106697333A (en) * 2017-01-12 2017-05-24 北京理工大学 Robustness analysis method for spacecraft orbit control strategy
CN107031868A (en) * 2017-03-23 2017-08-11 北京空间飞行器总体设计部 A kind of low rail remote sensing satellite autonomous Orbit control method
CN109630317A (en) * 2018-12-06 2019-04-16 上海空间推进研究所 Rail attitude control integration space propulsion system based on electrodynamic pump
CN112412660A (en) * 2020-12-03 2021-02-26 西安航天动力研究所 Space power system combining extrusion and electric pump auxiliary pressurization

Also Published As

Publication number Publication date
CN113944570B (en) 2023-02-03

Similar Documents

Publication Publication Date Title
CN110794863B (en) Heavy carrier rocket attitude control method capable of customizing control performance indexes
Lee et al. Optimal path planning of solar-powered UAV using gravitational potential energy
CN109911249B (en) Interstellar transfer limited thrust orbit-entering iterative guidance method for low thrust-weight ratio aircraft
CN105353621A (en) Fault mode thrust allocation method for geostationary orbit satellite electric thruster
CN112550770B (en) Rocket soft landing trajectory planning method based on convex optimization
CN101393458A (en) Longitudinal control method for high altitude lifting of aerospaceplane
CN112486196B (en) Aircraft rapid trajectory optimization method meeting strict time and position constraints
CN113525721B (en) Satellite orbit transformation method, device, electronic equipment and storage medium
CN112329136A (en) Carrier rocket online flight program reconstruction method based on balanced flight theory
Peng et al. Free return orbit design and characteristics analysis for manned lunar mission
CN112182772A (en) Rocket propulsion control method, device and storage medium
CN113581494B (en) GEO satellite electric hybrid propulsion orbital transfer method
CN109190155B (en) Hybrid continuous low-thrust track design method adopting electric propulsion/solar sail propulsion
CN112660426B (en) Rocket soft landing guidance method
CN113944570B (en) Track control method based on electric pump engine and computing equipment
CN112507461B (en) Starting method of engine at power soft landing stage of carrier rocket
CN106227968A (en) A kind of spacecraft sustainer Burnout optimization method
Wang et al. Six-DOF trajectory optimization for reusable launch vehicles via Gauss pseudospectral method
CN116301028B (en) Multi-constraint on-line flight trajectory planning middle section guiding method based on air suction hypersonic speed platform
CN104914873A (en) Coupling method for attitude and orbit control engine
CN114690793B (en) Sliding mode control-based reusable carrier rocket vertical soft landing guidance method
CN114781067A (en) ASV separation trajectory design method based on course differential motion model
CN115494727A (en) Carrier rocket orbit-entering trajectory planning method based on orbit prediction
CN112325711A (en) Carrier rocket orbit height keeping control method based on balanced flight theory
Wang et al. Trajectory optimization for reusable rocket landing

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant