Research Article  Open Access
Robust Linear Quadratic Regulator via Sliding Mode Guidance for Spacecraft Orbiting a Tumbling Asteroid
Abstract
Aiming to ensure the stability of the spacecraft with multiuncertainties and mitigate the threat of the initial actuator saturation, a Robust Linear Quadratic Regulator (RLQR) via sliding mode guidance (SMG) for orbiting a tumbling asteroid is proposed in this paper. The orbital motion of the spacecraft near a tumbling asteroid is modelled in the bodyfixed frame considering the sunrelative effects, and the orbiting control problem is formulated as a stabilization of a nonlinear timevarying system. RLQR based on the adaptive feedback linearization is proposed to stabilize the spacecraft orbiting with the uncertainties of the asteroid’s rotation and gravitational field. In order to avoid the initial actuator saturation, SMG is applied to generate the transition process trajectory of the closedloop system. The effectiveness of the proposed control scheme is verified by the simulations of orbiting the asteroid Toutatis 4179.
1. Introduction
The exploration of asteroids has become a hot topic of interests. Equilibria and periodic orbits near asteroids are usually highly unstable due to the irregular shapes and the complex rotation of the asteroids. Hence, the control of a spacecraft in closed proximity of an asteroid is among the most challenging technical problems in the exploration of asteroids [1].
Many researches have paid great attention to the orbital control of the spacecrafts near asteroids. Sawai et al. [2] presented a control based on the onedimensional altimetry measurements to stabilize a hovering spacecraft. Broschart and Scheeres [3] investigated the stability of realistic hovering control in the bodyfixed and the inertial reference frames, respectively. Then, Broschart and Scheeres [4] proposed the sufficient conditions for a deadband controller to bound spacecraft hovering motion in timeinvariant Lagrangian dynamical systems. Wie [5] presented the dynamic modeling and control analysis of multiple gravity tractors in halo orbits. Furfaro et al. [6] investigated a novel closedloop autonomous guidance law based on multiple sliding surfaces for the soft landing of the spacecraft on the designated point on the asteroid. Liu et al. [7] presented orbital control law for the spacecraft, which consists of PD controller and a nonsingular terminal sliding mode controller, to track the soft landing trajectory. Guelman [8] investigated a simple threedimensional guidance law for the orbit transfer to a quasicircular orbit about a rotating small celestial body using continuous thrust. In these previous contributions, asteroids are assumed to be a pure rotation configuration about the axis with its maximum inertia. However, the rotation of asteroids is very complex in practice and lots of them are timevarying [9]. The asteroid with nonprincipal axis rotational motion is socalled “tumbling asteroid” [10]. Nazari et al. [11] investigated the observer based bodyframe hovering control over a tumbling asteroid, which is based on the timevarying LQR or the LyapunovFloquet transformation and timeinvariant LQR, respectively. However, the rotation of the asteroid needs to be modelled in his control, which will cost a huge amount of telescope observing time to be determined [12], and the effect of the gravitational fields uncertainty is not analyzed. Hence, designing a robust control law for the spacecraft under the multiuncertainties of the asteroid’s gravitational field and rotation is urgent.
The initial actuator saturation problem also needs to be considered. The output of the controller may have a high peak at the start, whose magnitude depends on the size of the initial error. As the continuous adjustable thrust is small, the control output may rise over the limitation of the thrust and threaten the stability of the system [13]. Aiming to avoid this problem, arranging the transition process is a good idea to acquire an acceptable controlled quantity. In many researches of the active disturbance rejection controllers, Tracking Differentiator (TD) is applied to arrange the transition process to avoid much larger control outputs [14–17]. However, the stability of the TD is difficult to be approved and the magnitude of the control output via TD is hardly adjustable. The sliding mode control is an important improvement of the control theory [18]. Based on it, the sliding mode guidance (SMG) is confirmed due to its good performance in the orbital control of the spacecraft near asteroids [19]. We apply the SMG to arrange the transition process of the controller. By doing this, the max magnitude of the control output becomes adjustable and the threat of the initial actuator saturation problem is reduced.
In this paper, we proposed RLQR via the SMG for stabilizing the orbit of the spacecraft around a tumbling asteroid, which does not need to model the rotation of the asteroid. Firstly, the spacecraft orbital motion near a tumbling asteroid is modeled as a restricted threebody problem. RLQR based on the adaptive feedback linearization is proposed. The feedback linearization consists of a feedforward control, which is based on the spheric harmonic coefficient model of the asteroids gravitational field, and an adaptive compensator to ensure the robust stability against model uncertainties. Aiming to mitigate the threat of the initial actuator saturation, the SMG is applied to arrange the transition process of the RLQR. Simulations of orbiting the Toutatis 4769 are performed to verify the effectiveness of the proposed controller. The results of applying the RLQR with and without the SMG are compared to show the advantages.
The rest of the paper is organized in the following form. In Section 2, the orbital dynamic of spacecraft orbital motion is modeled and the problem formulation for control is proposed. In Section 3, the RLQR based on the adaptive feedback linearization is proposed. The SMG is applied to arrange the transition process of the proposed RLQR. In Section 4, simulations of orbiting the Toutatis 4769 are performed. Conclusions are drawn in Section 5.
2. Problem Formulation
For the relative orbital motion of a spacecraft near an asteroid, the dynamic model includes the nonspherical gravity field of the asteroid and the solar radiation pressure (SRP) [1]. In most previous researches, the orbital dynamic of the spacecraft was formulated into two regimes: the gravity dominated regime, in which the effects on the spacecraft from the sun were ignored, and the solar dominated regime, in which the gravity of the asteroid was out of consideration. Few researches considered both of the effects from the sun and the asteroid [20]. In this section, the orbital dynamic of the spacecraft near a tumbling asteroid is modeled as a restricted threebody problem with the joint perturbations.
2.1. Orbital Dynamic of the Spacecraft near Tumbling Asteroids
Before the modeling, two relative frames need to be defined. As shown in Figure 1, the inertial frame centered at the mass center of the asteroid. The axis is parallel to the sunline. The axis is coinciding with the direction of the velocity vector of the asteroid. The , , and axis compose the righthanded coordinate system. The bodyfixed frame fixes on asteroid with the origin coinciding with the mass center of the asteroid. The , , and are coinciding with the axis of the asteroids maximum, minimum, and intermediate moment of inertia, respectively.
The orbital motion of spacecraft in the inertial frame can be expressed aswhere is the orbital angular velocity vector of the asteroid, is the vector of the sun to the asteroid, is the vector of the sun to the probe, is the vector of the asteroid to the probe, , is the gravitational parameter of the sun, is the gravitational acceleration of the asteroid, is the solar radiation pressure on the probe, and is the vector of control.
Assume that the asteroid orbits the sun in a circle; (1) can be simplified as
Transform (2) into the bodyfixed frame; it can be expressed aswhere is the rotation vector of the asteroid and is the sunrelative effect on the spacecraft.
2.1.1. SunRelative Effects
Define , and the following equation can be established according to the assumption of asteroids circle orbit:
Apply the Taylor expansion; an approximate linearization of (4) can be expressed as
Considering that the main effect of the SRP arises from its firstorder perturbations [21], can be defined aswhere is a defined constant parameter of SRP.
Define the transformation matrix from the inertial frame to the bodyfixed frame as and the rotation vector as . Then, (6) can be rearranged aswhere
As usually , , and the elements of transformation matrix , , hence (9) can be simplified as
2.1.2. Asteroid’s Gravitational Field
The gravitational acceleration of the asteroid can be modelled by the gradient of the gravitational potential function of the asteroid as
Different from the regular large objects, the bodies of asteroids are very irregular. Several methods have been investigated to approach the gravitational field of asteroids, such as the Spherical Harmonic Expansion Model (SHEM) [22], the polyhedron model [23], and the ellipsoidal harmonic expansion model [24]. The secondorder SHEM, which models the asteroid’s perturbation potential with the most significant gravity coefficients and is widely used in the previous researches, is selected to formulate the nominal gravitational field in this paper. It can be expressed aswhere is the gravitational parameter of the asteroid, is the normalizing radius, and are the spherical harmonic gravity coefficient, and , , and are the radius, latitude, and longitude of the field point, respectively.
Regard the asteroid as a threeaxis ellipsoid with three main axes , and the spherical harmonic gravity coefficients and can be estimated by the three axes of the asteroid as
Then, the real gravity of the asteroid can be formulated aswhere is the uncertainty of the asteroids gravitational model.
Hypothesis 1. is bounded and the constant exists, which satisfy the following conditions:
2.2. Problem Formulation for Controller Design
The stable hovering in the inertial frame is equivalent to a stable orbiting in the bodyfixed frame. Hence, the stable orbiting condition in the bodyfixed frame is expressed as
Define the error vector as
Substitute (3) into (17); the time derivative of the error state is expressed as
Then, the stable orbiting control problem can be formulated as the stabilization of the nonlinear system equation (18).
The rotation of the asteroid is difficult to be modelled before the close proximity orbiting observation. Hence, the following assumptions need to be declared before the controller design:(i)The rotation vector is unmodeled, but its upper bound is known.(ii)The current state of the rotation can be observed, so that the error vector can be applied to the feedback.(iii)The position and the velocity of the spacecraft in the bodyfixed frame are provided by the navigation system.
3. Robust Linear Quadratic Regulator via Sliding Mode Guidance
3.1. Robust Linear Quadratic Regulator
The nonlinear system equation (18) can be transformed into (19). Consider
Because the rotation vector is unmodeled, the timevarying matrix is uncertain. However, the upper bound of is known. Hence, we can define the upper bound matrix of as
Design RLQR for the stabilization of the nonlinear system equation (19) as
is a feedforward control as
is a linear quadratic regulator, which minimizes the cost function as (25). It can be expressed aswhere and are positive weight matrices and is a positive symmetric matrix, which is the solution of the Riccati equation as
is an adaptive compensator, which can be expressed aswhere is the adaptive parameter and is an adjustable field in the error space that contains the origin.
Theorem 1. Considering the spacecraft orbiting problem which is formulated as the stabilization of a nonlinear system as (19), define the error state as (18) and apply the robust LQR as (24). If Hypothesis 1 is satisfied, the closedloop system is globally uniformly stable.
Proof. Define a Lyapunov function candidate aswhere is a constant, which make all satisfy the condition asThe time derivative of can be expressed aswhere the function is defined as and the function is defined as .
As , , and are positive matrices, for all , . According to the Lyapunov stability theorem, the nonlinear system as (18) under the control as (22) is uniformly stable. The bound of the steady state error depends on the adjustable field .
3.2. RLQR via the Sliding Mode Guidance
The proposed RLQR could ensure the robust stability of the orbiting. However, the proposed controller has the potential threat from the actuator saturation. Firstly, depending on the error state, the control output may have a high peak at the start. Secondly, the gain of the adaptive compensator may be too conservative, which is much larger than necessary. As the adjustable continuous thrusters only could provide a small thrust, the saturation of the thrust may threaten the safety of the spacecraft. To avoid the initial actuator saturation problem of the RLQR, the sliding mode guidance is applied to arrange the transition process of the controller.
3.2.1. The Sliding Mode Guidance
Define the outputs of the sliding mode guidance as and , and their initial values are
The sliding mode guidance is expressed aswhere is the sliding mode function as and and are selected parameters of the SMG.
Define a Lyapunov function candidate as
Then, the time derivative of is expressed as
Obviously, only when . According to the Lyapunov stability theorem, the SMG is globally asymptotically stable. When , and .
3.2.2. The RLQR Based on the Virtual Error Vector
Define the virtual error vector as
The time derivative of the virtual error can be expressed as
As the SMG is globally asymptotically stable and the matrix is bounded the constant exists, which satisfy the condition as
Then, design the RLQR based on the virtual error as
Theorem 2. Considering the spacecraft orbiting problem which is formulated as the stabilization of a nonlinear system as (18), define the sliding mode guidance as (32) and the virtual error vector as (17) and (35), and apply the robust LQR based on the virtual error as (38). If Hypothesis 1 is satisfied, the closedloop system is globally uniformly stable.
Proof. Define a Lyapunov function candidate aswhere is a constant, which make all satisfy the condition asThen, the time derivative of is represented asAs , , and are positive matrices, for all , . Hence, when , , for all , . According to the Lyapunov stability theorem, the RLQR via SMG could make the system as (18) uniformly stable. The bound of the steady state error is .
4. Simulations
We specifically use the numerical data of the asteroid Toutatis 4179 in the simulation. To make the simulation more realistic, the gravitational field of the asteroid is modeled by the polyhedron method [23] based on Hudson’s 33996face shape data [25], which can be expressed as (43). The secondorder SHEM is applied as the nominal model in the controller in the same time:where is the gravitational constant, is the constant density of the asteroid, is a vector from the field point to an arbitrary point on each edge, is a vector from the field point to an arbitrary point on each face, is a dyad defined in terms of the face and edge normal vectors associated with each edge, is the outer product of face normal vectors, is a logarithmic term expressing the potential of 1D straight wire, and is the solid angle subtended by a face when viewed from the field point.
According to the shape data, the maximum, intermediate, and minimum major axes of Toutatis 4179 are km. The rotation of Toutatis in two periods is shown in Figure 2, which starts at rad/s. The initial attitude of the asteroid expressed by 321 Euler angle is . The initial transformation matrix determined is by the initial attitude and the rotation of the asteroid. The gravitational parameters of Toutatis are shown in Table 1.

4.1. The Orbiting under the Control of the RLQR
The weight matrixes of RLQR are and . The adjustable field of the error state is . The initial position of the spacecraft in the bodyfixed frame is m. The initial velocity of the spacecraft in the bodyfixed frame is m/s.
As illustrated in Figure 3, the error vector under RLQR control converges to zero quickly. Then, the errors keep near to the zero axes during the remaining times of the simulation and the system becomes stable. The outputs of the RLQR are shown in Figure 4. In the first several seconds, the control acceleration has a high peak, and the adaptive gain grows quickly (as shown in Figure 5). Then, the control acceleration decreases and becomes very small that mainly restrains the effect of the sun and the gravity of the asteroid. The adaptive gain is almost stable at 0.2213. In the simulation, the control cost of the steadystable orbiting is about 0.0228 m/s per hour, which is acceptable for the deep space missions lasting months. Through these simulations, the effectiveness of the RLQR has been verified.
4.2. The Orbiting under the RLQR via the SMG
The parameters of the SMG are and . The adjustable field of the virtual error state is .
The error states under the control of the RLQR via SMG are shown in Figure 6. The convergence of the error states takes more time than the ones in Figure 3. The reduction of convergence speed is the price of the decrease of the control magnitude. As illustrated in Figure 7, the max value of the control acceleration is 0.046 m/s^{2}, which is much less than the one in Figure 4. What is more, the max magnitude of the control acceleration can be adjusted by the parameters and . Hence, the risk of the initial actuator saturation problem caused by the large initial error is mitigated. The gain of the adaptive compensator with the SMG is much smaller than the one without the SMG (see Figure 8). The control cost of the steadystable orbiting is similar to the one without the SMG.
Through the simulations in Figures 6–8, it was verified that the RLQR via the SMG could stabilize the orbiting and limit the magnitude of the control output. The threat of the initial actuator saturation has been mitigated.
4.3. The Proposed Controller under Different Magnitudes of Uncertainty
To analyze the effects of the gravity model’s uncertainties on the proposed algorithm, we did more simulations which apply no gravity model, the secondorder SHEM, the fourthorder SHEM, and the polyhedron model (no uncertainty), respectively. Adding the order of the SHEM to the fourth, the gravitational potential of the asteroid can be expressed as
The spheric harmonic coefficient can be determined as
In the simulation, the initial position of the spacecraft is km in the bodyfixed frame. The initial velocity of the spacecraft is m/s in the bodyfixed frame. The simulation results are shown as Figures 9–14.
As shown in Figures 9–11, the errors in three axes are converged near zero no matter which model applied. Also, the error curves controlled with different uncertainties have little difference. As shown in Figures 12–14, the control outputs are similar too.
The reason of these results may be shown in Figures 15–17. The asteroid’s gravitational acceleration is too small, which is in the magnitude of . Although we limited the control output, it is still much larger than the asteroid’s gravity.
As a conclusion, the uncertainty of the asteroid’s gravity has little effect on the proposed control algorithm. The proposed controller can be used in the orbiting around the asteroids which are not well mapped in terms of gravity field.
5. Conclusion
In this paper, RLQR via SMG is proposed to ensure the robust stability of the spacecraft orbiting a tumbling asteroid. Based on the adaptive feedback linearization, the controller is able to suppress the effects of the multiuncertainties in the asteroids rotation and gravity model. Applying the SMG to arrange the transition process, the max magnitude of the control output is limited and threatening of the initial actuator saturation of the controller is mitigated. Also, the control cost of the steadystable orbiting is acceptable for the deep space missions lasting months. As a consequence, the proposed control scheme can be properly selected to maintain spacecraft orbiting a tumbling asteroid and does not need to model the rotation of the asteroid.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This work is supported by the National Basic Research Program of China (973 Program) under Grant 2012CB720000.
References
 D. J. Scheeres, “Close proximity dynamics and control about asteroids,” in Proceedings of the American Control Conference (ACC '14), pp. 1584–1598, Portland, Ore, USA, June 2014. View at: Google Scholar
 S. Sawai, D. J. Scheeres, and S. B. Broschart, “Control of hovering spacecraft using altimetry,” Journal of Guidance, Control, and Dynamics, vol. 25, no. 4, pp. 786–795, 2002. View at: Publisher Site  Google Scholar
 S. B. Broschart and D. J. Scheeres, “Control of hovering spacecraft near small bodies: application to asteroid 25143 Itokawa,” Journal of Guidance, Control, and Dynamics, vol. 28, no. 2, pp. 343–354, 2005. View at: Publisher Site  Google Scholar
 S. B. Broschart and D. J. Scheeres, “Boundedness of spacecraft hovering under deadband control in timeinvariant systems,” Journal of Guidance, Control, and Dynamics, vol. 30, no. 2, pp. 601–610, 2007. View at: Publisher Site  Google Scholar
 B. Wie, “Dynamics and control of gravity tractor spacecraft for asteroid deflection,” Journal of Guidance, Control, and Dynamics, vol. 31, no. 5, pp. 1413–1423, 2008. View at: Publisher Site  Google Scholar
 R. Furfaro, D. Cersosimo, and D. R. Wibben, “Asteroid precision landing via multiple sliding surfaces guidance techniques,” Journal of Guidance, Control, and Dynamics, vol. 36, no. 4, pp. 1075–1092, 2013. View at: Publisher Site  Google Scholar
 K. Liu, F. Liu, S. Wang, and Y. Li, “Finitetime spacecraft's soft landing on asteroids using PD and nonsingular terminal sliding mode control,” Mathematical Problems in Engineering, vol. 2015, Article ID 510618, 10 pages, 2015. View at: Publisher Site  Google Scholar
 M. Guelman, “Closedloop control of close orbits around asteroids,” Journal of Guidance, Control, and Dynamics, vol. 38, no. 5, pp. 854–860, 2015. View at: Publisher Site  Google Scholar
 M. Kaasalainen, “Interpretation of lightcurves of precessing asteroids,” Astronomy and Astrophysics, vol. 376, no. 1, pp. 302–309, 2001. View at: Publisher Site  Google Scholar
 P. Pravec, A. W. Harris, P. Scheirich et al., “Tumbling asteroids,” Icarus, vol. 173, no. 1, pp. 108–131, 2005. View at: Publisher Site  Google Scholar
 M. Nazari, R. Wauson, T. Critz, E. A. Butcher, and D. J. Scheeres, “Observerbased bodyframe hovering control over a tumbling asteroid,” Acta Astronautica, vol. 102, pp. 124–139, 2014. View at: Publisher Site  Google Scholar
 P. Pravec, P. Scheirich, J. Ďurech et al., “The tumbling spin state of (99942) Apophis,” Icarus, vol. 233, pp. 48–60, 2014. View at: Publisher Site  Google Scholar
 Z. Zhu, Y. Xia, and M. Fu, “Adaptive sliding mode control for attitude stabilization with actuator saturation,” IEEE Transactions on Industrial Electronics, vol. 58, no. 10, pp. 4898–4907, 2011. View at: Publisher Site  Google Scholar
 G. Feng, Y.F. Liu, and L. P. Huang, “A new robust algorithm to improve the dynamic performance on the speed control of induction motor drive,” IEEE Transactions on Power Electronics, vol. 19, no. 6, pp. 1614–1627, 2004. View at: Publisher Site  Google Scholar
 P. Dong, G.J. Ye, J. Wu, J.M. Yang, and Y.R. Chen, “Autodisturbance rejection controller in the wind energy conversion system,” in Proceedings of the 4th International Power Electronics and Motion Control Conference (IPEMC '04), vol. 2, pp. 878–881, IEEE, Xi'an, China, August 2004. View at: Google Scholar
 J. H. Ruan, Z. W. Li, F. Y. Zhou, and Y. B. Li, “ADRC based ship tracking controller design and simulations,” in Proceedings of the IEEE International Conference on Automation and Logistics (ICAL '08), pp. 1763–1768, Qingdao, China, September 2008. View at: Publisher Site  Google Scholar
 B. Gao, J. Shao, and X. Yang, “A compound control strategy combining velocity compensation with ADRC of electrohydraulic position servo control system,” ISA Transactions, vol. 53, no. 6, pp. 1910–1918, 2014. View at: Publisher Site  Google Scholar
 K. D. Young, V. I. Utkin, and U. Ozguner, “A control engineer's guide to sliding mode control,” IEEE Transactions on Control Systems Technology, vol. 7, no. 3, pp. 328–342, 1999. View at: Publisher Site  Google Scholar
 Z. Zexu, W. Weidong, L. Litao et al., “Robust sliding mode guidance and control for soft landing on small bodies,” Journal of the Franklin Institute. Engineering and Applied Mathematics, vol. 349, no. 2, pp. 493–509, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 D. J. Scheeres, “Orbit mechanics about asteroids and comets,” Journal of Guidance, Control, and Dynamics, vol. 35, no. 3, pp. 987–997, 2012. View at: Publisher Site  Google Scholar
 J. W. McMahon and D. J. Scheeres, “New solar radiation pressure force model for navigation,” Journal of Guidance, Control, and Dynamics, vol. 33, no. 5, pp. 1418–1428, 2010. View at: Publisher Site  Google Scholar
 D. J. Scheeres, B. G. Williams, and J. K. Miller, “Evaluation of the dynamic environment of an asteroid: applications to 433 Eros,” Journal of Guidance, Control, and Dynamics, vol. 23, no. 3, pp. 466–475, 2000. View at: Publisher Site  Google Scholar
 R. A. Werner and D. J. Scheeres, “Exterior gravitation of a polyhedron derived and compared with harmonic and mascon gravitation representations of asteroid 4769 Castalia,” Celestial Mechanics and Dynamical Astronomy, vol. 65, no. 3, pp. 313–344, 1996. View at: Google Scholar
 R. Garmier and J.P. Barriot, “Ellipsoidal harmonic expansions of the gravitational potential: theory and application,” Celestial Mechanics and Dynamical Astronomy, vol. 79, no. 4, pp. 235–275, 2001. View at: Publisher Site  Google Scholar  MathSciNet
 R. S. Hudson, S. J. Ostro, and D. J. Scheeres, “Highresolution model of asteroid 4179 Toutatis,” Icarus, vol. 161, no. 2, pp. 346–355, 2003. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Peng Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.