The Student Room Group

Dynamics

Ok, so I can do q1 and 3 this week but not 2 :P

A smooth straight wire rotates with constant angular speed &#969; about the vertical axis through a fixed point O on the wire, and the angle between the wire and the upward vertical is constant and equal to &#952;, where 0<&#952;<½pi. A bead is free to slide on the wire. Show that z(t), the height of the bead above O, satisfies the equation

z:-(&#969;² sin² &#952; ) z=-gcos² &#952;

where z:=d²z/dt²

and find z(t) if z(0)=h and z.(0)=0. Deduce that if h>gcot²&#952;/&#969;², then
z(t)->&#8734; as t->&#8734; but if h<gcot²&#952;/&#969;², then z(t)->-&#8734; as t->&#8734;.

This is a big question it really is but I cant do any of it. I would be really grateful for help :smile:....even if you can only do one bit of it :redface:

Cheers :tsr:
Reply 1
lgs98jonee
Ok, so I can do q1 and 3 this week but not 2 :P

A smooth straight wire rotates with constant angular speed ? about the vertical axis through a fixed point O on the wire, and the angle between the wire and the upward vertical is constant and equal to ?, where 0<?<½pi. A bead is free to slide on the wire. Show that z(t), the height of the bead above O, satisfies the equation

z:-(?²sin²?)z=-gcos²?

where z:=d²z/dt²

and find z(t) if z(0)=h and z.(0)=0. Deduce that if h>gcot²?/?², then
z(t)->? as t->? but if h<gcot²?/?², then z(t)->-? as t->?.

This is a big question it really is but I cant do any of it. I would be really grateful for help :smile:....even if you can only do one bit of it :redface:

Cheers :tsr:


Smilie in wrong place
Reply 2
steve2005
Smilie in wrong place


sorted :biggrin:
Reply 3
Surely someone can help me :frown:
Reply 4
I've done the first bit (havent looked at the later bits yet):

The tricky bit was coming up force on the bead that would make it fly up the wire. The point of wire that the bead is sat on is accelerating towards the axis of rotation at r&#969;² m/s/s. Resolving this along the line of the wire gives r&#969;²sin&#952;. The bead isnt fixed to the wire so it is as if the wire is sliding out from under it (providing the bead is light enough) with the same acceleration. Part of the acceleration of the bead up the wire is therefore r&#969;²sin&#952;. The acceleration due to gravity (down the wire) resolves to gcos&#952;. The resultant acceleration up the wire is therefore given by

r&#969;²sin&#952;-gcos&#952;

This is an acceleration in the direction of the wire, and we want z:, the acceleration upwards, so we make a small modification by resolving into a vertical component (we're not interested in the horizontal):

z: = cos&#952;(r&#969;²sin&#952;-gcos&#952; )

All that remains is to replace that r, and simple trigonometry says that r=ztan&#952; (or r=zsin&#952;/cos&#952; ). Plugging that in:

z: = cos&#952;((zsin&#952;/cos&#952; )&#969;²sin&#952;-gcos&#952; )

and multiplying out:

z: = z&#969;²sin²&#952;-gcos²&#952;

which rearranges to the equation you have:

z: - z&#969;²sin²&#952; = -gcos²&#952;

I'll have a go at the other parts in a moment, my pizza just arrived. Sorry for the wordyness, that's just how I work when I do a mechanics question.
Reply 5
Is it something along the lines of....

if h>gcot²&#952;, then z:>0 so accelerating upwards
and if h<gcot²&#952;, then z:<0, so accelertating downwards

??
Reply 6
A difficult question (and quite refreshing, I'm on a gap year and trying to get back into this for upcoming interviews). Finding z(t) involves solving the second order differential equation you end up with in part one. The general solution for an equation of this type (as you probably know) is:

z(t)=Ae&#955;1 t+Be&#955;2 t+C
I'll call this (1)

Where &#955;1 and &#955;2 are the solutions of the auxilliary equation:

&#955;²-&#969;²sin²&#952;=0
(so &#955;1=&#969;sin&#952; and &#955;2=-&#969;sin&#952; )

The constants take a bit of work to determine but it's all fairly straightforward algebra: Differentiate (1) to get z.(t), and use z.(0)=0 to arrive at A=B. Then use z(0)=h in the (1) to get 2A+C=h. Thirdly, use the equation arrived at in the first part, put z(0)=h into it and conclude that z: (0)=h&#969;²sin²&#952;-gcos²&#952;. Differentiate z.(t) to get another form z: (t) and put z: (0)=h&#969;²sin²&#952;-gcos²&#952; to hopefully get your three constants A, B and C. Let me know if you need any of this expanding on.

You should arrive at A=B=½h-(gcos²&#952;/2&#969;²sin²&#952; ) or ½h-(gcot²&#952;/2&#969;²) and C=gcos²&#952;/&#969;²sin²&#952;

So your final z(t) is defined as

z=(½h-(gcot²&#952;/2&#969;²))(e&#969;tsin&#952;+e-&#969;tsin&#952;)+gcos²&#952;/&#969;²sin²&#952;
I'll call this (2)

I've checked that and I think it all works out, but I am known to slip up occasionally. I'll know for sure when I've done the rest of it, bear with me.

Edit: In fact, I think I have made a mistake. z.(t) needs to be a function involving h also and according to the above, it isn't. I'll try and find my error.
Edit2: I put z(0)=0 instead of z(0)=h, it's right now I think.
Edit3: Made another stupid error (left out a squared sign), I've corrected this post and now the very last part is fairly trivial:

Looking at (2) in pieces:

(e&#969;tsin&#952;+e-&#969;tsin&#952;)
This part gets more positive with increasing time, since the first half just increases exponentially and the other half gets increasingly tiny (and remains positive).

+gcos²&#952;/&#969;²sin²&#952;
This is a constant, negative because g is negative and the rest is squared. Although it is negative it becomes less significant as t increases since the magnitude of the first half of (2) increases rapidly.

(½h-(gcot²&#952;/2&#969;²)) or ½(h-(gcot²&#952;/&#969;²))
This is the crux of the matter, the sign of z as t becomes large rests here because it is constant and it multiplies the significant (and positive) part of (2). Obviously, if h>gcot²&#952;/&#969;² then this constant is positive and the whole thing will get increasingly positive (z tends to &#8734; ), similarly if h<gcot²&#952;/&#969;² then z tends to -&#8734;.