Biped Walking of Passive Dynamic Walker with IPMC Linear Actuator

Similar documents
The Effect of Semicircular Feet on Energy Dissipation by Heel-strike in Dynamic Biped Locomotion

Dynamics of Heel Strike in Bipedal Systems with Circular Feet

Enhanced IPMC actuation by thermal cycling

EAP - ElectroActive Polymers

Limit Cycle Walking on Ice

Tailoring actuation of ionic polymer-metal composites through cation combination

Operator-based Modeling for Nonlinear Ionic Polymer Metal Composite with Uncertainties

Controlled actuation of Nafion-based Ionic Polymer-metal Composites (IPMCs) with Ethylene Glycol as Solvent

An Earth Auger as Excavator for Planetary Underground Explorer Robot. Using Peristaltic Crawling

Reinforcement Learning of Potential Fields to achieve Limit-Cycle Walking

Analysis of Circular Robot Jumping by Body Deformation

Inverse Dynamics of Gel Robots made of Electro-Active Polymer Gel

High-Performance PEDOT:PSS/Single-Walled Carbon Nanotube/Ionic liquid Actuators Combining Electrostatic Double-Layer and Faradaic Capacitors

Characteristics Evaluation of PVC Gel Actuators

Novel Reaction Force Control Design Based on Biarticular Driving System Using Intrinsic Viscoelasticity of Muscle

Characterization of the Electromechanical Properties of Ionomeric Polymer-Metal Composite (IPMC)

151H0653H00L Bio-Inspired Motor Control

CHANGING SPEED IN A SIMPLE WALKER WITH DYNAMIC ONE-STEP TRANSITIONS

Inverse Kinematics of Gel Robots made of Electro-Active Polymer Gel

Evaluation of a surface acoustic wave motor with a multi-contact-point slider

Fast Low-Voltage Electroactive Actuators. Using Nanostructured Polymer Electrolytes

Design and Evaluation of a Gravity Compensation Mechanism for a Humanoid Robot

Study of flapping actuator modules using IPMC

Actuation of Artificial Muscle Based on IPMC by Electromyography (EMG) Signal

Rhythmic Robot Arm Control Using Oscillators

Identification of a Thin Flexible Object with Bipedal Gaits

Level-Ground Walk Based on Passive Dynamic Walking for a Biped Robot with Torso

Surface Acoustic Wave Linear Motor

Feasibility Study on IPMC Actuated Contractile Water Jet Thruster Body Contraction

A Novel Propulsion Method of Flexible Underwater Robots

Electro-statically stricted polymers (ESSP)

Dynamics in the dynamic walk of a quadruped robot. Hiroshi Kimura. University of Tohoku. Aramaki Aoba, Aoba-ku, Sendai 980, Japan

DEVELOPMENT OF JUMP ASSIST SYSTEM USING PNEUMATIC RUBBER MUSCLE

Integration simulation method concerning speed control of ultrasonic motor

Ion Concentration and Electromechanical Actuation Simulations of Ionic Polymer-Metal Composites

DEVELOPMENT OF A BAR-SHAPED ULTRASONIC MOTOR FOR MULTI-DEGREES OF FREEDOM MOTION

ARTISAN ( ) ARTISAN ( ) Human-Friendly Robot Design

Reading. Realistic Character Animation. Modeling Realistic Motion. Two Approaches

CONTROLLED SYMMETRIES AND PASSIVE WALKING

Bipedal Walking Gait with Variable Stiffness Knees

Angular Momentum Based Controller for Balancing an Inverted Double Pendulum

SHAPE MEMORY ALLOY ACTUATOR PROTECTED BY ROLLED FILM TUBE FOR ARTIFICIAL MUSCLE

Surface Acoustic Wave Atomizer with Pumping Effect

A consideration on position of center of ground reaction force in upright posture

Coordinating Feet in Bipedal Balance

Natural Frequency Analysis of Spring-Manipulator System for Force Generation Utilizing Mechanical Resonance

Robust Low Torque Biped Walking Using Differential Dynamic Programming With a Minimax Criterion

Derivation and Application of a Conserved Orbital Energy for the Inverted Pendulum Bipedal Walking Model

Balancing of an Inverted Pendulum with a SCARA Robot

Vibration Suppression Control of a Space Robot with Flexible Appendage based on Simple Dynamic Model*

Lecture Presentation Chapter 8 Equilibrium and Elasticity

Acknowledgements. Feedback Control of Bipedal Locomotion RAS Pioneer Award to Mark Spong. Videos and Papers

Prototyping Pneumatic Group Actuators Composed of Multiple Single-motion Elastic Tubes

Micro Pneumatic Curling Actuator - Nematode Actuator -

An Energy Circulation Driving Surface Acoustic Wave Motor

All questions are of equal value. No marks are subtracted for wrong answers.

Novel Actuating System Based on a Composite of Single-Walled Carbon Nanotubes and an Ionomeric Polymer

University of Twente

Chapter 12. Recall that when a spring is stretched a distance x, it will pull back with a force given by: F = -kx

Mechanical energy transfer by internal force during the swing phase of running

MCE503: Modeling and Simulation of Mechatronic Systems

Bringing the Compass-Gait Bipedal Walker to Three Dimensions

A passively safe cable driven upper limb rehabilitation exoskeleton

Semester I lab quiz Study Guide (Mechanics) Physics 135/163

Development of fish robot using stacked-type electrostatic actuators

Real-Time Implementation of a LQR-Based Controller for the Stabilization of a Double Inverted Pendulum

Efficient Swing-up of the Acrobot Using Continuous Torque and Impulsive Braking

Jumping via Robot Body Deformation - Mechanics and Mechanism for Higher Jumping

Design and Control of Variable Stiffness Actuation Systems

High speed and quick response precise linear stage system using V-shape transducer ultrasonic motors

Detumbling an Uncontrolled Satellite with Contactless Force by Using an Eddy Current Brake

Power Assist H Control of Shift Lever with Spring Connected Link

EXTENDED GRIPPING CONDITIONS OF ROCK CLIMBER-LIKE ROBOT FOR ASYMMETRIC GRIPPING CONFIGURATION IN MICROGRAVITY

Proposal of Step Climbing of Wheeled Robot Using Slip Ratio Control

Mathematical Model for Pressure-Deformation Relationship of Miniaturized McKibben Actuators. {ashwinkp,

Dynamic Reconfiguration Manipulability Analyses of Redundant Robot and Humanoid Walking

The Dynamic Postural Adjustment with the Quadratic Programming Method

Safe Joint Mechanism using Inclined Link with Springs for Collision Safety and Positioning Accuracy of a Robot Arm

AP Physics C! name I CM R 2. v cm K = m

Asymptotically Stable Gait Primitives for Planning Dynamic Bipedal Locomotion in Three Dimensions

Swinging-Up and Stabilization Control Based on Natural Frequency for Pendulum Systems

Friction Drive Simulation of a SAW Motor with Slider Surface Texture Variation

A Backstepping control strategy for constrained tendon driven robotic finger

A) 4.0 m/s B) 5.0 m/s C) 0 m/s D) 3.0 m/s E) 2.0 m/s. Ans: Q2.

Hiroshi Kaminaga 1, Satoshi Otsuki 1, and Yoshihiko Nakamura 1

for Articulated Robot Arms and Its Applications

Joint Torque Control for Backlash Compensation in Two-Inertia System

Positioning Servo Design Example

Decoupling Identification for Serial Two-link Robot Arm with Elastic Joints

Chapter 2 Surface Acoustic Wave Motor Modeling and Motion Control

General procedure for formulation of robot dynamics STEP 1 STEP 3. Module 9 : Robot Dynamics & controls

ieee transactions on ultrasonics, ferroelectrics, and frequency control, vol. 45, no. 5, september

COMPLIANT CONTROL FOR PHYSICAL HUMAN-ROBOT INTERACTION

Fig.1 Partially compliant eccentric slider crank linkage

PERFORMANCE OF HYDROTHERMAL PZT FILM ON HIGH INTENSITY OPERATION

The Effect of Viscosity and Ion Size Copley, Hubbard, Maisano on the Transduction of Ionic. Polymer Metal Composite

Active elastomer components based on dielectric elastomers

Fundamental properties of snake robot locomotion

Stable Limit Cycle Generation for Underactuated Mechanical Systems, Application: Inertia Wheel Inverted Pendulum

Experimental Investigation on the Power Generation Performance of Dielectric Elastomer Wave Power Generator Mounted on a Square Type Floating Body

Transcription:

Biped Walking of Passive Dynamic Walker with IPMC Linear Actuator N. Kamamichi 1,Y.Kaneda 1, M. Yamakita 1,3,K.Asaka 2,3 and Z. W. Luo 3 1 Tokyo Institute of Technology, 2-12-1 Meguro-ku Oh-okayama Tokyo, JAPAN 2 Special Division for Human Life Technology, AIST, 1-8-31 Midirigaoka Ikeda Osaka, JAPAN 3 Bio-Mimetic Control Research Center, RIKEN, Anagahora Shimoshidami Moriyama-ku Nagoya, JAPAN {nkama,ykaneda,yamakita}@ctrl.titech.ac.jp Abstract: We are developing a linear artificial muscle actuator using ionic polymer-metal composites (IPMC) which is an electro-active polymer that bends in response to electric stimuli and the goal of our study is to apply the actuator to robotic applications especially to a biped walking robot. In this paper, we will describe an empirical model of the actuator and demonstrate walking simulations of a small-sized biped walking robot called Wire Walker. Keywords: Passive Dynamic Walking, Artificial Muscle Actuator, IPMC, EAP 1. Introduction The study of humanoid robot and bipedal walking is one of the most active research area in robotics. Humanoid walking has been already demonstrated, however, energy-effective gait design and safe control against outside environment and human being are not confirmed completely yet. One of the solutions for energy-effective gait design is exploiting the natural dynamics of a walking robot. McGeer studied a simple unpowered walking machine which is well known to as passive dynamic walking 1). A bipedal robot with suitable parameters can walk down a gentle slope without any driving force utilizing only gravity effect and its physical dynamics, and generates a stable periodic gait automatically. The gait pattern of passive walker can be considered to be natural and energy-efficient. In our previous works 2, 3, 4), we have studied some application methods of gait synthesis for active walking on the level ground based on the passive dynamic walking. A control strategy based on passive dynamic walling has properties of automatic gait generation and realization of dynamic-based control. The validity of proposed methods had been investigated by numerical simulations and experiments. In the experimental verification, we considered that the geared motor was not able to exploit the natural dynamics fully, and soft actuation should be needed in order to exploit the natural dynamics. For biped walkers, soft and light actuators with back-driveability should be needed to realized energy-effective gaits by taking advantage of the natural dynamics and the gravity effect. Some researchers have applied soft actuators to biped walking robots, e.g., Direct Drive torque motor, MIT Leglab s Series Elastic Actuator 5), pneumatic McKibben muscle actuator. The study of M. Wisse et al. is especially excellent, they constructed a 2D autonomous walking machine: MIKE 6) which is actuated by McKibben muscle actuators, and demonstrated its stable and energy-effective walking on the level ground. In this paper, we apply a high polymer gel actuator, Ionic Polymer-Metal Composite(IPMC), to a biped walking robot. High polymer gel actuator is one of the candidates of artificial muscle actuators due to their softness and miniaturizability. For several decades, electro-active polymers(eap) 7), which respond to electric stimuli with shape change, received little attention because of their actuating limitation. During last ten years, development of EAP materials with large displacement and quick response changed the potential capability, and EAP received much attentions of engineer and researcher from many disciplines, e.g., robotics, medical service, toy industry. Ionic Polymer-Metal Composite(IPMC) actuator is one of the most promising EAP actuators for applications. The output force of IPMC actuator is still small, however, we proposed a new concept of unitization in order to realize desired characteristic 8). The structure of our proposed actuator can be easily constructed, and the actuator transform bending motion into linear motion. The elementary units can be connected in parallel and series in order to realize desired displacement and force. In this paper, we explain the structure of the actuator, and identify an empirical model of the actuator. Numerical simulation of a small-sized bipedal walking based on the model and basic experiment are demonstrated. 2. Ionic Polymer-Metal Composite (IPMC) Ionic Polymer-Metal Composite (IPMC) is an EAP material that bends in response to an electrical activation, and produced by chemically plating gold or platinum on a perfluorosulfonic acid membrane which is known as an ion-exchange membrane. When input voltage is applied to metal layers of the both surface, it bends in high speed (Fig. 1). The phenomenon of this motion was discovered by Oguro et al. in 1992 9). The charac-

elementary unit Figure 3: Basic concept of IPMC linear actuator Figure 1: Bending behavior of IPMC actuator Figure 2: Structure of IPMC actuator teristics of IPMC are as follows: Driving voltage is low. (1 2 [V]) Speed of response is fast. (> 100 [Hz]) It is durable and stable chemically. (It is possible to bend over 1 10 6 times) It is a flexible material. It moves in water and in wet condition. Miniaturization and weight saving is possible. It is silent. 0[V] 2.5[V] By exploiting the characteristics, IPMC actuators have been applied to robotic applications such as active catheter 10, 11, 12), a wiper for nanorover 13), a micropump 14), a fish-type underwater robot 12, 15), micromanipulator 16) and a distributed actuation device 17). It, however, also has disadvantages that actuation force is still small and that input voltage is restricted to the range that electrolysis of ionic polymer dose not occur. In order to improve the performance, developments of ionic polymer menbrance and plating method are required. 3. IPMC Linear Actuator In this section, a concept of linear actuator and the structure of its elementary unit is explained, and a performance of the actuator is analyzed from experimental results. load Laser Displacement Meter IPMC Actuator [current] Computer Figure 4: Experimental setup [position] [voltage] 3.1 Structure of Linear Actuator We propose an effective structure of actuator that transform bending motion into contracting motion. Fig. 2 shows the structure of the elementary unit. This elementary unit consists of four IPMC films and one side of the unit is formed from a pair of films which are connected by a flexible material. When input voltage is applied to electrodes of surface with anode outside, each membranes bend outside, then the actuator is constricted. An actuation force and displacement of each units are small, however, the elementary units are able to be connected in parallel and series as in Fig. 3, so the actuator can realize desired force and displacement. By shifting series of the elementary units by half pitch to avoid interferences as in Fig. 3, total actuator is made compact, and high power/volume and miniaturization are realized. 3.2 Characteristics of Actuator In order to check the characteristics of the actuator, we carried out fundamental experiments. Fig. 4 shows the experimental setup. In this experiment, one edge of the actuator is fixed on a board floated on the water surface to reduce effects of the weight of electrodes and ties. Displacement of the linear actuator was measured by a laser displacement meter. IPMC film which we used in this experiment is Nafion R 117(DuPont) plated with gold, and the amount of gold on the surface of membrane is about 10 [mg/cm 2 ]. The elementary unit was made from four IPMC films which were cut into ribbons with a width of 2 [mm] and length of 20 [mm], and the total length of actuator unit is about 40 [mm]. The maximum displacement was more than 10 [mm] under a step voltage

$ # "! (a) displacement (b) current $ # "! (a) displacement % % % (b) current % % % Figure 5: Step response with various inputs Figure 6: Step response with various loads of 2.5 [V], though the response of the actuator varies depending on its condition. Response in Step Voltage Fig. 5 shows the response in a step voltage without loads, where step input voltages of 1.5, 2.0, and 2.5 [V] were applied at 0 [s]. From (a), it can be seen that the peak value of displacement enlarges as input voltage increases. When the step voltage is applied, IPMC membrane bends toward anode side quickly and bends back gradually. The linear actuator has such a characteristic behavior. From (b), it is observed that the current increased sharply at the moment when the input voltage is applied, then it decreased exponentially. The peak value of current increases as input voltage does. Response with Loads Fig. 6 shows the response in step voltage with loads, step input voltage of 2.5 [V] was applied in loading. As the load is increased, displacement becomes small but the current does not almost change. Fig. 7 shows the response in step load. A constant load is added on after 1 [s]. Displacement becomes large in proportion to weight of loads. 3.3 Model Identification In order to know the capability of the linear actuator for robotic system, we identify the linear actuator as a mathematical model. IPMC actuators have been modeled in various ways in consideration of the physical and chemical phenomenon 18, 19), however, it is difficult to represent these models by systems of ordinary differential equation because of its complexity. In this paper, in consideration of model-based control, we identify the actuator as a linear time invariant model with a static nonlinearity from input-output data 20) using a subspace identification algorithm. % % % Figure 7: Response with various loads First, the model of the actuator is assumed to be represented by a system of Fig. 8. This model has two inputs and one output, and consists of two sub-system, P 1 (s) andp 2 (s), which are connected in series. P 1 (s) is a system with one input v, input voltage, and one output f 1 which is a force generated by electric stimuli. P 2 (s) is a system with one input force f 2 which is exerted to the actuator, and one output y, displacement of actuator. f 2 is assumed to be the difference between f 1 and f l. P (s) is defined as P 2 (s)p 1 (s). In order to conduct the simulations that the actuator is applied for a mechanical system, we need to identify both of sub-systems. It is difficult to add two inputs simultaneously, so we first assume that they are linear time invariant (LTI) systems, and identify them as following procedure: 1. Identification of P (s): Measure a response from input voltage v to displacement y, then compute the system P (s) = P 2 (s)p 1 (s) from input-output data using a subspace identification algorithm. 2. Identification of P 2 (s): Measure a response from load f l to displacement y, then compute the system P 2 (s) from input-output

f l v P 1 f 1 f 2 y P 2 b u h m h Figure 8: Block diagram m l m l g a θ 1 (x, y) θ 2 l = a b Figure 10: Experimental walking machine and its ideal model (a) displacement. u h # " f actuator r pulley (b) input voltage. Figure 9: Identification result with LTI model Figure 11: Actuator model of hip joint data using a subspace identification algorithm. 3. Computation of P 1 (s): Compute the system P 1 (s) asp 2 (s) 1 P (s). In procedure 1 and 2, we performed the system identification using N4SID function 21) in MATLAB. Fig. 9 shows a comparison between experimental result and simulated result using the identification model. It is shown that a quick transient of simulated result is nearly equal to the experimental one. But in part of slow decay, there is a little error between them. As the results of identification with linear approximate model, the characteristics of the actuator is captured, however, we have to identify the characteristics in consideration of nonlinear effects to obtain an appropriate model especially in large operating range. It is reasonable to assume that actuator has many nonlinear characteristics and one of the most influential factor is dynamics of load since the bending characteristics of the films highly depends on the load. In order to include the effect, we identified the actuator as Hammerstein model that consists of static nonlinearity and LTI system, where the static nonlinearity is followed by a LTI system in a cascade connection. Please see reference 8) for the detail. 4. Application to Biped Walker This section addresses an application of the actuator to a small-sized biped walking robot. 4.1 Model of Biped Robot In this paper, a simplest planar walking model, so-called compass-gait walker, is chosen as a target system. Fig. 10 shows the experimental walking machine and its ideal model, and Table 1 shows the values of the parameters. Because of the symmetric structure, the motion is constrained in a sagittal plane and the robot can be regarded as a two-legged robot. A heel-strike collision is assumed to be inelastic and without sliding in the simulation. The robot with suitable parameters can walk down a gentle slope without any driving force, and generates a stable periodic gait. Our biped model with the parameters has the nature of the passive dynamic walking. We consider the walking with actuators at hip joint on a slope or a level ground. To actuate the hip joint, contractive force of linear actuator transforms into torque by the structure as Fig. 11. We assume that one joint is actuated by a pair of linear actuators, then the hip joint of the robot is able to be rotated clockwise and counterclockwise as in this figure.

" " Table 1: Parameters of walking machine Actuator r m l leg mass 0.005 [kg] m h hip mass 0.010 [kg] total mass 0.020 [kg] a lower part of leg 0.05 [m] b upper part of leg 0.05 [m] l leg length 0.10 [m] r radius of hip pulley 0.004 [m] g gravity acceleration 9.81 [m/s 2 ] θ 1,θ 2 leg angle [rad] u h hip torque [N m] l m θ Figure 13: Experimental setup of swinging a pendulum #" % %! Figure 14: Experimental result (a) angular positions! # # " (b) input = = = ing cycle is so fast and the walking motion is generated with quick transient motion and small displacement of the actuator. Fig. 12 shows the simulated results of walking on a level ground. The number of units connected in parallel and series is set as 4 and 3, respectively, and the radius of pulley is set as r =0.004 [m]. In this simulation, we applied a square pulse as input voltage whose cycle is 0.48 [s] and whose amplitude is 2.5 [V]. From this results, it can be seen that the oneperiodic walking gait is generated and the walking cycle synchronizes with the cycle of the input signal. Figure 12: ground (c) stick diagram Simulated results of walking on a level 4.2 Numerical Simulation We conducted walking simulations using MATX 22) and we use LTI model as the actuator model since a walk- 4.3 Experiment We carry out the experiments of swinging a pendulum as preliminary tests for walking experiments. Fig. 14 shows the result. The weight of a small pendulum is 0.002 [kg] and a length is 0.1 [m]. We applied a square pulse as input voltages whose cycle is 0.5 [s] and whose amplitude is 3.0 [V]. In comparison with simulated results, the amplitude decreases slightly under the influence of frictions or modeling error, but it is confirmed that the motion is enough fast to realize a bipedal walking. 5. Conclusions We have discussed the development of linear actuator using IPMC materials in order to apply a soft actuator to biped walking robot. The elementary unit of the actuator was constructed, and its empirical model was identified from input-output data. Based on the empirical model, we demonstrated walking simulations of a

biped walking robot. Now we are planing to test the validity of the actuator using an experimental walking machine. In order to apply the artificial muscle actuator to a general robotic system, there exists a lot of problem such as limitation of output force, however, we think the mutual evolution of improvement of actuator technology and design of control system is important for further applications. References [1] T. McGeer, Passive dynamic walking, The Int. Journal of Robotics Research, Vol. 9, No. 2, pp. 62 82, 1990. [2] F. Asano et. al., A novel gait generation for biped walking robots based on mechanical energy constraint, Proc. of IEEE/RSJ Int. Conf. on Intelligent Robots and Systems(IROS), pp. 2637 2644, 2002. [3] F. Asano and M. Yamakita, Virtual gravity and coupling control for robotic gait synthesis, IEEE Trans. on Systems, Men and Cybernetics, Part A: Systems and Humans, Vol. 31, No. 6, pp. 737 745, 2001. [4] M. Yamakita et. al., Virtual coupling control for dynamic bipedal walking, Proc. of IEEE/RSJ Int. Conf. on Intelligent Robots and Systems(IROS), pp. 233 238, 2001. [5] D. W. Robinson et. al., Series elastic actuator development for a biomimetic walking robot, Proc. of IEEE/ASME International Conference on Advanced Intelligent Mechatronics(AIM), pp. 561 568, 1999. [6] M. Wisse and J. van Frankenhuyzen, Design and construction of MIKE; a 2D autonomous biped on passive dynamic walking, Proc. of Int. Symp. on Adaptive Morion of Animals and Machines(AMAM), 2003. [7] Y. Bar-Cohen, Electro-active polymers: current capabilities and challenges, Proc. of SPIE Int. Symp. on Smart Structures and Materials, EAPAD, 2002. [8] M. Yamakita et. al., Development of artificial muscle actuator using ionic polymer with its application to biped walking robots, Proc. of SPIE Int. Symposium on Smart Structures and Materials, EAPAD, 2003. [9] K. Oguro et al., Bending of an ion-conducting polymer film-electrode composite by an electric stimulus at low voltage, Journal of Micromachine Society, Vol. 5, pp. 27 30, 1992. (in Japanese) [10] S. Guo et al., Micro catheter system with active guide wire, Proc. of IEEE Int. Conf. on Robotics and Automation (ICRA), pp. 79 84, 1995. [11] S. Sawa et al., Polymer electrolyte actuator with gold electrodes, Proc. of SPIE Int. Symp. on Smart Structures and Materials, EAPAD, Vol. 3669, pp. 64 71, 1999. [12] EAMEX Corporation, http://www.eamex.co.jp/ [13] Y. Bar-Cohen et al., Challenges to the application of IPMC as actuators of planetary mechanisms, Proc. of SPIE Int. Symp. on Smart Structures and Materials, EAPAD, Vol. 3987, 2000. [14] S. Guo et al., Development of a new type of capsule micropump, Proc. of IEEE Int. Conf. on Robotics and Automation (ICRA), pp. 2171 2176, 1999. [15] S. Guo et al, Development of underwater microrobot using ICPF actuator, Proc. of IEEE Int. Conf. on Robotics and Automation (ICRA), pp. 1829 1834, 1998. [16] S. Tadokoro et al., Multi-DOF device for soft micromanipulation consisting of soft gel actuator elements, Proc. of IEEE Int. Conf. on Robotics and Automation (ICRA), pp. 2177 2182, 1999. [17] S. Tadokoro et al., Development of a distributed actuation device consisting of soft gel actuator elements, Proc. of IEEE Int. Conf. on Robotics and Automation (ICRA), pp. 2155 2160, 1998. [18] K. Asaka and K. Oguro, Bending of polyelectrolyte membrane platinum composites by electric stimuli Part II. Response kinetics, Journal of Electroanalytical Chemistry, 480, pp. 186 198, 2000. [19] S. Tadokoro et al., A dynamical model of ICPF actuator considering ion-induced lateral strain for molluskan robotics, Proc. of IEEE Int. Conf. on Robotics and Automation (ICRA), pp. 2010 2017, 2002. [20] K. Mallavarapu et al., Feedback control of the bending response of ionic polymer-metal composite actuators, Proc. of SPIE Int. Symposium on Smart Structures and Materials, EAPAD, Vol. 4329, pp. 301 310, 2001. [21] P. Van Overschee and B. De Moor, N4SID: Subspace algorithms for the identification of combined deterministic-stochastic systems, Automatica, 30(1), pp. 75 93, 1994. [22] M. Koga, Numerical Computation with MATX, Tokyo Denki University Press, 2000. (In Japanese)