Next Article in Journal
One-Step Dry-Etching Fabrication of Tunable Two-Hierarchical Nanostructures
Previous Article in Journal
Improvement of Single Event Transient Effects for a Novel AlGaN/GaN High Electron-Mobility Transistor with a P-GaN Buried Layer and a Locally Doped Barrier Layer
Previous Article in Special Issue
Intelligent Navigation of a Magnetic Microrobot with Model-Free Deep Reinforcement Learning in a Real-World Environment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Actuators with Two Double Gimbal Magnetically Suspended Control Moment Gyros for the Attitude Control of the Satellites

by
Romulus Lungu
1,2,
Alexandru-Nicolae Tudosie
1,*,
Mihai-Aureliu Lungu
1 and
Nicoleta-Claudia Crăciunoiu
1
1
Faculty of Electrical Engineering, University of Craiova, 200585 Craiova, Romania
2
International Academy of Astronautics (IAA), 75016 Paris, France
*
Author to whom correspondence should be addressed.
Micromachines 2024, 15(9), 1159; https://doi.org/10.3390/mi15091159
Submission received: 28 June 2024 / Revised: 6 September 2024 / Accepted: 10 September 2024 / Published: 16 September 2024

Abstract

:
The paper proposes a novel automatic control system for the attitude of the mini-satellites equipped with an actuator having N = 2 DGMSCMGs (Double Gimbal Magnetically Suspended Control Moment Gyros) in parallel and orthogonal configuration, as well as a DGMSCMG-type sensor for the measurement of the satellite absolute angular rate. The proportional-derivative controller, designed based on the Lyapunov-functions theory, elaborates the control law according to which the angular rates applied to the servo systems for the actuation of the DGMSCMGs gyroscopic gimbals are computed. The gimbal’s angular rates create gyroscopic couples acting on the satellite in order to control its attitude with respect to the local orbital frame. The new proposed control architecture was software implemented and validated, and the analysis of the obtained results proved the cancellation of the convergence errors and excellent angular rate precision.

1. Introduction

The CMG (Control Moment Gyros) type actuators are indispensable drive systems for the controllers of the mini-satellites (satellites under 500 kg), being suitable for fast rotational maneuvers [1,2,3,4,5]. Comparing to the mechanical CMGs, the control moment gyros having rotors with active magnetic bearings (AMB-rotors) and single or double gimbals (SGMSCMGs or DGMSCMGs) possesses the advantages of zero-friction (eliminating the lubrication necessity), of lower noise, of small vibrations, as well as of increased longevity [6,7,8]. There are two kinds of the control moment gyros with two gimbals: (1) magnetically suspended DGCMGs with a magnetically suspended rotor and a stator part that is base-fixed; (2) variable speed double-gimbal control moment gyros (VSDGCMGs) [9]. The DGMSCMGs generate two gyroscopic torques, so a higher cumulative total torque, thus reducing the satellite’s mass and dimensions. An architecture involving a DGMSCMG with AMB allows the decoupling of its translation dynamics from the rotation one, as well as from gimbals’ rotational dynamics; thus, it allows the dynamics’ decoupling, some observers being required for state observation [1,2,10,11,12,13].
An important issue treated by the specialized scientific papers [14,15] is the resonance vibration control; another issue is the one associated with the design of the gimbals’ driving servo systems [7,10,16,17]. To avoid singularities using the null motion, with or without stored energy control, one proposed SGCMG type actuators, as well as actuators based on two or three DGCMGs, in orthogonal or parallel configuration [16,18,19,20,21]. Various satellites’ attitude control systems, based on both gimbals’ angular rates, and gyro-motors angular rates changing control have also been studied [5,15,16,22]. The topic of the satellites’ attitude control using DGMSCMGs, proportional-derivative control laws, proportional-integrator control laws, optimal control laws (with H2, H, H2/H optimal criteria), gain scheduling control laws (handling the external disturbances such as the coupling gyroscopic effects), integral-terminal-sliding mode control laws, or adaptive control laws, with or without neural networks (NNs), was treated by a lot of scientific papers, such as [1,2,3,9,15,16,23,24,25,26]. The attitude control was improved by means of various control architectures involving the modification of the DGMSCMG gimbals’ angular rates and/or the permanent change of the gyro-motors’ own rotation angular rates [5,27].
There are a lot of drawbacks to the papers mentioned above. The dynamics considered in paper [2] do not consider the rotor-gimbal interactions, the control being not sufficiently precise. The sliding-mode control approach is used to control the satellites’ attitude and to suppress the vibrations in [28], but the drawbacks of this control method involve the mandatory decoupling between the channels associated with the gimbals and the employed coordinate transformation. The parameter perturbation robust control method was employed in [29] to control the dynamics of MSCMGs, the obtained control scheme reaching its critical stability when the rotor angular rate is increased. Paper [19] proposed neural network-based fault diagnosis architectures for the estimation and the fault isolation of the SGCMGs equipping spacecraft by considering multiple actuators (SGCMG type); the attitude control is accurately achieved, but the employed cluster type includes SGCMGs, without decoupling between the rotation and the translation dynamics.
A new dynamics is obtained in work [25] for classical double-gimbal CMG systems, its strong point being the solving of the mismatched problem; also, the unwanted effects of the disturbances are properly cancelled via a backstepping control-based approach. However, this robust control method involves some shortcomings [9]: (a) the dynamics does not consider the rotor dynamics; (b) the nonlinearities of the systems should be completely known or must be estimated with very good accuracy; (c) instead of the gimbal angular rates, the inputs of the controller are the gimbals’ angular displacements. Paper [10] proposed some adaptive control laws to control the translational and rotational dynamics of the AMB-rotor, as well as the control of the gimbals’ angular rates from the component of a DGMSCMG; as the main drawback of this work, we can state the non-consideration of the interactions between the DGMSCMG and the satellite, on which this is placed as actuator and as subsystem of the automatic attitude control architecture.
Different from the above-mentioned control architectures, the present paper aims to design a novel attitude control scheme for mini-satellites by using an actuator with N = 2 DGMSCMGs (in parallel/orthogonal configuration) and a DGMSCMG-type sensor for the control of the satellite absolute angular rate. The design of the attitude control architecture is based on the Lyapunov theory. An important element is the actuator, consisting of two DGMSCMGs and a DGMSCMG-type sensor, as part of the actuator-sensor-satellite subsystem. To design the DGMSCMG-based actuator and sensor, each gyroscopic rotor’s translation dynamics was decoupled from its rotation dynamics and from gimbals’ rotation dynamics, followed by the design of the adaptive control of all three subsystems (for the rotors’ translation control, for the rotors’ rotation control, as well as for the gimbals’ rotation control) by using the dynamic inversion concept and neural networks.
The paper is structured as follows: In Section 2, one presents the elements of the attitude calculation approach, based on the quaternion method [30]. In Section 3, one presents the dynamic models of the actuators with N = 2 DGMSCMGs, in parallel and orthogonal configurations; the used DGMSCMGs have the architecture described in work [31]. The next section introduces the dynamic model of the DGMSCMG-type sensor for the satellite’s absolute angular rate. Section 5 presents a PD-type control law for the satellite’s attitude control, based on a Lyapunov-type function [32]. In the sixth part of the paper, one presents the results for the performed simulations for both types of actuators (parallel and orthogonal configurations), based on the Simulink-Matlab models.

2. Attitude Definition and Calculation

The attitude of the satellite (S), with its tied frame denoted by O x b y b z b   , with respect to the local orbital frame O x 0 y 0 z 0 ( O y 0 axis oriented in the direction that connects the origin of the initial frame O i x i y i z i to the satellite’s frame origin, O z 0 axis oriented in the orbital plan’s anti-normal direction, and O x 0 axis completing the orthogonal frame) is given by the Euler angles θ ,   φ ,   Ψ (see Figure 1). Rotating successively θ , φ   , and Ψ   , one obtains the satellite’s frame rotation matrix A   , with respect to the local orbital frame, as presented in [17]. In order to avoid the singularities for large Euler angles or a large number of calculations, one uses quaternions for the rotation matrix expression, as in [30]:
q q 4 T = e 1 sin ϕ 2 e 2 sin ϕ 2 e 3 sin ϕ 2 cos ϕ 2 T ,   q = q 1 q 2 q 3 T = e sin ϕ 2 ,
q T q + q 4 2 = 1 ,
where ϕ is the satellite’s rotation angle about the Euler axis, e = e 1 e 2 e 3 T the unit vector on Euler’s axis direction, while e 1 , e 2 , and e 3 are the vector’s director cosines with respect to the local orbital frame.
Quaternion’s kinematics is defined by the differential equations [30]:
q ˙ = 1 2 ω s × q + 1 2 q 4 ω s ,
q 4 = 1 2 ω s T q ,
or
q ˙ = 1 2 q × + q 4 I 3   ω s ,
q 4 = 1 2 q ω s ,
where ω s = ω s x ω s y ω s z T is the satellite’s angular rate vector with respect to the local orbital frame (with its components ω s x ,   ω s y ,   ω s z relative to the axes of the satellite-tied frame O x b y b z b ), I 3 3 × 3 unit matrix,
ω s × = 0 ω s z ω s y ω s z 0 ω s x ω s y ω s x 0 ,   ω b × = 0 ω b z ω b y ω b z 0 ω b x ω b y ω b x 0 ,   q × = 0 q 3 q 2 q 3 0 q 1 q 2 q 2 0   ,
with ω b x ,   ω b y ,   ω b z the components of the angular rate vector, relative to the satellite’s axes.
According to Figure 2, one defines the quaternions
q = q 1 q 2 q 3 T ,   λ = λ 1 λ 2 λ 3 T ,   ε = ε 1 ε 2 ε 3 T ,
q T q + q 4 2 = 1   ,       λ T λ + λ 4 2 = 1 ,   ε T ε + ε 4 2 = 1 ,
while the rotation matrices are [2]
A = A q , q 4 = q 4 2 q T q I 3 + 2 q T q 2 q 4 q × ,
B = B λ , λ 4 = λ 4 2 λ T λ I 3 + 2 λ T λ 2 λ 4 λ × ,
C = C ε , ε 4 = ε 4 2 ε T ε I 3 + 2 ε T ε 2 ε 4 ε × ,
from which it results in the formula B = A C   .
The relation between the absolute angular rate ω b , the relative angular rate ω s , and the transport angular rate ω 0 is
ω b = ω s + A q , q 4 ω 0 .
If the local orbital frame is considered the reference frame, i.e., O x R y R z R O x 0 y 0 z 0 , then ω 0 = 0 0 ω 0 T and the attitude matrix determined by Equation (8) becomes
A = 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 2 q 1 q 2 q 3 q 4 q 1 2 + q 2 2 q 3 2 + q 4 2 2 q 2 q 3 + q 1 q 4 2 q 1 q 3 + q 2 q 4 2 q 2 q 3 q 1 q 4 q 1 2 q 2 2 + q 3 2 + q 4 2   .
Identifying the elements of the matrix A   ( θ ,   φ ,   Ψ ) in work [17] with those associated with matrix A in Equation (12), one obtains the Euler angles as follows:
θ = atan   a 21 a 22 ,   θ = atan   a 23 cos Ψ a 33 ,   Ψ = atan   a 13 a 33 .
In Figure 3, the block diagram of the satellite’s attitude calculation system is depicted. The components of the satellite’s absolute angular rate ω b x ,   ω b y ,   ω b z are measured by using rate gyros having their sensitivity axes along the axes of the satellite-tied frame, while the components of the satellite’s relative angular rate are computed by means of the block diagram depicted in Figure 3 and by using the formulas.
ω s x = ω b x + 2 ω 0 q 1 q 3 q 2 q 4   ,
ω s y = ω b y + 2 ω 0 q 2 q 4 + q 1 q 3   ,
ω s z = ω b z ω 0 q 1 2 + q 2 2 q 3 2 q 4 2 .

3. Dynamic Models of the Actuators with N = 2 DGMSCMGs

3.1. DGMSCMGs with Paralel Architecture

In the satellite’s attitude automatic control system, the subsystem consisting of two DGMSCMGs (in parallel/orthogonal configuration) acts as an actuator. Each one of these reacts at the angular rates ( ω i j and ω e j , where j is j = 1 for DGMSCMG 1 and j = 2 for DGMSCMG 2), computed by the attitude controller and applied to the inner and outer DGMSCMG’s gimbals by means of a servo system.
Figure 4a [31] depicts the simplified DGMSCMG structure (with j = 1   ,   2 ¯ ); for σ i 1 = σ e 1 =   = σ i 2 = σ e 2 = 0   , it is positioned with its inner gimbal axis (i) along the O x b axis, with its outer gimbal axis along the O y b axis, and with the gyroscopic rotor’s axis (r), i.e., the axis of its kinetic couple K 0   , along the O z b axis; this way, the outer gimbals are parallel.
Figure 4b depicts the gimbals’ angles and angular rates, the components of the satellite’s relative angular rate vector, as well as the gyroscopic rotors’ precession angular rates ( ϕ ˙ x r j and ϕ ˙ y r j   , with j = 1 , 2 ¯ ). Figure 4c depicts the graph of the interdependences between the rotor, the gimbals, and the base, which matches Figure 2. Having in mind that the DGSMCMGs are PD-type automatic control subsystems for the compensation of the linear/angular deviations and rates [10], the angular rate of the gyroscopic rotor can be written relative to the inner gimbal as: ω r j = α ˙ j + β ˙ j 0   ; therefore, the frames O x r j y r j z r j and O x i j y i j z i j are overlapped and, consequently, the rotation matrix of the rotor-tied frame with respect to the inner gimbal-tied frame becomes A i j r j I 3 .
According to Figure 4b,c, the absolute angular rates of the rotor and inner gimbals of the DGSMCMGs relative to the local orbital frame are the resultants of the relative angular rates σ ˙ i and σ ˙ c (with respect to the base-tied frame) and the transport angular rate (of the base relative to the local orbital frame) ω s . Thus, the coordinates of the outer gimbal (e) relative to the inner one (i) and the base (b) are
ϕ ˙ j = σ ˙ i j + σ ˙ e j + ω s = σ ˙ i j + σ ˙ e j + ω s x + ω s y + ω s z ,
σ ˙ i j = σ ˙ i j 0 0 T ,   σ ˙ e j = 0 σ ˙ e j 0 T ,   j = 1 ,   2 ¯ .
Projecting Equation (15) on the axes of the O x r j y r j z r j frames (overlapped over O x i j y i j z i j ) and considering Figure 4b, one yields
ϕ ˙ x r j ϕ ˙ x i j = σ ˙ i j + ω s x cos σ e j + ω s z sin σ e j   ,
ϕ ˙ y r j ϕ ˙ y i j = σ ˙ e j + ω s y cos σ i j ω s x sin σ e j ω s z cos σ e j sin σ i j   ,
where ω s x ,   ω s y , and ω s z are obtained via Equation (14).
According to Figure 4b, between the angular rates expressed with Equation (17) and the ones provided by the attitude controller, there are the next mathematical relations:
ω i j = ϕ ˙ x r j ϕ ˙ x i j ,   ω e j = ϕ ˙ y r j / cos σ i j ϕ ˙ y i j / cos σ i j ,   j = 1 , 2 ¯ .
By using Equations (17) and (18), one gets the relative angular rates applied to the servo systems driving the gyroscopic gimbals
σ ˙ i j e = ω i j ω s x cos σ e j + ω s z sin σ e j ,
σ ˙ e j e = ω e j ω s y + ω s x sin σ e j ω s z cos σ e j   tan σ i j   ,
where ω i j and ω e j are the components of the algebraic vectors ω i j and ω e j   , that, according to Figure 4b, have the following expressions:
ω i j = ω i j 0 0 T ,   ω e j = 0 ω e j 0 T ,   j = 1 , 2 ¯ .
By means of these forms, one can express the matrix
ω i j × = 0 0 0 0 0 ω i j 0 ω i j 0 ,   ω e j × = 0 0 ω e j 0 0 0 ω e j 0 0   .
Each of the two DGSMCMGs has the form of the one modeled and designed in paper [10]. Thus, the j-th DGSMCMG consists of four subsystems: (1) the automatic control subsystem (with a PID-type control law) for the control of the coordinates x r j and y r j (linear displacements) associated to the gyroscopic rotor; (2) the adaptive control subsystem for the control of the gyroscopic rotor’s angular displacements (consisting of a PD-type dynamic compensator, a linear observer, and a neural network modeling an adaptive control law component); (3) the servo system for the adaptive control of the gyroscopic gimbals’ angular rates σ ˙ i j and σ ˙ e j (also consisting of a PD-type dynamic compensator, a linear observer, and a neural network for the adaptive control law modeling). The gimbals are driven by electrical motors having their shafts’ axes co-linear with the gimbals’ axes (favorable arranged so that σ ˙ i j and σ ˙ e j are following the calculated angular rates σ ˙ i j e and σ ˙ e j e   , given by Equation (19)); (4) the subsystem modeling the interaction between the j-th DGSMCMG and the satellite.
The first three subsystems, presented in detail in work [10], have simplified forms, depicted in Figure 5a–c, while the last subsystem is presented in Figure 5d. For the first subsystem (Figure 5a), the command vector is u 1 j = i x j i y j T , where i x j and i y j are the correction currents, applied to the stator’s coils of the magnetic bearing (two coils in series on the gyroscopic rotor’s axes O x r and O y r ). For the second subsystem (the one in Figure 5b), the command vector is u 2 j = i α j i β j T , with i α j and i β j the correction currents applied to the same coils. For the third subsystem (the one depicted in Figure 5c), the command vector is u 3 j = i x i 1 i y e j T , where i x i 1 and i y e j are the currents applied to the command coils of the correction motors.
In Figure 5d, the term M k = M g represents the input of the satellite’s model, described by the equation:
I b ω ˙ b + ω b × I b ω b = M k + M p , M k = j = 1 2 M k j .
In the above formula, M p is the total disturbing moment (applied to the satellite), I b the satellite’s matrix of the inertia moments, considered together with the actuator ( j = 1 , 2 ¯ ) and the angular rate sensor j = 3 , given by the formula in [1]:
I b = J b + j = 1 3 A b e j T J e j A b e j + A b i j T I i j A b i j ,
where J b is the satellite’s matrix of inertia moments (without actuator and sensor), while J i j and J e j are the matrices of the inertia moments for the inner and outer gimbals of the j-th DGMSCMG, respectively; I i j = J i j + J r j and J r j are the matrices of inertia moments associated to the gyroscopic rotor, while A b i j = A i j b T and A b e j = A e j b T are the rotation matrix of the inner and outer gimbals, related to the satellite base. The components for j = 3 are for the DGMSCMG used as the angular rate sensor of the satellite’s base.
For computing the rotation matrices, one can use Figure 4b. The matrices that express the rotation of the outer gimbal relative to the base and the rotation of the inner gimbal with respect to the outer gimbal can be respectively determined as:
x e j y e j z e j = A b e j x b y b z b = cos σ e j 0 sin σ e j 0 1 0 sin σ e j 0 cos σ e j x b y b z b ,
x i j y i j z i j = A e j i j x e i y e i z e i = 1 0 0 0 cos σ i j sin σ i j 0 sin σ i j cos σ i j x b y b z b ,
A i j b A r j b = A b i j T = A e j i j A b e j T ,   A e j b = A b e j T ,   A i j e j = A e j i j T .
According to Figure 4c, the projections of the kinetic moment K 0 j = 0 0 K 0 T on the axes of the base-tied, on the axes of the inner gimbal-tied frame, and on the axes of the outer gimbal-tied frame have respectively the forms:
K r b = A r j b K 0 j ,   K r i j = A r j i j K 0 j ,   K r e j = A r j e j K 0 j ,
A r j b = A b r j T = A i j r j A e j i j A b e j T I 3 A e j i j A b e j T = A b i j T ,   A r j i j I 3 , A r j e j = A e j r j T = A i j r j A e j i j T I 3 A e j i j T = A e j i j T = A i j e j   ;
Thus, we write:
K r b = A i j b K 0 j , K r i j = K 0 j , K r e j = A r j e j K 0 j ,   K 0 j = K 0 , j = 1 , 2 ¯ .
The angular rates ω b , ω i j , ω e j generate the gyroscopic couples:
M g b = K r b × ω b = ω b K r b × = ω b A i j b K 0 j = M k b ,
M g i j = K r i j × ω i j = K 0 j × ω i j = M k i j ,
M g e j = K r e j × ω e j = A i j e j K 0 j × ω e j = M k e j ,   j = 1 ,   2 ¯ .
Projecting these torques on the axes of the O x b y b z b frame and summing them up, one obtains the following equivalent equations:
M k j = M g j = ω b × A i j b K 0 j A j j b K 0 j × ω i j A e j b A i j e j K 0 j × ω e j ,
M k j = M g j = ω b × A i j b K 0 j + A j j b ω i j × K 0 j × + A e j b ω e j × A i j e j K 0 j .
The total torque created by the two DGMSCMGs of the actuator has one of the following equivalent forms:
M k = M g = ω b × j = 1 2 A i j b K 0 j j = 1 2 A i j b K 0 j × ω i j j = 1 2 A e j b A i j e j K 0 j × ω e j ,
M k = M g = ω b × j = 1 2 A i j b K O j + j = 1 2 A i j b ω i j v K O j + j = 1 2 A e j b ω e j × A i j e j K O j .

3.2. DGMSCMGs with Orthogonal Architecture

The two DGMSCMGs of the actuator have the axes of the outer gimbals perpendicular; DGMSCMG 1 has the axis of the outer gimbal parallel to the O y b axis (see Figure 4b), while DGMSCMG 2 has the axis of its outer gimbal parallel to O z b axis. In Figure 6, we present the frames of the rotors, of the inner and outer gimbals, as well as the angular rates of the gyroscopic gimbals for both DGMSCMGs (in Figure 6a for DGMSCMG 1 and in Figure 6b—for DGMSCMG 2).
For DGMSCMG 1, Equations (15)–(25) remains valid, while for DGMSCMG 2, according to Figure 6b, one obtains:
ϕ ˙ y r 2 ϕ ˙ y i 2 = ω i 2 = σ ˙ i 2 + ω s x sin σ e 2 + ω s y cos σ e 2 ,
ϕ ˙ z r 2 ϕ ˙ z i 2 = ω e 2 cos σ i 2 = σ ˙ e 2 + ω s z   cos σ i 2 + ω s x cos σ e 2 + ω s y sin σ e 2   sin σ i 2   ,
which lead to new expressions for the calculated angular rates of the gimbals with respect to the base:
σ ˙ i 2 c = ω i 2 + ω s x sin σ e 2 ω s y cos σ e 2 ,
σ ˙ i 2 c = ω i 2 + ω s x sin σ e 2 ω s y cos σ e 2   .
The expressions of the kinetic moment vector, of the angular rate vector, and of the matrices ω i 2 × ,   ω e 2 × are
K 02 = K 0 0 0 T ,   σ ˙ i 2 = 0 σ ˙ i 2 0 T , σ ˙ e 2 = 0 0 σ ˙ e 2 T ,
ω i 2 = 0 ω i 2 0 T ,   ω e 2 = 0 0 ω e 2 T ,
ω i 2 × = 0 0 ω i 2 0 0 0 ω i 2 0 0 ,   ω e 2 × = 0 ω e 2 0 ω e 2 0 0 0 0 0 .
According to Figure 6b,
x e 2 y e 2 z e 2 = cos σ e 2 sin σ e 2 0 sin σ e 2 cos σ e 2 0 0 0 1 A b e 2   x b y b z b ,   x i 2 y i 2 z i 2 = cos σ i 2 0 sin σ i 2 0 1 0 sin σ i 2 0 cos σ i 2 A e 2 i 2   x e 2 y e 2 z e 2 .
Equations (25) and (39) are the same for j = 2, too. At the same time, one maintains the same calculation formulas (i.e., Equations (30)–(33), for j = 1, 2), and also Equation (23). Figure 7 presents the subsystems associated with the DGMSCMGs of the actuator in orthogonal configuration.

4. Dynamics of DGMSCMG 3 Sensor for Satellite Absolute Angular Rate Measurement

The disturbed satellite, rotating with ω s relative angular rate, should be returned to its initial position via an angular rate ω s , produced by the moment M k = M g , applied by the actuator to the satellite. M k depends only on ω b the absolute angular rate, in order to return the satellite (aperiodic according to Equation (22)) to its angular rate ω b = A q s , q 4 s ω 0 , where q s and q 4 s denote the stabilized values of the quaternion (33), M k = M k ω b , ω i , ω e , with ω i = ω i 1 T ω i 2 T T and ω e = ω e 1 T ω e 2 T T . Therefore, the satellite will rotate with a different angular rate ω ^ b ω b ; thus, an angular rate sensor (DGMSCMG 3) is required for the measurement of ω ^ b ; this sensor is connected to the satellite and to which no angular rates around the outer and inner gimbals’ axes are applied ω i 3 = ω e 3 = 0 . This sensor is placed in such a way that the axis of its outer gimbal is parallel to the axis of the outer gimbal of DGMSCMG 1.
For the sensor, the Equation (19) becomes (for i = 3):
σ ˙ i 3 c = ω s x cos σ e 3 + ω s z sin σ e 3 ,
σ ˙ e 3 c = ω s y + ω s x sin σ e 3 ω s z cos σ e 3   tan σ i 3 .
The sensor is modeled by the equations:
M k s = M g s = ω ^ b × A i 3 b K 03 = A i 3 b K 03   ω ^ b ,   K 03 = 0 0 K 0 T ,
ω ^ b = ω ^ b x ω ^ b y ω ^ b z ,   ω ^ b × = 0 ω ^ b z ω ^ b y ω ^ b z 0 ω ^ b x ω ^ b y ω ^ b x 0   .
Consequently, two torques are acting on the satellite: M k and M k s   , namely Δ M k =   = M k M k s . The subsystems of the DGMSCMG 3 sensor are depicted in Figure 8. As long as the angular rates ω e and ω i are null, the system depicted in Figure 9 is stabilized, i.e., M k = M k ω b M k s ω ^ b   , which means that Δ M k 0 and, thus, ω ^ b ω b = ( 11 ) A q s , q 4 s ω 0 and ω s 0 .

5. Satellite’s Attitude Control Using PD-Type Controller

To design the controller, a Lyapunov function is chosen as in [32]:
V = 1 2 ω s e T I b ω s e + 2 k p ln 1 + q e T q e ,   k p > 0 ,
ω s e = ω s d ω s ,   q e = M q d , q 4 d q q 4 T ,
where ω s d is the desired angular rate of the satellite relative to the local orbital frame, q d the desired quaternion of the satellite, q the quaternion expressing the satellite’s attitude relative to the local orbital frame, q e the error quaternion, while M q d , q 4 d is a matrix having the following form:
M q d , q 4 d = q 4 d q 3 d q 2 d q 1 d q 3 d q 4 d q 1 d q 2 d q 2 d q 1 d q 4 d q 3 d q 1 d q 2 d q 3 d q 4 d .
The error quaternion is the solution of the differential equation.
q ˙ e = F q e ω s e ,
where F q e and q e have the following forms:
F q e = 1 2 1 2 1 + q e T q e I 3 + q e × ,
q e = q e 1 q e 2 q e 3 ,   q e × = 0 q e 3 q e 2 q e 3 0 q e 1 q e 2 q e 1 0 .
The time derivative of the Lyapunov function is successively obtained as follows:
V ˙ = ω s e T I b ω ˙ s e + 4 k p 1 + q e T q e q e T q ˙ e = ω s e T I b ω ˙ s e + 4 k p q e T 1 + q e T q e F q e   ω s e = ω s e T I b ω ˙ s e + 4 k p q e T 1 + q e T q e 1 2 1 2 1 + q e T q e I 3 + q e x   ω s e = ω s e T I b ω ˙ s e + k p q e T ω s e + 2 k p q e T q e × 1 + q e T q e ω s e = ω s e T I b ω ˙ s e + k p ω s e T q e ;
because q e T q e x = 0   , one yields
V ˙ = ω s e T I b ω ˙ s e + k p ω s e T q e .
Imposing V ˙ = k d ω s e T ω s e , with k d > 0 , one obtains
V ˙ = ω s e T I b ω ˙ s e + k p ω s e T q e = k d ω s e T ω s e < 0 ,
i.e., the following inequality results:
ω s e T I b ω ˙ s e + k d ω s e + k p q e < 0 ;
V ˙ < 0 is the stability condition for the closed-loop system depicted in Figure 9.
Replacing ω ˙ s e with its expression ω ˙ s e = ω ˙ s d ω ˙ s in the following equation
I ω ˙ s e + k d ω s e + k p q e = 0   ,
one obtains
I b ω ˙ s d + k d ω s e + k p q e = I b ω ˙ s ,
where the term I b ω ˙ s comes from the satellite dynamics equation.
The term ω ˙ b is computed with the derivative of Equation (11), namely: ω ˙ b = ω ˙ s + A ˙ q , q 4   ω 0 =   = ω ˙ s ω s x A q , q 4 ω 0 = ( 12 ) ω ˙ s + ω 0 ω s x c o l 3 A ; introducing it into Equation (22), a new form is obtained:
I b ω ˙ s = ( 22 ) M k ω b × I b ω b ω 0 I b ω s × c o l 3 A + M p ,
where the term c o l 3 A represents the third column vector of matrix A, while M k has the form (32) or (33), i.e.,
M k = ω b × j = 1 2 A i j K 0 j j = 1 2 A i j b K 0 j × ω i j j = 1 2 A e j b A i j e j K 0 j × ω e j = ω b × j = 1 2 A i j b K 0 j + M k c ,
M k c = j = 1 2 A i j b K 0 j × ω i j j = 1 2 A e j b A i j e j K 0 j × ω e j .
Invoking Equations (55) and (56), one obtains
I b ω ˙ s = ω b × I b ω b + j = 1 2 A i j b K 0 j × ω 0 I b ω s × c o l 3 A + M p M k c .
One may use the following notations:
Q i j = A i j b K 0 j × ,   Q e j = A e j b A i j e j K 0 j × ,   Q j = Q i j Q e j ,   j = 1 , 2 ,
where Q i j and Q e j are (3 × 3) matrices, while Q j is a (3 × 6) matrix. With these notations, Equation (57) becomes
M k c = j = 1 2 Q i j ω i j j = 1 2 Q e j ω e j = j = 1 2 Q i j Q e j   ω i j T ω e j T T = j = 1 2 Q j ω g j ,
with
Q j = Q i j Q e j ,   ω g j = ω i j T ω e j T T ,   j = 1 ,   2 ¯ .
The expression (60) might be rewritten as follows:
M k c = Q 1 ω g 1 Q 2 ω g 2 = Q 1 Q 2   ω g 1 T ω g 2 T Τ = Q ω g ,
with
Q = Q 1 Q 2 ,   ω g = ω g 1 T ω g 2 T T = ω i T ω e T T ,   ω i = ω i 1 T ω i 2 T T ,   ω e = ω e 1 T ω e 2 T T .
By using the value of the torque M k c , computed by the attitude controller, as well as Equation (62), one obtains
ω g = ω i T ω e T T = Q + M k c ,   Q + = Q T Q 1 Q T ,
where Q + is the pseudo-inverse of the matrix Q .
Replacing Equation (58) into Equation (54) and using expression (62) for M k c , one obtains
I b ω ˙ s d + k d ω s e + k p q e + ω b × I b ω b + j = 1 2 A i j b K 0 j × + ω 0 I b ω s × c o l 3 A M p = M k c ,
or
M k c / M ^ k c = M k c ,
where M k c / is the total disturbance moment with expression
M k c / = ε = ω b × I b ω b + j = 1 2 A i j b K 0 j x + ω 0 I b ω s × c o l 3 A M p
and
M ^ k c = I b ω ˙ s d + k d ω s e + k p q e ,
This last expression is being modeled by the PD-type attitude controller.
Remark 1. 
The torque M k c is computed for two configurations (orthogonal and parallel), each time with respect to the angular rates (that must be applied to DGMSCMGs’ gimbals) and the rotation matrices associated to the gimbals (relative to the base). By means of the dynamic inversion (i.e., the pseudo-inverse of the matrix Q ) and the torque M k c (provided by the attitude controller), the angular rates ω i and ω e are computed and the imposed angular rates of the DGMSCMGs’ gimbals are obtained; these angular rates represent the inputs of the servo-systems acting the gimbals and producing the total torque M g = M k applied to the base in order to obtain its rotation.
The architecture of the satellite’s attitude controller is depicted in Figure 9. The reference model is a first-order command filter, having the transfer matrix H m s = 1 s + 1 I 3 . It provides both the desired angular rate ω s d and its time derivative ω ˙ s d ; the angular acceleration ω ˙ s d is considered in the control design in order to obtain ω ˙ s = ω ˙ s d .
Remark 2. 
The attitude control of the satellites, via the P.D.-type controller, is equivalent to the P.I.-type control of the angular rate ω s , relative to the local orbital frame and, consequently, equivalent to the P.I.-type control of the absolute angular rate ω b . The attitude controller calculates the imposed command torque M k c for the satellite’s orientation and its attitude’s stabilization, based on Lyapunov functions theory. For the satellite’s rotation, the M k torque must be applied; it should tend to M k c .

6. Numerical Simulations

One has performed a study concerning the dynamics of the system depicted in Figure 9. The numerical values for all three DGMSCMGs are those used in [10] (see Table 1), to which one has added the initial values of the parameters, which might be remarked directly from the graphics in Figure 10 and Figure 11—representing the dynamic characteristics of the system in Figure 9, determined using its Matlab-Simulink model.
The graphics in Figure 10a and in Figure 11a present the components of the satellite’s relative angular rate vector ω s , the components of the error relative angular rate vector ω s e , as well as the components of absolute angular rate vectors ω b and ω ^ b .
Figure 10b,c and Figure 11b,c present the characteristics expressing the satellite’s attitude with respect to the local orbital frame, namely the quaternion, the error quaternion, and the Euler angles, respectively. Figure 10d and Figure 11d depict the components of the torque vectors M ^ k c , M k c / , M k c with respect to the satellite-tied frame, as well as the components of the actuator’s and sensor’s output vectors M k , M k s , modeled by the systems given in Figure 5d and Figure 7d, respectively.
Figure 10e and Figure 11e present the characteristics, ω i j ,   ω e j ,   j = 1   , 2 ¯ , i.e., the angular rates applied to the gimbals of the actuator-type DGMSCMGs. In Figure 10f and in Figure 11f, there are graphically represented the dynamic characteristics of the models of both actuator’s DGMSCMGs, shown in a simplified manner in Figure 5a,b,c and Figure 7a,g, respectively, having their complete structure provided in work [10]. Thus, for each DGMSCMG of the actuators, one has represented:
  • the components of the vectors of the linear displacements associated to the rotors of the gyro-motors, with respect to the precession axes, i.e., y 1 j ,   j = 1 , 2 ¯ ;
  • the components of the error vectors y ˜ 1 j , with y ˜ 1 j = y 1 c y 1 j ,   j = 1 , 2 ¯ ;
  • the components of the pseudo-command vectors for the control of the gyroscopic rotor’s linear displacements, namely V 1 j ,   j = 1 , 2 ¯ ;
  • the components of the command vectors u 1 j ,   j = 1 , 2 ¯ , i.e., the currents applied to the correction stator coils of the gyroscopic rotors’ magnetic bearings;
  • the components of the vectors y 2 j ,   y ˙ 2 j ,   y ˜ 2 j ,   j = 1 , 2 ¯ , associated with the gyroscopic rotors, namely the vectors consisting of the precession angles, the precession angular rates, and the precession angles’ errors, respectively;
  • the components of the control vectors u 2 j ,   j = 1 , 2 ¯ , and the components of the pseudo-command vectors V 2 j ,   j = 1 , 2 ¯ , all of them associated with the subsystems for the adaptive control of the precession angles; the vectors u 2 j ,   j = 1 , 2 ¯ , contain the command currents applied to the same coils of the bearings;
  • the components of the vectors y 31 j c (calculated angular rates), y 31 j (rotation angular rates), y ˜ 31 j (error angular rates), and y 31 i j (rotation angles of the gimbals);
  • the components of the vectors V 3 j and u 3 j ,   j = 1 , 2 ¯ , associated to the command servo systems for the gimbals for the two DGMSCMGs; u 3 j contain the currents applied to the command coils of the correction motors on the inner and outer gimbals.
In the Figure 10g and Figure 11g, there are represented the dynamic characteristics of the DGMSCMG 3 gyro sensor, characteristics that are similar to those in Figure 10f and Figure 11f, obtained for the DGMSCMGs belonging to the actuator.
Comparing the dynamic characteristics of the attitude control system with N = 1 DGMSCMG to those of the attitude control system with N = 2 DGMSCMGs (in parallel/orthogonal configuration), one remark is that the architectures with actuators having N = 1 DGMSCMG have slightly higher overshoots, are more oscillating, and are a little bit slower.
Remark 3. 
Gyroscopic actuators with N DGMSCMGs (N = 1, 2, or 3), in orthogonal or parallel configuration, compared to similar DGVSCMG-based actuators [8,15,33,34], have the advantage of zero friction in their gyroscopic bearings, which makes unnecessary both the lubrication and the control of the rotors’ speed(s) and stored kinetic energy; thus, they have superior performances. However, DGMSCMGs must be equipped with systems for the automatic control of the gyroscopic rotors’ linear and angular displacements.
Remark 4. 
In some papers approaching the same topic, the control of DGMSCMGs’ dynamics is based on the decoupling of the rotors’ translation dynamics from their rotation dynamics. For the control of rotors’ and gimbals’ nonlinear rotation dynamics, one has used differential geometry’s principles [11,12]; in the present paper, the control is based on the coupling of the gyroscopic rotors’ rotation dynamics to the gimbals’ rotation dynamics. Comparing the dynamic and static characteristics of the described actuators’ structures (with N = 2 DGMSCMGs in orthogonal/parallel configurations) using the differential geometry’s principles [11,12] to those in the present paper (using the dynamic inversion concept and the adaptive control based on neural networks—a novel element in the field of gyro-systems), the superiority of the latter is noteworthy. This second control technique leads to superior and less oscillating dynamic characteristics, with smaller overshoot, convergence time, and static errors; that means superior performances, especially from the precision of the satellite’s attitude orientation/stabilization point of view.
An important novelty element in this paper (the use of a supplementary DGMSCMG-type sensor on the actuator’s feedback) contributes to the improvement of the above-mentioned performances. The gyroscopic torque generated by this sensor and applied to the satellite compensates the positioning error of the satellite due to the error of the satellite’s angular rate canceling during its stabilization.
Another advantage of the above-described control architecture is related to the adaptive control law, which uses a small number of sensors, namely sensors for gyroscopic rotors’ total linear displacements and for gimbals’ angular rates, respectively.

7. Conclusions

This work first introduced some calculation elements concerning the attitude control of the mini-satellites by using the quaternion theory and Euler angles. Then, there are presented some models for two architectures of the actuator, with N = 2 DGMSCMGs, in parallel and orthogonal configuration, additionally considering the third DGMSCMG as a sensor for the measurement of the satellite’s absolute angular rate ω b . One has calculated the total gyroscopic torques created by the two DGMSCMGs of the actuator relative to the angular rate ω b   , the angular rates ω i j ,   ω e j ,   j = 1 , 2 ¯ (computed by the attitude controller and applied to the servo-systems belonging to the DGMSCMGs for the actuation of the gimbals), as well as the kinetic torques induced to the satellite by the kinetic moments of the gyroscopic rotors and gimbals. A novel automatic attitude control system has been obtained, its main subsystem being a PD-type controller. For this adaptive control scheme, by means of complex simulations in the Matlab-Simulink environment, one has obtained the dynamic characteristics for both architectures (the one with two DGMSCMGs in parallel configuration and the one with two DGMSCMGs in orthogonal configuration).

Author Contributions

Conceptualization, R.L.; methodology, R.L., A.-N.T. and M.-A.L.; software, M.-A.L. and A.-N.T.; validation, M.-A.L. and A.-N.T.; formal analysis, R.L.; investigation R.L. and A.-N.T.; resources, R.L. and A.-N.T.; data curation, A.-N.T. and M.-A.L.; writing—original draft preparation, R.L., A.-N.T. and N.-C.C.; writing, review, and editing, A.-N.T. and N.-C.C.; visualization, R.L., A.-N.T. and M.-A.L.; supervision, R.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Zhon, J.; Zhon, D. Spacecraft Attitude Control with Double—Gimbaled Control Moment Gyroscopes. In Proceedings of the IEEE International Conference on Robotics and Biomimetics, Sanya, China, 5–8 December 2017. [Google Scholar]
  2. Zheng, S.; Han, B. Investigations of an Integrated Angular Velocity Measurement and Attitude Control System for Spacecraft Using Magnetically Suspended Double—Gimbal CMGs. Space Res. 2013, 51, 2216–2228. [Google Scholar] [CrossRef]
  3. Peng, C.; Fang, J.; Xu, S. Composite anti-disturbance controller for magnetically suspended control moment gyro subject to mismatched disturbances. Nonlinear Dyn. 2015, 79, 1563–1573. [Google Scholar] [CrossRef]
  4. Narita, M.; Chen, G.; Takami, I. Gain scheduling controller synthesis for active magnetic beaaring based on parameter dependent LMI with convexity condition. In Proceedings of the IEEE International Conference 2nd on Control and Robotics Engineering, Bangkok, Thailand, 1–3 April 2017; pp. 63–67. [Google Scholar]
  5. Higashiyama, D.; Shoji, Y.; Satoh, S.; Jikuya, I.; Yamada, K. Attitude control for spacecraft using pyramid-type variable-speed control moment gyros. Acta Astronaut. 2020, 173, 252–265. [Google Scholar] [CrossRef]
  6. Shelke, S. Controllability of Radial Magnetic Bearing. In Proceedings of the 3rd International Conference of Innovation Automation and Mechatronics Engineering—ICIAME, Vallabh Vidhyanagar, India, 5–6 February 2016; pp. 106–113. [Google Scholar]
  7. Han, B.; Ma, J.; Li, H. Reserch on Nelinear Friction Compensation of Harmonic Drive in Gimbal Servo—System of DGCMG. Int. J. Control Autom. Syst. 2016, 14, 779–786. [Google Scholar] [CrossRef]
  8. Sasaki, T.; Shimomura, T.; Schaub, H. Robust attitude control using a double-gimbal variable-speed control moment gyroscope. J. Spacecr. Rocket. 2018, 55, 1235–1247. [Google Scholar] [CrossRef]
  9. Lungu, R.; Lungu, M.; Efrim, C. Attitude adaptive control of satellites using double-gimbal magnetically suspended control moment gyros. Aerosp. Sci. Technol. 2022, 126, 107652. [Google Scholar] [CrossRef]
  10. Lungu, R.; Lungu, M.; Efrim, C. Adaptive control of DGMSCMG using dynamic inversion and neural networks. Adv. Space Res. 2021, 68, 3478–3494. [Google Scholar] [CrossRef]
  11. Xiaocen, C.; Maoyin, C. Precise Control of Magnetically Suspended Doble-Gimbal Control Moment Gyroscope Using Differential Geometric Decopling Method. Chin. J. Aeronaut. 2013, 26, 1017–1028. [Google Scholar]
  12. Lungu, R.; Efrim, C.; Zheng, Z.; Lungu, M. Dynamic Inversion Based Control of DGMSCMGs. Part2: Control Architecture Design and Validation. In Proceedings of the IEEE International Conference on Applied and Theoretical Electricity, Craiova, Romania, 27–29 May 2021. [Google Scholar]
  13. Cui, P.; Cui, J.; Yang, Q.; Zheng, S. The Coupling Characteristic Investigation of Double-Gimbal Magnetically Suspended Control Moment Gyro Used on Agile Maneuver Spacecraft. Math. Probl. Eng. 2015, 2015, 78145. [Google Scholar] [CrossRef]
  14. Li, H.; Yang, S.; Ren, H. Dinamic Decoupling Control of DGCMG Gimbal System via State Feedback Linearization. Mechatronics 2016, 36, 127–135. [Google Scholar] [CrossRef]
  15. Sasaki, T.; Shimomura, T.; Pullen, S.; Schaub, H. Attitude and vibration control with double-gimbal variable-speed control moment gyros. Acta Astronaut. 2018, 152, 740–751. [Google Scholar] [CrossRef]
  16. Lungu, M. Backstepping Control Method in Aerospace Engineering; Academica Greifswald Publisher: Rostock, Germany, 2022. [Google Scholar]
  17. Lungu, M.; Lungu, R. Adaptive neural network-based satellite attitude control by using the dynamic inversion technique and a VSCMG pyramidal cluster. Complexity 2019, 2019, 1645042. [Google Scholar] [CrossRef]
  18. Liu, F.; Gao, F.; Zhang, W.; Zhang, B.; He, J. The optimization design with minimum power for variable speed control moment gyroscopes with integrated power and attitude control. Aerosp. Sci. Technol. 2019, 88, 287–297. [Google Scholar] [CrossRef]
  19. Yuandong, L.I.; Qinglei, H.U.; Xiaodong, S.H. Neural network-based fault diagnosis for spacecraft with single-gimbal control moment gyros. Chin. J. Aeronaut. 2021, 35, 261–273. [Google Scholar]
  20. Sasaki, T.; Shimomura, T.; Kanata, S. LPV control and singularity avoidance of a spacecraft with DGCMGs. IFAC-PapersOnLine 2016, 49, 152–157. [Google Scholar] [CrossRef]
  21. Walker, J.M.; Culbertson, H.; Raitor, M.; Okamura, A. Haptic orientation guidance using two parallel double-gimbal control moment gyroscopes. IEEE Trans. Haptics 2017, 11, 267–278. [Google Scholar] [CrossRef]
  22. Leeghim, H.; Kim, D. Singularity-robust control moment gyro allocation strategy for spacecraft attitude control in the presence of disturbances. Aerosp. Sci. Technol. 2021, 119, 107178. [Google Scholar] [CrossRef]
  23. Lungu, R.; Lungu, M.; Ioan, M. Determination and control of the satellites’ attitude by using a pyramidal configuration of four control moment gyros. In Proceedings of the 12th IEEE International Conference on Informatics in Control, Automation and Robotics, Colmar, France, 21–23 July 2015; pp. 448–456. [Google Scholar]
  24. Isidori, A. Nonlinear Control Systems; Springer: Berlin/Heidelberg, Germany, 1995. [Google Scholar]
  25. Lungu, M. Control of double gimbal control moment gyro systems using the backstepping control method and a nonlinear disturbance observer. Acta Astronaut. 2021, 180, 639–649. [Google Scholar] [CrossRef]
  26. Acquatella, B.P.; Briese, L.E.; Schnepper, K. Guidance command generation and nonlinear dynamic inversion control for reusable launch vehicles. Acta Astronaut. 2020, 174, 334–346. [Google Scholar] [CrossRef]
  27. Ma, Z.Q.; Huang, P.F. Discrete-time sliding mode control for deployment of tethered space robot with only length and angle measurement. IEEE Trans. Aerosp. Electron. Syst. 2020, 56, 585–596. [Google Scholar] [CrossRef]
  28. Li, H.; Yu, J. Anti-disturbance control based on cascade ESO and sliding mode control for gimbal system of double gimbal CMG. IEEE Access 2019, 8, 5644–5654. [Google Scholar] [CrossRef]
  29. Han, B.; Chen, Y.; Zheng, S.; Li, M.; Shi, Y. Robust Control for a magnetically suspended control moment gyro with strong gyroscopic effects. In Proceedings of the IECON Annual Conference of the IEEE Industrial Electronics Society, Washington, DC, USA, 21–23 October 2018; pp. 2440–2446. [Google Scholar]
  30. Wie, B.; Lu, J. Feedback Control Logic for Spacecraft Eigenaxis Rotation Under Slew Rate and Control Constrains. J. Guid. Control Dyn. 1995, 18, 1372–1379. [Google Scholar] [CrossRef]
  31. Schweitzer, R. Active Magnetic Bearings—Chances and Limitations; International Center for Magnetic Bearings, EHT: Zurich, Switzerland, 2017. [Google Scholar]
  32. Yoon, H.; Tsiostras, P. Spacecraft adaptive attitude and power tracking with variable speed control moment giroscopes. J. Guid. Control Dyn. 2002, 25, 1081–1090. [Google Scholar] [CrossRef]
  33. Sasaki, T.; Shimomura, T. Convex optimization for power tracking of double-gimbal variable-speed control moment gyroscopes. J. Spacecr. Rocket. 2018, 55, 541–551. [Google Scholar] [CrossRef]
  34. Sasaki, T.; Shimomura, T. Fault-tolerant arhitecture of two parallel double-gimbal variable-speed control moment gyros. In Proceedings of the AIAA Guidance, Navigation andControl Conference, San Diego, CA, USA, 4–8 January 2016. [Google Scholar]
Figure 1. Inertial, local orbital, and satellite-tied frames.
Figure 1. Inertial, local orbital, and satellite-tied frames.
Micromachines 15 01159 g001
Figure 2. Frames, rotation matrix, and quaternion.
Figure 2. Frames, rotation matrix, and quaternion.
Micromachines 15 01159 g002
Figure 3. Block diagram of the satellite’s attitude calculation subsystem.
Figure 3. Block diagram of the satellite’s attitude calculation subsystem.
Micromachines 15 01159 g003
Figure 4. Base-, rotors-, and gimbal-tied frames and angular rates: (a) DGSMCMG’s structure; (b) gimbals’ angles and angular rates; (c) the graph of the interdependences between rotors, gimbals and base.
Figure 4. Base-, rotors-, and gimbal-tied frames and angular rates: (a) DGSMCMG’s structure; (b) gimbals’ angles and angular rates; (c) the graph of the interdependences between rotors, gimbals and base.
Micromachines 15 01159 g004
Figure 5. The subsystems of the j-th DGMSCMG (a) coordinates’ control subsystem; (b) gyroscopic rotors’ angular displacements control sub-system; (c) gyroscopic gimbals’ angular rates control servo-system and of the actuator with N = 2 DGMSCMGs in parallel configuration (d).
Figure 5. The subsystems of the j-th DGMSCMG (a) coordinates’ control subsystem; (b) gyroscopic rotors’ angular displacements control sub-system; (c) gyroscopic gimbals’ angular rates control servo-system and of the actuator with N = 2 DGMSCMGs in parallel configuration (d).
Micromachines 15 01159 g005
Figure 6. The frames of the DGMSCMG, the rotation angles and the angular rates (a) for DGMSCMG 1 and (b) for DGMSCMG 2.
Figure 6. The frames of the DGMSCMG, the rotation angles and the angular rates (a) for DGMSCMG 1 and (b) for DGMSCMG 2.
Micromachines 15 01159 g006
Figure 7. The subsystems of the DGMSCMGs’ models for an actuator with N = 2 DGMSCMGs in orthogonal configuration: for DGMSCMG 1 in (ad) [(a)—coordinates’ control subsystem; (b)—gyroscopic rotors’ angular displacements control sub-system; (c)—gyroscopic gimbals’ angular rates control servo-system; (d)—the subsystem modeling the interaction between DGSMCMG and the satellite] and for DGMSCMG 2 in (eh) [(e)—coordinates’ control subsystem; (f)—gyroscopic rotors’ angular displacements control sub-system; (g)—gyroscopic gimbals’ angular rates control servo-system; (h)—the subsystem modeling the interaction between DGSMCMG and the satellite].
Figure 7. The subsystems of the DGMSCMGs’ models for an actuator with N = 2 DGMSCMGs in orthogonal configuration: for DGMSCMG 1 in (ad) [(a)—coordinates’ control subsystem; (b)—gyroscopic rotors’ angular displacements control sub-system; (c)—gyroscopic gimbals’ angular rates control servo-system; (d)—the subsystem modeling the interaction between DGSMCMG and the satellite] and for DGMSCMG 2 in (eh) [(e)—coordinates’ control subsystem; (f)—gyroscopic rotors’ angular displacements control sub-system; (g)—gyroscopic gimbals’ angular rates control servo-system; (h)—the subsystem modeling the interaction between DGSMCMG and the satellite].
Micromachines 15 01159 g007
Figure 8. Subsystems of the DGMSCMG 3 used as sensor ) [(a)—coordinates’ control subsystem; (b)—gyroscopic rotors’ angular displacements control sub-system; (c)—gyroscopic gimbals’ angular rates control servo-system; (d)—the subsystem modeling the interaction between DGSMCMG and the satellite].
Figure 8. Subsystems of the DGMSCMG 3 used as sensor ) [(a)—coordinates’ control subsystem; (b)—gyroscopic rotors’ angular displacements control sub-system; (c)—gyroscopic gimbals’ angular rates control servo-system; (d)—the subsystem modeling the interaction between DGSMCMG and the satellite].
Micromachines 15 01159 g008
Figure 9. Automatic control system for satellite’s attitude control with PD controller, actuator with N = 2 DGMSCMGs, and sensor (DGMSCMG 3).
Figure 9. Automatic control system for satellite’s attitude control with PD controller, actuator with N = 2 DGMSCMGs, and sensor (DGMSCMG 3).
Micromachines 15 01159 g009
Figure 10. Dynamic characteristics of the system in Figure 9 for an actuator with N = 2 DGMSCMGs in parallel configuration: (a) the components of the satellite’s relative angular rate vectors, of the error relative angular rate vectors and of the absolute angular rate vectors; (b) satellite’s attitude with respect to the local orbital frame; (c) the Euler angles; (d) components of the torque vectors; (e) angular rates applied to the gimbals of the actuator-type DGMSCMG; (f) dynamic characteristics of the models of both actuator’s DGMSCMGs; (g) dynamic characteristics of the DGMSCMG 3 gyro sensor; (h) dynamic characteristics of the DGMSCMG 3 gyro sensor.
Figure 10. Dynamic characteristics of the system in Figure 9 for an actuator with N = 2 DGMSCMGs in parallel configuration: (a) the components of the satellite’s relative angular rate vectors, of the error relative angular rate vectors and of the absolute angular rate vectors; (b) satellite’s attitude with respect to the local orbital frame; (c) the Euler angles; (d) components of the torque vectors; (e) angular rates applied to the gimbals of the actuator-type DGMSCMG; (f) dynamic characteristics of the models of both actuator’s DGMSCMGs; (g) dynamic characteristics of the DGMSCMG 3 gyro sensor; (h) dynamic characteristics of the DGMSCMG 3 gyro sensor.
Micromachines 15 01159 g010aMicromachines 15 01159 g010bMicromachines 15 01159 g010cMicromachines 15 01159 g010d
Figure 11. Dynamic characteristics of the system in Figure 9 for an actuator with N = 2 DGMSCMGs in orthogonal configuration: (a) the components of the satellite’s relative angular rate vectors, of the error relative angular rate vectors and of the absolute angular rate vectors; (b) satellite’s attitude with respect to the local orbital frame; (c) the Euler angles; (d) components of the torque vectors; (e) angular rates applied to the gimbals of the actuator-type DGMSCMG; (f) dynamic characteristics of the models of both actuator’s DGMSCMGs; (g) dynamic characteristics of the DGMSCMG 3 gyro sensor; (h) dynamic characteristics of the DGMSCMG 3 gyro sensor.
Figure 11. Dynamic characteristics of the system in Figure 9 for an actuator with N = 2 DGMSCMGs in orthogonal configuration: (a) the components of the satellite’s relative angular rate vectors, of the error relative angular rate vectors and of the absolute angular rate vectors; (b) satellite’s attitude with respect to the local orbital frame; (c) the Euler angles; (d) components of the torque vectors; (e) angular rates applied to the gimbals of the actuator-type DGMSCMG; (f) dynamic characteristics of the models of both actuator’s DGMSCMGs; (g) dynamic characteristics of the DGMSCMG 3 gyro sensor; (h) dynamic characteristics of the DGMSCMG 3 gyro sensor.
Micromachines 15 01159 g011aMicromachines 15 01159 g011bMicromachines 15 01159 g011cMicromachines 15 01159 g011d
Table 1. Parameters of the DGMSCMG.
Table 1. Parameters of the DGMSCMG.
ParameterValueParameterValueParameterValueParameterValue
m [kg]2.8Jiz [kg·m2] 2 × 10 2 Jry [kg·m2] 45 × 10 3 kxr [N/A]0.21
lm [m] 4.1 × 10 2 Jey [kg·m2] 15 × 10 2 Jrz [kg·m2] 65 × 10 3 kyr [N/A]0.21
ls [m] 6.5 × 10 2 khx [N/m]0.8Jix [kg·m2] 2 × 10 2 kxi [Nm/A] 3 × 10 3
Jrx [kg·m2] 45 × 10 3 khy [N/m]0.8Jiy [kg·m2] 2 × 10 2 kye [Nm/A] 2 × 10 3
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

Lungu, R.; Tudosie, A.-N.; Lungu, M.-A.; Crăciunoiu, N.-C. Actuators with Two Double Gimbal Magnetically Suspended Control Moment Gyros for the Attitude Control of the Satellites. Micromachines 2024, 15, 1159. https://doi.org/10.3390/mi15091159

AMA Style

Lungu R, Tudosie A-N, Lungu M-A, Crăciunoiu N-C. Actuators with Two Double Gimbal Magnetically Suspended Control Moment Gyros for the Attitude Control of the Satellites. Micromachines. 2024; 15(9):1159. https://doi.org/10.3390/mi15091159

Chicago/Turabian Style

Lungu, Romulus, Alexandru-Nicolae Tudosie, Mihai-Aureliu Lungu, and Nicoleta-Claudia Crăciunoiu. 2024. "Actuators with Two Double Gimbal Magnetically Suspended Control Moment Gyros for the Attitude Control of the Satellites" Micromachines 15, no. 9: 1159. https://doi.org/10.3390/mi15091159

APA Style

Lungu, R., Tudosie, A. -N., Lungu, M. -A., & Crăciunoiu, N. -C. (2024). Actuators with Two Double Gimbal Magnetically Suspended Control Moment Gyros for the Attitude Control of the Satellites. Micromachines, 15(9), 1159. https://doi.org/10.3390/mi15091159

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