Abstract

The drag-free satellites, being space-borne ultrahigh precise measurement platforms, have played irreplaceable roles in a great number of space science missions such as navigation, earth science, fundamental physics, and astrophysics. Most of these missions have to be performed based on the satellites placed with double cube test-masses, which makes the satellite layout and control strategy be more complex. This paper investigates the orbit keeping control problem of a class of low Earth orbit drag-free satellites with double cube test masses. A disturbance observer-based composite control method is proposed, which consists of an extended sliding mode observer and the tube-based robust model predictive control approach. In this design, the observer is proposed to estimate the relative position and velocity of the satellite and the external space disturbance force. A tube-based robust model predictive control scheme is then developed to stabilize the satellite orbit control systems in the presence of actuator saturation, state constraints, and additive stochastic noises. Finally, a simulation example is presented to demonstrate the efficacy and superiority of the proposed orbit control method.

1. Introduction

In recent years, the drag-free satellites [1], being space-borne ultrahigh precise measurement platforms, have played irreplaceable roles in many space science missions, such as the test of equivalence principle [2], the measurement of the Earth gravity field [3], and the detection of gravitational waves [4]. The drag-free satellites possess many advantages; for example, they can provide autonomous precision orbit determination, map the static and time-varying components of the Earth’s mass distribution more accurately, deepen the understanding of the fundamental force of gravity, eventually open up a new window to the universe through the detection and observation of gravitational waves, and so forth.

The key technology of the drag-free satellite is the gravitational reference sensor (GRS), which insulates an internal free-floating test mass (TM, also called proof mass) from both external disturbances and disturbances caused by the spacecraft itself [5]. The drag-free satellites can be divided into two types [6]. The first one is the “accelerometer” drag-free mode, where an electrostatic accelerometer is used as the primary sensor and an electrostatic suspension actuator is paired to maintain the TM to be centred in its cage; therefore it can counter the disturbance forces acting on the spacecraft [711]. The second one is free-falling TM mode, in which the satellite provides indirect drag-free behaviour by tracking the movement of the free-falling TM in the cage [1215]. In particular, the structure of a satellite with two cube TMs is always regarded as the primary layout in these missions, which makes the GRS be more complex and the control system design work be more challenging. For example, as shown in Figure 1, the drag-free satellite containing two cube TMs will enter a low Earth orbit (LEO), which decays more rapidly due to the decelerating effects of the Earth’s atmosphere.

In the research area of control for drag-free satellites, a trade-off should be considered among system complexity, fuel conservation, cost of components, operations, performance, and so forth [16]. In fact, it is always difficult for the drag-free satellite to obtain control satisfactory performance due to the existence of external time-varying environmental disturbances and the unmodelled internal uncertainties. To solve this problem, various types of robust control algorithms have been proposed, such as model predictive control (MPC) [17], control [18, 19], -synthesis control [20], embedded model control (EMC) technique [21], and quantitative feedback theory [22]. However, these design results can only realize disturbance attenuation but not disturbance rejection. In order to improve closed-loop control performance, the disturbance observer-based control [2331] and active disturbance rejection control [3239] have been developed for drag-free satellite to pursue more ideal control system performance [25, 4050]. However, how to design a suitable observer to estimate the disturbances more precisely is still an open and attractive problem [51].

This paper investigates the orbit control problem for a class of LEO drag-free satellites with double cube TMs based on the relative position dynamics with state constraints, actuator saturation, and the additive stochastic disturbances. In this design, an extended sliding mode observer method is developed to estimate the system states. Then, a tube-based robust model predictive control (TRMPC) approach is presented [52] to stabilize the resulting orbit control systems. In particular, the TRMPC design approach is divided into the following two steps: (i) an offline evaluation of the constraints to ensure the uncertain future trajectories to lie in sequence of sets, which is called tubes, and (ii) an online MPC scheme designed for the nominal trajectories. As a result, the orbit control issue for the drag-free satellite is solved.

The remaining of this paper is organized as follows. Section 2 formulates the design problem under investigation; the design of disturbance observation-based composite control strategy is presented in Section 3. A demonstrated simulation example is provided in Section 4, and the paper is concluded in Section 5.

2. Problem Formulation

According to [22, 53], for a LEO satellite containing two TMs, as shown in Figure 2, the drag-free control strategy is defined as follows: TM1 is chosen as the gravitational reference, which flies freely in a pure gravitational orbit, and the satellite and TM2 are controlled to follow TM1; then the linearized relative position dynamics between the TMs and the satellite are given as follows:where means the satellite and means the TM2; , , and represent the relative position variables, the mass of the satellite or the TM2, and the control forces; is the sum of all kinds of disturbances, which for the satellite, the air drag, the solar radiation pressure, and the thruster quantization error and other stochastic noises is considered primarily, and, for the TM2, it mainly consists of electrostatic interference signals, actuator quantization error, and other stochastic noises [9, 40].where is orbit angular velocity of the TM1 and and are the damping factors between the TM and its condenser cage [54].

The desired positions and velocities of the satellite and TM2 are denoted as , , , and , respectively. The desired positions and are some constants dependent on the layout of the satellite. The relative velocities are defined as follows: and . Define the following error variables for case of presentation:

By substituting (3) into (1), one can obtain

Due to the limit volume of the capacitor cage, the relative error variables should satisfy and , which are the state constraints. For convenience, formations (4) should be translated into a state space model. By choosing and , , , , and , where and denote the unknown stochastic noises and unknown estimable disturbances, respectively. Then, the state space model is given as follows:where

Before designing the disturbance observer-based composite control strategy, the following assumptions are given about (t) and .

Assumption 1. The unknown external disturbance satisfies the following: and , where and are known real constants.

Assumption 2. The unknown stochastic noise is assumed to be discontinuous but bounded subject to and belongs to a bounded and convex subset containing the origin in its interior.

3. Design of Disturbance Observer-Based Composite Control

The control system composition of a LEO satellite with two cube TMs is shown in Figure 3. To achieve the “drag-free” goal and eliminate the effects of external disturbances (the environment disturbances may contain the atmosphere and the solar radiation pressure, the thruster quantization error, the electrostatic noises, the actuator quantization error, etc.), a disturbance observer-based composite control approach is proposed, which consists of an extended sliding mode observer and a tube-based robust model predictive control. The detailed design process is as follows. First, for system (5), the observer is proposed to estimate states for state feedback control and for active disturbance rejection control . Second, a tube-based robust model predictive control is adopted, which can cope with state constraints and actuator saturation and attenuate the effects of additive stochastic noises. For convenience, in the following discussion, the subscript in (5) is omitted.

3.1. Design of Extended Sliding Mode Observers

Motivated by the augmented strategy method in [34, 36, 55], we define the following extended vectors and matrices:

For system (7), consider the following continuous-time extended SMO [55]:where and is the discontinuous term designed as follows:where and , , and and are the maximal and minimal nonzero eigenvalues of matrix , respectively.

Remark 1. It is worth mentioning that the main function of the proposed extended SMO is to estimate the state and disturbances vectors simultaneously, which is a design basis of the subsequent composite control law.
Define ; then the error system is derived as follows:

Lemma 1 (see [56]). System (5) has the relative degree with respect to the unknown input (i.e., the system is strongly observable).

Lemma 2. If the matrix pair is detectable and the conditionholds, then the pair is detectable; that is, there exists an observer gain matrix such that is Hurwitz.

Proof. It can be derived that

Case 1. When , for the matrix in right side of equation (12), we havewhich means thatThat meansHence, if the pair is detectable, we haveThen, the following rank condition holds:

Case 2. When , if rank condition (11) holds, we haveFinally, for both cases, we can always imply , , , which implies that the pair is detectable.
To attenuate the influence of for error system (10), we define the prescribed performance index as follows:where is the weighted estimation error and is the weight matrix. The following main theorem provides the existence condition of the proposed observer (8).

Theorem 1. The solution of observer error system (10) is asymptotically stable with the prescribed performance index if there exist positive and definite matrix and matrix with appropriate dimensions, in minimizing , such that the following matrix constraints hold:Moreover, the observer gain is designed as .

Proof. For system (10), select the Lyapunov function as with being designed; it is easy to obtain thatIn the following analysis, let ; we consider two cases separately.

Case 3. ; in this case, we havewhere can be chosen such that .

Case 4. ; in this case, it can be derived thatBased on Cases 3 and 4, we can conclude thatTo minimize the effect of the disturbance on the estimation error in the sense of norm, we consider the following constraint:In light of Schur complement, it is easy to derive thatAccordingly, it can be seen that is ensured provided that holds. Therefore, the solution of error system (10) is asymptotically stable with the prescribed performance index as .

Remark 2. It should be pointed out that, in Theorem 1, the equality constraint can be rewritten asHence, we can introduce the following condition:where is a parameter to be designed. Then the design problem of observer gains L can be converted into the following minimization problem:

3.2. Design of Tube-Based Robust Model Predictive Control

For equations (5), is estimated and eliminated through the extended sliding mode observer; thus we have to compensate it by disturbance compensation feed-forward control. Translate formulations (5) without into the discrete model of the form in (30) as follows:

As mentioned before, the controller is designed as , whereis the feed-forward control part to compensate and is the model predictive control law to be designed. Define the estimation error as

Then, system (30) becomes

As proved in the observer design results in Theorem 1, we have with ; thus holds, where is a small constant and when .

Define ; it satisfies , which belongs to the bounded and convex subset . Hence, system (30) can be written in the following form:

The investigated system (30) is subject to hard constraints on both state and input vectors with the following form:where and are polytopes.

To solve this control problem, a robust MPC algorithm is considered [57] by repeatedly solving an optimal control problem, where the finite horizon quadratic cost to be minimized at the current time is

In (36), is the MPC prediction horizon, , , , and is the solution of the algebraic Riccati equation [57].

Due to the presence of the unknown disturbance , we rewrite the state vector of the system as the sum of a nominal part and an error part in the following form:where denotes the deviation of the real state with respect to the nominal one.

Design the following feedback policy (39) for system (30):where denotes the nominal input vector and the gain matrix should be selected such that is Schur-stable; then the corresponding nominal and error dynamics can be described, respectively, as follows:

Hence, the finite horizon optimal quadratic cost (36) can be redefined in terms of nominal state and control input asand the finite horizon optimal control problem can be reformulated as follows.

Definition 1. Given the nominal system dynamics (40), cost (42), and nominal constraints set , the nominal robust MPC finite horizon optimization problem can be described asThe solution of (43) is the optimal nominal control sequence and the first control action, that is, , represents the optimal control to be applied to system (40).
The proposed control law applied on the uncertain system (40), according to the control policy adopted, isThe composite closed-loop system then satisfiesFor the TRMPC approach, the matrix in the control policy (39) is designed to stabilize system (30). Consider the following closed-loop system:Hence, the satisfaction of the following condition aims to define the feedback gain that stabilizes the system.In order to robustly satisfy the mission constrains, they are tightened to allow the trajectories of the uncertain system, affected by disturbance, to lie in a tube centred on the nominal one, where each trajectory is related to a particular realization of the uncertainty at each time step . The derivations of the nominal state, input, and terminal constrains set are described according to the approach proposed in [57], such that the constraints in (39) of system (30) are satisfied for every realization of the disturbance sequence by suitable design of the tube.
We now define , where is the Minkowski sum and is the set multiplication; the uncertain set of the error is the minimal robust positive invariant set forThen the state and control input vector constraints in (35) are satisfied provided thatwhere denotes the Pontryagin set difference. It is obvious that the terminal constraint for system (30) at time instant is ensured if the normal system (40) satisfies the tighter constraintMoreover, these assertions only make sense if the disturbance set is sufficiently small to satisfy the following Assumption 3, as defined in [57].

Assumption 3. (Restricted disturbances for constraints satisfaction) and .
The next step is to define a robust positively invariant set for [19] to obtain the tighter constraints acting on the nominal system. Then the constraints are considered for the TRMPC problem. Once the uncertainty set is evaluated, an inner approximation of the nominal constraint set can be constructed. In this design, we adopt the following strategy presented in [57].

Algorithm 1. Computation of and (1)Define the linear state constraint as: ;(2)Construct the nominal state constraint inequality ;(3)Approximate the upper value of as ;(4)Choose a suitable and N such that , compute ;(5)Compute the nominal state constraint set ;(6)Compute the nominal control constraint set , where .Hence, the observer-based model predictive control strategy could be formally described by the following algorithm.

Algorithm 2. Disturbance observer-based model predictive control strategy.(1)Initialization: at time , set where denotes the current state.(2)At time , considering the current state (), based on the disturbance estimation from the observer (8), solve the nominal optimal control problem (43) to obtain the nominal control vector and the control input vector .(3)If the nominal optimal control problem (43) is infeasible, adopt safety/recovery procedure.(4)Apply the control to the system (45) and (46);(5)Calculate the estimation from the observer (8) as successor state of the system (30), and calculate the successor state of the nominal system (40).(6)Set () = (), set , and go to (2).

4. Simulation Results

In this section, a numerical simulation is carried out to verify the effectiveness of the extended sliding mode-based TRMPC approach. Suppose that the satellite is flying at a little eccentric low Earth orbit with altitude of 300 km. The parameters in (1) are given in Table 1 [54].

For the relative position motion between the TM1 and the satellite, define , where

The resolution and the maximum value of the thruster are  N and 0.4 N. Besides, a zero mean white noise with mean squared error being  N is added as the stochastic disturbance.

For the relative position motion between TM1 and TM2, define , where

The resolution and the maximum value of the thruster are  N and  N. Similarly, a zero mean white noise with mean squared error being  N is added as the stochastic disturbance.

The simulation results are shown in Figures 410 as follows. Define the estimated error (i = 1, 2), as shown in Figures 4 and 5, and can be estimated precisely by the ESMO. From Figures 611, the relative motion variables , , , and could be estimated well by the ESMO. In addition, the stable control accuracy of the composite control approach has achieved or , which shows that the developed extended sliding mode observer method and tube-based model predictive control law are effective.

5. Conclusions

This paper has considered the relative position control of the drag-free satellite with double cube test-masses in the presence of external disturbance, additive stochastic disturbances, actuator quantization error, actuator saturation, and state constraints. An extended sliding mode observer method is adopted to estimate the state vector and external disturbance, based on which a tube-based robust model predictive control scheme is developed. The designed control method can not only cope with the constraints of control and state but also attenuate the effect of additive stochastic noises. Future work will be focused on the consideration of relative attitude dynamics between the test-masses and the satellite.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (61833009, 11972130, and 61690212) and Guangdong Major Project of Basic and Applied Basic Research (Grant no. 2019B030302001).