2. Department of Electrical Engineering and Automation, Tianjin University, Tianjin 300072, China
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
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
Sliding mode control (SMC)
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 CRITERIAThe 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
$ \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) |
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
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
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.
$ \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) |
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
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
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
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
$ \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
$ \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
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
$ \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.
$ \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) |
$ \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.
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
$ \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) |
$ \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 SIMULATIONSSuppose 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
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
$ \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
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. |
As can be seen from the phase plane in
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
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.
[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 |