elevation, pitch and travel axis stabilization of 3dof
TRANSCRIPT
International Journal of Electrical and Computer Engineering (IJECE)
Vol. 10, No. 2, April 2020, pp. 1868~1884
ISSN: 2088-8708, DOI: 10.11591/ijece.v10i2.pp1868-1884 1868
Journal homepage: http://ijece.iaescore.com/index.php/IJECE
Elevation, pitch and travel axis stabilization of 3DOF helicopter
with hybrid control system by GA-LQR based PID controller
Ibrahim K. Mohammed, Abdulla I. Abdulla Department of Systems and Control Engineering, College of Electronics Engineering, Ninevah University, Iraq
Article Info ABSTRACT
Article history:
Received Apr 29, 2019
Revised Aug 28, 2019
Accepted Oct 9, 2019
This research work introduces an efficient hybrid control methodology
through combining the traditional proportional-integral-derivative (PID)
controller and linear quadratic regulator (LQR) optimal controlher.
The proposed hybrid control approach is adopted to design three degree of
freedom (3DOF) stabilizing system for helicopter. The gain parameters of
the classic PID controller are determined using the elements of the LQR
feedback gain matrix. The dynamic behaviour of the LQR based PID
controller, is modeled in state space form to enable utlizing state feedback
controller technique. The performance of the proposed LQR based LQR
controller is improved by using Genetic Algorithm optimization method
which are adopted to obtain optimum values for LQR controller gain
parameters. The LQR-PID hybrid controller is simulated using Matlab
environment and its performance is evaluated based on rise time, settling
time, overshoot and steady state error parameters to validate the proposed
3DOF helicopter balancing system. Based on GA tuning approach,
the simulation results suggest that the hybrid LQR-PID controller can be
effectively employed to stabilize the 3DOF helicopter system.
Keywords:
3DOF helicopter system
Genetic algorithm
Hybrid control
LQR controller
PID controller
Copyright Β© 2020 Institute of Advanced Engineering and Science.
All rights reserved.
Corresponding Author:
Ibrahim K. Mohammed,
Department of Systems and Control Engineering,
College of Electronics Engineeing,
Ninevah University, Mosul, Iraq
Email: [email protected]
1. INTRODUCTION
The A 3DOF helicopter model is a multi-input, multi-output system (MIMO) with strongly
nonlinear characteristics, used as experimental platform for verification of various control algorithms.
The helicopter is a complex mechanical scheme and has unstable dynamics which make using traditional
control methods to stabilize it is a big challenge. In control procedure of 3DOF helicopter system
the nonlinear dynamics is linearised about preselected operating conditions. Based on linearised system
model a classical proportional integral derivative PID controller technique is widely used due to its simplicity
and easy to realize [1]. However, using conventional PID technique is insufficient to guide the helicopter
through the demand trajectory as the regulation of the controller parameters is hard task even for
an experienced control engineers [1, 2]. In another control method, an optimal tracking approach uilizing
LQR controller to stabilize the helicopter system was proposed by [3]. In the proposed system,
a Genetic Algorithm (GA) is employed to optimize a classic LQR controller for position control of DC
motor. By using LQR controller a significant improvement in guidance of helicopter system is carried out.
Further investigations in control strategies for helicopter system are still the goal of many researchers. A new
control approach based on hybrid controller technique was proposed by many researchers to control
the 3DOF helicopter system. Fan and Joo [4] proposed a hybrid control method using fuzzy and PID
controller to balance the 3DOF helicopter system at desired equilibrium position. In this study, an efficient
Int J Elec & Comp Eng ISSN: 2088-8708
Elevation, pitch and travel axis stabilization of 3DOF helicopter with⦠(Ibrahim K. Mohammed)
1869
control method has been adopted to propose a hybrid control strategy using both conventional and intelligent
control algorithms. This hybrid control approach uses PID controller as a conventional control, while
the liner quadratic regulator (LQR) controller based on GA used as an intelligent control. The LQR optimal
controller is recommended to implement the intelligent system due to its high precision in movement
applications [5]. In this state feedback controller technique a tradeoff between the plant characteristics and
the input effort of the system can be achieved if the controller gain parameters are tuned propperly [6-8].
The LQR controller design can be improved by employing soft computing and optimization techniques,
which are adopted to achieve tuning proecss of the controller gain parameters, such as Differential Evolution
(DE) [9], Particle Swarm Inspired Evolutionary Algorithm (PS-EA) [10], Particle Swarm Optimisation
(PSO) [11], Genetic Algorithm (GA) [12], memetic algorithm (MA) [13], Imperialist Competitive Algorithm
(ICA) [14], Ant Colony Optimization (ACO) [15], Artificial Bee Colony (ABC) [16], and Artificial Immune
Systems (AIS) [17]. In order to validate the proposed helicopter model, the hybrid control system is
simulated by using the Matlab environment, and its balancing performance is then evaluated based on
stability criteria parameters which include rise and settling time, overshoot and control input. The purpose of
the proposed hybrid GA-LQR based PID controller is to design an efficient control system utilized to
stabilize a prototype 3DOF helicopter system at desired roll and pitch positions as well as angular travel
speed. During the current decade there has been a considerable interest by many control researchers in using
of hybrid control approach in various engineering applications. Vendoh and Jovitha [18] presented a hybrid
control system using LQR based PID controller for stabilization and trajectory tracking of magnetic levitation
system. Arbab et al. [19] proposed a hybrid control strategy using fuzzy based LQR controller for 3-DOF
helicopter system. Choudhary [1] design another hybrid control scheme using LQR based PID controller for
3-DOF helicopter system, however, the introduced hybrid control systems are not efficient as the parameters
of the LQR feedback gain matrix are obtained by applying trial and error procedure. The paper presents
a simple method for the approximation of PID controller gain parameters from feedback gain matrix of
the controller. The LQR controller gain is obtained from the weighting matrices Q and R which their
elements values are tuned effectively by the GA optimization method. The rest of the paper is organised as follows. Section 2 presents configuration and dynamics
modeling of the helicopter system. In section 3, controller techniques of helicopter hybrid control system are
introduced. Tuning method is described in section 4. Section 5 introduces helicopter control system design is
presented. In section 6 simulation results of the hybrid control system are introduced and followed by
conclusions and prospective research in section 7.
2. STRUCTURE AND DYNAMICS MODELING OF HELICOPTER
2.1. Helicopter stucture
The conceptual platform of 3-DOF helicopter system is shown in Figure 1. It consists of an arm
mounted on a base. The main body of the helicopter composed of propellers driven by two motors mounted
are the either ends of an short balance bar. The whole helicopter body is fixed on one end of the arm and
a balance block installed at the other end. The balance arm can rotate about travel axis as well as slope on an
elevation axis. The body of helicopter is free to roll about the pitch axis. The system is provided by encoders
mounted on these axis used to measure the travel motion of the arm and its elevation and pitch angle.
The propellers with motors can generate an elevation mechanical force proportional to the voltage power
supplied to the motors. This force can cause the helicopter body to lift off the ground. It is worth considering
that the purpose of using balance block is to reduce the voltage power supplied to the propellers motors.
2.2. Helicopter dynamics and mathematical modeling
In this study, the nonlinear dynamics of 3DOF helicopter system is modeled mathematically based
on developing the model of the system behavior for elevation, pitch and travel axis. The definition of
the symbols and nomecalture of the proposed helicopter system is included in Table 1.
2.2.1. Elevation axis model
The free body diagram of 3DOF helicopter system based on elevation axis is shown in Figure 2.
The movement of the elevation axis is governed by the following differential equations:
π½ππΜ = ππ + πππΉπcos(π) β ππ€,π β ππ,π (1)
where ππ is the centrifugal torque, which is function of (π), ππ,π is friction moment exerted on the helicopter
elevation axis, which is generated from combining joint friction and air resistance, and ππ€,π is the effective
gravitational torque due to the mass differential π = πβπ about the elevation axis.
ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 10, No. 2, April 2020 : 1868 - 1884
1870
Figure 1. 3DOF Helicopter System [20]
Figure 2. Schematic diagram of elevation axis model for 3DOF helicopter system
ππ€,π = πππβπ (2)
If the elevation angle π = 0 and the rotation angle of pitch axis π = 0, then the centrifugal torque (ππ) will
be zero (1) based on Eulerβs second law becomes as follows:
π½ππΜ = πππΉπ β ππ€,π + ππ,π (3)
π½ππΜ = ππ(πΉπ + πΉπ) β ππ€,π + ππ,π (4)
πΉπ = πΎππππ = π, π (5)
π½ππΜ = πΎπππ(ππ + ππ) β ππ€,π + ππ,π (6)
π½ππΜ = πΎπππππ β ππ€,π + ππ,π (7)
2.2.2. Pitch axis model
Consider the pitch schematic diagram of the system in Figure 3. It can be seen from the figure that
the main torque acting on the system pitch axis is produced from the thrust force generated by the properller
motors. When π β 0, the gravitational force will also generate a torque ππ€,π acts on the helicopter pitch axis.
ππ€ ,π + ππ ,π
π
ππ
ππ
Int J Elec & Comp Eng ISSN: 2088-8708
Elevation, pitch and travel axis stabilization of 3DOF helicopter with⦠(Ibrahim K. Mohammed)
1871
Figure 3. Schematic diagram of pitch axis model for 3DOF helicopter scheme
Table 1. Nomenclature of the helicopter system Nomenclature Definition
π, π, π
ππ
ππ
ππ
πβ
πΉπ
Rotation angle of elevation, pitch and travel axis respectively (degree)
Distance between the pitch axis and the helicopter center of gravity (m)
Distance from axle point to either propeller motor (m)
Distance from axle point to the balance block (m)
Distance between the pitch axis and the center of propellers (m)
Lift force of the properller motors (N)
πΉπ , πΉπ
ππ , ππ
Forces of fornt and back properller motors (N)
Votage applied to the front and back properller motor respectively (v)
ππ
ππ
π
ππΉ,π , ππΉ,π,ππΉ,π
ππ
ππ€,π ,ππ€,π
ππ , ππ, ππ
πΎπ
Mass of balance block (kg)
Mass of properller motors (kg)
Mass differential about the elevation axis (kg)
Torque component generated from combining the joint friction and air resistance for elevation, pitch and
rate axis respectively (N.m)
Centrifugal torque about the elevation axis (N.m)
Effective gravitational torque due to gravity about the elevation axis and pitch axis respectively (N.m)
Inertia moment of the system about elevation, pitch and travel axis respectively (kg.m2)
Force constant of the motor-propeller combination
The dynamics of the pitch axis can be modeled mathematically as follows:
π½ποΏ½ΜοΏ½ = πΉππβ β πΉππβ β ππ€,π β ππ,π (8)
Where ππ,πis friction moment exerted on the pitch axis.
ππ€,π = πβπππ sin(π) cos(π) (9)
Based on the assumption that the pitch angle π = 0, ππ€,π = 0 , then (8) becomes as follows:
π½ποΏ½ΜοΏ½ = πβ(πΉπ β πΉπ) β ππ,π (10)
π½ποΏ½ΜοΏ½ = πΎπππ(ππ β ππ) β ππ,π (11)
π½ποΏ½ΜοΏ½ = πΎπππππ β ππ,π (12)
2.2.3. Travel axis model The free body diagranm of the helicopter system dynamics based on travel axis is presented in
Figure 4. In this model, when π β 0, the main forces acting on the helicopter dynamics are the thrust forces
of properller motors (πΉπ, πΉπ). These forces have a component generates a torque on the travel axis. Assume
ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 10, No. 2, April 2020 : 1868 - 1884
1872
that the helicopter body has roll up by an angle π as shown in Figure 4. Then the dynamics of travel axis for
3DOF helicopter system is modeled as follows:
π½ποΏ½ΜοΏ½ = β(πΉπ + πΉπ)π ππ(π)ππ β ππ,π (13)
Figure 4. Schematic diagram of travel rate axis model for 3DOF helicopter scheme
The thrust foces of the two properller motors (πΉπ + πΉπ) are required to keep the helicopter in flight
case and is approximately π.
π½ποΏ½ΜοΏ½ = βππ ππ(π)ππ β ππ,π (14)
Where ππ,π is friction moment exerted on travel axis. As π approachs to zero, based on sinc function,
π ππ(π) = π, the above equation becomes as follows:
π½ποΏ½ΜοΏ½ = βππππ β ππ,π (15)
Based on the assumption that the coupling dynamics, gravitational torque (ππ€,π) and friction moment exerted
on elevation, pitch and travel axis are neglected, then the dynamics modeling (7), (12) and (15) for 3DOF
helicopter system can be simplified as in (16), (17) and (18) respectively [1].
πΜ =πΎπππ
π½πππ (16)
οΏ½ΜοΏ½ =πΎπππ
π½πππ (17)
οΏ½ΜοΏ½ =πππ
π½ππ (18)
2.3. System state space model
In order to design state feedback controller based on LQR technique for 3DOF helicopter system,
the dynamics model of the system should be formulated in state space form. In this study, the proposed
hybrid control algorithm is investigated for the purpose of control of pitch angle, elevation angle and travel
rate of 3DOF helicopter scheme by regulating the voltage suopplies to the front and back motors.
Let π₯(nx1) = [π₯1,π₯2, π₯3,π₯4, π₯5,π₯6, π₯7]π = [π, π, πΜ, π,Μ π, Κ, πΎ]π be the state vector of the system, the state
variables are chosen as the angles and rate and their corresponding angular velocities, and ΚΜ = π, οΏ½ΜοΏ½ = π.
The voltages supplied to the front and back propellers motors are considered the inputβs vector such that,
π’(π‘)(mx1) = [π’1,π’2]π
= [ππ, ππ]π and the elevation angle, pitch angle and travel rate are assumed
Int J Elec & Comp Eng ISSN: 2088-8708
Elevation, pitch and travel axis stabilization of 3DOF helicopter with⦠(Ibrahim K. Mohammed)
1873
the outputβs vector such that, π¦(π‘)(px1) = [π, π, π]π. Based on (13)-(15), choosing these state variables
yields the following system state space model:
π₯1Μ = π = π₯2
π₯2Μ = πΜ = π₯3
π₯3Μ = πΜ =πΎππππ½π
(ππ + ππ)
π₯4Μ = οΏ½ΜοΏ½ =πΎπππ
π½π(ππ β ππ) (16)
π₯5Μ = οΏ½ΜοΏ½ =ππππ½π
π₯2
π₯6Μ = ΚΜ = π₯1
π₯7Μ = οΏ½ΜοΏ½ = π₯4
The general state and output matrix equations describing the dynamic behavior of the linear-time-
invariant (LTI) helicopter system in state space form are as follows:
π₯(Μπ‘) = π΄π₯(π‘) + π΅π’(π‘) (17)
π¦(π‘) = πΆπ₯(π‘) + π·π’(π‘) (18)
Where π΄(nxn) is the system matrix, π΅(nxm) is the input matrix, πΆ(pxm) is the output matrix, and π·(mxp)
is feed forward matrix, for the designed system. Based on (16), the state space representation (17) and (18)
are rewritten as in (19) and (20) [1].
[ 0 0 1 0 0 0 00 0 0 1 0 0 00 0 0 0 0 0 00 0 0 0 0 0 0
0πππ
π½π0 0 0 0 0
1 0 0 0 0 0 00 0 0 0 1 0 0]
[ πππΜοΏ½ΜοΏ½πΚπΎ]
+
[
0 00 0
πΎπππ
π½π
πΎπππ
π½ππΎπππ
π½πβ
πΎππ
π½π
0 00 00 0 ]
[ππ
ππ] (19)
[πππ] = [
1 0 0 0 0 0 00 1 0 0 0 0 00 0 0 0 1 0 0
]
[ πππΜοΏ½ΜοΏ½πΚπΎ]
+ [0 00 00 0
] [ππ
ππ] (20)
[ πΜοΏ½ΜοΏ½πΜοΏ½ΜοΏ½οΏ½ΜοΏ½ΚΜοΏ½ΜοΏ½]
=
[ 0 0 1 0 0 0 00 0 0 1 0 0 00 0 0 0 0 0 00 0 0 0 0 0 00 2.0655 0 0 0 0 01 0 0 0 0 0 00 0 0 0 1 0 0]
[ πππΜοΏ½ΜοΏ½πΚπΎ]
+
[
0 00 0
5.8197 5.819763.9498 β63.9498
0 00 00 0 ]
[ππ
ππ] (21)
In this study, for purpose of control system design, the model of the system is formulated in state
space form using the physical parameters values stated in Table 2 [1]. Based on (19) and using
the parameters values in Table 1, the state equation of the system is given by (21).
ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 10, No. 2, April 2020 : 1868 - 1884
1874
Table 2. Values of physical parameters of 3DOF helicopter system Symbol Physical Unit Numerical Values
π½π kg.m2 1.8145 π½π kg.m2 1.8145 π½π kg.m2 0.0319 W N 4.2591 ππ m 0.88 ππ m 0.35
ππ m 0.17
πΎπ N/V 12
3. CONTROLLER TECHNIQUES
In this section, the theory of the PID and LQR controller techniques used for designing the proposed
helicopter control system are introduced. The gain parameters of the PID controller are obtained
approximately from the LQR gain matrix. Control equations of elevation angle, pitch angle, and travel rate
for helicopter system using PID controller are also presented.
3.1. PID controller
A PID is the most popular controller technique that is widely used in industrial applications due to
the simplicity of its structure and can be realized easily for various control problems as the gain parameters of
the controller are relatively independent [21][22]. Basically, the controller provides control command signals
π’(π‘)based on the error π(π‘) between the demand input and the actual output of the system. The continuous
time structure of the PID controller is as follows:
π’(π‘) = πΎππ(π‘) +πΎπ β« π(π)πππ‘
0+ πΎπ
ππ(π‘)
ππ‘ (22)
Where πΎπ, πΎπand πΎπ are the proportional, integral and differential components of the controller gain.
These controller gain parameters should be tuned properly to enable the output states of the system to
efficiently follow the desired input. If ππ, ππ and ππ are the desired elevation angle, pitch angle and travel
rate of the helicopter system, it can express the form of PID controllers used to meet the desired output states
as follows [1,23]. In this study, for elevation angle, the control equation is based on the following PID control
equation:
ππ = πΎππππ + πΎππποΏ½ΜοΏ½ + πΎππ β« ππππ‘ (23a)
ππ = πΎππ(ππ β π) β πΎπππΜ + πΎππ β«(ππ β π)ππ‘ (23b)
While the pitch angle is controlled by the following PD control equation:
ππ = πΎππππ + πΎππποΏ½ΜοΏ½ (24a)
ππ = πΎππ(ππ β π) β πΎπποΏ½ΜοΏ½ (24b)
The travel rate is governed by the following PI control equation:
ππ = πΎππππ + πΎππ β« ππππ‘ (25a)
ππ = πΎππ(ππ β π) + πΎππ β«(ππ β π)ππ‘ (25b)
where ππ = ππ β π, ππ = ππ β π, ππ = ππ β π, ποΏ½ΜοΏ½ = βπΜ and ποΏ½ΜοΏ½ = βοΏ½ΜοΏ½.
3.2. LQR controller
Linear quadratic regulator is a common optimal control technique which has been widely utilized in
various manipulator systems [24]. Using LQR technique in design control system requires all the states of
the system plant to be measurable as it bases on the full state feedback concept. Therefore, using LQR
controller to stabilize the 3DOF helicopter system based on the assumption that the states of the system are
considered measurable. LQR approach includes applying the optimal control effort:
Int J Elec & Comp Eng ISSN: 2088-8708
Elevation, pitch and travel axis stabilization of 3DOF helicopter with⦠(Ibrahim K. Mohammed)
1875
π’(π‘) = βπΎπ₯(π‘) (26)
Where πΎ is the state feedback gain matrix of LQR controller, that will enable the output states of the system
to follow the trajectories of reference input, while minimizing the following the cost function:
π½ = β« (π₯π(π‘)ππ₯(π‘)β
0β π’π(π‘)π π’(π‘))ππ‘ (27)
Where Q and R are referred to as weighting state and control matrices. The controller feedback gain matrix
can be determined by using below equation:
πΎ = π β1π΅ππ (28)
Where π is (nxn) matrix determined from the solution of the following Riccati matrix equation:
π΄ππ + ππ΄ β ππ΅π β1π΅ππ + π = 0 (29)
In this application, πΎ = [k11 k12 k13 k14 k15 k16 k17
k21 k22 k23 k24 k25 k26 k27] and π’(π‘) = [
π’1
π’2]. Based on this,
the control effort of the system stated in (26) can be written as follows:
[π’1
π’2] = β [
k11 k12 k13 k14 k15 k16 k17
k21 k22 k23 k24 k25 k26 k27]
[ πππΜοΏ½ΜοΏ½πΚπΎ]
(30)
For the purpose of simplicity of control problem the weighting matrices Q and R are chosen as the diagonal
matrices so that the cost function (27) can be reformulated as below:
π½ = β« (π11π₯12 + π22π₯2
2 + π33π₯32 + π44π₯4
2 + π55π₯52 + π66π₯6
2 + π77π₯72 + π11π’1
2 + π22π’22)
β
0ππ‘ (31)
Where π11, π22, π33, π44, π55, π66 and π77 denote the weighting elements of flight angles and their
corresponding angular velocities of the proposed 3DOF helicopter system respectively, while, π11and π22 are
the weighting elements of control inputs. It is worth considering that, in this study, based on the state and
control LQR weighting matrices Q and R, the feedback gain matrix (πΎ) can be calculated by utilizing
the Matlab command βlqrβ.
The optimal control approach LQR is highly recommended for stabilizing the 3DOF helicopter
system as it basically looks for a compromising between the best control performance and minimum control
effort. Based on LQR controller, an optimum tracking performance can be investigated by a proper setting of
the feedback controller gain matrix. To achieve this, the LQR controller is optimized by using GA tuning
method which is adopted to obtain an optimum elements values for Q and R weighting matrices.
4. TUNING METHOD
In this study, GA tuning approach is adopted to optimize the LQR gain matrix used to approximate
the gain parameters of PID controller for 3DOF helicopter system. GA is a global search optimization
technique bases on the strategy of genetics and natural selection [25][26]. This optimization method is
utilized to obtain an optimum global solution for more control and manipulating problems. The procedure of
GA approach includes of three basic stepes: Selection, Crossover and Mutation. Applying these main
operations creates new individuals which could be better than their parents. Based on the requirements of
desired response, the sequence of GA optimization technique is repeated for many iterations and finally stops
at generating optimum solution elements for the application problems [27]. The sequence of the GA tuning
method is presented in Figure 5 [3]. The steps of the GA loop are defined as follows:
1. Initial set of population.
2. Selecting individuals for mating.
3. Mating the population to create progeny.
4. Mutate progeny.
ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 10, No. 2, April 2020 : 1868 - 1884
1876
5. Inserting new generated individuals into population.
6. Are the system fitness function satisfied?
7. End search process for solution.
In this study, the aim of using GA optimization method is to tune the PID controller parameters by
minimizing a selected fitness function which will be discussed in detail in the section 4. The implementation
procedure of the GA tuning method begins with the definition step of the chromosome representation. Each
chromosome is represented by nine cells which correspond to the weight matrices elements of the LQR
controller as shown in Figure 6. By this representation it can adjust the LQR elements in order to achieve
the required performance. These cells are formed by real positive numbers and characterize the individual to
be evaluated [27].
Figure 5. Process loop of GA.
Figure 6. Definition of GA chromosome
5. HELICOPTER HYBRID CONTROLS YSTEM DESIGN
In this study, a control system using LQR based PID controller is designed to stabilize 3DOF
helicopter system. Based on step input, a hybrid controller is designed for the following desired performance
parameters: rise time (π‘π) less than 10 (ms), settling time (π‘π ) less than 30 (ms), maximum overshoot
percentage, (ππ) less than 5%. The block diagram of the proposed helicopter control system based on LQR
controller is presented in Figure 7.
Figure 7. LQR controller based on GA for 3DOF helicopter system
11q 22q 33q 44q 55q 66q 77q 11r 22r
Int J Elec & Comp Eng ISSN: 2088-8708
Elevation, pitch and travel axis stabilization of 3DOF helicopter with⦠(Ibrahim K. Mohammed)
1877
The control system is analysed mathematically and then simulated using Matlab/simulink tool to
validate the proposed LQR controller. Based on the desired performance parameters which include rise and
settling time, overshoot and error steady state parameters the fitness function of the control problem is
formulated as follows:
πΉ = 0.3π. π‘π + 0.3π. π‘π + 0.2π. π + 0.2π. ππ π (32)
Where, π is closed loop transfer function of the helicopter system, π. π‘π , π‘π , π, ππ π are the rise time, settling
time, maximum over shoot and error steady state of the closed-loop control system. It is worth considering
that the control input effort is considered in the evaluation process of the proposed stabilizing helicopter
system. In this study, the design of controller is effectively optimized by using GA tuning method which is
adopted to obtain optimum elements values for LQR weighting matrices Q and R. These optimized matrices
are using to calcutate optimum controller gain matrix by using (29) and (30). However, the gain matrix is
determined by using the Matlab command βlqrβ.
5.1. PID approximation
In this subsection, the gain parameters πΎπ, πΎπ , πΎπof the PID controller are calculated approximately
from the feedback gain matrix πΎ of the LQR controller based on GA optimization technique. Analyzing (30)
yields the following [4]:
[π’1
π’2] = β [
k11 k12 k13 k14 k15 k16 k17
k11 βk12 k13 βk14 βk15 k16 βk17]
[ πππΜοΏ½ΜοΏ½πΚπΎ]
(33)
5.1.1. Elevation control using PID controller Summing the rows of (33) results the following [1]:
π’1 + π’2 = β(2π11π + 2π13πΜ + 2π16Κ) = β(2π11π + 2π13πΜ + 2π16 β« π ππ‘) (34)
The above equation can be written as
ππ = β2π11(ππ β π) β 2π13πΜ β 2π16 β«(ππ β π) ππ‘ (35)
It is obvious that (23b) and (35) have the same structure, this means that the gain parameters of pitch PID
controller can be obtained from the gain elements of LQR controller. Thus, comparing (23b) with (35), yields
the following gain relationships:
πΎππ = 2π11
πΎππ = 2π13
πΎππ = 2π16
(36)
The block diagram of closed-loop control system for 3DOF helicopter system based on PID controller is
shown in Figure 8. Taking Laplace transform for elevation axis model (13) yields the following equation:
π½ππ(π ). π 2 = πΎππ1ππ (π ) (37)
The transfer function of the elevation axis plant is given by:
π(π )
ππ (π )=
πΎπππ
π½ππ 2 (38)
The transfer function of the PID controller is as follows:
ππ (π )
πΈπ(π )=
πΎπππ 2+πΎπππ +πΎππ
π (39)
ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 10, No. 2, April 2020 : 1868 - 1884
1878
where πΈπ(π ) = ππ(π ) β π(π ), the open loop transfer function of the elevation axis control πΊπ(π ) is given by:
πΊπ(π ) =π(π )
πΈπ(π )=
ππ (π )
πΈπ(π )
π(π )
ππ (π ) (40)
Based on (38) and (39), the open loop elevation transfer function becomes:
πΊπ(π ) =πΎπππ
π½ππ 2 πΎπππ 2+πΎπππ +πΎππ
π (41)
The closed loop transfer function for elevation angle control is as follows:
π(π )
ππ(π )=
πΊπ(π )
1+πΊπ(π )=
πΎππππΎπππ 2+πΎππππΎπππ +πΎππππΎππ
π½ππ 3+πΎππππΎπππ 2+πΎππππΎπππ +πΎππππΎππ (42)
Figure 8. Control system block diagram for helicopter elevation, pitch and travel axis using PID controller
5.1.2. Pitch control using PD controller
Similarly, the difference of the rows of (33) results in
π’1 β π’2 = β2π12π β 2π14οΏ½ΜοΏ½ β 2π15(π β ππ) β 2π17πΎ (43)
ππ = β2π12π β 2π14οΏ½ΜοΏ½ β 2π15(π β ππ) β 2π17 β«(π β ππ)ππ‘ (44)
Substitution (25b) in (24b) results,
ππ = βπΎπππ β πΎπποΏ½ΜοΏ½ + πΎπππΎππ(ππ β π) + πΎπππΎππ β«(ππ β π)ππ‘ (45)
It is clear that (44) and (45) have exactly the same structure. Then, by comparing these equations, it can
obtain the feedback gains parameters for the PID controller from the LQR gains parameters as follows:
πΎππ = 2π12
πΎππ = 2π14
πΎππ =π15
π12
πΎππ =π17
π12
(46)
Taking Laplace transform for pitch axis model (14) yields:
π½ππ(π )π 2 = πΎπππππ(π ) (47)
The transfer function for pitch axis model is given by:
Ο(s)
Vd(s)=
πΎπππ
π½πs2 (48)
Int J Elec & Comp Eng ISSN: 2088-8708
Elevation, pitch and travel axis stabilization of 3DOF helicopter with⦠(Ibrahim K. Mohammed)
1879
The transfer function of the PD controller is as follows:
ππ(π )
πΈπ(π )= πΎπππ + πΎππ (49)
where πΈπ(π ) = ππ(π ) β π(π ), the open loop transfer function of the pitch axis control is given by:
πΊπ(π ) =π(π )
πΈπ(π )=
π(π )
ππ(π )
ππ(π )
πΈπ(π )=
πΎπππ(πΎπππ +πΎππ)
π½ππ 2 (50)
The closed loop transfer function of pitch angle is given by:
π(π )
ππ(π )= β
πΎππππΎπππ +πΎππππΎππ
π½ππ 2+πΎππππΎπππ +πΎππππΎππ (51)
5.1.3. Travel control using PI controller
Taking Laplace transform for travel axis model (15) results:
π(π )π =πππ
π½ππ(π ) (52)
The transfer function for travel axis model is given by:
r(s)
Ο(s)=
πππ
π½πs (53)
The transfer function of the PI controller is as follows:
π(π )
πΈπ(π )= πΎππ +
πΎππ
π (54)
where πΈ(π ) = ππ(π ) β π(π ), the open loop transfer function of the travel axis control is given by:
πΊπ(π ) =r(s)
π(s)
ππ(π )
πΈπ(π )=
πππ(πΎπππ +πΎππ)
π½ππ 2 (55)
The closed loop transfer function for travel angle is as follows:
r(s)
rd(s)=
ππππΎπππ +ππππΎππ
π½ππ 2+ππππΎπππ +ππππΎππ (56)
6. CONTROLLER SIMULATION AND RESULTS
6.1. LQR controller
In order to validate the proposed helicopter stabilizing system, the LQR controller is analysed
mathematically using Matlab environment. Based on objective function (J) and using the Matlab command
βlqrβ the elements of the LQR weighting matrices Q, R are tuned using GA optimization method.
The parameters of the GA optimization approach chosen for tuning process of the 3DOF helicopter control
system are listed in Table 3. Converging elements of the LQR weighting matrices Q, R through iteration
based on GA optimization method are presented in Figure 9.
Table 3. Parameters of GA tuning method GA Parameter Value/Method
Population Size 20
Max No. of Generations 100
Selection Method Normalized Geometric Selection
Crossover Method Scattering
Mutation Method Uniform Mutation
ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 10, No. 2, April 2020 : 1868 - 1884
1880
Figure 9. Number of generation of GA-LQR parameters Q and R
Based on the fitness function, the optimal elements values for LQR weighting matrices Q, R
obtained based on the GA tuning approach are given by:
π =
[ 9.0425 0 0 0 0 0 0
0 0.7840 0 0 0 0 00 0 4.5093 0 0 0 00 0 0 0.2763 0 0 00 0 0 0 0.7077 0 00 0 0 0 0 0.9876 00 0 0 0 0 0 0.8004]
, R = [0.0483 0
0 0.0483]
By using the state matrix (π΄), input matrix (π΅), and weighting matrices (π, π ), the optimized feedback gain
matrix πΎ and closed-loop rictions of the system are determined using the Matlab command βlqrβ.
The calculated gain matrix πΎ is given by:
πΎ = [117.138 65.9975 68.4252 16.9323 50.6954 31.954 28.766117.138 β65.9975 68.4252 β16.9323 β50.6954 31.954 β28.766
]
The poles of the closed-loop control system for helicopter using optimized GA-LQR controller are listed in
Table 4. Based on the feedback gain matrix and using (26), the LQR control effort vector for
the 3DOF helicopter system is as follows:
Table 4. Poles of helicopter control system based on LQR technique Pole Value π1 -794.7130 π2 -1.3745
π3 -0.3405
π4 -2161.7
π5 -1.3623 + i1.0631
π6 -1.3623 - 1.0631i
π7 -1.1773
[π’1
π’2] = β [
117.138 65.998 68.425 16.932 50.695 31.954 28.766117.138 β65.998 68.425 β16.932 β50.695 31.954 β28.766
]
[ πππΜοΏ½ΜοΏ½πΚπΎ]
(57)
Int J Elec & Comp Eng ISSN: 2088-8708
Elevation, pitch and travel axis stabilization of 3DOF helicopter with⦠(Ibrahim K. Mohammed)
1881
6.2. PID controller
Based on (36), (46) and (57), the absolute values of PID, PD and PI gain parameters for elevation,
pitch and travel axis model respectively for helicopter system are listed in Table 5 [1]. Using the values in
Table 2 and 5, the closed-loop transfer function of elevation, pitch and travel axis (42), (51) and (56) become
as in (58), (59) and (60) respectively:
Table 5. Values of gain parameters for PID, PD and PI controllers PID Parameter Relationship Absolute Value
πΎππ 2π11 234.2765
πΎππ 2π13 136.8504
πΎππ 2π16 63.9086
πΎππ 2π12 131.9949
πΎππ 2π14 33.8645
πΎππ π15/π12 0.7681
πΎππ π17/π12 0.4359
π(π )
ππ(π )=
1445π 2+247π +674.9
1.8145π 3+1445π 2+2474π +674.9 (58)
π(π )
ππ(π )=
69.08π +269.3
0.0319π 2+69.08π +269.3 (59)
π(π )
ππ(π )=
2.878π +1.634
1.815π 2+2.878π +1.634 (60)
Based on bounded input signal, the elevation, pitch and travel axis model of 3DOF helicopter
system are unstable as they give unbounded outputs. The output time responses for elevation, pitch and travel
angle are illustrated in the Figure 10. In this study, in order to to achieve a stable output, a hybrid control
system using LQR based PID controller for 3DOF helicopter system is proposed to control the dynamic
behavior of the system. To validate the proposed helicopter stabilizing scheme, the controller is simulated
using Matlab programming tool. Three axis, elevation, pitch, travel rate, are considered in the simulation
process of the control system. The performance of the helicopter balancing system is evaluated under unit
step reference input using rise, settling time, overshoot and steady state error parameters for the elevation, pitch and travel angles to simulate the desired command given by the pilot.
6.2.1. Elevation axis model simulation
This section deals with the simulatin of LQR based PID controller used to control the position of
helicopter elevation model. Figure 11 presents tracking control curve of the demand input based on PID
controller for helicopter elevation angle. The simulation result shows that the controller succeeded to guide
the output state of the system through the desired input trajectory effectively with negligible overshoot, short
rise and settling time of 0.1 ms and 0.3 ms respectively.
6.2.2. Pitch axis model simulation
In this section, GA-LQR based PD controller is designed to control the dynamic model of helicopter
pitch angle. Based on the optimized PD gain parameters stated in Table 4, the output response of proposed
helicopter tracking system is illustrated in Figure 12. It is obvious that the controller forced the pitch angle
state to follow the desired trajectory effectively.
6.2.3. Travel axis model simulation
The control of travel rate for 3DOF helicopter system is governed GA-LQR based PI controller.
The output response of the PI tracking scheme using optimum gain parameters listed in Table 5 is shown in
Figure 13. It can be noted that the optimized hybrid controller enabled the output state to track the desired
input trajectory without overshoot, and shorter rise and settling time with minimal steady state tracking error.
Regarding the control effort, the curves of input signals supplied to the propeller motors for the proposed
3DOF helicopter system are shown in Figure 14. It can be seen from input response that the control inputs of
the helicopter control system were within acceptable values. Based on the Figures 11-13, it can say that
the control performance of optimized GA-LQR based PID, PD and PI controllers for helicopter elevation,
pitch and travel axis model respectively was acceptable through tracking the system output states for
the reference input efficiently. From the mini plots of Figures 11-13, the performance parameters of PID,
PD and PI controller for helicopter elevation, pitch and travel axis are listed in the Table 6. The results in
the table show the effectiveness of the optimized hybrid controllers for helicopter system application.
ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 10, No. 2, April 2020 : 1868 - 1884
1882
Figure 10. Open-loop response of helicopter system
Figure 11. Output response of closed-loop helicopter
control system for elevation angle state
Figure 12. Output response of closed-loop helicopter
control system for pitch angle state
Figure 13. Output response of closed-loop helicopter
control system for travel rate state
Figure 14. Conrol input of 3DOF helicopter control
system
Table 6. Performance parameters of helicopter controller system Controller Tr (sec) Ts (sec) Mp %
Elevation PID 0.343 0.535 1.1
Pitch PD 0.582 1.05 0
Travel PI 1.17 12.4 5.29
0 5 10 15 200
2000
4000
6000
8000
10000
12000
14000
Time (s)
An
gle
Elevation angle
Pitch angle
Travel angle
0 0.1 0.2 0.30
0.5
1
0 5 10 15 200
0.2
0.4
0.6
0.8
1
Time (s)
Ele
va
tio
n a
ng
le
0 0.02 0.040
0.5
1
0 5 10 15 200
0.2
0.4
0.6
0.8
1
Time (s)
Pitch
an
gle
0 0.02 0.040
0.5
1
0 2 4 6 8 10-1.5
-1
-0.5
0
0.5
1
Time (s)
Co
ntr
ol e
ffo
rt
u1
u2
0 0.02 0.04-1.5
-1
-0.5
0
0.5
1
0 5 10 15 200
0.2
0.4
0.6
0.8
1
1.2
Tra
ve
l ra
te
Time (s)
0 5 100
0.5
1
Int J Elec & Comp Eng ISSN: 2088-8708
Elevation, pitch and travel axis stabilization of 3DOF helicopter with⦠(Ibrahim K. Mohammed)
1883
7. CONCLUSION
In this study, an efficient hybrid control system has been developed for 3DOF helicopter system.
The dynamics of elevation, pitch and travel axis for helicopter system is modeled mathematically and then
formulated in state space form to enable utilizing state feedback controller technique. In the proposed
helicopter stabilizing scheme, a combination of a conventional PID control with LQR state feedback
controller is adopted to stabilize the elevation, pitch and travel axis of the helicopter scheme. The gain
parameters of the traditional PID controller are determined from the gain matrix of state feedback LQR
controller.
In this research, the LQR controller is optimized by using GA tuning technique. The GA
optimization method has been adopted to find optimum values for LQR gain matrix elements which are
utilized to find best PID gain parameters. The output response of the optimized helicopter control system has
been evaluated based on rise time, setting time, overshoot and steady state error parameters. The simulation
results have shown the effectiveness of the proposed GA-LQR based PID controller to stabilize the helicopter
system at desired values of elevation and pitch angle and travel parameters.
REFERENCES [1] S. K., Choudhary, βLQR Based PID Controller Design for 3-DOF Helicopter System,β International Journal of
Electrical and Information Engineering, vol.8, no.8, pp. 1498-1503, 2014.
[2] S. Franko, βLQR-Based Trajectory Control of Full Envelope, Autonomous Helicopter,β World Congress of
Engineering 2009, vol. I, London, 1-3, July 2009.
[3] N. Thomas and P. Poongodi, βPosition Control of DC Motor Using Genetic Algorithm Based PID Controller,β
World Congress of Engineering 2009, vol. II, London, 1-3 July, 2009.
[4] L. Fan, E. Joo, βDesign for Auto-tuning PID Controller Based on Genetic Algorithms,β IEEE 4th International
Conference on Industrial Electronics and Applications, 2009
[5] I. Mohammed and A. Ibrahim, βDesign of optimised linear quadratic regulator for capsule endoscopes based on
artificial bee colony tuning algorithmβ, International Journal for Engineering Modeling, vol. 31, no.12, pp. 77-98,
2018.
[6] K. Ogata, Modern Control Engineering. Prentice-Hall International Upper Saddle River, NJ, 1997.
[7] A. Tewari, Modern Control Design with Matlab and Simulink, John Willey and Sons, LTD; 2005.
[8] F. Hasbullah, W. Faris, βA Comparative analysis of LQR and fuzzy logic controller for active suspension using half
car model,β IEEE international Conference on Control Automation Robotics and Vision, pp. 2415-2420, 2010.
[9] H. Liouane, I. Chiha, A. Douik, and H. Messaoud. βProbabilistic differential evolution for optimal design of lqr
weighting matrices,β IEEE International Conference on Computational Intelligence for Measurement Systems and
Applications (CIMSA), pp. 18β23, 2012.
[10] D. Srinivasan, T. H. Seow, Evolutionary Computation, CEC β03, Canberra, Australia, vol. 4, pp. 2292β2297, 2003.
[11] X. Zeng, L. Jing, Z. Yao, and Y. Guo, βA PSO-based LQR controller for accelerator PWM power supply,β 2nd
International Conference on Mechatronics and Intelligent Materials, MIM 2012, May 18-19, pp. 71-75, 2012.
[12] M. B. Poodeh, S. Eshtehardiha, A. Kiyoumarsi, and M. Ataei, βOptimizing LQR and pole placement to
control buck converter by genetic algorithm,β International Conference on Control, Automation and Systems,
ICCAS '07, pp.2195-2200, 2007.
[13] J. Zhang, L. Zhang, and J. Xie, βApplication of memetic algorithm in control of linear inverted pendulum,β IEEE
International Conference on Cloud Computing and Intelligence Systems, pp. 103β107, 2011.
[14] E. Rakhshani, βIntelligent linear-quadratic optimal output feedback regulator for a deregulated automatic
generation control system, βElectric Power Components and Systems, vol. 40, no. 5, pp. 513β533, 2012.
[15] A. Jacknoon , M. A. Abido, βAnt Colony based LQR and PID tuned parameters for controlling Inverted Pendulum,
β IEEE International Conference on Communication, Control, Computing and Electronics Engineering
(ICCCCEE), pp. 19-25, 2017.
[16] B. Ata, R. Coban, βArtificial Bee Colony algorithm based Linear Quadratic optimal controller design for a nonliear
Inverted Pendulum,β International Journal of Intelligent System and Applications in Engineering, vol. 3 no. 1,
pp. 1-6, 2015.
[17] S. P. Ramaswamy, G. K Venayagamoorthy, and S. Balakrishnan, βOptimal control of class of non-linear plants
using artificial immune systems: Application of the clonal selection algorithm,β IEEE 22nd International
Symposium on Intelligent Control, pp. 249-254, 2007. [18] E. V. Kumar and J. Jerome, βLQR based optimal tuning of PID controller for trajectory tracking of Magnetic
Levitation System, β International Conference On DESIGN AND MANUFACTURING, Procedia engineering 24,
pp. 254-264, 2013.
[19] A. N. Khizer, et al., β3DoF Model Helicopter with Hybrid Control,β TELKOMNIKA Indonesian Journal of
Electrical Engineering, vol.12, no. 5, pp. 3863-3872, 2014.
[20] R. Akmeliawati and S. Raafat, βOptimized State Feedback Regulation of 3DOF Helicopter System via Extremum
Seeking,β 9th Asian Control Conference, pp. 1-6, 2013.
[21] I. Mohammed, B. Sharif and J. Neasham, βDesign and implementation of a magnetic levitation control system for
robotically actuated capsule endoscopes,β IEEE International Symposium on Robotic and Sensors Environments
ROSE2012, pp. 140-145, 2012.
ISSN: 2088-8708
Int J Elec & Comp Eng, Vol. 10, No. 2, April 2020 : 1868 - 1884
1884
[22] I. Mohammed, A. Abdulla, βFractional Order PID Controller Design for Speed Control DC Motor based on
Artificial Bee Colony Optimization, β International Journal of Computer Applications, vol. 179, no. 24, pp. 43-49,
2018.
[23] M. Bharathi, G. Kumar, βDesign Approach for Pitch Axis Stabilization Of 3-Dof Helicopter System an LQR
Controller,β International Journal of Advanced Research in Electrical, Electronics and Instrumentation
Engineering, vol. 1, no. 5, pp. 351-365, 2012.
[24] J. Fang, βThe LQR Controller Design of Two-Wheeled Self-Balancing Robot Based on the Particle Swarm
Optimization Algorithm,β Mathematical Problems in Engineering, vol. 2014, pp. 1-6, June 2014.
[25] C. Wongsathan, C. Sirima, βApplication of GA to Design LQR Controller for an Inverted Pendulum System,β
IEEE International Conference on Robotics and Biomimetics, pp. 951-954, 2008.
[26] M. J. Neath et al., βAn Optimal PID Controller for a Bidirectional Inductive Power Transfer System Using
Multiobjective Genetic Algorithm, β IEEE Transaction on Power Electronic, vol. 29, no. 3 pp. 1523-1531, 2014.
[27] A. Abdulla., I. Mohammed and A. Jasim, βRoll Control System Design Using Auto Tuning LQR Technique,β
International Journal of ngineering and Innovative Technology, vol.6, no. 2, pp. 11-22.