Abstract

Continuum manipulators are a class of robots with many degrees-of-freedom, leading to highly flexible motion with inherent compliance. These attributes make them well suited for manipulation tasks and physical interaction with the environment. A high impact yet challenging field for exploring continuum robot designs is free-floating underwater manipulation with a remotely operated vehicle (ROV). In this article, we propose a modular, reconfigurable, cable-driven continuum arm for free-floating underwater manipulation and present a corresponding kinematics, control, and computation framework. The mechanical design consists of a continuum arm, an actuation unit, and a waterproof enclosure. The kinematics model is introduced as two mappings between three spaces: the joint space, the configuration space, and the task space. The differential kinematics for each mapping is also derived. An electronics system design is proposed for underwater applications, including the communication framework between the topside computer (above surface), on-board computer, and manipulator mechatronics. Experimental validation is presented to demonstrate the robot’s underwater functionality, test the limits of its articulation, and evaluate the arm’s stiffness. Future work includes field testing with an ROV platform and development of advanced controls and planning for manipulation tasks.

1 Introduction

Continuum robots are an unconventional class of robots inspired by biological features such as an elephant’s trunk [13] or an octopus’s tentacles [47]. They have a very large number of degrees-of-freedom (DOFs), allowing them to achieve a nearly continuous bending shape. One unique feature of continuum robots is their inherent flexibility and compliance, making them ideal for safe physical interactions with the environment. The entire robot arm can be used as a manipulator to grasp objects in a whole arm caging style, or a gripper can be attached to the end of the arm to complete more versatile and precise tasks.

A significant challenge for continuum robots is the difficulty in modeling and control due to the high flexibility. Many researchers have investigated this problem and proposed solutions, including work by Chirikjian and Burdick to analyze the kinematics of hyper-redundant manipulators [811]. Gravagne and Walker proposed a kinematic model and dynamic model for planar motion of a continuum robot and analyzed the forces, compliance, and control [12,13].

Continuum robots have attracted interest primarily for biomedical applications because of their flexibility and compliance enabling safe interaction with the anatomy [1417]. These qualities also make them suitable for mobile robot manipulation. McMahan et al. demonstrated the field capabilities and manipulation of a pneumatic continuum manipulator attached to a mobile robot platform in Ref. [18]. Rone et al. used a continuum tail to stabilize a walking mobile robot [19,20]. However, there are no continuum manipulators suitable for underwater free-floating manipulation yet. Continuum robots for underwater applications have been limited to seabed locomotion and seabed mobile robots. For instance, Liu et al. used a continuum arm and gripper for underwater manipulation on an underwater walker [21]. Gong et al. implemented a soft manipulator for seafloor manipulation [22]. Cianchetti et al. developed an octopus-inspired robot that utilizes eight continuum arms for locomotion and manipulation [7].

Free-floating manipulation is a unique problem that utilizes the mobility of the ROV platform to achieve a much larger reachable workspace while maintaining fine manipulation capability, which opens up many new possible applications. Unfortunately, this also introduces many challenges for current robot designs, which has been addressed in literature [2325]. Often these designs use rigid-link serial arm structures that are heavy, expensive, and difficult to make waterproof, which limits their applications and increases the barrier for small business and hobbyists to adopt.

In this article, a modular open-source cable-driven continuum robot design is proposed. This robot is designed for free-floating underwater manipulation in conjunction with an underwater remotely operated vehicle (ROV). The continuum arm is made using a series of eight links connected via two DOFs universal joints. The robot utilizes Nitinol cables to actuate the arm and lead screw actuation to provide increased mechanical advantage. A waterproof enclosure and sleeve surrounds the robot arm and its actuation unit to allow the robot to be operated underwater. The enclosure is easy to attach and reconfigure to an ROV, allowing for custom configurations for specific tasks. The robot is fabricated using off-the-shelf components and 3D printed parts, which allows the design to be easily customized and reproduced at low cost. Important design decisions are discussed to support the robot’s utility for underwater manipulation and compatibility with an ROV.

To control the robot, a two-part kinematics model is adapted from Refs. [26,27] to accommodate our discrete universal joint design. The two-part model represents mappings among three kinematics spaces: the task space, the configuration space, and the joint space of the robot. These spaces represent the end link position, the arm shape, and the motor positions. The kinematics model maps the joint space to the configuration space and the configuration space to the task space. To derive the kinematic model, the arm is modeled as a series of rigid links connected by universal joints. While the arm design is technically discrete, not continuous, its high flexibility and degrees-of-freedom allow us to model and approximate it as a continuum robot and will be referred to throughout this paper as a continuum arm. The differential kinematics model is also derived.

There have been open-source hardware projects in the robotics community that have helped significantly lower the barrier of entry to hardware prototypes, such as the OpenHand project at Yale2 [28,29], the Soft Robotics Toolkit developed at Harvard and Trinity College Dublin3, and the FAB@HOME open-source 3D printer project started at Cornell.4 We are committed to contributing open-sourced mechanical, electronics, and software designs to the research community that enables customizability tailored to different users’ needs.

Main contributions of this work include four parts. First, the self-contained, sealed actuation unit allows the arm to be easily mounted in a number of configurations on a mobile robot platform, as demonstrated in Fig. 1(a). Second, this also allows the robot to be used for free-floating manipulation. To the best of our knowledge, we are the first to design and integrate such modular continuum robots for this application. However, the control and implementation of the entire free-floating system is outside the scope of this article and is planned for future works. Third, although the mechanism design utilizes discrete universal joints instead of continuum structures, we developed a hybrid kind of kinematics model that accommodates the continuum shape assumption and the discrete joint implementation, which could be useful for other similar designs in the future. Finally, the proposed design demonstrates a commitment to open-source usability. Components were chosen to be affordable and readily available. The CAD models are available for free download through our website, which can be used to 3D print a large portion of the components. The electronics design allows the entire robot to be viewed as modular USB devices, which is advantageous for easy integration with the electronics systems for a variety of on-board computer platforms. The addition of inertial measurement units (IMU)’s on each link allows for real-time estimation of the continuum arm’s position for ease of teleoperation.

Fig. 1
Conceptual CAD models of possible system configurations for integrating modular manipulators. Base BlueROV2 CAD Model from Blue Robotics, Inc. (a) Configuration 1: standard manipulation, (b) configuration 2: seafloor manipulation, and (c) configuration 3: crawling gait.
Fig. 1
Conceptual CAD models of possible system configurations for integrating modular manipulators. Base BlueROV2 CAD Model from Blue Robotics, Inc. (a) Configuration 1: standard manipulation, (b) configuration 2: seafloor manipulation, and (c) configuration 3: crawling gait.
Close modal

This article is organized as follows: in Sec. 2, the continuum robot design design, materials, and fabrication is detailed. In Sec. 3, the kinematics and differential kinematics models are derived. Section 4 discusses the control architecture and electronics hardware. Section 5 discusses tests to demonstrate the functionality of the proposed design. Finally, Sec. 6 outlines next steps in the design refinement and provides concluding remarks.

2 Mechanism Design

The continuum robot contains two major modules: the flexible continuum arm, which is driven by four Nitinol cables, and the actuation unit, which pulls the cables of the arm and houses the electronics. Both modules are enclosed by a waterproof case, which is composed of a cylindrical enclosure and a flexible bellows.

There are multiple methods of actuating continuum robots, including pneumatic actuation [18,21], shape memory alloy actuators [4,5], cable or tendon actuation [1,14,30,31], mechanical linkages [2,19,20,32], or some combination of multiple methods. A cable-driven mechanism is a good choice for an underwater continuum robot because it can be made compact and centralized. Combining with high mechanical advantage actuation such as lead screw, it can provide large amounts of force to resist strong external disturbances such as current, and its high back-drive impedance enables low power consumption to sustain stationary posture or configuration.

The continuum robot is designed to be attached to a free-floating ROV platform, such as the BlueROV2 by Blue Robotics, Inc.5 The modularity of the robot makes it possible to reconfigure the system mechanically for different tasks. Figure 1 shows example configurations: Fig. 1(a) shows two arms implemented in a forward configuration to complete fine manipulation tasks (e.g., object handling); Fig. 1(b) shows two arms implemented in a downward configuration to pick and place objects on the seafloor; and Fig. 1(c) shows four arms implemented in a forward configuration to crawl along a fence.6 The continuum robots can be reconfigured on other ROV platforms as well. In the future work, algorithms could be investigated to optimize the reconfiguration (e.g., number of arms, mounting positions and orientations) based on task specifications.

2.1 Continuum Arm Design.

The continuum arm uses a series of rigid links connected by universal joints that each allow rotations in two axes. Extension springs connect each link to provide structure stiffness of the arm. The robot could be made longer or shorter to accommodate a specific task by adding or subtracting links to the arm. The choice of design parameters of each link has a large effect on the kinematics calculations (e.g., the maximum required stroke length of the robot) as well as mechanics calculations (e.g., deformation under disturbance force). Table 1 summarizes the values of parameters used for the continuum robot. The parameter r is the radial distance from the center of the link to the cable routing point. The length L is the distance from the center of the link to the center of the universal joint.

Table 1

Continuum arm design parameters

ParameterSymbolValue
Radiusr21 mm
LengthL18 mm
Number of linksNl8
Number of cablesNc4
Maximum bending angleθmaxπ/2
Spring stiffnessks528.9 N/m
Unstretched spring lengthls25.4 mm
Diameter of Nitinol backbonesdb1.5 mm
ParameterSymbolValue
Radiusr21 mm
LengthL18 mm
Number of linksNl8
Number of cablesNc4
Maximum bending angleθmaxπ/2
Spring stiffnessks528.9 N/m
Unstretched spring lengthls25.4 mm
Diameter of Nitinol backbonesdb1.5 mm

Four Nitinol cables are attached to the final link via set screws. While the minimum number of pulling cables to actuate a continuum robot is two, we use four cables, and the actuation redundancy provides a larger wrench-feasible workspace in configuration space and generally reduces the peak tendon pulling load per individual cable. Large diameter cables were chosen to further increase the stiffness of the arm.

The continuum arm consists of eight 3D printed links that are connected via aluminum universal joints. These links have four routing holes for the Nitinol cable. Four stainless steel extension springs are attached between each link. The final link includes four threaded holes to house a set screw to fix each cable. Each link also has a flat area to accommodate an IMU for real-time shape sensing. This process is described in more detail in Sec. 5.

2.2 Actuation Unit Design.

The actuation unit uses Dynamixel XL430-W250-T servo motors produced by Robotis, Inc.7 to actuate four lead screws, providing the linear motion required to pull the cables. The lead screw mechanism is larger and heavier than alternative cable driving mechanisms such as pulleys, but provides higher mechanical advantage and therefore increases the maximum load of the robot. The added system weight is not a major concern since underwater vehicles generally have a high payload capacity due to the effects of buoyancy (Fig. 2).

Fig. 2
CAD models of the continuum robot arm design, the IMU sensor integration, and the actuation unit: (a) continuum arm and flexible sleeve, (b) IMU sensor integration, and (c) actuation unit exploded view
Fig. 2
CAD models of the continuum robot arm design, the IMU sensor integration, and the actuation unit: (a) continuum arm and flexible sleeve, (b) IMU sensor integration, and (c) actuation unit exploded view
Close modal
An important step in the design of the actuation unit was sizing the motors and lead screw pitch to withstand a load of up to 100 N per cable, which represents the force the robot might experience during a large current or other disturbance. The servo motors can provide up to 1.50 N·m of torque, so the lead screw travel distance must be chosen to achieve the desired maximum load. The relationship between pitch p, thread starts n, travel distance per revolution Dt (i.e., ratio of linear slide velocity to motor angular velocity), stall torque τs, and cable tension T is given by the equations:
(1)
(2)
(3)
Equation (2) yields a maximum travel distance of 15 mm/rev. The user can adjust the strength of the robot by simply changing the lead screws with a different thread; a smaller travel distance will allow the robot to resist larger loads, while a larger travel distance will allow it to move faster. For the proposed design, a 1/4 in.-12 ACME lead screw with four thread starts was chosen as a low-cost option with a large travel distance (8.382 mm/rev) while still withstanding the required loads. The lead screw stroke length Qmax is a function of the maximum desired bending angle given by the following equation, which is derived from the kinematic calculation of the motor actuation distance qi:
(4)

The actuation unit consists of three parallel 3D printed mounting plates (referred to as top, bottom, and electronics plates), which mount the continuum arm, motors, and electronics boards, respectively. The lead screws are attached to the motors using 3D printed couplings with a shoulder screw pinned through the lead screw to secure it and supported by bearings on the top plate. The top and bottom plates are aligned and fastened together using two stainless steel standoffs and two 3D printed supports, which accommodate the linear potentiometers. A 3D printed coupling connects each lead screw to a guide rail. To accommodate linear position feedback via a contact-based membrane potentiometer (ThinPot), the coupling also has a spring plunger attached to ensure firm and continuous contact with the potentiometer.

2.3 Waterproof Enclosure Design.

The actuation unit is housed in a cylindrical waterproof enclosure minimally modified from the Blue Robotics, Inc. 6 in. series watertight enclosure. It consists of an acrylic tube and a collar and flat plate at each end. O-rings are used to seal the connections between these components, and screws are used to secure the end plate to the collar. Sealed penetrators attached to the back plate allow communication cables to enter the enclosure without compromising the waterproof seal. The diameter of the enclosure was chosen to accommodate the motor layout, and the length was chosen to accommodate the required stroke length Qmax, calculated from the continuum arm design parameters as discussed later in this article.

The continuum arm is sealed via a flexible rubber bellows. Although originally intended for linear motion, the bellows also work well for the robot design by acting as a lightweight and extremely flexible sleeve. The diameter and the length of the bellows are chosen primarily to accommodate the continuum arm dimensions. Since the bellows are flexible, their length can stretch or compress to match the length and motion of the arm. The diameter also changes inversely proportional to the changes in the length; however, these changes are very small in comparison and can be neglected. The open ends of the bellows are closed via the large orange silicone collar and large plastic disk. Waterproof epoxy was put into all of the gaps in the joints to seal the bellows. The epoxy was also put between the flat collar and the end plate of the actuation unit to attach and seal the bellows to the actuation unit enclosure. The end plate and bellows are detachable from the collar together to access the inner mechanisms by unscrewing the fasteners, providing easy access for maintenance.

3 Mathematical Modeling

A continuum robot model is applied to the arm, which is actually a series of eight rigid links connected via two degrees-of-freedom universal joints. The configuration space variables are used to calculate the rotation matrices and displacements of each link relative to the robot’s base coordinate system. The sum of the distances of the cable routing points is used to calculate the motor position for the joint space, and the location of the end link is used to find the task space variables (the position and orientation of the end disk). The differential kinematics and Jacobian for the configuration space to joint space mapping and configuration space to task space mapping are also derived.

This section is broken up into two subsections. Section 3.1 details the kinematic modeling of the continuum arm, which defines the relationships between the robot shape, position and orientation of the end disk in the arm, and the linear position of the motor. Section 3.2 includes calculations of the Jacobian matrices for relating the velocities of these variables.

3.1 Kinematics.

The kinematic modeling of our design is facilitated via three sets of variables that are referred to as the task space variables, the configuration space variables, and the joint space variables. The configuration space variables indicate the in-plane bending angle θ and the bending plane direction angle ϕ of the manipulator as a continuum structure. The joint space variables refer to the actuation displacements for pulling cables, represented by q1 to qNc. The task space variables are defined as the position and orientation of the end disk represented by a vector xe and a rotation matrix Re, respectively. This leads to two kinematic mappings: the configuration to joint space mapping, and the configuration to task space mapping. Key assumptions for the kinematic model are as follows: (1) the continuum arm bends in a plane, (2) the bending angles between each two adjacent disks are the same, (3) the Nitinol cables are inextensible.

3.1.1 Configuration to Joint Space Mapping.

The configuration space ψ contains the in-plane bending angle θ and bending direction angle ϕ as illustrated in Fig. 3.
(5)
Fig. 3
Continuum robot mathematical model: (a) continuum robot configuration space, and (b) discrete rigid-link diagram
Fig. 3
Continuum robot mathematical model: (a) continuum robot configuration space, and (b) discrete rigid-link diagram
Close modal
The continuum arm is modeled as a series of transformations in between each disk consisting of an axial translation, a rotation about the two axes of the universal joint, and another axial translation. A unique coordinate frame is created for each disk with the origin located at the centroid of each disk:
(6)
where vector [a] denotes a skew-symmetric matrix formed by a according the convention:
(7)
where the matrix exponential expression (e.g., eϕ[z^]) represents a rotation of the given angle about the given axis.
The centroid of the ith link is given by:
(8)
The location of the ith joint (i.e., the center of the universal joint that connects two adjacent links) is given by:
(9)
The location of each routing point, denoted by pi,j, is calculated by a translation of r in the disk coordinate frame positive and negative X and Y axes:
(10)
(11)
The location is expressed in the previous disk coordinate for the convenience of calculating the distance of the jth tendon between the i and (i − 1)th disks, di,j:
(12)
The sum of the distances between consecutive routing points is calculated for each cable and gives the cable length lj for the given configuration:
(13)
The linear motor position qj can be calculated by finding the difference between this length and the neutral length l0,j, which is the length of the cables when the arm is straight. The linear motor position represents the position of the attachment point between the lead screw and the Nitinol actuation cables and is given by
(14)
(15)

3.1.2 Configuration to Task Space Mapping.

The 6 DOF task space model is given by the position and orientation of the final link, given by xe and Re, respectively. As mentioned earlier, this model assumes that the total bending angle θ is evenly distributed between each link. The pose is calculated iteratively using Eqs. (6), (8), and (9):
(16)
(17)

3.2 Differential Kinematics

3.2.1 Configuration to Joint Space Mapping.

The differential kinematics of the configuration to joint space is defined as follows:
(18)
where
(19)
The Jacobian is found by using the chain rule to find the derivatives of the norm pi,jpi1,j:
(20)
where
(21)
The derivative of pi,j is given by
(22)
(23)
(24)
where we denote ψ1 = θ, ψ2 = ϕ for indexing convenience and Γj was defined in Eq. (11).
The derivatives in Eqs. (22)(24) are ultimately dependent on the derivatives of the rotation matrices Ri, which are given by
(25)
(26)
where the matrices A and B are derived as follows:
(27)
(28)

It can be seen from Eqs. (25) and (26) that the derivatives of the rotation matrices are solved iteratively, similar to the kinematics presented earlier.

3.2.2 Configuration to Task Space Mapping.

The differential kinematics of the configuration to task space includes the linear velocity part and the angular velocity part. The linear velocity part is defined as follows:
(29)
(30)
The Jacobian is found by differentiating Eq. (16):
(31)
where the intermediate terms C, D, and E are defined as follows:
(32)
(33)
(34)
The angular velocity part is defined by:
(35)
(36)
where ω represents the angular velocity of the end disk expressed in the base frame.
The key to finding Jωψ is to utilize the full derivative formulation of a rotation matrix as explained in Refs. [33,34]:
(37)
where, as introduced in Eq. (7), [()] is the skew-symmetric matrix operator.
In addition to Eq. (37), in which we expressed the full derivative of the end disk orientation Re with the angular velocity, we can also derive the full derivative with the configuration space velocities:
(38)
By vectorizing Eqs. (37) and (38), respectively, and by equating the right-hand sides of them, we are able to find the expression for the Jacobian Jωψ. This nontrivial derivation is detailed in Appendix  A.
(39)
(40)
where the definitions and computations for (·)+, Vec(·), n^x, n^y, n^z, Reθ, and Reϕ are described in details in Appendix  A.

4 Electronics, Systems, and Controls

The kinematic model derived in Sec. 3 is implemented in robot operating system to achieve robot control. The communication is facilitated using USB between the topside control computer and the on-board electronics boards. The electronics hardware is prototyped using two microprocessor boards that communicate with the central computer and the auxiliary sensors. As shown in Fig. 4, the first microprocessor board is dedicated to motor control, while the second is dedicated to sensor readings and communications. The motor control board is connected to four Dynamixel XL-430-W250-T servo motors produced by Robotis Inc. (37). The sensor board is connected to four linear potentiometers that are used as additional position sensors. Inertia measurement units (IMUs) are attached to each of the disks of the continuum arm and are also connected to the sensor board.

Fig. 4
ROV integration schematic
Fig. 4
ROV integration schematic
Close modal

A user control interface was developed using a gamepad for the purpose of teleoperation. The command is captured using the joy stick, in the format of velocities in θ˙ and δ˙. The velocity command is sent to the kinematic model and then used to issue low-level commands for motor positions. As shown in Fig. 4, a tether is used to communicate with the ROV, and the robot electronics boards connect to the on-board central computer via USB. For testing purposes, these boards were connected to the topside computer directly. The ROV hardware and software schematic is outside the scope of this article.

5 Experimental Validation

Several tests were conducted to demonstrate the functionality of the robot, articulation, and stiffness. This section outlines the experiments performed and their outcomes.

5.1 Motion, Articulation, and Visualization.

Basic motion control of the robot is fundamental to its success. This was tested through joystick operation of the robot, utilizing the control framework discussed in Sec. 4. To minimize the effects of gravity, the robot was placed in a vertical orientation. A bending angle of up to 90 deg was demonstrated in a full 360 deg rotation. The arm can go from 0 deg to 90 deg in approximately 6 sec and complete a revolution in approximately 30 sec with the bending angle held at 90 deg while using a joystick controller. A fifth-order polynomial trajectory planner was also implemented to allow the robot to smoothly move from one configuration to another.

The functionality of the IMUs was also demonstrated. The orientation readings from each IMU is used to calculate the relative rotation between each link. These are then converted into unique joint angles, representing the bending angle along each of the two DOFs of the universal joints. The combination of these joint angles allows us to visualize the arm in real time using RVIZ and a URDF file created from the CAD model. An example of this visualization can be seen in Fig. 5. This figure demonstrates the accuracy of the visualization in representing the actual position of the arm over different angles.

Fig. 5
IMU visualization comparison
Fig. 5
IMU visualization comparison
Close modal

5.2 Underwater Functionality.

A fundamental goal of the robot is to be waterproof. The first step to determine this is to create and maintain a vacuum inside the enclosure. This allows us to test the seal without risking exposure or damage to the electrical components. A vacuum of approximately −400 mmHg was drawn and the pressure was monitored over time to ensure it did not drop, which would indicate a leak. The enclosure successfully maintained a vacuum for over fifteen minutes.

After the initial vacuum test, the robot was fully submerged in a tank to provide a more realistic test. The robot did not indicate any leaks when submerged; moreover, external power and communications were fully functional. The same basic joystick control test was performed underwater successfully, and a complete circular motion was achieved underwater. This setup is shown in Fig. 6. This figure shows the accuracy of the arm in achieving both a perfectly straight configuration and a maximum bending angle configuration while submerged.

Fig. 6
Submerged testing of the robot prototype: (a) straight configuration (θ = 0) and (b) bending configuration (θ = π/2)
Fig. 6
Submerged testing of the robot prototype: (a) straight configuration (θ = 0) and (b) bending configuration (θ = π/2)
Close modal

Preliminary water tests also highlighted some design shortcomings that could be overcome. One is that both the arm and actuation units are very buoyant due to the large empty volume inside the bellows and enclosure. To resolve this, small weights were added around the bellows as evenly distributed as possible. These were attached by stringing the weights on a wire and attaching a snap swivel at both ends of the wire. The snap swivels were then connected together to form a loop around the bellows with a small enough diameter to stay in place. Multiple loops were attached one at a time until the arm was able to remain straight when submerged, as shown in Fig. 6(a).

5.3 Arm Stiffness Test.

The success of the robot in resisting external disturbances, such as a force caused by a strong current, is also crucial for demonstrating the robot’s suitability for ROV manipulation. An experimental setup was designed to apply a known weight perpendicular to the end of the arm via a string. The setup contains a rail with a pulley that can be adjusted to change the angle and position of the string to keep the applied force perpendicular to the arm. The setup can be modified to flip the direction of the string to apply an outward radial force (against the direction of bending) or an inward radial force (with the direction of bending), as shown in Fig. 7. This experiment was repeated with increasing weight from 0 to 1.5 kg with arm bending angles ranging from 30 deg to 90 deg. For consistency, ϕ was kept at 0 deg. A tracking square was attached to the end link to track the deflection of the arm as the weight was applied by taking pictures with a mounted camera.

Fig. 7
Stiffness testing experimental setup: (a) inward bending experiment example and (b) arm stiffness test setup
Fig. 7
Stiffness testing experimental setup: (a) inward bending experiment example and (b) arm stiffness test setup
Close modal

The stiffness of the arm was analyzed by calculating the position and orientation of the tracking square over a series of frames for a specific bending angle. The algorithm detects Harris corners within a user-specified range, and the correct points corresponding to the tracking square are manually selected. The length of the square is known and used as a scale to calculate the exact displacements. Finally, the weight is plotted as a function of the calculated displacements and a linear polynomial curve is fit to the data to numerically find the stiffness, as shown in Fig. 8. The stiffness for each bending angle is summarized in Table 2.

Fig. 8
Experimental results of inward and outward radial loading
Fig. 8
Experimental results of inward and outward radial loading
Close modal
Table 2

Arm stiffness at various bending angles

Angle (deg)Spring stiffness inward (N/m)Spring stiffness outward (N/m)
3060.5720112.0680
4063.5241100.0550
5066.7533102.6783
6063.4751124.0438
7066.4339126.0441
8071.078689.0820
9071.679394.1832
Angle (deg)Spring stiffness inward (N/m)Spring stiffness outward (N/m)
3060.5720112.0680
4063.5241100.0550
5066.7533102.6783
6063.4751124.0438
7066.4339126.0441
8071.078689.0820
9071.679394.1832

The experimental data show that the stiffness of the arm is fairly linear for every angle and loading direction. In general, the stiffness increases as the bending angle increases. The outward loading case has a noticeably larger stiffness when compared to the inward loading case.

6 Conclusion

This article introduces the design of a modular, open-source continuum robot for underwater manipulation. The design is easy to reproduce and compatible with a variety of ROV platforms. This design consists of a tendon-driven continuum arm, lead screw actuation mechanism, and a waterproof enclosure to seal the entire robot. We derived a two-part kinematic and differential kinematic model mapping the joint space, configuration space, and task space model. We also developed the necessary electronics, controls, and communications. The current design and future modifications and customization will be released online.8

The operation of a functional prototype has been successfully demonstrated in testing preliminary functions, as shown in Fig. 6. The control in configuration space was achieved through user joystick operation, as discussed in Sec. 4. Underwater operation was tested by fully submerging the robot in a tank. The robot did not indicate any leaks when submerged; moreover, external power and communications were fully functional, allowing the robot to achieve the same range of motion underwater. Neutral buoyancy of the arm was achieved by incrementally adding ballast weights around the bellows until the arm was able to maintain its configuration without deflection.

While this testing has shown the robot to be fundamentally operational, there is still additional development and refinements to improve the robot’s performance and integrate it with a floating ROV platform. A rigorous analysis of the design space will be carried out to analyze how the robot can be optimized for metrics such as maximum stiffness, maximum supported external load, and others. The robot must still be fully integrated with the mobile platform, which includes creating a mounting design and redesigning the electronics system to facilitate power transmission and communications. Finally, advanced controls and planning must be implemented to accomplish specific tasks such as perching (in which the robot stabilizes itself using a pier or other underwater object) or object manipulation. Implementing these improvements will yield a robot that can overcome a wide array of difficult underwater manipulation problems.

Footnotes

Conflict of Interest

There are no conflicts of interest.

Data Availability Statement

The data and information that support the findings of this article are freely available.9

Acknowledgment

This research was supported in part by USDA-NIFA (Grant No. 2021-67022-35977).

Nomenclature

q =

joint space vector composed of linear motor position: q=[q1,,qNc]T

n =

thread starts

p =

thread pitch

r =

radius of continuum arm link

L =

length between link centroid and joint centroid

T =

cable tension

ji =

position vector of centroid of joint i

oi =

position vector of centroid of link i

pi,j =

position vector of routing point for cable j on link i

xe =

position vector of end link

J =

Jacobian matrix relating joint and configuration velocities

J =

Jacobian matrix relating task and configuration velocities

Re =

rotation matrix representing orientation of end link

db =

diameter of Nitinol cables

di,j =

distance between routing points for cable j between link i and i − 1

ks =

spring stiffness

lj =

length of cable j

l0,j =

neutral length of cable j

ls =

unstretched spring length

Dt =

thread travel distance

Nl =

number of links

Nc =

number of cables

Qmax =

maximum stroke length

Ri =

rotation matrix representing orientation of link i

θ =

in-plane bending angle

θmax =

maximum bending angle

τs =

motor stall torque

ϕ =

bending plane angle

ψ =

configuration space vector representing shape of continuum arm

ω =

angular velocity vector of end link

Appendix A: Derivation of Jωψ

We express the rotation matrix Re in Eq. (37) with the three unit vectors for derivation convenience:
(A1)
It allows us to rewrite Eq. (37) as follows:
(A2)
Let us use Vec( · ) to denote the vectorization of a matrix by stacking the columns of the matrix on top one another. We then can vectorize Eq. (A2) as follows:
(A3)
(A4)

Thereby, we have achieved the vectorization and rearrangement of Eq. (37), arriving at Eq. (A4). In addition, we note that matrix F has linearly independent columns and has a rank of 3, which suggests that it has a particular pseudoinverse as a left inverse.

Next, we consider the vectorization of Eq. (38):
(A5)
(A6)
where the two derivative matrices, Reθ and Reϕ, can be computed easily by recalling that Re represents the orientation of the last link (Re=RNl) and by using the iterative computation methods shown in Eqs. (25) and (26). Thereby, we have achieved the vectorization and rearrangement of Eq. (38), arriving at Eq. (A6).
Finally, by equating the right-hand sides of Eqs. (A4) and (A6), we obtain the expression of the Jacobian matrix:
(A7)
where (·)+ denotes the pseudoinverse matrix operation, which, in this case, is computed as follows:
(A8)

References

1.
Hannan
,
M.
, and
Walker
,
I.
,
2003
, “
Kinematics and the Implementation of an Elephant’s Trunk Manipulator and Other Continuum Style Robots
,”
J. Rob. Syst.
,
20
(
2
), pp.
45
63
.
2.
Liu
,
Y.
,
Ge
,
Z.
,
Yang
,
S.
,
Walker
,
I. D.
, and
Ju
,
Z.
,
2019
, “
Elephant’s Trunk Robot: An Extremely Versatile Under-Actuated Continuum Robot Driven by a Single Motor
,”
ASME J. Mech. Rob.
,
11
(
5
), p.
051008
.
3.
Yeshmukhametov
,
A.
,
Koganezawa
,
K.
, and
Yamamoto
,
Y.
,
2019
, “
A Novel Discrete Wire-Driven Continuum Robot Arm With Passive Sliding Disc: Design, Kinematics and Passive Tension Control
,”
Robotics
,
8
(
7
), p.
51
.
4.
Zheng
,
T.
,
Branson
,
D. T.
,
Kang
,
R.
,
Cianchetti
,
M.
,
Guglielmino
,
E.
,
Follador
,
M.
,
Medrano-Cerda
,
G. A.
,
Godage
,
I. S.
, and
Caldwell
,
D. G.
,
2012
, “
Dynamic Continuum Arm Model for Use With Underwater Robotic Manipulators Inspired by Octopus Vulgaris
,”
2012 IEEE International Conference on Robotics and Automation
,
Saint Paul, MN
,
May 14–18
, pp.
5289
5294
.
5.
Zheng
,
T.
,
Branson
,
D. T.
,
Guglielmino
,
E.
,
Kang
,
R.
,
Medrano Cerda
,
G. A.
,
Cianchetti
,
M.
,
Follador
,
M.
,
Godage
,
I. S.
, and
Caldwell
,
D. G.
,
2013
, “
Model Validation of an Octopus Inspired Continuum Robotic Arm for Use in Underwater Environments
,”
ASME J. Mech. Rob.
,
5
(
2
), p.
021004
.
6.
Cianchetti
,
M.
,
Follador
,
M.
,
Mazzolai
,
B.
,
Dario
,
P.
, and
Laschi
,
C.
,
2012
, “
Design and Development of a Soft Robotic Octopus Arm Exploiting Embodied Intelligence
,”
2012 IEEE International Conference on Robotics and Automation
,
Saint Paul, MN
,
May 14–18
, pp.
5271
5276
.
7.
Cianchetti
,
M.
,
Calisti
,
M.
,
Margheri
,
L.
,
Kuba
,
M.
, and
Laschi
,
C.
,
2015
, “
Bioinspired Locomotion and Grasping in Water: The Soft Eight-Arm Octopus Robot
,”
Bioinsp. Biom.
,
10
(3), p.
035003
.
8.
Chirikjian
,
G.
, and
Burdick
,
J.
,
1990
, “
An Obstacle Avoidance Algorithm for Hyper-Redundant Manipulators
,”
Proceedings IEEE International Conference on Robotics and Automation
,
Cincinnati, OH
,
May 13–18
, Vol.
1
, pp.
625
631
.
9.
Chirikjian
,
G.
,
1993
, “
General Methods for Computing Hyper-redundant Manipulator Inverse Kinematics
,”
Proceedings of 1993 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS ’93)
,
Yokohama, Japan
,
July 26–30
, Vol.
2
, pp.
1067
1073
.
10.
Chirikjian
,
G.
, and
Burdick
,
J.
,
1994
, “
A Modal Approach to Hyper-Redundant Manipulator Kinematics
,”
IEEE. Trans. Rob. Autom.
,
10
(
3
), pp.
343
354
.
11.
Chirikjian
,
G.
, and
Burdick
,
J.
,
1992
, “
Kinematically Optimal Hyper-Redundant Manipulator Configurations
,”
Proceedings 1992 IEEE International Conference on Robotics and Automation
,
Nice, France
,
May 12–14
, Vol.
1
, pp.
415
420
.
12.
Gravagne
,
I.
, and
Walker
,
I.
,
2002
, “
Manipulability, Force, and Compliance Analysis for Planar Continuum Manipulators
,”
IEEE. Trans. Rob. Autom.
,
18
(
3
), pp.
263
273
.
13.
Gravagne
,
I.
, and
Walker
,
I.
,
2002
, “
Uniform Regulation of a Multi-Section Continuum Manipulator
,”
Proceedings 2002 IEEE International Conference on Robotics and Automation (Cat. No.02CH37292)
,
Washington, DC
,
May 11–15
, Vol.
2
, pp.
1519
1524
.
14.
Simaan
,
N.
,
Taylor
,
R.
, and
Flint
,
P.
,
2004
, “
A Dexterous System for Laryngeal Surgery
,”
IEEE International Conference on Robotics and Automation, 2004
,
New Orleans, LA
,
Apr. 26–May 1
, Vol.
1
, pp.
351
357
.
15.
Conrad
,
B.
, and
Zinn
,
M.
,
2015
, “
Closed Loop Task Space Control of an Interleaved Continuum-Rigid Manipulator
,”
2015 IEEE International Conference on Robotics and Automation (ICRA)
,
Seattle, WA
,
May 26–30
, pp.
1743
1750
.
16.
Conrad
,
B. L.
, and
Zinn
,
M. R.
,
2017
, “
Interleaved Continuum-Rigid Manipulation: An Approach to Increase the Capability of Minimally Invasive Surgical Systems
,”
IEEE/ASME Trans. Mech.
,
22
(
1
), pp.
29
40
.
17.
Simaan
,
N.
,
Yasin
,
R. M.
, and
Wang
,
L.
,
2018
, “
Medical Technologies and Challenges of Robot-Assisted Minimally Invasive Intervention and Diagnostics
,”
Annual Rev. Control, Rob., Auto. Syst.
,
1
, pp.
465
490
.
18.
McMahan
,
W.
,
Chitrakaran
,
V.
,
Csencsits
,
M.
,
Dawson
,
D.
,
Walker
,
I.
,
Jones
,
B.
,
Pritts
,
M.
,
Dienno
,
D.
,
Grissom
,
M.
, and
Rahn
,
C.
,
2006
, “
Field Trials and Testing of the Octarm Continuum Manipulator
,”
Proceedings 2006 IEEE International Conference on Robotics and Automation ICRA 2006
,
Orlando, FL
,
May 15–19
, pp.
2336
2341
.
19.
Rone
,
W. S.
,
Saab
,
W.
, and
Ben-Tzvi
,
P.
,
2018
, “
Design, Modeling, and Integration of a Flexible Universal Spatial Robotic Tail
,”
ASME J. Mech. Rob.
,
10
(
4
), p.
041001
.
20.
Rone
,
W. S.
,
Saab
,
W.
,
Kumar
,
A.
, and
Ben-Tzvi
,
P.
,
2019
, “
Controller Design, Analysis, and Experimental Validation of a Robotic Serpentine Tail to Maneuver and Stabilize a Quadrupedal Robot
,”
ASME J. Dyn. Syst. Meas. Control.
,
141
(
8
), p.
081002
.
21.
Liu
,
J.
,
Iacoponi
,
S.
,
Laschi
,
C.
,
Wen
,
L.
, and
Calisti
,
M.
,
2020
, “
Underwater Mobile Manipulation: A Soft Arm on a Benthic Legged Robot
,”
IEEE Rob. Auto. Magaz.
,
27
(
4
), pp.
12
26
.
22.
Gong
,
Z.
,
Chen
,
B.
,
Liu
,
J.
,
Fang
,
X.
,
Liu
,
Z.
,
Wang
,
T.
, and
Wen
,
L.
,
2019
, “
An Opposite-Bending-and-Extension Soft Robotic Manipulator for Delicate Grasping in Shallow Water
,”
Front. Rob. AI
,
6
, p.
00026
.
23.
Simetti
,
E.
,
Casalino
,
G.
,
Torelli
,
S.
,
Sperindé
,
A.
, and
Turetta
,
A.
,
2014
, “
Underwater Floating Manipulation for Robotic Interventions
,”
IFAC Proc. Vol.
,
47
, pp.
3358
3363
.
24.
Khatib
,
O.
,
Yeh
,
X.
,
Brantner
,
G.
,
Soe
,
B.
,
Kim
,
B.
,
Ganguly
,
S.
,
Stuart
,
H.
,
Wang
,
S.
,
Cutkosky
,
M.
,
Edsinger
,
A.
,
Mullins
,
P.
,
Barham
,
M.
,
Voolstra
,
C. R.
,
Salama
,
K. N.
,
L’Hour
,
M.
, and
Creuze
,
V.
,
2016
, “
Ocean One: A Robotic Avatar for Oceanic Discovery
,”
IEEE Rob. Auto. Magaz.
,
23
(
4
), pp.
20
29
.
25.
Manley
,
J. E.
,
Halpin
,
S.
,
Radford
,
N.
, and
Ondler
,
M.
,
2018
, “
Aquanaut: A New Tool for Subsea Inspection and Intervention
,” OCEANS 2018 MTS/IEEE, Charleston, SC, Oct. 22–25, pp.
1
4
.
26.
Wang
,
L.
, and
Simaan
,
N.
,
2019
, “
Geometric Calibration of Continuum Robots: Joint Space and Equilibrium Shape Deviations
,”
IEEE Trans. Rob.
,
35
(
2
), pp.
387
402
.
27.
Simaan
,
N.
,
Xu
,
K.
,
Wei
,
W.
,
Kapoor
,
A.
,
Kazanzides
,
P.
,
Flint
,
P.
, and
Taylor
,
R.
,
2009
, “
Design and Integration of a Telerobotic System for Minimally Invasive Surgery of the Throat
,”
Inter. J. Rob. Res. Spec. Issue Med. Rob. (Special Issue Med. Rob.)
,
28
(
9
), pp.
1134
1153
.
28.
Ma
,
R. R.
,
Odhner
,
L. U.
, and
Dollar
,
A. M.
,
2013
, “
A Modular, Open-Source 3d Printed Underactuated Hand
,”
2013 IEEE International Conference on Robotics and Automation
,
Karlsruhe, Germany
,
May 6–10
,
IEEE
, pp.
2737
2743
.
29.
Ma
,
R.
, and
Dollar
,
A.
,
2017
, “
Yale Openhand Project: Optimizing Open-Source Hand Designs for Ease of Fabrication and Adoption
,”
IEEE Rob. Auto. Magaz.
,
24
(
1
), pp.
32
40
.
30.
Hirose
,
S.
, and
Ma
,
S.
,
1991
, “
Coupled Tendon-Driven Multijoint Manipulator
,”
Proceedings 1991 IEEE International Conference on Robotics and Automation
,
Sacramento, CA
,
Apr. 9–11
, Vol.
2
, pp.
1268
1275
.
31.
Na
,
Y.-m.
,
Lee
,
H.-s.
, and
Park
,
J.-k.
,
2020
, “
Fabrication and Experiment of an Automatic Continuum Robot System Using Image Recognition
,”
ASME J. Mech. Rob.
,
12
(
1
), p.
011017
.
32.
Liu
,
Y.
, and
Ben-Tzvi
,
P.
,
2020
, “
Design, Analysis, and Integration of a New Two-Degree-of-Freedom Articulated Multi-Link Robotic Tail Mechanism
,”
ASME J. Mech. Rob.
,
12
(
2
), p.
021101
.
33.
Paul
,
R. P.
,
1981
,
Robot Manipulators: Mathematics, Programming, and Control: The Computer Control of Robot Manipulators
,
MIT Press
,
Cambridge, MA
.
34.
Mooring
,
B. W.
,
Roth
,
Z. S.
, and
Driels
,
M. R.
,
1991
,
Fundamentals of Manipulator Calibration
,
Wiley-Interscience
,
Hoboken, NJ
.