[Paper Review] Modelling Developable Ribbons Using Ruling Bending Coordinates

Zherong Pan, Jin Huang+1 more|arXiv (Cornell University)|Mar 13, 2016

Advanced Materials and Mechanics17 references3 citations

TL;DR

This paper introduces a novel dynamic simulation method for developable ribbons using ruling bending coordinates, parameterizing ribbon shapes via crease angles and bending angles along a centerline. The approach ensures exact isometric developability, enables stable large-timestep integration via fully implicit optimization, and achieves high-fidelity simulation with minimal elements, outperforming FEM in stability and robustness to extreme deformation and collisions.

ABSTRACT

This paper presents a new method for modelling the dynamic behaviour of developable ribbons, two dimensional strips with much smaller width than length. Instead of approximating such surface with a general triangle mesh, we characterize it by a set of creases and bending angles across them. This representation allows the developability to be satisfied everywhere while still leaves enough degree of freedom to represent salient global deformation. We show how the potential and kinetic energies can be properly discretized in this configuration space and time integrated in a fully implicit manner. The result is a dynamic simulator with several desirable features: We can model non-trivial deformation using much fewer elements than conventional FEM method. It is stable under extreme deformation, external force or large timestep size. And we can readily handle various user constraints in Euclidean space.

Motivation & Objective

• To address the instability and locking issues in FEM-based simulations of developable ribbons due to approximate developability and artificial constraints.
• To develop a configuration space that exactly represents the true developable shape space, eliminating the need for additional constraints.
• To enable stable, large-timestep dynamic simulation of ribbons under extreme forces and complex constraints.
• To support flexible user constraints and robust collision handling in a unified optimization framework.
• To reduce computational cost and improve scalability compared to conventional FEM while maintaining geometric accuracy.

Proposed method

• Parameterize ribbon shape using a set of creases along the centerline and bending angles across them, ensuring exact isometric developability.
• Define a configuration space where world-space vertex positions are reconstructed analytically from ruling bending coordinates, avoiding mesh-based approximations.
• Discretize kinetic and potential energy using auxiliary variables for world-space positions, enabling analytical gradient computation.
• Formulate time integration as a single optimization problem in the configuration space, using fully implicit time stepping.
• Introduce a closeness metric $ E $ to enforce consistency between simulated and target configurations, improving stability under large timesteps.
• Leverage adjoint methods to accelerate Hessian and gradient evaluation, mitigating scalability issues for long ribbons.

Experimental results

Research questions

• RQ1 Can a configuration space be designed that exactly captures the developable shape space of ribbons, avoiding artificial constraints?
• RQ2 How can kinetic and potential energy be discretized in a way that supports stable, large-timestep integration in a reduced configuration space?
• RQ3 Can collision and contact handling be seamlessly integrated into a non-mesh-based, optimization-driven simulation framework?
• RQ4 How does the method compare in accuracy and stability to FEM-based approaches under extreme deformation and large timesteps?
• RQ5 What are the scalability limits of the method, and how can performance be improved for long ribbons?

Key findings

• The proposed method achieves geometric accuracy by ensuring all configurations are exactly developable, eliminating locking and stiffness issues common in FEM.
• The simulator remains stable under large timesteps ($ h = 0.01 $) and extreme external forces, as demonstrated in falling and impact simulations.
• The method correctly captures complex deformations such as centrifugal-induced triangular ruling patterns, which are missed by simplified lumped-mass kinetic energy models.
• Validation against FEM with remeshing shows high accuracy, with discrepancies visualized in Fig. 10, confirming the fidelity of the formulation.
• The use of adjoint methods reduces the computational cost of gradient evaluation from quadratic to manageable scaling, improving performance for long ribbons.
• Robust collision handling is achieved without mesh distortion, even in complex scenarios like double ribbon chains, where conventional solvers fail.