CN113627029A - Method for designing track near earth-moon triangle translational point - Google Patents
Method for designing track near earth-moon triangle translational point Download PDFInfo
- Publication number
- CN113627029A CN113627029A CN202110973579.XA CN202110973579A CN113627029A CN 113627029 A CN113627029 A CN 113627029A CN 202110973579 A CN202110973579 A CN 202110973579A CN 113627029 A CN113627029 A CN 113627029A
- Authority
- CN
- China
- Prior art keywords
- orbit
- moon
- earth
- under
- track
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 21
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 claims abstract description 19
- 238000013519 translation Methods 0.000 claims abstract description 16
- 238000013461 design Methods 0.000 claims abstract description 10
- 239000013598 vector Substances 0.000 claims abstract description 6
- 230000010354 integration Effects 0.000 claims description 13
- 238000006467 substitution reaction Methods 0.000 claims description 13
- 238000004091 panning Methods 0.000 claims description 10
- 208000037280 Trisomy Diseases 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 6
- 230000005484 gravity Effects 0.000 claims description 2
- 238000012216 screening Methods 0.000 claims description 2
- 239000000243 solution Substances 0.000 description 15
- 238000010586 diagram Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 238000012423 maintenance Methods 0.000 description 2
- 239000003637 basic solution Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/17—Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
-
- 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)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Operations Research (AREA)
- Complex Calculations (AREA)
Abstract
The invention discloses a method for designing a track near a moon triangle translation point, which comprises the following steps: (1) constructing a semi-analytical solution approximation of the motion near the moon trigonometric translation points L4 and L5 under the ephemeris model; (2) designing a task orbit which is maintained for a certain time near the earth-moon triangle translational points L4 and L5 under the ephemeris model; (3) and outputting a track design result, namely generating the position and the speed vector of the detector under the geocentric celestial coordinate system and the earth-moon instantaneous convergence coordinate system at certain time intervals. The invention can design the task track which is maintained for a long time from any time, the adopted semi-analytical solution has higher precision, and the invention is suitable for two triangular translation points and is easy for a track designer to directly call.
Description
Technical Field
The invention relates to the technical field of aerospace, in particular to a method for designing an orbit near a moon-earth triangle translation point.
Background
Theoretical research results under the existing ephemeris model are not complete, and a certain distance is reserved from the practical engineering application requirements. The simplified circular restrictive trisomy problem the orbit designed near the trigonometric pan point will usually diverge very quickly in a real lunar system. Taking the point L4 under the circular limiting trisomy problem as an example, under the simplified circular limiting trisomy problem model, the probe will be fixed at this point and remain unchanged, however, this orbit diverges quickly when operating in a real lunar system and is not suitable as a mission orbit. Briefly: and the track under the circular restrictive three-body problem model can not meet the requirement as the precision of the task track under the actual force model. The double-circle model of the earth-moon system is not dynamically self-consistent, and for collinear translational points, previous work has shown that the model causes distortion phenomena of certain dynamic structures. At present, no relevant research work exists on the triangle translational point, but the essential defect that the dynamics is not self-consistent still influences the applicability of the solution, so that the method has only theoretical research significance. The plane earth-moon-sun three-body model is more consistent with an actual physical system than a non-self-consistent double-circle model, but still causes a certain degree of distortion compared with the result under the ephemeris model. For the low-amplitude dynamics alternative orbit, the instability index given under the ephemeris model is 1 order of magnitude higher than the result under the planar earth-moon-sun trisomy model. This indicates that the instability of the small amplitude unstable orbit in the actual terrestrial system is worse than the result in the simplified planar model. Briefly: compared with the actual earth-moon system, the double circle model or the planar earth-moon-sun model has certain difference in dynamics, so that the designed orbit has qualitative difference compared with the actual operation orbit of the detector, and the long-time task orbit can be designed by taking the orbit as a basic solution and is limited.
Disclosure of Invention
The technical problem to be solved by the invention is to provide a method for designing a track near a moon-earth triangle translation point, which can design a task track maintained for a long time from any moment, has higher precision of a adopted semi-analytical solution, is suitable for two triangle translation points, and is easy for a track designer to directly call.
In order to solve the technical problem, the invention provides a method for designing a track near a moon-earth triangle translation point, which comprises the following steps:
(1) constructing a semi-analytical solution approximation of the motion near the moon trigonometric translation points L4 and L5 under the ephemeris model;
(2) designing a task orbit which is maintained for a certain time near the earth-moon triangle translational points L4 and L5 under the ephemeris model;
(3) and outputting a track design result, namely generating the position and the speed vector of the detector under the geocentric celestial coordinate system and the earth-moon instantaneous convergence coordinate system at certain time intervals.
Preferably, in the step (1), constructing a semi-analytic solution approximation of the motion near the earth-moon trigonometric translation points L4 and L5 under the ephemeris model specifically includes the following steps:
(11) establishing a motion equation under a geocentric celestial coordinate system;
(12) aiming at the small amplitude dynamics replacement orbit, the specific processing comprises the following steps:
(a) constructing a quasi-periodic orbit which is maintained for a certain time by adopting a numerical method of multi-point target shooting through a numerical iteration process, wherein the initial value of iteration is any orbit type near a triangular translation point of a circular restrictive three-body problem and needs small-amplitude motion, and before numerical improvement, the orbit is transferred from a convergence coordinate system X-Y-Z to a geocentric celestial coordinate system X-Y-Z;
(b) aiming at the quasi-periodic orbit constructed in the step (a), analyzing the frequencies by adopting a Fast Fourier Transform (FFT) mode, and only extracting specific frequency items and coefficients thereof aiming at the frequencies
Wherein ω is1~ω4Is four fundamental frequencies related to the relative movement of three, namely, sun, earth and moon, and the four fundamental frequencies are kept unchanged no matter what orbit is in the earth-moon system, and the extracted terms of the specific frequencies are represented together in the following form
Wherein
(1) Compared with the orbit under the circular restrictive trisomy problem, the formula is a better dynamic approximation of the real motion under the ephemeris model;
(c) repeating the steps (a) and (b) to gradually improve the frequency combination and the corresponding coefficient in the formula (1) by taking the formula as an initial guess value, and repeating the steps until the coefficient change amount corresponding to each frequency item is smaller than a preset value epsilon under a certain frequency item combination, namely, the coefficient change amount is satisfied
Wherein the superscript (m) represents the result of the mth iteration, and at this time, the small-amplitude dynamics substitution orbit is considered to be obtained under a certain precision condition;
(d) after the steps are carried out, a small-amplitude dynamics substitution track is obtained, and a semi-analysis description of the movement near the small-amplitude dynamics substitution track is given by analyzing the change rule of the deviation of the movement near the small-amplitude dynamics substitution track relative to the dynamics substitution track along with time;
(13) aiming at the large-amplitude dynamics alternative orbit, the specific processing comprises the following steps:
(a) firstly, the initial state quantity and the initial time t are measured under the geocentric convergence coordinate system x-y-z0Gridding in steps, i.e.
Wherein h is*Respectively representing the grid step lengths of corresponding quantities, and recording the set of all grid points as N0Wherein, an orbit is generated after any point is integrated, and a shorter integration time upper limit T is set0Go through set N0All tracks of (1), selected at integration time T0Then the points still can stay in the area near the triangle panning point, and the set of the points is recorded as N1It is obvious that
(b) For set N1Set a longer upper limit T of the integration time1Go through set N1All tracks of (1), selected at integration time T1Then the points still can stay in the area near the triangle panning point, and the set of the points is recorded as N2It is obvious that
(c) Repeating the step (b), gradually screening initial grid nodes by gradually increasing the upper limit of integration time, and finally obtaining nodes which still stay in the area near the triangular translational point within a million-year time scale as alternatives of a stable track;
(d) analyzing the frequency of the selected stable orbit through Fast Fourier Transform (FFT), and only keeping a term in the form of equation (1) as a semi-analytical solution approximation of the large-amplitude dynamics alternative orbit;
(e) like the case of the small-amplitude dynamics alternative orbit, a semi-analytical solution of the large-amplitude dynamics alternative orbit vicinity motion is constructed by analyzing the motion in its vicinity.
Preferably, in step (11), the relative motion is given by JPL DE 405 numerical ephemeris, taking into account the gravity of mass points of the earth, moon and sun; under the earth center celestial coordinate system, the motion equation of the small celestial body (detector) is
Wherein GmE,GmM,GmSGravitational constants of the earth, moon and sun, R, RM、RSRespectively are the position vectors of the small celestial body, the moon and the sun under the earth center celestial coordinate system.
Preferably, in the step (2), the step of designing the mission trajectory for maintaining a certain time around the earth-moon trigonometric translation points L4 and L5 under the ephemeris model specifically includes the following steps:
(21) aiming at the small amplitude dynamics replacement orbit, the specific processing comprises the following steps:
(a) selecting a proper time interval to generate a series of track nodes according to the constructed semi-analytical solution, wherein the total time interval of the whole track nodes needs to meet the task time requirement;
(b) combining possible task boundary constraints, adopting a multipoint target shooting method, and performing numerical improvement under an ephemeris model to generate a more accurate target orbit;
(22) and generating an initial orbit value and performing numerical integration on the initial epoch selected according to the constructed semi-analytical solution aiming at the large-amplitude dynamics alternative orbit so as to design a stable dynamics alternative orbit or an orbit nearby the stable dynamics alternative orbit.
The invention has the beneficial effects that: (1) the adopted ephemeris model is more accurate, so that the designed task orbit has higher precision and can meet the task precision requirement; (2) the high-precision semi-analytical solution is benefited, and a task track which is maintained for a long time can be designed from any time by combining a numerical method; (3) the adopted semi-analytical solution has higher precision and is popularized to two triangular translation points, and the solution is easy to be directly called by a rail designer through comparatively systematic arrangement.
Drawings
FIG. 1 is a schematic flow chart of the method of the present invention.
Fig. 2 is a schematic diagram of an orbit image of a task orbit near a small-amplitude dynamic substitution orbit near a point L4 in a geocentric celestial coordinate system.
Fig. 3 is a schematic diagram of an orbit image of a mission orbit near a small-amplitude dynamic substitution orbit near a point L4 under a geodetic instantaneous convergence coordinate system.
FIG. 4 is a schematic diagram of an orbit image of a large-amplitude dynamic alternative orbit near the point L5 in a geocentric celestial coordinate system.
Fig. 5 is a schematic diagram of a track image of a dynamic alternative track with large amplitude near the point L5 designed by the invention under a lunar instantaneous convergence coordinate system.
Detailed Description
As shown in fig. 1, a method for designing a track near a moon triangle panning point includes the following steps:
(1) constructing a semi-analytical solution approximation of the motion near the moon trigonometric translation points L4 and L5 under the ephemeris model;
(2) designing a task orbit which is maintained for a certain time near the earth-moon triangle translational points L4 and L5 under the ephemeris model;
(3) and outputting a track design result, namely generating the position and the speed vector of the detector under the geocentric celestial coordinate system and the earth-moon instantaneous convergence coordinate system at certain time intervals.
For the dynamic alternative orbit of L4 point or L5 point, whether small amplitude or large amplitude, it can be expressed as follows under the geocentric celestial coordinate system:
wherein
Frequency omega1~ω4Is related to the relative movement of the three parts of the sun, the earth and the moon and keeps unchanged.
If the trajectory is a dynamic substitution trajectory with small amplitude, the trajectory in the vicinity thereof (deviation from the trajectory with small amplitude) can be expressed as follows:
where α, β are the offset amplitude parameters in and out of the white-road plane of the relatively small amplitude dynamic surrogate orbit, θ1、θ2Respectively, the argument corresponding to the amplitude alpha, beta.
If the dynamic displacement orbit with large amplitude is adopted, the orbit (deviation of the orbit with relatively large amplitude) nearby can be expressed as follows:
where α, β are two amplitude parameters within the white track plane of the relatively large amplitude dynamic replacement trajectory for the trajectory, γ is an offset amplitude parameter outside the white track plane of the relatively large amplitude dynamic replacement trajectory for the trajectory, and θ1,θ2,θ3Respectively, the argument corresponding to the amplitude alpha, beta, gamma. The semi-analytical solution of the above form forms the basis for the overall track design.
Example 1:
assuming that the starting epoch of the mission trajectory is 2021 year 3 month 21 day 01 hour 01: 32.0 seconds (greenwich mean time), the maintenance time of the constructed trajectory is 10 years, and it is necessary to construct a small amplitude dynamics near L4 point in place of the trajectory near the trajectory, the amplitude parameter in the above equation is set to 0.005, and β is set to 0.005.
Selecting 915 track nodes in 10-year tracks with the time interval of 4 days according to the semi-analytical solution; and (4) carrying out numerical iteration to design a task orbit by adopting a multipoint shooting method under a complete ephemeris model.
And outputting the designed track node data for analysis and drawing. Fig. 2 and 3 show images of the designed trajectory (maintained for 10 years) in the geocentric celestial coordinate system and the earth-moon instantaneous convergence coordinate system (projection of the trajectory on x-y, x-z, y-z planes and image in three-dimensional space), respectively, according to the output data. In the geocentric celestial coordinate system, the height of the orbit near the L4 point is designed to be equivalent to the height of the lunar orbit, and the average angular velocity of motion is the same as the orbital angular velocity of the moon, so the trajectory appears as an orbit around the earth (fig. 2). Under the earth-moon instantaneous convergence coordinate system, the designed orbit near the point L4 needs to keep the characteristics of motion around the point L4 and does not diverge within the mission time (i.e. intersect with the earth-moon line or be more than 2 earth-moon distances away from the earth), so the trajectory is displayed as an orbit around the point L4 (FIG. 3). The third trajectory shown in fig. 2 is suitable for a space task staying near the point of march L4.
Example 2:
assuming that the start epoch of the mission trajectory is 0.0 minutes 15 minutes at 22 days 3 and 9 months 2025, the maintenance time of the constructed trajectory is 10 years, and a large-amplitude dynamic replacement trajectory itself near the point L5 needs to be constructed, the amplitude parameter in the above equation is set to 0.0, and β is set to 0.0. And according to a semi-analytical solution, giving the position and the speed of the initial epoch moment detector under the geocentric celestial coordinate system, and generating a 10-year task orbit by integration.
And outputting the designed track node data for analysis and drawing. Fig. 4 and 5 show images of the designed trajectory (maintained for 10 years) in the geocentric celestial coordinate system and the earth-moon instantaneous convergence coordinate system (projection of the trajectory on x-y, x-z, y-z planes and image in three-dimensional space), respectively, according to the output data. In the geocentric celestial coordinate system, the height of the orbit near the L5 point is designed to be equivalent to the height of the lunar orbit, and the average angular velocity of motion is the same as the orbital angular velocity of the moon, so the trajectory appears as an orbit around the earth (fig. 4). Under the earth-moon instantaneous convergence coordinate system, the designed orbit near the point L5 needs to keep the characteristics of motion around the point L5 and does not diverge within the mission time (i.e. intersect with the earth-moon line or be more than 2 earth-moon distances away from the earth), so the trajectory is displayed as an orbit around the point L5 (FIG. 5). The third trajectory shown in fig. 4 is suitable for a space task staying near the point of march L5.
Claims (4)
1. A method for designing a track near a moon-earth triangle translation point is characterized by comprising the following steps:
(1) constructing a semi-analytical solution approximation of the motion near the moon trigonometric translation points L4 and L5 under the ephemeris model;
(2) designing a task orbit which is maintained for a certain time near the earth-moon triangle translational points L4 and L5 under the ephemeris model;
(3) and outputting a track design result, namely generating the position and the speed vector of the detector under the geocentric celestial coordinate system and the earth-moon instantaneous convergence coordinate system at certain time intervals.
2. The method for designing an orbit near the moon triangle panning point as claimed in claim 1, wherein in the step (1), constructing a semi-analytic solution approximation of the motion near the moon triangle panning points L4 and L5 under the ephemeris model specifically comprises the following steps:
(11) establishing a motion equation under a geocentric celestial coordinate system;
(12) aiming at the small amplitude dynamics replacement orbit, the specific processing comprises the following steps:
(a) constructing a quasi-periodic orbit which is maintained for a certain time by adopting a numerical method of multi-point target shooting through a numerical iteration process, wherein the initial value of iteration is any orbit type near a triangular translation point of a circular restrictive three-body problem and needs small-amplitude motion, and before numerical improvement, the orbit is transferred from a convergence coordinate system X-Y-Z to a geocentric celestial coordinate system X-Y-Z;
(b) aiming at the quasi-periodic orbit constructed in the step (a), analyzing the frequencies by adopting a Fast Fourier Transform (FFT) mode, and only extracting specific frequency items and coefficients thereof aiming at the frequencies
Wherein ω is1~ω4Is four fundamental frequencies related to the relative movement of three, namely, sun, earth and moon, no matter what type of orbit in the earth-moon system, the four fundamental frequencies are kept unchanged, and the extracted terms of the specific frequencies are represented together in the following form
Wherein
(1) Compared with the orbit under the circular restrictive trisomy problem, the formula is a better dynamic approximation of the real motion under the ephemeris model;
(c) repeating the steps (a) and (b) to gradually improve the frequency combination and the corresponding coefficient in the formula (1) by taking the formula as an initial guess value, and repeating the steps until the coefficient change amount corresponding to each frequency item is smaller than a preset value epsilon under a certain frequency item combination, namely, the coefficient change amount is satisfied
Wherein the superscript (m) represents the result of the mth iteration, and at this time, the small-amplitude dynamics substitution orbit is considered to be obtained under a certain precision condition;
(d) after the steps are carried out, a small-amplitude dynamics substitution track is obtained, and a semi-analysis description of the movement near the small-amplitude dynamics substitution track is given by analyzing the change rule of the deviation of the movement near the small-amplitude dynamics substitution track relative to the dynamics substitution track along with time;
(13) aiming at the large-amplitude dynamics alternative orbit, the specific processing comprises the following steps:
(a) firstly, the initial state quantity and the initial time t are measured under the geocentric convergence coordinate system x-y-z0Gridding in steps, i.e.
Wherein h is*Respectively representing the grid step lengths of corresponding quantities, and recording the set of all grid points as N0Wherein, an orbit is generated after any point is integrated, and a shorter integration time upper limit T is set0Go through set N0All tracks of (1), selected at integration time T0Then the points still can stay in the area near the triangle panning point, and the set of the points is recorded as N1It is obvious that
(b) For set N1Set a longer upper limit T of the integration time1Go through set N1All tracks of (1), selected at integration time T1Then the points still can stay in the area near the triangle panning point, and the set of the points is recorded as N2It is obvious that
(c) Repeating the step (b), gradually screening initial grid nodes by gradually increasing the upper limit of integration time, and finally obtaining nodes which still stay in the area near the triangular translational point within a million-year time scale as alternatives of a stable track;
(d) analyzing the frequency of the selected stable orbit through Fast Fourier Transform (FFT), and only keeping a term in the form of equation (1) as a semi-analytical solution approximation of the large-amplitude dynamics alternative orbit;
(e) like the case of the small-amplitude dynamics alternative orbit, a semi-analytical solution of the large-amplitude dynamics alternative orbit vicinity motion is constructed by analyzing the motion in its vicinity.
3. The method for designing orbits near the Earth-moon trigonometric panning points as claimed in claim 2, wherein in step (11), the gravity of mass points of the Earth, the moon and the sun is taken into account, and the relative motion is given by JPL DE 405 numerical ephemeris; under the earth center celestial coordinate system, the motion equation of the small celestial body is
Wherein GmE,GmM,GmSGravitational constants of the earth, moon and sun, R, RM、RSRespectively are the position vectors of the small celestial body, the moon and the sun under the earth center celestial coordinate system.
4. The method for designing an orbit around the earth-moon trigonometric panning point of claim 1, wherein in the step (2), the step of designing the orbit of the task which is maintained for a certain time around the earth-moon trigonometric panning points L4 and L5 under the ephemeris model specifically comprises the following steps:
(21) aiming at the small amplitude dynamics replacement orbit, the specific processing comprises the following steps:
(a) selecting a proper time interval to generate a series of track nodes according to the constructed semi-analytical solution, wherein the total time interval of the whole track nodes needs to meet the task time requirement;
(b) combining possible task boundary constraints, adopting a multipoint target shooting method, and performing numerical improvement under an ephemeris model to generate a more accurate target orbit;
(22) and generating an initial orbit value and performing numerical integration on the initial epoch selected according to the constructed semi-analytical solution aiming at the large-amplitude dynamics alternative orbit so as to design a stable dynamics alternative orbit or an orbit nearby the stable dynamics alternative orbit.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110973579.XA CN113627029B (en) | 2021-08-24 | 2021-08-24 | Method for designing track near earth-moon triangle translational point |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110973579.XA CN113627029B (en) | 2021-08-24 | 2021-08-24 | Method for designing track near earth-moon triangle translational point |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113627029A true CN113627029A (en) | 2021-11-09 |
CN113627029B CN113627029B (en) | 2022-07-12 |
Family
ID=78387442
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110973579.XA Active CN113627029B (en) | 2021-08-24 | 2021-08-24 | Method for designing track near earth-moon triangle translational point |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113627029B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4616227A (en) * | 1983-07-15 | 1986-10-07 | Hitachi, Ltd. | Method of reconstructing synthetic aperture radar image |
CN105912819A (en) * | 2016-05-06 | 2016-08-31 | 北京理工大学 | Quick design method of earth-moon L1 Lagrange point transfer orbit |
CN106672266A (en) * | 2017-01-04 | 2017-05-17 | 北京理工大学 | Balance point Halo orbit phasing orbit transfer method taking time constraint into consideration |
CN106679674A (en) * | 2016-12-05 | 2017-05-17 | 北京理工大学 | Ephemeris model-based method for analyzing shade of earth-moon L2 point Halo orbit |
CN112257343A (en) * | 2020-10-22 | 2021-01-22 | 上海卫星工程研究所 | High-precision ground track repetitive track optimization method and system |
-
2021
- 2021-08-24 CN CN202110973579.XA patent/CN113627029B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4616227A (en) * | 1983-07-15 | 1986-10-07 | Hitachi, Ltd. | Method of reconstructing synthetic aperture radar image |
CN105912819A (en) * | 2016-05-06 | 2016-08-31 | 北京理工大学 | Quick design method of earth-moon L1 Lagrange point transfer orbit |
CN106679674A (en) * | 2016-12-05 | 2017-05-17 | 北京理工大学 | Ephemeris model-based method for analyzing shade of earth-moon L2 point Halo orbit |
CN106672266A (en) * | 2017-01-04 | 2017-05-17 | 北京理工大学 | Balance point Halo orbit phasing orbit transfer method taking time constraint into consideration |
CN112257343A (en) * | 2020-10-22 | 2021-01-22 | 上海卫星工程研究所 | High-precision ground track repetitive track optimization method and system |
Non-Patent Citations (8)
Title |
---|
YING-JING QIAN,ET AL.: "Periodic motion analysis around the libration points by polynomial expansion method in planar circular restricted three-body problem", 《NONLINEAR DYNAMICS》 * |
YING-JING QIAN,ET AL.: "Periodic motion analysis around the libration points by polynomial expansion method in planar circular restricted three-body problem", 《NONLINEAR DYNAMICS》, 5 December 2017 (2017-12-05), pages 39 - 54, XP036388755, DOI: 10.1007/s11071-017-3818-1 * |
喻圣贤: "深空探测中的轨道分析、设计与控制", 《中国优秀博硕士学位论文全文数据库(博士) 工程科技II辑》 * |
喻圣贤: "深空探测中的轨道分析、设计与控制", 《中国优秀博硕士学位论文全文数据库(博士) 工程科技II辑》, no. 8, 15 August 2013 (2013-08-15), pages 88 - 91 * |
张磊: "地月系平动点导航卫星星座设计与导航性能分析", 《中国优秀博硕士学位论文全文数据库(博士) 信息科技辑》 * |
张磊: "地月系平动点导航卫星星座设计与导航性能分析", 《中国优秀博硕士学位论文全文数据库(博士) 信息科技辑》, no. 4, 15 April 2017 (2017-04-15), pages 11 - 12 * |
高珊 等: "地月拉格朗日L2点中继星轨道分析与设计", 《深空探测学报》 * |
高珊 等: "地月拉格朗日L2点中继星轨道分析与设计", 《深空探测学报》, vol. 4, no. 2, 30 April 2017 (2017-04-30), pages 122 - 129 * |
Also Published As
Publication number | Publication date |
---|---|
CN113627029B (en) | 2022-07-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Pitjeva et al. | Proposals for the masses of the three largest asteroids, the Moon-Earth mass ratio and the Astronomical Unit | |
Watkins et al. | JPL Level-2 processing standards document for Level-2 product release 05 | |
Yuan | JPL level-2 processing standards document for level-2 product release 06 | |
Visser et al. | Energy integral method for gravity field determination from satellite orbit coordinates | |
Schwarz et al. | The use of FFT techniques in physical geodesy | |
CN110132286B (en) | X-ray pulsar navigation method considering spacecraft dynamic effect and system deviation | |
Abd Rabbou et al. | Integration of GPS precise point positioning and MEMS-based INS using unscented particle filter | |
Aleksandrova et al. | A new approach to calculation of secular frequencies in the dynamics of near-Earth objects in orbits with large eccentricities | |
Han et al. | The space–time references of BeiDou navigation satellite system | |
Ely | Mean element propagations using numerical averaging | |
CN113627029B (en) | Method for designing track near earth-moon triangle translational point | |
Pou et al. | High precision SEIS calibration for the InSight mission and its applications | |
Abd-Allah | Assessment of GOCE models in Egypt | |
Zagidullin et al. | Physical libration of the Moon: An extended problem | |
Pucacco et al. | Tidal effects on the LAGEOS–LARES satellites and the LARASE program | |
Xu | Sciences of Geodesy-II: innovations and future developments | |
Folta et al. | Astrodynamics convention and modeling reference for lunar, cislunar, and libration point orbits | |
Vallado et al. | Accurate orbit determination from short-arc dense observational data | |
Christensen-Dalsgaard | The Sun as a fundamental calibrator of stellar evolution | |
Vailado et al. | Accurate orbit determination from short-arc dense observational data | |
Bettadpur | Gravity recovery and climate experiment | |
Schaffrin et al. | A Bayes filter in Friedland form for INS/GPS vector gravimetry | |
Vijayan et al. | Nonlinear dynamic modeling of satellite relative motion with differential $\pmb {J} _ {2} $ and drag | |
Zhu et al. | A novel relative navigation method based on thrust on-line identification for tight formation keeping | |
Pitjeva | EPM—High-precision planetary ephemerides of IAA RAS for scientific research and astronavigation on the Earth and in space |
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 |