Development of a Swimming Humanoid robot for Research of Human Swimming

Similar documents
Available online at ScienceDirect. Procedia Engineering 112 (2015 )

Free Swimming of the Swimming Humanoid Robot for the Crawl Stroke

Simulation analysis of the influence of breathing on the performance in breaststroke

Development of a Simulation Model for Swimming with Diving Fins

Comparison of Active Drag Using the MRT-Method and the MAD-System in Front Crawl Swimming

Body Stabilization of PDW toward Humanoid Walking

Journal of Biomechanical Science and Engineering

Development of Fish type Robot based on the Analysis of Swimming Motion of Bluefin Tuna Comparison between Tuna-type Fin and Rectangular Fin -

A NEW GOLF-SWING ROBOT MODEL UTILIZING SHAFT ELASTICITY

Analysis of stroke technique using acceleration sensor IC in freestyle swimming

1. A tendency to roll or heel when turning (a known and typically constant disturbance) 2. Motion induced by surface waves of certain frequencies.

Swimming practical examination support materials

Joint Torque Evaluation of Lower Limbs in Bicycle Pedaling

BODY FORM INFLUENCES ON THE DRAG EXPERIENCED BY JUNIOR SWIMMERS. Australia, Perth, Australia

Kinematic Differences between Set- and Jump-Shot Motions in Basketball

CRITERIA OF BOW-DIVING PHENOMENA FOR PLANING CRAFT

Humanoid Robots and biped locomotion. Contact: Egidio Falotico

Dynamic analysis and motion measurement of ski turns using inertial and force sensors

Generation of Robot Motion Based on Measurement of Human Movement. Susumu Sakano 1, Satoru Shoji 1

Motion Control of a Bipedal Walking Robot

Level 1 Stroke Performance Criteria

ZMP Trajectory Generation for Reduced Trunk Motions of Biped Robots

A Study on the Human Impulse Characteristics Of the Standing Shooting Posture

Gait Analysis of a Little Biped Robot. Received May 2015; accepted July 2015

ZSTT Team Description Paper for Humanoid size League of Robocup 2017

Video recording setup

WATER HYDRAULIC SYSTEM FOR HIGH SPEED CYLINDER DRIVE

SHUFFLE TURN OF HUMANOID ROBOT SIMULATION BASED ON EMG MEASUREMENT

Timing 1. The hips and shoulders rotate at the end of the catch to assist in the acceleration of the stroke.

SIMULTANEOUS RECORDINGS OF VELOCITY AND VIDEO DURING SWIMMING

Designing Diving Beetle Inspired Underwater Robot(D.BeeBot)

Kungl Tekniska Högskolan

Centre for Autonomous Systems

Deck Configuration and Jurisdiction for S&T Judges and CJ s during the CJ Role Playing:

The Mechanics of Modern BREASTSTROKE Swimming Dr Ralph Richards

Synchronised Swimming. Skill Level 1

Characteristics of ball impact on curve shot in soccer

OPTIMAL TRAJECTORY GENERATION OF COMPASS-GAIT BIPED BASED ON PASSIVE DYNAMIC WALKING

Naval Special Warfare Combat Side Stroke Guide

COMPARISON STUDY BETWEEN THE EFFICIENY OF THE START TECHNIQUES IN THE ROMANIAN COMPETITIVE SWIMMING

Twisting techniques used in freestyle aerial skiing

TECHNOLOGY FOR DECREASING ACTIVE DRAG AT THE MAXIMAL SWIMMING VELOCITY

DEVELOPMENT OF A FULL-SIZED BIPEDAL HUMANOID ROBOT UTILIZING SPRING ASSISTED PARALLEL FOUR-BAR LINKAGES WITH SYNCHRONIZED ACTUATION

Design, Fabrication and Analysis of Microcontroller Based Bipedal Walking Robot Vaidyanathan.V.T 1 and Sivaramakrishnan.R 2

Computer Aided Drafting, Design and Manufacturing Volume 26, Number 2, June 2016, Page 53. The design of exoskeleton lower limbs rehabilitation robot

Biomechanical Comparison of Two Different Kicks in Soccer

Mobile Robots (Legged) (Take class notes)

Dynamically stepping over large obstacle utilizing PSO optimization in the B4LC system

Skippy: Reaching for the Performance Envelope

Simulation of the Hybtor Robot

Swimming Stroke Mechanics

WATER HYDRAULIC HIGH SPEED SOLENOID VALVE AND ITS APPLICATION

Return this cover page with the test.

Biomechanical Analysis of Body Movement During Skiing Over Bumps

Robot motion by simultaneously wheel and leg propulsion

+ t1 t2 moment-time curves


An investigation of kinematic and kinetic variables for the description of prosthetic gait using the ENOCH system

Walking with coffee: when and why coffee spills

Friction properties of the face of a hand-held tennis racket

Investigation of Bio-Kinematic Elements of Three Point Shoot in Basketball

BIOMECHANICAL MOVEMENT

Controlling Walking Behavior of Passive Dynamic Walker utilizing Passive Joint Compliance

To Appear in ASME Journal of Biomechanical Engineering

Synchronised Swimming. Skill Level 1

Emergent walking stop using 3-D ZMP modification criteria map for humanoid robot

CHAPTER IV FINITE ELEMENT ANALYSIS OF THE KNEE JOINT WITHOUT A MEDICAL IMPLANT

Impact Points and Their Effect on Trajectory in Soccer

Ball impact dynamics of knuckling shot in soccer

Sample Solution for Problem 1.a

ABSTRACT PATTERNS USING 3D-DYNAMIC MODELING. Kaustav Nandy, Master of Science, Department of Electrical And Computer Engineering

Testing of hull drag for a sailboat

Toward a Human-like Biped Robot with Compliant Legs

Basketball free-throw rebound motions

Abstract 1. INTRODUCTION

Biomechanics and Models of Locomotion

Key words: biomechanics, injury, technique, measurement, strength, evaluation

ITF Coaches Education Programme Biomechanics of the forehand stroke

Optimization of an off-road bicycle with four-bar linkage rear suspension

Development of Turtle-like Submergence Vehicle

Complex movement patterns of a bipedal walk

Swimming Breaststroke Checklist Marion Alexander, Yumeng Li, Adam Toffan, Biomechanics Lab, U of Manitoba

Tether-based Robot Locomotion Experiments in REX-J mission

REPORT DOCUMENTATION PAGE

Development of Technology to Estimate the Flow Field around Ship Hull Considering Wave Making and Propeller Rotating Effects

Available online at ScienceDirect. Procedia Engineering 112 (2015 ) 40 45

Aerodynamic Measures for the Vortex-induced Vibration of π-shape Composite Girder in Cable-stayed Bridge

INTERACTION OF STEP LENGTH AND STEP RATE DURING SPRINT RUNNING

A COMPARISON OF SELECTED BIOMECHANICAL PARAMETERS OF FRONT ROW SPIKE BETWEEN SHORT SET AND HIGH SET BALL

International Journal of Technical Research and Applications e-issn: , Volume 4, Issue 3 (May-June, 2016), PP.

Spring Locomotion Concepts. Roland Siegwart, Margarita Chli, Martin Rufli. ASL Autonomous Systems Lab. Autonomous Mobile Robots

Development and analysis of a novel add-on drive mechanism for motorizing a manual wheelchair and its step climbing strategy

Speed Control System Design in Bicycle Robot by Low Power Method. Abstract

PERCEPTIVE ROBOT MOVING IN 3D WORLD. D.E- Okhotsimsky, A.K. Platonov USSR

2014 CCAA STROKE and TURN TEST

Biomechanical analysis of the medalists in the 10,000 metres at the 2007 World Championships in Athletics

Wind turbine Varying blade length with wind speed

Freestyle A swimmer who can successfully stay on top of the water can have more efficiency because there is less surface area for the water to grab on

Design Of Humanoid s Lower Limb Model

ROSE-HULMAN INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering. Mini-project 3 Tennis ball launcher

Transcription:

Development of a Swimming Humanoid robot for Research of Human Swimming Changhyun Chung a and Motomu Nakashima b,* a Department of Mechanical and Environmental Informatics, Tokyo Institute of Technology, Japan b Department of Mechanical and Control Engineering, Tokyo Institute of Technology, Japan Abstract In this study, the swimming humanoid robot SWUMANOID for research of human swimming was developed. It was half the size of a real human and had the same body proportions and appearance; mass distribution was considered, as well. It was actuated by 2 motors, which were waterproofed and compactly designed. Swimming motions were determined based on human swimming movements. To realize the swimming arm stroke, a joint imitating the human s scapular retraction was installed and the methodology to realize the swimming motion with that joint was established. Finally, the efficiency of the design was validated by simulation and representatively, the crawl stroke was successfully realized in the water tank. Index Terms Swimming, Humanoid, Kinematics, SWUMANOID, SWUM I. INTRODUCTION Competitive swimming is one of the most popular sports. It has been played all over the world including in the Olympic Games. However, its mechanics have not been fully clarified yet since it is an extremely complicated phenomenon in which a complex human body moves unsteadily with many degrees-of-freedom in the three-dimensional water flow. To clarify the relationship between the swimming motion and the propulsive force, researchers have made various attempts to measure the fluid forces acting on a swimmer. The first approach was an experiment involving a human subject [1-3]. However, this method has problems with insufficient repeatability, physical fatigue of the subject, and difficulty in installing sensors in the subject. In a different way, some researchers used physical models instead of human subjects. In previous studies, Lauder and Dabnichki developed a full-scale mechanical arm to estimate propulsive forces. It had one degree of freedom (DOF) at the shoulder and was * Corresponding author: 2-12-1-I3-9 Ookayama, Meguro-ku, Tokyo, 152-8552, Japan, E-mail: motomu@mech.titech.ac.jp used to simulate a single plane underwater rotation at four elbow configurations [4]. Sidelnik and Young made a motor-driven arm model which could perform 2-DOF motion, and measured the fluid forces in the unsteady condition [5]. Nakashima and Takahashi conducted experiments using an underwater robot arm which had 5-DOF and could perform a swimming motion which was much more similar to that of a human than in previous studies [6, 7]. However, the effects of the trunk motion and interaction among four limbs were still ignored since the robot arm could represent only a single limb. In order to investigate such effects, a fullbody swimming humanoid robot will be necessary. In addition to this significance, the swimming humanoid robot has many possibilities in other applications such as education, entertainment, exploration, sea rescue and so on. In spite of its usefulness, only a few researchers studied the humanoid robot for underwater use. Oya and Suzuki proposed a small diver type humanoid robot to estimate water depth by the robot s behavior [8]. Li et al. developed a swimming robot using the flutter kick for underwater use [9]. Nakashima and Kobayashi developed smallsized swimming humanoids [1]. All of them were remodeled from a land-based robot for underwater use and had waterproof clothes similar to diving suits. Such thick clothes restricted the movement of the robot, making it difficult to perform the right movement. The former two robots moved in the water by the flutter kick only. The latter one swam with arms, although there were problems, such as that the waterproofing and the motion was unnatural due to the limitation of the degrees of freedom and the range of the joint angles. The objective of this study was to develop a full-body swimming humanoid robot for research of human swimming. The developed robot was named SWUMANOID (SWimming humanoid). It has a detailed human s body shape and 21-DOF for the joints in order to reproduce the complicated swimming motions with high fidelity. It also has a potential to be an exemplar of the swimming humanoid robots for the practical applications since it does not have any 19

waterproof clothes and has sufficient DOF as well as the range of motion. The paper is organized as follows: Section 2 describes the basic design of the robot, Section 3 describes the methodology to realize swimming motions, Section 4 describes the experiment of the crawl stroke motion, and Section 5 presents the conclusion and future work. II. BASIC DESIGN OF THE ROBOT A. Design concepts The design concepts of SWUMANOID are as follows: a) Waterproof design b) Imitation of human body geometry : Appearance, body proportions, center of mass and specific gravity c) Realization of swimming motions d) Half the size of a real human Since SWUMANOID is the swimming test model substituting a human subject, the design concepts a)-c) are obvious. The robot size was determined to be half the size of a real human in order to make development easier and to be more convenient when it is operated in experiments. Figure 1 shows SWUMANOID which was developed based on the design concepts a) to d), as well as the scanned body data of an elite swimmer. As shown in Table I, SWUMANOID is 89 [mm] in height and 5.5 [kg] in weight including the battery. The height was determined by the design concept d) and the weight was determined by the design concept b), which considered human s specific gravity. To realize the design concept c), SWUMANOID has a total of 21-DOF but the leg doesn t have sufficient degrees of freedom to perform the breaststroke kick, so in this study, realization of the breaststroke was excluded. B. Configuration of degrees of freedom Figure 2 shows the positions of the actuators and Items Dimensions Weight inc. battery Degree of freedom Microprocessor Battery Table 1 Specifications of SWUMANOID Specifications Height: 89 mm Width: 27 mm Depth: 1 mm 5.5 kg Total: 21 DOF Head (Free joint): 1 DOF Arm: 2 Arms x 6 DOF Waist: 2 DOF Leg: 2 Legs x 3 DOF ATmega 2561 (Atmel Corp.) Li-Po.1V 44 mha joints. All actuators were connected to the joints directly except for the head and ankle joints. The head joint was a free joint in order to keep the head down during swimming motions. The ankle joint was actuated using a four-bar linkage because there was no space for the installation of an actuator in the ankle. Each arm had 6-DOF. The elbow had 2-DOF and the shoulder had 4-DOF. Three of the shoulder joints were for the roll-pitch-yaw motion and one was for the scapular retraction (scapular joint in Figure 2). The reason for installation of the scapular joint is described in section 3.2. Each leg had 3-DOF for extension and flexion of the hip, knee and ankle joints. To produce more realistic motions, the ankle joint tilted about 2 degrees as shown in Figure 2. C. Actuators and control system The robot was driven by 2 servo actuators called DYNAMIXEL (RX-24F, Robotis Co.). They were linked by a daisy chain, so wiring was simple and we saved many steps in waterproofing the wires compared to normal RC servo motors. However, the motor was Fig. 1 Scanned body data of an elite swimmer compared with a developed swimming humanoid robot Fig. 2 DOF configuration

(a) DYNAMIXEL RX-24F (b) Modified motor (a) Waterproofing of the axis Fig. 4 Waterproofing of the motor (b) Waterproofing of the board were not enough space to install the o-ring and screws. Finally, the body cases produced by the rapid prototyping machine were coated with a waterproofing chemical. (c) Forearm assembly with a modified motor Fig. 3 The features of a developed compact waterproof motor too large to be installed inside of the robot cases as shown in Figure 3 (a). Therefore, we newly designed a motor case and gear housing to be compact and waterproof. Figure 3 (b) shows the modified motor. The one-sided output axis was centered, and unnecessary cases were removed. The modified motors were successfully installed in narrow parts, such as the longitudinal axis of the forearm as shown in Figure3 (c). The installed motors were controlled by a robot controller (CM-7, Robotis Co.) and it was possible to monitor conditions of the motor such as angles, temperature, voltage, torque, and so on through a personal computer. D. Waterproofing Waterproofing was one of the most important and essential parts to developing an underwater robot. In the previous research of an underwater humanoid robot, most robots wore waterproofing clothes and were sealed by silicone for the waterproofing [8-1]. It was an easy and simple method to waterproof a robot body, but this method was different from our research objective. First, it is difficult to keep a uniform body shape and it would exert a bad effect on the measurement of fluid force. Second, it makes it difficult to measure body volume and density. In this study, keeping a constant specific gravity was important not only for measurement of fluid force but also for a free swimming test in the future. In order to maintain a uniform body shape and similar specific gravity with that of a human, we utilized three types of waterproof methods as follows: First, we modified motor cases for the waterproofing. The rotational axis of the motor was sealed by a quad-ring and o-ring as shown in Figure 4 (a) and the control board of the motor was sealed by silicon as shown in Figure 4 (b), in case of leakage from the outer case. Second, the motors and electric wires were fixed and sealed by silicone and an epoxy bond because there E. Imitation of human body geometry Appearance, body proportions, center of mass, and specific gravity were considered to imitate a human subject. To mimic the real human body, the outer cases of the humanoid robot were made based on the scanned body data of a world-class elite male swimmer. Each case had a half length of real body parts and is hollow inside with 2~3 mm thickness to contain robot parts. After assembling, cases are sealed for waterproofing. They were built by a rapid prototyping machine with ABS material. The scanned body data of a swimmer and a photograph of the developed humanoid robot are shown in Figure 1. Every driving part of the robot was installed inside to keep the scanned body shape. Actuators were redesigned compactly to be installed inside. The controller and battery were also inside of the robot, so that the humanoid robot obtained a similar appearance to a human and could swim without any connection from the ground. The weight and volume of the robot s separate parts are listed in Table II. The total weight of the robot is approximately 5.5 kg and the specific gravity is about.82. In the experiment, the initial conditions, the center of gravity and the total weight of the robot was adjusted by adding more weights inside the robot. The robot could maintain proper buoyancy and balance without any other buoyancy equipments in the water, suggesting Table 2 Weight and volume of a humanoid robot Parts Weight(g) Added weight(g) Volume ( cm3 ) Specific gravity Head & Neck 264 59 1.85 Upper arm*2 57 72.79 Forearm*2 25 28.89 Breast 15 37 21.68 Waist 36 353 1.2 Thigh*2 56 682 14.89 Shank*2 67 7.96 Foot*2 82 6 2.71 Total 386 172 6753.82 1

a possibility of an autonomous underwater humanoid robot. III. METHODOLOGY TO REALIZE SWIMMING MOTIONS A. Overview of Swimming motion generation The method of swimming motion generation is schematically shown in Figure 5. First, we analyzed movies of model swimming by a human subject and determined joint angles of a simulation model using a swimming human simulator, SWUM []. Using SWUM, we could check the fluid forces acting on the robot as well as the joint torques when the robot performed swimming motions. In SWUM, the model of the swimmer s body was divided into 21 cylindrical segments, and one stroke cycle was divided into 18 motion frames. We calculated the desired positions and orientations of the robot arm, solving forward kinematics at each motion frame. Also, the joint angles of the robot arm were determined by applying the coordinate transformation and solving the inverse kinematics. On the other hand, leg motions were simply realized because the flutter kick and butterfly kick were expressed by two dimensional motions. The determined joint angles were simulated with a robot model and finally performed by the developed humanoid robot. B. Usage of scapular joint Arm motions during swimming are composed of the underwater stroke and recovery stroke. Swimmers perform these two strokes continuously. However, such a continuous stroke motion cannot be realized with a general robot arm which is constructed of a threerotation sequence of pitch, yaw and roll. This problem is schematically explained in Figure 6 (a). Stroke motion is mainly generated by the pitch joint (Shoulder Joint 1 in Figure 2). The angle of Shoulder Joint 1 is increased from to 18 degrees for the underwater stroke as shown in the left of Figure 6 (a). Then, it has to return to degrees during the recovery stroke as shown in the right of Figure 6 (a) in order to perform the next underwater stroke continuously. Therefore, the recovery stroke cannot become the over arm motion, as shown in the left of Figure 6 (b). In order to avoid this problem, the scapular joint was installed on the robot s arm based on the idea that the human scapular retraction motion makes it easy to perform an over arm motion. With the scapular joint, the recovery stroke can become the natural over arm motion as shown in the right of Figure 6 (b). In order to use the scapular joint systematically, we limited the range of the retraction motion to the over arm motions, such as the recovery motion in the crawl stroke. It was reasonable considering the motion of humans. Afterwards, we recalculate the position and orientation of the elbow coordinate considering the moved position of the shoulder joint by the scapular joint. From the recalculated goal position of the elbow, we could calculate the rest motor angles through the inverse kinematics. The details are described in the following sections. C. Inverse Kinematics Accompanying Scapular Joint The joint motion in the present study was determined based on standard six-beat crawl stroke of SWUM. Figure 7 shows the reference and link coordinate systems of the 6-DOF robot arm using the Denavit- Hartenberg convention. The values of the kinematic Fig. 5 Overview of swimming motion generation 2

(a) Rotation of Shoulder Joint 1 for underwater and recovery strokes Fig. 7 Coordinate system of the 6-DOF robot arm T elbow _ SWUM 21 12 22 32 13 23 33 p x p y p z 1 (2) (b) Comparison of recovery strokes without/with the scapular joint Fig. 6 Stroke motions of the robot parameters are listed in Table 3, where a, b and c are the lengths between the shoulder joints, and l 1, l 2 and l 3 are the arm s lengths. They correspond to the home position pictured in the link coordinate diagram of Figure 7. The position and orientation of the wrist with respect to the base frame is expressed as: b 1 2 3 Twrist T 1T 2T 3T 4T r r12 r13 px r21 r22 r23 py r r32 r33 pz 1 5 T T 4 5 w (1) where [r r 21 r ] T, [r 12 r 22 r 32 ] T and [r 13 r 23 r 33 ] T are the unit vectors of the frame attached to the wrist, and [p x p y p z ] T is the position vector of the origin of this frame with respect to the origin of the base frame (x b y b z b ) (Figure 7). The angle of the scapular joint in the over arm motion was determined as the upper arm angle with respect to the horizontal line, which is shown in Figure 8 (a). This angle can be calculated from the elbow transformation matrix. The transformation matrix is expressed as: where [r 13 r 23 r 33 ] T is the unit vector of z 3 in Figure 9 and it represents the upper arm s direction. Therefore, the angle of the upper arm with respect to the horizontal line or with respect to axis y sb in Figure 9 is obtained by followed equation: Table 3 Denavit-Hartenberge parameters of the arm i α i-1 a i-1 d i i 1-9 a b 1-9 2-9 c 2-9 3-9 l 1 3 4-9 l 2 4 5 9 5 6 l 3 a i the distance from Z i to Z i+1 measured along Xi; α i the angle from Z i to Z i+1 measured about Xi; d i the distance from X i-1 to X i measured along Zi; and i the angle from X i-1 to X i measured about Zi. (a) Determination of the angle of scapular joint from SWUM (b) Application of the angle of scapular joint in the robot model Fig. 8 Determination of the angle of scapular joint ( ) of the robot from the simulation model 3

T elbow _ SWUM 21 12 22 32 + 13 23 33 p x p y p z 1 c s + a c c c + a s 1 (5) Fig. 9 Relationship of SWUM s base frame and the robot s base frame 1 13 tan ( ) (3) Also, around the beginning and end of the overarm motion, we set the angle of scapular joint so that it varied gradually. As shown in Figure 9, the base frame of the simulation model {sb} and the base frame of the robot model {b} are different. The frame {sb} is located at the end of the shoulder and the frame {b} is located in the scapular joint because there is no scapular joint in the simulation model. Therefore, the original desired frame of the wrist, which was calculated from the simulation model, was derived from frame {sb}. For the robot, the origin of the base frame of the simulation model is the same as initial position of robot frame {2} in Figure 9, but the frame {2} moves depending on movement of at each motion frame. The difference between frame {sb} and frame {2} is given by equation (4). To simplify the notation, s i, c i, s i-9 and c i-9 are sin i, cos i, sin( i -9 ), and cos( i -9 ) respectively. cc1 9c2 9 + ss b 1 sc1 9c2 9 cs T 1T2T s1 9c2 9 2 9 2 9 c c s c 23 1 9 2 9 1 9 2 9 s s s s 1 9 2 9 + s c 2 9 c c 2 9 c s 1 9 s s 1 9 c 1 9 cs + ac cc + as 1 (4) We could not calculate the joint angles of the robot arm because the origin of frame {2} was not consistent with the origin of frame {sb} but if we moved the origin of frame {sb} to the origin of frame {2}, we could calculate the joint angles using the motion data of the simulation. Considering the difference, the desired position and orientation were modified as follows in equation (5): T elbow b 1 2 T 1T 2T3T As a result, the shoulder s position would have a small error because of the difference of the base frames. However, we could calculate the joint angles of the robot from the motion data of simulation, as described in the following section. To determine the arm joint angles for a given position and orientation of the wrist, we solved the inverse kinematics problems. The relationship between the elbow and the reference base frame is given by equation (5). The solution starts with the matrix equation (6): b 1 1 2 T Telbow 1T 2T3T L1 R1 (6) The equations L1(1,3)R1(1,3) and L1(3,3)R1(3,3) of the elements of the matrices L1 and R1 can be used to solve for 1. This yields: 1 33 1 9 tan ( ) (7) c + s 13 23 The equations L1(2,2)R1(2,2) and L1(2,1)R1(2,1) of the matrices L1 and R1 provide the solution for 3 : 1 s 12 cr22 3 tan ( ) (8) s + c 22 The equation L1(2,1)R1(2,1) and L1(2,3)R1(2,3) of the matrices L1 and R1 provide the solution for 2 : 1 s c 21 2 9 tan ( ) (9) ( s c ) c 1 1 1 b 1 2 T 1T T 13 23 T Twrist r12 r13 px 2 3 4 5 r21 r22 r23 py 3T 4T 5T wt r32 r33 pz 1 2 1 3 4 5 3T T 4T 5T wt L2 R2 3 (1) () 4

The equations L2(1,3)R2(1,3), L2(3,3)R2(1,3) and L2(2,1)R2(2,1) ), L2(2,2)R2(2,2) of the matrices L2 and R2 provide the solution for 4, and 5 : 1 c3 13 s3 33 4 tan ( ) 23 (12) 1 s3 c3 5 tan ( ) s c (13) 3 12 3 32 If multiple values were obtained as solutions, we determined the value considering the range of motion and natural pose of a real human. The joint angles of the shoulder were simulated in the crawl stroke, backstroke and butterfly stroke. In order to validate and confirm the efficiency of the scapular joint for the retraction motion ( ), we simulated two conditions. First, we did not use the scapular joint and performed stroke motions with a 5-DOF robot arm, according to the results of the conventional inverse kinematics. Next, we used the scapular joint as explained above, corresponding to Figure 6 (b). Figure 1 shows the joint angles of the shoulder. The graphs in Figure 1 (a)-(c) are without the scapular joint, while those in Figure 1 (d)-(f) are with the scapular joint. As shown in Figure 1 (a)-(c), when there is no scapular joint, the Shoulder Joints 1 and 2 switch their positions (time frame to 1 and 6 to 7 in crawl, 17 to and 9 to 1 in back, and 9 to 1 and 17 to in butterfly). On the other hand, as shown in Figure 1 (d)-(f), when the scapular joint moves, the motion of Shoulder Joint 1 was continuous and the range of motion of the Joint 2 became significantly smaller. Therefore, we realized swimming stroke motion sufficiently similar to the human s one using the scapular joint. Due to the simple structure of the robot shoulder, kinematics of the robot arm was analyzed clearly and the methodology to use the joint was established. IV. EXPERIMENT A. Experimental Conditions We realized the crawl stroke motions in the water representatively among the four styles of swimming because the arm stroke motion of the crawl is the most complicated and fundamental. The swimming humanoid robot was connected to the four rods as shown in Figure, and the rolling motion was generated by the driving mechanism [13]. We realized the crawl stroke motion on the water at 8.9 sec per cycle. The speed was quite slow and it was about 25% of our target speed because the objective of this experiment was to check that the robot could realize the stroke motion clearly. B. Results and Discussion Figure 12 shows the desired and generated position trajectories in the Cartesian space during the execution of the crawl stroke motion. The commanded joint angles were derived from the crawl stroke data of simulation through the proposed kinematics method. The robot was controlled in open loop because it was enough to produce correct stroke motion as shown in Fig. 12. Figure 13 shows a series of images during the crawl stroke. The swimming humanoid robot performed the crawl stroke naturally. The video of the experiment is shown in Movie S1. The results show the potential of the experimental system for measurement of fluid force during swimming using a swimming humanoid robot. Also, we confirmed that the robot could generate swimming motions using the proposed methodology. It means we can link the developed simulation model, SWUM and the robot model, SWUMANOID. We could simulate various underwater motions and realize with the robot as well using the proposed kinematic methodology. It could make it easier to develop underwater locomotion of the robot and could be a base of research of humanoid (a) Crawl stroke without scapular joint (b) Back stroke without scapular joint (c) Butterfly without scapular joint (d) Crawl stroke with scapular joint (e) Back stroke with scapular joint (f) Butterfly with scapular joint Fig. 1 Joint angles of the right shoulder 5

Fig. Connection between the swimming humanoid robot and the driving mechanism for the underwater test robot for underwater use. Note that some awkwardness of the motion can be seen in Movie S1. The reasons for this may be that the commanded joint angles were somewhat different from the original human s joint angles due to the measurement error, and that the 18 frames for one stroke cycle were not sufficient to represent the human s motion. V. CONCLUSION AND FUTURE WORK In this study, the swimming humanoid robot, named SWUMANOID, was developed. The appearance and body properties were considered for the research of human swimming. Also, a compact waterproofing motor was developed and installed in the humanoid robot. To realize the human swimming motion, we added a scapular joint for the retraction motion and also developed an inverse kinematics model to determine the angle of the retraction joint; this model can be used to determine the angle of other joints. The efficiency of the scapular joint was validated by simulation of the crawl, back and butterfly strokes, and by experimentation of the crawl stroke. Measurement of the propulsive force in various swimming motions, validation of the simulation, and a free swimming test will be our future research goal. VI. ACKNOWLEDGMENTS This work was supported by a Grant-in-Aid for Scientific Research ((B) 233229) from the Japan Society for the Promotion of Science. REFERENCES [1] A.P Hollander, G. De Groot, G.J van Ingen Schenau, H.M Toussaint, H. De Best, W. Peeters, A. Meulemans and A.W Schreurs, Measurement of active drag during crawl arm stroke swimming, Journal of Sports Sciences 4(1), 1986, pp.21-3. [2] S.V Kolmogorov, O.A Duplishcheva, Active drag, useful mechanical power output and hydrodynamic force coefficient in different swimming strokes at maximal velocity, Journal of Biomechanics 25(3), 1992, pp.1-8. (a) Angles of scapular joint ( ) (b) Angles of shoulder joint 1 ( 1) (c) Angles of shoulder joint 2 ( 2) (d) Angles of shoulder joint 3 ( 3) (e) Angles of elbow joint 1 ( 4) (f) Angles of elbow joint 2 ( 5) (g) Angles of hip joint (h) Angles of knee joint (i) Angles of ankle joint Fig. 12 Joint angles of the right arm and right leg 6

Fig. 13 Snapshots of SWUMANOID in swimming experiments [3] H. Takagi and R. Sanders, Measurement of propulsion by the hand during competitive swimming, The Engineering of Sport 4 (Eds. S. Ujihashi and S.J. Haake): Blackwell Publishing, 22, pp.6-637. [4] M.A Lauder, P. Dabnichki, Estimating propulsive forces-sink or swim? Journal of Biomechanics 38(1), 25, pp.1984-199. [5] N.O Sidelnik and B.W Young, Optimising the freestyle swimming stroke: the effect of finger spread, Sports Engineering 9, 26, pp.129-135. [6] M. Nakashima and A. Takahashi A, Clarification of Unsteady Fluid Forces Acting on Limbs in Swimming Using an Underwater Robot Arm (Development of an Underwater Robot Arm and Measurement of Fluid Forces), Journal of Fluid Science and Technology 7(1), 212, pp.1-3. [7] M. Nakashima and A. Takahashi, Clarification of Unsteady Fluid Forces Acting on Limbs in Swimming Using an Underwater Robot Arm. Journal of Fluid Science and Technology 7(1), 212, pp.4-128. [8] H. Oya and K. Suzuki, Behavior Based Water Depth Estimation for Diver Type Small Humanoid Robot, Proceedings of the 4th International Conference on Autonomous Robots and Agents, Feb. 1-12, 29, Wellington, New Zealand, pp.43-435. [9] Y. Li, E. Shimizu and M. Ito, Development of a humanoid robot for underwater use, Journal of Artificial Life and Robotics, Vol. 15, Issue 4, December 21, pp.551-554. [1] M. Nakashima and H. Kobayashi, Development of small-sized swimming humanoids, in The Engineering of Sport 7(Eds. M. Estivalet and P. Brisson). Vol. 2. Paris: Springer-Verlag, 28, pp. 33-39. [] M. Nakashima, K. Satou and Y. Miura, Development of swimming human simulation model considering rigid body dynamics and unsteady fluid force for whole body, Journal of Fluids Science and Technology 2(1), 27, pp.56-67. [12] http://www2.edu.ipa.go.jp/gz/h1swim/h1kn2/h1cr2 /h1cr22.mpg [13] M. Nakashima, Y. Ejiri, Measurement and Modeling of Unsteady Fluid Force Acting on the Trunk of a Swimmer Using a Swimming Mannequin Robot, Journal of Fluid Science and Technology 7(1), 212, pp.-24 7