Discrete formulation for the dynamics of rods deforming in space

Ana Casimiro, César Rodrigo

Research output: Contribution to journalArticlepeer-review


The movement of rods in an Euclidean space can be described as a field theory on a principal bundle. The dynamics of a rod is governed by partial differential equations that may have a variational origin. If the corresponding smooth Lagrangian density is invariant by some group of transformations, there exist the corresponding conserved Noether currents. Generally, numerical schemes dealing with PDEs fail to reflect these conservation properties. We describe the main ingredients needed to create, from the smooth Lagrangian density, a variational principle for discrete motions of a discrete rod, with the corresponding conserved Noether currents. We describe all geometrical objects in terms of elements on the linear Atiyah bundle using a reduced forward difference operator. We show how this introduces a discrete Lagrangian density that models the discrete dynamics of a discrete rod. The presented tools are general enough to represent a discretization of any variational theory in principal bundles, and its simplicity allows us to perform an iterative integration algorithm to compute the discrete rod evolution in time, starting from any predefined configurations of all discrete rod elements at initial times.

Original languageEnglish
Article number092901
JournalJournal Of Mathematical Physics
Issue number9
Publication statusPublished - 1 Sept 2019


Dive into the research topics of 'Discrete formulation for the dynamics of rods deforming in space'. Together they form a unique fingerprint.

Cite this