CN106446313A - Polar hovering satellite orbit-based system design method - Google Patents

Polar hovering satellite orbit-based system design method Download PDF

Info

Publication number
CN106446313A
CN106446313A CN201510487075.1A CN201510487075A CN106446313A CN 106446313 A CN106446313 A CN 106446313A CN 201510487075 A CN201510487075 A CN 201510487075A CN 106446313 A CN106446313 A CN 106446313A
Authority
CN
China
Prior art keywords
satellite
polar region
hovering
orbit
design
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
CN201510487075.1A
Other languages
Chinese (zh)
Other versions
CN106446313B (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.)
Beijing Institute of Spacecraft System Engineering
Original Assignee
Beijing Institute of Spacecraft System 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 Beijing Institute of Spacecraft System Engineering filed Critical Beijing Institute of Spacecraft System Engineering
Priority to CN201510487075.1A priority Critical patent/CN106446313B/en
Publication of CN106446313A publication Critical patent/CN106446313A/en
Application granted granted Critical
Publication of CN106446313B publication Critical patent/CN106446313B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

The invention provides a polar hovering satellite orbit-based system design method. The system design and transfer orbit design of a whole satellite are finished by performing design and characteristic analysis on a working orbit of the polar hovering satellite. The method comprises the steps of describing a dynamic motion of the polar hovering satellite and calculating control characteristics of the polar hovering satellite, thereby performing design optimization on the working orbit of the satellite to obtain optimization design parameters; calculating attributes such as mass and the like of subsystems as polar hovering satellites according to the optimization design parameters to obtain maximum mass of an effective load of the polar hovering satellite; and determining a transfer section orbit and a launch section orbit of the transfer orbit of the polar hovering satellite, and optimizing obtained orbit parameters, thereby designing an optimal transfer orbit. According to the method, the scheme is simple in principle and easy to solve; and a calculation result is of great guide significance for project realization and can be applied to a task design stage of the polar hovering satellite with demands of polar observation, communication and the like.

Description

Design method based on polar region hovering satellite orbit
Technical field
The invention belongs to spacecraft dynamics controls and application, it is related to a kind of dynamics Controlling of polar region hovering satellite and application Scheme and in particular to a kind of based on polar region hover satellite orbit design method, by polar region hover satellite work rail Road is designed completing system design and the transfer orbit design of whole star with specificity analysises.
Background technology
Polar region hovering satellite can provide the Continuous Observation to high latitude area especially arctic regions, can be to the sky of arctic regions Gas forecast, earth resource detection, radio communication and Space environment monitor etc. provide service, are the Arctic Ocean and Antarctic Continent exploitation Using providing space-based to support, there is huge application prospect and economic benefit.
Diver, J.M. propose the concept of " polar region hovering " earliest, and the satellite mechanical characteristic of fixing hovering height have been carried out point Analysis, has obtained in the minimum conclusion of perturbative force at about 380 earth radius of the earth.Colin R.McInnes and Matteo Ceriotti etc. devises the optimum polar region hovering track of fuel using the hybrid propulsion mode of electric propulsion and solar sail, has carried out fuel Consumption calculations and durability analysis, and calculate the subsystem quality of whole star, complete under particular job height and life conditions Payload mass calculates.
However, the domestic correlational study also not launching the hovering of relevant polar region.At present, also do not have any country or tissue emissions or Tested polar region hovering task satellite.In the prior art, the research of this respect is lacked from task layer in the face of polar region hovering satellite Transmitting, transfer, work overall process be designed, therefore, it is necessary to a kind of polar region of development hovering Orbit Design and system Design, is designed and specificity analysises by the working track of the satellite that hovered in polar region, so complete the design of whole star system and Can transfer orbit designs, be conducive to assessment prior art condition realize polar region hovering task.
Content of the invention
Carry out from the transmitting of polar region hovering satellite, transfer, work overall process faced by task layer for lacking in the research of prior art The present situation of design and problem, the present invention proposes a kind of polar region hovering satellite orbit and system design scheme, thus complete polar region hanging Stop Orbit Design, including working track design and transfer orbit design, and then complete satellite subsystem design and emission system Design.
Roughly, the orbit Design flow process of polar region hovering satellite of the present invention and system design flow particular content include following part.
(1) hovering working track is designed with circular Restricted three-body model and describes kinetics equation;The fixing hovering of analysis is high The mechanical characteristic of polar region hovering satellite under the conditions of degree;Working track design optimization problem is converted to optimal control problem;Using Excellent control algolithm solves the problems, such as the working track design optimization that hovers.
(2) satellite subsystem design according to the mechanical characteristic of working track design satellite respectively support the power of system, size, The attributes such as quality;According to satellite gross mass, calculate the biggest quality of payload;When requiring to payload, Satellite gross mass is iterated design.
(3) method using manifold design for the transfer orbit design of hovering reverse integration to track, finds out trajectory transfer orbit; Orbit parameter is changed between circular Restricted three-body model and two body Models;Design hovering satellite using Huo Man branch mode Launch trajectory;Using grid data service to design parameter optimization, the optimum transfer orbit of design.
Specifically, the invention provides a kind of design method of the satellite orbit that hovered based on polar region, defended by hovering to polar region The working track of star is designed completing system design and the transfer orbit design of whole star with specificity analysises.The method includes following Step:Step one, the dynamical motion of the satellite that hovered in polar region is described, and the control characteristic of the satellite that hovered in polar region is carried out Calculate, thus being designed to its working track optimizing to obtain optimal design parameter;Step 2, according to optimal design parameter, The attributes such as the quality of each subsystem to the satellite that hovers as polar region calculate, thus obtaining polar region hovering Satellite Payloads The biggest quality;And step 3, determine the transfer leg track of transfer orbit and the launch trajectory of polar region hovering satellite respectively, and The orbit parameter being obtained is optimized, thus designing the transfer orbit of optimum.
Step one executes:Using circular Restricted three-body model, the dynamical motion of the satellite that hovered in polar region is described; Under the conditions of fixing hovering height, the control characteristic of the satellite that hovered in polar region is analyzed;And using optimal control problem Mode, the working track of the satellite that hovered in polar region is designed optimizing.
In addition, step one also includes:Calculate optimal control problem using optimum control numerical optimization software kit, defended with polar region hovering The analysis result of the control characteristic of star is as the initial value of optimized algorithm, thus obtaining the working track of optimal-fuel and controlling acceleration Degree and the time history of controling power.
Polar region hover satellite dynamical motion equation be:
Wherein, dynamical motion defines in the rotated coordinate system, and zero points to ground in Solar-terrestrial system barycenter, x-axis for the sun Direction of bowl, z-axis is Solar-terrestrial system rotational angular velocity vector direction, and y-axis constitutes right hand rectangular coordinate system, and r is under rotating coordinate system The position vector of base hovering satellite, v is the velocity of base hovering satellite, and w is the angular velocity vector of rotating coordinate system, aT For controlling acceleration and its mould aTRepresent, m is the quality of base hovering satellite, IspFor electric propulsion engine specific impulse, g0For Sea level acceleration of gravity.
U is potential-energy function and form is:
Wherein, μ1、μ2It is respectively the gravitational constant of the sun and the earth, r1、r2It is respectively the position vector of the sun and the earth.
Optimal control problem at least includes:Object function, boundary condition and path constraints, wherein, object function is: J=-mf=-m (tf), boundary condition is:ry(t0)=0, m (t0)=m0, r (tf)=r (t0), v (tf)=v (t0), and path constraint Condition is:J represents target function value, and polar region hovering track has the periodicity of 1 year, t0For initial time and elect first point of Capricornus as, r (t0) for polar region hover satellite initial time position vector, v (t0) it is initial time Velocity, m (t0) for polar region hover satellite initial time quality, tfFor task service life, r (tf)、v(tf) hang for polar region Stop the position and speed vector in the last moment of satellite, m (tf) be the last moment quality, and (r2, x, r2, y, r2, z)T=r-r2, δeq=23.5 ° is ecliptic obliquity.
Step 2 executes:Mechanical characteristic according to working track come design polar region hover satellite subsystem attribute;According to The gross mass of polar region hovering satellite, calculates the biggest quality of payload system;And when payload system requires, right The gross mass of polar region hovering satellite is iterated, and wherein, attribute at least includes power, size, quality.
Specifically, step 2 includes:According to controling power time history, find the maximum of controling power;According to Optimal State variable, Obtain the mass change of polar region hovering satellite;Calculate fuel, fuel tank system, electric propulsion system, the matter of power-supply system respectively Amount, thus obtain the quality of payload system;And regard task situation, the first of polar region hovering satellite is recalculated by iteration Prothyl amount.
The transfer orbit of polar region hovering satellite is divided into transmitter section and transfer leg, and wherein, transmitter section is by rocket Upper Stage chemical propulsion Mode completes, and transfer leg is completed using the trajectory branch mode not applying pulse or continuous thrust.
Step 3 executes:Using manifold design, the track of polar region hovering satellite is inversely integrated, thus finding out transfer leg Track, and the orbit parameter of transfer leg track is changed between circular Restricted three-body model and two body Models;Using Huo Manzhuan Shifting mode come to design polar region hover satellite launch trajectory;And utilize grid data service, the orbit parameter being obtained is carried out Optimize, thus designing the transfer orbit of optimum.
In step 3, the process finding out transfer leg track includes:Set out with a point of the polar region hovering track of any time, Turn over predetermined angular with rotating coordinate system, according to the restricted three body dynamics models of circle to track reverse integration, thus obtaining pole The track manifold of ground hovering satellite;Block in the nearest position of the liftoff ball of manifold, obtain reverse time of integration now and herein Position and speed parameter;Position and speed parameter in rotating coordinate system in circular Restricted three-body model is obtained disome through Coordinate Conversion Position and speed under inertial coodinate system in model, and it is calculated corresponding Keplerian orbit six key element;And according to Kepler's rail The perigee altitude of road six element factor calculation transfer orbit and altitude of the apogee.
Therefore, the polar region hovering Orbit Design of the present invention and system design scheme principle are simply it is easy to solve, result of calculation There is very big directive significance for Project Realization, and the method for the present invention can be applicable to the polar region of the demands such as polar region observation, communication The hovering satellite task design phase.
Brief description
Fig. 1 is the schematic diagram of hovering working track in polar region involved in the present invention;
Fig. 2 is the schematic diagram of hovering satellite transfer orbit in polar region involved in the present invention;And
Fig. 3 is the figure cluster of the manifold of the reverse integration of polar region hovering satellite trajectory involved in the present invention.
Specific embodiment
It will be appreciated that hovering Orbit Design in polar region of the present invention is comprised the following steps with design method:
Step one, polar region hovering working track design.
Using circular Restricted three-body model (Circular Restricted Three-Body Problem, hereinafter referred to as CR3BP) The dynamical motion of description polar region hovering satellite;Under the conditions of fixing hovering height, the control characteristic of the satellite that hovered in polar region is carried out Calculate;Polar region hovering working track design optimization problem is described as the optimal control problem of following form:
Object function is J=-mf=-m (tf);
Kinetics equation is using circular Restricted three-body model;
Boundary condition is ry(t0)=0, m (t0)=m0, r (tf)=r (t0), v (tf)=v (t0);
Path constraints are:
Wherein, J represents target function value;Polar region hovering track has the periodicity of 1 year, selects initial time t0For Winter Solstice, r (t0) For spacecraft initial time position vector, v (t0) it is initial time velocity, m (t0) it is initial time quality, tfFor the task longevity Life, r (tf)、v(tf) for spacecraft end the moment position and speed vector, m (tf) be the last moment quality. (r2, x, r2, y, r2, z)T=r-r2, δeq=23.5 ° is ecliptic obliquity.
Then, calculate above-mentioned optimal control problem using optimum control numerical optimization software kit PSOPT, with fixing hovering height Result of calculation, as the initial value of optimized algorithm, obtains the working track of polar region hovering satellite optimal-fuel, obtains controlling acceleration simultaneously Degree and the time history of controling power.
Step 2, polar region hovering satellite subsystem design.
The controling power time history being obtained according to above-mentioned steps one, finds the maximum of T of controling powermax.Obtained according to step one Optimal State variable, obtains mass change Δ m, calculates fuel, fuel tank system, electric propulsion system, power-supply system respectively Quality.Finally, obtain the quality of payload system.If mission requirements require to payload mass, by repeatedly In generation, designs initial mass m of polar region hovering satellite of rerunning0.
Step 3, polar region hovering satellite transfer orbit design.
Polar region hovering satellite transfer orbit transfer leg mentality of designing be:With a little setting out of any time polar region hovering track, with rotation Turn coordinate system and turn over angle, θ and represent, according to the restricted three body dynamics models of circle to track reverse integration, obtain the rail of satellite Mark manifold.Block in the nearest position of the liftoff ball of manifold, obtaining the now reverse time of integration is tT, and position and speed ginseng herein Number.Position and speed parameter in rotating coordinate system in circular Restricted three-body model is obtained in two body motion model through Coordinate Conversion Position and speed under inertial coodinate system, is calculated corresponding Keplerian orbit six key element.According to Keplerian orbit six element factor calculation The perigee altitude of transfer orbit and altitude of the apogee.
Wherein, hovering satellite transfer orbit transmitter section in polar region is designed with the mode of Huo Man dipulse transfer, solves by LEO (Low Earth Orbit, hereinafter referred to as LEO) parking orbit is to the impulse magnitude Δ V of big ellipse target Orbit Transformation1、ΔV2. According to polar region hovering satellite operation track initial mass mtargetCalculate carrier rocket in LEO parking orbit and launch gross mass of entering the orbit mpark.
Being solved using grid-search algorithms makes mparkMinimum variable [tT, θ], complete the design of polar region hovering satellite transfer orbit.
The present invention is described in detail for 1-3 and specific embodiment below in conjunction with the accompanying drawings.
Hover task Orbit optimization design for polar region
The design of this part mainly divides several steps to complete, and is respectively:
(1) using circular Restricted three-body model, the polar region hovering spacecraft equation of motion is described as shown in equation below (1), Concrete such as Fig. 1.
The above-mentioned equation of motion defines in the rotated coordinate system, and rotating coordinate system is defined as follows:Zero in Solar-terrestrial system barycenter, X-axis refers to earthwards for the sun, and z-axis is Solar-terrestrial system rotational angular velocity vector direction, and y-axis constitutes right hand rectangular coordinate system. Wherein, r is the position vector of spacecraft under rotating coordinate system, and v is the velocity of spacecraft, and w is the angle speed of rotating coordinate system Degree vector, aTFor controlling acceleration, its mould aTRepresent, m is spacecraft mass, IspFor electric propulsion engine specific impulse, g0For Sea level acceleration of gravity.Shown in the form of potential-energy function U such as formula (2).
Wherein, μ1、μ2It is respectively the gravitational constant of the sun and the earth, r1、r2It is respectively the position vector of the sun and the earth.
(2) under the conditions of fixing hovering height, the control characteristic of the satellite that hovered in polar region is analyzed.
Under the conditions of fixing hovering height, position vector meets following constraint:
Wherein, d hovers highly for polar region, δeqFor ecliptic obliquity, ω is Solar-terrestrial system angular velocity of rotation, and μ is constant μ=3.0404 × 10-6, AU is astronomical unit.
Calculate any time tiControl acceleration:
Assume to control acceleration in [ti, ti+ Δ t] interval keep constant, mass change is
Repeat above formula (3), (4), the calculating of (5), can obtain working life fuel consumption with control acceleration when Between course.
(3) polar region hovering working track design optimization problem is described as optimal control problem.
Object function is:J=-mf=-m (tf) (6)
Kinetics equation adopts formula (1)
Boundary condition is:ry(t0)=0, m (t0)=m0, r (tf)=r (t0), v (tf)=v (t0) (7)
Path constraints are:
Wherein, J represents target function value;There is the periodicity of 1 year in polar region hovering track, select initial time t0For first point of Capricornus; r(t0) it is spacecraft initial time position vector, v (t0) it is initial time velocity, m (t0) it is initial time spacecraft mass; tfFor task service life, r (tf)、v(tf) for spacecraft end the moment position and speed vector, m (tf) be the last moment quality. (r2, x, r2, y, r2, z)T=r-r2;δeq=23.5 ° is ecliptic obliquity.
(4) optimum control numerical optimization software kit PSOPT is utilized to calculate above-mentioned optimal control problem, with above-mentioned steps (2) Result of calculation, as the initial value of optimized algorithm, obtains the working track of polar region hovering satellite optimal-fuel, obtains controlling acceleration simultaneously Degree and the time history of controling power.
Polar region hovering satellite subsystem design
According to " hover task Orbit optimization design for polar region " the controling power time history that process obtains, find the maximum of controling power Tmax.
According to " hover task Orbit optimization design for polar region " the Optimal State variable that process obtains, obtain the spacecraft matter in last moment Amount mf, then spacecraft fuel mass be:
mprop=m0-mf(9)
Calculate quality m of fuel tanktank
mtank=0.1mprop(10)
Calculate the maximum service rating P in electric propulsion system work processmax
Wherein, thruster efficiency etaSEP=0.7.
The quality of electric propulsion thruster and peak power P of electric propulsion systemmaxRelated:
mSEP=kSEPPmax(12)
Wherein, kSEP=20kg/kW.
It is mainly solar array quality for the power-supply system that electric propulsion system is powered to be about:
mPower=kSAPmax(13)
Wherein
The calculating of payload mass:
mpayload=m0-(mprop+mtank+nthrustersmSEP+mPower) (14)
Wherein, nthrustersNumber for electric propulsion engine, it is considered to Redundancy Design, typically takes nthrusters=2.
mpayloadIt is m for initial mass0Under the conditions of payload capacity quality, if to mpayloadRequire, then need to m0 It is iterated designing:
Subscript " ' " represent new design load.
Polar region hovering transfer orbit design
Polar region hovering transfer orbit can be divided into transmitter section and transfer leg.As shown in Fig. 2 wherein transmitter section is typically by rocket Upper Stage The mode of chemical propulsion completes, and transfer leg can take trajectory transfer to complete, and does not apply pulse or continuous thrust.
(1) transfer leg trajectory transfer design
Under the restricted circular three body dynamics models, sometime putting by hovering track, turned over rotating coordinate system Angle, θ represents, carries out track to hovering spacecraft and inversely integrates, can obtain a series of manifolds as shown in Figure 3.By manifold Block in the place nearest apart from the earth, the reverse time of integration is tT, obtain the preliminary orbit parameter of transfer leg (x1, y1, z1,).By the state vector (x in restricted circular three body dynamics model rotating coordinate systems1, y1, z1,) Be converted to the position and speed vector under two body Model inertial coodinate systems, be reconverted into Keplerian orbit six key element (a0, e0, i0, Ω0, w0, f0).
Calculate perigee and the altitude of the apogee of transfer orbit:
Wherein, REFor earth mean radiuss.
(2) launch trajectory design
Transmitter section considers disome kinetic model, using Huo Man transfer solving.
The track LEO parameter that known carrier rocket transmitting is entered the orbit, parking orbit orbit altitude is hpark, elliptic orbit 1 far Point height is hepo, calculate the eccentric ratio e of elliptic orbit 1t
Calculate the first subpulse Δ V of Huo Man transfer1
Wherein, μEFor earth gravitational field constant, Δ i is that target track is poor with the inclination angle of parking orbit, WΔiFor weight coefficient, have WΔi∈ [0,1].
Calculate the eccentricity of target track:
Calculate the second subpulse Δ V applying at elliptic orbit 1 apogee2
The gross mass that carrier rocket launches into LEO track is:
Wherein, mparkFor launching the spacecraft gross mass into target track, Δ Vtotal=Δ V1+ΔV2, mupFor rocket Upper Stage Quality, IspUFor the specific impulse of upper stage rocket engine, mtargetIt is the m calculating in step 20.
(3) with J1=mparkFor optimizing index, calculate optimized parameter (tT, θ, WΔi).
Using grid search method to J1It is optimized, each parameter area is respectively:
Step-size in search is taken to be respectively:
Obtain optimized parameter (tT *, θ*, WΔi *) after, according to (1), (2) calculate the track of transmitter section and transfer leg.
In sum, using the present invention, simply it is easy to solve, result of calculation has very big finger for Project Realization to solution principle Lead meaning, and can be applicable to the polar region hovering satellite task design phase of the demands such as polar region observation, communication.
Do not specify in the present invention and partly belong to techniques known.

Claims (10)

1. a kind of design method of the satellite orbit that hovered based on polar region, is carried out by the working track of the satellite that hovered in polar region Design and specificity analysises to complete the system design of whole star and transfer orbit designs it is characterised in that comprising the following steps:
Step one, the dynamical motion of the satellite that hovered in described polar region is described, and the control of the satellite that hovered in described polar region Characteristic is calculated, thus being designed to its working track optimizing to obtain optimal design parameter;
Step 2, according to described optimal design parameter, counts to the attribute of each subsystem as described polar region hovering satellite Calculate, thus obtaining the biggest quality of described polar region hovering Satellite Payloads;And
Step 3, determines the transfer leg track of transfer orbit and the launch trajectory of described polar region hovering satellite respectively, and to institute The orbit parameter obtaining is optimized, thus designing the transfer orbit of optimum,
Wherein, described attribute at least includes quality.
2. according to claim 1 based on polar region hover satellite orbit design method it is characterised in that in institute State execution in step one:
Using circular Restricted three-body model, the dynamical motion of the satellite that hovered in described polar region is described;
Under the conditions of fixing hovering height, the control characteristic of the satellite that hovered in described polar region is analyzed;And
By the way of optimal control problem, the working track of the satellite that hovered in described polar region is designed optimizing.
3. according to claim 2 based on polar region hover satellite orbit design method it is characterised in that described Step one also includes:
Calculate described optimal control problem using optimum control numerical optimization software kit, the control spy of the satellite that hovers with described polar region Property analysis result as optimized algorithm initial value, thus obtaining the working track of optimal-fuel and controlling acceleration and control The time history of power processed.
4. according to claim 3 based on polar region hover satellite orbit design method it is characterised in that
The hover dynamical motion equation of satellite of described polar region is: x . = r . v . m . = v - ▿ U - 2 ω × v + a T - mg T / I sp g 0 , Wherein, described Dynamical motion defines in the rotated coordinate system, and zero refers to earthwards in Solar-terrestrial system barycenter, x-axis for the sun, Z-axis is Solar-terrestrial system rotational angular velocity vector direction, and y-axis constitutes right hand rectangular coordinate system, and r is under described rotating coordinate system The position vector of described base hovering satellite, v is the velocity of described base hovering satellite, and w is described rotating coordinate system Angular velocity vector, aTFor controlling acceleration and its mould aTRepresent, m is the quality of described base hovering satellite, Isp For electric propulsion engine specific impulse, g0For sea level acceleration of gravity, U is potential-energy function and form is U = - μ 1 | | r - r 1 | | - μ 2 | | r - r 2 | | - ω 2 r x 2 + r y 2 2 , μ1、μ2It is respectively the gravitational constant of the sun and the earth, r1、r2Point Not Wei the sun and the earth position vector.
5. according to claim 4 based on polar region hover satellite orbit design method it is characterised in that described Optimal control problem at least includes:Object function, boundary condition and path constraints,
Wherein,
Described object function is:J=-mf=-m (tf),
Described boundary condition is:ry(t0)=0, m (t0)=m0, r (tf)=r (t0), v (tf)=v, (t0), and
Described path constraints are: tan - 1 ( - r 2 , y , r 2 , x ) - wt = 0 r 2 , x 2 + r 2 , y 2 - r 2 , z tan ( δ eq ) = 0 ,
J represents target function value, and polar region hovering track has the periodicity of 1 year, t0For initial time and elect first point of Capricornus as, r(t0) it is to hover the initial time position vector of satellite in described polar region, v (t0) it is initial time velocity, m (t0) it is described Hovering satellite in polar region is in the quality of initial time, tfFor task service life, r (tf), v, (tf) it is to hover the end of satellite in described polar region The position and speed vector in moment, m (tf) be the last moment quality, and (r2, x, r2, y, r2, z)T=r-r2, δeq=23.5 ° is Huang The red angle of cut.
6. according to claim 5 based on polar region hover satellite orbit design method it is characterised in that in institute State execution in step 2:
Mechanical characteristic according to described working track come design described polar region hover satellite subsystem attribute;
According to the gross mass of described polar region hovering satellite, calculate the biggest quality of described payload system;And
When described payload system requires, the gross mass of the satellite that hovered in described polar region is iterated,
Wherein, described attribute at least includes power, size, quality.
7. according to claim 6 based on polar region hover satellite orbit design method it is characterised in that described Step 2 includes:
According to described controling power time history, find the maximum of controling power;
According to described Optimal State variable, obtain the mass change of described polar region hovering satellite;
Calculate fuel, fuel tank system, electric propulsion system, the quality of power-supply system respectively, thus obtaining described effective load The quality of G system;And
Depending on task situation, recalculate the initial mass of described polar region hovering satellite by iteration.
8. according to claim 7 based on polar region hover satellite orbit design method it is characterised in that described The described transfer orbit of polar region hovering satellite is divided into transmitter section and transfer leg,
Wherein, described transmitter section is completed by the mode of rocket Upper Stage chemical propulsion, and described transfer leg is not using applying arteries and veins The trajectory branch mode of punching or continuous thrust completes.
9. according to claim 8 based on polar region hover satellite orbit design method it is characterised in that in institute State execution in step 3:
Using manifold design, the track of described polar region hovering satellite is inversely integrated, thus finding out described transfer leg track, And the orbit parameter of described transfer leg track is changed between described circle Restricted three-body model and two body Models;
Design the launch trajectory of described polar region hovering satellite using Huo Man branch mode;And
Using grid data service, the orbit parameter being obtained is optimized, thus designing the transfer orbit of optimum.
10. according to claim 9 based on polar region hover satellite orbit design method it is characterised in that In described step 3, the described process finding out described transfer leg track includes:
Set out with a point of the polar region hovering track of any time, turn over predetermined angular with rotating coordinate system, according to described circle The restricted three body dynamics models of shape inversely integrate to described track, thus obtaining the track manifold of described polar region hovering satellite;
Block in the nearest position of the liftoff ball of manifold, obtain the reverse time of integration now and position and speed parameter herein;
Position and speed parameter in rotating coordinate system in described circle Restricted three-body model is obtained described two through Coordinate Conversion Position and speed under inertial coodinate system in body Model, and it is calculated corresponding Keplerian orbit six key element;And
The perigee altitude of transfer orbit and altitude of the apogee according to Keplerian orbit six element factor calculation.
CN201510487075.1A 2015-08-10 2015-08-10 Design method based on polar region hovering satellite orbit Active CN106446313B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510487075.1A CN106446313B (en) 2015-08-10 2015-08-10 Design method based on polar region hovering satellite orbit

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510487075.1A CN106446313B (en) 2015-08-10 2015-08-10 Design method based on polar region hovering satellite orbit

Publications (2)

Publication Number Publication Date
CN106446313A true CN106446313A (en) 2017-02-22
CN106446313B CN106446313B (en) 2019-09-13

Family

ID=58092226

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510487075.1A Active CN106446313B (en) 2015-08-10 2015-08-10 Design method based on polar region hovering satellite orbit

Country Status (1)

Country Link
CN (1) CN106446313B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108959827A (en) * 2018-08-10 2018-12-07 哈尔滨工业大学 The design method of polar region suspension railway based on electronic sail
CN110758775A (en) * 2019-11-21 2020-02-07 北京理工大学 Multi-pulse area hovering method based on asteroid surface observation
CN111332498A (en) * 2020-01-10 2020-06-26 北京理工大学 Asteroid balance point hovering detection constant thrust threshold control method
CN113343442A (en) * 2021-05-21 2021-09-03 中国科学院国家空间科学中心 Method and system for solving fixed-time finite fuel multi-pulse transfer orbit

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103093096A (en) * 2013-01-15 2013-05-08 北京航空航天大学 Determination method and determination device of satellite orbit
CN105659867B (en) * 2010-07-02 2013-11-20 北京理工大学 The visible star of satellite navigation system and the DOP value prediction system of user-based security model task
CN104753124A (en) * 2014-12-11 2015-07-01 北京空间飞行器总体设计部 Simulation system for S4R power supply controller of spacecraft

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105659867B (en) * 2010-07-02 2013-11-20 北京理工大学 The visible star of satellite navigation system and the DOP value prediction system of user-based security model task
CN103093096A (en) * 2013-01-15 2013-05-08 北京航空航天大学 Determination method and determination device of satellite orbit
CN104753124A (en) * 2014-12-11 2015-07-01 北京空间飞行器总体设计部 Simulation system for S4R power supply controller of spacecraft

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
WANG ZHI 等: "Hybrid Re-entry Guidance for Reusable Launch Vehicle", 《"APISAT2014", 2014 ASIA-PACIFIC INTERNATIONAL SYMPOSIUM ON AEROSPACE TECHNOLOGY,APISAT2014》 *
林来兴 等: "卫星对空间目标悬停的轨道动力学与控制方法研究", 《中国空间科学技术》 *
饶殷睿 等: "航天器悬停构型设计与控制方法", 《航空学报》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108959827A (en) * 2018-08-10 2018-12-07 哈尔滨工业大学 The design method of polar region suspension railway based on electronic sail
CN108959827B (en) * 2018-08-10 2022-04-15 哈尔滨工业大学 Design method of polar region suspension track based on electric sail
CN110758775A (en) * 2019-11-21 2020-02-07 北京理工大学 Multi-pulse area hovering method based on asteroid surface observation
CN110758775B (en) * 2019-11-21 2021-01-08 北京理工大学 Multi-pulse area hovering method based on asteroid surface observation
CN111332498A (en) * 2020-01-10 2020-06-26 北京理工大学 Asteroid balance point hovering detection constant thrust threshold control method
CN111332498B (en) * 2020-01-10 2021-10-19 北京理工大学 Asteroid balance point hovering detection constant thrust threshold control method
CN113343442A (en) * 2021-05-21 2021-09-03 中国科学院国家空间科学中心 Method and system for solving fixed-time finite fuel multi-pulse transfer orbit

Also Published As

Publication number Publication date
CN106446313B (en) 2019-09-13

Similar Documents

Publication Publication Date Title
Junkins et al. Exploration of alternative state vector choices for low-thrust trajectory optimization
CN104133932B (en) A kind of satellite overall plan certainty annuity based on multidisciplinary optimization and implementation method
CN110276159B (en) Multi-model fusion-based satellite system multidisciplinary optimization method
CN106446313A (en) Polar hovering satellite orbit-based system design method
Klesh et al. Energy-optimal path planning for solar-powered aircraft in level flight
CN106021784A (en) Full-trajectory optimization design method based on two-layer optimization strategy
CN114967453A (en) Satellite east-west coordination state initial value estimation method based on neural network
CN110244767A (en) A kind of formation formation reconfiguration optimization using finite element method
CN115258196A (en) Low-orbit satellite constellation networking electric propulsion orbit transfer strategy optimization method and system
He et al. Solution set calculation of the Sun-perturbed optimal two-impulse trans-lunar orbits using continuation theory
CN101794336B (en) Gravitation-assist simulation method considering influence sphere of target celestial body
Ranieri et al. Indirect optimization of three-dimensional finite-burning interplanetary transfers including spiral dynamics
Lin et al. Entire flight trajectory design for temporary reconnaissance mission
Conway et al. Spacecraft trajectory optimization using direct transcription and nonlinear programming
Maddock et al. Vehicle and mission design of a future small payload launcher
CN112009727B (en) Optimal low-thrust transfer sectional design method for translation point orbit
Bucci et al. Optimal escape manifolds for cis-lunar halo orbits
Balesdent et al. Design and optimization of glide-back reusable launch vehicle architectures
Qingguo et al. A fast computational method for the landing footprints of space-to-ground vehicles
CN114186378A (en) Orbit-changing control method and device for spacecraft, electronic equipment and storage medium
Zhu et al. The Intelligent Trajectory Optimization of Multistage Rocket with Gauss Pseudo-Spectral Method.
O'Brien Using a robotic helicopter to fuel interest in and augment the human exploration of the planet Mars
Akhtar et al. Simulation-based optimization strategy for liquid fueled multi-stage space launch vehicle
Huang et al. Fuel-Optimal Orbit Transfer Trajectory Optimization Method for Satellites
Hernandez et al. Geosynchronous transfer orbit design: A practical approach

Legal Events

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