Numerical solution of a bending-torsion model for elastic rods

Aiming at simulating elastic rods, we discretize a rod model based on a general theory of hyperelasticity for inextensible and unshearable rods. After reviewing this model and discussing topological effects of periodic rods, we prove convergence of the discretized functionals and stability of a corresponding discrete flow. Our experiments numerically confirm thresholds, e.g., for Michell’s instability, and indicate a complex energy landscape, in particular in the presence of impermeability.