Mathematics and Statistics Vol. 10(1), pp. 25 - 34
DOI: 10.13189/ms.2022.100103
Reprint (PDF) (668Kb)

A Simulation of an Elastic Filament Using Kirchhoff Model

Saimir Tola 1, Alfred Daci 1,*, Gentian Zavalani 1,2
1 Faculty of Mathematical Engineering and Physical Engineering, Polytechnic University of Tirana, Tiranë, Albania
2 Center for Advanced Systems Understanding (CASUS), D-02826 Görlitz, Germany


This paper presents numerical simulations and comparisons between different approaches concerning elastic thin rods. Elastic rods are ideal for modeling the stretching, bending, and twisting deformations of such long and thin elastic materials. The static solution of Kirchhoff's equations [2] is produced using ODE45 solver where Kirchhoff and reference system equations are combined instantaneously. Solutions using formulations are based on Euler's elastica theory [1] which determines the deformed centerline of the rod by solving a boundary-value problem, on the Discreet Elastic Rod method using Bishop frame (DER) [5,6] which is based on discrete differential geometry, it starts with a discrete energy formulation and obtains the forces and equations of motion by taking the derivative of energies. Instead of discretizing smooth equations, DER solves discrete equations and obeys geometrical exactness. Using DER we measure torsion as the difference of angles between the material and the Bishop frame of the rod so that no additional degree of freedom is needed to represent the torsional behavior. We found excellent agreement between our Kirchhoff-based solution and numerical results obtained by the other methods. In our numerical results, we include simulation of the rod under the action of the terminal moment and illustrations of the gravity effects.

Kirchhoff Equation, ODE45, Fernet–Serret Frame, Bishop Frame, Discreet Elastic Rod

Cite This Paper in IEEE or APA Citation Styles
(a). IEEE Format:
[1] Saimir Tola , Alfred Daci , Gentian Zavalani , "A Simulation of an Elastic Filament Using Kirchhoff Model," Mathematics and Statistics, Vol. 10, No. 1, pp. 25 - 34, 2022. DOI: 10.13189/ms.2022.100103.

(b). APA Format:
Saimir Tola , Alfred Daci , Gentian Zavalani (2022). A Simulation of an Elastic Filament Using Kirchhoff Model. Mathematics and Statistics, 10(1), 25 - 34. DOI: 10.13189/ms.2022.100103.