General Mechanics/Motion Under Constraint

From Wikibooks, open books for an open world
Jump to navigation Jump to search

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 | edit source]

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 | edit source]

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
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
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 | edit source]

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.