Next Article in Journal
An Improved Arc Flow Model with Enhanced Bounds for Minimizing the Makespan in Identical Parallel Machine Scheduling
Next Article in Special Issue
Nonlinearly Parametrized Modeling and Adaptive Control for a Generic Hypersonic Vehicle
Previous Article in Journal
A New Theobromine-Based EGFRWT and EGFRT790M Inhibitor and Apoptosis Inducer: Design, Semi-Synthesis, Docking, DFT, MD Simulations, and In Vitro Studies
Previous Article in Special Issue
An Approximate Method of System Entropy in Discrete-Time Nonlinear Biological Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Observer-Based Predefined-Time Attitude Control for Spacecraft Subject to Loss of Actuator Effectiveness

College of Automation Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 211106, China
*
Author to whom correspondence should be addressed.
Processes 2022, 10(11), 2294; https://doi.org/10.3390/pr10112294
Submission received: 9 October 2022 / Revised: 27 October 2022 / Accepted: 1 November 2022 / Published: 4 November 2022
(This article belongs to the Special Issue Advances in Nonlinear and Stochastic System Control)

Abstract

:
The predefined-time tracking problem of spacecraft attitude systems with loss of actuator effectiveness and lumped disturbance including the external disturbance and system uncertainty is studied. In order to obtain the estimation of the actuator efficiency factor more quickly and accurately, a robust learning observer is designed. Based on the fault reconstruction information of the learning observer and the predefined-time stability lemma, a predefined-time tracking fault-tolerant control scheme is proposed for the faulty spacecraft attitude system. The stability of the learning observer and the whole control system is verified by the Lyapunov stability theory. Finally, the effectiveness and advantages of the proposed scheme are illustrated by simulation results including comparisons with existing works.

1. Introduction

The attitude control system is an important part of spacecraft control system. It affects the flight performance and attitude tracking property, and plays a crucial role in the safe flight mission of spacecraft [1]. In recent years, many scholars have deeply studied the design of controllers in the case of spacecraft external disturbances, system modeling uncertainties, and actuator faults [2].
The control methods of spacecraft mainly include adaptive control [3], sliding mode control [4], neural network control [5], model predictive control [6], fuzzy logic control [7] and dynamic inversion control [8], and so on. Sliding mode control is widely used in the field of aircraft control because of its strong robustness to external disturbances and uncertainties of the system. It should be pointed out that most of the current advanced controller designs can prove that the control system is asymptotically stable, that is, the derivative of Lyapunov function V ˙ ( x ) 0 , which theoretically indicates that the system can be stable in an infinite time. In the field of aircraft control, due to the requirements of spacecraft tasks, the system needs to respond quickly; as such, some scholars have proposed the finite time control method. Reference [9] designed a finite time tracking controller based on fast terminal sliding mode control technology and neural network for a spacecraft attitude control system with sensor failure and actuator saturation. In the stability proof of reference [9], the Lyapunov function satisfies V ˙ λ V α ( x ) , where 0 < α < 1 , the system will be stable within a finite time T r , i.e., T r V 1 α ( x 0 ) λ ( 1 α ) , x 0 is the initial value of the system state variable. For the rendezvous task of thrust vector spacecraft, a constrained optimal orbit finite time attitude controller design scheme is proposed in [10]. In reference [11], a finite time disturbance observer is designed for spacecraft attitude control system with actuator faults and mismatched disturbances, and the estimation is applied to the design scheme of a finite time sliding mode fault-tolerant controller. References [10,11] proved that the objective of stability is that the derivative of Lyapunov function satisfies V ˙ ( x ) ρ 1 V ( x ) ρ 2 V ρ ( x ) , where ρ 1 > 0 , ρ 2 > 0 , 0 < ρ < 1 , and it can prove that the state of the system converges in a finite time T f , and T f T 0 + 1 ρ 1 ( 1 ρ ) ρ 1 V 1 ρ ( x ( t 0 ) ) + ρ 2 ρ 2 . It can be seen from [9,10,11] that the proof of finite time convergence of the system is related to the initial value of the system. When the initial value of the system changes or uncertain, it is difficult to calculate the convergence time of the system, which limits the application of finite time control.
In order to solve the problem that the finite time stability is related to the initial value, a fixed time controller with globally fixed convergence time of the system has been proposed. Reference [12] designed a fixed time attitude controller using fixed time theory and improved inversion control strategy. A fixed time attitude tracking control scheme for rigid spacecraft is proposed in [13] based on the idea of adding an integrator and inversion control technology. The proof objective of [12,13] is to make Lyapunov function satisfy V ˙ ( x ( t ) ) k 1 V ρ 1 ( x ( t ) ) k 2 V ρ 2 ( x ( t ) ) , where 0 < k 1 , k 2 , ρ 1 < 1 , ρ 2 > 1 . If the above inequality is satisfied, it can be shown that the system is stable in a fixed time T, and the expression of T is T 1 k 1 ( 1 ρ 1 ) + 1 k 2 ( ρ 2 1 ) . It can be seen from the proof process of fixed time stability that the fixed time value of system convergence is related to the controller parameters of the system. In most cases, it is difficult to establish a direct relationship between stable fixed time and control parameters. Once the fixed time value of system stability is set, it is difficult to make the system stable through several parameter adjustments. In order to solve the above problems, Sanchez Torres et al. proposed a controller design scheme with a fixed time upper bound as an adjustable parameter, which is named as a predefined-time stable system [14]. The convergence time of the predefined-time stable system can be preset arbitrarily by designing the parameters of the controller, which is independent of the initial value of the system. At present, the control theory related to the predefined-time stability is mainly applied to fractional order systems [15], robotic control systems [16], space robots [17], surface ships [18], hypersonic vehicle attitude control systems [19], etc. To the best of our knowledge, the results about predefined-time fault-tolerant control of spacecraft attitude control systems with external disturbances, system uncertainties and actuator failures are limited, which remains challenging and motivates us to do this study.
Learning observer is an advanced fault diagnosis mechanism, which can estimate constant fault, time-varying fault, and periodic fault. Its basic principle is to use the fault information at the previous time and the state output error information of the system for iterative learning to obtain the estimated value of the fault information at the current time [20]. For the distributed system, reference [21] reconstructs the fault value of the system by using the fault information of the previous time and the residual information of the previous time. This kind of observer can be called iterative learning observer. References [22,23] use the fault information of the previous time and the estimated error value of the current time state to reconstruct the fault of the system. This kind of observer can be called recursive learning observer. At present, most learning observers are used to estimate additive faults, and there are few research results on multiplicative faults (such as actuator efficiency loss faults).
Based on the above analysis, this paper proposes a predefined-time control scheme based on robust learning observer for a spacecraft attitude control system with external disturbances, system uncertainty and actuator fault. Firstly, the mathematical model of spacecraft attitude control system is given and a learning observer is designed to estimate the efficiency loss fault of the actuator. By using the estimated value of the learning observer, a fault-tolerant controller with predefined-time is designed. Lyapunov function is used to prove that the attitude angle of the system can track the desired command signal within the predefined time. Finally, the numerical simulation proves the effectiveness and feasibility of the proposed control scheme. The contributions of this paper are summarized as follows:
  • A robust learning observer is proposed to reconstruct efficiency loss faults. Compared with the traditional iterative learning observer [20] and recursive learning observer [21,22], the learning estimation law proposed in this paper is composed of the fault estimation value and state estimation error of the previous time and the state estimation error of the current time, which can quickly and accurately reconstruct the effectiveness value of the actuators.
  • Using the reconstructed information of actuator efficiency factor, a fault-tolerant attitude controller is designed based on the lemma of predefined-time stability and sliding mode control theory, so that the attitude angle of the system can track the desired command signal within predefined time.
  • Compared with the finite-time and fixed-time attitude control systems, the convergence time of the predefined-time stable system can be preset arbitrarily by tuning a simple parameter, which is independent of the initial value of the system.
The paper is organized as follows. The spacecraft attitude model and control objective is formulated in Section 2. The design of learning observer is proposed in Section 3 while the predefined fault-tolerant controller is given in Section 4. The simulation results are shown in Section 5. Conclusions are made in the last section.

2. The Model of Spacecraft and Problem Formulation

In this section, the spacecraft kinematics described using modified Rodrigues parameters and dynamics of attitude loop can be described as
σ ˙ = G ( σ ) ω ( J + Δ J ) ω ˙ = ω × ( J + Δ J ) ω + τ + d e
where σ = [ σ 1 , σ 2 , σ 3 ] T 3 denotes the spacecraft attitude angle vector, ω = [ p , q , r ] T 3 denotes the attitude rate vector and is measurable; J = d i a g { J x , J y , J z } denotes the inertia matrix, Δ J denotes the inertia uncertainty; τ = [ τ 1 , τ 2 , τ 3 ] T denotes the control torque vector; d e denotes the external disturbance vector; G ( σ ) = 1 4 [ ( 1 σ T σ ) I 3 + 2 σ × + 2 σ σ T ] , and σ × and ω × denote cross product matrices which defined as follows
σ × = 0 σ 3 σ 2 σ 3 0 σ 1 σ 2 σ 1 0 , ω × = 0 r q r 0 p q p 0 .
The three-axis attitude control torque of spacecraft is realized by the actuator of the system. This paper considers the situation where the actuators experience a partial loss of actuator effectiveness. Then, the spacecraft attitude dynamics in Equation (1) can be modified as
J ω ˙ = ω × J ω + D E ( t ) u + d
where D 3 × m denotes the control distributed matrix, u = [ u 1 , u 2 , , u m ] T denotes the desired control torque allocating to the individual actuator, m indicates the number of actuators. E ( t ) = d i a g { e 1 , e 2 , , e m } m × m denotes the health condition of the actuators, which can be called the actuator effectiveness matrix, and e i [ 0 , 1 ] . When e i = 1 implies that the ith actuator is operating healthy without any fault. When e i = 0 implies completely failure of ith actuator. d = d e ω × Δ J ω Δ ω ˙ is can be regarded as the lumped disturbance vector.
Control objective. For the spacecraft attitude control system with external disturbance, system uncertainty and actuator efficiency loss fault, the efficiency factor in the system is reconstructed by designing a learning observer, and then a predefined-time fault-tolerant tracking controller is designed based on the predefined-time stability lemma and sliding mode theory, so that the attitude of the system can accurately track the command signal within the predefined time. The structure block diagram of the observer-based predefined-time attitude control for spacecraft is designed as shown in Figure 1.
The following lemma and assumptions are introduced, which will be used for deriving our main results in the sequel.
Lemma 1 
(Predefined-time stability [14]). For autonomous control systems x ˙ ( t ) = f ( x , t ) , if there is exists a Lyapunov function V ( x ) : n with V ( 0 ) = 0 , V ( x ) > 0 ( x 0 ) . Then, it can be concluded that the control system is globally predefined-time stable with respect to T in sense that x ( t ) = 0 for all t T and x ( t 0 ) 0 , provided that
V ˙ ( x ) 1 ρ T e x p [ V ρ ( x ) ] V 1 ρ ( x ) , x ( t ) 0
where 0 < ρ < 1 , T > 0 denotes predefined time value.
Assumption 1. 
The nonlinear function ω × J ω is uniformly bounded for ω 3 with a Lipschitz constant k > 0 , which could be formulated in the following
| | ω × J ω ω ^ × J ω ^ | | k | | ω ω ^ | | , ( ω , ω ^ 3 )
Assumption 2. 
The lumped disturbance vector d = [ d 1 , d 2 , d 3 ] T is bounded, but the upper bound is unknown. Namely, | d i | < d m ( i = 1 , 2 , 3 ) , where d m is the unknown positive scalar.
Assumption 3. 
Because of physical limitations on the actuator, the control action generated is limited by the saturation value, and for simplicity, we assume that all the actuator control input torques have the same constraint value u m , that is, | u i | u m for i = 1 , 2 , , m .

3. Design of Learning Observer

Inspired by the previous learning observer design [24], we designed an improved learning observer design scheme for the efficiency factor of the actuator. The actuator loss of effectiveness fault E ( t ) is a diagonal matrix, and E ( t ) u in Equation (2) can be written as
E ( t ) u = U e ( t )
where U = diag { u 1 , u 2 , , u m } and e ( t ) = [ e 1 , e 2 , , e m ] T .
The faulty dynamics equation of rigid spacecraft (2) can be transformed into
J ω ˙ ( t ) = ω × ( t ) J ω ( t ) + D U e ( t ) + d ( t )
An improved learning observer is designed as follows
J ω ^ ˙ ( t ) = ω ^ × ( t ) J ω ^ ( t ) + D U e ^ ( t ) + M ( ω ( t ) ω ^ ( t ) ) + n sgn ( ω ( t ) ω ^ ( t ) ) e ^ ( t ) = l e ^ ( t τ ) + H 1 ( ω ( t ) ω ^ ( t ) ) + H 2 ( ω ( t τ ) ω ^ ( t τ ) )
where ω ^ ( t ) 3 and e ^ ( t ) m denote the angular velocity estimation and loss effectiveness fault estimation, respectively. t τ denotes the last time value and τ is usually taken as the sampling time internal of the system. Signal e ^ ( t ) is also called learning observer input, which is updated or driven successively by the state estimation error at the past and current time, and previous learning observer input e ^ ( t τ ) . l is a positive constant and H 1 m × 3 and H 2 m × 3 are positive-definite gain matrices chosen by designer. n is a positive sign function gain. The parameter M 3 × 3 is a positive-definite matrix.
Remark 1. 
In order to ensure that our estimated fault value remains within an interval ( 0 , 1 ] , we improve our learning update law (6) as follows
e ^ ( t ) = Proj ( e m i n , 1 ] { l e ^ ( t τ ) + H 1 ( ω ( t ) ω ^ ( t ) ) + H 2 ( ω ( t τ ) ω ^ ( t τ ) ) }
where e m i n indicates the lower limit of the efficiency loss factor, which is usually set to a very small value, such as e m i n = 0.0001 . P r o j { · } stand for projection operator.
To evaluate the estimation performance of the robust learning observer, two new variables have been defined as follows ω ˜ ( t ) = ω ( t ) ω ^ ( t ) and e ˜ ( t ) = e ( t ) e ^ ( t ) . The observer error equation can be obtained by subtracting (7) from (6)
J ω ˜ ˙ ( t ) = ω × ( t ) J ω ( t ) + ω ^ × ( t ) J ω ^ ( t ) + D U e ˜ ( t ) M ω ˜ ( t ) + d ( t ) n sgn ( ω ˜ ( t ) )
In order to express the main results, the following assumptions and lemmas are required.
Assumption 4. 
Since the efficiency loss factor of the actuator satisfies the constraint 0 e i ( t ) 1 , the following inequality can be held as | | e ( t ) l e ( t τ ) | | m ( 1 + l ) v ¯ , where v ¯ is a positive constant.
Lemma 2. 
If learning update law e ^ ( t ) is defined in Equation (6), the following inequality holds:
e ˜ T ( t ) e ˜ ( t ) b 1 v T ( t ) v ( t ) + b 2 l 2 e ˜ ( t τ ) e ˜ ( t τ ) + b 3 ω ˜ T ( t ) H 1 T H 1 ω ˜ ( t ) + b 4 ω ˜ T ( t τ ) H 2 T H 2 ω ˜ ( t τ )
where v ( t ) = e ( t ) l e ( t τ ) and b i ( i = 1 , 2 , 3 ) is a number greater than 0, and its value will be given in the proof section.
Proof. 
The efficiency factor estimation error of the observer is
e ˜ ( t ) = e ( t ) l e ^ ( t τ ) H 1 ω ˜ ( t ) H 2 ω ˜ ( t τ ) = v ( t ) + l e ˜ ( t τ ) H 1 ω ˜ ( t ) H 2 ω ˜ ( t τ )
Next, we can get
e ˜ T ( t ) e ˜ ( t ) = v T ( t ) v ( t ) l 2 e ˜ T ( t τ ) ω ˜ ( t ) + ω ˜ T ( t ) H 1 T H 1 ω ˜ ( t ) + ω ˜ T ( t τ ) H 2 T H 2 ω ˜ ( t τ ) + 2 l e ˜ T ( t τ ) v ( t ) 2 ω ˜ T ( t ) H 1 T v ( t ) 2 ω ˜ T ( t τ ) H 2 T v ( t ) 2 l ω ˜ T ( t ) H 1 T e ˜ ( t τ ) 2 l e ˜ T ( t τ ) H 2 T ω ˜ ( t τ ) + 2 ω ˜ T ( t τ ) H 2 T H 1 ω ˜ ( t )
Using the well-known Young’s inequality, that is ± 2 x T y a i x T x + 1 a i y T y for all x , y 3 , we have, for arbitrary positive constants a i ( i = 1 , 2 , 3 , 4 , 5 , 6 ) ,
2 l e ˜ T ( t τ ) v ( t ) a 1 l 2 e ˜ T ( t τ ) e ˜ ( t τ ) + 1 a 1 v T ( t ) v ( t ) 2 ω ˜ T ( t ) H 1 T v ( t ) a 2 ω ˜ T ( t ) H 1 T H 1 ω ˜ ( t ) + 1 a 2 v T ( t ) v ( t ) 2 ω ˜ T ( t τ ) H 2 T v ( t ) a 3 ω ˜ T ( t τ ) H 2 T H 2 ω ˜ ( t τ ) + 1 a 3 v T ( t ) v ( t ) 2 l ω ˜ T ( t ) H 1 T ω ˜ ( t τ ) a 4 ω ˜ T ( t ) H 1 T H 1 ω ˜ ( t ) + l 2 a 4 e ˜ T ( t τ ) e ˜ ( t τ ) 2 l e ˜ T ( t τ ) H 2 ω ˜ ( t τ ) a 5 l 2 e ˜ T ( t τ ) e ˜ ( t τ ) + 1 a 5 ω ˜ T ( t τ ) H 2 T H 2 ω ˜ ( t τ ) 2 ω ˜ T ( t τ ) H 2 T H 1 ω ˜ ( t ) a 6 ω ˜ T ( t τ ) H 2 T H 2 ω ˜ ( t τ ) + 1 a 6 ω ˜ T ( t ) H 1 T H 1 ω ˜ ( t )
Combining above inequalities into (11) leads to
e ˜ T ( t ) e ( t ) ( 1 + 1 a 1 + 1 a 2 + 1 a 3 ) v T ( t ) v ( t ) + ( 1 + a 1 + 1 a 4 + a 5 ) l 2 e ˜ T ( t τ ) e ˜ ( t τ ) + ( 1 + a 2 + a 4 + 1 a 6 ) ω ˜ T ( t ) H 1 T H 1 ω ˜ ( t ) + ( 1 + a 3 + 1 a 5 + a 6 ) ω ˜ T ( t τ ) H 2 T H 2 ω ˜ ( t τ ) = b 1 v T ( t ) v ( t ) + b 2 l 2 e ˜ T ( t τ ) e ˜ ( t τ ) + b 3 ω ˜ T ( t ) H 1 T H 1 ω ˜ ( t ) + b 4 ω ˜ T ( t τ ) H 2 T H 2 ω ˜ ( t τ )
where b 1 = 1 + 1 / a 1 + 1 / a 2 + 1 / a 3 , b 2 = 1 + a 1 + 1 / a 4 + a 5 , b 3 = 1 + a 2 + a 4 + 1 / a 6 and b 4 = 1 + a 3 + 1 / a 5 + a 6 . This completes the proof. □
Theorem 1 (Main Result).
Consider estimation error Equation (9) satisfying Lemma 1. An improved learning observer is designed in the form of (7) with the observer gains selected as c 1 = λ m i n ( M ) a 7 2 k 1 b 3 | | H 1 | | 2 ( 1 + χ + u m a x 2 2 a 7 | | D | | 2 ) > 0 , c 2 = 1 b 2 l 2 ( 1 + χ + u m a x 2 2 a 7 | | D | | 2 ) 0 , c 3 = 1 b 4 | | H 2 | | 2 ( 1 + χ + u m a x 2 2 a 7 | | D | | 2 ) 0 and n d m , then efficiency factor estimation error e ˜ ( t ) and the state estimation error ω ˜ ( t ) will converge to a small interval containing the origin.
Proof. 
Select the Lyapunov function as follows
V o ( t ) = 1 2 ω ˜ T ( t ) J ω ˜ ( t ) + t τ t e ˜ ( s ) e ˜ ( s ) d s + t τ t ω ˜ ( s ) T ω ˜ ( s ) d s
Taking the time derivation of V o yield
V ˙ o ( t ) = ω ˜ T ( t ) J ω ˜ ˙ ( t ) + e ˜ T ( t ) e ˜ ( t ) e ˜ T ( t τ ) e ˜ ( t τ ) + ω ˜ T ( t ) ω ˜ ( t ) ω ˜ T ( t τ ) ω ˜ ( t τ )
Substituting (9) into V ˙ o , (16) can be further obtained as follows
V ˙ o ( t ) = ω ˜ T ( t ) ( ω ^ × ( t ) J ω ^ ( t ) ω × ( t ) J ω ( t ) ) + ω ˜ T ( t ) D U e ˜ ( t ) + ω ˜ T ( t ) [ d n s g n ( ω ˜ ( t ) ) ] ω ˜ T ( t ) M ω ˜ ( t ) + e ˜ T ( t ) e ˜ ( t ) e ˜ T ( t τ ) e ˜ ( t τ ) + ω ˜ T ( t ) ω ˜ ( t ) ω ˜ T ( t τ ) ω ˜ ( t τ )
Using the Young inequality principle, the following inequality holds
ω ˜ T ( t ) D U e ˜ ( t ) a 7 2 ω ˜ T ( t ) ω ˜ ( t ) + u m a x 2 2 a 7 | | D | | 2 e ˜ T ( t ) e ˜ ( t )
where a 7 is a positive constant. Substituting the inequality (18) into (17) leads to
V ˙ o ( t ) [ λ m i n ( M ) a 7 2 k 1 ] · | | ω ˜ ( t ) | | 2 χ | | e ˜ ( t ) | | 2 ( n d m ) | | ω ˜ ( t ) | | + ( 1 + u m a x 2 2 a 7 | | D | | 2 + χ ) · e ˜ ( T ) ( t ) e ˜ ( t ) | | e ˜ ( t τ ) | | 2 | | ω ˜ ( t τ ) | | 2
where χ is a positive constant.
Substituting Lemma 2 and Assumption 4 into Equation (18), one has
V ˙ o ( t ) [ λ m i n ( M ) a 7 2 k 1 b 3 | | H 1 | | 2 ( 1 + χ + u m a x 2 2 a 7 | | D | | 2 ) ] | | ω ˜ ( t ) | | 2 [ 1 b 2 l 2 ( 1 + χ + u m a x 2 2 a 7 | | D | | 2 ) ] | | e ˜ ( t τ ) | | 2 χ | | e ˜ ( t ) | | 2 [ 1 b 4 | | H 2 | | 2 ( 1 + χ + u m a x 2 2 a 7 | | D | | 2 ) ] | | ω ˜ ( t τ ) | | 2 ( n d m ) | | ω ˜ ( t ) | | + b 1 ( 1 + χ + u m a x 2 2 a 7 | | D | | 2 ) m ( 1 + l ) 2
Then, if the gain of the learning observer satisfies the following inequalities,
c 1 = λ m i n ( M ) a 7 2 k 1 b 3 | | H 1 | | 2 ( 1 + χ + u m a x 2 2 a 7 | | D | | 2 ) ] > 0 c 2 = 1 b 2 l 2 ( 1 + χ + u m a x 2 2 a 7 | | D | | 2 ) 0 c 3 = 1 b 4 | | H 2 | | 2 ( 1 + χ + u m a x 2 2 a 7 | | D | | 2 ) 0 c 4 = n d m 0
the derivation of V ( t ) respect to time yields to
V ˙ o ( t ) c 1 | | ω ˜ ( t ) | | 2 χ | | e ˜ ( t ) | | 2 + c 5
where c 5 = b 1 ( 1 + χ + u m a x 2 2 a 7 | | D | | 2 ) m ( 1 + l ) 2 . It is further obtained from inequality (21) that V ˙ o 0 when | | ω ˜ ( t ) | | > c 5 c 1 or | | e ˜ ( t ) | | > c 5 χ . According to theorem in [25], it can be concluded that the learning observer error system is uniformly ultimately bounded.
Correspondingly, it is concluded that the estimation errors of angular velocity and fault are ultimately uniformly bounded by
| | ω ˜ | | c 5 c 1 , | | e ˜ ( t ) | | c 5 χ
By adjusting the parameters to increase the value of c 1 and χ or decrease the value of c 5 , the estimation error of the observer can converge to arbitrarily small residual interval. Then, the proof of Theorem 1 is completed. □

4. Design of Predefined-Time Attitude Tracking Controller

In this section, a predefined-time attitude fault-tolerant tracking controller is designed by using the lemma of predefined-time stable and the reconstructed information E ^ ( t ) , which can be obtained by learning observer, where E ^ ( t ) = diag ( e ^ 1 , e ^ 2 , , e ^ m ) . E ˜ ( t ) m × m is defined as the estimation error matrix, which satisfies E ˜ ( t ) = E ( t ) E ^ ( t ) . From Theorem 1, e ˜ i converges to an arbitrarily small value, and hence the following inequality can be achieved | | D E ˜ ( t ) u + d | | | | D | | u m + d m d ¯ , where d ¯ is a unknown positive constant. The mathematical model expression of spacecraft attitude control system can be rewritten as
σ ˙ = G ( σ ) ω J ω ˙ = ω × J ω + D E ^ u + D E ˜ u + d
For the faulty attitude system (23), define two new error variables z 1 = σ σ d , z 2 = σ ˙ σ ˙ d , where σ d is the desired attitude command; therefore, the attitude tracking error dynamics can be expressed as
z ˙ 1 = z 2 z ˙ 2 = d G ( σ ) d t ω + G ( σ ) ω ˙ σ ¨ d
Define the Lyapunov function of attitude angle error V a = 1 2 z 1 T z 1 . Select a sliding mode variable as follows
s = z 2 + ξ
where ξ = 1 2 h 1 T 1 exp ( V a h 1 ) V a h 1 z 1 , 0 < h 1 < 1 , and T 1 is the predefined time value.
According to Lemma 1, the predefined-time attitude tracking controller can be expressed as
u = ( D E ^ ) { ω × J ω k s sgn ( s ) J G 1 ( d G ( σ ) d t ω σ ¨ d + ξ ˙ + 1 2 h 2 T 2 exp ( V b h 2 ) ( V b h 2 ) s ) }
where 0 < h 2 < 1 , T 2 is the predefined time value. E ˙ is the estimation of efficiency factor, which can be obtained from learning observer (7). ( D E ^ ) denotes the pseudo-inverse operation of the matrix D E ^ . k s is the parameter gain of the system and its value is taken as k s > d ¯ . V b is the Lyapunov function of sliding mode variable, V b = 1 2 s T s .
Theorem 2. 
For the faulty spacecraft attitude systems (23), by designing a fault-tolerant controller based on the estimation information from the learning observer (7), the control system can be stabilized in the predefined time T = T 1 + T 2 , and the attitude of the system can track the desired attitude angle command in the predefined time T.
Proof. 
Derivation of the Lyapunov functions V b on sliding surfaces variable, we can obtain
V ˙ b = s T s ˙ = s T [ d G ( σ ) d t ω σ ¨ d + ξ ˙ + G ( σ ) J 1 ( ω × J ω + D E ^ u + D E ˜ u + d ) ]
Substituting control law (26) into above equation gives
V ˙ b = s T [ G ( σ ) J 1 ( ( D E ˜ u + d ) k s sgn ( s ) ) 1 2 h 2 T 2 exp ( V b h 2 ) ( V b h 2 ) s ] 1 h 2 T 2 exp ( V b h 2 ) ( V b 1 h 2 )
According to Lemma 1, when t > T 2 , V b 0 and sliding mode variable s will converge to 0 in the predefined time T 2 .
When s = 0 , combining with (25), we have
z 2 = ξ = 1 2 h 1 T 1 exp ( V a h 1 ) V a h 1 z 1
Taking the derivative of V a with regard to time t, one has
V ˙ a = z 1 T z ˙ 1 = 1 2 h 1 T 1 z 1 T exp ( V a h 1 ) V a h 1 z 1 = 1 h 1 T 1 exp ( V a h 1 ) V a 1 h 1
According to Lemma 1, when s = 0 , tracking error vector z 1 will converge to 0 within the predefined time T 1 . In conclusion, the attitude system of the spacecraft will be stable within the predefined time T = T 1 + T 2 under the action of the controller, and the attitude angle of the system will track the reference command within the predefined time T. The proof of Theorem 2 is completed. □

5. Simulation

In order to verify the effectiveness and feasibility of the control scheme proposed in this paper, a learning observer (LO) and a predefined-time controller (PTC) are designed for a class of spacecraft operating in circular orbit. At the same time, they are compared with the existing adaptive sliding mode observer (ASMO) [26] and terminal sliding mode control (TSMC) schemes [10].
The initial value of attitude angle and angular rate of spacecraft attitude control system is set as σ ( 0 ) = [ 0.3 , 0.2 , 0.2 ] T and ω = [ 0 , 0 , 0 ] T . The learning time interval of the learning observer is τ = 0.01 . The external disturbance is set to d e = 0.01 ( | | ω | | 2 + 0.05 ) · [ sin ( 0.8 t ) , cos ( 0.5 t ) , cos ( 0.3 t ) ] T . The moment of inertia of the system is set to J = [ 36 , 0 , 0 ; 0 , 17 , 0 ; 0 , 0 , 26 ] T . In this section, we consider that the number of actuators is m = 3 , and they are installed orthogonally, that is, D is the identity matrix. The desired attitude angle command is σ d = [ 0 , 0 , 0 ] T . The design controller parameters and the observer values are given in Table 1.
In the simulation, the following two cases are considered.
Case 1. Proposed LO vs. adaptive sliding mode observer in [26].
In this case, the main purpose is demonstrating the validity of the proposed learning observer approach. In order to facilitate the analysis, We set different efficiency factors on the three actuators: constant or time-varying faults. Denote E ( t ) = diag { e 1 ( t ) , e 2 ( t ) , e 3 ( t ) } and the fault scenario is set as
e 1 ( t ) = 1 , t < 5 s 0.8 , t 5 s , e 2 ( t ) = 1 , t < 8 s 0.6 , t 8 s , e 3 ( t ) = 1 , t < 10 s 0.4 + 0.05 cos ( 0.25 t ) , t 10 s
Figure 2 shows the response curve of fault reconstruction by two observers. The two observers have good estimation performance for constant faults, but the adaptive sliding mode observer can not estimate time-varying faults well. Figure 3 shows the fault estimation error curves of the two observers. When no efficiency loss fault occurs in the system, the adaptive sliding mode observer is greatly affected by the comprehensive disturbance. The estimation error of the second channel e ˜ 2 fluctuates greatly. It takes 5 s for the adaptive sliding mode observer to make the estimation error converge to a small range. When the efficiency loss fault occurs, the learning observer can accurately estimate the efficiency factor in about 3 s. It can be seen from the Figure 4 that no matter what form of efficiency loss fault occurs, the learning observer can accurately reconstruct the fault factor within 3 s after the fault occurs. Compared with the learning observer, the adaptive sliding mode observer takes a longer time (about 8 s) to track the constant fault; however, for the time-varying fault, their estimation results are poor and there is a large estimation error. Figure 4 shows the estimation error curve of system state variables ω . Both observers have good estimation effect on the state variables of the system, and the estimation error converges to a very small interval. In conclusion, the learning observer proposed in this paper can quickly and accurately estimate the efficiency loss failure (constant or time-varying) of the actuator.
Case 2. Proposed LO-based predefined-time controller (PTC) vs. terminal sliding mode controller (TSMC) in [10].
For comparison, a terminal sliding mode fault tolerant controller is implemented in the attitude control system. Figure 5 demonstrates the attitude angle tracking performance and Figure 6 shows the response of system angular velocity. Figure 7 is the response of actuator torque under the two different control methods. It can be seen from Figure 5 that the system attitude angle changes gently under the predefined-time controller, and the attitude angle tracks the reference command within 20 s. Figure 6 and Figure 7 show that the output torque of the actuator directly affects the angular rate change of the system. When each actuator appears, the control torque change of the scheme proposed in this paper is relatively gentle, which will not cause a large range of changes in the angular rate curve. When no fault occurs (the first 5 s), the angular rate converges to the neighborhood of 0 in 3 s. At the time of failure, the attitude angular rate has only a small fluctuation, and then returns to the stable state within 2 s.
Overall, the results of the above two cases fully verify the stability and robustness of the learning observer-based predefined-time tracking control system under external disturbances, system uncertainties and loss of actuator effectiveness.

6. Conclusions

This paper proposes an observer-based fault-tolerant tracking control scheme for a spacecraft attitude system with loss of faults effectiveness. Firstly, a learning observer is designed to reconstruct the efficiency factor of the system actuator. It is theoretically proved that the estimation error of the observer can converge to any small interval. Based on the predefined-time control theory, a fault-tolerant tracking control scheme is designed using the fault reconstruction information of the learning observer. The numerical simulation is applied to compare with the existing observer and control method, which shows the advantages of the scheme proposed in this paper. This paper only considers the case that the actuator has efficiency loss failure, but we do not consider the case that the actuator has limited output and multiple faults. In addition, the control allocation problem of actuators will be the subject of future research.

Author Contributions

Conceptualization, T.C. and H.G.; methodology, T.C.; software, B.H.; validation, T.C., H.G. and B.H.; formal analysis, T.C. and B.H.; writing—original draft preparation, T.C.; writing—review and editing, T.C.; supervision, H.G.; project administration, H.G. and T.C.; funding acquisition, H.G. and T.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported in part by Postgraduate Research & Practice Innovation Program of Jiangsu Province (No. KYCX18_0299). The authors also gratefully acknowledge financial support from the program of China Scholarships Council (No. 201806830102).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank the anonymous reviewers for their constructive suggestions to improve the quality of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lee, D. Fault-tolerant finite-time controller for attitude tracking of rigid spacecraft using intermediate quaternion. IEEE Trans. Aerosp. Electron. Syst. 2021, 57, 540–553. [Google Scholar] [CrossRef]
  2. Hasan, M.N.; Haris, H.; Qin, S.Y. Vibration suppression and fault-tolerant attitude control for flexible spacecraft with actuator faults and malalignments. Aerosp. Sci. Technol. 2022, 120, 107290. [Google Scholar] [CrossRef]
  3. Xie, S.; Chen, Q. Adaptive nonsingular predefined time control for attitude stabilization of rigid spacecrafts. IEEE Trans. Circuits Syst. II Express Briefs 2022, 69, 189–193. [Google Scholar] [CrossRef]
  4. Tang, P.; Zhang, F.; Ye, J.; Lin, D. An integral TSMC-based adaptive fault-tolerant control for quadrotor with external disturbances and parametric uncertainties. Aerosp. Sci. Technol. 2021, 109, 106415. [Google Scholar] [CrossRef]
  5. Chen, Q.; Xie, S.; He, S. Neural-network-based adaptive singularity-free fixed-time attitude tracking control for spacecrafts. IEEE Trans. Cybern. 2021, 51, 5032–5045. [Google Scholar] [CrossRef]
  6. Chai, R.; Tsourdos, A.; Gao, H.; Xia, Y.; Chai, S. Dual-loop tube-based robust model predictive attitude tracking control for spacecraft with system constraints and additive disturbances. IEEE Trans. Ind. Electron. 2022, 69, 4022–4033. [Google Scholar] [CrossRef]
  7. Zhang, Z.; Shi, Y.; Yan, W. A novel attitude-tracking control for spacecraft networks with input delays. IEEE Trans. Control. Syst. Technol. 2020, 29, 1035–1047. [Google Scholar] [CrossRef]
  8. Wang, F.; Hou, M.; Cao, X.; Duan, G. Event-triggered backstepping control for attitude stabilization of spacecraft. J. Frankl. Inst. 2019, 356, 9474–9501. [Google Scholar] [CrossRef]
  9. Ye, D.; Xiao, Y.; Sun, Z.; Xiao, B. Neural Network Based Finite-Time Attitude Tracking Control of Spacecraft With Angular Velocity Sensor Failures and Actuator Saturation. IEEE Trans. Ind. Electron. 2022, 69, 4129–4136. [Google Scholar] [CrossRef]
  10. Zhang, J.; Biggs, J.D.; Ye, D.; Sun, Z. Finite-time attitude set-point tracking for thrust-vectoring spacecraft rendezvous. Aerosp. Sci. Technol. 2020, 96, 105588. [Google Scholar] [CrossRef]
  11. Guo, B.; Chen, Y. Adaptive fast sliding mode fault tolerant control integrated with disturbance observer for spacecraft attitude stabilization system. ISA Trans. 2019, 94, 1–9. [Google Scholar] [CrossRef] [PubMed]
  12. Du, H.; Zhang, J.; Wu, D.; Zhu, W.; Li, H.; Chu, Z. Fixed-time attitude stabilization for a rigid spacecraft. ISA Trans. 2020, 98, 263–270. [Google Scholar] [CrossRef] [PubMed]
  13. Zou, A.; Kumar, K.D.; de Ruiter, A.H.J. Fixed-time attitude tracking control for rigid spacecraft. Automatica 2020, 113, 108792. [Google Scholar] [CrossRef]
  14. Sanchez-Torres, J.D.; Gomez-Gutierrez, D.; López, E.; Loukianov, A.G. A class of predefined-time stable dynamical systems. IMA J. Math. Control. Inf. 2018, 35, 1–29. [Google Scholar] [CrossRef] [Green Version]
  15. Munoz-Vazquez, A.J.; Sanchez-Torres, J.D.; Defoort, M.; Boulaaras, S. Predefined-time convergence in fractional-order systems. Chaos Solitons Fractals 2021, 143, 110571. [Google Scholar] [CrossRef]
  16. Munoz-Vazquez, A.J.; Sanchez-Torres, J.D.; Jimenez-Rodriguez, E.; Loukianov, A.G. Predefined-Time Robust Stabilization of Robotic Manipulators. IEEE/ASME Trans. Mechatronics 2019, 24, 1033–1040. [Google Scholar] [CrossRef]
  17. Jin, R.; Geng, Y.; Chen, X. Predefined-time control for free-floating space robots in task space. J. Frankl. Inst. 2021, 358, 9542–9560. [Google Scholar] [CrossRef]
  18. Liang, C.; Ge, M.; Liu, Z.; Ling, G.; Liu, F. Predefined-time formation tracking control of networked marine surface vehicles. Control Eng. Pract. 2021, 107, 104682. [Google Scholar] [CrossRef]
  19. Li, Y.; Hou, M.; Liang, S.; Jiao, G. Predefined-time adaptive fault-tolerant control of hypersonic flight vehicles without overparameterization. Aerosp. Sci. Technol. 2020, 104, 105987. [Google Scholar] [CrossRef]
  20. Chen, W.; Saif, M. An iterative learning observer-based approach to fault detection and accommodation in nonlinear systems. In Proceedings of the 40th IEEE Conference on Decision and Control, Orlando, FL, USA, 4–7 December 2001; pp. 4469–4474. [Google Scholar]
  21. Chen, W.; Saif, M. Observer-based fault diagnosis of satellite system subject to time-varying thruster faults. J. Dyn. Syst. Meas. Control 2007, 129, 352–356. [Google Scholar] [CrossRef]
  22. Yao, L.; Qin, J.; Wang, H.; Jiang, B. Design of new fault diagnosis and fault tolerant control scheme for non-Gaussian singular stochastic distribution systems. Automatica 2012, 48, 2305–2313. [Google Scholar] [CrossRef]
  23. Xiao, B.; Hu, Q.L.; Friswell, M.I. Active fault-tolerant attitude control for flexible spacecraft with loss of actuator effectiveness. Int. J. Adapt. Control Signal Process. 2013, 27, 925–943. [Google Scholar] [CrossRef]
  24. Cao, T.; Gong, H.J.; Cheng, P.; Xue, Y.X. A novel learning observer-based fault-tolerant attitude control for rigid spacecraft. Aerosp. Sci. Technol. 2022, 128, 107751. [Google Scholar] [CrossRef]
  25. Hu, Q.L.; Zhang, X.X.; Niu, G.L. Observer-based fault tolerant control and experimental verification for rigid spacecraft. Aerosp. Sci. Technol. 2019, 92, 373–386. [Google Scholar] [CrossRef]
  26. Zhang, X.B.; Zhou, Z.P. Integrated fault estimation and fault tolerant attitude control for rigid spacecraft with multiple actuator faults and saturation. IET Control Theory Appl. 2019, 13, 2365–2375. [Google Scholar] [CrossRef]
Figure 1. The structure of the observer-based predefined-time attitude control strategy.
Figure 1. The structure of the observer-based predefined-time attitude control strategy.
Processes 10 02294 g001
Figure 2. The estimation of fault.
Figure 2. The estimation of fault.
Processes 10 02294 g002
Figure 3. The estimation error of fault.
Figure 3. The estimation error of fault.
Processes 10 02294 g003
Figure 4. The estimation error of angular rate.
Figure 4. The estimation error of angular rate.
Processes 10 02294 g004
Figure 5. The response of attitude angle tracking performance.
Figure 5. The response of attitude angle tracking performance.
Processes 10 02294 g005
Figure 6. The response of angle velocities performance.
Figure 6. The response of angle velocities performance.
Processes 10 02294 g006
Figure 7. The response of actuator torque.
Figure 7. The response of actuator torque.
Processes 10 02294 g007
Table 1. Parameters and initial states for simulation.
Table 1. Parameters and initial states for simulation.
ParametersValue
The initial value of observer ω ^ ( 0 ) = ω ( 0 ) , e ^ ( 0 ) = [ 1 ; 1 ; 1 ]
Observer gains l = 0.9 , n = 2.5 , e m i n = 0.0001 .
M = [ 80 , 0 , 0 ; 0 , 50 , 0 ; 0 , 0 , 60 ] ,
H 1 = [ 8 , 0 , 0 ; 0 , 15 , 0 ; 0 , 0 , 12 ; 10 , 0 , 1.2 ] ,
H 2 = [ 20 , 0 , 0 ; 0 , 3 , 0 ; 0 , 0 , 8 ; 1 , 0.5 , 6.5 ] .
Controller parameters h 1 = 0.5 , h 2 = 0.3 , k s = 5 ,
T 1 = 10 , T 2 = 10 .
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cao, T.; Gong, H.; Han, B. Observer-Based Predefined-Time Attitude Control for Spacecraft Subject to Loss of Actuator Effectiveness. Processes 2022, 10, 2294. https://doi.org/10.3390/pr10112294

AMA Style

Cao T, Gong H, Han B. Observer-Based Predefined-Time Attitude Control for Spacecraft Subject to Loss of Actuator Effectiveness. Processes. 2022; 10(11):2294. https://doi.org/10.3390/pr10112294

Chicago/Turabian Style

Cao, Teng, Huajun Gong, and Bing Han. 2022. "Observer-Based Predefined-Time Attitude Control for Spacecraft Subject to Loss of Actuator Effectiveness" Processes 10, no. 11: 2294. https://doi.org/10.3390/pr10112294

APA Style

Cao, T., Gong, H., & Han, B. (2022). Observer-Based Predefined-Time Attitude Control for Spacecraft Subject to Loss of Actuator Effectiveness. Processes, 10(11), 2294. https://doi.org/10.3390/pr10112294

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop