CN113734469A - Method and system for distributing thrust for maintaining position of electric thruster - Google Patents
Method and system for distributing thrust for maintaining position of electric thruster Download PDFInfo
- Publication number
- CN113734469A CN113734469A CN202111006812.3A CN202111006812A CN113734469A CN 113734469 A CN113734469 A CN 113734469A CN 202111006812 A CN202111006812 A CN 202111006812A CN 113734469 A CN113734469 A CN 113734469A
- Authority
- CN
- China
- Prior art keywords
- electric
- electric thruster
- thruster
- calculating
- thrust
- 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
- 238000000034 method Methods 0.000 title claims abstract description 47
- 239000013598 vector Substances 0.000 claims abstract description 62
- 230000007246 mechanism Effects 0.000 claims abstract description 32
- 238000004364 calculation method Methods 0.000 claims description 16
- 239000004576 sand Substances 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 abstract description 4
- 238000012423 maintenance Methods 0.000 abstract description 3
- 238000004088 simulation Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000446 fuel Substances 0.000 description 2
- 230000002035 prolonged effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/242—Orbits and trajectories
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
- B64G1/245—Attitude control algorithms for spacecraft attitude control
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/40—Arrangements or adaptations of propulsion systems
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
The invention discloses a method and a system for distributing thrust of position maintenance of an electric thruster, wherein the method comprises the following steps: determining normal pointing gears of the two electric thrusters according to the position holding drift rate control quantity; respectively calculating to obtain a target thrust direction [ X ] of the first electric thruster according to the determined normal pointing gears of the two electric thrustersn,Yn,Zn]And a boot time period Δ tnAnd a target thrust direction [ X ] of the second electric thrusters,Ys,Zs]And a boot time period Δ ts(ii) a Calculating to obtain the rotating speed omega of the first theoretical vector adjusting mechanismnAnd a second theoretical vector adjusting mechanism rotating speed omegas(ii) a Will [ X ]n,Yn,Zn]And Δ tnIs distributed to the first electric thruster to drive [ Xs,Ys,Zs]And Δ tsIs distributed to the second electric thruster to drive omeganIs assigned to a first vector adjustment mechanismωsIs assigned to the second vector adjustment mechanism. The method adopts a complete analysis method to calculate the starting time length and the speed increment of the electric thruster, fully considers the practical engineering constraint and can meet the position maintenance task requirement of the geostationary orbit satellite.
Description
Technical Field
The invention belongs to the technical field of satellite attitude and orbit control, and particularly relates to a method and a system for distributing position keeping thrust of an electric thruster.
Background
At present, the electric propulsion satellite in orbit at home and abroad can be used for carrying out east-west and north-south combined position keeping control by using four electric thrusters or can be used for carrying out north-south position keeping control by using two electric thrusters, and if the above propulsion mode is directly applied to the simply-configured electric propulsion satellite with the layout of two mechanical arms and two thrusters, the following defects are achieved: the thruster is started more times, and the service life of the thruster and the actual constraint of a fuel consumption mechanical arm type vector adjusting mechanism are possible.
Therefore, a position maintaining thrust distribution method considering the starting times and the limitation of a vector adjusting mechanism for the electric propulsion satellite with a simple configuration needs to be designed, so that the service life of the electric propulsion satellite with the stationary orbit is prolonged, and the reliability is improved.
Disclosure of Invention
The technical problem of the invention is solved: the method and the system for distributing the position keeping thrust of the electric thruster overcome the defects of the prior art, calculate the starting time and the speed increment of the electric thruster by adopting a complete analysis method, fully consider the practical engineering constraints and meet the position keeping task requirement of a stationary orbit satellite.
In order to solve the technical problem, the invention discloses a thrust distribution method for maintaining the position of an electric thruster, which comprises the following steps:
determining normal pointing gears of the two electric thrusters according to the position holding drift rate control quantity;
respectively calculating to obtain a target thrust direction [ X ] of the first electric thruster according to the determined normal pointing gears of the two electric thrustersn,Yn,Zn]And a boot time period Δ tnAnd a target thrust direction [ X ] of the second electric thrusters,Ys,Zs]And a boot time period Δ ts;
According to [ X ]n,Yn,Zn]And [ X ]s,Ys,Zs]And calculating to obtain the rotating speed omega of the first theoretical vector adjusting mechanismnAnd a second theoretical vector adjusting mechanism rotating speed omegas;
Will calculate the obtained [ Xn,Yn,Zn]And Δ tnDistributing the calculated [ X ] to the first electric thrusters,Ys,Zs]And Δ tsDistributing the calculated omega to a second electric thrusternAssigning the calculated ω to the first vector adjustment mechanismsIs assigned to the second vector adjustment mechanism.
In the above method for distributing position-maintaining thrust of an electric thruster, determining normal directional shift positions of two electric thrusters according to a position-maintaining drift rate control amount includes:
determining a unit vector [ X, Y, Z ] of a thrust direction under a track coordinate system; the origin of the orbit coordinate system is the center of mass of the satellite, the X axis points to the advancing direction of the satellite, the Z axis points to the ground, and the Y axis, the X axis and the Z axis meet the right-hand rule;
calculating to obtain the adjusting range C of the vector adjusting mechanismy=Y/Z;
For the regulation range CyN are equally divided to obtain n normal direction gears C which are arranged from small to largeys:Cys=Cy1、Cy2、…、Cyn(ii) a Wherein each normal direction is towards gear Cy1、Cy2、…、CynThe corresponding allowable ranges of the drift rate are as follows: [ Delta DL1,ΔDH1]、[ΔDL2,ΔDH2]、…、[ΔDLn,ΔDHn];
Acquiring a position holding drift rate control quantity delta D;
if Δ D ∈ [ Δ DL ]1,ΔDH1]Then, the normal direction gear of the two electric thrusters is determined to be Cys=Cy1;
If Δ D ∈ [ Δ DL ]2,ΔDH2]Then, the normal direction gear of the two electric thrusters is determined to be Cys=Cy2;
And so on;
if Δ D ∈ [ Δ DL ]n,ΔDHn]Then, the normal direction gear of the two electric thrusters is determined to be Cys=Cyn。
In the method for distributing the position maintaining thrust of the electric thrusters, the target thrust direction [ X ] of the first electric thruster is respectively calculated according to the determined normal directional gears of the two electric thrustersn,Yn,Zn]And a boot time period Δ tnAnd a target thrust direction [ X ] of the second electric thrusters,Ys,Zs]And a boot time period Δ tsThe method comprises the following steps:
acquiring eccentricity control vectors delta e in the x direction and the y direction under a track coordinate systemxAnd Δ eyAnd tilt control vectors Δ i in the x-and y-directionsxAnd Δ iy;
According to Δ ex、Δey、ΔixAnd Δ iyCalculating to obtain the first electric pushVelocity delta VN in the direction of force multiplier xxAnd a speed increment VS in the x-direction of the second electric thrusterxAnd the radial distribution identification flagV of the two electric thrusterszAnd the radial distribution difference Vz0;
According to flagVzAnd Vz0And calculating to obtain a speed increment VN of the first electric thruster in the y directionyAnd velocity increment VN in the z-directionzAnd a speed increment VS in the y-direction of the second electric thrusteryAnd z-direction velocity increment VSz;
According to VNy、VNz、VSyAnd VSzAnd calculating to obtain a target thrust direction [ X ] of the first electric thrustern,Yn,Zn]And a boot time period Δ tnAnd a target thrust direction [ X ] of the second electric thrusters,Ys,Zs]And a boot time period Δ ts。
In the above-described electric thruster position maintaining thrust force distributing method, VNxAnd VSxThe calculation flow of (2) is as follows:
according to the nominal orbital velocity V of the geostationary orbitsAnd the earth static orbit nominal orbit angular velocity omegaeAnd calculating to obtain a constant A, a constant B and an intermediate parameter aa:
then there are:
speed increment VN of first electric thruster in x directionxComprises the following steps:
speed increment VS of second electric thruster in x directionxComprises the following steps:
in the above-described electric thruster position holding thrust force distribution method, the flagVzAnd Vz0The calculation flow of (2) is as follows:
calculating to obtain an intermediate variable dVztemp:
According to dVztempDetermining flagVz:
If dVztempIf the value is less than or equal to 0, flagVz=1
If dVztempIf > 0, flagVz=-1
Calculating to obtain a radial distribution difference value Vz0:
Wherein, dVzIndicating the difference, dV, between the radial control quantities of the two electric thrustersz=abs(dVztemp),Vy0And the normal control quantity of the unilateral electric thruster is shown.
In the above-described electric thruster position maintaining thrust force distributing method, VNy、VNz、VSyAnd VSzThe calculation flow of (2) is as follows:
calculating to obtain an intermediate variable Vytemp:
When flagV is presentzWhen 1, there are:
VNy=Vytemp
VNz=Vz0+dVz
VSy=Vy0-Vytemp
VSz=Vz0
when flagV is presentzWhen-1, there is:
VNy=Vy0-Vytemp
VNz=Vz0
VSy=Vytemp
VSz=Vz0+dVz。
in the above-described electric thruster position holding thrust force distribution method, [ X [ X ] ]n,Yn,Zn]、Δtn、[Xs,Ys,Zs]And Δ tsThe calculation flow of (2) is as follows:
according to VNx、VNyAnd VNzCalculating to obtain [ X ]n,Yn,Zn]And Δ tn:
According to VSx、VSyAnd VSzCalculating to obtain [ X ]s,Ys,Zs]And Δ ts:
WhereinF represents the magnitude of the thrust of the electric thruster, VeRepresenting the electrical thruster specific impulse, mass representing the satellite mass.
In the above-described electric thruster position holding thrust force distribution method, the thrust force is distributed in accordance with [ X ]n,Yn,Zn]And [ X ]s,Ys,Zs]And calculating to obtain the rotating speed omega of the first theoretical vector adjusting mechanismnAnd a second theoretical vector adjusting mechanism rotating speed omegasThe method comprises the following steps:
determining a current heading [ X ] of a first electric thrustern0,Yn0,Zn0]And calculating to obtain [ X ]n,Yn,Zn]And [ X ]n0,Yn0,Zn0]Angle theta therebetweenn;
Determining a current heading [ X ] of a second electric thrusters0,Ys0,Zs0]And calculating to obtain [ X ]s,Ys,Zs]And [ X ]s0,Ys0,Zs0]Angle theta therebetweens;
Obtaining the rotating speed omega of the maximum theoretical vector adjusting mechanismmax;
When theta isn≥θsIn time, there are:
ωn=ωmax
when theta isn<θsIn time, there are:
ωs=ωmax。
in the method for distributing the thrust for maintaining the position of the electric thruster, a first vector adjusting mechanism is arranged on the south side of a satellite and connected with the first electric thruster for realizing the adjustment of the thrust direction in the x direction, the y direction and the z direction; the second vector adjusting mechanism is arranged on the north side of the satellite and connected with the second electric thruster, and is used for adjusting thrust direction in the x direction, the y direction and the z direction.
Correspondingly, the invention also discloses a position maintaining thrust distribution system of the electric thruster, which comprises:
the gear determining module is used for determining normal pointing gears of the two electric thrusters according to the position holding drift rate control quantity;
a first calculation module, configured to calculate, according to the determined normal pointing gears of the two electric thrusters, a target thrust direction [ X ] of the first electric thrustern,Yn,Zn]And a boot time period Δ tnAnd a target thrust direction [ X ] of the second electric thrusters,Ys,Zs]And a boot time period Δ ts;
A second calculation module for calculating according to [ Xn,Yn,Zn]And [ X ]s,Ys,Zs]And calculating to obtain the rotating speed omega of the first theoretical vector adjusting mechanismnAnd a second theoretical vector adjusting mechanism rotating speed omegas;
An allocation module for allocating the calculated [ Xn,Yn,Zn]And Δ tnDistributing the calculated [ X ] to the first electric thrusters,Ys,Zs]And Δ tsDistributing the calculated omega to a second electric thrusternAssigning the calculated ω to the first vector adjustment mechanismsIs assigned to the second vector adjustment mechanism.
The invention has the following advantages:
(1) the invention discloses a thrust distribution scheme for maintaining the position of an electric thruster, wherein each position protection control cycle of a south electric thruster and a north electric thruster is respectively started once, the starting time is short, and the service life of the electric thruster is prolonged.
(2) The invention discloses a thrust allocation scheme for maintaining the position of an electric thruster, which has less calculation amount and reduces the on-satellite calculation resource requirement.
(3) The invention discloses a thrust distribution scheme for maintaining the position of an electric thruster, which considers the moving range and the rotating speed bearing capacity of a vector adjusting mechanism and is safer and more reliable.
Drawings
FIG. 1 is a flow chart illustrating the steps of a method for distributing position maintaining thrust of an electric thruster in accordance with an embodiment of the present invention;
FIG. 2 is a schematic diagram of a layout of an electric propulsion satellite thruster according to an embodiment of the present invention;
FIG. 3 is a graph of a simulation example-geographical latitude and longitude curves in an embodiment of the present invention;
fig. 4 is a block diagram of a position maintaining thrust force distribution system of an electric thruster in an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the embodiments of the present invention will be described in detail with reference to the accompanying drawings.
The invention discloses a method for distributing position keeping thrust of an electric thruster, which is a simply configured method for distributing position keeping thrust of an electric propulsion satellite. Firstly, selecting a proper normal direction of the electric thruster according to the requirement of the position keeping drift rate control quantity; then, acquiring three-dimensional directions and starting-up time lengths of the two electric thrusters corresponding to the determined normal directions; and finally, distributing the theoretical rotating speeds of the two vector adjusting mechanisms according to the thrust distribution mode of the two electric thrusters. The invention adopts a complete analysis mode, solves the problem of the position maintaining control distribution of the electric thruster under the layout of the double mechanical arms and the electric thruster, reduces the fuel consumption and reduces the influence of the motion of the vector adjusting mechanism on the satellite.
As shown in fig. 1, in the present embodiment, the method for distributing position maintaining thrust of an electric thruster includes:
and step 101, determining normal pointing gears of the two electric thrusters according to the position holding drift rate control quantity.
In the present embodiment, the electric thruster position maintaining thrust force distribution method can be applied to, but is not limited to, a simple configuration of an electric propulsion satellite. As shown in fig. 2, the layout of the simple configuration electric propulsion satellite is as follows: two electric thrusters and two vector adjusting mechanisms. The first vector adjusting mechanism is arranged on the south side of the satellite, is connected with a first electric thruster and is used for realizing adjustment of thrust direction in the x direction, the y direction and the z direction; the second vector adjusting mechanism is arranged on the north side of the satellite and connected with the second electric thruster, and is used for adjusting thrust direction in the x direction, the y direction and the z direction. Namely, the first vector adjusting mechanism and the second vector adjusting mechanism are respectively and obliquely arranged on the north and south sides of the satellite back floor and are respectively used for realizing the adjustment of the three-dimensional direction of the first electric thruster and the second electric thruster. The method for distributing the position maintaining thrust of the electric thruster can be executed in one orbit control period, and ignition is successively carried out in the one position maintaining control period, so that the satellite can jointly adjust six orbits.
Preferably, a position hold drift rate control quantity, denoted as Δ D, may be determined. Wherein:
if Δ D ∈ [ Δ DL ]1,ΔDH1]Then, the normal direction gear of the two electric thrusters is determined to be Cys=Cy1;
If Δ D ∈ [ Δ DL ]2,ΔDH2]Then, the normal direction gear of the two electric thrusters is determined to be Cys=Cy2;
And so on;
if Δ D ∈ [ Δ DL ]n,ΔDHn]Then, the normal direction gear of the two electric thrusters is determined to be Cys=Cyn。
Further, the allowable range of the drift rate corresponding to each normal directional gear is determined as follows: determining unit vectors [ X, Y, Z ] of thrust directions under a orbital coordinate system](ii) a Calculating to obtain the adjusting range C of the vector adjusting mechanismyY/Z; for the regulation range CyN are equally divided to obtain n normal direction gears C which are arranged from small to largeys:Cys=Cy1、Cy2、…、Cyn(ii) a Each normal direction gear Cy1、Cy2、…、CynThe corresponding allowable ranges of the drift rate are as follows: [ Delta DL1,ΔDH1]、[ΔDL2,ΔDH2]、…、[ΔDLn,ΔDHn]. Wherein, the railThe origin of the orbit coordinate system is the center of mass of the satellite, the X axis points to the advancing direction of the satellite, the Z axis points to the ground, and the Y axis, the X axis and the Z axis meet the right-hand rule.
102, respectively calculating to obtain a target thrust direction [ X ] of the first electric thruster according to the determined normal pointing gears of the two electric thrustersn,Yn,Zn]And a boot time period Δ tnAnd a target thrust direction [ X ] of the second electric thrusters,Ys,Zs]And a boot time period Δ ts。
In this embodiment, the step 102 may specifically include the following sub-steps:
substep 1021, obtaining eccentricity control vector delta e in x direction and y direction of track coordinate systemxAnd Δ eyAnd tilt control vectors Δ i in the x-and y-directionsxAnd Δ iy。
Substep 1022, based on Δ ex、Δey、ΔixAnd Δ iyAnd calculating to obtain a speed increment VN of the first electric thruster in the x directionxAnd a speed increment VS in the x-direction of the second electric thrusterxAnd the radial distribution identification flagV of the two electric thrusterszAnd the radial distribution difference Vz0。
In the present exemplary embodiment, the speed increment VN in the x direction of the first electric thrusterxAnd a speed increment VS in the x-direction of the second electric thrusterxThe solving formulas of (1) are respectively as follows:
preferably, the nominal orbital velocity V can be determined according to the earth's stationary orbitsAnd the earth static orbit nominal orbit angular velocity omegaeAnd calculating to obtain a constant A, a constant B and an intermediate parameter aa:
further, the flag V is allocated radiallyzAnd the radial distribution difference Vz0The solution flow of (a) may be as follows:
first, an intermediate variable dV is calculatedztemp:
Then, according to the intermediate variable dVztempDetermining flagVz:
If dVztempIf the value is less than or equal to 0, flagVz=1
If dVztempIf > 0, flagVz=-1
Finally, calculating to obtain a radial distribution difference value Vz0:
Wherein, dVzIndicating the difference, dV, between the radial control quantities of the two electric thrustersz=abs(dVztemp),Vy0And the normal control quantity of the unilateral electric thruster is shown.
Substep 1023, based on flagVzAnd Vz0And calculating to obtain a speed increment VN of the first electric thruster in the y directionyAnd velocity increment VN in the z-directionzAnd a speed increment VS in the y-direction of the second electric thrusteryAnd z-direction velocity increment VSz。
In this embodiment, when flagVzWhen 1, there are: VNy=Vytemp、VNz=Vz0+dVz、VSy=Vy0-Vytemp、VSz=Vz0(ii) a When flagV is presentzWhen-1, there is: VNy=Vy0-Vytemp、VNz=Vz0、VSy=Vytemp、VSz=Vz0+dVz。
Preferably, for the intermediate variable VytempThe method comprises the following steps:
substep 1024, based on VNy、VNz、VSyAnd VSzAnd calculating to obtain a target thrust direction [ X ] of the first electric thrustern,Yn,Zn]And a boot time period Δ tnAnd a target thrust direction [ X ] of the second electric thrusters,Ys,Zs]And a boot time period Δ ts。
In this embodiment, VN may be first determinedx、VNyAnd VNzCalculating to obtain [ X ]n,Yn,Zn]And Δ tn:
Then according to VSx、VSyAnd VSzCalculating to obtain [ X ]s,Ys,Zs]And Δ ts:
Wherein F represents the thrust magnitude of the electric thruster, VeThe specific impulse of the electric thruster is shown,mass represents satellite quality.
In this embodiment, the current orientation [ X ] of the first electric thruster may be determined firstn0,Yn0,Zn0]And current orientation [ X ] of second electric thrusters0,Ys0,Zs0](ii) a And calculating to obtain [ Xn,Yn,Zn]And [ X ]n0,Yn0,Zn0]Angle theta therebetweenn,[Xs,Ys,Zs]And [ X ]s0,Ys0,Zs0]Angle theta therebetweens。
Then, the maximum theoretical vector adjusting mechanism rotating speed is obtained and recorded as omegamax。
When theta isn≥θsIn time, there are:
ωn=ωmax
when theta isn<θsIn time, there are:
ωs=ωmax
In conclusion, the method adopts a complete analysis method to calculate the starting time length and the speed increment of the electric thruster, fully considers the actual engineering constraint and can meet the requirement of the position maintenance task of the geostationary orbit satellite.
On the basis of the above embodiments, the following description will be made in conjunction with several specific examples.
For example, the thrust of a single electric thruster is 100mN, the mass of a satellite is 1380kg, and the fixed point longitude is 114E of east longitude; the position of the satellite keeps the control cycle to be 1 day; Δ D ═ 0.0043 °/day, Δ ex=0,Δey=-2.29e-5/day,Δix=0,Δiy=-0.0029°/day;Cy0.7. Then, a suitable normal pointing gear C can be selectedys=[1.7,1.0,0.7]The corresponding allowable ranges of the drift rate are as follows: [ -0.001,0.001]、[-0.003,0.003]、[-0.010,0.010]。
Based on the above known quantities in combination with the above steps, one can calculate: VNz=-0.237,VSx=0.011,dVztemp=-0.0014,dVz=abs(dVztemp)=0.0014,flagVz=1,Vytemp=0.079,Vz0=0.112。
Further, there are: VNz=0.114,VNy=0.079,VSz=0.112,VSy=0.078。
Then the process of the first step is carried out,
finally, according to the current direction, omega can be calculatedn=1°/s,ωs=0.54°/s。
Further, based on the above conditions: the thrust of a single electric thruster is 100mN, the mass of a satellite is 1380kg, and the fixed point longitude is 114 degrees E of east longitude; the position holding control period of the satellite is 1 day. The simulation duration is taken as 1 year. Then the simulation result as shown in fig. 3 can be obtained, and the thruster allocation method provided by the invention can keep the satellite within the 'dead zone' range of the fixed point longitude and latitude ± 0.5 °. Through simulation verification, the method can realize position keeping control electric thruster distribution under the simple configuration electric propulsion layout, and can be used for geosynchronous orbit position keeping control.
In conclusion, the thrust distribution method for maintaining the position of the electric thruster can distribute the thrust according to the accessibility of the vector adjusting mechanism and the rotation speed constraint, and realize the position maintaining control of the geosynchronous orbit.
On the basis of the above embodiment, as shown in fig. 4, the present invention also discloses an electric thruster position maintaining thrust distribution system, including: and a gear determining module 401, configured to determine normal pointing gears of the two electric thrusters according to the position holding drift rate control amount. A first calculating module 402, configured to respectively calculate, according to the determined normal directional gears of the two electric thrusters, a target thrust direction [ X ] of the first electric thrustern,Yn,Zn]And a boot time period Δ tnAnd a target thrust direction [ X ] of the second electric thrusters,Ys,Zs]And a boot time period Δ ts. A second calculation module 403 for calculating according to [ Xn,Yn,Zn]And [ X ]s,Ys,Zs]And calculating to obtain the rotating speed omega of the first theoretical vector adjusting mechanismnAnd a second theoretical vector adjusting mechanism rotating speed omegas. An assignment module 404 for assigning the calculated [ X ]n,Yn,Zn]And Δ tnDistributing the calculated [ X ] to the first electric thrusters,Ys,Zs]And Δ tsAssigned to a second electric thrusterWill calculate the resulting ωnAssigning the calculated ω to the first vector adjustment mechanismsIs assigned to the second vector adjustment mechanism.
For the system embodiment, since it corresponds to the method embodiment, the description is relatively simple, and for the relevant points, refer to the description of the method embodiment section.
Although the present invention has been described with reference to the preferred embodiments, it is not intended to limit the present invention, and those skilled in the art can make variations and modifications of the present invention without departing from the spirit and scope of the present invention by using the methods and technical contents disclosed above.
Those skilled in the art will appreciate that the invention may be practiced without these specific details.
Claims (10)
1. A method of distributing position maintaining thrust of an electric thruster, comprising:
determining normal pointing gears of the two electric thrusters according to the position holding drift rate control quantity;
respectively calculating to obtain a target thrust direction [ X ] of the first electric thruster according to the determined normal pointing gears of the two electric thrustersn,Yn,Zn]And a boot time period Δ tnAnd a target thrust direction [ X ] of the second electric thrusters,Ys,Zs]And a boot time period Δ ts;
According to [ X ]n,Yn,Zn]And [ X ]s,Ys,Zs]And calculating to obtain the rotating speed omega of the first theoretical vector adjusting mechanismnAnd a second theoretical vector adjusting mechanism rotating speed omegas;
Will calculate the obtained [ Xn,Yn,Zn]And Δ tnDistributing the calculated [ X ] to the first electric thrusters,Ys,Zs]And Δ tsDistributing the calculated omega to a second electric thrusternAssigning the calculated ω to the first vector adjustment mechanismsIs assigned to the second vector adjustment mechanism.
2. The method of claim 1, wherein determining the normally-oriented gears of the two electric thrusters based on the position hold drift rate control comprises:
determining a unit vector [ X, Y, Z ] of a thrust direction under a track coordinate system; the origin of the orbit coordinate system is the center of mass of the satellite, the X axis points to the advancing direction of the satellite, the Z axis points to the ground, and the Y axis, the X axis and the Z axis meet the right-hand rule;
calculating to obtain the adjusting range C of the vector adjusting mechanismy=Y/Z;
For the regulation range CyN are equally divided to obtain n normal direction gears C which are arranged from small to largeys:Cys=Cy1、Cy2、…、Cyn(ii) a Wherein each normal direction is towards gear Cy1、Cy2、…、CynThe corresponding allowable ranges of the drift rate are as follows: [ Delta DL1,ΔDH1]、[ΔDL2,ΔDH2]、…、[ΔDLn,ΔDHn];
Acquiring a position holding drift rate control quantity delta D;
if Δ D ∈ [ Δ DL ]1,ΔDH1]Then, the normal direction gear of the two electric thrusters is determined to be Cys=Cy1;
If Δ D ∈ [ Δ DL ]2,ΔDH2]Then, the normal direction gear of the two electric thrusters is determined to be Cys=Cy2;
And so on;
if Δ D ∈ [ Δ DL ]n,ΔDHn]Then, the normal direction gear of the two electric thrusters is determined to be Cys=Cyn。
3. The electrical thruster position of claim 2The method for distributing the maintained thrust is characterized in that a target thrust direction [ X ] of the first electric thruster is obtained through calculation according to the determined normal pointing gears of the two electric thrustersn,Yn,Zn]And a boot time period Δ tnAnd a target thrust direction [ X ] of the second electric thrusters,Ys,Zs]And a boot time period Δ tsThe method comprises the following steps:
acquiring eccentricity control vectors delta e in the x direction and the y direction under a track coordinate systemxAnd Δ eyAnd tilt control vectors Δ i in the x-and y-directionsxAnd Δ iy;
According to Δ ex、Δey、ΔixAnd Δ iyAnd calculating to obtain a speed increment VN of the first electric thruster in the x directionxAnd a speed increment VS in the x-direction of the second electric thrusterxAnd the radial distribution identification flagV of the two electric thrusterszAnd the radial distribution difference Vz0;
According to flagVzAnd Vz0And calculating to obtain a speed increment VN of the first electric thruster in the y directionyAnd velocity increment VN in the z-directionzAnd a speed increment VS in the y-direction of the second electric thrusteryAnd z-direction velocity increment VSz;
According to VNy、VNz、VSyAnd VSzAnd calculating to obtain a target thrust direction [ X ] of the first electric thrustern,Yn,Zn]And a boot time period Δ tnAnd a target thrust direction [ X ] of the second electric thrusters,Ys,Zs]And a boot time period Δ ts。
4. The electric thruster position maintaining thrust force distribution method of claim 3, wherein VNxAnd VSxThe calculation flow of (2) is as follows:
according to the nominal orbital velocity V of the geostationary orbitsAnd the earth static orbit nominal orbit angular velocity omegaeAnd calculating to obtain a constant A, a constant B and an intermediate parameter aa:
then there are:
speed increment VN of first electric thruster in x directionxComprises the following steps:
speed increment VS of second electric thruster in x directionxComprises the following steps:
5. the electric thruster position holding thrust force distribution method as claimed in claim 4, wherein the flagV is a flag VzAnd Vz0The calculation flow of (2) is as follows:
calculating to obtain an intermediate variable dVztemp:
According to dVztempDetermining flagVz:
If dVztempIf the value is less than or equal to 0, flagVz=1
If dVztempIf > 0, flagVz=-1
Calculating to obtain a radial scoreMatch difference value Vz0:
Wherein, dVzIndicating the difference, dV, between the radial control quantities of the two electric thrustersz=abs(dVztemp),Vy0And the normal control quantity of the unilateral electric thruster is shown.
6. The electric thruster position maintaining thrust force distribution method of claim 5, wherein VNy、VNz、VSyAnd VSzThe calculation flow of (2) is as follows:
calculating to obtain an intermediate variable Vytemp:
When flagV is presentzWhen 1, there are:
VNy=Vytemp
VNz=Vz0+dVz
VSy=Vy0-Vytemp
VSz=Vz0
when flagV is presentzWhen-1, there is:
VNy=Vy0-Vytemp
VNz=Vz0
VSy=Vytemp
VSz=Vz0+dVz。
7. the electric thruster position holding thrust force distribution method of claim 6, wherein [ X [ ]n,Yn,Zn]、Δtn、[Xs,Ys,Zs]And Δ tsThe calculation process is as followsThe following:
according to VNx、VNyAnd VNzCalculating to obtain [ X ]n,Yn,Zn]And Δ tn:
According to VSx、VSyAnd VSzCalculating to obtain [ X ]s,Ys,Zs]And Δ ts:
Wherein F represents the thrust magnitude of the electric thruster, VeRepresenting the electrical thruster specific impulse, mass representing the satellite mass.
8. The method of claim 7, wherein the method of assigning thrust force to position holding of an electrical thruster is based on [ X ]n,Yn,Zn]And [ X ]s,Ys,Zs]And calculating to obtain the rotating speed omega of the first theoretical vector adjusting mechanismnAnd a second theoretical vector adjusting mechanism rotating speed omegasThe method comprises the following steps:
determining a current heading [ X ] of a first electric thrustern0,Yn0,Zn0]And calculating to obtain [ X ]n,Yn,Zn]And [ X ]n0,Yn0,Zn0]Angle theta therebetweenn;
Determining a current finger of a second electric thrusterTo [ X ]s0,Ys0,Zs0]And calculating to obtain [ X ]s,Ys,Zs]And [ X ]s0,Ys0,Zs0]Angle theta therebetweens;
Obtaining the rotating speed omega of the maximum theoretical vector adjusting mechanismmax;
When theta isn≥θsIn time, there are:
ωn=ωmax
when theta isn<θsIn time, there are:
ωs=ωmax。
9. the method for distributing the position maintaining thrust of the electric thruster as claimed in claim 1, wherein the first vector adjusting mechanism is installed on the south side of the satellite and connected with the first electric thruster for realizing the adjustment of the thrust direction in the three directions of x, y and z; the second vector adjusting mechanism is arranged on the north side of the satellite and connected with the second electric thruster, and is used for adjusting thrust direction in the x direction, the y direction and the z direction.
10. An electric thruster position maintaining thrust distribution system, comprising:
the gear determining module is used for determining normal pointing gears of the two electric thrusters according to the position holding drift rate control quantity;
a first calculation module, configured to calculate, according to the determined normal pointing gears of the two electric thrusters, a target thrust direction [ X ] of the first electric thrustern,Yn,Zn]And a boot time period Δ tnAnd target push of the second electric thrusterDirection of force [ X ]s,Ys,Zs]And a boot time period Δ ts;
A second calculation module for calculating according to [ Xn,Yn,Zn]And [ X ]s,Ys,Zs]And calculating to obtain the rotating speed omega of the first theoretical vector adjusting mechanismnAnd a second theoretical vector adjusting mechanism rotating speed omegas;
An allocation module for allocating the calculated [ Xn,Yn,Zn]And Δ tnDistributing the calculated [ X ] to the first electric thrusters,Ys,Zs]And Δ tsDistributing the calculated omega to a second electric thrusternAssigning the calculated ω to the first vector adjustment mechanismsIs assigned to the second vector adjustment mechanism.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111006812.3A CN113734469B (en) | 2021-08-30 | 2021-08-30 | Method and system for distributing thrust for maintaining position of electric thruster |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111006812.3A CN113734469B (en) | 2021-08-30 | 2021-08-30 | Method and system for distributing thrust for maintaining position of electric thruster |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113734469A true CN113734469A (en) | 2021-12-03 |
CN113734469B CN113734469B (en) | 2023-02-03 |
Family
ID=78734024
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111006812.3A Active CN113734469B (en) | 2021-08-30 | 2021-08-30 | Method and system for distributing thrust for maintaining position of electric thruster |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113734469B (en) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6032904A (en) * | 1998-02-23 | 2000-03-07 | Space Systems/Loral, Inc. | Multiple usage thruster mounting configuration |
US6565043B1 (en) * | 2001-12-21 | 2003-05-20 | The Boeing Company | Redundant system for satellite inclination control with electric thrusters |
CN104477413A (en) * | 2014-11-06 | 2015-04-01 | 中国空间技术研究院 | Combined control method for performing station keeping and momentum wheel unloading based on electric thruster |
CN105353621A (en) * | 2015-11-30 | 2016-02-24 | 北京控制工程研究所 | Fault mode thrust allocation method for geostationary orbit satellite electric thruster |
CN105511490A (en) * | 2015-12-15 | 2016-04-20 | 北京理工大学 | Geostationary orbit satellite station keeping-angular momentum unloading combined control method |
CN106184818A (en) * | 2016-03-07 | 2016-12-07 | 北京理工大学 | Geostationary satellite Actuator failure operating mode lower railway position keeping method |
CN106394935A (en) * | 2016-10-31 | 2017-02-15 | 北京控制工程研究所 | Electric propulsion angular momentum unloading method taking thruster arc section losses into consideration |
CN107891997A (en) * | 2017-10-11 | 2018-04-10 | 北京控制工程研究所 | Taper layout electric propulsion satellite failure mode position keeps optimal thrust distribution method |
CN108490963A (en) * | 2018-02-08 | 2018-09-04 | 中国空间技术研究院 | Position keeping method under full electric propulsion satellite electric thruster fault mode and system |
CN108516106A (en) * | 2018-02-08 | 2018-09-11 | 中国空间技术研究院 | A kind of full electric propulsion Satellite Orbit Maneuver process angular momentum dumping method and system |
CN110254753A (en) * | 2019-06-04 | 2019-09-20 | 北京理工大学 | A kind of geostationary orbit satellite electric thruster and its layout optimization method |
-
2021
- 2021-08-30 CN CN202111006812.3A patent/CN113734469B/en active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6032904A (en) * | 1998-02-23 | 2000-03-07 | Space Systems/Loral, Inc. | Multiple usage thruster mounting configuration |
US6565043B1 (en) * | 2001-12-21 | 2003-05-20 | The Boeing Company | Redundant system for satellite inclination control with electric thrusters |
CN104477413A (en) * | 2014-11-06 | 2015-04-01 | 中国空间技术研究院 | Combined control method for performing station keeping and momentum wheel unloading based on electric thruster |
CN105353621A (en) * | 2015-11-30 | 2016-02-24 | 北京控制工程研究所 | Fault mode thrust allocation method for geostationary orbit satellite electric thruster |
CN105511490A (en) * | 2015-12-15 | 2016-04-20 | 北京理工大学 | Geostationary orbit satellite station keeping-angular momentum unloading combined control method |
CN106184818A (en) * | 2016-03-07 | 2016-12-07 | 北京理工大学 | Geostationary satellite Actuator failure operating mode lower railway position keeping method |
CN106394935A (en) * | 2016-10-31 | 2017-02-15 | 北京控制工程研究所 | Electric propulsion angular momentum unloading method taking thruster arc section losses into consideration |
CN107891997A (en) * | 2017-10-11 | 2018-04-10 | 北京控制工程研究所 | Taper layout electric propulsion satellite failure mode position keeps optimal thrust distribution method |
CN108490963A (en) * | 2018-02-08 | 2018-09-04 | 中国空间技术研究院 | Position keeping method under full electric propulsion satellite electric thruster fault mode and system |
CN108516106A (en) * | 2018-02-08 | 2018-09-11 | 中国空间技术研究院 | A kind of full electric propulsion Satellite Orbit Maneuver process angular momentum dumping method and system |
CN110254753A (en) * | 2019-06-04 | 2019-09-20 | 北京理工大学 | A kind of geostationary orbit satellite electric thruster and its layout optimization method |
Also Published As
Publication number | Publication date |
---|---|
CN113734469B (en) | 2023-02-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104015938B (en) | A kind of position keeping method of electric propulsion satellite | |
CN109552670B (en) | Application of low-thrust control in geostationary satellite orbit dip angle maintenance | |
CN110435930B (en) | Low-orbit optical satellite uniform-deceleration push-broom attitude planning method | |
EP3854698B1 (en) | Orientation control device, satellite, orientation control method, and program | |
US9764859B2 (en) | Efficient stationkeeping design for mixed fuel systems | |
CN109911249A (en) | The interspace transfer Finite Thrust of low thrust ratio aircraft is entered the orbit interative guidance method | |
CN108454886A (en) | A kind of electric propulsion system milli ox grade thrust on-orbit calibration method | |
CN110697086B (en) | Method for realizing satellite three-axis stable control by single control moment gyroscope and two flywheels | |
CN104729457A (en) | Method for determining position of sun relative to near-earth microsatellites | |
CN110096726B (en) | GEO satellite emergency transfer orbit rapid optimization design method based on lunar force borrowing | |
KR20170002287A (en) | Efficient stationkeeping design for mixed fuel systems in response to a failure of an electric thruster | |
JPH04349098A (en) | Method of controlling motion of the east direction of geostationary satellite | |
CN112607065A (en) | High-precision phase control method based on electric propulsion system | |
CN112130590B (en) | Satellite-borne antenna ground pointing determination method based on speed compensation under instantaneous inertial system | |
CN103955138B (en) | A kind of dynamic middle imaging satellite attitude control method based on increment type drift angle | |
CN111268177A (en) | Distributed closed-loop autonomous position maintaining control method for geostationary orbit satellite | |
CN113641182A (en) | High-precision aiming and pointing method and system for inter-satellite laser communication system | |
CN113734469B (en) | Method and system for distributing thrust for maintaining position of electric thruster | |
CN114889849A (en) | Estimation method for multi-constraint flying-over optimal speed of spacecraft | |
CN112945242B (en) | Method for autonomously planning optimal time and attitude of task on orbit by satellite | |
CN109190155B (en) | Hybrid continuous low-thrust track design method adopting electric propulsion/solar sail propulsion | |
CN111605735B (en) | Spacecraft three-phase control variable envelope angular momentum analysis method | |
US6860451B1 (en) | Spacecraft spin axis reorientation method | |
JP7068765B2 (en) | Satellite constellation formation system, satellite constellation formation method, satellite constellation formation program, and ground equipment | |
JP2004210032A (en) | Formation flying satellite |
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 |