<?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">ENG</journal-id><journal-title-group><journal-title>Engineering</journal-title></journal-title-group><issn pub-type="epub">1947-3931</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/eng.2014.62013</article-id><article-id pub-id-type="publisher-id">ENG-42938</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Engineering</subject></subj-group></article-categories><title-group><article-title>
 
 
  Turn Control of a Three-Dimensional Quasi-Passive Walking Robot by Utilizing a Mechanical Oscillator
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>ing</surname><given-names>Cao</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>Soichiro</surname><given-names>Suzuki</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Yohei</surname><given-names>Hoshino</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Department of Mechanical Engineering, Kitami Institute of Technology, Kitami, Japan</addr-line></aff><aff id="aff1"><addr-line>Graduate School of Engineering, Kitami Institute of Technology, Kitami, Japan</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>d1171508033@std.kitami-it.ac.jp(YH)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>28</day><month>01</month><year>2014</year></pub-date><volume>06</volume><issue>02</issue><fpage>93</fpage><lpage>99</lpage><history><date date-type="received"><day>May</day>	<month>17,</month>	<year>2013</year></date><date date-type="rev-recd"><day>June</day>	<month>17,</month>	<year>2013</year>	</date><date date-type="accepted"><day>June</day>	<month>25,</month>	<year>2013</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>
 
 
   A turn control strategy is proposed in order to improve environmental adaptability of a quasi-passive walking robot by utilizing a mechanical oscillator. The target trajectory of the fmechanical oscillator is determined by online planning of its period, phase, amplitude and angle of the central axis of oscillation. The motion of the mechanical oscillator is always entrained with the rocking motion of the robot based on forced entrainment in order to stabilize the robot. The turn radius can be controlled by adjusting the inclination angle of the central axis of the mechanical oscillator movement, and the control method is numerically and experimentally examined. Results show that the robot can turn with different radius and it is possible for the robot to walk in various environments. Finally, the gait of turn is compared with that of straight walking and analyzed in terms of mechanical work and energy. 
 
</p></abstract><kwd-group><kwd>Passive Walking; Turn Control; Mechanical Oscillator; Forced Entrainment</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Passive walking robots can walk on a shallow slope by utilizing gravity and inertia without any actuator or control [<xref ref-type="bibr" rid="scirp.42938-ref1">1</xref>]. Gait of passive walking robots can be more natural than gait of humanoid robots, and their energy-efficiency is also higher [<xref ref-type="bibr" rid="scirp.42938-ref2">2</xref>]. Research of passive walking not only contributes to the understanding of the mechanism of biped walking but also helps to improve the design and control of biped robots.</p><p>However, it is difficult to stabilize passive walking robots in different environments because passive walking is sensitive to variations of initial condition and slope angle [<xref ref-type="bibr" rid="scirp.42938-ref3">3</xref>]. Moreover, passive walking robots cannot walk on flat ground or an upward slope because the robots do not have actuators to recover mechanical energy. Addition of some control and actuation is therefore necessary to stabilize passive waking robots in complex environments. The walking robot with minimum control and actuation is called quasi-passive walking robot [2,4,5]. Quasi-passive walking robots are energy-efficient and retain the feature of passive walking gait. Collins et al. demonstrated that quasi-passive walking robots can walk on a flat ground with startling human-like gait with simple controls only, such as ankle push-off or hip actuation [<xref ref-type="bibr" rid="scirp.42938-ref2">2</xref>]. These robots can walk straight stably, but few quasi-passive walking robots can turn because they are symmetrically controlled and actuated.</p><p>A quasi-passive walking robot that can walk and turn on flat ground was demonstrated in our previous study [6-8]. It was experimentally demonstrated that synchronization of the period of lateral motion T<sub>L</sub> with the period of swing leg motion T<sub>S</sub> was a necessary condition for stable 3D passive walking [<xref ref-type="bibr" rid="scirp.42938-ref6">6</xref>]. In the next step, a mechanical oscillator actuated by a motor was mounted on a 3D passive walking robot to stabilize the robot [7,8]. The mechanical oscillator was always oscillated in the frontal plane to synchronize T<sub>L</sub> with T<sub>S</sub> by utilizing forced entrainment. The target trajectory of the mechanical oscillator is periodic, and thus is planned by adjusting the period, phase and amplitude of the target trajectory. The stabilization control method was numerically and experimentally demonstrated [7,8].</p><p>In variable environments, the ability to turn is a necessity for a biped robot in order to steer and avoid obstacles, as an example. In order to improve adaptability of the robot to changing environments, a turn control method was proposed and numerically examined [<xref ref-type="bibr" rid="scirp.42938-ref9">9</xref>]. However, the same method cannot appropriately apply to our experimental robot because the target trajectory of the mechanical oscillator becomes discontinuous when the stance leg changes. The power of the motor of the experimental robot is so limited that the mechanical oscillator cannot follow the discontinuous target trajectory and the experimental robot fails to turn.</p><p>In this study, we propose a novel turn control method by controlling the central axis of oscillation of the mechanical oscillator to enable the robot to turn stably on flat ground with turn radius controlled. This new method is examined experimentally and numerically. Since this method does not need to increase actuators or change the construction of the robot the turn control and stabilizing methods for straight walking can switch to each other directly.</p><p>Additionally, the gait of turn is compared with that of straight walking and analyzed in terms of mechanical work and energy.</p></sec><sec id="s2"><title>2. Experimental and Simulation Model</title><p>The experimental robot is made based on a passive walker, and the robot is composed of two straight legs, a trunk, a motor and a mechanical oscillator, as shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>(a). The legs are connected to hip axis by two passive joints, and the relative angles between the legs and the hip axis are measured by utilizing two rotary encoders, as shown in  <xref ref-type="fig" rid="fig1">Figure 1</xref>(b). The trunk is fixed to the hip axis, and the motor is mounted on the trunk to control the mechanical oscillator around the motor axis in the frontal plane. The batteries are fixed to the trunk to increase the moment of inertia about yaw axis and to lower the center of mass of the trunk. The center of mass of the trunk is lower than the hip axis because of the batteries, so the mechanical oscillator can keep upright. The robot’s foot sole is spherical and the center of the sphere is set higher than the center of mass of the robot to ensure that the robot is not only stable in standing posture but also robust against disturbance, initial conditions and path conditions in walking. The center of mass of the feet is adjusted backward by weights so that the swing leg can naturally swing forward even on a flat ground, as shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>(b).</p><p>The simulation model constructed on the Open Dynamics Engine (ODE) [<xref ref-type="bibr" rid="scirp.42938-ref10">10</xref>] has the same structure, mass distribution as the experimental model except for the feet, as shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>. The geometrical shape of the feet</p><p>is set to spherical the same as the experimental model, but the mass distribution of the feet is set cubic for simplicity. The world coordinate is OXYZ, and the roll, pitch and yaw are denoted by θ, γ and ψ, respectively. The double support phase is assumed to be instantaneous, and the motion that the swing foot reaches the ground is regarded as heel-strike, which is assumed to be inelastic and without sliding in simulation. In single support phase, the spherical stance foot purely rolls on the ground without slip, and the swing leg swings ahead like a pendulum. The friction in joints is set to zero in simulation.</p></sec><sec id="s3"><title>3. Algorithm of Turn Control</title><p>A simplified model of lateral motion in turn control is shown in <xref ref-type="fig" rid="fig3">Figure 3</xref>, where the trunk and legs are simplified to a block, with line AB as the central axis of the block, and line AC as the central axis of oscillation of the mechanical oscillator.</p><p>The roll angle of the lateral motion of the block is represented by θ, the inclination angle of the mechanical oscillator relative to line AB is represented by θ<sub>w</sub>, the inclination angle of the central axis of oscillation relative to line AB is represented by θ<sub>1</sub>, and the inclination angle of the mechanical oscillator relative to line AC is represented by θ<sub>2</sub>.</p><p>The target trajectory of θ<sub>w</sub>, θ<sub>1</sub> and θ<sub>2</sub> are represented by θ<sub>wt</sub>, θ<sub>1t</sub> and θ<sub>2t</sub>, respectively. The target trajectory θ<sub>wt</sub> is planned by θ<sub>1t</sub> and θ<sub>2t</sub>, because θ<sub>wt</sub> is equal to θ<sub>1t</sub> + θ<sub>2t</sub>. The turning radius is controlled by θ<sub>1t</sub>, and the robot is stabilized by θ<sub>2t</sub> based on the stabilization control method.</p><p>In the stabilization control method, the period of lateral motion T<sub>L</sub> is controlled and always synchronized with the period of swing leg motion T<sub>S</sub> by periodic oscillation of the mechanical oscillator. Following the method, in turn control θ<sub>2t</sub> is also periodic and thus is planned by control-</p><p>ling its period, amplitude and phase, respectively, as shown in <xref ref-type="fig" rid="fig4">Figure 4</xref>.</p><p>The period of θ<sub>2t</sub> is controlled on the basis of forced entrainment, which is an interesting phenomenon in nonlinear vibrations [<xref ref-type="bibr" rid="scirp.42938-ref11">11</xref>]. Forced van der Pol equation,</p><disp-formula id="scirp.42938-formula127966"><label>, (1)</label><graphic position="anchor" xlink:href="6-8102087\2df4fca3-a3f9-46c6-ae2a-9a5d395a55d1.jpg"  xlink:type="simple"/></disp-formula><p>is utilized to realize forced entrainment [<xref ref-type="bibr" rid="scirp.42938-ref7">7</xref>], where the roll angle θ of the robot is inputted into the Equation (1) as a periodic forcing function. The angular frequency of self-excitation of the Equation (1) is represented by Ω<sub>V</sub> and the angular frequency of θ is represented by ω. If Ω<sub>V</sub>@ω or the coefficient K is sufficiently large, the system indicates a phase-locking phenomenon and θ will entrain y. According to forced entrainment, the periods of y and <img src="6-8102087\09b4e070-6e65-4729-807d-425d213c5d62.jpg" /> are synchronized to the period of θ. The phase of <img src="6-8102087\d97ca25c-0467-40ef-aec7-fa894175d60e.jpg" />&#160;is the same as θ, but y shows a phase lead of π/2 relative to θ. The period of θ<sub>2t</sub> is controlled by y and <img src="6-8102087\041af939-6bae-4802-b4f2-c7af0e1edb6d.jpg" />, and thus the period of target trajectory is also synchronized with the period of lateral motion of the robot.</p><p>The amplitude β of the θ<sub>2t</sub> is determined by a proportional algorithm based on the stabilization control,</p><disp-formula id="scirp.42938-formula127967"><label>, (2)</label><graphic position="anchor" xlink:href="6-8102087\9d178e9c-8bae-467f-9a7b-38b248aedab4.jpg"  xlink:type="simple"/></disp-formula><p>where K<sub>P</sub> is the proportional gain, and α is a constant value obtained by preliminary simulation and determines the initial value of β.</p><p>According to y, <img src="6-8102087\42d5e9ea-33c1-46b8-8250-6747bef734ed.jpg" />, β, and phase difference φ, The target trajectory θ<sub>2t</sub> is determined as follows [<xref ref-type="bibr" rid="scirp.42938-ref7">7</xref>]:</p><disp-formula id="scirp.42938-formula127968"><label>, (3)</label><graphic position="anchor" xlink:href="6-8102087\946b2989-6c5a-463d-84c0-c9676920baa7.jpg"  xlink:type="simple"/></disp-formula><p>where C<sub>1</sub> and C<sub>2</sub> are the amplitudes of y and <img src="6-8102087\29628c6f-e633-4054-af5c-d706d2ee2b21.jpg" />. When β is positive, the phase difference φ between the target trajectory θ<sub>2t</sub> and the roll angle θ is set to 90˚ or −90˚, respectively, to increase or decrease T<sub>L</sub> most efficiently. When φ is set to 90˚, the phase difference is automatically selected as 90˚ or −90˚ according to the sign of β given by Equation (2).</p><p>The target trajectory of θ<sub>1</sub> is θ<sub>1t</sub>, which is planned to control the turning direction and turning radius r. When θ<sub>1</sub> is positive, the robot turns right, and when θ<sub>1</sub> is negative, the robot turns left. In order to investigate the relation between the turning radius r and θ<sub>1</sub>, r is measured when θ<sub>1</sub> is set to a constant value in ODE simulation, as shown in <xref ref-type="fig" rid="fig5">Figure 5</xref>. The vertical axis is r, and the horizontal axis is θ<sub>1</sub>, the data of which is symmetric about the vertical axis. The minimum turning radius is 0.53m when θ<sub>1</sub> is set to −60˚ or 60˚. If the absolute value of θ<sub>1</sub> is larger than 60˚, the robot cannot turn stably.</p><p>In order to control r, based on <xref ref-type="fig" rid="fig5">Figure 5</xref>, the relationship between θ<sub>1</sub> and r<sub>t</sub> is expressed by a function obtained by a curve fitting method, as</p><disp-formula id="scirp.42938-formula127969"><label>(4)</label><graphic position="anchor" xlink:href="6-8102087\10947d56-56f0-43b6-b3f8-2c53176fd1cf.jpg"  xlink:type="simple"/></disp-formula><p>The plus and minus signs “&#177;” in Equation (4) are used in left and right turn control, respectively.</p></sec><sec id="s4"><title>4. Simulation and Experiment</title><sec id="s4_1"><title>4.1. Turn Control</title><p>The turn control is realized in the ODE simulation. The robot walks straight on a flat ground for 3 seconds to stabilize the walking gait, then θ<sub>1</sub> is set to −60˚&#160;and the robot begins to turn right. The trajectory of the center of mass of the robot in right turn is shown in  <xref ref-type="fig" rid="fig6">Figure 6</xref>, where the horizontal and vertical axis are X and Y axis of the world coordinate, respectively. The turning radius is about 0.53 m, and the trajectory of the center of mass is snaky because the robot rolls in the lateral plane and the</p><p>center of mass of the robot moves between the left and right foot.</p><p>The period of right lateral motion T<sub>right</sub> is defined as the period of twice the time while θ is positive in one walking cycle, and the period of left lateral motion T<sub>left</sub> is defined as the period twice of the time while θ is minus in one walking cycle. In straight walking, T<sub>right</sub> and T<sub>left</sub> are the same because of the symmetric left and right lateral motion.</p><p>However, T<sub>right</sub> and T<sub>left</sub> become different in turn control by the inclination angle of the central axis of oscillation θ<sub>1</sub>. The relation of θ<sub>1</sub> to the period difference “T<sub>right</sub> - T<sub>left</sub>” is investigated in simulations and experiments, and the results are as shown in <xref ref-type="fig" rid="fig7">Figure 7</xref>, where the vertical axis is the period difference “T<sub>right</sub> - T<sub>left</sub>” and the horizontal axis is the inclination angle θ<sub>1</sub>. The period difference increases when θ<sub>1</sub> increases. The results of the simulation and experiments show that the method is effective to control the period difference “T<sub>right</sub> - T<sub>left</sub>”. According to <xref ref-type="fig" rid="fig5">Figure 5</xref> and <xref ref-type="fig" rid="fig7">Figure 7</xref>, if |T<sub>right</sub> - T<sub>left</sub>| is larger, the turn radius becomes smaller, because larger |T<sub>right</sub> -</p><p>T<sub>left</sub>| leads larger difference of strides between the left and right leg and makes the robot turn with smaller turn radius.</p><p>Turn is an important ability of the robot in variable environments, so the turn control method is examined with the experimental robot and simulation model walking on a path with different turn radiuses at bends on flat ground.</p><p>In the experiment the path has two bends and the turn radiuses are 0.75 m and 2.0 m, as shown in <xref ref-type="fig" rid="fig8">Figure 8</xref>(a), so the robot has to switch its turn radius to adapt to the changing environment. There are no external sensors on the robot to sense the changing environment, so θ<sub>1t</sub> is calculated in advance according to the turn radiuses. According to Equation (4), θ<sub>1t</sub> is calculated and adjusted to −58˚ and 35˚ in the experiments when the turn radiuses are set to 0.75 m and 2.0 m, respectively. The robot stably walks through the path in the experiments, as shown in <xref ref-type="fig" rid="fig8">Figure 8</xref>(b), which shows that the turn control makes it possible for the robot to walk stably even under different conditions.</p><p>In the simulation a more complicated path is built but control method is the same as the method in the experiment, in which θ<sub>1t</sub> is also calculated in advance according to the turn radiuses. The robot walks through the path stably, as shown in <xref ref-type="fig" rid="fig9">Figure 9</xref>, where the location of the robot is shown in the mini map of the path.</p></sec><sec id="s4_2"><title>4.2. Comparison of Turn and Straight Walking</title><p>In order to investigate the gait of the robot in turn control, the pitch angles of the right and left legs in right turn are compared with that in straight walking by simulation, as shown in  <xref ref-type="fig" rid="fig1">Figure 1</xref>0. The vertical axis is the pitch angle, and the horizontal axis is time. In straight walking the pitch angles of the two legs indicate similar wave shape, periods, and amplitudes yet opposite phase, and thus the time of left and right swing phase is the same, as shown in  <xref ref-type="fig" rid="fig1">Figure 1</xref>0(a). However, in right turn the pitch angles of the two legs indicate similar periods yet different wave shape and phases, and the time of the left swing phase is larger than the time of the right swing phase, as shown in  <xref ref-type="fig" rid="fig1">Figure 1</xref>0(b).</p><p>The lateral motion of the robot in right turn is compared with that in straight walking by using phase plane portraits, as shown in  <xref ref-type="fig" rid="fig1">Figure 1</xref>1. The vertical axis is the roll angle of the lateral motion θ, and the horizontal axis is the angular velocity of the lateral motion, with the black dot indicating the initial state, and the black arrow is the end of phase plane portraits. In straight walking the phase plane portraits are almost symmetric on the left and right half plane, and thus the lateral motion is symmetry in left and right stance phase, as shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>1(a). However, in right turn the phase plane portraits is asymmetric on the left and right half plane, as shown in</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>1(b), and the asymmetry comes from the inclination of the central axis of the mechanical oscillator movement. In right turn the roll angle of lateral motion in right stance phase is larger than that in left stance phase.</p><p>The unsymmetrical lateral motion of the robot can be understood from the viewpoint of mechanical work, so the positive and negative work performed by the motor on the robot in right turn are calculated and compared with that in straight walking in ODE simulation.</p><p>The positive and negative works in right turn are shown in  <xref ref-type="fig" rid="fig1">Figure 1</xref>2(a). The positive and negative works are represented by W<sub>p</sub> and W<sub>n</sub>, respectively. The left vertical axis shows the work performed by the motor on the robot, and the right vertical axis shows the roll angle of the robot θ. According to the roll angle θ, the right and left stance phase can be distinguished easily, because stance leg changes when θ is 0. In the right stance phase the motor performs more positive work than negative work, thus the mechanical energy, the amplitude of lateral motion of the robot and the time of right stance phase increases. In the left stance phase the motor performs more negative work than positive work, thus the mechanical energy, the amplitude of lateral motion of the robot and the period of right lateral motion decreases. Therefore, the period difference between the right and left lateral motion can be increased, as shown in <xref ref-type="fig" rid="fig7">Figure 7</xref>.</p><p>The positive and negative work in straight walking is shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>2(b). Although the mechanical oscil-</p><p>lator always accelerates and decelerates to periodically sway left and right in the frontal plane, positive work accounts for 91% of total mechanical work on average in both right and left stance phase. Therefore, the right and left lateral motion of the robot are symmetric, as shown in  <xref ref-type="fig" rid="fig1">Figure 1</xref>2(b).</p></sec></sec><sec id="s5"><title>5. Conclusion</title><p>In this research, a novel turn control method is proposed for a 3D quasi-passive walking robot, and the method is examined both numerically and experimentally. The turn radius of the robot can be controlled by the inclination angle of the central axis of oscillation of the mechanical oscillator. Based on the turn control, the robot successfully walks through a curved path in the experiment and simulation, which indicates that it is possible for the robot to walk in variable environments. In addition, the gait and mechanical work of turn control are investigated and compared with that in straight walking.</p><p>In our future work, external sensors will be used to detect changing environments in order to guide the robot to walk in complex environments. Moreover, the amplitude of the target trajectory of the mechanical oscillator will be calculated by utilizing the relation between work and energy in turn control.</p></sec><sec id="s6"><title>REFERENCES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.42938-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">T. McGeer, “Passive Dynamic Walking,” International Journal of Robotics Research, Vol. 9, No. 2, 1990, pp. 62-82.</mixed-citation></ref><ref id="scirp.42938-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">S. H. Collins, A. Ruina, R. Tedrake and M. Wisse, “Efficient Bipedal Robots Based on Passive Dynamic Walkers,” Science Magazine, Vol. 307, 2005, pp. 1082-1085.</mixed-citation></ref><ref id="scirp.42938-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">A. Goswami, B. Espiau and A. Keramane, “Limit Cycles and Their Stability in a Passive Bipedal Gait,” Proceedings of the IEEE International Conference on Robotics and Automation, Minneapolis, 1996, pp. 246-251.</mixed-citation></ref><ref id="scirp.42938-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">R. Tedrake, T. W. Zhang, M. F. Fong and H. S. Seung, “Actuating a Simple 3D Passive Dynamic Walker,” Proceedings of the IEEE on Robotics and Automation, Vol. 5, 2004, pp. 4656-4661.</mixed-citation></ref><ref id="scirp.42938-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">A. D. Kuo, “Stabilization of Lateral Motion in Passive Dynamic Walking,” International Journal of Robotics Research, Vol. 18, No. 9, 1999, pp. 917-930. http://dx.doi.org/10.1177/02783649922066655</mixed-citation></ref><ref id="scirp.42938-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">S. Suzuki and M. Hachiya, “Experimental Study on Stabilization of a Three-Dimensional Biped Passive Walking Robot,” Journal of the Society of Biomechanisms, Vol. 32, No. 4, 2008, pp. 239-246. http://dx.doi.org/10.3951/sobim.32.239</mixed-citation></ref><ref id="scirp.42938-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">M. Hachiya and S. Suzuki, “Stabilization of a Biped Quasi-Passive Walking Robot via Periodic Input,” Journal of the Society of Biomechanisms, Vol. 33, No. 1, 2009, pp. 57-63. http://dx.doi.org/10.3951/sobim.33.57</mixed-citation></ref><ref id="scirp.42938-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">S. Suzuki, M. Takada and Y. Iwakura, “Stability Control of a Three-Dimensional Passive Walker by Periodic Input Based on the Frequency Entrainment,” Journal of Robotics and Mechatronics, Vol. 23, No. 6, 2011, pp. 11001107.</mixed-citation></ref><ref id="scirp.42938-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">S. Suzuki, Y. Cao, M. Takada and K. Oi, “Climbing and Turning Control of a Biped Passive Walking Robot by Periodic Input Based on Frequency Entrainment,” Advanced Engineering Forum, Vol. 2, No. 3, 2011, pp. 4852.</mixed-citation></ref><ref id="scirp.42938-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">R. Smith, “Open Dynamics Engine v0.5 User Guide,” 2006. http://ode.org</mixed-citation></ref><ref id="scirp.42938-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">R. Rand, “Lecture Notes on Nonlinear Vibrations,” 2012. http://ecommons.library.cornell.edu/handle/1813/28989</mixed-citation></ref></ref-list></back></article>