No abstract

SUMMARY An angular momentum and energy‐conserving time integration algorithm for rigid body rotation is formulated in terms of the quaternion parameters and the corresponding four‐component conjugate momentum vector via Hamilton's equations. The introduction of an extended mass matrix leads to a symmetric set of eight state‐space equations of motion. The extra inertial parameter serves as a multiplier on the kinematic constraint, and it is demonstrated that convergence characteristics are improved by selecting this parameter somewhat larger than the inertial moments. External loads enter these equations via the set of momentum equations. Initially, the normalization of the quaternion array is introduced via a Lagrange multiplier. However, this Lagrange multiplier can be expressed explicitly in terms of the gradient of the external load potential, and elimination of the Lagrange multiplier from the final format leaves only an explicit projection applied to the external load potential gradient. An algorithm is developed by forming a finite increment of the Hamiltonian. This procedure identifies the proper selection of increments and mean values, and leads to an algorithm with conservation of momentum and energy. Implementation, conservation properties, and accuracy of the algorithm are illustrated by examples with a flying box and a spinning top. Copyright © 2012 John Wiley & Sons, Ltd.

No abstract

No abstract

SUMMARYA time integration algorithm is developed for the equations of motion of a flexible body in a rotating frame of reference. The equations are formulated in a hybrid state-space, formed by the local displacement components and the global velocity components. In the spatial discretization the local displacements and the global velocities are represented by the same shape functions. This leads to a simple generalization of the corresponding equations of motion in a stationary frame in which all inertial effects are represented via the classic global mass matrix. The formulation introduces two gyroscopic terms, while the centrifugal forces are represented implicitly via the hybrid state-space format. An angular momentum and energy conserving algorithm is developed, in which the angular velocity of the frame is represented by its mean value. A consistent algorithmic damping scheme is identified by applying the conservative algorithm to a decaying response, which is rendered stationary by an increasing exponential factor that compensates the decay. The algorithmic damping is implemented by introducing forward weighting of the mean values appearing in the algorithm. Numerical examples illustrate the simplicity and accuracy of the algorithm.