Making the relative motion of the frames more general
DOCX · 85.7 KB
Open DOCX file
Short note by Phil dated 7.20.12, written for his frames document. It generalizes where the rotation axis ω sits relative to frames S and S', derives de'n/dt = ω x e'n from an infinitesimal rotation, and poses the system of linear ODEs for finding R(θ) given ω(t). He tries an elimination, concludes there is probably no analytic solution, and sets it aside. It ends with two special cases: a tumbling camera platform and the Coriolis setup on the earth.
AI-written summary; may contain errors.
Extracted text (machine-read; may contain errors)
Making the relative motion of the frames more general? PhL 7.20.12
[ Here I generalize the location of the ω axis and then define my two special cases, all as now appear in frames doc.
[Also, I state the system of ODE's problem of how to calculate R(θ) given ω(t) which I don't think has an analytic solution even though only first order. You need to solve a system of 3 first order ODE's with variable coefficients. ]
Consider this drawing showing Frame S and Frame S' :
Frame S is aligned with the paper, with unit vector e3 coming directly out of the plane of paper.
The origin of frame S' lies on the end of a vector b whose tail lies on the origin of frame S. This vector b may be varying in time.
Frame S' is instantaneously rotating about some axis defined by vector ω = dφ/dt.
We make the following computation to learn how the basis vectors e'n are rotating relative to the axes of frame S (see Appendix A for the notation used here),
e'n(t+dt) = R(dφ) e'n(t) where R(dφ) = some small rotation about the ω vector axis
de'n = e'n(t+dt) - e'n(t) = [R(dφ) - 1] e'n(t) ≈ -idφ J e'n(t) = -idφk Jk e'n(t)
=> (de'n)i = -idφk (Jk)ij (e'n(t))j = - dφkkij (e'n(t))j = εikj dφk (e'n(t))j
=> de'n = dφ x e'n
=> (de'n/dt)S = ω(t) x e'n n = 1,2,3 where ω ≡ dφ/dt [ usual result]
Statement of the problem of finding θ(t) given ω(t)
This last equation is in fact a system of nine first-order, homogeneous, ordinary, linear differential equations in variable t in which there are nine unknown functions (e'n)i, where n = 1,2,3 and i = 1,2,3,
[(de'n/dt)S]i – εijkωj(t) (e'n)k = 0 n = 1,2,3 i = 1,2,3
In principle, given ω(t) and given some initial conditions, one could solve this to determine e'n(t) at any time t. This result could then be written in the form
e'n(t) = R[θ(t)] en
since at any time t the frame S' axes are just some rotation applied to those of frame S. That is to say, by solving the problem as just described, we would then know the parameter vector θ(t) and it would then be a function of ω(t) and the initial (boundary) conditions. Here is the problem in more conventional form
f'(n)i(t) + Σk Aik(t) f(n)k(t) = 0 Aik(t) = Σj εijkωj(t) = antisym
+ Σk Aik(t) fnk(t) = 0 n,i = 1..3
For some general ω(t), this is a non-trivial problem! Well, for each fixed n, it is a system of three equations. If we suppress the n, it becomes
+ Σk Aik(t) fk(t) = 0 n,i = 1..3
I think if you write this all out and play with it, it is solvable.
1 + A12u2 + A13u3 = 0 A12 = ε1j2ωj = ε132ω3 = - ω3
2 + A21u1 + A23u3 = 0 A13 = ε1j3ωj = ε123ω2 = ω2
3 + A31u1 + A32u2 = 0 A23 = ε2j3ωj = ε213ω1 = - ω1
1 - ω3u2 + ω2u3 = 0
2 + ω3u1 - ω1u3 = 0
3 - ω2u1 + ω1u2 = 0 - ω x u = 0
ω11 - ω1ω3u2 + ω1ω2u3 = 0
ω22 + ω2ω3u1 - ω2ω1u3 = 0
Add to get u3 eliminated.
ω11 + ω22 + ω2ω3u1 - ω1ω3u2 = 0
ω11 + ω22 + ω3[ω2u1 - ω1u2] = 0
Maybe define
v = ω2u1 - ω1u2
= 2 u1 – 1 u2 + ω21 - ω12
Well, I don't see how to get things separated here. Since Aij is not symmetric, the equation is probably not diagonalizable. We are off-topic, so let this be "till a rainy day". Web scan suggests that there is no analytic solution.
Special cases of ω positioning
The picture supports many different scenarios, here are two: [ my special cases ]
In one scenario, the ω axis always passes through the origin of frame S' so frame S' just tumbles about its origin in some manner determined by ω(t) while the frame S' origin simultaneously translates relative to frame S according to some b(t) . In this scenario, one might imagine that Frame S' is a self-propelled camera platform that is observing some activity in frame S.
In another scenario, suppose the ω axis passes through the origin of frame S, is aligned with e3, and ω is a constant. If the origin of frame S were put at the center of the earth and if b were a constant-length vector to the surface of the earth, and if perhaps e'1 pointed "east" when e'2 pointed "north" and e'3 pointed "up", the scenario would be suitable for studying Coriolis forces for physics on the surface of the earth.