Mathematical Problems in Engineering
Volume 2015 (2015), Article ID 278145, 11 pages
http://dx.doi.org/10.1155/2015/278145

Research Article

The Coupling Characteristic Investigation of Double-Gimbal Magnetically Suspended Control Moment Gyro Used on Agile Maneuver Spacecraft

Peiling Cui,1,2 Jian Cui,1,2 Qian Yang,¹ and Shiqiang Zheng1,2

¹School of Instrumentation Science and Optoelectronics Engineering, Beihang University, Beijing 100191, China
²Science and Technology on Inertial Laboratory, Beijing 100191, China

Received 4 November 2014; Revised 27 January 2015; Accepted 10 February 2015

Academic Editor: Paulo Batista Gonçalves

Copyright © 2015 Peiling Cui et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Double-gimbal magnetically suspended CMG is a novel attitude control actuator for the agile maneuver spacecraft. Taking the double-gimbal magnetically suspended control moment gyro used on agile maneuver spacecraft as the research object, the dynamic model of the magnetically suspended rotor, the inner gimbal, and the outer gimbal of double-gimbal magnetically suspended control moment gyro is built. The nonlinear coupling characteristic between the rotor, the gimbal, and the spacecraft is given. It can be seen that the motion of magnetically suspended rotor does not only rely on magnetic bearing force but also suffer from the influence of gimbal servo system and spacecraft motion. The coupling torque includes the gyro coupling torque and the inertial coupling torque. The work in this paper provides the foundation for further studies.

1. Introduction

Double-gimbal magnetically suspended control moment gyro (CMG) is a novel solution to realize high precision, long life, and attitude maneuver control of spacecraft by incorporating the advantages of magnetic bearing with double gimbals [1, 2]. Compared with single-gimbal magnetically suspended CMG, double-gimbal magnetically suspended CMG has more than one degree of freedom, and the singularity is smaller. Therefore, it is believed that double-gimbal magnetically suspended CMG can implement the rapid attitude maneuver with high efficiency [3, 4]. By the rotation of inner gimbal and outer gimbal, the angular momentum direction of magnetically suspended high-speed rotor is changed, and then the gyro torque is produced. For magnetic bearing, there are five more degrees of freedom than mechanical stator, including three translation motions and two radial rotation motions [5, 6].

Double-gimbal magnetically suspended CMG is composed of high-speed magnetically suspended rotor system, inner gimbal servo system, and outer gimbal servo system [7]. In the agile maneuver spacecraft using double-gimbal magnetically suspended CMG as actuator, there is strong coupling torque on the magnetic bearing resulting from the motion of spacecraft and gimbal, and then the high-speed rotor displacement is increased. When this phenomenon is serious, the rotor will collide with the protecting bearing, and the stability of magnetically suspended rotor will be reduced and even lost. Moreover, the spacecraft maneuver and the rotor radial motion will bring disturbance to the gimbal motion. The output torque precision of double-gimbal magnetically suspended control moment gyro will be influenced. Finally, the spacecraft attitude control precision will be reduced when the gimbal motion and rotor radial motion are large. It is important to build the dynamic model of magnetically suspended rotor, inner gimbal, outer gimbal, and spacecraft when using double-gimbal magnetically suspended control moment gyro as the spacecraft attitude control actuator.

For spacecraft dynamic modeling based on mechanical CMG, there is some research result. In [8], the mechanical CMG model is given, and the vector-form output torque of CMG in spacecraft body frame is provided. In [9], the dynamic model of one mechanical double-gimbal CMG is built using Lagrange equation, and the adaptive feedback control law is put forward. In [10], the inertial torque of inner gimbal and outer gimbal is considered, and the spacecraft attitude dynamic model with a cluster of double-gimbal CMGs is built. For spacecraft dynamic modeling based on magnetically suspended CMG, there is little research result. In [11], the spacecraft dynamic model using single-gimbal magnetically suspended CMG is built. It is aimed at large spacecraft, the spacecraft angular velocity is limited to 10⁻² rad/s, and the effect of spacecraft maneuver is not involved. As for double-gimbal magnetically suspended CMG, investigation mainly focuses on single part of the magnetically suspended CMG, such as rotor system [12, 13] and gimbal servo system [14]. Models are built separately, omitting the coupling characteristic between different parts. Besides, researchers do not consider maneuver of the spacecraft.

In this paper, taking double-gimbal magnetically suspended CMG used for attitude control of agile maneuver spacecraft as the research object, the dynamic model of magnetically suspended rotor, inner gimbal, and outer gimbal of double-gimbal magnetically suspended CMG is built. The nonlinear coupling characteristic between rotor, gimbal, and spacecraft is investigated.

2. The Working Principal of Double-Gimbal Magnetically Suspended CMG

Double-gimbal magnetically suspended CMG consists of three subsystems, including magnetically suspended high-speed rotor system, inner gimbal rate servo system, and outer gimbal rate servo system. The angular momentum is provided by magnetically suspended high-speed rotor. The torque is output by the direction changing of angular momentum, which will adjust the spacecraft attitude.

Double-gimbal magnetically suspended CMG prototype is shown in Figure 1. The gyro room and its subsidiary control system is called the magnetically suspended high-speed rotor system, including the magnetic bearing system and the motor for the high-speed rotor. The gyro room is equipped with the components, including the motor of rotor, axial magnetic bearing, radial magnetic bearing, the displacement sensor, and the protecting bearing. The rotor is driven by the motor, which will provide constant angle momentum. The rotor is suspended by the magnetic bearing, and it has five degrees of freedom. This kind of magnetic bearing is called five-degree-of-freedom magnetic bearing. The displacement of each degree of freedom is measured by the displacement sensor. If the rotor is not on the given zero position, the error signal is passed through the magnetic bearing controller. The amplifier will output the corresponding control current, which will drive the magnetic bearing to produce bearing force to return the rotor to the given position. The base is the mechanical structure that connects double-gimbal magnetically suspended CMG with spacecraft. The output torque of double-gimbal magnetically suspended CMG is imparted to spacecraft by the base.

Figure 1: Prototype of double-gimbal magnetically suspended CMG.

The magnetic bearing is an elasticity support with gap. When it is suspended normally, the magnetically suspended rotor has small rotation angle and translation motion relative to gyro room. Meanwhile, the angle motion of magnetically suspended high-speed rotor is interrelated with the rotation of inner gimbal and outer gimbal. The reference frames should be defined firstly.

At the initial time, it is defined that the gimbal angle, the gimbal rate, the base angle position, and the base angular velocity are all zero. For spacecraft dynamic modeling, the coordinate frames are defined, which are illustrated in Figure 2.

Figure 2: The illustration of coordinate frames.

is the inertial frame. The origin is earth center of mass. points to vernal equinox; points to the north pole.

is the orbit frame. The origin is the spacecraft center of mass. and represent the roll axis and the pitch axis, respectively. Both of them are located in the orbit plane. The roll axis points to the go-ahead direction. The pitch axis is vertical to orbit plane. represents the yaw axis, which points to earth’s core. This frame rotates with speed , and the rotation direction is contrary to the pitch axis.

is the spacecraft body frame. The origin is the spacecraft center of mass, , , and are fixed on the spacecraft, and they are the moments of inertial principal axis of the spacecraft. is the 3-1-2 Euler angle with respect to orbit frame.

is the zero position frame. At the zero position, it coincides with the outer gimbal frame. It is fixed with body frame.

is the outer gimbal frame. It is fixed with outer gimbal. The outer gimbal angle is . It is positive when the rotor rotation speed along direction is positive.

is the inner gimbal frame. It is fixed with inner gimbal. The inner gimbal angle is . It is positive when the rotor rotation speed along direction is positive.

is the stator frame. It is fixed with gyro room. This frame keeps still with respect to inner gimbal frame. and point to central directions of the radial magnetic bearing.

is the inner loop frame. It is fixed with the suppositional inner loop. This frame has three translation motions and two rotation degrees of freedom. is the center of rotor mass. and in Figure 3 are the suppositional inner loop and the dummy outer loop frames, respectively. The outer gimbal angle is the Euler angle of with respect to . The inter gimbal angle is the Euler angle of with respect to .

Figure 3: The coordinate transformation matrix from the stator frame to the rotor frame.

is the rotor frame. This frame is fixed with rotor and has one rotation degree of freedom along axis with respect to inner loop frame. The rotor rotation angle is . is the constant angular velocity, and it is expressed in rotor frame coordinate as .

is used to denote the coordinate transformation matrix from frame to frame . denotes the vector projected in coordinate . denotes the skew symmetric matrix of . On the condition of small angle transformation, then

3. The Dynamic Model of Double-Gimbal Magnetically Suspended CMG

3.1. The Dynamic Model of Magnetically Suspended Rotor

In order to obtain large rotor angular momentum, the rotor is designed as disc shape to improve the moment of inertia. The moment of inertia ratio between the pole and the equator is large. In order to perform the dynamic analysis and control system design, the magnetic bearing dynamic model of magnetically suspended CMG must be built. The main nomenclature used in this paper is listed in the Nomenclature section.

The magnetically suspended rotor has six degrees of freedom, including translation motion along , , and direction and rotation motion along , , and direction. The rotation along direction is controlled by the motor. The other five degrees of freedom are controlled by magnetic bearing, including the translation motion along , , and direction and the radial rotation along , direction. The exclusive force between the magnetically suspended rotor and the gyro room is the magnetic bearing force. In Figure 4, the radial bearing forces and can be synthesized into radial force . The radial bearing forces and can be synthesized into radial force . The radial torque in and direction is and , which will influence the radial translation and the rotation of magnetically suspended rotor. When the magnetic bearing rotates about the centroid, the translation and rotation motion about the rotor centroid can be considered separately. Firstly, the translation and rotation motion equations of the rotor must be built. By synthesizing the equations, the dynamic equation of the magnetically suspended rotor can be obtained.

Figure 4: The sketch map of magnetic bearing dynamic model.

The relative motion of magnetically suspended rotor in noninertial frame suffers from two kinds of force. One is the active control force of magnetic bearing and the other is the implicated inertial force and Coriolis inertial force that result from the implicated motion of inner gimbal and outer gimbal. The rotor dynamic equation in noninertial frame is . is the mass of magnetically suspended rotor. is the relative translation acceleration of magnetically suspended rotor. is the active force, namely, the bearing force. and are the implicated acceleration and Coriolis acceleration. The relative translation acceleration is greater than the implicated acceleration and Coriolis acceleration, so the implicated inertial force and Coriolis inertial force can be ignored. The relative translation motion model of the magnetically suspended rotor is simplified as , where , .

Similar with the translation motion of magnetically suspended rotor, the rotation of magnetically suspended rotor in inertial frame suffers from two kinds of torque. One is the active control torque of magnetic bearing, and the other is implicated inertial torque and Coriolis inertial torque that are produced from the implicated motion of the noninertial frame. Considering the ground experimental condition, the absolute angle motion of magnetically suspended rotor includes the following three kinds of motion.

(1)The rotation of rotor frame relative to magnetic bearing stator is described with Euler angles  ,  , and , which rotates about the inner loop rotation axis, the outer loop rotation axis, and the rotor spin axis.(2)The rotation of the magnetic bearing stator frame relative to inner gimbal frame is denoted by .(3)The rotation of outer gimbal frame relative to inertial frame is denoted by .

In the following, based on Newton-Euler method, Euler dynamic equation is used to derive the dynamic model of magnetically suspended rotor of double-gimbal magnetically suspended CMG. The reference frame is the magnetic bearing stator frame , and the movable frame is the inner loop frame .

The angle speed vector of magnetically suspended rotor includes six parts, namely, the rotation speed of rotor frame relative to stator frame , the rotation speed of stator frame relative to inner gimbal frame , the rotation speed of inner gimbal frame relative to outer gimbal frame , the rotation speed of outer gimbal frame relative to zero position frame , the rotation speed of zero position frame relative to body frame , and the rotation speed of body frame relative to inertial frame : . Because the projection gap is small, the maximal radial rotation angle of the Euler angles  ,   is less than 3 × 10⁻³ rad. It can be considered that , , , .

(1) Derive the angular velocity projection of the movable frame about each axis. The angular velocity of rotor frame relative to stator frame is . Its projection in inner loop frame is

Because stator frame is still relative to inner gimbal frame, . The angular velocity of inner gimbal relative to outer gimbal is . Its projection in inner loop frame is

The angular velocity of outer gimbal relative to zero position frame is . Its projection in inner loop frame is

The zero position frame is still relative to body frame, . The rotation speed of spacecraft relative to orbit frame is , which corresponds to roll, pitch, and yaw. The rotation speed of orbit frame is . So, the rotation speed of spacecraft can be described in frame as follows:

The spacecraft rotation speed in inner loop frame is

The absolute angular velocity of magnetically suspended rotor in inner loop frame is denoted by

Then

(2) The moment of inertia of magnetic bearing in is . Because the rotor is symmetric, its equator moment of inertia is , and its pole moment of inertia is . So the rotor moment of inertia relative to inner loop frame has the following relationship with :

The rotor angular momentum in inner loop frame is

Then

The projection of the angular velocity in inner loop frame can be simplified as

(3) In this paper, five-degree-of-freedom magnetic bearing is used. The axial magnetic bearing only provides translation force, and it cannot output torque to the rotor. Therefore, the total torque of the rotor is

(4) In magnetic bearing frame , the dynamic equation of the double-gimbal magnetically suspended CMG is

where , , and are the elements of . According to (2)~(14), the rotor dynamic equation is

Because the projection gap is small and the permitted maximal radial rotation angle is less than 3′′, the rotor motion is restricted in the projection gap. and are small, , , , so the above equation can be simplified as

When inner gimbal and outer gimbal are still on the ground, the dynamic equation of magnetically suspended rotor can be simplified as

3.2. The Dynamic Model of Inner Gimbal

The absolute angular velocity of inner gimbal is .

(1) Taking inner gimbal frame as moveable frame, the angular velocity of outer gimbal relative to zero position is . Its projection in inner gimbal frame is

The zero position frame is still relative to body frame, so . The spacecraft rotation speed in inner gimbal frame is denoted by

So the absolute velocity and the time derivative of inner gimbal can be described in inner gimbal frame as

(2) The torque that acts on inner gimbal includes two parts, namely, the magnetic bearing torque and the outer gimbal torque. Then

The projection of this torque in inner gimbal frame is

The torque acting on the inner gimbal by the outer gimbal is

where and are the constraint torque along and direction. is along direction. So the total torque of inner gimbal resulting from the magnetic bearing and outer gimbal is

(3) double-gimbal magnetically suspended CMGs are installed on spacecraft according to certain configurations, which can avoid singularity state. In inner gimbal frame , the dynamic model of inner gimbal of the () double-gimbal magnetically suspended CMG is

where the angular momentum of inner gimbal of the double-gimbal magnetically suspended CMG is . Then

The moment of inertia of inner gimbal is ; then

Substituting (20), (21), and (25) into the above equation, the dynamic model of inner gimbal is

3.3. The Dynamic Model of Outer Gimbal

Supposing that is the angular velocity of outer gimbal relative to zero position frame, is the angular velocity of zero position frame relative to spacecraft body frame, and is the implicated angular velocity of spacecraft body frame relative to inertial frame. The absolute angular velocity of outer gimbal is .

(1) Taking outer gimbal frame as movable frame, zero position frame is kept stable relative to body frame: .

The spacecraft rotating velocity in inner gimbal frame is described as

So the absolute angular velocity of outer gimbal and its time derivative can be described as

(2) Similar with inner gimbal, the torque acting on outer gimbal is the torque from spacecraft and the torque from inner gimbal. The projection of inner gimbal torque on outer gimbal is

The torque that the spacecraft is acting on outer gimbal is

where is the outer gimbal motor driving torque along direction. and are the constraint torque in which the base acts on outer gimbal along and direction. So the total torque acting on outer gimbal is

(3) In outer gimbal frame , the outer gimbal dynamic model of the double-gimbal magnetically suspended CMG is

where the angular momentum of outer gimbal of the double-gimbal magnetically suspended CMG is . Then .

The moment of inertia of outer gimbal is ; then

By substituting (31), (32), and (34) into the above equation, the dynamic equation of outer gimbal is

4. Agile Maneuver Spacecraft Dynamic Modelling Using Double-Gimbal Magnetically Suspended CMG as the Actuator

The torque in spacecraft body frame can be described as

Supposing that the moment of inertia of the spacecraft is , then

where is the vector form of the total torque that acts on spacecraft. is the disturbance, which can be described as

where is the constant rotor angular momentum. The outer disturbance torque can be denoted by

Because the rotor displacement relative to the gimbal is small, , . Equation (41) can be described as

where the total angular momentum of double-gimbal magnetically suspended CMGs is . Then

where is the inner gimbal angle and the outer gimbal angle of double-gimbal magnetically suspended CMGs. and are the gimbal rate of the double-gimbal magnetically suspended CMG. is the Jacobin matrix. The outer disturbance torque can be described as

Equations (16), (29), (38), and (40) form the dynamic model of double-gimbal magnetically suspended CMG when the spacecraft maneuvers. It can be seen that the relative motion of any component will result in the relative motion of other components. The motion of all components is coupling, which makes the dynamic model of spacecraft with double-gimbal magnetically suspended CMG complex.

5. The Coupling Characteristic in Double-Gimbal Magnetically Suspended CMG

Equations (16), (29), and (38) form the dynamic model of double-gimbal magnetically suspended CMG when spacecraft maneuvers. By writing them together, then

From the above equations, the following can be seen.

(1) Equation (46) is the nonlinear dynamic model of the relative motion between magnetically suspended rotor, inner gimbal, and outer gimbal of double-gimbal magnetically suspended CMG when spacecraft maneuvers. These models are built in inner loop frame, inner gimbal frame, and outer gimbal frame, respectively. The relative motion of the three components must be considered for analyzing the spacecraft stability and precision.

(2) The magnetic bearing suspending with gap is used in double-gimbal magnetically suspended CMG. The motion of magnetically suspended rotor does not only rely on magnetic bearing force but also suffer from the influence of gimbal servo system and spacecraft motion.

(3) Between the inner gimbal and the outer gimbal, there exists the dynamic coupling that results from the gyro effect. This dynamic coupling not only depends on the gimbal motion but also is related with the radial rotation. When spacecraft maneuvers, the gyro coupling torque of double-gimbal magnetically suspended CMG is proportional to spacecraft angular velocity and the cosine of the radial angle displacement. So this system is a nonlinear system. The disturbance influence of the strong coupling torque on the magnetically suspended rotor can make the magnetically suspended rotor unstable.

In (46), and can be obtained by using the dynamic equation of magnetically suspended rotor, namely, (16). , , and can be obtained by substituting them into the dynamic model of inner gimbal, namely, (29). And then, they are substituted into the dynamic equation of outer gimbal, namely, (38), to obtain , , and . By ignoring the high-order item, the dynamic model of magnetic bearing along -direction, -direction of the inner and outer gimbal can be obtained:

Based on the dynamic model of double-gimbal magnetically suspended CMG in (46), the coupling torque along -direction of magnetic bearing, -direction of magnetic bearing, the inner gimbal, and the outer gimbal can be obtained. The coupling torque can be divided into two portions, including the gyro coupling torque and the inertial coupling torque, where the gyro coupling torque is proportional to the spacecraft angular velocity. The gyro coupling torque is

where , , , and are the gyro coupling torque along -direction of magnetic bearing, -direction of magnetic bearing, the inner gimbal, and the outer gimbal, respectively. The inertial coupling torque is proportional to the angle acceleration of the spacecraft and gimbal:

where , , , and are the inertial coupling torque along -direction of magnetic bearing, -direction of magnetic bearing, the inner gimbal, and the outer gimbal, respectively.

6. Conclusion

Double-gimbal magnetically suspended CMG is a novel attitude control actuator for the agile maneuver spacecraft. The rotor has five degrees of freedom besides the rotation degree, including three translation motions and two radial motions. So the dynamic model of the double-gimbal magnetically suspended CMG is complex. Moreover, the maneuver of spacecraft will influence the stability and precision of high speed rotor.

In order to verify the interrelation in double-gimbal magnetically suspended CMG, its working principal and basic structure are given firstly. Then, the dynamic model of the magnetically suspended rotor, inter gimbal, outer gimbal, and spacecraft using double-gimbal magnetically suspended control moment gyro as actuator is built. From these models, it can be seen that the model of one component includes the coupling moment produced by other components, and the relative motion of one component will affect the motion of other components. The strong coupling between the components is investigated, which will provide the foundation for the further research.

Nomenclature

: | Radial force in direction
: | Radial force in direction
: | The radial torque in direction
: | The radial torque in direction
: | The moments of inertia of spacecraft
: | The moments of inertia of outer gimbal
: | The moments of inertia of inner gimbal
: | The moments of inertia of rotor
: | The rotation of rotor frame relative to magnetic bearing stator rotating about the inner loop rotation axis
: | The rotation of rotor frame relative to magnetic bearing stator rotating about outer loop rotation axis
: | The rotation of rotor frame relative to magnetic bearing stator rotating about the rotor spin axis
: | The rotation of the magnetic bearing stator frame relative to inner gimbal frame
: | The rotation of outer gimbal frame relative to inertial frame
: | The angle speed vector of magnetically suspended rotor
: | The rotation speed of rotor frame relative to stator frame
: | The rotation speed of stator frame relative to inner gimbal frame
: | The rotation speed of inner gimbal frame relative to outer gimbal frame
: | The rotation speed of outer gimbal frame relative to zero position frame
: | The rotation speed of zero position frame relative to body frame
: | The rotation speed of body frame relative to inertial frame.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This research has been supported by National Natural Science Foundation of China under Grant 61121003, the National Basic Research Program (973 Program) of China under Grant 2009CB72400101C, and National Civil Aerospace Pre-Research Project.

References

J. C. Fang, X. B. Xu, J. Q. Tang, and H. Liu, “Adaptive complete suppression of imbalance vibration in AMB systems using gain phase modifier,” Journal of Sound and Vibration, vol. 332, no. 24, pp. 6203–6215, 2013.

J. C. Fang, X. B. Xu, and J. J. Xie, “Active vibration control of rotor imbalance in active magnetic bearing systems,” Journal of Vibration and Control, vol. 21, no. 4, pp. 684–700, 2015.

R. J. Bauer, “Kinematics and dynamics of a double-gimbaled control moment gyroscope,” Mechanism and Machine Theory, vol. 37, no. 12, pp. 1513–1529, 2002.

B. Wie, “Singularity escape/avoidance steering logic for control moment gyro systems,” Journal of Guidance, Control, and Dynamics, vol. 28, no. 5, pp. 948–956, 2005.

J. C. Fang and Y. Ren, “Decoupling control of magnetically suspended rotor system in control moment gyros based on an inverse system method,” IEEE/ASME Transactions on Mechatronics, vol. 17, no. 6, pp. 1133–1144, 2012.

J. Fang and Y. Ren, “Self-adaptive phase-lead compensation based on unsymmetrical current sampling resistance network for magnetic bearing switching power amplifiers,” IEEE Transactions on Industrial Electronics, vol. 59, no. 2, pp. 1218–1227, 2012.

J. Fang, S. Zheng, and B. Han, “AMB vibration control for structural resonance of double-gimbal control moment gyro with high-speed magnetically suspended rotor,” IEEE/ASME Transactions on Mechatronics, vol. 18, no. 1, pp. 32–43, 2013.

C. J. Heiberg, “A practical approach to modeling Single-Gimbal control momentum gyroscopes in agile spacecraft,” in Proceedings of the Guidance, Navigation, and Control Conference and Exhibit, Denver, Colo, USA, August 2000.

J. Ahmed and D. S. Bernstein, “Adaptive control of double-gimbal control-moment gyro with unbalanced rotor,” Journal of Guidance, Control, and Dynamics, vol. 25, no. 1, pp. 105–115, 2002.

L. Tang and Y. Chen, “Modeling and dynamic analysis of double-gimbal control moment gyros,” Acta Aeronautica et Astronautica Sinica, vol. 29, no. 2, pp. 424–429, 2008.

H.-D. Duan and Z. Wu, “Dynamics and simulation of spacecraft with magnetically suspended control moment gyroscope,” Journal of System Simulation, vol. 21, no. 9, pp. 2783–2790, 2009 (Chinese).

Y. Ren and J. Fang, “High-stability and fast-response twisting motion control for the magnetically suspended rotor system in a control moment gyro,” IEEE/ASME Transactions on Mechatronics, vol. 18, no. 5, pp. 1625–1634, 2013.

M. Xiang and T. Wei, “Autobalancing of high-speed rotors suspended by magnetic bearings using LMS adaptive feedforward compensation,” Journal of Vibration and Control, vol. 20, no. 9, pp. 1428–1436, 2014.

H.-T. Li, J.-C. Fang, B.-C. Han, and T. Wei, “Study on the system disturbance rejection method used in the gimbal servo system of double gimbal magnetically suspended control moment gyro,” Journal of Astronautics, vol. 30, no. 6, pp. 2199–2205, 2009 (Chinese).