Reaction Wheel Installation Deviation Compensation for Overactuated Spacecraft with Finite-Time Attitude Control

A novel attitude tracking control scheme is presented for overactuated spacecraft to address the attitude stabilization problem in presence of reaction wheel installation deviation, external disturbance and uncertain mass of moment inertia. An adaptive sliding mode control technique is proposed to track the uncertainty. A Lyapunov-based analysis shows that the compensation control law can guarantee that the desired attitude trajectories are followed in finite-time.The key feature of the proposed control strategy is that it globally asymptotically stabilizes the system, even in the presence of reaction wheel installation deviation, external disturbances, and uncertain mass of moment inertia. The attitude track performance using the proposed finite-time compensation control is evaluated through a numerical example.


Introduction
In present, nearly all of the highly accurate slewing maneuvers necessitate the use of nonlinear differential equations for the kinematics and dynamics during the control system design [1].However, the attitude tracking problem is further complicated by the external disturbance and uncertain mass of moment inertia.To address these issues, there have been several important developments in the design of feedback control laws for spacecraft maneuvering.A number of control design approaches using adaptive control [2,3], sliding mode control [4][5][6][7],  ∞ [8,9], optimal control [10][11][12][13], and data driven control [14][15][16] have been proposed.However, few of them focus on the reaction wheel installation deviation that are of great theoretical and practical interest.In fact, the installation deviation is a widespread phenomenon, such as the actuator misalignment which is limited by the installation technique or generated by materials deforms the vehicle violent vibration during the launching process.In the area of actuator misalignment compensation, there currently exist few unified frameworks for the design of simple control structures.
Several solutions to actuator installation deviation have been presented in the literature [17][18][19][20].In [17], the authors presented a general adaptive tracking attitude controller design framework for a spacecraft subject to the actuator installation minor angle deviation.In [18], an adaptive attitude tracking method is proposed to compensate the actuator misalignment of nearly 15 degree.And in [19], a novel algorithm is employed precisely to estimate the information, such as installation angle of wheel and CMG alignments.And then the controller design can be on for the estimation information.Moreover, another recent paper in [20] proposed an adaptive control approach for satellite formation flying, in which backstepping technique is used to synthesize a controller to handle thrust magnitude error and misalignment.However, the torque is different between thruster and reaction wheel, one is literal and the other is time-variable.That is to say that this control strategy is not suitable for reaction wheel installation deviation compensation for overactuated spacecraft attitude control.
Treating the uncertain mass of moment inertia caused by the reaction wheel misalignment is another impossibly avoided problem.In practice, in order to ensure the reliability of on-orbit spacecraft operation, especially under high altitude sever external environment, overactuatation is widely employed to guarantee the control system reliability service.And finite-time is meanwhile necessary for time critical missions.As a result, more and more investigations also have focused on attitude control design with finite-time convergence.In [21][22][23], the finite-time control technique was applied to design an attitude controller.Feng et al. [24] proposed a terminal sliding mode controller to solve the singular problem for a second-order nonlinear dynamic system.A terminal sliding mode and the Chebyshev neural network were used in [25] to guarantee that the attitude manoeuvre was accomplished in finite time, even in the presence of an unknown inertia matrix, external disturbances, and control input constraints.Furthermore, two robust sliding mode controllers were proposed in [26] to realize attitude tracking in finite-time.Similar finite-time fault tolerant controllers for spacecrafts were investigated in [27][28][29].This work focuses on developing a control scheme to perform attitude compensation for an overactuated spacecraft with reaction wheel installation deviation, external disturbances, and uncertain inertia parameters.More specifically, the attitude tracking error is required to be zero in finite time.The proposed approach is illustrated in Figure 1.The compensation control module is added to the output of the nominal controller to compensate for the reaction wheel misalignment, disturbances, and uncertain moment of inertia.The proposed scheme solves a difficult problem of reliable and high accuracy attitude tracking control in finite time that rejects external disturbances and, at the same time, compensates for actuator misalignment and system uncertainties so that the control objective is met.
The remainder of this paper is organized as follows.In Section 2, we summarize the mathematical model for the rigid spacecraft attitude and control problem.A compensation control solution with the misalignment, disturbance, and mass moment of inertia is presented in Section 3. Simulation results are presented in Section 4. Some conclusions are given in Section 5.

Mathematical Model and Problem Formulation
The notation adopted throughout this paper is introduced as follows.The symbol ‖⋅‖ denotes the standard Euclidean norm or its induced norm; the symbol ‖ ⋅ ‖ ∞ denotes the infinite norm of a vector or matrix.For any given matrix A ∈ R × with full row rank, A † denotes its pseudoinverse.

Dynamic Model of Rigid
Spacecraft.Consider a rigid space system described by the following attitude kinematics and dynamics equations [30]: where  ∈ R 3 is the angular velocity of a body-fixed reference frame expressed in the body-fixed reference frame, J ∈ R 3×3 (positive and definite) is the total inertia matrix of the spacecraft, u = [ 1  2  3 ]  ∈ R 3 denotes the combined control torque produced by the actuators, and ∈ R 3 denotes the external disturbance torque from the environment, which is assumed to be unknown but bounded;  0 , q are the scalar and vector components of the unit quaternion, respectively, with q = [ 1  2  3 ]  ∈ R 3 , satisfying the constraint  2 0 + q  q = 1; I 3 represents the identity matrix with proper dimensions, and for ∀a = , a × denotes a skew-symmetric matrix, more precisely,

Reaction Wheel Configuration with Installation Deviation.
For orbiting spacecraft, loosely speaking, they have more than three reaction wheels aligned with the spacecraft body axes.However, in practice, the configuration of actuators will never be perfect; that is, to say, whether due to finite manufacturing tolerances or warping of the spacecraft structure during launch, some alignment errors can always exist.Thus, in this section, the reaction wheels misalignment is taken into consideration; the faulty dynamics can be described by where D ∈ R 3× denotes the actuator distribution matrix, ΔD denotes the actuator distribution matrix induced by misalignment, and  = [ 1 , . . .,   ]  ∈ R  denotes the actual output torque of the  reaction wheels.Due to the rotation of the payload or the existence of the flywheel installation deviation, the moment of inertia J is uncertain but positive definite symmetric matrices and record J = J 0 + ΔJ, where J 0 denotes the nominal rotational inertia and ΔJ denotes the uncertain rotational inertia.Here set 0 < ‖ΔJ‖ ≤ ‖J‖ ≤  max < ∞ and  max is a positive constant.

Attitude Tracking Model.
Assume that the desired attitude to be followed is described with a desired frame T with respect to I. It is specified by the desired unit quaternion The desired angular velocity is denoted by  d ∈ R 3 .Let the error quaternion Q e = ( 0 , e  ) ∈ R × R 3 denote the attitude between B and T, and let  e ∈ R 3 represent the corresponding error angular velocity.One has where R ∈ R 3×3 denote the corresponding rotation matrix that brings T onto B, and With ( 1)-( 5), the attitude tracking error dynamics is given by:

Control Objective.
The control objective of this work can be stated as considering the uncertain attitude tracking system ( 8)- (10) and design a control law to guarantee that the attitude tracking error converges to zero in finite-time, even in the presence of actuator misalignment, uncertain inertia matrix, and external disturbance d().
We present now the main results of this study.

Finite-Time Attitude Compensation Control
For the proposed control approach shown in Figure 1, the nominal control power and the compensation control effort are presented in this section.First, a finite-time sliding mode surface is proposed.Then, based on the finite-time sliding mode surface, a compensation controller is synthesized and added to the nominal controller to guarantee the global asymptotic stability of the resulting closed-loop attitude tracking system with finite-time convergence.
3.1.Finite-Time Sliding Mode Surface Design.We first introduce some lemmas which will be utilized in the subsequent control development and analysis.
Lemma 1 (see [31]).If  ∈ (0, 1), then the following inequality holds for any vector x = ( 1 ,  2 , . . .,   )  ∈ R  : Lemma 2 (see [32]).Suppose that V(x) is a  1 smooth positive-definite function such that where Then for any initial value (0) =  0 , it follows that V(x()) = 0 for all the time  ≥   1 , To this end, in this work, a sliding mode surface is introduced as Proof.From the sliding mode theory [33], it is known that once the state trajectories of the attitude tracking system reach the sliding surface, that is, s = 0, it follows that Consider a candidate Lyapunov function as Because the inequality which implies that V 1 = 0 if and only if e = 0. Thus,  1 is really a Lyapunov function such that the signal e will converge to zero and, accordingly,  0 tends to ±1 as  → ∞ by using the constraint in (18).Note that the equilibrium point ( 0 , e) = (−1, 0) is not a stable equilibrium point [34].Then, by Lemma 1, we obtain Because e = (−1, 0)  is not the stable equilibrium point, the signal e will converge to zero.Thus, lim  → ∞  0 () = 1 can be obtained from the constraint e  e +  2 0 = 1.There exists a finite time  ≥ 0 such that  0 () > 0 for  ≥ .Then, for  ≥ , one has Then, Using ( 21), ( 19) can be further bounded by ) Using 0.5 < ( + 1)/2 < 1 and Lemma 2, one has  1 () ≡ 0 for all  ≥  1 .According to definition of  1 () in (17),  0 () ≡ 1, e() ≡ 0, and  e () ≡ 0 for all  ≥  1 are concluded.Thereby, the proof is completed here.

Attitude Compensation Controller Design.
Considering the reaction wheel installation deviation and external disturbance, it is obtained from the sliding surface (14) that Because J is unknown but bounded, then J = 0 is established.Then, L can be represented to be bounded by ‖L‖ ≤  0 +  1 ‖‖ +  2 ‖‖ 2 , where   ,  = 1, 2, 3 are positive constants [35].
In order to facilitate analysis and proof, firstly, define  =  min (DD  ), 3‖ΔD‖ ∞ ‖D † ‖ ∞ =  < 1, and D † is the pseudoinverse of D. Now, we are ready to present the main result in Theorem 4.

Theorem 4.
Considering the uncertainty attitude tracking dynamics described by (5) with actuator misalignment ΔD and external disturbance torque d(), design an attitude compensation control law as where where  ∈ R + is control parameter and  1 is carefully chosen such that where π1 =  1 − π1 , k =   − k ,  = 3, 4, 5, and  1 ,  3 ,  4 ,  5 are the positive constants.
Remark 5.The controller (24) includes three parts:  adp () is used to compensate for system uncertainty caused by the external disturbance and moment inertia,  mis () is used to accommodate actuator misalignment, and  nom () is the nominal control.Theorem 6.Consider the attitude tracking system given by (5), (9), and (10).If the control scheme (24) is implemented, then the attitude tracking maneuver can be accomplished in a finite time   =  1 +  2 ; that is,  e () ≡ 0 and e() ≡ 0 are guaranteed for all the time  ≥   .
Proof.It is obtained from Theorem 4 that all the states of the attitude tracking system reach the sliding mode surface s() = 0 in finite-time  2 and maintain the motion state on the slide mode surface.Furthermore, from Theorem 3, it is obtained that once the system state reaches the slide mode surface (14) the system state can reach the equilibrium point ( 0 , e) = (1,0) in finite time  1 .Therefore, for any initial state Q(0) and (0), the desired attitude trajectory can be followed in a finite time   ; that is, e() ≡ 0,  0 () ≡ 1, and  e () ≡ 0 are achieved for all the time  ≥   .Thereby, the proof is completed here.

Reaction Wheel Configuration.
To demonstrate the effectiveness and performance of the proposed compensation control scheme, numerical simulations have been carried out using the rigid spacecraft system ( 3) and ( 6) in conjunction with the developed compensation control law (24).The spacecraft is activated by four reaction wheels with a limited control torque  max = 0.1 N⋅m.The configuration of those four actuators is shown in Figure 2.   = 35.26∘ and   = 45 ∘ are the nominal alignment angles,  = 1, 2, 3, 4. Δ  and Δ  are the misalignment angles.
With the configuration shown in Figure 2, the relation between the actual output torque of reaction wheel and the total torque acting on the spacecraft is to be calculated as Although the misalignment angles exist due to finitemanufacture technique and vehicle vibration, those angles Mathematical Problems in Engineering Δ  , Δ  ( = 1, 2, 3, 4) are small values.They can be approximated by cos Hence, (36) can be re written as where D and ΔD = (ΔD 1 , ΔD 2 , ΔD 3 , ΔD 4 ) are calculated as ) , ) .
We see in Figures 3-8 the controller managed to stabilize the origin equilibrium point in 30 seconds with great pointing accuracy.Indeed, since the knowledge of spacecraft inertia parameters was not required and an implicit integral item was incorporated in the control law design, external disturbance effect on the attitude control performance can be compensated efficiently, and also great robustness to system uncertainties, such as misalignment, can be guaranteed.
We can see in Figures 3 and 4 the time responses of angle velocity and attitude angle; the proposed control scheme surely realized the high precision stable control in the presence of external disturbance, uncertain moment of inertia, and reaction wheel misalignment, and the pointing accuracy is superior to 0.01 ∘ ; the attitude stable precision is superior to 0.001 ∘ /s.Meanwhile, from Figures 5 and 6   we can see that for compensating the misalignment and other uncertainties, the designed control command of control redundancy configuration for 4 reaction wheels  is allocated to the three-actual-output torque u, and then the purpose of compensation, attitude high precision control is realized.In addition, the finite-time control validity is shown in Figure 8.And from Figure 8, we can see that the spacecraft attitude control system status has realized the tracking control at   = 30.5.Thereinto, the spacecraft attitude has arrived at slide mode surface s at  1 = 25.9;afterwards, the statuses converge to equilibrium point at  2 = 5.4 under the normal control  nom ().The same validity of finite-time attitude compensation control strategy proposed in this paper can be further proved from the time response of the quaternion as shown in Figure 7.
From the above illustrated simulation results, it is shown that the proposed scheme can accomplish the attitude stabilization in finite-time in presence of time-varying external disturbances, uncertain inertial parameters, and even reaction wheel installation deviation.

Conclusions and Future Works
Considering the spacecraft issues about reaction misalignment, external disturbances, and parameters uncertainty, in this paper, a finite-time adaptive attitude compensation control has been proposed.A quantitative installation deviation angle analysis has been done and given out the value range of the reaction wheel misalignment angle.In the end the system stability and engineering practical value have been discussed from the perspective of theory and engineering.Numerical simulation of this novel control strategy was also presented to confirm the advantages and improvements over existing controllers.The case of actuator misalignment mentioned in Section 4 had only discussed for four reaction wheel configuration, but this compensation control scheme is suitable for more than that reaction wheel number.Moreover, the actuator faults have not been considered.The latter case should be as one of subjects for future research.Meanwhile, the method optimal control approach combined robust control [37,38] also can be applied in this field.

Figure 1 :
Figure 1: Structure of the attitude compensation controller.

Figure 2 :
Figure 2: Configuration of four reaction wheels.

Figure 8 :
Figure 8: Time response of sliding mode surface.