Advances in Engineering Research, volume 208

Proceedings of the International Conference on Innovation in Science and

Technology (ICIST 2020)

Preliminary Modelling of a CubeSat Attitude Control

System Using a Reaction Wheel

Linda Nur Afifa1,2 Tri Kuntoro Priyambodo*, Andi Dharmawan

1
Department of Information Technology Department of Computer Science and Electronics
Darma Persada University Gadjah Mada University
Jakarta, Indonesia Yogyakarta, Indonesia
2
Department of Computer Science and Electronics *[email protected], [email protected]
Gadjah Mada University
Yogyakarta, Indonesia
[email protected]

Abstract—Stabilization of satellite attitude is the most closely related to the motor's maximum speed and is usually
important thing that should always be maintained during orbit. limited by the reaction wheel's maximum voltage [8].
A model of satellite stillness controllers fairly stable against
interference and simple has been developed. This is evidenced by Due to these limitations, reaction wheels cannot perform
testing models to control motion on the Z-axis. A Proportional optimally, decreasing work performance [9]. Internal and
Integral (PI) controller has been selected as the actuator external interference, such as gravity torque, radiation torque,
controller. By tuning parameters P=423.6 and I=693.7, producing and environmental torque, may further limit the control
relatively good system characteristics, the maximum overshoot system's accuracy [10].
(MP) decreases by 3%, and the rise time improves (tr=3.1s), but
the settling time increases. By making a few modifications, this Following CubeSat's mission, modelling is required to
model can be used as a basis for transferring motion control to develop an appropriate ACS model and ensure its availability.
other axes—i.e., the X-axis (roll) and Y-axis (pitch). This article provides such modeling. According to the specified
reference, ACS controls satellite motion on three axes, i.e., X,
Keywords—CubeSat, ACS modeling, reaction wheel, PI Y, and Z axes. The control process is carried out sequentially,
controller depending on the sequence of satellite rotation. In this study,
the satellite rotation was modeled in the sequence of 323 or
I. INTRODUCTION ZYZ, i.e., rotated on the Z-axis, then the Y-axis, then to the Z-
axis. It is therefore modeled rotation control on the Z-axis.
A CubeSat is a satellite with a standard size of 10 x 10 x 10
cm3 and a maximum weight of 1.33 kg, categorized into The satellite attitude control modeling described in this
nanosatellite [1-3]. After launch, such a satellite must remain in paper begins by explaining how inertia and body frame are
a designated orbit and orientation. For this, an Attitude Control used. Then proceed with the equation's derivation for the
System (ACS) is necessary. satellite body's motion, actuators, and controllers. Furthermore,
the model equation is simulated using software to see its
ACS has two functions, first, to produce external torque to characteristics. By looking at the model's system
direct satellites, and second, to exchange momentum. The most characteristics, it can be concluded that the model is good or
widely used momentum exchange device is the reaction wheel. not, can be implemented or not.
Although this device has a fairly high accuracy level, it cannot
be very easy because it requires considerable power. This paper is organized into five sections. Section 1, the
introduction, explains the research background, Section
A reaction wheel is an actuator that works based on the describes related work, including previous research related to
principle of angular momentum. It employs a DC motor ACS modeling. Section 3, ACS modeling, contains a
mounted on a central flywheel, located at all three satellite axes derivation of mathematical equations for actuators, satellite
[4]. This rotating motor produces the torque required by the bodies, and control systems. Section 4 presents the simulation
satellite, and a change in angular velocity affects the motion and discussion, while Section 5 offers several conclusions.
and orientation of the satellite.
However, reaction wheels have a notable weakness, i.e., II. RELATED WORK
saturation [5], which occurs when a change in momentum The reaction wheel is one of the main components of ACS,
exceeds the storage capacity of the satellite [6,7]. Saturation is which is the source of disturbance [11], usually caused by

Copyright © 2021 The Authors. Published by Atlantis Press International B.V.

This is an open access article distributed under the CC BY-NC 4.0 license -http://creativecommons.org/licenses/by-nc/4.0/. 90
 Advances in Engineering Research, volume 208

flywheel imbalance, ripple torque, and friction to interfere with (3)

the accuracy and stability satellite [12].
Inertia frameworks, dynamics equations, and body frames
Reaction wheels are the most widely used actuators; it has are used as the basis for modeling the plant. The plant consists
highly accurate and produces precise torque, making them of the actuators and processes to be controlled. System
suitable for satellites that require high-speed maneuvering [13- performance is measured using the transfer function.
15]. Stabilization is an important issue that should be
considered at the time of ACS design [16].
In existing research [17,18], ACS has been modeled by
modeling the satellite body and the actuators used, i.e., reaction
wheels, and using the results to analyze system stability. The
control system often used to handle stabilization issues is PID
[19-21].
PID control systems are widely used to control the motion
of a system, such as vertical motion UAV [22], for VTOL
stabilization [23], to maintain quadrotor flight [24], namely by
tuning coefficients P, I, and D and can be optimized in value by
other methods such as Ant Colony, PSO, and others.
Fig. 2. Body framework.
III. ACS MODELING
In this section, the physical motion modeling of the satellite A. Satellite Body Motion Modeling
body, actuator, and controller will be described within the Satellite dynamics are defined by Equation 4.
context of ACS modeling. For modeling, several system
dynamics are required: the inertia framework, the motion (4)
frame, and the body frame. The inertia framework is labeled
. The point of origin is the center of the earth, Where is satellite inertia, is satellite dumping, and
is motor torque. By changing Equation 1 using the Laplace
the axis faces the vernal equinox, the axis indicates transform, output and input torque were compared, as
north, direction, and rotates using the right-hand rule, as expressed in Equation 5 and Equation 6.
described in figure 1.
(5)

(6)

B. Actuator Modeling
The actuator used to drive the satellite is the reaction wheel,
a DC motor [14,17,18]. DC motors are actuators that produce
angular acceleration on flywheels, thereby producing torque
and angular momentum that can change the satellite's
orientation. Change in momentum is generated by the reaction
wheel to move the satellite in the opposite direction.
Fig. 1. Inertia framework. DC motor modeling is grouped into three elements: the
electrical system, the mechanical system, and the rotational
The body frame describes the position of the object's body system, as expressed in Equation 7.
when rotating in three axes, namely the x-axis (roll), the Y-axis
(pitch), and the Z-axis (yaw), as described in Figure 2. The (7)
torque is generated on each axis, expressed by equation 1-3.
ACS controls satellite torque on roll, pitch and yaw axes. In is load torque, is the internal disturbance torque
this simulation, the torque to be controlled is the torque on the motor; in this initial study, the disturbance was ignored. By
Z-axis. using the Laplace transform, Equation 8 was obtained.

(1)
(8)
(2)

91
 Advances in Engineering Research, volume 208

C. Disturbance The analysis is subsequently used to examine the control

Frequently interfering with low-orbit satellites (LEO) are system's characteristics, using P and I control parameters of
1522 and 1880. The characteristics of the system are graphed
gravity gradient ( ), magnetic torque ( ), aerodynamic in Figure 4. From the figure, it can be seen that rise time starts
torque ( ), and solar radiation ( ) [1,25,26]. to increase at and that maximum overshoot
In this study, total disturbance torque was used as a (Mp=16.3%) exceeds the input value of 1. The system
reference; satellites should produce torque that is at least twice achieves a steady-state (settling time) at . The
as much as disturbance torque. In this study, the total amount characteristic values of this system are fully presented in Table
used was 6.850 10-7[Nm], the gravity gradient, magnetic 2.
torque, aerodynamic torque, and solar radiation.
The next experiment is used to tune the PI control
parameters. The P and I values used in this second experiment
IV. SIMULATION AND DISCUSSION were 423.6 and 693.7. Following these parameters, the system
Based on the motion modeling of satellite bodies, actuators, characteristics in Figure 5 are generated.
controllers, and interference, ACS systems were diagrammed
in a transfer function. To identify the ACS characteristics, the Figure 5 shows a rise time of , a maximum
parameter values used for each variable of the simulation are overshoot of , and a settling time of
given in Table 1. . The complete characteristic values of this
simulation can be seen in Table 3.
TABLE I. VARIABLE OF SIMULATION
Figure 4 and Figure 5 show the behavior of the system, as
Variable Simulasi Magnitude resulting from physical modeling. The output chart first rises to
Amplifier Gain 10
=1.56; after tuning, the control system parameters increases
Resistance 10 Ohm
to =3.1. The value is affected by satellite dumping.
Inductance 0.5 H
Inertia Motor 0.01 Kg.sqm
Inertia Satelit 2.5 Kg.sqm
Damping 1.17 Nms
DC motor constant 0.01 Nm/s

The results of this modeling of satellite motion, actuators,

and control systems were analyzed by looking at the
controller's characteristics and the overall output of the model
simulation. Figure 3 illustrates the satellite models' output (blue
lines) and actuators (red lines) by providing step input. In the
image, the settling time is reached at = 3.31, and the output is
close to the given input.

Fig. 4. System characteristics, P=1522, I=1880.

TABLE II. SYSTEM CHARACTERISTICS VALUE

Performance and Robustness
Tuned Block
Rise time 1.56 seconds 16.4 seconds
Settling time 7.81 seconds 29 seconds
Overshoot 16.3% 0%
Peak 1.16% 1
Gain margin 5.45 dB @4.37 32.2 dB @4.28
Fig. 3. System output. rad/s rad/s
Phase margin 60 deg @0.936 79.7 deg @ 0.108
Based on the simulation results, this system is stable, as rad/s rad/s
Closed-loop stability Stable Stable
indicated by the equation characteristics of the system being
negative (namely -0.5877 + 4.4278i, -0.5877 -4.4278i, -
0.3399 + 0.0000i, -0.1527 + 0.0000i).

92
 Advances in Engineering Research, volume 208

ACKNOWLEDGMENT
This research is supported by doctoral dissertation research
form The Directorate General of Higher Education, The
Indonesian Ministry of Education and Culture through contract
number 2086/UN1/DITLIT/DIT-LIT/PT/2020.

REFERENCES

[1] E. Oland and R. Schlanbusch, “Reaction Wheel Design for CubeSat,” in

4th International Conference on Recent Advances in Space
Technologies, Istanbul, 2009, pp. 778–783, 2009.
[2] The Tauri Group, “State of the Satellite Industry Report,” no.
September, p. 33, 2014. [Online]. Retrieved from:
http://www.sia.org/wp-content/uploads/2014/09/SSIR-September-2014-
Fig. 5. System characteristics P=423.6, I=693.7. Update.pdf
[3] T.K. Priyambodo, A.E. Putra, M. Asvial, and R.E. Poetro, “IiNUSAT-1 :
The 1 st Indonesian Inter-University Nano-Satellite for Research and
TABLE III. SYSTEM CHARACTERISTICS VALUE Education,” In 2014 IEEE International Conference on Aerospace
Performance and Robustness Electronics and Remote Sensing Technology, pp. 114–120, 2014.
Tuned Block [4] R. Votel and D. Sinclair, “Comparison of control moment gyros and
Rise time 3.1 seconds 16.4 seconds reaction wheels for small Earth-orbiting satellites,” in Proceedings of the
Settling time 10.5 seconds 29 seconds 26th Annual AIAA/USU Conference on Small Satellites, pp. 1–7, 2012.
Overshoot 10.5 % 0% [5] J. Chen, Z. Li, M.G. Gan, and G. Zhang, “Adaptive robust control for
Peak 1.11% 1 DC motors with input saturation,” IET Control Theory Appl., vol. 5, no.
Gain margin 16.2 dB @4.32 32.2 dB @4.28 16, pp. 1895–1905, Nov. 2011.
rad/s rad/s [6] J.L. Crassidis, Fundamentals of Spacecraft Attitude Determination and
Phase margin 60 deg @0.448 79.7 deg @ 0.108 Control. New York Heidelberg,Heidelberg Dordrecht London: Springer,
rad/s rad/s 2014.
Closed-loop stability Stable Stable [7] R. Fros, “How can I understand the saturation reaction wheel used in
Attitude Control,” 2014.
Furthermore, the maximum overshoot (MP) also decreased [8] K. Kong, H.C. Kniep, and M. Tomizuka, “Output Saturation in Electric
Motor Systems: Identification and Controller Design,” J. Dyn. Syst.
by 3%. is an important and necessary value when looking Meas. Control, vol. 132, no. 5, pp. 051002-1-051002–8, 2010.
to model a physical system. A good system will not produce [9] X. Cao and B. Wu, “Indirect adaptive control for attitude tracking of
much more output than a given input. In the control system, spacecraft with unknown reaction wheel friction,” Aerospace science
and technology, vol. 47, pp. 493–500, 2015.
this is represented by the maximum overshoot . In this [10] P. Guan, X. Liu, W. Zhang, and L. Xue, “The Direct Adaptive Fuzzy
simulation, the value was 10.5%. Robust Control for Satellite,” in Proceedings of 10th world Congress on
Intelligent Control and Automation, pp. 36–41, 2012.
Another important parameter to consider is settling time [11] S. Taniwaki and Y. Ohkami, “Experimental and Numerical Analysis of
Reaction Wheel Disturbances,” JSME Int. J. Ser. C, vol. 46, no. 2, pp.
( ), which indicates how quickly a system achieves its final 519–526, 2003.
value (usually 2% or 5% of the set value). In the second [12] M.F. Mehrjardi, H. Sanusi, M. Alauddin, and M. Ali, “Developing a
simulation, increased about 3s. This is an important record Proposed Satellite Reaction Wheel Model with Current Mode Control,”
at the time of modeling. In 2015 International Conference on Space Science and Communication
(IconSpace), pp. 0–3, 2015.
[13] Y. Yang, “Spacecraft Attitude and Reaction Wheel Desaturation
V. CONCLUSION Combined Control Method,” IEEE Transactions on Aerospace and
Electronic Systems, vol. 53, no. 1, 2017.
The success of modeling can be ascertained by looking at
[14] J. Li, M. Post, T. Wright, and R. Lee, “Design of Attitude Control
the system characteristics, i.e., rise time , settling time Systems for CubeSat-Class Nanosatellite,” J. Control Sci. Eng., vol.
( ), and maximum overshoot ( . Based on the simulation 2013, no. April, pp. 1–15, 2013.
results, the modeling achieved a better maximum overshoot [15] Y.T. Xing, K.S. Low, S. Member, and M.D. Pham, “Distributed Model
Predictive Control of Satellite Attitude Using Hybrid Reaction Wheels
( ) after tuning the controller parameters but required more and Magnetic Actuators,” In 2012 IEEE Symposium on Industrial
Electronics and Applications, pp. 230–235, 2012.
settling time ( ). It is necessary to review the physical
[16] A. Aydogan and O. Hasturk, “Adaptive LQR Stabilization Control of
modeling undertaken and the size of the variables given to the Reaction Wheel for Satellite Systems,” In 2016 14th International
system. Conference on Control, Automation, Robotics and Vision (ICARCV),
vol. 2016, no. November, pp. 13–15, 2016.

93
 Advances in Engineering Research, volume 208

[17] J. Carlos, “Attitude Control Model for CubeSats,” Iaa-Br-07-01, no. vertical moving control of UAV system,” in AIP Conference
February, pp. 1–22, 2016. Proceedings, vol. 1755, no. July, p. 170011, 2016.
[18] J. Li, M. Post, T. Wright, and R. Lee, “Design of Attitude Control [23] T.K. Priyambodo, A.E. Putra, and A. Dharmawan, “Optimizing control
Systems for CubeSat-Class Nanosatellite,” J. Control Sci. Eng., vol. based on ant colony logic for Quadrotor stabilization,” in 2015 IEEE
2013, pp. 1–15, 2013. International Conference on Aerospace Electronics and Remote Sensing
[19] P. Suganthi, S. Nagapavithra, and S. Umamaheswari, “Modeling and Technology (ICARES), pp. 1–4, Dec. 2015.
Simulation of Close Loop Speed control for BLDC Motor,” in Proc. [24] A. Dharmawan, A. Ashari, and A.E. Putra, “Quadrotor flight stability
IEEE Conference on Emerging Devices and Smart Systems (ICEDSS system with Routh stability and Lyapunov analysis,” in AIP Conference
2017), no. March, pp. 83–126, 2017. Proceedings, vol. 1755, no. July, p. 170007, 2016.
[20] V. Carrara and H.K. Kuga, “Torque and Speed Control Loops Of A [25] D. Hagen, Modelling and Controller Design considering Actuator
Reaction Wheel,” 11th Int. Conf. Vib. Probl., no. September, pp. 9–12, Dynamics. Norway: Narvik University College, 2006.
2013. [26] V. Izadi, M. Abedi, and H. Bolandi, “Verification of reaction wheel
[21] F. Mitin and A. Krivushov, “Application of Optimal Control Algorithm functional model in HIL test-bed,” in 2016 4th International Conference
for DC Motor,” in 29th DAAM International Symposium on Intelligent on Control, Instrumentation, and Automation (ICCIA), no. January, pp.
Manufacturing and Automation, pp. 0762–0766, 2018. 155–160, Jan. 2016.
[22] T.K. Priyambodo, A. Dharmawan, O.A. Dhewa, and N.A.S. Putro,
“Optimizing control based on fine tune PID using ant colony logic for

94