Chapter 7

**Hamilton's
Principle - Lagrangian and Hamiltonian Dynamics**

Many interesting physics systems describe systems of particles on which many forces are acting. Some of these forces are immediately obvious to the person studying the system since they are externally applied. Other forces are not immediately obvious, and are applied by the external constraints imposed on the system. These forces are often difficult to quantify, but the effect of these forces is easy to describe. Trying to describe such a system in terms of Newton's equations of motion is often difficult since it requires us to specify the total force. In this Chapter we will see that describing such a system by applying Hamilton's principle will allow us to determine the equation of motion for system for which we would not be able to derive these equations easily on the basis of Newton's laws. We should stress however, that Hamilton's principle does not provide us with a new physical theory, but it allows us to describe the existing theories in a new and elegant framework.

__Hamilton's
Principle__

The evolution of many physical systems involves the minimization of certain physical quantities. We already have encountered an example of such a system, namely the case of refraction where light will propagate in such a way that the total time of flight is minimized. This same principle can be used to explain the law of reflection: the angle of incidence is equal to the angle of reflection.

The minimization approach to physics was formalized in detail by Hamilton, and resulted in Hamilton's Principle which states:

"* Of all the possible paths
along which a dynamical system may more from one point to another within a
specified time interval (consistent with any constraints), the actual path
followed is that which minimizes the time integral of the difference between
the kinetic and potential energies.* "

We can express this principle in terms of the calculus of variations:

The quantity *T* - *U* is called the Lagrangian *L*.

Consider
first a single particle, moving in a conservative force field. For such a particle, the kinetic energy
*T* will just be a function of the
velocity of the particle, and the potential energy will just be a function of
the position of the particle. The
Lagrangian is thus also a function of the position and the velocity of the
particle. Hamilton's theorem
states that we need to minimize the Lagrangian and thus require that

In Chapter 6 we have developed the theory required to solve problems of this type and found that the Lagrangian must satisfy the following relation:

These last equations are called the **Lagrange equations of
motion. **Note that in order to generate these equations of
motion, we do not need to know the forces. Information about the forces is included in the details of
the kinetic and potential energy of the system.

Consider
the example of a plane pendulum.
For this system, there is only one coordinate we need to specify, namely
the polar angle *q*. The kinetic energy *T* of the pendulum is equal to

and the potential energy *U* is given by

The Lagrangian for this system is thus equal to

The equation of motion can now be determined and is found to be equal to

or

This equation is of course the same equation we can find by
applying Newton's force laws. I
this example, the only coordinate that was used was the polar angle *q*. Even though the pendulum is a 3-dimensional system, the
constraints imposed upon its motion reduced the number of degrees of freedom
from 3 to 1.

__Generalized
Coordinates__

If
we try to describe a system of *n*
particles, we need in general 3*n*
coordinates to specify the position of its components. If external constraints are imposed on
the system, the number of degrees of freedom may be less. If there are *m* constraints applied, the number of degrees of
freedom will be 3*n* - *m*. The
coordinates do not need to be the coordinates of a coordinate system, but can
be any set of quantities that completely specifies the state of the
system. The state of the system is
thus full specified by a point in the configuration space (which is a 3*n* - *m*
dimensional space). The time
evolution of the system can be described by a path in the configuration space.

The
generalized coordinates of a system are written as *q*_{1}, *q*_{2}, *q*_{3},
ÉÉ The generalized coordinates are
of course related to the physical coordinates of the particles being described:

where *i* = 1, 2, 3 and
*a* = 1, 2, É.., *n*. Since
the generalized coordinates in general will depend on time, we can also
introduce the generalized velocities.
The physical velocities will depends on the generalized velocities:

__Equations
of Motion in Generalized Coordinates__

Based on the introduction of the Lagrangian and generalized coordinates, we can rephrase Hamilton's principle in the following way:

"* Of all the possible paths
along which a dynamical system may more from one point to another in
configuration space within a specified time interval (consistent with any
constraints), the actual path followed is that which minimizes the time
integral of the Lagrangian function for the system.* "

Thus

and

When we use the Lagrange's equations to describe the evolution of a system, we must recognize that these equations are only correct of the following conditions are met:

1. the force acting on the system, except the forces of constraint, must be derivable from one or more potentials.

2. the equations of constraint must be relations that connect the coordinates of the particles, and may be time dependent (note: this means that they are independent of velocity).

Constraints that do not depend on velocity are called **holonomic
constraints**. There are two different types of holonomic constraints:

1.
fixed or **scleronomic constraints**: constraints that do not depend on time.

2.
moving or **rheonomic constraints**: constraints that depend on time.

__Example:
Problem 7.4__

A
particle moves in a plane under the influence of a force *f* = -*Ar*^{a}^{-1}
directed toward the origin; *A* and *a* are constants. Choose appropriate generalized
coordinates, and let the potential energy be zero at the origin. Find the Lagrangian equations of
motion. Is the angular momentum
about the origin conserved? Is the
total energy conserved?

Figure 1. Problem 7.4.

If we choose (*r*,*q*)
as the generalized coordinates, the kinetic energy of the particle is

(7.4.1)

Since the force is related to the potential by

(7.4.2)

we find

(7.4.3)

where we let *U*(*r* = 0) = 0.
Therefore, the Lagrangian becomes

(7.4.4)

LagrangeÕs equation for the coordinate *r* leads to

(7.4.5)

LagrangeÕs equation for the coordinate *q*
leads to

(7.4.6)

Since is identified as the angular momentum, (7.4.6) implies that angular momentum is conserved. Now, if we use l, we can write (7.4.5) as

(7.4.7)

Multiplying (7.4.7) by , we have

(7.4.8)

which is equivalent to

(7.4.9)

Therefore,

(7.4.10)

and the total energy is conserved.

__Example:
Problem 7.8__

Consider
a region of space divided by a plane.
The potential energy of a particle in region 1 is *U*_{1} and in region 2 it is *U*_{2}.
If a particle of mass *m*
and with speed *v*_{1} in
region 1 passes from region 1 to region 2 such that its path in region 1 makes
and angle *q*_{1}
with the normal to the plane of separation and an angle *q*_{2} with the normal when in region 2, show that

where *T*_{1} =
(1/2)*mv*_{1}^{2},

Figure 2. Problem 7.8.

Let us choose the (*x*,
*y*) coordinates so that the two
regions are divided by the *y*
axis:

If we consider the potential energy as a function of *x* as above, the Lagrangian of the particle is

(7.8.1)

Therefore, LagrangeÕs equations for the coordinates *x* and *y*
are

(7.8.2)

(7.8.3)

Using the relation

(7.8.4)

(7.8.2) becomes

(7.8.5)

Integrating (7.8.5) from any point in the region 1 to any point in the region 2, we find

(7.8.6)

(7.8.7)

or, equivalently,

(7.8.8)

Now, from (7.8.3) we have

and is constant. Therefore,

(7.8.9)

From (7.8.9) we have

(7.8.10)

Adding (7.8.8) and (7.8.10), we have

(7.8.11)

From (7.8.9) we also have

(7.8.12)

Substituting (7.8.11) into (7.8.12), we find

(7.8.13)

This problem is the mechanical analog of the refraction of light upon passing from a medium of a certain optical density into a medium with a different optical density.

__Lagrange's
Equations with Undetermined Multipliers__

We have seen already a number of examples were one could remove the equations of constraint by a suitable choice of coordinates. For example, when we looked at the motion of an object on the surface of a cylinder we could either:

1. Use a set of three coordinates to describe the motion, coupled with one equation of constraint.

2. Use a set of two coordinates (such at the azimuthal angle and the vertical position) to describe the motion, without an equation of constraint.

In this Section we will look at situations where the constraint depends on the velocity:

If the constraints can be expressed in a differential form,

we can directly incorporate it into the Lagrange equations:

It turns out that the forces of constraint can be determined
from the equations of constraints and the Lagrange multipliers *l*_{m}(*t*):

where *Q*_{j}
is the j^{th} component of the generalized force,
expressed in generalized coordinates.

The
use of Lagrange multiplier to determine the forces of constraints is nicely
illustrated in Example 7.9 in the textbook, where a disk rolling down an
inclined plane is being studied.
If the disk does not slip, we find that the distance along the plane *y* and the angle of rotation *q* are related, and the equation of
constraint is

The textbook explains in detail how the Lagrange equations are solved in this case, and I will not reproduce this here. The solution shows us that the Lagrange multiplier is given by

By combining the equation of constraint and the Lagrange multipliers we can determine the generalized forces of constraint:

and

Note
that these forces of constraint do not have to be all pure forces. The force of constraint associated with
*y* is the friction force between the disk
and the plane that is required to ensure that the disk rolls without
slipping. However, the force of
constraint associated with the angle *q*
is the torque of this friction force with respect to the center of the
disk. We need to note the
generalized force does not have to have the unit of force.

It
is also important to note that if we had chosen to solve the problem by
expressing the Lagrangian in terms of a single coordinate *y*, by eliminating the angle, we would not have
obtained any information about the forces of constraint. Although I have stressed that in many
cases, you can simplify the solution by the proper choice of coordinates such
that the equations of constraint are eliminated, in this case, the solution
will not provide any information about the forces of constraint.

__Example:
Problem 7.12__

A
particle of mass *m* rests on a smooth
plane. The plane is raised to an
inclination angle *q* at a
constant rate *a* (*q* = 0¡ at *t* = 0), causing the particle to move down the
plane. Determine the motion of the
particle.

Figure 3. Problem 7.12.

This problem is an example of a problem with a velocity-dependent constraint. However, if we can easily incorporate the constraint into the Lagrangian, we do not need to worry about constraint functions. In this example, we use our knowledge of the constraint immediately in our expression of the kinetic and the potential energy. Putting the origin of our coordinate system at the bottom of the plane we find

LagrangeÕs equation for *r* gives

or

(7.12.1)

The general solution is of the form *r* = *r*_{h}
+ *r*_{n} where *r*_{h} is the general solution of the
homogeneous equation and *r*_{n} is a particular solution of Eq.
(7.12.1). So

For *r*_{n},
try a solution of the form *r*_{p}
= *C* sin*a**t*. Then . Substituting
into (7.12.1) gives

So

We can determine *A*
and *B* from the initial
conditions:

(7.12.2)

(7.12.3)

Equation (7.12.2) implies:

Equation (7.12.3) implies:

Solving for *A* and *B* gives:

or

__ __

Although we have found an analytical solution to this
problem, we need to examine if the solution matches our expectation of the
motion of the mass *m*. The best way to do this is to plot a
graph of the motion of the mass in a Cartesian coordinate system. Consider the situation where *r*_{0} = 10 m. Figure 4 shows the trajectory of
the mass for two different value of *a*:
*a* = 0.1 rad/s and *a* = 0.03 rad/s.

__ __

Figure 4.
Solution of Problem 7.12 with *r*_{0}
= 10 m and *a* = 0.01 rad/s
(left) and *a* = 0.03 rad/s
(right).

__ __

__ __

__ __

__The
generalized momentum__

One of the big differences between the equations of motion obtained from the Lagrange equations and those obtained from Newton's equations is that in the latter case, the coordinate frame used is always a Cartesian coordinate frame. When we use the Lagrange equations we have the option to choice generalized coordinates that do not have to correspond to the coordinates of a Cartesian coordinate system.

The generalized coordinates are related to the Cartesian coordinates, and transformation rules allow use to carry out transformations between coordinate systems. The generalized forces of constraint are related to the Newtonian forces of constraint, as was illustrated in Example 7.9 in the textbook. The similarities between the Cartesian and the generalized parameters suggest it may also be useful to consider the concept of the generalized momentum.

In a Cartesian coordinate system we can easily determine the connection between the Lagrangian and the linear momentum. The Lagrangian is equal to

The Lagrange equation for this Lagrangian is given by

and we can rewrite this as

This last equation suggest that we define the generalized momentum of a particle in the following way:

It is obviously consistent with our definition of linear momentum in Cartesian coordinates.

Consider a particle moving in a two-dimensional plane and having its motion described in terms of spherical coordinates. The kinetic energy of the particle is equal to

Since there are two generalized coordinates we can determine two generalized momenta:

which is the linear momentum of the particle, and

which is the angular momentum of the particle. We thus see that two distinct concepts from our introductory courses emerge directly from our Lagrangian theory.

__Homogeneous
functions__

Consider
a homogeneous quadratic function *f* that
depends only on the products of the generalized velocities:

An example of such function would be the kinetic energy of a particle. Consider what happens when we differentiate this function with respect to one of the generalized velocities:

If we multiply this equation by *dq*_{l}/*dt* and sum over all values of l we obtain:

In general, if *f* is a
homogeneous function of the parameter *y*_{k}^{n}, then

__Conservation
of Energy__

If we consider a closed system, a system that does not interact with its environment, then we expect that the Lagrangian that describes this system does not depend explicitly on time. That is

Of course, this does not mean that *dL*/*dt* = 0
since

Using the Lagrange equations, we can rewrite this equation as

This equation can be written as

or

The constant *H* is called the Hamiltonian of the system and the
Hamiltonian is defined as

The Hamiltonian *H* is a conserved quantity for the system we are
currently considering. If we use
Cartesian coordinates we find that

In this case we find that the
Hamiltonian of the system is equal to the total energy of the system, and we
thus conclude that the total energy is conserved. The equality of *H* and
*E* is only satisfied if the
following conditions are met:

á The
potential *U* depends only on position,
and not on velocity.

á The transformation rules connecting Cartesian and generalized coordinates are independent of time.

The latter condition is not met in for example a moving coordinate system, and in such system, the Hamiltonian will not be equal to the total energy.

We thus conclude that if the Lagrangian of a system does not depend explicitly on time, the total energy of that system will be conserved.

__Example
Problem 7.22__

A
particle of mass *m* moves in one
dimension under the influence of a force *F*:

where *k* and *t* are positive constants. Compute the Lagrangian and Hamiltonian
functions. Compare the Hamiltonian
and the total energy and discuss the conservation of energy for the system.

The
potential energy *U* corresponding to this
force *F* must satisfy the relation

and *U* must thus be
equal to

Therefore, the Lagrangian is

The Hamiltonian is given by

so that

The Hamiltonian is equal to the total energy, *T* + *U*,
because the potential does not depend on velocity, but the total energy of the
system is not conserved because *H*
contains the time explicitly.

__Conservation
of Linear Momentum__

The Lagrangian should be unaffected by a translation of the entire system in space, assuming that space is homogeneous (which is one of the requirement of an inertial reference frame). Consider what happens when we carry out an infinitesimal displacement of the coordinate system along one of the coordinate axes. The change in the Lagrangian as a result of this displacement must be equal to zero:

We can rewrite this equation as

Since the displacement is arbitrary, this equation can only be correct if

Using the Lagrange equation this is equivalent to requiring

or

Assuming that the potential *U* does not depend on velocity we see that this
relation is equivalent to

The consequence of the independence of the Lagrangian under a translation of space is that linear momentum is conserved.

__Conservation
of Angular Momentum__

Space is an inertial reference frame is isotropic, which means that the properties of a system are unaffected by the orientation of the system. In this case we expect that the Lagrangian does not change when the coordinate axes are rotated through an infinitesimal angle. A rotation through such an angle produces the following change in the position vector:

The velocity vectors will change in the same way:

The Lagrangian should not change as a result of such a transformation. Thus we must require that

We thus conclude that

When we express the changes in terms of the rotation angle we obtain:

Since the angle of rotation was an arbitrary angle, this relation can only be satisfied if

or

The angular momentum of the system is thus conserved. This conserved quantity is a direct consequence of the invariance of the Lagrangian for infinitesimal rotations. We conclude that the important conserved quantities are a direct consequence of the properties of the space (and its symmetries).

__Canonical
Equations of Motion__

The Lagrangian we have discussed in this Chapter is a function of the generalized position and the generalized velocity. The equations of motion can also be expressed in terms of the generalized position and the generalized momentum. The generalized momentum is defined as

We can use the generalized momentum to rewrite the Lagrange equations of motion:

The Hamiltonian can also be expressed in terms of the generalized momentum

In general we will write the
Hamiltonian in terms of the generalized position and the generalized momentum
. The change in *H* due to small changes in time and in the generalized
position and momentum is equal to

The change in *H* can also be expressed in the following way:

After combining the last two equations we obtained the following relation:

or

Since the variations in time and the
generalized position and momenta are equal to independent, the coefficients of *dq*_{i}, *dp*_{i}, and *dt* must
be zero. Thus:

The first two equations are called **Hamilton's
equations of motion** or the **canonical
equations of motion**. Note:

á For each generalized coordinate there are two canonical equations of motion.

á For each generalized coordinate these is only one Lagrange equations of motion.

á Each canonical equation of motion is a first order differential equation.

á Each Lagrange equation of motion is a second order differential equation.

Although first order differential equations are in general easier to solve than second order differential equations, the Hamiltonian is often more difficult to construct than the Lagrangian since we must express the Hamiltonian in terms of the generalized position and the generalized momentum.

__Example:
Problem 7.38__

The
potential for an anharmonic oscillator is *U*
= *kx*^{2}/2 + *bx*^{4}/4 where *k* and *b*
are constants. Find Hamilton's equations
of motion.

The Hamiltonian of the system is

The Hamiltonian motion equations that follow this Hamiltonian are

__Example:
Problem 7.28.__

The
force *F* that is provided fixed the potential *U*:

The Lagrangian, expressed in polar coordinates, is thus equal to

In order to use Hamilton's equations of motion we must express the Hamiltonian in terms of the generalized position and momentum. The following relations can be used to do this:

Since the coordinate transformations are independent of *t*, and the potential energy is velocity-independent,
the Hamiltonian is the total energy.

HamiltonÕs equations of motion can now be found easily

Consider
a simple plane pendulum consisting of a mass *m* attached to a string of length *l*. After the pendulum is
set into motion, the length of the string is shortened at a constant rate:

The suspension point remains fixed. Compute the Lagrangian and Hamiltonian functions. Compare the Hamiltonian and the total energy, and discuss the conservation of energy for the system.

The kinetic energy and the potential energy of the system are expressed as

(7.24.1)

The Lagrangian is equal to

(7.24.2)

The Hamiltonian is

(7.24.3)

which is different from the total energy, *T* + *U*. The total energy is thus not conserved
in this system because work is done on the system and we have

(7.24.4)

**NOTE: WE WILL SKIP SECTIONS 7.12 AND 7.13.**

** **