Abstract
Due to the influence of various perturbations of space, satellites flying in formation cannot maintain specific configurations for long durations [1, 2]. In order to ensure that formation configurations are able to meet the requirements of space missions, it is important to maintain control of formation configurations. This is an urgent problem to be solved. The traditional control method for controlling formations is based on the average orbit element, and uses the assumption that the average orbit element deviation and the instantaneous orbit element deviation are approximately equal. However, the continuous control system is more difficult to achieve in engineering practice. Using a LQR (linear quadratic regulator) optimal control algorithm and SDRE (state-dependent Riccati equation) optimal control algorithm to maintain the formation flying [3, 4]. Through simulation, it was found that when using the SDRE controller in the system transition process time is shorter than when the LQR controller is used, and fuel consumption is less for the SDRE controller than for the LQR controller.
1 Introduction
With the development of small satellite technology, research on a number of small satellite co-operation to achieve a common function has become one of the hot issues in the field of space research. Compared with a traditional single large satellite [5], small satellites have the advantages of low quality, low cost and high reliability. When a satellite is broken, we can replace the bad satellite base on orbit reconstruction, thus extending the life of the entire system. As a result, small satellite networks have received wide attention. Since the beginning of the 1990s, the United States has launched a “ION-F”, “TechSat-21” and other research programs [5].
A small satellite is smooth between collaborative working, because between their formation keep close, so its development promotes the research of the relative motion theory of the spacecraft [6]. In formation flying, the traditional control method is based on the average orbit element, and uses the assumption that the average orbit element deviation and the instantaneous orbit element deviation are approximately equal.However, taking into account the actual project in the continuous control is more difficult.Based on the dynamics of formation flying, two controllers are designed: an LQR controller based on a linear model, and an SDRE controller based on a nonlinear model.
2 Formation flying dynamics
The dynamic model of a satellite in orbit can be obtainedaccording to Kepler’s equation:
From Figure 1, the equations of motion of the lead satellite can be obtained as [7]:
The equations of motion of the following satellite are, then:
In equation(3), f is the sum of perturbation and control forces exerted on the satellite [8], Hypothesis n represents the leader satellite’s orbital angular velocity,
And from (4) and (5) we can obtain:
When the distance between the leader satellite and the following satellite is far less than the distance between the center of the Earth and the center of the leader satellite, the nonlinear term in the equation can be simplified by
Making a series of linearization for nonlinear equation(6), then:
In equation (8),
When the leader satellite runs in a circular orbit, then equation(8) can be simplified as follows:
3 Design of LQR controller design
LQR (linear quadratic regulator) is one of the earliest and most mature state space design methods in modern control theory [9]. The optimal control law of state linear feedback is obtained by LQR, and it is easy to achieve closed-loop optimal control [10].
The linear dynamic model of satellite formation flying is rewritten as a state-space expression.
In equation(10), u represents the system control variables:
x indicates the system state variable:
A(x) is the state matrix of the system:
In the matrix,
B(x) is a control matrix:
In the matrix,
Assume the altitude of the leader-satellite is 800 km, the orbital radius of the leader-satellite is rc = 7.2 × 105 m, the orbital angular velocity is n = 0.045 rad/sec = 2.5 × 10−4 rad/sec. At this time the state space expression is:
In equation(11),
This indicates that the system can be controlled.
The objective function of LQR theory is:
In equation(13), xd is an ideal state, Q is used as the weight matrix of the error in the optimization process, and is a positive definite constant matrix of 6 × 6; R is the weight matrix of the control variables in the optimization process and is a positive semi-definite constant matrix of 3 × 3.
Solving Riccati equation:
P(x) can be solved.
Feedback matrix:
The control law of the linear system in the performance index is:
4 Design of SDRE controller design
The SDRE method is a kind of nonlinear control method. When applying this method, the nonlinear system dynamic equations must first be converted into SDC (State-dependent coefficients) forms [11].
For the nonlinear affine system:
SDC forms can be obtained by pseudo linearization of nonlinear affine systems.
In equation(18):
and
If the leader-satellite is running on a circular orbit, then
At this time the state space expression is:
In equation(19),
rc is the orbital radius of the leader-satellite, rc = 7.2 × 105 m; the orbital angular velocity is n = 0.045 deg/sec = 2.5 × 104 rad/sec; μ is the geocentric gravitational constant, μ = 3.986 × 1014 m3/s2.
In equation(20),
System controllability matrix:
This calculation shows that the system is controlled point by point.
The objective function of SDRE theory is:
In equation(23), xd is an ideal state, Q(x) is used as the weight matrix of the error in the optimization process, and R(x) is used as the weight matrix of the control variables in the optimization process. Q(x) and R(x) are functions of the system state x.
Solving Riccati equation:
P(x) > 0 can be solved.
The control law for the nonlinear system in the performance index is:
5 Numerical example and simulation results
Under ideal conditions, the initial relative position and velocity of the leader-satellite and the following-satellite formations are:
During a disturbance, the initial relative position and velocity error become:
So the initial relative position error and velocity error between the leader-satellite and the following-satellite are: x − xd = [ 20 30 −10 0.3 0.2 −0.1]
In order to facilitate the comparison of the LQR and SDRE methods, we do not consider the influencing factors of the percussion force. the weight matrix Q and R are respectively:
In the SDRE controller, the weight matrices Q and R are respectively:
Using MATLAB to get a variety of state simulations results are as follows:
6 Conclusions
Figures 2 to 4 compare the LQR controller algorithm and SDRE controller algorithm of the position error simulation. In 0-20 seconds, for the X-axis position error, the SDRE method tends to zero faster than the LQR. The curves for the Y- and Z-axes take at least 35 seconds to reach zero. Figures 5 to 7 compare the LQR and SDRE controller algorithms for the velocity error simulations. It can be seen that the SDRE method has less overshoot than the LQR method. From 0-7 seconds, for the X- and Y-axes, the SDRE method is faster than the LQR Method at tending to zero. From 7-35 seconds, the LQR method is faster than the SDRE method at tending to zero. But the Z-axis is the opposite. Figures 8 to 10 compare the LQR and SDRE controller algorithms of the control variables simulation . It can be seen that the range of control variables for the SDRE method is less than for the LQR method. At the same time, to reach steady state, the fuel consumed by the SDRE method is less than that of the LQR method.
In summary, the simulation results show that, in the system transition, both processing time and fuel consumption are lower for the SDRE controller than for the LQR controller.
References
[1] Lee K. H., Lunar orbiter trade study and conceptual design of onboard propulsion system, Journal of spacecraft and rockets, 2011, 48, 346–354.10.2514/1.50459Search in Google Scholar
[2] Perez-Chavela E., Tamayo C., Tamayo relative equilibria in the 4-vortex problem bifurcating from an equilateral triangle configuration, Applied Mathematics and Nonlinear Sciences, 2016, 1, 301–310.10.21042/AMNS.2016.1.00025Search in Google Scholar
[3] Lee J., Kim J.-S., Shim H., Disc margins of the discrete-time LQR and its application to consensus problem, International Journal of Systems Science: The Theory and Practice of Mathematical Modelling, Simulation, Optimization and Control in Relation to Biological, Economic, Industrial and Transportation Systems, 2012, 43, 1891–1900.10.1080/00207721.2011.555012Search in Google Scholar
[4] Liang Y.-W., Lin L.-G., Analysis of SDC matrices for success fully implementing the SDRE scheme, Automatica, 2013, 49, 3120–3124.10.1016/j.automatica.2013.07.026Search in Google Scholar
[5] Vaddi S., Sengupta P., Formation establishment and reconfiguration using impulsive control, Journal of Guidance, Control, and Dynamics, 2005, 28, 262–268.10.2514/1.6687Search in Google Scholar
[6] Guibout V., Scheeres D., Solving relative two point boundary value problems: Applications to spacecraft formation flight transfers, Journal of Guidance, Control, and Dynamics, 2004, 27, 693–704.10.2514/1.11164Search in Google Scholar
[7] Baker D.M.H., Head J.W., Neumann G.A., Smith D.E., Zuber M.T., The transition from complex craters to multi-ring basins on the moon: Quantitative geometric properties from lunar reconnaissance orbiter lunar orbiter laser altimeter (LOLA) data, Journal of geophysical research, Planets, 2012, 117, E00H16–1–E00H16–29.10.1029/2011JE004021Search in Google Scholar
[8] Abouelmagd E.I., Guirao J.L.G., On the perturbed restricted three-body problem, Applied Mathematics and Nonlinear Sciences, 2016, 1, 123–144.10.21042/AMNS.2016.1.00010Search in Google Scholar
[9] Shi L., Yuan Y., Chen J., Finite horizon LQR control with limited controller-system communication, IEEE Transactions on Automatic Control, 2013, 58, 1835–1841.10.1109/TAC.2012.2235722Search in Google Scholar
[10] Liu X., Wu Y., Zhang Y., Xiao S., A control method to make LQR robust: A planes cluster approaching mode, International Journal of Control, Automation, and Systems, 2014, 12, 302–308.10.1007/s12555-012-0435-0Search in Google Scholar
[11] Choi H.H., SDRE-based near optimal nonlinear controller design for unified chaotic systems, Nonlinear dynamics, 2012, 70, 2063–2070.10.1007/s11071-012-0598-5Search in Google Scholar
© 2017 Zhang Ke et al.
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License.