CN106542120B - In conjunction with the satellite three-axis attitude control method of magnetic torquer when flywheel drive lacking - Google Patents

In conjunction with the satellite three-axis attitude control method of magnetic torquer when flywheel drive lacking Download PDF

Info

Publication number
CN106542120B
CN106542120B CN201610877412.2A CN201610877412A CN106542120B CN 106542120 B CN106542120 B CN 106542120B CN 201610877412 A CN201610877412 A CN 201610877412A CN 106542120 B CN106542120 B CN 106542120B
Authority
CN
China
Prior art keywords
magnetic
flywheel
axis
control
torquer
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
CN201610877412.2A
Other languages
Chinese (zh)
Other versions
CN106542120A (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.)
Shanghai Aerospace Control Technology Institute
Original Assignee
Shanghai Aerospace Control Technology Institute
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 Shanghai Aerospace Control Technology Institute filed Critical Shanghai Aerospace Control Technology Institute
Priority to CN201610877412.2A priority Critical patent/CN106542120B/en
Publication of CN106542120A publication Critical patent/CN106542120A/en
Application granted granted Critical
Publication of CN106542120B publication Critical patent/CN106542120B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/32Guiding or controlling apparatus, e.g. for attitude control using earth's magnetic field
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/244Spacecraft control systems
    • B64G1/245Attitude control algorithms for spacecraft attitude control

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

It is in conjunction with the satellite three-axis attitude control method of magnetic torquer, step when flywheel drive lacking:(1) the geomagnetic field intensity vector B obtained according to magnetometer survey determines magnetic control matrix Γ (b), and can use situation according to actual flywheel, determines that flywheel controls matrix Kw;(2) three-axis attitude control magnetic current-order I is calculatedzk=Kic*Mzk;(3) three axis magnetic dumpings magnetic current-order I is calculatedxz=Kic*Mxz;(4) it calculates magnetic torquer and instructs Ick=Izk+Ixz;While magnetic torquer instructs progress gesture stability according to magnetic torquer, available counteraction flyback continues torque T based on instructioncCarry out gesture stability.The method of the present invention is when system is only left arbitrary two or even a flywheel is available, satellite three-axis attitude high-precision control is realized in conjunction with three-axis magnetorquer, suitable for the control of long-term absolute orientation or long-term Direct to the sun control, it is also applied for satellite any attitude Large angle maneuver control.

Description

In conjunction with the satellite three-axis attitude control method of magnetic torquer when flywheel drive lacking
Technical field
The invention belongs to Spacecraft Attitude Control fields, are related to satellite three-axis attitude controlling party when a kind of flywheel drive lacking Method.
Background technology
Modern Small Satellites do not carry gas ejecting system generally, counteraction flyback become moonlet high-performance pose stabilization control and The main executing agency of large angle attitude control.When the average value for the disturbance torque that satellite is subject to is not zero, flywheel Angular momentum can accumulate at any time, eventually lead to Speed of Reaction Wheels saturation.Therefore Modern Small Satellites also need configuration magnetic torquer to realize winged Wheel unloading.I.e. generally use counteraction flyback carries out gesture stability, using magnetic torquer cooperation flywheel work, completes flywheel and unloads It carries.
Previous satellite attitude control system, to meet lifetime of system and reliability requirement, usually configuration four or four Above counteraction flyback is designed by each flywheel installation direction, and when any one flywheel failure, system still has three or more Flywheel is available, remains to complete three-axis attitude (rolling, pitching, yaw-position) control task.But for Modern Small Satellites attitude control system System, to meet power consumption, quality and the limitation of cost, usually only only a few critical component realizes hardware redundancy, needs by event Software processing means under barrier pattern improve the reliability of system, extend task service life.For this purpose, its posture control system is usually only matched Three counteraction flybacks are set, are mounted on celestial body main shaft with orthogonal manner, no hardware backup, as shown in Figure 1.When any one or When arbitrary two flywheels break down, fly wheel system will be unable to provide three-axis attitude control moment, therefore attitude control system is moved back Typical flywheel underactuated control system (control input number is less than degree of freedom in system number) is turned to, no longer has three-axis attitude Control ability.Only it is 2 dimensions, to be extended to 3 in addition, since magnetic control torque is always in the plane of vertical magnetic field unit vector b Control moment is tieed up, flywheel or thruster need to be introduced.Therefore in the case of fly wheel system drive lacking, it must combine residue that can use flywheel Jointly controlled with magnetic torquer, could realize that three-axis attitude controls.
Invention content
Present invention solves the technical problem that being:It is combined when having overcome the deficiencies of the prior art and provide a kind of flywheel drive lacking The satellite three-axis attitude control method of magnetic torquer is remained to when system is only left arbitrary two or even a flywheel is available The high-precision control of satellite three-axis attitude is realized in conjunction with three-axis magnetorquer.
Technical solution of the invention is:In conjunction with the satellite three-axis attitude controlling party of magnetic torquer when flywheel drive lacking Method includes the following steps:
(1) the geomagnetic field intensity vector B obtained according to magnetometer survey determines magnetic control matrix Γ (b), and according to actual Flywheel can use situation, determine that flywheel controls matrix Kw,
Wherein,When i axis flywheels are available, ki=1, it is no Then when i axis flywheel failures, ki=0, i=x, y, z;
(2) three-axis attitude control magnetic current-order I is calculatedzk=Kic*Mzk, wherein kicFor magnetic torquer magnetic moment and electric current Conversion coefficient, Mzk=B × Λ-1(b)Tc/||B||2, Λ (b)=Γ (b)+Kw, TcThree axis provided for flywheel gesture stability algorithm Control moment instructs;
(3) three axis magnetic dumpings magnetic current-order I is calculatedxz=Kic*Mxz, wherein Mxz=-Ku(B×Δh)/||B||2, Δ H is the difference of current flywheel angular momentum and unloading target angular momentum, KuFor magnetic dumping coefficient matrix
Kux、Kuy、KuzRespectively three axis magnetic dumping coefficient of satellite;
(4) it calculates magnetic torquer and instructs Ick=Izk+Ixz, IckIt is three axis magnetic force away from device magnetic controlled current;Magnetic torquer according to While magnetic torquer instruction carries out gesture stability, available counteraction flyback continues torque T based on instructioncCarry out posture control System.
Magnetic torquer is instructed according to magnetic torquer maximum drive current range in the step (4) and is carried out at amplitude limit Reason, namely
Wherein ImaxFor separate unit magnetic torquer maximum operating currenbt.
The advantages of the present invention over the prior art are that:The prior art is matched by counteraction flyback standalone hardware redundancy It sets, to improve the reliability of system.The present invention can carry out satellite three-axis attitude control by residue with flywheel combination magnetic torquer, I.e. in the case of no flywheel hardware backup, the reliability of system is improved by software processing means.The method of the present invention with fly Attitude control method when wheel is all normal is compared, it is only necessary to change (the wheel control matrix K of constant parameter at onew), without control The switching of molding formula and control algolithm is easy to Project Realization because the method is succinct.Correspondingly, the method for the present invention is normal in flywheel When with flywheel exception, control mode switch need not be carried out, it is only necessary to the available flag of corresponding flywheel is provided in flywheel exception ?;This method is equally applicable to the case where counteraction flyback all works normally, and the geomagnetic torque generated at this time will be with flywheel Based on unloading effect, without being had an impact to gesture stability.
Description of the drawings
Fig. 1 is that typical flywheel installs configuration schematic diagram;
Fig. 2 is that typical magnetic torquer installs configuration schematic diagram;
Fig. 3 is the flow diagram of the method for the present invention.
Specific implementation mode
The method of the present invention can use situation and current geomagnetic fieldvector according to actual flywheel, by defining magnetic control matrix Γ (b) and wheel controls matrix Kw, it is calculated and instructs magnetic moment for three axis of gesture stability.By the magnetic current-order for gesture stability With the magnetic current-order summation for magnetic dumping, you can obtain final magnetoelectricity flow control instructions, at the same time, available flywheel Continue to execute the control instruction that respective wheel control algorithm provides.
Specifically by taking the installation configuration of the common flywheel and magnetic torquer of attached drawing 1 and attached drawing 2 as an example, system configures three flywheels altogether With three magnetic torquers, three flywheels are respectively along the orthogonal installation of tri- main shafts of body coordinate system ObXbYbZb, three magnetic torquers It is respectively parallel to the orthogonal installation of tri- main shafts of body coordinate system ObXbYbZb.ObXb corresponds to the axis of rolling, ObYb corresponds to pitch axis, ObZb corresponds to yaw axis.
As shown in figure 3, for the flow chart of the method for the present invention.The specific steps that the present invention uses:The first step is according to earth's magnetic field Vector operation magnetic control matrix Γ (b);And situation, setting wheel control matrix K can be used according to actual flywheeli;Second step combination posture Control instruction torque TcGesture stability magnetic current-order I is obtained with magnetic control matrix Γ (b)zk;Third step design flywheel magnetic dumping is calculated Method obtains unloading magnetic current-order Ixz;4th step passes through magnetic dumping current-order IxzWith gesture stability magnetic current-order IzkIt obtains Final magnetic torquer output (magnetoelectricity stream) instructs Ick, while magnetic torquer carries out gesture stability, available reaction flies Wheel continues to execute gesture stability instruction torque Tc.Specifically:
(1) magnetic control matrix is calculated
The unit vector b of geomagnetic fieldvector is calculated first
Wherein B is the geomagnetic field intensity vector that magnetometer survey obtains, the component B=at satellite body system ObXbYbZb [Bx By Bz], unit tesla;
Calculate magnetic control matrix Γ (b)
Wherein b×Indicate the multiplication cross matrix of unit of magnetic field strength vector
According to flywheel failure (available) situation, determine that flywheel controls matrix Kw
Wherein kx、ky、kzThe respectively available flag of Xb axis, Yb axis, Zb axis flywheels;When i (i=x or y or z) axis flywheel can Used time, ki=1, otherwise when i axis flywheel failures, take ki=0;Specific flywheel can use situation, can be examined according to in-orbit flywheel failure Note result determines in disconnected result or ground.
Magnetic control matrix characterizes the relationship between gesture stability instruction torque and magnetic torquer actual output torque;Wheel control square Battle array characterizes the available situation of each main shaft flywheel, can adapt to the case where only remaining flywheel of system can be used.
By taking the axis of rolling (one ObXb axis of attached drawing) flywheel failure as an example, flywheel controls matrix KwFor:
It should be noted that when any one or two flywheels of system spare are available, algorithm is still valid;But when system flies When wheel all fails, algorithm is unavailable.
(2) three-axis attitude control magnetic current-order is calculated
First calculate combination control matrix Λ (b), combination control matrix Λ (b) characterize flywheel and magnetic torquer composition Executive capability of the execution system to instruction torque.
Λ (b)=Γ (b)+Kw
The instruction magnetic moment M for gesture stability can be calculated in turnzk
Mzk=B × Λ-1(b)Tc/||B||2
Wherein TcFor the three axis control moments instruction that general flywheel gesture stability algorithm provides, WIE B are specifically referred to, WEISS H,ARAPOSTATHIS A.Quaternion feedback regulator for spacecraft eigenaxis rotations[J].Journal of Guidance,Control and Dynamics,1989,12(3):375-380。
And then instruction current I of the three-axis magnetorquer for gesture stability can be calculatedzk
Izk=Kic*Mzk
Wherein kicFor magnetic torquer magnetic moment and electric current conversion coefficient.Maximum functional electricity is corresponded to maximum output magnetic moment 20Am2 For the magnetic torquer for flowing 80mA, Kic=80/20.
Introduce the inverse of combinatorial matrix (magnetic control matrix+wheel control matrix), it is ensured that the reality of flywheel+magnetic torquer combination is defeated Gesture stability instruction torque can be tracked to greatest extent by going out torque.
(3) three axis flywheels unloading magnetic current-order is calculated
Magnetic dumping control algolithm can realize the momentum management to counteraction flyback, pass through celestial body angular momentum exchange so that Flywheel is operated near desired center rotating speed.
Mxz=-Ku(B×Δh)/||B||2
Wherein MxzIt is current flywheel angular momentum and this unloading target angular momentum to instruct magnetic moment, Δ h for three axis magnetic dumpings Difference, when target angular momentum is 0, Δ h is the current angular momentum Δ h=[h of three axis flywheelsx hy hz];KuFor magnetic dumping coefficient Matrix, when magnetic field unit is Gs, generally desirable Kux=Kuy=Kuz=10;
And then instruction current I of the three-axis magnetorquer for flywheel unloading can be calculatedxz
Ixz=Kic*Mxz
(4) magnetic torquer instruction (magnetoelectricity stream) is calculated
Ick=Izk+Ixz
IckIt is three axis magnetic force away from device magnetic controlled current, unit mA;When practical application, according to magnetic torquer maximum drive current model It encloses and amplitude limiting processing is carried out to instruction magnetoelectricity stream.
Wherein ImaxFor separate unit magnetic torquer maximum operating currenbt.
I.e. final magnetic torquer output (magnetoelectricity stream) instructs IckDirectly by gesture stability magnetic current-order IzkAnd magnetic dumping Current-order IxzSummation obtains.
While magnetic torquer carries out gesture stability, available counteraction flyback continues torque T based on instructioncCarry out appearance State controls.
In the case where counteraction flyback all works normally, this method is still valid;The earth magnetism that magnetic torquer generates at this time Based on torque will be acted on flywheel unloading, without being significantly affected to gesture stability generation.
It should be noted that magnetic torquer gesture stability algorithm is in the case where geomagnetic fieldvector is parallel with flywheel failure axis It will appear unusual, attitude control error caused to become larger;But inclined plane satellite is not deposited during satellite absolute orientation At the geomagnetic fieldvector situation parallel with failure axis (two axis of no magnetic field simultaneously zero passage), therefore flywheel+magnetic control of the present invention Method is suitable for the long-term three axis absolute orientation control of inclined plane satellite or long-term Direct to the sun control, while in motor-driven angular speed In the case of of less demanding, it is also applied for satellite any attitude Large angle maneuver control.
The content that description in the present invention is not described in detail belongs to the known technology of those skilled in the art.

Claims (2)

1. in conjunction with the satellite three-axis attitude control method of magnetic torquer when flywheel drive lacking, it is characterised in that include the following steps:
(1) the geomagnetic field intensity vector B obtained according to magnetometer survey determines magnetic control matrix Γ (b), and according to actual flywheel Situation can be used, determines that flywheel controls matrix Kw,
Wherein,When i axis flywheels are available, ki=1, otherwise when When i axis flywheel failures, ki=0, i=x, y, z;
(2) three-axis attitude control magnetic current-order I is calculatedzk=Kic*Mzk, wherein KicIt is converted for magnetic torquer magnetic moment and electric current Coefficient, Mzk=B × Λ-1(b)Tc/||B||2, Λ (b)=Γ (b)+Kw, TcThe three axis control provided for flywheel gesture stability algorithm Torque command;
(3) three axis magnetic dumpings magnetic current-order I is calculatedxz=Kic*Mxz, wherein Mxz=-Ku(B×Δh)/||B||2, Δ h is to work as The difference of preceding flywheel angular momentum and unloading target angular momentum, KuFor magnetic dumping coefficient matrix
Kux、Kuy、KuzRespectively three axis magnetic dumping coefficient of satellite;
(4) magnetic torquer magnetic controlled current I is calculatedck=Izk+Ixz, in magnetic torquer posture is carried out according to magnetic torquer magnetic controlled current While control, available counteraction flyback continues to be based on torque command TcCarry out gesture stability.
2. special in conjunction with the satellite three-axis attitude control method of magnetic torquer when flywheel drive lacking according to claim 1 Sign is:Amplitude limit is carried out to magnetic torquer magnetic controlled current according to magnetic torquer maximum drive current range in the step (4) Processing, namely
Wherein ImaxFor separate unit magnetic torquer maximum drive current.
CN201610877412.2A 2016-09-30 2016-09-30 In conjunction with the satellite three-axis attitude control method of magnetic torquer when flywheel drive lacking Active CN106542120B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610877412.2A CN106542120B (en) 2016-09-30 2016-09-30 In conjunction with the satellite three-axis attitude control method of magnetic torquer when flywheel drive lacking

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610877412.2A CN106542120B (en) 2016-09-30 2016-09-30 In conjunction with the satellite three-axis attitude control method of magnetic torquer when flywheel drive lacking

Publications (2)

Publication Number Publication Date
CN106542120A CN106542120A (en) 2017-03-29
CN106542120B true CN106542120B (en) 2018-11-02

Family

ID=58368239

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610877412.2A Active CN106542120B (en) 2016-09-30 2016-09-30 In conjunction with the satellite three-axis attitude control method of magnetic torquer when flywheel drive lacking

Country Status (1)

Country Link
CN (1) CN106542120B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3702289A1 (en) * 2019-03-01 2020-09-02 Centre National d'Etudes Spatiales Method for controlling the attitude of a satellite in low orbit

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109649690B (en) * 2018-12-27 2021-07-09 上海航天控制技术研究所 High-precision hierarchical saturation attitude control method and system for satellite based on flywheel configuration
CN109649693B (en) * 2019-01-21 2021-12-14 上海微小卫星工程中心 Pure magnetic control spinning sun-facing orientation method
CN109975880B (en) * 2019-04-04 2021-07-27 哈尔滨工业大学 Feature vector-based orientation method, device and system
CN110543181B (en) * 2019-08-26 2022-08-30 北京电子工程总体研究所 Under-actuated angle decoupling attitude control method and system
CN111762342B (en) * 2020-04-24 2022-04-12 北京控制工程研究所 Optimization method for under-actuated configuration of reaction wheel
CN112429277B (en) * 2020-11-25 2022-07-08 山东航天电子技术研究所 Method for restraining back electromotive force of magnetic torquer
CN112947116B (en) * 2021-02-09 2022-07-26 南京航空航天大学 Control method for joint execution mechanism of three-axis air bearing table
CN112937920A (en) * 2021-03-30 2021-06-11 湖南揽月机电科技有限公司 Multi-redundancy satellite intelligent attitude control assembly and working method thereof
CN112977890B (en) * 2021-04-06 2024-01-16 湖南揽月机电科技有限公司 Coil magnetic torquer in intelligent attitude control assembly
CN113335567B (en) * 2021-05-26 2022-09-23 航天科工空间工程发展有限公司 Wheel magnetic hybrid attitude control method and system for microsatellite
CN114476135B (en) * 2022-02-22 2024-02-09 中国空间技术研究院 Satellite bias momentum flywheel hot standby method

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3606636C1 (en) * 1986-02-28 1987-11-05 Messerschmitt Boelkow Blohm Method for determining geomagnetic field components with reference to a satellite-fixed coordinate system
CN102880183B (en) * 2012-10-15 2015-04-22 北京控制工程研究所 Control moment gyro angular momentum management method for yawing maneuvering
CN104176275B (en) * 2014-07-16 2016-05-04 北京控制工程研究所 A kind of rate damping method that uses momenttum wheel to combine with magnetic torquer
CN105644810B (en) * 2016-02-05 2017-09-12 上海微小卫星工程中心 A kind of opened loop control reaction wheel rotation speed change and discharging method

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3702289A1 (en) * 2019-03-01 2020-09-02 Centre National d'Etudes Spatiales Method for controlling the attitude of a satellite in low orbit
FR3093326A1 (en) * 2019-03-01 2020-09-04 Centre National d'Études Spatiales Method of controlling the attitude of a satellite in low orbit

Also Published As

Publication number Publication date
CN106542120A (en) 2017-03-29

Similar Documents

Publication Publication Date Title
CN106542120B (en) In conjunction with the satellite three-axis attitude control method of magnetic torquer when flywheel drive lacking
CN106660641B (en) Method for controlling the orbit of a satellite in earth orbit, satellite and system for controlling the orbit of such a satellite
KR101008176B1 (en) Maneuverability and controllability improvement using reaction wheel-based and thruster-based attitude controller simultaneously
CN109533396B (en) Satellite spin orientation method based on magnetic measurement and control
CN108069050B (en) Spacecraft initial attitude capture control method and system
Huang et al. Adaptive reaching law based sliding mode control for electromagnetic formation flight with input saturation
US9650159B2 (en) Method and device for electric satellite propulsion
CN107600463B (en) A kind of agility small satellite attitude fast reserve control method
Wu et al. Null motion strategy for spacecraft large angle agile maneuvering using hybrid actuators
Ye et al. Hybrid thrusters and reaction wheels strategy for large angle rapid reorientation with high precision
CN106094853B (en) A kind of control method of Upper Stage Orbit Transformation section vectored thrust
CN110697085B (en) Satellite control method combining double SGCMG and magnetic torquer
Jan et al. Attitude control system for ROCSAT-3 microsatellite: a conceptual design
Gui et al. Maneuver planning of a rigid spacecraft with two skew control moment gyros
Ao et al. Touchless attitude correction for satellite with constant magnetic moment
Mohammed et al. Magnetorquer control for orbital manoeuvre of low earth orbit microsatellite
US10934026B2 (en) Propulsion system with differential throttling of electric thrusters
CN206155823U (en) Be used for propulsive device of robot for space
JP4930072B2 (en) Satellite attitude control device
Kristiansen et al. A comparative study of actuator configurations for satellite attitude control
CN113247310B (en) Estimation method and system suitable for continuous attitude maneuver times of satellite
Li et al. Attitude control of staring-imaging satellite using Permanent Magnet momentum Exchange Sphere
Maclean et al. Computationally light attitude controls for resource limited nano-spacecraft
Jordaan et al. The attitude control of a tri-spin solar sail satellite
Kou et al. A single-axis pointing magnetic control method for microsatellite

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