CN113602534B - On-orbit calibration method for magnitude of micro electric propulsion thrust - Google Patents
On-orbit calibration method for magnitude of micro electric propulsion thrust Download PDFInfo
- Publication number
- CN113602534B CN113602534B CN202110715137.5A CN202110715137A CN113602534B CN 113602534 B CN113602534 B CN 113602534B CN 202110715137 A CN202110715137 A CN 202110715137A CN 113602534 B CN113602534 B CN 113602534B
- Authority
- CN
- China
- Prior art keywords
- orbit
- perturbation
- satellite
- thrust
- rail
- 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
- 238000000034 method Methods 0.000 title claims abstract description 71
- 230000008569 process Effects 0.000 claims abstract description 46
- 230000001133 acceleration Effects 0.000 claims abstract description 26
- 230000008859 change Effects 0.000 claims abstract description 18
- 230000001174 ascending effect Effects 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims abstract description 6
- 230000005484 gravity Effects 0.000 claims description 9
- 238000004458 analytical method Methods 0.000 claims description 5
- 150000001875 compounds Chemical class 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
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/40—Arrangements or adaptations of propulsion systems
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Remote Sensing (AREA)
- Aviation & Aerospace Engineering (AREA)
- Navigation (AREA)
Abstract
The invention discloses an on-orbit calibration method for the magnitude of micro electric propulsion thrust, which comprises the following steps: calculating the change of the semi-long axis in the rail lifting and lowering processes; theoretically analyzing the influence of earth spherical perturbation and atmospheric resistance perturbation in the rail ascending and descending processes; calculating a semi-major axis variable quantity formula obtained by calculation in the process of simultaneous rail lifting and rail lowering to obtain thrust acceleration; and then the electric thrust is obtained according to the thrust acceleration and the satellite mass. The thrust size does not need to be obtained through attitude or angular speed change conversion, namely the thruster does not need to be eccentrically installed to generate eccentric torque to change the attitude or the angular speed, the influence of spherical perturbation and atmospheric resistance perturbation which are relatively close to the magnitude order of the thrust of the miniature thruster on the calibration of the thrust size is eliminated, and the on-orbit calibration of the thrust size of a single thruster passing through the mass center can be realized.
Description
Technical Field
The invention belongs to the field of on-orbit calibration of electric propulsion thrust, and particularly relates to an on-orbit calibration method of micro electric propulsion thrust.
Background
The micro electric propulsion has the advantages of small volume, light weight, large specific impulse, high total impulse and the like, and is more and more widely applied in the field of commercial aerospace in recent years. However, the general thrust of the micro electric propulsion product is tens or hundreds of micro-newtons, and the temperature and pressure conditions consistent with the space environment are difficult to establish on the ground due to the limitation of the ground environment; and the accurate calibration of the thrust of the miniature electric propeller is difficult to realize due to the limitation of ground measurement means. In order to ensure the effectiveness of the in-orbit operation of the miniature electric thruster, the thrust of the miniature electric thruster needs to be calibrated in the in-orbit mode.
The existing thrust calibration methods comprise two methods, one is an orbit calibration method, and the orbit calibration method calculates to obtain the satellite speed increment according to the orbit change parameters; and calculating to obtain the thrust of the propulsion system according to the satellite speed increment and the electric propulsion working time. The method does not consider the influence of space perturbation force such as earth spherical perturbation and atmospheric resistance perturbation on thrust calibration. The other method is an attitude calibration method, in the method, a thruster which is not beyond a centroid is started to generate a control torque, and the thrust is calculated through the change of the satellite attitude and the angular speed. But the method is only suitable for installing a plurality of thrusters, each thruster is not installed through the mass center of the satellite, and the method cannot meet the on-orbit calibration of the thrust of a single thruster installed through the mass center.
Disclosure of Invention
Aiming at the defects of the existing micro electric thruster thrust magnitude on-orbit calibration technology, the invention provides an on-orbit calibration method of the micro electric thruster thrust magnitude, namely a lifting rail combined calibration method, which considers the influences of earth spherical perturbation and atmospheric resistance perturbation and can realize the on-orbit calibration of the thrust magnitude of a single mass center passing mounting thruster.
The invention is realized by adopting the following technical scheme: an on-orbit calibration method for the magnitude of micro electric propulsion thrust comprises the following steps:
step 1, calculating semimajor axis change in a rail lifting process;
a t is the thrust acceleration to which the satellite is subjected during the orbit raising process, a d1 For acceleration of lift rail drag, a e1 For the perturbed acceleration of the earth's spherical shape during the rail-lifting process, t 1 Total thrust application time for the rail lifting process, n 1 The total number of the running tracks of the satellite in the orbit raising process, T is the average orbit period in the orbit raising process, delta a is the variation of the semi-major axis of the orbit in the orbit raising process, a represents the semi-major axis of the average orbit, and mu =398600km 3 /s 2 Is an earth gravity parameter;
step 2, calculating the change of the semimajor axis in the rail descending process;
a t thrust acceleration to which the satellite is subjected during the course of falling into orbit, a d2 For falling rail resistance acceleration, a e2 For global perturbation acceleration, t, during the rail lowering process 2 For the total thrust application time in the rail lowering process, n 2 The total number of the running orbits of the satellite in the orbit descending process is T, the average orbit period in the orbit descending process is T, and delta a is the variation of the orbit semi-major axis in the orbit descending process;
step 3, theoretically analyzing the influence of earth spherical perturbation and atmospheric resistance perturbation in the rail ascending and descending processes;
1) For a low earth orbit satellite, the earth gravity perturbation mainly considers J 2 Determining that the J2 perturbation does not cause the change of the semi-major axis by integrating one circle along the satellite orbit;
2) During the process of raising and lowering the same orbit height, the atmospheric resistance suffered by the satellite is equal in magnitude and same in direction;
step 4, further according to the theoretical analysis of the step 3, equations (3) and (5) are combined to obtain the magnitude of the thrust acceleration;
step 5, calculating the electric thrust F according to the thrust acceleration and the satellite mass T The following are:
wherein m is the satellite mass.
Further, in step 3, a specific theoretical analysis process is as follows:
(1) The Gaussian perturbation equation of the track number shows that perturbation forces influencing the semi-major axis of the track are radial force and tangential force, and are shown as follows:
wherein a is the semi-major axis of the track, e is the eccentricity, and theta is the true paraxial point angle,
is the average velocity, F r 、F t 、F n The external forces of the satellite in the radial direction, the tangential direction and the normal direction are respectively considered, and J is mainly considered for the perturbation of the earth gravity of the satellite in the near-earth orbit 2 Influence of item perturbation, J 2 The components of the perturbation force in the radial direction, the tangential direction and the normal direction are as follows:
in the formula, J 2 ＝1.08263×10 -3 Is J 2 Coefficient of term R e =6378.14km for the radius of the earth, r = a (1-e)/(1 + ecos θ) for the centroid distance, u = ω + θ for the latitude argument, the perturbation force is substituted into the gaussian equation, and one revolution is integrated along the orbit to find J 2 Item perturbation does not cause a change in the semi-major axis, so J per complete turn during long run 2 Item perturbation all self-counteracts, incomplete rounds of J 2 The item perturbation influence is ignored;
(2) Because thin high-rise atmosphere exists on the near-earth orbit, the satellite runs under atmospheric resistance for a long time in orbit, the orbit height is gradually attenuated, and the perturbation of the atmospheric resistance and the aerodynamic coefficient C of the satellite surface d The atmospheric density rho, the windward area S and the satellite velocity v are related, and the equation is as follows:
because the space atmosphere model and the satellite surface aerodynamic coefficients cannot be accurately estimated, the atmospheric resistance received by the satellite cannot be accurately calculated, but the atmospheric resistance received by the satellite is considered to be equal in size and same in direction in the process of raising and lowering the same orbit height, and therefore the influence of the atmospheric resistance is counteracted through the simultaneous calculation of subtraction of the lifting orbit thrust equation.
Compared with the prior art, the invention has the advantages and positive effects that:
the calibration method provided by the scheme eliminates the influence of spherical perturbation and atmospheric resistance perturbation on the calibration of the thrust magnitude, wherein the magnitude order of the thrust magnitude of the micro thruster is relatively close to that of the spherical perturbation and the atmospheric resistance perturbation; and the on-orbit calibration of the thrust magnitude of a single through-center-of-mass mounting thruster can be realized.
Drawings
FIG. 1 is a schematic diagram of a combined calibration process of a lifting rail for micro electric propulsion thrust in an embodiment of the present invention;
FIG. 2 is a schematic diagram of a variation curve of a semi-major axis of a satellite orbit with time when a micro electric propulsion device is used to complete orbit lifting according to an embodiment of the present invention.
Detailed Description
In order to make the above objects, features and advantages of the present invention more clearly understood, the present invention will be further described with reference to the accompanying drawings and examples. In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention, however, the present invention may be practiced in other ways than those described herein, and thus, the present invention is not limited to the specific embodiments disclosed below.
The method for calibrating the magnitude of the micro electric propulsion thrust in the on-orbit requires a satellite to perform 2 stages of operation of ascending and descending the orbit, eliminates the influence of perturbation force by selecting the part with the consistent height of the orbit of the ascending and descending orbit, and realizes the accurate on-orbit calibration of the magnitude of the micro electric propulsion thrust.
The on-orbit calibration method for the magnitude of the micro electric propulsion thrust provided by the invention obtains the following calculation formula for the magnitude of the micro electric propulsion thrust:
in the formula, F T Representing micro-electric propulsion thrust, m representing satellite mass, n 1 Indicates the total orbit number n of the satellite in the orbit-ascending stage 2 Represents the total orbit number t of the satellite in the orbit reduction stage 1 Representing the total thrust application time, t, in the rail lifting phase 2 The total application time of the thrust in the rail descending stage is shown, delta a represents the variation of the semi-major axis of the track, a represents the semi-major axis of the average track, and mu =398600km 3 /s 2 Is the gravity parameter.
Specifically, as shown in fig. 1, a schematic diagram of a combined calibration flow of a lifting rail for micro electric propulsion thrust is shown, the objective of the invention is to accurately calibrate the magnitude of the micro electric propulsion thrust, and consider the influences of spherical perturbation and atmospheric resistance perturbation in the calibration process, and simultaneously realize the on-rail calibration of the magnitude of the thrust of a single over-center-of-mass mounted thruster, and the specific implementation steps of the combined calibration of the lifting rail for micro electric propulsion thrust are as follows:
step 1, calculating the semimajor axis change in the rail lifting process according to external force and time;
suppose the thrust acceleration of the satellite during the orbit raising process is a t Acceleration of resistance a d1 The earth's spherical perturbation acceleration is a e1 . Recording the total thrust application time t during the rail lifting process 1 Total number of orbits n of satellite 1 The average track period T and the variation Δ a of the track semimajor axis are as follows:
a t ·t 1 +a e1 ·n 1 ·T-a d1 ·n 1 ·T＝Δv (2)
in the formula (I), the compound is shown in the specification,
the velocity increment generated when the satellite is subjected to external force in the process of orbit rising.
Thus obtaining:
step 2, calculating the change of the semimajor axis in the rail descending process according to the external force and time;
assuming that the thrust acceleration received by the satellite in the process of orbit reduction is a t Acceleration of resistance a d2 The earth's spherical perturbation acceleration is a e2 . Recording the total application time t of the thrust during the rail descending process 2 Total number of orbits n of satellite operation 2 The average track period T and the variation Δ a of the track semi-major axis are
a t ·t 2 +a e2 ·n 2 ·T+a d2 ·n 2 ·T＝Δv (4)
Step 3, theoretically analyzing the influence of earth spherical perturbation and atmospheric resistance perturbation in the rail ascending and descending processes;
the Gaussian perturbation equation of the orbit root can know that perturbation forces influencing the semi-major axis of the orbit are radial force and tangential force, and the following is shown:
wherein a is the semi-major axis of the track, e is the eccentricity, theta is the true paraxial point angle,
to average speed, F r 、F t 、F n The external forces in the radial direction, the tangential direction and the normal direction are respectively applied to the satellite. For a low earth orbit satellite, the earth gravity perturbation mainly considers J 2 Influence of item perturbation, J 2 The components of the item perturbation force in the radial direction, the tangential direction and the normal direction are as follows:
in the formula, J 2 ＝1.082×63 -3 Is J 2 Coefficient of term, μ =398600km 3 /s 2 As a parameter of Earth's gravity, R e =6378.14km for earth radius, r = a (1-e)/(1 + ecos θ) for geocentric distance, and u = ω + θ for latitude argument. The perturbation force is substituted into the Gaussian equation, and the J can be found by integrating one circle along the track 2 Item perturbation does not cause variation in the semi-major axis. Thus, J for each complete revolution during long run 2 The item perturbation can be automatically counteracted, and J of incomplete circle 2 Item perturbation effects are negligible in small quantities.
Due to the thin high atmosphere on the low earth orbit, the satellite runs under atmospheric resistance for a long time in orbit, and the orbit height gradually attenuates. Atmospheric drag perturbation and satellite surface aerodynamic coefficient C d The atmospheric density ρ, the windward area S and the satellite velocity v are related, and the equation is as follows:
because the space atmosphere model and the surface aerodynamic coefficients of the satellite cannot be accurately estimated, the atmospheric resistance of the satellite cannot be accurately calculated. However, the atmospheric resistance experienced by the satellite during the same orbital altitude elevation and lowering can be considered equal and equal. Therefore, the influence of the atmospheric resistance can be offset by the subtraction simultaneous calculation of the lifting rail thrust equation.
Step 4, simultaneous equations (3) and (5) are used for solving the magnitude of the thrust acceleration;
and 5, calculating the electric thrust force according to the thrust acceleration and the satellite mass, wherein the electric thrust force comprises the following steps:
from the above formula, it can be seen that: the method does not need to obtain the thrust magnitude through attitude or angular speed change conversion, namely the thruster does not need to be eccentrically installed to generate eccentric torque to change the attitude or the angular speed, so that the method can realize the on-orbit calibration of the thrust magnitude of a single thruster passing through the center of mass.
The above method was verified by simulation as follows:
and setting the satellite mass to be 38kg and the electric thrust to be 500 mu N, and establishing a high-precision space environment ground simulation model. Firstly, the average semi-major axis of the initial orbit of the satellite is set to be 6878.5km, the thruster is opened for 150 tracks, each track is opened for 40 minutes, and the average semi-major axis of the orbit which can be reached by the satellite is known to be 6886.75km through simulation (see the curve ascending segment of fig. 2).
Then, the thruster is started at 146 tracks, each track is started for 40 minutes, and the simulation finds that the satellite can return to the initial track (see the curve descending segment of fig. 2) with the average semi-major axis of the track being 6878.5km, and the thrust of the micro electric propulsion is 498 mu N through calculation of the formula (1), and is basically consistent with the set value. The analysis can show that the proposed method can effectively calibrate the magnitude of the micro electric propulsion thrust.
The above description is only a preferred embodiment of the present invention, and not intended to limit the present invention in other forms, and any person skilled in the art may apply the above modifications or changes to the equivalent embodiments with equivalent changes, without departing from the technical spirit of the present invention, and any simple modification, equivalent change and change made to the above embodiments according to the technical spirit of the present invention still belong to the protection scope of the technical spirit of the present invention.
Claims (1)
1. An on-orbit calibration method for the magnitude of micro electric propulsion thrust is characterized by comprising the following steps:
step 1, calculating semimajor axis change in a rail lifting process;
a t is the thrust acceleration to which the satellite is subjected during the orbit raising process, a d1 For acceleration of lift rail drag, a e1 For the perturbed acceleration of the earth's spherical shape during the rail-lifting process, t 1 Total thrust application time for the rail lifting process, n 1 Is the total orbit number of the satellite in the orbit raising process, T is the average orbit period in the orbit raising process, delta a is the variation of the orbit semi-major axis in the orbit raising process,
denotes the mean orbit semimajor axis, μ =398600km 3 /s 2 Is an earth gravity parameter;
step 2, calculating the semimajor axis change in the rail descending process;
a t thrust and acceleration to which the satellite is subjected during the descent of the orbit, a d2 For falling rail resistance acceleration, a e2 For the acceleration of the earth's spherical perturbation during the rail lowering process, t 2 For the total thrust application time in the rail lowering process, n 2 The total number of the running orbits of the satellite in the orbit descending process is T, the average orbit period in the orbit descending process is T, and delta a is the variation of the orbit semi-major axis in the orbit descending process;
step 3, theoretically analyzing the influence of earth spherical perturbation and atmospheric resistance perturbation in the rail ascending and descending processes;
1) For a low earth orbit satellite, the earth gravity perturbation mainly considers J 2 Determining that the J2 perturbation does not cause the change of the semi-major axis by integrating one circle along the satellite orbit;
2) During the process of raising and lowering the same orbit height, the satellite receives the same atmospheric resistance in the same direction;
step 4, according to the theoretical analysis of the step 3, equations (3) and (5) are combined to obtain the magnitude of the thrust acceleration;
step 5, calculating the electric thrust F according to the thrust acceleration and the satellite mass T The following are:
wherein m is the satellite mass;
in the step 3, the specific theoretical analysis process is as follows:
(1) The Gaussian perturbation equation of the orbit root can know that perturbation forces influencing the semi-major axis of the orbit are radial force and tangential force, and the following is shown:
wherein a is the semi-major axis of the track, e is the eccentricity, theta is the true paraxial point angle,
is the average velocity, F r 、F t 、F n The external forces of the satellite in the radial direction, the tangential direction and the normal direction are respectively, and the earth gravity perturbation mainly considers J for the low earth orbit satellite 2 Influence of item perturbation, J 2 The components of the perturbation force in the radial direction, the tangential direction and the normal direction are as follows:
in the formula (I), the compound is shown in the specification,
is J 2 Coefficient of term, R e =6378.14km as earth radius, r = a (1-e)/(1 + e cos theta) as geocentric distance, u = ω + theta as latitude argument, the perturbation force is substituted into a Gaussian equation, and one circle is integrated along the track to find J 2 Item perturbation does not cause variation in the semi-major axis, so J is every complete revolution during long run 2 Item perturbation all self-counteracts, incomplete rounds of J 2 The item perturbation influence is ignored; (2) Because of the thin high atmosphere on the near-earth orbit, the satellite moves on orbit for a long time under the atmospheric resistance, the orbit height gradually attenuates, the atmospheric resistance perturbation and the aerodynamic coefficient C of the satellite surface d The atmospheric density ρ, the windward area S and the satellite velocity v are related, and the equation is as follows:
because the space atmosphere model and the satellite surface aerodynamic coefficient can not be accurately estimated, the atmospheric resistance received by the satellite can not be accurately calculated, but the atmospheric resistance received by the satellite is considered to be equal in size and same in direction in the process of raising and lowering the same orbit height, and therefore the influence of the atmospheric resistance is counteracted through the simultaneous calculation of the subtraction of the lifting orbit thrust equation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110715137.5A CN113602534B (en) | 2021-06-26 | 2021-06-26 | On-orbit calibration method for magnitude of micro electric propulsion thrust |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110715137.5A CN113602534B (en) | 2021-06-26 | 2021-06-26 | On-orbit calibration method for magnitude of micro electric propulsion thrust |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113602534A CN113602534A (en) | 2021-11-05 |
| CN113602534B true CN113602534B (en) | 2023-02-28 |
Family
ID=78303742
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110715137.5A Active CN113602534B (en) | 2021-06-26 | 2021-06-26 | On-orbit calibration method for magnitude of micro electric propulsion thrust |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113602534B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113815902B (en) * | 2021-11-24 | 2022-02-15 | 北京航天驭星科技有限公司 | Method for calibrating satellite orbit control effect |
| CN114771873B (en) * | 2022-03-24 | 2024-05-03 | 北京控制工程研究所 | Autonomous accurate maintenance method for ultra-low orbit satellite orbit |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103940431A (en) * | 2014-04-11 | 2014-07-23 | 北京空间飞行器总体设计部 | Tangential low-thrust in-orbit circular orbit calibration method based on (Global Navigation Satellite System) GNSS precise orbit determination |
| CN105651516A (en) * | 2014-11-11 | 2016-06-08 | 航天恒星科技有限公司 | Engine thrust calibration method based on GNSS observation value and calibration device |
| CN106697333A (en) * | 2017-01-12 | 2017-05-24 | 北京理工大学 | Robustness analysis method for spacecraft orbit control strategy |
| CN109190155A (en) * | 2018-07-25 | 2019-01-11 | 西北工业大学 | A kind of continuous low-thrust trajectory design method of mixing promoted using electric propulsion/solar sail |
| CN109870260A (en) * | 2019-02-27 | 2019-06-11 | 北京航空航天大学 | A kind of method of on-line measurement MEMS solid micro-thruster array thrust output |
| CN111994304A (en) * | 2020-08-31 | 2020-11-27 | 北京理工大学 | Low-thrust long-term position keeping method for geostationary orbit satellite |
| CN112298614A (en) * | 2020-09-18 | 2021-02-02 | 中国人民解放军战略支援部队航天工程大学 | Thrust on-orbit calibration test method |
| CN112393835A (en) * | 2020-11-03 | 2021-02-23 | 西北工业大学深圳研究院 | Small satellite on-orbit thrust calibration method based on extended Kalman filtering |
| CN112455725A (en) * | 2020-11-27 | 2021-03-09 | 山东航天电子技术研究所 | Method for transferring and converting pulse orbit transfer direction to limited thrust orbit |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU7357198A (en) * | 1997-03-25 | 1998-10-20 | Edward A. Belbruno | Low energy method for changing the inclinations of orbiting satellites using weak stability boundaries and a computer process for implementing same |
| GB2384058B (en) * | 2002-01-15 | 2005-11-30 | Rolls Royce Plc | Thrust correction |
| US20130002484A1 (en) * | 2011-07-03 | 2013-01-03 | Daniel A. Katz | Indoor navigation with gnss receivers |
| US10352784B2 (en) * | 2013-08-26 | 2019-07-16 | University Of Florida Research Foundation, Incorporated | Method and apparatus for measuring thrust |
| FR3014082B1 (en) * | 2013-11-29 | 2016-01-01 | Thales Sa | TUYER SYSTEM AND METHOD FOR ORBIT AND ATTITUDE CONTROL FOR GEOSTATIONARY SATELLITE |
| US10427806B2 (en) * | 2016-06-15 | 2019-10-01 | The Aerospace Corporation | Deployment and control algorithms for wheel cluster formations of satellites |
| CN106379555A (en) * | 2016-09-05 | 2017-02-08 | 北京理工大学 | Optimal orbital transfer method of low-earth-orbit satellite under limited thrust by taking J2 perturbation into consideration |
| CN111591469B (en) * | 2020-03-03 | 2021-12-07 | 航天科工空间工程发展有限公司 | Low-orbit constellation system phase keeping method, system, equipment and storage medium |
-
2021
- 2021-06-26 CN CN202110715137.5A patent/CN113602534B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103940431A (en) * | 2014-04-11 | 2014-07-23 | 北京空间飞行器总体设计部 | Tangential low-thrust in-orbit circular orbit calibration method based on (Global Navigation Satellite System) GNSS precise orbit determination |
| CN105651516A (en) * | 2014-11-11 | 2016-06-08 | 航天恒星科技有限公司 | Engine thrust calibration method based on GNSS observation value and calibration device |
| CN106697333A (en) * | 2017-01-12 | 2017-05-24 | 北京理工大学 | Robustness analysis method for spacecraft orbit control strategy |
| CN109190155A (en) * | 2018-07-25 | 2019-01-11 | 西北工业大学 | A kind of continuous low-thrust trajectory design method of mixing promoted using electric propulsion/solar sail |
| CN109870260A (en) * | 2019-02-27 | 2019-06-11 | 北京航空航天大学 | A kind of method of on-line measurement MEMS solid micro-thruster array thrust output |
| CN111994304A (en) * | 2020-08-31 | 2020-11-27 | 北京理工大学 | Low-thrust long-term position keeping method for geostationary orbit satellite |
| CN112298614A (en) * | 2020-09-18 | 2021-02-02 | 中国人民解放军战略支援部队航天工程大学 | Thrust on-orbit calibration test method |
| CN112393835A (en) * | 2020-11-03 | 2021-02-23 | 西北工业大学深圳研究院 | Small satellite on-orbit thrust calibration method based on extended Kalman filtering |
| CN112455725A (en) * | 2020-11-27 | 2021-03-09 | 山东航天电子技术研究所 | Method for transferring and converting pulse orbit transfer direction to limited thrust orbit |
Non-Patent Citations (3)
| Title |
|---|
| 全电推进卫星小推力变轨策略星上计算方法研究;王敏，李强，梁新刚，安然;《推进技术》;20190827(第01期);全文 * |
| 微纳卫星相对轨道机动控制技术研究;孙书剑;《中国博士学位论文全文数据库工程科技II辑》;20210115(第1期);全文 * |
| 激光微推力器作用下纳星轨道半长轴简化计算方法;周伟静，叶继飞，常浩，李南雷;《红外与激光工程》;20190425;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113602534A (en) | 2021-11-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113602534B (en) | On-orbit calibration method for magnitude of micro electric propulsion thrust | |
| Prudden et al. | Measuring wind with small unmanned aircraft systems | |
| CN108548542B (en) | Near-earth orbit determination method based on atmospheric resistance acceleration measurement | |
| CN109710961B (en) | High-altitude unmanned aerial vehicle limit rising data processing method based on GPS data | |
| CN105573337B (en) | A kind of braking Closed Loop Guidance method that leaves the right or normal track for meeting reentry angle and voyage constraint | |
| CN112298614B (en) | Thrust on-orbit calibration test method | |
| CN106021784B (en) | A kind of full track mark optimum design method based on bilevel optimization strategy | |
| CN113343442B (en) | Method and system for solving fixed-time finite fuel multi-pulse transfer orbit | |
| CN112666960B (en) | L1 augmentation self-adaption-based control method for rotary wing aircraft | |
| CN108082538B (en) | Multi-body system low-energy track capturing method considering initial and final constraints | |
| CN114684389B (en) | Lunar transfer window considering reentry constraint and accurate transfer track determining method | |
| CN108225323A (en) | Determine to settle in an area method, medium and the equipment on boundary based on deviation effects directional combination | |
| Bui et al. | Flight research of an aerospike nozzle using high power solid rockets | |
| CN112393835B (en) | Small satellite on-orbit thrust calibration method based on extended Kalman filtering | |
| CN104932266A (en) | Precision control method for entering section of lander based on feed-forward compensation | |
| CN109781374A (en) | A kind of method that real-time online quickly estimates aircraft thrust | |
| CN114234910A (en) | Inertia and ADS height fusion method based on air pressure reference self-adaptive correction | |
| Davis et al. | X-43A flight-test-determined aerodynamic force and moment characteristics at Mach 7.0 | |
| CN107839903B (en) | Method for estimating air bleeding time of sailboard in transfer orbit section of single-wing satellite | |
| CN103884333B (en) | A kind of survey of deep space independent navigation initial baseline catching method | |
| CN107804487A (en) | A kind of great-jump-forward based on the control of adaptive deviation, which reenters, returns to impact prediction method | |
| CN113418499A (en) | Method and system for resolving roll angle of rotary aircraft | |
| CN114771873B (en) | Autonomous accurate maintenance method for ultra-low orbit satellite orbit | |
| Spencer et al. | Mars pathfinder atmospheric entry reconstruction | |
| Moyano Cano | Quadrotor UAV for wind profile characterization |
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 |