- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
International Journal of Aerospace Engineering
Volume 2012 (2012), Article ID 231935, 9 pages
Guidance Stabilization of Satellites Using the Geomagnetic Field
1Center for Research and Development in Mathematics and Applications, University of Aveiro, Portugal
2School of Technology and Management, Polytechnic Institute of Viana do Castelo, Avenida do Atlântico, 4900-348 Viana do Castelo, Portugal
Received 21 October 2011; Revised 13 January 2012; Accepted 18 January 2012
Academic Editor: Kenneth M. Sobel
Copyright © 2012 Francisco Miranda. 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.
In the last years the small satellites have played an important role in the technological development. The attractive short period of design and low cost of them and the capacity to solve problems that are usually considered as problems to big and expensive spacecrafts lead us to study the control problem of these satellites. Active three-axis magnetic attitude stabilization of a low Earth orbit satellite is considered in this work. The control is created by interaction between the magnetic moment generated by magnetorquers mounted on the satellite body and the geomagnetic field. This problem is quite complex and difficult to solve. To overcome this difficulty guidance control is considered, where we use ε-strategies introduced by Pontryagin in the frame of differential games theory. Qualitative analysis and results of numerical simulation are presented.
The problem of attitude control of spacecraft has been widely studied in the last years. If the spacecraft is equipped with three independent actuators, a complete solution to the set point and tracking control problems is available. This problem was solved by different ways (see, e.g., [1–4]). If only two independent actuators are available, as discussed in detail in , the problem of attitude regulation is not solvable by means of continuous time-invariant control laws, whereas a time-varying control law, achieving local asymptotic nonexponential stability, was proposed in . Since magnetorquers are relatively reliable, lightweight, and energy efficient, they have become attractive for small and inexpensive satellites, but the above results are not directly applicable if the satellite is equipped with these magnetic coils as actuators. Such actuators operate on the basis of the interaction between the magnetic field generated by magnetic coils installed on the satellite body and the magnetic field of the Earth. This implies a spacecraft control significantly different. There are two classical uses of magnetic torque rods in attitude control. One is for momentum management of wheel-based systems . The other is for angular momentum and nutation control of spinning , momentum-biased , and dual-spin spacecraft . The present study is one of a growing number that considers active three-axis magnetic attitude stabilization of a low Earth orbit satellite considering a nontilted dipole model for the geomagnetic field. Magnetic coils have been extensively used since the early sixties as a simple and reliable technology to implement attitude control actuators in low Earth orbit satellites (see, e.g., [11–16] and the survey  and the references therein). The main difficulty in the implementation of this active three-axis magnetic system consists in the fact that the control torque can only be generated perpendicular to the geomagnetic field vector, which results in a time-varying nonlinear system. To solve this problem the most natural is to use the Lyapunov functions method combined with the Krasovskii-LaSalle theorem . The Lyapunov function is used to construct a stabilizer depending on current time and position of the system. Then, the Krasovskii-LaSalle theorem is applied to prove the asymptotic stability of the equilibrium position. However, a rigorous verification of the Krasovskii-LaSalle theorem conditions faces serious technical difficulties which can hardly be overcome even with the help of systems for symbolic computations. Other techniques were presented (see, e.g., [19, 20]). In , Wiśniewski presented a simple constant gain controller (CGC) via the Linear Quadratic Regulator (LQR) method and proposed other two attitude controllers using the periodic characteristic of the Earth magnetic field: infinite and finite horizon periodic controllers. In , Psiaki showed us how to design a class of stabilizing attitude controllers for nadir-pointing spacecraft using only magnetic actuation. Their control laws are designed using a new type of periodic linear quadratic regulator whose Riccati equation solution is approximated by a linear time-invariant solution for an averaged system. The resulting full-state feedback controller derives its periodicity from the time-varying control influence matrix which can be derived from on-board magnetometer measurements. Both authors considered the linear time-varying system taking the periodic nature of the magnetic field as an advantage to approximate the solution of the Riccati equation. A considerable amount of work has been dedicated to the problems of analysis and design of magnetic control laws in the linear case, that is, control laws for nominal operation of a satellite near its equilibrium attitude (see, e.g., [9, 17, 19–25]). Abdelrahman et al.  proposed an attitude control method of a satellite with only magnetic torquers as actuators via the state-dependent Riccati equation (SDRE) technique based on pseudolinear time-varying modeling of the spacecraft nonlinear dynamics. However, limited attention has been dedicated to global formulations of the problem. In [18, 27, 28], the attitude regulation problem for Earth pointing spacecraft is studied exploiting the (quasi) periodic behavior of the system, which used the standard passivity arguments to prove local asymptotic stabilizability of open-loop equilibria. A solution to the global stabilization problem by means of full (or partial) state feedback has been studied by Lovera and Astolfi [29–32]. Conditions for almost global attitude regulation are obtained in [29, 30] for the case of full state feedback, and local results are derived for the case of output (attitude only) feedback. The paper  shows us how stability conditions similar to those given in [29, 30] that can be derived for control laws achieving Earth pointing for magnetically actuated spacecraft, taking also into account the effect of gravity gradient torques. For this problem, an almost global stabilization result is given for the case of full state feedback, resorting to an adaptive PD-like state feedback control law. With respect to works [29, 30, 32] the results presented in  that can deal with a generic magnetically actuated satellite do not rely on restrictive assumptions on the controller parameters and guarantee that there are no trajectories of the closed-loop system along which the average controllability can be lost. These results do not rely on the (frequently adopted) periodicity assumption for the geomagnetic field along the considered orbit. The proposed framework for closed-loop stability analysis of magnetically controlled spacecraft can be also exploited to predict the effect of actuator faults on the behavior of the controlled satellite.
An alternative way to solve the active three-axis global magnetic attitude stabilization problem of a low Earth orbit satellite is suggested here. The method is based on guidance control , a more general control procedure. Considering an auxiliary global stabilization problem that we know an analytic solution, the method consists to use a trajectory of the globally stabilizable auxiliary system as a guide. We construct a local stabilizer that keeps the trajectory of the system in a neighborhood of a solution of the auxiliary system. In this way, the trajectory of the system tends to the equilibrium position. According to this approach and using -strategies introduced by Pontryagin in the frame of differential games theory , the stabilizing control is constructed as a function of time defined in a small time interval and not as a feedback. From the practical point of view, -strategy is similar to stabilizer which depends on the time and position only, because it usually is implemented as a generator of piecewise constant controls. The numerical implementation of this method is based on the construction of multistep reachability sets . The main contribution of this method is to solve a global stabilization problem applying local stabilization methods that avoid difficulties which we usually have in global stabilization problems, helping to overcome serious mathematical difficulties and being more effective in applications. This procedure also yields additional possibilities for the design of a stabilizer that eliminates the peak effect, that is, the large deviation of the solutions from the equilibrium position at the beginning of the stabilization process that we have when we construct stabilizers to obtain a high speed of damping of the control systems trajectories (studies of this effect are presented in [35–37]). The minimization problem used in this paper also eliminates this effect which represents a serious obstacle to guidance stabilization.
From all studies about the nonlinear controllability theory, the work  presents an application of this theory to the time-varying attitude dynamics of a magnetically actuated spacecraft in a Keplerian orbit in the geomagnetic field. The model considered in  does not include the gravitational torque and therefore the result can be applied to spherically symmetric satellites only. Considering a satellite subjected to gravity gradient torque, the controllability of the magnetic attitude of the satellite is studied here. Conditions of controllability when the satellite circular orbit is neither polar nor equatorial are obtained by local controllability theory.
It is assumed that the satellite moves along a circular orbit and the current value of the geomagnetic field with respect to the satellite axes is known due to three-axis magnetometer measurements. The position of the satellite mass center is also assumed to be known and the currents in the coils are control parameters.
2. Guidance Stabilization Process
Consider a control system where and are sufficiently smooth functions, that is, there exist all derivatives needed in our considerations. Let be such that for all , that is, is an equilibrium position. The global stabilization problem is to find a map such that and the equilibrium position of the differential equation is globally asymptotically stable. This problem is usually very hard to solve. One of the ways to avoid difficulties is to use a guidance control , a more general control procedure. Consider an auxiliary differential equation such that is an globally asymptotically stable equilibrium position. The trajectories of system (1), under the guidance control, are defined as solutions to the Cauchy problem: where the trajectory , that is considered to be a “guide’’, tends to the equilibrium position . Then we construct , , , , such that , , where is a solution to (1). To this end, we linearize system (1) along the trajectory and apply the methods developed to solve the stabilization problem for linear systems. This approach works only if is sufficiently small. For this reason, we construct the control solving the following minimization problem where , , , and . In this way, we generate an -strategy, that is, a sequence of controls defined on the intervals , , and so on, that minimizes the norm of the solution forcing this way the solution to be small for all . This implies that the trajectory of (1) follows the trajectory of the auxiliary differential equation. As , , then . These controls depend on and , . This problem has a solution if the control system , , is locally controllable.
Let , , be a measurable bounded control. The solution of the linearized system , where , is denoted by and the symbol is used for the set of all solutions. Define the reachability set . Using the following well-known result from the control theory (see ), we have a possibility to check if minimization problem (3) has a solution, that is, we can prove that there exists an -strategy such that , .
Theorem 1 (). The linearized system is controllable on the interval , that is, , if and only if, only the trivial solution to the adjoint differential equation , , satisfies the orthogonality condition , . The condition implies that control system (1) is controllable around the trajectory on the interval .
3. Application to the Satellite Stabilization Problem
Consider a satellite moving along a circular orbit. Introduce two Cartesian reference systems and . The system is the body reference system. The origin of this system, , coincides with the satellite mass center and the axes are directed along the principal inertia axes. The system is the orbital reference system. The origin coincides with the center of mass of the satellite. The axis is directed along the radius vector of the satellite mass center and the axis is perpendicular to the orbital plane. The angular position of the satellite with respect to the orbital system is described by three positive left hand rotations defined by the Euler angles . Coordinates of a vector in the satellite body system, , and in the orbital system, , satisfy the relation , where is an orthogonal 3-1-2 rotation matrix with . The components of the matrix: are Denote by , , the rows of the matrix . Obviously .
The satellite is equipped with three magnetic coils oriented along three orthogonal axes. Let be the vector of angular velocity in the body axes. The rotation of the satellite subjected to gravity gradient torque is described by the following equations: where(i)the tensor of inertia is a diagonal matrix with the diagonal elements , , , such that ,(ii)the vector of the geomagnetic field in the body axes is given by , where ,(iii) is the orbital inclination, that is, the angle between the equatorial and orbital planes,(iv) is the argument of latitude,(v)−1 is the angular velocity of the orbital motion,(vi)m is the radius of the orbit,(vii)Am2 is the Earth magnetic dipole moment,(viii)m3/c2 is the Earth gravitational parameter,(ix)H/m is the magnetic parameter,(x) describe the initial position of the mass center in the orbit,(xi)the vector has the components , , where stands for the current in the th coil, is the number of turns, and is the area of a loop.
The currents , , are control parameters. At the initial moment , the satellite has an angular position that corresponds to the matrix and the angular velocity . The problem is to find control laws that drive and angular velocity to and , respectively. Considering the Lyapunov function, known as the Jacobi integral, we see that this is a stable equilibrium position [18, 28].
Consider the auxiliary system: where for , , , , , , and has the same properties as .
Using the following lemma, the proof of Theorem 2 is trivial.
Lemma 3. The system where , has the asymptotically stable equilibrium position , , .
Proof. Considering the Lyapunov function and following the proof of theorem in , we get the result.
We can also see the proof of this lemma in .
Proof of Theorem 2. Using the equality: in (8), we have where Considering and using Lemma 3, we see that system (12)–(13) has the following asymptotically stable equilibrium position Using again equality (11) in (14)-(15), we obtain the result.
Auxiliary system (8)–(9) will be a guide for control system (6). Using the techniques presented in previous section, we numerically create a control such that system (6) follows auxiliary system (8)–(9). Finding , we get the currents , . Its only necessity to prove that the linear system of minimization problem (3) is locally controllable.
Considering , the corresponding adjoint system takes the form: where , , , , , If , then the first equation of the adjoint system implies , because . The orthogonality condition from Theorem 1 can be written in two equivalent forms and . Differentiating the latter, we obtain . Substituting this for in the first equation of the adjoint system, we have . Differentiating the first equation of the adjoint system and substituting the expression for obtained above, we get . Differentiating twice, we get . Substituting the above formulas for and its derivatives, we obtain , where are real analytic functions. From Theorem 1 we see that to prove the local controllability of the linearized system of (6) it suffices to show that on any interval of time. This happens, for example, if the vectors and are linearly independent. Thus we have the following result.
Lemma 4. If for almost all , then the linearized system of (6) is controllable.
Assume that the satellite circular orbit is neither polar nor equatorial. In this case, with the help of the symbolic computation system Mathematica, we prove that the conditions of Lemma 4 are satisfied for and is sufficiently small and for and in a neighborhood of the stable equilibrium position. See the two following examples to each one of these cases.
Example 5. Consider , and . Using the symbolic computation system Mathematica we obtain As for all , then the vectors and are linearly independent. This proves the conditions of Lemma 4.
Example 6. Now, consider any , and . Using the symbolic computation system Mathematica we obtain , where with Suppose that . As for all , then we have . From triangular inequality for moments of inertia, we have . Thus is equivalent to and is equivalent to From (24) and (25) we get Dividing both sides of (26) by , we have that is a contradiction. Thus we have to a finite number of points . This proves the conditions of Lemma 4.
4. Stabilization Algorithm and Numerical Simulation
Consider the following stabilization algorithm based on the previous results and a numerical simulation.
4.1. Stabilization Algorithm
The stabilization algorithm presented (see Algorithm 1) is a simple illustration of the previous method that we can implement to generate a sequence of controls defined on the intervals , , and so on. This sequence of controls can be used in a real situation in real time.
4.2. Numerical Simulation
To test the algorithm we consider satellite (6) and auxiliary system (8)–(9) having the same initial position in and , where and consider the following parameters: the orbital inclination , the tensor of inertia , the initial position of the mass center in the orbit , , , , , , and .
Implementing an -strategy in the form of three-step reachability set construction algorithm, that is, considering in the above algorithm, we obtain the results of the numerical simulation that are shown in Figures 1 and 2. This -strategy consists in the following. First, we divide the intervals of time , , and so on, on three subsets with the same size , , and using the Euler’s formula, we obtain approximations to which depend of piecewise constant controls . Finally, using an numerical algorithm to calculate the minimum of, we obtain these piecewise constant controls. Then we have and therefore we find the control laws , , that drive system (6) to the equilibrium position. Figure 1 represents the norm of the difference between the trajectory of satellite (6) and the trajectory of auxiliary system (8)–(9), and Figure 2 represents the norm of the difference between the trajectory of satellite (6) and the equilibrium position.
The problem of active three-axis magnetic attitude stabilization of a low Earth orbit satellite was studied. In this paper it was proposed an alternative approach to the existing methods. The method used in the stabilization of satellite is based on guidance control, a more general control procedure. Conditions of controllability were also obtained in this study. According to these conditions and using -strategies, the stabilizing control was constructed as a function of time defined in a small time interval and not as a feedback. The simulation results obtained using the algorithm proposed prove the feasibility of the approach presented here. This simple approach using local stabilization methods solves the complex and difficult global stabilization problem.
|:||The set of -by- real matrices|
|:||The set of real numbers|
|:||The usual -dimensional space of vectors , where ,|
|:||The Jacobian matrix of the function|
|:||The inner product of two vectors and in|
|:||The vector product of two 3-dimensional vectors and|
|:||The Euclidean norm of a vector|
|:||The matrix -Norm of a real matrix , ,|
|:||The unit ball in|
|:||The identity matrix|
|:||The transposed matrix of a real matrix|
|:||The inverse matrix of a real matrix|
|:||The interior of a set .|
This work was supported by the Portuguese Foundation for Science and Technology (FCT), the Portuguese Operational Programme for Competitiveness Factors (COMPETE), the Portuguese Strategic Reference Framework (QREN), and the European Regional Development Fund (FEDER).
- J. T.-Y. Wen and K. Kreutz-Delgado, “The attitude control problem,” IEEE Transactions on Automatic Control, vol. 36, no. 10, pp. 1148–1162, 1991.
- O.-E. Fjellstad and T. I. Fossen, “Comments on “The attitude control problem”,” IEEE Transactions on Automatic Control, vol. 39, no. 3, pp. 699–700, 1994.
- F. Caccavale and L. Villani, “Output feedback control for attitude tracking,” Systems and Control Letters, vol. 38, no. 2, pp. 91–98, 1999.
- M. R. Akella, “Rigid body attitude tracking without angular velocity feedback,” Systems and Control Letters, vol. 42, no. 4, pp. 321–326, 2001.
- C. I. Byrnes and A. Isidori, “On the attitude stabilization of rigid spacecraft,” Automatica, vol. 27, no. 1, pp. 87–95, 1991.
- P. Morin, C. Samson, J.-B. Pomet, and Z.-P. Jiang, “Time-varying feedback stabilization of the attitude of a rigid spacecraft with two controls,” Systems and Control Letters, vol. 25, no. 5, pp. 375–385, 1995.
- H. Iida and K. Ninomiya, “A new approach to magnetic angular momentum management for large scientific satellites,” NEC Research and Development, vol. 37, no. 1, pp. 60–77, 1996.
- L. D. D. Ferreira and J. J. Cruz, “Attitude and spin rate control of a spinning satellite using geomagnetic field,” Journal of Guidance, Control, and Dynamics, vol. 14, no. 1, pp. 216–218, 1991.
- M. E. Pittelkau, “Optimal periodic control for spacecraft pointing and attitude determination,” Journal of Guidance, Control, and Dynamics, vol. 16, no. 6, pp. 1078–1084, 1993.
- K. T. Alfriend, “Magnetic attitude control system for dual-spin satellites,” AIAA Journal, vol. 13, no. 6, pp. 817–822, 1975.
- J. S. White, F. H. Shigemoto, and K. Bourquin, “Satellite attitude control utilizing the Earth’s magnetic field,” Tech. Rep. D-1068, NASA, Washington, DC, USA, 1961.
- A. C. Stickler and K. T. Alfriend, “Elementary magnetic attitude control system,” Journal of Spacecraft and Rockets, vol. 13, no. 5, pp. 282–287, 1976.
- F. Martel, P. K. Pal, and M. L. Psiaki, “Active magnetic control system for gravity gradient stabilized spacecraft,” in Proceedings of the 2nd Annual AIAA/USU Conference on Small Satellites, Logan, Utah, USA, September 1988.
- K. L. Musser and W. L. Ebert, “Autonomous spacecraft attitude control using magnetic torquing only,” in Proceedings of the Flight Mechanics/Estimation Theory Symposium, pp. 23–38, NASA, May 1989.
- R. Wiśniewski, “Nonlinear control for satellite detumbling based on magnetic torquing,” in Proceedings of the 22nd (JSDE) Joint Service Data Exchange for Guidance, Navigation and Control, Scottsdale, Ariz, USA, October-November 1994.
- R. Wiśniewski and M. Blanke, “Three-axis satellite attitude control based on magnetic torquing,” in Proceedings of the 13th IFAC World Congress, San Francisco, Calif, USA, June 1996.
- E. Silani and M. Lovera, “Magnetic spacecraft attitude control: a survey and some new results,” Control Engineering Practice, vol. 13, no. 3, pp. 357–371, 2005.
- R. Wiśniewski and M. Blanke, “Fully magnetic attitude control for spacecraft subject to gravity gradient,” Automatica, vol. 35, no. 7, pp. 1201–1214, 1999.
- R. Wiśniewski, “Linear time-varying approach to satellite attitude control using only electromagnetic actuation,” Journal of Guidance, Control, and Dynamics, vol. 23, no. 4, pp. 640–647, 2000.
- M. L. Psiaki, “Magnetic torquer attitude control via asymptotic periodic linear quadratic regulation,” Journal of Guidance, Control, and Dynamics, vol. 24, no. 2, pp. 386–394, 2001.
- R. Wiśniewski and F. L. Markley, “Optimal magnetic attitude control,” in Proceedings of the 14th IFAC World Congress, Beijing, China, July 1999.
- M. Lovera, E. De Marchi, and S. Bittanti, “Periodic attitude control techniques for small satellites with magnetic actuators,” IEEE Transactions on Control Systems Technology, vol. 10, no. 1, pp. 90–95, 2002.
- H. Yan, I. M. Ross, and K. T. Alfriend, “Pseudospectral feedback control for three-axis magnetic attitude stabilization in elliptic orbits,” Journal of Guidance, Control, and Dynamics, vol. 30, no. 4, pp. 1107–1115, 2007.
- M. Corno and M. Lovera, “Spacecraft attitude dynamics and control in the presence of large magnetic residuals,” Control Engineering Practice, vol. 17, no. 4, pp. 456–468, 2009.
- T. Pulecchi, M. Lovera, and A. Varga, “Optimal discrete-time design of three-axis magnetic attitude control laws,” IEEE Transactions on Control Systems Technology, vol. 18, no. 3, pp. 714–722, 2010.
- M. Abdelrahman, I. Chang, and S.-Y. Park, “Magnetic torque attitude control of a satellite using the state-dependent Riccati equation technique,” International Journal of Non-Linear Mechanics, vol. 46, no. 5, pp. 758–771, 2011.
- C. Arduini and P. Baiocco, “Active magnetic damping attitude control for gravity gradient stabilized spacecraft,” Journal of Guidance, Control, and Dynamics, vol. 20, no. 1, pp. 117–122, 1997.
- C. J. Damaren, “Comments on “Fully magnetic attitude control for spacecraft subject to gravity gradient”,” Automatica, vol. 38, no. 12, p. 2189, 2002.
- M. Lovera and A. Astolfi, “Global attitude regulation using magnetic control,” in Proceedings of the 40th IEEE Conference on Decision and Control (CDC '01), pp. 4604–4609, Orlando, FL, USA, December 2001.
- M. Lovera and A. Astolfi, “Spacecraft attitude control using magnetic actuators,” Automatica, vol. 40, no. 8, pp. 1405–1414, 2004.
- M. Lovera and A. Astolfi, “Global magnetic attitude control of inertially pointing spacecraft,” Journal of Guidance, Control, and Dynamics, vol. 28, no. 5, pp. 1065–1067, 2005.
- M. Lovera and A. Astolfi, “Global magnetic attitude control of spacecraft in the presence of gravity gradient,” IEEE Transactions on Aerospace and Electronic Systems, vol. 42, no. 3, pp. 796–805, 2006.
- V. Bushenkov and G. Smirnov, Stabilization Problems with Constraints: Analysis and Computational Aspects, Gordon and Breach Science Publishers, Amsterdam, The Netherlands, 1997.
- L. S. Pontryagin, “Linear differential games. I,” Soviet Mathematics. Doklady, vol. 8, pp. 769–771, 1967.
- R. N. Izmailov, “The peak effect in stationary linear systems with scalar inputs and outputs,” Automation and Remote Control, vol. 48, no. 8, pp. 1018–1024, 1987.
- H. J. Sussmann and P. V. Kokotovic, “The peaking phenomenon and the global stabilization of nonlinear systems,” IEEE Transactions on Automatic Control, vol. 36, no. 4, pp. 424–440, 1991.
- G. Smirnov, V. Bushenkov, and F. Miranda, “Advances on the transient growth quantification in linear control systems,” International Journal of Applied Mathematics and Statistics, vol. 14, no. J09, pp. 82–92, 2009.
- S. P. Bhat, “Controllability of nonlinear time-varying systems: applications to spacecraft attitude control using magnetic actuation,” IEEE Transactions on Automatic Control, vol. 50, no. 11, pp. 1725–1735, 2005.
- G. Smirnov, Introduction to the Theory of Differential Inclusions, American Mathematical Society, Providence, RI, USA, 2002.
- I. V. Burkov, “Asymptotic stabilization of the position of a rigid body with fixed point without velocity measurements,” Systems and Control Letters, vol. 25, no. 3, pp. 205–209, 1995.
- S. Sastry, Nonlinear Systems: Analysis, Stability, and Control, Interdisciplinary Applied Mathematics, Vol. 10, Springer, Berlin, Germany, 1999.