Satellite Attitude Control System Design considering the Fuel Slosh Dynamics

The design of the satellite attitude control system (ACS) becomes more complex when the satellite structure has different type of components like, flexible solar panels, antennas, mechanical manipulators, and tanks with fuel. A crucial interaction can occur between the fuel slosh motion and the satellite rigid motion during translational and/or rotational manoeuvre since these interactions can change the satellite centre of mass position damaging the ACS pointing accuracy. Although, a well-designed controller can suppress such disturbances quickly, the controller error pointing may be limited by the minimum time necessary to suppress such disturbances thus affecting the satellite attitude acquisition. As a result, the design of the satellite controller needs to explore the limits between the conflicting requirements of performance and robustness. This paper investigates the effects of the interaction between the liquid motion (slosh) and the satellite dynamics in order to predict what the damage to the controller performance and robustness is. The fuel slosh dynamics is modelled by a pendulum which parameters are identified using the Kalman filter technique. This information is used to design the satellite controller by the linear quadratic regulator (LQR) and linear quadratic Gaussian (LQG) methods to perform a planar manoeuvre assuming thrusters are actuators.


Introduction
The problem of interaction between fluid and structure is important when one needs to study the dynamic behavior of offshore and marine structures and road and railroad containers partially filled with a fluid [1,2].In space missions the sloshing problem appears when the spacecraft is spinning and there is liquid inside it.An example is the damp device involving fluid as the damping material, whose motion can interact with flexible panels and/or manipulators [3].An interesting approach to analyze a rigid container mounted on flexible springs interacting with a perfect fluid including sloshing effects has been done by Lui and Lou [4].The success of a space mission can depend on taking into account the knowledge of the interaction between fluid motion (slosh) and structure dynamics since this interaction can damage the ACS pointing requirements.A space structure, like rockets, geosynchronous satellites, and the space station, usually contains liquid in tanks that can represent more than 40% of the initial mass of the system.As a result, the first step to design its ACS is to obtain a detailed dynamics model of the space structure.When the fuel tanks are only partially filled and suffer a transverse acceleration and/or rotational motion, large quantities of fuel move uncontrollably inside the tanks and generate the sloshing effects.Agrawal [5] has shown that the dynamics interaction between the fuel motion and the rigid and/or flexible body dynamics can result in some kind of control instability.For minimizing these effects the ACS must be designed using a robust control method in order to assure stability and good performance to achieve the attitude control system requirement [6].When a rigid-flexible satellite with fuel tanks inside is subjected to large angle manoeuvre its dynamics is only captured by complex nonlinear mathematical model.Besides, the remaining flexible and/or liquid vibration can introduce a tracking error resulting in a minimum attitude acquisition time.Souza [7] has done a detailed investigation of the influence of the nonlinearities introduced by the panel's flexibility into the ACS design.It was shown that system parameters variation can degrade the control system performance, indicating the necessity to improve 2 Shock and Vibration the ACS robustness.An experimental controller robustness and performance investigation was done by Conti and Souza [8], where the estimation of the platform inertia parameters was introduced as part of the platform ACS design.The problem of designing satellite nonlinear controller for rigid satellite has been done by Souza and Gonzales [9] using the state-dependent Riccati equation (SDRE) method which is able to deal with high nonlinear plants.Due to the complexity of modeling the fluid and/or flexible dynamic of the system it is common to use mechanical systems analogies that describe this dynamic.Besides, if one needs to know some physical parameters related with the slosh or the flexibility dynamics it is common to obtain them by experimental apparatus or some kind of estimating method such as Kalman filter [10].

Satellite with Sloshing Model
The phenomenon of sloshing is due to the movement of a free surface of a liquid that partially fills a compartment and this movement is oscillating.It depends on the shape of the tank, the acceleration of the gravity, and on the axial/rotational acceleration of the tank.As a representative of the behavior of the total weight of the system it is accepted that when the mass of the liquid oscillates the mass center of the rigid body also oscillates, thereby disturbing the rigid-flexible part of the vehicle under consideration.As an oscillating movement it is natural to consider the wave generated by the movement of the liquid as a stationary wave in all oscillation modes.Each mode of oscillation has a special feature of this phenomenon under study and one observes, in a quantitative sense, how much mass is displaced.Among all the modes that cause the greatest disruption in the system are the first and the second modes.Despite the fact that oscillation has lower frequency it is capable of resulting in violent shifting of the center of mass of the liquid creating an oscillation in the system as a role.The other oscillation modes act as a less aggressive and it may not even vary the position of its center of mass due to the symmetry of the wave which on average causes no displacement.Due to its complexity, the sloshing dynamics is usually represented by mechanical equivalents that describe and reproduce faithfully the actions and reactions due to forces and torques acting on the system.The main advantage of replacing the fluid model with an equivalent oscillating model [11] is simplifying the analysis of motion in the rigid body dynamics compared to the fluid dynamics equations.Due to the complexity of establishing an analytical model for the fluid moving freely within a closed tank, a simplified system is used, taking into account the following criteria [6]: (a) small displacements, (b) a rigid tank, and (c) no viscous, incompressible, and homogeneous liquid.Under these conditions the dynamics of the sloshing can be approximated by mechanical system consisting of a mass-spring or pendulum.Consider a rigid spacecraft moving in a fixed plane with a spherical fuel tank and including the lowest frequency slosh mode.Based on the Lagrange equation and the Rayleigh dissipation function one can model systems using the mechanical mass-spring and pendulum type system, respectively.Figure 1 shows a satellite model where slosh dynamics is represented by its pendulum analogous mechanical system, where the mass of the satellite and the moment of inertia, regardless of the fuel, are given by  and , respectively, and the mass equivalent of fuel and its inertia moment are given by   and   , respectively.The attitude control of the spacecraft is done by the force  and by the pitch moment .The constant thrust  is responsible for the orbital transfer of the spacecraft with respect to the inertial reference system (, , ) and it acts on the center of mass of spacecraft in the longitudinal axis.Also it is given the velocity of the center of the fuel tank   ,   and the attitude angle  of the spacecraft with respect to a fixed reference (, , ).Besides,  represents the linear velocity and  represents the angular velocity of the rigid body.The length of the pendulum is , the distance from satellite center of mass to the pendulum connected point is , and the angle of the pendulum with respect to the spacecraft longitudinal axis is , which is assumed in the equilibrium position ( = 0) about the reference axis .

The Satellite Equations of Motion
The satellite equations of motion can be derived using the Lagrange equations [10] given by d dt Details of the equations of motion derivation can be found in [12], where  is the Lagrangian of the system, the generalized coordinates are  and ,  is the Rayleigh dissipation function,   is the internal torque, and   is the external torque.Assume that ,   ,   , ,  are given by The position vector of the satellite mass center with respect to the inertial reference system (, , ) (see Figure 1) is Assuming the relations V  = ẋ + θ and V  = ż − θ the satellite velocity is given by The position of the mass of fuel with respect to the inertial reference system (, , ) (see Figure 1) is given by As a result, the mass of the fuel velocity is The Lagrangian of the entire system is given by Substituting ( 4) and ( 6) into (7), using the relations given by (2), and performing the derivations of (1), one obtains the satellite equations of motion given by Assuming the relations All equations derived previously are nonlinear.However, in order to design LQR and LQG controllers and estimate the sloshing parameters using the Kaman filter technique one has to get the linear set of equations of motion, which is obtained assuming that the system makes small movements around the point of equilibrium very close to zero; that is, (, θ , , ψ ) = (0, 0, 0, 0).Now, substituting (12) into ( 10) and (11) and assuming the linearization conditions, one has the satellite equation of motion given by where  * =   /( +   ),  * =   /( +   ), and  *  =   /( +   ).

Kalman Filter (KF)
The Kalman filter estimates the instantaneous state of a dynamic system from the reading of measurements.In other words the Kalman filter is the optimal solution of minimum variance, which means that the equations of the Kalman filter can be deduced from this premise [7].The Kalman filtering can be divided into two steps: time update and measurement update.
4.1.Time Update.This step propagates the states and the covariance of the time  −1 to   .For this it just integrates the following equations, with the boundary conditions x −1 = _ x −1 and P −1 = P−1 .Herein P is the covariance matrix and x is the vector state: where  is a nonlinear vector function of the state x and time .
where G is a matrix that adds noise to the system dynamic, Q is the process noise covariance, and F is the Jacobian matrix of .

Measurement
Update.This step updates the state and covariance to the time  due to the measurement y.Consider the following: where   is the filter gain,   is the Jacobian matrix of ℎ which in turn is a vectorial nonlinear function of the state, and y is the measurement vector.

Linear Quadratic Gaussian (LQG)
This method is, basically, the union with the LQR [13] and the Kalman filter.In the LQR method one assumes that all states are available to be feedback but in reality that is not true.Therefore, when there is any state that is not available one has to use the Kalman filter to estimate it in order to be a feedback.The separation principle ensures that each step of this process can be made independently of each other; one may first solve the LQR problem and then design the optimal estimator (Kalman filter), or vice versa, so that the global solution is always the same.Assume a plant described by the linear state equations given by ẋ () = Ax () + Bu () + Γ, where x represents the state vector, the matrix A is the state matrix, B is the input matrix, y is the output vector, C is the output matrix,  and  are white noise, and u is the control law.In the LQG method [14] the control law gain is obtained by the LQR method and it is given by where R is real symmetric positive definite matrix and P  is the symmetrical solution of the LQR Riccati equation given by Similarly the Kalman filter gain now is given by where V is real symmetric positive definite matrix and P  is the symmetrical solution matrix of the KF Riccati equation given by where P  = P   ≥ 0 and P  = P   ≥ 0 and Q, R, V, and W are weight matrices which can be regarded as setting parameters (tuning) that must be manipulated until they find one acceptable response to the system.
A necessary and sufficient condition to guarantee the existence of the K  and K  is if the system is completely controllable and observable.The LQG method is more realistic than the LQR method, since the states are not always available to be feedback or need to be measurable with the aid of sensors.The inclusion of the noise in the model which represents imperfections of the system is also one advantage of the LQG method.

Simulations and Results
The first simulation is the design of the control law using the LQR and LQG control theories, for the spacecraft with a partially filled tank, to account for the sloshing dynamics by the mechanically analog pendulum type.The pendulum physical parameters used in the simulations are as follows:  = 600 kg,   = 100 kg,  = 720 kg/m 2 ,   = 90 kg/m 2 ,  = 0.3 m,  = 0.3 m,  = 500 N, and  = 0.19 kg m 2 /s.The equations of motion that describe the dynamics of the system are given by ( 13) which need to be put in state space form.These equations describe the angular displacement of the spacecraft and the angular displacement of the pendulum and the initial conditions used are  = 2 ∘ , θ = 0.57 ∘ /s,  = 1 ∘ , and ψ = 0 ∘ .
Figures 2 and 3 show that the LQR control law performance is better than the LQG and the reason is because the first one considers that the sloshing variables are available to be feedback which is not true.
The performance of the LQG control law is damaged because the sloshing motion is controlled indirectly and    the sloshing state variables are not available to be feedback and they need to be estimated by the Kalman filter.Now one uses the Kalman filter to estimate the rod length of the pendulum using a database containing the time evolution of the variables  and θ .The rod length  is now considered as state variable so one has a system with five state variables which can be estimated.Using the same control LQR law, it is possible to analyze the control system performance during the estimation of the rod length.Besides, one inserts a white noise in the measurements  and θ , given by the noise process ]  = (0, 0.001).
Figure 4 shows that the estimated value of the rod length is  = 0.33 m, remembering that the "actual" value of the rod is 0.3 m. Figure 5 shows that plant of the system depends on the estimation of the rod, since at the time that the rod length is varying the position of the poles also moves.
Figures 6 and 7 show the LQR and the LQG controllers during the rod length estimation.One observes that the system is controlled in less than approximately 15 seconds and that the slosh motion is more oscillatory with the LQG controller.This is because the sloshing motion is controlled indirectly and the sloshing state variables are being estimated by the Kalman filter.

Figure 1 :
Figure 1: Satellite model with slosh dynamics pendulum analogous mechanical system.

Figure 2 :
Figure 2: Performance of the LQR control law.

Figure 3 :
Figure 3: Performance of the LQG control law.

Figure 4 :
Figure 4: The estimation of the rod length.

Figure 5 :
Figure 5: Poles position during the estimation.