Nov 13 to Dec. 31, 2016 - Over 6 weeks. As promised, here's my answer (first part anyway) to the problem posed in #1:
It will simplify to have the common transmission shaft linking the LHS motor and RHS generator to be a thin-walled tube of mean radius r, instantaneous angular velocity
ω, and subject to an axially uniform torque
τ generating a spatially uniform shear stress of magnitude
σ everywhere within the tube wall.
As shown in the figure, consider a thin annular section of the shaft wall - a ring - of width b, and in particular a small square shaped segment of it - wall element s of edge lengths a, two such edges oriented parallel to the shaft axis.
With the LHS sourced driving torque
τ directed as shown, antiparallel edge traction forces
f1,
f2, acting on s to form a CW couple, will be as shown. Assuming negligible shaft elastic energy storage capacity, power notionally flowing in to element s via
f1 is exactly compensated by an equal withdrawal via
f2. Being in shear equilibrium, an opposing ACW couple owing to axially oriented antiparallel edge forces
f3,
f4 must also be present. Which pair being normal to the peripheral velocity
v =
ω × r of s play no role in the notional power flows. In the lab frame S where system COM is stationary, static equilibrium requires at all times |
f1|=|
f2|=|
f3|=|
f4|, whether
τ or
ω are time varying or not (restricting evaluation to where wave propagation i.e axial phase delay is at all times negligible.)
Initial clue in #1 was to ponder
F = d
p/dt. It implies any shaft axial force will necessarily be directly proportional to time-rate-of-change of power flow dP/dt = d(
τ.ω)/dt from source to load. Hence zero whenever P =
τ.ω is constant not just zero. Suggesting examining just two complimentary cases:
1: Fixed
ω, ramp torque d
τ/dt,
2: Fixed
τ, ramp angular velocity d
ω/dt
Case 1 first. In lab frame S, time varying
τ thus
σ thus edge forces
f1,
f2,
f3,
f4 offer no evident avenue for an axial force as they all vary uniformly. But now move into the instantaneous proper frame S' of s, in S moving vertically up at
v =
ω × r as seen in the figure. In S' the relevant quantity is the temporal term of Lorentz boost:
t' = γ(t-
v.x/c²) - (1)
As only the pair
f3,
f4, have a separation distance b along
v, only they are affected by (1) such that a differential in magnitude exists between them in S'. Their equal purely time rate of change d
f/dt in S now having a spatial gradient added in S'. Setting t = 0, x = 0, for
f3 in S', gives
f3' =
f3'(t'=0) - (2)
f4' = -(
f3'(t'=0) + γd
f/dt(-
v.b/c²)) - (3)
δ
f' =
f4'+
f3' = γd
f/dt(
v.b/c²) - (4)
The γ factor in (4) is unimportant and cancels out on transforming back into frame S. However the crucial residual force δ
f, owing to nonsimultaneity does not and is responsible for 'rescuing' conservation of momentum. it also nicely acts against the nominal rate of change of power flow, as required. The force density is just (4) divided by the separation distance b, (and γ factored back out ) i.e.
|
ρ| = |δ
f/b| = |d
f/dt(v/c²)| - (5)
Which by simple inspection matches and cancels the rate of change of momentum transfer through element s heading from motor to load. Obviously since their is a balance across any given element of shaft, integrating over the whole will maintain that situation.
How about that. No recourse to heavy higher maths. I did leave out the connection to shear waves where body forces explicitly are proportional to the shear stress gradient, but the reader can pick that sort of thing up in e.g. eq'n. (1) here:
http://petrowiki.org/Compressional_and_shear_velocities
The interesting and subtle difference to this case is that the shear gradient is rotated orthogonal relative to the much more familiar shear wave situation - as needed.
Sorry but lots of trouble to deal with last few days threw my plans out, so no time to get case 2: dealt with before 2017 forum time arrives. While I'm recovering, maybe at last someone here can now figure out how to do that one. Back later....