Home / Math and Physics Files / Physics / Mechanics / new frames doc / Unresolved Frames Doc Folders / Hard Questions
Question 2 about L=Iw
DOCX · 28.9 KB
Open DOCX file
Phil's working notes (dated 1.11.15, with comments added 3.12.17) from a frames document, posed as a question set on L = Iω. He argues the equation can be evaluated in either frame and that no separate L' = I'ω' exists. He then takes body-frame axes along the principal axes so I' is diagonal, and derives Euler's equations, compared with Goldstein and Marion. He ends by reviewing Goldstein Chapter 5.
AI-written summary; may contain errors.
Extracted text (machine-read; may contain errors)
About L = I ω PhL 1.11.15
Question Set. We have L = I ω .
Is this equation "in a frame"? Is this the Frame S equation and there is a Frame S' equation that reads L' = I' ω' like F = ma and F' = ma' ?
For the case that Frame S' is the body frame, it seems that L' = 0 for the rigid body or any particle of it, since the body is at rest in this body frame. Then L' = I' ω' would read 0 = I' ω'. What does that mean? In one sense ω' = 0 in Frame S' since nothing is moving, but in another sense, Frame S' is moving relative to Frame S at some ω' and perhaps ω' = -ω.
Comment added 3.12.17. In the body frame, a particle has v' = 0 and p' = 0 so p' = mv' says 0 = 0. This does not mean that v = 0 for that particle in Frame S.
If L = I ω is a Frame S equation with Frame S being "the table on which the top is spinning" then how can you say I is diagonal? I think of the principle axes being the body Frame S' axes, maybe that is wrong.
Maybe L = I ω is not an equation that lives within a frame? -
So let's take it on one piece at a time.
Question 2a: We have L = I ω . Is this equation "in a frame"? Is this the Frame S equation and there is a Frame S' equation that reads L' = I' ω' like F = ma and F' = ma' ?
Conjecture A: Holding off on the question for a bit, this L = I ω is a vector equation with a matrix involved and like any other vector equation with a matrix it can be evaluated either in Frame S or in Frame S' Here is how that would work
(L)i = Iij(ω)j Frame S
(L)'i = I'ij(ω')j Frame S'
Comment added 3.12.17 What is ω here? In my Fig 1, ω is ω for the Frame S' relative to Frame S. In that Figure I show a particle of mass m. This particle could be doing any motion you want. But if the particle is at rest relative to Frame S', then I claim below 4.1.1 that such a particle also rotates at ω about the same ω axis, It has a a different rotation radius, but the same dφ in dt. So if this is a particle of a rigid body at rest in Frame S', then that particle is doing ω in Frame S about the same ω axis that describes the Frame S' axes. In Frame S', however, this particle is at rest relative to the Frame S' axes. The particle is static, there is no axis fixed in Frame S' about which the particle has some ω' which is non-zero.
Question: How does this fit in with my claim: If Frame S' rotates at ω relative to Frame S, then Frame S rotates at -ω relative to Frame S' ? These ω's are what Observers static in each Frame would see. I guess you could in some sense call one of these two ω's by the name ω' then you have ω' = -ω. But this is different from talking about the ω of a particle as measured in two frames? I think this is a bad use of the symbol ω' because ω describes the relation between two frames, it is not a property of one frame or the other frame. In the example above, I would reserve ω to mean that of a rigid body particle at rest in Frame S and then ω'= 0 for that particle. Then ω is a property of a particle about some axis. I think no matter now a particle is moving, at some time t there is some unique axis and it has some unique ω. I show this in Fig 8.7.5. When Frame S' is the body-frame, ω is always refer to the Frame S ang velox of a particle in that rigid body.
I think I could write this in Dirac notation and it would all make sense with an operator I.
(L)i = <ei| L> = <ei| I ω> = <ei| I | ω> = <ei| I |ej><ej| ω> = Iij(ω)j
Like any operator we will find that
I' = RIR-1
and since I is not R, we will NOT find that I' = I.
There is no non-zero vector called L' and no vector called ω'. We take the Passive View because the notation L' already has a separate meaning! It is L' = mr' x v'. The symbols r' and v' also have the "special meaning" and are "used up" symbols, so we cannot "create" new vectors such as L' ≡ RL.
In terms of the natural v and v' we can define p ≡ mv and p' ≡ mv'. The law F = ma in Frame S' becomes F'eff = ma' , however. Now consider L = I ω . We already have a meaning for L' , so it might not work to try do define L' ≡ I' ω'. Also, I don't think ω' had any meaning. The rotation of Frame S' seen from Frame S is ω, and going the other way it is -ω. It is like b, there is no b'. So my tentative answer to Question 2a is that no, there is no equation L' = I' ω' which is similar to p' ≡ mv'. The latter is a definition and I cannot define L' ≡ I' ω' . [ not quite sure of this ]
Question 2b : For the case that Frame S' is the body frame, it seems that L' = 0 for the rigid body or any particle of it, since the body is at rest in this body frame. Then L' = I' ω' would read 0 = I' ω'. What does that mean? In one sense ω' = 0 in Frame S' since nothing is moving, but in another sense, Frame S' is moving relative to Frame S at some ω' and perhaps ω' = -ω. [ bad use of symbol ω' ]
First of all, I guess I accept that L' ≡ sum of r' x p' (in body frame Frame S') really is 0 since all the p' are zero. Recall that,
L(c)S ≡ (r-c) x mvS = (r-c) x mv = (r-c) x p ≡ L(c)
L'(c')S' ≡ (r'-c') x mv'S' = (r'-c') x mv' = (r'-c') x p' ≡ L'(c') . (1.9.3)
For common origins special case 4 we get
L ≡ r-c) x mv = r x p // for each rigid body particle
L' ≡ r' x mv' = r' x p'
Later in Section 11 I show that
L(c) = L'(c') + m(r'-c') x [ (ω x r') + S] (11.2.14)
so in Special Case 4 we get
L = L' + mr' x (ω x r')
Then since we just said above that for rigid body Frame S' we have L' = 0, we get
L = mr' x (ω x r')
I think this agrees with Marion page 365 where α just labels particles, and he has no primes. In my version r' is a location in the body of a particle, whereas L is the Frame S natural L.
Now let's roll out that old vector identity and do the double cross product ourselves right here.
L = - mr' x (r' x ω) A x (A x C) = (AC)A - A2C
= -m [ (r'ω)r' - r'2ω ] // L is the Frame S L !
= m [ r'2ω - (r'ω)r' ] // Marion page 365 (12-18)
Li = m [ r'2(ω)i - (r')j(ω)j(r')i ] // Frame S evaluation
= m [ r'2(ω)jδij - (r')j(ω)j(r')i ]
= m [ r'2δij - (r')j(r')i ] (ω)j
= Iij (ω)j Iij = m[ r'2δij - (r')j(r')i] but sum over all particles!
(L)'i = m [ r'2(ω)'i - (r')'j(ω)'j(r')'i ] // Frame S' evaluation
= m [ r'2(ω)'jδij - (r')'j(ω)'j(r')'i ]
= m [ r'2δij - (r')'j(r')'i ] (ω)'j
= I'ij (ω)'j I'ij ≡ m[ r'2δij - (r')'j(r')'i] but sum over all particles!
Now components of vector r' in Frame S' could be called x'i for short.
Now it is this thing I'ij that we would like to be diagonal! How could we arrange that to happen?
Question 2c : How could you arrange that the above I'ij be diagonal? If you arrange so that the axes of Frame S' line up with the principal axes of the rotating object, then I'ij will be diagonal. For an object with a symmetry axis, that will always be a principal axis. But object need not have a symmetry axis. For ANY object , you could arrange to have Frame S' axes line up with the principal axis of that object. The idea here is that you can always diagonalize a symmetric matrix, I could probably write that up to make it very explicit, but maybe don't need to do that here.
So I had a wrong impression of things I think. I was imagining that Iij was non-diagonal in Frame S and diagonal in Frame S'. Probably that is true IF you select Frame S' as I just said.
So now ASSUME that Frame S' has its axes lined up with principal axes. Then we get
I'ij = I'iδij
and we then get
(L)'i = I'ij (ω)'j = I'iδij (ω)'j = I'i(ω)'i
Note that this is the Frame S L on the left!
Question 2d : How do I get from here to an equation of motion? I do know that
= N = ∂SL in Frame S which is inertial.
I also know from the G Rule that
= (∂SL) = (∂S'L) + ω x L
where (∂S'L) is "unnatural". So we then have
N = (∂S'L) + ω x L
Take a component in Frame S' :
(N)'i = [(∂tL)S']'i + [ω x L]'i
Now use commutation theorem which says
[(∂tL)S']'i = ∂t(L)'i // think this is right
Then you get
(N)'i = ∂t(L)'i + [ω x L]'i
NOW use the result from above
(L)'i = I'i(ω)'i // no implied sum on i
and then we have
(N)'i = ∂t[ I'i(ω)'i] + [ω x L]'i
= I'i∂t[(ω)'i] + [ω x L]'i
I would like now to show that
∂t[(ω)'i] = (∂tω)'i
This is the commutation theorem again which says you can interchange if both same. Then we get
(∂tω)'i = ()'i
Meanwhile, in the other term we have
[ω x L]'i = εijk(ω)'j(L)'k = εijk(ω)'j[ I'k(ω)'k ]
Then our equations are
(N)'i = I'i ()'i + εijk(ω)'j(ω)'k I'k
And we can now just write these out
(N)'1 = I'1 ()'1 + ε123(ω)'2(ω)'3 I'3 + m ε132(ω)'3(ω)'2 I'2
= I'1 ()'1 + (ω)'2(ω)'3 I'3 - m (ω)'3(ω)'2 I'2
= I'1 ()'1 + (ω)'2(ω)'3 [ I'3 - I'2 ]
The equations are cyclic, so we write them all at once
(N)'1 = I'1 ()'1 + (ω)'2(ω)'3 [ I'3 - I'2 ]
(N)'2 = I'2 ()'2 + (ω)'3(ω)'1 [ I'2 - I'3 ]
(N)'3 = I'2 ()'3 + (ω)'1(ω)'2 [ I'1 - I'1 ]
This finally is Goldstein p 158 A. We know of course in terms of Euler angles that
(ω)'1 = sinθsinψ + cosψ
(ω)'2 = sinθcosψ - sinψ
(ω)'3 = cosθ + // Frame S' (G.6.11)
I guess that to get ()'i you do all the messy derivatives!
()'1 = sinθsinψ + cosθ sinψ + sinθ cosψ + cosψ - sinψ
and so on. But I don't think you install these at this point, you just work with the three equations as shown.
(N)'1 = I'1 ()'1 + (ω)'2(ω)'3 [ I'3 - I'2 ]
(N)'2 = I'2 ()'2 + (ω)'3(ω)'1 [ I'2 - I'3 ]
(N)'3 = I'2 ()'3 + (ω)'1(ω)'2 [ I'1 - I'1 ]
Everything here is in Frame S' components.
Comment: I think I have finally cleared up these nasty equations about L = I ω.
Question 2d: If L = I ω is a Frame S equation with Frame S being "the table on which the top is spinning" then how can you say I is diagonal? I think of the principle axes being the body Frame S' axes, maybe that is wrong.
Maybe L = I ω is not an equation that lives within a frame?
Answer: Well, the equation is L = I ω , and we evaluate it in Frame S' components and this gives
(L)'i = I'ij (ω)'j where L is indeed the Frame S L . (L' = 0! ) Then it is I'ij which is diagonal, and yes that is in Frame S'.
L = I ω is an equation which can be evaluated in either Frame. it does not live within a Frame.
Continue 3.12.17. OK, I guess the above is all fine and I can write it up. But how then to you solve a problem like a top spinning on a table? Let's do a little Goldstein review on this right now. It is Chapter 5 that is of interest here.
Question 2e: What is the meaning of saying In = n I n as on Goldstein p 155?
Well think of n = niei then maybe In = <n | I | n> = nT I n = niIijnj = < niei | I | njej>
= ninj< ei | I | ej> = ninj Iij. So mechanically this is what n I n means. But what is the physics meaning? "moment of inertia about the n axis" I could just define it this way. Below 5.15 Gold shows that T = (1/2) ω I ω = (1/2) ω2 I = (1/2) ω2 Iω . So then the physics meaning is that Iω is the thing that makes the formula T = (1/2) ω2 Iω work, analogous to T = (1/2) v2 m for linear motion. That seems OK to me. So really n = as he states under 5-15 page 149. This fact is derived on page 149.
On page 157 Goldstein writes the Lagrangian in Frame S' (but no primes displayed) in 5-31, where these ωi are the body-frame guys I derive in Appendix G. I agree with page 157 A, looking at those body frame formulas. Gold then writes the Euler-Lagrange equation for coordinate ψ and gets (5-33). I obtain this result above without using the Lagrange equations. Gold says "clearly we can permute the indices" to get the other equations. I suppose that would correspond to doing the other two Euler Lagrange equations.
Gold rederives 5-34 using my method above on page 158. It is the next sections that I now want to study up on from my G notes and from his book. Start with Section 5.6.