I'm currently studying Goldstein's Classical Mechanics book and I can't get my head around his reasoning in section 2.4. (Extending Hamilton's principle to systems with constraints). I'd like to understand the example he gives. Here it comes:
Consider a smooth solid hemisphere of radius a placed with its flat side down and fastened to the Earth whose gravitational acceleration is g. Place a small mass M at the top of the hemisphere with an infinitesimal displacement off center so the mass slides down without friction. Choose coordinate x,y,z centered on the base of the hemisphere with z vertical and the x-z-plane containing the initial motion of the mass.
Let θ be the angle from the top of the sphere to the mass. The Lagrangian is L=12⋅M⋅(˙x2+˙y2+˙z2)−m⋅g⋅z. The initial conditions allow us to ignore the y coordinate, so the constraint equation is a−√x2+y2=0. Expressing the problem in terms of r2=x2+z2 and x/z=cos(θ), Lagrange's equations are Ma˙θ2−Mgcos(θ)+λ=0 and Ma2¨θ+Mgasin(θ)=0
Solve the second equation and then the first to obtain ˙θ2=−2gacos(θ)+2ga and λ=Mg(3cos(θ)−2) So λ is the magnitude of the force keeping the particle on the sphere and since λ=0 when θ=cos−1(23), the mass leaves the sphere at that angle.
I have the following questions:
Shouldn't it be x/z=tanθ?
Could it be that he's mixing up r and a? My guess is that from "Lagrange's equations are" it should say r instead of a. I get confused whether a is a system parameter or a Lagrangian multiplier.
Could you give me a) an explanation or b) a good read on why setting L′=L+λ⋅f gives us an analogue of Hamilton's principle on constraint systems? I don't understand Goldstein's derivation. (L is the original Lagrangian, f is the constraint and λ is the Lagrangian multiplier.)
Why can λ be thought of as the constraint force?
When I understand 3., I understand the example -- I reverse engineered the supposedly Lagrangian equations to see that L′ needs to be of form 12Mr2˙θ2−Mrgcos(θ)+λ⋅f with generalized coordinates θ and r. Then everything works out just fine.
Answer
- Yes, it should be x/z=tanθ, this is probably a typo.
- The constraint should be a−√x2+z2=0 for the argument to make sense. r is a coordinate which is variable but due to the constraint it will always be equal to a, so we can use a in the equations instead. (˙a=0).
- You know that a gradient of f is always perpendicular to the surface of constant f, so you can understand the extra term coming from ∂L′/∂xi as a force acting perpendicular to the surface of f=0 holding the particle on it. However, λ has to be solved so that the motion of the system is only along the constant surface. But imagine now a force equal to the solved λ∂f/∂xi - it would have the same effect as the constraint, so this is the force by which the constraint actually has to be acting to hold the particle/system. (This also answers 4.)
No comments:
Post a Comment