Numerical code for time dependent solution from The θ-formulation of the two-dimensional elastica: buckling and boundary layer theory

The numerical code used for this study is our own implementation (in C) of the method introduced by M. Gazzola et al., “Forward and inverse problems in the mechanics of soft filaments”, Royal Society open science 5, 171628 (2018). A python version of the code developed by the authors of this method is available on Github (https://protect-eu.mimecast.com/s/cnaJCy8A4IBJVLTZGHC0?domain=github.com).;Here we provide a typical script of our version of the code that allows to reproduce the numerical data presented in our study.; The file “cosseratFilament.c” that calls the functions in “functions.c”. The code can be compiled using the provided makefile with a standard C compiler (e.g gcc).; Here the script generates a peanut-like elastica as those studied in our paper. Starting from a circular initial condition, we constrain the distance d between two points of the elastica belonging to the same diameter of the initial circle. The imposed distance is reduced by small increment, (waiting enough time at each increment for the elastica to reach meachanical equilibrium) from the initial distance d=2R to d=0. During the simulation, the position of the vertices along the centerline and the orientation matrix of the elements of the rod are saved periodically.