Stress tensors, Riemannian metrics and the alternative descriptions in elasticity

N/

Symplectic-energy-momentum preserving variational integrators

The purpose of this paper is to develop variational integrators for conservative mechanical systems that are symplectic and energy and momentum conserving. To do this, a space–time view of variational integrators is employed and time step adaptation is used to impose the constraint of conservation of energy. Criteria for the solvability of the time steps and some numerical examples are given

Stability Analysis of a Rigid Body with Attached Geometrically Nonlinear Rod by the Energy-Momentum Method

This paper applies the energy-momentum method to the problem of nonlinear stability of relative equilibria of a rigid body with attached flexible appendage in a uniformly rotating state. The appendage is modeled as a geometrically exact rod which allows for finite bending, shearing and twist in three dimensions. Application of the energy-momentum method to this example depends crucially on a special choice of variables in terms of which the second variation block diagonalizes into blocks associated with rigid body modes and internal vibration modes respectively. The analysis yields a nonlinear stability result which states that relative equilibria are nonlinearly stable provided that; (i) the angular velocity is bounded above by the square root of the minimum eigenvalue of an associated linear operator and, (ii) the whole assemblage is rotating about the minimum axis of inertia

Hamiltonian reductions of the one-dimensional Vlasov equation using phase-space moments

We consider Hamiltonian closures of the Vlasov equation using the phase-space moments of the distribution function. We provide some conditions on the closures imposed by the Jacobi identity. We completely solve some families of examples. As a result, we show that imposing that the resulting reduced system preserves the Hamiltonian character of the parent model shapes its phase space by creating a set of Casimir invariants as a direct consequence of the Jacobi identity

Discrete Routh Reduction

This paper develops the theory of abelian Routh reduction for discrete mechanical systems and applies it to the variational integration of mechanical systems with abelian symmetry. The reduction of variational Runge-Kutta discretizations is considered, as well as the extent to which symmetry reduction and discretization commute. These reduced methods allow the direct simulation of dynamical features such as relative equilibria and relative periodic orbits that can be obscured or difficult to identify in the unreduced dynamics. The methods are demonstrated for the dynamics of an Earth orbiting satellite with a non-spherical $J_2$ correction, as well as the double spherical pendulum. The $J_2$ problem is interesting because in the unreduced picture, geometric phases inherent in the model and those due to numerical discretization can be hard to distinguish, but this issue does not appear in the reduced algorithm, where one can directly observe interesting dynamical structures in the reduced phase space (the cotangent bundle of shape space), in which the geometric phases have been removed. The main feature of the double spherical pendulum example is that it has a nontrivial magnetic term in its reduced symplectic form. Our method is still efficient as it can directly handle the essential non-canonical nature of the symplectic structure. In contrast, a traditional symplectic method for canonical systems could require repeated coordinate changes if one is evoking Darboux' theorem to transform the symplectic structure into canonical form, thereby incurring additional computational cost. Our method allows one to design reduced symplectic integrators in a natural way, despite the noncanonical nature of the symplectic structure.Comment: 24 pages, 7 figures, numerous minor improvements, references added, fixed typo

Frictional Collisions Off Sharp Objects

This work develops robust contact algorithms capable of dealing with multibody nonsmooth contact geometries for which neither normals nor gap functions can be defined. Such situations arise in the early stage of fragmentation when a number of angular fragments undergo complex collision sequences before eventually scattering. Such situations precludes the application of most contact algorithms proposed to date

Variational integrators, the Newmark scheme, and dissipative systems

Variational methods are a class of symplectic-momentum integrators for ODEs. Using these schemes, it is shown that the classical Newmark algorithm is structure preserving in a non-obvious way, thus explaining the observed numerical behavior. Modifications to variational methods to include forcing and dissipation are also proposed, extending the advantages of structure preserving integrators to non-conservative systems

A block diagonalization theorem in the energy-momentum method

We prove a geometric generalization of a block diagonalization theorem first found by the authors for rotating elastic rods. The result here is given in the general context of simple mechanical systems with a symmetry group acting by isometries on a configuration manifold. The result provides a choice of variables for linearized dynamics at a relative equilibrium which block diagonalizes the second variation of an augmented energy these variables effectively separate the rotational and internal vibrational modes. The second variation of the effective Hamiltonian is block diagonal. separating the modes completely. while the symplectic form has an off diagonal term which represents the dynamic interaction between these modes. Otherwise, the symplectic form is in a type of normal form. The result sets the stage for the development of useful criteria for bifurcation as well as the stability criteria found here. In addition, the techniques should apply to other systems as well, such as rotating fluid masses

Normalizing connections and the energy-momentum method

The block diagonalization method for determining the stability of relative equilibria is discussed from the point of view of connections. We construct connections whose horizontal and vertical decompositions simultaneosly put the second variation of the augmented Hamiltonian and the symplectic structure into normal form. The cotangent bundle reduction theorem provides the setting in which the results are obtained