Consistency of the Cartesian G Rules
DOCX · 47.1 KB
Open DOCX file
Working note by Phil dated 7.14.12, with later bracketed annotations. It computes velocity and acceleration of an ant on a rotating record in Cartesian components, first by differentiating v and then by the G rule a = ω x v + (dv/dt)S'. The two results seem to disagree, a paradox that the annotations resolve: the methods give different accelerations, aS versus the cross acceleration aS'.
AI-written summary; may contain errors. This description is approximate.
Extracted text (machine-read; may contain errors)
Consistency of the Cartesian G Rules PhL 7.14.12
Here I want to show that, if you stay in Cartesian coordinates, the four "G rules" for r, v, a and L are all consistent. I will do this in the "ant on phono record" context.
[ This is the first document (by creation date) that shows the (dr/dt)S notation, July 14. ] [ a keeper ]
[ Despite using the ∂S stuff correctly, at this time I was unaware of the fact that there are really 8 different accelerations and here is what happens below: (1) I compute a = aS by differentiation of v and I notice that there are terms present. (2) I compute the cross acceleration aS' and notice that there are no terms present. (3) I wrongly think of these two results as being the same acceleration, and I then have a paradox since one has stuff and the other does not ! ]
1. Start with the G rule for r:
[ I think this entire section is OK, though the result is not expressed in a useful manner. ]
v = (dr/dt)S = ω x r + (dr/dt)S' . [ correct !! ]
Our ant makes his voyage r(t) in frame S' which we write as
r(t) = r'1(t) e'1 + r'2(t) e'2 , [ ok if components mean (r)'i ]
where everything is S'-frame stuff. Then it is quite easy to compute
(dr/dt)S' = dt r(t) = '1(t) e'1 + '2(t) e'2 . [ since e'i constants; again ( )'i(t) components ]
Meanwhile,
ω x r = ω x r = [ωz] x [r'1(t) e'1 + r'2(t) e'2 ] = ω r'1 e'2 – ω r'2 e'1 . [ok]
We then end up with
v = (dr/dt)S = ω r'1 e'2 – ω r'2 e'1 + '1(t) e'1 + '2(t) e'2 [ok]
= ('1 – ω r'2) e'1 + ('2 + ω r'1 ) e'2 [ok, components must be (v)'i]
But now we want this in terms of frame S basis vectors,
e1' = R(φ) e1 = (e1' e1) e1 + (e1' e2) e2
e2' = R(φ) e2 = (e2' e1) e1 + (e2' e2) e2
e1' = R(φ) e1 = cosφ e1 + cos(π/2-φ) e2
e2' = R(φ) e2 = cos(π/2+φ) e1 + cosφ e2
e1' = R(φ) e1 = cosφ e1 + sinφ e2
e2' = R(φ) e2 = - sinφ e1 + cosφ e2
or
= = Rz(φ)T
Then we find
v = (dr/dt)S = ('1 – ω r'2) e'1 + ('2(t) + ω r'1 ) e'2
= ('1 – ω r'2) [cosφ e1 + sinφ e2]
+ ('2 + ω r'1 ) [- sinφ e1 + cosφ e2]
= [ ('1 – ω r'2)cosφ - ('2 + ω r'1 )sinφ ] e1
+ [ ('1 – ω r'2) sinφ + ('2 + ω r'1 ) cosφ] e2 [ ok unless algebra errors ]
Now let's try for more compact vector notation. How about
v = (dr/dt)S = (e'1, e'2)
(e'1, e'2) = (e1, e2) R(φ) = (e1, e2)
Then we get
v = (dr/dt)S = (e1, e2)
which agrees with the written-out result above. Notice there are no terms here.
Test Case: Suppose the ant is at rest in frame S'. This then says [(dr/dt)S'= 0 ]
v = (dr/dt)S = (e1, e2) = (e'1, e'2)
= -ωr'2 e'1 + ωr'1 e'2 [ would rather express this in terms of en ]
Now meanwhile
ω x r = [ωz] x [r'1(t) e'1 + r'2(t) e'2 ] = ω r'1 e'2 – ω r'2 e'1 = the above v.
In polar coordinates, r = r and we then get
ω x r = ωr = v
which seems right.
2. Differentiate v to get a
[ my general plan in this doc is to compare the "differentiation method" to the "G Rule method" to see if the two methods agree or not. So here we are doing the first method on v . ]
[ apart from possible algebra errors, I think this section is also OK.]
v = (dr/dt)S = (e1, e2)
dt = ω
dt =
// for access later: '1 - ω '2 - r'2 '2 + ω '1 + r'1
'1 – ω r'2 '2 + ω r'1
Then our result seems to be
a = (e1, e2) [ + ω ]
An impressive mess!
Write out as
a = a1 e1 + a2 e2
a1 = cosφ ('1 - ω '2 - r'2) - sinφ('2 + ω '1 + r'1) - ωsinφ ('1 – ω r'2) -ωcosφ ('2 + ω r'1)
a2 = sinφ ('1 - ω '2 - r'2) + cosφ('2 + ω '1 + r'1) + ωcosφ ('1 – ω r'2) -ωsinφ ('2 + ω r'1)
a1 = cosφ ('1 - ω '2 - r'2) - sinφ('2 + 2ω '1 + r'1) - ωsinφ ( – ω r'2) -ωcosφ ('2 + ω r'1)
a2 = sinφ ('1 - ω '2 - r'2) + cosφ('2 + 2ω '1 + r'1) + ωcosφ ( – ω r'2) -ωsinφ ('2 + ω r'1)
a1 = cosφ ('1 - 2ω '2 - r'2) - sinφ('2 + 2ω '1 + r'1) - ωsinφ ( – ω r'2) -ωcosφ (ω r'1)
a2 = sinφ ('1 - 2ω '2 - r'2) + cosφ('2 + 2ω '1 + r'1) + ωcosφ ( – ω r'2) -ωsinφ (ω r'1)
a1 = cosφ ('1 - 2ω '2 - r'2) - sinφ('2 + 2ω '1 + r'1) + ω2sinφ r'2 - ω2cosφ r'1
a2 = sinφ ('1 - 2ω '2 - r'2) + cosφ('2 + 2ω '1 + r'1) - ω2cosφ r'2 - ω2sinφ r'1
a1 = cosφ ('1 - 2ω '2 - r'2 - ω2 r'1) - sinφ('2 + 2ω '1 + r'1 – ω2r'2)
a2 = sinφ ('1 - 2ω '2 - r'2 - ω2 r'1) + cosφ('2 + 2ω '1 + r'1 – ω2 r'2)
=
Here is an alternate way to get the above result. Start again with
a = (e1, e2) [ + ω ]
Rewrite the second term as (you have to stare carefully to verify this, check all 4 items! )
-ω
Then we have
a = (e1, e2) [ - ω ]
= (e1, e2)
which replicates the result obtained above. Note the presence of terms!
3. Do the G rule for v:
We want to obtain the same result as above for a by using the G rule for v. That rule is
a = (dv/dt)S = ω x v + (dv/dt)S' . [ok]
So, ant makes this voyage [ normally would have used r', but OK to give voyage in terms of r ]
rant(t) = r'1(t) e'1 + r'2(t) e'2 [ r = (r)'i e'i, seems valid ]
vant(t) = '1(t) e'1 + '2(t) e'2 [ (dr/dt)S' = vS' = ()'i e'i ]
(dvant/dt)S' = '1(t) e'1 + '2(t) e'2 = (e'1, e'2) = (e1, e2)
[ ∂S'2r = ∂S' vS' = aS'S' = aS' = as shown on right above = one of the 6 cross accelerations! ]
and boom, we have the second [ first? ] term right away. There are no (or ω) factors anywhere. After all, this calculation is not supposed to know about the rotation.
STOP. Right away we have a problem Houston. Recall from above that we had
v = (dr/dt)S = (e1, e2)
which contains no terms. Thus ω x v contains no terms . And we just showed that (dv/dt)S' has no ω or terms. Therefore, our G Rule will give a result for a which has no terms. This then disagrees with the previous section's calculation of a. So we have now arrived at this paradox staying entirely with Cartesian coordinates, avoiding possible confusing extra issues with the polar coordinates.
[ A good paradox! Here is the resolution: the derivative method correctly computed a = aS , but in the G Rule method we computed the cross acceleration aS' which is a different object!!! ]
******************************************************************************
In the first term, we shall use the exact v from the previous problem, which was this
v = (dr/dt)S = (e1, e2)
Then,
ω x v = ω x (e1, e2)
= ω (e2, -e1)
and then our full result becomes
a = (dv/dt)S = ω (e2, -e1) + (e1, e2)
But now compare this with our result from Section 2 above
a = (e1, e2)
This Section 2 result includes factors , but our G Rule for v produces no such factors. So the results are in fact NOT consistent. Why is this? (My 3 week long problem continues, even in Cartesians)
[ see resolution above ]