1. Introduction
In addition to being subjected to the gravitational force of celestial bodies, a solar sail is also subjected to solar radiation pressure, which cannot vary arbitrarily in space. The solar sail spacecraft is a spacecraft that generates thrust using solar photons interacting with a high area ratio sail surface to achieve interplanetary travel [
1]. The solar sail spacecraft is actually a complex dynamical system in which orbital motion, attitude motion, and structural vibrations are coupled with each other. The standard for evaluating a solar sail dynamics model is that it should accurately reflect dynamical characteristics such as high flexibility, low stiffness, weak damping, and low fundamental frequency and modal density [
2]. The first successful application of the solar sail in the world was the IKAROS (Interplanetary Kite-craft Accelerated by Radiation of the Sun, IKAROS), which was launched in 2010 [
3].
At present, research on station-keeping control for solar sails is dominated by linear control methods. The basic principle of linear control is to first linearize the orbital dynamics equations, then use the reference orbit as the control command, and use various linear feedback control methods to design the station-keeping control law, such as linear state feedback control [
4,
5,
6], PID control [
7], and LQR control [
8,
9,
10,
11]. To update the control law for subsequent orbits, Moore and Ceriotti proposed a new method on the basis of the use of a Control Transition Matrix and linearization [
12]. Additionally, another major research approach for the station-keeping control of solar sails is auto disturbance rejection control (ADRC) [
13,
14,
15,
16,
17,
18,
19].
The solar sail spacecraft that utilizes hybrid small-thrust propulsion [
20,
21] is a new continuous small-thrust spacecraft combining the solar sail propulsion system and the solar electric propulsion system, which is suitable for complex orbital missions. On the one hand, it has the dual advantages of no energy consumption for solar sail propulsion and high efficiency for solar electric propulsion. On the other hand, it is one that overcomes the shortcomings of solar sails that cannot provide a propulsive part pointing in the direction of the sun. Simo and McInnes [
22] proposed a hybrid sail feedback linearized orbital controller, where they compensated for the nonlinear terms in the model with closed-loop feedback and designed a state error PD control law for the compensated linear system. Heiligers et al. [
23] designed a geosynchronous displaced orbit control method based on a hybrid propulsion solar sail, compared the fuel consumption of the hybrid propulsion solar sail and pulse control, and found that fuel consumption can be effectively reduced using a hybrid propulsion solar sail.
The newly launched Starlink satellite “V2 mini” utilizes argon Hall thrusters. The advent of this technology has made it possible to achieve hybrid orbit control without significantly increasing the mass of the solar sail spacecraft. The greatest advantage of argon electric thrusters over xenon or krypton Hall thrusters currently in use is their extremely low cost. As an example, 1 kg of high-purity xenon gas sells for tens of thousands of CNY, while the same quality of high-purity argon gas costs only a few CNY [
24].
The heliocentric displaced orbit is a non-Kepler periodic orbit displaced above the sun and formed by the mutual balance of spacecraft propulsion and gravity. Due to its special space position, it provides an ideal platform for deep space relay communication, sun–Earth observation, and other orbital missions [
25]. McInnes [
26] and Bookless [
27] summarized and analyzed the dynamic characteristics of different types of displaced orbits, and the results showed that the orbits under certain parameters are unstable, requiring station-keeping techniques to ensure the stable operation of a spacecraft in the target displaced orbit. Qian et al. [
28] adopted an LQR to design an orbit-keeping control law based on the linearized equation of state of the solar sail spacecraft near the reference heliocentric displaced orbit. Zhang et al. [
29] used the ADRC with low model dependence to design a station-keeping controller for a hybrid small-thrust spacecraft, and achieved good results. However, the parameter tuning is relatively cumbersome, and there is a significant overshoot, which means that the track position deviates significantly at one point. Chen et al. [
30] designed a control law based on a modified conditional integral sliding surface and combined it with an adaptive method to estimate uncertain parameters online.
Most of the above research on solar sail spacecraft in the heliocentric displaced orbit has linearized dynamic equations, resulting in high dependence on accurate models. Some studies have ignored the effects of model uncertainty due to modeling errors, complex deep space environments, etc. Meanwhile, single solar sail propulsion suffers from slow convergence speed and low control accuracy. And, complex conditions such as external disturbances, internal unmodeled dynamics, initial injection errors, and input saturation further increase the difficulty of the controller design. An adaptive high-performance sliding mode control strategy is proposed for the station-keeping problem under these complex conditions. Compared to the existing literature, the main contributions of this article include the following:
(1) A high-performance station-keeping controller, designed on the basis of an improved integral sliding mode and multivariate super-twisting sliding mode approaching law, has fast orbital tracking speed and high control accuracy;
(2) The weight coefficient matrix of the radial basis function (RBF) neural network is used as an uncertain parameter, and its adaptive law is designed using Lyapunov stability analysis. On this basis, an RBF neural network estimator is designed, which can estimate online and compensate for the combined disturbance acceleration constructed from external disturbances and internal unmodeled dynamics;
(3) Transforming the virtual control acceleration introduced into actual control variables and obtaining the optimal control variable with the objective of minimizing propellant consumption of the argon Hall thruster can solve the problem of low control accuracy caused by linearized dynamic equations.
(4) Although the use of hybrid low-thrust propulsion breaks the symmetry of the solar sail in the configuration, the proposed control strategy requires only a small amount of propellant consumption to effectively achieve the station-keeping control of the spacecraft.
2. Orbit Dynamic Model
2.1. Definition of Coordinate System and Attitude Angle
To derive the orbit dynamics model of a solar sail spacecraft utilizing hybrid small-thrust propulsion in a heliocentric displaced orbit, the following coordinate system is defined:
(1) Heliocentric displaced orbit coordinate system:
. As shown in
Figure 1, the origin is located at the center
of the sun. The
plane is parallel to the displaced orbital plane. The
-axis points to the vernal equinox of the J2000 ephemeris, and the
-axis points to the direction of the angular velocity of Earth’s revolution. The
-axis, the
-axis, and the
-axis form a right-handed coordinate system.
(2) Orbital coordinate system: . The origin is located at the center of the solar sail spacecraft. The -axis points in the direction of the position vector , and the -axis is perpendicular to the -axis in the plane of the displaced orbit and points in the direction of the spacecraft’s motion. The -axis, the -axis, and the -axis form a right-handed coordinate system. We assumed that the three axis unit vectors were , , and , respectively.
(3) Cylindrical coordinate system: . and denote the orbital radius and orbital height of the displaced orbit, respectively. is the angle between the projection of the position vector of the solar sail spacecraft in the plane and the -axis. We assumed that the unit vectors of the cylindrical coordinate system were , , and , respectively.
The transformation matrix of the orbital coordinate system to the heliocentric displaced orbit coordinate system is as follows:
where
;
.
Considering the special configuration of the solar sail spacecraft, new attitude angles are defined:
(1) Cone angle . Define to represent the unit normal vector of the sail surface. Then, the cone angle is the angle between the unit normal vector and the position vector . .
(2) Turning angle . The turning angle is the angle between the projection of the normal vector in the plane and the axis. .
2.2. Displaced Orbital Dynamics
In the system composed of the sun, Earth, and solar sail spacecraft, the heliocentric displaced orbit studied in this paper is relatively close to Earth and Earth’s gravity cannot be neglected. The displaced orbits studied in this paper are within the sphere of influence of a celestial body centered on the sun.
Considering the existence of model uncertainty and external disturbances in the displaced orbital motion system, the orbital dynamics equation is as follows:
where
represents the heliocentric gravitational constant;
is the propulsion acceleration under the action of solar radiation pressure;
is the electric propulsion acceleration,
;
is the acceleration introduced by the model uncertainty, mainly including the modeling error of solar radiation pressure propulsion caused by the special flexible configuration of the solar sail (sail membrane fold deformation, structural vibration, etc.); and
is the acceleration introduced by external disturbances, influenced by dynamic uncertainties such as solar storms and micrometeoroids in the environment.
Normalized units are used to simplify the analysis process. The Earth–sun distance is defined by unit length (1 AU), the solar mass is defined by unit mass (), and the angular velocity of Earth’s revolution is defined by unit angular velocity (). Then .
The solar radiation pressure propulsion acceleration is as follows:
The unit normal vector
of the sail surface can be expressed as follows:
It can also be expressed as a cylindrical coordinate, as follows:
There exists a geometrical relationship in the
coordinate system:
Taking the second order derivatives of
and
, respectively, gives the following:
Substituting Equations (5) and (6) into Equation (2), the orbital dynamics equations in the form of cylindrical coordinates are obtained as follows:
where
is the combined disturbance acceleration caused by model uncertainty and external disturbances,
.
Introducing the virtual control acceleration
, we can obtain the following:
where
3. High-Performance Station-Keeping Controller Design
For the dynamical model (10) of the heliocentric displaced orbit, a station-keeping controller based on an RBF neural network estimator, an improved integral sliding surface, and a multivariate super-twisting sliding mode approaching law are designed. The structure of the station-keeping control system is shown in
Figure 2.
3.1. RBF Neural Network Estimator Design
With excellent local nonlinear approximation, the RBF neural network is a neural network capable of approximating any continuous function with arbitrary accuracy. In practical engineering application systems, it is difficult to accurately measure the combined disturbance acceleration . The station-keeping control of the solar sail spacecraft is realized by designing an RBF neural network estimator to estimate and compensate for the combined disturbance acceleration in real time.
Assuming that the hidden layer has
neurons,
is a continuous nonlinear function. Denote
, and take the inputs to the neural network to be
.
is approximated by the RBF neural network:
where
represents the ideal RBF neural network weight;
denotes the approximation error, which is a very small real vector; and
is the radial basis function of the hidden layer.
The output of the Gaussian radial basis function of the
neuron of the hidden layer is as follows:
where
is the center vector of the Gaussian radial basis function, and
represents the width of the Gaussian basis function, which is a positive scalar.
is unknown, and the RBF neural network estimator
can be designed to fit it for real-time estimation.
Those marked with “^” indicate estimated values. To ensure the convergence of the errors and the real-time estimation, a reasonable adaptive law based on Lyapunov’s stability theory was designed for
, outlined in
Section 3.3.
According to Equations (13) and (15), the estimation error is as follows:
and
.
3.2. Control Law Design
For System (10), the cylindrical coordinate error
is defined as follows:
where
.
,
, and
denote nominal orbital radius, nominal orbital angle, and nominal displaced height, respectively.
The designed sliding surface
is as follows:
where
;
;
is the adjustment factor for the degree of the weakened integral (when
is a traditional integral sliding surface),
; and
is the thickness of the boundary layer,
. The saturation function is defined as follows:
When is inside the boundary layer, acts as a traditional integral. When is outside the boundary layer, the integral term is weakened. The extent to which it is weakened can be measured according to . It is clear that the extent of weakened demand satisfaction can be obtained through the regulating factor , . The integral term is weakened so that even under relatively large initial errors and disturbances, the control variable is not prone to jumping, the overshoot can be kept small, and the steady-state error is reduced.
Deriving Equation (18) and substituting Equations (10) and (13) into it, we obtain
The designed multivariate super-twisting sliding mode approaching law is as follows:
where
;
;
; and
.
Substituting Equation (21) into Equation (20), we can obtain the following:
The designed adaptive sliding mode control law for station-keeping control is as follows:
3.3. Stability Analysis
In the process of designing the controller, the adaptive law is obtained based on the Lyapunov function, which ensures the stability of the corresponding system. And, the stability of the control system, shown in
Figure 2, has not been proven via theoretical derivation. Therefore, in this section, the stability of the above closed-loop system will be analyzed based on the Lyapunov function.
Lemma 1 [31]. For a nonlinear system ,
assuming the existence of a positive definite function ,
this satisfies the following:
where ; ; ; ;
and the nonlinear system is finite-time stable.
The convergence time can be given by Theorem 1. For the orbit dynamics system (10) of a solar sail spacecraft utilizing hybrid small-thrust propulsion in a heliocentric displaced orbit, if its control law is designed as in Equation (23), and if the values of the parameters such as , , , , , and are appropriate, the control system can converge to the sliding surface in finite-time, and the orbital tracking error of the system can converge asymptotically to zero.
Proof of Theorem 1. The proof process will be divided into two steps. Step 1 proves that the system state converges to the sliding surface in finite-time, and Step 2 proves that after reaching the sliding surface, the orbital tracking error of the system converges asymptotically to zero.
Substituting Equation (23) into Equation (20), we can obtain the following:
where
is the estimation error of the RBF neural network estimator, assuming that it satisfies
and
.
Step 1. The Lyapunov function is chosen as follows:
Equation (26) can be written as , where
, .
satisfies the following:
Since
, it follows that
Deriving the derivative for
, we can obtain the following:
where
Let
; then, the adaptive law can be designed as follows:
where
represents the adjustable coefficient associated with the adaptive law,
;
represents the pseudo-converse operation. Thus,
If
and
are positive definite symmetric matrices, then
is a negative definite. It is necessary to satisfy
.
Using Equation (28), it can be derived that
where
From Lemma 1, if the values of the parameters such as , , , , , and are appropriate, the control system can converge to the sliding surface in finite-time.
Step 2. The system reaches the sliding surface,
, satisfying
; then,
When
tends to zero, there exists the following:
The Lyapunov function is chosen as follows:
The derivation of Equation (36) yields the following:
is equivalent to . According to LaSalle’s invariance principle, it can be known that () is the global asymptotically stable equilibrium point of Equation (34), i.e.,. The orbital tracking error of the system converges asymptotically to zero.
This completes Proof of Theorem 1. □
3.4. Control Variable Conversion
To achieve station-keeping control, the virtual control acceleration needs to be converted into solar sail cone angle
, turning angle
, and solar electric propulsion acceleration
after the following is obtained:
. From Equation (12), the solar electric propulsion acceleration components can be derived as follows:
where
,
,
, and
are real-time measurement values; the light pressure factor
can be found via Equation (39).
where
,
, and
are nominal parameters of the displaced orbit.
To more directly describe the cost of station-keeping control, real-time mass changes are applied here. The real-time mass change in the electric propellant satisfies the following:
where
is the electric thrust,
, and
;
is the specific impulse of electric propulsion; and
is the standard gravitational acceleration of Earth.
The propellant consumption can be reduced by solving the minimum value of electric propulsion acceleration using the following expression:
where
and
, respectively, are the cone angle and turning angle, which are optimal solutions for propellant consumption. On the basis of the optimal solutions, the electric thrust
can be obtained.