Next Article in Journal
E-Sail Optimal Trajectories to Heliostationary Points
Next Article in Special Issue
Numerical Analysis of Glauert Inflow Formula for Single-Rotor Helicopter in Steady-Level Flight below Stall-Flutter Limit
Previous Article in Journal
The Submerged Nozzle Damping Characteristics in Solid Rocket Motor
Previous Article in Special Issue
Influence of Spinner Shape on Droplet Impact over Rotating Spinners
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Attitude Dynamics of Spinning Magnetic LEO/VLEO Satellites

by
Vladimir S. Aslanov
1,* and
Dmitry A. Sizov
1,2
1
Department of Theoretical Mechanics, Samara National Research University, 34, Moskovskoye Shosse, Samara 443086, Russia
2
Department of Chemical Technology, Samara State Technical University, Syzran Branch, 45, Sovetskaya St., Syzran 446001, Russia
*
Author to whom correspondence should be addressed.
Aerospace 2023, 10(2), 192; https://doi.org/10.3390/aerospace10020192
Submission received: 29 December 2022 / Revised: 14 February 2023 / Accepted: 14 February 2023 / Published: 17 February 2023
(This article belongs to the Special Issue Advances in Aerospace Sciences and Technology III)

Abstract

:
With the growing popularity of small satellites, the interaction with the air in low and especially in very low Earth orbits becomes a significant resource for passive angular stabilisation. However, the possibility of spin motion remains a considerable challenge for missions involving aerodynamically stabilised satellites. The goal of this paper was to investigate the attitude motion of arbitrarily spinning satellites in LEO and VLEO under the action of aerodynamic, gravitational, and magnetic torques, taking into account the aerodynamic damping. Using an umbrella-shaped deployable satellite as an example, the study demonstrated that both regular and chaotic attitude regimes are possible in the attitude motion. The occurrence of chaos was verified by means of Poincaré sections. The results revealed that, to prevent chaotic motion, active attitude control and reliable deployment techniques for aerodynamically stabilised satellites are needed.

1. Introduction

If a satellite is in LEO or very low Earth orbit (VLEO) (typically characterised by altitudes of 80–450 km), its interaction with the atmosphere needs to be taken into account [1,2,3,4]. Aerodynamic torque resulting from this interaction can be used for angular stabilisation. Passive aerodynamic stabilisation is the simplest, since it does not require any power supply. This type of stabilisation has been studied since the late 1950s [5,6,7,8,9]. Currently, with the growing popularity of small satellites, which are more sensitive to the interaction with the air, the number of missions using partial or total aerodynamic stabilisation has increased substantially. Current research in this field is focused on atmosphere-breathing electric propulsion systems [10,11]; satellite shape optimisation and drag reduction [12,13,14]; using aerodynamic forces for formation flying control [15,16]; and creating variable shapes in order to change satellite equilibrium positions [17]. Another popular topic related to aerodynamic stabilisation is the design of re-entry capsules [18,19,20,21,22] and their systems [23].
In the general case, aerodynamically stabilised satellites may spin relative to an orbiting reference frame. A recent in-orbit test [19] of an inflatable satellite showed that this spinning may be caused by distortion of the aeroshell after its deployment. Regardless of the cause of spinning, the joint action of the environmental and gyroscopic torques produces several equilibrium angular positions, which can be either stable or unstable. The classical papers by Meirovitch, Wallace, and Frik [6], regarding spatial attitude motion, dealt only with two environmental torques, namely gravitational and aerodynamic, both depending only on the satellite’s orientation. More recent studies have considered other possible perturbations, such as orbit eccentricity [24], resisting medium [25,26], inertial torques due to elastic elements [27,28,29,30], and magnetic [31,32,33] and damping aerodynamic torques [2,30,33], demonstrating that, in the presence of unstable equilibrium positions, these perturbations may lead to chaos. Most of these papers, aimed at studying chaos in the attitude motion, were focused on the limiting case of planar pitch motion or did not simultaneously take into account magnetic and aerodynamic torques.
The goal of this paper was thus to investigate the attitude motion of arbitrarily spinning satellites in LEO and VLEO under the action of aerodynamic, gravitational, and magnetic torques, taking into account the aerodynamic damping. In order to achieve this goal, the unperturbed and perturbed motions are defined, the equations of motion were derived, the equilibrium positions were found, and a case study dealing with a deployable LEO/VLEO satellite was considered.
The paper has the following structure. In Section 2, the problem of satellite attitude motion under several environmental torques is formulated and the expressions for the torques and their potential functions are considered. Section 3 presents the definitions of the unperturbed and perturbed motions and the derivation of the equations of motion. In addition, the equilibrium positions for the unperturbed motion are found by analysing the dynamic potential. Section 4 is devoted to the case study of the attitude motion of a deployable satellite in LEO and VLEO. Finally, the conclusions are given in Section 5.

2. Problem Statement

The attitude motion of the satellite is considered under the following assumptions:
  • The orbit of the satellite remains circular.
  • The atmosphere is not rotating.
  • The density and temperature of the incident air stream are chosen using the Jacchia–Bowman 2008 Atmosphere Model [34].
  • The satellite is shaped as a body of revolution.
  • The centre of mass of the satellite lies on its axis of symmetry.
  • The transverse moments of inertia of the satellite are equal to each other.
The last two assumptions define the inertia matrix of the satellite:
J = A 0 0 0 C 0 0 0 C
where A is the longitudinal moment of inertia, i.e., about the axis of symmetry of the satellite, and C is the moment of inertia about the axes perpendicular to the axis of symmetry.

2.1. Coordinate Systems

The following coordinate systems or reference frames, denoted using a convenient notation proposed in [35], are used:
  • The orbital frame O is defined through an orthonormal right-hand set of unit vectors o ^ k , k = 1 , 2 , 3 with an origin at the satellite’s centre of mass (Figure 1). The o ^ 1 vector is tangential to the orbit in the flight direction, the o ^ 3 vector is directed along the radius vector from the centre of the Earth to the centre of mass of the satellite.
  • The body-fixed frame B is defined through a set b ^ k . These vectors coincide with the satellite’s principal axes of inertia, and the b ^ 1 vector lies along the axis of symmetry. The orientation of b ^ k relative to o ^ k is described by a symmetric ( 1 , 3 , 1 ) set of Euler angles corresponding to three successive positive rotations: first about Axis 1 by the precession angle ψ , then about Axis 3 by the nutation angle θ (angle of attack), 0 θ π , and finally, about Axis 1 by the spin angle φ , as is shown in Figure 1. The transformation matrix Θ between O and B is defined by
    v B = Θ v O
    where
    Θ = c θ c ψ s θ s θ s θ c φ s θ c θ c φ c ψ s φ s ψ c ψ s φ + c θ c φ s ψ s θ s φ c θ c ψ s φ c φ s ψ c φ c ψ c θ s φ s ψ ,
    c θ = cos θ , s θ = sin θ , etc.
  • The intermediate coordinate frames I and I are defined through unit vector sets i ^ k , i ^ k , respectively. The orientations of these frames relative to the orbital frame are determined by the above-mentioned rotations: a single rotation by the angle ψ for the i ^ k frame and two successive rotations by the angles ψ and θ for the i ^ k frame.
Note that the satellite’s attitude can be described not only by Euler angles, but also quaternions, Cayley–Klein parameters, etc. [35,36]. In this paper, Euler angles were used, since, in the authors’ opinion, they give the clearest physical picture of the considered problem.

2.2. Environmental Torques

The following environmental torques were considered: gravitational, magnetic, and aerodynamic, the last one consisting of restoring and damping torques. For the convenience of the analysis, it is useful to separate the conservative torques from non-conservative ones. Any torque M acting on the body is conservative, if it depends only on the body’s orientation and its curl is equal to zero [8]:
M φ θ M θ φ = 0 , M ψ φ M φ ψ = 0 , M θ ψ M ψ θ = 0
where
M ψ = M · o ^ 1 θ = 0 , φ = 0 , M θ = M · i ^ 3 φ = 0 , M φ = M · i ^ 1 .
If a torque satisfies Condition (4), it has a potential function, which can be found as
V = 0 ψ M ψ d ψ 0 θ M θ d θ 0 φ M φ d φ .
Further in this subsection, we will derive, if possible, the potential functions for each of the environmental torques considered.

2.2.1. Gravitational Torque

The torque due to the gravity gradient is defined as [7]
M G = 3 Ω 2 o ^ 3 × J o ^ 3
where Ω = μ G / R 0 3 is the mean motion, R 0 = R e + h is the orbit radius, R e is the mean radius of the Earth, h is the orbital altitude, and μ G is the gravitational parameter of the Earth. According to Equation (5), the projections of the vector M G on the (1,3,1) axes o ^ 1 , i ^ 3 , i ^ 1 (Figure 1) are, respectively,
M G ψ = 0 , M G θ = 3 2 C A Ω 2 sin 2 ψ sin 2 θ , M G φ = 0 .
The potential energy corresponding to the gravitational torque can be obtained by the substitution of Equation (8) into Equation (6). Using a well-known double-angle formula after the integration and omitting constant terms, we can write the potential energy in the form:
V G = 0 θ M θ d θ = 3 2 C A Ω 2 sin 2 θ sin 2 ψ .
which, for the chosen Euler angles, is symmetric with respect to the angles θ and ψ .

2.2.2. Magnetic Torque

The magnetic torque acting on a satellite is [37]
M M = μ × B
where μ is the vector of the satellite’s own magnetic moment:
μ = M 1 M 2 M 3 B = μ M ¯ 1 M ¯ 2 M ¯ 3 B ,
M ¯ k = M k / μ are dimensionless components of the satellite’s own magnetic moment, B is the Earth’s magnetic field, which is modelled as a field of a dipole placed in the centre of the Earth [7]:
B = μ 0 μ m 4 π R 0 3 3 o ^ 3 k e · o ^ 3 k e
where k e is the dipole direction assumed to be aligned with the Earth’s rotation axis, μ 0 is the magnetic permeability of free space, and μ m is the geomagnetic dipole moment. The magnetic field takes its simplest form if it is expressed in the orbital frame:
B = μ 0 μ m 4 π R 0 3 cos ν sin i cos i 2 sin i sin ν O
In Equation (13), i is the orbital inclination and ν is the true anomaly:
ν = Ω t .
Note, since the magnetic torque is time-dependent for non-equatorial orbits, it is not conservative in general and, thus, has no potential function.

2.2.3. Restoring Aerodynamic Torque

In this paper, the restoring aerodynamic torque is understood as a torque that depends only on the orientation of the body-fixed frame with respect to the flow vector. This orientation is described by two angles, namely the angle of attack θ and the spin angle φ . Thus, the restoring aerodynamic torque is independent of the precession angle ψ . In this case, the conditions (4) for the torque to be conservative are as follows [8]:
M ψ = c o n s t ,
M φ θ = M θ φ .
For all bodies of revolution with the centre of mass lying on the axis of geometric symmetry, the aerodynamic restoring torque about the longitudinal symmetry axis is equal to zero, so Conditions (15) and (16) become
M ψ = M φ = 0 ,
M θ φ = 0 .
Thus, for a body of revolution with the centre of mass on the axis of geometric symmetry, the restoring aerodynamic torque is fully described by the function M θ = f ( θ ) . This function can be written in the following form:
M θ = c r M ¯ θ
where
c r = ρ V 2 2 A r L r ,
V = Ω R 0 is the orbital speed, A r is the reference area, L r is the reference length, and ρ is the density of the atmosphere. The restoring aerodynamic torque coefficient M ¯ θ is a periodic odd function, so it can be approximated by a finite number of terms of a Fourier sine series:
M ¯ θ = i = 1 n b i sin i θ
where n is the number of harmonics. Substituting Equations (17), (19), and (21) into Equation (6), we obtain the potential energy corresponding to the restoring aerodynamic torque:
V A = 0 θ M θ d θ = c r i = 1 n 1 i b i cos i θ .

2.2.4. Damping Aerodynamic Torque

The second component of the aerodynamic torque, which depends not only on the orientation of the body, but also on its angular velocity, is often called the damping torque. In the body-fixed frame B, it can be expressed as [38]
M d = c d M ¯ 1 ω 1 ω 1 M ¯ 2 ω 2 ω 2 M ¯ 3 ω 3 ω 3 B
where
c d = ρ V 2 A r L r 2 ,
where M ¯ k ω k = M ¯ k ω k ( θ , ψ ) , k = 1 , 2 , 3 are the damping torque coefficients associated with the projections ω k of the angular velocity vector ω ,
ω = ω 1 ω 2 ω 3 = φ ˙ + ψ ˙ c θ + Ω c ψ s θ θ ˙ s φ ψ ˙ c φ s θ + Ω c θ c φ c ψ s φ s ψ θ ˙ c φ + ψ ˙ s θ s φ Ω c θ c ψ s φ + s ψ c φ B
on the body axes b ^ k [35]. The “longitudinal” damping torque coefficient M ¯ 1 ω 1 of a body of revolution is an even periodic function of the angle of attack alone, so it can be represented as a Fourier cosine series of a single variable. The “transverse” damping torque coefficients M ¯ 2 ω 2 , M ¯ 3 ω 3 are even periodic functions of both the angle of attack θ and the spin angle φ , so they need to be represented as 2D Fourier cosine series. Moreover, due to the axial symmetry of the body, M ¯ 3 ω 3 is equal to M ¯ 2 ω 2 , calculated for φ increased or decreased by π / 2 . Bearing all these facts in mind, one can write the damping torque coefficients as
M ¯ 1 ω 1 ( θ ) = i = 0 m a 1 , i cos i θ , M ¯ 2 ω 2 ( θ , φ ) = i = 0 p j = 0 r a 2 , i , j cos i θ cos j φ , M ¯ 3 ω 3 ( θ , φ ) = i = 0 p j = 0 r a 2 , i , j cos i θ cos j ( φ + π / 2 ) = M ¯ 2 ω 2 ( θ , φ + π / 2 )
where m, p, and r are the numbers of harmonics. Note that, since, by virtue of Equation (23), the damping aerodynamic torque depends on the angular velocity of the satellite, it has no potential function.

3. Equations of Motion

The attitude motion of aerodynamically stabilised satellites may be treated using the general principles of the analytical mechanics of a rigid body rotating about its centre of mass [36,39,40]. The equations of motion of the satellite are first derived in standard dimensional form using the Euler angles ψ , θ , and φ as generalised coordinates:
q = θ ψ φ T
considering q k as functions of time. Then, the equations are transformed into dimensionless form, using the true anomaly ν as the independent variable. This allows simplifying the expressions and investigating the behaviour of the system in an arbitrary space of parameters.
If some torques acting on the satellite are conservative, it is extremely useful for an intuitive analysis of possible satellite motions to define the perturbed and the unperturbed motions. The perturbed motion is defined as the attitude motion of the satellite under the action of all the environmental torques considered in the previous section. Then, the unperturbed motion is the motion of the satellite under the aerodynamic restoring torque and gravitational torque. Both of them are conservative, which justifies the choice of the Lagrangian approach for the derivation of the equations.

3.1. Unperturbed Motion

The kinetic energy of the satellite is
T = 1 2 A ω 1 2 + C ( ω 2 2 + ω 3 2 ) .
After the substitution of the components of the vector ω , defined by Equation (25), Equation (28) takes the form:
T = 1 2 A φ ˙ + ψ ˙ c θ + Ω c ψ s θ 2 + 1 2 C θ ˙ 2 + Ω c θ c ψ ψ ˙ s θ 2 2 Ω θ ˙ s ψ + Ω 2 s 2 ψ .
The potential energy is the sum of the gravitational and aerodynamic potential energies:
V = V G + V A
where V G and V A are defined by Equations (9) and (22), respectively. The Lagrangian of the system is
L = T V = 1 2 A φ ˙ + ψ ˙ c θ + Ω c ψ s θ 2 + 1 2 C θ ˙ 2 + Ω c θ c ψ ψ ˙ s θ 2 2 Ω θ ˙ s ψ + Ω 2 s 2 ψ V ( θ , ψ ) .
Since the generalised coordinate φ does not appear in the Lagrangian, the unperturbed motion has an integral in the form:
p φ = c o n s t
where
p φ = L φ ˙ = A φ ˙ + ψ ˙ c θ + Ω c ψ s θ .
Since the Lagrangian (31) is not an explicit function of time, the Hamiltonian of the system is also an integral of motion. It can be found as
H = k = 1 2 R q k q ˙ k R = 1 2 C θ ˙ 2 + ψ ˙ 2 s 2 θ 1 2 C Ω 2 c 2 θ c 2 ψ + s 2 ψ + p φ 2 2 A p φ Ω c ψ s θ + V ( θ , ψ ) = c o n s t
where R is the Routhian of the system:
R = L p φ φ ˙ .
The general form of the equations of the unperturbed motion is
d d t R q ˙ k R q k = 0 , k = 1 , 2 , p ˙ φ = 0 .
According to Equation (27), the first two of Equation (36) correspond to the non-cyclic coordinates ψ and θ , respectively, while the third of Equation (36) corresponds to the cyclic coordinate φ . The equations of unperturbed motion in dimensionless form can be obtained by solving Equation (36) for q ¨ k and dividing both sides of each equation by A Ω 2 :
γ ψ = 4 γ 3 s ψ c ψ + β θ s ψ csc θ 2 γ θ c ψ + ψ cot θ , γ θ = α i = 1 n b i sin i θ + 3 γ 1 s θ c θ s 2 ψ + β β γ + γ φ c ψ c θ ψ s θ + γ ψ c ψ , γ φ = γ 1 c θ θ c ψ 2 s 2 ψ + γ ψ s θ θ + s ψ + cot θ φ s ψ θ + ψ c θ 2 γ 1 θ + s ψ
where means differentiation with respect to the true anomaly, and α , β , and γ :
α = A r L ρ R 0 2 2 A ,
β = p φ A Ω = φ + ψ cos θ + sin θ cos ψ ,
γ = C / A
are dimensionless quantities.
In order to analyse possible satellite motions, it is convenient to introduce the dynamic potential. Let us represent the Hamiltonian (34) in the following form, ignoring the constant term p φ 2 / 2 A :
H = T * + U = c o n s t .
In Equation (41), T * = 1 2 C θ ˙ 2 + ψ ˙ 2 sin 2 θ and U is the dynamic potential:
U = V ( θ , ψ ) p φ Ω c ψ s θ 1 2 C Ω 2 c 2 θ c 2 ψ + s 2 ψ .
For the sake of simplicity, one can rewrite the dynamic potential in the dimensionless form:
U ¯ = U A Ω 2 = 3 2 ( 1 γ ) s 2 θ s 2 ψ + α i = 1 n 1 i b i cos i θ β s θ c ψ 1 2 γ c 2 θ c 2 ψ + s 2 ψ .
The dynamic equilibrium positions E of the satellite in its unperturbed motion characterised by q ˙ k = 0 , k = 1 , 2 are the roots of the equations:
q k U ¯ a t E = 0 , k = 1 , 2 .
Equation (44) describes the following three groups of equilibrium positions.
Group E 1 : ψ = 0 and θ is the root of the equation:
α M ¯ θ β cos θ + γ sin θ cos θ = 0 .
Group E 2 : ψ = ± π and θ is the root of the equation:
α M ¯ θ + β cos θ + γ sin θ cos θ = 0 .
Group E 3 : ψ and θ are the roots of the system:
β + ( 3 4 γ ) sin θ cos ψ = 0 , α M ¯ θ + 3 ( γ 1 ) sin θ cos θ = 0 .
The equilibrium is stable if U ¯ has a minimum at E and unstable if it has a maximum.

3.2. Perturbed Motion

The equations of motion for the perturbed system are
d d t L q ˙ k L q k = Q k , k = 1 , 2 , 3 .
The generalised forces due to aerodynamic damping and magnetic torques can be found as
Q k = q ˙ k M M + M d · ω
where M M , M d , and ω are defined by Equations (10), (23), and (25), respectively. The generalised forces (49) are given in their full form in Appendix A. The equations of perturbed motion in dimensionless form are
γ ψ = 4 γ 3 s ψ c ψ + β θ s ψ csc θ 2 γ θ c ψ + ψ cot θ + δ f 2 ψ , ψ , θ , θ , φ + ε g 2 ψ , θ , φ , ν , γ θ = α i = 1 n b i sin i θ + 3 γ 1 s θ c θ s 2 ψ + β β γ + γ φ c ψ c θ ψ s θ + γ ψ c ψ + δ f 1 ψ , ψ , θ , θ , φ + ε g 1 ψ , θ , φ , ν , γ φ = γ 1 c θ θ c ψ 2 s 2 ψ + γ ψ s θ θ + s ψ + cot θ φ s ψ θ + ψ c θ 2 γ 1 θ + s ψ + δ f 3 ψ , ψ , θ , θ , φ , φ + ε g 3 ψ , θ , φ , ν
where δ and ε are dimensionless parameters:
δ = A r L 2 ρ R 0 2 A , 0 < δ 1 ,
ε = μ 0 μ m μ 4 π A μ G .
In Equation (50), the terms δ f k and ε g k represent perturbations due to the aerodynamic damping torque and magnetic torque, respectively. If δ = ε 0 , Equation (50) transforms into the equations of unperturbed motion (37). It can be seen from Equation (52) that the dimensionless parameter ε depends on the satellite’s longitudinal moment of inertia and, by virtue of Equation (11), on the satellite’s own magnetic moment. In this paper, we assumed that the principal component of the magnetic dipole can be estimated in flight and cancelled using magnetopairs. The remaining magnetic torques, expressed by μ and M ¯ k , can thus be of the order of the other perturbations.

4. Case Study: Deployable Satellite

Consider a satellite similar to the one proposed in [21]. It consists of a cylindrical body, a hemispherical nose, and a deployable conic aerobrake. Further in this paper, we will refer to this satellite as the example satellite. Some of its parameters are summarised in Table 1. This paper deals only with the satellite’s final aerodynamic shape, i.e., after the deployment of the aerobrake (Figure 2).
Suppose that, during uneven deployment of the aerobrake, the satellite, for a certain amount of time, takes a skew shape for which M φ 0 . In this case, the interaction with the atmosphere will cause initial spin motion, which will remain after the deployment is completed, resulting in M φ = 0 . This scenario is quite possible, as confirmed by a similar situation with the re-Entry satellite with Gossamer aeroshell and GPS/Iridium (EGG) [19]. In order to numerically study the case, one needs first to determine the aerodynamic characteristics of the satellite, then find its equilibrium positions defined by Equations (45)–(47). The following subsections are devoted to these issues.

4.1. Aerodynamic Characteristics

The aerodynamic torque coefficients for the example satellite are calculated as follows [30]. Taking into account that, for a body with a reference length of about 1 m, the Knudsen number at altitudes above 130 km is larger than 10, which means that the flow is free molecular, then one can assume that the reflected air molecules’ speed distribution is Maxwellian and calculate the pressure and shear stress coefficients using the Schaaf and Chambre’s approach [41]. The numerical results were approximated by Fourier series according to Equations (21) and (26).
Figure 3 shows the calculated restoring aerodynamic torque coefficient and its Fourier series approximation (Equation (21)). It is worth noting that, in the vicinity of θ = π / 2 , there exists a region where M ¯ θ is positive. In some cases, this could provoke instability. One of these cases will be studied below. Figure 4 and Figure 5 depict the coefficients M ¯ 1 ω 1 and M ¯ 2 ω 2 of the damping torque, respectively, as well as their 2D Fourier series approximation (Equation (26)). The plot for coefficient M ¯ 3 ω 3 is not shown, since it can easily be obtained by shifting the M ¯ 2 ω 2 surface by π / 2 along the φ axis by virtue of the last of Equation (26). It can be seen in Figure 5 that, unlike the spin damping coefficient M ¯ 1 ω 1 , the other two damping torque coefficients slightly vary with φ . It is also interesting that the magnitudes of all damping torque coefficients for the considered shape are substantially higher when the angle of attack is small than when the symmetry axis of the satellite is nearly perpendicular to the flow. One may thus expect more intensive attitude motion damping when the satellite is close to its operational orientation.

4.2. Dynamic Potential and Equilibrium Positions in Unperturbed Motion

The dimensionless dynamic potential U ¯ of the unperturbed motion (Equation (43)) is determined by some parameters of the satellite, which were chosen above, namely its shape, the moments of inertia, etc., as well as the initial value of the dimensionless conjugated momentum β (Equation (39)) and the orbital altitude. Let us choose two orbital altitudes: 700 km (LEO) and 250 km (VLEO). At 700 km, the gravitational torque and aerodynamic restoring torque have roughly the same order of magnitude. In VLEO, the satellite’s interaction with the atmosphere is much stronger, so the aerodynamic restoring torque is much larger than the gravitational one. The initial conjugated momentum β depends on the initial conditions and was chosen arbitrarily, in order to represent the shapes of the dynamic potential, which are typical for the problem under consideration.
Figure 6 and Figure 7 represent the dimensionless dynamic potential for the chosen altitudes as 3D surfaces and the corresponding contour plots. The latter can also be understood as phase portraits of the unperturbed system. The stable and unstable equilibrium positions of the satellite are marked with black and white dots, respectively. Figure 6 shows that, for h = 700 km and β = 33 , 000 , the equilibrium positions from all three groups (Equations (45)–(47)) are present. The stable positions belong to Groups E 1 and E 3 . In the case of h = 250 km, β = 20 , illustrated by Figure 7, the situation is substantially different: only the positions from Groups E 1 and E 2 are present, and the stable ones belong to Group E 1 .

4.3. Simulations of Perturbed Motion

In order to illustrate different possible cases of the attitude motion of the example satellite, the Poincaré surfaces of the section were calculated via the numerical integration of the perturbed equations (Equation (50)) for several phase-space trajectories at two orbital altitudes chosen above. The components of the satellite’s own magnetic moment were taken equal to each other and to unity:
M ¯ k = 1 ,
k = 1 , 2 , 3 . The initial conditions were chosen in such a way that all trajectories at the same altitude had the same energy (Equation (41)). The simulations were performed for different absolute values of magnetic perturbation, which were determined by the dimensionless parameter ε (Equation (52)). For a better understanding of the influence of each torque acting on the example satellite, the orders of magnitude of the torques are given in Table 2.
Figure 8 and Figure 9 represent the Poincaré surfaces of the section projected onto the ( θ , ψ ) plane, as well as several typical phase trajectories coloured in red. If its own magnetic moment is absent ( ε = 0 ) , the Poincaré sections are consistent with the phase portrait of the unperturbed system, as shown in Figure 8a and Figure 9a, and all the section points corresponding to the chosen phase trajectories lie in the potential wells, i.e., the satellite oscillates in the vicinity of its stable equilibrium positions. Figure 8b and Figure 9b depict the Poincaré sections for rather weak magnetic perturbations. It can be seen that the system starts to demonstrate chaotic behaviour: the layers connecting different potential wells appear, and the phase trajectories that start in one well may go into another, though the satellite still oscillates near the same set of equilibrium positions, for which π / 2 ψ π / 2 . In the case of stronger magnetic perturbations, the chaotic behaviour becomes more prominent, and full rotations about the precession axis are possible, as is clearly indicated by Figure 8c and Figure 9c, showing the Poincaré sections with the precession angle ψ distributed in the interval from π to π .
Thus, the numerical simulations revealed the possibility of chaos in the attitude motion of aerodynamically stabilised spinning deployable satellites in the presence of magnetic perturbations. To cause the chaotic motion in VLEO, where the interaction with the air is stronger and the aerodynamic damping is significant, the magnitude of perturbations needs to be several orders higher than in LEO.

5. Conclusions

The attitude motion of arbitrarily spinning satellites in LEO and VLEO under the action of aerodynamic, gravitational, and magnetic torques was investigated, taking into account the aerodynamic damping. The separation of the equations of motion into the unperturbed and perturbation parts, proposed in the paper, facilitated an intuitive understanding of the satellite’s behaviour and allowed determining several regimes of angular oscillations corresponding to the wells of the dynamic potential surface. The numerical simulations showed that the magnetic perturbations may cause the chaotic attitude motion of the satellite, manifested in the transition between the potential wells. These results revealed the need for active attitude control to prevent chaotic motion in the case of unexpected spin motion, though this is challenging for small satellites due to their power limits. In this regard, another way of eliminating undesirable spin can be proposed, dealing with the design of reliable satellite deployment techniques, reducing the risk of resulting skew aerodynamic shapes.
The future work will focus on related problems not covered in this paper. In particular, it is interesting to control chaos using magnetorquers or movable aerodynamic surfaces in order to transform irregular attitude motions into periodic ones. Another area of research is to investigate the observed effects in a more general case, taking into account the perturbations caused by the evolution of orbital altitude or by the oscillations of flexible elements.

Author Contributions

Conceptualisation: V.S.A.; investigation: V.S.A. and D.A.S.; methodology: V.S.A. and D.A.S.; software: D.A.S.; writing—original draft: D.A.S.; writing—review and editing: V.S.A. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Russian Science Foundation (Project No. 19-19-00085).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Generalised Forces

The generalised forces (49) in the equations of the perturbed motion (50) can be written as
Q k = A Ω 2 Q ¯ k
where Q ¯ k are the dimensionless generalised forces:
Q ¯ 1 = δ M ¯ 2 ω 2 s 2 φ + M ¯ 3 ω 3 c 2 φ θ s ψ + M ¯ 3 ω 3 M ¯ 2 ω 2 s θ ψ c θ c ψ s φ c φ ε M ¯ 1 c θ + s θ M ¯ 3 s φ M ¯ 2 c φ c ψ c i + M ¯ 1 s θ + M ¯ 2 c φ M ¯ 3 s φ c θ c ν s i + 2 ε M ¯ 1 c θ + M ¯ 3 s φ M ¯ 2 c φ s θ s ν s ψ s i , Q ¯ 2 = δ M ¯ 2 ω 2 M ¯ 3 ω 3 s ψ θ s φ c φ + M ¯ 3 ω 3 ψ s 2 φ s θ + M ¯ 2 ω 2 ψ s θ c 2 φ s θ + δ M ¯ 1 ω 1 φ + M ¯ 1 ω 1 ψ c θ + M ¯ 1 ω 1 M ¯ 2 ω 2 c 2 φ M ¯ 3 ω 3 s 2 φ s θ c ψ c θ + ε M ¯ 1 s θ + M ¯ 2 c φ M ¯ 3 s φ c θ s ψ + M ¯ 2 s φ + M ¯ 3 c φ c ψ c i + 2 ε M ¯ 1 s θ c ψ + M ¯ 2 c φ M ¯ 3 s φ c θ c ψ M ¯ 2 s φ + M ¯ 3 c φ s ψ s ν s i , Q ¯ 3 = δ M ¯ 1 ω 1 β 2 ε M ¯ 2 s φ + M ¯ 3 c φ c θ s ψ + M ¯ 3 s φ M ¯ 2 c φ c ψ s ν s i + ε M ¯ 2 s φ + M ¯ 3 c φ c θ c ψ + M ¯ 2 c φ M ¯ 3 s φ s ψ c i M ¯ 2 s φ + M ¯ 3 c φ c ν s θ s i .
In Equation (A2), dimensionless quantities β , δ , and ε are defined by Equations (39), (51), and (52), respectively.

References

  1. Nurre, G.S. Effects of aerodynamic torque on an asymmetric, gravity-stabilized satellite. J. Spacecr. Rocket. 1968, 5, 1046–1050. [Google Scholar] [CrossRef]
  2. Beletskii, V.; Grushevskii, A. The evolution of the rotational motion of a satellite under the action of a dissipative aerodynamic moment. J. Appl. Math. Mech. 1994, 58, 11–19. [Google Scholar] [CrossRef]
  3. Kuznetsova, E.Y.; Sazonov, V.; Chebukov, S.Y. Evolution of the satellite rapid rotation under the action of the gravitational and aerodynamic torques. Mech. Solids 2000, 35, 1–10. [Google Scholar]
  4. Maslova, A.; Pirozhenko, A. Modeling of the aerodynamic moment acting upon a satellite. Cosm. Res. 2010, 48, 362–370. [Google Scholar] [CrossRef]
  5. DeBra, D.B. The effect of aerodynamic forces on satellite attitude. J. Astronaut. Sci 1959, 6, 40–45. [Google Scholar]
  6. Meirovitch, L.; Wallace, F., Jr. On the effect of aerodynamic and gravitational torques on the attitude stability of satellites. AIAA J. 1966, 4, 2196–2202. [Google Scholar] [CrossRef]
  7. Beletskii, V.V. Motion of an Artificial Satellite about its Center of Mass; NASA TT F-429; NASA: Washington, DC, USA, 1966.
  8. Frik, M.A. Attitude stability of satellites subjected to gravity gradient and aerodynamic torques. AIAA J. 1970, 8, 1780–1785. [Google Scholar] [CrossRef]
  9. Pande, K.; Venkatachalam, R. On optimal aerodynamic attitude control of spacecraft. Acta Astronaut. 1979, 6, 1351–1359. [Google Scholar] [CrossRef]
  10. Romano, F.; Espinosa-Orozco, J.; Pfeiffer, M.; Herdrich, G.; Crisp, N.; Roberts, P.; Holmes, B.; Edmondson, S.; Haigh, S.; Livadiotti, S.; et al. Intake design for an Atmosphere-Breathing Electric Propulsion System (ABEP). Acta Astronaut. 2021, 187, 225–235. [Google Scholar] [CrossRef]
  11. Vaidya, S.; Traub, C.; Romano, F.; Herdrich, G.; Chan, Y.A.; Fasoulas, S.; Roberts, P.; Crisp, N.; Edmondson, S.; Haigh, S.; et al. Development and analysis of novel mission scenarios based on Atmosphere-Breathing Electric Propulsion (ABEP). CEAS Space J. 2022, 14, 689–706. [Google Scholar] [CrossRef]
  12. Aslanov, V.; Sizov, D. 3U CubeSat aerodynamic design aimed to increase attitude stability and orbital lifetime. In Proceedings of the International Astronautical Congress, IAC, Naples, Italy, 12–14 October 2020. [Google Scholar]
  13. Otsu, H. Aerodynamic Characteristics of Re-Entry Capsules with Hyperbolic Contours. Aerospace 2021, 8, 287. [Google Scholar] [CrossRef]
  14. Jiang, Y.; Zhang, J.; Tian, P.; Liang, T.; Li, Z.; Wen, D. Aerodynamic drag analysis and reduction strategy for satellites in Very Low Earth Orbit. Aerosp. Sci. Technol. 2022, 132, 108077. [Google Scholar] [CrossRef]
  15. Ivanov, D.; Monakhova, U.; Guerman, A.; Ovchinnikov, M. Decentralized Control of Nanosatellite Tetrahedral Formation Flying Using Aerodynamic Forces. Aerospace 2021, 8, 199. [Google Scholar] [CrossRef]
  16. Traub, C.; Fasoulas, S.; Herdrich, G.H. A planning tool for optimal three-dimensional formation flight maneuvers of satellites in VLEO using aerodynamic lift and drag via yaw angle deviations. Acta Astronaut. 2022, 198, 135–151. [Google Scholar] [CrossRef]
  17. Miyamoto, K.; Chujo, T.; Watanabe, K.; Matunaga, S. Attitude dynamics of satellites with variable shape mechanisms using atmospheric drag torque and gravity gradient torque. Acta Astronaut. 2023, 202, 625–636. [Google Scholar] [CrossRef]
  18. Yamada, K.; Nagata, Y.; Abe, T.; Suzuki, K.; Imamura, O.; Akita, D. Suborbital reentry demonstration of inflatable flare-type thin-membrane aeroshell using a sounding rocket. J. Spacecr. Rocket. 2015, 52, 275–284. [Google Scholar] [CrossRef]
  19. Berthet, M.; Yamada, K.; Nagata, Y.; Suzuki, K. Feasibility assessment of passive stabilisation for a nanosatellite with aeroshell deployed by orbit-attitude-aerodynamics simulation platform. Acta Astronaut. 2020, 173, 266–278. [Google Scholar] [CrossRef]
  20. Alpert, H.; Miller, R.; Kazemba, C.; Swanson, G.; Cheatwood, F. LOFTID Heat Flux Gauge Calibration: What is Truth? In Proceedings of the International Planetary Probe Workshop, Hamelin, Germany, 19–23 April 2021. [Google Scholar]
  21. Zuppardi, G.; Savino, R.; Mongelluzzo, G. Aero-thermo-dynamic analysis of a low ballistic coefficient deployable capsule in Earth re-entry. Acta Astronaut. 2016, 127, 593–602. [Google Scholar] [CrossRef] [Green Version]
  22. Carná, S.R.; Bevilacqua, R. High fidelity model for the atmospheric re-entry of CubeSats equipped with the drag de-orbit device. Acta Astronaut. 2019, 156, 134–156. [Google Scholar] [CrossRef]
  23. Balaraman, P.; Stephen Joseph Raj, V.; Sreehari, V.M. Static and Dynamic Analysis of Re-Entry Vehicle Nose Structures Made of Different Functionally Graded Materials. Aerospace 2022, 9, 812. [Google Scholar] [CrossRef]
  24. Iñarrea, M.; Lanchares, V. Chaotic pitch motion of an asymmetric non-rigid spacecraft with viscous drag in circular orbit. Int. J. Non-Linear Mech. 2006, 41, 86–100. [Google Scholar] [CrossRef]
  25. Akulenko, L.; Leshchenko, D.; Rachinskaya, A. Evolution of the satellite fast rotation due to the gravitational torque in a dragging medium. Mech. Solids 2008, 43, 173–184. [Google Scholar] [CrossRef]
  26. Leshchenko, D.; Ershkov, S.; Kozachenko, T. Evolution of a heavy rigid body rotation under the action of unsteady restoring and perturbation torques. Nonlinear Dyn. 2021, 103, 1517–1528. [Google Scholar] [CrossRef]
  27. Yue, B.Z. Study on the chaotic dynamics in attitude maneuver of liquid-filled flexible spacecraft. AIAA J. 2011, 49, 2090–2099. [Google Scholar] [CrossRef]
  28. Chegini, M.; Sadati, H.; Salarieh, H. Chaos analysis in attitude dynamics of a flexible satellite. Nonlinear Dyn. 2018, 93, 1421–1438. [Google Scholar] [CrossRef]
  29. Aslanov, V.S. Chaotic attitude dynamics of a LEO satellite with flexible panels. Acta Astronaut. 2021, 180, 538–544. [Google Scholar] [CrossRef]
  30. Aslanov, V.S.; Sizov, D.A. Chaos in flexible CubeSat attitude motion due to aerodynamic instability. Acta Astronaut. 2021, 189, 310–320. [Google Scholar] [CrossRef]
  31. Iñarrea, M. Chaos and its control in the pitch motion of an asymmetric magnetic spacecraft in polar elliptic orbit. Chaos Solitons Fractals 2009, 40, 1637–1652. [Google Scholar] [CrossRef]
  32. Liu, Y.; Chen, L. Chaos in Spatial Attitude Motion of Spacecraft. In Chaos in Attitude Dynamics of Spacecraft; Springer: Berlin/Heidelberg, Germany, 2013; pp. 99–129. [Google Scholar]
  33. Aslanov, V.S.; Sizov, D.A. Chaotic pitch motion of an aerodynamically stabilized magnetic satellite in polar orbits. Chaos Solitons Fractals 2022, 164, 112718. [Google Scholar] [CrossRef]
  34. ISO 14222:2013; International Standard. Space Environment (Natural and Artificial)—Earth Upper Atmosphere. International Organization for Standardization: Geneva, Switzerland, 2013.
  35. Schaub, H.; Junkins, J.L. Analytical Mechanics of Space Systems; AIAA: Reston, VA, USA, 2009. [Google Scholar]
  36. Aslanov, V.S. Rigid Body Dynamics for Space Applications; Butterworth-Heinemann: Oxford, UK, 2017. [Google Scholar]
  37. Shrivastava, S.; Modi, V. Satellite attitude dynamics and control in the presence of environmental torques-A brief survey. J. Guid. Control. Dyn. 1983, 6, 461–471. [Google Scholar] [CrossRef]
  38. Gallais, P. Atmospheric Re-Entry Vehicle Mechanics; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2007. [Google Scholar]
  39. Chernousko, F.L.; Akulenko, L.D.; Leshchenko, D.D. Evolution of Motions of a Rigid Body about Its Center of Mass; Springer: Berlin/Heidelberg, Germany, 2017. [Google Scholar]
  40. Leshchenko, D.; Ershkov, S.; Kozachenko, T. Rotations of a Rigid Body Close to the Lagrange Case under the Action of Nonstationary Perturbation Torque. J. Appl. Comput. Mech. 2022, 8, 1023–1031. [Google Scholar]
  41. Schaaf, S.; Chambre, P. Flow of Rarefied Gases, High Speed Aerodynamics and Jet Propulsion. In Fundamentals of Gas Dynamics; Princeton University Press: Princeton, NJ, USA, 1958; Volume 3. [Google Scholar]
Figure 1. Coordinate frames and Euler angles.
Figure 1. Coordinate frames and Euler angles.
Aerospace 10 00192 g001
Figure 2. Example satellite.
Figure 2. Example satellite.
Aerospace 10 00192 g002
Figure 3. Restoring aerodynamic torque coefficient.
Figure 3. Restoring aerodynamic torque coefficient.
Aerospace 10 00192 g003
Figure 4. Coefficient M ¯ 1 ω 1 of the damping torque.
Figure 4. Coefficient M ¯ 1 ω 1 of the damping torque.
Aerospace 10 00192 g004
Figure 5. Coefficient M ¯ 2 ω 2 of the damping torque.
Figure 5. Coefficient M ¯ 2 ω 2 of the damping torque.
Aerospace 10 00192 g005
Figure 6. Dimensionless dynamic potential in LEO ( h = 700 km, β = 33 , 000 ). Black dots: stable equilibrium; white dots: unstable equilibrium.
Figure 6. Dimensionless dynamic potential in LEO ( h = 700 km, β = 33 , 000 ). Black dots: stable equilibrium; white dots: unstable equilibrium.
Aerospace 10 00192 g006
Figure 7. Dimensionless dynamic potential in VLEO ( h = 250 km, β = 20 ). Black dots: stable equilibrium; white dots: unstable equilibrium.
Figure 7. Dimensionless dynamic potential in VLEO ( h = 250 km, β = 20 ). Black dots: stable equilibrium; white dots: unstable equilibrium.
Aerospace 10 00192 g007
Figure 8. Poincaré sections ( h = 250 km, β = 33,000 ); (a): ε = 0 , (b): ε = 580 , (c): ε = 3900 . Red lines: typical phase trajectories.
Figure 8. Poincaré sections ( h = 250 km, β = 33,000 ); (a): ε = 0 , (b): ε = 580 , (c): ε = 3900 . Red lines: typical phase trajectories.
Aerospace 10 00192 g008
Figure 9. Poincaré sections ( h = 700 km, β = 20 ); (a): ε = 0 , (b): ε = 0.58 , (c): ε = 1.35 . Red lines: typical phase trajectories.
Figure 9. Poincaré sections ( h = 700 km, β = 20 ); (a): ε = 0 , (b): ε = 0.58 , (c): ε = 1.35 . Red lines: typical phase trajectories.
Aerospace 10 00192 g009
Table 1. Parameters of the example satellite.
Table 1. Parameters of the example satellite.
ParameterValue
Body length = Reference length L r 0.3 m
Body radius r0.15 m
Nose radius0.2 m
Reference area A r = π r 2 0.0707  m 2
Aerobrake half-cone angle45°
Aerobrake diameter1.1 m
Longitudinal moment of inertia A                   0.005 kg · m 2
Transverse moment of inertia C                   0.05 kg · m 2
Longitudinal shift of the centre of mass from the body’s geometric centre 0.125  m
Table 2. Orders of magnitude of the torques acting on the example satellite.
Table 2. Orders of magnitude of the torques acting on the example satellite.
h (km)Order of Magnitude (N·m)
GravitationalAerodynamic RestoringAerodynamic DampingMagnetic
250 10 8 10 4 10 12 ε = 580 : 10 6 ε = 3900 : 10 5
700 10 8 10 7 10 15 ε = 0.58 : 10 9 ε = 1.35 : 10 8
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Aslanov, V.S.; Sizov, D.A. Attitude Dynamics of Spinning Magnetic LEO/VLEO Satellites. Aerospace 2023, 10, 192. https://doi.org/10.3390/aerospace10020192

AMA Style

Aslanov VS, Sizov DA. Attitude Dynamics of Spinning Magnetic LEO/VLEO Satellites. Aerospace. 2023; 10(2):192. https://doi.org/10.3390/aerospace10020192

Chicago/Turabian Style

Aslanov, Vladimir S., and Dmitry A. Sizov. 2023. "Attitude Dynamics of Spinning Magnetic LEO/VLEO Satellites" Aerospace 10, no. 2: 192. https://doi.org/10.3390/aerospace10020192

APA Style

Aslanov, V. S., & Sizov, D. A. (2023). Attitude Dynamics of Spinning Magnetic LEO/VLEO Satellites. Aerospace, 10(2), 192. https://doi.org/10.3390/aerospace10020192

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop