Passive Stiffness Shaping in Cable-Suspended Aerial Manipulation via Movable Compliant Anchors
Authors: Antonio Franchi, Amr Afifi
Organizations: Robotics and Mechatronics Lab, Faculty of Electrical Engineering, Mathematics & Computer Science, University of Twente, Enschede, The Netherlands · Department of Computer, Control and Management Engineering, Sapienza University of Rome, 00185 Rome, Italy · Department of Information and Computing Sciences, Utrecht University, Utrecht, The Netherlands
Cable-suspended aerial manipulation offers a lightweight architecture for cooperative transportation and physical interaction, yet the passive mechanical response perceived at the load remains insufficiently understood and systematically exploited. This work interprets aerial vehicles as movable compliant anchors and develops a gravity-aware quasi-static theory for predicting and shaping the passive Cartesian stiffness of a suspended load. The formulation applies to an arbitrary number of aerial vehicles connected to a point load by taut, straight, inextensible cables. At a selected gravity-loaded equilibrium, aerial-anchor compliance and transverse cable geometric compliance combine in series within each leg, while the leg stiffnesses act in parallel on the load. For isotropic aerial-anchor behavior, each leg is exactly equivalent to a virtual unilateral elastic cable, revealing an axial--transverse stiffness decomposition governed by the equilibrium tension. These results define a nonlinear map from commanded-anchor configuration to passive load stiffness, whose differential enables local constraint-preserving shaping through anchor repositioning. A dynamic rigid-body validation framework with nonlinear vehicle control, elastic-damped tendons, and environmental contact is defined to assess when and to what extent the derived stiffness remains predictive beyond the assumptions of the analytical model.
Figures & tables
Fig. 1 : System architecture and passive-stiffness-shaping principle. Each taut cable connects the load point P to an actual aerial-anchor point Pa,i , while the finite closed-loop stiffness of aerial vehicle i relates this point to its commanded anchor Pd,i . Gravity determines the equilibrium cable directions and tensions. Repositioning the commanded anchors changes these operating-point quantities and thereby reshapes the passive Cartesian stiffness ellipsoid at the load. The three-vehicle drawing represents the general n -vehicle system.
Fig. 2 : Physical and virtual representations of one cable–anchor leg. The physical leg combines aerial-anchor compliance and the transverse geometric compliance of the taut inextensible cable in series. Under isotropic aerial-anchor stiffness, the leg is quasi-statically equivalent to a virtual unilateral elastic cable connecting the load directly to the commanded anchor point.
Fig. 3 : Axial–transverse eigenstructure of an isotropic cable–anchor leg. The Euclidean-identified stiffness operator has the simple eigenvalue ki along the cable direction and the eigenvalue γi , with algebraic and geometric multiplicity two, on the transverse plane. Increasing the operating-point tension changes the leg contribution from nearly rank one toward isotropic.
Fig. 4 : Gravity-aware stiffness pipeline and parallel intrinsic and coordinate descriptions. The aerial-anchor setpoints determine a selected gravity-loaded equilibrium branch and thereby the equilibrium cable directions and tensions. These quantities determine the leg and total stiffness maps. Their coordinate representations yield the stiffness-coordinate vector and its Jacobian for sensitivity analysis and control. Gravity affects stiffness through the selected equilibrium rather than through an additive stiffness term.
Fig. 5 : Constraint-preserving local stiffness-regulation loop.
Parameter
Value
System: number of quadrotors and cables
4
Quadrotor: mass
1.28kg
Quadrotor: principal inertia
diag(0.015,0.015,0.007)kgm2
Payload: mass
2.00kg
Payload: principal inertia
diag(0.05,0.05,0.05)kgm2
Translational control: position gain Kp
diag(12,12,14)N/m
TABLE I : MuJoCo model and controller parameters.
Phenomenon
Analytical representation
MuJoCo representation
Remaining simulation-to-reality gap
Payload
Point load with translational displacement
Free rigid body with mass, inertia, attitude, and four spatially separated attachment sites
Structural flexibility, uncertain inertial parameters, and unmodeled payload aerodynamics
Aerial vehicles
Movable anchors with prescribed linear Cartesian stiffness
Rigid quadrotors with nonlinear SE(3) tracking and finite translational gains
Motor dynamics, thrust uncertainty, battery effects, and aerodynamic interaction
Anchor compliance
Symmetric positive-definite linear stiffness map
Controller-induced displacement under cable loading
Gain variation, saturation, delays, estimation errors, and off-equilibrium nonlinearities
Cables
Taut, straight, massless, and inextensible
Elastic-damped spatial tendons with finite axial stiffness and rest length
Distributed mass, sagging, bending, drag, transverse vibration, and cable contact
Dynamics
Quasi-static evolution along a selected equilibrium branch
Coupled vehicle–payload transients with inertia and dissipation
Unmodeled high-frequency dynamics and hardware-dependent settling behavior
Contact
Excluded from the stiffness derivation
Optimization-based contact between the payload and environment
Surface compliance, friction uncertainty, impact dynamics, and geometry imperfections
TABLE II : Relationship among the analytical model, the MuJoCo implementation, and the corresponding physical system.
Fig. 6 : Parity-plot comparison across the selected stiffness and deformation quantities.
Fig. 7 : Stiffness tracking plots for 2 of the 4 cases used in zigzag sliding
Fig. 8 : Stiffness tracking plots for 2 of the 4 cases used in zigzag sliding
Fig. 9 : Cable tensions during stiffness regulation.
Fig. 10 : Commanded-anchor velocities during stiffness regulation.
Parameter
Value
Role
Payload regulation gain λp
1.0
Weight on payload-position error
Stiffness regulation gain λK
1.0
Weight on stiffness-tracking error
Payload objective weight wp
1.0
Relative weight of payload term in the QP cost
Stiffness objective weight wK
25.0
Relative weight of stiffness term in the QP cost
Velocity regularization εv
10−4
Quadratic penalty on anchor-velocity magnitude
Admissibility gain ηT
1.0
Gain in the minimum-tension inequality
TABLE III : Default parameters of the implemented QP anchor-setpoint controller.
Fig. 11 : Schematic of the zigzag guided-sliding task, including the longitudinal direction of progression and the lateral contact direction.
Fig. 12 : Longitudinal and lateral tracking errors during guided sliding.
Fig. 13 : Commanded-anchor formations used to generate the four passive-stiffness profiles for the guided-sliding comparison.
Department of Informatics, Bioengineering, Robotics and Systems Engineering, Università degli Studi di Genova, Via all’Opera Pia 13, 16145, Genoa, Italy