Next Article in Journal
Multifidelity Sensitivity Study of Subsonic Wing Flutter for Hybrid Approaches in Aircraft Multidisciplinary Design and Optimisation
Next Article in Special Issue
Decentralized Control of Nanosatellite Tetrahedral Formation Flying Using Aerodynamic Forces
Previous Article in Journal
Utilization of FPGA for Onboard Inference of Landmark Localization in CNN-Based Spacecraft Pose Estimation
Previous Article in Special Issue
Leader–Follower Synchronization of Uncertain Euler–Lagrange Dynamics with Input Constraints
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Satellite Motion Control System Parameters on Performance of Space Debris Capturing

1
Moscow Institute of Physics and Technology, State University, Institutsky Lane 9, Dolgoprudny, 141701 Moscow Region, Russia
2
Keldysh Institute of Applied Mathematics RAS, Miusskaya sq. 4, 125047 Moscow, Russia
*
Author to whom correspondence should be addressed.
Aerospace 2020, 7(11), 160; https://doi.org/10.3390/aerospace7110160
Submission received: 18 September 2020 / Revised: 27 October 2020 / Accepted: 4 November 2020 / Published: 6 November 2020
(This article belongs to the Special Issue Small Satellite Formation Flying Motion Control and Attitude Dynamics)

Abstract

:
Relative motion control problem for capturing the tumbling space debris object is considered. Onboard thrusters and reaction wheels are used as actuators. The nonlinear coupled relative translational and rotational equations of motion are derived. The SDRE-based control algorithm is applied to the problem. It is taken into account that the thrust vector has misalignment with satellite center of mass, and reaction wheels saturation affects the ability of the satellite to perform the docking maneuver to space debris. The acceptable range of a set of control system parameters for successful rendezvous and docking is studied using numerical simulations taking into account thruster discreteness, actuators constrains, and attitude motion of the tumbling space debris.

1. Introduction

Further exploration of near-Earth space in the short term will become impossible without solving the problem of space debris removal. A number of international projects are aimed at its solution [1,2]. Space debris removal approaches can be divided into two classes: passive and active [3,4,5]. In the case of passive approach the space debris are removed with the help of external natural forces, for example, aerodynamic drag in low Earth orbits [6,7], solar pressure [8], ionospheric drag [9] etc. Active removal of inactive spacecrafts and rocket stages often involves the use of special small satellites that can attach themselves on the space debris object or capture it with a manipulator [10], net [11,12] or harpoon with tether [13,14], and change its orbit using the on-board motion control system. Active debris removal implies autonomous relative motion control in order to achieve relative state vector required for capturing. An onboard propulsion is often considered for the translational motion control and reaction wheels are for the attitude maneuvers. However, small satellites have restrictions on mass, size, energy, on-board computing power, and the composition of the control system equipment, which complicates the on-board control algorithms for the main modes of motion at the mission design stage, taking into account the limited capabilities of the satellites. The reaction wheels may experience saturation which must be avoided, the saturation may be caused by thruster misalignment or by high angular velocity of the tumbling debris. These situations may lead to the failure of the object capturing. That is why it is required to study the influence of the control system parameters on rendezvous and capturing maneuver capabilities.
The problem of relative translational and attitude motion control is well studied and a big variety of control approaches is developed. For example, sliding control-based algorithms [15,16] are developed for the relative orbit-attitude tracking problem, for the rendezvous problem the swarm particle optimization algorithm is applied for the required trajectories generation [17], for the docking stage with non-cooperative object the majority of the proposed control algorithms are fuel-optimal or time-optimal [18,19,20,21,22]. The fuel-optimal trajectory generation algorithms are computationally intensive, its implementation in a real-time system is very challenging. Therefore, the optimal algorithms for calculating trajectories are often replaced by computationally simpler and faster non-optimal ones [23,24], however their performance is strongly dependent on initial conditions at the docking stage. As a compromise between two approaches a feedback control law developed for minimization of some defined cost-function can be applied to the problem. A linear quadratic regulator (LQR) is a well-known example of such an algorithm, though the relative motion equations are highly non-linear. To overcome this inconsistency the motion equations are linearized in the vicinity of the current state vector and LQR-like State-Dependent Riccati Equation-based (SDRE) control algorithm is applied [25,26,27]. In [28,29] a comparative study between SDRE and LQR is presented, and SDRE showed its advantages considering fuel consumption, rendezvous time, and trajectory accuracy. The SDRE-based algorithms are used to address various problems such as the position and attitude control of a single spacecraft [30] or relative motion control in satellite formation flying [29]. For the problem of space debris object capturing the kinematic coupling effect must be taken into account when the relative motion of not centers of mass of two bodies but motion between two defined body-fixed points is considered as in [31]. The paper [32] studies the application of the SDRE-based control for this type of relative motion equations. The main contribution of the current paper is in the study of the influence of the control system parameters on the performance of the SDRE-based control algorithm on the relative motion during the capturing, taking into account the reaction wheels saturation and thrusters misalignment.
This paper considers a spacecraft with thrusters installed on board to control the center of mass motion and it is equipped with reaction wheels for the attitude control [33]. The motion of a non-cooperative object relative to the spacecraft is considered known for the sake of simplicity. In real cases the target motion is estimated by relative motion determination system with some errors and uncertainties. The purpose of the work is to develop an algorithm for controlling both the center of mass motion and the angular motion to achieve the required relative position and attitude of the spacecraft relative to a non-cooperative object, which is necessary for the capturing. It is assumed that the object of space debris has an axis of dynamical symmetry and rotates freely under the influence of a gravitational torque and its tumbling motion is similar to nutation. To capture space debris it is necessary to match a desired point on the spacecraft reference frame with a point on the object surface [31,34]. As a capturing system can be considered a robotic arm capable of catching on an element of the space debris body.
SDRE control algorithm to achieve relative motion along a desired trajectory is developed in the paper [35]. The SDRE method requires linearization of the equations of motion in the vicinity of the current state [27], and the optimal controller coefficients are calculated as a result of solving the Riccati equation at each control cycle [36]. Nonlinear coupled translational and rotational motion equations of the spacecraft relative to the object are used. It is assumed that the value of the thrust is limited and discrete, and the thrust vector has a misalignment, which produces an additional disturbing torque and an external coupling effect [37]. Because of control limitations the successful capture of a space debris object is possible only in the region of acceptable values of the system parameters. The parameters include the initial conditions for the motion equations, the magnitude of the misalignment of the thrust, the value of the control constraints, the position of the capture point on the object, the parameters of the angular motion of the object. This paper is the continuation of authors pervious work [35], it considers a more generalized equations of motion and proposes a methodology for assessing the acceptable range of these parameters using a numerical study of the system motion.

2. Equations of Motions and Control Law

In this section a short introduction to equations of motion and control law is presented. The relative state vector consist of relative position and velocity between the centers of mass of a chaser (active satellite) and target (passive space debris) and relative attitude quaternion and angular velocity. Since for the capturing it is necessary to align two points on the surfaces of the satellites and in the object (the position of the capturing system of the chaser and capturing point of the target, see Figure 1), the relative rotational and translational motion equation are coupled.

2.1. Relative Rotational and Translational Equations of Motion

For modeling the dynamics of this problem the coupled sets of nonlinear rotational and translational relative equations of motion of spacecraft with respect to the non-cooperative object are used. The Equations (1) and (2) describe the changes of angular momentum of chaser H C and the target H T in inertial reference frame:
d H C d t | I = d H C d t | C + ω C × H C = N C + T C .
d H T d t | I = d H T d t | T + ω T × H T = N T
where | I , | C , and | T are derivative in the inertial, chaser-fixed, target-fixed reference frame respectively, ω C ,   ω T are chaser and target angular velocities, N C   and   N T are external torques acting on the spacecraft and debris respectively, and T C is the control torque. The chaser satellite is equipped with reaction wheels, so its angular momentum is calculated as follows:
H C = I C ω C + h W C H t = I T ω T
where h W C is the angular momentum produced by reaction wheels, I C and I T are chaser and debris inertia tensor. It is assumed that the chaser satellite is equipped with three reaction wheels, the angular momentum of each wheel is aligned with principal axis of the satellite. All the variables in Equations (1) and (2) are expressed in the corresponding body-fixed reference frame.
Since the target is passive and the chaser is equipped with control system, the relative rotational equation of motion should be expressed in target-fixed reference frame. By subtracting (2) from (1) and taking into account the transition matrix D ( q ) between the chaser reference frame to target reference frame the following relative angular motion equation is obtained:
I T ω ˙ T = I T D ( q ) I C 1 [ D ( q ) 1 ( ω T + ω T T ) × I C D ( q ) 1 ( ω T + ω T T ) D ( q ) 1 ( ω T + ω T T ) h W C h ˙ W C + T C + N C ] I T ω T T × ω T + [ ω T T × I T ω T T ]
Here ω T = D ( q ) ω C C ω T T is the relative angular velocity which is expressed in the target’s reference frame, superscript defines the reference frame where the value is expressed, “ C ” stands for chaser and “ T ” stands for non-cooperative target. The direction cosine matrix of rotation can be expressed using quaternion q that defines the transition from chaser to target reference frame.
The external torque N T acting on the target includes only gravitational torque, but N C also includes the torque N t h caused by thruster vector misalignment from the chaser center of mass.
N t h = r t h × F t h
where F t h is a thrust vector of the thruster, r t h is the radius vector from center of mass to the center of thrust application.
The kinematical equation can be written for quaternions as follow:
q ˙ = 1 2 Q ( q ) ω T
In (5) Q is given as:
Q = [ q 4 q 3 q 2 q 3 q 4 q 1 q 2 q 1 q 4 q 1 q 2 q 3 ]
In Figure 1 the local-vertical-local-horizontal coordinate system of target is shown as C T x T y T z T . C T ξ T η T ζ T is the target body-fixed reference frame and C C ξ C η C ζ C is the chaser body-fixed reference frame of the chaser.
The general nonlinear relative equations of motion centers of mass of two objects are given as follows:
x ¨ 2 ω O T y ˙ ω ˙ O T y ω O T 2 x = μ ( r T + x ) [ ( r T + x ) 2 + y 2 + z 2 ] 3 2 + μ r T 2 + a x
y ¨ + 2 ω O T x ˙ + ω ˙ O T x ω O T 2 y = μ y [ ( r T + x ) 2 + y 2 + z 2 ] 3 2 + a y
z ¨ = μ z [ ( r T + x ) 2 + y 2 + z 2 ] 3 2 + a z
where ρ 0 = [ x ; y ; z ] T is radius-vector of chaser center of mass in LHLV reference frame with target in the origin C T , r T is the value of the radius-vector of the target center of mass in inertial reference frame; a = [ a x   a y   a z ] T is the control acceleration, which is produced by thrusters; ω O T is the orbital angular velocity of the target; μ is the Earth gravitational constant. In the case of near circular orbits and r T ρ 0 these equations can be linearized and the well-known Clohessy-Wiltshire can be obtained [38]. To provide control acceleration it is assumed that the chaser satellite is equipped with six thrusters, a pair of thrusters with opposite directions for each principal axis of satellite.
To develop a control algorithm for capturing the target it is required to obtain the equations of relative motion of two points—the point of the chaser satellite where the capturing mechanism is placed, and the point of the target object for which the object can be captured. The coupling effect between translational and attitude motion should be taken into account. Consider two arbitrary points of the target r T j and of the chaser r C i (see Figure 1). For these two points the relative motion equations can be derived from the following relation for the relative vector ρ i j between these points:
ρ i j = ρ 0 + r C i r T j
All of the vectors are expressed in LVLH reference frame. Calculating the first and then second derivative the following equations can be obtained:
ρ ˙ i j = ρ ˙ 0 + ω × r C i
ρ ¨ i j = ρ ¨ 0 + ω ˙ × r C i + ω × ( ω × r C i )
Taking into account that ρ i j = [ x i j   y i j   z i j ] T Equation (12) can be rewritten in the following form:
x ¨ i j 2 ω O T y ˙ i j ω ˙ O T y i j 3 ω O T 2 x i j = a x + K 1 y ¨ i j + 2 ω O T x ˙ i j + ω ˙ O T x i j = a y + K 2 z ¨ i j + ω O T 2 z i j = a z + K 3
where K = [ K 1   K 2   K 3 ] T are terms which consist of ω and express the coupling between (13) and (4):
K 1 = [ ω y ( ω x y C i ω y x C i ) ω z ( ω z x C i ω x z C i ) ] + ω ˙ y z C i ω ˙ z y C i + 2 ω O T ( ω z x C i + ω x z C i ) + ω ˙ O T ( y C i + y T j ) + 3 ω O T 2 ( x C i + x T j ) + μ r T 2 μ r T r C 3
K 2 = [ ω z ( ω y z C i ω z y C i ) ω x ( ω x y C i ω y x C i ) ] + ω ˙ z x C i ω ˙ x z C i + 2 ω O T ( ω y z C i + ω z y C i ) + ω ˙ O T ( x C i + x T j ) + ω O T 2 ( y C i + y T j )
K 3 = [ ω x ( ω z x C i ω x z C i ) + ω y ( ω y z C i ω z y C i ) ] + ω ˙ x y C i ω ˙ y x C i + ω O T ( z C i z T j )
where r T j = [ x T j   y T j   z T j ] T is the radius-vector of the capturing point of the target and r C i = [ x C i   y C i   z C i ] T is the radius-vector of the capturing system of the chaser.
Note that these coupled equations are expressed in LVLH placed in the target. However, for the onboard algorithms is it necessary to develop the equations in body-fixed reference frame. This issue is addressed in the next chapter.

2.2. Modified Coupled Translational Equations of Motion

As mentioned above ρ ¨ i j is expressed in LVLH and it is required to be expressed in body reference frame of the target. The relation between radius-vector between the centers of mass ρ 0 in the LVLH reference frame and its value ρ ˜ 0 in target body-fixed reference frame is as follows:
ρ ˜ 0 = A ( q T ) ρ 0
where q T is quaternion of the target reference frame relative to LVLH, A ( q T ) is the direction cosine matrix. The kinematical relations for direction cosine matrix can be written as:
A ˙ = Ω A ,
where Ω is skew-symmetric matrix filled using target angular velocity relative to LVLH. Taking time derivative from (17) and substituting (18) one can obtain the following:
ρ ˜ ˙ 0 = Ω ρ ˜ 0 + A ρ ˙ 0 .
From this equation we get the following:
ρ ˙ 0 = A T ( ρ ˜ ˙ 0 Ω ρ ˜ 0 )
Taking derivative from (19) acceleration relationship can be found
ρ ˜ ¨ 0 = Ω ˙ ρ ˜ 0 Ω 2 ρ ˜ 0 + 2 Ω ρ ˜ ˙ 0 + A ρ ¨ 0
On the other hand we can rewrite Clohessy-Whiltshire equations [33] in vector-matrix form as:
ρ ¨ 0 + M ρ ˙ 0 + N ρ 0 = a + K
where
N = [ 3 ω O T 2 ω ˙ O T 0 ω ˙ O T 0 0 0 0 ω O T 2 ] M = [ 0 2 ω O T 0 2 ω O T 0 0 0 0 0 ]
The Equation (22) is written in LVLH. Rewrite both sides of the equation is in target body-fixed reference frame:
A T [ ρ ˜ ¨ 0 Ω ˙ ρ ˜ ˙ 0 + Ω 2 ρ ˜ 0 2 Ω ρ ˜ ˙ 0 ] + M A T ( ρ ˜ ˙ 0 Ω ρ ˜ 0 ) + N A T ρ ˜ 0 = A T a ˜ + A T K ˜
where K ˜ = A K , a ˜ = A a . From this equation one can rewrite it as follows:
ρ ˜ ¨ 0 = a ˜ + K ˜ + β ρ ˜ ˙ 0 + γ ρ ˜ 0
where
β = A [ M A T 2 A T Ω A Ω ˙ ] γ = A [ N A M A T Ω + A Ω 2 ]
So, the Equation (25) are expressed in body-fixed reference frame of target. Now it is possible to derive coupled equations of motion between two points expressed in the target body frame. Let the points of the target r T j and of the chaser r C i are written in target reference frame. Then, using (10) and substituting ρ ˜ 0 ,   ρ ˜ ˙ 0 ,   ρ ˜ ¨ 0 instead of ρ 0 , ρ ˙ 0 , ρ ¨ 0 in (12) the following equation is obtained:
ρ ˜ ¨ i j = β ρ ˜ ˙ i j + γ ρ ˜ i j + ( [ ω T ] x 2 β [ ω T ] x + [ ω ˙ T ] x ) r C i + γ r T j + a ˜ + K ˜
These equations of the relative motion of two arbitrary points of the target and the chaser is used for control algorithm development.

2.3. SDRE Control Algorithm

The common nonlinear motion equations of the considered system is as follows:
x ˙ = f ( x ( t ) ) + g ( x ( t ) , u ( t ) )
Here x n is the state vector; u m is the control vector, f (   ) , g (   ) are nonlinear smooth functions; for establishing the SDRE control algorithm for the dynamical system of type (28) the functional J of type (29) is considered:
J = 1 2 0 t f [ x ( t ) T Q   x ( t ) + u ( t ) T R   u ( t ) ] d t
where Q , R are positive definite weighting matrices. In (29) finite horizon t f is considered. Here the point x = 0 is assumed to be equilibrium of the system. In (29) the minimization of the infinite horizon cost function is required. SDRE method requires the linearization of the equation of motion in a neighborhood of the equilibrium. The optimal coefficients of the regulator are calculated as a result of solving Riccati equation at each time step. In this paper Q , R are considered as constant matrices. Next step is to linearize the nonlinear system. After linearization, the dynamical system has the following form:
x ˙ = A ( x ) x + B ( x ) u
where A ( x ) is the matrix of dynamic and B ( x ) is the control matrix. Similar to linear quadratic regulator a corresponding algebraic Riccati equation can be derived for nonlinear quadratic regulator. After forming the Hamilton function and using the maximum principle of Pontryagin [36], necessary conditions for optimality can be applied to obtain the optimal control law as (31)
u ( x ) = R 1 B T ( x ) P ( x ) x
The control function is similar to LQR, but unlike LQR here the coefficients are functions of state vector. The matrix P ( x ) is unique, symmetric, and positive-definite and can be achieved by solving following Riccati equation:
P ( x ) A ( x ) + A T ( x ) P ( x ) P ( x ) B ( x ) R 1 ( x ) B T ( x ) P ( x ) + Q ( x ) = 0
The closed loop system using this semi-optimal control (31):
x ˙ = ( A B ( x ) R 1 B T ( x ) P ( x ) ) x
The Riccati equation can be solved analytically using Hamiltonian matrix [39,40].

2.4. Control Application to the Problem of Capturing

In this paper the state vector x ( t ) consisting of 12 components is considered. The state vector is composed of three vector-part quaternion components of relative attitude between the target and the chaser [ q 1 , q 2 , q 3 ] T , the angular velocity vector [ ω x , ω y , ω z ] T , relative translational position, and velocity of docking points [ x ˜ i j , y ˜ i j , z ˜ i j , x ˜ ˙ i j , y ˜ ˙ i j , z ˜ ˙ i j ] T . The control vector u ( t ) contains six elements to provide a full feedback control for the coupled motion. The three first elements of vector u ( t ) are reaction wheels control vector h ˙ W C and thrusters force vector F . The matrix of dynamic A ( x ) and control matrix B ( x ) are as follows:
A = [ [ [ D T ( q ) ω T ] x q 4 D T 0 3 × 3 R o ] 0 6 × 6 0 6 × 6 [ 0 3 × 3 I 3 × 3 γ β ] ] .
B = [ 0 3 × 3 0 3 × 3 D I C 1 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 I 3 × 3 / m ] ,
where
R o = [ D I C 1 D T ω T ] I C D T + [ h W C ] x D I C 1 D T [ D I C 1 D T ω T T ] x I C D T [ ω T T ] x + [ I C D T ω T T ] x D I C 1 D T
Since the control is aimed to coincide attitude of the target and the chaser, then the docking points of the target r T j and of the chaser r C i should be chosen in order not to avoid the collision of them. From the mathematical point of view it means the constraint that can be written as scalar product of these two vectors must be negative, i.e., ( r T j , r C i ) < 0 . In the case when the capturing point of the target is not satisfied then we can apply an additional rotation to chaser with attitude matrix S , it means that the point will be calculated as S r C i .

3. Numerical Study

Consider a demonstration of the proposed docking algorithm application. At the initial time the parameters of the target objects is set according to the Table 1. The chaser satellite is initially has no angular velocity with respect to the LVLH, but the target has defined initial rotation. The relative initial conditions for the attitude motion and for the translational motion is also presented in the right part of the Table 1.
In this example the weight matrix Q is set to identity matrix and matrix R is block matrix of the following structure:
R = [ R r o t 0 3 x 3 0 3 x 3 R t r a n ]
The matrix R r o t equals to 160 I 6 x 6 and R t r a n = 160 I 6 x 6 . These parameters are chosen using trial and error approach in order not to meet reaction wheels saturation. It is assumed in the paper that R t r a n is fixed, but the influence of R r o t on the algorithm is studied. The thruster misalignment in this example is set to 3.5 mm. In this example the inertia tensor of the target is set to identity matrix, and the inertia tensor of the chaser is 2 I 3 x 3 . The maximal reaction wheel kinematic momentum is chosen as 1 Nms, which corresponds to reaction wheels mini-wheels produced by Honeywell company [41].
For the docking, it is necessary to control motion in a way that relative distance and velocity (both of rotational and translational motion) converge to zero. These relative position and velocity are shown correspondingly in Figure 2 and Figure 3. This relative translational motion is shown in Figure 4.
The attitude of the chaser and angular velocity with respect to the target is shown in Figure 5 and Figure 6. At the time of about 135 s quaternion parameters achieve values less than 10 5 (about 10 3 deg of the deviation) and it is assumed as the required accuracy for attitude tracking. At the same time the position errors are quite large and the more control effort from thrusters are required to obtain rendezvous. When the relative distance becomes less than 1 cm the docking is assumed to be accomplished. This time is longer and is 294 which is shown in Figure 2 and Figure 3.
Figure 7 shows the discrete thrust values that are used to control translational motion. The value of the discreet thrust is 0.5 N. Figure 8 shows angular momentum of reaction wheels that achieved their final constant values to hold to tracking attitude. Because of thruster misalignment and in order to track the target attitude the reaction wheels angular momentum is increased but still do not exceed the maximal values of 1 Nms.
Reaction wheels have operational limitation, exceeding these limits will lead to saturation, which is caused not only by the thruster misalignment but also by improper choices of control algorithm parameters such as weighting matrix R r o t . The inertia tensor of the target and its rotational velocity can cause saturation as well. Here the task is to determine the maximum values of angular velocity of the target, the misalignment of thrust vector, the weighting matrix at which the capturing maneuver will be successful, i.e., the required attitude and points positions deviation will be achieved. After numerical simulation it will be possible to determine the operational range of reaction wheels under critical variables to avoid the wheels saturation.
Consider a target with inertia tensor that is not identity matrix, but a inertia tensor of body with dynamical symmetry as the following:
I T = d i a g ( J , J , k J )
where J is inertia moment along principal axis, k is the ration coefficient for the inertia moment along axis of dynamical symmetry. According to physical constraint of the inertia moments k > 0.5 . It is assumed that the inertia tensor of the target is known. In practical cases it is estimated at the stage of space debris observation, for example using image processing technique [42]. The target with dynamical symmetry is chosen for consideration because in that case it is easier to study the influence of the inertial moment of the target on control algorithm performance. Moreover, a lot of debris such as upper stages of rockets and some inactive satellites are close to body with dynamical symmetry. Taking all the parameters from the previous simulation example except the target inertia tensor, the maximal value of the reaction wheels momentum is calculated for changing ratio coefficient k . The results are shown in Figure 9. It can be seen that the lowest maximal angular momentum is close to the case when k = 1 that correspond to spherical body. In the case when the body has a dynamical symmetry the more the ratio coefficient of matrix of inertia, the higher the maximal angular momentum of the reaction wheels. It can be explained by the increasing influence of the nutation motion of the body, the capturing point of the target position is harder to track by the reaction wheels of the chaser. To further study the influence of nutation motion on the control algorithm performance, the initial angular velocity is also varied. For the random ratio k from the interval k [ 0.5 , 3 ] and for the random initial angular velocity of the target (each component value is chosen from interval [ 0 , 24 ]   deg / s ) after simulation the maximum value of the reaction wheels momentum is obtained. Figure 10 shows these points for each simulation and corresponding values of maximum angular momentum which is required to be provided by reaction wheels. Almost all of the maximal values of the reaction wheels are less than 1 Nms, that is assumed to be the upper limit. However, in the case the angular velocity is higher that 24 deg/s the wheels saturation often occurs that causes the inappropriate errors for the capturing.
The reaction wheels saturation can be also caused by thruster misalignment. To study its effect on the capturing the random misalignment l [ 0 , 0.01 ] m is added to simulations additionally to the varying parameters k and angular velocity of the target. Moreover, the control parameter weighting matrix R r o t affects the required angular momentum of reaction wheel, consequently the changing elements of R r o t [ 120 , 200 ] are added as well to the simulation alongside other varying parameters. The results of multiple simulations are presented in the Figure 11 in the parameters ranges, where bubbles mean that there is no saturation and the capturing is successful, and the stars are for capturing the failure due to reaction wheels saturation. The smaller the bubbles the less the maximal angular momentum of the reactions wheels.
To interpret Figure 11 it will be helpful to present the box diagram of influential parameters. Figure 12 shows the dependence of maximum angular momentum of reaction wheels to the misalignment of thrusters from the center of mass. The 50% of the simulation results are inside the rectangular and each outside interval contains 25% of results, the mean value is depicted as horizontal line in the box. Figure 12 demonstrates a gradual increase in angular momentum which intuitively could be predicted. Starting from 8 mm of thruster misalignment some of the simulations results exceed the limit of the angular momentum. Figure 13 shows the effect of angular momentum of target on maximum angular momentum of reaction wheels. In the given scale the increment seems that the lower values are exponentially converging to the value of 0.9 Nms. Using Figure 11 with a fixed H T and l it is possible to choose the optimum value for control weighting matrix R r o t .
Thus, the developed relative coupled motion equations and SDRE-based control algorithm allow to study the influence of the most crucial parameters influence on the capturing performance. The multiple numerical study showed that for the presented case with defined parameters the successful capturing is possible when the thruster misalignment does not exceed 8 mm, the control weighting matrix R r o t elements are inside the interval [ 120 , 200 ] . Moreover, the successful capturing is limited to the target angular velocity of about 24 deg/s for each component and inertia tensor ratio k ( 0.5 , 2.8 ] . It should be noted that these results will be different for different initial relative vector and for other limit of the reaction wheels momentum but using the same methodology it is possible to obtain the allowable range for the successful capturing.

4. Conclusions

In this work an active space debris removal approach is considered. An algorithm for capturing a non-cooperative target in LEO orbit is presented. Moreover, a method to explore the boundaries and limitation of this suggested algorithm is developed. This method is used to determine the range of the acceptability of the algorithm with respect to the parameters such as targets angular momentum, and misalignment of the thrusters to avoid saturation of reaction wheels. For a debris with given moment on inertia and angular velocity this method allows to predict the possibility of tracking the target and capturing it. The multiple numerical study showed that for the presented case with defined parameters the successful capturing is possible when the thruster misalignment does not exceed 8 mm, the control weighting matrix R r o t elements are inside the interval [ 120 , 200 ] . Moreover, the successful capturing is limited to the target angular velocity of about 24 deg/s for each component and inertia tensor ratio k ( 0.5 , 2.8 ] . It should be noted that these results will be different for different initial relative vector and for other limit of the reaction wheels momentum but using the same methodology it is possible to obtain the allowable range for the successful capturing.

Author Contributions

M.A. conducted numerical study, D.I. proposed the main theory. All authors have read and agreed to the published version of the manuscript.

Funding

The work is supported by the Russian Foundation of Basic Research, grants NO 18-31-20014, 20-31-90072.

Conflicts of Interest

The authors declare no conflict of interest regarding the publication of this paper.

References

  1. Aglietti, G.S.; Taylor, B.; Fellowes, S.; Salmon, T.; Retat, I.; Hall, A.; Chabot, T.; Pisseloup, A.; Cox, C.; Zarkesh, A.; et al. The active space debris removal mission RemoveDebris. Part 2: In orbit operations. Acta Astronaut. 2020, 168, 310–322. [Google Scholar] [CrossRef] [Green Version]
  2. Hakima, H.; Bazzocchi, M.C.F.; Emami, M.R. A deorbiter CubeSat for active orbital debris removal. Adv. Space Res. 2018, 61, 2377–2392. [Google Scholar] [CrossRef]
  3. Shan, M.; Guo, J.; Gill, E. Review and comparison of active space debris capturing and removal methods. Prog. Aerosp. Sci. 2016, 80, 18–32. [Google Scholar] [CrossRef]
  4. Bonnal, C.; Ruault, J.M.; Desjean, M.C. Active debris removal: Recent progress and current trends. Acta Astronaut. 2013, 85, 51–60. [Google Scholar] [CrossRef]
  5. Mark, C.P.; Kamath, S. Review of Active Space Debris Removal Methods. Space Policy 2019, 47, 194–206. [Google Scholar] [CrossRef]
  6. Kristensen, A.S.; Ulriksen, M.D.; Damkilde, L. Self-Deployable Deorbiting Space Structure for Active Debris Removal. J. Spacecr. Rocket. 2017, 54, 323–326. [Google Scholar] [CrossRef]
  7. Dalla Vedova, F.; Morin, P.; Roux, T.; Brombin, R.; Piccinini, A.; Ramsden, N. Interfacing Sail Modules for Use with “Space Tugs”. Aerospace 2018, 5, 48. [Google Scholar] [CrossRef] [Green Version]
  8. Trofimov, S.; Ovchinnikov, M. Sail-Assisted End-of-Life Disposal of Low-Earth Orbit Satellites. J. Guid. Control Dyn. 2017, 40, 1796–1805. [Google Scholar] [CrossRef]
  9. Smith, B.G.A.; Capon, C.J.; Brown, M.; Boyce, R.R. Ionospheric drag for accelerated deorbit from upper low earth orbit. Acta Astronaut. 2020, 176, 520–530. [Google Scholar] [CrossRef]
  10. Felicetti, L.; Gasbarri, P.; Pisculli, A.; Sabatini, M.; Palmerini, G.B. Design of robotic manipulators for orbit removal of spent launchers’ stages. Acta Astronaut. 2016, 119, 118–130. [Google Scholar] [CrossRef]
  11. Benvenuto, R.; Lavagna, M.; Salvi, S. Multibody dynamics driving GNC and system design in tethered nets for active debris removal. Adv. Space Res. 2016, 58, 45–63. [Google Scholar] [CrossRef]
  12. Botta, E.M.; Miles, C.; Sharf, I. Simulation and tension control of a tether-actuated closing mechanism for net-based capture of space debris. Acta Astronaut. 2020, 174, 347–358. [Google Scholar] [CrossRef]
  13. Dudziak, R.; Tuttle, S.; Barraclough, S. Harpoon technology development for the active removal of space debris. Adv. Space Res. 2015, 56, 509–527. [Google Scholar] [CrossRef]
  14. Aslanov, V.; Yudintsev, V. Dynamics of large space debris removal using tethered space tug. Acta Astronaut. 2013, 91, 149–156. [Google Scholar] [CrossRef]
  15. Zhang, J.; Ye, D.; Biggs, J.D.; Sun, Z. Finite-time relative orbit-attitude tracking control for multi-spacecraft with collision avoidance and changing network topologies. Adv. Space Res. 2019, 63, 1161–1175. [Google Scholar] [CrossRef]
  16. Zhang, J.; Biggs, J.D.; Ye, D.; Sun, Z. Extended-State-Observer-Based Event-Triggered Orbit-Attitude Tracking for Low-Thrust Spacecraft. IEEE Trans. Aerosp. Electron. Syst. 2020, 56, 2872–2883. [Google Scholar] [CrossRef]
  17. Zhang, J.; Biggs, J.D.; Ye, D.; Sun, Z. Finite-time attitude set-point tracking for thrust-vectoring spacecraft rendezvous. Aerosp. Sci. Technol. 2020, 96, 105588. [Google Scholar] [CrossRef]
  18. Hartley, E.N.; Trodden, P.A.; Richards, A.G.; Maciejowski, J.M. Model predictive control system design and implementation for spacecraft rendezvous. Control Eng. Pract. 2012, 20, 695–713. [Google Scholar] [CrossRef] [Green Version]
  19. Breger, L.S.; How, J.P. Safe Trajectories for Autonomous Rendezvous of Spacecraft. J. Guid. Control Dyn. 2008, 31, 1478–1489. [Google Scholar] [CrossRef]
  20. Liu, X.; Lu, P. Solving Nonconvex Optimal Control Problems by Convex Optimization. J. Guid. Control Dyn. 2014, 37, 750–765. [Google Scholar] [CrossRef]
  21. Lu, P.; Liu, X. Autonomous Trajectory Planning for Rendezvous and Proximity Operations by Conic Optimization. J. Guid. Control Dyn. 2013, 36, 375–389. [Google Scholar] [CrossRef]
  22. Ventura, J.; Ciarcià, M.; Romano, M.; Walter, U. Fast and Near-Optimal Guidance for Docking to Uncontrolled Spacecraft. J. Guid. Control Dyn. 2017, 40, 3138–3154. [Google Scholar] [CrossRef]
  23. Sabatini, M.; Palmerini, G.B.; Gasbarri, P. A testbed for visual based navigation and control during space rendezvous operations. Acta Astronaut. 2015, 117, 184–196. [Google Scholar] [CrossRef]
  24. Ivanov, D.; Koptev, M.; Ovchinnikov, M.; Tkachev, S.; Proshunin, N.; Shachkov, M. Flexible microsatellite mock-up docking with non-cooperative target on planar air bearing test bed. Acta Astronaut. 2018, 153, 357–366. [Google Scholar] [CrossRef]
  25. Cloutier, J.R.; Stansbery, D.T. The capabilities and art of state-dependent Riccati equation-based design. In Proceedings of the 2002 American Control Conference (IEEE Cat. No.CH37301), Anchorage, AK, USA, 8–10 May 2002; Volume 1, pp. 86–91. [Google Scholar]
  26. Cloutier, J.R.; Cockburn, J.C. The state-dependent nonlinear regulator with state constraints. Proc. Am. Control Conf. 2001, 1, 390–395. [Google Scholar] [CrossRef]
  27. Çimen, T. Survey of State-Dependent Riccati Equation in Nonlinear Optimal Feedback Control Synthesis. J. Guid. Control Dyn. 2012, 35, 1025–1047. [Google Scholar] [CrossRef]
  28. Navabi, M.; Reza Akhloumadi, M. Nonlinear Optimal Control of Orbital Rendezvous Problem for Circular and Elliptical Target Orbit. Modares Mech. Eng. 2016, 15, 132–142. [Google Scholar]
  29. Felicetti, L.; Palmerini, G.B. A comparison among classical and SDRE techniques in formation flying orbital control. In Proceedings of the 2013 IEEE Aerospace Conference, Big Sky, MT, USA, 2–9 March 2013. [Google Scholar]
  30. Stansbery, D.T.; Cloutier, J.R. Position and attitude control of a spacecraft using the state-dependent Riccati equation technique. Proc. Am. Control Conf. 2000, 3, 1867–1871. [Google Scholar] [CrossRef]
  31. Segal, S.; Gurfil, P. Effect of Kinematic Rotation-Translation Coupling on Relative Spacecraft Translational Dynamics. J. Guid. Control Dyn. 2009, 32, 1045–1050. [Google Scholar] [CrossRef]
  32. Lee, D.; Bang, H.; Butcher, E.A.; Sanyal, A.K. Kinematically Coupled Relative Spacecraft Motion Control Using the State-Dependent Riccati Equation Method. J. Aerosp. Eng. 2015, 28, 04014099. [Google Scholar] [CrossRef]
  33. Navabi, M.; Akhloumadi, M.R. Nonlinear Optimal Control of Relative Rotational and Translational Motion of Spacecraft Rendezvous. J. Aerosp. Eng. 2017, 30, 04017038. [Google Scholar] [CrossRef]
  34. Alfriend, K.T. Spacecraft Formation Flying: Dynamics, Control and Navigation; Elsevier/Butterworth-Heinemann: Oxford, UK, 2010; ISBN 9780750685337. [Google Scholar]
  35. Akhloumadi, M.; Ivanov, D. Satellite relative motion SDRE-based control for capturing a noncooperative tumbling object. In Proceedings of the 9th International Conference on Recent Advances in Space Technologies, RAST 2019, Istanbul, Turkey, 11–14 June 2019; pp. 253–260. [Google Scholar]
  36. Kirk, D.E. Optimal Control Theory: An Introduction; Prentice-Hall: Upper Saddle River, NJ, USA, 1970; ISBN 9780486434841. [Google Scholar]
  37. Massari, M.; Zamaro, M. Application of SDRE technique to orbital and attitude control of spacecraft formation flying. Acta Astronaut. 2014, 94, 409–420. [Google Scholar] [CrossRef]
  38. Clohessy, W.H.; Wiltshire, R.S. Terminal Guidance System for Satellite Rendezvous. J. Astronaut. Sci. 1960, 27, 653–678. [Google Scholar] [CrossRef]
  39. Laub, A.J. Schur Method For Solving Algebraic Riccati Equations. In Proceedings of the IEEE Conference on Decision and Control, San Diego, CA, USA, 10–12 January 1979; pp. 60–65. [Google Scholar]
  40. Lancaster, P.; Rodman, L. Algebraic Riccati Equations; Clarendon Press: Oxford, UK, 1995; ISBN 9780198537953. [Google Scholar]
  41. Miniature Reaction Wheels Interface Connector Mounting Flange. Available online: www.honeywell.com/space/CEM (accessed on 27 October 2020).
  42. Ivanov, D.; Ovchinnikov, M.; Sakovich, M. Relative Pose and Inertia Determination of Unknown Satellite Using Monocular Vision. Int. J. Aerosp. Eng. 2018, 1–16. [Google Scholar] [CrossRef]
Figure 1. Reference frames of chaser satellite and target object.
Figure 1. Reference frames of chaser satellite and target object.
Aerospace 07 00160 g001
Figure 2. Relative position with respect to the target during the capturing maneuver.
Figure 2. Relative position with respect to the target during the capturing maneuver.
Aerospace 07 00160 g002
Figure 3. Relative velocity with respect to the target during the capturing maneuver.
Figure 3. Relative velocity with respect to the target during the capturing maneuver.
Aerospace 07 00160 g003
Figure 4. Relative motion of satellite to target.
Figure 4. Relative motion of satellite to target.
Aerospace 07 00160 g004
Figure 5. Relative quaternion with respect to the target.
Figure 5. Relative quaternion with respect to the target.
Aerospace 07 00160 g005
Figure 6. Relative angular velocity with respect to the target.
Figure 6. Relative angular velocity with respect to the target.
Aerospace 07 00160 g006
Figure 7. Thrust vector of satellite.
Figure 7. Thrust vector of satellite.
Aerospace 07 00160 g007
Figure 8. Angular momentum of reaction wheels.
Figure 8. Angular momentum of reaction wheels.
Aerospace 07 00160 g008
Figure 9. Maximum angular momentum with respect to inertial moment of the target.
Figure 9. Maximum angular momentum with respect to inertial moment of the target.
Aerospace 07 00160 g009
Figure 10. Maximum angular momentum of reaction wheels with respect to the target angular moment value.
Figure 10. Maximum angular momentum of reaction wheels with respect to the target angular moment value.
Aerospace 07 00160 g010
Figure 11. Random points and range of saturation of reaction wheels (bubbles means that there is no saturation and the capturing is successful, and the stars are for the capturing failure due to reaction wheels saturation).
Figure 11. Random points and range of saturation of reaction wheels (bubbles means that there is no saturation and the capturing is successful, and the stars are for the capturing failure due to reaction wheels saturation).
Aerospace 07 00160 g011
Figure 12. Box diagram of eccentricity of thrusters vector and angular momentum of reaction wheels.
Figure 12. Box diagram of eccentricity of thrusters vector and angular momentum of reaction wheels.
Aerospace 07 00160 g012
Figure 13. Box diagram of target’s angular momentum vector and angular momentum of reaction wheels.
Figure 13. Box diagram of target’s angular momentum vector and angular momentum of reaction wheels.
Aerospace 07 00160 g013
Table 1. Parameters and initial condition of the modelling.
Table 1. Parameters and initial condition of the modelling.
Orbital ParametersInitial Conditions
Altitude, km750 q ( t 0 ) = q T T ( t 0 ) = q C C ( t 0 ) [ 0 , 0 , 0 , 1 ] T
Eccentricity0.03 ω ( t 0 ) = ω T T ( t 0 ) ,   deg / s ω C C ( t 0 ) ,   deg / s [ 10 , 10 , 20 ]   T [ 0 , 0 , 0 ]   T
Inclination, deg70 ρ ( t 0 ) = [ x 0 , y 0 , z 0 ] T , m [ 50 , 27 , 100 ] T
Right ascension, deg50 ρ ˙ ( t 0 ) = r ˙ 0 ,   m / s [ 0 , 2 , 0 ] T
Argument of perigee, deg80 r C i ,   m [ 1 ,   1 ,   0 ] T
Initial true anomaly0 r T j , m [ 1 , 0 , 1 ] T
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Akhloumadi, M.; Ivanov, D. Influence of Satellite Motion Control System Parameters on Performance of Space Debris Capturing. Aerospace 2020, 7, 160. https://doi.org/10.3390/aerospace7110160

AMA Style

Akhloumadi M, Ivanov D. Influence of Satellite Motion Control System Parameters on Performance of Space Debris Capturing. Aerospace. 2020; 7(11):160. https://doi.org/10.3390/aerospace7110160

Chicago/Turabian Style

Akhloumadi, Mahdi, and Danil Ivanov. 2020. "Influence of Satellite Motion Control System Parameters on Performance of Space Debris Capturing" Aerospace 7, no. 11: 160. https://doi.org/10.3390/aerospace7110160

APA Style

Akhloumadi, M., & Ivanov, D. (2020). Influence of Satellite Motion Control System Parameters on Performance of Space Debris Capturing. Aerospace, 7(11), 160. https://doi.org/10.3390/aerospace7110160

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop