[luwang00]

This PR is ready for merging.

Feature or improvement description
This PR introduces slip supports to MoorDyn. A slip support can be attached to any Point type (free, fixed, coupled, or attached to a Body) and acts like a hoop that lines pass through. A slip support constrains line motion in the normal direction only (through stiffness and damping), while allowing the line to slip through the support along the tangent direction, resisted only by friction.

The formulation uses a smooth B-spline kernel to distribute loads across several neighboring nodes, giving a smooth load handoff as the line nodes slip past the support, without requiring remeshing. The normal constraint is enforced through stiffness and damping based on the normal offset and relative normal velocity between the line contact point and the slip support. In the tangential direction, the line experiences friction modeled with either regularized Coulomb friction with a Stribeck transition from static to kinetic friction, or the single-state LuGre model, which captures true sticking.

The contact point kinematics are also interpolated from neighboring line nodes using the same B-spline kernel, so the contact point lies on a smooth curve that can be offset from the polyline representation of the mooring line. The offset is larger for sharper bends and decreases with a smaller kernel width (a user input), at the cost of more ripple as nodes pass the support. Because the kernel width is specified relative to the segment length, refining the line discretization reduces the offset while keeping the same smoothness. Using the same kernel for interpolating the contact kinematics and spreading the loads ensures adjoint consistency, so the contact spring and damper act passively on the line.

With the current implementation, the force and velocity projections are consistent for the normal damping and the tangential friction. The node forces and the velocities they act against use the same per-node projections. As a result, damping and Coulomb friction are dissipative, and LuGre friction is passive. The spring force is different. It is computed from the gap at the contact point, normal to the line tangent there. When it is spread to the nodes, the component along each node's local tangent is dropped. Without that projection, the normal support would push the line along its own axis and cause spurious snatching. The trade-off is that the spring force is not exactly conservative. It cannot add energy through axial motion, but it is not the gradient of a potential either.

Each line can have multiple slip supports, and each Point can serve as a slip support for multiple lines. A free Point used as a slip support must have nonzero mass, and a very light Point may require a reduced, user-specified stiffness. If a line end moves past a slip support, the line is permanently released from that support.

Linearization is supported, and new output channels report the contact location, slip velocity, normal and friction forces, bristle deflection (LuGre), and the reaction force on the support Point.

Impacted areas of the software
MoorDyn, docs

Additional supporting information
Detailed user and theory guide for slip support can be found in
docs/source/user/moordyn/moordynf_features.rst
docs/source/user/moordyn/slip_support_theory.rst

Example results
The example below was simulated using MoorDyn driver. A slack line is draped over a slip support, the red cross. The left end of the line is anchored, and the right end follows prescribed motion. Initially, the right end is at rest. At t=20 s, the right end starts moving toward +x at a constant speed of 0.5 m/s until t = 132 s, at which point, the motion stops. During this process, the line is pulled to the right until taut, with the line nodes sliding past the slip support smoothly.

MD_SlideSupport_Demo1.mp4

In addition to the normal force at the slip support, friction based on the LuGre model is also included. The time series of the friction force is shown below, including a zoomed-in figure around 20 s when the sliding started.

The friction is zero initially when the system is at rest. This is a direct consequence of the dynamic initialization, which produces a friction-free configuration at t=0 by design. After the dynamic initialization, the bristle deflection state of the LuGre model is reset to zero, so that there is no elastic stiction force.

There is a spike in friction at about 20 s when the right end of the line starts to get pulled. For the most part of the time series, the line is sliding at about 0.5 m/s, well above the Stribeck velocity of 0.1 m/s defined in the input file. Therefore, the friction is capped at the dynamic friction limit that starts to grow rapidly after about 120 s. This is the point where the line starts to get pulled taut, and the normal force at the slip support rises rapidly. At 132 s, the pulling stops, and the line stiction at the slip support reestablishes. A spike in friction followed by oscillation and slow relaxation can be observed.

Focusing around 20 s when the line just starts to move, a spike in friction past the dynamic friction limit can be seen due to stiction. The max friction does not quite reach the static friction limit. This rate-dependent slip is a known behavior of the LuGre model, where a fast start can undercut the peak static friction. Right after the slipping starts, the friction can temporarily remain above the dynamic friction limit and takes a finite amount of time to drop down to the limit. This is all related to the friction lag effect.

The new slip support should also correctly interact with other existing MoorDyn features, such as line failure. In the example below, we continue to pull the right end of the chain after 132 s. The anchor point is set to detach after reaching 1e5 N tension. The line continues to move past the slip support after detaching from the anchor until it is fully released from the support and drops to the seabed.

MD_SlideSupport_Demo1.mp4

Note that in the two examples above, the slip support is attached to a fixed Point. It is also possible to use a free or coupled Point as slip support. In the latter case, the slip support becomes attached to the rest of the OpenFAST structure like a fairlead. This is shown in the example below. The horizontal cylinder floating on the free surface is modeled in HydroDyn and ElastoDyn. The two ends of the cylinder are both coupled to MoorDyn slip supports. The nonlinear hydrodynamic loading causes the cylinder to drift downstream, sliding along the two mooring lines that act like guide ropes.

HMember_MDSlipSupport.mp4

Generative AI usage
This development is heavily supported by AI both during formulation and implementation.

• The mathematical formulation has gone through careful human review.
• The code was not yet checked. Example cases have been developed to demonstrate that the results are plausible. The new codes are mostly isolated in the new module MoorDyn_SlipSupport.f90 with minimal interface with the existing MoorDyn codebase; therefore, it is extremely unlikely it will cause any issue with existing features.

Models/tools:

• Co-authored-by: Anthropic Claude claude@anthropic.com
• Assisted by: GitHub Copilot support@github.com

Test results, if applicable
The two examples above, a line over a slip support pulled taut and a floating horizontal cylinder with two guide ropes, are both added as new r-tests. The floating horizontal cylinder glue-code test also serves as an additional r-test for PR #3489.
The two new tests are md_SlipSupport_LuGre and Floating_Cylinder_Guide_Ropes.

• r-test branch merging required

[Copilot]

🟡 Changes recommended

Unresolved critical defects can inject energy through damping and friction, retain incorrect support releases after state restoration, and cause out-of-bounds output access.

6 open issues

• Reject support IDs below 1 to prevent out-of-bounds access · New
• Use consistent velocity projection and force spreading for damping and friction · New
• Restore support release status during correction iterations · New
• Recompute viscoelastic extensions after updating support geometry · New
• Match the full slip-support section header instead of any “slip” · New
• Preserve small nonzero friction inputs with scientific notation · New

What changed in this PR

Adds slip supports to MoorDyn-F, allowing lines to slide through Point-based guides with transverse restraint and friction.

Changes:

• Implements contact forces, Coulomb and LuGre friction, line release, outputs, and linearization integration.
• Adds Python input/output support and registers standalone and coupled regression cases.
• Documents inputs, usage, and theory, and updates both build systems.

File	Description
vs-build/​modules/​MoorDyn.vfproj	Includes the slip-support module in Visual Studio builds.
reg_tests/​CTestList.cmake	Registers two slip-support regression cases.
openfast_io/​openfast_io/​FAST_writer.py	Writes slip-support input tables.
openfast_io/​openfast_io/​FAST_reader.py	Reads slip-support input tables.
modules/​moordyn/​src/​MoorDyn.f90	Integrates parsing, initialization, forces, release, and linearization.
modules/​moordyn/​src/​MoorDyn_Types.f90	Adds slip-support types and copy/serialization routines.
modules/​moordyn/​src/​MoorDyn_SlipSupport.f90	Implements contact search, force distribution, friction, and release.
modules/​moordyn/​src/​MoorDyn_Registry.txt	Defines slip-support data and state fields.
modules/​moordyn/​src/​MoorDyn_Point.f90	Adds support reactions to Point forces.
modules/​moordyn/​src/​MoorDyn_Line.f90	Includes slip-support forces in line dynamics.
modules/​moordyn/​src/​MoorDyn_IO.f90	Adds slip-support output channels.
modules/​moordyn/​CMakeLists.txt	Includes the new module in CMake builds.
docs/​source/​user/​moordyn/​slip_support_theory.rst	Documents the contact and friction formulation.
docs/​source/​user/​moordyn/​moordynf_features.rst	Describes inputs, outputs, usage, and linearization.
docs/​source/​user/​moordyn/​index.rst	Links the new MoorDyn-F documentation.
docs/​source/​user/​api_change.rst	Records the optional slip-support input section.

_{Review effort: Balanced}