<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">WJET</journal-id><journal-title-group><journal-title>World Journal of Engineering and Technology</journal-title></journal-title-group><issn pub-type="epub">2331-4222</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/wjet.2014.24034</article-id><article-id pub-id-type="publisher-id">WJET-51808</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Chemistry&amp;Materials Science</subject><subject> Engineering</subject></subj-group></article-categories><title-group><article-title>
 
 
  Model Helicopter Control Using Body-Mounted Vibro-Tactile Transducers
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>ui</surname><given-names>Zhou</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Merhran</surname><given-names>Mehrandezh</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Raman</surname><given-names>Paranjape</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Industrial Systems Program, Faculty of Engineering, University of Regina, Regina, Canada</addr-line></aff><aff id="aff1"><addr-line>Electronic Systems Program, Faculty of Engineering, University of Regina, Regina, Canada</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>Erui198311@gmail.com(UZ)</email>;<email>Mehran.Mehrandezh@uregina.ca(MM)</email>;<email>Raman.Paranjape@uregina.ca(RP)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>29</day><month>09</month><year>2014</year></pub-date><volume>02</volume><issue>04</issue><fpage>322</fpage><lpage>330</lpage><history><date date-type="received"><day>27</day>	<month>September</month>	<year>2014</year></date><date date-type="rev-recd"><day>20</day>	<month>October</month>	<year>2014</year>	</date><date date-type="accepted"><day>10</day>	<month>November</month>	<year>2014</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  The sense of touch as a man-machine communication channel can be as acute as the sense of sight and sound. In some scenarios such as those seen in aerobatics, stunt flying, and combat flights, tactile sensors can even outperform the conventional non-contact sensors in terms of situation awareness. Fusion of tactile sensory information with those obtained via sight and sound can avoid diverting the user’s attention away from the operational task at hand as well. In this study, the performance of an operator, to servo control the motion of a 2-dof model helicopter with pitch/yaw maneuverability, subjected to an intuitive body-referenced arrangement of a cluster of vibro-tactile sensors is investigated. A blindfolded operator will then control the helicopter to a safe attraction zone via a joystick based on this tactile sensory information. A fine-tuned local controller would take over for the end-of-motion precise homing. This study can pave the way towards a systematic integration and characterization of tactile sensors in high performance weapon platforms with improved situation awareness in visually awkward maneuvers such as those seen in aerial combat scenarios.
 
</p></abstract><kwd-group><kwd>Haptics</kwd><kwd> Flight Control</kwd><kwd> Linear Quadratic Regulators (LQR)</kwd><kwd> Tactile Immersive Control</kwd><kwd> 2-Dof  Model Helicopter</kwd><kwd> Mechatronics</kwd><kwd> Real-Time Control</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Highly maneuverable flying machines such as helicopters are widely used in search and rescue, and surveillance due to their agility [<xref ref-type="bibr" rid="scirp.51808-ref1">1</xref>] . Control of these machines lends itself as a challenging Multi-Input Multi-Output (MIMO) problem due to the nonlinear dynamics and strong interaction between controlled variables. Incorporating haptic sensors into the closed-loop feedback control of flying machines has recently drawn a great deal of attention in both academia and the industry sector [<xref ref-type="bibr" rid="scirp.51808-ref2">2</xref>] . When employing body-referenced haptic sensors, the man-machine interaction can be as acute as that when using sense of sight and/or sound. In [<xref ref-type="bibr" rid="scirp.51808-ref3">3</xref>] , a haptic interface for UAV collision avoidance was used. It was reported that, with the limited visual information, the collision could be avoided using haptic clues. Literature pertinent to the control of highly maneuverable aircrafts is vast. In particular, the robust and optimal controls of helicopters have been studied for the past two decades, (e.g., Refs. [<xref ref-type="bibr" rid="scirp.51808-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.51808-ref4">4</xref>] ).</p><p>An effective human-in-the-loop hybrid control strategy based on sensory information obtained via an array of body-mounted haptic sensors (i.e., vibrating motors) is proposed. This method was successfully tested on a 2- dof model helicopter. The control strategy is summarized below.</p><p>At the first stage, a user attempts to navigate the model helicopter to an attraction zone through haptic sensory feedback by using a joystick. In the second stage, an autonomous nonlinear controller takes over making the helicopter hover at steady state with high precision. A simple switching algorithm is adopted for smooth transi- tion between the human- and computer-based controls.</p><p>The experimental model helicopter used in this study has two degrees of freedom, namely pitch and yaw (see <xref ref-type="fig" rid="fig1">Figure 1</xref>). The pitch propeller is driven by a voltage-controlled DC motor. Nose up is considered as positive pitch. Correspondingly, the yaw propeller is driven by a smaller voltage-controlled DC motor. Due to the gy- roscopic effect, the propeller causing pitch would also cause a load torque on the motor’s shaft which is in turn seen at the yaw axis. The model helicopter is pivoted on a fixed stand. Its center of gravity is located off the pi- vot point with a considerable offset in a way that the gravitation force will push the helicopter’s nose down all the time. Pitch and yaw angles are measured through optical encoders. The pitch angle is limited to &#177;40 degrees from the level configuration, while, thanks to the slip-ring employed at the pivot point, the yaw angle can vary, and being measured by an optical encoder, indefinitely.</p><p>The layout of the paper is as follows: The governing dynamics equations of the system are described in Sec- tion 2. The control strategy is explained in Section 3. The experimental setup and results are given in Sections 4 and 5, respectively. Section 6 concludes the paper and described the future works.</p></sec><sec id="s2"><title>2. Dynamics Model of the Helicopter</title><p><xref ref-type="fig" rid="fig1">Figure 1</xref> shows a schematic of the main components of the model helicopter. The governing dynamics equa- tions are as follows [<xref ref-type="bibr" rid="scirp.51808-ref5">5</xref>] :</p><disp-formula id="scirp.51808-formula1875"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-1560137x5.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.51808-formula1876"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-1560137x6.png"  xlink:type="simple"/></disp-formula><p>Definition of parameters in Equations (1) and (2) are given in <xref ref-type="table" rid="table1">Table 1</xref>.</p><p>The nonlinear governing dynamics given in Equations (1) and (2) can be linearized about the quiescent point</p><p>defined at<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x7.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x8.png" xlink:type="simple"/></inline-formula>, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x9.png" xlink:type="simple"/></inline-formula> as follows:</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> A schematic of the 2-dof model helicopter</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-1560137x10.png"/></fig><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Physical parameters of the system</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Symbol</th><th align="center" valign="middle" >Definition</th><th align="center" valign="middle" >Value</th><th align="center" valign="middle" >Unit</th></tr></thead><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x11.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >Moment of inertia (pitch direction)</td><td align="center" valign="middle" >0.0384</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x12.png" xlink:type="simple"/></inline-formula></td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x13.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >Moment of inertia (yaw direction)</td><td align="center" valign="middle" >0.0432</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x14.png" xlink:type="simple"/></inline-formula></td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x15.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >Mass of the helicopter</td><td align="center" valign="middle" >1.3872</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x16.png" xlink:type="simple"/></inline-formula></td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x17.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >Distance from center of mass to yaw axis</td><td align="center" valign="middle" >0.186</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x18.png" xlink:type="simple"/></inline-formula></td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x19.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >Thrust torque constant acting on pitch axis from pitch motor</td><td align="center" valign="middle" >0.204</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x20.png" xlink:type="simple"/></inline-formula></td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x21.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >Thrust torque constant acting on yaw axis from yaw motor</td><td align="center" valign="middle" >0.072</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x22.png" xlink:type="simple"/></inline-formula></td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x23.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >Thrust torque constant acting on pitch axis from yaw motor</td><td align="center" valign="middle" >0.0068</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x24.png" xlink:type="simple"/></inline-formula></td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x25.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >Thrust torque constant acting on yaw axis from pitch motor</td><td align="center" valign="middle" >0.0219</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x26.png" xlink:type="simple"/></inline-formula></td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x27.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >Power voltage acting on the pitch motor</td><td align="center" valign="middle" >Variable</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x28.png" xlink:type="simple"/></inline-formula></td></tr><tr><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x29.png" xlink:type="simple"/></inline-formula></td><td align="center" valign="middle" >Power voltage acting on the yaw motor</td><td align="center" valign="middle" >Variable</td><td align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x30.png" xlink:type="simple"/></inline-formula></td></tr></tbody></table></table-wrap><disp-formula id="scirp.51808-formula1877"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-1560137x31.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.51808-formula1878"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-1560137x32.png"  xlink:type="simple"/></disp-formula><p>The linearized dynamics equations can be presented in the state space, disregarding the gravitational torque</p><p>term, namely <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x33.png" xlink:type="simple"/></inline-formula> as follows:</p><disp-formula id="scirp.51808-formula1879"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-1560137x34.png"  xlink:type="simple"/></disp-formula><p>where the state and control vectors are: <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x35.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x36.png" xlink:type="simple"/></inline-formula> respectively.</p><p>As can be seen from Equation (1), the gravitational force generates a torque around the pivot point as</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x37.png" xlink:type="simple"/></inline-formula>. In order to compensate this torque, a non-linear term, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x38.png" xlink:type="simple"/></inline-formula>was fed-forward to the voltage</p><p>provided to the pitch DC motor as follows:</p><disp-formula id="scirp.51808-formula1880"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-1560137x39.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x40.png" xlink:type="simple"/></inline-formula> is:</p><disp-formula id="scirp.51808-formula1881"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-1560137x41.png"  xlink:type="simple"/></disp-formula><p>In Equation (7), <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x42.png" xlink:type="simple"/></inline-formula>denotes the feed-forward control gain (normally considered as 1) and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x43.png" xlink:type="simple"/></inline-formula> denotes the desired pitch angle. The gravitational force would have no effect on yaw motion. Therefore, the overall</p><p>pitch/yaw control input denoted by <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x44.png" xlink:type="simple"/></inline-formula> can be presented as:</p><disp-formula id="scirp.51808-formula1882"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-1560137x45.png"  xlink:type="simple"/></disp-formula></sec><sec id="s3"><title>3. The Proposed Control Strategy</title><p>The control strategy is carried out in two phases. In the first phase, a human operator attempts to bring the robot within an attraction zone using a joystick via haptic sensory information. We call this a human-in-the-loop con- trol strategy. In the second phase, a nonlinear controller takes over for the precise end-of-motion servoing pur- pose. A simple switching law is adopted for smooth transition between the two controllers. The control strate- gies and the switching law are described in the following sections.</p><sec id="s3_1"><title>3.1. Human-in-the-loop Control</title><p>The key to the successful integration of vibro-tactile sensory information into a control loop is to convey an or- ganized tactile sensation to the user via lightweight and compact devices mounted on user’s body without im- pairing his/her movement. In this context, an array of four small-size vibrating DC motors (denoted as front, back, left and right motors) was mounted on the user’s body (see <xref ref-type="fig" rid="fig2">Figure 2</xref>). The voltage provided to the qua- druplet vibro-motors, and correspondingly their vibration amplitudes, will change proportional to the error and also rate of change of error between the pitch and yaw angles and their desired values. Mathematically this can be written as:</p><disp-formula id="scirp.51808-formula1883"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-1560137x46.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x47.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x48.png" xlink:type="simple"/></inline-formula> denote the voltages provided to the two pairs (i.e., front-back and left-right) of the body- mounted vibro-motors corresponding to pitch and yaw errors, respectively, the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x49.png" xlink:type="simple"/></inline-formula><sub> </sub>denotes the calibration constant, the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x50.png" xlink:type="simple"/></inline-formula> denotes the Linear Quadratic Regulator (LQR) gain, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x51.png" xlink:type="simple"/></inline-formula> denotes the desired state of the helicopter. The arrangement of the voltages provided to the quadruplet vibro motors are as follows:</p><disp-formula id="scirp.51808-formula1884"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-1560137x52.png"  xlink:type="simple"/></disp-formula><p>In Equation (10), <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x53.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x54.png" xlink:type="simple"/></inline-formula><sub>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x54.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x55.png" xlink:type="simple"/></inline-formula></sub>, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x54.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x55.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x56.png" xlink:type="simple"/></inline-formula> denote the voltages provided to the front, back, right and left motors mounted on operator’s body as seen in <xref ref-type="fig" rid="fig2">Figure 2</xref>. In the next section, a description on calculating the <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x53.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x54.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x55.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x56.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x57.png" xlink:type="simple"/></inline-formula> is given.</p><sec id="s3_1_1"><title>3.1.1. The Linear Quadratic Regulator, an overview</title><p>In a Linear Quadratic Regulator (LQR), the control gain, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x58.png" xlink:type="simple"/></inline-formula>in a linear full-state feedback control law of the form <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x58.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x59.png" xlink:type="simple"/></inline-formula> is calculated by minimizing a quadratic objective function, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x58.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x59.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x60.png" xlink:type="simple"/></inline-formula>defined as:</p><disp-formula id="scirp.51808-formula1885"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-1560137x61.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x62.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x62.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x63.png" xlink:type="simple"/></inline-formula> are weighting matrices on states and the control variables, respectively. The following weight- ing matrices were adopted in our experiments:</p><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Four vibro-tactile sensors mounted on human’s body</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-1560137x64.png"/></fig><disp-formula id="scirp.51808-formula1886"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-1560137x65.png"  xlink:type="simple"/></disp-formula><p>The LQR control gain matrix, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x66.png" xlink:type="simple"/></inline-formula>was calculated as:</p><disp-formula id="scirp.51808-formula1887"><label>(13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-1560137x67.png"  xlink:type="simple"/></disp-formula></sec><sec id="s3_1_2"><title>3.1.2. The Switching Law</title><p>The user controls the motion of the model helicopter based on vibro-tactile sensation he/she receives through four body-mounted motors. A human-based control would shift to an autonomous control strategy after helicop- ter is brought to a safe zone (i.e., attraction zone). A weighted 2-norm of the state error as a performance index, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x68.png" xlink:type="simple"/></inline-formula>was used to trigger the switching between the human- and computer-based controllers. Mathematically it can be represented as:</p><disp-formula id="scirp.51808-formula1888"><label>, (14)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-1560137x69.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x70.png" xlink:type="simple"/></inline-formula> denotes the weighting factors, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x71.png" xlink:type="simple"/></inline-formula>denote the states, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x72.png" xlink:type="simple"/></inline-formula> denotes the desired value of states.</p><p>Weighting factors, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x73.png" xlink:type="simple"/></inline-formula>and the calibration constant, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x73.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x74.png" xlink:type="simple"/></inline-formula>were tuned offline during operator’s training based on his/her level of comfort when controlling the model helicopter. The overall structure of the proposed control strategy is depicted in <xref ref-type="fig" rid="fig3">Figure 3</xref>.</p></sec></sec><sec id="s3_2"><title>3.2. End-of-motion Servo Control</title><p>A full-state feedback control law similar to that in LQR was developed for the end-of-motion servoing. However, a loop-shaping optimization process was utilized to calculate the controller gains. The nonlinear model of the helicopter was assumed within the optimization process. Furthermore, the pitch/yaw motor saturation was taken into consideration through the optimization process as well. Geometric hard constrains on parameters such as settling time, rise time and the overshoot in system’s response were defined offline prior to optimization. The Response Optimization Toolbox from Math Work was used for this purpose [<xref ref-type="bibr" rid="scirp.51808-ref6">6</xref>] . As an example, <xref ref-type="fig" rid="fig4">Figure 4</xref> shows the simulated yaw response of the model helicopter when using the LQR design described in Section (3.1.1). <xref ref-type="fig" rid="fig5">Figure 5</xref> shows the system’s response based on the response optimization technique described earlier. As can be seen from <xref ref-type="fig" rid="fig5">Figure 5</xref>, a lesser overshoot, settling time, and rise time were achieved without saturating the yaw motor. <xref ref-type="fig" rid="fig6">Figure 6</xref> shows the geometric constrains imposed on the system’s response along with a number of ite- rations on calculated response trajectories. It is noteworthy that two integral terms on pitch and yaw angles were also added to the state vector further on to reduce the steady state error. The modified state vector <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x75.png" xlink:type="simple"/></inline-formula> was:</p><disp-formula id="scirp.51808-formula1889"><label>(15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-1560137x76.png"  xlink:type="simple"/></disp-formula><p>The optimal gain <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x77.png" xlink:type="simple"/></inline-formula> was calculated as:</p><disp-formula id="scirp.51808-formula1890"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/9-1560137x79.png"  xlink:type="simple"/></disp-formula></sec></sec><sec id="s4"><title>4. Experimental setup</title><p><xref ref-type="fig" rid="fig7">Figure 7</xref> shows the experimental apparatus. Two power modules were used to drive pitch and yaw motors of the model helicopter. Linear Pulse Width Modulated (PWM) servo amplifiers were also used to drive body-mounted vibrating motors. The motion of the helicopter can be controlled in real time via a Logitech joystick (shown in the figure) or completely autonomously. The Q8 board from Quanser was used for control and data acquisition in real time [<xref ref-type="bibr" rid="scirp.51808-ref7">7</xref>] . Wincon from Quanser along with a third-party real-time kernel provide a real-time operating environment within the Simulink. The Q8 board can support 8 analogue outputs and 8 quadratic encoder chan-</p><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> The proposed control structure</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-1560137x80.png"/></fig><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Simulated yaw response of the 2-dof model helicopter via LQR + Integrator</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-1560137x81.png"/></fig><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Simulated yaw response of the 2-dof model helicopter via loop- shaping and response optimization</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-1560137x82.png"/></fig><fig id="fig6"  position="float"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> User-defined geometric hard constraints on the 2-dof model heli- copter’s yaw response</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-1560137x83.png"/></fig><fig id="fig7"  position="float"><label><xref ref-type="fig" rid="fig7">Figure 7</xref></label><caption><title> The experimental apparatus</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-1560137x84.png"/></fig><p>nels. Two rotary optical encoders measure the pitch and yaw angles in real time. Pitch motion is limited to <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x85.png" xlink:type="simple"/></inline-formula> degrees from the level configuration. However, the yaw angle can change indefinitely thanks to the slip-ring mechanism used at the pivot point connecting the stand and the helicopter’s main body. The pitch motor can generate a larger thrust than that in the yaw motor. This extra power would be needed to compensate for the gravitational torque around the pivot due to the offset between the helicopter’s center of gravity and the hinge on the helicopter’s stand.</p></sec><sec id="s5"><title>5. Experimental Result</title><p><xref ref-type="fig" rid="fig8">Figure 8</xref> and <xref ref-type="fig" rid="fig9">Figure 9</xref> show two representative experimental results. The desired pitch and yaw angles were set at 0 and 90 degrees, respectively in both experiments. Pitch and yaw angles along with voltages provided to the pitch/yaw motors are shown in figures. The initial elapsed time for controllers to take effect in both experiments was set at 3 seconds approximately. In the first experiment shown in <xref ref-type="fig" rid="fig8">Figure 8</xref>, the human-in-the-loop control strategy is in place for 15 seconds following the initialization. The computer-based controller starts at 18 se- conds introducing some small overshoots in both pitch and yaw motions before bringing the helicopter to steady state.</p><fig id="fig8"  position="float"><label><xref ref-type="fig" rid="fig8">Figure 8</xref></label><caption><title> Experimental results #1 (original in color)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-1560137x86.png"/></fig><fig id="fig9"  position="float"><label><xref ref-type="fig" rid="fig9">Figure 9</xref></label><caption><title> Experimental result after changing parameters (original in color)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/9-1560137x87.png"/></fig><p>In the second experimentation, the calibration constant, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x88.png" xlink:type="simple"/></inline-formula>, and also the switching threshold, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x88.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/9-1560137x89.png" xlink:type="simple"/></inline-formula>were</p><p>changed based on operator’s learning and comfort level. A lesser overshoot was achieved but at the cost of an increased settling time. It is noteworthy that the sensitivity of the joystick was also decreased in the second ex- perimentation. In both experiments the operator could guide the model helicopter to a safe zone blindfolded. In both figures, the dashed line represents the time at which the computer-based controller takes effect.</p></sec><sec id="s6"><title>6. Conclusions and Future Work</title><p>A methodology to incorporate tactile sensory information into the feedback control loop of a 2-dof model heli- copter was addressed. It was shown through some experiments that body-referenced touch sensors can provide the human operator with sufficient information to guide a free-flying object even in the case that visual clues are obstructed. A blindfolded operator, after going through a training period, could navigate the model helicopter to a safe zone. The overall control strategy was divided into two main stages, namely the human-in-the-loop and the computer-based control. An LQR-based control strategy would suffice to generate proper control cues to the operator, subjected to a cluster of vibro-tactile motors mounted on his/her body. A computer-based control takes effect at the end of motion for precise homing.</p><p>Further characterization of the human learning curve for the more robust and smooth control of our model helicopter is under investigation. Incorporating visual clues when controlling the model helicopter in awkward maneuvers such as those found in stunt flying in form of a sensor fusion technique will be also investigated. In- corporating collision avoidance strategies into the control loop through tactile sensory information on a free fly- ing UAV is envisioned in a long term.</p></sec></body><back><ref-list><title>References</title><ref id="scirp.51808-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">McLean, D. (2006) Automatic Flight Control Systems. Measurement and Control, 36, 172-177. 
http://dx.doi.org/10.1177/002029400303600602</mixed-citation></ref><ref id="scirp.51808-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Yokokohji, Y., Tubouchi, T., Tanaka, A. and Yoshida, T. (2006) Guidelines for Human Interface Design of Rescue Robots. Proceedings of the SICE-ICASE IEEE International Conference, Busan, 18-21 October 2006, 3455-3460.</mixed-citation></ref><ref id="scirp.51808-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Lam, T.M., Mulder, M. and Paassen, V. (2007) Haptic Interface for UAV Collision Avoidance. International Journal of Aviation Psychology, 17, 167-195. http://dx.doi.org/10.1080/10508410701328649</mixed-citation></ref><ref id="scirp.51808-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Gwo-Ruey, Y. (2007) Robust-Optimal Control of a Nonlinear Two Degree-of-Freedom Helicopter. Proceedings of the 6th IEEE/ACIS International Conference on Computer and Information Science, Melbourne, 11-13 July 2007, 744-749.</mixed-citation></ref><ref id="scirp.51808-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Quanser Company (2014) Quanser 2 DOF Helicopter User Manual. 
http://users.encs.concordia.ca/~realtime/elec483/docs/2%20DOF%20Helicopter%20User%20Manual.pdf</mixed-citation></ref><ref id="scirp.51808-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Mathwork (2014) Simulink Response Optimization User’s Guide. 
http://www.mathworks.com/help/sldo/index.html#response-optimization</mixed-citation></ref><ref id="scirp.51808-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Quanser (2014) Q8 WinCon Support &amp; Installation Guide. http://www.quanser.com/Products/q8</mixed-citation></ref></ref-list></back></article>