Conference paper
Rigid Body Time Integration by Convected Base Vectors with Implicit Constraints
A conservative time integration algorithm based on a convected set of orthonormal base vectors is presented. The equations of motion are derived from an extended Hamiltonian formulation, combining the components of the three base vectors with a set of orthonormality constraints. The particular form of the kinetic energy used in the present formulation is deliberately chosen to correspond to a rigid body rotation, and the orthonormality constraints are introduced via the equivalent Green strain components of the base vectors.
The particular form of the extended inertia tensor used here implies a set of orthogonality relations between the base vector components and their conjugate momentum components. These orthogonality relations permit explicit elimination of the Lagrange multipliers associated with the constraints, leading to a projected form of the dynamic equation without explicit algebraic constraints.
The differential equations of motion are recast into discrete form using a suitable combination of mean values and increments, which is identified by considering a finite increment of the Hamiltonian. Examples illustrate the accuracy and conservation properties of the algorithm.
Language: | English |
---|---|
Year: | 2013 |
Proceedings: | 4th International Conference on Computational Methods in Structural Dynamics and Earthquake Engineering |
Types: | Conference paper |
ORCIDs: | Krenk, Steen and Nielsen, Martin Bjerre |