**Next:** Discretization **Up:** Mathematical Formulation **Previous:** Mathematical Formulation

## Constraint Dynamics

Numerous physical systems lead to a Hamiltonian system of ordinary differential equations:

where and are the momenta and positions of the particles, respectively. However, such a set of equations is altered when the motion of the bodies (in our case, unit-mass particles) is restricted by a constraint of the form

New equations of motion reflect the addition of a second ``force'' from the constraint acting upon the particles. For the general constrained system, determining the trajectories of motion now involves solving constrained dynamical equations of the form:

Here, is the Jacobian of the vector of constraints and is the vector of Lagrange multipliers associated with the constraints.
Solutions are said to lie on the solution manifold, , if the constraint is satisfied. If one differentiates the constraint with respect to time, another, ``hidden'' constraint is obtained. When the constraint represents a surface, the hidden constraint confines particle velocities to a plane tangent to the surface, providing a fourth and final equation of motion:

The position and momentum pairs which satisfy the two constraints are said to lie in the tangent bundle. The exact solution of the constrained system consists of position and momentum pairs that are members of this tangent bundle.
We can model the lattice (and traditional) Newton's Cradle using this formulation of constrained differential equations. In this model, the constraints describe the surfaces of spheres in three-dimensional space:

where and are the positions of the th ball and hook respectively. For the lattice Newton's Cradle, and are vectors in , and is the length of the pendulum rods. Introducing this constraint equation into equations (4-6), our system is reduced to:

**Next:** Discretization **Up:** Mathematical Formulation **Previous:** Mathematical Formulation

*1999-05-24*