CN113703487B - Small satellite formation configuration control method based on single electric push - Google Patents

Small satellite formation configuration control method based on single electric push Download PDF

Info

Publication number
CN113703487B
CN113703487B CN202111104268.6A CN202111104268A CN113703487B CN 113703487 B CN113703487 B CN 113703487B CN 202111104268 A CN202111104268 A CN 202111104268A CN 113703487 B CN113703487 B CN 113703487B
Authority
CN
China
Prior art keywords
formation configuration
plane
pulse
thrust
application
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.)
Active
Application number
CN202111104268.6A
Other languages
Chinese (zh)
Other versions
CN113703487A (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.)
Sun Yat Sen University
Original Assignee
Sun Yat Sen University
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 Sun Yat Sen University filed Critical Sun Yat Sen University
Priority to CN202111104268.6A priority Critical patent/CN113703487B/en
Publication of CN113703487A publication Critical patent/CN113703487A/en
Application granted granted Critical
Publication of CN113703487B publication Critical patent/CN113703487B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/10Simultaneous control of position or course in three dimensions
    • G05D1/101Simultaneous control of position or course in three dimensions specially adapted for aircraft
    • G05D1/104Simultaneous control of position or course in three dimensions specially adapted for aircraft involving a plurality of aircrafts, e.g. formation flying
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

The invention discloses a control method for a small satellite formation configuration based on single electric push, which comprises the following steps: calculating formation configuration parameter deviation; calculating the suboptimal three pulse size of the fuel and the corresponding pulse applying time, and the optimal single pulse size of the out-of-plane fuel and the corresponding pulse applying time; calculating the total duration of the application of the continuous low thrust at each pulse; setting continuous low-thrust application duration of each continuous low-thrust application arc section, and reserving gesture maneuvering time; maneuvering the satellite attitude to a target attitude that can produce a thrust in a predetermined direction during an attitude maneuver period; applying continuous low-thrust control at each set continuous low-thrust application arc segment, and maneuvering the satellite to a conventional attitude; and updating the current formation configuration parameters until the deviation of the formation configuration parameters is judged to be within a preset range. By using the invention, the control of the small satellite formation configuration with only a single electric propulsion, which is less in fuel consumption, is realized. The method can be widely applied to the field of satellite formation configuration control.

Description

Small satellite formation configuration control method based on single electric push
Technical Field
The invention relates to the field of satellite formation configuration control, in particular to a small satellite formation configuration control method based on single electric push.
Background
Compared with the traditional single large satellite, the small satellite formation control has the advantages of meeting task requirements, along with small size, low cost, high reliability and strong flexibility. In order to prolong the on-orbit life of the formation satellite, the total fuel consumption required by formation configuration control is usually used as an optimization target, and a pulse-based formation configuration reconstruction control optimization solution is obtained, but when a propulsion system of the formation system is electric propulsion or micro-cooling air propulsion, the suboptimal pulse solution is not applicable any more, and a formation continuous low-thrust configuration control method meeting both control precision requirements and suboptimal fuel consumption needs to be provided for the conditions of continuous low thrust and single electric propulsion.
Disclosure of Invention
In order to solve the technical problems, the invention aims to provide a small satellite formation configuration control method based on single electric propulsion, which realizes formation continuous low-thrust configuration control with suboptimal fuel consumption.
The technical scheme adopted by the invention is as follows: a control method for a small satellite formation configuration based on single electric push comprises the following steps:
s1, calculating formation configuration parameter deviation according to a current formation configuration parameter and a target formation configuration parameter;
s2, judging that the deviation of formation configuration parameters exceeds a preset threshold, and calculating the in-plane suboptimal three-pulse size, the corresponding suboptimal pulse application latitude amplitude angle, the out-of-plane optimal pulse size and the corresponding optimal pulse application latitude amplitude angle by taking less fuel consumption as a target according to the current formation configuration parameters;
s3, calculating the total duration of continuous small thrust application at each pulse according to the speed increment which can be provided in the unit time of the maximum thrust of the single electric propeller;
s4, considering the influence of the continuous small thrust application duration of each track on the formation configuration control time and the fuel consumption, setting the continuous small thrust application duration of each track and the attitude maneuver time which can generate the thrust in the required direction before and after each pulse application moment;
s5, maneuvering the satellite attitude to a target attitude capable of enabling a single electric propulsion to be adjusted to generate a preset control thrust direction in the attitude maneuvering time period;
s6, applying continuous low-thrust control at the set continuous low-thrust application position of each track, and maneuvering the satellite to a conventional posture after the application of the continuous low-thrust control is finished;
and S7, updating the current formation configuration parameters and returning to the step S1 until the deviation of the formation configuration parameters is judged to be within a preset range.
Further, the formation configuration parameter deviation comprises a difference between a formation auxiliary star and a main star semi-long axis, a relative eccentricity vector, a relative inclination vector and a formation around a flying center along a track direction offset distance.
Further, the calculation formula of the in-plane suboptimal three-pulse size is as follows:
Figure BDA0003271625370000021
Figure BDA0003271625370000022
Figure BDA0003271625370000023
in the above-mentioned method, the step of,
Figure BDA0003271625370000024
representing the first pulse size in plane, +.>
Figure BDA0003271625370000025
Representing in-plane second pulse size, < >>
Figure BDA0003271625370000026
Representing the third pulse size in plane, Δa m Indicating the expected suboptimal pulse control of the second half long axis change in ignition, Δλ sm Representing the amount of change in the constant offset along the track, Δe m The magnitude change of the E vector is represented, n represents the angular velocity of the primary star orbit, and T represents the period of the formation primary star orbit.
Further, the calculation formula of the latitude argument of the corresponding suboptimal pulse application is as follows:
Figure BDA0003271625370000027
Figure BDA0003271625370000028
Figure BDA0003271625370000029
in the above-mentioned method, the step of,
Figure BDA00032716253700000210
representing the corresponding latitude depression angle of the application of the first pulse in the plane, < >>
Figure BDA00032716253700000211
Indicating correspondence when applying in-plane second pulsesIs a latitude depression angle,/>
Figure BDA00032716253700000212
Representing the corresponding latitude depression angle of the application of the third pulse in the plane,/>
Figure BDA00032716253700000213
Component representing the desired change in the y-direction of the relative eccentricity vector,/->
Figure BDA00032716253700000214
Representing the component of the desired change in the x-direction of the relative eccentricity vector.
Further, the calculation formula of the out-of-plane optimal pulse size is as follows:
δv n1 =na|Δi target -Δi current |=na|Δi man |
in the above, δv n1 Represents the out-of-plane optimal pulse size, n represents the primary satellite orbit angular velocity, a represents the primary satellite orbit semi-major axis, Δi man Representing the relative tilt angle vector of the desired change.
Further, the calculation formula of the corresponding latitude amplitude of the optimal pulse application is as follows:
Figure BDA0003271625370000031
in the above, u n1 Representing the corresponding optimal pulse application latitude argument,
Figure BDA0003271625370000032
component representing the desired change in the y-direction of the relative tilt vector, is provided>
Figure BDA0003271625370000033
Representing the component of the desired change in the x-direction of the relative tilt vector.
Further, the calculation formula of the total duration of the continuous small thrust application at each pulse is as follows:
Figure BDA0003271625370000034
Figure BDA0003271625370000035
Figure BDA0003271625370000036
Figure BDA0003271625370000037
in the above-mentioned method, the step of,
Figure BDA0003271625370000038
indicates the duration of the force application required to achieve the first pulse size,/->
Figure BDA0003271625370000039
Indicates the duration of the force required to achieve the second pulse size,/->
Figure BDA00032716253700000310
Indicating the duration of the force applied, deltaT, required to achieve the third pulse size n1 Represents the duration of the applied force, deltaV, required to achieve an out-of-plane optimum pulse size max Representing the speed increase that can be provided per unit time of maximum thrust of a single electric propulsion.
Further, the formula for setting the continuous low thrust application duration per track before and after each pulse application time is as follows:
Figure BDA00032716253700000311
Figure BDA00032716253700000312
k∈1,2,3
Figure BDA00032716253700000313
Figure BDA00032716253700000314
in the above-mentioned method, the step of,
Figure BDA00032716253700000315
indicates the start time of each low thrust control in the plane, < > and the like>
Figure BDA00032716253700000316
Indicating the end time of each low thrust control in the plane, < > and the like>
Figure BDA00032716253700000317
Indicates the in-plane force application time length, (T) n1 ) start The out-of-plane single low thrust control start time is represented by (T) n1 ) end And the out-of-plane single small thrust control end time is indicated.
The method and the system have the beneficial effects that: according to the invention, the pulse control suboptimal solution is decomposed into continuous low-thrust formation configuration control near the suboptimal application point, and for the limitation of single electric pushing of a satellite, the direction control of the thrust is realized by regulating the gesture through a satellite flywheel at the suboptimal pulse, so that the aims of simultaneously meeting formation configuration control precision and saving fuel consumption are fulfilled.
Drawings
FIG. 1 is a flow chart of the steps of a control method for a formation configuration of a small satellite based on a single electric push.
Detailed Description
The invention will now be described in further detail with reference to the drawings and to specific examples. The step numbers in the following embodiments are set for convenience of illustration only, and the order between the steps is not limited in any way, and the execution order of the steps in the embodiments may be adaptively adjusted according to the understanding of those skilled in the art.
Referring to fig. 1, the invention provides a control method for a small satellite formation configuration based on single electric push.
In a first embodiment, in-plane formation configuration control:
s101, according to the current in-plane formation configuration parameter delta alpha current =(Δa current ,Δe current ,Δi current ,Δλ current ) T And a target in-plane formation configuration parameter Δα target =(Δa target ,Δe target ,Δi target ,Δλ target ) T Calculating the in-plane formation configuration parameter deviation delta alpha error =(Δa error ,Δe error ,Δi error ,Δλ error ) T Wherein Δa=a deputy -a chief To form the difference between the minor and major semi-major axes of the minor and major stars, Δe= (Δe) x ,Δe y ) T =(e d cosω d -e c cosω c ,e d sinω d -e c sinω c ) T As relative eccentricity vector, Δi= (Δi) x ,Δi y ) T =(i d -i c ,(Ω dc )sini c ) T As relative tilt vector, Δλ=Δu+ΔΩ cosi c Offset distances along the track around the center of flight for formation.
S102, judging the deviation delta alpha of formation configuration parameters in a plane error Exceeding a preset threshold Δα predefined According to the current in-plane formation parameter Deltaalpha current And calculating the suboptimal three-pulse size and the corresponding suboptimal pulse application latitude amplitude angle in the plane.
Specifically, based on the characteristic that in-plane and out-of-plane movements of the formation can be decoupled, the in-plane configuration parameter delta alpha is calculated for the formation in-plane =(Δa,Δe,Δλ) T Control target, calculating the magnitude of sub-optimal along-track three-pulse control quantity by taking fuel consumption as target
Figure BDA0003271625370000041
And a corresponding latitude amplitude u t1 ,u t2 ,u t3 The method comprises the steps of carrying out a first treatment on the surface of the Aiming at the control target of the formation out-of-plane configuration parameter delta i, calculating the optimal out-of-plane single pulse control quantity of fuel>
Figure BDA0003271625370000042
And a corresponding latitude amplitude u n1
The calculation formula is as follows:
Figure BDA0003271625370000051
Figure BDA0003271625370000052
Figure BDA0003271625370000053
wherein Δa m =Δa target -Δa current ,Δe m =|Δe m |=|Δe target -Δe current And n is the angular velocity of the orbit of the main star, and T is the orbit period of the main star.
The corresponding latitude amplitude of application is:
Figure BDA0003271625370000054
Figure BDA0003271625370000055
the corresponding time is +.>
Figure BDA0003271625370000056
S103, calculating the total duration of continuous small thrust application at each pulse according to the speed increment which can be provided in the unit time of the maximum thrust of the single electric propeller;
in particular, based on the use of continuously small thrust to achieve equivalent pulse controlThe thinking of the control effect is based on the maximum thrust F of a single electric propeller max Speed increment DeltaV capable of being provided in unit time (1 s) max Calculating the total duration DeltaT of the application of successive low thrust forces at each sub-optimal pulse 1 ,ΔT 2 ,ΔT 3 Wherein
Figure BDA0003271625370000057
m is the current satellite mass. />
Figure BDA0003271625370000058
To achieve a speed pulse with continuously low thrust +.>
Figure BDA0003271625370000059
The required continuous low thrust on-time. And so on to obtain the total duration of continuous small thrust application at the other two suboptimal pulses +.>
Figure BDA00032716253700000510
And->
Figure BDA00032716253700000511
S104, considering the influence of the continuous small thrust application duration of each track on the formation configuration control time and the fuel consumption, setting the continuous small thrust application duration of each track and the attitude maneuver time which can generate the thrust in the required direction before and after each pulse application moment;
specifically, considering that only one point of time is required for pulse application and continuous small thrust is required to be continuously applied for a certain period of time, in order to achieve approximately equivalent pulse control effect with continuous small thrust, it is preferable that each small thrust control start time in the plane is
Figure BDA00032716253700000512
Ending time of +.>
Figure BDA00032716253700000513
Wherein k is 1,2,3.
S105, single electric pushing of the satellite is adjusted to a control direction in a gesture maneuvering time period;
in particular, considering that the satellite has only a single electric propulsion and that the electric propulsion is fixed on the satellite, in order to provide a continuous small thrust along the track, it is necessary to apply at each time a continuous small thrust is applied
Figure BDA00032716253700000514
Previously, the satellite needs to be reserved from the normal attitude q normal Motorized to a target attitude enabling a single electric propulsion to generate thrust along the track +.>
Figure BDA0003271625370000061
Time of (2)>
Figure BDA0003271625370000062
S106, applying continuous low-thrust control at the set continuous low-thrust application position of each track, and maneuvering the satellite to a conventional posture after the application of the continuous low-thrust control is finished;
specifically, considering that the number of satellite attitude maneuvers is reduced as much as possible, 3m consecutive small thrust application arcs along the track may be set as a period in which the satellite attitude remains in a target attitude that can produce consecutive small thrust along the track. Outside the period, the satellite performs attitude maneuver and returns to the normal attitude q normal
And S107, updating the current in-plane formation configuration parameters and returning to the step S101 until the in-plane formation configuration parameter deviation is judged to be within a preset range.
Specifically, after a continuous low thrust control period, the configuration parameters in the current formation plane are updated
Figure BDA0003271625370000063
And updating the deviation of the configuration parameters in the current formation plane +.>
Figure BDA0003271625370000064
If->
Figure BDA0003271625370000065
Then the formation control based on the single electric push is ended, otherwise, the step S101 is returned to, and the iterative control is performed until +.>
Figure BDA0003271625370000066
I.e. the deviation of the configuration parameters in the formation plane is within a preset allowable range.
In a second embodiment, the out-of-plane formation configuration control:
s201, according to the current out-of-plane formation configuration parameter delta i current And a target out-of-plane formation configuration parameter Δi target Calculating out-of-plane formation configuration parameter deviation delta i error
S202, judging that the deviation of the out-of-plane formation configuration parameters exceeds a preset threshold, namely delta i error >Δi predefined Calculating a fuel optimal out-of-plane optimal pulse size based on the current out-of-plane formation parameters, targeting a minimum fuel consumption
Figure BDA0003271625370000067
And the corresponding optimal pulse application latitude argument u n1
The control quantity of single pulse outside the track plane is as follows:
δv n1 =na|Δi target -Δi current |=na|Δi man |
the corresponding latitude amplitude of application is:
Figure BDA0003271625370000068
the corresponding time is +.>
Figure BDA0003271625370000069
In the above formula, δi is the configuration out-of-plane reconstruction control I vector change amount.
S203, calculating the total duration of continuous small thrust application at each pulse according to the speed increment which can be provided in the unit time of the maximum thrust of the single electric propeller;
specifically, based on the thought of realizing the equivalent pulse control effect by adopting continuous small thrust, the method is as followsAccording to the maximum thrust F of a single electric propeller max Speed increment DeltaV capable of being provided in unit time (1 s) max Calculating the total duration DeltaT of the application of successive low thrust forces at each optimal out-of-plane pulse n1 Wherein
Figure BDA0003271625370000071
m is the current satellite mass. />
Figure BDA0003271625370000072
To achieve a speed pulse with continuously low thrust +.>
Figure BDA0003271625370000073
The required continuous low thrust on-time.
S204, considering the influence of the continuous small thrust application duration of each track on the formation configuration control time and the fuel consumption, setting the continuous small thrust application duration of each track and the attitude maneuver time which can generate the thrust in the required direction before and after each pulse application moment;
specifically, considering that only one point in time is required for pulse application and continuous small thrust is required to be continuously applied for a period of time, in order to achieve approximately equivalent pulse control effect with continuous small thrust, the out-of-plane single small thrust control start time may be
Figure BDA0003271625370000074
Ending time of +.>
Figure BDA0003271625370000075
S205, single electric pushing of the satellite is adjusted to a control direction in the attitude maneuver time period;
in particular, considering that the satellite has only a single electric push and that the electric push is fixed on the satellite, in order to provide a continuous small thrust out of plane, it is necessary to apply at each time a continuous small thrust is applied
Figure BDA0003271625370000076
Previously, the satellite needs to be reserved from the normal attitude q normal Motorized toTarget posture enabling single electric propulsion to generate out-of-plane thrust +.>
Figure BDA0003271625370000077
Time of (2)>
Figure BDA0003271625370000078
S206, applying continuous low-thrust control at the set continuous low-thrust application position of each track, and maneuvering the satellite to a conventional posture after the application of the continuous low-thrust control is finished;
specifically, considering that the number of satellite attitude maneuvers is reduced as much as possible, p out-of-plane continuous low-thrust application arcs may be set as a period, and the satellite attitude in the period is kept at a target attitude that can generate out-of-plane continuous low thrust. Outside the period, the satellite performs attitude maneuver and returns to the normal attitude q normal
S207, updating the current out-of-plane formation configuration parameters and returning to the step S201 until the deviation of the out-of-plane formation configuration parameters is judged to be within a preset range.
Specifically, after a continuous low thrust control period, the current formation out-of-plane configuration parameter Δi is updated current And update the current formation out-of-plane configuration parameter deviation deltai error If Δi error ≤(Δi error ) predefined If the formation control based on the single electric push is finished, otherwise, returning to step S201, and performing iterative control until delta i is reached error ≤(Δi error ) predefined I.e. the deviation of the formation out-of-plane configuration parameters is within a preset allowable range.
While the preferred embodiment of the present invention has been described in detail, the invention is not limited to the embodiment, and various equivalent modifications and substitutions can be made by those skilled in the art without departing from the spirit of the invention, and these modifications and substitutions are intended to be included in the scope of the present invention as defined in the appended claims.

Claims (8)

1. The control method for the formation configuration of the small satellite based on single electric push is characterized by comprising the following steps of:
s1, calculating formation configuration parameter deviation according to a current formation configuration parameter and a target formation configuration parameter;
s2, judging that the deviation of formation configuration parameters exceeds a preset threshold, and calculating the in-plane formation configuration control suboptimal three-pulse size, the corresponding suboptimal pulse application latitude amplitude angle, the out-of-plane formation configuration control optimal pulse size and the corresponding optimal pulse application latitude amplitude angle according to the current formation configuration parameters by taking less fuel consumption as a target;
s3, calculating the total duration of continuous small thrust application at each pulse according to the speed increment which can be provided in the unit time of the maximum thrust of the single electric propeller;
s4, considering the influence of the continuous small thrust application duration of each track on the formation configuration control time and the fuel consumption, setting the continuous small thrust application duration of each track and the attitude maneuver time which can generate the thrust in the required direction before and after each pulse application moment;
s5, maneuvering the satellite attitude to a target attitude capable of enabling a single electric propulsion to generate a preset thrust direction in an attitude maneuvering time period;
s6, applying continuous low-thrust control at the set continuous low-thrust application position of each track, and maneuvering the satellite to a conventional posture after the application of the continuous low-thrust control is finished;
and S7, updating the current formation configuration parameters and returning to the step S1 until the deviation of the formation configuration parameters is judged to be within a preset range.
2. The method for controlling the formation configuration of the small satellites based on single electric propulsion according to claim 1, wherein the deviation of the formation configuration parameters comprises the difference between the auxiliary satellites and the semi-long axes of the main satellites, the relative eccentricity vector, the relative inclination vector and the formation offset distance along the track around the center of flight.
3. The method for controlling the formation configuration of the small satellite based on single electric push according to claim 1, wherein the calculation formula of the in-plane formation configuration for controlling the suboptimal three pulse size is as follows:
Figure FDA0004180857570000011
Figure FDA0004180857570000012
Figure FDA0004180857570000013
in the above-mentioned method, the step of,
Figure FDA0004180857570000014
representing the first pulse size in plane, +.>
Figure FDA0004180857570000015
Representing in-plane second pulse size, < >>
Figure FDA0004180857570000016
Representing the third pulse size in plane, Δa m Indicating the expected suboptimal pulse control of the second half long axis change in ignition, Δλ sm Representing the amount of change in the constant offset along the track, Δe m The magnitude change of the E vector is represented, n represents the angular velocity of the primary star orbit, and T represents the period of the formation primary star orbit.
4. A control method for a formation configuration of a small satellite based on single electric propulsion according to claim 3, wherein the calculation formula of the latitude angle of application of the corresponding suboptimal pulse is as follows:
Figure FDA0004180857570000021
Figure FDA0004180857570000022
Figure FDA0004180857570000023
/>
in the above-mentioned method, the step of,
Figure FDA0004180857570000024
representing the corresponding latitude depression angle of the application of the first pulse in the plane, < >>
Figure FDA0004180857570000025
Representing the corresponding latitude depression angle of the application of the second pulse in the plane,/>
Figure FDA0004180857570000026
Representing the corresponding latitude depression angle of the application when the third pulse in the plane is applied,
Figure FDA0004180857570000027
component representing the desired change in the y-direction of the relative eccentricity vector,/->
Figure FDA0004180857570000028
Representing the component of the desired change in the x-direction of the relative eccentricity vector.
5. The method for controlling the formation configuration of the small satellite based on single electric push as claimed in claim 4, wherein the calculation formula of the optimal pulse size for controlling the formation configuration outside the plane is as follows:
δv n1 =na|Δi target -Δi current |=na|Δi man |
in the above, δv n1 Represents the out-of-plane optimal pulse size, n represents the primary satellite orbit angular velocity, a represents the primary satellite orbit semi-major axis, Δi man Representing the relative tilt vector, Δi, of the desired change target Representing the out-of-plane formation configuration parameters of the object, Δi current Representing current out-of-plane formation configuration parameters.
6. The control method of the small satellite formation configuration based on single electric push according to claim 5, wherein the calculation formula of the corresponding optimal pulse application latitude amplitude angle is as follows:
Figure FDA0004180857570000029
in the above, u n1 Representing the corresponding optimal pulse application latitude argument,
Figure FDA00041808575700000210
component representing the desired change in the y-direction of the relative tilt vector, is provided>
Figure FDA00041808575700000211
Representing the component of the desired change in the x-direction of the relative tilt vector.
7. The control method of the formation configuration of the small satellites based on single electric propulsion according to claim 6, wherein the calculation formula of the total duration of the continuous small thrust application at each pulse is as follows:
Figure FDA00041808575700000212
Figure FDA00041808575700000213
Figure FDA0004180857570000031
Figure FDA0004180857570000032
in the above-mentioned method, the step of,
Figure FDA0004180857570000033
indicates the duration of the force application required to achieve the first pulse size,/->
Figure FDA0004180857570000034
Indicates the duration of the force required to achieve the second pulse size,/->
Figure FDA0004180857570000035
Indicating the duration of the force applied, deltaT, required to achieve the third pulse size n1 Represents the duration of the applied force, deltaV, required to achieve an out-of-plane optimum pulse size max Representing the speed increase that can be provided per unit time of maximum thrust of a single electric propulsion.
8. The control method for a single electric propulsion-based minisatellite formation configuration according to claim 7, wherein the formula for setting the duration of each track of continuous small propulsion application before and after each pulse application time is as follows:
Figure FDA0004180857570000036
Figure FDA0004180857570000037
Figure FDA0004180857570000038
/>
Figure FDA0004180857570000039
in the above-mentioned method, the step of,
Figure FDA00041808575700000310
indicates the start time of each low thrust control in the plane, < > and the like>
Figure FDA00041808575700000311
Indicating the end time of each low thrust control in the plane, < > and the like>
Figure FDA00041808575700000312
Indicates the in-plane force application time length, (T) n1 ) start The out-of-plane single low thrust control start time is represented by (T) n1 ) end And the out-of-plane single small thrust control end time is indicated. />
CN202111104268.6A 2021-09-22 2021-09-22 Small satellite formation configuration control method based on single electric push Active CN113703487B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111104268.6A CN113703487B (en) 2021-09-22 2021-09-22 Small satellite formation configuration control method based on single electric push

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111104268.6A CN113703487B (en) 2021-09-22 2021-09-22 Small satellite formation configuration control method based on single electric push

Publications (2)

Publication Number Publication Date
CN113703487A CN113703487A (en) 2021-11-26
CN113703487B true CN113703487B (en) 2023-05-26

Family

ID=78661558

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111104268.6A Active CN113703487B (en) 2021-09-22 2021-09-22 Small satellite formation configuration control method based on single electric push

Country Status (1)

Country Link
CN (1) CN113703487B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114935947B (en) * 2022-07-21 2022-10-25 银河航天(北京)网络技术有限公司 Method for keeping control of satellite formation and electronic equipment
CN115309174A (en) * 2022-07-21 2022-11-08 银河航天(北京)网络技术有限公司 Method and electronic equipment for combined control of near-earth orbit satellite formation
CN116750210B (en) * 2023-07-12 2024-02-20 银河航天(北京)网络技术有限公司 Method for controlling orbit of satellite formation

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101941528A (en) * 2010-09-30 2011-01-12 哈尔滨工业大学 Flywheel based attitude maneuvering control device and method for successive approaching of satellite rounding instantaneous Euler shaft
CN106628261A (en) * 2016-10-20 2017-05-10 上海航天控制技术研究所 Satellite high-precision and high-stability attitude control method in electric propulsion position protecting process
CN108438255A (en) * 2018-03-14 2018-08-24 上海航天控制技术研究所 Satellite is diversion Formation Configuration initial method under a kind of engineering constraints

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6210749B2 (en) * 2013-06-18 2017-10-11 三菱電機株式会社 Spacecraft orbit control method and apparatus
US9334068B2 (en) * 2014-04-04 2016-05-10 NOA Inc. Unified orbit and attitude control for nanosatellites using pulsed ablative thrusters
CN104142686B (en) * 2014-07-16 2016-06-29 北京控制工程研究所 A kind of satellite Autonomous formation flight control method
CN106094529B (en) * 2016-07-15 2018-08-17 上海航天控制技术研究所 The in-orbit Auto-calibration method of thruster under formation task multiple-pulse control condition
CN109977576A (en) * 2019-04-03 2019-07-05 北京理工大学 A kind of constellation system multidisciplinary design optimization based on agent model
CN112607065B (en) * 2020-12-23 2023-04-18 长春工业大学 High-precision phase control method based on electric propulsion system
CN113343442B (en) * 2021-05-21 2022-01-25 中国科学院国家空间科学中心 Method and system for solving fixed-time finite fuel multi-pulse transfer orbit

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101941528A (en) * 2010-09-30 2011-01-12 哈尔滨工业大学 Flywheel based attitude maneuvering control device and method for successive approaching of satellite rounding instantaneous Euler shaft
CN106628261A (en) * 2016-10-20 2017-05-10 上海航天控制技术研究所 Satellite high-precision and high-stability attitude control method in electric propulsion position protecting process
CN108438255A (en) * 2018-03-14 2018-08-24 上海航天控制技术研究所 Satellite is diversion Formation Configuration initial method under a kind of engineering constraints

Also Published As

Publication number Publication date
CN113703487A (en) 2021-11-26

Similar Documents

Publication Publication Date Title
CN113703487B (en) Small satellite formation configuration control method based on single electric push
Xing et al. Nonlinear attitude state tracking control for spacecraft
Shen et al. Dynamic lateral entry guidance logic
US20100250031A1 (en) Propulsive guidance for atmospheric skip entry trajectories
CN109774974B (en) Track design method for space debris approaching
CN112947573B (en) Reentry guidance method for hypersonic aircraft under terminal time constraint
CN109774973B (en) Design method for ascending intersection orbit parameters of manned lunar surface lander
CN109911249B (en) Interstellar transfer limited thrust orbit-entering iterative guidance method for low thrust-weight ratio aircraft
CN106444822A (en) Space vector field guidance based stratospheric airship&#39;s trajectory tracking control method
CN107515611B (en) Superposition type hybrid sine maneuvering path planning method
CN109839116B (en) Intersection approach method with minimum speed increment
CN112256061A (en) Reentry guidance method for hypersonic aircraft under complex environment and task constraint
CN112572835B (en) Satellite in-orbit angular momentum management and control method with attitude switching function
CN113602532A (en) Solid carrier rocket in-orbit correction method
CN109190155B (en) Hybrid continuous low-thrust track design method adopting electric propulsion/solar sail propulsion
CN112180961B (en) Track tracking control method and system for all-state limited stratospheric airship
CN113835442A (en) Hypersonic glide aircraft linear pseudo-spectrum reentry guidance method and system
CN113343561A (en) Method and system for solving optimal moon fly-by transfer orbit
CN110162070B (en) Three-axis attitude motion trajectory planning system and method under constraint of free boundary of tail end
CN103818564A (en) Integrated control method adopting low-thrust spacecraft orbit maintaining and to-earth orientation posture keeping
CN114063645B (en) Eccentricity inclination angle vector-based inclination flying-around holding control effect evaluation method
Nan et al. Global 4-D trajectory optimization for spacecraft
CN112455725A (en) Method for transferring and converting pulse orbit transfer direction to limited thrust orbit
KR101920637B1 (en) method for controlling of satellite for large angle maneuver using satellite angular velocity
Ross et al. Bellman 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