Table of Contents Author Guidelines Submit a Manuscript
Mathematical Problems in Engineering
Volume 2013, Article ID 483457, 7 pages
Research Article

Averaging Methods for Design of Spacecraft Hysteresis Damper

1School of Technology and Management of Lamego, Avenida Visconde Guedes Teixeira, 5100-074 Lamego, Portugal
2Centre for Aerospace Science and Technologies, University of Beira Interior, Calçada Fonte do Lameiro, 6201-001 Covilhã, Portugal
3Scientific Area of Mathematics, ESTGV, Polytechnic Institute of Viseu, Campus Politécnico, 3504-510 Viseu, Portugal
4Centre of Physics, Department of Mathematics and Applications, University of Minho, Campus de Gualtar, 4710-057 Braga, Portugal

Received 29 April 2013; Accepted 22 May 2013

Academic Editor: Antonio F. Bertachini A. Prado

Copyright © 2013 Ricardo Gama et al. 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.


This work deals with averaging methods for dynamics of attitude stabilization systems. The operation of passive gravity-gradient attitude stabilization systems involving hysteresis rods is described by discontinuous differential equations. We apply recently developed averaging techniques for discontinuous system in order to simplify its analysis and to perform parameter optimization. The results obtained using this analytic method are compared with those of numerical optimization.

1. Introduction

Dampers that use magnetic hysteresis rods to dissipate the energy of undesired angular motions occurred during deployment or caused by perturbations are used in attitude control systems of small satellites since 1960s [1]. Mathematical modeling of such systems is quite a difficult task since the majority of existent hysteresis models result in differential equations with discontinuous right-hand side.

Analysis of dynamics for attitude control systems with magnetic hysteresis dampers and optimization of their parameters have been done in [2, 3], and the results of these studies have been implemented in real missions [4, 5]. However, these studies lack an accurate theoretical basis for application of averaging methods to such problems.

Recently, an adequate mathematical approach has been developed by the authors in [6]. Now we can address a complete mathematical theory for attitude stabilization systems with hysteresis.

Consider a differential equation describing a mechanical system with stabilizer. Here is a parameter. It is assumed that for all and ; that is, the velocity of the system near the origin is small. Here we do not assume that zero is an equilibrium position of system (1). The parameter should be chosen to optimize, in some sense, the behavior of the trajectories. The choice of this parameter can be based on various criteria. Obviously, it is impossible to construct a stabilizer optimal in all aspects. Consider, for example, a linear controllable system. The pole assignment theorem guarantees the existence of a linear feedback yielding a linear differential equation with any given set of eigenvalues, so one can choose a stabilizer with a very high damping speed. However, such a stabilizer is practically useless because of the so-called peak effect (see [7, 8]). Namely, there exists a large deviation of the solutions from the equilibrium position at the beginning of the stabilization process, whenever the module of the eigenvalues is big.

The aim of this paper is to develop effective analytical and numerical tools oriented to optimization of stabilizer parameters for passive attitude stabilization system with hysteresis rods.

Throughout this paper, we denote the set of real numbers by and the usual -dimensional space of vectors with components in by . We denote by the usual scalar product in and by a norm. By we denote the closed unit ball, that is, the set of vectors satisfying . The transpose of a matrix is denoted by . The set of positively definite symmetric -matrices is denoted by . If and are two subsets in and , we use the following notations: , . The convex hull and the closure of a subset are denoted by and , respectively. The Hausdorff distance between two sets is defined as We denote by the set of solutions to the Cauchy problem ,  ,  , and by the reachability set. We use also the notations and . The closed unit ball in the space of continuous functions with the uniform norm, , is denoted by . The set of locally integrable functions is denoted by . The upper limit of a set-valued map is given by where stands for the graph of the set valued map .

2. Statement of the Problem

Consider dynamics of a satellite in a circular geocentric orbit. The satellite is equipped with a gravity-gradient attitude control system that includes a number of magnetic hysteresis rods as a damper. The spacecraft's equations of motion [2] can be represented in the normalized form. Denote by the solution to the respective Cauchy problem where is a parameter from a compact set . Define the functions

Here are compact sets, are norms in , and . Consider the following mathematical programming problem:

Many problems of stabilization systems' parameters optimization can be written in this form (see [9]). For example, the minimization of the final deviation can be formalized as follows: and the minimization of the maximal deviation of trajectories satisfying certain restrictions at the final moment of time has the form:

3. Averaging for Discontinuous Systems

The averaging method is one of the most used methods to analyze differential equations of the form appearing in the study of nonlinear problems. The idea behind the averaging method is to replace the original equation by the averaged one:

This equation is simpler and has solutions close to the solutions of the original equation. A rigorous justification of the method is given by Bogolyubov's first theorem containing an estimate for the distance between the solutions of the exact and averaged systems on large time intervals [10]. The Samoilenko-Stanzhitskii theorem [11, Theorem  2] which is a generalization of Bogolyubov's second theorem shows that asymptotic stability of the zero equilibrium position of averaged (10) implies that the solutions to original (9) are close to zero on the infinite time interval.

For several models of systems with hysteresis, including the passive attitude stabilization systems, the function appearing in (9) is discontinuous (see, e.g., [3]), and the classical notion of solution and the classical averaging method cannot be used. For such systems Filippov proposed a generalized concept of solution, rewriting problem (9) as a differential inclusion where is an upper semicontinuous set-valued map obtained from by Filippov regularization [12, 13]. The use of this concept of solution makes it necessary to generalize the averaging method to differential inclusions. Many results extending Bogolyubov's first theorem to differential inclusions have been obtained (see, e.g., [14, 15]). In the case of Lipschitzian differential inclusions, the problem has been completely solved by Plotnikov [14]. Averaging results for inclusions with upper semicontinuous right-hand side have been obtained by Plotnikov [15] under conditions of Lipschitz continuity of the averaged inclusion and for inclusions with a piecewise Lipschitzian right-hand side. Recently [6], an averaged differential inclusion has been introduced allowing one to prove extensions of Bogolyubov's first theorem and of the Samoilenko-Stanzhitskii theorem for upper semi-continuous differential inclusions and, as a consequence, for discontinuous dynamical systems. Here we outline the main results from [6].

Let be a set-valued map. Set

We denote by the convex hull of the map and define the averaged differential inclusion as

Note that under Lipschitz condition this map coincides with if the limit exists in the sense of Hausdorff distance (see [6]).

Assume that the following conditions are satisfied:(C1), for all ; (C2) the set-valued map is upper semi-continuous; (C3) for any there exists measurable selection of , that is, there exists such that is measurable for all ; (C4) there exists a nonnegative such that for all ; (C5) there exists the limit

Under these conditions, the following version of Bogolyubov's first theorem is true.

Theorem 1. Let and let be a set-valued map satisfying conditions (C1)–(C5). Let be a compact set. Then for any there exists such that for any and any solution , there exists a solution satisfying


Next theorem is an extension of the Samoilenko-Stanzhitskii theorem.

Theorem 2. Let be a set-valued map satisfying conditions (C1)–(C5). Assume that is an asymptotically stable equilibrium position of the differential inclusion . Then for any , there exist and such that , whenever .

The last theorem shows that if the averaged inclusion has zero as its asymptotically stable equilibrium position, the trajectories of the original inclusion stay in the vicinity of the origin provided and are sufficiently small.

If the averaged inclusion has a special form, we can go further and make some conclusion on the detailed behaviour of the trajectories of the original system. Assume that the averaged inclusion has the form where , , the real parts of the matrix eigenvalues are negative for all , and the function is continuous for all . If is sufficiently small, then the set of solutions to the Lyapunov inequality (see [16]) for the matrix , is nonempty and compact for all . Let . Denote by the Euclidean norm defined by . There exist positive constants and satisfying whenever for some .

Theorem 3. Let ,  , and . There exists such that for all the condition implies the inequality .

This theorem shows that the behavior of the trajectory can be characterized in terms of the pair . The parameter is responsible for the damping speed of the process, while the form of the ellipsoid describes the amplitude of the deviation of the trajectory from the origin. The aim of parameter choosing can be formulated as follows: maximal value of and maximal sphericity of the ellipsoid . The latter property guarantees minimal overshooting of the damping process and, as a consequence, the largest region of applicability of the approximation obtained via averaging.

4. Choosing Passive Magnetic Stabilizer Parameters

The in-plane oscillations of a satellite moving along a polar circular orbit and equipped with a passive gravity-gradient attitude stabilization system with one hysteresis rod are described by the equation where is the pitch angle of spacecraft, is small parameter proportional to the rod's volume, and the force is given by

Here describes the projection of the geomagnetic field on the rod axis, is the hysteresis function, corresponds to the coercive force, and is the argument of latitude of the current point of the orbit [2]. The vector describes the orientation of the hysteresis rod in the satellite body. Equation (22) is equivalent to the system

After the change of variables one arrives at the system

Theorem 4. Assume that is an irrational number. Then the averaged system for (31) is where

, and .

Obviously we have Therefore, we see that the linearization of the averaged system always has a Lyapunov function , and the damping speed is determined by the value of . This means that the peak effect does not take place for the linearization of the averaged system.

To maximize the damping, one has to increase the total volume of the hysteresis material on board. However, it is wellknown that the efficiency of a damping rod is increased with the increase of the ratio between the rod's length and its cross-section dimension. Therefore, instead of one massive bar, the attitude control system should use several rather thin rods of the maximum length allowed by the spacecraft geometrical and system restrictions. On the other hand, to minimize the perturbation of the spacecraft angular motion caused by the damping system itself, the direction of total magnetic field in the rods should deviate as little as possible from the direction of the geomagnetic field at the current point of the orbit. Thus, in general case, one should use a system of three equal orthogonal hysteresis rods or a number of such systems. Here we consider in-plane satellite dynamics on a polar orbit, and for such purpose it suffices to analyze a pair of equal orthogonal rods.

Orientation of this pair of equal orthogonal rods can be defined as , , where and , and due to the system symmetry it is enough to study the interval . If the satellite is equipped with several identical hysteresis rods, the corresponding nonlinear system is

Here the terms , describe the interaction of the respective rod with the geomagnetic field. For a couple of equal orthogonal rods and for small deviation from the origin, the forces and are given by respectively. Here

For this case, the first approximation of the averaged system takes the form where

An easy calculation shows that the optimal value of is , so . In the next section we numerically analyze the validity of the previous analytical study.

5. Numerical Simulations

We approximate problem (6) by the following problem: where ,  ,  ,  , and is the Euler approximation for the solution . Note that this is a hard problem because of the discontinuity of the system. It is necessary to consider very fine partition of the time interval in order to get a good approximation of the solutions to the discontinuous differential equation. For smooth right-hand sides the number of points in the mesh can be significantly reduced (see [9]).

Let be small enough. We consider a satellite with two equal orthogonal hysteresis rods. A typical trajectory of the system is shown in Figure 1.

Figure 1: A typical trajectory for , , .

Note that the oscillations of the trajectory do not allow one to characterize its damping speed using the norm at the final moment of time. For this reason, we numerically solve the following problem: where is an interval corresponding to the period of oscillation of solutions. The minimization is done using multistart Nelder-Mead method. For , , , , , and , the fastest damping speed is observed for values of very close to , that is, to the expected value. Next we consider the problem

The results of minimization for , , , , , , and show that the problem has many local minima but with the values of the functional very close to ; that is, the nonlinear system with two orthogonal hysteresis rods has no overshooting for all values of . This allows one to conclude that the value is the best one. It guarantees high damping speed and does not cause peaking.

6. Conclusions

This paper is dedicated to the problem of parameter optimization for a gravitationally stabilized satellite with magnetic hysteresis damper. Its motion is described by differential equations with discontinuous right-hand side. The discontinuity is the principal obstacle in the application of the averaging method. Our recently obtained results on averaging of discontinuous systems are applied now to rigorously justify the use of this method for a satellite with hysteresis rods. We consider here the simplest case of in-plane oscillations on a polar circular orbit. Theorem 3 shows that the behavior of the system can be characterized in terms of Lyapunov function which can be chosen in order to guarantee the best properties of damping process. Further study is performed numerically, and the simulations are in excellent agreement with the analytical results, confirming also previous studies on hysteresis damping of satellite pitch oscillation in gravity-gradient mode.

Our results can also be applied to rigorously justify the use of averaging techniques in analysis of other engineering problems involving differential equations with discontinuous right-hand side.


This appendix contains the proofs of Theorems 3 and 4.

Proof of Theorem 3. Set . Then, using (21), from the Lyapunov inequality we get , whenever . There exists such that . Therefore, we have and . This ends the proof.

Proof of Theorem 4. Consider the function where
First, note that for any fixed pair the function is analytic. Therefore the integral (see Section 2) is a point. This implies that the averaged operator defined in (14) coincides with if the limit exists. Since is irrational number and can be considered as a -periodic function of the arguments and , we see that limit (A.3) does exist, and we have To evaluate this integral, we represent the derivative in the form where
Thus, we have
Observe that
Substituting (A.7) into (A.4) and integrating in order of and then in order of , we obtain the result.


This research is supported by the Portuguese Foundation for Science and Technologies (FCT), the Portuguese Operational Programme for Competitiveness Factors (COMPETE), the Portuguese Strategic Reference Framework (QREN), and the European Regional Development Fund (FEDER).


  1. R. E. Fischell and F. F. Mobley, “A system for passive gravity-gradient stabilization of earth satellites,” Guidance and Control, vol. 2, pp. 37–71, 1964. View at Google Scholar
  2. V. A. Sarychev, V. I. Penkov, M. Yu. Ovchinnikov, and A. D. Guerman, “Motions of a gravitationally stabilized satellite with hysteresis rods in a polar orbit,” Cosmic Research, vol. 26, no. 5, pp. 561–574, 1988. View at Google Scholar
  3. A. D. Guerman, M. Yu. Ovchinnikov, V. I. Pen'kov, and V. A. Sarychev, “Non-resonant motions of a satellite with hysteresis rods under conditions of gravity orientation,” Mechanics of Solids, vol. 24, pp. 1–11, 1989. View at Google Scholar
  4. M. Ovchinnikov, V. Pen'kov, O. Norberg, and S. Barabash, “Attitude control system for the first sweedish nanosatellite MUNIN,” Acta Astronautica, vol. 46, no. 2–6, pp. 319–326, 2000. View at Publisher · View at Google Scholar · View at Scopus
  5. F. Santoni and M. Zelli, “Passive magnetic attitude stabilization of the UNISAT-4 microsatellite,” Acta Astronautica, vol. 65, no. 5-6, pp. 792–803, 2009. View at Publisher · View at Google Scholar · View at Scopus
  6. R. Gama, A. Guerman, and G. Smirnov, “On the asymptotic stability of discontinuous systems analysed via the averaging method,” Nonlinear Analysis: Theory, Methods and Applications A, vol. 74, no. 4, pp. 1513–1522, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  7. R. N. Izmailov, “The peak effect in stationary linear systems with scalar inputs and outputs,” Automation and Remote Control, vol. 48, no. 8, part 1, pp. 1018–1024, 1987. View at Google Scholar · View at Scopus
  8. H. J. Sussmann and P. V. Kokotović, “The peaking phenomenon and the global stabilization of nonlinear systems,” IEEE Transactions on Automatic Control, vol. 36, no. 4, pp. 424–439, 1991. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  9. A. Guerman, A. Seabra, and G. Smirnov, “Optimization of parameters of asymptotically stable systems,” Mathematical Problems in Engineering, vol. 2011, Article ID 526167, 19 pages, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  10. N. N. Bogoliubov and Y. A. Mitropolsky, Asymptotic Methods in the Theory of Non-Linear Oscillations, Gordon and Breach, New York, NY, USA, 1961. View at MathSciNet
  11. A. M. Samoilenko and A. N. Stanzhitskii, “On averaging differential equations on an infinite interval,” Differential Equations, vol. 42, no. 4, pp. 476–482, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  12. A. F. Filippov, Differential Equations with Discontinuous Right-Hand Sides, vol. 18 of Mathematics and Its Applications, Kluwer Academic, New York, NY, USA, 1988. View at Zentralblatt MATH · View at MathSciNet
  13. G. V. Smirnov, Introduction to the Theory of Differential Inclusions, vol. 41 of AMS Graduate Studies in Mathematics, American Mathematical Society, Providence, RI, USA, 2002. View at MathSciNet
  14. V. A. Plotnikov, “Averaging method for differential inclusions and its application to optimal control problems,” Differential Equations, vol. 15, pp. 1013–1018, 1979. View at Google Scholar · View at MathSciNet
  15. V. A. Plotnikov, “Asymptotic methods in the theory of differential equations with discontinuous and multi-valued right-hand sides,” Ukranian Mathematical Journal, vol. 48, no. 11, pp. 1605–1616, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet
  16. A. M. Liapunov, Stability of Motion, Academic Press, New York, NY, USA, 1966. View at MathSciNet