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 , 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 .
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 . 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: (1)
represents a radial vector from the Earth’s center of mass to the satellite, μ is the gravitational constant, μ = GM, with M as the Earth’s mass.
The equations of motion of the following satellite are, then: (3)
In equation(3), f is the sum of perturbation and control forces exerted on the satellite , Hypothesis n represents the leader satellite’s orbital angular velocity, as the relative distance between the center of the leader satellite and the center of the following satellite. Then: (4) (5)
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 (7)
Making a series of linearization for nonlinear equation(6), then: (8)
In equation (8), α is the length of the semi major axis, θ represents the true anomaly and e indicates the eccentricity.
When the leader satellite runs in a circular orbit, then equation(8) can be simplified as follows: (9)
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 . The optimal control law of state linear feedback is obtained by LQR, and it is easy to achieve closed-loop optimal control .
The linear dynamic model of satellite formation flying is rewritten as a state-space expression. (10)
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: (11)
In equation(11), (12)
This indicates that the system can be controlled.
The objective function of LQR theory is: (13)
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: (14)
P(x) can be solved.
Feedback matrix: (15)
The control law of the linear system in the performance index is: (16)
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 .
For the nonlinear affine system: (17)
SDC forms can be obtained by pseudo linearization of nonlinear affine systems. (18)
If the leader-satellite is running on a circular orbit, then = 0.
At this time the state space expression is: (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. (20) (21)
System controllability matrix: (22)
This calculation shows that the system is controlled point by point.
The objective function of SDRE theory is: (23)
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: (24)
P(x) > 0 can be solved.
The control law for the nonlinear system in the performance index is: (25)
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:
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.
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. CrossrefGoogle Scholar
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. Google Scholar
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. CrossrefGoogle Scholar
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. Web of ScienceGoogle Scholar
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. CrossrefWeb of ScienceGoogle Scholar
About the article
Published Online: 2017-06-14