AUTONOMOUS SPIN AXIS ATTITUDE CONTROLLER FOR A SPINNING SPACECRAFT

Abstract: A system located onboard a geosynchronous satellite (10) for detecting spin axis precession and developing thruster control signals for restoring the spin axis (12) to the desired attitude. Fuel economy is enhanced by timing the attitude trim maneuvers to coincide with certain station-keeping maneuvers.

Description: AUTONOMOUS SPIN AXIS ATTITUDE CONTROLLER FOR A SPINNING SPACECRAFT Background of the Invention 1. Field of the Invention. This invention relates to attitude control systems for spacecraft and , more particularly , to such control systems which can he self-contained on a satellite in geosynchronous or other type of orbit . 2. Description of the Prior Art. In systems employing geostationary spinning satellites, as perhaps with certain other earth satellite systems , it becomes neces sary f rom time to time to _ determine the spin axis attitude and to execute attitude tr im maneuvers .'* These maneuvers are needed in order to compensate for the spin axis precession which is caused by solar torque. In predecessor systems , such a tt i tude t r im maneuvers have been commanded f rom ground contro l stations. In such systems, data which was pertinent to the determination of satellite spin axis attitude was collected by onboard sensor s and telemetered to the controlling ground station. There the data was processed to determine the appropriate attitude trim corrections . These were then transmitted to the satellite as command signals to perform the necessary attitude trim maneuvers . This performance of attitude trim correction in real time under ground control requires an operator in the control loop, is subject to anomalies from operator error, transmitted signal interference and the like, and imposes an operational burden on ground control systems. It would be desirable, therefore, to have the required attitude trim maneuvers determined and controlled by systems which are entirely self-contained on the satellite. Summary of the Invention in brief, arrangements in accordance with the present invention comprise spinning earth sensors and an onboard state estimator for attitude determination., plus thruster firing control logic to automatically perform attitude correction. In accordance with a particular feature of the present invention, special control logic is utilized to combine attitude control functions for the optimization of fuel usage. For its various trajectory control functions, the satellite is provided with both axial and radial jets, or thrusters. The axial thrusters, normally two in number, for redundancy, are directed parallel to the spin axis of the satellite and on opposite sides thereof. The radial thrusters, also normally two in number for redundancy, are mounted along a radius of the satellite. On a spinning satellite, the necessary thrusting maneuvers can be performed with only two jets, one axial jet which is displaced from the spin axis and one radial jet. For north-south station keeping, the axial jet is operated continuously. For east-west station keeping, the radial jet is pulsed. For precession maneuvers and nutation control, the axial jet is operated in a pulsed mode. If more than one axial or radial jet is provided, the pulses are directed to a selected one of the jets. The orientation of the satellite at all points in its orbit should be such that its spin axis is normal to its orbital plane which intersects the center of the earth. A pair of earth sensors are mounted within the satellite to view the earth periodically during satellite rotation through a window in the side of the satellite. The diametral plane of the spinning satellite should coincide with the satellite's orbital plane. However, as the satellite spin axis precesses so that it is no longer normal to the orbital plane and the diametral plane of the satellite deviates from the orbital plane, the view of the earth detected by one of the sensors departs from coincidence with that detected by the other sensor. For example, as the diametral plane of the satellite dips below the center of the earth, the chord of the earth which is detected by the north earth sensor as its field of view sweeps the earth will become larger while the chord of the earth sensed by the south earth sensor as its field of view sweeps the earth becomes smaller. On the opposite side of the earth, 180° around the satellite orbit, the opposite situation will obtain: the chord sensed by the south earth sensor will be larger than the chord sensed by the north earth sensor. The difference in sensed chord length from a reference value (normally 13.4° for a sensor cant angle of 5.5° from the spin plane) is used as a measure of satellite roll angle. This difference value is taken from a selected one of the two earth sensors, either north or south depending on which is free of sun interference. As an alternative, the difference in chord length between the two sensor measurements is a measure of the earth aspect angle and may be used as a measurement of satellite roll angle. The three-state estimator onboard the satellite estimates the roll, yaw, and earth sensor bias using earth horizon sensor measurements. The earth sensor provides the satellite roll angle measurements while the earth orbit kinematic coupling gives the satellite yaw angle information. Following a determination of the spin-axis attitude and the attitude trim maneuvers which are required, appropriate signals are generated to fire the thrusters at selected points and intervals to achieve the desired attitude corrections. In accordance with an 5 aspect of the invention, the timing of thruster firing is coordinated with orbital position for more efficient overall fuel usage. If, however, it is determined that the precession errors exceed a preset threshold, the thrusters are pulsed to reduce those errors to the 0 threshold range. Brief Description of the Drawing A better understanding of the present invention may be had from a consideration of the following detailed description, taken in conjunction with the accompanying 5 drawing in which: Fig. 1 is a schematic representation of a spinning satellite which may be controlled by the present invention; Fig. 2 is an illustrative- view indicating the 0 frames of reference of the satellite relative to the earth and to the sun; Fig. 3 is a schematic representation of a satellite illustrating the operation of the earth sensors; Fig.. . 4 is a sketch of earth sensor output 25. waveforms; Fig. 5 is a functional block diagram of the spin axis controller (SAC) of the present invention; Fig. 6 is a simplified schematic block diagram representing the spin axis controller; and 30 Fig. 7 is a plot illustrating spin axis correction of a satellite controlled in accordance with the present invention. Description of the Preferred Embodiment The present invention has been developed for the 35 Intelsat VI satellite and the preferred embodiment will be described in the context of that spacecraft. However, this is not to be taken as limiting the scope of the invention, which may have application to other spacecraft operating under .other conditions and orbits . Referring to Fig. 1, a satellite 10 is represented schematically, spinning about a longitudinal spin axis 12. The satellite 10 is indicated as equipped with an axial jet or thruster 14 , displaced from the spin axis 12 , and a radial jet or thruster 16. The axial jet 14 may be one of a pair of thrusters displaced on opposite sides of the spin axis 12. For spin axis attitude correction, only one thruster of an axial pair is pulsed at a time. Fig . 2 shows the two coordinate frames used in the analysis . The orbit reference f rame (Xg Yg z θ ) i- 3 oriented as follows : The orbit Xg-axis is in the the direction of flight; the orbit Yg-axis is normal to the orbital plane; and the orbit Zg-axis is pointed toward the earth along the radius vector f rom the satellite to the center of the earth. The sun fixed reference frame (X s Y s Z s ) is defined as follows: the Y s -axis coincides with the orbit Yg-axis ; the projection of the sun vector on the orbit plane is at the Z s -axis , and the X s -axis completes a r ight-handed orthogonal triad . Fig . 2 also g ives the coordinate transformation between these two coordinate systems where roll ( φ ) and yaw ( ψ) angles are defined as the rotation angles about the Xg-axis and Z g-axis , respectively; U is the rotation angle about the X s -axis , and V is about the Z s -axis . The two coordinate systems coincide at local midnight. The orbit coordinate frame rotates about the sun frame , Y s -axis , with orbit rate of ωg. For Intelsat VI the nutation damping time constant is much shor ter than the geostationary orbit period ; hence , it can be justif ied that the nutational dynamics can be neglected in the spacecraft dynamical model for SAC analysis . Fig . 3 is a diagram illustrating the way in which the orientation of the spin axis of the satellite 10 is determined from a scan of the earth as the satellite rotates. The satellite is provided with a pair of earth sensors (see Fig. 6) which are arranged to view the earth through a window, 18 n or 18 s , each time the windows face the earth. The spin plane of the satellite 10, indicated by the broken line 20, should coincide with the orbital plane which passes through the center of the earth 22. When properly oriented, the spin axis 12 will coincide with the orbit normal arrow 24. The north earth sensor field of view through the window 18 n extends by an angle <% above the line 20. Similarly, the field of view of the south earth sensor through window 18 s extends by an ' angle OQ below the broken line 20. If the spin plane deviates from an intersection with the center of the earth 22, one of the sensors will view the earth for a longer interval than if no deviation exists and the other sensor will view the earth for a shorter interval. This is because the field of view of one sensor is closer to. the equator than that of the other sensor, which is then closer to one of the poles. In the drawing of Fig. 3, in which the spin plane 20 is shown shifted in a southerly direction from the center of the earth 22, the north earth sensor scans a longer earth chord and the south earth sensor scans a shorter earth chord. Fig. 4 shows the variation of the sensed chord length for different attitude conditions ψg which corresponds to zero shift of the spin plane should equal a reference value Z re f which is stored in the processor memory. The measured deviation of sensed chord length from the reference value is a measure of the earth aspect angle which is used in determining the extent to which the spin plane, or the spin axis, deviates from the correct attitude. In synchronous orbit the external disturbances on the spacecraft are dominated by the solar radiation pressure. The solar torque can be predicted and expressed in terms of a Fourier series. With the secular term in the solar torque Fourier expression and the spacecraft angular momentum, the spacecraft precession rate due to solar torque can be modeled as ♦ - _ o The linearized kinematic rotation equation of the spacecraft in roll ( φ) and yaw ( ψ) in synchronous orbit are ♦ " » 0 * * A o sin "oc * w l A * « -u o φ + Ao cos m o _c * w, _. where W]_ and 2 are white Gaussian noises representing primarily the neglected higher order terms in the solar torque Fourier expression. The earth sensor bias variation has been modeled as a random walk with input noise sufficient to cause its 1<? value to grow at a certain rate. Although the mean square value of a random walk process grows linearly with time and only seems to be suitable for short term operation, it can be justified because the control loop time constant (in the order of 1 hour) is much shorter than the characteristic time of the random walk process. For the analysis and computer simulation to be followed, the noise inputs wj_, W2 can be modeled as σ i * N(0,1) ϊ 7 where T s ■ sampling period which represents the computer simulation time interval T 0 - random walk characteristic time σi * one sigma or expected value of the noise after time T 0 N(0,1) ■ Gaussian random noise generator with zero mean and σ« 1 δ The spinning earth horizon sensor measures the earth chord length (Z) which can be related to the spacecraft roll angle, φ , by the following approximation (Fig. 3): _ — φ where γ = the earth angular radius, as ~ sensor cant angle, Φ= the spacecraft roll angle measured from the orbit normal, and Z is the sensor output. The roll attitude measurement Z n can be obtained at iteration period n from the following equation: Z n = K(Z ' a a v vg_ -Z rβf ) ♦ Z b - αcos(2γ n ) - øsin(2γ n ) + ω z where 2 avg = average chord length Z re f = referenced chord length, a fixed number Z D = co mandable bias, for transient free sensor switching α, 3 - commandable parameters to compensate for sensor cyclic bias γ- jj = orbital angle at iteration period N ω z = sensor noise can be modeled as Gaussian white noise with zero mean and with standard deviation of σ z The sensor noise will be averaged out after M measurements taken in accordance with u m -i * σ * N(0, 1) There is no direct yaw angle measurement available with the spinning earth sensors; the yaw information can be obtained through the orbit kinematic coupling as described below. The spacecraft control system and the particular portion thereof relating to the spin axis attitude controller of the present invention are represented in the functional and schematic block diagrams of Figs. 5 and 6. In the functional block diagram of Fig. 5, a simple model of the spacecraft dynamics, represented within the block 30, is shown where the forcing function of the system is the angular drift rate about the sunline (V). The block 30 is shown containing the horizon sensor 32 which is coupled to provide an input to an earth chord processor 34, connected to be enabled by a counter 36. The counter 36 also controls a state estimate stage 38 which, with the earth chord processor 34, provides an output to a measurement error determination stage 40. The result is a correction of the state estimates, developed in correction stage 42 and applied to a coordinate transformation stage 44. The output of the stage 44 is applied to a maneuver control logic stage 46, controlled by a sidereal clock counter 48 and providing an output to a thruster command processor 50. The processor 50 is followed by an inverse transformation stage 52 and an estimate correction stage 54 which provides an output to a thruster driver 56. In the schematic block diagram of Fig. 6, the system is represented as comprising north and south earth sensors 70, 71 having separate field-of-view sensing elements 72, 74 and a sun sensor 76 having a sensing element 78. The earth sensor 70 provides an output to a first earth chord processor 82 comprising the earth sensor processing electronics (ESPE-A). A second earth chord processor 84 containing the earth sensor processing electronics (ESPE-B) is connected to the south earth sensor 71 for alternative earth chord sensing. There is also a sun sensor 76 with element 78 to provide a pulse for controlling the spin phase angle for thruster firing. The outputs of the blocks 82, 84 are applied to a microprocessor 90 which contains the spin axis controller algorithm corresponding to the functional diagram of Fig. 5. The output of the microprocessor 90 is coupled to a thruster command processor 92 which in turn is coupled to a spun driver unit 94, connected to control the axial thrusters 96A, 96B. Referring to the functional block diagram of Fig. 5, the state estimate propagation algorithm is executed every 10 minutes (corresponding to 2.5° of orbital motion). This calculation uses the previous estimates along with the specified drift rate to predict the present value of the states. The earth chord processor nominally uses 256 consecutive earth chord measurements to determine a low noise average chord width (Z aV g) which is then used in conjunction with the state predictions to determine a prediction error. This prediction error is multiplied by the estimator gains (K^K^Kt*,) and added to the respective predictions- to determine the new state estimates. A coordinate transformation of these state estimates from the body fixed frame (φ,ψ) to the sun fixed frame (U,V) is then performed before the maneuver logic is evaluated. The transformation of the state estimates to the sun fixed coordinates is to permit the use of constant axial thruster firing angles when performing spin axis attitude correction maneuvers. This eliminates the need for complex real time firing angle calculations, thus vastly simplifying the control procedures and processor demands. If attitude corrections are necessary, the magnitudes of the desired corrections in the sun fixed frame (U C ,V C ) are sent to the thruster command processor, and the appropriate axial thruster commands are executed if the SAC is enabled and no other thruster firings are in progress. The spin phase angle of the firings is controlled by a sun sensor pulse. Meanwhile, the corresponding corrections to the estimates of roll and yaw (Φfψ) are made and the results fed back to the state estimate propagation algorithm for use in the next cycle. The attached Appendix I lists the SAC commandable parameters. The ACP algorithm scaling is given in Appendix II. The SAC features two primary modes of operation (ETF, FTF), along with two special purpose executable sequences (IAT, QAT), which temporarily override the primary operating mode. The Error Threshold Firing Mode (ETF) is the simpler of the two operating modes. Attitude corrections are made whenever .a ground commandable threshold is exceeded. The threshold logic monitors the two attitude estimates in the sun fixed frame (U,V). When either of these estimates exceeds its threshold (U 0 or V 0 ) , a command to null the attitude in that particular direction is executed. For example, when V>V 0 , a V correction of -V 0 is performed. The thruster impulse is calibrated in orbit so that an integral number of thruster pulses will nominally result in the desired correction, -V 0 . This granularity equals the attitude correction resulting from a single thruster pulse and is a ground commandable parameter which is matched to actual thruster performance. ETF mode corrections can occur at any 2.5° point in the orbit. The only restrictions that prohibit corrections are the general restrictions for all SAC commanded corrections, such as SAC corrections will not be executed while other maneuvers are in progress or while the spacecraft is within an eclipse window (when the sun is not visible). The Fixed Time Firing Mode (FTF) involves a more complex logic which determines corrections according to both error magnitude and the spacecraft orbit location. Its objective is to restrict corrections to occur near the orbit node for the most efficient use of fuel. By making attitude corrections near the node, the small spacecraft velocity change which occurs will assist in north-south stationkeeping. Fig. 7 shows an example time history of the estimated spacecraft precession normal to the sunline during FTF mode operation. The attitude correction's attempt to maintain the precession within a ground specified band can be seen. Major corrections occur once a day at the node. For large precession rates where one correction per day is not sufficient, additional logic keeps the estimate within the band by a series of minor corrections made near the node in addition to the major correction at the node. The appropriate thresholds and timing are described as follows. If +/-θ 0 is allocated to be the spin axis attitude error budget, the delta angle (δ ) pe -axial thruster pulse and the estimator error (V ) are known. Then, the after node threshold (V a ) and the before node threshold (V****-*,) can be computed V b = abs (θ 0 - V) , V a ssign ( B )(V b - 2δ - ε y * δ) where B = estimated precession rate and εf = antenna flexure due to flexible motion coupling with nutation motion ( constant value) An extra δ angle was taken into account in computing V a s ince a s tep of δ atti tude correction in the angular momentum vector will induce a 2δ transien t motion in spacecraf t spin axis atti tude . The pa r ame te r M 0 designates the sidereal clock timing at which the FTF logic switches thresholds . The time paramete r for switching is def ined as the time from the orbit node in hours The use of M 0 in the firing logic is to ensure that all attitude corrections occur as near and syπimetric to the orbit node as possible. Thus, the FTF mode control logic (Fig.7) can be summarized: U maneuver ΔU C s U if U > ■ a UL at M. = 0 (orbit node) else lf^O V maneuver Δ V C = δ if SP * V > SP * V a and Q > M > 0 or Ivl > V b - nd M > M 0 else ΔV_ s V - V. + δ_ at M = 0 (orbit node) or ΔV c = 0 where δ p = SP*δ The percentage or the effectiveness fc^ of the SAC FTF mode attitude correction fuel consumption, which is also applicable to the north-south velocity change maneuver, is equal to cosine (ω 0 t 0 ) where ω 0 is the orbital rate and t 0 is the time from the orbit node crossing (i.e., ω 0 t 0 is orbital angle from the node). For an attitude trim occurring at t 0 = 0 hour (orbit node), ε = 100 percent, and at t 0 = 6 hours (antinode), ε= g percent. If an attitude trim occurs at t 0 greater than 6 hours, the fuel saving effectiveness becomes negative. Therefore, a SAC FTF mode attitude correction pulse firing at an orbital angle more than 90° away from the orbit node will cause the SAC fuel consumption to contribute to the north-south ΔV maneuver in a negative sense. Thus, it is desirable to design the FTF mode control logic suchthat t 0 is +/-6 hours from the orbit node crossing time. ' Two ground commanded control modes, Immediate Attitude Trim (IAT) and Quick Roll Estimate and Attitude Trim (QAT), are incorporated to null the spin axis attitude prior to the stationkeeping maneuvers. In IAT mode the SAC processor uses the latest onboard estimates Λ to execute an immediate attitude trim where U c = U and V c = v. In QAT mode the SAC initiates a quick roll attitude estimate using the average of 64 chord length measurements, rather than the normal 256 measurements in the ETF or FTF modes, and the sun frame components of the attitude error estimates are computed with yaw (ψ) equal to zero. Then, an immediate trim maneuver is executed in the same fashion as in IAT mode. The estimate of bias is not updated until 6 hours after QAT has occurred to eliminate the transient response. Several modes of SAC operation are commandable. In addition, the sidereal time, orbital angle, and state estimates may all be reset y ground command and nearly all the parameters discussed are programmable. The SAC mode command allows the ground to select ETF, FTF, QAT, or IAT operation of the SAC. The command is coded so that these mode commands are mutually exclusive. The default power-on operational mode is ETF. This mode of operation is also commanded by the processor at the completion of an IAT or QAT maneuver. The Spin Axis Controller logic diagram is attached hereto as Appendix III. An important function of the SAC mode command is the enable/disable of attitude corrections. If this bit is not set (i.e., attitude corrections are disabled), the estimators will continue to run, maneuver commands will be generated, but no thruster commands will be sent, and no thrusters will actually fire. The accumulated thruster pulse counts, SUMU and SUMV, are not changed if attitude corrections are disabled. To command thruster firings, the enable bit must be set and the platform despin controller (PDC) must be in platform pointing mode. Other bits in the SAC mode command control the selected input to the earth chord processor (ESPE-A or ESPE-B), the default or alternate gain selection for the estimators, the selection of the SAC spin rate to be a ground command constant or the auxiliary rate estimate from the PDC, the pulse width of commanded thruster maneuver (48 or 96 ms), and the selection of axial thruster (96A or 96B) to be used for attitude corrections. Default power-on conditions are: ESPE-A for earth chord processing, default estimator gains, SAC spin rate equal to a constant (30 rpm) , and a 96 ms pulse width with a zero thruster selection. Selection of thruster 96A or 96B is not mutually exclusive. If a zero thruster selection is sent, 96A will be used. If both thrusters are selected by ground command, 96B will be used. The variable reset command will reset the SAC estimates of roll, yaw, and bias to zero upon receipt of the command. A counter (36) which times out every 10 minutes is incremented each real time interrupt (RTI) of the processor. At time-out, the orbital angle,γ , is incremented by 2.5°, the estimate propagation occurs, and the earth chord processor 34 is reset if no maneuvers are currently in progress. The orbital angle, γ , is used in the coordinate transformation of roll and yaw to sun fixed coordinates. It may be set by ground command to any value between 0 and 360° with a resolution of 2.5°. When a new value of γ is set by ground command, the 10 minute timer 36 and the earth chord processor 34 are reset. All other SAC variables remain unchanged. It is possible to send a value of greater than 360° given the format of the serial command. The processor will use a value for γ which is the ground commanded value modulo 360°. A second counter 48 which times out every 5.61 minutes is also incremented each RTI of the processor. At time-out, the sidereal time (M) is incremented by its least significant bit (336.648 seconds), and the 5.61 minute counter is reset. The sidereal timer 48 continues to be incremented each 5.61 minutes until full scale is reached (23.94 hours) representing one sidereal day. It then sets the orbit node flag and begins to count again from zero. It may be set by ground command to any value between 0 and 23.94 hours. Upon receipt of such a command, the sidereal time will be set appropriately, and the 5.61 minute counter will be reset to zero. All other SAC variables remain unchanged. A more detailed exposition of the equations shown in the blocks 38, 40 and 42 of Fig. 5 is set fort as follows. As indicated in Fig. 5, the three state estimator 38, using earth sensor roll axis measurements, provides the state estimates for thruster control logic design. The three estimated states are roll angle (Φ ) , yaw angle (ψ), and sensor bias (b). Other variables (e-9-, the parameter A 0 in the precession rate model), which slowly vary compared to φ, ψ , or b, are treated as constants and may be updated periodically via ground command. The state estimator has the following simplified two state structure. Prediction N s *N-l * tt o l *N-l * o r sin ω o (N " l) l " *N = *N-l * "o^N-l * ^o : cos u o (N " l) τ H N = b N-1 Correction *N S" *N + K Φ (Z N * *N " 5 N> + Δ *N *N s " *N + κ ψ (Z N ~ *N " ->N> * Δ *N b N = b N + K b (Z N - T N - b N ) where l s estimator update time 5 . *N» *N» anc *- ^ s predic ed values of Φ, ψ, and b at sample N $^ » ϊ N , and b N = updated values of Φ, Φ, and b at sample N K φ , ^ and K b = estimator gains in Φ, *, and b respectively o Z - J S J - b-s j = measurement error at sample N Δ - S J and Δ$ N = control inputs in delta angles resulting from axial thruster firing This set of simplified estimator equations is applicable 5 when the estimator update time (t) is much shorter than the orbital period of 24 hours. The optimal estimator gains, K.(N), (N), and K D (N), can be determined by solving a Riccati equation. But the gains determined are time varying which is 0 difficult to implement in a microprocessor with limited memory size. Thus, it is preferable to seek a set of constant estimator gains, that provide sufficient speed of estimator error reduction and at the same time keep the estimator steady state error within acceptable bound. 5 Such a constant estimator gain matrix can be readily developed by methods known in the art. These may be selected to provide better spin axis attitude pointing in the steady state or by techniques which are more suitable for transient state operation, depending upon the choice of the user. With the use of digital simulation, it is possible to derive a set of constant estimator gains which may exhibit better spin axis pointing for the steady state operation and still provide sufficient speed for the estimator error reduction. The spin axis controller of the present invention is effective in controlling the trim maneuvers which are needed to maintain the spin axis of the spacecraft within acceptable limits of a preferred position. The spin axis controller accurately detects minor deviations of axis attitude from the orbit normal by measuring the difference between an earth chord observation from a selected earth sensor and a reference value. The. resulting measurement information is processed in accordance with the controller algorithm, including transformation of coordinate data to a sun fixed reference system, and is then utilized to control the firing of a selected axial thruster to correct the detected deviation. The timing of the firing of the thruster is selected to correspond with the rotational position and the orbital position of the satellite so that the firing of the selected thruster will act to both correct the spin axis attitude and adjust the orbital plane, if necessary, thereby economizing on the fuel needed to maintain the spacecraft in the proper attitude. To take care of possible contingencies, the spin axis controller is provided with a plurality of modes of operation. For larger precession errors which are in excess of a preset threshold level, the thruster can be commanded to fire in an error threshold firing mode so that the spin axis is drawn back within the threshold range. For more efficient overall fuel usage, the fixed time firing mode serves to restrict firings to occur near the orbital node crossing, such that the velocity impulses imparted by attitude corrections contribute favorably to the orbit inclination, i.e. north-south stationkeeping. Also, for minimizing the spin axis attitude errors prior to stationkeeping maneuvers, the spin axis controller includes two automatic sequences for nulling estimated attitude error: an immediate attitude trim mode in which the controller uses the latest onboard attitude estimate to execute an immediate attitude trim, and a quick roll estimate and attitude trim mode in which a quick roll estimate is utilized, followed by an immediate attitude trim maneuver. Either of the automatic sequences'can be executed from either onboard stored command or by direct ground command. These respective modes of operation of the spin axis controller are effective in assuring that the spacecraft can be controlled by the spin axis controller under the various operating conditions which may be encountered. Although there have been described above specific arrangements of an autonomous spin axis attitude controller in accordance with the invention for the purpose of illustrating the manner in which the invention may be used to advantage, it will be appreciated that the invention is not limited thereto. Accordingly, any and all modifications, variations or equivalent arrangements which may occur to those skilled in the art should be considered to be within the scope of the invention as defined in the annexed claims. APPENDIX I ∞MMANEABLE PARAMETERS Ccmmandable Default or Parameter Definitial Range Resolution Initial Value SPIN Spin rate, 150 to 270 0.0286 180.00 deg/sec Zβ Chord length 0 to 360 0.00549 0 bias, deg α,β Sensor cyclic +0.703125 0.00002145 0 bias, deg ΔP. Secular pre¬ +0.0027460 0.00002145 0 cession (per 10 min interval), deg Kφ Estimator gain, +1.0 1/256 0.5 roll Kψ Estimator gain. +2.0 1/128 0.5 yaw Kb Estimator gain, +1.0 1/256 0.148 sensor bias u 0 ,v 0 Firing thresholds +.00549° to 0.00002145 0.703 (ETF), deg +0.35156 Vb Firing threshold +.00549° to 0.00002145 0.703 (FTF), deg +0.35156 W Thruster pulse 48 or 96 width, ms Degree per +0.004° to pulse, deg +0.087890 0.00002145 0.02 N R Spin revolutions 255 between successive 1 15 pulses Reset value of Y 360 2.50 (orbit angle), deg M Reset value of 86,181,888 336.648 0 sidereal clock, sec Mo FTF threshold 43,090,944 336.648 43,090.944 switching, SAC Y ECLSTR Eclipse start 0 to 360 2.50 0 orbit angle, deg γ ECLENO Eclipse step 0 to 360 2.50 0 orbit angle, deg UBSTITUTE SHEET APPENDIX II ALGORITHM VARIABLE SCALING • Variab1• Range Resolu ion Bits Earth chord data 0 to 13 0.00002145 24 State esTiinates ).703125 0.00002145 16 Hβasuranβnts residual (e), deg Orbit angle, deg 0 to 360 2.5 8 Sidereal clock, sec 86,181.888 336.648 8 SUBSTITUTE SHEET APPENDIX III - 1 εcpMooε is SUBSTITUTE SHEET APPENDIX III - 2 SUBSTITUTE SHEET SUBSTITUTE SHEET

Claims: 24 WHAT IS CLAIMED IS: 1. An onboard control system for controlling the spin axis attitude of a spinning spacecraft comprising: earth sensing means for detecting deviation of the spin axis from a preferred attitude normal to the spacecraft orbit plane; signal processing means including a microprocessor for receiving signals from the sensing means indicative of detected deviation and developing corresponding thruster control signals; a pair of axial thrusters disposed on opposite sides of the spin axis for controlling spacecraft orbital position; a spun driver unit coupled to individually activate the firing of said thrusters in response to thruster control signals; and means for selectively applying thruster control signals to the spun driver unit to activate a selected thruster at a predetermined point in the rotation of said spacecraft which is effective to reduce said spin axis ' deviation. 25 2. The system of claim 1 wherein the earth sensing means comprise a first earth sensor having a field of view on one side of the spacecraft orbit plane and a second earth sensor having a field of view on the other side of said orbit plane, and wherein the signal processing means further include means for determining axis deviation from a comparison of signals from a selected earth sensor with a reference value corresponding to zero deviation. 26 3. The system of claim 1 wherein the earth sensing means comprise a first earth sensor having a field of view on one side of the spacecraft orbit plane and a second earth sensor having a field of view on the other side of said orbit plane, and wherein the signal processing means further include means for comparing signals from the first and second earth sensors to determine the degree of deviation. 27 4. The system of claim 2 further including a sun sensor for generating a sun pulse when the sensor scans the sun, and means for coordinating the thruster firing with the occurrence of the sun pulse. 28 5. The system of claim 4 wherein the thruster control signal applying means include means for applying said signals to generate a plurality of thruster pulses in the vicinity of the orbit node of the spacecraft. 29 6. The system of claim 5 wherein said plurality of thruster pulses are clustered on opposite sides of the orbit node and are symmetrical thereto. 30 7. The system of claim 1 wherein the signal processing means include means for averaging a predetermined plurality of signals from the sensing means and using the resulting average for developing thruster control signals. 31 8. The system of claim 7 further including means for interrupting said averaging means upon receipt of an activation signal from an associated ground control station in order to derive an average value from a number of sensing means signals which is less than said predetermined plurality and using said average value as a basis for nulling spin axis attitude deviation prior to stationkeeping maneuvers. 32 9. The system of claim 1 wherein the thrus ter control signal applying means include means for combining the thruster fir ings for spin axis correction with thruster firings for stationkeeping maneuvers. 33 10. The system of claim 1 further including means for compar ing the detected deviation signals with a predetermined threshold value and apply ing thrus ter control signals to activate a selected thruster only when said detected deviation signals exceed the threshold value. 34 11. The system of claim 10 further including means for limiting the thruster control signals to occur only near an orbit node. 35 12. The system of claim 10 wherein said threshold value is selectable in accordance with signals transmitted from an associated ground control station. 36 13. The system of claim 1 wherein the thruster control signal applying means further include means for activating a selected thruster in response to command signals from an associated ground control station to null spin axis deviation prior to the occurrence of stationkeeping maneuvers, based on onboard estimates of spacecraft reference data. 37 14. The method of controll ing the sp in axis attitude of a spinning spacecraft to reduce the deviation thereof from an orbit normal attitude comprising the steps of : sens ing the direct ion of the ear th f rom th e spacecraft to detect deviations of the spacecraft spin plane from coincidence with the spacecraf t orbital plane and deve loping s i gnals in accordance with detected deviation; generating signals for controlling the activation of spacecraf t thrusters in response to said detected deviation signals; and selectively applying the thruster control signals to activate a selected thruster at a predetermined point in the rotation of the spacecraft which is effective to reduce spin axis deviation. 38 15. The method of cla im 14 wherein the ear th sensing step includes scanning space with a first sensor having a field of view directed to one side of the spacecraft spin plane to develop a chord length signal corresponding to the sensing of the earth during a scan , and comparing said chord length signal with a reference signal corresponding to zero deviation in order to develop said detected deviation signals. 39 16. The method of claim 15 further including the steps of scanning space with a second sensor having a field of view on the opposite side of the spacecraft spin plane from that of the first sensor, and selecting chord length signals from one of the two sensors for comparing with the reference signal. 40 17. The method of claim 14 wherein the earth sensing step includes scanning space with first and second sensors having fields of view directed to opposide sides of the spacecraft spin plane to develop chord length signals from each sensor and comparing said chord length signals from the two sensors to develop said detected deviation signals. 41 18. The method of claim 15 further including the steps of scanning space to detect the sun and controlling the activation of the thrusters in accordance with the detection of the sun during said sun scan. 42 19. The method of claim 14 further including the s tep of l imi ting the application of thruster control pulses for spin axis attitude control to a period when the spacecraft is in the vicinity of an orbit node. 43 20. The method of claim 14 further including the step of apply ing the thruster control signals in combination with stationkeeping maneuvers. 44 21. The method of claim 15 further including the step of averaging a predetermined plurality of chord length signals and comparing the resulting average with said reference signal to determine deviation. 45 22. The method of claim 14 further including the step of comparing the detected deviation signals with a predetermined threshold value and activating a selected thruster only when said signals exceed the threshold value. 46 23. The method of claim 14 further including the step of activating a selected thruster a predetermined number of times in succession when the spacecraf t is in the vicinity of an orbit node. 47 24. The method of claim 14 wherein the activation of the selected thruster in the vicinity of an orbit node is controlled in accordance with the rotational phase of the spacecraf t to develop an element of orbit plane adjustment in addition to correction of the deviation of the spin axis. 48 25. The method of claim 14 in which the detected deviation signals are in body fixed frame coordinates and wherein the step of generating signals includes transforming the body fixed frame coordinates to sun fixed frame coordinates in order to control the thrusters from sun coordinate referenced signals. 49 26. The method of controlling the spin axis attitude of a spinning spacecraft to reduce the deviation thereof from an orbit normal attitude comprising the steps of: measuring the chord length of the earth within the field of view of a sensor rotating with the spacecraft and mounted at a predetermined cant angle relative to the orbit plane of the spacecraft; comparing signals derived from the measured chord length of the earth with a predetermined reference value to develop estimates of spin axis attitude error and sensor bias oriented in a first coordinate system which is referenced to the spacecraft fixed frame; transforming said developed estimates from said first coordinate system to a second coordinate system which is referenced to the sun fixed frame; activating a selected thruster at a predetermined position of the spacecraft in its orbit , about the earth; sensing the position of the sun during rotation of the spacecraft; and controlling the thruster activation to occur at a particular spin phase angle which is dependent upon the direction of the sun in said second coordinate system. 50 27. The method of claim 26 wherein the selected thruster activation at a predetermined position of the spacecraft occurs at an orbit node and is controlled to develop a major spin axis attitude correction in response to substantial spin axis attitude error measurements. 51 28. The method of claim 27 wherein the thruster activation comprises a series of minor axis attitude error corrections occurring during successive rotations of the spacecraft in the vicinity of said orbit node at spin phase angles determined in accordance with the sensed position of the sun. 52 29. The method of claim 26 further including the step of averaging a predetermined plurality of sensor signals measuring earth chord length to develop a composite signal for comparison with said predetermined reference value. 53 30. The method of claim 26 wherein the step of controlling thruster activation includes limiting the extent of thruster activation to that which is sufficient to maintain spin axis attitude error within a specified range. 54 31. The method of claim 29 wherein the step of averaging a predetermined plurality of earth chord length measurements includes interrupting said averaging at a predetermined point on receipt of a command from a ground station and using an average derived from a lesser number of chord length measurements for comparison with said predetermined reference value. 55 32. The method of claim 26 further including the step of providing an immediate attitude error correction upon receipt of a command signal from a ground station by using the most recently developed estimates of spin axis attitude error and sensor bias as transformed to the second coordinate system.

Inventors: CHEN CARL L, HUMMEL WILLIAM F JR, SLAFER LOREN I, SMAY JOHN W, CHEN, CARL, L., HUMMEL, WILLIAM, F., JR., SLAFER, LOREN, I., SMAY, JOHN, W.
Applicants: HUGHES AIRCRAFT CO, HUGHES AIRCRAFT COMPANY
Application year: 1986
Main classification: G05D1
IPC codes: B64G1, G05D1