CN111483619B - 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
CN111483619B
CN111483619B CN202010309856.2A CN202010309856A CN111483619B CN 111483619 B CN111483619 B CN 111483619B CN 202010309856 A CN202010309856 A CN 202010309856A CN 111483619 B CN111483619 B CN 111483619B
Authority
CN
China
Prior art keywords
spacecraft
earth
gravitational acceleration
geocentric
acceleration calculation
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
Application number
CN202010309856.2A
Other languages
Chinese (zh)
Other versions
CN111483619A (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 respectivelyGeocentric distance, geocentric latitude, and geocentric longitude;
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 connectionECEFTo calculate the groundGravitational acceleration of spacecraft under ball 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 gravitational acceleration a of the spacecraft under the earth-solid connectionECEFThe 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 acquires the spherical harmonic coefficient of the earth 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 sent to the spacecraft for gravitational accelerationA degree 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
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 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
The near-earth orbit dynamics model generation moduleAnd generating a near-earth orbit dynamic model according to the gravitational acceleration of the spacecraft under the earth inertial system to serve 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 controlIn the 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 above-mentionedFirst 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 spacecraft gravitational acceleration calculation module calculates a first auxiliary variable
Figure BDA0002457268780000099
The spacecraft gravitational acceleration calculation module is used for calculating the spherical harmonic of the gravitational field of the earthCoefficient of performance
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 (7)

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 FDA0002997333580000011
And
Figure FDA0002997333580000012
transformation matrix from earth inertia system to earth solid relation
Figure FDA0002997333580000013
And spacecraft in the earth's inertial systemPosition rECIAnd 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 FDA0002997333580000014
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 FDA0002997333580000015
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 FDA0002997333580000016
The spacecraft gravitational acceleration calculation module calculates the spherical harmonic coefficient according to the earth gravitational field
Figure FDA0002997333580000017
And
Figure FDA0002997333580000018
geocentric altitude r, geocentric latitude phi, geocentric longitude lambda and first auxiliary variable of spacecraft in earth-solid connection
Figure FDA0002997333580000019
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 FDA00029973335800000110
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 FDA00029973335800000111
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 FDA0002997333580000021
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 FDA0002997333580000022
wherein:
Figure FDA0002997333580000023
Figure FDA0002997333580000024
Figure FDA0002997333580000025
u is a spherical harmonic expansion formula of the earth gravitational potential; GM (GM)eIs the constant of the earth's gravity; delta0mIs a Kronecker operator;
second auxiliary variable cmAnd a third auxiliary variable smThe specific definition of (A) is as follows:
Figure FDA0002997333580000026
4. a spacecraft gravitational acceleration calculation method according to claim 3, wherein said spacecraft gravitational acceleration calculation module calculates a first auxiliary variable
Figure FDA0002997333580000027
The method comprises the following steps:
introducing said first auxiliary variable
Figure FDA0002997333580000028
The first auxiliary variable
Figure FDA0002997333580000029
And associated legendre function
Figure FDA00029973335800000210
The relationship of (a) to (b) is as follows:
Figure FDA0002997333580000031
the first auxiliary variable
Figure FDA0002997333580000032
Satisfies the following recursion relational expression:
Figure FDA0002997333580000033
Figure FDA0002997333580000034
5. a spacecraft gravitational acceleration calculation method according to claim 4, characterized in that said second auxiliary variable c is calculatedmAnd the third auxiliary variable smAnd the recurrence relation satisfies:
Figure FDA0002997333580000035
c0=1,s0=0
6. a spacecraft gravitational acceleration calculation method according to claim 1, characterized in that said earth gravitational field spherical harmonic coefficients
Figure FDA0002997333580000036
And
Figure FDA0002997333580000037
and obtaining by inquiring the earth gravity field model.
7. A spacecraft orbit control method, characterized by comprising the spacecraft gravitational acceleration calculation method according to any one of claims 1 to 6.
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 CN111483619A (en) 2020-08-04
CN111483619B true 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)

Families Citing this family (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

Family Cites Families (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
CN103091723B (en) * 2013-02-06 2015-05-13 中国科学院测量与地球物理研究所 Method of reducing influences of gravity satellite centroid adjustment errors to earth gravitational field accuracy
CN105739537B (en) * 2016-03-29 2018-06-19 北京理工大学 A kind of small feature loss surface attachment movement Active Control Method
CN106202640B (en) * 2016-06-28 2018-03-27 西北工业大学 Day ground three body gravitational fields in halo orbit spacecraft bias track design method
CN110826224A (en) * 2019-11-06 2020-02-21 北京理工大学 Method for determining spherical harmonic coefficient of small celestial body gravitational field based on gravitational acceleration

Also Published As

Publication number Publication date
CN111483619A (en) 2020-08-04

Similar Documents

Publication Publication Date Title
Wang et al. Estimation of information sharing error by dynamic deformation between inertial navigation systems
CN113819906A (en) Combined navigation robust filtering method based on statistical similarity measurement
US20080258971A1 (en) Method and System for Controlling the Direction of an Antenna Beam
KR100898169B1 (en) Initial alignment method of inertial navigation system
CN111323050A (en) Strapdown inertial navigation and Doppler combined system calibration method
CN105203104A (en) Gravity field modeling method suitable for high-precision inertial navigation system
CN111483619B (en) Spacecraft gravitational acceleration calculation method and orbit control method
CN103424127A (en) Method for transfer alignment of speed and specific force matching
CN106885577B (en) Autonomous orbit determination method for Lagrange navigation satellite
CN108592943A (en) A kind of inertial system coarse alignment computational methods based on OPREQ methods
CN115451952B (en) Multi-system integrated navigation method and device for fault detection and robust adaptive filtering
CN111552003A (en) Asteroid gravitational field full-autonomous measurement system and method based on ball satellite formation
CN112179334A (en) Star navigation method and system based on two-step Kalman filtering
CN104121930B (en) A kind of compensation method based on the MEMS gyro drift error adding table coupling
CN111207773A (en) Attitude unconstrained optimization solving method for bionic polarized light navigation
Liu et al. A method for SINS alignment with large initial misalignment angles based on Kalman filter with parameters resetting
CN112284412B (en) Ground static alignment method for avoiding precision reduction caused by singular Euler transformation
CN110044385B (en) Rapid transfer alignment method under condition of large misalignment angle
CN112393835B (en) Small satellite on-orbit thrust calibration method based on extended Kalman filtering
CN113108787B (en) Long-endurance inertial navigation/satellite global integrated navigation method
CN115031729A (en) SINS/DVL/USBL underwater tight combination navigation method and device and underwater carrier control equipment
KR101665375B1 (en) Navigation system and method
CN111238473B (en) Second-order damping method for height channel of inertial navigation system under geocentric geostationary coordinate system
CN115829879B (en) Attitude quaternion processing method, device and equipment for agile satellite
CN111896027B (en) Distance measuring sensor simulation modeling method considering topographic relief

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