## ON THE CONTROLLABILITY OF AN ORBITING SATELLITE MODEL WITH ELECTROMAGNETIC-ONLY ACTUATION

## Ye. Yevgenieva

Institute of Applied Mathematics and Mechanics, Sloviansk, Ukraine, Max Planck Institute for Dynamics of Complex Technical Systems, Magdeburg, Germany yevgeniia.yevgenieva@gmail.com

## A. Zuyev

Institute of Applied Mathematics and Mechanics, Sloviansk, Ukraine, Max Planck Institute for Dynamics of Complex Technical Systems, Magdeburg, Germany alexander.zuyev@gmail.com

## Ju. Kalosha

Institute of Applied Mathematics and Mechanics, Sloviansk, Ukraine, julykucher@gmail.com

July 2024

THIS PAPER IS DEDICATED TO THE 100 TH ANNIVERSARY OF PROFESSOR PAVEL V. KHARLAMOV.

## ABSTRACT

This paper presents sufficient conditions for small-time local controllability of a control-affine system that describes the rotational motion of a satellite in a circular orbit. The satellite is modeled as a rigid body subject to electromagnetic actuation. We focus on the underactuated scenario where the control torque is generated solely by magnetorquers. The main contributions of this work include proving small-time local controllability around the relative equilibrium under some natural assumptions on the mass distribution of the rigid body. This result is based on the Lie algebra rank condition and Sussmann's controllability condition. Furthermore, it is shown that the linearized system is not controllable in a neighborhood of the considered equilibrium.

MSC: 93B05, 93C10, 93C15, 70Q05.

Keywords: nonlinear control system · satellite with electromagnetic actuation · controllability · Lie brackets

## 1 Introduction

One of the fundamental problems in control theory is controllability: determining whether it is possible to drive the system from one state to another using a specified class of controls. Despite significant advances in this area, the study of controllability for general nonlinear systems remains a challenging problem [1, 10, 15]. Powerful methods for the controllability analysis of essentially nonlinear systems are extensively covered in [8]. A classical approach relies on the Lie algebraic methods [3]. For driftless control-affine systems, the controllabity conditions are completely characterized by the Chow-Rashevskii theorem [7, 21]. For wide classes of nonlinear systems with drift, some important sufficient controllability conditions have been developed in [4, 12-14, 17, 22-26].

This paper is devoted to the mathematical model of a rotating satellite with electromagnetic actuation. To our knowledge, there are no explicit controllability conditions available for this model as of now. It is mentioned in [28] that

for a satellite controlled by a set of magnetorquers, magnetic torquing lacks controllability in the direction of the local geomagnetic field. The controllability problem for a linear model of a spacecraft actuated by magnetic torques is addressed in [29].

In [3, 9], the Lie algebraic methods have been applied to describe reachability sets and controllability conditions for rotating rigid bodies. It is known that these kinds of systems also present significant challenges to the stabilization problem. Notably, it was shown in [5] that a mathematical model of an underactuated rigid spacecraft does not satisfy Brockett's necessary stabilizability condition.

The present work is motivated by the study [20] on a satellite model with electromagnetic-only actuation in a timevarying geomagnetic field and by our previous paper [16] on a satellite with an attached flexible boom. In this paper, we adapt the model considered in [16] to the case of rigid-only satellite and consider an underactuated situation where the control torque is generated by magnetorquers. For this mathematical model, we will show that the linear approximation is not controllable and prove nonlinear local controllability around a relative equilibrium under certain assumptions on the mass distribution.

## 2 Mathematical model

Consider a satellite which moves as a rigid body in a circular orbit around the Earth. We fix a right Cartesian frame Oxyz (body coordinate system - BCS) with unit vectors ( e , e 1 2 , e 3 ) assuming that O is the center of mass of the satellite. To describe the attitude of the satellite, we also consider the orbit coordinate system (OCS) as the right Cartesian frame Ox y z ′ ′ ′ with unit vectors ( η , η 1 2 , η 3 ) such that η 1 is orthogonal to the orbital plane, η 2 is tangent to the orbit, and η 3 points at the zenith (i.e. away from the Earth's center).

/negationslash

We denote by ω = ω e 1 1 + ω e 2 2 + ω e 3 3 the absolute angular velocity vector of the satellite and by ω 0 = const = 0 the orbital rate. Then, the absolute angular velocity of the OCS is -ω η 0 1 , and the relative angular velocity of the satellite with respect to the OCS is ω r = ω + ω η 0 1 . We represent the rotation of the BCS with respect to the OCS by the attutude quaternion with the vector part q = ( q 1 , q 2 , q 3 ) and the scalar part q 4 . The equations of motion of the considered model are described by Euler's equations with respect to the angular velocity components and by the kinematic equations in quaternion form [16, 19]:

$$\overset { \cdot } { I } \dot { \omega } + \omega \times \underset { \cdot } { I } \omega = \tau _ { c } + \underset { \cdot } { 3 } \omega _ { 0 } ^ { 2 } ( \eta _ { 3 } \times I \eta _ { 3 } ),$$

$$\dot { \ q } = \frac { 1 } { 2 } q _ { 4 } \omega _ { r } + \frac { 1 } { 2 } q \times \omega _ { r },$$

$$\overline { 2 }$$

$$\begin{smallmatrix} \psi & 2 ^ { \psi _ { 4 } \widetilde { \ } r \cdot } & 2 ^ { \psi _ { 4 } \widetilde { \ } r \cdot } \\ \dot { q } _ { 4 } = - \frac { 1 } { 2 } \left \langle q, \omega _ { r } \right \rangle. \end{smallmatrix}$$

We assume that Oxyz are the principal axes of inertia of the satellite, so that the inertia tensor is diagonal: I = diag ( I 1 , I 2 , I 3 ) . The term τ c in (1) is the torque generated by magnetorquers according to the following law:

/negationslash

$$\tau _ { c } = \beta u \times \eta _ { 1 }, \quad \beta = \text{const} \neq 0,$$

where u ∈ R 3 is the control input (magnetic moment generated by the satellite), and βη 1 is the Earth magnetic field vector. Formula (3) corresponds to the case of an equatorial orbit under the assumption of a constant geomagnetic field. The last term in the right-hand side of (1) describes the gravity gradient torque [27].

For a component-wise representation of the above equations, we express the unit vectors of the OCS by their coordinates in the BCS [20]:

$$\eta _ { 1 } = \begin{pmatrix} q _ { 1 } ^ { 2 } - q _ { 2 } ^ { 2 } - q _ { 3 } ^ { 2 } + q _ { 4 } ^ { 2 } \\ 2 ( q _ { 1 } q _ { 2 } - q _ { 3 } q _ { 4 } ) \\ 2 ( q _ { 1 } q _ { 3 } + q _ { 2 } q _ { 4 } ) \end{pmatrix}, \, \eta _ { 2 } = \begin{pmatrix} 2 ( q _ { 1 } q _ { 2 } + q _ { 3 } q _ { 4 } ) \\ q _ { 4 } ^ { 2 } - q _ { 1 } ^ { 2 } + q _ { 2 } ^ { 2 } - q _ { 3 } ^ { 2 } \\ 2 ( q _ { 2 } q _ { 3 } - q _ { 1 } q _ { 4 } ) \end{pmatrix}, \, \eta _ { 3 } = \begin{pmatrix} 2 ( q _ { 1 } q _ { 3 } - q _ { 2 } q _ { 4 } ) \\ 2 ( q _ { 2 } q _ { 3 } + q _ { 1 } q _ { 4 } ) \\ q _ { 3 } ^ { 2 } + q _ { 4 } ^ { 2 } - q _ { 1 } ^ { 2 } - q _ { 2 } ^ { 2 } \end{pmatrix}.$$

Formulas (3) and (4) allow presenting differential equations (1) and (2) with respect to ω t ( ) and q t ( ) as the following control system:

$$& \text{formulas} \left ( 3 \right ) \text{ and} \left ( 4 \right ) \text{allow presenting differentual equations} \left ( 1 \right ) \text{ and} \left ( 2 \right ) \text{with respect to} \, \omega \left ( t \right ) \text{ and} \, q \left ( t \right ) \text{ as the following} \\ & \text{col} \text{ system:} \\ & \dot { x } = f ( x, u ), \quad x = \left ( \omega _ { 1 }, \omega _ { 2 }, \omega _ { 3 }, q _ { 1 }, q _ { 2 }, \, q _ { 3 } \right ) ^ { \top } \in D \subset \mathbb { R } ^ { 6 }, \ u = \left ( u _ { 1 }, u _ { 2 }, u _ { 3 } \right ) ^ { \top } \in \mathbb { R } ^ { 3 }, \quad \text{(5)} \\ & \text{where} \, D = \left \{ \left ( \omega _ { 1 }, \omega _ { 2 }, \omega _ { 3 }, q _ { 1 }, \, q _ { 2 }, \, q _ { 3 } \right ) ^ { \top } \in \mathbb { R } ^ { 6 } \, \colon q _ { 1 } ^ { 2 } + q _ { 2 } ^ { 2 } + q _ { 3 } ^ { 2 } < 1 \right \} \text{ and} \\ & \text{f} ( x, u ) = \left ( \frac { \left ( \frac { \Lambda _ { 1 } - \Lambda _ { 1 } } { \Lambda _ { 1 } } \omega _ { 2 } \omega _ { 3 } + \frac { \omega _ { 2 } ^ { 2 } ( \Lambda _ { 1 } - \Lambda _ { 3 } } { \Lambda _ { 1 } } \Lambda _ { 2 } ^ { 2 } } \right ) ( q _ { 1 } q _ { 4 } + q _ { 2 } q _ { 3 } ) ( 2 q _ { 1 } ^ { 2 } + 2 q _ { 2 } ^ { 2 } - 1 ) + \frac { 1 } { \Lambda _ { 1 } } ( u _ { 3 } b _ { 3 } - u _ { 3 } b _ { 2 } ) \right ) \\ & \text{f} ( x, u ) = \left ( \frac { \left ( \frac { \Lambda _ { 1 } - \Lambda _ { 1 } } { \Lambda _ { 1 } } \omega _ { 3 } + \frac { \omega _ { 2 } ^ { 2 } ( \Lambda _ { 1 } - \Lambda _ { 1 } } { \Lambda _ { 1 } } \Lambda _ { 2 } ^ { 2 } } \right ) ( q _ { 1 } q _ { 3 } - q _ { 2 } q _ { 4 } ) ( 2 q _ { 1 } ^ { 2 } + 2 q _ { 2 } ^ { 2 } - 1 ) + \frac { 1 } { \Lambda _ { 1 } } ( u _ { 3 } b _ { 1 } - u _ { 1 } b _ { 3 } ) \right ) \\ & \text{f} ( x, u ) = \left ( \frac { \Lambda _ { 1 } - \Lambda _ { 1 } } { \Lambda _ { 1 } } \omega _ { 2 } + \frac { 1 2 \omega _ { 3 } ^ { 2 } ( \Lambda _ { 1 } - \Lambda _ { 1 } } { \Lambda _ { 1 } } ( q _ { 1 } q _ { 4 } + q _ { 2 } q _ { 3 } ) ( 2 q _ { 2 } q _ { 4 } - q _ { 1 } q _ { 3 } ) + \frac { 1 } { \Lambda _ { 1 } } ( u _ { 1 } b _ { 2 } - u _ { 2 } b _ { 1 } ) \right ) \\ & \quad \frac { \frac { \Lambda _ { 1 } } { \Lambda _ { 1 } } ( \left ( \omega _ { 1 } + \omega _ { 0 } ) q _ { 4 } - \omega _ { 2 } q _ { 3 } + \omega _ { 3 } q _ { 2 } ) } { \Lambda _ { 1 } } ( \Lambda _ { 1 } - \omega _ { 0 } ) q _ { 3 } + \omega _ { 2 } q _ { 4 } - \omega _ { 3 } q _ { 1 } ) \\ & \quad \frac { 1 } { \Lambda _ { 1 } } ( - ( \omega _ { 1 } - \omega _ { 0 } ) q _ { 2 } + \omega _ { 2 } q _ { 1 } + \omega _ { 3 } q _ { 4 } )$$

$$q _ { 4 } = \sqrt { 1 - q _ { 1 } ^ { 2 } - q _ { 2 } ^ { 2 } - q _ { 3 } ^ { 2 } }, \\ b _ { 1 } = \beta ( 1 - 2 q _ { 2 } ^ { 2 } - 2 q _ { 3 } ^ { 2 } ), \ b _ { 2 } = 2 \beta ( q _ { 1 } q _ { 2 } - q _ { 3 } q _ { 4 } ), \ b _ { 3 } = 2 \beta ( q _ { 1 } q _ { 3 } + q _ { 2 } q _ { 4 } ).$$

We do not consider the differential equation for ˙ q 4 in (5) because the attitude quaternion is normalized ( q 2 1 + q 2 2 + q 2 3 + q 2 4 = 1 ), so in our case, the component q 4 is defined by (6) for all x ∈ D .

Thus, we have derived the nonlinear control system (5) as a mathematical model of the considered orbiting satellite . This control system can be obtained from the mathematical model of [16] by neglecting the deformation of the flexible part of the satellite.

## 3 Main results

The considered system (5) admits an equilibrium ( x , u e e ) ∈ D × R 3 of the form

$$( x _ { e }, u _ { e } ) \, \colon \, \omega _ { 1 } = - \omega _ { 0 }, \omega _ { 2 } = \omega _ { 3 } = 0, \, q _ { 1 } = q _ { 2 } = q _ { 3 } = 0, \, u _ { 1 } = u _ { 2 } = u _ { 3 } = 0,$$

i.e. f ( x , u e e ) = 0 . This equilibrium corresponds to the circular motion of the satellite with the constant orbital rate ω 0 when the BCS aligns with the OCS.

To study the controllability property of system (5) near this equilibrium, we recall the notion of small-time local controllability. For a point ξ ∈ R n and ε &gt; 0 , we denote the ε -neighborhood of ξ by B ε ( ξ ) = { ˜ ξ ∈ R n : ‖ ˜ ξ - ‖ ξ &lt; ε } . All the norms in this paper are Euclidean.

Definition 1. System (5) is called small-time locally controllable at the equilibrium ( x , u e e ) if, for every ε &gt; 0 , there exists δ &gt; 0 such that, for every x (0) , x (1) ∈ B δ ( x e ) , there exists a function u ∈ L ∞ ([0 , ε ]; R 3 ) such that

$$\| u ( t ) - u _ { e } \| \leqslant \varepsilon \quad \text{for all } t \in [ 0, \varepsilon ],$$

and the solution x t ( ) of the corresponding Cauchy problem

$$\dot { x } ( t ) = f ( x ( t ), u ( t ) ), \quad x ( 0 ) = x ^ { ( 0 ) }$$

is such that

$$x ( \varepsilon ) = x ^ { ( 1 ) }.$$

Our main result reads as follows.

/negationslash

/negationslash

Theorem 1. System (5) is small-time locally controllable at the equilibrium ( x , u e e ) , provided that I 1 = I 2 + I 3 and I 2 = I 3 .

/negationslash

/negationslash

Remark. From the physical viewpoint, the restrictions of Theorem 1 mean that the mass distribution is not degenerate to a flat body case ( I 1 = I 2 + I 3 ), and the body is not symmetric with respect to its first principal axis of inertia ( I 2 = I 3 ).

The proof of Theorem 1 will be given in Section 4 using Lie algebraic sufficient controllability conditions due to H.J. Sussmann [24]. It should be emphasized that the linear approximation of system (5) at ( x , u e e ) is not controllable . To show this, we first identify the components of x as

$$x _ { 1 } = \omega _ { 1 }, \, x _ { 2 } = \omega _ { 2 }, \, x _ { 3 } = \omega _ { 3 }, \, x _ { 4 } = q _ { 1 }, \, x _ { 5 } = q _ { 2 }, \, x _ { 6 } = q _ { 3 },$$

$$& \text{$\mathcal{ }$ $v$ first identify the components of $x$ as} \\ & \text{$x_{1}=\omega_{1},\, x_{2}=\omega_{2},\, x_{3}=\omega_{3},\, x_{4}=q_{1},\, x_{5}=q_{2},\, x_{6}=q_{3},$} \\ & \text{and compute the Jacobian matrices} \frac { \mathcal { }$ of($x,u)\mathfrak { a } } { \partial x } = \left ( \frac { \mathfrak { d } f ( x, u ) } { \partial x _ { x } } \right ), \frac { \mathcal { }$ of($x,u)\mathfrak { a } } { \partial u } = \left ( \frac { \mathfrak { d } f ( x, u ) } { \partial u _ { y } } \right ) \text{at} \left ( x _ { e }, u _ { e } \colon \\ & \text{$A = \frac{ \mathfrak { d } f ( x, u ) } { \partial x } \right | _ { ( x, u ) = ( x _ { e }, u _ { e } ) } = \left ( \frac { 0 } { \mathfrak { d } f ( x, u ) } \left ( \frac { 1 } { \mathfrak { d } f ( x, u ) } \right ) \text{at} \left ( x _ { e }, u _ { e } \colon \\ & \text{$A = \frac{ \mathfrak { d } f ( x, u ) } { \mathfrak { d } f ( x, u ) = ( x _ { e }, u _ { e } ) } \right ) = \left ( \frac { 0 } { \mathfrak { d } f ( x, u ) } \left ( \frac { 1 } { \mathfrak { d } f ( x, u ) } \right ) \text{$\mathcal{ }$} \right ), \\ & \text{$B = \frac{ \mathfrak { d } f ( x, u ) } { \partial u } \right | _ { ( x, u ) = ( x _ { e }, u _ { e } ) } = \left ( \frac { 0 } { \mathfrak { d } f ( x, u ) } \right ) \text{$\mathcal{ }$} \right ). \\ & \text{$B = \frac{ \mathfrak { d } f ( x, u ) } { \mathfrak { d } u } \right ) \text{$\mathcal{ }$} \right ).$$

Then the linearization of (5) at the considered equilibrium reads as

$$\dot { \xi } = A \xi + B v, \ \xi \in \mathbb { R } ^ { 6 }, \ v \in \mathbb { R } ^ { 3 }.$$

The rank of the block matrix [ B,AB,A B,A B,A B,A B 2 3 4 5 ] is at most 4 ; therefore, Kalman's controllability condition for system (7) fails. This means that the linear test cannot be applied to the local controllability analysis of system (5) near the equilibrium ( x , u e e ) .

## 4 Small-time local controllability: Proof of Theorem 1

For further analysis, we rewrite system (5) in the following control-affine form:

$$\dot { x } = f _ { 0 } ( x ) + \sum _ { i = 1 } ^ { 3 } u _ { i } f _ { i } ( x ),$$

where the smooth vector fields f i : D → R 6 are

$$\text{mooth vector fields } f _ { i } \colon & D \to \mathbb { R } ^ { 6 } \text{ are} \\ f _ { 0 } ( x ) = & \left ( \frac { I _ { 2 } - I _ { 3 } } { I _ { 1 } } \omega _ { 2 } \omega _ { 3 } + \frac { 6 \omega _ { 6 } ^ { 2 } ( I _ { 2 } - I _ { 3 } ) } { I _ { 1 } } ( q _ { 1 } q _ { 4 } + q _ { 2 } q _ { 3 } ) ( 2 q _ { 1 } ^ { 2 } + 2 q _ { 2 } ^ { 2 } - 1 ) \\ \frac { I _ { 3 } - I _ { 1 } } { I _ { 2 } } \omega _ { 1 } \omega _ { 3 } + \frac { 6 \omega _ { 6 } ^ { 2 } ( I _ { 3 } - I _ { 1 } ) } { I _ { 2 } } ( q _ { 1 } q _ { 3 } - q _ { 2 } q _ { 4 } ) ( 2 q _ { 1 } ^ { 2 } + 2 q _ { 2 } ^ { 2 } - 1 ) \\ \frac { I _ { 1 } - I _ { 2 } } { I _ { 3 } } \omega _ { 1 } \omega _ { 2 } + \frac { 1 2 \omega _ { 6 } ^ { 2 } ( I _ { 1 } - I _ { 2 } ) } { I _ { 3 } } ( q _ { 1 } q _ { 4 } + q _ { 2 } q _ { 3 } ) ( q _ { 2 } q _ { 4 } - q _ { 1 } q _ { 3 } ) \\ \frac { 1 } { 2 } ( ( \omega _ { 1 } + \omega _ { 0 } ) q _ { 4 } - \omega _ { 2 } q _ { 3 } + \omega _ { 3 } q _ { 2 } ) \\ \frac { 1 } { 2 } ( ( \omega _ { 1 } - \omega _ { 0 } ) q _ { 3 } + \omega _ { 2 } q _ { 4 } - \omega _ { 3 } q _ { 1 } ) \\ \frac { 1 } { 2 } ( - ( \omega _ { 1 } - \omega _ { 0 } ) q _ { 2 } + \omega _ { 2 } q _ { 1 } + \omega _ { 3 } q _ { 4 } ) \\ / \quad \varrho \quad \ \ & \ \int _ { 2 \beta ( \alpha, \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \bar { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha } _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \ \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \alpha _ { \ \alpha _ { \ \alpha _ { \alpha _ { \ \alpha _ { \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \ | \alpha _ { \ \ | \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \ \ | \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \ \alpha _ { \ \ | \alpha _ { \ \ | \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \alpha _ { \ \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \ \ | \alpha _ { \alpha _ { \ \ | \alpha _ { \alpha _ { \ \ | \alpha _ { \ \alpha _ { \ \alpha _ { \ \ | \alpha _ { \ \alpha _ { \ \ | \alpha _ { \ \alpha _ { \ \alpha _ { \ \ | \alpha _ { \ \alpha _ { \ \alpha _ { \ \ | \alpha _ { \ \alpha _ { \ \alpha _ { \ \ | \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \ | \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \alpha _ { \alpha _ { \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \ \alpha _ { \alpha _ { \ \alpha _ { \ \alpha _ \ \alpha _ { \ \alpha _ \ \alpha _ { \ \alpha _ { \ \alpha _ \ \alpha _ \alpha _ \ \alpha _ \ \alpha _ \ \ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } ) } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } _ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } \; } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } }$$

$$f _ { 1 } ( x ) = \begin{pmatrix} 0 \\ - \frac { 2 \beta ( q _ { 1 } q _ { 3 } + q _ { 2 } q _ { 4 } ) } { 2 \beta ( q _ { 1 } q _ { 2 } - q _ { 3 } q _ { 4 } ) } \\ 0 \\ 0 \\ 0 \end{pmatrix}, \, f _ { 2 } ( x ) = \begin{pmatrix} 2 \beta ( q _ { 1 } q _ { 3 } + q _ { 2 } q _ { 4 } ) \\ 0 \\ \beta ( 2 q _ { 2 } + 2 q _ { 3 } - 1 ) \\ 0 \\ 0 \\ 0 \end{pmatrix}, \, f _ { 3 } ( x ) = \begin{pmatrix} 2 \beta ( q _ { 3 } q _ { 4 } - q _ { 1 } q _ { 2 } ) \\ \frac { I _ { 1 } } { \beta ( 1 - 2 q _ { 2 } - 2 q _ { 3 } ) } { I _ { 2 } } \\ 0 \\ 0 \\ 0 \end{pmatrix}, \\ \end{pmatrix}, \, f _ { 2 } ( x ) = \begin{pmatrix} 2 \beta ( q _ { 3 } q _ { 4 } - q _ { 1 } q _ { 2 } ) \\ \frac { I _ { 1 } } { \beta ( 1 - 2 q _ { 2 } - 2 q _ { 3 } ) } { I _ { 2 } } \\ 0 \\ 0 \\ 0 \end{pmatrix},$$

and the component q 4 is defined by (6).

We denote the above family of vector fields by F = { f , f 0 1 , ..., f m } , m = 3 . For f , f i j ∈ F , the Lie bracket of f i and f j is

$$[ f _ { i }, f _ { j } ] ( x ) = \frac { \partial f _ { j } ( x ) } { \partial x } f _ { i } ( x ) - \frac { \partial f _ { i } ( x ) } { \partial x } f _ { j } ( x ),$$

where ∂f i ( x ) ∂x is the Jacobian matrix of f i ( x ) . The Lie algebra of vector fields generated by F will be denoted by L F ( ) (for more clarifications, see [8, Section 3.2]). To prove Theorem 1, we will use Sussmann's controllability condition [24].

Definition 2. The control affine-system (8) satisfies the Lie algebra rank condition at ( x , u e e ) if

$$\{ g ( x _ { e } ) \colon g \in \mathcal { L } ( \mathcal { F } ) \} = \mathbb { R } ^ { 6 }.$$

Now we denote by Br ( F ) the set of iterated Lie brackets of F and, for an h ∈ Br ( F ) , we denote by δ i ( h ) the number of times that f i appears in h .

Definition 3. A control-affine system ˙ x = f 0 ( x ) + ∑ m i =1 u f i i ( x ) satisfies the Sussmann condition S θ ( ) for some θ ∈ [0 , + ∞ ) if, for any h ∈ Br ( F ) with δ 0 ( h ) odd and δ i ( h ) even for all i = 1 , m , the following condition holds:

$$h ( x _ { e } ) \in D _ { \theta } ( h ) \coloneqq \text{span} \{ g ( x _ { e } ) \colon g \in B r ( \mathcal { F } ) \text{ and } \delta ( g ) < \delta ( h ) \}, & & ( 1 0 )$$

where

$$\delta ( g ) \coloneqq \theta \delta _ { 0 } ( g ) + \sum _ { i = 1 } ^ { m } \delta _ { i } ( g ).$$

We refer to the following sufficient conditions for small-time local controllability due to H.J. Sussmann [24].

Theorem2. If system (8) satisfies the Lie algebra rank condition at an equilibrium ( x , u e e ) and the Sussmann condition S θ ( ) with some θ ∈ [0 , 1] , then system (8) is small-time locally controllable.

Proof of Theorem 1. We check the conditions of Theorem 2 by straightforward computation, in a similar way as it was done in [18, Theorem 2]. It is common to call h ∈ Br f ( ) a 'bad' bracket if it satisfies δ 0 ( h ) odd and δ i ( h ) is even for every i = 1 3 , , and a 'good' bracket if this verification fails (cf. [2,8]). The main task is to find six 'good' brackets g i which span R 6 at x = x e , and then analyze all possible cases for h in terms of condition (10).

As a set of 'good' brackets, we take

$$g _ { 1 } & = [ f _ { 0 }, [ f _ { 3 }, [ f _ { 0 }, [ f _ { 0 }, f _ { 2 } ] ] ] ], \, g _ { 2 } = f _ { 3 }, \, g _ { 3 } = f _ { 2 }, \\ g _ { 4 } & = [ f _ { 2 }, [ f _ { 0 }, [ f _ { 0 }, f _ { 3 } ] ] ], \, g _ { 5 } = [ f _ { 0 }, f _ { 3 } ], \, g _ { 6 } = [ f _ { 0 }, f _ { 2 } ].$$

The evaluation of these vector fields at the equilibrium yields

$$\text{vector fields at the equilibrium yields} \\ g _ { 1 } ( x _ { e } ) & = \frac { 3 \omega _ { 0 } ^ { 2 } \beta ^ { 2 } ( I _ { 1 } - I _ { 2 } - I _ { 3 } ) ( I _ { 2 } - I _ { 3 } ) } { I _ { 1 } ^ { 2 } I _ { 2 } I _ { 3 } } \left ( 1, 0, 0, 0, 0, 0 \right ) ^ { \top }, \\ g _ { 2 } ( x _ { e } ) & = \frac { \beta } { I _ { 2 } } \left ( 0, 1, 0, 0, 0, 0 \right ) ^ { \top }, \\ g _ { 3 } ( x _ { e } ) & = - \frac { \beta } { I _ { 3 } } \left ( 0, 0, 1, 0, 0, 0 \right ) ^ { \top }, \\ g _ { 4 } ( x _ { e } ) & = \frac { \beta ^ { 2 } ( I _ { 2 } + I _ { 3 } - I _ { 1 } ) } { 2 I _ { 1 } I _ { 2 } I _ { 3 } } \left ( 0, 0, 0, 1, 0, 0 \right ) ^ { \top }, \\ g _ { 5 } ( x _ { e } ) & = \left ( 0, 0, \frac { \omega _ { 0 } \beta ( I _ { 1 } - I _ { 2 } ) } { I _ { 2 } I _ { 3 } }, 0, - \frac { \beta } { 2 I _ { 2 } }, 0 \right ) ^ { \top }, \\ g _ { 6 } ( x _ { e } ) & = \left ( 0, \frac { \omega _ { 0 } \beta ( I _ { 1 } - I _ { 3 } ) } { I _ { 2 } I _ { 3 } }, 0, 0, 0, \frac { \beta } { 2 I _ { 3 } } \right ) ^ { \top }. \\ \text{Theorem 1, it is easy to see that}$$

Under the conditions of Theorem 1, it is easy to see that

$$\text{span} \{ g _ { 1 } ( x _ { e } ), g _ { 2 } ( x _ { e } ), g _ { 3 } ( x _ { e } ), g _ { 4 } ( x _ { e } ), g _ { 5 } ( x _ { e } ), g _ { 6 } ( x _ { e } ) \} = \mathbb { R } ^ { 6 },$$

which proves that the Lie algebra rank condition (9) holds.

Taking θ = 1 2 , we have

$$\delta ( g _ { i } ) \leqslant \frac { 7 } { 2 } \text{ \ for all } i = \overline { 1, 6 }.$$

Now for 'bad' brackets h ∈ Br f ( ) with δ h ( ) &gt; 7 2 , Sussmann's condition is obviously satisfied. If δ h ( ) /lessorequalslant 7 2 , then we have three possibilities:

- 1) δ 0 ( h ) /lessorequalslant 7 , δ i ( h ) = 0 for all i = 1 3 , . In this case, the corresponding iterated Lie brackets vanish at x e .
- 2) In the case where δ 0 ( h ) = 1 and δ i ( h ) = 2 for one of i ∈ { 1 2 3 , , } , we have

$$[ f _ { i }, [ f _ { i }, f _ { 0 } ] ] \equiv [ f _ { 0 }, [ f _ { i }, f _ { i } ] ] \equiv 0 \text{ \ for all } i = \overline { 1, 3 }.$$

- 3) In the case where δ 0 ( h ) = 3 and δ i ( h ) = 2 for one of i ∈ { 1 2 3 , , } , the only nonzero brackets at x e are:

$$h _ { 1 } = [ f _ { 2 }, [ f _ { 0 }, [ f _ { 0 }, f _ { 2 } ] ] ] ] = - \frac { \omega _ { 0 } \beta ^ { 2 } ( I _ { 1 } ^ { 2 } - ( I _ { 2 } + 2 I _ { 3 } ) I _ { 1 } + 2 I _ { 2 } I _ { 3 } + I _ { 3 } ^ { 2 } ) } { I _ { 1 } I _ { 2 } I _ { 3 } ^ { 2 } } \, ( 0, 0, 0, 1, 0, 0 ) ^ { \top },$$

$$h _ { 2 } = [ f _ { 3 }, [ f _ { 0 }, [ f _ { 0 }, [ f _ { 0 }, f _ { 3 } ] ] ] ] = - \frac { \omega _ { 0 } \beta ^ { 2 } ( I _ { 1 } ^ { 2 } - ( 2 I _ { 2 } + I _ { 3 } ) I _ { 1 } + I _ { 2 } ( I _ { 2 } + 2 I _ { 3 } ) ) } { I _ { 1 } I _ { 2 } ^ { 2 } I _ { 3 } } \left ( 0, 0, 0, 1, 0, 0 \right ) ^ { \top }.$$

As δ g ( 4 ) = 3 &lt; δ h ( i ) for i = 1 2 , , we state that h i ∈ D 1 2 ( h i ) for i = 1 2 , .

Hence, Sussmann's condition holds which concludes the proof of Theorem 1.

## 5 Conclusion

A particular outcome of the proof of Theorem 1 provides the set of 'good' Lie brackets (11) for control system (8). According to the available stabilizability results for essentially nonlinear control systems [11, 30], the structure of iterated brackets in the controllability condition is crucial for defining the stabilizing controllers. Control design issues for system (8) under the bracket-generating condition, with the brackets of the form (11), are left as prospective directions for future research.

## References

- [1] A. A. Agrachev. Is it possible to recognize local controllability in a finite number of differentiations? Open problems in mathematical systems and control theory , 15-18, 1999.
- [2] A. A. Agrachev. 'Good Lie Brackets' for Control Affine Systems. Journal of Dynamical and Control Systems , 30(5), 2024. https://doi.org/10.1007/s10883-023-09674-w.
- [3] A. A. Agrachev, Y. Sachkov. Control theory from the geometric viewpoint. Springer Science &amp; Business Media , 2013. https://doi.org/10.1137/0320017.
- [4] F. Bullo, N. E. Leonard, A. D. Lewis. Controllability and motion algorithms for underactuated Lagrangian systems on Lie groups. IEEE Transactions on Automatic Control , 45(8): 1437-1454, 2000. https://doi.org/10.1109/9.871753.
- [5] C. I. Byrnes. On Brockett's Necessary Condition for Stabilizability and the Topology of Liapunov Functions on R n . Communications in information and systems , 8(4):333-352, 2008.
- [6] Y. Chitour. A continuation method for motion-planning problems. ESAIM: Control, Optimisation and Calculus of Variations , 12(1):139-168, 2006. https://doi.org/10.1051/cocv:2005035.
- [7] W. L. Chow. Über Systeme von linearen partiellen Differentialgleichungen erster Ordnung. Mathematische Annalen , 117:98-105, 1940. https://doi.org/10.1007/bf01450011.
- [8] J.-M. Coron. Control and Nonlinearity . American Mathematical Society, USA, 2007.
- [9] P. Crouch. Spacecraft attitude control and stabilization: Applications of geometric control theory to rigid body models. IEEE Transactions on Automatic Control , 29(4):321-331, 1984.
- [10] L. Giraldi, P. Lissy, J. B. Pomet, C. Moreau. Necessary conditions for local controllability of a particular class of systems with two scalar controls. ESAIM: Control, Optimisation and Calculus of Variations , 30(4):1-37, 2024. https://doi.org/10.1051/cocv/2023073.
- [11] V. Grushkovskaya, A. Zuyev. Design of Stabilizing Feedback Controllers for High-Order Nonholonomic Systems. IEEE Control Systems Letters , 8:988-993, 2024. https://doi.org/10.1109/LCSYS.2024.3406931.
- [12] H. Hermes. On local and global controllability. SIAM Journal on Control , 12(2):252-261, 1974. https://doi.org/10.1137/0312019.
- [13] H. Hermes. Lie Algebras of Vector Fields and Local Approximation of Attainable Sets. SIAM Journal on Control and Optimization , 16(5):715-727, 1978. https://doi.org/10.1137/0316047.
- [14] H. Hermes. On local controllability. SIAM Journal on Control and Optimization , 20(2): 211-220, 1982. https://doi.org/10.1137/0320017.
- [15] Z. Ji, X. Zhang, D. Cheng. Global Controllability Criteria and Motion Planning of Regular Affine Systems With Drifts. IEEE Control Systems Letters , 7:2581-2586, 2023.
- [16] J. Kalosha, Y. Yevgenieva, A. Zuyev. Partial Stabilization of an Orbiting Satellite Model with a Flexible Attachment. Journal of Mathematical Sciences , 211(1): 111-119, 2024.
- [17] M. Kawski. High-order small-time local controllability. H.J. Sussmann (ed.), Nonlinear controllability and optimal control , 431-468, 1990. https://doi.org/10.1201/9780203745625.
- [18] E. Y. Kerai. Analysis of small time local controllability of the rigid body model. IFAC Proceedings Volumes , 28(8):597-602, 1995. IFAC Conference on System Structure and Control 1995, Nantes, France, 5-7 July 1995. https://doi.org/10.1016/S1474-6670(17)45529-7
- [19] F. L. Markley, J. L. Crassidis Fundamentals of spacecraft attitude determination and control , vol. 1286. Springer, 2014.
- [20] R. Misra, R. Wisniewski, A. Zuyev. Attitude Stabilization of a Satellite Having Only Electromagnetic Actuation Using Oscillating Controls. Aerospace , 9(8): 444, 2022. https://doi.org/10.3390/aerospace9080444.
- [21] P. K. Rashevskii. About connecting two points of complete non-holonomic space by admissible curve (in Russian). Uch. Zapiski Ped. Inst. Libknexta , 2:83-94, 1938.
- [22] G. Stefani. On the local controllability of a scalar-input control system. C.I. Byrnes and A. Lindquist (eds.), Theory and Applications of Nonlinear Control Systems , 1986.
- [23] H. J. Sussmann. Lie brackets and local controllability: a sufficient condition for scalar-input systems. SIAM J. Control Optim. , 21:686-713, 1983. https://doi.org/10.1137/0321042.

- [24] H. J. Sussmann. A general theorem on local controllability. SIAM Journal on Control and Optimization , 25(1):158-194, 1987.
- [25] H. J. Sussmann, V. Jurdjevic. Controllability of nonlinear systems. Journal of Differential Equations , 12:95-116, 1972.
- [26] H. J. Sussmann, V. Jurdjevic. Control systems on Lie groups. Journal of Differential equations , 12(2):313-329, 1972.
- [27] J. R. Wertz. Spacecraft attitude determination and control , vol. 73. Springer Science &amp; Business Media, 2012.
- [28] R. Wisniewski, M. Blanke. Three-Axis Satellite Attitude Control Based on Magnetic Torquing-Linear Optimal Approach. IFAC Proceedings Volumes , 29(8): 91-96, 1996. https://doi.org/10.1016/S1474-6670(17)43682-2.
- [29] Y. Yang. Controllability of spacecraft using only magnetic torques. IEEE Transactions on Aerospace and Electronic Systems , 52(2): 954-961, 2016. https://doi.org/10.1109/TAES.2015.150520.
- [30] A. Zuyev, V. Grushkovskaya. On stabilization of nonlinear systems with drift by time-varying feedback laws. In 2019 12th International Workshop on Robot Motion and Control (RoMoCo) , IEEE, 9-14, 2019. https://doi.org/10.1109/RoMoCo.2019.8787353
Yevgeniia Yevgenieva Alexander Zuyev Julia Kalosha
On the Controllability of an Orbiting Satellite Model with Electromagnetic-only Actuation
This paper presents sufficient conditions for small-time local
controllability of a control-affine system that describes the rotational motion
of a satellite in a circular orbit. The satellite is modeled as a rigid body
subject to electromagnetic actuation. We focus on the underactuated scenario
where the control torque is generated solely by magnetorquers. The main
contributions of this work include proving small-time local controllability
around the relative equilibrium under some natural assumptions on the mass
distribution of the rigid body. This result is based on the Lie algebra rank
condition and Sussmann's controllability condition. Furthermore, it is shown
that the linearized system is not controllable in a neighborhood of the
considered equilibrium.