WO2004032392A2 - Steering logic for control moment gyro system - Google Patents

Steering logic for control moment gyro system Download PDF

Info

Publication number
WO2004032392A2
WO2004032392A2 PCT/US2003/027312 US0327312W WO2004032392A2 WO 2004032392 A2 WO2004032392 A2 WO 2004032392A2 US 0327312 W US0327312 W US 0327312W WO 2004032392 A2 WO2004032392 A2 WO 2004032392A2
Authority
WO
WIPO (PCT)
Prior art keywords
cmg
gimbal
attitude control
further characterized
control according
Prior art date
Application number
PCT/US2003/027312
Other languages
French (fr)
Other versions
WO2004032392A3 (en
Inventor
Bong Wie
Original Assignee
Arizona Board Of Regents
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 Arizona Board Of Regents filed Critical Arizona Board Of Regents
Priority to AU2003294218A priority Critical patent/AU2003294218A1/en
Publication of WO2004032392A2 publication Critical patent/WO2004032392A2/en
Publication of WO2004032392A3 publication Critical patent/WO2004032392A3/en

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/28Guiding or controlling apparatus, e.g. for attitude control using inertia or gyro effect
    • B64G1/286Guiding or controlling apparatus, e.g. for attitude control using inertia or gyro effect using control momentum gyroscopes (CMGs)
    • 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

Definitions

  • This invention relates to control moment gyro (CMG) steering logic and more particularly to such logic that improves the escape from and/or avoids singularities. ,
  • Control moment gyros as applied to spacecraft attitude control and momentum management, have been extensively studied during the past three decades, Refs. 1 - 8, and more recently in Refs. 9 - 17. They have been successfully employed for a variety of space missions, such as the Skylab, the MIR station, and the International Space Station (ISS).
  • CMGs have never been used in commercial communications and imaging satellites because their higher torque capabilities have not been needed by most commercial and imaging satellites and also because CMGs are much more expensive and mechanically complex than reaction wheels.
  • a CMG contains a spinning rotor with large, constant angular momentum, but whose angular momentum vector direction can be changed with respect to the spacecraft by gimballing the spinning rotor.
  • Such a CMG is the Astrium Teldix DMG 15-45S.
  • the spinning rotor is mounted on a gimbal (or a set of gimbals), and torquing the gimbal results in a precessional, gyroscopic reaction torque orthogonal to both the rotor spin and gimbal axes.
  • the CMB is a torque amplification device because small gimbal torque input produces large control torque output on the spacecraft. Because the CMGs are capable of generating large control torques and storing large angular momentum over long periods of time, they have been employed for attitude control and momentum management of large space vehicles, such as the International Space Station (ISS).
  • ISS International Space Station
  • next-generation commercial imaging satellites may be equipped with CMGs because such satellites will require rapid rotational maneuverability for high- resolution images (Refs. 18 - 22).
  • a narrower field of vision will be required.
  • rapid slewing movement will be needed.
  • Pointing the entire spacecraft allows the body- fixed imaging system with a narrow field of view to achieve a higher definition and improves the resolution for its images.
  • the overall cost and effectiveness of such agile spacecraft is greatly affected by the average retargeting time.
  • mini-CMGs a low-cost attitude control system employing smaller and inexpensive CMGs, called mini-CMGs, is of current practical importance for developing future agile imaging spacecraft (Refs. 19 - 21) as well as small agile satellites (Ref. 22).
  • CMGs necessitates the development of CMG steering logic which generates the CMG gimbal rate commands in response to the CMG torque commands.
  • CMG steering logic which generates the CMG gimbal rate commands in response to the CMG torque commands.
  • CMG torque is available in all but one direction.
  • Approaches to the problem have been largely, but not entirely, successful.
  • the patented CMG steering logic of Refs. 13 and 14 is unable to escape the saturation singularities of certain CMG configurations, which can be problematic if CMG momentum desaturation is desired. Consequently, a new steering logic is desirable to overcome such a deficiency of the singularity-robust CMG steering logic. Furthermore, the new steering logic should provide an effective means of explicitly avoiding, instead of passing through, the internal elliptic singularities that are commonly encountered by most other pseudoinverse-based steering logic.
  • A -sinx, -c ⁇ cosx 2 sinx 3 c ⁇ cox ⁇ 4 [3] s ⁇ osx l s ⁇ cosx 2 s ⁇ cosx 3 s ⁇ cosx 4
  • Figs. 1(a) and (b) Two different types of the internal singularity for the pyramid array of four CMGs are illustrated in Figs. 1(a) and (b).
  • the hyperbolic singularity at (-90, 0, 90, 180) deg, which can be escaped by null motion, is shown in Fig. 1(b).
  • Fig. 1(b) Detailed analyses and discussions of CMG singularities can be found in Refs. 5 and 23.
  • 90 deg.
  • Fig. 1 is a 4-CMG configuration with two orthogonal pairs of scissored CMGs.
  • the Jacobian [3] becomes:
  • DGCMG double-gimbal control moment gyro
  • the rotor is suspended inside two gimbals and consequently the rotor momentum can be oriented on a sphere along any direction provided no restrictive gimbal stops.
  • DGCMGs For the different purposes of redundancy management and failure accommodation, several different arrangements of DGCMGs have been developed, such as three orthogonally mounted DGCMGs used in the Skylab and four parallel mounted DGCMGs employed for the International Space Station.
  • the total CMG momentum vector H ⁇ Hx,Hy,Hz) is expressed in the (x,y,z) axes as
  • the x-axis torque component is not a function of the outer gimbal /?,- , motions. Consequently, in Kennel's CMG steering law implemented on the International Space Station, the inner gimbal rate commands, .,. , are determined first for the commanded x-axis torque, then the outer gimbal rate commands, /?., for the commanded; - and z-axis torques.
  • Typical singularities of a system of four parallel double-gimbal CMGs are illustrated in Figs. 4(a) - (d).
  • the 4H saturation singularity is an elliptic singularity which cannot be escaped by null motion.
  • the 2H and OH singularities, shown in Fig. 4(b) and (c) respectively, are hyperbolic singularities which can be escaped by null motion.
  • a nonsingular configuration but with a zero momentum is also shown in Fig. 4(d).
  • the instantaneous torque vector, ⁇ generated by CMG gimbal motion x ⁇ can be defined as
  • x r x.
  • Most CMG steering laws determine the gimbal rate commands with some variant of the pseudoinverse of the form of the above expression [14].
  • A A ⁇ [AA ⁇ + 2I] _1 . [18]
  • the positive scalar ⁇ and the off-diagonal elements ⁇ i are to be properly selected such that A ⁇ ⁇ 0 for any nonzero constant ⁇ .
  • the singularity-robust inverse of the form (20) is based on the mixed, two-norm and weighted least-squares minimization although the resulting effect is somewhat similar to that of artificially misaligning the commanded control torque vector from the singular vector directions.
  • the singularity robust steering logic is based on the minimum two-norm, pseudoinverse solution, it does not explicitly avoid singularity encounters but it rather approaches and rapidly transits unavoidable singularities whenever needed.
  • the steering logic effectively generates deterministic dither signals when the system becomes near singular. Any internal singularities can be escaped for any nonzero constant torque commands using the singularity robust steering logic.
  • CMG steering logic with improved singularity robustness.
  • improved CMG steering logic that eliminates or avoids singularities such as occur in even the most robust prior CMG control is provided.
  • the CMG steering logic of this invention is based on the well-known mixed weighted two-norm and least-squares optimization solution with weighting matrices W and N.
  • W weighted pseudoinverse solution
  • N alone do not completely eliminate the deficiency; i.e. both W and N are required for completely resolving the singularity escape/avoidance problem.
  • the steering logic of the invention is simple to implement compared to many other explicit singularity-avoidance algorithms employing, for example, gradient methods, global search, optimization and/or null motion.
  • transient torque errors inherent to the steering logic of the invention during singularity escapes or transits will be acceptable in most practical situations in which precision pointing or tracking is not required during CMG momentum desaturation, CMG reconfiguration or spacecraft reorientation maneuvers.
  • Fig. 1 is a diagrammatic illustration of a CMG platform having a "pyramid" mounting arrangement of four single-gimbal CMGs;
  • Fig. 1(a) and (b) are vector diagrams of a 2H internal elliptical singularity (a) and 9H internal hyperbolic singularity (b) with the CMG arrangement of Fig. 1 ;
  • Fig. 2 is a diagrammatic illustration of the characteristics of two single gimbal CMGs with parallel gimbal axes
  • Fig. 2(a) is a diagrammatic illustration of a OH anti-parallel singularity of the CMG arrangement of Fig. 2;
  • Fig. 2(b) is a diagrammatic illustration of a 2H parallel singularity of a CMG arrangement of Fig 2;
  • Fig. 3 is a diagrammatic illustration of the characteristics of four parallel-mounted double-gimbal CMGs
  • Figs. 4(a) - (d) are diagrammatic representations of a four parallel-gimbal CMG attitude control in three singularity causing relationships (a) - (c) and in one relationship (d) not a singularity;
  • Fig. 5 is a series of plots vs. time of characteristics of a momentum desaturation singularity escape simulation for a two parallel single-gimbal CMG attitude control using the steering logic of the present invention
  • Fig. 10 is a series of plots vs. time of gimbal rate for singularity avoidance simulation with a gimbal rate limit of + 2 rad/sec;
  • Fig. 13 is a diagrammatic illustration of a three single-gimbal CMG array with parallel axes and Figs. 13(a) and (b) represent IH anti-parallel and 3H parallel singularities, respectively; and
  • Fig. 14 is a functional block diagram showing a control embodying the present invention to rotate a satellite in response to commanded rotation signal q c .
  • P is chosen such that A ⁇ P ⁇ ⁇ Ofor any torque command ⁇ , as follows:
  • ⁇ 0 sm( ⁇ t + ⁇ t ) [29] where the constant parameters, ⁇ 0 , ⁇ , ⁇ 0 , ⁇ , and ⁇ , need to be appropriately selected.
  • null motion can be generated at this hyperbolic singularity, this example demonstrates that the mere existence of null motion does not guarantee escape from a singularity (Refs. 5 and 23).
  • this type of internal singularity can be easily escaped by the singularity-robust steering logic of Refs. 13 and 14.
  • the commanded torque level as well as the CMG momentum is normalized as one without loss of generality. Thus it should be appropriately adjusted depending on the actual values of gimbal rate limit and CMG momentum magnitude.
  • the new steering logic also provides a simple yet effective way of passing through the OH internal hyperbolic singularity, but with a small transient torque error in u .
  • W ⁇ I is also needed to escape the saturation singularity, as shown in Fig. 8.
  • the x -axis torque component is not a function of the outer gimbal ⁇ t motions as can be seen in equation [10], Consequently the inner gimbal rate commands can be determined first for the commanded torque ⁇ x , then the outer gimbal rate commands for ( ⁇ y , ⁇ 2 ).
  • Kennel's CMG steering law distributes the CMG angular momentum vectors such that all inner gimbal angles are equal, which reduces the rate requirements on the outer gimbals as follows:
  • the desired satellite attitude q c at input 10 is presented as Euler angles, quaternions, Gibbs parameters, or some other convenient way of describing the attitude of a satellite.
  • the attitude enor 16 is used by the attitude control 18 to generate a desired body acceleration ⁇ c at output 20.
  • the modified pseudo inverse 36 is used to calculate the close loop gimbal rates ⁇ c at 42, according to the following:
  • the time varying off diagonal terms in P provide the escape of a control trajectory even when it is started a Moore-Penrose elliptical singularity.
  • the CMG anay 48 responds to the gimbal rate command 42.
  • the new gimbal angles ⁇ at 46 are used to generate at 44 the
  • the CMG array 48 movement produces torque h 50 on the satellite 52.
  • the satellite motion ⁇ 54 is detected by the attitude sensors 56.
  • Fig. 14 shows function blocks that may be implemented through hardware or software, preferably the latter in a computer based satellite control containing one or more signal processors programmed to produce output signals to control CMGs on the satellite as explained. Fundamentally the process is shown for a single signal path between two points, but it should be understood that single lines represent vector data which is 3 dimensional for the satellite attitude, attitude rate and torque, and a dimensional for the signals related to the n CMGs.
  • the terms "controller,” “control” and “control system” shall mean the calculating, signal producing and processing, and force, thrust or torque producing anangements that control spacecraft, robot or other object attitude in cooperation with a CMG anay.
  • controllers are known in the art and may include a dedicated control system computer or other data processing facility that calculates the control signals produced, the mechanical force, thrust or torque to be imparted using the steering logic programmed therein.

Abstract

A new steering logic based on a mixed weighted two-norm and least squares optimization solution overcomes a deficiency of being trapped in the momentum saturation singularities. The new steering logic utilizes an additional weighting matrix, (W), in conjunction with deterministic dither signals in a weighting matrix (V). The new steering logic also provides a simple means of avoiding troublesome, internal elliptic singularities which are commonly encountered by most other pseudoinverse-based steering logic. Various (CMG) systems, such as a typical pyramid array of four single-gimbal (CMGs), two and three parallel single-gimbal (CMG) configurations, and four parallel double-gimbal (CMGs) of the International Space Station are described as illustrative embodiments demonstrating the simplicity and effectiveness of the new steering logic.

Description

SINGULARITY ESCAPE/AVOIDANCE STEERING LOGIC FOR CONTROL
MOMENT GYRO SYSTEMS
Cross Reference to Related Application
This application claims priority from U.S. provisional patent applications Serial Nos. 60/407,334 filed August 28, 2002 and 60/493,923 filed August 8, 2003 in the name of Bong Wie entitled "Singularity Escape/Avoidance Steering Logic for Control Moment Gyro System," both of which are incorporated herein by reference.
Field of the Invention
This invention relates to control moment gyro (CMG) steering logic and more particularly to such logic that improves the escape from and/or avoids singularities. ,
Background
Control moment gyros (CMGs), as applied to spacecraft attitude control and momentum management, have been extensively studied during the past three decades, Refs. 1 - 8, and more recently in Refs. 9 - 17. They have been successfully employed for a variety of space missions, such as the Skylab, the MIR station, and the International Space Station (ISS). However, CMGs have never been used in commercial communications and imaging satellites because their higher torque capabilities have not been needed by most commercial and imaging satellites and also because CMGs are much more expensive and mechanically complex than reaction wheels. A CMG contains a spinning rotor with large, constant angular momentum, but whose angular momentum vector direction can be changed with respect to the spacecraft by gimballing the spinning rotor. Such a CMG is the Astrium Teldix DMG 15-45S. The spinning rotor is mounted on a gimbal (or a set of gimbals), and torquing the gimbal results in a precessional, gyroscopic reaction torque orthogonal to both the rotor spin and gimbal axes. The CMB is a torque amplification device because small gimbal torque input produces large control torque output on the spacecraft. Because the CMGs are capable of generating large control torques and storing large angular momentum over long periods of time, they have been employed for attitude control and momentum management of large space vehicles, such as the International Space Station (ISS). Four parallel mounted double-gimbal CMGs with a total weight of about 2400 LB and with a design life of 10 years are employed on the ISS. In particular, Kennel's steering law (Ref. 4) has been implemented on the double-gimbaled CMG system of the ISS.
Most of the next-generation commercial imaging satellites may be equipped with CMGs because such satellites will require rapid rotational maneuverability for high- resolution images (Refs. 18 - 22). In the interest of higher resolution images, a narrower field of vision will be required. To be able to image points of interest more widely distributed than points within the narrow angle of vision of the satellite, rapid slewing movement will be needed. Rather than sweeping the imaging system from side-to-side, the whole spacecraft body will turn rapidly. Pointing the entire spacecraft allows the body- fixed imaging system with a narrow field of view to achieve a higher definition and improves the resolution for its images. The overall cost and effectiveness of such agile spacecraft is greatly affected by the average retargeting time. Thus, the development of a low-cost attitude control system employing smaller and inexpensive CMGs, called mini-CMGs, is of current practical importance for developing future agile imaging spacecraft (Refs. 19 - 21) as well as small agile satellites (Ref. 22).
The use of CMGs necessitates the development of CMG steering logic which generates the CMG gimbal rate commands in response to the CMG torque commands. One of the principal difficulties in using CMGs for spacecraft attitude control and momentum management is the geometric singularity problem in which no control torque is generated for the commanded control torques along a particular direction. At such a singularity, CMG torque is available in all but one direction. The problem of overcoming singularities in CMG systems has previously been addressed. Approaches to the problem have been largely, but not entirely, successful.
The CMG singularity problem, studied previously by Margulies and Aubrun (Ref. 2) and Bedrossian et al. (Refs. 5 and 6), has been further examined recently in Ref. 23, incorporated herein by reference, to characterize and visualize the physical as well as mathematical nature of the singularities, singular momentum surfaces, and other singularity related problems. The steering logic that may be considered "baseline" steering logic for CMGs is pseudoinverse steering logic, discussed below. It is not free of singularities. A simple yet effective way of passing through, and also escaping from, any internal singularities, as applied to agile spacecraft pointing control, has been developed in Refs. 13 - 16, incorporated herein by reference. The CMG steering logic developed and patented by Wie, Bailey and Heiberg of Refs. 13 and 14 is mainly intended for typical reorientation maneuvers in which precision pointing or tracking is not required during reorientation maneuvers, and it fully utilizes the available CMG momentum space in the presence of any singularities. Although there are special missions in which prescribed attitude trajectories are to be "exactly" tracked in the presence of internal singularities, most practical cases will require a tradeoff between robust singularity transit/escape and the resulting, transient pointing errors.
Because the singularity-robust steering logic developed in Refs. 13 and 14 is based on the minimum two-norm, pseudoinverse solution, it does not explicitly avoid singularity encounters. Rather it approaches and rapidly transits unavoidable singularities whenever needed. It effectively generates deterministic dither signals when the system becomes near singular. Any internal singularities can be escaped for any nonzero constant torque commands using the singularity-robust steering logic.
However, the patented CMG steering logic of Refs. 13 and 14 is unable to escape the saturation singularities of certain CMG configurations, which can be problematic if CMG momentum desaturation is desired. Consequently, a new steering logic is desirable to overcome such a deficiency of the singularity-robust CMG steering logic. Furthermore, the new steering logic should provide an effective means of explicitly avoiding, instead of passing through, the internal elliptic singularities that are commonly encountered by most other pseudoinverse-based steering logic.
Control Moment Gyro Systems
There follows a summary of several representative CMG systems. Detailed descriptions of these systems can be found in Ref. 23. These CMG systems will be used in connection with a description of specific exemplary embodiments to demonstrate the simplicity and effectiveness of the new steering logic.
Pyramid Array of Four Single-Gimbal CMGs
For a typical pyramid mount of four single-gimbal CMGs with skew angle of β , shown in Fig. 1, the total CMG momentum vector is represented by H =
Figure imgf000004_0001
h2(x2)+ h3(x3)+ h4(x4)
[1]
Figure imgf000004_0002
where x. is the z'th gimbal angle, cβ ≡ cos/7, and sβ = sin ? . Then we obtain H = Ax [2] where x = (x, , x2 , x3 , x4 ) and A is the Jacobian matrix defined as
- c/Jcosx, sinx2 cβcosx3 -sinx4
A = -sinx, -cβcosx2 sinx3 cβcoxχ4 [3] sβ osxl sβcosx2 sβcosx3 sβcosx4
Equation (2) represents a linear mapping from x = (x1,x2,x3,x4) to H = (Hx,Hy,Hz). Two different types of the internal singularity for the pyramid array of four CMGs are illustrated in Figs. 1(a) and (b). The well known, troublesome elliptic singularity at (-90, 0, 90, 0) deg, which cannot be escaped by null motion, is illustrated in Fig. 1(a). The hyperbolic singularity at (-90, 0, 90, 180) deg, which can be escaped by null motion, is shown in Fig. 1(b). Detailed analyses and discussions of CMG singularities can be found in Refs. 5 and 23. When β = 90 deg. Fig. 1 is a 4-CMG configuration with two orthogonal pairs of scissored CMGs. The Jacobian [3] becomes:
0 sin 2 0 - sin;.
- sinxj 0 sin 3 0
A = cosx cosx2 cos 3 cosx
This special configuration is also of practical importance, as studied extensively in Refs. 1 and 23. Many other CMG configurations are some variants of this basic arrangement of two orthogonal pairs of two parallel CMGs.
Two and Three Parallel Single-Gimbal CMG Configurations
Two and three single-gimbal CMGs with parallel gimbal axes have been investigated in Refs. 1, 2 and 23 for two-axis control of a spacecraft. For a system of two CMGs without redundancy, the momentum vectors, Hj and H2 move in the (x,y) plane normal to the gimbal axis as shown in Fig. 2. For such "scissored" single-gimbal CMGs, the total CMG momentum vector is simply represented as
COSXj + cosx2
H = [4] smxj + smx2 where a constant unit momentum for each CMG is assumed.
Defining a new set of gimbal angles, and β , as follows: j + x2 X2 Xj a ,β [5]
2 " 2 where a is called the "rotation" angle and β the "scissor" angle (Ref. 1), one obtains
Figure imgf000006_0001
and
H = Ax [7] where H = [Hx, Hy \ x = (ά , β), and A is the Jacobian matrix defined as
Figure imgf000006_0002
For a system of three single-gimbal CMGs with parallel gimbal axes, the total CMG angular momentum vector is:
COSXj + cos 2 + cosx3
Η = sm x, + sin x2 + sm x3
and the Jacobian matrix for x = (x, , x2 , x3 ) becomes
-sm . smx2 - smx3
COSXj cosx2 cos 3
The singular momentum surfaces of this systems of three parallel single-gimbal CMGs are described by two circles as shown in Fig. 4. The internal hyperbolic singularity is shown in Fig. 4(a) and the external (saturation) elliptic singularity in Fig. 4(b). Detailed singularity analysis of this system can be found in Refs. 2 and 23.
Four Parallel Double-Gimbal CMGs
For a double-gimbal control moment gyro (DGCMG), the rotor is suspended inside two gimbals and consequently the rotor momentum can be oriented on a sphere along any direction provided no restrictive gimbal stops. For the different purposes of redundancy management and failure accommodation, several different arrangements of DGCMGs have been developed, such as three orthogonally mounted DGCMGs used in the Skylab and four parallel mounted DGCMGs employed for the International Space Station.
As shown by Kennel (Ref. 4), mounting of DGCMGs of unlimited outer gimbal angle freedom with all their outer gimbal axes parallel allows drastic simplification of the CMG steering law development in the redundancy management and failure accommodation and in the mounting hardware.
Such a parallel mounting arrangement of four double-gimbal CMGs with the inner and outer gimbal angles, at and βi of the z'th CMG, is shown diagra matically in Fig. 3.
The total CMG momentum vector H = {Hx,Hy,Hz) is expressed in the (x,y,z) axes as
Figure imgf000007_0001
where a constant unit momentum is assumed for each CMG. The time derivative of H becomes
∑ COSQ.,0.,
H = ∑(- sin a i cos βiάi - cos at sin βiβi [10] ∑(- sin i sin βiάi + cos α. cos /?, ?,- 1
Note that the x-axis torque component is not a function of the outer gimbal /?,- , motions. Consequently, in Kennel's CMG steering law implemented on the International Space Station, the inner gimbal rate commands, .,. , are determined first for the commanded x-axis torque, then the outer gimbal rate commands, /?., for the commanded; - and z-axis torques.
Typical singularities of a system of four parallel double-gimbal CMGs are illustrated in Figs. 4(a) - (d). The 4H saturation singularity is an elliptic singularity which cannot be escaped by null motion. The 2H and OH singularities, shown in Fig. 4(b) and (c) respectively, are hyperbolic singularities which can be escaped by null motion. A nonsingular configuration but with a zero momentum is also shown in Fig. 4(d).
The various CMG systems described in this section are used below to demonstrate the simplicity and effectiveness of the new singularity escape/avoidance logic according to this invention. Pseudoinverse Steering Logic
Ignoring the effect of spacecraft angular motion, the instantaneous torque vector, τ , generated by CMG gimbal motion x{ can be defined as
Figure imgf000008_0001
or, in matrix form, as
= ^ = ^ Xi = Ax [12] dt dx, i = l ' where A is a 3 x n Jacobian matrix and x = (x, ,...xn ) .
For the given control torque command , the gimbal rate command x , often referred to as the pseudoinverse steering logic, is then obtained as x = A+τ [13] where
A+ = AΓ (AAΓ )_1 [14]
This pseudoinverse is the minimum two-norm solution of the following constrained minimization problem: mini, 1, 2 x subject to Ax = τ [15] x where |x|| = xrx. Most CMG steering laws determine the gimbal rate commands with some variant of the pseudoinverse of the form of the above expression [14].
The pseudoinverse is a special case of the weighted minimum two-norm solution x = A+r where A+ = Q~x KT [AQ-X AT]X [16] of the following constrained minimization problem: minN 1, 2 x subject to Ax = τ [17] x "Q
2 where ||x|| = xrβx and Q = QT > 0. Below, the significance of choosing Q ≠ I, where I
is an identify matrix, is demonstrated. If rank (A) < m for certain sets of gimbal angles, or equivalently rank (AAr) < m, when A is an m x n matrix, the pseudoinverse does not exist, and it is said that the pseudoinverse steering logic encounters singular states. This singular situation occurs when all individual CMG torque output vectors are perpendicular to the commanded torque direction. Equivalently, the singular situation occurs when all individual CMG momentum vectors have extremal projections onto the commanded torque τ .
Since the pseudoinverse, A+ = Ar (AAΓ ) , is the minimum two-norm solution of gimbal rates subject to the constraint Ax = r , the pseudoinverse steering logic and all other pseudoinverse-based steering logic tend to leave inefficiently positioned CMGs alone causing the gimbal angles to eventually "hang-up" in anti-parallel singular arrangements. That is, they tend to steer the gimbals toward anti-parallel singular states. Despite this deficiency, the psuedoinverse steering logic, or some variant of pseudoinverse, is commonly employed for most CMG systems because of its simplicity for onboard, real-time implementation.
Nakamura, Y. and Hanafusa, H., in Ref. 24, describe a singularity robust steering logic in which:
A = AΓ [AAΓ + 2I]_1 . [18]
The singularity robust inverse steering logic, x = A*u , does not become singular; i.e., det [AAΓ + 0. This solves the large majority of singularity problems in CMG steering logic. However, x becomes zero if det (AAΓ J = 0 and if a control torque is commanded along the singular direction. The CMG system with the singularity robust inverse steering logic can then become trapped in the singular state.
Singularity Escape/Avoidance Steering Logic
A simple yet effective way of passing through and escaping any internal singularities was developed by Wie, Bailey and Heiberg (Refs. 13 - 16). Such a singularity-robust CMG steering logic was mainly intended for typical reorientation maneuvers in which precision pointing or tracking is not required during reorientation maneuvers, and it fully utilizes the available CMG momentum space in the presence of any singularities.
The singularity-robust steering logic of Refs. 13 - 14 has the following form: x = A#τ [19] where A = [ArA4 + Λ/]"1A7
[20] Ar[AAr + E]" and
Figure imgf000010_0001
The positive scalar λ and the off-diagonal elements εi are to be properly selected such that A τ ≠ 0 for any nonzero constant τ .
Note that there exists always a null vector of A# since rank(A# j< 3 for any λ and £,. when the Jacobian matrix A is singular. Consequently, a simple way of guaranteeing that A#r ≠ 0 for any nonzero constant τ command is to continuously modulate ε{ , for example, as follows: ε, = £0 sin(fitf + i) [22] where the amplitude εQ , the modulation frequency ω , and the phases φ , need to be appropriately selected. The scalar λ maybe adjusted as: λ = λ0 exp(- μ det(AAr )) [23] where λ0 and μ are constants to be properly selected. It is emphasized that the singularity-robust inverse of the form (20) is based on the mixed, two-norm and weighted least-squares minimization although the resulting effect is somewhat similar to that of artificially misaligning the commanded control torque vector from the singular vector directions. Because the singularity robust steering logic is based on the minimum two-norm, pseudoinverse solution, it does not explicitly avoid singularity encounters but it rather approaches and rapidly transits unavoidable singularities whenever needed. The steering logic effectively generates deterministic dither signals when the system becomes near singular. Any internal singularities can be escaped for any nonzero constant torque commands using the singularity robust steering logic.
However, a certain type of external saturation singularity cannot be escaped for any choice of λ and εi , as shown below and therefore a need exists for further improvement in CMG steering logic to escape and/or avoid all singularities, internal and external. Summary of the Invention
With this invention, there is provided CMG steering logic with improved singularity robustness. Particularly, improved CMG steering logic that eliminates or avoids singularities such as occur in even the most robust prior CMG control is provided. The CMG steering logic of this invention is based on the well-known mixed weighted two-norm and least-squares optimization solution with weighting matrices W and N. The proper use of W, in conjunction with the deterministic dither signals in N, eliminates a deficiency of the patented steering logic of Refs. 13 and 14. Either W (i.e. the weighted pseudoinverse solution) or N alone do not completely eliminate the deficiency; i.e. both W and N are required for completely resolving the singularity escape/avoidance problem.
Some internal elliptic singularities can be explicitly avoided or skirted by a proper selection of the weighting matrix W. However, the steering logic of the invention rather tends to approach and rapidly transit unavoidable singularities whenever necessary.
A trial-and-error, intuitive approach is used for choosing proper W's in several embodiments of the invention. The steering logic of the invention is simple to implement compared to many other explicit singularity-avoidance algorithms employing, for example, gradient methods, global search, optimization and/or null motion.
The transient torque errors inherent to the steering logic of the invention during singularity escapes or transits will be acceptable in most practical situations in which precision pointing or tracking is not required during CMG momentum desaturation, CMG reconfiguration or spacecraft reorientation maneuvers.
It will be seen that a new steering logic utilizing a mixed weighted two-norm and least squares optimization solution has been developed to overcome a deficiency of the previously patented steering logic. It can also be employed to explicitly avoid singularity encounters for which most other pseudoinverse-based steering logic either fail or must transit through. The practicality, effectiveness and simplicity of the new singularity escape/avoidance logic, even in the presence of gimbal-rate limit, are demonstrated for various CMG systems.
The above and further objects and advantages of the invention will be better understood from the following detailed description of preferred embodiments taken in consideration with the accompanying drawings. Brief Description of the Drawings
Fig. 1 is a diagrammatic illustration of a CMG platform having a "pyramid" mounting arrangement of four single-gimbal CMGs;
Fig. 1(a) and (b) are vector diagrams of a 2H internal elliptical singularity (a) and 9H internal hyperbolic singularity (b) with the CMG arrangement of Fig. 1 ;
Fig. 2 is a diagrammatic illustration of the characteristics of two single gimbal CMGs with parallel gimbal axes;
Fig. 2(a) is a diagrammatic illustration of a OH anti-parallel singularity of the CMG arrangement of Fig. 2; Fig. 2(b) is a diagrammatic illustration of a 2H parallel singularity of a CMG arrangement of Fig 2;
Fig. 3 is a diagrammatic illustration of the characteristics of four parallel-mounted double-gimbal CMGs;
Figs. 4(a) - (d) are diagrammatic representations of a four parallel-gimbal CMG attitude control in three singularity causing relationships (a) - (c) and in one relationship (d) not a singularity;
Fig. 5 is a series of plots vs. time of characteristics of a momentum desaturation singularity escape simulation for a two parallel single-gimbal CMG attitude control using the steering logic of the present invention; Fig. 6 is a series of plots vs. time of characteristics of a momentum desaturation singularity escape simulation for a three parallel single-gimbal CMG attitude control using the steering logic of the present invention with W = diag{ 1,2,3};
Fig. 7 is a series of plots vs. time of characteristics of a momentum desaturation singularity escape simulation for a four single-gimbal CMG (β = π/2) attitude control using the steering logic of the present invention with W = diag{l,2,3,4};
Fig. 8 is a series of plots vs. time of characteristics of a singularity transit simulation for a four single-gimbal CMG (β = 53.13 deg.) attitude control using the steering logic of the present invention with W = diag {1,1,1,1};
Fig. 9 is a series of plots vs. time of characteristics of a singularity avoidance simulation for a four parallel single-gimbal CMG (β = 53.13deg. ) attitude control using the steering logic of the present invention with W = diag {10"4,1, 1,1}; Fig. 10 is a series of plots vs. time of gimbal rate for singularity avoidance simulation with a gimbal rate limit of + 2 rad/sec;
Fig. 11 is a series of plots vs. time of characteristics of a 4H saturation singularity escape passing through the OH singularity of a four parallel double-gimbal CMG attitude control using the steering logic of the invention with W = diag {1,1 ,2,2, 1 , 1 ,2,2} ;
Fig. 12 is a series of plots vs. time of characteristics of a momentum desaturation singularity escape simulation for a four parallel double-gimbal CMG attitude control escaping the saturation singularity and also avoiding OH singularity using the steering logic of the present invention with W = diag{l,2,3,4,l,2,3,4}; Fig. 13 is a diagrammatic illustration of a three single-gimbal CMG array with parallel axes and Figs. 13(a) and (b) represent IH anti-parallel and 3H parallel singularities, respectively; and
Fig. 14 is a functional block diagram showing a control embodying the present invention to rotate a satellite in response to commanded rotation signal qc.
Detailed Description
In accordance with this invention, a new steering logic which is capable of escaping all types of singularities, is of substantially the form x = A*τ [24] where
Figure imgf000013_0001
s Q~lAT [AQ-1 AT + P"1 ]_1 [25]
≡ WAT[AWAT + v] where W ≡ Q~l and N ≡ P-1. For a 3 x n Jacobian matrix A, [ArEA + Q J is an n x n matrix and \AWAT + v\ is a 3 x 3 matrix. The singularity-robust inverse of the form [25] is the solution of the well-known, mixed two-norm and least-squares minimization problem:
Figure imgf000013_0002
where e = Ax - τ is the torque error vector. Because P and Q are positive definite matrices, ATPA = O and AWAT + V are always nonsingular. The weighting matrices, P and Q (equivalently, W and N), must be properly chosen: (i) to obtain acceptable levels of torque errors and gimbal rates, (ii) to escape any internal as well as external singularities, and (iii) to pass through singularities or avoid singularity encounters.
For example, P is chosen such that AτPτ ≠ Ofor any torque command τ , as follows:
Figure imgf000014_0001
and λ = λ0 exp(- μ det(AAr))
[28]
ε = ε0 sm(ωt + φt ) [29] where the constant parameters, λ0 ,μ,ε0, ω, and φ ,, need to be appropriately selected.
Because the condition ATPτ ≠ 0, which is a necessary condition, is not sufficient for escaping and/or avoiding singularities, the matrix Q also needs to be properly chosen, as follows:
Figure imgf000014_0002
Different values of W,- and/or nonzero off-diagonal elements are needed to be able to escape all types of singularities, including the external saturation singularities when CMG momentum desaturation is requested. Furthermore, a proper choice of W ≠ I provides an effective means for "explicitly" avoiding singularity encounters (to be demonstrated in the next section).
Several examples of demonstrating the significance of W ≠ I for escaping a certain type of external saturation singularities for which the previous singularity-robust steering logic of Wie, Bailey and Heiberg (Refs. 13 and 14) fails are presented in the discussion entitled "Exemplary Embodiments," below. Examples of avoiding an internal elliptic singularity for which the previous singularity-robust steering logic had to pass through are also presented. Exemplary Embodiments
In the following exemplary embodiments, the simplicity and effectiveness of the new steering logic will be demonstrated using various CMG systems, such as two and three parallel single-gimbal CMG configurations, a pyramid array of four single-gimbal CMGs and four parallel double-gimbal CMGs of the International Space Station.
Two Parallel Single-Gimbal CMGs
For a system of two parallel single-gimbal CMGs illustrated in Fig. 2 we have the CMG momentum vector H as
Figure imgf000015_0001
and the Jacobian matrix A as
Figure imgf000015_0002
where a is called the "rotation" angle and β the "scissor" angle.
For the internal anti-parallel singularity at x = ( ,β) = (θ,^/2), shown in Fig. 3(a), the Jacobian matrix becomes
Figure imgf000015_0003
with its singular vector u = (0,1); i.e., Aru-0. A null space vector of A such that Ax = 0 is: x =null(A) = (l,0).
Although null motion does exist from this singularity, this hyperbolic singularity cannot be escaped by null motion because the singular configuration ( (β = π 1 ) remains undisturbed during this null motion along the null manifold (or a degenerate trajectory). Although null motion can be generated at this hyperbolic singularity, this example demonstrates that the mere existence of null motion does not guarantee escape from a singularity (Refs. 5 and 23). However, this type of internal singularity can be easily escaped by the singularity-robust steering logic of Refs. 13 and 14.
For the saturation singularity at x = (a,β) = (θ,θ), shown in Fig. 3(b), the Jacobian matrix becomes 0 0
A = 2 0 with its singular vector u = (-1,0) and its null space vector x = null(A) = (0,1). This external elliptic singularity cannot be escaped by null motion because null motion does not exist in the vicinity of this momentum saturation singularity, as discussed in Refs. 5 and 23.
The singularity-robust inverse of Refs. 13 and 14 with λ - 0.1 and ε - 0.1 becomes
Figure imgf000016_0001
Note that x = [ά, /?) ≠ 0, but β - 0; i.e., the system remains singular while a changes. In fact, β = 0 for any torque command vector r and any εi . Consequently, the singularity- robust steering logic of Refs. 13 and 14 is unable to command a non-zero β for this special case.
However, this external singularity can be escaped by the new singularity escape/avoidance steering logic, proposed in this paper, of the form: x = A*τ where
Figure imgf000016_0002
As can be seen in Fig. 5, the momentum saturation singularity is escaped using the new steering logic with
Figure imgf000016_0003
Figure imgf000016_0004
where
£ = 0.1cost λ = 0.01 exp(- 10 det(AAr )) '
In Fig. 5, ψx,uy ) represent the actual torques generated by the CMG system for the commanded torques of ( τx , τy ) = (-1,0) for 0 < t < 3 sec. The CMG momentum is completely desaturated at t = 2 sec, but the desaturation torque is commanded until t - 3 sec to further explore any problem of encountering an internal singularity. (For all of the simulation results presented in this paper, the commanded torque level as well as the CMG momentum is normalized as one without loss of generality. Thus it should be appropriately adjusted depending on the actual values of gimbal rate limit and CMG momentum magnitude.) As can be seen in Fig. 5, the new steering logic also provides a simple yet effective way of passing through the OH internal hyperbolic singularity, but with a small transient torque error in u .
As discussed in Ref. 1, a steering logic based on the direct inverse of A can be obtained for this simple 2 x 2 problem, as follows: x = A~V [35] where
Figure imgf000017_0001
Then we obtain
sm.aτχ + cos aτy a [37] 2co&β Λ- λx
- cosaτx -s aτy β = [38] 2(sin/? + Λ2) where λ = λ0 exp(- μ det(AAr))sin β [39a] λ2 = λ0 exp(- μ det(AA ))cos β [39b] which are included in the steering law to avoid dividing by zero when the system becomes singular.
Three Parallel Single-Gimbal CMGs
For a system of three single-gimbal CMGs with parallel gimbal axes, W should be selected as Wl ≠ W2 ≠ W to escape the 3H saturation singularity at (xj , x2 , x3 ) = (0,0,0) with
xy) = (-l,0). Simulation results of (τxy) = ((-1,0) for 0 ≤ t ≤ 3 sec with W = diag {1,2,3} are shown in Fig. 6. Detailed singularity analyses of such a system of three parallel single-gimbal CMGs can be found in Ref. 23.
Pyramid Array of Four Single-Gimbal CMGs
First consider a special case (β = π/2) of the pyramid array of four single-gimbal CMGs. Such two orthogonal pairs of two parallel single-gimbal CMGs is of practical importance due to its simple arrangement of four single-gimbal CMGs for three-axis control applications (Refs. 1 and 23).
Similar to the case of three parallel single-gimbal CMGs, W≠ I is also needed to escape the saturation singularity, as shown in Fig. 8. For this simulation, W = diag {1,2,3,4} was used to escape the saturation singularity at (π/2,π/2,π/2,π/2 ) with (τxy2) =
(0,0,-1) for 0 < t ≤ A sec. It can be seen in Fig. 7 that an internal singularity is encountered at t =2 sec, but it is rapidly passed through.
In Fig. 8, simulation results are shown for a typical pyramid array of four single- gimbal CMGs (/? = 53.13 deg) with initial gimbal angles of (0,0,0,0) and a torque command of [τx , τy , τz ) = (1 ,0,0). For this simulation, W = diag {l,l,l,l} was used, i.e. the singularity robust steering logic of Refs. 13 and 14 was used for this simulation. It can be seen that the well-known internal elliptic singularity at (- /2,0,;r/2,θ) is encountered but it is successfully passed through. The inevitable transient torque errors during the singularity transit can be seen in Fig. 8. In Fig. 9, simulation results with W = diag |l0~4 ,1,1,1/ are shown for the previous case of a typical pyramid array of four single-gimbal CMGs (β = 53.13 deg) with initial gimbal angles of (0,0,0,0) and a torque command of \τx , τy , τz ) = (1 ,0,0). It can be seen that the well known internal elliptic singularity at - π 12,0,π 12,0) is not encountered. However, the inevitable transient torque errors caused by skirting such a troublesome, impassable elliptic singularity can also be seen in Fig. 9.
For the simulations shown in Figs. 8 and 9, the maximum gimbal rate of each CMG was explicitly imposed as ± 2 rad/sec. Such a gimbal rate saturation limit did not affect the singularity-avoidance performance of the proposed steering logic. The gimbal-rate time histories associated with the simulation shown in Fig. 9 are provided in Fig. 10.
Four Parallel Double-Gimbal CMGs
The new steering logic can also be employed for a system of four parallel double- gimbal CMGs described by Ax = τ where x = ( l, 2341234),
T = (τx , r, , τ. ) and A is a 3 x 8 Jacobian matrix which can be easily constructed from equation [10]. The saturation singularity escape capability of the new steering logic is demonstrated in Fig. 11 The simulation conditions for this case are: Initial gimbal angles: a( = β; = 0 for all i
Commanded torque: τxyz )= (θ,l, ) W = diag {1,1,2,2,1,1,2,2,}
Λ = 0.01exp(- det(AAr )) ε = O.lcost In Fig. 11 we have ax = a2 ≠ a3 = a4andβx = β2 ≠ β3 = β4. For this case with W = diag {1,1, 2,2, 1,1, 2,2} it escapes the saturation singularity but it passes through the OH singularity, as can be seen in Fig. 11. However, for W = I, this system remains singular at its singular momentum envelope, i.e. the patented steering logic of Refs. 13 and 14 is unable to escape the saturation singularity of this system of four double-gimbal CMGs even when CMG momentum desaturation is requested.
Furthermore, when we choose W = diag {1,2,3, 4, 1,2,3, 4}, it can be seen in Fig. 12 that it escapes the saturation singularity and that the OH singularity is completely avoided with no transient torque enors. The CMG momentum is completely desaturated at t = 4 sec to further explore any problem of encountering an internal singularity.
This example of a system of four double-gimbal CMGs demonstrates that it is feasible to completely avoid singularity encounters (with no transient torque enors) with a proper selection of W in conjunction with the deterministic dither signals generated in N.
For the parallel mounting arrangement of four double-gimbal CMGs, the x -axis torque component is not a function of the outer gimbal βt motions as can be seen in equation [10], Consequently the inner gimbal rate commands can be determined first for the commanded torque τx , then the outer gimbal rate commands for ( τy , τ2 ). In particular Kennel's CMG steering law distributes the CMG angular momentum vectors such that all inner gimbal angles are equal, which reduces the rate requirements on the outer gimbals as follows:
a, = ^^ + k( - a i ) [40]
^ COS Q.,. where k is called the inner gimbal angle distribution gain and a* is the desired inner gimbal angle for all CMGs chosen as
a - ∑<0°"'')»' . [41]
2_Jcos i The new steering logic can then be incorporated with Kennel's steering law to generate the outer gimbal rate commands β, using a reduced 2 x 4 Jacobian matrix without requiring
Kennel's virtual CMG for avoiding 2H-singularity encounters.
Simulation results also show that the new steering logic rapidly transits through the 2H singularity shown in Fig. 4(b).
In Fig. 14, an exemplary embodiment, like that of the incorporated-by-reference patent of Ref. 14, applied to satellite attitude control, the desired satellite attitude qc at input 10 is presented as Euler angles, quaternions, Gibbs parameters, or some other convenient way of describing the attitude of a satellite. At 14 it is compared with actual attitude qa at input 12 generated by the Inertial Measurement Unit (IMU) 56 or some other method for determining satellite attitude. The attitude enor 16 is used by the attitude control 18 to generate a desired body acceleration ώc at output 20. The torque command, hc at output 24 is calculated at 22 from the spacecraft inertial matrix, Js, nc = Js ώc. The modified pseudo inverse 36 is used to calculate the close loop gimbal rates δc at 42, according to the following:
Figure imgf000020_0001
Q and P being as defined above.
The time varying off diagonal terms in P provide the escape of a control trajectory even when it is started a Moore-Penrose elliptical singularity. The CMG anay 48 responds to the gimbal rate command 42. The new gimbal angles δ at 46 are used to generate at 44 the
1h Jacobean A = — [43]
2δ used in the pseudo inverse. The CMG array 48 movement produces torque h 50 on the satellite 52. The satellite motion ω 54 is detected by the attitude sensors 56.
It will be appreciated that Fig. 14 shows function blocks that may be implemented through hardware or software, preferably the latter in a computer based satellite control containing one or more signal processors programmed to produce output signals to control CMGs on the satellite as explained. Fundamentally the process is shown for a single signal path between two points, but it should be understood that single lines represent vector data which is 3 dimensional for the satellite attitude, attitude rate and torque, and a dimensional for the signals related to the n CMGs. As used herein the terms "controller," "control" and "control system" shall mean the calculating, signal producing and processing, and force, thrust or torque producing anangements that control spacecraft, robot or other object attitude in cooperation with a CMG anay. Such controllers are known in the art and may include a dedicated control system computer or other data processing facility that calculates the control signals produced, the mechanical force, thrust or torque to be imparted using the steering logic programmed therein.
Although prefened embodiments of the invention have been described in detail, it will be readily appreciated by those skilled in the art that further modifications, alterations and additions to the invention embodiments disclosed may be made without departure from the spirit and scope of the invention as set forth in the appended claims.
References
(I) Crenshaw, J., "2-SPEED, A Single-Gimbal Control Moment Gyro Attitude Control System" Paper No. 73-895, presented atAIAA Guidance and Control Conference, Key Biscayne, FL, August 20 - 22, 1973. (2) Margulies, G. and Aubrun, J. N., "Geometric Theory of Single-Gimbal Control Moment Gyro Systems," Journal of the Astronautical Sciences, Vol. 26, No. 2, 1978, pp. 159-191.
(3] Cornick, D. E., "Singularity Avoidance Control Laws for Single Gimbal Control Moment Gyros," Paper No. 79-1698, Proceedings of the AIAA Guidance and Control Conference, Boulder, CO, August 1979, pp. 20 - 33.
(4) Kennel, H. F., "Steering Law for Parallel Mounted Double-Gimb ailed Control Moment Gyros: Revision A," NASA TM-82390, January 1981.
(5) Bedrossian, N. S., Paradiso, J., Bergmann, E. V., and Rowell, D., "Redundant Single- Gimbal Control Moment Gyroscope Singularity Analysis," Journal of Guidance, Control and Dynamics, Vol. 13, No. 6, 1990, pp. 1096 - 1101.
(6) Bedrossian, N. S., Paradiso, J., Bergmann, E. V., and Rowell, D., "Steering Law Design for Redundant Single-Gimbal Control Moment Gyroscopes," Journal of Guidance, Control, and Dynamics, Vol. 13, No. 6, 1990, pp. 1083 - 1089 (also, N. S. Bedrossian's M.S. thesis at MIT, August 1987). (7) Vadali, S. R., Oh, H., and Walker, S., "Prefened Gimbal Angles for Single Gimbal Control Moment Gyroscopes," Journal of Guidance, Control, and Dynamics, Vol. 13, No. 6, 1990, pp. 1090 - 1095.
(8) Paradiso, J. A., "Global Steering of Single Gimballed Control Moment Gyroscopes Using a Direct Search," Journal of Guidance, Control, and Dynamics, Vol. 15, No. 5, 1992, pp. 1236 - 1244.
(9) Schaub, H. and Junkins, J. L., "Singularity Avoidance Using Null Motion and Variable-Speed Control Moment Gyros," Journal of Guidance, Control, and Dynamics, Vol. 23, No. 1, 2000, pp. 11 - 16.
(10) Ford, K. A. and Hall, C. D., "Singular Direction Avoidance Steering for Control Moment Gyros," Journal of Guidance, Control, and Dynamics, Vol. 23, No. 4, 2000, pp. 648
- 656.
(II) Heiberg, C. J., Bailey, D., and Wie, B., "Precision Spacecraft Pointing Using Single- Gimbal Control Moment Gyros with Disturbance," Journal of Guidance, Control, and Dynamics, Vol. 23, No. 1, 2000, pp. 77 - 85. (12) Wie, B., Space Vehicle Dynamics and Control, AIAA Education Series, AIAA, Washington, D.C., 1998, Chapter 7.
(13) Wie, B., Heiberg, C, and Bailey, D., "Singularity Robust Steering Logic for Redundant Single-Gimbal Control Moment Gyros," Journal of Guidance, Control, and Dynamics, Vol. 24, No. 5, 2001, pp. 865 - 872.
(14) Wie, B., Bailey, D., and Heiberg, C, "Robust Singularity Avoidance in Satellite Attitude Control," U.S. Patent No. 6,039,290, Issued March 21, 2000.
(15) Bailey, D., Heiberg, C. and Wie, B., "Continuous Attitude Control That Avoids CMG Anay Singularities," U.S. Patent No. 6,131,056, Issued October 10, 2000. (16) Wie, B., Heiberg, C, and Bailey, D., "Rapid Multi-Target Acquisition and Pointing Control of Agile Spacecraft," Journal of Guidance, Control, and Dynamics, Vol. 25, No. 1, 2002, pp. 96 - 104.
(17) Dominguez, J. and Wie, B., "Computation and Visualization of Control Moment Gyroscope Singularities," Paper No. 2002-4570, AIAA Guidance, Navigation, and Control Conference, Monterey, CA, August 5 - 8, 2002.
(18) Roser, X., Sghedoni, M., "Control Moment Gyroscopes (CMGs) and their Application in Future Scientific Missions," Proceedings of 3rd ESA International Conference on Spacecraft Guidance, Navigation, and Control Systems, ESTEC, Noordwijk, The Netherlands, Nov. 26 - 29, 1996, pp. 523 - 528. (19) Busseuil, J., Llibre, M., and Roser, X., "High Precision Mini-CMG's and their Spacecraft Applications," presented at AAS Guidance and Control-Advances in the Astronautical Sciences, Breckenridge, CO., February 1998.
(20) Salenc, C. and Roger, X., "AOCS for Agile Scientific Spacecraft with mini-CMGs," Proceedings of 4t l ESA International Conference on Spacecraft Guidance, Navigation, and Control Systems, ESTEC, Noordwijk, The Netherlands, Oct. 18 - 21, 1999, pp. 379 - 384.
(21) Defendini, A, Lagadec, K., Guay, P., and Griseri, G., "Low Cost CMG-based AOCS Designs," Proceedings of 4t l ESA International Conference on Spacecraft Guidance, Navigation, and Control Systems, ESTEC, Noordwijk, The Netherlands, Oct. 18 - 21, 1999, pp. 393 - 398. (22) Lappas, V. J., Steyn, W. J., and Underwood, C. I., "Attitude Control for Small Satellites Using Control Moment Gyros," presented at 52nd International Astronautical Federation Congress, Toulouse, France, Oct. 1 - 5, 2001. (23) Wie, B., "Singularity Analysis and Visualization of Single-Gimbal Control Moment Gyro Systems," AIAA Paper No. 2003-5658, AIAA Guidance, Navigation, and Control Conference, Austin, TX, August 11 - 14, 2003.
(24) Nakamura, Y. and Hanafusa, H., "Inverse Kinematic Solutions with Singularity Robustness for Robot Manipulator Control," Journal of Dynamic Systems, Measurement, and Control, Vol. 108, Sept. 1986, pp. 163 - 171.
Each of the foregoing references (1 - 24) is incorporated herein by reference.

Claims

I claim:
1. In a control moment gyro (CMG) attitude control of the kind comprismg an angular rate signal generator for each gyro in a CMG anay using pseudoinverse control law that uses a torque command signal and a Jacobian value for the angle of the gyro; the improvement characterized by a controller having:
(a) programming to calculate the Jacobian A from weighting matrices W and V, (b) programming enabling selection of terms of the matrix W, and (c) programming to compute a gimbal rate command x from a control torque command τ and the Jacobian developed from the matrices W and V.
2. The attitude control according to claim 1, further characterized by V being substantially of the form:
Figure imgf000025_0001
wherein λ = λ0 sin(ωt + φi ), and λQ,μ,ε0,ω and φt are selected constant parameters.
3. The attitude control accordmg to claim 1, further characterized by:
1 ε3 ε2
P_1 ≡ V = s_ 1 Sl > 0 ε2 εl 1 and wherein λ = λ0 exρ(- μ det(AAT )) , and εt = ε0 ~in(ωt + φt) , where λ0 ,μ,ε0, ω, and φt , are selected parameters.
4. The attitude control according to any one of claims 1, 2 or 3, further characterized by:
Figure imgf000026_0001
and different values of ωt and/or nonzero off-diagonal elements chosen to escape all types of singularities, including external saturation singularities.
5. The attitude control according to claim 4, further characterized by W≠ I, I being an identity matrix.
6. The attitude control according to claim 1, further characterized by the CMG anay comprising two parallel single-gimbal CMGs and wherein:
Figure imgf000026_0002
where
Figure imgf000026_0004
ε = 0.1 cost, and
Λ
Figure imgf000026_0003
7. The attitude control according to claim 1, further characterized by the CMG anay comprising three single-gimbal GMBs with parallel axes and W is selected as ψχ ≠ W2 ≠ W .
8. The attitude control according to claim 1, further characterized by the CMG anay comprising four single gimbal CMGs in a pyramid anay and:
W ≠ I, an identity matrice.
9. The attitude control according to claim 8, further characterized by substantially:
W = diag{l,2,3,4}.
1 10. The attitude control according to claim 8, further characterized by
2 substantially:
3 W = diag{l,l,l,l}.
1 11. The attitude control according to claim 8, further characterized by
2 substantially:
3 W = diag{10-4,1,U}.
1 12. The attitude control according to claim 1, further characterized by the CMG
2 anay comprising four double-gimbal CMGs and wherein:
3 x = (ax, 2,a3,a4x234), A τ = [τx , τy , τz ), and substantially 5 W = diag{l,l,2,2,l,l,2,2}.
1 13. The attitude control according to claim 1, further characterized by the CMG
2 anay comprising four double-gimbal CMGs and wherein:
3 x = (ax, 2, 3, 4x234), A τ = [τxy2 j, and substantially 5 W = diag{l,2,3,4,l,2,3,4}.
1 14. In an attitude control system of the kind that includes an anay of a plurality of
2 CMGs, means for generating an angular rate signal for each of the CMGs using
3 pseudoinverse control law that uses a torque command signal and a Jacobian value for the
4 angle of each gyro, the improvement characterized by:
5 (a) a singularity escaping and avoiding attitude controller developing
6 commands based upon a control logic and having;
7 (i) means for calculating x = A*τ,
8 where x = (x,,....xn), the Jacobian A =
Figure imgf000027_0001
where W ≡ Q"1 and V s P"1,
Figure imgf000028_0001
λ = λ0 exp(- μ det(AAT )) £• = ε0 sm(ωt + φt ), and cox λ λ λ λ co2 λ λ
Q'l ≡ W = > 0. λ λ co3 λ λ λ λ CO,
15. In a method of control moment gyro (CMG) attitude control using steering logic to develop steering inputs to the array and in which the steering logic is substantially of the form
x = AT , where x is the gimbal rate command, A is the Jacobian of CMG array angular momentum with respect to gimbal angle, τ is the control torque command; the improvement characterized by computing x = A# , where
Figure imgf000028_0002
r i- r- WAT[AWAT +v\ l where W ≡ Q~l and V ≡ P~ and W and P are weighting matrices chosen to avoid singularities.
16. The method of CMG attitude control according to claim 15, further characterized by V being substantially of the form
Figure imgf000028_0003
wherein λ = λQ ήxi(ωt + φι ), and λ0,μ,ε0,ω and φ. are selected constant parameters.
17. The method of CMG attitude control according to claim 15, further characterized by V being substantially of the form:
Figure imgf000029_0001
wherein λ = λQ exp(- μ det(AAr )), and ε i = εQ sin(ωt + φi).
18. The method of CMG attitude control according to any one of claims 15, 16 and 17, further characterized by W being substantially of the form:
Figure imgf000029_0002
and values of ω and/or nonzero off-diagonal elements being chosen to escape singularities.
19. The method of CMG attitude control according to claim 15, further characterized by W≠ I, I being an identify matrix.
20. The method of CMG attitude control according to claim 15, wherein the CMG anay comprising two parallel single-gimbal CMGs and further characterized by substantially:
Figure imgf000029_0003
where
Figure imgf000029_0004
ε = 0.\ cost, and
A = 0.01 exp (-10 det(AA T1)\)
21. The method of CMG attitude control according to claim 15, further characterized by the CMG anay comprising three single-gimbal GMBs with parallel axes and W is selected as WX ≠ W2 ≠ W3.
22. The method of CMG attitude control according to claim 15, further characterized by the CMG anay comprising four single gimbal CMGs in a pyramid anay and:
W ≠ I, an identity matrice.
23. The method of CMG attitude control according to claim 22, further characterized by substantially:
W = diag{ 1,2,3,4}.
24. The method of CMG attitude control according to claim 22, further characterized by substantially:
W = diag{l,l,l,l}.
25. The method of CMG attitude control according to claim 22, further characterized by substantially:
W = diag{10-4,U,1}.
26. The method of CMG attitude control according to claim 15, further characterized by the CMG anay comprising four double-gimbal CMGs and wherein substantially:
x = (ax,a2,a3,a4x234), τ = {τxy2 ), and substantially W = diag{l,l,2,2,l,l,2,2}.
27. The method of CMG attitude control according to claim 15, further characterized by the CMG anay comprising four double-gimbal CMGs and wherein substantially:
x = {ax,a2,a3,a4x234), τ =
Figure imgf000030_0001
and substantially W = diag{l,2,3,4,l,2,3,4}.
28. In a method of attitude confrol system of the kind that includes providing an anay of a plurality of CMGs, means for generating an angular rate signal for each of the CMGs using pseudoinverse control law that uses a torque command signal and a Jacobian value for the angle of each gyro, the improvement characterized by:
(a) developing singularity escaping and avoiding attitude commands based upon a control logic including;
(i) calculating x = A#τ, where x = (xj ,....x„ ), the Jacobian A# =
Figure imgf000031_0001
≡ WAT[AWAT +v\λ where W ≡ Q"1 and V ≡ P"1,
Figure imgf000031_0002
λ = λ0 exp(- μ det(AAT )), £,- = ε0 sin(ωt + φt ), and
Figure imgf000031_0003
PCT/US2003/027312 2002-08-28 2003-08-28 Steering logic for control moment gyro system WO2004032392A2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
AU2003294218A AU2003294218A1 (en) 2002-08-28 2003-08-28 Steering logic for control moment gyro system

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
US40733402P 2002-08-28 2002-08-28
US60/407,334 2002-08-28
US49392303P 2003-08-08 2003-08-08
US60/493,923 2003-08-08

Publications (2)

Publication Number Publication Date
WO2004032392A2 true WO2004032392A2 (en) 2004-04-15
WO2004032392A3 WO2004032392A3 (en) 2005-03-03

Family

ID=32073290

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2003/027312 WO2004032392A2 (en) 2002-08-28 2003-08-28 Steering logic for control moment gyro system

Country Status (3)

Country Link
US (1) US6917862B2 (en)
AU (1) AU2003294218A1 (en)
WO (1) WO2004032392A2 (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006020314A1 (en) * 2004-07-30 2006-02-23 Honeywell International Inc. Method and system for optimizing torque in a cmg array
WO2006085996A1 (en) * 2004-07-23 2006-08-17 Honeywell International Inc. Method and system for cmg array singularity avoidance
US7805226B2 (en) 2006-09-29 2010-09-28 Honeywell International Inc. Hierarchical strategy for singularity avoidance in arrays of control moment gyroscopes
US8589814B2 (en) 2010-04-16 2013-11-19 Honeywell International Inc. System and method for visual presentation of information in a process control system
CN104090489A (en) * 2014-07-02 2014-10-08 中国科学院长春光学精密机械与物理研究所 Flexible agile satellite attitude maneuver rolling optimization control method
US9191843B2 (en) 2013-06-12 2015-11-17 Honeywell International Inc. Apparatus and method for measuring and reporting redundant wireless connectivity over time
US9240164B2 (en) 2013-02-27 2016-01-19 Honeywell International Inc. Apparatus and method for providing a pan and zoom display for a representation of a process system
US9609524B2 (en) 2014-05-30 2017-03-28 Honeywell International Inc. Apparatus and method for planning and validating a wireless network
US9977413B2 (en) 2013-03-11 2018-05-22 Honeywell International Inc. Apparatus and method for managing open windows in a graphical display for a representation of a process system
CN109189086A (en) * 2018-07-27 2019-01-11 中国人民解放军战略支援部队航天工程大学 A kind of spacecraft attitude based on magnetic suspension control sensitivity gyro and vibration integrated control method

Families Citing this family (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2846107B1 (en) * 2002-10-21 2005-06-24 Cit Alcatel DEVICE FOR CONTROLLING ATTITUDE OF A GEOSTATIONARY SATELLITE
KR101118888B1 (en) * 2004-06-02 2012-03-19 록웰 콜린스 컨트롤 테크놀로지스, 인코포레이티드 Systems and methods for controlling dynamic systems
KR100552583B1 (en) 2004-11-16 2006-02-15 한국과학기술원 Control moment gyros having improved singularity problem and method for steering the same
US7370833B2 (en) * 2005-10-20 2008-05-13 Honeywell International Inc. Method and system for determining a singularity free momentum path
US7693619B2 (en) * 2005-11-30 2010-04-06 Honeywell International Inc. Method and system for controlling sets of collinear control moment gyroscopes with offset determination without attitude trajectory of spacecraft
KR100778098B1 (en) 2006-07-26 2007-11-22 한국항공우주연구원 Moment gyros'cluster for spacecraft three axis attitude control
US7627404B2 (en) * 2007-04-13 2009-12-01 The Boeing Company Singularity escape and avoidance using a virtual array rotation
JP5125542B2 (en) * 2008-01-22 2013-01-23 三菱電機株式会社 Attitude control device and attitude control method for artificial satellite
US8209070B2 (en) * 2008-12-17 2012-06-26 Honeywell International Inc. Methods and systems for efficiently orienting an agile vehicle using a gyroscope array
WO2010135421A2 (en) * 2009-05-19 2010-11-25 University Of Florida Research Foundation, Inc. Attitude control system for small satellites
US8346538B2 (en) 2010-05-25 2013-01-01 Honeywell International Inc. Methods and systems for reducing angular velocity using a gyroscope array
US9038958B1 (en) * 2012-05-29 2015-05-26 United States Of America As Represented By The Secretary Of The Navy Method and apparatus for contingency guidance of a CMG-actuated spacecraft
US8880246B1 (en) * 2012-08-22 2014-11-04 United States Of America As Represented By The Secretary Of The Navy Method and apparatus for determining spacecraft maneuvers
US9567112B1 (en) * 2013-06-27 2017-02-14 The United States Of America, As Represented By The Secretary Of The Navy Method and apparatus for singularity avoidance for control moment gyroscope (CMG) systems without using null motion
US9817371B2 (en) * 2014-10-22 2017-11-14 Honeywell International Inc. Single-event transient feedback disturbance suppression system
US10407187B2 (en) * 2015-06-22 2019-09-10 Mitsubishi Electric Corporation Device for controlling attitude of spacecraft and method for calculating CMG gimbal angle
FR3045847B1 (en) * 2015-12-17 2018-01-05 Centre National D'etudes Spatiales (Cnes) METHOD FOR CONTROLLING A REDUNDANT SYSTEM
WO2019204324A1 (en) * 2018-04-16 2019-10-24 The Regents Of The University Of California Linear and angular position stabilization and control of an underwater robotic system
KR102280790B1 (en) * 2019-12-20 2021-07-22 한국항공우주연구원 Satellite having reaction wheel having arrangement structure and direction of rotation for preventing zero-crossing
FR3107953B1 (en) * 2020-03-03 2022-03-11 Airbus Defence & Space Sas Method for acquiring images of a terrestrial area by a spacecraft

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5435503A (en) * 1993-08-27 1995-07-25 Loral Vought Systems Corp. Real time missile guidance system
US6681649B2 (en) * 2002-04-03 2004-01-27 Honeywell International Inc. Inertial control and measurement system
US6682019B2 (en) * 2002-04-04 2004-01-27 Honeywell International Inc. Minimum energy wheel configurations for energy storage and attitude control

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6131056A (en) * 1998-03-16 2000-10-10 Honeywell International Inc. Continuous attitude control that avoids CMG array singularities
US6039290A (en) 1998-03-16 2000-03-21 Honeywell Inc. Robust singularity avoidance in satellite attitude control
US6128556A (en) * 1998-03-16 2000-10-03 Honeywell International Inc. CMG control based on angular momentum to control satellite attitude
US6047927A (en) * 1998-03-16 2000-04-11 Honeywell Inc. Escaping singularities in a satellite attitude control
US6648274B1 (en) * 2002-04-12 2003-11-18 David A. Bailey Virtual reaction wheel array

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5435503A (en) * 1993-08-27 1995-07-25 Loral Vought Systems Corp. Real time missile guidance system
US6681649B2 (en) * 2002-04-03 2004-01-27 Honeywell International Inc. Inertial control and measurement system
US6682019B2 (en) * 2002-04-04 2004-01-27 Honeywell International Inc. Minimum energy wheel configurations for energy storage and attitude control

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006085996A1 (en) * 2004-07-23 2006-08-17 Honeywell International Inc. Method and system for cmg array singularity avoidance
US7246776B2 (en) 2004-07-23 2007-07-24 Honeywell International, Inc. Method and system for CMG array singularity avoidance
JP2008507446A (en) * 2004-07-23 2008-03-13 ハネウェル・インターナショナル・インコーポレーテッド Method and system for avoiding CMG array specificity
WO2006020314A1 (en) * 2004-07-30 2006-02-23 Honeywell International Inc. Method and system for optimizing torque in a cmg array
US7014150B2 (en) 2004-07-30 2006-03-21 Honeywell International Inc. Method and system for optimizing torque in a CMG array
US7805226B2 (en) 2006-09-29 2010-09-28 Honeywell International Inc. Hierarchical strategy for singularity avoidance in arrays of control moment gyroscopes
US8589814B2 (en) 2010-04-16 2013-11-19 Honeywell International Inc. System and method for visual presentation of information in a process control system
US9240164B2 (en) 2013-02-27 2016-01-19 Honeywell International Inc. Apparatus and method for providing a pan and zoom display for a representation of a process system
US9977413B2 (en) 2013-03-11 2018-05-22 Honeywell International Inc. Apparatus and method for managing open windows in a graphical display for a representation of a process system
US9191843B2 (en) 2013-06-12 2015-11-17 Honeywell International Inc. Apparatus and method for measuring and reporting redundant wireless connectivity over time
US9609524B2 (en) 2014-05-30 2017-03-28 Honeywell International Inc. Apparatus and method for planning and validating a wireless network
CN104090489A (en) * 2014-07-02 2014-10-08 中国科学院长春光学精密机械与物理研究所 Flexible agile satellite attitude maneuver rolling optimization control method
CN104090489B (en) * 2014-07-02 2016-12-07 中国科学院长春光学精密机械与物理研究所 A kind of flexible agile satellite attitude maneuvers rolling optimization control method
CN109189086A (en) * 2018-07-27 2019-01-11 中国人民解放军战略支援部队航天工程大学 A kind of spacecraft attitude based on magnetic suspension control sensitivity gyro and vibration integrated control method

Also Published As

Publication number Publication date
WO2004032392A3 (en) 2005-03-03
US6917862B2 (en) 2005-07-12
US20040111194A1 (en) 2004-06-10
AU2003294218A1 (en) 2004-04-23
AU2003294218A8 (en) 2004-04-23

Similar Documents

Publication Publication Date Title
US6917862B2 (en) Singularity escape/avoidance steering logic for control moment gyro systems
Wie Singularity escape/avoidance steering logic for control moment gyro systems
Wiśniewski et al. Fully magnetic attitude control for spacecraft subject to gravity gradient
Wisniewski Linear time-varying approach to satellite attitude control using only electromagnetic actuation
Chesi et al. Aerodynamic three-axis attitude stabilization of a spacecraft by center-of-mass shifting
US6360996B1 (en) Steering control for skewed scissors pair CMG clusters
Narkiewicz et al. Generic model of a satellite attitude control system
US7370833B2 (en) Method and system for determining a singularity free momentum path
Kusuda et al. Feedback control with nominal inputs for agile satellites using control moment gyros
Meng et al. Failure-tolerant control for small agile satellites using single-gimbal control moment gyros and magnetic torquers
Yoon Spacecraft attitude and power control using variable speed control moment gyros
Kojima et al. Spacecraft line of sight maneuver control using skew-arrayed two single-gimbal control moment gyros
Evain et al. Satellite attitude control with a six-control moment gyro cluster tested under microgravity conditions
Reijneveld et al. Attitude control system of the Delfi-n3Xt satellite
Karpenko et al. Flight implementation of pseudospectral optimal control for the TRACE space telescope
Karpenko et al. Implementation of shortest-time maneuvers for generic CMG steering laws
Courie et al. Simulation of satellite attitude control using Single Gimbal Control Moment Gyro (SGCMG) system
Hao et al. A practical rapid attitude maneuver control system using control moment gyros for microsatellite tsubame
Guo et al. Integrated control of attitude maneuver and vibration suppression using pyramid-type SGCMGs
Jian et al. Research on attitude control method of agile satellite based on variable structure control algorithm
Ashok et al. Cmg configuration and steering approach for spacecraft rapid maneuvers
Nanamori et al. Steering law of control moment gyros using optimization of initial gimbal angles for satellite attitude control
Lappas et al. Combined Singularity Avoidance for Redundant Control Moment Gyro (CMG) Clusters
Salleh et al. Three-axis attitude control performance of a small satellite using control moment gyroscope
Inumoh et al. Three-axis attitude control of a satellite in zero momentum mode using a tilted wheel methodology

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A2

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SY TJ TM TN TR TT TZ UA UG UZ VC VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A2

Designated state(s): GH GM KE LS MW MZ SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LU MC NL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
122 Ep: pct application non-entry in european phase
NENP Non-entry into the national phase

Ref country code: JP

WWW Wipo information: withdrawn in national office

Country of ref document: JP