  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.
