Method for station keeping control of flexible spacecraft using onboard gain scheduling scheme

Claim 2

Apparatus for compensating for thrusting disturbances in spacecraft, said spacecraft containing a plurality of thrusters commanded by thruster selection and timing logic means which determines the appropriate thrusters to be fired during a given sample period in a given maneuver and computes the duration of firing for the thrusters, said plurality of thrusters being capable of producing torques about axes of the spacecraft, said apparatus comprising: rate sensor means for measuring the spacecraft's body rate about its axes, estimator means, responsive to signal inputs from the rate sensor means for estimating angular acceleration, body rate, and body angle of said spacecraft, said estimator means utilizing time-varying gains which decline exponentially and are determined in an iterative fashion in computing said estimated angular acceleration and body rate; and controller means for generating acceleration control commands from angular acceleration, body rate, and body angle estimates provided by said estimator means, said controller means utilizing time-varying gains which decline exponentially and are determined in an iterative fashion in computing said acceleration control commands.

Claim 3

The apparatus of claim 1 further comprising means for generating feedforward signals indicating an actual rate change commanded during a previous sample period by thruster selection and timing logic about the spacecraft's axes for use by said estimator means for estimating angular acceleration, body rate and body angle of said spacecraft.

Claim 4

The apparatus of claim 1 wherein the gains used by said estimator means are chosen for a predetermined initial value and decline to a predetermined final, lower steady-state value.

Claim 5

The apparatus of claim 1 wherein the gains used by said controller means are chosen for a predetermined initial value and decline to a final, lower steady-state value.

Claim 6

The apparatus of claim 1 wherein said controller means further comprises filter means for compensating said acceleration control signal in order to achieve stability of structural modes of said spacecraft.

Claim 7

The apparatus of claim 5, wherein said filter means comprises a notch filter.

Claim 8

The apparatus of claim 1 wherein said time-varying estimator means gains and said time-varying controller means gains are determined in an iterative fashion by use of a difference equation.

Claim 9

The apparatus of claim 7 wherein the gains used by said estimator means are calculated in an iterative fashion by use of a homogeneous first-order difference equation.

Claim 10

The apparatus of claim 7 wherein the gains used by said controller means are calculated in an iterative fashion by use of a homogeneous first-order difference equation.

Claim 11

The apparatus of claim 7 wherein the gains used by said estimator means are calculated in an iterative fashion according to the relations K.sub.d (n)=K.sub.d (n-1)+Ω.sup.e T(K.sub.d.sup.f -K.sub.d (n-1)) and K.sub.r (n)=K.sub.r (n-1)+Ω.sup.e T(K.sub.r.sup.f -K.sub.r (n-1)) where K d (n) represents the gain used in computing the angular acceleration estimate at the nth iteration K r (n) represents the gain used in computing the body rate estimate at the nth iteration Ω e represents the inverse of a predetermined system time constant K d f represents a predetermined final value for K d K r f represents a predetermined final value for K r ; and T represents the chosen length of the sample period, and the resulting gains as a function of time satisfy K.sub.d (t)=K.sub.d.sup.f +(K.sub.d (0)-K.sub.d.sup.f)e.sup.-t/τ and K.sub.r (t)=K.sub.r.sup.f +(K.sub.r (0)-K.sub.r.sup.f)e.sup.-t/τ.

Claim 12

The apparatus of claim 7 wherein the gains used by said controller means are calculated in an iterative fashion according to the relations C.sub.d (n)=C.sub.d (n-1)+Ω.sup.c T(C.sub.d.sup.f -C.sub.d (n-1)), C.sub.r (n)=C.sub.r (n-1)+Ω.sup.c T(C.sub.r.sup.f -C.sub.r (n-1)), C.sub.p (n)=C.sub.p (n-1)+Ω.sup.c T(C.sub.p.sup.f -C.sub.p (n-1)), where C d (n) represents the gain applied to the angular acceleration estimate at the nth iteration C r (n) represents the gain applied to the body rate estimate at the nth iteration C p (n) represents the gain applied to the body angle estimate at the nth iteration Ω c represents the inverse of the system time constant C d f represents a predetermined final value for C d C r f represents a predetermined final value for C r C p f represents a predetermined final value for C p ; and T represents the chosen length of the sample period, and the resulting gains as a function of time satisfy C.sub.d =C.sub.d.sup.f +(C.sub.d (0)-C.sub.d.sup.f)e.sup.-t/τ, C.sub.p =C.sub.p.sup.f +(C.sub.p (0)-C.sub.p.sup.f)e.sup.-t/τ, and C.sub.r =C.sub.r.sup.f +(C.sub.r (0)-C.sub.r.sup.f)e.sup.-t/τ.

Claim 13

Apparatus for compensating for thrusting disturbances in spacecraft during stationkeeping operations, said spacecraft containing a plurality of thrusters commanded by thruster selection and timing logic means which determine the appropriate thrusters to be fired during a given sample period in a given maneuver and computes the duration of firing for the thrusters, said plurality of thrusters being capable of producing torques about each axis of the spacecraft, said apparatus comprising: rate sensor means for measuring the spacecraft's body rate about each of its three axes, means for generating feedforward signals indicating the actual rate change commanded during the previous sample period by the thruster selection and timing logic about each of the spacecraft's axes, estimator means, responsive to signal inputs from the spacecraft rate sensors and feedforward signals from the thruster selection and timing logic, for estimating disturbance torque, body rate, and body angle of said spacecraft, said estimator means using calculated time-varying gains calculated iteratively according to the relations K.sub.d (n)=K.sub.d (n-1)+Ω.sup.e T(K.sub.d.sup.f -K.sub.d (n-1)) and K.sub.r (n)=K.sub.r (n-1)+Ω.sup.e T(K.sub.r.sup.f -K.sub.r (n-1)) where K d (n) represents the gain used in computing the angular acceleration estimate at the nth iteration K r (n) represents the gain used in computing the body rate estimate at the nth iteration Ω e represents the inverse of a predetermined system time constant K d f represents a predetermined final value for K d K r f represents a predetermined final value for K r ; and T represents the chosen length of the sample period, and the resulting gains as a function of time satisfy K.sub.d (t)=K.sub.d.sup.f +(K.sub.d (0)-k.sub.d.sup.f)e.sup.-t/τ and K.sub.r (t)=K.sub.r.sup.f +(K.sub.r (0)-K.sub.r.sup.f)e.sup.-t/τ, where τ and Ω e are related by Ω e= 1/τ, resulting in rapid convergence of said rate and disturbance estimates and obviating the need for spacecraft prebias maneuvering while maintaining stability of the structural modes of said spacecraft; and controller means for generating acceleration control commands from disturbance torque, body rate, and body angle values provided by said estimator means, said controller means using calculated time-varying gains calculated in an iterative fashion according to the relations C.sub.d (n)=C.sub.d (n-1)+Ω.sup.c T(C.sub.d.sup.f -C.sub.d (n-1)), C.sub.r (n)=C.sub.r (n-1)+Ω.sup.c T(C.sub.r.sup.f -C.sub.r (n-1)), C.sub.p (n)=C.sub.p (n-1)+Ω.sup.c T(C.sub.p.sup.f -C.sub.p (n-1)), where C d (n) represents the gain applied to the angular acceleration estimate at the nth iteration C r (n) represents the gain applied to the body rate estimate at the nth iteration C p (n) represents the gain applied to the body angle estimate at the nth iteration Ω c represents the inverse of the system time constant C d f represents a predetermined final value for C d C r f represents a predetermined final value for C r C p f represents a predetermined final value for C p ; and T represents the chosen length of the sample period, and the resulting gains as a function of time satisfy C.sub.d =C.sub.d.sup.f +(C.sub.d (0)-C.sub.d.sup.f)e.sup.-t/τ, C.sub.p =C.sub.p.sup.f +(C.sub.p (0)-C.sub.p.sup.f)e.sup.-t/τ, and C.sub.r =C.sub.r.sup.f +(C.sub.r (0)-C.sub.r.sup.f)e.sup.-t/τ, where τ and Ω c are related by Ω c =1/τ, resulting in rapid convergence of said acceleration control estimate and obviating the need for prebias maneuvering, said controller means further comprising notch filter means for compensating said acceleration control signal in order to achieve stability of the structural modes of said spacecraft.

Claim 14

Apparatus for compensating for thrusting disturbances in spacecraft during stationkeeping operations, said spacecraft containing a plurality of thrusters commanded by thruster selection and timing logic means which determine the appropriate thrusters to be fired during a given sample period in a given maneuver and computes the duration of firing for the thrusters, said plurality of thrusters being capable of producing torques about each axis of the spacecraft, said apparatus comprising: rate sensor means for measuring the spacecraft's body rate about each of its three axes, means for generating feedforward signals indicating the actual rate change commanded during the previous sample period by the thruster selection and timing logic about each of the spacecraft's axes, estimator means, responsive to signal inputs from the spacecraft rate sensors and feedforward signals from the thruster selection and timing logic, for estimating disturbance torque, body rate, and body angle of said spacecraft, said estimator means using calculated time-varying gains calculated iteratively according to the relations K.sub.d (n)=K.sub.d (n-1)+Ω.sup.e T(K.sub.d.sup.f -K.sub.d (n-1)) and K.sub.r (n)=K.sub.r (n-1)+Ω.sup.e T(K.sub.r.sup.f -K.sub.r (n-1)) where K d (n) represents the gain used in computing the angular acceleration estimate at the nth iteration K r (n) represents the gain used in computing the body rate estimate at the nth iteration Ω e represents the inverse of a predetermined system time constant K d f represents a predetermined final value for K d K r f represents a predetermined final value for K r ; and T represents the chosen length of the sample period, and the resulting gains as a function of time satisfy K.sub.d (t)=K.sub.d.sup.f +(K.sub.d (0)-k.sub.d.sup.f)e.sup.-t/τ and K.sub.r (t)=K.sub.r.sup.f +(K.sub.r (0)-K.sub.r.sup.f)e.sup.-t/τ, where τ and Ω e are related by Ω e= 1/τ, resulting in rapid convergence of said rate and disturbance estimates and obviating the need for spacecraft prebias maneuvering while maintaining stability of the structural modes of said spacecraft; and controller means for generating acceleration control commands from disturbance torque, body rate, and body angle values provided by said estimator means, said controller means using calculated time-varying gains calculated in an iterative fashion according to the relations C.sub.d (n)=C.sub.d (n-1)+Ω.sup.c T(C.sub.d.sup.f -C.sub.d (n-1)), C.sub.r (n)=C.sub.r (n-1)+Ω.sup.c T(C.sub.r.sup.f -C.sub.r (n-1)), C.sub.p (n)=C.sub.p (n-1)+Ω.sup.c T(C.sub.p.sup.f -C.sub.p (n-1)), where C d (n) represents the gain applied to the angular acceleration estimate at the nth iteration C r (n) represents the gain applied to the body rate estimate at the nth iteration C p (n) represents the gain applied to the body angle estimate at the nth iteration Ω c represents the inverse of the system time constant C d f represents a predetermined final value for C d C r f represents a predetermined final value for C r C p f represents a predetermined final value for C p ; and T represents the chosen length of the sample period, and the resulting gains as a function of time satisfy C.sub.d =C.sub.d.sup.f +(C.sub.d (0)-C.sub.d.sup.f)e.sup.-t/τ, C.sub.p =C.sub.p.sup.f +(C.sub.p (0)-C.sub.p.sup.f)e.sup.-t/τ, and C.sub.r =C.sub.r.sup.f +(C.sub.r (0)-C.sub.r.sup.f)e.sup.-t/τ, where τ and Ω c are related by Ω c =1/τ, resulting in rapid convergence of said acceleration control estimate and obviating the need for prebias maneuvering, said controller means further comprising notch filter means for compensating said acceleration control signal in order to achieve stability of the structural modes of said spacecraft.

Claim 17

The method of claim 13 further comprising the step of generating feedforward signals indicating an actual body rate change commanded during a previous sample period about the spacecraft's axes for use in the step of estimating the angular acceleration, the body rate and the body angle of the spacecraft.

Claim 18

The method of claim 13 wherein the step of estimating comprises calculating the gains in an iterative fashion by use of a homogenous first-order difference equation.

Claim 19

The method of claim 13 wherein the step of generating comprises calculating the gains in an iterative fashion by use of a homogenous first-order difference equation.

Claim 20

The method of any one of claim 13 wherein the step of estimating comprises calculating the gains in an iterative fashion according to the relations K.sub.d (n)=K.sub.d (n-1)+Ω.sup.e T(K.sub.d.sup.f -K.sub.d (n-1)) and K.sub.r (n)=K.sub.r (n-1)+Ω.sup.e T(K.sub.r.sup.f -K.sub.r (n-1)) where K d (n) represents the gain used in computing the angular acceleration estimate at the nth iteration, K r (n) represents the gain used in computing the body rate estimate at the nth iteration Ω e represents the inverse of a predetermined system time constant K d f represents a predetermined final value for K d K r f represents a predetermined final value for K r ; and T represents the chosen length of the sample period, and the resulting gains as a function of time which satisfy K.sub.d (t)=K.sub.d.sup.f +(K.sub.d (0)-K.sub.d.sup.f)e.sup.-t/τ and K.sub.r (t)=K.sub.r f+(K.sub.r (0)-K.sub.r.sup.f)e.sup.-t/τ.

Claim 21

The method of claim 13 wherein the step of generating comprises calculating the gains in an iterative fashion according to the relations C.sub.d (n)=C.sub.d (n-1)+Ω.sup.c T(C.sub.d.sup.f -C.sub.d (n-1)), C.sub.r (n)=C.sub.r (n-1)+Ω.sup.c T(C.sub.r.sup.f -C.sub.r (n-1)), C.sub.p (n)=C.sub.p (n-1)+Ω.sup.c T(C.sub.p.sup.f -C.sub.p (n-1)), where C d (n) represents the gain applied to the angular acceleration estimate at the nth iteration C r (n) represents the gain applied to the body rate estimate at the nth iteration C p (n) represents the gain applied to the body angle estimate at the nth iteration Ω c represents the inverse of the system time constant C d f represents a predetermined final value for C d C r f represents a predetermined final value for C r C p f represents a predetermined final value for C p ; and T represents the chosen length of the sample period, and the resulting gains as a function of time satisfy C.sub.d =C.sub.d.sup.f +(C.sub.d (0)-C.sub.d.sup.f)e.sup.-t/τ, C.sub.p =C.sub.p.sup.f +(C.sub.p (0)-C.sub.p.sup.f)e.sup.-t/τ, and C.sub.r =C.sub.r.sup.f +(C.sub.r (0)-C.sub.r.sup.f)e.sup.-t/τ.

Claim 22

The method of claim 13 further comprising the step of compensating the acceleration control signal in order to achieve stability of structural modes of the spacecraft.

Claim 23

The method of claim 19, wherein the step of compensating comprises using a notch filter.