Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 257193, 8 pages

http://dx.doi.org/10.1155/2013/257193

## Design of the Microsatellite Attitude Control System Using the Mixed Method via LMI Optimization

Instituto Nacional de Pesquisas Espaciais (INPE), Avenida dos Astronautas 1758, 12227-101 So Jos Dos Campos, SP, Brazil

Received 5 April 2013; Revised 4 July 2013; Accepted 11 July 2013

Academic Editor: Maria Zanardi

Copyright © 2013 Erberson Rodrigues Pinheiro and Luiz Carlos Gadelha de Souza. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Due to the space missions limited budget, small satellite cluster or constellation would be an economical choice. From risk-sharing viewpoint, a number of smaller satellites have a significant reliability advantage over a bigger one. Generally, one satellite is subject to two types of uncertainties: structured uncertainty that represents some satellite parameter variation and the unstructured uncertainty that represents some kind of the satellite model error. On the other hand, the Satellite Attitude Control (SAC) design becomes more vulnerable to uncertainty disturbances like model error and moment-of-inertia variation as the satellite has great decrease in size and weight. This is the case for a microsatellite with mass less than 100 kg where the ACS performance and robustness become very sensitive to both kinds of uncertainties. As a result, the design of the SAC has to deal with both types of uncertainties which is associated with the drawback between controller performance and robustness. The purpose of this work is to model a microsatellite taking into account the uncertainties and to perform the Control System Design based on the mixed methodology via LMI optimization.

#### 1. Introduction

Microsatellites play an important role in space missions, such as position location, Earth observation, atmospheric data collection, space science, and communication. Some spacecrafts used to observation need high-accuracy performance on pointing requirement, so it is necessary to apply a three-axis attitude control, leading a multivariable control system [1]. In the face of disturbance and uncertainty, it is necessary to design a robust control for analysis and synthesis of attitude control system. Examples of satellite robust control system design using multiobjective and nonlinear approaches can be found in [2, 3], respectively. Low orbit spacecrafts are under a more strong influence of gravity gradient torque, aerodynamic torque, and magnetic torque. Some equipments on the microsatellite like cameras, telescopes, and solar array can move causing change on moment of inertia. Microsatellites with mass less than 100 kg are more sensitive to moment of inertia variation and disturbances like external torques [4]. In this work, we will be using a kind of robust control called mixed control. This combination was introduced by Bernstein and Haddad [4]; where the idea was to minimize an norm of a transfer function subjected to a constraint given by a norm of another transfer function. In the paper [5] was considered the state and the output feedback of the mixed control to solve the non-linear Riccati equation in a convex optimization context. As for vibration control of rigid-flexible satellite, an alternative approach is to use piezoelectric shunt damping technique as has been done in [6]. In this work, one uses the mixed control via the LMI approach [7] to design an attitude control of a microsatellite subjected to an external disturbances and with uncertainty in the moment of inertia.

#### 2. Microsatellite Attitude Dynamics

It is defined as a body-fixed reference frame with its origin located in the center of mass of a microsatellite and is given the unit vectors being along the principal axes. The Euler equations of a microsatellite are given by [7] where , , and are the principal moments of inertia, , , and are the body-axis components of angular velocity, , , and are the nonmodelled external torques, , , and are the components of gravity gradient torques that will be inserted into the equations, and , , and are the control torques.

It is necessary to consider another reference system A local-vertical local-horizontal (LVLH) with its origin at the center of mass of the microsatellite. The LVLH frame has the following unitary vectors , with in the direction of the microsatellite velocity in the orbital plane, pointing to the Earth, and normal to the orbit plane.

To describe the orientation of the body-fixed frame with respect of LVLH frame in terms of Euler angles, the following coordinate transformation is used:

The angular velocity of the body-fixed frame relative to the LVLH is given by The angular velocity of the body fixed frame relative to the inertial frame fixed in the Earth center becomes where is the orbital frequency of the microsatellite.

For small attitude deviation from LVLH orientation, the following linearized attitude dynamics can be obtained as follows:

##### 2.1. Gravity Gradient Torque

In the space, the gravitational field is not uniform, so the variation in the gravitational field over the body yields the gravitational torque of the center of mass of the body. On the assumption that the microsatellite center of mass is in a Keplerian circular orbit and the Earth is spherical, the gravity gradient torque along the body axes becomes [8] Inserting the gravity gradient equations (7) into the Euler equations (1) and making a linearization, one has These equations are the Euler equations for the microsatellite, from which one observes that the pitch axis is decoupled from the row and yaw axes.

#### 3. Structured Uncertainty

To represent the system uncertainty it will be used Linear Fractional Transformation (LFT) [9]. As shown in Figure 1, using the LFT procedure the block transfer function from the perturbation signal to error signal is given by

The plant of the system can be represented by the block which is given by Considering that there is uncertainty in the principal moment of inertia of the microsatellite, it can be expressed as a nominal value plus a perturbation [7] given by where is the variation and is the normalized uncertainty. Inserting it into the Euler equations, one has the dynamic equation with uncertainty. This uncertainty can be pulled out of the system, and it can be considered as a disturbance.

Let us do this calculation, initially, for pitch axis of the microsatellite which is decouple, is the disturbing torque. As result, the equation of motion is given by This equation of motion can be put in the block diagram as shown in Figure 2.

Using the Linear Fractional Transformation (LFT), the first block of the Figure 2 can be represented by Comparing this equation with the LFT (9), the first block will be Doing the same procedure for the second and the third block is given by As a result, the new block diagram for the pitch axis taking into account the uncertainty is as shown in Figure 3.

This representation helps to understand how the uncertainty acts in the system and how it can be lumped out of the system like a perturbation. Usually, the uncertainty is incorporated in the generalized plant in the diagonal form as shown in Figure 4.

The generalized plant represented in Figure 4 can be given by

Taking into account all derivations performed up to now and the relationship between the input and the output of the microsatellite system, one obtains the generalized plant in matrix form given by

The mixed design aprroach consists first to minimizing the perturbation effect of moment of inertia uncertainty (structured) by the norm. The closed loop system will be robust and stable for any uncertainty satisfying the relation . Specifically the lower the value of ; the closed loop system remains stable for a large uncertainty, that is, a large range of inertia moment variation.

On the other hand, the external disturbance uncertainty (unstructured) will be minimized by the norm. As a result, the closed loop system remains stable for external perturbation, which is associated with good performance, for example, quick time response and small overshot.

Therefore, in order to include both kinds of uncertainties so as the controller designed presents good robustness and adequate performance, the new generalized plant must have both signals that will be minimized. Here, the generalized plant in matrix form is given by

In that case, one has the inputs , the outputs , the states , and the perturbation . As a result, the new state space model that includes both kinds of design requirements is given by

#### 4. The Mixed Controller Theory

Figure 5 shows the block diagram of the mixed control approach [5], where the above block represents the uncertainty; the medium block is the generalized plant, and the below block is the controller to designed.

The mixed controller design is a multiobjective control problem where the goal is to minimize the norm in order to improve performance subjected to the minimization of norm to guaranty robustness requirement [9] which can be expressed by

As shown in [8], the mixed control problem is equivalent to minimizing Trace Tr(Q) over the matrices , , and satisfying the Linear Matrix Inequalities [9] given by

Assuming that the LMIs (22) and (23) have solutions , , and , the mixed controller is given by which results in the following two expressions:

#### 5. Simulations and Results

The information about the microsatellite used in the simulation is given by Tables 1 and 2

In the mixed controller design, a key point is to find the appropriate values of . Besides, one must keep in mind that for one has the pure control problem, and for one has the pure control problem. Here, just for simulations propose one decides to begin using to perform a comparative study. As for the uncertainty, one assumes that the variation on the moment of inertia is about 10%. In order to obtain the maximum and the minimum uncertainty variations, one considers two kinds of plants given by Variation in the inertia moment for pitch angle: Plant with uncertainty 1: , and . Plant with uncertainty 2: , and . Variation in the inertia moment for pitch angle for roll and yaw angles: Plant with uncertainty 1: , and . Plant with uncertainty 2: , and . Initial conditions:

In Figures 6, 7, and 8 the dashed line represents the plant with uncertainty variation, and the continues line represents the nominal plant without uncertainty. They show the Euler angles (roll, pitch, and yaw) control for pure control (-blue line); the controller is very robust with respect to uncertainty, because there is no difference between the nominal plant and the plant with uncertainty. On the other hand, it is noted that the difference between the nominal plant and the plant with uncertainty increases for pure control (-red line).

Figures 9, 10, and 11 show that for pure control (blue line), the controller signal has bigger overshot than for pure control (red line), which represents the drawback between its robustness and performance. Considering that the microsatellite actuator must have small torque, the mixed control with is a good choice to design the controller, because it is not so slow like pure , and the control signal is not so strong like pure control.

Spacecraft is subjected to small disturbances on the space, and these disturbances can be persistent. In the case of a low orbit, the microsatellite is more subject to disturbances due to the gravity gradient, magnetic, and aerodynamic torques. The gravity gradient torque is included into the equations, the magnetic torque is cyclic and can be approximated by sinusoids with different frequencies, and the aerodynamic torque is cumulative and can be approximated to a step. As a result, one has assumed that these torques can be represented by the following equation:

Figures 12, 13, and 14 show the Euler angles (roll, pitch, yaw) for the pure control (blue line), the mixed control (-black line), and the pure control (red line). One observes that the pure control has the best capacity of attenuation with respect to a sinusoidal disturbance. The pure control is not robust with respect to a sinusoidal disturbance. In order to have a good balance control between robustness and performance one must choose some values of such that the mixed control provide a good performance even with the perturbation of the external torques. Again, the best mixed control controller value for is 2.

#### 6. Conclusion

This paper presents a microsatellite model taking into account the uncertainties and the design of the satellite control system based on the mixed methodology via LMI optimization. This control technique is used to design the microsatellite attitude control system in the face of environmental disturbance (unstructured uncertainty) and moment of inertia variation (structured uncertainty). It is well know that the controller provides robust stability with respect to structured uncertainty while the controller provides good performance with respect to unstructured uncertainty. Here, one investigates the conjunction of both methods in order to improve the performance and robustness of the SAC system. To do this, one assumes that the microsatellite is subjected to uncertainty in the moment of inertia variation of about 10%, and environmental disturbances were approximated to sinusoidal function plus a step function. The simulations have shown that the controller has presented the best robustness and performance than the controller with respect to uncertainty due to inertia moment variation and due to external disturbance. However, in all simulations the controller signal was bigger than the controller, which can cause bigger overshot and can saturate the actuator, once the microsatellite usually needs a small actuator. As a result, the way to achieve robustness stability and good performance was to design the controller using the mixed control, because in this procedure one can choose an adequate value for the tuning parameter so as one can have robust control and with low control signal.

#### References

- X. C. Méndez Cubillos and L. C. G. de Souza, “Using of H-infinity control method in attitude control system of rigid-flexible satellite,”
*Mathematical Problems in Engineering*, vol. 2009, Article ID 173145, 9 pages, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - I. Mainenti, L. C. G. Souza, and F. L. Souza, “Design of a nonlinear controller for a rigid-flexible satellite using multi-objective generalized extremal optimization with real codification,”
*Shock and Vibration*, vol. 19, pp. 1–10, 2012. View at Google Scholar - L. C. G. Souza and R. G. Gonzalez, “Application of the state-dependent riccati equation and kalman filter techniques to the design of a satellite control system,”
*Shock and Vibration*, vol. 19, pp. 22–28, 2012. View at Google Scholar - D. S. Bernstein and W. M. Haddad, “LQG control with an
*H*_{∞}performance bound: a Riccati equation approach,”*IEEE Transactions on Automatic Control*, vol. 34, no. 3, pp. 293–305, 1989. View at Publisher · View at Google Scholar · View at Scopus - P. P. Khargonekar and M. A. Rotea, “Mixed
*H*_{2}/*H*_{∞}control: a convex optimization approach,”*IEEE Transactions on Automatic Control*, vol. 36, no. 7, pp. 824–837, 1991. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - T. P. Sales, D. A. Rade, and L. C. G. Souza, “Passive vibration control of flexible spacecraft using shunted piezoelectric transducers,”
*Aerospace Science and Technology*, vol. 1, pp. 12–26, 2013. View at Publisher · View at Google Scholar - C.-D. Yang and Y.-P. Sun, “Mixed
*H*_{2}/*H*_{∞}state-feedback design for microsatellite attitude control,”*Control Engineering Practice*, vol. 10, no. 9, pp. 951–970, 2002. View at Publisher · View at Google Scholar · View at Scopus - J. R. Wertz and W. J. Larson,
*Space Mission Analysis and Design*, Microcosm Press, Hawthorne, Calif, USA, 1989. - S. Boyd, L. El Ghaoui, E. Feron, and V. Balakrishnan,
*Linear Matrix Inequalities in System and Control Theory*, SIAM Studies in Applied Mathematics, SIAM, Philadelphia, Pa, USA, 1994. View at MathSciNet