Authors: Pietro Pustina, Ebrahim Shahabi, Daniel Feliu-Talegon, Alessandro De Luca, Cosimo Della Santina
Organizations: Department of Computer, Control and Management Engineering, Sapienza University of Rome, Rome, Italy · Department of Cognitive Robotics, Delft University of Technology, 2628 CN Delft, The Netherlands · Institute of Robotics and Mechatronics, German Aerospace Center (DLR), 82234 Oberpfaffenhofen, Germany
Controlling the shape of a continuum soft robot typically requires an accurate dynamic model and actuation of all degrees of freedom. We show that regulating only the actuated coordinates, through collocated shape control, achieves provably stable convergence of those coordinates and, under an explicit compatibility condition, of the entire robot shape. While collocated control is a cornerstone of high-performance motion control in rigid robotics, extending this formulation to continuum soft robots has remained challenging due to the complexity of their dynamics. We present the first general framework for collocated control of continuum soft robots and derive a unified family of controllers, including PD, PID, PsatID, and their counterparts with compensation and cancellation components. The framework unifies existing approaches while introducing new controller designs. In particular, we develop three classes of PD and PID like regulators with local, semi-global, and global stability guarantees, and provide rigorous convergence analyses for each. Extensive experimental validation demonstrates the effectiveness of the proposed methods across different model discretizations and controller parameters. The resulting framework provides practical design guidelines for selecting and implementing controllers with known stability guarantees, without requiring a complete dynamic model of the robot
Figures & tables
Figure 1: Drive a soft robot’s shape by actuating only where you can - collocatedcontrol - with guarantees that scale from model-free practical regulation to model-based global stability guarantees. This approach paves the way for shape regulation in soft robots with stability guarantees, without requiring a full dynamic model.
Figure 2: Block diagram of the collocated shape-regulation framework, incorporating the coordinate transformation, h(q) , a collocated shape regulator, and a decomposition of the system dynamics into actuated and unactuated components.
Figure 3: Summary of the proposed controller families, their performance guarantees, and implementation costs.
Name
Expression of um
Cancellation
ga(θa,θu)+ka(θa,θu)
Gravity cancellation and collocated elasticity compensation
ga(θa,θu)+ka(θad,θu)
Collocated gravity compensation and elasticity cancellation
ga(θad,θu)+ka(θa,θu)
Collocated gravity and elasticity compensation
ga(θad,θu)+ka(θad,θu)
Table 2: Possible model-based terms um for ( 19 ) guaranteeing stability according to Theorem 2 .
Figure 4: Possible unactuated equilibria θud1 and θud2 compatible with the desired θad value of the actuation coordinates.
FAMILY
MEASURES
MODEL
INTEGRAL
CLASS
GUARANTEE
THEOREM
PD
θa
none
—
—
bounded / practical
—
PD + feedforward
θa
Fcd
—
under class
local / global ∗
T1
PD / PID + cancellation
full θ
full potential forces
optional
under class
semi-global ∗
T2–3
P-satI-D
θa
none
saturated integral
elastically dominated
global ∗
T4
Table 3: Summary of controller families and their guarantees. Select the row that matches the available measurements and model knowledge, then read across to identify the applicable controller structure and its corresponding performance guarantee. ∗ Guarantee holds under the stated class assumption; (Sec. 4 ).
Controller
Provably stable
u1=KPθ~a−KDθ˙a
yes
u2=KPθ~a−KDθ˙a+ga(θd)+ka(θd)
yes
u3=KPθ~a−KDθ˙a+ga(θ)+ka(θd)
yes
u4=KPθ~a−KDθ˙a+ga(θ)+ka(θ)
yes
u5=KPθ~a−KDθ˙a+ga(θa,θud)+ka(θd)
no
u6=KPθ~a−KDθ˙a+ga(θad,θu)+ka(θd)
yes
Table 4: Regulators considered in the experiments. For each controller, we report if some form of stability has been proven.
Figure 5: The experimental setup consists of a slender continuum soft arm of Eco-Flex 20 silicone actuated by a pair of tendons tensioned by Dynamixel XM430-W210-T motors. Eight Optitrack Prime x 13 cameras measure the position and orientation of three rings of markers placed along the robot body.
Figure 6: Robot kinematic models considered in the experiments.
Table 5: Curvature models used in the experiments. The variable s denotes the curvilinear abscissa.
Figure 7: Time evolution of the actuation coordinate and the corresponding reference (black dashed lines) for the PD regulators under the PC discretization with 4 DoF, a proportional gain kP=0.5\penalty[\mathrm{N}\text{,}\mathrm{m}] and no mass at the tip.
Figure 8: Time evolution of the actuation coordinate and the corresponding reference (black dashed lines) for the PID/P-satI-D regulators under the PC discretization with 4 DoF, a proportional gain kP=0.5\penalty[\mathrm{N}\text{,}\mathrm{m}] and no mass at the tip.
Figure 9: Time evolution of three controllers (PD ( u1 ), PD+ ( u2 ), and P-satI-D+ ( u9 )) for two different target shapes using PC discretization with 4 DoF and a proportional gain kP=0.5\penalty[\mathrm{N}\text{,}\mathrm{m}] .
Appendix figures & tables13 assets
Supplementary material from the paper’s appendix.
Appendix
Figure 10: (a) Measured and identified step response of the open-loop system, and (b) estimated closed-loop response of the linear plant P(s) under a PID regulator with kP=0.294\penalty[\mathrm{N}\text{,}\mathrm{m}] , kD=0.0313\penalty[\mathrm{N}\text{,}\mathrm{m}\text{,}\mathrm{s}] , and kI=0.689\penalty[\mathrm{N}\text{,}\mathrm{m}\text{,}{\mathrm{s}}^{-1}] .
Figure 11: Block scheme of the closed-loop system.
Figure 12: RMSE, absolute value of the steady state error, and its percentage for the PD regulators. The solid lines represent the mean value for each metric of the PD regulator. For each controller, gain, and metric, the plots show the median, lower and upper quartiles, as well as the minimum and maximum values, and with a colored dot the mean value.
Figure 13: Overshoot, settling time and transient time for the PD regulators. The solid lines represent the mean value for each metric of the PD regulator. For each controller, gain, and metric, the plots show the median, lower and upper quartiles, as well as the minimum and maximum values, and with a colored dot the mean value.
Figure 14: RMSE, absolute value of the steady state error and its percentage for the PID regulators. The solid lines represent the mean value for each metric of the PID regulator. For each controller, gain, and metric, the plots show the median, lower and upper quartiles, as well as the minimum and maximum values, and with a colored dot the mean value.
Figure 15: Overshoot, settling time and transient time for the PID regulators. The solid lines represent the mean value for each metric of the PID regulator. For each controller, gain, and metric, the plots show the median, lower and upper quartiles, as well as the minimum and maximum values, and with a colored dot the mean value.
Figure 16: Mean RMSE, absolute and percentage steady-state errors, overshoot, settling time, and transient time for the PD regulator (left side of the figure) and the PID regulator (right side of the figure). The lines represent the metric values averaged across all model discretizations.
Figure 17: Mean of the RMSE [left], absolute value of the steady state error [center] and its percentage [right] for the analyzed regulators with varying proportional gain (in [Nm]).
Figure 18: Mean value of the overshoot [left], settling time [center] and transient time [right] for the analyzed regulators with varying proportional gain (in [Nm]).
Figure 19: Mean of the RMSE [left], absolute value of the steady state error [center] and its percentage [right] respectively) [center] and its percentage [right] with varying proportional gain (in [Nm]) in the presence of a payload of 0.09\penalty[\mathrm{kg}] attached to the robot tip.
Figure 20: Mean value of the overshoot [left], settling time [center] and transient time [right] for the analyzed regulators with varying proportional gain (in [Nm]) in the presence of a payload of 0.09\penalty[\mathrm{kg}] attached to the robot tip.
Figure 21: Steady state and RMS values of the Cartesian, angular, and task errors for the considered discretizations.
Figure 22: Snapshots of the robot and corresponding dynamic model for the closed-loop system under the P-satI-D controller with kP=0.5\penalty[\mathrm{N}\text{,}\mathrm{m}] .