CN111483619A - Spacecraft gravitational acceleration calculation method and orbit control method - Google Patents

Spacecraft gravitational acceleration calculation method and orbit control method Download PDF

Info

Publication number
CN111483619A
CN111483619A CN202010309856.2A CN202010309856A CN111483619A CN 111483619 A CN111483619 A CN 111483619A CN 202010309856 A CN202010309856 A CN 202010309856A CN 111483619 A CN111483619 A CN 111483619A
Authority
CN
China
Prior art keywords
spacecraft
earth
gravitational acceleration
acceleration calculation
geocentric
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
CN202010309856.2A
Other languages
Chinese (zh)
Other versions
CN111483619B (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.)
Shanghai Engineering Center for Microsatellites
Innovation Academy for Microsatellites of CAS
Original Assignee
Shanghai Engineering Center for Microsatellites
Innovation Academy for Microsatellites of CAS
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 Shanghai Engineering Center for Microsatellites, Innovation Academy for Microsatellites of CAS filed Critical Shanghai Engineering Center for Microsatellites
Priority to CN202010309856.2A priority Critical patent/CN111483619B/en
Publication of CN111483619A publication Critical patent/CN111483619A/en
Application granted granted Critical
Publication of CN111483619B publication Critical patent/CN111483619B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/242Orbits and trajectories

Abstract

The invention provides a spacecraft gravitational acceleration calculation method and an orbit control method, and provides a spacecraft gravitational acceleration calculation method for eliminating singularity, aiming at the situation that the spacecraft gravitational acceleration calculated by a traditional method is singular. The dynamic model is the basis of spacecraft orbit design and control, wherein the earth gravity is an important part for forming the near-earth orbit dynamic model. When the gravity acceleration of the spacecraft is calculated by using the traditional calculation method, a singular condition exists when the geocentric latitude of the spacecraft is close to 90 degrees. In order to enhance the robustness of the algorithm and avoid the occurrence of singular conditions, a new calculation method is provided.

Description

Spacecraft gravitational acceleration calculation method and orbit control method
Technical Field
The invention relates to the technical field of spacecraft power, in particular to a spacecraft gravitational acceleration calculation method and a spacecraft orbit control method.
Background
The dynamic model is the basis of spacecraft orbit design and control, wherein the earth gravity is an important part for forming the near-earth orbit dynamic model. Since the earth is an ellipsoid, the gravitational force of the earth has non-spherical gravitational components in addition to the central gravitational force. The gravitational force is typically determined using a spherical harmonic model.
The spherical harmonic expansion of the gravitational potential of the earth is
Figure BDA0002457268780000011
Wherein: GM (GM)eIs the constant of the earth's gravity; reIs the earth mean radius; r, phi and lambda are geocentric distance, geocentric latitude and geocentric longitude respectively;
Figure BDA0002457268780000012
and
Figure BDA0002457268780000013
is a completely normalized spherical harmonic coefficient;
Figure BDA0002457268780000014
is a fully normalized associative legendre function.
The acceleration of the earth gravity can be obtained by calculating the first derivative of the earth gravity position to the position, but the spherical harmonic expansion formula of the earth gravity position takes the geocentric height r, the geocentric longitude phi and the geocentric latitude lambda under the earth-fixed relation as parameters, and usually an orbit integral equation is described by rectangular coordinates under the earth-fixed relation, so that the derivative of the earth gravity position to the position of the spacecraft under the earth-fixed relation needs to be calculated firstly, and then the derivative is converted into the earth-fixed relation. Then
Figure BDA0002457268780000015
Wherein:
Figure BDA0002457268780000021
Figure BDA0002457268780000022
Figure BDA0002457268780000023
and
Figure BDA0002457268780000024
the recurrence relation of (c) is as follows:
Figure BDA0002457268780000025
Figure BDA0002457268780000026
Figure BDA0002457268780000027
ijis Kronecker operator.
According to the expression and the spherical harmonic coefficient of the earth gravity field model, the earth gravity acceleration under the earth-solid relation can be calculated, but when the geocentric latitude of the spacecraft is close to 90 degrees, the expression is singular, so that a novel calculation method is provided for enhancing the robustness of the algorithm and avoiding the occurrence of the singular condition.
Disclosure of Invention
The invention aims to provide a spacecraft gravitational acceleration calculation method and an orbit control method, and aims to solve the problem that spacecraft gravitational acceleration under earth-fixed connection calculated by the existing method is singular.
In order to solve the technical problem, the invention provides a spacecraft gravitational acceleration calculation method, which relates to a parameter acquisition module, a spacecraft gravitational acceleration calculation module and a near-earth orbit dynamics model generation module, wherein:
the parameter acquisition module acquires the spherical harmonic coefficient of the earth gravitational field
Figure BDA0002457268780000031
And
Figure BDA0002457268780000032
transformation matrix from earth inertia system to earth solid relation
Figure BDA0002457268780000033
And the position r of the spacecraft in the inertial system of the earthECIAnd sending the acceleration to the spacecraft gravitational acceleration calculation module;
the spacecraft gravitational acceleration calculation module converts the earth inertial system into the earth-solid system according to the transformation matrix
Figure BDA0002457268780000034
And the position r of the spacecraft in the inertial system of the earthECIAnd calculating the position of the spacecraft under the earth-fixed connection:
Figure BDA0002457268780000035
the spacecraft gravitational acceleration calculation module calculates the position r according to the earth contactECEFCalculating the geocentric height r, the geocentric latitude phi and the geocentric longitude lambda of the spacecraft under earth-fixed connection;
calculating a first auxiliary variable by the spacecraft gravitational acceleration calculation module
Figure BDA0002457268780000036
The spacecraft gravitational acceleration calculation module calculates the spherical harmonic coefficient according to the earth gravitational field
Figure BDA0002457268780000037
And
Figure BDA0002457268780000038
geocentric altitude r, geocentric latitude phi, geocentric longitude lambda and first auxiliary variable of spacecraft in earth-solid connection
Figure BDA0002457268780000039
Calculating the gravitational acceleration a of the spacecraft under the earth-solid connectionECEF
The spacecraft gravitational acceleration calculation module converts the earth inertial system into the earth-solid system according to the transformation matrix
Figure BDA00024572687800000310
And the gravitational acceleration a of the spacecraft under the earth-solid connectionECEFAnd calculating the gravitational acceleration of the spacecraft under the earth inertial system
Figure BDA00024572687800000311
And generating a near-earth orbit dynamic model by the near-earth orbit dynamic model generating module according to the gravitational acceleration of the spacecraft under the earth inertial system to serve as orbit determination navigation data of the spacecraft.
Optionally, in the method for calculating gravitational acceleration of spacecraft, the spacecraft gravitational acceleration calculation module calculates the gravitational acceleration of spacecraft according to a position r under earth contactECEFCalculating the geocentric altitude r, the geocentric latitude phi and the geocentric longitude lambda of the spacecraft under the earth-fixed connection comprises the following steps:
Figure BDA00024572687800000312
wherein x, y and z are positions r of the spacecraft in earth-fixed connectionECEFThree-dimensional coordinates of (1).
Optionally, in the method for calculating gravitational acceleration of spacecraft, the calculation module for gravitational acceleration of spacecraft may calculate the height of the earth's center of the spacecraft in earth-fixed connectionrThe geocentric latitude phi and the geocentric longitude lambda are calculated to calculate the earth-solid relationGravitational acceleration a of spacecraftECEFThe method comprises the following steps:
Figure BDA0002457268780000041
wherein:
Figure BDA0002457268780000042
Figure BDA0002457268780000043
Figure BDA0002457268780000044
optionally, in the method for calculating gravitational acceleration of spacecraft, the spacecraft gravitational acceleration calculation module calculates a first auxiliary variable
Figure BDA0002457268780000045
The method comprises the following steps:
introducing said first auxiliary variable
Figure BDA0002457268780000046
The first auxiliary variable
Figure BDA0002457268780000047
And associated legendre function
Figure BDA0002457268780000048
The relationship of (a) to (b) is as follows:
Figure BDA0002457268780000049
the first auxiliary variable
Figure BDA00024572687800000410
Satisfies the following recursion relational expression:
Figure BDA00024572687800000411
Figure BDA00024572687800000412
optionally, in the method for calculating gravitational acceleration of spacecraft, the second auxiliary variable cmAnd a third auxiliary variable smThe specific definition of (A) is as follows:
Figure BDA0002457268780000051
optionally, in the method for calculating gravitational acceleration of spacecraft, the second auxiliary variable c is calculatedmAnd the third auxiliary variable smAnd the recurrence relation satisfies:
Figure BDA0002457268780000052
optionally, in the method for calculating gravitational acceleration of spacecraft, the earth gravitational field spherical harmonic coefficient
Figure BDA0002457268780000053
And
Figure BDA0002457268780000054
and obtaining by inquiring the earth gravity field model.
The invention also provides a spacecraft orbit control method, which comprises the spacecraft gravitational acceleration calculation method.
In the spacecraft gravitational acceleration computing method and the orbit control method provided by the invention, an earth gravitational field spherical harmonic coefficient, a conversion matrix from an earth inertia system to an earth solid relation and a position of a spacecraft under the earth inertia system are obtained by a parameter obtaining module and are sent to a spacecraft gravitational acceleration computing module, the spacecraft gravitational acceleration computing module computes the position of the spacecraft under the earth solid relation according to the conversion matrix from the earth inertia system to the earth solid relation and the position of the spacecraft under the earth inertia system, the spacecraft gravitational acceleration computing module computes the geocentric altitude, the geocentric latitude and the geocentric longitude of the spacecraft under the earth solid relation according to the position under the earth solid relation, the spacecraft gravitational acceleration computing module computes the geocentric altitude, the geocentric latitude and the geocentric longitude of the spacecraft under the earth solid relation according to the spacecraft gravitational field spherical harmonic coefficient, the geocentric altitude, the geocentric latitude and the geocentric longitude, the spacecraft gravitational acceleration calculation method based on the earth-fixed relation can be used for calculating the spacecraft gravitational acceleration under the earth-fixed relation according to a transformation matrix from the earth inertial system to the earth-fixed relation and the spacecraft gravitational acceleration under the earth-fixed relation.
Drawings
Fig. 1 is a schematic flow chart of a spacecraft gravitational acceleration calculation method according to an embodiment of the present invention.
Detailed Description
The spacecraft gravitational acceleration calculation method and the orbit control method provided by the invention are further described in detail below with reference to the accompanying drawings and specific embodiments. Advantages and features of the present invention will become apparent from the following description and from the claims. It is to be noted that the drawings are in a very simplified form and are not to precise scale, which is merely for the purpose of facilitating and distinctly claiming the embodiments of the present invention.
The core idea of the invention is to provide a spacecraft gravitational acceleration calculation method and an orbit control method, so as to solve the problem that the spacecraft gravitational acceleration calculated by the existing method under earth-fixed connection is singular.
In order to realize the above thought, the invention provides a spacecraft gravitational acceleration calculation method and an orbit control method, wherein the spacecraft gravitational acceleration calculation method relates to a parameter acquisition module and a spacecraft gravitational acceleration calculation module, wherein: the parameter acquisition module acquiresSpherical harmonic coefficient of spherical gravitational field
Figure BDA0002457268780000061
And
Figure BDA0002457268780000062
transformation matrix from earth inertia system to earth solid relation
Figure BDA0002457268780000063
And the position r of the spacecraft in the inertial system of the earthECIAnd sending the acceleration to the spacecraft gravitational acceleration calculation module; the spacecraft gravitational acceleration calculation module converts the earth inertial system into the earth-solid system according to a conversion matrix
Figure BDA0002457268780000064
And the position r of the spacecraft in the inertial system of the earthECIAnd calculating the position of the spacecraft under the earth-fixed connection:
Figure BDA0002457268780000065
the spacecraft gravitational acceleration calculation module is used for calculating the gravitational acceleration of the spacecraft according to the position r under the earth contactECEFCalculating the geocentric height r, the geocentric latitude phi and the geocentric longitude lambda of the spacecraft under earth-fixed connection; the spacecraft gravitational acceleration calculation module calculates a first auxiliary variable
Figure BDA0002457268780000066
The spacecraft gravitational acceleration calculation module calculates the spherical harmonic coefficient according to the spacecraft gravitational field
Figure BDA0002457268780000067
And
Figure BDA0002457268780000068
geocentric altitude r, geocentric latitude phi, geocentric longitude lambda and first auxiliary variable of spacecraft in earth-solid connection
Figure BDA0002457268780000069
Spacecraft under earth-fixed relation calculationAcceleration of gravity aECEF(ii) a The spacecraft gravitational acceleration calculation module converts the earth inertial system into the earth-solid system according to a conversion matrix
Figure BDA00024572687800000610
And the gravitational acceleration a of the spacecraft under the earth-solid connectionECEFAnd calculating the gravitational acceleration of the spacecraft under the earth inertial system
Figure BDA00024572687800000611
And the near-earth orbit dynamics model generation module generates a near-earth orbit dynamics model according to the gravitational acceleration of the spacecraft under the earth inertial system, and the near-earth orbit dynamics model is used as orbit determination navigation data of the spacecraft.
< example one >
The embodiment provides a spacecraft gravitational acceleration calculation method, as shown in fig. 1, the spacecraft gravitational acceleration calculation method relates to a spacecraft orbit control system, that is, the spacecraft gravitational acceleration calculation method is completed by the spacecraft orbit control system, the spacecraft orbit control system includes a parameter acquisition module and a spacecraft gravitational acceleration calculation module, wherein: the parameter acquisition module acquires the spherical harmonic coefficient of the earth gravitational field
Figure BDA00024572687800000612
And
Figure BDA00024572687800000613
transformation matrix from earth inertia system to earth solid relation
Figure BDA0002457268780000071
And the position r of the spacecraft in the inertial system of the earthECIAnd sending the acceleration to the spacecraft gravitational acceleration calculation module; the spacecraft gravitational acceleration calculation module converts the earth inertial system into the earth-solid system according to a conversion matrix
Figure BDA0002457268780000072
And the position r of the spacecraft in the inertial system of the earthECIAnd calculating the position of the spacecraft under the earth-fixed connection:
Figure BDA0002457268780000073
the spacecraft gravitational acceleration calculation module is used for calculating the gravitational acceleration of the spacecraft according to the position r under the earth contactECEFCalculating the geocentric height r, the geocentric latitude phi and the geocentric longitude lambda of the spacecraft under earth-fixed connection; the spacecraft gravitational acceleration calculation module calculates a first auxiliary variable
Figure BDA0002457268780000074
The spacecraft gravitational acceleration calculation module calculates the spherical harmonic coefficient according to the earth gravitational field
Figure BDA0002457268780000075
And
Figure BDA0002457268780000076
geocentric altitude r, geocentric latitude phi, geocentric longitude lambda and first auxiliary variable of spacecraft in earth-solid connection
Figure BDA0002457268780000077
Calculating the gravitational acceleration a of the spacecraft under the earth-solid connectionECEF(ii) a The spacecraft gravitational acceleration calculation module converts the earth inertial system into the earth-solid system according to a conversion matrix
Figure BDA0002457268780000078
And the gravitational acceleration a of the spacecraft under the earth-solid connectionECEFAnd calculating the gravitational acceleration of the spacecraft under the earth inertial system
Figure BDA0002457268780000079
And the near-earth orbit dynamics model generation module generates a near-earth orbit dynamics model according to the gravitational acceleration of the spacecraft under the earth inertial system, and the near-earth orbit dynamics model is used as orbit determination navigation data of the spacecraft.
Specifically, in the spacecraft orbit control system, the spacecraft gravitational acceleration calculation module calculates the position r according to the earth contactECEFCalculating the geocentric altitude r, the geocentric latitude phi and the geocentric longitude lambda of the spacecraft under the earth-fixed connection comprises the following steps:
Figure BDA00024572687800000710
wherein x, y and z are positions r of the spacecraft in earth-fixed connectionECEFThree-dimensional coordinates of (1).
Further, in the spacecraft orbit control system, the spacecraft gravitational acceleration calculation module calculates the spacecraft gravitational acceleration a under earth-solid connection according to the geocentric altitude r, the geocentric latitude phi and the geocentric longitude lambda of the spacecraft under earth-solid connectionECEFThe method comprises the following steps:
Figure BDA0002457268780000081
wherein:
Figure BDA0002457268780000082
Figure BDA0002457268780000083
Figure BDA0002457268780000084
in addition, in the spacecraft orbit control system, the spacecraft gravitational acceleration calculation module calculates a first auxiliary variable
Figure BDA0002457268780000085
The method comprises the following steps: introducing said first auxiliary variable
Figure BDA0002457268780000086
The first auxiliary variable
Figure BDA0002457268780000087
And associated legendre function
Figure BDA0002457268780000088
The relationship of (a) to (b) is as follows:
Figure BDA0002457268780000089
the first auxiliary variable
Figure BDA00024572687800000810
Satisfies the following recursion relational expression:
Figure BDA00024572687800000811
Figure BDA00024572687800000812
in particular, in the spacecraft orbit control system, the second auxiliary variable cmAnd a third auxiliary variable smThe specific definition is as follows:
Figure BDA00024572687800000813
further, in the spacecraft orbit control system, the second auxiliary variable c is calculatedmAnd the third auxiliary variable smAnd the recurrence relation satisfies:
Figure BDA0002457268780000091
in addition, in the spacecraft orbit control system, the earth gravitational field spherical harmonic coefficient
Figure BDA0002457268780000092
And
Figure BDA0002457268780000093
and obtaining by inquiring the earth gravity field model.
In summary, the above embodiments have described in detail different configurations of the spacecraft gravitational acceleration calculation method, and it is understood that the present invention includes, but is not limited to, the configurations listed in the above embodiments, and any modifications made on the configurations provided in the above embodiments are within the scope of the present invention. One skilled in the art can take the contents of the above embodiments to take a counter-measure.
< example two >
The embodiment provides a spacecraft orbit control method based on a spacecraft orbit control system in the previous embodiment, where the spacecraft orbit control method includes a spacecraft gravitational acceleration calculation method in the previous embodiment, and specifically includes: the parameter acquisition module acquires the spherical harmonic coefficient of the earth gravitational field
Figure BDA0002457268780000094
And
Figure BDA0002457268780000095
transformation matrix from earth inertia system to earth solid relation
Figure BDA0002457268780000096
And the position r of the spacecraft in the inertial system of the earthECIAnd sending the acceleration to a spacecraft gravitational acceleration calculation module; the spacecraft gravitational acceleration calculation module converts the earth inertial system into the earth-solid system according to a conversion matrix
Figure BDA0002457268780000097
And the position r of the spacecraft in the inertial system of the earthECIAnd calculating the position of the spacecraft under the earth-fixed connection:
Figure BDA0002457268780000098
the spacecraft gravitational acceleration calculation module is used for calculating the gravitational acceleration of the spacecraft according to the position r under the earth contactECEFCalculating the geocentric height r, the geocentric latitude phi and the geocentric longitude lambda of the spacecraft under earth-fixed connection; the above-mentionedThe spacecraft gravitational acceleration calculation module calculates a first auxiliary variable
Figure BDA0002457268780000099
The spacecraft gravitational acceleration calculation module calculates the spherical harmonic coefficient according to the earth gravitational field
Figure BDA00024572687800000910
And
Figure BDA00024572687800000911
earth center height of spacecraft in earth-fixed connectionrGeocentric latitude phi, geocentric longitude lambda and first auxiliary variable
Figure BDA00024572687800000912
Calculating the gravitational acceleration a of the spacecraft under the earth-solid connectionECEF(ii) a The spacecraft gravitational acceleration calculation module converts the earth inertial system into the earth-solid system according to a conversion matrix
Figure BDA00024572687800000913
And the gravitational acceleration a of the spacecraft under the earth-solid connectionECEFAnd calculating the gravitational acceleration of the spacecraft under the earth inertial system
Figure BDA00024572687800000914
And the near-earth orbit dynamics model generation module generates a near-earth orbit dynamics model according to the gravitational acceleration of the spacecraft under the earth inertial system and uses the near-earth orbit dynamics model as orbit determination navigation data of the spacecraft.
In the spacecraft gravitational acceleration computing method and the orbit control method provided by the invention, an earth gravitational field spherical harmonic coefficient, a conversion matrix from an earth inertia system to an earth solid relation and a position of a spacecraft under the earth inertia system are obtained by a parameter obtaining module and are sent to a spacecraft gravitational acceleration computing module, the spacecraft gravitational acceleration computing module computes the position of the spacecraft under the earth solid relation according to the conversion matrix from the earth inertia system to the earth solid relation and the position of the spacecraft under the earth inertia system, the spacecraft gravitational acceleration computing module computes the geocentric altitude, the geocentric latitude and the geocentric longitude of the spacecraft under the earth solid relation according to the position under the earth solid relation, the spacecraft gravitational acceleration computing module computes the geocentric altitude, the geocentric latitude and the geocentric longitude of the spacecraft under the earth solid relation according to the earth gravitational field spherical harmonic coefficient, the spacecraft gravitational acceleration calculation method based on the earth-fixed relation can be used for calculating the spacecraft gravitational acceleration under the earth-fixed relation according to a transformation matrix from the earth inertial system to the earth-fixed relation and the spacecraft gravitational acceleration under the earth-fixed relation.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is relatively simple because the system corresponds to the method disclosed by the embodiment, and the relevant points can be referred to the method part for description.
The above description is only for the purpose of describing the preferred embodiments of the present invention, and is not intended to limit the scope of the present invention, and any variations and modifications made by those skilled in the art based on the above disclosure are within the scope of the appended claims.

Claims (8)

1. A spacecraft gravitational acceleration calculation method is characterized by involving a parameter acquisition module, a spacecraft gravitational acceleration calculation module and a near-earth orbit dynamics model generation module, wherein:
the parameter acquisition module acquires the spherical harmonic coefficient of the earth gravitational field
Figure FDA0002457268770000011
And
Figure FDA0002457268770000012
transformation matrix from earth inertia system to earth solid relation
Figure FDA0002457268770000013
And the position r of the spacecraft in the inertial system of the earthECIAnd sending the acceleration to the spacecraft gravitational acceleration calculation module;
the spacecraft gravitational acceleration calculation module converts the earth inertial system into the earth-solid system according to the transformation matrix
Figure FDA0002457268770000014
And the position r of the spacecraft in the inertial system of the earthECIAnd calculating the position of the spacecraft under the earth-fixed connection:
Figure FDA0002457268770000015
the spacecraft gravitational acceleration calculation module calculates the position r according to the earth contactECEFCalculating the geocentric height r, the geocentric latitude phi and the geocentric longitude lambda of the spacecraft under earth-fixed connection;
calculating a first auxiliary variable by the spacecraft gravitational acceleration calculation module
Figure FDA0002457268770000016
The spacecraft gravitational acceleration calculation module calculates the spherical harmonic coefficient according to the earth gravitational field
Figure FDA0002457268770000017
And
Figure FDA0002457268770000018
geocentric altitude r, geocentric latitude phi, geocentric longitude lambda and first auxiliary variable of spacecraft in earth-solid connection
Figure FDA0002457268770000019
Calculating the navigation under earth-fixed connectionAcceleration of gravity of celestial body aECEF
The spacecraft gravitational acceleration calculation module converts the earth inertial system into the earth-solid system according to the transformation matrix
Figure FDA00024572687700000110
And the gravitational acceleration a of the spacecraft under the earth-solid connectionECEFAnd calculating the gravitational acceleration of the spacecraft under the earth inertial system
Figure FDA00024572687700000111
And generating a near-earth orbit dynamic model by the near-earth orbit dynamic model generating module according to the gravitational acceleration of the spacecraft under the earth inertial system to serve as orbit determination navigation data of the spacecraft.
2. A spacecraft gravitational acceleration calculation method according to claim 1, wherein said spacecraft gravitational acceleration calculation module is based on a position r under earth contactECEFCalculating the geocentric altitude r, the geocentric latitude phi and the geocentric longitude lambda of the spacecraft under the earth-fixed connection comprises the following steps:
Figure FDA0002457268770000021
wherein x, y and z are positions r of the spacecraft in earth-fixed connectionECEFThree-dimensional coordinates of (1).
3. A spacecraft gravitational acceleration computing method according to claim 2, wherein said spacecraft gravitational acceleration computing module computes spacecraft gravitational acceleration a in earth-solid relation according to geocentric altitude r, geocentric latitude Φ and geocentric longitude λ of spacecraft in earth-solid relationECEFThe method comprises the following steps:
Figure FDA0002457268770000022
wherein:
Figure FDA0002457268770000023
Figure FDA0002457268770000024
Figure FDA0002457268770000025
4. a spacecraft gravitational acceleration calculation method according to claim 3, wherein said spacecraft gravitational acceleration calculation module calculates a first auxiliary variable
Figure FDA0002457268770000026
The method comprises the following steps:
introducing said first auxiliary variable
Figure FDA0002457268770000027
The first auxiliary variable
Figure FDA0002457268770000028
And associated legendre function
Figure FDA0002457268770000029
The relationship of (a) to (b) is as follows:
Figure FDA00024572687700000210
the first auxiliary variable
Figure FDA00024572687700000211
Satisfies the following recursion relational expression:
Figure FDA0002457268770000031
Figure FDA0002457268770000032
5. a spacecraft gravitational acceleration calculation method according to claim 4, characterized in that the second auxiliary variable cmAnd a third auxiliary variable smThe specific definition of (A) is as follows:
Figure FDA0002457268770000033
6. a spacecraft gravitational acceleration calculation method according to claim 5, characterized in that said second auxiliary variable c is calculatedmAnd the third auxiliary variable smAnd the recurrence relation satisfies:
Figure FDA0002457268770000034
7. a spacecraft gravitational acceleration calculation method according to claim 1, characterized in that said earth gravitational field spherical harmonic coefficients
Figure FDA0002457268770000035
And
Figure FDA0002457268770000036
and obtaining by inquiring the earth gravity field model.
8. A spacecraft orbit control method, characterized by comprising the spacecraft gravitational acceleration calculation method according to any one of claims 1 to 7.
CN202010309856.2A 2020-04-20 2020-04-20 Spacecraft gravitational acceleration calculation method and orbit control method Active CN111483619B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010309856.2A CN111483619B (en) 2020-04-20 2020-04-20 Spacecraft gravitational acceleration calculation method and orbit control method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010309856.2A CN111483619B (en) 2020-04-20 2020-04-20 Spacecraft gravitational acceleration calculation method and orbit control method

Publications (2)

Publication Number Publication Date
CN111483619A true CN111483619A (en) 2020-08-04
CN111483619B CN111483619B (en) 2021-07-23

Family

ID=71811915

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010309856.2A Active CN111483619B (en) 2020-04-20 2020-04-20 Spacecraft gravitational acceleration calculation method and orbit control method

Country Status (1)

Country Link
CN (1) CN111483619B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2768557C1 (en) * 2021-03-23 2022-03-24 Федеральное Государственное Унитарное Предприятие "Всероссийский Научно-Исследовательский Институт Физико-Технических И Радиотехнических Измерений" (Фгуп "Вниифтри") Method for measuring the gravitational acceleration of a space vehicle

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5393017A (en) * 1993-01-25 1995-02-28 Lockheed Missiles & Space Company, Inc. Technique for dispensing earth satellites into multi-planar orbits
CN103091723A (en) * 2013-02-06 2013-05-08 中国科学院测量与地球物理研究所 Method of reducing influences of gravity satellite centroid adjustment errors to earth gravitational field accuracy
CN105739537A (en) * 2016-03-29 2016-07-06 北京理工大学 Active control method for adhesion motion on small celestial body surface
CN106202640A (en) * 2016-06-28 2016-12-07 西北工业大学 Day ground three body gravitational fields in halo orbit spacecraft bias track method for designing
CN110826224A (en) * 2019-11-06 2020-02-21 北京理工大学 Method for determining spherical harmonic coefficient of small celestial body gravitational field based on gravitational acceleration

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5393017A (en) * 1993-01-25 1995-02-28 Lockheed Missiles & Space Company, Inc. Technique for dispensing earth satellites into multi-planar orbits
CN103091723A (en) * 2013-02-06 2013-05-08 中国科学院测量与地球物理研究所 Method of reducing influences of gravity satellite centroid adjustment errors to earth gravitational field accuracy
CN105739537A (en) * 2016-03-29 2016-07-06 北京理工大学 Active control method for adhesion motion on small celestial body surface
CN106202640A (en) * 2016-06-28 2016-12-07 西北工业大学 Day ground three body gravitational fields in halo orbit spacecraft bias track method for designing
CN110826224A (en) * 2019-11-06 2020-02-21 北京理工大学 Method for determining spherical harmonic coefficient of small celestial body gravitational field based on gravitational acceleration

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
柳培清: "GPS卫星精密定轨中的地球引力加速度研究", 《测绘科学》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2768557C1 (en) * 2021-03-23 2022-03-24 Федеральное Государственное Унитарное Предприятие "Всероссийский Научно-Исследовательский Институт Физико-Технических И Радиотехнических Измерений" (Фгуп "Вниифтри") Method for measuring the gravitational acceleration of a space vehicle

Also Published As

Publication number Publication date
CN111483619B (en) 2021-07-23

Similar Documents

Publication Publication Date Title
CN111323050B (en) Strapdown inertial navigation and Doppler combined system calibration method
CN106871928B (en) Strap-down inertial navigation initial alignment method based on lie group filtering
Wang et al. Estimation of information sharing error by dynamic deformation between inertial navigation systems
Pascoal et al. Navigation system design using time-varying complementary filters
CN104165640A (en) Near-space missile-borne strap-down inertial navigation system transfer alignment method based on star sensor
CN113819906A (en) Combined navigation robust filtering method based on statistical similarity measurement
CN102944241B (en) Spacecraft relative attitude determining method based on multicell liner differential inclusion
CN105203104A (en) Gravity field modeling method suitable for high-precision inertial navigation system
CN108427427A (en) A kind of spacecraft is to earth&#39;s surface oriented object posture angle computational methods
CN108225325B (en) Polar region transverse navigation method based on virtual ball model
CN111366984B (en) Method for determining gravitational field model based on gravity satellite inter-satellite laser ranging system
CN106885577B (en) Autonomous orbit determination method for Lagrange navigation satellite
CN111380518A (en) SINS/USBL tight combination navigation positioning method introducing radial velocity
CN108592943A (en) A kind of inertial system coarse alignment computational methods based on OPREQ methods
CN111552003A (en) Asteroid gravitational field full-autonomous measurement system and method based on ball satellite formation
CN103424127A (en) Method for transfer alignment of speed and specific force matching
CN111483619B (en) Spacecraft gravitational acceleration calculation method and orbit control method
Ding et al. A new contour-based combined matching algorithm for underwater terrain-aided strapdown inertial navigation system
CN112179334A (en) Star navigation method and system based on two-step Kalman filtering
CN111207773A (en) Attitude unconstrained optimization solving method for bionic polarized light navigation
CN110554443A (en) Method for determining earth gravity field based on carrier phase observed value and point acceleration method
CN112393835B (en) Small satellite on-orbit thrust calibration method based on extended Kalman filtering
CN109799704A (en) Flexible spacecraft attitude and structure vibration coupling control method
CN113108787B (en) Long-endurance inertial navigation/satellite global integrated navigation method
CN106643726B (en) Unified inertial navigation resolving 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