Abstract
This paper proposes a method to design the robust parametric control for autonomous rendezvous of spacecrafts with the inertial information with uncertainty. We consider model uncertainty of traditional C-W equation to formulate the dynamic model of the relative motion. Based on eigenstructure assignment and model reference theory, a concise control law for spacecraft rendezvous is proposed which could be fixed through solving an optimization problem. The cost function considers the stabilization of the system and other performances. Simulation results illustrate the robustness and effectiveness of the proposed control.
1. Introduction
With further exploration into the space, a set of complex missions is in the space development agenda such as large-scale structure assembling, sending and picking up astronauts, and repairing, saving, and docking, orbital propellant resupply based on the autonomous rendezvous technology [1]. Due to the essential position, many scholars have been focusing on the control problem during rendezvous and some results enlightened deeper research. In the approximately circular orbit, C-W equations [2], derived by Clohessy and Wiltshire, have been widely applied for the depiction of the relative motion between neighboring spacecrafts. The early stage of control design based on C-W equation revealed a number of open-loop methods such as V-bar, R-bar, dual impulsive, and multiple impulsive [3]. With the benefits of control theory flourishing, plenty of advanced control methods are used to solve the rendezvous problems such as using artificial potential function in [4], sliding mode control in [5], adaptive control in [6], and H-infinity theory in [7].
Though the C-W equation supplies an explicit description of the relative motion for spacecrafts, there is an obstacle when applied in reality that the real-time angle velocity of the target spacecraft could not be obtained accurately as result of detection errors and perturbation from environment. This parameter uncertainty affects the control force and system stability directly. It is necessary to investigate the uncertain model for spacecraft rendezvous not depending on accurate value of real-time angle velocity. The traditional robust control method could deal with parametric uncertainty to recognize rendezvous but some expecting system characters are hard to be included during the control design.
In this paper, the spacecraft rendezvous problem with uncertain parameter would be solved by robust parametric method which allows freedom to improve system performance. The robust control integrates eigenstructure assignment and model reference theory to propose a concise control law for spacecraft rendezvous which takes into consideration the system performance such as the control constraints and fuel saving. In the rest of this paper, a relative motion model with uncertainty for the spacecraft rendezvous is to be established; the design of the robust parametric control law follows; besides, we apply the robust parametric control for an example to illustrate the effectiveness of this design approach.
2. Problem Formation
2.1. Equations of Motion
The coordinate frame for the two spacecrafts rendezvous is based on the target spacecraft orbit, described in Figure 1. We set the original point at the target’s mass center; , , and indicate along-track, the radial, and out of plane components of the position vector of the chaser satellite in the target satellite’s local-vertical-local-horizontal (LVLH) frame, respectively.
The spacecraft rendezvous in the circle orbit would obey the C-W equations where , , and stand for the relative position between the chase spacecraft and the target spacecraft; represents the average angle velocity of the target spacecraft; , , and stand for the control acceleration on each axis.
According to the equation, the state and control vector can be described as and output vector can be Then, we get where and represents the matrix with the values of all elements equal to zero; represents the unit matrix.
2.2. Problem Description
The classical C-W equations need accurate angle velocity simultaneously which is difficult to obtain due to the detection error. Therefore, we consider the uncertain item to the angle velocity to make the system model closer to reality.
When the angle velocity changes are the system model can be described as where
The object of the designing control law is to recognize where is the output of the system and represent the reference relative position between chase spacecraft and target spacecraft. Meanwhile, the uncertainty brings trouble to the stability of the system which would be taken into consideration during designing the control law.
3. Design of Robust Parametric Control
The design of the control law aims at reaching the reference point of the chase spacecraft and keeping the closed loop system stable. It could be separated into two parts as stabilization controller and trajectory tracking controller.
3.1. Trajectory Tracking Controller
To begin with, we would design the tracking controller based on the model reference theory. Lemma 1 supplies theoretical evidence for the linear tracking problems referred to [8].
Lemma 1. For the system, if the stabilization feedback control law exists, the control law following the form as would obtain the result of tracking reference signal, which means that where feedforward control law could be calculated from the following equation: and , could be calculated as
According to Lemma 1, the rendezvous system could track the reference position when the feedback control law stabilizes the system. Then, the critical task of designed controller is to find a robust stabilization control law . Regarding the eigenstructure assignment of linear system, some useful results would be utilized in the later part which are from [8].
Lemma 2. Suppose , , and is controllable. , , are a set of complex numbers, which are symmetric about the real axis. Then, the matrices and satisfying are given by where , , are arbitrary vectors which satisfy and and are right comprime polynomial matrices satisfying
For the rendezvous system in this paper, we could calculate according to Lemma 2 as
Lemma 2 supplies a concise parametric formula for state feedback law in which the poles of the closed-loop system are included. Proper poles would not only guarantee the system stabilization but also enhance system characters through optimization in some specific fields. Besides, the parametric method offers all kinds of freedom to design the control system with the free parametric vectors , , which enable us to adjust these parameters for system stabilization.
3.2. Stabilization Controller
Using the control law the closed-loop system can be described as where When is a nondefective matrix and the closed-loop system owns the required poles , the sufficient condition for the system stabilization with the uncertainty item is [9] where is a symmetric positive definite solution of the following: Lemma 3 provides the parametric expression for based on the eigenstucture of the system.
Lemma 3. The solution to (23) has the following parametric representation: where and , , , are respectively the eigenvalues and corresponding eigenvectors of .
For a better stabilization system, we hope to minimize which is the result of some adjustments for the parameters and .
3.3. Optimization of Control Law
We have established the connection between the system characters and the parameters and through the design of the control law. Therefore, the design problem for the rendezvous system can be converted into the following nonlinear optimization problem:
where , , , and specify the desired areas of the closed-loop eigenvalues.
The performance index is chosen as follows: where denotes the initial state of system; , , and are the weighting factors. The first part of (27) is chosen due to the consideration of the input constraint. The second item of (27) takes into consideration fuel consumption. The last part of (27) is used for global stability of the rendezvous system.
The optimization discussed above could be solved resorting to the optimization tool in MATLAB for its convenience. Then, the poles of the system and free parametric vectors would be fixed to calculate the feedback matrix for the robust control.
4. Numerical Simulations
In this section, our control law designed through the method proposed above would be tested by an example of spacecrafts in the final approaching in rendezvous mission. With the assumption that the target is in the geosynchronous orbit, we set the standard angle velocity rad/s. Suppose the initial state vector is and the desired final state is where is final time of the rendezvous mission. The weighting factors in the performance index are Specify the desired closed-loop eigenvalue regions as
By solving the optimization problem (26a), (26b), and (26c), we gain the poles of the system and free parametric vectors , , According to Lemma 2 and (15), we get
Then, the stabilization control law to (15) is
Assume the uncertainty which leads to Meanwhile, the closed-loop poles and parametric vectors could be used to calculate as It is obvious that the rendezvous process could reach the desired state with the control law when the inequality (22) has satisfied. The rendezvous trajectory and the relative position of the two spacecraft are showed in Figures 2 and 4 and the effectiveness could be proved simultaneously. Due to the proper optimization function, the control inputs have been constrained to which can be seen in Figure 3. The motion in every axis direction changes smoothly so that the simulation system gets closer to the real engineering requirement showed in Figure 3.
5. Conclusion
This paper has proposed a method to design the robust control law for spacecraft rendezvous in the final approach subject to parameter uncertainty in near circle orbit. Based on the eigenstructure assignment and model reference theory, the control law is constructed with the closed-loop poles and design freedom. Through solving an optimization problem, we obtain the poles and parametric vectors to calculate the control law which has been proved useful by simulation.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by National Natural Science Foundation of China (no. 61074111), Natural Science Foundation of Jilin Province, China (no. 20140520061JH), and Educational Commission Foundation of Jilin Province, China (no. 2013122).