CN109649689B - Limited thrust orbital transfer gravity loss calculation method and thrust calculation device - Google Patents
Limited thrust orbital transfer gravity loss calculation method and thrust calculation device Download PDFInfo
- Publication number
- CN109649689B CN109649689B CN201811495695.XA CN201811495695A CN109649689B CN 109649689 B CN109649689 B CN 109649689B CN 201811495695 A CN201811495695 A CN 201811495695A CN 109649689 B CN109649689 B CN 109649689B
- Authority
- CN
- China
- Prior art keywords
- thrust
- orbital transfer
- loss
- satellite
- general assembly
- 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.)
- Active
Links
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/242—Orbits and trajectories
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Toys (AREA)
- Testing Of Engines (AREA)
Abstract
A gravity loss result can be obtained only according to the length of an ignition arc segment, the thrust magnitude of an engine and initial and tail end orbit information do not need to be obtained, and complicated orbit control model integral operation is not needed. The method is oriented to the field of relative navigation tests, can quickly calculate the propellant added when the satellite adopting the near-circular orbit changes orbit, and has strong reference value for the relative navigation test remote guidance section orbit changing task.
Description
Technical Field
The invention relates to a method and a device for calculating the limited thrust orbital transfer gravity loss, and belongs to the technical field of satellite orbital transfer.
Background
The relative navigation test generally involves two or more satellites, and due to the limitation of relative navigation instruments, the distance between the satellites used for the relative navigation test needs to satisfy certain constraints. To achieve this constraint, the relative navigation masters cannot avoid performing a orbital transfer maneuver to arrive at a configuration that can be used for relative navigation trials.
At present, orbital transfer engines adopted by satellites are limited in thrust, and ignition orbital transfer is mostly carried out by adopting an inertial space orientation method. In the traditional engine gravity loss calculation method, integration operation is mostly carried out on an orbit control model to obtain the speed increment of orbital transfer, and in addition, the speed increment of pulse-form orbital transfer needs to be obtained, so that the gravity loss of satellite orbital transfer is obtained. The calculation method needs to know information such as the thrust magnitude and the initial and final end orbit, and cannot calculate on the premise of no information.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: the method and the device for calculating the orbital transfer gravity loss of the limited thrust overcome the defects of the prior art, obtain a simple relation between the gravity loss and the ignition arc length, and can be applied to quick calculation of the orbital transfer gravity loss of the near-circular orbit.
The purpose of the invention is realized by the following technical scheme:
a method for calculating the gravity loss of limited thrust orbital transfer is suitable for a near-circular orbit satellite and adopts inertial directional orbital transfer, and the gravity loss eta of the limited thrust orbital transfer is as follows:
In the method for calculating the limited thrust orbital transfer gravity loss, the orbital transfer engine thrust F of the near-circular orbit satelliteGeneral assemblyComprises the following steps:
in the formula (I), the compound is shown in the specification,the thrust for effective orbital transfer is generated in the thrust of the engine.
The above-mentioned limited thrust becomesMethod for calculating orbital gravity loss, propellant loss M of near-circular orbit satelliteLoss of powerComprises the following steps:
Mloss of power=η×MGeneral assembly
In the formula, MGeneral assemblyThe total propellant demand for the satellite orbital transfer.
According to the method for calculating the limited thrust orbital transfer gravity loss, the orbital eccentricity of the near-circular orbit satellite is not more than 0.0025, and the orbital height is not more than 2000 km.
A thrust computing device can compute the thrust F of a orbital transfer engine of a near-circular orbit satelliteGeneral assemblyThrust F of orbital transfer engine of said satelliteGeneral assemblyComprises the following steps:
wherein
Wherein eta is the limited thrust orbital transfer gravity loss,generating effective orbital transfer thrust for engine thrust;is the ignition arc segment length.
The thrust computing device, the propellant loss M of the near-circular orbit satelliteLoss of powerComprises the following steps:
Mloss of power=η×MGeneral assembly
In the formula, MGeneral assemblyThe total propellant demand for the satellite orbital transfer.
According to the thrust calculation device, the orbital eccentricity of the near-circular orbit satellite is not more than 0.0025, and the orbital height is not more than 2000 km.
Compared with the prior art, the invention has the following beneficial effects:
(1) compared with the prior art, the method can obtain the gravity loss result only according to the length of the ignition arc segment, does not need to know the thrust magnitude of the engine and the initial and tail end orbit information, and does not need complicated orbit control model integral operation;
(2) the method is oriented to the field of relative navigation tests, can quickly calculate the propellant added when the satellite adopting the near-circular orbit changes orbit, and has extremely high reference value for the relative navigation test remote guidance section orbit changing task;
(3) the method has good calculation real-time performance and is beneficial to realizing the rapid calculation of the limited thrust of the orbital transfer engine.
Drawings
Fig. 1 is a position relationship diagram of the ignition process of the satellite according to the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, embodiments of the present invention will be described in detail with reference to the accompanying drawings.
A method for calculating the loss of gravity by changing orbit with limited thrust is suitable for changing orbit of a near-circular orbit satellite by inertial orientation, the orbit eccentricity of the near-circular orbit satellite is not more than 0.0025, and the orbit height is not more than 2000 km. The limited thrust orbital transfer gravity loss eta is as follows:
Generating effective orbital transfer thrust according to limited thrust orbital transfer gravity loss eta and engine thrustOrbital transfer engine thrust F capable of obtaining satellite with near-circular orbitGeneral assemblyComprises the following steps:
according to the limited thrust orbital transfer gravity loss eta and the total propellant demand M of satellite orbital transferGeneral assemblyCan obtain the propellant loss M of the near-circular orbit satelliteLoss of powerComprises the following steps:
Mloss of power=η×MGeneral assembly
In the formula, MGeneral assemblyThe total propellant demand for the satellite orbital transfer.
A thrust computing device can compute the thrust F of a orbital transfer engine of a near-circular orbit satelliteGeneral assemblyThe orbit eccentricity ratio of the near-circular orbit satellite is not more than 0.0025, and the orbit height is not more than 2000 km. Thrust F of orbital transfer engine of near-circular orbit satelliteGeneral assemblyComprises the following steps:
wherein
Wherein eta is the limited thrust orbital transfer gravity loss,generating effective orbital transfer thrust for engine thrust;is the ignition arc segment length.
Example (b):
a method for calculating the limited thrust orbital transfer gravity loss is disclosed, and the specific formula is deduced as follows.
The gravity loss is defined by the formula:
in the formula, eta is the gravity loss of the satellite orbit changing; fGeneral assemblyThe unit N is the thrust of the orbital transfer engine;the thrust value of the effective orbital transfer generated in the thrust of the engine is N.
According to the characteristics of the near-circular orbit, see fig. 1. The effective thrust magnitude is known as:
in the formula, Δ T is the total ignition time of the satellite in units of s; alpha is the included angle between the effective thrust and the total thrust in the orbital transfer process, and is the variable quantity, namely unit rad.
According to the above two formulas, the total thrust FGeneral assemblyFixed, the formula for the gravity loss calculation can be replaced by:
according to the characteristics of the circular orbit, the following relationship holds in the ignition process:
α=θ
in the formula, θ is represented by unit rad in FIG. 1.
Therefore, the temperature of the molten metal is controlled,
by integral transformation, can obtain
In the formula (I), the compound is shown in the specification,is the ignition arc length, in units rad.
In conclusion, the calculation formula of the gravity loss of the orbital transfer of the near-circular orbit satellite is obtained by the following steps:
the formula is only related to the propulsion arc length of the transition.
Generating effective orbital transfer thrust according to limited thrust orbital transfer gravity loss eta and engine thrustOrbital transfer engine thrust F capable of obtaining satellite with near-circular orbitGeneral assemblyComprises the following steps:
according to the limited thrust orbital transfer gravity loss eta and the total propellant demand M of satellite orbital transferGeneral assemblyCan obtain the propellant loss M of the near-circular orbit satelliteLoss of powerComprises the following steps:
Mloss of power=η×MGeneral assembly
In the formula, MGeneral assemblyThe total propellant demand for the satellite orbital transfer.
Those skilled in the art will appreciate that those matters not described in detail in the present specification are well known in the art.
Claims (7)
1. A method for calculating the gravity loss of limited thrust orbital transfer is suitable for a near-circular orbit satellite to adopt inertial directional orbital transfer and is characterized in that: the limited thrust orbital transfer gravity loss eta is as follows:
2. The method for calculating the limited thrust orbital transfer gravity loss according to claim 1, wherein the method comprises the following steps: orbital transfer engine thrust F of near-circular orbit satelliteGeneral assemblyComprises the following steps:
3. The method for calculating the limited thrust orbital transfer gravity loss according to claim 1, wherein the method comprises the following steps: propellant loss M of near-circular orbit satelliteLoss of powerComprises the following steps:
Mloss of power=η×MGeneral assembly
In the formula, MGeneral assemblyThe total propellant demand for the satellite orbital transfer.
4. The method for calculating the limited thrust orbital transfer gravity loss according to claim 1, wherein the method comprises the following steps: the orbit eccentricity ratio of the near-circular orbit satellite is not more than 0.0025, and the orbit height is not more than 2000 km.
5. A thrust force calculation device characterized in that: the device is used for calculating the orbital transfer engine thrust F of the near-circular orbit satelliteGeneral assemblyThrust F of orbital transfer engine of said satelliteGeneral assemblyComprises the following steps:
wherein
6. The thrust force calculation device according to claim 5, wherein: propellant loss M of near-circular orbit satelliteLoss of powerComprises the following steps:
Mloss of power=η×MGeneral assembly
In the formula, MGeneral assemblyThe total propellant demand for the satellite orbital transfer.
7. The thrust force calculation device according to claim 5, wherein: the orbit eccentricity ratio of the near-circular orbit satellite is not more than 0.0025, and the orbit height is not more than 2000 km.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811495695.XA CN109649689B (en) | 2018-12-07 | 2018-12-07 | Limited thrust orbital transfer gravity loss calculation method and thrust calculation device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811495695.XA CN109649689B (en) | 2018-12-07 | 2018-12-07 | Limited thrust orbital transfer gravity loss calculation method and thrust calculation device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109649689A CN109649689A (en) | 2019-04-19 |
CN109649689B true CN109649689B (en) | 2021-10-01 |
Family
ID=66113061
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811495695.XA Active CN109649689B (en) | 2018-12-07 | 2018-12-07 | Limited thrust orbital transfer gravity loss calculation method and thrust calculation device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109649689B (en) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5222023A (en) * | 1991-04-02 | 1993-06-22 | Space Systems/Loral, Inc. | Compensated transition for spacecraft attitude control |
CN101186236A (en) * | 2007-12-26 | 2008-05-28 | 北京控制工程研究所 | Orbit changing method for reducing space craft gravity loss |
US7918420B2 (en) * | 2007-07-17 | 2011-04-05 | The Boeing Company | System and methods for simultaneous momentum dumping and orbit control |
CN102874418A (en) * | 2012-10-24 | 2013-01-16 | 北京空间飞行器总体设计部 | Method for improving orbit-transferring safety of inclined orbit satellite |
CN103412563A (en) * | 2013-07-25 | 2013-11-27 | 北京控制工程研究所 | Method for estimating orbit control comprehensive efficiency and propellant consumption |
US9114890B1 (en) * | 2012-08-10 | 2015-08-25 | Lockheed Martin Corporation | Thruster orbit control method and configuration |
CN105620792A (en) * | 2016-02-05 | 2016-06-01 | 上海微小卫星工程中心 | Method for controlling attitude and orbit of satellite by adopting obliquely-arranged thrusters |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2001234479A1 (en) * | 2000-01-19 | 2001-07-31 | Kline And Walker, Llc | Protected accountable primary focal node interface |
US7967255B2 (en) * | 2006-07-27 | 2011-06-28 | Raytheon Company | Autonomous space flight system and planetary lander for executing a discrete landing sequence to remove unknown navigation error, perform hazard avoidance and relocate the lander and method |
EP2673195B1 (en) * | 2011-04-01 | 2017-02-01 | Geryon Space Technologies, Limited | Multi-body dynamics method of generating fuel efficient transfer orbits for spacecraft |
US10005568B2 (en) * | 2015-11-13 | 2018-06-26 | The Boeing Company | Energy efficient satellite maneuvering |
CN106114909B (en) * | 2016-06-23 | 2018-08-31 | 中国空间技术研究院 | A kind of Satellite Orbit Maneuver propellant waste computational methods |
-
2018
- 2018-12-07 CN CN201811495695.XA patent/CN109649689B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5222023A (en) * | 1991-04-02 | 1993-06-22 | Space Systems/Loral, Inc. | Compensated transition for spacecraft attitude control |
US7918420B2 (en) * | 2007-07-17 | 2011-04-05 | The Boeing Company | System and methods for simultaneous momentum dumping and orbit control |
CN101186236A (en) * | 2007-12-26 | 2008-05-28 | 北京控制工程研究所 | Orbit changing method for reducing space craft gravity loss |
US9114890B1 (en) * | 2012-08-10 | 2015-08-25 | Lockheed Martin Corporation | Thruster orbit control method and configuration |
CN102874418A (en) * | 2012-10-24 | 2013-01-16 | 北京空间飞行器总体设计部 | Method for improving orbit-transferring safety of inclined orbit satellite |
CN103412563A (en) * | 2013-07-25 | 2013-11-27 | 北京控制工程研究所 | Method for estimating orbit control comprehensive efficiency and propellant consumption |
CN105620792A (en) * | 2016-02-05 | 2016-06-01 | 上海微小卫星工程中心 | Method for controlling attitude and orbit of satellite by adopting obliquely-arranged thrusters |
Non-Patent Citations (2)
Title |
---|
A novel single thruster control strategy for spacecraft attitude stabilization;Godard;《ACTA ASTRONAUTICA》;20130615;全文 * |
卫星轨道圆化的点火控制策略;荆武兴;《宇航学报》;19971031;第18卷(第4期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109649689A (en) | 2019-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2021036778A1 (en) | Method for directly planning reentry trajectory in altitude-velocity profile | |
CN108759565B (en) | A kind of carrier rocket grade return phase precise guidance method based on virtual proportional guidance | |
CN107609267B (en) | A kind of moon Finite Thrust repeatedly captures track implementation method | |
CN103486905B (en) | Determining method for terminal guidance shift-exchange conditions of reenter vehicle | |
CN113602532B (en) | Solid carrier rocket in-orbit correction method | |
CN105501467A (en) | Carrier rocket exhaustion shutdown control method based on zero-range orientation | |
CN103412563A (en) | Method for estimating orbit control comprehensive efficiency and propellant consumption | |
CN104932517A (en) | Underwater unmanned aircraft track tracking dynamic surface control optimization method under environment interference | |
CN109269504B (en) | Attitude maneuver path planning method with terminal constraint | |
CN108511908A (en) | A kind of satellite antenna automatic following control system and method inhibiting function with phase | |
CN111591470A (en) | Aircraft accurate soft landing closed-loop guidance method adapting to thrust adjustable mode | |
CN102288177A (en) | Strapdown system speed calculating method based on angular speed output | |
CN106813663B (en) | Inertial navigation data and satellite navigation data synchronization method | |
CN109649689B (en) | Limited thrust orbital transfer gravity loss calculation method and thrust calculation device | |
CN109484675B (en) | Spacecraft orbit-in control method by utilizing space vector matching | |
CN111504140A (en) | Remote solid rocket instantaneous orbit landing point closed-loop optimal feedback guidance method | |
CN102749846B (en) | Design method of double parallel configuration VSDGCMGs singularity avoidance steering law | |
CN106021835B (en) | A kind of flight path design method towards optimal scouting | |
CN104375512B (en) | Spacecraft maneuvering path optimization method based on spectral analysis | |
CN113447025B (en) | Inertial navigation high-precision attitude angle resolving method and system based on Krafft angle | |
CN104848857A (en) | Method for automatically distributing accuracy indexes of ballistic missile inertia measurement system | |
CN110044384B (en) | Double-step filtering method suitable for vehicle-mounted transfer alignment | |
CN112455725A (en) | Method for transferring and converting pulse orbit transfer direction to limited thrust orbit | |
CN113447024A (en) | Inertial navigation attitude angle resolving method and system based on extended Krafft angle | |
CN106021726B (en) | It is a kind of transmitting observation satellite to Frozen Orbit eccentricity control method |
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 |