Received 11 May 2016; accepted 20 May 2016; published 29 August 2016
Earlier some research was done on helicopter and its control in foreign countries  -  . The cost of helicopter model is very high where developing countries could not afford. To encourage the research on control system in academic institution in developing countries, LCADRH is designed and developed for the students. The helicopter is controlled in two degrees of freedom (2-DOF). In future, control oriented design related to 3-DOF can be extended  -  . The system identification of nonlinear system like helicopter plays a major role for its precise control subject to system dynamics and uncertainty  -  . Research work is carried out on grey box system identification for model helicopter in  . The small scale unmanned rotorcraft system is modeled using system identification for the flight control in  . Prediction error method of LPV models is done in  . LQR control was successful for aerospace engineers for the flight control provides good stability and control   . LQR controller design for helicopter control was done in 2-DOF   . Helicopter control is highly tedious due to the presence of nonlinear dynamics and uncertainty. Research work on helicopter control is necessary due to the demand of precise control in emergency conditions like earth quake, fire accidents. Nonlinear control of 2-DOF helicopter research is carried out in  -  .
In this paper, the mathematical model of LCADRH is identified using system identification. The main components of LCADRH is designed and integrated to get a helicopter prototype. The effect of nonlinearity and inherent coupling existence are observed in closed loop position control. Therefore, proposed LCADRH model can be used for providing hands-on training for undergraduate and postgraduate students in the area of modelling and control of nonlinear system. The block diagram of proposed system is shown in Figure 1. The LQR controller is designed for the closed loop position control of LCADRH.
2. Design of Mathematical Model
There are many methods available for system identification of non-linear MIMO systems using measured input-output data such as black box modelling, grey box modelling etc. In this paper, we design an indigenous archetype dual rotor helicopter and its dynamics are modelled by grey box modelling. Using the system identification tool in MATLAB an idnlgrey object of this model is created and the unknown idnlgrey model parameters and initial states using measured data are estimated by prediction error method (PEM). The PEM is a parametric estimation method in which the system parameters are found out providing the initial set of values. Here the initial values are set by physical interpretation and presumption.
The objective of the PEM is to minimize the prediction error by the prediction parameter which minimizes the variation of prediction error,
is the prediction error,
is the model output,
Figure 1. Block diagram of Low cost archetype dual rotor helicopter system for academic research.
is the predicted estimate of the output,
is the vector containing all unknown parameters,
is the vector containing predicted parameters,
is the number of data samples.
The PEM uses optimization to minimize the cost function. Initially the state space matrixes A, B contains the coefficients with initial value of fixed parameters and some fixed values. The estimating algorithm estimates the parameters in A and B matrix. The model is estimated using PEM with so far determined parameters as the initial values in MATLAB. To validate this method, a simulation model is used to generate the data for the system identification.
3. LQR Controller Design
While designing controllers for complex systems with stringent performance requirements, LQR method is used. It is a powerful method to find the best controller that minimizes cost. It is mainly used where the system dynamics are given by a set of linear equations and the cost is given by a quadratic function. Here the weighting factors are supplied by a human. The settings of the controller are found by using a mathematical algorithm that minimizes the cost function. Q and R are the two matrices that parameterize the cost function. They weight the state vector and the system input respectively. LQR method tries to achieve optimal control by solving the algebraic riccatti equation.
The LQR controller is designed for linear state space model of LCADRH,
The states of LCADRH are,
is pitch angle,
is yaw angle,
is angular pitch velocity,
is angular yaw velocity,
and the input of LCADRH are,
is pitch motor voltage
is yaw motor voltage
satisfies the following cost function,
where Q and R are the weighting matrices that to be designed such that subject to the system dynamics.
4. Low Cost Design of Dual Rotor Helicopter Prototype in 2-DOF
The main components of LCADRH are 12V DC servo motor, 6V DC servo motor, optical pitch encoder, optical yaw encoder, hollow shaft slip ring and motor power circuit. The LCADRH shown in Figure 2 is designed and integrated for the control of elevation about the pitch axis and side to side motions about yaw axis. 12V DC servo motor drives the pitch propeller for elevation and 6V DC servo motor drives the yaw propeller for side to
Figure 2. Low cost designed two rotor helicopter model for academic research.
side motions. Optical pitch encoder and optical yaw encoder provides the digital position feedback information. The hollow shaft slip ring is used to transfer power of two motors and optical pitch encoder signal without tangling of wires.
System Specification and Cost Details
The main components of LCADRH are shown in Figure 3. Table 1 summarizes the specification of components and angle range of pitch and yaw of LCADRH. The Table 2 shows the overall cost detail of LCADRH. The total cost for the complete design of LCADRH is $894.51. The designed helicopter model in low cost is beneficial for the students in academic institution to test and explore modern identification and control methodologies.
5. Results Discussion
5.1. Closed-Loop Position Control of LCADRH
The position control of LCADRH in open loop is tedious due to high nonlinearity and strong coupling effect. In open loop, LCADRH tend to change yaw position when pitch nose goes up and down and vice versa due to the existence of coupling effect between pitch and yaw actuators. The LQR closed loop position control of LCADRH is done under four cases to analyze the coupling effect that exists between the pitch and yaw actuators. In case I, the pitch angle is varied by giving step input of 10 degree whereas yaw angle is constant at 0 degree. In case II, yaw angle is varied by giving step input of 30 degree while the pitch angle is constant at 0 degree. In case III, the elevation of LCADRH is raised by setting pitch angle to square input with the amplitude of 5 degree and frequency of 0.05 Hz whereas the yaw angle is constant at 0 degree. In case IV, side to side motion is increased by setting the yaw angle to square input of 30 degree and frequency of 0.05 Hz and the pitch angle is constant at 0 degree. The objective of this closed-loop position control is to compare the measured closed-loop response with simulated response. The LQR controller gives encouraging results in reducing coupling effect, structural vibration and oscillation.
5.2. Case I―Pitch Step Input
To analyze the coupling impact of LCADRH MIMO system in LQR closed loop position control, the first control variable is varied by keeping second control variable as constant at 0 degree. The simulated pitch and yaw position response tracks the reference with less steady state error. The measured pitch angle is 10.2 degree at t = 1.478 s. The steady state error of measured pitch angle is very less at time t = 20 s to t = 30 s as
(a) (b)(c) (d)
Figure 3. Main components of LCADRH (a) Pitch encoder; (b) Yaw encoder; (c) Pitch and Yaw motor; (d) Slip ring.
Table 1. LCADRH specification.
Table 2. Cost detail of LCADRH.
shown in Figure 4(a). Moreover deviation of measured pitch angle is also very less as shown in Figure 4(a). The difference between the simulated and measured pitch angle is 0.5 degree. The experimental result of pitch angle tracks the desired pitch angle consistently. By conducting this experiment, it is observed that there is less coupling impact of elevation motion of LCADRH on yaw axis as shown in Figure 4(a). Figure 4(a) depicts
Figure 4. LCADRH MIMO system outputs (a) LQR controlled closed loop position response under pitch step input (b) and (c) Motor input voltages.
LCADRH rotates about the yaw axis by 3.56 degree at t = 0.622 s due to the rise of elevation about pitch axis by 10 degree. This is due to the generation of rotary force by pitch motor about the yaw axis. From t = 10 s to t = 50 s, the measured yaw angle is constant at 3.516 degree. It is observed from the results that the implementation of LQR controller for the closed loop position control of LCADRH gives better result in tracking the desired pitch angle. From the experimental result of pitch and yaw angle, it is observed that structural vibration is purged. Therefore LCADRH is designed well at low cost.
The Figure 4(b) and Figure 4(c) shows the input voltage of pitch and yaw motor. Initially when the elevation of LCADRH is raised about an pitch axis at time t = 0 s, the pitch motor voltage is 7.334 V. LCADRH goes down about an pitch axis by the decrease of yaw motor input voltage from −1.959 V at t = 0 s to −9.056 V at t = 11.99 s. The measured pitch and yaw motor input voltage in Figure 4(b) and Figure 4(c) shows less variation.
5.3. Case II―Yaw Step Input
In this case, coupling existence between the pitch and yaw actuators of LCADRH is observed by giving step input of 30 degree about yaw axis. The simulated pitch and yaw angle possess consistent tracking of desired pitch and yaw angle as shown in Figure 5(a). The pitch motor simulated voltage of 10 V causes the overshoot on yaw axis by angle 5.28 degree from reference. Then eventually simulated yaw angle tracks the reference with the error of 0.07 degree from t = 30 s to t = 50 s. As shown in Figure 5(a), initially the measured pitch angle has overshoot of 9.58 degree when the yaw motor input voltage is 7.034 V. This yaw motor voltage is the reason for the generation of torque on pitch axis. Experimental result depicts then the pitch angle eventually stabilizes at t=6.186s. The measured pitch angle is 0 degree at t = 10 s to t = 50 s. Experimental result shows that the absence of coupling impact of yaw actuator on pitch axis. The 10 V applied to the pitch motor at t = 0.1 s generates a torque on yaw axis. The measured yaw angle tracks the desired yaw angle with the error of 2.43 degree which is very less.
The Figure 5(b) and Figure 5(c) shows the pitch an yaw motor input voltages. As illustrated in Figure 5(b), motion about the pitch axis tracks the reference exactly when pitch motor voltage is approximately −1.635 V.
Figure 5. LCADRH MIMO system outputs (a) LQR controlled closed loop position response under yaw step input (b) and (c) Motor input voltages.
This shows the absence of coupling effect of yaw actuator on pitch axis in LQR closed loop position control of LCADRH.
5.4. Case III―Pitch Square Input
In case III, square input of 5 degree with a periodicity of 20 seconds is given as desired set point for pitch to analyze the forward and backward movement of LCADRH and its corresponding coupling impact on yaw. The simulated result of pitch angle and yaw angle follows the desired set point of pitch and yaw angle. As shown in Figure 6(a), experimental result illustrates that the rotation of LCADRH about yaw axis is by 5.537 degree at t = 20.62 s due to the rise of elevation about pitch axis. When the desired pitch angle is raised from −5 degree to 5 degree and decreased from 5 degree to −5 degree, measured yaw angle rotates clockwise and anticlockwise by 1 degree. This measure of side to side motion of LCADRH proves the inherent presence of coupling existence between pitch and yaw actuators. Initially, the negative potential voltage is applied to the yaw motor as shown in Figure 6(c).
5.5. Case IV―Yaw Square Input
In case IV, square input of 30 degree with a periodicity of 20 seconds is given as desired set point for yaw to analyze the side to side motions of LCADRH and its corresponding coupling impact on pitch actuator. Figure 7 shows that the simulated pitch and yaw response tracks desired pitch and yaw angle with less steady state error. For a given desired yaw angle, the body of LCADRH rotated clockwise and in an anticlockwise direction with less steady state error. The coupling impact of yaw actuator on pitch actuator is reduced by the LQR position controller.
The LCADRH is designed and developed for the academic research on control system experiments in educational institution. The real time measurement and control of helicopter is provided by integrating with data
Figure 6. LCADRH MIMO system outputs (a) LQR controlled closed loop position response under pitch square input (b) and (c) Motor input voltages.
Figure 7. LCADRH MIMO system outputs (a) LQR controlled closed loop position response under yaw square input (b) and (c) Motor input voltages.
acquisition board and sensors to provide hardware features, software support capability. The mathematical modelling of LCADRH is identified from its real input output data using grey box modelling. LQR controller is designed for LCADRH for the pitch and yaw position control. The coupling effect, nonlinearity and system dynamics of LCADRH is analyzed using experimental setup by giving step and square wave input. The simulation results of pitch and yaw angle follows the desired angle with less steady state error. The experimental results of pitch and yaw position shows the absence of structural vibration, oscillation and measurement noise. The designed LQR position controller controls elevation and side to side motion with consistent tracking of desired pitch and yaw angle and it reduces the existence of coupling effect between pitch actuator and yaw actuator.
Because of its cost effective nature, the proposed LCADRH model can be developed in third world countries giving the under graduate and post graduate students an affordable opportunity to learn modelling, simulation, and control oriented experiments in a non linear coupled MIMO system.
The authors gratefully acknowledge Dr. S. Baskar, Professor, EEE, Thiagarajar College of Engineering, Madurai, for his meticulous guidance, encouragement and his valuable discussions. One of the authors of this paper (P.S. Manoharan) acknowledges Science and Engineering Research Board, Department of Science and Technology (DST), India for sanctioning the funding under Fast Track Young Scientist Scheme, vide sanction number SERB/F/2056/11-12 dated 15.02.2012.
 Ahmad, S.M., Chipperfield, A.J. and Tokhi, O. (2000) Dynamic Modeling and Optimal Control of a Twin Rotor MIMO System. IEEE Conference on National Aerospace and Electronics, 391-398.
 Liu, H., Lu, G. and Zhong, Y.S. (2013) Robust LQR Attitude Control of a 3-DOF Laboratory Helicopter for Aggressive Maneuvers. IEEE Transactions on Industrial Electronics, 60, 4627-4636.
 Mettler, B., Tischler, M.B. and Kanade, T. (2002) System Identification Modeling of a Small Scale Unmanned Rotorcraft for Flight Control Design. Journal of the American Helicopter Society, 47, 50-63.
 Zhao, Y., Huangb, B., Su, H.Y. and Chu, J. (2012) Prediction Error Method for Identification of LPV Models. Journal of Process Control, 22, 180-193.
 Vinodh Kumar, E. and Jerome, J. (2013) Robust LQR Controller Design for Stabilizing and Trajectory Tracking of Inverted Pendulum. Procedia Engineering, 64, 169-178.
 Yu, G.-R. and Liu, H.-T. (2005) Sliding Mode Control of a Two-Degree-of-Freedom Helicopter via Linear Quadratic Regulator. Systems, Man and Cybernetics, 4, 3299-3304.
 Ramalakshmi, A.P.S. and Manoharan, P.S. (2012) Non-Linear Modeling and PID Control of Twin Rotor MIMO System. Proceedings of IEEE Conference on Advanced Communication Control and Computing Technologies, 366-369.
 Juang, J.-G., Liu, W.-K. and Lin, R.-W. (2011) A Hybrid Intelligent Controller for a Twin Rotor MIMO System and Its Hardware Implementation. ISA Transactions, 50, 609-619.