Acta Polytechnica CTU Proceedings doi:10.14311/APP.2017.7.0053 Acta Polytechnica CTU Proceedings 7:53–57, 2017 © Czech Technical University in Prague, 2017 available online at http://ojs.cvut.cz/ojs/index.php/app RIGHTING OF CHINESE MITTEN CRABS (ERIOCHEIR SINENSIS) AND THEIR MODELS Tobias Riphausa,∗, Florian Hoffmannb, Susanna Labischb a International Degree Course in Biomimetics, City University of Applied Science Bremen, Germany b Biomimetics-Innovation-Centre (B-I-C), City University of Applied Science Bremen, Germany ∗ corresponding author: T. Riphaus via slabisch@bionik.hs-bremen.de Abstract. The usage of unmanned underwater vehicles for marine tasks is continuously growing and bioinspired stabilizing systems shall help them to gain and keep a stable position during work. Therefore the righting maneuver of E. sinensis has been studied. These crabs are able to perform a 180◦-rotation with an angular velocity of 4.30 s−1 when falling underwater from a supine starting position. High-speed particle image velocimetry has shown, that propulsive forces with a peak of 0.021 ± 0.001 N were produced by the hind legs to initiate and stop the rotation. In a numerical multibody simulation a constant force of 0.009 N acting for 0.2 s leads to the same rotation. In order to prove this mechanism, it was implemented into a robotic system. Its mean density of 1.15 g/cm3 deviates not more than 4% from the biological and numerical models. It can complete a 180◦-turn within 1.03 ± 0.12 s with a rotational velocity of up to 4.25 s−1. Keywords: biomimetics, righting maneuver, robot stabilization. 1. Introduction Since the first success of unmanned underwater ve- hicles (UUVs) their field of operation is constantly spreading. In spite of all the technology, the opti- mization of their stabilizing behavior is still a crucial topic in contemporary research. The reason why even modern UUVs are facing problems in holding stable positions, to carry out their tasks or turn around quickly once they have lost their preferred position, is the kind of drive that is mainly used. Propeller and pump-jets are advantageous for high-speed loco- motion, but too expensive for the use in a seldom occurring event. To gain new ideas for solving such problems more and more research groups are investigating biological models. First successful terrestrial applications are the amphibious robot RoboCrab by Krummel et al., which adapts the skills of horseshoe crabs within the surf zone, or the four-wheeled autonomous robot by Chang- Siu et al., which is capable of mid-air reorienting during free fall modeled on geckoes (Hemidactylus platyurus) [1][2]. But the usage of biological models for improving aquatic robots is poorly documented. Thus, the re- search group DFU of the Department of Biomimetics in Bremen started investigations on different benthic marine animals. The aim was to search for a strategy that allows a quick turning or efficient reorientation manoeuvre which uses neither a propeller nor is driven by a permanent rotation. Especially promising seemed the crustacean order Decapoda. Therefore the species Eriocheir sinensis has been studied in terms of their mid-water righting behavior [3][4]. Besides the mor- phological biomechanical analysis, also the necessary Figure 1. Individual of E. sinensis (A, [3] modi- fied), its numerical model (B) [4] and the adapted robot O-Typus (C) with an insight into its electronic devices (D) [5]. hydrodynamic forces have to be determined to show, that such a strategy is transferable to a technical application. 2. Material and Methods 2.1. Eriocheir sinensis - biological model In total 52 rotations of 10 individuals of E. sinensis were analysed. The rotations are differentiated in the kinematic approach with 39 analysed rotations (28/11, initialising/stopping rotation) and for the PIV (Particle Image Velocimetry) analysis with 13 analysed rotations (9/4, initialising/stopping rotation). 53 http://dx.doi.org/10.14311/APP.2017.7.0053 http://ojs.cvut.cz/ojs/index.php/app T. Riphaus, F. Hoffmann, S. Labisch Acta Polytechnica CTU Proceedings 2.1.1. Kinematic analysis To observe the natural righting behavior of E. sinensis (Figure 1A), the crabs were held in a supine starting position at the height of 0.60 m within a freshwater tank with a ground area of 0.6 × 0.6 m and a high of 0.8 m and released by starting device. The body movement was high speed recorded (250 fps) to be able to manually track two points on the body and the tip of the 5th leg. Furthermore falling speed, body angle and leg trajectories were plotted and the speed of the leg tip calculated. The leg induced force (Fd) during the strike was calculated by Fd = 1 2 ρ v2 c cd A (1) where ρ = water density, vc = leg centroid velocity, A = area of 5th legs and cd = drag coefficient from static force balance measurements. Centroid velocity was calculated by: vc = vt Lbc Lbt (2) where vc = velocity of centroid, vt = velocity of leg tip, Lbc = length leg base to leg centroid and Lbt = length leg base to leg tip. Lbt and Lbc were experimental determined on 15 legs of 10 dead E. sinensis. Centroid x/y coordinates were calculated with an integrated function of ImageJ from the leg area. The resulting factor for Lbc / Lbt was calculated as 0.44 ± 0.016. 2.1.2. Piv measurements and force calculation The tank of the kinematic study was also used for the PIV measurements. Crabs were released by hand or the automated release fixture. The laser sheet was placed vertically in the middle of the tank illu- minating a surface of approx. 0.3 x 0.3 m. Flow was visualized by Polyamid Seeding particles (ø 57 µm, density: 1.03 g/cm3). Images were taken by a high speed camera (Photron APX-RS) at an interframing time of 3.75 ms (266 fps). Crabs were released at different distances to the lasersheet depending on the leg length. We used falling sequences where the last leg was cutting the lasersheet at positions from 50 to 80% of the leg. By measuring the bound circulation forces around the leg the generated propulsion force was calculated by using the Kutta-Joukowski theorem. Circulation Γ was measured from the beginning to the end of the leg stroke in order to get a temporal resolution of the circulation. In PIVLab [6] under MATLAB (The MathWorks Inc.) the PIV images were analysed and the spanwise circulation (circulation around the spanwise axis) was derived by a loop integral of the tangent velocity: Γ = ∮ Svort vt dSvort = ∫ Avort ωZ dAvort (3) where Svort = circular path around the vortex core, vt = tangential velocity, Avort = area of the vortex core and ωZ = spanwise vorticity. The vortex ring momentum calculation model was used to estimate an average propulsive force [7][8]. The average propulsive force of the leg stroke was herewith calculated by F = 2 Γ ρ ν b (4) where Γ = circulation, ρ = fluid density, ν = fluid velocity and b = leg length. The force results were multiplied by two, to get the force of both legs. 2.2. Simulated Model For the numerical simulation of the physical behavior the software ADAMS (rel. 2013.1) was used. It allows to describe the limb movement as a motion of mechani- cally joint rigid bodies, builds the relevant differential equations of motion and solves them numerically. Basis for the multi-body-analysis is experimental gained data of an individual crab showing the dimen- sions, the mass, the density and the center of mass (COM) of each limb. All body-parts are connected with rigid joints. In order to simulate the correct falling speed in water the gravity was reduced to 1.09 m/s2. During movement the crabs hold their front legs close to the carapax and stretching out only the hind legs. Therefore the main body can easily be simplified to a cylinder. With a diameter of 6.00 cm and a height of 2.75 cm for carapax and front legs such a cylinder will have an average density of 1.12 g/cm3 which leads to a mass of 87 g. This data fits the measured crab very well. The hind leg pair is modelled with the merus/ischium, the propodus/carpus and the dactylus in their measured length and density. The parts of the leg are conical so that the COM of these parts matches also (Figure 1B). The acting forces are applied on the joint between dactylus and propodus with a constant value for a given duration (thick arrows in Fig 1B). 2.3. O-Typus - Technical Model 2.3.1. Setup To confirm whether the observed mechanism is cor- rectly understood, the basic model is implemented in the purpose-built robot O-Typus (Figure 1C+D) [5]. Its chassis is congruent with the carapace of an individual of E. sinensis and scaled 2:1, to house drive units and the sensor system. The actuating elements consist of two waterproof HK15328D servo motors (max torque: 12,8 kg.cm) with a 6V external power supply (2250 mA). Their control is based on the PWM-signals from an Ar- duino Uno. The position tracking is enabled by the MPU6050 accelerometer. The real time operation sends the current Yaw-/ Pitch-/ Roll-Angles (YPR) to the Arduino Uno to log the data (~27 Hz). 54 vol. 7/2017 Righting of Chinese Mitten Crabs (Eriocheir sinensis) and their Models Figure 2. A. Stereotypical turning maneuvers ob- served for E. sinensis after a virtually moment-free start ([9], modified). ∆t = 0.2 s. B. Positions of body axis, right 5th leg and COM for one individual (E. s.). Origin of ordinates on caudal end of carapace at starting point. ∆t = 0.02 s (0.1 s for body axis). 2.3.2. Experiments The performed tests are adapted by the setup of the experiments carried out in section 2.1.1. Therefore, the robot starts in a supine position at the top of the tank. To initialize the turning, the two legs will strike simultaneously. The appropriate kinematic pa- rameters were determined by an iteratively change of the striking angle θ at the maximum angular strik- ing velocity ω of 4.76 s−1. The acting forces were kinematically analyzed by equation 1. 3. Results Falling E. sinensis show a body rotation along the transversal axis during free fall (see Figure 2A). The rotation is initiated directly after release by one stroke cycle of the 5th leg pair followed by a passive phase and stopped again with the 5th leg pair by a counterstrike. The turn is completed in 0.76 ± 0.15 s which leads to a mean angular velocity of ω = 4.3 s−1. In Figure 2B the trace of one individual turning maneuver is given in sagittal view. Besides the reorientation of the body axis the trajectory of the 5th leg tip is given to show the efficient strike and counterstrike. The moving location of the COM testifies the complex interaction of the acting forces. 3.1. Force calculation of 5th leg pair by kinematic approach and PIV Forces which are induced by moving the hind leg pair in phase 2a (see Figure 3B) have a maximum at 33% of the stroke with a mean force of 0.018 ± 0.004 N and a maximum peak of 0.023 ± 0.002 N. The total momentum during phase 2a is 0.0036 Ns. For phase 3a the maximum is at 35% of leg stroke with a mean force 0.015 ± 0.003 N and a maximum peak of 0.018 ± 0.002 N. Figure 3. Results from experiment and multibody- simulation for A: body angle and B: induced forces. Experimental forces are mean values of the force calcu- lations. Body angle of experimental data (blue rhomb), n=8 and multibody-simulation (red line) are showing the same trend in time and angle. The calculated forces are 1.5 to 2 times higher (PIV: 2a n=9, 3a n=3, Kinematic: 2a n=28, 3a n=11) than the forces used in the simulation. The PIV data shows that propulsive forces were produced by the hind legs to initiate and stop the rotation. Peak forces to initiate the rotation (45% of stroke cycle) of 0.021 N (± 0.001 N) (n=9) lead to a total momentum of 0.0021 Ns. 3.2. Simulation The simulation shows that a constant force of 0.0045 N acting on each hind leg respectively for 0.2 s leads to a successful rotation within time. This means that the total momentum of 0.0018 Ns causes a constant angular velocity of 3.99 s−1 in the middle part of the turning maneuver, part 2b in Figure 3. 3.3. Technical Model The experiments confirmed the assumption, that a stroke with two legs is sufficient for a 180◦-turn (see Figure 4B). The correspondent kinematic parameters are θ = 2.44 and ω = 4.76 s−1. The complete turn was herewith possible within 1.00 ± 0.12 s (n = 6). An exemplary run is shown in Figure 4B. The turn is divisible into the same five phases like E. sinensis. During the initializing (1) the chassis is descending in a supine position( till I/ 2). The striking phase (2a) includes the active leg movement and thereby the ac- celeration (I-II/ 2-3). The following gliding phase (2b) describes the passive turning based on inertia (II-III/ 55 T. Riphaus, F. Hoffmann, S. Labisch Acta Polytechnica CTU Proceedings Figure 4. Turning movement of O-Typus. B. observed turning behavior with five characteristic Phases (Initializing PI=̂1, Striking PS=̂2a, Gliding PG=̂2b, Counterstriking PC=̂3a, Descending PD=̂3b); ø∆t = 0,36 s. A. Course of the body axis and trajec- tory of the left leg’s tip. Leg positions at the end of PI (I), end of PS (II) and the target position (III) in grey. IV shows the maximal overwinding; ∆t = 0,2 s (0,1 s during passive turning (II-III)). C. Angle of attack (green) and angular velocity (red) of three runs. 3-5). When the horizontal target position is reached, the counterstriking phase (3a) starts with a counter movement of the legs, to stop the rotation (III-IV/ 5-6). Thereafter the chassis descends (3b) in a stable up-right position. Additionally Figure 4A shows the trajectory of one striking leg. Point I poses the initiation of the stroke, Point II the end of the stroke, Point III the end of the 180◦-turn and thereby the initiation of the counter- strike and Point IV the eventual slight overwinding. The graphs of Figure 4C are describing the shifting of the angle of attack and its angular velocity. The same phases can be located here. The stroke starts at 2000 ms and leads over from phase 1 to 2a. The maximal angular velocity of 4.06 s−1 is reached at the end of 2a. Subsequently ω starts decreasing within phase 2b. After approximately 1150 ms the turn is completed and results in phase 3a, which prevents overwinding. Afterwards the body descends stable with minor wobbling. The force which is induced by the moving leg pair in phase 2a has a mean value of 0.0175 N. This propulsive force is acting over 0.51 s, which results in a total generated momentum of 0.00898 Ns. Figure 5. Comparison of characteristic turning phases of (A) E. sinensis, (B) O-Typus and (C) the simulated model. Phases: 1 (yellow), 2a (blue), 2b (green), 3a (orange) and 3b (grey). 4. Conclusions 4.1. Comparison of E. sinensis with its Models A qualitative comparison of the righting behavior of E. sinensis with the turning characteristics of the adapted robot and the multibody simulation is shown in Figure 5. O-Typus shows like the biological model a complete rotation after a single stroke. Its rotation could be divided into the exact same phases like E. sinensis. The correspondence can be seen by comparing the phases. The duration of some phases are partially varying. A considerable difference between A and B can be seen within phase 2b. According to this, the larger O-Typus needs significantly more time to complete the passive turn. A full comparison of the designed and measured parameters is shown in Table 1. While the numerical model could be simulated in the original scale, the robot had to be enlarged to fit all devices. However the mean density was adapted. The measured propulsive forces of E. sinensis deviate among the two methods. While the PIV data corresponds closely to the values gained by the simulation, the kinematic approach matches nicely with the acting forces on O-Typus. Due to the considerable bigger dimensions of O-Typus, the total momentum to generate a similar turn is about 3 times of E. sinensis value. Herewith O-Typus reaches a mean angular velocity of 3.93 ± 0.30 s−1, which deviates no more than 8.6% of its model. 4.2. Transferability to Aquatic Robots As the motivation for this study showed, there is a desperate need in new stabilizing systems for technical aquatic applications. The observed and tested righting behavior of E. sinensis revealed a high efficiency and therefor a potential approach to solve the problem. In order to verify, whether this new stabilizing system is an optimization compared to current methods, further 56 vol. 7/2017 Righting of Chinese Mitten Crabs (Eriocheir sinensis) and their Models Parameter E. sinensis Simulation Robot Kinematic PIV Scaling 100% 100% 200% Mass [g] 88.65 87 220 Mean Density [g/cm3] 1.12 1.12 1.15 Mean Force [N] 0.018 0.011 0.009 0.0175 ± 0.004 ± 0.003 Momentum [Ns] 0.0036 0.0021 0.0018 0.00898 Mean Angular vel. [s−1] 4.30 3.99 3.93 ± 0.30 Table 1. Summary of relevant parameters regarding the turning of E. sinensis, its numerical model and its technical application. Values of the measured force are mean values during the striking phase 2a. The angular velocity is the averaged value of the whole turn. Figure 6. Smart-E. Modern thruster-driven AUV, which is declared to be highly efficient [10]. studies have to be carried out. Thus, two different efficiencies have to be compared. On one hand, the power consumption has to be evaluated. At the current state of research no state- ment is possible, as there is no data about robotic devices, driven by thruster, available to compare it with. It would be worth to use an existing modern Robot like the AUV Smart-E [10], to compare the two mechanisms at the same chassis. On the other hand, time efficiency has to be consid- ered. Most remotely operated vehicles are featuring slow turning capabilities. Unfortunately non pub- lished the specific turning kinematics. However, the maximum angular velocity of 4.25 s−1 reached by O- Typus is considerably (visually) higher than of most found thruster-driven devices. Acknowledgements Special thanks are referring to the Department of Biomimetics in Bremen, which made these studies possible. References [1] G. Krummel, K. N. Kaipa, S. K. Gupta. A horseshoe crab inspired surf zone robot with righting capabilities. In Proc. of the ASME 2014 Int. Design Eng. Tech. Conf. & Comput. and Inform. in Eng. Conf. 2014. doi:10.1115/DETC2014-34679. [2] E. Chang-Siu, T. Libby, M. Tomizuka, R. J. Full. A Lizard-inspired active Tail enables rapid Maneuvers and dynamic Stabilization in a terrestrial Robot. IEEE/RSJ Int Conf of Intelligent Robots and Syst 2011. doi:10.1109/IROS.2011.6094658. [3] P. Bausch, F. Hoffmann, A. B. Kesel. Im freien Fall unter Wasser - Lagestabilisierung bei Wollhandkrabben. Eine experimentelle Analyse. In Bionik: Patente aus der Natur (ed. A.B. Kesel und D. Zehren). 155-161, Bionik-Innovations-Centrum, Bremen, 2015. [4] F. Hoffmann, S. Labisch, R. Sonntag, et al. Analysis of underwater Righting Reflex of the Chinese mitten crab (Eriocheir sinensis) using PIV and Numerical Multibody Simulation. In Proceedings of the Annual Main Meeting of the Society for Experimental Biology, Prag, [A 11.68]. 2015. [5] T. Riphaus, S. Labisch. Konstruktion des biologisch inspirierten Roboters O-Typus zur Implementierung und Untersuchung des Stabilisationsverhaltens von Decapoden. In Proceedings of the 7. IEEE Germany Student Conference, Magdeburg. 2016. [6] W. Thielicke, E. J. Stamhuis. PIVlab - towards User- friendly, affordable and accurate Digital Particle Image Velocimetry in MATLAB., 2014. J.OpenRes.Softw.,2,. doi:http://doi.org/10.5334/jors.bl. [7] J. M. V. Rayner. A new Approach to Animal Flight Mechanics. J Exp Biol 80:17–54, 1979. [8] G. R. Spedding, J. M. V. Rayner, C. J. Pennycuick. Momentum and Energy in the Wake of a Pigeon (Columba livia) in Slow Flight. J Exp Biol 111:81–102, 1984. [9] P. Bausch, P. Berger, F. Hoffmann, A. B. Kesel. How to land on your Legs - Underwater Righting Reflexes in Crabs and Crayfish. In Proceedings of the Annual Main Meeting of the Society for Experimental Biology, Prag, [A 11.67]. 2015. [10] B. Meyer, K. Ehler, C. Osterloh, E. Maehle. An autonomous omnidirectional Underwater Robot. Journal of Behavioral Robotics pp. 204–210, 2013. 57 http://dx.doi.org/10.1115/DETC2014-34679 http://dx.doi.org/10.1109/IROS.2011.6094658 J. Open Res. Softw., 2,. http://dx.doi.org/http://doi.org/10.5334/jors.bl Acta Polytechnica CTU Proceedings 7:53–57, 2017 1 Introduction 2 Material and Methods 2.1 Eriocheir sinensis - biological model 2.1.1 Kinematic analysis 2.1.2 Piv measurements and force calculation 2.2 Simulated Model 2.3 O-Typus - Technical Model 2.3.1 Setup 2.3.2 Experiments 3 Results 3.1 Force calculation of 5th leg pair by kinematic approach and PIV 3.2 Simulation 3.3 Technical Model 4 Conclusions 4.1 Comparison of E. sinensis with its Models 4.2 Transferability to Aquatic Robots Acknowledgements References