Active Disturbance Rejection Station-Keeping Control of Unstable Orbits around Collinear Libration Points

An active disturbance rejection station-keeping control scheme is derived and analyzed for station-keeping missions of spacecraft along a class of unstable periodic orbits near collinear libration points of the Sun-Earth system. It is an error driven, rather than model-based control law, essentially accounting for the independence of model accuracy and linearization. An extended state observer is designed to estimate the states in real time by setting an extended state, that is, the sum of unmodeled dynamic and external disturbance. This total disturbance is compensated by a nonlinear state error feedback controller based on the extended state observer. A nonlinear tracking differentiator is designed to obtain the velocity of the spacecraft since only position signals are available. In addition, the system contradiction between rapid response and overshoot can be effectively solved via arranging the transient process in tracking differentiator. Simulation results illustrate that the proposedmethod is adequate for station-keeping of unstable Halo orbits in the presence of system uncertainties, initial injection errors, solar radiation pressure, and perturbations of the eccentric nature of the Earth’s orbit. It is also shown that the closed-loop control system performance is improved significantly using our method comparing with the general LQR method.


Introduction
In recent years there has been an increasing interest in libration points missions.The libration points, which are normally called equilibrium points or Lagrangian points, correspond to regions in space where the centrifugal forces and the gravitational forces from the Sun and the Earth cancel each other.The existences of periodic orbits and quasi-periodic orbits in the vicinity of collinear libration points have been proved and analyzed rigorously in celestial mechanics [1,2].These orbits offer potentially valuable opportunities for investigations concerning solar and heliospheric effects on planetary environment and highly precise visible light telescopes.
However, these libration point orbits are inherently unstable but controllable.Thus, additional control force is needed for a spacecraft to remain close to the nominal orbit.The challenges of station-keeping control emerge from highly accuracy, low computation burden, and minimal fuel cost control requirements under the condition of spacecraft dynamic uncertainties, unmodeled perturbations, and initial orbit injection errors [3][4][5][6].Hence, station-keeping control for libration point orbit missions is of virtual importance but with a great deal of difficulties.
The study of station-keeping control on libration point orbits has become a popular research topic ever since the problem was firstly proposed.A vast majority of the stationkeeping control methods are designed based on LTI model via local linearization at the libration points due to the high nonlinearity of the dynamic equation of libration point orbit.The nonlinear dynamic equation is obtained utilizing a Clohessy-Wiltshire-(CW-) like reference frame, which is 2 Mathematical Problems in Engineering widely used for Earth-centered spacecraft dynamic analysis [7,8].Breakwell et al. [9,10] firstly introduced classical optimal control strategies for Halo orbit missions.Erickson and Glass [11] specially analyzed the ISEE-3 mission to make this approach into implementation.Cielaszyk and Wie [12] employed a new LQR control method combined with a disturbance accommodating controller for Lissajous and Halo orbits maintenance based on LTI model.Di Giamberardino and Monaco [13] designed a nonlinear controller based on LTI model to solve the problem of tracking Halo orbit about the  2 point.However, the LTI model includes only 1st, 2nd, or 3rd order term of the gravitational force and is only effective in the neighborhood of the libration point.As a matter of fact, the higher the amplitude of the libration point orbit, the greater the influence of nonlinear factors.
To improve the modeling accuracy, LTV method is employed instead of the LTI approach.Gurfil and Kasdin [14] expended the LTI model to LTV model and designed a timevarying controller for formation flying in the Sun-Earth system.Barcolona Group proposed the so-called Floquet Mode Control [15][16][17], which employs nonlinear techniques to compute the invariant manifolds and determine the discrete impulsive maneuvers to cancel deviations along the unstable component based on LTV model.Howell and Pernicka [18] developed a target-point approach to maintain the spacecraft within some torus about the nominal Halo orbit.Furthermore, several other control strategies are developed to overcome the disadvantage of searching for an optimal controller for both Floquet Mode approach and target-point approach.Kulkarni et al. [3] extended the traditional  ∞ framework to periodic discrete LTV systems for stabilization of spacecraft flight in Halo orbits.Wang et al. [6] presented a nonlinear controller based on polynomial eigenstructure assignment of LTV model for the control of Sun-Earth  2 point stationkeeping and formation flying, without considering the system uncertainties.
Besides the LTI and LTV model-based method, existing works also involve methodologies directly developed from original nonlinear dynamical equation of motion.Rahmani et al. [19] solved the problem of Halo orbit control using the optimal control theory and the variation of the extreme technique.Xin et al. [4] used a suboptimal control technique (the  −  technique) to complete the mission of multiple spacecraft formation flying in deep space about the  2 point.Marchand and Howell [20] employed the feedback linearization approach for formation flight in the vicinity of libration points.Gurfil et al. [21] presented a novel nonlinear adaptive neural control methodology for deep-space formation flying.Bai and Junkins [22] proposed a modified Chebyshev-Picard interation method for station-keeping of  2 Halo orbits of Earth-Moon system.These nonlinear control methods have improved the control performance.However, the robustness under nonlinear system uncertainties and the request for low computation burden are unreasonably neglected, which should be considered before implementation.
Consequently, this paper derives an active disturbance rejection station-keeping control method considering mainly the following issues.
(1) Proposing an error driven, rather than model-based control law which takes into account system uncertainties, unmodeled disturbance, and orbit injection errors to achieve better robustness.
(2) Considering output feedback from practical viewpoint rather than full information, since only position signals of spacecraft can be measured.
(3) Better station-keeping performance as well as simpler computation burden.
A new nonlinear station-keeping control law based on active disturbance rejection control (ADRC) method, which refers to the so-called active disturbance rejection stationkeeping control (ADRSC) method in this paper, is proposed and analyzed.ADRC was firstly proposed by Han [23,24], and Sun [25] and has been successfully applied in various industrial processes such as vehicle flight control, ship tracking control, robot control, and power plant control.
The remainder of this paper is organized as follows.Section 2 presents the equations of motion and periodic reference orbits about the Sun-Earth libration point.Section 3 presents an introduction to the theory of ADRC synthesis.Section 4 presents an ADRSC method for nonlinear stationkeeping of periodic orbits around collinear libration point.Section 5 carries out simulation results to validate the effectiveness of the new method.Finally, Section 6 concludes the paper.

Equation of Motion
In this section, the dynamic models of spacecraft based on the circular restricted three-body problem (CR3BP) are established.Additionally, the relative Halo reference orbits derived from LP map method are presented.

Dynamics of CR3BP.
The CR3BP is one of the most common nonlinear models which investigate the relative motion around the libration points.As is shown in Figure 1, the infinitesimal mass  of spacecraft compared to the mass  1 of the Sun and the mass  2 of the Earth is assumed.It is further assumed that the two primary bodies rotate about their barycenter in a circular orbit under a constant angular velocity .A rotating coordinate system (, , , ) is defined with the origin set at the barycenter of the Sun-Earth system.The -axis is directed from the Sun towards the Earth.The -axis is perpendicular to the plane of rotation and is positive when pointing upward.The -axis completes the set to yield a right-hand reference frame.Normalization is performed by defining the distance  between the Sun and the Earth as the unit AU of length, the time of 1/ as the unit of time TU, and the total mass of the Sun and the Earth as the unit of mass.Normalization is obtained by setting  =  2 /( 1 +  2 ),  2 = , and located at (−, 0, 0) and  2 is located at (1 − , 0, 0).The wellknown equations of motion for CR3BP can be written in the dimensionless form [1] where the dot represents time derivative in the rotating frame, and the distance  1 = √ ( + ) 2 +  2 +  2 and  2 = √ ( − (1 − )) 2 +  2 +  2 .

Periodic Reference Orbits.
We can find libration points of the Sun-Earth system, denoted by  1 - 5 , by setting all derivatives in (1) to zero.For the Sun-Earth system,  1 is located at  = 0.989986055.It has been proved that the three collinear points,  1 ,  2 , and  3 , are unstable but controllable [1].And a family of periodic orbits and quasi-periodic orbits exist around the collinear points [2].Periodic solutions of the nonlinear equations of motion can be constructed using the method of successive approximations in conjunction with a technique similar to the Lindstedt-Poincare method [26].These periodic solutions are firstly presented by Richardson [27,28] and have been widely quoted and used.To obtain the periodic solutions around  1 , a new rotating coordinate system ( 1 , , , ) is defined with the origin set at  1 and the distance  between the Earth and  1 as the new unit of length, which is also shown in Figure 1.The relationship between the (, , , ) reference frame and the ( 1 , , , ) reference frame is as follows: where  = 0.010010905 [27].
The equations of the Halo orbit to the third order are given by where   (),   () and   () are the coordinates along , , and  axes.The Halo orbit designed for ISEE-3 mission [27] is selected as one of the target station-keeping orbits.The other reference orbit is selected with higher vertical amplitude.The values of the various constants in (2) of both Halo orbits are given in Table 1.

Introduction of ADRC
Due to complexity of modern systems, more attention has been paid on data-driven control scheme recently [29,30].
The ADRC is a typical data-driven control.It inherits from PID, using the error driven, rather than model-based control law to eliminate errors.Han's ADRC consists of three parts, a nonlinear tracking differentiator (TD) [31], which is used to arrange the ideal transient process of the system, an extended state observer (ESO) [32], which could estimate the disturbance from the system output and then the ADRC compensates the disturbances according to estimated values, and a nonlinear state error feedback (NLSEF) [24], which is used to get the control input of the system.Consider the following system: where  is the output variable,  is the control input,  is magnification factor, and () is the external disturbance; ( 1 ,  2 , (), ) includes three parts: modeling dynamics, uncertain dynamics, and disturbance.The ADRC approach makes an effort to compensate for the unknown dynamics and the external disturbances in real time without an explicit mathematical expression.

Mathematical Problems in Engineering
where V() denotes the control objective,  is speed factor and decides tracking speed, ℎ 0 is filtering factor, and fhan where  1 ,  2 , and  3 are observer gains,  is the error, and fal(, , ) is as follows: 3.3.NLSEF.NLSEF generates control voltage  for system by using the errors between the output of ESO and TD A nonlinear combination of errors signal can be constructed as [24]  0 = −fhan ( 1 ,  2 , , ℎ 1 ) , where the nonlinear coefficient  is selected as 0 <  < 1, and  is the proportional coefficients.
The controller is designed as

ADRSC for Station-Keeping
As is mentioned in Section 2, the reference Halo orbit around collinear libration orbit is inherently unstable, and station-keeping control approach is designed under model uncertainties and disturbances.In this paper, we extend the ADRC for the ADRSC.The structure of ADRSC algorithm is presented in Figure 2. ADRSC will be proposed and analyzed in this section from the following three aspects: discrete TD, discrete ESO, and discrete NLSEF.where  = 1, 2, 3 denote the , , and  axes component, respectively, V  () denotes the reference orbit.TD I provides transition process for nominal orbit V  and differential trajectory of set position value, that is, V 1 and its differential V 2 .

TD II.
It is designed for flight position tracking and estimating of the velocity of the spacecraft where  = 1, 2, 3 denote the , , and  axes component, respectively.For instance, state vector [V 11 V 12 ]  denotes the estimation of [ Ẋ]  which is the position and velocity vector of spacecraft trajectory along  axes.V  () denotes the flight trajectory of spacecraft.And  0 is the tracking speed parameter and ℎ 0 is the filter factor, which makes an effort of filtering.While the integration step is fixed, increasing the filtering factor will make the effort of filter better.The outputs of TD II, which are the spacecraft flight position and velocity, provide the inputs of ESO for feedback as described in Figure 2.

Discrete ESO for Estimating the Nonlinear Dynamics
Model Uncertainties and Unmodeled Disturbance.The equations of station-keeping of Halo orbits using ADRSC can be generally defined as  represents the uncertain dynamics vector, respectively.Note that (1) will be simulated as the real flight dynamic mode.Hence, robust performance will be presented under such formation of dynamic model uncertainties.
It is necessary to point out that a distinct improvement can be obtained using ADRSC for station-keeping since the dynamic model does not need to be expressively known which is different with the aforementioned studies.Define F(X, Ẋ) = g(X, Ẋ) + w.In fact, in the context of feedback control, F(X, Ẋ) is something to be overcome by the control signal, and it is therefore denoted as the "total disturbance." At this point, F(X, Ẋ) is extended as an additional state variable; that is, where  1 ,  2 , and  3 denote the corresponding axis components of F(X, Ẋ), respectively, and let Ḟ (X, Ẋ) = a(X, Ẋ), where a(X, Ẋ) is unknown.One can rewrite (14) as follows: where  1 ,  2 , and  3 represent the position, velocity, and   (X, Ẋ), respectively,  = 1, 2, 3.Then, one can use the following discrete nonlinear observer { 1 ,  2 ,  3 } to estimate state vector { 1 ,  2 ,  3 }: where {  ,   , ℎ  } are the NLSEF parameters described in (10). 1 and  2 are the state errors described in (9). 3 denotes unmodeled dynamics and disturbances   (X, Ẋ) which is estimated from (17).
The parameter ℎ  of NLSEF, which is also named filter factor of fhan, makes great contribution in NLSEF control performance.With the integration step fixed, ℎ  must be selected properly.It is noted that larger errors corresponding to larger ℎ  and smaller errors corresponding to smaller ℎ  will make a better performance.
In this paper, a parameter self-turning approach is firstly proposed for ℎ  selection as follows: where the parameters , ℎ 0 , and  of self-turning ℎ  can be easily chosen in practice, as will be presented in the simulation.

Disturbing Forces.
In any actual mission, the perturbation factors and the injection errors, coupled with the inherent unstable nature of the Halo orbits around the collinear libration points, will cause a spacecraft to drift from the periodic reference orbit [33].A robust performance in the presence of disturbance can be obtained by ADRSC, and the model of disturbance is not necessary for ADRSC as long as the disturbance is bounded.Physically speaking, deepspace disturbances are always bounded [21].To illustrate the capability of ADRSC to reject unknown disturbances, here we introduce the perturbative forces due to the eccentric nature of the Earth's orbit and the solar radiation pressure (SRP) disturbance.

The Perturbative Forces due to the
where   is the position vector of the Earth at the pericenter of the Earth's elliptical orbit from the libration point.  is the position vector similar to a circular orbit.

SRP Disturbance.
In the deep-space mission, SRP is another disturbance to account for.Here, we adopt a widely used model [34].According to this model, the disturbance acceleration to SRP is where  is a parameter that depends on the coefficient of reflectivity, the area, and mass of the spacecraft, the solar flux, and the speed of light.
_ r 1 = r 1 /‖r 1 ‖, and n is the attitude vector of the spacecraft in the rotating frame.To simplify (21), we assume that n = _ r 1 , which yields

Simulation Scenario and Numerical Values.
To compare the performance of ADRSC with linearization method, the simulation of an LQR controller [1] based on LTI system from local linearization around  1 point of nonlinear system equation ( 1) is also carried out.In order to compare the performance of ADRSC between low and large amplitude orbits, we select two orbits for simulation, one orbit for the ISEE-3 mission with vertical displacement 110,000 km and the other of 800,000 km.The control forces for station-keeping are provided by an ionic engine with the maximum thrust 60 mN which is currently commercially available [35].The mass of the spacecraft is designed as  = 500 kg.
The disturbing forces described in 4.4 are considered in the simulation.The eccentricity of Earth orbit  = 0.016675.And the SRP parameter is calculated assuming that the solar flux is   = 1358 W/m 2 , the speed of light is  = 3 × 10 8 m/s, the cross-sectional area is  = 3.5 m 2 , and the coefficient of reflectivity is  = 0.6, which results in  =   (1 + )/() = 5 × 10 −8 m/s 2 [21].It is also assumed that the initial position injection errors along the , , and  directions are 1000 km and the velocity injection errors are 1 m/s.The MATLAB Simulink software is used for simulation.The time step ℎ is set as 0.0001 TU for better accuracy of the spacecraft.The simulation time is chosen as 1800 days, almost 10 periods of Halo motion about  1 libration point.
In the simulation, considering the system is not full state feedback, the TD is designed to track the spacecraft flight position and obtain the velocity as well.Figure 3 illustrates the TD II tracking result of the spacecraft flight.In Figure 3, FP and FV denote the spacecraft flight position and velocity; TP and TV denote the TD tracking position and velocity, respectively.Figure 3 presents the accurately tracking of position signal and velocity signal along the , , and  axes, which will be used for feedback control.The position errors along , , and  directions can be kept within 1 km, which is very small compared with the orbit altitude of the ISEE-3 mission (110,000 km).Also, as shown in Figure 5(c) the control forces begin at 60 mN to drive the relative position errors to zero, but quickly reduce to the steady progress in 5 days.Hence, these results illustrate the effectiveness of the ADRSC of unstable Halo orbits near collinear libration points.

ISEE-3 Halo Orbit (𝑍
In order to describe the results of ADRSC more accurately when the spacecraft moves steady on the Halo orbit, here  we introduce the popular evaluation parameters for stationkeeping [3,33], the velocity increments Δ  , Δ  , and Δ  (in unit of m/s/T), and the mean absolute value of the position errors   ,   ,   (in unit of km) as follows [33]: where  is the orbitalperiod.
velocity increment and the mean absolute value of the position errors of ISEE-3 station-keeping with orbit injection errors and disturbance are given in Table 3.It can be seen that the mean absolute position error is 0.389 km, which illustrates the higher precision characteristic for the spacecraft station-keeping mission with ADRSC, comparing with that of LQR controller 2.965 km shown in Table 3.Meanwhile, the velocity increment of ADRSC under disturbing force and injection error is 36.567m/s/T which is approximately equal to that of LQR controller 36.705m/s/T.Note that there is a trade-off between accuracy of the spacecraft to follow the Halo orbit and the thrust usage [3].The earlier approaches, which use impulsive maneuvers, are designed for station-keeping missions and do not have stringent orbittracking requirements.While high accuracy orbit stationkeeping is definitely required for current missions, which use continuous low-thrust propulsion.

Summary.
From the above analyses and control results between ADRSC and LQR, we can find the following.(1) The ADRSC, which is an error driven control method, is adequate for the station-keeping control of unstable orbits without any knowledge about the spacecraft dynamic model.
(2) The ADRSC presents praiseworthy station-keeping performance with orbit injection errors as well as unmodeled disturbances such as the SRP and the perturbative forces due to the eccentric nature of the Earth's orbit using an ionic engine with maximum thrust 60 mN.
(3) The ADRSC approach has better station-keeping control ability and higher orbit maintenance accuracy compared with LQR controller despite of the level of both the control methods thrust usage.
(4) The ADRSC has better robust performance compared with the general station-keeping method depending on spacecraft dynamic model.
There has been a great deal of commendable research on unstable libration point orbits station-keeping as mentioned before.A major point of departure between ADRSC method and earlier approaches to station-keeping is that ADRSC is an error driven rather than model-based control law, which can inherently get across the dependency on model accuracy as well as the drawbacks of linearization.With the combination of TD, ESO, and NLSEF, the unmodeled disturbances as well as the unmodeled system dynamic can be compensated in real time.Thus, fast response-time requirement and high accuracy of orbit maintenance requirement can be satisfied by ADRSC.
It is important to keep in mind not only the tracking accuracy but also the robustness of the station-keeping controller, as the space environment cannot be accurately modeled as well as the internal and external disturbance.ADRSC extends the unmodeled spacecraft dynamic and the disturbance as a state, which can be estimated from

4. 1 .
Discrete TD for the Tracking of the Reference Orbit and the Spacecraft Trajectory.Assuming that only position signals of spacecraft can be detected and used for feedback, there are two different kinds of discrete TDs designed for stationkeeping control.4.1.1.TD I.It is designed for reference Halo orbit tracking

4 . 3 .
Discrete NLSEF for Station-Keeping.Since the states, unmodeled dynamics, and disturbances have been estimated by ESO presented in Section 4.2, state errors between the output of ESO and TD I are combined for feedback control law, large errors corresponding to low gains and small errors corresponding to high gains

Figure 5 :
Figure 5: The relative position errors, velocity errors, and control input of ADRSC control of ISEE-3.
Figures 7 and 8 illustrate the ADRSC results for Halo orbit ( = 800, 000 km) with orbit injection errors as well as disturbing force.From Figures 8(a) and 8(b), one can find that ADRSC controller can drive the relative position errors and velocity errors along , , and  steadily to zero within 7.5 days, compared with LQR controller 25 days in Figures 9(a) and 9(b).During the steady periodic, the position errors are kept 2 km which illustrate high accurately maintained ability of ADRSC of high Halo orbit.

Figure 6 :
Figure 6: The relative position errors, velocity errors, and control input of LQR control of ISEE-3.

Table 1 :
[31]values of the various constants in(2)for ISEE-3 ( = 110,000 km) orbit and Halo orbit ( = 800,000 km) near  1 libration point of the Sun-Earth system.For ISEE-3 ( = 110,000 km) orbit and for Halo orbit ( = 800,000 km).The TD has the ability to track the given input reference signal with quick response and no overshoot by providing transition process for expected input V and differential trajectory of set value, that is, V 1 , and its differential V 2 .One feasible second-order TD can be designed as[31]

Table 2 :
Initial orbits parameters of Halo missions without injection errors.

Table 3 :
The velocity increments and the mean absolute value of the position errors of ADRSC and LQR control.