In this article, a three-axis attitude manoeuvre spacecraft consisting of a central rigid body and a rotating solar array is studied. The rotating solar array is considered a disturbance to the spacecraft. In the design of the controller, the coupled terms and the rotating solar array are considered a disturbance. The improved extended state observer is proposed by combing the sliding mode observer with the originally extended state observer to estimate the disturbance. The sliding mode control method is adopted to adjust the attitude of the spacecraft. Numerical simulations are presented to demonstrate the outstanding performance of the present observer.

#### Keywords

- rigid spacecraft
- rotating solar array
- disturbance
- improved extended state observer
- sliding mode control

#### MSC 2010

- ???

Solar arrays on the spacecraft track and face the Sun to provide the power to the payload. The rotation of the solar array has effects on the attitude of the spacecraft, which causes a series of problems about multi-body dynamics and the accuracy of the attitude control.

There exists the relative movement between the central platform and the solar array due to the rotation of the solar array. The multi-body dynamics of such structure have attracted extensive attention of scholars. Gasbarri et al. [1] devoted themselves to study the structural uncertainty of a very large space structure. The model was established by the finite element method (FEM), and a nonlinear control scheme was adopted to solve the robustness of the system. In their another paper, the model was reconstructing and a new control concept called ‘control-oriented modelization’ was proposed to solve the problem of the model's uncertainties and disturbances [2]. Sattar and Wei [3] established an accurate mathematical model for kinematic and coupled micro vibrations of a rigid aluminium plate driven with external excitation torque. Iwata [4] proposed a composite control scheme for compensating and rejecting the disturbance of the Advanced Land Observing Satellite (ALOS). Based on the adaptive-robust control law, Azadi et al. [5] used piezoelectric layers as actuators to suppress the residual vibration of the flexible appendages during the spacecraft attitude manoeuvring in space. Lu and Liu [6] designed the _{∞} optimal control to improve the robustness of the spacecraft with the rotating solar array, where the singularity of the system was also considered. Besides the above, the sliding mode control (SMC) method is widely used as a robust approach to solve the disturbance and uncertainty. Hu et al. [7] combined the variable structure output feedback and input shaping and considered adaptive law to design control law. The control scheme not only guaranteed the stability of the closed-loop system but also overcame the parameter uncertainty and external disturbances. Hu [8] studied the SMC scheme for the three-axis stabilised flexible spacecraft. Cao et al. [9] combined SMC with input shaping techniques to reduce the vibration of the driving solar arrays. The controller accomplished asymptotic attitude manoeuvring in the closed-loop system and performed good robustness. For the rigid spacecraft with the disturbance, the SMC with the adaptive algorithm was adopted to track the attitude and compensate the disturbance [10].

In this article, the flexible vibration of the appendage is not considered. Instead, the multi-rigid-body motion of the spacecraft with disturbance is studied. To estimate and compensate the uncertainty and disturbance, more and more researchers pay attention to the new control algorithm that involves the application of the state observer. The disturbance observer-based control (DOBC) is a promising approach, which can be used to eliminate the disturbance and improve the robustness of the system. For the spacecraft with an attitude adjustment and Sun-tracking missions, a feed-forward control based on a disturbance observer was designed [11]. Chen and Kun-Yung [12] proposed the adaptive sliding mode tracking scheme to control the two-link manipulator robot system, where the disturbance observer is used to stabilise and compensate the unknown time-varying disturbances. Lee [13] proposed a composite control technique with a nonlinear disturbance observer for the attitude tracking of the rigid spacecraft, where the uncertain inertia parameters were compensated through feedforward and the robustness performance of the spacecraft was enhanced. In addition to the above, there are two other state observers widely used, namely, sliding mode observer (SMO) [14,15,16,17] and extended state observer (ESO) [18,19,20,21].

The SMOs are widely used for finite-time convergence, robustness with respect to uncertainties and the possibility of uncertainty estimation [15]. The nonlinear sliding mode observers are adopted on the robot manipulators estimation [14]. An important algorithm for control and observation using SMO is the so-called super twisting algorithm described by the differential inclusion [16]. The second-order SMC was used to track the constant reference and had been validated on the unstable non-minimum phase inverted pendulum system [22]. Moreno and Osorio [17] gave the rigorous and detailed proof using the Lyapunov approach for the convergence of the SMO. The synchronous motor was controlled by the adaptive sliding mode method, where the estimation of the rotor speed was obtained by the SMO [23].

The ESO was first proposed together with the concept of the active disturbance rejection controller by Professor J.Q. Han. It was an effective way to estimate the so-called uncertainty from the unknown model and perturbation [21]. The key part of the ESO is that the uncertainty and disturbance of the system are considered as a new order state. Then, the estimations of all states including the added new one can be obtained by a specific nonlinear function, where the gains of the function need to be selected appropriately. The bandwidth and stability of the linear ESO (LESO) are determined by the gains, which can be ensured by allocating the poles rationally [24, 25]. The impact of the observed bandwidth changes on the performance of the LESO was discussed in detail by Godbole et al. [26]. However, the LESO was out of its depth in dealing with the nonlinear system. Then, the ESO with the nonlinear function must be adopted. The stability and convergence of nonlinear ESO were given by using the piece-wise smooth Lyapunov function [18]. Huang and Han [19] and Huang [20] proposed the self-stable-region (SSR) theory. This theory could improve the ESO convergence performance and observation accuracy. Wang et al. [27] added an output signal filter to the ESO to eliminate measure noise. Then, the ADRC and feedback compensation were presented together in dealing with the uncertainty of the lateral thrust and aerodynamics blended control system [28]. Li et al. [29] adopted the ESO to estimate all states of the attitude stabilisation system for the spacecraft. Besides, the uniformly ultimately bounded stability of the ESO was assured by rigorous Lyapunov analysis. For the attitude tracking of the rigid spacecraft, the SMC scheme was designed to force the state variables to converge to the reference states [30]. In this case, the ESO was indispensable to estimate the disturbance. Zhu et al. [31] used the ESO to estimate the target acceleration and achieve the missile interception. The chattering was proved to be decreased by combing the SMC with the ESO. The third-order discrete ESO was employed on the antenna pointing system [32]. It made each channel decouple and linearise successfully.

In this article, a three-axis attitude motion of the spacecraft with a single rotating solar array is studied. The SMC method is designed to achieve the attitude adjustment. The improved ESO (IESO) is proposed to estimate the disturbance and uncertainty, where the sliding mode observer is adopted to replace the nonlinear function

As shown in Figure 1, the satellite consists of a central rigid platform and one solar array. The attitude manoeuvring of the spacecraft is in a small range and the solar array is driven in a fairly low velocity, which means these two motions do not excite the elastic deformation of the solar array. Thus, the whole spacecraft is considered as rigid. The reference coordinate system _{0}_{0}_{0} is set on the core of the earth. The body coordinate system _{1}_{1}_{1} is fixed on the satellite. The origin _{2}_{2}_{2} is embedded on the solar array that the origin

Let _{C}_{x}_{y}_{z}^{T} denotes the vector of the angular velocities of the spacecraft and _{A}_{A}^{T} denotes the vector of the angular velocities of the solar array. The dynamic equations of the spacecraft with rotating solar array are expressed as follows [6, 10, 13, 30]:
_{t} = _{C} + _{CA}_{b}_{AC} is the rotary inertia of the whole spacecraft with respect to the body-fixed coordinate system, and is the rotary inertia of the solar arrays with respect to the rotation coordinate system.
_{C}_{A}_{d}_{CA} is the coordinate transformation matrix from the rotation coordinate frame to the body-fixed coordinate frame. It is expressed as follows:

Euler angle is used to denote the attitude of the satellite usually. Three coordinate rotations in sequence can describe any rotation. The satellite rotates in a sequence of _{0} is the orbital angular velocity. In this article, the spacecraft attitude manoeuvre is in a small range, the variations of the attitude angles are small enough that the linearisation form is as follows:

Thus for simplicity, Eq. (5) can be rewritten in the short form:
^{T}, _{0}_{0}, _{0}^{T} and _{x}_{y}_{z}^{T}.

The control system for the spacecraft with the rotating solar array consists of double closed loops. As shown in Figure 2, the outer loop is for the attitude motion system with an angular position regulator while the inner loop is the attitude dynamic system. In Figure 2, the attitude dynamic system employs an ESO to estimate the state variable. The sliding mode controller is adopted to adjust the attitude. The attitude angular _{d}

If we set _{x}_{y}_{z}_{A}_{A}^{T}, _{A}_{A}^{T}, the state space form of Eq. (1) can be written as
_{1}, τ_{2}, τ_{3}, τ_{A}]^{T}.

The symbols in the above matrices are expressed as

Some terms in Eq. (7) are the products of two or more angle velocity vectors, so they can be ignored in the calculation as higher-order small quantities. Then the simplified equation is obtained as follows:

According to Eq. (5), three channels are decoupled and the rotation of the solar array has effects on the pitching axis directly. Taking the pitching channel as an example, the controller is designed in this part. The state variables are defined as follows:
_{1} is the state of the system while _{2} is the extended state.

Based on the concept of ESO [21], second-order ESO is defined as follows:
_{11} and _{12} are the tracking state of _{1} and _{2}, respectively. _{11} is the output error between the observed value and state output of the inner-loop. The double number of the subscript is defined to distinguish from the ESO in the solar array driving system, which will be mentioned later. The nonlinear function

To improve the performance of the observer, this article proposed the IESO. It is combined with the SMO as follows:
_{i}_{11} of IESO is the estimation of the state of the dynamic system. _{12} uses a similar form of the SMO to estimate the disturbance instead of the higher-order state.

The error dynamic equations are obtained by Eqs (10) and (13):
_{12} = _{12} − _{2}.

The stability analysis in the same form as Eq. (14) has been discussed in detail by Moreno and Osorio [16, 17]; thus some useful conclusions are cited directly in this article.

[16]. _{1} ≥ 0 _{2} ≥ 0. _{i} can be selected large enough so that the system is strongly globally asymptotically stable, where the origin is an equilibrium point. Then k_{i} can be given by

Where

The controller and observer for the other two channels can be designed according to the above-mentioned steps.

As shown in Figure 3, if the controller and the ESO operate well, then the inner loop can be considered as a pathway. The attitude motion system is presented as an error feedback control system.

For the attitude motion formulation in Eq. (6), the ^{*} of this system; the attitude motion is rewritten as follows:
^{*} is designed as
_{0} = _{0},_{0},_{0}], _{1}(•), _{2}(_{3}(^{T}.

Thus, the nonlinear angle position error feedback law is designed as

Without loss of generality, taking the pitching channel as an example to discuss the stability of the regulator, the input error for the outer-loop system is

As the desired attitude angle velocity

To show the stability of the ESO, Eq. (21) can be linearised as follows:
_{0} ≠ 0, setting the Lyapunov function as
_{1} > 0 by first mean value theorem for integration. Differentiating _{1} with respect to time and considering Eq. (24), one obtains

Only when _{0} = 0, one has

As the solar array is restricted to rotating around the _{A}_{A}_{B}_{A}_{B}_{m}_{0} is the viscous friction coefficient. _{r}

It can be seen that there exists the attitude motion coupled term

In Figure 4, the second-order ESO is employed on the feedback path. The estimation of the rotating velocity _{21} is transformed into a motor current signal _{real}_{real}_{ref}_{m}_{d}_{A}_{22} is obtained as follows:

In this article, the attitude manoeuvre is achieved by the sliding mode method. With the help of the IESO, the perturbation can be eliminated as much as possible. The sliding surface for the three-axis manoeuvre is selected as follows:
_{1},_{2},_{3}]^{T}. _{1},_{2},_{3}] is a designed parameter to stabilise the system. should meet the Hurwitz conditions that the system is stable as long as the trajectory of the system reaches the sliding surface in finite time. _{1} and _{2} are defined as
_{d}_{1} and _{2} are the vectors of the estimations obtained by the ESO.

The exponential rate reaching law is applied as follows:

Combing with Eq. (31), the derivative of _{i} = _{i}_{i}_{i}

The control law is obtained by Eqs (32) and (33) as follows:

To discuss the stability of the control law, a candidate Lyapunov function with the linear sliding surface is given as

Taking the derivatives of Eq. (35) and considering the Eqs (33) and (34) yield

It is obvious that
_{i}_{i}_{i}

The numerical simulation of the proposed control scheme to the rigid spacecraft is given by the MATLAB software. The rotary inertia of the spacecraft consists of _{C} and _{b} as follows:

The parameters of the control scheme for the spacecraft and the solar array system driving system are summarised in Tables 1 and 2, respectively.

Control parameters selected for numerical analysis

SMC | |||

_{1} = | _{s} | ||

_{2}= | |||

_{1} = _{1}, _{2}= _{2}, | Error regulator | _{0} = | |

_{3} = _{4} = |

Solar array driving system selected of numerical analysis

_{m} | 0.3 | Viscous friction coefficient _{0} | 26 |

_{r} | 300 | Proportional coefficient _{p} | 10 |

0.2 | Integration coefficient _{i} | 300 | |

68 | Desired angular velocity _{d}^{°}/ | 0.06 |

In the beginning of the manoeuvre, the attitude angle is assumed to be [_{0}, _{0}, _{0}]^{T} = [5, −2.1]^{T} and the angle position of the solar array is zero relative to the central rigid body.

For the analysis of the validity and advantages of the IESO, Figures 5 and 6 show the observer errors obtained from Eqs (11) and (13) under the SMC law. _{1i}(_{2i} represents the estimated errors of the disturbance, respectively. As shown in Figures 5 and 6, both the cases estimate the states successfully. The error of each observer reaches a stable stage of <3 s. The right columns of Figures 5 and 6 are the partially enlarged image of the observer errors at the stable stage. They demonstrate the high precision of the ESO was guaranteed. Obtained from the observed errors, it is obvious that the ones which adopted the improved ESO converge to zero have little fluctuation. The ranges of the errors obtained by the original ESOs are, and |_{13}| < 1, respectively. However, the ones of the IESOs are |_{11}| < 0.01, |_{11}| < 0.02 and |_{13}| < 0.01, respectively. The reduction of the amplitude is an order of magnitude. From Figure 6, even though the error ranges for the disturbance observations of two kinds of the ESO are roughly similar, the fluctuation frequency of the IESO is much lower than the original ESO. The relative stable estimations are always easy to be observed and make use of in engineering.

Figures 7 and 8 show the inertial angular velocities and the Euler attitude angles, respectively. The right columns of these two figures are the partially enlarged drawing of the corresponding steady stage. From the left columns of Figures 7 and 8, both the inertial angular velocity and the Euler attitude angle under two kinds of ESO reach the commands <5 s. However, the oscillation performance is much different from these two observers. From the partially enlarged drawing in Figure 7, the amplitude with original ESO is nearly 10 times more than that under improved ESO. Correspondingly, the curves of three attitude angles have an analogous trend. The amplitudes of the Euler angles obtained from IESO are < 2 × 10^{−4}, which is much lower than the reference ones (obtained from ESO) about 1 × 10^{−2}. Even though the steady-state values are not zero, the precision is sufficient to meet the engineering requirement.

Figure 9 shows the rotation of the solar array under the driving torque, where the two control schemes are considered. Figure 9(a) is the variation of the angle position while Figure 9(b) is the angle velocity. At the initial stage, the angle is negative, which represents the solar array rotates in the opposite direction to the objective. This phenomenon depends on the coordinate frame selected for the observation. The attitude driving torque is large enough, and there exists a lag in the rotation of the solar array while offsetting its driving torque. However, with the attitude motion of the spacecraft turning to stable, the solar array driving performs well and the velocity tends to steady. Similar to the variation of attitude angles, the angle velocity of the solar array obtained by the IESO is more stable than that of the ESO. It demonstrates the better performance of the IESO over the original ESO from the perspective of solar array.

In this article, a three-axis attitude manoeuvre spacecraft with a rotating solar array is studied from the attitude dynamic point of view. Considering the rigid solar array, the multi-body dynamic model is obtained.

To decouple the system, the coupled terms of the dynamical equation are considered the disturbance. The SMC scheme with second-order ESO is adopted to adjust the attitude and reduce the disturbance of the spacecraft. The scheme helps the spacecraft to compensate the disturbance effectively and reach the target rapidly. From the simulation results, the IESO presents good performance to reduce the disturbance and achieve the attitude manoeuvre. The IESO has higher accuracy than the original ESO on the premise of ensuring the convergence time. Moreover, the solar array driving system is also benefitted from the IESO. The fluctuation of the driving is so little that the attitude adjustment can be smooth. In addition, the good robustness of the IESO is also presented by the results that the inertia moments variation within a certain range.

However, it should be mentioned that the solar array driving is simplified by neglect some electromagnetic effects. If they have been fully considered, the complex driving scheme may be a subject for future researches.

#### Control parameters selected for numerical analysis

SMC | |||

_{1} = | _{s} | ||

_{2}= | |||

_{1} = _{1}, _{2}= _{2}, | Error regulator | _{0} = | |

_{3} = _{4} = |

#### Solar array driving system selected of numerical analysis

_{m} | 0.3 | Viscous friction coefficient _{0} | 26 |

_{r} | 300 | Proportional coefficient _{p} | 10 |

0.2 | Integration coefficient _{i} | 300 | |

68 | Desired angular velocity _{d}^{°}/ | 0.06 |

Regarding new wave distributions of the non-linear integro-partial Ito differential and fifth-order integrable equations Nonlinear Mathematical Modelling of Bone Damage and Remodelling Behaviour in Human Femur Value Creation of Real Estate Company Spin-off Property Service Company Listing Entrepreneur's Passion and Entrepreneurial Opportunity Identification: A Moderated Mediation Effect Model Applications of the extended rational sine-cosine and sinh-cosh techniques to some nonlinear complex models arising in mathematical physics Study on the Classification of Forestry Infrastructure from the Perspective of Supply Based on the Classical Quartering Method A Modified Iterative Method for Solving Nonlinear Functional Equation New Principles of Non-Linear Integral Inequalities on Time Scales Has the belt and road initiative boosted the resident consumption in cities along the domestic route? – evidence from credit card consumption Analysis of the agglomeration of Chinese manufacturing industries and its effect on economic growth in different regions after entering the new normal Study on the social impact Assessment of Primary Land Development: Empirical Analysis of Public Opinion Survey on New Town Development in Pinggu District of Beijing Possible Relations between Brightest Central Galaxies and Their Host Galaxies Clusters and Groups Attitude control for the rigid spacecraft with the improved extended state observer An empirical investigation of physical literacy-based adolescent health promotion MHD 3-dimensional nanofluid flow induced by a power-law stretching sheet with thermal radiation, heat and mass fluxes The research of power allocation algorithm with lower computational complexity for non-orthogonal multiple access Research on the normalisation method of logging curves: taking XJ Oilfield as an example A Method of Directly Defining the inverse Mapping for a HIV infection of CD4+ T-cells model On the interaction of species capable of explosive growth Research on Evaluation of Intercultural Competence of Civil Aviation College Students Based on Language Operator Combustion stability control of gasoline compression ignition (GCI) under low-load conditions: A review Research on the Psychological Distribution Delay of Artificial Neural Network Based on the Analysis of Differential Equation by Inequality Expansion and Contraction Method The Comprehensive Diagnostic Method Combining Rough Sets and Evidence Theory Study on Establishment and Improvement Strategy of Aviation Equipment Design of software-defined network experimental teaching scheme based on virtualised Environment Research on Financial Risk Early Warning of Listed Companies Based on Stochastic Effect Mode System dynamics model of output of ball mill The Model of Sugar Metabolism and Exercise Energy Expenditure Based on Fractional Linear Regression Equation Constructing Artistic Surface Modeling Design Based on Nonlinear Over-limit Interpolation Equation Optimal allocation of microgrid using a differential multi-agent multi-objective evolution algorithm About one method of calculation in the arbitrary curvilinear basis of the Laplace operator and curl from the vector function Numerical Simulation Analysis Mathematics of Fluid Mechanics for Semiconductor Circuit Breaker Cartesian space robot manipulator clamping movement in ROS simulation and experiment Effects of internal/external EGR and combustion phase on gasoline compression ignition at low-load condition Research of urban waterfront space planning and design based on children-friendly idea Characteristics of Mathematical Statistics Model of Student Emotion in College Physical Education Human Body Movement Coupling Model in Physical Education Class in the Educational Mathematical Equation of Reasonable Exercise Course Sensitivity Analysis of the Waterproof Performance of Elastic Rubber Gasket in Shield Tunnel Impact of Web Page House Listing Cues on Internet Rental Research on management and control strategy of E-bikes based on attribute reduction method A study of aerial courtyard of super high-rise building based on optimisation of space structure Exact solutions of (2 + 1)-Ablowitz-Kaup-Newell-Segur equation