Abstract
Despite recent major advances in computational power and control algorithms, the stable and robust control of a bipedal robot is still a challenging issue due to the complexity and high nonlinearity of robot dynamics. To address the issue an efficient and powerful alternative based on a biologically inspired control framework employing neural oscillators is proposed and tested. In a numerical test the virtual force controller combined with the neural oscillator of a humanoid robot generated rhythmic control signals and stable bipedal locomotion when coupled with proper impedance components. The entrainment nature inherent to neural oscillators also achieved stable and robust walking even in the presence of unexpected disturbances, in that the centre of mass (COM) was successfully kept in phase with the zero moment point (ZMP) input trajectory. The efficiency of the proposed control scheme is discussed alongside simulation results.
1. Introduction
Keeping the stability and robustness of bipedal locomotion has recently gained a lot of attention in the humanoid robot community since considerable efforts are required to deal with both the high nonlinearity of multibody robot dynamics and the unexpected disturbances from new environments. Modern robot control based on motion planning based approaches, i.e., the zero moment point criterion [1, 2], the three-dimensional linear inverted pendulum model [3] and virtual model-based control [4], have tried to achieve stability and robustness. For motion planning a fairly accurate robot model with a large number of degrees of freedom is required to generate dynamically admissible motion patterns. However, an efficient pattern generation of such a large dynamic system is still quite challenging withstanding recent advances in computational power. Thus, from a practical viewpoint, a biologically inspired control approach that does not need a reference motion may be a good alternative to enhancing pattern generation competence.
Human and animal locomotion inherently exhibit stable rhythmic movements, adapted to the natural frequency of their body dynamics, in spite of variety in their sensors and actuators. The neural circuits of oscillators on the spinal cord [5] known as Central Pattern Generators (CPGs) contribute to such efficient motor movement and novel stability properties in biological models. Moreover, neural oscillators can entrain to sensory feedback, which plays a key role in adapting locomotion to unexpected forces from new environment. Hence, biologically inspired motion control based on the neural oscillator has drawn much attention as an attempt to realize natural motion in bipedal humanoid robots.
The neural oscillator for robotic applications was first formalized in mathematical form by Matsuoka [6,7]. In the paper the neurons were proven to generate rhythmic output and the necessary and sufficient conditions for self-oscillations were provided. Through Matsuoka's model, Taga et al. showed that the neural oscillators of a bipedal robot could entrain to the characteristics of various environments, so as to generate stable bipedal locomotion by the feedback of sensory signals from robot joints to the neural oscillators [8, 9]. In a 2D simulation robot control became robust against perturbation and the robot was able to climb an upward slope [10]. These methodologies were later applied to a 3D simulation (Miyakoshi et al. [11]), to dynamic quadrupedal walking (Fukuoka et al. [12]) and to rhythmic robotic arm movement (Williamson [13]).
The entrainment property of the neural oscillator was later tested on real systems [15-17] such as the robot walking on uneven terrains. Even though it enabled adaptive bipedal locomotion in a 3D environment, the neural oscillator achieved global stability mostly through its phase locking properties, which sometimes cause large deviations in motion from the reference one [18]. In addition, since previous research mostly focused on cyclic gait motion, such as bipedal walking, direct application of the neural oscillator based control schemes are not easy in robot practices, which need to control various robot behaviours that are non-cyclic. Therefore, it would be desirable for a user to be able to handle overall robot motion, i.e., ZMP trajectory, to fulfil a mission objective, while allowing autonomous adaptation of robots for stability purposes. On the other hand, a new intuitive stability control scheme based on virtual model control for legged locomotion was proposed by Pratt [14]. Even though this algorithm is extended for rough terrain walking, the approach is still limited to 2D.
For practical purposes new methods in this study are designed to maintain the stability autonomously without sophisticated stability controllers with a predetermined ZMP trajectory input so that humanoid robots could show arbitrary whole body motions as well as various bipedal motions with minimal interruption by a user. The new methods that couple the neural oscillator and the virtual components apply the interaction between neural oscillators and virtual mechanical components to humanoid locomotion control and produce periodic swing motions, which allow the adaptivity of humanoid robots to unknown disturbances through entrainment of their frequency and phase.
In addition, the methodologies on how to generate human-like motions while properly maintaining stability of humanoid robots will be discussed. In the proposed self-stabilizing controller, the motion of the inverted pendulum is often regarded as a supporting leg in the stance phase of humanoid locomotion [19]. We therefore develop a new control method for sustaining and enlarging the stability of an inverted pendulum in the pitching and rolling direction. The proposed controller contributes to whole body motion generation and adaptive bipedal locomotion of humanoid robots even under unknown disturbances and uneven terrains. Based on this self-stabilizing controller, we focus on the implementation of the new bipedal walking strategy. This is verified through simulations and experiments designed to test whether the proposed approach yields a robust and efficient control of rhythmic locomotion
This paper is organized as follows: Section 2 briefly introduces the mathematical form and inherent properties of Matsuoka's neural oscillator. In Sections 3 to 5, we address basic concepts for controlling bipedal robots and how to implement these in humanoid robots with basic tests on dynamic stability through simulations. More details of the bipedal stability on the proposed control model under unknown disturbances are discussed in Section 5. Section 6 presents the proposed bipedal walking strategy with simulations and experiments. Conclusions are drawn in Section 7.
2. Rhythmic Movement Generation
A. Matsuoka's Neural Oscillator
Our work is motivated by studies on and aspects of biologically inspired locomotion control employing oscillators. Especially, the fact that the basic motor pattern generated by the Central Pattern Generator (CPG) of the inner body of a human or animal is usually modified by sensory signals from motor information to deal with environmental disturbances. Similarly to the sensory system of a human or animal, the neural oscillators are entrained with external stimuli at a sustained frequency. They show stability against perturbations through global entrainment between the neuro-musculo-skeletal systems and the ground [8]. Thus, neural oscillators have been applied to the CPG of humanoid robots with rhythmic motions [9, 10]. The oscillators provide robust performance in a wide variety of rhythmic tasks when they are implemented to a system such as a robotic arm. The reason is that the oscillators use sensory signals about the joint state to adapt the frequency and phase of the joint motion regardless of the references corresponding to changes in the environment.
Matsuoka's neural oscillator consists of two simulated neurons arranged in mutual inhibition, as shown in Fig. 1 [5, 6]. If gains are properly tuned, the system exhibits limit cycle behaviours. The trajectory of a stable limit cycle can be derived analytically, describing the firing rate of a neuron with self-inhibition. The neural oscillator is represented by a set of nonlinear coupled differential equations given as

Schematic diagram of the NEURAL OSCILLATOR
where
B. Coupling Neural Oscillator to Mechanical Systems
This subsection addresses a new control method that exploits the natural dynamics of the oscillator coupled with the dynamic system that closely interacts with the environment. This method enables the robot to adapt to changing conditions. For simplicity, we employed a general 2nd order mechanical system connected to the neural oscillator, as seen in the lower system in Fig. 3. The desired torque input at the

Mechanical system model coupled to the neural oscillator
where
Generally, it is known that Matsuoka's neural oscillator exhibits the following properties: the natural frequency of the output signal increases in proportion to 1/
If we appropriately tune the parameters of the neural oscillator, the oscillator exhibits stable limit cycle behaviours. In Fig. 1, the gain

Simulation result on entrainment property of the neural oscillator.
3. The Control Model for Bipedal Stability
In humanoid locomotion, the pitching motion should be performed under the stable single support phase of the rolling motion. Now we explain how to attain the stable single support phase corresponding to the periodic bipedal locomotion. For the stable rolling motion in the frontal plane, we consider an inverted pendulum model coupled with such a virtual mechanical component as a spring and damper and the neural oscillator, as seen in Fig. 5 (a) for generating an appropriate rolling motion. The coupled model enables the inverted pendulum to stably move in a frontal plane according to a desired ZMP (Zero Moment Point) trajectory sustaining the stability.
Assuming that θ, the angle between the vertical axis and the pendulum in Fig. 5 (a), is small enough and linearized near 0, the dynamic equation of the coupled inverted pendulum is given by
where
If the desired ZMP trajectory,
where

(a) The proposed control model for stable inverted pendulum motion (b) Simulation result with respect to rolling motion of humanoid robot (c) The rolling COM motion of humanoid robot
As illustrated in Fig. 5 (b) and (c), the humanoid robot exhibits a stable rolling motion, therefore satisfying the desired ZMP. Assuming variation in environmental conditions, the robustness of the proposed control method is verified with the simulation result in Fig. (6). In the simulation results in (a) and (b) in Fig. 6, when the coupled model is implemented to the pendulum model (see the Fig. 6 (b)), the COM motion is sustained stably even with an allowable change of a ZMP input.
Parameters of the neural oscillator and control model

Simulation result with respect to the COM motion of humanoid robot in terms of an arbitrary ZMP input without (a) and with the coupled model (b)
4. Implementation and Basic Tests of The Coupled Model
A. Implementation
This subsection describes how to embody the bipedal locomotion of the humanoid robot using the proposed control scheme to guarantee stability on the single and double support phases. The humanoid robot in Fig. 7, developed by KIST, is used to test and evaluate the proposed controller. It has 35 degrees of freedom (DOF): 6 in each arm, 6 in each leg, 4 in each hand, 2 in its neck and 1 in its waist. It's 150cm high and weighs about 67kg. The robot is equipped with a stereo camera, a microphone, four force/torque sensors and a pose sensor.

KIST Humanoid robot
We propose a control method that uses several virtual mechanical components such as springs and dampers and couples them with neural oscillators, as seen in Fig. 8.

The coupled model with virtual springs and dampers
A virtual spring-damper is defined from the hip to the ankle for each leg and could generate the virtual vertical forces
Considering the right and left legs as separate systems, the dynamic equations of each leg can be written in the following form
where M is the 3 × 3 inertia matrix,
where
where
B. Basic Control Tests of the Coupled Model
To test the performance of the proposed control methods, a dynamic humanoid motion was simulated using a dynamic robot engine, SIMSTUDIO (SimLab, KOREA). Prior to a complete cycle of bipedal locomotion, tests on basic bipedal motions were simulated. From the simulation results the adaptive properties of the neural oscillator were investigated.
The rolling and pitching of humanoid are two basic bipedal motions when ZMP trajectory is given as a function over time. The key dynamic parameters and snapshots of a whole robot body are presented in Figs. 9–12.

Input ZMP profile (dashed line), the output of the COM position (solid line) and the output of the neural oscillator in the rolling motion (dotted line)

Snap shots on the stable rolling motion of a humanoid robot

Input ZMP profile (dashed line), the output of the COM position (solid line) and the output of the neural oscillator in the pitching motion (dash-dotted line)

Snap shots on the stable pitching motion of a humanoid robot
The dashed lines in Figs. 9 and 11 indicate the input ZMP trajectory over time. It is shown that the COM trajectory presented in dashed lines properly oscillates in response to the arbitrarily defined ZMP input, so the humanoid can maintain stability. This implies that the proposed control method in Section 3 works properly. For example, let's investigate the outputs of the neural oscillators presented in dotted lines in Figs. 9 and 11. The COM motion of the humanoid robot is considered as the sensory signal of the neural oscillator. Then the outputs of neural oscillator entrain to that of the humanoid robot and drive the humanoid robot corresponding to the sensory input appropriately as a torque input. Consequently, we note that this leads the adaptive motion of humanoid robots to maintain bipedal stability even under an unknown disturbance. In addition, the ZMP planning in the proposed control method is possible in terms of properly adjusting the gains of the impedance model and the neural oscillator described in Eq. (4). In Figs. 9 and 11, the simulations were performed with the consideration that the real ZMP is larger than the movement of the COM within the range of the ZMP criterion.
5. Adaptive Motion under unknown Disturbance
In this section, the inherent adaptive feature of neural oscillators is tested through the proposed balance controller. The proposed controller is based on a neural oscillator coupled with the virtual components under two states of the humanoid robot under unknown disturbances. Fig. 13 illustrates how to apply an external acceleration as a disturbance to the MAHRU robot. The size and direction of the red s indicates the magnitude and direction of the external acceleration as unknown disturbances,. The unknown external acceleration is induced in the COM position of the MAHRU robot periodically as time passes.

Schematic figure on an unknown harmonic disturbance applied to the humanoid robot
From the results shown in Figs. 14 and 15, it can be observed that the humanoid robot not coupled to the neural oscillator and virtual components tips over if there is an unknown disturbance. In the graphs in Figs. 15 and 17, the red thin, the blue dashed and the grey thick lines indicate the calculated ZMP, the input of an external acceleration as an unknown disturbance and the COM position, respectively. Also, note that the black dotted line is the output of the neural oscillator and the red dash-dotted line is the desired ZMP input. The COM and the ZMP plots of Fig. 15 represent that the humanoid robot became unstable. Here the positive and negative signs of the moving distance imply the left and right direction of the humanoid robot. In addition, it can be verified from Figs. 16 and 17 that the COM follows the external disturbance well regardless of the motion of the humanoid robot. This is caused by the neural oscillator coupled with the COM of the humanoid robot, since the neural oscillator entrains the sensory signal fed from the COM motion of the humanoid robot.

Snap shots on the motion of the MAHRU robot under an unknown sinusoidal disturbance of 1.9G without the neural oscillator; G is set with 9.8 m/s2

The outputs of the COM position and the ZMP of the humanoid robot; The input of the external acceleration

Snap shots on the motion of the MAHRU robot under an unknown sinusoidal disturbance of 1.9G with the neural oscillator based controller. G is set with 9.8 m/s2

The outputs of the COM position and the ZMP of the humanoid robot; The input of the external acceleration; The output of the neural oscillator
The proposed control method was implemented on the real humanoid robot to verify the validity of the proposed motion adaptor with the self-stabilizing bipedal controller. The kinematic and dynamic properties of a robot such as the width of the feet, the area of each foot and the length and mass of each leg are different from those of a human. So, considering the calculated COM, we adjust the ZMP measured from the human motion so that the bipedal stability of the humanoid robot is maintained within the stable ZMP range. This is done by multiplication by a suitable scale factor. The modified ZMP is input into the self-stabilizing bipedal controller. Then the stable COM corresponding to the ZMP input is autonomously generated in terms of the self-stabilizing balance controller, as illustrated in Fig. 5 (c). Under the dancing motion of the humanoid robot, to employ the motion data obtained from the motion conversion method, unknown disturbances composed of two dumbbells weighing 10kg are applied directly to the humanoid robot on the frontal and sagittal sides, as seen in Fig. 18. It is observed from the experiment that the humanoid robot is able to maintain bipedal stability. Hence, even though the external disturbances of a sinusoidal form effect to the humanoid robot, the humanoid robot coupled to the neural oscillator and virtual components can stably exhibit a novel adaptive motion corresponding to an unknown disturbance.

Snap shots on the whole body motion and balancing test of humanoid robot under unknown disturbances, Dotted circles indicate the dumbbells that produce unknown disturbances
6. Bipedal Locomotion Strategy
The authors finally propose the appropriate bipedal walking strategy of humanoid robots for our control model. Our objective is to embody the simple and integrated algorithm for generating various bipedal motions of humanoid robots. We also aim to have robust and adaptive properties under an unknown disturbance even in an unpredictable terrain. The proposed algorithm is illustrated in Fig. 19. Comparing the bipedal walking strategy with the whole body motion generation, the only difference between them is that the input is not a reference ZMP but a rhythmic pattern generated by the neural oscillator. However, a desired ZMP pattern can be input into the self-stabilizing controller for bipedal walking as illustrated in Section 3. In experiments with respect to the bipedal walking of the humanoid robot, a simple sinusoidal ZMP pattern was used.

The proposed bipedal locomotion strategy
The output of the COM motion generated in the rolling direction dominates the entire bipedal walking motion, such as walking period, velocity, stride distance etc. The output of the rolling COM motion causes the COM motion in the pitching direction because the neural oscillator that produces the pitching motion is internally connected to the rolling oscillator. In this way, the rolling motion and pitching motion of the COM are harmonized periodically. Also the output of the COM motion in the lateral plane determines the trajectory of the swing motion. The sequence in the locomotion phases is shown in Fig. 23. If the COM is inhibited by an unknown disturbance, such as pushing or pulling forces, the neural oscillator generates the adapting motion properly by means of sensory feedback. In addition, if the swinging foot doesn't contact the ground thoroughly, the swing motion and the single support phase of individual legs are sustained. On the contrary, during the initial phase, if the swinging leg absorbs much of the impact energy in its landing on the ground, the module, in order to distinguish between distinct walking phases, selects a certain phase according to the ZMP and COM position. The initial phase begins when the swing leg hits the ground and lands on it and ends when the walking phase is changed into the double support phase. By this walking strategy, the swing and stance motions can be autonomously changed according to the altered COM motion. Ultimately the COM of a humanoid robot will be appropriately controlled, sustaining the bipedal stability according to an arbitrary ZMP input. Figs. 20 and 21 show the stable walking of MAHRU robot on flat terrain, controlled in terms of the proposed walking algorithm. In order to verify the validity of the proposed walking algorithm under unknown environmental changes, the height of the ground was altered aperiodically, irrespective of walking phases during bipedal walking. As shown in Fig. 22, when a swinging leg lands on the ground, the environmental change with regard to the height of the ground is applied. At first, while the left leg supports the whole weight of the robot, the height of the left leg changes in landing on the ground in about 7s. Therefore, the condition is applied to each walking phase from 7s to 11s in sequence. In Fig. 22, the red arrows indicate the direction that the 1.5cm high plate is inserted between the foot of the robot and the ground. Even though an unknown environmental condition exists by changing the height of the ground during bipedal walking it was verified from the experimental results of Fig. 21 that the proposed self-stabilizing controller and walking algorithm enabled a humanoid robot to exhibit efficient adaptable motions sustaining bipedal stability. However, most of the existing problems associated with unknown terrain as well as additional uncertainty in dynamic characteristics remain largely undetermined so far. This work will be extended to solve such a difficult problem in the future.

Simulation Result on MAHRU flat terrain walking

ZMP profile on flat terrain walking; the red thin line and the black dotted line indicate the reference ZMP input and the ZMP output, respectively.

Snap shots on bipedal walking of the humanoid robot under unknown changes of the ground 1.5cm in height.
7. Conclusion
To attain a stable periodic locomotion in biped humanoid robots, we proposed a new control architecture consisting of neural oscillators and virtual components. This causes the COM position of the biped robot to be controlled to follow the time varying desired ZMP input, thereby sustaining the bipedal stability. As a consequence, a humanoid robot can maintain bipedal stability regardless of single or double phases. For this, we simply let the virtual components and the neural oscillator generate the Cartesian forces. We can determine the required joint torques incorporating the Jacobian of the three link inverted pendulum. It is also noted that the appropriate swing and stance motions were generated according to the rolling and pitching motion because the neural oscillators between the dynamic models of the rolling and pitching direction are connected internally in phase. Since the stable rolling motion of the COM was induced by the ZMP reference input, which properly generated the pitching motion of the COM in accordance with the rolling motion, stable walking could be achieved with this novel walking strategy. Extensive simulations and experiments were carried out to verify the effectiveness of the proposed method. More experiments are currently under way for further evaluation of the proposed control method.
