General Mechanics/Motion Under Constraint

So far we've tacitly assumed we can just calculate force as a function of position, set up the ODE

and start solving.

It is not always so simple.

Often, we have to deal with motion under constraint; a bead sliding on a wire, a ball rolling without slipping, a weight dangling from a string.

There has to be some force keeping the bead on the wire, but we don't know what it is in advance, only what it does. This isn't enough information for us to write down the ODE.

We need a way of solving the problem without knowing the forces in advance.

How easy this is depends on the type of constraint.

  • If the constraint is an inequality, as with the weight on the string, there is no straightforward analytical method.
  • If the constraint can be written as a set of differential equations, and those equations can't be integrated in advance, there is an analytical method, but it is beyond the scope of this book. A ball rolling without slipping falls into this category.
  • If the constraint can be written as a set of algebraic equations, and frictional forces are negligible, there is a straightforward method that solves the problem.

Generalised coordinates edit

Suppose we have a system of n particles satisfying k constraints of the form

 

then we can use the constraints to eliminate k of the 3n coordinates of the particles, giving us a new set of 3n-k independent generalised coordinates; q1, q2, … q3n-'k.

Unlike the components of the position vectors, these new coordinates will not all be lengths, and will not typically form vectors. They may often be angles.

We now need to work out what Newton's laws will look like in the generalised coordinates.

Derivation edit

The first step is to eliminate the forces of constraint.

We will need to consider a virtual displacement. This is an infinitesimal displacement made, while holding the forces and constraints constant. It is not the same as the infinitesimal displacement made during an infinitesimal time, since the forces and constraints may change during that time.

We write the total force on particle i as
 
the sum of the externally applied forces and the forces of constraint.
Newton's second law states
 
We take the dot product of this with the virtual displacement of particle i and sum over all particles.
 
We now assume that the forces of constraint are perpendicular to the virtual displacement. This assumption is generally true in the absence of friction; e.g, the force of constraint that keeps a ball on a surface is normal to the surface.
This assumption is called D'Alembert's principle. Using it we can eliminate the forces of constraint from the problem, giving
 
or
 
The left hand side of this equation is called the virtual work.

Now we must change to the generalised co-ordinate system.

We write
 
Using the chain rule gives
 
and
 
Note that this implies
 
The virtual work is, dropping the superscript,
 
where the Qj are the components of the generalised force.

We now manipulate the right hand side of (1) into a form comparable with this last equation

The right hand side term is
 
The terms in the coefficient of qj can be rearranged
 
on substituting in equation (2) from above

On taking a close look at this last equation, we see a resemblance to the total kinetic energy,

 
We now further rearrange to get an expression explicitly involving T.
 

Putting this last expression into (1) along with the generalised force give

 

Since the δqj, unlike the δri, are independent, this last equation can only be true if all the coefficients vanish.

That is we must have

 

These are the equations of motion for the system, in a general set of coordinates for which all constraints are automatically satisfied.

For example, suppose we have a cylinder, mass m, radius a, rolling without slipping on a flat plane.

The kinetic energy of the cylinder is

 

using the results from Rigid Bodies, where x is the axis in the plane perpendicular to the axis of the cylinder, and θ is the angle of rotation.

Rolling without slipping implies

 

so we get

 

The cylinder has the same kinetic energy as if its mass were 20% greater. If there is no torque on the cylinder then Qx=Fx, and the cylinder behaves in every respect as though it were a 20% larger point mass.

To use (3) more generally, we need an expression for the Qj

Suppose, as is often the case, that

 

then, by definition

 

so, equating coefficients, the generalised force is

 

Putting this generalised force into (3) gives

 

since V has been assumed independent of the velocities.

In fact, this last equation will still be true for some velocity dependent forces, most notably magnetism, for a suitable definition of V, but we won't prove this here.

We call the T-V the Lagrangian, L, and write

 

We call these equations Lagrange's equations. They are useful whenever Cartesian co-ordinates are inconvenient, including motion under constraint.

Example edit

Suppose we have two identical point masses, m, connected by a string, length a. The string is threaded through an hole in a flat table so that the upper mass is moving in a horizontal plane without friction, and the lower mass is always vertically below the hole. The distance of the upper mass from the hole is r.


The position of the mass on the table is best described using polar coordinates, (r,θ). Its kinetic energy is then

 

The velocity of the lower mass is d(a-r)/dt=-dr/dt, so the total kinetic energy is

 

The potential energy is

 

where g is the gravitational acceleration.

This means

 

and the equations of motion are

 

The first of these equations says that the angular momentum is constant, as expected since there is no torque on the particles. If we call this constant angular momentum l then we can write

 

and the second equation of motion becomes

 

Clearly, if initially

 

then the lower ball will be pulled out of the hole, at which point these equations of motions cease to apply. They only hold when 0≤ra, a constraint which is not easily tamable.

Notice, we have not needed to calculate the tension in the string, which is the force of constraint in this problem.