Mixed finite element discretization of intrinsic geometrically exact beams for multiplier-free multibody dynamics
Abstract
The Reissner-Simo and Hodges models are two equivalent continuous descriptions of finite-strain beam dynamics. The Reissner-Simo formulation uses displacements and rotations, while the Hodges formulation is intrinsic and avoids both variables. Although equivalent in theory, the two approaches behave differently after discretization and offer distinct numerical advantages. In this work, we develop a structure-preserving discretization of the intrinsic planar beam formulation. Because the intrinsic equations involve linear differential operators, both kinematic and dynamic boundary conditions can be imposed naturally using mixed finite elements. The resulting formulation also enables multibody systems to be assembled without algebraic constraints, avoiding the stiff differential-algebraic equations typically introduced by kinematic constraints. We demonstrate the approach on different examples, also showing that closed kinematic loops can be modeled without algebraic constraints. The resulting interconnected systems retain a port-Hamiltonian structure, with all nonlinearities confined to the interconnection operator. This structure allows exact energy preservation when combined with implicit midpoint time integration. Furthermore the scheme appears to require less Newton iterations compared to existing energy preserving scheme.
Related
- Input-to-State Stability Framework for Fully Distributed Primal-Dual Dynamics for Quadratic GNEPs Without Multiplier Consensus
- Fischer-Burmeister-based Multiplier Dynamics in Inequality-Constrained Optimization
- Euler Time-Discretization of Funnel-Controlled Dynamics
- Large time-step discretization of adiabatic quantum dynamics
- Latent Ornstein–Uhlenbeck-Driven Log-ARCH Models: Continuous-Time Stochastic Dynamics, Exact Discretization, and Long-Span QML Theory
BibTeX
@article{brugnoli2026mixed,
title = {{Mixed finite element discretization of intrinsic geometrically exact beams for multiplier-free multibody dynamics}},
author = {Andrea Brugnoli and Phi Kinon and Francesco Sanfedino and Peter Betsch and Olivier Bauchau},
journal = {Mechanism and Machine Theory},
year = {2026},
month = oct,
doi = {10.1016/j.mechmachtheory.2026.106654},
url = {https://hal.science/hal-05805375v1},
note = {CC BY-NC-ND 4.0},
}
RIS
TY - JOUR
TI - Mixed finite element discretization of intrinsic geometrically exact beams for multiplier-free multibody dynamics
AU - Andrea Brugnoli
AU - Phi Kinon
AU - Francesco Sanfedino
AU - Peter Betsch
AU - Olivier Bauchau
PY - 2026
DA - 2026/10/19
JO - Mechanism and Machine Theory
DO - 10.1016/j.mechmachtheory.2026.106654
UR - https://hal.science/hal-05805375v1
AB - The Reissner-Simo and Hodges models are two equivalent continuous descriptions of finite-strain beam dynamics. The Reissner-Simo formulation uses displacements and rotations, while the Hodges formulation is intrinsic and avoids both variables. Although equivalent in theory, the two approaches behave differently after discretization and offer distinct numerical advantages. In this work, we develop a structure-preserving discretization of the intrinsic planar beam formulation. Because the intrinsic equations involve linear differential operators, both kinematic and dynamic boundary conditions can be imposed naturally using mixed finite elements. The resulting formulation also enables multibody systems to be assembled without algebraic constraints, avoiding the stiff differential-algebraic equations typically introduced by kinematic constraints. We demonstrate the approach on different examples, also showing that closed kinematic loops can be modeled without algebraic constraints. The resulting interconnected systems retain a port-Hamiltonian structure, with all nonlinearities confined to the interconnection operator. This structure allows exact energy preservation when combined with implicit midpoint time integration. Furthermore the scheme appears to require less Newton iterations compared to existing energy preserving scheme.
SN - 0094-114X
ER -