In order to deal with the low convergence rate of the standard sliding mode in satellite attitude control, a novel variable structure sliding mode is constructed in this paper by designing the update law of the sliding mode parameter. By implementing this method, the advantage such as simple structure and strong robustness of the standard sliding mode are maintained and the system convergence rate is largely improved. The fixed sliding mode parameter is modified, and the update law is designed. When the system state is away from the sliding mode surface, the parameter is fixed, and when the system state approaches the sliding mode surface, the parameter begins to update. The constraint on control torque and angular velocity is taken into consideration, and the constraint on control parameters is given to ensure that the system state do not exceed its upper bound. System stability is proved by the Lyapunov stability theorem, and the superiority of the proposed controller is demonstrated by numerical simulation.
National Natural Science Foundation of China619032891. Introduction
In the satellite attitude control issue, sliding mode control has drawn much attention because of its simple structure and clear physical meaning. Researchers have done a lot of work focusing on sliding mode controllers for satellite attitude control. Xiao [1, 2] and Hu [3, 4] have done some work for satellite attitude control, and some typical sliding modes are proposed in their work. The focus of their work is system robustness for model uncertainty, unknown disturbance, and actuator failure. Wu [5, 6] also proposed some methods dealing with the chattering issue in sliding mode controllers. Most of these works mentioned above pointed out that strong robustness is an advantage of the standard sliding mode controller, but the exponential convergence rate (i.e., low convergence rate issue) near the system equilibrium point is its main drawback.
In order to improve the system convergence rate governed by sliding mode controllers, researchers developed some fast convergence control methods. Liang [7, 8] and Wang [9, 10] developed finite-time controllers for satellite attitude control. The focus of their work is the structure of the finite-time sliding mode, and some typical finite-time sliding modes are proposed. It is pointed out that the finite-time controller could largely improve the system convergence rate compared with the standard sliding mode but the cost is the complex structure of controllers which causes the loss of strong robustness. Wu [11] developed a finite-time controller based on the fast terminal sliding mode surface. The terminal sliding mode proposed by Dr. Wu in Ref [11] with finite time convergence rate is combined with standard sliding mode and a fractional order sliding mode. Hence, the system could have expected performance. Li [12] proposed a finite-time controller with a three-stage structure. A maneuver stage with constant angular velocity is constructed, and the angular velocity capability is fully utilized. Generally, finite-time control is an effective method to improve the system convergence rate but the cost is its complex structure and also the loss of strong robustness.
In order to have both fast convergence rate and strong robustness, researchers have done some modifications on the standard sliding mode. Jovan [13] proposed a simple structure sliding mode controller with the time varying parameter. System robustness is largely improved by the construction of the parameter update law. Although his main focus is on system robustness but not the convergence rate, the time varying sliding mode method proposed in his work is insightful. Verbin [14] and Lappas [15] also have done some work for satellite attitude control. They pointed that the key to improve the system convergence rate is to design the angular velocity trajectory properly. The braking curve for angular velocity is proposed in their work, and system performance is largely improved. Li [16–18] proposed some methods to design the angular velocity trajectory, and the system convergence rate is largely improved when compared with traditional methods. Sun [19–21] developed an improved direct instantaneous torque control based on the adaptive terminal sliding mode for a segmented-rotor SRM. A novel reaching law is proposed to reduce torque ripple and startup time. The modified sliding mode controllers are proposed, and system performance is largely improved. These works mentioned above could largely improve the system convergence rate, but the structure of the controller and sliding mode is still relatively complex and is hard for engineering practice.
Noting that the control capability is limited in engineering practice, excessive angular velocity would bring some unexpected characters (such as the failure of the star sensor, drop of attitude determination accuracy, large flexible deformation, and fuel consumption for reaction wheel damping); hence, it is necessary to consider control torque and angular velocity constraint when designing controllers. Hence, the goal of this paper could be concluded as follows: develop a sliding mode controller with simple structure and strong robustness under system control torque and angular velocity constraint. Also, the highlight and the main contribution of this paper could be concluded as follows: the standard sliding mode is modified by implementing a dynamic parameter update law; hence, the system convergence rate is improved, and this is a new method by the design sliding mode with expecting characters.
The structure of this paper could be listed as follows: the first part introduces the necessity of this paper, and this second part gives the dynamic and kinetic model used in this paper. The 3–5 parts are the main contributions of this paper, 3rd part gives the controller, 4th part discusses the constraint on control parameters, and 5th part gives the numerical simulation and its analysis. The 6th part concludes this paper.
2. Dynamic and Kinetic Model
The dynamic model of the satellite could be written as follows:(1)Jω˙+ω×Jω=u+d,where ω is a three-dimensional vector which describes system angular velocity, J is the inertia matrix of the satellite which is a 3×3 symmetric matrix, d is three-dimensional unknown disturbance torque, and u is a three-dimensional vector and the control torque generated by actuators such as thruster and reaction wheel, and the design of u is the main focus of this paper. Product matrix r× of vector r is defined as(2)r×=0−r3r2r30−r1−r2r10.
Generally, inertia matrix J could not be known accurately, and it could be described as follows:(3)J=J^+J˜,where J^ is the inertia matrix best estimate and J˜ is the error matrix.
Product matrix has an important property which will be used in the later part that the eigenvalues of r× satisfie(4)λr×=0,±r,λmaxr×=r.
The kinetic model based on attitude quaternion could be written as follows:(5)q=q0qvTT=cosφ2eTsinφ2T,q0=cosφ2,qv=esinφ2,q˙0=−12qvTω,q˙v=12q0I3+qv×ω=12Fω,where q0 and qv are the scalar part and vector part of attitude quaternion, respectively, with q0≤1,qvTqv+q02=1 and e and φ are the Euler axis and Euler angle, respectively, and the definition is given in [12], and the eigenvalue of F satisfies(6)λF=q0,1,λmaxF=1.
In order to simplify the text, the maximum and minimum eigenvalue of matrix A is described as λMA and λmA.
3. Variable Structure Sliding Mode Controller
In the satellite attitude control field, the standard sliding mode surface could be written as follows:(7)s=ω+kqv,k>0.
When the system state reaches this sliding mode surface, we have(8)ω=−kqv,and(9)q˙0=−12ωTqv=12kqvTqv=12k1−q02.
Based on equation (9), it is obvious that the system convergence rate is exponential. It is worth noticing that the reverse of angular velocity and attitude quaternion is an important property of this sliding mode surface. This property implies that when the norm of angular velocity is fixed, this sliding mode surface has the best efficiency on angular velocity, and because of this property, the standard sliding mode has been developed by researchers.
However, the main drawback of the standard sliding mode surface is its low convergence rate when the system state approaches its equilibrium point. The sliding mode parameter k is fixed, and the norm of angular velocity is proportional to attitude quaternion, and when approaching the equilibrium point, the norm of angular velocity descends drastically, and this caused the low convergence rate of the standard sliding mode. If better convergence rate is needed, a common method is to enlarge the parameter, but this method would cause some unexpected characters such as severe vibration, large control torque, and angular velocity.
In order to solve the low convergence rate issue and avoid those drawbacks mentioned above, a variable structure sliding mode is proposed in this paper as follows:(10)s=ω+kqv,k0>0,k˙=0,s≥ε,12pq031−q0,s<ε,where k0 is the initial value of sliding mode parameter k and is a positive scalar, p is a positive scalar that needs to be designed properly, and ε is a small positive scalar.
It is obvious that sliding mode surface (10) is constructed as two stages: when the system state is away from the sliding mode surface i.e., s≥ε, the proposed sliding mode is totally same as the standard sliding mode surface, and the parameter is fixed; when the system state reaches the sliding mode surface, i.e., s<ε, the angular velocity has been reversed to attitude quaternion and the sliding mode parameter begins to update. Based on the update law, it could be found that its derivative is positive scalar, and this means that the system could maneuver along the sliding mode surface with a small initial parameter, and the sliding mode parameter is enlarged when the system state approaches its equilibrium point. By this method, the norm of angular velocity would not drop so drastically as of the standard sliding mode; hence, the system convergence rate would be improved. It is worth noticing that on sliding mode surface (10), the angular velocity is reversed to attitude quaternion, and this means the advantage of the sliding mode could be maintained. Also, it could be found that when the system converges to its equilibrium point, the derivative of sliding mode parameters also tends to zero; hence, update law (10) does not have the singularity issue.
The attitude controller based on sliding mode (10) could be written as follows:(11)u=−ρ1kss+ω×J^ω−12kJ^Fω−l1sgns,s≥ε,−ρ2kss+ω×J^ω−12kJ^Fω−12pq031−q0J^qv−l2sgns,s<ε,where(12)l1=d¯+λω2+12kω,l2=d¯+λω2+12kω+12pq031−q0qv.
In equation (11), ks is a positive scalar, λ is a positive scalar and satisfies λ≥λMJ˜, and ρi is the control gain factor defined as follows:(13)ρi=−kssTri+kssTri2+ks2s2u¯2−ri2ks2s2,u≥u¯,1,u<u¯,r1=ω×J^ω−12kJ^Fω−l1sgns,r2=ω×J^ω−12kJ^Fω−12pq031−q0J^qv−l2sgns,where u¯ is the norm upper bound of control torque. The goal of ρi is to ensure that the control torque does not exceed its upper bound. Based on its definition, it could be found that when the demand torque does not exceed its upper bound, the factor does not change the controller structure, and when the demand torque exceeds its upper bound, ρi could reduce the proportion term.
The next step is to prove system stability under controller (11). Select the Lyapunov function as follows:(14)Vs=12sTJs.
When the system state does not reach the sliding mode, i.e., s≥ε, calculate the derivative of (14) and substitute controller (11), and it could be obtained that(15)V˙s=sTJs˙=sTJω˙+ksTJq˙v,=sTu+d−ω×Jω+12ksTJFω,=−ρ1kssTs+sTω×J^ω−12ksTJ^Fω−l1sTsgns+sTd−sTω×Jω+12ksTJFω,=−ρ1kssTs−sTω×J˜ω+12ksTJ˜Fω+sTd−d¯+λω2+12kωsTsgns,≤−ρ1kssTs≤0.
When the system state reaches the sliding mode, i.e., s<ε, calculate the derivative of (14) and substitute controller (11), and it could be obtained that(16)V˙s=sTJs˙=sTJω˙+ksTJq˙v+k˙sTJqv,=sTu+d−ω×Jω+12ksTJFω+12pq031−q0sTJqv,=−ρ2kssTs+sTω×J^ω−12ksTJ^Fω−12pq031−q0sTJ^qv−l2sTsgns+sTd−sTω×Jω+12ksTJFω+12pq031−q0sTJqv,=−ρ2kssTs+sTd−sTω×J˜ω+12ksTJ˜Fω+12pq031−q0sTJ˜qv−d¯+λω2+12kω+12pq031−q0qvsTsgns,≤−ρ2kssTs≤0.
Hence, the sliding mode state s is globally asymptotically stable, and based on the discussion above, systems (1) and (5) are globally asymptotically stable on sliding mode surface (10); hence, systems (1) and (5) governed by controller (11) are globally asymptotically stable. System stability has been proved.
4. System State Constraint
Generally, in order to get better convergence rate, the initial value of sliding mode parameter k0 and parameter p should be selected as large as possible. However, it is obvious that larger control parameter would cause larger control torque and angular velocity, and system performance would be harmed, and the system state would diverge under some extreme cases. Hence, it is necessary to discuss the relation between angular velocity upper bound ω¯, control torque upper bound u¯, and control parameters.
Firstly, the constraint of control torque u¯ on control parameters is given. Noting that when the system state does not reach the sliding mode, the proportion term kss could be reduced by the gain factor ρi; hence, the first term of controller (11) does not have the control torque saturation issue, and the key to satisfy the control torque constraint is the later terms. Noting that when s≥ε,(17)r1=ω×J^ω−12kJ^Fω−l1sgns,≈ω×Jω+12kJFω,≤λMJω¯2+12k0ω¯,where ω¯ is the norm upper bound of angular velocity; hence, if the following inequality is satisfied, the control torque constraint could be satisfied:(18)λMJω¯2+12k0ω¯≤u¯.
When the system state reaches the sliding mode surface, the angular velocity satisfies(19)ω=−kqv.
Define auxiliary variable z as follows:(22)z=λMJ32k2+118pqv.
It is obvious that z is the norm upper bound of r2; hence, if z≤u¯ is satisfied, r2≤u¯ could be ensured. Calculate the derivative of (22) and it could be obtained that(23)z˙=λMJ3kk˙qv+32k2+118pdqvdt,=12λMJkq0qv3pq021−q0−32k2−118p.
Hence, if z reaches its maximum value, the following inequality holds:(24)k2=2pq021−q0−127p.
Since the following inequality holds, the control torque constraint could be ensured:(26)328λMJp≤u¯.
Above all, the constraint of control torque on control parameters could be concluded as follows:(27)λMJω¯2+12k0ω¯≤u¯,328λMJp≤u¯.
The former inequality is the constraint on the initial value of the sliding mode parameter, and the later inequality is the constraint on the update parameter. If they both are satisfied, system control torque would not exceed its norm upper bound during the whole control process.
The next step is to discuss the constraint of angular velocity on control parameters. Noting that, before reaching the sliding mode surface, the sliding mode structure is totally same as the standard sliding mode, and the method on the standard sliding mode could be implemented. Noting that, under this condition system could be treated as a two-stage linear system, and the overshot of the system state generally is usually less than 100%; hence, by selecting control parameters as follows, the constraint of angular velocity could be ensured:(28)k0≤12ω¯.
When it reaches the sliding mode surface, calculate the derivative of the angular velocity norm:(29)dωdt=k˙qv+kdqvdt,=12pq031−q0qv−12k2q0qv,=12q0qvpq021−q0−k2.
When the norm of angular velocity reaches its maximum value,(30)k2=pq021−q0.
Hence, we have(31)ω=kqv≤pq021−q01−q02≤2pq021−q02≤2p4.
Hence, the constraint of angular velocity on control parameters could be concluded as follows:(32)k0≤12ω¯,2p4≤ω¯.
The former inequality is the constraint on the initial value of the sliding mode parameter, and the later inequality is the constraint on the update parameter. If they are both satisfied, system angular velocity would not exceed its norm upper bound during the whole control process.
5. Numerical Simulation and Analysis
Set system configurations as follows:(33)J=diag302520kg·m2,J^=diag222421kg·m2,λ=2,tsample=0.5s,d¯=0.001N·m,u¯=0.5N·m,ω¯=0.1rads,ωt0=−0.03−0.040.05Trads,qt0=0663322T.
In order to demonstrate the superiority of the proposed controller, the standard sliding mode controller is compared as follows:(34)s=ω+kqv,u=−kss+ω×Jω−12kJFω−d¯sgns.
Case 1.
Standard sliding mode controller with the large parameter.
Generally, larger control parameters in (34) could bring better a convergence rate. In order to select control parameters properly, the method to determine control parameters in [12] is referred in this paper, and following parameters are selected:(35)ks=2,k=0.25.
The simulation results are shown as follows:
Based on Figures 1 and 2, it could be found that system convergence time is about 90 s; however, Figure 3 demonstrates that the initial control torque is relative large, and there exists a 10 s period with control torque exceeding the system norm upper bound. Also, based on Figure 3, it could be found that there exists a 15 s period with angular velocity exceeding the system upper bound. This group of simulation demonstrates that the standard sliding mode could have a good convergence rate by selecting the large parameter but the cost is excessive control torque and angular velocity.
Curve of angular velocity.
Curve of attitude quaternion.
Curve of angular velocity and control torque norm.
Case 2.
Standard sliding mode controller with the small parameter.
Based on the simulation results above, it could be found that system angular velocity exceeds its upper bound; hence, following parameters are selected:(36)ks=2,k=0.1.
The simulation results are shown as follows:
In this group of simulation, Figures 4 and 5 demonstrate that system convergence time is about 150 s, and Figure 6 demonstrates that system control torque and angular velocity are both lower than the system upper bound. Comparing this group with Case 1, it could be easily found that if control torque and angular velocity constraints are satisfied by selecting the small parameter, but the convergence rate of the standard sliding mode is relatively low.
Curve of angular velocity.
Curve of attitude quaternion.
Curve of angular velocity and control torque norm.
Case 3.
Dynamic sliding mode proposed in this paper.
Based on previous simulation, it could be found that system performance is mainly determined by k, and ks determines the initial control torque. Also, based on the update law, it could be found that larger p and r should be selected and the constraints should be satisfied. Hence, the following parameters are selected:(37)ks=2,k0=0.05,p=0.5,r=0.15.
Based on equations (27) and (32), it could be found that the constraints on control parameters are all satisfied with equation (36).
The simulation results are shown as follows (Table 1).
Based on Figures 7 and 8, it could be found that system convergence time is about 80 s which is at the same level with Case 1. Comparing simulation results of Case 3 with that of Case 1 and Case 2, it could be easily found that the sliding mode parameter in Case 3 is the same as Case 2 (k = 0.05), and the convergence time in Case 3 is about half of Case 2. And, the sliding mode parameter in Case 3 is 20% of Case 1 (k(0) = 0.05 in Case 3 and k = 0.25 in Case 1), and the convergence time is at the same level. This demonstrates that, by implementing the update law of the sliding mode parameter, system convergence time could be largely improved.
Based on Figure 9, it could be found that system control torque and angular velocity are lower than the system upper bound during the whole control torque; however, in Case 1, this property does not hold. This demonstrates that the proposed method could improve the system convergence rate, while the demand control torque and angular velocity are relatively small compared with the standard sliding mode. Also, this demonstrates the effectiveness of the constraints on control parameters (27) and (32) proposed in this paper.
Based on Figure 10, it could be found that sliding mode parameter k is norm upper bounded, and the dynamic sliding mode proposed in this paper does not have the singularity issue. Also, it could be found that the basic structure of the standard sliding mode is maintained; hence, the system has strong robustness to model uncertainty and unknown disturbance torque.
Above all, it could be concluded that the proposed method could improve the system convergence rate without causing control torque and angular velocity saturation issue. The method proposed in this paper could deal with the low convergence rate issue of the standard sliding mode and maintain its advantages.
Simulation results comparison.
Convergence time (s)
Max control torque
Max angular velocity
Conclusion
Standard sliding mode with larger parameter
90
1.3 Nm
0.18 rad/s
Best convergence rate, but control torque and angular velocity exceed system upper bound
Standard sliding mode with larger parameter
150
0.27 Nm
0.08 rad/s
Torque and velocity constraints are satisfied, but convergence rate is the slowest
Dynamic sliding mode proposed in this paper
80
0.5 Nm
0.08 rad/s
Same convergence rate as Case 1, and torque and velocity constraints are satisfied
Curve of angular velocity.
Curve of attitude quaternion.
Curve of angular velocity and control torque norm.
Curve of sliding mode parameter k.
6. Conclusion
In this paper, a dynamic sliding mode with the time varying sliding mode parameter is proposed for satellite attitude control based on the standard sliding mode. The focus of the proposed method is to deal with the low convergence rate issue and maintain other advantages of the standard sliding mode. It is also proved that the designing of angular velocity trajectory properly is the key to improve convergence rate and system performance. The results of this paper are listed as follows:
The reverse of angular velocity and attitude quaternion is the key to improve the system convergence rate, and the low convergence rate is caused by the descend of angular velocity in the standard sliding mode. By enlarging the norm of angular velocity and maintaining the reverse property, the system convergence rate could be improved.
By implementing the update law of the sliding mode parameter, the system could have better performance compared with the fixed parameter sliding mode. System control torque and velocity are utilized more effectively by the method of variable structure sliding mode. Moreover, by designing the update law properly, the system could both have good convergence rate and strong robustness.
System control torque and angular velocity constraints could be satisfied by selecting control parameters properly, and the method proposed in this paper is suitable for engineering practice.
The dynamic sliding mode method proposed in this paper is enlightening for the sliding mode controller design, and later works could focus on the dynamic sliding mode with finite-time or fixed-time convergence rate.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest regarding the publication of this paper..
Acknowledgments
This work was supported partially by National Natural Science Foundation of China (Project no. 61903289). The authors greatly appreciate the above financial support.
XiaoB.YinS.A structure simple controller for satellite attitude tracking maneuver201764214361446XiaoB.YinS.GaoH. J.Tracking control of robotic manipulators with uncertain kinematics and dynamics2016631064396449HuQ. L.DuC. L.XieL. H.WangY. Y.Discrete-time sliding mode control with time-varying surface for hard disk drives2009171175183HuQ. L.ZhangY. M.HuoX.XiaoB.Adaptive integral-type sliding mode control for spacecraft attitude maneuvering under actuator stuck failures20112413245WuS. N.WangR.WuZ. G.Robust attitude maneuver control of spacecraft with reaction wheel low-speed friction compensation201543218213WuS. N.SunZ. W.DengH.Robust sliding mode controller design for globally fast attitude tracking of target spacecraft20117520872098LiangH. Z.WangJ. Y.SunZ. W.Robust decentralized coordinated attitude control of spacecraft formation2011695-6280288LiangH. Z.SunZ. W.WangJ. Y.Finite-time attitude synchronization controllers design for spacecraft formations via behavior-based approach20131117371753WangJ. Y.LiangH. Z.SunZ. W.Dual-quaternion-based finite-time control for spacecraft tracking in six degrees of freedom20133528545WangJ. Y.LiangH. Z.SunZ. W.ZhangS. J.Finite-time control for spacecraft formation with dual-number-based description2012353950962WuS. N.RadiceG.SunZ. W.Robust finite-time control for flexible spacecraft attitude maneuver2014271185190LiY.YeD.SunZ.Robust finite time control algorithm for satellite attitude control201768465710.1016/j.ast.2017.05.0142-s2.0-85019394775JovanD. B.SaiM. L.Robust tracking control design for spacecraft under control input saturation2004274627633VerbinD.LappasV. J.Rapid rotational maneuvering of rigid satellites with hybrid actuators configuration2013362532547VerbinD.Rapid rotational maneuvering of rigid satellites with reaction wheels20133615381544LiY.YeD.SunZ.Time efficient sliding mode controller based on bang-bang logic for satellite attitude control20187534235210.1016/j.ast.2017.12.0422-s2.0-85042058173LiY.YeD.Near time-optimal controller based on analytical trajectory planning algorithm for satellite attitude maneuver201984497509LiY.SunZ. W.YeD.Time efficient robust PID plus controller for satellite attitude stabilization control considering angular velocity and control torque constraint201730510.1061/(asce)as.1943-5525.00007432-s2.0-85018904116SunX. D.FengL. Y.DiaoK.YangZ. B.An improved direct instantaneous torque control based on adaptive terminal sliding mode for a segmented-rotor SRM202099SunX. D.CaoJ. H.LeiG.GuoY. G.A composite sliding mode control for SPMSM drives based on a new hybrid reaching law with disturbance compensation202199SunX. D.WuJ. L.LeiG.GuoY. G.Torque ripple reduction of SRM drive using improved direct torque control with sliding mode controller and observer202099