Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Chiral Organophosphorus Pharmaceuticals: Properties and Application
Next Article in Special Issue
Assessment of Non-Coplanar Maneuver Parameters and Perturbing Accelerations Using the Minimal Number of Observations
Previous Article in Journal
Results of Hyperbolic Ricci Solitons
Previous Article in Special Issue
Heuristic Technique for the Search of Interception Trajectories to Asteroids with the Use of Solar Sails
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Station-Keeping Control Strategy for a Symmetrical Spacecraft Utilizing Hybrid Low-Thrust Propulsion in the Heliocentric Displaced Orbit

1
School of Astronautics, Harbin Institute of Technology, Harbin 150001, China
2
Aerospace System Engineering Shanghai, Shanghai 201109, China
3
Shanghai Institute of Satellite Engineering, Shanghai 201109, China
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(8), 1549; https://doi.org/10.3390/sym15081549
Submission received: 14 June 2023 / Revised: 23 July 2023 / Accepted: 4 August 2023 / Published: 6 August 2023
(This article belongs to the Special Issue Advances in Mechanics and Control II)

Abstract

:
The solar sail spacecraft utilizing a hybrid approach of solar sail and solar electric propulsion in the heliocentric displaced orbit is affected by external disturbances, internal unmodeled dynamics, initial injection errors, and input saturation. To solve the station-keeping control problem under such complex conditions, an adaptive control strategy is proposed. First, the dynamical equations of the spacecraft utilizing hybrid low-thrust propulsion in the cylindrical coordinate system are derived. Second, the combined disturbance acceleration introduced by external disturbances and internal unmodeled dynamics is constructed, and a radial basis function neural network estimator is designed to estimate it online in real time. Third, an adaptive high-performance station-keeping controller based on an improved integral sliding surface and multivariate super-twisting sliding mode approaching law is designed. Then, stability analysis is conducted using Lyapunov theory, adaptive laws are designed, and the introduced virtual control accelerations are converted into actual control variables. Finally, simulations are conducted under different simulation conditions based on the disturbance sources. The results show that although the use of hybrid low-thrust propulsion breaks the symmetry of the solar sail in configuration, the proposed control strategy can effectively achieve the station-keeping and disturbance estimation of the spacecraft with only a small amount of propellant consumed and position tracking errors up to decimeters.

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: o - x y z . As shown in Figure 1, the origin is located at the center o of the sun. The x o y plane is parallel to the displaced orbital plane. The o x -axis points to the vernal equinox of the J2000 ephemeris, and the o z -axis points to the direction of the angular velocity of Earth’s revolution. The o x -axis, the o y -axis, and the o z -axis form a right-handed coordinate system.
(2) Orbital coordinate system: o s s - x s s y s s z s s . The origin is located at the center o s s of the solar sail spacecraft. The o s s x s s -axis points in the direction of the position vector r s , and the o s s y s s -axis is perpendicular to the o s s x s s -axis in the plane of the displaced orbit and points in the direction of the spacecraft’s motion. The o s s x s s -axis, the o s s y s s -axis, and the o s s z s s -axis form a right-handed coordinate system. We assumed that the three axis unit vectors were i s , j s , and k s , respectively.
(3) Cylindrical coordinate system: ρ s , θ s , z s . ρ s and z s denote the orbital radius and orbital height of the displaced orbit, respectively. θ s is the angle between the projection of the position vector r s of the solar sail spacecraft in the x o y plane and the o x -axis. We assumed that the unit vectors of the cylindrical coordinate system were ρ ^ , θ ^ , and z ^ , respectively.
The transformation matrix of the orbital coordinate system to the heliocentric displaced orbit coordinate system is as follows:
C = cos φ s sin φ s 0 0 0 1 sin φ s cos φ s 0
where cos φ s = ρ s r s ; sin φ s = z s r s .
Considering the special configuration of the solar sail spacecraft, new attitude angles are defined:
(1) Cone angle α s . Define n to represent the unit normal vector of the sail surface. Then, the cone angle α s is the angle between the unit normal vector n and the position vector r s . α s π / 2 π / 2 .
(2) Turning angle δ s . The turning angle δ s is the angle between the projection of the normal vector n in the y s s o s s z s s plane and the o s s z s s axis. δ s π π .

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:
r ¨ s = μ s r s 3 r s + a s + a e p + a Δ f + a d
where μ s represents the heliocentric gravitational constant; a s is the propulsion acceleration under the action of solar radiation pressure; a e p is the electric propulsion acceleration, a e p = [ a e p , x a e p , y a e p , z ] ; a Δ f 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 a d 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 ( M s = 1 ), and the angular velocity of Earth’s revolution is defined by unit angular velocity ( ω e = 1 ). Then μ s = 1 .
The solar radiation pressure propulsion acceleration is as follows:
a s = β s 1 r s 2 n cos 2 α s
The unit normal vector n of the sail surface can be expressed as follows:
n = cos α s sin α s sin δ s sin α s cos δ s i s j s k s
It can also be expressed as a cylindrical coordinate, as follows:
n = cos α s cos φ s sin α s cos δ s sin φ s ρ ^ + sin α s sin δ s θ ^ + cos α s sin φ s sin α s cos δ s cos φ s z ^
There exists a geometrical relationship in the o - x y z coordinate system:
r s = ρ s + z s
Taking the second order derivatives of ρ s and z s , respectively, gives the following:
d 2 ρ s d t 2 = d 2 ρ s d t 2 ρ s d θ s d t 2 ρ ^ + 2 d ρ s d t d θ s d t + ρ s d 2 θ s d t 2 θ ^
d 2 z s d t 2 = d 2 z s d t 2 z ^
Substituting Equations (5) and (6) into Equation (2), the orbital dynamics equations in the form of cylindrical coordinates are obtained as follows:
ρ ¨ s = ρ s θ ˙ s 2 ρ s r s 3 + β s cos 3 α s r s 3 ρ s sin α s cos 2 α s cos δ s r s 3 z s + a e p , x cos θ s + a e p , y sin θ s + a f d , 1 θ ¨ s = 2 ρ s ρ ˙ s θ ˙ s + β s sin α s cos 2 α s sin δ s ρ s r s 2 a e p , x sin θ s ρ s + a e p , y cos θ s ρ s + a f d , 2 z ¨ s = z s r s 3 + β s cos 3 α s r s 3 z s + β s sin α s cos 2 α s cos δ s r s 3 ρ s + a e p , z + a f d , 3
where a f d is the combined disturbance acceleration caused by model uncertainty and external disturbances, a f d = a Δ f + a d = a f d , 1 a f d , 2 a f d , 3 T .
Introducing the virtual control acceleration u = u 1 u 2 u 3 T , we can obtain the following:
ρ ¨ s θ ¨ s z ¨ s = a ρ s + a f d , 1 a θ s + a f d , 2 a z s + a f d , 3 + u 1 u 2 u 3
where
a ρ s = ρ s θ ˙ s 2 ρ s r s 3 a θ s = 2 ρ s ρ ˙ s θ ˙ s a z s = z s r s 3
u 1 = β s cos 3 α s r s 3 ρ s sin α s cos 2 α s cos δ s r s 3 z s + a e p , x cos θ s + a e p , y sin θ s u 2 = β s sin α s cos 2 α s sin δ s ρ s r s 2 a e p , x sin θ s ρ s + a e p , y cos θ s ρ s u 3 = β s cos 3 α s r s 3 z s + β s sin α s cos 2 α s cos δ s r s 3 ρ s + a e p , z

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 a f d . 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 n neurons, a f d is a continuous nonlinear function. Denote Ξ = [ ρ s θ s z s ] T , and take the inputs to the neural network to be x = [ Ξ ; Ξ ˙ ] . a f d , i ( x ) ( i = 1 , 2 , 3 ) is approximated by the RBF neural network:
a f d , i ( x ) = w i T χ i ( x ) + ν i
where w i n represents the ideal RBF neural network weight; ν i denotes the approximation error, which is a very small real vector; and χ i ( x ) n is the radial basis function of the hidden layer.
The output of the Gaussian radial basis function of the j th   ( j = 1 , 2 , , n ) neuron of the hidden layer is as follows:
χ i j ( x ) = exp ( x c i j ) T ( x c i j ) 2 σ i j 2
where c i j is the center vector of the Gaussian radial basis function, and σ i j represents the width of the Gaussian basis function, which is a positive scalar.
a f d , i is unknown, and the RBF neural network estimator a ^ f d , i ( x ) can be designed to fit it for real-time estimation.
a ^ f d , i ( x ) = w ^ i T χ i ( x )
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 w ^ i , outlined in Section 3.3.
According to Equations (13) and (15), the estimation error is as follows:
e a = a f d , i ( x ) a ^ f d , i ( x ) = w i T w ^ i T χ i ( x ) + ν i
and e ˙ a = w ^ ˙ i T χ i ( x ) .

3.2. Control Law Design

For System (10), the cylindrical coordinate error e = [ e 1 e 2 e 3 ] T is defined as follows:
e = Ξ Ξ d
where Ξ d = [ ρ s , d θ s , d z s , d ] T . ρ s , d , θ s , d , and z s , d denote nominal orbital radius, nominal orbital angle, and nominal displaced height, respectively.
The designed sliding surface s i ( i = 1 , 2 , 3 ) is as follows:
s i = e ˙ i + k 0 i e i + k 1 i w i w ˙ i = k 1 i w i + ε i 1 + τ i sat s i τ i s i e ˙ i / k 0 i
where k 0 i > 0 ; k 1 i > 0 ; τ i is the adjustment factor for the degree of the weakened integral (when τ i = 1 is a traditional integral sliding surface), τ i > 1 ; and ε i is the thickness of the boundary layer, ε i > 0 . The saturation function is defined as follows:
sat s i = 1 s i > ε i s i / ε i s i ε i 1 s i < ε i
Then,
w ˙ i = e i ( s i ε i ) ( 1 + τ i ) / k 0 i s i > ε i e i s i ε i e i + ( s i ε i ) ( 1 + τ i ) / k 0 i s i < ε i
When s i is inside the boundary layer, w ˙ i acts as a traditional integral. When s i is outside the boundary layer, the integral term w ˙ i is weakened. The extent to which it is weakened can be measured according to ( s i ε i ) ( 1 + τ i ) / k 0 i . It is clear that the extent of weakened demand satisfaction can be obtained through the regulating factor τ i , τ i [ 1 , k 0 i e i / ( s i ε i ) 1 ) . 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
s ˙ 1 s ˙ 2 s ˙ 3 = a ρ s + w 1 T χ 1 ( x ) + ν 1 ρ ˙ s , d + k 01 e 1 + k 11 w 1 a θ s + w 2 T χ 2 ( x ) + ν 2 θ ˙ s , d + k 02 e 2 + k 12 w 2 a z s + w 3 T χ 3 ( x ) + ν 3 z ˙ s , d + k 03 e 3 + k 13 w 3 + u 1 u 2 u 3
The designed multivariate super-twisting sliding mode approaching law is as follows:
s ˙ i = k 2 i s i 1 / 2 sign s i k 3 i s i + μ i μ ˙ i = k 4 i sign s i k 5 i s i
where k 2 i > 0 ; k 3 i > 0 ; k 4 i > 0 ; and k 5 i > 0 .
Substituting Equation (21) into Equation (20), we can obtain the following:
u 1 u 2 u 3 = a ρ s + w 1 T χ 1 ( x ) + ν 1 ρ ˙ s , d + k 01 e 1 + k 11 w 1 a θ s + w 2 T χ 2 ( x ) + ν 2 θ ˙ s , d + k 02 e 2 + k 12 w 2 a z s + w 3 T χ 3 ( x ) + ν 3 z ˙ s , d + k 03 e 3 + k 13 w 3 + k 21 s 1 1 / 2 sign s 1 k 31 s 1 + μ 1 k 22 s 2 1 / 2 sign s 2 k 32 s 2 + μ 2 k 23 s 3 1 / 2 sign s 3 k 33 s 3 + μ 3
The designed adaptive sliding mode control law for station-keeping control is as follows:
u 1 u 2 u 3 = a ρ s + w ^ 1 T χ 1 ( x ) ρ ˙ s , d + k 01 e 1 + k 11 w 1 a θ s + w ^ 2 T χ 2 ( x ) θ ˙ s , d + k 02 e 2 + k 12 w 2 a z s + w ^ 3 T χ 3 ( x ) z ˙ s , d + k 03 e 3 + k 13 w 3 + k 21 s 1 1 / 2 sign s 1 k 31 s 1 + μ 1 k 22 s 2 1 / 2 sign s 2 k 32 s 2 + μ 2 k 23 s 3 1 / 2 sign s 3 k 33 s 3 + μ 3

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  x ˙ = f x t , assuming the existence of a positive definite function  V s x R , this satisfies the following:
V ˙ s x κ 1 V s x κ 2 V s γ x
where κ 1 > 0 ; κ 2 > 0 ; γ ( 0 , 1 ) ; V 0 = V s x ( 0 ) ; and the nonlinear system x ˙ = f x t is finite-time stable.
The convergence time can be given by
t s 1 κ 1 1 γ ln κ 1 V 0 1 γ + κ 2 κ 2
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  k 0 i , k 1 i , k 2 i , k 3 i , k 4 i , and k 5 i  are appropriate, the control system can converge to the sliding surface in finite-time, and the orbital tracking error  e  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 e of the system converges asymptotically to zero.
Substituting Equation (23) into Equation (20), we can obtain the following:
s ˙ i = k 2 i s i 1 / 2 sign s i k 3 i s i + μ i + e a , i μ ˙ i = k 4 i sign s i k 5 i s i
where e a , i is the estimation error of the RBF neural network estimator, assuming that it satisfies e a , i ζ i s i 1 / 2 and e a = [ e a , 1 e a , 2 e a , 3 ] T .
Step 1. The Lyapunov function is chosen as follows:
V s i = 2 k 4 i s i + k 5 i s i 2 + 1 2 μ i 2 + 1 2 k 2 i s i 1 / 2 sign s i + k 2 i s i μ i 2 + 1 2 γ e a , i 2
Equation (26) can be written as V s i = η i T Ψ i η i + 1 2 γ e a , i 2 , where
η i = s i 1 / 2 sign s i s i μ i T , Ψ i = 1 2 4 k 4 i + k 2 i 2 k 2 i k 3 i k 2 i k 2 i k 3 i 2 k 5 i + k 3 i 2 k 3 i k 2 i k 3 i 2 .
V s i satisfies the following:
λ min Ψ i η i 2 2 V s i λ max Ψ i η i 2 2 + ζ i 2 2 γ s i
Since s i 1 / 2 η i 2 , it follows that
λ min Ψ i η i 2 2 V s i λ max Ψ i + ζ i 2 2 γ η i 2 2
Deriving the derivative for V s i , we can obtain the following:
V ˙ s i = 1 s i 1 / 2 η i T Θ i η i η i T Λ i η i + e a , i κ 1 , i η i + 1 s i 1 / 2 κ 2 , i η i + 1 γ e a , i e ˙ a , i = 1 s i 1 / 2 η i T Θ i η i η i T Λ i η i + e a , i κ 1 , i η i + 1 s i 1 / 2 κ 2 , i η i 1 γ w ^ ˙ i T χ i
where
Θ i = k 2 i 2 2 k 4 i + k 2 i 2 0 k 2 i 0 2 k 5 i + 5 k 3 i 2 3 k 3 i k 2 i 3 k 3 i 1
Λ i = k 3 i k 4 i + 2 k 2 i 2 0 0 0 k 5 i + k 3 i 2 k 3 i 0 k 3 i 1
κ 1 , i = 3 2 k 2 i k 3 i 2 k 5 i + k 3 i 2 k 3 i ,   κ 2 , i = 2 k 4 i + k 2 i 2 2 0 k 2 i 2
Let κ 1 , i η i + 1 s i 1 / 2 κ 2 , i η i 1 γ w ^ ˙ i T χ i = 0 ; then, the adaptive law can be designed as follows:
w ^ ˙ i = γ κ 1 i + κ 2 i s i 1 / 2 η i χ i T
where γ represents the adjustable coefficient associated with the adaptive law, γ > 0 ; (   ) represents the pseudo-converse operation. Thus,
V ˙ s i = 1 s i 1 / 2 η i T Θ i η i η i T Λ i η i
If Θ i and Λ i are positive definite symmetric matrices, then V ˙ s i is a negative definite. It is necessary to satisfy 4 k 4 i k 5 i > 8 k 4 i + 9 k 2 i 2 k 3 i 2 .
V ˙ s i 1 s i 1 / 2 λ min Θ i η i 2 2 λ min Λ i η i 2 2
Using Equation (28), it can be derived that
V ˙ s i γ 1 i V s i 1 / 2 γ 2 i V s i
where
γ 1 i = λ min 1 / 2 Ψ i λ min Θ i λ max Ψ i + ζ i 2 2 γ ,   γ 2 i = λ min Λ i λ max Ψ i + ζ i 2 2 γ
From Lemma 1, if the values of the parameters such as k 0 i , k 1 i , k 2 i , k 3 i , k 4 i , and k 5 i are appropriate, the control system can converge to the sliding surface in finite-time.
Step 2. The system reaches the sliding surface, s i = 0 , satisfying s i ε i ; then,
e ˙ i + k 0 i e i + k 1 i e i d τ = 0
When s ˙ i tends to zero, there exists the following:
e ¨ i + k 0 i e ˙ i + k 1 i e i = 0
The Lyapunov function is chosen as follows:
V e = 1 2 e ˙ T e ˙ + 1 2 i = 1 3 k 1 i e i 2 = 1 2 i = 1 3 e ˙ i 2 + 1 2 i = 1 3 k 1 i e i 2
The derivation of Equation (36) yields the following:
V ˙ e = i = 1 3 e ˙ i e ¨ i + k 1 i e i = i = 1 3 e ˙ i e ¨ i + k 1 i e i = i = 1 3 e ˙ i k 0 i e ˙ i = i = 1 3 k 0 i e ˙ i 2 0
V ˙ e 0 is equivalent to e ˙ i 0 . According to LaSalle’s invariance principle, it can be known that ( e i = 0 , e ˙ i = 0 ) is the global asymptotically stable equilibrium point of Equation (34), i.e., lim t e i t = lim t e ˙ i t = 0   ( i = 1 , 2 , 3 ) . The orbital tracking error e 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 α s , turning angle δ s , and solar electric propulsion acceleration a e p after the following is obtained: a e p = a e p , x 2 + a e p , y 2 + a e p , z 2 1 / 2 . From Equation (12), the solar electric propulsion acceleration components can be derived as follows:
a e p , x = u 1 cos θ s u 2 ρ s sin θ s + β s sin α s cos 2 α s sin δ s ρ s r s 2 sin θ s β s cos 3 α s r s 3 ρ s + sin α s cos 2 α s cos δ s r s 3 z s cos θ s a e p , y = u 1 sin θ s + u 2 ρ s cos θ s β s sin α s cos 2 α s sin δ s ρ s r s 2 cos θ s β s cos 3 α s r s 3 ρ s + sin α s cos 2 α s cos δ s r s 3 z s sin θ s a e p , z = u 3 β s cos 3 α s r s 3 z s β s sin α s cos 2 α s cos δ s r s 3 ρ s
where ρ s , θ s , z s , and r s are real-time measurement values; the light pressure factor β s can be found via Equation (39).
β s = 1 + z s , d ρ s , d 2 1 / 2 z s , d ρ s , d 2 + 1 ω s , d 2 ρ s , d 2 + z s , d 2 3 / 2 2 3 / 2 z s , d ρ s , d 2 + 1 ω s , d 2 ρ s , d 2 + z s , d 2 3 / 2 2
where ρ s , d , z s , d , and ω s , d 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:
m ˙ = T I s p g 0
where T is the electric thrust, T = [ T x T y T z ] T , and T = m a e p ; I s p is the specific impulse of electric propulsion; and g 0 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:
α , δ = arg min α s [ π / 2 , π / 2 ] δ s [ π , π ] a e p
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 T can be obtained.

4. Simulation Results and Discussion

4.1. Simulation Conditions

The studied nominal heliocentric displaced orbit is located above the ecliptic plane. The nominal orbit always has the same phase as Earth, with a period of 1 year. The dimensionless parameters of the heliocentric displaced orbit are shown in Table 1. From Equation (39), the light pressure factor can be calculated as 0.7356, which is less than 1 and can be realized for engineering applications.
It is assumed that the initial mass of the solar sail spacecraft is 500 kg. The combination of a slider and RSB was used as the propellant-free actuator, and the argon Hall thruster was used as the solar electric thruster. The relevant parameters of the argon Hall thruster are shown in Table 2.
The relevant parameters of the adaptive sliding mode controller for station-keeping are shown in Table 3.
We assumed that the orbital radius error and the displaced height error of the solar sail spacecraft in the heliocentric displaced orbit at the initial injection were 10 , 000   km and 10 , 000   km , respectively. The simulation step was set to 0.000001 TU (1 TU ≈ 58 d) and the simulation time was set to 365 d.
Referring to the simulation conditions in the literature [29,30], the dimensionless external disturbance acceleration a d and the internal unmodeled dynamics a Δ f were set as follows:
a d 1 = 0.00001 sin 80 t + 100 a d 2 = 0.00001 cos 60 t 100 a d 3 = 0.00001 sin 80 t + 1 ,   a Δ f 1 = 0.01 ρ s cos 5 θ s 0.1 / r e 3 a Δ f 2 = 0.01 sin 2 θ s / ρ s r e 3 a Δ f 3 = 0.01 z s sin 3 θ s + 0.1 / r e 3
where r e denotes the position vector of the geocenter with respect to the solar sail spacecraft.

4.2. Simulation Results and Discussion under Different Sources of Disturbance

Simulations under different sources of perturbations will be studied in the following sections. To further verify the effectiveness of the adaptive sliding mode station-keeping control strategy proposed in this paper, the orbit position tracking error results were compared with the literature [28,29,30].

4.2.1. When the Solar Sail Spacecraft Is Mainly Subjected to External Disturbances

The simulation results of the station-keeping control are shown in Figure 3, Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8.
Figure 3 shows the response curve of the improved integral sliding surface. The sliding surface eventually converges to near zero. Figure 4 describes the trajectory of the solar sail spacecraft under the action of the proposed high-performance station-keeping control strategy. The black solid line represents the reference orbit and the yellow dashed line represents the actual trajectory. It can be seen from Figure 4 that the two overlap well, indicating that the station-keeping control strategy proposed in this paper has good accuracy. Figure 5 shows the response curves of orbital position errors and angular velocity. The orbital radius error and displaced height error converge to within 0.8 m after 20 d, and the orbital angular velocity converges to the nominal angular velocity after 8 days.
Figure 6 shows the estimation results of the RBF neural network estimator for the external disturbance acceleration. It can be seen that, during the spacecraft’s approach to the nominal orbit, the estimator’s performance in estimating disturbance acceleration is not ideal. However, when the orbital angular velocity converges to the nominal angular velocity, the estimation of a d 2 is enhanced. When the orbital position tends to be stable, the estimator’s estimate of each component of the disturbance acceleration can better approximate the actual value.
Figure 7 shows the results of the control variables in the case of minimum propellant consumption. From Figure 7, it can be seen that the cone angle and turning angle stabilize to 20.464 ° and 0.001 ° , respectively, after 15 days, and the propulsion force of the argon Hall electric thruster stabilizes within 0.02 mN after 15 days.
Figure 8 describes the mass variation in the solar sail spacecraft under different propulsion methods. M1 represents the mass variation when a hybrid propulsion method is used; M2 represents the mass variation when the argon Hall electric thruster is operating at maximum power; and M3 represents the mass variation when only electric thrusters are used for propulsion. Throughout one year, M1 decreases by 0.46 kg and M2 decreases by 218.67 kg, while the amount of change in M2 is much smaller than M3, which means that a single electric thruster cannot accomplish the station-keeping control mission. Therefore, hybrid propulsion is currently the most efficient way to achieve the station-keeping control of a solar sail spacecraft.

4.2.2. When the Solar Sail Spacecraft Is Subjected to Both External Disturbances and Internal Unmodeled Dynamic

Figure 9 shows the response curves of orbital position errors and angular velocity. Similar to the results in Section 4.2.1, the orbital radius error and displaced height error converge to within 0.8 m after 20 d, and the orbital angular velocity converges to the nominal angular velocity after 8 days.
Figure 10 shows the estimation results of the RBF neural network estimator for the combined disturbance acceleration. Similar to the results in Section 4.2.1, the estimation in the transition stage is unsatisfactory, while in the steady state, the estimated value is able to approximate the actual value well.
Figure 11 shows the results of the control variables in the case of minimum propellant consumption. As can be seen from Figure 11, the cone and turning angles oscillate and change during the one-year running period, and, in particular, the turning angle has a large range of variation, which means that the orientation of the sail surface is always changing. From a comparison with Section 4.2.1, it was found that the unmodeled dynamics are the main reason for the obvious change in attitude angles in Figure 11. This is due to the relatively large order of magnitude of the unmodeled dynamics. The thrust of the argon Hall thruster changes correspondingly with the change in attitude angles.
Figure 12 describes the mass variation in a solar sail spacecraft. After one year, M1 decreases by 26.69 kg. Considering argon prices, mass requirements for solar sail spacecraft, and launch costs, it can be found that hybrid propulsion has more obvious advantages for long-period missions such as heliocentric displaced orbits. In addition to this method of minimizing propellant consumption, a single electric thruster can be used to achieve station-keeping control under a given attitude angle. Of course, this method requires higher mass of the propellant to be consumed.

4.2.3. Comparison of Simulation Results

Since the simulation conditions are consistent, the results of the comparison between the proposed station-keeping control strategy and the other three control strategies in terms of orbit position errors will be given directly here, as shown in Table 4. As can be seen from Table 4, the adaptive sliding mode controller proposed in this paper has faster station tracking speed and higher control accuracy compared with existing controllers.
From the above simulation results, it is clear that the proposed station-keeping control strategy has the following characteristics: (1) the tracking position errors can reach the order of decimeters, close to zero; (2) the steady-state errors of orbit tracking under different sources of disturbances are in the same order of magnitude, and are nearly the same size; and (3) the values of the relevant parameters of the station-keeping controller under different sources of disturbances are consistent.

5. Conclusions

In this paper, an adaptive high-performance sliding mode control strategy was designed to solve the station-keeping control problem of the heliocentric displaced orbit. An RBF neural network estimator was designed, which achieved online estimation and compensation of the combined disturbance acceleration of each channel; an adaptive orbit-keeping controller based on the improved integral sliding surface and multivariate super-twisting sliding mode convergence law was designed, which achieved fast tracking of the nominal orbit with high control accuracy. The simulation results showed that the proposed station-keeping control strategy has strong robustness and a fast orbit tracking speed, with position tracking errors up to decimeters.

Author Contributions

Conceptualization, T.Z. and R.M.; methodology, T.Z.; software, T.Z.; validation, T.Z., Y.Z. and Z.L.; formal analysis, T.Z.; investigation, T.Z. and Y.Z.; resources, T.Z., Z.Z., B.L. and C.Y.; data curation, T.Z. and Z.L.; writing—original draft preparation, T.Z.; writing—review and editing, R.M., Y.Z., Z.L., Z.Z., B.L. and C.Y.; visualization, T.Z.; supervision, R.M., Z.Z., B.L. and C.Y.; project administration, R.M.; funding acquisition, R.M., Z.Z., B.L. and C.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Shanghai Institute of Satellite Engineering (SISE) and by the pre-research on civilian space technology, grant number D010301. The APC was funded by the Institute of Aircraft Systems Engineering.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author.

Acknowledgments

We would like to thank SISE for funding; the relevant engineering data that SISE provided assisted us significantly.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yang, J. Thinking on the Development of China’s Cause. Important Issues 2017, 11, 1–2. [Google Scholar] [CrossRef]
  2. Yin, J.; Zhang, R.; Zhang, X. Orbital Dynamic Characters and Orbital Maintenances of Pole-sitter Spacecraft. Spacecr. Eng. 2018, 27, 26–34. [Google Scholar]
  3. Ullery, D.C.; Soleymani, S.; Heaton, A.; Orphee, J.; Johnson, L.; Sood, R.; Kung, P.; Kim, S.M. Strong Solar Radiation Forces from Anomalously Reflecting Metasurfaces for Solar Sail Attitude Control. Sci. Rep. 2018, 8, 10026. [Google Scholar] [CrossRef] [PubMed]
  4. Baoyin, H.; McInnes, C.R. Solar sail orbits at artificial Sun-Earth libration points. J. Guid. Control Dyn. 2005, 28, 1328–1331. [Google Scholar] [CrossRef]
  5. Zhang, J.; Cao, X.; Liu, W.; She, W.; Wei, J. Orbit keeping and control of solar sail spacecraft with RCD on artificial libration points. Flight Dyn. 2021, 39, 77–81. [Google Scholar]
  6. Bianchi, C.; Niccolai, L.; Mengali, G.; Quarta, A.A. Collinear artificial equilibrium point maintenance with a wrinkled solar sail. Aerosp. Sci. Technol. 2021, 119, 107150. [Google Scholar] [CrossRef]
  7. Simo, J.; McInnes, C.R. Displaced solar sail orbits: Dynamics and applications. Adv. Astronaut. Sci. 2010, 136, 1803–1816. [Google Scholar]
  8. Ceriotti, M.; McInnes, C.R. Hybrid solar sail and solar electric propulsion for novel Earth observation missions. Acta Astronaut. 2011, 69, 809–821. [Google Scholar] [CrossRef] [Green Version]
  9. Li, Z.; Yuan, C.; Yu, H.; Zuo, C. Station-Keeping Control of Solar Sail Spacecraft Planet-Centred Displaced Orbit. J. Spacecr. Rockets 2019, 45, 33–37. [Google Scholar]
  10. Niccolai, L.; Mengali, G.; Quarta, A.A.; Caruso, A. Feedback control law of solar sail with variable surface reflectivity at Sun-Earth collinear equilibrium points. Aerosp. Sci. Technol. 2020, 106, 106144. [Google Scholar] [CrossRef]
  11. Moore, L.; Ceriotti, M.; McInnes, C.R. Station-keeping for a solar sail during lander/probe deployment using feedback control. Acta Astronaut. 2022, 201, 182–197. [Google Scholar] [CrossRef]
  12. Moore, L.; Ceriotti, M. Solar sails for perturbation relief: Application to asteroids. Adv. Space Res. 2021, 67, 3027–3044. [Google Scholar] [CrossRef]
  13. Gong, S.; Li, J.; Simo, J. Orbital motions of a solar sail around the L2 Earth–Moon libration point. J. Guid. Control Dyn. 2014, 37, 1349–1356. [Google Scholar] [CrossRef] [Green Version]
  14. Yuan, J.; Gao, C.; Zhang, J. Periodic orbits of solar sail equipped with reflectance control device in Earth–Moon system. Astrophys. Space Sci. 2018, 363, 23. [Google Scholar] [CrossRef]
  15. Zhu, M. Dynamics and Control of Solar Sail Spacecraft. Ph.D. Thesis, University of Science and Technology of China, Heifei, China, 2016. [Google Scholar]
  16. Lou, Z.; Zhang, K.; Wang, Y.; Gao, Q. Active disturbance rejection station-keeping control for solar-sail libration-point orbits. J. Guid. Control Dyn. 2016, 39, 1913–1917. [Google Scholar] [CrossRef]
  17. Lou, Z. Dynamics and Control of Continuous Low-Thrust Spacecraft on Libration Points. Ph.D. Thesis, University of Science and Technology of China, Heifei, China, 2019. [Google Scholar]
  18. Huang, J.; Biggs, J.D.; Bai, Y.; Cui, N. Integrated guidance and control for solar sail station-keeping with optical degradation. Adv. Space Res. 2021, 67, 2823–2833. [Google Scholar] [CrossRef]
  19. Huang, J.; Bai, Y.; Biggs, J.D.; Cui, N. Indirect Estimation Based Active Disturbance Rejection Control for Solar Sail Station-Keeping. Int. J. Aeronaut. Space Sci. 2021, 22, 1493–1503. [Google Scholar] [CrossRef]
  20. Heiligers, J.; McInnes, C.R.; Biggs, J.D.; Ceriotti, M. Displaced geostationary orbits using hybrid low-hrust propulsion. Acta Astronaut. 2012, 71, 51–67. [Google Scholar] [CrossRef] [Green Version]
  21. Simo, J.; McInnes, C.R. Designing displaced lunar orbits using low-thrust propulsion. J. Guid. Control Dyn. 2010, 33, 259–265. [Google Scholar] [CrossRef] [Green Version]
  22. Simo, J.; McInnes, C.R. Feedback stabilization of displaced periodic orbits: Application to binary asteroids. Acta Astronaut. 2014, 96, 106–115. [Google Scholar] [CrossRef] [Green Version]
  23. Heiligers, J.; Ceriotti, M.; McInnes, C.R.; Biggs, J.D. Displaced geostationary orbit design using hybrid sail propulsion. J. Guid. Control Dyn. 2012, 34, 1852–1866. [Google Scholar] [CrossRef]
  24. SpaceX just Announced Their Hall Thruster Parameters, actually Ar Propellant. Available online: https://m.163.com/dy/article/HUMMT1QL05119RIN.html (accessed on 22 April 2023).
  25. McKay, R.; Macdonald, M.; Biggs, J.; McInnes, C. Survey of highly non-Keplerian orbits with low-thrust propulsion. J. Guid. Control Dyn. 2011, 34, 645–666. [Google Scholar] [CrossRef] [Green Version]
  26. McInnes, C.R. Passive control of displaced solar sail orbits. J. Guid. Control Dyn. 1998, 21, 975–982. [Google Scholar] [CrossRef]
  27. Bookless, J.; McInnes, C. Dynamics and control of displaced periodic orbits using solar-sail propulsion. J. Guid. Control Dyn. 2006, 29, 527–537. [Google Scholar] [CrossRef]
  28. Qian, H.; Zheng, J.; Yu, X.; Gao, L. Dynamics and control of displaced orbits for solar sail spacecraft. Chin. J. Space Sci. 2013, 33, 458–464. [Google Scholar] [CrossRef]
  29. Zhang, K.; Lou, Z.; Wang, Y.; Chen, S. Station-keeping control of spacecraft using hybrid low-thrust propulsion in heliocentric displaced orbits. Acta Aeronaut. Astronaut. Sin. 2015, 36, 3910–3918. [Google Scholar]
  30. Chen, Y.; Qi, R.; Zhang, J.; Wang, H. High-performance sliding mode control for orbit keeping of spacecraft using hybrid low-thrust propulsion. Acta Aeronaut. Astronaut. Sin. 2019, 40, 322827. [Google Scholar]
  31. Yu, S.; Yu, X.; Shirinzadeh, B.; Man, Z. Continuous finite-time control for robotic manipulators with terminal sliding mode. Automatica 2005, 41, 1957–1964. [Google Scholar] [CrossRef]
Figure 1. Displaced orbit and orbital coordinate system of the solar sail spacecraft.
Figure 1. Displaced orbit and orbital coordinate system of the solar sail spacecraft.
Symmetry 15 01549 g001
Figure 2. Block diagram of the station-keeping control system for a solar sail spacecraft.
Figure 2. Block diagram of the station-keeping control system for a solar sail spacecraft.
Symmetry 15 01549 g002
Figure 3. Change curves of improved integral sliding surfaces.
Figure 3. Change curves of improved integral sliding surfaces.
Symmetry 15 01549 g003
Figure 4. Reference orbit and controlled trajectory.
Figure 4. Reference orbit and controlled trajectory.
Symmetry 15 01549 g004
Figure 5. Change curves of position errors and angular velocity: (a) change curves of orbital radius error and displaced height error; (b) change curve of angular velocity.
Figure 5. Change curves of position errors and angular velocity: (a) change curves of orbital radius error and displaced height error; (b) change curve of angular velocity.
Symmetry 15 01549 g005
Figure 6. Change curves of external disturbances and their estimations.
Figure 6. Change curves of external disturbances and their estimations.
Symmetry 15 01549 g006
Figure 7. Station-keeping control variables: (a) change curves of attitude angles; (b) change curves of propulsion force components of argon Hall thruster.
Figure 7. Station-keeping control variables: (a) change curves of attitude angles; (b) change curves of propulsion force components of argon Hall thruster.
Symmetry 15 01549 g007
Figure 8. Mass variation in a spacecraft under different propulsion methods.
Figure 8. Mass variation in a spacecraft under different propulsion methods.
Symmetry 15 01549 g008
Figure 9. Change curves of position errors and angular velocity: (a) change curves of orbital radius error and displaced height error; (b) change curve of angular velocity.
Figure 9. Change curves of position errors and angular velocity: (a) change curves of orbital radius error and displaced height error; (b) change curve of angular velocity.
Symmetry 15 01549 g009
Figure 10. Change curves of external disturbances and their estimations.
Figure 10. Change curves of external disturbances and their estimations.
Symmetry 15 01549 g010
Figure 11. Station-keeping control variables: (a) change curves of attitude angles; (b) change curves of propulsion force components of argon Hall thruster.
Figure 11. Station-keeping control variables: (a) change curves of attitude angles; (b) change curves of propulsion force components of argon Hall thruster.
Symmetry 15 01549 g011
Figure 12. Mass variation in a spacecraft.
Figure 12. Mass variation in a spacecraft.
Symmetry 15 01549 g012
Table 1. Parameters of nominal heliocentric displaced orbit.
Table 1. Parameters of nominal heliocentric displaced orbit.
Orbital RadiusDisplaced HeightAngular Velocity
0.7AU0.4AU1
Table 2. Parameters of argon Hall thruster.
Table 2. Parameters of argon Hall thruster.
PropellantThrust (mN)Specific Impulse (s)Total EfficiencyPower (kW)Mass (kg)
Argon170250050%4.22.1
Table 3. Parameters of adaptive sliding mode controller.
Table 3. Parameters of adaptive sliding mode controller.
Parameters of the Sliding SurfaceParameters of the
Approaching Law
Parameters of the Adaptive Law
k 0 i = 165 ; k 1 i = 0.0003 ;
ε i = 0.00001 ; τ i = 0.98
k 2 i = 1 ; k 4 i = 0.001 γ = 0.001
Table 4. Comparison of orbit tracking position errors.
Table 4. Comparison of orbit tracking position errors.
Maximum Thrust of Electric Thruster (mN)Convergence Time (d)Orbital Radius Tracking Error (km)Displaced Height Tracking Error (km)
In this article17020-0.00060.0008
In the literature [28]23040-0.00120.0016
In the literature [29]23090Within 0.05Within 0.05
In the literature [30]230730-30.0006-198.9652
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, T.; Mu, R.; Zhou, Y.; Liao, Z.; Zhang, Z.; Liao, B.; Yao, C. A Station-Keeping Control Strategy for a Symmetrical Spacecraft Utilizing Hybrid Low-Thrust Propulsion in the Heliocentric Displaced Orbit. Symmetry 2023, 15, 1549. https://doi.org/10.3390/sym15081549

AMA Style

Zhang T, Mu R, Zhou Y, Liao Z, Zhang Z, Liao B, Yao C. A Station-Keeping Control Strategy for a Symmetrical Spacecraft Utilizing Hybrid Low-Thrust Propulsion in the Heliocentric Displaced Orbit. Symmetry. 2023; 15(8):1549. https://doi.org/10.3390/sym15081549

Chicago/Turabian Style

Zhang, Tengfei, Rongjun Mu, Yilin Zhou, Zizheng Liao, Zhewei Zhang, Bo Liao, and Chuang Yao. 2023. "A Station-Keeping Control Strategy for a Symmetrical Spacecraft Utilizing Hybrid Low-Thrust Propulsion in the Heliocentric Displaced Orbit" Symmetry 15, no. 8: 1549. https://doi.org/10.3390/sym15081549

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop