<?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">JFCMV</journal-id><journal-title-group><journal-title>Journal of Flow Control, Measurement &amp; Visualization</journal-title></journal-title-group><issn pub-type="epub">2329-3322</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jfcmv.2014.24016</article-id><article-id pub-id-type="publisher-id">JFCMV-50486</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Biomedical&amp;Life Sciences</subject><subject> Engineering</subject></subj-group></article-categories><title-group><article-title>
 
 
  Precise Positioning of Pneumatic Artificial Muscle Systems
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>haofei</surname><given-names>Wang</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Kaiji</surname><given-names>Sato</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>Toshiharu</surname><given-names>Kagawa</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Interdisciplinary Graduate School of Science and Engineering, Tokyo Institute of Technology, Yokohama, Japan</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>sato.k.ah@m.titech.ac.jp(KS)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>18</day><month>09</month><year>2014</year></pub-date><volume>02</volume><issue>04</issue><fpage>138</fpage><lpage>153</lpage><history><date date-type="received"><day>7</day>	<month>June</month>	<year>2014</year></date><date date-type="rev-recd"><day>12</day>	<month>September</month>	<year>2014</year>	</date><date date-type="accepted"><day>22</day>	<month>September</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>
 
 
  Pneumatic artificial muscles (PAMs) currently possess a high power-to-weight ratio, a high power-to-volume ratio, and a high degree of safety. They have therefore been applied to many power assist devices and positioning mechanisms such as bionic robots, welfare devices, and parallel manipulators. However, the significant nonlinear characteristics of PAM mechanisms limit their positioning accuracies. The accuracies are generally lower than 5 μm, which preclude the PAM from precision systems. Nevertheless, enhancing the positioning accuracy is desired to extend the application fields of PAMs. This study aims to clarify a practical controller design method to achieve the precise positioning of PAM systems. As the first step of this research, a linear motion mechanism with a pair of McKibben PAMs was constructed and a conventional dynamic model for this system is introduced. The dynamic model is used to explain the basic characteristics of the PAM mechanism and discuss the necessary characteristics for precise positioning. Then open-loop step and sinusoidal responses of the PAM mechanism were examined by experimental and simulated results. Next, for precise positioning, the practical controller design procedure is discussed and determined based on the measured open-loop responses. The proposed controller design procedure can be easily implemented into PAM mechanisms without an exact dynamic model. The positioning performance of such a system was experimentally evaluated. The experimental results show that although the positioning accuracy depends on the target position, the positioning error is lower than 1 μm even in the worst case and the positioning resolution can be set to 0.5 μm.
 
</p></abstract><kwd-group><kwd>Precision</kwd><kwd> Positioning</kwd><kwd> McKibben Pneumatic Muscle</kwd><kwd> Pneumatic Artificial Muscle</kwd><kwd> Accuracy</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>McKibben pneumatic artificial muscles, which are generally called pneumatic artificial muscles (PAMs), consist of an inflatable rubber tube sheathed by a braided mesh. It generates a pulling force via pressurized air and contracts in the axial direction while expanding in the radial direction. PAMs possess many advantages over traditional motor actuators, such as a high power-to-weight ratio, a high power-to-volume ratio, a high degree of safety [<xref ref-type="bibr" rid="scirp.50486-ref1">1</xref>] , and stick-slip-free operation [<xref ref-type="bibr" rid="scirp.50486-ref2">2</xref>] . They have therefore been applied in bionic robots [<xref ref-type="bibr" rid="scirp.50486-ref3">3</xref>] [<xref ref-type="bibr" rid="scirp.50486-ref4">4</xref>] , welfare devices [<xref ref-type="bibr" rid="scirp.50486-ref5">5</xref>] -[<xref ref-type="bibr" rid="scirp.50486-ref7">7</xref>] , and parallel manipulators [<xref ref-type="bibr" rid="scirp.50486-ref8">8</xref>] [<xref ref-type="bibr" rid="scirp.50486-ref9">9</xref>] . In particular, PAMs attract a great deal of attention in welfare devices. Owing to its high safety, low weight, and powerful output, it has been applied as a power assist device for rehabilitation applications such as an upper limb exoskeletal rehabilitation robot [<xref ref-type="bibr" rid="scirp.50486-ref5">5</xref>] and a knee-ankle-foot orthosis [<xref ref-type="bibr" rid="scirp.50486-ref6">6</xref>] .</p><p>However, PAMs have significant nonlinearity [<xref ref-type="bibr" rid="scirp.50486-ref10">10</xref>] , creep phenomenon [<xref ref-type="bibr" rid="scirp.50486-ref11">11</xref>] , and hysteresis [<xref ref-type="bibr" rid="scirp.50486-ref12">12</xref>] . These provide the low controllability and make it difficult to move finely. Thus precise motion control of PAMs is an important and unsolved problem. Precise motion control of the PAM provides the potential to extend its application owing to its many advantages. For example, the precise motion control of PAMs can make precise assembly robots powerful and lightweight and precise assembly robots are expected to become safer for the assemblers who work with the robots. In addition, the precise motion of PAMs enables nursing care robots to provide disabled people with delicate support such as injections and repairing clothes. A neurosurgical robot using hydraulic and pneumatic cylinders for precise positioning has been reported in [<xref ref-type="bibr" rid="scirp.50486-ref13">13</xref>] . Since PAMs have superior characteristics compared to hydraulic and pneumatic cylinders in the application of welfare devices, the precise motion of the PAM may be may also be used for a neurosurgical robot in the future.</p><p>In recent years, a variety of control methods have been developed to achieve excellent control performance in PAM mechanisms. These methods can be divided into three categories: conventional Proportional Integral Derivative (PID) control, nonlinear model-based control, and intelligent control. Conventional PID controllers are typically easy to design. However, they cannot sufficiently compensate for the nonlinearity of PAMs and lead to poor accuracy. The positioning errors by PID controllers are on the level of 0.1˚ for a rotational PAM mechanism [<xref ref-type="bibr" rid="scirp.50486-ref14">14</xref>] [<xref ref-type="bibr" rid="scirp.50486-ref15">15</xref>] and 200 &#181;m for a linear PAM mechanism [<xref ref-type="bibr" rid="scirp.50486-ref16">16</xref>] .</p><p>The nonlinear model-based control is another major approach for PAM mechanisms. Tondu and Lopez proposed a sliding mode controller using a nonlinear length-tension model of PAMs [<xref ref-type="bibr" rid="scirp.50486-ref17">17</xref>] . Even though some modified nonlinear model-based controllers with different models were subsequently presented, the positioning accuracy was on order of a micrometer because of inaccurate modeling [<xref ref-type="bibr" rid="scirp.50486-ref18">18</xref>] -[<xref ref-type="bibr" rid="scirp.50486-ref20">20</xref>] . Li and Kawashima proposed a nonlinear dynamic model where the force generated by a PAM is a function of the operating pressure and contraction length of the PAM [<xref ref-type="bibr" rid="scirp.50486-ref21">21</xref>] . The model also contains a viscous damping coefficient. They applied this nonlinear model to their PAM control system and reported that the highest positioning accuracy was 6 &#181;m.</p><p>Although the nonlinear model-based control provides a positioning accuracy that is higher than PID control for PAM mechanisms, its design procedures require the modeling of mechanisms and a sufficient knowledge of control theory. This makes it difficult to design for engineers who are unfamiliar with the controller design and makes the controller impracticable for real-world applications.</p><p>Intelligent control has also been widely used for PAM mechanisms [<xref ref-type="bibr" rid="scirp.50486-ref22">22</xref>] -[<xref ref-type="bibr" rid="scirp.50486-ref25">25</xref>] . The most accurate positioning result by this approach was reported in [<xref ref-type="bibr" rid="scirp.50486-ref25">25</xref>] where the researchers proposed an intelligent switching control method to adjust the gains of the PID controller using neural networks and achieved a rotational positioning error of 0.05˚. The problems of imprecise positioning and complex controller design still exist in intelligent control for PAM mechanisms due to the necessity of modeling mechanisms and sufficient knowledge of the intelligent algorithm.</p><p>This research aims to provide a practical controller design method for the precise motion for PAM mechanisms. In this paper, a practical and effective controller for the precise positioning of a PAM mechanism with a pair of PAMs is proposed and its usefulness is demonstrated experimentally. Characteristics of the PAM mechanism are examined by simulation using a conventional dynamic model and by experiment using an experimental setup. Comparison results between the simulated and experimental results are used to discuss the usefulness of the dynamic model in macroscopic analysis and the problems in controller design for precise positioning. Based on the discussion, we propose a practical controller that can be easily designed for precise positioning of the PAM mechanism, and the positioning performances are evaluated experimentally. In contrast to the nonlinear model-based controllers and the intelligent controllers, the proposed controller achieved the sub-micro- meter order positioning accuracy.</p></sec><sec id="s2"><title>2. Characteristics of the PAM Mechanism</title><sec id="s2_1"><title>2.1. Experimental Setup</title><p><xref ref-type="fig" rid="fig1">Figure 1</xref> shows the experimental setup in this research. An antagonistic structure PAM mechanism is used as a</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Experimental setup</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2760053x5.png"/></fig><p>controlled mechanism. This structure is widely used as the basic actuating element in many applications of PAMs [<xref ref-type="bibr" rid="scirp.50486-ref3">3</xref>] [<xref ref-type="bibr" rid="scirp.50486-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.50486-ref7">7</xref>] . The two PAMs (FESTO DMSP-10-100N-RM-CM) generate pulling forces to push and pull the mover (1.75 kg) along the moving direction in a working range of −1 mm to 1 mm. Air is injected from the pressure supply (0.5 MPa, absolute) and controlled by a 5-port 3-way proportional servo valve (FESTO MPYE- 5-1/8LF-010-B). The pressures in the two PAMs are measured using two pressure sensors (SMC PSE540A-01, resolution: 0.0012 MPa), which are only used for observing the behavior of the mechanism. A linear encoder (GSI Mercury II 5800, resolution: 100 nm) is used as a single feedback sensor in this research. The data pro- cessing unit is used as a controller hardware. The controller is implemented at a sampling rate of 10 kHz.</p></sec><sec id="s2_2"><title>2.2. Conventional Dynamic Model</title><p><xref ref-type="fig" rid="fig2">Figure 2</xref> shows a schematic diagram of a conventional dynamic model of PAM system. This dynamic model is used to examine basic characteristics of the PAM mechanism and discuss essential characteristics for controller design of precise positioning.</p><p>In <xref ref-type="fig" rid="fig2">Figure 2</xref>, the input voltage of servo valve <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x6.png" xlink:type="simple"/></inline-formula> determines the effective opening area of the ports which connect to PAM 1, PAM 2 (<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x7.png" xlink:type="simple"/></inline-formula>and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x8.png" xlink:type="simple"/></inline-formula>) and atmosphere (<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x9.png" xlink:type="simple"/></inline-formula>and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x10.png" xlink:type="simple"/></inline-formula>). The values of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x11.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x12.png" xlink:type="simple"/></inline-formula> are refer to [<xref ref-type="bibr" rid="scirp.50486-ref26">26</xref>] and the values of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x13.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x14.png" xlink:type="simple"/></inline-formula> are decided by fitting the step responses of simulated results to experimental results. The flow dynamic block shown in <xref ref-type="fig" rid="fig2">Figure 2</xref> is modeled as the process of a flow of an ideal gas through a converging nozzle [<xref ref-type="bibr" rid="scirp.50486-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.50486-ref7">7</xref>] [<xref ref-type="bibr" rid="scirp.50486-ref27">27</xref>] . This process is expressed as a nonlinear function for calculation of mass flow rate. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x15.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x16.png" xlink:type="simple"/></inline-formula> are the mass flow rates when air flowing into PAM 1 and PAM 2 respectively; <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x17.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x18.png" xlink:type="simple"/></inline-formula> are the mass flow rates when air leaking out from PAMs. In this process, the air flow is divided into two regimes which are choked and unchoked flow through the orifice. Based on the differentiation of state equation of ideal gas and the volume estimation [<xref ref-type="bibr" rid="scirp.50486-ref27">27</xref>] - [<xref ref-type="bibr" rid="scirp.50486-ref29">29</xref>] , the change rate of pressure in each PAM <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x19.png" xlink:type="simple"/></inline-formula> can be calculated [<xref ref-type="bibr" rid="scirp.50486-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.50486-ref7">7</xref>] [<xref ref-type="bibr" rid="scirp.50486-ref27">27</xref>] [<xref ref-type="bibr" rid="scirp.50486-ref29">29</xref>] .</p><p>The force generated by a single PAM can be described as a nonlinear function of the contraction length <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x20.png" xlink:type="simple"/></inline-formula> and the air pressure <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x21.png" xlink:type="simple"/></inline-formula> in the PAM [<xref ref-type="bibr" rid="scirp.50486-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.50486-ref27">27</xref>] [<xref ref-type="bibr" rid="scirp.50486-ref29">29</xref>] . In this research, the nonlinear force is identified as the expression shown as following [<xref ref-type="bibr" rid="scirp.50486-ref21">21</xref>] .</p><disp-formula id="scirp.50486-formula177"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2760053x22.png"  xlink:type="simple"/></disp-formula><p>where the contraction length <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x23.png" xlink:type="simple"/></inline-formula> is the difference between the original length of PAM <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x24.png" xlink:type="simple"/></inline-formula> and the length after contraction<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x25.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x26.png" xlink:type="simple"/></inline-formula>is the atmosphere pressure, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x27.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x28.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x29.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x30.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x24.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x31.png" xlink:type="simple"/></inline-formula>are identified coefficients.</p><p>Since the PAMs are constructed in an opposing pair configuration in this research, the motion equation of the PAM mechanism can be described as following [<xref ref-type="bibr" rid="scirp.50486-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.50486-ref27">27</xref>] .</p><disp-formula id="scirp.50486-formula178"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2760053x32.png"  xlink:type="simple"/></disp-formula><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Schematic diagram of the dynamic model</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2760053x33.png"/></fig><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x34.png" xlink:type="simple"/></inline-formula> is the mass of the mover, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x35.png" xlink:type="simple"/></inline-formula>is the position of the mover (<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x36.png" xlink:type="simple"/></inline-formula>,<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x37.png" xlink:type="simple"/></inline-formula>),<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x38.png" xlink:type="simple"/></inline-formula>is the damping coefficient and is tuned as <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x39.png" xlink:type="simple"/></inline-formula> optimally. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x40.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x41.png" xlink:type="simple"/></inline-formula> are forces generated by two PAMs. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x41.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x42.png" xlink:type="simple"/></inline-formula>is the friction force which is represented by a dynamic friction force and a static friction force. The mover is considered to be at a stop in the simulation using the friction force model when the absolute velocity of the mover is smaller than 10<sup>−9</sup> m/s. In this case, the friction force acts as the static friction force shown in <xref ref-type="fig" rid="fig3">Figure 3</xref>(a) where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x38.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x41.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x43.png" xlink:type="simple"/></inline-formula> is the applied force to the mover. In the other cases, the friction force acts as the dynamic friction force shown in <xref ref-type="fig" rid="fig3">Figure 3</xref>(b).</p><p>From <xref ref-type="fig" rid="fig2">Figure 2</xref> and Equation (2), the general nonlinear dynamic model of the PAM mechanism from the input voltage <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x44.png" xlink:type="simple"/></inline-formula> to the mover position <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x44.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2760053x45.png" xlink:type="simple"/></inline-formula> can be considered as a 3<sup>rd</sup>-order system [<xref ref-type="bibr" rid="scirp.50486-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.50486-ref27">27</xref>] .</p></sec><sec id="s2_3"><title>2.3. Open-Loop Responses of the PAM Mechanism</title><p>In this section, simulated open-loop responses using the dynamic model and experimental open-loop responses using the actual PAM mechanism are compared. The model parameters were identified from the results of fundamental experiments. Accordingly, the validity of the dynamic model is confirmed. In addition, as necessary characteristics in the controller design for precise positioning, the open-loop step responses which is the basic action in positioning control and the open-loop sinusoidal responses which can be summarized as gain characteristics in control system are investigated.</p><sec id="s2_3_1"><title>2.3.1. Open-Loop Step Responses</title><p>The open-loop step responses of the PAM mechanism were investigated by experimental and simulated results. <xref ref-type="fig" rid="fig4">Figure 4</xref> shows the experimental and simulated time responses of the pressure difference between PAM 1 and PAM 2 (P<sub>1</sub> ‒ P<sub>2</sub>) and the mover position x to the stepwise increasing and decreasing input voltages u. Although a single PAM can contract about 25 mm, the working range of PAMs is smaller than +/−3 mm, because the two PAMs were assembled in the condition of no initial contraction and the experimental setup was constructed in a push-pull structure [<xref ref-type="bibr" rid="scirp.50486-ref30">30</xref>] . In this research, furthermore the pressure supply (compressor) and the servo valve decrease the working range to +/−1.5 mm. Moreover, according to the time responses in <xref ref-type="fig" rid="fig4">Figure 4</xref>, we can observe that the input voltage u out of the range from 1.2 V to −1.2 V saturates the pressure and the mover position.</p><p>From these compared results, it can be observed that the simulated results tally with the experimental ones. Besides, no vibration in transient responses can be found both in the simulated and the experimental results. This indicates that step input signals can give smooth and no vibration step responses of the open-loop PAM mechanism and use for rough positioning.</p><p>Based on the simulated and experimental results shown in <xref ref-type="fig" rid="fig4">Figure 4</xref>, the static relationship between the position x and the input voltage u and the relationship between the pressure difference (P<sub>1</sub> − P<sub>2</sub>) and the input voltage u are summarized in <xref ref-type="fig" rid="fig5">Figure 5</xref>(a) and <xref ref-type="fig" rid="fig5">Figure 5</xref>(b) respectively. Both the simulated and the experimental results represent strong nonlinearities, so they should be compensated and can be used as compensators in the controller design shown in Section 3.</p><fig-group id="fig3"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Friction force in the dynamic model. a) Dynamic friction; b) Static friction.</title></caption><fig id ="fig3_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2760053x47.png"/></fig><fig id ="fig3_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2760053x46.png"/></fig></fig-group><fig-group id="fig4"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Open-loop step responses of the PAM mechanism.</title></caption><fig id ="fig4_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2760053x48.png"/></fig></fig-group></sec><sec id="s2_3_2"><title>2.3.2. Sinusoidal Responses</title><p><xref ref-type="fig" rid="fig6">Figure 6</xref> shows the experimental and simulated time responses to a sinusoidal input (amplitude = 1 (V), frequency = 1 (Hz), offset = 0 (V)). It can be found that the simulated result broadly corresponds with the experimental result when the mover moves in a large range.</p><p>During precise positioning, the position of mover is adjusted finely and small vibration is reduced around each target position, so the input-output relationship of vibration in various target position should be clarified. <xref ref-type="fig" rid="fig7">Figure 7</xref>(a) shows the experimental time responses of the pressure difference (P<sub>1</sub> ‒ P<sub>2</sub>) and the position x to three sinusoidal voltage inputs u. The input offset values, frequencies, and amplitudes are shown in the figure. <xref ref-type="fig" rid="fig7">Figure 7</xref>(b) shows the zoomed time response of the pressure difference (P<sub>1</sub> ‒ P<sub>2</sub>) and the position x. In the experiment, the offset of the sinusoidal voltage inputs u was 0.5 V, the vibration frequency was 5 Hz and the amplitude was 0.01 V. According to this figure, the simulated result does not match the experimental result microscopically. The dynamic model cannot represent some fine motion of the PAM mechanism accurately. This indicates that it is hard to estimate positioning accuracy precisely by simulation using the dynamic model.</p><fig-group id="fig5"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Static relationships between input voltage, mover position, and pressure difference. (a) Static relationship between mover position and input voltage; (b) Static relationship between pressure difference and input voltage.</title></caption><fig id ="fig5_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2760053x49.png"/></fig></fig-group><fig id="fig6"  position="float"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> Open-loop time responses of the PAM mechanism to a sinusoidal input in simulation and experiment</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2760053x50.png"/></fig><fig-group id="fig7"><label><xref ref-type="fig" rid="fig7">Figure 7</xref></label><caption><title> Open-loop experimental sinusoidal responses of the PAM mechanism. (a) Open-loop experimental sinusoidal responses of three different input voltage; (b) Open-loop experimental and simulated sinusoidal responses.</title></caption><fig id ="fig7_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2760053x51.png"/></fig></fig-group><p>Moreover, using different sinusoidal inputs, similar open-loop experiments were performed to examine the vibration amplitude relationship between the position x and the pressure difference (P<sub>1</sub> ‒ P<sub>2</sub>). The summarized results are shown in <xref ref-type="fig" rid="fig8">Figure 8</xref>. This vibration amplitude relationship depends on the amplitude of the mover’s vibration rather than the mover position and the vibration frequency. The amplitude relationship can be approximated as the least square fitted line shown in <xref ref-type="fig" rid="fig8">Figure 8</xref> and described by Equation (3).</p><disp-formula id="scirp.50486-formula179"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2760053x52.png"  xlink:type="simple"/></disp-formula><p>A dead zone can be found in the relationship between the vibration amplitude of the mover position x and that of the pressure difference (P<sub>1</sub> ‒ P<sub>2</sub>), and this dead zone should be compensated in the controller design. The simulated results also are shown in <xref ref-type="fig" rid="fig8">Figure 8</xref>. Although the dead zone is similar, the slopes of the vibration amplitude relationship between the position x and the pressure difference (P<sub>1</sub> ‒ P<sub>2</sub>) are different between experimental and simulated results.</p><p>According to the comparison between experimental and simulated results, this dynamic model is useful macroscopically, but not microscopically. The actual microscopic responses are important in the controller design for precise positioning. However it is hard to use more complex dynamic model in the controller design. The dynamic model is already complex enough to hesitate to use for the controller design. Thus the design procedure without the dynamic model is desired.</p></sec></sec></sec><sec id="s3"><title>3. Controller Design for Precise Positioning</title><sec id="s3_1"><title>3.1. Design Concept</title><p>This study focuses on point-to-point precise positioning. A suitable controller and its design procedure are proposed for the precise positioning of PAM mechanisms.</p><p>As discussed in Section 2.3, the macroscopic characteristics of the dynamic model match the ones of actual PAM mechanism, but not the microscopic characteristics of the dynamic model, such as vibration characteristic in fine motion which is important in precise positioning. Furthermore, the PAM mechanism has significant nonlinear characteristics and it is difficult to represent the characteristics by a bode diagram. Besides, even though great effort and much time need to be made to identify an accurate dynamic model and make some adaptive elements for control, the controller parameters in some adaptive elements based on the model still should be tuned by trials and errors for good performance of positioning.</p><p>On the other hand, the PAM mechanism has characteristics which make it relatively easy to design the controller by trials and errors. It can be represented by a nonlinear 3<sup>rd</sup>-order system as discussed in Section 2.2 and has high damping characteristics macroscopically as shown in <xref ref-type="fig" rid="fig4">Figure 4</xref>. It is generally practical to adjust the parameters of a PID controller for the 3<sup>rd</sup>-order systems by trials and errors. Significant nonlinear characteristics often become obstructed the adjustment. However the static relationships in <xref ref-type="fig" rid="fig5">Figure 5</xref>(a) and the amplitude relationship in <xref ref-type="fig" rid="fig8">Figure 8</xref>, which are useful to compensate significant nonlinear characteristics, can be easily-mea- sured and used without dynamic models. Therefore, instead of the accurate model of the PAM mechanism, the</p><fig id="fig8"  position="float"><label><xref ref-type="fig" rid="fig8">Figure 8</xref></label><caption><title> Amplitude relationship between position and pressure dif- ference</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2760053x53.png"/></fig><p>measured input-output relationships are used for the controller design in this study. Meanwhile, this method can make the controller easy to design for engineers without a deep understanding of the knowledge of fluid mechanics and control engineering.</p><p>The positioning action of the proposed controller is divided into two phases: a reaching phase in the transient state and a following phase in the steady state. In the reaching phase, an open-loop controller is used to drive the mover to the target position. From the response in <xref ref-type="fig" rid="fig4">Figure 4</xref>, the controller can be expected to provide no or small overshoot responses. A closed-loop controller is active in the following phase and determines the precise positioning in the steady state. Thus the closed-loop controller is designed in the steady-state condition and the behavior in the transient state does not need to be accounted for in the controller design. The controller for the reaching phase is determined without the need to account for stability. The positioning accuracy basically depends on the characteristics of the control system in the following phase, but not the reaching phase. These features make it easy to design the controller.</p><p>According to the controller design concept shown above, the controller structure can be represented by <xref ref-type="fig" rid="fig9">Figure 9</xref>. The “Controller switch algorithm” switches between the open-loop controller and closed-loop controllers, and the flowchart shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>0 shows the concrete switching condition where T<sub>s</sub> is the sampling time and i is the computing step of the control system. During the positioning action in the reaching phase, the mover is driven using the open-loop controller. At this time, the input voltage to the PAM mechanism is chosen such that it is sufficient to drive the mover to the target position x<sub>r</sub>. In the following phase, the closed-loop controller works on reducing the steady-state error and the residual vibration.</p><fig id="fig9"  position="float"><label><xref ref-type="fig" rid="fig9">Figure 9</xref></label><caption><title> Basic controller structure</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2760053x54.png"/></fig><fig id="fig10"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>0</label><caption><title> Flowchart of controller switch algorithm</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2760053x55.png"/></fig></sec><sec id="s3_2"><title>3.2. Design Procedure</title><p>Based on the controller design concept introduced above, the controller is designed as follows.</p><p>Step 1: Decision of the open-loop controller</p><p>The open-loop controller is used to provide a voltage to drive the mover to the target position in the reaching phase. When the target position is known, the input voltage for driving the mover to the target position can be determined from the relationship shown in <xref ref-type="fig" rid="fig5">Figure 5</xref>(a). This corresponds to the “Feed-forward element S” in <xref ref-type="fig" rid="fig1">Figure 1</xref>1.</p><p>Step 2: Decision of the closed-loop controller</p><p>In the closed-loop controller, the “Feed-forward element S” provides an enough voltage which enables the mover to move to the target position and stay there roughly. It has been reported that a feed-forward element can be obtained from model inversion or estimated by an online numerical differentiator [<xref ref-type="bibr" rid="scirp.50486-ref31">31</xref>] [<xref ref-type="bibr" rid="scirp.50486-ref32">32</xref>] . Our current work focuses on achieving an easy controller design procedure for PAM mechanisms for engineers who are unfamiliar with the controller design. Thus the feed-forward element experimentally-measured by open-loop responses is adopted in the controller design.</p><p>1) The PID compensator design</p><p>The PD elements are calibrated to reduce the residual vibration and their outputs are connected to the input of “Compensator D” which shows the vibration amplitude characteristic. The integrator is calibrated to reduce the steady-state error. Unlike the outputs of PD elements, the output of integrator is a signal with DC component which is used to reduce the steady-state error rather than the amplitude of residual vibration. Therefore, the integrator output is not connected to the input of “Compensator D”. In addition to the calibration of K<sub>p</sub>, K<sub>i</sub> and K<sub>d</sub>, the parameter T<sub>d</sub> of the low-pass filter should be determined as well. In some cases, the low-pass filter deteriorates the dynamic characteristics. Accordingly, in order to suppress the negative effect of the low-pass filter, a cut-off frequency is selected for the low-pass filter f<sub>d</sub> = 300 Hz (f<sub>d</sub> = 1/T<sub>d</sub>) that is ten times the maximum value of residual vibration frequency f<sub>v</sub><sub>max</sub> (f<sub>v</sub><sub>max</sub> = 0.1 f<sub>d</sub>). Because the residual vibration often has a frequency near the edge of or outside the bandwidth of the control system, the designed low-pass filter is expected to exhibit little or no effect on the response of the control system.</p><p>2) Compensation of the dynamic nonlinearity</p><p>The mechanism has a nonlinear characteristic between the input voltage to the PAM mechanism and the mover displacement, which should be compensated for to allow precise positioning. The amplitude relationship shown in <xref ref-type="fig" rid="fig8">Figure 8</xref> is used as one of the control elements. <xref ref-type="fig" rid="fig1">Figure 1</xref>2 shows the control element for canceling the nonlinear characteristic, which corresponds to the “Compensator D” in <xref ref-type="fig" rid="fig1">Figure 1</xref>1. The relationship shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>2 is expressed in Equation (4). The C<sub>PV</sub> in <xref ref-type="fig" rid="fig1">Figure 1</xref>1 is the linearized gain from pressure to input voltage shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>3. The nonlinearity between the input and the output of closed-loop controller is compensated for by the “Compensator D” and C<sub>PV</sub>.</p><disp-formula id="scirp.50486-formula180"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2760053x56.png"  xlink:type="simple"/></disp-formula><p>Step 3: Fine-tuning of the acceleration feedback element</p><p>The acceleration feedback element is added to further reduce the residual vibration. Since the acceleration is essentially proportional to the pressure, the output of the acceleration feedback element is connected to the node of the pressure signal.</p><fig id="fig11"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>1</label><caption><title> Block diagram of control system for precise positioning of PAM mechanisms</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2760053x57.png"/></fig><fig id="fig12"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>2</label><caption><title> “Compensator D” obtained from the amplitude relationship between position and pre- ssure difference</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2760053x58.png"/></fig><fig id="fig13"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>3</label><caption><title> Linearized gain between the pres- sure and the input voltage C<sub>PV</sub></title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2760053x59.png"/></fig><p>In summary, the following measured relationships were determined to be applied in the controller design: the static relationship between the mover position and the valve voltage was used as “Feed-forward element S”, the static relationship between the pressure difference and the valve voltage was linearized as a gain of “C<sub>PV</sub>”. The vibration amplitude relationship between the mover position and the pressure difference (P<sub>1</sub> - P<sub>2</sub>) was applied to construct “Compensator D”. These compensators of the proposed control system can be determined from open- loop responses and the parameters can be tuned easily without any model parameters. Therefore, compared to the controllers in prior studies, the current controller can easily be implemented in PAM mechanisms in practical applications.</p></sec></sec><sec id="s4"><title>4. Positioning Performance</title><p>Positioning experiments were carried out to validate the usefulness of the proposed controller. The duration of every positioning experiment was 5 s and the step change begins at the 1 s mark.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>4 shows positioning responses of the proposed control system to a 1 mm step input. For comparison, the experimental and simulated responses of a conventional PID control system are also shown. In the PID controller design, the proportional gain and the integral gain were adjusted and the derivative gain was then tuned to reduce the residual vibration as much as possible. The gains used in the PID controller and proposed controller are shown in <xref ref-type="table" rid="table1">Table 1</xref>. A comparison of the results indicates that the proposed controller can reduce the residual vibration amplitude to 0.8 μm whereas the residual vibration amplitude for the PID controller was 2.8 μm. On the other hand, the simulated response of PID control system does not show the significant residual vibration and is not exact enough to discuss the precise positioning performance.</p><p>Moreover, in order to explain the effectiveness of “Compensator D”, the step response by the controller with a linearized “Compensator D” is shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>4 as well. The controller is the same as the result that in the proposed controller, but the “Compensator D” is replaced with a linearized “Compensator D”. It can be found that the vibration amplitude of the proposed control system in steady-state is smaller than that of the control system with a linearized “Compensator D”. The result indicates that the nonlinear “Compensator D” of the proposed controller is more effective to reduce the residual vibration in precise positioning.</p><p>To examine the influence of the nonlinearity of the PAM, the positioning performance was evaluated at different target positions. In the positive direction, positioning experiments with x<sub>r</sub> = 0.3 mm, 0.5 mm, 0.7 mm and 1 mm were performed and their results are shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>5. In the negative direction, positioning performance when x<sub>r</sub> = −0.5 mm and −1 mm were tested and the results are shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>6. The proposed controller was also evaluated when the target position was the original position x<sub>r</sub> = 0 mm and the result is shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>7. In the experiment, the initial position was 1 μm.</p><p>The experimental positioning results show that the positioning accuracy decreases with an increase of the distance between the target position and the original position. However the positioning errors are all less than 1 μm in the working range of the PAM mechanism. <xref ref-type="table" rid="table2">Table 2</xref> summarizes the averaged maximum steady-state errors and their standard deviations during 3 s to 5 s of 20 positioning experiments. According to the results in <xref ref-type="table" rid="table2">Table 2</xref>, sub-micrometer order positioning is achieved by the proposed controller. The positioning accuracy and the residual vibration frequency depend on the mover position although the reason is not clear. However, according to the positioning results, the resonant frequency changes with the mover position. The reason of this problem is not clear yet, but it will be clarified as the subject of next step of this research.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>8 shows the resolution test results in different positions. In these figures, <xref ref-type="fig" rid="fig1">Figure 1</xref>8(c) shows the resolution test result near x<sub>r</sub> = 0 mm where we can find that the vibration frequency of the mover is lower than those in Figures 18(a)-(b) <xref ref-type="fig" rid="fig1">Figure 1</xref>8(d) and <xref ref-type="fig" rid="fig1">Figure 1</xref>8(e), but similar to <xref ref-type="fig" rid="fig1">Figure 1</xref>7. According to these responses, a 0.5 μm resolution can be specified.</p><fig id="fig14"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>4</label><caption><title> A comparison of the positioning results of different controllers</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2760053x60.png"/></fig><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Controller gains of PID element in two control systems</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Gain</th><th align="center" valign="middle" >PID controller</th><th align="center" valign="middle" >Proposed controller</th></tr></thead><tr><td align="center" valign="middle" >K<sub>p</sub></td><td align="center" valign="middle" >1.2</td><td align="center" valign="middle" >4</td></tr><tr><td align="center" valign="middle" >K<sub>i</sub></td><td align="center" valign="middle" >30</td><td align="center" valign="middle" >40</td></tr><tr><td align="center" valign="middle" >K<sub>d</sub></td><td align="center" valign="middle" >0.005</td><td align="center" valign="middle" >0.15</td></tr></tbody></table></table-wrap><fig-group id="fig15"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>5</label><caption><title> Positioning results for positive target positions. (a) x<sub>r</sub> = 0.3 mm; (b) x<sub>r</sub> = 0.5 mm; (c) x<sub>r</sub> = 0.7 mm; (d) x<sub>r</sub> = 1 mm.</title></caption><fig id ="fig15_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2760053x61.png"/></fig><fig id ="fig15_2"><label>(c)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2760053x62.png"/></fig></fig-group><fig-group id="fig16"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>6</label><caption><title> Positioning results for negative target positions. (a) x<sub>r</sub> = −0.5 mm; (b) x<sub>r</sub> = −1 mm.</title></caption><fig id ="fig16_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2760053x63.png"/></fig></fig-group></sec><sec id="s5"><title>5. Conclusions</title><p>The practical and effective controller and its design procedure for precise positioning of PAM mechanisms have been proposed and its usefulness has been evaluated experimentally. PAM systems have nonlinear characteris-</p><fig id="fig17"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref>7</label><caption><title> Positioning results for target position x<sub>r</sub> = 0 mm</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2760053x64.png"/></fig><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Average and standard deviation of maximum steady-state errors in 20 positioning experiments</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Target position (mm)</th><th align="center" valign="middle" >Average (μm)</th><th align="center" valign="middle" >Standard deviation (μm)</th></tr></thead><tr><td align="center" valign="middle" >−1</td><td align="center" valign="middle" >0.880</td><td align="center" valign="middle" >0.062</td></tr><tr><td align="center" valign="middle" >−0.5</td><td align="center" valign="middle" >0.625</td><td align="center" valign="middle" >0.043</td></tr><tr><td align="center" valign="middle" >0</td><td align="center" valign="middle" >0.655</td><td align="center" valign="middle" >0.059</td></tr><tr><td align="center" valign="middle" >0.3</td><td align="center" valign="middle" >0.485</td><td align="center" valign="middle" >0.035</td></tr><tr><td align="center" valign="middle" >0.5</td><td align="center" valign="middle" >0.640</td><td align="center" valign="middle" >0.049</td></tr><tr><td align="center" valign="middle" >0.7</td><td align="center" valign="middle" >0.800</td><td align="center" valign="middle" >0.032</td></tr><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0.850</td><td align="center" valign="middle" >0.045</td></tr></tbody></table></table-wrap><p>tics such as a friction force and a nonlinear air flow dynamic. These characteristics provide the low controllability and make them difficult to move finely even if a high precision sensor is used.</p><p>The introduced conventional nonlinear dynamic model cannot show the same responses as the experimental mechanism microscopically although it provides the reasonable responses macroscopically. The practical controller design procedure needs three measured elements but not the model identification. The experimentally-measured elements were used in the controller design directly. This feature makes the controller design easy for engineers who are unfamiliar with control engineering. The positioning performance was evaluated experimentally. Experimental results demonstrated the effectiveness of the proposed controller with sub-micrometer positioning accuracy and 0.5 μm positioning resolution of the PAM mechanism. Therefore, it can be concluded that the proposed controller with a few nonlinear elements can be easily designed and provides high-precision positioning of PAM mechanisms. According to the positioning results, the positioning accuracy and the residual vibration frequency depend on the mover position. The reason of this problem is not clear yet, but it will be examined as a future work of this research.</p><p>The sampling rate of 10 kHz was used because the controller structure in this study was simple and good controller hardware was available. In this paper, we did not focus the reduction of the sampling rate. However the large difference between the sampling rate and the frequency of the residual vibration suggests the possibility to reduce it and make lower cost and lower-power controller hardwares useful. Based on the proposed controller design procedure, even though the change of the characteristics makes re-measurement of the input-output relationships necessary for good positioning performance, it is easy to measure the relationships. Furthermore, the compensation of the nonlinear characteristics enables the PD gains to increase and the increase can lead the robustness to the disturbance force. Therefore, the proposed controller design procedure is really reliable on industrial point of view and makes the total positioning system robust.</p><disp-formula id="scirp.50486-formula181"><graphic  xlink:href="http://html.scirp.org/file/2-2760053x65.png"  xlink:type="simple"/></disp-formula><p>(a) (b)</p><disp-formula id="scirp.50486-formula182"><graphic  xlink:href="http://html.scirp.org/file/2-2760053x67.png"  xlink:type="simple"/></disp-formula><p>(c) (d)</p><disp-formula id="scirp.50486-formula183"><graphic  xlink:href="http://html.scirp.org/file/2-2760053x69.png"  xlink:type="simple"/></disp-formula><p>(e)</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>8. Resolution test results for different positions. (a) Resolution test results near 0.5 mm; (b) Resolution test results near 1 mm; (c) Resolution test results near 0 mm; (d) Resolution test results near −0.5 mm; (e) Resolution test results near −1 mm.</p><p>Although this paper focused on precise positioning, precise tracking control is important in this research. Thus an important future work is to examine a precise tracking control method and its design procedure for PAM systems based on the proposed controller and achieve high tracking performance of PAM systems.</p></sec></body><back><ref-list><title>References</title><ref id="scirp.50486-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Reynolds, D.B., Repperger, D.W., Phillips, C.A. and Bandry, G. (2003) Modeling the Dynamic Characteristics of Pneumatic Muscle. Annals of Biomedical Engineering, 31, 310-317.&lt;br /&gt; http://dx.doi.org/10.1114/1.1554921</mixed-citation></ref><ref id="scirp.50486-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Jouppila, V., Gadsden, S.A. and Ellman, A. (2014) Experimental Comparisons of Sliding Mode Controlled Pneumatic Muscle and Cylinder Actuators. Journal of Dynamic Systems, Measurement, and Control, 136, Article ID: 044503. 
http://dx.doi.org/10.1115/1.4026873</mixed-citation></ref><ref id="scirp.50486-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Hosoda, K., Takuma, T., Nakamoto, A. and Hayashi, S. (2008) Biped Robot Design Powered by Antagonistic Pneumatic Actuators for Multi-Modal Locomotion. Robotics and Autonomous Systems, 56, 46-53. 
http://dx.doi.org/10.1016/j.robot.2007.09.010</mixed-citation></ref><ref id="scirp.50486-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Narioka, K. and Hosoda, K. (2011) Motor Development of a Pneumatic Musculoskeletal Infant Robot. Proceedings of the 2011 IEEE International Conference on Robotics and Automation, Shanghai, 9-13 May 2011, 963-968. 
http://dx.doi.org/10.1109/icra.2011.5980416</mixed-citation></ref><ref id="scirp.50486-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Balasubramanian, S., Ward, J., Sugar, T. and He, J. (2007) Characterization of the Dynamic Properties of PneumaticMuscle Actuators. Proceedings of the 2007 IEEE 10th International Conference on Rehabilitation Robotics, Noordwijk, 12-15 June 2007, 764-770. http://dx.doi.org/10.1109/icorr.2007.4428511</mixed-citation></ref><ref id="scirp.50486-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Wong, Z., Teng, C. and Chong, Y. (2012) Power Assisted Pnumatic-Based Knee-Ankle-Foot-Orthosis for Rehabilitation. Proceedings of the 2012 IEEE EMBS Conference on Biomedical Engineering and Sciences, Langkawi, 17-19 December 2012, 300-304. http://dx.doi.org/10.1109/iecbes.2012.6498058</mixed-citation></ref><ref id="scirp.50486-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Hussain, S., Xie, S.Q. and Jamwal, P.K. (2013) Robust Nonlinear Control of an Intrinsically Compliant Robotic Gait Training Orthosis. IEEE Transactions on Systems, Man, and Cybernetics: Systems, 43, 655-665. 
http://dx.doi.org/10.1109/tsmca.2012.2207111</mixed-citation></ref><ref id="scirp.50486-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Zhu, X., Tao, G., Yao, B. and Cao, J. (2008) Adaptive Robust Posture Control of a Parallel Manipulator Driven by Pneumatic Muscles. Automatica, 44, 2248-2257. &lt;br /&gt;http://dx.doi.org/10.1016/j.automatica.2008.01.015</mixed-citation></ref><ref id="scirp.50486-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Zhu, X., Tao, G., Yao, B. and Cao, J. (2009) Integrated Direct/Indirect Adaptive Robust Posture Trajectory Tracking Control of a Parallel Manipulator Driven by Pneumatic Muscles. IEEE Transactions on Control Systems Technology, 17, 576-588. http://dx.doi.org/10.1109/tcst.2008.2001715</mixed-citation></ref><ref id="scirp.50486-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Wickramatunge, K.C. and Leephakpreeda, T. (2010) Study on Mechanical Behaviors of Pneumatic Artificial Muscle. International Journal of Engineering Science, 48, 188-198. &lt;br /&gt;http://dx.doi.org/10.1016/j.ijengsci.2009.08.001</mixed-citation></ref><ref id="scirp.50486-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Minh, T.V., Kamers, B., Ramon, H. and Brussel, H.V. (2012) Modeling and Control of a Pneumatic Artificial Muscle Manipulator Joint-Part I: Modeling of a Pneumatic Artificial Muscle Manipulator Joint with Accounting for Creep Effect. Mechatronics, 22, 923-933. http://dx.doi.org/10.1016/j.mechatronics.2012.06.002</mixed-citation></ref><ref id="scirp.50486-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Vo-Minh, T., Tjahjowidodo, T., Ramon, H. and Brussel, H.V. (2011) A New Approach to Modeling Hysteresis in a Pneumatic Artificial Muscle Using the Maxwell-Slip Model. IEEE/ASME Transactions on Mechatronics, 16, 177-186. 
http://dx.doi.org/10.1109/tmech.2009.2038373</mixed-citation></ref><ref id="scirp.50486-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Raoufi, C., Goldenberg, A.A. and Kucharczyk, W. (2008) A New Hydraulically/Pneumatically Actuated MR-Compatible Robot for MRI-Guided Neurosurgery. Proceedings of the 2nd International Conference on Bioinformatics and Biomedical Engineering, Shanghai, 16-18 May 2008, 2232-3335. &lt;br /&gt;http://dx.doi.org/10.1109/icbbe.2008.889</mixed-citation></ref><ref id="scirp.50486-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Pujana-Arrese, A., Mendizabal, A., Arenas, J., Prestamero, R. and Landaluze, J. (2010) Modelling in Modelica and Position Control of a 1-DoF Set-Up Poweredby Pneumatic Muscles. Mechatronics, 20, 535-552. 
http://dx.doi.org/10.1016/j.mechatronics.2010.05.002</mixed-citation></ref><ref id="scirp.50486-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Li, X., He, F., Xia, H. and Ting, G. (2012) Implicit Generalized Predictive Control of Hip-Joint Rehabilitation Training Device Driven by Pneumatic Muscle Actuator. Applied Mechanics and Materials, 138-139, 273-278. 
http://dx.doi.org/10.4028/www.scientific.net/amm.138-139.273</mixed-citation></ref><ref id="scirp.50486-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Andrikopoulos, G., Nikolakopoulos, G. and Manesis, S. (2013) Non-Linear Control of Pneumatic Artificial Muscles. Proceeding of the 21st Mediterranean Conference on Control &amp; Automation, Chania, 25-28 June 2013, 729-734.  
http://dx.doi.org/10.1109/med.2013.6608804</mixed-citation></ref><ref id="scirp.50486-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Tondu, B. and Lopez, P. (2000) Modeling and Control of McKibben Artificial Muscle Robot Actuators. IEEE Control Systems, 20, 15-38. http://dx.doi.org/10.1109/37.833638</mixed-citation></ref><ref id="scirp.50486-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Xing, K., Huang, J., Wang, Y., Wu, J., Xu, Q. and He, J. (2010) Tracking Control of Pneumatic Artificial Muscle Actuators Based on Sliding Mode and Non-linear Disturbance Observer. Control Theory &amp; Applications, 4, 2058-2070. 
http://dx.doi.org/10.1049/iet-cta.2009.0555</mixed-citation></ref><ref id="scirp.50486-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Yang, L. and Lilly, H.L. (2003) Sliding Mode Tracking for Pneumatic Muscle Actuators in Bicep/Tricep Pair Configuration. Proceedings of the 2003 American Control Conference, Denver, 4-6 June 2003, 4669-4674. 
http://dx.doi.org/10.1109/acc.2003.1242460</mixed-citation></ref><ref id="scirp.50486-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Choi, T. and Lee, J. (2010) Control of Manipulator Using Pneumatic Muscles for Enhanced Safety. IEEE Transactions on Industrial Electronics, 57, 2815-2825. http://dx.doi.org/10.1109/tie.2009.2036632</mixed-citation></ref><ref id="scirp.50486-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Li, H., Kawashima, K., Tadano, K., Ganguly, S. and Nakano, S. (2013) Achieving Haptic Perception in Forceps’ Manipulator Using Pneumatic Artificial Muscle. IEEE/ASME Transactions on Mechatronics, 18, 74-85. 
http://dx.doi.org/10.1109/tmech.2011.2163415</mixed-citation></ref><ref id="scirp.50486-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Lilly, J.H. (2003) Adaptive Tracking for Pneumatic Muscle Actuators in Bicep and Tricep Configurations. IEEE Transactions on Neural Systems and Rehabilitation Engineering, 11, 333-339. &lt;br /&gt; 
http://dx.doi.org/10.1109/TNSRE.2003.816870</mixed-citation></ref><ref id="scirp.50486-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Chang, M., Liou, J. and Chen, M. (2011) T-S Fuzzy Model-Based Tracking Control of a One-Dimensional Manipulator Actuated by Pneumatic Artificial Muscles. Control Engineering Practice, 19, 1442-1449. 
http://dx.doi.org/10.1016/j.conengprac.2011.08.002</mixed-citation></ref><ref id="scirp.50486-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Leephakpreeda, T. (2011) Fuzzy Logic Based PWM Control and Neural Controlled-Variable Estimation of Pneumatic Artificial Muscle Actuators. Expert Systems with Applications, 38, 7837-7850.&lt;br /&gt; 
http://dx.doi.org/10.1016/j.eswa.2010.12.120</mixed-citation></ref><ref id="scirp.50486-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Ahn, K.K., Thanh, T.D.C. and Ahn, Y.K. (2005) Intelligent Switching Control of Pneumatic Artificial Muscle Manipulator. JSME International Journal Series C Mechanical Systems, Machine Elements and Manufacturing, 48, 657-667.  
http://dx.doi.org/10.1299/jsmec.48.657</mixed-citation></ref><ref id="scirp.50486-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Sato, K. and Sano, Y. (2014) Practical and Intuitive Controller Design Method for Precision Positioning of a Pneumatic Cylinder Actuator Stage. Precision Engineering, 38, 703-710. &lt;br /&gt;http://dx.doi.org/10.1016/j.precisioneng.2014.03.006</mixed-citation></ref><ref id="scirp.50486-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">Shen, X. (2010) Nonlinear Model-Based Control of Pneumatic Artificial Muscle Servo Systems. Control Engineering Practice, 18, 311-317. http://dx.doi.org/10.1016/j.conengprac.2009.11.010</mixed-citation></ref><ref id="scirp.50486-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">Chou, C. and Hannaford, B. (1996) Measurement and Modeling of McKibben Pneumatic Artificial Muscles. IEEE Transactions on Robotics and Automation, 12, 90-102. http://dx.doi.org/10.1109/70.481753</mixed-citation></ref><ref id="scirp.50486-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">Ganguly, S., Garg, A., Pasricha, A. and Dwivedy, S.K. (2012) Control of Pneumatic Artificial Muscle System through Experimental Modelling. Mechatronics, 22, 1135-1147. &lt;br /&gt;http://dx.doi.org/10.1016/j.mechatronics.2012.09.010</mixed-citation></ref><ref id="scirp.50486-ref30"><label>30</label><mixed-citation publication-type="other" xlink:type="simple">Festo (2010) Fluidic Muscle DMSP/MAS. Festo Brochure. &lt;br /&gt;   
https://www.festo.com/cat/ja_jp/data/doc_engb/PDF/EN/DMSP-MAS_EN.PDF</mixed-citation></ref><ref id="scirp.50486-ref31"><label>31</label><mixed-citation publication-type="other" xlink:type="simple">Fliess, M., Levine, J., Martin, P. and Rouchon, P. (1999) A Lie-Backlund Approach to Equivalence and Flatness of Nonlinear Systems. IEEE Transactions on Automatic Control, 44, 922-937. &lt;br /&gt;http://dx.doi.org/10.1109/9.763209</mixed-citation></ref><ref id="scirp.50486-ref32"><label>32</label><mixed-citation publication-type="other" xlink:type="simple">Fliess, M. and Join, C. (2008) Intelligent PID Controllers. 16th Mediterranean Conference on Control and Automation, Ajaccio, 25-27 June 2008, 326-331.</mixed-citation></ref></ref-list></back></article>