CN112084581A - Spacecraft low-thrust perturbation intersection trajectory optimization method and system - Google Patents

Spacecraft low-thrust perturbation intersection trajectory optimization method and system Download PDF

Info

Publication number
CN112084581A
CN112084581A CN202011013479.4A CN202011013479A CN112084581A CN 112084581 A CN112084581 A CN 112084581A CN 202011013479 A CN202011013479 A CN 202011013479A CN 112084581 A CN112084581 A CN 112084581A
Authority
CN
China
Prior art keywords
thrust
spacecraft
pulse
increment
small
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
CN202011013479.4A
Other languages
Chinese (zh)
Other versions
CN112084581B (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.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN202011013479.4A priority Critical patent/CN112084581B/en
Publication of CN112084581A publication Critical patent/CN112084581A/en
Application granted granted Critical
Publication of CN112084581B publication Critical patent/CN112084581B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • 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
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/04—Constraint-based CAD
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • 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)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Navigation (AREA)

Abstract

The invention discloses a low-thrust perturbation intersection track optimization method and a system, which comprise the following steps: the orbit number of the spacecraft and the target at the starting time and the intersection time is given, and four-pulse speed increment is calculated; assuming that the switching strategy of the small thrust is on-off-on, estimating two starting time lengths of the small thrust; according to the starting time length of the two small thrusts, the middle point moment of the small thrusts is taken as the equivalent pulse moment, and a new four-pulse speed increment is recalculated; outputting the starting time length with small thrust until the pulse increment size change calculated in the two previous and later times is smaller than a preset value; the estimated small-thrust startup duration is used as a constraint and input into an indirect method optimization model for solving, and the optimal control rate, the transfer track and the mass variation are obtained; and (4) calculating the increment percentage of the speed increment and the pulse speed increment corresponding to the optimal control rate, if the increment percentage is greater than the threshold value, solving the step 5 again, and if the increment percentage is less than the threshold value, outputting the optimal control rate and the transfer track.

Description

Spacecraft low-thrust perturbation intersection trajectory optimization method and system
Technical Field
The invention belongs to the technical field of space navigation control, and particularly relates to a method and a system for optimizing a low-thrust perturbation intersection track of a spacecraft.
Background
On the near-earth orbit, the earth non-spherical gravitational perturbation, the atmospheric resistance and the like are the most main perturbation items. When the small-thrust intersection trajectory optimization is carried out during the space debris clearing task, the influence of the perturbation term must be considered, but the corresponding orbit integral calculation amount is obviously increased compared with a two-body analysis dynamic model. Moreover, when the number of transfer turns is large, the calculation amount is greatly increased, and the efficiency of optimizing the rendezvous trajectory is reduced. Therefore, a new intersection trajectory planning method needs to be found, and planning efficiency is improved.
Disclosure of Invention
The invention aims to solve the technical problem of how to improve the optimization efficiency of a low-thrust intersection track when the influence of a perturbation item is considered, and provides a method and a system for quickly optimizing the low-thrust intersection track.
In order to solve the problem, the technical scheme adopted by the invention is as follows:
a spacecraft low thrust perturbation intersection trajectory optimization method comprises the following steps:
step 1: given spacecraft and target at a starting time t0And the number of tracks at the time of the crossing, tfThe track transfer time length is delta t, and four-pulse velocity increment delta v is calculated1,Δv2,Δv3,Δv4The first two pulses in the first and the last two pulses in the last cycle, and the pulse increment of the first cycle is recorded as Deltav0And the last pulse increment is recorded as Deltavf，
Δv0＝Δv1+Δv2
Δvf＝Δv3+Δv4
Step 2: assuming that the low-thrust switching strategy is on-off-on, the first start-up of the low-thrust switching strategy corresponds to a pulse speed increment of Δ v0The second start-up corresponds to a pulse velocity increment of Δ vfAccording to the mass and the thrust of the spacecraft, two equivalent starting time lengths of the small thrust are estimated to be respectively
Figure BDA0002698307390000011
Figure BDA0002698307390000012
Wherein k is the arc segment loss coefficient of the maneuvering effect when the impulse is expanded into small thrust; f is the thrust of the spacecraft, and m is the mass of the spacecraft;
and step 3: according to the estimated equivalent starting time length of the two small thrusts, taking the middle point moment of the small thrusts as the equivalent pulse moment, and taking the middle point moment of the small thrusts as the equivalent pulse moment
Figure BDA0002698307390000021
And
Figure BDA0002698307390000022
equivalent transfer duration of
Figure BDA0002698307390000023
Returning to the step 1 to recalculate the new four-pulse speed increment to obtain the pulse increment delta v of the first circle0' and pulse increment of the last revolution Δ vf'；
And 4, step 4: repeating the steps 1, 2 and 3 until the pulse increment size change calculated in two times is less than a preset value, | (delta v)0'+Δvf')-(Δv0+Δvf) If yes, outputting corresponding small-thrust starting time;
Δ v obtained at this time0'+Δvf' is the speed increment required to meet the estimated small thrust, at this time with the Δ v0'+Δvf' relative toCorresponding delta t1And Δ t2Starting up time length for two estimated small thrusters;
and 5: the estimated small-thrust starting-up time is used as a constraint and input into an optimization model of an indirect method small-thrust intersection track for solving, and the optimal control rate, the transfer track in the aircraft transfer process and the mass variation are obtained;
step 6: calculating the speed increment delta v corresponding to the optimal control rate according to the mass variation in the aircraft transfer processoptThe percentage of increase from the pulse velocity increment in step 4,
Figure BDA0002698307390000024
Figure BDA0002698307390000025
if the optimal control rate is greater than the preset threshold value, the step 5 is solved again, and if the optimal control rate is less than the preset threshold value, the optimal control rate and the transfer track are output, wherein IspIs specific impulse of engine, g is acceleration of gravity, m0Mass m at the initial moment of the spacecraftfAnd (4) meeting the quality of the spacecraft.
The invention also provides a system for quickly optimizing the low-thrust perturbation intersection trajectory of the spacecraft, which comprises a memory and a processor, wherein the memory stores the method for quickly optimizing the low-thrust perturbation intersection trajectory of the spacecraft, and the processor realizes the steps of the method when operating the method for quickly optimizing the low-thrust perturbation intersection trajectory of the spacecraft.
Compared with the prior art, the invention has the following beneficial effects:
according to the method and the system for quickly optimizing the low-thrust perturbation rendezvous trajectory of the spacecraft, the low-thrust rendezvous pulse velocity increment is calculated according to the given initial orbit position and the target orbit position of the spacecraft, and then the low-thrust rendezvous pulse velocity increment is used as an initial value and an optimality criterion for solving and is substituted into an indirect method optimization model to obtain the optimal low-thrust transfer trajectory, so that the problem that the calculation efficiency is low due to the fact that the dynamics recursion calculation amount is large when the transfer time is long in the traditional indirect method low-thrust rendezvous optimization algorithm is solved.
Drawings
FIG. 1 is a flow chart of the system of the present invention;
FIG. 2 is a schematic diagram of the power-on and power-off conditions of low thrust under the optimal solution conditions in the embodiment;
FIG. 3 is a trajectory diagram of the spacecraft under optimal solution conditions in the embodiment;
fig. 4 is a schematic diagram of the startup and shutdown conditions of a small thrust under the condition of a locally optimal solution obtained by using a common indirect method trajectory optimization algorithm in the comparative embodiment.
Detailed Description
As shown in fig. 1, a specific embodiment of a method for quickly optimizing a small thrust perturbation intersection trajectory of a spacecraft according to the present invention is provided, in this embodiment, the number of orbits of the spacecraft and a target at a reference time is shown in table 1:
TABLE 1 reference time spacecraft and target orbit root
Figure BDA0002698307390000031
The initial mass of the spacecraft is 755.32kg, the thrust is 0.5N, and the specific impulse is 1000 s. To be calculated at a relative reference time t0Track transfer began on day 154.1, tf157.8 days of intersection, 3.7 days of Δ t, the optimum control rate for the required small thrust.
Step 1: given spacecraft and target at a starting time t0And the number of tracks at the time of the crossing, tfThe track transfer time is delta t, and the increment of the four-pulse speed is calculated to be delta v1,Δv2,Δv3,Δv4The first two pulses are in the first circle and the last two pulses are in the last circle, and the increment of the first circle is recorded as delta v0And the last pulse increment is recorded as Deltavf，
Figure BDA0002698307390000041
The method for calculating the four-pulse velocity increment comprises a four-pulse perturbation orbit intersection optimization algorithm or a pulse velocity increment estimation method. The pulse velocity increment estimation method is used in the present embodiment.
[Δv1,Δv2,Δv3,Δv4]＝f(Element0,Elementf,tf,Δt) (2)
Wherein Element0And ElementfRespectively a starting orbit and a target orbit of the spacecraft, tfAt the time of the crossing, Δ t is the track transfer duration, t0＝tfAnd- Δ t is the departure time, and in this embodiment, a pulse velocity increment estimation method in a rapid estimation method of the optimal velocity increment of long-time rail transit under J2 perturbation, which is disclosed by the patent number CN110789739A, is used as the function f. Will tfSubstituting the number of orbits into the function f to obtain Δ v, wherein the number of the orbits is 157.8 days, the number of the Δ t is 3.7 days0＝28.4m/s，Δvf＝14.7m/s。
Step 2: assuming that the low-thrust switching strategy is on-off-on, the first start-up of the low-thrust switching strategy corresponds to a pulse speed increment of Δ v0The second start-up corresponds to a pulse velocity increment of Δ vfAccording to the mass and the thrust of the spacecraft, two equivalent starting time lengths of the small thrust are estimated to be respectively
Figure BDA0002698307390000042
Wherein k is the arc segment loss coefficient of the maneuvering effect when the impulse is expanded into small thrust; f is the thrust of the spacecraft, m is the mass of the spacecraft, k is a coefficient less than 1, and k is 0.7 in this embodiment; will delta v0＝28.4m/s，ΔvfSubstituting 14.7m/s into equation 3 yields Δ t1＝0.71，Δt2Day 0.37. Δ t1For the first pulse speed increment at start-up Δ v0Corresponding low thrust equivalent boot-up duration, Δ t2For the second on-time pulse rateIncrement Δ vfThe corresponding small thrust is equivalent to the starting time.
And step 3: according to the estimated small-thrust starting-up time length, taking the small-thrust midpoint moment as an equivalent pulse moment, namely taking the small-thrust midpoint moment as
Figure BDA0002698307390000043
And
Figure BDA0002698307390000044
the equivalent transfer duration is
Figure BDA0002698307390000045
Returning to the step 1 to recalculate the new four-pulse speed increment to obtain the pulse increment delta v of the first circle0' and pulse increment of the last revolution Δ vf'；
In this embodiment, t isf(iii) substitution of formula (1) for 157.8-0.37/2-157.615 days and Δ t-3.7-0.71/2-0.37/2-3.16 days to give Δ ν0'＝32.07m/s，Δvf'＝16.19m/s。
And 4, step 4: repeating the steps 1, 2 and 3 until the change of the pulse speed increment calculated in the two times is smaller than a preset value, and outputting the equivalent starting time with low thrust;
|(Δv0'+Δvf')-(Δv0+Δvf)|＜ (4)
Δ v obtained at this time0'+Δvf' i.e. the pulse velocity increment required for the estimated small thrust intersection, with this Δ v0'+Δvf' corresponding Deltat1And Δ t2The estimated small-thrust equivalent starting time length is obtained;
in this example, 0.1m/s was taken, and after 3 iterations, the final Δ t was obtained1Day 0.82,. DELTA.t20.41 days,. DELTA.v0＝32.69m/s，Δvf＝16.51m/s。
And 5: and (4) taking the estimated small-thrust startup duration as a constraint input to an optimization model for optimizing a small-thrust intersection track by an indirect method for solving to obtain the optimal control rate, the position speed and the mass variation in the transfer process.
The optimization model for optimizing the low-thrust intersection track by an indirect method refers to the following steps:
let the thrust of the spacecraft be F and the initial mass be m0Taking into account the earth's non-spherical gravitational field J2Term, in normalized position r ═ X Y Z]TAnd velocity v ═ vX vY vZ]TThe finite thrust spacecraft dynamics equation for the state quantity is:
Figure BDA0002698307390000051
Figure BDA0002698307390000052
Figure BDA0002698307390000053
Figure BDA0002698307390000054
the fuel consumption model is:
Figure BDA0002698307390000055
wherein the basic dimension units RUnitIs set as the radius of the earth, TUnitSet as 600s, mUnitSet as initial mass m of spacecraft0Then the normalization of other units can be obtained by dimension analysis; alpha ═ alphaX αY αZ]TIs a unit vector of thrust direction, alphaX、αY、αZRespectively representing unit vectors of thrust of the spacecraft in X, Y, Z three-axis coordinate directions, X, Y, Z representing coordinate positions of the spacecraft on three axes in a space coordinate system, vX、vY、vZIndicating that the speed of the spacecraft is X, Y, Z IIIComponent in the direction of axis coordinate, u ∈ [0,1 ]]And u is the amplitude of the thrust force,
Figure BDA0002698307390000061
for the maximum amplitude of the normalized thrust force,
Figure BDA0002698307390000062
to specific fuel consumption, IspIs the engine specific impulse, g is the gravity acceleration, mu is the earth gravity constant, r is the modulus of the vector r; reRepresenting the radius of the earth.
For each section of track intersection, the initial time and the final time of the intersection are fixed, the number of the tracks at the corresponding moment is converted into the position speed, and the boundary value conditions are as follows
Figure BDA0002698307390000063
Wherein, m (t)0) I.e. the mass m of the spacecraft at the initial moment0；
The optimal thrust amplitude u (t) and the optimal thrust direction α (t) need to be solved so that the fuel consumption index is minimal:
Figure BDA0002698307390000064
using the principle of minima, first writing the Hamiltonian
Figure BDA0002698307390000065
Wherein
Figure BDA0002698307390000066
To minimize the Hamiltonian, according to (6)
Figure BDA0002698307390000067
Of the expression (b), the optimum thrust direction being such thatvA takes the minimum value, i.e.
Figure BDA0002698307390000068
At the same time, the optimum thrust amplitude u should be such that
Figure BDA0002698307390000071
Taking the minimum value, namely:
Figure BDA0002698307390000072
where ρ is defined as the switching function:
Figure BDA0002698307390000073
according to the principle of minima, the covariates have the following relationships:
Figure BDA0002698307390000074
unfolding to obtain:
Figure BDA0002698307390000075
Figure BDA0002698307390000076
Figure BDA0002698307390000077
Figure BDA0002698307390000078
at the same time, as the terminal quality is not constrained, then
λm(tf)＝0 (17)
If an initial covariant lambda is givenr(t0)，λv(t0)，λm(t0) And initial state r (t)0)，v(t0)，m(t0) Then the state and the covariates of the terminal can be calculated by differential equations. Optimal control problem is equivalent to solving the initial covariate λr(t0)，λv(t0)，λm(t0) Such that the end state and covariates satisfy the targeting function phi (z):
φ(z)＝[r(tf)-rf,v(tf)-vf,λm(tf)]T＝0 (18)
that is, the target function has 7 equations in total, and the solution variable is z ═ λr(t0),λv(t0),λm(t0)]Also 7 dimensions, corresponding to 7 covariates at the initial instant.
The number of the spacecraft orbits at the departure time and the target orbits at the rendezvous time is respectively converted into unitized position velocity r0＝[0.425410721120000 0.540281353102000 -0.895982743973000]，v0＝[0.467147840028000 0.316148727267000 0.413824368484000]And rf＝[-0.167451008984000 0.0699372107820000 -1.10700302437600]，vf＝[0.497209586529000 0.493710947528000 -0.0434604813910000]。m0Is the initial mass. I.e. known r0,v0,m0Solving for λr(t0),λv(t0),λm(t0) Let the system of differential equations (6), (7), (15) at tfThe time satisfies the terminal condition expression (18).
According to Δ t1And Δ t2Constraining thrust amplitude u (t) at [ t0+Δt1,tf-Δt2]The interval is always 0, and a target practice is performed by randomly guessing the covariates using a published nonlinear equation set solution and optimization toolkit Minpack-1.
Solving to obtain the optimal solution
[λr(t0),λv(t0),λm(t0)]＝[0.283555391709000，-0.0590384733860000，-0.336635442397000，0.625749386435000，-0.148776401651000，0.280129358394000，0.00319735301800000]。
Step 6: calculating the corresponding speed increment delta v of the optimal track according to the mass variation in the aircraft transfer processoptThe percentage of increase from the pulse velocity increment in step 4,
according to the optimal solution [ lambda ]r(t0),λv(t0),λm(t0)]And r (t)0)，v(t0)，m(t0) Then, t can be obtained by integrating equations (6) and (7)0To tfThe optimal track between the position r, the speed v and the change relation of the mass m along with the time is obtained according to the value of the mass m along with the change of the time0And mfIs m (t)0) And m (t)f) Equivalent velocity increment Deltav corresponding to optimal trackopt：
Figure BDA0002698307390000091
Calculating to obtain the equivalent velocity increment delta v corresponding to the optimal trackopt＝46.3m/s，
Δ v estimated in step 40'+Δvf'＝32.69+16.51＝49.2
Then
Figure BDA0002698307390000092
In this embodiment, the preset threshold is 10%, and since the threshold is less than 10%, the solution obtained in step 5 is the optimal solution, and the optimal solution is output.
The thrust on-off condition of the optimal solution is shown in figure 2, the track is shown in figure 3, and the fuel consumption is 3.4 kg. Optimization takes less than 2 minutes. If the solving strategy designed by the invention is not used, the local optimal solution can be obtained with a higher probability by a common indirect method track optimization algorithm, the thrust on-off condition of the local optimal solution solved at a certain time is shown in figure 4, and the corresponding fuel consumption is 6.4 kg.
The embodiment proves that the method has high optimization speed, can avoid local optimal solution, and solves the problem of quick optimization of the small-thrust perturbation intersection track.
The invention also provides a system for quickly optimizing the low-thrust perturbation intersection trajectory of the spacecraft, which comprises a memory and a processor, wherein the memory stores the method for quickly optimizing the low-thrust perturbation intersection trajectory of the spacecraft, and the processor realizes the steps of the method when operating the method for quickly optimizing the low-thrust perturbation intersection trajectory of the spacecraft.
The above-mentioned embodiments only express a certain implementation mode of the present application, and the description thereof is more specific and detailed, but not construed as limiting the scope of the invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the concept of the present application, which falls within the scope of protection of the present application. Therefore, the protection scope of the present patent shall be subject to the appended claims.

Claims (5)

1. A spacecraft low-thrust perturbation intersection trajectory optimization method is characterized by comprising the following steps:
step 1: given spacecraft and target at a starting time t0And the number of tracks at the time of the crossing, tfThe track transfer time length is delta t, and four-pulse velocity increment delta v is calculated1,Δv2,Δv3,Δv4The first two pulses in the first and the last two pulses in the last cycle, and the pulse increment of the first cycle is recorded as Deltav0And the last pulse increment is recorded as Deltavf，
Δv0＝Δv1+Δv2
Δvf＝Δv3+Δv4
Step 2: assuming the low thrust switching strategy is on-off-on, the first power-on of the low thrust switching strategyCorresponding pulse velocity increment of Δ v0The second start-up corresponds to a pulse velocity increment of Δ vfAccording to the mass and the thrust of the spacecraft, two equivalent starting time lengths of the small thrust are estimated to be respectively
Figure FDA0002698307380000011
Figure FDA0002698307380000012
Wherein k is the arc segment loss coefficient of the maneuvering effect when the impulse is expanded into small thrust; f is the thrust of the spacecraft, and m is the mass of the spacecraft;
and step 3: according to the estimated starting time length of the two equivalent small thrusts, taking the middle point moment of the small thrusts as the equivalent pulse moment, and taking the middle point moment of the small thrusts as
Figure FDA0002698307380000013
And
Figure FDA0002698307380000014
equivalent transfer duration of
Figure FDA0002698307380000015
Returning to the step 1 to recalculate the new four-pulse speed increment to obtain the pulse increment delta v of the first circle0' and pulse increment of the last revolution Δ vf'；
And 4, step 4: repeating the steps 1, 2 and 3 until the change of the pulse velocity increment of the two previous and next calculations is less than a preset value, | (delta v)0'+Δvf')-(Δv0+Δvf) I <, outputting corresponding small-thrust equivalent starting time;
Δ v obtained at this time0'+Δvf' i.e. the pulse velocity increment required for the estimated small thrust intersection, with this Δ v0'+Δvf' corresponding Deltat1And Δ t2To estimateThe counted low thrust is equivalent to the starting time;
and 5: taking the estimated low-thrust equivalent starting time as a constraint input into an optimization model of an indirect low-thrust intersection track to solve, and obtaining the optimal control rate, the transfer track in the aircraft transfer process and the mass variation;
step 6: calculating the corresponding speed increment delta v of the optimal track according to the mass variation in the aircraft transfer processoptThe percentage of increase from the pulse velocity increment in step 4,
Figure FDA0002698307380000021
Figure FDA0002698307380000022
if the optimal control rate is greater than the preset threshold value, the step 5 is solved again, and if the optimal control rate is less than the preset threshold value, the optimal control rate and the transfer track are output, wherein IspIs specific impulse of engine, g is acceleration of gravity, m0Mass m at the initial moment of the spacecraftfAnd (4) meeting the quality of the spacecraft.
2. The method of claim 1, wherein: the optimization model of the indirect low-thrust intersection trajectory in the step 5 refers to:
let the thrust of the spacecraft be F and the initial mass be m0Taking into account the earth's non-spherical gravitational field J2Term, in normalized position r ═ X Y Z]TAnd velocity v ═ vX vY vZ]TThe finite thrust spacecraft dynamics equation for the state quantity is:
Figure FDA0002698307380000023
Figure FDA0002698307380000024
Figure FDA0002698307380000025
Figure FDA0002698307380000026
the fuel consumption model is:
Figure FDA0002698307380000027
wherein the basic dimension units RUnitIs set as the radius of the earth, TUnitSet as 600s, mUnitSet as initial mass m of spacecraft0Then the normalization of other units can be obtained by dimension analysis; alpha ═ alphaX αY αZ]TIs a unit vector of thrust direction, alphaX、αY、αZRespectively representing unit vectors of thrust of the spacecraft in X, Y, Z three-axis coordinate directions, X, Y, Z representing coordinate positions of the spacecraft on three axes in a space coordinate system, vX、vY、vZRepresents the component of the spacecraft velocity in the X, Y, Z three-axis coordinate direction, and u is equal to 0,1]And u is the amplitude of the thrust force,
Figure FDA0002698307380000031
for the maximum amplitude of the normalized thrust force,
Figure FDA0002698307380000032
mu is the Earth's gravitational constant, R is the modulus of the vector R, R is the specific fuel consumptioneRepresents the radius of the earth;
for each section of track intersection, the initial time and the final time of the intersection are fixed, the number of the tracks at the corresponding moment is converted into the position speed, and the boundary value conditions are as follows
Figure FDA0002698307380000033
The optimal thrust amplitude u (t) and thrust direction α (t) need to be solved so that the fuel consumption index is minimal:
Figure FDA0002698307380000034
using the principle of minima, first writing the Hamiltonian
Figure FDA0002698307380000035
Wherein
Figure FDA0002698307380000036
To minimize the Hamiltonian, according to (6)
Figure FDA0002698307380000039
Of the expression (b), the optimum thrust direction being such thatvA takes the minimum value, i.e.
Figure FDA0002698307380000037
At the same time, the optimum thrust amplitude u should be such that
Figure FDA0002698307380000038
Taking the minimum value, namely:
Figure FDA0002698307380000041
where ρ is defined as the switching function:
Figure FDA0002698307380000042
according to the principle of minima, the covariates have the following relationships:
Figure FDA0002698307380000043
unfolding to obtain:
Figure FDA0002698307380000044
Figure FDA0002698307380000045
Figure FDA0002698307380000046
wherein
Figure FDA0002698307380000047
At the same time, as the terminal quality is not constrained, then
λm(tf)＝0 (17)
If an initial covariant lambda is givenr(t0)，λv(t0)，λm(t0) And initial state r (t)0)，v(t0)，m(t0) Then the state and the covariates of the terminal can be obtained by calculation through a differential equation, and the optimal control problem is equivalent to solving the initial covariates lambdar(t0)，λv(t0)，λm(t0) So that the end state and the covariant state changeThe quantities satisfy the targeting function phi (z):
φ(z)＝[r(tf)-rf,v(tf)-vf,λm(tf)]T＝0 (18)
that is, the target function has 7 equations in total, and the solution variable is z ═ λr(t0),λv(t0),λm(t0)]Also 7 dimensions, corresponding to 7 covariates at the initial instant.
3. The method of claim 2, wherein the method of solving the optimization model in step 5 is:
according to Δ t1And Δ t2Constraining thrust amplitude u (t) at [ t0+Δt1,tf-Δt2]The interval is always 0, and a target practice is performed by randomly guessing the covariates using a published nonlinear equation set solution and optimization toolkit Minpack-1.
4. The method of claim 1, wherein: the calculation method of the pulse velocity increment in step 1 is an optimal four-pulse transfer optimization method or a four-pulse transfer estimation method.
5. The utility model provides a spacecraft low thrust perturbation intersection orbit optimization system which characterized in that: the method comprises a memory and a processor, wherein the memory stores a method for quickly optimizing the low-thrust perturbation intersection trajectory of the spacecraft, and the processor realizes the steps of the method in any one of claims 1 to 4 when running the method for quickly optimizing the low-thrust perturbation intersection trajectory of the spacecraft.
CN202011013479.4A 2020-09-24 2020-09-24 Spacecraft low-thrust perturbation intersection trajectory optimization method and system Active CN112084581B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011013479.4A CN112084581B (en) 2020-09-24 2020-09-24 Spacecraft low-thrust perturbation intersection trajectory optimization method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011013479.4A CN112084581B (en) 2020-09-24 2020-09-24 Spacecraft low-thrust perturbation intersection trajectory optimization method and system

Publications (2)

Publication Number Publication Date
CN112084581A true CN112084581A (en) 2020-12-15
CN112084581B CN112084581B (en) 2022-08-19

Family

ID=73739688

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011013479.4A Active CN112084581B (en) 2020-09-24 2020-09-24 Spacecraft low-thrust perturbation intersection trajectory optimization method and system

Country Status (1)

Country Link
CN (1) CN112084581B (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113353289A (en) * 2021-04-25 2021-09-07 北京控制工程研究所 Autonomous driving and separating method and device for space game and storage medium
CN114275191A (en) * 2021-02-26 2022-04-05 北京空间飞行器总体设计部 Rail-controlled speed increment estimation method suitable for lunar rail intersection docking task
CN114384803A (en) * 2022-01-12 2022-04-22 北京理工大学 Low-thrust optimal orbital transfer method considering terrain shadow constraint
CN114384806A (en) * 2022-01-12 2022-04-22 北京理工大学 Segmented optimization method for multi-circle orbital transfer of electric propulsion spacecraft under multi-perturbation earth shadow constraint
CN114462256A (en) * 2022-04-11 2022-05-10 中国人民解放军国防科技大学 Method, device, equipment and medium for determining non-cooperative low-thrust maneuvering target track
CN114485620A (en) * 2022-01-29 2022-05-13 中国科学院国家空间科学中心 Orbital dynamics fused asteroid detector autonomous visual positioning system and method
CN114491819A (en) * 2022-02-16 2022-05-13 哈尔滨逐宇航天科技有限责任公司 Rocket carrying capacity intelligent solving method based on speed loss calculation

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110789739A (en) * 2019-11-08 2020-02-14 中国人民解放军国防科技大学 Method for quickly estimating optimal speed increment of long-time rail crossing under J2 perturbation
CN111268176A (en) * 2020-01-17 2020-06-12 中国人民解放军国防科技大学 Perturbation track four-pulse intersection rapid optimization method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110789739A (en) * 2019-11-08 2020-02-14 中国人民解放军国防科技大学 Method for quickly estimating optimal speed increment of long-time rail crossing under J2 perturbation
CN111268176A (en) * 2020-01-17 2020-06-12 中国人民解放军国防科技大学 Perturbation track four-pulse intersection rapid optimization method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
周军等: "考虑地球扁率J2摄动影响的异面椭圆轨道多冲量最优交会", 《宇航学报》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114275191A (en) * 2021-02-26 2022-04-05 北京空间飞行器总体设计部 Rail-controlled speed increment estimation method suitable for lunar rail intersection docking task
CN113353289A (en) * 2021-04-25 2021-09-07 北京控制工程研究所 Autonomous driving and separating method and device for space game and storage medium
CN113353289B (en) * 2021-04-25 2022-12-13 北京控制工程研究所 Autonomous driving and separating method and device for space game and storage medium
CN114384803A (en) * 2022-01-12 2022-04-22 北京理工大学 Low-thrust optimal orbital transfer method considering terrain shadow constraint
CN114384806A (en) * 2022-01-12 2022-04-22 北京理工大学 Segmented optimization method for multi-circle orbital transfer of electric propulsion spacecraft under multi-perturbation earth shadow constraint
CN114485620A (en) * 2022-01-29 2022-05-13 中国科学院国家空间科学中心 Orbital dynamics fused asteroid detector autonomous visual positioning system and method
CN114485620B (en) * 2022-01-29 2023-07-28 中国科学院国家空间科学中心 Autonomous visual positioning system and method for asteroid detector fused with orbit dynamics
CN114491819A (en) * 2022-02-16 2022-05-13 哈尔滨逐宇航天科技有限责任公司 Rocket carrying capacity intelligent solving method based on speed loss calculation
CN114491819B (en) * 2022-02-16 2022-08-02 哈尔滨逐宇航天科技有限责任公司 Rocket carrying capacity intelligent solving method based on speed loss calculation
CN114462256A (en) * 2022-04-11 2022-05-10 中国人民解放军国防科技大学 Method, device, equipment and medium for determining non-cooperative low-thrust maneuvering target track
CN114462256B (en) * 2022-04-11 2022-06-21 中国人民解放军国防科技大学 Method, device, equipment and medium for determining non-cooperative low-thrust maneuvering target track

Also Published As

Publication number Publication date
CN112084581B (en) 2022-08-19

Similar Documents

Publication Publication Date Title
CN112084581B (en) Spacecraft low-thrust perturbation intersection trajectory optimization method and system
Bombardelli et al. Asymptotic solution for the two-body problem with constant tangential thrust acceleration
Petukhov Method of continuation for optimization of interplanetary low-thrust trajectories
CN110032768B (en) Four-pulse orbit intersection optimization method using accurate dynamic model
WO2017005052A1 (en) Optimization and design method for gradient segmentation of intervals of spacecraft pulse rendezvous trajectory
CN111444603B (en) Method for rapidly planning shortest time off-orbit trajectory of recoverable spacecraft
Petukhov Application of the angular independent variable and its regularizing transformation in the problems of optimizing low-thrust trajectories
Luo et al. Optimization of multiple-impulse minimum-time rendezvous with impulse constraints using a hybrid genetic algorithm
CN115841051B (en) Method and storage medium based on space vehicle orbit estimation
CN113110561A (en) Random distribution robust optimization-based SMPC algorithm for maintaining formation of satellite
CN111127295B (en) SGP4 orbit model integrated parallel method based on GPU
Benedikter et al. Convex optimization of launch vehicle ascent trajectory with heat-flux and splash-down constraints
Sauer Jr Solar electric performance for medlite and delta class planetary missions
Wang et al. Six-DOF trajectory optimization for reusable launch vehicles via Gauss pseudospectral method
Yang et al. An aerodynamic shape optimization study to maximize the range of a guided missile
Petropoulos et al. Low-thrust transfers using primer vector theory and a second-order penalty method
CN107315357B (en) Approximate modeling method of rod-shaped flexible mechanism for adjusting attitude of chip star
Zhang et al. Semi-analytical approach for computing near-optimal low-thrust transfers to geosynchronous orbit
Gao et al. Optimization of low-thrust many-revolution transfers and Lyapunov-based guidance
Songyan et al. Robust spline-line energy management guidance algorithm with multiple constraints and uncertainties for solid rocket ascending
CN112257259B (en) Method and system for estimating whole-course trajectory of ballistic missile based on improved autonomous multiple models
CN110569625B (en) Method and system for quickly generating deep-space low-thrust flying-over orbit
CN111178589A (en) Iterative guidance method and system based on gray prediction
Nosratolahi et al. Design of satellite’s combined attitude and thermal control system equipped with FMC actuators
Sun et al. SGCMG non-singularity steering based on adaptive Gauss pseudospectral method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant