CN113602531A - Method and system for generating GEO orbital transfer strategy of assembly under abnormal separation condition - Google Patents

Method and system for generating GEO orbital transfer strategy of assembly under abnormal separation condition Download PDF

Info

Publication number
CN113602531A
CN113602531A CN202110758134.XA CN202110758134A CN113602531A CN 113602531 A CN113602531 A CN 113602531A CN 202110758134 A CN202110758134 A CN 202110758134A CN 113602531 A CN113602531 A CN 113602531A
Authority
CN
China
Prior art keywords
transfer
orbit
orbital
abnormal
target
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
CN202110758134.XA
Other languages
Chinese (zh)
Other versions
CN113602531B (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 Institute of Satellite Engineering
Original Assignee
Shanghai Institute of Satellite Engineering
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 Institute of Satellite Engineering filed Critical Shanghai Institute of Satellite Engineering
Priority to CN202110758134.XA priority Critical patent/CN113602531B/en
Publication of CN113602531A publication Critical patent/CN113602531A/en
Application granted granted Critical
Publication of CN113602531B publication Critical patent/CN113602531B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B64—AIRCRAFT; AVIATION; COSMONAUTICS
    • B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00—Cosmonautic vehicles
    • B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/242—Orbits and trajectories
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B64—AIRCRAFT; AVIATION; COSMONAUTICS
    • B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00—Cosmonautic vehicles
    • B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/242—Orbits and trajectories
    • B64G1/2427—Transfer orbits
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B64—AIRCRAFT; AVIATION; COSMONAUTICS
    • B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00—Cosmonautic vehicles
    • B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/40—Arrangements or adaptations of propulsion systems
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Navigation (AREA)

Abstract

The invention provides a method and a system for generating a GEO orbital transfer strategy of an assembly under the condition of abnormal separation, which comprises the following steps: determining an orbit parameter of a separation point of the assembly and a platform parameter related to the generation of an orbit transfer strategy of the assembly according to the state of the abnormal separation moment between the assembly and the carrier or between the assemblies; determining a design constraint condition of a track transfer strategy according to the platform parameters; according to the constraint conditions, carrying out optimal orbital transfer times estimation on the basis of the principle that the speed increment required by orbital transfer is minimum; establishing a mathematical model for describing the single orbital transfer state change and the connection among multiple orbital transfers according to the estimation result; and generating an orbit transfer strategy meeting design constraints by using the mathematical model and taking the GEO orbit as a target through iterative optimization. The method overcomes the defects of large calculated amount, low calculating speed and the like of the traditional limited thrust complex design method under the condition of abnormal separation of the combined spacecraft, and has certain engineering practicability.

Description

Method and system for generating GEO orbital transfer strategy of assembly under abnormal separation condition
Technical Field
The invention relates to astronavigation aircraft orbital dynamics, in particular to a method and a system for generating a GEO orbital transfer strategy of an assembly under the condition of abnormal separation.
Background
Geostationary orbit spacecraft (GEO) has the characteristic of being relatively static to the earth, and is widely applied to the fields of communication, navigation, relay and the like. Compared with the traditional single spacecraft, the combined spacecraft has the advantages that each cabin section carries the propulsion system, so that on one hand, the combined spacecraft has better robustness when dealing with the abnormal conditions of the active section and the multiple orbit transfer sections, and on the other hand, the orbit transfer strategy needs to be designed for the combined body or each cabin section of the combined body rapidly under the condition that the combined body is separated abnormally.
At present, a plurality of transfer orbit strategy design schemes are applied at home and abroad, but the transfer orbit strategy design schemes are mostly complex integral optimization based on limited thrust or a segmented orbit transfer method based on a hybrid propulsion system. In the Chinese patent of "a method for transferring and converting pulse orbital transfer to limited thrust orbital transfer" (patent document: CN112455725A), such as Liujun Yao, Zhao Jian Wei, Zongshi, etc., firstly, the position of a pulse ignition point is determined, then limited thrust integral optimization is carried out at the position, and the ignition direction is corrected. In the Chinese patent 'small geostationary satellite orbit transfer method and system' (patent document: CN111891396A), Linbao army, Jiang Guo Wei, Fanyuan and the like, a method is provided, wherein a chemical propulsion system is firstly utilized to lift the altitude of a spacecraft at the near site, an electric propulsion system is utilized to adjust the altitude, the inclination angle and the eccentricity ratio of the spacecraft at the near site, and finally the chemical propulsion system is utilized to capture the spacecraft at a fixed point; on one hand, the altitude, the inclination angle and the eccentricity of the near place are adjusted in stages, the fuel consumption is large, and the on-orbit service life of the spacecraft is influenced.
In summary, the design optimization of the orbital transfer strategy method needs to be developed according to the requirement for rapidly generating the GEO orbit transfer strategy of the assembly spacecraft under the abnormal separation condition.
Disclosure of Invention
Aiming at the defects in the prior art, the invention aims to provide a method and a system for generating a GEO orbital transfer strategy of an assembly under the condition of abnormal separation.
The invention provides a method for generating a GEO orbital transfer strategy of an assembly under the condition of abnormal separation, which comprises the following steps:
step A: determining an orbit parameter of a separation point of the assembly and a platform parameter related to the generation of an orbit transfer strategy of the assembly according to the state of the abnormal separation moment between the assembly and the carrier or between the assemblies;
and B: determining a design constraint condition of a track transfer strategy according to the platform parameters;
and C: according to the constraint conditions, carrying out optimal orbital transfer times estimation on the basis of the principle that the speed increment required by orbital transfer is minimum;
step D: establishing a mathematical model for describing the single orbital transfer state change and the connection among multiple orbital transfers according to the estimation result;
step E: and generating an orbit transfer strategy meeting design constraints by using the mathematical model and taking the GEO orbit as a target through iterative optimization.
Preferably, the step a includes:
step S2.1: the combination body comprises a plurality of cabin sections, each cabin section carries a propulsion system, and the abnormal separation comprises a launching section abnormal and a transfer section abnormal, wherein the launching section abnormal and the transfer section abnormal are separated in advance, the combination body cannot be sent to a preset orbit by a carrier, and the transfer section abnormal are separated in advance; firstly, the orbit state at the moment of abnormal separation is determined, the general form of GEO orbit is adopted for description, and the semimajor axis a0Angle of inclination i0Eccentricity e0And the groundLongitude λ0Geographic latitude η0；
Step S2.2: defining the relevant platform parameters of the combined spacecraft, main satellite or main satellite and propulsion pod, including the separation moment weight m0M 'of available Fuel residual quantity'0Thrust F of engine0Specific impulse Isp of engine0。
Preferably, in the step B:
three design constraints are summarized by combining actual engineering: and the ground measurement and control conditions restrict the geographical longitude span and the measurement and control duration of the ignition section, the characteristics of the thruster, the loss of the arc section and the safety protection requirements restrict the single longest ignition time and restrict the available fuel of the rail transfer object.
Preferably, in the step C, in order to determine the dimension of the track transfer optimization parameter, the optimal number of times of track transfer needs to be estimated according to the following steps on the principle that the increment of the speed required by the transfer is minimum:
step S4.1: and calculating the minimum velocity increment delta v of the spacecraft from the abnormal separation moment state to the GEO target orbit according to the following formula:
Figure BDA0003147993000000031
Δv＝vn-v0
wherein ,r0Designing the earth center distance of the abnormal separation moment of the object for the transfer orbit, and calculating the earth center distance according to the number of the orbits; a is0Is a semi-major axis of the track; mu is an earth gravity constant; v. of0Designing the speed of the abnormal separation moment of the object for the transfer orbit; v. ofnIs the target track speed; a isnIs a target track semi-major axis;
step S4.2: calculating fuel consumption delta m corresponding to the minimum velocity increment of the orbital transfer according to a rocket formula:
Figure BDA0003147993000000032
wherein g is 9.80665 m-s2Is the acceleration of the earth's gravity;
step S4.3: calculating the second flow dm of the propellant according to the specific impulse of the engine, and calculating the total track change time t by combining the fuel consumption:
Figure BDA0003147993000000033
Figure BDA0003147993000000034
step S4.4: single maximum ignition time T combined with thrustermaxCalculating to obtain an optimal orbital transfer frequency estimated value N:
Figure BDA0003147993000000035
wherein [ x ] is an eave function and represents the minimum integer which is more than or equal to x.
Preferably, the step D includes:
step S5.1: according to the variable to be optimized for each ignition: semi-major axis a of orbital transfer targetkBefore track change, offset circle QkK is 1, …, N, calculating pre-ignition parameters:
track period Tk：
Figure BDA0003147993000000036
Drift rate of longitude
Figure BDA0003147993000000037
Figure BDA0003147993000000038
wherein ,ωe＝7.2921×10-5rad/s, and pi is the circumference ratio.
Point of intersection of the riseGeographic longitude λk：
Figure BDA0003147993000000041
Step S5.2: according to the semi-major axis a of the target at each orbital transferkAnd k is 1, …, N, and the single ignition point and the target track parameter are calculated as follows:
calculating the target track velocity v by the activity formulak：
Figure BDA0003147993000000042
wherein ,rkDesigning the ground center distance of the object at the ignition moment for the transfer orbit, and calculating the ground center distance according to the number of the orbits;
orbital transfer velocity increment Δ vk：
Figure BDA0003147993000000043
Wherein, α and βkCalculated according to the following formula
Figure BDA0003147993000000044
βk＝π-ik-1-α
Track-changing back inclination angle ik：
Figure BDA0003147993000000045
Orbital transfer fuel consumption Δ mk：
Figure BDA0003147993000000046
Ignition duration tk：
Figure BDA0003147993000000047
Preferably, the step E includes:
step S6.1: converting the GEO orbit transfer strategy solving problem into a multivariable, multi-target and multi-constraint optimizing problem, wherein an optimization model is described as follows:
an objective function:
Figure BDA0003147993000000048
wherein ,
Figure BDA0003147993000000051
representing a track target penalty function;
Figure BDA0003147993000000052
representing the deviation between the jth orbital parameter orbital transfer final value and a target orbit; Δ Orb ═ (Δ a, Δ i, Δ e, Δ λ, Δ η), and represents the deviation of the tracking result from the target track; gamma represents a weight coefficient;
constraint conditions are as follows:
Figure BDA0003147993000000053
step S6.2: selecting a multivariable, multi-target and multi-constraint optimization algorithm for iterative optimization aiming at the optimization model;
step S6.3: and outputting an optimization result, and determining a GEO orbit transfer strategy of the assembly spacecraft under the abnormal separation condition.
The invention provides a system for generating a GEO orbital transfer strategy of an assembly under the condition of abnormal separation, which comprises:
a module A: determining an orbit parameter of a separation point of the assembly and a platform parameter related to the generation of an orbit transfer strategy of the assembly according to the state of the abnormal separation moment between the assembly and the carrier or between the assemblies;
and a module B: determining a design constraint condition of a track transfer strategy according to the platform parameters;
and a module C: according to the constraint conditions, carrying out optimal orbital transfer times estimation on the basis of the principle that the speed increment required by orbital transfer is minimum;
a module D: establishing a mathematical model for describing the single orbital transfer state change and the connection among multiple orbital transfers according to the estimation result;
and a module E: and generating an orbit transfer strategy meeting design constraints by using the mathematical model and taking the GEO orbit as a target through iterative optimization.
Preferably, the module a comprises:
module S2.1: the combination body comprises a plurality of cabin sections, each cabin section carries a propulsion system, and the abnormal separation comprises a launching section abnormal and a transfer section abnormal, wherein the launching section abnormal and the transfer section abnormal are separated in advance, the combination body cannot be sent to a preset orbit by a carrier, and the transfer section abnormal are separated in advance; firstly, the orbit state at the moment of abnormal separation is determined, the general form of GEO orbit is adopted for description, and the semimajor axis a0Angle of inclination i0Eccentricity e0Geographic longitude λ0Geographic latitude η0；
Module S2.2: defining the relevant platform parameters of the combined spacecraft, main satellite or main satellite and propulsion pod, including the separation moment weight m0M 'of available Fuel residual quantity'0Thrust F of engine0Specific impulse Isp of engine0。
In the module B:
three design constraints are summarized by combining actual engineering: and the ground measurement and control conditions restrict the geographical longitude span and the measurement and control duration of the ignition section, the characteristics of the thruster, the loss of the arc section and the safety protection requirements restrict the single longest ignition time and restrict the available fuel of the rail transfer object.
Preferably, in the module C, in order to determine the optimal dimension of the track transfer parameter, the optimal number of times of track transfer needs to be estimated according to the following module based on the principle that the increment of the speed required by the transfer is minimum:
module S4.1: and calculating the minimum velocity increment delta v of the spacecraft from the abnormal separation moment state to the GEO target orbit according to the following formula:
Figure BDA0003147993000000061
Δv＝vn-v0
wherein ,r0Designing the earth center distance of the abnormal separation moment of the object for the transfer orbit, and calculating the earth center distance according to the number of the orbits; a is0Is a semi-major axis of the track; mu is an earth gravity constant; v. of0Designing the speed of the abnormal separation moment of the object for the transfer orbit; v. ofnIs the target track speed; a isnIs a target track semi-major axis;
module S4.2: calculating fuel consumption delta m corresponding to the minimum velocity increment of the orbital transfer according to a rocket formula:
Figure BDA0003147993000000062
wherein g is 9.80665m/s2Is the acceleration of the earth's gravity;
module S4.3: calculating the second flow dm of the propellant according to the specific impulse of the engine, and calculating the total track change time t by combining the fuel consumption:
Figure BDA0003147993000000063
Figure BDA0003147993000000064
module S4.4: single maximum ignition time T combined with thrustermaxCalculating to obtain an optimal orbital transfer frequency estimated value N:
Figure BDA0003147993000000065
wherein [ x ] is an eave function and represents the minimum integer which is more than or equal to x.
Preferably, the module D comprises:
module S5.1: according to the variable to be optimized for each ignition: semi-major axis a of orbital transfer targetkBefore track change, offset circle QkK is 1, …, N, calculating pre-ignition parameters:
track period Tk：
Figure BDA0003147993000000071
Drift rate of longitude
Figure BDA0003147993000000072
Figure BDA0003147993000000073
wherein ,ωe＝7.2921×10-5rad/s, and pi is the circumference ratio.
Ascending node geographic longitude λk：
Figure BDA0003147993000000074
Module S5.2: according to the semi-major axis a of the target at each orbital transferkAnd k is 1, …, N, and the single ignition point and the target track parameter are calculated as follows:
calculating the target track velocity v by the activity formulak：
Figure BDA0003147993000000075
wherein ,rkDesigning the ground center distance of the object at the ignition moment for the transfer orbit, and calculating the ground center distance according to the number of the orbits;
orbital transfer velocity increment Δ vk：
Figure BDA0003147993000000076
Wherein, α and βkCalculated according to the following formula
Figure BDA0003147993000000077
βk＝π-ik-1-α
Track-changing back inclination angle ik：
Figure BDA0003147993000000078
Orbital transfer fuel consumption Δ mk：
Figure BDA0003147993000000079
Ignition duration tk：
Figure BDA00031479930000000710
Preferably, said module E comprises:
module S6.1: converting the GEO orbit transfer strategy solving problem into a multivariable, multi-target and multi-constraint optimizing problem, wherein an optimization model is described as follows:
an objective function:
Figure BDA0003147993000000081
wherein ,
Figure BDA0003147993000000082
representing a track target penalty function;
Figure BDA0003147993000000083
denotes the jthDeviation between the track parameter orbital transfer final value and a target track; Δ Orb ═ (Δ a, Δ i, Δ e, Δ λ, Δ η), and represents the deviation of the tracking result from the target track; gamma represents a weight coefficient;
constraint conditions are as follows:
Figure BDA0003147993000000084
module S6.2: selecting a multivariable, multi-target and multi-constraint optimization algorithm for iterative optimization aiming at the optimization model;
module S6.3: and outputting an optimization result, and determining a GEO orbit transfer strategy of the assembly spacecraft under the abnormal separation condition.
Compared with the prior art, the invention has the following beneficial effects:
the method overcomes the defects of large calculated amount, low calculating speed and the like of the traditional limited thrust complex design method under the condition of abnormal separation of the combined spacecraft, and has certain engineering practicability.
Drawings
Other features, objects and advantages of the invention will become more apparent upon reading of the detailed description of non-limiting embodiments with reference to the following drawings:
FIG. 1 is a schematic block diagram of a method for rapidly generating an assembly GEO orbital transfer strategy under an abnormal separation condition;
FIG. 2 is a schematic view of a configuration of an assembled spacecraft (two sections);
FIG. 3 is a schematic view of an abnormally separated (out of tolerance) orbit of an assembled spacecraft from a carrier;
FIG. 4 is a diagram of a GEO orbit transfer strategy simulation of an assembly spacecraft based on abnormal separation conditions.
Detailed Description
The present invention will be described in detail with reference to specific examples. The following examples will assist those skilled in the art in further understanding the invention, but are not intended to limit the invention in any way. It should be noted that it would be obvious to those skilled in the art that various changes and modifications can be made without departing from the spirit of the invention. All falling within the scope of the present invention.
As shown in the attached figure 1, the invention provides a method for quickly generating a GEO (geosynchronous orbit) orbit-changing strategy of an assembly under the condition of abnormal separation. The method specifically comprises the following steps:
step A: and determining the track parameters of the separation point of the assembly and the platform parameters related to the generation of the track transfer strategy of the assembly according to the abnormal separation time state between the assembly and the carrier or between the assemblies.
And B: and determining the design constraint conditions of the track transfer strategy according to factors such as ground measurement and control requirements, single ignition capability and the like.
And C: and (4) carrying out optimal orbital transfer times estimation on the basis of the minimum speed increment required by orbital transfer.
Step D: as shown in fig. 4, a mathematical model describing the relationship between the state change of single orbital transfer and the multiple orbital transfers is established.
Step E: and with the GEO orbit as a target, quickly generating an orbit transfer strategy meeting design constraints through iterative optimization.
The step A comprises the following steps:
step S2.1: as shown in fig. 2, a combed spacecraft typically includes a plurality of sections, and each section carries a propulsion system. As shown in fig. 3, the abnormal separation includes an abnormal launching section in which the vehicle fails to send the assembled spacecraft to a predetermined orbit, i.e., the separation in advance, and an abnormal transfer section in which a cabin section in the assembled spacecraft fails to send the assembled spacecraft to a quasi-geosynchronous orbit, i.e., the separation in advance. Without loss of generality, it is described in the general form of a GEO orbit, i.e. semimajor axis a0Angle of inclination i0Eccentricity e0Geographic longitude λ0Geographic latitude η0；
Step S2.2: in addition, there is a need to specify relevant platform parameters for the transfer strategy design object (combined spacecraft, main satellite or main satellite and propulsion pod), typically including the separation time moment weight m0M 'of available Fuel residual quantity'0Thrust F of engine0Specific impulse Isp of engine0；
And step B, summarizing design constraints in 3 aspects by combining actual engineering requirements: the ground measurement and control conditions are used for constraining the geographical longitude span and the measurement and control duration of the ignition section, the characteristics of a thruster, the loss of an arc section and the safety protection requirements on the longest ignition time of a single ignition and the available fuel of a rail transfer object;
in step C, in order to determine the optimal dimension of the track transfer parameter, the optimal number of times of track transfer needs to be estimated according to the following steps and the principle that the increment of the speed required by transfer is minimum:
step S4.1: and calculating the minimum velocity increment delta v of the spacecraft from the abnormal separation moment state to the GEO target orbit according to the following formula:
Figure BDA0003147993000000101
Δv＝vn-v0
wherein ,r0Designing the ground center distance of the abnormal separation moment of the object for the transfer orbit, and calculating the ground center distance according to the number of the orbits; v. of0Designing the speed of the abnormal separation moment of the object for the transfer orbit; v. ofnIs the target track speed; a isnThe GEO orbit is 42164km which is the semimajor axis of the target orbit;
step S4.2: calculating fuel consumption delta m corresponding to the minimum velocity increment of the orbital transfer according to a rocket formula:
Figure BDA0003147993000000102
wherein g is 9.80665m/s2Is the acceleration of the earth's gravity;
step S4.3: calculating the second flow dm of the propellant according to the specific impulse of the engine, and calculating the total track change time t by combining the fuel consumption:
Figure BDA0003147993000000103
Figure BDA0003147993000000104
step S4.4: single maximum ignition time T combined with thrustermaxCalculating to obtain an optimal orbital transfer frequency estimated value N:
Figure BDA0003147993000000105
wherein [ x ] is an eave function and represents the minimum integer which is more than or equal to x.
The step D comprises the following steps:
step S5.1: according to the variable to be optimized for each ignition, i.e. the semi-major axis a of the target of the orbital transferkBefore track change, offset circle QkK — 1, …, N calculates the pre-ignition parameters:
track period Tk：
Figure BDA0003147993000000106
Drift rate of longitude
Figure BDA0003147993000000107
Figure BDA0003147993000000111
wherein ,ωe＝7.2921×10-5rad/s
Ascending node geographic longitude λk：
Figure BDA0003147993000000112
Step S5.2: according to the semi-major axis a of the target at each orbital transferkAnd k is 1, …, N, and the single ignition point and the target track parameter are calculated as follows:
calculating the target orbit according to the activity formulaVelocity vk：
Figure BDA0003147993000000113
wherein ,rkDesigning the ground center distance of the object at the ignition moment for the transfer orbit, and calculating the ground center distance according to the number of the orbits;
orbital transfer velocity increment Δ vk：
Figure BDA0003147993000000114
Wherein, α and βkCalculated according to the following formula
Figure BDA0003147993000000115
βk＝π-ik-1-α
Track-changing back inclination angle ik：
Figure BDA0003147993000000116
Orbital transfer fuel consumption Δ mk：
Figure BDA0003147993000000117
Ignition duration tk：
Figure BDA0003147993000000118
The step E comprises the following steps:
step S6.1: converting the GEO orbit transfer strategy solving problem into a multivariable, multi-target and multi-constraint optimizing problem, wherein an optimization model is described as follows:
an objective function:
Figure BDA0003147993000000121
wherein ,
Figure BDA0003147993000000122
representing a track target penalty function;
Figure BDA0003147993000000123
representing the deviation between the jth orbital parameter orbital transfer final value and a target orbit; Δ Orb ═ (Δ a, Δ i, Δ e, Δ λ, Δ η), and represents the deviation of the tracking result from the target track; gamma represents a weight coefficient, and a larger positive integer is taken for carrying out large-weight punishment on the parameters which do not meet the precision index;
constraint conditions are as follows:
Figure BDA0003147993000000124
step S6.2: selecting the existing multivariate, multi-target and multi-constraint optimization algorithm for iterative optimization aiming at the optimization model;
step S6.3: and outputting an optimization result, and determining a GEO orbit transfer strategy of the assembly spacecraft under the abnormal separation condition.
In this embodiment, assuming that an anomaly occurs at the time of separation of the vehicle and the combined spacecraft, the orbit inclination angle is deviated at the time of separation due to the fault of the vehicle, and the separation orbit parameter is determined according to step a:
Orb0＝(a0,i0,e0,λ0,η0)＝(24473.64km,30°,0.731215,-9.21°,30°)
determining platform parameters related to the combination body and the track transfer as follows: weight m at the moment of separation03863kg of available residual fuel quantity m'01820kg, engine thrust F0350N, engine specific impulse Isp0＝315s；
Determining track transfer strategy design constraints according to step B: shortest measurement and control duration tfireMore than or equal to 30min, ground measurement and control geographic longitude range [30 DEG E, 170 DEG E]Maximum time t of single ignition of engineengine4050s or less, and most fuel is available for rail transfer
Figure BDA0003147993000000125
According to step C, the GEO target orbit parameter Orb is combinednThe number of times of tracking was estimated (42164km,30 °,0,60 °,0 °):
Figure BDA0003147993000000126
and finally, according to the optimization model established in the step E, combining the step D to update parameters of single orbital transfer state change and multiple orbital transfer, and generating a track transfer strategy which meets design constraints through iterative optimization, wherein the track transfer strategy is shown in the table 1.
TABLE 1 iterative optimization solution for orbital transfer strategy
Figure BDA0003147993000000131
Figure BDA0003147993000000141
The invention also provides a system for generating the GEO orbital transfer strategy of the combination under the condition of abnormal separation, which comprises the following steps:
a module A: and determining the track parameters of the separation point of the assembly and the platform parameters related to the generation of the track transfer strategy of the assembly according to the abnormal separation time state between the assembly and the carrier or between the assemblies.
And a module B: and determining the design constraint conditions of the track transfer strategy according to the platform parameters.
And a module C: and according to the constraint conditions, carrying out optimal orbital transfer times estimation on the basis of the principle that the speed increment required by orbital transfer is minimum.
A module D: and establishing a mathematical model for describing the single orbital transfer state change and the connection among multiple orbital transfers according to the estimation result.
And a module E: and generating an orbit transfer strategy meeting design constraints by using the mathematical model and taking the GEO orbit as a target through iterative optimization.
Those skilled in the art will appreciate that, in addition to implementing the system and its various devices, modules, units provided by the present invention as pure computer readable program code, the system and its various devices, modules, units provided by the present invention can be fully implemented by logically programming method steps in the form of logic gates, switches, application specific integrated circuits, programmable logic controllers, embedded microcontrollers and the like. Therefore, the system and various devices, modules and units thereof provided by the invention can be regarded as a hardware component, and the devices, modules and units included in the system for realizing various functions can also be regarded as structures in the hardware component; means, modules, units for performing the various functions may also be regarded as structures within both software modules and hardware components for performing the method.
The foregoing description of specific embodiments of the present invention has been presented. It is to be understood that the present invention is not limited to the specific embodiments described above, and that various changes or modifications may be made by one skilled in the art within the scope of the appended claims without departing from the spirit of the invention. The embodiments and features of the embodiments of the present application may be combined with each other arbitrarily without conflict.

Claims (10)

1. A method for generating a GEO orbital transfer strategy of an assembly under the condition of abnormal separation is characterized by comprising the following steps:
step A: determining an orbit parameter of a separation point of the assembly and a platform parameter related to the generation of an orbit transfer strategy of the assembly according to the state of the abnormal separation moment between the assembly and the carrier or between the assemblies;
and B: determining a design constraint condition of a track transfer strategy according to the platform parameters;
and C: according to the constraint conditions, carrying out optimal orbital transfer times estimation on the basis of the principle that the speed increment required by orbital transfer is minimum;
step D: establishing a mathematical model for describing the single orbital transfer state change and the connection among multiple orbital transfers according to the estimation result;
step E: and generating an orbit transfer strategy meeting design constraints by using the mathematical model and taking the GEO orbit as a target through iterative optimization.
2. The method for generating the GEO orbital transfer strategy of the combined body under the abnormal separation condition according to claim 1, wherein the step a comprises:
step S2.1: the combination body comprises a plurality of cabin sections, each cabin section carries a propulsion system, and the abnormal separation comprises a launching section abnormal and a transfer section abnormal, wherein the launching section abnormal and the transfer section abnormal are separated in advance, the combination body cannot be sent to a preset orbit by a carrier, and the transfer section abnormal are separated in advance; firstly, the orbit state at the moment of abnormal separation is determined, the general form of GEO orbit is adopted for description, and the semimajor axis a0Angle of inclination i0Eccentricity e0Geographic longitude λ0Geographic latitude η0；
Step S2.2: defining the relevant platform parameters of the combined spacecraft, main satellite or main satellite and propulsion pod, including the separation moment weight m0M 'of available Fuel residual quantity'0Thrust F of engine0Specific impulse Isp of engine0。
3. The method for generating the GEO-orbital transfer strategy of the assembly under the abnormal separation condition according to claim 2, wherein in the step C, in order to determine the optimal orbital transfer parameter dimension, the optimal orbital transfer times are estimated according to the following steps based on the principle that the increment of the speed required by the transfer is minimum:
step S4.1: and calculating the minimum velocity increment delta v of the spacecraft from the abnormal separation moment state to the GEO target orbit according to the following formula:
Figure FDA0003147992990000011
Δv＝vn-v0
wherein ,r0Designing the earth center distance of the abnormal separation moment of the object for the transfer orbit, and calculating the earth center distance according to the number of the orbits; a is0Is a semi-major axis of the track; mu is an earth gravity constant; v. of0Designing the speed of the abnormal separation moment of the object for the transfer orbit; v. ofnIs the target track speed; a isnIs a target track semi-major axis;
step S4.2: calculating fuel consumption delta m corresponding to the minimum velocity increment of the orbital transfer according to a rocket formula:
Figure FDA0003147992990000021
wherein g is 9.80665m/s2Is the acceleration of the earth's gravity;
step S4.3: calculating the second flow dm of the propellant according to the specific impulse of the engine, and calculating the total track change time t by combining the fuel consumption:
Figure FDA0003147992990000022
Figure FDA0003147992990000023
step S4.4: single maximum ignition time T combined with thrustermaxCalculating to obtain an optimal orbital transfer frequency estimated value N:
Figure FDA0003147992990000024
wherein [ x ] is an eave function and represents the minimum integer which is more than or equal to x.
4. The method for generating the GEO orbital transfer strategy of the assembly under the abnormal separation condition according to claim 3, wherein the step D comprises:
step S5.1: according to the variable to be optimized for each ignition: semi-major axis a of orbital transfer targetkBefore track change, offset circle QkK is 1, …, N, calculating pre-ignition parameters:
track period Tk：
Figure FDA0003147992990000025
Drift rate of longitude
Figure FDA0003147992990000026
Figure FDA0003147992990000027
wherein ,ωe＝7.2921×10-5rad/s, and pi is the circumference ratio.
Ascending node geographic longitude λk：
Figure FDA0003147992990000028
Step S5.2: according to the semi-major axis a of the target at each orbital transferkAnd k is 1, …, N, and the single ignition point and the target track parameter are calculated as follows:
calculating the target track velocity v by the activity formulak：
Figure FDA0003147992990000031
wherein ,rkDesigning the ground center distance of the object at the ignition moment for the transfer orbit, and calculating the ground center distance according to the number of the orbits;
orbital transfer velocity increment Δ vk：
Figure FDA0003147992990000032
Wherein, α and βkCalculated according to the following formula
Figure FDA0003147992990000033
βk＝π-ik-1-α
Track-changing back inclination angle ik：
Figure FDA0003147992990000034
Orbital transfer fuel consumption Δ mk：
Figure FDA0003147992990000035
Ignition duration tk：
Figure FDA0003147992990000036
5. The method for generating the GEO orbital transfer strategy of the combined body under the abnormal separation condition according to claim 4, wherein the step E comprises:
step S6.1: converting the GEO orbit transfer strategy solving problem into a multivariable, multi-target and multi-constraint optimizing problem, wherein an optimization model is described as follows:
an objective function:
Figure FDA0003147992990000037
wherein ,
Figure FDA0003147992990000038
representing a track target penalty function;
Figure FDA0003147992990000039
representing the deviation between the jth orbital parameter orbital transfer final value and a target orbit; Δ Orb ═ (Δ a, Δ i, Δ e, Δ λ, Δ η), and represents the deviation of the tracking result from the target track; gamma represents a weight coefficient;
constraint conditions are as follows:
Figure FDA0003147992990000041
step S6.2: selecting a multivariable, multi-target and multi-constraint optimization algorithm for iterative optimization aiming at the optimization model;
step S6.3: and outputting an optimization result, and determining a GEO orbit transfer strategy of the assembly spacecraft under the abnormal separation condition.
6. An assembly GEO orbital transfer strategy generation system under the condition of abnormal separation is characterized by comprising:
a module A: determining an orbit parameter of a separation point of the assembly and a platform parameter related to the generation of an orbit transfer strategy of the assembly according to the state of the abnormal separation moment between the assembly and the carrier or between the assemblies;
and a module B: determining a design constraint condition of a track transfer strategy according to the platform parameters;
and a module C: according to the constraint conditions, carrying out optimal orbital transfer times estimation on the basis of the principle that the speed increment required by orbital transfer is minimum;
a module D: establishing a mathematical model for describing the single orbital transfer state change and the connection among multiple orbital transfers according to the estimation result;
and a module E: and generating an orbit transfer strategy meeting design constraints by using the mathematical model and taking the GEO orbit as a target through iterative optimization.
7. The system for generating a GEO-orbital transfer strategy for an abnormal separation situation of claim 6, wherein the module a comprises:
module S2.1: the combination body comprises a plurality of cabin sections, each cabin section carries a propulsion system, and the abnormal separation comprises a launching section abnormal and a transfer section abnormal, wherein the launching section abnormal and the transfer section abnormal are separated in advance, the combination body cannot be sent to a preset orbit by a carrier, and the transfer section abnormal are separated in advance; firstly, the orbit state at the moment of abnormal separation is determined, the general form of GEO orbit is adopted for description, and the semimajor axis a0Angle of inclination i0Eccentricity e0Geographic longitude λ0Geographic latitude η0；
Module S2.2: defining the relevant platform parameters of the combined spacecraft, main satellite or main satellite and propulsion pod, including the separation moment weight m0M 'of available Fuel residual quantity'0Thrust F of engine0Specific impulse Isp of engine0。
8. The system for generating the GEO-orbital transfer strategy for the combination under the abnormal separation condition of claim 7, wherein in the module C, in order to determine the optimal orbital transfer parameter dimension, the optimal orbital transfer times estimation is performed according to the following modules on the principle that the increment of the speed required by the transfer is minimum:
module S4.1: and calculating the minimum velocity increment delta v of the spacecraft from the abnormal separation moment state to the GEO target orbit according to the following formula:
Figure FDA0003147992990000051
Δv＝vn-v0
wherein ,r0Designing the earth center distance of the abnormal separation moment of the object for the transfer orbit, and calculating the earth center distance according to the number of the orbits; a is0Is a semi-major axis of the track; mu is an earth gravity constant; v. of0Designing the speed of the abnormal separation moment of the object for the transfer orbit; v. ofnIs the target track speed; a isnIs a target track semi-major axis;
module S4.2: calculating fuel consumption delta m corresponding to the minimum velocity increment of the orbital transfer according to a rocket formula:
Figure FDA0003147992990000052
wherein g is 9.80665m/s2Is the acceleration of the earth's gravity;
module S4.3: calculating the second flow dm of the propellant according to the specific impulse of the engine, and calculating the total track change time t by combining the fuel consumption:
Figure FDA0003147992990000053
Figure FDA0003147992990000054
module S4.4: single maximum ignition time T combined with thrustermaxCalculating to obtain an optimal orbital transfer frequency estimated value N:
Figure FDA0003147992990000055
wherein [ x ] is an eave function and represents the minimum integer which is more than or equal to x.
9. The system for generating a GEO-orbital transfer strategy for an abnormal separation situation of claim 8, wherein the module D comprises:
module S5.1: according to the variable to be optimized for each ignition: semi-major axis a of orbital transfer targetkBefore track change, offset circle QkK is 1, …, N, calculating pre-ignition parameters:
track period Tk：
Figure FDA0003147992990000056
Drift rate of longitude
Figure FDA0003147992990000057
Figure FDA0003147992990000058
wherein ,ωe＝7.2921×10-5rad/s, and pi is the circumference ratio.
Ascending node geographic longitude λk：
Figure FDA0003147992990000061
Module S5.2: according to the semi-major axis a of the target at each orbital transferkAnd k is 1, …, N, and the single ignition point and the target track parameter are calculated as follows:
calculating the target track velocity v by the activity formulak：
Figure FDA0003147992990000062
wherein ,rkDesigning the ground center distance of the object at the ignition moment for the transfer orbit, and calculating the ground center distance according to the number of the orbits;
orbital transfer velocity increment Δ vk：
Figure FDA0003147992990000063
Wherein, α and βkCalculated according to the following formula
Figure FDA0003147992990000064
βk＝π-ik-1-α
Track-changing back inclination angle ik：
Figure FDA0003147992990000065
Orbital transfer fuel consumption Δ mk：
Figure FDA0003147992990000066
Ignition duration tk：
Figure FDA0003147992990000067
10. The system for generating a GEO-orbital transfer strategy for an abnormal separation situation of claim 10, wherein the module E comprises:
module S6.1: converting the GEO orbit transfer strategy solving problem into a multivariable, multi-target and multi-constraint optimizing problem, wherein an optimization model is described as follows:
an objective function:
Figure FDA0003147992990000071
wherein ,
Figure FDA0003147992990000072
representing a track target penalty function;
Figure FDA0003147992990000073
representing the deviation between the jth orbital parameter orbital transfer final value and a target orbit; Δ Orb ═ (Δ a, Δ i, Δ e, Δ λ, Δ η), and represents the deviation of the tracking result from the target track; gamma represents a weight coefficient;
constraint conditions are as follows:
Figure FDA0003147992990000074
module S6.2: selecting a multivariable, multi-target and multi-constraint optimization algorithm for iterative optimization aiming at the optimization model;
module S6.3: and outputting an optimization result, and determining a GEO orbit transfer strategy of the assembly spacecraft under the abnormal separation condition.
CN202110758134.XA 2021-07-05 2021-07-05 Method and system for generating combined GEO (generic object oriented) orbit strategy under abnormal separation condition Active CN113602531B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110758134.XA CN113602531B (en) 2021-07-05 2021-07-05 Method and system for generating combined GEO (generic object oriented) orbit strategy under abnormal separation condition

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110758134.XA CN113602531B (en) 2021-07-05 2021-07-05 Method and system for generating combined GEO (generic object oriented) orbit strategy under abnormal separation condition

Publications (2)

Publication Number Publication Date
CN113602531A true CN113602531A (en) 2021-11-05
CN113602531B CN113602531B (en) 2023-05-12

Family

ID=78337265

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110758134.XA Active CN113602531B (en) 2021-07-05 2021-07-05 Method and system for generating combined GEO (generic object oriented) orbit strategy under abnormal separation condition

Country Status (1)

Country Link
CN (1) CN113602531B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102424116A (en) * 2011-12-08 2012-04-25 中国空间技术研究院 Method for optimizing orbital transfer strategy of geostationary orbit satellite
EP2990338A1 (en) * 2014-08-28 2016-03-02 The Boeing Company Satellite transfer orbit search methods
CN108216687A (en) * 2017-12-25 2018-06-29 中国空间技术研究院 GEO satellite based on particle cluster algorithm becomes rail policy calculation method, system and medium
CN109625323A (en) * 2018-11-09 2019-04-16 中国科学院空间应用工程与技术中心 A kind of satellite chemical propulsion orbit changing method and system
US20200055617A1 (en) * 2018-08-17 2020-02-20 Mitsubishi Electric Research Laboratories, Inc. System and Method of Tracking a Spacecraft Trajectory for Orbital Transfer

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102424116A (en) * 2011-12-08 2012-04-25 中国空间技术研究院 Method for optimizing orbital transfer strategy of geostationary orbit satellite
EP2990338A1 (en) * 2014-08-28 2016-03-02 The Boeing Company Satellite transfer orbit search methods
CN108216687A (en) * 2017-12-25 2018-06-29 中国空间技术研究院 GEO satellite based on particle cluster algorithm becomes rail policy calculation method, system and medium
US20200055617A1 (en) * 2018-08-17 2020-02-20 Mitsubishi Electric Research Laboratories, Inc. System and Method of Tracking a Spacecraft Trajectory for Orbital Transfer
CN109625323A (en) * 2018-11-09 2019-04-16 中国科学院空间应用工程与技术中心 A kind of satellite chemical propulsion orbit changing method and system

Also Published As

Publication number Publication date
CN113602531B (en) 2023-05-12

Similar Documents

Publication Publication Date Title
Polk et al. Demonstration of the NSTAR ion propulsion system on the Deep Space One mission
CN109625323A (en) A kind of satellite chemical propulsion orbit changing method and system
DK2586711T3 (en) Method and system for control of a unit of at least two satellites adapted to provide a service
McAdams et al. Trajectory design and maneuver strategy for the MESSENGER mission to Mercury
Kos et al. Altair descent and ascent reference trajectory design and initial dispersion analyses
Martin et al. Saturn V guidance, navigation, and targeting.
Lafleur et al. Low-earth-orbit constellation phasing using miniaturized low-thrust propulsion systems
CN113602531A (en) Method and system for generating GEO orbital transfer strategy of assembly under abnormal separation condition
Cupples et al. Application of Solar Electric Propulsion to a Comet Surface Sample Return Mission
Milligan et al. SMART-1 electric propulsion: an operational perspective
Patha et al. Guidance, energy management, and control of a fixed-impulse solid-rocket vehicle during orbit transfer
Abilleira 2011 Mars Science Laboratory Mission Design Overview
CN112393648A (en) Balance flight theoretical method for autonomous control under rocket thrust failure mode
Garner et al. In-flight operation of the Dawn ion propulsion system through orbit capture at Vesta
McAdams et al. MESSENGER–Six Primary Maneuvers, Six Planetary Flybys, and 6.6 Years to Mercury Orbit
Konstantinov et al. Spacecraft Station-keeping on the Molniya Orbit Us...
CN117208231B (en) GEO satellite minimum orbit height calculation method based on satellite propellant constraint
CN114460952B (en) Double-star cooperative orbit transfer method and system for initializing elliptical orbit flight accompanying configuration
Evdokimov et al. Landing a descent module on the Vostochnyi launch site after returning from the Moon
Gong Application of Celestial Mechanics Theory in Spacecraft Orbit Design
Blumer A future concept of coordinated orbit control of colocated geostationary satellites
Ruschmann et al. Efficient Station-Keeping For Cluster Flight
Kaplan All-electric thruster control of a geostationary communications satellite
Morantea et al. Low-Thrust Trajectory Optimization and Autonomy Analysis for a Medium-Earth-Orbit Constellation Deployment
Vaz et al. Sub-optimal orbital maneuvers for artificial satellites

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