Marine systems simulation
Rigid body cable
+ Collaboration diagram for Rigid body cable:

Cable model with bending/torsional stiffness and winch-functionality. The cable is discretized as "6 degrees of freedom" rigid cylinders which are connected with axial and angular constraints. For high bending stiffness the cable will work as a truss/beam.

The constraints are regularized with a compliant variation of the Baumgarte method. In effect the cylinders can be thought of as connected with axial and angular springs like a traditional spring-mass-damper model, but without some of the numerical instability problems associated with high stiffness cables.

The behavior of the springs are defined with three parameters; \(\alpha\), \(\beta\) and \(\epsilon\). \(\beta\) controls the rate at which the constraint seeks towards equilibrium, \(\alpha\) is the damping factor, and \(\epsilon\) is the constraint compliance.

The effective spring stiffness is a combination of \(\beta\), \(\epsilon\), and the mathematical formulation of the constraint function.

For the axial spring the effective restoring force is: \( F = \beta^2/\epsilon*dL \) where \(dL\) is the distance between two successive cylinder endpoints.

For the angular spring the effective restoring moment is: \(M = \beta^2/\epsilon*sin(\theta/2)\) where theta is the difference in orientation angle between two successive cylinders. Note that the bending and the torsional stiffness can be set independently.

\(\alpha\) less than \(\beta\) means the system is underdamped, \(\alpha\) equal to \(\beta\) means critical damping, and \(\alpha\) larger than Beta means overdamped. The effective damping scales with \(\alpha/\epsilon\)

A fairly robust technique for adjusting the parameters when using an explicit integration scheme is as follows:

The cable model is equipped with virtual winch functionality which makes it possible to retract/extend the cable from the end points. The cable is manipulated by the retracted length and retraction speed at the end points A and B. The retracted length us used to change the cable length, while the retracted speed is used to satisfy the dynamics of the internal cable equations at the end points. The current extended length of the cable is \( L_{cable} = L_{total} - L_{retracted_A} - L_{retracted_B}\)

The definition of retracted length and retracted speed vs. cable length

Input ports

Name Width Description
posA 3 position of end-connection point A
posB 3 position of end-connection point B
velA 3 velocity of end-connection point A
velB 3 velocity of end-connection point B
retractedLengthA 1 length of cable retracted on side A
retractedLengthB 1 length of cable retracted on side B
retractedSpeedA 1 retraction rate of cable retracted on side A
retractedSpeedB 1 retraction rate of cable retracted on side B

Output ports

Name Width Description
forceA 3 Reaction force on side A
forceB 3 Reaction force on side B

Configuration parameters

Name Width Description
Length 1 total length of the cable [m]
numElements 1 number of rigid elements in the cable [#]
AlphaN 1 axial damping factor [s^-1]
BetaN 1 axial restoring factor [s^-1]
EpsilonN 1 axial compliance factor [kg^-1]
AlphaM 1 bending damping factor [s^-1]
BetaM 1 bending restoring factor [s^-1]
EpsilonM 1 bending compliance factor [kg^-1*m^-2]
AlphaT 1 torsional damping factor [s^-1]
BetaT 1 torsional restoring factor [s^-1]
EpsilonT 1 torsional compliance factor [kg^-1*m^-2]
Weight 1 weight of the cable [kg/m]
Radius 1 radius of the cable. [m]

NB: Copied from segment beneath for safe-keeping

Parameters
ParameterDefaultComment
Length-total length of the cable [m]
numElements-number of rigid elements in the cable [#]
AlphaN-axial damping factor [s^-1]
BetaN-axial restoring factor [s^-1]
EpsilonN-axial compliance factor [kg^-1]
AlphaM-bending damping factor [s^-1]
BetaM-bending restoring factor [s^-1]
EpsilonM-bending compliance factor [kg^-1*m^-2]
AlphaT-torsional damping factor [s^-1]
BetaT-torsional restoring factor [s^-1]
EpsilonT-torsional compliance factor [kg^-1*m^-2]
Weight-weight of the cable [kg/m]
Radius-radius of the cable. [m]