Advanced
Satellite Attitude Control with a Modified Iterative Learning Law for the Decrease in the Effectiveness of the Actuator
Satellite Attitude Control with a Modified Iterative Learning Law for the Decrease in the Effectiveness of the Actuator
International Journal of Aeronautical and Space Sciences. 2010. Jun, 11(2): 87-97
Copyright ©2010, The Korean Society for Aeronautical Space Science
This is an Open Access article distributed under the terms of the Creative CommonsAttribution Non-Commercial License(http://creativecommons.org/licenses/by-nc/3.0/)which permits unrestricted non-commercial use, distribution, andreproduction in any medium, provided the original work is properly cited.
  • Published : June 15, 2010
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
HoJin Lee
Youdan Kim
ydkim@snu.ac.kr
Hee Seob Kim

Abstract
A fault tolerant satellite attitude control scheme with a modified iterative learning law is proposed for dealing with actuator faults. The actuator fault is modeled to reflect the degradation of actuation effectiveness, and the solar arrayinduced disturbance is considered as an external disturbance. To estimate the magnitudes of the actuator fault and the external disturbance, a modified iterative learning law using only the information associated with the state error is applied. Stability analysis is performed to obtain the gain matrices of the modified iterative learning law using the Lyapunov theorem. The proposed fault tolerant control scheme is applied to the rest-to-rest maneuver of a large satellite system, and numerical simulations are performed to verify the performance of the proposed scheme.
Keywords
1. Introduction
Thousands of satellites are now in operation for various purposes such as communication, navigation, military service, weather forecasting, and terrestrial and astronomical observations. Since satellite launches usually entail much cost and time, the social and economic cost that arises from satellite failure is critical. For satellite attitude control, reaction wheels, thrusters, and control momentum gyros are widely used as actuators. Therefore, the attitude control performance can be affected by actuator failure. If actuator failure occurs, the mission of the satellite cannot be accomplished or will be limited. Recently, the demand for satellite fault-tolerant control systems is increasing; therefore, considerable research has been undertaken that utilizes various control methods such as adaptive control, neural networks, sliding mode control, and so on (Henry, 2008; Jiang et al., 2008; Talebi and Patel, 2006; Tehrani et al., 2005; Wu and Saif, 2005). Tehrani used an estimator, based on neural networks, for fault diagnosis of the reaction wheel in a satellite attitude control system (Tehrani et al., 2005). Wu and Saif applied a sliding mode controller, based on neural networks, for satellite fault diagnosis (Wu and Saif, 2005). Henry used H-infinity and H_filter based schemes for the fault diagnosis of microscope satellite thrusters (Henry, 2008).
Several estimation and control schemes have been applied to deal with the failures of satellite systems. Especially, the iterative learning law is a fault estimation scheme for detecting and estimating faults. The iterative learning law using previous information can be applied to estimate both constant faults and time-varying faults (Chen and Saif, 2001, 2007).
In this paper, a fault tolerant control scheme with a modified iterative learning law is proposed to deal with the decreased effectiveness of satellite actuators. The modified iterative learning law using only information that is associated with the state error is adopted to estimate the influence of the actuator fault and the external disturbance. Lyapunov stability analysis is performed to obtain a stable controller. The performance of the proposed satellite attitude fault tolerant control scheme is verified by numerical simulations.
This paper is organized as follows. The second section describes the satellite system and the rest-to-rest maneuver of the satellite. In the third section, a fault tolerant controller with a modified iterative learning law is designed, and stability analysis is performed. In the fourth section, numerical simulations are performed to verify the performance of the proposed fault tolerant satellite control scheme. Finally, conclusions are drawn in the fifth section.
2. Satellite Dynamics and Mission
- 2.1 Satellite Dynamics
The satellite system considered in this paper is specified as follows (Chobotov, 1991).
Lager Image
Lager Image
In the above, q 0 , q 1 , q 2 , q 3 are the quaternion variables, ω x , ω y , ω z are the angular rates, the matrix I s includes the elements of the moment of inertia, T x , T y , T z are control input elements, and w dx , w dy , w dz are the external disturbances.
The quaternion vector q that is related to the Euler angle is defined as follows.
Lager Image
In Eq. (2), ( I s ω) x , ( I s ω) y and ( I s ω) z are defined as follows.
Lager Image
The dynamic equations of Eqs. (1) and (2) can be rewritten as follows.
Lager Image
Lager Image
Here, x 1 =[ q 0 q 1 q 2 q 3 ] T , x 2 =[ω x ω y ω z ] T ,
Lager Image
Lager Image
Note that B = I 3×3 for a healthy condition of the actuator. When one of the actuators partially loses its actuation effectiveness, Eq. (6) can be written as follows:
Lager Image
where the actuation effectiveness matrix, B f = diag { B fi },(∀=x~z), is a matrix that represents the fault effect of the actuator. The degradation of the actuation effectiveness in a real situation cannot be clearly modeled. Therefore, the chattering effect with high-frequency variation due to the degradation of the actuation effectiveness also should be considered. In this paper, the effects of the actuator fault, B fi are modeled as follows.
Lager Image
In Eq. (8), a fi is the magnitude of the actuation effectiveness, and the second term is the high-frequency variation due to the fault with amplitude b fi and frequency ω fi . Usually, the magnitude of b fi is very small when compared to that of a fi . For simplicity, let us introduce a matrix, C f , as follows
Lager Image
- 2.2 Rest-to-Rest Maneuver
The satellite usually performs various maneuvers such as rest-to-rest and despin. In this paper, the rest-to-rest maneuver is considered. The angle of maneuver of the principal axis of the satellite can be represented as follows.
Lager Image
In Eq. (10), f s ( . ) is a smooth approximation of the signum function for the maneuver.
For the rest-to-rest maneuver, the following boundary conditions should be satisfied.
Lager Image
In Eq. (11), θ f ? θ 0 is the maneuver attitude angle, t f is the target maneuver time, and T max is the maximum available torque. In this paper, the following function, f s ( . ), is adopted for near-minimum time control (Junkins and Kim, 1993) to prevent unnecessary activation of the flexible modes of the solar panels.
Lager Image
In Eq. (12), and the torque-shaping parameter, Δt=α t , α, is in the range of 0 <; α≤0.25.
The target maneuver time, t f , for each axis can be obtained by using Eqs. (10-12) as follows.
Lager Image
Note that each axis's desired reference trajectory for the satellite rest-to-rest maneuver can be obtained by using Eqs. (10, 12, 13).
3. Fault Tolerant Control Scheme with a Modified Iterative Learning Law
In this section, a fault tolerant control scheme with a modified iterative learning law is designed to deal with the decrease in the effectiveness of the actuator and the external disturbance.
First, let us design the sliding mode controller to make the satellite attitude angles track the reference trajectory in Eq. (10) using the pseudo control input. Note that the state vector, x 2 , in Eq. (5) can be considered as the pseudo control input. This pseudo control input becomes the reference trajectory of the angular rate, x 2 .
Using the quaternion vector and the desired trajectory, the sliding surface is defined as:
Lager Image
where the reference trajectory,
Lager Image
, is constructed to perform the rest-to-rest maneuver. Since the condition, q 0 2 + q 1 2 + q 2 2 + q 3 2 =1, should be satisfied, q 0 can follow q 0d as the sliding surface s approaches zero asymptotically.
Differentiating Eq. (14) with respect to time yields:
Lager Image
where
Lager Image
To satisfy the reaching condition of the sliding surface, the following Lyapunov candidate function is considered.
Lager Image
Differentiating Eq. (16) with respect to time and substituting Eq. (15) into the resulting equation yields:
Lager Image
The pseudo control input vector for the attitude angle can be obtained by the sliding mode control method as:
Lager Image
where K is the control gain matrix, and δ is the boundary layer thickness vector.
The pseudo control input of Eq. (18) is used as the reference trajectory of the angular rate, x 2d .
Substituting Eq. (18) into Eq. (17) gives:
Lager Image
Finally, it is concluded that the stability of the reference trajectory tracking controller is guaranteed.
Now, let us design a fault tolerant control scheme to deal with the decrease in effectiveness of the actuator and the external disturbance. Using the reference angular rate trajectory, x 2d , which is obtained from Eq. (18), and the angular rate state vector, x 2 , the angular rate state error vector is defined as:
Lager Image
Differentiating Eq. (20) with respect to time yields:
Lager Image
where T df ( t )= I s -1[ w d ( t )- C f u ( t )] is the influence of the actuator fault and the external disturbance.
Let us select the fault tolerant control input using Eq. (21) as:
Lager Image
where (F ^ ) is the estimated fault signal that can be obtained using the modified iterative learning law to compensate the influence of the actuator fault and the external disturbance. Usually, faults that occur in the system are unknown. To deal with this problem, the modified iterative learning law is used.
The general iterative learning law is updated by both the previous information and the state estimation error (Chen and Saif, 2001, 2007). Note that the modified iterative learning law proposed in this paper is updated using only the information associated with the state error; therefore, the influence of the actuator fault and the external disturbance can be estimated as follows:
Lager Image
where the parameter, τ, is an updating interval.
The gain matrices, L 1 and L 2 in Eq. (23), are determined using Lyapunov stability analysis. Substituting Eq. (22) into Eq. (21) and using Eq. (23) yield:
Lager Image
Let us consider the following Lyapunov candidate function.
Lager Image
Differentiating Eq. (25) with respect to time and substituting Eq. (24) into the resulting equation yield:
Lager Image
The following inequalities can be obtained (Yan et al., 1998).
Lager Image
Lager Image
The following equation also can be obtained using Eq. (24).
Lager Image
Similarly, the following inequalities can be obtained:
Lager Image
Lager Image
and
Lager Image
As a result, the inequality, Eq. (29), can be obtained using Eqs. (30-32) as follows:
Lager Image
where
Lager Image
and
Lager Image
By assuming
Lager Image
and
Lager Image
, the following can be derived:
Lager Image
where D max =max( D x , D y , D z ) and T max =max( T x , T y , T z ).
By substituting Eqs. (27, 28, 33, 34) into Eq. (26), the following equation can be obtained.
Lager Image
To satisfy Lyapunov stability, the following relation should be satisfied for the gain matrices, L 1 and L2.
Lager Image
Lager Image
From the above equation, the following inequality is obtained for using the gain matrix, L2:
Lager Image
where
Lager Image
Therefore, the gain matrix, L2, should be selected to satisfy the following relation.
Lager Image
In summary, the gain matrices, L1 and L2, must be chosen using Eqs. (36, 39), and the fault estimate signal, (수식삽입), of Eq. (23) is updated using the gain matrices and the information associated with the state error.
Remark : The state error can be used for fault monitoring purposes. When there are no faults in the system, the state error should be zero or close to zero. On the other hand, an increase in the state error would point to the occurrence of a fault. After the fault occurs, the estimated fault signal would learn about the fault and the state error will again be driven to zero or a value close to zero. This means that the iterative learning law can learn and update by the previous and the present state error information.
4. Numerical Simulation
Numerical simulation has been performed to verify the proposed fault tolerant control scheme. The Hubble space telescope is considered as the satellite system (Thienel and Sanner, 2007; Wie et al., 1993). The inertia, I s , is given as follows.
Lager Image
With regard to the external disturbance, the disturbance that is induced by the solar array is considered as follows (Wie et al., 1993).
Lager Image
The attitude (quaternion) reference trajectory is for the rest-to-rest maneuver. The maximum value of the available torque, T max , of the actuator is 0.82 Nm, and the torqueshaping α is selected as 0.25 for a near-minimum smooth maneuver. The initial attitude quaternion is chosen as , the initial reference angular velocity, ω d (0) is taken as zero, and the initial angular velocity is chosen as (deg/s)(Thienel and Sanner, 2007). The initial and target attitude angles are chosen as 0° and 40°, respectively, about all axes. The target attitude quaternion is calculated as . The final target maneuver time, t f , which is calculated by Eq. (13), is about 1,886 seconds. This final time is taken as the maximum target maneuver time at each axis. The reference trajectory of the rest-to-rest maneuver with respect to the quaternion values is shown in Fig. 1 .
Lager Image
Reference trajectory of the rest-to-rest maneuver (quaternion).
- 4.1 Case 1 ? Sequential Occurrence of Actuator Faults
In this case, the actuator faults are considered as follows.
Lager Image
In Eq. (8), the magnitude of the actuation effectiveness with regard to the actuator fault is chosen as 0.7, a fi = while b fi and the frequency, ω fi are chosen as 0.02 a fi and 5 Hz, respectively. Figure 2 shows the considered change in the fault effectiveness.
The gain matrix is chosen as K =0.1 I 3×3 , and the updating interval, τ, is chosen as 1 second. The parameters of the modified iterative learning law are chosen as ρ=0.1, γ 1 2 =0.1 and γ 3 4 5 =10. The upper bound on the external disturbance, D max , is chosen as 0.4 Nm.
The simulation results are shown in Figs. 3 - 6 . It can be seen from Figs. 3 and 4 that the state values, after the occurrence of the actuator faults, follow the reference trajectories very well. Figure 5 shows that the actuator faults can be estimated by the modified iterative learning law. As shown in Fig. 6 , a high-frequency control torque input is used for the y-axis. The control torque input response of the y-axis, from 0 to 50 seconds, is shown in Fig. 6(b) ; it can be seen that a high-frequency control input is required to deal with the disturbance that is induced by the solar array. As seen in the simulation results, the fault tolerant control scheme with the modified iterative learning law deals with the actuator faults and the external disturbance.
Lager Image
Model of the fault effect (B) (Case 1).
Lager Image
Attitude quaternion time responses (Case 1).
Lager Image
Angular rate time responses (Case 1).
Lager Image
Fault effect and disturbance estimation (Case 1).
Lager Image
Control torque input time response (Case 1).
- 4.2 Case 2 - Simultaneous Occurrence of Actuator Faults
In the second case, numerical simulation is performed for actuator faults that occur at same time. The actuator faults are considered to occur at t = t f - 800 sec about all axes. All the design parameters are identical to those in Case 1. Figure 7 shows the considered change in the fault effectiveness. The simulation results are shown in Figs. 8 - 11 . It can also be seen that the proposed fault tolerant control scheme with the modified iterative learning law copes well with the simultaneous occurrence of actuator faults.
Lager Image
Model of the fault effect (B) (Case 2).
Lager Image
Attitude quaternion time responses (Case 2).
Lager Image
Angular rate time responses (Case 2).
Lager Image
Fault effect and disturbance estimation (Case 2)
Lager Image
Control torque input time response (Case 2).
5. Conclusions
A satellite attitude control scheme is proposed to deal with actuator faults and external disturbances. The reference attitude angle of the satellite is designed using a near minimum-time maneuver in the rest-to-rest maneuver. A fault tolerant control scheme for the satellite attitude system with a modified iterative learning law is proposed to deal with the degradation of the actuation effectiveness and the external disturbance due to the vibration of the solar array. The modified iterative learning law is considered to estimate the unknown influence of the actuator fault and the external disturbance. Note that only information that is related to the state error is used. The fault tolerant control scheme is applied to a large satellite system, and the performance of the proposed satellite attitude control scheme is verified by numerical simulation. The proposed algorithm can be applied in an attitude control system to improve the reliability of satellite systems.
Acknowledgements
This study has been supported by the Korea Aerospace Research Institute (KARI) under the KOMPSAT-3 Development Program that is funded by the Ministry of Education, Science, and Technology (MEST) of the Republic of Korea.
References
Chen W , Saif M 2001 An iterative learning observerbased approach to fault detection and accomodation in nonlinear systems Proceedings of the IEEE Conference on Decision and Control Orlando FL 4469 - 4474
Chen W , Saif M 2007 Observer-based fault diagnosis of satellite systems subject to time-varying thruster faults Journal of Dynamic Systems Measurement and Control Transactions of the ASME 129 352 - 356    DOI : 10.1115/1.2719773
Chobotov V. A 1991 Spacecraft Attitude Dynamics and Control Original ed Krieger Malabar FL
Henry D 2008 Fault diagnosis of microscope satellite thrusters using H-infinity/H_filters Journal of Guidance Control and Dynamics 31 699 - 711    DOI : 10.2514/1.31003
Jiang T , Khorasani K , Tafazoli S 2008 Parameter estimation-based fault detection isolation and recovery for nonlinear satellite models IEEE Transactions on Control Systems Technology 16 799 - 808    DOI : 10.1109/TCST.2007.906317
Junkins J. L , Kim Y 1993 Introduction to Dynamics and Control of Flexible Structures American Institute of Aeronautics and Astronautics Washington DC
Talebi H. A , Patel R. V 2006 An intelligent fault detection and recovery scheme for reaction wheel actuator of satellite attitude control systems An intelligent fault detection and recovery scheme for reaction wheel actuator of satellite attitude control systems Munich Germany 3282 - 3287    DOI : 10.1109/CACSD-CCA-ISIC.2006.4777164
Tehrani E. S , Khorasani K , Tafazoli S 2005 Dynamic neural network-based estimator for fault diagnosis in reaction wheel actuator of satellite attitude control system Proceedings of the International Joint Conference on Neural Networks Montreal QC 2347 - 2352    DOI : 10.1109/IJCNN.2005.1556268
Wie B , Liu Q , Bauer F 1993 Classical and robust H(infinity) control redesign for the hubble space telescope Journal of Guidance Control and Dynamics 16 1069 - 1077    DOI : 10.2514/1.20591
Yan X. G , Wang J. J , Lu X. Y , Zhang S. Y 1998 Decentralized output feedback robust stabilization for a class of nonlinear interconnected systems with similarity IEEE Transactions on Automatic Control 43 294 - 299    DOI : 10.2514/3.21129