EP3643621B1 - Satellite attitude control system using eigen vector, non-linear dynamic inversion, and feedforward control - Google Patents

Satellite attitude control system using eigen vector, non-linear dynamic inversion, and feedforward control Download PDF

Info

Publication number
EP3643621B1
EP3643621B1 EP19203064.1A EP19203064A EP3643621B1 EP 3643621 B1 EP3643621 B1 EP 3643621B1 EP 19203064 A EP19203064 A EP 19203064A EP 3643621 B1 EP3643621 B1 EP 3643621B1
Authority
EP
European Patent Office
Prior art keywords
satellite
orientation
reaction wheel
control system
eigen vector
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
EP19203064.1A
Other languages
German (de)
French (fr)
Other versions
EP3643621A1 (en
Inventor
John Benton Derrick II
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.)
General Atomics Corp
Original Assignee
General Atomics Corp
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 General Atomics Corp filed Critical General Atomics Corp
Priority to SI201930327T priority Critical patent/SI3643621T1/en
Publication of EP3643621A1 publication Critical patent/EP3643621A1/en
Application granted granted Critical
Publication of EP3643621B1 publication Critical patent/EP3643621B1/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/08Control of attitude, i.e. control of roll, pitch, or yaw
    • G05D1/0808Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
    • 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/28Guiding or controlling apparatus, e.g. for attitude control using inertia or gyro effect
    • B64G1/283Guiding or controlling apparatus, e.g. for attitude control using inertia or gyro effect using reaction wheels
    • 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
    • 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

Definitions

  • the present invention relates generally to orienting a satellite, and more specifically to orienting a satellite using eigen vector rotation and non-linear dynamic inversion.
  • Satellites may be controlled by a system known as an attitude determination and control system (ADACS).
  • ADACS attitude determination and control system
  • an ADACS system may control the attitude, or orientation, of a satellite using reaction wheels that turn according to software algorithms executed on a computer within the satellite.
  • a satellite reaction wheel system may be characterized by dynamics that are highly non-linear.
  • Some control algorithms rely on linear control techniques such as a proportional-integral-derivative (PID) control to control the attitude of the satellite.
  • PID proportional-integral-derivative
  • Linear control systems applied to non-linear systems can exhibit oscillations, overshoots, and even instability.
  • ADACS systems respond to pointing commands in a manner that resembles a yaw, pitch, and roll sequence. This method of maneuvering may be inefficient in terms of energy use and time.
  • WO 2017/159156 A1 describes a method and system for spacecraft orientation control using an inner-loop control determining first control inputs for momentum exchange devices to control an orientation of the spacecraft and an outer-loop control determining second control inputs for thrusters of the spacecraft to concurrently control a pose of the spacecraft and a momentum stored by the momentum exchange devices of the spacecraft.
  • the outer-loop control determines the second control inputs using a model of dynamics of the spacecraft including dynamics of the inner-loop control, such that the outer-loop control accounts for effects of actuation of the momentum exchange devices according to the first control inputs determined by the inner-loop control.
  • the thrusters and the momentum exchange devices are controlled according the first and second control inputs.
  • US 9 745 082 B2 describes an attitude control system and method for a satellite based on four single degree-of-freedom control moment gyroscopes with variable speed flywheels (or reaction wheels) in a pyramid configuration, combined with path and endpoint constraint time-optimal control.
  • EP 0 926 066 A1 describes an attitude control system and method of a spacecraft in a state control device of a moving body constituted by a navigation dynamics, an actuator for driving the navigation dynamics, first controlling means for controlling the actuator in PID control in response to a first output signal outputted from the navigation dynamics and adding means for outputting a control signal for controlling the actuator in feedforward control in response to outside noise by adding an estimated value of the outside noise to a control signal outputted from the first controlling means.
  • the method of satellite orientation control includes applying with a processing device a satellite orientation control system, the satellite orientation control system comprising a double feedback loop system, wherein an outer loop receives a desired orientation of the satellite and an estimated orientation of the satellite as inputs, determines an eigen vector to rotate the satellite from one orientation to another based on the desired orientation and the estimated orientation and outputs the determined eigen vector or desired body rotational rates determined based on the determined eigen vector, and an inner loop receives the eigen vector or the desired body rotational rates from the outer loop as an input and executes a non-linear dynamic inversion algorithm based on the eigen vector or the desired body rotational rates to output an output signal to at least one reaction wheel of the satellite, rotating the at least one reaction wheel in response to the output signal, and orienting the satellite based upon the rotation of the at least one reaction wheel.
  • the system comprises a satellite including a pointing command generator, a navigation system and at least one reaction wheel and includes a satellite orientation control system communicatively connected to the pointing command generator, the navigation system and the at least one reaction wheel, the satellite orientation control system comprising a double feedback loop system including an outer loop and an inner loop, wherein the outer loop is configured to receive a desired orientation of the satellite from the pointing command generator and an estimated orientation of the satellite from the navigation system as inputs, determine an eigen vector to rotate the satellite from one orientation to another based on the desired orientation and the estimated orientation and output the determined eigen vector or desired body rotational rates determined based on the determined eigen vector to the inner loop, and wherein the inner loop is configured to receive the eigen vector or the desired body rotational rates from the outer loop as an input and execute a non-linear dynamic inversion algorithm based on the eigen vector or the desired body rotation
  • the outer loop further comprises executing a satellite orientation error command based on the estimated orientation of the satellite and the desired orientation of the satellite.
  • the eigen vector command executes to decompose the satellite orientation error command into a scalar component and a vector component.
  • the outer loop further comprises a feed-forward control system.
  • the feed-forward control system controls for timing errors between rotation of the at least one reaction wheel and pointing at a target.
  • the feed-forward control system accounts for motion of tracking dishes, antennae, cameras, robotic arms, and solar arrays associated with the satellite. In some examples of the method and system described above, the feed-forward control system further comprises receiving the desired orientation as an input.
  • the inner loop further comprises receiving measured satellite rotation rates from a navigation system associated with the satellite.
  • the inner loop further comprises receiving a measured reaction wheel speed of the at least one reaction wheel from a rotation wheel tachometer associated with the satellite. In some examples of the method and system described above, the inner loop further comprises determining a desired rotational acceleration for the at least one reaction wheel. In some examples of the method and system described above, the inner loop further comprises receiving a mass moment of inertia tensor of at least one reaction wheel, a rotation axis vector of at least one reaction wheel, and a mass moment of inertia tensor of the satellite.
  • the present disclosure provides a satellite control system that exhibits improved stability and increased efficiency by implementing a non-linear dynamic inversion inner-loop control algorithm coupled with an eigen vector outer-loop control algorithm.
  • the attitude determination and control system (ADACS) system may operate using commands to rotate directly about an eigen vector (i.e., to go from one orientation directly to another).
  • the outer-loop control system is augmented with a feed-forward control element to enhance pointing accuracy when tracking moving targets.
  • FIG. 1 shows an example of a satellite 100 in accordance with aspects of the present disclosure.
  • Satellite 100 includes a pointing command generator 105, navigation system 110, reaction wheels 115, a reaction wheel tachometer 120, and an attitude control system 125 (which may also be referred to as a satellite control system).
  • the attitude control system 125 may provide for stability and robustness over a wide range of operational modes and on-orbit conditions; smooth, quick, and accurate control system response to attitude commands; energy efficient operation to maximize battery life; ease of control system 125 setup; and ability to expand the control system's 125 use to large and complex satellites 100.
  • the system 125 achieves these results by implementing an eigen vector outer-loop control algorithm, combined with a non-linear dynamic inversion inner-loop control algorithm. These results represent an improvement on existing satellite 100 control technology, but they are not a comprehensive list of features or advantages of the present system.
  • the outer loop of the attitude control system 125 utilizes an eigen vector outer-loop control algorithm. That is, between any two orientations of the satellite 100 attitude, there exists a rotation axis, called the eigen vector, which leads the satellite 100 directly from one orientation to the other. By rotating along this axis, the time and energy needed to complete the attitude maneuver may be minimized.
  • the inner loop of the attitude control system 125 utilizes non-linear dynamic inversion control.
  • the system accounts for stored momentum in the reaction wheels 115 and the satellite 100 which produces a smooth, uncoupled, linear response on a highly non-linear system.
  • the reaction wheels 115 for the satellite 100 may have a top speed of 10,000 RPM.
  • the satellite 100 may become uncontrollable once the reaction wheels 115 have 4,000 RPM of stored momentum due to the non-linear coupling of the satellite 100 control axes.
  • the dynamic inversion control for the attitude control system 125 may be designed to work on satellites 100 with multiple reaction wheels 115 (e.g., more than 3) and with the reaction wheel spin axes orientated in any direction.
  • This aspect of the control system 125 may also be used on advanced satellites 100 with redundant and skewed-axis reaction wheels 115.
  • the system disclosed herein also takes into account cross coupling between the satellite 100 control axes and off-axis reaction wheels 115.
  • the dynamic inversion control for the attitude control system 125 may operate based on the following parameters: the number of reaction wheels 115, the mass moment of inertia of each reaction wheel, the mass moment of inertia tensor of the satellite 100 (including the reaction wheels 115), the coordinates of the unit vectors pointing along each reaction wheel's spin axis, the torque and speed limit of each reaction wheel, the satellite 100 attitude slew rate limit, and the closed-loop bandwidth of the rotational rate and attitude control loops.
  • a control system that only operates on proportional error feedback will result in the satellite 100 attitude lagging the desired orientation when the command is changing.
  • Changing commands occur when the satellite 100 is attempting to track a moving target such as a point on the Earth's surface or to be aligned with the nadir orientation.
  • the lagging response can result in significant steady-state pointing errors (e.g., approximately 10 degrees).
  • Feed-forward control can accomplish the goal of near zero steady-state error while not inducing unwanted oscillations and overshoot. The only trade-off is that these algorithms are somewhat more difficult to develop.
  • feed-forward algorithm is designed for pointing vector control (e.g. earth position pointing)
  • a second feed-forward algorithm is designed for fully constrained attitude commands such as nadir pointing. It is desired to have a feed-forward control system that will perfectly track the attitude commands in the absence of modeling errors and disturbances.
  • the vector pointing feed-forward control algorithm tracks the rate of change of the pointing vector command by using cross products and dot product derivatives.
  • the output of this algorithm is a body-frame rotational rate command that is sent directly into the inner-loop rotational rate control system. Simulation studies demonstrate that the Earth pointing error is reduced from approximately 10 degrees to less than 0.2 degrees with the use of this algorithm.
  • the fully-constrained attitude command feed-forward control algorithm first computes the time derivative of the quaternion attitude command. Quaternion math (i.e. products, conjugates, etc.) is used to convert the quaternion and quaternion derivative into a body-frame rotational rate command that is also sent directly into the inner-loop rotational rate control system. As was the case for the vector pointing feed-forward control algorithm, steady-state pointing error in the nadir pointing mode is virtually eliminated.
  • the control system inner loop may take the measured body rotational rates, the desired body rotational rates, and the measured reaction wheel speeds as inputs to produce desired reaction wheel rotational accelerations for the reaction wheels 115.
  • Pointing command generator 105 may be an example of, or include aspects of, the corresponding elements described with reference to FIG. 2 .
  • Navigation system 110 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 2 and 4 .
  • Reaction wheels 115 rotate in response to the output signal generated by the outer loop. Reaction wheels 115 also orient the satellite 100 based upon the rotation of the at least one reaction wheel 115. Reaction wheels 115 and the reaction wheel tachometer 120 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 2 and 4 .
  • Attitude control system 125 may be an example of, or include aspects of, the corresponding elements described with reference to FIG. 2 .
  • FIG. 2 shows an example of an attitude control system 220 with dynamic inversion in accordance with aspects of the present disclosure.
  • the example shown includes a pointing command generator 200, a navigation system 205, reaction wheels 210, a reaction wheel tachometer 215, and the attitude control system 220.
  • the pointing command generator 200 provides a desired attitude for input to the attitude control system 220.
  • the navigation system 205 provides an estimated attitude and may provide measured body rotational rates as inputs to the attitude control system 220.
  • the reaction wheel tachometer 215 may monitor the reaction wheels 210 and provide measured reaction wheel speeds for input to the attitude control system 220.
  • a control system outer loop 225 takes an estimated attitude and a desired attitude as inputs and produces desired body rotational rates.
  • a control system inner loop 230 takes the measured body rotational rates, the desired body rotational rates, and the measured reaction wheel speeds as input to produce desired reaction wheel rotational accelerations for the reaction wheels 210.
  • Pointing command generator 200 may be an example of, or include aspects of, the corresponding elements described with reference to FIG. 1 .
  • Navigation system 205 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 1 and 4 .
  • Reaction wheels 210 and reaction wheel tachometer 215 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 1 and 4 .
  • Attitude control system 220 may be an example of, or include aspects of, the corresponding elements described with reference to FIG. 1 .
  • Attitude control system 220 may include outer loop 225 and inner loop 230.
  • Outer loop 225 determines an eigen vector to rotate the satellite from one orientation to another.
  • the outer loop further comprises receiving a desired orientation of the satellite and an estimated orientation of the satellite as inputs.
  • the outer loop further comprises executing a satellite orientation error command based on the estimated orientation of the satellite and the desired orientation of the satellite.
  • the eigen vector command executes to decompose the satellite orientation error command into a scalar component and a vector component.
  • the outer loop further comprises a feed-forward control system.
  • the feed-forward control system controls for timing errors between rotation of the at least one reaction wheel and pointing at a target.
  • the feed-forward control system accounts for motion of tracking dishes, antennae, cameras, robotic arms, and solar arrays associated with the satellite.
  • the feed-forward control system further comprises receiving the desired orientation as an input.
  • Outer loop 225 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 3 and 4 .
  • Inner loop 230 executes a non-linear dynamic inversion algorithm to output a signal to at least one reaction wheel 210 of a satellite.
  • the inner loop comprises receiving the eigen vector determined from the outer loop as an input.
  • the inner loop further comprises receiving measured satellite rotation rates from a navigation system 205 associated with the satellite.
  • the inner loop further comprises receiving a measured reaction wheel speed of the at least one reaction wheel from a rotation wheel tachometer associated with the satellite.
  • the inner loop further comprises determining a desired rotational acceleration for the at least one reaction wheel. In some examples, the inner loop further comprises receiving a mass moment of inertia tensor of at least one reaction wheel, a rotation axis vector of at least one reaction wheel, and a mass moment of inertia tensor of the satellite.
  • Inner loop 230 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 3 and 4 .
  • FIG. 3 shows an example of an outer loop 300 of an attitude control system, such as attitude control system 220 of FIG. 2 , in accordance with aspects of the present disclosure.
  • the example shown includes outer loop 300 and inner loop 330.
  • the control system outer loop 300 takes an estimated attitude 305 and a desired attitude 310 as inputs. Each of these attitudes may be represented by a quaternion.
  • the estimated attitude 305 and the desired attitude 310 may be combined to form an orientation error command, which may be passed to a decomposition function 315.
  • the decomposition function 315 may generate a vector component (i.e., an eigen vector) and a scalar component.
  • the desired attitude 310 may be passed to a feed-forward control 320.
  • the scalar component may be passed to a limiter 325.
  • the output of the limiter 325 may be combined with the vector component, and this combined product may then be combined with the output of the feed forward control 320 to produce a desired body rotational rate.
  • the desired body rotational rate may be passed out of the control system outer loop 300 to a control system inner loop 330.
  • Outer loop 300 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 2 and 4 .
  • Outer loop 300 may include estimated attitude 305, desired attitude 310, decomposition function 315, feed forward control 320, and limiter 325.
  • Inner loop 330 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 2 and 4 .
  • FIG. 4 shows an example of an inner loop 420 of an attitude control system, such as attitude control system 220 of FIG. 2 , in accordance with aspects of the present disclosure.
  • the example shown includes an outer loop 400, a navigation system 405, a reaction wheel tachometer 410, reaction wheels 415, and the inner loop 420.
  • a desired body rotational rate may be received from the control system outer loop 400.
  • the navigation system 405 may provide one or more measured body rotational rates.
  • the body rotational rates may be combined with the desired body rotational rate to produce a desired body rotational acceleration.
  • the reaction wheel tachometer 410 may provide measured wheel speeds.
  • a satellite moment of inertia (MOI) tensor 425 may represent the satellite mass moment of inertia.
  • a wheel MOI matrix 430 may be derived from the reaction wheel mass moment of inertia and rotation axes.
  • a cross product function 435 may take as input the measured body rotational rates, and a combination of the measured body rotational rates with the satellite MOI tensor 425 and the measured wheel speeds. The cross product function 435 may output an Euler moment.
  • An inverse function 440 may generate an inverse of the wheel MOI matrix 430.
  • the satellite MOI tensor 425 may be combined with the desired body rotational acceleration to produce a moment command, which may then be combined with the Euler moment.
  • the result may be combined with the inverted matrix 430 and passed to a limiter 445.
  • the limiter 445 may then output desired reaction wheel rotational accelerations and pass them to the reaction wheels 415.
  • Outer loop 400 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 2 and 3 .
  • reaction wheel tachometer 410 and reaction wheels 415 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 1 and 2 .
  • Inner loop 420 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 2 and 3 .
  • FIG. 5 shows an example of a process for orienting a satellite in accordance with aspects of the present disclosure.
  • these operations may be performed by a processor executing a set of codes to control functional elements of a system. Additionally or alternatively, the processes may be performed using special-purpose hardware. Generally, these operations may be performed according to the methods and processes described in accordance with aspects of the present disclosure. For example, the operations may be composed of various substeps, or may be performed in conjunction with other operations described herein.
  • a system determines an eigen vector to rotate a satellite from one orientation to another.
  • the operations of this step refer to, or are performed by, an outer loop as described with reference to FIGs. 2-4 .
  • a system executes a non-linear dynamic inversion algorithm to output a signal to at least one reaction wheel of the satellite.
  • the operations of this step refer to, or are performed by, an inner loop as described with reference to FIGs. 2-4 .
  • a system rotates the at least one reaction wheel in response to the output signal.
  • the operations of this step refer to, or are performed by, a reaction wheels as described with reference to FIGs. 1 , 2 , and 4 .
  • a system orients the satellite based upon the rotation of the at least one reaction wheel.
  • the operations of this step refer to, or are performed by, a reaction wheels as described with reference to FIGs. 1 , 2 , and 4 .
  • FIG. 6 shows an example of an attitude control system response graph 600 in accordance with aspects of the present disclosure.
  • the attitude control system response graph 600 represents an example of the improvement in satellite orientation times for the system disclosed herein when compared to alternative systems that do not implement the disclosed features.
  • Control system response graph 600 may include vertical axis 605, horizontal axis 610, first sun pointing command 615, first earth surface pointing command 620, second sun pointing command 625, first improved response 630, first comparison response 635, second improved response 640, second comparison response 645, third improved response 650, and third comparison response 655.
  • the control system response graph 600 represents a simulation study of the attitude control system disclosed herein (the improved responses) to existing technology (the comparison responses). Other simulations (not shown) compared the attitude response of a satellite with the wheels set initially to zero speed and another case with the wheels set initially to 5,000 RPM. The attitude response with 5,000 RPM initial wheel speed for the system disclosed herein was identical to the response with zero initial wheel speed. The stored momentum in the wheels was completely accounted for and did not result in unwanted coupling, overshoot, oscillation, or instability.
  • the vertical axis 605 may represent a satellite pointing error, measured in degrees.
  • the horizontal axis 610 may represent time, in seconds.
  • the first improved response 630 shows a reduced response time for the first sun pointing command 615 in comparison to the first comparison response 635.
  • the second improved response 640 shows a reduced response time for the first earth surface pointing command 620 in comparison to the second comparison response 645.
  • the third improved response 650 shows a reduced response time for the second sun pointing command 625 in comparison to the third comparison response 655.
  • modules may be implemented as a hardware circuit comprising custom very large scale integration (VLSI) circuits or gate arrays, off-the-shelf semiconductors such as logic chips, transistors, or other discrete components.
  • VLSI very large scale integration
  • a module may also be implemented in programmable hardware devices such as field programmable gate arrays, programmable array logic, programmable logic devices or the like.
  • Modules may also be implemented in software for execution by various types of processors.
  • An identified module of executable code may, for instance, comprise one or more physical or logical blocks of computer instructions that may, for instance, be organized as an object, procedure, or function. Nevertheless, the executables of an identified module need not be physically located together, but may comprise disparate instructions stored in different locations which, when joined logically together, comprise the module and achieve the stated purpose for the module.
  • a module of executable code could be a single instruction, or many instructions, and may even be distributed over several different code segments, among different programs, and across several memory devices.
  • operational data may be identified and illustrated herein within modules, and may be embodied in any suitable form and organized within any suitable type of data structure. The operational data may be collected as a single data set, or may be distributed over different locations including over different storage devices, and may exist, at least partially, merely as electronic signals on a system or network.

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Automation & Control Theory (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Mobile Radio Communication Systems (AREA)

Description

    BACKGROUND OF THE INVENTION 1. Field of the Invention
  • The present invention relates generally to orienting a satellite, and more specifically to orienting a satellite using eigen vector rotation and non-linear dynamic inversion.
  • 2. Discussion of the Related Art
  • Various systems and processes are known in the art for orienting a satellite. Satellites may be controlled by a system known as an attitude determination and control system (ADACS). For example, an ADACS system may control the attitude, or orientation, of a satellite using reaction wheels that turn according to software algorithms executed on a computer within the satellite. In some cases, a satellite reaction wheel system may be characterized by dynamics that are highly non-linear.
  • Some control algorithms rely on linear control techniques such as a proportional-integral-derivative (PID) control to control the attitude of the satellite. Linear control systems applied to non-linear systems can exhibit oscillations, overshoots, and even instability. Additionally, some ADACS systems respond to pointing commands in a manner that resembles a yaw, pitch, and roll sequence. This method of maneuvering may be inefficient in terms of energy use and time.
  • WO 2017/159156 A1 describes a method and system for spacecraft orientation control using an inner-loop control determining first control inputs for momentum exchange devices to control an orientation of the spacecraft and an outer-loop control determining second control inputs for thrusters of the spacecraft to concurrently control a pose of the spacecraft and a momentum stored by the momentum exchange devices of the spacecraft. The outer-loop control determines the second control inputs using a model of dynamics of the spacecraft including dynamics of the inner-loop control, such that the outer-loop control accounts for effects of actuation of the momentum exchange devices according to the first control inputs determined by the inner-loop control. The thrusters and the momentum exchange devices are controlled according the first and second control inputs.
  • US 9 745 082 B2 describes an attitude control system and method for a satellite based on four single degree-of-freedom control moment gyroscopes with variable speed flywheels (or reaction wheels) in a pyramid configuration, combined with path and endpoint constraint time-optimal control.
  • EP 0 926 066 A1 describes an attitude control system and method of a spacecraft in a state control device of a moving body constituted by a navigation dynamics, an actuator for driving the navigation dynamics, first controlling means for controlling the actuator in PID control in response to a first output signal outputted from the navigation dynamics and adding means for outputting a control signal for controlling the actuator in feedforward control in response to outside noise by adding an estimated value of the outside noise to a control signal outputted from the first controlling means.
  • SUMMARY
  • A method of orienting a satellite using eigen vector rotation and non-linear dynamic inversion is described. The method of satellite orientation control according to the present invention includes applying with a processing device a satellite orientation control system, the satellite orientation control system comprising a double feedback loop system, wherein an outer loop receives a desired orientation of the satellite and an estimated orientation of the satellite as inputs, determines an eigen vector to rotate the satellite from one orientation to another based on the desired orientation and the estimated orientation and outputs the determined eigen vector or desired body rotational rates determined based on the determined eigen vector, and an inner loop receives the eigen vector or the desired body rotational rates from the outer loop as an input and executes a non-linear dynamic inversion algorithm based on the eigen vector or the desired body rotational rates to output an output signal to at least one reaction wheel of the satellite, rotating the at least one reaction wheel in response to the output signal, and orienting the satellite based upon the rotation of the at least one reaction wheel.
  • A system arranged for orienting a satellite using eigen vector rotation and non-linear dynamic inversion is described. The system according to the present invention comprises a satellite including a pointing command generator, a navigation system and at least one reaction wheel and includes a satellite orientation control system communicatively connected to the pointing command generator, the navigation system and the at least one reaction wheel, the satellite orientation control system comprising a double feedback loop system including an outer loop and an inner loop, wherein the outer loop is configured to receive a desired orientation of the satellite from the pointing command generator and an estimated orientation of the satellite from the navigation system as inputs, determine an eigen vector to rotate the satellite from one orientation to another based on the desired orientation and the estimated orientation and output the determined eigen vector or desired body rotational rates determined based on the determined eigen vector to the inner loop, and wherein the inner loop is configured to receive the eigen vector or the desired body rotational rates from the outer loop as an input and execute a non-linear dynamic inversion algorithm based on the eigen vector or the desired body rotational rates to determine an output signal for at least one reaction wheel of the satellite; a memory; a processor device in communication with the memory and configured to apply the satellite orientation control system to output the output signal to the at least one reaction wheel of the satellite to cause the at least one reaction wheel to rotate in response to the output signal and orient the satellite based upon the rotation of the at least one reaction wheel.
  • In some examples of the method and system described above, the outer loop further comprises executing a satellite orientation error command based on the estimated orientation of the satellite and the desired orientation of the satellite.
  • In some examples of the method and system described above, the eigen vector command executes to decompose the satellite orientation error command into a scalar component and a vector component.
  • In some examples of the method and system described above, the outer loop further comprises a feed-forward control system.
  • In some examples of the method and system described above, the feed-forward control system controls for timing errors between rotation of the at least one reaction wheel and pointing at a target.
  • In some examples of the method and system described above, the feed-forward control system accounts for motion of tracking dishes, antennae, cameras, robotic arms, and solar arrays associated with the satellite. In some examples of the method and system described above, the feed-forward control system further comprises receiving the desired orientation as an input.
  • In some examples of the method and system described above, the inner loop further comprises receiving measured satellite rotation rates from a navigation system associated with the satellite.
  • In some examples of the method and system described above, the inner loop further comprises receiving a measured reaction wheel speed of the at least one reaction wheel from a rotation wheel tachometer associated with the satellite. In some examples of the method and system described above, the inner loop further comprises determining a desired rotational acceleration for the at least one reaction wheel. In some examples of the method and system described above, the inner loop further comprises receiving a mass moment of inertia tensor of at least one reaction wheel, a rotation axis vector of at least one reaction wheel, and a mass moment of inertia tensor of the satellite.
  • BRIEF DESCRIPTION OF THE DRAWINGS
    • FIG. 1 shows an example of a satellite in accordance with aspects of the present disclosure.
    • FIG. 2 shows an example of an attitude control system with dynamic inversion in accordance with aspects of the present disclosure.
    • FIG. 3 shows an example of an outer loop of the attitude control system of FIG. 2 in accordance with aspects of the present disclosure.
    • FIG. 4 shows an example of an inner loop of the attitude control system of FIG. 2 in accordance with aspects of the present disclosure.
    • FIG. 5 shows an example of a process for orienting a satellite in accordance with aspects of the present disclosure.
    • FIG. 6 shows an example of an attitude control system response graph in accordance with aspects of the present disclosure.
    DETAILED DESCRIPTION
  • The following description is not to be taken in a limiting sense but is made merely for the purpose of describing the general principles of exemplary embodiments. The scope of the invention should be determined with reference to the claims.
  • Reference throughout this specification to "one embodiment," "an embodiment," or similar language means that a particular feature, structure, or characteristic described in connection with the embodiment is included in at least one embodiment of the present invention. Thus, appearances of the phrases "in one embodiment," "in an embodiment," and similar language throughout this specification may, but do not necessarily, all refer to the same embodiment.
  • Furthermore, the described features, structures, or characteristics of the invention may be combined in any suitable manner in one or more embodiments. In the following description, numerous specific details are provided, such as examples of programming, software modules, user selections, network transactions, database queries, database structures, hardware modules, hardware circuits, hardware chips, etc., to provide a thorough understanding of embodiments of the invention.
  • In other instances, well-known structures, materials, or operations are not shown or described in detail to avoid obscuring aspects of the invention.
  • The present disclosure provides a satellite control system that exhibits improved stability and increased efficiency by implementing a non-linear dynamic inversion inner-loop control algorithm coupled with an eigen vector outer-loop control algorithm. Thus, the attitude determination and control system (ADACS) system may operate using commands to rotate directly about an eigen vector (i.e., to go from one orientation directly to another). Additionally, the outer-loop control system is augmented with a feed-forward control element to enhance pointing accuracy when tracking moving targets.
  • FIG. 1 shows an example of a satellite 100 in accordance with aspects of the present disclosure. Satellite 100 includes a pointing command generator 105, navigation system 110, reaction wheels 115, a reaction wheel tachometer 120, and an attitude control system 125 (which may also be referred to as a satellite control system).
  • The attitude control system 125 may provide for stability and robustness over a wide range of operational modes and on-orbit conditions; smooth, quick, and accurate control system response to attitude commands; energy efficient operation to maximize battery life; ease of control system 125 setup; and ability to expand the control system's 125 use to large and complex satellites 100. The system 125 achieves these results by implementing an eigen vector outer-loop control algorithm, combined with a non-linear dynamic inversion inner-loop control algorithm. These results represent an improvement on existing satellite 100 control technology, but they are not a comprehensive list of features or advantages of the present system.
  • First, the outer loop of the attitude control system 125 utilizes an eigen vector outer-loop control algorithm. That is, between any two orientations of the satellite 100 attitude, there exists a rotation axis, called the eigen vector, which leads the satellite 100 directly from one orientation to the other. By rotating along this axis, the time and energy needed to complete the attitude maneuver may be minimized.
  • Secondly, the inner loop of the attitude control system 125 utilizes non-linear dynamic inversion control. The system accounts for stored momentum in the reaction wheels 115 and the satellite 100 which produces a smooth, uncoupled, linear response on a highly non-linear system. As an example, the reaction wheels 115 for the satellite 100 may have a top speed of 10,000 RPM. However, the satellite 100 may become uncontrollable once the reaction wheels 115 have 4,000 RPM of stored momentum due to the non-linear coupling of the satellite 100 control axes.
  • Also, the dynamic inversion control for the attitude control system 125 may be designed to work on satellites 100 with multiple reaction wheels 115 (e.g., more than 3) and with the reaction wheel spin axes orientated in any direction. This aspect of the control system 125 may also be used on advanced satellites 100 with redundant and skewed-axis reaction wheels 115. The system disclosed herein also takes into account cross coupling between the satellite 100 control axes and off-axis reaction wheels 115.
  • The dynamic inversion control for the attitude control system 125 may operate based on the following parameters: the number of reaction wheels 115, the mass moment of inertia of each reaction wheel, the mass moment of inertia tensor of the satellite 100 (including the reaction wheels 115), the coordinates of the unit vectors pointing along each reaction wheel's spin axis, the torque and speed limit of each reaction wheel, the satellite 100 attitude slew rate limit, and the closed-loop bandwidth of the rotational rate and attitude control loops.
  • A control system that only operates on proportional error feedback will result in the satellite 100 attitude lagging the desired orientation when the command is changing. Changing commands occur when the satellite 100 is attempting to track a moving target such as a point on the Earth's surface or to be aligned with the nadir orientation. The lagging response can result in significant steady-state pointing errors (e.g., approximately 10 degrees).
  • This steady-state error can be reduced or eliminated by adding an integral control algorithm. However, the attitude of the satellite 100 is two integrals away from the reaction wheel torque commands, and there is little natural damping of the system. These two aspects make integral control unappealing because any integral control will result in oscillations and overshoot of the attitude commands.
  • Feed-forward control can accomplish the goal of near zero steady-state error while not inducing unwanted oscillations and overshoot. The only trade-off is that these algorithms are somewhat more difficult to develop. Despite the additional developmental work, two different feed-forward control algorithms solve the lagging response problem. One feed-forward algorithm is designed for pointing vector control (e.g. earth position pointing), and a second feed-forward algorithm is designed for fully constrained attitude commands such as nadir pointing. It is desired to have a feed-forward control system that will perfectly track the attitude commands in the absence of modeling errors and disturbances.
  • The vector pointing feed-forward control algorithm tracks the rate of change of the pointing vector command by using cross products and dot product derivatives. The output of this algorithm is a body-frame rotational rate command that is sent directly into the inner-loop rotational rate control system. Simulation studies demonstrate that the Earth pointing error is reduced from approximately 10 degrees to less than 0.2 degrees with the use of this algorithm.
  • The fully-constrained attitude command feed-forward control algorithm first computes the time derivative of the quaternion attitude command. Quaternion math (i.e. products, conjugates, etc.) is used to convert the quaternion and quaternion derivative into a body-frame rotational rate command that is also sent directly into the inner-loop rotational rate control system. As was the case for the vector pointing feed-forward control algorithm, steady-state pointing error in the nadir pointing mode is virtually eliminated.
  • The control system inner loop may take the measured body rotational rates, the desired body rotational rates, and the measured reaction wheel speeds as inputs to produce desired reaction wheel rotational accelerations for the reaction wheels 115.
  • Pointing command generator 105 may be an example of, or include aspects of, the corresponding elements described with reference to FIG. 2. Navigation system 110 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 2 and 4.
  • Reaction wheels 115 rotate in response to the output signal generated by the outer loop. Reaction wheels 115 also orient the satellite 100 based upon the rotation of the at least one reaction wheel 115. Reaction wheels 115 and the reaction wheel tachometer 120 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 2 and 4.
  • Attitude control system 125 may be an example of, or include aspects of, the corresponding elements described with reference to FIG. 2.
  • FIG. 2 shows an example of an attitude control system 220 with dynamic inversion in accordance with aspects of the present disclosure. The example shown includes a pointing command generator 200, a navigation system 205, reaction wheels 210, a reaction wheel tachometer 215, and the attitude control system 220.
  • The pointing command generator 200 provides a desired attitude for input to the attitude control system 220. The navigation system 205 provides an estimated attitude and may provide measured body rotational rates as inputs to the attitude control system 220. The reaction wheel tachometer 215 may monitor the reaction wheels 210 and provide measured reaction wheel speeds for input to the attitude control system 220.
  • A control system outer loop 225 takes an estimated attitude and a desired attitude as inputs and produces desired body rotational rates. A control system inner loop 230 takes the measured body rotational rates, the desired body rotational rates, and the measured reaction wheel speeds as input to produce desired reaction wheel rotational accelerations for the reaction wheels 210.
  • Pointing command generator 200 may be an example of, or include aspects of, the corresponding elements described with reference to FIG. 1. Navigation system 205 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 1 and 4.
  • Reaction wheels 210 and reaction wheel tachometer 215 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 1 and 4.
  • Attitude control system 220 may be an example of, or include aspects of, the corresponding elements described with reference to FIG. 1. Attitude control system 220 may include outer loop 225 and inner loop 230.
  • Outer loop 225 determines an eigen vector to rotate the satellite from one orientation to another. The outer loop further comprises receiving a desired orientation of the satellite and an estimated orientation of the satellite as inputs. In some examples, the outer loop further comprises executing a satellite orientation error command based on the estimated orientation of the satellite and the desired orientation of the satellite. In some examples, the eigen vector command executes to decompose the satellite orientation error command into a scalar component and a vector component.
  • In some examples, the outer loop further comprises a feed-forward control system. In some examples, the feed-forward control system controls for timing errors between rotation of the at least one reaction wheel and pointing at a target. In some examples, the feed-forward control system accounts for motion of tracking dishes, antennae, cameras, robotic arms, and solar arrays associated with the satellite. In some examples, the feed-forward control system further comprises receiving the desired orientation as an input.
  • Outer loop 225 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 3 and 4.
  • Inner loop 230 executes a non-linear dynamic inversion algorithm to output a signal to at least one reaction wheel 210 of a satellite. The inner loop comprises receiving the eigen vector determined from the outer loop as an input. In some examples, the inner loop further comprises receiving measured satellite rotation rates from a navigation system 205 associated with the satellite. In some examples, the inner loop further comprises receiving a measured reaction wheel speed of the at least one reaction wheel from a rotation wheel tachometer associated with the satellite.
  • In some examples, the inner loop further comprises determining a desired rotational acceleration for the at least one reaction wheel. In some examples, the inner loop further comprises receiving a mass moment of inertia tensor of at least one reaction wheel, a rotation axis vector of at least one reaction wheel, and a mass moment of inertia tensor of the satellite.
  • Inner loop 230 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 3 and 4.
  • FIG. 3 shows an example of an outer loop 300 of an attitude control system, such as attitude control system 220 of FIG. 2, in accordance with aspects of the present disclosure. The example shown includes outer loop 300 and inner loop 330.
  • The control system outer loop 300 takes an estimated attitude 305 and a desired attitude 310 as inputs. Each of these attitudes may be represented by a quaternion. The estimated attitude 305 and the desired attitude 310 may be combined to form an orientation error command, which may be passed to a decomposition function 315. The decomposition function 315 may generate a vector component (i.e., an eigen vector) and a scalar component. The desired attitude 310 may be passed to a feed-forward control 320.
  • The scalar component may be passed to a limiter 325. The output of the limiter 325 may be combined with the vector component, and this combined product may then be combined with the output of the feed forward control 320 to produce a desired body rotational rate. The desired body rotational rate may be passed out of the control system outer loop 300 to a control system inner loop 330.
  • Outer loop 300 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 2 and 4. Outer loop 300 may include estimated attitude 305, desired attitude 310, decomposition function 315, feed forward control 320, and limiter 325.
  • Inner loop 330 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 2 and 4.
  • FIG. 4 shows an example of an inner loop 420 of an attitude control system, such as attitude control system 220 of FIG. 2, in accordance with aspects of the present disclosure. The example shown includes an outer loop 400, a navigation system 405, a reaction wheel tachometer 410, reaction wheels 415, and the inner loop 420.
  • A desired body rotational rate may be received from the control system outer loop 400. The navigation system 405 may provide one or more measured body rotational rates. The body rotational rates may be combined with the desired body rotational rate to produce a desired body rotational acceleration. The reaction wheel tachometer 410 may provide measured wheel speeds.
  • A satellite moment of inertia (MOI) tensor 425 may represent the satellite mass moment of inertia. A wheel MOI matrix 430 may be derived from the reaction wheel mass moment of inertia and rotation axes. A cross product function 435 may take as input the measured body rotational rates, and a combination of the measured body rotational rates with the satellite MOI tensor 425 and the measured wheel speeds. The cross product function 435 may output an Euler moment. An inverse function 440 may generate an inverse of the wheel MOI matrix 430.
  • The satellite MOI tensor 425 may be combined with the desired body rotational acceleration to produce a moment command, which may then be combined with the Euler moment. The result may be combined with the inverted matrix 430 and passed to a limiter 445. The limiter 445 may then output desired reaction wheel rotational accelerations and pass them to the reaction wheels 415.
  • In one example, the vector of reaction wheel angular acceleration may be given by: α w = I w 1 I b α b des w b x I b w b + I w w w
    Figure imgb0001
    where I w = I 1 u 1 x I 2 u 2 x I 3 u 3 x I n u nx I 1 u 1 y I 2 u 2 y I 3 u 3 y I n u ny I 1 u 1 x I 2 u 2 z I 3 u 3 z I n u nz ,
    Figure imgb0002
    w b = 0 w bz w by w bz 0 w bx w by w bx 0 ,
    Figure imgb0003
    α b des = K I w b des w b ,
    Figure imgb0004
    And where Ib is the body mass moment of inertia (including the reaction wheels 415), w b
    Figure imgb0005
    is the body frame angular rate vector with respect to inertial space, α b
    Figure imgb0006
    is the body frame acceleration vector with respect to inertial space, Iw is the wheel mass moment of inertia matrix, w w
    Figure imgb0007
    is the vector of reaction wheel angular rotation, and α w
    Figure imgb0008
    is the vector of reaction wheel 415 angular acceleration.
  • Outer loop 400 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 2 and 3.
  • Navigation system 405, reaction wheel tachometer 410 and reaction wheels 415 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 1 and 2.
  • Inner loop 420 may be an example of, or include aspects of, the corresponding elements described with reference to FIGs. 2 and 3.
  • FIG. 5 shows an example of a process for orienting a satellite in accordance with aspects of the present disclosure. In some examples, these operations may be performed by a processor executing a set of codes to control functional elements of a system. Additionally or alternatively, the processes may be performed using special-purpose hardware. Generally, these operations may be performed according to the methods and processes described in accordance with aspects of the present disclosure. For example, the operations may be composed of various substeps, or may be performed in conjunction with other operations described herein.
  • At step 500, a system determines an eigen vector to rotate a satellite from one orientation to another. The operations of this step refer to, or are performed by, an outer loop as described with reference to FIGs. 2-4.
  • At step 505, a system executes a non-linear dynamic inversion algorithm to output a signal to at least one reaction wheel of the satellite. The operations of this step refer to, or are performed by, an inner loop as described with reference to FIGs. 2-4.
  • At step 510, a system rotates the at least one reaction wheel in response to the output signal. The operations of this step refer to, or are performed by, a reaction wheels as described with reference to FIGs. 1, 2, and 4.
  • At step 515, a system orients the satellite based upon the rotation of the at least one reaction wheel. The operations of this step refer to, or are performed by, a reaction wheels as described with reference to FIGs. 1, 2, and 4.
  • FIG. 6 shows an example of an attitude control system response graph 600 in accordance with aspects of the present disclosure. The attitude control system response graph 600 represents an example of the improvement in satellite orientation times for the system disclosed herein when compared to alternative systems that do not implement the disclosed features.
  • Control system response graph 600 may include vertical axis 605, horizontal axis 610, first sun pointing command 615, first earth surface pointing command 620, second sun pointing command 625, first improved response 630, first comparison response 635, second improved response 640, second comparison response 645, third improved response 650, and third comparison response 655.
  • The control system response graph 600 represents a simulation study of the attitude control system disclosed herein (the improved responses) to existing technology (the comparison responses). Other simulations (not shown) compared the attitude response of a satellite with the wheels set initially to zero speed and another case with the wheels set initially to 5,000 RPM. The attitude response with 5,000 RPM initial wheel speed for the system disclosed herein was identical to the response with zero initial wheel speed. The stored momentum in the wheels was completely accounted for and did not result in unwanted coupling, overshoot, oscillation, or instability.
  • The vertical axis 605 may represent a satellite pointing error, measured in degrees. The horizontal axis 610 may represent time, in seconds. The first improved response 630 shows a reduced response time for the first sun pointing command 615 in comparison to the first comparison response 635. The second improved response 640 shows a reduced response time for the first earth surface pointing command 620 in comparison to the second comparison response 645. The third improved response 650 shows a reduced response time for the second sun pointing command 625 in comparison to the third comparison response 655.
  • Some of the functional units described in this specification have been labeled as modules, or components, to more particularly emphasize their implementation independence. For example, a module may be implemented as a hardware circuit comprising custom very large scale integration (VLSI) circuits or gate arrays, off-the-shelf semiconductors such as logic chips, transistors, or other discrete components. A module may also be implemented in programmable hardware devices such as field programmable gate arrays, programmable array logic, programmable logic devices or the like.
  • Modules may also be implemented in software for execution by various types of processors. An identified module of executable code may, for instance, comprise one or more physical or logical blocks of computer instructions that may, for instance, be organized as an object, procedure, or function. Nevertheless, the executables of an identified module need not be physically located together, but may comprise disparate instructions stored in different locations which, when joined logically together, comprise the module and achieve the stated purpose for the module.
  • Indeed, a module of executable code could be a single instruction, or many instructions, and may even be distributed over several different code segments, among different programs, and across several memory devices. Similarly, operational data may be identified and illustrated herein within modules, and may be embodied in any suitable form and organized within any suitable type of data structure. The operational data may be collected as a single data set, or may be distributed over different locations including over different storage devices, and may exist, at least partially, merely as electronic signals on a system or network.
  • While the invention herein disclosed has been described by means of specific embodiments, examples and applications thereof, numerous modifications and variations could be made thereto by those skilled in the art without departing from the scope of the invention set forth in the claims.

Claims (12)

  1. A system comprising:
    a satellite (100) including a pointing command generator (105, 200), a navigation system (110, 205) and at least one reaction wheel (115, 210, 415), and
    a satellite orientation control system (125, 220) communicatively connected to the pointing command generator (105, 200), the navigation system (110, 205) and the at least one reaction wheel (115, 210, 415), the satellite orientation control system (125, 220) comprising a double feedback loop system including an outer loop (225, 300, 400) and an inner loop (230, 330, 420),
    wherein the outer loop (225, 300, 400) is configured to receive a desired orientation (310) of the satellite (100) from the pointing command generator (105, 200) and an estimated orientation (305) of the satellite (100) from the navigation system (110, 205) as inputs, determine an eigen vector to rotate the satellite (100) from one orientation to another based on the desired orientation (310) and the estimated orientation (305) and output the determined eigen vector or desired body rotational rates determined based on the determined eigen vector to the inner loop (230, 330, 420), and
    wherein the inner loop (230, 330, 420) is configured to receive the eigen vector or the desired body rotational rates from the outer loop (225, 300, 400) as an input and execute a non-linear dynamic inversion algorithm based on the eigen vector or the desired body rotational rates to determine an output signal for the at least one reaction wheel (115, 210, 415) of the satellite (100) ;
    a memory; and
    a processor device in communication with the memory and configured to:
    apply the satellite orientation control system (125, 220) to output the output signal to the at least one reaction wheel (115, 210, 415) of the satellite (100) to cause the at least one reaction wheel (115, 210, 415) to rotate in response to the output signal and orient the satellite (100) based upon the rotation of the at least one reaction wheel (115, 210, 415).
  2. The system of claim 1, wherein:
    the outer loop (225, 300, 400) is further configured to execute a satellite orientation error command based on the estimated orientation (305) of the satellite (100) and the desired orientation (310) of the satellite (100).
  3. The system of claim 2, wherein:
    the eigen vector command, when executed, is configured to decompose the satellite orientation error command into a scalar component and a vector component.
  4. The system of claim 1, wherein:
    the outer loop (225, 300, 400) further comprises a feed-forward control system (320).
  5. The system of claim 4, wherein:
    the feed-forward control system (320) is configured to control for timing errors between rotation of the at least one reaction wheel (115, 210, 415) and pointing at a target.
  6. The system of claim 4, further comprising at least one of tracking dishes, antennae, cameras, robotic arms and solar arrays, wherein:
    the feed-forward control system (320) is configured to account for motion of tracking dishes, antennae, cameras, robotic arms, and solar arrays associated with the satellite.
  7. The system of claim 4, wherein:
    the feed-forward control system (320) is configured to receive the desired orientation (310) as an input.
  8. The system of claim 1, wherein:
    the inner loop (230, 330, 420) is further configured to receive measured satellite rotation rates from a navigation system (110, 205, 405) associated with the satellite (100).
  9. The system of claim 1, comprising a rotation wheel tachometer, wherein:
    the inner loop (230, 330, 420) is further configured to receive a measured reaction wheel speed of the at least one reaction wheel (115, 210) from the rotation wheel tachometer (120, 410) associated with the satellite (100).
  10. The system of claim 1, wherein:
    the inner loop (230, 330, 420) is further configured to determine a desired rotational acceleration for the at least one reaction wheel.
  11. The system of claim 1, wherein:
    the inner loop (230, 330, 420) is further configured to receive a mass moment of inertia tensor (430) of at least one reaction wheel (115, 210), a rotation axis vector of at least one reaction wheel (115, 210), and a mass moment of inertia tensor (425) of the satellite (100).
  12. A method of satellite orientation control, the method comprising:
    applying with a processing device a satellite orientation control system (125, 220), the satellite orientation control system (125, 220) comprising a double feedback loop system, wherein:
    an outer loop (225, 300, 400) receives a desired orientation (310) of the satellite (100) and an estimated orientation (305) of the satellite (100) as inputs, determines an eigen vector to rotate the satellite (100) from one orientation to another based on the desired orientation (310) and the estimated orientation (305) and outputs the determined eigen vector or desired body rotational rates determined based on the determined eigen vector, and
    an inner loop (230, 330, 420) receives the eigen vector or the desired body rotational rates from the outer loop (225, 300, 400) as an input and executes a non-linear dynamic inversion algorithm based on the eigen vector or the desired body rotational rates to output an output signal to at least one reaction wheel (115, 210, 415) of the satellite (100);
    rotating the at least one reaction wheel (115, 210, 415) in response to the output signal; and
    orienting the satellite (100) based upon the rotation of the at least one reaction wheel (115, 210, 415) .
EP19203064.1A 2018-10-25 2019-10-14 Satellite attitude control system using eigen vector, non-linear dynamic inversion, and feedforward control Active EP3643621B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
SI201930327T SI3643621T1 (en) 2018-10-25 2019-10-14 Satellite attitude control system using eigen vector, non-linear dynamic inversion, and feedforward control

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US16/170,157 US11279501B2 (en) 2018-10-25 2018-10-25 Satellite attitude control system using eigen vector, non-linear dynamic inversion, and feedforward control

Publications (2)

Publication Number Publication Date
EP3643621A1 EP3643621A1 (en) 2020-04-29
EP3643621B1 true EP3643621B1 (en) 2022-06-29

Family

ID=68342468

Family Applications (1)

Application Number Title Priority Date Filing Date
EP19203064.1A Active EP3643621B1 (en) 2018-10-25 2019-10-14 Satellite attitude control system using eigen vector, non-linear dynamic inversion, and feedforward control

Country Status (13)

Country Link
US (1) US11279501B2 (en)
EP (1) EP3643621B1 (en)
JP (2) JP7149247B2 (en)
CN (1) CN111103889B (en)
AU (1) AU2019240665B2 (en)
CA (1) CA3058338C (en)
DK (1) DK3643621T3 (en)
ES (1) ES2926461T3 (en)
HU (1) HUE059973T2 (en)
LT (1) LT3643621T (en)
PL (1) PL3643621T3 (en)
PT (1) PT3643621T (en)
SI (1) SI3643621T1 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102018118480A1 (en) * 2018-07-31 2020-02-06 Rockwell Collins Deutschland Gmbh Control system and method for regulating a swirl wheel device for stabilizing a spacecraft

Family Cites Families (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3460096A (en) * 1966-07-14 1969-08-05 Roger L Barron Self-organizing control system
DE3417661A1 (en) * 1983-05-13 1984-11-15 Mitsubishi Denki K.K., Tokio/Tokyo System for controlling the orientation of an artificial satellite
US4916622A (en) * 1988-06-16 1990-04-10 General Electric Company Attitude control system
FR2655167B1 (en) * 1989-11-29 1992-04-03 Aerospatiale METHOD OF CONTROLLING ATTITUDE IN ROLL AND LACET OF A SATELLITE.
US5100084A (en) * 1990-04-16 1992-03-31 Space Systems/Loral, Inc. Method and apparatus for inclined orbit attitude control for momentum bias spacecraft
US5611505A (en) * 1994-11-18 1997-03-18 Hughes Electronics Spacecraft energy storage, attitude steering and momentum management system
US5667171A (en) * 1995-04-28 1997-09-16 Hughes Aircraft Company Satellite spin axis stabilization using a single degree of freedom transverse momentum storage device
JP3185738B2 (en) 1997-12-25 2001-07-11 日本電気株式会社 Moving object state control apparatus and state control method
CA2338459A1 (en) 1998-07-23 2000-02-03 Douglas A. Staley System and method for spacecraft attitude control
US6463365B1 (en) * 2000-02-01 2002-10-08 Raytheon Company System and method for controlling the attitude of a space craft
JP2001260996A (en) 2000-03-23 2001-09-26 Toshiba Corp Attitude control device for spacecraft
US7630869B2 (en) * 2003-05-27 2009-12-08 University Of Washington Method for predicting vibrational characteristics of rotating structures
US7014150B2 (en) * 2004-07-30 2006-03-21 Honeywell International Inc. Method and system for optimizing torque in a CMG array
FR2918765B1 (en) * 2007-07-09 2009-10-02 Sagem Electronique Sa METHOD FOR DETERMINING A SERVO ERROR IN A CONTINUOUS LOOP OF A PSEUDO-RANDOM CODE.
CN100565405C (en) * 2008-09-12 2009-12-02 航天东方红卫星有限公司 A kind of spacecraft attitude control system of handling the unusual avoidance of rule
US8688296B2 (en) * 2008-11-17 2014-04-01 David A. Bailey Method for maximum data collection with a control moment gyroscope controlled satellite
CN101734379B (en) * 2009-12-22 2012-11-14 北京航空航天大学 FPGA-based highly-integrated high-precision control system for micro flywheel
US8918236B2 (en) * 2011-06-24 2014-12-23 Honeywell International Inc. Methods and systems for adjusting attitude using reaction wheels
CN102343985B (en) * 2011-07-08 2013-07-24 北京航空航天大学 Satellite time optimal posture maneuvering method with reaction flywheel
US9296474B1 (en) * 2012-08-06 2016-03-29 The United States of America as represented by the Administrator of the National Aeronautics & Space Administration (NASA) Control systems with normalized and covariance adaptation by optimal control modification
CN103674032B (en) * 2012-09-04 2016-02-24 西安电子科技大学 Merge the autonomous navigation of satellite system and method for pulsar radiation vector timing observation
US9745082B2 (en) 2015-06-02 2017-08-29 The Charles Stark Draper Laboratory, Inc. Rapid slew and settle systems for small satellites
US10005568B2 (en) * 2015-11-13 2018-06-26 The Boeing Company Energy efficient satellite maneuvering
US10180686B2 (en) * 2016-03-17 2019-01-15 Mitsubishi Electric Research Laboratories, Inc. Concurrent station keeping, attitude control, and momentum management of spacecraft
CN110235071B (en) 2016-11-10 2023-02-17 俄亥俄大学 Automatic car guiding and trajectory tracking
CN108327927B (en) * 2018-01-17 2020-11-06 浙江大学 Reaction wheel set self-adaptive moment distribution control method based on microsatellite
US10647449B2 (en) * 2018-05-30 2020-05-12 The Boeing Company Indirect self-imaging systems and methods

Also Published As

Publication number Publication date
NZ757832A (en) 2021-05-28
CA3058338A1 (en) 2020-04-25
CN111103889B (en) 2024-04-26
CA3058338C (en) 2022-02-15
AU2019240665A1 (en) 2020-05-14
AU2019240665B2 (en) 2021-06-03
US20200130869A1 (en) 2020-04-30
US11279501B2 (en) 2022-03-22
EP3643621A1 (en) 2020-04-29
JP2022121605A (en) 2022-08-19
SI3643621T1 (en) 2022-11-30
LT3643621T (en) 2022-10-10
JP7149247B2 (en) 2022-10-06
JP2020066427A (en) 2020-04-30
DK3643621T3 (en) 2022-09-05
ES2926461T3 (en) 2022-10-26
CN111103889A (en) 2020-05-05
PL3643621T3 (en) 2022-11-21
HUE059973T2 (en) 2023-01-28
PT3643621T (en) 2022-09-14

Similar Documents

Publication Publication Date Title
Wang et al. Coordinated stabilization of tumbling targets using tethered space manipulators
Bedrossian et al. Zero-propellant maneuver guidance
EP0453096B1 (en) Apparatus for inclined orbit attitude control of a satellite in an inclined orbit
CN107850900A (en) Quick rotation and installing system for grapefruit satellite
Nanos et al. On the dynamics and control of free-floating space manipulator systems in the presence of angular momentum
US20220097872A1 (en) Attitude control system and method
Matouk et al. Second-order sliding mode for Position and Attitude tracking control of Quadcopter UAV: Super-Twisting Algorithm
Sun et al. Adaptive control of space proximity missions with constrained relative states, faults and saturation
Saif Feedback linearization control of quadrotor with tiltable rotors under wind gusts
EP3643621B1 (en) Satellite attitude control system using eigen vector, non-linear dynamic inversion, and feedforward control
US6471161B1 (en) Satellite attitude control system
Aghili Time-optimal detumbling control of spacecraft
Kouhi et al. Robust control of a spin-stabilized spacecraft via a 1DoF gimbaled-thruster and two reaction wheels
Dizadji et al. 3-Axis attitude control of satellite using adaptive direct fuzzy controller
Yoon Spacecraft attitude and power control using variable speed control moment gyros
Sun et al. Tether attachment point stabilization of noncooperative debris captured by a tethered space system
Zhang The output torque estimation of MCMG for agile satellites
Jin Attitude control of underactuated and momentum-biased satellite using state-dependent Riccati equation method
Shen et al. Constrained attitude control of agile spacecraft using CMGs
NZ758060B2 (en) Transcatherer atrial sealing skirt, anchor, and tether
Karpenko et al. Attitude guidance: reducing the cost of rotating from ‘A’to ‘B’
Lee et al. Non-certainty-equivalent adaptive satellite attitude control using solar radiation pressure
Elkhayat et al. Backstepping control using single gimbal control moment gyro with parameter uncertainties
NZ758060A (en) Patient interface and method for making same
Cherd et al. Optimal control of a tetrahedral configured reaction wheels for Quaternion rotation based on rigid satellite model

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE APPLICATION HAS BEEN PUBLISHED

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

AX Request for extension of the european patent

Extension state: BA ME

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE

17P Request for examination filed

Effective date: 20200626

RBV Designated contracting states (corrected)

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: EXAMINATION IS IN PROGRESS

17Q First examination report despatched

Effective date: 20210430

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: GRANT OF PATENT IS INTENDED

RIC1 Information provided on ipc code assigned before grant

Ipc: B64G 1/28 20060101ALI20220103BHEP

Ipc: B64G 1/24 20060101AFI20220103BHEP

INTG Intention to grant announced

Effective date: 20220125

RIN1 Information on inventor provided before grant (corrected)

Inventor name: DERRICK II, JOHN BENTON

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE PATENT HAS BEEN GRANTED

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

REG Reference to a national code

Ref country code: AT

Ref legal event code: REF

Ref document number: 1501205

Country of ref document: AT

Kind code of ref document: T

Effective date: 20220715

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: DE

Ref legal event code: R096

Ref document number: 602019016379

Country of ref document: DE

REG Reference to a national code

Ref country code: DK

Ref legal event code: T3

Effective date: 20220830

REG Reference to a national code

Ref country code: RO

Ref legal event code: EPE

REG Reference to a national code

Ref country code: FI

Ref legal event code: FGE

REG Reference to a national code

Ref country code: SE

Ref legal event code: TRGR

REG Reference to a national code

Ref country code: PT

Ref legal event code: SC4A

Ref document number: 3643621

Country of ref document: PT

Date of ref document: 20220914

Kind code of ref document: T

Free format text: AVAILABILITY OF NATIONAL TRANSLATION

Effective date: 20220908

REG Reference to a national code

Ref country code: NL

Ref legal event code: FP

REG Reference to a national code

Ref country code: NO

Ref legal event code: T2

Effective date: 20220629

REG Reference to a national code

Ref country code: ES

Ref legal event code: FG2A

Ref document number: 2926461

Country of ref document: ES

Kind code of ref document: T3

Effective date: 20221026

Ref country code: GR

Ref legal event code: EP

Ref document number: 20220401907

Country of ref document: GR

Effective date: 20221010

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: HR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220629

REG Reference to a national code

Ref country code: EE

Ref legal event code: FG4A

Ref document number: E022682

Country of ref document: EE

Effective date: 20220913

REG Reference to a national code

Ref country code: SK

Ref legal event code: T3

Ref document number: E 40500

Country of ref document: SK

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: RS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220629

REG Reference to a national code

Ref country code: HU

Ref legal event code: AG4A

Ref document number: E059973

Country of ref document: HU

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SM

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220629

REG Reference to a national code

Ref country code: DE

Ref legal event code: R097

Ref document number: 602019016379

Country of ref document: DE

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: AL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220629

PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

26N No opposition filed

Effective date: 20230330

P01 Opt-out of the competence of the unified patent court (upc) registered

Effective date: 20230804

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: RO

Payment date: 20230925

Year of fee payment: 5

Ref country code: NL

Payment date: 20230915

Year of fee payment: 5

Ref country code: MC

Payment date: 20230928

Year of fee payment: 5

Ref country code: IT

Payment date: 20230913

Year of fee payment: 5

Ref country code: IE

Payment date: 20230912

Year of fee payment: 5

Ref country code: GB

Payment date: 20230831

Year of fee payment: 5

Ref country code: EE

Payment date: 20230905

Year of fee payment: 5

Ref country code: CZ

Payment date: 20230919

Year of fee payment: 5

Ref country code: BG

Payment date: 20230913

Year of fee payment: 5

Ref country code: LU

Payment date: 20230927

Year of fee payment: 5

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: SK

Payment date: 20230912

Year of fee payment: 5

Ref country code: SE

Payment date: 20230830

Year of fee payment: 5

Ref country code: PL

Payment date: 20230905

Year of fee payment: 5

Ref country code: GR

Payment date: 20230913

Year of fee payment: 5

Ref country code: FR

Payment date: 20230911

Year of fee payment: 5

Ref country code: BE

Payment date: 20230918

Year of fee payment: 5

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: LV

Payment date: 20230905

Year of fee payment: 5

Ref country code: LT

Payment date: 20230922

Year of fee payment: 5

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: ES

Payment date: 20231106

Year of fee payment: 5

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: IS

Payment date: 20231013

Year of fee payment: 5

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: TR

Payment date: 20231012

Year of fee payment: 5

Ref country code: SI

Payment date: 20230913

Year of fee payment: 5

Ref country code: PT

Payment date: 20231003

Year of fee payment: 5

Ref country code: NO

Payment date: 20231010

Year of fee payment: 5

Ref country code: MT

Payment date: 20231016

Year of fee payment: 5

Ref country code: HU

Payment date: 20230926

Year of fee payment: 5

Ref country code: FI

Payment date: 20231011

Year of fee payment: 5

Ref country code: DK

Payment date: 20231016

Year of fee payment: 5

Ref country code: DE

Payment date: 20230830

Year of fee payment: 5

Ref country code: CH

Payment date: 20231102

Year of fee payment: 5

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: CY

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220629

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220629

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: AL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20220629

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: AL

Payment date: 20231010

Year of fee payment: 5

PGRI Patent reinstated in contracting state [announced from national office to epo]

Ref country code: AL

Effective date: 20220830