CN104317303A - Method for controlling formation maintenance or flying-around withdrawal of spacecraft - Google Patents
Method for controlling formation maintenance or flying-around withdrawal of spacecraft Download PDFInfo
- Publication number
- CN104317303A CN104317303A CN201410532652.XA CN201410532652A CN104317303A CN 104317303 A CN104317303 A CN 104317303A CN 201410532652 A CN201410532652 A CN 201410532652A CN 104317303 A CN104317303 A CN 104317303A
- Authority
- CN
- China
- Prior art keywords
- space vehicle
- passive space
- delta
- axis
- spacecraft
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
The invention discloses a method for controlling formation maintenance or flying-around withdrawal of spacecraft. In two adjacent orbit periods of the target spacecraft, integration is performed on the x-coordinate component of the coordinate system of the main control spacecraft relative to the target spacecraft, a difference quotient of the integration result relative to the square of each corresponding period is obtained, and therefore, a period average drift velocity can be obtained; an expected period average drift velocity after the completion of formation maintenance control or flying-around withdrawal is set according to requirements, the period average drift velocity before control is subtracted from the expected period average drift velocity to obtain a period average drift velocity increment, and the horizontal components of double pulses of the formation maintenance control and single pulse of the flying-around withdrawal control both are obtained by virtue of calculation based on the increment; as a result, the formation maintenance control has the characteristics of independence, low frequentness and the like, and the lying-around withdrawal is independent, and safe and quick.
Description
Technical field
The invention belongs to Spacecraft Formation Flying guidance, control technology field, be specifically related to the maintenance of a kind of Spacecraft formation or the control method of withdrawing of being diversion.
Background technology
In Spacecraft Formation Flying field, generally exist with flying, leading relative motion forms such as flying and be diversion between two spacecrafts.Master control spacecraft surrounding target spacecraft does oval relative motion and is called and is diversion, and master control spacecraft is followed and done relative motion after passive space vehicle, is called with flying, and is called that neck flies above.And master control spacecraft enters with flying or lead the state of flying from the state of being diversion for passive space vehicle, being then called is diversion withdraws.The motion of the master control spacecraft relative target spacecraft orbits controlling maintained in the specific interval range of passive space vehicle orbital coordinate system x-axis is called, and formation maintains control.
In formation flight engineering reality, master control spacecraft always drifts about along near-circular orbit passive space vehicle orbital coordinate system+x or-x direction relative to the relative motion configuration of passive space vehicle.In order to reduce as far as possible with the relative motion state that flies, leads to fly or be diversion maintenance frequency, guarantee to maintain safe spacing between two spacecrafts, accurate description relative motion complete cycle configuration has great engineering practical value in the average drift movement velocity in passive space vehicle orbital coordinate system x-axis direction.The free movement analytic solution of existing CW equation comprise a CW drift velocity parameter, but due to CW equation be in formation flight process master control spacecraft relative to the approximate description of the relative motion of passive space vehicle, the relative practical situation error of CW drift velocity is comparatively large, is not enough to the drift velocity of accurate description and forecast relative motion complete cycle configuration.In formation flight field, a relative drift rate parameter is there is in the Relative Orbit Elements theory that development in recent years is got up, make difference by master control spacecraft mean orbit angular velocity relative to passive space vehicle mean orbit angular velocity to obtain, identical with the conceptual nature being called as mean longitude degree drift rate used in the retentive control of high rail satellite position, it is the accurate description of the drift motion trend to relative motion complete cycle configuration, just immeasurable, but because mean orbit angular velocity is difficult to Obtaining Accurate accurately, so this parameter is also inconvenient in Practical Project uses.
Summary of the invention
Technical matters to be solved by this invention is to provide a kind of Spacecraft formation and maintains or the control method of withdrawing of being diversion, and makes formation maintain control and has the features such as autonomous and low-frequency degree, and making to be diversion, it is autonomous and safe and efficient to withdraw.
Technical scheme of the present invention is:
A kind of Spacecraft formation maintains control method, and step is as follows:
(1) when expecting that in follow-up a period of time, orbits controlling all can not occur for master control spacecraft and passive space vehicle, a sampling initial time t carrying out cycle average drift velocity calculating is specified
00;
(2) cycle average drift velocity is carried out as follows
calculating:
X in formula
rt () is the x-axis location components of master control spacecraft in the orbital coordinate system of passive space vehicle corresponding to sampling instant t, obtain according to measurement; T is the passive space vehicle orbital period; t
0for the integration initial time of previous orbital period in adjacent two orbital periods, t
0value is:
(3) judge whether moment t meets t>t
00+ 2T, when meeting t>t
00during+2T, obtain the cycle average drift velocity result that first convergence is available, and proceed to step (4); Otherwise, return step (2);
(4) judge whether moment t arrives to form into columns to maintain and control Startup time t
beforeKeepif, due in t
beforeKeep, obtain this moment t
beforeKeepcorresponding cycle average drift velocity
proceed to step (5); If there is no due in t
beforeKeep, then step (2) is returned;
(5) according to cycle average drift velocity
cycle average drift velocity after having controlled with the maintenance expected
computation period average drift velocity increment
(6) according to cycle average drift velocity increment
calculating forms into columns maintains in level of control dipulse the formation component Δ v realizing the long adjustment of relative major semi-axis
a:
(7) the relative eccentric ratio vector before controlling to implement is maintained according to forming into columns
target relative eccentric ratio vector after having controlled with the maintenance expected
calculate relative eccentric ratio delta vector
wherein
and determine relative eccentric ratio delta vector
at passive space vehicle nodal coordinate system Ox
ny
nz
nmiddle x
naxle and y
nthe coordinate components δ Δ e of axle
xwith δ Δ e
y;
(8) according to formula
calculate the amplitude δ Δ e of relative eccentric ratio delta vector; The formation component Δ v forming into columns and maintain and realize relative eccentric ratio adjustment in level of control dipulse is calculated according to the amplitude δ Δ e of relative eccentric ratio delta vector
e:
In formula
for the flat major semi-axis of passive space vehicle is long, n
targetfor the track angular speed of passive space vehicle;
(9) according to Δ v
awith Δ v
edetermine the horizontal dipulse size maintaining and control of forming into columns, computing formula is:
Wherein, Δ v
x1be the size of the first pulse, Δ v
x2it is the size of the second pulse;
(10) according to δ Δ e
x, δ Δ e
ydetermine the track argument of horizontal dipulse application point on master control spacecraft:
u
x1=atan2(δΔe
y,δΔe
x)
u
x2=atan2(-δΔe
y,-δΔe
x)
Wherein, u
x1for the track argument of described first pulse application point on master control spacecraft, u
x2for the track argument of described second pulse application point on master control spacecraft, atan2 (
★,
★) the expansion codomain that is suitable for for engineering is to the arctan function of [-π, π] scope;
(11) at track argument u
x1place applies the first pulse Δ v to master control spacecraft
x1, at track argument u
x2place applies the second pulse Δ v to master control spacecraft
x2, completing forms into columns maintains control.
The orbital coordinate system of passive space vehicle is defined as: initial point C is positioned at passive space vehicle barycenter, and z-axis points to direction, the earth's core by passive space vehicle barycenter, and x-axis is vertical with z-axis and point to passive space vehicle heading, and y-axis and z-axis and x-axis meet right-hand rule.
When the motion of master control spacecraft relative target spacecraft exceeds the specific interval range of passive space vehicle orbital coordinate system x-axis, think that moment t arrives to form into columns to maintain instantly and control Startup time t
beforeKeep, otherwise, think that moment t not yet arrives to form into columns to maintain instantly and control Startup time t
beforeKeep.
Passive space vehicle nodal coordinate system is defined as: initial point Ο is the earth's core, axle x
npassive space vehicle ascending node of orbit is pointed to, axle z from Ο
nalong passive space vehicle orbital angular momentum vector direction, axle y
nwith z
naxle and x
naxle meets right-hand rule.
A kind of spacecraft is diversion and withdraws control method, and step is as follows:
(1) when expecting that in follow-up a period of time, orbits controlling all can not occur for master control spacecraft and passive space vehicle, a sampling initial time t carrying out cycle average drift velocity calculating is specified
00;
(2) cycle average drift velocity is carried out as follows
calculating:
X in formula
rt () is the x-axis location components of master control spacecraft in the orbital coordinate system of passive space vehicle corresponding to sampling instant t, obtain according to measurement; T is the passive space vehicle orbital period; t
0for the integration initial time of previous orbital period in adjacent two orbital periods, t
0value is:
(3) judge whether moment t meets t>t
00+ 2T, when meeting t>t
00during+2T, obtain the cycle average drift velocity result that first convergence is available, and proceed to step (4); Otherwise, return step (2);
(4) judge whether moment t arrives the departure time t that is diversion specified instantly
beforeWithdrawif, due in t
before Withdraw, obtain this moment t
before Withdrawcorresponding cycle average drift velocity
proceed to step (5); If there is no due in t
before?
withdraw, then step (2) is returned;
(5) basis
with the cycle average drift velocity after having withdrawn of being diversion expected
computation period average drift velocity increment
(6) according to cycle average drift velocity increment
calculating is diversion and is withdrawn pulse
at the component Δ v of passive space vehicle orbital coordinate system x-axis
x:
Setting is diversion and is withdrawn pulse
at the component Δ v of passive space vehicle orbital coordinate system z-axis
zfor:
In formula
before withdrawing for being diversion, master control spacecraft is at the z-axis relative velocity component relative to passive space vehicle withdrawing pulse action point place that is diversion;
According to Δ v
xwith Δ v
zobtain being diversion and withdraw pulse
size delta v
xzfor:
According to Δ v
xwith Δ v
zcalculate this to be diversion and to withdraw pulse
relative to the pitching angle theta of passive space vehicle orbital coordinate system+x-axis
withdraw, wherein θ
withdraw=atan2 (-Δ v
z, Δ v
x); And specify this to be diversion to withdraw pulse
it is 0 ° relative to the crab angle of passive space vehicle orbital coordinate system+x-axis;
(7) withdraw pulse action point place and apply to be diversion and to withdraw pulse being diversion
completing is diversion withdraws control.
The orbital coordinate system of passive space vehicle is defined as: initial point C is positioned at passive space vehicle barycenter, and z-axis points to direction, the earth's core by passive space vehicle barycenter, and x-axis is vertical with z-axis and point to passive space vehicle heading, and y-axis and z-axis and x-axis meet right-hand rule.
When preparing to withdraw from passive space vehicle front, being diversion and withdrawing pulse action point selection at the oval some place crossing with passive space vehicle orbital coordinate system+x-axis of being diversion; When preparing to withdraw from passive space vehicle rear, being diversion and withdrawing pulse action point selection at the oval some place crossing with passive space vehicle orbital coordinate system-x-axis of being diversion.
The present invention's advantage is compared with prior art: for obtaining the situation of master control spacecraft relative to passive space vehicle station-keeping data continuously, directly at near-circular orbit passive space vehicle x coordinate direction, get the relative measurement data of two complete cycles, by after cycle summation again to summed result about passive space vehicle orbital period square do difference coefficient, provide the parameter that is referred to as cycle average drift velocity.It is firm that this parameter directly has speed amount, describes the average drift movement velocity of relative motion complete cycle configuration exactly.Set after formation retentive control completes as required or be diversion withdraw after the cycle average drift velocity expected, before this expectation drift velocity deducts control, cycle average drift velocity obtains cycle average drift velocity increment, dipulse and the horizontal component of withdrawing control monopulse of being diversion of formation retentive control obtain by this incremental computations, thus make formation retentive control have the features such as autonomous and low-frequency degree, making to be diversion, it is autonomous and safe and efficient to withdraw.
Accompanying drawing explanation
Fig. 1 is the FB(flow block) of the inventive method;
Fig. 2 is passive space vehicle orbital coordinate system definition schematic diagram;
Fig. 3 is passive space vehicle nodal coordinate system definition schematic diagram;
Fig. 4 is that master control Spacecraft formation maintains control dipulse application point schematic diagram;
Fig. 5 is be diversion the pulse action point and pulse orientation schematic diagram withdrawn.
Embodiment
As shown in Figure 1, Spacecraft formation of the present invention maintains control method, and step is as follows:
(1) when expecting that in follow-up a period of time, orbits controlling all can not occur for master control spacecraft and passive space vehicle, a sampling initial time t carrying out cycle average drift velocity calculating is specified
00.
(2) cycle average drift velocity is carried out as follows
calculating:
X in formula
rt () is the x-axis location components of master control spacecraft in the orbital coordinate system of passive space vehicle corresponding to sampling instant t, obtain according to measurement; As shown in Figure 2, the orbital coordinate system of passive space vehicle is defined as: initial point C is positioned at passive space vehicle barycenter, z-axis points to direction, the earth's core by passive space vehicle barycenter, and x-axis is vertical with z-axis and point to passive space vehicle heading, and y-axis and z-axis and x-axis meet right-hand rule; T is the passive space vehicle orbital period; t
0for the integration initial time of previous orbital period in adjacent two orbital periods, t
0value is:
I.e. t
0at t≤t
00when+2T, equal value is t
00, but t
0from t>t
00first point of+2T starts to slide on a timeline value, thus the slip on a timeline of average drift velocity performance period calculates.
(3) judge whether moment t meets t>t
00+ 2T, when meeting t>t
00during+2T, obtain the cycle average drift velocity result that first convergence is available, and proceed to step (4); Otherwise, return step (2).
(4) judge whether moment t arrives to form into columns to maintain instantly and control Startup time t
beforeKeepif, due in t
beforeKeep, obtain this moment t
beforeKeepcorresponding cycle average drift velocity
proceed to step (5); If there is no due in t
beforeKeep, then step (2) is returned;
When the motion of master control spacecraft relative target spacecraft exceeds the specific interval range of passive space vehicle orbital coordinate system x-axis, think that moment t arrives to form into columns to maintain instantly and control Startup time t
beforeKeep, otherwise, think that moment t not yet arrives to form into columns to maintain instantly and control Startup time t
beforeKeep.
(5) according to cycle average drift velocity
cycle average drift velocity after having controlled with the maintenance expected
calculate the cycle average drift velocity increment of forming into columns and maintaining and controlling
Reversal periods average drift velocity after the retentive control expected completes is
according to form into columns with fly or be diversion configuration along passive space vehicle orbital coordinate system x-axis the circle time of specifying wish drift about distance determine, be known quantity.
(6) according to cycle average drift velocity increment
calculating forms into columns maintains in level of control dipulse the formation component Δ v realizing the long adjustment of relative major semi-axis
a:
(7) the relative eccentric ratio vector before controlling to implement is maintained according to forming into columns
target relative eccentric ratio vector after having controlled with the maintenance expected
calculate relative eccentric ratio delta vector
And determine relative eccentric ratio delta vector
at passive space vehicle nodal coordinate system Ox
ny
nz
nmiddle x
naxle and y
nthe coordinate components δ Δ e of axle
xwith δ Δ e
y, namely at passive space vehicle nodal coordinate system Ox
ny
nz
nin have:
As shown in Figure 3, passive space vehicle nodal coordinate system Ox
ny
nz
nbe defined as: initial point Ο is the earth's core, axle x
npassive space vehicle ascending node of orbit is pointed to, axle z from Ο
nalong passive space vehicle orbital angular momentum vector direction, axle y
nwith z
naxle and x
naxle meets right-hand rule.
(8) according to formula
calculate the amplitude δ Δ e of relative eccentric ratio delta vector; The formation component Δ v forming into columns and maintain and realize relative eccentric ratio adjustment in level of control dipulse is calculated according to the amplitude δ Δ e of relative eccentric ratio delta vector
e:
In formula
for the flat major semi-axis of passive space vehicle is long, n
targetfor the track angular speed of passive space vehicle.
(9) according to Δ v
awith Δ v
edetermine the horizontal dipulse size maintaining and control of forming into columns, computing formula is:
Wherein, Δ v
x1be the size of the first pulse, Δ v
x2it is the size of the second pulse;
As Δ v
x1or Δ v
x2when being greater than 0, corresponding pulse action direction is consistent with passive space vehicle orbital coordinate system+x direction; As Δ v
x1or Δ v
x2when being less than 0, corresponding pulse action direction is contrary with passive space vehicle orbital coordinate system-x direction.
(10) according to δ Δ e
x, δ Δ e
ydetermine the track argument of horizontal dipulse application point on master control spacecraft:
u
x1=atan2(δΔe
y,δΔe
x)
u
x2=atan2(-δΔe
y,-δΔe
x)
Wherein, u
x1for the track argument of described first pulse application point on master control spacecraft, u
x2for the track argument of described second pulse application point on master control spacecraft, see Fig. 4; Atan2 (
★,
★) the expansion codomain that is suitable for for engineering is to the arctan function of [-π, π] scope;
(11) at track argument u
x1place applies the first pulse Δ v to master control spacecraft
x1, at track argument u
x2place applies the second pulse Δ v to master control spacecraft
x2, completing forms into columns maintains control.
As shown in Figure 1, spacecraft of the present invention is diversion and withdraws control method, and step is as follows:
(1) when expecting that in follow-up a period of time, orbits controlling all can not occur for master control spacecraft and passive space vehicle, a sampling initial time t carrying out cycle average drift velocity calculating is specified
00;
(2) cycle average drift velocity is carried out as follows
calculating:
X in formula
rt () is the x-axis location components of master control spacecraft in the orbital coordinate system of passive space vehicle corresponding to sampling instant t, obtain according to measurement; As shown in Figure 2, the orbital coordinate system of passive space vehicle is defined as: initial point C is positioned at passive space vehicle barycenter, z-axis points to direction, the earth's core by passive space vehicle barycenter, and x-axis is vertical with z-axis and point to passive space vehicle heading, and y-axis and z-axis and x-axis meet right-hand rule; T is the passive space vehicle orbital period; t
0for the integration initial time of previous orbital period in adjacent two orbital periods, t
0value is:
I.e. t
0at t≤t
00when+2T, equal value is t
00, but t
0from t>t
00first point of+2T starts to slide on a timeline value, thus the slip on a timeline of average drift velocity performance period calculates;
(3) judge whether moment t meets t>t
00+ 2T, when meeting t>t
00during+2T, obtain the cycle average drift velocity result that first convergence is available, and proceed to step (4); Otherwise, return step (2);
(4) judge whether moment t arrives the departure time t that is diversion instantly
beforeWithdrawif, due in t
beforeWithdraw, obtain this moment t
beforeWithdrawcorresponding cycle average drift velocity
proceed to step (5); If there is no due in t
beforeWithdrawthen return step (2);
Be diversion departure time t
beforeWithdrawfor the known quantity of specifying;
(5) basis
with the cycle average drift velocity after having withdrawn of being diversion expected
calculating is diversion withdraws the cycle average drift velocity increment of control
Wherein, the cycle average drift velocity after having withdrawn of being diversion expected
for the known quantity of specifying;
(6) according to cycle average drift velocity increment
calculating is diversion and is withdrawn pulse
at the component Δ v of passive space vehicle orbital coordinate system x-axis
x:
Setting is diversion and is withdrawn pulse
at the component Δ v of passive space vehicle orbital coordinate system z-axis
zfor:
In formula
before withdrawing for being diversion, master control spacecraft is withdrawing the z-axis relative velocity component relative to passive space vehicle at pulse action point place; As shown in Figure 5, when preparing to withdraw from passive space vehicle front, being diversion and withdrawing pulse action point selection at the oval some place crossing with passive space vehicle orbital coordinate system+x-axis of being diversion; When preparing to withdraw from passive space vehicle rear, being diversion and withdrawing pulse action point selection at the oval some place crossing with passive space vehicle orbital coordinate system-x-axis of being diversion;
According to Δ v
xwith Δ v
zobtain being diversion and withdraw pulse
size delta v
xzfor:
See Fig. 5, according to Δ v
xwith Δ v
zcalculate this to be diversion and to withdraw pulse
relative to the pitching angle theta of passive space vehicle orbital coordinate system+x-axis
withdraw:
θ
Withdraw=atan2(-Δv
z,Δv
x)
Atan2 (
★,
★) the expansion codomain that is suitable for for engineering is to the arctan function of [-π, π] scope; And specify this to be diversion to withdraw pulse
it is 0 ° relative to the crab angle of passive space vehicle orbital coordinate system+x-axis;
(7) apply to be diversion and to withdraw pulse withdrawing pulse action point place
completing is diversion withdraws control.
The content be not described in detail in instructions of the present invention belongs to the known technology of those skilled in the art.
Claims (7)
1. Spacecraft formation maintains a control method, it is characterized in that step is as follows:
(1) when expecting that in follow-up a period of time, orbits controlling all can not occur for master control spacecraft and passive space vehicle, a sampling initial time t carrying out cycle average drift velocity calculating is specified
00;
(2) cycle average drift velocity is carried out as follows
calculating:
X in formula
rt () is the x-axis location components of master control spacecraft in the orbital coordinate system of passive space vehicle corresponding to sampling instant t, obtain according to measurement; T is the passive space vehicle orbital period; t
0for the integration initial time of previous orbital period in adjacent two orbital periods, t
0value is:
(3) judge whether moment t meets t>t
00+ 2T, when meeting t>t
00during+2T, obtain the cycle average drift velocity result that first convergence is available, and proceed to step (4); Otherwise, return step (2);
(4) judge whether moment t arrives to form into columns to maintain and control Startup time t
beforeKeepif, due in t
beforeKeep, obtain this moment t
beforeKeepcorresponding cycle average drift velocity
proceed to step (5); If there is no due in t
beforeKeep, then step (2) is returned;
(5) according to cycle average drift velocity
cycle average drift velocity after having controlled with the maintenance expected
computation period average drift velocity increment
(6) according to cycle average drift velocity increment
calculating forms into columns maintains in level of control dipulse the formation component Δ v realizing the long adjustment of relative major semi-axis
a:
(7) the relative eccentric ratio vector before controlling to implement is maintained according to forming into columns
target relative eccentric ratio vector after having controlled with the maintenance expected
calculate relative eccentric ratio delta vector
wherein
and determine relative eccentric ratio delta vector
at passive space vehicle nodal coordinate system Ox
ny
nz
nmiddle x
naxle and y
nthe coordinate components δ Δ e of axle
xwith δ Δ e
y;
(8) according to formula
calculate the amplitude δ Δ e of relative eccentric ratio delta vector; The formation component Δ v forming into columns and maintain and realize relative eccentric ratio adjustment in level of control dipulse is calculated according to the amplitude δ Δ e of relative eccentric ratio delta vector
e:
In formula
for the flat major semi-axis of passive space vehicle is long, n
targetfor the track angular speed of passive space vehicle;
(9) according to Δ v
awith Δ v
edetermine the horizontal dipulse size maintaining and control of forming into columns, computing formula is:
Wherein, Δ v
x1be the size of the first pulse, Δ v
x2it is the size of the second pulse;
(10) according to δ Δ e
x, δ Δ e
ydetermine the track argument of horizontal dipulse application point on master control spacecraft:
u
x1=atan2(δΔe
y,δΔe
x)
u
x2=atan2(-δΔe
y,-δΔe
x)
Wherein, u
x1for the track argument of described first pulse application point on master control spacecraft, u
x2for the track argument of described second pulse application point on master control spacecraft, the expansion codomain that atan2 (*, *) is suitable for for engineering is to the arctan function of [-π, π] scope;
(11) at track argument u
x1place applies the first pulse Δ v to master control spacecraft
x1, at track argument u
x2place applies the second pulse Δ v to master control spacecraft
x2, completing forms into columns maintains control.
2. Spacecraft formation according to claim 1 maintains control method, it is characterized in that, the orbital coordinate system of passive space vehicle is defined as: initial point C is positioned at passive space vehicle barycenter, z-axis points to direction, the earth's core by passive space vehicle barycenter, x-axis is vertical with z-axis and point to passive space vehicle heading, and y-axis and z-axis and x-axis meet right-hand rule.
3. Spacecraft formation according to claim 1 maintains control method, it is characterized in that, when the motion of master control spacecraft relative target spacecraft exceeds the specific interval range of passive space vehicle orbital coordinate system x-axis, think that moment t arrives to form into columns to maintain instantly and control Startup time t
beforeKeep, otherwise, think that moment t not yet arrives to form into columns to maintain instantly and control Startup time t
beforeKeep.
4. a kind of Spacecraft formation according to claim 1 maintains control method, and it is characterized in that, passive space vehicle nodal coordinate system is defined as: initial point Ο is the earth's core, axle x
npassive space vehicle ascending node of orbit is pointed to, axle z from Ο
nalong passive space vehicle orbital angular momentum vector direction, axle y
nwith z
naxle and x
naxle meets right-hand rule.
5. spacecraft is diversion and withdraws a control method, it is characterized in that step is as follows:
(1) when expecting that in follow-up a period of time, orbits controlling all can not occur for master control spacecraft and passive space vehicle, a sampling initial time t carrying out cycle average drift velocity calculating is specified
00;
(2) cycle average drift velocity is carried out as follows
calculating:
X in formula
rt () is the x-axis location components of master control spacecraft in the orbital coordinate system of passive space vehicle corresponding to sampling instant t, obtain according to measurement; T is the passive space vehicle orbital period; t
0for the integration initial time of previous orbital period in adjacent two orbital periods, t
0value is:
(3) judge whether moment t meets t>t
00+ 2T, when meeting t>t
00during+2T, obtain the cycle average drift velocity result that first convergence is available, and proceed to step (4); Otherwise, return step (2);
(4) judge whether moment t arrives the departure time t that is diversion specified instantly
beforeWithdrawif, due in t
beforeWithdraw, obtain this moment t
beforeWithdrawcorresponding cycle average drift velocity
proceed to step (5); If there is no due in t
beforeWithdraw, then step (2) is returned;
(5) basis
with the cycle average drift velocity after having withdrawn of being diversion expected
computation period average drift velocity increment
(6) according to cycle average drift velocity increment
calculating is diversion and is withdrawn pulse
at the component Δ v of passive space vehicle orbital coordinate system x-axis
x:
Setting is diversion and is withdrawn pulse
at the component Δ v of passive space vehicle orbital coordinate system z-axis
zfor:
In formula
before withdrawing for being diversion, master control spacecraft is at the z-axis relative velocity component relative to passive space vehicle withdrawing pulse action point place that is diversion;
According to Δ v
xwith Δ v
zobtain being diversion and withdraw pulse
size delta v
xzfor:
According to Δ v
xwith Δ v
zcalculate this to be diversion and to withdraw pulse
relative to the pitching angle theta of passive space vehicle orbital coordinate system+x-axis
withdraw, wherein θ
withdraw=atan2 (-Δ v
z, Δ v
x); And specify this to be diversion to withdraw pulse
it is 0 ° relative to the crab angle of passive space vehicle orbital coordinate system+x-axis;
(7) withdraw pulse action point place and apply to be diversion and to withdraw pulse being diversion
completing is diversion withdraws control.
6. spacecraft according to claim 5 is diversion and withdraws control method, it is characterized in that, the orbital coordinate system of passive space vehicle is defined as: initial point C is positioned at passive space vehicle barycenter, z-axis points to direction, the earth's core by passive space vehicle barycenter, x-axis is vertical with z-axis and point to passive space vehicle heading, and y-axis and z-axis and x-axis meet right-hand rule.
7. spacecraft according to claim 5 is diversion and withdraws control method, it is characterized in that, when preparing to withdraw from passive space vehicle front, being diversion and withdrawing pulse action point selection at the oval some place crossing with passive space vehicle orbital coordinate system+x-axis of being diversion; When preparing to withdraw from passive space vehicle rear, being diversion and withdrawing pulse action point selection at the oval some place crossing with passive space vehicle orbital coordinate system-x-axis of being diversion.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410532652.XA CN104317303B (en) | 2014-10-10 | 2014-10-10 | A kind of Spacecraft formation maintains or the control method withdrawn of being diversion |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410532652.XA CN104317303B (en) | 2014-10-10 | 2014-10-10 | A kind of Spacecraft formation maintains or the control method withdrawn of being diversion |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104317303A true CN104317303A (en) | 2015-01-28 |
CN104317303B CN104317303B (en) | 2016-09-21 |
Family
ID=52372545
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410532652.XA Active CN104317303B (en) | 2014-10-10 | 2014-10-10 | A kind of Spacecraft formation maintains or the control method withdrawn of being diversion |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104317303B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104765374A (en) * | 2015-02-10 | 2015-07-08 | 北京控制工程研究所 | High-orbit natural-flying-around-track correcting method |
CN114063645A (en) * | 2021-11-12 | 2022-02-18 | 中国西安卫星测控中心 | Eccentricity dip angle vector-based method for evaluating tilt fly-around holding control effect |
CN114872935A (en) * | 2022-06-09 | 2022-08-09 | 北京航天驭星科技有限公司 | Orbital transfer method for controlling approach and flying around of two spacecrafts on near-circular orbit |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6341249B1 (en) * | 1999-02-11 | 2002-01-22 | Guang Qian Xing | Autonomous unified on-board orbit and attitude control system for satellites |
CN101520511B (en) * | 2009-03-13 | 2011-05-04 | 北京航空航天大学 | Method for formation configuration of distributed satellites with synthetic aperture radars |
JP6210749B2 (en) * | 2013-06-18 | 2017-10-11 | 三菱電機株式会社 | Spacecraft orbit control method and apparatus |
CN103345256B (en) * | 2013-06-28 | 2016-04-13 | 上海新跃仪表厂 | Based on relative accompanying flying and the position transfer control method of Relative Orbit Elements |
-
2014
- 2014-10-10 CN CN201410532652.XA patent/CN104317303B/en active Active
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104765374A (en) * | 2015-02-10 | 2015-07-08 | 北京控制工程研究所 | High-orbit natural-flying-around-track correcting method |
CN104765374B (en) * | 2015-02-10 | 2017-05-10 | 北京控制工程研究所 | High-orbit natural-flying-around-track correcting method |
CN114063645A (en) * | 2021-11-12 | 2022-02-18 | 中国西安卫星测控中心 | Eccentricity dip angle vector-based method for evaluating tilt fly-around holding control effect |
CN114872935A (en) * | 2022-06-09 | 2022-08-09 | 北京航天驭星科技有限公司 | Orbital transfer method for controlling approach and flying around of two spacecrafts on near-circular orbit |
Also Published As
Publication number | Publication date |
---|---|
CN104317303B (en) | 2016-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104142686B (en) | A kind of satellite Autonomous formation flight control method | |
CN103728980B (en) | The control method of spacecraft relative orbit | |
CN103941271B (en) | A kind of GPS/SINS hypercompact combination navigation method of time-space difference | |
CN104076819B (en) | The boundary control method of satellite bounded accompanying flying under a kind of round reference orbit | |
CN102878997B (en) | Quick high accuracy Extrapolation method on a kind of star of highly eccentric orbit | |
CN105182985A (en) | Hypersonic flight vehicle dive segment full amount integration guidance control method | |
CN103268067B (en) | A kind of satellite based on intending hypercomplex number and plan hypercomplex number kinematical equation points to tracking and controlling method | |
CN103438890B (en) | Based on the planetary power descending branch air navigation aid of TDS and image measurement | |
CN105116914B (en) | A kind of stratospheric airship analytic modell analytical model predicted path tracking and controlling method | |
CN110632935B (en) | Autonomous control method for formation satellite flying around | |
CN102508280B (en) | Method for assisting double-antenna measuring unit in determining integer ambiguity and heading by gyroscope | |
CN104443432A (en) | Autonomous track transfer guidance method for satellite limit thrust coplane circular orbit | |
CN104317303A (en) | Method for controlling formation maintenance or flying-around withdrawal of spacecraft | |
CN102944241A (en) | Spacecraft relative attitude determining method based on multicell liner differential inclusion | |
CN105157724A (en) | Transfer alignment time delay estimation and compensation method based on velocity plus attitude matching | |
CN105241319B (en) | A kind of guided cartridge of spin at a high speed real-time alignment methods in the air | |
CN103033197B (en) | Method for correcting zero drift of MEMS gyroscope | |
CN104076818A (en) | Space rendezvous system gain scheduling control method with linearization errors taken into consideration | |
CN101038170B (en) | Method for online estimating inertia/satellite combined guidance system data asynchronous time | |
CN103950555A (en) | High-precision keeping and control method for relative positions with ultra-close distance | |
CN103245948B (en) | Image match navigation method for double-area image formation synthetic aperture radars | |
CN102997932A (en) | Method for eliminating shaking impact of rotation table during high precision inertial navigation system calibration | |
Wang et al. | Novel in-flight coarse alignment of low-cost strapdown inertial navigation system for unmanned aerial vehicle applications | |
CN103900561A (en) | Middle-distance and long-distance approximate position estimation method of space tethered robot | |
CN106054227A (en) | Pseudo range difference single-star high-dynamic positioning method in inertial navigation assistance |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |