IEEE/CAA Journal of Automatica Sinica  2016, Vol.3 Issue (1): 42-50   PDF    
Characteristic Model-based Discrete-time Sliding Mode Control for Spacecraft with Variable Tilt of Flexible Structures
Lei Chen, Yan Yan, Chaoxu Mu, Changyin Sun     
1. School of Automation, Southeast University, Nanjing 210096, China;
2. Department of Electrical Engineering and Automation, Tianjin University, Tianjin 300072, China
Abstract: In this paper, the finite-time attitude tracking control problem for the spacecrafts with variable tilt of flexible appendages in the conditions of exogenous disturbances and inertia uncertainties is addressed. First the characteristic modeling method is applied to the problem of the spacecraft modeling. Second, a novel adaptive sliding mode surface is designed based on the characteristic model. Furthermore, a discrete-time sliding mode control (DTSMC) law, which makes the tracking error converge into a predefined bound in finite time, is proposed by employing the parameters of characteristic model associated with the sliding mode surface to provide better performances, robustness, faster response, and higher control precision. The designed DTSMC includes the adaptive control architecture and is chattering-free. Finally, digital simulations of a sun synchronous orbit satellite (SSOS) are presented to illustrate effectiveness of the control strategies as well as to verify the practical feasibility of the rapid maneuver mission.
Key words: Spacecraft     attitude maneuver     characteristic model     sliding mode control     finite time    
Ⅰ. INTRODUCTION

Many future space missions will involve spacecrafts with flexible appendages such as large antennas,sun shields,solar arrays,and solar sails. These complex flexible nonlinear systems (CFNSs) are the systems with both nonlinear and flexible characteristics. Compared with the rigid body spacecrafts without flexible appendages,the design and implementation of attitude and orbit control systems (AOCSs) for the spacecrafts with such structures involve significant difficulties,because of the combination of uncertainty on several vehicle parameters,the frequency and damping of each flexible mode,and the potentially significant coupling between axes[1, 2, 3]. At the early stage,researchers developed the low-dimensional models for the CFNSs by ignoring the impacts of flexible modalities or considering flexible modalities as external disturbances. To improve the performances of the close-loop systems, several nonlinear control techniques based on those models such as compound control[4],$H_\infty$ state-feedback control[5], geometric control[6],etc.,have been proposed for solving the attitude tracking control problem of spacecrafts. Furthermore,a finite-time control theory[7] has been designed to provide fast convergence. However,in terms of attitude tracking tasks,large angle or fast angular rate is generally considered in the rotation process. Since these results do not consider flexible characteristics of the CFNSs,the vibration in the flexible structures introduced by the orbiting attitude slewing operation may degrade the attitude pointing accuracy. This is the reason why the flexible appendages cannot be ignored anymore,especially when,in some cases,they are also required to rotate[8].

Currently,the spacecraft dynamic model is described as a relationship between the 6 DOF vector of acceleration at the composite center of the mass. Barring the small angle and angular rate assumptions,the commonly adopted linearized model is not suitable to design attitude control laws. Moveover,according to the modal analysis theory,the order of the flexible appendages is usually enormously high leading to that it is impossible to depict the CFNSs with accurate mathematical models. Generally,the modal truncation method is used to reduce the order. However,to cut off the order means to lose some information of the CFNSs,which leads to some unpredictable challenges. Taking the Hubble telescope launched in 1990 for an example[9],a PID controller was designed with this method. But in space,frequency sympathetic vibration caused by disturbance occurred,which made the system performance worse.~~Several~~recent~~studies have considered the application of linear robustness analysis techniques to this problem,but significant difficulties were also encountered with promising results reported,such as the extremely complicated problem of computing tight upper and lower bounds on singular value[10]. One important reason for the difficulties in modeling and control is that the present modeling and control theory is based on accurate dynamic analysis as well as modeling and control are considered separately. As a result,the current modeling method needs to be modified by an engineering modeling method,i.e., a characteristic modeling method[9]. The theory of characteristic modeling,as a kind of data-based model-free approach for control system design,requires the analysis of the dynamic and control performance of the system instead of an accurate plant dynamic analysis. Based on this methodology,a golden section adaptive control law[11],which aims at designing an engineering-oriented adaptive controller by using only a few parameters,has been successfully used in various industries.

Sliding mode control (SMC)[12] has been widely studied and used for many years due to its simplicity and robustness in system variations and disturbances. The idea of the SMC is the peculiar sliding motion introduced by discontinuous control,which enables the system trajectory to reach and stay in a prescribed switching manifold indefinitely. It is also well known that SMC is an efficient method to deal with nonlinear systems with exogenous disturbances and parametric uncertainties[13]. Since modern industrial control systems are digitally implemented[14],which means that the industrial SMC would be executed in discrete time, the need for research in DTSMC is evident. However,even though the long history of SMC development and main results have been reported, only a few of them among available results are devoted to discrete-systems[15]. It should be pointed out that DTSMC cannot be obtained from its continuous counterpart by means of simple equivalence,theoretically. Although the investigation of the discretization behaviors of terminal SMC systems[16] revealed the underpinning generation mechanism periodic orbits,which will help develop effective discrete-time terminal SMC with less chattering and better robustness,it has not been gotten fully used for actual controlled plants.

In this paper,by introducing a novel sliding mode surface based on the characteristic modeling,a finite-time DTSMC method is proposed for the tracking control problem of the rigid spacecraft with flexible structures. Firstly,according to the characteristic modeling method,the characteristic model of the CFNSs is established and the parameters of the model are obtained with identification algorithm. Secondly,a characteristic model-based DTSMC law is designed and Lypapunov stability law is applied to tackle the problem how states can reach the bounded sliding surface in finite time. Finally,the proposed controller is used for the attitude tracking control problem of the CFNSs.

This paper is organized as follows: Section II describes the 6 DOF spacecraft simulation model,AOCS and characteristic modeling criteria used in this study. The main researches are presented in Section III,in which the characteristic model-based DTSMC is proposed to achieve the highly accurate tracking performance. Simulation results are provided in Section IV. Finally,some conclusions are presented in Section V.

Ⅱ. SPACECRAFT MODEL AND INVOLVED CRITERIA

A. Spacecraft Attitude Dynamics and Kinematics

The main body model is obtained thanks to Euler/Netwon equations applied to the rigid body. Considering the worst case,the corresponding uncertainty around their nominal value is 50 %. The dynamic model of the rigid body is

$ \begin{align} \boldsymbol{I}_s \dot{\boldsymbol{\omega}}_s + \boldsymbol{\omega}_s^{\times} \boldsymbol{I}_s \boldsymbol{\omega}_s + \boldsymbol{F}_{sa}\ddot{\boldsymbol{\eta}}_{sa} = \boldsymbol{T}_T, \end{align} $ (1)

where $\boldsymbol{\omega}_s = \left[\omega_x,\ \omega_y,\ \omega_z \right]^{\rm T} \in R ^{3 \times 1}$ is the angular velocity, $\boldsymbol{I}_s$ $\in$ $R ^{3 \times 3}$ is the rotational inertia of the spacecraft,$\boldsymbol{T}_T \in R ^{3 \times 1}$ is the aggregate external torque,$\boldsymbol{\eta}_{sa} \in R^{m \times 1}$ is the vibration modal coordinates array of the flexible appendage,$\boldsymbol{F}_{sa}$ $\in$ $R^{3 \times m}$ is the coupling coefficient matrix between the flexible appendage and the rigid body,and $(\cdot)^{\times}$ denotes a $3 \times 3$ skew-symmetric matrix,that is

$ \begin{align}\notag \boldsymbol{a}^{\times} = \left[\begin{array}{ccc}0& -a_3&a_2\\ a_3&0&-a_1\\ -a_2& a_1&0\end{array} \right]. \end{align} $

The dynamic model of the flexible appendage cantilevered on the main body is commonly described by the so-called cantilever hybrid model[17]. This model gives a relationship between the 6 DOF acceleration vector and the 6 DOF forces vector applied by the main body to the appendage:

$ \begin{align} \ddot{\boldsymbol{\eta}}_{sa} + 2 \xi_{sa} \boldsymbol{\omega}_{sa} \dot{\boldsymbol{\eta}}_{sa} +\boldsymbol{\omega}_{sa}^2 \boldsymbol{\eta}_{sa} + \boldsymbol{F}_{sa} ^{\rm T} \dot{\boldsymbol{\omega}}_s = \boldsymbol{{0}}, \end{align} $ (2)

where $\boldsymbol{\omega}_{sa} \in R^{m \times m}$is the modal frequency matrix of the flexible appendage vibration,and $\xi_{sa}$ is the vibration damping ratio of appendage.

Quaternion is generally used in the onboard attitude presentation, which is defined by $q_0=\cos(\frac{\alpha}{2})$,$\boldsymbol{q}=$ $[\gamma_1 \sin(\frac{\alpha}{2})$, $\gamma_2\sin(\frac{\alpha}{2})$, $\gamma_3\sin(\frac{\alpha}{2})]^{\rm T}$,where $[\gamma_1,\gamma_2,\gamma_3]$ is the principle axis from the current attitude to the reference attitude and $\alpha$ is the principle angle. The kinematics model is established referring to the aforementioned frame and quaternion:

$ \begin{align} \dot{\boldsymbol{q}} = \frac{1}{2} \left( \boldsymbol{q}^{\times} + q_0 \boldsymbol{I} \right) \boldsymbol{\omega}_s,\ \ \ \dot{q}_0= -\frac{1}{2} \boldsymbol{q}^{\rm T} \boldsymbol{\omega}_s. \end{align} $ (3)

Considering the tilt of the cantilevered appendage,the varying angle $\theta_{sa}$ affects all terms of the spacecraft dynamic model. Due to this problem,the coupling matrix between the flexible appendage and the rigid body is time varying as follows:

$ \begin{align} \boldsymbol{F}_{sa}=\left[ \begin{array}{ccc} \cos(\theta_{sa}) & 0 & -\sin(\theta_{sa}) \\ 0 & 1 & 0\\ \sin(\theta_{sa}) & 0 & \cos(\theta_{sa}) \end{array} \right] \boldsymbol{F}_{s0}. \end{align} $ (4)

B. Characteristic Modeling

Modeling based on plant dynamic characteristics and performance requirements rather than only accurate plant dynamic analysis is the key idea of the characteristic modeling. Unlike other intelligent modeling methods such as the T-S fuzzy modeling[18],the characteristic modeling compresses all information of the system into several characteristic parameters. A characteristic model has the following features[19].

1) A plant characteristic model is equivalent to its practical plant in output for the same input,i.e.,the output error can maintain within a permitted range in a dynamic process as well as their outputs are equal in the steady state.

2) The order and form of a characteristic model mainly depends on control performance requirements.

3) Compared with an original dynamic equation,the structure of a characteristic model should be simpler,easier,and more convenient for realization in engineering.

4) Unlike the reduced-order model of a high-order system,a characteristic model compresses all the information of the high-order system into several characteristic parameters,which means no information is lost. Generally,a characteristic model is represented by a slowly time-varying difference equation.

Consider the nonlinear system

$ \begin{align} \dot{x}(t) = f \left( x,\dot{x},...,x^{(n)},u,\dot{u}, ...,u^{(m)}\right), \end{align} $ (5)

where $x$ and $u$ denote the state and input of the system, respectively. Choosing

$ \begin{align} \begin{split} &x = x_1,~~ \dot{x} = x_2,~~ ...,~~ x^{(n)} = x_{n+1},\\ &u = u_1,~~ \dot{u} = u_2,~~ ...,~~ u^{(m)} = u_{m+1}, \end{split} \end{align} $ (6)

then (5) can be rewritten as

$ \begin{align} \dot{x}_1(t) = f( x_1,\ ...,\ x_{n+1},\ u_1,\ ...,\ u_{m+1}). \end{align} $ (7)

Assume that the properties of the nonlinear system (7) are as follows[20]:

a) There is only a single input and a single output.

b) The power of $u(t)$ is 1.

c) If $x_i=0$ and $u_i=0$,we have $f(\cdot)=0$.

d) $f(\cdot)=0$ is continuously differentiable with all variables, and partial derivative values are bounded.

e) $| f ( x(t+\triangle t),u(t+\triangle t) ) | - | f( x(t),u(t) ) | < M \triangle t$,where the constant $M > 0$ and $\triangle t$is sampling time.

f) The states and control value are bounded.

Based on (7),the following lemma will be used in the derivation of the main results.

Lemma 1[20]. For any nonlinear system that can be described as (7),if assumptions a)-d) are satisfied and sampling time $\triangle t$ satisfies certain conditions,the desired signal $r$ and its derivative value are bounded,the error characteristic model of the system can be expressed with the following 2-order difference equation:

$ \begin{align} e(k+1)=&\ f_1(k)e(k)+f_2(k)e(k-1)\notag\\ & +g_0(k)u(k)+g_1(k)u(k-1), \end{align} $ (8)

where $g_0(k)=-g_1(k)+ {O} (\triangle t)$,$ {O} (\triangle t)$ represents the high-order infinitesimal term of the sampling time,$e(k)=x(k)-r(k)$,and $u(k)$ is the bounded sampling control. If the system is stable and assumptions e)-f) are satisfied,then

1) $f_1(k)$,$f_2(k)$,$g_0(k)$,and $g_1(k)$ are slowly time varying.

2) The ranges of these coefficients can be determined beforehand.

3) In dynamic process,under the same input,selecting suitable sampling period $\triangle t$ can make sure that the output error between the characteristic model and the controlled plant keeps within a permitted limit.

4) In steady state,both outputs are equal.

For the minimum-phase system,in general the error characteristic model is chosen as follows:

$ \begin{align} e(k+1)=f_1(k)e(k)+f_2(k)e(k-1)+g_0(k)u(k). \end{align} $ (9)

Ⅲ. SPACECRAFT ATTITUDE CONTROL

A. The Characteristic Model of the Spacecraft

Due to the attitude kinematics and dynamics of (1),(2),and (3), defining $\boldsymbol{x}=\left[q_0,\boldsymbol{q}^{\rm T},\boldsymbol{\eta}_{sa}^{\rm T}, \dot{\boldsymbol{\eta}}_{sa}^{\rm T},\boldsymbol{\omega}^{\rm T}_s\right]^{\rm T}$ and $\boldsymbol{u}=\boldsymbol{T}_T$,the attitude kinematics and dynamics of the CFNSs can be rewritten as:

$ \begin{align} \boldsymbol{A}_1\dot{\boldsymbol{x}}=\boldsymbol{A}_2(\boldsymbol{x})+\boldsymbol{Bu}, \end{align} $ (10)

where

$ \begin{align} {A}_1=\left[\begin{array}{ccccc} 1 & \boldsymbol{O}_{1\times3} & \boldsymbol{O}_{1\times5} & \boldsymbol{O}_{1\times m} & \boldsymbol{O}_{1\times3}\\ \boldsymbol{O}_{3\times1} & \boldsymbol{I}_{3\times3} & \boldsymbol{O}_{3\times m} & \boldsymbol{O}_{3\times m} & \boldsymbol{O}_{3\times3} \\ \boldsymbol{O}_{m\times3} & \boldsymbol{O}_{m\times3} & \boldsymbol{I}_{m\times m} & \boldsymbol{O}_{m\times m} & \boldsymbol{O}_{m\times3} \\ \boldsymbol{O}_{m\times1} & \boldsymbol{O}_{m\times3} & \boldsymbol{O}_{m\times m} & \boldsymbol{I}_{m\times m} & \boldsymbol{F}_{sa~m \times3}^{\rm T} \\ \boldsymbol{O}_{3\times1} & \boldsymbol{O}_{3\times3} & \boldsymbol{O}_{3\times m} & \boldsymbol{F}_{sa~3\times m} & \boldsymbol{I}_{s~3\times3} \end{array} \right],\notag \\ \boldsymbol{A}_{2}(\boldsymbol{x})=\left[\begin{array}{c} -\frac{1}{2}\boldsymbol{q}^{\rm T}\boldsymbol{\omega}_s\\ \frac{1}{2} \left(\boldsymbol{q}^{\times}+q_0\boldsymbol{I}\right)\boldsymbol{\omega}_s\\ \dot{\boldsymbol{\eta}}_{sa}\\ -2\xi_{sa}\boldsymbol{\omega}_{sa}\dot{\boldsymbol{\eta}}_{sa}-\boldsymbol{\omega}^2_{sa}\boldsymbol{\eta}_{sa}\\ -\boldsymbol{\omega}^{\times}_{s}\boldsymbol{I}_s\boldsymbol{\omega}_{s} \end{array} \right], ~~ \boldsymbol{B}=\left[\begin{array}{c} \boldsymbol{O}_{1\times3}\\ \boldsymbol{O}_{3\times3}\\ \boldsymbol{O}_{m\times3}\\ \boldsymbol{O}_{m\times3}\\ \boldsymbol{I}_{1\times3}\end{array} \right]. \notag \end{align} $

It is clear that (10) can be rewritten as the form of (7) and the assumptions a)-f) are satisfied,so characteristic model of the angular velocity can be described in the form of (8) after it is discretized. In order to verify the possible intervals of the characteristic parameters of the CFNSs,simulation is done with several kinds of control signals. Suppose that the control input $u$ has the following four types:

1) Step signal: $u(k)=10$.

2) 0.2 Hz sinusoid signal: $ u(k)=10\sin(0.4k \pi \triangle t)$.

3) 0.2 Hz square wave signal:$ u(k)=10 \textrm{sgn} \left[\sin (0.4k \pi \triangle t)\right].$

4) White noise.

A SSOS,which has one degree freedom single solar panel driven by a constant velocity to point to the sun,is applied to the simulation in open loop. The main parameters of the flexible satellite are presented in Table Ⅰ.

Table Ⅰ
DEFINITION OF THE SSOS

Here we use the recursive least-square method to estimate parameters. Set the forgetting factor $\lambda_s=0.97$. The recursive least-square method is as follows:

$ \begin{align}\label{equa11} \begin{cases} \boldsymbol{K}(k) = \dfrac{\boldsymbol{P}(k-1) \boldsymbol{\phi}(k-1)}{\lambda_s + \boldsymbol{\phi}^{\rm T}(k-1) \boldsymbol{P}(k-1) \boldsymbol{\phi}(k-1)}, \\[3mm] \hat{\boldsymbol{\theta}}(k) = \hat{\boldsymbol{\theta}}(k-1) + \boldsymbol{K}(k) \left( y(k) - \boldsymbol{\phi}^{\rm T}(k-1) \hat{\boldsymbol{\theta}}(k-1) \right) ,\\[2mm] \boldsymbol{P} (k) = \dfrac{1}{\lambda_s} \left( \boldsymbol{I} - \boldsymbol{K}(k) \boldsymbol{\phi}^{\rm T}(k-1) \right) \boldsymbol{P}(k-1) , \end{cases} \end{align} $ (11)

where $\hat{\boldsymbol{\theta}}(k-1) = [\hat{f}_1(k),\ \hat{f}_2(k),\ \hat{g}_0(k),\ \hat{g}_1(k)]^{\rm T}$ and $\boldsymbol{\phi}(k-1) = [y(k-1),\ y(k-2),\ u(k-1),\ u(k-2)]^{\rm T}$.

The simulation results in the steady state are listed in Table Ⅱ, where errors are the difference between the output of the characteristic model (8) and the dynamic model (10).

Table Ⅱ
SIMULATION RESULTS WITH THE INPUTS

The simulation results show it is effective to replace the spacecraft model with the corresponding characteristic model. Moreover,to provide the simplicity of design,the characteristic model is chosen as follows:

$ \begin{align} e_i(k+1)=f_{i1}(k)e_i(k)+f_{i2}(k)e_i(k-1)+g_{i0}(k)u_i(k), \end{align} $ (12)

where $i=1,2,3$ represents the roll,pitch,and yaw axis, respectively,$r_i(k)$ is the desired output,and $e_i(k)=y_i(k)-r_i(k)$.

$ \begin{align}\label{equa13} \boldsymbol{x}_i(k+1)=\left[\begin{array}{cc}0 & 1\\ f_{i2}(k) & f_{i1}(k) \end{array}\right]\boldsymbol{x}_i(k) +\left[\begin{array}{c}0 \\ g_{i0}(k)\end{array}\right]u_i(k). \end{align} $ (13)

Furthermore,with the introduction of new state variables $\boldsymbol{x}_i (k)$ $=$ $[e_i(k-1),e_i(k)]^{\rm T}$,the overall system of the spacecraft can be rewritten as

B. Sliding Mode Control for the Spacecraft

In this paper,a novel sliding mode surface is designed as follows:

$ \begin{align} &s_i (k) = \left[\begin{array}{cc}cs_{i2}(k-1) & cs_{i1}(k-1)\end{array}\right]\boldsymbol{x}_i(k),\notag\\ &cs_{i1}(k-1)=\left(L_1\hat{f}_{i1}(k-1)\right)^\tau,\notag\\ &cs_{i2}(k-1)=\left(L_2\hat{f}_{i2}(k-1)\right)^\tau, \end{align} $ (14)

where $L_1=0.382$ and $L_2=0.618$ are golden section coefficients, $\tau>1$ is the adjustable parameter,and $\hat{f}_{i1}(k)$ and $\hat{f}_{i2}(k)$ are the estimated values of the corresponding coefficients in (12). The coefficients are estimated by the gradient projection algorithm as follows[20]:

$ \begin{align}\label{equa15} \begin{cases} \hat{\boldsymbol{\theta}}_{\rm in} = \hat{\boldsymbol{\theta}}_{i}(k-1)\notag\\[0mm] \quad\quad\ \ + \dfrac{\lambda_{e1}\boldsymbol{\phi}_i(k-1)\left(e_i(k)-\boldsymbol{\phi}_i(k-1)^{\rm T}\hat{\boldsymbol{\theta}}_{i}(k-1)\right)}{\lambda_{e2}+\boldsymbol{\phi}_i(k-1)^{\rm T}\boldsymbol{\phi}_i(k-1)},\\[3mm] \hat{\boldsymbol{\theta}}_{i}(k) = {\rm Pro} \left(\hat{\boldsymbol{\theta}}_{in}(k-1)\right), \end{cases} \end{align} $

where $\hat{\boldsymbol{\theta}}_{i}(k) = [\hat{f}_{i1}(k), \hat{f}_{i2}(k),\hat{g}_{i0}(k)]^{\rm T}$ and $\boldsymbol{\phi}_i(k-1) = [e_i(k-1),e_i(k-2),u_i(k-1)]^{\rm T}$. The positive constants $\lambda_{e1}$ and $\lambda_{e2}$ satisfy $0 < \lambda_{e1} < 1$ and $\lambda_{e2}>0$,respectively. ${\rm Pro}(x)$ represents the orthogonal projection from $x$ to the bounded closed convex set $\Xi$. According to Table Ⅱ,we have

$ \begin{align} \Xi= \begin{cases} 1.9\leq\hat{f}_{i1}\leq2,\ -1\leq\hat{f}_{i2}\leq-0.9,\\ \hat{g}_{i0}\in\left(g_{i0\min},g_{i0\max}\right). \\ \end{cases} \end{align} $ (15)

Usually,the projection operator plays a crucial role in the projection-type adaptive mapping scheme[21]. However it is noticeable that the orthogonal projection algorithm only works when the estimated values leave the given intervals which have been tested.

By (13) and (14),we obtain the equivalent control of DTSMC law as follows:

$ \begin{align} u_{eqi}(k)=-\left[\begin{array}{cc}cs_{i2}(k) & cs_{i1}(k) \end{array}\right]\left[\begin{array}{cc}0 & 1 \\ \hat{f}_{i2}(k) & \hat{f}_{i1}(k)\end{array}\right]\boldsymbol{x}_i(k). \end{align} $ (16)

The variable structure control of DTSMC law is designed as

$ \begin{align} u_{vsi}(k)=s_i(k)-\rho_is_i^{\frac{p_i}{q_i}}-\vartheta_i\textrm{sgn}(s_i(k)), \end{align} $ (17)

where $p_i$ and$q_i$ are both odd positive integers with $p_i < q_i$ and $\rho_i$ and $\vartheta_i$ are positive constants. In order to avoid the singularity,the parameters $p_i$ and $q_i$ should be chosen carefully[22]. $\textrm{sgn}(\cdot)$denotes the sign function,that is

$ \begin{align}\notag \textrm{sgn}(x)= \begin{cases}1,& x>0,\\0,& x=0,\\-1,& x < 0.\end{cases} \end{align} $

The control input for the spacecraft is proposed as

$ \begin{align} u_i(k)=\frac{u_{eqi}(k)+u_{vsi}(k)}{cs_{i1}(k)\hat{g}_{i0}(k)}. \end{align} $ (18)

Then we will analyze the behavior of sliding mode surface and the output tracking precision using the characteristic model-based DTSMC law. Before moving on,we give the following lemma and definition.

Definition 1[23]. The Lypaunov function is chosen to be $V(k)$ $=$ $s^2(k)/2$. The sliding manifold $s(k)=0$ is attractive if $V(k$ $+$ $1)-V(k) < 0$. This condition is satisfied if

$ \begin{align} \left(s(k+1)-s(k)\right)\textrm{sgn}((s(k)) < 0,\notag\\ \left(s(k+1)+s(k)\right)\textrm{sgn}((s(k))>0. \end{align} $ (19)

Lemma 2[24]. Consider the scalar dynamical system

$ \begin{align} z(k+1) = z(k) -l z(k)^\alpha +j(k), \end{align} $ (20)

where $l>0$ and $0 < \alpha < 1$ is a ratio of odd integers. If $| j(k)|$ $\leq$ $\gamma$,$\gamma>0$,there is a finite number $K^*>0$ such that

$ \begin{align} | z(k)| \leq \psi(\alpha)\cdot \max\left\{\left(\frac{\gamma}{l}\right)^{\frac{1}{\alpha}},l^{\frac{1}{1-\alpha}}\right\}, \end{align} $ (21)

where function $\psi(\alpha)$ is defined as

$ \begin{align}\notag \psi(\alpha)=1+\alpha^{\frac{\alpha}{1-\alpha}}-\alpha^{\frac{1}{1-\alpha}}. \end{align} $

Based on Definition 1 and Lemma 2,we have the following theorem.

Theorem 1. If the error characteristic model of system (10) can be expressed with (12),the sliding mode surface is defined as (14),and the control law is designed as (19),then

a) The overall system (13) is stable with

$ \begin{align} |s_i(k)|>\frac{\rho_i|s_i(k)|^{\frac{p_i}{q_i}}+\vartheta_i}{2}. \end{align} $ (22)

b) There is a finite number $K^*_1>0$ such that

$ \begin{align} | s_i(k)| \leq \psi\left(\frac{p_i}{q_i}\right)\cdot \max \left\{\left(\frac{\vartheta_i}{\rho_i}\right)^{\frac{q_i}{p_i}},\rho_i^{\frac{q_i}{q_i-p_i}}\right\}, \ \forall k\geq K^*_1. \end{align} $ (24)

c) There is a finite number $K^*_2>0$ such that

$ \begin{align} | e_i(k) | \leq \frac{smax+| \left(L_2f_{i2}(\infty)\right)^\tau|\cdot| e_i(0)|}{| \left(L_1f_{i1}(\infty)\right)^\tau |},\ \forall k\geq K^*_2, \end{align} $ (24)

where

$ \begin{align}\notag smax=\psi\left(\frac{p_i}{q_i}\right)\cdot \max \left\{\left(\frac{\vartheta_i}{\rho_i}\right)^{\frac{q_i}{p_i}},\rho_i^{\frac{q_i}{q_i-p_i}}\right\}, \end{align} $

$f_{i1}(\infty)$ and $f_{i2}(\infty)$ can be found in Table Ⅱ, and $e_i(0)$ is the initial error.

Proof. Choosing Lypapunov function $V_i(k)=s_i^2(k)/2$,the behavior of $s_i(k)$ can be prescribed as

$ \begin{align} s_i(k+1) = \left[\begin{array}{cc}cs_{i2}(k)&cs_{i1}(k)\end{array}\right]\boldsymbol{x}_i(k+1). \end{align} $ (25)

Substituting (13) and (19) into (26) yields

$ \begin{align} s_i(k+1) = s_i(k)-\rho_is^{\frac{p_i}{q_i}}_i(k)-\vartheta_i\textrm{sgn}(s_i(k)), \end{align} $ (26)

and then

$ \begin{align} \left(s_i(k+1)-s_i(k)\right)\textrm{sgn}(s_i(k))=-\rho_i | s_i(k) | ^{\frac{p_i}{q_i}},\notag \\ -\vartheta_i+| s_i(k) | < 0,\notag\\[2mm] \left(s_i(k+1)+s_i(k)\right)\textrm{sgn}(s_i(k))= -\rho_i | s_i(k) | ^{\frac{p_i}{q_i}},\notag \\ -\vartheta_i+2| s_i(k) |>0. \end{align} $ (27)

Hence the control law (19) makes (13) move to the predefined sliding mode surface. Observing $s_i(k)=0$ in sliding phase,it follows that

$ \begin{align} e_i(k)&=-\left(\frac{cs_{i2}(k-1)}{cs_{i1}(k-1)}\right)^\tau\notag\\ &= -\left(\frac{L_2\hat{f}_{i2}(k-1)}{L_1\hat{f}_{i1}(k-1)}\right)^\tau e_i(k-1). \end{align} $ (28)

Substituting (16) into (29),we have $0.7258\leq L_1\hat{f}_{i1}\leq 0.764$ and $-0.618\leq L_2\hat{f}_{i2}\leq -0.5562$,which results in $| (\frac{cs_{i2}(k-1)}{cs_{i1}(k-1)})^\tau |$ $ < $ $1$. Thus the overall system is stable.

Then it follows from (27),considering $|-\vartheta_i\textrm{sgn}(s_i(k))|=\vartheta_i$,based on Lemma 2, there exists a finite number $K^*_1>0$ such that when $\forall k\geq K^*_1$ there is

$ \begin{align}\notag | s_i(k)| \leq \psi\left(\frac{p_i}{q_i}\right)\cdot \max \left\{\left(\frac{\vartheta_i}{\rho_i}\right)^{\frac{q_i}{p_i}},\rho_i^{\frac{q_i}{q_i-p_i}}\right\}. \end{align} $

Substituting (19) into (12) yields

$ \begin{align} &e_i(k+1)=\notag\\ &\ \ \ \frac{\left(s_i(k)-\rho_is_i^{\frac{p_i}{q_i}}(k)-\vartheta_i\textrm{sgn}(s_i(k)) -\left(L_2\hat{f}_{i2}(k)\right)^\tau e_i(k)\right)}{\left(L_1\hat{f}_{i1}(k)\right)^\tau}. \end{align} $ (29)

Noticing that $| (L_2\hat{f}_{i2}(k))^\tau | < 1$ which results in $| (L_2\hat{f}_{i2}(k))^\tau |$ $\cdot$ $| e_i(k)|$ $ < $ $| (L_2f_{i2}{(\infty)})^\tau| \cdot| e_i(0)|$,it follows from Lemma 2 that there is a finite number $K^*_2>0$ such that when $\forall k\geq K^*_2$ there is

$ \begin{array}{l} |{e_i}(k)|\\ \le \frac{{\left( {\psi (\frac{{{p_i}}}{{{q_i}}}) \cdot \max \left\{ {{{(\frac{{{\vartheta _i}}}{{{\rho _i}}})}^{\frac{{{q_i}}}{{{p_i}}}}},\rho _i^{\frac{{{q_i}}}{{{q_i} - {p_i}}}}} \right\} + |{{\left( {{L_2}{f_{i2}}(\infty )} \right)}^\tau }| \cdot |{e_i}(0)|} \right)}}{{\left| {{{\left( {{L_1}{{\hat f}_{i1}}(k)} \right)}^\tau }} \right|}}\\ \le \frac{{\left( {\psi (\frac{{{p_i}}}{{{q_i}}}) \cdot \max \left\{ {{{(\frac{{{\vartheta _i}}}{{{\rho _i}}})}^{\frac{{{q_i}}}{{{p_i}}}}},\rho _i^{\frac{{{q_i}}}{{{q_i} - {p_i}}}}} \right\} + |{{\left( {{L_2}{f_{i2}}(\infty )} \right)}^\tau }||{e_i}(0)|} \right)}}{{|{{\left( {{L_1}{f_{i1}}(\infty )} \right)}^\tau }|}}. \end{array} $ (30)

Remark 1. According to (29),when the error reaches its sliding surface,it will converge to zero with a speed determined by chosen $\tau$,which lies in $\left[0.728^\tau,0.85^\tau \right]$. Usually,we choose $2\leq \tau \leq 10$.

Remark 2. Actually,in practice,according to Table Ⅱ and (16),we can obtain $f_{i1}(\infty)\approx\lim\nolimits_{k \rightarrow K^*_2}\hat{f}_{i1}(k)\approx2$ and $f_{i2}(\infty)$ $\approx$ $\lim\nolimits_{k \rightarrow K^*_2}\hat{f}_{i2}(k)\approx-1$,which lead to

$ \begin{align} | e_i(k) | \leq\frac{smax+ 0.618^\tau | e_i(0) |}{0.764^\tau},\ \ \ \forall k\geq K^*_2. \end{align} $ (31)

Theorem 1 demonstrates that for the closed system (10) with the characteristic model (12),under the DTSMC law (19),the sliding mode surface and the error will converge to its relevant bound in finite time and the predefined bound mainly depends on the adjustable parameters of the DTSMC law.

Ⅳ. NUMERICAL SIMULATIONS

Suppose that the SSOS is required to rapidly maneuver considering the tilt of the cantilevered appendage. Simulations are conducted for demonstrating the performance of the proposed DTSMC law. The closed system logic structure is shown in Fig.1.

Download:
Fig. 1. Control logic block diagram for the spacecraft.

Four different sinusoidal-wave disturbances as in (33) are added to the spacecraft randomly. The external disturbances used in the simulation are far worse than those typically observed in practice[25].

$ \begin{align}\label{equa33} d_1(t)=[&0.10\sin(0.4t),0.05\cos(0.5t),\notag\\ &0.08\cos(0.7t)]^{\rm T} {\rm (N\cdot m)},\notag\\[1mm] d_2(t)=[&0.06\cos(0.4t),\notag\\ &0.10\sin(0.5t),0.05\sin(0.7t)]^{\rm T} {\rm (N\cdot m)},\notag\\[1mm] d_3(t)=[&0.08\sin(0.4t+\frac{\pi}{4}),\ \notag\\ &0.06\cos(0.5t+\frac{\pi}{4}),\ 0.07\cos(0.7t+\frac{\pi}{4})]^{\rm T} {\rm (N\cdot m)},\notag\\[1mm] d_4(t)=[&0.08\cos(0.4t+\frac{\pi}{4}),\ \notag\\ &0.08\sin(0.5t+\frac{\pi}{4}),\ 0.10\sin(0.7t+\frac{\pi}{4})]^{\rm T} {\rm (N\cdot m)}. \end{align} $ (32)

The desired attitude velocity command is described by using the bang-coast-bang method and is chosen as

$ \begin{align} \label{equa34} r_i(t)= \begin{cases} 0.115t,& 0 < t\leq10 {\rm (^\circ/s)},\\ 1.15,& 10 < t\leq30 {\rm (^\circ/s)},\\ 1.15-0.015t,& 30 < t\leq40 {\rm (^\circ/s)} ,\\ 0,& t>40 {\rm (^\circ/s)}. \end{cases} \end{align} $ (33)

For the sliding mode,the parameter is $m=2$. The parameters of the controller are tuned as $p_i=1$,$q_i=5$,$\rho_i=1$,and $\vartheta_i =1\times 10^{-4}$. The simulation results are shown in Figs. 2-7. The simulation results reveal that converging times of three-axis are very close. Maneuvering angles are $43^\circ$, $20^\circ$,and $43^\circ$,respectively. The rapid maneuvering performance and point accuracy shown in Figs. 3 and 4 are significantly improved compared to the robust fuzzy $H_{\infty}$ controller[1]. Table III shows the simulation results. The interference of the flexible appendage vibration is suppressed. It follows from Figs. 5 and 6 that tracking effects are achieved in finite time. The control input is shown in Fig.7. It illustrates that chattering is avoided because of the bounded DTSMC implementation.

Download:
Fig. 2. Quaternion responses.

Download:
Fig. 3. Attitude angles.

Download:
Fig. 4. Angular rates.

Download:
Fig. 5. Tracking errors.

Download:
Fig. 6. Partial enlargement of tracking errors.

Download:
Fig. 7. Control inputs.

Table Ⅲ
THE RESULTS OF COMPARISON WITH [1]
 

As can be seen from the phase plane in Figs. 8-10,the characteristic model-based sliding mode surface and the controller demonstrate a fast convergence and strong robustness.

Download:
Fig. 8. Phase plane of the first channel.

Download:
Fig. 9. Phase plane of the second channel.

Download:
Fig. 10. Phase plane of the third channel.

Then the rotational inertia of satellite is considered as $1.5{\pmb I}_s$ and $0.5{\pmb I}_s$,respectively. The simulation results in Table Ⅳ illuminate the proposed controller can provide a satisfactory tracking result.

Table Ⅳ
SIMULATION RESULTS
Ⅴ. CONCLUSION

The requirement of attitude tracking in AOCSs is a challenging and important objective for many practical spacecraft missions. In this paper,a finite-time DTSMC law has been proposed to achieve this objective for the spacecraft with variable tilt of flexible structures using the characteristic modeling method. The designed methodology has used the characteristic model,which includes all the information of the spacecraft,instead of the dynamic model to avoid the enormously complex cantilever hybrid model. Based on such characteristic model, an adaptive sliding mode surface,which was composed by the parameters of the characteristic model and the relative angular velocity errors of the spacecraft,has been designed to construct the finite-time DTSMC law. Numerical simulations of a SSOS has been performed to validate the effectiveness of the proposed control law in the presence of severe model uncertainties and disturbances. Simulation results of the control performance have demonstrated that the spacecraft converges to the desired attitude and angular velocity in a short time. Furthermore,the vibration caused by the maneuver of the appendage has been suppressed as well as the chattering phenomenon caused by the switching control of the SMC has been reduced.

References
[1] Chen J X, Sun F C, Sun Y G, Yu L Y. Modeling and controller design for complex flexible nonlinear systems via a fuzzy singularly perturbed approach. Information Sciences, 2014, 255:187-203
[2] Sun Chang-Yin, Mu Chao-Xu, Yu Yao. Some control problems for near space hypersonic vehicles. Acta Automatica Sinica, 2013, 39(11):1901-1913(in Chinese)
[3] Wang W, Menon P P, Bates D G, Bennani S. Robustness analysis of attitude and orbit control systems for flexible satellites. IET Control Theory and Applications, 2010, 4(12):2958-2970
[4] Xia Y Q, Fu M Y. Compound Control Methodology for Flight Vehicles. Berlin:Springer, 2013. 55-79
[5] Gao H J, Yang X B, Shi P. Multi-objective robust H control of spacecraft rendezvous. IEEE Transactions on Control Systems Technology, 2009, 17(4):794-802
[6] Lee T. Exponential stability of an attitude tracking control system on SO(3) for large-angle rotational maneuver. System and Control Letters, 2012, 61(1):231-237
[7] Sun C Y, Mu C X, Zhang R M. Terminal Sliding Mode Control Method for Hypersonic Vehicles. Beijing:Science Press, 2014. 156-182
[8] Guy N, Alazard D, Cumer C, Charbonnel C. Dynamic modeling and analysis of spacecraft with variable tilt of flexible appendages. Journal of Dynamic Systems, Measurement, and Control, 2014, 136(2):021020
[9] Wu H X, Liu Y W, Liu Z H, Xie Y C. Characteristic modeling and the control of flexible structure. Science in China Series:Information Sciences, 2001, 44(4):278-291
[10] Valentin-Charbonnel C, Duc G, Le Ballois S. Low-order robust attitude control of an earth observation satellite. Control Engineering Practice, 1999, 7(4):493-506
[11] Meng B, Wu H X, Lin Z L, Li G. Characteristic model based control of the X-34 reusable launch vehicle in its climbing phase. Science in China Series F:Information Sciences, 2009, 52(11):2216-2225
[12] Yang K D, Utkin V I, Ozguner U. A control engineer's guide to sliding mode control. IEEE Transactions on Control Systems Technology, 1999, 7(3):328-342
[13] Yu X H, Wang B, Li X J. Computer-controlled variable structure systems:the state-of-the-art. IEEE Transactions on Industrial Informatics, 2012, 8(2):197-205
[14] Song R Z, Xiao W D, Sun C Y. Optimal tracking control for a class of unknown discrete-time systems with actuator saturation via data-based ADP algorithm. Acta Automatica Sinica, 2013, 39(9):1413-1420
[15] Galias Z, Yu X H. Study of periodic solutions in discretized twodimensional sliding-mode control systems. IEEE Transactions on Circuits and Systems II:Express Briefs, 2011, 58(6):381-385
[16] Galias Z, Yu X H. Dynamical behaviors of discretized second-order terminal sliding-mode control systems. IEEE Transactions on Circuits and Systems II:Express Briefs, 2012, 59(9):597-601
[17] Chen L, Yan Y, Sun C Y. Nonlinear characteristic model-based SMC and its application to flexible. In:Proceedings of the 19th IFAC World Congress. Cape Town, South Africa:IFAC, 2014. 4595-4600
[18] Yang F S, Zhang H G, Wang Y C. An enhanced input-delay approach to sampled-data stabilization of T-S fuzzy systems via mixed convex combination. Nonlinear Dynamics, 2014, 75(3):501-512
[19] Wu H X, Hu J, Xie Y C. Characteristic model-based all-coefficient adaptive control method and its applications. IEEE Transactions on Systems, Man, and Cybernetics, Part C:Applications and Reviews, 2007, 37(2):213-221
[20] Wu H X, Hu J, Xie Y C. Characteristic Model-Based Intelligent Adaptive Control. Beijing:China Science and Technology Press, 2009. 113-158
[21] Yang F S, Zhang H G, Jiang B, Liu X D. Adaptive reconfigurable control of systems with time-varying delay against unknown actuator faults. International Journal of Adaptive Control and Signal Processing, 2014, 28(11):1206-1226
[22] Mu Chao-Xu, Yu Xing-Huo, Sun Chang-Yin. Phase trajectory and transient analysis for nonsingular terminal sliding mode control systems. Acta Automatica Sinica, 2013, 39(6):902-908(in Chinese)
[23] Sarpturk S Z, Istefanopulos Y, Kaynak O. On the stability of discretetime sliding mode control systems. IEEE Transactions on Automatic Control, 1987, 32(10):930-932
[24] Li S H, Du H B, Yu X H. Discrete-time terminal sliding mode control systems based on Euler's discretization. IEEE Transactions on Automatic Control, 2014, 59(2):546-552
[25] Wu B L, Wang D W, Poh E K. Decentralized sliding-mode control for attitude synchronization in spacecraft formation. International Journal of Robust and Nonlinear Control, 2013, 23(1):1183-1197