Check my solution to this dynamics problem.

[mfig]

Hello,

I found this problem while preparing for the qualifying exams and thought I would tackle it. I have no idea if it is correct, but I would like someone to tell me if there is an easier way to solve this than I used. Perhaps some Lagrangian or Hamiltonian approach? Thanks.


The situation is as follows. A solid disk of radius R and mass M_{w} is resting on its edge on a horizontal plane surface near a rigid vertical pole. Free to slide frictionlessly on the pole is a collar of mass M_{c}. Connecting the collar to a frictionless pin at the center of the disk is a thin beam of length L and mass M_{b}. The angle between the beam and the plane is \theta_{0}. Assuming the system is released, find an expression for the velocity of the collar when the angle is θ.

To solve this I used conservation of mechanical energy.
The initial energy is all potential, and the energy at θ (once the system is released) is a mixture of kinetic and potential.

E_{{{\it before}}} = E_{{{\it pot.} {\it\space collar \space initial}}}+E_{{{\it pot.} {<br /> \it beam \space initial}}}

E_{{{\it pot.}{\it \space collar \space initial}}}=M_{{c}}\cdot g \cdot L\cdot \sin \left( <br /> \theta_{{0}} \right)

E_{{{\it pot.} {\it beam \space initial}}}=\frac{M_{{b}}\cdot g \cdot L\cdot \sin<br /> \left( \theta_{{0}} \right)}{2}<br />

Now we find the energies of the components when the system is released.

E_{{{\it kin.} {\it \space wheel}}}=\frac{{\it I}_{{{\it \space wheel}}}\cdot{\space <br /> {\omega}^{2}_{{{\it wheel}}}}}{2}+\frac{M_{{w}}\cdot{{V}^{2}_{{c. m. {\it \space wheel}<br /> }}}}{2}<br />

\omega_{{{\it wheel}}}={\frac {V_{{c. m. {\it \space wheel}}}}{R}}


{\it I}_{{{\it wheel}}}=\frac{M_{{w}}\cdot{R}^{2}}{2}


V_{{c. m. {\it \space wheel}}}={\frac {V_{{c}}}{\tan \left( \theta<br /> \right) }}<br />

E_{{{\it kin.} {\it \space collar}}}=\frac{M_{{c}}{V_{{c}}}^{2}}{2}

E_{{{\it pot.} {\it \space collar}}}=M_{{c}}\cdot g \cdot L\sin \left( \theta \right)

E_{{{\it kin.}{\it \space beam}}}=\frac{{\it I}_{{c. m. {\it <br /> \space beam}}}\cdot {{\omega}^{2}_{{{\it beam}}}}}{2} +\frac{M_{{b}}\cdot {{V}^{2}_{{c\cdot m\cdot {<br /> \it beam}}}}<br /> }{2}

E_{{{\it pot.} {\it \space beam}}}=\frac{M_{{b}}\cdot g \cdot L\cdot \sin \left( \theta<br /> \right)}{2}<br />

{\it I}_{{c. m. {\it \space beam}}}=\frac{M_{{b}}{L}^{2}}{12}

\omega_{{{\it beam}}}={\frac {2\cdot V_{{c}}\cdot \cos \left( \theta \right) }{L<br /> }}<br />

V_{{c. m.{\it \space beam}}}=\frac{V_{{c}}}{2}\cdot \sqrt {1+ \frac{1}{\left( {\tan}^{2}<br /> \left( \theta \right) \right) }}<br />

So putting this all together, we get an expression for the total energy of the system after release.

E_{{{\it after}}}=E_{{{\it kin.}{\it \space wheel}}}+E_{{{\it kin.} <br /> {\it \space collar}}}+E_{{{\it pot.} {\it \space collar}}}+E_{{{\it kin.} {<br /> \it \space beam}}}+E_{{{\it pot.}{\it \space beam}}}<br />


Equating this with the initial energy, and solving for V_{c} yields:

V_{c} = \sqrt { \frac{M_{{c}}\cdot g\cdot L\cdot \left( \sin \left( \beta \right) - \sin \left( \theta \right) \right)+\frac{M_{{b}}\cdot g\cdot L\cdot \left(\sin<br /> \left( \beta \right) - \sin<br /> \left( \theta \right) \right)}{2}}{{\frac {3\cdot M_{{w}}}{4\cdot<br /> {\tan} ^{2} \left( \theta \right)}}+\frac{M_{{c}}}{2}+\frac{M<br /> _{{b}}\cdot {\cos} ^{2} \left( \theta \right) }{6}+\frac{M_{{b}}<br /> \cdot \left( 1+ \frac{1}{{\tan}^{2} \left( \theta \right) } \right) <br /> }{8}}} <br />


For definiteness, I use:

M_{w} = 21
M_{c} = 33
M_{b} = 11
\theta_{0}=63°
\theta=32°
L = 3.11

And get V_{c} = 2.593

 

[haruspex]

I believe you have the ratio of velocities of wheel and collar inverted. Haven't checked any further.
(It wasn't obvious how the beam is attached to collar and wheel, but judging from the equations they seem to be a joints that move freely in a vertical plane; and I guess the pillar is in the plane of the disc.)

 

[mfig]

Thanks, haruspex.

I will double check that.

The beam is only there to connect the collar, which slides freely up and down the vertical pole, to the center of the wheel. The wheel stays on the ground and rolls as the collar drops down the pole. So the collar falls, pushing the wheel away from the pole. The distance between the collar and the center of the wheel is constant (the length of the beam). Does that clear it up?

Anyone else have a look? Thanks!



Here is an image I drew with Paint.


http://speedypic.net/di-7G0C.png

 

[mfig]

I added an image. Hopefully this helps clear up anything that wasn't clear from my description.

 

[mfig]

Hello,

I think haruspex is correct about the velocity ratios being inverted. The correct relation is:

V_{{c. m. {\it \space wheel}}}=V_{c}\cdot \tan \left( \theta \right)

So I fixed this and now get this solution. Note only one term was affected.V_{c} = \sqrt { \frac{M_{{c}}\cdot g\cdot L\cdot \left( \sin \left( \beta \right) - \sin \left( \theta \right) \right)+\frac{M_{{b}}\cdot g\cdot L\cdot \left(\sin \left( \beta \right) - \sin \left( \theta \right) \right)}{2}}{{\frac {3\cdot M_{{w}}\cdot {\tan} ^{2} \left( \theta \right)}{4}}+\frac{M_{{c}}}{2}+\frac{M _{{b}}\cdot {\cos} ^{2} \left( \theta \right) }{6}+\frac{M_{{b}} \cdot \left( 1+ \frac{1}{{\tan}^{2} \left( \theta \right) } \right) }{8}}}And using the same numbers as before (I realized I didn't list g=9.81), I get V_{c}=3.833

Can anyone think of another way to solve this that is simpler? Even a hint will help. Thanks.

 

[haruspex]

Looks to me that the inversion of the velocity ratio also affects the final term in your last line, making it collapse to M_b sec²(θ)/8.
Similarly, for ω_b, I get V_c sec(θ)/L: V_c = (d/dt)(L sin(θ)) = L cos(θ) ω_b
It can be confusing trying to work out these velocity relationships. Safest is to write down the distance relationship then differentiate.

 

[mfig]

Now I am confused, haruspex. I used this picture to find the expression for the angular velocity of the beam.
I was trying to find the component of the collar's velocity perpendicular to the beam, as shown, then use

\omega=\frac{V_{\bot}}{R}

R here is L/2 and V is the velocity of the collar. Why is this wrong?

 

[haruspex]

It's wrong because you've left out the fact that the disc's movement also contributes to the rotation of the beam.

 

[mfig]

Whoa. You are blowing my mind here. I assumed that the perpendicular components of the velocities of both ends of the beam would be the same. You are telling me that the perpendicular velocities of the two ends of the beam are not the same, even though the beam is rigid? And that we simply add each contribution to get ω?

How is it even possible to have two ends of a rigid beam turning it at different rates, perpendicular to the axis of rotation, without bending the beam?

Even then I still get a factor of 2 that you don't have.

Wheel contribution: \frac{V_{wheel} \cdot sin \left( \theta \right)}{\frac{L}{2}}

Collar contribution: \frac{V_{collar} \cdot cos \left( \theta \right)}{\frac{L}{2}}

Now the sum is:

\frac{V_{wheel} \cdot sin \left( \theta \right)}{\frac{L}{2}} + \frac{V_{collar} \cdot cos \left( \theta \right)}{\frac{L}{2}}

or:

\frac{V_{collar} \cdot tan \left( \theta \right) \cdot sin \left( \theta \right)}{\frac{L}{2}} + \frac{V_{collar} \cdot cos \left( \theta \right)}{\frac{L}{2}}

which simplifies to:

\frac{2 \cdot V_{collar}}{L\cdot cos \left( \theta \right)}

As you had. I have a factor of 2 that you don't. Where did I go wrong this time?

 

[haruspex]

I assumed that the perpendicular components of the velocities of both ends of the beam would be the same.

If that were the case then the instantaneous centre of rotation would be the beam's centre, but in general that will not be so.
Consider a plank standing upright against a wall. If you slide the bottom end away from the wall, the instantaneous centre of rotation will at first be the top end of the plank. Nearly all the linear movement at this time is at the bottom end, and the rotation rate will be that movement speed divided by the plank length.
When the plank has dropped to 45 degrees, the two ends are moving at the same linear rates, v say. The instantaneous centre of rotation is now the centre of the plank.
More generally, at angle θ to the horizontal, the height of the top end of the plank is L sin(θ), and the distance of the bottom from the wall is L cos(θ). If the rotation rate (dθ/dt) is ω then the two speeds are L cos(θ) ω and L sin(θ) ω respectively (chain rule).

we simply add each contribution to get ω?

Not exactly. As you can see from the plank/wall discussion, it's more complicated than that. With the two endpoint movements being at right-angles to each other, it's more like a root-sum-square (divided by plank length).
If the two endpoint movements were parallel then it would be the difference that led to the rotation (i.e. adding magnitudes if in opposite directions). In the trivial case of the two endpoints moving with identical velocity vectors there is of course no rotation.

 

[mfig]

What you say makes good sense. I sincerely appreciate your patience with me on this. I want to re-learn this stuff but it has been a long time.

So the instantaneous point of rotation will not be the center. That means we cannot just add the contributions. But what method can we use to get to ω? You say it is " more like a root-sum-square" but what does that mean? Is it a root-sum-square or only like one? You got a definite result, so please show me how you arrived there. Did you explicitly find an expression for the instantaneous center of rotation?

Thanks again!

 

[haruspex]

You could take one endpoint as the frame of reference. The relative movement of the other end is just the difference of the two movements as vectors. Since the beam is constant in length, this is necessarily perpendicular to the beam, so you can then just take its magnitude and divide by the beam length. You can see this trivially gives the right answer when the movements are parallel. In the present case, the movements being perpendicular to each other, the resultant has root-sum-square magnitude.